潛水-微承壓含水層水文地質(zhì)參數(shù)求解方法研究

任柳妹,楊軍耀,呂 路

(太原理工大學(xué)水利科學(xué)與工程,山西太原030024)

潛水-微承壓含水層水文地質(zhì)參數(shù)求解方法研究

任柳妹,楊軍耀,呂 路

(太原理工大學(xué)水利科學(xué)與工程,山西太原030024)

潛水-微承壓含水層水文地質(zhì)參數(shù)是正確評價山西臨汾澇洰河生態(tài)建設(shè)工程,確定水庫滲漏量的重要依據(jù)。以澇河河谷中段C6典型抽水試驗為基礎(chǔ),基于含水層試驗(Aquifer Test)專業(yè)軟件,分析多種方法下獲取的潛水-微承壓含水層水文地質(zhì)參數(shù)。結(jié)果表明,考慮潛水重力釋水的Boulton法和水位恢復(fù)法求得的水文地質(zhì)參數(shù)穩(wěn)定可靠,可為后期水文地質(zhì)數(shù)值模擬提供基礎(chǔ)數(shù)據(jù)。

水文地質(zhì)參數(shù)；潛水-微承壓含水層；含水層試驗；Boulton法；澇洰河生態(tài)建設(shè)工程

0 引 言

水文地質(zhì)參數(shù)是水文地質(zhì)條件中反映含水層或透水層水文地質(zhì)性能的指標(biāo),是進行地下水資源評價、含水層污染風(fēng)險評價以及地下水?dāng)?shù)值模擬和溶質(zhì)運移模擬的前提條件[1],其精確度直接影響后期水文工作?？紤]人工求參的繁瑣以及人為性,利用計算機求解水文地質(zhì)參數(shù)已得到廣泛應(yīng)用和認(rèn)可[2- 3],如MATLAB、1stOpt通用優(yōu)化計算平臺、Lingo和GAMS等[4- 8]。但考慮到此類優(yōu)化和數(shù)值算法涉及計算機編程,有相當(dāng)難度。而借助成熟的計算機平臺不僅可以節(jié)省大量資金和工作量,還具有精度高、可比性好的優(yōu)點,為了解研究區(qū)含水層水力特征、地下水?dāng)?shù)值模型和水資源評價提供重要基礎(chǔ)。

含水層試驗(Aquifer Test 4.2)是加拿大斯倫貝謝水務(wù)公司發(fā)布的最新專門用于抽水試驗和微水試驗資料分析、數(shù)據(jù)處理的圖形化分析和報告軟件,包含多種分析模型[9-10],主要有Theis、Hantush、Theis with Jacob Correction、Neuman、Papadopulos&Cooper、Cooper&Jacob、Double Porosity 等模型,能確定多種類型含水層的水文地質(zhì)參數(shù),如承壓含水層、潛水含水層、越流含水層和基巖裂隙含水層等,并具有水位預(yù)測、井群干擾降深計算、含水層補給邊界條件和含水介質(zhì)性質(zhì)判斷等功能。

本文以澇洰河生態(tài)建設(shè)工程澇河段抽水試驗為例,在分析研究了該地區(qū)潛水-微承壓含水層結(jié)構(gòu)特征的基礎(chǔ)上,應(yīng)用Aquifer Test抽水試驗軟件中考慮潛水重力釋水的Boulton法、Theis配線法、考慮弱透水層彈性釋水的Hantush with Storage法和水位恢復(fù)法這4種不同方法,探索求解潛水-微承壓水水文地質(zhì)參數(shù)的適用方法。

1 研究區(qū)概況

1.1 地形地貌

澇洰河生態(tài)建設(shè)工程位于山西省臨汾市堯都區(qū)境內(nèi),西至汾河,分布在澇河及洰河兩岸500～1 000 m不等寬度范圍內(nèi),占地面積22.4 km2。主要河流為汾河及其支流澇河、洰河,為東西向汾河匯集的汾河水系。

研究區(qū)位于臨汾盆地中部、汾河?xùn)|部的沖積平原區(qū),分布有河谷及河漫灘,以及汾河一、二、三級階地等多種微地貌單元。東高西低,澇河、洰河的河槽處形成低谷；澇洰河交匯處以西的區(qū)域地勢相對平坦,澇河與洰河之間的大片區(qū)域也相對平坦；澇河與洰河交匯處以東的區(qū)域河床兩側(cè)分布陡坎,高差變化范圍廣且坡度極大,河床內(nèi)不平坦,分布中等和較高坡；位于澇河與洰河之間的狹長溝谷區(qū)兩側(cè)也分布陡坡和極陡坡。研究區(qū)地質(zhì)地貌見圖1。

