李勁彬

摘 要：地下水作為重要的環(huán)境資源，其安全受到人們的廣泛關(guān)注和高度重視。地下水?dāng)?shù)值模擬作為評價地下水水量和水質(zhì)的重要技術(shù)手段，因而也受到了相關(guān)科研人員的廣泛關(guān)注。文章在分析地下水?dāng)?shù)值模擬發(fā)展歷程的基礎(chǔ)上，簡單介紹了國內(nèi)外幾種使用較為頻繁的地下水?dāng)?shù)值模擬軟件，如GMS，Visual MODFLOW，F(xiàn)EFLOW等，闡述了地下水?dāng)?shù)值模擬的發(fā)展趨勢。

關(guān)鍵詞：地下水；數(shù)值模擬；GMS；Visual MODFLOW；FEFLOW

中圖分類號：TV1 文獻標(biāo)志碼：A 文章編號：2095-2945（2018）07-0044-02

Abstract： Groundwater is an important environmental resource， and its safety has been widely concerned about and highly valued. Numerical simulation of groundwater is an important technical means to evaluate the quantity and quality of groundwater. On the basis of analyzing the development course of groundwater numerical simulation， this paper briefly introduces several kinds of software of groundwater numerical simulation， such as GMS， Visual MODFLOW， FEFLOW， etc.， which are used frequently at home and abroad， and the development trend of numerical simulation of groundwater is also described.

Keywords： groundwater； numerical simulation； GMS； Visual MODFLOW； FEFLOW

引言

地下水作為一種十分寶貴的自然資源，同時也作為城市生活用水、工業(yè)用水、農(nóng)業(yè)用水等重要的供水水源，與人類生活息息相關(guān)。近些年，由于城鎮(zhèn)化發(fā)展進程的加快，地下水長期處于超采狀態(tài)[1]，形成區(qū)域降落漏斗[2]。而且由于建設(shè)項目的增多，地下水水質(zhì)也時刻受到威脅[3]。因此如何科學(xué)合理利用好地下水水資源，關(guān)系到我國國民經(jīng)濟是否能夠健康平穩(wěn)可持續(xù)發(fā)展。通過數(shù)值模擬方法[4]，可以可視化分析研究地下水流和溶質(zhì)運移問題，其方法直觀、有效、靈活，因而越來越受到領(lǐng)域內(nèi)學(xué)者的重視以及廣泛應(yīng)用。目前其主要應(yīng)用于地下水資源評價、預(yù)報、管理以及地下水中溶質(zhì)運移擴散研究當(dāng)中[4]，并且取得了一定的研究成果。

1 數(shù)值模擬概述

數(shù)值模擬法是以描述地下水流運動微分方程、邊界條件和初始條件為基礎(chǔ)，建立地下水水位隨時間連續(xù)變化的函數(shù)數(shù)學(xué)模型，并在計算機上采用離散化方法對所建立的數(shù)學(xué)模型進行求解，所求得的解是一種近似解。數(shù)值模擬法在研究地下水流動態(tài)方面的優(yōu)點就在于其不僅可以有效地反映復(fù)雜地質(zhì)條件下地下水運動情況，而且還可以解決地下水水質(zhì)以及熱量運移等問題。常用的數(shù)值模擬方法包括有限差分法[5]、有限單元法[6]、有限分析法[7]、邊界單元法[8]等，其中最常用的方法為有限差分法和有限單元法[9]。有限差分法于20世紀(jì)60年代開始逐漸應(yīng)用于地下水流模擬計算，其簡單、靈活、通用性強，計算機上容易實現(xiàn)。后者由我國數(shù)學(xué)家馮康創(chuàng)建，并于20世紀(jì)60年代后期開始應(yīng)用于地下水?dāng)?shù)值模擬計算之中[10]。

2 地下水?dāng)?shù)值模擬發(fā)展簡史

人類對地下水的開發(fā)利用有著悠久的歷史，從大約5600年前中國先祖?zhèn)兺ㄟ^木質(zhì)結(jié)構(gòu)井取水開始，到1856年法國著名水利工程師達西（H. Darcy）提出揭示地下水滲流運動規(guī)律的達西定律，之后很長時間水文地質(zhì)關(guān)于地下水的計算主要以解析法為主，到了20世紀(jì)50～60年代，離散介質(zhì)點網(wǎng)絡(luò)模擬成為地下水主要的計算方法[11]。1956年斯圖爾曼（R. W. Stallman）首次將數(shù)值法應(yīng)用于水文地質(zhì)計算，隨著華爾頓（W. C. Walton）將電子計算機引入地下水?dāng)?shù)值模擬，便利了水文地質(zhì)數(shù)學(xué)問題的求解，促進了水文地質(zhì)問題的定量求解，加速了數(shù)值模擬的發(fā)展。至此，越來越多與地下水流相關(guān)的各類問題都通過數(shù)值模擬方法求解。

