李　舒,陳元芳,李致家

(1. 河海大學水文水資源學院,江蘇 南京　210098; 2. 黃河水利科學研究院,河南 鄭州　450003)

?

井工礦開采對窟野河水資源的影響

李舒1,2,陳元芳1,李致家1

(1. 河海大學水文水資源學院,江蘇 南京210098; 2. 黃河水利科學研究院,河南 鄭州450003)

為了定量研究流域內井工礦開采對水資源的影響,采用Mann-Kendall-Pettitt方法檢測窟野河1966—2009年年徑流的突變點,利用水文模型(SWAT和SWAT-VISUAL MODFLOW)定量計算井工礦開采對窟野河流域1997—2009年徑流及2009年地下水的影響。結果表明,井工礦開采是該時期徑流銳減及地下水疏干的主要原因,其對年徑流減少的貢獻量為24.20 mm,占總減水量的64.16%,對地下水的影響量為4.09×108m3。

井工礦開采;徑流;地下水;Mann-Kendall-Pettitt方法;SWAT模型;SWAT-VISUAL MODFLOW耦合模型;窟野河流域

窟野河流域內蘊藏著豐富的優(yōu)質煤炭資源。1998年原煤產量約為0.2174×108t,2009年激增到2.03×108t。井工礦開采引發(fā)河川基流減少,地下水位下降[1,2]。隨著遙感和地理信息技術的發(fā)展,利用水文模型可以量化研究大多數(shù)人類活動對徑流及地下水的影響[1,3]。但國內外利用水文模型定量評估窟野河流域內井工礦開采對徑流及地下水的影響研究很少,本文針對這個問題研究了1979—2009年井工礦開采對窟野河徑流及地下水的影響。

1　研究區(qū)域與數(shù)據(jù)

圖1　窟野河水文站、雨量站及氣象站分布Fig. 1　Distributions of hydrological, rainfall, and meteorological stations in Kuye River Basin

窟野河全長242 km,流域面積約8 706 km2,位于東經(jīng)109°28′～110°45′、北緯38°22′～39°50′。本研究使用的數(shù)據(jù)有:2000年分辨率30 m×30 m的ASTER GDEM,1966—2009年的日序列水文數(shù)據(jù)(圖1),1∶100 000的1996年、2000年土地利用數(shù)據(jù)集,1∶1 000 000土壤類型數(shù)據(jù)集,全國地質資料館地質資料。

2　方　　法

2.1Mann-Kendall非參數(shù)秩次相關檢驗法與Pettitt變異點檢驗法相結合

Rouge等[4]將Mann-Kendall趨勢檢驗法[5-6]與Pettitt變異點檢測方法[7]相結合,提出一種能夠檢測水文序列漸變點與突變點的方法(以下簡稱MK-P方法)，并檢測了美國1 217個水文氣象站1910—2009年年平均降雨和溫度的漸變點和突變點。

2.2人類活動對徑流直接影響和間接影響的分離方法

人類活動對徑流產生的直接影響(井工礦開采)和間接影響是共存的,為了便于分離,假設它們是相互獨立的。參照水文模型模擬分離法,推導出人類活動對徑流的直接影響分離公式:

(1)

式中:ΔWT——徑流變化總量;ΔWC+HID——氣候變化和人類活動(間接影響)影響量之和;ΔWHD——人類活動直接影響量;WC+HD+HID——氣候變化和人類活動影響時期的觀測徑流量;WC+HID——氣候變化和人類活動(間接影響)影響時期的觀測徑流量;WCN+HIDN——氣候變化和人類活動影響時期的徑流量還原為氣候變化和人類活動(間接影響)影響時期的徑流量。

2.3水文模型

2.3.1SWAT模型參數(shù)率定及驗證方法

模型參數(shù)率定采用SUFI-2 算法[8],并用P和R進行不確定性分析(P、R為判斷因子)。當P≥70%,同時R<1時,模擬結果不確定性較小[9]。Nash-Suttcliffe Efficiency系數(shù)(NSE)、確定性系數(shù)R2和相對誤差RE這3個指標可以衡量模擬效果,若NSE>0.5,R2>0.5,RE在±25%范圍內,則模型模擬效果較好[10]。

圖2　窟野河地下水含水層劃分示意圖Fig. 2　Sketch of groundwater aquifer division in Kuye River Basin

