地下水數值模擬研究

李勁彬

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

關鍵詞：地下水；數值模擬；GMS；Visual MODFLOW；FEFLOW

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

1 數值模擬概述

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

2 地下水數值模擬發(fā)展簡史

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

3 地下水模擬軟件概述

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

3.1 GMS

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

3.2 Visual MODFLOW

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

3.3 FEFLOW

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

4 國內研究進展

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

5 發(fā)展趨勢

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

參考文獻：

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

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

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

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

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

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

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

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

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

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

[11]張人權，粱杏，靳孟貴，等.水文地質學基礎（第六版）[M].北京：地質出版社，2011：1-13.

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

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

[14]郭曉東，田輝，張梅桂，等.我國地下水數值模擬軟件應用進展[J].地下水，2010，32（4）：5-7.

[15]薛禹群.中國地下水數值模擬的現狀與展望[J].高校地質學報，2010，16（01）：1-6.

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

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

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

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

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

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

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

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