GEOFLUIDS Geofluids 1468 - 8123 1468 - 8115 Hindawi 10.1155 / 2019/3920414 3920414 研究文章 试验研究基于渗透率演化在煤矿的突水风险和预防回填的讨论 1 https://orcid.org/0000 - 0003 - 1774 - 448 x 鸡胸 1 https://orcid.org/0000 - 0001 - 6831 - 598 x 维情 1 境况不佳的 1 https://orcid.org/0000 - 0001 - 8742 - 0144 2 Lo Russo 斯特凡诺 1 国家重点实验室的煤炭资源与安全开采 矿业学院 中国矿业大学&技术 徐州221116 中国 cumt.edu.cn 2 交通运输工程学院 淮阴理工学院 淮安223003 中国 hyit.edu.cn 2019年 30. 12 2019年 2019年 10 07年 2019年 14 11 2019年 25 11 2019年 30. 12 2019年 2019年 版权©2019李濛et al。 这是一个开放的文章在知识共享归属许可下发布的,它允许无限制的使用,分布和繁殖在任何媒介,提供最初的工作是正确的引用。

煤炭开采引起的骨折经常发生在地质地层直到发生故障时,提供渠道水流。因此,有必要研究岩石的渗透率演化下的负载。钻孔取样是在基岩层进行下一个含水层,和砂岩的渗透率演化标本在不同围压下岩石力学测试中测试实验室。结果表明,增加围压渗透率逐渐降低,而峰值强度增加而增加的压力。发生在最小和最大渗透率砂岩标本受到弹性形变和应变软化阶段,分别。失败,和最大磁导率,这些砂岩标本不同时发生。为了防止流道形成由于岩石裂缝的发展和失败,回填采空区的方法,提出了回填裂缝发展的影响也进行了讨论。

江苏省高校自然科学基础 KJB440004 中国矿业大学和技术 SKLCRSM19X001 国家杰出青年科学基金 51725403 国家博士后项目创新人才 BX20180361
1。介绍

在煤炭开采过程中,地层运动发生在矿井地面压力的影响进一步产生骨折( 1- - - - - - 3]。在这种背景下,输水渠道形成一次骨折相互连接,这样水的含水层流动开采面板中,造成突水( 4, 5]。近年来,有许多突水事故一个含水层下采煤引起的,导致经济损失和威胁着煤矿的安全生产 6, 7]。此外,骨折不断被开发和扩展,而渗透率根据不同采动应力的影响也不尽相同。因此,有必要探索在岩石渗透率的演化变形故障提供实验数据揭示了煤层在地下含水层的突水机理。

研究人员进行了大量的煤和岩石的渗透率实验研究[ 8- - - - - - 12和数值方法 13- - - - - - 15]。杆菌等。 16)提出了一个新的实用技术方法的评价水力传导率和在复杂的裂隙岩体隧道涌水量,用来评估水流入隧道计划沿着新建的高速公路项目。Bukharov和伊 17)提出了一个方法确定骨折在硬岩石的渗透率在大坝渗流问题,和Zhang et al。 18]研究脆性岩石的变形与渗透率之间的关系。Parisa [ 19]估计岩石的渗透率通过使用神经网络模型;王等人。 20.]研究了透气性的特点运用低渗透性岩石三轴试验,时提供指导建设地下储油设施。郭et al。 21)表明,岩石的渗透率与温度和压力呈正相关在研究温度对岩石的渗透性的影响。孟和李 22)研究了高质量的煤的渗透率的变化。Głowacki和Selvadurai 23)测试的渗透性石灰石在不同压力和提出了一个经验空间渗透性演化的模型。研究在渗流过程中岩石渗透率演化可以大大有助于揭示岩石的渗流机理。

此外,煤炭开采造成突水发生,通常是由于地层破碎煤炭开采期间,这将形成地下水的渗流通道。对于这个问题,提出了一些方法来控制在煤矿突水;同时,现场取得了良好的应用效果。李等人。 24)发现了一个突水机制适合特定的地质和水文地质条件,应用灌浆控制措施。为了防止煤矿的突水,吴et al。 25)提出了脆弱性指数方法通过耦合层次分析法和地理信息系统评价的突水风险。此外,煤层下组合结构,特别是在水体可以通过适当的规划优化提取煤的恢复和系统层控制调查( 26]。然而,上述控制方法不能保证煤炭开采过程中地层的完整性,这将破坏生态环境和地下水。与此同时,一些研究可以结合岩石渗流和突水的机理,迄今缺乏关注的在基岩含水层的渗透性演化。

