概率成像技术在重力位场中的应用研究
|
摘要
重力数据处理反演解释方法技术长期以来都存在许多问题,而概率成像技术在重力勘探中能够解决复杂地质问题,且方法简单,相对容易实现,运算速度较快,计算解释结果相对稳定准确,在部分地区亦不需要很多条件参数约束,同时更不需要在重力反演中给出初始模型以及拟合迭代的具体过程,只需要把滑动窗口内的相关函数扫描点逐点进行计算即可,适合在大规模数据成像中使用。通过研究布格异常与基于位场分离的梯度成像,给出了多种组合的球体模型成像结果,成像方法可显示出异常地质体的空间分布形态,提高重力数据横向与垂向分辨率。
Inversion interpretation of gravity data processing has many problems for a long time.However,probability tomography is simple,fast-calculating and easy to implement,and it can solve complex geological problems occurred in gravity exploration.The results of calculation and interpretation are accurate.It does not require lots of conditional parameter constraints in some regions,and initial model and related processes of fitting iterations are not required in gravity inversion.We just need to calculate the scanned points of correlation function within the sliding window one by one.This method is quite useful in large-scale data imaging.This paper studied the Bouguer anomaly and gradient imaging based on potential field separation and achieved the imaging results of various combinations of sphere models.The imaging method explained the spatial relationship of abnormal geological bodies and improved the lateral and vertical resolution of gravity data.
Inversion interpretation of gravity data processing has many problems for a long time.However,probability tomography is simple,fast-calculating and easy to implement,and it can solve complex geological problems occurred in gravity exploration.The results of calculation and interpretation are accurate.It does not require lots of conditional parameter constraints in some regions,and initial model and related processes of fitting iterations are not required in gravity inversion.We just need to calculate the scanned points of correlation function within the sliding window one by one.This method is quite useful in large-scale data imaging.This paper studied the Bouguer anomaly and gradient imaging based on potential field separation and achieved the imaging results of various combinations of sphere models.The imaging method explained the spatial relationship of abnormal geological bodies and improved the lateral and vertical resolution of gravity data.
引文
[1]OLDENBURG D W.The inversion and interpretation of gravity anomalies[J].Geophysics,1974,39(4):526-536.
[2]GREEN W R.Inversion of Gravity Profiles by Use of a Backus-Gilbert Approach[J].Geophysics,1975,40(5):763-772.
[3]LAST B J,KUBIK K.Compact gravity inversion[J].Geophysics,1983,48(48):713-721.
[4]管志宁,侯俊胜,黄临平,等.重磁异常反演的拟BP神经网络方法及其应用[J].地球物理学报,1998(02):242-251.
[5]ZHDANOV M S,ELLIS R,MUKHERJEE S.Three-dimensional regularized focusing inversion of gravity gradient tensor component data[J].Geophysics,2004,69(4):925-937.
[6]VASCO D W.Inversion of airborne gravity gradient data,southwestern Oklahoma[J].Geophysics,1991,56(1):90-101.
[7]PORTNIAGUINE O,ZHDANOV M S.Focusing geophysical inversion images[J].Geophysics,1999,64(3):874-887.
[8]FEDI M,RAPOLLA A.3-D inversion of gravity and magnetic data with depth resolution[J].Geophysics,1999,64(2):452-460.
[9]BEAR G W.Linear inversion of gravity data for 3-D density distributions[J].Geophysics,1995,60(5):1354-1364.
[10]BARBOSA V C F.Generalized compact gravity inversion[J].Geophysics,1994,7(59):57-68.
[11]MAURIELLO P,PATELLA D.Localization of magnetic sources underground by a data adaptive tomographic scanner[J/OL].Physics,2005.[2018-03-20].DOI:10.2528/PIERM08050504.
[12]MAURIELLO P,PATELLA D.Imaging polar and dipolar sources of geophysical anomalies by probability tomography.Part II:application to the vesuvius volcanic area[J/OL].Physics,2006.[2018-03-20].http://www.docin.com/p-1694540523.html.
[13]ALAIA R,PATELLA D,MAURIELLO P.Imaging multipole gravity anomaly sources by 3D probability tomography[J].Journal of Geophysics&Engineering,2009,6(3):298-310.
[14]IULIANO T,MAURIELLO P,PATELLA D.Looking inside Mount Vesuvius by potential fields integrated probability tomographies[J].Journal of Volcanology&Geothermal Research,2002,113(3):363-378.
[15]CHEN Z,MENG X,GUO L,et al.GICUDA:A parallel program for 3Dcorrelation imaging of large scale gravity and gravity gradiometry data on graphics processing units with CUDA[J].Computers&Geosciences,2012,46(3):119-128.
[16]孟小红,刘国峰,陈召曦,等.基于剩余异常相关成像的重磁物性反演方法[J].地球物理学报,2012,55(1):304-309.
[17]郭良辉,孟小红,石磊,等.重力和重力梯度数据三维[J].地球物理学报,2009,52(4),1098-1106.
[18]许令周.概率成像方法对台阶模型的成像[J].西南民族大学学报(自然科学版),2007,33(2):364-367.
[19]许令周.重力场中概率成像方法对脉状矿体的识别[J].内蒙古师范大学学报(自然科学汉文版),2007(4):463-465,469.
[20]马国庆,杜晓娟,李丽丽.改进的位场相关成像方法[J].地球科学(中国地质大学学报),2013,38(5):1121-1127.
[21]PATELLA D.Introduction to ground surface self-potential tomography[J].Geophysical Prospecting,1997a,45,653-681.
[22]PATELLA D.Self-potential global tomography including topographic effects[J].Geophysical Prospecting,1997b,45,843-863.
[23]MAURIELLO P,MONNA D,PATELLA D.3D geoelectric tomography and archaeological applications[J].Geophysical Prospecting,1998,46:543-570.
[24]MAURIELLO P,PATELLA D.Resistivity anomaly imaging by probability tomography[J].Geophysical Prospecting,1999a,47:411-429.
[25]MAURIELLO P,PATELLA D.Principles of probability tomography for natural-source electromagnetic induction fields,Geophysics,1999b,64,1403-1417.
[26]阮帅,王绪本,高永才,等.位场正则化下延在自然电位概率成像中的应用[J].物探化探计算技术,2006(4):332-336,285.
[27]王绪本,高永才.基于正则化向下延拓的重磁概率成像方法研究[C]//2008年重磁数据处理解释应用研讨会论文集.成都:成都理工大学,2008:49-54.
[28]黄金明.三维局部重磁场源全方位成像可视化研究[D].北京:中国地质大学(北京),2003.
[29]MAURIELLO P,PATELLA D.Gravity probability tomography:a new tool for buried mass distribution imaging[J].Geophysical Prospecting,2001,49:1-12.
[30]吉日嘎拉图,鲁光银.基于线性回归的重力数据三维RMS成像[J].物探化探计算技术,2014,36(2):185-189.
[31]杨旭明,王克斌.一种层析成像的正演和反演方法[J].物探化探计算技术,2002,24(4):337-340.
[32]闫浩飞,刘国峰.一种重力异常概率成像的扩展计算[J].地球物理学进展,2014,29(4):1837-1842.
[33]徐龙强,刘志义.物探化探数据成像预处理[J].物探化探计算技术,1991,13(1):41-46.
[34]路利春,文博,程斌,等.成像技术在矿区重力勘探中的应用研究[J].地质与资源,2018(1):99-102,88.
[35]许令周,关继腾,房文静.减小概率成像分辨宽度的一种方法[J].物探化探计算技术,2005,27(2):138-140,93.
[2]GREEN W R.Inversion of Gravity Profiles by Use of a Backus-Gilbert Approach[J].Geophysics,1975,40(5):763-772.
[3]LAST B J,KUBIK K.Compact gravity inversion[J].Geophysics,1983,48(48):713-721.
[4]管志宁,侯俊胜,黄临平,等.重磁异常反演的拟BP神经网络方法及其应用[J].地球物理学报,1998(02):242-251.
[5]ZHDANOV M S,ELLIS R,MUKHERJEE S.Three-dimensional regularized focusing inversion of gravity gradient tensor component data[J].Geophysics,2004,69(4):925-937.
[6]VASCO D W.Inversion of airborne gravity gradient data,southwestern Oklahoma[J].Geophysics,1991,56(1):90-101.
[7]PORTNIAGUINE O,ZHDANOV M S.Focusing geophysical inversion images[J].Geophysics,1999,64(3):874-887.
[8]FEDI M,RAPOLLA A.3-D inversion of gravity and magnetic data with depth resolution[J].Geophysics,1999,64(2):452-460.
[9]BEAR G W.Linear inversion of gravity data for 3-D density distributions[J].Geophysics,1995,60(5):1354-1364.
[10]BARBOSA V C F.Generalized compact gravity inversion[J].Geophysics,1994,7(59):57-68.
[11]MAURIELLO P,PATELLA D.Localization of magnetic sources underground by a data adaptive tomographic scanner[J/OL].Physics,2005.[2018-03-20].DOI:10.2528/PIERM08050504.
[12]MAURIELLO P,PATELLA D.Imaging polar and dipolar sources of geophysical anomalies by probability tomography.Part II:application to the vesuvius volcanic area[J/OL].Physics,2006.[2018-03-20].http://www.docin.com/p-1694540523.html.
[13]ALAIA R,PATELLA D,MAURIELLO P.Imaging multipole gravity anomaly sources by 3D probability tomography[J].Journal of Geophysics&Engineering,2009,6(3):298-310.
[14]IULIANO T,MAURIELLO P,PATELLA D.Looking inside Mount Vesuvius by potential fields integrated probability tomographies[J].Journal of Volcanology&Geothermal Research,2002,113(3):363-378.
[15]CHEN Z,MENG X,GUO L,et al.GICUDA:A parallel program for 3Dcorrelation imaging of large scale gravity and gravity gradiometry data on graphics processing units with CUDA[J].Computers&Geosciences,2012,46(3):119-128.
[16]孟小红,刘国峰,陈召曦,等.基于剩余异常相关成像的重磁物性反演方法[J].地球物理学报,2012,55(1):304-309.
[17]郭良辉,孟小红,石磊,等.重力和重力梯度数据三维[J].地球物理学报,2009,52(4),1098-1106.
[18]许令周.概率成像方法对台阶模型的成像[J].西南民族大学学报(自然科学版),2007,33(2):364-367.
[19]许令周.重力场中概率成像方法对脉状矿体的识别[J].内蒙古师范大学学报(自然科学汉文版),2007(4):463-465,469.
[20]马国庆,杜晓娟,李丽丽.改进的位场相关成像方法[J].地球科学(中国地质大学学报),2013,38(5):1121-1127.
[21]PATELLA D.Introduction to ground surface self-potential tomography[J].Geophysical Prospecting,1997a,45,653-681.
[22]PATELLA D.Self-potential global tomography including topographic effects[J].Geophysical Prospecting,1997b,45,843-863.
[23]MAURIELLO P,MONNA D,PATELLA D.3D geoelectric tomography and archaeological applications[J].Geophysical Prospecting,1998,46:543-570.
[24]MAURIELLO P,PATELLA D.Resistivity anomaly imaging by probability tomography[J].Geophysical Prospecting,1999a,47:411-429.
[25]MAURIELLO P,PATELLA D.Principles of probability tomography for natural-source electromagnetic induction fields,Geophysics,1999b,64,1403-1417.
[26]阮帅,王绪本,高永才,等.位场正则化下延在自然电位概率成像中的应用[J].物探化探计算技术,2006(4):332-336,285.
[27]王绪本,高永才.基于正则化向下延拓的重磁概率成像方法研究[C]//2008年重磁数据处理解释应用研讨会论文集.成都:成都理工大学,2008:49-54.
[28]黄金明.三维局部重磁场源全方位成像可视化研究[D].北京:中国地质大学(北京),2003.
[29]MAURIELLO P,PATELLA D.Gravity probability tomography:a new tool for buried mass distribution imaging[J].Geophysical Prospecting,2001,49:1-12.
[30]吉日嘎拉图,鲁光银.基于线性回归的重力数据三维RMS成像[J].物探化探计算技术,2014,36(2):185-189.
[31]杨旭明,王克斌.一种层析成像的正演和反演方法[J].物探化探计算技术,2002,24(4):337-340.
[32]闫浩飞,刘国峰.一种重力异常概率成像的扩展计算[J].地球物理学进展,2014,29(4):1837-1842.
[33]徐龙强,刘志义.物探化探数据成像预处理[J].物探化探计算技术,1991,13(1):41-46.
[34]路利春,文博,程斌,等.成像技术在矿区重力勘探中的应用研究[J].地质与资源,2018(1):99-102,88.
[35]许令周,关继腾,房文静.减小概率成像分辨宽度的一种方法[J].物探化探计算技术,2005,27(2):138-140,93.