Paleostress Reconstruction from 3D Seismic Data and Slip Tendency in the Northern Slope Area of the Bongor Basin, Southwestern Chad

Hakro Ahmed Suhail^1*, Xintao Yuan^2,3, Jianlin Li¹, Hong Liu¹, Yuan Liu¹, Cong Ma¹

¹Geophysical Research Institute, BGP/CNPC, Hebei, China.
²CNPC International (Chad) Co., Ltd., N’Djamena, Chad.
³Research Institute of Petroleum Exploration and Development, Beijing, China.

DOI: 10.4236/ojg.2023.135024                

Abstract

Paleostress plays a significant role in controlling the formation, accumulation, and distribution of reservoirs, and this could be an important factor in controlling the production of hydrocarbons from the unconventional reservoirs. In this study, we will use 3D seismic reflection data to perform the slip-tendency-based stress inversion to determine the stress field in the basement of the northern slope area in the Bongor Basin. The dataset for this technique is easily available in the oil and gas companies. The stress inversion results from the basement of the northern slope area of Bongor basin show that the maximum principal stress axis (σ1) is oriented vertically, the intermediate principal stress axis (σ2) is oriented N18° and the minimum principal stress axis (σ3) is oriented N105°, and σ2/σ1 = 0.60 and σ3/σ1 = 0.29. The findings of this paper provide significant information to understand the fault reactivation at the critical stage of hydrocarbon accumulation and the regional tectonic evolution.

Keywords

Seismic Reflection Data, Slip Tendency, Bongor Basin, Stress Inversion

1. Introduction

Faults and fractures are significant factors in subsurface fluid current in complex geological structures [1] . Fault reactivation is a hot topic in the development of reservoirs [2] , and natural fractures are employed to anticipate the favorable zone for oil and gas exploration [3] [4] [5] [6] [7] . In general, both faults and fractures are considered to be fractures [8] . During the tectonic evolution, rocks may have been experienced a wide range of stress magnitude and orientations, the resulting in rock failure and fractures; thus, the reconstruction of the paleostress plays an important role in precise fracture prediction. In addition, paleostress reconstruction is a crucial approach for studying a region’s tectonic evolution [9] - [21] . Paleostress field investigation is useful to explain fault reactivation [22] [23] , and the timing and designs of fault leakage [22] [23] [24] [25] . Nowadays, a couple of methods have been proposed for estimating Paleostress analysis, such as earthquake focal mechanisms [26] , fault-slip data inversion [27] - [34] . However, these techniques are based on the outcrops and oriented core data. Such data is not easily available in the petroleum companies. Geologists are frequently concerned about the subsurface geological structures, and the data can be easily obtained from the 3D seismic interpretation, which usually does not relate to the field outcrop data.

In this paper, we will use the stress inversion method based on 3D seismic interpretation data, and slip tendency analysis. The method was used to reconstruct the Paleostress tensor in the basement of the northern Slope area of the Bongor Basin. The inversion results can be used to examine regional tectonic evolutions and are an important source of information for determining fault reactivity at the critical stage in the hydrocarbon accumulation.

2. Background of the Study Area

The Bongor Basin is located in the southwestern Chad, north of the middle section of the Central African Shear Zone (CASZ) (Figure 1(a)). On the other hand, this basin is an intracontinental rift basin that formed throughout the Mesozoic and Cenozoic under the influence of the (CASZ). The formation and evolution of this basin are closely related to the evolution of the African Plate, especially the shear zone in Central and West Africa, with typical passive rift characteristics [35] [36] . The Bongor basin covers an area of approximately 18,000 km². From east to west, this Basin extends for ~280 km long and 40 - 80 km wide, and is divided into four structural units based on the Bouguer gravity, aeromagnetic, and seismic data: the northern slope, central sub-basin, southern sub-basin, and southern uplift (Figure 1(b)) [37] . The central sub-basin is further separated into four depressions Annona Depression, Cola Depression, Mango Depression, and Moul Depression.

The basement of the Bongor basin is composed of Precambrian metamorphic rock, similar to the basement of the Sudan Rift Basin in nature [38] [39] [40] , on which tens of thousands of meters of Mesozoic and Cenozoic continental clastic rock strata are deposited, including the Lower Cretaceous, Tertiary, and Quaternary (Figure 2). According to the characteristics of gravity, aeromagnetic and seismic reflection wave groups and the distribution of regional unconformity, the sedimentary cover in the basin is further divided into two major structural layers: the Lower Cretaceous and the Cenozoic. These structural layers reflect that the tectonic evolution of the Bongor Basin underwent major rift activity,

Figure 1. Location map of study area (b) Simplified structural map of the Bongor basin. Modified after [37] . Oil/gas fields: 1-Ronier; 2-Mimosa; 3-Baobab; 4-Phoenix; 5-Phoenix S; 6-Raphia S-8; 7-Raphia S-6; 8-Daniela; 9-Lanea; 10-Lanea E; 11-Pavetta; 12-Moul; 13-Mango; 14-Vitex; 15-Delo. WARS: West African Rift System; CARS: Central African Rift System; EARS: East African Rift System; CASZ: Central African Shear Zone.

strong inversion, and a subsequent depression subsidence stage. The lower Cretaceous is the main period for the development of the rift, during which the main source and reservoir rock of the basin were deposited. During the lower Cretaceous, Bongor Basin underwent rapid subsidence in the early (Prosopis and Mimosa Formation sedimentary periods) and slow subsidence in the late (Kabula and Baobab Formation sedimentary periods) stages. About 85 - 80 Ma ago, the Santonian compression event was caused by the collision of the African and the Eurasian plate.

Figure 2. Regional stratigraphy of the Bongor basin showing the petroleum system and tectonic phases. Modified after [47] .

The reversal event, especially in the Bongor Basin is very strong, with the overall uplift of the basin being strong in the west and weak in the east, resulting in a lack of Upper Cretaceous in the basin. Entering the Cenozoic era, the regional stress field underwent significant changes, exhibiting strong extension in the West African Rift System and Sudan region, while thermal subsidence was dominant in the Bongor Basin. Correspondingly, the Termit Basin of the West African Rift System experienced strong subsidence, with a thickness of 3000 meters in the Tertiary strata, while the Bongor Basin had a maximum thickness of fewer than 500 meters in the Tertiary sediments. Therefore, the reservoir formation conditions of the Bongor Basin are similar to the basins of the West African Rift (WAR) and Sudan. There are significant differences between the Muglad Basin and the Melut Basin, as well as the Doba Basin in southern Chad.

The Bongor Basin is dominated by half-graben structures. Main boundary faults are generally striking in the ENE or WSW (Figure 3) and the dip angle of boundary faults and major faults varies from 30˚ to 50˚ but is most prominent at 45˚ in the SE and NW. The structural style in the Bongor basin is the large-scale domino fault blocks formed during the lower Cretaceous expansion and the transpressional inversion anticlines that formed on the hanging wall. The extensional topography has been intact; however, there has been an uplift of around 1500 m, with an 8% shortening ratio for the entire basin and, 20% for the half-grabens of the northern slope. A locally or regionally inverted rift basin would be a good way to describe the Bongor basin [41] [42] with various inversion styles and intensities in different regions of the basin.

3. Methodology

3.1. Slip Tendency Analysis

The idea behind slip tendency analysis is that the resolved shear stress and normal stress on a surface are excellent predictors of both the likelihood and the slip direction on the surface [43] . The approach has been used successfully to define fault slip and fault directions [44] [45] [46] . High slip-tendency oriented fractures are frequently more effective flow conduits than low slip-tendency orientation fractures [1] [43] [47] [48] . Stress anisotropy has the greatest impact when a

Figure 3. Seismic section illustrating the fault characteristic in the study area, location of seismic section shown in Figure 4.

Figure 4. Structural map of the basement of the northern slope of Bongor basin.

crack or fracture’s effective stress conditions are close to those needed for slip [49] [50] [51] . However, preferential fluid flow in fault and fracture routes becomes more obvious when the differential stress increases and the region of both faults and fractures with higher slip propensity increases [1] [43] [52] [53] .

In general, Equation (1) defines the slip tendency ( $Ts$ ) of a fault as the shear stress ( $\tau$ ) to the normal stress ( $\sigma n$ ) on the surface [43] :

$Ts=\frac{\tau }{\sigma n}$ (1)

In this study, we derived the equation of extensional regimes. For example, there are three principal stresses σ1 > σ2 > σ3 in space, the Anderson model suggests, there is a major stress axis that is perpendicular to the ground since the Earth’s surface is a principal stress plane.

In the extensional deformation regimes, the maximum principal stress axis (σ1) is oriented vertically, and the intermediate principal stress axis (σ2) and minimum principal stress axis (σ3) are horizontal. The estimate of stress tensor (maximum a posteriori (MAP) is also known as the maximum likelihood and least squares estimations. The relative magnitude of the principal stresses is for the numerical best-fit stress tensor.

These parameters are determined using a 3D stress computer program developed by the Southwest Research Institute (SWRI). The stress regime is defined by the nature of the vertical stress axes:

1) Extensional stress regime when σ2 is the maximum horizontal stress axis and σ1 is vertical.

2) Strike-slip stress regime when σ1 is the maximum horizontal stress and σ2 is vertical.

3) Compressional stress regime when σ1 is the maximum horizontal stress σ3 is vertical.

3.2. Slip Tendency-Based Stress Inversion

The slip tendency values can be used to calculate the likelihood and the fault slip direction on the fault surface. A higher slip tendency indicates the fault is more likely to reactivate, and that there will be more sliding along the fault with respect to time. Therefore, faults with higher slip tendency values accure a larger fault displacement than faults with lower slip tendency. However, there are a few exceptions, there is a fault with a higher slip tendency and slips during the end of the deformation history. As a result, a wide range of fault displacement values could be found in the faults with higher slip tendency. Therefore, the slip tendency can be used as a proxy for fault displacement [54] .

The stress inversion method characterized the degree of agreement between the computed slip tendency values of the stress tensor and the measured displacement values for a set of observed surfaces. In order to assess a candidate stress tensor’s ability to explain the pattern of displacement on the surface of different variable [43] .

4. Results

4.1. Data

The northern slope of the Bongor Basin covers an area of ~3200 km². 3D data for the whole area is available. The original data has been well processed; the seismic volume’s linear spacing is 25 m, the primary frequency is 15 - 30 Hz, and it has been loaded to a workstation and interpreted using GeoEast v.3.3.1. The structural maps of the basement (Figure 4) obtained from the time-depth conversion of 3D seismic interpretation data.

The geometric parameters of fault strike, dip angle and displacement in the basement were measured through 3D seismic interpretation. A total of n = 1123 fault data points were collected from the basement (Figure 4 and Table S1).

Fault Displacement was calculated from the seismic sections that were orthogonal to the direction of the adjacent fault. The horizontal and vertical distances between the footwall and hanging wall were used to compute the fault displacement.

A stereographic projection of faults with various displacement values is shown in Figure 5. The stereographic projection shows that the strike of the fault with larger displacement is in the NE-SW direction, but we can find some lower displacement values in the NW-SE direction, which could be the impact of the pre-existing fault structure in the basement. In general, high-angle faults typically accumulated higher displacement and are easier to reactivate. However, faults with the lower displacement are distributed in the SE-NW (Figure 5).

Figure 5. Fault displacement of basement. Stereographic projection (lower hemisphere, Schmidt projection) with fault poles.

4.2. Stress Inversion Results

A total of n = 1123 fault points were measured from the basement. The fault geometric parameters, i.e. strike, dip, and displacement which were obtained from the 3D seismic interpretation data to acquire the stress state condition under the slip tendency. In order to analyze this particular dataset, no artificial data are required. Table 1, describes the stress tensor’s maximum a posteriori (MAP) estimation.

The maximum and minimum slip tendency values in the data were recorded in 3D stress. Slip tendency results show in (Figure 6). For each fault point’s the fault displacement and slip tendency values were normalized [53] . The slip tendency values ranged from 0 to 1 allowing a direct comparison of the slip tendency and fault displacement. The Paleostress inversion determines the slip tendency that best corresponds to the fault displacement. The final stress inversion results show that the maximum principal stress axis σ1 is oriented vertically, the intermediate principal stress axis σ2 is oriented N18˚ and minimum principal stress axis σ3 is oriented N105˚, and σ2/σ1 = 0.60 and σ3/σ1 = 0.29.

Table 1. MAP Stress tensor estimate of basement.

Figure 6. Principal stress orientation and slip tendency plot of the basement of the northern slope area of the Bongor basin.

5. Discussion

In order to keep the production from the unconventional reservoir. Paleostress plays an important role in controlling and distributing of fractured reservoirs. Nowadays, couples of methods to invert Paleostress, but these methods are still dependent on the outcrop or oriented core data. Those data are generally not available in the oil and gas companies. The slip tendency-based stress inversion approach was previously demonstrated mainly by outcrop data [34] but it hasn’t been used with the seismic reflection data. Previously 3D seismic reflection data was used by [55] to perform the slip tendency based stress inversion through the slip tendency algorithm presented by [54] to identify the stress field of the Xicaogu area in the Bohai Bay Basin. However, we used the 3D Stress software package by the SWRI to perform the Slip tendency-based stress inversion to infer the stress field in the basement of the northern slope area in the Bongor basin. This method is based on the slip tendency and fault displacement and the dataset for this method can be easily collected from the oil and gas companies.

The stress inversion was accomplished based on the fault geometric data obtained from the seismic interpretation data. The fault surface orientation and displacement in the basement of the northern slope area of the Bongor basin were examined. The result shows that the fault displacements of each fault point exist in different areas of the stereographic projection and the faults with higher displacement values, steeper dip angle, and high slip tendency values are easier to reactivate than those faults with low slip tendency and displacement. This relationship provides a fault displacement that can be used as a proxy for slip tendency.

The stress state that best matched the slip tendency values and the fault displacement of the observed surface, the stress inversion was achieved based on the fault geometric data obtained from the 3D seismic interpretation data. The Paleostress inversion results show that the slip tendency has been found especially in the NE-SW direction.

Different stages of evolution were thought to have occurred in the Bongor Basin: a basement-forming period, and a strong setting period in the lower Cretaceous Prosopis and Mimosa formations. A Stable setting period in the lower Cretaceous Kabula, Ronier, and Baobab formation, and thermal precipitation during the Cenozoic. The Cretaceous rifting and Neogene post-rifting thermal subsidence were thought to have occurred in the Bongor Basin. The basin strike-slip regime is characterized by the west-northwest, east-southeast striking faults; these faults primarily controlled the deposition of lower Cretaceous sediments in the Bongor Basin, indicating that the faults may have been related to pre-existing basement faults and reactivated during the lower Cretaceous. This phenomenon might be connected to the basement faults reactivating as a result of the strike-slip movement of the Central African Shear Zone. The cumulative fault displacement during the tectonic or stress event was different, but the slip could accumulate on randomly oriented, potentially pre-existing faults.

The slip tendency-based stress inversion method opens a new and intriguing view point on paleostress analysis. Some modern software, such as Trap Tester, Move Midland Valley, can also be used to estimate the paleostress. This method is currently undergoing additional development and testing.

6. Conclusion

Based on the fault geometrics and displacement obtained from the 3D seismic interpretation data, stress inversion was then obtained by determining the stress state that best matched with the slip tendency values and the displacement of measures surface. The Paleostress inversion results show that the maximum principal stress axis (σ1) is oriented vertically, the intermediate principal stress axis (σ2) is oriented N18˚, the minimum principal stress axis and, (σ3) is oriented N105˚, and σ2/σ1 = 0.60 and σ3/σ1 = 0.29. The results show a strong correlation, indicating that most of the faults trend towards the northeast and southwest on the northern slope of the Bongor Basin, and could be related to pre-existing basement faults that were reactivated in the lower Cretaceous because of the strike-slip movement of the Central African Shear Zone.

Acknowledgements

This research work was supported by the China National Petroleum Corporation (No. CNPCIC-CT-2002-008).

Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Author Contributions

H.S.; X.Y.; contribute to conception and design of the paper. J.L.; contributes to program, H.L.; analyze the fault kinematics. Y.L.; C.M.; contribute to seismic interpretation. All authors contributed to manuscript revision, read and approved the submitted version.

Appendix

Table S1. Details of faults data collection points from basement. n = 1123 fault data points was selected to measure the strike, dip, horizontal slip, vertical fault throw and displacement. For location, see the Figure 4.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1]	Ferrill, D.A., Winterle, J., Wittmeyer, G., Sims, D., Colton, S., Armstrong, A. and Morris, A. (1999) Stressed Rock Strains Groundwater at Yucca Mountain, Nevada. GSA Today, 9, 1-7.
[2]	Moeck, I., Kwiatek, G. and Zimmermann, G. (2009) Slip Tendency Analysis, Fault Reactivation Potential and Induced Seismicity in a Deep Geothermal Reservoir. Journal of Structural Geology, 31, 1174-1182. https://doi.org/10.1016/j.jsg.2009.06.012
[3]	Magoon, L.B. and Dow, W. (1994) The Petroleum System—From Source to Trap. American Association of Petroleum Geologists, Tulsa. https://doi.org/10.1306/M60585
[4]	Jia, C., Zheng, M. and Zhang, Y. (2012) Unconventional Hydrocarbon Resources in China and the Prospect of Exploration and Development. Petroleum Exploration and Development, 39, 139-146. https://doi.org/10.1016/S1876-3804(12)60026-3
[5]	Wang, W., Pang, X., Chen, Z., Chen, D., Zheng, T., Luo, B., Li, J. and Yu, R. (2019) Quantitative Prediction of Oil and Gas Prospects of the Sinian-Lower Paleozoic in the Sichuan Basin in Central China. Energy, 174, 861-872. https://doi.org/10.1016/j.energy.2019.03.018
[6]	Gong, L., Wang, J., Gao, S., Fu, X., Liu, B., Miao, F., Zhou, X. and Meng, Q. (2021) Characterization, Controlling Factors and Evolution of Fracture Effectiveness in Shale Oil Reservoirs. Journal of Petroleum Science and Engineering, 203, Article ID: 108655. https://doi.org/10.1016/j.petrol.2021.108655
[7]	Zhang, L., Pang, X., Pang, H., Huo, X., Ma, K. and Huang, S. (2022) Hydrocarbon Accumulation Model Based on Threshold Combination Control and Favorable Zone Prediction for the Lower Enping Formation, Southern Lufeng Sag. Advances in Geo-Energy Research, 6, 438-450. https://doi.org/10.46690/ager.2022.05.08
[8]	Li, H., Yu, F., Wang, M., Wang, Y. and Liu, Y. (2022) Quantitative Prediction of Structural Fractures in the Paleocene Lower Wenchang Formation Reservoir of the Lufeng Depression. Advances in Geo-Energy Research, 6, 375-387. https://doi.org/10.46690/ager.2022.05.03
[9]	Jamison, W.R. (1992) Stress Controls on Fold Thrust Style. In: McClay, K.R., Ed., Thrust Tectonics, Springer, Berlin, 155-164. https://doi.org/10.1007/978-94-011-3066-0_13
[10]	Van Der Pluijm, B.A., Craddock, J.P., Graham, B.R. and Harris, J.H. (1997) Paleostress in Cratonic North America: Implications for Deformation of Continental Interiors. Science, 277, 794-796. https://doi.org/10.1126/science.277.5327.794
[11]	Zain Eldeen, U., Delvaux, U. and Jacobs, P. (2002) Tectonic Evolution in the Wadi Araba Segment of the Dead Sea Rift, South-West Jordan. EGU Stephan Mueller Special Publication Series, 2, 63-80. https://doi.org/10.5194/smsps-2-63-2002
[12]	Rowland, S.M., Duebendorfer, E.M. and Schiefelbein, I.M. (2007) Structural Analysis & Synthesis: A Laboratory Course in Structural Geology. 4th Edition, Wiley-Blackwell, Hoboken.
[13]	Soumaya, A., Ben Ayed, N., Rajabi, M., Meghraoui, M., Delvaux, D., Kadri, A., Ziegler, M., Maouche, S. and Braham, A. (2018) Active Faulting Geometry and Stress Pattern near Complex Strike-Slip Systems along the Maghreb Region: Constraints on Active Convergence in the Western Mediterranean. Tectonics, 37, 3148-3173. https://doi.org/10.1029/2018TC004983
[14]	Dasgupta, A., Mukherjee, S., Vanik, N., Chatterjee, R. and Pal, S.K. (2023) Paleostress Analyses and Oblique Rift Kinematics of Barmer Rift Basin, Western Rajasthan, India. Marine and Petroleum Geology. (Submitted)
[15]	Kumar, A., Shaikh, M.A., Singh, S., Singh. T., Mukherjee, S. and Singh, S. (2022) Active Morphogenic Faulting and Paleostress Analyses from the Central Nahan Salient, NW Siwalik Himalaya. International Journal of Earth Sciences, 111, 1251-1267. https://doi.org/10.1007/s00531-022-02176-3
[16]	Biswas, T., Bose, N., Dutta, D. and Mukherjee, S. (2022) Arc-Parallel Shears in Collisional Orogens: Global Review and Paleostress Analyses from the NW Lesser Himalayan Sequence (Garhwal Region, Uttarakhand, India). Marine and Petroleum Geology, 138, Article ID: 105530. https://doi.org/10.1016/j.marpetgeo.2022.105530
[17]	Maurya, D.M., Shaikh, M.A. and Mukherjee, S. (2021) Comment on Structural Attributes and Paleostress Analysis of Quaternary Landforms along the Vigodi Fault (VF) in Western Kachchh Region. Quaternary International, 601, 143-147. https://doi.org/10.1016/j.quaint.2021.04.029
[18]	Misra, A.A., Bhattacharya, G., Mukherjee, S. and Bose, N. (2014) Near N-S Paleo-Extension in the Western Deccan Region in India: Does It Link Strike-Slip Tectonics with India-Seychelles Rifting? International Journal of Earth Sciences, 103, 1645-1680. https://doi.org/10.1007/s00531-014-1021-x
[19]	Dutta, D., Biswas, T. and Mukherjee, S. (2019) Arc-Parallel Compression in the NW Himalaya: Evidence from Structural and Palaeostress Studies of Brittle Deformation from the Clasts of the Upper Siwalik, Uttarakhand, India. Journal of Earth System Science, 128, 125. https://doi.org/10.1007/s12040-019-1138-1
[20]	Goswami, T. K., Gogoi, M., Mahanta, B.N., Mukherjee, S., Saikia, H., Shaikh, M.A., Kalita, P., Baral, U. and Sarmah, R.K. (2022) Brittle Tectonics in the Western Arunachal Frontal Fold Belt, India: Change in Stress Regime from Pre-Collisional Extension to Collisional Compression. Geological Journal, 57, 5019-5038. https://doi.org/10.1002/gj.4393
[21]	Shaikh, M.A., Maurya, D.M., Mukherjee, S., Vanik, N., Padmalal, A. and Chamyal, L. (2020) Tectonic Evolution of the Intra-Uplift Vigodi-Gugriana-Khirasra-Netra Fault System in the Seismically Active Kachchh Rift Basin, India: Implications for the Western Continental Margin of the Indian Plate. Journal of Structural Geology, 140, Article ID: 104124. https://doi.org/10.1016/j.jsg.2020.104124
[22]	Gartell, A.P. and Lisk, M. (2005) Potential New Method for Paleostress Estimation by Combining Three-Dimensional Fault Restoration and Fault Slip Inversion Techniques: First Test on the Skua Field, Timor Sea. In: Boult, P. and Kaldi, J., Eds., Evaluating Fault and Cap Rock Seals, AAPG Hedberg Series, No. 2, American Association of Petroleum Geologists, Tulsa, 23-26.
[23]	Kulikowski, D. and Amrouch, K. (2018) 4D Modelling of Fault Reactivation Using Complete Paleostress Tensors from the Cooper-Eromanga Basin, Australia. Australian Journal of Earth Sciences, 65, 661-681. https://doi.org/10.1080/08120099.2018.1465472
[24]	Du Rouchet, J. (1981) Stress Fields, a Key to Oil Migration. American Association of Petroleum Geologists Bulletin, 65, 74-85. https://doi.org/10.1306/2F919774-16CE-11D7-8645000102C1865D
[25]	Henk, A. (2005) Pre-Drilling Prediction of the Tectonic Stress Field with Geomechanical Models. First Break, 23, 53-57. https://doi.org/10.3997/1365-2397.2005021
[26]	Gephart, J.W. and Forsyth, D.W. (1984) An Improved Method for Determining the Regional Stress Tensor Using Earthquake Focal Mechanism Data: Application to the San Fernando Earthquake Sequence. Journal of Geophysical Research, 89, 9305-9320. https://doi.org/10.1029/JB089iB11p09305
[27]	Etchecopar, A., Vasseur, G. and Daignieres, M. (1981) An Inverse Problem in Microtectonics for the Determination of Stress Tensors from Fault Striation Analysis. Journal of Structural Geology, 3, 51-65. https://doi.org/10.1016/0191-8141(81)90056-0
[28]	Angelier, J. (1984) Tectonic Analysis of Fault Slip Data Sets. Journal of Geophysical Research, 89, 5835-5848. https://doi.org/10.1029/JB089iB07p05835
[29]	Célérier, B. (1988) How Much Does Slip on a Reactivated Fault Plane Constrain the Stress Tensor? Tectonics, 7, 1257-1278. https://doi.org/10.1029/TC007i006p01257
[30]	Yamaji, A. (2000) The Multiple Inverse Method: A New Technique to Separate Stresses from Heterogeneous Fault-Slip Data. Journal of Structural Geology, 22, 441-452. https://doi.org/10.1016/S0191-8141(99)00163-7
[31]	Tezuka, K. and Niitsuma, H. (2000) Stress Estimated Using Microseismic Clusters and Its Relationship to the Fracture System of the Hijiori Hot Dry Rock Reservoir. Engineering Geology, 56, 47-62. https://doi.org/10.1016/S0013-7952(99)00133-7
[32]	Zalohar, J. and Vrabec, M. (2010) Kinematics and Dynamics of Fault Reactivation: The Cosserat Approach. Journal of Structural Geology, 32, 15-27. https://doi.org/10.1016/j.jsg.2009.06.008
[33]	Hansen, J.A. (2013) Direct Inversion of Stress, Strain or Strain Rate Including Vorticity: A Linear Method of Homogenous Fault-Slip Data Inversion Independent of Adopted Hypothesis. Journal of Structural Geology, 51, 3-13. https://doi.org/10.1016/j.jsg.2013.03.014
[34]	Tiwari, S.K., Beniest, A. and Biswal, T.K. (2020) Extension Driven Brittle Exhumation of the Lower-Middle Crustal Rocks, a Paleostress Reconstruction of the Neoproterozoic Ambaji Granulite, NW India. Journal of Asian Earth Sciences, 195, Article ID: 104341. https://doi.org/10.1016/j.jseaes.2020.104341
[35]	Guiraud, R. and Maurin, J.C. (1992) Early Cretaceous Rifts of Western and Central Africa: An Overview. Tectonophysics, 213, 153-168. https://doi.org/10.1016/0040-1951(92)90256-6
[36]	Genik, G.J. (1992) Regional Framework, Structural and Petroleum Aspects of Rift Basins in Niger, Chad and the Central African Republic (C.A.R.). In: Ziegler, P.A., Ed., Geodynamics of Rifting, Elsevier, Amsterdam, 169-185. https://doi.org/10.1016/B978-0-444-89912-5.50036-3
[37]	Dou, L., Xiao, K., Hu, Y., Song, H., Cheng, D. and Du, Y. (2011) Petroleum Geology and a Model of Hydrocarbon Accumulations in the Bongor Basin, the Republic of Chad. Acta Petrolei Sinica, 32, 379-386.
[38]	Tong, X.G., Dou, L., Tian, Z.J., Pan, X.H. and Zhu, X.D. (2004) Geological Mode and Hydrocarbon Accumulation Mode in Muglad Passive Rift Basin of Sudan. Acta Petrolei Sinica, 25, 19-24.
[39]	Tong, X.G., Xu, Z., Shi, B., Dou, L. and Xiao, K. (2006) Petroleum Geologic Property and Reservoir-Forming Pattern of Melut Basin in Sudan. Acta Petrolei Sinica, 27, 1-5.
[40]	Dou, L.R., Pan, X.H., Tian, Z.J., Xiao, K.Y. and Zhang, Z.W. (2006) Hydrocarbon Formation and Distribution of Rift Basins in Sudan—A Comparative Analysis of Them with Rift Basins in East China. Petroleum Exploration and Development, 33, 255-261.
[41]	Macgregor, D.S. (1995) Hydrocarbon Habitat and Classification of Inverted Rift Basins. Geological Society Special Publication, 88, 83-93. https://doi.org/10.1144/GSL.SP.1995.088.01.06
[42]	Chen, L., Ji, H., Dou, L., Du, Y., Xu, Z., Zhang, L., Yang, X. and Fu, S. (2018) The Characteristics of Source Rock and Hydrocarbon Charging Time of Precambrian Granite Reservoirs in the Bongor Basin, Chad. Marine and Petroleum Geology, 97, 323-338. https://doi.org/10.1016/j.marpetgeo.2018.06.003
[43]	Morris, A. Ferrill, D.A. and Henderson, D.B. (1996) Slip-Tendency Analysis and Fault Reactivation. Geology, 24, 275-278. https://doi.org/10.1130/0091-7613(1996)024<0275:STAAFR>2.3.CO;2
[44]	Streit, J.E. and Hillis, R.R. (2004) Estimating Fault Stability and Sustainable Fluid Pressures for Underground Storage of CO2 in Porous Rock. Energy, 29, 1445-1456. https://doi.org/10.1016/j.energy.2004.03.078
[45]	Lisle, R.J. and Srivastava, D.C. (2004) Test of the Frictional Reactivation Theory for Faults and Validity of Fault-Slip Analysis. Geology, 32, 569-572. https://doi.org/10.1130/G20408.1
[46]	Collettini, C. and Trippetta, F. (2007) A Slip Tendency Analysis to Test Mechanical and Structural Control on Aftershock Rupture Planes. Earth and Planetary Science Letters, 255, 402-413. https://doi.org/10.1016/j.epsl.2007.01.001
[47]	Barton, C.A., Zoback, M.D. and Moos, D. (1995) Fluid Flow along Potentially Active Faults in Crystalline Rock. Geology, 23, 683-686. https://doi.org/10.1130/0091-7613(1995)023<0683:FFAPAF>2.3.CO;2
[48]	Sibson, R.H. (2000) Fluid Involvement in Normal Faulting. Journal of Geodynamics, 29, 469-499. https://doi.org/10.1016/S0264-3707(99)00042-3
[49]	Stock, J.M., Healy, J.H., Hickman, S.H. and Zoback, M.D. (1985) Hydraulic Fracturing Stress Measurements at Yucca Mountain, Nevada, and Relationship to the Regional Stress Field (USA). Journal of Geophysical Research, 90, 8691-8706. https://doi.org/10.1029/JB090iB10p08691
[50]	Wiprut, D. and Zoback, M.D. (2002) Fault Reactivation, Leakage Potential, and Hydrocarbon Column Heights in the Northern North Sea. Norweigian Petroleum Society Special Publications, 11, 203-219. https://doi.org/10.1016/S0928-8937(02)80016-9
[51]	Zhang, X., Koutsabeloulis, N. and Heffer, K. (2007) Hydromechanical Modeling of Critically Stressed and Faulted Reservoirs. AAPG Bulletin, 91, 31-50. https://doi.org/10.1306/08030605136
[52]	Zoback, M., Barton, C., Finkbeiner, T. and Dholakia, S. (1996) Evidence for Fluid Flow along Critically-Stressed Faults in Crystalline and Sedimentary Rock. In: Jones, G., Fisher, Q. and Knipe, R., Eds., Faulting, Faults Sealing and Fluid Flow in Hydrocarbon Reservoirs, University of Leeds, London, 47-48.
[53]	Takatoshi, I. and Kazuo, H. (2003) Role of Stress-Controlled Flow Pathways in HDR Geothermal Reservoirs. Pure and Applied Geophysics, 60, 1103-1124. https://doi.org/10.1007/978-3-0348-8083-1_19
[54]	McFarland, J.M., Morris, A.P. and Ferrill, D.A. (2012) Stress Inversion Using Slip Tendency. Computers and Geosciences, 41, 40-46. https://doi.org/10.1016/j.cageo.2011.08.004
[55]	Ping, G., Liu, X., Li, M., Zhang, X., Gao, Y., Wang, S. and Liu, X. (2022) Using Seismic Data and Slip Tendency to Estimate Paleostress: A Case Study from Xicaogu Area of Bohai Bay Basin, China. Frontiers in Earth Science, 10, Article ID: 812874. https://doi.org/10.3389/feart.2022.812874