呂 敏

（遼寧省丹東水文局，遼寧丹東118000）

?

基于二維非穩(wěn)態(tài)模型的地下水資源預(yù)測研究

呂 敏

（遼寧省丹東水文局，遼寧丹東118000）

摘要：本文根據(jù)國內(nèi)外地下水資源預(yù)測的研究進(jìn)展，結(jié)合灌區(qū)的地質(zhì)條件，運用二維非穩(wěn)態(tài)模型對灌區(qū)地下水資源分布情況進(jìn)行了預(yù)測，經(jīng)與實測值對比，表明二維非穩(wěn)態(tài)模型的計算精度較高，可以準(zhǔn)確地反映出地下水資源分布情況，研究結(jié)果可為地下水資源預(yù)測提供參考。

關(guān)鍵詞：地下水；預(yù)測；非穩(wěn)態(tài)模型；灌區(qū)

地下水掠奪式開采導(dǎo)致了一系列環(huán)境問題，因此地下水保護(hù)迫在眉睫。地下水系統(tǒng)一般包括含水系統(tǒng)和地下徑流。由于地下水是不可見的，且受地質(zhì)條件影響遠(yuǎn)大于地表水，因此地下水資源預(yù)測比地表水預(yù)測難得多。

從國內(nèi)外的研究情況看，地下水測量主要經(jīng)歷了原始測量、傳統(tǒng)測量、數(shù)字測量三個階段［3-4］。目前，一般使用GPS定位和探測儀相結(jié)合的方法對地下水進(jìn)行測量。隨著數(shù)值模擬和水動力學(xué)的發(fā)展，計算流體力學(xué)逐漸成為學(xué)者們的研究熱點。

本文采用的二維非穩(wěn)態(tài)地下水資源預(yù)測模型由連續(xù)性方程、運動方程、控制方程組成［5］。

1　地下水資源分布

研究區(qū)域地下水主要為潛水和承壓水，潛水埋藏深度一般是2～10m，承壓水埋藏深度大多在90m以下。由于承壓水埋藏深度遠(yuǎn)大于潛水，其開采成本往往是潛水的數(shù)倍乃至數(shù)十倍，因此目前人類大量開采的灌溉水多為潛水，本文的地下水資源預(yù)測的主要對象也為潛水。

在研究區(qū)域設(shè)置了38個地下水資源測量點，每隔15天進(jìn)行一次水位測量。根據(jù)現(xiàn)有資料及1996年、2007年、2014年的地下水資源埋藏深度數(shù)據(jù)，將現(xiàn)有數(shù)據(jù)進(jìn)行最小二乘法擬合，并進(jìn)行數(shù)理統(tǒng)計，發(fā)現(xiàn)這些測量數(shù)據(jù)均服從正態(tài)分布，1996年的測試點中大部分在擬合曲線誤差范圍內(nèi)，僅有個別測量點偏離幅度較大，對于這些點，應(yīng)當(dāng)剔除。3組的測量數(shù)據(jù)中誤差最小的是2014年，運用F1uent軟件對2014年的研究區(qū)域的水位進(jìn)行計算，得出了水位分布云圖，見圖1。

圖1　研究區(qū)域水位分布云圖

從圖1可以看出，藍(lán)色區(qū)域只分布在研究區(qū)的周邊，而研究區(qū)的中央水位呈現(xiàn)淺綠色，說明研究區(qū)域中央水位很低，四周水位較高，形成了顯著的地下水漏斗。對比三組數(shù)據(jù)的漏斗面積，發(fā)現(xiàn)漏斗有生長的趨勢，如不采取合理措施，將導(dǎo)致嚴(yán)重后果。

2　地下水資源預(yù)測

2.1 數(shù)學(xué)模型

地下水資源預(yù)測研究建立在二維非穩(wěn)態(tài)數(shù)學(xué)模型基礎(chǔ)上，通過F1uent軟件對地下水資源分布情況進(jìn)行計算，得出地下水資源的動態(tài)變化。二維非穩(wěn)態(tài)數(shù)學(xué)模型是描述整個研究區(qū)域水資源分布的計算式集合，本文建立的水文地質(zhì)數(shù)學(xué)模型為［6-7］：式中：K—滲透系數(shù)，無量綱；h—含水層厚度，m；n—外法線方向，無量綱；q—出含水層流量，kg/s；S—儲水系數(shù)，無量綱。

二維非穩(wěn)態(tài)數(shù)學(xué)模型在研究過程對一些影響不大的因素進(jìn)行了簡化，降低計算難度，縮短了計算時間。對于含水層，該模型認(rèn)為潛水的埋藏深度隨時間動態(tài)變化，潛水流動方向與地表地形一致，將潛水層當(dāng)作各項異性的二維非穩(wěn)態(tài)流動。

2.2 地下水預(yù)測結(jié)果

