■ 河北建設(shè)勘察研究院有限公司 吳立彬

地下水是重要的水資源，對(duì)人類(lèi)生存至關(guān)重要。在礦產(chǎn)資源開(kāi)采過(guò)程中，重金屬元素和有害離子等污染物質(zhì)會(huì)遷移到水文地質(zhì)環(huán)境，影響地下水質(zhì)量，造成地下水污染[1]。如果不能采取有效措施對(duì)水質(zhì)污染提前預(yù)警，勢(shì)必會(huì)嚴(yán)重威脅人類(lèi)健康，因此需要對(duì)水文地質(zhì)勘探地下水污染進(jìn)行預(yù)測(cè)。

學(xué)術(shù)界對(duì)于污染預(yù)測(cè)模型已經(jīng)積累一定研究成果。符東等[2]優(yōu)化篩選評(píng)價(jià)因子，構(gòu)建BP神經(jīng)網(wǎng)絡(luò)，并對(duì)主要污染因子進(jìn)行預(yù)測(cè)，得到了監(jiān)測(cè)斷面的水質(zhì)標(biāo)準(zhǔn)；侯澤宇等[3]利用學(xué)習(xí)機(jī)替代模型，解決了含水層參數(shù)的同步反演識(shí)別問(wèn)題，大幅度減少了模擬的計(jì)算負(fù)荷，縮短了識(shí)別時(shí)間。但上述模型在預(yù)測(cè)準(zhǔn)確度上還存在優(yōu)化空間?，F(xiàn)有預(yù)測(cè)模型以地下水污染的靜態(tài)分析為依據(jù)，忽略了各影響因子之間的相互聯(lián)系。而模糊綜合評(píng)價(jià)克服了各因子之間相互獨(dú)立的局限性，可以較為全面地預(yù)測(cè)地下水污染情況，評(píng)估水資源質(zhì)量。本研究基于模糊綜合評(píng)價(jià)，構(gòu)建水文地質(zhì)勘探地下水污染預(yù)測(cè)模型，以期提升預(yù)測(cè)質(zhì)量。

1.基于模糊綜合評(píng)價(jià)的水文地質(zhì)勘探地下水污染預(yù)測(cè)模型

1.1 選取參考指標(biāo)計(jì)算地下水污染物通量

首先選取影響地下水質(zhì)的參考指標(biāo)。選取指標(biāo)時(shí)，應(yīng)綜合地下水污染問(wèn)題建立綜合性的評(píng)價(jià)體系。選取的參考指標(biāo)要具有通用的概念，能夠與其他環(huán)境或區(qū)域的水文地質(zhì)勘察情況進(jìn)行橫向比較，在時(shí)間序列上也具有縱向可比性[4]，其中包含變化與非變化指標(biāo)，能夠全面反映地下水質(zhì)變化情況。參考指標(biāo)需要擁有側(cè)重點(diǎn)和層次性，使地下水污染預(yù)測(cè)更具科學(xué)性。綜合考慮技術(shù)支持能力和研究工作量，參考相關(guān)質(zhì)量監(jiān)測(cè)標(biāo)準(zhǔn)，將擬測(cè)定指標(biāo)分為基本指標(biāo)和特征指標(biāo)，構(gòu)建參考指標(biāo)體系如下：設(shè)定地下水質(zhì)量為目標(biāo)層，基本指標(biāo)和特征指標(biāo)為系統(tǒng)層，二者分別對(duì)應(yīng)具體指標(biāo)層?；局笜?biāo)對(duì)應(yīng)的指標(biāo)層包括PH值、氧化還原電位、總硬度、溶解性總固體、揮發(fā)酚、LAS、電導(dǎo)率、MnO4-指數(shù)、溶解氧、重金屬離子濃度、硝酸鹽、氮氰化物等；特征指標(biāo)對(duì)應(yīng)的指標(biāo)層為有色金屬礦，包括銅、鋅、鉛、汞、氰化物和石油類(lèi)等。在選取地下水污染物指標(biāo)的基礎(chǔ)上，計(jì)算各污染物通量[5]。污染物通量是特定時(shí)間內(nèi)通過(guò)監(jiān)測(cè)斷面的污染物總量，不僅受到點(diǎn)源污染影響，還與地下水流量等因素的影響，其計(jì)算結(jié)果直接反映各污染源對(duì)水體的貢獻(xiàn)量，展示污染分布特征，是地下水污染預(yù)測(cè)的先決條件[6]。計(jì)算方法包括斷面瞬時(shí)通量與時(shí)段通量計(jì)算，本文主要利用污染物時(shí)段通量進(jìn)行預(yù)測(cè)。污染物時(shí)段通量的計(jì)算公式為：

