ORIGINAL RESEARCH article

Front. Earth Sci., 02 June 2022

Sec. Structural Geology and Tectonics

Volume 10 - 2022 | https://doi.org/10.3389/feart.2022.918683

Fracture Identification of Deep Dolomite Reservoir Based on R/S-FD Analysis: A Case Study of the Cambrian Sinian Reservoirs in the Sandaoqiao Gas Field, Northern Tarim Basin

  • QM

    Qingxiu Meng 1,2,3

  • WD

    Wenlong Ding 1,2,3*

  • XD

    Xindong Diao 4

  • PH

    Pengyuan Han 1,2,3

  • HW

    Huanhuan Wang 1,2,3

  • ZX

    Zikang Xiao 1,2,3

  • 1. School of Energy Resources, China University of Geosciences, Beijing, China

  • 2. Key Laboratory of Marine Reservoir Evolution and Hydrocarbon Enrichment Mechanism, Ministry of Education, China University of Geosciences, Beijing, China

  • 3. Beijing Key Laboratory of Unconventional Natural Gas Geological Evaluation and Development Engineering, China University of Geosciences, Beijing, China

  • 4. Sinopec Northwest Oilfield Company, Urumqi, China

Article metrics

View details

1

Citations

2,2k

Views

574

Downloads

Abstract

The pre-Mesozoic Sandaoqiao gas field, northern Tarim Basin, Western China, hosts a typical deep fractured dolomite reservoir. Taking this area as an example, this study carried out identification research on the fracture development section of the Cambrian Sinian deep (5,700–6200 m) thick (60–80 m) fractured dolomite reservoirs. The core scale identification and numerical analysis method of heavy gauge finite difference (R/S-FD) are used to identify the fracture development section. In the process of applying this method to the study area, by comparing the fracture development identification results of continuous and complete coring sections, three logging curves with high fracture sensitivity (AC/DEN/CAL) are selected. The adjusted R/S-FD analysis method can effectively identify the fracture development section of the fractured dolomite reservoir. Among them, five fracture development sections were identified in well Q1, with an average thickness of 6.8 m. The fracture development section is in good consistency with the reservoir interpretation section of conventional logging. Well, Q101 identified 11 fracture development sections with an average thickness of 2.5 m. The results show that the gas logging section lags behind the fracture development section, mainly corresponding to the lower part of the identified fracture development section. The thickness and distribution of the longitudinal fracture development section of the two wells are obviously different. The research shows that the R/S-FD method can well identify the fracture development section of a fractured dolomite reservoir.

Introduction

Carbonate reservoirs are one of the most important reservoir types in the world (Zhu et al., 2020). Dolomite reservoir is one of the common carbonate reservoir types (Antonellini and Mollema, 2000; Ma et al., 2008; Guo et al., 2020). Compared with limestone, this kind of reservoir facies is more brittle and has weak dissolution (Zheng H. et al., 2007; Zheng R. C. et al., 2007; Davies and Smith, 2011). The main reservoir space of the Cambrian Sinian reservoir in the sandaoqiao gas field is mainly fractured, and a typical fractured dolomite reservoir is developed. Structural fractures directly affect the quality and production of dolomite reservoirs.

Fractures in reservoirs are the key to high and stable production of oil and gas reservoirs, especially for tight sandstone, shale, and carbonate reservoirs (Chen et al., 2021; Lai et al., 2020; Song et al., 2015; Zhou et al., 2020; Liu et al., 2022; Zhou et al., 2022). The development degree and distribution characteristics of fractures directly affect the physical properties and production of oil and gas reservoirs. In recent years, fractures have always been the focus of reservoir research (Luo et al., 2012; Yin and Wu., 2020; Zhao et al., 2021; Wang et al., 2022). Identifying the characteristics of fracture development sections in underground reservoirs vertically is one of the research contents (Prioul and Jocker, 2009; Santos et al., 2015; Fernández-Ibanez et al., 2018; He et al., 2020; Wang et al., 2021; Yin et al., 2020).

There are many ways to identify the fracture development section longitudinally (Hong et al., 2020; Lan et al., 2021; Zheng et al., 2020; Li, 2022). These methods include geological outcrop, core, logging data, and seismic data (Aghli et al., 2016; Bates et al., 1999; Lai et al., 2017; Liu et al., 2019; Yin et al., 2018a; Li H. et al., 2020; Li, 2022). Considering the recognition accuracy, recognition range, representativeness, cost, and other factors, the use of conventional logging curves is more practical and economical (Miranda and Andrade, 1999; Zhao et al., 2011; Ding et al., 2012; Ding et al., 2013; Xue et al., 2014; Yang et al., 2017; Yin et al., 2018b; Afshari et al., 2018). R/S method uses several conventional logging curves sensitive to fracture response to identify fracture development sections (Mandelbrot and Wallis, 1969; Li et al., 2018; Xiao et al., 2019a; Xiao et al., 2019b; Li et al., 2019; Zhao et al., 2019). The finite difference method improves the accuracy and effect of the R/S method. R/S-FD method is used to identify the longitudinal fracture development section of a single well, and its application in shale and tight sandstone reservoirs is relatively mature (Wang et al., 2018; Xiao et al., 2019a; Li et al., 2019). This study applies this method to a deep fractured dolomite reservoir for the first time. The results have also been verified by other reservoir interpretation methods.

Geological Setting

The Tarim Basin is located in western China and is an important developing oil- and gas-producing area. The Sandaoqiao gas field is located in the northern Tarim Basin (Figure 1), in the western part of the Yakela Faulted Uplift (Han et al., 2016; Yang et al., 2018). The Kuqa depression is located to the north, and the Halahatang depression is located to the south. The Sandaoqiao gas field is generally located in the compression-derived structural slope of the far margin uplift of the foreland basin. The Sinian system of the Cambrian system was severely eroded with the formation of a structural compression fold, forming a buried hill drape structure (Li Y. et al., 2020; Wang et al., 2021), with a reservoir characterized by a denuded residual buried hill (Figure 1). The Cambrian Sinian strata are pinched out on the Yakela Faulted-Uplift oriented in the NW direction (Han et al., 2015). The Carboniferous strata overlying the Sandaoqiao gas field are compact and stable, with good caprock conditions, which prevented the surface water from leaching and dissolving Cambrian Sinian strata during the regional long-term denudation period. These conditions have formed a condensate gas reservoir producing area with a Carboniferous cap rock and developed structural fractures.

FIGURE 1

In the study area, the depth of the Cambrian Sinian system is 5,700–6,200 m, which is deep-ultradeep. The thickness of the Cambrian Sinian strata is uniform. From the late Paleozoic to Mesozoic, this stratum experienced multiple superimposed weak-strong thrusting fault events, accompanied by strong weathering and denudation. The Mesozoic Cretaceous strata cover formed a denudation residual buried hill drape structure. The lithology is characterized by the development of deeply fractured dolomite reservoirs (micritic and powdery). According to the measured physical property data, the porosity and permeability of dolomite reservoirs in the study area are relatively low. Structural fracture is the main reservoir space. The properties of the structural fractures are considered the key factors in determining the reservoir space and seepage performance.

Characteristics of Tectonic Stress

In the geological history period, the study area was mainly affected by NW strong compressive stress and pre-Mesozoic strong denudation, forming a buried hill drape structure with strata pinching out in the NW direction.

According to the analysis results of imaging logging and other logging data, the azimuth of the fast shear wave is NE, and the induced fractures thus strike NE. It is comprehensively inferred that the horizontal principal stress direction is NE. The NE trending fault strike (Figure 2) is roughly consistent with the current horizontal principal stress direction.

FIGURE 2

Fracture Characteristics

According to the seismic data, the main faults in the study area are NE-trending, and most of them terminate in the upper Cambrian strata. They are mainly high-angle and oblique parallel reverse faults with lower-order faults.

The fracture development of the core sections of wells Q1 and Q101 in fine crystalline powder dolomite has resulted in mainly high-angle and half-filled fractures (Figure 3), and conjugate fractures are found locally. The total length of the observed core is 50.3m, and the average coring rate is 76%. We analyzed and photographed all coring sections. Dissolution pores are rare.

FIGURE 3

The cast thin sections show that the dolomite dissolution in this section is very weak (Figure 4). We observed and analyzed a total of 16 cast thin sections, 7 ordinary thin sections, and 10 SEM analysis reports. The typical fracture development parts in the study area are selected for sampling. It is clear that the fractured dolomite reservoir is the main feature in this section of the Sandaoqiao gas field. The scanning electron microscopy (SEM) results showed that there were microcracks in some parts of the section. This dolomite is generally subhedral heteromorphic, without obvious dissolution.

FIGURE 4

Methodology

R/S-FD Method

The numerical analysis method used in this study is based on the R/S method. In R/S analysis, R is the range, that is, the difference between the maximum cumulative deviation and the minimum cumulative deviation, representing the complexity of the time series; S is the standard deviation, that is, the square root of the change, representing the average trend of the time series (Hurst. 1951; Pang and North, 1996). The ratio of the range to the standard deviation (i.e., the rescaling range R/S) represents a dimensionless time series and the fluctuation intensity (Xiao et al., 2019a; Yang et al., 2020).

The concave section of the R(n)/S(n) curve can be regarded as the location of fracture development. To avoid recognition error, the R/S-FD method uses the finite difference method to calculate the second derivative of these discrete data (Xiao et al., 2019a). The zone with a positive second derivative can be identified as a fracture development section. F is the multiplication of the K value of the multigroup logging curve, which represents the accumulation of the fracture development probability of multigroup logging data identification.

The basic data of the R/S-FD method is a conventional logging curve. Firstly, the crack sensitivity curve suitable for the study area is selected by various methods. On the basis of referring to several commonly used fracture sensitivity curves, this study first uses the fracture development of the coring section to restrict the curve selection and selects three fracture sensitivity curves (DEN, AC, and CAL) to calculate the K value. Considering that the sampling interval of logging data is mostly 0.1 or 0.125 m, it is of certain indicative significance to optimize the fracture sensitivity curve by using whether there are fractures in the core section with a certain length (0.3–0.5 m). Then calculate the R/S value of the logging data of the study interval, and carry out the finite difference calculation (Figure 5). Finally, the calculated F value is used to identify the fracture development section.

FIGURE 5

Data Selection and Calculation

The following factors are considered in the selection of the numerical analysis method. The pre-Mesozoic Cambrian Sinian system in the study area is the main body of a large set of dolomites with single and continuous lithologic structures. Unlike tight sandstone/shale sections, the dolomite reservoir in the study area does not include lithology mutation surfaces such as sand-mud interbedding surfaces. The reservoir as a whole is a fractured dolomite reservoir with undeveloped pores (based on the 50.3 m coring data of four wells in the Sandaoqiao gas field and imaging logging data, the pores are undeveloped). The interference of holes is eliminated. The logging curve is complete, and the quality of the logging core data is good. The combination of R/S logging and core analysis is suitable.

Taking well Q1 as an example, the double logarithm curve of R/S value and serial number is made, and seven kinds of logging curves are selected for this analysis (Figure 6). The sensitivity of different good logs to fracture development varies. At the same time, the difference in lithofacies between study areas will make the analysis effect of different logging curves different.

FIGURE 6

Seven kinds of logging curves are selected in the analysis of wells Q1 and Q101: the CAL, AC, and DEN curves, which are widely used in R/S analysis, and other control curves RS, RD, SP, and GR. After comparing and analyzing the calculation results of the confirmatory fracture development section and logging curve in the coring section, CAL, AC and DEN are selected as the research objects (Table 1). That is, the F value used in this study is a product of the CAL, AC, and DEN K values (K_CAL/K_AC/K_DEN). In the fracture development section, the identification result with good applicability is obtained.

TABLE 1

Well NameCoring depth section/mCore length/mFracture developmentRecognition (corresponding ѻ/Non-corresponding ×)
ACDENCALRSRDSPGR
Q15,741.1–5,741.40.3YesѺѻѻѻѻ×ѻ
Q15,745.6–5,7460.4YesѺѻѻѻѻѻѻ
Q15,747.3–5,747.70.4YesѺѻѻѻѻ×ѻ
Q15,749–5,749.30.3YesѺѻѻ××××
Q15,753.7–5,7540.3YesѺѻѻ×××ѻ
Q15,754.5–5,754.80.3YesѺѻѻ×××ѻ
Q15,756.1–5,756.40.3YesѺѻѻѻѻѻѻ
Q15,757.8–5,758.20.4YesѺѻѻ××ѻѻ
Q15,759.7–5,7600.3YesѺѻѻѻѻѻѻ
Q15,763.8–5,764.10.3YesѺѻѻ××××
Q1015,757.8–5,758.10.3YesѺѻѻѻѻѻ×
Q1015,760.6–5,761.10.5YesѺѻѻѻѻѻѻ

Recognition effect of various logging curves on the fracture development depth of the coring section.

Table 1 shows that the commonly used fracture sensitivity curves AC, DEN, and CAL are consistent in the identification of fracture development section in the core section. The other four kinds of curves show a relatively low recognition rate for the fracture development section.

The calculation method used in this study is adjusted on the basis of the R/S-FD method. Based on the particularity of the dolomite value, the value range is adjusted, and the K value calculated by the logging curve at this depth is smaller (this may be related to the particularity of dolomite logging curves, with a small number of especially large or small values). Because K represents the fluctuation in logging data, we homogenize the interval of k > 0 to 1-2, which makes the probability of fracture development accumulate positively. If k = 0, then F = 0; if K≠0, then F is the product of the K values of the CAL, AC, and DEN. Thus, the fracture development section can be clearly identified.

Results

Characteristics of the Fracture Development Section

The dolomite section of well Q1 has developed thick fractures (a 60 m section), which can be divided into five sections, with an average thickness of 6.8 m and an average fracture spacing of 5.6 m (Figure 7). The fracture development section of well Q1 is relatively thicker and denser. The fractures of the core of well Q1 are high-angle and half-filled. The vertical distribution of the fracture development section in well Q101 is relatively uniform (an 80 m section), which can be divided into 11 sections, and the thickness of these fractures is thinner than those in well Q1, with an average length of 2.5 m and an average spacing of 4.1 m (Figure 8). Relatively speaking, the fracture development section of well Q101 is thinner and more dispersed.

FIGURE 7

FIGURE 8

Correlation Between the Fracture Development Section and Actual Production Data

The R/S-FD method is used to analyze the fractured dolomite reservoir of the Cambrian Sinian system in the Sandaoqiao area. The results show that the comprehensive identification results of multiple logging curves are in good correspondence with the actual observation results of cored sections, traditional numerical logging interpretations, and gas logging displays.

The fracture development section identified by well Q1 shows a good correspondence with the reservoir section interpreted by traditional logging data. The identified fracture development section corresponds to logging interpretation reservoir section Nos. 11–15 (Figure 7). The lower gas logging section of the dolomite reservoir in the Q101 well has a good correspondence with the identified fracture development section, corresponding to gas logging section Nos. 19–25 and in the lower part of the identified fracture development section, lagging behind the fracture development section (Figure 8). This correspondence shows that the R/S-FD method has certain applicability in fractured dolomite reservoirs.

Discussion

Fracture Types

The fracture types in the fracture development section identified by the R/S-FD method are noteworthy information (Xiao et al., 2019a). Through the identification results of Q1 and Q101 wells, combined with core, imaging logging, and other data, it is considered that most of the identified fractures are medium high angle fractures, semi filled and unfilled. Part of the fractures generated under the compression of the early formation are completely filled in the later stage. Because the filling material is similar to the composition of the reservoir, it is less abnormal in density logging and caliper logging. The identification effect of this method on this kind of fracture is general.

Controlling Factors

The difference between the two wells in the characteristics of longitudinal fracture development sections is affected by many factors, mainly including lithology, structural position, structural stage, and so on (Zhao et al., 2019). In addition, compared with the two wells, well Q1 shows higher and more stable oil and gas production. This difference may be due to the greater thickness of the average fracture development section of well Q1.

Method Applicability

The R/S-FD method also shows its limitations in the application of different strata in different areas (Xiao et al., 2019a; Li et al., 2019; Yang et al., 2020). Firstly, the stratum thickness used to identify the fracture development section should not be too thick. With the increase of thickness, reservoir heterogeneity is also increasing, which will make the anomalies identified by the R/S method contain a variety of other information. Secondly, the lithology of the study interval should be unified as far as possible. The sudden change in lithology will bring complex changes to the logging curve, which will reduce the accuracy of R/S recognition. In the application of the dolomite stratum, it is generally considered that a thickness less than 100 m is more suitable. In the fracture identification of the 300 m reservoir section of well Q102, the results show that the spacing of fracture development sections is too large, and there is no good corresponding result in the comparison with the coring section and production means.

Conclusion

1. In this study, the adjusted R/S-FD method is applied to fractured dolomite reservoirs in the deep layer (5,700–6200 m) of the Cambrian Sinian system (60–80 m), and good recognition results are obtained. In the process of applying this method to the study area, by comparing the fracture development identification results of continuous and complete coring sections, three logging curves with high fracture sensitivity (AC/DEN/CAL) are selected.

2. The adjusted R/S-FD analysis method can effectively identify the fracture development section of the fractured dolomite reservoir. Among them, five fracture development sections were identified in well Q1, with an average thickness of 6.8 m. The fracture development section is in good consistency with the reservoir interpretation section of conventional logging. Well, Q101 identified 11 fracture development sections with an average thickness of 2.5 m. The results show that the gas logging section lags behind the fracture development section, mainly corresponding to the lower part of the identified fracture development section. The difference between the two wells in the characteristics of longitudinal fracture development sections is affected by many factors, mainly including lithology, structural position, structural stage, and so on.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.

Author contributions

WD and QM are responsible for the idea and writing of this study and XD, PH, HW, and ZX are responsible for the logging interpretation.

Funding

This research was supported by the National Natural Science Foundation of China (Grant No. 42072173) and the National Sciences and Technology Major Project of China (2016ZX05046-003-001).

Acknowledgments

The authors thank the staff of all of the laboratories that cooperated in performing the tests and analyses. We are also grateful to the reviewers, whose comments improved the quality of this manuscript.

Conflict of interest

Author XD was employed by the Sinopec Northwest Oilfield Company.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

  • 1

    AfshariM. J.SomogyváriM.ValleyB.JalaliM.LoewS.BayerP. (2018). Fracture Network Characterization Using Stress-Based Tomography. J. Geophys. Res. Solid Earth123, 9324–9340. 10.1029/2018JB016438

  • 2

    AghliG.SoleimaniB.Moussavi-HaramiR.MohammadianR. (2016). Fractured Zones Detection Using Conventional Petrophysical Logs by Differentiation Method and its Correlation with Image Logs. J. Petroleum Sci. Eng.142, 152–162. 10.1016/j.petrol.2016.02.002

  • 3

    AntonelliniM.MollemaP. N. (2000). A Natural Analog for a Fractured and Faulted Reservoir in Dolomite: Triassic Sella Group, Northern italy. Aapg Bull.84 (3), 314–344. 10.1306/c9ebcddd-1735-11d7-8645000102c1865d

  • 4

    BatesC. R.LynnH. B.SimonM. (1999). The Study of a Naturally Fractured Gas Reservoir Using Seismic Techniques. AAPG Bull.83 (9), 1392–1406. 10.1306/e4fd41c5-1732-11d7-8645000102c1865d

  • 5

    ChenG. B.LiT.YangL.ZhangG. H.LiJ. W.DongH. J. (2021). Mechanical Properties and Failure Mechanism of Combined Bodies with Different Coal-Rock Ratios and Combinations. J. Min. Strata Control Eng.3 (2), 023522. 10.13532/j.jmsce.cn10-1638/td.20210108.001

  • 6

    DaviesG. R.SmithL. B. (2006). Structurally Controlled Hydrothermal Dolomite Reservoir Facies: an Overview. Bulletin90 (11), 1641–1690. 10.1306/05220605164

  • 7

    DingW.LiC.LiC.XuC.JiuK.ZengW.et al (2012). Fracture Development in Shale and its Relationship to Gas Accumulation. Geosci. Front.3 (001), 97–105. 10.1016/j.gsf.2011.10.001

  • 8

    DingW.ZhuD.CaiJ.GongM.ChenF. (2013). Analysis of the Developmental Characteristics and Major Regulating Factors of Fractures in Marine-Continental Transitional Shale-Gas Reservoirs: A Case Study of the Carboniferous-Permian Strata in the Southeastern Ordos Basin, Central China. Mar. Petroleum Geol.45, 121–133. 10.1016/j.marpetgeo.2013.04.022

  • 9

    Fernández-IbáñezF.DeGraffJ. M.IbrayevF. (2018). Integrating Borehole Image Logs with Core: A Method to Enhance Subsurface Fracture Characterization. Bulletin102 (06), 1067–1090. 10.1306/0726171609317002

  • 10

    GuoR.ZhangS.BaiX.WangK.SunX.LiuX. (2020). Hydrothermal Dolomite Reservoirs in a Fault System and the Factors Controlling Reservoir Formation-A Case Study of Lower Paleozoic Carbonate Reservoirs in the Gucheng Area, Tarim Basin. Mar. Petroleum Geol.120, 104506. 10.1016/j.marpetgeo.2020.104506

  • 11

    HanQ.LiZ.WangC.YangZ.-C.YanL.MengQ.-L.et al (2016). Characteristics and Age of a Proterozoic Buried Hill Reservoir in the Northern Shaya Uplift of the Tarim Basin. Geol. China43 (2), 486–499. (in Chinese with English abstract). 10.12029/gc20160210

  • 12

    HanQ.LiZ.YangZ.YanL.ShiY.SuJ. (2015). Characteristics of Hydrocarbon Accumulation in the Sandaoqiao Area of Tarim Basin. Petroleum Geol. Recovery Effic.22 (006), 14–20. (in Chinese with English abstract). 10.3969/j.issn.1009-9603.2015.06.003

  • 13

    HeX.ZhangP.HeG.GaoY.LiuM.ZhangY.et al (2020). Evaluation of Sweet Spots and Horizontal-Well-Design Technology for Shale Gas in the Basin-Margin Transition Zone of Southeastern Chongqing, SW China. Energy Geosci.1 (3–4), 134–146. 10.1016/j.engeos.2020.06.004

  • 14

    HongD.CaoJ.WuT.DangS.HuW.YaoS. (2020). Authigenic Clay Minerals and Calcite Dissolution Influence Reservoir Quality in Tight Sandstones: Insights from the Central Junggar Basin, NW China. Energy Geosci.1 (1–2), 8–19. 10.1016/j.engeos.2020.03.001

  • 15

    HurstH. E. (1951). Long Term Storage Capacity of Reservoirs. Am. Soc. Civ. Eng. Tans116, 776–808. 10.1061/taceat.0006518

  • 16

    LaiJ.ChenK.XinY.WuX.ChenX.YangK.et al (2021). Fracture Characterization and Detection in the Deep Cambrian Dolostones in the Tarim Basin, china: Insights from Borehole Image and Sonic Logs. J. Petroleum Sci. Eng.196, 107659. 10.1016/j.petrol.2020.107659

  • 17

    LaiJ.FanZ.WangZ.ChenJ.ZhouZ.et al (2017). Fracture Detection in Oil-Based Drilling Mud Using a Combination of Borehole Image and Sonic Logs. Mar. Petroleum Geol.84, 195–214. 10.1016/j.marpetgeo.2017.03.035

  • 18

    LanS. R.SongD. Z.LiZ. L.LiuY. (2021). Experimental Study on Acoustic Emission Characteristics of Fault Slip Process Based on Damage Factor. J. Min. Strata Control Eng.3 (3), 033024. 10.13532/j.jmsce.cn10-1638/td.20210510.002

  • 19

    LiA.DingW.LuoK.XiaoZ.WangR.YinS.et al (2019). Application of R/S Analysis in Fracture Identification of Shale Reservoir of the Lower Cambrian Niutitang Formation in Northern Guizhou Province, South China. Geol. J.55, 4008–4020. 10.1002/gj.3648

  • 20

    LiH.QinQ.ZhangB.GeX.HuX.FanC.et al (2020). Tectonic Fracture Formation and Distribution in Ultradeep Marine Carbonate Gas Reservoirs: A Case Study of the Maokou Formation in the Jiulongshan Gas Field, Sichuan Basin, Southwest China. Energy fuels.34 (11), 14132–14146. 10.1021/acs.energyfuels.0c03327

  • 21

    LiH. (2022). Research Progress on Evaluation Methods and Factors Influencing Shale Brittleness: A Review. Energy Rep.8, 4344–4358. 10.1016/j.egyr.2022.03.120

  • 22

    LiJ. Z.LaubachS. E.GaleJ. F. W.MarrettR. A. (2018). Quantifying Opening-Mode Fracture Spatial Organization in Horizontal Wellbore Image Logs, Core and Outcrop: Application to Upper Cretaceous Frontier Formation Tight Gas Sandstones, usa. J. Struct. Geol.108, 137–156. 10.1016/j.jsg.2017.07.005

  • 23

    LiY.ZhouD.WangW.JiangT.XueZ. (2020). Development of Unconventional Gas and Technologies Adopted in China. Energy Geosci.1 (1–2), 55–68. 10.1016/j.engeos.2020.04.004

  • 24

    LiuJ.DingW.XiaoZ.DaiJ. (2019). Advances in Comprehensive Characterization and Prediction of Reservoir Fractures. Prog. Geophys.34, 2283–2300. (in Chinese with English abstract). 10.6038/pg2019CC0290

  • 25

    LiuZ.ZhangY.ZhangY. (2022). Influencing Factor Analysis on the Fractured Tight Sandstone Gas Reservoir Characteristics: A Case Study of Bozi 3 Gas Reservoir in the Tarim Basin. Front. Earth Sci10, 881934. 10.3389/feart.2022.881934

  • 26

    LuoX.TangL.XieD. (2012). Mesozoic Bottom Boundary Unconformity and its Significance for Hydrocarbon Exploration in the Yakela Faulted uplift,Tarim Basin. Oil Gas Geol.33 (1), 30–36. (in Chinese with English abstract). 10.11743/ogg20120104

  • 27

    MaY.GuoT.ZhaoX.CaiX. (2008). The Formation Mechanism of High-Quality Dolomite Reservoir in the Deep of Puguang Gas Field. Sci. China Ser. D-Earth Sci.51 (1), 53–64. 10.1007/s11430-008-5008-y

  • 28

    MandelbrotB. B.WallisJ. R. (1969). Robustness of the Rescaled Range R/s in the Measurement of Noncyclic Long Run Statistical Dependence. Water Resour. Res.5 (5), 967–988. 10.1029/WR005i005p00967

  • 29

    MirandaJ. G. V.AndradeR. F. S. (1999). Rescaled Range Analysis of Pluviometric Records in Northeast brazil. Theor. Appl. Climatol.63 (1-2), 79–88. 10.1007/s007040050094

  • 30

    PangJ.NorthC. P. (1996). Fractals and Their Applicability in Geological Wireline Log Analysis. J. Pet. Geol.19 (3), 339–350. 10.1111/j.1747-5457.1996.tb00438.x

  • 31

    PrioulR.JockerJ. (2009). Fracture Characterization at Multiple Scales Using Borehole Images, Sonic Logs, and Walkaround Vertical Seismic Profile. Bulletin93, 1503–1516. 10.1306/08250909019

  • 32

    SantosR. F. V. C.MirandaT. S.BarbosaJ. A.GomesI. F.MatosG. C.GaleJ. F. W.et al (2015). Characterization of Natural Fracture Systems: Analysis of Uncertainty Effects in Linear Scanline Results. Bulletin99 (12), 2203–2219. 10.1306/05211514104

  • 33

    SongD.WangT.-G.LiH. (2015). Geochemical Characteristics and Origin of the Crude Oils and Condensates from Yakela Faulted-Uplift, Tarim Basin. J. Petroleum Sci. Eng.133, 602–611. 10.1016/j.petrol.2015.07.007

  • 34

    WangS.WangG.LaiJ.LiD.LiuS.ChenX.et al (2020). Logging Identification and Evaluation of Vertical Zonation of Buried Hill in Cambrian Dolomite Reservoir: A Study of Yingmai-Yaha Buried Hill Structural Belt, Northern Tarim Basin. J. Petroleum Sci. Eng.195, 107758. 10.1016/j.petrol.2020.107758

  • 35

    WangX.DingW.CuiL.WangR.HeJ.LiA.et al (2018). The Developmental Characteristics of Natural Fractures and Their Significance for Reservoirs in the Cambrian Niutitang Marine Shale of the Sangzhi Block, Southern china. J. Petroleum Sci. Eng.165, 831–841. 10.1016/j.petrol.2018.02.042

  • 36

    WangZ.XuK.ZhangH. (2022). Fracture Effectiveness Evaluation of Ultra-deep Tight Sandstone Reservoirs: A Case Study of the Keshen Gas Field, Tarim Basin, West China. Front. Earth Sci.10, 883479. 10.3389/feart.2022.883479

  • 37

    XiaoZ.DingW.HaoS.TaleghaniA. D.WangX.ZhouX.et al (2019b). Quantitative Analysis of Tight Sandstone Reservoir Heterogeneity Based on Rescaled Range Analysis and Empirical Mode Decomposition: A Case Study of the Chang 7 Reservoir in the Dingbian Oilfield. J. Petroleum Sci. Eng.182, 106326. 10.1016/j.petrol.2019.106326

  • 38

    XiaoZ.DingW.LiuJ.TianM.YinS.ZhouX.et al (2019a). A Fracture Identification Method for Low-Permeability Sandstone Based on R/S Analysis and the Finite Difference Method: A Case Study from the Chang 6 Reservoir in Huaqing Oilfield, Ordos Basin. J. Petroleum Sci. Eng.174, 1169–1178. 10.1016/j.petrol.2018.12.017

  • 39

    XueY.ChengL.MouJ.ZhaoW. (2014). A New Fracture Prediction Method by Combining Genetic Algorithm with Neural Network in Low-Permeability Reservoirs. J. Petroleum Sci. Eng.121 (2), 159–166. 10.1016/j.petrol.2014.06.033

  • 40

    YangH.PanH.WuA.LuoM.KonatéA. A.MengQ. (2017). Application of Well Logs Integration and Wavelet Transform to Improve Fracture Zones Detection in Metamorphic Rocks. J. Petroleum Sci. Eng.157, 716–723. 10.1016/j.petrol.2017.07.057

  • 41

    YangR.DingW.LiuJ.ZhaoZ.LiS.XiaoZ. (2020). Fracture Identification of Ordovician Carbonate Reservoir Based on R/S Analysis in the North of the Shunbei No. 5 Fault Zone, Tarim Basin. Interpretation8, T907–T916. 10.1190/int-2019-0284.1

  • 42

    YangY.TangL.DiaoX.XieD. (2018). Differential Deformation and its Control Mechanism of Fault Structures in Yakela Fault-Salient, Tarim Basin. Oil Gas Geol.039 (001), 89–97. (in Chinese with English abstract). 10.11743/ogg20180109

  • 43

    YinS.DongL.YangX.WangR. (2020). Experimental Investigation of the Petrophysical Properties, Minerals, Elements and Pore Structures in Tight Sandstones. J. Nat. Gas Sci. Eng.76, 103189–103214. 10.1016/j.jngse.2020.103189

  • 44

    YinS.LvD.DingW. (2018a). New Method for Assessing Microfracture Stress Sensitivity in Tight Sandstone Reservoirs Based on Acoustic Experiments. Int. J. Geomechanics18 (4), 1–16. 10.1061/(asce)gm.1943-5622.0001100

  • 45

    YinS.WuZ. (2020). Geomechanical Simulation of Low-Order Fracture of Tight Sandstone. Mar. Petroleum Geol.117, 104359–104416. 10.1016/j.marpetgeo.2020.104359

  • 46

    YinS.ZhaoJ.WuZ.DingW. (2018b). Strain Energy Density Distribution of a Tight Gas Sandstone Reservoir in a Low-Amplitude Tectonic Zone and its Effect on Gas Well Productivity: a 3d Fem Study. J. Petroleum Sci. Eng.170, 89–104. 10.1016/j.petrol.2018.06.057

  • 47

    ZhaoG.DingW.SunY.WangX.TianL.LiuJ.et al (2020). Fracture Development Characteristics and Controlling Factors for Reservoirs in the Lower Silurian Longmaxi Formation Marine Shale of the Sangzhi Block, Hunan Province, China. J. Petroleum Sci. Eng.184, 106470. 10.1016/j.petrol.2019.106470

  • 48

    ZhaoK. K.JiangP. F.FengY. J.SunX. D.ChengL. X.ZhengJ. W. (2021). Investigation of the Characteristics of Hydraulic Fracture Initiation by Using Maximum Tangential Stress Criterion. J. Min. Strata Control Eng.3 (2), 023520. 10.13532/j.jmsce.cn10-1638/td.20201217.001

  • 49

    ZhengH.WuM.WuX. (2007). Oil-gas Exploration Prospect of Dolomite Reservoir in the Lower Paleozoic of Tarim Basin. Acta Pet. Sin.28 (2), 1–8. 10.3321/j.issn:0253-2697.2007.02.001

  • 50

    ZhengH.ZhangJ.QiY. (2020). Geology and Geomechanics of Hydraulic Fracturing in the Marcellus Shale Gas Play and Their Potential Applications to the Fuling Shale Gas Development. Energy Geosci.1 (1–2), 36–46. 10.1016/j.engeos.2020.05.002

  • 51

    ZhengR. C.ZhongG, H. U.FengQ. P., (2007). Genesis of Dolomite Reservoir of the Changxing Formation of Upper Permian,northeast Sichuan Basin. J. Mineralogy Petrology, 27(4), 78–84. 10.1016/S1872-5813(08)60001-8

  • 52

    ZhongyingZ.YulinW.GuangdiL.XiangS. (2011). An Rescaled Range Analysis on the Characteristics of Coal Seam Development in the Eastern Depression of the Liaohe Basin. Min. Sci. Technol. (China)21 (02), 223–227. 10.1016/j.mstc.2011.02.011

  • 53

    ZhouX.WangR.DuZ. (2022). Characteristics and Main Controlling Factors of Fractures within Highly-Evolved Marine Shale Reservoir in Strong Deformation Zone. Front. Earth Sci.10, 832104. 10.3389/feart.2022.832104

  • 54

    ZhouY.YangF.JiY.ZhouX.ZhangC. (2020). Characteristics and Controlling Factors of Dolomite Karst Reservoirs of the Sinian Dengying Formation, Central Sichuan Basin, Southwestern china. Precambrian Res.343, 105708. 10.1016/j.precamres.2020.105708

  • 55

    ZhuG.MilkovA. V.LiJ.XueN.ChenY.HuJ.et al (2021). Deepest Oil in Asia: Characteristics of Petroleum System in the Tarim Basin, china. J. Petroleum Sci. Eng.199, 108246. 10.1016/j.petrol.2020.108246

Summary

Keywords

rescaled range analysis, R/S-FD method, fracture identification, fractured dolomite reservoir, sandaoqiao gas field

Citation

Meng Q, Ding W, Diao X, Han P, Wang H and Xiao Z (2022) Fracture Identification of Deep Dolomite Reservoir Based on R/S-FD Analysis: A Case Study of the Cambrian Sinian Reservoirs in the Sandaoqiao Gas Field, Northern Tarim Basin. Front. Earth Sci. 10:918683. doi: 10.3389/feart.2022.918683

Received

12 April 2022

Accepted

29 April 2022

Published

02 June 2022

Volume

10 - 2022

Edited by

Wei Ju, China University of Mining and Technology, China

Reviewed by

Kun Zhang, Henan Polytechnic University, China

Jianhua He, Chengdu University of Technology, China

Updates

Copyright

*Correspondence: Wenlong Ding,

This article was submitted to Structural Geology and Tectonics, a section of the journal Frontiers in Earth Science

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics