興隆莊采煤塌陷區地下水補給規律分析

馮忠倫,袁娜,王雅欣,徐健,楊丹,林洪孝

興隆莊采煤塌陷區地下水補給規律分析

馮忠倫,袁娜,王雅欣,徐健,楊丹,林洪孝*

山東農業大學水利土木工程學院,山東泰安271018

本文在分析降水與地下水埋深關系的基礎上，采用地下水位波動法研究興隆莊采煤塌陷區蓄水與地下水的補給規律，得出研究區地下水補給量的變化規律為從東北向西南遞減，與地下水流向基本一致；補給量隨降水量增大而增大，但在降水量增大到一定程度后變化趨于平穩，與地下水埋深的關系則為先增后減；將塌陷前后地下水補給量進行對比，表明塌陷后地下水對塌陷區蓄水有較大補給作用。此研究成果可為塌陷區蓄水與地下水的綜合開發利用提供重要參考。

地下水補給;水位波動;塌陷區;分析

地下水補給是地下水系統的基本組成部分[1]，準確評價地下水的補給量或補給資源是分析水文循環規律、合理制定水資源規劃和地下水可持續開采方案的基礎，對于實現當地經濟社會的可持續發展具有重要的戰略意義[2]。興隆莊煤礦位于濟寧曲阜市苑莊鎮境內，由于煤礦資源多年開采，塌陷區面積達到4 km2，蓄水容量約2057.8×104m3，對生態環境造成極大破壞。目前主要的治理技術包括土壤基質改良、植被恢復、土壤動物和微生物應用等[3]，開發利用途徑有農林魚禽生態恢復模式、動態塌陷區林（果）、草、牧生態恢復與重建模式、集約化農業生態恢復模式，而對于塌陷區水資源的開發利用研究較少。本文在研究塌陷區蓄水與地下水相互補給關系的基礎上，提出向工業供水[4]，為綜合治理和利用采煤塌陷區提供重要參考依據。

1 研究區概況

本文研究的興隆莊采煤塌陷區位于濟寧曲阜市苑莊鎮境內，地處北以京滬鐵路到西郭家村一線為界，南以程家莊站西約1km處京滬鐵路到中疃村一線為界，東邊平行于京滬鐵路，北起于西郭家村，南止于中疃村，北邊以京滬鐵路西約2 km為界，東西寬約4 km，南北長約6 km，形狀為西北到東南向的狹長方型，面積約為92.233 km2；其中興隆莊采煤形成的6個較大塌陷區，塌陷面積達到4～5 km2，水深1～12 m不等，蓄水容積約為2057.8×104m3，見圖1所示。經多年塌陷，目前已成相對穩定狀態。

圖1 興隆莊采煤塌陷區位置分布圖Fig.1 Location of Xinglongzhuang coal mining subsidence area

2 研究方法

由于地下水補給的復雜性，確定地下水補給的方法有多種，其中包含水均衡法、基流分割法、流域數值模擬法、零通量法、蒸滲儀法、達西定律、氯元素守恒法、放射性同位素法、地下水位波動法[5-6]。本文采用地下水位波動法研究興隆莊采煤塌陷區的地下水補給規律，認為非承壓含水層的水位變化主要歸因于地下水的補給，其計算公式如下[7]：

式中：R—地下水補給量，mm；yS—給水度；hΔ—tΔ時間內地下水位變化幅度，m。

公式（1）假定地下水進入到含水層后，立即進入到貯存量中，并進行補給重新分配，如蒸發和側向徑流[5]，該方法可以用來估算長期的地下水補給量（季度或年度）[7]。本文中關于yS和hΔ的確定方法如下：

2.1hΔ的確定

h Δ是指地下水位最高點與前期地下水位下降曲線在該點的對應值之差，即地下水變化幅度。地下水位下降曲線是假設在沒有降雨的情況下地下水位的變化軌跡。目前確定hΔ的主要方法有圖形外推法、主回歸曲線法（Master Recession Curve）和RISE程序法[8]。本文采用圖形外推法計算地下水的變化幅度。

2.2yS的確定

Sy是地下水資源評價、動態預報、水文地質計算、農田排水及滲流計算中一個十分重要的參數[9]。其值直接影響到補給量的計算。給水度的測定計算方法主要有原狀土取樣釋水試驗測定、非穩定流抽水試驗、水量平衡法和數值解法[10]。其中以抽水實驗應用最為廣泛，本文以抽水試驗計算該值。

