基于积分方程法的任意形状源多次激发三维电磁场反演
|
摘要
利用积分方程和阻尼最小二乘法,对任意形状源多个位置激发下三维异常体的电磁场进行了反演计算。计算中将任意形状源条件下多组激发点和接收点的电磁场数据统一考虑,求出了此种条件下的雅克比矩阵,反演得到了地下异常体的电阻率分布。模型试算结果表明,反演收敛速度快,对初值依赖性小,结果准确可靠。
3D abnormal body's electromagnetic field with arbitrary shape source multiple locations excitation can be inverted using integral equation and damped least squares.Multiple groups of electromagnetic field data in different excitation and receiving point to be uniform consideration in inversion,a matrix is obtained,distribution of resistivity of underground abnormal body is achieved.Model test shows that inversion of fast convergence speed,less dependent on the initial value,the result is accurate and reliable.
3D abnormal body's electromagnetic field with arbitrary shape source multiple locations excitation can be inverted using integral equation and damped least squares.Multiple groups of electromagnetic field data in different excitation and receiving point to be uniform consideration in inversion,a matrix is obtained,distribution of resistivity of underground abnormal body is achieved.Model test shows that inversion of fast convergence speed,less dependent on the initial value,the result is accurate and reliable.
引文
[1]HOHMANN G W.Three dimensional induced polarization and electromagnetic modeling[J].Geophysics1975,40(2):309-324.
[2]WANNAMAKER P E,HOHMANN G W,SANFILIPO W A.Electromagnetic modeling of three dimensional bodies in layered earths using integral equations[J].Geophysics 1984,49(1):60-74.
[3]朱伯芳.有限单元法原理与应用[M].北京:水利出版社,1979.
[4]徐世浙.地球物理中的有限单元法[M].北京:科学出版社,1994.
[5]周熙襄.电法勘探数值模拟技术[M].成都:四川科学技术出版社,1986.
[6]李建平,李桐林,张亚东.层状介质任意形状回线源瞬变电磁场正反演研究[J].物探与化探,2012,36(2):256-259.
[7]林昌洪,谭捍东,舒晴,等.可控源音频大地电磁三维共轭梯度反演研究[J].地球物理学报,2012,55(11):3829-3839.
[8]MACKIE R L,MADDEN T R.Three-dimensional magnetotelluric inversion using conjugate gradients[J].Geophys,J.Int,1993,115:215-229.
[9]NEWMAN G A,ALUMBAUGH D L.Three-dimensional magnetotelluric inversion using non-linear conjugate gradients[J].Geophys,J.Int,2000,140:410-424.
[10]胡祖志,胡祥云,何展翔.大地电磁非线性共轭梯度拟三维反演[J].地球物理学报,2006,49(4):1226-1234.
[11]LIN C H,TAN H D,TONG T.Three-dimension conjugate gradient inversion of magnetotelluric full information data[J].Applied Geophysics,2011,8(1):1-10.
[12]刘斌,李术才,李树忱,等.基于不等式约束的最小二乘法三维电阻率反演及其算法优化[J].地球物理学报,2012,55(1):260-268.
[13]宛新林,席道瑛,高尔根,等.用改进的光滑约束最小二乘正交分解法实现电阻率三维反演[J].地球物理学报,2005,48(2):439-444.
[14]SASAKI Y.3-D resistivity inversion using the finite-element method[J].Geophysics,1994,59(11):1839-1848.
[15]EATON P A.3Delectromagnetic inversion using integral equation[J].Geophysical prospecting,1989,37:407-426.
[16]PETRICK W R JR,SILL W R,WARD S H.Threedimensional resistivity inversion using alpha centers[J].Geophysics,1981,46(8):1148-1163.
[17]EL-QADY G,USHIJIMA K.Inversion of DC resistivity data using neural networks[J].Geophysical Prospecting,2001,49(4):417-430.
[2]WANNAMAKER P E,HOHMANN G W,SANFILIPO W A.Electromagnetic modeling of three dimensional bodies in layered earths using integral equations[J].Geophysics 1984,49(1):60-74.
[3]朱伯芳.有限单元法原理与应用[M].北京:水利出版社,1979.
[4]徐世浙.地球物理中的有限单元法[M].北京:科学出版社,1994.
[5]周熙襄.电法勘探数值模拟技术[M].成都:四川科学技术出版社,1986.
[6]李建平,李桐林,张亚东.层状介质任意形状回线源瞬变电磁场正反演研究[J].物探与化探,2012,36(2):256-259.
[7]林昌洪,谭捍东,舒晴,等.可控源音频大地电磁三维共轭梯度反演研究[J].地球物理学报,2012,55(11):3829-3839.
[8]MACKIE R L,MADDEN T R.Three-dimensional magnetotelluric inversion using conjugate gradients[J].Geophys,J.Int,1993,115:215-229.
[9]NEWMAN G A,ALUMBAUGH D L.Three-dimensional magnetotelluric inversion using non-linear conjugate gradients[J].Geophys,J.Int,2000,140:410-424.
[10]胡祖志,胡祥云,何展翔.大地电磁非线性共轭梯度拟三维反演[J].地球物理学报,2006,49(4):1226-1234.
[11]LIN C H,TAN H D,TONG T.Three-dimension conjugate gradient inversion of magnetotelluric full information data[J].Applied Geophysics,2011,8(1):1-10.
[12]刘斌,李术才,李树忱,等.基于不等式约束的最小二乘法三维电阻率反演及其算法优化[J].地球物理学报,2012,55(1):260-268.
[13]宛新林,席道瑛,高尔根,等.用改进的光滑约束最小二乘正交分解法实现电阻率三维反演[J].地球物理学报,2005,48(2):439-444.
[14]SASAKI Y.3-D resistivity inversion using the finite-element method[J].Geophysics,1994,59(11):1839-1848.
[15]EATON P A.3Delectromagnetic inversion using integral equation[J].Geophysical prospecting,1989,37:407-426.
[16]PETRICK W R JR,SILL W R,WARD S H.Threedimensional resistivity inversion using alpha centers[J].Geophysics,1981,46(8):1148-1163.
[17]EL-QADY G,USHIJIMA K.Inversion of DC resistivity data using neural networks[J].Geophysical Prospecting,2001,49(4):417-430.