<?xml version="1.0" encoding="us-ascii"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Earth Sci.</journal-id>
<journal-title>Frontiers in Earth Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Earth Sci.</abbrev-journal-title>
<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">1401455</article-id>
<article-id pub-id-type="doi">10.3389/feart.2024.1401455</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Earth Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Production characteristics and influencing factors of coalbed methane wells: a case study of the high-ranking coal seam in the southeastern Qinshui Basin, China</article-title>
<alt-title alt-title-type="left-running-head">Chen 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.2024.1401455">10.3389/feart.2024.1401455</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Chen</surname>
<given-names>Xiaolong</given-names>
</name>
<uri xlink:href="https://loop.frontiersin.org/people/2686908/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Gao</surname>
<given-names>Yufei</given-names>
</name>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Wang</surname>
<given-names>Yaqing</given-names>
</name>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
</contrib-group>
<aff>
<institution>China National Offshore Oil Corporation Research Institute Co., Ltd.</institution>, <addr-line>Beijing</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1447721/overview">Hongjian Zhu</ext-link>, Yanshan University, China</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2399920/overview">Qitao Zhang</ext-link>, The Pennsylvania State University (PSU), United States</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2677347/overview">Yanwei Wang</ext-link>, China University of Geosciences, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Yaqing Wang, <email>644646437@qq.com</email>
</corresp>
</author-notes>
<pub-date pub-type="epub">
<day>22</day>
<month>05</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>12</volume>
<elocation-id>1401455</elocation-id>
<history>
<date date-type="received">
<day>15</day>
<month>03</month>
<year>2024</year>
</date>
<date date-type="accepted">
<day>26</day>
<month>04</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 Chen, Gao and Wang.</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Chen, Gao and Wang</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). 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.</p>
</license>
</permissions>
<abstract>
<p>This study focuses on coalbed methane (CBM) wells in high-ranking coal seam as the research subject. Considering the influence of effective stress and matrix shrinkage, a comprehensive permeability calculation model for CBM reservoirs is established. Based on this model, the variations in pressure and permeability during well production are quantified. By integrating static geological parameters, a finely classified classification of CBM wells is achieved using self-organizing map (SOM) neural network. Subsequently, an analysis of production dynamic characteristics and productivity differences among different types of CBM wells is performed, followed by providing drainage optimization suggestions. The results of SOM analysis show that 7,000 m<sup>3</sup>/d and 1,500 m<sup>3</sup>/d can be used as the production boundaries for the wells with different productivity in Block P. The daily gas production of exceptional well exceeds 7,000 m<sup>3</sup>/d, and the permeability remains relatively stable throughout the drainage process of this well. The daily gas production of the potential well ranges from 1,500 to 7,000 m<sup>3</sup>/d, and the permeability exhibits a significant decrease during the drainage process. The daily gas production of Inefficient well is consistently below 1,500 m<sup>3</sup>/d with moderate permeability variation. In addition to well location and structural geology, production variability is also influenced by the matching of reservoir conditions and drainage systems. This is primarily manifested in discontinuous drainage systems and rapid decline in bottom hole pressure (BHP) during early production. The analysis of drainage parameters indicates that in order to achieve optimal production from CBM wells, the BHP should exhibit an initial rapid decline followed by a slowly decrease during the early production period, with an average pressure drop ranging from 0.005 to 0.02 MPa/d. The research findings can offer technical guidance for the future advancement of CBM in the P Block.</p>
</abstract>
<kwd-group>
<kwd>coalbed methane</kwd>
<kwd>permeability</kwd>
<kwd>production characteristics</kwd>
<kwd>high-ranking coal seam</kwd>
<kwd>Qinshui Basin</kwd>
</kwd-group>
<contract-sponsor id="cn001">China Postdoctoral Science Foundation<named-content content-type="fundref-id">10.13039/501100002858</named-content>
</contract-sponsor>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Geoscience and Society</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>The utilization of clean energy is regarded as one of the efficacious measures to alleviate global climate change (<xref ref-type="bibr" rid="B7">Dickey, 2006</xref>; <xref ref-type="bibr" rid="B42">Zhu et al., 2021</xref>). CBM, characterized by its non-toxic combustion, absence of particulate matter emissions, and lower carbon dioxide release compared to coal, oil, or wood, is acknowledged as a clean fuel (<xref ref-type="bibr" rid="B5">Chen et al., 2017</xref>; <xref ref-type="bibr" rid="B29">Peng et al., 2017</xref>; <xref ref-type="bibr" rid="B16">Kong et al., 2022</xref>; <xref ref-type="bibr" rid="B25">Ni et al., 2023</xref>). Over the past 2 decades, CBM has emerged as a significant energy resource and is anticipated to play a pivotal role in meeting future global energy demands (<xref ref-type="bibr" rid="B24">Moore, 2012</xref>; <xref ref-type="bibr" rid="B22">Lou et al., 2013</xref>; <xref ref-type="bibr" rid="B8">Fu et al., 2017</xref>; <xref ref-type="bibr" rid="B11">Gao et al., 2023a</xref>; <xref ref-type="bibr" rid="B13">Ifrene et al., 2023</xref>; <xref ref-type="bibr" rid="B38">Wang et al., 2024a</xref>). China, as the current global leader in coal consumption and production, can benefit from the development of CBM resources. This not only aids in reducing natural greenhouse gas emissions and maximizing energy utilization but also enhances safety measures in coal mining operations by mitigating risks associated with underground disasters such as fires caused by spontaneous combustion of natural gas and coal (<xref ref-type="bibr" rid="B1">Bao et al., 2020</xref>; <xref ref-type="bibr" rid="B17">Kong et al., 2021</xref>). The CBM reserves in China are distributed among high-rank, middle-rank, and low-rank coal seams, with each accounting for approximately one-third of the total resources. However, the focus of CBM development is primarily on the high-rank and middle-rank seams. A prime example of a high-rank CBM field is located in the southern Qinshui Basin of Shanxi Province (<xref ref-type="bibr" rid="B2">Bustin, 1997</xref>; <xref ref-type="bibr" rid="B4">Chen et al., 2016</xref>; <xref ref-type="bibr" rid="B10">Gao et al., 2023b</xref>; <xref ref-type="bibr" rid="B31">Qiu et al., 2023</xref>). In recent years, scholars have conducted comprehensive research on the geological background, physical properties of CBM reservoirs, resource evaluation, reservoir formation mechanisms, and production potential in the Qinshui Basin. Li et al. demonstrated that the tectonic subsidence and weak hydrodynamics in the southern block of the Qinshui Basin created favorable constraints for the occurrence of CBM (<xref ref-type="bibr" rid="B20">Li et al., 2018</xref>). Liu et al. discovered that during the formation of CBM reservoirs, the spatial distribution of CBM composition exhibits a gradual change in different locations, and there is relatively weak interference between CBM wells in the southern block of Qinshui Basin (<xref ref-type="bibr" rid="B21">Liu et al., 2018</xref>). Lu et al. conducted a study on the production characteristics and key factors of numerous gas Wells in the southern Qinshui Basin using variable correlation analysis and grey correlation theory. They concluded that hydraulic fracturing is an effective method to ensure high production of CBM wells, and identified gas content and permeability as two crucial factors influencing gas production of CBM wells (<xref ref-type="bibr" rid="B23">Lv et al., 2012</xref>). Jiang et al. conducted a study on the characteristics and gas production potential of the No. 15 coal seam in the northern block of Qinshui Basin, revealing that the No. 15 coal seam exhibits significant thickness, high adsorption capacity, moderate gas content, low critical desorption pressure, and substantial gas production potential (<xref ref-type="bibr" rid="B14">Jiang et al., 2022</xref>). At present, the commercial development of CBM from high-rank coal reservoirs in Qinshui Basin has achieved remarkable results (<xref ref-type="bibr" rid="B26">Nie et al., 2000</xref>; <xref ref-type="bibr" rid="B30">Qin et al., 2005</xref>; <xref ref-type="bibr" rid="B6">Chen et al., 2015</xref>; <xref ref-type="bibr" rid="B37">Wang et al., 2024b</xref>).</p>
<p>The P block is situated in the southeastern region of the Qinshui Basin, where CBM exploration activities commenced in the 1970s. This area represents a typical site for CBM development in China, with two primary target coal seams: the No. 3 coal seam of the Shanxi Formation and the No. 15 coal seam of the Taiyuan Formation. The CBM well of No. 3 coal in Block P has a prolonged production history, and the overall gas production of this coal seam has exhibited a noticeable declining trend. Therefore, it is imperative to develop No. 15 coal seam in order to enhance the gas production of the block (<xref ref-type="bibr" rid="B41">Yuanyuan et al., 2012</xref>). However, due to the different reservoir characteristics, there are significant differences in the gas production rules of No. 3 and No. 15 coal seam CBM wells, and the development method of No. 15 coal seam cannot be fully used for reference by No. 3 coal seam. Therefore, it is the key to develop efficient development methods to clarify the production characteristics and influencing factors of CBM wells in No. 15 coal seam.</p>
<p>Permeability is a key parameter that determines fluid migration ability, and its influence on the production dynamics of CBM wells cannot be ignored. In the process of CBM production, the permeability of the reservoir will change significantly due to the strong influence of effective stress, matrix contraction and gas slippage (<xref ref-type="bibr" rid="B3">Cai et al., 2014</xref>). Scholars have established multiple permeability models to evaluate the dynamic changes of reservoir permeability. These models can be divided into three categories: stress-based models (<xref ref-type="bibr" rid="B27">Palmer and Mansoori, 1998a</xref>), strain-based models (<xref ref-type="bibr" rid="B33">Shi and Durucan, 2005a</xref>) and production data-based models (<xref ref-type="bibr" rid="B19">Lai et al., 2013</xref>). However, these models consider fewer factors and have limited applicability. Fine classification of CBM wells is the prerequisite for accurate analysis of production dynamics. Most previous studies have classified wells based on geostatic parameters, without quantifying the influence of dynamic characteristic parameters (dynamic changes in permeability and pressure) on well classification. Shang et al. classified CBM wells into four types according to gas production and analyzed the forming reasons. The results showed that coal structure was the basic factor affecting the productivity of CBM wells, and gas content and permeability were the direct geological factors controlling the production of CBM wells (<xref ref-type="bibr" rid="B32">Sang et al., 2009</xref>). Yi et al. classified CBM wells according to the type of reservoir pressure drop, and divided them into three types: rapid drop type, medium-term stable type and slow drop type. Among them, the rapid drop type is more favorable for high and stable production of CBM (<xref ref-type="bibr" rid="B40">Yi et al., 2019</xref>). Hou et al. classified CBM wells into four categories based on the static coal reservoir properties and the dynamic average daily gas production data of a single well, and systematically analyzed the gas production characteristics of each type of well (<xref ref-type="bibr" rid="B12">Hou et al., 2018</xref>). At present, most of the studies on CBM well classification do not consider dynamic parameters such as permeability and average reservoir pressure. In addition, the classification results of CBM wells are greatly affected by subjective factors. Therefore, it is the key point of the future CBM development evaluation in Block P to objectively classify the CBM wells according to the dynamic and static parameters, and comprehensively analyze the production characteristics and control factors of different types of CBM wells.</p>
<p>This study focuses on the No. 15 coal seam in the P block as a research subject, investigating the impact of various static factors on CBM production. Based on this, a neural network method is employed to finely classify CBM wells by considering both static and dynamic parameters. The production dynamics and productivity differences among different types of CBM wells are analyzed, providing CBM drainage optimization suggestions that offer technical support for future development of No. 15 coal seam.</p>
</sec>
<sec id="s2">
<title>2 Regional geology and reservoir characteristics</title>
<p>The P block is situated in the southeastern part of the Qinshui Basin. The structural configuration within the block is characterized by simplicity, presenting a predominantly west-dipping monoclinic structure overall, with only two faults observed in the northeastern region (refer to <xref ref-type="fig" rid="F1">Figure 1</xref>). Secondary folds are primarily developed within the block, exhibiting a general orientation along the north-south axis on a planar scale. These folds possess distinctive features including broad and gentle characteristics, symmetrical wings, and shallow dip angles. The fold structures in the east-west direction consist of Zhengcun anticline, Panhe syncline, Shigou anticline, Huojiashan syncline, Mashancun anticline, etc. The strata exhibit gentle inclination angles ranging from 5&#xb0; to 8&#xb0; on average, with a maximum dip angle of only 15&#xb0;. Generally speaking, the synclines tend to have relatively high gas content compared to the anticlines.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Location of the P block in China.</p>
</caption>
<graphic xlink:href="feart-12-1401455-g001.tif"/>
</fig>
<p>In the P block, the primary coal-bearing strata comprise the Upper Carboniferous Taiyuan Formation and the Lower Permian Shanxi Formation (refer to <xref ref-type="fig" rid="F2">Figure 2</xref>). The Taiyuan Formation predominantly accumulated in a marine sedimentary environment, consisting of sandstone, mudstone, coal, and limestone. Conversely, the Shanxi Formation primarily deposited in a terrestrial sedimentary setting composed of sandstone, mudstone, and coal. The coal seam within the block comprises the No. 3 coal seam of Shanxi Formation and the No. 15 coal seam of Taiyuan Formation, with a current focus on exploring and developing CBM in the latter. The No. 15 coal seam consists of sediments from barrier coast and delta systems. Furthermore, within the paleo-sedimentary environment, there is well-developed distributary channeling and inter-distributary bay formation, leading to bifurcation or even disappearance of certain coal seams within the block.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Stratigraphic column of the Permo-Carboniferous coal-bearing strata in P block.</p>
</caption>
<graphic xlink:href="feart-12-1401455-g002.tif"/>
</fig>
<p>The No. 15 coal seam is buried at depths ranging from 350 m to 700 m, with an average burial depth of 510 m. The thickness of the coal seam remains consistent in the horizontal plane, varying between 2.68 m and 7.78 m, with an average thickness of 4.14 m. The coal body structure primarily consists of intact structural coal and fragmented coal, exhibiting a permeability range of 0.08 to 5.707 mD, with an average permeability value of 1.84 mD. Additionally, the raw coal exhibits an ash yield ranging from 9.15% to 15.54%, averaging at around 12.17%, indicating its classification as low and medium ash coal. The desorption test results indicate that No. 15 coal exhibits a high gas content, ranging from 12.57 to 27.71 m<sup>3</sup>/t, with an average of 20.91 m<sup>3</sup>/t, and the gas saturation exceeds 90%. Additionally, the reservoir pressure ranges from 1.82 to 2.93 MPa, averaging at 2.21 MPa. The overall development geological conditions of No. 15 coal are more suitable for the development of CBM.</p>
</sec>
<sec id="s3">
<title>3 Mathematical model of dynamic permeability of CBM well</title>
<p>The consideration of permeability as a parameter in dynamic analysis is crucial; however, it is impractical to measure the permeability of CBM wells in real-time during production. Therefore, this section initially derives the calculation model for determining the dynamic permeability of CBM wells. In the actual production process of CBM, reservoir pressure is closely correlated with permeability. Therefore, the dynamic characteristics of coal reservoir pressure are investigated through analysis of the material balance equation. The material balance equation for CBM can be expressed as follows:<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>A</mml:mi>
<mml:mi>h</mml:mi>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>&#x3c1;</mml:mi>
<mml:mi>B</mml:mi>
</mml:msub>
<mml:mi>A</mml:mi>
<mml:mi>h</mml:mi>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:mfrac>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>A</mml:mi>
<mml:mi>h</mml:mi>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>g</mml:mi>
</mml:msub>
</mml:mfrac>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3c1;</mml:mi>
<mml:mi>B</mml:mi>
</mml:msub>
<mml:mi>A</mml:mi>
<mml:mi>h</mml:mi>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:mfrac>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>where, <italic>G</italic>
<sub>
<italic>p</italic>
</sub> is cumulative gas production, m<sup>3</sup>; <italic>A</italic> is well control area, m<sup>2</sup>; <italic>h</italic> is coal seam thickness, m; <inline-formula id="inf1">
<mml:math id="m2">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is initial porosity, %; <italic>S</italic>
<sub>
<italic>wi</italic>
</sub> is initial water saturation, %; <italic>B</italic>
<sub>
<italic>gi</italic>
</sub> is initial gas volume coefficient, m<sup>3</sup>/m<sup>3</sup>; <inline-formula id="inf2">
<mml:math id="m3">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c1;</mml:mi>
<mml:mi>B</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is coal density, kg/m<sup>3</sup>; <italic>V</italic>
<sub>
<italic>L</italic>
</sub> is langmuir volume, m<sup>3</sup>/t; <italic>p</italic>
<sub>
<italic>i</italic>
</sub> is initial reservoir pressure, MPa; <italic>p</italic>
<sub>
<italic>L</italic>
</sub> is langmuir pressure, MPa; <inline-formula id="inf3">
<mml:math id="m4">
<mml:mrow>
<mml:mi>&#x3d5;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is porosity, %; <italic>S</italic>
<sub>
<italic>w</italic>
</sub> is water saturation, %; <italic>B</italic>
<sub>
<italic>g</italic>
</sub> is gas volume coefficient, m<sup>3</sup>/m<sup>3</sup>; <italic>p</italic> is reservoir pressure, MPa.</p>
<p>The material balance equation for water can be expressed as follows:<disp-formula id="e2">
<mml:math id="m5">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>A</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>&#x3d5;</mml:mi>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>A</mml:mi>
<mml:mi>h</mml:mi>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>where, <italic>W</italic>
<sub>
<italic>B</italic>
</sub> is gas volume coefficient, m<sup>3</sup>/m<sup>3</sup>; c<sub>w</sub> is water compression coefficient, MPa<sup>&#x2212;1</sup>; <italic>W</italic>
<sub>
<italic>e</italic>
</sub> is water invasion volume, m<sup>3</sup>; <italic>W</italic>
<sub>
<italic>p</italic>
</sub> is cumulative water production, m<sup>3</sup>.</p>
<p>W<sub>e</sub> &#x3d; 0, Eq. <xref ref-type="disp-formula" rid="e2">2</xref> is simplified as:<disp-formula id="e3">
<mml:math id="m6">
<mml:mrow>
<mml:mi>A</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>&#x3d5;</mml:mi>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>A</mml:mi>
<mml:mi>h</mml:mi>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
</p>
<p>The porosity of coal exhibits a self-adjusting effect, which varies in response to the pressure exerted on the coal reservoir. According to the P&#x26;M model (<xref ref-type="bibr" rid="B34">Shi et al., 2021</xref>), at the current reservoir pressure, the ratio of porosity to initial porosity can be expressed as follows:<disp-formula id="e4">
<mml:math id="m7">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x3d5;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b5;</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>K</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>M</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>where, <italic>c</italic>
<sub>
<italic>m</italic>
</sub> is P&#x26;M model coefficient, MPa<sup>&#x2212;1</sup>; <italic>&#x3b5;</italic>
<sub>
<italic>l</italic>
</sub> is langmuir volumetric strain fitting parameter; <italic>K</italic> is volume modulus, MPa; <italic>M</italic> is unaxial modulus, MPa.<disp-formula id="e5">
<mml:math id="m8">
<mml:mrow>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>g</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>M</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>K</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>M</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>&#x3b3;</mml:mi>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>where, <italic>g</italic> is a correction factor between 0 and 1; <inline-formula id="inf4">
<mml:math id="m9">
<mml:mrow>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is a number factor between 0 and 1; <italic>&#x3b3;</italic> is coal matrix compression coefficient, MPa<sup>&#x2212;1</sup>.<disp-formula id="e6">
<mml:math id="m10">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>M</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>E</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>v</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>
<disp-formula id="e7">
<mml:math id="m11">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>K</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>M</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</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:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>v</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>where, <italic>E</italic> is young&#x2019;s modulus, MPa; <italic>&#x3bd;</italic> is poisson ratio.</p>
<p>Combine Eqs <xref ref-type="disp-formula" rid="e4">4</xref>&#x2013;<xref ref-type="disp-formula" rid="e7">7</xref> to obtain the equation of porosity under the influence of effective stress and matrix shrinkage:<disp-formula id="e8">
<mml:math id="m12">
<mml:mrow>
<mml:mtable class="align" columnalign="left">
<mml:mtr>
<mml:mtd columnalign="right">
<mml:mrow>
<mml:mi>&#x3d5;</mml:mi>
</mml:mrow>
<mml:mspace width="0.3em"/>
</mml:mtd>
<mml:mtd columnalign="left">
<mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mspace width="0.3em"/>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="}" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>v</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>3</mml:mn>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>&#x3b3;</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="right"/>
<mml:mtd columnalign="left">
<mml:mspace width="1em"/>
<mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mspace width="0.2em"/>
<mml:msub>
<mml:mi>&#x3b5;</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>v</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>3</mml:mn>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>
</p>
<p>The iterative solution of Eqs <xref ref-type="disp-formula" rid="e1">1</xref>, <xref ref-type="disp-formula" rid="e3">3</xref>, <xref ref-type="disp-formula" rid="e8">8</xref> enables the calculation of coal reservoir pressure during the production process, thereby providing fundamental parameters for subsequent permeability calculations.</p>
<p>For an unsaturated coal reservoir, when the reservoir pressure exceeds the critical desorption pressure, CBM cannot be released from the coal matrix, resulting in pore saturation with coalbed water. This stage is characterized by single-phase water flow production. During this stage, the permeability of the coal reservoir is solely influenced by effective stress. The mathematical model describing changes in coal permeability during this stage is as follows (<xref ref-type="bibr" rid="B6">Chen et al., 2015</xref>):<disp-formula id="e9">
<mml:math id="m13">
<mml:mrow>
<mml:mi>K</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi>I</mml:mi>
</mml:msub>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>f</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>v</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>where, <italic>K</italic>
<sub>
<italic>i</italic>
</sub> is initial permeability, mD; <italic>K</italic> is dynamic permeability of single-phase flow stage, mD.</p>
<p>When the reservoir pressure is below the critical desorption pressure, CBM is released from the coal matrix due to the decrease in pressure. The permeability of the reservoir is not only affected by effective stress but also by matrix shrinkage. According to the classical P&#x26;M model (<xref ref-type="bibr" rid="B34">Shi et al., 2021</xref>), the relationship between porosity and permeability can be expressed as follows:<disp-formula id="e10">
<mml:math id="m14">
<mml:mrow>
<mml:mi>K</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msup>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x3d5;</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>P</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>3</mml:mn>
</mml:msup>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>where, <italic>K</italic>
<sub>
<italic>cd</italic>
</sub> is the permeability at critical desorption pressure, mD; <inline-formula id="inf5">
<mml:math id="m15">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the prosity at critical desorption pressure</p>
<p>Combine Eqs <xref ref-type="disp-formula" rid="e8">8</xref>, <xref ref-type="disp-formula" rid="e10">10</xref> to obtain the equation of permeability under the influence of effective stress and matrix shrinkage:<disp-formula id="e11">
<mml:math id="m16">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>w</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msup>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="}" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>v</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>3</mml:mn>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x3b3;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b5;</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>v</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>3</mml:mn>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>3</mml:mn>
</mml:msup>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>where, <italic>K</italic>
<sub>
<italic>now</italic>
</sub> is the permeability below critical desorption pressure, mD.</p>
<p>The permeability calculation method is as follows: the change in reservoir pressure over time during CBM well production can be determined using Eqs <xref ref-type="disp-formula" rid="e1">1</xref>, <xref ref-type="disp-formula" rid="e3">3</xref>, <xref ref-type="disp-formula" rid="e8">8</xref>. Equation <xref ref-type="disp-formula" rid="e9">9</xref> is employed to calculate permeability when the reservoir pressure exceeds the critical desorption pressure, while Eq. <xref ref-type="disp-formula" rid="e11">11</xref> is utilized for calculating permeability when the reservoir pressure falls below the critical pressure.</p>
</sec>
<sec sec-type="results|discussion" id="s4">
<title>4 Results and discussion</title>
<sec id="s4-1">
<title>4.1 Classification of CBM wells based on SOM neural network</title>
<p>The fine classification of CBM wells necessitates a comprehensive consideration of both the static reservoir parameters and the dynamic drainage parameters. <xref ref-type="fig" rid="F3">Figure 3</xref> is the scatter diagram of the correlation between static parameters and gas production. The permeability, ratio of critical desorption pressure to reservoir pressure (RCR), and gas content are ultimately selected as the key static parameters for well type classification due to their relatively high correlation with gas production. Considering that the key dynamic parameters of CBM wells in existing studies usually include gas production, pressure and permeability (<xref ref-type="bibr" rid="B41">Yuanyuan et al., 2012</xref>; <xref ref-type="bibr" rid="B22">Lou et al., 2013</xref>; <xref ref-type="bibr" rid="B6">Chen et al., 2015</xref>; <xref ref-type="bibr" rid="B8">Fu et al., 2017</xref>; <xref ref-type="bibr" rid="B1">Bao et al., 2020</xref>), this paper selects average daily gas production, average reservoir pressure drop, and change amplitude of permeability as dynamic parameters. The average reservoir pressure drop and dynamic permeability are derived from the mathematical model established in the preceding section. After the parameters of CBM well classification are determined, the intelligent algorithm is used to classify the CBM wells in detail. To minimize the influence of subjective factors in the classification process, this study employs a SOM neural network for automatic clustering. The key feature of this algorithm lies in its ability to determine the optimal number of clusters based on data characteristics, thereby revealing the underlying essential features concealed within the data (<xref ref-type="bibr" rid="B15">Kalteh et al., 2008</xref>).</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Correlation between reservoir static parameters and gas production.</p>
</caption>
<graphic xlink:href="feart-12-1401455-g003.tif"/>
</fig>
<p>The SOM neural network is employed for clustering 71 CBM wells. <xref ref-type="table" rid="T1">Table 1</xref> presents the data required for SOM clustering of some CBM wells. After repeated training, it has been determined that a competition layer consisting of a 10&#xd7;10 network structure (i.e., 100 neurons in total) yields favorable training effects and clustering outcomes, while also exhibiting a rapid learning rate. Consequently, the competition layer (output layer) of the SOM neural network is configured with 100 neurons, whereas the input layer comprises six neurons.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Partial parameter data of CBM wells.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Gas production/(m<sup>3</sup>/d)</th>
<th align="center" style="color:#222A35">Percentage of pressure drop/%</th>
<th align="center" style="color:#222A35">Percentage change of permeability/%</th>
<th align="center" style="color:#222A35">Permeability/mD</th>
<th align="center">RCR</th>
<th align="center" style="color:#222A35">CBM content/(m<sup>3</sup>/t)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center" style="color:#222A35">17,515.18</td>
<td align="center" style="color:#222A35">70.00</td>
<td align="center" style="color:#222A35">2.10</td>
<td align="center" style="color:#222A35">1.60</td>
<td align="center" style="color:#222A35">0.69</td>
<td align="center" style="color:#222A35">23.50</td>
</tr>
<tr>
<td align="center" style="color:#222A35">11,185.40</td>
<td align="center" style="color:#222A35">72.00</td>
<td align="center" style="color:#222A35">5.90</td>
<td align="center" style="color:#222A35">1.10</td>
<td align="center" style="color:#222A35">0.77</td>
<td align="center" style="color:#222A35">21.30</td>
</tr>
<tr>
<td align="center" style="color:#222A35">9,623.57</td>
<td align="center" style="color:#222A35">89.00</td>
<td align="center" style="color:#222A35">3.50</td>
<td align="center" style="color:#222A35">1.50</td>
<td align="center" style="color:#222A35">0.67</td>
<td align="center" style="color:#222A35">22.40</td>
</tr>
<tr>
<td align="center" style="color:#222A35">7,563.00</td>
<td align="center" style="color:#222A35">83.00</td>
<td align="center" style="color:#222A35">3.50</td>
<td align="center" style="color:#222A35">1.30</td>
<td align="center" style="color:#222A35">0.75</td>
<td align="center" style="color:#222A35">24.50</td>
</tr>
<tr>
<td align="center" style="color:#222A35">5,813.68</td>
<td align="center" style="color:#222A35">33.30</td>
<td align="center" style="color:#222A35">20.90</td>
<td align="center" style="color:#222A35">0.95</td>
<td align="center" style="color:#222A35">0.62</td>
<td align="center" style="color:#222A35">21.2</td>
</tr>
<tr>
<td align="center" style="color:#222A35">3,840.53</td>
<td align="center" style="color:#222A35">27.60</td>
<td align="center" style="color:#222A35">30.20</td>
<td align="center" style="color:#222A35">0.88</td>
<td align="center" style="color:#222A35">0.70</td>
<td align="center" style="color:#222A35">20.20</td>
</tr>
<tr>
<td align="center" style="color:#222A35">1,897.08</td>
<td align="center" style="color:#222A35">27.60</td>
<td align="center" style="color:#222A35">19.00</td>
<td align="center" style="color:#222A35">0.87</td>
<td align="center" style="color:#222A35">0.69</td>
<td align="center" style="color:#222A35">17.10</td>
</tr>
<tr>
<td align="center" style="color:#222A35">807.51</td>
<td align="center" style="color:#222A35">7.00</td>
<td align="center" style="color:#222A35">3.90</td>
<td align="center" style="color:#222A35">0.59</td>
<td align="center" style="color:#222A35">0.54</td>
<td align="center" style="color:#222A35">16.30</td>
</tr>
<tr>
<td align="center" style="color:#222A35">503.35</td>
<td align="center" style="color:#222A35">11.40</td>
<td align="center" style="color:#222A35">6.00</td>
<td align="center" style="color:#222A35">0.27</td>
<td align="center" style="color:#222A35">0.63</td>
<td align="center" style="color:#222A35">18.80</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>
<xref ref-type="fig" rid="F4">Figure 4</xref> illustrates the weight connections among network neurons post-training, with each small hexagon representing a neuron. The color intensity of the connection between two hexagons indicates their proximity, reflecting the similarity in features of two input vectors. The hexagon with the deepest shade of black is assigned a weight of 0, while the boundaries of each well type region are defined by neighboring neurons exhibiting progressively darker colors. Within the depicted figure, these dark neurons effectively partition the entire network map into three distinct segments, thereby illustrating how the SOM neural network classifies samples into three distinct types.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Distance distribution of neurons after training.</p>
</caption>
<graphic xlink:href="feart-12-1401455-g004.tif"/>
</fig>
<p>The statistical analysis is conducted on the parameter distribution of three types of wells. <xref ref-type="table" rid="T2">Table 2</xref> presents the classification of well types based on SOM clustering results. Based on the characteristics of various well parameters, CBM wells are categorized into three groups: exceptional wells, potential wells, and inefficient wells. According to the distribution characteristics of dynamic and static parameters of the three wells, it can be observed that the exceptional well exhibits favorable geological conditions and a rational drainage system, thereby fully realizing its gas production potential. Conversely, the potential well possesses advantageous geological conditions but lacks an optimal drainage system, resulting in an underdeveloped gas production potential. Implementing certain measures could potentially transform them into exceptional wells. Lastly, the inefficient well is characterized by unfavorable geological conditions and limited gas production capacity. The exceptional wells are predominantly distributed in close proximity to the central syncline structure, while the potential wells exhibit a wide distribution pattern. The majority of inefficient wells are concentrated in the northeastern and southern regions of the block.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>The classification of well types based on SOM clustering results.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th rowspan="2" align="center">Well type</th>
<th colspan="3" align="center" style="color:#222A35">Dynamic parameter</th>
<th colspan="3" align="center" style="color:#222A35">Static parameter</th>
</tr>
<tr>
<th align="center" style="color:#222A35">Gas/production (m<sup>3</sup>/d)</th>
<th align="center" style="color:#222A35">Percentage of pressure drop/%</th>
<th align="center" style="color:#222A35">Percentage change of permeability/%</th>
<th align="center" style="color:#222A35">Permeability/mD</th>
<th align="center">RCR</th>
<th align="center" style="color:#222A35">CBM content/(m<sup>3</sup>/t)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center" style="color:#222A35">Exceptional well</td>
<td align="center" style="color:#222A35">&#x3e;7,000</td>
<td align="center" style="color:#222A35">60.1&#x2013;93.9</td>
<td align="center" style="color:#222A35">3.2&#x2013;10.5</td>
<td align="center" style="color:#222A35">&#x3e;0.9</td>
<td align="center" style="color:#222A35">&#x3e;0.7</td>
<td align="center" style="color:#222A35">&#x3e;20</td>
</tr>
<tr>
<td align="center" style="color:#222A35">Potential well</td>
<td align="center" style="color:#222A35">1,500&#x2013;7,000</td>
<td align="center" style="color:#222A35">22.3&#x2013;35.2</td>
<td align="center" style="color:#222A35">19.3&#x2013;30.2</td>
<td align="center" style="color:#222A35">&#x3e;0.8</td>
<td align="center" style="color:#222A35">&#x3e;0.5</td>
<td align="center" style="color:#222A35">&#x3e;20</td>
</tr>
<tr>
<td align="center" style="color:#222A35">Inefficient well</td>
<td align="center" style="color:#222A35">25&#x2013;1,500</td>
<td align="center" style="color:#222A35">4.9&#x2013;15.5</td>
<td align="center" style="color:#222A35">2.8&#x2013;8.3</td>
<td align="center" style="color:#222A35">0.5&#x2013;0.9</td>
<td align="center" style="color:#222A35">0.5&#x2013;0.7</td>
<td align="center" style="color:#222A35">14&#x2013;24</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s4-2">
<title>4.2 Production characteristics of CBM wells</title>
<p>The production curves of three exceptional wells are depicted in <xref ref-type="fig" rid="F5">Figure 5</xref>. In order to enhance the clarity of the production curve morphology for different CBM wells, we have simplified the production curves for each well. The production curves of the three wells differed before they reached steady production. Well Z1 and Z3 were already producing gas in the early stages of production, and well Z1 reached its peak gas production quickly, while well Z3 needed a longer time to reach the peak gas production. Well Z2 only produced water in the early production stage and needed a longer time to reach its peak gas production rate. The difference of the curve shape is related to the geological structure and drainage system of the well. Well Z1 and Z3 are located in the syncline, with good gas preservation conditions and high CBM content. Under the reasonable production system, there is enough natural gas supply in the area controlled by a single well, with an average daily gas production of more than 15,000 m<sup>3</sup>/d. The Z2 well is situated along the anticline axis, exhibiting well-developed interbedded layers and fractures with favorable permeability. However, the presence of certain vertical fractures facilitates gas leakage, resulting in limited gas content and prolonged drainage water and pressure reduction period. Despite a reasonable production regime, well Z2 still exhibits lower daily gas production compared to wells Z1 and Z3. The three typical wells exhibit distinct curve morphologies; however, they all demonstrate consistent quantitative characteristics in terms of dynamic parameters. During the 700-day production period, the average daily gas production of the three wells exceeded 10,000 m<sup>3</sup>/d; the average decline in reservoir pressure surpassed 60%; the water production curve exhibited a single peak, and water production was exceptionally low during the intermediate stage of production, with an average daily water production rate below 1.0 m<sup>3</sup>/d. As water production decreases, the influence of effective stress on reservoir permeability damage gradually diminishes, leading to a balance between damage and recovery processes, resulting in stable or slightly reduced permeability (<xref ref-type="bibr" rid="B9">Fu et al., 2002</xref>).</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>The production curves of three exceptional wells.</p>
</caption>
<graphic xlink:href="feart-12-1401455-g005.tif"/>
</fig>
<p>The production curves of the three potential wells are depicted in <xref ref-type="fig" rid="F6">Figure 6</xref>. While the production curves characteristics of Q1, Q2, and Q3 exhibit similarities with Z1, Z2, and Z3, respectively, notable distinctions can be observed in terms of their production dynamic parameters. The average daily gas production of the three potential wells falls significantly short of that of the exceptional well, with Q2 exhibiting the highest average daily gas production at only 6,285 m<sup>3</sup>/d. The water production of the three wells exhibits a substantial increase during the early and middle periods, followed by a gradual decrease in the later stages; however, the decline rate is minimal, ensuring a consistent level of water production. The average decline in reservoir pressure is lower compared to that of exceptional wells. Due to the high water production and low gas production characteristics, the impact of effective stress on reservoir permeability surpasses that of matrix shrinkage, resulting in a significant decrease in reservoir permeability. Some potential wells exhibit comparable geological conditions to the exceptional wells, yet they yield significantly disparate gas production owing to unreasonable production systems, which will be comprehensively analyzed in the subsequent section.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>The production curves of three potential wells.</p>
</caption>
<graphic xlink:href="feart-12-1401455-g006.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="F7">Figure 7</xref> shows the typical production curve of inefficient wells. The average daily gas production of most inefficient wells is less than 1,000 m<sup>3</sup>/d, but the water production is relatively high, with a multi-peak water production curve, and the average daily water production is between 1.15 and 16 m<sup>3</sup>/d. Inefficient wells have a slow decline in reservoir pressure during production, resulting in low gas production. Although the well is characterized by high water production and low gas production, the overall liquid output is comparatively smaller while experiencing a lesser decline in reservoir permeability.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>The production curves of three inefficient wells.</p>
</caption>
<graphic xlink:href="feart-12-1401455-g007.tif"/>
</fig>
</sec>
<sec id="s4-3">
<title>4.3 Productivity differences among different well types</title>
<sec id="s4-3-1">
<title>4.3.1 Analysis of causes for inefficient wells</title>
<p>The production capacity of gas wells is closely correlated with their structural position, leading to varying productivity among gas wells located in different structural positions. Due to the effect of north-south compressive stress, numerous fractures are formed along the axis of the anticline within the study block. The upper overburden exhibits thinness and vertical fractures have developed, facilitating gas migration. Consequently, the gas content within the anticline axis is limited, resulting in low production. On the contrary, the gas preservation conditions in the syncline axis and the flank are good, the CBM is enriched, and the gas content is large. The average and maximum gas production of gas wells are generally higher than those of wells in the axis of the anticline. The inefficient wells are predominantly distributed within the Zhengcun anticline, distant from the Panhe syncline, exhibiting low reservoir permeability, gas content, and RCR, resulting in inadequate gas production capacity.</p>
<p>In the process of CBM production, the generation of coal powder is inevitable. The optimal amount of coal powder facilitates the improvement of reservoir permeability and enhances production capacity of gas wells. However, excessive coal powder can cause reservoir damage and hinder production capacity of gas wells (<xref ref-type="bibr" rid="B18">Kravchenko and Lyerly, 2018</xref>). Some low-efficiency wells in the No. 15 coal seam are seriously affected by coal powder, and the well D5 has a high amount of coal powder blocking the reservoir, resulting in a very low gas production in the early stage. In order to solve this problem, measures of washing coal powder were carried out, and the average daily gas production before washing coal powder was 103.3 m<sup>3</sup>/d, and the average daily gas production after washing coal powder was 446.5 m<sup>3</sup>/d (see <xref ref-type="fig" rid="F8">Figure 8</xref>).</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>The influence of coal power on the production curve of well D5.</p>
</caption>
<graphic xlink:href="feart-12-1401455-g008.tif"/>
</fig>
<p>The continuity of production directly affects the uninterrupted of gas-water flow and the stability of depressurization and desorption in coal reservoirs. The discontinuous production, on the one hand, leads to the rapid rise of reservoir pressure near the wellbore, which stops the seepage action and causes coal powder accumulation to block the seepage channel, reducing the permeability of coal seam; on the other hand, the stagnant flow of free CBM leads to the gradual accumulation of small bubbles in the micro-cracks of coal reservoir into large bubbles, which block the seepage channel and produce Jamin&#x2019;s effect, increasing the gas flow resistance. The blockage of the seepage channels of coal seam caused by multiple switching wells will seriously affect the production in the later stage, and even the gas wells completely failed (<xref ref-type="bibr" rid="B28">Palmer and Mansoori, 1998b</xref>). The hazards of discontinuous production are illustrated by the example of well D6, whose production process is shown in <xref ref-type="fig" rid="F9">Figure 9</xref>. The well underwent a normal drainage process during the initial stage of production; however, due to power and mechanical failures, it experienced 28 shutdowns throughout the production period. These intermittent shutdowns significantly impacted production continuity, leading to a substantial rise in the dynamic liquid level of the gas well and greatly affecting its gas production capacity, resulting in an average daily gas production of only 80.8 m<sup>3</sup>/d. <xref ref-type="fig" rid="F10">Figure 10</xref> illustrates the correlation between average daily gas production and shutdown frequency of inefficient wells. Statistical analysis indicates that a lower number of shutdowns is associated with higher gas production and greater potential for gas production.</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>The influence of discontinuous production system on the production curve of well D6.</p>
</caption>
<graphic xlink:href="feart-12-1401455-g009.tif"/>
</fig>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>The correlation between average daily gas production and shutdown frequency of inefficient wells.</p>
</caption>
<graphic xlink:href="feart-12-1401455-g010.tif"/>
</fig>
</sec>
<sec id="s4-3-2">
<title>4.3.2 Analysis of the productivity differences between exceptional wells and potential wells</title>
<p>Both exceptional wells and potential wells have good geological conditions. The main reason for the difference in productivity between the two is the different matching degree between drainage system and reservoir conditions. The CBM reservoir in the study area is an undersaturated reservoir, necessitating decompression and pressure reduction during the development process. The initial rate of pressure reduction significantly influences both overall pressure reduction and the expansion of the pressure reduction funnel.</p>
<p>
<xref ref-type="fig" rid="F11">Figure 11</xref> shows the production curves of the exceptional well Z4 and the potential well Q4. The two wells are adjacent to each other and have the same geological conditions, but the productivity differences are very obvious. The average daily gas production of well Z4 is 14061 m<sup>3</sup>/d, accompanied by an average daily water production of 1.65 m<sup>3</sup>/d. In contrast, well Q4 exhibits a significantly lower average daily gas production of only 1,606 m<sup>3</sup>/d, coupled with an average daily water production of 0.91 m<sup>3</sup>/d. This disparity can be primarily attributed to the differential pressure drop rates experienced during the early stage of drainage. The pressure drop rate in Well Z4 before peak production was 0.01 MPa/d, with a 76% reduction in bottom flow pressure, while the pressure drop rate in Well Q4 before peak production was 0.076 MPa/d, with a 70% reduction in bottom flow pressure. When the disparity in bottom flow pressure drop is insignificant, the pressure drop rate of Q4 well exhibits an approximate 9 times faster than that of Z4 well. This phenomenon potentially leads to an elevation in the pressure differential between overlying rock and pore spaces within the rock, resulting in coal matrix compression, closure of minuscule fissures within the rock, and gradual reduction in permeability&#x2014;thus inducing a pressure-sensitive effect. On the other hand, this phenomenon will result in an excessive growth rate of gas production, thereby accelerating the fluid flow velocity within the coal reservoir fissures. The rapid fluid flow carries a substantial amount of coal powder through the seepage channel, leading to blockages in the coal seam fissures due to the velocity-sensitive effect (<xref ref-type="bibr" rid="B35">Shi and Durucan, 2005b</xref>). The internal permeability of the reservoir is compromised under the combined influence of pressure-sensitive and velocity-sensitive effects, thereby impeding the outward expansion rate of the pressure drop funnel. Consequently, this limited the extent of desorption zone, ultimately leading to a decline in gas production from the potential well and an incomplete realization of its gas production potential.</p>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>The production curves of the exceptional well Z4 and the potential well Q4.</p>
</caption>
<graphic xlink:href="feart-12-1401455-g011.tif"/>
</fig>
</sec>
<sec id="s4-3-3">
<title>4.3.3 Optimization method of pressure drop rate</title>
<p>The rate of drainage pressure drop in the early stage of CBM well production should not be excessively high, as this may lead to a reduction in desorption area and a significant decline in gas production potential. The production characteristics observed in Q4 well align with this particular pressure drop trend. The system of drainage with a large pressure drop velocity first and a small one later is beneficial to improve the gas production effect. The production characteristics of Well Z4 align with the observed pressure drop path, thereby enhancing its production performance. Subsequently, the underlying internal mechanism responsible for this phenomenon is discussed.</p>
<p>The seepage states under three different pressure drops are illustrated in <xref ref-type="fig" rid="F12">Figure 12</xref>. The black solid line a represents the stage when the BHP reaches the critical desorption pressure. The yellow dotted line represents the pressure distribution after rapidly reducing the BHP from the state a, and the purple dotted line represents the pressure distribution after slowly reducing the BHP from the state a.</p>
<fig id="F12" position="float">
<label>FIGURE 12</label>
<caption>
<p>The seepage states under three different pressure drops (<xref ref-type="bibr" rid="B39">Yan et al., 2021</xref>).</p>
</caption>
<graphic xlink:href="feart-12-1401455-g012.tif"/>
</fig>
<p>The difference between the BHP and the critical desorption pressure is constrained for a, resulting in a limited amount of gas desorbed into the micro-fracture system. The majority of gas is predominantly stored on the surface of the micro-fracture as minute bubbles, which lack flowability but contribute to a reduction in effective flow area and weaken the flow ability of the fluid. However, the flow resistance in this case is consistently negligible, and the primary consequence of reducing BHP is an augmentation in gas production rather than an escalation in resistance near the bottom hole. Consequently, the reduction of BHP in this scenario exerts a positive influence on the expansion of desorption zone.</p>
<p>In scenario b, the decline in BHP continues at an accelerated rate, facilitating a sustained period of rapid production growth and reaches the peak in a short period of time (refer to <xref ref-type="fig" rid="F13">Figure 13</xref>). However, it is important to note that the rapid decrease in BHP results in a significant release of desorption gas around the bottom hole, leading to gas-water two-phase seepage at the bottom hole, which will sharply increase the flow resistance in the near-well area, resulting in the difficulty of fluid flowing to the bottom hole from the distal coal seam (<xref ref-type="bibr" rid="B36">Wang et al., 2023</xref>). Therefore, the pressure drop is difficult to be transmitted to the distal end, resulting in the slow expansion of the desorption area, which is consistent with the production characteristics of Q4 well. The scenario c can be interpreted as a gradual and controlled decline in BHP following the relatively large BHP decline for a certain duration, namely scenario a. On one hand, this pressure drop mode hinders the increase of gas saturation in the near-wellbore region, mitigates the occurrence of gas lock phenomenon in the near-wellbore area, and facilitates smooth fluid flow to reach the bottom of the well, thereby promoting desorption area expansion. On the other hand, the BHP is still slightly decreased, the production pressure difference is constantly increased, and the average daily gas production will also increase, which will reduce the gas saturation of the coal seam to a certain extent, thus reducing the flow resistance of the coal seam (refer to <xref ref-type="fig" rid="F13">Figure 13</xref>).</p>
<fig id="F13" position="float">
<label>FIGURE 13</label>
<caption>
<p>The seepage states under three different pressure.</p>
</caption>
<graphic xlink:href="feart-12-1401455-g013.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="F14">Figures 14</xref>, <xref ref-type="fig" rid="F15">15</xref> represent the reasonable value of the average pressure drop rate in the early stage of drainage in the study area analyzed from a statistical perspective. The correlation analysis reveals a logarithmic relationship between the pressure drop rate and the peak gas production volume during the initial production period. As the pressure drop rate gradually increases, there is a decrease in peak gas and water production. In the initial stage of gas well production, when the pressure drop rate exceeds 0.005 MPa/d, most wells exhibit a peak gas production of less than 10,000 m<sup>3</sup>/d; whereas when the pressure drop rate falls below 0.005 MPa/d, the peak gas production is generally greater than 13,000 m<sup>3</sup>/d. The well Z12, which exhibits the highest gas production in the P block, demonstrates an exceptionally minimal pressure drop during production, measuring only 0.00067 MPa/d. Moreover, it achieves a peak gas production rate of 52,316 m<sup>3</sup>/d and maintains a steady gas production rate of 43,000 m<sup>3</sup>/d. The initial decline rate of drainage pressure in most high production water wells is typically below 0.02 MPa/d, whereas wells with a decline rate exceeding 0.02 MPa/d exhibit limited peak water production, usually less than 10 m<sup>3</sup>/d. Therefore, to ensure stable and continuous gas production, the pressure drop rate should not be too fast in the early stage of production, and the average pressure drop rate should be controlled within 0.005&#x2013;0.02 MPa/d.</p>
<fig id="F14" position="float">
<label>FIGURE 14</label>
<caption>
<p>Correlation between pressure drop rate and peak gas production in early production stage.</p>
</caption>
<graphic xlink:href="feart-12-1401455-g014.tif"/>
</fig>
<fig id="F15" position="float">
<label>FIGURE 15</label>
<caption>
<p>Correlation between pressure drop rate and peak water production in early production stage.</p>
</caption>
<graphic xlink:href="feart-12-1401455-g015.tif"/>
</fig>
</sec>
</sec>
</sec>
<sec sec-type="conclusion" id="s5">
<title>5 Conclusion</title>
<p>
<list list-type="simple">
<list-item>
<p>(1) Utilizing static and dynamic parameters, the SOM neural network is employed for the classification of CBM wells. The CBM wells in P Block&#x2019;s No. 15 coal seam are categorized into three types: exceptional wells, potential wells, and inefficient wells. The exceptional well exhibits favorable geological conditions and a rational drainage system, thereby fully realizing its gas production potential. Conversely, the potential well possesses advantageous geological conditions but lacks an optimal drainage system, resulting in an underdeveloped gas production potential. Implementing certain measures could potentially transform them into exceptional wells. Lastly, the inefficient well is characterized by unfavorable geological conditions and limited gas production capacity.</p>
</list-item>
<list-item>
<p>(2) Poor reservoir conditions, high content of coal powder and discontinuous drainage system are the main causes of inefficient wells. Both exceptional wells and potential wells have good geological conditions. The main reason for the difference in productivity between the two is the different matching degree between drainage system and reservoir conditions. The pressure decline rate of most potential wells is excessively rapid during the early stage of drainage, leading to pressure-sensitive, velocity-sensitive and jiamin effect effects on coal reservoirs, resulting in damage to the reservoir and a gradual decrease in permeability that ultimately impacts gas production during later stages.</p>
</list-item>
<list-item>
<p>(3) A logarithmic relationship between the pressure drop rate and the peak gas production volume during the initial production period. In the initial stage of gas well production, when the pressure drop rate exceeds 0.005 MPa/d, most wells exhibit a peak gas production of less than 10,000 m<sup>3</sup>/d; the initial decline rate of drainage pressure in most high production water wells is typically below 0.02 MPa/d, whereas wells with a decline rate exceeding 0.02 MPa/d exhibit limited peak water production. Therefore, it is suggested that the average pressure drop rate in the early stage of drainage should be controlled within 0.005&#x2013;0.02 MPa/d.</p>
</list-item>
</list>
</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s6">
<title>Data availability statement</title>
<p>The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.</p>
</sec>
<sec id="s7">
<title>Author contributions</title>
<p>XC: Writing&#x2013;review and editing, Writing&#x2013;original draft, Visualization, Validation, Supervision, Resources, Project administration, Methodology, Investigation, Funding acquisition, Formal analysis, Data curation, Conceptualization. YG: Formal Analysis, Methodology, Supervision, Writing&#x2013;review and editing, Investigation. YM: Methodology, Resources, Writing&#x2013;original draft, Writing&#x2013;review and editing.</p>
</sec>
<sec sec-type="funding-information" id="s8">
<title>Funding</title>
<p>The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was supported by the China Postdoctoral Science Foundation (2023M733916).</p>
</sec>
<ack>
<p>We would like to thank China Postdoctoral Science Foundation (2023M733916) for their support of this work. The authors also thank the reviewers for their valuable comments.</p>
</ack>
<sec sec-type="COI-statement" id="s9">
<title>Conflict of interest</title>
<p>Authors XC, YG, and YW were employed by China National Offshore Oil Corporation Research Institute Co. Ltd.</p>
</sec>
<sec sec-type="disclaimer" id="s10">
<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">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Ju</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Yin</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Xiong</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Qi</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Influence of reservoir properties on the methane adsorption capacity and fractal features of coal and shale in the upper Permian coal measures of the South Sichuan coalfield, China</article-title>. <source>Energy Explor. Exploitation</source> <volume>38</volume> (<issue>1</issue>), <fpage>57</fpage>&#x2013;<lpage>78</lpage>. <pub-id pub-id-type="doi">10.1177/0144598719877527</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bustin</surname>
<given-names>R. M.</given-names>
</name>
</person-group> (<year>1997</year>). <article-title>Importance of fabric and composition on the stress sensitivity of permeability in some coals, northern Sydney basin, Australia: relevance to coalbed methane exploitation</article-title>. <source>AAPG Bull.</source> <volume>81</volume> (<issue>11</issue>), <fpage>1894</fpage>&#x2013;<lpage>1908</lpage>. <pub-id pub-id-type="doi">10.1306/3b05c65e-172a-11d7-8645000102c1865d</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cai</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Pan</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Connell</surname>
<given-names>L. D.</given-names>
</name>
<etal/>
</person-group> (<year>2014</year>). <article-title>Effects of pressure and temperature on gas diffusion and flow for primary and enhanced coalbed methane recovery</article-title>. <source>Energy Explor. Exploitation</source> <volume>32</volume> (<issue>4</issue>), <fpage>601</fpage>&#x2013;<lpage>619</lpage>. <pub-id pub-id-type="doi">10.1260/0144-5987.32.4.601</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Tao</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>W.</given-names>
</name>
<etal/>
</person-group> (<year>2016</year>). <article-title>Discussion about &#x201c;critical depth&#x201d; of deep coalbed methane in Zhengzhuang area, Qinshui Basin</article-title>. <source>J. China Coal Soc.</source> <volume>41</volume> (<issue>12</issue>), <fpage>3069</fpage>&#x2013;<lpage>3075</lpage>.</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Tao</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>J.</given-names>
</name>
<etal/>
</person-group> (<year>2017</year>). <article-title>
<italic>In-situ</italic> stress measurements and stress distribution characteristics of coal reservoirs in major coalfields in China: implication for coalbed methane (CBM) development</article-title>. <source>Int. J. Coal Geol.</source> <volume>182</volume>, <fpage>66</fpage>&#x2013;<lpage>84</lpage>. <pub-id pub-id-type="doi">10.1016/j.coal.2017.09.009</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>Y. X.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>D. M.</given-names>
</name>
<name>
<surname>Yao</surname>
<given-names>Y. B.</given-names>
</name>
<name>
<surname>Cai</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Dynamic permeability change during coalbed methane production and its controlling factors</article-title>. <source>J. Nat. Gas Sci. Eng.</source> <volume>25</volume>, <fpage>335</fpage>&#x2013;<lpage>346</lpage>. <pub-id pub-id-type="doi">10.1016/j.jngse.2015.05.018</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Dickey</surname>
<given-names>M. W.</given-names>
</name>
</person-group> (<year>2006</year>). &#x201c;<article-title>Economic pumping technology for coalbed methane (CBM), stripper oil, and shallow gas well deliquification</article-title>,&#x201d; in <conf-name>SPE Eastern Regional Meeting</conf-name>, <conf-loc>Canton, Ohio, USA</conf-loc>, <conf-date>October, 2006</conf-date>.</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Tao</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>S.</given-names>
</name>
<etal/>
</person-group> (<year>2017</year>). <article-title>Characteristics of pore structure and fractal dimension of low-rank coal: a case study of Lower Jurassic Xishanyao coal in the southern Junggar Basin, NW China</article-title>. <source>Fuel</source> <volume>193</volume>, <fpage>254</fpage>&#x2013;<lpage>264</lpage>. <pub-id pub-id-type="doi">10.1016/j.fuel.2016.11.069</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Qin</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>G.</given-names>
</name>
<etal/>
</person-group>(<year>2002</year>). <article-title>Experimental research of influence of coal matrix shrinkage on permeability</article-title>. <source>Journal of China University of Mining</source> &#x26; <source>Technology</source>, <volume>31(2)</volume>: <fpage>129</fpage>&#x2013;<lpage>131</lpage>. <pub-id pub-id-type="doi">10.3321/j.issn:1000-1964.2002.02.005</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gao</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Kong</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Peng</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2023b</year>). <article-title>Pore-scale mechanisms and hysteresis effect during multi-cycle injection and production process in underground hydrogen storage reservoir</article-title>. <source>Energy</source> <volume>283</volume>, <fpage>129007</fpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2023.129007</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gao</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Cheng</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Luo</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2023a</year>). <article-title>Structural evolution and characterization of organic-rich shale from macroscopic to microscopic resolution: the significance of tectonic activity</article-title>. <source>Adv. Geo-Energy Res.</source> <volume>10</volume> (<issue>2</issue>), <fpage>84</fpage>&#x2013;<lpage>90</lpage>. <pub-id pub-id-type="doi">10.46690/ager.2023.11.03</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hou</surname>
<given-names>C. H.</given-names>
</name>
<name>
<surname>Shi</surname>
<given-names>J. T.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Peng</surname>
<given-names>L. S.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Chang</surname>
<given-names>Y. C.</given-names>
</name>
<etal/>
</person-group> (<year>2018</year>). <article-title>Classification and evaluation method of coalbed methane wells in M block</article-title>. <source>Nat. Gas. Ind.</source> <volume>38</volume> (<issue>S1</issue>), <fpage>74</fpage>&#x2013;<lpage>79</lpage>.</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ifrene</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Irofti</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Ni</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Egenhoff</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Pothana</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>New insights into fracture porosity estimations using machine learning and advanced logging tools</article-title>. <source>Fuels</source> <volume>4</volume> (<issue>3</issue>), <fpage>333</fpage>&#x2013;<lpage>353</lpage>. <pub-id pub-id-type="doi">10.3390/fuels4030021</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jiang</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Duan</surname>
<given-names>Y.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <article-title>Reservoir characteristics and gas production potential of deep coalbed methane: insights from the no. 15 coal seam in shouyang block, Qinshui Basin, China</article-title>. <source>Unconv. Resour.</source> <volume>2</volume>, <fpage>12</fpage>&#x2013;<lpage>20</lpage>. <pub-id pub-id-type="doi">10.1016/j.uncres.2022.06.001</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kalteh</surname>
<given-names>A. M.</given-names>
</name>
<name>
<surname>Hjorth</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Berndtsson</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Review of the self-organizing map (SOM) approach in water resources: analysis, modelling and application</article-title>. <source>Environ. Model. Softw.</source> <volume>23</volume> (<issue>7</issue>), <fpage>835</fpage>&#x2013;<lpage>845</lpage>. <pub-id pub-id-type="doi">10.1016/j.envsoft.2007.10.001</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kong</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Gao</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Lian</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Characteristics of gas-oil contact and mobilization limit during gas-assisted gravity drainage process</article-title>. <source>Adv. Geo-Energy Res.</source> <volume>6</volume> (<issue>2</issue>), <fpage>169</fpage>&#x2013;<lpage>176</lpage>. <pub-id pub-id-type="doi">10.46690/ager.2022.02.08</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kong</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Gao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Sarma</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Experimental investigation of immiscible water-alternating-gas injection in ultra-high water-cut stage reservoir</article-title>. <source>Adv. Geo-Energy Res.</source> <volume>5</volume> (<issue>2</issue>), <fpage>139</fpage>&#x2013;<lpage>152</lpage>. <pub-id pub-id-type="doi">10.46690/ager.2021.02.04</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kravchenko</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Lyerly</surname>
<given-names>H. K.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>The impact of coal-powered electrical plants and coal ash impoundments on the health of residential communities</article-title>. <source>N. C. Med. J.</source> <volume>79</volume> (<issue>5</issue>), <fpage>289</fpage>&#x2013;<lpage>300</lpage>. <pub-id pub-id-type="doi">10.18043/ncm.79.5.289</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lai</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Fu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>Z. h.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>A drainage data-based calculation method for coalbed permeability</article-title>. <source>J. Geophys. Eng.</source> <volume>10</volume> (<issue>6</issue>), <fpage>065005</fpage>. <pub-id pub-id-type="doi">10.1088/1742-2132/10/6/065005</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Ranjith</surname>
<given-names>P. G.</given-names>
</name>
<name>
<surname>Cai</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Geological controls on variable gas concentrations: a case study of the northern Gujiao Block, northwestern Qinshui Basin, China</article-title>. <source>Mar. Petroleum Geol.</source> <volume>92</volume>, <fpage>582</fpage>&#x2013;<lpage>596</lpage>. <pub-id pub-id-type="doi">10.1016/j.marpetgeo.2017.11.022</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Sang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Xue</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Lan</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Ren</surname>
<given-names>B.</given-names>
</name>
<etal/>
</person-group> (<year>2018</year>). <article-title>Evolution and geochemical characteristics of gas phase fluid and its response to inter-well interference during multi-well drainage of coalbed methane</article-title>. <source>J. Petroleum Sci. Eng.</source> <volume>162</volume>, <fpage>491</fpage>&#x2013;<lpage>501</lpage>. <pub-id pub-id-type="doi">10.1016/j.petrol.2017.10.007</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lou</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Song</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Dong</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Productivity equation of fractured well in CBM reservoirs</article-title>. <source>J. Nat. Gas Sci. Eng.</source> <volume>11</volume>, <fpage>39</fpage>&#x2013;<lpage>45</lpage>. <pub-id pub-id-type="doi">10.1016/j.jngse.2012.12.001</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lv</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Luo</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Production characteristics and the key factors in high-rank coalbed methane fields: a case study on the Fanzhuang Block, Southern Qinshui Basin, China</article-title>. <source>Int. J. Coal Geol.</source> <volume>96</volume>, <fpage>93</fpage>&#x2013;<lpage>108</lpage>. <pub-id pub-id-type="doi">10.1016/j.coal.2012.03.009</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Moore</surname>
<given-names>T. A.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Coalbed methane: a review</article-title>. <source>Int. J. Coal Geol.</source> <volume>101</volume>, <fpage>36</fpage>&#x2013;<lpage>81</lpage>. <pub-id pub-id-type="doi">10.1016/j.coal.2012.05.011</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ni</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Ling</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Afari</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Different ways to approach shale reservoirs&#x2019; CO<sub>2</sub> storage potential in America</article-title>. <source>Heliyon</source> <volume>9</volume>, <fpage>e18458</fpage>. <pub-id pub-id-type="doi">10.1016/j.heliyon.2023.e18458</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nie</surname>
<given-names>B. S.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>X. Q.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>E. Y.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>Mechanism and modes of gas diffusion in coal seams</article-title>. <source>China Saf. Sci. J.</source> <volume>10</volume> (<issue>6</issue>), <fpage>24</fpage>&#x2013;<lpage>28</lpage>. <pub-id pub-id-type="doi">10.3969/j.issn.1003-3033.2000.06.006</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Palmer</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Mansoori</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>1998a</year>). <article-title>How permeability depends on stress and pore pressure in coalbeds: a new model</article-title>. <source>SPE Reserv. Eval. Eng.</source> <volume>1</volume> (<issue>06</issue>), <fpage>539</fpage>&#x2013;<lpage>544</lpage>. <pub-id pub-id-type="doi">10.2118/52607-pa</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Palmer</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Mansoori</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>1998b</year>). <article-title>How permeability depends on stress and pore pressure in coalbeds: a new model</article-title>. <source>SPE Reserv. Eval. Eng.</source> <volume>1</volume> (<issue>06</issue>), <fpage>539</fpage>&#x2013;<lpage>544</lpage>. <pub-id pub-id-type="doi">10.2118/52607-pa</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Peng</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zou</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>G.</given-names>
</name>
<etal/>
</person-group> (<year>2017</year>). <article-title>Factors affecting coalbed methane (CBM) well productivity in the Shizhuangnan block of southern Qinshui basin, North China: investigation by geophysical log, experiment and production data</article-title>. <source>Fuel</source> <volume>191</volume>, <fpage>427</fpage>&#x2013;<lpage>441</lpage>. <pub-id pub-id-type="doi">10.1016/j.fuel.2016.11.071</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Qin</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Song</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Fu</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>The mechanism of the flowing ground water impacting on coalbed gas content</article-title>. <source>Chin. Sci. Bull.</source> <volume>50</volume>, <fpage>118</fpage>&#x2013;<lpage>123</lpage>. <pub-id pub-id-type="doi">10.1007/bf03184093</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Qiu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Song</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Response law and indicator selection of seismic wave velocity for coal seam outburst risk</article-title>. <source>Adv. Geo-Energy Res.</source> <volume>9</volume> (<issue>3</issue>), <fpage>198</fpage>&#x2013;<lpage>210</lpage>. <pub-id pub-id-type="doi">10.46690/ager.2023.09.07</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>M. x.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Geological controls over coal-bed methane well production in southern Qinshui basin</article-title>. <source>Procedia Earth Planet. Sci.</source> <volume>1</volume> (<issue>1</issue>), <fpage>917</fpage>&#x2013;<lpage>922</lpage>. <pub-id pub-id-type="doi">10.1016/j.proeps.2009.09.142</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shi</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Durucan</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2005a</year>). <article-title>A model for changes in coalbed permeability during primary and enhanced methane recovery</article-title>. <source>SPE Reserv. Eval. Eng.</source> <volume>8</volume> (<issue>04</issue>), <fpage>291</fpage>&#x2013;<lpage>299</lpage>. <pub-id pub-id-type="doi">10.2118/87230-pa</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shi</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Fang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Lu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Hou</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>X.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>A new coal reservoir permeability model considering the influence of pulverized coal blockage and its application</article-title>. <source>Nat. Gas. Ind. B</source> <volume>8</volume> (<issue>1</issue>), <fpage>67</fpage>&#x2013;<lpage>78</lpage>. <pub-id pub-id-type="doi">10.1016/j.ngib.2020.06.004</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shi</surname>
<given-names>J. Q.</given-names>
</name>
<name>
<surname>Durucan</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2005b</year>). <article-title>A model for changes in coalbed permeability during primary and enhanced methane recovery</article-title>. <source>SPE Reserv. Eval. Eng.</source> <volume>8</volume> (<issue>04</issue>), <fpage>291</fpage>&#x2013;<lpage>299</lpage>. <pub-id pub-id-type="doi">10.2118/87230-pa</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Dai</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Shen</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Soltanian</surname>
<given-names>M. R.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>An integrated multi-scale model for CO2 transport and storage in shale reservoirs</article-title>. <source>Appl. Energy</source> <volume>331</volume>, <fpage>120444</fpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2022.120444</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2024b</year>). <article-title>Multiscale model for hydrogen transport and storage in shale reservoirs</article-title>. <source>SPE J.</source>, <fpage>1</fpage>&#x2013;<lpage>27</lpage>. <pub-id pub-id-type="doi">10.2118/219472-pa</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>Y. W.</given-names>
</name>
<name>
<surname>Dai</surname>
<given-names>Z. X.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>G. S.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Xia</surname>
<given-names>Y. Z.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>Y. H.</given-names>
</name>
</person-group> (<year>2024a</year>). <article-title>A hybrid physics-informed data-driven neural network for CO2 storage in depleted shale reservoirs</article-title>. <source>Petroleum Sci.</source> <volume>21</volume> (<issue>1</issue>), <fpage>286</fpage>&#x2013;<lpage>301</lpage>. <pub-id pub-id-type="doi">10.1016/j.petsci.2023.08.032</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yan</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Guan</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Q.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>A prediction model for pressure propagation and production boundary during coalbed methane development</article-title>. <source>Energy &#x26; Fuels</source> <volume>35</volume> (<issue>2</issue>), <fpage>1219</fpage>&#x2013;<lpage>1233</lpage>. <pub-id pub-id-type="doi">10.1021/acs.energyfuels.0c03354</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yi</surname>
<given-names>Y. X.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>S. H.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>S. H.</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>X. L.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>K. F.</given-names>
</name>
<name>
<surname>Dang</surname>
<given-names>F.</given-names>
</name>
<etal/>
</person-group> (<year>2019</year>). <article-title>Analysis on the type of reservoir pressure drop and drainage control of coalbed methane well in the southern block of Shizhuang</article-title>. <source>Coal Geol. Explor.</source> <volume>47</volume> (<issue>05</issue>), <fpage>118</fpage>&#x2013;<lpage>126</lpage>.</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yuanyuan</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Xianjie</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Caifeng</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Production characteristics and affecting factors of high&#x2014;mid rank coalbed methane wells: taking Fanzhuang and Hancheng mining areas as examples</article-title>. <source>Petroleum Explor. Dev.</source> <volume>39</volume> (<issue>6</issue>), <fpage>773</fpage>&#x2013;<lpage>777</lpage>. <pub-id pub-id-type="doi">10.1016/s1876-3804(12)60103-7</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Ju</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Bu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>M.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>Multi-scale multi-dimensional characterization of clay-hosted pore networks of shale using FIBSEM, TEM, and X-ray micro-tomography: implications for methane storage and migration</article-title>. <source>Appl. Clay Sci.</source> <volume>213</volume>, <fpage>106239</fpage>. <pub-id pub-id-type="doi">10.1016/j.clay.2021.106239</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>