GEOFLUIDSgydF4y2Ba GeofluidsgydF4y2Ba 1468 - 8123gydF4y2Ba 1468 - 8115gydF4y2Ba HindawigydF4y2Ba 10.1155 / 2018/8127207gydF4y2Ba 8127207gydF4y2Ba 研究文章gydF4y2Ba 气体渗流模型研究对减少气体上保护煤层开采过程中积累的风险gydF4y2Ba http://orcid.org/0000 - 0001 - 8673 - 3077gydF4y2Ba 张gydF4y2Ba 坎昆gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba http://orcid.org/0000 - 0001 - 6561 - 1011gydF4y2Ba 张gydF4y2Ba LeigydF4y2Ba 3gydF4y2Ba 李gydF4y2Ba MingxuegydF4y2Ba 4gydF4y2Ba http://orcid.org/0000 - 0001 - 9847 - 851 xgydF4y2Ba 王gydF4y2Ba 程ydF4y2Ba 2gydF4y2Ba 5gydF4y2Ba 王gydF4y2Ba 温gydF4y2Ba 1gydF4y2Ba 北京重点实验室精确Intergrown能源和资源的开采gydF4y2Ba 中国矿业大学科技(北京)gydF4y2Ba 北京100083年gydF4y2Ba 中国gydF4y2Ba cumtb.edu.cngydF4y2Ba 2gydF4y2Ba 国家重点实验室的煤炭资源安全开采(CUMT)gydF4y2Ba 徐州gydF4y2Ba 江苏221116年gydF4y2Ba 中国gydF4y2Ba cumt.edu.cngydF4y2Ba 3gydF4y2Ba 深煤炭资源教育部重点实验室gydF4y2Ba 矿业学院gydF4y2Ba 中国矿业大学&技术gydF4y2Ba 徐州gydF4y2Ba 江苏221116年gydF4y2Ba 中国gydF4y2Ba cumt.edu.cngydF4y2Ba 4gydF4y2Ba 物理学院的gydF4y2Ba 中国矿业大学&技术gydF4y2Ba 徐州gydF4y2Ba 江苏221116年gydF4y2Ba 中国gydF4y2Ba cumt.edu.cngydF4y2Ba 5gydF4y2Ba 矿业学院gydF4y2Ba 贵州大学gydF4y2Ba 贵阳550025gydF4y2Ba 中国gydF4y2Ba gzu.edu.cngydF4y2Ba 2018年gydF4y2Ba 23gydF4y2Ba 12gydF4y2Ba 2018年gydF4y2Ba 2018年gydF4y2Ba 10gydF4y2Ba 05年gydF4y2Ba 2018年gydF4y2Ba 06gydF4y2Ba 09年gydF4y2Ba 2018年gydF4y2Ba 23gydF4y2Ba 12gydF4y2Ba 2018年gydF4y2Ba 2018年gydF4y2Ba 版权©2018年坎昆Zhang et al。gydF4y2Ba 这是一个开放的文章在知识共享归属许可下发布的,它允许无限制的使用,分布和繁殖在任何媒介,提供最初的工作是正确的引用。gydF4y2Ba

保护煤层开采(PCSM)是一个最重要的缓解措施的地区爆发深部煤层开采过程中,在中国拥有高爆发的风险。在PCSM过程中,甲烷浓度超过极限的现象通常发生在工作面。是至关重要的理解从保护煤层工作面气体排放影响因素(PSWF)和获得一个方程来确定气体排放的数量。气体渗流模型(GSM)开发过程中模拟气体排放上PCSM。在这项研究中,制定一个方程来确定气体的数量而眠的缝成PSWF保护。这些方程已经开发利用扩散菲克第二定律和达西流定律。渗透率和压力之间的关系描述了处于弹塑性状态,和围岩的力学影响。它可以得出从GSM的初始气体压力保护焊缝,夹层岩石的特点,和通风PSWF压力是主要的影响因素从保护煤层的解吸气体排放到PSWF。发达GSM测试计算气体排放数量从PCSM过程利用一个煤矿的实际地质条件数据,中国位于汉城。结果显示与实地测量获得伟大的协议,这是由结合通风空气沼气量的拟合曲线PSWF夹层间距。 A loss coefficient ( δgydF4y2Ba)为1.012×10gydF4y2Ba−3gydF4y2Ba在这项研究中获得的。gydF4y2Ba

北京市自然科学基金gydF4y2Ba 8184082gydF4y2Ba 中国国家自然科学基金gydF4y2Ba 51704274gydF4y2Ba 51502339gydF4y2Ba 51704278gydF4y2Ba 一流的专业施工贵州省的本科学校gydF4y2Ba SJZY2017006gydF4y2Ba 中国矿业大学和技术gydF4y2Ba SKLCRSM16KF 06gydF4y2Ba 年轻的科学家和技术人员增长贵州教育部门的项目gydF4y2Ba 钱他肯塔基州字[2018]115gydF4y2Ba
1。介绍gydF4y2Ba

煤与瓦斯突出事故是一个复杂的动态的现象,在中国经常发生在煤矿。因为它的突然出现和强烈的破坏性,事故危及煤矿生产和人员的安全。它也会导致巨大的经济损失矿业公司(gydF4y2Ba 1gydF4y2Ba- - - - - -gydF4y2Ba 3gydF4y2Ba]。为了解决这个问题,许多研究已经在全国范围内进行的。这些研究阐明了在煤与瓦斯突出机制;此外,爆发风险指标和缓解措施提出了在这些研究[gydF4y2Ba 4gydF4y2Ba- - - - - -gydF4y2Ba 6gydF4y2Ba]。优先保护煤层开采(PCSM)是用来释放煤层压力,提高渗透率。这是在中国的一个主要突出缓解措施,主要是因为气体含量及其相关的爆发风险相对较低。通过提取出保护焊缝,相关的压力释放后,它大大增加了保护煤层的渗透性。减压是天然气开采增强的结果(gydF4y2Ba 7gydF4y2Ba- - - - - -gydF4y2Ba 9gydF4y2Ba]。在保护煤层开采(PCSM),甲烷浓度超过极限的现象通常发生在保护煤层工作面(PSWF)。因为seam提供了保护气体减压,大量的气体和流入PSWF眠。这是甲烷浓度超过限制的主要原因在PCSM的过程。保证安全开采和提高甲烷采矿、压力释放气体分离技术是目前使用gydF4y2Ba 10gydF4y2Ba- - - - - -gydF4y2Ba 12gydF4y2Ba]。目前,更多的是注意保护煤层的卸压效果,但因素导致甲烷浓度超过极限时忽略PCSM [gydF4y2Ba 13gydF4y2Ba,gydF4y2Ba 14gydF4y2Ba]。定性的方法主要用于估计压力释放气体渗流,这并不适用于瓦斯抽放和控制(gydF4y2Ba 15gydF4y2Ba- - - - - -gydF4y2Ba 17gydF4y2Ba]。除此之外,几乎没有研究压力释放气体渗流的主要影响因素。它也不利于提取和控制压力释放气体。因此,压力释放气体的渗流特征的定量分析具有重要意义的布局瓦斯抽放钻孔和保护煤层的开采参数。gydF4y2Ba

本文总结了先前的研究在煤和岩石渗透率。应力-应变关系、压力变化,和失败的地板已经决定PCSM的过程。根据评审结果,GSM及其方程建立了用于计算气体的数量而眠PCSM期间保护缝到保护煤层。确定的影响因素也通过分析渗流过程。最后,对模型进行了验证与数据从一个原位通风空气中的甲烷PSWF煤矿,位于汉城,中国。结果表明,现有的GSM完全赞同实地测量获得。gydF4y2Ba

2。之前的研究在渗透率的演化规律gydF4y2Ba

自从达西定律的命题,岩石渗透率的研究一直是一个有吸引力的和正在进行的研究领域。经过多年的开发和研究,建立了渗透率模型的两种主要类型:porosity-permeability模型和stress-permeability模型。这些模型由Palmer-Mansoori(出版社)和灰色模型,分别(gydF4y2Ba 18gydF4y2Ba]。在煤的渗透率模型提出的灰色(gydF4y2Ba 19gydF4y2Ba),地面和扩张/收缩应力引起的吸附/解吸气体是全面调查。gydF4y2Ba (1)gydF4y2Ba σgydF4y2Ba hgydF4y2Ba egydF4y2Ba −gydF4y2Ba σgydF4y2Ba hgydF4y2Ba 0gydF4y2Ba egydF4y2Ba =gydF4y2Ba −gydF4y2Ba vgydF4y2Ba 1gydF4y2Ba −gydF4y2Ba vgydF4y2Ba pgydF4y2Ba −gydF4y2Ba pgydF4y2Ba 0gydF4y2Ba +gydF4y2Ba EgydF4y2Ba 1gydF4y2Ba −gydF4y2Ba vgydF4y2Ba ΔgydF4y2Ba εgydF4y2Ba 年代gydF4y2Ba ΔgydF4y2Ba pgydF4y2Ba 年代gydF4y2Ba ΔgydF4y2Ba pgydF4y2Ba 年代gydF4y2Ba 。gydF4y2Ba

在这里,gydF4y2Ba σgydF4y2Ba hgydF4y2Ba egydF4y2Ba 是有效的水平应力,它是表达MPa;gydF4y2Ba σgydF4y2Ba hgydF4y2Ba 0gydF4y2Ba egydF4y2Ba 是初始有效的水平应力,这是表达MPa;gydF4y2Ba ΔgydF4y2Ba pgydF4y2Ba 年代gydF4y2Ba 吸附平衡压力的变化,表达了MPa;gydF4y2Ba ΔgydF4y2Ba εgydF4y2Ba 年代gydF4y2Ba 单元的应变变化,这是由于吸附平衡压力的变化;gydF4y2Ba vgydF4y2Ba 是泊松比;gydF4y2Ba pgydF4y2Ba 气体压力和gydF4y2Ba pgydF4y2Ba 0gydF4y2Ba 初始气体压力,MPa表达;和gydF4y2Ba EgydF4y2Ba 是平均绩点的弹性模量。gydF4y2Ba

