李金亮

（鹤壁水文水资源勘测局 河南鹤壁 458000）

在研究水资源污染现象中发现［1-3］，该现象的形成原因较为复杂，主要是由于在一个城市相同的时间段或不同的城市相同的时间段内，灾害因子的表现形式不同［4］，同时，水资源污染事件的发生具有不确定性、阶段性、突发性等特点［5］。处理水资源污染事件前，需要工作人员进行大量的水文、水质、地质工程的勘查，掌握水环境基本特征。为了实现优化水资源污染实时监测工作，水环境监督管理单位提出了基于水环境过程变化动态模拟的污染模型，并明确了将数字化技术作为支撑执行此项工作。但是，该方法构建的模型为水环境二维模型，存在一定的局限性，且无法及时感知突发性污染问题［6］。

1 GIS水文建模的应急监测方法

1.1 建立水资源空间数据库

为实现应急监测水资源污染，引进MySQL 关系数据库［7］，通过获取水资源空间数据的方式，构建水资源环境空间数据库。在此过程中，将空间数据以离散数据或集成数据的方式存储在关系数据库中，并使用空间数据决策引擎工具（ArcSDE 处理工具）导入与存储数据。同时，通过Arcview 工具，可视化处理地质信息与地理环境信息，在存储终端设置不同的空间模块，存储空间数据与属性数据［8］。在此过程中，通过VIEW管理采集或获取的空间数据，将不同类型的地理信息按照实体类型划分为多个信息结构并按照类别表示。

在构建数据库时，需要建立属性数据库，管理各类属性数据信息，以表的形式，结合关系模型，实现对各项属性数据信息的连接。针对每一个图层，均建立相应主题的属性列表，同时，针对各类属性数据，建立相应的检索关键词，并利用Link工具，通过对相同关键词的获取，实现其与主题属性之间的连接。为了确保两种数据库的相互协调［9］，单独存储数据，并通过Arc-SDE 的应用程序接口，实现与空间数据库各数据信息相连，确保监测过程中各类数据能够更加快速地完成传输、查询和输出。

1.2 GIS水文建模的污染源定位

在建立水资源空间数据库后，引入GIS 水文建模［10-11］，定位水资源具体污染和污染源位置。图1为基于GIS水文建模的水资源污染源定位流程示意图。

根据图1可知，通过GIS水文建模，展示地级市县、县级市界等区域内水资源环境的空间信息，将收集的空间数据和属性数据显示在相应水资源环境划分区域内，同时，在模型中明确标记水资源周围环境中的排污口、排污单位等。结合监测输入的数据及参数，预测监测区域内的水资源未来水质条件，并动态模拟单污染源或多污染源的叠加影响效果［12］。在ArcView 窗口中，针对模拟后的水环境水质进行可视化展现。

图1 基于G IS 水文建模的水资源污染源定位流程示意图

1.3 动态数据的水污染扩散分析与应急监测

结合动态数据分析水污染扩散情况，通过扩散约束条件应急监测水质。在分析水污染扩散时，结合经典水动力学方程和水质扩散算法，模拟真实的水环境条件下水资源的运动。水资源在流动过程中具有典型的连续性，因此，通过下述公式实现对水资源运动的表达：

式（1）中，h表示水资源的实际深度，t表示监测过程中水资源的流动时间，u表示水资源流动过程中在水平方向上的平均流速，v表示水资源流动过程中在竖直方向上的平均流速，x表示水环境周围水平方向风速，y表示水环境周围竖直方向风速。

在此基础上，针对水资源污染扩散的分析，需要引入动量方程和非持久性污染物迁移方程，基于有限体积法，实现对其离散处理。为了确保监测结果的精度，在DEM网格中，采用SIMPLEC法校正水深。在分析过程中，由于水资源持续流动，因此，结合干湿法，处理水资源的动态边界问题。同时，划分水资源中污染等级，将水资源污染等级按照从重到轻依次划分为6 个等级，分别为I等级、II等级、III等级、IV等级、V等级和V劣等级。其中，I～III 等级表明在水资源中污染物没有或含量极少，不会威胁到水资源的正常使用和运行；IV等级为水资源中含有少量污染物，并且已经出现影响水资源正常使用的情况；V等级为无法使用的水资源；V 劣等级是指水资源中存在大量污染物，并且需要紧急治理，确保水资源的正常使用。6 个等级分别设置不同的阈值，根据上述扩散分析得出模拟预测结果，与各个阈值范围进行比较，实现快速分类水资源污染等级，并给出不同等级污染的应急保护、治理措施，降低水资源污染的影响，以此实现对水资源污染应急监测。