式（1）中，F(xiàn)表示污染物總量；t表示預(yù)測(cè)時(shí)段；D（t）表示瞬時(shí)流量；P（t）表示瞬時(shí)濃度。此次研究使用時(shí)段平均濃度與水量的乘積，計(jì)算監(jiān)測(cè)斷面的污染物通量。利用半變異函數(shù)計(jì)算污染物通量在空間上的連續(xù)性變化，具體計(jì)算過(guò)程為：

式（2）中，v（d）表示污染物半變異值；d表示空間分割距離；n（d）表示距離d的點(diǎn)對(duì)數(shù)；S（xi）表示空間變量?；谏鲜鲇?jì)算方法，得到污染物空間分布特征。以氯化物和氨氮為典型無(wú)機(jī)污染物；石油類(lèi)和苯為典型有機(jī)污染物，對(duì)典型污染物進(jìn)行異常值處理，可得到最優(yōu)擬合模型，如圖1所示。

圖1 典型污染物的最優(yōu)擬合

根據(jù)污染物異常值的最優(yōu)擬合結(jié)果，得到空間分布差異特征，以便對(duì)研究區(qū)域污染物進(jìn)行差異性分析。通過(guò)選取影響地下水質(zhì)的參考指標(biāo)，計(jì)算污染物通量，確定污染物濃度差異與分布范圍，并以此為基礎(chǔ)進(jìn)行下一階段的模型設(shè)計(jì)。

1.2 利用模糊綜合評(píng)價(jià)建立指標(biāo)關(guān)系矩陣

現(xiàn)有模型局限于靜態(tài)分析，考慮到對(duì)地質(zhì)勘探地下水長(zhǎng)期動(dòng)態(tài)結(jié)果進(jìn)行評(píng)價(jià)，因此利用模糊綜合評(píng)價(jià)法對(duì)研究斷面的水質(zhì)進(jìn)行評(píng)價(jià)[7]。模糊綜合評(píng)價(jià)法將主觀與客觀因素相結(jié)合，通過(guò)設(shè)定的指標(biāo)體系和賦值權(quán)重，計(jì)算各因子對(duì)指標(biāo)的隸屬度，分析結(jié)果向量，進(jìn)而得到區(qū)域地下水質(zhì)量等級(jí)及排序[8]。首先建立因子集和評(píng)價(jià)集，以參考指標(biāo)建立的評(píng)價(jià)因子集U可表示為：

式（3）中，um表示地下水質(zhì)量指標(biāo)；m表示選取指標(biāo)總個(gè)數(shù)。根據(jù)地下水質(zhì)量標(biāo)準(zhǔn)，確定水質(zhì)標(biāo)準(zhǔn)評(píng)價(jià)集V為：

式（10）中，axy表示第x個(gè)指標(biāo)因子第y類(lèi)標(biāo)準(zhǔn)的隸屬度。正向指標(biāo)評(píng)價(jià)集的隸屬度函數(shù)計(jì)算公式與上述公式相反。由于各指標(biāo)因子對(duì)地下水質(zhì)量的影響程度不同，需要對(duì)各個(gè)因子賦予不同權(quán)重[9]。權(quán)重反映指標(biāo)因子的重要程度，本文結(jié)合層次分析與熵權(quán)法，對(duì)權(quán)重進(jìn)行賦值。層次分析中，對(duì)選定的評(píng)價(jià)指標(biāo)進(jìn)行數(shù)據(jù)處理，對(duì)于正向指標(biāo)存在：

式（11）中，b表示指標(biāo)標(biāo)度值；mc表示該指標(biāo)實(shí)測(cè)濃度；表示濃度標(biāo)準(zhǔn)平均值。對(duì)于逆向指標(biāo)，指標(biāo)標(biāo)度值為正向評(píng)價(jià)因子的倒數(shù)[10]。構(gòu)造判斷矩陣，得到特征向量，歸一化處理后得到權(quán)重系數(shù)。熵權(quán)法中，利用評(píng)價(jià)指標(biāo)實(shí)測(cè)值構(gòu)造判斷矩陣，進(jìn)行歸一化處理，對(duì)于正向指標(biāo)存在：

式（12）中，h表示指標(biāo)歸一化值；samax與samin分別表示指標(biāo)中最滿(mǎn)意和不滿(mǎn)意的值。對(duì)于逆向指標(biāo)存在：

