This paper develops the boundary element method, the authors employ two-layered earth Green 's functions as the weighting functions of residual and derive boundary integral equations. The forward problems of point...This paper develops the boundary element method, the authors employ two-layered earth Green 's functions as the weighting functions of residual and derive boundary integral equations. The forward problems of point sources on 2 - D and 3-D structures with an influencing cover are solved by this method. The results show that this method markedly improves the original boundary element method. The features of the improved method are greater numerical accuracy and much smaller systems of equations and thus considerable savings for the storage capacity of computers, allowing us to solve the above problems with only ordinary microcomputers. The results in this paper extend the scope of applying the boundary element method while using electrical methods for geophysical prospecting.展开更多
As a means of quantitative interpretation,forward calculations of the global lithospheric magnetic field in the Spherical Harmonic(SH)domain have been widely used to reveal geophysical,lithological,and geothermal vari...As a means of quantitative interpretation,forward calculations of the global lithospheric magnetic field in the Spherical Harmonic(SH)domain have been widely used to reveal geophysical,lithological,and geothermal variations in the lithosphere.Traditional approaches either do not consider the non-axial dipolar terms of the inducing field and its radial variation or do so by means of complicated formulae.Moreover,existing methods treat the magnetic lithosphere either as an infinitesimally thin layer or as a radially uniform spherical shell of constant thickness.Here,we present alternative forward formulae that account for an arbitrarily high maximum degree of the inducing field and for a magnetic lithosphere of variable thickness.Our simulations based on these formulae suggest that the satellite magnetic anomaly field is sensitive to the non-axial dipolar terms of the inducing field but not to its radial variation.Therefore,in forward and inverse calculations of satellite magnetic anomaly data,the non-axial dipolar terms of the inducing field should not be ignored.Furthermore,our results show that the satellite magnetic anomaly field is sensitive to variability in the lateral thickness of the magnetized shell.In particular,we show that for a given vertically integrated susceptibility distribution,underestimating the thickness of the magnetic layer overestimates the induced magnetic field.This discovery bridges the greatest part of the alleged gap between the susceptibility values measured from rock samples and the susceptibility values required to match the observed magnetic field signal.We expect the formulae and conclusions of this study to be a valuable tool for the quantitative interpretation of the Earth's global lithospheric magnetic field,through an inverse or forward modelling approach.展开更多
On the basis of the results of improved analytical expression of computation of gravity anomalies due to a homogeneous polyhedral body composed of polygonal facets, and applying the forward theory with the coordinate ...On the basis of the results of improved analytical expression of computation of gravity anomalies due to a homogeneous polyhedral body composed of polygonal facets, and applying the forward theory with the coordinate transformation of vectors and tensors, we deduced both the analytical expressions for gravity gradient tensors and for magnetic anomalies of a polygon, and obtained new analytical expressions for computing vertical gradients of gravity anomalies and vertical component of magnetic anomalies caused by a polyhedral body. And also we developed explicitly the complete unified expressions for the calculation of gravity anomalies, gravity gradient, and magnetic anomalies due to the homogeneous polyhedron. Furthermore, we deduced new analytical expressions for computing vertical gradients of gravity anomalies due to a finite rectangular prism by applying the newly obtained expressions for gravity gradient tensors due to a polyhedral target body. Comparison with forward calculation of models shows the correctness of these new expressions. It will reduce forward calculation time of gravity-magnetic anomalies and improve computational efficiency by applying our unified expressions for joint forward modeling of gravity-magnetic anomalies due to homogeneous polyhedral bodies.展开更多
The pavement layered structures are composed of surface layer,road base and multi-layered soil foundation.They can be undermined over time by repeated vehicle loads.In this study,a hybrid numerical method which can ev...The pavement layered structures are composed of surface layer,road base and multi-layered soil foundation.They can be undermined over time by repeated vehicle loads.In this study,a hybrid numerical method which can evaluate the displacement responses of pavement structures under dynamic falling weight deflectometer(FWD)loads.The proposed method consists of two parts:(a)the dynamic stiffness matrices of the points at the surface in the frequency domain which is based on the domain-transformation and dual vector form equation,and(b)interpolates the dynamic stiffness matrices by a continues rational function of frequency.The mixed variables formulation(MVF)can treat multiple degree of freedom systems with considering the coupling term between degree of freedoms.The accuracy of the developed method has been demonstrated by comparison between the proposed method and published results from the other method.Then the proposed method can be applied as a forward calculation technique to emulate the falling weight deflectometer test for multi-layered pavement structures.展开更多
The electromagnetic induction method is widely used to measure sea ice thickness. Based on the electrical properties of sea ice and seawater, the method measures the apparent conductivity, which represents the conduct...The electromagnetic induction method is widely used to measure sea ice thickness. Based on the electrical properties of sea ice and seawater, the method measures the apparent conductivity, which represents the conductivity of the half-space, and calculates the thickness of the sea ice. During the fourth Chinese National Arctic Research Expedition in summer 2010, an integrated electromagnetic induction system was set up on the icebreaker R/V XUE LONG to measure sea ice thickness along the ship＇s tracks to the north of the Chukchi Sea. The conductivities of sea ice, seawater, and brine were measured and a simple forward model was used to explain the effect of changes in those conductivities on the apparent conductivity over a horizontal layered structure. The results of this analysis indicated that when using the electromagnetic induction method to measure sea ice thickness, the conductivity of sea ice could be neglected and the conductivity of seawater could be treated as a constant. The ice distribution results derived from the electromagnetic induction method showed that the typical sea ice thickness was 160 cm and 90 cm during the outbound and the return legs of the voyage, respectively.展开更多
In the interaction computation for 3D gravity and magnetic anomalies due to arbitrarily shaped homogenous magnetized polyhedron model composed of triangular facets, there are many difficult points, such as mass comput...In the interaction computation for 3D gravity and magnetic anomalies due to arbitrarily shaped homogenous magnetized polyhedron model composed of triangular facets, there are many difficult points, such as mass computing, absence of a mature computer technique in 3D geological body modeling, inconvenient human-computer interaction, hard program coding, etc.. Based on the formulae of the magnetic field due to horizontal regular bodies, and by applying forward theory with the three-dimensional Cartesian coordinate system transformation, the forward problems of magnetic anomalies and gradient tensors for arbitrary slantwise regular bodies were solved. It is shown that the magnetic calculating expressions of the arbitrary posture regular body are corrected by comparing results with the homogeneous polyhedral body model outcome data. Furthermore, in the same condition, the former significantly reduced forward time. Applying a new forward method of regular body expressions in arbitrary posture, developed software for interaction computation between the 3D geological body model and magnetic field has advantages of fast calculation speed, easy manipulation, etc..展开更多
Initial alignment is the precondition for strapdown inertial navigation system(SINS)to navigate.Its two important indexes are accuracy and rapidity,the accuracy of the initial alignment is directly related to the work...Initial alignment is the precondition for strapdown inertial navigation system(SINS)to navigate.Its two important indexes are accuracy and rapidity,the accuracy of the initial alignment is directly related to the working accuracy of SINS,but in self-alignment,the two indexes are often contradictory.In view of the limitations of conventional data processing algorithms,a novel method of compass alignment based on stored data and repeated navigation calculation for SINS is proposed.By means of data storage,the same data is used in different stages of the initial alignment,which is beneficial to shorten the initial alignment time and improve the alignment accuracy.In order to verify the correctness of the compass algorithm based on stored data and repeated navigation calculation,the simulation experiment was done.In summary,when the computer performance is sufficiently high,the compass alignment method based on the stored data and the forward and reverse navigation calculation can effectively improve the alignment speed and improve the alignment accuracy.展开更多
文摘This paper develops the boundary element method, the authors employ two-layered earth Green 's functions as the weighting functions of residual and derive boundary integral equations. The forward problems of point sources on 2 - D and 3-D structures with an influencing cover are solved by this method. The results show that this method markedly improves the original boundary element method. The features of the improved method are greater numerical accuracy and much smaller systems of equations and thus considerable savings for the storage capacity of computers, allowing us to solve the above problems with only ordinary microcomputers. The results in this paper extend the scope of applying the boundary element method while using electrical methods for geophysical prospecting.
基金supported by the National Natural Science Foundation of China(Grant Nos.42250103 and 42174090)the Opening Fund of Key Laboratory of Geological Survey and Evaluation of Ministry of Education(Grant No.GLAB2023ZR02)the Ministry of Science and Technology(MOST)Special Fund from the State Key Laboratory of Geological Processes and Mineral Resources(Grant No.MSFGPMR2022-4)。
文摘As a means of quantitative interpretation,forward calculations of the global lithospheric magnetic field in the Spherical Harmonic(SH)domain have been widely used to reveal geophysical,lithological,and geothermal variations in the lithosphere.Traditional approaches either do not consider the non-axial dipolar terms of the inducing field and its radial variation or do so by means of complicated formulae.Moreover,existing methods treat the magnetic lithosphere either as an infinitesimally thin layer or as a radially uniform spherical shell of constant thickness.Here,we present alternative forward formulae that account for an arbitrarily high maximum degree of the inducing field and for a magnetic lithosphere of variable thickness.Our simulations based on these formulae suggest that the satellite magnetic anomaly field is sensitive to the non-axial dipolar terms of the inducing field but not to its radial variation.Therefore,in forward and inverse calculations of satellite magnetic anomaly data,the non-axial dipolar terms of the inducing field should not be ignored.Furthermore,our results show that the satellite magnetic anomaly field is sensitive to variability in the lateral thickness of the magnetized shell.In particular,we show that for a given vertically integrated susceptibility distribution,underestimating the thickness of the magnetic layer overestimates the induced magnetic field.This discovery bridges the greatest part of the alleged gap between the susceptibility values measured from rock samples and the susceptibility values required to match the observed magnetic field signal.We expect the formulae and conclusions of this study to be a valuable tool for the quantitative interpretation of the Earth's global lithospheric magnetic field,through an inverse or forward modelling approach.
基金This paper is supported by the National Natural Science Foundation of China (No.40374039)Program for New Century Excellent Talents in University (No. NCET-04-0726)the Focused Subject Program of Beijing (No. XK104910598).
文摘On the basis of the results of improved analytical expression of computation of gravity anomalies due to a homogeneous polyhedral body composed of polygonal facets, and applying the forward theory with the coordinate transformation of vectors and tensors, we deduced both the analytical expressions for gravity gradient tensors and for magnetic anomalies of a polygon, and obtained new analytical expressions for computing vertical gradients of gravity anomalies and vertical component of magnetic anomalies caused by a polyhedral body. And also we developed explicitly the complete unified expressions for the calculation of gravity anomalies, gravity gradient, and magnetic anomalies due to the homogeneous polyhedron. Furthermore, we deduced new analytical expressions for computing vertical gradients of gravity anomalies due to a finite rectangular prism by applying the newly obtained expressions for gravity gradient tensors due to a polyhedral target body. Comparison with forward calculation of models shows the correctness of these new expressions. It will reduce forward calculation time of gravity-magnetic anomalies and improve computational efficiency by applying our unified expressions for joint forward modeling of gravity-magnetic anomalies due to homogeneous polyhedral bodies.
基金The authors are grateful for the financial support of the National Key Research and Development Program of China(No.2016YFC0802400)the National Natural Science Foundation of China under Grant No.(51508203,51678536,41404096)+2 种基金the Outstanding Young Talent Research Fund of Zhengzhou University(1621323001)Program for Innovative Research Team(in Science and Technology)in University of Henan Province(18IRTSTHN007)the Program for Science and Technology Innovation Talents in Universities of Henan Province(Grant No.19HASTIT043),and the authors extend their sincere gratitude.
文摘The pavement layered structures are composed of surface layer,road base and multi-layered soil foundation.They can be undermined over time by repeated vehicle loads.In this study,a hybrid numerical method which can evaluate the displacement responses of pavement structures under dynamic falling weight deflectometer(FWD)loads.The proposed method consists of two parts:(a)the dynamic stiffness matrices of the points at the surface in the frequency domain which is based on the domain-transformation and dual vector form equation,and(b)interpolates the dynamic stiffness matrices by a continues rational function of frequency.The mixed variables formulation(MVF)can treat multiple degree of freedom systems with considering the coupling term between degree of freedoms.The accuracy of the developed method has been demonstrated by comparison between the proposed method and published results from the other method.Then the proposed method can be applied as a forward calculation technique to emulate the falling weight deflectometer test for multi-layered pavement structures.
基金supported by the National Natural Science Foundation of China(Grant no.41006116)the National Major Scientific Research Project(Grant no.2013CBA01804)the Chinese Polar Environmental Comprehensive Investigation and Assessment Programs(Grant no.CHINARE-2015-02-02)
文摘The electromagnetic induction method is widely used to measure sea ice thickness. Based on the electrical properties of sea ice and seawater, the method measures the apparent conductivity, which represents the conductivity of the half-space, and calculates the thickness of the sea ice. During the fourth Chinese National Arctic Research Expedition in summer 2010, an integrated electromagnetic induction system was set up on the icebreaker R/V XUE LONG to measure sea ice thickness along the ship＇s tracks to the north of the Chukchi Sea. The conductivities of sea ice, seawater, and brine were measured and a simple forward model was used to explain the effect of changes in those conductivities on the apparent conductivity over a horizontal layered structure. The results of this analysis indicated that when using the electromagnetic induction method to measure sea ice thickness, the conductivity of sea ice could be neglected and the conductivity of seawater could be treated as a constant. The ice distribution results derived from the electromagnetic induction method showed that the typical sea ice thickness was 160 cm and 90 cm during the outbound and the return legs of the voyage, respectively.
基金supported by the National High Technology Research and Development Program of China(Nos. 2007AA06Z134, 2006AA06A201)the Program of In-troducing Talents of Discipline to Universities (No.B07011)Schlumberger-China University of Geosciences Education Fund (No.SLBX0807)
文摘In the interaction computation for 3D gravity and magnetic anomalies due to arbitrarily shaped homogenous magnetized polyhedron model composed of triangular facets, there are many difficult points, such as mass computing, absence of a mature computer technique in 3D geological body modeling, inconvenient human-computer interaction, hard program coding, etc.. Based on the formulae of the magnetic field due to horizontal regular bodies, and by applying forward theory with the three-dimensional Cartesian coordinate system transformation, the forward problems of magnetic anomalies and gradient tensors for arbitrary slantwise regular bodies were solved. It is shown that the magnetic calculating expressions of the arbitrary posture regular body are corrected by comparing results with the homogeneous polyhedral body model outcome data. Furthermore, in the same condition, the former significantly reduced forward time. Applying a new forward method of regular body expressions in arbitrary posture, developed software for interaction computation between the 3D geological body model and magnetic field has advantages of fast calculation speed, easy manipulation, etc..
基金This work was supported by the National Nature Science Foundation of China(Grant No.5200110367)Natural Science Foundation of Jiangsu Province(Grant No.SBK2020043219)+1 种基金Scientific Research Foundation of the Higher Education Institutions of Jiangsu Province(Grant No.19KJB510052)NUPTSF(Grant No.NY219023).
文摘Initial alignment is the precondition for strapdown inertial navigation system(SINS)to navigate.Its two important indexes are accuracy and rapidity,the accuracy of the initial alignment is directly related to the working accuracy of SINS,but in self-alignment,the two indexes are often contradictory.In view of the limitations of conventional data processing algorithms,a novel method of compass alignment based on stored data and repeated navigation calculation for SINS is proposed.By means of data storage,the same data is used in different stages of the initial alignment,which is beneficial to shorten the initial alignment time and improve the alignment accuracy.In order to verify the correctness of the compass algorithm based on stored data and repeated navigation calculation,the simulation experiment was done.In summary,when the computer performance is sufficiently high,the compass alignment method based on the stored data and the forward and reverse navigation calculation can effectively improve the alignment speed and improve the alignment accuracy.
