# Production Data Analysis and Practical Applications in the Sulige Tight Gas Reservoir, Ordos Basin, China

Authors: Minhua Cheng; Wen Xue; Meng Zhao; Guoting Wang; Bo Ning; Guang Ji; Lihua Cheng; Yingying Tian
Journal: Geofluids
 (2021)

---

## Abstract

Successful exploitation of tight sandstone gas is one of the important means to ensure the “increasing reserves and production” of the oil and gas initiative and also one of the important ways to ensure national energy security. To further improve the accuracy of historical matching of field data such as gas production and bottom-hole pressure during the production process of this type of gas reservoir, in this study, a new expression of wellbore pressure for the uniform flow of vertical fractured wells in Laplace space based on the point sink function model of vertical fractures in tight sandstone gas reservoirs is constructed. This innovation is based on a typical production data analysis plot of the Blasingame type that uses the numerical inversion decoupling mathematical equation. After analyzing the pressure and pressure derivative characteristics of each flow stage in the typical curves, a new technique of type-curve matching was proposed. In order to verify the correctness of the model and the application value of the field, based on the previous production data of Sulige Gas Field in China, a new set of production data diagnostic chart of tight sandstone gas reservoir was formed. A case analysis showed that the application of the production data analysis method and data diagnosis plot in the field accurately evaluated the development effect of the tight sandstone gas reservoirs, clarified the scale of effective sand bodies, and provided technical support for optimizing and improving the well pattern and realizing the efficient development of gas fields.


---

## Body

## 1. Introduction

Tight gas is widely distributed in China, and the favorable area is32×104 km2 and is primarily distributed in Ordos, Sichuan, Tarim, Songliao, and other sedimentary basins. These areas have become important fields for oil and gas production and storage [1, 2]. According to the fourth China oil and gas resources census, China’s tight gas geological resources reached 21.86×1012 m3, and the recoverable resources were 12.14×1012 m3 [3]. The Sulige gas field is a typical representative of a tight sandstone gas reservoir in China. The primary production layers are the Upper Paleozoic Shihezi Formation and the Shanxi Formation. These have the characteristics of many gas-bearing strata, a thin single-layer thickness, poor reservoir physical properties, strong heterogeneity, a low formation pressure coefficient, and a low reserve abundance [4]. To achieve the scale-benefit development of tight gas reservoirs, it is necessary to implement fracturing measures for reservoirs to minimize the gas flow distance and reduce the flow resistance, so as to maximize the effective scale utilization. With the deepening of development, the well type put into production in Sulige gas field has gradually changed from a single-well to a multiwell cluster deployment. However, it is still dominated by vertical well development, accounting for more than 80% of the total number of wells [5]. The reliability of the productivity evaluation results of the fractured vertical wells directly affects the stable production, development technology optimization, and surface engineering construction of the Sulige gas field.The traditional production data analysis methods have a rich history experience. For example, the early Arps production decline analysis method has been widely used as the primary means to analyze production declines due to its simple application and use of less actual development data. However, gas wells in tight gas reservoirs often have the characteristics of unstable working systems and large differences in production characteristics at different production stages. These led to many unfavorable conditions in the application of the Arps method, such as data dispersion, more singular points, a low fitting accuracy of the test data, and difficulty in identifying the correlation coefficient close to the decline type under the different types of decline. The modern production decline analysis method is based on the theory of unstable flow and combines the advantages of the reservoir engineering method and the modern well test analysis to establish a new type of simulation chart. It uses the daily production dynamic data to quantitatively analyze the seepage characteristics of oil and gas wells, determine reservoir parameters, and calculate well-controlled reserves. This method primarily includes the empirical method represented by the Arps method [6], the classical analysis method represented by the Fetkovich method [7], the modern double logarithmic curve fitting method represented by Blasingame et al.’s method [8], and the reservoir engineering method represented by the FMB method [9]. This method and the entire life cycle of the modern well test analysis have become the primary methods of reservoir dynamic description, and the corresponding analysis software has been widely used in various oilfields at home and abroad. The decline analysis models of vertical wells, vertical fractured wells and horizontal wells in composite, and double permeable and triple media have been established successively. But some researchers have confused the concepts of Blasingame et al.’s material balance (pseudo) time with Fetkovich’s dimensionless time, causing errors in the typical graphs. In addition, the production integral curve and the production integral derivative curve of the vertical fractured wells in tight gas reservoirs and the fitting calculation of effective permeability and the effective fracture half-length are rarely provided in detail.Therefore, this study redefines the material balance (pseudo) time function, transforms the dimensionless pressure reciprocal under the condition of bottom-hole production into the dimensionless production under the condition of bottom-hole pressure, and discusses the productivity dynamics of vertical fractured wells under the condition of unstable flow in a closed formation. After integrating the point source function of the vertical fractured wells, the typical curves of the Blasingame-type production data analysis are obtained using the Laplace space change and numerical inversion calculation, and a new production data analysis method for vertical fractured wells in low permeability tight gas reservoirs is formed. Combined with a large number of field production data research work, by the fitting and regression of actual production data under different reservoir conditions, the problems caused by poor quality of data during each flow stage of the gas wells are summarized. A set of data diagnosis plot suitable for tight sandstone gas reservoirs is formed so as to evaluate the development effect of a vertical fracture single well and a gas zone more accurately.


