東莊水庫巖溶滲漏三維數值模擬

曾 峰，蔡金龍，萬偉鋒，苗 旺

(黃河勘測規劃設計有限公司，河南鄭州 450003)

東莊水庫巖溶滲漏三維數值模擬

曾 峰，蔡金龍，萬偉鋒，苗 旺

(黃河勘測規劃設計有限公司，河南鄭州 450003)

東莊水庫2.7 km長碳酸鹽巖庫段的巖溶滲漏問題是制約工程建設的關鍵技術問題。為了分析蓄水后庫水的滲漏情況，并掌握不同防滲處理措施的效果，在總結分析東莊水庫水文地質特征的基礎上，通過建立地下水三維數值模型對水庫滲漏量進行模擬，得出不同防滲處理措施下的水庫滲漏量均屬于輕微滲漏，為防滲帷幕的設計工作提供必要的參考。

東莊水庫;滲漏量;防滲處理;三維數值模型

陜西省涇河東莊水庫的建設歷經60余年的風雨，前期勘察工作幾經中斷，其中各方研究和爭論的焦點就是2.7 km碳酸鹽巖庫段的巖溶滲漏問題，該問題是決定東莊水庫能否成壩的關鍵工程地質問題。經過勘察發現，碳酸鹽巖庫段地表存在溶洞、溶隙、溶孔等溶蝕特征，地下一定深度范圍內已存在巖溶現象，且地下水位低于河水位30～45 m，形成懸托河，在水庫建成蓄水后，存在巖溶滲漏問題。但關于水庫滲漏的部位、滲漏方式、滲漏途徑以及滲漏量卻一直沒有定論。本文就是在分析和總結項目建議書階段的成果和認識的基礎上，對巖溶滲漏問題進行進一步研究，探討影響庫水滲漏的邊界條件、左岸分水嶺等問題，對水庫可能發生滲漏部位、方式、途徑進行分析，并通過建立地下水三維滲漏模型，分析估算不同防滲處理措施的水庫滲漏量，提出可能的防滲處理措施。

1 區域水文地質概況

東莊水庫地處渭北中部，地勢整體上北高南低，自北向南呈階梯狀排列，地形差異巨大，最高海拔1599 m(鉆天嶺)，最低海拔約430 m(涇河峽谷出口)。構造條件復雜，地層分布受到構造的控制作用明顯。區內主要河流有4條，自西向東分別是:漆水河、泔河、涇河以及冶峪河。

以乾縣—富平斷裂為界，渭北巖溶可分為兩個大區:一是沿東西向展布的碳酸鹽巖斷續出露的北部山區;二是山前深埋巖溶區。其中在北部山區，受構造和當地河流排泄基準的控制，自西向東分布著5個泉域，分別是周公廟泉域(主泉高程850 m)、龍巖寺泉域(主泉高程547 m，已干涸)、篩珠洞泉域(主泉實測高程452 m)、袁家坡泉域(主泉高程377 m)和瀵泉泉域(主泉高程350 m)。東莊水庫即位于篩珠洞泉域巖溶地下水子系統內［1］。

2 篩珠洞泉域子系統邊界條件

篩珠洞泉域子系統呈西窄東寬的楔形展布，面積381.85 km2，北部以老龍山斷層北側的砂頁巖為隔水邊界，西南部以唐王陵向斜核部砂泥質碎屑巖為相對隔水邊界，東南邊界以北東走向的張家山斷裂為潛流滲水排泄邊界，地下水在此處跨越斷層帶向山前深部巖溶水系統潛流排泄(圖1)。系統內地下水總體上由西向東、東南徑流［2］。

3 碳酸鹽巖庫段巖溶水文地質特征

碳酸鹽巖庫段為老龍山斷層—東莊壩址2.7 km長河段，屬于涇河深切峽谷，岸坡陡峭，河谷斷面呈“V”字型，出露的地層為奧陶系下統冶里—亮甲山組(O1y-l)、中統馬家溝群(O2m)及第四系松散堆積層(Q)。根據地表巖溶發育特征、地質構造、地層巖性、薄片鑒定和實測地質剖面資料，并結合礦物成分分析成果，將2.7 km碳酸鹽巖庫段巖溶發育程度分為A(中等發育區)、B(弱發育區)、C(弱發育區)三個區，見圖2。

