基于松弛变量的不等式约束三维电阻率并行反演方法研究

doi:10.13722/j.cnki.jrme.2014.1731

Abstract
Figure/Table
References (0)
Related Citation (15)

Download: PDF (12653 KB) HTML (1 KB)
Export: BibTeX | EndNote (RIS)

Abstract Multiple solutions may be resulted from the inversion of 3D electrical resistivity and the calculation is time-consuming，which restricts the application of resistivity detection. An idea is presented to solve this problem，with the inequality constraints to the inversion equations and the parallel algorithms being applied. Based on the traditional smooth constraints，the relaxed variables were introduced to apply the inequality constraints to the inverse equations，which carry the upper and lower limits information of subsurface media. A new objective function was obtained using the primal dual interior point method. The resistivity values in the function were defined within the inequality constrains range，and the parameters were optimized within the feasible region defined by the constraints. The method theoretically suppressed the multiplicity of inverse solution. The parallel calculation algorithm for partial derivative matrix and the parallel algorithm for Cholesky decomposition to the overall coefficient matrix were designed，which speeds the inversion more than fifty percent. Based on the above research，the parallel 3D electrical resistivity inversion method with inequality constraint based on relax variables was developed，and numerical tests and engineering application were carried out. Results showed that the above method made full use of the inequality constrains and removed the false anomalies，suppressed the multiplicity and improved the accuracy and calculation efficiency.

Key words： numerical simulation 3D electrical resistivity inversion inequality constraint slack variables primal dual interior point method parallel inversion

	Service

	E-mail this article
	Add to my bookshelf
	Add to citation manager
	E-mail Alert
	RSS
	Articles by authors
	LI Shucai
	WANG Chuanwu
	NIE Lichao
	LIU Bin
	CHEN Lei
	LIU Zhengyu
	TIAN Mingzhen

Cite this article:

LI Shucai,WANG Chuanwu,NIE Lichao, et al. Parallel 3D electrical resistivity inversion method with#br# inequality constraint based on slack variables[J]. , 2016, 35(6): 1122-1132.

URL:

https://rockmech.whrsm.ac.cn/EN/10.13722/j.cnki.jrme.2014.1731 OR https://rockmech.whrsm.ac.cn/EN/Y2016/V35/I6/1122

[1] SASAKI Y. 3D resistivity inversion using the finite element method[J]. Geophysics，1994，59(11)：1 839–1 848.
[2] 黄俊革，阮百尧，鲍光淑. 基于有限单元法的三维地电断面电阻率反演[J]. 中南大学学报：自然科学版，2004，35(2)：295–299. (HUANG Junge，RUAN Baiyao，BAO Guangshu. Resistivity inversion on 3-D section based on FEM[J]. Journal of Central South University：Science and Technology，2004，35(2)：295–299.(in Chinese))
[3] 刘斌，李术才，聂利超，等. 基于自适应加权光滑约束与PCG算法的三维电阻率探测反演成像[J]. 岩土工程学报，2012，34(9)：1 646–1 653.(LIU Bin，LI Shucai，NIE Lichao，et al. Study on inversion method of 3D resistivity detection using adaptive-weighted smooth constraint and PCG algorithm[J]. Chinese Journal of Geotechnical Engineering，2012，34(9)：1 646–1 653.(in Chinese))
[4] OLDENBURG D W，LI Y. Inversion of induced polarization data[J]. Geophysics，1994，59(9)：1 327–1 341.
[5] LI Y，OLDENBURG D W. 3D inversion of induced polarization data[J]. Geophysics，2000，65(6)：1 931–1 945.
[6] KAIPIO J P，KOLEHMAINEN V，VAUKONEN M，et al. Inverse problems with structural prior information[J]. Inverse Problems，1999，15(3)：713–729.
[7] 刘斌，聂利超，李术才，等. 三维电阻率空间结构约束反演成像方法[J]. 岩石力学与工程学报，2012，31(11)：2 258–2 268.(LIU Bin，NIE Lichao，LI Shucai，et al. 3D electrical resistivity inversion tomography with spatial structure constraint[J]. Chinese Journal of Rock Mechanics and Engineering，2012，31(11)：2 258–2 268.(in Chinese))
[8] 刘斌，李术才，聂利超，等. 矿井突水灾变过程电阻率约束反演成像实时监测模拟研究[J]. 煤炭学报，2012，37(10)：1 722–1 731. (LIU Bin，LI Shucai，NIE Lichao，et al. Research on simulation of mine water inrush real-time monitoring of using electrical resistivity constrained inversion imaging method[J]. Journal of China Coal Society，2012，37(10)：1 722–1 731.(in Chinese))
[9] 刘斌，李术才，李树忱，等. 基于不等式约束的最小二乘法三维电阻率反演及其算法优化[J]. 地球物理学报，2012，55(1)：260–268. (LIU Bin，LI Shucai，LI Shuchen，et al. 3D electrical resistivity inversion with least-squares method based on inequality constraint and its computation efficiency optimization[J]. Chinese Journal of Geophysics，2012，55(1)：260–268.(in Chinese))
[10] WEI H，SASAKI H，KUBOKAWA J，et al. An interior point nonlinear programming for optimal power flow problems with a novel data structure[J]. IEEE Transactions on Power Systems，1998，13(3)：870–877.
[11] 范宏，韦化. 基于扰动 KKT 条件的原始–对偶内点法和分支定界法的最优潮流研究[J]. 电力自动化设备，2004，24(5)：5–9. (FAN Hong，WEI Hua. Study on optimal power flow based on primal-dual interior point algorithm under perturbed KKT conditions and branch-and-bound method[J]. Electric Power Automation Equipment，2004，24(5)：5–9.(in Chinese))
[12] LIN C G，TAN H D，TONG T. Parallel rapid relaxation of 3D magneto telluric data[J]. Applied Geophysics，2009，6(1)：77–83.
[13] 陈召曦，孟小红，郭良辉，等. 基于GPU并行的重力、重力梯度三维正演快速计算及反演策略[J]. 地球物理学报，2012，55(12)：4 069–4 077.(CHEN Zhaoxi，MENG Xiaohong，GUO Lianghui，et al. Three dimension fast forward modeling and the inversion strategy for large scale gravity and gravimetry data based on GPU[J]. Chinese Journal of Geophysics，2012，55(12)：4 069–4 077.(in Chinese))
[14] 刘劲松，刘福田，刘俊，等. 地震层析成像LSQR算法的并行化[J]. 地球物理学报，2006，49(2)：540–545.(LIU Jinsong，LIU Futian， LIU Jun，et al. Parallel LSQR algorithms used in seismic tomography[J]. Chinese Journal of Geophysics，2006，49(2)：540–545.(in Chinese))