程林森，雷少剛，卞正富

(中國礦業(yè)大學環(huán)境與測繪學院，江蘇 徐州　221008)

?

半干旱區(qū)煤炭開采對土壤含水量的影響

程林森，雷少剛①，卞正富

(中國礦業(yè)大學環(huán)境與測繪學院，江蘇 徐州221008)

摘要：土壤含水量是影響干旱區(qū)植被生長的一個重要因素，研究由于煤礦開采引起的土壤水時序變化有利于礦區(qū)植被的重建。以大柳塔礦井52303工作面為例，采用剖面水分速測儀對開采前后壓縮區(qū)10～100 cm深處土壤含水量的變化進行監(jiān)測。通過對監(jiān)測結果的分析發(fā)現(xiàn)，未受采煤影響時，土壤含水量較平穩(wěn)，不隨溫度變化而改變；受采煤影響初期，土壤含水量變化不明顯，存在滯后性，滯后時間為地下開采到達測點正下方后4～5 d；之后，在測點達到最大沉降量過程中發(fā)現(xiàn)，下沉對土體的擾動作用導致土壤含水量出現(xiàn)先上升后下降現(xiàn)象，且地下開采對土壤水擾動程度由淺至深逐漸減弱。沉陷引起的土體擾動導致土壤粒徑減小，容重增加，孔隙比降低，使得土壤持水能力增強，這是導致在受開采干擾初期(04-03—04-10)土壤含水量短暫上升的重要原因，而裂縫產生以及雨水補給能力的降低是導致后期土壤含水量降低的主要原因。

關鍵詞：半干旱區(qū)；煤炭開采；土壤水；時序監(jiān)測

煤炭是我國重要的能源來源,在我國一次性能源消費結構中所占比例一直是70%左右[1]。目前,我國東部礦區(qū)煤炭資源趨于枯竭,國家煤炭資源開發(fā)重心正向西北地區(qū)轉移,隨著向西部進一步擴張,在干旱缺水的地區(qū)集中大規(guī)模井工開采將會導致環(huán)境問題[2],使得某些環(huán)境要素進一步惡化;煤炭開采對地表土壤水的影響是對區(qū)域生態(tài)環(huán)境影響的一個重要表現(xiàn),對其進行研究有助于生態(tài)脆弱區(qū)植被的重建。采煤對土壤水影響主要是擾動土體,改變土體理化性質[3]和土壤持水能力[4],從而引起土壤水的變化;而西部礦區(qū)井工開采具有高強度和高產特點,且埋深較淺,對地表土體擾動尤為顯著。已有研究表明,受采煤影響,塌陷區(qū)域內土壤容重、孔隙比和機械組成等物理性質發(fā)生不同程度變化[5],這也導致土壤滲透系數(shù)的改變;而地下水位降低,潛水補給能力降低,地表裂縫產生導致地表蒸發(fā)量增加,都影響土壤含水率[6];同時由于井工開采方式會導致覆巖從下至上發(fā)生冒落、裂隙和彎曲下沉,使地表產生大量裂縫[7],這也增大了土壤水蒸發(fā)面積,進一步降低土壤含水量[8]。薛豐昌等[9]通過對礦區(qū)預采區(qū)、采動區(qū)和采空區(qū)土壤含水率的調查比較,發(fā)現(xiàn)采礦影響下神東礦區(qū)域土壤含水率存在一定空間分異性,與未塌陷區(qū)相比,塌陷區(qū)土壤含水量降低;受坡位影響,采煤對土壤水影響程度也存在差異[10]。國外對采煤塌陷區(qū)土壤水的研究認為,采煤沉陷對包氣帶水分影響十分嚴重,分析了煤炭開采對土壤水力學參數(shù)及其空間變異性的影響[11-13],同時對采后塌陷區(qū)土壤結構變化進行研究[14]。從上述研究可以看出,礦區(qū)土壤水分變化研究主要集中在采前和采后變化差異上,關于土壤水分時序變化方面的研究較缺乏,而國外對于土壤水的連續(xù)監(jiān)測研究已有先例,但主要集中在流域范圍內。COSH等[15-16]通過對流域范圍內表層3～7 cm處土壤水的長時間連續(xù)監(jiān)測,采用時序穩(wěn)定分析來評估典型樣點,以獲得流域范圍內土壤水平均值;CHO等[17]通過土壤水的連續(xù)監(jiān)測來研究氣象因素對土壤水的影響,研究發(fā)現(xiàn)降雨為主要影響因素。因此,對于西部礦區(qū)土壤水的研究目前尚缺乏塌陷區(qū)不同區(qū)域(中性區(qū)、拉伸區(qū)和壓縮區(qū))以及時序監(jiān)測結果分析。筆者將對塌陷區(qū)內特定區(qū)域采前、采中和采后的土壤水進行時序監(jiān)測,以研究土壤水受采煤影響整個過程中的連續(xù)變化情況,以期為工作面上方土壤水的時序變化研究提供參考,進而有利于整個礦區(qū)的植被重建。