## 2. Production Data Analysis

### 2.1. Material Balance (Pseudo) Time

For gas wells, when pseudo-steady-state seepage occurs in closed reservoirs, the variable flow problem can be transformed into a constant flow problem by solving the material balance (pseudo) time [10, 11]. This method not only can deal with a change in the natural gas physical properties with pressure but also clearly show the flow characteristics of the different production stages [12]. When the formation temperature is constant, the gas flow formula can be obtained using the gas state equation, the material balance equation, and the gas compression coefficient [13]:
(1)qg=−GZipipavgcgavgZavgdpavgdt.After the deformation of (1), the variable integral is separated, and the two sides are multiplied by μgicti/qgt to obtain:
(2)μgictiqgt∫0tqgμgavgcgavgdt=−GctiqgtμgiZipi∫pipavgpavgμgavgZavgdpavg.Assumingcti=cgavg, the gas-material equilibrium (pseudo) time can be obtained from Russell et al.’s pseudo pressure [14] and Meunier et al.’s pseudo time [15]:
(3)ta=μgictiqgt∫0tqgτμgavgcgavgdτ=Gctiqgtmpi−mpavg,(4)mpi−mpwfqgt=1Gctita+μgiBgi2πkghln4πreD2eγCA=1Gctita+bpss.The results of the liquid percolation can be directly used in the gas percolation problem. After the deformation of Equation (3), the gas variable flow inflow performance relationships in the same form as the liquid percolation proposed in references [16, 17] can be obtained:


### 2.2. Mathematical Production Model of the Fractured Well

