陳凌顏， 滕 凱

(1.齊齊哈爾市河道管理處，黑龍江 齊齊哈爾 161006； 2.齊齊哈爾市水務局,黑龍江 齊齊哈爾 161006)

?

利用非穩(wěn)定流求解潛水水文地質參數的降深比值6點法

陳凌顏1， 滕 凱2

(1.齊齊哈爾市河道管理處，黑龍江 齊齊哈爾 161006； 2.齊齊哈爾市水務局,黑龍江 齊齊哈爾 161006)

利用Boulton井流公式求解潛水非穩(wěn)定流條件下的水文地質參數，是涉及4個未知數的超越方程,無法通過常規(guī)的解析法直接獲解.傳統(tǒng)的圖解法不但受圖表束縛,而且求解過程繁復、成果精度存在人為誤差.文中采用優(yōu)化擬合方法,在工程適用參數范圍內,用較為簡單的函數完成了由圖表給出的井函數數值關系的替代,并利用降深比值關系,通過選取水位降深曲線上的6個點,經整理獲得了僅含1個未知量的表達式,經簡單試算即可完成參數求解,結果更接近工程實際,且計算過程簡捷,不依賴圖表,便于實際工程應用.

越流補給;水位降深比值;水文地質參數;優(yōu)化擬合;簡化計算

潛水是地下水資源的重要組成部分,具有埋藏深度小、補給條件好、開發(fā)利用成本低等優(yōu)點,在地下水資源的開發(fā)利用中占有重要地位.因此,研究地下潛水水文地質參數的獲取方法,對定量評價地下潛水水資源量及制訂防污染計劃意義重大.因利用Boulton井流公式求解潛水非穩(wěn)定流條件下水文地質參數為涉及4個未知數(即:含水層導水系數T,m2/h；含水層彈性釋水系數S;給水度μ;越流因數D,m)的超越方程,無法通過解析法直接完成參數求解.目前,工程上普遍采用圖解法[1],即降深-時間量板法[2]、標準曲線比對法和直線圖解法[3-4].而現有利用計算機軟件實現的智能算法均派生于Theis公式[5-7],并不能直接應用于Boulton井流公式.由于標準曲線比對法需分別完成對抽水初期及后期降深曲線與標準曲線的手動擬合,實際工作不但受圖表束縛,而且求解成果因受標準曲線的密度及比對過程中人為因素的影響,其精度難免存在誤差[8];直線圖解法是在假定抽水后期潛水滯后重力給水對水位降深沒有影響的情況下,降深與時間的對數可用直線關系表示,而實際工程并非完全如此,另一方面,由于利用所取直線的斜率確定的切線位置不唯一,所得參數也很難一致,有時甚至差距很大[9].因此,研究提出一種計算方法簡單、求解成果精度可靠的潛水非穩(wěn)定流情況下水文地質參數的獲取方法,對進一步提高求解潛水含水層水文地質參數的工作效率和成果精度具有重要意義.

筆者采用優(yōu)化擬合方法[10-12],通過對潛水非穩(wěn)定流井函數標準曲線的線型分析,借助現有表格給出的井函數數值對應關系,通過逐次逼近擬合計算,獲得了潛水非穩(wěn)定流井流函數的簡化替代公式.并利用該公式通過降深比值關系,經數學推導獲得了通過選取抽水降深關系曲線上的6個點直接完成相關水文地質參數求解的計算公式,所得參數結果更接近工程實際,計算方法簡單直接,可在實際工程中推廣應用.

1 公式的建立

1.1 Boulton模型

根據Boulton理論,在潛水含水層中一完整井以定流量Q進行連續(xù)抽水，如圖1所示,其含水層的水位變化可用如下公式表示[13]:

(1)

(2)

(3)

(4)

(5)

(6)

圖1 潛水含水層示意圖

1.2 降深比值求解法

由于利用潛水非穩(wěn)定流抽水資料確定水文地質參數需分別根據抽水初期和后期的計算結果才能最終確定,因此,下面分別就抽水初期和后期的參數求解方法進行研究.

1.2.1 抽水初期

在抽水初期(抽水主孔或觀測孔)的降深曲線上分別選取(t1,s1)、(t2,s2)和(t3,s3),則由式(1)可得到以下方程組:

(7)

(8)

(9)

