用有限差分模型反求高升水源各項補排量

湯玉福

(遼寧省營口水文局，遼寧 營口 115000)

?

用有限差分模型反求高升水源各項補排量

湯玉福

(遼寧省營口水文局，遼寧 營口 115000)

摘要：根據(jù)含水層的結(jié)構(gòu)變化、富水性特征、地下水等水位線的疏密程度，將計算區(qū)按三角形單元進行剖分；在計算區(qū)域剖分和時間剖分基礎(chǔ)上，按均衡原理，建立不規(guī)則網(wǎng)格的有限差分模型；利用計算區(qū)的長期觀測資料，對模型進行識別和驗證，并反求水文地質(zhì)參數(shù)和各項補排量。

關(guān)鍵詞：高升水源地；地質(zhì)構(gòu)造；差分模型；識別；驗證；水文地質(zhì)參數(shù)

1水源地現(xiàn)狀

高升水源地處下遼河中部平原,位于市區(qū)東北部高升鎮(zhèn)北,繞陽河?xùn)|岸一帶，是盤錦市區(qū)主要的供水水源，日取水量為6萬m3/d，位于盤錦市東北方向，距市區(qū)20km，水源地分布在高升鎮(zhèn)、喜彬鄉(xiāng)、大荒農(nóng)場及遼河油田高升采區(qū)，總面積217.5km2。

1.1地形地貌

水源地所在區(qū)域內(nèi)地貌形態(tài)的形成是以地殼升降運動為主，也是新構(gòu)造運動、基底構(gòu)造斷裂、古氣候變化、河流等內(nèi)外動力地質(zhì)作用的綜合影響的結(jié)果。由于地殼不斷沉降，故形成了以堆積地形為主的地貌形態(tài)，呈沖積平原的地貌景觀。地勢由東北向西南微傾斜，地面標(biāo)高為4～8m，其上分布有零星固定和半固定沙丘，表面分布亞黏土和粉細砂。

1.2地質(zhì)構(gòu)造

水源地所在區(qū)域在地質(zhì)構(gòu)造上屬于下遼河凹陷一部分，侏羅紀以前與華北地臺為統(tǒng)一整體，三疊紀以后，由于構(gòu)造運動形成了一些零星分割的盆地，隨著構(gòu)造運動的加強，形成控制盆地形成和發(fā)展的多字型斷裂構(gòu)造，同時伴隨有劇烈的巖漿活動。喜山運動時，盆地發(fā)生強烈分異。早第三紀，在多字型斷裂構(gòu)造的控制下，隨著盆地大幅度斷裂式下沉，發(fā)生強烈分異作用，形成一系列緊密相間的隆起和拗陷，水源所在區(qū)域就位于其中的盤山拗陷。

1.3氣象水文

水源地所在區(qū)域?qū)倥瘻貛Т箨懶园霛駶櫦撅L(fēng)氣候，四季分明，春季少雨多風(fēng)，夏季炎熱多雨，秋季氣高涼爽，冬季寒冷少雪。多年平均降水量為660.2 mm，但年內(nèi)時空分布不均，春季(3—5月)平均降水量96.5 mm，占全年降水量的15.5%；夏季(6—8月)平均降水量392.1 mm，占全年降水量的62.9%；秋季(9—11月)平均降水量121.0 mm，占全年降水量的19.4%；冬季(12—2月)平均降水量13.6 mm，占全年降水量的2.2%。年降水量最大為934.5 mm(1986年)，年降水量最小為363.4 mm(1989年)多年平均蒸發(fā)量為1705.0 mm，多年平均相對濕度為66%。水源地干流總長283km，流域面積9 950 km2，區(qū)域內(nèi)設(shè)有王回窩堡水文站，多年平均流量為8.47 m3/s。

1.4水文地質(zhì)及地下水的補給、徑流和排泄條件

地下水主要賦存于第四系松散層中，含水層厚度為70～122 m，巖性為中砂、細砂、間夾薄層亞黏土但不連續(xù)，故多層含水層的水力聯(lián)系較為密切，其底板為下更新統(tǒng)下遼河組淤泥質(zhì)亞黏土，埋深為70～130m。此外在第四系之下，還廣泛分布第三系孔隙-裂隙水。水源區(qū)域內(nèi)地下水以垂向補給為主，側(cè)向補給為輔，前者以大氣降雨面狀入滲補給和河渠滲漏線狀補給為主，以灌溉滲漏的局部面狀補給和開采條件下的激發(fā)補給為輔。地下水總體流向是由北東向南西流，然后向遼東灣排泄，由于水力坡度較小，因此徑流條件極差。區(qū)域內(nèi)地勢低平坦，地下水位埋深淺(1～2m)，局部直接出露地表，形成沼澤。因此，蒸發(fā)是地下水主要排泄途徑。其次是向河流排泄。

2有限差分法的基本思想

