<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.3 20210610//EN" "JATS-journalpublishing1-3-mathml3.dtd">
<article article-type="brief-report" dtd-version="1.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:ali="http://www.niso.org/schemas/ali/1.0/" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Earth Sci.</journal-id>
<journal-title-group>
<journal-title>Frontiers in Earth Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Earth Sci.</abbrev-journal-title>
</journal-title-group>
<issn pub-type="epub">2296-6463</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1733210</article-id>
<article-id pub-id-type="doi">10.3389/feart.2026.1733210</article-id>
<article-version article-version-type="Version of Record" vocab="NISO-RP-8-2008"/>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Brief Research Report</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>Study of the hydraulic cavitation mechanism and aperture optimization of drilling in high-gas coal seams</article-title>
<alt-title alt-title-type="left-running-head">Li et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/feart.2026.1733210">10.3389/feart.2026.1733210</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Li</surname>
<given-names>Hongtao</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/3398302"/>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; original draft" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-original-draft/">Writing - original draft</role>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Xu</surname>
<given-names>Jinhai</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &#x26; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/">Writing - review and editing</role>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhang</surname>
<given-names>Zelin</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &#x26; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/">Writing - review and editing</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Data curation" vocab-term-identifier="https://credit.niso.org/contributor-roles/data-curation/">Data curation</role>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Li</surname>
<given-names>Wenxian</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Methodology" vocab-term-identifier="https://credit.niso.org/contributor-roles/methodology/">Methodology</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &#x26; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/">Writing - review and editing</role>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Fan</surname>
<given-names>Yongping</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &#x26; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/">Writing - review and editing</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Supervision" vocab-term-identifier="https://credit.niso.org/contributor-roles/supervision/">Supervision</role>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Jin</surname>
<given-names>Zhiyong</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &#x26; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/">Writing - review and editing</role>
</contrib>
</contrib-group>
<aff id="aff1">
<label>1</label>
<institution>School of Mines, State Key Laboratory for Fine Exploration and intelligent Development of Coal Resources, China University of Mining &#x26; Technology</institution>, <city>Xuzhou</city>, <state>Jiangsu</state>, <country country="CN">China</country>
</aff>
<aff id="aff2">
<label>2</label>
<institution>Zuoquan Wulinhou Mining Co., Ltd., Shanxi Lu&#x2019;an Group</institution>, <city>Jinzhong</city>, <country country="CN">China</country>
</aff>
<author-notes>
<corresp id="c001">
<label>&#x2a;</label>Correspondence: Jinhai Xu, <email xlink:href="mailto:441560931@qq.com">441560931@qq.com</email>
</corresp>
</author-notes>
<pub-date publication-format="electronic" date-type="pub" iso-8601-date="2026-02-24">
<day>24</day>
<month>02</month>
<year>2026</year>
</pub-date>
<pub-date publication-format="electronic" date-type="collection">
<year>2026</year>
</pub-date>
<volume>14</volume>
<elocation-id>1733210</elocation-id>
<history>
<date date-type="received">
<day>27</day>
<month>10</month>
<year>2025</year>
</date>
<date date-type="rev-recd">
<day>05</day>
<month>12</month>
<year>2025</year>
</date>
<date date-type="accepted">
<day>08</day>
<month>01</month>
<year>2026</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2026 Li, Xu, Zhang, Li, Fan and Jin.</copyright-statement>
<copyright-year>2026</copyright-year>
<copyright-holder>Li, Xu, Zhang, Li, Fan and Jin</copyright-holder>
<license>
<ali:license_ref start_date="2026-02-24">https://creativecommons.org/licenses/by/4.0/</ali:license_ref>
<license-p>This is an open-access article distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="https://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution License (CC BY)</ext-link>. The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</license-p>
</license>
</permissions>
<abstract>
<p>To enhance gas pre-extraction efficiency in low-permeability, high-gas coal seams and eliminate the danger of coal and gas outbursts at the Wulihou Coal Mine, this study analyzed the influence of the three stress zones surrounding cavitation holes (employed in alternating cavitation gas extraction technology) on gas extraction, based on the stress distribution characteristics around these holes. Two indices&#x2014;the pressure relief zone area ratio and the plastic zone area ratio&#x2014;were proposed to evaluate the pressure relief effect induced by the cavitation process. The effectiveness of cavitation holes for pressure relief and permeability enhancement was investigated using FLAC3D numerical simulation. Results indicate that as the cavitation hole size increases, the area of the associated stress relief zone expands, with its spatial extent best characterized by an exponential function. Similarly, the area of the plastic zone surrounding the cavitation hole also expands and is well-described by a cubic function. Both the pressure relief zone area ratio and the plastic zone area ratio exhibit an initial increase followed by a decrease with increasing cavitation hole radius, reaching their peak values within the radius range of [0.4, 0.5] m. Therefore, the optimal radius for cavitation holes is recommended to be between 0.4 m and 0.5 m. This research provides valuable reference for evaluating the pressure relief and permeability improvement effects, as well as for optimizing the aperture of hydraulic cavitation boreholes in high-gas coal seams.</p>
</abstract>
<kwd-group>
<kwd>cavity radius</kwd>
<kwd>gas extraction</kwd>
<kwd>hydraulic cavitation</kwd>
<kwd>soft and low-permeability coal seam</kwd>
<kwd>unloading-induced permeability increase</kwd>
</kwd-group>
<funding-group>
<funding-statement>The author(s) declared that financial support was not received for this work and/or its publication.</funding-statement>
</funding-group>
<counts>
<fig-count count="11"/>
<table-count count="1"/>
<equation-count count="16"/>
<ref-count count="34"/>
<page-count count="00"/>
</counts>
<custom-meta-group>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Geohazards and Georisks</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<label>1</label>
<title>Introduction</title>
<p>Coal is China&#x27;s dominant energy source, accounting for approximately 57.7% of China&#x27;s energy consumption in 2019 (<xref ref-type="bibr" rid="B10">Li, 2014</xref>; <xref ref-type="bibr" rid="B18">Lou et al., 2024</xref>; <xref ref-type="bibr" rid="B20">Lu et al., 2011</xref>). In the process of mine production, coal and gas outbursts and gas explosion accidents often cause heavy casualties and economic losses, which restrict the development of coal mines. Gas extraction is an effective method to eliminate gas accidents (<xref ref-type="bibr" rid="B33">Zhou et al., 2020</xref>; <xref ref-type="bibr" rid="B19">Lou et al., 2025</xref>). However, the permeability of coal bodies is typically low, ranging from (1&#x2013;100) &#xd7; 10<sup>&#x2013;3</sup> &#x3bc;m<sup>2</sup>, which poses challenges and inefficiencies in gas extraction (<xref ref-type="bibr" rid="B32">Zhang et al., 2023</xref>; <xref ref-type="bibr" rid="B2">Chen H et al., 2017</xref>; <xref ref-type="bibr" rid="B8">Huang et al., 2014</xref>). To enhance gas extraction efficiency effectively, corresponding measures to improve coal seam permeability are essential (<xref ref-type="bibr" rid="B5">Gao et al., 2015</xref>; <xref ref-type="bibr" rid="B21">Lu et al., 2017</xref>; <xref ref-type="bibr" rid="B4">Chen L et al., 2017</xref>). These include CO<sub>2</sub> or N<sub>2</sub> gas-phase fracturing (<xref ref-type="bibr" rid="B1">Busch and Gensterblum, 2011</xref>; <xref ref-type="bibr" rid="B28">Yang et al., 2022</xref>; <xref ref-type="bibr" rid="B34">Zhou et al., 2023</xref>; <xref ref-type="bibr" rid="B3">Chen et al., 2020</xref>), high-pressure water jet techniques, hydraulic fracturing, and hydraulic slotting, among others (<xref ref-type="bibr" rid="B9">Kong et al., 2016</xref>; <xref ref-type="bibr" rid="B7">He et al., 2021</xref>; <xref ref-type="bibr" rid="B29">Zhang et al., 2017</xref>; <xref ref-type="bibr" rid="B17">Liu et al., 2025</xref>; <xref ref-type="bibr" rid="B30">Zhang H et al., 2019</xref>). In engineering applications, the low-permeability characteristic of Chinese coal seams means gas drainage using these technologies still falls short of meeting actual safe production requirements (<xref ref-type="bibr" rid="B5">Gao et al., 2015</xref>; <xref ref-type="bibr" rid="B12">Li et al., 2024</xref>; <xref ref-type="bibr" rid="B13">Li et al., 2025</xref>).</p>
<p>
<xref ref-type="bibr" rid="B29">Zhang et al. (2017)</xref> proposed a novel in-seam borehole hydraulic flushing (ISBHF) technology for gas extraction. The consistency between field test results and simulation outcomes validated the reliability of the corresponding numerical model. After applying this new technology, the number of gas extraction boreholes at the heading face decreased by 50%, while the gas extraction concentration and flow rate increased by 10-fold and 6-fold, respectively. Consequently, the gas extraction duration was shortened by 65%, and the monthly excavation length of coal mine roadways doubled. For structural coal seams, <xref ref-type="bibr" rid="B31">Zhang R et al. (2019)</xref> conducted research on bedding hydraulic cavitation to strengthen gas drainage. The results indicated that this technology reduced the number of boreholes by 66.7% and increased the gas extraction rate and concentration by 1.33 times and 3 times, respectively, demonstrating significant advantages over traditional extraction methods.</p>
<p>Beyond hydraulic flushing and cavitation technologies, <xref ref-type="bibr" rid="B14">Lin et al. (2015)</xref> developed a cross-borehole hydraulic slotting technique. Field application data showed that after half a month of gas extraction, the average gas extraction concentration of slotted boreholes reached 26%, compared to only 7% for conventional boreholes&#x2014;meaning the concentration of slotted boreholes was approximately 3.7 times that of conventional ones, highlighting the remarkable effectiveness of slotting technology in improving extraction concentration. <xref ref-type="bibr" rid="B30">Zhang H et al. (2019)</xref> presented a hydraulic-flushing gas drainage technology. This technology efficiently flushed out more tectonic coal from thick coal seams, and the large-diameter cavities formed in the tectonic coal sub-layers significantly expanded the scope of stress relief and permeability enhancement, thereby remarkably improving gas extraction efficiency. To address the core demands of technological optimization and construction cost reduction, <xref ref-type="bibr" rid="B3">Chen et al. (2020)</xref> established a multi-physical coupling model. Through this model, they determined the optimal spacing of hydraulic flushing boreholes (HFB), providing precise parameter support for field tests and filling the gap in quantitative design within this field. However, existing hydraulic enhanced extraction technologies all have obvious limitations, a consensus widely recognized by numerous scholars (<xref ref-type="bibr" rid="B11">Li et al., 2023</xref>; <xref ref-type="bibr" rid="B16">Liu et al., 2023</xref>; <xref ref-type="bibr" rid="B23">Wang et al., 2021</xref>). For instance, hydraulic slotting technology exhibits poor sustainability in permeability enhancement effects. Under high <italic>in situ</italic> stress conditions, the fractures generated by slotting are prone to closure, leading to a rapid decline in extraction efficiency. Hydraulic fracturing technology, on the other hand, tends to cause local stress concentration in coal masses, which in turn restricts the permeability improvement effect. More critically, both technologies demonstrate poor performance in gas extraction from soft tectonic coal seams. Soft tectonic coal inherently has low strength and disorganized fracture development. Fractures from hydraulic slotting are liable to collapse and close, while hydraulic fracturing may result in overall coal fragmentation that clogs boreholes. This has become a key bottleneck limiting gas control in soft coal seams.</p>
<p>In contrast, borehole hydraulic cavitation has emerged as an innovative pressure relief and permeability enhancement technology, attracting growing attention and extensive research in recent years. The core process of this technology involves three steps: first, drilling conventional boreholes into the coal seam; second, fracturing the coal surrounding the boreholes using high-pressure water jets; and finally, flushing the fragmented coal out of the boreholes to form large-diameter cavities. Compared with competing technologies such as hydraulic slotting and fracturing, its unique advantage lies in achieving long-term stress relief by forming physical cavities rather than relying on fracture structures that are prone to closure. This provides a novel solution to the challenges of gas extraction in soft tectonic coal seams.</p>
<p>Gas cavitation plays a crucial role in the safe and efficient extraction of coalbed methane. Although previous studies have predominantly focused on seepage characteristics, gas cavitation exerts a significant impact on coal mass by altering its stress distribution: through stress relief and concentration effects, this stress change promotes gas desorption from coal matrix pores, increasing coal permeability, and thereby achieves effective gas extraction. Meanwhile, gas cavitation damages the coal surrounding the cavitation hole, leading to the formation of a plastic zone in the affected area. Notably, gas cavitation disrupts the original stress balance of the coal mass, resulting in the formation of a two-zone annular structure within the coal, which is composed of a stress concentration zone and an original stress zone. As the core disturbed region, the stress concentration zone undergoes the redistribution of radial and tangential stresses, accompanied by corresponding variations in coal permeability, and can be further subdivided into two sub-zones. In the fluidized zone, tangential compressive stress drives the opening of coal fractures, leading to a sharp increase in permeability; in the compaction zone, radial compressive stress dominates the coal compaction process, resulting in a temporary reduction in permeability. In contrast, the original stress zone remains in an undisturbed state, with its physical&#x2013;mechanical properties and seepage characteristics maintaining their initial conditions. Numerical simulations have been conducted to analyze the stress distribution around the cavitation hole and the formation characteristics of the plastic zone. Relevant indicators for evaluating the effectiveness of gas cavitation have been proposed based on field practice, which provide a solid theoretical foundation for the practical application of gas cavitation technology in coalbed methane extraction.</p>
<p>Existing research has not yet clearly studied the mechanism of borehole hydraulic cavitation. Therefore, based on the redistribution of coal stress caused by the excavation of the hole, this article attempts to explain the mechanism of pressure relief and permeability enhancement of the hole from the perspective of stress. In this article, based on the simulation of engineering practice, limited by the field parameters, only simplified models are used instead of complex models. The innovation is that it is not limited to the conventional gas seepage analysis, but instead offers a new way to analyze the stress change and coal damage in the process of gas cavitation from the mechanical failure mechanism.</p>
<p>Initially, the mechanism of the distressed zone and permeability increase is examined based on the stress distribution law surrounding the cavitation hole. Subsequently, the finite-difference FLAC3D numerical analysis software is utilized in the context of the gas extraction project at the Wulihou Coal Mine in Shanxi, China, to investigate the optimal radius of the cavitation hole. This investigation is conducted by comparing and analyzing variations in the area of the pressure relief zone, the pressure relief area ratio, the plastic zone area, and the plastic zone area ratio for different cavitation hole radii. On this basis, combined with practical engineering, a set of quantitative indicators has been carefully established to accurately evaluate the effect of gas cavitation, injecting new vitality into engineering practice.</p>
</sec>
<sec id="s2">
<label>2</label>
<title>Project overview and engineering geological condition</title>
<p>The Wulihou Coal Mine is situated in Zuoquan County, Shanxi Province. The coal seam has a thickness of approximately 6 m, with an absolute gas outflow of approximately 18 cubic meters per minute and a relative gas outflow of 11 cubic meters per ton, indicating it is a high-gas mine. The maximum principal stress of the lower coal group in Wulihou Coal Mine is horizontal stress, with an average value of 18.5 MPa and a direction angle of N31.65 E &#x223c; N36.53 E; the minimum principal stress is still horizontal stress, with an average value of 10.5 MPa and a direction angle of N121.32 E &#x223c; N125.32 E; the intermediate principal stress is vertical stress, with an average value of 11.5 MPa and a direction angle of is N214.95&#xb0; E &#x223c; N219.83&#xb0; E. The lower coal seam has a large dip angle, and the strength of the roof and floor mudstone and coal body is small; they are unstable rock strata. Furthermore, the coal seams exhibit limited permeability and diminished compressive strength, complicating gas extraction and heightening the associated risks of coal and gas outbursts.</p>
<p>To guarantee the secure excavation of the belt highway and the track roadway, alternating cavitation technology is employed for gas extraction, as seen in <xref ref-type="fig" rid="F1">Figure 1</xref>. The technical principle is as follows: first, ordinary boreholes are drilled into the coal seam in the belt bottom drainage roadway, subsequently employing high-pressure water to generate current as the energy source to construct a series of cylindrical cave-building chambers within the coal seam. This causes the coal body to unload and increase permeability, thereby facilitating the extraction of coal seam gas. To increase production efficiency and decrease costs, the interval is alternately utilized to position standard drilling holes and cavitation holes. Specifically, a row of acupoint holes is situated between two rows of standard drilling. However, existing research is insufficient to understand the pressure relief and permeability enhancement effects around the cavitation hole in the coal body. The radius (R) of the cavitation hole may influence the coal unloading and permeability increase, necessitating more investigation.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Gas extraction schematic diagram of interval alternating cavitation hole-making technology.</p>
</caption>
<graphic xlink:href="feart-14-1733210-g001.tif">
<alt-text content-type="machine-generated">Diagram showing a coal seam area with labeled elements including cavitation hole, common borehole, belt roadway, track roadway, and belt bottom drainage roadway, with red lines representing borehole directions connecting to various points in the coal seam.</alt-text>
</graphic>
</fig>
</sec>
<sec id="s3">
<label>3</label>
<title>Mechanism of pressure relief and permeability enhancement by cavitation holes</title>
<p>The primary coal structure is composed of pores, cracks, coal matrix blocks, and gas, as shown in <xref ref-type="fig" rid="F2">Figure 2</xref>. It is assumed that the coal matrix, pores, cracks, and gas are subjected to the original rock stress (&#x3c3;<sub>1</sub>, &#x3c3;<sub>2</sub>, &#x3c3;<sub>3</sub>) before the cave is formed; that is, in the original rock stress state. Then the coal body is caved, the original rock stress balance of the coal body is broken, and the coal body stress is redistributed around the cavitation. After the stress equilibrates, three zones will be formed around the cavitation hole: a stress-decreasing zone (referred to as the stress-relaxation area), the stress concentration zone, and the initial stress area. Within the stress-decreasing zone (the pressure relief shown in <xref ref-type="fig" rid="F2">Figure 2</xref>), the coal body expands the acupoints around the cavitation hole, driving the pores and cracks to deform. The pores and cracks also expand under the effect of the gas pressure (such as the gas pressure in <xref ref-type="fig" rid="F2">Figure 2</xref>) inside the hole, propagating the pores and cracks, which leads to an increase in the permeability of the coal body and a decrease in the pressure of the reservoir. As a result of the changes, the gas in the coal body desorbs, diffuses, and seeps into the cavitation. When the stress exceeds the strength limit of the coal body in the stress concentration area, the coal matrix block around the cavitation hole is damaged, and new cracks sprout and connect with the original pores and cracks, which increases the permeability of the coal body. The gas in the coal body diffuses and seeps into the cavitation hole. Within the primary rock stress zone area, the coal body is not affected by the disturbance of the cavitation hole, and it has no effect on the gas in the coal body.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Gas extraction schematic diagram.</p>
</caption>
<graphic xlink:href="feart-14-1733210-g002.tif">
<alt-text content-type="machine-generated">Diagram illustrating coal matrix with labeled pores, fractures, cavitation boundary, and gas pressure. Multiple arrows indicate principal stresses (sigma one, sigma two, sigma three), stress relief, and directions of force application within the coal structure.</alt-text>
</graphic>
</fig>
<p>The above analysis shows that the essence of pressure relief is to form a stress reduction zone and a stress concentration area around the cavitation hole, which causes internal changes in the coal body and increases the permeability of the coal body. For this reason, the issue will be analyzed from the perspective of stress distribution. It is assumed that the coal rock mass, overburden rock, and floor are isotropic linear elastic bodies (<xref ref-type="bibr" rid="B22">Ma et al., 2023</xref>; <xref ref-type="bibr" rid="B26">Xue et al., 2013</xref>). There is an excavation in the infinite elastic body to create acupuncture holes, simplifying the space problem into a plane strain problem. The force diagram of the cavitation hole is shown in <xref ref-type="fig" rid="F3">Figure 3</xref>. The vertical and horizontal original rock stresses are p and q, respectively; the radius of the cavitation hole is a; the angle of horizontal axis from the horizontal t to the right to the right time to the unit body is <inline-formula id="inf1">
<mml:math id="m1">
<mml:mrow>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>; the radial stress is <inline-formula id="inf2">
<mml:math id="m2">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, the ring stress is <inline-formula id="inf3">
<mml:math id="m3">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, the shear stress is <inline-formula id="inf4">
<mml:math id="m4">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c4;</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, and the side pressure coefficient k &#x3d; p/q. The stress on the cavitation hole is the same as that on the horizontal to the right horizontal axis. According to the theory of elasticity, the stresses acting on any tiny cell are (<xref ref-type="disp-formula" rid="e1">Equation 1</xref>).<disp-formula id="e1">
<mml:math id="m5">
<mml:mrow>
<mml:mfenced open="" close="}" separators="&#x7c;">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mi>a</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>4</mml:mn>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mi>a</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>3</mml:mn>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mi>a</mml:mi>
<mml:mn>4</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>4</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mi>a</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>3</mml:mn>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mi>a</mml:mi>
<mml:mn>4</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>4</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c4;</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mi>a</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>3</mml:mn>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mi>a</mml:mi>
<mml:mn>4</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>4</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>sin</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>
</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Force diagram of the cavitation hole.</p>
</caption>
<graphic xlink:href="feart-14-1733210-g003.tif">
<alt-text content-type="machine-generated">Diagram showing a square domain with side length H, subjected to vertical load p and horizontal load q. A circular hole of radius a is centered with polar coordinates r and &#x3B8; labeled. Stresses &#x3C3;r, &#x3C3;&#x3B8;, and shear stress &#x3C4;r&#x3B8; are indicated at a differential element near the hole.</alt-text>
</graphic>
</fig>
<p>Based on the vertical stress of the primary stress state on most global coal bodies, the vertical stress of the original rock stress is greater than the horizontal stress, that is, <inline-formula id="inf5">
<mml:math id="m6">
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>. The distribution of cutting stress <inline-formula id="inf6">
<mml:math id="m7">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> along the periphery of the hole is shown in <xref ref-type="fig" rid="F3">Figure 3</xref>. The shape of the circular stress distribution curve is related to the value of K. This curve shape has a value of <inline-formula id="inf7">
<mml:math id="m8">
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>. If <inline-formula id="inf8">
<mml:math id="m9">
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>, the top (B-B line direction) of the cavitation hole will occur. According to the principle of pressure as &#x201c;&#x2b;&#x201d; and &#x201c;&#x2212;&#x201d; according to rock mechanics, the ring stress distribution curve at the top will be moved. The shape of the acupoint is roughly similar. Therefore, in order not to lose the general nature, the shape of the hoop stress <inline-formula id="inf9">
<mml:math id="m10">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> distribution curve in <xref ref-type="fig" rid="F4">Figure 4</xref> is analyzed.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Distribution curve of hoop stress <italic>&#x3c3;&#x3b8;</italic> around the cavitation hole.</p>
</caption>
<graphic xlink:href="feart-14-1733210-g004.tif">
<alt-text content-type="machine-generated">Diagram shows a central circle labeled as excavation boundary with a surrounding irregular curve labeled as hoop stress distribution curve, intersected by dashed AA and solid BB axes, and labeled arrows for reference.</alt-text>
</graphic>
</fig>
<p>
<xref ref-type="fig" rid="F4">Figure 4</xref> shows that the hoop stress <inline-formula id="inf10">
<mml:math id="m11">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> around the periphery of the cavitation hole has a stress reduction in the direction of the B-B line; that is, it is the stress of the force release or decompression effect. Therefore, it is proposed to analyze the decompression effect of the cavitation hole in the B-B line by the hoop stress in the B-B line direction <inline-formula id="inf11">
<mml:math id="m12">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> In addition, in order to quantify the scope of influence of the decompression effect, the decompression area ratio is used as a judgment index, which is calculated by the following formula:<disp-formula id="e2">
<mml:math id="m13">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>
</p>
<p>where &#x3b7;<sub>s</sub> denotes the pressure relief area ratio; S<sub>r</sub> represents the area where the hoop stress &#x3c3;<sub>&#x3b8;</sub> is diminished following cavitation; S signifies the area of the cavitation hole.</p>
<p>The hoop stress around the cavitation hole is <inline-formula id="inf12">
<mml:math id="m14">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>. There is an increase in stress in the direction of the A-A line; that is, it is a stress concentration, and the stress should be the force concentration that will lead to shear damage of the coal matrix, sprouting new cracks, and increasing the permeability of the coal body. In order to quantify the scope of influence of the increased permeability of the coal body due to the stress concentration, the area ratio of the plastic zone is used as a judgment index, and its calculation formula is as follows:<disp-formula id="e3">
<mml:math id="m15">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
</p>
<p>where &#x3b7;<sub>p</sub> denotes the plastic zone area ratio; S<sub>FP</sub> represents the plastic zone area after cavitation.</p>
</sec>
<sec id="s4">
<label>4</label>
<title>Numerical simulation</title>
<sec id="s4-1">
<label>4.1</label>
<title>Model building</title>
<p>To analyze the effect of unloading pressure and increasing penetration of different cavitation hole diameters, FLAC3D was utilized to create a three-dimensional numerical grid model, as illustrated in <xref ref-type="fig" rid="F5">Figure 5</xref>. The dimensions of the model are 12 m &#xd7; 10 m &#xd7; 6 m, with a stress boundary condition applied at the top and a displacement boundary condition at the bottom, front, back, left, and right sides. The coal body utilizes the Moore&#x2013;Cullen model, and its mechanical parameters are presented in <xref ref-type="table" rid="T1">Table 1</xref>.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Drill hole arrangement for gas extraction. <bold>(a)</bold> Coal body model. <bold>(b)</bold> Common borehole and cavitation holes in the perspective view of the coal body model.</p>
</caption>
<graphic xlink:href="feart-14-1733210-g005.tif">
<alt-text content-type="machine-generated">Panel a shows a three-dimensional turquoise cube composed of a grid mesh. Panel b displays a transparent rectangular prism labeled with dimensions twelve meters by ten meters by six meters, highlighting colored tubes identified as common borehole and cavitation hole.</alt-text>
</graphic>
</fig>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Mechanical parameters of the coal body (<xref ref-type="bibr" rid="B24">Wang et al., 2024</xref>; <xref ref-type="bibr" rid="B6">Guo et al., 2023</xref>).</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">
<inline-formula id="inf13">
<mml:math id="m16">
<mml:mrow>
<mml:mi>&#x3c1;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>/(kg <inline-formula id="inf14">
<mml:math id="m17">
<mml:mrow>
<mml:mo>&#xb7;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> m<sup>&#x2212;3</sup>)</th>
<th align="center">
<italic>K</italic>/GPa</th>
<th align="center">
<italic>G</italic>/GPa</th>
<th align="center">
<inline-formula id="inf15">
<mml:math id="m18">
<mml:mrow>
<mml:mi>&#x3c6;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>/(<inline-formula id="inf16">
<mml:math id="m19">
<mml:mrow>
<mml:mo>&#xb0;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula>)</th>
<th align="center">
<italic>C</italic>
<sub>
<italic>m</italic>
</sub>/MPa</th>
<th align="center">
<inline-formula id="inf17">
<mml:math id="m20">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>/MPa</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">1410</td>
<td align="left">2.0</td>
<td align="left">0.88</td>
<td align="left">20</td>
<td align="left">0.5</td>
<td align="left">0.2</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>
<inline-formula id="inf18">
<mml:math id="m21">
<mml:mrow>
<mml:mi mathvariant="bold-italic">&#x3c1;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>, density; <italic>K</italic>, bulk modulus; <italic>G</italic>, shear modulus; <inline-formula id="inf19">
<mml:math id="m22">
<mml:mrow>
<mml:mi mathvariant="bold-italic">&#x3c6;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>, friction angle; <italic>C</italic>
<sub>
<italic>m</italic>
</sub>, cohesion; <inline-formula id="inf20">
<mml:math id="m23">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">&#x3c3;</mml:mi>
<mml:mi mathvariant="bold-italic">t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, tensile strength.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>The numerical simulation methodology is delineated as follows: establish the cavitation hole radii at 0.3 m, 0.4 m, 0.5 m, 0.6 m, and 0.7 m, respectively, to investigate the effects of varying cavitation hole diameters on the plastic zone, maximum principal stress, and shear strain increment of the surrounding coal body. The distance between the cavitation hole and the standard drill holes is 3 m, and the radius of the standard drill holes is 0.1 m.</p>
<p>When the cavitation hole has a radius of 0.7 m, the top and bottom boundaries are 3 m from the hole&#x2019;s center, approximately 4.29 times the radius. For the numerical model, the boundary extent shall be specified as six times the radius, pursuant to <xref ref-type="disp-formula" rid="e4">Equations 5</xref>&#x2013;<xref ref-type="disp-formula" rid="e8">8</xref> (<xref ref-type="bibr" rid="B15">Liu and Yu, 2012</xref>; <xref ref-type="bibr" rid="B25">Xu, 2016</xref>). Beyond this range, rock mass variation is approximately 3%, which is negligible. In this case, with a 4.29&#xd7; radius boundary, the rock mass variation is approximately 5%, still within an acceptable range. This choice minimizes the influence of external rock masses while avoiding a significant increase in grid count and computational load, thus balancing resource use and meeting engineering analysis requirements.<disp-formula id="e4">
<mml:math id="m24">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x3b3;</mml:mi>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msubsup>
<mml:mi>r</mml:mi>
<mml:mn>0</mml:mn>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
<disp-formula id="e5">
<mml:math id="m25">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x3b3;</mml:mi>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msubsup>
<mml:mi>r</mml:mi>
<mml:mn>0</mml:mn>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>
<disp-formula id="e6">
<mml:math id="m26">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c4;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>
<disp-formula id="e7">
<mml:math id="m27">
<mml:mrow>
<mml:mi>u</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x3b3;</mml:mi>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:msubsup>
<mml:mi>r</mml:mi>
<mml:mn>0</mml:mn>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>G</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
<disp-formula id="e8">
<mml:math id="m28">
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>
</p>
<p>In this study, the radius of conventional boreholes is 0.1 m, so the mesh size is set to 0.1 m during modeling, resulting in a mesh count of 278,414. Expanding the boundary would significantly increase the mesh count. For these reasons, the current model boundary is selected.</p>
</sec>
<sec id="s4-2">
<label>4.2</label>
<title>Analysis of the pressure relief effects of the cavitation hole</title>
<sec id="s4-2-1">
<label>4.2.1</label>
<title>Stress relief zone analysis</title>
<p>The distribution of maximum principal stresses for different radii of cavitation holes is shown in <xref ref-type="fig" rid="F6">Figure 6</xref>, where the negative numbers indicate rock mass compression. <xref ref-type="fig" rid="F6">Figure 6a</xref> shows that when the radius of the cavitation hole is 0.3 m, the minimum value of the maximum principal stress at the top of the cavitation hole is 1.28 MPa. Compared with the stress in the original rock of 2.39 MPa, it is reduced by 1.11 MPa, and the maximum stress-release rate is 46%. <xref ref-type="fig" rid="F6">Figure 6b</xref> shows that when the radius of the cavitation hole is 0.4 m, the stress-release area in the vertical direction of the cavitation hole is in the shape of a flat long ellipse, and its area is larger than that of the cavitation hole with a radius of 0.3 m, while the stress concentration area in the horizontal direction is in the shape of a circle, and its area is larger than that of the cavitation hole with a radius of 0.3 m. The minimum value of the maximum principal stress at the top of the cavitation hole is 1.05 MPa, which is 1.34 MPa less than the original rock stress, and the maximum stress-release rate is 56%. <xref ref-type="fig" rid="F6">Figure 6c</xref> shows that when the radius of the cavitation hole is 0.5 m, the stress-release area in the vertical direction of the cavitation hole is in the shape of a long flat ellipse, and its area is larger than that of the cavitation hole with a radius of 0.4 m, while the stress concentration area in the horizontal direction is in the shape of a goldfish, and its area is larger than that of the cavitation hole with a radius of 0.4 m. The minimum value of the maximum principal stress at the top of the cavitation hole is 0.876 MPa, which is 1.11 MPa less than the original rock stress, and the maximum stress-release rate is 63%. <xref ref-type="fig" rid="F6">Figure 6d</xref> shows that when the radius of the cavitation hole is 0.6 m, the stress relief area in the vertical direction of the cavitation hole is in the shape of a long flat ellipse, and its area is larger than that of the cavitation hole. The radius of the cavitation hole is 0.5 m larger, and the stress concentration area in the horizontal direction is in the shape of a goldfish, and its area is half that of the cavitation hole. The diameter of 0.5 m is larger. The minimum value of the maximum principal stress at the top of the cavitation hole is 0.849 MPa, which is 1.54 MPa less than the original rock stress, and the maximum stress-release rate is 64%. <xref ref-type="fig" rid="F6">Figure 6e</xref> shows that when the radius of the cavitation hole is 0.7 m, the stress relief area in the vertical direction of the cavitation hole is in the shape of a long flat ellipse, and its area is larger than that of the cavitation hole with a radius of 0.7 m. When the radius is 0.6 m, the stress relief area decreases, whereas the horizontal stress-concentration zone exhibits a goldfish-like shape, with an area larger than that corresponding to a cavitation hole radius of 0.6 mm. The minimum value of the maximum principal stress at the top of the cavitation hole is 0.77 MPa, which is 1.62 MPa less than the original rock stress, and the maximum stress-release rate is 68%.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Distribution of maximum principal stresses for different radii of cavitation holes: <bold>(a)</bold> 0.3 m, <bold>(b)</bold> 0.4 m, <bold>(c)</bold> 0.5 m, <bold>(d)</bold> 0.6 m, and <bold>(e)</bold> 0.7 m.</p>
</caption>
<graphic xlink:href="feart-14-1733210-g006.tif">
<alt-text content-type="machine-generated">Five color contour plots (labeled a to e) display the variation of principal stress &#x3C3;&#x2081; in megapascals across a coal body of fixed dimensions, each with a progressively larger white circular void near the center, indicating different borehole sizes. Each plot uses a color scale from red (higher stress values) to blue (lower stress values) to visualize stress distribution, with axes labeled &#x22;coal body height in meters&#x22; and &#x22;length of the coal body in meters.&#x22;</alt-text>
</graphic>
</fig>
<p>From the shape of the principal stress distribution, when the radius of the cavitation hole changes from 0.3 m to 0.7 m, the stress-release area in the vertical direction of the cavitation hole is flattened with an elongated ellipse shape. However, the stress concentration area in the horizontal direction varies. When the radius of the cavitation hole is 0.3 m or 0.4 m, the stress concentration area assumes a circular shape. When the diameter of the cavitation hole measures 0.5 m, 0.6 m, or 0.7 m, the stress concentration region in the horizontal axis exhibits a goldfish-like morphology. As the radius expands, the zones of stress release and horizontal stress concentration in the vertical direction progressively enlarge.</p>
<p>
<xref ref-type="fig" rid="F6">Figures 6a&#x2013;e</xref> show that, compared to the stress in the original rock, the maximum principal stress in the vertical direction adjacent to the cavitation hole diminishes, while the maximum principal stress in the horizontal direction escalates. As the hole diameter increases, the maximum principal stress in the vertical direction near the cavitation hole continues to decrease significantly, resulting in an expanding stress-release area. The reason is that as the radius of the cavitation cavity increases, the disruption to the coal body intensifies, allowing for greater expansion of the coal body and thus enlarging the stress-release region, so facilitating gas seepage.</p>
<p>The area of the pressure relief zone corresponding to the radius of the cavitation hole is extracted, and a plot of the area of the pressure relief zone against the radius of the cavitation hole is illustrated in <xref ref-type="fig" rid="F7">Figure 7</xref>, which shows that the area of the pressure relief zone progressively enlarges with the radius of the cavitation hole. When the radius of the cavitation hole is &#x2265; 0.4 m, the increase is approximately linear, suggesting that the radius of the cavitation hole significantly influences the area of the pressure relief zone.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Graph of the variation of the area of the pressure relief zone with the radius of the cavitation hole.</p>
</caption>
<graphic xlink:href="feart-14-1733210-g007.tif">
<alt-text content-type="machine-generated">Line chart with black square markers shows area of pressure relief zone in meters on the vertical axis and cavitation hole radius in meters on the horizontal axis. Data points display a positive, nonlinear correlation.</alt-text>
</graphic>
</fig>
<p>To better describe the relationship between the area of the pressure relief zone and the radius of the cavitation hole, two reasonable fitting model assumptions are proposed for the variation curve of the pressure relief zone area with the cavitation hole radius. Two more reasonable fitting model assumptions can be made for the curves of the radius change plots (<xref ref-type="disp-formula" rid="e9">Equations 9</xref>&#x2013;<xref ref-type="disp-formula" rid="e12">12</xref>):<list list-type="order">
<list-item>
<label>1.</label>
<p>Parabolic model. The curve of the model is</p>
</list-item>
</list>
<disp-formula id="e9">
<mml:math id="m29">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:msup>
<mml:mi>R</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mi>R</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>
</p>
<p>where <italic>a</italic>
<sub>3</sub>, <italic>b</italic>
<sub>3</sub>, and <italic>c</italic>
<sub>3</sub> are parameters to be determined in relation to the area of the pressure relief zone.<list list-type="order">
<list-item>
<label>2.</label>
<p>Exponential function model. The curve of the model is</p>
</list-item>
</list>
<disp-formula id="e10">
<mml:math id="m30">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:msup>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:mi>R</mml:mi>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>where <italic>a</italic>
<sub>4</sub>, <italic>b</italic>
<sub>4</sub>, and <italic>c</italic>
<sub>4</sub> are parameters to be determined with respect to the area of the pressure relief zone.</p>
<p>Based on the extracted values for the radius of the cavitation hole and the area of the pressure relief zone, the fitted curves are obtained, respectively,<disp-formula id="e11">
<mml:math id="m31">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>11.9</mml:mn>
<mml:msup>
<mml:mi>R</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>23.2</mml:mn>
<mml:mi>R</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>5.2</mml:mn>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>
<disp-formula id="e12">
<mml:math id="m32">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1.4</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>15</mml:mn>
<mml:mo>&#xd7;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>.</mml:mo>
<mml:msup>
<mml:mn>1</mml:mn>
<mml:mi>R</mml:mi>
</mml:msup>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>
</p>
<p>The coefficients of determination of the fit are, respectively, <inline-formula id="inf21">
<mml:math id="m33">
<mml:mrow>
<mml:msubsup>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0.9920</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>; <inline-formula id="inf22">
<mml:math id="m34">
<mml:mrow>
<mml:msubsup>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0.9942</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<p>The above analysis shows that the relationship between the area of the pressure relief zone of the cavitation hole and the radius of the cavitation hole is better fitted by the exponential function model, and the exponential function model can better describe the change of the area of the pressure relief zone of the cavitation hole with the change of the radius of the cavitation hole.</p>
</sec>
<sec id="s4-2-2">
<label>4.2.2</label>
<title>Plastic zone analysis</title>
<p>Tensile failure occurs when the tensile stress <inline-formula id="inf23">
<mml:math id="m35">
<mml:mrow>
<mml:mi>&#x3c3;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> exceeds the tensile strength <inline-formula id="inf24">
<mml:math id="m36">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> of coal, while shear failure occurs when the shear stress <inline-formula id="inf25">
<mml:math id="m37">
<mml:mrow>
<mml:mi>&#x3c4;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> exceeds the shear strength <inline-formula id="inf26">
<mml:math id="m38">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c4;</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> of coal. This kind of strength failure can be intuitively characterized by the distribution characteristics of the plastic zone. <xref ref-type="fig" rid="F8">Figure 8</xref> shows the distribution state of the plastic zone under different cavitation hole radii, where &#x201c;none&#x201d; indicates no yielding, &#x201c;shear&#x201d; indicates shear yielding, &#x201c;tension&#x201d; indicates tensile yielding, &#x201c;n&#x201d; indicates entering yielding at the current step, and &#x201c;p&#x201d; indicates entering yielding at the previous step. The core value of this figure lies in clearly presenting the entire process of coal fracture initiation and propagation dominated by tangential stress, and the difference in tangential pressure distribution directly determines the range, morphology, and failure type of fracture opening.</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Plastic zone distribution map of cavitation holes with different radii: <bold>(a)</bold> 0.3 m, <bold>(b)</bold> 0.4 m, <bold>(c)</bold> 0.5 m, <bold>(d)</bold> 0.6 m, and <bold>(e)</bold> 0.7 m.</p>
</caption>
<graphic xlink:href="feart-14-1733210-g008.tif">
<alt-text content-type="machine-generated">Five vertically stacked simulation graphics labeled a to e display a circular central void with colored zones representing different mechanical states, explained in each panel's legend. As the sequence progresses, colored zones&#x2014;red, magenta, blue, and green&#x2014;expand outward from the central void, indicating changes in shear and tension states within a cyan background.</alt-text>
</graphic>
</fig>
<p>
<xref ref-type="fig" rid="F8">Figure 8a</xref> shows that the top and bottom of the cavitation hole exhibit tensile failure characteristics, while the two sides show shear failure. The essence of this phenomenon is the uneven distribution of tangential stress at different positions around the hole: at the top and bottom of the hole, tangential tensile stress dominates. When this stress exceeds the tensile strength <inline-formula id="inf27">
<mml:math id="m39">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> of coal, it directly drives the opening and extension of primary micro-fractures, forming tensile failure zones. On both sides of the hole, tangential compressive stress interacts with radial stress to generate shear stress, triggering the development of shear fractures. Although tangential compressive stress does not directly lead to fracture opening at this time, it provides mechanical conditions for the formation of shear fractures by changing the stress-field distribution. Note that the range of the plastic zone under this radius is relatively small, indicating that the influence depth of tangential stress is limited and the degree of fracture opening is weak, which is highly consistent with the subsequent result that the stress-release rate is only 46%.</p>
<p>A comparison of <xref ref-type="fig" rid="F8">Figures 8a,b</xref> with the increase in cavitation hole radius shows that the plastic zone around the hole expands significantly, and the plastic zone near ordinary boreholes also increases synchronously. The top of the cavitation hole changes from simple tensile failure to combined tension&#x2013;shear failure, while the two sides are still dominated by shear failure. The core driving factor for this change is the reconstruction of the tangential stress field: the increase in hole radius breaks the original stress balance around the hole in a wider range, further increases the peak value of tangential stress, and expands its influence radius outward. In the top area, the enhancement of tangential tensile stress not only intensifies the opening degree of existing fractures but also promotes the initiation of new tensile fractures. At the same time, the concentration of tangential compressive stress on the sides causes the interweaving of shear fractures and tensile fractures, forming combined tension&#x2013;shear failure zones. The increase in the plastic zone near ordinary boreholes is direct evidence of the outward propagation of the tangential stress field generated by the cavitation hole, proving that tangential stress can drive fractures to extend across boreholes and build connected channels for gas flow. At this time, the effect of tangential tensile stress of coal on fracture opening has initially shown advantages, corresponding to the increase of stress-release rate to 56% under this radius, which further verifies the positive correlation between tangential stress, fracture opening, and stress relief.</p>
<p>In <xref ref-type="fig" rid="F8">Figure 8c</xref>, the plastic zones at the haunch, shoulder, top, and bottom of the cavitation hole are significantly expanded compared with those when the radius is 0.4 m, while the plastic zones on both sides remain basically unchanged. This differentiated expansion characteristic is highly consistent with the distribution law of tangential stress: the haunch and shoulder positions are the secondary concentration areas of tangential stress. With the increase in hole radius, the tangential compressive stress in this area further accumulates, leading to the continuous increase of shear stress and extending the depth of shear fractures. Meanwhile, the coverage range of tangential tensile stress at the top and bottom expands, making tensile fractures form longer connected zones. From the perspective of mechanical mechanism, the gradient change of tangential stress promotes fractures to present a differentiated development mode of &#x201c;tensile opening at the top, shear extension at the shoulders.&#x201d; The reason why the plastic zones on both sides are stable is that their tangential stress is in a relatively balanced state without obvious increase or decrease, so the degree of fracture opening does not change significantly. The morphological change of the plastic zone under this radius also provides plastic mechanical support for the subsequent stress-release rate reaching 63%; that is, the expansion of the fracture opening range directly improves the efficiency of coal stress relief.</p>
<p>In <xref ref-type="fig" rid="F8">Figures 8d,e</xref>, the plastic zones at the haunch, shoulder, top, and bottom of the cavitation hole continue to expand, while the plastic zones on both sides still show no obvious changes, and the characteristic of combined tension&#x2013;shear failure becomes more prominent. This phenomenon indicates that when the cavitation hole radius exceeds 0.5 m, although the peak value of tangential stress continues to increase, the concentrated area of stress distribution does not change, remaining concentrated at the top, bottom, and shoulder positions. At this time, the tangential tensile stress is large enough to not only fully open the primary fractures but also promotes the connection of different fracture zones to form a network fracture system. The continuous action of tangential compressive stress on the sides refines the shear fractures and further improves the degree of coal fragmentation.</p>
<p>The plastic zone area corresponding to the radius of the cavitation hole is extracted, and the change in the plastic zone area is plotted against the radius of the cavitation hole. This is shown in <xref ref-type="fig" rid="F9">Figure 9</xref>, which indicates that the plastic zone area increases gradually with the radius of the cavitation hole. At &#x2265;0.5 m radius, the increase is approximately linear, indicating that the radius of the cavitation hole has a significant effect on the plastic zone area.</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>Graph of the plastic zone area versus the radius of the cavitation hole.</p>
</caption>
<graphic xlink:href="feart-14-1733210-g009.tif">
<alt-text content-type="machine-generated">Line graph showing the increase in plastic zone area, measured in meters, as cavitation hole radius, also in meters, increases from 0.3 to 0.7. Four data points are plotted with a curved line connecting them.</alt-text>
</graphic>
</fig>
<p>Based on the plot of the plastic zone area versus the radius of the cavitation hole, two more reasonable fits can be proposed for the curves modeling assumptions (<xref ref-type="disp-formula" rid="e13">Equations 13</xref>&#x2013;<xref ref-type="disp-formula" rid="e16">16</xref>):<list list-type="order">
<list-item>
<label>1.</label>
<p>Linear model. The model curve is</p>
</list-item>
</list>
<disp-formula id="e13">
<mml:math id="m40">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mi>R</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(13)</label>
</disp-formula>
</p>
<p>where <italic>a</italic>
<sub>1</sub> and <italic>b</italic>
<sub>1</sub> are parameters to be determined in relation to the plastic zone area.<list list-type="order">
<list-item>
<label>2.</label>
<p>3-function model. The model curve is</p>
</list-item>
</list>
<disp-formula id="e14">
<mml:math id="m41">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:msup>
<mml:mi>R</mml:mi>
<mml:mn>3</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:msup>
<mml:mi>R</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mi>R</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>d</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(14)</label>
</disp-formula>where <italic>a</italic>
<sub>2</sub>, <italic>b</italic>
<sub>2</sub>, <italic>c</italic>
<sub>2</sub>, and <italic>d</italic>
<sub>2</sub> are parameters to be determined in relation to the area of the pressure relief zone.</p>
<p>Based on the extracted values of the radius of the cavitation hole and the plastic zone surface, the fitted curves can be obtained, respectively,<disp-formula id="e15">
<mml:math id="m42">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>9.5</mml:mn>
<mml:mi>R</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>2.2</mml:mn>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(15)</label>
</disp-formula>
<disp-formula id="e16">
<mml:math id="m43">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>116.4</mml:mn>
<mml:msup>
<mml:mi>R</mml:mi>
<mml:mn>3</mml:mn>
</mml:msup>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>177.2</mml:mn>
<mml:msup>
<mml:mi>R</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>95.4</mml:mn>
<mml:mi>R</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>15.4</mml:mn>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(16)</label>
</disp-formula>
</p>
<p>The coefficients of determination of the fit are, respectively, <inline-formula id="inf28">
<mml:math id="m44">
<mml:mrow>
<mml:msubsup>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0.9636</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>; <inline-formula id="inf29">
<mml:math id="m45">
<mml:mrow>
<mml:msubsup>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0.9796</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<p>The above analysis reveals that the relationship between the plastic zone area of the cavitation hole and the radius of the cavitation hole is better fitted by the cubic function model, and the cubic function model can better describe the change of the plastic zone area of the cavitation hole with the radius of the cavitation hole.</p>
</sec>
</sec>
<sec id="s4-3">
<label>4.3</label>
<title>Comparison analysis of pressure relief and permeability enhancement effects of cavitation holes with different radii</title>
<p>The above analysis reveals that both the area of the pressure relief zone and the plastic zone area increase with the increase of the radius of the cavitation hole. It shows that the larger the radius of the cavitation hole is, the better the effects of pressure relief and permeability enhancement are. However, in engineering applications, the larger the radius of the cavitation hole, the higher the cost. To understand the relationship between the improved effect of the distressed zone and penetration enhancement and the cost of drilling a larger cavitation hole, we analyzed two aspects of the distressed zone area ratio and the plastic zone area ratio, which are defined in <xref ref-type="disp-formula" rid="e2">Equation 2</xref> and <xref ref-type="disp-formula" rid="e3">Equation 3</xref>, respectively.</p>
<p>According to <xref ref-type="disp-formula" rid="e2">Equation 2</xref>, pressure relief area ratio was calculated for the corresponding cavitation hole radius and plotted with the cavitation hole radius, as shown in <xref ref-type="fig" rid="F10">Figure 10</xref>, which reveals that the pressure relief area ratio increases and then decreases with the increase of the hole radius, and the maximum value of the pressure relief area ratio is between [0.4, 0.5] m for the hole radius.</p>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>Graph of the variation in the area ratio of the pressure relief zone with the radius of the cavitation hole.</p>
</caption>
<graphic xlink:href="feart-14-1733210-g010.tif">
<alt-text content-type="machine-generated">Line graph showing area ratio of pressure relief zone on the vertical axis and cavitation hole radius in meters on the horizontal axis. The plotted curve rises steeply, peaks near a radius of 0.45 meters at an area ratio of about 4.7, and then gradually declines.</alt-text>
</graphic>
</fig>
<p>The plastic zone area ratios for the corresponding cavitation hole radii are calculated according to <xref ref-type="disp-formula" rid="e3">Equation 3</xref> and plotted as a function of the variation of hole radius, as shown in <xref ref-type="fig" rid="F11">Figure 11</xref>, which reveals that the area ratio of the plastic zone increases and then decreases with the increase of the radius of the cavitation hole. The maximum value of the area ratio of the plastic zone is at the radius of the cavitation hole of approximately 0.4 m. The area ratio of the plastic zone increases and then decreases with the increase of the radius of the cavitation hole.</p>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>Graph of the variation in the plastic zone area ratio with the radius of the cavitation hole. The above analyses show that the optimal cavitation hole radius is between 0.4 m and 0.5 m.</p>
</caption>
<graphic xlink:href="feart-14-1733210-g011.tif">
<alt-text content-type="machine-generated">Line graph showing the relationship between cavitation hole radius in meters on the x-axis and plastic zone area ratio on the y-axis. Data points rise to a peak near 0.42 meters, then gradually decline and level off.</alt-text>
</graphic>
</fig>
</sec>
</sec>
<sec id="s5">
<label>5</label>
<title>Engineering applications</title>
<p>The No. 15 coal seam in Baiyangling Coal Mine (<xref ref-type="bibr" rid="B27">Yan, 2022</xref>) has an average thickness of 5 m, and the 15108 track drift area is a low-permeability and difficult-to-drain coal seam. The original gas content ranges from 4.52 m<sup>3</sup>/t to 6.5 m<sup>3</sup>/t, the inverse maximum gas pressure is 0.45 MPa, and the maximum initial gas desorption velocity &#x394;p is 38.5 mmHg. The maximum firmness coefficient f is 1.22, the coal mass failure types are mainly Class III and IV, and the permeability coefficient is 0.079 m<sup>2</sup>/(MPa<sup>2</sup>&#xb7;d). The cavitation hole radius has a significant impact on drainage effect and borehole stability: when the radius is 0.4 m, the permeability of the plastic zone around the hole reaches 300 times that of the elastic zone, the gas diffusion channel is unobstructed, and the gas pressure around the hole drops below 0.1 MPa after 150 days of drainage. The initial single-hole gas flow rate is nearly twice that of the 0.2 m radius, the stable drainage concentration is 21.49% (with 7 m borehole spacing), which is 3.36 times higher than that of non-cavitation holes (6.39%&#x2013;8.64%), and the flow attenuation coefficient is 0.015 d<sup>&#x2212;1</sup>, only 1/4.27 of that of non-cavitation holes (0.118 d<sup>&#x2212;1</sup>). For the 0.2 m radius, due to the limited plastic zone and insufficient permeability improvement, local high-gas areas still exist after 150 days of drainage. For the 0.6 m radius, the plastic zone doubles, the hole collapse risk is more than three times that of the 0.4 m radius, the borehole spacing must be reduced to 5 m, the construction cost increases by more than 30%, and additional investment in hole collapse control is required.</p>
<p>The optimal parameters are 0.4 m cavitation hole radius &#x2b; 7 m borehole spacing. At this time, there is no blank zone in the drainage area, and the total drainage volume is 1.41 times higher than that of the 5 m spacing. This practice verifies that when the cavitation hole radius is selected between 0.4 m and 0.5 m, it not only establishes a parameter benchmark for the drainage of low-permeability coal seams but also achieves the dual goals of drainage efficiency and borehole stability. Its applicability has been effectively verified under different coal quality conditions.</p>
</sec>
<sec sec-type="conclusion" id="s6">
<label>6</label>
<title>Conclusion</title>
<p>
<list list-type="order">
<list-item>
<p>Based on the stress distribution characteristics around the cavitation hole, three distinct stress zones are identified: the stress relief area, the stress concentration region, and the original rock stress area. Both the stress relief area and stress concentration region increase gas permeability within the coal mass, while the original rock stress area does not affect coal gas permeability.</p>
</list-item>
<list-item>
<p>The area of the stress relief zone and the area of the stress-plastic zone around the cavitation hole vary with the cavitation hole radius according to an exponential function model and a cubic function model, respectively. Their functional models are <inline-formula id="inf30">
<mml:math id="m46">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1.4</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>15</mml:mn>
<mml:mo>&#xd7;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>.</mml:mo>
<mml:msup>
<mml:mn>1</mml:mn>
<mml:mi>R</mml:mi>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula>; <inline-formula id="inf31">
<mml:math id="m47">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:mi>P</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>116.4</mml:mn>
<mml:msup>
<mml:mi>R</mml:mi>
<mml:mn>3</mml:mn>
</mml:msup>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>177.2</mml:mn>
<mml:msup>
<mml:mi>R</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>95.4</mml:mn>
<mml:mi>R</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>15.4</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
</list-item>
<list-item>
<p>Through FLAC3D numerical simulation and stress distribution mechanism analysis, the optimal cavitation hole radius for achieving peak pressure relief and permeability enhancement lies within the range of 0.4&#x2013;0.5 m. The scope determined in this study is based on scientific reasoning and argumentation, which ensures its reliability. Moreover, preliminary validation through numerical calculations supports this conclusion. Subsequent field applications will further verify the scientific nature and effectiveness of this scope.</p>
</list-item>
<list-item>
<p>The pressure relief effect of cavitation holes is closely related to the surrounding tangential stress distribution. When the lateral pressure coefficient (k &#x3c; 1), the tangential stress exhibits stress reduction along the B-B line and stress concentration along the A-A line, which correspond to the pressure relief effect and the formation of the plastic zone, respectively. Both effects can be quantified by the pressure relief zone area ratio and the plastic zone area ratio. When the radius of the cavitation hole ranges from 0.4 m to 0.5 m, the tangential stress change achieves the optimal effect, and the relevant indicators of the pressure relief zone and the plastic zone reach their peak values. If the radius is too small, the disturbance to the coal mass will be insufficient; if the radius is too large, the pressure relief efficiency will decrease, and the hole wall is prone to collapse. Meanwhile, the redistribution of tangential stress not only reduces the stress value but also changes the stress state of the coal mass. It promotes gas desorption in the pressure relief zone and unblocks the gas seepage channels in the stress concentration zone, synergistically improving the gas extraction efficiency.</p>
</list-item>
</list>
</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s7">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.</p>
</sec>
<sec sec-type="author-contributions" id="s8">
<title>Author contributions</title>
<p>HL: Writing &#x2013; original draft. JX: Writing &#x2013; review and editing. ZZ: Writing &#x2013; review and editing, Data curation. WL: Methodology, Writing &#x2013; review and editing. YF: Writing &#x2013; review and editing, Supervision. ZJ: Writing &#x2013; review and editing.</p>
</sec>
<ack>
<title>Acknowledgements</title>
<p>The authors thank all the athletes for their participation in this study.</p>
</ack>
<sec sec-type="COI-statement" id="s10">
<title>Conflict of interest</title>
<p>Authors HL, ZZ, WL, YF, and ZJ were employed by Zuoquan Wulinhou Mining Co., Ltd.</p>
<p>The remaining author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="ai-statement" id="s11">
<title>Generative AI statement</title>
<p>The author(s) declared that generative AI was not used in the creation of this manuscript.</p>
<p>Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.</p>
</sec>
<sec sec-type="disclaimer" id="s12">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Busch</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Gensterblum</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>CBM and CO<sub>2</sub>-ECBM related sorption processes in coal: a review</article-title>. <source>Int. J. Coal Geol.</source> <volume>87</volume> (<issue>2</issue>), <fpage>49</fpage>&#x2013;<lpage>71</lpage>. <pub-id pub-id-type="doi">10.1016/j.coal.2011.04.011</pub-id>
</mixed-citation>
</ref>
<ref id="B2">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Increasing permeability of coal seams using the phase energy of liquid carbon dioxide</article-title>. <source>J. CO<sub>2</sub> Util.</source> <volume>19</volume>, <fpage>112</fpage>&#x2013;<lpage>199</lpage>. <pub-id pub-id-type="doi">10.1016/j.jcou.2017.03.010</pub-id>
</mixed-citation>
</ref>
<ref id="B3">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Xie</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Qin</surname>
<given-names>B.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Increased permeability and coal and gas outburst prevention using hydraulic flushing technology with cross-seam borehole</article-title>. <source>J. Nat. Gas Sci. Eng.</source> <volume>73</volume>, <fpage>103067</fpage>. <pub-id pub-id-type="doi">10.1016/j.jngse.2019.103067</pub-id>
</mixed-citation>
</ref>
<ref id="B4">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Feng</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Experimental research on timeliness of promoting gas drainage by N<sub>2</sub> injection in coal bed</article-title>. <source>J. Northeast. Univ. Nat. Sci.</source> <volume>38</volume>, <fpage>26</fpage>&#x2013;<lpage>30</lpage>. <pub-id pub-id-type="doi">10.12068/j.issn.1005-3026.2017.07.024</pub-id>
</mixed-citation>
</ref>
<ref id="B5">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Lin</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Pang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Drilling large diameter cross-measure boreholes to improve gas drainage in highly gassy soft coal seams</article-title>. <source>J. Nat. Gas Sci. Eng.</source> <volume>26</volume>, <fpage>193</fpage>&#x2013;<lpage>204</lpage>. <pub-id pub-id-type="doi">10.1016/j.jngse.2015.05.035</pub-id>
</mixed-citation>
</ref>
<ref id="B6">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Guo</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Numerical analysis of the stress-unloading and permeability-inereasing characteristies and influencing factors of regional in seam borehole hydraulic cavitation</article-title>. <source>Jounal Xi&#x2019;an Univesity Sci. Technol.</source> <volume>43</volume> (<issue>2</issue>), <fpage>292</fpage>&#x2013;<lpage>300</lpage>. <pub-id pub-id-type="doi">10.13800/j.cnki.xakjdxxb.2023.0209</pub-id>
</mixed-citation>
</ref>
<ref id="B7">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>He</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Shen</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Ba</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Xiong</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Numerical research on the dynamic rock-breaking process of impact drilling with multi-nozzle water jets</article-title>. <source>J. Petroleum Sci. Eng.</source> <volume>207</volume>, <fpage>109145</fpage>. <pub-id pub-id-type="doi">10.1016/j.petrol.2021.109145</pub-id>
</mixed-citation>
</ref>
<ref id="B8">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Fan</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Aminian</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Optimal scheduling for enhanced coal bed methane production through CO<sub>2</sub> injection</article-title>. <source>Appl. Energy</source> <volume>113</volume>, <fpage>1475</fpage>&#x2013;<lpage>1483</lpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2013.08.074</pub-id>
</mixed-citation>
</ref>
<ref id="B9">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kong</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Feng</surname>
<given-names>J.</given-names>
</name>
<etal/>
</person-group> (<year>2016</year>). <article-title>Coupled analysis about multi-factors to the effective influence radius of hydraulic flushing: application of response surface methodology</article-title>. <source>J. Nat. Gas Sci. Eng.</source> <volume>32</volume>, <fpage>538</fpage>&#x2013;<lpage>548</lpage>. <pub-id pub-id-type="doi">10.1016/j.jngse.2016.04.043</pub-id>
</mixed-citation>
</ref>
<ref id="B10">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Mining thin sub-layer as self-protective coal seam to reduce the danger of coal and gas outburst</article-title>. <source>Nat. Hazards</source> <volume>71</volume> (<issue>1</issue>), <fpage>41</fpage>&#x2013;<lpage>52</lpage>. <pub-id pub-id-type="doi">10.1007/s11069-013-0898-1</pub-id>
</mixed-citation>
</ref>
<ref id="B11">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Fu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>Experimental study on compressive behavior and failure characteristics of imitation steel fiber concrete under uniaxial load</article-title>. <source>Constr. Build. Mater.</source> <volume>399</volume>, <fpage>132599</fpage>. <pub-id pub-id-type="doi">10.1016/j.conbuildmat.2023.132599</pub-id>
</mixed-citation>
</ref>
<ref id="B12">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Zhai</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
<etal/>
</person-group> (<year>2024</year>). <article-title>Roadway portal and self-moving hydraulic support for rockburst prevention in coal mine and its application</article-title>. <source>Phys. Fluids</source> <volume>36</volume>, <fpage>124136</fpage>. <pub-id pub-id-type="doi">10.1063/5.0243798</pub-id>
</mixed-citation>
</ref>
<ref id="B13">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Cao</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2025</year>). <article-title>Characteristics and trends of coal mine safety development</article-title>. <source>Energy Sources, Part A Recovery, Util. Environ. Eff.</source> <volume>47</volume>, <fpage>2316</fpage>&#x2013;<lpage>2334</lpage>. <pub-id pub-id-type="doi">10.1080/15567036.2020.1852339</pub-id>
</mixed-citation>
</ref>
<ref id="B14">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lin</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zou</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>C.</given-names>
</name>
<etal/>
</person-group> (<year>2015</year>). <article-title>Cross-borehole hydraulic slotting technique for preventing and controlling coal and gas outbursts during coal roadway excavation</article-title>. <source>J. Nat. Gas Sci. Eng.</source> <volume>26</volume>, <fpage>518</fpage>&#x2013;<lpage>525</lpage>. <pub-id pub-id-type="doi">10.1016/j.jngse.2015.06.035</pub-id>
</mixed-citation>
</ref>
<ref id="B15">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Elastoplastic stress and displacement analyticalsolutions to deep-buried circular tunnels consideringintermediate principal stress and dilatancy</article-title>. <source>Eng. Mech.</source> <volume>29</volume> (<issue>08</issue>), <fpage>289</fpage>&#x2013;<lpage>296</lpage>. <pub-id pub-id-type="doi">10.6052/j.issn.1000-4750.2011.08.0511</pub-id>
</mixed-citation>
</ref>
<ref id="B16">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>D.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>Experimental study of effect of liquid nitrogen cold soaking on coal pore structure and fractal characteristics</article-title>. <source>Energy</source> <volume>275</volume>, <fpage>127470</fpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2023.127470</pub-id>
</mixed-citation>
</ref>
<ref id="B17">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2025</year>). <article-title>Investigations on the mechanism of the microstructural evolution of different coal ranks under liquid nitrogen cold soaking</article-title>. <source>Energy Sources, Part A Recovery, Util. Environ. Eff.</source> <volume>47</volume> (<issue>1</issue>), <fpage>2595</fpage>&#x2013;<lpage>2612</lpage>. <pub-id pub-id-type="doi">10.1080/15567036.2020.1841856</pub-id>
</mixed-citation>
</ref>
<ref id="B18">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lou</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Kang</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Wei</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Yue</surname>
<given-names>J.</given-names>
</name>
<etal/>
</person-group> (<year>2024</year>). <article-title>Plugging methods for underground gas extraction boreholes in coal seams: a review of processes, challenges and strategies</article-title>. <source>Gas Sci. Eng.</source> <volume>122</volume>, <fpage>205225</fpage>. <pub-id pub-id-type="doi">10.1016/j.jgsce.2024.205225</pub-id>
</mixed-citation>
</ref>
<ref id="B19">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lou</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Yao</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Qin</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>Z.</given-names>
</name>
<etal/>
</person-group> (<year>2025</year>). <article-title>A novel dynamic filling material for plugging fractures around underground gas extraction boreholes: experimental and engineering performances</article-title>. <source>Energy</source> <volume>314</volume>, <fpage>134202</fpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2024.134202</pub-id>
</mixed-citation>
</ref>
<ref id="B20">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lu</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Improving the gate road development rate and reducing outburst occurrences using the waterjet technique in high gas content outburst-prone soft coal seam</article-title>. <source>Int. J. Rock Mech. Min. Sci.</source> <volume>48</volume> (<issue>8</issue>), <fpage>1271</fpage>&#x2013;<lpage>1282</lpage>. <pub-id pub-id-type="doi">10.1016/j.ijrmms.2011.09.003</pub-id>
</mixed-citation>
</ref>
<ref id="B21">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Ge</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Xia</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Progress on the hydraulic measures for grid slotting and fracking to enhance coal seam permeability</article-title>. <source>Int. J. Min. Sci. Technol.</source> <volume>27</volume> (<issue>5</issue>), <fpage>867</fpage>&#x2013;<lpage>871</lpage>. <pub-id pub-id-type="doi">10.1016/j.ijmst.2017.07.011</pub-id>
</mixed-citation>
</ref>
<ref id="B22">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ma</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Theoretical analysis of water inrush from goaf in floor and mining stability for upward repeated mining</article-title>. <source>J. China Coal Soc.</source> <volume>48</volume>, <fpage>37</fpage>&#x2013;<lpage>46</lpage>. <pub-id pub-id-type="doi">10.13225/j.cnki.jccs.2022.0109</pub-id>
</mixed-citation>
</ref>
<ref id="B23">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Liao</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Chu</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Study on mechanism of permeability improvement for gas drainage by cross-seam cavitation borehole</article-title>. <source>Coal Sci. Technol.</source> <volume>49</volume> (<issue>5</issue>), <fpage>75</fpage>&#x2013;<lpage>82</lpage>. <pub-id pub-id-type="doi">10.13199/j.cnki.cst.2021.05.010</pub-id>
</mixed-citation>
</ref>
<ref id="B24">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Analysis of trapezoid fracture parameters and fracture evolution of overburden in shallow coal seam mining</article-title>. <source>J. Henan Polytechnic University(natural Science)</source> <volume>43</volume> (<issue>3</issue>), <fpage>50</fpage>&#x2013;<lpage>59</lpage>. <pub-id pub-id-type="doi">10.16186/j.cnki.1673-9787.2023060054</pub-id>
</mixed-citation>
</ref>
<ref id="B25">
<mixed-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Xu</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2016</year>). <source>Elasticity</source>. <edition>5th edition</edition>. <publisher-loc>Beijing</publisher-loc>: <publisher-name>Higher Education Press</publisher-name>.</mixed-citation>
</ref>
<ref id="B26">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xue</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Mechanism of deformation-induced damage and gas permeability enhancement of coal under typical mining layouts</article-title>. <source>Chin. J. Geotechnical Eng.</source> <volume>35</volume> (<issue>2</issue>), <fpage>328</fpage>&#x2013;<lpage>336</lpage>. <pub-id pub-id-type="doi">10.11779/CJGE201302016</pub-id>
</mixed-citation>
</ref>
<ref id="B27">
<mixed-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Yan</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2022</year>). <source>Study on technology of increasing flow andpromoting pumping of water jet cavitation in parallel drilling along layer</source>. <publisher-name>Henan Polytechnic University</publisher-name>. <pub-id pub-id-type="doi">10.13301/j.cnki.ct.2023.01.037</pub-id>
</mixed-citation>
</ref>
<ref id="B28">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Du</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Jin</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Gong</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>N<sub>2</sub> injection to enhance coal seam gas drainage (N<sub>2</sub>-ECGD): insights from underground field trial investigation</article-title>. <source>Energy</source> <volume>239</volume>, <fpage>122247</fpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2021.122247</pub-id>
</mixed-citation>
</ref>
<ref id="B29">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Cheng</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Yuan</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Dong</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>L.</given-names>
</name>
<etal/>
</person-group> (<year>2017</year>). <article-title>A novel in-seam borehole hydraulic flushing gas extraction technology in the heading face: enhanced permeability mechanism, gas flow characteristics, and application</article-title>. <source>J. Nat. Gas Sci. Eng.</source> <volume>46</volume>, <fpage>498</fpage>&#x2013;<lpage>514</lpage>. <pub-id pub-id-type="doi">10.1016/j.jngse.2017.08.022</pub-id>
</mixed-citation>
</ref>
<ref id="B30">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Cheng</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Yuan</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Pan</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Hydraulic flushing in soft coal sublayer: gas extraction enhancement mechanism and field application</article-title>. <source>Energy Sci. and Eng.</source> <volume>7</volume> (<issue>5</issue>), <fpage>1970</fpage>&#x2013;<lpage>1973</lpage>. <pub-id pub-id-type="doi">10.1002/ese3.405</pub-id>
</mixed-citation>
</ref>
<ref id="B31">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Cheng</surname>
<given-names>Y.-P.</given-names>
</name>
<name>
<surname>Yuan</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Enhancement of gas drainage efficiency in a special thick coal seam through hydraulic flushing</article-title>. <source>Int. J. Rock Mech. Min. Sci.</source> <volume>124</volume>, <fpage>104085</fpage>. <pub-id pub-id-type="doi">10.1016/j.ijrmms.2019.104085</pub-id>
</mixed-citation>
</ref>
<ref id="B32">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Kong</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Peng</surname>
<given-names>S.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>Laboratory experiments of CO<sub>2</sub>-enhanced coalbed methane recovery considering CO<sub>2</sub> sequestration in a coal seam</article-title>. <source>Energy</source> <volume>262</volume>, <fpage>125473</fpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2022.125473</pub-id>
</mixed-citation>
</ref>
<ref id="B33">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhou</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Elsworth</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Fan</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Airflow disturbance induced by coal mine outburst shock waves: a case study of a gas outburst disaster in China</article-title>. <source>Int. J. Rock Mech. Min. Sci.</source> <volume>128</volume>, <fpage>104262</fpage>. <pub-id pub-id-type="doi">10.1016/j.ijrmms.2020.104262</pub-id>
</mixed-citation>
</ref>
<ref id="B34">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhou</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Bai</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Luo</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Effect of damage zone around borehole on carbon dioxide injection promoted gas extraction in soft and low-permeability coal seam</article-title>. <source>Front. Earth Sci.</source> <volume>17</volume> (<issue>3</issue>), <fpage>776</fpage>&#x2013;<lpage>787</lpage>. <pub-id pub-id-type="doi">10.1007/s11707-022-1036-8</pub-id>
</mixed-citation>
</ref>
</ref-list>
<fn-group>
<fn fn-type="custom" custom-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1414047/overview">Xiaohu Yang</ext-link>, Xi&#x2019;an Jiaotong University, China</p>
</fn>
<fn fn-type="custom" custom-type="reviewed-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1983084/overview">Haowei Yao</ext-link>, Zhengzhou University of Light Industry, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/3271557/overview">Zhiyuan Shen</ext-link>, Shenyang Jianzhu University, China</p>
</fn>
</fn-group>
</back>
</article>