3 地下水模擬軟件概述

地下水?dāng)?shù)值模擬軟件隨著計算機的引入而快速發(fā)展，為了方便、快捷、直觀解決地下水流相關(guān)問題，國外科研團隊研究開發(fā)了眾多地下水?dāng)?shù)值模擬軟件。其中主要有GMS[12]、Visual MODFLOW[12]以及FEFLOW[13]等。

3.1 GMS

Groundwater Modeling System，簡稱GMS，即地下水模擬系統(tǒng)，由美國Brigham Young環(huán)境模擬研究中心與美國軍隊排水工程試驗工作站在已有地下水模擬軟件的基礎(chǔ)上共同合作研究開發(fā)，該軟件以概念模型和網(wǎng)格模型為基礎(chǔ)， 軟件界面綜合了MODFLOW、SEAM3D、MODPATH、MT3DMS等可視化三維地下水模擬軟件包，前處理功能優(yōu)越，建模便捷直觀，圖形圖像處理功能強大，是國內(nèi)外目前功能最為完善的地下水模擬軟件。

3.2 Visual MODFLOW

Visual MODFLOW是一款在美國地質(zhì)勘探局（USGS）軟件MODFLOW的基礎(chǔ)之上研究開發(fā)的地下水流模擬軟件，由加拿大Waterloo水文地質(zhì)公司研究開發(fā)，其綜合了已有的 MT3D、MODPATH、Win PEST以及RT3D等地下水?dāng)?shù)值模型，并將現(xiàn)代可視化技術(shù)應(yīng)用其中，使之成為目前國際上最具影響力的地下水流和溶質(zhì)運移可視化模擬軟件。該軟件在我國目前主要應(yīng)用于地下水資源評價、地下水水位動態(tài)預(yù)測、建設(shè)項目選址、地下水污染物運移、煤礦涌水量預(yù)測等方面。endprint

3.3 FEFLOW

FEFLOW軟件在1979年由德國WASY水資源規(guī)劃系統(tǒng)研究所研究開發(fā)，其作為有限元地下水?dāng)?shù)值模擬的代表，有利于用戶進行模型建立、分析復(fù)雜三維地質(zhì)體地下水流動態(tài)以及溶質(zhì)運移情況，該軟件是迄今為止功能最齊全的地下水水量和水質(zhì)計算模擬軟件[14]。其主要優(yōu)勢表現(xiàn)在模擬非飽和帶與飽和帶地下水流及其溫度分布、應(yīng)用ArcGIS和FEFLOW進行地下水空間分布模擬等方面。

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

我國在地下水?dāng)?shù)值模擬方面的研究起步相對較晚，大概開始于20世紀(jì)70年代[15]。經(jīng)過國內(nèi)專家學(xué)者20多年的不懈努力，我國地下水?dāng)?shù)值模擬研究飛速發(fā)展，應(yīng)用遍及地下水相關(guān)的各個領(lǐng)域，并且針對應(yīng)用過程中存在的問題，科研工作者不斷進行理論和方法上的創(chuàng)新，拓寬研究方向和加深研究深度，現(xiàn)如今我國的地下水?dāng)?shù)值模擬研究已經(jīng)接近國際化水平。國際上出現(xiàn)的地下水?dāng)?shù)值模型我國已基本具有，如地下水污染模型[16]、壩體滲漏模型[17]、地下水資源管理模型[18]、熱運移模型[19]、溶質(zhì)運移模型[20]、地面沉降模型[21]等。近些年，隨著地理信息系統(tǒng)技術(shù)不斷快速發(fā)展，其與地下水?dāng)?shù)值模型的整合也越來越受到國內(nèi)科研人員的重視，整合后的模型將具有整體化、自動化、實時化等優(yōu)點。

5 發(fā)展趨勢

隨著世界經(jīng)濟的發(fā)展和科學(xué)技術(shù)的不斷進步，地下水資源需求量不斷增加，與地下水相關(guān)各類問題的范圍與復(fù)雜性也日益增大，人類活動對地下水水位以及水質(zhì)的影響也日益加劇，因此合理制定開采方案并進行環(huán)境影響預(yù)報顯得尤為重要，地下水?dāng)?shù)值模擬發(fā)揮著越來越重要的作用。但由于實際水文地質(zhì)條件復(fù)雜，野外試驗數(shù)據(jù)缺乏，學(xué)科交叉應(yīng)用不強等原因，從而導(dǎo)致數(shù)值模擬的精度不高。近年來，隨著計算機技術(shù)的飛速發(fā)展以及3S技術(shù)（遙感技術(shù)、地理信息系統(tǒng)技術(shù)、全球定位系統(tǒng)）的不斷成熟，地下水?dāng)?shù)值模擬軟件與其結(jié)合必將對模擬精度的提高帶來質(zhì)的飛躍。其次，隨著人類對地下水賦存環(huán)境及其演變規(guī)律認識的不斷深入，多學(xué)科多方法的耦合研究，將使得地下水?dāng)?shù)值模擬研究方法具有更廣闊的前景[22]。

參考文獻：

[1]陳飛，侯杰，于麗麗，等.全國地下水超采治理分析[J].水利規(guī)劃與設(shè)計，2016（11）：3-7.

[2]陳宇良，都基眾.地下水降落漏斗的形成與控制[J].工程勘察，2007（06）：23-26+73.

[3]龔星，陳植華，孫璐.地下水環(huán)境影響評價若干關(guān)鍵問題探討[J].安全與環(huán)境工程，2013，20（02）：95-99.

[4]孫從軍，韓振波，趙振，等.地下水?dāng)?shù)值模擬的研究與應(yīng)用進展[J].環(huán)境工程，2013，31（05）：9-13+17.

[5]張宏仁.解地下水滲流問題的有限差分法（一）[J].水文地質(zhì)工程地質(zhì)，1980（01）：53-57.

[6]阮淼森.有限單元法在地下水資源管理中計算水位動態(tài)的應(yīng)用[J].工程勘察，1982（03）：40-43.

[7]周志芳，李艷.復(fù)雜巖體地下水運動問題的有限分析法[J].水科學(xué)進展，1997（03）：36-42.

[8]邊際.邊界單元法在地下水非穩(wěn)定流計算中的初步應(yīng)用[J].水利學(xué)報，1984（04）：20-30.

[9]童彥釗，吳平，韓強強.地下水?dāng)?shù)值模擬軟件簡述[J].科技信息，2011（29）：480+459.

[10]王文科.求解地下水非穩(wěn)定流問題的局部坐標(biāo)有限分析法[J].長春地質(zhì)學(xué)院學(xué)報，1995，25（4）：415-422.

[11]張人權(quán)，粱杏，靳孟貴，等.水文地質(zhì)學(xué)基礎(chǔ)（第六版）[M].北京：地質(zhì)出版社，2011：1-13.

[12]王慶永，賈忠華，劉曉峰，等.VisualMODFLOW及其在地下水模擬中的應(yīng)用[J].水資源與水工程學(xué)報，2007（05）：90-92.

[13]賀國平，邵景力，崔亞莉，等.FEFLOW在地下水流模擬方面的應(yīng)用[J].成都理工大學(xué)學(xué)報（自然科學(xué)版），2003（04）：356-361.

[14]郭曉東，田輝，張梅桂，等.我國地下水?dāng)?shù)值模擬軟件應(yīng)用進展[J].地下水，2010，32（4）：5-7.

[15]薛禹群.中國地下水?dāng)?shù)值模擬的現(xiàn)狀與展望[J].高校地質(zhì)學(xué)報，2010，16（01）：1-6.

[16]薛禹群，吳吉春，謝春紅.越流含水層系統(tǒng)地下水污染數(shù)值模擬[J].地質(zhì)學(xué)報，1997，71（2）：186-192.

[17]李恩羊.滲灌條件下土壤水分運動的數(shù)學(xué)模擬[J].水利學(xué)報，1982（4）：1-10.

[18]邵景力，魏加華，崔亞莉，等.用遺傳算法求解地下水資源管理模型[J].地球科學(xué)——中國地質(zhì)大學(xué)學(xué)報，1998，23（5）：532-536.

[19]黎明，李國敏，楊遼，等.陜西渭北東部巖溶熱水運移數(shù)值模型研究[J].中國科學(xué).D輯：地球科學(xué)，2006（S2）：33-38.

[20]孫克讓，徐紹輝，朱學(xué)愚.地下水污染模型的特征有限元解法

──以佳木斯為例[J].太原工業(yè)大學(xué)學(xué)報，1994（04）：21-27.

[21]郭占榮，曲煥林，崔小東，等.天津市地面沉降數(shù)值模擬研究[J].地球?qū)W報，1998（04）：80-87.

[22]郝治福，康紹忠.地下水系統(tǒng)數(shù)值模擬的研究現(xiàn)狀和發(fā)展趨勢[J].水利水電科技進展，2006（01）：77-81.endprint