1. 引言
GMS是目前国际上最先进的综合性地下水模拟软件包，其软件模块多，功能全，几乎可以模拟与地下水有关的所有水流和溶质运移问题 [1] [2] [3] 。随着科学技术的快速发展，GMS版本和功能也进行了升级和完善，补充许多新的应用程序，不断完善各模块的功能，其模拟运行的精确度也不断提高 [4] 。目前，国内外已有一些关于运用GMS模拟地下水溶质运移的研究 [5] [6] [7] 。
地下水水质的分析与研究是为了满足用水户不同的水质需求 [8] 。对存在的地下水水质问题寻求可能的解决方法与对策，从而实现水资源的持续利用 [9] [10] 。根据反映地下水物理化学特性的多个指标进行地下水水质评评价，通过数学方法的计算，得到最终适用的地下水水质等级评价结果 [11] 。地理信息技术等更为区域地下水水质演变与地下水水质演变成因分析的实现提供了可能 [12] [13] 。监测地下水质动态，建立水质资料数据库，预报地下水中污染物是国内外发展趋势 [14] [15] 。
本文以铜街–曼家寨矿区东部废石场为例，在研究场区水文地质条件的基础上，运用GMS中的MODFLOW模块建立三维数值模型，进行地下水流场数值模拟。针对废石淋滤液下渗情景下对地下水环境造成的影响进行预测分析，选用MT3D模块进行污染物溶质运移数值模拟，研究污染物在地下环境中的浓度变化情况。通过地下水环境影响数值模拟，预测地下水可能受到污染的情况。
2. 研究区概况
2.1. 研究区概况
铜街–曼家寨矿区位于云南省文山壮族、苗族自治州马关县都龙镇，设计开采标高为+1420~885米，开采方式为露天开采。矿山现有一个废石场，根据废石浸出毒性检测结果，废石属第I类一般工业固体废物，废石场位于曼家寨采矿场东部，占地面积115.74 hm2，库容8700万m3。
2.2. 区域水文地质条件
区域含水层组主要有岩溶–裂隙水、裂隙水含水层组两大类：
1) 碳酸岩盐岩溶–裂隙水含水层组
分布于都龙西部，面积31.14 km2，可划分为三个富水等级：∈3t富水性强、∈3x富水性中等、∈2l富水性较弱。
2) 裂隙水含水层组
分布于都龙以东及北部地区，面积68.86 km2，可划分为三个富水等级：∈2t富水性较弱、∈1ch富水性弱、花岗岩(
、
)富水性强。
区内隔水层为田蓬组第一段、第五段(∈2t1、∈2t5)、新鲜花岗岩(
、
)等地层。它们将各含水层组隔离，使测区内形成一个具多层结构的含水层系统。
2.3. 水文地质单元划分
本项目将研究区的水文地质单元划分为6个水文地质块段，分别命名为四台坡水文地质块段(I1)、南加水文地质块段(I2)、南当厂大沟文地质块段(I3)、铜街–曼家寨水文地质块段(II1)、烟山水文地质块段(II2)、马鞍山–大寨水文地质单元(III)。如图1所示。东部废石场位于铜街–曼家寨水文地质块段(II1)与四台坡水文地质块段(I1)两个水文地质单元地块之间(图1)。
Figure 1. Area map of hydrogeological units in the survey area
图1. 调查区水文地质单元分区图
3. 地下水流数值模拟和污染物溶质运移数值模拟
3.1. 水文地质概念模型
使用GMS软件进行数值模拟，其概念模型的建立主要包括研究区的划定及概化、边界条件的确定及水文地质参数的赋值(见3.2)。
3.1.1. 网格剖分
模拟区范围为13.8 km2，平面上将模拟区剖分为861,813个网格，垂向上分为2层，其中模型有效网格270,738个，每个网格长宽均为10 m，模拟区平面和空间三维模型如图2所示。
3.1.2. 模拟流场及初始、边界条件
为了更为准确概化模型区水文地质条件，模型在确定边界时遵循尽可能采用自然边界的原则。预测区北东部马关–都龙断裂，南部辣子寨、四台坡次级分水岭均概化为零流量边界，南部万龙山尾矿库概化为定水头边界，西部F2断裂概化为弱透水边界，次级分水岭概化为零流量边界。
3.2. 水文地质参数确定
3.2.1. 渗透系数
本次模拟将模拟区含水层进行了概化，将目标含水层区域渗透系数主要划分为六个分区，分区结果如图3所示。同时对导水断裂F1周围采取加大渗透系数处理的办法，结果见表1。
Figure 2. 3D model of the simulated area space
图2. 模拟区空间三维模型
Table 1. Percentage table of target aquifer permeability coefficient in evaluation area
表1. 评价区目标含水层渗透系数分区表
Figure 3. Schematic diagram of the permeability coefficient of the study area
图3. 研究区渗透系数分区示意图
3.2.2. 降雨入渗系数
项目区域多年降雨年平均降水量1351.1 mm。区域降雨入渗系数分布如图4所示。
3.2.3. 弥散系数
开展野外弥散试验并参考经验值在模型中进行适当修正，纵向弥散系数为2 m2/d，横向弥散系数为0.2 m2/d。
3.3. 模型的验证识别
本次地下水位分布基本和地形变化一致，受区内地形地貌、富水岩组分布及构造发育等条件控制，水力坡度从北部山顶各山脊线沿山坡向谷底及溪谷下游逐渐减少，受区域地表水控制最终汇入南家河，这与当地水文地质条件及实地调查结果相符，流场基本合理。本次模拟利用废石场场内及拦挡坝处勘查期间七个钻孔水位观测孔，用以验证其模拟计算水位与实际工程地勘中的实际观测水位之间的匹配性。
模拟区潜水含水层观测水位和计算水位拟合情况如图5所示。
3.4. 模拟因子选择和污染源强设定
本次源强依据废石场废石浸溶试验结果，选取浸溶实验浸出且有质量标准、占标率较高的氟化物作为预测特征因子。氟化物的初始浓度为3.77 mg/L。利用验证后的模型来预测不同时间下各预测因子的浓度。
由于东部废石场目前已经堆存了10年，本次模拟主要考虑有三个时间段对地下水环境的影响：东部废石场堆存10年的污染物模拟、往后继续堆存至服务期满(18年)的污染物模拟、封场后5年(23年)污染物模拟。
3.5. 污染物浓度变化趋势
如图6~8显示了废石场渗滤液中氟化物随水流运移扩散的污染晕10年、18年、23年时浓度分布情况。
Figure 4. Distribution of rainfall infiltration recharge in the study area
图4. 研究区降雨入渗补给分布图
Figure 5. Observed hole water level drilling fit map
图5. 观测孔水位钻孔拟合图
Figure 6. Distribution of fluoride contamination in groundwater during 10 years of storage
图6. 堆存10年时地下水中氟化物污染晕分布图
Figure 7. Distribution of fluoride contamination in groundwater during 18 years of storage
图7. 堆存18年时地下水中氟化物污染晕分布图
Figure 8. Distribution of fluoride contamination in groundwater during 23 years of storage
图8. 堆存23年时地下水中氟化物污染晕分布图
通过模拟结果可以看出，现有东部废石场内废石经雨水林滤产生的淋溶液下渗对下伏潜水含水层已经造成影响，并存在超标现象，随着时间的推移、污染物的扩散，超标面积逐渐增大，同时影响范围逐渐扩大。按照东部废石场堆存10年的时间计算，现有氟化物污染面积为2,016,557 m2，超标面积为2,079,100 m2，氟化物污染物已经进入了万龙山尾矿库中。堆存15年后，污染范围会进一步扩大，氟化物污染面积为2,889,522 m2，超标面积为2,801,635 m2，由于受地形和边界条件的控制，23年时氟化物超标面积变为2,856,687 m2，污染物进一步扩散至万龙山尾矿库中。
4. 总结
1) 运用GMS软件模拟的结果可以得出废石场渗滤液中的污染物随地下水溶质运移情况，通过模型校正调参，其模拟效果较好，反映了该区域地下水的一些基本情况，不仅在理论上是合理的，而且在实际模拟过程也是可行的。
2) 预测结果表明：在设定的污染运移条件下，污染物进入地下水，含水层中污染物浓度逐渐升高，污染范围不断扩大，对场地地下水水质产生污染。根据东部废石场模拟预测分析，东部废石场内污染物将持续扩散至万龙山尾矿库内，对下游的泉点产生影响小。超标范围集中在企业的用地范围内，万龙山尾矿库水旱季全部回用，雨季部分外排，废石场淋滤水对地下水影响可接受。
3) 本研究区的废石属第I类一般工业固体废物，根据《一般工业固体废物贮存、处置场污染控制标准》(GB 18599-2001)并未要求对贮存、处置场进行防渗处理。但《地下水质量标准》(GB/T14848-2017)实施后，由于地下水标准要求严格，第I类一般工业固体废物堆场地下水预测可能会出现超标范围，因此在地下水敏感地区，需要依据查明地下水污染分布特征及迁移转化特征及污染区水文地质条件，采取相应的防渗措施。
NOTES
*通讯作者。