Prediction of Intraformational Remaining Oil Distribution Based on Reservoir Heterogeneity: Application to the J-Field

Zhang, Jie; Fang, Feifei; Wang, Jie; Tian, Yajie; Mo, Fei; Li, Qi; Li, Sainan; Huang, Xiaoliang; Yang, Yi

doi:https://doi.org/10.1155/2021/8870274

Advances in Civil Engineering

On this page

Abstract Introduction Conclusions Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2021 | Article ID 8870274 | https://doi.org/10.1155/2021/8870274

Prediction of Intraformational Remaining Oil Distribution Based on Reservoir Heterogeneity: Application to the J-Field

Jie Zhang,¹Feifei Fang,²Jie Wang,²Yajie Tian,³Fei Mo,²Qi Li,¹Sainan Li,²Xiaoliang Huang,²and Yi Yang¹

Academic Editor: Xuemei Liu

Received30 Jun 2020

Revised28 Jan 2021

Accepted12 Mar 2021

Published31 Mar 2021

Abstract

At high water cut stage, the study of remaining oil distribution in water-flooding reservoir is the basis of implementing potential-tapping measures and enhancing oil recovery. At present, most of the oilfields in China have entered the stage of ultrahigh water cut. The reserves of the oilfields are highly developed, the situation of water flooding is extremely complex, and it is difficult to predict the distribution of the remaining oil, which seriously restricts the adjustment of the production measures, tapping the potential and improving the ultimate recovery rate. In view of aforementioned difficulties, this study puts forward a research approach to predict remaining oil distribution based on reservoir heterogeneity, which can quantitatively characterize reservoir heterogeneity. In order to avoid the drawback that a single parameter cannot fully describe the characteristics of pore structure, the composite index of pore structure (SQRT(K/Φ)) is introduced to study the pore microstructure. The composite index of pore structure is used to predict the distribution of remaining oil in the formation, and the results are basically consistent with those calculated by numerical simulation. It is concluded that the larger the fractal dimension of the composite index of pore structure is, the stronger the heterogeneity of reservoir is; the smaller the composite index of pore structure is, the smaller the recovery degree is. The composite index of pore structure is used to analyze and predict the distribution of remaining oil in the layer, which provides a new direction for the prediction method of remaining oil.

1. Introduction

At present, the development of the main oilfields in China mostly enters the stage of high water cut stage, the distribution of underground fluid is becoming more and more complex, and the exploitation is becoming increasingly challenging. Therefore, it is of great significance to know exactly the remaining oil reserves and their distribution for the adjustment of the production measures, potential tapping, and ultimate recovery enhancement of the oilfield, and it is also an important means to realize the development strategy of “water control and oil stabilization” [1].

In recent years, although the method of remaining oil description has developed rapidly, most of them are qualitative description except numerical simulation method, and the further development of this technology is hindered by the degree of popularization, the depth of research, and the high cost of testing [2]. Present research on remaining oil mainly focuses on the following three aspects: (1) remaining oil distribution and characteristic study [3], (2) measuring and monitoring study on remaining oil saturation (ROS) [4], and (3) potential-tapping technology study [5].

Because of the complexity of oil-water relationship of the reservoir in the high water cut stage, it is more difficult to study the distribution of remaining oil. It is inevitable to have great limitations to predict the distribution of remaining oil based on a single subject. Therefore, only the application of multidisciplinary theory, method, and technology can accurately predict the distribution of remaining oil.

The manifested evidences suggest that interlayer heterogeneity, interior-layer heterogeneity, and plane heterogeneity have an effect on the remaining oil distribution [6]. Because of the difference of heterogeneity in the reservoir layer, the water drive effect of different microfacies in the vertical direction is different, and the distribution of remaining oil is uneven and concentrated in the strongly heterogeneous region. Influenced by the heterogeneity of reservoir layers, the remaining oil is usually concentrated in the single layer with large permeability range and poor physical properties. In the plane, restricted by the distribution, scale, connectivity, and spatial combination of sand bodies, the microfacies with relatively poor physical properties and large permeability range in the reservoir are also the main enrichment areas and potential exploration targets of remaining oil.

Webber puts forward the classification system of quantitative recognition of heterogeneity in the stage of oilfield evaluation and development [7]. Considering heterogeneity and its influence on fluid flow, reservoir heterogeneity is divided into seven types: closed or semiclosed or unclosed fault, unit boundary, permeable layer in unit, interlayer in unit, lamination or cross bedding, micro-heterogeneity, and closed or open fracture. Daws and Prosser proposed to divide reservoir heterogeneity into five grades from large to small: sedimentary strata, sedimentary facies, boundary/bedding, lamination, and diagenesis [8]. According to the characteristics of continental reservoirs and production practice, Qiu et al. put forward a classification scheme to study reservoir heterogeneity from interlayers and interior layers, plane, and microlevel [9]. In addition, there are some effective methods for studying reservoir heterogeneity that uses dynamic-based indicators such as tracers and local grid refinement and downscaling [10–13].

