期刊論文地下水流模擬應(yīng)用

　　【摘 要】介紹了基于有限元原理的FEFLOW軟件,它的應(yīng)用領(lǐng)域覆蓋了水量、水質(zhì)與溫度的模擬,軟件提供地理信息系統(tǒng)數(shù)據(jù)接口,能自動產(chǎn)生空間有限單元網(wǎng)格。FEFLOW程序包中配備了快速精確的數(shù)值算法來控制和優(yōu)化求解過程,結(jié)果輸出具有先進的圖形視覺化效果,并能很好地處理相關(guān)的水文地質(zhì)條件。舉例說明了FEFLOW軟件在建立黃河下游影響帶三維地下水流模擬模型中的應(yīng)用。

　　【關(guān)鍵詞】有限元原 理水文地質(zhì)條件 影響帶 模擬模型

　　利用數(shù)值模擬軟件對地下水流和溶質(zhì)運移等問題進行模擬的方法以其有效性、靈活性和相對廉價性逐漸成為地下水研究領(lǐng)域的一種不可缺少的重要方法,得到了廣泛的應(yīng)用。模型計算指的是進行地下水流動或溶質(zhì)運移正反演計算,常用的方法主要是有限差分法、有限元法、邊界元法等[1]。有限差分法是一種古典的數(shù)值計算方法,其基本思想是:用滲流區(qū)內(nèi)有限個離散本論文由

　　點的集合代替連續(xù)的滲流區(qū),在這些離散點上用差商近似地代替微商,將微分方程及其定解條件化為以未知函數(shù)在離散點上的近似值為未知量的代數(shù)方程(稱之為差分方程),然后求解差分方程,從而得到微分方程的解在離散點上的近似值[2]。美國地質(zhì)調(diào)查局在20世紀(jì)80年代開發(fā)的MODFLOW (Modular three-dimensional finite-difference ground-water flow model)為有限差分的典型代表。該軟件包括水井、補給、河流、溝渠、蒸發(fā)蒸騰和通用水頭邊界六個子程序包,分別用來處理相關(guān)的水文地質(zhì)條件[3]。自MODFLOW問世以來,不斷有新的子程序包被開發(fā)出來。例如用來模擬水位下降引起地面沉降的子程序包[4],用來模擬水平流動障礙(Horizontal flow-barrier)的子程序包等[5]。新子程序包的加入,使MODFLOW的應(yīng)用范圍不斷擴大。實踐證明,只要合理使用,它也可以用來解決許多地下水在裂隙介質(zhì)中的流動問題。經(jīng)過合理的線性化處理,MODFLOW還可以用于解決空氣在土壤中的運動問題[6]。如果將它與其他用于溶質(zhì)運移模擬的程序結(jié)合起來,可以用來模擬諸如海水入侵等地下水密度為變量的問題[7]。有限單元法是求解偏微分方程定解問題的另一種有效的數(shù)值方法,它的基本思想早在20世紀(jì)40年代就已經(jīng)提出。有限單元法成功地解決了眾多領(lǐng)域的許多計算問題,例如結(jié)構(gòu)力學(xué)、結(jié)構(gòu)工程學(xué)中的應(yīng)力分析,結(jié)構(gòu)穩(wěn)定性問題;土力學(xué)、巖石力學(xué)中的應(yīng)力變形與穩(wěn)定性分析等;水利工程的液體流動、水工結(jié)構(gòu)和壩體分析以及地下水流動、水動力彌散問題等。20世紀(jì)70年代末,德國WASY水資源規(guī)劃和系統(tǒng)研究所開發(fā)了基于有限單元法的FEFLOW ( Finite element subsurface FLOWsystem)軟件,它是迄今為止功能最為齊全的地下水模擬軟件包之一。該軟件包具有圖形人機對話、地理信息系統(tǒng)數(shù)據(jù)接口、自動產(chǎn)生空間各種有限單元網(wǎng)格、空間參數(shù)區(qū)域化及快速精確的數(shù)值算法和先進的圖形視覺化技術(shù)等特點。在FEFLOW系統(tǒng)中,用戶可以很方便迅速地產(chǎn)生空間有限單元網(wǎng)格,設(shè)置模型的參數(shù)和定義邊界條件,運行數(shù)值模擬以及實時圖形顯示結(jié)果與成圖。