As shown in the model in Figure1, a vertical fractured well with a fracture half-length (xf) was located in the center of the homogeneous formation with a drainage radius (re). The fluid was assumed to be single-phase compressible, regardless of the capillary force and gravity, and when isothermal Darcy percolation occurred, the production was constant.Figure 1
Vertical fracture well model of the circular confined reservoir.Without considering the pressure loss of fluid in the fracture, the fracture can be assumed as a straight line sink without considering the pressure loss. The point source in the straight line sink was abstracted as a center. The dimensionless seepage control equation of the point source is as follows [12, 18]
(5)∂2pD∂rD2+1rD∂pD∂rD=∂pD∂tD,pDrD,0=0,∂pDreD,tD∂rD=0,limrD⟶0rD∂pD∂rD=−1.By redefining the dimensionless parameters for the vertical fractured wells in tight gas reservoirs, the following was obtained:(6)pD=2πkghmpi−mpwfqgμgiBgi,tD=kgtaφμgictixf2,tD=kgtaφμgictixf2,rD=rxf,reD=rexf,xD=xxf,yD=yxf.The expression of the wellbore pressure was obtained using the Laplace transform:(7)p~D=1sK1reDsI1reDsI0rDs+1sK0rDs.The study of the vertical fractured wells in conventional reservoirs showed that the wellbore pressure curves of the uniform flow and infinite conductivity vertical fractures were very similar, and the former was smoother than the latter [19, 20]. The bottom-hole pressure formula of the infinite-conductivity fracture well in the Laplace space was obtained using the pressure superposition of Equation (7):
(8)p~DxD,yD,αD,s=12s∫−11K1reDsI1reDsI0xD−αD2+yD2s+K0xD−αD2+yD2sdαD.For cracks with uniform flow rates, the usual treatment method is to integrally average the bottom-hole pressure equation of the infinite-conductivity fracture well, namely, Equation (8) [21]. The specific method is to let yD=0 and integrate xD along the crack and obtain the following:
(9)p~wDxD,s=1ssπ2−Ki12s+K12s−12s+K1reDsI1reDsIi12s−I12s.Equation (9) is the dimensionless bottom-hole pressure expression of the well with a uniform flow vertical fracture in the Laplace space.The Blasingame-type production data analysis curve of the dimensionless production rate can be obtained using the inversion of the dimensionless bottom-hole pressure reciprocal of the vertical fractured well with uniform flow. Then, the typical plot can be formed. The rate of the constant bottom-hole pressure and the pressure of the constant rate in the Laplace space satisfies the following:(10)p~wDsq~wDs=1s2.According to the study, the dimensionless production under the condition of the bottom-hole constant pressure under the material balance (pseudo) time was consistent with the reciprocal expression of the dimensionless wellbore pressure under the condition of the bottom-hole constant production [22]. The production equation obtained from Equations (9) and (10) is as follows:
(11)qwDtaD=1pwDtaD=1L−1p~wDs.The numerical inversion calculation of Equation (11) was conducted using the method in Reference [23], and the solution of Equation (11) in the Laplace relative space was obtained. The above solutions diverged in the pseudo-steady-state stage. In order to normalize the curves in the pseudo -steady-state stage, the dimensionless production rate, qwD of Equation (11) was revised to qdD, and the dimensionless time, tD, was revised to tdD. The specific results are as follows:
(12)qdD=lnreD2+34qwD,(13)tdD=2reD2lnreD2+34−1taD.Substituting Equations (12) and (13) into Equation (11), the dimensionless production decline curve of the Blasingame type was obtained (Figure 2).Figure 2
The log-log plot of the dimensionless production versus the dimensionless material balance pseudo time.As shown in Figure2, the decline curve of the Blasingame was a set of curves with different dimensionless reservoir drainage radius (reD) values during the early unsteady flow process. When reaching the boundary of the closed formation, this set of curves controlled by the dimensionless reservoir drainage radius converges into a harmonic decline curve. The decline curve can be divided into three flow stages: the AB section is an unstable flow stage in the early stage of the well opening. The production data of this stage showed that the fluid was not affected by the boundary, which was equivalent to flowing in an infinite formation. The boundary had no effect on the flow. The dimensionless production was primarily controlled by the dimensionless reservoir drainage radius. When the dimensionless reservoir drainage radius increases gradually, the dimensionless production decline curve will move downward. According to the curve of this period, the relevant information of near well zone can be obtained, such as the effective permeability and skin factor. There were few corresponding production data points in the BC section, primarily due to the short flow period, which was the transition part from the unstable flow stage to the pseudo-steady flow stage. The CD section is a pseudo-steady-state flow stage. The pressure disturbance completely spreads to the boundary; that is, the pressure drop propagates to the boundary and has an impact on the flow. The production of the CD section is completely derived from the expansion of rock and fluid due to pressure drop, and the pressure drop velocity at each point was the same. According to the analysis of this stage, the required parameters, such as the well control dynamic reserves and the reservoir drainage radius, can be calculated and deduced.


### 2.3. Production Integral and Integral Derivative Curve

The dimensionless production decline curve plot can make the derivative curve smoother and easier to judge through the application of the derivative form after the integration of the dimensionless production. As can be seen in the actual derivation process, since the output integral results are sensitive to the error of the early data points, a small error of the early data points will make the derivative curve produce a relatively large cumulative error. According to the field data and the data obtained in the laboratory, the obtained plot can be used to fit the actual data. Reservoir-related physical properties can be obtained by the dimensionless production decline curve, the dimensionless production integral curve, and the derivative curve fitting.When the bottom-hole pressure is known, the integral average equation satisfies the following:(14)qidDtD=1tD∫0tDi+11pwDvdv,(15)It=1t∫0tftdt.Based on the time corresponding to the timest1,t2,⋯,tN, Equation (15) is discretely decomposed and is applied to the trapezoidal method. The resulting form satisfies the following equation:
(16)I2tN=∫t1tNftdt=∑i=1N−1∫titi+1ftdt.A linear interpolation on the interval (ti,fi) and (ti+1,fi+1) is able to obtain the Blasingame-type dimensionless production rate integral curve and the dimensionless production rate integral curve of the derivation (Figure 3).Figure 3
The log-log plot dimensionless production integral and derivative versus dimensionless material balance pseudo time.


## 2.1. Material Balance (Pseudo) Time

