利用降深比值求解越補(bǔ)水文地質(zhì)參數(shù)簡化方法

滕 凱

(齊齊哈爾市水務(wù)局，黑龍江省齊齊哈爾市 161006)

文章編號：1006—2610(2015)01—0016—04

利用降深比值求解越補(bǔ)水文地質(zhì)參數(shù)簡化方法

滕 凱

(齊齊哈爾市水務(wù)局，黑龍江省齊齊哈爾市 161006)

越流補(bǔ)給條件下水文地質(zhì)參數(shù)求解涉及3個未知數(shù)的超越方程，無法通過常規(guī)的解析法直接獲解。目前的圖解法及試算法或是成果精度差、或是計算過程繁復(fù)，不便實(shí)際應(yīng)用。采用優(yōu)化擬合方法，在工程適應(yīng)參數(shù)范圍內(nèi)，用較為簡單的函數(shù)完成了由圖表給出的井函數(shù)數(shù)值關(guān)系的替代，并利用水位降深比值關(guān)系經(jīng)整理獲得了僅含1個未知量的表達(dá)式，經(jīng)簡單試算即可完成參數(shù)求解，計算過程簡捷，不依賴圖表，便于實(shí)際工程應(yīng)用。

越流補(bǔ)給；水位降深比值；水文地質(zhì)參數(shù)；優(yōu)化擬合；簡化計算

1 問題的提出

地下水是水資源的重要組成部分，在國民經(jīng)濟(jì)及社會發(fā)展中占據(jù)主導(dǎo)地位。因此，作為地下水資源量評價及含水層污染風(fēng)險評估的主要依據(jù)——水文地質(zhì)參數(shù)，一直以來備受水文地質(zhì)工作者的廣泛關(guān)注。由于構(gòu)造單一各向均質(zhì)的含水層在實(shí)際工程中并不多見，因此，研究具有越流補(bǔ)給條件下的水文地質(zhì)參數(shù)的求解更具實(shí)際意義。因越流補(bǔ)給水文地質(zhì)參數(shù)的求解涉及3個未知數(shù)(即：含水層導(dǎo)水系數(shù)T，m2/h：含水層彈性釋水系數(shù)S；越流因數(shù)B，m)，且屬超越方程求解問題，用常規(guī)的解析法不能直接完成求解。目前，用來進(jìn)行有越流補(bǔ)給條件下水文地質(zhì)參數(shù)求解的方法主要為圖解法[1-2](標(biāo)準(zhǔn)曲線法，拐點(diǎn)法)和試算法[3-4](降深比值法)，由于圖解法在圖形繪制與數(shù)據(jù)讀取上存在人為誤差，且實(shí)際工作受圖表束縛；試算法由于重復(fù)計算及查表插值，求解過程比較繁復(fù)。隨著計算機(jī)技術(shù)的進(jìn)一步普及，肖常來、李培月等先后提出了全程曲線擬合法[5]及數(shù)值模擬法[6]，但因全程曲線擬合法涉及建立目標(biāo)函數(shù)及約束條件，計算中必需對井流公式的解析解進(jìn)行近似化；數(shù)值模擬法需綜合考慮含水層結(jié)構(gòu)、地下水動態(tài)、水文信息等諸多因素，涉及12個方面較為復(fù)雜的技術(shù)工作，而且調(diào)參過程工作量大，計算時間長。由于這兩種方法均需借助相應(yīng)的計算軟件及程序，又要處理上述較為繁雜的技術(shù)問題，對應(yīng)用者的基本技能要求較高，實(shí)際推廣應(yīng)用尚存在問題。因此，研究提出一種計算方法簡單、求解成果精度可靠的越流補(bǔ)給條件下水文地質(zhì)參數(shù)的求解方法，是目前該領(lǐng)域十分關(guān)切和正在探索的主要內(nèi)容之一。本文采用優(yōu)化擬合方法，通過對井函數(shù)標(biāo)準(zhǔn)曲線的線型分析，借助現(xiàn)有表格給出的井函數(shù)數(shù)值對應(yīng)關(guān)系，通過逐次逼近擬合計算，獲得了井流函數(shù)的簡化替代公式，并利用該公式依據(jù)抽水降深資料直接完成了相關(guān)水文地質(zhì)參數(shù)的求解計算，計算方法簡單直接，可在實(shí)際工程中推廣應(yīng)用。

