劉志濤，周群道，楊建華

(山東省魯北地質(zhì)工程勘察院，山東 德州　253072)

?

地下水溶質(zhì)運移數(shù)值法和解析法預(yù)測結(jié)果對比分析
——以沾化電廠為例

劉志濤，周群道，楊建華

(山東省魯北地質(zhì)工程勘察院，山東 德州253072)

摘要：數(shù)值法和解析法是當(dāng)前解決地下水流和溶質(zhì)運移問題最常用的兩種方法，雖然數(shù)值法以其廣泛的適用性和較高的仿真性等優(yōu)點取得了越來越普遍的應(yīng)用，但解析法也以其簡單易用等特點一直成為首選方法之一。該文以沾化電廠地下水環(huán)評項目為例，分別采用兩種方法對污染物在長期和短期滲漏情況下對地下水所造成的污染情況進行了預(yù)測，并在對兩種方法預(yù)測結(jié)果的差異進行分析的基礎(chǔ)上，論述了解析法在該區(qū)的適宜性，為解析法在該區(qū)的應(yīng)用提供了參考。

關(guān)鍵詞：地下水；溶質(zhì)運移；數(shù)值法；解析法

引文格式：劉志濤，周群道，楊建華.地下水溶質(zhì)運移數(shù)值法和解析法預(yù)測結(jié)果對比分析——以沾化電廠為例[J].山東國土資源，2016,32(7)：78-82. LIU Zhitao, ZHOU Qundao, YANG Jianhua. Comparative Analysis on Prediction Results by Using Numerical Method and Analytical Method of Solute Transport in Groundwater——Taking Zhanhua Power Plant As an Example[J].Shandong Land and Resources, 2016,32(7)：78-82.

0引言

數(shù)值法和解析法是目前進行地下水水流運動和溶質(zhì)運移預(yù)測最常用的方法，尤其是隨著計算機科學(xué)技術(shù)的飛速發(fā)展，數(shù)值法在很多領(lǐng)域都得到了越來越廣泛的應(yīng)用，并使得許多過去難以解決的復(fù)雜地下水運動問題在一定程度上得到了解決，同時也提高了對相關(guān)領(lǐng)域的研究水平[1]。但建立地下水?dāng)?shù)值模型對水文地質(zhì)資料的精度和深度要求比解析法要高很多，且由于模型識別驗證需要一定的時間，使得項目持續(xù)時間較長。

在實際地下水環(huán)境影響評價工作中，由于各地水文地質(zhì)工作的研究程度有著很大不同，同時受到項目周期、資金等因素的影響，解析法仍有著廣泛的應(yīng)用。該文以濱州沾化電廠為例，分別用數(shù)值法和解析法進行了溶質(zhì)運移預(yù)測，并對預(yù)測結(jié)果進行了對比分析，以探索解析法在該區(qū)的適宜性。

1工作區(qū)水文地質(zhì)特征

工作區(qū)處于濱州市沾化北部，屬魯西北平原松散巖類水文地質(zhì)區(qū)、濱海沖積海積平原咸水水文地質(zhì)亞區(qū)[2]，區(qū)域地層呈水平層狀分布。區(qū)內(nèi)地下水類型為第四系松散巖類孔隙水，按水力特征又分為淺層孔隙潛水-微承壓水和中、深層承壓水[3]，其中淺層地下水是該次工作主要研究對象。

(1)淺層潛水-微承壓水。工作區(qū)內(nèi)淺層潛水-微承壓水全部為咸水，埋藏深度小于60m，含水層巖性主要是粉砂，共2～3層，累計厚度15m左右，單井涌水量200m3/d左右，富水性較差。地下水埋深1.5～4.3m，由南向北遞增，年變化幅度一般在1.5m以內(nèi)，地下水水位動態(tài)主要受氣象因素的影響，水位動態(tài)曲線與大氣降水的分配過程基本一致。地下水水質(zhì)普遍較差，礦化度均大于10g/L，水化學(xué)類型以Cl-Na型為主。工作區(qū)地下水水位埋深較淺，以垂直方向的補給、排泄為主，大氣降水為其主要補給源，蒸發(fā)式主要排泄途徑。地下水流向自南向北，徑流極為緩慢，水力梯度在0.6×10-3左右。

