本发明涉及水环境保护领域，更具体地涉及一种地下水水质综合预警方法，为地下水污染防控提供科学依据。
背景技术：
地下水污染预警基于预警理论，通过明确水质污染的现状、发展程度以及危害程度明确预警的级别，并且及时给出有效的警报，警示人们对地下水的开采利用所带来的影响，进一步规范约束人类的行为，同时为政府部门制定保护地下水资源的管理方针奠定基础，从而通过有效的防治措施控制地下水的演变方向。
目前国外通过预警技术方法解决地下水污染问题的应用尚不多见，对于地下水污染预警多着眼于水质状态预警，即通过进行水质现状监测实现预警的目标。这种方法缺乏一定的预防性，更加趋向于鉴别和控制水质恶化的结果。我国地下水污染预警主要进行的是趋势预警，发展多集中在预测模型方面。总的来说，目前地下水污染预警的技术方法还不够完善，有待于深入的探索。
技术实现要素：
基于上述预警研究中存在的问题，本发明的目的在于提供一种综合考虑地下水水质现状预警和水质趋势预警的预警方法。
为实现上述目的，本发明提供了一种地下水水质综合预警方法，包括：
(1)监测待预警区域地下水监测井中现状条件下的各污染物浓度，采用改进的内梅罗综合评价指数法计算水质综合评价指数；通过所述水质综合评价指数与水质状态预警等级划分标准的对比，确定水质状态预警级别。
(2)通过地下水水流和溶质运移模拟软件进行各污染物的迁移模拟，得出各污染物的浓度随时间的变化，据此计算水质变化速率；通过水质变化速率与水质趋势预警等级划分标准的对比，确定研究区水质趋势预警级别；
(3)综合考虑所述待预警区域的地下水水质状态预警警度和地下水水质趋势预警警度，采用组合矩阵法建立地下水水质预警多因素综合评价技术方法，确定所述待预警区域的地下水水质的预警警度。
基于上述技术方案可知，本发明的预警方法具有以下优点：
(1)本发明方法综合考虑了地下水水质现状情况、水质趋势变化，同时考虑了地下水中多种组分，建立了地下水水质综合多因素预警评价方法；
(2)本发明采用地下水溶质运移模型进行水质趋势预警，结果准确程度高，符合研究区实际的水文地质情况，特别适用于实际复杂的水文地质条件；
(3)本发明为了解决忽视实测浓度不大但对水质影响极大的评价因子和过分凸显污染指数最大值情况，采用改进的内梅罗综合评价指数法进行地下水水质评价。
附图说明
图1为本发明的地下水水质综合预警方法的流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本发明作进一步的详细说明。
本发明公开了一种地下水水质综合预警方法，包括基于地下水水质现状监测数据采用改进的内梅罗综合评价指数法进行水质状态预警，采用visualmodflow软件建立的水流模型和溶质运移模型进行水质趋势预警，采用矩阵法叠加水质状态预警结果和水质趋势预警结果确定水质综合预警警度。本发明的方法能够有效进行地下水水质综合预警，为政府部门开展地下水污染防治工作提供可靠的科学依据。
更具体地，本发明的地下水水质综合预警技术方法，包括以下步骤：
(1)改进的内梅罗综合评价指数法评价水质
为了评价地下水水质现状和未来不同时间节点上的地下水水质趋势，采用改进的内梅罗综合评价指数法的计算公式和步骤如下：
①各项评价因子的评价
根据地下水水质监测结果中各项污染组分的实测浓度与地下水质量标准，计算各项评价因子的污染指数fi。计算中将《地下水质量标准》(gb14848-93)iii类标准值作为基准。计算公式如下：
式中：ci——第i种评价因子的实测浓度，i＝1，2，...，n；
sij——第i种评价因子的第j类标准浓度，j＝1，2，...，m。
②水质综合评价
按照下述公式进行综合评价分值f的计算：
式中：f——水质综合评价指数；
f′max———n种评价因子分值fi的最大值(修正)。
其中：
式(3)和(4)中：wi——第i种评价因子的权重值；
s最大——第i种评价因子的最大标准值；
si——各评价因子的标准值；
n——评价因子的个数。
式(5)和(6)中：fmax——n种评价因子分值fi的最大值；
fw——k项污染因子分值fi的平均值，k为n项中
③根据f值按照以下对应关系划分地下水质量级别。
当f＜0.8时，所述地下水质量级别为优良；
当f介于0.80-2.50之间时，所述地下水质量级别为良好；
当f介于2.50-4.25之间时，所述地下水质量级别为较好；
当f介于4.25-7.2之间时，所述地下水质量级别为较差；
当f＞7.20时，所述地下水质量级别为极差。
(2)确定水质状态预警级别
通过水质综合评价指数f与水质状态预警等级划分标准的对比确定水质状态预警的级别，具体对应关系为：
当f＜0.8时，所述水质状态预警等级划分为零级，警报颜色为绿色，说明地下水潜水面处未受到污染，地下水适用于各种用途；
当f介于0.80-2.50之间时，所述水质状态预警等级划分为一级，警报颜色为蓝色，说明地下水潜水面处受到污染较轻，地下水主要适用于集中式饮用水水源和工、农业用水；
当f介于2.50-4.25之间时，所述水质状态预警等级划分为二级，警报颜色为黄色，说明地下水潜水面处受到污染中等，地下水适用于工、农业用水，适当处理后可作生活饮用水；
当f介于4.25-7.2之间时，所述水质状态预警等级划分为三级警报颜色为橙色，说明地下水潜水面处受到污染较重，地下水不宜饮用，其他用水可根据使用目的选用；
当f＞7.20时，所述水质状态预警等级划分为四级，警报颜色为红色，说明地下水潜水面处受到污染非常严重，地下水不宜使用。
(3)确定污染物浓度随时间的变化趋势
采用地下水水质和溶质运移模拟软件visualmodflow进行污染物在地下水中的迁移模拟，得出污染物浓度随时间的变化，具体步骤包括：
①收集获取研究区水文地质资料(含水层厚度、含水层顶底板等高线、含水层等水位线、渗透系数、给水度、孔隙度等)和污染源资料(污染源位置、特征污染物种类、排放量和浓度、排放去向等)；
②建立水文地质概念模型，包括研究区域范围、概化研究区域水文地质条件、边界条件、水文地质参数、源汇项、水文地质内部结构以及地下水流场等；
③建立地下水水流模型；
地下水水流模型的数学表达式如下：
式中：kx、ky、kz-渗透系数在x、y和z方向上的分量，m/d；
h-水头，m；
w-汇源项，l/d；
ss-孔隙介质的贮水率，1/m；
t-时间，d。
水流模型的初始条件方程为：
h(x，y，z，t)＝h0(x，y，z)；x，y，z∈ω式(8)
式中：h0-已知初始水头；
ω-研究区域的范围。
水流模型有三种类型的边界条件，即为水头边界、流量边界以及混合边界，各类边界的基本方程如下所示。
水头边界基本方程：
式中：γ1-研究区域的第一类边界(水头边界)；
h1(x，y，z，t)-沿着γ1的已知水头值，会随着时间而变。
流量边界基本方程：
式中：k-含水层的渗透系数，m/d；
γ2-研究区域的第二类边界(流量边界)；
q(x，y，z)-沿着γ2的已知流量值。
混合边界基本方程：
式中：γ3-研究区域的第三类边界(混合边界)；
q(x，y，z)-沿着γ3的已知流量值，m3/d；
k’-河床的渗透系数，m/d；
h0-河流水位，m；
b’-河床厚度，m。
④建立地下水溶质运移模型；
地下水溶质运移模型的数学表达式如下：
式中：c-地下水污染物浓度，kg/m3；
θ-介质的孔隙率，无单位；
t-时间，d；
xi、xj-沿直角坐标轴的距离，m；
d-水动力弥散系数，m2/d；
v-地下水实际流速，m/d；
q-单位体积含水层流量，l/d；
cs-源或汇水流中污染物浓度，kg/m3；
∑r-化学反应项，kg/(m3·d)。
⑤联合运行水流模型和溶质运移模型，预测污染物在预测时段内不同时间节点上的污染物浓度。
(4)计算水质变化速率
根据不同时间节点上的水质综合评价指数计算水质变化速率，具体公式见：
式中：r-水质变化速率；
e(t)-未来某一时刻的水质综合评价指数，水质综合评价指数通过改进的内梅罗综合评价指数法求得；
e(t)-已知某一时刻的水质综合评价指数；
δt-时段个数，根据经验数值，以每2年为一个时段。
(5)确定水质趋势预警级别
通过水质变化速率r与水质趋势预警等级划分标准的对比，确定研究区地下水水质趋势预警级别，具体对应关系如下：
当r≤0时，水质趋势预警等级划分为零级，警报颜色为绿色，说明地下水水质稳定或有好转趋势；
当0＜r≤0.1时，水质趋势预警等级划分为一级，警报颜色为蓝色，说明地下水水质轻度恶化趋势；
当0.1＜r≤0.25时，水质趋势预警等级划分为二级，警报颜色为黄色，说明地下水水质中度恶化趋势；
当0.25＜r≤0.4时，水质趋势预警等级划分为三级，警报颜色为橙色，说明地下水水质重度恶化趋势；
当r＞0.4时，水质趋势预警等级划分为四级，警报颜色为红色，说明地下水水质有突变趋势。
(6)确定水质预警综合级别
采用组合矩阵法确定地下水水质预警综合警度，具体方法见下表。如，当水质趋势预警级别为一级，水质状态预警级别为二级时，水质综合预警等级为一级。
实施案例
下面结合说明书附图和实施案例，对本发明作进一步具体阐释说明。
某畜禽养殖场可能对周边第四系孔隙潜水含水层水质造成一定程度的污染，因此以地下水流下游的水井g1作为评价对象，通过本发明的地下水水综合预警方法确定水井的预警综合等级。经过地下水监测得知，g1中cod为13.1mg/l，氨氮为0.068mg/l，硝酸盐氮为0.125mg/l。
(1)改进的内梅罗综合评价指数法评价地下水水质
①各项污染组分的评价
将《地下水质量标准》(gb14848-93)iii类标准值作为基准，根据式(1)计算，fcod＝0.655，f氨氮＝0.34，f硝酸盐氮＝0.125。
②水质综合评价
根据式(4)计算，cod、氨氮和硝酸盐氮的权重w分别为0.0098、0.9804、0.0098。根据式(3)计算，
(2)确定水质状态预警级别
f＝0.2450，通过水质综合评价指数与水质状态预警等级划分标准的对比确定水质状态预警的级别可知，g1预警级别为零级，警报颜色为绿色。
(3)确定污染物浓度随时间的变化趋势
利用visualmodflow软件，进行地下水中cod、氨氮以及硝酸盐氮的运移模拟，基于式(7)-式(12)计算预测时间为2a时g1中cod、氨氮以及硝酸盐氮的浓度分别为15.52mg/l、0.154mg/l、0.223mg/l，预测时间为4a时g1中cod、氨氮以及硝酸盐氮的浓度分别为18.02mg/l、0.218mg/l、0.358mg/l。
(4)计算水质变化速率
预测时间为2a时，g1的fcod＝0.776，f氨氮＝0.77，f硝酸盐氮＝0.0112，
根据式(13)计算，r＝0.2416。
(5)确定水质趋势预警级别
根据水质变化速率与预警级别之间的关系表确定g1的水质趋势预警级别为二级，警报颜色为黄色，说明潜水含水层有中度恶化趋势。
(6)确定水质预警综合级别
水质现状预警级别为零级，水质趋势预警级别为二级，根据组合矩阵法确定地下水水质预警综合警度，水质预警综合级别为一级。警报颜色为蓝色，需要发布蓝色预警。
该水井地下水水质现状综合指数为0.2450，水质现状优良，4年后地下水水质综合评价指数为0.8115，水质有中度恶化趋势。可见，该水井地下水面临一定程度的污染风险，污染程度一般，须关注地下水水质的变化情况，并制定相应的污染防治措施以应对养殖基地污水渗漏带来的污染。另外，该井中氨氮浓度4年后达到0.218mg/l，超过水质质量iii类标准，需要提前采取相应的污染防治措施。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。