[发明专利]一种区域浅层地下水特殊脆弱性的评价方法

权利要求书

1.一种地下水污染防控方法，其特征在于，包括：

获取区域浅层地下水特殊脆弱性的评价指标，所述评价指标包括净补给量、土壤介质、包气带介质、地下水水流速度和土地利用类型；

对各评价指标进行加权求和得到地下水特殊脆弱性指数；

根据地下水特殊脆弱性指数划分地下水特殊脆弱性等级，得到研究区地下水脆弱性评价分区图；

根据研究区地下水脆弱性评价分区图进行地下水污染防控；

其中，所述获取区域浅层地下水特殊脆弱性的评价指标，所述评价指标包括净补给量、土壤介质、包气带介质、地下水水流速度和土地利用类型，包括：

对区内浅层地下水的水均衡条件分析，利用下式计算净补给量：

Q_净补＝Q_降+Q_侧补+Q_河+Q_灌-Q_蒸-Q_开-Q_侧排-Q_基流

其中：Q_降表示降水入渗补给量；Q_侧补表示地下水侧向补给量；Q_河表示河流入渗补给量；Q_灌表示灌溉水入渗补给量；Q_蒸表示蒸发排泄量；Q_开表示人工开采排泄量；Q_侧排表示地下水侧向排泄量；Q_基流表示河流排泄地下水水量；

利用下式计算降水入渗补给量：

Q_降＝αP_降F_总/10

其中，α表示大气降水入渗系数；P_降表示大气降水量；F_总表示评价区面积；

利用下式计算灌溉水回渗补给量：

Q_水田＝βYF_灌/100+βP_灌F_灌/10

Q_菜田＝βYF_灌/100

其中，Q_水田表示水田回渗量；Q_菜田表示菜田回渗量；β表示灌溉水回渗系数；Y表示灌溉定额；F_灌表示灌溉面积；P_灌表示灌溉期大气降水利用量；

利用下式计算侧向径流补给/排泄量：

Q_侧补/Q_侧排＝365BKiH/10⁴

其中，B表示计算断面宽度；K表示渗透系数；i表示水力梯度；H表示含水层厚度；

利用下式计算河流入渗补给量：

Q_河＝L×γ×t

其中，Q_河表示河流入渗补给量；L表示河流补给段长度；γ表示河流入渗补给强度；t表示河流输水时间；

利用下式计算蒸发排泄量：

Q_蒸＝F_蒸E(1-h/l)²

其中，Q_蒸表示蒸发排泄量；E表示蒸发皿观测年蒸发量；F_蒸表示蒸发面积；h表示水位埋深；l表示潜水蒸发极限深度；

利用下式计算河流排泄地下水水量：

Q_基流＝M·F_河

其中，Q_基流表示河流排泄地下水水量；M表示河流平均基流模数；F_河表示流域面积；

其中，所述对各评价指标进行加权求和得到地下水特殊脆弱性指数，包括:

利用下式对各评价指标进行加权求：

RI＝R_WR_r+S_WS_r+I_WI_r+V_WV_r+L_WL_r

其中，RI为地下水特殊脆弱性指数；R为净补给量；S为土壤介质；I为包气带介质；V为地下水水流速度；L为土地利用类型；

其中，所述根据地下水特殊脆弱性指数划分地下水特殊脆弱性等级，得到研究区地下水脆弱性评价分区图，包括：

选取栅格单元为基本的评价单元，利用下式计算适宜栅格单元大小：

G_s＝7.49+0.0006S-2.0*10^-9S²+2.9*10^-15S³

式中，Gs是适宜栅格大小；S为原始等高线数据精度的分母；

基于ArcGIS栅格空间分析对评价指标的评分分区图按各个评价指标的相对权重值进行图层叠加，计算研究区内每个栅格单元上的地下水特殊脆弱性指数，根据地下水特殊脆弱性指数取值按照几何间距法将研究区划分为5个脆弱性分区，最终得到地下水特殊脆弱性分区图。