有限差分法是指用泰勒技術(shù)展開式將變量的導(dǎo)數(shù)寫成變量，在不同時間或空間點值的差分形式的方法，按時間步長和空間步長將時間和空間區(qū)域剖分成若干網(wǎng)格，用未知函數(shù)在網(wǎng)格結(jié)(節(jié))點上的值所構(gòu)成的差分近似代替所用偏微分方程中出現(xiàn)的各階導(dǎo)數(shù)，從而把表示變量連續(xù)變化關(guān)系的偏微分方程離散為有限個代數(shù)方程，然后解此線性代數(shù)方程組，以求出溶質(zhì)在各網(wǎng)格結(jié)(節(jié))點上不同時刻的濃度。有限差分法一般分以下幾個步驟：①剖分滲流區(qū)，確定離散點；②建立水動力彌散問題的差分方程組；③求解差分方程組。

3模型的建立

首先根據(jù)含水層的結(jié)構(gòu)變化、富水性特征、地下水等水位線的疏密程度，將計算區(qū)按三角形單元進行剖分，共剖分為230個三角形單元，含136個節(jié)點，其中一類邊界點10個，二類邊界點30個，內(nèi)節(jié)點96個，待求水位節(jié)點126個，然后根據(jù)含水層的結(jié)構(gòu)和水文地質(zhì)參數(shù)以及富水性，將全區(qū)域劃分為15個水文地質(zhì)參數(shù)區(qū)，參數(shù)分區(qū)及剖分圖見圖1。在計算區(qū)域剖分和時間剖分的基礎(chǔ)上，按均衡原理確定不規(guī)則網(wǎng)格的有限差分模型，公式為：

(1)

為便于求解，采用超松弛迭代法，其通式為：

則：hin+1=hin+ω△hi

(2)

式中：hin為第n次迭代的水位值，m；n為迭代次數(shù)；ω為超松弛因子，采用0.85。

圖1　參數(shù)分區(qū)及剖分圖

4模型的識別與驗證

利用計算區(qū)的長期觀測資料，對已建立的數(shù)學(xué)模型進行識別和驗證，目的在于檢驗數(shù)學(xué)模型的適用性，并反求水文地質(zhì)參數(shù)和各項補給量。

首先將各項補排量置零，運行模型求得的水位變化與實測水位的變化進行比較，如果兩者基本一致，說明已建立的潛水?dāng)?shù)學(xué)模型符合實際，同時檢驗了邊界性質(zhì)。求參方法是采用正演擬合水位，間接求水文地質(zhì)參數(shù)，然后用反演擬合水量，間接檢驗參數(shù)。根據(jù)各時段的補給量、排泄量和參數(shù)初值，運行模型求解水位。如模擬計算水位和實測對比水位相差較大，則根據(jù)參數(shù)的變化范圍和參數(shù)與水位的相關(guān)性，再試一組參數(shù)值，直到水位擬合較好為止。這時選用的一組參數(shù)即為所求參數(shù)。然后在不改變參數(shù)的情況下，以實測水位為依據(jù)，反求補給量，如與實測補給量相近，證明該組參數(shù)比較可靠。

5成果比較

5.1水位擬合情況

表1　水位擬合精度統(tǒng)計表

由表1可以看出，計算水位和實測水位之間絕對誤差<0.5m的節(jié)點占計算總節(jié)點數(shù)的85%以上。而絕對誤差0.9~1.0m的節(jié)點，第1時段有3個，第2時段有4個，分別占計算節(jié)點數(shù)的2%和3%。由此可見擬合較好精度較高，說明建立的數(shù)學(xué)模型可靠。

5.2水文地質(zhì)參數(shù)

表2　水文地質(zhì)參數(shù)表

由表2可見，所求的水文地質(zhì)參數(shù)基本符合評價區(qū)的含水層巖性情況和變化規(guī)律，由北向南滲透系數(shù)漸減。并與抽水試驗求得的參數(shù)接近。

5.3區(qū)域內(nèi)各項資源量計算

表3　區(qū)域內(nèi)各項資源量計算

通過數(shù)值模擬求得計算區(qū)的各項資源量與其他常規(guī)方法計算結(jié)果相近，總補給量比總排泄量多254.56m3/a，地下水處于正平衡，地下水位上升，與2009年的地下水動態(tài)資料吻合。

參考文獻：

[1]李曉昭,廖作鴻.有限差分方法在期權(quán)定價中的應(yīng)用[J].科技情報開發(fā)與經(jīng)濟,2004,14(04)：83-84.

Using Finite Difference Model to Reverse Various Complementary Displacement of Gaosheng Water Source

TANG Yu-fu

(Yingkou Hydrological Bureau Liaoning Province，Yingkou 115000，China)

Abstract：According to the structural changes of aquifer，water rich characteristics and density of groundwater level line，the calculation area was divided in line with triangular elements.On the basis of calculation area and time division，a finite difference model of irregular net was established based on the balance principle.The model was identified and confirmed using the long-term observation date in calculation area and to reverse the hydrogeological parameters and various complementary displacements.

Key words：Gaosheng water head site；geological structure；difference model；identification；confirmation；hydrogeological parameter

[作者簡介]湯玉福(1962-)，男，云南大理人，高級工程師，從事水文水資源管理工作。

[收稿日期]2015-06-28

中圖分類號：P641

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

文章編號：1007-7596(2015)10-0025-03