Abstract

Hydraulic fracturing technology has become a key technology for the development of low-permeability/tight oil and gas reservoirs. The evaluation on the postfracturing effect is imperative to the formulation and implementation of the fracturing and development plan. Based on the characteristics of the flow in fracture network after a large-scale hydraulic fracturing, a numerical method for evaluating the effect of fracturing in vertical well was established. This study conducts postfracturing effect evaluations to block C Oilfield’s wells that underwent conventional fracturing and volumetric fracturing, respectively, proposes the definition of fracture network conductivity and its relationship with cumulative production, and analyzes the fracturing construction parameters. The results suggest that the conventional fracturing can only form a single fracture instead of a stimulated reservoir volume (SRV) region. However, the volumetric fracturing transformation can form a complex fracture network system and SRV region and meanwhile bring obvious increase in the production. The effective time lasts for a longer period, and the increase of average daily oil is 2.2 times more than that of conventional fracturing. Additionally, with the progress of the production, the SRV area within the core region of the volume transformation gradually decreased from 6664.84 m2 to 4414.45 m2; the SRV area of the outer region decreased from 7913.5 m2 to 5391.3 m2. As the progress develops, the equivalent permeability and the area of the fracture gradually decrease as the fracturing effect gradually weakens, and so does the conductivity of the network decreasing exponentially; a good correlation is observed between the conductivity of the fracture network, the cumulative production, and fracturing construction parameters, which can serve as the evaluation parameters for the fracturing effects and the basis for fracturing productivity prediction and provide a guidance for fracturing optimization design.

1. Introduction

As the technology of unconventional resource exploration and development advances, tight reservoirs and their production contribution have increased year by year and become a substantial part in China’s current and future crude oil reserves [13]. The development of tight reservoirs has become the key role in the progress of enhancing national energy source security and realizing energy self-sufficiency. The tight oil and gas reservoirs are characterized by low permeability, poor pore throat structure, and no natural production [46]. The conventional development method can barely achieve profitability and thus requires volume fracturing which refers to the generation of fracture network by artificial approaches, improving the reservoir seepage capacity and the contact area between the fracture and the reservoir [711]. The evaluation of production after volume fracturing has always been the focus for petroleum engineers.

Methods evaluating fracturing effects can be distinguished into direct methods and indirect methods. The former one mainly includes microseismic, inclinometer, well temperature logging, and tracer distributed fiber. However, this method is limited by the monitoring method and time conditions and therefore can only detect the fracture effect at a certain point in time and cannot perform long-term monitoring [12]. At the same time, the indirect method does not restrict to the monitoring method and the required data is easy to obtain. Therefore, compared to the direct method, indirect approach is more widely used in the petroleum industry. Perrine proposed a pressure method for multiphase flow test [13]. He pointed out that if the flow coefficient and compression coefficient of a single-phase fluid are replaced by the total fluidity of the multiphase fluid and the comprehensive compression coefficient, the well test formula derived in the case of single-phase flow can also be used for a multiphase flow test well. Martin pointed out that through theoretical analysis, the essence of Perrine’s pressure method is to ignore the saturation gradient in the reservoir [14]. Based on the method’s well test analysis, the integrated flow coefficient, skin factor, and average formation pressure can be obtained. The production index of the well can also be obtained by analogy with the steady flow of the single-phase flow. Serra et al. proposed a method for calculating a relative permeability curve using unstable pressure drop or pressure recovery data [15]. Li et al. proposed a new production data analysis method to address the abrupt change issue based on the virtual equivalent time. The mathematical model with a variable flowrate provides new definitions of normalized pseudopressure and pseudotime with consideration of the pressure-dependent permeability for multifractured horizontal wells in shale gas reservoirs [16]. Hatzignatiou et al. studied the effect of wellbore storage on multiphase flow system pressure [17]. Zhang and Yang adopted a semianalytical technique to examine stress-sensitive effect on the transient pressure responses of a multiple-fractured horizontal well in an unconventional reservoir with an arbitrary shape; sensitivity analysis on pressure responses and their corresponding derivative curves has been performed with respect to the convergence skin effect, boundary shape, maximum distance, and minimum fracture conductivity [18]. Xu et al. adopted a semianalytical pseudo-3-D geomechanical model to describe induced fracture network. The hydraulically stimulated volume is represented by a horizontally expanding ellipse containing a simplified fracture network consisting of two sets of vertical planar fractures perpendicular to one another. This model provides a mathematically equivalent description of the process of hydraulic fracture propagation and the characteristics of the induced fractures [19].