1研究區(qū)概況

大柳塔礦區(qū)位于陜西省神木縣西北約52.5 km處,北鄰哈拉溝,南與朱蓋塔礦井接壤,東以悖牛川為界,西以烏蘭木倫河為界。該區(qū)地處陜北高原與毛烏素沙漠東南緣的接合部位,屬半干旱大陸性季風氣候區(qū),冬季漫長寒冷,夏季短暫炎熱,風沙頻繁,暴雨集中,晝夜溫差懸殊。多年氣溫為-28.1～38.9 ℃,降水量為251.3～646.5 mm,多集中在7—9月,平均蒸發(fā)量為1 788.4 mm。所選取的52303工作面地面位置位于大柳塔礦井田的東南區(qū)域三盤區(qū),工作面長4 547.6 m,寬301 m,煤層平均厚度為6.94 m,煤層傾角為1°～3°,平均埋深為235.0 m,屬于典型的超大工作面開采。地表大部為第四系松散沉積物覆蓋,平均厚度為30.0 m,上覆基巖主要由粉砂巖和細砂巖組成,平均厚度為205.0 m,地表土壤類型主要為沙土、栗鈣土和黃土。

2研究方法

2.1監(jiān)測點的選擇

選取大柳塔礦52303工作面為研究對象,根據(jù)開采進度表在掘進方向前方還未受開采影響的區(qū)域選取監(jiān)測點,測點位于工作面上方,地形較平坦,土壤類型為硬梁地,質地均勻,采用單點連續(xù)監(jiān)測,以時間替代空間,采前土壤水的相對穩(wěn)定與采中、采后土壤水時序變化形成對照。通過對該區(qū)100 cm深度內土壤水的連續(xù)動態(tài)監(jiān)測,了解該區(qū)域在受開采影響不同階段土壤水的變化情況。

2.2數(shù)據(jù)的獲取

采用英國PR2(Profile Probe Type PR2)土壤剖面水分速測儀進行連續(xù)監(jiān)測,每30 min監(jiān)測1次,分別監(jiān)測10、20、30、40、60和100 cm深處土壤體積含水量的變化;監(jiān)測時間為2014-03-16—2014-06-03,即監(jiān)測點受地下開采影響前至沉陷相對穩(wěn)定后一段時間,獲得監(jiān)測點在受到地下開采影響前后10～100 cm深度土壤含水量數(shù)據(jù)。同時利用附近已有的2個地表沉陷觀測站,觀測站走向分別為沿工作面和與工作面方向相垂直,采用上海華測X91型接收機(標稱精度為5 mm+1 km)進行實時動態(tài)差分觀測,監(jiān)測地表下沉情況,從而分析不同下沉階段土壤水的變化。土壤取樣采用環(huán)刀法,在土壤水監(jiān)測點周圍取3個樣點,每個樣點取4個環(huán)刀,用玻璃片將環(huán)刀兩頭密封后帶回實驗室分析。

3監(jiān)測結果

