Multiwave inversion: A key step for depth model building — Examples from the Sultanate of Oman

2021 ◽  
Vol 40 (8) ◽  
pp. 610-618
Author(s):  
Daniela Donno ◽  
Mohammad Sheryar Farooqui ◽  
Mostafa Khalil ◽  
David McCarthy ◽  
Didier Solyga ◽  
...  

The near surface in the Middle East, particularly in the Sultanate of Oman, is characterized by very shallow high-velocity carbonates and anhydrites interleaved by slow-velocity clastic layers, resulting in sharp velocity inversions in the first few hundred meters below the surface. In addition, the surface is characterized by features such as unconsolidated materials within dry riverbeds (known as “wadis”), small jebels, and sand dunes, which cause distortions in the underlying shallow and deeper seismic images. This work presents the building of a near-surface model by using multiwave inversion that jointly inverts information from P-wave first breaks and surface-wave dispersion curves. The use of surface waves in combination with first breaks captures the lateral and vertical velocity variations, especially in the shallowest parts of the near surface. This paper focuses on the analysis of two drawbacks of this technology: the picking of the input data information, which can be cumbersome and time consuming, and the limited penetration depth of surface waves at the typical frequencies of active data. To overcome these issues, an innovative workflow is proposed that combines the use of an unsupervised machine learning technique to guide the pick extraction phase and the reconstruction of ultra-low-frequency surface waves (0.5 to 1.5 Hz) through an interferometry process using information from natural and ambient noise. Deeper near-surface P- and S-wave velocity models can be obtained with multiwave inversion using these ultra-low frequencies. The integration of a near-surface model into the velocity model building workflow brings a major improvement in depth imaging from shallow to deep structures, as demonstrated on two data sets from the Sultanate of Oman.

2019 ◽  
Vol 218 (3) ◽  
pp. 1873-1891 ◽  
Author(s):  
Farbod Khosro Anjom ◽  
Daniela Teodor ◽  
Cesare Comina ◽  
Romain Brossier ◽  
Jean Virieux ◽  
...  

SUMMARY The analysis of surface wave dispersion curves (DCs) is widely used for near-surface S-wave velocity (VS) reconstruction. However, a comprehensive characterization of the near-surface requires also the estimation of P-wave velocity (VP). We focus on the estimation of both VS and VP models from surface waves using a direct data transform approach. We estimate a relationship between the wavelength of the fundamental mode of surface waves and the investigation depth and we use it to directly transform the DCs into VS and VP models in laterally varying sites. We apply the workflow to a real data set acquired on a known test site. The accuracy of such reconstruction is validated by a waveform comparison between field data and synthetic data obtained by performing elastic numerical simulations on the estimated VP and VS models. The uncertainties on the estimated velocity models are also computed.


2018 ◽  
Vol 6 (4) ◽  
pp. SM27-SM37 ◽  
Author(s):  
Jing Li ◽  
Kai Lu ◽  
Sherif Hanafy ◽  
Gerard Schuster

Two robust imaging technologies are reviewed that provide subsurface geologic information in challenging environments. The first one is wave-equation dispersion (WD) inversion of surface waves and guided waves (GW) for the shear-velocity (S-wave) and compressional-velocity (P-wave) models, respectively. The other method is traveltime inversion for the velocity model, in which supervirtual refraction interferometry (SVI) is used to enhance the signal-to-noise ratio of far-offset refractions. We have determined the benefits and liabilities of both methods with synthetic seismograms and field data. The benefits of WD are that (1) there is no layered-medium assumption, as there is in conventional inversion of dispersion curves. This means that 2D or 3D velocity models can be accurately estimated from data recorded by seismic surveys over rugged topography, and (2) WD mostly avoids getting stuck in local minima. The liability is that WD for surface waves is almost as expensive as full-waveform inversion (FWI) and, for Rayleigh waves, only recovers the S-velocity distribution to a depth no deeper than approximately 1/2 to 1/3 wavelength of the lowest-frequency surface wave. The limitation for GW is that, for now, it can estimate the P-velocity model by inverting the dispersion curves from GW propagating in near-surface low-velocity zones. Also, WD often requires user intervention to pick reliable dispersion curves. For SVI, the offset of usable refractions can be more than doubled, so that traveltime tomography can be used to estimate a much deeper model of the P-velocity distribution. This can provide a more effective starting velocity model for FWI. The liability is that SVI assumes head-wave first arrivals, not those from strong diving waves.


