2. 断层数据解读与预处理
在断层缓冲区分析中,断层线数据的质量直接决定后续缓冲距离、矿点归属判断和成矿规律解释的可信度。很多分析结果“看起来不对”,根源往往不在缓冲区参数,而在进入分析前的断层数据本身:几何类型混杂、坐标系缺失、属性字段含义不清、悬挂节点与自相交未处理。本节围绕一条典型的断层线要素类,讲清楚如何读懂它、检查它、修正它,使其成为可直接参与缓冲区运算的干净数据。
一、先读懂数据:几何类型与属性字段
打开断层要素类后,第一件事是确认几何类型。ArcGIS Pro 中要素类几何类型分为点、线、面、多点和注记,断层通常为线(Polyline)。但要注意“线”内部还有差别:单部件线(Singlepart)与多部件线(Multipart)。一条断层若由若干断续段组成,可能被存成一个多部件要素,此时对整条要素做缓冲会得到一条连续带,而实际断层在段间存在缺口。判断方法是在属性表中查看要素数量与地图上可见段数是否一致,或用“多部件转单部件”工具拆分后对比。几何类型判断错误,后续缓冲距离和相交统计都会失真。
属性字段方面,断层数据常见字段包括断层编号(如 Fault_ID)、断层名称、断层性质(正断层、逆断层、走滑断层)、走向、倾角、长度、数据来源与精度等级。其中断层性质影响缓冲策略:活动断裂常需更大缓冲半径,而一般断层可按固定距离处理;走向字段可用于校验几何方向是否与地质认识一致。需要特别检查字段类型:编号若被存为数值型,前导零会丢失;长度字段若为文本型,则无法直接参与计算。建议在预处理阶段统一为:编号用文本型,长度、走向、倾角用双精度型,性质用文本型并建立值域约束。
二、数据检查:从坐标系到几何异常
坐标系检查是预处理的第一道关口。在 ArcGIS Pro 中查看要素类属性,确认是否定义了投影坐标系。若显示“未知”或仅有地理坐标系(如 GCS_WGS_1984,单位为度),则不能直接以“米”为单位设置缓冲距离。此时应使用“投影”工具转换为适合研究区的投影坐标系,例如我国中东部常用 CGCS2000 3 度带高斯-克吕格投影,或按研究区选择的 UTM 带。转换时需明确源坐标系与目标坐标系,并保留原始数据备份,避免反复转换引入误差。
几何检查重点看三类问题:一是悬挂节点,即断层线端点未与其他线相接,形成“断头”;二是自相交,线要素自身交叉成环;三是重复线,同一断层被数字化两次。ArcGIS Pro 的“检查几何”工具可批量识别空几何、自相交、短线段等问题,结果写入几何检查表,逐条定位修复。对于悬挂节点,若属于数字化遗漏,可用“延伸线”工具按容差延伸至相邻线;若本身是断层端点,则不应强行连接,需在属性中标注端点性质。
三、投影转换:让缓冲距离有真实单位
缓冲区分析的本质是“以距离为半径生成面”,距离必须有明确的长度单位。若断层线处于地理坐标系,其坐标单位为度,1 度经度在不同纬度对应的地面距离差异很大,直接按“0.01”设置缓冲会得到随纬度变化的不一致结果。因此,进入缓冲区之前必须完成投影转换。操作路径为:数据管理工具 → 投影与变换 → 投影,选择输入要素类,指定输出坐标系(如 CGCS2000_3_Degree_GK_Zone_39),并确认地理变换参数(若涉及不同基准面)。转换完成后,在要素类属性中核对线性单位是否为米,并用“计算几何”工具重新计算断层长度字段,与转换前按度计算的结果对比,确认量级合理。
需要提醒的是,投影转换不是“越转越好”。若研究区跨两个 3 度带,应选择覆盖全区的中央经线,或改用 Albers 等面积投影以兼顾面积统计。转换后应保留原始数据副本,并在元数据中记录源坐标系、目标坐标系和转换时间,便于回溯。
四、拓扑修复:建立规则并逐条处理
拓扑修复建议在文件地理数据库的要素数据集内进行。新建拓扑,添加断层线要素类,设置规则:不能有悬挂节点(针对应连接的断层段)、不能自相交、不能有重复线。验证拓扑后,错误 inspector 会列出每处异常及涉及要素。处理原则是“先判断地质合理性,再决定修复方式”:若悬挂节点是数字化遗漏,用延伸或捕捉连接;若自相交是线段折返,用“重塑”或节点编辑修正;若重复线确属同一断层,保留精度较高的一条并删除另一条,同时在属性中记录处理说明。
修复完成后再次验证拓扑,直至错误数为零。随后用“要素类转要素类”或“复制要素”生成一份干净数据,命名建议体现处理状态,例如 Fault_Line_Projected_Clean,避免与原始数据混淆。此时数据已具备三个条件:几何类型明确、坐标系为投影米制、拓扑无错误,可以进入缓冲区距离设定与矿点空间关系分析。