(2)中深層承壓水。工作區(qū)內(nèi)中深層承壓水皆為咸水，埋深在60～200m范圍內(nèi)，含水層巖性主要為粉細(xì)砂，富水性差，單井涌水量一般小于500m3/d，水化學(xué)類型主要為Cl·SO4-Na型水，礦化度在19g/L左右。該含水層與淺層地下水之間存在較厚的粘土層，隔水性能相對較好，兩含水層之間水力聯(lián)系較弱。

(3)深層承壓水。深層承壓水頂板埋深一般在200m左右，含水層巖性多為粉細(xì)砂，富水性中等，單井涌水量500～1000m3/d，礦化度均大于3g/L，水化學(xué)類型多為氯化物-鈉型。深層承壓水埋藏深度較大，含水層層次多，累計厚度亦較大，與上部含水層之間存在以粘土為主的穩(wěn)定隔水層。

工作區(qū)內(nèi)，淺、中、深各含水層組之間普遍存在厚度較大、分布廣泛、隔水性能相對較好的粉質(zhì)粘土或粘土層，且各含水層的水質(zhì)、水位有著明顯差別，反映了淺、中、深3個含水層組之間地下水水力聯(lián)系微弱[4]。

2數(shù)值法預(yù)測

2.1數(shù)值模型的建立

在平面上，模擬區(qū)范圍主要根據(jù)工作區(qū)水文地質(zhì)條件并結(jié)合地下水流場確定，東部邊界與等水位線基本垂直，南、北部邊界與等水位線大致水平，西部以河流為界，面積約90km2(圖1)。在垂向上，該次模型概化為兩層。第一層主要巖性為粉土、粉質(zhì)粘土，底界面埋深3.4～5.5m，第二層以粉砂為主，底界面埋深19.4～24.5m，是主要含水層。第二層以下有累計厚度超過15m的粉質(zhì)粘土層，其垂向滲透系數(shù)在(1.27～8.11)×10-8cm/s之間，具有良好的隔水性能，因此將其概化為模型的不透水底部邊界。垂向上水量交換主要為大氣降水入滲補給和潛水蒸發(fā)排泄。

圖1　模擬區(qū)范圍圖

地下水水流模型采用非穩(wěn)定三維滲流數(shù)學(xué)模型[5-6]，具體可以描述為：

式中：H—水頭；K—滲透系數(shù)；Ss—貯水率；W—源匯項；G—研究域；?！芯坑虻倪吔?。

溶質(zhì)運移模型采用不考慮吸附、揮發(fā)以及生物化學(xué)反應(yīng)情況下的三維非穩(wěn)定流數(shù)學(xué)模型[7-8]，具體可以描述為：

式中：C為濃度；Dij為水動力彌散系數(shù)張量的分量；W為源匯項；Cw為隨W注入的污染物濃度；Rd為滯留因子；n為孔隙度；ui為水流實際速度分量；其他符號同前。

模擬軟件采用Visaul Modflow 4.2，水文地質(zhì)參數(shù)主要根據(jù)該次工作野外試驗以及以往該區(qū)所開展的地質(zhì)勘查成果初步確定，然后通過擬合流場進行調(diào)參，調(diào)整后的參數(shù)取值見表1。由于模型建立、識別和驗證不是該文的主要研究對象，這里不再對其過程進行詳細(xì)論述。

表1　水文地質(zhì)參數(shù)取值

2.2地下水污染預(yù)測情景設(shè)定

項目污染源多種多樣，但根據(jù)污染源滲漏情況，一般可分為長期滲漏和短期滲漏，該文分別選擇一種情況進行預(yù)測。

(1)長期滲漏。該次對廠區(qū)污水處理站池底出現(xiàn)裂縫，使得污水發(fā)生長期微量滲漏的情況進行預(yù)測。假定池底出現(xiàn)長5m，寬5cm的裂縫，根據(jù)包氣帶土層滲透系數(shù)計算得通過裂縫滲漏的污水量為0.06m3/d。特征污染物為COD，濃度為200mg/L。

(2)短期滲漏。該次對柴油儲罐發(fā)生爆炸破壞地面防滲層，導(dǎo)致石油類污染物短期滲漏的情況進行預(yù)測。假定爆炸破壞地表防滲結(jié)構(gòu)面積為25m2，事故發(fā)生后污染物24小時清理完畢，同樣根據(jù)包氣帶土層滲透系數(shù)計算得通過破壞層的污染物體積為5.86m3/d。特征污染物為輕柴油(石油類)，濃度為6800mg/L。

2.3預(yù)測結(jié)果