2 公式的建立

2.1 基本計算公式

在抽水影響范圍內(nèi)，如含水層為多層、各層近似水平、均質(zhì)等厚、各向同性并側(cè)向無限延伸，且有上復(fù)或下臥的弱透水層在越流含水層抽水時能產(chǎn)生越流補(bǔ)給，此時的地下水滲流運(yùn)動即為第一類越流補(bǔ)給，見圖1所示。其越流系統(tǒng)井流問題的基本求解方程為：

圖1 越補(bǔ)含水層示意圖

(1)

其 中：

(2)

(3)

(4)

(5)

當(dāng)抽水主孔或觀測孔在t1、t2和t3時刻的水位降深分別為s1、s2和s3，則由式(1)可得到以下方程：

(6)

(7)

(8)

(9)

(10)

式中：s12、s23分別為與抽水歷時t1、t2和t2、t3相對應(yīng)的降深比值。

利用式(9)、(10)即可通過試算法求得參數(shù)α及B值，進(jìn)而完成其他相關(guān)水文參數(shù)的計算。

2.2 簡化公式的建立

圖關(guān)系曲線圖

經(jīng)逐次逼近擬合[8]可得到式(2)的最優(yōu)替代函數(shù)：

(11)

其 中：A=1.619u1.2-3.3u0.8+1.944u0.4+0.187

(12)C=4.287u0.7-1.411u1.05-3.941u0.35+1.034

(13)

D=exp(0.415u1.2+0.568u0.48+2.055u0.24-1.862)

(14)

式中：A、C、D均為與參數(shù)α及t有關(guān)的中間變量。

將式(11)～(14)分別代入式(9)、(10)中，經(jīng)進(jìn)一步整理即可獲得求解參數(shù)α的計算公式：

(15)

其 中:

(16)

(17)

(18)

(19)

(20)

(21)

再由式(22)求得彈性釋水系數(shù)S為：

(22)

2.3 精度分析

圖3 式(11)替代式(2)誤差包絡(luò)線圖

(i=1,2,3 ,…,n)

(23)

式中：zi為擬合替代相對誤差，%；i為擬合計算中所選取的第i個數(shù)據(jù)。

由圖3可見，在工程適用參數(shù)范圍內(nèi)，利用本文近似井函數(shù)式(11)替代式(2)的最大相對正、負(fù)誤差分別為4.58%和-4.64%，完全可以滿足實(shí)際工程的設(shè)計精度要求。

3 算 例

選文獻(xiàn)[4]例3-6：某河階地，上部為潛水層，其下為2m厚的弱透水亞砂土，再下為1.5m厚的中、粗砂層(承壓)。水源地提取承壓水，為求取參數(shù)以T32號孔做非穩(wěn)定流抽水試驗(yàn)，距它197m處有T31號觀測孔，抽水流量Q=69.1m3/h，觀測孔水位降深值如表1。

表2 觀測孔水位降深資料表

解：根據(jù)表1繪制s～t關(guān)系曲線(見圖4)，并在s～t關(guān)系曲線上取：t1=25min,s1=0.383m;t2=90min,s2=0.645m;t3=300min,s3=0.768m。

圖4 s～t關(guān)系曲線圖

依據(jù)上述已知參數(shù)，并設(shè)α=2 894 500m2/d，由式(12)～(14)及式(16)～(18)即可分別求得：

將上述參數(shù)代入式(15)可求得式左、右分別為：

當(dāng)在s～t關(guān)系曲線上取另一組數(shù)據(jù)時，即：

t1=25min, s1=0.383m

t2=240min, s2=0.757m

t3=480min, s3=0.794m

