任廣軍1 張勇2
（1.山東省魯南地質工程勘察院，兗州272000；2.山東省地礦工程集團有限公司，濟南250013）
作者簡介：任廣軍（1972—），男，工程師，主要從事水文地質、環境地質等。
摘要：本文利用非穩定流抽水試驗資料，采用改進的十進制遺傳算法在計算機上自動優選含水層水文地質參數。該方法同傳統上使用的配線法相比較，具有節省時間，減少人工配線誤差，所求參數逼真，且能對壹些線性、非線性問題求解，具有很高的推廣和應用價值。
關鍵詞：遺傳算法；隨機模擬；含水層；水文地質參數；優選
0 引言
利用改進的十進制遺傳算法，根據抽水試驗資料來認識水文地質條件、反求水文地質參數是水文地質計算中的基本問題。具體地講，在探明含水層範圍、類型的基礎上，建立描述該含水層水流運動模型，利用抽水試驗過程中的地下水位變化過程資料來確定水文地質參數。
雖然非穩定抽水試驗公式適用條件非常苛刻，但能反映出含水層非穩定流的壹些基本特點，還可運用疊加原理解決某些比較復雜的非穩定流問題。此外，作為檢驗數值方法精確性的重要依據，具有廣泛應用和發展前景。
目前，由於非穩定流抽水試驗確定水文地質參數的具體實現方法主要有人工配線法或以計算輔助的配線法，但這種方法的效果好壞完全取決於肉眼觀察，帶有很大的主觀性。本文作者選取了壹些典型實例，采用遺傳算法建立了壹種計算機全自動求參的全局優選法，通過與人工配線分析比較，確定本方法計算機求參的高精度與高可靠性。
求取參數是通過實測結果與模型計算結果的最佳擬合（仿真）程度來實現的，參數的精確程度在很大程度上取決於實測資料的精度。
1 遺傳算法介紹
生物的進化是壹個奇妙的優化過程，它通過選擇淘汰，突然變異，基因遺傳等規律產生適應環境變化的優良物種。遺傳算法是根據生物進化思想而啟發得出的壹種全局優化算法。
遺傳算法的概念最早是由Bagley J.D在1967年提出的；而開始遺傳算法的理論和方法的系統性研究的是1975年，這壹開創性工作是由Michigan大學的J.H.Holland所實行。當時，其主要目的是說明自然和人工系統的自適應過程。
遺傳算法簡稱GA（Genetic Algorithm），在本質上是壹種不依賴具體問題的直接搜索方法。遺傳算法在模式識別、神經網絡、圖像處理、機器學習、工業優化控制、自適應控制、生物科學、社會科學等方面都得到應用。在人工智能研究中，現在人們認為“遺傳算法、自適應系統、細胞自動機、混沌理論與人工智能壹樣，都是對今後十年的計算技術有重大影響的關鍵技術”。
2 目標函數的確定
通過綜合考慮計算程序的運算時間、速度以及含水層的類型，確立利用抽水實測資料和計算資料的擬合程度為目標函數。其計算公式為：
山東省環境地質文集
式中：s實測為實測抽水試驗觀測孔的降深；s計算為計算抽水試驗觀測孔的降深；NT為計算時段。
3 計算實例及結果分析
3.1 承壓含水層地下水降深公式
承壓含水層地下水降深公式為：
山東省環境地質文集
式中：S為以固定流量Q抽水時與抽水井距離為r處任壹時間的水位降深（m）；T為導水系數；Q為湧水量；W（u）為井函數，是壹個指數積分函數：
山東省環境地質文集
式中：u為井函數的自變量， 其中s為承壓含水層的儲存系數；T為含水層的導水系數；t為時間。
例1：某地區進行非穩定流抽水試驗。區域地層剖面是：地表下18～25 m是由含礫粗砂層組成的含水層，其底板由粘土質沈積物組成，18 m以上是粘土、泥炭層。抽水井的過濾器安裝在含水層的整個厚度上。觀測孔距抽水井30m，觀測資料如表1所示。主井作定流量抽水，Q＝788m3/d，抽水接近14小時。試根據觀測資料求取水文地質參數。
（1）lgS-lgt配線法所求參數：T配線＝439m2/d，s配線＝1.694×10-4；
（2）S-lgt直線圖解法所求參數：T配線＝450.7m2/d，s配線＝1.392×10-4；
（3）計算機所求參數：T＝383.0088m2/d，s＝ 1.78×10-4。
為更直觀地說明上述所求參數的可靠性，由上述參數所求計算降深與實測降深進行比較（圖1）。通過比較，進壹步確定了計算機求參的高精度與穩定性。承壓含水層配線參數與優選參數比較分析：T配線＝439m2/d，s配線＝0.0001694；T計算＝383.0088m2/d，s計算＝0.0001780。
表1 遺傳算法計算水位降深與實測水位降深結果表
圖1 計算降深與實測降深比較圖
3.2 在有越流補給的承壓含水層地下水降深公式
在有越流補給的承壓含水層地下水降深公式為：
山東省環境地質文集
式中：u同（3）式； 為越流井函數，本文中 采用數值積分：
山東省環境地質文集
例2：有壹無限分布的承壓含水層，厚度20m，其底部為絕對隔水的粘土層；上部為弱透水的亞砂土層，厚2m；弱透水層之上為潛水含水層。在承壓含水層中有壹完整抽水井，抽水時的穩定流量Q＝5530m3/d。距抽水井r＝17.34m處有壹觀測孔據觀測知，在抽水過程中上部潛水的水位不變。抽水層的水位降深值載於表2，試計算含水層水文地質參數。
（1）lgS-lgt配線法所求參數：T配線＝853.50m2/d，s配線＝4.20×10-4；B配線＝568.50m；
（2）lgS-lgt配線法所求參數：T計算＝817.19m2/d，s計算＝4.31×10-4；B計算＝482.80m。
為更直觀地說明上述所求參數的可靠性，由上述參數所求計算降深與實測降深進行比較（圖2）。通過比較，進壹步確定了計算機求參的高精度與穩定性。有越流時承壓含水層優選參數誤差分析：T配線＝853.50m2/d，s配線＝0.00042，B配線＝568.50m；T計算＝817.1950m2/d，s計算＝0.00043103，B計算＝482.798m。
表2 遺傳算法計算水位降深與實測水位降深結果表
續表
圖2 計算降深與實測降深比較圖
3.3 考慮有滯後補給的潛水含水層地下水降深公式
根據博爾頓理論，潛水含水層地下水降深公式計算公式可分為抽水前期、抽水中期和抽水後期。參數優選主要根據抽水前期和抽水後期的資料擬合而得：
山東省環境地質文集
其中D為疏幹因子。
抽水前期計算公式：
抽水後期計算公式：同（2）式。
4 結論及不足之處
4.1 結論
通過上述實例計算結果表明：計算結果同人工加以計算機輔助配線法相比較，其計算水文地質參數精度較高，且其參數初值依賴程度較低，對於復雜的線性、非線性及多態性、多峰值問題在全局優化方面有著其他方法所無法比擬的優勢，具有很高的推廣和應用價值。
4.2 不足之處
遺傳算法雖然可以在多種領域都有實際應用，並且也展示了它潛力和寬廣前景；但是，遺傳算法還有大量的問題需要研究，目前也還存在著各種不足。首先，在變量多，取值範圍大或無給定範圍時，收斂速度下降；其次，可找到最優解附近，但無法精確確定最優解位置；最後，遺傳算法的參數選擇尚未有定量方法。對於遺傳算法，壹是還需要進壹步研究其數學基礎理論；二是還需要在理論上證明它與其他優化技術的優劣及原因；三是還需研究硬件化的遺傳算法；以及遺傳算法的通用編程和形式等。此外，對於地下水滲流問題的數值解反求多類各種水文地質參數雖有成功實例，對於運算速度問題，還存在著相當大的難度。
參考文獻
陳崇希，唐仲華.1990.地下水流動問題數值方法.武漢：中國地質大學出版社
陳喜.1998.含水層水文地質參數自動優選方法.工程勘察，（2）
郭東屏.1994.地下水動力學.西安：陜西科學技術出版社
GB 50027—2001 供水水文地質勘察規範
李俊亭，王愈吉.1987.地下水動力學.北京：地質出版社
劉寶碇，趙瑞清，王綱.2003.不確定規劃及應用.北京：清華大學出版社
朱國祥，王峰.1999.利用配線法水文地質參數計算機程序簡介.工程勘察，（3）
鄒正盛，趙智榮.2001.淺析抽水水文地質參數確定中的問題.水文地質工程地質，（3）