高靖勛

摘 要：地下水是水資源的重要組成部分，對于社會的經(jīng)濟發(fā)展和人們的日常生活來說都是不可或缺的。該文結合地下水數(shù)值模擬的特點，分析了數(shù)值模擬方法在實際工程應用和地下水污染控制與修復等方面的優(yōu)勢，介紹了國內(nèi)外地下水數(shù)值模擬及相關軟件的研究進展。

關鍵詞：地下水 數(shù)值模擬 模擬軟件 發(fā)展趨勢

中圖分類號：P641.7 文獻標識碼：A 文章編號：1672-3791（2019）01（c）-0088-03

水是生命之源，地下水作為水資源的主要組成部分之一，在社會的經(jīng)濟發(fā)展和人們的日常生活中都發(fā)揮著重要的作用。但隨著經(jīng)濟的高速發(fā)展和社會的不斷進步，人們對地下水的需求量節(jié)節(jié)攀升。受工農(nóng)業(yè)、生活等人類活動的影響，地下水污染問題也日趨嚴重。如何經(jīng)濟、高效、合理地處理實際工程中出現(xiàn)的水文地質(zhì)問題、進行地下水污染修復成為了當下備受關注的問題。

地下水數(shù)值模擬是一種以計算機、軟件為載體，利用地下水數(shù)值模型模擬計算不同條件下地下水流、溶質(zhì)運移情況的技術方法，具有經(jīng)濟有效、靈活性強等特點[1]，在地下水研究領域應用廣泛。進行數(shù)值模擬對地下水資源保護、地下水污染評估、污染修復方案優(yōu)化具有非常重要的指導作用。

1 國外研究進展

法國水力學家達西（H.Darcy）進行了大量均勻沙柱滲流實驗，于1856年提出了達西定律[2]，奠定了地下水定量計算的基礎。之后，裘布依以Darcy定律為基礎，與1863年提出了地下水一維穩(wěn)定流和向井二維穩(wěn)定流[3]。1935年，在數(shù)學家C.I.盧賓的輔助下，泰斯（C.V.Theis）推導出了描述地下水非穩(wěn)定流動的泰斯公式，奠定了非穩(wěn)定流理論的基礎。

伴隨著計算機技術的發(fā)展，地下水定量計算也從解析法、物理模擬法進入了數(shù)值模擬階段。1935年初，斯圖爾曼等（R.W.Stallman）[4]率先在地下水研究領域應用了數(shù)值模擬方法。數(shù)值法利用數(shù)學方法可以將不規(guī)則的連續(xù)地下水系統(tǒng)離散化，得到離散解，更加準確地刻畫含水層的非均質(zhì)性和各項異性、入滲、蒸發(fā)等[5]。數(shù)值模擬方法主要包括有限差分法、有限元法、邊界元法、積分有限差分法等，目前比較常用的是有限元法和有限差分法。

1964年，Tyson和Weber在模擬美國加利福尼亞州地下水流時運用了有限差分法。之后，Bedehoeft[6]于1969應用有限差分法進行了地下水數(shù)值模擬，使該方法的影響力進一步擴大。

1965年，有限元法被引入地下水數(shù)值模擬領域[7]，這種方法對邊界條件比較復雜的水文地質(zhì)區(qū)適應性更好，占用內(nèi)存小，精度較高。三維等參有限元法由Capta等提出，Huyakom在1986年應用有限元法進行了孔隙、裂隙介質(zhì)中的地下水數(shù)值模擬[8]。Numan認為，使用有限元法處理非穩(wěn)定流問題時，不同時間步長的選擇可能產(chǎn)生不穩(wěn)定的解，對此，Wood等[9]進行了相關研究，確定了有限元計算二維流的時間步長條件。

1998年，Scheibe等[10]發(fā)表了對三維流及溶質(zhì)運移模擬模型尺度問題的研究結果。2000年，Porter等[11]發(fā)表了對于有限差分法在地質(zhì)學、地球物理學、水文學綜合性應用的看法。為了有效解決地下水位預測中預測結果不確定這一問題，Li等[12]在2003年提出了隨機地下水模型。

2 國內(nèi)研究進展

同國外相比，中國在地下水數(shù)值模擬領域的研究起步較晚。1974年，劉明新、陳忠祥[13]率先引入數(shù)值模擬方法處理地層滲流問題。進入20世紀90年代后，地下水數(shù)值模擬開始在我國快速發(fā)展起來[14]。薛禹群[15]認為，多尺度有限元法比傳統(tǒng)的有限元法更加適合處理一般的非均質(zhì)多孔介質(zhì)條件。

隨著計算機科學的高速發(fā)展，計算機運算速度和運算能力不斷提高，地下水數(shù)值模擬技術在我國的應用也越來越廣泛。盧文喜[16]認為，在進行邊界條件預報時，需要考慮人為和自然因素的變化以及相鄰地區(qū)水流的收支。王浩然等[17]建立了有限元-邊界元耦合模型，避免了人為邊界引發(fā)的流場失真。杜新強[18]等分析了怎樣推求三維流模擬中弱透水層的初始水位。

數(shù)學模型是對地下水流、含水層系統(tǒng)、污染物運移以及它們之間可能存在的物理、化學、生物反應的概化，因此會對相關條件進行許多假設。為了在當前計算能力可解模型的前提下使模擬結果更加科學準確，避免地下水污染修復系統(tǒng)設計失效，處理地下水系統(tǒng)參數(shù)的固有不確定性，就變成了研究過程中的關鍵環(huán)節(jié)。

水文地質(zhì)參數(shù)的確定方法，往往會對模擬結果產(chǎn)生很大的影響，我國許多水文地質(zhì)學家都對此進行了參數(shù)最優(yōu)解方面的研究。目前，地下水數(shù)值模擬也在各個相關領域被廣泛應用[14]。孫澤清[19]選取云南某煉油廠進行了石油泄漏污染物運移的模擬和預測。張浩佳等[20]應用GSFLOW模型對開都河流域進行了地表水-地下水耦合模擬，為地表水-地下水相交互用評估和資源統(tǒng)籌提供了科學依據(jù)。

3 地下水數(shù)值模擬軟件

在計算機技術的推動下，出現(xiàn)了各種地下水數(shù)值模擬軟件。從初始的DOS版本軟件如MODFLOW、PEST、MT3DMS[21]到如今的可視化軟件如GMS、Visual MODFLOW、FEFLOW[8]等，地下水數(shù)值模擬的準確性和效率不斷提高。GMS、Visual MODFLOW也是目前在科研院所和生產(chǎn)單位中應用較廣的模擬軟件，它們對穩(wěn)定、非穩(wěn)定流和溶質(zhì)運移的模擬結果都比較可靠。

4 存在問題

地下水數(shù)值模型可以對水文地質(zhì)問題進行定量刻畫，因此利用地下水數(shù)值模擬技術得到解的準確度通常也比較高。但是，準確的數(shù)值模擬結果必定是建立在充足的水文地質(zhì)資料和水文地質(zhì)條件的合理概化之上的。目前，實際項目中，水文地質(zhì)資料的缺失十分常見，這無疑給模擬邊界的圈定和機理分析帶來了困難，同時也降低了模擬結果的精度。而且，如果過分依賴計算機軟件模擬而忽略了實際水文地質(zhì)條件的研究，那模擬結果的準確性就完全無從談起了。

5 前景展望

地下水作為水資源的重要組成部分，其污染控制與修復工作對環(huán)境保護和可持續(xù)發(fā)展有著重要的意義[22]。隨著信息技術的發(fā)展和信息獲取手段的不斷擴充，地下水模型的精度必定會逐步提高，地下水數(shù)值模擬技術與數(shù)學、化學、物探等學科的交叉聯(lián)動必定也會越來越多，這也會為各種耦合模型的開發(fā)和地下水數(shù)值模擬技術的發(fā)展注入源源不斷的活力。