利用Visua1MODFLOW軟件對研究區(qū)域地下水資源分布情況進(jìn)行求解。首先根據(jù)GPS信息，繪制研究區(qū)域的數(shù)據(jù)地圖，并依據(jù)此建立有限元分析模型，隨后采用Gambit軟件對地圖進(jìn)行網(wǎng)格劃分，本研究區(qū)域共劃分為3200個計算網(wǎng)格。影響地下水資源分布的因素主要有含水層結(jié)構(gòu)和水文地質(zhì)參數(shù)，這些采用經(jīng)驗系數(shù)。根據(jù)歷史數(shù)據(jù)采用時間列序分析法進(jìn)行水資源預(yù)測，得到了2016年和2024年的地下水等水位線圖，見圖2和圖3。

圖2　2016年的地下水等水位線圖

圖3　2024年的地下水等水位線圖

從圖2和圖3看，2024年研究區(qū)域地下水流場的變化趨勢與2016年類似，呈西北向東南流動。水位最高處位于研究區(qū)域上游，地下水流動方向與地表地形趨勢基本相同。由于研究區(qū)域上游地下水補給量大，因此地下水位相對較高，而下游地下水位較低，并且有進(jìn)一步降低的趨勢。

從2014～2024年間，整個研究區(qū)域共出現(xiàn)了3個顯著的地下水漏斗，且漏斗面積逐年增加，最大漏斗中心點水位竟低于外圍點6m。

對于出現(xiàn)此種狀況的原因，主要有以下幾點。

（1）研究區(qū)域的植物結(jié)構(gòu)出現(xiàn)了明顯變化；

（2）灌區(qū)管理制度變化，由于地下水開采較為方便，符合實時灌溉的需求，因此地下水開采量劇增；

（3）近年來，研究區(qū)域降水量減少，加之上游水庫蓄水，導(dǎo)致下游來流較少，地下水補給量減少；

（4）干渠水量減少，地下水渠系滲漏補給和田間灌溉補給降低。

3　結(jié)語

由于地下水是不可見的，且受地質(zhì)條件影響遠(yuǎn)大于地表水，因此地下水資源預(yù)測比地表水預(yù)測難得多。從國內(nèi)外的研究情況看，一般使用GPS定位和探測儀相結(jié)合的方法對地下水進(jìn)行測量。隨著數(shù)值模擬和水動力學(xué)的發(fā)展，計算流體力學(xué)逐漸成為學(xué)者們的研究熱點。根據(jù)國內(nèi)外地下水資源預(yù)測的研究進(jìn)展，結(jié)合灌區(qū)的地質(zhì)條件，對灌區(qū)地下水資源分布情況進(jìn)行了預(yù)測。并得出以下結(jié)論。

（1）二維非穩(wěn)態(tài)模型的計算精度較高，可以準(zhǔn)確的反映出地下水資源分布情況；

（2）研究區(qū)域上游地下水補給量大，因此地下水位相對較高，而下游地下水位較低，并且有進(jìn)一步降低的趨勢；

（3）給出了地下水分布不均的原因，希望可以為今后地下水資源預(yù)測提供指導(dǎo)。

參考文獻(xiàn)

［1］付宏淵，邱祥，曾鈴，等.基于多變量時間序列（GAR）模型的洞庭湖區(qū)域地下水資源量預(yù)測［J］.長沙理工大學(xué)學(xué)報，2014 （03）：40-47.

［2］邢譯心，鮑新華，吳永東，等.基于Visua1MODFLOW的尚志市水源地地下水資源預(yù)測與開采利用［J］.水電能源科學(xué)，2015 （02）：42-45.

［3］王文科，楊澤元，程東會，等.面向生態(tài)的干旱半干旱地區(qū)區(qū)域地下水資源評價的方法體系［J］.吉林大學(xué)學(xué)報（地球科學(xué)版），2011（01）：159-167.

［4］魏勝利，蘇偉杰.分布式地下水模擬模型MODFLOW介紹與應(yīng)用［J］.水利規(guī)劃與設(shè)計，2011（03）：30-32.

［5］周浩.大小凌河地區(qū)地下水資源可利用量預(yù)測模擬研究［J］.水利技術(shù)監(jiān)督，2015（02）：16-20.

［6］邵景力，趙宗壯，崔亞莉，等.華北平原地下水流模擬及地下水資源評價［J］.資源科學(xué)，2009（03）：361-367.

［7］柳育.淺析崆峒區(qū)北部塬區(qū)地下水資源利用存在的問題及建議［J］.水利規(guī)劃與設(shè)計，2012（04）：26-28.

［8］鄭剛.塔里木河下游沙漠邊緣綠洲地下水特征及其變化分析——以喀拉米吉鎮(zhèn)綠洲為例［J］.水利技術(shù)監(jiān)督，2013（06）：7-12.

作者簡介：呂 敏（1963年—），女，工程師。

收稿日期：2015-06-05

DOI：10.3969 /j.issn.1672-2469.2016.02.010

中圖分類號：P641

文獻(xiàn)標(biāo)識碼：A

文章編號：1672-2469（2016）02-0027-03