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.展开更多
Seismic traveltime tomographic inversion has played an important role in detecting the internal structure of the solid earth. We use a set of blocks to approximate geologically complex media that cannot be well descri...Seismic traveltime tomographic inversion has played an important role in detecting the internal structure of the solid earth. We use a set of blocks to approximate geologically complex media that cannot be well described by layered models or cells. The geological body is described as an aggregate of arbitrarily shaped blocks,which are separated by triangulated interfaces. We can describe the media as homogenous or heterogeneous in each block. We define the velocities at the given rectangle grid points for each block,and the heterogeneous velocities in each block can be calculated by a linear interpolation algorithm. The parameters of the velocity grid positions are independent of the model parameterization,which is advantageous in the joint inversion of the velocities and the node depths of an interface. We implement a segmentally iterative ray tracer to calculate traveltimes in the 3D heterogeneous block models.The damped least squares method is employed in seismic traveltime inversion,which includes the partial derivatives of traveltime with respect to the depths of nodes in the triangulated interfaces and velocities defined in rectangular grids. The numerical tests indicate that the node depths of a triangulated interface and homogeneous velocity distributions can be well inverted in a stratified model.展开更多
The Emeishan large igneous province(hereafter named by its acronym ELIP) is the first accepted large igneous region in China.The current study tries to reconstruct the density structure of the crust in this region. Fo...The Emeishan large igneous province(hereafter named by its acronym ELIP) is the first accepted large igneous region in China.The current study tries to reconstruct the density structure of the crust in this region. For this purpose, we conducted the gravity survey along an 800-km-long profile, which stretched laterally along the latitude 27°N from Lijiang(Yunnan province) to Guiyang(Guizhou province). The fieldwork included 338 gravity measurements distributed from the inner zone to the outer zone of the mantle plume head.After a series of gravity reductions, we calculated the Bouguer gravity anomaly and then constructed the density model for ELIP by iterative forward modeling from an initial density model tightly constrained by wide-angle seismic reflection data. The topography of the Moho, here physically interpreted as a density discontinuity of ~0.4 g·cm^(–3), gradually rises from the inner zone(~50 km deep) to the outer zone(~40 km), describes a thicker crust in the inner zone than in any other segment of the profile and largely reproduces the shape of the Bouguer gravity anomaly curve. Both the Bouguer gravity and the density structure show significant differences with respect to the inner zone and the other two zones of ELIP according to the commonly accepted partition of the Emeishan area. A thicker and denser middle-lower crust seems to be the main feature of the western section of the profile, which is likely related to its mafic magmatic composition due to magmatic underplating of the Permian mantle plume.展开更多
Compared with the surface,the deep environment has the advantages of allowing“super-quiet and ultra-clean”-geophysical field observation with low vibration noise and little electromagnetic interference,which are con...Compared with the surface,the deep environment has the advantages of allowing“super-quiet and ultra-clean”-geophysical field observation with low vibration noise and little electromagnetic interference,which are conducive to the realization of long-term and high-precision observation of multi-physical fields,thus enabling the solution of a series of geoscience problems.In the Panyidong Coal Mine,where there are extensive underground tunnels at the depth of 848 m below sea level,we carried out the first deep-underground geophysical observations,including radioactivity,gravity,magnetic,magnetotelluric,background vibration and six-component seismic observations.We concluded from these measurements that(1)the background of deep subsurface gravity noise in the long-period frequency band less than 2 Hz is nearly two orders of magnitude weaker than that in the surface observation environment;(2)the underground electric field is obviously weaker than the surface electric field,and the relatively high frequency of the underground field,greater than 1 Hz,is more than two orders of magnitude weaker than that of the surface electric field;the east-west magnetic field underground is approximately the same as that at the surface;the relatively high-frequency north-south magnetic field underground,below 10 Hz,is at least one order of magnitude lower than that at the surface,showing that the underground has a clean electromagnetic environment;(3)in addition to the highfrequency and single-frequency noises introduced by underground human activities,the deep underground space has a significantly lower background vibration noise than the surface,which is very beneficial to the detection of weak earthquake and gravity signals;and(4)the underground roadway support system built with ferromagnetic material interferes the geomagnetic field.We also found that for deep observation in the“ultra-quiet and ultra-clean”environment,the existing geophysical equipment and observation technology have problems of poor adaptability and insufficient precision as well as data cleaning problems,such as the effective separation of the signal and noise of deep observation data.It is also urgent to interpret and comprehensively utilize these high-precision multi-physics observation data.展开更多
基金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.
基金supported financially by the Ministry of Science and Technology of China(2011CB808904)the National Natural Science Foundation of China(Nos.41021063,41174075,41004034,41174043,and 41274090)
文摘Seismic traveltime tomographic inversion has played an important role in detecting the internal structure of the solid earth. We use a set of blocks to approximate geologically complex media that cannot be well described by layered models or cells. The geological body is described as an aggregate of arbitrarily shaped blocks,which are separated by triangulated interfaces. We can describe the media as homogenous or heterogeneous in each block. We define the velocities at the given rectangle grid points for each block,and the heterogeneous velocities in each block can be calculated by a linear interpolation algorithm. The parameters of the velocity grid positions are independent of the model parameterization,which is advantageous in the joint inversion of the velocities and the node depths of an interface. We implement a segmentally iterative ray tracer to calculate traveltimes in the 3D heterogeneous block models.The damped least squares method is employed in seismic traveltime inversion,which includes the partial derivatives of traveltime with respect to the depths of nodes in the triangulated interfaces and velocities defined in rectangular grids. The numerical tests indicate that the node depths of a triangulated interface and homogeneous velocity distributions can be well inverted in a stratified model.
基金supported by the National Key Research and Development Project of China (grant 2016YFC0600302)the National Natural Science Foundation of China (grant 41774100)National Basic Research Program of China (973 Program, grant 2011CB808904)
文摘The Emeishan large igneous province(hereafter named by its acronym ELIP) is the first accepted large igneous region in China.The current study tries to reconstruct the density structure of the crust in this region. For this purpose, we conducted the gravity survey along an 800-km-long profile, which stretched laterally along the latitude 27°N from Lijiang(Yunnan province) to Guiyang(Guizhou province). The fieldwork included 338 gravity measurements distributed from the inner zone to the outer zone of the mantle plume head.After a series of gravity reductions, we calculated the Bouguer gravity anomaly and then constructed the density model for ELIP by iterative forward modeling from an initial density model tightly constrained by wide-angle seismic reflection data. The topography of the Moho, here physically interpreted as a density discontinuity of ~0.4 g·cm^(–3), gradually rises from the inner zone(~50 km deep) to the outer zone(~40 km), describes a thicker crust in the inner zone than in any other segment of the profile and largely reproduces the shape of the Bouguer gravity anomaly curve. Both the Bouguer gravity and the density structure show significant differences with respect to the inner zone and the other two zones of ELIP according to the commonly accepted partition of the Emeishan area. A thicker and denser middle-lower crust seems to be the main feature of the western section of the profile, which is likely related to its mafic magmatic composition due to magmatic underplating of the Permian mantle plume.
基金The research was supported by the National Natural Science Foundation of China （Grant No. 41130419, 41240027, 41374061, 41204048） and China Postdoctoral Science Foundation （Grant No. 2012M510533）.
基金supported by the National Natural Science Foundation of China (Grant Nos.62127815,42150201,U1839208)
文摘Compared with the surface,the deep environment has the advantages of allowing“super-quiet and ultra-clean”-geophysical field observation with low vibration noise and little electromagnetic interference,which are conducive to the realization of long-term and high-precision observation of multi-physical fields,thus enabling the solution of a series of geoscience problems.In the Panyidong Coal Mine,where there are extensive underground tunnels at the depth of 848 m below sea level,we carried out the first deep-underground geophysical observations,including radioactivity,gravity,magnetic,magnetotelluric,background vibration and six-component seismic observations.We concluded from these measurements that(1)the background of deep subsurface gravity noise in the long-period frequency band less than 2 Hz is nearly two orders of magnitude weaker than that in the surface observation environment;(2)the underground electric field is obviously weaker than the surface electric field,and the relatively high frequency of the underground field,greater than 1 Hz,is more than two orders of magnitude weaker than that of the surface electric field;the east-west magnetic field underground is approximately the same as that at the surface;the relatively high-frequency north-south magnetic field underground,below 10 Hz,is at least one order of magnitude lower than that at the surface,showing that the underground has a clean electromagnetic environment;(3)in addition to the highfrequency and single-frequency noises introduced by underground human activities,the deep underground space has a significantly lower background vibration noise than the surface,which is very beneficial to the detection of weak earthquake and gravity signals;and(4)the underground roadway support system built with ferromagnetic material interferes the geomagnetic field.We also found that for deep observation in the“ultra-quiet and ultra-clean”environment,the existing geophysical equipment and observation technology have problems of poor adaptability and insufficient precision as well as data cleaning problems,such as the effective separation of the signal and noise of deep observation data.It is also urgent to interpret and comprehensively utilize these high-precision multi-physics observation data.
