Home Application of the wavelet transform and Hilbert–Huang transform in stratigraphic sequence division of Jurassic Shaximiao Formation in Southwest Sichuan Basin
Article Open Access

Application of the wavelet transform and Hilbert–Huang transform in stratigraphic sequence division of Jurassic Shaximiao Formation in Southwest Sichuan Basin

  • Zheng Li , Jingchun Tian , Laixing Cai and Tian Yang EMAIL logo
Published/Copyright: September 3, 2024
Become an author with De Gruyter Brill

Abstract

In the southwestern Sichuan Basin, the Jurassic Shaximiao Formation encompasses a multitude of working areas, displaying intricate sedimentary traits. Traditional methods of stratigraphic division based on sequence suffer from inherent subjectivity and limitations. This study employs a combined mathematical approach to use the wavelet transform (WT) and the Hilbert–Huang transform (HHT). It decomposes the natural gamma ray (GR) logging curve into energy spectrum plots and wavelet coefficients at different scales, high and low frequency signals at different frequencies, and a set of intrinsic mode function components and residual functions. The study conducted a detailed stratigraphic division of the Jurassic Shaximiao Formation in the southwestern Sichuan Basin using these methods. The WT offers greater resolution for the periodic changes in the base level, whereas the HHT demonstrates a superior correlation with the positions of stratigraphic interfaces. The combined utilization of the continuous wavelet transform, the discrete wavelet transform, and HHT methods has demonstrated encouraging outcomes in the stratigraphic division of the Jurassic Shaximiao Formation. These methods have been shown to enhance the accuracy of stratigraphic division and to reduce the influence of subjective factors. This study presents new insights and approaches for geological data processing, offering significant theoretical and practical implications and novel technical means for oil and gas exploration and development.

1 Introduction

In recent years, China has made significant progress in exploring and developing of tight sandstone oil and gas resources. The proportion of tight sandstone oil and gas reserves has been increasing, with major deposits found in basins such as the Sichuan Basin, Ordos Basin, and Tarim Basin [1,2,3,4]. The Shaximiao Formation in the Middle Jurassic of the Sichuan Basin has become an important exploration target in recent years [5]. The daily gas production is 136 × 104 m3, which has great exploration potential and is the main replacement force of tight sandstone gas reservoirs in the Xujiahe Formation [6]. And the southwestern region of the Sichuan Basin has significant terrestrial tight gas resources that have not been fully explored, indicating high exploration potential [7,8]. Accurate stratigraphic partitioning is critical for successful oil and gas exploration, and can even determine the success or failure of gas field development [9]. Stratigraphic division in conventional methods is typically performed using data such as outcrops, core samples, and well logs. However, this process can be influenced by subjective human factors [10,11,12,13]. In the southwestern Sichuan Basin, the Shaximiao Formation covers a large area spanning multiple working zones. However, no unified sequence stratigraphy work has been conducted, leading to numerous stratigraphic division schemes. Some scholars propose that the Shaximiao Formation can be subdivided into three third-order sequences and five fourth-order sequences, while others propose that the formation can be divided into Upper Shaximiao Formation (J2s) and Lower Shaximiao Formation (J2x), or alternatively be divided into two sections: Shaximiao Formation Section 1 (J2s1) and Shaximiao Formation Section 2 (J2s2) based on lithological characteristics and thickness [14,15,16]. This has hindered a comprehensive understanding of the internal complex sand bodies, which has been severely impeding the progress of exploration and development [17]. Therefore, it is necessary to conduct more refined, high-resolution sequence stratigraphy work to establish a precise foundation for subsequent oil and gas exploration and development.

Logging curves can effectively reflect information about cycles of different periods in sedimentary strata and serve as compelling evidence for performing sequence stratigraphic division [18]. Related high-resolution sequence stratigraphic classification methods, such as the Fourier transform [19,20], the wavelet transform (WT) [21,22,23], and the Hilbert–Huang transform (HHT) [24,25], can be used in intelligent identification tasks within sequence stratigraphy. These methods can help reduce the impact of the aforementioned limitations on sequence stratigraphic work. The WT can identify both time and frequency information within signals, extracting the necessary information from well logging data. It provides high resolution for the periodic variations in sequence stratigraphic cycles on the base level [26]. The earliest application of the WT in stratigraphic partitioning dates back to 2003 [27]. Since then, an increasing number of scholars have applied the WT to sequence stratigraphic division work, including both the continuous wavelet transform (CWT) and the discrete wavelet transform (DWT) [28,29,30,31]. The HHT is a tool used to deal with complex, nonlinear, and nonstationary signals [32,33,34]. It has good identification capabilities for small-scale sedimentary bodies and concealed geological features, especially providing higher resolution for delineating the boundaries of different geological bodies [35,36]. Using of computational techniques to identify and compare reference surface cycles in well logging data can significantly reduce workload and human error [37].

The sedimentary characteristics of the Shaximiao Formation in the southwestern Sichuan Basin are complex, characterized by the development of shallow-water delta-river sedimentary systems. Sediments of shallow-water deltas and fluvial facies exhibit striking similarities in their characteristics and distribution. The complexity of the sedimentary environment, in conjunction with the extensive distribution area, presents a significant challenge to the conduct of stratigraphic work on target layers within the study area. It is evident that the stratigraphic division criteria established through conventional methods are unable to fully address this issue. Therefore, this study combines the advantages of the WT and HHT, conducts comparative analyses, and addresses the limitations of individual methods to accurately identify hidden sedimentary cyclical information within well logging curves. The approach enables more precise high-resolution sequence stratigraphic work.

2 Geological setting

The Sichuan Basin is situated in southwestern China (Figure 1a). It is a cratonic basin with multiple-cycle superimposed sedimentation, covering an area of approximately 18,000 km2. Ending the Late Triassic period, the Sichuan Basin transitioned into the Jurassic Red Basin stage. The Micangshan-Dabashan uplift occurred in the northeastern part of the basin during the deposition of the Shaximiao Formation, and served as the primary source of sediment during this period [38] (Figure 1b). The Sichuan Basin underwent marine sedimentation from the Ediacaran to the Middle Triassic of the Mesozoic Era, and it was mainly composed of carbonate rocks. During the Late Middle Triassic period, the Upper Yangtze region experienced uplift due to the Indosinian movement. This caused the cessation of marine transgression and led to erosion of the basin interior. By the Late Triassic, the sedimentary environment had shifted from shallow marine platforms to inland lake basins, marking a transition from marine to continental facies deposition. During the Jurassic, Cretaceous, and Paleogene periods, sedimentation was mainly continental, consisting of clastic rocks [39]. In the study area, the Middle Jurassic Lianggaoshan Formation is underlain by the Shaximiao Formation. The sedimentary environment of the first member of the Shaximiao Formation is similar to that of the Lianggaoshan Formation, characterized by the development of a deltaic-fluvial depositional system. Towards the end of the first member of the Shaximiao Formation, there is a brief period of littoral-lake sedimentation. Fluvial depositional environments dominate the second member of the Shaximiao Formation. The formations are mainly composed of sandstone, siltstone, and mudstone, with a thickness of approximately 1,200–2,000 m [40,41] (Figure 1c).

Figure 1 
               (a) Map of the Sichuan Basin’s location. (b) Location map of the study area in southwest Sichuan Basin. Location of the studied wells H, YQ, and ZT is shown on the map by filled red circles. (c) Schematic diagram of the regional stratigraphy of the study area in Jurassic Shaximiao Formation in southwest Sichuan Basin.
Figure 1

(a) Map of the Sichuan Basin’s location. (b) Location map of the study area in southwest Sichuan Basin. Location of the studied wells H, YQ, and ZT is shown on the map by filled red circles. (c) Schematic diagram of the regional stratigraphy of the study area in Jurassic Shaximiao Formation in southwest Sichuan Basin.

3 Methodology

3.1 Wavelet Transform

The WT is a sophisticated technique that builds on the Fourier Transform. Its fundamental concept involves breaking down signals into wavelet function series, allowing for the projection of signals in various time-frequency spaces. This approach addresses the limitation of poor time domain resolution in the Fourier Transform and demonstrates excellent localization characteristics in both time and frequency domains simultaneously. This makes it particularly suitable for processing well logging curves [42,43]. Through the WT, well logging curves can be transformed into one-dimensional wavelet coefficient and two-dimensional spectrograms, which reveal hidden information in the curves. This makes it useful for geological interpretation. The WT includes the CWT and the DWT.

The CWT possesses means the properties of time-frequency localization and adjustable time windows, enabling more precise determination of the time-frequency characteristics of well logging data. This capability facilitates the identification of different scales of sedimentary cycles based on these variations, allowing for more efficient sequence stratigraphic works [44,45]. The specific definition of the CWT is defined as follows:

(1) W f ( a , b ) = a 1 / 2 R f ( t ) ψ t b a d t ,

where “a” represents the scale parameter, “b” represents the time shift parameter, f(t) represents the well logging signal sequence, and “ψ” is the wavelet function. The WT adjusts the sampling step size in the time domain for different frequencies by varying “a” and “b.” The initial step involves comparing the selected wavelet ψ with the initial portion of the input signal f(t). This is followed by computing the correlation between the chosen wavelet ψ and the input signal f(t), shifting it to the right, and repeating this process until the entire input signal f(t) has been covered. Finally, by scaling (stretching) the wavelet, the process is repeated to accomplish continuous WT analysis. The CWT has high frequency resolution but poor time resolution at low frequencies. Conversely, the CWT has high time resolution at high frequencies and low frequency resolution. This characteristic allows for the analysis of slow changes in low-frequency signals and rapid changes in high-frequency signals, which is the basis for signal analysis using the CWT.

Because the scaling and translation coefficients of the CWT are independent of each other, there is some similarity between the different wavelet functions obtained by scaling and translation. However, due to the independence between these two coefficients, redundant signals can be generated. The DWT was introduced to reduce redundancy.

The DWT decomposes a signal into a set of mutually orthogonal discrete and shifted signals, obtaining corresponding high frequency and low frequency wavelet reconstruction coefficients. The low-frequency component tends to match the inherent characteristics of the signal itself, while the high-frequency component is often contaminated by noise. By decomposing the low-frequency component further, we can obtain different levels of wavelet reconstruction coefficients, which allows for a more precise signal analysis [42,46]. The specific definition of the DWT is defined as follows:

(2) Ψ j , k ( t ) = a 0 j 2 Ψ ( a 0 j t k b 0 ) ,

where Ψ j , k ( t ) represents the discrete wavelet function, a 0 and b 0 are discretized scale and shift parameters, respectively, with a = a 0 j and b = ka 0 j b 0, j ∈ Z, a 0 is always greater than 1, and k is an integer. First, the selected wavelet function ψ is matched with the characteristics of the input signal f(t). Second, the chosen wavelet function ψ is used as the core of the filter to perform convolution on the input signal f(t), decomposing the signal into low-frequency components (approximation coefficients) and high-frequency components (detail coefficients). Subsequently, the convolution result is down sampled in order to reduce the data volume while preserving the essential information, thereby generating a set of approximation coefficients and a set of detail coefficients representing signal components at different frequency ranges. Finally, the aforementioned steps are repeated, with the low-frequency part (approximation coefficients) being further decomposed until the desired decomposition level is achieved or the maximum depth of decomposition is reached. The inverse WT is employed to reconstruct the decomposed approximation coefficients and detail coefficients into an approximate value of the original signal, thereby enabling the analysis of signal characteristics at different scales.

3.2 Hilbert-Huang transform

The HHT differs significantly from traditional signal processing methods and is a time-frequency analysis technique that possesses unique advantages in high-resolution analysis of nonlinear signals. This method decomposes the original signal into intrinsic mode functions (IMF) of different frequencies using empirical mode decomposition (EMD). The HT is then applied to these IMFs to generate instantaneous frequencies that reflect the abrupt changes in the original signal. This approach effectively addresses issues such as signal distortion [47,48,49]. The specific definition of the HHT is defined as follows:

(3) H ( t ) = j = 1 n e j ( t ) + C n ( t ) ,

where e j ( t ) is defined as the IMF components at different scales, representing different frequency components that change according to the variations in the original signal. C n ( t ) is defined as the residual component function, indicating the trend of signal changes, where n denotes the level of signal decomposition. First, compute the local maxima and minima of the input signal f(t) and use these to construct the upper and lower envelopes of the original signal. The upper envelope is formed by connecting line segments between adjacent local maxima, while the lower envelope is formed by connecting line segments between adjacent local minima. This process yields an approximation of the upper and lower bounds of the original signal. Second, the upper and lower envelopes are added together and divided by two to obtain the mean function, which describes the overall trend of the original signal. The mean function is then subtracted from the original signal to obtain the residual signal, which reflects the local oscillatory component between the original signal and the mean function. Subsequently, it is necessary to ascertain whether the residual signal satisfies the convergence criteria of the IMF. This entails determining whether the residual signal is a monotonic or quasi-static function and whether the number of extremum points and zero crossings of the residual signal is equal or differs by no more than one. In the event that the aforementioned criteria are not met, it is necessary to repeat the aforementioned steps. Conversely, if the criteria are met, the residual signal must be defined as one of the IMFs and saved. Finally, the aforementioned steps should be repeated until the maximum decomposition limit is reached, with the aim of gradually capturing different frequency components within the signal.

3.3 Specific steps

The above signal processing methods can be implemented using MATLAB software’s Wavelet Toolbox, Main Menu toolbox, and HHT custom function code program. To perform time–frequency analysis on well logging curves, it is necessary to organize the well logging data of the target formation into text format. These data can then be processed using MATLAB software for analysis. Well logging curves can accurately reflect the superimposed oscillations of sedimentary cycles of different periods. By applying the WT and the HHT to well logging curves, it is possible to identify the periodic oscillatory characteristics of sedimentary cycles. These oscillations can then be matched with sequence interfaces of different orders, aiding in conducting more refined sequence stratigraphic work.

4 Application to the Shaximiao Formation

In the field of the HT analysis research, matching pursuit algorithm is used for selection of the most appropriate wavelet for signal processing [42]. Examining different wavelets revealed that Daubechies (db) wavelet is associated with the lowest mean squared error in signal approximation. The db wavelet was selected for analysis of well logging data based on the unique geological features of the research area. In principle, signals can be decomposed to an infinite degree after being input into the HHT [50]. However, following an analysis of the data from various well locations within the research area, it was determined that the most effective approach would be to employ seven IMFs and one residual signal (R) in order to conduct high-resolution stratigraphic work.

The Shaximiao Formation in the southwestern region of the Sichuan Basin can be divided into two third-order base level descents/rises and five fourth-order base level descents/rises. This corresponds to two third-order sequences and five fourth-order sequences, respectively. The GR logging curve of well YQ was analyzed using MATLAB software. A comprehensive interpretation of the upper boundary of the Shaximiao Formation (SB1), the lowest boundary of the Shaximiao Formation (SB2), secondary stratigraphic boundary (sb1, sb2, sb3), the maximum flooding surface (MFS1), and the secondary flooding surface (mfs1, mfs2, mfs3, and mfs4) was conducted for high-resolution sequence stratigraphy.

The GR curve of well YQ was subjected to the CWT using the db6 wavelet as the mother wavelet. This resulted in the corresponding time–frequency spectrogram and wavelet coefficient reconstruction curve (Figure 2). The time–frequency spectrogram shows regions of increasing color intensity from light to dark, indicating the increasing energy of sedimentary water bodies and representing the cyclicity of sedimentary formations. The accommodation space to sedimentation rate (A/S) ratio can effectively reflect the variations in the cycles of the reference surfaces [51]. When the A/S ratio is smaller, the sand bodies exhibit obvious overlapping and cutting characteristics, and the vertical accretion is stronger. Conversely, when the A/S ratio is large, isolated sandstone surrounded by mudstone often develops, resulting in poor lateral connectivity. Typically, a complete sequence boundary (SB) comprises a rising surface (A/S > 1) and a falling surface (A/S < 1). Regarding the time–frequency spectrogram, there is a clear transition of the strong energy cluster from the SB to the weak energy cluster at the MFS, followed by a transition back from the weak energy cluster at the MFS to the strong energy cluster at the SB, forming a complete stratigraphic cycle. Based on Figure 2, it is evident that SB1, SB2, sb2, and sb3 correspond well on the time–frequency spectrogram. Using sb2 as an example, there is a distinct mutation between strong and weak energy clusters, whereas the correspondence for sb1 is not evident. MFS1, mfs1, mfs3, and mfs4 also show good correspondence in the time-frequency spectrogram, while mfs2 does not correspond to a stronger energy cluster. The GR logging curve underwent denoising processing at three different frequencies: a = 32 (high frequency), a = 64 (medium frequency), and a = 128 (low frequency). The resulting wavelet reconstruction coefficient curves were obtained. Most interfaces correspond well with the abnormal regions of the wavelet reconstruction coefficient curves, providing valuable reference information. When using the CWT for sequence stratigraphy identification, it is important to note that higher frequencies reveal more details. In the CWT, not all areas of abrupt change in the time–frequency spectrogram and wavelet reconstruction coefficients can be explained as SBs and MFCs. Most of these correspond to specific areas of SBs. Therefore, combining geological knowledge and other evidence is necessary to identify and delineate stratigraphic SBs.

Figure 2 
               The results of applying CWT on the GR log of well YQ to interpret the stratigraphic cycles of the Jurassic Shaximiao Formation. MFS is the maximum flooding surface and SB is the sequence boundary.
Figure 2

The results of applying CWT on the GR log of well YQ to interpret the stratigraphic cycles of the Jurassic Shaximiao Formation. MFS is the maximum flooding surface and SB is the sequence boundary.

Based on Figure 3, continuing with the selection of db wavelet, the GR logging curve is decomposed into 11 high-frequency to low-frequency detail signals (d1–d11) and one approximation signal (a11). From d1–d11, it can be observed that there is a good correspondence between the anomaly regions and the MFSs. Using mfs2 as an example, a corresponding response at d5 and d7 is observed when the base level rises (A/S>1), while its response in the CWT is less prominent (Figure 2). In contrast, sb1 shows a clear correspondence at d5, d7, and d11, indicating a complete process of the base level transitioning from rising to falling.

Figure 3 
               The results of applying DWT on the GR log of well YQ to interpret the stratigraphic cycles of the Jurassic Shaximiao Formation. MFS is the maximum flooding surface and SB is the sequence boundary.
Figure 3

The results of applying DWT on the GR log of well YQ to interpret the stratigraphic cycles of the Jurassic Shaximiao Formation. MFS is the maximum flooding surface and SB is the sequence boundary.

As decomposition levels increase, the signal’s noise is gradually filtered out, allowing for better identification of MFSs. The a11 curve shows a similar fluctuation trend compared to the original GR well logging data. After removing significant noise, it displays clearer and more intuitive peaks and troughs.

Generally, the a11 signal’s troughs correspond to the opposing sides of the base level descent (A/S < 1) and base level ascent (A/S > 1). In contrast, the peaks correspond to the transition surfaces from base level ascent (A/S > 1) to base level descent (A/S < 1). Overall, the use of DWT provides a more intuitive recognition of the transition surfaces (MFSs) between base level ascent (A/S > 1) and base level descent (A/S < 1) compared to the response of the CWT. However, it also exhibits some degree of offset. For instance, taking MFS1 as an example, there is a noticeable downward shift. Similar to the CWT, the DWT also limits the identification of sequence boundaries and flooding surfaces to a relatively small range, rather than providing precise responses for accurate interface positioning.

The results of processing the GR logging curve using HHT are shown in Figure 4. The GR logging curve of well YQ was decomposed into seven IMFs (IMF1–IMF7) and one residual function (R7) through EMD. A good correspondence between various sequence boundaries and maximum flooding surfaces was observed in IMF3 and IMF7. Using MFS1, mfs1, mfs2, and mfs4 as examples, there is a significant peak response from baseline rise (A/S > 1) to baseline fall (A/S < 1) in IMF3 and IMF7. Additionally, in the IMF3 curve, mfs3 also exhibits a significant peak response. Similarly, on the opposite side of the baseline fall (A/S < 1) and baseline rise (A/S > 1), SB1, SB2, sb1, sb2, and sb3 have reasonable interpretations in the IMF3 curve. The IMF7 and R7 curves closely match the variations in the original GR log curve. The peak-valley positions of the MFS1 correspond precisely to the peaks and valleys in IMF7 and R7, completely coupled with the position of the MFS1. In summary, the use of the HHT for identifying sequence boundaries and maximum flooding surfaces results in more precise point locations compared to the range reflected by the DWT and CWT. For instance, in the case of MFS1, the peak value in IMF7 is significantly higher than the anomaly response shown at a = 128 in the CWT and at d11 in the DWT.

Figure 4 
               The results of applying HHT on the GR log of well YQ to interpret the stratigraphic cycles of the Jurassic Shaximiao Formation. MFS is the maximum flooding surface and SB is the sequence boundary.
Figure 4

The results of applying HHT on the GR log of well YQ to interpret the stratigraphic cycles of the Jurassic Shaximiao Formation. MFS is the maximum flooding surface and SB is the sequence boundary.

5 Discussion

The structure of the fluvial facies sequence is usually influenced by changes in base level and accommodation space, resulting in a dynamic process. High-resolution sequence stratigraphy, which utilizes the HT and HHT, can provide more detailed and comprehensive stratigraphic information. This serves as an essential method for finely dividing strata. The text uses GR logging curves as a starting point and applies the CWT, DWT, and HHT to process them. This process aims to remove extraneous interference signals from the GR logging curves, resulting in more precise identification of sequence boundaries and flooding surfaces.

As shown in Figure 5, the application of the CWT, DWT, and the HHT to the GR logging curves of well YQ reveals significant anomaly responses for most of the sequence boundaries and flooding surfaces. This approach minimizes the impact of subjective factors in sequence stratigraphy analysis. In the study area, the Suining Formation comprises “flood-overlake” strata dominated by meander river deposits and is located at the upper part of the top boundary (SB1) of the Shaximiao Formation [52]. The section below the lower boundary (SB2) corresponds to the Lianggaoshan Formation, which is primarily characterized by the development of delta-lake sedimentary systems [53,54]. A stable thick layer of sandstone deposited at SB1 and SB2 serves as a marker bed for stratigraphic division. The interval between 1,220 and 1,320 m can be identified as the SB1 interface if stratigraphic divisions are based solely on well logging curves. At depths of 2,340 and 2,370 m, potential indications of the SB2 interface exist. However, for a more precise delineation of SB1 and SB2, it is necessary to consider the anomalous response characteristics of low-frequency signals such as a = 32, d7, d11, IMF7, and R7. The time–frequency spectrogram does not provide significant indications due to the divergence of boundary energy.

Figure 5 
               The results of applying CWT, DWT, HHT on the GR log of well YQ to interpret the stratigraphic cycles of the Jurassic Shaximiao Formation. MFS is the maximum flooding surface and SB is the sequence boundary.
Figure 5

The results of applying CWT, DWT, HHT on the GR log of well YQ to interpret the stratigraphic cycles of the Jurassic Shaximiao Formation. MFS is the maximum flooding surface and SB is the sequence boundary.

The lower part of the Shaximiao Formation continues the sedimentary systems of the Lianggaoshan Formation period, which are predominantly characterized by shallow-water delta-fluvial facies. It experienced a brief period of shallow lake sedimentation in the later stages, resulting in a stable set of conchostracans shales [38]. These shales represent the MFS1 of this formation. Based on the well logging curves, MFS1 can be identified at 2,138–2,197 m depths. However, based on the analysis of the low-frequency signal IMF7, it is suggested that MFS1 should be positioned at 2,170 m.

The upper part of the study area, the Shaximiao Formation is characterized by typical deposits of fluvial deposition [55]. The thickness of the sand bodies increases significantly, allowing for the identification of four fourth-order sequences. Each of these sequences represents a trend where the curvature of meandering river channels initially increases and then decreases. On the one hand, sb1, sb2, and sb3 indicate positions where the curvature of river channels is relatively low, with a significant development of sandstone. On the other hand, mfs1, mfs2, mfs3, and mfs4 represent locations where the curvature of river channels is relatively high, with a fairly significant development of mudstone. These fourth-order sequence boundaries and their sub-level flooding surfaces are characterized by multiple similar points in well logging curves, which greatly increases uncertainty. However, analyzing the changing trends in the color of energy clusters in the time–frequency spectrogram, as well as high-frequency signals such as d5, d6, and IMF3, can provide a more precise identification of these boundaries.

The CWT maps the GR logging curve into the WT domain, allowing the creation of time–frequency spectrograms and wavelet reconstruction coefficient curves. This visual approach, supplemented by denoising curves, aids in identifying the features of the GR logging curve, facilitating the recognition of the cyclical periods of different sequences. The DWT and HHT are used to decompose the GR logging curve into fluctuation signals of different scales and frequencies, which reduces irrelevant noise. This process enhances the clarity of the amplitude in the original logging curve, enabling a more efficient revelation of the cyclical characteristics of sequence stratigraphy. After undergoing processing with the CWT and DWT, the GR logging curve can capture the cyclical signals of different stratigraphic sequences more clearly. Typically, interfaces can be constrained within specific depth ranges. However, the HHT provides more accurate identification of anomalous points by enhancing the curve characteristics of anomalous regions to determine the specific locations of interfaces.

However, when applying the WT and the HHT for high-resolution stratigraphy work, it is necessary to ensure that the input signal is sampled at a high rate. Therefore, before processing the GR curve, it is necessary to combine other methods to address this issue. Second, there may be instances where stratigraphic boundaries are eroded or older layers are overlaid by newer ones. In some well locations (Figure 6), the response characteristics of SB1 and SB2 in the CWT, DWT, and HHT are not distinct. In such cases, confirmation can only be obtained through other sources of data.

Figure 6 
               Comparative analysis of stratigraphic well Logs for Wells H, YQ, and ZT in the study area based on CWT, DWT and HHT of GR well log curves.
Figure 6

Comparative analysis of stratigraphic well Logs for Wells H, YQ, and ZT in the study area based on CWT, DWT and HHT of GR well log curves.

The CWT, DWT, and HHT can be used to identify sequence boundaries. However, well logging curves have inherent uncertainties, and anomalously low or high-value areas could be misinterpreted as sequence stratigraphic interfaces, adding to the overall uncertainty. Therefore, it is essential to first select logging curves that exhibit high sensitivity to sequence stratigraphic cycles. A combination of multiple logging curve processing methods should be employed to determine the stratigraphic division. The different advantages of each method should be leveraged by integrating and optimizing the results obtained from various methods. Various well logging curve processing methods can distinguish high-frequency signals corresponding to fourth-order cycles, while medium to low-frequency signals correspond to third-order cycles. This allows for the acquisition of high-resolution signals at different frequencies. High-resolution sequence stratigraphy can be achieved by combining and corroborating multiple methods.

It is evident that the stratigraphic division method described in this study is derived from the distinctive sedimentary characteristics of the Shaximiao Formation in southwestern Sichuan Basin. When dividing shallow-water delta-river facies, the primary considerations are the periodicity of sedimentary units and variations in accommodation space (A/S ratio) [14]. The utilization of the CWT, DWT, and HHT for the processing of the GR curve allows for the accurate identification of its periodic variations and anomalous responses, thus resulting in more precise stratigraphic division results. However, when it comes to stratigraphic division for different sedimentary systems, such as for gravity flow deposition, considerations shift towards the sedimentary body type, structural features, and the dynamic nature of the depositional environment [56]. The particle size and bedding characteristics of gravity flow deposition differ significantly from those of shallow-water delta-river facies sedimentation. The methods employed in this study may not be readily apparent in identifying these features, and the responses of different sequence boundaries may vary. Consequently, the aforementioned methods may not be suitable for the analysis of more complex sedimentary systems. A case-by-case approach is therefore required. The research methodology of this study primarily focuses on stratigraphic division work for shallow-water delta-river facies sedimentary systems. The method combining the CWT, DWT, and HHT has been applied to the study area of two additional wells, Well H and Well ZT. It was observed that similar sequence boundaries and flooding surfaces exhibit comparable response characteristics in these wells (Figure 5). In the southwestern Sichuan region, the Jurassic Shaximiao Formation is mainly supplied by sediment sources that are directed toward the northeast. In contrast, secondary sediment sources are present in the southwest direction. The stratigraphic thickness exhibits a distribution pattern of being thinner in the south and west and thicker in the north and east. The northeastern part of the region typically features thicker strata ranging between 1,000 and 1,600 m, while the southwestern part tends to have thinner strata ranging between 500 and 800 m. Sequence stratigraphy division is conducted by combining the three methods. The study utilized the CWT, DWT, and HHT to analyze sequence stratigraphy by comparing it with continuous well stratigraphic profiles (Figure 6). The results show that the overall trend of stratigraphic thickness variation is consistent with the distribution of stratigraphic thickness in the study area, further confirming the reliability of this comprehensive analysis approach.

Following the application of the CWT and DWT, it can be employed for the extraction of features from geological data, thereby enabling the identification of the frequency domain characteristics of underground structures. The HHT can be employed to extract both vibration modes and resonance characteristics from geological data. By leveraging these analytical techniques in combination, they can guide both imaging and signal recovery processes. The application of the method is not limited to stratigraphic division. There has been significant progress in the characterization of sand bodies [57], the interpretation of seismic profiles [58], the identification of fractures in well logging data [59], and the evaluation of reservoirs [29,60], among other areas. The application of these methods in geology is bound to increase with the advancement of technology, further enhancing the efficiency and accuracy of geological data processing and promoting rapid development in the fields of geology and geophysics.

6 Conclusion

This study presents a comprehensive study of the stratigraphic division of the Shaximiao Formation in southwestern Sichuan, which is subject to complex stratigraphic conditions. The research employs a range of high-resolution stratigraphy methods, including the CWT, DWT, and HHT, in order to gain a detailed understanding of the formation’s stratigraphic characteristics. The principal findings are as follows:

  1. The CWT, DWT, and HHT are more accurate in revealing the periodicity of stratigraphic changes and the anomalous response characteristics of well logging curves than traditional methods, overcoming the subjectivity and limitations of these traditional methods.

  2. Any individual method can be employed to guide stratigraphic work. However, by combining the CWT, DWT, and HHT, the precision of stratigraphic division can be further enhanced, providing favorable support for the stratigraphic division of the Shaximiao Formation in southwestern Sichuan where geological conditions are complex.

These findings are of significant importance not only for geological exploration and resource evaluation, but also for the methodology of geology. Nevertheless, despite the significant outcomes achieved, there are still some issues that require further investigation. These include the question of the universal applicability of the methods and the potential for further optimization to enhance efficiency. Future research could investigate the potential for combining other high-resolution geological methods with those proposed in this study, with a view to enhance the accuracy and reliability of stratigraphic division.

Acknowledgements

We would like to thank the constructive comments by the editors and the anonymous reviewers.

  1. Funding information: This work is funded by the Scientific Research and Technological Development Project of PetroChina: A comprehensive evaluation and optimization of the conditions for tight gas accumulation in the southwest Sichuan Province (JS2021-038).

  2. Author contributions: Zheng Li: conceptualization, data analysis, methodology, plotting, and writing – original draft. Tian Yang: supervision, validation, writing – review and editing. Jingchun Tian: supervision, validation, and writing – review and editing. Laixing Cai: supervision, validation, and writing – review and editing.

  3. Conflict of interest: All authors state no conflict of interest.

  4. Data availability statement: The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.

References

[1] Dai J, Ni Y, Hu G, Huang S, Liao F, Yu C, et al. Stable carbon and hydrogen isotopes of gases from the large tight gas fields in China. Sci China Earth Sci. 2014;57(1):88–103. 10.1007/s11430-013-4701-7.Search in Google Scholar

[2] Dai J, Ni Y, Wu X. Tight gas in China and its significance in exploration and exploitation. Pet Explor Dev. 2012;39(3):277–84. 10.1016/S1876-3804(12)60043-3.Search in Google Scholar

[3] Li Y, Chen S, Lu J, Xiao Z, He Q, Su K, et al. Main controlling factors of near-source and inter bedded-accumulation tight sandstone gas from coal-bearing strata: A case study from natural gas of Xujiahe Formation. Cent Sichuan Basin Nat Gas Geosci. 2019;30(6):798–808. 10.11764/j.issn.1672-1926.2019.02.001. (in Chinese with English abstract).Search in Google Scholar

[4] Zou C, Tao S, Han W, Zhao Z, Ma W, Li C, et al. Geological and geochemical characteristics and exploration prospect of coal‐derived tight sandstone gas in China: Case study of the Ordos, Sichuan, and Tarim Basins. Acta Geol Sin‐Engl. 2018;92(4):1609–26. 10.1111/1755-6724.13647.Search in Google Scholar

[5] Chen S, Yang Y, Qiu L, Wang X, Habilaxim E. Source of quartz cement and its impact on reservoir quality in Jurassic Shaximiao Formation in central Sichuan Basin, China. Mar Pet Geol. 2024;159:106543. 10.1016/j.marpetgeo.2023.106543.Search in Google Scholar

[6] Pan H, Jiang Y, Guo G, Yang C, Deng H, Zhu X, et al. Pore-throat structure characteristics and fluid mobility analysis of tight sandstone reservoirs in Shaximiao Formation, Central Sichuan. Geol J. 2023;58(11):4243–56. 10.1002/gj.4734.Search in Google Scholar

[7] Wu X, Huang S, Liao F, Li Z. Carbon isotopic characteristics of Jurassic alkane gases in the Sichuan Basin, China. Energy Explor Exploit. 2010;28(1):25–36. 10.1260/0144-5987.28.1.25.Search in Google Scholar

[8] Zhao Z, Tang D, Wang X, Chen SL. Discussion on main controlling factors of natural gas enrichment and high yield in tight sandstone gas reservoirs: Case study of Xujiahe Formation in Sichuan Basin. Nat Gas Geosci. 2019;30(7):963–72. 10.11764/j.issn.1672-1926.2019.03.014.Search in Google Scholar

[9] Amodu A, Oyetade OP, Fadiya SL, Fowora O. Sequence stratigraphic analysis and hydrocarbon prospectivity of AMO Field, deep offshore Niger Delta, Nigeria. Eng Geol. 2022;3(1):80–93. 10.1016/j.engeos.2021.11.003.Search in Google Scholar

[10] Allen J. Studies in fluviatile sedimentation: bars, bar-complexes and sandstone sheets (low-sinuosity braided streams) in the Brownstones (L. Devonian), Welsh Borders. Sediment Geol. 1983;33(4):237–93. 10.1016/0037-0738(83)90076-3.Search in Google Scholar

[11] Bohacs KM, Carroll AR, Neal JE, Mankiewicz PJ, Gierlowski-Kordesch EH, Kelts KR. Lake-basin type, source potential, and hydrocarbon character: an integrated sequence-stratigraphic-geochemical framework. Tulsa, Oklahoma, USA: Lake basins through space and time: AAPG Studies in Geology. Vol. 46, 2000. p. 3–34. 10.1306/St46706C1.Search in Google Scholar

[12] Catuneanu O. Principles of sequence stratigraphy. Newnes; 2022.10.1016/B978-0-444-53353-1.00006-7Search in Google Scholar

[13] Kumar P, Foufoula Georgiou E. Wavelet analysis for geophysical applications. Rev Geophys. 1997;35(4):385–412. 10.1029/97RG00427.Search in Google Scholar

[14] Xiao F, Wei T, Wang X, Guan X, Wu C, Hong H. Research on the sequence stratigraphy of the Shaximiao Formation in Chuan-zhong-Chuanxi area. Sichuan Basin Nat Gas Geosci. 2020;31(9):1216–24. 10.11764/j.issn.1672-1926.2020.04.023. (in Chinese with English abstract).Search in Google Scholar

[15] Zhang X, Deng H, Fu M, Ling C, Xu Z, Wu D. Controls of architecture under the constraints of a multilevel interface on physical property heterogeneities from the meandering river of the Shaximiao formation in Western Sichuan. Acta Sedimentol Sin. 2022;42(4):1384–400. 10.14027/j.issn.1000-0550.2022.105. (in Chinese with English abstract).Search in Google Scholar

[16] Song L, Liu S, Zeng Q, Zhou D, Tang D, Wang J. Genetic mechanism of relatively high-quality reservoirs of middle Jurassic Shaximiao formation tight sandstone in the transition zone between the central and western Sichuan Basin. J Jilin Univ (Earth Sci Ed). 2023;54(2):371–88. 10.13278/j.cnki.jjuese.20220208. (in Chinese with English abstract).Search in Google Scholar

[17] Zhu P, Ma T, Wang X, Li X, Dong Y, Yang W, et al. Wavelet transform coupled with Fischer plots for sequence stratigraphy: A case study in the Linxing area, Ordos Basin, China. J Pet Sci Eng. 2023;231:212306. 10.1016/j.geoen.2023.212306.Search in Google Scholar

[18] Li X, Fan Y, Deng S, Wang T. Automatic demarcation of sequence stratigraphy using the method of well logging multiscale data fusion. Pet Explor Dev. 2009;36(2):221–7. 1000-0747(2009)02-0221-07. (in Chinese with English abstract).Search in Google Scholar

[19] Pan S, Hsieh B, Lu M, Lin Z. Identification of stratigraphic formation interfaces using wavelet and Fourier transforms. Comput Geosci-UK. 2008;34(1):77–92. 10.1016/j.cageo.2007.01.002.Search in Google Scholar

[20] Weedon GP. Time-series analysis and cyclostratigraphy: examining stratigraphic records of environmental cycles. Cambridge, UK: Cambridge University Press; 2003.10.1017/CBO9780511535482Search in Google Scholar

[21] Chen S, Liu P, Tang D, Tao S, Zhang T. Identification of thin-layer coal texture using geophysical logging data: Investigation by Wavelet Transform and Linear Discrimination Analysis. Int J Coal Geol. 2021;239:103727. 10.1016/j.coal.2021.103727.Search in Google Scholar

[22] Wang R, Xie J, Ran A, Wang S, Wang J, Hu X, et al. Comparison of INPEFA technology and wavelet transform in sequence stratigraphic division of mixed reservoir: A case study of lower Es3 of KL oilfield in Laizhouwan Sag. J Pet Explor Prod Technol. 2022;12(12):3213–25. 10.1007/s13202-022-01523-z.Search in Google Scholar

[23] Zhang S, Yan J, Cai J, Zhu X, Hu Q, Wang M, et al. Fracture characteristics and logging identification of lacustrine shale in the Jiyang Depression, Bohai Bay Basin, Eastern China. Mar Pet Geol. 2021;132:105192. 10.1016/j.marpetgeo.2021.105192.Search in Google Scholar

[24] Wen Z, Zhu X, Jin M, Xu J. Demarcating stratigraphic sequence based on Hilbert-Huang transform of well-logging data. Well Logging Technol. 2007;31(6):533–6. 1004-1338(2007)06-0533-04. (in Chinese with English abstract).Search in Google Scholar

[25] Xu J, Wang G, Liu L, Zhang C, Xue Z. Application of Hilbert-Huang transform to study classification and correlation of isochrone substratum. Shiyou Diqiu Wuli Kantan (Oil Geophys Prospecting). 2010;45(3):423–30. 10.13810/j.cnki.issn.1000-7210.2010.03.02. (in Chinese with English abstract).Search in Google Scholar

[26] Rivera N, Ray S, Chan A, Jensen J, editors. Well-log feature extraction using wavelets and genetic algorithms. Proc Am Assoc Petr Geol ann meeting. Houston, Texas: 2002.Search in Google Scholar

[27] Climent LPF, Associes RR, Lescar F. Sequence stratigraphy applied to log interpretation: Improving methodology by means of signal processing techniques and outcrop calibration. Barcelona, Spain: AAPG. 2003.Search in Google Scholar

[28] Karacan CÖ, Olea RA. Inference of strata separation and gas emission paths in longwall overburden using continuous wavelet transform of well logs and geostatistical simulation. J Appl Geophys. 2014;105:147–58. 10.1016/j.jappgeo.2014.03.019.Search in Google Scholar

[29] Perez-Muñoz T, Velasco-Hernandez J, Hernandez-Martinez E. Wavelet transform analysis for lithological characteristics identification in siliciclastic oil fields. J Appl Geophys. 2013;98:298–308. 10.1016/j.jappgeo.2013.09.010.Search in Google Scholar

[30] Schultz SK, MacEachern JA, Catuneanu O, Dashtgard SE, Díaz N. High-resolution sequence stratigraphic framework for the late Albian Viking Formation in central Alberta. Mar Pet Geol. 2022;139:105627. 10.1016/j.marpetgeo.2022.105627.Search in Google Scholar

[31] Zhang Q, Zhang F, Liu J, Wang X, Chen Q, Zhao L, et al. A method for identifying the thin layer using the wavelet transform of density logging data. J Pet Sci Eng. 2018;160:433–41. 10.1016/j.petrol.2017.10.048.Search in Google Scholar

[32] Huang NE, Shen Z, Long SR, Wu MC, Shih HH, Zheng Q, et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc R Soc Lond Ser A: Math Phys Eng Sci. 1998;454(1971):903–95. 10.1098/rspa.1998.0193.Search in Google Scholar

[33] Huang NE. Computer implemented empirical mode decomposition method, apparatus, and article of manufacture for two-dimensional signals; 2001.Search in Google Scholar

[34] Huang NE, Wu MC, Long SR, Shen SS, Qu W, Gloersen P, et al. A confidence limit for the empirical mode decomposition and Hilbert spectral analysis. Proc R Soc Lond Ser A: Math Phys Eng Sci. 2003;459(2037):2317–45. org/10.1098/rspa.2003.1123.Search in Google Scholar

[35] Amoura S, Gaci S, Barbosa S, Farfour M, Bounif MA. Investigation of lithological heterogeneities from velocity logs using EMD-Hölder technique combined with multifractal analysis and unsupervised statistical methods. J Pet Sci Eng. 2022;208:109588. 10.1016/j.petrol.2021.109588.Search in Google Scholar

[36] Liang Y, Gu H, Yao Z. The application of improved Hilbert-Huang transform in reservoir prediction. Geophys Prospect Pet. 2016;55(4):606–15. 10.3969/j.issn.1000-1441.2016.04.016.Search in Google Scholar

[37] Roveri N, Carcaterra A. Damage detection in structures under traveling loads by Hilbert–Huang transform. Mech Syst Signal Process. 2012;28:128–44. 10.1016/j.ymssp.2011.06.018.Search in Google Scholar

[38] He Q, Yang T, Cai L, Ren Q, Dai J. Porosity and permeability correction of laumontite-rich reservoirs in the first member of the Shaximiao Formation in the Central Sichuan Basin. Mar Pet Geol. 2023;150:106115. 10.1016/j.marpetgeo.2023.106115.Search in Google Scholar

[39] Zhang X, Fu M, Deng H, Li Z, Zhao S, Gluyas JG, et al. The differential diagenesis controls on the physical properties of lithofacies in sandstone reservoirs from the Jurassic Shaximiao Formation, western Sichuan depression, China. J Pet Sci Eng. 2020;193:107413. 10.1016/j.petrol.2020.107413.Search in Google Scholar

[40] Tan X, Jiang W, Luo L, Gluyas J, Song L, Liu J, et al. Variations of sedimentary environment under cyclical aridification and impacts on eodiagenesis of tight sandstones from the late Middle Jurassic Shaximiao Formation in Central Sichuan Basin. Mar Pet Geol. 2024;161:106699. 10.1016/j.marpetgeo.2024.106699.Search in Google Scholar

[41] Zhang X, Deng H, Li T, Xu Z, Fu M, Ling C, et al. The conditions and modelling of hydrocarbon accumulation in tight sandstone reservoirs: A case study from the Jurassic Shaximiao formation of western Sichuan Basin, China. J Pet Sci Eng. 2023;225:211702. 10.1016/j.geoen.2023.211702.Search in Google Scholar

[42] Kadkhodaie A, Rezaee R. Intelligent sequence stratigraphy through a wavelet-based decomposition of well log data. J Nat Gas Sci Eng. 2017;40:38–50. 10.1016/j.jngse.2017.02.010.Search in Google Scholar

[43] Nazeer A, Abbasi SA, Solangi SH. Sedimentary facies interpretation of Gamma Ray (GR) log as basic well logs in Central and Lower Indus Basin of Pakistan. Geod Geodyn. 2016;7(6):432–43. 10.1016/j.geog.2016.06.006.Search in Google Scholar

[44] Antonini M, Barlaud M, Mathieu P, Daubechies I. Image coding using wavelet transform. IEEE Trans Image Process. 1992;1:20–5.10.1109/83.136597Search in Google Scholar PubMed

[45] Bolton EW, Maasch KA, Lilly JM. A wavelet analysis of Plio‐Pleistocene climate indicators: A new view of periodicity evolution. Geophys Res Lett. 1995;22(20):2753–6. 10.1029/95GL02799.Search in Google Scholar

[46] Daubechies I. Orthonormal bases of compactly supported wavelets. Commun Pure Appl Math. 1988;41(7):909–96. org/10.1002/cpa.3160410705.Search in Google Scholar

[47] Huang NE, Shen Z, Long SR, Wu MC, Shih HH, Zheng Q, et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc R Soc London Ser A Math Phys Eng Sci. 1998;454(1971):903–95. 10.1098/rspa.1998.0193.Search in Google Scholar

[48] Huang NE. Hilbert-Huang transform and its applications. World Scientific; 2014.10.1142/8804Search in Google Scholar

[49] Xu J, Wang G, Liu L, Zhang C, Xue Z. Application of Hilbert-Huang transform to study classification and correlation of isochrone substratum. Shiyou Diqiu Wuli Kantan (Oil Geophys Prospect). 2010;45(3):423–30. 10.13810/j.cnki.issn.1000-7210.2010.03.024. (in Chinese with English abstract).Search in Google Scholar

[50] Zhang J, Wu S, Hu G, Fan T, Yu B, Lin P, et al. Sea-level control on the submarine fan architecture in a deepwater sequence of the Niger Delta Basin. Mar Pet Geol. 2018;94:179–97. 10.1016/j.marpetgeo.2018.04.002.Search in Google Scholar

[51] Deng H, Wu H, Wang N, Timothy AC. Division of fluvial sequence stratigraphy: An example from the Lower Cretaceous Fuyu oil-bearing layer, the Songliao Basin. Oil Gas Geol. 2007;28(5):621–7. 0253-9985(2007)05-0621-07.Search in Google Scholar

[52] Chen X, Ji Y, Yang K, Liu J, Huang F, Zhao C, et al. Flood overlake sedimentary characteristics of the Suining formation (Lower Jurassic) in Western Si chuan depression. Acta Sedimentol Sin. 2014;32(5):912–20. 10.14027/j.cnki.cjxb.2014.05.035. (in Chinese with English abstract).Search in Google Scholar

[53] Dawei C, Zhijie Z, Haitao H, Shaomin Z, Chunyu Q, Xuanjun Y, et al. Sequence structure, sedimentary evolution and their controlling factors of the Jurassic Lianggaoshan Formation in the East Sichuan Basin, SW China. Pet Explor Dev. 2023;50(2):262–72. 10.1016/S1876-3804(23)60388-X.Search in Google Scholar

[54] Yi J, Zhang S, Cai L, Chen S, Luo X, Yu J, et al. Strata and sedimentary filling characteristics of the Lower Jurassic Lianggaoshan Formation and its hydrocarbon exploration in eastern Sichuan Basin. J Jilin Univ (Earth Sci Ed). 2022;52(3):795–815. 10.13278/j.cnki.Jjuese.20210429. (in Chinese with English abstract).Search in Google Scholar

[55] Yu W, Yang T, Cai L, Li X, He Q. Paleoenvironment and paleoclimate evolution during the depositional period of the Middle Jurassic Shaximiao Formation in the central Sichuan basin: A case tudy of well Yongqian 1. Acta Geol Sin. 2023. 10.19762/j.cnki.dizhixuebao.2023283. (in Chinese with English abstract).Search in Google Scholar

[56] Yang T, Cao Y, Wang Y, Cai L, Liu H, Jin J. Sedimentary characteristics and depositional model of hyperpycnites in the gentle slope of a lacustrine rift basin: A case study from the third member of the Eocene Shahejie Formation, Bonan Sag, Bohai Bay Basin, Eastern China. Basin Res. 2023;35(4):1590–618. 10.1111/bre.12766.Search in Google Scholar

[57] Xu Z, Zhang B, Li F, Cao G, Liu Y. Well-log decomposition using variational mode decomposition in assisting the sequence stratigraphy analysis of a conglomerate reservoir. GeoPhysics. 2018;83(4):B221–8. 10.1190/GEO2017-0817.1.Search in Google Scholar

[58] Cooper M. Structural style and hydrocarbon prospectivity in fold and thrust belts: A global review. London, UK: Geological Society of London; 2007.10.1144/GSL.SP.2007.272.01.23Search in Google Scholar

[59] Tokhmechi B, Memarian H, Rasouli V, Noubari HA, Moshiri B. Fracture detection from water saturation log data using a Fourier–wavelet approach. J Pet Sci Eng. 2009;69(1-2):129–38. 10.1016/j.petrol.2009.08.005.Search in Google Scholar

[60] Chandrasekhar E, Eswara Rao V. Wavelet analysis of geophysical well-log data of Bombay Offshore Basin, India. Math Geosci. 2012;44(8):901–28. 10.1007/s11004-012-9423-4.Search in Google Scholar

Received: 2024-03-03
Revised: 2024-04-23
Accepted: 2024-05-01
Published Online: 2024-09-03

© 2024 the author(s), published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Regular Articles
  2. Theoretical magnetotelluric response of stratiform earth consisting of alternative homogeneous and transitional layers
  3. The research of common drought indexes for the application to the drought monitoring in the region of Jin Sha river
  4. Evolutionary game analysis of government, businesses, and consumers in high-standard farmland low-carbon construction
  5. On the use of low-frequency passive seismic as a direct hydrocarbon indicator: A case study at Banyubang oil field, Indonesia
  6. Water transportation planning in connection with extreme weather conditions; case study – Port of Novi Sad, Serbia
  7. Zircon U–Pb ages of the Paleozoic volcaniclastic strata in the Junggar Basin, NW China
  8. Monitoring of mangrove forests vegetation based on optical versus microwave data: A case study western coast of Saudi Arabia
  9. Microfacies analysis of marine shale: A case study of the shales of the Wufeng–Longmaxi formation in the western Chongqing, Sichuan Basin, China
  10. Multisource remote sensing image fusion processing in plateau seismic region feature information extraction and application analysis – An example of the Menyuan Ms6.9 earthquake on January 8, 2022
  11. Identification of magnetic mineralogy and paleo-flow direction of the Miocene-quaternary volcanic products in the north of Lake Van, Eastern Turkey
  12. Impact of fully rotating steel casing bored pile on adjacent tunnels
  13. Adolescents’ consumption intentions toward leisure tourism in high-risk leisure environments in riverine areas
  14. Petrogenesis of Jurassic granitic rocks in South China Block: Implications for events related to subduction of Paleo-Pacific plate
  15. Differences in urban daytime and night block vitality based on mobile phone signaling data: A case study of Kunming’s urban district
  16. Random forest and artificial neural network-based tsunami forests classification using data fusion of Sentinel-2 and Airbus Vision-1 satellites: A case study of Garhi Chandan, Pakistan
  17. Integrated geophysical approach for detection and size-geometry characterization of a multiscale karst system in carbonate units, semiarid Brazil
  18. Spatial and temporal changes in ecosystem services value and analysis of driving factors in the Yangtze River Delta Region
  19. Deep fault sliding rates for Ka-Ping block of Xinjiang based on repeating earthquakes
  20. Improved deep learning segmentation of outdoor point clouds with different sampling strategies and using intensities
  21. Platform margin belt structure and sedimentation characteristics of Changxing Formation reefs on both sides of the Kaijiang-Liangping trough, eastern Sichuan Basin, China
  22. Enhancing attapulgite and cement-modified loess for effective landfill lining: A study on seepage prevention and Cu/Pb ion adsorption
  23. Flood risk assessment, a case study in an arid environment of Southeast Morocco
  24. Lower limits of physical properties and classification evaluation criteria of the tight reservoir in the Ahe Formation in the Dibei Area of the Kuqa depression
  25. Evaluation of Viaducts’ contribution to road network accessibility in the Yunnan–Guizhou area based on the node deletion method
  26. Permian tectonic switch of the southern Central Asian Orogenic Belt: Constraints from magmatism in the southern Alxa region, NW China
  27. Element geochemical differences in lower Cambrian black shales with hydrothermal sedimentation in the Yangtze block, South China
  28. Three-dimensional finite-memory quasi-Newton inversion of the magnetotelluric based on unstructured grids
  29. Obliquity-paced summer monsoon from the Shilou red clay section on the eastern Chinese Loess Plateau
  30. Classification and logging identification of reservoir space near the upper Ordovician pinch-out line in Tahe Oilfield
  31. Ultra-deep channel sand body target recognition method based on improved deep learning under UAV cluster
  32. New formula to determine flyrock distance on sedimentary rocks with low strength
  33. Assessing the ecological security of tourism in Northeast China
  34. Effective reservoir identification and sweet spot prediction in Chang 8 Member tight oil reservoirs in Huanjiang area, Ordos Basin
  35. Detecting heterogeneity of spatial accessibility to sports facilities for adolescents at fine scale: A case study in Changsha, China
  36. Effects of freeze–thaw cycles on soil nutrients by soft rock and sand remodeling
  37. Vibration prediction with a method based on the absorption property of blast-induced seismic waves: A case study
  38. A new look at the geodynamic development of the Ediacaran–early Cambrian forearc basalts of the Tannuola-Khamsara Island Arc (Central Asia, Russia): Conclusions from geological, geochemical, and Nd-isotope data
  39. Spatio-temporal analysis of the driving factors of urban land use expansion in China: A study of the Yangtze River Delta region
  40. Selection of Euler deconvolution solutions using the enhanced horizontal gradient and stable vertical differentiation
  41. Phase change of the Ordovician hydrocarbon in the Tarim Basin: A case study from the Halahatang–Shunbei area
  42. Using interpretative structure model and analytical network process for optimum site selection of airport locations in Delta Egypt
  43. Geochemistry of magnetite from Fe-skarn deposits along the central Loei Fold Belt, Thailand
  44. Functional typology of settlements in the Srem region, Serbia
  45. Hunger Games Search for the elucidation of gravity anomalies with application to geothermal energy investigations and volcanic activity studies
  46. Addressing incomplete tile phenomena in image tiling: Introducing the grid six-intersection model
  47. Evaluation and control model for resilience of water resource building system based on fuzzy comprehensive evaluation method and its application
  48. MIF and AHP methods for delineation of groundwater potential zones using remote sensing and GIS techniques in Tirunelveli, Tenkasi District, India
  49. New database for the estimation of dynamic coefficient of friction of snow
  50. Measuring urban growth dynamics: A study in Hue city, Vietnam
  51. Comparative models of support-vector machine, multilayer perceptron, and decision tree ‎predication approaches for landslide ‎susceptibility analysis
  52. Experimental study on the influence of clay content on the shear strength of silty soil and mechanism analysis
  53. Geosite assessment as a contribution to the sustainable development of Babušnica, Serbia
  54. Using fuzzy analytical hierarchy process for road transportation services management based on remote sensing and GIS technology
  55. Accumulation mechanism of multi-type unconventional oil and gas reservoirs in Northern China: Taking Hari Sag of the Yin’e Basin as an example
  56. TOC prediction of source rocks based on the convolutional neural network and logging curves – A case study of Pinghu Formation in Xihu Sag
  57. A method for fast detection of wind farms from remote sensing images using deep learning and geospatial analysis
  58. Spatial distribution and driving factors of karst rocky desertification in Southwest China based on GIS and geodetector
  59. Physicochemical and mineralogical composition studies of clays from Share and Tshonga areas, Northern Bida Basin, Nigeria: Implications for Geophagia
  60. Geochemical sedimentary records of eutrophication and environmental change in Chaohu Lake, East China
  61. Research progress of freeze–thaw rock using bibliometric analysis
  62. Mixed irrigation affects the composition and diversity of the soil bacterial community
  63. Examining the swelling potential of cohesive soils with high plasticity according to their index properties using GIS
  64. Geological genesis and identification of high-porosity and low-permeability sandstones in the Cretaceous Bashkirchik Formation, northern Tarim Basin
  65. Usability of PPGIS tools exemplified by geodiscussion – a tool for public participation in shaping public space
  66. Efficient development technology of Upper Paleozoic Lower Shihezi tight sandstone gas reservoir in northeastern Ordos Basin
  67. Assessment of soil resources of agricultural landscapes in Turkestan region of the Republic of Kazakhstan based on agrochemical indexes
  68. Evaluating the impact of DEM interpolation algorithms on relief index for soil resource management
  69. Petrogenetic relationship between plutonic and subvolcanic rocks in the Jurassic Shuikoushan complex, South China
  70. A novel workflow for shale lithology identification – A case study in the Gulong Depression, Songliao Basin, China
  71. Characteristics and main controlling factors of dolomite reservoirs in Fei-3 Member of Feixianguan Formation of Lower Triassic, Puguang area
  72. Impact of high-speed railway network on county-level accessibility and economic linkage in Jiangxi Province, China: A spatio-temporal data analysis
  73. Estimation model of wild fractional vegetation cover based on RGB vegetation index and its application
  74. Lithofacies, petrography, and geochemistry of the Lamphun oceanic plate stratigraphy: As a record of the subduction history of Paleo-Tethys in Chiang Mai-Chiang Rai Suture Zone of Thailand
  75. Structural features and tectonic activity of the Weihe Fault, central China
  76. Application of the wavelet transform and Hilbert–Huang transform in stratigraphic sequence division of Jurassic Shaximiao Formation in Southwest Sichuan Basin
  77. Structural detachment influences the shale gas preservation in the Wufeng-Longmaxi Formation, Northern Guizhou Province
  78. Distribution law of Chang 7 Member tight oil in the western Ordos Basin based on geological, logging and numerical simulation techniques
  79. Evaluation of alteration in the geothermal province west of Cappadocia, Türkiye: Mineralogical, petrographical, geochemical, and remote sensing data
  80. Numerical modeling of site response at large strains with simplified nonlinear models: Application to Lotung seismic array
  81. Quantitative characterization of granite failure intensity under dynamic disturbance from energy standpoint
  82. Characteristics of debris flow dynamics and prediction of the hazardous area in Bangou Village, Yanqing District, Beijing, China
  83. Rockfall mapping and susceptibility evaluation based on UAV high-resolution imagery and support vector machine method
  84. Statistical comparison analysis of different real-time kinematic methods for the development of photogrammetric products: CORS-RTK, CORS-RTK + PPK, RTK-DRTK2, and RTK + DRTK2 + GCP
  85. Hydrogeological mapping of fracture networks using earth observation data to improve rainfall–runoff modeling in arid mountains, Saudi Arabia
  86. Petrography and geochemistry of pegmatite and leucogranite of Ntega-Marangara area, Burundi, in relation to rare metal mineralisation
  87. Prediction of formation fracture pressure based on reinforcement learning and XGBoost
  88. Hazard zonation for potential earthquake-induced landslide in the eastern East Kunlun fault zone
  89. Monitoring water infiltration in multiple layers of sandstone coal mining model with cracks using ERT
  90. Study of the patterns of ice lake variation and the factors influencing these changes in the western Nyingchi area
  91. Productive conservation at the landslide prone area under the threat of rapid land cover changes
  92. Sedimentary processes and patterns in deposits corresponding to freshwater lake-facies of hyperpycnal flow – An experimental study based on flume depositional simulations
  93. Study on time-dependent injectability evaluation of mudstone considering the self-healing effect
  94. Detection of objects with diverse geometric shapes in GPR images using deep-learning methods
  95. Behavior of trace metals in sedimentary cores from marine and lacustrine environments in Algeria
  96. Spatiotemporal variation pattern and spatial coupling relationship between NDVI and LST in Mu Us Sandy Land
  97. Formation mechanism and oil-bearing properties of gravity flow sand body of Chang 63 sub-member of Yanchang Formation in Huaqing area, Ordos Basin
  98. Diagenesis of marine-continental transitional shale from the Upper Permian Longtan Formation in southern Sichuan Basin, China
  99. Vertical high-velocity structures and seismic activity in western Shandong Rise, China: Case study inspired by double-difference seismic tomography
  100. Spatial coupling relationship between metamorphic core complex and gold deposits: Constraints from geophysical electromagnetics
  101. Disparities in the geospatial allocation of public facilities from the perspective of living circles
  102. Research on spatial correlation structure of war heritage based on field theory. A case study of Jinzhai County, China
  103. Formation mechanisms of Qiaoba-Zhongdu Danxia landforms in southwestern Sichuan Province, China
  104. Magnetic data interpretation: Implication for structure and hydrocarbon potentiality at Delta Wadi Diit, Southeastern Egypt
  105. Deeply buried clastic rock diagenesis evolution mechanism of Dongdaohaizi sag in the center of Junggar fault basin, Northwest China
  106. Application of LS-RAPID to simulate the motion of two contrasting landslides triggered by earthquakes
  107. The new insight of tectonic setting in Sunda–Banda transition zone using tomography seismic. Case study: 7.1 M deep earthquake 29 August 2023
  108. The critical role of c and φ in ensuring stability: A study on rockfill dams
  109. Evidence of late quaternary activity of the Weining-Shuicheng Fault in Guizhou, China
  110. Extreme hydroclimatic events and response of vegetation in the eastern QTP since 10 ka
  111. Spatial–temporal effect of sea–land gradient on landscape pattern and ecological risk in the coastal zone: A case study of Dalian City
  112. Study on the influence mechanism of land use on carbon storage under multiple scenarios: A case study of Wenzhou
  113. A new method for identifying reservoir fluid properties based on well logging data: A case study from PL block of Bohai Bay Basin, North China
  114. Comparison between thermal models across the Middle Magdalena Valley, Eastern Cordillera, and Eastern Llanos basins in Colombia
  115. Mineralogical and elemental analysis of Kazakh coals from three mines: Preliminary insights from mode of occurrence to environmental impacts
  116. Chlorite-induced porosity evolution in multi-source tight sandstone reservoirs: A case study of the Shaximiao Formation in western Sichuan Basin
  117. Predicting stability factors for rotational failures in earth slopes and embankments using artificial intelligence techniques
  118. Origin of Late Cretaceous A-type granitoids in South China: Response to the rollback and retreat of the Paleo-Pacific plate
  119. Modification of dolomitization on reservoir spaces in reef–shoal complex: A case study of Permian Changxing Formation, Sichuan Basin, SW China
  120. Geological characteristics of the Daduhe gold belt, western Sichuan, China: Implications for exploration
  121. Rock physics model for deep coal-bed methane reservoir based on equivalent medium theory: A case study of Carboniferous-Permian in Eastern Ordos Basin
  122. Enhancing the total-field magnetic anomaly using the normalized source strength
  123. Shear wave velocity profiling of Riyadh City, Saudi Arabia, utilizing the multi-channel analysis of surface waves method
  124. Effect of coal facies on pore structure heterogeneity of coal measures: Quantitative characterization and comparative study
  125. Inversion method of organic matter content of different types of soils in black soil area based on hyperspectral indices
  126. Detection of seepage zones in artificial levees: A case study at the Körös River, Hungary
  127. Tight sandstone fluid detection technology based on multi-wave seismic data
  128. Characteristics and control techniques of soft rock tunnel lining cracks in high geo-stress environments: Case study of Wushaoling tunnel group
  129. Influence of pore structure characteristics on the Permian Shan-1 reservoir in Longdong, Southwest Ordos Basin, China
  130. Study on sedimentary model of Shanxi Formation – Lower Shihezi Formation in Da 17 well area of Daniudi gas field, Ordos Basin
  131. Multi-scenario territorial spatial simulation and dynamic changes: A case study of Jilin Province in China from 1985 to 2030
  132. Review Articles
  133. Major ascidian species with negative impacts on bivalve aquaculture: Current knowledge and future research aims
  134. Prediction and assessment of meteorological drought in southwest China using long short-term memory model
  135. Communication
  136. Essential questions in earth and geosciences according to large language models
  137. Erratum
  138. Erratum to “Random forest and artificial neural network-based tsunami forests classification using data fusion of Sentinel-2 and Airbus Vision-1 satellites: A case study of Garhi Chandan, Pakistan”
  139. Special Issue: Natural Resources and Environmental Risks: Towards a Sustainable Future - Part I
  140. Spatial-temporal and trend analysis of traffic accidents in AP Vojvodina (North Serbia)
  141. Exploring environmental awareness, knowledge, and safety: A comparative study among students in Montenegro and North Macedonia
  142. Determinants influencing tourists’ willingness to visit Türkiye – Impact of earthquake hazards on Serbian visitors’ preferences
  143. Application of remote sensing in monitoring land degradation: A case study of Stanari municipality (Bosnia and Herzegovina)
  144. Optimizing agricultural land use: A GIS-based assessment of suitability in the Sana River Basin, Bosnia and Herzegovina
  145. Assessing risk-prone areas in the Kratovska Reka catchment (North Macedonia) by integrating advanced geospatial analytics and flash flood potential index
  146. Analysis of the intensity of erosive processes and state of vegetation cover in the zone of influence of the Kolubara Mining Basin
  147. GIS-based spatial modeling of landslide susceptibility using BWM-LSI: A case study – city of Smederevo (Serbia)
  148. Geospatial modeling of wildfire susceptibility on a national scale in Montenegro: A comparative evaluation of F-AHP and FR methodologies
  149. Geosite assessment as the first step for the development of canyoning activities in North Montenegro
  150. Urban geoheritage and degradation risk assessment of the Sokograd fortress (Sokobanja, Eastern Serbia)
  151. Multi-hazard modeling of erosion and landslide susceptibility at the national scale in the example of North Macedonia
  152. Understanding seismic hazard resilience in Montenegro: A qualitative analysis of community preparedness and response capabilities
  153. Forest soil CO2 emission in Quercus robur level II monitoring site
  154. Characterization of glomalin proteins in soil: A potential indicator of erosion intensity
  155. Power of Terroir: Case study of Grašac at the Fruška Gora wine region (North Serbia)
  156. Special Issue: Geospatial and Environmental Dynamics - Part I
  157. Qualitative insights into cultural heritage protection in Serbia: Addressing legal and institutional gaps for disaster risk resilience
Downloaded on 8.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/geo-2022-0657/html
Scroll to top button