3 計算結果及分析

3.1基礎資料的獲取

本文選取研究區6個塌陷區附近的10處地下水位觀測站（鄒城M1140010、鄒城M1110590、兗州M1111550、兗州M1111640、曲阜M112143A、曲阜M112154A、曲阜M112200A、曲阜M112205A、曲阜M112215A、曲阜M112216A）2004～2008年的地下水觀測資料，測量頻率為5 d1次；同時選取研究區內及臨近菠蘿樹、尼山水庫、書院、西諏、兗州水文站2004～2008年的降雨觀測資料進行統計分析，其地下水位與降水量的對應動態變化情況見圖2，地下水位的變化隨降雨量的變化呈現出周期性的波動，地下水位與降雨量聯系密切。

圖2 各觀測井地下水位與降雨量關系Fig.2 Water level in wells and average precipitation for study area

3.2補給量計算

地下水的變化幅度通過圖形外推法獲取，給水度采用《山東省曲阜市苑莊水源地供水水文地質詳查》抽水試驗的結果，取040=.Sy，由式（1）計算得到年補給量，計算結果見表1。

表1 補給量計算Table 1 Calculated result of groundwater recharge

研究區年平均補給量為147.39 mm，其中補給量最大的是曲阜M112216A觀測井，補給量為260.85mm，最小的是鄒城M1140010觀測井，補給量為83.09mm，比較結果見圖3。補給量的變化規律為東北向西南方向遞減，與地下水流動方向基本一致，見圖4。

圖3 各觀測井計算結果比較Fig.3 Comparison of recharge in the observation wells

圖4 多年平均地下水位等值線Fig.4 Contour lines of average groundwater table

3.3地下水補給影響因素分析

3.3.1 補給量與降水量的關系各觀測井的地下水補給來源主要為大氣降水，通過分析各觀測井降水量與補給量的相關關系見圖5。可以看出，總體上，隨著降水量增大，補給量增大，但在降水量達到950mm時，其變化趨勢減弱，尤其是圖中圈出來的3個點，在降水量較大時其補給量較小。

圖5 降雨量與補給量相關關系Fig.5 Correlation between rainfall and recharge

3.3.2 補給量與地下水位埋深的關系各觀測井的地下水位埋深差異較大，見圖6。地下水位埋深最淺的是鄒城M1140010觀測井，只有2.4 m，埋深最深的是鄒城M1110590觀測井，為14.1 m。通過分析地下水位埋深與補給量相關關系圖見圖7，可以看出，地下水的補給量與地下水位埋深之間的關系為先上升后下降，在水位埋深4 m以上的區域，地下水容易接受補給，其補給量隨地下水位埋深的增大而增大，這是因為越靠近包氣帶上部區域越容易受外界干擾，部分降雨在未達到飽水帶時就通過土壤蒸發或植物散發消耗掉，從而導致能到達飽水帶的補給量小。但在水位埋深4 m以下的區域，其補給量隨地下水位埋深的增大而減小，這是由于在到達潛水面之前水量大部分都留在了包氣帶，包氣帶越厚，能到達潛水面的水量就越少。

圖6 觀測井地下水位埋深Fig.6 Ground water level depth of observation wells

圖7 各觀測井水位埋深與補給量關系Fig.7 Correlation between groundwater level depth and recharge

3.4塌陷前后補給量對比

選取與塌陷區水力聯系密切的5個觀測井（兗州M1111550、兗州M1111640、曲阜M112143A、曲阜M112200A、曲阜M112205A）塌陷前后的地下水觀測資料進行計算分析。由于興隆莊煤礦于1981年投產運行，因此選取1975～1979年的地下水觀測數據計算塌陷前地下水補給量，由式（1）計算得到年補給量，計算結果見表2。塌陷后的地下水觀測資料仍選取2004～2008年的數據進行計算，計算結果見表1。

表2 塌陷前補給量計算Table 2 Calculated result of groundwater recharge before the collapse

圖8 塌陷前后地下水位關系比較Fig.8 Comparison of groundwater levels before and after the collapse

圖9 塌陷前后補給量比較Fig.9 Comparison of recharge before and after the collapse

