Modeling and Simulation of Micro/NanoScale Oil and Gas Migration in Tight Rocks and Shales
View this Special IssueResearch Article  Open Access
Qi Chen, Shaojun Wang, Dan Zhu, Guoxuan Ren, Yuan Zhang, Jinghong Hu, "A Comprehensive Model for Estimating Stimulated Reservoir Volume Based on Flowback Data in Shale Gas Reservoirs", Geofluids, vol. 2020, Article ID 8886988, 14 pages, 2020. https://doi.org/10.1155/2020/8886988
A Comprehensive Model for Estimating Stimulated Reservoir Volume Based on Flowback Data in Shale Gas Reservoirs
Abstract
Stimulated reservoir volume (SRV) which is generated by horizontal drilling with multistage hydraulic fracturing governs the production in the shale gas reservoirs. Although microseismic data has been used to estimate the SRV, it is highpriced and sometimes overestimated. Additionally, the effect of stress sensitivity on SRV is not considered in abnormal overpressure areas. Thus, the objective of this work is to characterize subsurface fracture networks with stress sensitivity of permeability through the shale gas well production data of the early flowback stage. The flowback regions are first identified with the flowback data of two shale gas wells in South China. Then, we measured the permeability stress sensitivity of the core after fracturing, coupled to the dynamic relative permeability (DRP) calculation to obtain an accurate and simple DRP curve. After that, a comprehensive model is built considering dynamic twophase relative permeability function and stress sensitivity. Finally, we compared the calculated results with the microseismic data. The results show that the proposed model could reasonably predict the SRV using the flowback data after fracturing. Additionally, compared with the microseismic data, the stress sensitivity should be included, especially in the abnormal overpressure block. It is believed that this mathematical model is accurate and useful. The work provides an efficient approach to estimate stimulated reservoir volume in the shale gas reservoirs.
1. Introduction
The development of shale gas has gained increasing attention with the decline of the production from conventional reservoirs [1, 2]. Due to the ultralow permeability of shales (from ) [3], a combination of horizontal drilling and multistage hydraulic fracturing has been widely used to SRV, which increases effective contact significantly to improve gas production from an individual well. Therefore, characterization of the stimulated reservoir volume appears to be of vital importance to manage and predict shale reservoir performance [4, 5].
While microseismic data has been extensively used to validate the stimulated reservoir volume, the SRV obtained directly from the microseismic data is generally overestimated. For example, SRV would include large unstimulated void regions by the convexity assumption [6–9]. Also, implementation of microseismic techniques in a field sometimes is not costeffective, which limits a wide application in fields. Therefore, using early production data after hydraulic fracturing together with early flowback data of hydraulic fluids is fast becoming a key instrument in evaluating SRV [10–13], namely, transient analysis of flowback data.
Twophase flow appears to be a classic and representative flow behavior during the flowback process after hydraulic fracturing [2, 14]. For example, a large number of production data obtained from wells in the Horn River shale show that twophase flow was observed for every single well during the flowback process after hydraulic fracturing [15, 16]. Clarkson [17] and WilliamsKovacs and Clarkson [18] also discovered twophase flow within 48 hours of well opening and production in the Barnett and Marcellus shales. Similarly, the twophase flow stage occurs in the Silurian Longmaxi Formation shale gas wells in South China after the well opening [19]. And an even greater advance in the research of production data patterns is that after processing the production data from shale gas wells in the Horn River Basin, Adefidipe et al. [20, 21] and Xu et al. [22, 23] divided the production data into two stages as early gas production (EGP) and late gas production (LGP) for the “V” shape of the relationship between the gaswater ratio (GWR) and the cumulative gas production. They believed that the EGP stage was special and worth studying because the LGP stage already had a matrix gas supply.
To take advantage of the twophase flow characteristics during the early gas production (EGP) stage, research studies have endeavored to estimate the SRV by coupling rock physics with various types of reservoir models. For example, Ezulike et al. [24–26] integrated an explicitly determined dynamic relative permeability (DRP) function in linear dualporosity models to explore the twophase flowback data. The function extends the singlephase water flow equation to gaswater twophase flow which was calculated using cumulative production data of water and hydrocarbon compounds. Clarkson et al. [17, 27, 28] initially used an improved method for coalbed methane (CBM) flow to analyze multiphase flows. They assumed that the flow of gas and water through the generated fracture network was similar to the simultaneous flow of gas and water during longterm production from fractured coal reservoirs. Later, they proposed a large number of twophase flow models based on this assumption to quantitatively calculate multiflow data, including tight oil and tight gas, to help analyze the flow characteristics of shale gas wells which are analogous to analyzing twophase (gas and water) CBM production. And in the followup research work [29–32], they verified the analysis results from the above methods by integrating various techniques, including a modified material balance equation (MBE) applied to beforebreakthrough (BBT) (before gas breakthrough) singlephase ratetransient analysis (RTA) and afterbreakthrough (ABT) (after gas breakthrough) multiphase ratetransient analysis (RTA).
While the aforementioned models provide insights in calculating SRV through multiflow data quantitatively, challenges of applying these models to the field remain because the Langmuir volume needs to be increased in the calculation process, which will lead to some extent uncertainty and inaccuracy. To overcome the complexity of applying these models in fields, Clarkson et al. [17, 28] and Fu et al. [33] established a linear matching curve relation to describe the straight line relationship between ratenormalized pressure (RNP) and material balance time (MBT). This is because fracture permeability and fracture volume can be obtained on the basis of the characteristics of the fitting curve. Also, Alkouh et al. [34] modified a loglog plot of RNP and RNP derivative against MBT for both the oil and gas cases and provided a straightforward approach to calculate the key fracture parameters covering fracture permeability and fracture halflength. However, the influence of fracture compressibility during the flowback process is not taken into consideration in the models, which likely leads to estimation error of the total compressibility.
However, due to some defects in the hypothesis of DRP function, a nonlinear phenomenon occurs when the model is applied to some field cases. And the phenomenon of the nonzero intercept leads researchers to question the physical properties of these problems. It may be due to the inaccuracies caused by stress sensitivity in some areas. In particular, there is abnormal overpressure in shale gas reservoirs in southern China [35], so the influence of stress sensitivity, especially for the twophase relative permeability, has to be taken into account [36]. Moreover, all the researchers did not make a production data verification and comparison with the microseismic data.
We thus aimed to develop a stressdependent twophase relative permeability, which can be used to extend the existing singlephase model. Also, we proposed a new multiphase model to estimate SRV and fracture characteristics. Moreover, to verify our model, we compared the results with the microseismic data.
In this work, we first processed the production data of shale gas wells and obtain the gaswater ratio characteristics. Then, we calculated the simple and practical twophase dynamic relative permeability (DRP) based on the production data collected in the field belonging to shale gas wells in southern Sichuan. Subsequently, the stressstrain curves of permeability are obtained by using the laboratory experiments. Finally, the material balance equation and diffusion equation are used to process the production data to derive the model, thus calculating fracture parameters. We compared the calculation results with the microseismic data to verify the rationality and accuracy of the calculation results.
2. Data Preparation
2.1. Field Data from a Typical Shale Gas Reservoir in Southern China
In this work, we selected the flowback data from two production wells located in southern China: the shale gas formations of the Lower Silurianare, a typical marine shale gas reservoir. A large number of micronnanoscale pores are developed in the mineral grains and organic matter [37]. These pores are the main channels for shale gas storage and migration, which play an important role in improving shale gas storage performance and providing good storage space for marine shale gas accumulation and enrichment [38]. The depth of the shale gas reservoir is over 2000 m at subsurface with a formation pressure coefficient up to 2.0. The thickness of the entire shale reservoir is about 100 m with porosity ranging from 2% to 4% and matrix permeability of 0.001 mD. The production at the early stage was about with a stable production around at .
The eastwest anticline structure belt with few faults is arranged in the left echelon row as the main geological structure of this area [39]. This drilling platform is located at the end of the anticline structure. Although two compressoshear faults exist around the well pad, the overall upper lateral sealing performance leads to an effective accumulation and preservation of shale gas [37, 38]. Figure 1 shows the layout of the horizontal well group. Due to the low quartz content in the highquality shale segment at the bottom of the Longmaxi Formation, natural fractures are not developed throughout the reservoir [40]. In the whole horizontal well group, Well 1, Well 4, and Well 5 are located in one formation, while Well 2, Well 3, Well 6, and Well 7 belong to another formation.
In this work, we used data from Well 3 and Well 4 to test our model for the following two reasons. Firstly, the two wells were drilled through the shale gas reservoir in different directions and are located in two different formations. Secondly, microseismic monitoring is conducted in the two wells together with a comprehensive drilling, completion data, and well test data. Figure 2 shows gas and water production rates of the two wells. The two wells both experienced gas breakthrough immediately after well opening and production, meaning that there was no singlephase flow. The difference is that Well 3 experienced a precipitous drop in production at one point (80 hrs), while Well 4 experienced a more gentle fluctuation, possibly due to different formation conditions. We are sorry that there is no further research on the causes of this phenomenon in this paper on account of insufficient highfrequency production data or complete adjacent test well data.
(a)
(b)
2.2. Producing GasWater Ratios (GWR)
In order to explore the trend of GWR of production data, this section simply processes and analyzes the production data to observe whether there is an immediate gas breakthrough after the production of shale gas wells [17, 27], which means that SRV can be calculated by establishing relevant models based on production data in this area.
This is largely because the effective fracture network system is saturated with both the gas and water phases after two shutins. The gas source here is assumed to be from three aspects: (1) the originally existing initial gas in the active natural fracture, (2) the gas displaced by fracturing fluid under the influence of the strong countercurrent water imbibitions into the shale matrix due to the huge pressure difference during the first shutin period, and (3) the gas accumulation that resulted from spontaneous imbibitions of fracturing fluid during the second shutin period [41, 42].
Similarly, before establishing the model, we also need to use the gaswater yield ratio to gain the Vshaped trend. Thus, we processed the initial twophase production data of Well 3 and Well 4 and finally obtained the trend of the gaswater ratio (GWR). Figure 3 shows a Vshaped trend GWR over time as a feature of early twophase flowback data. Although the inflection points in the GWR curves of the two wells occur at different times due to different formation conditions and production systems, the same Vshaped trend still appears obviously. Ghanbari et al. [15] and Abbasi [16] also found the Vshaped GWR in the wells drilled in the Horn River shales in Canada.
(a)
(b)
The gaswater ratio decreases and then increases, and we can study it from the most basic theory. Assume that the gaswater twophase flow satisfies Darcy’s law: where is the fluid flow rate under pressure , m^{3}/s; is the fluid permeability, μm^{2}; is the fluid flow area, m^{2}; is the fluid viscosity, mP·s; is the pressure, MPa; and is the fluid flow length, m.
After ignoring the capillary force in the fracture network system, the gaswater ratio becomes where and are the flow rates of the gas phase and water phase, respectively, m^{3}/d; and are the viscosity of the gas phase and water phase under pressure, respectively, mP·s; and are the relative permeability of the gas phase and water phase under their respective saturation, μm^{2}; and are the saturation of the gas phase and water phase, dimensionless; and and are the pressure of the gas phase and water phase, respectively, MPa.
After the well opening for production, the gas viscosity decreases with the decrease of pressure, while the water viscosity remains relatively unchanged. According to the change of the gaswater ratio over time, the ratio of gaswater relative permeability decreases. Therefore, on the basis of the relationship between permeability and saturation, it can be inferred that the saturation ratio has a corresponding variation trend, meaning that the initial gas saturation () is not zero, which also proves the existence of initial free gas in the fracture before well opening. When the wellbore storage effect disappears, the gas flows from the matrix to the fracture network, resulting in the replenishment of gas in the fracture system. Accordingly, the gas saturation increases, so the gaswater relative permeability increases and the gaswater ratio decreases. This is why the slope of the GWR curve drops first and then rises.
In general, the negative slope on the GWR diagnostic plot is called early gas production (EGP) and the rise of the GWR diagnostic plot is called late gas production (LGP). It is generally believed that the EGP region is the stage of wellbore storage effects and both the gas and water productions come from the effective fracture connected with the horizontal well. Also, the increase of water relative permeability is greater than the reduction of gas viscosity. The LGP phase is the result of the matrix gas transfer to the fracture network after wellbore effects becoming negligible [10, 11, 17, 25–27]. According to the above analysis, we have figured out that in the EGP stage, the ground flow of flowback is only from the effective fracture system, while after the LGP stage beginning, the matrix gas initiates participation in the flow. Therefore, we will establish an analysis model for the EGP phase to calculate the SRV of shale gas wells after multistage fracturing.
3. Methods
3.1. Material Balance Equation (MBE) and Diffusion Equation for the Fracture System
In order to facilitate the establishment of the mathematical model for the EGP phase, the fracture network around the shale gas fractured well is simplified into the SRV region composed of the matrix system and fracture system. As shown in Figure 4, the effective fracture system is made up of a fracture section and a matrix section. The fracture section consists of artificial hydraulic fractures, the secondary fractures generated by fracturing, and active natural fractures. The matrix section is the shale matrix connected to the fracture section. It is worth noting that the desorption and adsorption processes were ignored in the material balance equation (MBE) [17, 43].
This model simplifies complex, active natural and secondary fractures, as well as artificial hydraulic fractures, into a simple fracture system. The length of the artificial hydraulic fracture is used as the width of the entire stimulated reservoir volume, and the length of the horizontal wellbore is used as the length. In the whole stimulated reservoir area, the height of the major fracture, including the matrix part connected to the fracture system, participates in the flow. In the equivalent fracture system, the fracture is saturated with fracturing fluid (water phase) and natural gas (gas phase). It is assumed that no matrix gas is involved in the flow during the early gas production (EGP).
It is assumed that the fracture system can be approximated as a homogeneous/closed/tank system. And the fluid flow from the fracture to the horizontal well is assumed to be linear. The mechanism driving the gaswater flow includes two aspects: (1) fracture closure and (2) expansion of the fluid (gaswater phase). Kuchuk et al. [44] also considered the assumption that fracture closure and expansion of water and free gas in the fracture system drive the accumulation of natural gas and water at the surface suggesting the occurrence of pseudosteadystate flow in fractured reservoirs.
In general, for the purpose of facilitating the establishment and solution of the model, we made the following assumptions: (1) capillary pressure in fracture systems is ignored; (2) the gas from the matrix is negligible in the EGP stage; (3) the fracture system is approximated as a homogeneous/closed system; (4) Darcy’s law applies to fluid flow; and (5) the effective fracture system is saturated with fracturing fluid (water phase) and natural gas (gas phase) initially.
Effective compressibility terms are defined using simplified gas material balance equations (MBEs) [45–48]. In this work, the material balance equation (MBE) of the gas phase in the fracture network is simplified. The effects of gas and water (fracturing fluid) expansion and fracture closure are classified as the “effective compressibility” term given by where and are the cumulative gas production and initial volume of gas in the fracture, respectively, m^{3}; and are the gas formation volume factor and gas formation volume factor at initial conditions, respectively, m^{3}/m^{3}; and are the initial gas phase saturation and initial water phase saturation, respectively, dimensionless; and are the gas compressibility and water compressibility, respectively, Pa^{1}; and are the cumulative water production and initial volume of water in the fracture, respectively, m^{3}; and are the volume of effective fractures and volume of effective fractures at initial conditions, m^{3}; and is the fracture pressure, Pa.
is analogous to the total compressibility term proposed for conventional multiphase well testing. This is represented as a function of measurable flowback parameters, including cumulative gas production and cumulative water production. Each term in its expression outlines different driving mechanisms in the fracture system: (1) expansion of the gas phase, (2) expansion of liquid phase, and (3) fracture closure. Among them, is similar to the formation compressibility coefficient in traditional conventional material balance analysis [49]. When dealing with fractures, it refers to the inverse of the fracture stiffness determined by the elastic or strain energy required to keep the hydraulic fracture open [50].
Ignoring the gas flow rate from the matrix system to the fracture, the gas phase material balance equation is where and are the density of gas at surface conditions and reservoir conditions, respectively, kg/m^{3}; and is the volume of gas in the fracture system at any time, m^{3}.
According to the relationship between gas volume , fracture volume , and water volume in the fracture system, the definition of gas compressibility is obtained as where , , and are the fracture volume, gas volume, and water volume in the fracture system, respectively, m^{3}.
By substituting the gas production and water production , then using the chain rule, the total volume of fractures is derived: where is the water formation volume factor at initial conditions, m^{3}/m^{3}.
In the early gas production (EGP), it is assumed that , the final material balance equation, can be obtained by defining a compressibility term and the total fluid flow rate by simplifying the above equation is where is total producing rate, m^{3}/d.
Given that singlephase, steadystate flow can be described using the continuity equation and Darcy’s law, the singlephase gas diffusion equation in the fracture system is given by Zhang and Winter [51]: where is the gas compressibility factor, dimensionless; is the porosity for the fracture system, dimensionless; and is the gradient operator.
Define the pseudopressure and pseudotime functions [46, 52] as where is the relative permeability of the gas at a given time, μm^{2}; and is the total effective compressibility, Pa^{1}.
Then, the governing equation of singlephase gas flow in the fracture system is
3.2. Coupling Stress Sensitivity with Dynamic Relative Permeability Function
In Equation (10), is the variable function of singlephase gas with time. By introducing the dynamic relative permeability (DRP) function , the diffusion equation of singlephase gas is transformed into the diffusion equation of the twophase system [53].
The method in this paper is similar to that in Ezulike and Dehghanpour’s study [53], but due to the geological conditions of formations belonging to southern Sichuan being different from those of North America, the dynamic relative permeability (DRP) function of the gas phase is obtained by using the method of fitting field data. In the previous study, we have assumed that the fracture system is saturated with gas and water. The relationship between cumulative gas and water production and time can be reported through field data. Therefore, the gaswater saturation in the fracture system can be established as a function of gaswater production with time. Given the characteristics of abnormal overpressure in this area, the stress sensitivity in this area cannot be ignored [36, 54, 55]. Therefore, the stresssensitive model was considered in the dynamic relative permeability function. The specific process is shown in Figure 5.
In Figure 6, the relative permeability of the cores after fracturing in the research area is tested through laboratory experiments. The curve of pressure change with time during the flowback period is shown in Figure 7. Stressstrain experiments were performed to describe the curves between dimensionless permeability and effective stress by stabilizing the confining pressure and reducing the internal pressure. This is more in line with the actual formation flowback pressure changes. The early flowback stage of shale gas wells after well opening is significantly short; therefore, it can be assumed that the formation pressure has not changed much, but the inner pressure begins to decrease after the fluid flows out of the cores.
(a)
(b)
Take the stress sensitivity experiments of cores belonging to the stimulated area of Well 3 as an example. Both the matrix and the fractures are most likely contained in these cores. We had required as many experiments as possible to obtain the stresssensitive test data in the stimulated region after fracturing as accurately as possible. However, subject to the insufficient samples or the different distances of samples from the horizontal well, these data maybe cannot represent the whole area in fact but still have considerable reference and research value.
The specific experimental procedures are as follows: (1) the initial confining pressure was set as the original formation pressure which is 38 MPa, and the internal pressure was 23 MPa; (2) the confining pressure was increased to 58 MPa slowly which is the formation pressure before well opening for flowback, and the internal pressure was increased to 43 MPa at the same time to keep the effective stress constant; and (3) the internal pressure was reduced to different pressure points to increase the effective stress, and the gas permeability of the sample was measured after each pressure point was stabilized.
Dimensionless permeability is defined as where is the relative permeability of the gas at a given time, μm^{2}; is the initial gas relative permeability, μm^{2}; and are are stress sensitivity coefficients by experiments, dimensionless.
The experimental results are reported in Figure 8. According to the calculation method in Figure 5, the results in Figure 9 were finally obtained. Then, according to the fitting curve, the function of gasphase relative permeability changing with time is
(a)
(b)
Since the above equation is obtained based on flowback data and core experiment relative permeability, and are called the flowback data coefficient and flowback data intercept, respectively, which are determined by the data fitting curve.
3.3. TwoPhase Flowback Model for Early Gas Production (EGP)
By substituting Equation (13) into Equation (10), the governing equation of singlephase gas flow in the fracture system is changed into the twophase flow model. This means that when gaswater twophase flow occurs in the fracture system, the new real pseudotime function is transformed into
In order to establish the relationship between pseudopressure function and pseudotime function, we make the following transformation:
And is derived by the final material balance equation (Equation (7)); in a similar way, we can gain and by the defined pseudopressure (Equation (9)) and new pseudotime equation (Equation (14)), respectively. Finally, Equation (15) becomes
Define the equivalent gas rate as [48]
Substituting in Equation (16) using the real gas law, then combine it with the equivalent gas rate (Equation (17)):
By substituting (Equation (18)) in Equation (11), we can establish the twophase diffusion equation:
Then, the following inner and outer boundary conditions can be solved: , , and , .
Define the fracture storage coefficient as
Then, the following relationship is gained:
Substitute Equation (20) into Equation (18) and integrate to obtain
By combining Equation (22) with Equation (15), we can modify the final twophase flow model of the early gas production stage as follows:
Theoretically, a plot of the ratenormalized pseudopressure (RNP) vs. the pseudotime should yield a straight line relationship. With the pseudotime function as an independent variable and the left side of the equation as a dependent variable, the slope and intercept can be obtained according to the fitting curve, and the relationship between the equivalent fracture porosity/halflength of the effective fracture system and effective fracture system permeability can be described as follows:
3.4. Analysis Procedure
We propose the following analysis procedure: (1)Obtain and process water and gas flowback data to explore a Vshaped gaswater ratio trend (see Figures 2 and 3 and Equation (2))(2)The early flowback period (EGP) is distinguished from the late flowback period (LGP) according to the Vshaped trend of the GWR curve (see Figure 3)(3)Conduct a simpler fracture network system model for the EGP (see Figure 4)(4)Calculate effective compressibility by Equation (3) to deduce the twophase material balance equation (MBE) (Equation (7)) for the fracture system(5)Define the pseudopressure and pseudotime functions (Equations (9) and (10)) to deduce the gas diffusion equation (Equation (11)) of the fracture system(6)Plot the gaswater relative permeability curve (see Figure 5) and the curves of dimensionless permeability with effective stress (see Figure 8) by experiments(7)Calculate gas DRP coupling stress sensitivity (see Figure 9) following the steps by Figure 6 to transform the gas diffusion equation (Equation (11)) into the twophase diffusion equation (Equation (19))(8)Modify the final twophase flow mode (Equation (23)) by combining Equation (22) with Equation (15)(9)Plot ratenormalized pressure change with pseudotime (see Figure 10)(10)Calculate SRV by Equation (24) and effective fracture system permeability by Equation (25)(11)Verify the analytical model against microseismic data (see Figure 11)
(a)
(b)
(a)
(b)
The flowback data we need to obtain include production rates and pressure and cumulative production data profiles. Then, we got a Vshaped trend in the gaswater ratio curve by processing the flowback data. Thus, an analytical model was established for the EGP stage. Finally, the calculated results were compared with the microseismic data to validate the mathematical results.
4. Results and Discussions
4.1. Stimulated Reservoir Volume
We apply the analytical model presented above to analyze the flowback data of Well 3 and Well 4. However, there are several issues that need to be addressed and discussed: (1)It is difficult to gain an appropriate initial gas saturation of the fracture system from actual field data. Unlike conventional numerical simulation, this parameter is unknown in actual field data. After the fracturing operation is completed, the effective fracture system connecting the wellbore is filled with fracturing fluid approximately. Thus, after the well is opened for flowback operation and before production, the volume of the recovered fracturing fluid under the ground is filled with gas renewedly. Therefore, in this paper, we take a reasonable value which is the recovery percentage of the total injection amount of fracturing fluid as the initial fracture system gas saturation in calculation. According to Xu’s paper [48], there is a little significant difference between the actual value of initial gas saturation and the recovery percentage. However, this was still not rigorous enough. Then, other possible values are substituted into the calculation by us to get the new SRV. The results showed that the impact of this value on SRV estimation was not obvious within the error range, and the difference was less than 1% compared with the microseismic data
Due to that, fracturing fluid may leak off into the existing inactive natural fractures and into the matrix during injection; the fluid flowback percentage decreases, but the initial gas saturation in the fracture system increases. Hence, when a significant portion of the injected fluid volume does not contribute to create fracture volume, this value of the initial gas saturation for calculation is actually low. One direction of future work is to consider using the results to iteratively optimize the initial gas saturation value. (2)The value of the fracture closure term in the total compressibility cannot be accurately expressed. Fortunately, it is found that this value has no obvious influence on the results in the calculation process. Since this value refers to the inverse of fracture stiffness when dealing with fractures, it is considered to give a reasonable value in the subsequent work from the perspective of rock fracture mechanics
The final model calculation results are shown in Figure 10. Since microseismic detection technology was used to estimate the final hydraulic fracture stimulated effect in the two wells. We compared the calculated results with the microseismic data and finally found that the difference between the calculated results of the analytical model and the microseismic data was minor, as shown in Figure 11. This reason is that quantities of fractures that do not connect effectively with fracture systems are counted in the final result when microseismic detection techniques are used. Those “dead” fractures skew the microseismic results.
4.2. Effective Fracture System Permeability of Stimulated Reservoir
As shown in Table 1, after obtaining the slope of ratenormalized pressure change with pseudotime (see Figure 10), SRV was gained by Equation (17). On the basis of the average length of hydraulic fractures provided by the microseismic data of the drilling platform where Well 3 and Well 4 are located and the average porosity of the cores after fracturing obtained through experiments, the permeability of the effective fracture system in the stimulated region is calculated by using Equation (18). The results are shown in Table 1, where the absolute deviation refers to the difference between the SRV calculated by the analytical model and estimated by the microseismic data, and the relative deviation refers to the ratio of the absolute deviation to the SRV estimated by microseismic data. Numerical errors in porosity of the effective fracture system are responsible for the error of permeability.