圖1 篩珠洞泉域子系統邊界條件示意圖Fig.1 The sketch showing the boundary conditions of subsystem in Shaizhudong spring area

圖2 碳酸鹽巖庫段巖溶發育程度分區圖Fig.2 The zoning map for karst development degree in reservoir segments of carbonate rock

庫壩區地下水屬于碳酸鹽巖巖溶裂隙水，無可靠的不透水層或相對隔水層。庫壩區地下水流向總體與河流走向一致，向南、南東方向徑流。據鉆孔揭露，近河兩岸地下水位低于河水位，涇河在此段為懸托型河谷。在前期各階段的工作中，曾推測在庫區左岸鉆天嶺一帶有地下水的分水嶺，但一直沒有充足的證據。項建補充階段在庫區左岸的鉆天嶺附近布置的兩個新鉆孔(YRZK01和YRZK05)地下水位分別為737 m和723 m(2013年7月)，初步判定在鉆天嶺附近存在地下水分水嶺［3］(圖3)。

圖3 篩珠洞泉域地下水等水位線圖Fig.3 The sketch showing the groundwater contour of Shaizhudong spring area

4 水庫滲漏途徑和滲漏形式分析

4.1 滲漏途徑

右岸遠處五峰山一帶存在高于水庫正常高水位的地下水位，地下水整體自西向東徑流，水庫蓄水后，不會向西產生永久滲漏。其中右岸A、B區產生的滲漏量有限;C區為庫首區厚層、巨厚層灰巖，蓄水后由于水力梯度相比蓄水前陡增，存在滲漏問題，是防滲的重點。

左岸鉆天嶺附近存在地下水分水嶺，但分水嶺水位低于正常蓄水位且分水嶺的分布范圍尚不明確，因此左岸仍存在向東滲漏的可能，根據現有資料分析，蓄水后向東滲漏的可能性較小。

根據庫壩區地下水位觀測資料，在巨厚層灰巖段(C區)，尤其是靠近河床部位，地下水位與河水有一定的水力聯系，因而，存在壩基和兩岸壩肩繞壩滲漏問題，也是繞壩滲漏的主要部位。

東莊壩址左岸約7.0 km河段為一河灣地段，弦長約3.5 km。河灣地塊地層巖性為奧陶系中統(O2)灰巖，無穩定連續的隔水層，因此，水庫蓄水后在庫首段可能存在河灣滲漏問題。

4.2 滲漏形式

根據地表地質測繪、鉆孔、平硐揭露的巖溶現象，并結合壓水試驗、示蹤試驗等成果分析，水庫巖溶滲漏形式以溶隙型為主，局部淺表層風化卸荷帶、斷層破碎帶可能存在順縫隙的脈管型滲漏。

5 水庫滲漏量三維數值模擬

5.1 水文地質概念模型

本次數值模擬的計算范圍為整個篩珠洞泉域子系統，呈狹長三角形狀，位于老龍山斷層以南，唐王陵向斜的北翼，模型總面積約380 km2，總厚度約150～1250 m，頂部為現狀地表，主要接受大氣降雨的面狀補給以及河流的線狀補給，底部以透水性較弱的奧陶系碳酸鹽巖頂面為界，可視為隔水邊界;北部以老龍山斷層為界，老龍山斷層南側為奧陶系的碳酸鹽巖，北側為二疊系和三疊系的碎屑巖，由于構造作用，巖層近直立，構成隔水邊界;西南部為唐王陵向斜核部碎屑巖組成的隔水邊界;根據地下水水化學與同位素資料推測有來自龍巖寺泉域地下水的徑流補給，巖溶區西南邊界為流量邊界;東南邊界為北東走向的張家山斷裂，在張家山斷層南側碳酸鹽巖陷落并與松散層對接，在宋家山以東構成流量邊界;篩珠洞泉和風箱道泉作為溢出邊界處理。

5.2 地下水滲流數學模型

本次研究的重點是模擬水庫蓄水后碳酸鹽巖庫段的滲漏特征，故采用穩定流進行模擬。根據上述水文地質概念模型，可以建立下列與之相適應的地下水運動三維穩定流數學模型:

式中:H，Hr為地下水位標高(m)，涇河水位標高(m); Kxx、Kyy、Kzz為含水層各向異性主方向滲透系數(m/d); Kr為河床淤積層垂向滲透系數(m/d);μs為潛水含水層給水度;W為源匯項(m3/d);h為點(x，y，z)在t時刻的水頭值;h0為初始水位標高(m);q，Qr為變流量邊界流量(m3/d·m2)，河水入滲量或地下水溢出量(m3/d); A，Mr為涇河入滲與地下水溢出面積(m2)，河床淤積層厚度(m);ε為潛水面垂向交換量(入為正、出為負) (m3/d·m2);x，y，z，t為坐標變量(m)，時間變量(d); Γ1，Γ2為一類邊界，二類邊界;n，Ω為二類邊界外法線方向，計算區范圍。

5.3 地下水滲流數值模型的建立

5.3.1 網格剖分

本次三維數值模擬計算采用的是Visual Modflow 4.1軟件系統，MODEFLOW計算采用三維有限差分方法，六面體單元。模型中的網格間距為200 m，在碳酸鹽巖庫段加密為50 m間距，垂直方向上分四個模擬層，模型頂板高程按1∶5萬數字地形圖確定，模型底板高程為280 m。模擬層的厚度依地形高度變化。

5.3.2 參數識別

本次模擬主要在以往采集的水位資料的基礎上，增加本次新測量地下水位點，用模擬區控制點的實測水位與計算水位吻合程度來對模型進行校正，通過不斷調整水文地質參數及源匯項，使得區內控制點的實測和計算水位達到擬合要求，從而完成模型的校正。圖4即為校正后模型計算出的地下水天然流場。

圖4 模型計算的地下水天然流場及觀測點分布圖Fig.4 Groundwater natural flow field calculated by model and the distribution of observation points

各觀測孔實測水位與計算水位擬合情況見圖5。

圖5 各觀測孔水位擬合結果圖Fig.5 The fitting results for the water level in observation holes

根據模型識別結果，碳酸鹽巖裸露區降水入滲系數取0.2，碳酸鹽巖裸露河谷區降水入滲系數取0.15，黃土覆蓋區降水入滲系數取0.016(圖6)。工程區多年平均降水量為560 mm(降水強度為1.5×10-3m/d)，降水入滲補給強度分別為3.2×10-4m/d、2.3×10-4m/d和2.5×10-5m/d，第三層參數分區見圖7。

圖6 研究區降雨入滲系數分區Fig.6 Zoning for infiltration coefficients of rainfall in study area

圖7 第三層滲透系數K分區(m/d)Fig.7 The zoning of infiltration coefficient K in third layer(m/d)

5.4 滲漏量數值模擬預測

5.4.1 模擬工況

為了分析不同帷幕條件下水庫滲漏量的變化情況，本次模擬設定以下幾種不同工況進行水庫滲漏量的預測:

工況1 壩址區和左、右岸均設置防滲帷幕，即全封閉方案。

工況2 壩址區設置防滲帷幕，左岸A、B和C區均設置防滲帷幕。

工況3 僅在壩址區設置防滲帷幕，并向壩肩兩岸延伸接弱透水(3 Lu)巖體。

工況4 僅在壩址區設置防滲帷幕。

5.4.2 模擬計算結果

僅在壩基設置防滲帷幕，計算蓄水后滲漏量與地下水流場。當壩基帷幕底界高程為380 m時，蓄水后等水頭線分布見圖8。

圖8 工況0蓄水后地下水位等水位線分布Fig.8 The distribution of the groundwater contour after anti-seepageand impoundment in the dam foundation

工況1 當壩基帷幕底界高程為380 m時，蓄水后地下水位等水位線分布見圖9。

圖9 工況1蓄水后地下水位等水位線分布Fig.9 The distribution of the groundwater contour after impoundment in first condition

工況2 當壩基帷幕底界高程為380 m時，蓄水后地下水位等水位線分布如圖10。

圖10 工況2蓄水后地下水位等水位線分布Fig.10 The distribution of the groundwater contour after impoundment in second condition

