根据研究区水文地质条件可将研究区的水文地质概念模型概化为具有非均质、各向异性三维地下水非稳定渗流系统，概化模型边界范围见图2。模型范围基于ArcGIS平台流域分析结合水文地质条件圈定。模型西部以地表河流沙河为界，全长约5.7 km，模型北部、东部和南部均以地下水分水岭为界，全长约28.85 km，模拟区总面积约32.35 km²，基本构成一个较完整的水文地质单元模型顶面接受降雨补给，地表水体为排泄边界，底面与下伏地层无水量交换。

模型结构通过勘察资料及水文地质条件确定。根据评价区水文地质条件分析，确定评价区边界条件如下：

① 四周边界

北部、东部及南部均由丘陵形成的地表分水岭作为隔水边界，西部以沙河为边界，定为第一类边界(给定水头边界)。

② 上边界为降水补给、蒸发排泄边界。

③ 下边界以微风化基岩作为底部相对隔水边界。模型垂向自上往下依次分为1层：一层为：基岩风化裂隙水，模型空间结构见图3。

图2. 概念模型边界示意图

图3. 概念模型空间结构图

根据前文概化的水文地质概念模型，可用以下地下水渗流偏微分方程及其定解条件来表示：

∂ ∂ x ( K x x ∂ H ∂ x ) + ∂ ∂ y ( K y y ∂ H ∂ y ) + ∂ ∂ z ( K z z ∂ H ∂ z ) + w = μ s ∂ H ∂ t (1)

H ( x , y , z , 0 ) = H 0 , ( x , y , z ) ∈ Ω (2)

K ∂ H ∂ n | S 2 = q ( x , y , z , t ) , ( x , y , z ) ∈ S 2 (3)

H ( x , y , z , t ) = H 1 , ( x , y , z ) ∈ S 1 (4)

式中，Ω：地下水渗流区域，根据概化模型，Ω为6.65 km²；H₀：初始地下水位，H₀的取值后文给出，量纲为L；H₁：指定水位，量纲为L；S₁：第一类边界水头边界；S₂：第二类边界流量边界；μ_s：单位储水系数，量纲为L⁻¹；K_xx，K_yy，K_zz：分别为x、y、z主方向的渗透系数，量纲为LT⁻¹；w：源汇项，包括蒸发，降雨入渗补给，量纲为T⁻¹；q(x，y，z，t)：表示在边界不同位置上不同时间的流量，量纲为L³T⁻¹；

∂ H ∂ n ：表示水力梯度在边界法线上的分量。

结合本次概化的水文地质模型，式(2) 中由于为零流量边界，故q取为0；北部竹皮河等效定义为定水头边界，故式(4)中 H₁取值为竹皮河水头值65 m，南部东宝水库等效定义为定水头边界，故该处H₁为87 m。

将划定的模拟区离散为不规则三角剖分网格。模型平面上共剖分8180个节点，模型总共15,952个单元。根据勘探钻探资料结合评价区地层资料，利用克里金插值法得到各片高程，输入FEFLOW后得到评价区三维地质模型如图3，其中结点数18,526个，有限单元数18,077个。

根据抽水试验、渗水试验等水文地质试验取得本次模拟的水文地质参数，校验模型需对模型进行上百次运算，参数校验除基于勘察资料外，主要采用“试错法”调整。最终拟合结果的合理性依据：研究区流场形态、模型收敛性、水位观测孔拟合数量。根据模型识别验证，参数取值见表1。

表1. 评价区水文地质初始参数取值表

本次预测评价中初始流场采取的技术方法是将模拟区参数分区及初始参数取值表输入模型，经过稳定流计算得到评价区稳定流条件下的天然流场，然后根据实际观测水位对天然流场进行参数校正，得到校正后的地下水初始流场及水文地质参数(图4)，其中图4展示了校正后的地下水初始流场。通过图4可以看到，拟合点基本均匀分布在标准线附近，反应了模拟结果与实际测量值之间总体趋势是一致的，所有值基本都分布在置信区间(5 m)内，可以认为模拟值与实际值拟合情况较好，图5可以看到初始流场的模拟值与实际观测值在总体趋势上趋于一致，拟合效果较好。

通过上述拟合数据对比，可以说明本次建立的地下水稳定流模型基本符合评价区实际水文地质条件，基本反映了地下水系统的流场特征，故可以利用该模型得到的水头值作为非稳定流的初始流场并以此为基础对建设区地下水环境影响进行预测评。

由于水动力弥散尺度效应的存在为模拟和预测地下水中溶质在介质中的运移规律带来了困难，本次弥散度的选取在弥散试验的基础上结合了Geihar [11] 等(1992)对世界范围内所收集的59个大区域弥散资料进行的整理分析成果综合确定。最终确定的溶质运移模型参数如表2所示。

情景类型：非正常工况(事故条件下)、人工防渗措施失效。泄漏源强类型：间断源强。垃圾渗滤液收集室位于厂区中部垃圾卸料大厅下方，其长度与垃圾卸料大厅相的长度等，总容积1180 m³，按高度为3 m计算。渗滤液收集室总面积为393 m²，事故条件人工防渗失效，则天然渗透系数取1.015 × 10⁻⁴ cm/s，按30年模拟。泄漏点：厂区中部垃圾渗滤液收集室。泄漏面积：393 m²。泄漏时间：假设垃圾渗滤液收集室渗漏发生事故导致污染物渗入地下，从发生泄漏到泄漏污染物处理完毕不再发生污染的时间长为5天。泄漏量：渗入地下体积172.3 m³。泄漏浓度：污水中NH₃-N：15,000 mg/L。

为了研究污染物随时间迁移规律，本次预测年限设定为发生泄露100天、1000天和30年，本次评价的特征污染物参照《GB /T14848-93》中相应高锰酸盐指数的相应质量等级标准浓度值(0.2 mg/L)，预测过程中浓度超过上述标准值的区域即为超标区域。

由图6可知，在30年模拟期中，事故工况且人工防渗失效下污染物下渗后直接进入基岩风化裂隙水中，受地下水流向控制逐步向北部迁移扩散，超标污染晕在平面上逐年扩散，但并未超出厂区边界。

图4. 初始流场水位拟合散点分布图

图5. 初始流场水位拟合折线图

表2. 溶质运移模型参数表

图6展示了模型运行100天、1000天和30年三个时段下地下水中污染物的迁移扩散情况。表3针对三个典型时间段，统计了污染晕的运移距离、最高浓度、污染面积和垂向迁移距离等模拟结果。

在平面上污染晕整体运移方向在流向控制下向北部迁移，污染晕最高浓度逐步降低，面积呈先扩大后缩小的趋势。三个时间段污染晕最高浓度分别约为915.2 mg/L、335.6 mg/L、77.6 mg/L，从污染区边界算起，其最大迁移距离分别约为25.4 m、31.4 m、66.9 m，污染面积约为8800.0 m²、9839.5 m²、6694.9 m²，在30年的模拟期内污染物未迁移出厂区边界。

图6. 污染物迁移结果图

表3. 溶质运移结果表

在垂向上污染晕逐步向下迁移，三个时段中最大垂向迁移距离分别约为63.0 m、35.0 m、7.0 m。

综上所述，事故工况且防渗失效下，30年运行期间污染物对基岩风化裂隙水含水层造成了一定的污染，但污染物不会迁移出厂边界，对下游村庄和沙河基本无影响。