Power-gradient velocity model

Geophysics ◽  
2009 ◽  
Vol 74 (5) ◽  
pp. U13-U33 ◽  
Author(s):  
Alexey Stovas

A proposed power-gradient velocity model incorporates several well-known velocity models as special cases. The model covers a wide range of possible velocity distributions and has four parameters, giving more flexibility in velocity-model manipulation. For this nonlinear model, the kinematic characteristics — offset-traveltime parameteric equations, traveltime parameters, relative geometric spreading, and propagator phase — are computed. The characteristics are investigated with respect to a parameter responsible for nonlinearity of velocity distribution. The inversion of travel-time parameters was studied in three- and four-parameter frameworks.

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.


1990 ◽  
Vol 80 (2) ◽  
pp. 395-410 ◽  
Author(s):  
Glenn D. Nelson ◽  
John E. Vidale

Abstract We present a new method for locating earthquakes in a region with arbitrarily complex three-dimensional velocity structure, called QUAKE3D. Our method searches a gridded volume and finds the global minimum travel-time residual location within the volume. Any minimization criterion may be employed. The L1 criterion, which minimizes the sum of the absolute values of travel-time residuals, is especially useful when the station coverage is sparse and is more robust than the L2 criterion (which minimizes the RMS sum) employed by most earthquake location programs. On a UNIX workstation with 8 Mbytes memory, travel-time grids of size 150 by 150 by 50 are reasonably employed, with the actual geographic coverage dependent on the grid spacing. Location precision is finer than the grid spacing. Earthquake recordings at six stations in Bear Valley are located as an example, using various layered and laterally varying velocity models. Locations with QUAKE3D are nearly identical to HYPOINVERSE locations when the same flat-layered velocity model is used. For the examples presented, the computation time per event is approximately 4 times slower than HYPOINVERSE, but the computation time for QUAKE3D is dependent only on the grid size and number of stations, and independent of the velocity model complexity. Using QUAKE3D with a laterally varying velocity model results in locations that are physically more plausible and statistically more precise. Compared to flat-layered solutions, the earthquakes are more closely aligned with the surface fault trace, are more uniform in depth distribution, and the event and station travel-time residuals are much smaller. Hypocentral error bars computed by QUAKE3D are more realistic in that the trade-off of depth versus origin time is implicit in our error estimation, but ignored by HYPOINVERSE.


Solid Earth ◽  
2018 ◽  
Vol 9 (6) ◽  
pp. 1487-1505 ◽  
Author(s):  
Claudia Werner ◽  
Erik H. Saenger

Abstract. Time reverse imaging (TRI) is evolving into a standard technique for locating and characterising seismic events. In recent years, TRI has been employed for a wide range of applications from the lab scale, to the field scale and up to the global scale. No identification of events or their onset times is necessary when locating events with TRI; therefore, it is especially suited for locating quasi-simultaneous events and events with a low signal-to-noise ratio. However, in contrast to more regularly applied localisation methods, the prerequisites for applying TRI are not sufficiently known.To investigate the significance of station distributions, complex velocity models and signal-to-noise ratios with respect to location accuracy, numerous simulations were performed using a finite difference code to propagate elastic waves through three-dimensional models. Synthetic seismograms were reversed in time and reinserted into the model. The time-reversed wave field back propagates through the model and, in theory, focuses at the source location. This focusing was visualised using imaging conditions. Additionally, artificial focusing spots were removed using an illumination map specific to the set-up. Successful locations were sorted into four categories depending on their reliability. Consequently, individual simulation set-ups could be evaluated by their ability to produce reliable source locations.Optimal inter-station distances, minimum apertures, relations between the array and source locations, heterogeneities of inter-station distances and the total number of stations were investigated for different source depths and source types. Additionally, the accuracy of the locations was analysed when using a complex velocity model or a low signal-to-noise ratio.Finally, an array in southern California was investigated regarding its ability to locate seismic events at specific target depths while using the actual velocity model for that region. In addition, the success rate with recorded data was estimated.Knowledge about the prerequisites for using TRI enables the estimation of success rates for a given problem. Furthermore, it reduces the time needed to adjust stations to achieve more reliable locations and provides a foundation for designing arrays for applying TRI.


2020 ◽  
Vol 10 (12) ◽  
pp. 4391
Author(s):  
Yasir Bashir ◽  
Nordiana Mohd Muztaza ◽  
Seyed Yaser Moussavi Alashloo ◽  
Syed Haroon Ali ◽  
Deva Prasad Ghosh

Fractured imaging is an important target for oil and gas exploration, as these images are heterogeneous and have contain low-impedance contrast, which indicate the complexity in a geological structure. These small-scale discontinuities, such as fractures and faults, present themselves in seismic data in the form of diffracted waves. Generally, seismic data contain both reflected and diffracted events because of the physical phenomena in the subsurface and due to the recording system. Seismic diffractions are produced once the acoustic impedance contrast appears, including faults, fractures, channels, rough edges of structures, and karst sections. In this study, a double square root (DSR) equation is used for modeling of the diffraction hyperbola with different velocities and depths of point diffraction to elaborate the diffraction hyperbolic pattern. Further, we study the diffraction separation methods and the effects of the velocity analysis methods (semblance vs. hybrid travel time) for velocity model building for imaging. As a proof of concept, we apply our research work on a steep dipping fault model, which demonstrates the possibility of separating seismic diffractions using dip frequency filtering (DFF) in the frequency–wavenumber (F-K) domain. The imaging is performed using two different velocity models, namely the semblance and hybrid travel time (HTT) analysis methods. The HTT method provides the optimum results for imaging of complex structures and imaging below shadow zones.


Author(s):  
Ying Liu ◽  
Huajian Yao ◽  
Haijiang Zhang ◽  
Hongjian Fang

Abstract Southwest China, located at the southeastern margin of the Tibetan plateau, plays an important role for the plateau growth and its material extrusion. It has complicated tectonic environment and strong seismic activities including the 2008 Wenchuan great earthquake. Numerous geophysical studies have been conducted in southwest China. However, a community velocity model (CVM) in this region is still not available, which makes it difficult to have a consistent catalog of earthquake locations and focal mechanisms and a consistent velocity model for simulating strong ground motions and evaluating earthquake hazards. In this study, we aim at building a high-resolution CVM (both VP and VS) of the crust and uppermost mantle in southwest China along with earthquake locations by joint inversion of body- and surface-wave travel-time data. In total, we have assembled 386,958 P- and 372,662 S-wave first arrival times and nearly 8100 Rayleigh-wave dispersion curves in the period band of 5–50 s. A multigrid strategy is adopted in the joint inversion. A coarser horizontal grid interval of 0.5° is first used and then a finer grid interval of 0.25° is used with initial models interpolated from the coarser-grid inverted velocity models. The spatial resolution of both VP and VS models can reach up to 0.5° horizontally and 10 km vertically according to the checkerboard tests. The comparisons of our inverted VP and VS models with those from other studies show general consistency in large-scale features. The inverted models are further validated by P-wave arrival times from active sources and Rayleigh-wave data. In general, our velocity models show two low-velocity zones in the middle-lower crust and a prominent high-velocity region in between them. Our new models have been served as the first version of the CVM in southwest China (SWChinaCVM-1.0) for future studies.


Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE291-VE302 ◽  
Author(s):  
Stefan Dümmong ◽  
Kristina Meier ◽  
Dirk Gajewski ◽  
Christian Hübscher

Velocity-model determination during seismic data processing is crucial for any kind of depth imaging. We compared two approaches of grid tomography: prestack stereotomography and normal-incidence-point (NIP) wave tomography. Whereas NIP wave tomography is based on wavefield attributes obtained during the common reflection surface stack and thus on the underlying hyperbolic second-order traveltime approximation, prestack stereotomography describes traveltimes by local slopes (i.e., linearly) in the prestack data domain. To analyze the impact of the different traveltime approximations and the different input-data domains on velocity model building, we applied two implementations of these techniques to two profiles of a field marine data set from the Levante Basin, eastern Mediterranean. Because ofthe presence of a thick, tabular mobile unit of the Messinian evaporites, strong vertical and lateral velocity contrasts had been expected. The velocity models revealed the reconstruction of high-velocity contrasts by grid tomographic methods is limited because of the smooth description of the velocity distribution. The lateral resolution of velocities obtained from prestack stereotomography appears to be better than those from NIP wave tomography, which is related to the difference in the approximation of traveltimes, the determination of input data, and the description of the velocity distribution. Other differences are caused mainly by different implementations of the inversion schemes. Nevertheless, both algorithms provide suitable models for high-quality depth imaging, whereas most of the reflections are fairly flat in CIGs.


Author(s):  
Hasbi Ash Shiddiqi ◽  
Pa Pa Tun ◽  
Lars Ottemöller

ABSTRACT Earthquake monitoring in Myanmar has improved in recent years because of an increased number of seismic stations. This provides a good quality dataset to derive a minimum 1D velocity model and local magnitude (ML) scale for the Myanmar region, which will improve the earthquake location and magnitude estimates in this region. We combined and reprocessed earthquake catalogs from the Department of Meteorology and Hydrology of Myanmar and the International Seismological Centre. Additional waveform data from various sources were processed as well. A total of 419 earthquakes were selected based on azimuthal gap, minimum number of stations, and root mean square travel‐time residuals. A set of initial seismic velocity models was derived from various seismic velocity models. These models were randomly perturbed and used as initial models in a coupled hypocenter and 1D seismic velocity inversion procedure. We compared the average mean travel‐time residuals from the initial and inverted models. The best final model showed an improvement of location standard errors compared to the old model. Furthermore, the local magnitude scale inversion for the Myanmar region was performed using 194 earthquakes having a minimum of two amplitude observations. The following ML scale was obtained ML=logA(nm)+1.485×logR(km)+0.00118×R(km)−2.77+S. This scale is valid for hypocentral distance up to 1000 km and magnitudes up to ML 6.2.