The heterogeneity of oil and gas reservoir is not only the main factor to control oil and gas accumulation but also the basis of remaining oil distribution research [14, 15]. With the increasingly severe exploitation situation of major oil fields and the deepening understanding of the importance of remaining oil exploration, the research on reservoir heterogeneity has been deepened and developed, and the research contents and methods have been improved and enriched.

According to the present situation of oil and gas field development in China, the relationship between reservoir heterogeneity and remaining oil distribution is discussed by studying reservoir heterogeneity. We introduce the parameter of the composite index of pore structure and apply it to the study of the heterogeneity in the macroscopic layer. The composite index of pore structure can accurately represent not only the characteristics of micropore structure but also the degree of heterogeneity in the reservoir, which is a kind of parameter connecting the micro- and macrodescription. Based on this, the composite index of the pore structure combined with the numerical simulation method is proposed to analyze and predict the distribution of remaining oil.

2. Quantitative Characterization of Reservoir Heterogeneity

According to the reservoir classification scheme put forward by Qiu and Chen [9], reservoir heterogeneity is classified into interlayer heterogeneity, plane heterogeneity, interior-layer heterogeneity, and pore heterogeneity. Interlayer heterogeneity includes the rhythmicity of rock granularity, the sequence of bedding structure, the permeability difference and the location of high permeability section, the distribution frequency and size of discontinuous thin shale in the layer, the distribution of other impermeable interlayer, and the ratio of vertical permeability. Pore heterogeneity mainly refers to the heterogeneity of micropore structure, including sand body pore, throat size and its uniformity, configuration, and connectivity of pore throat. In this paper, the quantitative characterization of heterogeneity of the study area is described from a macroperspective.

2.1. Geological Characteristics

The study area is about 25 km east of “K” City, which is located in a reverse fault zone in the northwest margin of “J” basin. The north, southeast, and west are, respectively, bounded by faults. The faults in the north act as the function of blocking the reservoir. The study area is a structural-lithological reservoir blocked by faults (Figure 1).

In general, the study area is a NW-SE structural reservoir blocked by faults, with a dip angle of 2–30, and the depth of the middle part of the reservoir is 1340 m. The block is formed by three faults.

Sedimentary microfacies can be divided into three subfacies: fan root, fan subfacies, and fan margin. Among them, fan root subfacies mainly develop braided channel microfacies, fan subfacies mainly develop braided channel and slow flow sand microfacies, and fan margin mainly develops overland mud microfacies.

The “K” formation in the study area is unconformity deposited on the carboniferous reservoir, and the upper layer is successively deposited by the Baijiantan formation (T₃b), Badaowan Formation in Jurassic System (J₁b), Sangonghe Formation (J₁s), Xishanyao Formation (J₂x), Toutunhe Formation (J₂t), Qigu formation (J₃q), and Tugulu group (K1t) of Cretaceous System. The study area is K₇ sand formation, and K₇ sand consists of four layers K₇¹, K₇², K₇³, and K₇⁴and is subdivided into nine layers K₇¹, K₇²⁻¹, K₇²⁻², K₇²⁻³, K₇³⁻¹, K₇³⁻², K₇³⁻³, K₇⁴⁻¹, and K₇⁴⁻². The main oil layers are K₇², K₇³, and K₇⁴.

2.2. Study on the Heterogeneity of Permeability in Layers

The heterogeneity in a single sand layer refers to the change of its internal vertical reservoir properties, including the change of permeability in the vertical layer. The variation range of porosity of K₇ sand group in the study area is mainly from 12%–16%, and the variation range of permeability is 8 × 10⁻³ μm²–17 × 10⁻³ μm². According to the results of stratigraphic division, K₇ sand consists of four layers K₇¹, K₇², K₇³, and K₇⁴ and is subdivided into nine layers K₇¹, K₇²⁻¹, K₇²⁻², K₇²⁻³, K₇³⁻¹, K₇³⁻², K₇³⁻³, K₇⁴⁻¹, and K₇⁴⁻².

The permeability heterogeneity parameters of each sublayer are shown in Table 1.

The analysis of permeability heterogeneity shows that porosity and permeability have little change, and the heterogeneity in the layer is weak. In general, the average permeability difference of each small layer in K₇ sand formation is small, and the heterogeneity is not strong [2].

