<?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">864527</article-id>
<article-id pub-id-type="doi">10.3389/fenrg.2022.864527</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>Research on Refracturing Technology of Horizontal Wells Based on Dynamic Drainage Volume</article-title>
<alt-title alt-title-type="left-running-head">Dong et al.</alt-title>
<alt-title alt-title-type="right-running-head">Refracturing Technology of Horizontal Wells</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Dong</surname>
<given-names>Qi</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>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Li</surname>
<given-names>Jianshan</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Yang</surname>
<given-names>Lian</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Wang</surname>
<given-names>Fei</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhao</surname>
<given-names>Kun</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1655358/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Huo</surname>
<given-names>Fei</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Oil and Gas Technology Research Institute ChangQing Oilfield Company</institution>, <addr-line>Xi&#x2019;an</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>State Engineering Laboratory for Exploration and Development of Low Permeability Oil and Gas Fields</institution>, <addr-line>Xi&#x2019;an</addr-line>, <country>China</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>Key Laboratory of Unconventional Oil and Gas Development</institution>, <institution>Ministry of Education</institution>, <addr-line>Qingdao</addr-line>, <country>China</country>
</aff>
<aff id="aff4">
<sup>4</sup>
<institution>School of Petroleum Engineering</institution>, <institution>China University of Petroleum</institution>, <addr-line>Qingdao</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/1389365/overview">Qi Zhang</ext-link>, China University of Geosciences Wuhan, 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/1666605/overview">Mingshan Zhang</ext-link>, University of Alberta, Canada</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1468723/overview">Debin Kong</ext-link>, University of Science and Technology Beijing, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1650286/overview">Shiyuan Zhan</ext-link>, Chengdu University of Technology, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Qi Dong, <email>dqi_cq@petrochina.com.cn</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Advanced Clean Fuel Technologies, a section of the journal Frontiers in Energy Research</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>03</day>
<month>05</month>
<year>2022</year>
</pub-date>
<pub-date pub-type="collection">
<year>2022</year>
</pub-date>
<volume>10</volume>
<elocation-id>864527</elocation-id>
<history>
<date date-type="received">
<day>28</day>
<month>01</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>12</day>
<month>04</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2022 Dong, Li, Yang, Wang, Zhao and Huo.</copyright-statement>
<copyright-year>2022</copyright-year>
<copyright-holder>Dong, Li, Yang, Wang, Zhao and Huo</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>Refracturing is an important method to improve the productivity of low-yield and low-efficiency horizontal wells. Due to the influence of geological factors and primary fracturing, it is of great importance to select suitable candidate wells and determine reasonable refracturing timing. Considering that production and pressure are the first-hand data available in the fields, they can be used to evaluate the dynamic drainage volume of fractured horizontal wells and to predict refracturing timing combined with fracturing parameters. In this study, the stimulation potential of refracturing is evaluated by introducing the dynamic drainage volume, and the refracturing timing discrimination model is established. Based on this model, the refracturing timing charts of different initial fracturing scales are drawn. In addition, the reservoir numerical simulation method is used to optimize the method of horizontal well refracturing. The results show that it is advisable to repeat fracturing when the change rate of dynamic drainage volume growth rate is less than 10%. The refracturing methods of restoring the conductivity of fractures or sealing the original fractures and refracturing the new fractures have good development effects. The research is of great significance to the design of horizontal well refracturing.</p>
</abstract>
<kwd-group>
<kwd>horizontal well</kwd>
<kwd>refracturing</kwd>
<kwd>dynamic drainage volume</kwd>
<kwd>refracturing timing</kwd>
<kwd>increasing production potential</kwd>
</kwd-group>
</article-meta>
</front>
<body>
<sec id="s1">
<title>Introduction</title>
<p>Horizontal well refracturing is one of the major technologies in unconventional oil and gas development, which can achieve greater utilization of single well reserves and further enhance recovery (<xref ref-type="bibr" rid="B4">Brown and Economides, 1992</xref>; <xref ref-type="bibr" rid="B6">Dou and Zhou, 2003</xref>). Due to the failure of fracture initiation, sand filling failure and other problems in the production of primary fractured horizontal Wells, fractures may lose their conductivity. Therefore, refracturing is required to restore or improve the production capacity of horizontal wells (<xref ref-type="bibr" rid="B10">Giger, 1985</xref>; <xref ref-type="bibr" rid="B19">Vincent, 2011</xref>; <xref ref-type="bibr" rid="B17">Ren et al., 2015</xref>; <xref ref-type="bibr" rid="B20">Wan et al., 2019</xref>). At present, a large number of wells in the early development stage have been refractured (<xref ref-type="bibr" rid="B18">Siebrits et al., 2000</xref>; <xref ref-type="bibr" rid="B7">French et al., 2014</xref>; <xref ref-type="bibr" rid="B8">Fu et al., 2017</xref>; <xref ref-type="bibr" rid="B2">Athavale et al., 2019</xref>; <xref ref-type="bibr" rid="B9">Garza et al., 2019</xref>), and considerable results have been achieved. The research on related technology types and operation schemes have gradually attracted the attention of scholars and experts around the wold. A large number of studies have been carried out at home and abroad on repeated fracturing well selection methods, production increase potential, and production increase measures (<xref ref-type="bibr" rid="B24">Yang et al., 2016</xref>; <xref ref-type="bibr" rid="B1">Aniemena and Kohshour, 2017</xref>; <xref ref-type="bibr" rid="B23">Xia et al., 2017</xref>; <xref ref-type="bibr" rid="B21">Wang, 2018</xref>; <xref ref-type="bibr" rid="B11">Ib&#xe1;&#xf1;ez et al., 2020</xref>; <xref ref-type="bibr" rid="B13">Kong et al., 2020</xref>; <xref ref-type="bibr" rid="B15">Lu et al., 2020</xref>; <xref ref-type="bibr" rid="B12">Kong et al., 2021</xref>; <xref ref-type="bibr" rid="B27">Zhang and Sun, 2021</xref>), However, relatively few studies have been conducted on the refracturing timing and the modification method of refracturing. At the same time, many scholars have carried out research on reservoir capacity characterization and dynamic prediction (<xref ref-type="bibr" rid="B14">Lee et al., 2003</xref>; <xref ref-type="bibr" rid="B22">Williams et al., 2004</xref>; <xref ref-type="bibr" rid="B3">Blasingame and Rushing, 2005</xref>; <xref ref-type="bibr" rid="B5">Cheng et al., 2009</xref>), There are mainly four methods:type-curve analysis,flow-regime analysis, numerical simulation and empir-ical methods. Despite the diversity of models, they always require knowledge of the reservoir that may not be known, originally, such as the shape of the drainage volume.Moghanloo and Yuan (<xref ref-type="bibr" rid="B16">Moghanloo et al., 2015</xref>) put forward the concept of dynamic drainage volume (DDV) for the first time. Although the model is not suitable for any reservoir, such as a fractured reservoir, but the data required by the model is easy to obtain. DDV refers to a continuously expanding oil drainage volume formed when the pressure fluctuation begins to spread in the whole reservoir when the production well is produced. With the progress of production, the boundary of DDV expands outward. Using DDV to evaluate the productivity of fractured horizontal wells does not need any prior knowledge of complex reservoir characteristics, and considering the advantages of dynamic change process. Therefore, the results are more in line with reality (<xref ref-type="bibr" rid="B26">Yuan et al., 2016</xref>; <xref ref-type="bibr" rid="B28">Zheng et al., 2017</xref>).</p>
<p>Based on the material balance theory, we deduce the calculation formulas of the DDV of horizontal wells. Using the diffusion law of DDV, the stimulation potential of wells is calculated, and the relationship between stimulation potential and refracturing timing is established. Moreover, the charts of refracturing timing with different initial pressure scales are drawn. Finally, the reconstruction mode of horizontal well refracturing is optimized by numerical simulation.</p>
</sec>
<sec id="s2">
<title>DDV Model of Depleted Horizontal Wells</title>
<p>When the horizontal wells are exhausted in the initial stage of production, it depends on the natural energy of the reservoir for production. The oil drainage volume in the reservoir reconstruction area directly reflects the productivity of the horizontal wells. Therefore, based on the material balance equation, the equations of DDV is established, the DDV is evaluated by using the pressure and production data, and the relationship between the refracturing timing and scale is established.</p>
<sec id="s2-1">
<title>Material Balance Equations</title>
<p>Assuming that the oil-water two-phase flow in the reservoir is quasi-steady state with constant temperature, and the fluid and rock in the reservoir are slightly compressible, the material balance equation is<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:mtext>&#xa0;&#xa0;&#xa0;&#xa0;</mml:mtext>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>where <inline-formula id="inf1">
<mml:math id="m2">
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the water production rate, m<sup>3</sup>/day; <inline-formula id="inf2">
<mml:math id="m3">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is water saturation; <inline-formula id="inf3">
<mml:math id="m4">
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the pore volume of oil drainage, m<sup>3</sup>; <inline-formula id="inf4">
<mml:math id="m5">
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the volume coefficient of formation water; <inline-formula id="inf5">
<mml:math id="m6">
<mml:mi>t</mml:mi>
</mml:math>
</inline-formula> is the time, day.<disp-formula id="e2">
<mml:math id="m7">
<mml:mrow>
<mml:mtext>&#xa0;&#xa0;&#xa0;&#xa0;&#xa0;&#xa0;&#xa0;</mml:mtext>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>o</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>o</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>where <inline-formula id="inf6">
<mml:math id="m8">
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the oil production rate, m<sup>3</sup>/day; <inline-formula id="inf7">
<mml:math id="m9">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>o</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is oil saturation; <inline-formula id="inf8">
<mml:math id="m10">
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>o</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the volume coefficient of formation oil. Combined with <xref ref-type="disp-formula" rid="e1">Eqs 1</xref>, <xref ref-type="disp-formula" rid="e2">2</xref>, we can get<disp-formula id="e3">
<mml:math id="m11">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>/</mml:mo>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>P</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>P</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
<disp-formula id="e4">
<mml:math id="m12">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>O</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>o</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>o</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>o</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>o</mml:mi>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>/</mml:mo>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>o</mml:mi>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>P</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>P</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
</p>
<p>
<xref ref-type="disp-formula" rid="e3">Eq. 3</xref> &#xd7; <inline-formula id="inf9">
<mml:math id="m13">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>W</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>, <xref ref-type="disp-formula" rid="e4">Eq. 4</xref> &#xd7; <inline-formula id="inf10">
<mml:math id="m14">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>o</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>, we can get<disp-formula id="e5">
<mml:math id="m15">
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>W</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>o</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>W</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>O</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>W</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>O</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>P</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>P</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>W</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>W</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>O</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>O</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>P</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>where: <inline-formula id="inf11">
<mml:math id="m16">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>W</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>O</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>; <inline-formula id="inf12">
<mml:math id="m17">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>O</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the isothermal compressibility coefficient of formation oil, MPa<sup>&#x2212;1</sup>; <inline-formula id="inf13">
<mml:math id="m18">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>W</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the isothermal compressibility coefficient of formation water, MPa<sup>&#x2212;1</sup>; <inline-formula id="inf14">
<mml:math id="m19">
<mml:mi>P</mml:mi>
</mml:math>
</inline-formula> is reservoir pressure, MPa. The comprehensive compressibility coefficient of the reservoir can be expressed as<disp-formula id="e6">
<mml:math id="m20">
<mml:mrow>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>P</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>o</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>o</mml:mi>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>o</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>P</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>o</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>P</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>
</p>
<p>
<xref ref-type="disp-formula" rid="e5">Eq. 5</xref> can be simplified as<disp-formula id="e7">
<mml:math id="m21">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>P</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>o</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
</p>
<p>According to <xref ref-type="disp-formula" rid="e7">Eq. 7</xref>
<disp-formula id="e8">
<mml:math id="m22">
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>P</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>o</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>
</p>
<p>Integrate the <xref ref-type="disp-formula" rid="e8">Eq. 8</xref> to obtain<disp-formula id="e9">
<mml:math id="m23">
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x222b;</mml:mo>
</mml:mstyle>
<mml:mn>0</mml:mn>
<mml:mi>t</mml:mi>
</mml:munderover>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mi>d</mml:mi>
<mml:mi>P</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x222b;</mml:mo>
</mml:mstyle>
<mml:mn>0</mml:mn>
<mml:mi>t</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>o</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>
<disp-formula id="e10">
<mml:math id="m24">
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x222b;</mml:mo>
</mml:mstyle>
<mml:mn>0</mml:mn>
<mml:mi>t</mml:mi>
</mml:munderover>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mi>d</mml:mi>
<mml:mi>P</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="true">&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mrow>
<mml:mi>W</mml:mi>
<mml:mi>P</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>O</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="true">&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mrow>
<mml:mi>O</mml:mi>
<mml:mi>P</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>where <inline-formula id="inf15">
<mml:math id="m25">
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="true">&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula> is the horizontal average volume coefficient of the formation; <inline-formula id="inf16">
<mml:math id="m26">
<mml:mrow>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mrow>
<mml:mi>W</mml:mi>
<mml:mi>P</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is cumulative water yield, m<sup>3</sup>; <inline-formula id="inf17">
<mml:math id="m27">
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>O</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="true">&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula> the average volume coefficient of formation oil; <inline-formula id="inf18">
<mml:math id="m28">
<mml:mrow>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mrow>
<mml:mi>O</mml:mi>
<mml:mi>P</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the cumulative oil production, m<sup>3</sup>.</p>
</sec>
<sec id="s2-2">
<title>DDV calculation of Oil</title>
<p>With the continuous exploitation of oil and gas reservoir, the reservoir pressure continues to spread outward, and the boundary of DDV continues to expand. This problem can be regarded as a dynamic boundary problem. It is considered that under the limitation of pressure range, the oil drain volume changes continuously and finally tends to the extreme value. The DDV of oil is defined as<disp-formula id="e11">
<mml:math id="m29">
<mml:mrow>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mrow>
<mml:msup>
<mml:mi>P</mml:mi>
<mml:mi>&#x221e;</mml:mi>
</mml:msup>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>P</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>P</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>P</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mover accent="true">
<mml:mo>&#x2192;</mml:mo>
<mml:mrow>
<mml:mi>y</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mover>
</mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>P</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mover accent="true">
<mml:mo>&#x2192;</mml:mo>
<mml:mrow>
<mml:mi>y</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mover>
</mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mrow>
<mml:msup>
<mml:mi>P</mml:mi>
<mml:mi>&#x221e;</mml:mi>
</mml:msup>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>Where <inline-formula id="inf19">
<mml:math id="m30">
<mml:mi>&#x3bb;</mml:mi>
</mml:math>
</inline-formula> is the rate at which the DDV expands; <inline-formula id="inf20">
<mml:math id="m31">
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> can be obtained by integrating <inline-formula id="inf21">
<mml:math id="m32">
<mml:mi>&#x3bb;</mml:mi>
</mml:math>
</inline-formula>; <inline-formula id="inf22">
<mml:math id="m33">
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mrow>
<mml:msup>
<mml:mi>P</mml:mi>
<mml:mi>&#x221e;</mml:mi>
</mml:msup>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the final pore volume of oil drainage, m<sup>3</sup>; <inline-formula id="inf23">
<mml:math id="m34">
<mml:mrow>
<mml:mover accent="true">
<mml:mi>P</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula> is the average formation pressure, MPa; <inline-formula id="inf24">
<mml:math id="m35">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the abandoned formation pressure, MPa. By substituting <xref ref-type="disp-formula" rid="e10">Eq. 10</xref> into <xref ref-type="disp-formula" rid="e11">Eq. 11</xref>, the DDV of oil can be obtained<disp-formula id="e12">
<mml:math id="m36">
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x222b;</mml:mo>
</mml:mstyle>
<mml:mn>0</mml:mn>
<mml:mi>t</mml:mi>
</mml:munderover>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mrow>
<mml:msup>
<mml:mi>P</mml:mi>
<mml:mi>&#x221e;</mml:mi>
</mml:msup>
</mml:mrow>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>P</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>P</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mi>d</mml:mi>
<mml:mi>P</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="true">&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mrow>
<mml:mi>W</mml:mi>
<mml:mi>P</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>O</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="true">&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mrow>
<mml:mi>O</mml:mi>
<mml:mi>P</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>
</p>
</sec>
<sec id="s2-3">
<title>Calculation Process</title>
<p>Firstly, the average pressure of the studied reservoir is estimated to obtain the values of <inline-formula id="inf25">
<mml:math id="m37">
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf26">
<mml:math id="m38">
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>O</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf27">
<mml:math id="m39">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>O</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf28">
<mml:math id="m40">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf29">
<mml:math id="m41">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>f</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>. Then the average physical parameter values between the original reservoir pressure and the average pressure are calculated by using the empirical formula and substituted into <xref ref-type="disp-formula" rid="e10">Eq. 10</xref>. Estimating <inline-formula id="inf30">
<mml:math id="m42">
<mml:mi>&#x3bb;</mml:mi>
</mml:math>
</inline-formula> values, use an iterative process based on varying l to find the drainage pore volume and average reservoir pressure at any time in the production history of the well, when the error range is met, the DDV of oil and the corresponding average pressure value are output (<xref ref-type="fig" rid="F1">Figure 1</xref>).</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Dynamic drainage volume calculation flow chart.</p>
</caption>
<graphic xlink:href="fenrg-10-864527-g001.tif"/>
</fig>
</sec>
</sec>
<sec id="s3">
<title>Determination Method of Refracturing Timing Based on DDV</title>
<p>Because it is difficult to obtain a large number of production data of fractured horizontal wells, this paper set different fracture parameters and established the numerical simulation model of multi-stage fractured horizontal wells under depletion development mode. Based on the reservoir physical parameters and production data under different initial pressure modes, the variation law of DDV is obtained by using the numerical simulation model. The relationship between initial pressure scale and refracturing timing is also analyzed. The numerical simulation model of a typical ultra-low permeability reservoir in Changqing Oilfield is established: the number of grids is 81 &#xd7; 131&#xd7;3. The grid unit size is 20&#xa0;m, the average porosity is 11%, and the average permeability is 0.5 &#xd7; 10<sup>&#x2013;3</sup>&#xa0;&#x3bc;m<sup>2</sup>. The fluid physical parameters are listed in <xref ref-type="table" rid="T1">Table 1</xref>, the relative permeability curve is shown in <xref ref-type="fig" rid="F2">Figure 2</xref> and the initial fracturing parameters are listed in <xref ref-type="table" rid="T2">Table 2</xref>.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Parameters of depleted development liquid model.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Parameter</th>
<th align="left">Value</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Initial Reservior Pressure/MPa</td>
<td align="char" char=".">14.71</td>
</tr>
<tr>
<td align="left">Reservior Oil Viscosity/(mPa&#xb7;s)</td>
<td align="char" char=".">1.27</td>
</tr>
<tr>
<td align="left">Surface Oil Density/(g/cm<sup>3</sup>)</td>
<td align="char" char=".">0.78</td>
</tr>
<tr>
<td align="left">Oil Volume Factor/(m<sup>3</sup>/m<sup>3</sup>)</td>
<td align="char" char=".">1.24</td>
</tr>
<tr>
<td align="left">Original Gas-Oil ratio/(m<sup>3</sup>/m<sup>3</sup>)</td>
<td align="char" char=".">103.30</td>
</tr>
<tr>
<td align="left">Oil Compressibility Factor/MPa<sup>&#x2212;1</sup>
</td>
<td align="char" char=".">0.001493</td>
</tr>
<tr>
<td align="left">Rock Compressibility Factor/MPa<sup>&#x2212;1</sup>
</td>
<td align="char" char=".">0.000435</td>
</tr>
<tr>
<td align="left">Water Compressibility Factor/MPa<sup>&#x2212;1</sup>
</td>
<td align="char" char=".">0.000501</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Relative permeability curve of depleted development model.</p>
</caption>
<graphic xlink:href="fenrg-10-864527-g002.tif"/>
</fig>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Initial fracturing parameters of depleted development well.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Fracture parameter name</th>
<th colspan="3" align="center">Value</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Fracture Spacing/m</td>
<td align="char" char=".">40</td>
<td align="char" char=".">80</td>
<td align="char" char=".">160</td>
</tr>
<tr>
<td align="left">Fracture Half-length/m</td>
<td align="char" char=".">100</td>
<td align="char" char=".">240</td>
<td align="char" char=".">320</td>
</tr>
<tr>
<td align="left">Fracture Conductivity/(10<sup>&#x2212;3</sup>&#xa0;&#x3bc;m<sup>2</sup> m)</td>
<td align="char" char=".">20</td>
<td align="char" char=".">120</td>
<td align="char" char=".">240</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>There are 27 fracturing parameter combinations in <xref ref-type="table" rid="T2">Table 2</xref>, which are shown in <xref ref-type="table" rid="T3">Table 3</xref>.</p>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>Initial fracturing scheme of depleted development.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th rowspan="2" align="left">Scheme</th>
<th align="center">Fracture</th>
<th align="center">Fracture</th>
<th align="center">Fracture conductivit</th>
</tr>
<tr>
<th align="center">Spacing/m</th>
<th align="center">Half-Length/m</th>
<th align="center">/(10<sup>&#x2212;3</sup>&#xa0;&#x3bc;m<sup>2</sup> m)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">1</td>
<td align="char" char=".">40</td>
<td align="char" char=".">100</td>
<td align="char" char=".">20</td>
</tr>
<tr>
<td align="left">2</td>
<td align="char" char=".">40</td>
<td align="char" char=".">100</td>
<td align="char" char=".">120</td>
</tr>
<tr>
<td align="left">3</td>
<td align="char" char=".">40</td>
<td align="char" char=".">100</td>
<td align="char" char=".">240</td>
</tr>
<tr>
<td align="left">4</td>
<td align="char" char=".">40</td>
<td align="char" char=".">240</td>
<td align="char" char=".">20</td>
</tr>
<tr>
<td align="left">5</td>
<td align="char" char=".">40</td>
<td align="char" char=".">240</td>
<td align="char" char=".">120</td>
</tr>
<tr>
<td align="left">6</td>
<td align="char" char=".">40</td>
<td align="char" char=".">240</td>
<td align="char" char=".">240</td>
</tr>
<tr>
<td align="left">7</td>
<td align="char" char=".">40</td>
<td align="char" char=".">320</td>
<td align="char" char=".">20</td>
</tr>
<tr>
<td align="left">8</td>
<td align="char" char=".">40</td>
<td align="char" char=".">320</td>
<td align="char" char=".">120</td>
</tr>
<tr>
<td align="left">9</td>
<td align="char" char=".">40</td>
<td align="char" char=".">320</td>
<td align="char" char=".">240</td>
</tr>
<tr>
<td align="left">10</td>
<td align="char" char=".">80</td>
<td align="char" char=".">100</td>
<td align="char" char=".">20</td>
</tr>
<tr>
<td align="left">11</td>
<td align="char" char=".">80</td>
<td align="char" char=".">100</td>
<td align="char" char=".">120</td>
</tr>
<tr>
<td align="left">12</td>
<td align="char" char=".">80</td>
<td align="char" char=".">100</td>
<td align="char" char=".">240</td>
</tr>
<tr>
<td align="left">13</td>
<td align="char" char=".">80</td>
<td align="char" char=".">240</td>
<td align="char" char=".">20</td>
</tr>
<tr>
<td align="left">14</td>
<td align="char" char=".">80</td>
<td align="char" char=".">240</td>
<td align="char" char=".">120</td>
</tr>
<tr>
<td align="left">15</td>
<td align="char" char=".">80</td>
<td align="char" char=".">240</td>
<td align="char" char=".">240</td>
</tr>
<tr>
<td align="left">16</td>
<td align="char" char=".">80</td>
<td align="char" char=".">320</td>
<td align="char" char=".">20</td>
</tr>
<tr>
<td align="left">17</td>
<td align="char" char=".">80</td>
<td align="char" char=".">320</td>
<td align="char" char=".">120</td>
</tr>
<tr>
<td align="left">18</td>
<td align="char" char=".">80</td>
<td align="char" char=".">320</td>
<td align="char" char=".">240</td>
</tr>
<tr>
<td align="left">19</td>
<td align="char" char=".">160</td>
<td align="char" char=".">100</td>
<td align="char" char=".">20</td>
</tr>
<tr>
<td align="left">20</td>
<td align="char" char=".">160</td>
<td align="char" char=".">100</td>
<td align="char" char=".">120</td>
</tr>
<tr>
<td align="left">21</td>
<td align="char" char=".">160</td>
<td align="char" char=".">100</td>
<td align="char" char=".">240</td>
</tr>
<tr>
<td align="left">22</td>
<td align="char" char=".">160</td>
<td align="char" char=".">240</td>
<td align="char" char=".">20</td>
</tr>
<tr>
<td align="left">23</td>
<td align="char" char=".">160</td>
<td align="char" char=".">240</td>
<td align="char" char=".">120</td>
</tr>
<tr>
<td align="left">24</td>
<td align="char" char=".">160</td>
<td align="char" char=".">240</td>
<td align="char" char=".">240</td>
</tr>
<tr>
<td align="left">25</td>
<td align="char" char=".">160</td>
<td align="char" char=".">320</td>
<td align="char" char=".">20</td>
</tr>
<tr>
<td align="left">26</td>
<td align="char" char=".">160</td>
<td align="char" char=".">320</td>
<td align="char" char=".">120</td>
</tr>
<tr>
<td align="left">27</td>
<td align="char" char=".">160</td>
<td align="char" char=".">320</td>
<td align="char" char=".">240</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The cumulative oil production and cumulative water production under different development modes are calculated and brought into the material balance equation. The production is normalized by considering the physical properties of the reservoir, and finally the variation value of the DDV with time is obtained. <xref ref-type="fig" rid="F3">Figure 3</xref> shows the variation of DDV with time when the interval of fracturing section is 40&#xa0;m, the half length of fracture is 100&#xa0;m, and the conductivity of initial fracturing fracture is 20 &#xd7; 10<sup>&#x2013;3</sup>&#xa0;&#x3bc;m<sup>2</sup> m. <xref ref-type="fig" rid="F4">Figure 4</xref> shows the change rate of DDV growth rate. The trend of drainage pore volume shows an asymptotic features, which is consistent with the previous studies (<xref ref-type="bibr" rid="B25">Yin et al., 2015</xref>).Due to the limitation of reservoir pressure, the DDV finally tends to a limit value.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Dynamic drainage volume variation over time.</p>
</caption>
<graphic xlink:href="fenrg-10-864527-g003.tif"/>
</fig>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Dynamic drainage volume growth rate.</p>
</caption>
<graphic xlink:href="fenrg-10-864527-g004.tif"/>
</fig>
<p>See <xref ref-type="fig" rid="F5">Figure 5</xref> for the flow chart of determining the refracturing timing by using the DDV method. According to the cumulative water production and cumulative oil production under different initial fracture distribution, the DDV change law curve with time is obtained by using <xref ref-type="disp-formula" rid="e12">Eq. 12</xref> respectively. The curve slope is then used to calculate the growth rate and change rate of the growth rate of DDV. When the change rate of the growth rate of DDV is less than 10%, the dynamic drainage boundary propagates slowly, the reservoir is in the development state of low production and low efficiency, and the output of the initially fractured oil well has not reached the expected effect. However, at this time, a large amount of residual oil still exists in the reservoir and between fractures. The horizontal well needs to be refractured to restore the oil well productivity. Through this method, we can obtain the refracturing timing of different initial fracturing scales under 27 groups of schemes, as shown in <xref ref-type="fig" rid="F6">Figure 6</xref>.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Flow chart of determining refracturing timing by dynamic drainage volume method.</p>
</caption>
<graphic xlink:href="fenrg-10-864527-g005.tif"/>
</fig>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Refracturing timing sequence under different initial fracturing modes.</p>
</caption>
<graphic xlink:href="fenrg-10-864527-g006.tif"/>
</fig>
<p>The refracturing timing under different initial pressure scales is classified, and the chart shown in <xref ref-type="fig" rid="F7">Figure 7</xref>. The horizontal axis represents the interval and the vertical axis represents the half length of the fracture. There are 9 grids in the chart and each grid is divided into 3 layers according to the fracture conductivity. The fracture conductivity from bottom to top is 20 &#xd7; 10<sup>&#x2013;3</sup>&#xa0;<italic>&#x3bc;</italic>m<sup>2</sup>&#xb7;m, 120 &#xd7; 10<sup>&#x2013;3</sup>&#xa0;<italic>&#x3bc;</italic>m<sup>2</sup>&#xb7;m and 240 &#xd7; 10<sup>&#x2013;3</sup>&#xa0;<italic>&#x3bc;</italic>m<sup>2</sup>&#xb7;m, the chart is divided into 27 grids, representing 27 different schemes. We divided the grids into four categories: Class I is horizontal wells with poor initial development effect, which are production wells with low conductivity and large interval. It is recommended to carry out repeated fracturing in the first 20&#xa0;months of the initial stage of production; Class II and III are horizontal wells with good development effect at the initial stage of production and low production after a period of time due to insufficient fracture length. It is recommended to extend the fracture length by half; Class IV is horizontal wells with good fracture length and conductivity and less influence of interval spacing. It is recommended to carry out refracturing finally. In conclusion, it is suggested that the sequence of refracturing for horizontal wells with different initial pressure scales should be followed: first, refracturing for horizontal wells of class I to restore the conductivity of fractures and add new fractures; Secondly, horizontal wells of class II and III are refractured to prolong the fracture length and restore the fracture conductivity; Finally, horizontal wells of class IV are refractured.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Chart of refracturing timing at different primary fracturing scales.</p>
</caption>
<graphic xlink:href="fenrg-10-864527-g007.tif"/>
</fig>
</sec>
<sec id="s4">
<title>Study on Optimization of Horizontal Refracturing Reconstruction Mode</title>
<p>In the process of reservoir reconstruction, two types of refracturing are usually used: original fracture refracturing and new fracture refracturing. The main ways of refracturing the original fracture are: prolonging the old fracture, turning the old fracture and restoring the conductivity of the original fracture. There are mainly two ways to refracture new fractures: retaining the original fracture and sealing the original fracture. Choosing different refracturing reconstruction methods will lead to great changes in the seepage law and stimulation mechanism of horizontal wells after refracturing.</p>
<sec id="s4-1">
<title>Reconstruction of Fractured Horizontal Well With original Fracture Refracturing</title>
<sec id="s4-1-1">
<title>Restore the Original Fracture Conductivity</title>
<p>By restoring the conductivity of the failed fracture, the horizontal well can obtain higher productivity. Due to the crushing and migration of proppant, the fracture conductivity changes from time to time and decreases rapidly in the initial production stage. In view of this situation, it is recommended to use high-strength proppant during refracturing, as shown in <xref ref-type="fig" rid="F8">Figure 8</xref>. The fracture conductivity before refracturing is 25 &#xd7; 10<sup>&#x2013;3</sup>&#xa0;&#x3bc;m<sup>2</sup>&#xa0;m, the recovery conductivity after refracturing is 300 &#xd7; 10<sup>&#x2013;3</sup>&#xa0;&#x3bc;m<sup>2</sup>&#xa0;m. Through the established numerical simulation model, the interval between segments is 80&#xa0;m and the half length of fractures is 120&#xa0;m.</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Variation diagram of conductivity of recovery fracture.</p>
</caption>
<graphic xlink:href="fenrg-10-864527-g008.tif"/>
</fig>
<p>In the simulation process, the horizontal well is normally produced for three years under the original fracture conditions, and refractureing production is carried out for ten years. <xref ref-type="fig" rid="F9">Figure 9</xref> shows the comparison of pressure and remaining oil distribution before and after restoring conductivity of failed fractures. By comparing and analyzing the simulation results, it can be seen that in the initial stage of development, the fracture conductivity is reduced due to fracture failure, and the pressure propagation near the horizontal well is difficult. The distribution of reservoir pressure directly affects the failure development effect of the horizontal well. After refracturing, the fracture conductivity is restored and the &#x201c;pollution&#x201d; in the fracture is relieved. The pressure can diffuse to the far well zone along the fracture, thus improve the effective production degree of low-permeability reservoir in the far well zone, and expand the DDV of horizontal wells.</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>Pressure distribution before and after recovery of fracture conductivity. <bold>(A)</bold> Pressure distribution of fracturing failure horizontal well before refracturing. <bold>(B)</bold> Pressure distribution after restoring conductivity. <bold>(C)</bold> Remaining oil distribution after 10&#xa0;years without refracturing. <bold>(D)</bold> Remaining oil distribution after 10&#xa0;years after restoring conductivity.</p>
</caption>
<graphic xlink:href="fenrg-10-864527-g009.tif"/>
</fig>
</sec>
<sec id="s4-1-2">
<title>Increase the Original Fracture Length</title>
<p>Based on the basic model of numerical simulation, See <xref ref-type="table" rid="T4">Table 4</xref> for parameter values of cracks. In the simulation process, similarily, the horizontal well is normally produced for three years under the original fracture conditions, and refractureing production is carried out for ten years. <xref ref-type="fig" rid="F10">Figure 10</xref> shows the comparison of pressure and residual oil distribution before and after increasing the original fracture length. From the simulation results, it can be seen that the depleted development production depends on the reconstrution of the reservoir. Extending the original fracture length expands the fracturing scale of the horizontal well, increases the effective reconstruction volume of the reservoir, expands the spread range of the pressure wave, and significantly increases the oil drainage volume of the horizontal well. Therefore, the far well reservoir fluid can flow into the wellbore and improve the production of the horizontal well.</p>
<table-wrap id="T4" position="float">
<label>TABLE 4</label>
<caption>
<p>Increase the original fracture length Fracture parameter table.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Fracture parameter name</th>
<th align="center">Value</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Fracture Spacing/m</td>
<td align="char" char=".">80</td>
</tr>
<tr>
<td align="left">Primary Fracture Half-length/m</td>
<td align="char" char=".">120</td>
</tr>
<tr>
<td align="left">Refracturing Fracture Half-length/m</td>
<td align="char" char=".">200</td>
</tr>
<tr>
<td align="left">Fracture Conductivity/(10<sup>&#x2212;3</sup>&#xa0;&#x3bc;m<sup>2</sup> m)</td>
<td align="char" char=".">100</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>Pressure distribution before and after increasing original fracture length and remaining oil distribution in 10&#xa0;years of production. <bold>(A)</bold> Pressure distribution before refracturing <bold>(B)</bold> Pressure distribution after increasing the orginally fracture length. <bold>(C)</bold> Remaining oil distribution in 10 years without refracturing. <bold>(D)</bold> Remaining oil distribution after ten years after increasing the orginally fracture length.</p>
</caption>
<graphic xlink:href="fenrg-10-864527-g010.tif"/>
</fig>
</sec>
<sec id="s4-1-3">
<title>Reorientation of Original Fractures</title>
<p>With the process of depletion development, the <italic>in-situ</italic> stress field near the hydraulic fracture of initial fracturing will change, and the tensile stress parallel to the original fracture will be less than the vertical tensile stress. When the stress difference caused by pore pressure is greater than the initial maximum and minimum horizontal stress difference, the maximum stress direction of the original stress field will change to the minimum stress direction. After that, the refracturing measures are implemented to add new fractures. Under the action of the stress field, the fractures will turn in the stress shadow area near the horizontal well, and then the fracture direction will be parallel to the original fractures and continue to expand to form new fractures, as shown in <xref ref-type="fig" rid="F11">Figure 11</xref>.</p>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>Diagram of fracture diversion in horizontal well refracturing.</p>
</caption>
<graphic xlink:href="fenrg-10-864527-g011.tif"/>
</fig>
<p>Using the reservoir numerical simulation method, ignoring the actual complex fracture network and the complex variation law of <italic>in-situ</italic> stress deflection area. The refracturing measures are implemented after three years of production, and the original fracture deflects 45&#xb0; near the well zone. <xref ref-type="fig" rid="F12">Figure 12</xref> shows the residual oil and pressure distribution before and after fracture reorientation. After three years of production, the <italic>in-situ</italic> stress field of the reservoir changes, and the fracture continues to expand along the direction of the maximum principal stress after turning. Therefore, refracturing effectively expands the propagation range of pressure wave. At the same time, the fracture reorientation increases the contact area between the fracture and the reservoir. The steering fracture can extend to the far well area, which can expand the spread range of the fracture and significantly increase the oil drainage volume of the horizontal well.</p>
<fig id="F12" position="float">
<label>FIGURE 12</label>
<caption>
<p>Pressure distribution and remaining oil distribution before and after fracture turning. <bold>(A)</bold> Pressure distribution before fracture reorientation. <bold>(B)</bold> Pressure distribution after fracture reorientation. <bold>(C)</bold> Remaining oil distribution before fracture reorientation. <bold>(D)</bold> Remaining oil distribution after fracture reorientation.</p>
</caption>
<graphic xlink:href="fenrg-10-864527-g012.tif"/>
</fig>
<p>According to the numerical simulation test scheme designed above, the model adopts constant bottom hole flow pressure production. <xref ref-type="fig" rid="F13">Figure 13</xref> shows the change results of cumulative oil increase in ten years of refracturing production in the test scheme. By analyzing the variation law of cumulative oil increment in depleted development horizontal wells, it can be concluded that the influence of refracturing ways on their effect is as follows. The growth of fracture half length and fracture reorientation effectively increase the scope of reservoir reconstruction, so as to increase the oil drainage volume of horizontal wells and obtain better refracturing effect. Restoring the conductivity of the original fracture dredges the seepage channel of crude oil in the reservoir. The crude oil in the fracture is easier to flow to the wellbore and increases the oil production. The fracture conductivity has a greater impact on the re fracturing effect than the half length and steering of the fracture.</p>
<fig id="F13" position="float">
<label>FIGURE 13</label>
<caption>
<p>Change of cumulative oil increment in ten years of test scheme.</p>
</caption>
<graphic xlink:href="fenrg-10-864527-g013.tif"/>
</fig>
</sec>
</sec>
<sec id="s4-2">
<title>Reconstruction of Fractured Horizontal Well With Refracturing New Fractures</title>
<p>Using the numerical simulation method, the initial fracture half length is set as 120&#xa0;m and the fracture conductivity is 100 &#xd7; 10<sup>&#x2013;3</sup>&#xa0;&#x3bc;m&#xa0;m, the interval before refracturing is 160&#xa0;m, and that after refracturing new fractures is 80&#xa0;m. The horizontal well normally produces for three years under the condition of original fractures, and then refracturing is implemented. <xref ref-type="fig" rid="F14">Figure 14</xref> shows the comparison of pressure distribution and remaining oil distribution of different refracturing methods. <xref ref-type="fig" rid="F15">Figure 15</xref> shows the results of ten-year cumulative oil increase under two refracturing modes.</p>
<fig id="F14" position="float">
<label>FIGURE 14</label>
<caption>
<p>Distribution diagram of pressure and remaining oil under different refracturing modes. <bold>(A)</bold> Pressure distribution refracturing before <bold>(B)</bold> Pressure distribution after remaining gold fracture and refracturing new fracture <bold>(C)</bold> Pressure distribution after plugging old fracture and refracturing new fracture <bold>(D)</bold> Remaining oil distribution before refracturing <bold>(E)</bold> Remaining oil distribution after remaining old fracture and refracturing new fracture <bold>(F)</bold> Remaining oil distribution after plugging old fracture and refracturing new fracture</p>
</caption>
<graphic xlink:href="fenrg-10-864527-g014.tif"/>
</fig>
<fig id="F15" position="float">
<label>FIGURE 15</label>
<caption>
<p>Variation diagram of ten-year cumulative oil increase under compound pressure under different transformation modes.</p>
</caption>
<graphic xlink:href="fenrg-10-864527-g015.tif"/>
</fig>
<p>By analyzing the variation law of seepage field and production, it is found that after refracturing, new fractures are added, fracture spacing is reduced, crude oil seepage resistance is reduced, pressure sweep range is expanded, and crude oil flow channel is increased. Therefore the crude oil production is improved. Moreover, sealing the original fracture and refracturing the new fracture changes the pressure distribution between the fractures, reduces the crude oil seepage resistance, and effectively initates the remaining oil around the original fracture. Compared with retaining the original fracture and refracturing the new fracture, it can obtain better development effect.</p>
</sec>
</sec>
<sec id="s5">
<title>Field Application</title>
<p>One horizontal well in a typical ultra-low permeability reservoir in Changqing Oilfield was put into operation in October 2012. After 6&#xa0;years of initial fracturing production, it was refractured, 19 old fractures were blocked, and 23 new fractures were refractured between the old fractures. After the measures, the daily oil production is 15.9&#xa0;t and the daily oil increase is 13.3&#xa0;t, as shown in <xref ref-type="fig" rid="F16">Figure 16</xref>. At present, it has been produced for nearly one year, with a daily oil production of 14.29&#xa0;t, a daily oil increase of 11.61&#xa0;t, a cumulative oil increase of 3442&#xa0;t and a cumulative oil production of 4407&#xa0;t.</p>
<fig id="F16" position="float">
<label>FIGURE 16</label>
<caption>
<p>Production variation law of a well before and after refracturing.</p>
</caption>
<graphic xlink:href="fenrg-10-864527-g016.tif"/>
</fig>
</sec>
<sec sec-type="conclusion" id="s6">
<title>Conclusion</title>
<p>In this paper, by introducing the calculation method and process of DDV, the discrimination model of refracturing timing is established. The refracturing timing chart under different initial pressure scale is formed, and the different reconstruction methods of refracturing are optimized and designed by numerical simulation method. The results show that.<list list-type="simple">
<list-item>
<p>1) With the depletion of the reservoir, the reservoir pressure continues to spread outward, and the boundary of the reservoir drainage volume continues to expand. However, the drainage volume finally tends to the extreme value under the limitation of the pressure range. The increase rate of DDV of horizontal wells shows a trend of first fast and then slow. When the change rate of growth rate is less than 10%, it can not meet the expected production requirements and needs to be refractured.</p>
</list-item>
<list-item>
<p>2) Considering the horizontal fracture distribution of horizontal wells, it is suggested that the sequence of refracturing of horizontal wells with different initial pressure scales is as follows: firstly, refracturing is carried out for horizontal wells with few fractures, poor fracture conductivity and poor initial development effect, so as to restore the fracture conductivity and infill new fractures; Secondly, the horizontal wells with insufficient fracture half length and poor development effect are refractured to prolong the fracture length and restore the fracture conductivity, so as to increase the reservoir reconstruction volume and obtain better development effect.</p>
</list-item>
<list-item>
<p>3) Through the implementation of refracturing reconstruction measures, the reconstruction effect of original fracture refracturing is as follows: restore the original fracture conductivity &#x3e; increase the original fracture length &#x3e; turn the original fracture. The development effect of plugging the original fractures and refracturing the new fractures is better than that of retaining the original fracturing and refracturing the new fracturing.</p>
</list-item>
</list>
</p>
</sec>
</body>
<back>
<sec id="s7">
<title>Data Availability Statement</title>
<p>The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s8">
<title>Author Contributions</title>
<p>QD: Conceptualization, Methodology, Formal Analysis, Writing&#x2014;Original Draft; JL: Data Curation, Writing&#x2014;Original Draft; LY: Visualization, Investigation; FW: Resources, Supervision; KZ: Visualization; FH:Writing&#x2014;Review and Editing.</p>
</sec>
<sec id="s9">
<title>Funding</title>
<p>This study was supported by the National major special project &#x201c;tackling key problems and field test of tight oil and gas multi-layer and multi-stage fracturing technology in Changqing Oilfield&#x201d; (2020CGCGZ004).</p>
</sec>
<sec sec-type="COI-statement" id="s10">
<title>Conflict of Interest</title>
<p>QD, JL, LY, and FW were employed by the company Oil &#x26;and Gas Technology Research Institute ChangQing Oilfield Company.</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="s11">
<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="confproc">
<person-group person-group-type="author">
<name>
<surname>Aniemena</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Kohshour</surname>
<given-names>I. O.</given-names>
</name>
</person-group> (<year>2017</year>). &#x201c;<article-title>A Fast Semi-analytical Method for Refracturing Candidate Selection and Performance Estimation of Shale Wells</article-title>,&#x201d; in <conf-name>Unconventional Resources Technology Conference</conf-name>, <conf-loc>Austin, Texas, USA</conf-loc>, <conf-date>24-26 July 2017</conf-date>. <comment>paper 2693452-MS presented at the</comment>. <pub-id pub-id-type="doi">10.15530/URTEC-2017-2693452</pub-id> </citation>
</ref>
<ref id="B2">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Athavale</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Onafowokan</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Kasdano</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Lakshminarayanan</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Sousa</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2019</year>). &#x201c;<article-title>Refracturing Technology Restores Profitability of Older Unconventional Wells in Southeast New Mexico</article-title>,&#x201d; in <conf-name>SPE Hydraulic Fracturing Technology Conference and Exhibition</conf-name> (<publisher-loc>Woodlands, TX</publisher-loc>: <publisher-name>OnePetro</publisher-name>). <pub-id pub-id-type="doi">10.2118/194372-MS</pub-id> </citation>
</ref>
<ref id="B3">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Blasingame</surname>
<given-names>T. A.</given-names>
</name>
<name>
<surname>Rushing</surname>
<given-names>J. A.</given-names>
</name>
</person-group> (<year>2005</year>). &#x201c;<article-title>A Production-Based Method for Direct Estimation of Gas in Place and Reserves</article-title>,&#x201d; in <conf-name>SPE Eastern Regional Meeting</conf-name>. <pub-id pub-id-type="doi">10.2118/98042-MS</pub-id> </citation>
</ref>
<ref id="B4">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Brown</surname>
<given-names>J. E.</given-names>
</name>
<name>
<surname>Economides</surname>
<given-names>M. J.</given-names>
</name>
</person-group> (<year>1992</year>). &#x201c;<article-title>An Analysis of Horizontally Fractured Horizontal Wells</article-title>,&#x201d; in <conf-name>SPE Rocky Mountain Regional Meeting</conf-name>, <conf-loc>Casper,Wyoming</conf-loc>, <conf-date>18-21 May 1992</conf-date>. <comment>paper 24322-MS presented at the</comment>. <pub-id pub-id-type="doi">10.2118/24322-MS</pub-id> </citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cheng</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Lee</surname>
<given-names>W. J.</given-names>
</name>
<name>
<surname>McVay</surname>
<given-names>D. A.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>A New Approach for Reliable Estimation of Hydraulic Fracture Properties Using Elliptical Flow Data in Tight Gas Wells</article-title>. <source>SPE Reserv. Eval. Eng.</source> <volume>12</volume> (<issue>02</issue>), <fpage>254</fpage>&#x2013;<lpage>262</lpage>. <pub-id pub-id-type="doi">10.2118/105767-PA</pub-id> </citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dou</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Remainined Oil Distribution in Complicated Faulted Blocks and its Accessory Potential Exploitation Measures</article-title>. <source>Petroleum Explor. Dev.</source> <volume>30</volume> (<issue>5</issue>), <fpage>90</fpage>&#x2013;<lpage>93</lpage>. <pub-id pub-id-type="doi">10.3321/j.issn:1000-0747.2003.05.027</pub-id> </citation>
</ref>
<ref id="B7">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>French</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Rodgerson</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Feik</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2014</year>). &#x201c;<article-title>Re-fracturing Horizontal Shale Wells: Case History of a Woodford Shale Pilot Project</article-title>,&#x201d; in <conf-name>SPE Hydraulic Fracturing Technology Conference</conf-name> (<publisher-loc>Woodlands, TX</publisher-loc>: <publisher-name>OnePetro</publisher-name>). <pub-id pub-id-type="doi">10.2118/168607-MS</pub-id> </citation>
</ref>
<ref id="B8">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Fu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Xiu</surname>
<given-names>N.</given-names>
</name>
</person-group> (<year>2017</year>). &#x201c;<article-title>Large-scale Experimental Investigation of Refracture Reorientation and Field Evaluation in Sand Reservoir</article-title>,&#x201d; in <conf-name>51st US Rock Mechanics/Geomechanics Symposium</conf-name> (<publisher-loc>San Francisco, CA</publisher-loc>: <publisher-name>OnePetro</publisher-name>). </citation>
</ref>
<ref id="B9">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Garza</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Baumbach</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Prosser</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Pettigrew</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Elvig</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2019</year>). &#x201c;<article-title>An Eagle Ford Case Study: Improving an Infill Well Completion through Optimized Refracturing Treatment of the Offset Parent Wells</article-title>,&#x201d; in <conf-name>SPE Hydraulic Fracturing Technology Conference and Exhibition</conf-name> (<publisher-loc>Bahrain</publisher-loc>: <publisher-name>OnePetro</publisher-name>). <pub-id pub-id-type="doi">10.2118/194374-MS</pub-id> </citation>
</ref>
<ref id="B10">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Giger</surname>
<given-names>F. M.</given-names>
</name>
</person-group> (<year>1985</year>). &#x201c;<article-title>Horizontal Wells Production Techniques in Heterogeneous Reservoirs</article-title>,&#x201d; in <conf-name>SPE Middle East Oil Technical Conference and Exhibition</conf-name>, <conf-date>11-14 March 1985</conf-date>. <comment>paper 13710-MS presented at the</comment>. <pub-id pub-id-type="doi">10.2118/13710-MS</pub-id> </citation>
</ref>
<ref id="B11">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Ib&#xe1;&#xf1;ez</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Bezerra de Melo</surname>
<given-names>R. C.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Dong</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Alvarellos</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Lakshmikantha</surname>
<given-names>M. R.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). &#x201c;<article-title>Improved Refracturing Evaluation Workflow: An Integrated Methodology Incorporating Advanced Geomechanical Modeling</article-title>,&#x201d; in <conf-name>SPE Canada Unconventional Resources Conference</conf-name> (<publisher-loc>Virtual</publisher-loc>: <publisher-name>OnePetro</publisher-name>). <pub-id pub-id-type="doi">10.2118/199997-MS</pub-id> </citation>
</ref>
<ref id="B12">
<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="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kong</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Lian</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Pore-scale Investigation of Immiscible Gas-Assisted Gravity Drainage</article-title>. <source>Phys. Fluids</source> <volume>32</volume> (<issue>12</issue>), <fpage>122004</fpage>. <pub-id pub-id-type="doi">10.1063/5.0033027</pub-id> </citation>
</ref>
<ref id="B14">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Lee</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Rollins</surname>
<given-names>J. B.</given-names>
</name>
<name>
<surname>Spivey</surname>
<given-names>J. P.</given-names>
</name>
</person-group> (<year>2003</year>). <source>Pressure Transient Testing (eBooks)[M]</source>. <publisher-loc>Houston, TX</publisher-loc>: <publisher-name>Society of Petroleum Engineers</publisher-name>, <fpage>9</fpage>. </citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lu</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Su</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhan</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Almrabat</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Modeling for Reorientation and Potential of Enhanced Oil Recovery in Refracturing</article-title>. <source>Adv. Geo-Energ. Res.</source> <volume>4</volume> (<issue>1</issue>), <fpage>20</fpage>&#x2013;<lpage>28</lpage>. <pub-id pub-id-type="doi">10.26804/ager.2020.01.03</pub-id> </citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Moghanloo</surname>
<given-names>R. G.</given-names>
</name>
<name>
<surname>Yuan</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Ingrahama</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Krampf</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Arrowooda</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Dadmohammadi</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Applying Macroscopic Material Balance to Evaluate Interplay between Dynamic Drainage Volume and Well Performance in Tight Formations</article-title>. <source>J. Nat. Gas Sci. Eng.</source> <volume>27</volume>, <fpage>466</fpage>&#x2013;<lpage>478</lpage>. <pub-id pub-id-type="doi">10.1016/j.jngse.2015.07.047</pub-id> </citation>
</ref>
<ref id="B17">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Ren</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Deng</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Lin</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2015</year>). &#x201c;<article-title>Refracturing Strategy for Sevensprings Oilfield in China Based on Geological and Engineering Characteristics</article-title>,&#x201d; in <conf-name>SPE/IATMI Asia Pacific Oil &#x26; Gas Conference and Exhibition</conf-name>, <conf-loc>Nusa Dua, Bali, Indonesia</conf-loc>, <conf-date>20-22 October 2015</conf-date>. <comment>paper 176163-MS presented at the</comment>. <pub-id pub-id-type="doi">10.2118/176163-MS</pub-id> </citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Siebrits</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Elbel</surname>
<given-names>J. L.</given-names>
</name>
<name>
<surname>Hoover</surname>
<given-names>R. S.</given-names>
</name>
<name>
<surname>Diyashev</surname>
<given-names>I. R.</given-names>
</name>
<name>
<surname>Griffin</surname>
<given-names>L. G.</given-names>
</name>
<name>
<surname>Demetrius</surname>
<given-names>S. L.</given-names>
</name>
<etal/>
</person-group> (<year>2000</year>). &#x201c;<article-title>Refracture Reorientation Enhances Gas Production in Barnett Shale Tight Gas Wells</article-title>,&#x201d; in <conf-name>Paper Presented at the SPE Annual Technical Conference and Exhibition</conf-name> (<publisher-loc>Dallas, TX</publisher-loc>: <publisher-name>OnePetro</publisher-name>). <pub-id pub-id-type="doi">10.2118/63030-MS</pub-id> </citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Vincent</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Restimulation of Unconventional Reservoirs: When Are Refracs Beneficial?</article-title> <source>J. Can. Petroleum Technol.</source> <volume>50</volume> (<issue>05</issue>), <fpage>36</fpage>&#x2013;<lpage>52</lpage>. <pub-id pub-id-type="doi">10.2118/136757-PA</pub-id> </citation>
</ref>
<ref id="B20">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Wan</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Rasouli</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Damjanac</surname>
<given-names>B.</given-names>
</name>
</person-group> (<year>2019</year>). &#x201c;<article-title>Numerical Simulation of Integrated Hydraulic Fracturing, Production and Refracturing Treatments in the Bakken Formation</article-title>,&#x201d; in <conf-name>53rd US Rock Mechanics/Geomechanics Symposium</conf-name> (<publisher-loc>New York, NY</publisher-loc>: <publisher-name>OnePetro</publisher-name>). </citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>Wei.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Application of Cluster Analysis in Selective Well Selection of Horizontal Well Repeated Fracturing[J]</article-title>. <source>Chem. Eng. Equip.</source> <volume>47</volume> (<issue>3</issue>), <fpage>76</fpage>&#x2013;<lpage>78</lpage>. </citation>
</ref>
<ref id="B22">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Williams</surname>
<given-names>G. J. J.</given-names>
</name>
<name>
<surname>Mansfield</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>MacDonald</surname>
<given-names>D. G.</given-names>
</name>
<name>
<surname>Bush</surname>
<given-names>M. D.</given-names>
</name>
</person-group> (<year>2004</year>). &#x201c;<article-title>Top-Down Reservoir Modelling</article-title>,&#x201d; in <conf-name>SPE annual technical conference and exhibition</conf-name>. <pub-id pub-id-type="doi">10.2118/89974-MS</pub-id> </citation>
</ref>
<ref id="B23">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Xia</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Mondal</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Fonseca</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2017</year>). &#x201c;<article-title>Understanding Stress Reorientation Process in Shale Gas Play and its Impact on Refracturing Time Window</article-title>,&#x201d; in <conf-name>51st US Rock Mechanics/Geomechanics Symposium</conf-name> (<publisher-loc>San Francisco, CA</publisher-loc>: <publisher-name>OnePetro</publisher-name>). </citation>
</ref>
<ref id="B24">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Yang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Xue</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Datta-Gupta</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>King</surname>
<given-names>M. J.</given-names>
</name>
</person-group> (<year>2016</year>). &#x201c;<article-title>Rapid Refracturing Candidate Selection in Shale Reservoirs Using Drainage Volume and Instantaneous Recovery Ratio</article-title>,&#x201d; in <conf-name>SPE/AAPG/SEG Unconventional Resources Technology Conference</conf-name>, <conf-loc>San Antonio, Texas, USA</conf-loc>, <conf-date>1-3 August 2016</conf-date>. <comment>paper 2459368-MS presented at the</comment>. <pub-id pub-id-type="doi">10.15530/URTEC-2016-2459368</pub-id> </citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yin</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Xie</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Datta-Gupta</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Hill</surname>
<given-names>A. D.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Improved Characterization and Performance Prediction of Shale Gas Wells by Integrating Stimulated Reservoir Volume and Dynamic Production Data</article-title>. <source>J. Petroleum Sci. Eng.</source> <volume>127</volume>, <fpage>124</fpage>&#x2013;<lpage>136</lpage>. <pub-id pub-id-type="doi">10.1016/j.petrol.2015.01.030</pub-id> </citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yuan</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Moghanloo</surname>
<given-names>R. G.</given-names>
</name>
<name>
<surname>Shariff</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Integrated Investigation of Dynamic Drainage Volume and Inflow Performance Relationship (Transient IPR) to Optimize Multistage Fractured Horizontal Wells in Tight/shale Formations[J]</article-title>. <source>J. Energy Resour. Technol.</source> <volume>138</volume> (<issue>5</issue>), <fpage>052901</fpage>. <pub-id pub-id-type="doi">10.1115/1.4032237</pub-id> </citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>An Exploratory Multi-Scale Framework to Reservoir Digital Twin</article-title>. <source>Adv. Geo-Energy Res.</source> <volume>5</volume> (<issue>3</issue>), <fpage>239</fpage>&#x2013;<lpage>251</lpage>. <pub-id pub-id-type="doi">10.46690/ager.2021.03.02</pub-id> </citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zheng</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Yuan</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Moghanloo</surname>
<given-names>R. G.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Analytical Modeling Dynamic Drainage Volume for Transient Flow towards Multi-Stage Fractured Wells in Composite Shale Reservoirs</article-title>. <source>J. Petroleum Sci. Eng.</source> <volume>149</volume>, <fpage>756</fpage>&#x2013;<lpage>764</lpage>. <pub-id pub-id-type="doi">10.1016/j.petrol.2016.11.023</pub-id> </citation>
</ref>
</ref-list>
</back>
</article>