Traditionally, airborne time-domain electromagnetic （ATEM） data are inverted to derive the earth model by iteration. However, the data are often highly correlated among channels and consequently cause ill-posed and ...Traditionally, airborne time-domain electromagnetic （ATEM） data are inverted to derive the earth model by iteration. However, the data are often highly correlated among channels and consequently cause ill-posed and over-determined problems in the inversion. The correlation complicates the mapping relation between the ATEM data and the earth parameters and thus increases the inversion complexity. To obviate this, we adopt principal component analysis to transform ATEM data into orthogonal principal components （PCs） to reduce the correlations and the data dimensionality and simultaneously suppress the unrelated noise. In this paper, we use an artificial neural network （ANN） to approach the PCs mapping relation with the earth model parameters, avoiding the calculation of Jacobian derivatives. The PC-based ANN algorithm is applied to synthetic data for layered models compared with data-based ANN for airborne time-domain electromagnetic inversion. The results demonstrate the PC-based ANN advantages of simpler network structure, less training steps, and better inversion results over data-based ANN, especially for contaminated data. Furthermore, the PC-based ANN algorithm effectiveness is examined by the inversion of the pseudo 2D model and comparison with data-based ANN and Zhody＇s methods. The results indicate that PC-based ANN inversion can achieve a better agreement with the true model and also proved that PC-based ANN is feasible to invert large ATEM datasets.展开更多
Electrically anisotropic strata are abundant in nature, so their study can help our data interpretation and our understanding of the processes of geodynamics. However, current data processing generally assumes isotrop...Electrically anisotropic strata are abundant in nature, so their study can help our data interpretation and our understanding of the processes of geodynamics. However, current data processing generally assumes isotropic conditions when surveying anisotropic structures, which may cause discrepancies between reality and electromagnetic data interpretation. Moreover, the anisotropic interpretation of the time-domain airborne electromagnetic （TDAEM） method is still confined to one dimensional （1D） cases, and the corresponding three-dimensional （3D） numerical simulations are still in development. In this study, we expanded the 3D TDAEM modeling of arbitrarily anisotropic media. First, through coordinate rotation of isotropic conductivity, we obtained the conductivity tensor of an arbitrary anisotropic rock. Next, we incorporated this into Maxwell＇s equations, using a regular hexahedral grid of vector finite elements to subdivide the solution area. A direct solver software package provided the solution for the sparse linear equations that resulted. Analytical solutions were used to verify the accuracy and feasibility of the algorithm. The proven model was then applied to analyze the effects of arbitrary anisotropy in 3D TDAEM via the distribution of responses and amplitude changes, which revealed that different anisotropy situations strongly affected the responses of TDAEM.展开更多
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.展开更多
Airborne electromagnetic transient method enjoys the advantages of high-efficiency and the high resolution of electromagnetic anomalies,especially suitable for mining detection around goaf areas and deep exploration o...Airborne electromagnetic transient method enjoys the advantages of high-efficiency and the high resolution of electromagnetic anomalies,especially suitable for mining detection around goaf areas and deep exploration of minerals.In this paper,we calculated the full-wave airborne transient electromagnetic data,according to the result of numerical research,the advantage of switch-off time response in electromagnetic detection was proofed via experiments.Firstly,based on the full-wave airborne transient electromagnetic system developed by Jilin University(JLU-ATEMI),we proposed a method to compute the full-waveform electromagnetic(EM)data of 3D model using the FDTD approach and convolution algorithm,and verify the calculation by the response of homogenous half-space.Then,through comparison of switch-off-time response and off-time response,we studied the effect of ramp time on anomaly detection.Finally,we arranged two experimental electromagnetic detection,the results indicated that the switch-off-time response can reveal the shallow target more effectively,and the full-waveform airborne electromagnetic system is an effective technique for shallow target detection.展开更多
We develop a new computational method for modeling and inverting frequency domain airborne electromagnetic(EM)data.Our method is based on the contraction integral equation method for forward EM modeling and on inversi...We develop a new computational method for modeling and inverting frequency domain airborne electromagnetic(EM)data.Our method is based on the contraction integral equation method for forward EM modeling and on inversion using the localized quasi-linear(LQL)approximation followed by the rigorous inversion,if necessary.The LQL inversion serves to provide a fast image of the target.These results are checked by a rigorous update of the domain electric field,allowing a more accurate calculation of the predicted data.If the accuracy is poorer than desired,rigorous inversion follows,using the resulting conductivity distribution and electric field from LQL as a starting model.The rigorous inversion iteratively solves the field and domain equations,converting the non-linear inversion into a series of linear inversions.We test this method on synthetic and field data.The results of the inversion are very encouraging with respect to both the speed and the accuracy of the algorithm,showing this is a useful tool for airborne EM interpretation.展开更多
基金supported by the National Natural Science Foundation of China (Grant No. 40974039)High-Tech Research and Development Program of China (Grant No.2006AA06205)Leading Strategic Project of Science and Technology, Chinese Academy of Sciences (XDA08020500)
文摘Traditionally, airborne time-domain electromagnetic （ATEM） data are inverted to derive the earth model by iteration. However, the data are often highly correlated among channels and consequently cause ill-posed and over-determined problems in the inversion. The correlation complicates the mapping relation between the ATEM data and the earth parameters and thus increases the inversion complexity. To obviate this, we adopt principal component analysis to transform ATEM data into orthogonal principal components （PCs） to reduce the correlations and the data dimensionality and simultaneously suppress the unrelated noise. In this paper, we use an artificial neural network （ANN） to approach the PCs mapping relation with the earth model parameters, avoiding the calculation of Jacobian derivatives. The PC-based ANN algorithm is applied to synthetic data for layered models compared with data-based ANN for airborne time-domain electromagnetic inversion. The results demonstrate the PC-based ANN advantages of simpler network structure, less training steps, and better inversion results over data-based ANN, especially for contaminated data. Furthermore, the PC-based ANN algorithm effectiveness is examined by the inversion of the pseudo 2D model and comparison with data-based ANN and Zhody＇s methods. The results indicate that PC-based ANN inversion can achieve a better agreement with the true model and also proved that PC-based ANN is feasible to invert large ATEM datasets.
基金financially supported by National Nonprofit institute Research Grant of IGGE(Nos.AS2017J06,AS2017Y04,and AS2016J10)Survey on coastal area for airborne magnetic method of UNV in Jiangsu(No.DD20160151-03)+3 种基金Key National Research Project of China(No.2017YFC0601900)Key Program of National Natural Science Foundation of China(No.41530320)Natural Science Foundation(No.41274121)China Natural Science Foundation for Young Scientists(No.41404093)
文摘Electrically anisotropic strata are abundant in nature, so their study can help our data interpretation and our understanding of the processes of geodynamics. However, current data processing generally assumes isotropic conditions when surveying anisotropic structures, which may cause discrepancies between reality and electromagnetic data interpretation. Moreover, the anisotropic interpretation of the time-domain airborne electromagnetic （TDAEM） method is still confined to one dimensional （1D） cases, and the corresponding three-dimensional （3D） numerical simulations are still in development. In this study, we expanded the 3D TDAEM modeling of arbitrarily anisotropic media. First, through coordinate rotation of isotropic conductivity, we obtained the conductivity tensor of an arbitrary anisotropic rock. Next, we incorporated this into Maxwell＇s equations, using a regular hexahedral grid of vector finite elements to subdivide the solution area. A direct solver software package provided the solution for the sparse linear equations that resulted. Analytical solutions were used to verify the accuracy and feasibility of the algorithm. The proven model was then applied to analyze the effects of arbitrary anisotropy in 3D TDAEM via the distribution of responses and amplitude changes, which revealed that different anisotropy situations strongly affected the responses of TDAEM.
基金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.
基金Project(41674109) supported by the National Natural Science Foundation of China
文摘Airborne electromagnetic transient method enjoys the advantages of high-efficiency and the high resolution of electromagnetic anomalies,especially suitable for mining detection around goaf areas and deep exploration of minerals.In this paper,we calculated the full-wave airborne transient electromagnetic data,according to the result of numerical research,the advantage of switch-off time response in electromagnetic detection was proofed via experiments.Firstly,based on the full-wave airborne transient electromagnetic system developed by Jilin University(JLU-ATEMI),we proposed a method to compute the full-waveform electromagnetic(EM)data of 3D model using the FDTD approach and convolution algorithm,and verify the calculation by the response of homogenous half-space.Then,through comparison of switch-off-time response and off-time response,we studied the effect of ramp time on anomaly detection.Finally,we arranged two experimental electromagnetic detection,the results indicated that the switch-off-time response can reveal the shallow target more effectively,and the full-waveform airborne electromagnetic system is an effective technique for shallow target detection.
文摘We develop a new computational method for modeling and inverting frequency domain airborne electromagnetic(EM)data.Our method is based on the contraction integral equation method for forward EM modeling and on inversion using the localized quasi-linear(LQL)approximation followed by the rigorous inversion,if necessary.The LQL inversion serves to provide a fast image of the target.These results are checked by a rigorous update of the domain electric field,allowing a more accurate calculation of the predicted data.If the accuracy is poorer than desired,rigorous inversion follows,using the resulting conductivity distribution and electric field from LQL as a starting model.The rigorous inversion iteratively solves the field and domain equations,converting the non-linear inversion into a series of linear inversions.We test this method on synthetic and field data.The results of the inversion are very encouraging with respect to both the speed and the accuracy of the algorithm,showing this is a useful tool for airborne EM interpretation.
