<?xml version="1.0" encoding="UTF-8"?>
<!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. Energy Res.</journal-id>
<journal-title>Frontiers in Energy Research</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Energy Res.</abbrev-journal-title>
<issn pub-type="epub">2296-598X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">763052</article-id>
<article-id pub-id-type="doi">10.3389/fenrg.2021.763052</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Energy Research</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Experimental Study on the Fracture Behaviors Effect of Stress Conditions on the Fracture of Shale by Super-Critical Carbon Dioxide</article-title>
<alt-title alt-title-type="left-running-head">Liu et&#x20;al.</alt-title>
<alt-title alt-title-type="right-running-head">Carbon Sequestration</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Liu</surname>
<given-names>Guojun</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1431813/overview"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Zhao</surname>
<given-names>Yuan</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1461309/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Cheng</surname>
<given-names>Yugang</given-names>
</name>
<xref ref-type="aff" rid="aff5">
<sup>5</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1">
<label>
<sup>1</sup>
</label>Hunan Engineering Research Center of Structural Safety and Disaster Prevention for Urban Underground Infrastructure, Hunan City University, <addr-line>Yiyang</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<label>
<sup>2</sup>
</label>College of Civil Engineering and Mechanics, Xiangtan University, <addr-line>Xiangtan</addr-line>, <country>China</country>
</aff>
<aff id="aff3">
<label>
<sup>3</sup>
</label>Sinohydro Bureau 8 Co. LTD., POWERCHINA, <addr-line>Changsha</addr-line>, <country>China</country>
</aff>
<aff id="aff4">
<label>
<sup>4</sup>
</label>Department of Hydraulic Engineering, Tsinghua University, <addr-line>Beijing</addr-line>, <country>China</country>
</aff>
<aff id="aff5">
<label>
<sup>5</sup>
</label>State Key Laboratory of Mountain Bridge and Tunnel Engineering, Chongqing Jiaotong University, <addr-line>Chongqing</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/1256438/overview">Guowen Xu</ext-link>, Colorado School of Mines, United&#x20;States</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/1459776/overview">Delei Shang</ext-link>, Shenzhen University, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1036480/overview">Yong Li</ext-link>, Chongqing University, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Guojun Liu, <email>liuguojun@hncu.edu.cn</email>; Yuan Zhao, <email>zhaoyuan@cqu.edu.cn</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Carbon Capture, Utilization and Storage, a section of the journal Frontiers in Energy Research</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>10</day>
<month>12</month>
<year>2021</year>
</pub-date>
<pub-date pub-type="collection">
<year>2021</year>
</pub-date>
<volume>9</volume>
<elocation-id>763052</elocation-id>
<history>
<date date-type="received">
<day>23</day>
<month>08</month>
<year>2021</year>
</date>
<date date-type="accepted">
<day>29</day>
<month>10</month>
<year>2021</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2021 Liu, Zhao and Cheng.</copyright-statement>
<copyright-year>2021</copyright-year>
<copyright-holder>Liu, Zhao and Cheng</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&#x20;terms.</p>
</license>
</permissions>
<abstract>
<p>This paper examines the fracture propagation problems of supercritical carbon fracturing in low permeability shale. Acoustic emission monitoring and computerized tomography (CT) scanning methods were used to study the influence of initial stress ratios on crack initiation and propagation crack in fracturing experiments. The results show that crack initiation pressure and crack morphology are very different under different stress conditions. Under the condition of constant confining pressure, when the initial stress ratio &#x3bb; &#x3d; 1, cracks are mainly in a horizontal direction; while for an initial stress ratio of &#x3bb; &#x3c; 1, cracks are mainly in a vertical direction. With the decrease of &#x3bb;, crack initiation pressure, reopening pressure, and fracturing liquid volume also decrease, and crack propagation is not as obvious. According to CT scanning results, the crack propagation direction is the same as the maximum principal stress, and fewer cracks are initiated with a smaller initial stress ratio. Based on the acoustic emission characteristics, the fracturing process (including crack initiation, propagation, and closure), can be divided into three stages: 1) the pressure accumulation in the wellbore, 2) Pump Closure; and 3) crack reopening. This study provides the basis for a reasonable selection of shale gas fracturing formation and geo-sequestration of greenhouse gas&#x20;CO<sub>2</sub>.</p>
</abstract>
<kwd-group>
<kwd>shale gas</kwd>
<kwd>super-critical carbon dioxide</kwd>
<kwd>carbon sequestration</kwd>
<kwd>mechanical characteristics</kwd>
<kwd>acoustic emission</kwd>
</kwd-group>
<contract-sponsor id="cn001">Natural Science Foundation of Hunan Province<named-content content-type="fundref-id">10.13039/501100004735</named-content>
</contract-sponsor>
<contract-sponsor id="cn002">Education Department of Hunan Province<named-content content-type="fundref-id">10.13039/100009377</named-content>
</contract-sponsor>
</article-meta>
</front>
<body>
<sec id="s1">
<title>Introduction</title>
<p>China is very rich in shale gas, with resources of 3.608 &#xd7; 1013&#xa0;m<sup>3</sup>, which has great potential to develop the shale gas industry (<xref ref-type="bibr" rid="B19">Wang, 2011</xref>). However, as the shale matrix has characteristics of ultra-low permeability, it is necessary to explore the relationship between natural fractures and artificial fracturing in forming a fracture network and a high permeable path to increase production. At present, hydraulic fracturing technology is the main method that greatly improves the yield of shale gas (<xref ref-type="bibr" rid="B10">Jin, et&#x20;al., 2012</xref>; <xref ref-type="bibr" rid="B11">Jin and Shah, 2013</xref>); however, the process consumes a large amount of water during hydraulic fracturing (<xref ref-type="bibr" rid="B18">Vengosh et&#x20;al., 2014</xref>; <xref ref-type="bibr" rid="B1">Brantley et&#x20;al., 2014</xref>) Furthermore, chemical additives in the fracturing fluids may cause serious pollution to the reservoir and surrounding environment as hydraulic fracturing in China could cause swelling clays to obstruct gas channels (<xref ref-type="bibr" rid="B17">Small et&#x20;al., 2014</xref>; <xref ref-type="bibr" rid="B6">Huang et&#x20;al., 2015</xref>), which can easily decrease permeability and gas production (<xref ref-type="bibr" rid="B13">Lai et&#x20;al., 2005</xref>) Additionally, most favorable exploration areas are short of water, meaning non-aqueous fracturing technology is of great significance for the exploration and development of shale gas in China. In recent years, some scholars have put forward supercritical CO<sub>2</sub> (when temperature and pressure are above 304.1&#xa0;K and 7.38 MPa, CO<sub>2</sub> will reach a supercritical state) as a new fracturing fluid for shale gas development (<xref ref-type="bibr" rid="B9">Ishida et&#x20;al., 2012</xref>; <xref ref-type="bibr" rid="B2">Chen L. Q. et&#x20;al., 2015</xref>; <xref ref-type="bibr" rid="B14">Middleton et&#x20;al., 2015</xref>)<sup>.</sup> Supercritical CO<sub>2</sub> is beneficial for improving the permeability of the reservoir and using supercritical CO<sub>2</sub>, as the fracturing fluid reduces reservoir damage and can obtain higher oil and gas recovery (<xref ref-type="bibr" rid="B4">Gupta et&#x20;al., 2005</xref>; <xref ref-type="bibr" rid="B12">Kolle, 2000</xref>). Domestic research also shows that supercritical CO<sub>2</sub> can make clay mineral dehydration shrink the shale gas reservoir as porosity increases, and effectively protects the reservoir. In addition, the adsorption capacity of shale to CO<sub>2</sub> is much larger than CH<sub>4</sub>, which can effectively replace CH<sub>4</sub> in shale and improve the recovery of shale gas, while finally achieving CO<sub>2</sub> geological sequestration (<xref ref-type="bibr" rid="B20">Xie et&#x20;al.,.2014</xref>). It is therefore of great scientific significance to study the mechanical characteristics and failure mechanisms of shale gas in the supercritical CO<sub>2</sub> fracturing process, to find a suitable method for the development of shale gas reservoirs in China.</p>
<p>
<xref ref-type="bibr" rid="B8">Hubbert and Willis (1972)</xref> first proposed crack initiation criteria for longitudinal fractures and suggested that fracture pressure was related to horizontal principal stress, rock tensile strength, and pore pressure. They outline that fracture pressure has nothing to do with the size, elasticity parameters of the rock, or vertical principal stress. <xref ref-type="bibr" rid="B5">Haimson and Fairhurst (1967)</xref> considered the diafiltration effect of the rock to improve the rupture pressure calculation formula, observing that the formation of cracks mainly depends on the stress state on the borehole wall (<xref ref-type="bibr" rid="B7">Zubkov et al., 2007</xref>; <xref ref-type="bibr" rid="B15">Schmidt and Zoback, 1989</xref>) considered the effect of porosity and Poisson&#x2019;s ratio on initiation pressure, which improved the Hubbert and Willis criteria, and proposed the initiation criteria of permeable rock formation and non-permeable rock formation. <xref ref-type="bibr" rid="B21">Yang et&#x20;al. (2014)</xref> pointed out the importance of supercritical CO<sub>2</sub> compressibility by studying the pore pressure of supercritical CO<sub>2</sub> wave propagation in saturated pore media. <xref ref-type="bibr" rid="B3">Chen Y. et&#x20;al. (2015)</xref> considered the characteristics of low viscosity, high compressibility, and the low-pressure rate of supercritical CO<sub>2</sub> fluid, and the supercharging rate model of supercritical CO<sub>2</sub> fluid in wellbore was established, thus concluding that the fracturing pressure of supercritical CO<sub>2</sub> is 75.5% lower than that of conventional hydraulic fracturing.</p>
<p>Current research on supercritical CO<sub>2</sub> fracturing is relatively rare. The deformation and failure characteristics of shale and the mechanism of crack propagation at the coupling of the original rock stress field and fluid field are key problems that need to be solved urgently in the process of shale gas exploitation. The original stress state, rock tensile strength, and development degree of natural crack are important factors that influence fracture morphology. Therefore, a three-axis fracturing indoor experiment was carried out to simulate supercritical CO<sub>2</sub> fracturing shale under different stress conditions, and acoustic emission detection technique and CT scanning were combined to investigate initial fracturing condition and analyze the process of crack formation, crack propagation, and crack morphological changes of fracturing shale under different stress conditions. The results of this research provide important parameters to reveal the mechanism that <italic>in situ</italic> stress states influence the law of shale crack propagation.</p>
</sec>
<sec id="s2">
<title>Test Equipment and Solutions</title>
<sec id="s2-1">
<title>Specimen Preparation</title>
<p>The shale specimens were taken from an outcrop of the lower Silurian Longmaxi formation, located in the Changning region of the Sichuan Basin. The total organic carbon (TOC) content was 4.18% and the vitrinite reflectance (Ro) was 2.36%, the uniaxial tensile strength was 9.2&#xa0;MPa, the elastic modulus was 16.0&#xa0;GPa, and Poisson&#x2019;s ratio was 0.31. The shale was mainly composed of quartz and illite, with the addition of calcite, feldspar, dolomite, chlorite, pyrite, and other minerals. The average value of the total amounts of brittle minerals was about 70.54%, which is similar to typical shale gas basins in North America.</p>
<p>To reduce the dispersion of the samples of the test results, they were drilled from the same piece of shale. The drilling direction was perpendicular to the bedding plane, and samples with smooth ends and no surface cracks were selected for processing into cylinder specimen (100&#xa0;mm in diameter and 200&#xa0;mm in length, axial and radial pressure were applied respectively), as shown in <xref ref-type="fig" rid="F1">Figure&#x20;1</xref>. Before the fracturing experiment, an axial central borehole with an 8-mm diameter was drilled to a depth of 140&#xa0;mm. An injection pipe was inserted into the borehole to ensure the tightness of the system during the experiment.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Diagram of drilling hole and preparation process of fracturing sample.</p>
</caption>
<graphic xlink:href="fenrg-09-763052-g001.tif"/>
</fig>
<p>The deviation of height and diameter in samples was controlled at &#xb1;0.3&#xa0;mm, and the deviation of flatness was controlled at &#xb1;0.02&#xa0;mm in the process. The shale samples were tested by a sonic wave before the fracturing experiment, the samples with relatively closer wave velocity were selected. The basic physical parameters of experimental shale samples are shown in <xref ref-type="table" rid="T1">Table&#x20;1</xref>.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Basic physical parameters of experimental samples.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Sample number</th>
<th align="center">Density (g.cm<sup>&#x2212;3</sup>)</th>
<th align="center">Porosity (%)</th>
<th align="center">Longitudinal velocity (m&#x387;s<sup>&#x2212;1</sup>)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">1(S<sub>1</sub>)</td>
<td align="char" char=".">2.59</td>
<td align="char" char=".">3.5569</td>
<td align="char" char=".">4,580</td>
</tr>
<tr>
<td align="left">2(S<sub>2</sub>)</td>
<td align="char" char=".">2.56</td>
<td align="char" char=".">3.6578</td>
<td align="char" char=".">4,540</td>
</tr>
<tr>
<td align="left">3(S<sub>3</sub>)</td>
<td align="char" char=".">2.52</td>
<td align="char" char=".">3.6925</td>
<td align="char" char=".">4,500</td>
</tr>
<tr>
<td align="left">4(S<sub>4</sub>)</td>
<td align="char" char=".">2.58</td>
<td align="char" char=".">3.5869</td>
<td align="char" char=".">4,600</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s2-2">
<title>Test Equipment</title>
<p>The supercritical CO<sub>2</sub> fracturing shale and increased permeability experimental device were independently developed by Chongqing University (<xref ref-type="fig" rid="F2">Figure&#x20;2A</xref>). The experimental device was specifically designed for supercritical CO<sub>2</sub> fracturing and was composed of six subsystems: an air supply system, CO<sub>2</sub> injection system, triaxial loading and control system, temperature control system (oil bath<bold>),</bold> acoustic emission monitoring system, and a data collecting system.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Schematic diagram of supercritical CO<sub>2</sub> fracturing shale test system. <bold>(A)</bold> High precision plunger pump. <bold>(B)</bold> Multi-field coupling experiment machine. <bold>(C)</bold> PCI-2 Acoustic emission. <bold>(D)</bold> CT (SOMTOM Scope).</p>
</caption>
<graphic xlink:href="fenrg-09-763052-g002.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="F2">Figure&#x20;2</xref> shows the test system. This test was developed to study the fracture law of shale caused by supercritical CO<sub>2</sub> and explore the influence law of supercritical CO<sub>2</sub> on shale initial cracking pressure, furthering research on self-developing supercritical CO<sub>2</sub> cracking using an increasing permeability test device (<xref ref-type="bibr" rid="B24">Zhou et&#x20;al., 2016</xref>), for shale gas reservoirs. The system is composed of six parts: Gas source supply system, CO<sub>2</sub> boosting control system, triaxial servo loading system, temperature control system, acoustic emission monitoring system, and data acquisition system. It can simulate an experiment of fracture propagation laws in the process of supercritical CO<sub>2</sub> fracturing shale under the coupling of multiple fields and collect the stress-strain relationship and acoustic emission signals in real time during the process of the shale fracture.<list list-type="simple">
<list-item>
<p>1) <xref ref-type="fig" rid="F2">Figure&#x20;2A</xref> is a High precision plunger pump, which is mainly used to provide high-pressure CO<sub>2</sub> fluid. The plunger pump can provide continuous, stable and pulse free flow, and can carry out constant pressure and constant flow operation. The maximum pressure is 50&#xa0;MPa and the maximum capacity is 266&#xa0;ml. Supercritical CO<sub>2</sub> is a critical state in which the temperature and pressure of CO<sub>2</sub> meet the requirements of 31.04&#xb0;C and 7.38&#xa0;mpa respectively. The plunger pump can fully meet the pressure control range and pressure mode required by the experimental process. The plunger pump adopts a constant pressure mode to&#x20;load.</p>
</list-item>
<list-item>
<p>2) <xref ref-type="fig" rid="F2">Figure&#x20;2B</xref> is a multi field coupling physical simulation test machine, which can realize the simulation test of real shale reservoir temperature, pressure, and stress field. After all the three-axis servo loading devices are lifted to the constant temperature oil bath, the stress is loaded when the temperature is constant to the set temperature. The pressure field mainly considers supercritical CO<sub>2</sub> fluid, and the gas inlet and outlet ends are all produced. The pressure sensor is used to monitor the change of fluid pressure during the experiment. The maximum load axial pressure of the system is 100&#xa0;MPa, the maximum confining pressure is 30&#xa0;MPa, the maximum fluid pressure can be provided at 20&#xa0;MPa, the maximum axial displacement is 60&#xa0;mm, the maximum circumferential deformation is 10&#xa0;mm, and the axial confining pressure adopts hydraulic loading, which can meet the requirements of the experiment. The transformation and upgrading of the test system provide a new means for the study of shale mechanical properties and fluid seepage law under the condition of multi field coupling, which can promote more in-depth study on the geological storage of carbon dioxide and the efficient exploitation of shale&#x20;gas.</p>
</list-item>
<list-item>
<p>3) <xref ref-type="fig" rid="F2">Figure&#x20;2C</xref> is a PCI-2 acoustic emission system, which has the advantages of computerization, full digitalization, multi-channel, low noise, low threshold value, ultra fast processing speed, and reliable stability. The 8-channel disp acoustic emission positioning system is used to monitor the fracture initiation and extension information, and the acoustic emission data, the event distribution location in the three-dimensional acoustic emission positioning, and the rock stress-strain curve are used to analyze the mechanical behavior of shale fracture propagation.</p>
</list-item>
<list-item>
<p>4) <xref ref-type="fig" rid="F2">Figure&#x20;2D</xref> is the medical CT produced by SIEMENS, which can perform non-destructive testing of&#x20;rocks.</p>
</list-item>
</list>
</p>
</sec>
<sec id="s2-3">
<title>Test Scheme</title>
<p>To reach the supercritical state of CO<sub>2</sub>, the experiment temperature was kept constant at 308&#xa0;K. To simulate the influence of different stress conditions on failure modes and initiation pressures of crack, horizontal stress remained constant, the initial stress ratio (the ratio of horizontal stress and vertical stress) was changed by changing vertical stress during the experiment. The confining pressure (horizontal stress) was constant at 10&#xa0;MPa, and the axial (vertical stress) was set at 10, 20, 30, and 40&#xa0;MPa, respectively. In the process of supercritical CO<sub>2</sub> fracturing, acoustic emission signals were monitored simultaneously, and the experimental scheme is shown in <xref ref-type="table" rid="T2">Table&#x20;2</xref>.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>The test parameters of supercritical CO<sub>2</sub> fracturing experiment.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Sample number</th>
<th align="center">Vertical stress (MPa)</th>
<th align="center">Horizontal stress (MPa)</th>
<th align="center">Initial stress ratio (&#x3bb;)</th>
<th align="center">Injection rate of CO<sub>2</sub> (ml<sup>.</sup>min<sup>&#x2212;1</sup>)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">1(S<sub>1</sub>)</td>
<td align="char" char=".">10</td>
<td align="char" char=".">10</td>
<td align="char" char=".">1</td>
<td align="char" char=".">20</td>
</tr>
<tr>
<td align="left">2(S<sub>2</sub>)</td>
<td align="char" char=".">20</td>
<td align="char" char=".">10</td>
<td align="char" char=".">0.5</td>
<td align="char" char=".">20</td>
</tr>
<tr>
<td align="left">3(S<sub>3</sub>)</td>
<td align="char" char=".">30</td>
<td align="char" char=".">10</td>
<td align="char" char=".">0.33</td>
<td align="char" char=".">20</td>
</tr>
<tr>
<td align="left">4(S<sub>4</sub>)</td>
<td align="char" char=".">40</td>
<td align="char" char=".">10</td>
<td align="char" char=".">0.25</td>
<td align="char" char=".">20</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Fracturing and re-fracturing experiments were carried out to simulate the fracturing and re-fracturing process. Therefore, supercritical CO<sub>2</sub> injection was divided into two stages, the injection rate of CO<sub>2</sub> fluid was constant at 20&#xa0;ml. min<sup>&#x2013;&#x2212;1</sup> at each stage. First, CO<sub>2</sub> was injected continuously until the pressure was no longer elevated, and after the completion of fracturing, when the pressure inside the specimen decreased, the pump pressure stopped immediately. Second, when the internal pressure of the crack achieved balance, CO<sub>2</sub> fluid continued to be injected.</p>
<p>Before and after the fracturing test, the specimens were photographed using a high definition digital camera and CT scanning to obtain surface and internal crack characteristics of the samples before and after fracturing. The CT scans of each specimen, along with samples in a vertical direction (200&#xa0;mm) of 0.75&#xa0;mm thickness to form a&#x20;layer.</p>
</sec>
</sec>
<sec id="s3">
<title>Analysis of Experimental Results</title>
<sec id="s3-1">
<title>Analysis of Pump Pressure Curve</title>
<p>The acoustic emission count-time curve and pump pressure-time curve during shale fracturing experiments are shown in <xref ref-type="fig" rid="F3">Figure&#x20;3</xref>. The overall trend of the pump pressure curve and acoustic emission curve of the four specimens was similar. The process can be divided into three stages:</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>The relationship curves between acoustic emission count-time and pump pressure-time of samples. <bold>(A)</bold> The relationship curve between acoustic emission count-time and pump pressure-time of S<sub>1</sub>. <bold>(B)</bold> The relationship curve between acoustic emission count-time and pump pressure-time of S<sub>2</sub>. <bold>(C)</bold> The relationship curve between acoustic emission count-time and pump pressure-time of S<sub>3</sub>. <bold>(D)</bold> The relationship curve between acoustic emission count-time and pump pressure-time of S<sub>4</sub>
</p>
</caption>
<graphic xlink:href="fenrg-09-763052-g003.tif"/>
</fig>
</sec>
<sec id="s3-2">
<title>The Stage of Pressure Accumulation in Wellbore</title>
<p>The initial pressure of failure of shale can be defined as peak pressure at this stage. Carbon dioxide fluid, an effective fracturing fluid, first gathered in an artificial cavity pocket with the gradual injection of fluid, the pressure of carbon dioxide in the fracturing chamber grew sharply, before the carbon dioxide fluid started to squeeze into the pore and fissure system around the wellbore. When the pressure reached a certain value, the pressure of fracturing fluid growth rate became relatively slow due to the diffusion of pressure in cracks and pores. The small amount of acoustic emission signals produced in this phase indicate that the native cracks of shale began to open gradually under the effect of carbon dioxide fluid along with the formation of micro cracks, but shale had not yet formed macroscopic cracks through this process. Carbon dioxide fluid pressure still rose and its increase rate began to increase faster after the pore and fissure space was full of carbon dioxide with a continuous injection of fluid until the shale failed, and pump pressure vs. time curve reached the first peak value. A large number of acoustic emission counts were found during the peak value of pressure, which indicates that fractures were generated and extended dramatically, and transfixion cracks were formed.</p>
</sec>
<sec id="s3-3">
<title>Stage of Pump Closure</title>
<p>The pressure of fracturing fluid decreased rapidly after the closure of the fracturing pump, as crack coalescence formed in the shale specimens and the seepage channel had been generated, which resulted in the faster spread of the fluid inside the crack and pressure dropping quickly in the fracture until finally, the pressure remained at a certain value. The acoustic emission signal was weak and the counts were relatively small (or even none), which showed that hardly any new cracks were generated at this stage and mainly show the closure of the generated fracture.</p>
</sec>
<sec id="s3-4">
<title>Crack Reopening Stage</title>
<p>The fracturing fluid of supercritical carbon dioxide was continually injected into the opened cracks by opening the plunger pump after the pressure of the fracturing chamber reached the balance value. The pressure of the fracture increased dramatically at the same time, the fractured cracks reopened after the fracture stop closure, and crack pressure reached the second peak value, which reached Pr. The AE counts of repeated fracturing reduced significantly compared to the initial fracturing, but it still had a strong acoustic emission signal, which indicates that the original cracks continue to propagate and extend at the second fracturing process and many secondary cracks were generated.</p>
<p>The results of supercritical CO<sub>2</sub> fluid fracturing shale are shown in <xref ref-type="table" rid="T3">Table&#x20;3</xref>. In combination with <xref ref-type="table" rid="T1">Table&#x20;1</xref> and <xref ref-type="table" rid="T2">Table&#x20;2</xref>, it can be concluded that the initiation pressure of shale is always higher than horizontal stress; however, the stress of repeated fracturing is close to horizontal stress (10&#xa0;MPa). With the decrease of &#x3bb;, the initiation pressure of shale, the stress of repeated fracturing, and the amount of fracturing fluid decreased. Among them, the initiation pressure of Specimen 1 was highest (13.7&#xa0;MPa), and the stress of repeated fracturing on Specimen 4 was the smallest (9.4&#xa0;MPa).</p>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>Experimental results of shale specimens.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Sample number</th>
<th align="center">Breaking time (s)</th>
<th align="center">Initiation pressure (MPa)</th>
<th align="center">Reopening pressure (MPa)</th>
<th align="center">Consumed fluid volume (ml)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">1(S<sub>1</sub>)</td>
<td align="char" char=".">1,025</td>
<td align="char" char=".">13.7</td>
<td align="char" char=".">10.6</td>
<td align="char" char=".">1,192</td>
</tr>
<tr>
<td align="left">2(S<sub>2</sub>)</td>
<td align="char" char=".">1,000</td>
<td align="char" char=".">11.9</td>
<td align="char" char=".">10.5</td>
<td align="char" char=".">1,022</td>
</tr>
<tr>
<td align="left">3(S<sub>3</sub>)</td>
<td align="char" char=".">836</td>
<td align="char" char=".">11.3</td>
<td align="char" char=".">9.8</td>
<td align="char" char=".">825</td>
</tr>
<tr>
<td align="left">4(S<sub>4</sub>)</td>
<td align="char" char=".">985</td>
<td align="char" char=".">10.7</td>
<td align="char" char=".">9.4</td>
<td align="char" char=".">810</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>According to first strength theory, during the supercritical CO<sub>2</sub> fluid fracturing process when the fluid pressure exceeds the stress required for the cracking of the rock at the hole wall, cracks are generated at the wall of the hole, and the occurrence of cracks is related to the type of fluid, the size of ground stress and the physical and mechanical properties of the rock. In this paper, with a constant confining pressure and the greater the vertical stress, the smaller the initial stress ratio and the lower the stress of the rock cracking at the hole wall, which results in a decrease in the initiation pressure with the decrease of the initial stress ratio (<xref ref-type="bibr" rid="B23">Zhao et&#x20;al., 2015</xref>). The relationship between the initiation pressure and the initial stress ratio of the shale specimens is shown in <xref ref-type="fig" rid="F4">Figure&#x20;4</xref>, which shows that the initiation pressure increases with the increase of the initial stress ratio and that there is a linear relationship between the initiation pressure and the initial stress ratio, which can be fitted to y &#x3d; 3.7x &#x2b; 10 (<italic>R</italic>
<sup>2</sup> &#x3d; 0.998).</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>The fitting relationship between crack initiation pressure and initial stress&#x20;ratio.</p>
</caption>
<graphic xlink:href="fenrg-09-763052-g004.tif"/>
</fig>
<p>When the shale samples initial pressure relief, the samples will reach the second pressure peak, continuing to inject CO<sub>2</sub>. The reopening pressure can be obtained. As shown in <xref ref-type="table" rid="T3">Table&#x20;3</xref>. The reopening pressure of the cracks are lower than the initial pressure of the crack, which is due to damage to the shale specimen and forms a certain crack network during the initial fracturing as the strength of the shale itself was reduced. In the second fracturing, the crack continued to propagate based on the existing crack. This can also be verified in the CT scanning image where the initial stress ratio was different, as the number of fractures formed and the crack volume space were different in terms of fracturing. The different crack volumes led to a difference in the amount of CO<sub>2</sub> fluid entering the crack pore space, so different initial stress ratios in addition to the impact of initiation pressure also lead to a difference in the amount of fracturing&#x20;fluid.</p>
</sec>
<sec id="s3-5">
<title>CT Scanning and the Crack Characteristics of Shale Samples</title>
<p>From the obtained CT scanning images of the shale samples, each specimen (S<sub>1</sub>&#x2013;S<sub>4</sub>) was selected from two groups of CT scanning images in the vertical and parallel direction of the axial section before and after fracturing, with a total of 16 images shown in <xref ref-type="fig" rid="F5">Figure&#x20;5</xref>.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>CT scanning before and after the fracturing of four samples. <bold>(A)</bold> CT scanning image before and after fracturing of sample 1. <bold>(B)</bold> CT scanning image before and after fracturing of sample 2. <bold>(C)</bold> CT scanning image before and after fracturing of sample 3. <bold>(D)</bold> CT scanning image before and after fracturing of sample 4.</p>
</caption>
<graphic xlink:href="fenrg-09-763052-g005.tif"/>
</fig>
<p>The CT scanning image of the initial stress ratio at 1 is shown in <xref ref-type="fig" rid="F5">Figure&#x20;5A</xref>. Some new cracks were created at the level of the vertical wellbore axis, while transverse cracks were found at the level of the parallel wellbore axis, which suggests that cracks mainly extend from the center of the wellbore to the horizontal direction of the wellbore, at the same vertical stress and horizontal stress (<xref ref-type="fig" rid="F5">Figure&#x20;5A</xref>).</p>
<p>In <xref ref-type="fig" rid="F5">Figures 5B&#x2013;D</xref>, there were no obvious cracks in the shale specimens and the shales kept integrity before fracturing. In comparison to <xref ref-type="fig" rid="F5">Figure&#x20;5A</xref>, some cracks were detected in both the vertical and parallel wellbore sections after fracturing in (<xref ref-type="fig" rid="F5">Figures 5B&#x2013;D</xref>), and cracks and fractures were similar to the end face of the specimens, which indicates that cracks spread from top and bottom, and the cracks in the specimens were connected. Furthermore, it was noted that the cracks mainly occurred along the direction of maximum principal stress, thus, cracks in the specimens after fracturing confirmed that supercritical CO<sub>2</sub> fracturing can form a complex network of cracks in shale specimens.</p>
<p>To better understand the process of supercritical CO<sub>2</sub> fracturing, the mechanism of the fracturing is described. From <xref ref-type="fig" rid="F5">Figures 5B&#x2013;D</xref>, we can see that cracks were first produced at the wall of the hole, and tensile failure occurred along the direction of maximum principal stress and accompanied by shear failure. The mechanical failure of the fracturing may relate to the developed stratification of shale and the characteristics of supercritical CO<sub>2</sub> with low viscosity and high diffusivity. Compared to water, the low viscosity and high diffusion of supercritical CO<sub>2</sub> make it easier to enter the micropore of the shale and make the cracks extend and turn along the weak bedding surface. These research findings are consistent with those of <xref ref-type="bibr" rid="B2">Chen L. Q. et&#x20;al. (2015)</xref>, who obersved that, due to the special properties of supercritical CO<sub>2</sub> fracturing, it forms a complex fracturing networks more easily. In the process of crack propagation&#x2014;due to the increase of free space&#x2014;the phase change of CO<sub>2</sub> from supercritical states to gas states was likely to have occurred. In addition, the compressibility of gaseous CO<sub>2</sub> is much larger than supercritical CO<sub>2,</sub> so the phase change will produce gas explosions steeply, and is favorable for the production of complex fracture networks in shale specimens.</p>
<p>As shown in <xref ref-type="fig" rid="F5">Figures 5B&#x2013;D</xref>, it is interesting to note that new cracks at the same height of the shale specimens varied from different initial stress ratios. Specimen 2 had the largest number of cracks and the most complex fracture networks, while there were fewer cracks and fractures in Specimen 4, which had monotonous fracture networks. Under different stress state conditions, the crack initiation pressure and propagation direction varies with stress states and eventually forms different crack distributions and quantities, which show that the complexity of fracturing networks is affected by the initial stress ratio. In conclusion, these results show that cracks extend along the direction of maximum principal stress, and the cracks in shale form the main crack parallel to the direction of maximum principal stress after fracturing by supercritical CO<sub>2</sub>&#x20;fluid.</p>
</sec>
</sec>
<sec sec-type="discussion" id="s4">
<title>Discussion</title>
<p>There are many theories relating to obtaining the pressure of crack initiation. One of the most classic formulae was proposed by <xref ref-type="bibr" rid="B8">Hubbert and Willis (1972)</xref>. In this paper, according to the first strength theory, the initial rupture usually occurs (<xref ref-type="bibr" rid="B23">Zhao et&#x20;al., 2015</xref>; <xref ref-type="bibr" rid="B22">Yang et&#x20;al., 2019</xref>) at the hole wall. The tangential pressure at the hole wall may be defined as:<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>&#x3c6;</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mi>&#x3bb;</mml:mi>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:mn>2</mml:mn>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mi>&#x3bb;</mml:mi>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3c6;</mml:mi>
<mml:mo>-</mml:mo>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>Where <inline-formula id="inf1">
<mml:math id="m2">
<mml:mi>P</mml:mi>
</mml:math>
</inline-formula> is the injection pressure; <inline-formula id="inf2">
<mml:math id="m3">
<mml:mi>&#x3bb;</mml:mi>
</mml:math>
</inline-formula> is the initial stress ratio coefficient; <inline-formula id="inf3">
<mml:math id="m4">
<mml:mi>&#x3b2;</mml:mi>
</mml:math>
</inline-formula> is the seepage coefficient; <inline-formula id="inf4">
<mml:math id="m5">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mtext>c</mml:mtext>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the horizontal stress.</p>
<p>The tensile strength of shale is far less than its compressive strength and the conditions for tensile failure or shear failure, so long as the tangential pressure reaches the tensile strength.<disp-formula id="e2">
<mml:math id="m6">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>&#x3c6;</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>-</mml:mo>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mtext>t</mml:mtext>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>Where <inline-formula id="inf5">
<mml:math id="m7">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mtext>t</mml:mtext>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the tensile strength of the shale, the value of tensile stress is negative. The <xref ref-type="disp-formula" rid="e1">Eq. 1</xref> taken into the <xref ref-type="disp-formula" rid="e2">Eq. 2</xref> can be obtained:<disp-formula id="e3">
<mml:math id="m8">
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mi>&#x3bb;</mml:mi>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mtext>c</mml:mtext>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:mn>2</mml:mn>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mi>&#x3bb;</mml:mi>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mtext>c</mml:mtext>
</mml:msub>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3c6;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mtext>t</mml:mtext>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
</p>
<p>Because <inline-formula id="inf6">
<mml:math id="m9">
<mml:mrow>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3c6;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> &#x2208; [&#x2212;1, 1], when shale is most likely to damage at <inline-formula id="inf7">
<mml:math id="m10">
<mml:mrow>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3c6;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; 1 or &#x3d; &#x2212;1 (<inline-formula id="inf8">
<mml:math id="m11">
<mml:mi>&#x3c6;</mml:mi>
</mml:math>
</inline-formula> &#x3d; 0&#xb0;or <inline-formula id="inf9">
<mml:math id="m12">
<mml:mi>&#x3c6;</mml:mi>
</mml:math>
</inline-formula> &#x3d; 90&#xb0;), the initiation pressure <inline-formula id="inf10">
<mml:math id="m13">
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mtext>b</mml:mtext>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> of shale can be obtained:<disp-formula id="e4">
<mml:math id="m14">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mtext>b</mml:mtext>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mtext>min</mml:mtext>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mtext>t</mml:mtext>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mn>3</mml:mn>
<mml:mi>&#x3bb;</mml:mi>
</mml:mfrac>
<mml:mo>-</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mtext>c</mml:mtext>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mtext>t</mml:mtext>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>3</mml:mn>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mi>&#x3bb;</mml:mi>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mtext>c</mml:mtext>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>}</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
</p>
<p>When &#x03BB; &#x3d; 1, there are:<disp-formula id="e5">
<mml:math id="m15">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mtext>b</mml:mtext>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mtext>t</mml:mtext>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>2</mml:mn>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mtext>c</mml:mtext>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>
</p>
<p>It can be obtained according to this formula: <inline-formula id="inf11">
<mml:math id="m16">
<mml:mrow>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3c6;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; 1 (<inline-formula id="inf12">
<mml:math id="m17">
<mml:mi>&#x3c6;</mml:mi>
</mml:math>
</inline-formula> &#x3d; 0&#xb0;) or <inline-formula id="inf13">
<mml:math id="m18">
<mml:mrow>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3c6;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; &#x2212;1 (<inline-formula id="inf14">
<mml:math id="m19">
<mml:mi>&#x3c6;</mml:mi>
</mml:math>
</inline-formula> &#x3d; 90&#xb0;). At this point, there are two possible directions for crack initiation, one is horizontal and the other is vertical, the experimental results belong to the theory.</p>
<p>When <inline-formula id="inf15">
<mml:math id="m20">
<mml:mi>&#x3bb;</mml:mi>
</mml:math>
</inline-formula> &#x3c; 1, there are:<disp-formula id="e6">
<mml:math id="m21">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mtext>t</mml:mtext>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>3</mml:mn>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mi>&#x3bb;</mml:mi>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mtext>c</mml:mtext>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3c;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x3c3;</mml:mi>
<mml:mtext>t</mml:mtext>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mn>3</mml:mn>
<mml:mi>&#x3bb;</mml:mi>
</mml:mfrac>
<mml:mo>-</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mtext>c</mml:mtext>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>
<disp-formula id="e7">
<mml:math id="m22">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mtext>b</mml:mtext>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mtext>t</mml:mtext>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>3</mml:mn>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mi>&#x3bb;</mml:mi>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mtext>c</mml:mtext>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
</p>
<p>At the same time, the fracture initiation will occur in the vertical direction of the fracture hole and can be concluded from the <xref ref-type="disp-formula" rid="e7">Eq. 7</xref>. The smaller the initial stress ratio is, the&#x20;smaller the initiation pressure. At the same time, the initial stress ratio is different, the stress distribution in the wellbore is different, and the location of the crack will be different so that the direction of crack propagation will not be the same, so the crack shape will be different. In short, the reliability of the experimental results is further verified by the theory.</p>
<p>The crack initiation orientation is affected by the maximum principal stress characteristics, the cracks are firstly initiated in the wellbore. The occurrence of expansion along with the maximum principal stress, eventually form a main crack in the vertical direction of minimum principal stress. Stress state will influence crack initiation orientation, crack morphology, crack distribution and crack initiation pressure. The study of the extraction and analysis of rock fracturing crack propagation through supercritical CO<sub>2</sub> fracturing shale of low permeability under different triaxial stress conditions can provide a reference for reasonably selecting the fracturing layer of the shale gas reservoir and design of CO<sub>2</sub> geological storage.</p>
<p>
<xref ref-type="bibr" rid="B16">Shang, et&#x20;al. (2018)</xref> used the local asymmetric fracturing stimulation method with small scale high-pressure fluids such as carbon dioxide, liquid nitrogen, water, and jet cutting, etc. to create slots and perforations for the aid of fracture initiation and orientation. The study results show that the asymmetric fracturing with slot and perforation orientation significantly increased fracture complexity. The research in this paper, in terms of the initial stress ratio and bedding effect of the crack propagation morphology and development process of rock, is consistent with the conclusions of <xref ref-type="bibr" rid="B16">Shang et&#x20;al. (2018)</xref>.</p>
<p>Compared with hydraulic fracturing, supercritical CO<sub>2</sub> has low viscosity, meaning that it easily enters the very small pores and fractures in the shale specimen. The crack propagation process is not only more important for the microstructure and micro defects of shale specimens, but also more sensitive and helpful for the penetration of micro defects. The low viscosity characteristics of supercritical CO<sub>2</sub> are helpful to crack propagation. In addition, during the process of crack propagation, the supercritical CO<sub>2</sub> may also change from a supercritical state to a gaseous state. The compressibility of CO<sub>2</sub> gas is so much greater than the supercritical state, that a sharp expansion of gas is produced during the process of phase transition. This is favorable to a weak surface or crack in the shale specimen. Therefore, the supercritical CO<sub>2</sub> fracturing of shale simulation test method and test results established in this paper provides a new technical path for the development of anhydrous fracturing of shale&#x20;gas.</p>
</sec>
<sec sec-type="conclusion" id="s5">
<title>Conclusion</title>
<p>Supercritical dioxide carbon fracturing experiments were carried by advanced multi-filed equipment under the condition of different stress states. Acoustic emission signal monitoring and crack CT scanning were used to analyze the fracturing crack propagation. The conclusions of this study are as follows:<list list-type="simple">
<list-item>
<p>1) When vertical stress is equal to horizontal stress, fracturing cracks propagate along the horizontal direction. When the vertical stress is greater than the horizontal stress, fracturing cracks propagate along the vertical extension. Supercritical CO<sub>2</sub> fracturing can form a main crack in triaxial stress conditions in shale, with the main crack expansion occurring along the direction of maximum principal stress expansion.</p>
</list-item>
<list-item>
<p>2) Shale crack initiation pressure is always greater than horizontal stress, and reopening pressure is lower than initiation pressure. The initial stress ratio influences crack initiation and propagation pressure. Stress state can influence crack initiation orientation, crack morphology, crack distribution and crack initiation pressure. With the decrease of the stress ratio, the crack becomes less obvious.</p>
</list-item>
<list-item>
<p>3) The supercritical CO<sub>2</sub> fracturing process can be divided into three representative stages: (a) the stage of pressure accumulation in the wellbore; (b) the stage of pump closure; and (c) the crack reopening stage. At the crack reopening stage, the crack reopens again and continues to extend to form a complex fracture network. Acoustic emission signal intensity directly reflects shale crack initiation, propagation, and the coalescence process.</p>
</list-item>
</list>
</p>
</sec>
</body>
<back>
<sec id="s6">
<title>Data Availability Statement</title>
<p>The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.</p>
</sec>
<sec id="s7">
<title>Author Contributions</title>
<p>GL: Conceptualization, Writing - Original Draft, Funding acquisition, YZ: Resources, Formal analysis, Visualization, YC: Investigation, Visualization, Funding acquisition.</p>
</sec>
<sec id="s8">
<title>Funding</title>
<p>This work was supported by the National Natural Science Foundation of China (No. 51904049); the China Postdoctoral Science Foundation (2021M693750); the Natural Science Foundation of Chongqing (General Program, No. cstc2019jcyj-msxmX0702); the Program for Natural Science Foundation of Hunan Province (No. 2021JJ40028), Research Program of Department of Education in Hunan Provincial (No. 19C1781).</p>
</sec>
<sec sec-type="COI-statement" id="s9">
<title>Conflict of Interest</title>
<p>YZ was employed by the company Sinohydro Bureau 8 Co. LTD.POWERCHINA.</p>
<p>The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="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>Brantley</surname>
<given-names>S. L.</given-names>
</name>
<name>
<surname>Yoxtheimer</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Arjmand</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Grieve</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Vidic</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Pollak</surname>
<given-names>J.</given-names>
</name>
<etal/>
</person-group> (<year>2014</year>). <article-title>Water Resource Impacts During Unconventional Shale Gas Development: the Pennsylvania Experience</article-title>. <source>Int. J.&#x20;Coal Geology.</source> <volume>126</volume>, <fpage>140</fpage>&#x2013;<lpage>156</lpage>. <pub-id pub-id-type="doi">10.1016/j.coal.2013.12.017</pub-id> </citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>L. Q.</given-names>
</name>
<name>
<surname>Tian</surname>
<given-names>S. C.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>G. S.</given-names>
</name>
<name>
<surname>Fan</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2015a</year>). <article-title>Initiation Pressure Models for Supercritical CO2 Fracturing and Sensitivity Analysis</article-title>. <source>J.&#x20;Rock Soil Mech.</source> <volume>36</volume>, <fpage>125</fpage>&#x2013;<lpage>131</lpage>. <pub-id pub-id-type="doi">10.16285/j.rsm.2015.S2.016</pub-id> </citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Nagaya</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Ishida</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2015b</year>). <article-title>Observations of Fractures Induced by Hydraulic Fracturing in Anisotropic Granite</article-title>. <source>Rock Mech. Rock Eng.</source> <volume>48</volume> (<issue>4</issue>), <fpage>1455</fpage>&#x2013;<lpage>1461</lpage>. <pub-id pub-id-type="doi">10.1007/s00603-015-0727-9</pub-id> </citation>
</ref>
<ref id="B4">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Gupta</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Gupta</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Langlinais</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2005</year>). <source>Feasibility of Supercritical Carbon Dioxide as a Drilling Fluid for Deep Underbalanced Drilling Operation</source>. <publisher-loc>Dallas, Texas</publisher-loc>: <publisher-name>Society of Petroleum Engineers</publisher-name>. <pub-id pub-id-type="doi">10.2118/96992-MS</pub-id> </citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Haimson</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>FaIrhurst</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>1967</year>). <article-title>Initiation and Extension of Hydraulic Fractures in Rocks</article-title>. <source>Soc. Pet. Eng. J.</source> <volume>7</volume> (<issue>6</issue>), <fpage>310</fpage>&#x2013;<lpage>318</lpage>. <pub-id pub-id-type="doi">10.2118/1710-pa</pub-id> </citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huang</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Lu</surname>
<given-names>Y. Y.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>J.&#x20;R.</given-names>
</name>
<name>
<surname>Ao</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>L. W.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>On the Erosion of Shale Impacted by Supercritical Carbon Dioxide jet.[J]</article-title>. <source>J.&#x20;Rock Mech. Eng.</source> <volume>34</volume> (<issue>4</issue>), <fpage>787</fpage>&#x2013;<lpage>794</lpage>. <pub-id pub-id-type="doi">10.13722/j.cnki.jrme.2015.04.016</pub-id> </citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hubbert</surname>
<given-names>M. K.</given-names>
</name>
<name>
<surname>Willis</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>1972</year>). <article-title>Mechanics of Hydraulic Fracturing [j]</article-title>. <source>Trans. AIME.</source> <volume>18</volume> (<issue>1</issue>), <fpage>369</fpage>&#x2013;<lpage>390</lpage>. <pub-id pub-id-type="doi">10.1080/14786435708241195</pub-id> </citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ishida</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Aoyagi</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Niwa</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Acoustic Emission Monitoring of Hydraulic Fracturing Laboratory Experiment With Supercritical and Liquid CO<sub>2</sub>[J]</article-title>. <source>Geophys. Res. Lett.</source> <volume>39</volume> (<issue>16</issue>). <pub-id pub-id-type="doi">10.1029/2012gl052788</pub-id> </citation>
</ref>
<ref id="B10">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Jin</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Shah</surname>
<given-names>S. N.</given-names>
</name>
<name>
<surname>Mao</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2012</year>). &#x201C;<article-title>Hydraulic Fracturing Model Based on Nonlinear Fracture Mechanics: Theory and Simulation</article-title>,&#x201D; in <conf-name>SPE Annual Technical Conference and Exhibition</conf-name>, <conf-loc>San Antonio, Texas</conf-loc>. </citation>
</ref>
<ref id="B11">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Jin</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Shah</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2013</year>). &#x201C;<article-title>Fracture Propagation Direction and its Application in Hydraulic Fracturing</article-title>,&#x201D; <conf-name>SPE Hydraulic Fracturing Technology Conference</conf-name>, <conf-loc>Texas</conf-loc>. <pub-id pub-id-type="doi">10.2118/163832-MS</pub-id> </citation>
</ref>
<ref id="B12">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Kolle</surname>
<given-names>J.&#x20;J.</given-names>
</name>
</person-group> (<year>2000</year>). &#x201C;<article-title>Coiled-Tubing Drilling With Supercritical Carbon Dioxide</article-title>,&#x201D; in <conf-name>SPE/Petroleum Society of CIM International Conference on Horizontal Well Technology</conf-name>, <conf-loc>Calgary</conf-loc>. <pub-id pub-id-type="doi">10.2118/65534-MS</pub-id> </citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lai</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Ye</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>LiuXiangjunYang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>In House Study on Water Locking Damage of Tight Sand Gas Reservoirs With Low Permeability [J]</article-title>. <source>Nat. Gas Industry.</source> <volume>25</volume> (<issue>4</issue>), <fpage>125</fpage>&#x2013;<lpage>127</lpage>. <pub-id pub-id-type="doi">10.1016/j.molcatb.2005.02.001</pub-id> </citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Middleton</surname>
<given-names>R. S.</given-names>
</name>
<name>
<surname>Carey</surname>
<given-names>J.&#x20;W.</given-names>
</name>
<name>
<surname>Currier</surname>
<given-names>R. P.</given-names>
</name>
<name>
<surname>Hyman</surname>
<given-names>J.&#x20;D.</given-names>
</name>
<name>
<surname>Kang</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Karra</surname>
<given-names>S.</given-names>
</name>
<etal/>
</person-group> (<year>2015</year>). <article-title>Shale Gas and Non-Aqueous Fracturing Fluids: Opportunities and Challenges for Supercritical CO2</article-title>. <source>Appl. Energ.</source> <volume>147</volume>, <fpage>500</fpage>&#x2013;<lpage>509</lpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2015.03.023</pub-id> </citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Schmidt</surname>
<given-names>D. R.</given-names>
</name>
<name>
<surname>Zoback</surname>
<given-names>M. D.</given-names>
</name>
</person-group> (<year>1989</year>). <article-title>Poroelasticity Effects in the Determination of Minimum Horizontal Principal Stress in Hydraulic Fracturing Test&#x2014;A Proposed Breakdown Equation Employing a Modified Effective Stress Relation for Tensile Failure[J]</article-title>. <source>Int. J.&#x20;Rock Mech. Mining Sci. Geomechanics Abstr.</source> <volume>26</volume> (<issue>6</issue>), <fpage>499</fpage>&#x2013;<lpage>506</lpage>. <pub-id pub-id-type="doi">10.1016/0148-9062(89)91427-7</pub-id> </citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shang</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Yin</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Deng</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Kang</surname>
<given-names>X.</given-names>
</name>
<etal/>
</person-group> (<year>2018</year>). <article-title>Local Asymmetric Fracturing to Construct Complex Fracture Network in Tight Porous Reservoirs During Subsurface Coal Mining: An Experimental Study</article-title>. <source>J.&#x20;Nat. Gas Sci. Eng.</source> <volume>59</volume>, <fpage>343</fpage>&#x2013;<lpage>353</lpage>. <pub-id pub-id-type="doi">10.1016/j.jngse.2018.09.005</pub-id> </citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Small</surname>
<given-names>M. J.</given-names>
</name>
<name>
<surname>Stern</surname>
<given-names>P. C.</given-names>
</name>
<name>
<surname>Bomberg</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Christopherson</surname>
<given-names>S. M.</given-names>
</name>
<name>
<surname>Goldstein</surname>
<given-names>B. D.</given-names>
</name>
<name>
<surname>Israel</surname>
<given-names>A. L.</given-names>
</name>
<etal/>
</person-group> (<year>2014</year>). <article-title>Risks and Risk Governance in Unconventional Shale Gas Development</article-title>. <source>Environ. Sci. Technol.</source> <volume>48</volume> (<issue>15</issue>), <fpage>8289</fpage>&#x2013;<lpage>8297</lpage>. <pub-id pub-id-type="doi">10.1021/es502111u</pub-id> </citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Vengosh</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Jackson</surname>
<given-names>R. B.</given-names>
</name>
<name>
<surname>Warner</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Darrah</surname>
<given-names>T. H.</given-names>
</name>
<name>
<surname>Kondash</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>A Critical Review of the Risks to Water Resources From Unconventional Shale Gas Development and Hydraulic Fracturing in the united states</article-title>. <source>Environ. Sci. Technol.</source> <volume>48</volume> (<issue>15</issue>), <fpage>8334</fpage>&#x2013;<lpage>8348</lpage>. <pub-id pub-id-type="doi">10.1021/es405118y</pub-id> </citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Feasibility Analysis on Shale Gas Exploitation With Supercritical Co2</article-title>. <source>Pet. Drilling Tech.</source> <volume>39</volume> (<issue>3</issue>), <fpage>30</fpage>&#x2013;<lpage>35</lpage>. <pub-id pub-id-type="doi">10.3969/j.issn.1001-0890.2011.03.005</pub-id> </citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xie</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Fang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Shi</surname>
<given-names>L.</given-names>
</name>
<etal/>
</person-group> (<year>2014</year>). <article-title>Carbon Geological Utilization and Storage in China: Current Status and Perspectives</article-title>. <source>Acta Geotech.</source> <volume>9</volume> (<issue>1</issue>), <fpage>7</fpage>&#x2013;<lpage>27</lpage>. <pub-id pub-id-type="doi">10.1007/s11440-013-0277-9</pub-id> </citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname>
<given-names>D. X.</given-names>
</name>
<name>
<surname>Qi</surname>
<given-names>L. I.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Numerical Analysis of Propagation of Pore Pressure Waves in Compressible Fluid Saturated Porous Media</article-title>. <source>Rock Soil Mech.</source> <volume>35</volume> (<issue>7</issue>), <fpage>2047</fpage>&#x2013;<lpage>2056</lpage>. <pub-id pub-id-type="doi">10.16285/j.rsm.2014.07.001</pub-id> </citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhang.</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Du</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Shang</surname>
<given-names>D.</given-names>
</name>
<etal/>
</person-group> (<year>2019</year>). <article-title>Supercritical Co2 Fracturing With Different Drilling Depths in Shale</article-title>. <source>Energy Sourc. A: Recov., Util., Environ. Eff.</source> <volume>18</volume> (<issue>3</issue>), <fpage>1</fpage>&#x2013;<lpage>20</lpage>. <pub-id pub-id-type="doi">10.1080/15567036.2019.1673850</pub-id> </citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Cao</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Qin</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>The Analysis of Antireflection Range in Coal Seam Hydraulic Fracturing</article-title>. <source>J.&#x20;Mining Saf. Eng.</source> <volume>32</volume> (<issue>4</issue>), <fpage>644</fpage>&#x2013;<lpage>650</lpage>. <pub-id pub-id-type="doi">10.13545/j.cnki.jmse.2015.04.020</pub-id> </citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhou</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Xian</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>D.</given-names>
</name>
<etal/>
</person-group> (<year>2016</year>). <article-title>Supercritical Carbon Dioxide Fracturing in Shale and the Coupled Effects on the Permeability of Fractured Shale: An Experimental Study</article-title>. <source>J.&#x20;Nat. Gas Sci. Eng.</source> <volume>36</volume>, <fpage>369</fpage>&#x2013;<lpage>377</lpage>. <pub-id pub-id-type="doi">10.1016/j.jngse.2016.10.005</pub-id> </citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zubkov</surname>
<given-names>V. V.</given-names>
</name>
<name>
<surname>Koshelev</surname>
<given-names>V. F.</given-names>
</name>
<name>
<surname>Lin&#x2019;kov</surname>
<given-names>A. M.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Numerical Modeling of Hydraulic Fracture Initiation and Development</article-title>. <source>J. Min. Sci.</source> <volume>43</volume>, <fpage>40</fpage>&#x2013;<lpage>56</lpage>. <pub-id pub-id-type="doi">10.1007/s10913-007-0006-6</pub-id> </citation>
</ref>
</ref-list>
</back>
</article>