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

任柳妹,楊軍耀,呂 路

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

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

任柳妹,楊軍耀,呂 路

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

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

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

0 引 言

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

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

本文以澇洰河生態(tài)建設(shè)工程澇河段抽水試驗(yàn)為例,在分析研究了該地區(qū)潛水-微承壓含水層結(jié)構(gòu)特征的基礎(chǔ)上,應(yīng)用Aquifer Test抽水試驗(yàn)軟件中考慮潛水重力釋水的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ū),分布有河谷及河漫灘,以及汾河一、二、三級(jí)階地等多種微地貌單元。東高西低,澇河、洰河的河槽處形成低谷；澇洰河交匯處以西的區(qū)域地勢(shì)相對(duì)平坦,澇河與洰河之間的大片區(qū)域也相對(duì)平坦；澇河與洰河交匯處以東的區(qū)域河床兩側(cè)分布陡坎,高差變化范圍廣且坡度極大,河床內(nèi)不平坦,分布中等和較高坡；位于澇河與洰河之間的狹長(zhǎng)溝谷區(qū)兩側(cè)也分布陡坡和極陡坡。研究區(qū)地質(zhì)地貌見(jiàn)圖1。

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

1.2 水文地質(zhì)條件

澇洰河生態(tài)建設(shè)工程研究區(qū)共分3個(gè)水文地質(zhì)單元或分區(qū)。即澇洰河交匯處至汾河入口段澇河左右岸為Ⅰ區(qū)；澇洰河河間地塊區(qū)、澇河右岸階地區(qū)與洰河左岸階地區(qū)為Ⅱ區(qū)；澇洰河交匯處至汾河入口、澇河河谷區(qū)與洰河河谷區(qū)為Ⅲ區(qū)。其中,澇河右岸階地區(qū)域0～30 m范圍內(nèi)沒(méi)有連續(xù)穩(wěn)定的含水層,局部夾粉砂、細(xì)砂透鏡體。淺層地下水賦存于亞砂土中,粉質(zhì)粘土或亞粘土為相對(duì)隔水層,洰河左岸階地區(qū)域與澇河右岸階地區(qū)含水介質(zhì)特征基本相同,富水性差,供水意義不大。澇洰河河間地塊區(qū)域0～50 m范圍內(nèi)沒(méi)有連續(xù)穩(wěn)定的含水層,局部夾細(xì)砂、粗砂透鏡體,淺層地下水賦存于黃土中,亞砂土及局部細(xì)砂、粗砂透鏡體構(gòu)成相對(duì)潛水含水層,下部粉質(zhì)粘土或亞粘土為相對(duì)隔水層,富水性差。研究區(qū)主要為沿河流走向的連續(xù)含水層,主要分布在澇洰河交匯處至汾河入口澇河河谷區(qū)、澇洰河交匯處上游河谷區(qū)、洰河河谷區(qū),南北向無(wú)水力聯(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)層概化為潛水-微承壓含水層,含水層厚度全部利用。地下水接受降水入滲補(bǔ)給、側(cè)向徑流補(bǔ)給,灌溉回歸入滲補(bǔ)給量較小。地下水排泄以人工開(kāi)采、側(cè)向徑流排泄為主,蒸發(fā)次之,水流向由東向西。澇河河谷段水文地質(zhì)剖面見(jiàn)圖2。

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

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

2 抽水試驗(yàn)

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

表1 各抽水試驗(yàn)主孔信息

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

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

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

以澇河河谷段的抽水試驗(yàn)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ì)砂承壓含水層,這樣的一個(gè)雙層含水層結(jié)構(gòu)厚度記為4.5 m。觀測(cè)井C6-1距離其18.73 m。

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

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

圖3 Boulton法擬合曲線

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

圖4 Theis配線法擬合曲線

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

圖5 Hantush with Storage法擬合曲線

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

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

3.2 不同方法對(duì)比分析

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

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

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

圖7 地下水滲流場(chǎng)等水位線對(duì)比(單位：m)

圖7 地下水滲流場(chǎng)等水位線對(duì)比(單位：m)

4 結(jié) 語(yǔ)

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

(1)利用Aquifer Test計(jì)算機(jī)平臺(tái),分別用考慮潛水重力釋水的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種方法計(jì)算結(jié)果范圍內(nèi)調(diào)參,可大大節(jié)省時(shí)間和工作量。

[1]李培月. 非穩(wěn)定流抽水試驗(yàn)確定越流承壓含水層水文地質(zhì)參數(shù)方法對(duì)比研究[D]. 西安： 長(zhǎng)安大學(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]江思珉, 朱國(guó)榮, 施小清, 等. 水文地質(zhì)參數(shù)反演的Hooke-Jeeves粒子群混合算法[J]. 水科學(xué)進(jìn)展, 2010, 9(5)： 606- 612．

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

[5]程先云, 鄭凡東, 楊浩, 等. Theis井函數(shù)計(jì)算方法及井模型參數(shù)優(yōu)化計(jì)算研究[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 計(jì)算水文地質(zhì)參數(shù)的對(duì)比分析[J]. 地下水, 2015, 37(4)： 165- 167．

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

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

[12]陳慶秋. 考慮潛水層重力釋水滯后作用條件下潛水—微承壓含水層完整單井非穩(wěn)定流的解析解[J]. 科學(xué)通報(bào), 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

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

任柳妹(1991—),女,安徽合肥人,碩士研究生,研究方向?yàn)樗呐c水資源；楊軍耀(通訊作者)．

P641.2(225)

A

0559- 9342(2017)03- 0038- 06