Joint inversion method of 3D electrical resistivity detection based on inequality constraints
Liu Bin,Liu Zheng-Yu,Song Jie,Nie Li-Chao,Wang Chuan-Wu,Chen Lei
DOI: https://doi.org/10.6038/cjg20170232
2017-01-01
Abstract:Non-uniqueness is one inherent problem in 3D electrical resistivity linear and non-linear inversion. The linear inversion method has a rapid calculation speed, but its result is strongly dependent on the initial model, which may result in falling into the local minimum. The nonlinear inversion method is independent on the initial model, but its search and inversion efficiency is extremely low. And there is no literature about 3D resistivity nonlinear inversion at present. To solve these problems, based on the complementary strengths between linear and nonlinear inversion method, this study proposed a joint inversion method combining the least square method (a linear inversion method) and improved genetic algorithm (GA) method (a nonlinear inversion method). At first, the resistivity varying range is considered as an inequality constraint, which is incorporated into inversion equation. Then the definition method of looser inequality constraint and stricter inequality constraint is given, which is deduced from drilling holes. And the least square linear inversion method and GA nonlinear inversion method based on inequality constraints are put forward separately. Secondly, the GA method is improved in aspects of mutation direction control, production of initial population, by which searching direction, evolutionary efficiency and population diversity are improved greatly. Then, the joint inversion method and its implementation scheme are presented. The improved GA is used in the first inversion-stage. The advantages of its low dependence on the initial model are utilized. When the space near the optimized solution can be searched, the optimal individual should be exported. Then the least square method is used in the second inversion-stage. The optimal individual obtained by GA is taken as the initial model and efficient local linear search is executed in the space near the optimized solution. In this way, 3D imaging of geo-electric structure can be achieved. At last, verification on synthetic data and real examples has been carried out. Comparing with the traditional inversion method of least square shows that this joint inversion method has significant advantages in suppressing non-uniqueness, getting rid of the initial model dependence and improving inversion effect.