CN105809561B - 一种含瓦斯煤弹性应变储能的计算方法 - Google Patents
一种含瓦斯煤弹性应变储能的计算方法 Download PDFInfo
- Publication number
- CN105809561B CN105809561B CN201610185876.7A CN201610185876A CN105809561B CN 105809561 B CN105809561 B CN 105809561B CN 201610185876 A CN201610185876 A CN 201610185876A CN 105809561 B CN105809561 B CN 105809561B
- Authority
- CN
- China
- Prior art keywords
- methane gas
- containing methane
- coal containing
- coal
- strain energy
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 title claims abstract description 146
- 239000003245 coal Substances 0.000 title claims abstract description 111
- 238000000205 computational method Methods 0.000 title claims abstract description 20
- 238000004146 energy storage Methods 0.000 title claims abstract description 13
- 238000000034 method Methods 0.000 claims abstract description 30
- 238000012360 testing method Methods 0.000 claims abstract description 14
- 238000004364 calculation method Methods 0.000 claims abstract description 11
- 238000002474 experimental method Methods 0.000 claims description 8
- 238000011160 research Methods 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 abstract description 2
- 239000011435 rock Substances 0.000 description 12
- 230000005489 elastic deformation Effects 0.000 description 5
- 230000006835 compression Effects 0.000 description 4
- 238000007906 compression Methods 0.000 description 4
- 238000005065 mining Methods 0.000 description 4
- 230000002265 prevention Effects 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 2
- 125000004122 cyclic group Chemical group 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000010339 dilation Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Marine Sciences & Fisheries (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Civil Engineering (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Structural Engineering (AREA)
- Mining & Mineral Resources (AREA)
- Architecture (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开一种含瓦斯煤弹性应变储能的计算方法,涉及煤岩力学特性研究技术领域,包括步骤1、获得含瓦斯煤应力‑体积应变曲线;步骤2、分析含瓦斯煤应力‑体积应变曲线,计算含瓦斯煤储存的弹性应变能;步骤3、进一步对含瓦斯煤弹性应变能计算新方法进行可靠性验证,以含瓦斯煤力学特性为依据,能够科学划分含瓦斯煤加载时弹性变形区段和塑性变形区段,扩容临界点前为弹性变形区段,扩容临界点后为塑性变形区段,有效解决了应力‑轴向应变曲线无法准确全面计算弹性应变能的难题,含瓦斯煤加载试验操作简单,无需反复加卸载试验,具有广泛的实用性。
Description
技术领域
本发明涉及煤岩力学特性研究技术领域,具体涉及一种含瓦斯煤弹性应变储能的计算方法。
背景技术
随着开采深度增加,深部煤炭开采动力灾害发生的频度和强度增大,严重威胁煤矿安全高效开采。深部煤炭开采动力灾害主要分为冲击地压和煤与瓦斯突出两类。随着开采向深部延伸,煤层瓦斯压力和瓦斯含量增大,深部煤层开采发生的动力灾害很难分清是冲击地压还是煤与瓦斯突出。目前,有关冲击地压和煤与瓦斯突出的防治,基于能量角度研究成果较丰富,也得到了国内外专家学者的认同。从能量的角度研究煤层动力灾害,认为煤与瓦斯动力灾害发生时,煤岩体内释放的能量主要来自外力所做的功和瓦斯内能,外力对煤岩体所做的功主要以弹性变形能的形式储存。这些能量在动力灾害发生过程中主要转化为煤体破碎耗散能和抛出时的动能。工作面煤体开采扰动而诱发的能量累积和释放过程是动力灾害发生的主要原因和内驱力。因此,对煤岩体内能量积聚及其转化过程的准确掌握和预测是有效防控煤与瓦斯动力灾害的关键。尤其是煤岩体储存的弹性应变能的计算准确与否直接关系煤与瓦斯动力灾害防控效果。
目前,煤体弹性应变能的计算主要是依据冲击地压能量判据的能量冲击指标法、弹性应变能指标法和岩爆能量指标法中对煤岩体弹性应变能的计算方法。能量冲击指标法等主要以煤体全应力-轴向应变曲线为基础,认为煤体抗压峰值强度前应力-轴向应变曲线投影围成的面积近似为煤体储存的弹性应变能。众所周知,煤岩体受外力作用至峰值强度时不但会发生弹性变形也发生塑性变形,而目前将煤体抗压峰值强度前应力-轴向应变曲线投影围成的面积定为煤体储存的弹性应变能,主要是忽略了峰值强度前部分塑性变形,这一计算方法对于不含瓦斯的岩石受力后储存的弹性应变能的计算基本能满足试验和实际要求,而对深部含瓦斯煤层,忽略部分塑性变形则会导致弹性应变能的计算结果存在很大偏差,使得煤与瓦斯动力灾害预判准确性大大降低。而若采用弹性应变能指数法计算弹性应变能,需开展加卸载实验,然而加卸载实验仅能开展峰前强度的实验(如峰值强度的70%、80%等),无法开展峰值强度的加卸载实验,由于循环加卸载试验条件限制,同样会导致弹性变形能计算结果存在偏差。因此,迫切需要对含瓦斯煤层受力作用后储存的弹性应变能的准确计算方法进行研究。
经研究发现,受高采动应力作用,含瓦斯煤发生扩容力学现象,考虑煤岩体发生扩容后应变变形主要以塑性变形为主,能够划分弹性变形区段和塑性变形区段,扩容临界点前为弹性变形区段,扩容临界点后为塑性变形区段,该方法以发生扩容临界点为界限,依据含瓦斯煤应力-体积应变曲线投影到横坐标围成的阴影部分面积计算弹性应变能,有效解决了应力-轴向应变曲线或循环加卸载试验无法准确全面计算弹性应变能的难题,为含瓦斯煤动力灾害有效防治提供了保障。
发明内容
针对现有技术的不足,本发明公开一种含瓦斯煤弹性应变储能的计算方法,能够解决现有忽略抗压峰值强度前塑性变形,导致弹性应变能的计算结果存在很大偏差的问题。
为实现以上目的,本发明通过以下技术方案予以实现:
一种含瓦斯煤弹性应变储能的计算方法,包括以下步骤:
S1:获得含瓦斯煤应力-体积应变曲线,取不同地质赋存条件煤样,测试含瓦斯煤应力-体积应变曲线,分析含瓦斯煤体积应变特征;进行三轴加载试验,其试验采用荷载加载方式;
S2:分析含瓦斯煤应力-体积应变曲线,计算含瓦斯煤储存的弹性应变能,针对含瓦斯煤应力-体积应变曲线特性,基于含瓦斯煤扩容力学现象,以含瓦斯煤扩容临界点为界,划分含瓦斯煤加载时弹性变形区段和塑性变形区段,扩容临界点前为弹性变形区段,扩容临界点后为塑性变形区段,以弹性变形区段应力-体积应变曲线投影到横坐标围成的阴影部分面积计算含瓦斯煤弹性应变能;
S3:进一步对含瓦斯煤弹性应变能计算新方法进行可靠性验证。
优选的,所述荷载加载方式为三轴应力加载方式。
优选的,所述步骤S1中,需在含瓦斯煤试样上的不同方向设置应变传感器,从而观测含瓦斯煤应力-体积应变曲线。
优选的,所述含瓦斯煤弹性应变能的计算方法可靠性依据不含瓦斯煤的弹性应变能指标法进行验证。
优选的,所述荷载加载过程的一次加压压差为0.005MPa/s。
本发明公开一种含瓦斯煤弹性应变储能的计算方法,以含瓦斯煤力学特性为依据,能够划分含瓦斯煤加载时弹性变形区段和塑性变形区段,扩容临界点前为弹性变形区段,扩容临界点后为塑性变形区段,有效解决了应力-轴向应变曲线无法准确全面计算弹性应变能的难题,含瓦斯煤加载试验操作简单,无需反复加卸载试验,具有广泛的实用性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明中含瓦斯煤应力-体积应变曲线图;
图2是本发明中常用弹性变形能计算示意图;
图3是本发明中大理岩应力-应变曲线;
图4是本发明的流程图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1至图4所示,本发明实施例所述一种含瓦斯煤弹性应变能的计算新方法,包括:
步骤1、获得含瓦斯煤应力-体积应变曲线
取不同地质赋存条件煤样,测试含瓦斯煤应力-体积应变曲线,分析含瓦斯煤体积应变特征:
进行三轴加载试验:利用MTS-816型岩石伺服气固耦合试验系统测试含瓦斯煤的应力-体积应变曲线,其中采用荷载加载方式,荷载加载过程的一次加压压差为0.005MPa/s;
测试含瓦斯煤三轴全应力-应变曲线时,在含瓦斯煤试样上的不同方向设置应变传感器,测定含瓦斯煤体积应变,从而观测含瓦斯煤体积应变由压缩转变为非线性膨胀拐点,即扩容临界点;
步骤2、分析含瓦斯煤应力-体积应变曲线,计算含瓦斯煤储存的弹性应变能;
含瓦斯煤应力-体积应变曲线,如图1所示。煤体加载时,其体积应变经历先压缩后膨胀的扩容过程,认为含瓦斯煤发生扩容后应变变形主要表现为塑性变形,故以含瓦斯煤扩容临界点(A点)为界,以加载开始扩容段(0-A段)体积应变曲线投影围成的面积计算弹性应变能。因此,含瓦斯煤储存的弹性应变能可用体积应变曲线与横坐标围成的阴影部分面积SOAB表示,故其储存弹性变形能We为:
式中:εv为体积应变值,εB为扩容临界点对应的体积应变值,σi为差应力。
由图1分析可知,阴影部分面积SOAB是曲边三角形,因此,SOAB面积采用定积分的原理进行计算,即将SOAB曲边形分成n个小曲边梯形,而每个小曲边梯形可以近似看成小矩形,将所有小矩形的面积之和作为曲边三角形面积,故实际计算可按式(2)进行。
We=∑σiεΔvi (2)
式中:εΔvi为σi时对应的单位体积应变值,σi为差应力。
步骤3、进一步对含瓦斯煤弹性应变能计算方法进行可靠性验证
目前,针对煤岩体加载过程中储存弹性应变能计算方法应用较多是:能量冲击性指标法、理想弹塑性模型的弹性能和弹性应变能指数法,如图2所示。
选取对锦屏大理岩具有岩爆倾向性岩石进行能量计算验证,如图3所示。
根据图3的轴向应变和体积应变曲线进行弹性能量计算,如表1所示。
表1大理岩应力–应变曲线弹性应变能计算结果
*注:实际储存弹性应变能以弹性应变能指数2计算
由表1分析可知,大理岩不同加载速度加载时,以弹性应变能指数法计算的实际储存弹性应变能分别为130.46焦耳、145.82焦耳,而以应力-体积应变曲线计算的弹性应变能分别为124.95焦耳、139.55焦耳。两种方法计算结果很接近,误差约为4.22%。由此表明,以扩容临界点为界的体积应变曲线计算弹性能的新方法的正确性、合理性。
综上所述本发明公开一种含瓦斯煤弹性应变储能的计算方法,以含瓦斯煤体力学特性为依据,能够划分含瓦斯煤加载时弹性变形区段和塑性变形区段,扩容临界点前为弹性变形区段,扩容临界点后为塑性变形区段,有效解决了应力-轴向应变曲线无法准确全面计算弹性应变能的难题含瓦斯煤加载试验操作简单,无需反复加卸载试验,具有广泛的实用性。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (5)
1.一种含瓦斯煤弹性应变储能的计算方法,其特征在于,包括以下步骤:
S1:获得含瓦斯煤应力-体积应变曲线,取不同地质赋存条件煤样,测试含瓦斯煤应力-体积应变曲线,分析含瓦斯煤体积应变特征;进行三轴加载试验,其试验采用荷载加载方式;
S2:分析含瓦斯煤应力-体积应变曲线,计算含瓦斯煤储存的弹性应变能,针对含瓦斯煤应力-体积应变曲线特性,基于含瓦斯煤扩容力学现象,以含瓦斯煤扩容临界点为界,划分含瓦斯煤加载时弹性变形区段和塑性变形区段,扩容临界点前为弹性变形区段,扩容临界点后为塑性变形区段,以弹性变形区段应力-体积应变曲线投影到横坐标围成的阴影部分面积计算含瓦斯煤弹性应变能;
S3:进一步对含瓦斯煤弹性应变能计算新方法进行可靠性验证。
2.如权利要求1所述含瓦斯煤弹性应变储能的计算方法,其特征在于:所述荷载加载方式为三轴应力加载方式。
3.如权利要求1所述含瓦斯煤弹性应变储能的计算方法,其特征在于:所述步骤S1中,需在含瓦斯煤试样上的不同方向设置应变传感器,从而观测含瓦斯煤应力-体积应变曲线。
4.如权利要求1所述含瓦斯煤弹性应变储能的计算方法,其特征在于:所述含瓦斯煤弹性应变能的计算方法可靠性依据不含瓦斯煤的弹性应变能指标法进行验证。
5.如权利要求1所述含瓦斯煤弹性应变储能的计算方法,其特征在于:所述荷载加载过程的一次加压压差为0.005MPa/s。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610185876.7A CN105809561B (zh) | 2016-03-25 | 2016-03-25 | 一种含瓦斯煤弹性应变储能的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610185876.7A CN105809561B (zh) | 2016-03-25 | 2016-03-25 | 一种含瓦斯煤弹性应变储能的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105809561A CN105809561A (zh) | 2016-07-27 |
CN105809561B true CN105809561B (zh) | 2018-07-03 |
Family
ID=56455037
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610185876.7A Active CN105809561B (zh) | 2016-03-25 | 2016-03-25 | 一种含瓦斯煤弹性应变储能的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105809561B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107991195A (zh) * | 2017-11-27 | 2018-05-04 | 中南大学 | 一种在三点弯曲断裂试验中测定岩石材料峰值荷载点处弹性变形能的方法 |
CN108051296A (zh) * | 2017-11-27 | 2018-05-18 | 中南大学 | 一种在复合型断裂试验中测定岩石材料峰值荷载点处弹性变形能的方法 |
CN107831069A (zh) * | 2017-11-27 | 2018-03-23 | 中南大学 | 一种在点荷载试验中测定岩石材料峰值荷载点处弹性变形能的方法 |
CN108170941B (zh) * | 2017-12-26 | 2021-07-27 | 东北大学 | 一种模具钢锻造过程等向性的预测方法 |
CN109253928B (zh) * | 2018-10-22 | 2021-05-28 | 山东科技大学 | 一种冲击地压启动能量阈值确定方法及冲击地压预测方法 |
CN112434418B (zh) * | 2020-11-20 | 2022-09-09 | 山东科技大学 | 一种开采引起冲击动能估算方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103454139A (zh) * | 2013-09-06 | 2013-12-18 | 安徽理工大学 | 含瓦斯煤岩体扩容关键影响因素及重要度确定方法 |
-
2016
- 2016-03-25 CN CN201610185876.7A patent/CN105809561B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103454139A (zh) * | 2013-09-06 | 2013-12-18 | 安徽理工大学 | 含瓦斯煤岩体扩容关键影响因素及重要度确定方法 |
Non-Patent Citations (2)
Title |
---|
主应力轴旋转对压实黄土动变形特性的影响;杨利国等;《工程地质学报》;20100131(第1期);参见第3.1节以及图3 * |
深部高瓦斯工作面控制煤体扩容主动防控方法;谢广祥等;《中国矿业大学学报》;20140531;第43卷(第3期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN105809561A (zh) | 2016-07-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105809561B (zh) | 一种含瓦斯煤弹性应变储能的计算方法 | |
CN105866373B (zh) | 一种煤与瓦斯动力灾害的新型判别方法 | |
Wang et al. | Caving mechanisms of loose top-coal in longwall top-coal caving mining method | |
Mitri, HS*, Tang, B.* & Simon | FE modelling of mining-induced energy release and storage rates | |
Brown | Progress and challenges in some areas of deep mining | |
Hu et al. | Model-based simulation of the synergistic effects of blast and fragmentation on a concrete wall using the MPM | |
Mase et al. | Non-linear site response analysis of soil sites in Northern Thailand during the M w 6.8 Tarlay earthquake | |
Hao et al. | A novel technology for enhancing coalbed methane extraction: Hydraulic cavitating assisted fracturing | |
CN107060761B (zh) | 一种近水平厚煤层综采采场基本顶初次破断距计算方法 | |
CN109100790A (zh) | 一种人工裂缝的模拟方法及装置 | |
Zhao et al. | Load transfer mechanism and reinforcement effect of segmentally yieldable anchorage in weakly consolidated soft rock | |
Yong et al. | Dynamic effect and control of key strata break of immediate roof in fully mechanized mining with large mining height | |
Kaiser et al. | Overcoming challenges of rock mass characterization for underground construction in deep mines | |
CN106815445A (zh) | 一种浅埋煤层区段煤柱载荷计算模型及其计算方法 | |
CN109918697A (zh) | 一种矿井冲击地压强度的预测方法及装置 | |
CN115410663A (zh) | 动态冲击/接触弹塑性大变形断裂分析显式相场物质点法 | |
CN104832163B (zh) | 一种煤矿地下开采中冲击地压危险性的监测方法 | |
CN106326636A (zh) | 一种基于可释放弹性应变能的岩爆碎块弹射速度预测方法 | |
CN106682757A (zh) | 一种井壁坍塌压力确定方法和装置 | |
CN111425180A (zh) | 一种ucg-ccs覆岩裂隙发育高度预测方法 | |
Hu et al. | Effect of Loading Rates on Mechanical Characteristics and Rock Burst Tendency of Coal‐Rock Combined Samples | |
Li et al. | Mechanism research on crack propagation in coal induced by CO2 phase-transition fracturing under different lateral compression coefficients | |
Wang et al. | Study on high and steep slope stability of surface mine based on RFPA-SRM | |
CN116167147B (zh) | 基于多层感知机算法的煤岩冲击倾向性直接指数评价方法 | |
Shen et al. | Research on the influence of parallel joints on the extension length and morphology of blasting cracks |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |