煤炭是当今世界社会运行的3大支柱能源之一，在我国一次性能源的生产和消耗中，煤炭占70%以上^[1]，且煤炭的能源支柱作用短期内不会改变^[2]。大规模的煤炭开采改变了矿区水文地质条件，造成大面积塌陷，出现地裂缝、塌陷坑、塌陷洞等，严重破坏土体结构，影响土壤水分布特征^[3]，水土保持功能降低，改变原有地形地貌，加剧降雨形成的水流对土壤的侵蚀^[4]，人为造成水土流失。神府-东胜煤田是我国重要的煤炭基地，近年来煤炭资源的过度开采，给当地生态环境造成了严重影响，其中采煤塌陷造成的水土流失尤为严重^[5-6]。部分学者^[7-10]通过对比采煤塌陷区和非塌陷区的土壤含水量，分析了塌陷区土壤含水量分布特征，但少见针对塌陷区内裂缝对土壤含水量影响的研究。笔者以神府-东胜煤田石圪台采煤塌陷区为研究对象，分析不同裂缝宽度下的土壤水特征，并对采煤塌陷裂缝区与非裂缝区的土壤含水量空间变异特征进行分析，以期了解采煤塌陷区裂缝对土壤水的影响，为矿区土地复垦及生态环境治理提供基础理论依据。

研究区位于神府-东胜煤田石圪台矿区^[11]。石圪台煤矿东西长10 km，南北宽8.5 km，矿区面积65.25 km²，地质储量8.35亿t，其中可开采量6.62亿t。矿区处于陕北黄土高原北部和毛乌素沙漠东南缘，植被覆盖度低，土壤以风沙土为主。地貌单元主要分为2大类：北部风积沙区，沙丘连绵，地形相对平缓；南部丘陵沟壑区，梁峁相间分布，地形东高西低^[12]。矿区内主要水系为乌兰木伦河，年径流量0.38亿~5.25亿m³。气候属北温带干旱半干旱大陆性气候，四季冷热多变，日夜温差悬殊，多年平均降雨量约为400 mm，蒸发量约为2 000 mm^[13]。

2016年11月在石圪台煤矿塌陷区选取100 m×100 m区域作为试验区，区内包括裂缝区、有植被覆盖的非裂缝区和无植被覆盖的非裂缝区。在实验区内每隔20 m作为1个步长均匀布点，布点36个(见图 1)。非裂缝区布点26个，其中有植被区7个，无植被区19个，裂缝区布点10个。根据塌陷裂缝边界及植被区与无植被区分界情况，补充10个测点(B1~B10)。在塌陷区不同宽度裂缝处(8、4、3、2和1 cm)和非裂缝区各取样1个，共6个进行土壤机械组成及含水量分析。在每个实验点0~60 cm深度内的0、10、20、30、40、50和60 cm等7个土层分别用环刀取样，每层重复取样3个，最后实验室测定结果取平均值。共采取土壤样品1 092个置于铝盒中，现场用精度为0.01 g的电子秤称量密封并记录数据。

图 1(Fig. 1)

图 1 实验区采样布点图 Fig. 1 Sampling layout map in the experimental area

土壤含水量测定采用烘干法^[12]；采用传统统计学分析土壤含水量的垂向变化特征^[14]；利用P-P检验图法^[15]，通过SPSS 20对数据进行正态性检验。在此基础上，利用地统计学中的Kriging插值法，通过Surfer11.0软件进行土壤含水量的平面变化特征分析。

不同宽度裂缝下不同深度土层土壤含水量测定结果见表 1。对照点与各裂缝点土壤含水量在土层深度上变化的变异系数值范围为0.33~0.46，均呈中等变异程度^[16](0.1 < C_V < 1)。各裂缝点不同深度土层的含水量均值全部小于对照点，且随着裂缝宽度增加，土壤含水量呈下降趋势，其均值依次降低为2.96%、2.78%、2.24%、1.45%和1.18%，分别减少12.68%、17.99%、33.92%、57.22%和65.19%。说明裂缝的存在对土壤含水量产生明显影响，且随着裂缝宽度的增加，影响作用加强，土壤含水量减小幅度不断增大。当地表裂缝宽为8 cm时，其土壤含水量均值仅为对照点的35%。对于相同深度的土层，裂缝宽度越大，土壤含水量越小。将对照点与不同宽度裂缝的同一深度土层的含水量差异进行显著性检验分析，发现在0~60 cm范围，各土层的含水量差异的显著性水平P值在0.003~0.045之间，均小于0.05，说明土壤裂缝宽度对土壤含水量的影响显著。

表 1(Tab. 1)

表 1 不同宽度裂缝条件下土壤含水量统计值 Tab. 1 Statistical value of soil moisture content under different fracture width % 裂缝宽度 Crack width/cm 土层深度Depth/cm 平均值 Average value 变异系数 Coefficient of variation 0 10 20 30 40 50 60 1 1.48 1.85 2.65 3.22 3.47 3.94 4.08 2.96 0.34 2 1.34 2.02 2.26 3.06 3.19 3.67 3.92 2.78 0.34 3 0.64 1.23 1.82 2.45 3.01 3.11 3.41 2.24 0.46 4 0.60 1.15 1.36 1.73 1.54 1.68 2.07 1.45 0.33 8 0.55 0.88 1.02 1.43 1.27 1.34 1.76 1.18 0.34 对照点Contrast point 1.79 1.94 2.87 3.89 3.91 4.34 4.98 3.39 0.36 P 0.003 0.045 0.007 0.003 0.014 0.016 0.004

表 1 不同宽度裂缝条件下土壤含水量统计值 Tab. 1 Statistical value of soil moisture content under different fracture width

从垂向分布看，对照点及裂缝宽为1 cm、2 cm和3 cm的测点，土壤含水量随土层深度增加而不断增大。而当裂缝宽为4 cm和8 cm时，地面以下30 cm内，土壤含水量较小，30~60 cm范围，各土层间土壤含水量变幅较小，且明显小于裂缝宽度为1 cm、2 cm和3 cm的相应深度的土壤含水量。说明当裂缝宽度较大时，30 cm以下土壤含水量受裂缝影响较大。

由于土壤结构的非均质性以及含水量受土壤结构、质地、降雨等因素的影响，土壤含水量在垂向上存在变化特征。采煤塌陷区土壤塌陷和大量裂缝的出现，破坏土体结构，进而影响土壤的孔隙度，对土壤含水量的分布产生较大影响^[17]。从表 2和表 3裂缝区和非裂缝区0~60 cm不同深度下的土壤含水量特征值可以看出，裂缝区土壤含水量在0~30 cm范围随着深度增加增幅显著，30 cm以下，增幅较小。最大值均处于60 cm处，最小值处于土壤表层(0 cm)。裂缝区和非裂缝区土壤含水量的均值和中位值变化趋势基本一致，随着土层深度的增加而变大。地表(0 cm)土壤含水量明显低于地下各层(10~60 cm)，主要因为地表土壤，在冬季(11月)降雨补给少，同时实验区地处干旱半干旱地区蒸发作用强烈所致；而地表以下(10~60 cm)的土壤层，历经雨季后水分下渗到土壤深层，地表以下土壤中蒸发作用减弱，水分蒸发量降低。同一土层深度，裂缝区土壤含水量小于非裂缝区，是由于裂缝的存在破坏土壤原生结构，出现了水分优先路径，使得水分更快向土壤深处运动，降低了土壤持水能力，并且裂缝增加土壤与外界的接触面积，使得蒸发量增加，导致裂缝区土壤含水量偏低。

表 2(Tab. 2)

表 2 裂缝区土壤含水量特征值统计 Tab. 2 Eigenvalue of soil moisture content in crack area 特征值Eigenvalue 深度Depth/cm 0 10 20 30 40 50 60 变化范围Range of change/% 0.43~1.05 1.12~1.86 1.86~3.35 1.51~3.39 1.04~3.78 1.47~4.48 1.03~4.63 平均值Average value/% 0.57 1.43 2.39 2.47 2.62 2.56 2.74 中位值Median value/% 0.58 1.39 2.32 2.46 2.59 2.41 2.92 标准差Standard deviation 0.31 0.27 0.44 0.53 0.85 0.83 1.04 离散系数Discrete coefficient 0.54 0.19 0.18 0.21 0.32 0.33 0.38

表 2 裂缝区土壤含水量特征值统计 Tab. 2 Eigenvalue of soil moisture content in crack area

表 3(Tab. 3)

