Production Analysis of Tight Sandstone Reservoir in Consideration of Stress-Sensitive Permeability

Yinghao Shen1, *, Xinyu Yang2, Yuelei Zhang3
1 Unconventional Natural Gas Research Institute, China University of Petroleum, Beijing, P.R. China
2 CNPC Bohai Drilling Engineering Company Limited, Tianjin, P.R. China
3 CNODC Brasil Peteoleo e Gas Ltda, Rio de Janeiro, Brazil

Article Metrics

CrossRef Citations:
Total Statistics:

Full-Text HTML Views: 2526
Abstract HTML Views: 1558
PDF Downloads: 732
ePub Downloads: 302
Total Views/Downloads: 5118
Unique Statistics:

Full-Text HTML Views: 1226
Abstract HTML Views: 799
PDF Downloads: 490
ePub Downloads: 199
Total Views/Downloads: 2714

© 2017 Shen et al.

open-access license: This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 International Public License (CC-BY 4.0), a copy of which is available at: This license permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

* Address correspondence to this author the Unconventional Natural Gas Research Institite, China University of Petroleum, 18th Fuxue Road, Beijing, P.R. China; Tel: +86 18210117563; E-mail:



Tight sandstone reservoirs play an important role in the oil industry. The permeability of tight sandstone reservoir generally has stronger stress sensitivity than that of conventional reservoir because of the latter’s poor physical properties. However, the production analysis of tight sandstone reservoir did not fully considered the stress-sensitive permeability yet.


This paper proposed a production analysis method considering the stress- sensitive permeability.


This paper firtstly investigated the stress sensitivity characteristics and the effect of stress-sensitive permeability on a tight reservoir. Decline-type curves that consider stress-sensitive permeability are then established, and a systematic analysis method was built for the production analysis to obtain the single-well controlled dynamic reserves and reservoir physical properties.


A field analysis was performed in combination with Block Yuan-284 of Changqing Oilfield. Results show that with the reduction of reservoir pressure, stress sensitivity leads to the decline in reservoir permeability and the increase in seepage resistance, thus reducing the actual single-well controlled reserve and radius.


By utilizing the analysis method based on the decline curves, we can effectively predict the single-well controlled dynamic reserves of such reservoirs and evaluate the characteristic parameters of reservoirs.

Keywords: Unconventional reservoir, Stress sensitive, Tight sandstone, Production analysis, Decline-type curves.


Unconventional oil and gas resources have become an important aspect of current oil and gas production [1, 2], and stress sensitivity is an important feature of unconventional oil and gas reservoirs, such as tight gas reservoir, tight oil reservoir, and shale gas [3-8]. These reservoirs have poor physical properties, and the effect of stress sensitivity is also relatively more significant. Production analysis using decline-type curves is efficient in obtaining single-well controlled reserve and reservoir physical properties. This method is more advanced than the volumetric method because of the strong heterogeneity, very complex pore structure, and certain particularity in the seepage mechanism of unconventional reservoir [9-13]. The accurate evaluation and prediction of dynamic reserves is the basis of making a reasonable development program, and it plays a vital role in the reasonable and efficient exploitation of unconventional oil and gas resources.

The current production analysis methods do not consider the effect of stress sensitivity, and thus their calculation results produce a large deviation. This article brings stress-sensitive permeability into decline-type curves to obtain a more accurate production analysis, including single-well controlled dynamic reserves, to make the calculation result correspond with the actual situation of reservoir, and to effectively guide the development of stress-sensitive oil–gas reservoirs.


Pore pressure reduces constantly in the production of oil and gas, and thus effective stress increases accordingly. The corresponding relation is expressed in the following equation [14]:


where S is the total stress, σ is the effective stress (matrix stress and pressure between particles), a is the stress sensitivity correction coefficient, and is the fluid pressure or pore pressure.

Formula 1 shows that under the condition of unchanged total stress, the change in fluid pressure leads to the change in effective stress, which will correspondingly result in the change in rock physical property.

A permeability stress sensitivity experiment was conducted for the ultra-low permeability sandstone from Block Yuan-284 of Changqing Oilfield. The core belongs to the Triassic Chang-6 Reservoir Group of Ordos Basin. The main fragment component of the sandstone is feldspar lithic sandstone, which is fine-medium grain and very fine-fine grain. The stress sensitivity characteristics of sandstone were analyzed in the research area. The confining pressure was maintained, and the displacement pressure was changed to simulate the process in which the formation rock was compressed during exploitation. The stress sensitivity experiment was conducted using eight cores. Fig. (1) presents the curves of k/ki changing with pore pressure, where k/ki is the ratio of permeability and initial permeability under different pore pressure states. Fig. (1) also illustrates the experiment results. We can fit the results with the exponential law curves in which the fitting compliance degrees are all above 0.9.

Fig. (1). Permeability stress-sensitive experiment results for chang-6 tight sandstones.


Blasingame’s decline-type curve [15, 16] is an important method developed on the basis of Fetkovich’s decline-type curve [17-19] for production data analysis to evaluate single-well controlled dynamic reserves and reservoir parameters. However, to the best of our knowledge, this method currently does not consider the effect of stress sensitivity. In this study, we established a physical and mathematical model that considers stress-sensitive permeability to obtain a new series of decline-type curves that allow fr the stress-sensitivity of permeability.

Assuming the presence of a vertical well with a circular reservoir (Fig. 2), the formation and the fluid satisfy the following conditions:

Fig. (2). Schematic Diagram of a Vertical Well in a Circular Closed Formation.

  1. The reservoir is homogeneous and has a uniform thickness of h.
  2. The original formation pressure is Pi.
  3. The oil well runs throughout the whole reservoir with a wellbore radius of rw and constant production rate q.
  4. The outer boundary of the reservoir is closed, and the radius of the outer boundary is re.
  5. Permeability presents an exponential change law with pore pressure.
  6. The formation fluid is Newtonian fluid.
  7. The effect of capillary force, temperature, and gravity is ignored.

The basic flow equation of a homogeneous reservoir is


When the permeability’s change in pressure follows an exponential law, the permeability modulus is effective in characterizing the permeability variation law. The relation of permeability change with effective stress can be expressed as follows [20]:


where γ is the permeability modulus of 1 MPa.

The integral for Equation (3) is


where ki is the reservoir permeability under the original formation pressure, mD.

Equation (4) is the expression of the exponential change law of permeability with pressure.

By bringing (4) into (2), we obtain the following eq:


The above formula is the basic percolation differential equation of a tight sandstone reservoir that considers the stress-sensitive effect of permeability. The following dimensionless variables are introduced:




The following dimensionless permeability modulus is defined as:


According to flow mechanics theory, the mathematical model of a homogeneous closed reservoir that considers stress sensitivity can be established.

The flow control equation is:


The initial conditions are expressed as:


The internal boundary are expressed as:



The external boundary conditions are expressed as:


To calculate the above model, we used the perturbation technique for calculation and introduced the transformation relation as follows:


The following perturbation techniques are applied as follows:





Given the small dimensionless permeability modulus, the zero-order perturbation solution is used. Thus, Equations (10–14) are expressed as follows:


ηOD(rD,0) = 0 (21)




By conducting the Laplace transform for Equation (20), its Laplace spatial solution can be obtainedas:




In the above functions, KO and K1 are a zero-order and a first-order modified Bessel function, respectively, g is the Laplace Variable, L-1 is the Laplace inverse transform, and O(γD), ηwD is the allowance above zero-order solution.

By transforming the above analytical solution into a real spatial solution through the Stehfest inversion algorithm, we obtained PwD~tD. A typical decline-type curve mainly includes the decline-type curve dimensionless production proposed by Fetkovich, qDd, the dimensionless material balance time tDd proposed by Blasingame, and the dimensionless production integral and its derivative proposed by McCray, qDdi and qDid.

The decline-type curve dimensionless production and time are defined as follows:



Here, rwe=rw·e-s.

The auxiliary drawing function generally includes the dimensionless production integral and its derivative proposed by McCray. It is introduced for practical application and is conducive to improving the efficiency and accuracy of field engineers to explain and analyze production data. In actual use, the method of fitting three groups of data simultaneously is generally adopted to obtain the best fitting effect.

The following is the dimensionless production integral:


The following is the dimensionless production integral derivative:


By drawing qDd tDd, qDdi tDd and qDdid tDd on the log-log plot, decline-type curves for vertical wells in homogeneous reservoirs of stress-sensitive reservoir are formed.

Fig. (3) presents the decline-type curve under different dimensionless permeability moduli γD when the dimensionless reservoir radius reD is 5000, and the values of γD are 0, 0.0015, 0.003, and 0.0045 respectively. The difference between curves in the unsteady stage in the figure is small, and it is mainly in the pseudo steady stage. The curves are shown when γD is 0, 0.0015, 0.003, and 0.0045 from right to left.

Fig. (3). Decline-type Curves for Vertical Wells in Stress-sensitive Reservoirs (Different Permeability Moduli).

Fig. (4) illustrates the decline-type curves under different dimensionless reservoir radii reD when the dimensionless permeability modulus γD is 0.003 and the values of reD are 200, 500, 2000, and 5000, respectively. The difference exists both in the unsteady stage and in the quasi-steady stage. The curves are shown when reD is 200, 500, 2000, and 5000 from top to bottom.

Figs. (3) and (4) summarize the following basic laws of the decline-type curves for vertical wells in stress-sensitive reservoirs:

Fig. (4). Decline-type Curves for Vertical Wells in Stress-sensitive Reservoirs (Different Reservoir Radii).

  1. Stress-sensitive permeability affects the decline-type curves, mainly those existing in the pseudo steady stage.
  2. Stress-sensitive permeability not only affects the production curve patterns but also the production integral and production integral derivative curve to a certain extent.
  3. After the curve reaches the pseudo steady state when stress sensitivity occurs, the production curve slope of the pseudo steady state for the conventional reservoir does not follow -1 any longer. The permeability is reduced because of the effect of stress sensitivity, and thus production decreases within the scope of discussion time.
  4. Under the same reservoir radius, the lowering degree of the curve increases with the increase in permeability modulus. Therefore, when other conditions are the same, the effect on production and the declining range of production becomes greater with the increase in stress-sensitive degree.
  5. Under the same dimensionless permeability modulus, the larger the dimensionless radius of the reservoir is, the greater is the lowering degree. That is, the effective transformation of the near-well-bore area can reduce the effect of stress sensitivity on the production.


An instance analysis was conducted by selecting a well in Block Yuan 284 that belongs to Ordos Basin. Block Yuan 284’s purpose of horizon is the Chang-6 formation. The effective thickness of this well was 10.94 m, with porosity being 0.13, irreducible water saturation being 43%, and well radius being 0.1 m. Its production history is illustrated in Fig. (5). The history of this well was studied first. No other operation measures wereadopted after the fracturing production of this well. Its production and pressure data were monitored daily that are representative of relatively high accuracy and good correlation (Fig. 5). The water injection enhancement productivity effect was poor, pressure and production declined with time, and the oil well remained at the depletion development state.

This well does not undergo unstable pressure testing after being put into production. Therefore, we mainly used its production data to judge the model as shown in Fig. (6).

The basic model characteristics are as follows:

  • Pseudo-steady characteristicsare observed . The tightness of reservoir leads to poor flowability and a very slow effect of the water injection well. The dynamic response characteristics indicate that this well does not take effect, and thus the oil well gradually achieves the quasi-steady state with no peripheral pressure response characteristics.
  • Hydraulic fracture characteristicshave been observed . The flow characteristics of the non-steady stage do not conform to the characteristics of the hydraulic fractured well. Although fracturing production is made for this well, the role of fractures in this well is to improve the near-well-bore permeability because of poor reservoir permeability. The physical property of the reservoir periphery is relatively poor, and it is a common feature of this well area that can be observed from the unstable pressure testing curve of similar wells in this well area.
Fig. (5). Production History of the Target Well.

Fig. (6). Curves for the Normalized Productivity Index Function vs. Material Balance Time.

To sum up, the basic model of this well is a vertical well within a stress-sensitive closed reservoir. This well can be the model of a vertical well within a stress-sensitive closed reservoir. Through the exponential curve fitting, we obtained its permeability modulus of 0.0097 (1/MPa). The dimensionless permeability modulus of 0.0031 was calculated by Equation (7), to select a reasonable decline chart series from it.

To calculate the material balance time, the material balance time of the oil well is defined as follows:


The calculated normalized production, production integral, and its derivative are defined as follows:

Normalized production


Production integral (defined as the average production of well at a certain time)


Production integral derivative


A log-log plot was drawn for the curve of three production functions versus material balance time, with its coordinate scale, the same as that of the typical decline chart, and a typical decline chart with γD equal to 0.003 was used. Through fitting, the production curve can be better fitted with a typical decline-type curve with a dimensionless reservoir radius reD equal to 200 (Fig. 7).

Fig. (7). Fitting Curve of the Typical Decline-type Chart.

The following fitting parameters can be obtained through fitting:


(tDd)MP=1.0, (tc)MP=6000(h)

(qDd)MP=1.0, (qO/ p)MP=0.38(m3/ (dMPa))

The typical curve analysis and calculation result are as follows:

Single-well controlled dynamic reserves


Single-well controlled dynamic reserves expressed by quality


Reservoir discharge area


Equivalent control radius


Original effective permeability


Skin coefficient

s = -ln[(re /reD)/ rw] = -2.01. (41)

After the model and the parameters were determined, the production history was required to be fit to verify the correctness of the model and the parameters. The well is a vertical well model within a stress-sensitive closed reservoir, and the basic parameters for fitting have been determined through the chart. The history fitting curve is presented in Fig. (8).

Fig. (8). Production History Fitting Curve.

