2. 地层与断层数据基础

在ArcGIS Pro中做地层与断层的空间叠加分析,第一步不是打开工具,而是把地质对象"翻译"成GIS能理解的数据模型。地层界面、断层线、断层带这三类对象在野外是三维地质实体,在GIS里却必须落到点、线、面、栅格或TIN等具体载体上。数据表达选错了,后面的叠加、相交、缓冲区分析都会失真。本节围绕"对象—数据模型—适用场景—组织方式"这条主线,把基础打牢。

三类地质对象的数据表达

地层界面本质是一个地质分界面,最自然的表达是面要素或栅格表面。当作平面分布图使用时,用面要素(Polygon)表示某地层出露范围或顶面投影范围,属性表记录地层代号、岩性、时代;当作起伏面使用时,用栅格(DEM或地层顶面高程栅格)或TIN表达界面高程,每个像元值即该点界面标高。判断依据很简单:若分析关注"某地层在平面上覆盖哪里",用面;若关注"界面埋深、厚度、起伏",用栅格或TIN。

断层线是断层面与地表(或某一基准面)的交线,标准表达是线要素(Polyline)。属性上至少要有断层编号、性质(正断层、逆断层、平移断层)、走向、倾角、倾向、落差等字段。断层线是叠加分析中最常用的对象,因为它能直接与地层面相距、相交、做缓冲区。

断层带是断层两侧岩层破碎影响的范围,具有宽度,因此用面要素(Polygon)或双线表达。断层带面可以由断层线按破碎带宽度做双侧缓冲区生成,也可以直接数字化实测边界。它适合做"断层影响范围"类分析,例如评估某地层在断层带内的分布比例。

地质对象 → 数据模型 → 适用场景 地层界面 地质分界面 有范围、有起伏 断层线 断层面与基准面交线 有走向、倾角、落差 断层带 破碎影响范围 有宽度、可双侧扩展 面要素 Polygon 栅格 / TIN 表面 线要素 Polyline 面要素 Polygon 线 + 缓冲区生成 平面分布分析 地层出露范围、面积统计、与断层相交 → 用面要素,属性含地层代号、岩性 界面起伏与厚度分析 埋深、等厚图、顶面高程、体积估算 → 用栅格/TIN,像元值即界面标高 断层与地层空间关系 相交、相距、切割、缓冲区影响范围 → 断层线 + 地层面,做空间叠加 断层影响带评估 破碎带内某地层占比、避让范围 → 断层带面,与地层面叠加统计

矢量与栅格的适用场景判断

矢量与栅格不是谁更先进的问题,而是各自适配不同的分析问题。矢量(点、线、面)保留精确边界和拓扑关系,适合做边界类、关系类分析:地层与断层是否相交、断层线穿过哪些地层、断层带覆盖的地层面积比例。这类问题用"相交(Intersect)""空间连接(Spatial Join)""按位置选择"就能得到干净结果,且属性表可直接统计。

栅格适合做连续场、逐像元运算:地层顶面高程、厚度等值面、断层距离场。例如要计算"距断层500米内的地层厚度均值",可先由断层线生成欧氏距离栅格,再与地层厚度栅格做栅格计算器运算。栅格的优势是每个像元都有值,便于连续统计;代价是精度受像元大小控制,边界会有锯齿。

实操中常见的组合是:矢量定关系,栅格算数值。先用断层线与地层面做相交,得到被断层切割的地层分段;再把地层顶面栅格与断层距离栅格叠加,评估断层附近界面起伏变化。两者互为补充,而不是二选一。

在ArcGIS Pro中的组织方法

数据组织直接决定后续叠加分析是否顺畅。建议按"工程—地理数据库—要素数据集—要素类/栅格数据集"的层级管理,而不是把Shapefile散落在文件夹里。具体做法:新建一个文件地理数据库(如 Geology.gdb),在其中建立要素数据集 FaultStratum,统一指定投影坐标系(例如项目区适用的3度带高斯-克吕格投影),再把断层线、断层带、地层界面面要素放入同一数据集。同一数据集内坐标系一致,做相交、缓冲区时不会因坐标系不一致而报错或产生偏移。

栅格数据单独放在 Raster 文件夹或镶嵌数据集(Mosaic Dataset)中,命名体现内容与分辨率,例如 Stratum_Top_10m、Fault_Dist_10m。关键字段要提前设计好:断层线至少含 FaultID、Type、Strike、Dip、Throw;地层面含 StratumCode、Lithology、Age。字段类型用文本或短整型,避免中文别名与字段名混用导致工具识别异常。

组织检查要点:① 所有参与叠加的矢量是否在同一投影坐标系;② 断层线是否为单部件线、无自相交;③ 地层面是否存在缝隙或重叠;④ 栅格像元大小是否与分析尺度匹配(如1:1万填图宜用5~10米栅格)。这四项在叠加前处理干净,能避免大部分拓扑错误和面积统计偏差。

一个可复现的组织示例

假设某项目区有1:5万地质图,需分析断层对某地层的影响。示例参数如下(属教学假设,非实测):断层线要素类 Fault_Line,共42条;地层界面面要素类 Stratum_A,投影为CGCS2000 3度带;断层破碎带宽度按性质取值,正断层取80米、逆断层取120米。操作上,先对 Fault_Line 按 Type 字段筛选,分别做双侧缓冲区生成断层带面;再用"相交"工具将断层带面与 Stratum_A 叠加,输出要素类保留双方属性,即可用 StratumCode 分组统计断层带内该地层的面积与占比。判断结果是否合理,可检查叠加后总面积是否等于原断层带面积(允许微小容差),若明显偏大,多半是面存在重叠或坐标系不一致。

把这三类对象的数据表达、矢量栅格的适用边界、以及地理数据库中的组织规范理清之后,后续的相交、缓冲区、栅格计算才有可靠输入。数据基础的质量,往往比工具参数更决定叠加分析的结论是否可信。