Open access peer-reviewed chapter

Gas Well Testing

Written By

Freddy Humberto Escobar

Submitted: 22 September 2015 Reviewed: 26 January 2017 Published: 02 August 2017

DOI: 10.5772/67620

From the Edited Volume

Advances in Natural Gas Emerging Technologies

Edited by Hamid A. Al-Megren and Rashid H. Altamimi

Chapter metrics overview

1,965 Chapter Downloads

View Full Metrics

Abstract

Modeling liquid flow for well test interpretation considers constant values of both density and compressibility within the range of dealt pressures. This assumption does not apply for gas flow case in which the gas compressibility factor is also included for a better mathematical representation. The gas flow equation is normally linearized to allow the liquid diffusivity solution to satisfy gas flow behavior. Depending upon the viscosity-compressibility product, three treatments are considered for the linearization: square of pressure squared, pseudopressure, or linear pressure. When wellbore storage conditions are insignificant, drawdown tests are best analyzed using the pseudopressure function. Besides, since the viscosity-compressibility product is highly sensitive in gas flow; then, pseudotime best captures the gas thermodynamics. Buildup pressure tests, for example, require linearization of both pseudotime and pseudopressure. The conventional straight-line method has been customarily used for well test interpretation. Its disadvantages are the accuracy in determining of the starting and ending of a given flow regime and the lack of verification. This is not the case of the Tiab’s Direct Synthesis technique (TDS) which is indifferently applied to either drawdown or buildup tests and is based on features and intersection points found of the pressure and pressure derivative log-log plot.

Keywords

  • TDS technique
  • pseudotime
  • pseudopressure
  • rapid flow
  • viscosity
  • rate transient analysis
  • pressure transient analysis

1. Introduction

Contrary to liquids, a gas is highly compressible and much less viscous. In general, gas viscosity is about a 100 times lower than the least viscous crude oil. It is important, however, to try to provide the same mathematical treatment to oil and gas hydrocarbons, so interpretation methodologies can easily be applied in a more practical way. Then, the gas flow equation is normally linearized to allow the liquid diffusivity solution to satisfy the gas behavior when analyzing transient test data of gas reservoirs. Depending on the values of reservoir pressure, viscosity, and gas compressibility factor, the gas flow behavior can be treated as a function of either pressure to the second power or linear pressure with a region which does not correspond to any of these and it is better represented by a synthetic function call pseudopressure. Pseudopressure is a function that integrates pressure, density, and compressibility factor. The gas system’s total compressibility highly depends on gas compressibility which for ideal gases changes inversely with the pressure. Then, another artificial function referred as pseudotime is included to further understand the transient behavior of gas flow in porous media. For instance, when wellbore storage conditions are insignificant, drawdown tests are best analyzed using the pseudopressure function. On the other hand, buildup pressure tests require linearization of both pseudotime and pseudopressure.

This chapter will be devoted to provide both fundamental of gas flow in porous media as well as interpretation of pressure and rate data in gas reservoirs. The use of the oil flow equations and interpretation techniques is carefully extended for gas flow so that reservoir permeability, skin factor, and reservoir area can be easily estimated from a gas pressure or gas rate test by using conventional analysis and characteristic points found on the pressure derivative plot (TDS technique). Conventional analysis—the oldest pressure transient test interpretation technique—is based upon understanding the flow behavior in a given reservoir geometry, so the pressure versus time function is plotted in such way that a linear trend can be obtained. Both slope and intercept of such linear tendency are used to characterize the reservoir. Conventional analysis has two main drawbacks: (1) difficulty of finding a given flow regime and (2) absence of parameter verification. On the other hand, TDS technique—is strongly based on the log-log plot of pressure and pressure derivative versus time curves which provide the best way for flow regime identification; then, it uses the “fingerprints” or characteristic points found in such plot which are entered in practical and direct analytical equations to easily find reservoir parameters. Moreover, the same parameters can be obtained from different sources for verification purposes. Such is the case, for instance, of the reservoir area in elongated systems which can be estimated five times.

The chapter will include both interpretation techniques TDS and conventional in two cases: (1) infinite and (2) finite reservoirs. Channels or elongated systems in which reservoir hemilinear, parabolic or linear flow regimes developed once radial flow regime vanishes are reported in Refs [8, 13, 14]. This formation of linear flow regime normally occurs in fluvial deposits (channels), sand lens, parallel faulting, terrace faulting, and carbonate reefs. Then, such systems are worth of transient pressure analysis characterization. Latest researches on the determination of drainage area in constant-pressure-bounded systems using either conventional analysis or TDS technique are also reported by Escobar et al. [10].

It is convenient to mention some other important aspects concerning gas well testing which have appeared recently. The first case is the transient rate analysis in hydraulically fractured wells which was presented by [19] for both oil and gas wells. The traditional model for elliptical flow included the reservoir area as a variable. Handling the interpretation using TDS Technique may be little difficult for unexperienced interpreters. Therefore, [20] introduced a model excluding the reservoir drainage area and avoiding the necessity of developing pseudosteady-state regime. When a naturally fractured reservoir is subjected to hydraulic fracturing, the interpretation should be performed according to the presented by [21]. [35] presented the pressure behavior of finite-conductivity fractured wells in gas composite systems. As far a horizontal wells, the recent works by [23] and [24] included off-centered wells for transient-rate or transient-pressure cases, respectively. [29] presented a study of production performance of horizontal wells when rapid flow conditions are given.

Practical exercises will provide in the chapter provide a better understanding and applicability of the interpretation techniques.

The purpose of this chapter is two folded: (1) to present the governing equation for gas flow used in well test interpretation and (2) to use both conventional and TDS Techniques as valuable tools for well test interpretation in both transient rate and transient pressure analysis. Some detailed examples will be given for demonstration purposes.

Advertisement

2. Transient pressure analysis

Transient pressure analysis is performed measuring the bottom-hole pressure while the flow rate is kept constant.

2.1. Fluid flow equations

The gas diffusivity equation in oil-field units is given by:

1rr(Pμ(P)Z(P)rPr)=φ0.0002637t(PZ(P))E1

Which can be modified to respond for three-phase flow (oil, water, and gas):

1rr(rPr)=φct0.0002637λtPtE2

where, the total compressibility, ct, and total mobility, λt, are given by:

ctcgsg+coSo+cwSw+cfE3
λt=kgμg+koμo+kwμwE4

As can be inferred from Eq. (3), the total compressibility varies significantly when dealing with monophasic gas flow since the gas compressibility varies along with the pressure. Agarwal [1] introduced the pseudotime function to alleviate such problem. This function accounts for the time dependence of gas viscosity and total system compressibility:

ta=treftdtμ(t)ct(t)E5

Pseudotime is better defined as a function of pressure as a new function given in hr psi/cp:

ta(P)=PrefP(dt/dP)μ(P)ct(P)dPE6

Notice that μ and ct are now pressure-dependent properties.

As expressed by Eq. (1), viscosity and gas compressibility factor are strong functions of pressure; then, to account for gas flow behavior, Al-Hussainy et al. [2] introduced the pseudopressure function which basically includes the variation of gas viscosity and compressibility into a single function which is given by:

m(P)=2PrefPPμ(P)Z(P)dPE7

After replacing Eqs. (6) and (7) into Eq. (1), it yields:

1rr(rm(P)r)=φ0.0002637km(P)ta(P)E8

Contrary to liquid well testing, rapid gas flow has a strong influence on well testing, [32]. As the flow rate increases, so does the skin factor, then:

sa=s+DqE9

Eq. (9) shows that the apparent skin factor is a function of the mechanical skin factor—which is assumed to be constant during the test— and the product of the flow rate with the turbulence factor or non-Darcy term. This implies that two flow test ought to be run at different flow rates to find mechanical skin factor and the turbulence factor from:

(sa)1=s+Dq1E10
(sa)2=s+Dq2E11

Solving the simultaneous equations:

D=(sa)1(sa)2q1q2E12
s=(sa)1(sa)1(sa)2q1q2q1E13

where, the skin factors 1 and 2 are estimated from each pressure test. However, there is a need of estimating the turbulence factor by empirical correlations for buildup cases or when a single test exists. Then, the non-Darcy flow coefficient is defined by [26]:

D=2.222×1015γgkhβμgrwhp2E14

The above equation is also applied to partially completed or partially penetrated wells. hp is the length of the perforated interval. For fully perforated wells, hp = h.

Parameter β is called turbulence factor or inertial factor can be found by correlations. The correlation proposed by Geertsma [21] is given by:

β=4.851×104φ5.5kE15

The consideration on the skin factor effect on gas testing was recognized by Fligelman et al. [25] who provided correction charts to account for apparent skin factor values.

2.2. Conventional analysis

The solution to the transient diffusivity equation, Eq. (8), is given by:

m(P)D(1,tDa)=12Ei(14tDa)E16

The dimensionless parameters used in this chapter are given below. The rigorous dimensionless time is:

tD=0.0002637ktφ(μct)rw2E17

Including the pseudotime function, ta(P), the dimensionless pseudotime is:

tDa=(0.0002637kφrw2)ta(P)E18

Notice that the viscosity-compressibility product is not seen in Eq. (16) since they are included in the pseudotime function. However, if we multiply and, then, divide by (μct)i, a similar equation to the general dimensionless time expression will be obtained.

tDa=(0.0002637kφ(μct)irw2)[(μct)i×ta(P)]E19

The dimensionless pseudopressure and pseudopressure derivatives are:

t*Δm(P)D'=hk[t*Δm(P)']1422.52qscTE19a
m(P)D=hk[m(Pi)m(P)]1422.52qTE20
ta(P)*Δm(P)D'=hk[ta(P)*Δm(P)']1422.52qscTE21

And the dimensionless wellbore storage coefficient is given by:

CD=(0.8935φhctrw2)CE22

The dimensionless radii are given:

rD=rrwE23
rDe=rerwE24

For practical purposes, Eq. (16) will end up in a semilog behavior of pseudopressure drops against time. After replacing the respective dimensionless quantities into the mentioned straight-line semilog expression, it is obtained [4]:

m(Pi)m(Pwf)=1.422×106qTkh[1.1513log(kt1688φ(μgct)irw2)+s'+D]E25
m(Pi)m(Pwf)=1.422×106qTkh[1.1513log(kta(P)1688φrw2)+s'+D]E26

The above equations are applied during transient or radial flow regime. They are used to find reservoir transmissibility and apparent skin factor from the slope and intercept, respectively, of a semilog plot of well-flowing pressure versus time. After applying the superposition principle, the above equations for the buildup case are converted into:

m(Pi)m(Pwf)=1.422×106qTkhlog(tp+ΔtΔt)E27
m(Pi)m(Pwf)=1.422×106qTkhlog(ta(P)p+Δta(P)Δta(P))E28

From a semilog plot of pseudopressure versus time (or pseudotime), its slope allows calculating the reservoir permeability and the intercept is used to find the pseudoskin factor, respectively:

k=1637.74qTmhE29
s'=[m(Pi)m(P1hr)mlog(kφ(μct)irw2)3.227+0.8686]E30

Notice that for the pseudotime case, (µct)i product in the above equation will be set as the unity. The gas pseudoskin factor is estimated for the buildup case as:

s'=[m(P1hr)m(Pwf)mlog(kφ(μct)irw2)3.227+0.8686]E31

The governing dimensionless pressure equation during pseudosteady-state period is given by [28]:

m(P)D=2tDreD2+lnrD0.75+s'E32

By replacing the dimensionless quantities, changing the log base, the above equation leads to:

m(Pi)m(Pwf)=0.2395qTtAhφ+3263qTkh[log0.472rerw+s'2.303]E33

A Cartesian plot of m(Pwf) versus time or pseudotime during pseudosteady state will yield a straight line in which slope, m*, is useful to find the well drainage area:

A=0.23395(5.615)qTφhm*E34

Such deliverability tests as backpressure, isochronal, modified isochronal, and flow after flow are conducted for the purpose of determining the flow exponent n (n = 1 is considered turbulent flow and 0.5 < n < 1 is considered to be rapid flow) and the performance coefficients. They assumed that stabilization is reached during the testing which is not true in most of the cases. Then, they are not included in this chapter but can be found in Chapter 4 of Ref. [4].

2.3. TDS technique

Tiab [33] proposed a revolutionary technique which is very useful to interpret pressure tests using characteristics points found on the pressure and pressure derivative versus time log-log plot. He obtained practical analytical solutions for the determination of reservoir parameters.

m(P)Dr=(7.029×104khqT)(m(Pi)m(Pwf)(tn)qn)=12(lntD+0.80907+2s)E35

From a log-log plot of pseudopressure and pseudopressure derivative against pseudotime, Figure 1, several main characteristics are outlined:

Figure 1.

Log-log plot of pseudopressure and pseudopressure derivative versus pseudotime. After Ref. [7].

1. The early unit-slope line originated by wellbore storage is described by the following equation:

m(P)D=tDaCDE36

Replacing the dimensionless parameters in Eq. (36), a new equation to estimate the wellbore storage coefficient is obtained:

C=(0.419qTct)(ta(P)Δm(P))E37

2. The intersection of the early unit-slope line with the radial horizontal straight line gives:

(tDaCD)i=0.5E38

From this, an equation to estimate either permeability or wellbore storage is obtained once the dimensionless parameters are replaced.

As presented by Tiab [33], the governing equation for the well pressure behavior during radial flow reformulated by Escobar et al. [7] in terms of pseudofunctions is expressed by:

ta(P)i=1695ctCkhE39

3. According to Ref. [28], another form of Eq. (35) is obtained when wellbore storage and skin factor are included:

m(P)Dr=12{ln(tDaCD)r+0.80907+ln(CDe2s)}E40

From the above equation, the derivative of pseudopressure with respect to the natural log of tDa/CD is given by:

[tDaCDm(P)'D]r=0.5E41

From Eq. (21), the dimensionless pseudopressure derivative with respect to the natural log of log tDa/CD gives:

[tDaCDm(P)'D]r=[7.029×104khqT][ta(P)*m(P)']E42

Combination of Eqs. (41) and (42) will result into an equation to estimate permeability:

k=711.26qTh[ta(P)*Δm(P)']rE43

3. Dividing Eq. (40) by Eq. (41), replacing the dimensionless quantities and, then, solving for the pseudoskin factor will yield:

s'=0.5[[Δm(P)]r[ta(P)*Δm(P)']rln(k(ta(P))rφrw2)+7.4316]E44

Finally, the pressure derivative during the pseudosteady-state flow regime of closed systems is governed by:

tDa*m(P)'D=2πtDaE45

The intersection point of the above straight line and the radial flow regime straight line is:

taDARPi=14πE46

After substituting the dimensionless pseudotime function into Eq. (46), a new equation for the well drainage area is presented:

A=kta(P)rpi301.77φE47

Further applications of gas well test can be found in the literature. Escobar et al. [12] introduced the mathematical expressions for interpretation of pressure tests using the pseudopressure and pseudopressure derivative as a function of pseudotime for hydraulically fractured wells and naturally fractured (heterogeneous) formations. Fligelman [30] presented an interpretation methodology using TDS technique for finite-conductivity fractured wells. They used pseudopressure and rigorous time. In 2012, Escobar et al. [16] implemented the transient pressure analysis on gas fractured wells in bi-zonal reservoirs. Moncada et al. [31] extended the TDS for oil and gas flow for partially completed and partially penetrated wells. As far as horizontal wells, it is worth to mention the work performed in Refs. [11] and [15] on homogeneous and naturally fractured reservoirs.

2.4. Example 1

Chaudhry [4] presented a reservoir limit test for a gas reservoir (example 5-2 of Ref. [4]). However, once the pressure derivative was taken to the test data, no late pseudosteady state regime was observed. Then, the input data given below were used to simulate a pressure test given in Table 1.

Sg = 70%Sw = 30%q = 6184 MSCF/D
h = 41 ftk = 44 mdBg = 0.00102 ft3/STB
rw = 0.4271 ftφ = 10.04%ct = 0.0002561 psi−1
ωg = 0.0992 md/cpγg = 0.732Pcr = 380.16 psia
Tcr = 645.06 RT = 710 Rre = 2200 ft (349 Ac)
m(Pi) = 340920304.2 psi2/cp

t, hrP, psit, hrP, psi
039651.27133677.2527
0.0013960.6291.60053670.3779
0.0023956.43132.01483663.602
0.0033952.37742.53653656.8388
0.0043948.45163.19333650.1477
0.0053944.64314.02023643.5144
0.0063940.94385.01073637.2126
0.0073937.34696.01073632.0323
0.0083933.84667.01073627.6664
0.0093930.43828.01073623.8936
0.01133922.82779.01073620.5718
0.01433913.840210.01073617.6047
0.0183903.257312.01073612.4479
0.02263891.132521.01073596.5482
0.02853877.500630.01073586.4447
0.03583862.481639.01073579.0231
0.04513846.205448.01073573.1536
0.05683829.259257.01073568.2967
0.07153812.132966.01073564.1453
0.093795.233575.01073560.4805
0.11333779.268684.01073557.2312
0.12713771.759493.01073554.3136
0.163757.7512102.01073551.6672
0.20153745.1803179.51073535.5164
0.25373734.0438269.51073523.5847
0.31933724.1258359.51073514.1663
0.4023715.1658449.51073505.545
0.50613706.8112539.51073497.2804
0.63723698.9669629.51073489.1576
0.80213691.489719.51073481.0959
1.00983684.2741809.51073473.056
ta(P), psi-hr/cpΔm(P), psi2/cpta(P)* Δm(P)′, psi2/cpta(P), psi-hr/cpΔm(P), psi2/cpta(P)* Δm(P)′, psi2/cp
00045802.4141870739.38417533667.77
37.7745270218.7615265698.015957569.76481838960.00317954428.54
75.5167518758.6126520866.104172365.78511815682.59818369258.28
113.2278751325.0471767304.190690970.24671799003.72618783428.06
150.9089968476.59021005966.456114363.69231786548.91819193297.17
188.5611172014.9041237504.102143779.30511774358.75619599754.75
226.1851362245.2141462410.249178978.26671766967.69919986008.75
263.78171544967.1391681099.321214477.4741764620.95220303606.64
301.3521713052.7971893921.256249946.70711752899.59420571329.92
338.89661876882.1932101169.292285390.45871717299.58720802733.84
426.29232225474.8942563944.683320812.05231692510.18221006509.26
536.13682605002.3023110507.871356214.04731681851.99921188553.27
674.15633004035.2183754179.95426967.52861725356.9821505009.16
847.52553404366.4314491744.317744772.17061759676.34122480949.29
1065.2353781729.3485321159.0891061858.0721799660.97423101515.46
1338.55114101068.9156235191.6451378463.161769298.72823557564.54
1681.59344335396.3637226023.9631694705.5111759287.15423918367.9
2112.07684462372.6668257991.5142010656.3171770713.33524217012.8
2652.25584463511.7599301323.1232326363.1011736601.94424472345.81
3330.09954300906.00510331215.572641858.7651709368.01224697796.42
4180.80194056795.1211304325.652957168.151701504.5924897723.65
4684.38083917486.95111762181.873272311.761711704.51225077261.23
5880.95673569692.65812616571.143587306.1051744436.56525240114.57
7383.62183191507.54813383609.836295436.7941821464.18326234973.81
9271.2022873762.47714063387.959432920.2822036512.44126970809.17
11642.90042591188.80914668995.8112564795.272355209.64427552034.83
14623.5032366500.19915216145.54156919652740929.32328084374.29
18369.90032197999.34715726460.218814758.133212985.0328594968.02
23079.35182076483.3916205683.1821933329.443695740.20529097073.99
28999.92391985415.32516662774.4925047755.954195533.85529595671.37
36443.53741918507.60517104058.5928158083.334703995.7230093194.56

Table 1.

Pressure, pseudopressure, time, and pseudotime data for example 1.

Estimate permeability, skin factor, and drainage area by both conventional analysis and TDS technique.

2.4.1. Solution by conventional analysis

Figure 2 presents a semilog pressure of pseudopressure versus pseudotime. The slope and intercept of the radial flow regime straight line in such plot are given below:

m=3995147.42 (psi2/cp)/(log hrpsi/cp)E84
m(P)1hr=342125555.5 psi2/cpE85

Figure 2.

Semilog plot for example 1.

Use of Eqs. (27) and (28) allows finding reservoir permeability and pseudoskin factor, respectively:

k=1637.74qTmh=1637.74(6184)(710)(3995147.42)(41)=43.45mdE8000
s'=[340920304.25342125555.53995147.42log(43.45(0.1004)(0.42712))3.227+0.8686]=0.5172E8001

To find the well drainage area, the Cartesian plot given in Figure 3 was built. Its slope, m* = 0.0914 (psi2/cp)/(hr−psi/cp), is plugged into Eq. (34):

A=0.23395qTφhm*=0.23395(6184)(710)(0.1004)(41)(0.0914)(43560)=336.4AcE9003

Figure 3.

Cartesian plot for example 1.

2.4.2. Solution by TDS technique

Figure 4 presents the pseudopressure and pressure derivative versus pseudotime log-log plot in which wellbore storage, radial flow regime, and late pseudosteady-state regimes are clearly observed. The following characteristic points were read from Figure 4:

ta(P)r = 1694705.5 psi hr/cpm(P)r = 23918367.9 psi2/cp
ta(P)*∆m(P)’r = 1735066.96 psi2/cpta(P)rpi = 10113641.48 psi hr/cp

Figure 4.

Pseudopressure drop and pseudopressure derivative versus time log-log plot for example 1.

Permeability and pseudoskin factor are respectively estimated from Eqs. (42) and (44),

k=711.26qTh[ta(P)*Δm(P)']r=711.26(6184)(710)(41)(1735066.96)=43.9mdE9005
s'=0.5[23918367.91735066.96ln(43.9(1694705.5)(0.1004)(0.42712))+7.4316]=0.454E9008

and well drainage area is found with Eq. (47):

A=kta(P)rpi301.77φ=(43.9)(10113641.48)301.77(0.1004)=336.4AcE9009

Finally, the inertial factor and the non-Darcy flow coefficient are estimated with Eqs. (14) and (15):

β=4.851×104φ5.5k=4.851×104(0.1004)43.9=2265091235.63ft1E9100
D=2.222×1015(0.732)(41)(43.9)(2265091235.63)(0.0992)(0.4271)(412)=9×105D/MscfE9101

The true skin factor is found with Eq. (9):

sa=s+Dq=0.454+9×105(6184)=1.42E9011

It can be seen that the simulated parameters closely match the results obtained from the examples.

Advertisement

3. Transient rate analysis

Transient rate analysis is performed by recording the continuous changing flow rate under a constant bottom-hole pressure condition. This procedure is normally achieved in very low gas formations and shale gas systems.

3.1. Basic flow and dimensional equations

The Laplace domain, the rate of solution for a well producing against a constant bottom-hole well-flowing pressure was given by [34]:

qD=1uK0(u)E48

The solution for a bounded reservoir was presented by [5]:

q¯D=I1(reDu)K1(u)K1(reDu)I1(u)u[I0(u)K1(reDu)+K0(u)I1(reDu)]E49

For considerable longer times, Ref. [27] showed that the qD function in Eq. (48) may be approximated by:

1qD=12[lntD+0.80907]E50

where the dimensionless reciprocal rate and reciprocal rate derivative are given by:

1/qD=khΔm(P)1422.52TqE51
tD*(1/qD)'=kh[t*Δm(P)']1422.52TqE52

Including pseudoskin effects in Eq. (49),

1qD=12[lntD+0.80907+2s']E53

3.2. Conventional analysis

After replacing the dimensionless quantities and changing the logarithm base, it yields:

1q=1.422×106qTkhΔm(P)[1.1513log(kta(P)1688φrw2)+s']E54

As for the case of pressure transient analysis, from a semilog plot of pseudopressure versus time (or pseudotime), its slope allows calculating the reservoir permeability and the intercept is used to find the pseudoskin factor, respectively:

k=1637.74TmhΔm(P)E55
s'=[(1/q)1hrmlog(kφ(μct)irw2)3.227+0.8686]E56

Considering approximation for large time to the analytical Laplace inversion of Eq. (49), the following expression is obtained:

qD=1lnreD0.75exp[2tDreD2(lnreD0.75)]E57

For tDtDpss, this flow period is known as the exponential decline period. tDpss is the time required for the development of true pseudosteady state at the producing well for constant rate production case. Eq. (57) concerns only the circular reservoir. The solution can be generalized for other reservoir shapes by using the Dietz shape factor [6], CA,

qD=2ln(4ADγCA)exp[4πtDADln(4ADγCA)]E58

where, AD (dimensionless area) and reD (dimensionless radius) are given by:

AD=Arw2E59
reD=rerwes=rerweffE60

Eq. (58) suggests that a plot of log(q) versus time will yield a straight line with negative slope Mdecline,

Mdecline=2(0.0002637)kreD2(lnreD0.75)φμctrw2E61

and intercept at (t = 0):

qint=khΔm(P)1637.74Bμ(lnreD0.75)E62

The reservoir area can be determined by solving the Eq. (62) for reD:

reD=exp(1637.74BμkhΔm(P)(lnreD0.75)qint+0.75)E63

3.3. TDS technique

Escobar et al. [9] extended the TDS Technique for gas well in homogeneous and naturally fractured formations using rigorous time. The equations they presented for wellbore storage coefficient and permeability are given below:

C=0.4196TqtNμΔm(P)N=0.4198TμΔm(P)[tt*(1/q)']NE64
k=711.5817ThΔm(P)[t×(1/q)']rE65

Using a procedure similar to the pressure transient case, Escobar et al. [9] found an expression to estimate the pseudoskin factor:

s'=0.5{(1/q)r[t×(1/q)']r-ln(ktrφμctrw2)+7.43}E66

For the estimation of reservoir area, Escobar et al. [9] also presented an equation that uses the starting time of the pseudosteady-state period, tspss.

re2=(0.0015ktspssφμct)1/2E67

As treated in pressure transient analysis, Eq. (41), the reciprocal rate derivative takes a value of 0.5 during radial flow. The intercept of this with the reciprocal rate derivative of Eq. (57) will provide:

tDrpi=12reD2[ln(reD)-0.75]E68

in which numerical solution gives:

reD=1.0292tDrpi0.4627E69

After replacing the dimensionless quantities, we obtain:

re=22.727×10-3rweff(kφμctrweff2)0.4627trpi0.4627E70

Refs. [13] and [14] presented rate transient analysis for long homogeneous and naturally fractured oil reservoirs using TDS technique and conventional analysis, respectively. Equations can be easily translated to gas flow.

3.4. Example 2

Escobar et al. [9] presented an example for a homogeneous bounded reservoir. Figure 5 and Table 2 present the reciprocal rate and reciprocal rate derivative versus rigorous time for this exercise. Other relevant data for this example are given below:

Figure 5.

Reciprocal rate and reciprocal rate derivative for example 2—homogeneous bounded reservoir. After Ref. [9].

h = 80 ftk = 25 mdrw = 0.3 ft
φ = 25%ct = 0.00187 psi−1μg = 0.0122 md/cp
γg = 0.85T = 670 Rre = 30 ft (0.065 Ac)
Δm(P) = 340920304.2 psi2/cpΔP = 580 psi

t, hr1/q, MSCF1t*(1/q)′, MSCF1t, hr1/q, MSCF1t*(1/q)′, MSCF1
3.11E-0443189.24959179751.10714.72E-0217360.95346135561.954
3.89E-0440834.91569162199.43976.09E-0216811.48278135417.1412
4.67E-0439060.52002163534.76457.58E-0216365.45213135450.8332
5.45E-0437655.94848160065.48770.09476764315935.44127134654.7122
7.01E-0435537.19664157149.87770.12217558815464.07533130033.9989
8.72E-0433849.9472154708.18150.15207516515058.4209120446.6934
1.07E-0332361.46026152593.87760.17699147914766.60107110524.7186
1.35E-0330871.32659150543.04990.22483080214266.9873291414.11215
1.66E-0329601.36073148839.15110.28462995613699.3497871972.81156
1.97E-0328617.55908147555.50610.35439563513079.4416855862.58323
2.36E-0327646.46771146325.50090.45007428112279.6627141400.85561
2.80E-0326785.03603145257.08680.56967258811348.6554530197.97654
3.42E-0325823.21027144105.05980.70920394610349.6439222079.63438
4.04E-0325068.84401143224.78960.9005612389119.4179715300.5905
4.73E-0324398.17164142461.49541.1397578527783.90359110316.16432
5.54E-0323753.19888141745.66071.3151687026930.7213057965.441894
6.63E-0323056.85049140996.52141.5224724356043.0905946009.327954
7.87E-0322423.58559140330.6251.8015351525027.2929964239.056471
9.12E-0321908.28203139803.03532.1204639714077.3216322935.619806
1.15E-0221141.51732139045.84552.6307500812923.3109681712.191843
1.49E-0220332.18915138280.54263.284554161917.14852909.5234175
1.86E-0219683.45322137688.68614.2413406181040.993708385.7970016
2.34E-0219065.59817137131.47595.261912839543.216581162.6885651
3.02E-0218406.42719136508.25326.569520997230.482916841.88651512
3.77E-0217872.96521135981.94479.1209515527.980820760.548915189

Table 2.

Reciprocal rate, reciprocal rate derivative versus time data for example 2.

Find reservoir permeability, skin factor, and drainage radius for this example using the TDS Technique.

3.4.1. Solution

The following characteristic points were read from Figure 5:

tr = 4.72 × 10−2 hr[t*(1/q)’]r = 7.293 × 10−6 D/Mscf
(1/q)r = 5.76 × 10−5 D/Mscftrpi = 0.06 hr

Eqs. (65), (66), and (70) are used to obtain permeability, skin factor, and drainage.

k=711.5817ThΔm(P)[t×(1/q)']r=711.5817(670)(80)(30976300)(7.293×106)=26.37mdE8300
s'=0.5{(5.76×105)(7.293×106)-ln((25)(0.0472)(0.25)(0.0122)(0.00187)(0.3)2)+7.43}=0.68E8301
re=22.727×10-3(0.3)((25)(0.25)(0.0122)(0.00187)(0.3)2)0.4627(0.6)0.4627=19.5ftE8003

Notice that the results closely match the permeability and external reservoir radius as presented by Ref. [9].

Finally, it is worth to mention that nowadays, conventional shale-gas reservoirs have become very attractive in the oil industry. Then, their characterization via well test analysis is very important. Shale-gas reservoir is normally tested under constant well-flowing pressure conditions—transient rate analysis—then, the recent studies performed in Refs. [17] and [22] should be read. If such wells are tested under constant rate conditions—pressure transient analysis—then the reader should refer to the works by Bernal et al. [3] and Escobar et al. [18].

Advertisement

Nomenclature

AWell drainage area, ft2 and Ac
BVolumetric factor, rb/MSCF
CWellbore storage coefficient, bbl/psi
ctTotal compressibility, 1/psi
DTurbulent flow factor, Mscf/D
hFormation thickness, ft
hpPerforated interval, ft
I0, I1Bessel function
kPermeability, md
K0,K1Bessel function
mSemilog slope
m*Cartesian slope
m(P)Pseudopropressure function, psi2/cp
MdeclineSlope of plot of log(q) versus time
nFlow exponent
PPressure, psi
PDDimensionless pressure
PwfWell-flowing pressure, psi
qGas flow rate, MSCF
1/qReciprocal of the flow rate, D/Mscf
rRadius, ft
reExternal reservoir radius, ft
rwRadio del pozo, ft
rweffEffective wellbore radius, rwe−s, ft
s’Apparent or pseudoskin factor
saTotal skin factor
tTime, hr
tD*PD’Dimensionless pressure derivative
tDpssExponential decline period
t*(1/q)’Reciprocal rate derivative, D/Mscf
tD*(1/qD)’Dimensionless reciprocal rate derivative
tpHorner or producing time
tpssExponential decline period, hr
tspssTime to initiate pseudosteady state, hr
uArgument for a Bessel function
ZGas supercompressibility factor
αTurbulence factor or inertial factor
ΔChange, drop
φPorosity, fraction
γEuler’s constant—1.781 or e0.5772
γgGas gravity
λMobility, md/cp
μViscosity, cp
1 hrOne hour
crCondition at critical point
DADimensionless referred to drainage area
DaDimensionless referred to pseudotime
DDimensionless
DeDimensionless referred to external
eExternal
effEffective
gGas
iInitial or intercept
pssPseudosteady state
rRadial flow
refReference
rpiIntercept radial-pseudosteady
tTotal
ta(P)Pseudotime, psi-hr/cp
wWell

References

  1. 1. Agarwal, G., 1979. Real gas pseudo-time a new function for pressure buildup analysis of MHF gas wells. In 54th Technical Conference and Exhibition of the Society of Petroleum Engineers of AIME held in Las Vegas, NV, 23–26 September 1973.
  2. 2. Al-Hussainy, R., Ramey, H.J. Jr., and Crawford, P.B. 1966. The flow of real gases through porous media. Journal of Petroleum Technology, Transactions AIME. 18:624–636.
  3. 3. Bernal, K.M., Escobar, F.H., and Ghisays-Ruiz, A. 2014. Pressure and pressure derivative analysis for hydraulically-fractured shale formations using the concept of induced permeability field. Journal of Engineering and Applied Sciences. 9(10):1952–1958. ISSN 1819-6608.
  4. 4. Chaudhry, A.U. 2003. Gas well testing handbook. Gulf Professional Publishing, Burlington, MA, USA, 887 p.
  5. 5. Da Prat, G., Cinco-Ley, H., and Ramey, H. 1981, June 1. Decline curve analysis using type curves for two-porosity systems. Society of Petroleum Engineers. 21:354–362. doi:10.2118/9292-PA.
  6. 6. Earlougher, R. C. 1971,October 1. Estimating drainage shapes from reservoir limit tests. Society of Petroleum Engineers. 23. pp. 1266–1268. doi:10.2118/3357-PA.
  7. 7. Escobar, F.H., Lopez, A.M., and Cantillo, J.H. 2007, Decemeber. Effect of the pseudotime function on gas reservoir drainage area determination. CT&F—Ciencia, Tecnología and Futuro. 3(3):113–124. ISSN 0122-5383.
  8. 8. Escobar, F.H., Hernández, Y.A., and Hernández, C.M. 2007. Pressure transient analysis for long homogeneous reservoirs using TDS technique. Journal of Petroleum Science and Engineering. 58(1–2):68–82. ISSN 0920-4105.
  9. 9. Escobar, F.H., Sanchez, J.A., and Cantillo, J.H. 2008, December. Rate transient analysis for homogeneous and heterogeneous gas reservoirs using the TDS technique. CT&F—Ciencia, Tecnología y Futuro. 4(4):45–59.
  10. 10. Escobar, F.H., Hernandez, Y.A., and Tiab, D. 2010, June.Determination of reservoir drainage area for constant-pressure systems using well test data. CT&F—Ciencia, Tecnología y Futuro. 4(1):51–72. ISSN 0122-5383.
  11. 11. Escobar, F.H., Muñoz, Y.E.M., and Cerquera, W.M. 2011, September. Pressure and pressure derivate analysis vs. pseudotime for a horizontal gas well in a naturally fractured reservoir using the TDS technique. Entornos Journal. Issue (24):39–54.
  12. 12. Escobar, F.H., Martinez, L.Y., Méndez, L.J., and Bonilla, L.F. 2012, March. Pseudotime application to hydraulically fractured vertical gas wells and heterogeneous gas reservoirs using the TDS technique. Journal of Engineering and Applied Sciences. 7(3):260–271.
  13. 13. Escobar, F.H., Rojas, M.M., and Bonilla, L.F. 2012, March. Transient-rate analysis for long homogeneous and naturally fractured reservoir by the TDS technique. Journal of Engineering and Applied Sciences. 7(3):353–370. ISSN 1819-6608.
  14. 14. Escobar, F.H., Rojas, M.M., and Cantillo, J.H. 2012, April. Straight-line conventional transient rate analysis for long homogeneous and heterogeneous reservoirs. Dyna. 79(172):153–163. ISSN 0012-7353.
  15. 15. Escobar, F.H., Muñoz, Y.E.M., and Cerquera, W.M. 2012. Pseudotime function effect on reservoir width determination in homogeneous and naturally fractured gas reservoir drained by horizontal wells. Entornos Journal. Issue (24):221–231.
  16. 16. Escobar, F.H., Zhao, Y.L., and Zhang, L.H. 2014. Interpretation of pressure tests in hydraulically-fractured wells in bi-zonal gas reservoirs. Ingeniería e Investigación. 34(4):76–84.
  17. 17. Escobar, F.H., Montenegro, L.M., and Bernal, K.M. 2014. Transient-rate analysis for hydraulically-fractured gas shale wells using the concept of induced permeability field. Journal of Engineering and Applied Sciences. 9(8):1244–1254.
  18. 18. Escobar, F.H., Bernal, K.M., and Olaya-Marin, G. 2014, August. Pressure and pressure derivative analysis for fractured horizontal wells in unconventional shale reservoirs using dual-porosity models in the stimulated reservoir volume. Journal of Engineering and Applied Sciences. 9(12):2650–2669.
  19. 19. Escobar, F.H., Castro, J.R. and Mosquera, J.S. 2014, May. Rate-Transient Analysis for Hydraulically Fractured Vertical Oil and Gas Wells. Journal of Engineering and Applied Sciences. 9(5):739–749. ISSN 1819-6608.
  20. 20. Escobar, F.H., Ghisays-Ruiz, A. and Bonilla, L.F. 2014, September. New Model for Elliptical Flow Regime in Hydraulically-Fractured Vertical Wells in Homogeneous and Naturally-Fractured Systems. Journal of Engineering and Applied Sciences. 9(9):1629–1636. ISSN 1819-6608.
  21. 21. Escobar, F.H., Zhao, Y.L. and Fahes, M. 2015, July. Characterization of the naturally fractured reservoir parameters in infinite-conductivity hydraulically-fractured vertical wells by transient pressure analysis. Journal of Engineering and Applied Sciences. 10(12):5352–5362. ISSN 1819-6608.
  22. 22. Escobar, F.H., Rojas, J.D., and Ghisays-Ruiz, A. 2015, January. Transient-rate analysis hydraulically-fractured horizontal wells in naturally-fractured shale gas reservoirs. Journal of Engineering and Applied Sciences. 10(1):102–114. ISSN 1819-6608.
  23. 23. Escobar, F.H., Pabón, O.D., Cortes, N.M., Hernández, C.M. 2016, September. Rate-transient analysis for off-centered horizontal wells in homogeneous anisotropic hydrocarbon reservoirs with closed and open boundaries. Journal of Engineering and Applied Sciences. 11(17):10470–10486. ISSN 1819-6608.
  24. 24. Escobar, F.H., Cortes, N.M., Pabón, O.D., Hernández, C.M. 2016, September. Pressure-transient analysis for off-centered horizontal wells in homogeneous anisotropic reservoirs with closed and open boundaries. Journal of Engineering and Applied Sciences. 11(17):10156–10171. ISSN 1819-6608.
  25. 25. Fligelman, H., Cinco-Ley, H., and Ramey, H.J., Jr. 1981, March. Drawdown testing for high velocity gas flow. Paper SPE 9044 presented at the 1981 California Regional Meeting in Bakersfield, CA.
  26. 26. Geertsma, J. 1974, October 1. Estimating the coefficient of inertial resistance in fluid flow through porous media. Society of Petroleum Engineers Journal. 14:445–450. doi:10.2118/4706-PA.
  27. 27. Jacob, C.E., and Lohman, S. W. 1952, August. Non-steady flow to a well of constant drawdown in an extensive aquifer. Transactions American Geophysical Union. 559–569.
  28. 28. Jones, P. 1962, June 1. Reservoir limit test on gas wells. Journal of Petroleum Technology. 14:613–619. doi:10.2118/24-PA.
  29. 29. Lu, J., Li, S., Rahma, M.M. and Escobar, F.H. Escobar, F.H. and Zhang, C.P. 2016, August. Production Performance of Horizontal Gas Wells Associated with Non-Darcy Flow. 11(15):9428–9435. ISSN 1819-6608.
  30. 30. Nunez, W., Tiab, D., and Escobar, F.H. 2003, January,1. Transient pressure analysis for a vertical gas well intersected by a finite-conductivity fracture. Society of Petroleum Engineers. doi:10.2118/80915-MS.
  31. 31. Moncada, K., Tiab, D., Escobar, F.H., Montealegre-M, M., Chacon, A., Zamora, R.A., and Nese, S.L. 2005, December. Determination of vertical and horizontal permeabilities for vertical oil and gas wells with partial completion and partial penetration using pressure and pressure derivative plots without type-curve matching. CT&F—Ciencia, Tecnología y Futuro. 2(6):77–95.
  32. 32. Ramey, H.J., Jr. 1965, February. Non-Darcy flow and wellbore storage effects in pressure buildup and drawdown of gas wells. Journal of Petroleum Technology. 17:223.
  33. 33. Tiab, D. 1995. Analysis of pressure and pressure derivative without type-curve matching: 1-skin and wellbore storage. Journal of Petroleum Science and Engineering. 1995;12:171–181.
  34. 34. Van Everdingen, A.F., and Hurst, W. 1949, December. The application of the Laplace transformation to flow problems in reservoirs. Society of Petroleum Engineers. 1. doi:10.2118/949305-G.
  35. 35. Zhao, Y.L. Escobar, F.H., Hernandez, C.M., Zhang, C.P. 2016, August. Performance analysis of a vertical well with a finite-conductivity fracture in gas composite reservoirs. Journal of Engineering and Applied Sciences. 11(15):8992–9003. ISSN 1819-6608.

Written By

Freddy Humberto Escobar

Submitted: 22 September 2015 Reviewed: 26 January 2017 Published: 02 August 2017