計(jì)算各指標(biāo)熵值E，表示為：

式（15）中，Ej表示各指標(biāo)熵值；g表示指標(biāo)個(gè)數(shù)。計(jì)算層次分析與熵權(quán)法得到權(quán)重的算術(shù)平均值，即為評(píng)價(jià)指標(biāo)的權(quán)重，構(gòu)建因子權(quán)重集W，與模糊矩陣A相乘，得到綜合指標(biāo)。綜合指標(biāo)反映全部隸屬度信息，更能真實(shí)反映實(shí)際情況[11]。基于上述模糊綜合評(píng)價(jià)過(guò)程，在各指標(biāo)間建立有效連接，以此實(shí)現(xiàn)對(duì)地下水質(zhì)的綜合評(píng)價(jià)。

1.3 構(gòu)建地下水污染預(yù)測(cè)模型

為實(shí)現(xiàn)污染物運(yùn)移數(shù)值模擬，建立地下水污染預(yù)測(cè)模型，通過(guò)模擬結(jié)果預(yù)測(cè)和分析污染物的濃度變化和運(yùn)移路徑，進(jìn)而提出有效控制措施，減少水質(zhì)污染情況。首先根據(jù)研究區(qū)地質(zhì)條件，確定模擬邊界。水平邊界包括第一類(lèi)定水頭邊界和第二類(lèi)定流量邊界[12]。垂直邊界主要為向上的潛水面和底部隔水邊界[13]。根據(jù)地質(zhì)條件，建立地下水流運(yùn)動(dòng)控制方程，具體過(guò)程如下：

式（16）中，x，y，z表示三維坐標(biāo)軸；λ1，λ2，λ3，分別表示x，y，z方向的滲透系數(shù)；WL表示地下水位；o表示源匯項(xiàng)；表示含水層給水度；t表示時(shí)間?？紤]對(duì)流和彌散作用對(duì)污染物遷移預(yù)測(cè)的影響，建立污染物遷移控制方程：

式（17）中，k表示地層有效孔隙度；c表示污染物濃度；表示水動(dòng)力彌散系數(shù)；表示實(shí)際流速。由于該模型為三維非均質(zhì)模型，且均有各向非穩(wěn)定性，需要確定初始和邊界的定解條件[14]。模擬地下水流動(dòng)的定解條件，要確定初始水位和邊界條件；模擬污染物遷移的定解條件，要確定該區(qū)域初始濃度、濃度邊界分布和連續(xù)注入濃度[15]。通過(guò)上述過(guò)程的數(shù)值模擬，完成對(duì)污染物濃度變化和運(yùn)移路徑的預(yù)測(cè)，至此在模擬綜合評(píng)價(jià)方法的輔助下，實(shí)現(xiàn)對(duì)水文地質(zhì)勘探地下水污染預(yù)測(cè)模型的構(gòu)建。

2.實(shí)驗(yàn)研究

2.1 實(shí)驗(yàn)準(zhǔn)備

將某水文地質(zhì)勘探區(qū)域作為實(shí)驗(yàn)測(cè)試對(duì)象，對(duì)該區(qū)域的地下水污染情況進(jìn)行預(yù)測(cè)。該區(qū)域從上到下分別為人工堆積層、堆積物和天然地層。人工堆積層包括細(xì)砂和素填土，包含塊石和碎石土；堆積物包括細(xì)砂、粉砂和粉質(zhì)黏土；天然地層為強(qiáng)風(fēng)化和全風(fēng)化的花崗巖。該區(qū)域地下水的主要補(bǔ)給來(lái)源為大氣降水，方式以垂向滲入為主；水平徑流緩慢，通過(guò)排水管道、壩體排水和蒸發(fā)進(jìn)行地下水排泄。選取5個(gè)測(cè)試點(diǎn)，對(duì)地下水污染情況進(jìn)行監(jiān)測(cè)。根據(jù)各監(jiān)測(cè)井污染物指標(biāo)的組合賦權(quán)權(quán)重，計(jì)算各監(jiān)測(cè)井評(píng)價(jià)指標(biāo)的污染權(quán)重，具體結(jié)果見(jiàn)表1。

根據(jù)表1結(jié)果，可得到各測(cè)試點(diǎn)的主要污染指標(biāo)及其權(quán)重，利用本文構(gòu)建的污染預(yù)測(cè)模型確定各測(cè)試點(diǎn)的水質(zhì)污染預(yù)測(cè)等級(jí)。

表1 各測(cè)試點(diǎn)評(píng)價(jià)指標(biāo)污染權(quán)重

