重力梯度欧拉反褶积及其在文顿盐丘的应用
|
摘要
在矿产勘察与区域地质调查中,梯度数据由于精度高而得到广泛应用。基于三种张量分量的重力梯度数据,本文提出一种新的梯度欧拉反褶积算法,避免了常规计算中由于构造指数选取不当和张量分量转换产生的误差。建立单长方体和不同深度双长方体的理论模型,验证了该方法能够较准确地反演地质体的埋深和水平位置。将此方法应用于文顿盐丘的实测重力梯度数据,获得了盖岩顶部的形态与分布等信息,并探测到新的小规模地质体。该结果与已知资料相符,进一步丰富了该地区的地下构造信息。
In the mineral exploration and regional geological survey works,the gradiometry data is widely used because of its high precision.A new gradiometry Euler deconvolution is proposed in this paper combining three kinds of tensors of gravity gradiometry data,which can avoid errors from improperly selecting structural index and tensor transformation in calculations.Building theoretical models of single prism and dual prisms in different depths verifies that the proposed method can accurately inverse geological body depth and horizontal location.This method is applied to gradiometry data in Vinton Dome.The shape and distribution of cap rock top are depicted and new small-scale geological bodies are detected.The results agree very well with the known data,further enriches the underground structure information of this area.
In the mineral exploration and regional geological survey works,the gradiometry data is widely used because of its high precision.A new gradiometry Euler deconvolution is proposed in this paper combining three kinds of tensors of gravity gradiometry data,which can avoid errors from improperly selecting structural index and tensor transformation in calculations.Building theoretical models of single prism and dual prisms in different depths verifies that the proposed method can accurately inverse geological body depth and horizontal location.This method is applied to gradiometry data in Vinton Dome.The shape and distribution of cap rock top are depicted and new small-scale geological bodies are detected.The results agree very well with the known data,further enriches the underground structure information of this area.
引文
[1] 耿美霞,杨庆节.应用RBF神经网络反演二维重力密度分布[J].石油地球物理勘探,2013,48(4):651-657.GENG Meixia,YANG Qingjie.2-D density inversion with the RBF neural network method[J].Oil Geophysical Prospecting,2013,48(4):651-657.
[2] 林宝泽,肖锋,王明常.二维重力数据径向反演及应用[J].石油地球物理勘探,2018,53(2):403-409.LIN Baoze,XIAO Feng,WANG Mingchang.2D gra-vity data radial inversion[J].Oil Geophysical Prospecting,2018,53(2):403-409.
[3] Thompson D T.EULDPH:A new technique for ma-king computer-assisted depth estimates from magnetic data[J].Geophysics,1982,47(1):31-37.
[4] 周文纳,杜晓娟,李吉焱.重力数据的平面全张量梯度角度边界识别方法[J].石油地球物理勘探,2013,48(6):1009-1015.ZHOU Wenna,DU Xiaojuan,LI Jiyan.Gravity data edge identification using plane full tensor gradient angle[J].Oil Geophysical Prospecting,2013,48(6):1009-1015.
[5] 王彦国,吴姿颖,邓居智,等.基于幂次平均的离散归一化总梯度法[J].石油地球物理勘探,2018,53(6):1351-1364.WANG Yanguo,WU Ziying,DENG Juzhi,et al.Discretization of normalized total gradient method based on power mean[J].Oil Geophysical Prospecting,2018,53(6):1351-1364.
[6] Zhang C,Mushayandebvu M F,Reid A B,et al.Euler deconvolution of gravity tensor gradient data[J].Geophysics,2000,65(2):512-520.
[7] Ma G.Combination of horizontal gradient ratio and Euler methods for the interpretation of potential field data[J].Geophysics,2013,78(5):J53-J60.
[8] Ma G,Huang D,Liu C.Application of balanced edge detection filters to estimate the location parameters of the causative sources using potential field data[J].Journal of Applied Geophysics,2013,99:18-23.
[9] 马国庆,杜晓娟,李丽丽,等.梯度反褶积法及其在矿产勘探中的应用[J].吉林大学学报:地球科学版,2013,43 (1):259-266.MA Guoqing,DU Xiaojuan,LI Lili,et al.Gradient deconvolution method and its application in mineral exploration[J].Journal of Jilin University (Earth Science Edition),2013,43(1):259-266.
[10] Guo C C,Xiong S Q,Xue D J,et al.Improved Euler method for the interpretation of potential data based on the ratio of the vertical first derivative to analytic signal[J].Applied Geophysics,2014,11(3):331-339.
[11] Ma G.The application of extended Euler deconvolution method in the interpretation of potential field data[J].Journal of Applied Geophysics,2014,107:188-194.
[12] Zhou W,Nan Z,Li J.Self-constrained Euler deconvolution using potential field data of different altitudes[J].Pure & Applied Geophysics,2016,173(6):2073-2085.
[13] 周文月.重力及梯度数据联合欧拉反褶积方法研究[D].吉林长春:吉林大学,2016.
[14] Coker M O,Bhattacharya J P,Marfurt K J.Fracture patterns within mudstones on the flanks of a salt dome:syneresis or slumping?[C].Gulf Coast Association of Geological Societies,2007,57:125-137.
[15] Ennen C,Hall S.Structural mapping of the Vinton salt dome,Louisiana,using gravity gradiometry data[C].SEG Technical Program Expanded Abstracts,2011,30:830-835.
[16] Oliveira Jr V C,Barbosa V C F.3-D radial gravity gradient inversion[J].Geophysical Journal International,2013,195(2):883-902.
[17] Hou Z,Huang D,Wei X.Fast inversion of probability tomography with gravity gradiometry data based on hybrid parallel programming[J].Journal of Applied Geophysics,2016,124:27-38.
[18] Hou Z,Wei X,Huang D.Fast density inversion solution for full tensor gravity gradiometry data[J].Pure & Applied Geophysics,2016,173(2):509-523.
[19] 秦朋波,黄大年.重力和重力梯度数据联合聚焦反演方法[J].地球物理学报,2016,59(6):2203-2224.QIN Pengbo,HUANG Danian.Integrated gravity and gravity gradient data focusing inversion[J].Chinese Journal of Geophysics,2016,59(6):2203-2224.
[20] 王浩然,陈超,杜劲松.重力梯度张量数据的三维反演方法与应用[J].石油地球物理勘探,2013,48(3):474-481.WANG Haoran,CHEN Chao,DU Jinsong.3-D inversion of gravity gradient tensor data and its applications[J].Oil Geophysical Prospecting,2013,48(3):474-481.
[2] 林宝泽,肖锋,王明常.二维重力数据径向反演及应用[J].石油地球物理勘探,2018,53(2):403-409.LIN Baoze,XIAO Feng,WANG Mingchang.2D gra-vity data radial inversion[J].Oil Geophysical Prospecting,2018,53(2):403-409.
[3] Thompson D T.EULDPH:A new technique for ma-king computer-assisted depth estimates from magnetic data[J].Geophysics,1982,47(1):31-37.
[4] 周文纳,杜晓娟,李吉焱.重力数据的平面全张量梯度角度边界识别方法[J].石油地球物理勘探,2013,48(6):1009-1015.ZHOU Wenna,DU Xiaojuan,LI Jiyan.Gravity data edge identification using plane full tensor gradient angle[J].Oil Geophysical Prospecting,2013,48(6):1009-1015.
[5] 王彦国,吴姿颖,邓居智,等.基于幂次平均的离散归一化总梯度法[J].石油地球物理勘探,2018,53(6):1351-1364.WANG Yanguo,WU Ziying,DENG Juzhi,et al.Discretization of normalized total gradient method based on power mean[J].Oil Geophysical Prospecting,2018,53(6):1351-1364.
[6] Zhang C,Mushayandebvu M F,Reid A B,et al.Euler deconvolution of gravity tensor gradient data[J].Geophysics,2000,65(2):512-520.
[7] Ma G.Combination of horizontal gradient ratio and Euler methods for the interpretation of potential field data[J].Geophysics,2013,78(5):J53-J60.
[8] Ma G,Huang D,Liu C.Application of balanced edge detection filters to estimate the location parameters of the causative sources using potential field data[J].Journal of Applied Geophysics,2013,99:18-23.
[9] 马国庆,杜晓娟,李丽丽,等.梯度反褶积法及其在矿产勘探中的应用[J].吉林大学学报:地球科学版,2013,43 (1):259-266.MA Guoqing,DU Xiaojuan,LI Lili,et al.Gradient deconvolution method and its application in mineral exploration[J].Journal of Jilin University (Earth Science Edition),2013,43(1):259-266.
[10] Guo C C,Xiong S Q,Xue D J,et al.Improved Euler method for the interpretation of potential data based on the ratio of the vertical first derivative to analytic signal[J].Applied Geophysics,2014,11(3):331-339.
[11] Ma G.The application of extended Euler deconvolution method in the interpretation of potential field data[J].Journal of Applied Geophysics,2014,107:188-194.
[12] Zhou W,Nan Z,Li J.Self-constrained Euler deconvolution using potential field data of different altitudes[J].Pure & Applied Geophysics,2016,173(6):2073-2085.
[13] 周文月.重力及梯度数据联合欧拉反褶积方法研究[D].吉林长春:吉林大学,2016.
[14] Coker M O,Bhattacharya J P,Marfurt K J.Fracture patterns within mudstones on the flanks of a salt dome:syneresis or slumping?[C].Gulf Coast Association of Geological Societies,2007,57:125-137.
[15] Ennen C,Hall S.Structural mapping of the Vinton salt dome,Louisiana,using gravity gradiometry data[C].SEG Technical Program Expanded Abstracts,2011,30:830-835.
[16] Oliveira Jr V C,Barbosa V C F.3-D radial gravity gradient inversion[J].Geophysical Journal International,2013,195(2):883-902.
[17] Hou Z,Huang D,Wei X.Fast inversion of probability tomography with gravity gradiometry data based on hybrid parallel programming[J].Journal of Applied Geophysics,2016,124:27-38.
[18] Hou Z,Wei X,Huang D.Fast density inversion solution for full tensor gravity gradiometry data[J].Pure & Applied Geophysics,2016,173(2):509-523.
[19] 秦朋波,黄大年.重力和重力梯度数据联合聚焦反演方法[J].地球物理学报,2016,59(6):2203-2224.QIN Pengbo,HUANG Danian.Integrated gravity and gravity gradient data focusing inversion[J].Chinese Journal of Geophysics,2016,59(6):2203-2224.
[20] 王浩然,陈超,杜劲松.重力梯度张量数据的三维反演方法与应用[J].石油地球物理勘探,2013,48(3):474-481.WANG Haoran,CHEN Chao,DU Jinsong.3-D inversion of gravity gradient tensor data and its applications[J].Oil Geophysical Prospecting,2013,48(3):474-481.