Application of 9-component S-wave 3D seismic data to study sedimentary facies and reservoirs in a biogasbearing area:A case study on the Pleistocene Qigequan Formation in Taidong area,Sanhu Depression,Qaidam Basin,NW China

To solve the problems in restoring sedimentary facies and predicting reservoirs in loose gas-bearing sediment,based on seismic sedimentologic analysis of the first 9-component S-wave 3D seismic dataset of China,a fourth-order isochronous stratigraphic framework was set up and then sedimentary facies and reservoirs in the Pleistocene Qigequan Formation in Taidong area of Qaidam Basin were studied by seismic geomorphology and seismic lithology.The study method and thought are as following.Firstly,techniques of phase rotation,frequency decomposition and fusion,and stratal slicing were applied to the 9-component S-wave seismic data to restore sedimentary facies of major marker beds based on sedimentary models reflected by satellite images.Then,techniques of seismic attribute extraction,principal component analysis,and random fitting were applied to calculate the reservoir thickness and physical parameters of a key sandbody,and the results are satisfactory and confirmed by blind testing wells.Study results reveal that the dominant sedimentary facies in the Qigequan Formation within the study area are delta front and shallow lake.The RGB fused slices indicate that there are two cycles with three sets of underwater distributary channel systems in one period.Among them,sandstones in the distributary channels of middle-low Qigequan Formation are thick and broad with superior physical properties,which are favorable reservoirs.The reservoir permeability is also affected by diagenesis.Distributary channel sandstone reservoirs extend further to the west of Sebei-1 gas field,which provides a basis to expand exploration to the western peripheral area.

A Method of Improving Seismic Data Resolution：Comprehensive Inversion of Well logging and Seismic Data

Comprehensive inversion of logging and seismic data presented in this paper is a method to improve seismic data resolution. It involves using ample high-frequency information and complete low-frequency information of known logging to make up for the lack of limited bandwidth of practical seismic recording, obtaining an approximate reflection coefficient sequence (or wave impedance) of high resolution by iterative inversion and providing more reliable seismic evidence for further lithologic interpretation and lateral tracking, correlation and prediction of thin reservoir. The comprehensive inversion can be realized in the following steps: (1) to establish an initial model of higher resolution; (2) to obtain wavelets, and (3) to constrain iterative inversion. The key to this inversion lies in building an initial model. It is assumed from our experience that when the initial model is properly given, iterative inversion can be quickly converged to the ideal result.

Integration of Well Logs and Seismic Data for Prospects Evaluation of an X Field, Onshore Niger Delta, Nigeria

Hydrocarbon reservoir beds have been delineated using direct hydrocarbon indicator on seismic sections as well as well logs data in X field, Onshore Niger Delta. The research methodology involved horizon interpretation to produce sub-surface structure map. Geophysical well log signatures were employed in identifying hydrocarbon bearing sand. The well-to-seismic tie revealed that the reservoir tied directly with hydrocarbon indicator (bright spot) on the seismic sections. The major structure responsible for the hydrocarbon entrapment is anticline. The crest of the anticline from the depth structural map occurs at 3450 metres.

Separation of P. and SV-wavefields from multi-componen seismic data

In multi-component seismic exploration, the horizontal and vertical components both contain P- and SV-waves. The P- and SV-wavefields in a seismic record can be separated by their horizontal and vertical displacements when upgoing P- and SV-waves arrive at the sea floor. If the sea floor P wave velocity, S wave velocity, and density are known, the separation can be achieved in ther-p domain. The separated wavefields are then transformed to the time domain. A method of separating P- and SV-wavefields is presented in this paper and used to effectively separate P- and SV-wavefields in synthetic and real data. The application to real data shows that this method is feasible and effective. It also can be used for free surface data.

Robust high frequency seismic bandwidth extension with a deep neural network trained using synthetic data

Geophysicists interpreting seismic reflection data aim for the highest resolution possible as this facilitates the interpretation and discrimination of subtle geological features.Various deterministic methods based on Wiener filtering exist to increase the temporal frequency bandwidth and compress the seismic wavelet in a process called spectral shaping.Auto-encoder neural networks with convolutional layers have been applied to this problem,with encouraging results,but the problem of generalization to unseen data remains.Most published works have used supervised learning with training data constructed from field seismic data or synthetic seismic data generated based on measured well logs or based on seismic wavefield modelling.This leads to satisfactory results on datasets similar to the training data but requires re-training of the networks for unseen data with different characteristics.In this work seek to improve the generalization,not by experimenting with network architecture(we use a conventional U-net with some small modifications),but by adopting a different approach to creating the training data for the supervised learning process.Although the network is important,at this stage of development we see more improvement in prediction results by altering the design of the training data than by architectural changes.The approach we take is to create synthetic training data consisting of simple geometric shapes convolved with a seismic wavelet.We created a very diverse training dataset consisting of 9000 seismic images with between 5 and 300 seismic events resembling seismic reflections that have geophysically motived perturbations in terms of shape and character.The 2D U-net we have trained can boost robustly and recursively the dominant frequency by 50%.We demonstrate this on unseen field data with different bandwidths and signal-to-noise ratios.Additionally,this 2D U-net can handle non-stationary wavelets and overlapping events of different bandwidth without creating excessive ringing.It is also robust in the presence of noise.The significance of this result is that it simplifies the effort of bandwidth extension and demonstrates the usefulness of auto-encoder neural network for geophysical data processing.

The Analysis of Seismic Data Structure and Oil and Gas Prediction

In this paper, a new concept called numerical structure of seismic data is introduced and the difference between numerical structure and numerical value of seismic data is explained. Our study shows that the numerical seismic structure is closely related to oil and gas-bearing reservoir, so it is very useful for a geologist or a geophysicist to precisely interpret the oil-bearing layers from the seismic data. This technology can be applied to any exploration or production stage. The new method has been tested on a series of exploratory or development wells and proved to be reliable in China. Hydrocarbon-detection with this new method for 39 exploration wells on 25 structures indi- cates a success ratio of over 80 percent. The new method of hydrocarbon prediction can be applied for: (1) depositional environment of reservoirs with marine fades, delta, or non-marine fades (including fluvial facies, lacustrine fades); (2) sedimentary rocks of reservoirs that are non-marine clastic rocks and carbonate rock; and (3) burial depths range from 300 m to 7000 m, and the minimum thickness of these reservoirs is over 8 m (main frequency is about 50 Hz).

Time-domain identification of dynamic properties of layered soil by using extended Kalman filter and recorded seismic data

A novel time-domain identification technique is developed for the seismic response analysis of soil-structure interaction.A two-degree-of-freedom (2DOF) model with eight lumped parameters is adopted to model the frequency- dependent behavior of soils.For layered soil,the equivalent eight parameters of the 2DOF model are identified by the extended Kalman filter (EKF) method using recorded seismic data.The polynomial approximations for derivation of state estimators are applied in the EKF procedure.A realistic identification example is given for the layered-soil of a building site in Anchorage,Alaska in the United States.Results of the example demonstrate the feasibility and practicality of the proposed identification technique.The 2DOF soil model and the identification technique can be used for nonlinear response analysis of soil-structure interaction in the time-domain for layered or complex soil conditions.The identified parameters can be stored in a database tor use in other similar soil conditions,lfa universal database that covers information related to most soil conditions is developed in the thture,engineers could conveniently perform time history analyses of soil-structural interaction.

Resolving co- and early post-seismic slip variations of the 2021 MW 7.4 Madoi earthquake in east Bayan Har block with a block-wide distributed deformation mode from satellite synthetic aperture radar data

On 21 May 2021(UTC),an MW 7.4 earthquake jolted the east Bayan Har block in the Tibetan Plateau.The earthquake received widespread attention as it is the largest event in the Tibetan Plateau and its surroundings since the 2008 Wenchuan earthquake,and especially in proximity to the seismic gaps on the east Kunlun fault.Here we use satellite interferometric synthetic aperture radar data and subpixel offset observations along the range directions to characterize the coseismic deformation of the earthquake.Range offset displacements depict clear surface ruptures with a total length of~170 km involving two possible activated fault segments in the earthquake.Coseismic modeling results indicate that the earthquake was dominated by left-lateral strike-slip motions of up to 7 m within the top 12 km of the crust.The well-resolved slip variations are characterized by five major slip patches along strike and 64%of shallow slip deficit,suggesting a young seismogenic structure.Spatial-temporal changes of the postseismic deformation are mapped from early 6-day and 24-day InSAR observations,and are well explained by time-dependent afterslip models.Analysis of Global Navigation Satellite System(GNSS)velocity profiles and strain rates suggests that the eastward extrusion of plateau is diffusely distributed across the east Bayan Har block,but exhibits significant lateral heterogeneities,as evidenced by magnetotelluric observations.The block-wide distributed deformation of the east Bayan Har block along with the significant co-and post-seismic stress loadings from the Madoi earthquake imply high seismic risks along regional faults,especially the Tuosuo Lake and Maqên-Maqu segments of the Kunlun fault that are known as seismic gaps.

A novel technique for automatic seismic data processing using both integral and local feature of seismograms

A novel technique for automatic seismic data processing using both integral and local feature of seismograms was presented in this paper. Here, the term integral feature of seismograms refers to feature which may depict the shape of the whole seismograms. However, unlike some previous efforts which completely abandon the DIAL approach, i.e., signal detection, phase identifi- cation, association, and event localization, and seek to use envelope cross-correlation to detect seismic events directly, our technique keeps following the DIAL approach, but in addition to detect signals corresponding to individual seismic phases, it also detects continuous wave-trains and explores their feature for phase-type identification and signal association. More concrete ideas about how to define wave-trains and combine them with various detections, as well as how to measure and utilize their feature in the seismic data processing were expatiated in the paper. This approach has been applied to the routine data processing by us for years, and test results for a 16 days＇ period using data from the Xinjiang seismic station network were presented. The automatic processing results have fairly low false and missed event rate simultaneously, showing that the new technique has good application prospects for improvement of the automatic seismic data processing.

Seismic Signal and Data Analysis of Rock Media with Vertical Anisotropy

This paper is concerned with anisotropic effects on seismic data and signal analysis for transversely isotropic rock media with vertical anisotropy. It is understood that these effects are significant in many practical applications, e.g. earthquake forecasting, materials exploration inside the Earth’s crust, as well as various practical works in oil industry. Under the framework of the most accepted anisotropic media model (i.e. VTI media, transverse isotropy with a vertical axis symmetry), with applications of a set of available anisotropic rock parameters for sandstone and shale, we have performed numerical calculations of the anisotropic effects. We show that for rocks with strong anisotropy, the induced relative depth error can be significantly large. Nevertheless, with an improved understanding of the seismic-signal propagation and proper data processing, the error can be reduced, which in turn may enhance the probability of forecasting accurately the various wave propagations inside the Earth’s crust, e.g. correctly forecasting the incoming earthquakes from the center of the Earth.

Sedimentary Microfacies and Porosity Modeling of Deep-Water Sandy Debris Flows by Combining Sedimentary Patterns with Seismic Data： An Example from Unit I of Gas Field A, South China Sea

Sandy debris flow deposits are present in Unit I during Miocene of Gas Field A in the Baiyun Depression of the South China Sea. The paucity of well data and the great variability of the sedimentary microfacies make it difficult to identify and predict the distribution patterns of the main gas reservoir, and have seriously hindered further exploration and development of the gas field. Therefore, making full use of the available seismic data is extremely important for predicting the spatial distribution of sedimentary microfacies when constructing three-dimensional reservoir models. A suitable reservoir modeling strategy or workflow controlled by sedimentary microfacies and seismic data has been developed. Five types of seismic attributes were selected to correlate with the sand percentage, and the root mean square （RMS） amplitude performed the best. The relation between the RMS amplitude and the sand percentage was used to construct a reservoir sand distribution map. Three types of main sedimentary microfacies were identified： debris channels, fan lobes, and natural levees. Using constraints from the sedimentary microfacies boundaries, a sedimentary microfacies model was constructed using the sequential indicator and assigned value simulation methods. Finally, reservoir models of physical properties for sandy debris flow deposits controlled by sedimentary microfacies and seismic inversion data were established. Property cutoff values were adopted because the sedimentary microfacies and the reservoir properties from well-logging interpretation are intrinsically different. Selection of appropriate reservoir property cutoffs is a key step in reservoir modeling when using simulation methods based on sedimentary microfacies control. When the abnormal data are truncated and the reservoir properties probability distribution fits a normal distribution, microfacies-controlled reservoir property models are more reliable than those obtained from the sequence Gauss simulation method. The cutoffs for effective porosity of the debris channel, fan lobe, and natural levee facies were 0.2, 0.09, and 0.12, respectively; the corresponding average effective porosities were 0.24, 0.13, and 0.15. The proposed modeling method makes full use of seismic attributes and seismic inversion data, and also makes the property data of single-well depositional microfacies more conformable to a normal distribution with geological significance. Thus, the method allows use of more reliable input data when we construct a model of a sandy debris flow.

Velocity structure in the South Yellow Sea basin based on first-arrival tomography of wide-angle seismic data and its geological implications

The South Yellow Sea basin is filled with Mesozoic-Cenozoic continental sediments overlying pre-Palaeozoic and Mesozoic-Palaeozoic marine sediments.Conventional multi-channel seismic data cannot describe the velocity structure of the marine residual basin in detail,leading to the lack of a deeper understanding of the distribution and lithology owing to strong energy shielding on the top interface of marine sediments.In this study,we present seismic tomography data from ocean bottom seismographs that describe the NEE-trending velocity distributions of the basin.The results indicate that strong velocity variations occur at shallow crustal levels.Horizontal velocity bodies show good correlation with surface geological features,and multi-layer features exist in the vertical velocity framework(depth:0–10 km).The analyses of the velocity model,gravity data,magnetic data,multichannel seismic profiles,and drilling data showed that high-velocity anomalies(>6.5 km/s)of small(thickness:1–2 km)and large(thickness:>5 km)scales were caused by igneous complexes in the multi-layer structure,which were active during the Palaeogene.Possible locations of good Mesozoic and Palaeozoic marine strata are limited to the Central Uplift and the western part of the Northern Depression along the wide-angle ocean bottom seismograph array.Following the Indosinian movement,a strong compression existed in the Northern Depression during the extensional phase that caused the formation of folds in the middle of the survey line.This study is useful for reconstructing the regional tectonic evolution and delineating the distribution of the marine residual basin in the South Yellow Sea basin.

Clastic compaction unit classification based on clay content and integrated compaction recovery using well and seismic data

Compaction correction is a key part of paleogeomorphic recovery methods. Yet, the influence of lithology on the porosity evolution is not usually taken into account. Present methods merely classify the lithologies as sandstone and mudstone to undertake separate porositydepth compaction modeling. However, using just two lithologies is an oversimplification that cannot represent the compaction history. In such schemes, the precision of the compaction recovery is inadequate. To improve the precision of compaction recovery, a depth compaction model has been proposed that involves both porosity and clay content. A clastic lithological compaction unit classification method, based on clay content, has been designed to identify lithological boundaries and establish sets of compaction units. Also, on the basis of the clastic compaction unit classification, two methods of compaction recovery that integrate well and seismic data are employed to extrapolate well-based compaction information outward along seismic lines and recover the paleo-topography of the clastic strata in the region. The examples presented here show that a better understanding of paleo-geomorphology can be gained by applying the proposed compaction recovery technology.

Determination of combined effect of ampli- fication and attenuation of soft rock site us- ing digital seismic data

In the paper, for the application of stochastic simulation of ground motion, we put forward a method to determine ″the combined effect of amplification and attenuation″ (combined effect for short) of soft rock site by using digital seismic data of moderate and small earthquakes. Our approach aims at solving the problem of the combined effect of soft rock site, which is difficult to determine in most regions of China because fewer measures were done for S-wave velocity structure. The combined effect of soft rock site can be determined by using the approach recom- mended by us. An example is given to discuss the practical application of the method.

The movement age of hidden fault and analysis on width of its effect zone from shallow seismic sounding and drilling data

Field experimental seismic sounding permitted us to obtain optimal shallow seismic reflection sounding parameters.In process of data processing,we obtained a high-qualitative shallow seismic reflection sounding profile by using the techniques such as filtering,edition surgical blanking,prediction deconvolution,fitting static correlation of first arrival time,and velocity analysis.Comprehensive analysis on the information of reflection wave groups along the seismic sounding profile and the stratigraphic and neogeochronological data obtained from many drills near the sounding line reveals that the upper termination of the detected fault zone is located at depth of 75~80 m,in the Middle Pleistocene deposits dated to be about 220 ka BP.The continuity,discontinuity,increasing and decreasing amount of reflection wave groups and change of their configurations,in combination with geological columns of drills,permitted us to know that the width of upper termination of the fault zone is 100 m.It can be inferred from the variation of number of reflection wave groups along the profile that the scarp of hidden fault is 200 m wide and the fault is a synsedimentary active fault in the Early Pleistocene and the early stage of Middle Pleistocene.No tectonic movement,which offset the covering deposits,had occurred since the late stage of Middle Pleistocene.

LATE QUATERNARY SEISMIC SEQUENCES AND SEDIMENTATION ALONG THE COAST OFF TAISHAN,GUANGDONG,AS REVEALED BY HIGH RESOLUTION SEISMIC DATA

Seismic stratigraphic analysis of 150 km of high resolution seismic profiles (GEOPULSE and 3.5 kHz) obtained along the muddy coast off Taishan, Guangdong provided valuable information on the offshore Quaternary stratigraphy, tectonics, sedimentary evolution, and environmental and engineering geology of this area.The following late Quaternary stratigraphic sequences were mapped. (1) Pleistocene alluvial silty sand and shallow marine clay; (2) coarse fluvial channel sand; (3) nearshore and bay clayey sediments. The incision of fluvial channels occurred in two stages: during the Late Pleistocene and at the beginning of the Holocene. The Pleistocene fluvial channels have been filled by coarse sand overlain by more than 10 m of late Pleistocene clay and 20-30 m of Holocene clay.The coastal area is dissected by two major fault systems: the first group strikes ENE-WSW and controls the evolution of the adjacent Quaternary basin, while the second one trends NNE-SSW. These fault systems must be taken into account in any coastal engineering considerations. In addition, both granitic and metamorphic basements were recognized.

Integration of Geological, Geophysical and Seismological Data for Seismic Hazard Assessment Using Spatial Matching Index

Probabilistic seismic hazard assessment (PSHA) takes into account as much data as possible for defining the initial seismic source zone model. In response to this, an algorithm has been developed for integration of geological, geophysical and seismological data through a spatial index showing the presence or absence of a potential seismic source feature in the input data. The spatial matching index (SMI) is calculated to define the coincidence of independent data showing any indications for existence of a fault structure. It is applied for hazard assessment of Bulgaria through quantification of the seismic potential of 416 square blocks, 20 × 20 km in size covering the entire territory of Bulgaria and extended by 20 km outside of the country borders. All operations are carried out in GIS environment using its capabilities to work with different types of georeferenced spatial data. Results show that the highest seismic potential (largest SMI) is observed in 56 block elements (13% of the territory) clearly delineating cores of the source zones. Partial match is registered in 98 block elements when one of the features is missing. Not any evidence for earthquake occurrence is predicted by our calculation in 117 elements, comprising 28% of the examined area. The quantitative parameter for spatial data integration which is obtained in the present research may be used to analyze information regardless of its type and purpose.

Combination of structure tensor and tilt angle in the edge detection of potential field data

Edge detection is a commonly requested task in the interpretation of potential field data. Different methods have different results for varied depths and shapes of geological bodies. In this paper,we propose using the combination of structure tensor and tilt angle to detect the edges of the sources,which can display the edges of shallow and deep bodies simultaneously. Through tests on synthetic potential field data,it is obvious that the proposed edge detection methods can display the sources edges more clearly and precisely,compared with other commonly used methods. The application on real potential field data shows similar result,obtaining the edges of layers and faults clearly. In addition,another advantage of the new method is its insensitivity to noise.

Structure and Temporal Variability of Mediterranean Water in Hydrological and Marine Seismic Data South of Portimao Canyon (Gulf of Cadiz), from 1999 to 2002

Hydrological and marine seismic data, collected in the Gulf of Cadiz (respectively in July 1999, 2000, 2001 and 2002, and in April 2000 and 2001) are analysed to reveal the various structures of Mediterranean Water (MW). Both the hydrological and seismic data clearly identify the MW undercurrents on the Iberian slope, detached MW eddies (meddies and a cyclone) and smaller fragments of MW (filaments and small eddies). Seismic reflectivity and synthetic reflectivity computed from hydrology, indicate that strong acoustic reflectors, associated with 8 - 64 m thick homogeneous water layers, are found above and below meddies and filaments, around the MW undercurrents, but mostly in the lower part of cyclones and below submesoscale eddies. Reflectors are also observed in the near surface layers where thermohaline contrasts are quite pronounced. The successful use of seismic data to locate submesoscale MW structures, superior to that of hydrology, is related to the improved horizontal resolution.

Tectonics of the Solomon Sea Basin from Vertical Gravity Gradient and Seismic Data

The Solomon Sea Basin is a Cenozoic back-arc spreading basin within the convergence system of the Pacific and Indo-Australian plates.Against the background of subduction polarity reversal,the current Solomon Sea Basin gradually formed a rhombic morphology with the subduction of the basin along the New Britain Trench and the Trobriand Trough.By analyzing the vertical gravity gradient,natural earthquake and seismic reflection data,this study determines the structural characteristics of the Solomon Sea Basin.It was found that the tectonics of the basin are characterized by the original expansion structure within the central part in addition to the structure induced by the latest subduction along the basin margin.The original spreading structure of the basin presented an east–west linear graben and horst controlled by normal faults during the basin expansion period.As a result of the subduction and slab-pull of the Solomon Sea Basin,extensional structure belts parallel to the New Britain Trench formed along the basin margin.