3.1測點地表沉陷監(jiān)測結果

主要通過對走向觀測站的監(jiān)測結果(圖1)來判定測點位置的地表沉陷情況。

土壤水測點位置在4號站點。

圖1顯示,測點位置在03-26之前就已經(jīng)受到開采影響而開始下沉,04-05還未達到最大沉降量,之后處于下沉活躍期,04-26、05-08和06-11測點位置沉降量基本不變,說明已處于相對穩(wěn)沉狀態(tài)。

3.2土壤水監(jiān)測結果

2014-03-16—2014-06-03土壤含水量的變化趨勢見圖2。根據(jù)地表沉陷觀測結果(圖1)以及該工作面推進進度表可知,03-26之前監(jiān)測點已經(jīng)受到開采影響而開始下沉,03-29前后地下開采達到測點正下方,由圖2可知,04-03之前土壤含水量變化較平穩(wěn);之后的一段時間,受開采和降雨影響,土壤含水量變化幅度較大。因此,可分04-03之前的土壤含水量未擾動期和之后的擾動期對監(jiān)測結果進行采前和采后對比分析。

圖2　不同深度土壤含水量變化趨勢

4分析與討論

4.1未擾動期土壤水分變化分析

土壤含水量未受地下開采影響的時間為03-16—04-03,將這段時間監(jiān)測數(shù)據(jù)與氣象數(shù)據(jù)進行對比分析,研究氣象條件對土壤含水量的影響。根據(jù)氣象數(shù)據(jù),此段時間沒有降雨,所以不存在降雨對土壤水監(jiān)測數(shù)據(jù)的影響,因此分析溫度與土壤含水量變化的關系(圖3)。由圖3可知,該段時間內不同深度土壤含水量變化不隨溫度的變化而發(fā)生改變,因此溫度變化對土壤水影響不大。根據(jù)開采進度表推算地下開采于03-29到達監(jiān)測點正下方;根據(jù)觀測站觀測結果可知,該點在03-26前后就受到地下開采的影響而下沉,直到04-03不同深度土壤含水量標準差都很小(表1),均<0.01 m3·m-3,且變異系數(shù)均較小,表明地表下沉初期對土壤含水量影響不明顯,具有一定滯后性,滯后時間為地下開采到達監(jiān)測點正下方之后4～5 d。這主要是由于地表下沉初期沉陷速度相對較緩慢,對土體擾動較小,導致土壤含水量變化不明顯。

圖3　03-16—04-03期間溫度與土壤含水量變化

4.2擾動期土壤水分變化分析

土壤含水量變化在04-03之后趨于明顯(圖2),對于該段時間內不同深度土壤水變化進行相關性分析(表2)發(fā)現(xiàn),除60和100 cm深度外,各相鄰層土壤含水量變化都具有較強的正相關性,說明各層土壤含水量在受地下開采和降雨影響下的變化具有一定關聯(lián)性;而100 cm深處由于受降雨影響小,且受開采影響后土壤水較早進入穩(wěn)定期,因此和其他各層之間存在負相關性。

表1未擾動期不同深度土壤含水量變化幅度

Table 1Magnitude of the change in soil moisture content with soil depth during the undisturbed period

深度/cm標準差/(m3·m-3)變異系數(shù)/%100.0026.4200.0031.1300.0031.0400.0031.0600.0030.71000.0061.2

表2不同深度土壤含水量相關性分析

Table 2Correlation analysis of soil moisture content with soil depth

深度/cm不同深度(cm)土壤含水量相關系數(shù)20304060100100.747*0.476*0.507*0.062*-0.387*200.858*0.770*0.499*-0.403*300.722*0.744*-0.501*400.553*-0.274*60-0.506*

*表示在α=0.01水平上相關顯著(雙側檢驗)。