塌陷后，各觀測井的地下水位較塌陷前均有不同程度的下降，見圖8。塌陷前該區域的年平均補給量為114.5 mm，塌陷后年平均補給量為130.0 mm，部分觀測井塌陷后的補給量大于塌陷前的補給量，見圖9。這一現象主要是由于降雨量的變化引起的，塌陷前年平均降雨量為656.6 mm，塌陷后年平均降雨量為854.52 mm。若按照降雨量的變化關系對補給量進行同倍比放大，則塌陷前該區域的年平均補給量為148.8 mm，表明在相同降雨情況下，塌陷后地下水的補給量小于塌陷前的地下水補給量。而該區域淺層地下水與深層地下水聯系較弱，不存在越流補給的情況，因此，塌陷后地下水補給量較塌陷前減少只能是地下水側向補給塌陷區蓄水和塌陷區蓄水增大了蒸發量造成的。

4 結論

(1)研究區年平均補給量為147.39mm，其中補給量最大的是曲阜M112216A觀測井，補給量為260.85 mm，最小的是鄒城M1140010觀測井，補給量為83.09 mm。補給量的變化規律為從東北向西南方向遞減，與地下水流動方向基本一致。

(2)研究區補給量與降水關系密切，隨降水量增大而增大，但在降水量達到950mm之后，變化趨勢減弱。

(3)補給量隨地下水埋深的關系為先增后減，埋深在小于4 m的區域，補給量隨埋深的增大而增大，但埋深大于4 m的區域，補給量隨埋深的增大而減小。

(4)塌陷后地下水的補給量比塌陷前要小，表明塌陷區周圍地下水對塌陷區蓄水有較大的側向補給。

[1]SANFORD W.Recharge and groundwater models:An overview[J].Hydrogeology Journal,2002,10(1):110-120

[2]汪丙國.地下水補給評價方法研究[D].北京:中國地質大學,2008

[3]劉飛,陸林.采煤塌陷區的生態恢復研究進展[J].自然資源學報,2009(04):612-620

[4]林洪孝.兗礦國宏化工有限公司曲阜苑莊地下水泗河水興隆莊塌陷區蓄水聯合調度開發利用技術研究[R].山東:山東農業大學,2011

[5]尹立河.基于多種方法的地下水補給研究[D].北京:中國地質大學,2011

[6]張志才.降雨入滲補給地下水研究[D].南京:河海大學,2006

[7]Richard WH,Peter GC.Using groundwater levels to estimate recharge[J].Hydrogeology Journal,2002,10:91-109

[8]Geoffrey ND,Richard WH,David LL,et al.Comparison of local to regional-scale estimates of ground-water recharge in Minnesota,USA[J].Hydrogeology Journal,2007,334:231-249

[9]李洋,褚立孔,蒲治國.確定含水層給水度新方法[J].江蘇地質,2006(04):290-293

[10]劉學軍,揚維仁.給水度測定方法研究[J].地下水,2003(04):221-223-229

Analysis on the Ground water Supply Regular in Subsidence Area in Xinglongzhuang Coal Mining

FENG Zhong-lun,YUAN Na,WANG Ya-xin,XU Jian,YANG Dan,LIN Hong-xiao*
Water Conservancy and Civil Engineering College,Shandong Agricultural University,Taian271018,China

Based on the analysis of the relationship between the rainfall and groundwater depth on the groundwater recharge, the water-table fluctuation method was adopted in this paper to estimate the groundwater recharge around Xinglongzhuang coal mining subsidence area.This paper reached a conclusion that the changes of groundwater recharge amount decreased from northeast to southwest,and basically the same as the groundwater flow.The groundwater recharge generally increased with rainfall,however,as the rainfall increases to a certain extent,the change trend will not increase obviously,but the relationship with the buried depth of groundwater is first increased and then decreased.Meanwhile,the comparison of groundwater recharge before and after the collapse shows that the subsidence area will be recharged by groundwater after the collapse.The research results can provide an important reference for the comprehensive development and utilization of water in the subsidence area.

Groundwater recharge;water-table fluctuation;subsidence area;analysis

TV213.9

A

1000-2324(2014)04-0576-05

2013-06-04

2013-06-24

國家自然科學基金資助項目(41202174);科技部國際科技合作與交流計劃項目(2007DFB70200);高等學校博士學科點專項科研基金資助項目(20123702120020)

馮忠倫(1989-),男,碩士研究生,主要從事水資源開發利用與管理方面的研究.

*通訊作者:Author for correspondence.E-mail:linhx@sdau.edu.cn