Generally speaking, the larger the V_k, T_k, and J_k, the stronger the reservoir heterogeneity, otherwise weaker. The classification standard of the heterogeneity of clastic rock reservoir in China mainly adopts the standard issued by China National Petroleum Corporation (CNPC), as shown in Table 2.

According to Table 2, we judged that the heterogeneity type in K₇ sand group is a relatively homogeneous type. Permeability contour map of K₇ sand formation is shown in Figure 2. It can be seen that the plane permeability is also relatively uniform, mainly concentrated in 10∼18 mD.

2.3. Study on Heterogeneity in Layers with Pore Structure Index

The size and distribution of pore radius can be roughly reflected by the pressure and saturation of mercury withdrawal [16, 17]. It is not only a function of pore throat radius distribution and pore volume but also a function of porosity, permeability, and saturation. The quantitative characteristic parameters of capillary pressure curve, which can characterize the pore structure of reservoir, mainly include displacement pressure, saturation median pressure and median radius, maximum mercury saturation, and efficiency of mercury withdrawal, which are the important basis to reflect the microheterogeneity of the reservoir. Mercury injection data of K₇¹, K₇², K₇³, and K₇⁴ are collected, as shown in Figure 3.

(a)

(b)

(c)

(d)

According to the mercury injection data of each layer, the curves of the K₇⁴ layer are similar, which proves that the pore structure of the rock samples taken from the K₇⁴ layer is similar. After analysis and comparison, we conclude that the heterogeneity of four small layers is gradually weakened.

The composite index of the pore structure in four layers of K₇ sand group is used to study the reservoir heterogeneity, which can truly reflect the pore structure characteristics of samples. Using geological model data to calculate the frequency distribution of composite index of pore structure of each layer, it can be seen from Figure 4 that the heterogeneity in layer K₇¹ is strong, and the index is widely changed between 0.5∼2.5 as can be seen from Figure 4(a). The heterogeneity in layer K₇⁴ is relatively weak, and the index is mainly changed between 0.5∼1.5 as can be seen from Figure 4(d).

(a)

(b)

(c)

(d)

It is generally believed that the fractal dimension of SQRT (K/Φ) in three-dimensional Euclidean space is mostly in the range of 2–3 [18, 19]. The closer the fractal dimension is to 2, the weaker the heterogeneity is and the better the reservoir capability is. The closer the fractal dimension is to 3, the rougher the pore throat surface is, the worse the reservoir capability is, and the stronger the heterogeneity of the reservoir is [20]. We used the core data to analyze the fractal dimension of the composite index of pore structure SQRT (K/Φ) in the layer, as shown in Figure 5. The fractal dimension of the K₇¹ layer is 2.88, that of the K₇² layer is 2.74, that of the K₇³ layer is 2.73, and that of the K₇⁴ layer is 2.38. The fractal dimension of the K₇¹ layer is the largest, so the heterogeneity in the layer is relatively the strongest. At the same time, the fractal dimension of the K₇⁴ layer is the smallest, and the heterogeneity in the layer is also the smallest. We conclude that, from the K₇¹ layer to K₇⁴ layer, the heterogeneity in the reservoir is gradually weakened.

(a)

(b)

(c)

(d)

The intraformational heterogeneity obtained by using the frequency distribution of the composite index of the pore structure calculated by the geological model data is consistent with the intraformational heterogeneity obtained by using the core data to analyze the interlayer’s composite index of pore structure. The larger the fractal dimension of SQRT (K/Φ) is, the stronger the heterogeneity is. The composite index of the pore structure can well reflect the heterogeneity in the layer.

3. Study on the Distribution Law of Remaining Oil in the Layer Based on Heterogeneity

The macroheterogeneity of the reservoir is mainly reflected in the geometry of sand body, reservoir lithology, physical property, permeability heterogeneity, and interlayer distribution, which have different controls over the distribution of remaining oil. The main factors controlling the distribution of remaining oil in this study area are as follows.

The heterogeneity is the key factor that affects the sweep volume of vertical injection water and the distribution of remaining oil. In the positive rhythm sand body, the remaining oil is mainly distributed in the upper-middle section with low permeability and weak washing degree. In the reverse rhythm sand body, the swept volume of injected water is large, the displacement of injected water is uniform, and the upper and lower parts are all the potential areas of remaining oil. It is generally believed that water is always injected along the relatively high permeability section, resulting in high water drive degree and low remaining oil saturation.

