Analysis of ecological water level recovery degree under coal mining subsidence disturbance in ecologically fragile area
-
摘要:
生态脆弱区生态水位变异程度对区域生态地质环境保护至关重要。煤层采动下生态水位势必发生变异，当前研究集中于采动潜水渗漏下生态水位下降，极少考虑潜水不渗漏/采煤沉陷扰动下生态水位恢复程度。为此，采用地下水动力学、地下水位实测、数理统计、综合分析等方法，提出以“判别煤层采动下潜水渗漏状态→实测采煤工作面生态水位变化→建立生态水位恢复程度井流解析模型→预测不同生态水位恢复时间→对比不同生态水位恢复时间解析值与实测值”为思路的解析方法，研究生态脆弱区采煤沉陷扰动下生态水位恢复程度。结果表明：①基于“关键层位置+薄板理论+土拱效应+下行裂隙”建立了采动覆岩–土结构下导水裂隙发育高度计算方法，克服了现有经验公式未考虑土层效应导致的预测精度不足，结合残余隔水层厚度与潜水渗漏状态，实现了采动潜水渗漏状态判别；②采煤沉陷扰动下实测生态水位呈现“先迅速下降–缓慢回升–趋于稳定”的变化规律，但采后生态水位不能完全恢复至采前状态；③建立了采煤沉陷扰动下生态水位恢复程度井流解析模型，预计生态水位3个不同恢复程度的恢复时间解析值，对比于实测值，两者时间误差均小于10%；④从地表地形地貌、大气降雨补给、潜水含水层补径排、矿区井下疏放水等角度探讨了采后生态水位未完全恢复至采前状态原因。
Abstract:The variation degree of ecological water level in ecologically fragile areas is very important for ecological geological environment protection. Ecological water level inevitably undergoes variation under coal seam mining. The current research focuses on the decline of ecological water level under mining induced groundwater leakage, with little consideration given to the ecological water level restoration degree under non leakage of groundwater or disturbance of coal mining subsidence. For this reason, methods such as groundwater dynamics, groundwater level measurement, mathematical statistics, and comprehensive analysis are adopted, and an analytical method was proposed based on the idea of “distinguishing the leakage status of groundwater under coal seam mining → measuring the ecological water level changes in the coal mining face → establishing an ecological water level restoration degree well flow analysis model → predicting different ecological water level restoration times → comparing the analytical and measured values of different ecological water level restoration times” to study the degree of ecological water level restoration under the disturbance of coal mining subsidence in ecologically fragile areas. The results are as follows. ① Based on “key strata layer location + thin plate theory + soil arch effect + descending fissure”, a calculation method for the water-conducting fissure development height under overburden bedrock and soil structure is established, which overcomes the insufficient prediction accuracy caused by the existing empirical formula not considering the soil layer effect. Furthermore, combines the relationship between the thickness of the residual water-retaining layer and the aquifer seepage state, the aquifer seepage state is confirmed. ② Under the disturbance of coal mining subsidence, the measured ecological water level shows a change law of “rapid decline first, slow recovery, and stabilization”, but the ecological water level after coal mining can’t completely recover to the premining state. ③ An analytical model of the recovery degree of ecological water level under the disturbance of coal mining subsidence is established. The analytical values of the recovery time of three different recovery degrees of ecological water level are estimated, and the time errors of both are less than 10% compared with the measured values. ④ The reasons for the incomplete recovery of postmining ecological water level to pre-mining state are discussed from the perspectives of surface topography, atmospheric rainfall recharge, recharge and drainage of phreatic aquifers, and underground water drainage in mining areas.
-
0. 引 言
西北生态脆弱区具有显著的“富煤、贫水、弱生态”特征。第四系萨拉乌苏组潜水是该区域惟一具有供水和生态意义水资源。潜水生态水位与区域植被生态、地质环境密切相关[1]。煤炭资源大规模开发，势必扰动生态水位变异，诱发原本脆弱的生态环境进一步恶化。因此，研究生态脆弱区煤层采动下生态水位变异程度对矿区生态地质环境保护、绿色矿山建设具有重要意义。
保水采煤的目标是在防治采场突水的同时，维持具有供水意义和生态价值的含水层稳定或将生态水位变化控制在合理范围内。为此，煤层采动下生态水位变异研究取得了系列成果。现场实测方面，结合地下水传感器监测钻孔潜水水位，分析了潜水位下降与煤层开采强度的关系[2]、潜水位波动与地表沉陷的耦合关系[3]、潜水位下降对干旱矿区植被蒸腾[4]及河岸区地下水蒸散发量的影响[5]、地下水位与干旱区生态系统关系[6]、保水采煤环境工程地质模式及保水采煤类型区域研究[7]。钻孔水位监测较难实现大区域全方位、小尺度、精细化监测，为此，结合数值模拟，分析了采矿活动和降雨事件对潜水位的影响规律[8–10]、预测了矿井涌水量及评价了矿区地下水资源[11–13]、研究了高强度煤层采动下最佳生态水位将至警戒水位的生态效应[14]，弥补了钻孔水位监测的不足。此外，基于偏差信息和非线性自回归神经网络混合模型[15]、创新趋势分析和Mann-Kendall方法[16]、频谱法和校正方法、时间序列分析法[17]等数学理论方法，在选择多个影响因素的基础上，预测采矿活动对地下水位的影响，但缺乏水位变化的物理意义。采矿活动影响下地下水演变机理方面，Abdullah Karaman等假设将采煤工作面看作移动抽水井，模拟了地面沉降对地下水位的影响[18]；随后采用类似曲线分析计算法建立了采矿速率、水力扩散系数与潜水水位的耦合关系[19]，探究了采矿活动影响下地下水位的演变规律。杨倩等[20]阐述了采动裂隙导水和变形压力作用下的承压水位变化机制，建立了地下水稳定运动、非稳定运动2种井流方程的覆岩承压含水层水位采动变化数学模型。李涛等[21]利用Theis公式与相似模拟试验相结合的方法，建立了煤层采动对地下水位恢复的井流预测模型。
当前研究集中于煤层采动潜水渗漏条件下生态水位如何下降，而往往忽略潜水不渗漏/采煤沉陷扰动下生态水位恢复演变规律。此外，采动生态水位演变往往依靠现场实测，水位恢复程度预测解析解不明确。基于此，笔者首先基于“关键层位置+薄板理论+土拱效应+下行裂隙”建立煤层采动下覆岩–土隔水层厚度计算方法，判别潜水渗漏状态；其次，结合2个工作面的潜水生态水位实测数据，研究煤层采动下潜水生态水位变化规律；然后，基于井流模型建立采煤沉陷扰动下生态水位恢复程度解析解；最后，对比分析生态水位恢复实测值与解析值，并探讨了生态水位未完全恢复的原因。
1. 煤层采动下潜水渗漏状态判别
1.1 覆岩残余隔水层厚度的确定方法
随着我国东部煤炭资源逐渐枯竭，煤炭开采重心已转移至西北生态脆弱区。西部地区主采侏罗系煤层，地表被第四系沙层所覆盖，其下的第四系萨拉乌苏组潜水含水层具有重要的生态意义。第四系潜水隔水层主要是新近系保德组红土与中更新统离石组黄土，西部地区地层综合柱状图如图1所示。
![]()
图 1 西部矿区地层综合柱状图Figure 1. Comprehensive column chart of strata in the western mining area煤层采动后，形成覆岩导水裂隙带。目前煤层开采导水裂隙带发育高度(简称“导高”)预计的方法，主要针对我国东部石炭–二叠系煤层开采实测值总结。东部矿区主要以深井开采为主，而西部矿区开采主要为浅埋厚煤层，具有大采高、大采深的特点，上覆潜水隔水土层主要以黄土、红土为主，导水裂隙带易贯穿基岩进入(部分)土层，与东部矿区煤层埋藏以及开采条件大有不同，先前总结的经验公式难以实现西部矿区导高的准确预测[22–24]。
基于上述分析，针对西部矿区特征，结合已有研究成果，基于“关键层位置+薄板理论+土拱效应+下行裂隙”建立了覆岩–土结构下采动导高计算模型(图2)。首先，当覆岩当关键层位置距开采煤层大于7～10倍煤层采厚距离时，导高等于该关键层距开采煤层的高度[25]。反之，基于薄板理论计算主关键层上部基岩破断情况，即：通过对比薄板的极限挠度值与岩层下部自由空间高度的大小，判断各基岩岩层是否全部断裂，确定发育在基岩中的导高。若上部基岩岩层全部断裂，则需要判断导水裂隙带是否继续向上发育至土层。为此，基于普氏理论和岩体极限平衡理论确定土体破坏临界高度[26]。如果土层的厚度足够大，当土层被破坏到临界高度时，就会形成一个自然平衡拱，此时土层就不会继续垮落，导水裂隙带发育高度即为临界高度与基岩厚度之和。如果土层没有形成稳定拱，考虑采煤沉陷下行裂隙深度是否发育至基岩顶界面[27]，若两者贯通，导水裂隙带贯穿地表；若没有发育至基岩则导水裂隙带发育到土–岩交界面。综上，确定覆岩–土结构下采动导水裂隙带发育高度，为下述潜水渗漏状态判别提供理论依据。
![]()
图 2 覆岩–土结构下采动导水裂隙带高度计算流程Figure 2. Flow of calculating the height of water fracture zone under overburden bedrock and soil structure1.2 考虑残余隔水层厚度的潜水渗漏状态判别
在确定覆岩–土结构下采动导高的基础上，确定覆岩残余隔水层厚度及其类型。为保证煤矿开采安全，通常根据松散层、黏性土层厚度以及采厚等因素留设防水或防砂安全煤岩柱，但是对于考虑残余隔水层厚度及其岩性组合的潜水渗漏状态无法判别。
为此，对于覆岩残余隔水层厚度与渗漏状态的确定，笔者在分析潜水渗漏模式阈值和残余隔水层厚度组合线性变化规律的基础上，确定潜水渗漏与否的阈值：红土、黄土和基岩厚度分别为22、36和120 m[26]。李涛等[28]采用水–电相似模拟得出：离石黄土42.6 m或保德红土21.0 m为潜水不发生渗漏的最小厚度。综上，对比分析残余隔水层厚度与渗漏阈值(表1)，确定潜水的渗漏状态。
表 1 潜水渗漏模式判断标准Table 1. Criteria for judging aquifer leakage modes渗漏情况 判断方法1 判断方法2 潜水不发生渗漏 $ \dfrac{x}{{22}} + \dfrac{y}{{36}} + \dfrac{{\textit{z}}}{{120}} - 1 \geqslant 0 $ $ \dfrac{x}{{42.6}} + \dfrac{y}{{21}} - 1 \geqslant 0 $ 潜水发生渗漏 $ \dfrac{x}{{22}} + \dfrac{y}{{36}} + \dfrac{{\textit{z}}}{{120}} - 1 < 0 $ $ \dfrac{x}{{42.6}} + \dfrac{y}{{21}} - 1 < 0 $ 注：x为残余红土厚度，m；y为残余黄土厚度，m；z为残余基岩厚度，m。 2. 采动沉陷扰动下生态水位实测
2.1 煤层采动潜水渗漏状态
在西北生态脆弱区金鸡滩煤矿108工作面和小保当煤矿01工作面分别开展采煤扰动下生态水位实测。为判断潜水渗漏状态，在108工作面布置了KY1、KY2两个钻孔实测煤层采动导高分别为170.55、178.45 m[14](图3)。在钻孔KY1中，残余基岩和土层的厚度分别为9.26、39.7 m。在钻孔KY2中，残余基岩和土层的厚度分别为19.16、26.8 m。对于01工作面，布置了钻孔S5测试导高为157.39 m[3](图4)，残余隔水基岩和土层的厚度分别为74.07、65 m。
此外，结合覆岩残余隔水层厚度的确定方法和工作面综合钻孔柱状图，确定01、108工作面覆岩主关键层与开采煤层的距离均大于7～10倍煤层采厚，基于此，导高等于该关键层距开采煤层的高度分别为：179.13、191.76 m。结合工作面综合柱状图，108工作面残余基岩和土层的厚度分别为10.58、33.25 m；01工作面残余基岩和土层的厚度分别为52.33、65 m。
结合上述潜水渗漏状态参与隔水层厚度阈值，确定01、108工作面煤层采动过程中潜水均未发生渗漏。
2.2 煤层采动下生态水位实测
为了监测工作面整个采煤期间潜水位变化，在108工作面设置了水位监测钻孔KY8和KY10(图3)。监测钻孔KY8位于工作面倾斜方向的中心，距离工作面终采线170 m；监测钻孔KY10位于108工作面运输巷上方，距离工作面终采线117 m (图3)。01工作面布置S5和S13三个水文监测钻孔(图4)。同时，监测各个钻孔孔口的地面沉降值。
KY8钻孔于2017–12–06开始监测，此时采煤工作面煤壁与监测钻孔KY8的水平距离为3.4 m。根据图5可将KY8水位及地面沉降变化分为4个阶段，分别为1—水位快速下降阶段；2—受地面迅速沉降影响水位短暂平稳阶段；3—地面沉降幅度减弱、水位缓慢上升阶段；4—地面沉降幅度进一步减小直至沉降稳定、水位逐渐稳定阶段。工作面过孔后24 d左右水位下降至最低点，采煤进尺远离钻孔并运动到终采线后钻孔处沉降幅度逐渐减小并趋于稳定，水位逐渐恢复。
01工作面S5钻孔于2019–04–10开始观测，潜水位随回采均表现出“先下降后回升”的规律，如图6所示。采煤工作面过孔后25 d左右各水文孔达到最低水位，与地表下沉具有同时性，即地表开始下沉时，潜水位开始下降，下降幅度基本一致。采后138 d左右采煤进尺远离钻孔并运动到终采线后钻孔处沉降幅度逐渐减小并趋于稳定，水位逐渐恢复[3]。
2.3 生态水位演变规律
根据01和108工作面采动潜水变化实测可知，对于某一观测点来说，地面沉降与水位下降同步，地面沉降的幅度将影响潜水位的下降速度：地表沉降活跃期，地表沉降速度急剧加快，甚至大于水位下降幅度，潜水位迅速下降；当地表沉降进入活跃阶段后期，地面沉降幅度逐渐减弱，潜水位逐渐稳定；地面沉降衰退阶段，地表沉降趋于稳定，潜水位易在大气降水和潜水侧向补给作用下逐渐回升，进入回升阶段。
因此，在煤层采动潜水不渗漏/采煤沉陷扰动条件下，生态水位呈现“迅速下降—缓慢回升—趋于稳定”的演化规律，但采后潜水水位通常不能完全恢复至采前水位。
3. 采煤沉陷扰动下生态水位恢复程度井流模型解析
为模拟测孔的潜水位在采煤后恢复所需要的时间，建立满足Theis公式的单井数学模型。将模拟恢复时间与上述实际测量恢复时间进行对比，研究采煤对水位波动的影响规律。
3.1 井流模型的建立
随着采煤进尺经过监测点并运动到停采线水位下降，将此过程的瞬时水位看作一口虚拟抽水井(图7)，以定流量Q在tp时间内抽水造成潜水位降落漏斗；而潜水位的恢复发生在采煤结束地面沉降稳定停止虚拟抽水后，将这一过程看作虚拟注水。因此潜水位降深s′与时间t的关系如图8所示，图中t′为水位恢复的时间，sr′为水位恢复期的修正剩余降深[29-30]。
3.2 井流模型的解析
数学模型的建立采用定流量的Theis公式
$$ s' = \frac{Q}{{4\pi T}}W(u) $$ (1) $$ W(u) = \int_u^\infty {\frac{{{{\mathrm{e}}^{ - y}}}}{y}} {\mathrm{d}}y $$ (2) $$ u = \frac{{{r^2}\mu }}{{4Tt}} $$ (3) $$ W(u) = \ln \frac{{2.25Tt}}{{{r^2}\mu }} $$ (4) 式(3)、(4)为Theis公式的简化形式。式中，s′为虚拟抽水井的降深，m；r为钻孔到虚拟抽水井的距离，m；u为含水层的给水度； T为含水层导水系数，m2/d，T=KM，K为含水层的渗透系数，m/d，M为初始含水层厚度，m；t为虚拟抽水开始到计算时刻的时间，d ；Q为虚拟抽水井的定流量，m3/d；W(u)为井函数；y为积分变量。
首先计算虚拟抽水流量Q，选取模型中非对称的两任意观测点，已知其虚拟抽水结束后的水位s1's2'，假定虚拟抽水井位置，预设r1 r2，结合井式(1)得到式(5)，即可以求出Q值，再把求出的Q值代入Jacob式(6)，得到近似的虚拟抽水进行的时间tp。
$$ K = \frac{{0.318\;3Q(\ln \;{r_2} - \ln \;{r_1})}}{{(2M - {s_1}' - {s_2}')({s_1}' - {s_2}')}} $$ (5) $$ s' = \frac{{0.183Q}}{T}\lg \frac{{2.25T{t_{\mathrm{p}}}}}{{{r^2}\mu }} $$ (6) 根据计算出的Q值和tp值代入恢复期的Thies式(7)，其等价于同流量的虚拟抽水和虚拟注水的叠加，并代入任意降深可以求出对应的恢复所用时间t'。
$$ {s_{\mathrm{r}}}' = \frac{Q}{{4\pi T}}\left(\ln \frac{{t' + {t_{\mathrm{p}}}}}{{t'}} - \frac{{{r^2}\mu }}{{4Tt'}}\right) $$ (7) 根据井流模型解析生态水位波动过程的流程如图9所示。
4. 采煤沉陷扰动下生态水位恢复程度预测
4.1 生态水位恢复实测值与预测值对比分析
根据两工作面已有的工程地质、水文地质参数，结合2.2节中的水位监测数据，使用3.2节中所述的方法，能够对108工作面和01工作面的生态水位恢复程度进行预测解析。
01工作面中，$ K = 3.5\;{\mathrm{m/d}}，M = 10\;{\mathrm{m}}，T = KM = 35\;{{\mathrm{m}}}^{2}/{\mathrm{d}}，\mu =0.27 $。选取2.2节中测量出的S5和S13钻孔的水位降深以及其距离虚拟井的水平距离$ s'_1=2.595\;{\mathrm{m}}，s'_2=1.679\;{\mathrm{m}}，{r}_{1}=120\;{\mathrm{m}}，{r}_{2}= 128.5\;{\mathrm{m}} $为观测数据代入式(5)计算出虚拟抽水井的流量Q=2 314.479 m3/d。将流量Q代入抽水阶段的Jacob式(6)得出虚拟抽水时间$ {t}_{p}=26.94\;{\mathrm{d}} $，S5观测井实际抽水时间为27 d (图10中的I阶段)。S5水位恢复到第1个峰值时的水位降深为1.635 m，将虚拟抽水井流量Q和虚拟抽水时间tp代入式(7)可以得出虚拟注水时间t'=27.37 d，实际恢复时间为26 d (图10中的II–1阶段)；恢复到最高水位时的水位降深为1.081 m，虚拟注水时间t'=64.26 d，实际恢复时间为63 d (图10中的II–2阶段)。
![]()
图 10 S5钻孔水位波动预测时间与实际时间对比Figure 10. Comparison between predicted time and actual time of water level fluctuation in S5 borehole108工作面KY8钻孔的生态水位恢复过程解析[30]，得出KY8观测井虚拟抽水时间$ {t}_{{\mathrm{p}}}=243.988\;{\mathrm{d}} $，实际抽水时间为244 d(图11中的I阶段)。KY8观测井虚拟水位恢复到60%时，得$ t'_{60\%}=23.715\;{\mathrm{d}} $；KY8观测井实际恢复时间为20 d左右(图11中的II–1阶段)；恢复到110%时$ s'_{1\tau }=2.267\;{\mathrm{m}} $， $ t'_{110\%}= 42.127\;{\mathrm{d}} $；实际为40 d左右(图11中的II–2阶段)；恢复到200%时，$ t'_{200\%}=78.178\;{\mathrm{d}} $；实际为80 d左右(图11中的II–3阶段)。
![]()
图 11 KY8钻孔水位波动预测时间与实际时间对比Figure 11. Comparison between predicted time and actual time of water level fluctuation in KY8 borehole以上模拟过程的计算结果与实测结果之间存在一定误差，主要来源于计算过程中对井公式的简化带来的误差，以及开采过程中地面沉降带来的上覆隔水层的变化会对结果造成影响。但是误差较小，因此可以通过建立井流模型与实际测量数据相结合来估计采煤后的水位恢复时间。
4.2 生态水位未完全恢复原因
综合现场监测和解析模型生态水位变化数据可知：潜水不渗漏情况下采动生态水位下降主要由采煤沉降导致。随着地面沉降速度下降并趋于0，生态水位会有一定程度恢复，但通常在相当长时间内不能完全恢复至采前生态水位。笔者认为主要由地表地形地貌、大气降雨补给、潜水含水层补径排、矿区井下疏放水等综合因素诱发采后生态水位未能完全恢复。
其中，地形地貌控制着潜水水位的分布特征，潜水面形态与地表形态基本一致，采煤扰动下地表地形、微地貌均会发生变化，尤其在采煤沉陷盆地范围内，生态水位也相应发生变化。若煤层采动区域内潜水含水层分布局限，侧向补给缺乏[3]，或大气降水补给条件较差，无法弥补由于沉陷引起的生态水位下降程度就会导致采煤沉陷扰动下生态水位无法完全恢复。为了防治矿井水害，井下通常采取钻孔疏放水的措施，易诱发区域地下水系统中各含水层出现不同程度的补给排泄，进而会导致生态水位未能恢复至采前状态。
5. 结 论
1)针对西部矿区特征，基于关键层位置+薄板理论+土拱效应+下行裂隙建立了采动覆岩–土结构下导水裂隙带高度计算模型，结合覆岩残余基岩和土层阈值，构建了煤层采动下潜水渗漏状态判别方法。
2)在煤层采动潜水不渗漏/采煤沉陷扰动条件下，生态水位呈现迅速下降→缓慢回升→趋于稳定的演化规律，但采后潜水水位往往不能完全恢复至采前水位。
3)建立了采煤沉陷扰动下生态水位恢复程度井流解析模型，对比实测结果表明：生态水位不同恢复程度的恢复时间解析值与实测值误差均小于10%。
-
![]()
图 1 西部矿区地层综合柱状图
Figure 1. Comprehensive column chart of strata in the western mining area
![]()
图 2 覆岩–土结构下采动导水裂隙带高度计算流程
Figure 2. Flow of calculating the height of water fracture zone under overburden bedrock and soil structure
![]()
图 10 S5钻孔水位波动预测时间与实际时间对比
Figure 10. Comparison between predicted time and actual time of water level fluctuation in S5 borehole
![]()
图 11 KY8钻孔水位波动预测时间与实际时间对比
Figure 11. Comparison between predicted time and actual time of water level fluctuation in KY8 borehole
表 1 潜水渗漏模式判断标准
Table 1 Criteria for judging aquifer leakage modes
渗漏情况 判断方法1 判断方法2 潜水不发生渗漏 $ \dfrac{x}{{22}} + \dfrac{y}{{36}} + \dfrac{{\textit{z}}}{{120}} - 1 \geqslant 0 $ $ \dfrac{x}{{42.6}} + \dfrac{y}{{21}} - 1 \geqslant 0 $ 潜水发生渗漏 $ \dfrac{x}{{22}} + \dfrac{y}{{36}} + \dfrac{{\textit{z}}}{{120}} - 1 < 0 $ $ \dfrac{x}{{42.6}} + \dfrac{y}{{21}} - 1 < 0 $ 注：x为残余红土厚度，m；y为残余黄土厚度，m；z为残余基岩厚度，m。 
下载: 导出CSV
-
[1] 池明波. 我国西北矿区水资源承载力评价与科学开采规模决策[D]. 徐州: 中国矿业大学,2019: 5–6. CHI Mingbo. Evaluation of water resources carrying cpacity and decision-making of scientific coal mining scale in Northwest Mining Area of China:a case study in Yining Mining Area[D]. Xuzhou:China University of Mining and Technology,2019: 5–6.
[2] 范立民,向茂西,彭 捷,等. 西部生态脆弱矿区地下水对高强度采煤的响应[J]. 煤炭学报,2016,41(11):2672−2678. doi: 10.13225/j.cnki.jccs.2016.0243 FAN Liming,XIANG Maoxi,PENG Jie, et al. Groundwater response to intensive mining in ecologically fragile area[J]. Journal of China Coal Society,2016,41(11):2672−2678. doi: 10.13225/j.cnki.jccs.2016.0243
[3] 侯恩科,谢晓深,王双明,等. 中深埋厚煤层开采地下水位动态变化规律及形成机制[J]. 煤炭学报,2021,46(5):1404−1416. HOU Enke,XIE Xiaoshen,WANG Shuangming, et al. Dynamic law and mechanism of groundwater induced by medium-deep buried and thick coal seam mining[J]. Journal of China Coal Society,2021,46(5):1404−1416.
[4] 王 皓,周振方,杨 建,等. 蒙陕接壤区典型煤层开采地下水系统扰动的定量表征[J]. 煤炭科学技术,2020,51(7):83−93. doi: 10.1007/s10230-020-00714-6 WANG Hao,ZHOU Zhenfang,YANG Jian, et al. Quantitative characterization of the disturbance of groundwater system in typical coal seam mining in contiguous area of Inner Mongolia and Shaanxi[J]. Coal Science and Technology,2020,51(7):83−93. doi: 10.1007/s10230-020-00714-6
[5] YUE Weifeng,WANG Tiejun,TRENTON E. Franz, et al. Spatiotemporal patterns of water table fluctuations and evapotranspiration induced by riparian vegetation in a semiarid area[J]. Water Resources Research,2016,52(3):1948−1960. doi: 10.1002/2015WR017546
[6] CHEN Weichi,LI Wenping,YANG Zhi, et al. Analysis of mining-induced variation of the water table and potential benefits for ecological vegetation:a case study of Jinjitan coal mine in Yushenfu mining area,China[J]. Hydrogeology Journal,2021,29(4):1−17.
[7] 李文平,王启庆,刘士亮,等. 生态脆弱区保水采煤矿井(区)等级类型[J]. 煤炭学报,2019,44(3):718−726. doi: 10.13225/j.cnki.jccs.2018.6031 LI Wenping,WANG Qiqing,LIU Shiliang, et al. Grade types of water-preserved coal mining coalmines in ecologically fragile area[J]. Journal of China Coal Society,2019,44(3):718−726. doi: 10.13225/j.cnki.jccs.2018.6031
[8] XUE Sen,LIU Yu,LIU Shiliang, et al. Numerical simulation for groundwater distribution after mining in Zhuanlongwan mining area based on visual MODFLOW[J]. Environmental Earth Sciences,2018,77:400. doi: 10.1007/s12665-018-7575-3
[9] 刘中培,李 鑫,陈 莹,等. 基于MODFLOW与气候模式的矿区地下水流模拟[J]. 中国农村水利水电,2021,5:113−117,124. doi: 10.3969/j.issn.1007-2284.2021.05.020 LIU Zhongpei,LI Xin,CHEN Ying, et al. Groundwater flow simulation in mining area based on MODFLOW and climate model[J]. China Rural Water Resources and Hydropower,2021,5:113−117,124. doi: 10.3969/j.issn.1007-2284.2021.05.020
[10] DONG Donglin,SUN Wenjie,XI Sha. Optimización de la Capacidad de Drenaje de una Mina Usando FEFLOW para la Veta de Carbón N° 14 de la Mina Linnancang en China[J]. Mine Water and the Environment,2012,31(4):353−360. doi: 10.1007/s10230-012-0205-5
[11] HUANG Dan,LIU Zaibin,WANG Wenke. Evaluating the impaction of coal mining on ordovician karst water through statistical methods. Water,2018,10(10):1409.
[12] FETHI Lachaal,AMMAR Mlayah,MOURAD Bédir, et al. Implementation of a 3–D groundwater flow model in a semi-arid region using MODFLOW and GIS tools:The Zéramdine–Béni Hassen Miocene aquifer system(east-central Tunisia)[J]. Computers and Geosciences,2012,48:187−198. doi: 10.1016/j.cageo.2012.05.007
[13] 卢帆远. 露天矿采场底板含水层的水位预测及治理优化方法研究[D]. 徐州:中国矿业大学,2022: 41–64. LU Fanyuan. Research on water level prediction and control optimization method of open pit stope floor aquifer[D] Xuzhou:China University of Mining and Technology,2022: 41–64.
[14] 刘 瑜. 陕北侏罗系煤层开采导水裂缝带动态演化规律研究及应用[D]. 徐州:中国矿业大学,2018: 134–162. LIU Yu. Dynamic Evolution and Application of Water Conducting Fractured Zone during Extraction of Jurassic Coal Seams in Northern Shaanxi [D]. Xuzhou:China University of Mining and Technology,2018: 134–162.
[15] ZHANG Wenrui,LIU Tingxi,DUAN Limin, et al. Forecasting groundwater level of karst aquifer in a large mining area using partial mutual information and NARX hybrid model[J]. Environmental Research,2022,213:113747. doi: 10.1016/j.envres.2022.113747
[16] QU Shen,WANG Guangcai,SHI Zheming,et al. Temporal changes of hydraulic properties of overburden aquifer induced by longwall mining in Ningtiaota coalfield,northwest China[J]. Journal of Hydrology,2020,582.
[17] ZHAO Dekang,WU Qiang,CUI Fangpeng,et al. Using random forest for the risk assessment of coal- floor water inrush in Panjiayao COAL Mine,northern China. Hydrogeology Journal,2018,26(7):2327–2340.
[18] KARAMAN Abdullah. Type-curve analysis of water-level changes induced by a longwall mine[J]. Journal of Environmental Geology,2001,40(7):897–901.
[19] KARAMAN Abdullah,SEYIDALI S Akhiev,PHILIP J Carpenter. A new method of analysis of water-level response to a moving boundary of a longwall mine[J]. Water Resources Research,1999,35(4):1001–1010.
[20] 杨 倩. 采动影响下煤层顶板承压含水层水位变化及覆岩变形与破坏规律研究[D]. 阜新: 辽宁工程技术大学, 2013: 23–45. YANG Qian. Study on water level change in confined aquifer and deformationand failure lawofoverlying strata caused by coal mining [D]. Fuxin: Liaoning Technical University, 2013: 23–45.
[21] 李 涛,李文平,孙亚军. 半干旱矿区近浅埋煤层开采潜水位恢复预测[J]. 中国矿业大学学报,2011,40(6):894–900. LI Tao,LI Wenping,SUN Yajun,Prediction of groundwater table recovery affected by approximate shallowburied coal mining in semiarid region [J]. Journal of China University of Mining and Technology,2011,40(6):894–900.
[22] 赵德深,徐孟林,夏洪春. 基于熵权–层次分析的导水裂缝带高度预测研究[J]. 采矿与岩层控制工程学报,2013,18(110):8−10. ZHAO Deshen,XU Menglin,XIA Hongchun. Prediction of water-induced fissure zone heoght based on entropy weight-AHP[J]. Coal Mining Technology,2013,18(110):8−10.
[23] 胡小娟,李文平,曹丁涛,等. 综采导水裂隙带多因素影响指标研究与高度预计[J]. 煤炭学报,2012,37(4):613−620. doi: 10.13225/j.cnki.jccs.2012.04.026 HU Xiaojuan,LI Wenping,CAO Dingtao, et al. Index of multiple factors and expected height of fully mechanized water flowing fractured zone[J]. Journal of China Coal Society,2012,37(4):613−620. doi: 10.13225/j.cnki.jccs.2012.04.026
[24] 国家矿山安全监察局. 建筑物、水体、铁路及主要井巷煤柱留设与压煤开采规范[M]. 北京:煤炭工业出版社,2017. [25] 许家林,朱卫兵,王晓振. 基于关键层位置的导水裂隙带高度预计方法[J]. 煤炭学报,2012,37(5):762−769. doi: 10.13225/j.cnki.jccs.2012.05.002 XU Jialin,ZHU Weibing,WANG Xiaozhen. New method to predict the height of fractured water-conducting zone by location of key strata[J]. Journal of China Coal Society,2012,37(5):762−769. doi: 10.13225/j.cnki.jccs.2012.05.002
[26] 刘士亮. 陕北侏罗系煤田开采环境工程地质模式研究[D]. 徐州:中国矿业大学,2019: 26–90. LIU Shiliang. Environmental engineering geological patterns relating to mining activities in jurassic coalfield of northern Shaanxi [D]. Xuzhou:China University of Mining and Technology,2019: 26–90.
[27] 黄庆享. 榆神府矿区覆岩隔水岩组的隔水性与保水开采研究[C]//中国煤炭学会开采专业委员会. 煤炭开采新理论与新技术——中国煤炭学会开采专业委员会2010年学术年会论文集. 中国矿业大学出版社,2010:11–18. HUANG Qingxiang. Research on water-retaining and water-retaining mining of water-retaining rock formations in overlying rocks in Yushenfu Mining Area [C]// Mining Committee of China COAL Society. New Theory and Technology of Coal mining –- Proceedings of the,2010 Academic Annual Meeting of the Mining Committee of China COAL Society. China University of Mining and Technology Press,2010:11–18.
[28] 李 涛,王苏健,韩 磊,等. 生态脆弱矿区松散含水层下采煤保护土层合理厚度[J]. 煤炭学报,2017,42(1):98–105. LI Tao,WANG Sujian,HAN Lei,Reasonable thickness of protected loess under loose aquifer in ecologicallfragile mining area [J]. Journal of China Coal Society,2017,42(1):98–105.
[29] 李 涛. 陕北煤炭大规模开采含隔水层结构变异及水资源动态研究[D]. 徐州:中国矿业大学,2012.117–142. LI Tao. Study on aquifer-aquifuge structure variation and water resources dynamic evolution with large scale mining in northern Shaanxi Province [D]. Xuzhou:China University of Mining and Technology,2012.117–142.
[30] LIU Shiliang,DAI Song,ZHANG Wenhui,et al. Impacts of underground coal mining on phreatic water level variation in arid and semiarid mining areas:a case study from the Yushenfu mining area,China[J]. Environmental Earth Sciences,2022.81(9):81.
-
期刊类型引用(1)
1. 刘士亮，张文慧，王傲，郑雨生，王山林，毛德强，孙中光. 风沙滩地区煤层采动下潜水位恢复时间影响因素敏感性分析. 中国矿业大学学报. 2024(05): 993-1005 . 
百度学术
其他类型引用(4)
计量
- 文章访问数: 61
- HTML全文浏览量: 8
- PDF下载量: 30
- 被引次数: 5
