Fast Inversion of the Earthquake Rupture Processes with Complicated Velocity Structure: An Application to the Earthquake of 2017 Mw 6.5 Jiuzhaigou, China

2021 ◽  
pp. 2150030
Author(s):  
Jiemin Wang ◽  
Haitao Yin ◽  
Zhijun Feng ◽  
Pifeng Ma ◽  
Liang Wang

Due to the limitation of seismic station coverage or the network transport interrupted when the earthquake occurred, an accurate seismic shakemap may not be released to the public quickly. When the near-source observed waveforms for the intensity prediction technology used are incomplete, we synthesize the seismic waveform into observation waveforms. An accurate seismic rupture process is necessary to synthesize virtual station observations. So, we should release the rupture process as soon as possible after a large earthquake. Most large earthquakes occur at the junction of two or three tectonic terranes. With violent tectonic movements, fault basins and uplift zones are distributed on the edge of the plateau. With complex structural conditions, the 1D layered half-space velocity structure model could not meet the requirement of earthquake rupture process inversion. It takes much time to calculate 3D Green’s function with a 3D velocity model for the complete waveform inversion of the earthquake rupture process. To rapidly invert the rupture process as accurately as possible, according to the geological conditions of the station, we calculated several Green’s function libraries in advance. We extracted Green’s functions from these libraries for each site based on the sites’ coordinates once an earthquake occurs. The time we spend in extracting Green’s functions from several Green libraries equals that we spend in extracting Green’s functions from one single library. The applicability of this method was tested in the 2017 Jiuzhaigou M6.5 earthquake with complex structural conditions in the mountain uplift zone. With our model, the time we spent in calculating the rupture process was almost the same as that we spent with the 1D velocity structure model, which was far less than that we could have spent in calculating 3D Green’s function. The degree of fitting between the synthetic data and the observation data of our model was much higher than the fitting of the 1D velocity model, which means that the earthquake rupture process we determined was more reliable.

Geophysics ◽  
2019 ◽  
Vol 84 (6) ◽  
pp. Q57-Q72
Author(s):  
Yanadet Sripanich ◽  
Ivan Vasconcelos ◽  
Kees Wapenaar

The Marchenko method retrieves Green’s functions between the acquisition surface and any arbitrary point in the medium. The process generally involves solving an inversion starting with an initial focusing function, e.g., a direct-wave Green’s function from the desired subsurface position, typically obtained using an approximate velocity model. We have formulated the Marchenko method in the time-imaging domain. In that domain, we recognize that the traveltime of the direct-wave Green’s function is related to the Cheop’s traveltime pyramid commonly used in time-domain processing, which in turn can be readily obtained from the local slopes of the common-midpoint gathers. This observation allows us to substitute the velocity-model-based initial focusing operator with that from a data-driven slope estimation process. Moreover, we found that working in the time-imaging domain allows for the specification of the desired subsurface position in terms of vertical time, which is connected to the Cartesian depth position via the time-to-depth conversion. Our results suggest that the prior velocity model is only required when specifying the position in depth, but this requirement can be circumvented by making use of the time-imaging domain within its usual assumptions (e.g., mild lateral heterogeneity). Provided that those assumptions are satisfied, the estimated Green’s functions from the proposed method have comparable quality to those obtained with the knowledge of a prior velocity model.


Author(s):  
Norman J. Morgenstern Horing

Chapter 09 Nonequilibrium Green’s functions (NEGF), including coupled-correlated (C) single- and multi-particle Green’s functions, are defined as averages weighted with the time-development operator U(t0+τ,t0). Linear conductivity is exhibited as a two-particle equilibrium Green’s function (Kubo-type formulation). Admitting particle sources (S:η,η+) and non-conservation of number, the non-equilibrium multi-particle Green’s functions are constructed with numbers of creation and annihilation operators that may differ, and they may be derived as variational derivatives with respect to sources η,η+ of a generating functional eW=TrU(t0+τ,t0)CS/TrU(t0+τ,t0)C. (In the non-interacting case this yields the n-particle Green’s function as a permanent/determinant of single-particle Green’s functions.) These variational relations yield a symmetric set of multi-particle Green’s function equations. Cumulants and the Linked Cluster Theorem are discussed and the Random Phase Approximation (RPA) is derived variationally. Schwinger’s variational differential formulation of perturbation theories for the Green’s function, self-energy, vertex operator, and also shielded potential perturbation theory, are reviewed. The Langreth Algebra arises from analytic continuation of integration of products of Green’s functions in imaginary time to the real-time axis with time-ordering along the integration contour in the complex time plane. An account of the Generalized Kadanoff-Baym Ansatz is presented.


Author(s):  
Norman J. Morgenstern Horing