For gas wells, when pseudo-steady-state seepage occurs in closed reservoirs, the variable flow problem can be transformed into a constant flow problem by solving the material balance (pseudo) time [10, 11]. This method not only can deal with a change in the natural gas physical properties with pressure but also clearly show the flow characteristics of the different production stages [12]. When the formation temperature is constant, the gas flow formula can be obtained using the gas state equation, the material balance equation, and the gas compression coefficient [13]:
(1)qg=−GZipipavgcgavgZavgdpavgdt.After the deformation of (1), the variable integral is separated, and the two sides are multiplied by μgicti/qgt to obtain:
(2)μgictiqgt∫0tqgμgavgcgavgdt=−GctiqgtμgiZipi∫pipavgpavgμgavgZavgdpavg.Assumingcti=cgavg, the gas-material equilibrium (pseudo) time can be obtained from Russell et al.’s pseudo pressure [14] and Meunier et al.’s pseudo time [15]:
(3)ta=μgictiqgt∫0tqgτμgavgcgavgdτ=Gctiqgtmpi−mpavg,(4)mpi−mpwfqgt=1Gctita+μgiBgi2πkghln4πreD2eγCA=1Gctita+bpss.The results of the liquid percolation can be directly used in the gas percolation problem. After the deformation of Equation (3), the gas variable flow inflow performance relationships in the same form as the liquid percolation proposed in references [16, 17] can be obtained:


## 2.2. Mathematical Production Model of the Fractured Well

As shown in the model in Figure1, a vertical fractured well with a fracture half-length (xf) was located in the center of the homogeneous formation with a drainage radius (re). The fluid was assumed to be single-phase compressible, regardless of the capillary force and gravity, and when isothermal Darcy percolation occurred, the production was constant.Figure 1
Vertical fracture well model of the circular confined reservoir.Without considering the pressure loss of fluid in the fracture, the fracture can be assumed as a straight line sink without considering the pressure loss. The point source in the straight line sink was abstracted as a center. The dimensionless seepage control equation of the point source is as follows [12, 18]
(5)∂2pD∂rD2+1rD∂pD∂rD=∂pD∂tD,pDrD,0=0,∂pDreD,tD∂rD=0,limrD⟶0rD∂pD∂rD=−1.By redefining the dimensionless parameters for the vertical fractured wells in tight gas reservoirs, the following was obtained:(6)pD=2πkghmpi−mpwfqgμgiBgi,tD=kgtaφμgictixf2,tD=kgtaφμgictixf2,rD=rxf,reD=rexf,xD=xxf,yD=yxf.The expression of the wellbore pressure was obtained using the Laplace transform:(7)p~D=1sK1reDsI1reDsI0rDs+1sK0rDs.The study of the vertical fractured wells in conventional reservoirs showed that the wellbore pressure curves of the uniform flow and infinite conductivity vertical fractures were very similar, and the former was smoother than the latter [19, 20]. The bottom-hole pressure formula of the infinite-conductivity fracture well in the Laplace space was obtained using the pressure superposition of Equation (7):
(8)p~DxD,yD,αD,s=12s∫−11K1reDsI1reDsI0xD−αD2+yD2s+K0xD−αD2+yD2sdαD.For cracks with uniform flow rates, the usual treatment method is to integrally average the bottom-hole pressure equation of the infinite-conductivity fracture well, namely, Equation (8) [21]. The specific method is to let yD=0 and integrate xD along the crack and obtain the following:
(9)p~wDxD,s=1ssπ2−Ki12s+K12s−12s+K1reDsI1reDsIi12s−I12s.Equation (9) is the dimensionless bottom-hole pressure expression of the well with a uniform flow vertical fracture in the Laplace space.The Blasingame-type production data analysis curve of the dimensionless production rate can be obtained using the inversion of the dimensionless bottom-hole pressure reciprocal of the vertical fractured well with uniform flow. Then, the typical plot can be formed. The rate of the constant bottom-hole pressure and the pressure of the constant rate in the Laplace space satisfies the following:(10)p~wDsq~wDs=1s2.According to the study, the dimensionless production under the condition of the bottom-hole constant pressure under the material balance (pseudo) time was consistent with the reciprocal expression of the dimensionless wellbore pressure under the condition of the bottom-hole constant production [22]. The production equation obtained from Equations (9) and (10) is as follows:
(11)qwDtaD=1pwDtaD=1L−1p~wDs.The numerical inversion calculation of Equation (11) was conducted using the method in Reference [23], and the solution of Equation (11) in the Laplace relative space was obtained. The above solutions diverged in the pseudo-steady-state stage. In order to normalize the curves in the pseudo -steady-state stage, the dimensionless production rate, qwD of Equation (11) was revised to qdD, and the dimensionless time, tD, was revised to tdD. The specific results are as follows:
(12)qdD=lnreD2+34qwD,(13)tdD=2reD2lnreD2+34−1taD.Substituting Equations (12) and (13) into Equation (11), the dimensionless production decline curve of the Blasingame type was obtained (Figure 2).Figure 2
The log-log plot of the dimensionless production versus the dimensionless material balance pseudo time.As shown in Figure2, the decline curve of the Blasingame was a set of curves with different dimensionless reservoir drainage radius (reD) values during the early unsteady flow process. When reaching the boundary of the closed formation, this set of curves controlled by the dimensionless reservoir drainage radius converges into a harmonic decline curve. The decline curve can be divided into three flow stages: the AB section is an unstable flow stage in the early stage of the well opening. The production data of this stage showed that the fluid was not affected by the boundary, which was equivalent to flowing in an infinite formation. The boundary had no effect on the flow. The dimensionless production was primarily controlled by the dimensionless reservoir drainage radius. When the dimensionless reservoir drainage radius increases gradually, the dimensionless production decline curve will move downward. According to the curve of this period, the relevant information of near well zone can be obtained, such as the effective permeability and skin factor. There were few corresponding production data points in the BC section, primarily due to the short flow period, which was the transition part from the unstable flow stage to the pseudo-steady flow stage. The CD section is a pseudo-steady-state flow stage. The pressure disturbance completely spreads to the boundary; that is, the pressure drop propagates to the boundary and has an impact on the flow. The production of the CD section is completely derived from the expansion of rock and fluid due to pressure drop, and the pressure drop velocity at each point was the same. According to the analysis of this stage, the required parameters, such as the well control dynamic reserves and the reservoir drainage radius, can be calculated and deduced.