2.3.2SWAT-VISUAL MODFLOW耦合模型參數(shù)率定及驗證方法

根據(jù)SWAT模型中的HRU和VISUAL MODFLOW中的有限差分網(wǎng)格特點構建HRU-cell的交互界面[11]。利用PEST算法[12]率定參數(shù),用RMSE和NRMSE指標檢驗SWAT-VISUAL MODFLOW耦合模型(以下簡稱S-VM模型)的模擬效果。NRMSE<10%表示模擬效果相當好,20%≥NRMSE≥10%表示模擬效果較好,30%≥NRMSE≥20%表示模擬效果一般,NRMSE>30%表示模擬效果差[13]。

2.4VM模型的設置

模型含水層劃分及邊界條件見圖2、圖3。窟野河的溫家川控制水文站監(jiān)測的日水位作為定水頭邊界值。應力期、時間步長、時間單位分別是月、月、日。初始導水率采用1 km×1 km網(wǎng)格大小的Kriging插值法[14]計算,水頭統(tǒng)一假設為30 m。其中,薩拉烏蘇組含水層、燒變巖含水層、侏羅系含水層和白堊系含水層的初始底部深度分別為30 m、150 m、200 m和750～1 580 m。

圖3　窟野河地下水邊界Fig. 3　Groundwater boundary in Kuye River Basin

孔號見煤深度冒落帶厚度裂隙帶厚度冒裂帶厚度83259.6917.7399.98117.71S432.1614.5469.0683.60S3163.2011.5553.8265.37

2.5井工礦情景設置

本次研究引入“三帶”理論[15]作為井工礦開采情景設置準則。根據(jù)鉆孔位置(圖3)假設了3條裂縫(長5 km、寬1 km)。針對窟野河煤層上覆巖層特性選取式(2)(3)計算井工礦開采下的定水頭邊界:

(2)

(3)

式中:H1——冒落帶高度;Hη——導水裂隙帶高度;M——煤層采厚或厚度。

由表1和實地考察可以判斷采空區(qū)造成的裂隙到地表,由此裂縫的定水頭設定為煤礦埋深高度。

根據(jù)MK-P突變點分析,以及1979年后大規(guī)模實施水土保持措施,1996年后煤礦大規(guī)模開采,將1979—2009年分為2個時期:第一時期(1979—1996),第二時期(1997—2009)。研究過程見圖4。

圖4　研究過程示意圖Fig. 4　Sketch of methodological steps

3　結　　果

3.1井工礦開采對徑流影響

按照第2節(jié)介紹的方法計算井工礦開采對窟野河徑流的影響,模擬結果見圖5。第一時期率定期和驗證期模型模擬效果判定值計算結果如下:NSE和R2分別達到了0.76、0.60以及0.76、0.61,RE分別為6.1%,8.0%。P和R分別達到了0.73、0.61以及0.7、0.46。圖6顯示了第一時期和第二時期的年徑流觀測值和模擬值。井工礦開采對徑流減少的貢獻量為24.20 mm,占總減水量的64.16%。

3.2井工礦開采對地下水的影響

窟野河2009年S-VM地下水模型模擬結果見圖7。A情景(井工況開采情景)下模型率定期和驗證期的RMSE和NRMSE分別為0.27、0.1及0.02%、0.009%。驗證期RMSE和NRMSE值均小于率定期,說明該耦合模型適用于模擬流域井工礦開采條件下地下水的變化情況。圖8顯示了2009年12月A情景、B情景下第一含水層水頭情況(B情景即無井工況開采情景),白色部分表示該區(qū)域的地下水已疏干。A情景的含水層疏干面積約是B情景的4倍,由此可以推斷出井工礦開采加劇了地下水的疏干。

圖5　窟野河1979—1996年月徑流觀測值與模擬值Fig. 5　Observed and simulated values of monthly runoff in Kuye River from 1979 to 1996

圖6　窟野河1979—2009年年徑流觀測值與模擬值Fig. 6　Observed and simulated values of annual runoff in Kuye River from 1979 to 2009

圖7　2009年3—10月地下水位模擬值與實測值擬合結果Fig. 7　Fitting of observed and simulated groundwater levels in March to October 2009

圖8　2009年12月地下水水位示意圖Fig. 8　Sketch of groundwater level in December 2009

通過水均衡(表2)及面積加權法計算出2009年井工礦開采對地下水的總影響量為4.09×108m3,與吳喜軍[1]利用統(tǒng)計學方法得出的結果一致。B情景下,含水層在一定范圍內出現(xiàn)了疏干情況,說明本次研究設置的S-VM模型參數(shù)也反映了人工打井取水對地下水的影響。以上結果表明本次研究假設的A、B情景具有一定的合理性。

表2　研究區(qū)地下水均衡結果Table 2　Groundwater balance in study area

注:Q地表-地下為地表水轉化為地下水量,Q地下-地表為地下水轉化為地表水量。

4　討　　論

4.1MK-P方法在窟野河流域應用要點

本次研究使用了MK-P方法檢測計算了窟野河年徑流的突變點。在使用該方法時,首先需要確定一個閾值,該閾值是區(qū)分平穩(wěn)變化點和突變點的臨界值。根據(jù)Rouge等[4]將該方法應用于美國水文站的結果,確定了閾值為0.4,計算得到的窟野河年徑流突變點與呂新等[2]的研究成果一致,并且該流域在1979年和1996年后分別進行了大規(guī)模的水保措施和煤礦開采活動,說明該閾值的取值是合理的。在不同的流域使用該方法時應注意閾值的取值范圍會在0.1～0.5之間變化,需要根據(jù)流域徑流的實際情況選取。

4.2不確定性

本次研究結果的不確定性主要有以下幾個原因:(a)SWAT模型中使用的數(shù)據(jù)集尺度為日,但窟野河流域0.5 h內的強降雨出現(xiàn)次數(shù)較多;(b)SWAT模型參數(shù)的不確定性;(c)沒有考慮流域內水工建筑物及不同的農業(yè)耕種方式給徑流帶來的影響;(d)在VM模型建立中,用Kriging插值法計算導水率時采用較低的分辨率;(e)人為假設的模型邊界與實際情況有所不同。本次研究采用P、R和NRMSE定量分析這些不確定性對計算結果的影響,其結果表明SWAT模型和S-VM模型模擬結果的不確定性都在可接受范圍內,因此本研究計算出的井工礦開采對窟野河徑流及地下水的影響量是合理的。

4.3結果的分析及展望

研究中通過一系列概化方法計算了井工礦對徑流及地下水水量的影響,利用SWAT模型計算了井工礦對徑流的影響,其結果較好。在設置井工礦開采情景時,根據(jù)實地調研情況,概化了3條地裂縫,取得較好的模擬結果。

井工礦開采會使當?shù)厮牡刭|情況發(fā)生很復雜的變化,如煤層開采后形成采空區(qū)。這些采空區(qū)在一定條件下(如強降雨等),使地表出現(xiàn)塌陷區(qū),甚至出現(xiàn)“天坑”,進而對含水層產生擾動影響,以至于改變含水層的儲水率、導水率等地質參數(shù),最終對徑流及地下水產生影響。此外,窟野河流域分布著很多礦區(qū),不同的開采方法也會對水文地質情況產生不同的影響。

5　結　　語

本次研究首先利用MK-P突變點檢測方法檢測了窟野河流域年徑流量的突變點(1979年和1996年)。采用SWAT模型和S-VM模型分別定量計算了第二時期井工礦開采對窟野河徑流的影響量以及2009年井工礦開采對窟野河地下水的影響量。計算結果表明,井工礦開采對年徑流減少的貢獻量為24.20 mm,占總減水量的64.16%,其對地下水的影響量為4.09×108m3。井工礦開采是窟野河徑流減少及松散巖類孔隙水含水層疏干的主要原因。量化流域井工礦開采對徑流及地下水的影響是一項復雜的研究,準確理解井工礦開采影響流域水循環(huán)機理是利用耦合模型模擬井工礦開采對徑流及地下水影響這一物理過程的基礎。

[1] 吳喜軍. 煤炭開采地區(qū)河道徑流變化與生態(tài)基流研究:以陜北窟野河流域為例[D].陜西:西安理工大學, 2013.

[2] 呂新, 王雙明, 楊澤元,等. 神府東勝礦區(qū)煤炭開采對水資源的影響機制:以窟野河流域為例[J]. 煤田地質與勘探, 2014, 42(2): 54-57. (LYU Xin, WANG Shuangming, YANG Zeyuan, et al. Influence of coal mining on water resources: a case study in Kuye River basin[J]. Coal Geology & Exploration, 2014, 42(2): 54-57. (in Chinese))

[3] DONG W, CUI B S, LIU Z H, et al. Relative effects of human activities and climate change on the river runoff in an arid basin in northwest China[J]. Hydrological Processes, 2014, 28(18): 4854-4864.

[4] ROUGé C, GE Y, CAI X M. Detecting gradual and abrupt changes in hydrological records[J]. Advances in Water Resources, 2013, 53(1): 33-44.

[5] MANN H. Non-parametric tests against trend[J]. Econometrica, 1945, 13: 245-259.

[6] KENDALL M, GIBBONS J D. Rank correlation methods[M]. Griffin: Oxford University Press, 1975.

[7] PETTITT A N. A non-parametric approach to the change point problem[J]. Journal of the Royal Statistical Society, 1979, 28(2): 126-135.

[8] ABBASPOUR K C, JOHNSON A C, van GENUCHTEN M T H. Estimating uncertain flow and transport parameters using a sequential uncertainty fitting procedure[J]. Vadose Zone Journal, 2004, 3: 1340-1352.

[9] ABBASPOUR K C, YANG J, MAXIMOV I, et al. Modelling hydrology and water quality in the pre-alpine/alpine Thur watershed using SWAT[J]. Journal of Hydrology, 2007, 333(2/3/4): 413-430.

[10] MORIASI D N, ARNOLD J G, van LIEW M W, et al. Model evaluation guidelines for systematic quantification of accuracy in watershed simulations[J]. Transactions of the ASABE, 2006, 50(3): 885-900.

[11] 張雪剛,毛媛媛,董家瑞,等.SWAT模型與MODFLOW模型的耦合計算及應用[J].水資源保護,2010, 26(3):49-52.(ZHANG Xuegang, MAO Yuanyuan, DONG Jiarui, et al. A coupled model simulation and application of SWAT-MODFLOW[J]. Water Resources Protection, 2010, 26(3): 49-52. (in Chinese))

[12] ADAM S, TRACY N, PETER M. Natural recharge estimation and uncertainty analysis of an adjudicated groundwater basin using a regional-scale flow and subsidence model (Antelope Valley, California, USA)[J]. Hydrogeology Journal, 2015, 23(6): 1267-1291.

[13] MICHELE R, NICOLA L, ZINA F. Evaluation and application of the OILCROP-SUN model for sunflower in southern Italy[J]. Agricultural Systems, 2003, 78(1): 17-30.

[14] NOEL C. Spatial prediction and ordinary kriging[J]. Mathematical Geology, 1998, 20(4): 405-421.

[15] 煤炭科學研究院北京開采研究所. 煤礦地表移動與覆巖破壞規(guī)律及其應用[M].北京:煤炭工業(yè)出版社,1981.

Impacts of underground mining on water resources of Kuye River

LI Shu1,2, CHEN Yuanfang1, LI Zhijia1

(1.CollegeofHydrologyandWaterResources,HohaiUniversity,Nanjing210098,China;2.YellowRiverInstituteofHydraulicResearch,Zhengzhou450003,China)

In order to quantitatively study the impacts of underground mining on the water resources of the Kuye River Basin, the Mann-Kendall-Pettitt method was used to detect the abrupt change point of the annual runoff over the period from 1966 to 2009. The impacts of underground mining on runoff of the Kuye River Basin over the period from 1997 to 2009 and groundwater in the year 2009 were quantitatively estimated using the SWAT and SWAT-VISUAL MODFLOW models. The results show that the sharp decrease in runoff and the depletion of groundwater during the study period are mainly attributed to underground mining, which is responsible for 24.20 mm of annual runoff reduction, accounting for 64.16% of the total reduction, and 4.09×108m3of the groundwater affected.

underground mining; runoff; groundwater; Mann-Kendall-Pettitt method; SWAT model; SWAT-VISUAL MODFLOW coupling model; Kuye River Basin

1000-1980(2016)04-0347-06

10.3876/j.issn.1000-1980.2016.04.011

2015-08-31

國家自然科學基金(41130639, 51179045)；水利部公益性行業(yè)科研專項(201501022, 201301068)

李舒(1986—)，男，河南鄭州人，博士研究生，主要從事水文水資源研究。E-mail: lididshu30@163.com

P339

A