所示(gydF4y2Ba 2gydF4y2Ba),Seidle et al。gydF4y2Ba 20.gydF4y2Ba提出了一个更简洁的表达。gydF4y2Ba (2)gydF4y2Ba kgydF4y2Ba =gydF4y2Ba kgydF4y2Ba 0gydF4y2Ba egydF4y2Ba −gydF4y2Ba 3gydF4y2Ba cgydF4y2Ba fgydF4y2Ba σgydF4y2Ba −gydF4y2Ba σgydF4y2Ba 0gydF4y2Ba 。gydF4y2Ba

在此,gydF4y2Ba cgydF4y2Ba fgydF4y2Ba 是夹板体积压缩系数对有效应力的变化正常夹,gydF4y2Ba kgydF4y2Ba 0gydF4y2Ba 在初始有效应力是初始渗透率夹板吗gydF4y2Ba σgydF4y2Ba 0gydF4y2Ba ,gydF4y2Ba kgydF4y2Ba 是夹板渗透率有效应力gydF4y2Ba σgydF4y2Ba 。gydF4y2Ba

方程(gydF4y2Ba 2gydF4y2Ba)被用于一些实验,进行拟合回归方程。所示(gydF4y2Ba 3gydF4y2Ba),任和爱德华兹(gydF4y2Ba 21gydF4y2Ba)提出,渗透率和压力关系与完整岩石有关。gydF4y2Ba (3)gydF4y2Ba kgydF4y2Ba =gydF4y2Ba kgydF4y2Ba 0gydF4y2Ba egydF4y2Ba −gydF4y2Ba 0.25gydF4y2Ba σgydF4y2Ba −gydF4y2Ba σgydF4y2Ba 0gydF4y2Ba 。gydF4y2Ba

一来et al。gydF4y2Ba 22gydF4y2Ba)总结了岩石断裂表面的渗透性之间的关系和主应力。方程(gydF4y2Ba 4gydF4y2Ba)用于计算断裂表面的渗透性。gydF4y2Ba (4)gydF4y2Ba kgydF4y2Ba =gydF4y2Ba kgydF4y2Ba 1gydF4y2Ba σgydF4y2Ba 1gydF4y2Ba +gydF4y2Ba σgydF4y2Ba 3gydF4y2Ba 2gydF4y2Ba 米gydF4y2Ba 。gydF4y2Ba

在这里,gydF4y2Ba σgydF4y2Ba 1gydF4y2Ba 和gydF4y2Ba σgydF4y2Ba 3gydF4y2Ba 最小和最大的压力,分别表示的MPa;gydF4y2Ba kgydF4y2Ba 1gydF4y2Ba 内在渗透表现在米gydF4y2Ba2gydF4y2Ba当gydF4y2Ba σgydF4y2Ba 1gydF4y2Ba +gydF4y2Ba σgydF4y2Ba 3gydF4y2Ba /gydF4y2Ba 2gydF4y2Ba =gydF4y2Ba 1gydF4y2Ba 米gydF4y2Ba PgydF4y2Ba 一个gydF4y2Ba ;和gydF4y2Ba 米gydF4y2Ba 是一个材料参数回归分析发现。在一项由一来et al .,−0.8616和2.613×10的值gydF4y2Ba−13gydF4y2Ba米gydF4y2Ba2gydF4y2Ba确定了gydF4y2Ba 米gydF4y2Ba 和gydF4y2Ba kgydF4y2Ba 1gydF4y2Ba 分别gydF4y2Ba 。gydF4y2Ba(gydF4y2Ba 22gydF4y2Ba]。gydF4y2Ba

帕默提出的最具代表性的渗透率模型和MansoorigydF4y2Ba 23gydF4y2Ba和帕默et al。gydF4y2Ba 24gydF4y2Ba)是表示为gydF4y2Ba (5)gydF4y2Ba φgydF4y2Ba φgydF4y2Ba 0gydF4y2Ba =gydF4y2Ba 1gydF4y2Ba +gydF4y2Ba cgydF4y2Ba 米gydF4y2Ba φgydF4y2Ba 0gydF4y2Ba pgydF4y2Ba −gydF4y2Ba pgydF4y2Ba 0gydF4y2Ba +gydF4y2Ba εgydF4y2Ba 1gydF4y2Ba φgydF4y2Ba 0gydF4y2Ba KgydF4y2Ba 米gydF4y2Ba −gydF4y2Ba 1gydF4y2Ba pgydF4y2Ba pgydF4y2Ba εgydF4y2Ba +gydF4y2Ba pgydF4y2Ba −gydF4y2Ba pgydF4y2Ba 0gydF4y2Ba pgydF4y2Ba εgydF4y2Ba +gydF4y2Ba pgydF4y2Ba 0gydF4y2Ba ,gydF4y2Ba cgydF4y2Ba 米gydF4y2Ba =gydF4y2Ba ggydF4y2Ba 米gydF4y2Ba −gydF4y2Ba γgydF4y2Ba KgydF4y2Ba 米gydF4y2Ba +gydF4y2Ba fgydF4y2Ba −gydF4y2Ba 1gydF4y2Ba ,gydF4y2Ba kgydF4y2Ba kgydF4y2Ba 0gydF4y2Ba =gydF4y2Ba φgydF4y2Ba φgydF4y2Ba 0gydF4y2Ba 3gydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba φgydF4y2Ba 孔隙度;gydF4y2Ba φgydF4y2Ba 0gydF4y2Ba 初始孔隙度;gydF4y2Ba KgydF4y2Ba 是GPa的体积弹性模量;gydF4y2Ba 米gydF4y2Ba 是GPa的约束轴向弹性模量;gydF4y2Ba γgydF4y2Ba 是固体压缩,它也被称为粮食压缩性;gydF4y2Ba pgydF4y2Ba 年代gydF4y2Ba 标准条件下的大气压力,表达MPa;和gydF4y2Ba fgydF4y2Ba 分数的值在0和1之间变化。值为0和1建议垂直和水平楔子,分别。gydF4y2Ba

原则上,porosity-permeability或stress-permeability模型被定义为有效应力的变化,从而改变煤和岩石断裂和变形的内部结构矩阵,甚至导致失败。因而岩体的渗透性的影响。这也是评价的主要理论原则渗透增强和保护煤层开采的卸压。程等。gydF4y2Ba 25gydF4y2Ba)提出了一个理论模型渗透率演化和减压的煤炭,其中包含气体PCSM深处。现有的地应力在深缝,有效应力变化,被认为是,它渗透直接或间接控制。此外,现有的骨折以及新的裂缝进一步扩展,由PCSM的过程。煤炭矩阵变形对渗透率变化的影响可以忽略。通过假设有效应力、气体吸附/解吸和应变为变量,建立了提高煤渗透率理论模型(gydF4y2Ba 25gydF4y2Ba]:gydF4y2Ba (6)gydF4y2Ba kgydF4y2Ba kgydF4y2Ba 0gydF4y2Ba =gydF4y2Ba 1gydF4y2Ba +gydF4y2Ba εgydF4y2Ba VgydF4y2Ba 0gydF4y2Ba −gydF4y2Ba εgydF4y2Ba VgydF4y2Ba +gydF4y2Ba fgydF4y2Ba 米gydF4y2Ba εgydF4y2Ba 米gydF4y2Ba 0gydF4y2Ba −gydF4y2Ba εgydF4y2Ba 米gydF4y2Ba φgydF4y2Ba 0gydF4y2Ba 3gydF4y2Ba ×gydF4y2Ba 1gydF4y2Ba −gydF4y2Ba φgydF4y2Ba 0gydF4y2Ba 1gydF4y2Ba −gydF4y2Ba φgydF4y2Ba 0gydF4y2Ba −gydF4y2Ba εgydF4y2Ba VgydF4y2Ba 0gydF4y2Ba −gydF4y2Ba εgydF4y2Ba VgydF4y2Ba +gydF4y2Ba fgydF4y2Ba 米gydF4y2Ba εgydF4y2Ba 米gydF4y2Ba 0gydF4y2Ba −gydF4y2Ba εgydF4y2Ba 米gydF4y2Ba 2gydF4y2Ba 。gydF4y2Ba

在这里,gydF4y2Ba εgydF4y2Ba 是应变,下标“0”指的是起始点,下标“m”是指煤炭矩阵,下标“V”是指体积,然后呢gydF4y2Ba fgydF4y2Ba 米gydF4y2Ba 是影响因素,范围在0和1之间,和煤的解吸和收缩矩阵断裂应变的影响。gydF4y2Ba

因为渗透率和压力之间的关系是不同的断裂和完整的煤岩体,岩体的失败保护和保护之间的接缝被认为是在任何研究中,这是在PCSM与GSM相关联。gydF4y2Ba

3所示。气体渗流模型上保护煤层开采gydF4y2Ba

气体渗流模型(GSM)与上层PCSM成立,主要取决于气体源和渗流媒体。气体来源主要是依赖于受保护的缝,和渗流媒体很大程度上与岩体层间的关联。如图gydF4y2Ba 1gydF4y2Ba,王gydF4y2Ba 26gydF4y2Ba)确定各种类型的煤层气利用数值模拟迁移。在保护煤层开采过程中,大量的气体被而眠的底部保护焊缝,提供减压以及煤和岩石的扩张。这是天然气的主要来源,流入PSWF通过层间的岩体。gydF4y2Ba

气体传输机制在水库(gydF4y2Ba 26gydF4y2Ba]。gydF4y2Ba

不同的失败与岩体层间的关联机制,取决于其岩性及其距离保护的地板接缝。这个实验相应的透气性也各不相同。有两种类型的层间的岩体,即弹性完整和塑料破碎岩体。流过程可分为三个阶段从受保护的缝到PSWF:(1)气体解吸发生通过提供保护煤层卸压,(2)气体通过完整的岩体渗流阶段发生,和3)气体通过裂隙岩体渗流阶段发生。为了简化开发模型,以下假设进行分析:gydF4y2Ba

整个模型系统有相同的恒定的温度gydF4y2Ba

水和其他气体之间的相互作用被忽视,即。,只有气体建模gydF4y2Ba

气体流沿着最短路径断裂的岩石gydF4y2Ba

天然气被视为理想气体;动态粘度和扩散系数视为常量值gydF4y2Ba

研究表明,有三个主要的几个州原位煤层气,即吸收气体,游离气和溶解气。的过程中,煤层的卸压和扩张,释放气体主要是吸收与煤层的形成。解吸过程是相反的吸收过程的机制,这是由朗缪尔方程(gydF4y2Ba 27gydF4y2Ba]。在气体解吸阶段,气体扩散到完整岩体的弹性岩石孔隙相对较小。扩散的过程遵循一个不稳定和非平衡吸附模型,称为菲克第二定律。gydF4y2Ba (7)gydF4y2Ba 问gydF4y2Ba =gydF4y2Ba DgydF4y2Ba ×gydF4y2Ba cgydF4y2Ba −gydF4y2Ba cgydF4y2Ba 0gydF4y2Ba HgydF4y2Ba ×gydF4y2Ba 一个gydF4y2Ba 。gydF4y2Ba

在此,gydF4y2Ba 问gydF4y2Ba 在单位时间内气体扩散的数量,表现在米gydF4y2Ba3gydF4y2Ba/ s;gydF4y2Ba DgydF4y2Ba 扩散系数;gydF4y2Ba cgydF4y2Ba −gydF4y2Ba cgydF4y2Ba 0gydF4y2Ba 浓度差;gydF4y2Ba HgydF4y2Ba 夹板的距离,这是表示在米;和gydF4y2Ba 一个gydF4y2Ba 是孔隙面积,m表示gydF4y2Ba2gydF4y2Ba。正常的气体是假定服从理想气体的状态方程:gydF4y2Ba (8)gydF4y2Ba PgydF4y2Ba VgydF4y2Ba =gydF4y2Ba NgydF4y2Ba RgydF4y2Ba TgydF4y2Ba 。gydF4y2Ba

在此,gydF4y2Ba PgydF4y2Ba 理想气体的压力,这是表达MPa;gydF4y2Ba VgydF4y2Ba 是理想气体体积,用mgydF4y2Ba3gydF4y2Ba;gydF4y2Ba NgydF4y2Ba 是气体的物质的量,这是摩尔表示;gydF4y2Ba TgydF4y2Ba 理想气体的热力学温度,用开尔文表达;gydF4y2Ba RgydF4y2Ba 理想气体常数。gydF4y2Ba

方程gydF4y2Ba 9gydF4y2Ba演示了气体的浓度之间的关系和衬底的气体量:gydF4y2Ba (9)gydF4y2Ba cgydF4y2Ba =gydF4y2Ba NgydF4y2Ba VgydF4y2Ba 。gydF4y2Ba

保护煤层的瓦斯压力表示gydF4y2Ba PgydF4y2Ba 1gydF4y2Ba ,距离表示gydF4y2Ba HgydF4y2Ba 。气体压力的完整岩体被认为等于大气压力(gydF4y2Ba PgydF4y2Ba 0gydF4y2Ba )。方程(gydF4y2Ba 10gydF4y2Ba)是用于计算气体量(gydF4y2Ba 问gydF4y2Ba 1gydF4y2Ba ),扩散到完整岩体:gydF4y2Ba (10)gydF4y2Ba 问gydF4y2Ba 1gydF4y2Ba =gydF4y2Ba DgydF4y2Ba ×gydF4y2Ba pgydF4y2Ba 1gydF4y2Ba −gydF4y2Ba pgydF4y2Ba 0gydF4y2Ba HgydF4y2Ba RgydF4y2Ba TgydF4y2Ba ×gydF4y2Ba 一个gydF4y2Ba 。gydF4y2Ba

当气体扩散到完整的岩体,它进入完好无损,裂隙岩体的渗流阶段。方程(gydF4y2Ba 11gydF4y2Ba在岩体)表明,气体渗流遵循达西定律。gydF4y2Ba (11)gydF4y2Ba 问gydF4y2Ba =gydF4y2Ba KgydF4y2Ba ×gydF4y2Ba PgydF4y2Ba −gydF4y2Ba PgydF4y2Ba 0gydF4y2Ba μgydF4y2Ba lgydF4y2Ba ×gydF4y2Ba 一个gydF4y2Ba 。gydF4y2Ba

在这里,gydF4y2Ba 问gydF4y2Ba 气体渗流的数量,这是表示在米gydF4y2Ba3gydF4y2Ba/ s;gydF4y2Ba KgydF4y2Ba 渗透率的岩石,用mgydF4y2Ba2−gydF4y2Ba;gydF4y2Ba PgydF4y2Ba −gydF4y2Ba PgydF4y2Ba 0gydF4y2Ba 是压差,在Pa表达;gydF4y2Ba lgydF4y2Ba 渗流的距离,这是表示在米;gydF4y2Ba 一个gydF4y2Ba 渗流面积,表现在米gydF4y2Ba2gydF4y2Ba;和gydF4y2Ba μgydF4y2Ba 动态粘度系数,表示在MPa·s。方程(gydF4y2Ba 12gydF4y2Ba)是用来计算气体压力后完整岩体的阶段。所示(gydF4y2Ba 12gydF4y2Ba的垂直长度),完整的岩体称为gydF4y2Ba lgydF4y2Ba 1gydF4y2Ba 。据流守恒原则,受影响的岩体的渗透性是众所周知的gydF4y2Ba kgydF4y2Ba 2gydF4y2Ba 。gydF4y2Ba (12)gydF4y2Ba pgydF4y2Ba zgydF4y2Ba =gydF4y2Ba pgydF4y2Ba 1gydF4y2Ba −gydF4y2Ba DgydF4y2Ba μgydF4y2Ba lgydF4y2Ba 1gydF4y2Ba pgydF4y2Ba 1gydF4y2Ba −gydF4y2Ba pgydF4y2Ba 0gydF4y2Ba kgydF4y2Ba 2gydF4y2Ba HgydF4y2Ba RgydF4y2Ba TgydF4y2Ba 。gydF4y2Ba

气体流(gydF4y2Ba 问gydF4y2Ba 2gydF4y2Ba )保护煤层裂隙岩体PSWF,利用计算(gydF4y2Ba 13gydF4y2Ba)。岩石渗透率叫做gydF4y2Ba kgydF4y2Ba 3gydF4y2Ba ,gydF4y2Ba lgydF4y2Ba 2gydF4y2Ba 被称为垂直长度,PSWF通风压力吗gydF4y2Ba pgydF4y2Ba 2gydF4y2Ba 。gydF4y2Ba (13)gydF4y2Ba 问gydF4y2Ba 2gydF4y2Ba =gydF4y2Ba pgydF4y2Ba 1gydF4y2Ba −gydF4y2Ba pgydF4y2Ba 2gydF4y2Ba μgydF4y2Ba lgydF4y2Ba 2gydF4y2Ba −gydF4y2Ba DgydF4y2Ba μgydF4y2Ba lgydF4y2Ba 1gydF4y2Ba pgydF4y2Ba 1gydF4y2Ba −gydF4y2Ba pgydF4y2Ba 0gydF4y2Ba kgydF4y2Ba 2gydF4y2Ba HgydF4y2Ba RgydF4y2Ba TgydF4y2Ba kgydF4y2Ba 3gydF4y2Ba 一个gydF4y2Ba 。gydF4y2Ba

PCSM后,(gydF4y2Ba 14gydF4y2Ba)是用于计算最大和最小原则强调对裂隙岩体在任何时候(gydF4y2Ba 28gydF4y2Ba]。gydF4y2Ba (14)gydF4y2Ba σgydF4y2Ba 1gydF4y2Ba ,gydF4y2Ba σgydF4y2Ba 3gydF4y2Ba =gydF4y2Ba 问gydF4y2Ba 1gydF4y2Ba πgydF4y2Ba θgydF4y2Ba ±gydF4y2Ba 1gydF4y2Ba πgydF4y2Ba 问gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba +gydF4y2Ba 问gydF4y2Ba 2gydF4y2Ba 2gydF4y2Ba 罪gydF4y2Ba θgydF4y2Ba +gydF4y2Ba 问gydF4y2Ba 2gydF4y2Ba θgydF4y2Ba +gydF4y2Ba γgydF4y2Ba lgydF4y2Ba 因为gydF4y2Ba αgydF4y2Ba 。gydF4y2Ba

在此,gydF4y2Ba 问gydF4y2Ba 1gydF4y2Ba 的垂直压力作用于地板上,这是表达MPa;gydF4y2Ba 问gydF4y2Ba 2gydF4y2Ba 是水平应力,它是表达MPa;和gydF4y2Ba αgydF4y2Ba 煤层的角度,这是表达的度。如图gydF4y2Ba 2gydF4y2Ba,gydF4y2Ba θgydF4y2Ba =gydF4y2Ba θgydF4y2Ba 1gydF4y2Ba −gydF4y2Ba θgydF4y2Ba 2gydF4y2Ba ,在那里gydF4y2Ba θgydF4y2Ba 1gydF4y2Ba 和gydF4y2Ba θgydF4y2Ba 2gydF4y2Ba 角点之间的连接和地板上,和他们也表示度;gydF4y2Ba lgydF4y2Ba 点之间的垂直距离是M和保护煤层的地板,这是表示在米;和gydF4y2Ba γgydF4y2Ba 岩石密度,MN / m表示gydF4y2Ba3gydF4y2Ba。gydF4y2Ba