參考文獻

[1] 李凡，李家科，馬越，等.地下水數(shù)值模擬研究與應用進展[J].水資源與水工程學報，2018（1）：99-104.

[2] 張人權.水文地質(zhì)學基礎[M].北京：地質(zhì)出版社，2011.

[3] 薛禹群，謝春紅.地下水數(shù)值模擬[M].北京：地質(zhì)出版社，2007：1-5.

[4] Stallman RW. Numerical analysis of regional water levels to define aquifer hydrology[J].EOS Transactions American Geophysical Union，1956，37（4）：451-460.

[5] 何建東.延安市新區(qū)造地工程地下水流數(shù)值模擬研究[D].合肥工業(yè)大學，2015.

[6] Bredehoeft JD. Finite difference approximations to the equations of Ground-Water flow[J].Water Resources Research，1969，5（2）：531-534.

[7] Pinder GF，Bredehoeft JD.Application of digital computer for aquifer evaluation[J].Water Resources Research，1968，4（5）：1069-1093.

[8] 侯嘉維.馬海盆地地下水數(shù)值模擬與資源評價[D].吉林大學，2016.

[9] Wood WL. A note on how to avoid spurious oscillation in the finite-element solution of the unsaturated flow equation[J].Journal of Hydrology，1996，176（1）：205-218.

[10] Scheibe T， Yabusaki S. Scaling of flow and transport behavior in heterogeneous groundwater systems[J].Advances in Water Resources，1998，22（3）：223-238.

[11] Porter DW， Bruce PG， Walter FJ， et al. Data fusion modeling for groundwater systems[J]. Journal of Contaminant Hydrology，2000，42（2）：303-335.

[12] LI Shu-guang， Mclaughlin D， LIAO Hua-sheng. A computationally practical method for stochastic groundwater modeling[J]. Advances in Water Resources， 2003，26（11）：1137-1148.

[13] 劉明新，陳鐘祥.單重和雙重孔隙介質(zhì)二相二維滲流的一個新的數(shù)值模擬方法[J].力學學報，1984（3）：225-233.

[14] 韓小強.內(nèi)蒙古黃旗海盆地三維地質(zhì)建模及地下水數(shù)值模擬[D].南京大學，2018.

[15] 魏林宏，束龍倉，郝振純.地下水流數(shù)值模擬的研究現(xiàn)狀和發(fā)展趨勢[J].重慶大學學報：自然科學版，2000（S1）：50-52.

[16] 盧文喜.地下水運動數(shù)值模擬過程中邊界條件問題探討[J].水利學報，2003（3）：33-36.

[17] 王浩然，朱國榮，趙金熙.基于區(qū)域分解法的地下水有限元與邊界元耦合模型——淄博市王旺莊水源地地下水數(shù)值模擬[J].地質(zhì)論評，2003（1）：48-52.

[18] 杜新強，常翠，冶雪艷.地下水流三維數(shù)值模擬中弱透水層初始水位的推求方法分析[J].黑龍江水專學報，2008，35（4）：94-97.

[19] 孫澤清.云南某煉油廠地下水污染物運移數(shù)值模擬研究[D].華東理工大學，2017.

[20] 張浩佳，吳劍鋒，林錦，等.GSFLOW在干旱區(qū)地表水與地下水耦合模擬中的應用[J].南京大學學報：自然科學版，2015，51（3）：596-603.

[21] 丁繼紅，周德亮，馬生忠.國外地下水模擬軟件的發(fā)展現(xiàn)狀與趨勢[J].勘察科學技術，2002（1）：37-42.

[22] 劉志陽.地下水污染修復技術綜述[J].環(huán)境與發(fā)展，2016，28（2）：1-4.