因此,本研究旨在调查岩层的水渗流机理,提出有效的控制方法。钻孔取样在基岩含水层测试进行不同围压下岩样的渗透率演化,揭示在开采煤层突水的原因。通过岩石渗流实验研究,岩石的渗透率演化和破坏过程之间的关系可以获得准确。最后,回填控制方法对预防煤矿water-inrush-induced灾害提出和讨论。

2。实验的细节 2.1。实验材料

砂岩标本被从基岩层回填面板煤矿在安徽省,中国。的煤矿矿井含水层明显影响,因为第四纪松散地层的底部含水层直接覆盖基岩层,下面没有隔水层。此外,煤层之间的最小距离回填面板和含水层为21.8 m。煤层是含水层的影响,不能发掘传统的采矿方法。含水层的厚度范围从6.12到39.19米的平均20.7米,和它包含的材料包括砂砾、粗燕麦粉、粘土质砾石,粗砂,中沙和粘土砂。眼前的屋顶是粉砂岩和细砂岩组成的高相关的湿和干燥的优势和平均厚度5.8米。核心进行了钻井回填面板测试基岩的渗流特点,提供实验数据,旨在防止未来突水含水层下采煤引起的风险。砂岩的矿物主要有长石、石英、高岭石、伊利石和绿泥石。它平均密度为2570公斤/米3。圆柱形标本测量直径50毫米,和100毫米高度,从大部分标本。

利用扫描电子显微镜(SEM),砂岩的微观结构260 x标本后,1400 x 2400 x放大了,如图 1。大多数矿物粒子表现为不规则的块,它们是相互联系的;然而,有矿物粒子微观缺陷包括微孔隙和微裂隙。微孔隙和微裂隙表现出较低的抗拉和抗压强度与完整岩石相比,他们被认为是形成原始不变的流体渗流通道,允许标本( 27]。

扫描电镜获得的微观结构的砂岩。(一)260 x;(b) 1400 x;(c) 2400 x。

2.2。实验装置和测试方法

的渗流特征砂岩标本通过使用计算机控制的MTS815.03测试系统进行测试。最大轴向载荷是4600 kN。如图 2,上下两端的每个标本受限多孔钢板,以确保统一的水压力施加在整个试样横截面和保证均匀液体渗透到标本的表面。上部和下部的水压力分别对上、下多孔板的标本。垂直孔分别是在多孔板的中心,这样水形成的流道。塑料绝缘套管、热收缩塑料套管,用于密封标本。

原理图的测试岩石标本的渗流特征。

显示了渗流实验的基本机制如下:轴向压力<我n李ne-formula> P 1 ,围压<我n李ne-formula> P 2 ,孔隙压力<我n李ne-formula> P 3 实施(保持<我n李ne-formula> P 3 < P 2 ),然后孔隙压力<我n李ne-formula> P 4 在低端的标本,所以渗透压的差异(减少<我n李ne-formula> P = P 3 P 4 )可以建立整个试样的两端。因此,水被迫标本。

研究压力对砂岩的渗流特征的影响在不同埋藏深度,砂岩标本的围压力下的渗流特征4,6,8,分别测定10 MPa。此外,轴向应力对标本,然后,渗透<我n李ne-formula> k 的样品是通过使用下列公式计算 28, 29日]: (1) k = 1 5 n = 1 n 526年 × 10 6 × lg Δ P 1 Δ P , 在哪里<我n李ne-formula> n ,<我n李ne-formula> Δ P 1 ,<我n李ne-formula> Δ P 数据采集是指行号,渗透压力之间的差异<我n李ne-formula> 1 th和<我n李ne-formula> th线,分别。

在进行渗流实验之前,有必要提前标本完全饱和渗流过程是块状如果标本部分饱和;此外,圆柱形砂岩标本密封测试期间,或标本内的水与油混合在三轴细胞本身。

3所示。实验结果分析