Lei established the concept of numerical well test model for the first time, considering the wellbore storage effect and skin effect [20]. Puchyr first proposed the concept of numerical well test model, which assumes that the fluid is an incompressible single-phase fluid under single well conditions and forms a corresponding reservoir numerical simulator [21]. Archer through the study of numerical well test models found that Green’s function is the best choice for numerical interpretation models [22].

Meyer and Bazan aimed at the formation of complex fracture network structures with high conductivity after volume fracturing in oil wells based on the material balance and momentum conservation equations; the discrete fracture network is proposed. However, when using this model to calculate the fracture network parameters, it is necessary to combine the fracturing construction parameters and the in situ stress parameters, and the seepage field law of the well production stage after fracturing cannot be obtained [23]. Xu et al. implement a discrete fracture model (DFM) and use flexible Delaunay triangulation to represent individual fractures. Numerical simulations are conducted to investigate the influence of rock deformation on a reservoir with complex fracture network [24].

However, for tight reservoirs, it is difficult to derive analytical solutions for seepage equations under complex boundary conditions and for numerical well testing, previous works used simple grids that show certain limitation regardless of whether being structural or unstructured.

In this paper, based on the mathematical model of numerical well testing, we used the hybrid meshing method with radial mesh and constrained PEBI mesh to evaluate the fracturing effect based on the production data. This method was applied to the vertical wells undergoing either conventional fracturing or volumetric fracturing in a block of the C Oilfield, China, and evaluated and compared their effects after fracturing.

2. Model Description

2.1. Physical Model

Volumetric fracturing not only forms the main fractures but also creates a reformed zone composed of secondary fractures. After the fracturing transformation, the fracturing area becomes complicated and needs to be simplified in a numerical solution. Based on the black oil model, the following basic assumptions are proposed:

(1) the fluid flow in the reservoir satisfies the Darcy flow; (2) there are only oil and water phases in the model. The oil and water components are nonmiscible to each other, and there is no mass exchange between the two components; (3) the fluid in the reservoir is in a constant temperature and thermodynamic equilibrium; (4) considering the pressure sensitivity, both the fluid and the rock are compressible; (5) reservoir has heterogeneity and anisotropy; and (6) capillary force is negligible. The reservoir modification volume (SRV) formed by volumetric fracturing transformation is represented by “.” In this paper, the main fractures are viewed as a grid with infinite conductivity; the transformation is distinguished into the core area (main fracture affected area) and the outer area (secondary fracture affected area), and the equivalent permeability of the corresponding area is increased (Figure 1). The porosity and permeability of the core zone are relatively greater, while those of the outer zone are less. The porosity and permeability of the unmodified zone are the basic parameters of the reservoir.

2.2. Mathematical Model

Based on the above assumption, the substitution of Darcy’s law into the continuity equation can be derived from the following governing equation.

The oil phase equation is

The water phase equation is where is the mesh permeability, μm2; and are the relative permeability of oil and water, dimensionless; is viscosity, Pa·s; is the formation volume coefficient; is pressure, Pa; is density, kg/m3; is gravitational acceleration, m/s2; is vertical coordinate; is reservoir porosity; and are oil and water saturation, dimensionless; and are the flowrate of oil and water for ground conditions, respectively, m3/day; subscripts o and w represent oil phase and water phase, respectively.

The above two equations have four variables: , , , and . To solve the equation, two additional equations need to be added. Assume only oil and water exist in the reservoir. where is the capillary force of the oil phase and the water phase and a function of the water saturation and can be obtained from the oil-water phase flow experiment. The two variables in the governing equation can be eliminated by equations (3) and (4), and only two variables and are left to be solved.

Oil production equation

Water production equation where is the distance from the grid adjacent to the fractures, m; is the fracture width, m; is the thickness of the grid, m; is the wellbore pollution coefficient, dimensionless; is the wellbore pressure, Pa; and is the pressure of the grid adjacent to fracture, Pa.

It is obvious from equations (5) and (6) that there is a certain relationship between oil production and water production:

The well is producing at constant total liquid production , where oil production and water production can be obtained via formula (7). Considering the storage effect of the well, we can get the following equation: where is the wellbore storage coefficient, m3/MPa.

