基于全张量磁梯度方向Theta法进行边界识别的研究与应用
|
摘要
边界识别是位场(重磁)数据解释的重要内容,其对圈定异常范围,划分断裂构造等具有重大意义.Theta法是众多边界识别方法中一种应用较为广泛的方法,其为总水平导数与解析信号定义的比值函数.Theta的极大值对应着地质体的边界,具有能均衡不同埋深地质体边界的优势,但其受斜磁化的影响较为严重.全张量磁梯度数据具有高精度、高分辨率、多参量的优点,且其受斜磁化的影响较小,成为近年来位场数据处理与解释的研究热点.为了减少Theta法受斜磁化的影响,本文从Theta法的基本公式出发,推导出了由全张量磁梯度数据定义的新的Theta法表达式,本文称之为全张量磁梯度方向Theta法.通过模型试验,验证了全张量磁梯度方向Theta法的理论可行性,其边界识别的效果不仅十分清晰,而且受高纬度斜磁化情况下的影响较小.将该方法应用到大兴安岭北段某地区实测航磁数据的解释中,获得了异常体边界的大致分布.
Edge detection is an important part of the interpretation of the potential field(gravity and magnetic) data,which is of great significance for delineating the anomaly range and dividing the fracture structure. The Theta method is a widely used method in many edge detection methods,which is a ratio function defined by the total horizontal derivative and the analytical signal. The maximum value of Theta corresponds to the boundary of the geological body,and has the advantage of balancing the boundaries of different buried depths,but it is more affected by the oblique magnetization. The full tensor magnetic gradient data has the advantages of high precision,high resolution and multi-parameters,and its influence by oblique magnetization is small,which has become a research hotspot of data processing and interpretation of potential field in recent years. In order to reduce the influence of the Theta method on the oblique magnetization,this paper derives the new Theta expression defined by the full tensor magnetic gradient data from the basic formula of Theta method. This paper is called the full tensor magnetic gradient direction Theta method. Through the model test,the theoretical feasibility of the full tensor magnetic gradient direction Theta method is verified. The effect of edge detection is not only clear,but also less affected by high latitude oblique magnetization. Applying the method to the interpretation of the measured aeromagnetic data in a certain area in the northern section of the Greater Khingan Mountains,the approximate distribution of the anomalous body boundaries is obtained.
Edge detection is an important part of the interpretation of the potential field(gravity and magnetic) data,which is of great significance for delineating the anomaly range and dividing the fracture structure. The Theta method is a widely used method in many edge detection methods,which is a ratio function defined by the total horizontal derivative and the analytical signal. The maximum value of Theta corresponds to the boundary of the geological body,and has the advantage of balancing the boundaries of different buried depths,but it is more affected by the oblique magnetization. The full tensor magnetic gradient data has the advantages of high precision,high resolution and multi-parameters,and its influence by oblique magnetization is small,which has become a research hotspot of data processing and interpretation of potential field in recent years. In order to reduce the influence of the Theta method on the oblique magnetization,this paper derives the new Theta expression defined by the full tensor magnetic gradient data from the basic formula of Theta method. This paper is called the full tensor magnetic gradient direction Theta method. Through the model test,the theoretical feasibility of the full tensor magnetic gradient direction Theta method is verified. The effect of edge detection is not only clear,but also less affected by high latitude oblique magnetization. Applying the method to the interpretation of the measured aeromagnetic data in a certain area in the northern section of the Greater Khingan Mountains,the approximate distribution of the anomalous body boundaries is obtained.
引文
Codell L,Grauch V J S.1985.Mapping basement magnetization zones from aeromagnetic data in the San Juan Basin,New Mexico [M].// Hinze W J ed.The Utility of Regional Gravity and Magnetic Anomaly Maps.Tulsa,OK:Society of Exploration Geophysics,181-197.
Cooper G R J,Cowan D R.2006.Enhancing potential field data using filters based on the local phase [J].Computers & Geosciences,32(10):1585-1591.
Cordell L.1979.Gravimetric expression of graben faulting in Santa Fe Country and the Espanola Basin,New Mexico [M].// New Mexico Geological Society Guidebook,30th Annual Fall Field Conference.New Mexico:New Mexico Geological.Society,59- 64.
Evjen H M.1936.The place of the vertical gradient in gravitational interpretations [J].Geophysics,1(1):127-136.
Guan Z N.2005.Geomagnetic Field and Magnetic Exploration [M].Beijing:Geological Publishing House.
Hood P,McClure D J.1965.Gradient measurements in ground magnetic prospecting [J].Geophysics,30(3):403- 410.
Hood P J,Teskey D J.1989.Aeromagnetic gradiometer program of the Geological Survey of Canada [J].Geophysics,54(8):1012-1022.
Hsu S K,Sibuet J C,Shyu C T.1996.High-resolution detection of geologic boundaries from potential field anomalies:an enhanced analytic signal technique [J].Geophysics,61(2):373-386.
Ma G Q,Du X J,Li L L.2011.Edge detection of potential field data using correlation coefficients of horizontal and vertical derivatives [J].Journal of Jilin University (Earth Science Edition) (in Chinese),41(S1):345-348,doi:10.13278/j.cnki.jjuese.2011.s1.010.
Ma G Q,Huang D N,Yu P,et al.2012.Application of improved balancing filters to edge identification of potential field data [J].Chinese Journal of Geophysics (in Chinese),55(12):4288- 4295,doi:10.6038/j.issn.0001-5733.2012.12.040.
Ma G Q,Li L L.2012.Edge detection in potential fields with the normalized total horizontal derivative [J].Computer & Geosciences,41:83-87.
Ma G Q,Li L L,Du X J.2012.Boundary detection and interpretation by magnetic gradient tensor data [J].Oil Geophysical Prospecting (in Chinese),47(5):815-821,doi:10.13810/j.cnki.issn.1000-7210.2012.05.020.
Miller H G,Singh V.1994.Potential field tilt—a new concept for location of potential field sources [J].Journal of Applied Geophysics,32(2-3):213-217.
Nabighian M N.1972.The analytic signal of two-dimensional magnetic bodies with polygonal cross-section:Its properties and use for automated anomaly interpretation [J].Geophysics,37(3):507-517.
Nabighian M N.1984.Toward a three-dimensional automatic interpretation of potential field data via generalized Hilbert transforms:Fundamental relations [J].Geophysics,49(6):780-786.
Roest W R,Verhoef J,Pilkington M.1992.Magnetic interpretation using the 3-D analytic signal [J].Geophysics,57(1):116-125.
She H Q,Li J W,Xiang A P,et al.2012.U-Pb ages of the Zircons from primary rocks in middle-northern Daxinganling and its implications to geotectonic evolution [J].Acta petrologica Sinica (in Chinese),28(2):571-594.
Shu Q,Ma G Q,Liu C,et al.2018.Balanced edge detection and depth imaging technique for the interpretation of full-tensor magnetic gradient data [J].Chinese J.Geophys.(in Chinese),61(4):1539-1548,doi:10.6038/cjg2018L0726.
Sun A,Guo H,Tian M Y,et al.2017.Application of aeromagnetic three-component survey based on the Euler deconvolution method [J].Chinese Journal of Geophysics (in Chinese),60(11):4491- 4505,doi:10.6038/cjg20171133.
Verduzco B,Fairhead J D,Green C M,et al.2004.New insights into magnetic derivatives for structural mapping [J].The Leading Edge,23(2):116-119.
Wang L F,Guo C C,Xue D J,et al.2016.Analytic signals of magnetic gradient tensor and their application to estimate source location [J].Progress in Geophysics (in Chinese),31(3):1164-1172,doi:10.6038/pg20160333.
Wijns C,Perez C,Kowalczyk P.2005.Theta map:Edge detection in magnetic data [J].Geophysics,70(4):L39-L43.
管志宁.2005.地磁场与磁力勘探[M].北京:地质出版社.
马国庆,杜晓娟,李丽丽.2011.利用水平与垂直导数的相关系数进行位场数据的边界识别[J].吉林大学学报(地球科学版),41(S1):345-348,doi:10.13278/j.cnki.jjuese.2011.s1.010.
马国庆,黄大年,于平,等.2012.改进的均衡滤波器在位场数据边界识别中的应用[J].地球物理学报,55(12):4288- 4295,doi:10.6038/j.issn.0001-5733.2012.12.040.
马国庆,李丽丽,杜晓娟.2012.磁张量数据的边界识别和解释方法[J].石油地球物理勘探,47(5):815-821,doi:10.13810/j.cnki.issn.1000-7210.2012.05.020.
佘宏全,李进文,向安平,等.2012.大兴安岭中北段原岩锆石U-Pb测年及其与区域构造演化关系[J].岩石学报,28(2):571-594.
舒晴,马国庆,刘财,等.2018.全张量磁梯度数据解释的均衡边界识别及深度成像技术[J].地球物理学报,61(4):1539-1548,doi:10.6038/cjg2018L0726.
孙昂,郭华,田明阳,等.2017.基于欧拉反褶积方法的航磁三分量应用研究[J].地球物理学报,60(11):4491- 4505,doi:10.6038/cjg20171133.
王林飞,郭灿灿,薛典军,等.2016.磁梯度张量解析信号分析法及其在场源位置识别中的应用[J].地球物理学进展,31(3):1164-1172,doi:10.6038/pg20160333.
Cooper G R J,Cowan D R.2006.Enhancing potential field data using filters based on the local phase [J].Computers & Geosciences,32(10):1585-1591.
Cordell L.1979.Gravimetric expression of graben faulting in Santa Fe Country and the Espanola Basin,New Mexico [M].// New Mexico Geological Society Guidebook,30th Annual Fall Field Conference.New Mexico:New Mexico Geological.Society,59- 64.
Evjen H M.1936.The place of the vertical gradient in gravitational interpretations [J].Geophysics,1(1):127-136.
Guan Z N.2005.Geomagnetic Field and Magnetic Exploration [M].Beijing:Geological Publishing House.
Hood P,McClure D J.1965.Gradient measurements in ground magnetic prospecting [J].Geophysics,30(3):403- 410.
Hood P J,Teskey D J.1989.Aeromagnetic gradiometer program of the Geological Survey of Canada [J].Geophysics,54(8):1012-1022.
Hsu S K,Sibuet J C,Shyu C T.1996.High-resolution detection of geologic boundaries from potential field anomalies:an enhanced analytic signal technique [J].Geophysics,61(2):373-386.
Ma G Q,Du X J,Li L L.2011.Edge detection of potential field data using correlation coefficients of horizontal and vertical derivatives [J].Journal of Jilin University (Earth Science Edition) (in Chinese),41(S1):345-348,doi:10.13278/j.cnki.jjuese.2011.s1.010.
Ma G Q,Huang D N,Yu P,et al.2012.Application of improved balancing filters to edge identification of potential field data [J].Chinese Journal of Geophysics (in Chinese),55(12):4288- 4295,doi:10.6038/j.issn.0001-5733.2012.12.040.
Ma G Q,Li L L.2012.Edge detection in potential fields with the normalized total horizontal derivative [J].Computer & Geosciences,41:83-87.
Ma G Q,Li L L,Du X J.2012.Boundary detection and interpretation by magnetic gradient tensor data [J].Oil Geophysical Prospecting (in Chinese),47(5):815-821,doi:10.13810/j.cnki.issn.1000-7210.2012.05.020.
Miller H G,Singh V.1994.Potential field tilt—a new concept for location of potential field sources [J].Journal of Applied Geophysics,32(2-3):213-217.
Nabighian M N.1972.The analytic signal of two-dimensional magnetic bodies with polygonal cross-section:Its properties and use for automated anomaly interpretation [J].Geophysics,37(3):507-517.
Nabighian M N.1984.Toward a three-dimensional automatic interpretation of potential field data via generalized Hilbert transforms:Fundamental relations [J].Geophysics,49(6):780-786.
Roest W R,Verhoef J,Pilkington M.1992.Magnetic interpretation using the 3-D analytic signal [J].Geophysics,57(1):116-125.
She H Q,Li J W,Xiang A P,et al.2012.U-Pb ages of the Zircons from primary rocks in middle-northern Daxinganling and its implications to geotectonic evolution [J].Acta petrologica Sinica (in Chinese),28(2):571-594.
Shu Q,Ma G Q,Liu C,et al.2018.Balanced edge detection and depth imaging technique for the interpretation of full-tensor magnetic gradient data [J].Chinese J.Geophys.(in Chinese),61(4):1539-1548,doi:10.6038/cjg2018L0726.
Sun A,Guo H,Tian M Y,et al.2017.Application of aeromagnetic three-component survey based on the Euler deconvolution method [J].Chinese Journal of Geophysics (in Chinese),60(11):4491- 4505,doi:10.6038/cjg20171133.
Verduzco B,Fairhead J D,Green C M,et al.2004.New insights into magnetic derivatives for structural mapping [J].The Leading Edge,23(2):116-119.
Wang L F,Guo C C,Xue D J,et al.2016.Analytic signals of magnetic gradient tensor and their application to estimate source location [J].Progress in Geophysics (in Chinese),31(3):1164-1172,doi:10.6038/pg20160333.
Wijns C,Perez C,Kowalczyk P.2005.Theta map:Edge detection in magnetic data [J].Geophysics,70(4):L39-L43.
管志宁.2005.地磁场与磁力勘探[M].北京:地质出版社.
马国庆,杜晓娟,李丽丽.2011.利用水平与垂直导数的相关系数进行位场数据的边界识别[J].吉林大学学报(地球科学版),41(S1):345-348,doi:10.13278/j.cnki.jjuese.2011.s1.010.
马国庆,黄大年,于平,等.2012.改进的均衡滤波器在位场数据边界识别中的应用[J].地球物理学报,55(12):4288- 4295,doi:10.6038/j.issn.0001-5733.2012.12.040.
马国庆,李丽丽,杜晓娟.2012.磁张量数据的边界识别和解释方法[J].石油地球物理勘探,47(5):815-821,doi:10.13810/j.cnki.issn.1000-7210.2012.05.020.
佘宏全,李进文,向安平,等.2012.大兴安岭中北段原岩锆石U-Pb测年及其与区域构造演化关系[J].岩石学报,28(2):571-594.
舒晴,马国庆,刘财,等.2018.全张量磁梯度数据解释的均衡边界识别及深度成像技术[J].地球物理学报,61(4):1539-1548,doi:10.6038/cjg2018L0726.
孙昂,郭华,田明阳,等.2017.基于欧拉反褶积方法的航磁三分量应用研究[J].地球物理学报,60(11):4491- 4505,doi:10.6038/cjg20171133.
王林飞,郭灿灿,薛典军,等.2016.磁梯度张量解析信号分析法及其在场源位置识别中的应用[J].地球物理学进展,31(3):1164-1172,doi:10.6038/pg20160333.
© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号 地址：北京市海淀区学院路29号 邮编：100083 电话：办公室：(+86 10)66554848；文献借阅、咨询服务、科技查新：66554700 |