The intraformational heterogeneity obtained by using the frequency distribution of the composite index of the pore structure calculated by the geological model data is consistent with the intraformational heterogeneity obtained by using the core data to analyze the interlayer’s composite index of the pore structure. From the K₇¹ layer to K₇⁴ layer, the heterogeneity in the reservoir is gradually weakened. From the oil saturation of the K₇² and K₇⁴ layer (Figure 6), the displacement effect of the K₇⁴ layer is better than that of the K₇² layer.

(a)

(b)

It can be seen from Table 3 that the recovery degree of the K₇⁴ layer is 19.18% and that of the K₇² layer is 27.20%.

The frequency distribution of the composite index of pore structure of K₇² and K₇⁴ are shown in Figure 7.

(a)

(b)

The displacement effect of the K₇⁴ layer is better than that of the K₇² layer, but the recovery degree of the K₇⁴ layer is lower than that of the K₇² layer. This is because, although the heterogeneity of the K₇⁴ layer is weaker than that of the K₇² layer, the composite index of the pore structure of the K₇⁴ layer is mainly concentrated in 0.5–1.5, while the composite index of the pore structure of the K₇² layer is widely distributed in the range from 0.5–3. The composite index of the pore structure of the K₇⁴ layer is relatively smaller than the K₇² layer. As we know, the smaller the composite index of pore structure is, the lower the recovery degree is, so the recovery degree of K₇⁴ layer is relatively smaller.

It can be seen that, in the process of oil displacement, the heterogeneity in the layer affects the distribution of remaining oil, and the remaining oil in the layer with small composite index of the pore structure will be relatively enriched. The composite index of the pore structure can be used to predict the distribution of remaining oil.

4. Conclusions

In this study, the composite index of the pore structure (SQRT(K/Φ)) is introduced to study the pore microstructure of the reservoir, which can be used to quantitatively characterize the heterogeneity of the reservoir, and this parameter (SQRT(K/Φ)) can comprehensively describe the characteristics of the micropore structure well. The interlayer distribution of remaining oil can be predicted by using the composite index of the pore structure (SQRT(K/Φ)), and the result is basically consistent with that calculated by the numerical simulation method. It is found that the larger the fractal dimension of the composite index of the pore structure is, the stronger the heterogeneity of reservoir is; the smaller the composite index of the pore structure is, the smaller the recovery degree is. The composite index of the pore structure can be used to analyze and predict the interlayer distribution of remaining oil, which provides a new direction for the prediction method of remaining oil, but we have to also acknowledge that our proposed composite index appears to work for the particular field, and it is not universal and may not for other formations.

In this study, K₇ sand formation of an oilfield in China is taken as an example, and the composite index of the pore structure (SQRT(K/Φ)) is used to analyze the intraformational heterogeneity and predict the distribution of remaining oil in the formation. The conclusions are(1)The intraformational heterogeneity is related to the composite index of the pore structure. The larger the fractal dimension of the composite index of the pore structure is, the stronger the intraformational heterogeneity of the reservoir is; the smaller the fractal dimension of the composite index of the pore structure is, the weaker the intraformational heterogeneity of the reservoir is.(2)The interlayer distribution of remaining oil can be predicted by using the composite index of pore structure (SQRT(K/Φ)), and the result is basically consistent with that calculated by the numerical simulation method. The smaller the index is, the lower the recovery degree is, and the remaining oil is also relatively richer. This also explains the abnormal phenomenon that the displacement effect of the K₇⁴ layer is better, but the recovery degree is smaller.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

All authors declare that there are no conflicts of interest regarding the publication of this article.

Acknowledgments

The work was financially supported by the National Science and Technology Major Project of China (Grant no. 2017ZX05013-001).

References

