CC BY
84
DOI: https://doi.Org/10.18599/grs.2018.1.39-43
The results interpretation of thermogasdynamic studies of vertical gas wells incomplete in terms of the reservoir
penetration degree
M.N. Shamsiev
Institute of Mechanics and Engineering, FRC Kazan Science Center of the Russian Academy of Sciences, Kazan, Russian Federation
E-mail: mshamsiev@imm.knc.ru
Abstract. A method is proposed for interpreting thermogasdynamic studies of vertical gas wells that are incomplete in terms of the reservoir penetration degree on the basis of inverse problem theory. The inverse problem has the aim to determine the reservoir parameters for nonisothermal filtration of a real gas to a vertical well in an anisotropic reservoir. In this case, the values of the pressure and temperature at the well bottom, recorded by deep instruments, are assumed to be known. The solution of the inverse task is to minimize the functional. The iterative sequence for minimizing the functional is based on the LevenbergMarquardt method. The convergence and stability of the iterative process for various input information have been studied on specific model examples. The effect of reservoir anisotropy on the pressure and temperature changes at the bottom of the well is studied. It is shown that if the reservoir is not completely penetrated by the results of pressure and temperature measurements at the bottom of the well after its start-up, anisotropy of the reservoir can be estimated. It should be noted that when studying thermodynamic processes in the vicinity of a well, which penetrates thick layers, it is necessary to take into account not only the heat exchange of the reservoir with the surrounding rocks, but also the geothermal temperature gradient.
Keywords: anisotropy, thermogasdynamic studies, incomplete well
Recommended citation: Shamsiev M.N. (2018). The results interpretation of thermogasdynamic studies of vertical gas wells incomplete in terms of the reservoir penetration degree. Georesursy = Georesources, 20(1), pp. 39-43. DOI: https://doi.org/10.18599/grs.2018.L39-43
Introduction
Most of the gas fields have a layered structure, due to the peculiarities of the sedimentation process. In layered reservoirs, the filtration properties in the plane of the layers differ from the filtration properties in the direction perpendicular to the layers, i.e. such reservoirs are anisotropic. In real gas reservoirs, anisotropy can be caused by fracturing, layering, the presence of various inclusions, which lead to unequal properties of the medium in different directions.
The heterogeneity of the formation in the vertical and horizontal directions is characterized by the anisotropy parameter. The parameter of reservoir anisotropy is of decisive importance in predicting the technological regime of wells operation, exposing reservoirs with bottom water, multi-layer deposits, etc. In work (Gritsenko et al., 1995), a grapho-analytical method is suggested for estimating anisotropy based on the results of gas-hydrodynamic studies of vertical wells that are incomplete in terms of the reservoir penetration degree. The incompleteness of the bottom entails the
© 2018 M.N. Shamsiev.Published by Georesursy LLC This is an open access article under the CC BY 4.0 license (https://creativecommons.org/licenses/by/4.0/)
appearance of additional filtering resistances arising in the bottomhole zone and a decrease in production rates as a result of the deviation in the liquid flow geometry from the flat-radial flow (Basniev et al., 1993; Gritsenko et al., 1995; Korotaev et al., 1991; Masket, 1949; Shchurov, 1983). In this regard, consideration of the characteristics of the inflow to incomplete wells is of great practical importance.
In this paper, we propose a method for interpreting thermogasdynamic studies of vertical gas wells that are incomplete in terms of the reservoir penetration degree on the basis of the inverse problem theory. The effect of anisotropy on the pressure and temperature changes at the bottom of the well is investigated. It is shown that the results of pressure and temperature measurements at the bottom of the well after its start-up can be used to evaluate the anisotropy of the reservoir.
Nonisothermal gas filtration in an anisotropic reservoir
Thermogasdynamic methods for studying gas wells are based on the processes of pressure and temperature redistribution in the reservoir when they are put into operation and after a stop. The process of non-isothermal
¡SCIENTIFIC AND TECHNICAL JOURNAL
www.geors.ru GEDRESURSY
filtration of a real gas to a vertical well in an anisotropic reservoir is described by a system of differential equations:
* 11. r dr
B
m— dt
r ^ P_
Tç.
rK p dp^
fj.Tç dr ) dz\/j. Tç dz
re{rw,Rk), z e (0,L), t> 0,
(1)
dz
+
\dp 'dT -s^l ,Kdp 'dT dp) 'dt)
dr _ dr c dr_ T H dz _dz G dz_
„ dT 1 d ( . &T\ d (, dT Q — =--— + A
dt r dry dr J dz
p RTq
rs(rw>Rk), zs(0,L), t> 0, with initial
p(r,z,0) = p0, rz{rw,Rk),ze{0,L), (3)
T{r,z,0) = T0+(L-z)G, re (rw,Rk ),ze (0,L), (4) and boundary conditions:
(2)
k PTst rdP Cw S H Ps,Tc, dr Lw dt
2n J
L-K
z e [L-Lw,l], 1 dt p RTç
f \
P_
yTÇ;
dz = Q,.
(5)
' k dp 'dT -aÔ-P\ \
Kßdr Jir G dr /
mrj
dp dr
pRTç
= 0, z e [0,L-Lw),
0,ze[L-Lw,L],
dT dr
= 0, ze[0,L-Lw),
p(rk,z,t) = p0, z e [0,Z,], t> 0, T(rk,z,t) = T0 + (L-z)G, z<=[0,L], t> 0,
= 0, r e (rw,Rk), t>0,
dp 'dp
dz. z=0 j z=L
(6)
(7)
(8)
(9) (10)
(11)
dz
df dz
z=0 =L
z=0 z=L
V 7lt ^ T\Zz=L+Kf) 0, re(rw,4 t> 0
HHH
(12)
where L - the thickness of the reservoir, Lw - the depth of the reservoir penetration, m - the gas viscosity, m - the reservoir porosity, p0, T0 - the pressure and temperature at the reservoir boundary, pst, Tst - the standard pressure and temperature, rw, Rk - the radii of the well and external boundary, Q - the well production rate, z - the gas supercompressebility, k, kz - the reservoir permeability along the directions of r and z axes, n - the adiabatic expansion factor, s - the Joule-Thomson coefficient, Cl=mCpp/RTç + {l-m)C2p2 - the volumetric heat
capacity of the reservoir, C C2 - the specific heat capacity of the gas and medium, p2 - the medium density, R - the gas constant, 12 - the thermal conductivity of the gas and medium, Cw - wellbore storage coefficient, G - geothermal gradient. The first expression of condition (12) simulates the process of heat exchange on the roof and the bottom of the gas reservoir by rocks, i.e. non-thermal insulated reservoir. Conditions (5) -(8) simulate the process of heat and mass transfer to a vertical well that is incomplete in the degree of reservoir penetration.
Zone of thermal compensation above the roof and below the bottom of the formation is calculated by the formula:
K-^t
C2p2
The wellbore storage coefficient is calculated by the formula (Korotaev et al., 1991; Khairullin et al., 2013, Shamsiev, Talipova, 2015). The gas supercompressibility coefficient Z is calculated by the formula of Gurevich-Latonov (Bondarev et al., 1988). The adiabatic expansion coefficient n and the JouleThomson coefficient e are calculated by the formulas (Korotaev et al., 1991).
Equations (1) - (12) belong to the class of quasilinear parabolic equations. The greatest difficulty is the numerical solution of equation (2), which simultaneously describes conductive and convective heat transfer, including the Joule-Thomson effect, as well as a decrease in the gas temperature due to its adiabatic expansion. For the numerical solution of the system (1) - (12), the finite difference method is applied. The solution area is covered by a non-uniform grid, which thickens at the approach to the well.
Results of the calculations
In the case of model examples, influence the reservoir anisotropy on pressure and temperature curves at the bottom of the well is studied. A model reservoir with the following data is considered: H = 20 m, R, = 250 m, r = 0.1 m, pn = 20 MPa, Tn = 300 K,
k ' w ' r 0 ' 0
Tst = 293 K, m = 0.012 mPa s, m = 0.2, kr = 0.01 mkm2, Cp =2093 J/ kg K, C2 = 920 J/ kg K, r2 = 2700 kg / m3, RP= 520 J/kg K, 11 = 1.52 W/m K, 12 = 1.9 W/m K, G = 0.01 K/m, Q = 500 thousand m3/day, operation time of the well t = 5 days.
exp J
From the results obtained, it follows that the process of heat exchange of the reservoir with surrounding rocks affects the temperature field and has little effect on the pressure field in the reservoir. It should be noted that when studying thermodynamic processes in the vicinity of a well, which penetrates thick layers, it is necessary to take into account not only the heat exchange of the reservoir with the surrounding rocks, but also the geothermal
SCIENTl FK AND TECHNICAL JOURNAL
'4M GEDRESDURCES www.geors.ru
temperature gradient. Comparison of the numerical solution of the system (1) - (12) with the numerical solution of the one-dimensional task (Shamsiev, Talipova, 2015) at G = 0, kr = kz is shown in Figure 1. Deviations of the final sections of the derivatives of the bottomhole temperature curves are observed. This is due to the process of heat exchange with surrounding rocks. In the case of a thermally insulated reservoir, numerical solutions of one-dimensional and two-dimensional tasks almost coincide.
Figure 2 shows curves of bottomhole pressure, temperature and their derivatives depending on the
Figure 1. The bottomhole temperature curve and its derivative. The reservoir is completely penetrated and not thermally insulated, G = 0. k = k. 1 - numerical solution of system (1) - (12), 2 - numerical solution of one-dimensional task.
reservoir anisotropy after the well has been put into operation without full penetration. The calculations were carried out at к= kz = 0.01 цт2; кг = 0.01 цт2, к = 0.02 цт2; кг = 0^.01 Цт2, к = 0.005 цт2. When launching a vertical well with a constant production rate and incornplete penetration of the reservoir, the initial radial, spherical and late radial flows are observed. The initial radial flow is masked by the presence of the influence of the borehole volume (Figure 2). The negative slope of the derivative curve characterizes the presence of a spherical flow. The less the reservoir is penetrated, the longer the duration of the spherical flow. The straight section on the pressure derivative curve is the diagnostic sign of the late radial flow (Figure 2). The deviations of the final sections of the pressure and temperature curves characterize the influence of the reservoir boundary.
The results of calculations showed that the less the reservoir is penetrated, the rnore sensitive the pressure and ternperature curves to variations in the perrneability coefficients к, к.
r z
Evaluation of the reservoir anisotropy
The results interpretation of therrnogasdynarnic studies of gas wells is based on the solution of the inverse problern. As the initial inforrnation, the data on the change in bottornhole pressure and ternperature, recorded by deep instrurnents after the start-up of the well, are used. The airn of the inverse task is to deterrnine the permeability coefficients к, кг and the porosity of the reservoir m, when the process of non-isothermal filtration of a real gas to a vertical well incornplete in the degree of reservoir penetration is described by the systern of equations (1) - (12).
Figure 2. Curves of pressure (a), temperature (b) and their derivatives. The reservoir is penetrated by 50%. 1 - k =k, 2 - k<k, 3 - k>k
r z r z
SCIENTIFIC AND TECHNICAL JOURNAL
www.geors.ru GEDRESURSY
The following initial information is considered known:
p{rw,L,t) = tft), T(rw,L,t) = ¥{t), (13)
where f(t), y(t) - observed values of pressure and temperature at the bottom of the well.
The solution of the inverse problem (1) - (12) and
(13) is to minimization the functional (Khairullin et al., 2013, Shamsiev, Badertdinova, 2012; Shamsiev, Talipova, 2015):
F(o)= ]^f)-p{rw,L,tf +[¥{t)-T{rw,L,t)f}dtl (14) o
where a = (k kz, m), 0 < a.< a.< b. (a,, b = const), X -weighting parameter.
The iterative sequence for minimizing the functional
(14) is based on the Levenberg-Marquardt method. The convergence and stability of the iterative process with different input information have been studied using standard madel examples. The iterative process is considered to be completed when one of the specified precision is reached (10-6 - for the functional, 10-6 -for the gradient, 10-6 - for the argument) or when performing a specified number of iterations (N.er = 40). For exact values of the initial information, the iterative process of minimizing the functional (14) converges in 6-8 iterations. To study the stability in the model curves, changes in bottomhole pressure and temperature randomly introduced errors, where </>Si{t) = </>{t)+ ad^ , vs2{t) = v{t)+ (oS2, where §1 = 0.05 Mna, S2= 0.05 K, co - a random variable distributed along a uniform law on the interval [-1, 1]. With the perturbed initial data, the iterative process of minimizing the functional (14) converges for 10-15 iterations. The results of the calculations show that the proposed method is stable with respect to the errors of the initial information.
Figure 3 shows one of the characteristic calculations of the convergence of the iterative process of minimizing the functional (13) with the perturbed initial data
Figure 3. Convergence of the iterative process. 1 - к, 2 - к, 3 - m
Figure 4. Curves of temperature changes. The reservoir is penetrated by 50%. 1 - perturbed, 2 - calculated curves
(Figure 4, curve-1), where a - true pararneters, a abs -calculated pararneters. The calculated ternperature curves are shown in Figure 4 (curve-2). The iterative process converges in 12 iterations.
Frorn the results obtained, it follows that when the reservoir is partially penetrated by the results of therrnogasdynarnic studies of gas wells, anisotropy of the reservoir can be estirnated.
References
Basniev K.S., Kochina I.N., Maksimov V.M. (1993). Podzemnaya gidrornekhanika [Underground hydrornechanics]. Moscow: Nedra, 303 p. (In Russ.)
Bondarev E.A., Vasil'ev V.I., Voevodin A.F., Pavlov N.N., Shadrina A.P. (1988). Terrnogidrodinarnika sistern dobychi i transporta gaza [Therrnohydrodynarnics of gas production and transportation systerns]. Moscow: Nauka, 270 p. (In Russ.)
Gritsenko A.I., Aliev Z.S., Errnilov O.M., Rernizov V.V., Zotov G.A. (1995). Rukovodstvo po issledovaniyu skvazhin [Guidelines for the study of wells]. Moscow: Nedra, 523 p. (In Russ.)
Khayrullin M.Kh., Shamsiev M.N., Kazunin D.V. (2013). Interpretatsiya rezul'tatov gazogidrodinarnicheskikh issledovaniy vertikal'nykh skvazhin s uchetorn vliyaniya ob'erna stvola [Interpretation of the results of gas-hydrodynamic studies of vertical wells taking into account the influence of the borehole volurne]. Gazovaya promyshlennost = Gas Industry, 2, pp. 16-17. (In Russ.)
Khayrullin M.Kh., Shamsiev M.N., Morozov P.E., Abdullin A.I., Salirn'yanov I.T., Gadil'shina V.T. (2013). Chislennoe reshenie pryarnykh i obratnykh zadach teplornassoperenosa v neftyanykh plastakh [Nurnerical solution of direct and inverse problerns of heat and rnass transfer in oil reservoirs]. VestnikKazanskogo tekhnologicheskogo universiteta [Bulletin of Kazan Technological University], 24(16), pp. 125-129. (In Russ.)
Korotaev Yu.P., Krivoshein B.L., Novakovskiy V.N. (1991). Termodinamika gazopromyslovykh sistem [Thermodynamics of gas field systerns]. Moscow: Nedra, 275 p. (In Russ.)
Masket M. (1949). Techenie odnorodnoy zhidkosti v poristoy srede [Homogeneous fluid flow in a porous medium]. Moscow: Gostoptekhizdat, 628 p. (In Russ.)
Shamsiev M.N., Badertdinova E.R. (2012). Otsenka fil'tratsionnykh i teplofizicheskikh parametrov neftyanogo plasta po rezul'tatam izmereniy temperatury na zaboe skvazhiny [Evaluation of the filtration and thermophysical pararneters of the oil reservoir using bottorn hole ternperature rneasurernents].
Izmeritel'naya tekhnika = Measurement techniques, 3, pp. 45-47. (In Russ.)
Sharnsiev M.N., Talipova A.A. (2015). Terrnogidrodinarnicheskie issledovaniya vertikal'nykh gazovykh skvazhin [Therrnohydrodynarnic studies of vertical gas wells]. Geologiya, geofizika i razrabotka neftyanykh i gazovykh mestorozhdeniy = Geology, geophysics and development of oil and gas fields, 6, pp. 43-46. (In Russ.)
Shchurov V.I. (1983). Tekhnologiya i tekhnika dobychi nefti [Technology and technical equiprnent for oil production]. Moscow: Nedra, 510 p. (In Russ.)
SCIENTIFIC AND TECHNICAL JOURNAL
GEDRESDURCES
www.geors.ru
About the Author
Marat N. Shamsiev - DSc (Engineering), Leading Researcher, Laboratory of Underground Hydrodynamics
Institute of Mechanics and Engineering, FRC Kazan Scientific Center of the Russian Academy of Sciences
2/31, Lobachevsky St., Kazan, 420111, Russian Federation
Tel: +7 843 2929345
E-mail: mshamsiev@imm.knc.ru
Manuscript received 18 December 2017;
Accepted 27February 2018; Published 30 November 2017
SCIENTIFIC AND TECHNICAL JOURNAL
www.geors.ru GEDRESURSY M.
