Structure characteristics of microbial community in acid mine drainage in Shandi River Basin
-
摘要:
为了揭示不同氧化还原环境条件下酸性矿井水微生物群落结构特征及其与水环境相互作用关系，选取娘子关泉域内典型废弃煤矿区山底河流域，采集不同氧化还原环境的AMD水样，用于水化学同位素分析和微生物16S rRNA基因V4、V5区测序。水化学同位素分析得出AMD水化学特征主要受pH和围岩水岩作用影响，好氧和厌氧环境下AMD的pH与Eh呈负相关；缺氧厌氧环境下AMD的pH与 Fe3+、Eh呈负相关。AMD中Ca、Mg离子则主要来源于围岩和煤中方解石、白云石的溶解；K、Na主要来源于长石等硅酸岩在酸性条件下的溶解；${\mathrm{SO}}_4^{2-} $、Fe含量主要来源含煤地层中黄铁矿氧化；AMD的硫同位素演化特征为好氧环境下AMD的δ34S–${\mathrm{SO}}_4^{2-} $较低，缺氧厌氧环境下AMD的δ34S–${\mathrm{SO}}_4^{2-} $较高。山底流域AMD微生物多样性影响分析，得出pH越低，样品的细菌丰度和多样性越低；温度越高样品的细菌丰度和多样性越高。通过AMD的微生物群落与环境因子的相关性分析，得出柳沟AMD的优势菌属为Ferrovum、Gallionella和Ferritrophicum。庙沟露天矿AMD的优势菌属为Ferrovum、Gallionella和Acidithiobacillus。榆林垴钻井AMD的优势菌属为Ferrovum和Gallionella。山底河流域AMD的微生物群落受pH、温度和降水的影响较大。
Abstract:In order to reveal the structural characteristics of microbial communities in acid mine drainage(AMD) under different redox environmental conditions and the interaction between microbial communities and groundwater environment, Select the typical abandoned coal mining area in the Shandi River basin of the Niangziguan Spring area, and collect AMD samples from different redox environments for Hydrochemical characteristics and isotope analysis and the high-throughput sequencing of the V4, V5 region of microbial 16S rRNA. The analysis of Hydrochemical characteristics and isotope indicates that the hydrochemical characteristics of AMD are mainly influenced by pH and surrounding rock. The analysis of hydrochemical isotopes shows that the hydrochemical characteristics of AMD are mainly influenced by pH and rock water interaction. Under aerobic and anaerobic environments, the pH of AMD is negatively correlated with Eh; The pH of AMD in anaerobic and hypoxic environments is negatively correlated with Fe3+and Eh.The Ca and Mg ions in AMD mainly come from the dissolution of calcite and dolomite in the surrounding rock and coal. The K and Na ions in AMD mainly come from the dissolution of silicate rocks such as feldspar under acidic conditions; The main source of ${{\mathrm{SO}}_4^{2-}} $ and Fe content in AMD is the oxidation of pyrite in coal-bearing strata; The sulfur isotope evolution of AMD is characterized by lower δ34S–${{\mathrm{SO}}_4^{2-}}$ values in aerobic environments and higher δ34S–${{\mathrm{SO}}_4^{2-}}$ values in anaerobic and hypoxic environments. The analysis of the microbial diversity impact of AMD in the bottom of the mountain basin shows that the lower the pH, the lower the bacterial abundance and diversity of the sample; The higher the temperature, the higher the bacterial abundance and diversity of the sample. Through the correlation analysis between the microbial community and environmental factors of AMD, it was found that the dominant bacterial genera in the Liugou are Ferrovum, Gallionella and Ferritrophicum. The dominant genera of AMD in the Miaogou open-pit mine are Ferrovum, Gallionella and Acidithiobacillus. The dominant genera of AMD in Yulinnao drilling are Ferrovum and Gallionella. The microbial community of AMD in the Shandi River Basin is greatly affected by pH, temperature and precipitation.
-
0. 引 言
煤岩体结构改造是解决煤矿众多技术难题的共性核心科学问题[1]。煤矿中的坚硬厚及特厚煤层综放开采、冲击矿压、低渗透性煤层瓦斯抽采、煤与瓦斯突出等是制约矿井安全高效生产的技术难题，而煤岩体致裂手段是实现煤岩性质弱化、卸压增透[2]的有效方法。依据冲击应力产生的形式与特征，常用的致裂方式有炸药爆破、水压致裂和二氧化碳相变爆破。传统炸药爆破因其施工要求高，炸药材料受到严格管控，可控性、预见性较低等特点，使用的条件越发受到严格限制，如《煤矿安全规程》禁止在采动影响区内使用炸药爆破；水力压裂则存在形成裂缝单一[3]、水资源需求大、保压困难，工作面易污染等问题；二氧化碳相变爆破作为一种非炸药破岩技术，在煤岩体结构改造中起到重要补充作用。
二氧化碳相变爆破是一种非炸药低扰动破岩技术，它将二氧化碳密闭于专用爆破管内，通过高热快速激发，将液态二氧化碳迅速转化为气态，瞬间产生强大的膨胀能量，快速释放的高压气体在岩体内产生拉应力场，从而达到致裂岩体的效果。该技术起源于1914年欧美的CARDOX技术，最开始应用于高瓦斯煤层的开采。在美国、英国和加拿大被广泛应用，曾被美国矿业局誉为当时煤矿领域最著名的技术发明[4-5]。据统计，在20世纪40年代，相变爆破生产了美国25%的井工煤炭产量，年爆破量曾经达到了280万次之多[6]。随着机械化采煤技术及相关综采设备的不断发展，该技术逐渐退出煤矿行业，逐渐转向城区土石爆破和水泥行业清堵等方面，现阶段共计100多个国家应用该项技术[7]。相变爆破技术在20世纪90年代引入我国，21世纪初在我国高瓦斯矿井采煤工作面进行采煤试验，经过相关实践发现了其应用后，巷道内瓦斯浓度增加，从而发现了该项技术应用于矿井瓦斯治理的可行性[8]。2012年来，该项技术被广泛系统地应用于矿井瓦斯综合治理[9]，主要包括卸压掘进、卸压均化瓦斯场、卸压增透抽采瓦斯、顶板安全致裂、井筒揭煤和防治冲击矿压等多个方面[10]，在不同矿区均得到有效应用推广[7, 11]。2014年后从地下矿山拓展到露天矿山，应用于矿山、道路、基坑等露天岩石爆破工程[12-13]。
尽管相变爆破在煤矿中的应用得到不断的发展，但其相关的基础理论尚处于起步阶段，相关的生产规范还没有形成，大规模的推广应用仍旧存在一定困难。笔者依据二氧化碳相变爆破技术特征，确定了相变爆破致裂机理，分析了相变爆破影响因素，探究了相变爆破应用场景，并提出二氧化碳相变爆破致裂技术研究展望，旨在推进相变爆破的进一步推广应用。
1. 二氧化碳相变爆破技术特征
1.1 相变爆破原理
二氧化碳相变爆破是利用二氧化碳相变时产生的物理膨胀特性，将冲击能量直接作用在周围介质中，它主要由排气头、定压卸能片、储液管、发热管和注液头组成（图1）[14]。爆破产生原理如图2所示，将液态二氧化碳充装进管内，利用发热管对二氧化碳进行加热，随着储液管内温度不断增加，二氧化碳发生汽化，体积急剧膨胀，压力也随之增加，当储液管内的压力大于定压卸能片的强度时，定压卸能片破裂，二氧化碳从排气头中喷出，周围介质受到二氧化碳冲击作用产生破裂，同时高压气体能够侵入裂隙，使之继续扩展[15]。
![]()
图 1 二氧化碳相变爆破系统示意Figure 1. Schematic diagram of carbon dioxide phase change blasting system![]()
图 2 二氧化碳相变爆破原理T—温度；P—压强；n—物质的量；V—体积；P0—泄能片强度Figure 2. Principle of carbon dioxide phase change explosion1.2 相变爆破器材与装备
根据储液管的形式和泄气形式，相变爆破致裂管可分为重复使用致裂管和一次性致裂管（图3）。重复使用致裂管采用管壁较厚（10 mm以上）的钢管作为储液管、管内的高压气体通过管底卸压片破坏释放；一次性致裂管管壁较薄（1.5 mm），采用管板工艺加固，致裂管侧面的焊缝为二氧化碳释放通道[16]，二氧化碳气体是从致裂管侧面释放而不是从致裂管底部释放，这种气体释放方式，也有助于消除“飞管”问题，同时使得气体膨胀压力更加均匀地作用在炮孔壁上。虽然这两种形式结构不同，但其工作原理是相似的[17]。
![]()
图 3 不同相变爆破装置结构及射流特征Figure 3. Jet characteristics of different phase change blasting devices传统相变爆破多采用重复使用的致裂管，其结构如图4所示。注液头内有一个充液阀和一个点火头，充液阀用于每次爆破后重新充装二氧化碳，点火头电极与外部导线连接后可对发热管加热。发热管是相变爆破过程中主要加热元件，内装化学物质，爆破时，点火头电极通电后产生的热能使化学物质发生迅速反应，释放大量的热量，使二氧化碳迅速汽化。定压卸能片在相变爆破过程中主要用于控制爆破压力，一旦汽化的二氧化碳压力超过爆破片强度，爆破片就发生破坏，气体从而得到释放，一般采用剪切型爆破片，这爆破片结构简单、安装方便，为保证密封效果，在爆破片和爆破筒主体之间加入起密封作用的垫片。排气头是二氧化碳是释放通道，设有排气孔，用于控制气体释放方式，依据其需求不同，可设置成不同形状的排气孔。储液管是相变爆破装置的主体，用来盛装高压的二氧化碳。
1.3 相变爆破能量
相变爆破的能量密度比密实的化学炸药小，对于相变爆破能量的准确计算对相变爆破效果的评估至关重要，能够量化研究该技术爆破能力。关于相变爆破能量的计算有采用压缩气体和水蒸气容器爆炸能量、采用射流速度计算爆炸能量、采用爆容计算爆破能量、采用爆破振动等效爆炸能量和采用真实状态方程计算爆破能量方式，不同计算方式确定的相变爆破能量特征见表1。
表 1 相变爆破能量计算方法Table 1. Calculation method for phase change blasting energy计算原理 计算公式 单位质量CO2爆破能量/kJ TNT当量/g 评价 压缩气体和水蒸气
容器爆炸能量[18-20]$E=\dfrac{P_1 V}{K-1}\left[1-\left(\dfrac{P_2}{P_1}\right)^{\dfrac{K-1}{K}}\right] $ 385.1～ 1256.4 90.6～295.6 使用广泛，参数易于获取，计算简单，爆炸能量计算结果较小，未考虑温度对爆炸能量的影响，容积对爆炸能量影响较大 射流速度计算爆炸能量[6, 21] $E=\dfrac{1}{2} v^2+\dfrac{k}{k-1} \dfrac{p_2}{\rho}\left(\dfrac{p_1}{p_2}\right)^{1 / k} $ 562.8～641.3 132.2～150.6 公式没有考虑到体积和质量对爆炸能量的影响，因此存在较大误差 爆容计算爆破能量[22] $E=\dfrac{Q W_{{\mathrm{T N T}}}}{Q_{{\mathrm{T N T}}}} $ 1687.3 397.0 爆容分析受到炸药配比等因素的影响，
结果不一定会十分精确爆破振动等效爆炸能量[23] $E=\dfrac{R M S}{R M S_{{\mathrm{T N T}}}} W_{{\mathrm{T N T}}} $ 1615 ～2827.5 380.0～430.0 从爆破振动间接测试，需现场爆破
试验测试，工程量大真实状态方式计算能量[24] $E=\displaystyle\int_1^2 p \mathrm{~d} V \cong-\Delta A=A_1-A_2=R\left(\phi_1 T_1-\phi_2 T_2\right) $ 327.3～ 1250 77.0～294.1 CO2热力学特性真实状态方程，计算精度高，计算复杂 注：E为相变爆破能量；P1为致裂管内绝对压力；P2为外部环境绝对压力；V为致裂管容积；K为CO2的绝热指数；v为气体出气口速度；ρ为气体密度；Q为相变爆破的爆容；QTNT为炸药的爆容；WTNT为单位TNT炸药爆破能量；RMS为相变爆破振动波的均方跟振幅；RMSTNT为TNT炸药爆炸振动波的均方跟振幅；A为亥姆霍兹自由能；R为气体常数；T为系统温度；φ为无量纲亥姆霍兹能量。 相变爆破影响因素众多，由于不同型号致裂管生产标准不统一，理论计算仅供参考[25]，现统计不同相变爆破方式下爆破能量大小，绘制出单根相变爆破能量分布和单位质量二氧化碳相变爆破能量分布如图5所示[15]，由于计算标准和试验方法的差异导致相变爆破能量数值存在差异，但呈现整体分布特征，单根致裂管爆破能量集中在500～1 000 kJ，部分相变爆破能量计算数值较高，在2 500 kJ左右。经过对不同相变爆破能量统计可得，单根致裂管爆破能量平均值为1 018.3 kJ，中位数为904 kJ。单位质量二氧化碳相变爆破能量平均值为825.5 kJ，中位数为618.3 kJ。较传统炸药爆破，相变爆破属于一种低能量致裂方式[26]。
2. 二氧化碳相变爆破致裂机理
2.1 二氧化碳相变射流传播特征
二氧化碳相变过后产生的射流是致裂煤岩体的根本动力[27-28]。致裂管二氧化碳射流与岩体作用特征如图6所示。相变致裂管与围岩之间多为不耦合结构，二氧化碳从致裂管射流过后在致裂管与钻孔孔壁之间衰减传播，之后与孔壁碰撞形成冲击压力并在岩体中激发应力波[29]，钻孔孔壁在应力波作用下产生初始裂隙，随后高压二氧化碳进入初始裂隙中对裂隙进一步起到扩展作用。
![]()
图 6 二氧化碳射流与岩体作用特征Figure 6. Characteristics of interaction between carbon dioxide jet and rock mass为确定相变爆破引发钻孔孔壁应力波特征，做出如下假设：钻孔中间充填物质为空气，因其作用较弱，将其按照真空处理；致裂管采用钢材制成，将其视为刚体；整个过程不考虑温度场的影响，二氧化碳在间隙发生等熵膨胀，遇到孔壁后产生冲击应力，在岩体中激发应力波。
高压二氧化碳与钻孔孔壁发生碰撞时，产生透射和反射现象，导致压力会增加n倍，因此在不耦合结构下，二氧化碳在钻孔孔壁上产生的初始压力[29]为
$$ {P_2} = n{P_{\mathrm{m}}}{\left( {\frac{{{d_{\mathrm{g}}}}}{{{d_0}}}} \right)^6} $$ (1) 式中：P2和Pm分别为孔壁压力和致裂管破裂压力；dg和d0分别为致裂管和致裂孔直径；n为压力增大倍数。
二氧化碳在碰撞孔壁后压力会迅速衰减，在钻孔周围形成准静态应力场，并且进入应力波产生的初始裂隙中并使裂纹扩展。钻孔中形成的初始准静态压力P3[30]为
$$ {P_3} = K{P_{\mathrm{m}}} $$ (2) 式中：K为二氧化碳气体衰减系数。
相变爆破的本质是二氧化碳爆破射流产生的动态应力波作用和准静态高压气体作用[31-32]。相变爆破作为一种物理爆破，利用二氧化碳相态变化产生的膨胀能量，其应力波传播速度小（仅为炸药爆破的1/10左右），冲击应力峰值也远小于炸药爆破[33]。水压致裂则是利用钻孔中高压水的作用对介质进行压力，其作用时间更长，压力也更低，可以将其视为准静态过程。不同致裂方式的压力特征见表2，其压力时间历程曲线如图7所示[34]。依据不同冲击应力特征及所处的动态应变率范围，将炸药爆破致裂、相变爆破致裂和水压致裂分别确定为高应力快速冲击、中等冲击和准静态冲击。相变爆破能量偏低，产生应力波峰值小，加载速率低，形成的气楔作用时长大，高压作用持续时间长，属于中等加载速率长持时载荷[35-36]。
表 2 不同致裂方式技术参数Table 2. Technical parameters of different methods致裂类型 升压时间/s 峰值压力/MPa 加载速率/(MPa·s−1) 总时长/s 炸药爆破 10−7 ＞104 ＞108 10−6 相变爆破 10−3 102 102～106 10−2 水压 102 10 ＜10−1 104 ![]()
图 7 不同致裂方式压力时间历程曲线Figure 7. Pressure time history curve of different fracturing methods2.2 中等加载速率下多重起裂机理
致裂钻孔在动态应力波作用下会产生初始断裂。致裂钻孔周围介质受到外部原始地应力作用和钻孔内部冲击压力P。如图8所示，钻孔周围介质受力可以通过2种状受力特征叠加获得[14]。
$$ \left\{ \begin{gathered} {\sigma _r} = \frac{{{R^2}}}{{{r^2}}}P + \frac{{\gamma H}}{2}\left( {1 + \lambda } \right)\left( {1 - \frac{{{R^2}}}{{{r^2}}}} \right) - \frac{{\gamma H}}{2}\left( {1 - \lambda } \right)\left( {1 - 4\frac{{{R^2}}}{{{r^2}}} + 3\frac{{{R^4}}}{{{r^4}}}} \right)\cos \;2\theta \\ {\sigma _\theta } = \frac{{{R^2}}}{{{r^2}}}P - \frac{{\gamma H}}{2}\left( {1 + \lambda } \right)\left( {1 + \frac{{{R^2}}}{{{r^2}}}} \right) - \frac{{\gamma H}}{2}\left( {1 - \lambda } \right)\left( {1 + 3\frac{{{R^4}}}{{{r^4}}}} \right)\cos\; 2\theta {\text{ }} \\ \end{gathered} \right. $$ (3) 式中：σr为径向应力；σθ为切向应力；θ为该点的水平夹角；λ为水平侧压系数；H为岩体埋深；γ为岩层平均容重；R为致裂孔半径；P为冲击压力；r为该点距离钻孔中心距离。
根据式3绘制出静水压力条件下钻孔周围受力特征如图9所示。冲击应力和地应力作用产生的径向应力方向相同起到叠加作用，致裂孔周围介质承受较大的压应力作用；冲击应力和地应力作用产生的切向应力作用相反起到相互削弱作用，但由于冲击应力数值较大、衰减较快，致裂周围较近处受到拉应力作用，随着距离增加，拉应力数值逐渐降低并逐渐转换为压应力。
由于煤岩体动态抗压强度远大于抗拉强度，相变爆破和地应力产生的应力峰值小于煤岩体抗压强度。钻孔周边介质受到径向压应力作用和切向拉应力作用，当应力超过其抗拉强度时，在钻孔周围发生拉断裂。
根据Mott理论[10]，当断裂产生过后，会在断裂周围产生卸载波并向四周传播，卸载波的传播能够有效地释放应力，也就是说，卸载波传播过的区域即使受到冲击应力的作用也会产生卸载，不会产生新的断裂，断裂仅能产生在卸载波未影响的区域。
取钻孔周围一小段微元体分析（图10），该微元体在应力作用下产生断裂，当应变大于断裂应变（εc）时产生断裂，发生断裂后应力得到充分卸载到0。
$$ \left\{ \begin{gathered} \sigma = {\sigma _{\mathrm{c}}}{\text{ }}\left( {\varepsilon < {\varepsilon _{\text{c}}}} \right) \\ \sigma = 0{\text{ }}\left( {\varepsilon \geqslant {\varepsilon _{\text{c}}}} \right) \\ \end{gathered} \right. $$ (4) 相变爆破卸载波传播速度与材料本身性质和运动参数相关，对于相同材料的不同冲击特征作用下钻孔周围介质起裂特征如图11所示。依据冲击破裂起裂机理和最终状态，将不同冲击形成的破坏分为3条不同的路径：①高应力快速冲击（Ⅰ→Ⅱ）；②准静态作用（Ⅰ→Ⅲ→Ⅳ）；③中等冲击（Ⅰ→Ⅲ→Ⅴ→Ⅵ→Ⅶ→Ⅷ）：
![]()
图 11 不同冲击作用下致起裂机理Figure 11. Crack initiation mechanism under different impact effects①产生不同形式断裂的核心是冲击波传播速率和应力峰值大小的差异。当冲击应力数值大、速率快时，应力作用直接超过其抗压强度，煤体产生压缩破坏，同时冲击速率快，在卸载波还未完全形成时，煤体已经完全破坏，卸载波无法形成有效卸载作用（Ⅰ→Ⅱ），在钻孔周围形成粉碎破裂区。
②随着冲击应力峰值降低，达不到煤体动态抗压强度，但产生的切向应力大于抗拉强度时，煤体在冲击应力形成的切向拉应力作用下产生初始断裂（A、B），初始断裂产生过后，卸载波随之从断裂处向四周传播（Ⅰ→Ⅲ）；如果冲击应力速率较慢，卸载波已经传至整个区域，煤体得到有效卸载，在钻孔内应力作用下也不会产生新的断裂，最终形成单条裂隙（Ⅰ→Ⅲ→Ⅳ）。
③如果冲击速率介于上述2种条件之间时，在产生初次断裂（Ⅰ→Ⅲ）之后，冲击波仍会持续对周围介质产生作用，冲击波传播快于卸载波传播，在卸载波未到达的区域由于拉应力的作用产生新的断裂（C、D）点（Ⅰ→Ⅲ→Ⅴ）；新的断裂继续产生卸载波（Ⅰ→Ⅲ→Ⅴ→Ⅵ）；如果卸载波还未完全扩展到整个区域时，仍旧会产生新的断裂（E、F）点（Ⅰ→Ⅲ→Ⅴ→Ⅵ→Ⅶ）；直至卸载波传播至整个区域。在这种周而复始的卸载波和冲击波相互作用下，煤体最终形成多重断裂形态，在钻孔周围形成多个起裂点。
二氧化碳相变爆破应力峰值较传统炸药爆破小，较水压致裂大；作用时长较传统炸药爆破大，较水压致裂小。采用不同冲击波类型探讨了不同冲击特征下煤体裂隙分布特征，利用数值模拟方式确定了不同冲击载荷下裂隙分布特征，并与相关研究[23]对比，得到如图12所示的不同冲击条件下煤体裂隙分布特征。数值模拟得出的结果和理论分析、实验得到效果具有一致性。
在高应力快速冲击致裂的传统爆破作用下，爆破产生的高强度冲击波远大于煤体抗压强度，在钻孔周围煤体中形成压缩粉碎区，爆破产生能量未能充分利用形成多裂隙状态；在相变爆破致裂条件下，冲击应力峰值降低，作用时长增加，钻孔周围的煤体产生拉破坏形成断裂，由于冲击作用时间较长，在钻孔周围形成多条裂隙，并有部分裂隙在应力持续作用下形成贯穿煤体的主裂隙；在水压致裂条件下，煤体受到准静态作用，仅在拉应力的作用产生单条裂隙。
2.3 高压气体多重裂隙扩展
二氧化碳相变射流产生的动态应力波会在致裂孔孔壁上产生初始的多重断裂，随后高压二氧化碳楔入，初始的多重裂隙会进一步扩展[37]。动态应力波形成的初始裂隙可以简化为楔形，多重楔形裂纹的扩展模型如图13所示，由于在高压气体作用下裂纹传播速度较小，将多重裂纹的传播过程视为准静态过程。
根据弹性理论，对称双翼裂纹扩展的裂隙宽度w[38]为
$$ \begin{split} w\left( \theta \right) = &\frac{{4\left( {1 - v} \right){L_{\mathrm{f}}}}}{{\pi G}}\int_\theta ^1 {\left[ {\int_0^{{\theta _2}} {\frac{{\left( {p\left( {{\theta _1}} \right) - \sigma } \right)}}{{\sqrt {\theta _2^2 - \theta _1^2} }}{f_{\mathrm{N}}}\left( {{\theta _1}} \right)f\left( {{\theta _1},\frac{{{L_{\mathrm{f}}}}}{{{r_{\mathrm{w}}}}}} \right){\mathrm{d}}{\theta _1}} } \right]} \\&\quad {\text{ }}\frac{{\theta _2^{}}}{{\sqrt {\theta _2^2 - \theta _{}^2} }}{f_{\mathrm{N}}}\left( {{\theta _2}} \right)f\left( {{\theta _2},\frac{{{L_{\mathrm{f}}}}}{{{r_{\mathrm{w}}}}}} \right){\mathrm{d}}{\theta _2} \\[-1pt] \end{split}$$ (5) 式中：Lf为裂纹长度；θ为裂纹上某一点的相对位置；θ1，θ2分别为裂纹扩展瞬间的相对长度和相对微段长度；v为泊松比；G为剪切模量；rw为钻孔半径；f (θ, Lf/rw)为考虑裂纹位置与裂纹形状、钻孔尺寸的权重系数，可以表示为
$$ f\left( {\theta ,\frac{{{L_{\mathrm{f}}}}}{{{r_{\mathrm{w}}}}}} \right) = 1 + 0.3\left( {1 - \theta } \right){\left( {\frac{{{r_{\mathrm{w}}}}}{{{r_{\mathrm{w}}} + {L_{\mathrm{f}}}}}} \right)^2} $$ (6) 为实现相变爆破致裂多重裂隙的计算模型，在上述对称双翼裂纹模型上通过引入考虑裂纹数目的因子[39]：
$$ {f_{\mathrm{N}}} = \frac{{{f_{{\mathrm{N}}{\theta } }} + {f_{{\mathrm{N}}{\theta} }}\dfrac{{{L_{\mathrm{f}}}{N_{\mathrm{f}}}}}{{\pi {r_{\mathrm{w}}}}}}}{{{f_{{\mathrm{N}}{\theta} }} + \dfrac{{{L_{\mathrm{f}}}{N_{\mathrm{f}}}}}{{\pi {r_{\mathrm{w}}}}}}} $$ (7) $$ {f_{{\mathrm{N}}{\theta} }} = 1 + \frac{\pi }{2}\left( {\frac{{2\sqrt {{N_{\mathrm{f}}} - 1} }}{{{N_{\mathrm{f}}}}} - 1} \right)\sqrt {\left( {1 - {\theta ^2}} \right)} $$ (8) 式中：Lf为裂纹数目。
该因子同时考虑到了裂纹和钻孔形状的影响，当裂纹长度相对于裂纹之间的距离较小时，该因子是统一的，而当裂纹长度较大时，该因子接近fNθ，与裂纹数目相关。
根据断裂力学理论和叠加原理，对称双翼型裂纹的应力强度因子可以由地应力作用下裂纹尖端的应力强度因子和高压气体作用下裂纹尖端的应力强度因子叠加得到：
$$ {K_{\mathrm{I}}} = 2\sqrt {\frac{{{L_{\mathrm{f}}}}}{\pi }} \int_0^1 {\frac{{\left[ {p\left( {\theta ,t} \right) - \sigma } \right]}}{{\sqrt {\left( {1 - {\theta ^2}} \right)} }}} f\left( {\theta ,\frac{{{L_{\mathrm{f}}}}}{{{r_{\mathrm{w}}}}}} \right){\mathrm{d}}\theta $$ (9) 同时在引入关于裂纹数目的权重因子f1Nθ，得到多重裂隙扩展作用下裂纹应力强度因子：
$$ {f_{1{\mathrm{N}}{\theta} }} = \left( {1 + \frac{\pi }{4}} \right)\frac{{\sqrt {{N_{\mathrm{f}}} - 1} }}{{{N_{\mathrm{f}}}}} $$ (10) 因此，二氧化碳相变爆破应力气体协同致裂过程如图14所示。二氧化碳相变会产生动态应力波和准静态的高压气体，由于动态应力波处于中等加载速率范围内，会在致裂孔孔壁形成多重初始裂隙，高压二氧化碳随后贯入到裂纹内，在气体压力函数和裂纹尖端应力强度因子相互作用下，驱动着裂纹不断向外扩展，最终到达破碎煤岩体的目的。
![]()
图 14 相变爆破应力气体协同致裂过程Figure 14. Phase change blasting stress gas coupling cracking process3. 相变爆破致裂影响因素
3.1 泄能方向
不同的二氧化碳相变爆破致裂装备会形成如图15所示的底部泄能和侧壁泄能等不同的方式。采用数值模拟方式研究不同泄能方式损伤破坏特征如图16所示[40]。
![]()
图 16 不同泄能方式损伤破坏特征Figure 16. Damage and failure characteristics of different energy release methods不同泄能方式形成的损伤破坏区域均成椭圆形分布，最大的损伤破坏区域半径为1.1 m。底部泄能处形成水平方向为长轴，垂直方向为短轴的椭圆形，在煤岩体底部作用均匀，而远离卸压区域的位置损伤较小，当煤岩体受到侧壁泄能相变爆破影响时，由于致裂管周围不均匀释放，在爆破孔壁周围形成非对称分布，在靠近泄能槽处煤岩体损伤破坏区域更大。
3.2 煤岩体性质
煤岩体性质会受到所处的地质构造和形成原因的影响，这会对相变爆破致裂效果产生影响。以煤岩体弹性模量、泊松比和抗压强度3种因素（分别为A1、B1、C1）建立正交分析，选取三水平，各个因素水平见表3。
表 3 煤岩体性质因素水平Table 3. Level table of coal properties factors水平 1 2 3 因素A1/GPa 1.0 3.0 5.0 因素B1 0.1 0.25 0.4 因素C1/MPa 5 15 25 根据正交表设置9种不同的方案，通过数值模拟方式研究不同方案下模型损伤破坏体积变化特征，分析各个因素的极差数据如图17所示[41]。随着因素A1弹性模量的增加，相变爆破破坏体积先增加后减小；随着因素B1泊松比的增加，相变爆破破坏体积逐步增加，但增加到一定值时增速减缓；随着因素C1抗压强度的增加，相变爆破破坏体积先增加后减小，对于较小强度煤体，由于相变爆破会形成大量的破碎区消耗能量，煤体破坏效果反而会稍微降低。
对于不同因素引起的数据差异进行方差分析，计算各个因素的误差见表4，其中F比是对应因素均方和与误差均方和的比值，F的比值越大，表明该因素对试验效果的影响程度越大。因此，对于爆破效果的影响因素由大到小分别为抗压强度、泊松比和弹性模型。对于不同性质煤岩体进行相变爆破致裂方案设计时，应动态调整相应方案。
表 4 煤岩体性质方差分析Table 4. Variance analysis of coal properties来源 平方和S 自由度 均方和V F比 因素A1 0.959 2 0.480 0.435 因素B1 3.041 2 1.521 1.378 因素C1 4.439 2 2.220 2.011 误差E 2.207 2 1.104 — 3.3 相变爆破参数
针对相变爆破参数，同样的采用正交分析研究二氧化碳充装量、泄能片爆破压力和致裂孔间距等3个因素（分别为A2、B2、C2）对致裂效果的影响，各个因素水平表见表5。
表 5 相变爆破参数因素水平Table 5. Level table of phase change blasting parameters水平 1 2 3 因素A2/kg 1.5 1.8 2.1 因素B2/MPa 120 150 180 因素C2/m 6 8 10 相变爆破参数因素影响下失效比例指标极差分析如图18所示[41]，由图可知，随着因素A2充装量的增加，煤体失效比例指标增加；随着因素B2爆破压力增加，煤体失效比例指标呈现先增加后减小的趋势，过高的爆破压力引起的能量高度集中反而不利于相变爆破致裂效果；随着因素C2致裂孔间距的增加，煤体失效比例指标减小，相变爆破致裂效果与致裂孔间距呈现负相关关系。
对相变爆破各参数进行方差分析见表6，对于相变爆破二氧化碳充装量、爆破压力和致裂孔间距3个因素，相变爆破致裂孔间距对于爆破参数影响最大，其次是二氧化碳充装量，爆破压力对于相变爆破效果影响最小。在相变爆破参数设计时，应着重调整相变爆破致裂孔间距以适应不同条件。
表 6 相变爆破参数方差分析Table 6. Variance analysis of phase change blasting parameters来源 平方和 自由度 均方和 F比 因素A2 3.3×10−4 2 1.65×10−4 3.11×10−3 因素B2 8.4×10−5 2 4.2×10−5 8.49×10−4 因素C2 7.1×10−3 2 3.55×10−3 6.70×10−2 误差E 0.106 2 0.053 — 3.4 其他影响因素
除了二氧化碳泄能口方向、煤岩体性质和爆破参数外，相变爆破还会受到初始地应力、钻孔布置参数和钻孔切槽特性等因素的影响。
在地下采矿过程中相变爆破通常受到初始地应力的影响，且随着开采深度的不断增加，地应力影响逐渐增加。因初始地应力会与爆破产生应力波产生叠加作用，会影响动态应力波的传播规律，进而影响裂纹传播特征。初始地应力产生的围压效应会对相变爆破致裂产生的裂隙起到抑制作用[42]，但地应力形成的主应力差会引导裂纹方向发生偏转[35, 43]。
多排钻孔进行相变爆破时，常用的布孔方式有矩形布孔和菱形布孔2种不同方式。在钻孔孔距和排距相同的条件下，菱形布孔较矩形布孔方式裂隙贯通范围更广[35, 44]。采用菱形布孔进行相变爆破后瓦斯抽采达标时间小于矩形布孔，致裂增透效果更优[45]。为进一步提升相变爆破效果，可以考虑通过设置控制孔（空孔）的方式提升相变爆破效果。控制孔可以为相变爆破提供辅助自由面，应力波在控制孔处叠加反射，加速裂纹的扩展[46]。但相变爆破依靠高压二氧化碳致裂岩样，当控制孔布置过多时，高压气体会沿着裂隙通过控制孔逸出，降低相变爆破致裂效果[47]。
为控制相变爆破裂纹扩展方向，相关学者结合切槽定向爆破方式，提出了二氧化碳切槽定向致裂技术[48]。对圆形致裂孔进行切槽后，切槽处会产生应力集中，同时降低了槽尖起裂应力，使得裂纹沿槽线方向起裂，对裂纹扩展方向起到了引导的作用。
4. 相变爆破致裂工程应用
4.1 高瓦斯煤体相变爆破致裂增透
二氧化碳相变爆破过程会形成多重裂隙，不仅会促进瓦斯的有效解析，而且能够为瓦斯提供有效的流动通道。依据图12所做的相变爆破裂隙特征，建立多重裂隙煤体瓦斯渗流特征如图19所示[14]。相变爆破过后裂隙区域煤体渗透率远大于原始煤体渗透率，相变爆破形成的多重裂隙形成瓦斯流动的有效通道，因其数目多，分布范围广，能够对整个煤体中的瓦斯起到卸压增透作用。
![]()
图 19 多重裂隙煤体瓦斯渗流特性Figure 19. Gas seepage characteristics of coal bodies with multiple fractures潞安集团高河煤矿主采3号煤层，原始瓦斯含量为11.16 m3/t，采用二氧化碳相变爆破进行卸压增透提高抽采效率，通过实施二氧化碳相变爆破预裂过后瓦斯抽采效果如图20所示，压裂钻孔瓦斯抽采浓度提高3倍以上，煤层透气性系数整体提高10～20倍，提高了瓦斯抽采效率，保障了矿井安全生产。
![]()
图 20 相变爆破致裂增透效果Figure 20. Permeability enhancement of phase change blasting fracturing4.2 相变爆破预裂提升块煤率
块煤率是煤矿井下开采的煤炭资源成块状的比率，是衡量煤炭品质的一个重要指标。在我国西北部赋存的浅埋煤层中，由于煤体本身硬度较大，综采时煤体不易从煤壁剥落[49]，导致采煤机效率低、截割速度慢、块煤率降低。采用人工干预手段降低煤体硬度、增加煤体裂隙发育程度，能够有效地提升工作面块煤率，降低采煤设备的比能耗和提高设备的工作效率。
采用PFC数值模拟软件，研究煤体预裂前后截割特征如图21所示[50]。煤体预裂过后，经过滚筒截割时，原本形成的大块煤进一步破碎为小块煤，形成的块煤间尺寸差异性降低，截割力分布更加均匀。
榆神矿区40105工作面回采4-2煤层，煤层硬度为3.0，在回采巷道两侧平行于工作面布置2排钻孔进行二氧化碳相变爆破，二氧化碳相变爆破过后，块煤粒级变化如图22所示[51]。相变爆破预裂过后粒度大于100 mm的煤块占比减小，而粒度小于100 mm块体粒度增大，与数值模拟结论一致。现场应用效果表明相变爆破能够有效提升回采的块煤率。
![]()
图 22 相变爆破预裂煤体块度变化Figure 22. Change in block size of pre cracked coal body in phase change blasting4.3 工作面相变爆破预裂顶板
工作面回采过程中上隅角悬顶不仅会引起瓦斯积聚，而且大面积悬顶的突然垮落会引起其他次生灾害的发生[52]。采用二氧化碳相变爆破的方式预裂顶板，能够在顶板中产生裂隙，在工作面回采过后上隅角区域会随之垮落。沙坪煤矿9204工作面相变爆破预裂效果如图23所示。通过钻孔窥视仪对顶板钻孔进行窥视，致裂前煤层顶板完整、裂隙不发育，致裂过后，致裂孔内出现多条裂隙。工作面推进过预裂区域后采空区顶板垮落较好。
![]()
图 23 相变爆破预裂顶板效果对比Figure 23. Comparison of effect of phase change blasting on pre-splitting the roof由于矿压较小，顶板坚硬完整，综放工作面在初采期间底板难以垮落。为了避免工作面悬顶大面积的突然垮落，减少初采期间放煤损失，在磁窑沟煤矿进行相变爆破的预裂试验。在初采开切眼位置及端头两巷布置相变爆破致裂孔，对工作面基本顶进行预裂，最终效果对比如图24所示。采用相变爆破致裂过后，工作面基本顶内产生贯通裂隙，破坏了顶板完整结构，随着开采的不断推进，工作面采空区全部垮落、正常放煤推进距离和初次来压步距均减小了50%左右。
![]()
图 24 相变爆破致裂初采放顶效果Figure 24. Effect of phase change blasting on initial mining and topping5. 相变爆破致裂研究展望
对于二氧化碳相变爆破致裂机理的研究得到了初步进展，同时在高瓦斯煤层卸压增透、块煤率提升及工作面顶板预裂等场景中得到成功应用。未来在相变爆破致裂机理方面应在致裂多速率研究、多尺度分析和多物理场耦合方面等方面进一步研究。同时开展延时相变爆破技术研究，拓宽相变爆破应用场景。
1）进一步明晰煤岩体多速率作用下起裂机理。各种致裂手段的本质是煤岩体在不同应力作用下破坏断裂特征，煤岩体裂隙分布特征与不同致裂手段形成的加载速率相关（图25）。进一步研究不同应变率作用下断裂特征与冲击应力的关系，能够为致裂方式的选择提供更科学理论基础。采用岩石动力学测试方法构建煤岩体在相变爆破中等速率条件下的本构关系，开展不同加载速率下中心圆孔受力破坏特征试验，建立煤岩体裂纹数目、破碎块度与加载速率之间的定量关系。
![]()
图 25 不同加载速率冲击破碎特征Figure 25. Impact fragmentation characteristics at different loading rates2）进行相变爆破致裂煤岩体多尺度分析。二氧化碳相变爆破致裂在宏观上表现为裂隙的扩展，在微观层面上，二氧化碳相变作用对于孔隙的影响也是不可忽略的重要因素[53]，另外，二氧化碳会与部分矿物发生反应，改变煤岩体矿物特征，也会影响煤岩体性质的变化，对于相变爆破致裂效果分析仍需要建立宏微观统一模型来评价。通过图像识别、微震监测反演等方式确定宏观裂隙的扩展特征，采用扫描电镜、CT扫描、核磁共振、低温氮吸附、压汞实验等方法建立致裂前后孔隙变化特征，通过XRD、FTIR、GCMS等仪器测定煤岩体矿物成分及微晶结构变化特征[3]，如图26所示。最终利用分形的方法建立不同表征方法跨尺度之间的联系，形成相变爆破致裂煤岩体多尺度模型。
3）研究相变爆破致裂过程中多物理场耦合特征。二氧化碳相变爆破过程中不仅涉及到煤岩体固体力学物理场，还包括二氧化碳流场变化特征[54]，同时相变过程会引发温度场相应变化[55]，二氧化碳与矿物的反应[56]涉及化学场，整个相变爆破致裂应从热−流−固−化等多场耦合特征进行进一步分析。另外，可以进一步引入相场损伤的方式来表明裂纹扩展特征。通过研究4个物理场之间的耦合关系（图27），建立多物理场PDE方程，采用COMSOL等多场耦合软件，研究不同物理场在相变爆破过程中演化规律，确定不同因素对于致裂效果影响大小，从不同物理场角度针对性提出提升相变爆破致裂效果方法。
![]()
图 27 相变爆破致裂多场耦合特征Figure 27. Multi field coupling characteristics of phase change blasting fracturing4）开展延时相变爆破致裂技术研究。目前，相变爆破多采用小规模、一次性爆破方式，单次破岩体量较小，随着电子雷管的开发与普及，使得大规模延时相变爆破技术具备了可行性[57]，能够达到提升爆破效果和降低爆破振动的双重目的，进一步拓宽相变爆破技术应用范围。煤岩体中原生静态裂隙和爆破产生的动态裂隙之间也会产生相互影响（图28），不同阶段的裂隙一方面会形成后续爆破的自由面，另一方面也可能成为二氧化碳泄气通道，需要研究裂隙与相变爆破致裂的互馈影响关系，并确定延时相变爆破网络时空演化规律，明确裂隙网络的演化特征及其与相变爆破参数的相关性，最终形成基于爆破参数合理设计的相变爆破裂隙网络调控方法。
![]()
图 28 延时相变爆破裂纹扩展研究Figure 28. Research of delayed phase change blasting crack propagation5）扩展二氧化碳相变爆破工程应用范围。相变爆破致裂虽然在一些工程中得到了成功应用，但多依据工程经验进行施工，相关技术规范还未形成，相变爆破致裂优势没有得到充分发挥。在相变爆破设备装置、参数设计规范和智能管控等方面还需要进一步研究，开发适用于不同工况的致裂管体和释放喷嘴，研究相变爆破直接定向致裂设备，形成相变爆破孔网布置和爆破参数设计规范，建立相变爆破致裂全过程智能管理监控系统，能够为二氧化碳相变爆破进一步扩展应用提供指导。
6. 结 论
1）二氧化碳相变爆破能量的本质是二氧化碳相变时物理膨胀。根据相变爆破装置可分为重复使用致裂管和一次性致裂管，二者工作原理相似，泄能方式不同。二氧化碳相变爆破能量采用不同能量计算方式均表明，相较于传统炸药爆破，相变爆破属于一种低能量爆破方式。
2）二氧化碳相变射流能够产生动态应力波和高压气体双重作用，煤岩体在中等冲击作用下，受到拉应力破坏产生径向初始断裂，并在冲击波和卸载波综合作用下形成多重起裂特征，高压气体在多重裂隙中进一步扩展，驱动裂隙向外扩展，相变爆破裂隙在应力气体协同作用下扩展演化。
3）泄能方向、煤岩体性质、相变爆破参数、初始地应力、钻孔布置参数和钻孔切槽特性等因素均会对二氧化碳相变爆破致裂效果产生影响。泄能方向对煤岩体破坏起到直接作用，引发非对称损伤破坏，煤岩体性质对于爆破效果的影响因素由大到小分别为抗压强度、泊松比和弹性模型，相变爆破参数对于爆破效果的影响因素由大到小分别为致裂孔间距、充装量和爆破压力，初始地应力、钻孔布置参数和钻孔切槽特性等影响裂纹发育扩展特征。
4）二氧化碳相变爆破致裂在各种工程中得到了成功应用。揭示了相变爆破多重裂隙渗流特征，多重裂隙形成了瓦斯流动的主要通道，对低渗煤体卸压增透能够起到有效作用；对比了预裂前后煤体截割特征，预裂爆破能够减少特大块煤体的产生，形成更加均匀的块煤；相变爆破预裂顶板产生的裂隙能够促进顶板的垮落。
5）未来对于相变爆破致裂基础理论的研究应主要集中在煤岩体多速率下起裂机理、相变爆破致裂多尺度分析、致裂过程中多物理场耦合特征和延时相变爆破技术，进一步揭示二氧化碳相变爆破致裂机理，拓宽二氧化碳相变爆破致裂工程应用场景。
-
![]()
图 1 山底河流域水文地质图和采样点
Figure 1. Hydrogeological map and sampling sites of Shandi River Basin
![]()
图 4 山底河流域AMD中 Ca2+与Mg2+关系
Figure 4. Relationship between Ca2+ and Mg2+ in AMD in Shanxi River Basin
![]()
图 5 山底河流域${\mathrm{SO}}_4^{2-} $质量浓度和pH随δ34S–${\mathrm{SO}}_4^{2-} $的变化
Figure 5. Changes of ${\mathrm{SO}}_4^{2-} $ and pH with δ34S–${\mathrm{SO}}_4^{2-} $ in Shanxi River Basin
![]()
图 6 样品α多样性的Chao1和Shannon指数曲线
Figure 6. Chao1 and Shannon index curves of Alpha diversity of samples
![]()
图 7 柳沟AMD样品优势菌属及丰度变化
Figure 7. Changes of dominant genera and abundance of bacteria in AMD samples from Liugou
![]()
图 8 庙沟露天矿AMD样品优势菌属及丰度变化
Figure 8. Changes of dominant genera and abundance of bacteria in AMD samples from Miaogou open-pit mine
![]()
图 9 小沟露天矿AMD样品优势菌属及丰度变化
Figure 9. Changes of dominant genera and abundance of bacteria in AMD samples from Xiaogou open-pit mine
![]()
图 10 榆林垴钻井AMD样品优势菌属相对丰度变化
Figure 10. Changes of dominant genera and abundance of bacteria in AMD samples from Yulinnao drilling
![]()
图 11 榆林垴钻井AMD样品微生物群落的CCA分析
Figure 11. CCA analysis of microbial communities of AMD samples from Yulinnao drilling
![]()
图 12 柳沟AMD样品微生物群落的CCA分析
Figure 12. CCA analysis of microbial communities of AMD samples from Liugou
![]()
图 13 庙沟露天矿AMD样品微生物群落的CCA分析
Figure 13. CCA analysis of microbial communities of AMD samples from Miaogou open pit mine
![]()
图 14 小沟露天矿AMD样品微生物群落的RDA分析
Figure 14. RDA analysis of microbial communities of AMD samples from Xiaogou open-pit mines
表 1 山底河流域AMD取样点基本情况
Table 1 Basic situation of AMD sampling sites in Shandi River Basin
样品 样品数量(组) 监测点类别 采空区煤层 氧化还原环境 水位标高/m 裂隙水 3 石炭二叠裂隙水 未开采 935.0～1 002.5 AMD01 12 采空区矿井溢出点 15号 缺氧 827.8 AMD02 12 露天矿坑水 12号 好氧 857.6 AMD03 12 露天矿矿渣淋滤水 15号 好氧 875.4 AMD04 12 采空区矿井水 15号 厌氧或缺氧 863.1 
下载: 导出CSV
-
[1] BALDI F,CLARK T,POLLACK S S,et al. Leaching of pyrites of various reactivities by thiobacillus ferrooxidans[J]. Applied and Environmental Microbiology,1992,58(6):1853−1856. doi: 10.1128/aem.58.6.1853-1856.1992
[2] SASAKI K,TSUNEKAWA M,OHTSUKA T,et al. The role of sulfur-oxidizing bacteria Thiobacillus thiooxidans in pyrite weathering[J]. Colloids and Surfaces A:Physicochemical and Engineering Aspects,1998,133(3):269−278.
[3] 盛益之. 酸性矿井水环境中Fe(Ⅱ)的氧化和微生物多样性研究[D]. 北京:中国地质大学(北京),2016. SHENG Yizhi. Low-pH Fe(Ⅱ) oxidation and microbial diversity in acid mine drainage environments[D]. Beijing:China University of Geosciences,2016.
[4] ZHANG X,NIU J J,LIANG Y L,et al. Metagenome-scale analysis yields insights into the structure and function of microbial communities in a copper bioleaching heap[J]. BMC Genetics,2016,17(1):1−12. doi: 10.1186/s12881-015-0265-z
[5] HAO C B,WEI P F,PEI L X,et al. Significant seasonal variations of microbial community in an acid mine drainage lake in Anhui Province,China[J]. Environmental Pollution,2017,223:507−516. doi: 10.1016/j.envpol.2017.01.052
[6] ETTAMIMI S,CARLIER J D,COX C J,et al. A meta-taxonomic investigation of the prokaryotic diversity of water bodies impacted by acid mine drainage from the São Domingos mine in southern Portugal[J]. Extremophiles,2019,23(6):821−834. doi: 10.1007/s00792-019-01136-1
[7] MÉNDEZ-GARCÍA C,PELÁEZ A I,MESA V,et al. Microbial diversity and metabolic networks in acid mine drainage habitats[J]. Frontiers in Microbiology,2015,6:475.
[8] 孙亚军,张莉,徐智敏,等. 煤矿区矿井水水质形成与演化的多场作用机制及研究进展[J]. 煤炭学报,2022,47(1):423−437. SUN Yajun,ZHANG Li,XU Zhimin,et al. Multi-field action mechanism and research progress of coal mine water quality formation and evolution[J]. Journal of China Coal Society,2022,47(1):423−437.
[9] 肖升木,谢学辉,柳建设,等. 矿物培养影响微生物群落的ERIC研究[J]. 高校地质学报,2007,13(4):662−668. doi: 10.3969/j.issn.1006-7493.2007.04.008 XIAO Shengmu,XIE Xuehui,LIU Jianshe,et al. Study on the changes of microbial ecology with ERIC method after cultivation by different mineral resources[J]. Geological Journal of China Universities,2007,13(4):662−668. doi: 10.3969/j.issn.1006-7493.2007.04.008
[10] BAKER B J,LUTZ M A,DAWSON S C,et al. Metabolically active eukaryotic communities in extremely acidic mine drainage[J]. Applied and Environmental Microbiology,2004,70(10):6264−6271. doi: 10.1128/AEM.70.10.6264-6271.2004
[11] 陈迪. 高硫煤废弃矿井微生物群落演替规律及铁硫代谢基因的功能预测[D]. 徐州:中国矿业大学,2020. CHEN Di. Succession rule of microbial community and functional prediction of iron-sulfur metabolic genes in high-sulfur abandoned coal mines[D]. Xuzhou:China University of Mining and Technology,2020.
[12] DRUSCHEL G K,BAKER B J,GIHRING T M,et al. Acid mine drainage biogeochemistry at iron mountain,California[J]. Geochemical Transactions,2004,5(2):13. doi: 10.1186/1467-4866-5-13
[13] LIU Y C,DONG Q,WU C,et al. Study of the succession of microbial communities for sulfur cycle response to ecological factors change in sediment of sewage system[J]. Environmental Science and Pollution Research,2015,22(12):9250−9259. doi: 10.1007/s11356-014-3934-0
[14] 吕保义,谢建云,郑喻,等. 稀土尾矿库周边地下水微生物的群落多样性研究[J]. 环境工程,2015,33(S1):101−104,116. LYU Baoyi ,XIE Jianyun ,ZHENG Yu,et al. Study on microbial biodiversity in the groundwater around tailings of rare earth[J]. Environmental Engineering,2015,33(S1):101−104,116.
[15] 唐春雷,梁永平,晋华,等. 山底河流域煤矿酸性矿井水野外监测[J]. 中国岩溶,2022,41(4):522−531. doi: 10.11932/karst20220402 TANG Chunlei,LIANG Yongping,JIN Hua,et al. Overview of field monitoring for acid mine water system of the coal mine in Shandi river basin[J]. Carsologica Sinica,2022,41(4):522−531. doi: 10.11932/karst20220402
[16] 石维芝,赵春红,梁永平,等. 煤矿酸性“老窑水”低Ca/Mg成因机制[J]. 中国岩溶,2022,41(4):511−521. doi: 10.11932/karst2022y10 SHI Weizhi ,ZHAO Chunhong ,LIANG Yongping ,et al. Genetic mechanism analysis of low Ca/Mg value of acid goaf water in coal mine drainage[J]. Carsologica Sinica,2022,41(4):511−521. doi: 10.11932/karst2022y10
[17] 唐春雷,梁永平,晋华,等. 山西娘子关泉群及其水的来源[J]. 中国岩溶,2022,41(2):174−182. doi: 10.11932/karst20220201 TANG Chunlei,LIANG Yongping,JIN Hua,et al. Niangziguan spring group in Shanxi Province and its water source[J]. Carsologica Sinica,2022,41(2):174−182. doi: 10.11932/karst20220201
[18] 唐春雷,申豪勇,赵春红,等. 古堆泉域岩溶地下水水化学特征及成因[J]. 环境科学,2023,44(9):4874−4883. TANG Chunlei,SHEN Haoyong,ZHAO Chunhong,et al. Hydrochemical characteristics and formation causes of ground karst water systems in Gudui spring catchment[J]. Environmental Science,2023,44(9):4874−4883.
[19] 赵春红,梁永平,卢海平,等. 娘子关泉域岩溶水 ${{\mathrm{SO}}_4^{2-} }$、δ34S特征及其环境意义[J]. 中国岩溶,2019,38(6):867−875. doi: 10.11932/karst20190604 ZHAO Chunhong,LIANG Yongping,LU Haiping,et al. Chemical characteristics and environmental significance of $ {{\mathrm{SO}}_4^{2-} } $ and sulfur isotope in the karst watershed of the Niangziguan spring,Shanxi Province[J]. Carsologica Sinica,2019,38(6):867−875. doi: 10.11932/karst20190604
[20] 梁永平,申豪勇,赵春红,等. 对中国北方岩溶水研究方向的思考与实践[J]. 中国岩溶,2021,40(3):363−380. LIANG Yongping,SHEN Haoyong,ZHAO Chunhong,et al. Thinking and practice on the research direction of karst water in northern China[J]. Carsologica Sinica,2021,40(3):363−380.
[21] 梁永平,赵春红,唐春雷,等. 山西娘子关泉水及污染成因再分析[J]. 中国岩溶,2017,36(5):633−640. LIANG Yongping,ZHAO Chunhong,TANG Chunlei,et al. Reanalysis of spring water and its pollution causes of the Niangziguan spring in Shanxi[J]. Carsologica Sinica,2017,36(5):633−640.
[22] 胡晓兵,方健聪,翟虎威,等. 硫氧同位素在识别辛安泉域岩溶水 ${{\mathrm{SO}}_4^{2-}} $来源中的应用[J]. 地质科技通报,2022,41(5):333−340. HU Xiaobing ,FANG Jiancong ,ZHAI Huwei ,et al. Application of sulfur and oxygen isotopes in identifying the source of sulfate in karst water from Xinan spring area[J]. Bulletin of Geological Science and Technology,2022,41(5):333−340.
[23] REN K,ZENG J,LIANG J P,et al. Impacts of acid mine drainage on karst aquifers:Evidence from hydrogeochemistry,stable sulfur and oxygen isotopes[J]. Science of The Total Environment,2021,761:143223. doi: 10.1016/j.scitotenv.2020.143223
[24] 唐春雷,梁永平,王维泰,等. 龙子祠泉域岩溶水水化学−同位素特征[J]. 桂林理工大学学报,2017,37(1):53−58. doi: 10.3969/j.issn.1674-9057.2017.01.007 TANG Chunlei,LIANG Yongping,WANG Weitai,et al. Hydrogeochemical and isotopic characteristics of the karst groundwater systems in Longzici spring basin[J]. Journal of Guilin University of Technology,2017,37(1):53−58. doi: 10.3969/j.issn.1674-9057.2017.01.007
[25] 唐春雷,赵春红,申豪勇,等. 娘子关泉群水化学特征及成因[J]. 环境科学,2021,42(3):1416−1423. TANG Chunlei,ZHAO Chunhong,SHEN Haoyong,et al. Chemical characteristics and causes of groups water in Niangziguan spring[J]. Environmental Science,2021,42(3):1416−1423.
[26] ULLRICH S R,POEHLEIN A,LEVICÁN G,et al. Iron targeted transcriptome study draws attention to novel redox protein candidates involved in ferrous iron oxidation in “Ferrovum” sp. JA12[J]. Research in Microbiology,2018,169(10):618−627. doi: 10.1016/j.resmic.2018.05.009
[27] LIN H,TANG Y L,DONG Y B,et al. Characterization of heavy metal migration,the microbial community,and potential bioremediating Genera in a waste-rock pile field of the largest copper mine in Asia[J]. Journal of Cleaner Production,2022,351:131569. doi: 10.1016/j.jclepro.2022.131569
[28] WATANABE T,KATAYANAGI N,AGBISIT R,et al. Influence of alternate wetting and drying water-saving irrigation practice on the dynamics of Gallionella-related iron-oxidizing bacterial community in paddy field soil[J]. Soil Biology and Biochemistry,2021,152:108064. doi: 10.1016/j.soilbio.2020.108064
[29] LI Y Y,LIU L,WANG H J. Mixotrophic denitrification for enhancing nitrogen removal of municipal tailwater:Contribution of heterotrophic/sulfur autotrophic denitrification and bacterial community[J]. Science of The Total Environment,2022,814:151940. doi: 10.1016/j.scitotenv.2021.151940
[30] WANG Y F,WU G X,ZHENG X N,et al. Synergistic ammonia and nitrate removal in a novel pyrite-driven autotrophic denitrification biofilter[J]. Bioresource Technology,2022,355:127223. doi: 10.1016/j.biortech.2022.127223
[31] LIN Z Q,PANG S M,ZHOU Z,et al. Novel pathway of acephate degradation by the microbial consortium ZQ01 and its potential for environmental bioremediation[J]. Journal of Hazardous Materials,2022,426:127841. doi: 10.1016/j.jhazmat.2021.127841
[32] MURAVYOV M,PANYUSHKINA A,FOMCHENKO N. Effect of copper/nickel ratio on the efficiency of biobeneficiation of bulk copper-nickel sulfide concentrates[J]. Minerals Engineering,2022,182:107586. doi: 10.1016/j.mineng.2022.107586
[33] BAO Y P,GUO C L,WANG H,et al. Fe- and S-metabolizing microbial communities dominate an AMD-contaminated river ecosystem and play important roles in Fe and S cycling[J]. Geomicrobiology Journal,2017,34(8):695−705. doi: 10.1080/01490451.2016.1243596
[34] ORMEÑO-ORRILLO E,MARTÍNEZ-ROMERO E. A genomotaxonomy view of the Bradyrhizobium genus[J]. Frontiers in Microbiology,2019,10:1334. doi: 10.3389/fmicb.2019.01334
[35] LI A H,LIU H C,HOU W G,et al. Pseudorhodobacter sinensis sp. nov. and Pseudorhodobacter aquaticus sp. nov. ,isolated from crater lakes[J]. International Journal of Systematic and Evolutionary Microbiology,2016,66(8):2819−2824. doi: 10.1099/ijsem.0.001061
[36] YAN Z Z,MENG H J,ZHANG Q Q,et al. Effects of cadmium and flooding on the formation of iron plaques,the rhizosphere bacterial community structure,and root exudates in Kandelia obovata seedlings[J]. Science of The Total Environment,2022,851:158190. doi: 10.1016/j.scitotenv.2022.158190
[37] CLYDE E J,CHAMPAGNE P,JAMIESON H E,et al. The use of a passive treatment system for the mitigation of acid mine drainage at the Williams Brothers Mine (California):Pilot-scale study[J]. Journal of Cleaner Production,2016,130:116−125. doi: 10.1016/j.jclepro.2016.03.145
[38] OTHMAN A,SULAIMAN A,SULAIMAN S K. Carbide lime in acid mine drainage treatment[J]. Journal of Water Process Engineering,2017,15:31−36. doi: 10.1016/j.jwpe.2016.06.006
[39] NAICKER K,CUKROWSKA E,MCCARTHY T S. Acid mine drainage arising from gold mining activity in Johannesburg,South Africa and environs[J]. Environmental Pollution,2003,122(1):29−40. doi: 10.1016/S0269-7491(02)00281-6
[40] 赵峰华. 煤矿酸性水地球化学[M]. 北京:煤炭工业出版社,2005. [41] 岳梅,赵峰华,任德贻. 煤矿酸性水水化学特征及其环境地球化学信息研究[J]. 煤田地质与勘探,2004,32(3):46−49. doi: 10.3969/j.issn.1001-1986.2004.03.016 YUE Mei,ZHAO Fenghua,REN Deyi. The environment geochemistry information of the coal mine acid mining drainage[J]. Coal Geology & Exploration,2004,32(3):46−49. doi: 10.3969/j.issn.1001-1986.2004.03.016
[42] SIMATE G S,NDLOVU S. Acid mine drainage:Challenges and opportunities[J]. Journal of Environmental Chemical Engineering,2014,2(3):1785−1803. doi: 10.1016/j.jece.2014.07.021
[43] RODRÍGUEZ-GALÁN M,BAENA-MORENO F M,VÁZQUEZ S,et al. Remediation of acid mine drainage[J]. Environmental Chemistry Letters,2019,17(4):1529−1538. doi: 10.1007/s10311-019-00894-w
[44] LUO Z H,LI Q,CHEN N,et al. Genome-resolved metagenomics reveals depth-related patterns of microbial community structure and functions in a highly stratified,AMD overlaying mine tailings[J]. Journal of Hazardous Materials,2023,447:130774. doi: 10.1016/j.jhazmat.2023.130774
[45] YANG D,FAN R,GREET C,et al. Microfluidic screening to study acid mine drainage[J]. Environmental Science & Technology,2020,54(21):14000−14006.
[46] HUANG Y,LI X T,JIANG Z,et al. Key factors governing microbial community in extremely acidic mine drainage (pH<3)[J]. Frontiers in Microbiology,2021,12:761579. doi: 10.3389/fmicb.2021.761579
[47] SHE Z X,PAN X,WANG J,et al. Vertical environmental gradient drives prokaryotic microbial community assembly and species coexistence in a stratified acid mine drainage lake[J]. Water Research,2021,206:117739. doi: 10.1016/j.watres.2021.117739
[48] CHEN J J,LIU Y L,DIEP P,et al. Genomic analysis of a newly isolated acidithiobacillus ferridurans JAGS strain reveals its adaptation to acid mine drainage[J]. Minerals,2021,11(1):74. doi: 10.3390/min11010074
[49] JIN D C,WANG X M,LIU L L,et al. A novel approach for treating acid mine drainage through forming schwertmannite driven by a mixed culture of Acidiphilium multivorum and Acidithiobacillus ferrooxidans prior to lime neutralization[J]. Journal of Hazardous Materials,2020,400:123108. doi: 10.1016/j.jhazmat.2020.123108
[50] BOND P L,DRUSCHEL G K,BANFIELD J F. Comparison of acid mine drainage microbial communities in physically and geochemically distinct ecosystems[J]. Applied and Environmental Microbiology,2000,66(11):4962−4971. doi: 10.1128/AEM.66.11.4962-4971.2000
[51] 高旭波,潘振东,龚培俐,等. 微生物诱导碳酸盐岩沉淀过程及作用机理[J]. 中国岩溶,2022,41(3):441−452. doi: 10.11932/karst20220311 GAO Xubo ,PAN Zhendong ,GONG Peili ,et al. Process and mechanism of microbial induced carbonate precipitation[J]. Carsologica Sinica,2022,41(3):441−452. doi: 10.11932/karst20220311
[52] KUANG J L,HUANG L N,CHEN L X,et al. Contemporary environmental variation determines microbial diversity patterns in acid mine drainage[J]. The ISME Journal,2012,7(5):1038−1050.
[53] WANG M M,WANG X N,ZHOU S N,et al. Strong succession in prokaryotic association networks and community assembly mechanisms in an acid mine drainage-impacted riverine ecosystem[J]. Water Research,2023,243:120343. doi: 10.1016/j.watres.2023.120343
计量
- 文章访问数: 30
- HTML全文浏览量: 14
- PDF下载量: 11