## 2.3. Production Integral and Integral Derivative Curve

The dimensionless production decline curve plot can make the derivative curve smoother and easier to judge through the application of the derivative form after the integration of the dimensionless production. As can be seen in the actual derivation process, since the output integral results are sensitive to the error of the early data points, a small error of the early data points will make the derivative curve produce a relatively large cumulative error. According to the field data and the data obtained in the laboratory, the obtained plot can be used to fit the actual data. Reservoir-related physical properties can be obtained by the dimensionless production decline curve, the dimensionless production integral curve, and the derivative curve fitting.When the bottom-hole pressure is known, the integral average equation satisfies the following:(14)qidDtD=1tD∫0tDi+11pwDvdv,(15)It=1t∫0tftdt.Based on the time corresponding to the timest1,t2,⋯,tN, Equation (15) is discretely decomposed and is applied to the trapezoidal method. The resulting form satisfies the following equation:
(16)I2tN=∫t1tNftdt=∑i=1N−1∫titi+1ftdt.A linear interpolation on the interval (ti,fi) and (ti+1,fi+1) is able to obtain the Blasingame-type dimensionless production rate integral curve and the dimensionless production rate integral curve of the derivation (Figure 3).Figure 3
The log-log plot dimensionless production integral and derivative versus dimensionless material balance pseudo time.


## 3. Production Data Diagnosis

It was found in the field research that the quality of production data was crucial for the analysis results. Some production data will have some “illusions,” and these “illusions” may be due to data acquisition instrument problems or improper practical operation methods, which are inevitable [24]. These “illusions” should be removed when processing the data. In common cases, such as the poor quality of the wellhead pressure into bottom-hole pressure, the early production data analysis will be affected by the so-called well washing operation.The number of vertical fractured vertical wells in the Sulige gas field is huge, and because of the gas well’s own factors or external reasons, the actual production data curve fluctuates more. Hence, it is difficult to intuitively see the quality of the gas well production data. A set of typical curves is an effective diagnostic tool. By noise reduction processing of production data, the processed production data curve can be compared with the typical curve chart, and the real interpretation of data can be quickly and directly sought. Based on a large number of production data analyses, this study comprehensively summarizes several common problems encountered in the comparison of production data of vertical fractured gas wells in the Sulige gas field using typical curve charts (Figure4). In addition, a set of diagnostic plots for production data of the tight sandstone gas reservoirs are established (Figure 5 and Table 1). The diagnostic plot can help researchers accurately find problems in the data analysis, improve the data quality, and accurately predict the future of gas wells, the reservoir information, and the characteristic parameters of gas wells.Figure 4
Common problems of the match data of the type curve (data was obtained from the Sulige gas field).Figure 5
Production data diagnosis plate of the tight sandstone gas reservoir.Table 1
Common problems in the production data diagnosis of the tight sandstone gas reservoir.
ProblemCauseI Unsteady state① Flat early dataContaminated well② High early dataInitial reservoir pressure is low③ Normalized rate increases with timeWell washing④ Step-shaped early dataInitial reservoir pressure is high? Large tubing sizeII Transitional stage⑤ Early arrival of boundary flowBottom hole overstock liquidIII Pseudo-steady state⑥ Data below the 45-degree slope lineBottom hole overstock liquidPressure loss caused by wellbore instabilityMultiwell interference (Energy leakage)⑦ Not unifiedUnstable production (Workover, etc.)In the actual field application, the visual data curve and the typical curves are fitted and analyzed from three different flow stages. If the curve trend of the data appears similar to the curve in the diagnostic plot in this study, that is, this part of the data is different from the trend of the typical curve chart, then the engineer can be reminded that the data of this portion has quality problems. Compared with the problems in different flow stages, as shown in Table1, the causes of the problems are discovered, and the corresponding solutions are selected to improve the data quality and improve the accuracy of the gas well analysis.The use of the data diagnosis plot of the tight sandstone gas reservoir to detect the quality of dynamic data from a gas well is not only intuitive, but also convenient. Prior to the evaluation of a single well, the data diagnostic plot can be compared, the problems can be discovered in time, and the reasons can be discovered to ensure the accuracy of the gas well evaluation. Prior to the analysis and evaluation of the block, gas wells with good production data or the gas wells with improved data quality can be selected for priority analysis to help engineers grasp the scale of effective sand bodies and develop the overall block more efficiently.