(10)

(11)

式中:s12，s23分別為與抽水歷時t1，t2和t2，t3相對應的降深比值.

圖2 W-uα關系曲線

由圖2可見,W-uα具有較好的指數函數關系.依據函數曲線的這一特征,經過對多組具有類似曲線形式的備選函數的擬合比選,并以標準剩余差最小為目標函數[14],即:

(12)

(13)

其中:

(14)

(15)

(16)

式中:A，B，C均為與參數a及t有關的中間變量.

將式(13)—(16)分別代入式(10)、(11)中,經整理即可獲得求解參數a的計算公式，即

(17)

其中:

A21=A2-A1,A32=A3-A2,

(18)

B21=B2-B1,B32=B3-B2,

(19)

C21=C2-C1,C32=C3-C2.

(20)

(21)

(22)

(23)

再由式(5)求得彈性釋水系數S，即

(24)

1.2.2 抽水后期

采用與利用抽水降深初期資料求解水文地質參數相同的方法,即可完成利用抽水降深后期資料求解水文地質參數的計算.

(25)

其中:

(26)

(27)

(28)

式中:A′，B′，C′均為與參數ay及t有關的中間變量.

(29)

其中:

(30)

(31)

(32)

(33)

(34)

(35)

再由式(6)求得給水度μ，即

(36)

進而可求得潛水含水層的總給水度μ0及給水度特征比值η,即

(37)

2 精度分析

(38)

式中:zi為擬合替代相對誤差，%;i為擬合計算中所選取的第i個數據精度比較點.

圖3 式(13)替代誤差包絡線

圖4 式(25)替代誤差包絡線

由圖3和圖4可見,在工程適用參數范圍內,利用本文近似井函數式(13)(或式(25))替代文獻[2]附表2中理論值的最大相對正、負誤差分別為5.21%和-5.37%,完全可以滿足實際工程的設計精度要求.

3 算 例

選文獻[2]算例:在某地第四紀含水層中進行抽水試驗,含水層厚7.75 m,由含礫石的中砂組成.隔水底板為砂巖.抽水井為完整井,井徑為0.32 m,抽水持續(xù)時間為48 h 50 min,抽水流量Q=53 m3/h,抽水井至觀測孔的距離r=10 m.抽水試驗資料見表1,試計算該含水層水文地質參數.

根據表1繪制s-t關系曲線，如圖5所示,并對曲線進行認真修正,剔除個別點因測量或設備因素產生的誤差影響.

分別在s-t關系曲線的初期段和后期段上選取:t1=0.75 min,s1=0.137 m;t2=2.00 min,s2=0.185 m;t3=10.00 min,s3=0.221 m;t4=670.00 min,s4=0.373 m;t5=1 450.00 min,s5=0.429 m;t6=2 650.00 min,s6=0.475 m.

表1 觀測孔實測資料

對抽水初期階段,依據上述已知參數,并設a=695 100 m2/d,由式(14)—(16)及式(18)—(20)即可分別求得:A1=-1.912 389,A2=-1.951 335,A3=-1.971 105,B1=1.179 148,B2=0.638 819,B3=-0.215 999,C1=0.594 548,C2=1.130 126,C3=1.673 282,A21=-0.038 946,A32=-0.019 771,B21=-0.540 328,B32=-0.854 818,C21=0.535 577,C32=0.543 156,s12=0.740 541,s23=0.837 104.

圖5 s-t關系曲線

將上述參數代入式(17)可求得等號左、右兩端分別為:

由上述計算結果可得該潛水含水層導水系數T的平均值為

總給水度及給水度特征比值為

μ0=μ+S=0.060 43;

4 結 語

依據抽水降深比值關系,對以圖表形式給出的潛水非穩(wěn)定流井函數值進行函數替代,實現了利用s-t曲線上的6個點完成4個水文參數的求解目標,與現行的其他方法比較,具有以下特點.

1)公式的表達形式比較簡單、計算過程簡捷直觀.由于計算所利用的6個點均取自實際抽水降深曲線,因此,所獲結果與其他方法比較，更接近工程實際.

2)本文求解方法不受圖表束縛及標準曲線族密度影響,避免了圖表取值的人為誤差及反復進行擬合曲線比對的舍取之繁,便于實際推廣應用.