Fig. (8) indicates that the fitting degree is high, and thus the adopted model and parameters have high credibility and can correctly reflect the actual status of the oil well.


  1. The permeability of Chang-6 tight sandstone in Block Yuan-284 of Changqing Oilfield had significant stress sensitivity characteristics, and thus stress-sensitive permeability needed to be considered in the analysis of reservoir engineering to improve the accuracy and reliability of analysis.
  2. This article introduced the stress sensitivity of permeability, established the seepage control equation in consideration of stress-sensitive permeability, and calculated it. Thus, the typical decline chart under stress sensitivity was obtained.
  3. The typical decline chart indicated that stress sensitivity decreased the actual single-well controlled dynamic reserves and the single-well controlled radius. The reason for this phenomenon is that, with the reduction of reservoir pressure, stress sensitivity leads to the lowering of reservoir permeability and the increase in seepage resistance.
  4. An sample application was performed according to the established typical decline chart in consideration of the stress sensitivity of permeability and in combination with a certain well in Block Yuan 284 of Changqing Oilfield. The analysis method based on this chart can effectively predict the dynamic reserves and evaluate the characteristic parameters of the reservoirs.


The authors confirm that this article content has no conflict of interest.


This work is supported by the National Natural Science Foundation of China (Grant No. 51604287, 51490652).


[1] Bentley RW. Global oil & gas depletion: an overview. Energy Policy 2002; 30(3): 189-205.
[2] Sorrell S, Speirs J, Bentley JR, Brandt A, Miller R. Global oil depletion: A review of the evidence. Energy Policy 2010; 38(9): 5290-5.
[3] Holditch SA. Tight gas sands. SPE Journal of Petroleum Technology 2006; 58(6): 86-93.
[4] Jones FO, Owens WW. A laboratory study of low-permeability gas sands. SPE Journal of Petroleum Technology 1980; 32(9): 1631-40.
[5] Vairogs J, Hearn CL, Dareing DW, Rhoades VW. Effect of rock stress on gas production from low-permeability reservoirs. SPE Journal of Petroleum Technology 1971; 23(9): 1161-7.
[6] Wang L, Liu H, Zhen S, Liu X, Li J, Song Y. Research on quantitative interpretation for stress sensitivity of low-permeability reservoir. Acta Petrolei Sinica 2009; 30(01): 96-9.
[7] Jiao C, He S, Xie Q, et al. Experiment for stress sensitivity of ultra-low permeability sandstone reservoir. Acta Petrolei Sinica 2011; 32(3): 489-94.
[8] Zhu L, Ji Y, Yang T, Li X. Analysis on production of coal bed methane considering the change in permeability of coal rock. Open Petroleum Engineering Journal 2016; 9: 289-98.
[9] Wang S. Research on Predication Method for Dynamic Reserves of Low-permeability Gas Reservoir. Beijing: China University of Petroleum (Beijing) 2009.
[10] Wang W, Shen P, Ma X, Fan L, Tang J. Determination of dynamic reserves for heterogeneous complicated low-permeability gas reservoir. Natural Gas Industry 2004; 24(07): 80-2.
[11] Hao Y, Xu X. Correct calculation of dynamic reserves for low-permeability gas reservoir. Acta Petrolei Sinica 2002; 29(5): 66-8.
[12] Ilk D, Rushing JA, Blasingame TA. Integration of production analysis and rate-time analysis via parametric correlations -- theoretical considerations and practical applications Proceeding of SPE Hydraulic Fracturing Technology Conference.
[13] Cipolla CL, Lolon E, Erdle J, Tathed V. Modeling well performance in Shale-Gas Reservoirs Proceeding of SPE/EAGE Reservoir Characterization and Simulation Conference, Society of Petroleum, Engineers, Abu Dhabi, UAE, 2009.
[14] Dobynin VM. Effect of overburden pressure on some properties of sandstones. SPE Journal 1962; 2(4): 360-6.
[15] Blasingame TA, McCray TL, Lee WJ. Decline curve analysis for variable pressure drop/variable flowrate systems In: Proceeding of SPE Gas Technology Symposium. Houston, Texas: Society of Petroleum Engineers 1991.
[16] Palacio J, Blasingame T. Decline curve analysis using type curves: analysis of gas well production data In: Proceeding of SPE Joint Rocky Mountain Regional and Low Permeability Reservoirs Symposium. Denver, USA: Society of Petroleum Engineers 1993.
[17] Fetkovich MJ. Decline curve analysis using type curves. SPE Journal of Petroleum Technology 1980; 32(06): 1065-77.
[18] Fetkovich M, Vienot M, Bradley M, Kiesow UG. Decline curve analysis using type curves: case histories. SPE Formation Evaluation 1987; 2(4): 637-56.
[19] Yin D, Wang D, Zhang C, Duan Y. Shale gas productivity predicting model and analysis of influence factor. Open Petroleum Engineering Journal 2015; 8: 203-7.
[20] Pedrosa OA. Pressure transient response in stress-sensitive formations Proceeding of SPE California Regional Meeting.