## 4. Applications in the Field

### 4.1. Single-Well Application

Well V0 in the Sulige gas field is a tight sandstone gas well that was produced after fracturing (artificial fracture 92 m at depth 637 m) in November 2008. The target layer is the reservoir of the Shan 1 member that is primarily composed of grey white-grey medium and coarse sandstone and grey-black fine sandstone. The debris content in the debris particles is high, and the rock types are primarily lithic quartz sandstone and quartz sandstone. The porosity of the reservoir is primarily 4.0–10.0%, with an average of 6.1%, and the permeability is primarily 0.100–1.000 mD, with an average of 0.265 mD. This is typical for a tight sandstone reservoir. In order to avoid a long time shut-in pressure recovery test, the modern production data analysis of well V0 was studied. The gas reservoir description data, production dynamic curve, and fitting curve are shown in Table2 and Figures 6 and 7, respectively.Table 2
Reservoir description parameters of well V0.
ParameterValueParameterValueNet pay (m)15Geological reserve (104 m3)3727.73Shaft radius (m)0.067Total porosity (%)9.5Reservoir temperature (K)383.15Initial viscosity of gas (mPa·s)0.0036Initial water saturation (%)68Initial volume factor0.0021Figure 6
Casing pressure and daily production curve of gas well V0.Figure 7
Typical plot matching result of well V0.The biggest difference between the production data analysis method of the vertical fractured wells in tight gas reservoirs and ordinary vertical wells is that the fracture half-length analysis is increased, and the interpretation results of this study also found this increased skin factor.The specific fitting and research steps were as follows:(i)
The material balance (pseudo) times at each production time point were obtained(ii)
Moving the processed production data points on the Blasingame-type dimensionless yields integral and integral derivative chart. The yield integral curve was first fitted because the curve was smoother than the yield decline curve. Then, the production decline curve and the integral derivative curve were used for the fitting. Finally, the time step and dimensionless discharge radius were adjusted to select the best fitting effect to obtain the fitting parameter values(iii)
The real parameter information of well V0 was obtained from the fitting parameters and well control parameters (Table3)Table 3
Interpretation parameters of well V0.
ParameterValueParameterValueOGIP (m3)2853.33Drainage radius (m)143.53Effective permeability (mD)0.24Fracture half-length (m)28.71Drainage area (km2)0.64Skin factor2.6According to the analysis of modern production data, the following understanding of well V0 was obtained. Compared with the diagnostic chart of tight sandstone, the production data qualities of gas wells were better. It can be seen from Figure7 that when the material balance time was longer than 0.7 days, the fitting rate between the actual data and the typical data was greater than 90%. When the material balance fitting time was shorter than 0.7 days, the fitting rate between the actual data and the typical curve was below 60%. The data noise was primarily caused by factors such as well pollution and inaccurate logging of the early production data. Although there was a certain gap between the actual well-controlled reserves and the previous gas reservoir description, the quality of the data was relatively high. Hence, the interpretation results can be judged accurately. The effective permeability was consistent with the overall permeability of the Shan 1 reservoir, and the fitting effect was good. The fracturing effect was far from the design goal, and the near wellbore skin pollution was serious.It is recommended to modify the geological modeling carrier and numerical simulation data based on the parameters with the good fitting effects (fracture half-length, formation permeability, and well-controlled reserves). The core of the target layer wellbore was taken, and a conventional gas reservoir test and analysis were conducted to study the seepage law of fluid in the porous media. Combined with the microseismic monitoring of well v0, secondary or multiple repeated fracturing was conducted. It is suggested that the length of the artificial fracture be extended to the outer boundary of the reservoir so as to form an effective area for the expanding volume fracturing and further increase the dominant channel for fluid seepage. In view of the serious skin pollution, it is recommended to implement an acid plugging removal operation, which plays an auxiliary role in improving the fluid seepage capacity.


