四川地区流动重力资料的位场 分离与异常特征提取

李大虎, 丁志峰, 梁明剑, 李军, 苏琴

李大虎, 丁志峰, 梁明剑, 李军, 苏琴. 2014: 四川地区流动重力资料的位场 分离与异常特征提取. 地震学报, 36(2): 261-274. DOI: 10.3969/j.issn.0253-3782.2014.02.011
引用本文: 李大虎, 丁志峰, 梁明剑, 李军, 苏琴. 2014: 四川地区流动重力资料的位场 分离与异常特征提取. 地震学报, 36(2): 261-274. DOI: 10.3969/j.issn.0253-3782.2014.02.011
Li Dahu, Ding Zhifeng, Liang Mingjian, Li Jun, Su Qin. 2014: Field separation and anomaly feature extraction of mobile gravity data in Sichuan area. Acta Seismologica Sinica, 36(2): 261-274. DOI: 10.3969/j.issn.0253-3782.2014.02.011
Citation: Li Dahu, Ding Zhifeng, Liang Mingjian, Li Jun, Su Qin. 2014: Field separation and anomaly feature extraction of mobile gravity data in Sichuan area. Acta Seismologica Sinica, 36(2): 261-274. DOI: 10.3969/j.issn.0253-3782.2014.02.011

四川地区流动重力资料的位场 分离与异常特征提取

基金项目: 中国地震局地震科技星火计划项目(XH12041Y)和四川省地震局地震科技专项(LY1312)联合资助.
详细信息
    通讯作者:

    李大虎, e-mail: lixiang2006@sina.com

  • 中图分类号: P312.1

Field separation and anomaly feature extraction of mobile gravity data in Sichuan area

  • 摘要: 用于监测四川地区重力场随时间变化的流动重力测网于2008年5月12日汶川MS8.0地震后进行了全面改造, 形成了新的闭合环线并对整网点位每年实施2次重力观测工作. 该文以2010—2012年获取的5期流动重力资料为基础, 分别采用小波多尺度分解和三维视密度反演等位场分离方法, 对测区内重力场的动态变化进行多层次分解及异常特征提取, 以探求四川测区范围内重力异常动态变化特征以及地震前后不同深度范围内密度变化情况与重力场响应之间的关系. 结果表明, 对同一期资料进行小波分解得到不同阶次的细节图像, 揭示了不同深度处重力异常动态变化特征; 分解结果同时也表明了四川测区范围内重力场的变化形态与川内特有的断裂构造格架和深部结构环境及过渡带特征密切相关. 通过对不同时期的差分重力异常资料进行三维视密度反演对比研究发现, 在地震前后四川测区范围内由浅到深密度变化呈现出不同程度的差异, 震区附近出现较为明显的横向密度变化特征且浅层物质密度局部性的异常变化差异较为突出, 测区范围内中上地壳深度层密度反演切片显示出的趋势性差异基本上反映了地震前后测区内地下物质的动态变化趋势, 因此, 深刻认识不同深度范围内重力场的动态演化特征, 对研究壳内物质密度变化以及与地震的孕育和发展有关的深部介质环境具有积极的意义.
    Abstract: The mobile gravity measuring network for monitoring the Sichuan area gravitational field varying along with the time has made an overall transformation after the Wenchuan MS8.0 earthquakes on 12 May 2008, forming a new closed loop with gravity observation to all of the point locations twice every year. Taking five-issue mobile gravity data in the period of 2010—2012 as the foundation, this paper detects dynamic change of gravitational field in multi-scale decomposition and extracts anomaly features by using wavelet multi-scale decomposition and the 3-D apparent density inversion, so as to search for dynamic change characteristics of gravity anomaly and the relationship between the range of density variation and the gravity field response within different depths before and after the earthquake in Sichuan measuring area. The results show that different order details chart got from wavelet decomposition on the same stage gravity data reveals dynamic change characteristics of gravity anomaly in the different depth. Decomposition results also show that changes of gravitational field are closely related with fracture structure framework, deep structure and characteristics of transitional zone in Sichuan measuring area. Comparison of 3-D apparent density inversions based on the differential gravity anomaly data of different periods suggests density change is different from shallow to deep depth before and after an earthquake. Variation of transverse density characteristics is obviously presented nearby the seismic region, and anomalous change difference in the density of shallow matter is especially prominent. The apparent density slice maps within the depth of the middle-upper crust of the testing area shows trend differences, which basically reflects the dynamic change process of underground substances before and after the earthquake. Therefore deep understanding of gravity dynamic evolution characteristics in the different depth bounds is of positive significance to study the matter density change in the crust and deep medium environment related to the pregnancy and development of the earthquake.
  • 图  1   四川地区流动重力测网及构造略图

    Figure  1.   Sketch of mobile gravity measuring network and tectonic setting in Sichuan region

    图  2   1109期小波多尺度分解结果(单位: 10-8 m/s2)

    (a) 重力场变化等值线图; (b) 一阶细节; (c) 二阶细节; (d) 三阶细节; (e) 四阶细节; (f) 四阶逼近F1: 龙门山断裂;F2: 抚边河断裂;F3: 鲜水河断裂; F4: 玉农希断裂;F5: 理塘断裂; F6: 安宁河断裂;F7: 则木河断裂;F8: 小江断裂; F9: 大凉山断裂

    Figure  2.   Wavelet multi-scale decomposition results of 1 109-issue gravity measurements (unit: 10-8 m/s2)

    (a) Contour map of gravitational field variation; (b) The first-order details; (c) The second-order details;(d) The third-order details; (e) The fourth-order details; (f) The fourth-order approximation F1: Longmenshan fault; F2: Fubianhe fault; F3: Xianshuihe fault; F4: Yunongxi fault; F5: Litang fault; F6: Anninghe fault; F7: Zemuhe fault; F8: Xiaojiang fault; F9: Daliangshan fault

    图  3   2010—2012年重力场变化等值线图(单位: 10-8 m/s2)

    (a) 2010-08—2011-03; (b) 2011-03—09; (c) 2011-09—2012-04; (d) 2012-04—09

    Figure  3.   Contour map of gravitational field variation (in unit of 10-8 m/s-2) during 2010—2012

    (a) 2010-08—2011-03; (b) 2011-03—09; (c) 2011-09—2012-04; (d) 2012-04—09

    图  4   2010—2012年观测区不同深度层视密度反演结果(单位: kg/m3)

    (a)—(c) 2010-08—2011-03; (d)—(f) 2011-03—09; (g)—(i) 2011-09—2012-04; (j)—(l) 2012-04—09

    Figure  4.   Apparent density inversion results of the layers at the depths of 5,10 and 20 km during the 2010—2012

    (a)—(c) 2010-08—2011-03; (d)—(f) 2011-03—09; (g)—(i) 2011-09—2012-04; (j)—(l) 2012-04—09

    表  1   2010—2012年四川地区重力测量情况

    Table  1   Gravity measur ements of Sichuan region during 2010—2012

    下载: 导出CSV
  • 陈石 , 王谦身, 祝意青, 蒋长胜, 卢红艳, 王武星. 2011. 汶川MS8.0震前区域重力场变化与震质中研究[J]. 地球物理学进展, 26 (4): 1147-1156.

    Chen S, Wang Q S, Zhu Y Q, Jiang C S, Lu H Y, Wang W X. 2011. Regional gravity variation before Wenchuan MS8.0 earthquake and epicentroid research[J]. Progress in Geophysics, 26 (4): 1147-1156 (in Chinese).

    陈运泰, 顾浩鼎, 卢造勋. 1980. 1975年海城地震和1976年唐山地震前后重力变化[J]. 地震学报, 2 (1): 21-3.

    Chen Y T, Gu H D, Lu Z X. 1980. Variations of gravity before and after the Haicheng earthquake, 1975 and the Tang-shan earthquake, 1976[J]. Acta Seismologica Sinica, 2 (1): 21-31 (in Chinese).

    刁博, 王家林, 程顺有. 2007. 重力异常小波多分辨分析分解阶次的确定[J]. 地球科学: 中国地质大学学报, 32 (4): 564-568.

    Diao B, Wang J L, Cheng S Y. 2007. The confirmation of decomposition level in wavelet multi-resolution analysis for gravity anomalies[J]. Earth Science: Journal of China University of Geosciences, 32 (4): 564-568 (in Chinese).

    顾功叙, 郭宗汾, 刘克人, 郑金涵, 卢红艳, 刘端法. 1997. 中国京津唐张地区时间上连续的重力变化与地震的孕育和发生[J]. 科学通报, 42 (18): 1919-1930.

    Gu G X, Guo Z F, Liu K R, Zheng J H, Lu H Y, Lu DF. 1997. Time continuity of gravity variation in the BJTZ region and earthquake development and occurrence[J]. Chinese Science Bulletin, 42 (18): 1919-1930 (in Chinese).

    李大虎, 何强, 李才明, 屈进红. 2011. MATLAB与C++混合编程实现航磁异常提取的小波分析方法研究[J]. 地震研究, 34 (2): 121-126.

    Li D H, He Q, Li C M, Qu J H. 2011. Research on wavelet analysis method of magnetic anomaly extraction based on MATLAB and C++ mixed programming[J]. Journal of Seismological Research, 34 (2): 121-126 (in Chinese).

    李辉, 申重阳, 孙少安, 王晓权, 项爱民, 刘少明. 2009. 中国大陆近期重力场动态变化图像[J]. 大地测量与地球动力学, 29 (3): 1-10.

    Li H, Shen C Y, Sun S A, Wang X Q, Xiang A M, Liu S M. 2009. Dynamic gravity change in recent years in China Continent[J]. Journal of Geodesy and Geodynamics, 29 (3): 1-10 (in Chinese).

    李振海, 汪海洪. 2010. 重力数据网格化方法比较[J]. 大地测量与地球动力学, 30 (1): 140-144.

    Li Z H, Wang H H. 2010. Comparison among methods for gravity data gridding[J]. Journal of Geodesy and Geodynamics, 30 (1): 140-144 (in Chinese).

    刘少明, 申重阳, 孙少安, 陈步云. 2005. 三峡库首区重复重力资料的小波分解与构造活动[J]. 大地测量与地球动力学, 25 (1): 13-18.

    Liu S M, Shen C Y, Sun S A, Chen B Y. 2005. Wavelet decomposition of repeated gravity data and tectonic activities in front area of Three Gorges reservoir[J]. Journal of Geodesy and Geodynamics, 25 (1): 13-18 (in Chinese).

    马宏生, 张国民, 周龙泉, 刘杰, 邵志刚, 夏红. 2008. 川滇地区中小震重新定位于速度结构的联合反演研究[J]. 地震, 28 (2): 29-38.

    Ma H S, Zhang G M, Zhou L Q, Liu J, Shao Z G, Xia H. 2008. Simultaneous inversion of small earthquake relocation and velocity structure in Sichuan-Yunnan area[J]. Earthquake, 28 (2): 29-38 (in Chinese).

    申重阳, 李辉, 付广裕. 2003. 丽江7.0级地震重力前兆模式研究[J]. 地震学报, 25 (2): 163-171.

    Shen C Y, Li H, Fu G Y. 2003. Study on a gravity precursor mode of Lijiang earthquake with MS7.0 [J]. Acta Seismologica Sinica, 25 (2): 163-171 (in Chinese).

    申重阳, 李辉, 孙少安. 2009. 重力场动态变化与汶川MS8.0地震孕育过程[J]. 地球物理学报, 52 (10): 2547-2557.

    Shen C Y, Li H, Su S A. 2009. Dynamic variations of gravity and the preparation process of the Wenchuan MS8.0 earthquake[J]. Chinese Journal of Geophysics, 52 (10): 2547-2557 (in Chinese).

    孙少安, 项爱民, 李辉. 2002. 长江三峡库首区重力场变化与构造活动[J]. 大地测量与地球动力学, 22 (3): 56-59.

    Sun S A, Xiang A M, Li H. 2002. Changes of gravity field and tectonic activity in front region of the Three Gorges project[J]. Journal of Geodesy and Geodynamics, 22 (3): 56-59 (in Chinese).

    唐荣昌, 韩渭宾. 1993. 四川活动断裂与地震[M]. 北京: 地震出版社: 33-38.

    Tang R C, Han W B.1993. Sichuan Active Fault and Earthquake[M]. Beijing: Seismologica Press: 33-38 (in Chinese).

    文百红, 程方道. 1990. 用于划分磁异常的新方法: 插值切割法[J]. 中南矿冶学院学报, 21 (3): 229-334.

    Wen B H, Chen F D.1990. A new interpolating cut method for identifying reglonal and local fields of magnetic anomaly [J]. Journal of Central-South Institute of Mining and Metallurgy, 21 (3): 229-334 (in Chinese).

    徐世浙, 余海龙, 李海侠, 田钢, 杨金玉. 2009. 基于位场分离与延拓的视密度反演[J]. 地球物理学报, 52 (6): 1592-1598.

    Xu S Z, Yu H L, Li H X, Tian G, Yang J Y. 2009. The inversion of apparent density based on the separation and continuation of potential field[J]. Chinese Journal of Geophysics, 52 (6): 1592-1598 (in Chinese).

    玄松柏, 申重阳, 李辉. 2010. 丽江7.0级地震前后地壳密度动态变化三维反演[J]. 大地测量与地球动力学, 30 (1): 14-20.

    Xuan S B, Shen C Y, Li H. 2010. On 3D inversion for crustal density change before and after Lijiang MS7.0 earthquake[J]. Journal of Geodesy and Geodynamics, 30 (1): 14-20 (in Chinese).

    杨文采, 施志群, 侯遵泽. 2001. 离散小波变换与重力异常多重分解[J]. 地球物理学报, 44 (4): 534-541.

    Yang W C, Shi Z Q, Hou Z Z. 2001. Discrete wavelet transform for multiple decompo sition of gravity anoma lies[J]. Chinese Journal of Geophysics, 44 (4): 534-541 (in Chinese).

    朱艾斓, 徐锡伟, 周永胜, 尹京苑, 甘卫军, 陈桂华. 2005. 川西地区小震重新定位及其活动构造意义[J]. 地球物理学报, 48 (3): 629-636

    Zhu A L, Xu X W, Zhou Y S, Yin J Y, Gan W J, Chen G H. 2005. Relocating of small earthquakes in western Sichuan, China and its implications for active tectonics[J]. Chinese Journal of Geophysics, 48 (3): 629-636 (in Chinese).

    朱桂芝, 祝意青, 唐文, 程宏宾. 2004. 河西地区重力动态变化二维离散多尺度分解[J]. 大地测量与地球动力学, 24 (1): 46-49.

    Zhu G Z, Zhu Y Q, Tang W, Cheng H B. 2004. Multiscale decomposition of 2D discrete wavelet of dynamic gravity variation in Hexi region[J]. Journal of Geodesy and Geodynamics, 24 (1): 46-49 (in Chinese).

    祝意青, 徐云马, 吕弋培, 李铁明. 2009. 龙门山断裂带重力变化与汶川8.0级地震关系研究[J]. 地球物理学报, 52 (10): 2538-2546.

    Zhu Y Q, Xu Y M, Lü Y P, Li T M. 2009. Relations between gravity variation of Longmenshan fault zone and Wenchuan M8.0 earthquake[J]. Chinese Journal of Geophysics, 52 (10): 2538-2546 (in Chinese).

    祝意青, 苏琴, 梁伟锋, 刘芳. 2010. 川西地区重力变化与汶川地震[J]. 大地测量与地球动力学, 30 (2): 1-5.

    Zhu Y Q, Su Q, Liang W F, Liu F. 2010. Gravity variation of western Sichuan and Wenchuan earthquake [J]. Journal of Geodesy and Geodynamics, 30 (2): 1-5 (in Chinese).

    Boschetti F, Hornby P, Horowitz F G. 2001. Wavelet based inversion of gravity data[J]. Exploration Geophysics, 32 (1): 48-55.

    Burchfiel B C. 2003. New technology: New geological challenges[J]. Geol Soc Amer Today, 14 (2): 4-10.

    Cooper G. 2004. The stable downward continuation of potential field data[J]. Exploration Geophysics, 35 (2): 260- 265.

    Fedi M, Florio G. 2002. A stable downward continuation by using the ISVD method[J]. Geophys J Int, 151 (1): 146-156.

    Møller M J, Olsen H, Ploug C, Strykowski G, Hjorth H. 2007. Gravity field separation and mapping of buried quaternary valleys in Lolland, Denmark using old geophysical data[J]. J Geodyn, 43 (2): 330-337.

图(4)  /  表(1)
计量
  • 文章访问数:  694
  • HTML全文浏览量:  297
  • PDF下载量:  21
  • 被引次数: 0
出版历程
  • 收稿日期:  2013-01-18
  • 修回日期:  2013-12-04
  • 发布日期:  2014-02-28

目录

    /

    返回文章
    返回