圖1 研究區(qū)地質(zhì)地貌

1.2 水文地質(zhì)條件

澇洰河生態(tài)建設(shè)工程研究區(qū)共分3個水文地質(zhì)單元或分區(qū)。即澇洰河交匯處至汾河入口段澇河左右岸為Ⅰ區(qū)；澇洰河河間地塊區(qū)、澇河右岸階地區(qū)與洰河左岸階地區(qū)為Ⅱ區(qū)；澇洰河交匯處至汾河入口、澇河河谷區(qū)與洰河河谷區(qū)為Ⅲ區(qū)。其中,澇河右岸階地區(qū)域0～30 m范圍內(nèi)沒有連續(xù)穩(wěn)定的含水層,局部夾粉砂、細(xì)砂透鏡體。淺層地下水賦存于亞砂土中,粉質(zhì)粘土或亞粘土為相對隔水層,洰河左岸階地區(qū)域與澇河右岸階地區(qū)含水介質(zhì)特征基本相同,富水性差,供水意義不大。澇洰河河間地塊區(qū)域0～50 m范圍內(nèi)沒有連續(xù)穩(wěn)定的含水層,局部夾細(xì)砂、粗砂透鏡體,淺層地下水賦存于黃土中,亞砂土及局部細(xì)砂、粗砂透鏡體構(gòu)成相對潛水含水層,下部粉質(zhì)粘土或亞粘土為相對隔水層,富水性差。研究區(qū)主要為沿河流走向的連續(xù)含水層,主要分布在澇洰河交匯處至汾河入口澇河河谷區(qū)、澇洰河交匯處上游河谷區(qū)、洰河河谷區(qū),南北向無水力聯(lián)系。

上部潛水含水層為粗砂、礫石,厚1.2～11.3 m,由東向西逐漸變厚,東部一定范圍內(nèi)表層覆蓋0.5～3.5 m的粉土,使其含水層局部具有微承壓性；中間弱透水層以粉土、粉質(zhì)粘土為主,厚8.9～12.5 m,局部夾粗砂、細(xì)砂透鏡體；下部為微承壓含水層,以中粗砂為主,砂質(zhì)粒徑由東向西逐漸變粗,靠近汾河處為粗砂、礫石,含水層厚2.4～14.3 m,頂板埋深15.2～20.8 m,底板埋深21.3～24.7 m,含水層連續(xù)性好,分布穩(wěn)定。將研究目標(biāo)層概化為潛水-微承壓含水層,含水層厚度全部利用。地下水接受降水入滲補給、側(cè)向徑流補給,灌溉回歸入滲補給量較小。地下水排泄以人工開采、側(cè)向徑流排泄為主,蒸發(fā)次之,水流向由東向西。澇河河谷段水文地質(zhì)剖面見圖2。

圖2 澇河河谷段水文地質(zhì)剖面(孔深：m；高程：m)

表2 各觀測孔觀測結(jié)果

2 抽水試驗

在入汾河口布置H13抽水主孔,澇洰河近交匯段布置G12抽水主孔,澇河河谷中段和上段布置C6和H1抽水主孔,洰河河谷中段和上段布置G9-1和G1抽水主孔(見圖1)。各抽水試驗主孔信息見表1。

表1 各抽水試驗主孔信息

采用1次降深,帶1個觀測孔的抽水試驗及水位恢復(fù)試驗。以穩(wěn)定流抽水,按非穩(wěn)定流要求觀測,觀測時間為抽水開始或停泵后均按1、2、3、4、6、8、10、15、20、25、30、40、50、60、80、100、120 min各觀測1次,以后每隔30 min觀測1次,同步觀測觀測孔的水位埋深,用電測水位儀觀測水位,誤差不超過±5 mm,用水表觀測流量,誤差不超過±0.1 m3,抽水穩(wěn)定延續(xù)時間不少于8 h。各觀測孔觀測結(jié)果見表2。

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

研究區(qū)抽水試驗求取的微承壓含水層水文地質(zhì)參數(shù)有滲透系數(shù)k、給水度μ、貯水系數(shù)μ*等。水位下降將引起2部分水的排出,一部分為在上部潛水面下降部位引起的重力排水,另一部分為下部飽水部分由于水頭下降引起的彈性釋水[11]。潛水層重力釋水具有普遍的滯后作用。

以澇河河谷段的抽水試驗C6數(shù)據(jù)資料為基礎(chǔ),利用Aquifer Test軟件求解微承壓含水層的水文地質(zhì)參數(shù)。抽水井C6穩(wěn)定抽水量為0.91 L/s,穩(wěn)定水位降深為5.497 m,上部為厚3.1 m的卵石潛水層,下部為厚1.4 m的細(xì)砂承壓含水層,這樣的一個雙層含水層結(jié)構(gòu)厚度記為4.5 m。觀測井C6-1距離其18.73 m。