通過查找氣象數(shù)據(jù)(04-10、04-14、04-15、04-25、05-01、05-08和05-09降水量分別為3.8、0.1、0.2、13.8、3.4、0.1和12.0 mm),各層出現(xiàn)土壤含水量突變值是由于降水影響所致。由圖2可知,在04-09之前未受到降水影響的情況下各層都出現(xiàn)土壤含水量上升階段,這說明在受地下開采擾動地表沉陷活躍期間,各層土壤含水量并未立即下降,而是在初期先有個上升過程,在無外來降水補給情況下,其補給源只能來源于土壤含水量較高的底層,且補給量大于損失量,土壤含水量上升表明這一階段土壤持水能力上升。而從擾動程度來看,通過對各層土壤含水量變化程度(10、20、30、40、60和100 cm深處土壤含水量變異系數(shù)分別為25.3%、8.5%、7.9%、2.4%、3.0%和1.8%)進行分析可知,隨著深度的增加,土壤含水量變異系數(shù)在減小,說明開采對表層土壤水的擾動要大于深層。

土壤持水性能主要受土壤總孔隙度、毛管孔隙度、土壤容重、土壤有機質和土壤顆粒組成的影響[18]。土壤孔隙度降低會導致容重增加,從而改變土壤結構,影響土壤透氣性、入滲性能和持水能力等[19]。土壤有機質影響土壤的比表面積[20],并且自身結構疏松多孔,有利于水分保持。土壤顆粒組成在一定程度上決定了土壤的結構和性質,對土壤持水性和保肥供肥能力均有明顯影響[21]。由上述相關性分析可知,相鄰土壤層之間含水量的變化具有較強相關性,為了分析上述現(xiàn)象,且考慮到采煤對表層土壤水影響較大,因此取開采前后測點周圍10～20 cm深處土樣進行分析。通過對采樣結果的分析發(fā)現(xiàn),土壤粒徑減小(表3),主要是由于沉陷所致土體之間相互應力作用,土壤顆粒重新組合,這在一定程度上導致采后土壤孔隙比小于采前,從而使容重增加,而有機質含量變化不顯著(P>0.05)(表4)。已有研究表明,土壤孔隙大小不同,所起作用也不同,團粒內部毛管孔隙(小孔隙)能保持水分,而團粒間非毛管孔隙(大孔隙)則能保持通氣。當土壤受壓縮作用時,團粒間大孔隙容積降低,中等大小孔隙容積卻有所增加[22]。馮杰等[23]的研究表明,受擾動后原狀土的中等孔隙發(fā)育,在低吸力段持水性較好。因此,土體在受煤礦開采沉陷影響時,擠壓作用使得土壤團聚體間大孔隙減小,中等孔隙發(fā)育,土壤含水量處于上升階段。而隨著沉陷加劇,由于土體非連續(xù)移動,地表出現(xiàn)裂縫,增加了土壤水蒸發(fā)面積,從而使土壤含水量下降。

表3開采前后不同粒徑土壤占比的變化

Table 3Composition of soil particle size before mining and after mining

粒徑/mm不同粒徑土壤占比/%采前采后<0.005 3.7 8.80.005～<0.07539.436.70.075～<0.2557.154.60.25～0.50.50.3

表4開采前后土壤孔隙比、容重和有機質含量的變化

Table 4Void ratio, soil bulk density and organic matter content before and after coal mining

時間容重/(g·cm-3)天然孔隙比w(有機質)/%采前1.300.9541采后1.450.8242

從整個過程來看,土壤含水量受降水影響較大的時間為04-10、04-25和05-09,對這3個降水過程中土壤含水量變化速率(即每小時土壤含水量變化程度)進行分析。由于60 cm深處只在02-25受到降水影響,100 cm深處在整個監(jiān)測過程中受降水影響較小,因此選取10、20、30和40 cm深處進行分析,結果見表5。由表5可知,10 cm深處土壤含水量變化速率變化較小,基本平穩(wěn),主要是由于0～10 cm深處為沙土,雨水下滲速度較快所致。而其他3層在這3次降水過程中土壤含水量變化速率在減小,說明這3層土壤滲透速度降低。這主要是由于在地表下沉過程中擠壓作用使得土壤團聚體間孔隙減小,土壤導水率降低,同時降水隨裂縫流失,從而引起土壤含水量變化速率減小,這也說明雨水對土壤水的補給能力下降。

表53次降雨過程土壤含水量的變化速率

Table 5Change rates of soil moisture content during the course of three rainfall events

時間(月-日)不同深度(cm)土壤含水量變化速率/(m3·m-3·h-1)1020304004-102.018.211.226.004-251.88.29.013.205-091.63.96.26.6

5結論

未擾動期各層土壤含水量變化幅度不大,溫度對土壤水的影響較小,地下開采對測點的影響在初期不明顯,對土壤水影響存在一定滯后性,根據(jù)監(jiān)測結果分析,其影響在地下開采到達測點正下方之后4～5 d趨于明顯。

在開采擾動期間,在測點達到最大沉降量的過程中,壓縮區(qū)各層土壤水受不同程度的影響,且表層土壤水受影響程度大于深層;由于沉陷引起的土體擾動導致土壤粒徑減小,容重增加,孔隙比降低,使得土壤持水能力增強,這是導致在受開采擾動初期(04-03—04-10)土壤含水量短暫上升的重要原因,而裂縫產生以及雨水補給能力的降低是導致后期土壤含水量降低的主要原因。

參考文獻：

[1]卞正富.我國煤礦區(qū)土地復墾與生態(tài)重建研究[J].資源產業(yè),2005,7(2):18-24.

[2]錢鳴高.煤炭產能擴張引發(fā)中西部環(huán)境隱憂[J].資源環(huán)境與發(fā)展,2011(4):4-7.

[3]雷少剛,卞正富.西部干旱區(qū)煤炭開采環(huán)境影響研究進展[J].生態(tài)學報,2014,34(11):2837-2843.

[4]趙紅梅,張發(fā)旺,宋亞新,等.神府東勝礦區(qū)不同塌陷階段土壤水分變化特征[J].南水北調與水利科技,2008,6(3):92-96.

[5]王健,武飛,高永,等.風沙土機械組成、容重和孔隙度對采煤塌陷的響應[J].內蒙古農業(yè)大學學報(自然科學版),2006,27(4):37-41.

[6]卞正富,雷少剛,常魯群,等.基于遙感影像的荒漠化礦區(qū)土壤含水率的影響因素分析[J].煤炭學報,2009,39(4):520-525.

[7]張發(fā)旺,侯新偉,韓占濤,等.采煤塌陷對土壤質量的影響效應及保護技術[J].地理與地理信息科學,2003,19(3):67-70.

[8]魏江生,賀曉,胡春元,等.干旱半干旱地區(qū)采煤塌陷對沙質土壤水分特性的影響[J].干旱區(qū)資源與環(huán)境,2006,20(5):84-88.

[9]薛豐昌,卞正富.神東礦區(qū)采礦對土壤含水率影響分析[J].煤炭科學技術,2007,35(9):83-85.

[10]呂晶潔,胡春元,賀曉.采煤塌陷對固定沙丘土壤水分動態(tài)的影響研究[J].干旱區(qū)資源與環(huán)境,2005,19(7):152-156.

[11]BUCZKO U,GERKE H H.Estimating Spatial Distributions of Hydraulic Parameters for a Two-Scale Structured Heterogeneous Lignitic Mine Soil[J].Journal of Hydrology,2005,312(1):109-124.

[12]ANTONOPOULOS V Z,WYSEURE G C L.Modeling of Water and Nitrogen Dynamics on an Undisturbed Soil and a Restored Soil After Open-Cast Mining[J].Agricultural Water Management,1998,37(1):21-40.

[13]ZIPPER C,BALFOUR W,ROTH R,etal.Domestic Water Supply Impacts by Underground Coal Mining in Virginia,USA[J].Environmental Geology,1997,29(1/2):84-93.

[14]DECK O,ALHEIB M,HOMAND F.Taking the Soil-Structure Interaction Into Account in Assessing the Loading of a Structure in a Mining Subsidence Area[J].Engineering Structure,2003,25(4):435-448.