Geophysics ◽  
1951 ◽  
Vol 16 (1) ◽  
pp. 63-80 ◽  
Author(s):  
Milton B. Dobrin

A non‐mathematical summary is presented of the published theories and observations on dispersion, i.e., variation of velocity with frequency, in surface waves from earthquakes and in waterborne waves from shallow‐water explosions. Two further instances are cited in which dispersion theory has been used in analyzing seismic data. In the seismic refraction survey of Bikini Atoll, information on the first 400 feet of sediments below the lagoon bottom could not be obtained from ground wave first arrival times because shot‐detector distances were too great. Dispersion in the water waves, however, gave data on speed variations in the bottom sediments which made possible inferences on the recent geological history of the atoll. Recent systematic observations on ground roll from explosions in shot holes have shown dispersion in the surface waves which is similar in many ways to that observed in Rayleigh waves from distant earthquakes. Classical wave theory attributes Rayleigh wave dispersion to the modification of the waves by a surface layer. In the case of earthquakes, this layer is the earth’s crust. In the case of waves from shot‐holes, it is the low‐speed weathered zone. A comparison of observed ground roll dispersion with theory shows qualitative agreement, but it brings out discrepancies attributable to the fact that neither the theory for liquids nor for conventional solids applies exactly to unconsolidated near‐surface rocks. Additional experimental and theoretical study of this type of surface wave dispersion may provide useful information on the properties of the surface zone and add to our knowledge of the mechanism by which ground roll is generated in seismic shooting.


1996 ◽  
Vol 86 (6) ◽  
pp. 1704-1713 ◽  
Author(s):  
R. D. Catchings ◽  
W. H. K. Lee

Abstract The 17 January 1994, Northridge, California, earthquake produced strong ground shaking at the Cedar Hills Nursery (referred to here as the Tarzana site) within the city of Tarzana, California, approximately 6 km from the epicenter of the mainshock. Although the Tarzana site is on a hill and is a rock site, accelerations of approximately 1.78 g horizontally and 1.2 g vertically at the Tarzana site are among the highest ever instrumentally recorded for an earthquake. To investigate possible site effects at the Tarzana site, we used explosive-source seismic refraction data to determine the shallow (<70 m) P-and S-wave velocity structure. Our seismic velocity models for the Tarzana site indicate that the local velocity structure may have contributed significantly to the observed shaking. P-wave velocities range from 0.9 to 1.65 km/sec, and S-wave velocities range from 0.20 and 0.6 km/sec for the upper 70 m. We also found evidence for a local S-wave low-velocity zone (LVZ) beneath the top of the hill. The LVZ underlies a CDMG strong-motion recording site at depths between 25 and 60 m below ground surface (BGS). Our velocity model is consistent with the near-surface (<30 m) P- and S-wave velocities and Poisson's ratios measured in a nearby (<30 m) borehole. High Poisson's ratios (0.477 to 0.494) and S-wave attenuation within the LVZ suggest that the LVZ may be composed of highly saturated shales of the Modelo Formation. Because the lateral dimensions of the LVZ approximately correspond to the areas of strongest shaking, we suggest that the highly saturated zone may have contributed to localized strong shaking. Rock sites are generally considered to be ideal locations for site response in urban areas; however, localized, highly saturated rock sites may be a hazard in urban areas that requires further investigation.


2018 ◽  
Vol 58 (2) ◽  
pp. 884
Author(s):  
Lianping Zhang ◽  
Haryo Trihutomo ◽  
Yuelian Gong ◽  
Bee Jik Lim ◽  
Alexander Karvelas

The Schlumberger Multiclient Exmouth 3D survey was acquired over the Exmouth sub-basin, North West Shelf Australia and covers 12 600 km2. One of the primary objectives of this survey was to produce a wide coverage of high quality imaging with advanced processing technology within an agreed turnaround time. The complexity of the overburden was one of the imaging challenges that impacted the structuration and image quality at the reservoir level. Unlike traditional full-waveform inversion (FWI) workflow, here, FWI was introduced early in the workflow in parallel with acquisition and preprocessing to produce a reliable near surface velocity model from a smooth starting model. FWI derived an accurate and detailed near surface model, which subsequently benefitted the common image point (CIP) tomography model updates through to the deeper intervals. The objective was to complete the FWI model update for the overburden concurrently with the demultiple stages hence reflection time CIP tomography could start with a reasonably good velocity model upon completion of the demultiple process.


Author(s):  
Gleb S. Chernyshov ◽  
◽  
Anton A. Duchkov ◽  
Aleksander A. Nikitin ◽  
Ivan Yu. Kulakov ◽  
...  

The problem of tomographic inversion is non–unique and requires regularization to solve it in a stable manner. It is highly non–trivial to choose between various regularization approaches or tune the regularization parameters themselves. We study the influence of one particular regularization parameter on the resolution and accuracy the tomographic inversion for the near–surface model building. We propose another regularization parameter, which allows to increase the accuracy of model building.


Geophysics ◽  
2012 ◽  
Vol 77 (1) ◽  
pp. V11-V20 ◽  
Author(s):  
Roohollah Askari ◽  
Robert J. Ferguson

Wavenumber, group velocity, phase velocity, and frequency-dependent attenuation characterize the propagation of surface waves in dispersive, attenuating media. We use a mathematical model based on the generalized [Formula: see text] transform to simultaneously estimate these characteristic parameters for later use in joint inversion for near-surface shear wave velocity. We use a scaling factor in the generalized S transform to enable the application of the method in a highly dispersive medium. We introduce a cost function in the [Formula: see text]-domain to estimate an optimum value for the scaling factor. We also use the cost function to generalize the application of the method for noisy data, especially data with a low signal-to-noise ratio at low frequencies. In that case, the estimated wavenumber is perturbed. As a solution, we estimate wavenumber perturbation by minimizing the cost function, using Simulated Annealing. We use synthetic and real data to show the efficiency of the method for the estimation of the propagation parameters of highly dispersive and noisy media.


Geophysics ◽  
2005 ◽  
Vol 70 (2) ◽  
pp. V21-V29 ◽  
Author(s):  
Xander H. Campman ◽  
Kasper van Wijk ◽  
John A. Scales ◽  
Gérard C. Herman

When traveling through a complex overburden, upcoming seismic body waves can be disturbed by scattering from local heterogeneities. Currently, surface-consistent static and amplitude corrections correct for rapid variations in arrival times and amplitudes of a reflector, but these methods impose strong assumptions on the near-surface model. Observations on synthetic and laboratory experiments of near-surface scattering with densely sampled data suggest that removing noise from near-receiver scattering requires multichannel approaches rather than single-channel, near-surface corrections. In this paper we develop a wavefield-based imaging method to suppress surface waves scattered directly beneath the receivers. Using an integral-equation formulation, we account for near-surface heterogeneities by a surface impedance function. This impedance function is used to model scattered surface waves, excited by upcoming wavefronts. The final step in our algorithm is to subtract the scattered surface waves. We successfully apply this method to laboratory data of scattered surface waves, excited and monitored with a noncontacting acquisition system.


Sign in / Sign up

Export Citation Format

Share Document