工況3 當壩基帷幕底界高程為380 m時，蓄水后地下水位等水位線分布見圖11。

工況4 當壩基帷幕底界高程為380 m時，蓄水后地下水位等水位線分布見圖12。

從表1中可以看出，僅壩基采取防滲措施，水庫年滲漏量達到了0.65×108m3/a，為涇河多年平均入庫徑流量(12.1×108m3/a)的5.40%，滲漏量較大，須采取工程防滲處理措施。在庫壩區分別采用不同防滲措施的其它4種工況條件下(幕底高程380 m，不考慮淤積)，水庫年滲漏量分別為1923×104m3/a、2018× 104m3/a、1974×104m3/a和3339×104m3/a，占多年平均徑流量的1.59%～2.76%，庫壩區采取防滲處理后的水庫滲漏量減小明顯，均屬于輕微滲漏。

圖11 工況3蓄水后地下水位等水位線分布Fig.11 The distribution of the groundwater contour after impoundment in third condition

圖12 工況4蓄水后地下水位等水位線分布Fig.12 The distribution of the groundwater contour after impoundment in fourth condition

表1 各工況滲漏量對比分析一覽表Table 1 The list for comparative analysis on leakage in different conditions

6 結論

(1)渭北地區巖溶地下水分為北部山區巖溶水和南部山前深埋巖溶水，東莊水庫處于北部山區篩珠洞泉域巖溶地下水子系統(I1)中沙坡斷層以北水文地質單元(I1

1)內。(2)東莊水庫庫壩區左岸的鉆天嶺附近鉆孔水位720～740m，存在地下水分水嶺。

(3)水庫主要的滲漏途徑為近壩段繞壩滲漏、庫首河灣滲漏以及A、B區通過C區向下游的繞滲。滲漏形式以溶隙型為主，局部淺表層風化卸荷帶、構造破碎帶可能存在順縫隙的脈管型滲漏。

(4)在庫壩區分別采用不同防滲措施的4種工況條件下，水庫年滲漏量分別為1923×104m3/a、2018 ×104m3/a、1974×104m3/a和3339×104m3/a，占多年平均徑流量的1.59%～2.76%，庫壩區采取防滲處理后的水庫滲漏量減小明顯，均屬于輕微滲漏。

［1］ 萬偉鋒，鄒劍峰，曾峰，等.陜西省涇河東莊水利樞紐工程可行性研究階段巖溶滲漏專題研究報告［R］.鄭州:黃河勘測規劃設計有限公司，2014.

［2］ 周益民，王泉偉，萬偉鋒，等.陜西省涇河東莊水利樞紐工程項目建議書補充報告［R］.鄭州:黃河勘測規劃設計有限公司，2014.

［3］ 萬偉鋒，王泉偉，鄒劍峰，等.東莊水庫巖溶滲漏幾個關鍵問題的探討［J］.人民黃河，2015，37(2):99-103.

(責任編輯:于繼紅)

Three-dimensional Simulation of Karst Leakage of Dongzhuang Reservoir

ZENG Feng，CAI Jinlong，WAN Weifeng，MIAO Wan
(Yellow River Engineering Consulting Co.，Ltd.，Zhengzhou，Henan 450003)

The karst leakage of 2.7km carbonate rock is one of the key technological problems that restrict the construction of Dongzhuang Reservoir.In order to analyze the leakage of reservoir and the effect of different anti-seepage treatment，on the basis of summary and analysis of hydrological geology characteristics，a three-dimensional groundwater numerical model of leakage was established，it confirm that reservoir leakage of different anti-seepage treatment are slight，which provide the necessary reference for the design of impervious curtains.

Dongzhuang reservoir;leakage;anti-seepage treatment;three-dimensional simulation

TV697.3+2

:A

:1671-1211(2015)05-0538-05

10.16536/j.cnki.issn.1671-1211.201505004

2015-06-23;改回日期:2015-07-28

曾峰 (1982-)，男，高級工程師，博士，地質工程專業，從事水文地質及環境地質方面的研究。E-mail:zengfeng1982 @126.com

數字出版網址:http://www.cnki.net/kcms/detail/42.1736.X.20150818.0841.002.html數字出版日期:2015-08-18 08:41