Rate Decline of Acid Fracturing Stimulated Well in Bi-Zone Composite Carbonate Gas Reservoirs

: This paper develops a model of the multi-wing ﬁnite-conductivity fractures considering stress sensitivity for low-permeability bi-zone composite gas reservoirs. A new semi-analytical solution in the Laplace domain is presented. The main solution includes the theory of source function, Laplace integral transformation, perturbation technique, and Stehfest numerical inversion. Wellbore pressure is obtained by coupling solutions of reservoirs and fractures. The results showed that the pressure and derivative curves generated by this model include a bi-linear ﬂow stage. The model was validated by comparing its results with Wang’s results and the commercial well-test simulator; the results showed excellent agreement. This model illustrated the seepage characteristic of acid fracturing stimulated wells during refracturing treatment and how they are inﬂuenced by reservoir and hydraulic fractures parameters (asymmetrical factor and fractures distribution, etc.). The model is suitable to solve the solution of arbitrary-angle hydraulic fracture in refracturing and helpful to understand the transient production rate characteristic of the multi-wing fracturing well.


Introduction
Hydraulic fracturing is a common and effective practice in the production of lowpermeability gas reservoirs.At present, low gas and gas prices have aroused an interest in enhancing gas and gas production.Compared with drilling a new well, acid fracturing stimulated wells were considered a more budget-friendly way of increasing energy production [1,2].In refracturing, though the initial hydraulic fracture can have a significant influence on the propagation of a subsequent hydraulic fracture.The new fractures do not route at the same path as the initial fracture, owing to the change in stress anisotropy [3,4].
The pressure dynamic characteristics are widely used in order to evaluate the fracturecharacteristic parameters, which is thought as the most effective diagnostic method in gas and gas development.After hydraulic fracturing on the vertical well, there are often many secondary fractures distributing complexly around the wellbore, and there is an angle between the fracture and the fracture, and the length of every hydraulic fracture is not equal.The existence of complex fractures can effectively reduce the resistance that the fluid flows into the hydraulic fractures [5,6].
A large number of simulations and studies have been displayed [7,8] that the change of fracture parameters, especially fracture conductivity, will greatly affect well productivity.Dejam et al. [9] used the combination of the Laplace transform (LT) and finite Fourier Cosine transforms (FFCT) to solve the diffusion coefficient equation and proposed a semi-analytical solution for the transient pressure in vertical Wells in bounded double-pore formations.It showed that compared with the well location and reservoir shape of hydraulic fracturing, initiate at some angles along the direction of maximum stress, and various factors could have an effect on the direction and length of fractures.For the study of pressure dynamics on complex fractures formed during refracturing, the new seepage mathematical model of the reservoir and fractures about multi-wing fractures was established.The semi-analytical solution of the models is obtained by using the Laplace transform [22][23][24][25][26][27].
The analytical model of transient flow established by Shan et al. showed that orthogonal fractures generated during the refracturing process will produce orthogonal linear flow in the fractured reservoir volume.The angle between the wellbore and fracture can be estimated by using the pressure transient data analysis technique combined with the double logarithm curve diagnostic diagram [28].
Using the nodal analysis technique, Wu et al. [29] established a set of diffusion coefficient equations describing hydraulically redirected fracture flow.Then, the point source solution is coupled with the discrete fracture solution to obtain the semi-analytical solution.Five typical flow patterns can be observed in the transient response of this model.At the same time, the flux along the redirected fracture is calculated, and the flux distribution at different times is analyzed.On this basis, the influence curves of key parameters such as principal fracture angle, redirected fracture angle, permeability anisotropy, fracture conductivity, and fracture length ratio on transient seepage behavior are established.For rectangular anisotropic reservoirs, the influence of principal fracture angle, permeability anisotropy, and fracture conductivity on the curve type is concentrated in the early and middle stages, and the redirected fracture angle and fracture length ratio mainly affect the linear flow pattern.
Luo et al. [30] proposed a new semi-analytical method in the Laplace domain by coupling the Fredholm integral equation of the fracture with the reservoir equation and studying the unstable pressure characteristics of the "Z" shaped fracture.
Despite the above summary of the literature, the majority of efforts related to complex fractures' physical modeling or pressure behavior analysis do not take reservoir and fluid properties in different regions into account.Furthermore, the transient production rate per performance for low-permeability bi-zone reservoirs is rarely discussed.
In this research, the physical model of refractured wells with finite-conductivity fractures is derived in Section 2.Then, in Section 3, a mathematical model describing multi-wing fractures (MWFs) with bi-zone reservoirs accounting is derived.In Section 4, the mathematical model is solved through the Stehfest numerical inversion.The pressure and derivative curves are predicted.Results, discussions, and conclusions are provided in Sections 5 and 6, respectively.

Physical Model
In practical refracturing measurement, owing to different stress distributions around the wellbore and initial hydraulic fracture, MWFs with arbitrary angles besides initial hydraulic fracture are formed and produced at a constant rate q sc in a low-permeability gas reservoir (as shown in Figure 1).However, for convenience of research, it is assumed that width of all hydraulic fractures is equal and each fracture is penetrated fully.
The model presented in this article assumed that MWFs are with M fractures and with arbitrary angle between fracture and fracture.The well is located in center of bi-zone composite reservoirs with inner radius r m , and all fractures do not penetrate inner region.The infinite lateral boundary closed-upper and closed-bottom boundaries are chosen as the outer boundary in this model.

•
The fluid in the reservoir conforms to Darcy's law and law of isothermal percolation;

•
The upper and bottom boundaries of the reservoir are closed, and the fluid is incompressible;

•
Stress sensitivity of permeability in tight reservoirs is taken into account; • Negligible gravity and capillary effect are studied.

Line-Source Model in Bi-Zone Composite Reservoirs Considering Stress Sensitivity
Bi-zone composite low-permeability reservoirs with impermeable top and bottom boundaries have a line source.It is assumed that the fluid flow in the reservoirs results from an instantaneous pressure drop created at t = 0.The pressure waves propagate outward around the line source, which can make the fluids of inner and outer regions flow into the wellbore (as shown in Figure 2).
The majority of experimental studies demonstrate that during the creation of lowpermeability reservoirs, the permeability reduction brought on by stress sensitivity cannot be ignored.The common permeability change with pressure can be calculated with the following expression [31][32][33].It can be seen from Equation (1) that reservoir permeability is a function relating to pressure.Therefore, it presents a strong non-linearity in establishing the differential equation of low-permeability reservoir, which brings great difficulty to solving the solution of this model.Other basic assumptions are as follows:

•
The length of each hydraulic fracture is not equal, considering existence of angle (θ) between fracture and horizontal direction.The length (L f ), width (W f ), and permeability (k f ) of the different hydraulic fractures are different;

•
The fluid in the reservoir conforms to Darcy's law and law of isothermal percolation; • The upper and bottom boundaries of the reservoir are closed, and the fluid is incompressible;

•
Stress sensitivity of permeability in tight reservoirs is taken into account;

•
Negligible gravity and capillary effect are studied.

Line-Source Model in Bi-Zone Composite Reservoirs Considering Stress Sensitivity
Bi-zone composite low-permeability reservoirs with impermeable top and bottom boundaries have a line source.It is assumed that the fluid flow in the reservoirs results from an instantaneous pressure drop created at t = 0.The pressure waves propagate outward around the line source, which can make the fluids of inner and outer regions flow into the wellbore (as shown in Figure 2).According to corresponding study by Pedrosa [34] and Klkanl and Pedrosa [35], and combining equations of motion, equations of state, and continuous differential equations, the governing differential equations to describe bi-zone composite low-permeability reservoir of the inner and outer regions are established.According to definitions of dimensionless variables in Table A1, the dimensionless mathematical model of vertical linear The majority of experimental studies demonstrate that during the creation of lowpermeability reservoirs, the permeability reduction brought on by stress sensitivity cannot be ignored.The common permeability change with pressure can be calculated with the following expression [31][32][33].
It can be seen from Equation (1) that reservoir permeability is a function relating to pressure.Therefore, it presents a strong non-linearity in establishing the differential equation of low-permeability reservoir, which brings great difficulty to solving the solution of this model.
According to corresponding study by Pedrosa [34] and Klkanl and Pedrosa [35], and combining equations of motion, equations of state, and continuous differential equations, the governing differential equations to describe bi-zone composite low-permeability reservoir of the inner and outer regions are established.According to definitions of dimensionless variables in Table A1, the dimensionless mathematical model of vertical linear source with bi-zone composite low-permeability reservoir.Based on Laplace's integral transform in Equation ( 2), the dimensionless mathematical models in Laplace domain are obtained.
(1) Governing differential equation of inner region can be written as (2) Governing differential equation of outer region can be written as (3) Inner boundary condition can be written as lim (4) Infinite lateral boundary condition can be written as (5) Interface condition between the inner and the outer region can be written as The line-source model represented by Equations ( 4)-( 8) is strongly nonlinear.To successfully obtain transient pseudo-pressure responses caused by a line source in bi-zone composite low-permeability reservoir, perturbation technique [34] is used widely.When dimensionless permeability modulus is very small, zero-order approximate solution of the Energies 2023, 16, 2954 6 of 19 perturbation technique can satisfy the requirements of engineering precision.Zero-order approximate perturbation is as follows: Substitute Equation (10) into Equations ( 3)- (8), respectively.According to zero-order approximate perturbation given by Equations ( 9) and ( 3)-( 8) can be written as According to the theory of Bessel's function, we can get the general solution of governing differential equation of inner and outer region after perturbation as follows: ) where ε 0 = √ s; ε 1 = √ η 12 s.According to the properties of Bessel's function, With Equations ( 13)-( 16) and Equations ( 11) and ( 12) inputted into inner and outer boundary conditions and interface condition after perturbation, the coefficient of Equations (11) and (12) are obtained, respectively.
Energies 2023, 16, 2954 Substituting Equations ( 17) and (18) into Equation (11), pseudo-pressure is obtained at arbitrary position.However, it is noted that r D represents the dimensionless distance between the line source and any arbitrary point in low-permeability reservoirs in Equations ( 11) and ( 12) when line source is located at the origin.
If line source is located at (x w , y w ) instead of the origin, r D can be written as follows.

Surface-Source Model in Bi-Zone Composite Reservoirs Considering Stress Sensitivity
In order to obtain the uniform rate surface solution of the bi-zone composite lowpermeability reservoir, we can integrate Equation ( 13) to obtain the surface-source solution.
Based on assumption that hydraulic fracture is symmetric and parallel with the axis, the mathematical model solution is given to a single hydraulic fracture.However, if there is a certain angle between the hydraulic fractures and the coordinate axis or if the lengths of each hydraulic fracture are not equal, we can use the coordinate rotation and transformation method to distribute the hydraulic fractures along the changed coordinates.A schematic diagram of the physical model of the coordinate transformation is shown in Figure 3.
Substituting Equations ( 17) and (18) into Equation (11), pseudo-pressure is obtaine at arbitrary position.However, it is noted that rD represents the dimensionless distan between the line source and any arbitrary point in low-permeability reservoirs in Equ tions (11) and ( 12) when line source is located at the origin.(2

Surface-Source Model in Bi-Zone Composite Reservoirs Considering Stress Sensitivity
In order to obtain the uniform rate surface solution of the bi-zone composite low permeability reservoir, we can integrate Equation ( 13) to obtain the surface-source sol tion.Based on assumption that hydraulic fracture is symmetric and parallel with the ax the mathematical model solution is given to a single hydraulic fracture.However, if the is a certain angle between the hydraulic fractures and the coordinate axis or if the lengt of each hydraulic fracture are not equal, we can use the coordinate rotation and transfo mation method to distribute the hydraulic fractures along the changed coordinates.schematic diagram of the physical model of the coordinate transformation is shown Figure 3.The total pressure drop is obtained by the pressure drop superposition principle.A cording to previous research [36], the uniform flow surface solution of the ith hydrau fracture in the polar coordinates for the low-permeability composite reservoir can be wr ten as follows: The total pressure drop is obtained by the pressure drop superposition principle.According to previous research [36], the uniform flow surface solution of the ith hydraulic fracture in the polar coordinates for the low-permeability composite reservoir can be written as follows: According to the pressure drop superposition principle, the total pseudo-pressure response of M hydraulic fractures can be expressed as Energies 2023, 16, 2954 For a multi-wing well intercepting M hydraulic fractures, the assumption of constant total flow rate can yield the following equation: Based on the study of Cinco-Ley and Meng, we establish the following mathematical model to describe fluid flow in hydraulic fractures.Each hydraulic fracture is assumed as a three-dimensional rectangle with finite conductivity, and the volume of each hydraulic fracture is L f × W f × h.End point of each hydraulic fracture is impermeable.Compared with length of hydraulic fracture, width of hydraulic fracture is neglected.After the rotation of the coordinates, we established a mathematical model of low-permeability reservoirs with finite-conductivity hydraulic fracture in the x − y coordinate system.
With consideration of stress-sensitivity, the solution of ith fracture is obtained [25].

Coupling of Reservoir Model and Hydraulic Fracture Model
In order to obtain wellbore pseudo-pressure of low-permeability composite reservoir with multi-wing finite-conductivity fracture, combine Equations ( 23) and ( 26) to yield the final expression of dimensionless bottom-hole pseudo-pressure.
At the surface of the ith fracture, the reservoir pseudo-pressure and rate should be equal to those of the hydraulic fracture: Combining Equations ( 23) and ( 26), pseudo-pressure and rate at the surface of the ith fracture satisfy Equations ( 27) and (28).
However, the right-side portion of Equation ( 29) is only Fredholm integral equation.Because the rate of fracture is uneven during constant production, it is more difficult for us to solve wellbore pseudo-pressure of finite-conductivity fractures directly.We can divide fractures into N segments and assume that rate of each segment is equal (as is shown in Figure 4).According to the study of the Cinco-Ley [10], Equation ( 29) can be rewritten by discretizing single-wing fracture as follows.
Although we can only get M × N linear algebraic equations by Equation ( 30), M × N + 1 is still unknown.Therefore, if we want to get wellbore of bi-zone composite low-permeability reservoirs, a linear algebraic equation must be given.We can rewrite Equation (25) as following equation: Combining Equations ( 30) and (31), M × N + 1 linear algebraic equations, which can solve wellbore of bi-zone composite low-permeability reservoirs, are gotten.
When M = 2, the asymmetry factor is defined by using the length of the 2 wings.
Combining Duhamel's principle and pseudo-pressure drop superposition principle used in well-test analysis, wellbore pseudo-pressure is obtained considering skin effect and wellbore storage.
Finally, wellbore pseudo-pressure of multi-wing, formed by refracturing old hydraulic fracture, is obtained for bi-zone composite low-permeability reservoirs considering stress sensitivity.
The dimensionless wellbore flow rate equation for the constant pseudo-pressure production in low-permeability gas reservoirs can be determined by the relationship between the dimensionless pseudo-pressure and the rate in the Laplace domain [37]: According to the study of the Cinco-Ley [10], Equation ( 29) can be rewritten by discretizing single-wing fracture as follows.
Although we can only get M × N linear algebraic equations by Equation (30), M × N + 1 is still unknown.Therefore, if we want to get wellbore of bi-zone composite low-permeability reservoirs, a linear algebraic equation must be given.We can rewrite Equation (25) as following equation: Combining Equations ( 30) and ( 31), M × N + 1 linear algebraic equations, which can solve wellbore of bi-zone composite low-permeability reservoirs, are gotten.
When M = 2, the asymmetry factor is defined by using the length of the 2 wings.
Combining Duhamel's principle and pseudo-pressure drop superposition principle used in well-test analysis, wellbore pseudo-pressure is obtained considering skin effect and wellbore storage.
Finally, wellbore pseudo-pressure of multi-wing, formed by refracturing old hydraulic fracture, is obtained for bi-zone composite low-permeability reservoirs considering stress sensitivity.
The dimensionless wellbore flow rate equation for the constant pseudo-pressure production in low-permeability gas reservoirs can be determined by the relationship between the dimensionless pseudo-pressure and the rate in the Laplace domain [37]:

Mathematical Model Solution
Stehfest numerical inversion is used to obtain the pressure and rate of the real-time domain [38].In this section, our main purpose is to compare the solution obtained in our simple model with the commercial well-test simulator.The model without consideration of stress sensitivity in this article can be simply asymmetry finite-conductivity hydraulic fracture and finite-conductivity symmetry hydraulic fracture with bi-zone composite reservoir separately [39,40].Martirosyan, A.V. and Ilyushin, Y.V. presented the laboratory installation and mathematical model by changing the temperature field due to the similitude of the mathematical apparatus and the hydrodynamic processes behavior [41].Ahmadi et al. selected low-permeability carbonate plugs and applied nanoparticles in the porous media [42,43].If M = 2, L f1 = L f2 = 40, θ 1 = 0, θ 2 = 180 • in our model, it can be simplified as a conventional finite-conductivity symmetry hydraulic fracture with the bi-zone composite reservoir.To verify the credibility of the model presented in this article, we compared it with the commercial well-test simulator.If M = 2, L f1 = 20, L f2 = 60, θ 1 = 0, θ 2 = 180 • , and M 12 = η 12 = 1 in our model, we can again verify the correction of our solution by comparing it with Wang's model calculated by the asymmetry factor.
If our model can be simplified, it can be used to simulate a conventional finiteconductivity symmetry hydraulic fracture with a bi-zone composite reservoir.We can use a commercial well-test simulator to validate our model.Figure 5 shows that the results of our model are comparable to those of commercial well-test simulators for bi-linear flow, linear flow, and radial flow wells.The comparisons in Figure 5 show that the results calculated from our model are consistent with those obtained from the commercial simulator, confirming the credibility of the model presented in this article.

Mathematical Model Solution
Stehfest numerical inversion is used to obtain the pressure and rate of the real-time domain [38].In this section, our main purpose is to compare the solution obtained in our simple model with the commercial well-test simulator.The model without consideration of stress sensitivity in this article can be simply asymmetry finite-conductivity hydraulic fracture and finite-conductivity symmetry hydraulic fracture with bi-zone composite reservoir separately [39,40].Martirosyan, A.V. and Ilyushin, Y.V. presented the laboratory installation and mathematical model by changing the temperature field due to the similitude of the mathematical apparatus and the hydrodynamic processes behavior [41].Ahmadi et al. selected low-permeability carbonate plugs and applied nanoparticles in the porous media [42,43].If M = 2, Lf1 = Lf2 = 40, θ1 = 0, θ2 = 180° in our model, it can be simplified as a conventional finite-conductivity symmetry hydraulic fracture with the bi-zone composite reservoir.To verify the credibility of the model presented in this article, we compared it with the commercial well-test simulator.If M = 2, Lf1 = 20, Lf2 = 60, θ1 = 0, θ2 = 180°, and M12 = η12 = 1 in our model, we can again verify the correction of our solution by comparing it with Wang's model calculated by the asymmetry factor.
If our model can be simplified, it can be used to simulate a conventional finite-conductivity symmetry hydraulic fracture with a bi-zone composite reservoir.We can use a commercial well-test simulator to validate our model.Figure 5 shows that the results of our model are comparable to those of commercial well-test simulators for bi-linear flow, linear flow, and radial flow wells.The comparisons in Figure 5 show that the results calculated from our model are consistent with those obtained from the commercial simulator, confirming the credibility of the model presented in this article.In addition, we can verify our model and solutions by comparison with Wang's model mentioned above.Wang studied the pseudo-pressure dynamic characteristics of asymmetric hydraulic fracture by the means that well deviate from the center of the hydraulic fracture.If the angle between the hydraulic fracture and the x-axis is 180° and another angle between the hydraulic fracture and the x-axis is 0°, and the length of 2 hydraulic fractures is not equal in our model, the result of the model presented in this article should be the same as Wang's model.It is important to confirm the accuracy of our model under conductivity since the conductivity of the hydraulic fracture has a significant impact on the dimensionless pseudo-pressure and pseudo-pressure derivative during the early stages.Figure 6 shows that the results obtained in this article, under conductivity 1 and 10, are consistent with that of Wang's model, which can verify the accuracy of our model again.In addition, we can verify our model and solutions by comparison with Wang's model mentioned above.Wang studied the pseudo-pressure dynamic characteristics of asymmetric hydraulic fracture by the means that well deviate from the center of the hydraulic fracture.If the angle between the hydraulic fracture and the x-axis is 180 • and another angle between the hydraulic fracture and the x-axis is 0 • , and the length of 2 hydraulic fractures is not equal in our model, the result of the model presented in this article should be the same as Wang's model.It is important to confirm the accuracy of our model under conductivity since the conductivity of the hydraulic fracture has a significant impact on the dimensionless pseudo-pressure and pseudo-pressure derivative during the early stages.Figure 6 shows that the results obtained in this article, under conductivity 1 and 10, are consistent with that of Wang's model, which can verify the accuracy of our model again.

Results and Discussion
In order to study the flow regimes of low-permeability bi-zone composite reservo more graphically, type curves of pseudo-pressure response and production rate per mance are illustrated in Figures 7 and 8.According to the dimensionless pseudo-pressu derivate characteristic, pseudo-pressure response curves of low-permeability bi-zo composite reservoirs considering stress sensitivity are divided into sex flow periods.I noted that the red line presents the pseudo-pressure response curve of bi-zone compo reservoirs without considering stress sensitivity.
Period Ⅰ is the well storage and skin effect stage that occurred during the early peri both dimensionless pseudo-pressure and pseudo-pressure-derivative curves coincide a show as unit-slope straight lines during pure well storage.Curves of production rate a derivative coincide and exhibit a −1-slope line during pure well storage.
Period Ⅱ is the bi-linear flow stage corresponding to simultaneous linear flow in b hydraulic fractures and reservoirs.The pseudo-pressure-derivative curve exhibits a qu ter-slope line in an ideal condition.The production rate curve decreases, and the deri tive curve is shown as a level straight line whose value corresponds to conductivity.
Period Ⅲ is the early linear flow stage, in which the pseudo-pressure derivative cu appears with a half-slope.During this period, the fluid flow mainly appears as a lin flow perpendicular to the hydraulics fracture surfaces.The production rate curve ke decreasing, and the derivative curve also decreases compared with the bi-linear fl stage.
Period Ⅳ is the radial flow stage of the inner region.The pseudo-pressure wave d not propagate to the inner boundary.The fluid flow mainly appears as a radial fl around the whole MWF system during this stage.The pseudo-pressure derivative cu behaves as a 0.5-value horizontal line.Compared with the linear flow stage, the prod tion rate derivative curve drops significantly.
Period Ⅴ is the transition flow stage between the inner region and outer region.Period Ⅵ is the radial flow stage of the outer region.The pseudo-pressure wave pr agates to the inner boundary.The fluid of the outer region flow mainly appears as a rad flow around the whole MWFs system during this stage.Owing to different M12 value the inner and outer region, the pseudo-pressure derivative curve behaves as a 0.5 M value horizontal line without consideration of stress sensitivity.If stress sensitivity is c sidered, pseudo-pressure and pseudo-pressure derivative curves manifest as a line

Results and Discussion
In order to study the flow regimes of low-permeability bi-zone composite reservoirs more graphically, type curves of pseudo-pressure response and production rate performance are illustrated in Figures 7 and 8.According to the dimensionless pseudo-pressurederivate characteristic, pseudo-pressure response curves of low-permeability bi-zone composite reservoirs considering stress sensitivity are divided into sex flow periods.It is noted that the red line presents the pseudo-pressure response curve of bi-zone composite reservoirs without considering stress sensitivity.Period I is the well storage and skin effect stage that occurred during the early period; both dimensionless pseudo-pressure and pseudo-pressure-derivative curves coincide show as unit-slope straight lines during pure well storage.Curves of production rate and derivative coincide and exhibit a −1-slope line during pure well storage.Transient production rate curves influenced by the fracture number and dimensionless hydraulic fracture conductivity are shown in Figures 9 and 10.Fluid flowing into the wellbore through hydraulic fractures results in a pseudo-pressure loss, which is reflected by hydraulic fracture conductivity.Due to the tiny pseudo-pressure drop loss and short duration of bi-linear flow caused by large conductivity, the transient production rate and derivate curves diminish as the conductivity increases (as shown in Figure 9).More fractures can both boost output when there is constant pseudo-pressure and decrease pseudopressure loss when fluid is flowing into the wellbore from the reservoir.Figure 10 shows that fracture number only affects early-and intermediate-time transient production rate and derivate curves.The higher the fracture number is, the larger the transient production rate and derivate curves are during the early and intermediate time (as shown in Figure 10).Period II is the bi-linear flow stage corresponding to simultaneous linear flow in both hydraulic fractures and reservoirs.The pseudo-pressure-derivative curve exhibits a quarter-slope line in an ideal condition.The production rate curve decreases, and the derivative curve is shown as a level straight line whose value corresponds to conductivity.
Period III is the early linear flow stage, in which the pseudo-pressure derivative curve appears with a half-slope.During this period, the fluid flow mainly appears as a linear flow perpendicular to the hydraulics fracture surfaces.The production rate curve keeps decreasing, and the derivative curve also decreases compared with the bi-linear flow stage.
Period IV is the radial flow stage of the inner region.The pseudo-pressure wave does not propagate to the inner boundary.The fluid flow mainly appears as a radial flow around the whole MWF system during this stage.The pseudo-pressure derivative curve behaves as a 0.5-value horizontal line.Compared with the linear flow stage, the production rate derivative curve drops significantly.
Period V is the transition flow stage between the inner region and outer region.Period VI is the radial flow stage of the outer region.The pseudo-pressure wave propagates to the inner boundary.The fluid of the outer region flow mainly appears as a radial flow around the whole MWFs system during this stage.Owing to different M 12 values of the inner and outer region, the pseudo-pressure derivative curve behaves as a 0.5 M 12 -value horizontal line without consideration of stress sensitivity.If stress sensitivity is considered, pseudo-pressure and pseudo-pressure derivative curves manifest as a line upward trend.Stress sensitivity has little influence on the production rate derivative curve, but large stress sensitivity can lead to a small production rate curve in the log-log plot.
Transient production rate curves influenced by the fracture number and dimensionless hydraulic fracture conductivity are shown in Figures 9 and 10.Fluid flowing into the wellbore through hydraulic fractures results in a pseudo-pressure loss, which is reflected by hydraulic fracture conductivity.Due to the tiny pseudo-pressure drop loss and short duration of bi-linear flow caused by large conductivity, the transient production rate and derivate curves diminish as the conductivity increases (as shown in Figure 9).More fractures can both boost output when there is constant pseudo-pressure and decrease pseudo-pressure loss when fluid is flowing into the wellbore from the reservoir.Figure 10 shows that fracture number only affects early-and intermediate-time transient production rate and derivate curves.The higher the fracture number is, the larger the transient production rate and derivate curves are during the early and intermediate time (as shown in Figure 10).Figures 11 and 12 describe transient production rate curves affected by the mobi ratio (inner region to the outer region) and inner radius.The mobility ratio represents reservoir and fluid's physical property difference between the inner and outer regio When the mobility ratio is greater than 1, it indicates that the mobility of the inner is be than that of the outer and vice versa.A smaller mobility ratio can lead to increasing tr sient production rate curves during the radial flow stage of the outer region (as shown Figure 11).The inner radius has an influence on not only the duration of the inner rad flow but also the transient production rate curves of the outer radial flow.A small in radius corresponds to a short duration of inner radial flow and small transient product rate curves (as shown in Figure 12).Figures 11 and 12 describe transient production rate curves affected by the mobi ratio (inner region to the outer region) and inner radius.The mobility ratio represents reservoir and fluid's physical property difference between the inner and outer regio When the mobility ratio is greater than 1, it indicates that the mobility of the inner is be than that of the outer and vice versa.A smaller mobility ratio can lead to increasing tr sient production rate curves during the radial flow stage of the outer region (as shown Figure 11).The inner radius has an influence on not only the duration of the inner rad flow but also the transient production rate curves of the outer radial flow.A small in radius corresponds to a short duration of inner radial flow and small transient product rate curves (as shown in Figure 12).When the mobility ratio is greater than 1, it indicates that the mobility of the inner is better than that of the outer and vice versa.A smaller mobility ratio can lead to increasing transient production rate curves during the radial flow stage of the outer region (as shown in Figure 11).The inner radius has an influence on not only the duration of the inner radial flow but also the transient production rate curves of the outer radial flow.A small inner radius corresponds to a short duration of inner radial flow and small transient production rate curves (as shown in Figure 12).Figures 13 and 14 describe transient production rate curves affected by the fractur distribution and asymmetry factor.It is obvious that the effects of fracture distributio appear during the bi-linear and linear stages.When fracture clusters are distributed at a even angle (Model 3 in Figure 13), the contact area between the reservoir and the fractur is fully utilized, which contributes to an increase in production rate while pseudo-pres sure is kept constant.Smaller transient production rate curves during the early stage ar correlated with a more even fracture distribution (i.e., from Model 1 to Model 3 in Figur 13) (as shown in Figure 13).It is obvious that effects of fracture cluster asymmetry facto appear during the bi-linear and linear stages.The larger the fracture cluster asymmetr factor leads to smaller transient production rate curves during the bi-linear and linea stages (as shown in Figure 14).Figures 13 and 14 describe transient production rate curves affected by the fractur distribution and asymmetry factor.It is obvious that the effects of fracture distributio appear during the bi-linear and linear stages.When fracture clusters are distributed at a even angle (Model 3 in Figure 13), the contact area between the reservoir and the fractur is fully utilized, which contributes to an increase in production rate while pseudo-pres sure is kept constant.Smaller transient production rate curves during the early stage ar correlated with a more even fracture distribution (i.e., from Model 1 to Model 3 in Figur 13) (as shown in Figure 13).It is obvious that effects of fracture cluster asymmetry facto appear during the bi-linear and linear stages.The larger the fracture cluster asymmetr factor leads to smaller transient production rate curves during the bi-linear and linea stages (as shown in Figure 14).Figures 13 and 14 describe transient production rate curves affected by the fracture distribution and asymmetry factor.It is obvious that the effects of fracture distribution appear during the bi-linear and linear stages.When fracture clusters are distributed at an even angle (Model 3 in Figure 13), the contact area between the reservoir and the fracture is fully utilized, which contributes to an increase in production rate while pseudo-pressure is kept constant.Smaller transient production rate curves during the early stage are correlated with a more even fracture distribution (i.e., from Model 1 to Model 3 in Figure 13) (as shown in Figure 13).It is obvious that effects of fracture cluster asymmetry factor appear during the bi-linear and linear stages.The larger the fracture cluster asymmetry factor leads to smaller transient production rate curves during the bi-linear and linear stages (as shown in Figure 14).

Conclusions
In this paper, the semi-analytical solution of a finite-conductivity acid fracturing stimulated well formed during refracturing measurement is presented in lateral infinite for bi-zonal reservoirs.This paper considers the influence of stress sensitivity.Dimensionless pseudo-pressure and pseudo-pressure derivative curves under constant production and transient production rate and derivative curves under constant pseudo-pressure are drawn in a log-log plot.Based on the work presented in this paper, the following conclusions are obtained.
(1) Complicated fractures around the wellbore during refracturing measurement can be modeled using a modified multiple fractures model.Based on the theory of point function, when stress sensitivity is considered, the result of the multiple fractures model describing complicated fractures can be solved analytically in the Laplace domain.(2) We compared the simplified model obtained from this paper with the result calculated by Wang for asymmetry fracture without stress sensitivity and the result calculated by the commercial well-test simulator for symmetric fracture of bi-zonal

Conclusions
In this paper, the semi-analytical solution of a finite-conductivity acid fracturing stimulated well formed during refracturing measurement is presented in lateral infinite for bi-zonal reservoirs.This paper considers the influence of stress sensitivity.Dimensionless pseudo-pressure and pseudo-pressure derivative curves under constant production and transient production rate and derivative curves under constant pseudo-pressure are drawn in a log-log plot.Based on the work presented in this paper, the following conclusions are obtained.
(1) Complicated fractures around the wellbore during refracturing measurement can be modeled using a modified multiple fractures model.Based on the theory of point function, when stress sensitivity is considered, the result of the multiple fractures model describing complicated fractures can be solved analytically in the Laplace domain.(2) We compared the simplified model obtained from this paper with the result calculated by Wang for asymmetry fracture without stress sensitivity and the result calculated by the commercial well-test simulator for symmetric fracture of bi-zonal

Conclusions
In this paper, the semi-analytical solution of a finite-conductivity acid fracturing stimulated well formed during refracturing measurement is presented in lateral infinite for bi-zonal reservoirs.This paper considers the influence of stress sensitivity.Dimensionless pseudo-pressure and pseudo-pressure derivative curves under constant production and transient production rate and derivative curves under constant pseudo-pressure are drawn in a log-log plot.Based on the work presented in this paper, the following conclusions are obtained.
(1) Complicated fractures around the wellbore during refracturing measurement can be modeled using a modified multiple fractures model.Based on the theory of point function, when stress sensitivity is considered, the result of the multiple fractures model describing complicated fractures can be solved analytically in the Laplace domain.(2) We compared the simplified model obtained from this paper with the result calculated by Wang for asymmetry fracture without stress sensitivity and the result calculated by the commercial well-test simulator for symmetric fracture of bi-zonal composite reservoirs without stress sensitivity, respectively.The results showed excellent agreement.
(3) The log-log typical curves can be generated using a solution of this model, which mainly includes the bi-linear flow stage (quarter-slope portion), followed by the linear flow stage (quarter-slope portion), then the radial flow stage of the stimulated region (0.5 value) and radial flow stage of the un-stimulated region (0.5 M 12 value, without stress sensitivity).(4) The model illustrated how the transient production rate curves are influenced by reservoir and hydraulic fracture parameters.Reasonable fracture distribution can effectively decrease the pseudo-pressure loss of the early stage, and the more uneven the fracture distribution along the angle is, the lower the pseudo-pressure curve is; the smaller asymmetric factor leads to larger pseudo-pressure loss and the unobvious bi-linear characteristic.

Figure 1 .
Figure 1.Schematic of the physical model of low-permeability bi-zone composite reservoir.

Figure 1 .
Figure 1.Schematic of the physical model of low-permeability bi-zone composite reservoir.

Figure 2 .
Figure 2. Schematic of continuous line source for low-permeability bi-zonal composite reservoir.
is located at (xw, yw) instead of the origin, rD can be written as follows

Figure 3 .
Figure 3. Schematic of coordinate conversion and rotation.

Figure 3 .
Figure 3. Schematic of coordinate conversion and rotation.

Figure 4 .
Figure 4. Schematic of discretization of the hydraulic fracture.

Figure 4 .
Figure 4. Schematic of discretization of the hydraulic fracture.

Figure 5 .
Figure 5. Model verification with a commercial well test.

Figure 5 .
Figure 5. Model verification with a commercial well test.

Figure 8 .
Figure 8. Transient production rate curves for bi-zone low-permeability reservoirs.

Figure 9 .
Figure 9. Transient production rate curves affected by dimensionless hydraulic fracture conduc ity.

Figure 10 .
Figure 10.Transient production rate curves affected by fracture number.

Figure 10 .
Figure 10.Transient production rate curves affected by fracture number.

Figure 10 .
Figure 10.Transient production rate curves affected by fracture number.

Figures 11 and 12
Figures 11 and 12 describe transient production rate curves affected by the mobility ratio (inner region to the outer region) and inner radius.The mobility ratio represents the reservoir and fluid's physical property difference between the inner and outer regions.When the mobility ratio is greater than 1, it indicates that the mobility of the inner is better than that of the outer and vice versa.A smaller mobility ratio can lead to increasing transient production rate curves during the radial flow stage of the outer region (as shown in Figure11).The inner radius has an influence on not only the duration of the inner radial flow but also the transient production rate curves of the outer radial flow.A small inner radius corresponds to a short duration of inner radial flow and small transient production rate curves (as shown in Figure12).

Figure 11 .
Figure 11.Transient production rate curves affected by mobility ratio.

Figure 12 .
Figure 12.Transient production rate curves affected by inner radius.

Figure 12 .
Figure 12.Transient production rate curves affected by inner radius.

Figure 12 .
Figure 12.Transient production rate curves affected by inner radius.

Figure 13 .
Figure 13.Transient production rate curves affected by fracture distribution.

Figure 14 .
Figure 14.Transient production rate curves affected by asymmetry factor.

Figure 13 .
Figure 13.Transient production rate curves affected by fracture distribution.

Figure 14 .
Figure 14.Transient production rate curves affected by asymmetry factor.

Figure 14 .
Figure 14.Transient production rate curves affected by asymmetry factor.
Transient production rate curves affected by fracture distribution.

Table A1 .
Angle between the fracture and the horizontal axis, degree r m Radius of inner region, m x, y, z Directional coordinates, m ψ 0D1 Zero-order perturbation deformation function of inner region ξ 0DF Zero-order perturbation deformation function of hydraulic fracture θ Angle between the fracture and the horizontal axis, degree r m Radius of inner region, m x w , y w xand y-coordinates of line source, m M Hydraulic fractures number, integer N Grid number of hydraulic fractures dispersed, integer q Strength of continuous line source, m 2 /s q f Strength of continuous fracture line source, m 2 /s Node i-th hydraulic fractures is the j-th grid r mi,jMiddle node i-th hydraulic fractures is the j-th gridK v (x)The second kind modified Bessel function I v (x) The first kind modified Bessel function s Definitions of Dimensionless Variables.2D = πk 1 hT sc p sc Tq sc (ψ e − ψ 2 ) = πk 1 hT sc p sc Tq sc (ψ e − ψ f ) 1 φ 1 C t1 /k 1 µ 2 φ 2 C t2 /k 2 µ p sc Tq sc πk 1 hT sc γ