If the SRV estimated by microseismic data is taken as a reference for comparison, the specific calculation process is as follows: where and are the absolute deviation and the relative deviation between the SRV calculated by the analytical model and estimated by the microseismic data, respectively, m^{3}; and and are the SRV calculated by the analytical model and the SRV estimated by the microseismic data, m^{3}.
4.3. Analytical Model without Considering Stress Sensitivity
If the stress sensitivity effect is not considered in the process of the DRP calculation, the deviation of the DRP curve will occur, as shown in Figure 12. And the curve with a large error was substituted into the step shown in Section 3.4 to get the result in Figure 13. As a result, both the slope and the intercept of the RNP curve had changed. We obtained the new slope and intercept of the curve and calculated a new SRV according to Equations (17) and (18). Then, the calculation results and corresponding deviations as shown in Figure 14 were obtained.
(a)
(b)
(a)
(b)
It can be seen that the relative deviation calculated without considering the stress sensitivity effect is greater than that calculated with considering the stress sensitivity effect. Moreover, because the real SRV has been overestimated by the microseismic data, the calculation results of the model without considering the stress sensitivity effect are higher than those of the microseismic data, which further indicates that the calculation results without considering the stress sensitivity effect in the highpressure area will have a large error, which cannot be ignored. We also calculated the permeability of the effective fracture system, and the effective permeability of Well 3 and Well 4 was 0.432 mD and 0.283 mD, respectively. It can be seen that the calculation results are 43.50% and 15.75% more than those considering the stress sensitivity effect.
5. Conclusions
In this paper, a comprehensive model is developed to predict SRV in the shale gas reservoirs using the flowback data after fracturing. The stress sensitivity is included as well. The results of SRV are compared with the microseismic data. The following conclusions can be drawn: (1)The flowback data show that the gaswater ratio is Vshaped, i.e., the early descending stage and the late ascending stage in this field, which can be used to estimate the volume of the effective fracture system(2)The stress sensitivity is a key factor affecting the permeability of the effective fracture system as well as the SRV in the shale gas reservoirs. Once it is neglected, the estimation will be overestimated(3)Stress sensitivity is taken into account to forecast SRV in this typical block. Results show that the relative deviation of stimulated reservoir volume calculated by this proposed model and the microseismic data is less than 10%, indicating that this method could provide reasonable prediction
Data Availability
The test data used to support the findings of this study are included within the article. Readers can obtain data supporting the research results from the test data table in the paper.
Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this paper.
Acknowledgments
This research was funded by the National Science and Technology Major Project of China (2017ZX05009005), Fundamental Research Funds for the Central Universities (2652018209), and National Natural Science Foundation of China (51804282).
References
 G. Wang, A. Jia, Y. Wei, and C. Xiao, “Transient pressure analysis for multifractured horizontal well with the use of multilinear flow model in shale gas reservoir,” Geofluids, vol. 2020, Article ID 8348205, 20 pages, 2020. View at: Publisher Site  Google Scholar
 W. Zhou, R. Banerjee, B. D. Poe, J. Spath, and M. Thambynayagam, “Semianalytical production simulation of complex hydraulicfracture networks,” SPE Journal, vol. 19, no. 1, pp. 6–18, 2014. View at: Google Scholar
 C. E. Neuzil, “How permeable are clays and shales?” Water Resources Research, vol. 30, no. 2, pp. 145–150, 1994. View at: Publisher Site  Google Scholar
 R. S. Jones Jr., B. Pownall, and J. Franke, “Estimating reservoir pressure from early flowback data,” in Unconventional Resources Technology Conference, pp. 25–27, Denver, Colorado, August 2014. View at: Google Scholar
 B. Zanganeh, M. Soroush, J. D. WilliamsKovacs, and C. R. Clarkson, “Parameters affecting load recovery and oil breakthrough time after hydraulic fracturing in tight oil wells,” in SPE/CSUR Unconventional Resources Conference, pp. 20–22, Calgary, Alberta, October 2015. View at: Google Scholar
 A. Mohammadnejad, R. F. Shelley, L. V. Lehman, K. Shah, D. Gusain, and M. T. Conway, “Development of the brittle shale fracture network model,” in SPE Hydraulic Fracturing Technology Conference, pp. 4–6, The Woodlands, Texas, USA, February 2013. View at: Google Scholar
 J. Xu, C. Guo, W. Teng, M. Wei, and R. Jiang, “Production performance analysis of tight oil/gas reservoirs considering stimulated reservoir volume using elliptical flow,” Journal of Natural Gas Science and Engineering, vol. 26, pp. 827–839, 2015. View at: Publisher Site  Google Scholar
 M. Wu, M. Ding, J. Yao, S. Xu, L. Li, and X. Li, “Pressure transient analysis of multiple fractured horizontal well in composite shale gas reservoirs by boundary element method,” Journal of Petroleum Science and Engineering, vol. 162, pp. 84–101, 2018. View at: Publisher Site  Google Scholar
 Z. Chen, X. Liao, X. Zhao, W. Yu, and K. Sepehrnoori, “A workflow based on a semianalytical model to estimate the properties of stimulated reservoir volume of tightoil wells,” Journal of Petroleum Science and Engineering, vol. 178, pp. 892–903, 2019. View at: Publisher Site  Google Scholar
 D. Ilk, D. M. Anderson, G. W. J. Stotts, L. Mattar, and T. A. Blasingame, “Productiondata analysis, pitfalls, diagnostics,” SPE Resrervoir Evaluation & Engineering, vol. 13, no. 3, pp. 538–552, 2013. View at: Google Scholar
 D. Ilk, S. M. Currie, D. Symmons, J. A. Rushing, N. J. Broussard, and T. A. Blasingame, “A comprehensive workflow for early analysis and interpretation of flowback data from wells in tight gas/shale reservoir systems,” in SPE Annual Technical Conference and Exhibition, pp. 19–22, Florence, Italy, September 2010. View at: Google Scholar
 M. A. Abbasi, H. Dehghanpour, and R. V. Hawkes, “Flowback analysis for fracture characterization,” in SPE Canadian Unconventional Resources Conference, Calgary, Alberta, November 2012. View at: Google Scholar
 M. A. Abbasi, D. O. Ezulike, H. Dehghanpour, and R. V. Hawkes, “A comparative study of flowback rate and pressure transient behavior in multifractured horizontal wells completed in tight gas and oil reservoirs,” Journal of Natural Gas Science and Engineering, vol. 17, no. 2, pp. 82–93, 2014. View at: Publisher Site  Google Scholar
 Y. Zhang and C. EhligEconomides, “Accounting for remaining injected fracturing fluid in shale gas wells,” in SPE/AAPG/SEG Unconventional Resources Technology Conference, pp. 25–27, Denver, Colorado, USA, August 2014. View at: Google Scholar
 E. Ghanbari, D. Bearinger, M. A. Abbasi, and H. Dehghanpour, “Flowback volumetric and chemical analysis for evaluating load recovery and its impacton on earlytime production,” in SPE Unconventional Resources Conference Canada, pp. 5–7, Calgary, Alberta, Canada, November 2013. View at: Google Scholar
 M. A. Abbasi, A comparative study of flowback rate and pressure transient behavior in multifractured horizontal wells, [Ph.D. thesis], University of Alberta, Edmonton, Alberta, Canada, 2013.
 C. R. Clarkson, “Modeling twophase flowback of multifractured horizontal wells completed in shale,” in SPE Canadian Unconventional Resources Conference, Calgary, Alberta, Canada, November 2012. View at: Google Scholar
 J. D. WilliamsKovacs and C. R. Clarkson, “Modeling of twophase flowback from multifractured horizontal tight gas wells stimulated with nitrogen energized frac fluid,” in SPE Unconventional Resource Conference Canada, Calgary, Alberta, November 2013. View at: Google Scholar
 Z. Tao, L. Xiangfang, Y. Lifeng, L. Jing, W. Yonghui, and F. Dong, “Effects of shutin timing on flowback rate and productivity of shale gas wells,” Natural Gas Industry, vol. 37, no. 8, pp. 48–60, 2017. View at: Google Scholar
 O. A. Adefidipe, Y. Xu, H. Dehghanpour, and C. J. Virues, “Immediate gas production from shale gas wells: a twophase flowback model,” in SPE Unconventional Resources ConferenceUSA, The Woodlands, Texas, April 2014. View at: Google Scholar
 O. A. Adefidipe, Y. Xu, H. Dehghanpour, and C. J. Virues, “Estimating effective fracture volume from earlytime production data: a material balance approach,” in SPE/CSUR Unconventional Resources ConferenceCanada, Calgary, Alberta, October 2014. View at: Google Scholar
 Y. Xu, O. A. Adefidipe, and H. Dehghanpour, “Volumetric analysis of twophase flowback fata for fracture characterization,” in SPE Western Regional Meeting, Garden Grove, California, April 2015. View at: Google Scholar
 Y. Xu, Y. Fu, D. O. Ezulike, H. Dehghanpour, and C. J. Virues, “Modeling twophase flowback data using an open tank model,” in SPE/CSUR Unconventional Resources Conference, Calgary, Alberta, October 2015. View at: Google Scholar
 D. O. Ezulike, H. Dehghanpour, and R. V. Hawkes, “Understanding flowback as a transient twophase displacement process: an extension of the linear dualporosity model,” in CPR Unconventional Resources ConferenceCanada, Calgary, Alberta, November 2013. View at: Google Scholar
 D. O. Ezulike and H. Dehghanpour, “A workflow for flowback data analysis–creating value out of chaos,” in 2014 Unconventional Resources Technology Conference, Denver, Colorado, August 2014. View at: Google Scholar
 D. O. Ezulike, H. Dehghanpour, and C. J. Virues, “A flowbackguided approach for production data analysis on tight reservoirs,” in 2014 SPE/CSUR Unconventional Resources ConferenceCanada, Calgary, Alberta, October 2014. View at: Google Scholar
 C. R. Clarkson and J. D. WilliamsKovacs, “A new method for modeling multiphase flowback of multifractures horizontal tight oil wells to determine hydraulic fracture properties,” in SPE Annual Technical Conference and Exhibition, New Orleans, Louisiana, October 2013. View at: Google Scholar
 C. R. Clarkson, F. Qanbari, and J. D. WilliamsKovacs, “Innovative use of ratetransient analysis methods to obtain hydraulicfracture properties for lowpermeability reservoirs exhibiting multiphase flow,” Leading Edge, vol. 33, no. 10, pp. 1108–1122, 2014. View at: Publisher Site  Google Scholar
 J. D. WilliamsKovacs and C. R. Clarkson, “Stochastic modeling of multiphase flowback from multifractured horizontal tight oil wells,” in SPE Unconventional Resource Conference–Canada, Calgary, Alberta, November 2013. View at: Google Scholar
 J. D. WilliamsKovacs and C. R. Clarkson, “Stochastic modeling of twophase flowback of multifractured horizontal wells to estimate hydraulic fracture properties and forecast production,” in SPE Unconventional Resource Conference USA, The Woodlands, Texas, April 2013. View at: Google Scholar
 J. D. WilliamsKovacs and C. R. Clarkson, “A modified approach for modeling twophase flowback from multifractured horizontal shale gas wells,” in Unconventional Resources Technology Conference, San Antonio, Texas, July 2015. View at: Google Scholar
 J. D. WilliamsKovacs and C. R. Clarkson, “A modified approach for modeling twophase flowback from multifractured horizontal shale gas wells,” Journal of Natural Gas Science and Engineering, vol. 30, pp. 127–147, 2016. View at: Publisher Site  Google Scholar
 Y. Fu, D. O. Ezulike, and H. Dehghanpour, “Estimating effective fracture porevolume from early singlephase flowback data and relating it to fracture design parameters,” in SPE/CSUR Unconventional Resources Conference, Calgary, Alberta, October 2015. View at: Google Scholar
 A. B. Alkouh, S. McKetta, and R. A. Wattenbarger, “Estimation of effectivefracture volume using waterflowback and production data for shalegas wells,” Journal of Canadian Petroleum Technology, vol. 53, no. 5, pp. 293–303, 2014. View at: Google Scholar
 X. Wang, Y. Zhu, and C. Fu, “Experimental investigation of the stressdependent permeability in the Longmaxi Formation shale,” Journal of Petroleum Science and Engineering, vol. 175, pp. 932–947, 2019. View at: Publisher Site  Google Scholar
 X. Yan, J. Sun, and D. Liu, “Numerical simulation of shale gas multiscale seepage mechanismcoupled stress sensitivity,” Journal of Chemistry, vol. 2019, Article ID 7387234, 13 pages, 2019. View at: Publisher Site  Google Scholar
 W. Liu, J. Liu, M. Cai, C. Luo, X. Shi, and J. Zhang, “Pore evolution characteristic of shale in the Longmaxi Formation, Sichuan Basin,” Petroleum Research, vol. 2, no. 4, pp. 291–300, 2017. View at: Publisher Site  Google Scholar
 S. Chen, Y. Zhu, S. Chen, Y. Han, and C. Fu, “Hydrocarbon generation and shale gas accumulation in the Longmaxi Formation, Southern Sichuan Basin, China,” Marine and Petroleum Geology, vol. 86, pp. 248–258, 2017. View at: Google Scholar
 C. Wang, B. Zhang, Y. Lu et al., “Lithofacies distribution characteristics and its controlling factors of shale in Wufeng FormationMember 1 of Longmaxi Formation in the Jiaoshiba area,” Petroleum Research, vol. 3, no. 4, pp. 306–319, 2018. View at: Publisher Site  Google Scholar
 T. Dong, S. He, M. Chen et al., “Quartz types and origins in the paleozoic WufengLongmaxi Formations, Eastern Sichuan Basin, China: implications for porosity preservation in shale reservoirs,” Marine and Petroleum Geology, vol. 106, pp. 62–73, 2019. View at: Publisher Site  Google Scholar
 Y. Cheng, “Impact of water dynamics in fractures on the performance of hydraulically fractured wells in gasshale reservoirs,” Journal of Canadian Petroleum Technology, vol. 51, no. 2, pp. 143–151, 2013. View at: Google Scholar
 J. Zhang, A. Kamenov, A. D. Hill, and D. Zhu, “Laboratory measurement of hydraulicfracture conductivities in the Barnett shale,” SPE Production & Operations, vol. 29, no. 3, pp. 216–227, 2014. View at: Publisher Site  Google Scholar
 L. Fan, J. W. Thompson, and J. R. Robinson, “Understanding gas production mechanism and effectiveness of well stimulated in the Haynesville shale through reservoir simulation,” in SPE Unconventional Resources & International Petroleum Conference, Calgary, Alberta, October 2010. View at: Google Scholar
 F. J. Kuchuk, D. Biryukov, and T. Fitzpatrick, “Rate transient and decline curve analyses for continuously (dualporosity) and discretely naturally fractured reservoirs,” in SPE Annual Technical Conference and Exhibition, Amsterdam, Netherlands, October 2014. View at: Google Scholar
 N. M. A. Rahman, L. Kok, and K. Zaoral, “A new method for computing pseudotime for real gas flow using the material balance equation,” Journal of Canadian Petroleum Technology, vol. 45, no. 10, pp. 36–44, 2006. View at: Google Scholar
 S. Moghadam, O. Jeje, and L. Mattar, “Advanced gas material balance in simplified format,” Journal of Canadian Petroleum Technology, vol. 50, no. 1, pp. 90–98, 2013. View at: Google Scholar
 V. K. Singh, “Overview of material balance equation (MBE) in shale gas & nonconventional reservoir,” in SPE Middle East Oil and Gas Show and Conference, Manama, Bahrain, March 2013. View at: Google Scholar
 Y. Xu, O. A. Adefidipe, and H. Dehghanpour, “A flowing material balance equation for twophase flowback analysis,” Journal of Petroleum Science & Engineering, vol. 142, pp. 170–185, 2016. View at: Publisher Site  Google Scholar
 J. C. Martin, “Simplified equations of flow in gas drive reservoirs and the theoretical foundation of multiphase pressure build up analyses,” Transactions of the AIME, vol. 216, no. 1, pp. 321–323, 2013. View at: Google Scholar
 D. P. Craig, Analytical modeling of a fractureinjection/falloff sequenceand the development of a refracture candidate diagnostic test, [Ph.D. thesis], Texas A&M University, College Station, Texas, 2006.
 D. Zhang and C. Winter, “Momentequation approach to single phase fluid flow in heterogeneous reservoirs,” SPE Journal, vol. 4, no. 2, pp. 118–127, 1999. View at: Google Scholar
 S. H. Tabatabaie, L. Mattar, and M. PooladiDarvish, “Pseudotime calculationin low permeability gas reservoirs,” in SPE Unconventional Resources Conference Canada, Calgary, Alberta, Canada, November 2013. View at: Google Scholar
 D. O. Ezulike and H. Dehghanpour, “Modelling flowback as a transient twophase depletion process,” Journal of Natural Gas Science and Engineering, vol. 19, pp. 258–278, 2014. View at: Publisher Site  Google Scholar
 Y. Zhao, Z. Wang, X. Qin, J. Li, and H. Yang, “Stressdependent permeability of coal fracture networks: a numerical study with lattice Boltzmann method,” Journal of Petroleum Science and Engineering, vol. 173, pp. 1053–1064, 2019. View at: Publisher Site  Google Scholar
 Z. Ru, K. An, and J. Hu, “The impact of sulfur precipitation in developing a sour gas reservoir with pressuresensitive effects,” Advances in GeoEnergy Research, vol. 3, no. 3, pp. 268–276, 2019. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Qi Chen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.