由上述計算成果可見，雖然在s～t關(guān)系曲線上所選取的2組觀測數(shù)據(jù)變化較大，但所求水文地質(zhì)參數(shù)卻極為接近，說明本文方法具有較好的可靠性及實(shí)用性。

4 結(jié) 語

(1) 在工程適用范圍內(nèi)，本文提出的井函數(shù)替代公式表達(dá)形式簡單、擬合精度高。利用該公式進(jìn)行第一類越流補(bǔ)給水文地質(zhì)參數(shù)求解不受圖表束縛，避免了圖表取值的人為誤差，而且計算過程簡捷，便于實(shí)際推廣應(yīng)用。

(2) 為有效提高計算成果的精度，考慮水位降深觀測成果可能存在誤差，在利用本文方法計算前，可根據(jù)s～t的觀測成果完成s～t曲線繪制，并通過s～t曲線修正相關(guān)水位降深值 。

(3) 利用s～t曲線進(jìn)行(t1,s1)、(t2,s2)及(t3,s3)3個點(diǎn)選取時，為使各點(diǎn)具有更好的代表性，應(yīng)盡量將點(diǎn)(t1,s1)選在s～t曲線斜率較大區(qū)段(抽水前期區(qū)段)，點(diǎn)(t2,s2)應(yīng)選在s～t曲線斜率變化較大區(qū)段(曲率半徑較小區(qū)段)，點(diǎn)(t3,s3)應(yīng)選在s～t曲線斜率較小區(qū)段(抽水后期時段)。

(4) 采用與本文類似的方法也可完成第二類越流補(bǔ)給水文地質(zhì)參數(shù)的簡化求解計算問題，本文不詳述。

[1] 薛禹群,吳吉春.地下水動力學(xué)[M].3版.北京：地質(zhì)出版社，2010.

[2] 《供水水文地質(zhì)手冊》編寫組. 供水水文地質(zhì)手冊[M].北京：地質(zhì)出版社，1990.

[3] 劉兆昌,朱昆. 供水水文地質(zhì)[M].北京：地質(zhì)出版社，1979.

[4] 苑蓮菊,李振栓,楊展,等.工程滲流力學(xué)及應(yīng)用[M].北京：中國建材工業(yè)出版社，2001,7.

[5] 肖長來,梁秀娟,崔建銘,等.確定含水層參數(shù)的全程曲線擬合法[J].吉林大學(xué)學(xué)報(地球科學(xué)版)：2005，35(6)：751-755.

[6] 李培月.非穩(wěn)定流抽水試驗(yàn)確定越流承壓含水層水文地質(zhì)參數(shù)方法對比研究[D].西安：長安大學(xué)，2011.

[7] 王慧文.偏最小二乘回歸法及其應(yīng)用[M].北京：國防工業(yè)出版社,1999.

[8] 閻鳳文.測量數(shù)據(jù)處理方法[M].北京：原子能出版社,1988.

Simplified Method for Calculation of Hydrogeological Parameters of Leakage Recharge by Water Drawdown Ratio

TENG Kai

(Qiqihar Water Authority, Qiqihar 161006，China)

Three unknown numbers are involved to derive the hydrogeological parameters in condition of leakage recharge. So, the conventional method cannot be applied to directly get them. Both graphic method and trial method are not practical because of the poor precision of the results or the complicated calculation process. In range of the applicable parameters, by the optimized fitting method, the simple function is applied instead of the value relationship of the well function given by graph. Furthermore, an expression with one unknown number is derived by application of the water drawdown ratio. Through the simple trial calculation, parameters are derived. Its calculation process is simple, without relying on graph, and convenient to apply in practice.

leakage recharge; water drawdown ratio; hydrogeological parameter; optimized fitting; simplified calculation

2014-06-30

滕凱(1957- )，男，黑龍江齊齊哈爾市人，高級工程師，齊齊哈爾市水務(wù)局總工程師，主要從事防災(zāi)減災(zāi)及水利工程優(yōu)化設(shè)計研究.

P641.2

A

10.3969/j.issn.1006-2610.2015.01.005