Equations (1), (2), (5), and (6), together with closed boundaries and initial conditions, form a well-posed equation group. Normally, Cartesian grids are commonly used in reservoir simulations. However, in reservoirs with complex geological conditions, the boundary shape of the reservoir cannot be accurately described, and each grid is not guaranteed to be effective. The use of PEBI (perpendicular bisection) grids here can overcome the shortcomings of Cartesian grids, which are widely used in the field of numerical simulation, especially in numerical well testing [2527]. By fully implicitly linearizing the equations, the well-posed equations can be numerically solved [28].

3. Production Data Interpretation Method

According to the historical data of development and production, the bottom hole pressure can be calculated at different development times; then, the designated zone model (inner zone, outer zone, and unreconstructed zone) is used to fit and conduct the sensitivity analysis to analyze the primary and secondary fracture penetration parameters such as the permeability, fracture’s half-length, the reformed area, and the bedrock’s permeability. The following information indicates a vertical well whose 45-day postfractured historical production data were used for curve fitting. The following table shows the basic parameters of the reservoir.

Any change in the bottom hole pressure will be observed via tracking bedrock’s permeability, fracture’s half-length, and the main and secondary fracture zone permeability. The process of curve fitting and results are as follows: (1)Establish a model, including well type (horizontal well, fracturing vertical well, and fracturing horizontal well), sedimentary facies distribution, boundary properties, and faults(2)Enter the reservoir and layers’ properties, including the top and bottom boundary types, permeability of porous media, compressibility, field depth, thickness, saturation distribution, and porosity distribution(3)Enter the properties of the well including wellbore coordinates, fracture half-length, fracture orientation, wellbore reservoir coefficient, and skin factor(4)Enter the flowrate and bottom hole pressure data during the curve fitting. The well produces at constant flowrate. The calculated bottom hole pressure will match the measured bottom hole pressure, which is often converted from casing pressure and dynamic liquid level(5)Predict a set of fracturing effect evaluation parameters: the area of the core zone and the outer zone, the permeability, the reservoir permeability distribution, and the fracture half-length. During parameter adjustments, the parameters such as the area and permeability of the core area and the external area need to be adjusted repeatedly(6)Match the calculated bottom hole pressure with the measured bottom hole pressure to obtain the fitting result. Determine whether the fitting result meets the accuracy requirement. If it does, the designated parameter is the fracturing effect evaluation parameter; if not, return to step (5).

The matrix permeability mainly affects the flow at a distant end. When the permeability of the matrix decreases, the flow capacity in a distance is weakened, which also reduces the flow capacity of the matrix to the fractures. Therefore, the permeability of the matrix affects the shape of the entire pressure curve (Figure 2(a)). The half-length of the fractures mainly affects the front part’s shape of the pressure curve, and the latter half is parallel to each other; the fractures only affect the shape of the pressure curve in the early stage, because the fracture has a short half-length and the well’s productivity index decreases, causing rapid drop in pressure (Figure 2(b)). The foremost section of the pressure curve almost overlaps, indicating that the inner zone permeability has a small effect on the pressure relative to the main fractures. After this section, the pressure curve begins to be nearly parallel, which means that the flow characteristics at the later stage cause no significant impact, and the main fracture permeability mainly affects the flow in well’s vicinity (Figure 2(c)). The secondary fracture zone has less impact on the early pressure curve, indicating that the second outer zone mainly affects the medium-term pressure (Figure 2(d)).

4. Application and Evaluation

4.1. Fracturing Effect Evaluation

The method is applied to three typical vertical wells adjacent to the same block in the C Oilfield. Due to low production and lack of theoretical understanding, these three wells underwent conventional fracturing and reforming measures in the early stage, but soon, the stimulation effect disappeared, and production decreased. In recent years, the theory of “large-scale hydraulic fracturing” has been put forward and has been successful in field tests. Therefore, these three wells were subjected to volumetric hydraulic fracturing at the same time. The basic parameters of the reservoir are shown in Table 1. The fracturing parameters are shown in Table 2.

The C1 well was fractured and put into operation in August 2007. The initial daily production liquid was 4.12 m3, the daily oil production was 2.93 t, the water content was 14.1%, and the hydrodynamic surface was 1494 m. In April 2013, the liquid production decreased. Since then, the daily output was 0.89 m3, the daily oil production was 0.81 t, and the water content was 9.38%. After the temporary plugging and fracturing in September 2013, the daily liquid production was 2.75 m3, the oil production was 1.35 t, and the output increased. The water cut is stable. The fracturing production data and the bottomhole flow pressure data are collated to obtain Figure 3.