2018 ◽  
Author(s):  
Claudia Werner ◽  
Erik H. Saenger

Abstract. Time Reverse Imaging (TRI) is evolving into a standard technique for localizing and characterizing seismic events. In recent years, TRI has been applied to a wide range of applications from the lab scale over the field scale up to the global scale. No identification of events and their onset times is necessary when localizing events with TRI. Therefore, it is especially suited for localizing quasi-simultaneous events and events with a low signal-to-noise ratio. However, in contrast to more regularly applied localization methods, the prerequisites for applying TRI are not sufficiently known. To investigate the significance of station distributions, complex velocity models and signal-to-noise ratios for the localization quality, numerous simulations were performed using a finite difference code to propagate elastic waves through three-dimensional models. Synthetic seismograms were reversed in time and re-inserted into the model. The time-reversed wavefield backpropagates through the model and, in theory, focuses at the source location. This focusing was visualized using imaging conditions. Additionally, artificial focusing spots were removed with an illumination map specific to the setup. Successful localizations were sorted into four categories depending on their reliability. Consequently, individual simulation setups could be evaluated by their ability to produce reliable localizations. Optimal inter-station distances, minimum apertures, relations between array and source location, heterogeneities of inter-station distances and total number of stations were investigated for different source depth as well as source types. Additionally, the quality of the localization was analysed when using a complex velocity model or a low signal-to-noise ratio. Finally, an array in Southern California was investigated for its ability to localize seismic events in specific target depths while using the actual velocity model for that region. In addition, the success rate with recorded data was estimated. Knowledge about the prerequisites for using TRI enables the estimation of success rates for a given problem. Furthermore, it reduces the time needed for adjusting stations to achieve more reliable localizations and provides a foundation for designing arrays for applying TRI.


2021 ◽  
Vol 225 (2) ◽  
pp. 1020-1031
Author(s):  
Huachen Yang ◽  
Jianzhong Zhang ◽  
Kai Ren ◽  
Changbo Wang

SUMMARY A non-iterative first-arrival traveltime inversion method (NFTI) is proposed for building smooth velocity models using seismic diving waves observed on irregular surface. The new ray and traveltime equations of diving waves propagating in smooth media with undulant observation surface are deduced. According to the proposed ray and traveltime equations, an analytical formula for determining the location of the diving-wave turning points is then derived. Taking the influence of rough topography on first-arrival traveltimes into account, the new equations for calculating the velocities at turning points are established. Based on these equations, a method is proposed to construct subsurface velocity models from the observation surface downward to the bottom using the first-arrival traveltimes in common offset gathers. Tests on smooth velocity models with rugged topography verify the validity of the established equations, and the superiority of the proposed NFTI. The limitation of the proposed method is shown by an abruptly-varying velocity model example. Finally, the NFTI is applied to solve the static correction problem of the field seismic data acquired in a mountain area in the western China. The results confirm the effectivity of the proposed NFTI.


2021 ◽  
Author(s):  
Francesca D’Ajello Caracciolo ◽  
Rodolfo Console

AbstractA set of four magnitude Ml ≥ 3.0 earthquakes including the magnitude Ml = 3.7 mainshock of the seismic sequence hitting the Lake Constance, Southern Germany, area in July–August 2019 was studied by means of bulletin and waveform data collected from 86 seismic stations of the Central Europe-Alpine region. The first single-event locations obtained using a uniform 1-D velocity model, and both fixed and free depths, showed residuals of the order of up ± 2.0 s, systematically affecting stations located in different areas of the study region. Namely, German stations to the northeast of the epicenters and French stations to the west exhibit negative residuals, while Italian stations located to the southeast are characterized by similarly large positive residuals. As a consequence, the epicentral coordinates were affected by a significant bias of the order of 4–5 km to the NNE. The locations were repeated applying a method that uses different velocity models for three groups of stations situated in different geological environments, obtaining more accurate locations. Moreover, the application of two methods of relative locations and joint hypocentral determination, without improving the absolute location of the master event, has shown that the sources of the four considered events are separated by distances of the order of one km both in horizontal coordinates and in depths. A particular attention has been paid to the geographical positions of the seismic stations used in the locations and their relationship with the known crustal features, such as the Moho depth and velocity anomalies in the studied region. Significant correlations between the observed travel time residuals and the crustal structure were obtained.


Sign in / Sign up

Export Citation Format

Share Document