煤层气体扩散和渗流。gydF4y2Ba

permeability-stress模型被选中来决定完好无损,裂隙岩体的渗透性。通过使用实验室数据安装(gydF4y2Ba 3gydF4y2Ba)和(gydF4y2Ba 4gydF4y2Ba),最短路径选择,即,gydF4y2Ba θgydF4y2Ba =gydF4y2Ba 0gydF4y2Ba 。假设完整岩体的应力gydF4y2Ba γgydF4y2Ba lgydF4y2Ba ,原岩体应力gydF4y2Ba γgydF4y2Ba hgydF4y2Ba ,gydF4y2Ba hgydF4y2Ba 是保护煤层的埋藏深度。这种假设是基于在实际生产过程中,大多数气体释放保护缝不扩散到PSWF。通过不断推进工作面,保护煤层的气体压力逐渐减小。方程(gydF4y2Ba 15gydF4y2Ba)必须乘以系数的丧失gydF4y2Ba δgydF4y2Ba 。用参数为(gydF4y2Ba 13gydF4y2Ba),(gydF4y2Ba 15gydF4y2Ba)方法如下:gydF4y2Ba (15)gydF4y2Ba 问gydF4y2Ba 2gydF4y2Ba =gydF4y2Ba 一个gydF4y2Ba δgydF4y2Ba pgydF4y2Ba 1gydF4y2Ba −gydF4y2Ba DgydF4y2Ba μgydF4y2Ba lgydF4y2Ba 1gydF4y2Ba pgydF4y2Ba 1gydF4y2Ba −gydF4y2Ba pgydF4y2Ba 0gydF4y2Ba /gydF4y2Ba kgydF4y2Ba 0gydF4y2Ba HgydF4y2Ba RgydF4y2Ba TgydF4y2Ba ∫gydF4y2Ba 0gydF4y2Ba lgydF4y2Ba 1gydF4y2Ba egydF4y2Ba −gydF4y2Ba 0.25gydF4y2Ba γgydF4y2Ba lgydF4y2Ba 1gydF4y2Ba +gydF4y2Ba lgydF4y2Ba 2gydF4y2Ba −gydF4y2Ba hgydF4y2Ba −gydF4y2Ba xgydF4y2Ba dgydF4y2Ba xgydF4y2Ba −gydF4y2Ba pgydF4y2Ba 2gydF4y2Ba ∫gydF4y2Ba 0gydF4y2Ba lgydF4y2Ba 2gydF4y2Ba kgydF4y2Ba 1gydF4y2Ba γgydF4y2Ba xgydF4y2Ba /gydF4y2Ba 因为gydF4y2Ba αgydF4y2Ba 米gydF4y2Ba dgydF4y2Ba xgydF4y2Ba μgydF4y2Ba lgydF4y2Ba 2gydF4y2Ba =gydF4y2Ba pgydF4y2Ba 1gydF4y2Ba −gydF4y2Ba DgydF4y2Ba γgydF4y2Ba μgydF4y2Ba lgydF4y2Ba 1gydF4y2Ba pgydF4y2Ba 1gydF4y2Ba −gydF4y2Ba pgydF4y2Ba 0gydF4y2Ba /gydF4y2Ba 4gydF4y2Ba kgydF4y2Ba 0gydF4y2Ba HgydF4y2Ba RgydF4y2Ba TgydF4y2Ba egydF4y2Ba −gydF4y2Ba 0.25gydF4y2Ba γgydF4y2Ba lgydF4y2Ba 2gydF4y2Ba −gydF4y2Ba hgydF4y2Ba −gydF4y2Ba egydF4y2Ba −gydF4y2Ba 0.25gydF4y2Ba γgydF4y2Ba lgydF4y2Ba 1gydF4y2Ba +gydF4y2Ba lgydF4y2Ba 2gydF4y2Ba −gydF4y2Ba hgydF4y2Ba −gydF4y2Ba pgydF4y2Ba 2gydF4y2Ba 一个gydF4y2Ba δgydF4y2Ba kgydF4y2Ba 1gydF4y2Ba γgydF4y2Ba 米gydF4y2Ba lgydF4y2Ba 2gydF4y2Ba 米gydF4y2Ba 米gydF4y2Ba +gydF4y2Ba 1gydF4y2Ba μgydF4y2Ba 因为gydF4y2Ba 米gydF4y2Ba αgydF4y2Ba 。gydF4y2Ba

图gydF4y2Ba 2gydF4y2Ba总结了整个上层PCSM的GSM。gydF4y2Ba

方程(gydF4y2Ba 15gydF4y2Ba)表明,有三个主要因素相关的气体量,这可能会导致气体流入PSWF从受保护的缝,即gydF4y2Ba

seam的初始气体压力保护gydF4y2Ba

夹层岩体层间的间距,支离破碎的百分比和完整的岩体,骨折和完整岩体的渗透系数gydF4y2Ba

PSWF通风压力和渗流区域的地板上gydF4y2Ba

4所示。地板上失败的深度分析gydF4y2Ba 4.1。理论计算模型层裂缝深度发展gydF4y2Ba

分析的最重要的部分是确定地板失败深度根据GSM保护焊缝,在引入部分gydF4y2Ba 3gydF4y2Ba。计算地板的方法失败的深度由理论计算和经验回归方程。gydF4y2Ba

煤炭开采过程中过重的负荷转移到煤工作面前的身体紧凑的下降在采空区岩石形成支承压力峰值区域。在工作面前,地板岩石在压缩状态,经历压缩位移。在背后的采空区工作面,地板岩石在减压区,它经历了膨胀位移。岩石剪切断裂通常发生在地板由于剪切变形,这是位于压缩和扩张的边境地区。传统的工作面推进过程,在地板上有三个阶段的压力:压力升级之前采矿、应力降开采后,和压力恢复。这些阶段的应力反复发生与工作面推进。在周围的地板区域工作面,岩体应力分布在地板上被开采的影响。它分为四个水平区域:(1)原岩应力区,(2)压缩区域,(3)扩张地区,和(4)应力恢复区。地板层的分布分为以下四个垂直区:(i)破碎区,(2)新受损区,(3)最初受损区域,(iv)断裂区(gydF4y2Ba 29日gydF4y2Ba]。图gydF4y2Ba 3gydF4y2Ba显示了地板的地区,受支承压力的影响。gydF4y2Ba

地板上失败的深度由支承压力引起的。gydF4y2Ba

针对摩尔-库仑强度理论,根据弹塑性力学和塑性区成立于地面岩石塑性变形。在一定范围的工作面楼当轴承压力达到或超过极限值的压力。gydF4y2Ba

支承压力达到最大容量时,岩石的一部分失败完全通过加入塑料区岩体的支承压力的周边地区。与岩体的塑性变形,形成一个连续的滑动面采空区方向。gydF4y2Ba

地板的最大破坏深度,这是由于应力集中在轴承极限压力条件下的岩石上,可以通过使用(gydF4y2Ba 16gydF4y2Ba)。地层的极限承载力是针对摩尔-库仑失败准则计算(gydF4y2Ba 30.gydF4y2Ba]。gydF4y2Ba

塑料的宽度面积计算使用(gydF4y2Ba 16gydF4y2Ba)。方程(gydF4y2Ba 17gydF4y2Ba)是用于确定地板的最大破坏深度。gydF4y2Ba (16)gydF4y2Ba lgydF4y2Ba ′gydF4y2Ba =gydF4y2Ba 一个gydF4y2Ba ×gydF4y2Ba 1gydF4y2Ba 1gydF4y2Ba −gydF4y2Ba ngydF4y2Ba 米gydF4y2Ba γgydF4y2Ba HgydF4y2Ba /gydF4y2Ba 2gydF4y2Ba σgydF4y2Ba 米gydF4y2Ba +gydF4y2Ba ngydF4y2Ba 米gydF4y2Ba γgydF4y2Ba HgydF4y2Ba 2gydF4y2Ba −gydF4y2Ba 1gydF4y2Ba 或gydF4y2Ba lgydF4y2Ba ′gydF4y2Ba =gydF4y2Ba 0.015gydF4y2Ba HgydF4y2Ba ,gydF4y2Ba (17)gydF4y2Ba hgydF4y2Ba =gydF4y2Ba lgydF4y2Ba ′gydF4y2Ba 因为gydF4y2Ba βgydF4y2Ba 2gydF4y2Ba 因为gydF4y2Ba πgydF4y2Ba /gydF4y2Ba 4gydF4y2Ba +gydF4y2Ba βgydF4y2Ba /gydF4y2Ba 2gydF4y2Ba ×gydF4y2Ba egydF4y2Ba πgydF4y2Ba /gydF4y2Ba 4gydF4y2Ba +gydF4y2Ba βgydF4y2Ba /gydF4y2Ba 2gydF4y2Ba ×gydF4y2Ba 棕褐色gydF4y2Ba βgydF4y2Ba 。gydF4y2Ba

在这里,gydF4y2Ba 一个gydF4y2Ba 表示米,是工作面顶板控制距离的两倍;gydF4y2Ba ngydF4y2Ba 米gydF4y2Ba 的最大应力集中系数;gydF4y2Ba σgydF4y2Ba 米gydF4y2Ba 煤炭的极限应力的身体优势,这是表达MPa;gydF4y2Ba HgydF4y2Ba 是垂直距离工作面地面,表现在米;和gydF4y2Ba βgydF4y2Ba 是地板层的内摩擦角,它是表达度。gydF4y2Ba

方程(gydF4y2Ba 18gydF4y2Ba失败的)用于理论计算深度的地板gydF4y2Ba 31日gydF4y2Ba]。gydF4y2Ba (18)gydF4y2Ba hgydF4y2Ba =gydF4y2Ba 1.57gydF4y2Ba γgydF4y2Ba 2gydF4y2Ba HgydF4y2Ba 2gydF4y2Ba lgydF4y2Ba 4gydF4y2Ba RgydF4y2Ba cgydF4y2Ba 2gydF4y2Ba 。gydF4y2Ba

在这里,gydF4y2Ba hgydF4y2Ba 是破碎区和深度gydF4y2Ba lgydF4y2Ba 工作面长度,这两种米表示,然后呢gydF4y2Ba RgydF4y2Ba cgydF4y2Ba 岩石的单轴抗压强度,是MPa表示。gydF4y2Ba

在以前的研究中,进行统计分析,以确定故障的深度地板各煤矿在煤炭开采过程中。相应的线性回归方程是由使用多元线性和非线性统计分析方法。一些典型的方程如下gydF4y2Ba 31日gydF4y2Ba- - - - - -gydF4y2Ba 33gydF4y2Ba]:gydF4y2Ba (19)gydF4y2Ba hgydF4y2Ba =gydF4y2Ba 0.0113gydF4y2Ba HgydF4y2Ba +gydF4y2Ba 6.25gydF4y2Ba lngydF4y2Ba lgydF4y2Ba 40gydF4y2Ba +gydF4y2Ba 2.52gydF4y2Ba lngydF4y2Ba 米gydF4y2Ba hgydF4y2Ba 1.48gydF4y2Ba ,gydF4y2Ba (20)gydF4y2Ba hgydF4y2Ba =gydF4y2Ba 0.0111gydF4y2Ba lgydF4y2Ba +gydF4y2Ba 0.006gydF4y2Ba HgydF4y2Ba +gydF4y2Ba 4.541gydF4y2Ba fgydF4y2Ba cgydF4y2Ba −gydF4y2Ba 0.009gydF4y2Ba αgydF4y2Ba −gydF4y2Ba 2.4gydF4y2Ba ,gydF4y2Ba (21)gydF4y2Ba hgydF4y2Ba =gydF4y2Ba −gydF4y2Ba 2.0234gydF4y2Ba +gydF4y2Ba 1.48gydF4y2Ba ×gydF4y2Ba 10gydF4y2Ba −gydF4y2Ba 26gydF4y2Ba HgydF4y2Ba 9gydF4y2Ba +gydF4y2Ba 0.1016gydF4y2Ba lgydF4y2Ba +gydF4y2Ba 0.1913gydF4y2Ba αgydF4y2Ba +gydF4y2Ba 1.0637gydF4y2Ba 米gydF4y2Ba hgydF4y2Ba −gydF4y2Ba 5.5536gydF4y2Ba DgydF4y2Ba +gydF4y2Ba 7.5070gydF4y2Ba 我gydF4y2Ba 。gydF4y2Ba

在这里,gydF4y2Ba 米gydF4y2Ba hgydF4y2Ba 煤层的开采高度,表达在米;gydF4y2Ba DgydF4y2Ba 地板的失败的抵抗能力;和gydF4y2Ba 我gydF4y2Ba 的相关因素存在断层破碎带或失败情况。利用上述方程估算地板失败深度,完整的长度和裂隙岩体在GSM。gydF4y2Ba

对于一个案例研究在城市某煤矿,山西,中国,验证程序和相关计算下面描述。gydF4y2Ba

实际的两个PSWFs用于验证模型。沿走向的距离是1142米,180米的PSWF4216斜宽度。PSWF4215,沿走向的距离是1158米和200米的宽度。这两个工作面临位于2号煤层;煤层的厚度约为0 - 2.5,平均厚度为1米。2号煤层的深度是430米。间距的2号3号煤层煤层和底层(受保护的缝)的范围9.8 - -20.7米,平均间距为15.5米。这个保护煤层的瓦斯含量不高,和0.55 MPa的压力。低了3号煤层突出危险。当3号煤层的相对气体排放上升到27.82米gydF4y2Ba3gydF4y2Ba/ t,气体压力成为2.03 MPa和低渗透煤层。表gydF4y2Ba 1gydF4y2Ba显示了该工作面围岩岩体的力学参数。图gydF4y2Ba 4gydF4y2Ba显示了工作面钻孔图。gydF4y2Ba

岩石的物理力学参数和煤炭。gydF4y2Ba

序列号gydF4y2Ba 岩性gydF4y2Ba 厚度(m)gydF4y2Ba 密度(公斤·mgydF4y2Ba−3gydF4y2Ba)gydF4y2Ba 体积弹性模量(GPa)gydF4y2Ba 剪切模量(GPa)gydF4y2Ba 凝聚力(MPa)gydF4y2Ba 内摩擦角(°)gydF4y2Ba 扩展强度(MPa)gydF4y2Ba
1gydF4y2Ba 细砂岩gydF4y2Ba 8.01gydF4y2Ba 2500年gydF4y2Ba 3所示。0gydF4y2Ba 1。5gydF4y2Ba 2.0gydF4y2Ba 35gydF4y2Ba 1。5gydF4y2Ba
2gydF4y2Ba 粉砂岩gydF4y2Ba 2.99gydF4y2Ba 2400年gydF4y2Ba 2.6gydF4y2Ba 1。6gydF4y2Ba 2.2gydF4y2Ba 36gydF4y2Ba 1。5gydF4y2Ba
3gydF4y2Ba 煤炭行gydF4y2Ba 0.49gydF4y2Ba 1450年gydF4y2Ba 1。8gydF4y2Ba 1。0gydF4y2Ba 1。5gydF4y2Ba 30.gydF4y2Ba 1。0gydF4y2Ba
4gydF4y2Ba 粉砂岩gydF4y2Ba 3.67gydF4y2Ba 2500年gydF4y2Ba 2.6gydF4y2Ba 1。6gydF4y2Ba 2.2gydF4y2Ba 36gydF4y2Ba 1。5gydF4y2Ba
5gydF4y2Ba 细砂岩gydF4y2Ba 4.23gydF4y2Ba 2500年gydF4y2Ba 2.9gydF4y2Ba 1。9gydF4y2Ba 2.8gydF4y2Ba 39gydF4y2Ba 1。8gydF4y2Ba
6gydF4y2Ba 2 #煤层gydF4y2Ba 1.00gydF4y2Ba 1450年gydF4y2Ba 1。8gydF4y2Ba 1。0gydF4y2Ba 1。5gydF4y2Ba 30.gydF4y2Ba 1。0gydF4y2Ba
7gydF4y2Ba 砂质泥岩gydF4y2Ba 2.58gydF4y2Ba 2300年gydF4y2Ba 2.0gydF4y2Ba 1。3gydF4y2Ba 2.0gydF4y2Ba 38gydF4y2Ba 1。5gydF4y2Ba
8gydF4y2Ba 粉砂岩gydF4y2Ba 4.92gydF4y2Ba 2400年gydF4y2Ba 2.6gydF4y2Ba 1。6gydF4y2Ba 2.2gydF4y2Ba 36gydF4y2Ba 1。5gydF4y2Ba
9gydF4y2Ba Middle-fine砂岩gydF4y2Ba 8.11gydF4y2Ba 2600年gydF4y2Ba 4.0gydF4y2Ba 2.4gydF4y2Ba 2.8gydF4y2Ba 38gydF4y2Ba 2.8gydF4y2Ba
10gydF4y2Ba 粉砂岩gydF4y2Ba 1.98gydF4y2Ba 2400年gydF4y2Ba 2.6gydF4y2Ba 1。6gydF4y2Ba 2.2gydF4y2Ba 36gydF4y2Ba 1。5gydF4y2Ba
11gydF4y2Ba 3 #煤层gydF4y2Ba 6.09gydF4y2Ba 1400年gydF4y2Ba 1。8gydF4y2Ba 1gydF4y2Ba 1。5gydF4y2Ba 30.gydF4y2Ba 1。0gydF4y2Ba
12gydF4y2Ba 砂质泥岩gydF4y2Ba 2.97gydF4y2Ba 2400年gydF4y2Ba 2.0gydF4y2Ba 1。3gydF4y2Ba 1。5gydF4y2Ba 30.gydF4y2Ba 1。2gydF4y2Ba
13gydF4y2Ba 细砂岩gydF4y2Ba 5.03gydF4y2Ba 2500年gydF4y2Ba 2.8gydF4y2Ba 1。8gydF4y2Ba 2.9gydF4y2Ba 36gydF4y2Ba 1。7gydF4y2Ba
14gydF4y2Ba Middle-fine砂岩gydF4y2Ba 45.64gydF4y2Ba 2500年gydF4y2Ba 3所示。0gydF4y2Ba 1。5gydF4y2Ba 2.0gydF4y2Ba 35gydF4y2Ba 1。5gydF4y2Ba

工作面钻孔图。gydF4y2Ba

为了验证的准确性发达GSM保护焊缝,地板上失败的深度挖掘过程中首先确定工作面。使用理论方程(gydF4y2Ba 16gydF4y2Ba)- (gydF4y2Ba 18gydF4y2Ba从矿业(获得)和数据gydF4y2Ba HgydF4y2Ba 是430,gydF4y2Ba βgydF4y2Ba 是40°,gydF4y2Ba lgydF4y2Ba 是180,gydF4y2Ba γgydF4y2Ba 0.025公斤/厘米吗gydF4y2Ba3gydF4y2Ba,gydF4y2Ba RgydF4y2Ba cgydF4y2Ba 21.59 MPa),失败楼深处被估计为15.2米和17.52米。利用回归方程gydF4y2Ba 19gydF4y2Ba- - - - - -gydF4y2Ba 21gydF4y2Ba和获得的值(gydF4y2Ba 米gydF4y2Ba hgydF4y2Ba 是1.2米,gydF4y2Ba αgydF4y2Ba 是0°,gydF4y2Ba DgydF4y2Ba 是0.2,gydF4y2Ba 我gydF4y2Ba 是0,gydF4y2Ba fgydF4y2Ba 3),地板失败深度估计为13.73米,15.80米和16.44米。通过考虑理论和回归方程的结果,计算层深度失败被发现在13.73到17.52米,平均深度15.625失败。这个值被发现非常接近实际的平均层间间距。进一步确定失败的地板上,通过使用模拟数值分析进行了数值模拟软件。gydF4y2Ba

4.2。数值模拟地板裂缝发展深度gydF4y2Ba

莫尔-库仑本构模型被用来模拟,提出了岩石力学参数表gydF4y2Ba 1gydF4y2Ba。提出了一种模型的示意图如图gydF4y2Ba 5gydF4y2Ba。模型的长度是300米,它有一个宽80米。保护煤层的屋顶用JSET被划分为一个分层结构。结构破损没有观察到在地板上地层之间的防护和保护,是类似于保护煤层的屋顶。地层和岩石参数有不同的厚度。如图gydF4y2Ba 5gydF4y2Ba、地板层被分成一个多边形与全等方面利用2 d泰森多边形法多边形生成器。在这个数值模型,在模拟变形和刚性的多边形。然而,单个多边形不能打破。,一个llfractures must follow polygonal boundaries. The cemented contacts between two polygons can rupture through shear or tension depending on the stress state and properties of the contact surface. Therefore, fracture development can be realistically simulated through the contact rupture initiation, growth, and nucleation. The force-displacement relationship for a contact is determined by normal ( KgydF4y2Ba ngydF4y2Ba )和剪切(gydF4y2Ba KgydF4y2Ba 年代gydF4y2Ba )的接触刚度。接触性质、凝聚力、抗拉强度和摩擦角定义接触强度。骨折时创建界面的应力水平超过一个阈值在紧张或剪切。接触本构模型如图gydF4y2Ba 6gydF4y2Ba。gydF4y2Ba

原理图的数值模型。gydF4y2Ba

本构模型的接触模拟:(一)块之间的正常和剪切刚度;(b)本构行为的剪切和紧张(gydF4y2Ba 我gydF4y2Ba =gydF4y2Ba 年代gydF4y2Ba ,gydF4y2Ba ngydF4y2Ba)[gydF4y2Ba 34gydF4y2Ba]。gydF4y2Ba

PSWF开始挖掘50米距离左边边缘模型,挖掘5米100米距离的每一步。图gydF4y2Ba 7gydF4y2Ba说明的工作面推进地板裂缝发展逐步深入而裂缝发展的范围不断扩大。当PSWF达到30米的推进距离,防护和保护之间的岩体接缝连接骨折但骨折是稀疏的发展。开发变得更加密集推进距离时35米。推进40米的距离,受保护的缝了骨折。这导致气体使解除吸附seam和流入PSWF保护。gydF4y2Ba

地板裂缝深度图的不同推进距离。gydF4y2Ba

基于数值模拟的结果和理论计算,推断,地板断裂发生在seam PCSM过程中保护。因此,可以认为,没有完整的岩体之间的防护和保护。这种假设是审议在接下来的计算。gydF4y2Ba

5。气体渗流模型验证gydF4y2Ba

最具有挑战性的物理测量的参数在GSM,虽然一些参数在相同地质条件下保持不变。然而,它可以从上述分析得出结论,层间的岩体裂隙岩体。方程(gydF4y2Ba 15gydF4y2Ba)可以简化为gydF4y2Ba (22)gydF4y2Ba 问gydF4y2Ba 2gydF4y2Ba =gydF4y2Ba pgydF4y2Ba 1gydF4y2Ba −gydF4y2Ba pgydF4y2Ba 2gydF4y2Ba 一个gydF4y2Ba δgydF4y2Ba kgydF4y2Ba 1gydF4y2Ba γgydF4y2Ba lgydF4y2Ba 2gydF4y2Ba 米gydF4y2Ba 米gydF4y2Ba +gydF4y2Ba 1gydF4y2Ba μgydF4y2Ba 。gydF4y2Ba

根据实际数据考虑PSWF4216在汉城,中国,与天然气相关的参数如下:gydF4y2Ba PgydF4y2Ba 1gydF4y2Ba =gydF4y2Ba 2.03gydF4y2Ba MPa(保护煤层的气体压力),gydF4y2Ba PgydF4y2Ba 2gydF4y2Ba =gydF4y2Ba 0.1gydF4y2Ba MPa (PSWF的气体压力),gydF4y2Ba γgydF4y2Ba =gydF4y2Ba 0.025gydF4y2Ba 米gydF4y2Ba NgydF4y2Ba /gydF4y2Ba 米gydF4y2Ba 3gydF4y2Ba 和面积(斜宽度PSWF×地板宽度)的工作面是180米gydF4y2Ba2gydF4y2Ba和200米gydF4y2Ba2gydF4y2Ba分别为PSWF4216 PSWF4215。的价值gydF4y2Ba μgydF4y2Ba 天然气是1.75×10gydF4y2Ba−11gydF4y2BaMPa·s。的gydF4y2Ba 米gydF4y2Ba 和gydF4y2Ba kgydF4y2Ba 1gydF4y2Ba 被分配为−0.8616和2.613×10gydF4y2Ba−13gydF4y2Ba米gydF4y2Ba2gydF4y2Ba分别介绍了在文献[gydF4y2Ba 22gydF4y2Ba]。用上面的数据(gydF4y2Ba 22gydF4y2Ba),保护脸的气体排放计算通过使用gydF4y2Ba (23)gydF4y2Ba 问gydF4y2Ba ggydF4y2Ba 米gydF4y2Ba 3gydF4y2Ba /gydF4y2Ba 最小值gydF4y2Ba =gydF4y2Ba 60gydF4y2Ba 问gydF4y2Ba 2gydF4y2Ba =gydF4y2Ba 53986.01gydF4y2Ba ×gydF4y2Ba δgydF4y2Ba lgydF4y2Ba −gydF4y2Ba 0.8616gydF4y2Ba 。gydF4y2Ba

保护和保护之间的接缝,改变值与层间间距得到的瓦斯抽放钻孔。这是放置在地板上的空气返回网关。在工作面推进距离30米,地板骨折PSWF可以发展成缝的保护。在工作面通风空气中甲烷的数据也被后进步30米的距离。如图gydF4y2Ba 8gydF4y2Ba、散点图和拟合曲线的通风空气中甲烷与层间间距的PSWF4216获得相应的层间间距数据。gydF4y2Ba

散点图的通风空气对PSWF4216甲烷。gydF4y2Ba

损失系数(gydF4y2Ba δgydF4y2Ba )计算通过使用通风空气中甲烷的散点图(图gydF4y2Ba 8gydF4y2Ba)。相关系数(gydF4y2Ba RgydF4y2Ba 2gydF4y2Ba )拟合方程的0.9022,这显示了一个可接受的工程实测数据的拟合情况。这些发现证实了(gydF4y2Ba 22gydF4y2Ba)是用于计算保护煤层气体排放数量,对应与现场测量精确。此外,损失系数(gydF4y2Ba δgydF4y2Ba )为1.012×10gydF4y2Ba−3gydF4y2Ba获得了。为了进一步验证的正确性GSM,损失系数(gydF4y2Ba δgydF4y2Ba )为1.012×10gydF4y2Ba−3gydF4y2Ba和参数的PSWF4215(只有PSWF4215的面积是不同于PSWF4216)。被替换成(gydF4y2Ba 22gydF4y2Ba)。图gydF4y2Ba 9gydF4y2Ba显示了理论曲线,它是根据(gydF4y2Ba 22gydF4y2BaPSWF4215)。图gydF4y2Ba 9gydF4y2Ba显示测量数据与理论曲线相匹配,这也进一步验证了GSM的可行性和有效性。gydF4y2Ba

测量数据和PSWF4215的理论曲线。gydF4y2Ba

6。讨论gydF4y2Ba

深煤层开采过程中,PCSM已经成为一个最重要的地区爆发的措施缓解固有的爆发风险很高。在PCSM过程中,甲烷浓度超过极限的现象通常发生在工作面。是非常重要的理解PSWF的气体排放影响因素和开发一个方程计算其数量。在这项研究中,一个GSM开发上层PCSM的过程。使用GSM的方程、气体的数量而眠保护缝到保护煤层估计。菲克第二定律的扩散和达西流定律被用于估计。此外,渗透率和压力之间的关系研究处于弹塑性状态。围岩的力学也被调查。有四个假设和简化几个用于GSM的建立;它极大地简化了计算复杂度,减少了选择参数的工程计算。 As shown in ( 15gydF4y2Ba),它更简单,更容易使用GSM。该模型不需要花很多的时间和金钱为额外的实地测试,并没有额外的GSM的不确定性因素。这个模型对于工程应用,具有独特的优势和效率高。然而,温度、气体成分、含水率因素假设和简化管理。他们也影响煤和岩石的渗透率,降低了模型的精度。在未来的工作中,上述因素被认为是在GSM。gydF4y2Ba

7所示。结论gydF4y2Ba

PCSM的过程中,围岩渗透性分为两种类型,即失败塑料和弹性岩体。岩石渗透率和压力之间的关系也总结下这两种演化类型,GSM之间的防护和保护提供理论依据。最终,从保护煤层气体排放分类分三个阶段,即(1)气体从保护煤层卸压解吸阶段,(2)气体渗流阶段完整的岩体,和(3)气体在裂隙岩体渗流阶段(图gydF4y2Ba 10gydF4y2Ba)。从这些方程,可以得出的初始气体压力保护焊缝,夹层岩石的特点,和通风压力的工作面瓦斯涌出量的主要影响因素,从保护煤层PSWF流淌。gydF4y2Ba

气体渗流模型的计算过程。gydF4y2Ba

保护煤层的实际地质条件研究了煤矿,位于汉城,中国。GSM是用来计算从保护煤层气体排放数量。建模结果精确对应的拟合曲线与现场测量通风空气PSWF4216甲烷数量和PSWF4215层间间距和损耗系数(gydF4y2Ba δgydF4y2Ba )为1.012×10gydF4y2Ba−3gydF4y2Ba。值得注意的是,实际情况提出了简化的模型。这些条件是基于上述假设。在未来的研究中,从采空区气体排放及其影响因素确定。此外,一些这些方程中参数的调整分析的基础上额外的字段数据测量的工作面PCSM的过程。未来的研究将帮助我们进一步降低气体积累风险PCSM的过程。gydF4y2Ba

数据可用性gydF4y2Ba

使用的数据来支持本研究的发现可以从相应的作者。gydF4y2Ba

的利益冲突gydF4y2Ba

作者宣称没有利益冲突。gydF4y2Ba

确认gydF4y2Ba

金融支持这项工作得到了北京自然科学基金(8184082);中国国家自然科学基金(51704274号,51502339,,51704278);一流的本科院校的专业建设贵州省(SJZY2017006);国家重点实验室开放项目的煤炭资源安全开采,CUMT (SKLCRSM16KF 06);贵州教育的年轻科学家和技术人员发展项目部门(钱他肯塔基州子[2018]115号);和曰气杰出学者项目,中国矿业大学科技,北京。gydF4y2Ba

