Prestack Depth Migration Research Articles

Seismic Exploration Methods for Structural Studies and for Active Fault Characterization: A Review

In this paper, seismic exploration methods are reviewed with a particular emphasis on the use of the reflection seismology to investigate the subsurface structures and characterize active faults. The paper provides a descriptive overview, intended for a non-specialist audience, of the methods and of their recent developments aimed at improving the resolution, accuracy, and computational efficiency of seismic imaging. Techniques such as seismic ray tomography, full-waveform inversion and pre-stack depth migration are briefly introduced, highlighting their potential applications in structural geology studies. The main seismic attributes that have become increasingly important in the interpretation of faults and fractures are also presented, along with some examples of application. Finally, some case studies of active fault characterization are discussed. From these examples, the crucial role played nowadays by the seismic exploration methods for structural studies and for active fault characterization is evident.

High-precision and high-efficiency first-arrival slope tomography via eikonal solvers and the adjoint-state method

Abstract First-arrival slope tomography (FAST) introduces first-arrival slopes, corresponding to the horizontal components of the slowness vectors at the receiver and source positions to supplement first-arrival traveltime for better guiding ray propagation in the media until the best match is achieved with the observed data. FAST can recover the velocity model with higher resolution and precision than first-arrival traveltime tomography (FATT) but is computationally intensive. In this context, we propose an improved approach, referred to as high-precision and high-efficiency first-arrival slope tomography (HFAST). HFAST redefines one of the slopes using the reciprocity principle and simultaneously employs the first-arrival traveltime and slopes to ensure high-quality model building. On the other hand, HFAST extracts calculated data and derives the gradient of the misfit function from the solutions of relatively limited forward and inverse problems, resulting in a low computational cost. The cost of HFAST is proportional to the minimum between the receivers and sources, whereas the cost of FAST is scaled to the sum of the receivers and sources. Numerical experiments involving the checkerboard and SEAM II Foothill models demonstrate that HFAST can achieve a higher inversion precision than FATT, especially in the recovery of small-scale anomalies and the presence of velocity reversal. Moreover, HFAST is more computationally efficient than FAST and suitable for managing large data sets. Therefore, HFAST can be regarded as a valuable supplement to current first-arrival-based model building methods and has the potential to be applied in static corrections, prestack depth migration and waveform inversion in the future.

Wide-angle seismic imaging of divergent and transform segments of the Pará-Maranhão-Barreirinhas-Ceará margin, NW Brazil

The structure of the North-East equatorial Brazilian margin was investigated during the MAGIC (Margins of brAzil, Ghana and Ivory Coast) seismic wide-angle experiment. This study focuses on the MC5 profile, that spans NW-SE 720 km in length, from the São Paulo Double Fracture Zone to the Barreirinhas margin and continental Borborema province. Its main objective is to understand the fundamental processes which lead to the thinning and finally to the breakup of the continental crust in a specific context of a divergent pull-apart system with two strike-slip borders. The experiment was devised to obtain the 2D structure along this profiles from joint pre-stack depth migration of streamer data and travel-time inversion by forward modeling of 43 Ocean Bottom Seismometers and 21 Land Seismic Stations records. Along the MC5 wide-angle transect, 4 major sectors are identified: 1) the São Paulo Double Fracture Zone presenting a 4.5 km thick volcano-sedimentary Basin on top of a 5.5 km thick basement; 2) a volcanic alignment and intermediate domain SE-ward, formed by the 4.5 km thick Basin III; 3) the 7.5 km thick Basin II, and the 5.5 km thick Basin I composing the continental slope and shelf. While all the offshore basement remains about 6 km thick in the deep-sea domains, acoustic velocity evolves from two-layer 4.8–6 km/s and 6.1–6.8 km/s beneath Basin III to two-layer high velocity 6.1–6.8 km/s and 7.2–7.4 km/s beneath Basin II and I, The necking zone, forming the Parnaiba Platform and associated Piaui-Camocim and Ceará Basins, is 50 km wide; 4) the Médio Coreaú and Ceará Central thrust belt, where the un-thinned continental crust thickness reaches 32 km. Finally, a schematic kinematic reconstruction that satisfies these observation is argumented.

Characteristic reflector-based wave equation reflection traveltime inversion

Reconstructing the deep background velocity model is crucial for imaging subsurface structures. The reflection traveltime inversion (RTI) is capable of recovering the kinematic information in the velocity model under a limited seismic observation aperture. The estimation of the reflection traveltime residuals between the synthetic and observed seismic waves is the key to the RTI. However, it is difficult to accurately obtain the rapidly changing reflection traveltime shifts in the prestack data domain. Estimating the traveltime difference is especially difficult where the discrepancy between the noisy observed data and demigrated synthetic data includes a waveform distortion, inconsistent reflected events, and so on. To overcome this problem, we have developed a new reflection traveltime residual estimation method. In this method, after identifying a characteristic reflector in the prestack depth migration section with an initial velocity, the demigrated and observed seismic reflections in the data domain corresponding to the characteristic reflector share the same traveltime at the zero offset. Thus, their traveltime residual can be automatically tracked by their crosscorrelation with a modified dynamic programming algorithm based on a zero-offset control point, which can avoid traveltime picking or windowing of relevant reflection arrivals in the prestack seismic data domain. Then, the estimated reflection traveltime residual is applied to the RTI. To improve the inversion accuracy, we have further developed a layer-stripping velocity model-building workflow. Synthetic and field data tests demonstrate that the proposed method can obtain a robust reflection traveltime residual and further achieve a correct low-to-intermediate wavenumber velocity model update.

A hypothesis of the obliquely subducted Gagua Ridge below Nanao Basin based on uplifted structures in Ryukyu forearc

At the southernmost part of the Ryukyu subduction zone, six long-offset multi-channel seismic profiles were collected across three forearc basins and the southern Ryukyu accretionary prism during the TAIGER experiment in 2009 and the TAICRUST project in 1995. These profiles were reprocessed to generate pre-stack depth migration (PSDM) sections. In addition, two velocity-interface models were obtained by reanalyzing active source data recorded from 28 ocean-bottom seismometers during the same TAIGER experiment, in consideration of the PSDM sections and previous tomography models. Due to the northwest convergence of the Philippine Sea Plate (PSP), it is suggested that the Gagua Ridge may have been obliquely subducting northwestward beneath the Ryukyu prism and below the Nanao Basin. The PSDM sections and the velocity-interface models indicate the subducted Gagua Ridge causing the uplift of the sedimentary basement and the lower crustal structure below the Nanao Basin. The sedimentary and crustal structures near ~ 122.5° E beneath the Nanao Basin were also uplifted where shallow earthquakes had occurred by the oblique subduction of the Gagua Ridge depicted in the isopach map of the crust. Furthermore, a recent earthquake (Mw ~ 6.0) occurred near the northeast coast of Taiwan in 2018 at a depth of approximately 12 km below the Hsincheng Ridge. Our study suggests that this earthquake was caused by a thrust fault near the décollement, which might have been formed by the subduction of the PSP.

Imaging exhumed continental and proto-oceanic crusts in the Camamu triple junction, Brazil

During the SALSA experiment, in 2014, twelve combined wide-angle refraction and coincident multi-channel seismic profiles were acquired in the Jequitinhonha-Almada-Camamu, Jacuípe, and Sergipe-Alagoas basins, NE Brazil. Profile SL09 images the Almada-Camamu basin and the São Francisco craton, with 18 four-channel ocean-bottom seismometers and 22 land stations. The datasets were forward modelled and combined with pre-stack depth migration to increase the horizontal resolution of the velocity models.Our results show that sediment thickness varies between 3.8 km in the oceanward part of the profile, 4.3 km in the Almada basin and 6.5 km in the Camamu basin. Crustal thickness at the north-western edge of the profile is of around 40 km, with velocity gradients indicating a continental origin.The Camamu basin, which corresponds to the triple junction between the aborted N–S oriented Tucano rift, the SW-NE oriented Jacuipe-Sergipe-Alagoas branch, and the N–S Jequitinhonha-Almada branch, presents two crustal layers: a very thin upper layer, about 1.5 km thick, which increases seawards to 3 km in the Almada basin, and an higher velocity (HV) layer (6.8–7.2 km/s) about 4 km thick. This lower layer gradually disappears in the Almada basin. At the south-eastern edge of the profile, the resolution is lower but the thickness of the crust seems to increase up to 5 km. Deep wide-angle reflections indicate upper mantle stratification.Crustal organisation and P-wave propagation velocities in the Almada and Camamu basins indicate a transitional crust domain of exhumed continental crust affinity. In the Camamu triple junction and beneath this thin exhumed continental crust, the HV layer may probably reflect intruded materials. No exhumed upper mantle is observed along the entire profile. The easternmost part of the profile may correspond to a proto-oceanic crust. Typical oceanic crust is never imaged along the 260 km-long offshore profile.

Multi-window spectrum analysis of Búzios oil field PSDM seismic data applied to frequency attenuation study

The intrinsic attenuation is often quantified using the inverse quality factor Q. Methods to estimate Q are mainly based on time domain or frequency domain and sensitive to amplitude and/or frequency band variations. In this sense, noise, signal artifacts or processing steps can alter the amplitude, frequency content and the waveform, what may lead to inaccurate calculations of Q. A way to better understand the frequency band variation in seismic data and verify patterns or problems occurring in the signal frequency band is performing layered or window-by-window spectrum analysis. This work describes a multi-window spectrum analysis method used to better understand the frequency attenuation process as the seismic wave travels through the sediments from the sea bottom down to the reservoir depths in a 3D Kirchhoff Pre-Stack Depth Migrated (PSDM) seismic data from Búzios oil field. To perform the spectral analysis in frequency domain, the PSDM volume was converted to time domain using the same interval velocity volume used to depth-migrate the data. The variation of frequency band as the signal propagates through the sediment package was analyzed dividing the studied crossline sections in 13 layers of 200 ms each, from sea bottom down to the base of pre-salt reservoir formations. The results suggest that most part of higher frequency loss has been occurred in post-salt package. This may be an indication of significant viscoelastic behavior or strong stratigraphic filtering of this interval. Sequentially, higher frequency content is lost slowly from salt top down to the last layer located inside pre-salt formation. Very slight frequency spectrum variation below salt top indicates elastic behavior of the salt formation and underling pre-salt formation in this interval. The spectrum analysis method proposed here has shown to be useful and efficient to identify intervals of different attenuation patterns in seismic data. Additionally, the frequency band study produced by our spectrum analysis can be used to quality-control Q-factor estimation computed by different methods.

Depth migration of crustal-scale seismic reflection profiles: A case study in the Bohai Bay basin

The application of the pre-stack depth migration (PSDM) method using deep seismic reflection data is still a challenge. In this paper, a deep seismic reflection dataset collected from the Bohai Bay Basin is used to test the applicability of the pre-stack depth migration method. After improving the signal-to-noise ratio of the prestack gathers and building an accurate velocity model capturing lateral and vertical variations, a high-quality pre-stack depth migration profile is obtained. This profile shows detailed crustal structures and an undulating Moho. The depth difference between the east and west sides of the Moho imaged in the profile is about 6 km. The Moho is disrupted by several discontinuities. Compared with conventionally stacked time-migration profiles, the pre-stack depth migration profile provides a higher-resolution subsurface image with accurate depth information. Combining geological and geophysical data, the pre-stack depth migration profile reveals geologic structures responsible for crust thinning during the destruction of the North China Craton. The profile also provides more accurate depth domain information for further study of crust-mantle interaction in this area.

Seismic Imaging of Mineral Exploration Targets: Evaluation of Ray- vs. Wave-Equation-Based Pre-Stack Depth Migrations for Crooked 2D Profiles

Seismic imaging is now a well-established method in mineral exploration with many successful case studies. Seismic data are usually imaged in the time domain (post-stack or pre-stack time migration), but recently pre-stack depth imaging has shown clear advantages for irregular/sparse acquisitions and very complex targets. Here, we evaluate the effectiveness of both ray-based and wave-equation-based pre-stack depth imaging methodologies applied to crooked-line 2D seismic reflection profiles. Seismic data were acquired in the Kylylahti mining area in eastern Finland over severely folded, faulted and subvertical Kylylahti structure, and associated mineralization. We performed 3D ray-based imaging, i.e., industry-standard Kirchhoff migration and its improved version (coherency migration, CM), and wave-equation-based migration, i.e., reverse time migration (RTM) using a velocity model built from first-arrival traveltime tomography. Upon comparing the three different migrations against available geological data and models, it appeared that CM provided the least noisy and well-focused image, but failed to image the internal reflectivity of the Kylylahti formation. RTM was the only method that produced geologically plausible reflections inside the Kylylahti formation including a direct image of the previously known shallow massive sulfide mineralization.

Imaging of thrust structures under Colombian foothills landscape through Marchenko approach

The Marchenko imaging approach successfully deals with multiple scattering during seismic migration, obtaining clear images in structurally complex geological zones below at a specific depth. Through the Marchenko equation, the method retrieves Green's function of a virtual source in profundity using reflections recorded on the surface and a smoothed velocity model. According to the theorem of reciprocity, Green's function represents the medium's response to a source on the surface recorded by a virtual receiver in-depth and decomposes it into a down-going wavefield and an up-going wavefield. Both wave fields' deconvolution provides the reflection response of virtual sources and virtual receivers repositioned in the subsurface, mimicking redatumming. The result is an image free of multiple internal reflections from the area below virtual both sources and receivers. After successfully testing with synthetic seismograms generated in a plane-horizontal layers model, we applied the Marchenko algorithm to seismic data of the Upper Magdalena Valley project using an available velocity model. Compared to an image supplied by a Kirchhoff pre-stack depth migration algorithm, the new image correctly positions the seismic events avoiding misleading interpretations in locating complex structures.

Offshore subsurface characterization enabled by fiber-optic distributed acoustic sensing (DAS): An East China Sea 3D VSP survey example

Recent advances in distributed acoustic sensing (DAS) technology have allowed more intense measurements of subsurface and environment events, providing improved geohazard monitoring and subsurface characterization. This study discussed the subsurface evaluation enabled by the DAS-vertical seismic profiling (DAS-VSP) system in the East China Sea, China. Taking advantage of a continuous recording of the vibrational wavefields through the strain deformation of the fiber deployed along the borehole, DAS-VSP is considered an emerging and promising alternative borehole acquisition method. It provides a wider-spectrum range of recordings in a cheaper, denser, and more continuous tense compared to conventional geophones. We explored the 3D DAS-VSP signal processing and imaging strategy. Based on a set of 3D DAS-VSP data of exceptional quality, this study reviewed the common processing challenges and practical solutions for de-noise, de-ghosting, de-multiple, and wavefield separation arising from the DAS acquisition mechanism in the offshore VSP scenario. High-quality down-going multiples were separated and imaged using the pre-stack Gaussian beam depth migration in addition to the primary reflection wavefield, providing significant additional illumination to support the subsurface evaluation. The current results validated the efficiency of the DAS-VSP survey and encouraged better geological interpretation.

Adaptive Focused Beam Prestack Depth Migration Under the Condition of Rugged Topography

The complex topography is a challenging issue in onshore seismic exploration. The rugged terrain and lateral change of near-surface velocity pose a significant obstacle to the accurate imaging of seismic data. The adaptive focused beam migration method retains the good applicability of the ray methods for calculating the seismic wavefield under complex surface conditions. It can effectively solve the contradiction between the imaging accuracy of deep and shallow strata in traditional Gaussian beam migration. We extend the adaptive focused beam migration approach to the deep domain imaging of seismic data under complex surface conditions. Firstly, the basic principles of the adaptive focused beam are reviewed. Then, Green’s functions of the seismic source and the receiving point of rugged topography are characterized by the adaptive focused beam, and an adaptive focused beam prestack depth migration method based on cross-correlation imaging is proposed. The full-wave-arrival imaging strategy is applied to image all wave arrivals of subsurface imaging points. A single input seismic trace is adopted for imaging, which can directly emit the focused beam from the receiving point of rugged topography for wave field continuation, thus avoiding multiple focusing. As a result, the applicability of the migration approach to complex surface conditions was improved, and the imaging accuracy of the migration method was also effectively enhanced. The numerical model migration test of different rugged topography conditions and tectonic forms verified that the proposed method was an effective prestack depth migration applicable to accurately imaging seismic data under the rugged topography condition.

From FWI to ultra-high-resolution imaging

In recent years, the development of time-lag full-waveform inversion (FWI) has enabled the use of the full wavefield (primary reflections, diving waves, and their multiples and ghosts) in the inversion process. With this advancement, it is possible to obtain a very detailed velocity model, ultimately reaching the point of deriving from the velocity a migration-like reflectivity image called the FWI image. When the FWI maximum frequency is increased, high-resolution velocity models are obtained, revealing superior reservoir information compared to conventional imaging results. Two case studies will be discussed in this paper. The first is in the Greater Castberg area where the 150 Hz FWI image greatly surpassed the Q Kirchhoff prestack depth migration image from the water-bottom level down to the reservoir (located at a depth of about 1.5 km). The second case study is over the Nordkapp Basin. The use of the full wavefield for the shallow ultra-high-resolution (UHR) FWI image (run at 200 Hz) revealed reverse faulting and pockmark details that were invisible with Kirchhoff prestack depth migration and reverse time migration. By using additional information present in multiples, ghosts, and diving waves, a spatial resolution down to 3 m was achieved. This made it possible to image very thin features without the need for a dedicated high-resolution acquisition design. The current UHR FWI image workflow provides velocity and reflectivity information in the near surface that is important in identifying optimal locations for various purposes such as well placement and wind-turbine installation.

Superior resolution through multiparameter FWI imaging: A new philosophy in seismic processing and imaging

Seismic processing and imaging workflows have been refined over many decades to attenuate aspects of the recorded wavefield which would be improperly mapped into the image domain by legacy migration algorithms such as Kirchhoff prestack depth migration. These workflows, which include techniques such as deghosting, designature, demultiple, and regularization, have become increasingly complex and time-consuming due to the sequential fashion in which they must be tested and applied. The single-scattering (primary-only) preprocessed data are then migrated and used in extensive model building workflows, including reflection residual moveout tomography, to refine low-frequency subsurface models. Obtaining optimal results at each stage requires subjective assessment of a wide range of parameter tests. Results can be highly variable, with different decisions resulting in very different outcomes. Such workflows mean that projects may take many months or even years. Full-waveform inversion (FWI) imaging offers an alternative philosophy to this conventional approach. FWI imaging is a least-squares multiscattering algorithm that uses the raw field data (transmitted and reflected arrivals as well as their multiples and ghosts) to determine many different subsurface parameters, including reflectivity. Because this approach uses the full wavefield, the subsurface is sampled more completely during the inversion. Here, we demonstrate the application of a novel multiparameter FWI imaging technique to generate high-resolution amplitude variation with angle reflectivity simultaneously with other model parameters, such as velocity and anisotropy, directly from the raw field data. Given that these results are obtained faster than the conventional workflow with a higher resolution, improved illumination, and reduced noise, we highlight the potential of multiparameter FWI imaging to supersede the conventional workflow.

Hydrocarbon accumulation conditions and exploration targets of the Sinian Dengying Formation in the Southeastern Sichuan Basin

No breakthrough has been made in petroleum exploration of the Sinian Dengying Formation in the southeastern Sichuan Basin as the major factors controlling the hydrocarbon accumulation of the Dengying Formation have not been well established, including the development period and distribution of the platform margin, the effectiveness of the source-reservoir combination and the reliability of the Cambrian pre-salt structural trap. Based on the interpretation of seismic data and hydrocarbon-source correlation, the distribution of the platform margin of the Sinian Dengying Formation in the southeastern Sichuan Basin was mapped. Moreover, the gas source rocks of the Dengying Formation were clarified according to geochemical characteristics of fluid inclusions, and the pre-salt structures were further confirmed by seismic forward modeling and experiments of prestack depth migration. Results showed that the sedimentary facies of the 4th member of Dengying Formation in the southeastern Sichuan Basin was platform margin of gentle slope type. The platform margin reservoir of the Dengying Formation and the Lower Cambrian source rocks constitute a favorable source - reservoir combination. The carbon isotope characteristics of gas in fluid inclusions indicate that some of the oil/gas in the Dengying Formation were contributed by the Doushantuo Formation source rocks. As the difference in seismic velocities between carbonate rocks and gypsum/salt beds may cause a reduction in the amplitude of pre-salt structures, the low-amplitude structures identified on current seismic profiles are underestimated instead of overestimated. Accordingly, the Sinian structural traps are confirmed to develop in the deep pre-salt sequences in the low-steep anticline belt, southeastern Sichuan Basin. On the basis of the above new understanding, we propose two types of hydrocarbon migration modes for the Sinian Dengying Formation in southeastern Sichuan Basin, including lateral migration of hydrocarbons generated from the Cambrian Qiongzusi Formation and vertical migration of hydrocarbons derived from the Sinian Doushantuo Formation. Lithologic traps at the platform margin and pre-salt structural traps within the platform are favorable targets for natural gas exploration.

Depth-domain angle and depth variant seismic wavelets extraction for prestack seismic inversion

Many prestack depth migration methods have been developed and widely used to generate depth-domain seismic images, resulting in a need for depth-domain prestack seismic inversion of the subsurface elastic properties for reservoir characterization. Time-domain inversion techniques often are used after the depth-domain data set is transformed to the time domain. We provide a new technique for directly inverting depth-domain prestack seismic data for subsurface elastic properties: P-impedance ([Formula: see text]), S-impedance ([Formula: see text]), and density ([Formula: see text]) in the depth domain. The proposed depth-domain workflow eliminates the need for time-depth/depth-time conversion, making it efficient and effective. Using a depth-wavenumber decomposition approach, the suggested workflow first extracts a collection of depth and angle varying wavelets to handle the potential nonstationarity of the depth-domain prestack seismic data. The extracted depth and angle variant wavelets are used in a basis pursuit inversion, which enhances the resolution of the inversion results. The workflow is tested on the Wenan 3D field data set, demonstrating its viability for practical applications.

Basin-scale estimates of thermal stress and expelled petroleum from Mesozoic-Cenozoic potential source rocks, southern Gulf of Mexico

The Campeche and Yucatan salt basins remain two of the least explored areas of the Gulf of Mexico basin. This study uses a grid of 23,600 line-km pre-stack depth migrated (PSDM) 2D seismic reflection profiles, shipborne gravity data, and open-source geologic information to model the thermal stress of four potential source intervals (Oxfordian, Tithonian-centered, Cenomanian-Turonian, lower Miocene). We performed map-based, and 1D thermal modeling along two margin-perpendicular transects, each consisting of five pseudo-wells tied to the regional grid of seismic reflection data. Our modeling takes into account thermal stress variations related to the depth of base lithosphere, crustal type and thickness, paleo-water depth, Jurassic salt thickness, and the transient heat flow effects related to recent clastic sedimentation. We predict that deeply-buried, salt-related minibasins along the outer marginal trough are mature for petroleum expulsion with deeply-buried Mesozoic source rocks within the oil window from late Paleogene to early Neogene time. The ‘lag time’ required for vertical oil migration explains why oil maturation occurred during the late Paleogene to early Neogene but active oil seeps are observed today at the sea surface. We predict that oil is present in subsurface traps in the deepest part of the outer marginal trough and we calculate that the Oxfordian source interval has expelled a cumulative 20 million bbl of oil equivalent [BOE]/km2 and that the Tithonian-centered source interval has expelled 67 million bbl of oil equivalent [BOE]/km2.

Seismic Survey in Lesser Himalayan Thrust Belt, Western Nepal

Two hundred km of 2D seismic survey was carried out at the Lesser Himalayan Thrust Belts in Dailekh district, western Nepal. The main motivation is to elucidate the geologic relationship between the known oil and gas seeps, subsurface structure, and stratigraphy in the area. This is a challenging task which is from its extreme structural and geological complexity such as thrust faulting, tight folding, steep dip layers, and strong lateral variations in seismic velocity. Seismic data were acquired with SERCEL 428XL system and processed by GEOEAST computer software. In order to increase signal-to-noise ratio (SNR), suppress interference, and search for optimum acquisition parameters, a series of comparative tests on the different charge depth and size, group interval, CDP fold, geophone array, and single high-sensitivity geophone were conducted. We also tested 2S3L (two lines shooting and three lines receiving) wide line profiling. The results indicate that single hole with charge depth of 12 m, 4-16 kg charge size (less charge size for the densely populated areas), single high-sensitivity geophone, and 1S2L wide line profiling with 132 folds are the optimum acquisition parameters. On the basis of comparative process experiment, data processing workflow consisting of data preparation, prestack denoising, amplitude compensation, deconvolution, tomography static correction, velocity analysis, residual static correction, CRS stack, poststack migration, prestack time migration (PSTM), and prestack depth migration (PSDM) was selected. Maybe affected by problem of conflicting dip in complex media, CRS stack section does not show satisfactory geological characteristics. PSTM profile has moderate signal-to-noise (S/N) ratio; the shallow, medium, and deep continuous reflections can be observed in section. More details of the geological structures can be observed in PSDM section, especially in medium and shallow layers (less than 3000 ms or 4000 m), but PSDM method is more expensive and highly time consuming than that of CRS stack and PSTM. So, the PSTM section can be reasonably used for geological interpretation. By reference to field mapping, thrust characteristics, and MT data, the final interpretation to the PSTM section identified the interfaces of 6 geological units (Paleoproterozoic Nabhisthan Fm., Paleoproterozoic Dubidanda Fm., Neogene to Late Cretaceous Surkhet group, Late Carboneferous to Early Cretaeous Gondwana group, Mesoproterozoic Upper Lakharpata group, and Lower Lakharpata group) and delineated Main Boundary Thrust (MBT), Ramgarh Thrust (RMT), Padukasthan Thrust (PT), and Dailekh Thrust (DT). The bottom of Surkhet group which is our top target zone is about 4250 meters deep.

Elastic anisotropic finite-difference full-wave modeling and imaging of 2D tilted transversely isotropic (TTI) media

Imaging below the tilted transversely isotropic (TTI) media pose serious problems and distortions of the subsurface geological targets of interest for hydrocarbon exploration, which is very common in fold-thrust belts or active tectonic areas like the Rocky Mountain foothills of the Canada or foothills of the Himalaya. To model and image below the TTI media having symmetry axis orthogonal to dipping anisotropic layers, a 2D TTI thrust sheet is considered for elastic anisotropic finite-difference full-wave modeling and generation of synthetic seismic data using the robust staggered-grid scheme of wave propagation. Free-surface boundary conditions at the top and absorbing boundary conditions for the other three sides of the model have been imposed to reduce undesirable edge effects and suppress dispersions. The synthetic seismic data generated for the 2D TTI thrust model are migrated using the anisotropic Kirchhoff pre-stack time (PSTM) and depth migration (PSDM) technique followed by the migration velocity analysis (MVA) algorithm. The traveltime computations are made using robust eikonal equations with fourth-order elastic anisotropic finite-difference solutions to handle large tilts (ν) along with anisotropic ray-tracing to image below the TTI thrust sheet. The MVA algorithm adopted also estimates the anisotropic parameters ε and δ accurately within the acceptable error limit through successive iterations of inversion and parameter update, which results in sufficient flattening of the reflection events in the common image gathers (CIGs). This is a powerful and stable methodology to handle large tilt (ν = 60°) of the steeply dipping reflectors in complex geological structure like the TTI thrust model without much distortion of CIGs for which the symmetry axis of each TTI block is set orthogonal to the bottom of the reflectors having different tilts.

Seismic reflection waveform inversion based on Gauss–Newton optimization

Abstract Owing to the lack of long-offset signals and effective low-frequency components, it is extremely difficult for classical full-waveform inversion (FWI) to obtain the long-wavelength components of medium and deep elastic parameter models. Considering ray theory, reflection traveltime tomography is restricted by resolution, and the reflection waveform inversion (RWI) method of the wave equation, and its practicability, has received greater attention. However, previous studies into RWI have primarily used first-order optimization methods, in which convergence and resolution still need substantial improvement. Based on the second-order optimization inversion theory, we deduced the reflection wave-sensitive kernel, objective function gradient and Hessian operator of background and perturbation model. Additionally, we revealed the de-fuzzification effect of the Hessian matrix on the gradient of inversion and the working mechanism of background velocity updating on direction optimization and inversion resolution improvement. Experiments on the SEG/EAGE overthrust model showed that an RWI with the Gauss–Newton method using an approximate Hessian matrix significantly improved convergence and wide-spectrum modeling capability compared with the conjugate gradient method. On the streamer seismic data of the East China Sea, the second-order optimization RWI method surpassed the commonly used reflection traveltime tomography based on pre-stack depth migration. In addition, the RTM of inversion model can depict the complex fault system in the Changjiang sag with a high resolution, which improves the imaging quality of the deep basement.