Multiparticle thermodynamic Green’s functions, defined in terms of grand canonical ensemble averages of time-ordered products of creation and annihilation operators, are interpreted as tracing the amplitude for time-developing correlated interacting particle motions taking place in the background of a thermal ensemble. Under equilibrium conditions, time-translational invariance permits the one-particle thermal Green’s function to be represented in terms of a single frequency, leading to a Lehmann spectral representation whose frequency poles describe the energy spectrum. This Green’s function has finite values for both t>t′ and t<t′ (unlike retarded Green’s functions), and the two parts G1> and G1< (respectively) obey a simple proportionality relation that facilitates the introduction of a spectral weight function: It is also interpreted in terms of a periodicity/antiperiodicity property of a modified Green’s function in imaginary time capable of a Fourier series representation with imaginary (Matsubara) frequencies. The analytic continuation from imaginary time to real time is discussed, as are related commutator/anticommutator functions, also retarded/advanced Green’s functions, and the spectral weight sum rule is derived. Statistical thermodynamic information is shown to be embedded in physical features of the one- and two-particle thermodynamic Green’s functions.


Author(s):  
Klaus Morawetz

The method of the equation of motion is used to derive the Martin–Schwinger hierarchy for the nonequilibrium Green’s functions. The formal closure of the hierarchy is reached by using the selfenergy which provides a recipe for how to construct selfenergies from approximations of the two-particle Green’s function. The Langreth–Wilkins rules for a diagrammatic technique are shown to be equivalent to the weakening of initial correlations. The quantum transport equations are derived in the general form of Kadanoff and Baym equations. The information contained in the Green’s function is discussed. In equilibrium this leads to the Matsubara diagrammatic technique.


1994 ◽  
Vol 03 (02) ◽  
pp. 523-589 ◽  
Author(s):  
T.T.S. KUO ◽  
YIHARN TZENG

We present an elementary and fairly detailed review of several Green’s function methods for treating nuclear and other many-body systems. We first treat the single-particle Green’s function, by way of which some details concerning linked diagram expansion, rules for evaluating Green’s function diagrams and solution of the Dyson’s integral equation for Green’s function are exhibited. The particle-particle hole-hole (pphh) Green’s function is then considered, and a specific time-blocking technique is discussed. This technique enables us to have a one-frequency Dyson’s equation for the pphh and similarly for other Green’s functions, thus considerably facilitating their calculation. A third type of Green’s function considered is the particle-hole Green’s function. RPA and high order RPA are treated, along with examples for setting up particle-hole RPA equations. A general method for deriving a model-space Dyson’s equation for Green’s functions is discussed. We also discuss a method for determining the normalization of Green’s function transition amplitudes based on its vertex function. Some applications of Green’s function methods to nuclear structure and recent deep inelastic lepton-nucleus scattering are addressed.


1999 ◽  
Author(s):  
Paul E. Barbone

Abstract We derive a one-way wave equation representation of the “free space” Green’s function for an inhomogeneous medium. Our representation results from an asymptotic expansion in inverse powers of the wavenumber. Our representation takes account of losses due to scattering in all directions, even though only one-way operators are used.


1987 ◽  
Vol 30 (1) ◽  
pp. 28-35 ◽  
Author(s):  
P. W. Eloe

AbstractLet G(x,s) be the Green's function for the boundary value problem y(n) = 0, Ty = 0, where Ty = 0 represents boundary conditions at two points. The signs of G(x,s) and certain of its partial derivatives with respect to x are determined for two classes of boundary value problems. The results are also carried over to analogous classes of boundary value problems for difference equations.


2019 ◽  
Vol 3 (2) ◽  
pp. 36
Author(s):  
Ifan Johnston ◽  
Vassili Kolokoltsov

We look at estimates for the Green’s function of time-fractional evolution equations of the form D 0 + * ν u = L u , where D 0 + * ν is a Caputo-type time-fractional derivative, depending on a Lévy kernel ν with variable coefficients, which is comparable to y - 1 - β for β ∈ ( 0 , 1 ) , and L is an operator acting on the spatial variable. First, we obtain global two-sided estimates for the Green’s function of D 0 β u = L u in the case that L is a second order elliptic operator in divergence form. Secondly, we obtain global upper bounds for the Green’s function of D 0 β u = Ψ ( - i ∇ ) u where Ψ is a pseudo-differential operator with constant coefficients that is homogeneous of order α . Thirdly, we obtain local two-sided estimates for the Green’s function of D 0 β u = L u where L is a more general non-degenerate second order elliptic operator. Finally we look at the case of stable-like operator, extending the second result from a constant coefficient to variable coefficients. In each case, we also estimate the spatial derivatives of the Green’s functions. To obtain these bounds we use a particular form of the Mittag-Leffler functions, which allow us to use directly known estimates for the Green’s functions associated with L and Ψ , as well as estimates for stable densities. These estimates then allow us to estimate the solutions to a wide class of problems of the form D 0 ( ν , t ) u = L u , where D ( ν , t ) is a Caputo-type operator with variable coefficients.


Sign in / Sign up

Export Citation Format

Share Document