A two-stage CO-PSO minimum structure inversion using CUDA for extracting IP information from MT data

Abstract

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.

Keywords

Cauchy oscillation particle swarm optimization / magnetotelluric sounding / nonlinear inversion / induced polarization (IP) information extraction / compute unified distributed architecture (CUDA)

References

[1] HeJ-shanDual frequency induced polarization method [M], 2006, Beijing, Higher Education Press

[2] WuH-r, WangS-ming. The feasibility of nature source induced polarization [J]. Geophysical and Geochemical Exploration, 1978, 1: 62-64

[3] WareGTheoretical and field investigations of telluric currents and induced polarization [D], 1974, Berkeley, University of California at Berkeley

[4] MuraliS, RaoI B R, BhimasankaramV L S. Comparison of anomalous effects determined using telluric fields and time domain IP technique (test results) [J]. Exploration Geophysics, 1980, 11(2): 45-46

[5] GasperikovaE, MorrisonH F. Mapping of induced polarization using natural fields [J]. Geophysics, 2001, 66(1): 137-147

[6] GhorbaniA, CamerlynckC, FlorschN. CR1Dinv: A matlab program to invert 1D spectral induced polarization data for the Cole–Cole model including electromagnetic effects [J]. Computers & Geosciences, 2009, 35(2): 255-266

[7] YueA-p, DiQ-y, WangM-yue. 1-D forward modeling of the CSAMT signal incorporating IP effect [J]. Chinese Journal of Geophysics, 2009, 52(7): 1937-1946

[8] HeZ-x, HuZ-z, LuoW-f, WangC-fu. Mapping reservoirs based on resistivity and induced polarization derived from continuous 3D magnetotelluric profiling: Case study from Qaidam basin, China [J]. Geophysics, 2010, 75(1): B25-B33

[9] YuC-t, LiuH-f, ZhangX-jun. The analysis on IP signals in TEM response based on SVD [J]. Applied Geophysics, 2013, 10(1): 79-87

[10] TangR, XuP, XiangY, ZhangXu. The sensitivity analysis of different induced polarization models used in magnetotelluric method [J]. Acta Geodaetica et Geophysica, 2014, 49(2): 225-233

[11] ShawR, SrivastavaS. Particle swarm optimization: A new tool to invert geophysical data [J]. Geophysics, 2007, 72(2): F75-F83

[12] YuanS-y, WangS-x, TianNan. Swarm intelligence optimization and its application in geophysical data inversion [J]. Applied Geophysics, 2009, 6(2): 166-174

[13] JiangF-b, DaiQ-w, DongLi. An ICPSO-RBFNN nonlinear inversion for electrical resistivity imaging [J]. Journal of Central South University, 2016, 23(8): 2129-2138

[14] DiasC AA non-grounded method for measuring electrical induced polarization and conductivity [D], 1968, Berkeley, University of California

[15] DiasC A. Developments in a model to describe lowfrequency electrical polarization of rocks [J]. Geophysics, 2000, 65(2): 437-451

[16] KennedyJ, EberhartR. Particle swarm optimization [C]. IEEE International Conference on Neural Networks, 199519421948

[17] Nabavi-KeriziS H, AbadiM, KabirE. A PSO-based weighting method for linear combination of neural networks [J]. Computers & Electrical Engineering, 2010, 36(5): 886-894

[18] MussiL, DaolioF, CagnoniS. Evaluation of parallel particle swarm optimization algorithms within the CUDA™ architecture [J]. Information Sciences, 2011, 181(20): 4642-4657

[19] OuyangA-j, TangZ, ZhouX, XuY-m, PanG, LiK-qin. Parallel hybrid PSO with CUDA for lD heat conduction equation [J]. Computers & Fluids, 2015, 110(30): 198-210

[20] DaliN, BouamamaS. GPU-PSO: Parallel particle swarm optimization approaches on graphical processing unit for constraint reasoning: case of Max-CSPs [J]. Procedia Computer Science, 2015, 60(1): 1070-1080

[21] UgolottiR, NashedY S G, MesejoP L, MussiL, CagoniS. Particle swarm optimization and differential evolution for model-based object detection [J]. Applied Soft Computing, 2013, 13(6): 3092-3105

[22] JiangF-b, DaiQ-w, DongLi. Ultra-high density resistivity nonlinear inversion based on principal component-regularized ELM [J]. Chinese Journal of Geophysics-Chinese Edition, 2015, 58(9): 3356-3369

[23] ConstableS C, ParkerR L, ConstableC G. Occam’s inversion: A practical algorithm for generating smooth models from electromagnetic sounding data [J]. Geophysics, 1987, 52(3): 289-300

[24] JiangF-b, DaiQ-w, DongLi. Nonlinear inversion of electrical resistivity imaging using pruning Bayesian neural networks [J]. Applied Geophysics, 2016, 13(2): 267-278

[25] LiuM-l, LiuX-n, WuM-xin. Integrating spectral indices with environmental parameters for estimating heavy metal concentrations in rice using a dynamic fuzzy neural-network model [J]. Computers & Geosciences, 2011, 37(10): 1642-1652

[26] JiangF-b, DongL, DaiQ-w, DavidC N. Using wavelet packet denoising and ANFIS networks based on COSFLA optimization for electrical resistivity imaging inversion [J]. Fuzzy Sets and Systems, 2018, 337: 93-112

[27] LiG, TangJ-t, XiaoX, LiJ, ZhuH-j, ZhouC, YanF-bao. Near-source noise suppression of AMT by compressive sensing and mathematical morphology filtering [J]. Applied Geophysics, 2017, 14(4): 581-589

[28] JiangF-b, DongL, DaiQ-wei. Electrical resistivity imaging inversion: An ISFLA trained kernel principal component wavelet neural network approach [J]. Neural Networks, 2018