Analysis of Seismic Response Characteristics of Coal-Bearing Strata in Pinghu Formation of Xihu Depression ()
1. Introduction
Seismic forward modeling is an essential part of seismic reservoir prediction research. It helps to better understand the geophysical characteristics and geological significance of complex propagation media in strata. Fu et al. (2010) used wave equation forward numerical simulation to study the seismic reflection structure characteristics of bioherms. Liu (2014) conducted forward modeling to study the influence of coal seam thickness on reflection wave energy and frequency, laying a theoretical foundation for interpreting coal thickness changes using seismic attribute technology. Li et al. (2018) established models based on drilling data for different types of sand-mudstone combinations and conducted forward modeling to study the internal seismic facies structure characteristics of deep-water fans, matching them with actual seismic facies identification. Niu et al. (2022) used pre-stack reflection low-value anomalies of sand bodies and post-stack reflection low-value anomalies of coal-mud interbeds to identify lithologic traps in coal-bearing strata in the Pingbei area of the Xihu Depression. Zhang et al. (2024) used seismic forward modeling to study the relationship between sand thickness and seismic attributes in tight sandstone reservoirs. Wang et al. (2017) used seismic forward and inverse modeling combined with petrophysical analysis to study tight reservoirs in the Loess Plateau region. Mirkamali et al. (2023) used seismic forward modeling combined with seismic petrophysics, geological model construction, and 2D finite-difference elastic modeling to investigate and interpret thin layers in carbonate reservoirs in southwestern Iran. Johansen conducted seismic forward modeling of the Kvalhovden outcrop in Spitsbergen, Norway, exploring the resolution gap between outcrop, reservoir, and seismic scale observations (Johansen et al., 2023). Zhang et al. (2020) conducted detailed seismic forward modeling analysis of coastal oilfields, exploring the relationship between seismic wave penetration and branching.
In the study of seismic reflection configuration characteristics, there is limited research on the influence of coal seams on the seismic reflection configuration of sand bodies. The coal-bearing strata of the Pinghu Formation in the Pingbei area of the Xihu Depression are the main source rocks in the Xihu Depression. Due to the widespread distribution of coal-bearing strata, the large number of coal seams, thin single-layer thickness, and the significant influence of coal seams on the identification of sedimentary microfacies in seismic profiles, there are many controversies regarding the sedimentary facies of the Pinghu Formation. Xiong et al. (2007) believe that the spatial distribution of coal seams in the Pinghu Formation is characterized by being numerous, scattered, thin, widespread, and uniform. Zhang & Xu (2012) conducted detailed research on the sequence stratigraphy and sedimentary facies of the Pinghu Formation in the Xihu Depression, concluding that one of the lithological characteristics of the Pinghu Formation is the presence of numerous thin and scattered coal seams in both mudstone and sandstone. Zhou et al. (2016) believe that the tidal flat and tide-controlled delta environments in the Pinghu Formation of the Xihu Depression have better coal accumulation effects than the braided river delta, and the main factors influencing the matching relationship between the regional sedimentary system and the coal-rich zone are the coal accumulation environment and mode. Shen (2018) believes that the main coal accumulation units in the Pinghu Formation are the supratidal and intertidal zones, with coal seams characterized by thin single-layer thickness and large cumulative thickness. He also points out that the uneven distribution of coal seams in the Pinghu Formation is due to changes in accommodation space during different tectonic evolution stages. Wang et al. (2021) believes that the coal seams in the Pinghu Formation of the Xihu Depression developed in tide-controlled delta and river-controlled delta peat swamp environments, characterized by multiple layers and thin thickness.
Through previous research, it has been found that the coal accumulation environment of the Paleogene Pinghu Formation in the Xihu Depression is mainly tidal flats and tide-controlled delta plain peat swamps (Zhou et al., 2016). Based on well logging data and Tesseral software, this paper constructs theoretical geological models for the coal-bearing strata of the Paleogene Pinghu Formation in the Xihu Depression. Through seismic forward modeling, the influence of coal seam thickness and spatial distribution on the seismic reflection characteristics of sand bodies is analyzed. Variables such as the distance between sand and coal layers as well as coal seam thickness are integrated into the model. The results are then compared and analyzed with actual seismic profiles and lithological data. Ultimately, conclusions are drawn, providing a basis for the study of regional sedimentary systems and the configuration of typical seismic profiles.
2. Regional Geological Background and Coal Accumulation
Environment Characteristics
2.1. Regional Geological Background
The Xihu Depression is one of the most oil- and gas-rich depressions in China’s offshore areas with significant exploration potential. It is located in the northeastern part of the East China Sea Basin, belonging to the eastern depression belt of the East China Sea Shelf Basin, with a water depth of about 100 m and an area of 4.27 × 104 km2 (Xu et al., 2020). To the east is the Diaoyu Islands Fold Belt, and to the west is the Central Uplift (Figure 1) (Zhang et al., 2011). The Xihu Depression began to deposit in the late Cretaceous, with obvious east-fault and west-overlap characteristics, forming a half-graben basin. It includes three main structural units from east to west: the western slope zone, the central sag-reversal zone, and the eastern fault-step zone. The strata of the Xihu Depression, from bottom to top, include the Paleocene, Eocene Pinghu Formation, Oligocene Huagang Formation, Miocene Longjing Formation, Yuquan Formation, and Liulang Formation, Pliocene Santan Formation, and Quaternary Donghai Group. The Pinghu Formation is the main target layer for oil and gas exploration (Zhang et al., 2012).
![]()
Figure 1. Location of research area and brief table of stratigraphic sequence.
2.2. Coal Accumulation Environment Characteristics
The main coal-bearing strata in the Xihu Depression are distributed in the Pinghu Formation. The main coal accumulation environments of the Pinghu Formation are:
1) Braided river delta plain peat swamp: The braided river delta plain is a terrestrial deposit with well-developed distributary channels, which strongly erode muddy deposits, making peat accumulation difficult and significantly affecting coal accumulation. The lower delta plain is an underwater deposit, where underwater distributary bays are similar to tide-controlled delta distributary bays, easily developing swamp environments conducive to coal seam formation (Figure 2(a)).
2) Tidal flat peat swamp: Tidal flat peat swamps develop in the supratidal and intertidal zones. Due to the periodic nature of tides, parts of the intertidal and supratidal zones periodically emerge above water, and under suitable climatic conditions, these areas easily develop peat swamps and form coal seams. Due to the periodic influence, the coal seams formed in this environment are characterized by thin single-layer thickness and multiple layers (Li et al., 2013) (Figure 2(b)).
3) Tide-controlled delta plain peat swamp: Due to coastline advancement, tide-controlled delta plain distributary bays are easily filled with sediments and become shallower. Under these conditions, plants and peat accumulation easily develop into swamp environments (Guo et al., 1998). Tide-controlled delta distributary bays, after being filled with sediments and becoming shallower, are suitable for peat swamp development and coal seam formation. Abandoned tidal sandbars and abandoned distributary channels in the tide-controlled delta plain are also important coal accumulation environments (Figure 2(c)).
In general, the coal accumulation environment characteristics of the Pinghu Formation in the Xihu Depression are: stable structure, suitable water depth, stable water body, stable climate, conducive to plant growth and peat accumulation, and easy development of peat swamps.
Figure 2. Contour map of coal series thickness in the research area.
3. Basic Principles of Acoustic Forward Modeling
Seismic forward modeling, or seismic wavefield simulation, simulates seismic wave signals in real geological structures based on the geophysical properties (velocity, density, etc.) of different lithologies in well logging data, obtaining seismic records close to real ones. This forward model is based on the Tesseral-2D velocity-depth model, which effectively simulates the 2D wave propagation effects of seismic energy in real situations, without considering the hardness of solid media, representing an ideal fluid medium situation. By summarizing velocity and density parameters from well logging data and assigning them to corresponding lithological bodies in the 2D model, the numerical solution of the wave equation is used to synthesize seismic records, accurately reflecting seismic attributes such as amplitude, phase, and frequency during seismic wave propagation.
The wave equation in seismic forward modeling mainly considers the dynamic properties of seismic waves. The propagation law of seismic waves in underground media can be approximately represented by the following acoustic wave equation (Xie et al., 2015):
(1)
(2)
(3)
(4)
(u: wavefield value; v: wave velocity at point (x, y); S(t): source function; l: grid interval distance; ∆t: unit sampling step).
By combining Equations (1)-(4), the following is obtained:
(5)
When the source function, grid interval distance, unit sampling step, and wave velocity are determined, Equation (5) can be used to recursively calculate the wavefield value at any sampling point at any time in the established 2D geological model. After processing, the final synthetic seismic forward record can be obtained.
4. Forward Model Establishment and Result Analysis
4.1. Forward Model Parameter Selection
Based on well logging data from the study area, two combinations of “sand-coal-sand” and “sand-coal interbed” were summarized for sand, coal, and mud in the study area (Figure 3). By statistically analyzing data from 15 wells in the study area, combined with the high-velocity coal seam characteristics of the study area, the P-wave velocity, density, and wave impedance of sandstone, mudstone, and coal seams in the Pinghu Formation were obtained, laying the foundation for the establishment of seismic forward geological models (Table 1).
Table 1. Statistical table of geophysical parameters.
Lithology |
Velocity (m/s) |
Density (g/cm3) |
Acoustic Impedance (10⁷ m/s∙kg/m3) |
Sandstone |
3907 |
2.373 |
9273 |
Mudstone |
3628 |
2.654 |
9630 |
Coal Seam |
4118 |
2.608 |
10,742 |
Figure 3. Combination form of sand and coal. (a) Sand-coal-sand; (b) Sand coal interlayer.
4.2. Forward Model Establishment
Based on the above geophysical parameter statistics, well logging facies lithology configuration, sand-coal distance, and sand, coal, and mud combination types, four theoretical forward models were established. The horizontal width of the models is 1000 m. To explore the influence of coal seams on seismic reflection characteristics in different combination types, Model 1 (Figure 5(a)) was established with lithology settings from bottom to top as follows: first layer mudstone, second layer sandstone, third layer coal seam (thickness gradually increasing from 0 - 5 m), fourth layer sandstone, fifth layer mudstone; the horizontal 0 - 400 m is set as no coal seam, and the horizontal 400 - 1000 m is set as the “sand-coal-sand” combination (Figure 3(a)). Model 2 (Figure 5(b)) was established with lithology settings from bottom to top as follows: first layer mudstone, second layer sandstone, third layer coal seam (thickness gradually increasing from 0 - 5 m), fourth layer sandstone, fifth layer coal seam (thickness gradually increasing from 0 - 5 m), sixth layer sandstone, seventh layer mudstone; the horizontal 0 - 400 m is set as no coal seam, and the horizontal 400 - 1000 m is set as the “sand-coal interbed” combination (Figure 3(b)). To explore the influence of the distance between the underlying coal seam and sandstone on seismic reflection characteristics in different combination types, Model 3 (Figure 5(c)) was established with lithology settings from bottom to top as follows: first layer mudstone, second layer sandstone, third layer mudstone, fourth layer coal seam with a thickness of 5 m and a distance from the overlying sandstone gradually increasing from 0 - 20 m, fifth layer mudstone. Model 4 (Figure 5(d)) was established with lithology settings as follows: first layer mudstone, second layer coal seam with a thickness of 5 m and a distance from the underlying sandstone gradually increasing from 0 - 20 m, third layer mudstone, fourth layer sandstone, fifth layer mudstone. The relationship between coal seam thickness and horizontal distance in the four models (Figure 4(a)), and the relationship between coal seam distance to overlying sandstone and horizontal distance in Model 3 and Model 4 (Figure 4(b)) are as follows.
![]()
Figure 4. The relationship between coal seam thickness, coal seam position, and lateral distance.
4.3. Forward Simulation Analysis
Seismic forward simulations were conducted using zero-phase Ricker wavelets with main frequencies of 30 Hz, 40 Hz, 50 Hz, 60 Hz, and 70 Hz. After setting the source parameters, spatial range parameters, and time sampling interval parameters, forward simulations were performed using the acoustic wave equation. After velocity analysis and Kirchhoff pre-stack migration processing, the final depth-domain synthetic seismic profiles based on Model 1, Model 2, and Model 3 were obtained. Figure 5(a) shows the seismic forward profile based on Model 1 with a wavelet main frequency of 60 Hz, Figure 5(b) shows the seismic forward profile based on Model 2 with a wavelet main frequency of 50 Hz, Figure 5(c) shows the seismic forward profile based on Model 3 with a wavelet main frequency of 60 Hz, and Figure 5(d) shows the seismic forward profile based on Model 4 with a wavelet main frequency of 50 Hz. In Figure 7, AA’ and BB’ correspond to the reflection characteristics of the forward simulation results of Model 1 and Model 2 in the original seismic profile, with their geographical locations shown in Figure 1. Model 3 and Model 4 aim to study the influence of the relative distance between coal seams and overlying or underlying sandstone on seismic configuration characteristics.
Influence of Coal Seam Thickness on Seismic Reflection Configuration
Analyzing coal seam thickness variations based on different models. For Model 1 (Figure 5(a)), the sandstone-coal combination is “sandstone-coal-sandstone”. When the coal seam thickness is 0 - 1 m, sandstone-coal interface reflection characteristics emerge, with minimal enhancement in seismic reflection amplitude. As the coal seam thickness increases, strong amplitudes appear in the sandstone-coal interface seismic reflections, but no phase inversion of the seismic signal occurs. The top surface of the coal seam reflects as a wave crest, while the bottom surface reflects as a wave trough. Throughout this process, the reflection characteristics of the top and bottom interfaces of the sandstone-mudstone remain unchanged (Figure 6(a)), and its true seismic reflection characteristics correspond to Profile AA’ (Figure 7).
Figure 5. Theoretical forward model.
Figure 6. Forward modeling results.
Figure 7. Original seismic profile.
For Model 2 (Figure 5(b)), the sandstone-coal combination is “interbedded sandstone and coal”. Within the horizontal range of 400 - 570 m, corresponding to a coal seam thickness of 0 - 1.4 m, sandstone-coal interface reflection characteristics appear, with minimal enhancement in seismic reflection amplitude. The amplitude at the interface between the upper coal seam and sandstone exhibits a wave crest, while the lower one exhibits a wave trough. The reflection amplitudes of both the top and bottom interfaces of the sandstone-mudstone change from wave crests to wave troughs, with a slight enhancement in amplitude. Within the horizontal range of 570 - 670 m, corresponding to a coal seam thickness of 1.4 - 2.25 m, due to the thin coal seams and the interference between upper and lower coal seams, the sandstone between coal seams exhibits no reflection characteristics. Meanwhile, the reflection characteristics of the upper coal seam are also affected by the wave reflection of the lower coal seam, resulting in weakened amplitude. The sandstone-coal reflection interface still exhibits a wave crest, while the reflection characteristic of the sandstone-coal bottom interface changes from a wave crest to a wave trough. Both the amplitudes of the top and bottom interfaces of the sandstone-mudstone exhibit wave troughs. Within the horizontal range of 670 - 830 m, corresponding to a coal seam thickness of 2.25 - 3.6 m, the amplitude at the reflection interface between the upper coal seam and sandstone increases and exhibits wave crest reflection. The amplitude of the top interface of the sandstone-mudstone increases, while the amplitude of the bottom interface decreases, with both top and bottom interface amplitudes exhibiting wave troughs. Within the horizontal range of 830 - 930 m, corresponding to a coal seam thickness of 3.6 - 4.4 m, the amplitudes of the top and bottom interfaces of the upper sandstone-coal decrease, while the reflection amplitude intensity at the top and bottom interfaces of the lower sandstone-coal shows no significant change. The amplitude of the top interface of the sandstone-coal changes from a wave trough to a wave crest, and the lower one changes from a wave crest to a wave trough. Within the horizontal range of 930 - 1000 m, corresponding to a coal seam thickness of 4.4 - 5 m, the amplitudes of the top and bottom interfaces of the upper sandstone-coal increase. The amplitude intensity of the top interface of the lower sandstone-coal shows no significant change, exhibiting a change from a wave crest to a wave trough, while the amplitude intensity of the bottom interface decreases. The amplitude intensities of both the top and bottom interfaces of the sandstone-mudstone increase (Figure 6(b)), and its true seismic reflection characteristics correspond to Profile BB’ (Figure 7).
5. Conclusion and Insights
Based on the geological features of coal measures in the Pinghu Formation of Pingbei, Xihu Sag, four forward models were built and analyzed using acoustic wave equations. The following insights were obtained:
1) In the “sandstone-coal-sandstone” combination, as coal seam thickness increases, the reflection amplitude intensity of sandstone-mudstone interfaces augments without polarity reversal.
2) In the “interbedded sandstone and coal” scenario, increasing coal seam thickness boosts reflection amplitude intensity at sandstone-mudstone interfaces. Multiple polarity reversals occur when coal thickness ranges from 0 - 5 m, attributed to interference from multiple coal seams.
3) The influence of coal seams on sandstone reflections diminishes as their distance increases. Between 0 - 10 m, reflection amplitude at sandstone-mudstone interfaces intensifies without polarity reversal. Between 10 - 20 m, amplitude decreases and polarity reversal occurs.
Acknowledgements
Thanks to the school and colleagues for their help, as well as the strong support for me.