2. 矿区水系数据获取与预处理

在矿区水文地质与构造解译工作中，水系是最直观、最容易与地形和构造建立联系的一类要素。它既反映地表径流的汇聚路径，也常常沿断裂、节理或岩性接触带发育，因此水系数据的质量直接决定后续叠加分析的可信度。本节要完成的任务是把三类来源不同、格式不同的水系信息整合为统一、可用的矢量数据：一是从DEM自动提取的水系，二是对已有地形图或影像进行矢量化得到的河流与沟渠，三是野外调查记录的水系点。三者最终要在同一坐标系下相互印证，并通过完整性检查。

一、从DEM提取水系：理解算法与参数

从DEM提取水系的核心是“填洼—流向—汇流累积—阈值提取”这条链路。以30 m分辨率的SRTM或ASTER GDEM为例，首先进行洼地填充（Fill Sinks），避免因数据噪声形成虚假闭合洼地导致流向中断。随后计算D8流向，即每个栅格流向其8邻域中坡度最陡的方向。再基于流向计算汇流累积量（Flow Accumulation），每个像元的值代表其上游汇水像元数。最后设定阈值，将累积量大于阈值的像元提取为河网栅格，再栅格转矢量得到线状水系。

阈值的选择直接决定河网密度。以30 m栅格为例，若阈值设为1000个像元，约相当于0.9 km²的汇水面积，提取结果偏向主干河道；若设为200个像元（约0.18 km²），则能捕捉更多冲沟，但也会引入大量伪河道。矿区地形常受采动、排土场和道路改造影响，建议先用较大阈值提取主干，再用较小阈值补充沟渠，并与影像比对后取舍。需要强调的是，DEM提取的水系是地形驱动的“理论汇流路径”，并不等同于实际有水的河道，尤其在干旱矿区，很多提取线只是季节性冲沟。

二、河流与沟渠的矢量化

对于已有地形图、遥感影像或矿区排水系统图的对象，矢量化是获取真实河道与人工沟渠的主要手段。矢量化时应遵循“先分类、后分层”的原则：天然河流、季节性冲沟、人工排水沟分别放在不同图层，并在线属性中记录类型、宽度等级、流向和是否常年有水。线要素的绘制方向应统一为从上游到下游，便于后续计算流向和长度统计。对于宽度较大的河流，可先绘中心线，再在需要时缓冲为面；对于矿区排水沟，应尽量沿沟壁中心绘制，并在转折处增加节点以保证几何贴合。

矢量化过程中最常见的质量问题是悬挂线、自相交和重复线。悬挂线指端点未与其他线或边界连接的孤立线段，会造成网络拓扑断裂；自相交会使后续网络分析失败；重复线则会导致长度统计翻倍。建议在编辑时开启捕捉（ snapping），容差设为略小于最小沟渠宽度，例如0.5 m；完成后再用拓扑检查工具逐项排查。

三、野外调查水系点的整理

野外调查通常以GPS或手机记录泉点、井口、地表水取样点、断头沟、渗水点等位置，附带编号、坐标、高程、流量描述和照片。整理时首先要统一字段结构，例如“点号、类型、X、Y、高程、流量描述、调查日期、备注”。其次要检查坐标列是否被表格软件误读为文本或科学计数法，导入QGIS时应明确指定X字段为经度、Y字段为纬度，并选择与记录一致的坐标系（野外手持GPS多为WGS 84）。若记录的是度分秒格式，需先转换为十进制度。

野外点与提取水系、矢量化水系之间应建立空间对应关系。例如泉点应落在沟谷线上或其附近，若偏离超过合理距离（如50 m），需复核是GPS误差还是点位记录错误。这一步是后续“水系—构造”关联分析的基础。

四、统一坐标系与完整性检查

三类数据往往坐标系不一致：DEM多为地理坐标（WGS 84），矢量化底图可能是地方投影坐标，野外GPS为WGS 84。统一坐标系时应选择矿区常用的投影坐标系，例如CGCS2000 3度带高斯-克吕格投影，以保证长度和面积量算准确。在QGIS中可通过“重投影图层”或导出时指定目标CRS完成，注意不要仅修改图层显示CRS而不真正重投影。统一后应检查所有图层的范围是否重叠、单位是否为米。

完整性检查建议覆盖以下方面：几何完整性（空几何、自相交、悬挂线）、属性完整性（类型、流向字段是否为空）、空间完整性（水系是否覆盖整个矿区、野外点是否全部落入研究区）、逻辑完整性（河流流向是否与DEM高程递减方向一致）。下面这张图概括了三类水系数据从获取、处理到统一入库并接受完整性检查的整体关系。

五、一个可复用的检查清单

为便于实际执行，可将上述要点整理为一张检查表，在数据入库前逐项确认。下表字段均围绕本节任务设置，可直接作为作业记录使用。

完成上述步骤后，得到的是三类相互独立又彼此印证的水系数据：DEM提取线提供地形汇流骨架，矢量化线提供真实河道与人工沟渠，野外点提供实地验证。它们在同一投影下叠加，既可用于检查提取结果的合理性，也为后续与断层、岩性界线等构造信息叠加打下基础。需要提醒的是，矿区水系受人为改造影响大，任何自动提取结果都应视为待验证的初稿，最终以影像判读和野外核实为准。