2 对比实验

在引入GIS 水文建模技术后，提出了一种全新的水资源污染应急监测方法。为了进一步验证研究方法的合理性和可行性，选择新的监测方法作为实验组，将传统水资源污染监测中的基于分布式水文模型的监测方法作为对照组，将两种不同的监测方法应用到相同的实验环境中，并针对相同的监测对象进行对比实验研究。

选择在具有模型计算功能的软件中输入真实的污染物地点、污染排放量、水资源污染浓度等数据，模拟水资源中污染物扩散情况。通过该软件，获取坐标文件和污染浓度文件，并且坐标文件的信息对应污染物浓度文件的信息。在模拟水资源污染扩散的过程中，分别利用提出的监测方法和基于分布式水文模型的监测方法对其进行监测。表1为模拟的水资源环境水质条件记录表。

表1中，IV、V和V劣均为水资源的水质等级，污染程度从小到大排列依次为IV＞V＞V劣，其中，IV等级表明水资源污染程度较轻，而V 劣等级表明水资源污染程度最严重。从表1中记录的各个区域水资源环境水质条件可知，在该水资源环境范围内，区域A、区域C和区域E均在不同月份出现严重污染的情况。

表1 模拟的水资源环境水质条件记录表

2.1 监测精度对比

为验证设计监测方法在实际应用中的监测精度，将表1记录的内容作为基础，利用本文监测方法和基于分布式水文模型的监测方法，测定5 个区域内的污染物浓度，并将得出的实验结果绘制成表2所示。

从表2的数据可知，实验组的监测结果更接近于实际测量结果，最大误差仅为0.04mg/L，其他误差均在0.01mg/L，并且监测结果能够为水资源污染等级的划分提供重要依据。同时，在实验过程中发现，对照组监测方法在测定5 个区域的水资源污染物浓度时，区域A、区域C和区域E的监测结果与其他两个区域相比略接近实测值，说明对照组监测方法适用于监测存在重度污染的水资源环境的污染浓度，并且监测的结果均没有达到实验组的监测精度。因此，实验初步证明，本文提出的基于GIS水文建模的监测方法在实际应用中能够实现高精度监测水资源污染浓度，提高了监测精度，并且监测结果能够为水资源水质等级划分提供重要数据依据。

表2 实验组与对照组水资源污染浓度监测结果记录表（单位：mg / L）

2.2 监测时效性对比

为验证该方法是否具备快速响应的应急能力，对比实验组和对照组两种监测方法的时效性，监测5 个区域的水资源污染浓度变化，通过人为的方式，改变水资源污染扩散模拟过程中污染物的浓度。分别利用两种监测方法对水资源环境中污染物浓度瞬间变化时的平均响应时间作为评价指标，指标值越小，表明方法的性能越好。瞬时变化响应时间的计算公式为：

式（2）中，表示水资源环境污染浓度瞬间变化时的监测方法的平均响应时间，t1表示监测方法响应时刻，t2表示监测方法发现浓度异常变化时刻，n表示监测方法从发现浓度异常变化时刻到作出响应时刻之间监测次数。根据上述公式，计算得出实验组和对照组两种监测方法在水资源环境污染浓度瞬间变化时的平均响应时间，并将计算得出的结果绘制成图2所示。

从图2中两条曲线可知，实验组监测方法的平均响应时间在40ms 以下，对照组监测方法的平均响应时间超过50ms，两者相比，实验组比对照组低了10ms以上，实验组的平均响应时间更短，因此，本文监测方法在实际应用中具备更好的时效性。本文基于GIS水文建模的监测方法与基于分布式水文模型的监测方法相比，前者的监测精度和监测时效性更强，其能够为水资源保护和水资源治理措施的提出提供可靠的数据依据。

图2 实验组和对照组监测时效性对比图

3 结语

由于人类社会活动的不规范、对土地资源与环境资源的无节制开垦，全球生态环境建设呈现极端变化趋势，水资源被污染的面积越来越大，该现象对于人类生活所造成的负面影响呈现一种增加趋势。因此，针对上述问题，本文在现有成果的基础上引进GIS技术，采用对水资源进行水文建模的方式，研究具有针对性的水资源污染的应急监测方法，通过该方法，实现对突发性污染事件的及时感知与处理，为环境治理和保护提供可靠依据。