[15]COSH M H,JACKSON T J,MORAN S,etal.Temporal Persistence and Stability of Surface Soil Moisture in a Semi-Arid Watershed[J].Remote Sensing of Environment,2008,112(2):304-313.

[16]COSH M H,JACKSON T J,STARKS P,etal.Temporal Stability of Surface Soil Moisture in the Little Washita River Watershed and Its Applications in Satellite Soil Moisture Product Validation[J].Journal of Hydrology,2006,323(1):168-177.

[17]CHO E,CHOI M.Regional Scale Spatio-Temporal Variability of Soil Moisture and Its Relationship With Meteorological Factors Over the Korean Peninsula[J].Journal of Hydrology,2014,516:317-329.

[18]劉效東,喬玉娜,周國逸.土壤有機質對土壤水分保持及其有效性的控制作用[J].植物生態(tài)學報,2011,35(12):1209-1218.

[19]鄭紀勇,邵明安,張興昌.黃土區(qū)坡面表層風沙土容重和飽和導水率空間變異特征[J].水土保持學報,2004,18(3):53-56.

[20]鐘國輝,田發(fā)益，旺姆，等.西藏主要農區(qū)農田土壤肥力狀況研究[J].土壤,2005,37(5):523-529.

[21]黃昌勇,徐建明.土壤學[M].3版.北京:中國農業(yè)出版社,2012:56-70.

[22]李志洪,趙蘭坡,竇森.土壤學[M].北京:化學工業(yè)出版社,2005:85-92.

[23]馮杰,尚熳廷,劉佩貴.大孔隙土壤與均質土壤水分特征曲線比較研究[J].土壤通報,2009,40(5):1006-1009.

(責任編輯： 李祥敏)

Effect of Coal Mining on Soil Moisture Content in Semi-Arid Areas.

CHENGLin-sen,LEIShao-gang,BIANZheng-fu

(School of Environment Science and Spatial Informatics， China University of Mining and Technology, Xuzhou 221008, China)

Abstract:Soil moisture content is an important factor affecting vegetation growth in semi-arid areas. The study on effects of coal mining on soil moisture dynamics may help reestablish vegetation in the mining area.Profile soil moisture probes were used to monitor changes in soil moisture content (SMC) in the 10 to 100 cm soil layer of the compression area located at the 52303 working face of the Daliuta coal mine before and after the mining operation. Analysis of the monitoring data shows that before the mining operation, SMC remained quite stable and didn′t change much with temperature. In the initial period of the coaling operation SMC changed, but the change was not significant and lagged temporally behind the process of mining. The time lag reached 4-5 days, that is to say, the change in SMC became apparent 4 to 5 days later after the mining operation reached right beneath the monitoring site. Then subsidence began to develop till maximum. During the course it disturbed the soil mass, causing changes in SMC, which rose first and then fell. The disturbance of mining to SMC waned gradually with soil depth in the soil profile from top to bottom. The disturbance of the soil mass caused by subsidence also led to decrease in soil particle size, increase in bulk density and decline of void ratio, thus enhancing the water holding capacity of the soil mass, which explains why SMC rises temporarily in the initial period of the mining operation. Meanwhile the occurrence of ground fissures and reduced water supply from rainfall are the two major causes that lower SMC in the late period of the operation.

Key words:semi-arid area;coal mining;soil moisture;time sequential monitoring

作者簡介：程林森(1990—),男,安徽桐城人,碩士生,研究方向為礦山生態(tài)修復。E-mail： chenglinsen@126.com

DOI：10.11934/j.issn.1673-4831.2016.02.008

中圖分類號：X825;S152.7

文獻標志碼：A

文章編號：1673-4831(2016)02-0219-05

通信作者①E-mail: lsgang@126.com

基金項目：國家科技基礎性工作專項(2014FY110800);教育部新世紀優(yōu)秀人才支持計劃(NCET-12-0956)

收稿日期：2015-11-10