The fracturing production data and the converted bottom hole pressure data are compared, where the production data is relatively stable and the pressure is fitted by the interval of the adjacent well interference. A total of two sections are fitted to observe the fracturing effect dynamic changes.

The first period of time is the production data for 24 days from September 19 to October 12, 2013. After curve fitting, the permeability around the fractures was 10 mD and the permeability of the unmodified region was 0.3 mD. The fracture half-yield is 103 m, and the average liquid production is 2.33 m3, which is 85% of the initial level. The figure below shows the permeability distribution.

It can be seen from Figure 4(a) that the SRV is not formed during the fracturing, and only a single main fracture is formed. Figure 4(b) is a fitted pressure history diagram.

As shown in Figure 5(a), when the flowback of fracturing fluid finished and the well was producing at early stages, the pressure in the well’s vicinity was as high as 7.76 MPa due to the existence of crude oil around it. As the well kept producing, the bottom hole pressure field caused by the production gradually spread to the reservoir’s boundary, making the pressure in the vicinity of well and fracture decline while the boundary pressure gradually increases. By the 5th day (Figure 5(b)), the pressure field’s effect had spread to the boundary, but the boundary pressure had not reached its maximum value, and the pressure drop around the well and the fractures remained small. As the production reaches the eighth day (Figure 5(c)), the pressure effect continued to expand toward the boundary while the influence of the fracturing fluid began to disappear and formed obvious contrast to pressure drop around well and fractures. When the production went into the 19th day (as shown in Figure 5(d)), the pressure around the wellbore and the fractures was smaller than the one at the boundary, forming a depressurization cone produced by pressure drawdown.

The second period of time was 38 days of production data ranging from December 17, 2013, to January 23, 2014. By curve fitting the fracture half-length to 62 m, the effective permeability reduced to 2 mD, and the average daily liquid production reduced to 0.96 m3, which was 35% of the initial value. The stimulation effect from fracturing reduced quickly. Hence, it is obvious that the half-length of the fractures has a great impact on the production, and as the fractures gradually close, the fracture conductivity decreases and so does and the production.

The remaining time periods of the three wells went through curve fitting according to the same method mentioned above, and the fitting results are shown in Table 3.

It can be seen from Table 3 that Mitsui’s large-scale hydraulic fracturing is more effective than conventional fracturing in terms of stimulation results. Conventional fracturing in Well C1 could only form a single fracture due to the small volume of liquid and proppant. As the development progresses, fracture gradually closed the equivalent permeability and gradually reduced from 10 mD to 6 mD; the half-length of the fracture rapidly declined with decreased production. At last, as the fracture closed, the equivalent permeability dropped to 0.8 mD. The stimulation effect caused by fracturing basically disappeared, and the production dropped to the without-fracturing level. In the stage of second large-scale hydraulic fracturing, due to the significant volume of liquid and proppant on the ground and a small ratio, a large-scale volume reconstruction area was formed underground. The fracture half-length explained was 110 m, and the permeability of the core area was 15 mD. The area was 5747.6 m2; the permeability of the external reconstruction area was 1 mD, and the reconstruction area was 8471.8 m2. The amount of fluid entering the ground during conventional fracturing in Wells C3 and C4 has increased, forming a smaller fracture network. As the process developed, the fractures gradually closed and the stimulated effect disappeared quickly. In large-scale hydraulic fracturing process, the amount of liquid and the volume of proppant entering formation of Well C4 were much higher than those of Well C3, while the permeability and reformed area of the inner zone and the permeability of the outer zone were slightly better than those of Well C3. It can be observed that the amount of liquid and proppant entering formation, when reaching a certain level, would actually reduce the stimulation effect.

4.2. The Calculation Method and Influence Factor Analysis of the Stimulated Reservoir Conductivity (SRC)

The main purpose of large-scale hydraulic fracturing is to form seepage channels with high conductivity within the stimulated reservoir, thereby reducing the resistance of oil and gas flowing from the formation to the wellbore. In terms of measuring stimulated reservoir’s conductivity, it has been widely recognized that the single fracture conductivity characterization method is not suitable for fracture network conductivity. Therefore, here, the product of the equivalent permeability of the fracture network and the corresponding area was used. Its definition is expressed as follows: where is the area of the stimulated reservoir, m2, and is the corresponding equivalent permeability, mD.

It can be seen from Figure 6 that the SRC shows an exponential decrease with time, but the stimulated reservoir conductivity in the external area (SRCE) decreases faster. This characteristic is consistent with the declining law of the conductivity of a single fracture, which also exemplifies the SRC formula proposed in this paper. The external area was mainly composed of microfractures produced by fracturing. Along with the development process, these microfractures were more vulnerable to formation stress sensitivity and therefore could close more easily with severer conductivity decline.

Figure 7 is a diagram of the relationship between the stimulated reservoir conductivity of the internal area (SRCI) and SRCE with 1.5 years of cumulative production after fracturing. According to the diagram, the cumulative production has a strong correlation with the fracture network conductivity of the inner and outer zones, where the correlation coefficients are 0.83 and 0.88, respectively. This result has proved that the fracture network conductivity defined by parameters inversely computed from this paper’s numerical method is capable of reflecting the large-scale hydraulic fracturing effect and can serve in accordance with oil and gas production forecast.

In order to explore the influence caused by fracturing operating parameters on the SRC, the correlation of the first stage of SRCI and SRCE with the breakdown pressure, injection volume, proppant volume, proppant concentration, reservoir thickness, and fracture half-length was analyzed. The analysis results are shown in Figure 8.

It can be seen from Figure 8 that the first stage SRCI has no correlation with breakdown pressure and reservoir thickness but has a weak correlation with proppant volume and a strong negative correlation with proppant concentration with a correlation coefficient of 0.95; strong positive correlations with injection volume and fracture half-length are observed, where the correlation coefficients are 0.88 and 0.99, respectively. The first stage SRCE has no correlation with the thickness of the oil layer but has a weak correlation with the fracture pressure, proppant volume, and fracture half-length and a strong positive correlation with the injected volume, with a correlation coefficient of 0.83; a strong negative correlation with sand ratio exists, and the correlation coefficient is 0.93.

Through the above analysis, there is a good correlation between the SRC and the cumulative production, which reflects an important parameter for the effect of large-scale hydraulic fracturing. The SRC is also correlated to the injection volume, proppant volume, and proppant concentration, and therefore, it can provide guidance for the optimization of large-scale hydraulic fracturing; at the same time, it can also influence the geological characteristics of the reservoir according to the productivity of oil and gas wells and the effect of fracturing engineering. Reevaluation would improve the spatial distribution characteristics of reservoirs and the accuracy of sweet spot prediction. In the actual application process, for oil and gas wells in the same area, the relationship between the evaluation results of large-scale hydraulic fracturing engineering reconstruction effects and productivity is counted first. Then, through actual examples, the productivity of oil and gas wells after large-scale hydraulic fracturing is compared and analyzed; the stimulation effect is also affected by geological conditions such as reservoir porosity, permeability, and oil abundance. Finally, the factors that lead to productivity differences between wells are determined, such that the factors that influence the well productivity can be distinguished and reservoir’s spatial distribution can be confirmed, where the evaluation of characteristics and formation body location is more accurate.

5. Conclusion

(1)The effect of large-scale hydraulic fracturing is better than that of conventional fracturing. The fracture network area and equivalent permeability of the internal and external regions formed by the transformation are significantly better than conventional fracturing does; after large-scale hydraulic fracturing, the effect of fracturing reformation in different production stages is different, and the permeability and the area of fracturing reformation area gradually decrease with time(2)The conductivity of the internal and external fracture network shows an exponential decrease with time, and the conductivity of the external fracture network decreases faster, which is consistent with the decrease of the conductivity of a single fracture, exemplifying and verifying the formula proposed in this paper for the fracture network’s conductivity(3)There is a good correlation between fracture network conductivity and cumulative production, which can be used to reflect the effect of large-scale hydraulic fracturing and can be used as the basis for production prediction of oil and gas wells after fracturing. At the same time, a good correlation also exists between SRC and the injection volume, proppant volume, and proppant concentration, which can provide guidance for the optimization of large-scale hydraulic fracturing design and influence the spatial distribution characteristics and sweet spot of the reservoir, providing a basis for judgment

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declared that they have no conflicts of interest in this work.

Acknowledgments

We gratefully acknowledge the financial supports from the National Science and Technology Major Project (2017ZX05013-001).