The airborne electromagnetic (AEM) method has a high sampling rate and survey flexibility. However, traditional numerical modeling approaches must use high-resolution physical grids to guarantee modeling accuracy, e...The airborne electromagnetic (AEM) method has a high sampling rate and survey flexibility. However, traditional numerical modeling approaches must use high-resolution physical grids to guarantee modeling accuracy, especially for complex geological structures such as anisotropic earth. This can lead to huge computational costs. To solve this problem, we propose a spectral-element (SE) method for 3D AEM anisotropic modeling, which combines the advantages of spectral and finite-element methods. Thus, the SE method has accuracy as high as that of the spectral method and the ability to model complex geology inherited from the finite-element method. The SE method can improve the modeling accuracy within discrete grids and reduce the dependence of modeling results on the grids. This helps achieve high-accuracy anisotropic AEM modeling. We first introduced a rotating tensor of anisotropic conductivity to Maxwell's equations and described the electrical field via SE basis functions based on GLL interpolation polynomials. We used the Galerkin weighted residual method to establish the linear equation system for the SE method, and we took a vertical magnetic dipole as the transmission source for our AEM modeling. We then applied fourth-order SE calculations with coarse physical grids to check the accuracy of our modeling results against a 1D semi-analytical solution for an anisotropic half-space model and verified the high accuracy of the SE. Moreover, we conducted AEM modeling for different anisotropic 3D abnormal bodies using two physical grid scales and three orders of SE to obtain the convergence conditions for different anisotropic abnormal bodies. Finally, we studied the identification of anisotropy for single anisotropic abnormal bodies, anisotropic surrounding rock, and single anisotropic abnormal body embedded in an anisotropic surrounding rock. This approach will play a key role in the inversion and interpretation of AEM data collected in regions with anisotropic geology.展开更多
Frequency-domain airborne electromagnetics is a proven geophysical exploration method.Presently,the interpretation is mainly based on resistivity-depth imaging and onedimensional layered inversion;nevertheless,it is d...Frequency-domain airborne electromagnetics is a proven geophysical exploration method.Presently,the interpretation is mainly based on resistivity-depth imaging and onedimensional layered inversion;nevertheless,it is difficult to obtain satisfactory results for two- or three-dimensional complex earth structures using 1D methods.3D forward modeling and inversion can be used but are hampered by computational limitations because of the large number of data.Thus,we developed a 2.5D frequency-domain airborne electromagnetic forward modeling and inversion algorithm.To eliminate the source singularities in the numerical simulations,we split the fields into primary and secondary fields.The primary fields are calculated using homogeneous or layered models with analytical solutions,and the secondary(scattered) fields are solved by the finite-element method.The linear system of equations is solved by using the large-scale sparse matrix parallel direct solver,which greatly improves the computational efficiency.The inversion algorithm was based on damping leastsquares and singular value decomposition and combined the pseudo forward modeling and reciprocity principle to compute the Jacobian matrix.Synthetic and field data were used to test the effectiveness of the proposed method.展开更多
Numerical simulation in the frequency-space domain has inherent advantages, such as: it is possible to simulate wave propagation from multiple sources simultaneously; there are no cumulative errors; only the interest...Numerical simulation in the frequency-space domain has inherent advantages, such as: it is possible to simulate wave propagation from multiple sources simultaneously; there are no cumulative errors; only the interesting frequencies can be selected; and it is more suitable for wave propagation in viscoelastic media. The only obstacle to using the method is the requirement of huge computer storage. We extend the compressed format for storing the coefficient matrix. It can reduce the required computer storage dramatically. We get the optimal coefficients by least-squares method to suppress the numerical dispersion and adopt the perfectly matched layer (PML) boundary conditions to eliminate the artificial boundary reflections. Using larger grid intervals decreases computer storage requirements and provides high computational efficiency. Numerical experiments demonstrate that these means are economic and effective, providing a good basis for elastic wave imaging and inversion.展开更多
Non-spreading nature of Bessel spatiotemporal wavepackets is theoretically and experimentally investigated and orders of magnitude improvement in the spatiotemporal spreading has been demonstrated.The spatiotemporal c...Non-spreading nature of Bessel spatiotemporal wavepackets is theoretically and experimentally investigated and orders of magnitude improvement in the spatiotemporal spreading has been demonstrated.The spatiotemporal confinement provided by the Bessel spatiotemporal wavepacket is further exploited to transport transverse orbital angular momentum through embedding spatiotemporal optical vortex into the Bessel spatiotemporal wavepacket, constructing a new type of wavepacket: Bessel spatiotemporal optical vortex. Both numerical and experimental results demonstrate that spatiotemporal vortex structure can be well maintained and confined through much longer propagation. High order spatiotemporal optical vortices can also be better confined in the spatiotemporal domain and prevented from further breaking up, overcoming a potential major obstacle for future applications of spatiotemporal vortex.展开更多
Motivated by the recent work that the periodicity of a black hole is responsible for the area spectrum,we exclusively utilize the period of motion of an outgoing wave,which is shown to be related to the vibrational fr...Motivated by the recent work that the periodicity of a black hole is responsible for the area spectrum,we exclusively utilize the period of motion of an outgoing wave,which is shown to be related to the vibrational frequency of the perturbed black hole,to study area spectra of a non-rotating BTZ black hole and a rotating BTZ black hole.It is found that the area spectra and entropy spectra for both space times are equally spaced.In addition,we find that though the entropy spectra of the 3-dimensional BTZ black holes take the same form as those of the 4-dimensional black holes,the area spectra depend on the dimension of space times.Our result confirms that the entropy spectrum of a black hole is more fundamental than the area spectrum.展开更多
基金financially supported by the Key Program of National Natural Science Foundation of China(No.41530320)China Natural Science Foundation for Young Scientists(No.41404093)+1 种基金Key National Research Project of China(Nos2016YFC0303100 and 2017YFC0601900)China Natural Science Foundation(No.41774125)
文摘The airborne electromagnetic (AEM) method has a high sampling rate and survey flexibility. However, traditional numerical modeling approaches must use high-resolution physical grids to guarantee modeling accuracy, especially for complex geological structures such as anisotropic earth. This can lead to huge computational costs. To solve this problem, we propose a spectral-element (SE) method for 3D AEM anisotropic modeling, which combines the advantages of spectral and finite-element methods. Thus, the SE method has accuracy as high as that of the spectral method and the ability to model complex geology inherited from the finite-element method. The SE method can improve the modeling accuracy within discrete grids and reduce the dependence of modeling results on the grids. This helps achieve high-accuracy anisotropic AEM modeling. We first introduced a rotating tensor of anisotropic conductivity to Maxwell's equations and described the electrical field via SE basis functions based on GLL interpolation polynomials. We used the Galerkin weighted residual method to establish the linear equation system for the SE method, and we took a vertical magnetic dipole as the transmission source for our AEM modeling. We then applied fourth-order SE calculations with coarse physical grids to check the accuracy of our modeling results against a 1D semi-analytical solution for an anisotropic half-space model and verified the high accuracy of the SE. Moreover, we conducted AEM modeling for different anisotropic 3D abnormal bodies using two physical grid scales and three orders of SE to obtain the convergence conditions for different anisotropic abnormal bodies. Finally, we studied the identification of anisotropy for single anisotropic abnormal bodies, anisotropic surrounding rock, and single anisotropic abnormal body embedded in an anisotropic surrounding rock. This approach will play a key role in the inversion and interpretation of AEM data collected in regions with anisotropic geology.
基金supported by the Doctoral Fund Project of the Ministry of Education(No.20130061110060 class tutors)the National Natural Science Foundation of China(No.41504083)National Basic Research Program of China(973Program)(No.2013CB429805)
文摘Frequency-domain airborne electromagnetics is a proven geophysical exploration method.Presently,the interpretation is mainly based on resistivity-depth imaging and onedimensional layered inversion;nevertheless,it is difficult to obtain satisfactory results for two- or three-dimensional complex earth structures using 1D methods.3D forward modeling and inversion can be used but are hampered by computational limitations because of the large number of data.Thus,we developed a 2.5D frequency-domain airborne electromagnetic forward modeling and inversion algorithm.To eliminate the source singularities in the numerical simulations,we split the fields into primary and secondary fields.The primary fields are calculated using homogeneous or layered models with analytical solutions,and the secondary(scattered) fields are solved by the finite-element method.The linear system of equations is solved by using the large-scale sparse matrix parallel direct solver,which greatly improves the computational efficiency.The inversion algorithm was based on damping leastsquares and singular value decomposition and combined the pseudo forward modeling and reciprocity principle to compute the Jacobian matrix.Synthetic and field data were used to test the effectiveness of the proposed method.
基金supported by the 863 Program (Grant no.2006AA09Z323)the 973 Program (Grant No.2006CB202402)
文摘Numerical simulation in the frequency-space domain has inherent advantages, such as: it is possible to simulate wave propagation from multiple sources simultaneously; there are no cumulative errors; only the interesting frequencies can be selected; and it is more suitable for wave propagation in viscoelastic media. The only obstacle to using the method is the requirement of huge computer storage. We extend the compressed format for storing the coefficient matrix. It can reduce the required computer storage dramatically. We get the optimal coefficients by least-squares method to suppress the numerical dispersion and adopt the perfectly matched layer (PML) boundary conditions to eliminate the artificial boundary reflections. Using larger grid intervals decreases computer storage requirements and provides high computational efficiency. Numerical experiments demonstrate that these means are economic and effective, providing a good basis for elastic wave imaging and inversion.
文摘Non-spreading nature of Bessel spatiotemporal wavepackets is theoretically and experimentally investigated and orders of magnitude improvement in the spatiotemporal spreading has been demonstrated.The spatiotemporal confinement provided by the Bessel spatiotemporal wavepacket is further exploited to transport transverse orbital angular momentum through embedding spatiotemporal optical vortex into the Bessel spatiotemporal wavepacket, constructing a new type of wavepacket: Bessel spatiotemporal optical vortex. Both numerical and experimental results demonstrate that spatiotemporal vortex structure can be well maintained and confined through much longer propagation. High order spatiotemporal optical vortices can also be better confined in the spatiotemporal domain and prevented from further breaking up, overcoming a potential major obstacle for future applications of spatiotemporal vortex.
基金supported by the National Natural Science Foundation of China (Grant Nos.10773002,10875012 and 11175019)the Team Research Program of Hubei University for Nationalities (Grant No.MY2011T006)Beijing Postdoctoral Research Foundation (Grant No.71006015201201)
文摘Motivated by the recent work that the periodicity of a black hole is responsible for the area spectrum,we exclusively utilize the period of motion of an outgoing wave,which is shown to be related to the vibrational frequency of the perturbed black hole,to study area spectra of a non-rotating BTZ black hole and a rotating BTZ black hole.It is found that the area spectra and entropy spectra for both space times are equally spaced.In addition,we find that though the entropy spectra of the 3-dimensional BTZ black holes take the same form as those of the 4-dimensional black holes,the area spectra depend on the dimension of space times.Our result confirms that the entropy spectrum of a black hole is more fundamental than the area spectrum.