表 3 非裂缝区土壤含水量特征值统计 Tab. 3 Eigenvalue of soil moisture content in non-crack area 特征值Eigenvalue 深度Depth/cm 0 10 20 30 40 50 60 变化范围Range of change/% 1.12~4.39 2.29~7.09 2.29~8.44 3.17~8.67 3.06~8.39 3.59~6.48 4.17~11.41 平均值Average value/% 1.92 3.56 4.39 4.63 4.51 4.98 5.71 中位值Median value/% 1.53 2.98 3.86 3.86 4.37 4.87 5.25 标准差Standard deviation 0.85 1.42 1.51 1.63 1.11 0.79 1.52 离散系数Discrete coefficient 0.44 0.39 0.34 0.35 0.24 0.16 0.27

表 3 非裂缝区土壤含水量特征值统计 Tab. 3 Eigenvalue of soil moisture content in non-crack area

离散系数是指标准差与平均数的比例^[18]，它反映单位均值上的离散程度，离散系数越大，说明其离散程度越大。对于土壤含水量垂向分布的离散程度，采用离散系数来度量。裂缝区土壤含水量不同深度的离散系数在0.18~0.54之间，非裂缝区处于0.16~0.44之间。裂缝区与非裂缝区土壤含水量垂向不同深度的离散程度均为中等变异性。受自然因素与人为因素的影响，裂缝区与非裂缝区在土壤表层(0 cm)离散系数均为最大。

地统计学是以半方差函数为工具，研究在空间分布上既有随机性又有相关性的自然现象的科学。P-P图是正态分布检验的一种方法，本研究利用SPSS20对土壤含水量进行P-P图检验分析。从图 2可以看出，各个深度土壤含水量P-P检验图近似为直线型，可以看作正态分布。对数据进行半方差函数拟合，分析模型误差^[19-20]。从表 5拟合模型的参数误差可知满足分析要求。利用拟合模型通过Surfer 11.0中的Kriging插值对土壤含水量空间分布进行分析。实验区0~60 cm不同深度土层含水量空间分布见图 3~9。

图 2(Fig. 2)

图 2 不同深度土壤含水量P-P正态检验图 Fig. 2 P-P normal distribution test diagram of soil moisture content at different depths

表 4(Tab. 4)

表 4 不同深度土壤含水量半方差函数拟合模型的参数误差 Tab. 4 Parameter errors from semi-variogram functional model of soil moisture content at different depths 深度 Depth/cm 均方根预测误差 Predicted error of root mean square 平均标准误差 Mean standard error 标准化平均误差 Standardized mean error 标准化均方根预测误差 Predicted error of normalized mean square root 0 0.374 5 0.372 9 0.012 8 0.997 8 10 0.124 7 0.108 2 0.007 6 1.083 8 20 0.283 2 0.279 5 0.002 1 1.065 1 30 0.304 1 0.286 7 0.012 6 0.983 5 40 0.178 2 0.152 9 0.007 5 0.974 4 50 0.302 5 0.344 3 0.017 2 1.029 0 60 0.388 5 0.377 0 0.021 3 0.936 4

表 4 不同深度土壤含水量半方差函数拟合模型的参数误差 Tab. 4 Parameter errors from semi-variogram functional model of soil moisture content at different depths

图 3(Fig. 3)

图 3 0 cm土壤含水量等值线图 Fig. 3 Contour map of soil moisture content at the depth of 0 cm

图 4(Fig. 4)

图 4 10 cm土壤含水量等值线图 Fig. 4 Contour map of soil moisture content at the depth of 10 cm

图 5(Fig. 5)

图 5 20 cm土壤含水量等值线图 Fig. 5 Contour map of soil moisture content at the depth of 20 cm

图 6(Fig. 6)

图 6 30 cm土壤含水量等值线图 Fig. 6 Contour map of soil moisture content at the depth of 30 cm

图 7(Fig. 7)

图 7 40 cm土壤含水量等值线图 Fig. 7 Contour map of soil moisture content at the depth of 40 cm

图 8(Fig. 8)

图 8 50 cm土壤含水量等值线图 Fig. 8 Contour map of soil moisture content at the depth of 50 cm

图 9(Fig. 9)

图 9 60 cm土壤含水量等值线图 Fig. 9 Contour map of soil moisture content at the depth of 60 cm

从表层(0 cm)土壤含水量等值线图(图 3)可以看出，该层土壤含水量总体较低。裂缝区土壤含水量范围在0.43%~1.05%之间，平均0.57%，离散系数0.54，属于中等变异强度。非裂缝区土壤含水量范围在1.12%~4.39%之间，平均1.92%，离散系数0.44，属于中等变异, 变异性小于裂缝区(表 2和表 3)。在实验区分别存在2个土壤含水量高值区和低值区，高值区1个位于实验区西北角，1个位于东部偏南，低值区1个位于试验区中部，1个位于东北角。结合实验区布点情况，可以看出，土壤含水量高值区均位于非裂缝区，同时位于西北角的高值区土壤含水量相对较高且等值线较密，低值区1个位于裂缝宽度较大的区域，另1个位于无植被区但靠近宽度较大的裂缝处。说明采煤塌陷裂缝使得土壤表层含水量减少，且宽度越大处，含水量越低，植被根系对土壤水具有一定的固结作用。

从10 cm深度土壤含水量分布等值线图(图 4)可以看出，该层土壤含水量总体上比地表(0 cm)处高。裂缝区与非裂缝区土壤含水量均属于中等变异, 非裂缝区变异性大于裂缝区。分布有3个土壤含水量高值区，全部位于非裂缝区，且等值线分布较密，说明在高值区含水量平面变幅较大；而土壤含水量低值区位于裂缝区及无植被区接近裂缝区，其中最低值位于裂缝宽度较大处。

20 cm土层含水量分布等值线如图 5所示。该层土壤含水量总体上比10 cm土层高。塌陷区土壤含水量变异性属于中等变异，但裂缝区土壤含水量变异性小于非裂缝区。位于植被覆盖区的土壤含水量高值区，等值线较为密集，说明含水量平面变幅较大，而位于塌陷区的低值区等值线较为稀疏，说明含水量平面变幅较小。同时，可以看到土壤含水量高值区相对于10 cm分布图向北移动，但最高区依然位于非裂缝区的植被覆盖区，说明植物对增加土壤含水量具有较大作用。

30 cm土层含水量分布等值线见图 6。非裂缝区土壤含水量与20 cm土层基本一致，裂缝区土壤含水量有较小幅度的增大。裂缝区与非裂缝区的土壤含水量变异性均属于中等变异。

土壤含水量高值区依然位于非裂缝区，且最高区位于植被覆盖区，最低区位于裂缝区与无植被区的结合处，且靠近宽度较大的裂缝。

40 cm土层含水量分布等值线见图 7。非裂缝区土壤含水量与30 cm土层基本一致，裂缝区土壤含水量有较小的增幅。裂缝区与非裂缝区土壤含水量均属于中等变异。该层土壤含水量最小区位于采煤塌陷裂缝区中裂缝宽度较大处，裂缝宽度达4 cm。说明裂缝宽度对土壤含水量的影响较为强烈，且裂缝宽度越大，垂直影响深度越深。

50 cm土层含水量分布等值线如图 8所示。与上层相比，该层土壤含水量均值依然增加。裂缝区与非裂缝区土壤含水量均属于中等变异。3个高值区均位于非裂缝区，低值区位于宽裂缝区或靠近宽裂缝的无植被区。

60 cm土层含水量分布等值线见图 9。裂缝区与非裂缝区土壤含水量均属于中等变异。土壤含水量最大值为11.41%，极差为10.38%，土壤含水量平均值为4.55%。与0~50 cm土层对比，土壤含水量最大值、极差和平均值均达到了最大，主要因为该层土壤深度较深，受蒸发蒸腾作用和人为影响较小，土壤含水量损失较少。土壤含水量高值区在植被覆盖区。2处含水量低值区均位于裂缝区中裂缝宽度较大或裂缝较发育的位置，土壤含水量损失严重。

1) 裂缝的存在会使相同深度土层的土壤含水量减少，且随着裂缝宽度的增加，土壤含水量减小幅度不断增大。当裂缝宽度>3 cm时，30 cm以下土壤含水量受裂缝影响较大。

2) 在干旱半干旱地区强烈的蒸发作用和冬季降雨补给较少的情况下，地表(0 cm)土壤含水量明显偏低；而地表以下(10~60 cm)的土壤层含水量明显高于地表(0 cm)。

3) 裂缝的存在降低了土壤持水能力；并且裂缝的存在，增加了土壤与外界的接触面积，增加蒸发量，导致了同一土层深度裂缝区土壤含水量明显小于非裂缝区。

4) 采煤塌陷区裂缝对土壤含水量的平面分布有强烈的影响。0~60 cm各土层土壤含水量高值区总是处在无裂缝区的植被覆盖区，低值区处在裂缝宽度较大的区域或无植被区靠近宽裂缝处。