轴向压力和渗透率之间的关系与砂岩在不同围压下的轴向应变图所示 3:这是渗透率之间的相似和不同围压下轴向应变。在不同围压下,主微孔隙和微裂隙在砂岩致密化和关闭阶段标本达到弹性极限压力不断增加。此外,轴向应力和应变之间的线性关系:弹性变形之后,新裂缝砂岩样品开发和扩展,发生塑性变形。砂岩试件达到峰值强度后,压力迅速下降,表现出应变软化。此后,砂岩标本进一步扩张,直到达到他们的剩余强度。

之间的关系压力、渗透率和不同围压下轴向应变。(一)4 MPa;(b) 6 MPa;(c) 8 MPa;(d) 10 MPa。

由于主要的微孔隙和微裂隙砂岩标本成为强化和更加紧密的关闭压力增加在弹性变形阶段,渗透迅速下降。在塑性变形阶段,新的骨折快速开发和扩展和渗透砂岩试件达到极限强度后逐渐增加。然后是渗透率砂岩标本受损后达到最大值。此外,新裂缝砂岩标本中不断发展并最终失败发生,它可以不断地提供渠道水流。

3表明,最小和最大渗透率砂岩标本受到时发生弹性变形阶段和应变软化阶段。此外,发生最大渗透率砂岩试件达到峰值强度后,暗示砂岩标本的失败不同时出现在达到最大渗透率。渗透率突然达到最大后砂岩标本完全损坏。在煤矿,一个渠道水流形成由于突变让水流面板,以便煤矿受到water-inrush-induced灾难,威胁着煤矿的安全开采。

最大渗透率之间的关系和砂岩的峰值强度与围压图所示 4。逐渐的最大渗透率随围压增加而降低,峰值强度随着围压上升。当围压从4 MPa提高到10 MPa,标本的最大渗透率下降<我n李ne-formula> 48.7 × 10 7 μ2来<我n李ne-formula> 17.1 × 10 7 μ2在峰值强度从38.9 MPa提高到98.2 MPa。

最大磁导率之间的关系和砂岩的峰值强度和围压。

4所示。讨论

根据压力和渗透率之间的关系,从变形到破坏过程可分为六个阶段(图 5)[ 30., 31日]。六个阶段包括微裂隙的闭合(OA),弹性变形阶段(AB)、裂缝发育和扩展阶段(BC),快速发展的不稳定骨折(CD),应变软化阶段(DE)和残余强度阶段(EF)。

之间的关系压力、渗透率和压力在这些砂岩标本。

OA,主裂缝进一步被压缩,关闭这样的空隙体积试样整体降低,渗透率下降。在弹性变形阶段(AB),应变逐渐随压力增加而降低。此外,有一个近似的线性轴向应力和应变之间的关系和标本的渗透率逐渐增加。在裂缝的发展和扩展(BC),标本是经塑性变形后达到弹性极限和新骨折和关节出现压力增大,导致渗透率迅速增加。之后,在不稳定骨折快速发展阶段(CD)。在这个阶段,骨折和关节标本的快速开发和扩展后的试样达到屈服极限,直到标本完全破坏,通透性显著增加。在应变软化阶段(DE),标本内的骨折和关节逐渐损坏标本后达到峰值强度和应力逐渐降低到一个稳定的价值(剩余强度)增加压力。此外,渗透率迅速升至最大标本后达到峰值强度。在阶段(EF),标本受到应变软化后,有过压力和渗透率变化增加压力。

在煤矿突水图所示 6。煤炭开采后,立即屋顶坍塌,形成一个屈服区上覆地层压力的影响下由于损失降低煤层的支持。相比之下,主屋顶逐渐走向的下部煤层裂缝逐渐扩展和内部应力达到极限强度时破裂。此外,沿着上覆地层破裂发展形成破碎区。屈服和断裂的区域形成了输水断裂区。一旦输水含水层裂隙带发展,流入的水通过连接板骨折引起的煤炭开采,造成上述突水。

在煤矿突水含水层诱导。

从实验结果可以看出部分 3的上覆地层逐渐沉没在煤层开采后,骨折在地层的影响下逐步发展并最终破裂压力引起的煤矿。随着地层在每一层被破坏,断裂形成的网格,因此多孔。含水层的水可以流到矿业委员会通过骨折和突水随之而来。

通过以上分析,可以推测的是控制地层的裂缝网格是必要的,以免形成一个连接含水层,并防止water-inrush-induced灾害。因此,采空区的面板可以通过使用回填材料回填煤层开采后,由于煤炭开采空间可以回填,回填材料(作为主要支持)被用来承担大量的表层,有效地限制了沉降的屋顶( 32- - - - - - 34]。因此,上覆地层主要是表示为弯曲下沉和骨折在有限的领域,它有效地降低输水裂隙带的高度,保证下煤层的安全开采含水层( 35]。

已经提出,固体废物(煤矸石、粉煤灰和矿渣)回填到采空区控制地层移动和裂缝发展( 36, 37];固体充填开采的方法示意图如图 7。调查勘探的发展在上覆地层断裂区固体充填开采,开采高度和压缩比的关系输水裂隙带的高度是通过回归使用SPSS统计分析软件 35, 38]: (2) H = 31.96 + 2.72 34.56 φ , 在哪里<我n李ne-formula> H 指输水裂隙带的高度,<我n李ne-formula> φ 是压缩比,<我n李ne-formula> 是矿业的高度。

该方法中固体充填开采煤矿的示意图。

而压缩率代表的比例的最终沉降矿业的屋顶高度,压缩比越大,回填效果就越好。 (3) φ = e , 在哪里<我n李ne-formula> e 指的是屋顶的最终沉降。

根据公式( 2),挖掘高度和压缩比的关系的高度的输水派生断裂区,如图 8

输水裂隙带的高度的关系挖掘高度和压缩比。

如图 8开采高度越低,压缩比越大,降低输水裂隙带的高度。因此,对于那些煤矿地质条件特殊,压缩比是一个关键因素在控制输水裂隙带的高度。此外,压缩比越大,越好回填效果和上覆地层的沉降越少,因此,输水区骨折不到达含水层,从而减少突水的风险。

为了监控输水裂隙带的高度在回填采矿面板中,两个监测钻孔回填上面钻孔和安装面板。钻井液是用来监测地层的裂缝发展。两个钻孔的深度1和钻孔2是295.0米,和最终的洞是在地板上煤层的层。挖掘回填面板后,地层位移和变形相对温和,和骨折的发展高度约为10.0米。同时,输水裂隙带的高度没有发展到含水层。监测结果表明,采用回填方法可以有效地减少开发的高度输水区骨折在地层下实现安全开采的含水层。

5。结论

本文旨在研究岩石的渗透率变化在加载和揭示的机制在煤矿突水含水层。为了实现这一目标,在基岩层钻孔取样进行了含水层的渗透性演化下砂岩标本在不同围压下了。

由于主要的微孔隙和微裂隙砂岩标本进一步被强化和关闭随着压力的增加在弹性变形阶段,渗透迅速下降。在塑性变形阶段,新的骨折快速开发和扩展:渗透率逐渐增加砂岩试件达到极限强度后,然后,渗透率突然达到最大值后砂岩标本完全损坏。此外,最小和最大渗透率砂岩标本时发生弹性形变和应变软化阶段。此外,最大渗透率砂岩试件达到峰值强度后出现,这意味着失败的砂岩标本没有发生的同时最大渗透率。

渗透率突然达到最大值后砂岩标本完全损坏。在煤矿,一个渠道水流形成由于突变使得水流面板和造成突水。回填的方法可以用来减少煤炭开采引起的裂缝发展和防止输水骨折区到达含水层,从而减少突水的风险。

数据可用性

使用的数据来支持本研究的结果包括在本文中。

的利益冲突

作者宣称没有利益冲突。

确认

金融支持这项工作由全国博士后项目提供创新人才(BX20180361)、国家杰出青年科学基金(51725403),国家重点实验室的独立研究项目的煤炭资源安全开采,CUMT (SKLCRSM19X001)和江苏省高校自然科学基金(KJB440004)。

J。 w·P。 数值模拟的气体在采空区岩层采动裂隙网络迁移到 国际矿业科技杂志》上 2017年 27 4 681年 685年 10.1016 / j.ijmst.2017.05.015 2 - s2.0 - 85019932715 Palchik V。 形成由于长壁开采上覆岩层的断裂区 环境地质 2003年 44 1 28 38 10.1007 / s00254 - 002 - 0732 - 7 J。 Lv M。 Q。 Y。 K。 煤微晶结构变化与甲烷吸附/解吸 燃料 2019年 239年 13 23 10.1016 / j.fuel.2018.10.155 2 - s2.0 - 85056197099 D。 H。 Y。 力学行为的煤层渗透岩溶陷落柱:采动地下水侵入的风险 自然灾害 2015年 75年 3 2137年 2151年 10.1007 / s11069 - 014 - 1416 - 9 2 - s2.0 - 84920082107 年代。 W。 Y。 太阳 W。 D。 实验模拟煤矿断层突水通道进化的地板上 矿井水和环境 2017年 36 3 443年 451年 10.1007 / s10230 - 017 - 0433 - 9 2 - s2.0 - 85013996544 j . C。 b . H。 中国煤炭开采含水层下:一个案例研究 国际岩石力学和采矿科学杂志》上 2004年 41 4 629年 639年 10.1016 / j.ijrmms.2003.01.005 2 - s2.0 - 1842736088 J。 年代。 R。 Y。 矿山突水灾害场景分析使用贝叶斯网络 安全科学 2016年 89年 231年 239年 10.1016 / j.ssci.2016.06.013 2 - s2.0 - 84976894786 H。 t Y。 解释rock-electric使用孔隙网络模型和渗流特征 石油科学与工程》杂志上 2019年 180年 1 10 10.1016 / j.petrol.2019.05.005 2 - s2.0 - 85065618443 Vinsot 一个。 Linard Y。 Lundy M。 Necib 年代。 Wechner 年代。 见解稀释过程的化学渗透水从水井callovo-oxfordian泥质岩石 Procedia地球和行星科学 2013年 7 871年 874年 10.1016 / j.proeps.2013.03.093 J。 J。 X。 G。 D。 实验研究的装载煤炭在真三轴条件下的渗流特征 岩石力学和岩石工程 2019年 52 8 2815年 2833年 10.1007 / s00603 - 018 - 1720 - x 2 - s2.0 - 85059592124 K。 x B。 z W。 a . H。 渗流的数值模拟结构面及其可行性 中南大学学报 2013年 20. 5 1326年 1331年 10.1007 / s11771 - 013 - 1619 - 3 2 - s2.0 - 84888415215 Z。 J。 Q。 B。 M。 妞妞 Q。 各向异性特征,煤级煤骨折富康矿区,中国 燃料 2018年 211年 182年 193年 10.1016 / j.fuel.2017.09.067 2 - s2.0 - 85029584747 Sweetenham m·G。 麦克斯韦 r·M。 桑蒂 p . M。 评估的时间和震级precipitation-induced渗透到无聊的隧道,通过断裂的岩石 隧道与地下空间技术 2017年 65年 62年 75年 10.1016 / j.tust.2017.02.003 2 - s2.0 - 85014182955 Z。 J。 Q。 妞妞 Q。 J。 H。 X。 各向异性渗透率的变化煤煤在富康矿区不同有效应力下,中国 燃料 2018年 234年 1481年 1497年 10.1016 / j.fuel.2018.08.013 2 - s2.0 - 85051546308 X X。 d Q。 岩体非饱和渗流分析的数值程序包含裂缝网络和排水洞 《水文 2019年 574年 23 34 10.1016 / j.jhydrol.2019.04.014 2 - s2.0 - 85064223338 杆菌 N。 Pranzini G。 Alfi 一个。 Boerio V。 评价岩体渗透率张量和预测隧道流入通过geostructural调查和有限元渗流分析 工程地质 2008年 101年 3 - 4 174年 184年 10.1016 / j.enggeo.2008.05.002 2 - s2.0 - 53549120204 Bukharov r·K。 伊尔 F。 测定渗透率在坚硬的岩石裂缝 水工构造物 1991年 25 8 500年 502年 10.1007 / BF01424131 2 - s2.0 - 34249925289 R。 Z。 太阳 Q。 年代。 变形机制和渗透率之间的关系对脆性岩石 自然灾害 2013年 66年 2 1179年 1187年 10.1007 / s11069 - 012 - 0543 - 4 2 - s2.0 - 84874213646 Parisa B。 委员会为岩石渗透率预测神经网络模型 应用地球物理学杂志 2014年 104年 5 142年 148年 H。 W。 J。 Skoczylas F。 气体渗透性能的低渗透性岩石三轴压缩试验的过程 材料的信件 2014年 116年 2 386年 388年 10.1016 / j.matlet.2013.11.061 2 - s2.0 - 84890050144 X。 G。 Y。 Y。 T。 研究温度对岩石渗透率敏感性的影响 石油科学与工程》杂志上 2017年 156年 616年 622年 10.1016 / j.petrol.2017.06.045 2 - s2.0 - 85021047793 Z。 G。 实验研究级煤的渗透率在不同压力及其影响因素 工程地质 2013年 162年 14 108年 117年 10.1016 / j.enggeo.2013.04.013 2 - s2.0 - 84879503336 Głowacki 一个。 Selvadurai 答:p S。 压力引起的渗透率变化在印第安纳州石灰岩 工程地质 2016年 215年 122年 130年 10.1016 / j.enggeo.2016.10.015 2 - s2.0 - 84998881003 B。 h . B。 J·J。 突水机理由灌浆和控制Chensilou我衡量一个案例研究,中国 阿拉伯地球科学杂志》 2017年 10 1 10 Q。 Y。 D。 W。 底板突水预测:基于gis技术的应用ahp脆弱指数方法Donghuantuo煤矿,中国 岩石力学和岩石工程 2011年 44 5 591年 600年 10.1007 / s00603 - 011 - 0146 - 5 2 - s2.0 - 80053912634 Gandhe 一个。 Venkateswarlu V。 古普塔 r . N。 原位覆盖调查失败的特厚煤层长壁放顶煤开采过程 岩石力学和岩石工程 2005年 78年 5 155年 162年 美国问。 y . H。 Y Y。 W。 问:L。 实验研究砂岩渗流行为的物质有不同的气体压力 Mechanica学报 2015年 31日 6 837年 844年 10.1007 / s10409 - 015 - 0432 - 7 2 - s2.0 - 84951287087 石屏 l Yushou l l Zhenye W。 帮派 Z。 Permeability-strain方程对应的完整应力——应变路径Yinzhuang砂岩 国际岩石力学和采矿科学和地质力学学报文摘 1994年 31日 4 383年 391年 10.1016 / 0148 - 9062 (94)90906 - 7 2 - s2.0 - 0028576135 年代。 Z。 H。 C。 J。 试验研究不同围压下砂岩的孔隙率和渗透率 中国岩石力学与工程学报 2003年 22 5 742年 746年 z . P。 b . Y。 l·W。 z Y。 G。 b . T。 煤矿deformation-failure和渗透率的影响煤层地板上 煤炭地质与勘查 2012年 40 2 39 43 z . P。 J。 x C。 y D。 C。 在煤矿采空区岩体渗透率的计算模型及其数值模拟分析 中国煤炭学会杂志》上 2016年 41 8 1997年 2005年 x J。 g . L。 Q。 耦合的离散element-finite差分法分析沉降控制完全机械化固体回填采矿 环境地球科学 2016年 75年 8 10.1007 / s12665 - 016 - 5517 - 5 2 - s2.0 - 84963649308 破车 M。 Witthaus H。 的研究和应用进展与充填采煤 国际矿业科技杂志》上 2013年 23 1 7 12 10.1016 / j.ijmst.2013.01.002 2 - s2.0 - 84876953792 M。 J。 首歌 W。 日尔曼 d . M。 回收的碎矸石回填材料在煤矿:颗粒大小对压实的影响行为 环境科学与污染研究 2019年 26 9 8789年 8797年 10.1007 / s11356 - 019 - 04379 - 9 2 - s2.0 - 85061033684 30712211 M。 j . X。 x J。 N。 Q。 水保护基于固体充填采矿法在水下轴承地层及其应用 中国煤炭学会杂志》上 2017年 42 1 127年 133年 B。 F。 一个实验调查粒状煤矸石回填材料的压实特性修改绑定 环境地球科学 2018年 77年 7 10.1007 / s12665 - 018 - 7474 - 7 2 - s2.0 - 85045002309 J。 Y。 Z。 J。 H。 Y。 特征的宏观变形和粒子破碎粉碎煤矸石颗粒材料在循环荷载:在固体回填煤矿 粉技术 2019年 343年 159年 169年 10.1016 / j.powtec.2018.11.049 2 - s2.0 - 85056624505 J。 H。 X。 F。 预测勘探区域的高度高于固体充填开采的开采面板 矿井水和环境 2014年 33 4 317年 326年 10.1007 / s10230 - 014 - 0310 - 8 2 - s2.0 - 84911806644