3)為有效提高計算結果的精度,考慮水位降深觀測成果可能存在誤差,在利用本文方法計算前,可根據s-t的觀測成果完成s-t曲線繪制,并通過s-t曲線修正相關水位降深值s.

4)利用s-t曲線進行(t1，s1)—(t6，s6)的6個點選取時,為使各點具有更好的代表性,應盡量將點(t1，s1)—(t3，s3)的3個點選在s-t曲線斜率較大區(qū)段(抽水前期區(qū)段),將(t4，s4)—(t6，s6)的3個點選在s-t曲線的后區(qū)段(抽水后期時段).

5)采用與本文類似的方法也可完成類似情況下的水文地質參數的簡化求解計算問題,本文不詳述.

[1]段乃金.解析法求潛水含水層水文地質參數的對比分析[D].西安:長安大學,2010.

[2]供水水文地質手冊編寫組.供水水文地質手冊[M].北京:地質出版社,1990.

[3]陳崇希,林敏.地下水動力學[M].武漢:中國地質大學出版社,1999.

[4]張元禧,施鑫源.地下水水文學[M].北京:水利水電出版社,1998.

[5]高瑞忠,朝倫巴根,賈德彬,等．基于非穩(wěn)定流抽水試驗的BP神經網絡確定含水層參數研究[J]．沈陽農業(yè)大學學報,2004,35(5-6)：510-512.

[6]魏連偉,鄒景力,張建立,等．模擬退火算法反演水文地質參數算例研究[J]．吉林大學學報:地球科學版,2004,34(4):612-614．

[7]郭建青,李彥,王洪勝,等．確定含水層參數的混沌序列優(yōu)化算法[J]．中國農村水利水電,2006(12):26-29．

[8]李偉,趙燕容,朱旭芬,等.自動配線法求水文地質參數的優(yōu)化路徑法研究[J].勘察科學技術,2013(2):6-10.

[9]李炳森.利用抽水試驗降深比值求水文地質參數[J].煤炭技術,2005,24(8):74-75.

[10]滕凱.標準門洞形過水斷面臨界水深的簡化計算[J]．華北水利水電學院學報,2012,33(5):1-3.

[11]滕凱.馬蹄形斷面隧洞正常水深的簡化計算法[J]．華北水利水電學院學報,2013,34(5):31-33，76.

[12]滕凱,張麗偉.弧底梯形渠道收縮水深計算方法的簡化[J]．華北水利水電學院學報,2013,34(3):55-58.

[13]薛禹群,吳吉春.地下水動力學[M].3版.北京:地質出版社,2010.

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

[15]閻鳳文.測量數據處理方法[M].北京:原子能出版社,1988.

(責任編輯：喬翠平)

Six Point Method of Drawdown Ratio for Solving Hydrogeological Parameters of Phreatic Water under Unsteady Flow Condition

CHEN Lingyan1, TENG Kai2

(1.Tsitsihar River Management Office, Tsitsihar 161006, China; 2.Tsitsihar Water Authority, Tsitsihar 161006, China)

During solving hydrogeological parameters of phreatic water under unsteady flow Condition by Boulton well-flow formulas, the transcendental equations containing four unknown numbers are involved, which are not directly solved by the conventional analytical method. Not only the traditional descriptive geometry solution is bounded by charts, but also the solution procedure is complicated, and the precision of results is affected by personality errors. In this article, the optimal fitting method was used to finish the replacement of magnitude relation of well function given by charts with simpler functions, the expressions containing one unknown number were obtained through utilizing the relation of drawdown ratio and selecting six points in drawdown curves, the parameter solving was completed through simple and tentative calculation, and the solutions were much more approximate to the practice, the solution procedure was simple, did not rely on charts, and was convenient for practical application.

leakage recharge; drawdown ratio; hydrogeological parameters; optimum fitting; simplified calculation

2014-11-08

齊齊哈爾市科技攻關項目(NYGG201301).

陳凌顏(1982—)，女，黑龍江肇源人，工程師，主要從事工程建設管理方面的研究. 滕 凱(1957—)，男，黑龍江齊齊哈爾人，高級工程師，主要從事水利防災減災及工程優(yōu)化設計方面的研究.

10.3969/j.issn.1002-5634.2015.01.002

TV211.12

A

1002-5634(2015)01-0006-06