The precise determination of earthquake location is the fundamental basis in seismological community,and is crucial for analyzing seismic activity and performing seismic tomography.First arrivals are generally used to...The precise determination of earthquake location is the fundamental basis in seismological community,and is crucial for analyzing seismic activity and performing seismic tomography.First arrivals are generally used to practically determine earthquake locations.However,first-arrival traveltimes are not sensitive to focal depths.Moreover,they cannot accurately constrain focal depths.To improve the accuracy,researchers have analyzed the depth phases of earthquake locations.The traveltimes of depth phases are sensitive to focal depths,and the joint inversion of depth phases and direct phases can be implemented to potentially obtain accurate earthquake locations.Generally,researchers can determine earthquake locations in layered models.Because layered models can only represent the first-order feature of subsurface structures,the advantages of joint inversion are not fully explored if layered models are used.To resolve the issue of current joint inversions,we use the traveltimes of three seismic phases to determine earthquake locations in heterogeneous models.The three seismic phases used in this study are the first P-,sPg-and PmP-waves.We calculate the traveltimes of the three seismic phases by solving an eikonal equation with an upwind difference scheme and use the traveltimes to determine earthquake locations.To verify the accuracy of the earthquake location method by the inversion of three seismic phases,we take the 2021 M_(S)6.4 Yangbi,Yunnan earthquake as an example and locate this earthquake using synthetic and real seismic data.Numerical tests demonstrate that the eikonal equation-based earthquake location method,which involves the inversion of multiple phase arrivals,can effectively improve earthquake location accuracy.展开更多
3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic m...3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.展开更多
We present a compact upwind second order scheme for computing the viscosity solution of the Eikonal equation. This new scheme is based on： 1. the numerical observation that classical first order monotone upwind sche...We present a compact upwind second order scheme for computing the viscosity solution of the Eikonal equation. This new scheme is based on： 1. the numerical observation that classical first order monotone upwind schemes for the Eikonal equation yield numerical upwind gradient which is also first order accurate up to singularities; 2. a remark that partial information on the second derivatives of the solution is known and given in the structure of the Eikonal equation and can be used to reduce the size of the stencil. We implement the second order scheme as a correction to the well known sweeping method but it should be applicable to any first order monotone upwind scheme. Care is needed to choose the appropriate stencils to avoid instabilities.展开更多
We develop a second-order continuousfinite element method for solving the static Eikonal equation.It is based on the vanishing viscosity approach with a homotopy method for solving the discretized nonlinear system.Mor...We develop a second-order continuousfinite element method for solving the static Eikonal equation.It is based on the vanishing viscosity approach with a homotopy method for solving the discretized nonlinear system.More specifically,the homotopy method is utilized to decrease the viscosity coefficient gradually,while Newton’s method is applied to compute the solution for each viscosity coefficient.Newton’s method alone converges for just big enough viscosity coefficients on very coarse grids and for simple 1D examples,but the proposed method is much more robust and guarantees the convergence of the nonlinear solver for all viscosity coefficients and for all examples over all grids.Numerical experiments from 1D to 3D are presented to confirm the second-order convergence and the effectiveness of the proposed method on both structured or unstructured meshes.展开更多
In this paper, a reactive dynamic user equilibrium model is extended to simulate two groups of pedestrians traveling on crossing paths in a continuous walking facility. Each group makes path choices to minimize the tr...In this paper, a reactive dynamic user equilibrium model is extended to simulate two groups of pedestrians traveling on crossing paths in a continuous walking facility. Each group makes path choices to minimize the travel cost to its destination in a reactive manner based on instantaneous information. The model consists of a conservation law equation coupled with an Eikonal-type equation for each group. The velocity-density relationship of pedestrian movement is obtained via an experimental method. The model is solved using a finite volume method for the conservation law equation and a fast-marching method for the Eikonal-type equation on unstructured grids. The numerical results verify the rationality of the model and the validity of the numerical method. Based on this continuum model, a number of results, e.g., the formation of strips or moving clusters composed of pedestrians walking to the same destination, are also observed.展开更多
PFMM(perspective fast marching method)是一种有效解决透视投影下从明暗恢复形状(SFS)问题的方法,但是适应条件受限,且对初始数据的精度较为敏感。本文通过对Eikonal方程系数的分析,提出了在透视投影下基于自适应Eikonal方程的PFMM,...PFMM(perspective fast marching method)是一种有效解决透视投影下从明暗恢复形状(SFS)问题的方法,但是适应条件受限,且对初始数据的精度较为敏感。本文通过对Eikonal方程系数的分析,提出了在透视投影下基于自适应Eikonal方程的PFMM,解决了PFMM对初始数据过于依赖的问题,是PFMM的推广。对合成图像的实验表明本文算法比PFMM精度更高,对透视投影下SFS问题可以得到比较好的结果。展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.42174111 and42064004)the Beijing Natural Science Foundation(Grant No.8222033)the Ningxia Science and Technology Leading Talent Training Program(Grant No.2022GKLRLX04)。
文摘The precise determination of earthquake location is the fundamental basis in seismological community,and is crucial for analyzing seismic activity and performing seismic tomography.First arrivals are generally used to practically determine earthquake locations.However,first-arrival traveltimes are not sensitive to focal depths.Moreover,they cannot accurately constrain focal depths.To improve the accuracy,researchers have analyzed the depth phases of earthquake locations.The traveltimes of depth phases are sensitive to focal depths,and the joint inversion of depth phases and direct phases can be implemented to potentially obtain accurate earthquake locations.Generally,researchers can determine earthquake locations in layered models.Because layered models can only represent the first-order feature of subsurface structures,the advantages of joint inversion are not fully explored if layered models are used.To resolve the issue of current joint inversions,we use the traveltimes of three seismic phases to determine earthquake locations in heterogeneous models.The three seismic phases used in this study are the first P-,sPg-and PmP-waves.We calculate the traveltimes of the three seismic phases by solving an eikonal equation with an upwind difference scheme and use the traveltimes to determine earthquake locations.To verify the accuracy of the earthquake location method by the inversion of three seismic phases,we take the 2021 M_(S)6.4 Yangbi,Yunnan earthquake as an example and locate this earthquake using synthetic and real seismic data.Numerical tests demonstrate that the eikonal equation-based earthquake location method,which involves the inversion of multiple phase arrivals,can effectively improve earthquake location accuracy.
基金The authors thank the funds supported by the China National Nuclear Corporation under Grants Nos.WUQNYC2101 and WUHTLM2101-04National Natural Science Foundation of China(42074132,42274154).
文摘3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.
基金partially supported by ONR Grant N00014-02-1-0090ARO MURI Grant W911NF-07-1-0185NSF Grant DMS0811254
文摘We present a compact upwind second order scheme for computing the viscosity solution of the Eikonal equation. This new scheme is based on： 1. the numerical observation that classical first order monotone upwind schemes for the Eikonal equation yield numerical upwind gradient which is also first order accurate up to singularities; 2. a remark that partial information on the second derivatives of the solution is known and given in the structure of the Eikonal equation and can be used to reduce the size of the stencil. We implement the second order scheme as a correction to the well known sweeping method but it should be applicable to any first order monotone upwind scheme. Care is needed to choose the appropriate stencils to avoid instabilities.
基金supported by Natural Science Foundation of Jiangsu Province(Nos.KFR21026,PAF20042)National Natural Science Foundation of China(Nos.GBA20029,GCA20004)+2 种基金Science Challenge Project(No.TZ2018002)National Science and Technology Major Project(No.J2019-II-0007-0027)WH is supported by NSF DMS-1818769.
文摘We develop a second-order continuousfinite element method for solving the static Eikonal equation.It is based on the vanishing viscosity approach with a homotopy method for solving the discretized nonlinear system.More specifically,the homotopy method is utilized to decrease the viscosity coefficient gradually,while Newton’s method is applied to compute the solution for each viscosity coefficient.Newton’s method alone converges for just big enough viscosity coefficients on very coarse grids and for simple 1D examples,but the proposed method is much more robust and guarantees the convergence of the nonlinear solver for all viscosity coefficients and for all examples over all grids.Numerical experiments from 1D to 3D are presented to confirm the second-order convergence and the effectiveness of the proposed method on both structured or unstructured meshes.
基金supported by grants from the Research Grants Council of the Hong Kong Special Administrative Region, China (HKU 7183/06E)the University of Hong Kong (10207394)the National Natural Science Foundation of China (70629001 and 10771134)
文摘In this paper, a reactive dynamic user equilibrium model is extended to simulate two groups of pedestrians traveling on crossing paths in a continuous walking facility. Each group makes path choices to minimize the travel cost to its destination in a reactive manner based on instantaneous information. The model consists of a conservation law equation coupled with an Eikonal-type equation for each group. The velocity-density relationship of pedestrian movement is obtained via an experimental method. The model is solved using a finite volume method for the conservation law equation and a fast-marching method for the Eikonal-type equation on unstructured grids. The numerical results verify the rationality of the model and the validity of the numerical method. Based on this continuum model, a number of results, e.g., the formation of strips or moving clusters composed of pedestrians walking to the same destination, are also observed.
文摘PFMM(perspective fast marching method)是一种有效解决透视投影下从明暗恢复形状(SFS)问题的方法,但是适应条件受限,且对初始数据的精度较为敏感。本文通过对Eikonal方程系数的分析,提出了在透视投影下基于自适应Eikonal方程的PFMM,解决了PFMM对初始数据过于依赖的问题,是PFMM的推广。对合成图像的实验表明本文算法比PFMM精度更高,对透视投影下SFS问题可以得到比较好的结果。
