Finite‐difference EM modeling and imaging in a medium with a vertical axis of symmetry

1999 ◽  
Author(s):  
Efthimios Tartaras ◽  
Michael Zhdanov
Keyword(s):  
Finite Difference ◽  
Vertical Axis ◽  
Axis Of Symmetry

scholarly journals Theoretical Gravity Anomalies of Glaciers Having Parabolic Cross-Sections

1964 ◽  
Vol 5 (38) ◽  
pp. 255-257 ◽  
Author(s):  
Charles E. Corbató
Keyword(s):  
Cross Sections ◽  
Lower Boundary ◽  
Gravity Anomalies ◽  
Vertical Axis ◽  
Two Dimensional ◽  
Axis Of Symmetry ◽  
Upper Boundary

AbstractEquations and a graph are presented for calculating gravity anomalies on a two-dimensional glacier model having a horizontal upper boundary and a lower boundary which is a parabola with a vertical axis of symmetry.

scholarly journals An efficient wavefield inversion for transversely isotropic media with a vertical axis of symmetry

Geophysics ◽  
2020 ◽  
Vol 85 (3) ◽  
pp. R195-R206 ◽  
Author(s):  
Chao Song ◽  
Tariq Alkhalifah
Keyword(s):  
High Resolution ◽  
Optimization Problem ◽  
Source Function ◽  
Transversely Isotropic ◽  
Vertical Axis ◽  
Trade Off ◽  
Transversely Isotropic Media ◽  
Highly Nonlinear ◽  
Isotropic Media ◽  
Axis Of Symmetry

Conventional full-waveform inversion (FWI) aims at retrieving a high-resolution velocity model directly from the wavefields measured at the sensor locations resulting in a highly nonlinear optimization problem. Due to the high nonlinearity of FWI (manifested in one form in the cycle-skipping problem), it is easy to fall into local minima. Considering that the earth is truly anisotropic, a multiparameter inversion imposes additional challenges in exacerbating the null-space problem and the parameter trade-off issue. We have formulated an optimization problem to reconstruct the wavefield in an efficient matter with background models by using an enhanced source function (which includes secondary sources) in combination with fitting the data. In this two-term optimization problem to fit the wavefield to the data and to the background wave equation, the inversion for the wavefield is linear. Because we keep the modeling operator stationary within each frequency, we only need one matrix inversion per frequency. The inversion for the anisotropic parameters is handled in a separate optimization using the wavefield and the enhanced source function. Because the velocity is the dominant parameter controlling the wave propagation, it is updated first. Thus, this reduces undesired updates for anisotropic parameters due to the velocity update leakage. We find the effectiveness of this approach in reducing parameter trade-off with a distinct Gaussian anomaly model. We find that in using the parameterization [Formula: see text], and [Formula: see text] to describe the transversely isotropic media with a vertical axis of symmetry model in the inversion, we end up with high resolution and minimal trade-off compared to conventional parameterizations for the anisotropic Marmousi model. Application on 2D real data also indicates the validity of our method.

scholarly journals Frequency domain multiparameter acoustic inversion for transversely isotropic media with a vertical axis of symmetry

Geophysics ◽  
2019 ◽  
Vol 84 (1) ◽  
pp. C1-C14 ◽  
Author(s):  
Ramzi Djebbi ◽  
Tariq Alkhalifah
Keyword(s):  
Frequency Domain ◽  
Transversely Isotropic ◽  
Vertical Axis ◽  
The Other ◽  
Sensitivity Analyses ◽  
Data Set ◽  
Trade Off ◽  
Transversely Isotropic Media ◽  
Isotropic Media ◽  
Axis Of Symmetry

Multiparameter full-waveform inversion for transversely isotropic media with a vertical axis of symmetry (VTI) suffers from the trade-off between the parameters. The trade-off results in the leakage of one parameter’s update into the other. It affects the accuracy and convergence of the inversion. The sensitivity analyses suggested a parameterization using the horizontal velocity [Formula: see text], Thomsen’s parameter [Formula: see text], and the anelliptic parameter [Formula: see text] to reduce the trade-off for surface recorded seismic data. We aim to invert for this parameterization using the scattering integral (SI) method. The available Born sensitivity kernels, within this approach, can be used to calculate additional inversion information. We mainly compute the diagonal of the approximate Hessian, used as a conjugate-gradient preconditioner, and the gradients’ step lengths. We consider modeling in the frequency domain. The large computational cost of the SI method can be avoided with direct Helmholtz equation solvers. We applied our method to the VTI Marmousi II model for various inversion strategies. We found that we can invert the [Formula: see text] accurately. For the [Formula: see text] parameter, only the short wavelengths are well-recovered. On the other hand, the [Formula: see text] parameter impact is weak on the inversion results and can be fixed. However, a good background [Formula: see text], with accurate long wavelengths, is needed to correctly invert for [Formula: see text]. Furthermore, we invert a real data set acquired by CGG from offshore Australia. We simultaneously invert all three parameters using our inversion approach. The velocity model is improved, and additional layers are recovered. We confirm the accuracy of the results by comparing them with well-log information, as well as looking at the data and angle gathers.

The effect of transverse isotropy on isotropic traveltime inversion of vertical seismic profile data—A modeling study

Geophysics ◽  
1990 ◽  
Vol 55 (9) ◽  
pp. 1235-1241 ◽  
Author(s):  
Jan Douma
Keyword(s):  
Transverse Isotropy ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Vertical Axis ◽  
Seismic Profile ◽  
Isotropic Layer ◽  
Low Velocity ◽  
Traveltime Inversion ◽  
Axis Of Symmetry ◽  
Transversely Isotropic Layer

Traveltime inversion of multioffset VSP data can be used to determine the depths of the interfaces in layered media. Many inversion schemes, however, assume isotropy and consequently may introduce erroneous structures for anisotropic media. Synthetic traveltime data are computed for layered anisotropic media and inverted assuming isotropic layers. Only the interfaces between these layers are inverted. For a medium consisting of a horizontal isotropic low‐velocity layer on top of a transversely isotropic layer with a horizontal axis of symmetry (e.g., anisotropy due to aligned vertical cracks), 2-D isotropic inversion results in an anticline. For a given axis of symmetry the form of this anticline depends on the azimuth of the source‐borehole direction. The inversion result is a syncline (in 3-D a “bowl” structure), regardless of the azimuth of the source‐borehole direction for a vertical axis of symmetry of the transversely isotropic layer (e.g., anisotropy due to horizontal bedding).

Numerical modeling of EM fields over local anomalies with a vertical axis of symmetry

1990 ◽  
Vol 60 (1-4) ◽  
pp. 53-61 ◽  
Author(s):  
Michael S. Zhdanov ◽  
Vjacheslav V. Spichak ◽  
Leonid Yu. Zaslavsky
Keyword(s):  
Numerical Modeling ◽  
Vertical Axis ◽  
Axis Of Symmetry ◽  
Local Anomalies

Wave forces on vertical axisymmetric bodies

1975 ◽  
Vol 67 (2) ◽  
pp. 369-376 ◽  
Author(s):  
Jared L. Black
Keyword(s):  
Water Surface ◽  
Pressure Field ◽  
Vertical Axis ◽  
Computational Effort ◽  
Circular Cylinders ◽  
Wave Forces ◽  
Local Pressure ◽  
Integral Equation Formulation ◽  
Axis Of Symmetry ◽  
Axisymmetric Bodies

An integral-equation formulation is used to calculate wave forces on bodies having a vertical axis of symmetry. The development enables one to calculate the forces without completely determining the local pressure field, thus offering a considerable reduction of computational effort. Numerical results are presented for a hemisphere at the water surface and vertical circular cylinders.

scholarly journals A numerical study of flow around an impulsively started circular cylinder by a deterministic vortex method

1991 ◽  
Vol 233 ◽  
pp. 243-263 ◽  
Author(s):  
Chien-Cheng Chang ◽  
Ruey-Ling Chern
Keyword(s):  
Finite Difference ◽  
Circular Cylinder ◽  
Numerical Study ◽  
Vortex Method ◽  
Reynolds Numbers ◽  
Numerical Results ◽  
High Reynolds Numbers ◽  
Small Times ◽  
Axis Of Symmetry ◽  
High Reynolds

Impulsively started flow around a circular cylinder at various Reynolds numbers is studied by a deterministic hybrid vortex method. The key feature of the method consists in solving the viscous vorticity equation by interlacing a finite-difference method for diffusion and a vortex-in-cell method for convection. The vorticity is updated along the surface of the cylinder to satisfy the no-slip condition. The present method is basically different from previous applications of vortex methods, which are primarily in the context of random vortex algorithms. The Reynolds numbers of the flows under investigation range from 300 to 106. Numerical results are compared with analytical solutions at small times, and compared with finite-difference solutions and flow visualization results at relatively long times. Satisfactory agreement is found in the evolutions of the separation angles, wake lengths, surface pressure and drag coefficients, streamline patterns, and some velocities on the axis of symmetry behind the circular cylinder. The present hybrid vortex method is highly stable and suffers from little numerical diffusivity, yielding convincing numerical results for unsteady vortical flows at moderately high Reynolds numbers.

Sign in / Sign up

Export Citation Format

Share Document