　　1 FEFLOW簡介從FEFLOW的問世起一直到現(xiàn)在,在理論研究和對實際問題的處理上,它經(jīng)過了不斷的發(fā)展、修改、提高,在它的發(fā)展過程中,FEFLOW經(jīng)過了大量的測試和檢驗,它成功地解決了一系列與地下水有關(guān)的實質(zhì)性問題,如判斷污染物遷移途徑、追溯污染物的來源,地?zé)岬哪M,海水入侵預(yù)測等。

　　1.1 FEFLOW的應(yīng)用領(lǐng)域水量模擬:模擬水源地開采或者油田注水對區(qū)域地下水流場的影響、模擬水庫放水或者河流斷流時,河道沿線地下水流場的變化等;水質(zhì)模擬:模擬污染物在地下水中的遷移過程及其時間空間分布模式[8]、模擬沿海地區(qū)抽取地下水引起的海水入侵等[9];溫度模擬:模擬非飽和帶以及飽和帶溫度場的分布[10]。

　　1.2 系統(tǒng)輸入特點通過標(biāo)準(zhǔn)數(shù)據(jù)輸入接口,用戶既能直接利用已有的GIS空間多邊形數(shù)據(jù)生成有限單元網(wǎng)格,也可以基于地圖用鼠標(biāo)設(shè)計,網(wǎng)格的數(shù)目可以自己指定,可以方便地調(diào)整網(wǎng)格的幾何形狀,增加和放疏網(wǎng)格大小等。用戶可以視具體情況定義第一、第二和第三類邊界,而且還可以對邊界條件增加特定的限制條件;能方便地定義復(fù)合含水層中的分層開采井和混合開采井以及注水井。所有邊界條件及其限制條件、滲透系數(shù)、補排量既可設(shè)置為常數(shù),也能定義為隨時間變化的函數(shù)。FEFLOW提供了克里格(Kriging)、阿基瑪(Akima)和距離反比加權(quán)(IDW)三種方法,對離散的空間抽樣數(shù)據(jù)進行內(nèi)插或外推。輸入數(shù)據(jù)格式既可以是ASCⅡ碼文件,也可以是GIS地理信息系統(tǒng)文件,如ArcView的shp格式。

　　1.3 FEFLOW配備了若干先進的數(shù)值求解法來控制和優(yōu)化求解過程 快速直接求解法,如PCG,BICGSTAB,CGS,GMRES以及帶預(yù)處理的再啟動ORTHOMIN法;靈活多變的up-wind技術(shù),如流線up-wind,奇值捕捉法Shock capturing,以減少數(shù)值彌散;皮卡和牛頓迭代法求解非線性流場問題;自動調(diào)節(jié)模擬時間步長;模擬污染物遷移過程包括對流、水動力彌散、線性及非線性吸附、一階化學(xué)非平衡反應(yīng);為非飽和帶模擬提供了多種參數(shù)模型,如指數(shù)式、Van Genuchten式和多種形式的Richard方程;垂向滑動網(wǎng)格(BASD)技術(shù)處理自由表面含水系統(tǒng)以及非飽和帶的模擬[11];開放性外部程序接口,以便用戶在FEFLOW系統(tǒng)中連接和使用自己的程序模塊。1.4 系統(tǒng)結(jié)果輸出及顯示FEFLOW提供了其他任何地下水模擬軟件都無法比擬的、豐富實用的圖形顯示和數(shù)據(jù)結(jié)果分析工具。其先進的圖形視覺化及數(shù)據(jù)分析技術(shù)表現(xiàn)在:有限單元網(wǎng)、邊界條件和模型參數(shù)的三維可視化及顯示;標(biāo)量數(shù)據(jù)的三維彩色(透明或灰度)等勢面顯示以及其二維平面彩色或等值線顯示;三維地下水流徑追蹤,流動時間及流速動畫顯示(包括二維平面、剖面投影或二維平面追蹤);三維體截段的空間顯示和三維交叉剖面組的空間顯示;三維圖形的任意旋轉(zhuǎn),二維、三維圖形的放大或縮小;總體和局部水量平衡分析(包括任意幾何多邊形內(nèi)的水流通量分析);計算和圖形顯示通過各種邊界條件的水通量、物質(zhì)通量以及其在特定時間區(qū)間內(nèi)的積分量。

　　2 FEFLOW中特殊水文地質(zhì)問題的處理

　　2.1 地層的不連續(xù)在FEFLOW中,三維模型的空間結(jié)構(gòu)是由許多層(layer)組成,這些層由片(slice)隔開,每兩片中就是一層(圖1)。網(wǎng)格、邊界條件、初始條件是對應(yīng)于片,而參數(shù)(給水度、滲透系數(shù)等)是對應(yīng)于層。FEFLOW只能處理延伸到整個研究區(qū)域的片和層,所以,理論上FEFLOW不能處理地層不連續(xù)等特殊情況。

　　2.1.1 地層的缺失地層由于地殼的抬升而受強烈的剝蝕,含水巖層出露地面(圖1),或者含水層缺失。FEFLOW在處理這種情況時,仍然認(rèn)為地層是連續(xù)到整個研究區(qū)域;只是在地層缺失的地方,認(rèn)為厚度是無限小,默認(rèn)是0.01 m,其參數(shù)參考相鄰層的參數(shù)值。

　　2.1.2 斷層斷層是地層發(fā)生錯位,是一種特殊的水文地質(zhì)條件。到目前為止,FEFLOW還沒有很好的方法來處理斷層、跌水現(xiàn)象,WASY公司也一直設(shè)法解決這一難題。常用的處理方法是把斷層部分作為一個單獨的參數(shù)分區(qū),根據(jù)斷層兩邊的水力聯(lián)系程度來確定含水層的參數(shù)(主要是滲透系數(shù))。

　　2.1.3 小面積強透水帶的處理當(dāng)遇到天窗或者滲透系數(shù)非常大的區(qū)域時(這些區(qū)域的面積往往比較小),可以把它隔離出來當(dāng)成一單獨的層(圖2)。在這一特殊的層上,強透水帶位置的厚度、滲透系數(shù)等參數(shù)按照實際賦值;除此之外的地方,參數(shù)考慮相鄰層相應(yīng)的參數(shù)。

　　2.2 第三類邊界混合邊界河流在FEFLOW中是作為第三類邊界處理。河流通過底積層和含水層發(fā)生水力聯(lián)系(圖3)。根據(jù)達西定律有當(dāng)動態(tài)地模擬一條河流時,需要知道河流在典型水文年的水位動態(tài)曲線。例如已知兩個水文站的水位歷時曲線,采用一維線性插值方法,可以得到沿主流道任意結(jié)點的水位歷時曲線。然而這種插值法用來模擬短期內(nèi)洪水對地下水的影響時就會出現(xiàn)較大的誤差,因為它認(rèn)為上游和下游水位變化是同步的,沒有考慮時間滯后的影響;也就是說上游和下游的水位同時到達歷時曲線的波峰和波谷(圖4),這種結(jié)果并不是所需要的。為了解決上述問題,FEFLOW提供了一種考慮時間滯后的一維線性插值方法,通過這種插值(圖5),可以非常準(zhǔn)確地模擬一場暴雨后,沿河附近區(qū)域地下水位的變化。

　　小編推薦優(yōu)秀電子期刊　《中國新通信》是國家級期刊嗎

　　《中國新通信》雜志是由工業(yè)和信息化部主管，電子工業(yè)出版社主辦，《中國新通信》雜志社編輯出版的信息通信技術(shù)專業(yè)期刊，原《中國數(shù)據(jù)通信》(月刊)。自1999年創(chuàng)刊以來，受到了廣大讀者和作者的關(guān)心和支持，2005年被科學(xué)技術(shù)協(xié)會評為優(yōu)秀期刊，于2006年1月改名為《中國新通信》。