We have developed a Liouville partial-differential-equation (PDE)-based method for computing complex-valued eikonals in real phase space in the multivalued sense in attenuating media with frequency-independent qualify factors, where the new method computes the real and imaginary parts of the complex-valued eikonal in two steps by solving Liouville equations in real phase space. Because the earth is composed of attenuating materials, seismic waves usually attenuate so that seismic data processing calls for properly treating the resulting energy losses and phase distortions of wave propagation. In the regime of high-frequency asymptotics, the complex-valued eikonal is one essential ingredient for describing wave propagation in attenuating media because this unique quantity summarizes two wave properties into one function: Its real part describes the wave kinematics and its imaginary part captures the effects of phase dispersion and amplitude attenuation. Because some popular ordinary-differential-equation (ODE)-based ray-tracing methods for computing complex-valued eikonals in real space distribute the eikonal function irregularly in real space, we are motivated to develop PDE-based Eulerian methods for computing such complex-valued eikonals in real space on regular meshes. Therefore, we solved novel paraxial Liouville PDEs in real phase space so that we can compute the real and imaginary parts of the complex-valued eikonal in the multivalued sense on regular meshes. We call the resulting method the Liouville PDE method for complex-valued multivalued eikonals in attenuating media; moreover, this new method provides a unified framework for Eulerianizing several popular approximate real-space ray-tracing methods for complex-valued eikonals, such as viscoacoustic ray tracing, real viscoelastic ray tracing, and real elastic ray tracing. In addition, we also provide Liouville PDE formulations for computing multivalued ray amplitudes in a weakly viscoacoustic medium. Numerical examples, including a synthetic gas-cloud model, demonstrate that our methods yield highly accurate complex-valued eikonals in the multivalued sense.
more »
« less
Eulerian partial-differential-equation methods for complex-valued eikonals in attenuating media
Seismic waves in earth media usually undergo attenuation, causing energy losses and phase distortions. In the regime of high-frequency asymptotics, a complex-valued eikonal is an essential ingredient for describing wave propagation in attenuating media, where the real and imaginary parts of the eikonal function capture dispersion effects and amplitude attenuation of seismic waves, respectively. Conventionally, such a complex-valued eikonal is mainly computed either by tracing rays exactly in complex space or by tracing rays approximately in real space so that the resulting eikonal is distributed irregularly in real space. However, seismic data processing methods, such as prestack depth migration and tomography, usually require uniformly distributed complex-valued eikonals. Therefore, we have developed a unified framework to Eulerianize several popular approximate real-space ray-tracing methods for complex-valued eikonals so that the real and imaginary parts of the eikonal function satisfy the classic real-space eikonal equation and a novel real-space advection equation, respectively, and we dub the resulting method the Eulerian partial-differential-equation method. We further develop highly efficient high-order methods to solve these two equations by using the factorization idea and the Lax-Friedrichs weighted essentially nonoscillatory schemes. Numerical examples demonstrate that our method yields highly accurate complex-valued eikonals, analogous to those from ray-tracing methods. Our methods can be useful for migration and tomography in attenuating media.
more »
« less
- Award ID(s):
- 2012046
- PAR ID:
- 10278867
- Date Published:
- Journal Name:
- GEOPHYSICS
- Volume:
- 86
- Issue:
- 4
- ISSN:
- 0016-8033
- Page Range / eLocation ID:
- T179 to T192
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
First-arrival traveltime tomography is an essential method for obtaining near-surface velocity models. The adjoint-state first-arrival traveltime tomography is appealing due to its straightforward implementation, low computational cost, and low memory consumption. Because solving the point-source isotropic eikonal equation by either ray tracers or eikonal solvers intrinsically corresponds to emanating discrete rays from the source point, the resulting traveltime gradient is singular at the source point, and we denote such a singular pattern the imprint of ray-illumination. Because the adjoint-state equation propagates traveltime residuals back to the source point according to the negative traveltime gradient, the resulting adjoint state will inherit such an imprint of ray-illumination, leading to singular gradient-descent directions when updating the velocity model in the adjoint-state traveltime tomography. To mitigate this imprint, we solve the adjoint-state equation twice but with different boundary conditions: one being taken to be regular data residuals and the other taken to be ones uniformly, so that we are able to use the latter adjoint state to normalize the regular adjoint state and we further use the normalized quantity to serve as the gradient direction to update the velocity model; we call this process ray-illumination compensation. To overcome the issue of limited aperture, we have developed a spatially varying regularization method to stabilize the new gradient direction. A synthetic example demonstrates that our method is able to mitigate the imprint of ray-illumination, remove the footprint effect near source points, and provide uniform velocity updates along raypaths. A complex example extracted from the Marmousi2 model and a migration example illustrate that the new method accurately recovers the velocity model and that an offset-dependent inversion strategy can further improve the quality of recovered velocity models.more » « less
-
more » « less
SUMMARY We introduce a scheme for probabilistic hypocentre inversion with Stein variational inference. Our approach uses a differentiable forward model in the form of a physics informed neural network, which we train to solve the Eikonal equation. This allows for rapid approximation of the posterior by iteratively optimizing a collection of particles against a kernelized Stein discrepancy. We show that the method is well-equipped to handle highly multimodal posterior distributions, which are common in hypocentral inverse problems. A suite of experiments is performed to examine the influence of the various hyperparameters. Once trained, the method is valid for any seismic network geometry within the study area without the need to build traveltime tables. We show that the computational demands scale efficiently with the number of differential times, making it ideal for large-N sensing technologies like Distributed Acoustic Sensing. The techniques outlined in this manuscript have considerable implications beyond just ray tracing procedures, with the work flow applicable to other fields with computationally expensive inversion procedures such as full waveform inversion.
-
more » « less
Abstract To complement velocity distributions, seismic attenuation provides additional important information on fluid properties of hydrocarbon reservoirs in exploration seismology, as well as temperature distributions, partial melting, and water content within the crust and mantle in earthquake seismology. Full waveform inversion (FWI), as one of the state‐of‐the‐art seismic imaging techniques, can produce high‐resolution constraints for subsurface (an)elastic parameters by minimizing the difference between observed and predicted seismograms. Traditional waveform inversion for attenuation is commonly based on the standard‐linear‐solid (SLS) wave equation, in which case the quality factor (
Q ) has to be converted to stress and strain relaxation times. When using multiple attenuation mechanisms in the SLS method, it is difficult to directly estimate these relaxation time parameters. Based on a time domain complex‐valued viscoacoustic wave equation, we present an FWI framework for simultaneously estimating subsurfaceP wave velocity and attenuation distributions. BecauseQ is explicitly incorporated into the viscoacoustic wave equation, we directly deriveP wave velocity andQ sensitivity kernels using the adjoint‐state method and simultaneously estimate their subsurface distributions. By analyzing the Gauss‐Newton Hessian, we observe strong interparameter crosstalk, especially the leakage from velocity toQ . We approximate the Hessian inverse using a preconditioned L‐BFGS method in viscoacoustic FWI, which enables us to successfully reduce interparameter crosstalk and produce accurate velocity and attenuation models. Numerical examples demonstrate the feasibility and robustness of the proposed method for simultaneously mapping complex velocity andQ distributions in the subsurface. -
more » « less
Abstract We revisit the inverse source problem in a two dimensional absorbing and scattering medium and present a direct reconstruction method, which does not require iterative solvability of the forward problem, using measurements of the radiating flux at the boundary. The attenuation and scattering coefficients are known and the unknown source is isotropic. The approach is based on the Cauchy problem for a Beltrami-like equation for the sequence valued maps, and extends the original ideas of Bukhgeim from the non-scattering to the scattering media. We demonstrate the feasibility of the method in a numerical experiment in which the scattering is modeled by the two dimensional Henyey–Greenstein kernel with parameters meaningful in optical tomography.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1190/GEO2020-0659.1
