CC BY
29
ê Dmitry A. Martyushev
Improving the geological and hydrodynamic model of a carbonate oil object.
UDC 622.276.6
Improving the geological and hydrodynamic model of a carbonate oil object by taking into account the permeability anisotropy parameter
Dmitry A. MARTYUSHEV
Perm National Research Polytechnic University, Perm, Russia
Significant share of the developed oil assets related to carbonate complex-built objects has formidably increased in Russia, including the Perm Region. Reliable knowledge of the parameters for the cavern-pore type of the reservoir allows clarifying the existing geological and hydrodynamic models (GHM), selecting a rational development system, regulating the development processes and providing optimal geological and technical measures for this formation. In the construction and adaptation of GHM for oil fields, especially those related to complex-built carbonate reservoirs, knowledge of both horizontal and vertical permeability (anisotropy parameter) is important. When creating GHM of carbonate objects in Perm Region deposits, vertical permeability is often taken to be zero, although this is far from being the case. Determining the vertical permeability (anisotropy parameter), its dynamics when changing the formation and bottomhole pressures and using it in GHM is an urgent task that will improve the quality and reliability of using digital models to calculate and predict the oil production process.
Article describes the methodology for determining permeability anisotropy according to the interpretation of hydrodynamic investigations of wells. Proposed methodology for determining the anisotropy parameter processed the results of more than 200 studies conducted on production and injection wells of the Famennian deposit at the Gagarinskoye field. For each lithological-facies zone, dependence of the permeability anisotropy index on the bot-tomhole pressure is constructed. To predict and evaluate the effectiveness of the applied geological and technical measures and technological development indicators, author modified the geological and hydrodynamic model taking into account the obtained dependencies on the change in the anisotropy parameter.
Using a modified hydrodynamic model, it was possible to significantly improve the adaptation of both production and injection wells. Thus, the quality and reliability of the digital model of the Famennian deposit at the Gagarin-skoye field for calculating and predicting the oil production process has improved.
Key words: vertical permeability; horizontal permeability; permeability anisotropy parameter; complex-built carbonate deposit; modified geological and hydrodynamic model
How to cite this article: Martyushev D.A. Improving the geological and hydrodynamic model of a carbonate oil object by taking into account the permeability anisotropy parameter. Journal of Mining Institute. 2020. Vol. 243, p. 313-318. DOI: 10.31897/PMI.2020.3.313
Introduction. In the Perm Region, the share of developed oil assets related to complex carbonate objects has significantly increased. Complex carbonate reservoirs are characterized by different types of cavernosity (pores, cracks, cavities), significant inhomogeneity, fluid flows between types of cavernosity, as well as the difference between lateral and vertical permeability, and these factors significantly affect the productivity of production wells and an overall oil recovery factor [1, 11, 14, 16]. Reliable knowledge on the parameters of the cavern-pore type of the reservoir allows clarifying the existing geological and hydrodynamic models (GHM), select a rational development system, adjust development processes and provide optimal geological and technical measures (GTM) for the given formation [5, 10, 13].
In the construction and adaptation of oil fields' GHM, especially those related to complex carbonate reservoirs, knowledge of both horizontal and vertical permeability (anisotropy parameter) is important [10, 14]. When creating GHM of carbonate objects in the Perm Region's deposits, vertical permeability is often taken to be zero, although this is far from being the case (it is known that cav-ernosity can exist in the mudstone formation, although permeability is almost zero in laboratory study) [10]. Knowledge of vertical permeability under conditions of a complex carbonate reservoir is of particular importance with the active role of formation water. Oil production may be high for some period, but then it will decrease when a bottom water cone appears, which is often formed during irrational development of reserves and unreasonable technological regimes of production wells' operation. Process of cones' formation is much more intense in cavern reservoirs with vertical
ê Dmitry A. Martyushev
Improving the geological and hydrodynamic model of a carbonate oil object.
cracks, since they extend far down into the underlying formations, forming paths for creating water cones [6, 15, 17]. Thus, determination of vertical permeability (anisotropy parameter) and its use in GHM is a relevant task that will improve the quality and reliability of digital models' application for calculating and predicting the oil production process [3, 4, 8, 18].
Methods for determining the anisotropy parameter (vertical permeability). Permeability anisotropy parameter in carbonate reservoirs can be determined by geophysical, field and hydrodynamic investigations (HDI) of wells and formations [7]. Examining of rock sample and geophysical study cannot provide accurate knowledge on the features of the formation in connection with the investigation peculiarities. The most developed methods for determining the parameter of permeability anisotropy are the methods of well interference testing and tracing of marked substances, however, they are quite time consuming and expensive. Method of tracing indicators can be reliably implemented at a later stage of field development with intensive fluid watering. Particular attention should be paid to the interpretation of data from hydrodynamic investigations of wells (PRC/LRC (level recovery curve)/PSC (pressure stabilization curve) as the most common method of monitoring field development, which allows solving a huge range of production and scientific problems [7, 12].
Determination of permeability anisotropy parameter according to hydrodynamic investigations of wells. In [9] authors described a method for determining vertical permeability and permeability anisotropy parameter. Essence of the technique is as follows: initial pressure recovery curve (PRC) is reconstructed in the coordinates Pbot(t) - lg(t). Final part of the PRC is considered -final section is distinguished and the slope P is determined. PRC is also built in coordinates Pbot(t) - 1/t05, on which the straight section is selected so that the time corresponding to its end is less than the time corresponding to the beginning of the straight section in coordinates Pbot(t) - lg(t), and the slope of the section o is determined. And further, if drilled (working) thickness hdr is known, values o and P and the coordinates of the last point of the line constructed in the coordinates Pbot(t) - 1/t0 5, following properties are sequentially determined: total formation thickness, vertical permeability Pv, vertical piezoconductivity, horizontal permeability Pg, and consequently the anisotropy parameter.
Relevance and value of these calculations is obvious, since, on the one hand, informational content of the well research by pressure recovery is enhanced, and on the other hand, it is possible to obtain very valuable information about the anisotropy of the formation without the cost of special investigations.
Using this technique, results of more than 200 investigations conducted on the production and injection wells in the Famennian deposit at the Gagarinskoye field were processed. Carbonate oil deposit in the Famennian deposits at the Gagarinskoye field is characterized by double cavernosity, and the geological model for the deposit formation with a successive change in the following lithological-facies sedimentation conditions can be considered typical: reef slope, lower and upper rear drags, bioherm core [7]. Deposits of the upper rear drags, related to the central interreef part of the formation are characterized by the largest capacitive space. Facies of the lower rear drag, bio-herm core, and reef slope are confined to lower relief sections of the geological section. During sedimentation, a greater amount of micritic material was carried to them, which reduced their capacitive characteristics. On a number of deposits, rifogenic formations, which are usually characterized by low permeability and porosity properties, are also isolated into a separate lithological-facies zone. Complex structure of the reservoir, caused by sedimentation conditions, led to alternation in the section and on the reservoirs area of different capacities for each stratigraphic range [2, 7]. For each lithological-facies zone, dependence of the permeability anisotropy index on the bottomhole pressure was constructed (Fig. 1). Obtained dependencies can be used to predict the permeability anisotropy parameter in each lithological-facies zone, which may affect the selection of the optimal technological mode of the well operation and the choice of geological and technical measures.
314 -
Journal of Mining Institute. 2020. Vol. 243. P. 313-318 • Oil and Gas
ê Dmitry A. Martyushev
Improving the geological and hydrodynamic model of a carbonate oil object.
To predict and evaluate the effectiveness of the applied geological and technical measures and technological development indicators, author of the article modified the GHM taking into account the obtained dependencies for the change in the permeability anisotropy parameter. Modified GHM, compared with the base model, which does not take into account the anisotropy index, will allow reliable evaluation of the technological development indicators and well operation modes after the completion of GTM.
Improving the geological and hydrodynamic model taking into account changes in the anisotropy index. There are many ways to account for cavernosity in hydrodynamic models. The most laborious, but at the same time physical, is the creation of a double medium model. Obvious advantage of using this GHM is taking into account the complex structure of the cavities during multiphase filtration of the formation fluids. However, the process of building the model and its use is associated with significant difficulties, such as the ambiguity of determining the parameters of the cavern medium and the model's overloading with a slowdown in the calculation speed as a result of the permeability and porosity properties' distribution between two cavities, as well as the interflows between them. Therefore, such models are rarely used. Another method of GHM is based on the use of seams with low porosity and high permeability, such as the area of wells. Non-adjacent junctions are also used, in which the fluid can be transferred almost instantly from one point in the formation to another. Indirect modeling of the filtration in the cavern-pore reservoir can be done using relative phase permeability curves (RPP). There are methods to account for changes in reservoir permeability depending on reservoir pressure dynamics. A similar method was chosen to take into account the cavern component of the cavity space of the Famennian deposit at the Gagarinskoye field. It gave satisfactory indicators with rational time costs for modifications. Modeling was performed in a hydrodynamic simulator "Tempest" ver. 8.3.1 by Roxar.
Preliminary vertical permeability was numerically adjusted according to the results of the interpretation of hydrodynamic investigations (HDI). Then, the process of modifying the hydrodynamic model took place by searching for permeability multipliers in several stages:
1. Converting the obtained dependence of the anisotropy index on the ratio of the current bot-tomhole pressure to the initial reservoir pressure into the dependence of the anisotropy index on the current reservoir pressure. Dependencies will be the same.
2. Dependence of permeability on reservoir pressure is constructed as a result of the hydrody-namic investigations' analysis, conducted under transient conditions. Permeability obtained in the investigations will be equal to the lateral permeability Kx = Ky.
3. With a sufficiently small pressure step (for example, 5 bar) from atmospheric to reservoir pressure (discreteness can be reduced above this threshold), permeability values along the lateral are obtained.
4. Value of horizontal permeability is substituted into the dependence of the anisotropy index on the current reservoir pressure. Result is the distribution of vertical permeability with regard to reservoir pressure.
It is assumed that at the initial formation pressure the permeabilities are not changed, and when the pressure changes in one direction or another, they increase or decrease. Since the distribution of
1 2 3
Fig. 1. Dependence of anisotropy index on bottomhole pressure for lithological-facies zones of the Famennian deposit at Gagarinskoye field 1 - facies of bioherm core, y = 34.719e~1037x, R2 = 0.5255; 2 - facies of the upper rear drag, y = 0.0758x~2 32", R2 = 0.6562; 3 - facies of the lower rear drag, y = 28.606e~" llx, R2 = 0.6
Dmitry A. Martyushev
Improving the geological and hydrodynamic model of a carbonate oil object.
S I
w s
s-i 3
s m
in
O 3
P S
140 120 100 80 60 40 20
100 200 300 Formation pressure, bar
400
2000 1800 1600 1400 1200 1000 800 600 400 200 0
3
permeability in the area of each well can vary significantly, it is not the permeability values that are needed for use, but their multipliers (with Pfr = Pfrin the multiplier is equal to one). Process of modifying the hydrodynamic model by searching for permeability factors is shown in Fig.2.
Analyzing the obtained graphs, it can be concluded that with a significant decrease in formation pressure, the role of vertical filtration becomes significant. This procedure was performed for each litho-logical-facies zone of the development object under consideration. Thus, there are dependencies of lateral and vertical permeability changes for each zone. In contrast to the common method using the KVSP keyword, the use of the KVPX, KVPY, and KVPZ keywords allows for multidirectional distribution of changes in permeability with regard to pressure, which was done. Appearance of the modified geological and hydrodynamic model is presented on the example of a cube of altered regions (Fig.3).
When working with the modified model, it was possible to improve adaptation both for production (Fig.4) and injection wells (Fig.5).
Thus, the calculations using the modified hydrodynamic model in well 429 allowed recreating the dynamics of the bottomhole pressure, adjusting the formation pressure trend and significantly improve the convergence of the water cut parameter. In well 71, the necessary historical repression was achieved, and the entire volume of injected liquid flowed into the formation. This did not work on the original model and the bottomhole pressure was on the shelf for some periods of time (the auto-hydraulic fracturing effect was simulated).
Fig.2. Dependence of permeability multipliers on formation pressure for facies of the upper rear drag
1 - ratio of multipliers Kz/Kx; 2 - Kz; 3 - Kx
Fig. 3. Appearance of the modified geological and hydrodynamic model of the Famennian deposit
at the Gagarinskoye field
Dmitry A. Martyushev
Improving the geological and hydrodynamic model of a carbonate oil object...
a
b
Fluid flow rate
Fluid flow rate of the original model Real fluid flow rate Oil flow rate
Oil flow rate of the original model Real oil flow rate
in in OO CO OO 0 in CO CO CO
o <N CD <N o <N o <N o <N CD <N CD <N CD <N CD <N CD <N O CD <N CD <N CD <N CD <N CD <N CD <N CD <N CD <N CD <N CD <N CD <N
< Q £ < Q £ < Q t < Q < Q t < Q t ¡Ï < Q t < Q
_ Total fluid production
— Total fluid production of original model
— Water cut
■ - Water cut of original model
■ ■ Total oil production
■ - Total oil production of original model
— Real water cut
100
80 60 40 20 0
1E-7 Êf 8E-8 3 6E-8 § 4E-8 ^ 2E-8 0E0
< < Q < < P Injection
Injection of original model Real injection Total injection
Total injection of original model Real total injection
* <
150 100 50 0
in in 22 22
CD <N CD <N CD <N CD <N CD <N CD <N CD <N CD <N CD <N CD <N o
£ (Ü fü & y (Ü ë
U < — Pfr < Q < < < < u
Pfr of original model
Real Pbot
Pbot
■ Pbot of original model Real Pfr
Fig.4. Comparison of adaptation results for production well 429 on a modified and original model: a - flow rates; b - total production and water cut; c - injection; d- pressure
b
Tf \Q OO
Q\
Q\ Q\ Q\
Tf KO OO
C^ C^ o^ o^
200
15 0
S sa
100 o
50
0
500 400 300 200 100 0
^r OO
O, ^ C^ C^
Fig.5. Comparison of adaptation results for injection well 71 on a modified and original model: a - flow rates; b - total production and water cut; c - injection; d- pressure See the legend in Fig.4.
a
c
a
d
c
Conclusion. This article considers an urgent problem aimed at determining the permeability anisotropy parameter (vertical permeability) and taking it into account in existing geological and hydrodynamic models (using the Famennian deposit at the Gagarinskoye field as an example). For each lithological-facies zone of the Famennian deposit at the Gagarinskoye field, the dependence of the change in the permeability anisotropy parameter on the bottomhole pressure is obtained. It has been established that the values of the anisotropy index in various lithological-facies zones vary
ê Dmitry A. Martyushev
Improving the geological and hydrodynamic model of a carbonate oil object.
over a wide range. Using the obtained dependencies for estimating the permeability anisotropy parameter, the current geological and hydrodynamic model of the Famennian deposit at the Gagarin-skoye field was modified. Using a modified hydrodynamic model, it was possible to significantly improve the adaptation of both production and injection wells. Thus, the quality and reliability of the digital model of the Famennian deposit at the Gagarinskoye field for calculating and predicting the oil production process has improved.
Using the technique described in the article, it is possible to evaluate the change in the permeability anisotropy parameter of other carbonate objects at Russian oil fields and further modify the existing geological and hydrodynamic models.
REFERENCES
1. Abrosimov A.A., Shelyago E.V., Yazynina I.V. Justification of representative data volume of porosity and permeability properties for obtaining statistically reliable petrophysical connections. Zapiski Gornogo instituta. 2018. Vol. 233, p. 487-491. DOI: 10.31897/PMI.2018.5.487
2. Begma D.S., Belkina V.A. Lithological-facies features of the Upper Jurassic field T. structure. Izvestiya Tomskogopolitekh-nicheskogo universiteta. Inzhiniring georesursov. 2017. Vol. 328. N 1, p. 109-122 (in Russian).
3. Bozhenyuk N.N., Belkina V.A., Strekalov A.V. Geological model of the Vikulov deposits, taking into account the analysis of reservoir connectivity and data on horizontal wells. Izvestiya Tomskogo politekhnicheskogo universiteta. Inzhiniring georesursov.
2017. Vol. 329. N 4, p. 30-44 (in Russian).
4. Galkin V.I., Ponomareva I.N., Repina V.A. Study of oil recovery in reservoirs of various types of cavities using multivariate statistical analysis. Vestnik Permskogo natsionalnogo issledovatelskogo politekhnicheskogo universiteta. Geologiya. Neftegazovoe i gornoe delo. 2016. Vol. 15. N 19, p. 145-154. DOI: 10.15593/2224-9923/2016.19.5 (in Russian).
5. Kashnikov Yu.A., Gladyshev S.V., Razyapov R.K., Kontorovich A.A., Krasilnikova N.B. Hydrodynamic modeling of the development priority area at the Yurubcheno-Tokhomskoye field taking into account the geomechanical effect of crack closure. Neftyanoe khozyaistvo. 2011. N 4, p. 104-107 (in Russian).
6. Kudryashova D.A. Use of probable-statistical methods to determine the sources of watering of candidate wells for water isolation operations (on the example of the Visean object at the Perm Territory field). Vestnik Permskogo natsionalnogo issledovatel-skogo politekhnicheskogo universiteta. Geologiya. Neftegazovoe i gornoe delo. 2018. Vol. 17. N 1, p. 26-36 (in Russian).
7. Martyushev D.A., Zaitsev R.A. Influence of petrophysical parameters of reefogenic carbonate reservoirs at oil fields of the Tournais-Famennian deposits of the Upper Prikamye on the productivity of production wells. Izvestiya Tomskogo politekhnicheskogo universiteta. Inzhiniring georesursov. 2019. Vol. 330. N 11, p. 77-85. DOI: 10.18799/24131830/ 2019/11/2350 (in Russian).
8. Sayakhutdinov A.I., Kondrateva N.R., Gallyamova D.Ch., Karachurina E.V., Bulatova E.S. Hydrodynamic modeling technique for complex carbonate cavern reservoirs using the example of the Mancharovskoye field. Neftyanoe khozyaistvo. 2014. N 2, p. 114-115 (in Russian).
9. Cherepanov S.S., Martyushev D.A., Ponomareva I.N., Khizhnyak G.P. Estimation of the permeability anisotropy of carbonate reservoirs by pressure recovery curves. Neftyanoe khozyaistvo. 2013. N 4, p. 60-61 (in Russian).
10. Repina V.A., Galkin V.I., Galkin S.V. Complex petrophysical correction in the adaptation of geological hydrodynamic models (on the example of Visean pool of Gondyrev oil field). Zapiski Gornogo instituta. 2018. Vol. 231, p. 268-274. DOI: 10.25515/PMI.2018.3.268
11. Khramchenkov M.G., Korolev E.A. Dynamics of the cracks development in oil-saturated carbonate formations of the Bashkir layer in the Republic of Tatarstan. Neftyanoe khozyaistvo. 2017. N 4, p. 54-57 (in Russian).
12. Tsagan-Mandzhiev T.N. Improving the reliability of determining the vertical permeability of the formation according to hydrodynamic investigations. Gazovayapromyshlennost. 2012. N 5, p. 19-23 (in Russian).
13. Cherepanov S.S., Chumakov G.N., Galkin S.V. Possibilities of accounting for caverns of reservoirs during geological and hydrodynamic modeling of deposits development with water flooding. Neftepromyslovoe delo. 2016. N 8, p. 5-8 (in Russian).
14. Guerriero V., Mazzoli S., Iannace A., Vitale S., Strauss C. A permeability model for naturally fractured carbonate reservoirs. Marine and petroleum geology. 2013. Vol. 40, p. 115-134. DOI: 10.1016/j.marpetgeo.2012.11.002
15. Litvinenko V. Advancement of geomechanics and geodynamics at the mineral ore mining and underground space development. Geomechanics and Geodynamics of Rock Masses. International European Rock Mechanics Symposium. EUROCK
2018. Saint Petersburg, Russian Federation, 22 May 2018. Taylor and Francis Group, London, UK. 2018. Vol.1, p. 3-16.
16. Menezes F.F. Anisotropy of volume change and permeability evolution of hard sandstones under triaxial stress conditions. Journal of Petroleum Science and Engineering. 2019. Vol. 174, p. 921-939. DOI: 10.1016/j.petrol.2018.11.079
17. Noirot M., Massonnat G., Jourde H. On the use of Wireline Formation testing (WFT) data: 2. Consequences of permeability anisotropy and heterogeneity on the WFT responses inferred flow modeling. Journal of Petroleum Science and Engineering. 2015. Vol. 133. P. 776-784. DOI: 10.1016/j.petrol.2013.08.055
18. Zhang W., Wang Q. Permeability anisotropy and gas slippage of shales from the Sichuan Basin in South China. International Journal of Coal Geology. 2018. Vol. 194. p. 22-32. DOI: 10.1016/j.coal.2018.05.004
Author Dmitry A. Martyushev, Candidate of Engineering Sciences, Associate Professor, martyushevd@inbox.ru (Perm National Research Polytechnic University, Perm, Russia). The paper was received on 17 September, 2019. The paper was accepted for publication on 27 December, 2019.
318 -
Journal of Mining Institute. 2020. Vol. 243. P. 313-318 • Oil and Gas