該次模擬，當(dāng)預(yù)測點濃度大于相應(yīng)污染物的標(biāo)準(zhǔn)濃度時視為超標(biāo)，標(biāo)準(zhǔn)濃度根據(jù)《地下水質(zhì)量標(biāo)準(zhǔn)》(GB/T14848-93)中的三類標(biāo)準(zhǔn)或《生活飲用水衛(wèi)生標(biāo)準(zhǔn)》(GB5749-2006)確定，COD和石油類分別為3.0mg/L和0.3mg/L。當(dāng)模擬濃度小于于相應(yīng)污染物的檢出限時，視為對地下水無影響，COD和石油類對應(yīng)的檢出限分別為0.1mg/L和0.005mg/L。根據(jù)模擬預(yù)測結(jié)果，在1年末、5年末、10年末其超標(biāo)范圍、影響范圍、最大運移距離見表2和圖2。

3解析法預(yù)測

3.1污染物遷移解析解

同樣對長期滲漏和短期滲漏兩種情況進行預(yù)測，污染場景假設(shè)同前。

表2　地下水污染數(shù)值法預(yù)測結(jié)果

圖2　地下水污染數(shù)值法預(yù)測結(jié)果圖(10年末)

(1)點源定通量連續(xù)注入污染物二維遷移問題解析解。假設(shè)含水層在x，y平面上無限分布且含水層介質(zhì)均質(zhì)各向同性，有點源污染物在含水層厚度上連續(xù)定質(zhì)量注入，地下水沿x方向均勻等速流動，不考慮吸附解吸和化學(xué)反應(yīng)作用[9-10]，其解析解為：

(2)點源瞬時注入污染物二維遷移問題的解析解。假設(shè)含水層在x，y平面上無限分布且含水層介質(zhì)均質(zhì)各向同性，初始時刻在整個含水層厚度上瞬時注入一定量的污染物，地下水沿x方向均勻等速流動，不考慮吸附解吸和化學(xué)反應(yīng)作用，其解析解為：

式中：m為瞬時注入污染物質(zhì)量；其他符號同前。

3.2預(yù)測結(jié)果

各參數(shù)取值與數(shù)值法相同，另外地下水流速按u=k·I/n計算，水力梯度I根據(jù)枯、豐水期等水位線圖確定，最終計算得滲漏點處水流速度約為5.38×10-3m/d。預(yù)測結(jié)果見表3和圖3。

表3　地下水污染解析法預(yù)測結(jié)果

圖3　地下水污染解析法預(yù)測結(jié)果圖(10年末)

4不同預(yù)測方法的預(yù)測結(jié)果分析

根據(jù)預(yù)測結(jié)果，對于解析法，污染物超標(biāo)范圍和影響范圍呈規(guī)則的橢圓形，而數(shù)值法則呈不規(guī)則橢圓形，主要是因為解析法中假定地下水為均勻穩(wěn)定流，而數(shù)值法中地下水流場基本符合地下水的實際運動特征。通過對兩種預(yù)測方法影響范圍和超標(biāo)范圍面積的比較，數(shù)值法均小于解析法，主要是由于數(shù)值法模擬了污染物從地面下滲到含水層中然后再不斷擴散這一整個過程，而解析法中則假定污染物直接在含水層中注入，不考慮下滲過程，也就是說在垂向上，深度越大，數(shù)值法對污染面積的預(yù)測結(jié)果越小，而解析法預(yù)測結(jié)果不隨深度變化。

兩種方法對最大運移距離的預(yù)測結(jié)果較為接近，主要是因為工作區(qū)水力梯度很小且變化不大，使得地下水徑流緩慢且較穩(wěn)定，含水層單一且厚度也相對較穩(wěn)定，另外包氣帶厚度較薄，下滲過程對預(yù)測結(jié)果的影響較小。但如果條件改變，利用解析法求解時就必須考慮這些影響因素，否則預(yù)測結(jié)果將與實際差別較大，比如當(dāng)包氣帶厚度較大或垂向滲透性較差時，污染物下滲至含水層的時間可能較長，若不考慮下滲過程，則最大運移距離預(yù)測結(jié)果可能要大很多，尤其是當(dāng)?shù)叵滤魉佥^大時；另外，在多含水層系統(tǒng)且水力聯(lián)系較密切時，含水層厚度的確定也應(yīng)考慮下滲過程，根據(jù)預(yù)測時段內(nèi)可能影響到的深度確定，否則對預(yù)測結(jié)果的影響很大。

由于解析法引入了更多的各種假定，使得與地下水真實狀態(tài)相差較大，而數(shù)值法則能較好地根據(jù)實際調(diào)查資料建立模型，更為接近地下水實際特征，因此從其理論基礎(chǔ)上來講，數(shù)值法預(yù)測結(jié)果要優(yōu)于解析法。雖然解析法對影響面積和超標(biāo)范圍的預(yù)測結(jié)果和數(shù)值法相比有一定的差異，但未造成質(zhì)的影響，比如對于儲油罐滲漏情況下，數(shù)值模擬法所預(yù)測的橢圓形影響范圍的長軸和短軸分別為148m和57m，解析法則為182m和58m，相差并不是很大，兩種方法對最大運移距離的預(yù)測結(jié)果很接近，因此在該區(qū)采用解析法是適宜的，以其預(yù)測結(jié)果作為污染物防治的依據(jù)也是完全可行的。

5結(jié)論

從理論基礎(chǔ)上講，數(shù)值模擬法預(yù)測結(jié)果優(yōu)于解析法。通過對本例數(shù)值法和解析法預(yù)測結(jié)果的對比分析，其最大距離預(yù)測結(jié)果較為接近，影響面積和超標(biāo)范圍的預(yù)測結(jié)果雖有一定的差異，但未對結(jié)果造成質(zhì)的影響，據(jù)此作為地下水污染防治的依據(jù)也是完全可行的。在對造成預(yù)測結(jié)果差異的原因進行分析的基礎(chǔ)上，認(rèn)為在魯北平原區(qū)，水文地質(zhì)條件相對簡單，同一地段地層分布相對較穩(wěn)定，包氣帶厚度一般較小，可忽略污染物下滲過程，地下水徑流緩慢且較穩(wěn)定，在此條件下采用解析法是適宜的。

參考文獻：

[1]沈媛媛，蔣云鐘，雷曉輝等.地下水?dāng)?shù)值模型在中國的應(yīng)用現(xiàn)狀及發(fā)展趨勢[J].中國水利水電科學(xué)研究院學(xué)報，2009，7(1):57-58.

[2]康鳳新，徐軍祥，張中祥.山東省地下水資源及其潛力評價[J].山東國土資源，2010，26(8):4-6.

[3]劉桂儀.魯北平原深層地下水基本特征與水環(huán)境問題[J].山東國土資源，2001，17(5):43-44.

[4]張中祥，羅斐，張海林.魯西北平原深層地下水動力場時空演化[J].山東國土資源，2011，27(10):32-33.

[5]郭東屏.地下水動力學(xué)[M].陜西：陜西科學(xué)技術(shù)出版社，1994：62-87.

[6]郭東屏，張石峰.滲流理論基礎(chǔ)[M].陜西：陜西科學(xué)技術(shù)出版社，1994：27-32.

[7]薛禹群，謝春紅.地下水?dāng)?shù)值模擬[M].北京：科學(xué)出版社，2007：18-28.

[8]孫訥正.地下水污染-數(shù)值模型和數(shù)值方法[M].北京：地質(zhì)出版社，1989：10-18.

[9]王洪濤.多孔介質(zhì)污染物遷移動力學(xué)[M].北京：高等教育出版社，2008：246-253.

[10]陳崇希，李國敏.地下水溶質(zhì)運移理論及模型[M].北京：中國地質(zhì)大學(xué)出版社，1995：95-102.

收稿日期：2015-09-12；

修訂日期：2015-11-19；編輯：曹麗麗

作者簡介：劉志濤(1982—)，男，河南鞏義人，工程師，主要從事水工環(huán)地質(zhì)工作；E-mail：strawer@126.com

中圖分類號：X523

文獻標(biāo)識碼：B

Comparative Analysis on Prediction Results by Using Numerical Method and Analytical Method of Solute Transport in Groundwater——Taking Zhanhua Power Plant As an Example

LIU Zhitao, ZHOU Qundao, YANG Jianhua

(Lubei Geo-engineering Exploration Institute,Shandong Dezhou 253072, China)

Abstract:At present, numerical method and analytical method are the methods most commonly used for solving groundwater problems. Although numerical method has been widely used for its wide applicability and higher simulation, the analytical method has been one of the first choice for its simple and easy to use. Taking Zhanhua power plant as an example, the pollution caused by long-term and short-term leakage of pollutants have been predicated in this paper. The application suitability of the analytical method in the study area has been discussed based on analyzing the differences between the predictions results gained by using two methods. It will provide reference for the application of this method.

Key words:Groundwater；solute transport；numerical method；analytical method