刘gydF4y2Ba y W。gydF4y2Ba 李gydF4y2Ba g F。gydF4y2Ba 可靠性研究保护层开采和减压瓦斯抽放gydF4y2Ba 《采矿与安全工程gydF4y2Ba 2013年gydF4y2Ba 3gydF4y2Ba 426年gydF4y2Ba 431年gydF4y2Ba 刘gydF4y2Ba 一个。gydF4y2Ba 你gydF4y2Ba 年代。gydF4y2Ba 王gydF4y2Ba F。gydF4y2Ba 元gydF4y2Ba Y。gydF4y2Ba 张gydF4y2Ba C。gydF4y2Ba 张gydF4y2Ba Y。gydF4y2Ba 数值模拟的上部和下部的卸压规则保护烟煤thin-protective-seam矿业gydF4y2Ba 灾难的进步gydF4y2Ba 2013年gydF4y2Ba 6gydF4y2Ba 16gydF4y2Ba 25gydF4y2Ba 傅gydF4y2Ba j . H。gydF4y2Ba 程gydF4y2Ba y . P。gydF4y2Ba 中国的煤与瓦斯突出情况和控制对策gydF4y2Ba 《采矿与安全工程gydF4y2Ba 2007年gydF4y2Ba 24gydF4y2Ba 3gydF4y2Ba 253年gydF4y2Ba 259年gydF4y2Ba 曹国伟gydF4y2Ba w·A。gydF4y2Ba 李大钊gydF4y2Ba s . O。gydF4y2Ba Xue-shenggydF4y2Ba d . U。gydF4y2Ba 煤与瓦斯突出的预测基于距离判别分析方法及其应用gydF4y2Ba 《采矿与安全工程gydF4y2Ba 2009年gydF4y2Ba 4gydF4y2Ba 470年gydF4y2Ba 447年gydF4y2Ba 汉gydF4y2Ba J。gydF4y2Ba 张gydF4y2Ba h·W。gydF4y2Ba 首歌gydF4y2Ba w·H。gydF4y2Ba 王gydF4y2Ba Z。gydF4y2Ba 地应力矿区煤与瓦斯突出的领域gydF4y2Ba 中国岩石力学与工程学报gydF4y2Ba 2008年gydF4y2Ba 2gydF4y2Ba 3852年gydF4y2Ba 3859年gydF4y2Ba 张gydF4y2Ba C。gydF4y2Ba 你gydF4y2Ba 年代。gydF4y2Ba 张gydF4y2Ba lgydF4y2Ba 程ydF4y2Ba M。gydF4y2Ba 研究渗流方向对渗透率的影响压力测试gydF4y2Ba 阿拉伯科学与工程》杂志上gydF4y2Ba 2016年gydF4y2Ba 41gydF4y2Ba 11gydF4y2Ba 4583年gydF4y2Ba 4596年gydF4y2Ba 10.1007 / s13369 - 016 - 2215 - 2gydF4y2Ba 2 - s2.0 - 84991608360gydF4y2Ba 张gydF4y2Ba C。gydF4y2Ba 你gydF4y2Ba 年代。gydF4y2Ba 白gydF4y2Ba Q。gydF4y2Ba 杨gydF4y2Ba G。gydF4y2Ba 张gydF4y2Ba lgydF4y2Ba 评估减压开采指标基于表面气体通气孔在煤矿长壁开采数据gydF4y2Ba 天然气的科学与工程》杂志上gydF4y2Ba 2015年gydF4y2Ba 24gydF4y2Ba 431年gydF4y2Ba 440年gydF4y2Ba 10.1016 / j.jngse.2015.04.012gydF4y2Ba 2 - s2.0 - 84927770159gydF4y2Ba 张gydF4y2Ba C。gydF4y2Ba 你gydF4y2Ba 年代。gydF4y2Ba 张gydF4y2Ba lgydF4y2Ba 白gydF4y2Ba Q。gydF4y2Ba 元gydF4y2Ba Y。gydF4y2Ba 王gydF4y2Ba F。gydF4y2Ba 的方法确定采空区渗透率的演化规律及其在煤矿长壁分布gydF4y2Ba 地球物理与工程杂志》上gydF4y2Ba 2016年gydF4y2Ba 13gydF4y2Ba 2gydF4y2Ba 181年gydF4y2Ba 193年gydF4y2Ba 10.1088 / 1742 - 2132/13/2/181gydF4y2Ba 2 - s2.0 - 84962787526gydF4y2Ba 王gydF4y2Ba F。gydF4y2Ba 任gydF4y2Ba T。gydF4y2Ba 你gydF4y2Ba 年代。gydF4y2Ba 亨格福特gydF4y2Ba F。gydF4y2Ba 阿齐兹gydF4y2Ba N。gydF4y2Ba 实现温室气体减排的地下深孔定向钻进技术在中国煤矿gydF4y2Ba 国际期刊的温室气体控制gydF4y2Ba 2012年gydF4y2Ba 11gydF4y2Ba 6gydF4y2Ba 290年gydF4y2Ba 303年gydF4y2Ba 10.1016 / j.ijggc.2012.09.006gydF4y2Ba 2 - s2.0 - 84867236911gydF4y2Ba KaracangydF4y2Ba c . O。gydF4y2Ba 鲁伊斯gydF4y2Ba f。gydF4y2Ba 象牙海岸gydF4y2Ba M。gydF4y2Ba 菲普斯gydF4y2Ba 年代。gydF4y2Ba 煤矿甲烷:捕捉和利用实践的回顾与矿业安全与温室气体减排效益gydF4y2Ba 国际煤炭地质杂志》上gydF4y2Ba 2011年gydF4y2Ba 86年gydF4y2Ba 2 - 3gydF4y2Ba 121年gydF4y2Ba 156年gydF4y2Ba 10.1016 / j.coal.2011.02.009gydF4y2Ba 2 - s2.0 - 79954570672gydF4y2Ba 张gydF4y2Ba C。gydF4y2Ba 你gydF4y2Ba 年代。gydF4y2Ba 程ydF4y2Ba M。gydF4y2Ba 张gydF4y2Ba lgydF4y2Ba 减压和甲烷生产压力释放的性能在长壁开采天然气开采技术gydF4y2Ba 地球物理与工程杂志》上gydF4y2Ba 2016年gydF4y2Ba 14gydF4y2Ba 1gydF4y2Ba 77年gydF4y2Ba 89年gydF4y2Ba 10.1088 / 1742 - 2140/14/1/77gydF4y2Ba 2 - s2.0 - 85010634310gydF4y2Ba LvgydF4y2Ba Y。gydF4y2Ba 试验研究岩石和气体流量的煤炭开采煤层周围gydF4y2Ba 国际矿业科技杂志》上gydF4y2Ba 2012年gydF4y2Ba 22gydF4y2Ba 4gydF4y2Ba 499年gydF4y2Ba 502年gydF4y2Ba 10.1016 / j.ijmst.2012.01.010gydF4y2Ba 2 - s2.0 - 84865227918gydF4y2Ba 道gydF4y2Ba y Q。gydF4y2Ba 徐gydF4y2Ba J。gydF4y2Ba 李gydF4y2Ba s . C。gydF4y2Ba 唐gydF4y2Ba x J。gydF4y2Ba 程gydF4y2Ba m·J。gydF4y2Ba 综合控制技术在工作面甲烷gydF4y2Ba 重庆大学学报gydF4y2Ba 2008年gydF4y2Ba 9gydF4y2Ba 1068年gydF4y2Ba 1073年gydF4y2Ba 张gydF4y2Ba lgydF4y2Ba 任gydF4y2Ba T。gydF4y2Ba 阿齐兹gydF4y2Ba N。gydF4y2Ba 温度和湿度的影响煤吸附特性的烟煤悉尼盆地,澳大利亚gydF4y2Ba 国际期刊的石油、天然气和煤炭技术gydF4y2Ba 2014年gydF4y2Ba 8gydF4y2Ba 1gydF4y2Ba 62年gydF4y2Ba 78年gydF4y2Ba 10.1504 / IJOGCT.2014.064429gydF4y2Ba 2 - s2.0 - 84906672146gydF4y2Ba 张gydF4y2Ba lgydF4y2Ba 任gydF4y2Ba T。gydF4y2Ba 阿齐兹gydF4y2Ba N。gydF4y2Ba 你gydF4y2Ba 年代。gydF4y2Ba 三轴渗透测试和微观结构的研究hard-to-drain煤炭从悉尼盆地,澳大利亚gydF4y2Ba 国际期刊的石油、天然气和煤炭技术gydF4y2Ba 2014年gydF4y2Ba 8gydF4y2Ba 4gydF4y2Ba 432年gydF4y2Ba 448年gydF4y2Ba 10.1504 / IJOGCT.2014.066311gydF4y2Ba 2 - s2.0 - 84918498647gydF4y2Ba 杨gydF4y2Ba W。gydF4y2Ba 林gydF4y2Ba b . Q。gydF4y2Ba 杨ydF4y2Ba Q。gydF4y2Ba 翟gydF4y2Ba C。gydF4y2Ba 长壁开采采场的应力再分配和多层煤层瓦斯治理gydF4y2Ba 国际岩石力学和采矿科学杂志》上gydF4y2Ba 2014年gydF4y2Ba 72年gydF4y2Ba 8gydF4y2Ba 15gydF4y2Ba 10.1016 / j.ijrmms.2014.08.009gydF4y2Ba 2 - s2.0 - 84910646978gydF4y2Ba PalchikgydF4y2Ba V。gydF4y2Ba 使用高斯分布估计的采空区瓦斯抽放效率gydF4y2Ba 数学地质gydF4y2Ba 2002年gydF4y2Ba 34gydF4y2Ba 6gydF4y2Ba 743年gydF4y2Ba 765年gydF4y2Ba 10.1023 /:1019857326863gydF4y2Ba 2 - s2.0 - 0038255616gydF4y2Ba 帕尔默gydF4y2Ba 我。gydF4y2Ba 煤炭:渗透率变化分析建模gydF4y2Ba 国际煤炭地质杂志》上gydF4y2Ba 2009年gydF4y2Ba 77年gydF4y2Ba 1 - 2gydF4y2Ba 119年gydF4y2Ba 126年gydF4y2Ba 10.1016 / j.coal.2008.09.006gydF4y2Ba 2 - s2.0 - 57649159089gydF4y2Ba 灰色的gydF4y2Ba 我。gydF4y2Ba 在煤层油藏工程:一部分时物理过程在煤层的气体存储和运动gydF4y2Ba SPE油藏工程gydF4y2Ba 1987年gydF4y2Ba 2gydF4y2Ba 1gydF4y2Ba 28gydF4y2Ba 34gydF4y2Ba 10.2118 / 12514 - pagydF4y2Ba SeidlegydF4y2Ba j . P。gydF4y2Ba JeansonnegydF4y2Ba d . J。gydF4y2Ba 埃里克森gydF4y2Ba d . J。gydF4y2Ba 火柴的应用几何应力依赖煤的渗透率gydF4y2Ba 学报SPE落基山区域会议gydF4y2Ba 1992年5月gydF4y2Ba 怀俄明gydF4y2Ba 433年gydF4y2Ba 444年gydF4y2Ba 10.2523 / 24361 - msgydF4y2Ba 任gydF4y2Ba t . X。gydF4y2Ba 爱德华兹gydF4y2Ba j·S。gydF4y2Ba 德gydF4y2Ba s E。gydF4y2Ba 采空区气体模型技术最大化甲烷捕捉从表面采空区井gydF4y2Ba 2002年gydF4y2Ba 矿井通风gydF4y2Ba 一来gydF4y2Ba d . N。gydF4y2Ba 朗兹gydF4y2Ba i S。gydF4y2Ba 金曼gydF4y2Ba s W。gydF4y2Ba 耶茨gydF4y2Ba C。gydF4y2Ba 失业gydF4y2Ba 年代。gydF4y2Ba 岩土工程因素对气体流动的影响在英国煤矿长壁板gydF4y2Ba 国际岩石力学和采矿科学杂志》上gydF4y2Ba 2006年gydF4y2Ba 43gydF4y2Ba 3gydF4y2Ba 369年gydF4y2Ba 387年gydF4y2Ba 10.1016 / j.ijrmms.2005.07.006gydF4y2Ba 2 - s2.0 - 30944447627gydF4y2Ba 帕尔默gydF4y2Ba 我。gydF4y2Ba MansoorigydF4y2Ba J。gydF4y2Ba 如何在煤层渗透性取决于应力和孔隙压力:一个新的模型吗gydF4y2Ba SPE储层评价和工程gydF4y2Ba 1998年gydF4y2Ba 1gydF4y2Ba 6gydF4y2Ba 539年gydF4y2Ba 544年gydF4y2Ba 10.2118 / 52607 - pagydF4y2Ba 帕尔默gydF4y2Ba 我。gydF4y2Ba MavorgydF4y2Ba M。gydF4y2Ba 甘特gydF4y2Ba B。gydF4y2Ba 在生产和注入煤层渗透率变化gydF4y2Ba 国际煤层气研讨会,阿拉巴马大学gydF4y2Ba 2007年gydF4y2Ba 阿拉巴马州的塔斯卡卢萨gydF4y2Ba 程gydF4y2Ba y . P。gydF4y2Ba 刘gydF4y2Ba h . Y。gydF4y2Ba 郭gydF4y2Ba p K。gydF4y2Ba 锅gydF4y2Ba r·K。gydF4y2Ba 王gydF4y2Ba lgydF4y2Ba 的理论模型和演化特征mining-enhanced更深的瓦斯煤层渗透率gydF4y2Ba 中国煤炭学会杂志》上gydF4y2Ba 2014年gydF4y2Ba 39gydF4y2Ba 8gydF4y2Ba 1650年gydF4y2Ba 1658年gydF4y2Ba 王gydF4y2Ba g·R。gydF4y2Ba 同步流的数值模拟通过双重孔隙度煤层甲烷和水的脱气过程中,[博士。论文)gydF4y2Ba 1985年gydF4y2Ba 宾夕法尼亚州立大学gydF4y2Ba 史gydF4y2Ba J。gydF4y2Ba 李gydF4y2Ba X。gydF4y2Ba 徐gydF4y2Ba B。gydF4y2Ba 杜gydF4y2Ba X。gydF4y2Ba 李gydF4y2Ba Y。gydF4y2Ba 温gydF4y2Ba 年代。gydF4y2Ba 张gydF4y2Ba D。gydF4y2Ba 回顾desorption-diffusion-flow煤层气的典范gydF4y2Ba Scientia中央研究院自然史,Mechanica & AstronomicagydF4y2Ba 2013年gydF4y2Ba 43gydF4y2Ba 12gydF4y2Ba 1548年gydF4y2Ba 1557年gydF4y2Ba 10.1360 / 132013 - 154gydF4y2Ba 太阳gydF4y2Ba J。gydF4y2Ba 故障特征的地板“three-zone”煤层沿倾斜方向gydF4y2Ba 《采矿与安全工程gydF4y2Ba 2014年gydF4y2Ba 1gydF4y2Ba 115年gydF4y2Ba 121年gydF4y2Ba 史gydF4y2Ba l Q。gydF4y2Ba 汉gydF4y2Ba J。gydF4y2Ba 理论和实践的分裂four-zone煤炭矿区的地板上gydF4y2Ba 中国矿业大学与技术杂志》上gydF4y2Ba 2005年gydF4y2Ba 34gydF4y2Ba 1gydF4y2Ba 16gydF4y2Ba 23gydF4y2Ba 钱gydF4y2Ba m·G。gydF4y2Ba 史gydF4y2Ba p W。gydF4y2Ba 我的压力和地层控制gydF4y2Ba 2003年gydF4y2Ba 徐州gydF4y2Ba 中国矿业大学和科技新闻gydF4y2Ba 史gydF4y2Ba l Q。gydF4y2Ba 徐gydF4y2Ba d . J。gydF4y2Ba 邱gydF4y2Ba M。gydF4y2Ba 京gydF4y2Ba X。gydF4y2Ba 太阳gydF4y2Ba H . H。gydF4y2Ba 改进的公式对损坏的地板在工作区域的深度gydF4y2Ba 中国煤炭学会杂志》上gydF4y2Ba 2014年gydF4y2Ba 38gydF4y2Ba 2gydF4y2Ba 299年gydF4y2Ba 303年gydF4y2Ba 关gydF4y2Ba y . B。gydF4y2Ba 李gydF4y2Ba h . M。gydF4y2Ba 陆gydF4y2Ba j . C。gydF4y2Ba 研究的。9煤层地板在西安Dewang煤矿的断裂规律gydF4y2Ba 中国煤炭学会杂志》上gydF4y2Ba 2003年gydF4y2Ba 28gydF4y2Ba 2gydF4y2Ba 121年gydF4y2Ba 125年gydF4y2Ba 张gydF4y2Ba Y。gydF4y2Ba 矿山压力的感应功能的损伤地板gydF4y2Ba 太原理工大学学报gydF4y2Ba 2002年gydF4y2Ba 33gydF4y2Ba 3gydF4y2Ba 252年gydF4y2Ba 256年gydF4y2Ba LisjakgydF4y2Ba 一个。gydF4y2Ba GrasselligydF4y2Ba G。gydF4y2Ba 回顾离散不连续岩体的破裂过程的建模技术gydF4y2Ba 《岩石力学与岩土工程gydF4y2Ba 2014年gydF4y2Ba 6gydF4y2Ba 4gydF4y2Ba 301年gydF4y2Ba 314年gydF4y2Ba 10.1016 / j.jrmge.2013.12.007gydF4y2Ba 2 - s2.0 - 84925233852gydF4y2Ba