УДК 532.546
ГИДРОДИНАМИЧЕСКОЕ ИССЛЕДОВАНИЕ СКВАЖИН C МАГИСТРАЛЬНОЙ ТЕХНОГЕННОЙ ТРЕЩИНОЙ ГИДРОРАЗРЫВА ПЛАСТА: ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ И АНАЛИЗ ПРОМЫСЛОВЫХ ДАННЫХ
INTERFERENCE TEST TO FRACTURED INJECTION WELLS: MATHEMATICAL MODEL AND FIELD CASE
Г. Ф. Асалхузина, А. Я. Давлетбаев, Р. И. Нуриев
G. F. Asalkhuzina, A. Ya. Davletbaev, R. I. Nuriev
Башкирский государственный университет, г. Уфа
ООО «ПФК «ГИС-Нефтесервис», г. Уфа
Ключевые слова: гидропрослушивание; математическое моделирование Key words: interference test; mathematical modeling
Hydraulic fracturing is the state-of-the-art development technique to optimize the productivity of low permeability hydrocarbon reservoirs and is performed both in production and injection wells. Depending on the proppant volumes placed during a hydraulic fracturing operation, fracture half-lengths of 150-200 m and even further are desired in order to tap oil and gas resource plays. Development of hydraulic fractures mainly occurs along the axes of the regional stress field assuming fracturing operations are conducted in areas lacking a significant alteration of pressure and temperature, i.e. new exploration wells [1].
Modern pumping units which are employed for hydraulic fracturing operations achieve wellhead pressures of up to 20 MPa during the injection of fracturing fluid into the reservoir. Spontaneous hydraulic-fracturing growth in the horizon of interest is initiated as soon as injection pressure exceeds yield strength of the particular reservoir rock [2-8]. Real-time data acquisition of wellhead pressures from active and observation wells and the results of interference testing confirm that in some cases hydraulic fractures may achieve a length of up to 1 000 m [8].
S. Ekie et al. [9] present a mathematical model for the fracture formation for two vertical wells — one active well and one observation well. Conducting pulse-tests for the tracking of registration response pressures in several vertical observation wells can be used to determine the orientation of the fracture originating from the active fracture well. N. Mousli et al. [10], D. Meehan et al. [11] presented analytical solutions of modeling interference test for wells with parallel. D. Tiab and E. Abobise [12] discussed the design of pulse testing of a fractured active well and vertical observation well (without fractures), presenting correlations for fracture orientation, average permeability and other parameters of the reservoir on the response amplitude and the response time in the observation well. In the same year K. Cooper and R. Collins [13] applied of a mathematical model which is based on interference test of a two-well system with a fractured (one fracture) and another non-fracture well that was used to interpret the results of a field case and to estimate underlying parameters of the system. The article of H. Najurieta et al. [14] presents a two-dimensional model for heterogeneous anisotropic reservoirs; demonstrating the mathematical modeling for all wells. It can be used to calculate 2D pressure maps, transmissivity and diffusivity between active and observation wells estimated. Reservoir case studies with conductive faults were also presented in the article [15].
Contrary to the abovementioned papers we elaborate a case that consists of the active fractured well and the observation fractured well (fractures are oriented along the regional stress) either having a low-permeability porous bridge between the fracture tips, or two wells which are being directly connected through intersecting fracture planes without a porous media bridge. Furthermore, we want to highlight the results of interfe-
rence test from the work [8] by interpreting the chosen data model and providing an estimation of the fracture and reservoir characteristics. Formulating the abovementioned problem, injection pressure is assumed to surpass the fracture closure pressure, i. e. its geometry (height, length, width) is treated constant.
Mathematical model. In our approach we consider a symmetric geometry of two hydraulic fractured vertical wells each with one fracture as shown in Figure 1.
The wells are located in the same row and have same orientation along the lines of the regional stress. Fractures are symmetrical and parallel to the Ox axis, xfj and xf2
are the length of the active well and the observation well respectively, both fractures have same width wf and permeability kf. L is the distance between the boundaries of
layer and fractures. In our model the distance d between the fracture tips acts as a so called porous medium bridge.
Fig. 1. Model geometry with distance d > 0 between the fracture tips We assume that the reservoir is isotropic, i. e. k = k , and formation thickness h
* > mx my '
is constant, coinciding with the fracture height. Top and bottom of the reservoir are no-flow boundaries. Here x and y is the distance along the x axis and along y axis respectively, p is the reservoir pressure, which is described by the Laplace equation. With respect to the view from the problem's symmetry Ox axis or along the fracture planes, half of the computational domain is considered. Fluid flow in the fractures behaves only one-dimensional. We suppose that viscosity of reservoir liquid differs a little from viscosity of the injected liquid. Wherein injected agent is water.
Distribution of pressure within the fractures varies according to their length and is constant in each vertical section, as described by equation (1)
d2pf <Pf Wf, dpf
dx2
dt
L < x < L + 2Xf 1, L + 2Xf 1 + d < x < L + 2xf 1 + d + 2xf 2,
(1)
0 < y < wf ,
where is in the formation the pressure is given by equation (2)
d2pm , 52 pm _ L ^iCmt dpm
dx2
dy2
к.
dt
(2)
The inflow of liquid at the boundary of reservoir and fractures in the active well (L < x < L + 2xf j, 0 < y < wf /2) described by equation (3)
(
qA = -2 •
( xf i •h) •
k dp
m г m
f 1 И dy
(ь
+(wfj2•h) •
y=wf! 2
к dp m г m + К dPm \
И dx x=L И dx x=L+2 xfi j
(3)
and in the observation well (L + 2xf 1 + d < x < L +2xf1 + 2xf2, 0 < y < Wf/2). described by equation (4)
(
q0 =-2 •
(xf2•h)•
к dp
m г m
И dy
\
+(wfj2•h)'
к dp
m * m
И dx
y=Wfl2
к dp
, m * m
И dx
x=L+2 Xfi + d+2 Xf 2 j j
(4)
Boundary conditions at time = 0 in the «fracture — formation» and at the outer boundary of the reservoir suggest a constant initial pressure
Pm 1=0 Pf\t=0 Pi
'\y=L
' lx=2L+2Xfi +2Xf2+d
(5)
At the boundary of the «fracture — formation» pressures and fluid flow are equalized by equation (6)
p/Lv=pm
<y- /2
kf dpf
И dy
к dp
m г m
И dy
(6)
The condition of symmetry across the Ox axis, i.e. along fractures, is given in equation (7)
dp
dx
dP„
y=0
dx
= 0.
(7)
y=o
The boundary conditions in the wells are determined by setting the pressure or flow rate of injected fluid as follows (equation 8)
pi
dpf
=L+xf i PwA
y=0
^аИ
dx x=l+xf 1 2kfwfh y=0
(8)
Here pm and pf represent the pressure, km and kf are the permeability, (m and < are the porosity, and cmt and cft are the total compressibility of the fracture and
w
№
y
У
formation (matrix) respectively. p. is initial reservoir pressure, pwA is pressure in the
well, QwA is the flow rate at the active well.
Effect of the porous medium bridge between the fractures on the pressure response: d > 0. The system of equations (1)-(4) and the corresponding boundary conditions (5)-(8) for the geometry, as shown in Figure 1, were solved using the finite difference method for the Newton iteration scheme on a non-uniform rectangular difference grid [16]. The accuracy of the approximation was tested against an analytical solution for the case of wells with a single vertical fracture of finite conductivity [17]. An example of calculation was performed using the following parameters of the system:
km = 1-10"15 m2; 0,3 mPa-s; L =500 m; h = 21,23 m; <j,m = 0,17; ^ =0,414;
c = 3,6687-10"9 1/Pa; cf = 9,4845-10"9 1/Pa; d = 1, 3, 5, 10, 30, 100, 200, 500, 800 m;
mt > ft* ' ???????? >
xf 1 = 900 - d m; xf 2=100 m; kfwy -
10-10-12 m2^m;
Pi
= 27^10b Pa.
Subsequently, we began modeling of interference test to study the effect of a porous medium between the fractures. In the active wells we prescribed a variety of cycles: constant pressure on the injection cycle p/|I=L+^ = p, and a rate equal to 0 at the stop
X=L+Xf
y=0
cycle
dp.
dx
= 0 , thereby gaining reproducible responses of the pressure in obser-
y=o xf
vation wells. Figure 2 shows the pressure evolution in the active well (Fig. 2 a) and the pressure response in observation wells for different values of the distance between the fracture tips: d =100, 200, 500, 800 m (Fig. 2 b) and d = 3, 5, 10, 30 m (Fig. 2 c). The outcome shows that for d > 30 m a poor pressure response would be seen in the observation well. Diagnosing the pressure response of field measurements might become a difficult task when data are noisy in case of low accuracy measurement equipment. However, if d < 30 m (Fig. 2 c), a clear signal from of reaction pressure can be obtained. Furthermore, the computed pressure amplitude in the observation well is increasing when the distance between the fracture tips becomes more narrow.
Fig. 2. Pressure change in the active well (a) and pressure response in the observation well (b, c): 1 — d = 100 m;
2 — d = 200 m; 3 — d = 500 m; 4 — d = 800 m (b); 5 — d = 3 m; 6 — d = 5 m; 7 — d = 10 m; 8 — d = 30 m (c). The smaller the distance d, the larger the amplitude in the pressure response become
The pressure profiles along the Ox axis with fractures at d = 800 m (Fig. 3 a), d = 30 m
(Fig. 3 b), d =1 m (Fig. 3 c) at time 1 and 30 days is shown in Figure. 3. In the case of d = 800 m (Fig. 3 a) the pressure is increasing only along the fractures in the active
well, the pressure in observation wells remains equal to the initial reservoir pressure p..
Comparing the profiles with d = 30 m (Fig. 3 b) and d =1 m (Fig. 3 c) we can infer that the greater the distance between the fractures, the greater the pressure loss between the active and the observation wells.
Figure 4 shows the pressure distribution in the «wells — fractures — formation» system in 2D for a range of d = 500 m (Fig. 4 a), d = 30 m (Fig. 4 b) and d = 1 m (Fig. 4 c) at a time of 30 days. It is seen that for d = 500 m high pressure located just around the active well (along the fracture and around the reservoir) and therefore the response of pressure is not diagnosed in the observation well. However, for d = 1 m and d = 30 m high pressure regimes are distributed along the two wells and their attached fracture.
The presence of a small porous medium bridges having low filtration properties leads to a significant loss of pressure response between the active and the observation wells.
The field case which are presented by A. Davletbaev et al. [8] show that pressure is responding almost instantaneous between the active and observation wells. Furthermore, the pressure difference between the wells is much smaller than in the results which are shown in Figures 2-4. Thus, the case of two fractures without the porous medium bridge must be considered for the interpretation of the interference test with d = 0.
The case without the porous medium bridge between the fractures: d = 0. For the interpretation of the field test we have chosen the case of a single fracture intersecting both the active and the observation well. In this geometry the fracture has a length 2xf 1 + 2xf 2 and an absent porous medium bridge between the fractures, i. e.
d = 0. In this system the response time of pressure and differential pressure between the active and observation wells is given by the dimensionless group for fracture conductivity.
Fig. 3. Pressure distribution along the fracture
for d = 800 m (a); d = 30 m (b); d = 1 m (c) at different times: 1 — t = 1 day; 2 — t = 30 days
x, m
Fig. 4. Pressure distribution in the «wells — fractures — formation»
for a variety of d =500 m (a); d =30 m (b); d =1 m (c) at time t = 30 days
F =
1 cd
kfwf
km(Xf 1 + Xf 2 )
f 2 y
In Figure 5 we plot sensitivity curves for different pressures in the active and observation well for the calculations performed with a range of fracture conductivities, i. e. Fcd = 10 (Fig. 5 a), FCD = 25
(Fig. 5 b), Fcd = 100 (Fig. 5 c). The Figure emphasizes that higher fracture conductivities result in a better match for both pressure changes and absolute values of pressures in the active and the observation wells.
Field case interpretation. The mathematical model, which is presented by A. Davletbaev et al. [8], was used to interpret interference-test pressure data for wells pair (active well and observation well) and to acquire refined fracture parameters. History of injected water volumes into the active well was included in the mathematical model. A genetic algorithm is applied to obtain the best combination of the measured pressure and the computed (theoretical) pressure in the observation well. The following model parameters were used in the interpretation of measured data: u = 0,33 mPa-s;L = 200 m; h = 9,4 m; 6f = S = 0,16; c, = c = 9,48-10-
~ ' ' ' ' ' f Tm^'jt mt '
9 1/Pa; d = 0 m; xf j = 846 m; xf 2 = 100 m, the other parameters (km, pt, FCD) were
determined by simulation.
Figure 6 shows the best combination of calculated and measured curves resulting after 100 iterations performed. The best combination was obtained for a matrix rock permeability to the water phase k = 0,293T0-15 m2 (if phase relative permeability of wa-
Fig. 5. Pressure changes in the active (1) and the observation wells (2)
for Fcd = 10 (a); Fcd = 25 (b); Fcd = 100 (c)
ter k = 0,2, then absolute matrix permeability k = 1,465T0-
dimensionless fracture conductivity FCD
p. = 33,13 MPa.
of 58,5, and a reservoir
m ), resulting in a pressure
Fig. 6. Pressure readings in the active and the observation wells during the field interference test (measured and calculated curves)
Conclusions
• By modeling the injection of water into the «well — fracture — formation» geometry we demonstrate effects for several arrangements of a porous medium bridge and fracture conductivity on the distribution of the pressure field and the pressure response between the active and the observation wells. In presence of a low permeability porous medium bridge between the active and the observation wells we found that fast pressure pulses in the observation well cannot be produced, which could been also seen in field tests. Therefore, the field pulse test was simulated by two fractures without porous medium bridge.
• According to the hydraulic fracturing report the initial dimensionless fracture conductivity and half-length in the active well was 14,4, and 130 m respectively. However, the analysis of a field pulse test identifies a dimensionless fracture conductivity of 58,5 and fracture half-length of 946 m. Thus, the injection above the formation fracturing pressure leads to spontaneous growth of fracture length and increase of the conductivity in injection wells.
The authors express their gratitude to Georg Seidl for his help.
References
1. Selection of the optimal development system for deposits with low-permeability reservoirs / V. A. Baikov [et al.]. -Available at: http://www.ogbus.ru/authors/Baikov/Baikov_2.pdf.
2. S^by J., Pj0mdal H. P., Van den Hoek P. Managed Induced Fracturing Improves Waterflood Performance in South Oman // International Petroleum Technology Conference. - Doha, Qafar (November, 2005).
3. Chavez J. C., Carruthers J., McCurdy P. Water Flooding Efficiency in a Scenario of Multiple Induced Fractures, an Applied Geomechanical Study // SPE 97526: SPE International Improved Oil Recovery Conference in Asia Pacific. - Kuala Lumpur, Malaysia (5-6 December, 2005).
4. Van den Hoek P. J., Volchkov D., Burgos G., Masfry R. A. Application of New Fall-Off Test Interpretation Methodology to Fractured Water Injection Wells Offshore Sakhalin // SPE 102304: SPE Russian Oil and Gas Technical Conference and Exhibition. - Moscow, Russia (3-6 October, 2006).
5. Danymic Induced Fractures in Waterfloods and EOR / P. J. Van den Hoek [et al.] // SPE 115204: SPE Russian Oil and Gas Technical Conference and Exhibition. - Moscow, Russia (29-30 October, 2008).
6. Santarelli F. J., Havm0ller O., Naumann M. Geomechanical Aspects of 15 Years Water Injection on a Field Complex: An Analysis of the Past to Plan the Future // SPE 112944: SPE North Africa Technical Conference and Exhibition. -Marrakech, Morocco (12-14 March, 2008).
7. Testing and Analysis of Wells Producing Commingled Layers in Priobskoye Field / A. Davletbaev [et al.] // SPE 117411: SPE Russian Oil & Gas Technical Conference and Exhibition. - Moscow, Russia (28-30 October, 2008).
8. Multi-Layer Steady-State Injection Test with Higher Bottomhole Pressure than the Formation Fracturing Pressure / A. Davletbaev [et al.] // SPE 136199: SPE Russian Oil & Gas Technical Conference and Exhibition. - Moscow, Russia (26-28 October, 2010).
9. Ekie S., Hadinoto N., Raghavan R. Pulse-testing of vertically fractured wells // SPE 6751: 52nd Annual Meeting of the Society of Petroleum Engineers in Denver, CO. - 1977.
10. The influence of vertical fractures intercepting active and observation wells on influence tests / N. A. Mousli [et al.] // Petrol. Eng. Journal. - SPE 9346. -1982. - Vol. 22. - P. 933-944.
11. Meehan D. N., Horne R. N., Ramey Jr. H. J. Interference testing of finite conductivity hydraulically fractured wells // Soc. of Petrol. Eng. - SPE 19784. - 1989. - P. 137-152.
12. Tiab D., Abobise E.O. Determining Fracture Orientation From Pulse Testing // SPE Formation Evaluation. - SPE 11027 - 1989. - Vol. 4, № 3. - P. 459^66.
13. Cooper K. J., Collins R. E. Applications of transient pressure interference tests to fractured and non-fractured injection wells // Soc. of Petrol. Eng. - SPE 19795. - 1989. - P. 327-339.
14. Transmissivity and Diffusivity Mapping From Interference Test Data: A Field Example / H. Najurieta [et al.] // SPE Formation Evaluation .- SPE 22701. -1995. - Vol. 10, № 3. - P. 180-185.
15. Interference testing in reservoirs with conductive faults or fractures / M. Abbaszadeh [et al.] // Soc. of Petrol. Eng. -SPE 38941. - 1997. - P. 153-163.
16. Aziz K., Settari A. Petroleum reservoir simulation // York: Elsevier Applied Science Publishers. - 1979.
17. Cinco-Ley H., Samaniego V. F., Dominguez A. N. Transient pressure behavior for a well with a finite-conductivity vertical fracture // Soc. of Petrol. Eng. - SPE 6014. - 1978 - P. 253-264.
Сведения об авторах
Давлетбаев Альфред Ядгарович, к. ф.-м. н., доцент кафедры прикладной физики, Башкирский государственный университет, г. Уфа, тел. 89174812699, е-mail: davletbaevay@rambler. ru
Асалхузина Гузяль Фаритовна, аспирант кафедры прикладной физики, Башкирский государственный университет, г. Уфа, тел. 89178053644, e-mail: kafedrapf@rambler. ru
Нуриев Рустам Илдусович, начальник отдела, производственно-коммерческая фирма «ГИС-Нефтесервис», г. Уфа, тел. 89191478793, е-mail: nurievrustam@rambler. ru
Information about the authors
Davletbaev A. Ya., Candidate of Physics and Mathematics, Associate Professor at the Department of Applied Physics, Bashkir State University, Ufa, phone: 89174812699, е-mail: davletbaevay@r ambler. ru
Asalkhuzina G. F., Postgraduate at the Department of Applied Physics, Bashkir State University, Ufa, phone: 89178053644, e-mail: kafedrapf@rambler.ru
Nuriev R. I., Head of the Department, LLC «Production and Commercial Company «GIS-Oilservice», Ufa, phone: 89191478793, е-mail: nurievrustam@rambler.ru