2.2 地下水污染預(yù)測(cè)效果評(píng)估

對(duì)該研究區(qū)域的5個(gè)測(cè)試點(diǎn)進(jìn)行連續(xù)一年的地下水污染預(yù)測(cè)，并記錄實(shí)際污染物濃度值，測(cè)定實(shí)際水質(zhì)污染情況。將此次構(gòu)建模型與傳統(tǒng)模型的預(yù)測(cè)結(jié)果進(jìn)行對(duì)比，評(píng)估構(gòu)建的預(yù)測(cè)模型的實(shí)際應(yīng)用效果。利用模型得到的水質(zhì)評(píng)價(jià)等級(jí)，衡量預(yù)測(cè)效果，具體結(jié)果見(jiàn)圖2。

圖2 不同模型的水質(zhì)預(yù)測(cè)等級(jí)對(duì)比

根據(jù)圖2預(yù)測(cè)對(duì)比結(jié)果，本文構(gòu)建的預(yù)測(cè)模型得到的水質(zhì)評(píng)價(jià)等級(jí)與實(shí)際測(cè)定結(jié)果比較接近，而現(xiàn)有預(yù)測(cè)模型的預(yù)測(cè)結(jié)果與實(shí)際污染情況出現(xiàn)較大的偏差，說(shuō)明本文構(gòu)建的模型能夠可靠反映地下水質(zhì)狀況，與實(shí)際污染狀況相符合，比現(xiàn)有模型更具實(shí)用性，預(yù)測(cè)結(jié)果更加合理，驗(yàn)證了本文模型的可靠性。在地下水質(zhì)等級(jí)評(píng)估的基礎(chǔ)上，進(jìn)一步對(duì)污染物濃度進(jìn)行預(yù)測(cè)，計(jì)算實(shí)際測(cè)量值與模型預(yù)測(cè)值，結(jié)果見(jiàn)表2。

表2 地下水污染物濃度實(shí)際測(cè)量值與模型預(yù)測(cè)值

計(jì)算污染物濃度實(shí)際值與預(yù)測(cè)值的相對(duì)誤差，衡量預(yù)測(cè)模型精度，結(jié)果如圖3所示。

根據(jù)圖3對(duì)比結(jié)果，本文構(gòu)建模型的地下水污染物濃度預(yù)測(cè)誤差在3%—5%范圍內(nèi)，相比較現(xiàn)有預(yù)測(cè)模型而言，預(yù)測(cè)誤差降低了8%—10%，說(shuō)明預(yù)測(cè)準(zhǔn)確度得到提升。因此，在實(shí)際應(yīng)用中，與實(shí)際測(cè)量數(shù)據(jù)的結(jié)果相符合，預(yù)測(cè)效果更好 。

圖3 不同模型的污染物濃度預(yù)測(cè)誤差

3.結(jié)語(yǔ)

隨著礦產(chǎn)資源的開(kāi)采，有害元素遷移到水文地質(zhì)環(huán)境，地下水污染程度嚴(yán)重。現(xiàn)有污染預(yù)測(cè)模型局限于靜態(tài)分析，忽略各影響因子的相互聯(lián)系，導(dǎo)致預(yù)測(cè)準(zhǔn)確度較低。本文基于模糊綜合評(píng)價(jià)構(gòu)建水文地質(zhì)勘探地下水污染預(yù)測(cè)模型，選取參考指標(biāo)，計(jì)算地下水污染物通量，確定濃度差異；利用模糊綜合評(píng)價(jià)，建立指標(biāo)關(guān)系矩陣，得到質(zhì)量等級(jí)及排序，加強(qiáng)因子關(guān)聯(lián)性；構(gòu)建地下水污染預(yù)測(cè)模型，完成污染物運(yùn)移數(shù)值模擬。實(shí)驗(yàn)結(jié)果表明，構(gòu)建的預(yù)測(cè)模型能夠可靠反映地下水質(zhì)狀況，污染物濃度預(yù)測(cè)誤差在3%—5%范圍內(nèi)，與現(xiàn)有模型相比，準(zhǔn)確度提高8%—10%，更有利于預(yù)測(cè)地下水環(huán)境狀況。然而，由于水文地質(zhì)勘探環(huán)境較為復(fù)雜，對(duì)預(yù)測(cè)模型模擬效果存在干擾，此次建立的模型還存在不足之處。后續(xù)研究可對(duì)污染物展開(kāi)溯源分析，進(jìn)一步優(yōu)化污染預(yù)測(cè)結(jié)果，提高模型應(yīng)用效果。