Abstract

The distribution of lower-order faults affects the development of oil and gas and the distribution of remaining oil, which is also the key to the development of fault block reservoirs in eastern China. The lower-order faults are characterized by short extension, small fault displacement, and difficult identification by traditional seismic interpretation methods. Therefore, it is necessary to predict the development law of faults in combination with the paleostress field during the fault formation period to improve the accuracy of fault seismic interpretation. In this paper, on the basis of reservoir structural analysis and rock mechanics experiments, the geological and mechanical model of the Sanduo period is established to predict the simulation of the paleostress field during the development of lower-order faults. The development of lower-order faults in rift basins is mainly controlled by the minimum principal stress and stress difference, and the shear stress in the profile controls the tendency of faults. Controlled by the paleostress field, the occurrence of lower-order faults in the Weizhuang area, Gaoyou Sag, and Subei Basin is very complex, and the strike of the fault is distributed as a tension shear broom on the plane. The prediction results can provide a reference for the excavation of remaining oil and the development of oil and gas.

1. Introduction

Complex fault block oilfields are widely distributed in China, with complex internal fault systems. Faults of different orders control the formation and distribution of oil and gas in fault block oilfields [13]. Higher-order faults control the structural trend and occurrence of the fault block, thus controlling the formation and accumulation of oil and gas. The existence of lower-order faults further complicates the fault block oilfield and is an important controlling factor for remaining oil enrichment [47]. With the exploration and development of oil fields, the excavation of remaining oil has become the primary content of increasing production and tapping potential of oil fields. Therefore, it is of great theoretical and practical significance to study the development law of lower-order faults. Any fault relative to the fault with a higher order can be called a lower-order fault, which specifically refers to the fourth- and fifth-order faults in oilfield development [1, 4, 7, 8]. The lower-order faults are characterized by short extension and small fault displacement, basically not controlling the accumulation of oil and gas, and it is difficult to identify accurately in seismic data, thus becoming a key point in the fine description of high water cut reservoirs. Because its formation is mainly controlled by the local tectonic stress field, analyzing the local regional stress field becomes an important method to predict lower-order faults [1, 810].

At present, some new seismic technologies, such as 3D seismic fine interpretation, fine coherence analysis, and seismic attribute technology, are mainly used to realize the identification and description of lower-order faults [3, 4, 7]. On the basis of identification and description, methods such as the structural background method, fault combination analysis method, structural stress field analysis method, structural physical simulation method, and rock stratum curvature method are mainly used to predict the development law of lower-order faults, thus guiding the implementation of local structures [1, 8]. Shale deformation can be determined through physical simulation and numerical simulation in combination with synchronous CT and thin section observations [1113]. Geomechanical modeling is an effective method for predicting the large-scale deformation of shale [1318]. According to 3D structural maps and rock mechanics experiments, a geomechanical model can be established to simulate the paleostress field in different periods [1, 2, 19]. In this paper, on the basis of 3D seismic interpretation, the authors apply the numerical simulation method of the tectonic stress field to simulate the distribution of paleotectonic stress and then predict the plane and profile development laws of lower-order faults from the perspective of stress. Finally, the latest structural map is used to reveal the genetic mechanism of lower-order faults.

2. Geological Setting and Stratigraphy

2.1. Geological Setting

The Weizhuang area is located west of the Hanliu fault zone in the Gaoyou Sag, Subei Basin (Figure 1), which is the stress release area of the Hanliu fault zone. The distribution of faults in this area is extremely complex, and it is a narrow fault nose structure group that is uplifted from south to north. The main oil-bearing series are the 2nd member (E1f2) and the 1st member of the Funing Formation (E1f1). It has been nearly 30 years since the region was put into development in 1996. At present, all of them have entered the stage of a medium-high water cut, and development is becoming increasingly difficult. The two challenges faced by oilfield development are as follows: the degree of oil recovery is increasing, the reserves are insufficient, and stable production requires the replacement of new fault blocks. However, due to the poor seismic data, it is difficult to implement the structure. Second, after multiple rounds of adjustment and densification, the well pattern in the main part is relatively complete, with little room for further improvement [8, 20, 21].

The region was developed in 1996 and has been developed for nearly 30 years. At present, all of them have entered the stage of high water content, and the difficulty of development has become increasingly difficult. The two challenges facing oilfield development are as follows: first, the degree of oil recovery is increasing, reserves are insufficient, and stable production requires new fault blocks to replace, but due to poor seismic data, it is difficult to implement the structure. Second, after multiple rounds of adjustment and densification, the well pattern in the main part is relatively complete, with little room for further improvement [8, 20, 21].

2.2. Stratigraphy

Gaoyou Sag has nearly 7000 m thick strata from the Upper Cretaceous to the Quaternary, which is the sag with the largest sedimentary thickness and the most fully developed strata in the Subei Basin. The Paleogene strata are the main oil generation strata and are mainly composed of continental clastic deposits (Figure 2) [8, 20, 21]. The Hanliu fault zone is developed in the sedimentary cover of the Gaoyou Sag. The Mesozoic Cenozoic strata cut through from bottom to top are the Paleogene Funing Formation, Dainan Formation, Sanduo Formation, Neogene Yancheng Formation, and Quaternary Dongtai Formation. The target horizon studied in this paper is the second member of the Funing Formation (E1f2). The lithology is mainly gray black and dark gray mudstone mixed with marlstone, thin limestone, dolomite, and a small amount of shale in the upper part, and light gray siltstone and gray black mudstone are interbedded in varying thickness in the lower part, with an apparent thickness of 33~300 m.

2.3. Tectonic Evolution

From the first to fourth members of the Funing Formation, the fault activity gradually increased during the sedimentary period, and the number of active faults reached a maximum in the Dainan period; during the sedimentation period of the Sanduo Formation, the fault activity peaked, and the number of active faults in the deep depression zone decreased. Based on comprehensive analysis of stratum thickness, structural evolution profile, fault activity law, etc., the structural evolution of the Hanliu fault zone can be divided into four stages [8]: the sedimentary period of the first and third members of the Funing Formation is the stage controlled solely by the Zhen-1 fault (Figure 1); the sedimentary period from the fourth member of the Funing Formation to the second member of the Dainan Formation is the formation and development stage of the Hanliu fault; during the sedimentation period of the Sanduo Formation, the Hanliu fault zone was complicated and controlled by two faults (Neogene-Quaternary depression development period). During the sedimentation of the Funing Formation, the fault strike is generally close to the E-W direction, indicating that the direction of the minimum principal stress in the Funing Formation is close to the N-S direction. The stress direction of the sedimentary period of the Sanduo Formation is basically consistent with the orientation and stress direction of active faults during the sedimentary period of the Funing Formation.

Influenced by the interaction between the paleostress field and the faults in the surrounding areas, according to the observation and statistical results of core fractures, the fractures in the study area are mainly near east-west and northeast fractures. High-angle fractures and vertical fractures account for more than 80% of the total fractures. The true opening of macro fractures is mainly distributed between 0.3 mm and 0.6 mm, the linear fracture density is small, and the development degree of fractures between wells varies greatly. Most structural fractures are filled with calcite, and some oil shows and a small amount of argillaceous filling. Tensile, shear, and tensile-shear joints are developed (Figure 3).

3. Data and Methods

3.1. Numerical Simulation of the Stress Field

Numerical simulation is an effective method for analyzing tectonic stress fields, and finite element simulation is a more commonly used method [19, 22, 23]. The finite element method is a numerical solution method to approximate the general continuum problem. Its basic idea is to simplify the studied continuum into a discrete model composed of finite elements and then apply the computer to solve the numerical solution. A geological body is discretized into finite continuous elements, which are connected by nodes, and each element is given with its actual rock mechanical parameters. We convert the solution of the continuous field function in the study area to the solution of field function values at a finite number of discrete points (nodes), and the basic variables are displacement, strain, and stress [19, 2326]. According to the boundary stress conditions and node equilibrium conditions, establish and solve the equations with node displacement as the unknown quantity and the overall stiffness matrix as the coefficient. The interpolation function is used to obtain the displacement on each node, and then, the stress and strain values in each element are calculated [2732]. These elements are combined to calculate the tectonic stress field of the whole geological body. With the increase in the number of elements, the model is closer to the actual geological body, and the results of numerical simulation need to be measured and evaluated by certain inspection standards to judge the effectiveness and accuracy of simulation. However, the paleostress field is a stress field in a certain historical period. At present, there is no effective way to evaluate them, and there is no mature and reliable test standard for its results. In general, successful paleostress field simulation should at least reach the following two points: the simulated stress field should correspond well with the tectonic intensity of the simulated area, and the simulation results can be used to explain the nature and intensity of tectonic activity [8, 33]; the simulated stress value should conform to the objective geological laws, and the relationship between the stress value and the rock strength, burial depth, structural characteristics, and other factors should be fully considered to ensure the validity and accuracy of the simulation results of the stress field [3436].

3.2. Calculation Method of Rock Strain Energy

According to the theory of elasticity, strain energy accumulates in the solid during deformation, and the strain energy in the solid can be measured by the strain energy density [37, 38]:

In equation (1), , , and are the principal stresses corresponding to the element and , , and are the corresponding stress directions. Both theory and experiments have proven that sandstone, siltstone, and shale show strong brittleness. The maximum strain energy criterion for brittle rock fracture considers that when the release rate of elastic strain energy accumulated inside is equal to the energy required to produce the surface of the fracture body per unit area, the brittle material will fracture [39, 40].

4. Results

4.1. Geologic Model

Through seismic interpretation, the fracture map and structural relief of shale tops in the study area are obtained. According to the bottom structure map of the second member of the Funing Formation, the model is simplified, and the preexisting high-order faults are selected to establish the geological model of the Weizhuang area (Figure 4). The main strike of the fault in the study area is NEE, and the dip angle of the fault is 45°. Thickness is a key factor affecting rock fracture [41] and is closely related to the activity and distribution of high-order faults [42]. According to the logging data, the model is determined as the sand mud interbed model, which is divided into five layers, and the thickness of the model is 150 m. The top and bottom sandstone layers are 40 m, the middle sandstone layer is 30 m, and the thickness of the two sets of mudstone layers is 20 m.

4.2. Geomechanical Model

The mechanical properties of different tectonic units are different. The strength of the general fault zone is weaker than that of the normal sedimentary strata [25, 28, 30]. The strength of the new strata in the sag area is greater than that of the old strata in the bulge area. According to the actual situation of the fault zone, the fault is assigned a mechanical parameter, and the rest of the sedimentary stratum is assigned a unified value of the mechanical parameter of the stratum. The size of the rock mechanical parameter is referred to in some experimental tests (Table 1). On the basis of loading mechanical parameters, the Solid45 element in ANSYS [23] is used for grid division, and 25986 elements are divided, with 16299 nodes participating in the operation (Figure 5). The Solid45 element is 8 nodes, each element is a cube, and nodes are 8 vertices on the cube. The external force acting on the cube element shall be loaded through the node. The external force does not act on the element directly but deforms the cube element through the node. The element is applicable to the numerical simulation of the stress field of a layered structure.

According to the analysis of paleostress and strain, the model loading is constantly changed. Finally, it is determined to apply 9.8 m/s2 gravity acceleration to the model during the Sanduo Period. Taking the rotation direction of the plane of the preexisting high-order fault in this period as the macroconstraint condition and meeting the requirements of finite element analysis and calculation, a tensile stress of 6 MPa is applied to the north-south boundary, a tensile stress of 2 MPa is applied to the east-west boundary, and a shear stress of 2 MPa is applied to the east boundary of the model to ensure the convergence of the calculation results and apply east-west linear constraints to the model frame (Figure 5). In addition, the direction displacement constraint of the bottom of the model is set to zero.

4.3. Numerical Simulation Results of the Stress Field

Using ANSYS software, the maximum principal stress distribution diagram and the minimum principal stress distribution diagram of the study area are obtained after solving. As shown in Figure 6, due to the strike-slip effect of the fault, the minimum principal stress is high in the hanging wall of the fault and low in the southeastern and southern parts of the study area. The intermediate principal stress is the nearly east-west principal stress. Due to the strike-slip action of the fault, the minimum principal stress is high in the hanging wall near the fault and low in the southwestern and southern parts of the study area.

5. Discussion

5.1. Development Degree of Lower-Order Faults

As shown in Figure 7, the difference between the maximum principal stress and minimum principal stress is high in the footwall of the fault, and the stress is 3-8 MPa; at the hanging wall of the fault, i.e., the southwest and south of the study area, the stress is low, and the stress is 0-1 MPa.

As shown in Figure 8, the distribution of strain energy is similar to the distribution of stress difference, showing high values in the footwall and low values in the hanging wall of the fault. In the area of high strain energy, the probability of rock fracture is higher; in contrast, the probability of rock fracture is small.

5.2. Development Law of the Lower-Order Fault Tendency

Figure 9 shows the distribution relationship between the lower-order fault tendency and the profile shear stress in the Weizhuang area, where negative values represent dextrorsal shear stress and south-dipping fault development areas; a positive value represents a sinistral direction and a north-dipping fault development area. The middle and eastern parts of the study area are generally inclined to the south, while the southwestern part is generally inclined to the north. The simulation results are basically consistent with the trend of lower-order faults, indicating that the shear stress of the profile controls the lower-order fault tendency.

5.3. Strike Development Law of Lower-Order Faults

As shown in Figure 10, the footwall of the Hanliu fault is of sinistral shear stress, with NNE-NE trending faults developed. The strike of low-order faults is nearly perpendicular to the minimum principal stress. The hanging wall of the fault is dextral shear stress, developing NE-trending faults. In the southwestern part of the study area, the footwall of the fault is dextral shear stress, developing NNE-trending faults.

As shown in Figure 11, the minimum horizontal principal stress in the Weizhuang area is near the north-south tensile stress. Because the rotation direction of plane shear stress near the fault changes, the stress near the fault deflects in different directions. At the footwall of the fault, the stress direction is NNE and nearly north-south at the footwall of the fault. Near the fault, the horizontal minimum principal stress direction rotates to a certain extent, and the stress direction turns parallel to the fault strike.

The method proposed by [1] is used to quantitatively predict the strike of small faults in this area. The intermediate principal stress is mainly used to control the strike of faults. In combination with the direction of paleostress, the strike of lower-order faults in the Weizhuang area is quantitatively predicted through the transformation of the stress coordinate system and geodetic coordinate system (Figure 12). Affected by the Hanliu fault, the Weizhuang area is located in the stress release area of the Hanliu fault. The strike of the lower-order fault is very complex, and it is distributed in a tension shear broom shape on the plane. The direction gradually changes from NNW to EW and NE. In the southwestern part of the study area, the strike of the lower-order fault changes to nearly NS.

6. Conclusions

Geomechanical modeling is an effective means for numerical simulation of paleostress fields and quantitative prediction of lower-order faults. According to the strike of active faults in different periods, the direction of paleostress and fault system can be determined, and the geomechanical model in different periods can be established to effectively predict the paleostress field in different periods.

The development degree of lower-order faults is mainly controlled by the minimum principal stress, strain energy and stress difference. The occurrence of faults is closely related to the stress type and planar and section shear stresses. The shear stress of the profile controls the tendency of the fault.

In the stress release zone near the fracture, the main stress direction is easy to deflect. Influenced by the Hanliu fault, the study area is located in the stress release area of the Hanliu fault, and the strike of lower-order faults is distributed in a tension shear broom shape on the plane. The direction of stress near the fault turns parallel to the fault strike.

Data Availability

The data that support the findings of this study are available from the corresponding author, Jingshou Liu.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was supported by the project ZR2020QD035 of the Shandong Provincial Natural Science Foundation, the “CUG Scholar” Scientific Research Funds at China University of Geosciences (Wuhan) (project no. 2022046), and the Open Research Fund of State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Grant No. SKLGME021023.