C. Li, X. L. Yang, and G. L. Zhadan, “Technical countermeasures for water control and oil stabilization in high water-containing exploitation period of grit frame reservoir,” Journal of Xinjiang Petroleum Institute, vol. 16, no. 3, pp. 56–59, 2004.
View at: Google Scholar
H. Q. Tan, Q. Fu, and L. X. Li, “Application of streamline numerical simulation method to study of remaining oil distribution in high water-cut stage oilfield,” Journal of Chengdu University of Technology (Science & Technology Edition), vol. 44, no. 1, pp. 30–35, 2017.
View at: Google Scholar
S. B. Zhang, G. M. Yao, and Z. Q. Dong, “Study on the characteristics of remaining oil saturation distribution,” Fault-Block Oil & Gas Field, vol. 7, no. 3, pp. 19–22, 2000.
View at: Google Scholar
A. A. Al-harbi, D. P. Schmitt, and S. X. Mark-Ma, “Toward quantitative remaining oil saturation (ROS): determination challenges and techniques,” in Proceedings of theSPE Annual Technical Conference and Exhibition, Denver, Colorado, USA, October 2011.
View at: Google Scholar
Z.-z. Han, R.-c. Yang, A.-p. Fan, Q.-c. Chen, and Y.-t. Shao, “Remaining oil distribution in Ng33 bottom water reservoir of Lin 2-6 fault-block in Huimin depression and potential tapping in horizontal well,” Mining Science and Technology (China), vol. 19, no. 1, pp. 102–107, 2009.
View at: Publisher Site | Google Scholar
J. Li, Y. Liu, Y. Gao, B. Cheng, F. Meng, and H. Xu, “Effects of microscopic pore structure heterogeneity on the distribution and morphology of remaining oil,” Petroleum Exploration and Development, vol. 45, no. 6, pp. 1112–1122, 2018.
View at: Publisher Site | Google Scholar
K. J. Webber, “How heterogeneity affects oil recovery,” in Reservoir Characterization, L. W. Lake and H. B. Carroll, Eds., pp. 487–544, Academic Press, Orlando, FL, USA, 1986.
View at: Google Scholar
J. A. Daws and D. J. Prosser, “Scales of permeability heterogeneity within the Brent group,” Journal of Petroleum Geology, vol. 15, no. 4, pp. 397–417, 1992.
View at: Publisher Site | Google Scholar
Y. N. Qiu and Z. Q. Chen, Reservoir Description, Petroleum Industry Press, Beijing, China, 1996.
G. M. Shook and K. M. Mitchell, A Robust Measure of Heterogeneity for Ranking Earth Models: The F PHI Curve and Dynamic Lorenz Coefficient, Society of Petroleum Engineers, Muscat, Oman, 2009.
B. Rashid, A. H. H. Muggeridge, A. Bal, and G. Williams, “Quantifying the impact of permeability heterogeneity on secondary-recovery performance,” SPE Journal, vol. 17, no. 02, pp. 455–468, 2012.
View at: Publisher Site | Google Scholar
V. A. Torrealba, H. Hoteit, and A. Chawathe, “Improving chemical EOR simulations and reducing the subsurface uncertainty using downscaling conditioned to tracer data,” in Proceedings of the SPE Annual Technical Conference and Exhibition, Society of Petroleum Engineers, San Antonio, Tx, USA, October 2017.
View at: Publisher Site | Google Scholar
H. Hoteit and A. Chawathé, Making Field-Scale Chemical EOR Simulations a Practical Reality Using Dynamic Gridding, Society of Petroleum Engineers, Muscat, Oman, 2014.
H. N. Li, H. M. Xu, and N. Xu, “The heterogeneity model of low permeability reservoir and remaining oil distribution-a case study of Dujiatai oil zone in Qi 9-Huan 50 block, the west liaohe depression,” Petroleum Geology & Experiment, vol. 28, no. 4, pp. 404–408, 2006.
View at: Google Scholar
W. Lin, S. C. Xiong, and Y. Liu, “Spontaneous imbibition in tight porous media with different wettability: pore-scale simulation,” Physics of Fluids, vol. 32, no. 11, 2021.
View at: Publisher Site | Google Scholar
X. L. Zhao, Z. M. Yang, and W. Lin, “Study on pore structures of tight sandstone reservoirs based on nitrogen adsorption, high-pressure mercury intrusion, and rate-controlled mercury intrusion,” Journal of Energy Resources Technology, vol. 141, no. 11, Article ID 112903, 2019.
View at: Publisher Site | Google Scholar
X. L. Zhao, Z. M. Yang, and Z. Y. Wang, “An experimental study on stress sensitivity of tight sandstones with different microfractures,” Advances in Civil Engineering, vol. 2020, Article ID 1865464, 10 pages, 2020.
View at: Publisher Site | Google Scholar
X. Zhao, Z. Yang, W. Lin et al., “Fractal study on pore structure of tight sandstone based on full-scale map,” International Journal of Oil, Gas and Coal Technology, vol. 22, no. 2, pp. 123–139, 2019.
View at: Publisher Site | Google Scholar
W. Lin, X. Li, Z. Yang et al., “A new improved threshold segmentation method for scanning images of reservoir rocks considering pore fractal characteristics,” Fractals, vol. 26, no. 02, Article ID 1840003, 2018.
View at: Publisher Site | Google Scholar
W. Lin, Z. Yang, X. Li et al., “A method to select representative rock samples for digital core modeling,” Fractals, vol. 25, no. 04, Article ID 1740013, 2017.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Jie Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies