基于地质体空间位置优化约束的航空重力梯度数据三维物性反演
|
摘要
重力数据的物性反演面临着严重的多解性问题,降低多解性的有效手段是加入约束条件.而边界识别、深度估计及成像方法可获取地质体的水平位置、深度范围等几何参数信息,本文将基于数据本身挖掘的地质体几何参数信息约束到物性反演中,以降低反演的多解性.通过引入基于深度信息的深度加权函数及基于水平位置的水平梯度加权函数建立优化约束条件,有效地提高了反演结果的横向及纵向分辨率.重力梯度数据包含更多的地质体空间特征信息,将优化约束反演方法应用到全张量数据的反演中,模型试验表明本文方法反演结果与理论模型更加吻合.最后对美国路易斯安那州文顿盐丘实测航空重力梯度数据的应用表明,本文方法在其他地球物理、地质资料不足的情况下获得更可靠的反演结果.
The inversion of gravity data for physical properties faces a serious non-uniqueness problem.The effective means to reduce such non-uniqueness is to add more constraints.The edge detecting,depth estimation and imaging methods can obtain the geometric parameters such as the horizontal position and depth range of the geological body.In this paper,the geological body geometric parameter information based on the data itself is constrained to the physical property inversion to reduce the non-uniqueness of inversion.By introducing the depth information-based depth weighting function and the horizontal position-based horizontal gradient weighting functionto establish optimizing constraints,the horizontal and vertical resolution of inversion results are effectively improved.The gravity gradient data contains more information on spatial features of the geological body,and the optimizing constrained inversion method is applied to the inversion of full tensor data.Tests on a synthetic model show that the inversion results of this method are more consistent with the theoretical model.Finally,the application of this method to the airborne gravity gradient data in the Vinton Salt in Louisiana,USA shows that this approach permits to obtain more reliable inversion results in the case with insufficient other geophysical and geological data.
The inversion of gravity data for physical properties faces a serious non-uniqueness problem.The effective means to reduce such non-uniqueness is to add more constraints.The edge detecting,depth estimation and imaging methods can obtain the geometric parameters such as the horizontal position and depth range of the geological body.In this paper,the geological body geometric parameter information based on the data itself is constrained to the physical property inversion to reduce the non-uniqueness of inversion.By introducing the depth information-based depth weighting function and the horizontal position-based horizontal gradient weighting functionto establish optimizing constraints,the horizontal and vertical resolution of inversion results are effectively improved.The gravity gradient data contains more information on spatial features of the geological body,and the optimizing constrained inversion method is applied to the inversion of full tensor data.Tests on a synthetic model show that the inversion results of this method are more consistent with the theoretical model.Finally,the application of this method to the airborne gravity gradient data in the Vinton Salt in Louisiana,USA shows that this approach permits to obtain more reliable inversion results in the case with insufficient other geophysical and geological data.
引文
Barbosa V C F,Silva J B C,Medeiros W E.1999.Stable inversion of gravity anomalies of sedimentary basins with nonsmooth basement reliefs and arbitrary density contrast variations.Geophysics,64(3):754-764.
Carter-McAuslan A,Lelièvre P G,Farquharson C G.2015.Astudy of fuzzy c-means coupling for joint inversion,using seismic tomography and gravity data test scenarios.Geophysics,80(1):W1-W15.
Cella F,Fedi M.2012.Inversion of potential field data using the structural index as weighting function rate decay.Geophysical Prospecting,60(2):313-336.
Commer M.2011.Three-dimensional gravity modelling and focusing inversion using rectangular meshes.Geophysical Prospecting,59(5):966-979.
Commer M,Newman G A,Williams K H,et al.2011.3Dinducedpolarization data inversion for complex resistivity.Geophysics,76(3):F157-F171.
Fedi M,Abbas M A.2013.A fast interpretation of self-potential data using the depth from extreme points method.Geophysics,78(2):E107-E116.
Fregoso E,Gallardo L A,García-Abdeslem J.2015.Structural joint inversion coupled with Euler deconvolution of isolated gravity and magnetic anomalies.Geophysics,80(2):G67-G79.
Geng M X,Huang D N,Yang Q J,et al.2014.3Dinversion of airborne gravity-gradiometry data using cokriging.Geophysics,79(4):G37-G47.
Giraud J,Pakyuz-Charrier E,Jessell M,et al.2017.Uncertainty reduction through geologically conditioned petrophysical constraints in joint inversion.Geophysics,82(6):ID19-ID34.
Guillen A,Menichetti V.1984.Gravity and magnetic inversion with minimization of a specific functional.Geophysics,49(8):1354-1360.
Guo M,Deng J Z,Chen X,et al.2018.Bayesian joint inversion of magnetotelluric and gravity based on petrophysical Constraints.Progress in Geophysics(in Chinese),33(5):1897-1902,doi:10.6038/pg2018BB0461.
Ialongo S,Fedi M,Florio G.2014.Invariant models in the inversion of gravity and magnetic fields and their derivatives.Journal of Applied Geophysics,110:51-62.
Last B J,Kubik K.1983.Compact gravity inversion.Geophysics,48(6):713-721.
Lelièvre P G,Oldenburg D W.2009.A comprehensive study of including structural orientation information in geophysical inversions.Geophysical Journal International,178(2):623-637.
Li Y G,Oldenburg D W.1996.3-D inversion of magnetic data.Geophysics,61(2):394-408.
Li Y G,Oldenburg D W.1998.3-D inversion of gravity data.Geophysics,63(1):109-119.
Liu Y P,Wang Z W,Du X J,et al.2013.3Dconstrained inversion of gravity data based on the Extrapolation Tikhonov regularization algorithm.Chinese Journal of Geophysics(in Chinese),56(5):1650-1659.
Lo Re D,Florio G,Ferranti L,et al.2016.Self-constrained inversion of microgravity data along a segment of the Irpinia fault.Journal of Applied Geophysics,124:148-154.
Newman G A,Alumbaugh D L.2000.Three-dimensional magnetotelluric inversion using non-linear conjugate gradients.Geophysical Journal International,140(2):410-424.
Oliveira V C Jr,Barbosa V C F.2013.3-D radial gravity gradient inversion.Geophysical Journal International,195(2):883-902.
Paoletti V,Ialongo S,Florio G,et al.2013.Self-constrained inversion of potential fields.Geophysical Journal International,195(2):854-869.
Portniaguine O,Zhdanov M S.1999.Focusing geophysical inversion images.Geophysics,64(3):874-887.
Silva J B C,Medeiros W E,Barbosa V C F.2001.Potential-field inversion:choosing the appropriate technique to solve a geologic problem.Geophysics,66(2):511-520.
Sun S D,Chen C.2016.A self-constrained inversion of magnetic data based on correlation method.Journal of Applied Geophysics,135:8-16.
Yin C C,Sun S Y,Gao X H,et al.2018.3Djoint inversion of magnetotelluric and gravity data based on local correlation constraints.Chinese Journal of Geophysics,61(1):358-367,doi:10.6038/cjg2018K0765.
Zhang Q,Yu P,Zhang D L,et al.2018.Preliminary inversion of the potential field model based on DEXP.Progress in Geophysics(in Chinese),33(5):2076-2082,doi:10.6038/pg2018BB0558.
Zhang Y,Yan J G,Li F,et al.2015.A new bound constraints method for 3-D potential field data inversion using Lagrangian multipliers.Geophysical Journal International,201(1):267-275.
Zhou S,Huang D N.2015.Determining the depth of certain gravity sources without a priori specification of their structural index.Journal of Applied Geophysics,122:172-180.
Zhou S,Huang D N,Jiao J.2016.A filter to detect edge of potential field data based on three-dimensional structural tensors.Chinese Journal of Geophysics(in Chinese),59(10):3847-3858,doi:10.6038/cjg20161028.
Zhou S,Jiang J B,Lu P Y.2017.Source parameters estimation by the normalized downward continuation of gravity gradient data.Journal of Applied Geophysics,147:81-90.
郭曼,邓居智,陈晓,等.2018.基于岩石物性约束的大地电磁与重力贝叶斯联合反演.地球物理学进展,33(5):1897-1902,doi:10.6038/pg2018BB0461.
刘银萍,王祝文,杜晓娟等.2013.基于Extrapolation Tikhonov正则化算法的重力数据三维约束反演.地球物理学报,56(5):1650-1659.
殷长春,孙思源,高秀鹤等.2018.基于局部相关性约束的三维大地电磁数据和重力数据的联合反演.地球物理学报,61(1):358-367,doi:10.6038/cjg2018K0765.
张琦,于平,张代磊,等.2018.基于DEXP的位场模型初步反演.地球物理学进展,33(5):2076-2082,doi:10.6038/pg2018BB0558.
周帅,黄大年,焦健.2016.基于三维构造张量的位场边界识别滤波器.地球物理学报,59(10):3847-3858.
Carter-McAuslan A,Lelièvre P G,Farquharson C G.2015.Astudy of fuzzy c-means coupling for joint inversion,using seismic tomography and gravity data test scenarios.Geophysics,80(1):W1-W15.
Cella F,Fedi M.2012.Inversion of potential field data using the structural index as weighting function rate decay.Geophysical Prospecting,60(2):313-336.
Commer M.2011.Three-dimensional gravity modelling and focusing inversion using rectangular meshes.Geophysical Prospecting,59(5):966-979.
Commer M,Newman G A,Williams K H,et al.2011.3Dinducedpolarization data inversion for complex resistivity.Geophysics,76(3):F157-F171.
Fedi M,Abbas M A.2013.A fast interpretation of self-potential data using the depth from extreme points method.Geophysics,78(2):E107-E116.
Fregoso E,Gallardo L A,García-Abdeslem J.2015.Structural joint inversion coupled with Euler deconvolution of isolated gravity and magnetic anomalies.Geophysics,80(2):G67-G79.
Geng M X,Huang D N,Yang Q J,et al.2014.3Dinversion of airborne gravity-gradiometry data using cokriging.Geophysics,79(4):G37-G47.
Giraud J,Pakyuz-Charrier E,Jessell M,et al.2017.Uncertainty reduction through geologically conditioned petrophysical constraints in joint inversion.Geophysics,82(6):ID19-ID34.
Guillen A,Menichetti V.1984.Gravity and magnetic inversion with minimization of a specific functional.Geophysics,49(8):1354-1360.
Guo M,Deng J Z,Chen X,et al.2018.Bayesian joint inversion of magnetotelluric and gravity based on petrophysical Constraints.Progress in Geophysics(in Chinese),33(5):1897-1902,doi:10.6038/pg2018BB0461.
Ialongo S,Fedi M,Florio G.2014.Invariant models in the inversion of gravity and magnetic fields and their derivatives.Journal of Applied Geophysics,110:51-62.
Last B J,Kubik K.1983.Compact gravity inversion.Geophysics,48(6):713-721.
Lelièvre P G,Oldenburg D W.2009.A comprehensive study of including structural orientation information in geophysical inversions.Geophysical Journal International,178(2):623-637.
Li Y G,Oldenburg D W.1996.3-D inversion of magnetic data.Geophysics,61(2):394-408.
Li Y G,Oldenburg D W.1998.3-D inversion of gravity data.Geophysics,63(1):109-119.
Liu Y P,Wang Z W,Du X J,et al.2013.3Dconstrained inversion of gravity data based on the Extrapolation Tikhonov regularization algorithm.Chinese Journal of Geophysics(in Chinese),56(5):1650-1659.
Lo Re D,Florio G,Ferranti L,et al.2016.Self-constrained inversion of microgravity data along a segment of the Irpinia fault.Journal of Applied Geophysics,124:148-154.
Newman G A,Alumbaugh D L.2000.Three-dimensional magnetotelluric inversion using non-linear conjugate gradients.Geophysical Journal International,140(2):410-424.
Oliveira V C Jr,Barbosa V C F.2013.3-D radial gravity gradient inversion.Geophysical Journal International,195(2):883-902.
Paoletti V,Ialongo S,Florio G,et al.2013.Self-constrained inversion of potential fields.Geophysical Journal International,195(2):854-869.
Portniaguine O,Zhdanov M S.1999.Focusing geophysical inversion images.Geophysics,64(3):874-887.
Silva J B C,Medeiros W E,Barbosa V C F.2001.Potential-field inversion:choosing the appropriate technique to solve a geologic problem.Geophysics,66(2):511-520.
Sun S D,Chen C.2016.A self-constrained inversion of magnetic data based on correlation method.Journal of Applied Geophysics,135:8-16.
Yin C C,Sun S Y,Gao X H,et al.2018.3Djoint inversion of magnetotelluric and gravity data based on local correlation constraints.Chinese Journal of Geophysics,61(1):358-367,doi:10.6038/cjg2018K0765.
Zhang Q,Yu P,Zhang D L,et al.2018.Preliminary inversion of the potential field model based on DEXP.Progress in Geophysics(in Chinese),33(5):2076-2082,doi:10.6038/pg2018BB0558.
Zhang Y,Yan J G,Li F,et al.2015.A new bound constraints method for 3-D potential field data inversion using Lagrangian multipliers.Geophysical Journal International,201(1):267-275.
Zhou S,Huang D N.2015.Determining the depth of certain gravity sources without a priori specification of their structural index.Journal of Applied Geophysics,122:172-180.
Zhou S,Huang D N,Jiao J.2016.A filter to detect edge of potential field data based on three-dimensional structural tensors.Chinese Journal of Geophysics(in Chinese),59(10):3847-3858,doi:10.6038/cjg20161028.
Zhou S,Jiang J B,Lu P Y.2017.Source parameters estimation by the normalized downward continuation of gravity gradient data.Journal of Applied Geophysics,147:81-90.
郭曼,邓居智,陈晓,等.2018.基于岩石物性约束的大地电磁与重力贝叶斯联合反演.地球物理学进展,33(5):1897-1902,doi:10.6038/pg2018BB0461.
刘银萍,王祝文,杜晓娟等.2013.基于Extrapolation Tikhonov正则化算法的重力数据三维约束反演.地球物理学报,56(5):1650-1659.
殷长春,孙思源,高秀鹤等.2018.基于局部相关性约束的三维大地电磁数据和重力数据的联合反演.地球物理学报,61(1):358-367,doi:10.6038/cjg2018K0765.
张琦,于平,张代磊,等.2018.基于DEXP的位场模型初步反演.地球物理学进展,33(5):2076-2082,doi:10.6038/pg2018BB0558.
周帅,黄大年,焦健.2016.基于三维构造张量的位场边界识别滤波器.地球物理学报,59(10):3847-3858.