移动阅读
马雄德1,2，黄金廷3，李吉祥4，宁世雄4
(1.长安大学 环境科学与工程学院，陕西 西安 710054; 2.长安大学 旱区地下水文与生态效应教育部重点实验室，陕西 西安 710054; 3.中国地质调查局西安地质调查中心,陕西 西安 710054; 4.西安科技大学 地质与环境学院，陕西 西安 710054)
摘 要:分布在干旱半干旱地区的植被，由于降水不足以维持其长期生存，需要地下水提供部分或全部水源，因而对地下水有一定的依赖性。煤层开采破坏含水层后地下水位会大幅降落，这在一定程度上会给依赖地下水植被造成水分胁迫，进而控制生态系统演化过程。针对榆神矿区煤层开采引起地下水位变化的基本特征，提出了生态安全约束下矿区地下水位控制阈值的确定方法。研究表明，植物根系与地下水毛细上升带保持接触时，植物就可以吸收利用地下水，因此本文将最大毛细上升高度与根系长度之和作为植被利用地下水的最大临界埋深。在毛管流理论指导下，以颗粒排列方式与孔隙直径大小的关系建立了最大毛细上升高度计算公式，并给出了通过颗粒级配曲线确定最大毛细上升高度的方法。据此计算的毛乌素风积沙最大毛细上升高度的取值区间为0.7～2.0 m，进一步确定了榆神矿区生态安全约束下的矿区水位控制下限为4.0 m。在此基础上，以2016年地下水流场基准，以水位埋深4.0 m为界将榆神矿区划分为生态约束区和无约束区。水位埋深小于4.0 m的区域植被对地下水依赖程度高，属于生态约束区，煤层开采造成地下水位下降极易使植被遭受水分胁迫，因而是矿区生态环境保护的重点。研究成果阐明了榆神矿区生态环境及地下水位对煤层开发的限制条件，为进一步推进保水采煤技术的发展奠定了理论基础。
关键词:榆神矿区;面向生态;最大毛细高度;水位阈值;保水采煤
移动阅读
马雄德，黄金廷，李吉祥，等.面向生态的矿区地下水位阈限研究[J].煤炭学报,2019,44(3):675-680.doi:10.13225/j.cnki.jccs.2018.6050
MA Xiongde，HUANG Jinting，LI Jixiang，et al.Groundwater level threshold under the constrain of ecology security in mining area[J].Journal of China Coal Society,2019,44(3):675-680.doi:10.13225/j.cnki.jccs.2018.6050
中图分类号:TD823;TD88
文献标志码:A
文章编号:0253-9993(2019)03-0675-06
收稿日期:2018
1225
修回日期:20190220
责任编辑:韩晋平
基金项目:国家自然科学基金资助项目(41672250)
作者简介:马雄德(1978—)，男，青海互助人，高级工程师，博士。E-mail:hgmxd@chd.edu.com
通讯作者:黄金廷(1979—)，男，辽宁建昌人，高级工程师，博士。E-mail:59280431@qq.com
MA Xiongde1,2，HUANG Jinting3，LI Jixiang4，NING Shixiong4
(1.College of Environmental Science and Engineering,Chang’an University,Xi’an 710054,China; 2.Key Laboratory of Subsurface Hydrology and Ecological Effects in Arid Region,Ministry of Education,Chang’an University,Xi’an 710054,China; 3.Xi’an Center of Geological Survery,CGS,Xi’an 710054,China; 4.College of Geology and Environment,Xi’an University of Science and Technology,Xi’an 710054，China)
Abstract:Vegetation distributed in arid and semi-arid regions,due to insufficient precipitation to maintain long-term survival,requires groundwater to provide part or all of its water demand,and thus is phreatophyte.When the underground mining destroys the aquifer,the groundwater level will drop significantly,which will cause water stress to phreatophyte,and control the ecosystem evolution process.In this paper,the authors investigated the basic characteristics of groundwater level change caused by coal seam mining in Yushen mining area,and proposed a method to determine groundwater level threshold under ecological constraints in mining area.As long as the roots are in contact with the groundwater capillary rise zone,the groundwater can be utilized by plants.Therefore,the sum of the maximum capillary rise height and the root length is regarded as the maximum critical depth of the groundwater.Under the guidance of capillary flow theory,the formula to calculate maximum capillary rise height is established by the relationship between particle arrangement and pore diameter,and the method of determining the maximum capillary rise height by particle gradation curve is given.According to the calculation,the maximum capillary rise height of the Mu Us aeolian sand is 0.7-2.0 m and the lower limit of the water level under the ecological constraints in the mine area is 4.0 m.Based on the groundwater flow field in 2016,the Yushen mining area is divided into ecologically restricted areas and unconstrained areas with a water depth of 4.0 m.Vegetation with a water depth of less than 4.0 m is highly dependent on groundwater and belongs to the ecological confinement zone.The decline of groundwater level caused by coal seam mining is very likely to cause water stress on vegetation,which is the focus of ecological environment protection in mining areas.This illuminates the limitations of ecology and water level on coal mining in Yushen mining area,which lays a theoretical foundation for further promoting water-preserved coal mining technology.
Key words:Yushen coal mine area;ecology orientated;maximum capillary height;groundwater level threshold;water-preserved coal mining
长期以来，人们非常关注湿地或河岸带地下水与植被生长之间的关系，而忽略了在干旱环境中也存在植被依赖地下水的情况。在地下水浅埋区，地下水为植被生态系统提供生长所需的部分或全部水源，物种组成和自然生态过程主要由地下水控制，因而被称作依赖地下水的生态系统(GDEs)[1]。
榆神矿区地处毛乌素沙漠南缘，属于我国西北干旱半干旱地区，年降水量不足400 mm，仅能维持耐旱植物长期生存，植被稀疏，种类单调，生物量低。尽管降水量稀少，但广泛分布的风积沙十分有利于大气降水入渗，并在地形及水文地质条件的控制下，缓慢向毛乌素沙地东南缘径流，使沙漠滩地区地下水埋藏较浅，常分布一些河岸带植被和湿地植被。由于降水量不能满足这些植物长期生长所需水分，因而需要吸收利用地下水[2-3],也属于GDEs。在长期自然选择中，为了能够吸收各个深度上的水分，这类植被往往会形成两个根系密度较高的分布区。一个处于浅层土壤中，可以吸收储存在土壤中的大气降水入渗的水分，另一个位于地下水最大毛细上升高度区附近，便于吸收利用地下水[4]。
研究发现，在潜水位处于一定范围时，地下水会通过毛细作用改变根系层土壤的含水率，进而影响植被长势、作物产量等。VIDELA MENSEGUE等[5](2015)发现，地下水位在1.5～2.0 m波动时，地下水对大豆生长需水的贡献量大约为30%。BABAJIMOPOULOS等[6](2007)发现在野外条件下，玉米根系吸收的地下水达3.6 mm/d，占总蒸腾量的18%。HURST 等[7](2004)分析了甘蔗作物吸收利用地下水的特征，发现地下水位大于1 m时就需要灌溉来补充作物需水，灌溉量取决于作物根系(长度和密度)，这说明浅层地下水的毛细上升为生长季植被提供了所需水源[8]。因而地下水位降低时，势必会抑制毛细上升对根系层土壤含水率的调节作用，对植被或作物产生不利影响[9]。
在榆神矿区东部开采区由于煤层较厚，埋藏浅，煤层开采过程中出现了地下水位下降的普遍问题[10-11]。地下水位下降后毛细上升不能到达植被根系层，从而引起了植被退化[12-13]。范立民等[10](2016)研究认为，1995—2014年，榆神府矿区有758.9 km2的区域地下水位下降超过8 m。这些矿区的煤层埋深普遍小于100 m，煤层开采后导水裂隙带直接导通含水层，使地下水位发生剧烈变化。对于湿生植被和中生植被而言，开采造成的水位下降极大地减弱了地下水对植被的涵养作用[14]，使生态安全遭受威胁，因而亟需建立矿区地下水位控制阈限。
因此我们定义面向生态的地下水位阈限的含义为,在地下水与植被关系密切的区域，煤层开采引起地下水位下降，其埋深不宜超过地下水的最大毛细高度与植被根系长度之和，否则植被会遭受干旱胁迫。笔者建立了采用土壤颗粒级配曲线确定最大毛细高度模型，结合榆神矿区沙柳根系调查结果，计算确定了榆神矿区地下水位控制下限，为保水采煤技术提供了理论依据。
非饱和土中的毛细上升现象是指地下水向水位以上某一高度运动的物理现象，其作用力主要来自土粒之间气-水交界面上表面张力产生的负孔隙水压力。地下水的毛细上升高度估算一般是将土壤颗粒等效成一系列毛管束，采用Young-Laplace方程进行计算。对于等效直径为d的毛细管，假定毛细水上升高度为hc，稳定之后，水柱产生的重力与作用在弯液面周围的表面张力Ts的垂直投影建立垂向力平衡方程，可得到
(1)
式中，ρw为水的密度;g为重力加速度;α为接触角。
于是毛细管中水的上升高度(即毛细上升高度)可表示为
(2)
在常温下(20 ℃)，水的密度ρw=1 g/cm3，重力加速度g=980 cm/s2，表面张力Ts=72.75 mN/m[14]，假定水充分浸湿砂粒表面，接触角α=0。代入式(2)，可得到毛细上升高度与等效毛细管直径的关系:
(3)
这说明充分静止的非饱和土壤中毛细上升高度与等效毛细管直径成反比关系。
一般情况下，土是由一系列粒径不同的土颗粒，按照一定的级配和复杂的排列方式组成的，所以精确计算毛细管的直径极其困难。为了获得简单易用的毛细上升高度计算解析公式，可以假定土壤骨架由一系列等径球状颗粒组成，则有2种理想的土颗粒排列形式，即立方体排列和四面体排列(图1)。相对而言，立方体形式排列的土颗粒最为松散，具有较大的孔隙直径(dsc)，而四面体形式排列的土颗粒最为紧密，具有较小的孔隙直径(dth)，由此可以分别估算出土粒不同组合方式下毛细上升高度的取值区间。
图1 颗粒排列方式与孔隙直径大小关系
Fig.1 Relationship curve between particle arrangement and pore diameter
图1(a)所示的立方体形式排列的土颗粒直径与最小孔隙直径的关系为
(4)
图1(b)所示的四面体形式排列的土颗粒直径与最小孔隙直径的关系为
(5)
在自然界，由颗粒组成的骨架不可能有这样规则的排列形式，式(4)和式(5)大致上构成了孔隙直径的一个参考区间。将式(4),(5)分别代入式(3)，则可得到理想化土体内毛细上升高度的阈值区间:
(6)
当然，土是由一定级配的土粒构成，不存在完全等粒的土体骨架，但可以将不同粒径构成的连续的毛细水带概化为不同直径的管束。这样，在平衡状态下，不同直径的毛细管具有不同的上升高度，直径越小的毛细管中水上升的高度越大。PECK(1974)[15]认为毛细上升高度与粒径累计曲线上累计含量为10%对应的粒径(D10)成反比。LANE和WASHBURN(1964)[16]通过实测8种不同的土毛细上升高度，建立了毛细上升高度与D10之间的线性关系。这充分说明，10%的细粒直径足够描述最小连续毛细管的有效直径[17]，因而能代表该土壤毛细上升高度的最大值。因此将式(6)转化为
(7)
由式(7)可知，在确定了土中累计含量为10%的细粒直径后，便可获得在分选程度、颗粒排列、形状等因素影响下的毛细上升高度的取值范围。再加上根系分布长度，即可确定基于生态安全的水位下限值。
粒径累计曲线的形状决定了粒径分布，也就是级配特性。笔者在榆神矿区选取了11个点采样，采集了0～200 cm处土壤样品，在实验室进行颗粒分析。试样的粗细程度和颗粒级配用筛分方法测定，用级配区表示砂的级配。筛分是用一套孔径为0.5，0.25，0.075，0.005 mm的标准筛，将200 g干砂由粗到细依次过筛，称量各筛上的筛余量(g)，计算各筛上的分计筛余率(%)，再计算累计筛余率(%)。为了便于分析，颗粒级配曲线横坐标的筛孔尺寸采用对数，纵坐标的通过百分率为常坐标。11组实验结果为:粒径主要分布在0.075～0.250 mm，含量高达86.7%以上，粒径分布较为集中。颗粒级配较为相近，差异性不大，级配曲线较陡，如图2所示。由图2可知，10%的细粒直径(D10)为0.01 cm。
图2 风积沙级配曲线
Fig.2 Grading curve for aeolian sand
在榆神矿区选择5个地方开挖了树龄在3～5 a的沙柳根系。沙柳根系长度最长1.3～2.4 m，平均2.2 m。在垂向上，沙柳根系密度出现了两个高值区，第1个高值区位于地表下0～35 cm，这使沙柳能够最大限度地吸收降水入渗后的土壤水及养分。第2个高值区位于地下水位以上20～40 cm，这允许沙柳通过毛细上升吸收利用地下水。柳琳秀[18](2015)也发现了类似的规律，说明沙柳生长与地下水关系密切[19]。
为了能够吸收不同深度上的水分，旱区的深根植物一般会发展出特殊的根系系统。在特定区域要么根系密度非常大，要么根系水力传导率较高，其作用和不同环境条中水分利用率的提高是分不开的。在实验室条件[20]和野外条件下[21]，地下水位以上根系密度增加这一现象都可以观察到，因此地下水对植被水分的供应是不可忽视的。
采用前文所提出的式(7)，在获得土壤颗分的基础上，就可以确定基于生态的水位控制阈值。由级配曲线确定10%的细粒直径(D10)为0.01 cm，代入式(7)可得毛细上升高度的取值区间为0.7～2.0 m,沙柳根系长度按平均值取2.2 m,于是面向生态保护的矿区水位埋深阈值区间确定为2.9～4.2 m。由此可知，对于沙柳而言，煤层开采过程中水位埋深不能超过4.2 m这个阈限，否则会对植被生长造成不良影响。
采用前述采样地采集的样品进行了毛细上升高度试验。将不同深度上取得的样品晾干后过0.5 mm的筛去除杂质，按取样顺序装入长200 cm内径2 cm的厚壁玻璃管中，并固定在支架上，其下端插入水槽内(定水位)。当玻璃管中毛细水上升速率减小到0.2 cm/10 h(玻璃管上刻度的精确读数为0.2 cm)，认为毛细高度达到了最大值，试验结束。实验结果表明，所测定的毛细上升最大高度范围在66～116 cm，最大毛细高度为116 cm(图2)。测定的毛细高度在式(7)预测取值区间内，但由于毛细上升高度除了与粒径有关外，还与诸多因素相关，如空隙度、结构和水的特征，包括水温，矿化度，黏滞度等，毛细高度最大值并没有达到土壤颗粒按理想四面体形式排列所计算的2.0 m。
图3 毛细上升高度随时间变化曲线
Fig.3 Capillary rise height over time
蒸散发会降低土壤水的势能，使根系吸水变得更加困难。当土壤比较湿润时，土壤水势较高，容易被植物根系吸收。而较为干燥的土壤水势较低，土壤水被束缚在颗粒周围，很难被根系吸收。因而，当土壤水势达到一定阈值时，植物会产生水分胁迫现象。植物遭受水分胁迫的程度采用水分胁迫系数(Ks)来表示：
(8)
式中，ETa为实际腾发量;Kc为作物系数;ETo为潜在腾发量。
采用文献[12]建立的数学模型计算两种情景下实际腾发量(ETa)，即天然水位(变化幅度在1.5～4.0 m)和固定水位(8.0 m)。模型中根系长度确定为2.0 m。潜在腾发量ETo采用修正的Penman-Monteith公式计算[22]。
图4表达了不同水位条件下植物遭受干旱胁迫的程度。从图4可以看出，以天然水位作为计算情景时，除了2007年外，计算的不同年份干旱胁迫系数约等于1，说明在天然条件下植物不遭受干旱胁迫。当地下水位限定在8.0 m后，情况发生了显著的变化。理论上地下水毛细上升最大2.0 m，模型中给定的植物根系长度2.0 m，因此，地下水对于植物需水基本没有贡献，这使干旱胁迫系数降低到0.2左右，植物遭受干旱胁迫，进而会出现病态的生理生化反应。
图4 不同水位下水分胁迫变化规律
Fig.4 Variation of water stress at different water levels
对于依赖地下水的生态系统(GDEs)，植物是否能够吸收利用地下水，除了与地下水位及毛细高度相关外，还与植被根系长度息息相关。物种不同，植物根系发育长度也不同，浅根植物需要更高的潜水位，而深根植物则可以提升更深的地下水。因此在区域上界定一个对所有植被都适用的水位控制阈限是不现实的，但是可以针对制定优势植被制定水位控制阈值。
SCHENK和JACKSON(2002)[23]发现所有物种95%以上的根系都集中在2.0 m以浅的土壤中。ROMERO-SALTOS等(2005)[24]也通过示踪试验验证了埋深2 m以外的根系几乎不吸水。吴宏伟(2017)[25]通过木质部导管直径确定的根系内部毛细上升高度最大值为1.86 m。因此，为了探讨一般规律，将榆神矿区所有植被可能吸收利用地下水的根系长度界定为2.0 m。再加上最大毛细上升高度后，确定植被能够吸收利用地下水的阈值为4.0 m。这与文献[19]采用NDVI指数与地下水位叠加后建立的统计关系相符。
以2016年地下水流场基准，在ArcGIS软件中圈定榆神矿区地下水位埋深小于4.0 m的区域，这是植被与地下水关系密切的区域，也是开采造成地下水位下降后会对植被产生干旱胁迫的区域，是矿区生态环境保护的重点(图5)。
图5 榆神矿区生态约束分区
Fig.5 Distribution and regionalization of ecology restrict on coal mining
生态约束区主要位于秃尾河以西及榆溪河流域，植被与地下水关系密切，是典型的GDEs。在煤层开采过程中，如果地下水位埋深下降超过4 m，植被生长将会遭受干旱胁迫，进而会产生一系列生理响应[26]。
无约束区，地下水对植被生长所需水分的贡献有限，在这些区域开采造成地下水位下降后，植被响应也有限。尽管如此，无约束区煤层开采如果造成了区域地下水位下降[27]，也会袭夺临近流域的地下水，从而波及到生态约束区。
(1)在地下水与植被关系密切的区域，煤层开采引起地下水位下降，其埋深不宜超过地下水的最大毛细高度与植被根系长度之和，否则植被会遭受干旱胁迫。文中建立了采用土壤颗粒级配曲线确定最大毛细高度模型，结合榆神矿区沙柳根系调查结果，计算确定了榆神矿区地下水位控制下限为4.2 m。
(2)由于毛细上升高度除了与粒径有关外，还与诸多因素相关，如空隙度、结构和水的特征，包括水温，矿化度，黏滞度等，实测的毛细高度最大值不能达到土壤颗粒按理想四面体形式排列所计算的2.0 m。但根据前人的研究，将植物根系长度进行概化后确定为2.0 m，由此确定地下水位埋深阈限为4.0 m是相对安全的。
(3)以地下水位埋深4.0 m为界，榆神矿区开采过程中存在生态约束的区域位于秃尾河以西及榆溪河流域。该区域地下水为埋深普遍小于4.0 m，地下水与植被关系密切，在后期开采中要控制地下水位下降幅度。秃尾河左岸的广大区域，依赖地下水的植被对煤层开采的限制较小。
参考文献
[1] BOULTON A J.Chances and challenges in the conservation of groundwaters and their dependent ecosystems[J].Aquatic Conservation Marine & Freshwater Ecosystems,2005,15(4):319-323.
[2] LOHEIDE S P,BUTLER J J,GORELICK S M.Estimation of groundwater consumption by phreatophytes using diurnal water table fluctuations:A saturated-unsaturated flow assessment[J].Water Resources Research,2005,41(7):372-380.
[3] ORELLANA F,VERMA P,II S P L,et al.Monitoring and modeling water-vegetation interactions in groundwater-dependent ecosystems[J].Reviews of Geophysics,2012,50(3):RG3003.
[4] BUTLER J J,KLUITENBERG G J,WHITTEMORE D O,et al.A field investigation of phreatophyte-induced fluctuations in the water table[J].Water Resources Research,2007,43(430):299-309.
[5] MENSEGUE H R V,DEGIOANNI A J,CISNEROS J M.Estimating shallow water table contribution to soybean water use in Argentina[J].European Scientific Journal,2015,11(14):1857-7881.
[6] BABAJIMOPOULOS C,PANORAS A,GEORGOUSSIS H,et al.Contribution to irrigation from shallow water table under field conditions[J].Agricultural Water Management,2007,92(3):205-210.
[7] HURST C A,THORBURN P J.Sugarcane water use from shallow water tables:Implications for improving irrigation water use efficiency[J].Agricultural Water Management,2004,65(1):1-19.
[8] WU Y,LIU T,PAREDES P,et al.Water use by a groundwater dependent maize in a semi-arid region of Inner Mongolia:Evapotranspiration partitioning and capillary rise[J].Agricultural Water Management,2015,152:222-232.
[9] KROES J,SUPIT I,DAM J V,et al.Impact of capillary rise and recirculation on crop yields[J].Hydrology & Earth System Sciences Discussions,2018:1-31.
[10] 范立民,马雄德,冀瑞君.西部生态脆弱矿区保水采煤研究与实践进展[J].煤炭学报,2015,40(8):1711-1717.
FAN Limin,MA Xiongde,JI Ruijun.The progress of research and engineering practice of water-preserved coal mining in western eco-environment frangible area[J].Journal of China Coal Society,2015,40(8):1711-1717.
[11] 范立民.保水采煤的科学内涵[J].煤炭学报,2017,42(1):27-35.
FAN Limin.Scientific connotation of water-preserved mining[J].Journal of China Coal Society,2017,42(1):27-35.
[12] 马雄德,范立民,严戈,等.植被对矿区地下水位变化响应研究[J].煤炭学报,2017,42(1):44-49.
MA Xiongde,FAN Limin,YAN Ge,et al.Vegetation responses to groundwater level change in mining area[J].Journal of China Coal Society,2017,42(1):44-49.
[13] 卞正富,雷少刚,常鲁群,等.基于遥感影像的荒漠化矿区土壤含水率的影响因素分析[J].煤炭学报,2009,34(4):520-525.
BIAN Zhengfu,LEI Shaogang,CHANG Luqun,et al.Affecting factors analysis of soil moisture for arid mining area based on TM images[J].Journal of China Coal Society,2009,34(4):520-525.
[14] VARGAFTIK N B,VOLKOV B N,VOLJAK L D.International tables of the surface tension of water[J].Journal of Physical & Chemical Reference Data,1983,12(3):817-820.
[15] PECK R B,HANSON W E,THORNBURN T H.Foundation engineering(2nd Edition)[J].Soil Mechanics & Foundation Engineering,1974,32(3):73-75.
[16] LANE K S,WASHBURN D E,KRYNINE D P.Capillarity tests by capillarimeter and by soil filled tubes[J].Highway Research Board Proceedings,1947,26.
[17] LIU Q,YASUFUKU N,MIAO J,et al.An approach for quick estimation of maximum height of capillary rise[J].Soils & Foundations,2014,54(6):1241-1245.
[18] 柳琳秀.毛乌素沙地三种植物根系垂直分布研究[D].包头:内蒙古大学,2015.
LIU Linxiu.Study on root vertical distribution of three kinds of Shrubs in Mu us sandy land[D].Baotou:Inner Mongolia University,2015.
[19] 马雄德,范立民,张晓团,等.基于植被地下水关系的保水采煤研究[J].煤炭学报,2017,42(5):1277-1283.
MA Xiongde,FAN Limin,ZHANG Xiaotuan,et al.Water-preserved mining based on relationship between vegetation and groundwater[J].Journal of China Coal Society,2017,42(5):1277-1283.
[20] JARRELL W M,VIRGINIA R A.Response of mesquite to nitrate and salinity in a simulated phreatic environment:Water use,dry matter and mineral nutrient accumulation[J].Plant & Soil,1990,125(2):185-196.
[21] FALKINER R A,NAMBIAR E K S,POLGLASE P J,et al.Root distribution of Eucalyptus grandis,and Corymbia maculata,in degraded saline soils of south-eastern Australia[J].Agroforestry Systems,2006,67(3):279-291.
[22] EVEN M J.How does groundwater subsidy of vegetation change as a function of landscape position and soil profile characteristics at the Ciha Fen(Johnson County,IA,USA)[D].Iowa:The University of Iowa,2014.
[23] SCHENK H J,JACKSON R B.The Global Biogeography of Roots[J].Ecological Monographs,2002,72(3):311-328.
[24] ROMERO-SALTOS H,STERNBERG L M,NEPSTAD D.Rainfall exclusion in an eastern amazonian forest alters soil water movement and depth of water uptake[J].American Journal of Botany,2005,92(3):443-455.
[25] 吴宏伟.大气-植被-土体相互作用:理论与机理[J].岩土工程学报,2017,39(1):1-47.
WU Hongwei.Atmosphere-plant-soil interactions theories and mechanisms[J].Chinese Journal of Geotechnical Engineering.2017,39(1):1-47
[26] FAN Limin,MA Xiongde.A review on investigation of water-preserved coal mining in western China[J].International Journal of Coal Science & Technology,2018,5(4):411-416.
[27] 冀瑞君,彭苏萍,范立民,等.神府矿区采煤对地下水循环的影响——以窟野河中下游流域为例[J].煤炭学报,2015,40(4):938-943.
JI Ruijun,PENG Suping,FAN Limin,et al.Effect of coal exploitation on groundwater circulation in the Shenfu mine area:An example from middle and lower reaches of the Kuye River basin[J].Journal of China Coal Society,2015,40(4):938-943.