The ratio between vertical and radial amplitudes of Rayleigh waves(hereafter,the Rayleigh wave ZH ratio)is an important parameter used to constrain structures beneath seismic stations.Some previous studies have explor...The ratio between vertical and radial amplitudes of Rayleigh waves(hereafter,the Rayleigh wave ZH ratio)is an important parameter used to constrain structures beneath seismic stations.Some previous studies have explored crust and upper mantle structures by joint inversion of the Rayleigh wave ZH ratio and surface wave dispersion.However,all these studies have used a 1-D depth sensitivity kernel,and this kernel may lack precision when the structure varies a great deal laterally.Here,we present a systematic investigation of the two-dimensional(2-D)Rayleigh wave ZH ratio kernel based on the adjoint-wavefield method and perform two synthetic tests using the new kernel.The 2-D ZH ratio kernel is consistent with the traditional 1-D sensitivity kernel but has an asymmetric pattern with a preferred orientation toward the source.The predominant effect caused by heterogeneity can clearly be seen from kernels calculated from models with 2-D heterogeneities,which confirms the necessity of using the new 2-D kernel in some complex regions.Inversion tests using synthetic data show that the 2-D ZH ratio kernel has the potential to resolve small anomalies as well as complex lateral structures.展开更多
In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation ...In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation in inversion. A feasible way to avoid the excessive storage demand is to reconstruct the source wavefield backward in time by storing the entire history of the wavefield in perfectly matched layers. In this paper, we make full use of the elementwise global property of the Laplace operator of the spectral element method (SEM) and propose an efficient source wavefield reconstruction method at the cost of storing the wavefield history only at single boundary layer nodes. Numerical experiments indicate that the accuracy of the proposed method is identical to that of the conventional method and is independent of the order of the Lagrange polynomials, the element type, and the temporal discretization method. In contrast, the memory-saving ratios of the conventional method versus our method is at least N when using either quadrilateral or hexahedron elements, respectively, where N is the order of the Lagrange polynomials used in the SEM. A higher memorysaving ratio is achieved with triangular elements versus quadrilaterals. The new method is applied to reverse time migration by considering the Marmousi model as a benchmark. Numerical results demonstrate that the method is able to provide the same result as the conventional method but with about 1/25 times lower storage demand. With the proposed wavefield reconstruction method, the storage demand is dramatically reduced;therefore, in-core memory storage is feasible even for large-scale three-dimensional adjoint inversion problems.展开更多
基金This study was funded by the National Key R&D Program of China(2016YFC0600301,2018YFC1503400)the National Natural Science Foundation of China(41790464)+1 种基金Natural Science Foundation of Jiangsu Province of China(BK20190499)the Fundamental Research Funds for the Central Universities(2019B0071428).
文摘The ratio between vertical and radial amplitudes of Rayleigh waves(hereafter,the Rayleigh wave ZH ratio)is an important parameter used to constrain structures beneath seismic stations.Some previous studies have explored crust and upper mantle structures by joint inversion of the Rayleigh wave ZH ratio and surface wave dispersion.However,all these studies have used a 1-D depth sensitivity kernel,and this kernel may lack precision when the structure varies a great deal laterally.Here,we present a systematic investigation of the two-dimensional(2-D)Rayleigh wave ZH ratio kernel based on the adjoint-wavefield method and perform two synthetic tests using the new kernel.The 2-D ZH ratio kernel is consistent with the traditional 1-D sensitivity kernel but has an asymmetric pattern with a preferred orientation toward the source.The predominant effect caused by heterogeneity can clearly be seen from kernels calculated from models with 2-D heterogeneities,which confirms the necessity of using the new 2-D kernel in some complex regions.Inversion tests using synthetic data show that the 2-D ZH ratio kernel has the potential to resolve small anomalies as well as complex lateral structures.
基金financial support for this work contributed by the National Key Research and Development Program of China (grant numbers 2016YFC0600101 and 2016YFC 0600201)the National Natural Science Foundation of China (grant numbers 41874065, 41604076, 41674102, 41674095, 41522401, 41574082, and 41774097)
文摘In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation in inversion. A feasible way to avoid the excessive storage demand is to reconstruct the source wavefield backward in time by storing the entire history of the wavefield in perfectly matched layers. In this paper, we make full use of the elementwise global property of the Laplace operator of the spectral element method (SEM) and propose an efficient source wavefield reconstruction method at the cost of storing the wavefield history only at single boundary layer nodes. Numerical experiments indicate that the accuracy of the proposed method is identical to that of the conventional method and is independent of the order of the Lagrange polynomials, the element type, and the temporal discretization method. In contrast, the memory-saving ratios of the conventional method versus our method is at least N when using either quadrilateral or hexahedron elements, respectively, where N is the order of the Lagrange polynomials used in the SEM. A higher memorysaving ratio is achieved with triangular elements versus quadrilaterals. The new method is applied to reverse time migration by considering the Marmousi model as a benchmark. Numerical results demonstrate that the method is able to provide the same result as the conventional method but with about 1/25 times lower storage demand. With the proposed wavefield reconstruction method, the storage demand is dramatically reduced;therefore, in-core memory storage is feasible even for large-scale three-dimensional adjoint inversion problems.
