Traditional two-dimensional（2D） complex resistivity forward modeling is based on Poisson＇s equation but spectral induced polarization（SIP） data are the coproducts of the induced polarization（IP） and the electro...Traditional two-dimensional（2D） complex resistivity forward modeling is based on Poisson＇s equation but spectral induced polarization（SIP） data are the coproducts of the induced polarization（IP） and the electromagnetic induction（EMI） effects.This is especially true under high frequencies,where the EMI effect can exceed the IP effect.2D inversion that only considers the IP effect reduces the reliability of the inversion data.In this paper,we derive differential equations using Maxwell＇s equations.With the introduction of the Cole-Cole model,we use the finite-element method to conduct2 D SIP forward modeling that considers the EMI and IP effects simultaneously.The data-space Occam method,in which different constraints to the model smoothness and parametric boundaries are introduced,is then used to simultaneously obtain the four parameters of the Cole-Cole model using multi-array electric field data.This approach not only improves the stability of the inversion but also significantly reduces the solution ambiguity.To improve the computational efficiency,message passing interface programming was used to accelerate the 2D SIP forward modeling and inversion.Synthetic datasets were tested using both serial and parallel algorithms,and the tests suggest that the proposed parallel algorithm is robust and efficient.展开更多
Due to the induced polarization(IP)eff ect,the sign reversal often occurs in timedomain airborne electromagnetic(AEM)data.The inversions that do not consider IP eff ect cannot recover the true umderground electrical s...Due to the induced polarization(IP)eff ect,the sign reversal often occurs in timedomain airborne electromagnetic(AEM)data.The inversions that do not consider IP eff ect cannot recover the true umderground electrical structures.In view of the fact that there are many parameters of airborne induced polarization data in time domain,and the sensitivity diff erence between parameters is large,which brings challenges to the stability and accuracy of the inversion.In this paper,we propose an inversion mehtod for time-domain AEM data with IP effect based on the Pearson correlation constraints.This method uses the Pearson correlation coeffi cient in statistics to characterize the correlation between the resistivity and the chargeability and constructs the Pearson correlation constraints for inverting the objective function to reduce the non uniqueness of inversion.To verify the eff ectiveness of this method,we perform both Occam’s inversion and Pearson correlation constrained inversion on the synthetic data.The experiments show that the Pearson correlation constrained inverison is more accurate and stable than the Occam’s inversion.Finally,we carried out the inversion to a survey dataset with and without IP eff ect.The results show that the data misfit and the continuity of the inverted section are greatly improved when the IP eff ect is considered.展开更多
The study of induced polarization （IP） information extraction from magnetotelluric （MT） sounding data is of great and practical significance to the exploitation of deep mineral, oil and gas resources. The linear i...The study of induced polarization （IP） information extraction from magnetotelluric （MT） sounding data is of great and practical significance to the exploitation of deep mineral, oil and gas resources. The linear inversion method, which has been given priority in previous research on the IP information extraction method, has three main problems as follows： 1） dependency on the initial model, 2） easily falling into the local minimum, and 3） serious non-uniqueness of solutions. Taking the nonlinearity and nonconvexity of IP information extraction into consideration, a two-stage CO-PSO minimum structure inversion method using compute unified distributed architecture （CUDA） is proposed. On one hand, a novel Cauchy oscillation particle swarm optimization （CO-PSO） algorithm is applied to extract nonlinear IP information from MT sounding data, which is implemented as a parallel algorithm within CUDA computing architecture; on the other hand, the impact of the polarizability on the observation data is strengthened by introducing a second stage inversion process, and the regularization parameter is applied in the fitness function of PSO algorithm to solve the problem of multi-solution in inversion. The inversion simulation results of polarization layers in different strata of various geoelectric models show that the smooth models of resistivity and IP parameters can be obtained by the proposed algorithm, the results of which are relatively stable and accurate. The experiment results added with noise indicate that this method is robust to Gaussian white noise. Compared with the traditional PSO and GA algorithm, the proposed algorithm has more efficiency and better inversion results.展开更多
基金jointly sponsored by the National Natural Science Foundation of China(Grant No.41374078)the Geological Survey Projects of the Ministry of Land and Resources of China(Grant Nos.12120113086100 and 12120113101300)Beijing Higher Education Young Elite Teacher Project
文摘Traditional two-dimensional（2D） complex resistivity forward modeling is based on Poisson＇s equation but spectral induced polarization（SIP） data are the coproducts of the induced polarization（IP） and the electromagnetic induction（EMI） effects.This is especially true under high frequencies,where the EMI effect can exceed the IP effect.2D inversion that only considers the IP effect reduces the reliability of the inversion data.In this paper,we derive differential equations using Maxwell＇s equations.With the introduction of the Cole-Cole model,we use the finite-element method to conduct2 D SIP forward modeling that considers the EMI and IP effects simultaneously.The data-space Occam method,in which different constraints to the model smoothness and parametric boundaries are introduced,is then used to simultaneously obtain the four parameters of the Cole-Cole model using multi-array electric field data.This approach not only improves the stability of the inversion but also significantly reduces the solution ambiguity.To improve the computational efficiency,message passing interface programming was used to accelerate the 2D SIP forward modeling and inversion.Synthetic datasets were tested using both serial and parallel algorithms,and the tests suggest that the proposed parallel algorithm is robust and efficient.
基金This paper was fi nancially supported by the National Natural Science Foundation of China(Nos.42030806,41774125,41904104,41804098)the Pioneer Project of Chinese Academy of Sciences(No.XDA14020102).
文摘Due to the induced polarization(IP)eff ect,the sign reversal often occurs in timedomain airborne electromagnetic(AEM)data.The inversions that do not consider IP eff ect cannot recover the true umderground electrical structures.In view of the fact that there are many parameters of airborne induced polarization data in time domain,and the sensitivity diff erence between parameters is large,which brings challenges to the stability and accuracy of the inversion.In this paper,we propose an inversion mehtod for time-domain AEM data with IP effect based on the Pearson correlation constraints.This method uses the Pearson correlation coeffi cient in statistics to characterize the correlation between the resistivity and the chargeability and constructs the Pearson correlation constraints for inverting the objective function to reduce the non uniqueness of inversion.To verify the eff ectiveness of this method,we perform both Occam’s inversion and Pearson correlation constrained inversion on the synthetic data.The experiments show that the Pearson correlation constrained inverison is more accurate and stable than the Occam’s inversion.Finally,we carried out the inversion to a survey dataset with and without IP eff ect.The results show that the data misfit and the continuity of the inverted section are greatly improved when the IP eff ect is considered.
基金Projects(41604117,41204054)supported by the National Natural Science Foundation of ChinaProjects(20110490149,2015M580700)supported by the Research Fund for the Doctoral Program of Higher Education,China+1 种基金Project(2015zzts064)supported by the Fundamental Research Funds for the Central Universities,ChinaProject(16B147)supported by the Scientific Research Fund of Hunan Provincial Education Department,China
文摘The study of induced polarization （IP） information extraction from magnetotelluric （MT） sounding data is of great and practical significance to the exploitation of deep mineral, oil and gas resources. The linear inversion method, which has been given priority in previous research on the IP information extraction method, has three main problems as follows： 1） dependency on the initial model, 2） easily falling into the local minimum, and 3） serious non-uniqueness of solutions. Taking the nonlinearity and nonconvexity of IP information extraction into consideration, a two-stage CO-PSO minimum structure inversion method using compute unified distributed architecture （CUDA） is proposed. On one hand, a novel Cauchy oscillation particle swarm optimization （CO-PSO） algorithm is applied to extract nonlinear IP information from MT sounding data, which is implemented as a parallel algorithm within CUDA computing architecture; on the other hand, the impact of the polarizability on the observation data is strengthened by introducing a second stage inversion process, and the regularization parameter is applied in the fitness function of PSO algorithm to solve the problem of multi-solution in inversion. The inversion simulation results of polarization layers in different strata of various geoelectric models show that the smooth models of resistivity and IP parameters can be obtained by the proposed algorithm, the results of which are relatively stable and accurate. The experiment results added with noise indicate that this method is robust to Gaussian white noise. Compared with the traditional PSO and GA algorithm, the proposed algorithm has more efficiency and better inversion results.