### 4.2. Block Applications

Clarifying the scale of the effective sand body and reasonably designing the well spacing are important means to improve the production degree and recovery efficiency of tight gas reservoirs. This production data analysis method can determine the scale of effective sand bodies from a dynamic perspective and use the dynamic data of gas wells to analyze the controlled reserves and production range of gas wells, thus providing a basis for optimizing well spacing.Gas wells in tight gas reservoirs are typically put into production after fracturing. The typical dimensionless production analysis chart deduced in this study was compared with the actual production data of gas wells to fit the dimensionless production, the dimensionless production integral, and the dimensionless production derivative under dimensionless time under different ratios of the gas discharge radius to the fracture half-length. Combined with the production data diagnosis chart, the gas wells with better production data quality were selected to predict the fracture half-length and discharge radius. The Sulige gas field X block gas well evaluation was selected in this study. After the plate fitting analysis, the discharge radius was 130–390 m, with an average of 242 m, and the sand body width in this block was preliminarily estimated to be approximately 484 (Table4). Combined with the well logging data of the block, it was concluded that the size of the core beach and the distribution characteristics of the channel in the block were in good agreement with the conclusion obtained by the production data analysis method. Therefore, this conclusion provides strong evidence for the demonstration of the sand body size and guides the improvement and optimization of the well pattern in the block.Table 4
Evaluation of the dynamic drainage radius for the Sulige X block wells.
WellV1V2V3V4V5V6V7V8V9V10V11V12V13AverageDrainage radius/m340384362318256183161145134163279207210242Based on the dynamic understanding of the above sand body scale, the reasonable well pattern density suitable for this area is determined, which effectively guides the well pattern adjustment in this area. The production performance shows that after optimizing the well pattern, the interference between wells is low and the gas well production effect is better.


## 4.1. Single-Well Application

Well V0 in the Sulige gas field is a tight sandstone gas well that was produced after fracturing (artificial fracture 92 m at depth 637 m) in November 2008. The target layer is the reservoir of the Shan 1 member that is primarily composed of grey white-grey medium and coarse sandstone and grey-black fine sandstone. The debris content in the debris particles is high, and the rock types are primarily lithic quartz sandstone and quartz sandstone. The porosity of the reservoir is primarily 4.0–10.0%, with an average of 6.1%, and the permeability is primarily 0.100–1.000 mD, with an average of 0.265 mD. This is typical for a tight sandstone reservoir. In order to avoid a long time shut-in pressure recovery test, the modern production data analysis of well V0 was studied. The gas reservoir description data, production dynamic curve, and fitting curve are shown in Table2 and Figures 6 and 7, respectively.Table 2
Reservoir description parameters of well V0.
ParameterValueParameterValueNet pay (m)15Geological reserve (104 m3)3727.73Shaft radius (m)0.067Total porosity (%)9.5Reservoir temperature (K)383.15Initial viscosity of gas (mPa·s)0.0036Initial water saturation (%)68Initial volume factor0.0021Figure 6
Casing pressure and daily production curve of gas well V0.Figure 7
Typical plot matching result of well V0.The biggest difference between the production data analysis method of the vertical fractured wells in tight gas reservoirs and ordinary vertical wells is that the fracture half-length analysis is increased, and the interpretation results of this study also found this increased skin factor.The specific fitting and research steps were as follows:(i)
The material balance (pseudo) times at each production time point were obtained(ii)
Moving the processed production data points on the Blasingame-type dimensionless yields integral and integral derivative chart. The yield integral curve was first fitted because the curve was smoother than the yield decline curve. Then, the production decline curve and the integral derivative curve were used for the fitting. Finally, the time step and dimensionless discharge radius were adjusted to select the best fitting effect to obtain the fitting parameter values(iii)
The real parameter information of well V0 was obtained from the fitting parameters and well control parameters (Table3)Table 3
Interpretation parameters of well V0.
ParameterValueParameterValueOGIP (m3)2853.33Drainage radius (m)143.53Effective permeability (mD)0.24Fracture half-length (m)28.71Drainage area (km2)0.64Skin factor2.6According to the analysis of modern production data, the following understanding of well V0 was obtained. Compared with the diagnostic chart of tight sandstone, the production data qualities of gas wells were better. It can be seen from Figure7 that when the material balance time was longer than 0.7 days, the fitting rate between the actual data and the typical data was greater than 90%. When the material balance fitting time was shorter than 0.7 days, the fitting rate between the actual data and the typical curve was below 60%. The data noise was primarily caused by factors such as well pollution and inaccurate logging of the early production data. Although there was a certain gap between the actual well-controlled reserves and the previous gas reservoir description, the quality of the data was relatively high. Hence, the interpretation results can be judged accurately. The effective permeability was consistent with the overall permeability of the Shan 1 reservoir, and the fitting effect was good. The fracturing effect was far from the design goal, and the near wellbore skin pollution was serious.It is recommended to modify the geological modeling carrier and numerical simulation data based on the parameters with the good fitting effects (fracture half-length, formation permeability, and well-controlled reserves). The core of the target layer wellbore was taken, and a conventional gas reservoir test and analysis were conducted to study the seepage law of fluid in the porous media. Combined with the microseismic monitoring of well v0, secondary or multiple repeated fracturing was conducted. It is suggested that the length of the artificial fracture be extended to the outer boundary of the reservoir so as to form an effective area for the expanding volume fracturing and further increase the dominant channel for fluid seepage. In view of the serious skin pollution, it is recommended to implement an acid plugging removal operation, which plays an auxiliary role in improving the fluid seepage capacity.