3.1 Aquifer Test多方法求解參數(shù)

(1)考慮潛水重力釋水的Boulton法。將主孔C6和觀測井C6-1的抽水及恢復(fù)全過程試驗數(shù)據(jù)分別輸入相應(yīng)的位置,用Analysis選取Boulton分析。選擇軟件自動配準(zhǔn)后,運用參數(shù)控制功能可達到更好的擬合效果。Boulton法擬合曲線見圖3。導(dǎo)水系數(shù)T=26 m2/d,給水度μ=5.36×10-4,貯水系數(shù)μ*=4.96×10-5,μ/μ*=17.8。在澇河河谷中段,由于承壓含水層水分的運動滯后于地下水水位的降落速度,潛水給水能力略大于微承壓含水層,說明河流沖積平原地區(qū)潛水重力釋水不可忽視[12-13]。

圖3 Boulton法擬合曲線

(2)Theis配線法。選擇Theis配線法分析,自動擬合后,基于專業(yè)知識判斷,利用參數(shù)控制功能得到理想擬合曲線和參數(shù)。Theis配線法擬合曲線見圖4。導(dǎo)水系數(shù)T=32 m2/d,貯水系數(shù)μ*=6.89×10-5。

圖4 Theis配線法擬合曲線

(3)考慮弱透水層彈性釋水的Hantush with Storage法。Hantush with Storage法擬合曲線見圖5。導(dǎo)水系數(shù)T=26 m2/d, 貯水系數(shù)μ*=4.8×10-5。弱透水層彈性釋水系數(shù)比含水層小幾個數(shù)量級,且越流作用極為遲緩。本研究區(qū)潛水-微承壓含水層結(jié)構(gòu)中弱透水層厚度不是特別大,含水層抽水時弱透水層補給量可忽略不計。

圖5 Hantush with Storage法擬合曲線

圖6 水位恢復(fù)法擬合曲線

(4)水位恢復(fù)法。水位恢復(fù)法擬合曲線見圖6。抽水結(jié)束后,含水層中的地下水位開始逐漸回升,水位的變化不再受外界因素的干擾,恢復(fù)變化比較平穩(wěn),也就比較容易準(zhǔn)確觀測。因此,通過分析計算觀測井水位恢復(fù)階段的試驗數(shù)據(jù),能夠得到較為準(zhǔn)確可靠的水文地質(zhì)參數(shù)。導(dǎo)水系數(shù)T=30.8 m2/d,貯水系數(shù)μ*=5.24×10-5。

3.2 不同方法對比分析

考慮潛水重力釋水、弱透水層越流補給等因素,分別運用以上4種不同方法,求取澇河河谷中段潛水-微承壓含水層水文地質(zhì)參數(shù),分析含水層水力特性。計算結(jié)果對比見表3。從表3可知,Theis配線法忽略了潛水含水層重力釋水滯后效應(yīng),計算結(jié)果偏大；Hantush with Storage法由于考慮弱透水層越流補給,所求滲透系數(shù)較小。

表3 不同方法水文地質(zhì)參數(shù)對比

應(yīng)用GMS模擬澇洰河地下水流場,選取本次抽(注)水試驗,依次將以上4種方法所求的水文地質(zhì)參數(shù)作為初始值,計算結(jié)果與實測水位的多次擬合結(jié)果顯示,Boulton法和水位恢復(fù)法所求參數(shù)的計算水位與實測水位擬合效果較Theis配線法和Hantush with Storage法要好,絕對誤差小于0.1 m的占80%,說明計算得到的滲流場與實際觀測資料基本吻合。地下水滲流場等水位線對比見圖7。

圖7 地下水滲流場等水位線對比(單位：m)

圖7 地下水滲流場等水位線對比(單位：m)

4 結(jié) 語

本文通過分析臨汾澇洰河研究區(qū)的水文地質(zhì)條件,基于Aquifer Test專業(yè)軟件采用不同方法獲取潛水-微承壓含水層相關(guān)水文地質(zhì)參數(shù),并在此基礎(chǔ)上應(yīng)用GMS模擬地下水位評價比較各方法,得出以下結(jié)論：

(1)利用Aquifer Test計算機平臺,分別用考慮潛水重力釋水的Boulton法和Theis配線法、考慮弱透水層彈性釋水的Hantush with Storage法和水位恢復(fù)法4種求參分析方法獲取潛水-微承壓含水層水文地質(zhì)參數(shù)：導(dǎo)水系數(shù)T=28.7 m2/d、滲透系數(shù)k=6.38 m/d、貯水系數(shù)μ*=5.47×10-5、給水度μ=5.36×10-4。

(2)考慮潛水重力釋水的Boulton法和水位恢復(fù)法求得的水文地質(zhì)參數(shù)穩(wěn)定可靠,Theis配線法、考慮弱透水層彈性釋水的Hantush with Storage法求取參數(shù)值偏大。后續(xù)建立水文地質(zhì)模型,在4種方法計算結(jié)果范圍內(nèi)調(diào)參,可大大節(jié)省時間和工作量。

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

[2]SAMUEL M P, JHA M K, et al.Estimation of Aquifer Parameters from Pumping Test Data by Genetic Algorithm Optimization Technique[J]. Journal of Irrigation and Drainage Engineering, 2003, 129(5): 348- 359．

[3]江思珉, 朱國榮, 施小清, 等. 水文地質(zhì)參數(shù)反演的Hooke-Jeeves粒子群混合算法[J]. 水科學(xué)進展, 2010, 9(5)： 606- 612．

[4]師永麗, 張永波, 張志祥. 基于Aquifer Test與matlab耦合的水文地質(zhì)參數(shù)求解[J]. 水電能源科學(xué), 2015, 33(8)： 63- 66．

[5]程先云, 鄭凡東, 楊浩, 等. Theis井函數(shù)計算方法及井模型參數(shù)優(yōu)化計算研究[J]. 水文, 2015, 35(3)： 8- 13．

[6]曹先仙, 趙峰科. 基于GMS 的水文地質(zhì)參數(shù)優(yōu)化[J]. 勘察科學(xué)技術(shù), 2015(5)： 26- 29．

[7]吳春勇, 趙軍庭, 陳社明, 等. 水源地裂隙含水層水文地質(zhì)參數(shù)的確定[J]. 節(jié)水灌溉, 2013(5)： 43- 46．

[8]SU Chen, CHEN Zongyu, CHEN Jiang, et al. Mechanics of Aquitard Drainage by Aquifer-System Compaction and Its Implications for Water-Management in the North China Plain[J]. Journal of Earth Science, 2014, 25(3)： 598- 604．[9]任改娟, 楊立順, 回廣榮, 等. 傳統(tǒng)公式法和Aquifer Test 計算水文地質(zhì)參數(shù)的對比分析[J]. 地下水, 2015, 37(4)： 165- 167．

[10]薛曉飛, 張永波. 截滲墻對地下水影響的數(shù)學(xué)模型與評價[J]. 水力發(fā)電, 2015, 41(1)： 15- 23．

[11]薛禹群. 地下水動力學(xué)[M]. 北京： 地質(zhì)出版社, 1997．

[12]陳慶秋. 考慮潛水層重力釋水滯后作用條件下潛水—微承壓含水層完整單井非穩(wěn)定流的解析解[J]. 科學(xué)通報, 1993, 38(10)： 930- 932．

[13]APPIAH-ADJEI E K, SHU Longcang, Adjei K A, et al. Evaluation of unconfined aquifer parameters from flow to partially penetrating wells in Tailan River basin, China[J]. Environmental Earth Sciences, 2013, 69(3)： 799- 809．

(責(zé)任編輯 楊 健)

Research on Hydrogeological Parameter Solving Method of Phreatic-feeble confined Aquifer

REN Liumei, YANG Junyao, Lü Lu

(College of Water Resources Science and Engineering, Taiyuan University of Technology, Taiyuan 030024, Shanxi, China)

The hydrogeological parameters of phreatic-feeble confined aquifer are important indexes to accurately evaluate Lao-Ju River Ecological Construction Project in Linfen, Shanxi, and the key basis to confirm reservoir leakage. On the basis of C6 hole pumping test in the middle of Laohe River, the hydrogeological parameters of phreatic-feeble confined aquifer obtained by various methods are analyzed by Aquifer Test software. The result shows that the parameters calculated by the Boulton method which considering phreatic and reliable aquifer gravity water release and the water level recovery method are stable and reliable, which can provide basic data for later hydrogeological numerical simulation．

hydrogeological parameter; phreatic-feeble confined aquifer; Aquifer Test; Boulton method; Lao-Ju River Ecological Construction Project

2015- 10- 14

國家自然科學(xué)基金項目(41502243)；山西省自然科學(xué)基金資助項目(2015021169)

任柳妹(1991—),女,安徽合肥人,碩士研究生,研究方向為水文與水資源；楊軍耀(通訊作者)．

P641.2(225)

A

0559- 9342(2017)03- 0038- 06