## 4.2. Block Applications

Clarifying the scale of the effective sand body and reasonably designing the well spacing are important means to improve the production degree and recovery efficiency of tight gas reservoirs. This production data analysis method can determine the scale of effective sand bodies from a dynamic perspective and use the dynamic data of gas wells to analyze the controlled reserves and production range of gas wells, thus providing a basis for optimizing well spacing.Gas wells in tight gas reservoirs are typically put into production after fracturing. The typical dimensionless production analysis chart deduced in this study was compared with the actual production data of gas wells to fit the dimensionless production, the dimensionless production integral, and the dimensionless production derivative under dimensionless time under different ratios of the gas discharge radius to the fracture half-length. Combined with the production data diagnosis chart, the gas wells with better production data quality were selected to predict the fracture half-length and discharge radius. The Sulige gas field X block gas well evaluation was selected in this study. After the plate fitting analysis, the discharge radius was 130–390 m, with an average of 242 m, and the sand body width in this block was preliminarily estimated to be approximately 484 (Table4). Combined with the well logging data of the block, it was concluded that the size of the core beach and the distribution characteristics of the channel in the block were in good agreement with the conclusion obtained by the production data analysis method. Therefore, this conclusion provides strong evidence for the demonstration of the sand body size and guides the improvement and optimization of the well pattern in the block.Table 4
Evaluation of the dynamic drainage radius for the Sulige X block wells.
WellV1V2V3V4V5V6V7V8V9V10V11V12V13AverageDrainage radius/m340384362318256183161145134163279207210242Based on the dynamic understanding of the above sand body scale, the reasonable well pattern density suitable for this area is determined, which effectively guides the well pattern adjustment in this area. The production performance shows that after optimizing the well pattern, the interference between wells is low and the gas well production effect is better.


## 5. Conclusion

The difficulty of analyzing production data using typical curves is the establishment of typical curves and the diagnosis of production data. In this study, combined with the seepage characteristics of gas wells after fracturing in the Sulige gas field and based on the unstable pressure solution and the introduction of the material balance (pseudo) time, the typical Blasingame-type production data chart of vertical fractured wells in tight gas reservoirs was established. This was combined with field experience, and a set of diagnosis charts of the formation data of the tight sandstone gas reservoir was combined with the field experience.(1)
It was seen from the typical curves that the production decline curve could be divided into three flow stages: the first stage was an unstable decline section, and the curve did not coincide. Stage 2 was the transition section, where the curve gradually overlapped. Stage 3 was a pseudo-steady decreasing section with a slope of -1, which was close to the Arps harmonic decreasing trend(2)
Based on a large amount of field research work, a set of diagnostic curves for the production data of the tight sandstone gas reservoirs was established. The diagnostic curves can help engineers quickly discover the causes of the poor quality production data and select the corresponding solutions to improve the data quality and improve the accuracy of the gas well analysis(3)
The specific analysis steps of production data fitting were provided in detail. Using an analysis of the production data, such as unstable pressure and production, the parameters, such as the well control reserves, the effective permeability, and the skin factor were obtained. The calculation results can be used as a basis for the quantitative evaluation of gas well development effects(4)
This study expounded on the use of gas well dynamic data to analyze the controlled reserves and the production range of gas wells in the block. In addition, gas wells with good data quality or gas wells that can improve data quality for analysis were selected so as to obtain an accurate and effective sand body size of the tight sandstone. This will not only provide technical support for the reasonable optimization of the well pattern but also provide suggestions for adjustments of the gas reservoir development plan



---

*Source: 5549725-2021-07-01.xml*