文章快速检索     高级检索
  大地测量与地球动力学  2019, Vol. 39 Issue (3): 307-312  DOI: 10.14075/j.jgg.2019.03.017

引用本文  

毛经伦, 李辉, 祝意青, 等. 利用EGM2008地球重力场模型研究川西地区水平构造应力场[J]. 大地测量与地球动力学, 2019, 39(3): 307-312.
MAO Jinglun, LI Hui, ZHU Yiqing, et al. Research on Horizontal Tectonic Stress Field Based on EGM2008 in Western Sichuan[J]. Journal of Geodesy and Geodynamics, 2019, 39(3): 307-312.

项目来源

国家自然科学基金(41374026)。

Foundation support

National Natural Science Foundation of China, No. 41374026.

通讯作者

李辉,研究员,主要研究方向为地球重力场及其变化观测, E-mail:lihuieq@163.com

Corresponding author

LI Hui, researcher, majors in observation and research of the gravity field and its changing, E-mail: lihuieq@163.com.

第一作者简介

毛经伦,硕士生,主要从事地震重力测量、时变重力及其与强震关系的研究,E-mail:371857582@qq.com

About the first author

MAO Jinglun, postgraduate, majors in seismic gravimetry, time-varying gravity and its relation to strong earthquakes, E-mail:371857582@qq.com.

文章历史

收稿日期:2018-04-17
利用EGM2008地球重力场模型研究川西地区水平构造应力场
毛经伦1,2     李辉1     祝意青2     李忠亚1     
1. 中国地震局地震研究所地震大地测量重点实验室,武汉市洪山侧路40号,430071;
2. 中国地震局第二监测中心,西安市西影路316号,710054
摘要:利用EGM2008地球重力场模型数据计算川西地区的水平构造应力场。结果显示,采用重力场数据可以快速获得水平构造应力场图像;龙门山断裂带及鲜水河-安宁河断裂带水平构造应力聚集程度较高,构造应力值较大,反映出这些地区地质构造运动活跃;水平构造应力方向从巴颜喀拉块体到华南地块再到川滇块体呈现顺时针方向旋转;构造应力场随着深度的增加呈逐渐减小的趋势,地震活动频率也相应降低;构造应力场等值线密集区构造运动活跃、地震活动频发,表明地震活动与构造应力聚集分布存在一定的对应关系。
关键词EGM2008重力场构造应力场川西地区

地球重力场研究是间接研究构造应力场的有效方法[1-6]。国内学者研究建立重力异常及垂线偏差转换水平构造应力场的方法[2],计算出不同区域的水平构造应力场[2-4];提出动态重力场转换动态水平构造应力场的理论方法[5];运用有限元求解华北地区重力势能差产生的水平构造应力[6]。研究表明,通过地应力、静岩应力、构造应力的内在联系可以实现地球重力场数据向应力场数据的转换[3]。包含龙门山断裂带、鲜水河-安宁河断裂带等区域的川西地区,是研究青藏高原向东南扩展变形及其与华南地块相互作用的关键地带,也是探讨汶川MS8.0地震和芦山MS7.0地震构造变形机制、开展地震风险评估的关键地域[7-8]。本文选取构造活动强烈的川西地区重力观测数据,研究探讨水平构造应力场的转换。通过对川西地区(26.5°~32°N,100°~105°E)的静态构造应力场进行分析,探究构造应力场背景及构造运动与地震孕育之间的联系。

1 重力场与构造应力场的转换

应用希尔伯特三维位场转换公式[2],由重力异常gz解出重力水平分量gxgy

$ \left\{ \begin{array}{l} {g_x} = \\ - \frac{1}{{2\pi }}\int_{ - \infty }^{ + \infty } {\int_{ - \infty }^{ + \infty } {\frac{{(x - \xi ){g_z}}}{{{{[{{(x - \xi )}^2} + {{(y - \eta )}^2} + {H^2}]}^{3/2}}}}{\rm{d}}\xi {\rm{d}}\eta } } \\ {g_y} = \\ - \frac{1}{{2\pi }}\int_{ - \infty }^{ + \infty } {\int_{ - \infty }^{ + \infty } {\frac{{(y - \eta ){g_z}}}{{{{[{{(x - \xi )}^2} + {{(y - \eta )}^2} + {H^2}]}^{3/2}}}}{\rm{d}}\xi {\rm{d}}\eta } } \end{array} \right. $ (1)

式中,(xy)为计算点坐标,(ζη)为整个区域的流动坐标,H为空间延拓高度。

由式(1)得出重力总水平梯度的大小与方向:

$ \alpha = \arctan \frac{{{g_y}}}{{{g_x}}} $ (2)
$ {g_L} = {({g_x}^2 + {g_y}^2)^{1/2}} $ (3)

其中,重力总水平梯度的方向近似于水平最大主应力方向。

水平构造应力计算公式为[2]

$ \Delta {\sigma _{{\rm{xx}}}} = \frac{{g{\rho _{x, y}}{g_L}}}{{4\pi f{\rho _m}}} $ (4)

式中,ρm为地幔密度,取值3 300 kg/m3f为万有引力常数,取值6.67×10-11N·m2/kg2g为重力加速度,取值9.8 m/s2ρx, y为经过均衡改正的单位柱体密度:

$ {\rho _{x, y}} = {\rho _c} - {\sigma _{x, y}} + {(\sigma )_{x, y}} $ (5)

δx, y为地形补偿密度:

$ {\sigma _{x, y}} = \frac{{{\rho _c}{H_{x, y}}}}{D} $ (6)

(δ)x, y为山根补偿密度:

$ {(\sigma )_{x, y}} = \frac{{({\rho _m} - {\rho _c})({d_{x, y}} - D)}}{D} $ (7)

式中,ρc为地壳平均密度,取值2 670 kg/m3Hx, y为单位柱体地形的平均海拔高;dx, y为单位柱体的山根平均深度;D为补偿面深度,Airy模型取30 km。

2 川西地区构造应力场计算 2.1 数据来源

本文利用EGM2008地球重力场模型[9]数据计算静态构造应力场。该模型由地面区域平均重力数据、卫星测高数据、GRACE卫星重力数据等重力资料与全球高分辨率地形数据构成,可提供2 159次、2 190阶的球谐系数,空间分辨率为5′,即9 km分辨率。本文所使用的EGM2008地球重力场模型来源于斯克里普斯海洋研究所(SIO)2016年发布的V24.1重力模型和2014年发布的V18.1地形模型,数据空间分辨率为1′。其中,SIO V24.1重力模型由EGM2008的陆地重力异常数据及卫星测高的海域重力异常数据构成,新版EGM2008地球重力场模型对大陆重力场改动较小,主要对海洋重力场进行了优化[10],能较好地反映我国大陆因地形起伏引起的重力场变化[11];SIO V18.1地形模型由ETOPO1陆地地形数据与海洋测深及卫星测高所得到的海底地形数据构成[10],利用SIO V18.1地形模型绘制的川西地区地形图如图 1所示。

图 1 SIO V18.1川西地区地形等值线 Fig. 1 Topography contour map in western Sichuan depends on SIO V18.1
2.2 川西地区地壳浅层水平构造应力场

静态构造应力场计算时,积分区域经离散化、格网化后以2°×2°范围为积分区域进行计算。数据采用Kriging法插值格网。在重力场与构造应力场转换时,将选取范围外延2°,以消除边界效应。通过式(3)计算川西地区重力总水平梯度,利用式(4)计算川西地区地壳浅层水平构造应力场,绘制出川西地区地壳浅层水平构造应力矢量图(图 2)和地壳浅层水平构造应力等值线图(图 3)。图 2中箭矢长度表示水平构造应力大小,箭矢所指方向为箭矢起点处最大水平主应力的迹线方向。箭矢颜色越趋于红色表示水平构造应力越大,越趋于蓝色表示水平构造应力越小。红色折线表示活动断层,黑色折线表示块体分界线[12]图 3中越趋于红色表示水平构造应力场等值线越密集,构造应力值越大,越趋于蓝色表示水平构造应力场等值线越稀疏,构造应力值越小。红色五角星标示汶川MS8.0地震和芦山MS7.0地震震中位置。川西地区水平构造应力在0~17 MPa范围内。

图 2 川西地区地壳浅层水平构造应力矢量 Fig. 2 Ground horizontal tectonic stress field vector map in western Sichuan

图 3 川西地区地壳浅层水平构造应力等值线 Fig. 3 Ground horizontal tectonic stress field contour map in western Sichuan
2.3 川西地区构造应力场随深部变化结果

根据重力延拓法计算川西地区10 km、20 km和50 km深度处的重力总水平梯度,再由式(4)计算出不同深度处的水平构造应力场如图 4图 5图 6所示。蓝色五角星所在位置为相应深度水平构造应力最大值处,黑色菱形为汶川MS8.0地震和芦山MS7.0地震震中位置。

图 4 川西地区10 km深处水平构造应力等值线 Fig. 4 Horizontal tectonic stress field contour map of 10 km deep in western Sichuan

图 5 川西地区20 km深处水平构造应力等值线 Fig. 5 Horizontal tectonic stress field contour map of 20 km deep in western Sichuan

图 6 川西地区50 km深处水平构造应力等值线 Fig. 6 Horizontal tectonic stress field contour map of 50 km deep in western Sichuan
3 结果分析与讨论

图 2可见,川西地区构造应力比较复杂,北部龙门山地区的构造应力总体呈东南方向,东南部的构造应力总体呈近东西向。在龙门山断裂带区域(30°~32°N,102.5°~105°E),构造应力较大且聚集明显,表现为巴颜喀拉块体向东南方向运动时受到华南地块的强烈阻碍。在鲜水河-安宁河断裂带区域(29°~30°N,101.5°~102.5°E),构造应力主要表现为近南北向,表明该地区受到北部块体的挤压。

此外,龙门山应力区目前处于近东西向的挤压构造应力环境中。在横向上,水平构造应力由巴颜喀拉块体的北东-南西向到龙门山块体的近东西向,再到华南地块的北西西-南东东向,表现为顺时针方向旋转;在纵向上,水平构造应力由华南地块的北西西-南东东向,到川滇块体的近东-西向,亦表现为顺时针方向旋转。应力方向在不同块体间的转变,反映出青藏高原东缘巴颜喀拉块体在南东向水平运动时对龙门山块体产生了挤压,同时由于华南地块特别是四川盆地的强烈阻挡,导致龙门山推覆构造带出现明显的地壳缩短和山体隆升。该结果与通过震源机制解获得的最大水平主应力优势方向较为一致[13]

图 3可见,川西地区水平构造应力变化范围为0.050~16.958 MPa,这一结果与郭飞霄等[4]利用重力数据反演川滇地区构造应力场的结果相一致。应力集中或构造运动活跃的区域,等值线密集程度越高。龙门山断裂带(30°~32°N,102.5°~105°E)、鲜水河-安宁河断裂带(29°~30°N,101.5°~102.5°E)等构造强活动区,对应于等值线图中红色区域;等值线稀疏或无等值线出现的地带为非构造应力聚集区,表明相应地带的构造运动不活跃,反映了华南地块尤其是四川盆地(28.5°~31.5°N,103°~105°E)为构造运动较稳定区域;等值线图中马尔康、道孚、峨边、雅江等地区附近的等值线密集程度介于二者之间,为构造运动相对活跃区。与孟文等[13]获得的龙门山断裂带构造应力实测数据1.74~25.69 MPa相比,本文所得结果略小,可能是由于本文的构造应力仅来源于重力场数据反演,虽然体现了实际构造应力的大部分信息,但仍有部分构造应力信息缺失。但是,构造应力分布与实测数据较为一致,构造应力最大值均出现在康定附近,龙门山南段在较大尺度范围内地壳浅层构造应力强度较高,龙门山中段构造应力强度有所降低。

图 2图 3可以看出,川西地区恰好处于地壳活跃地带,历史上也曾发生过多次强烈地震,如汶川MS8.0地震和芦山MS7.0地震。由图可见,龙门山断裂带构造应力明显聚集,且构造应力数值高于其他地区,地壳浅层构造应力最大值接近17 MPa,说明该地区构造活动强烈。龙门山应力聚集区目前正处于近东西向的挤压构造应力环境之中,与该区域逆冲断层为主的强烈构造活动情况相一致。GPS速度场资料显示,汶川地震以前,龙门山块体运动速度由藏东块体到龙门山块体衰减了49%[14]。汶川地震以后,巴颜喀拉地块与川滇块体相对于华南地块的运动速率进一步降低[15],而该地区的水平构造应力相对于周边地区呈急剧增强态势,说明此区域断层在漫长时期(至少为万年尺度以上)是强活跃的,目前出现强闭锁现象。鲜水河-安宁河断裂带与龙门山断裂带交汇处,尤其在康定南部有较大的应力聚集,构造应力方向主要沿断裂带方向分布,构造活动强烈,且主要为走滑断层。川西地区南部有一定程度的构造应力聚集,构造活动较明显。华南地块尤其是四川盆地,应力分布较平缓,应力聚集程度较低,反映出该地区构造运动相对稳定。

图 4图 5图 6所示,各深度处构造应力最大值分别为:地壳浅层16.958 MPa、10 km深度11.388 MPa、20 km深度8.108 MPa、50 km深度4.484 MPa,构造应力大小随着深度增加呈递减趋势, 该结果与通过重力场算得的构造应力分布结果一致[3-4]。此外,10 km、20 km和50 km深度的水平构造应力最大值所在位置均处于汶川MS8.0地震和芦山MS7.0地震震中附近,说明该地区不同深度构造活动均比较活跃,汶川MS8.0地震和芦山MS7.0地震在一定程度上均印证了这一特征。

图 7为2009-01~2017-09不同深度处地震统计图,信息来源于中国地震台网中心地震数据管理与服务系统川西地区4级以上地震的统计资料。综合图 4~图 7可得,随着深度的增加构造应力逐渐变小,构造运动相对稳定,地震活动频率降低。

图 7 2009-01~2017-09不同震源深度的地震统计 Fig. 7 Seismic statistics of different source depths from January 2009 to September 2017

图 8给出了川西地区地壳浅层构造应力场与地震分布情况,信息来源于中国地震信息网1900~2017年川西地区6级以上地震的统计资料。图中红色圆点标示地震分布,圆点越大表示震级越大。从图 8可以看出,在构造应力场等值线相对密集的地区构造运动活跃、地震活动频发,表明地震活动与构造应力聚集分布存在一定的对应关系。龙门山断裂带与鲜水河-安宁河断裂带属水平构造应力较大的地区,中强、强地震分布较为稀疏,可能会导致断层闭锁而造成能量积累,仍有强震发生的可能。需要注意的是,图中标示的地震几乎均发生于应力突变区域,可见应力分布不均匀以及应力突变可能是导致地震频发的一个重要因素。

图 8 川西地区地壳浅层水平构造应力等值线与6级以上地震分布 Fig. 8 Distribution of ground horizontal tectonic stress field contour lines and the earthquake more than 6 degree in western Sichuan
4 结语

1) 利用重力异常数据能够快速获取水平构造应力场。获得的水平构造应力场能够基本反映川西地区地壳构造活动状况,川西地区构造活动特征可通过水平构造应力大小及方向较好地诠释。然而,通过重力资料获得的水平构造应力场虽能反映大部分水平构造应力特征,但仍有部分水平构造应力信息缺失。

2) 以静态水平构造应力场作为背景构造应力场,结合流动重力观测资料、GPS观测位移资料、地震波观测资料等,进一步开展动态水平构造应力场研究,能够更深入地探究构造应力场及重力场与地震孕育的关系,为研究断层运动机理、探索孕震过程作好铺垫。

参考文献
[1]
Zang A, Stephansson O. Stress Field of the Earth's Crust[M]. London: Springer, 2010 (0)
[2]
游永雄. 重力场转换区域构造应力场的研究[J]. 地球物理学报, 1994, 37(增1): 259-271 (You Yongxiong. A Research on the Conversion of Gravity to Tectonic Stress Fields[J]. Acta Geophysica Sinica, 1994, 37(S1): 259-271) (0)
[3]
曹学伟.基于大地测量资料的区域构造应力场反演及模型研究[D].青岛: 山东科技大学, 2010 (Cao Xuewei. Study on Inversion of Regional Tectonic Stress Field and Models Based on Geodesy Data[D]. Qingdao: Shandong University of Science and Technology, 2010) (0)
[4]
郭飞霄, 肖云, 苗岳旺. 利用EGM2008模型和地面实测重力资料反演川滇地区构造应力场[J]. 大地测量与地球动力学, 2015, 35(3): 445-448 (Guo Feixiao, Xiao Yun, Miao Yuewang. Inversion of Tectonic Stress Field in Sichuan-Yunnan Area Using EGM2008 Model and Ground Gravity Observation[J]. Journal of Geodesy and Geodynamics, 2015, 35(3): 445-448) (0)
[5]
向文, 李辉. 重力场与构造应力场内在关系的理论研究[J]. 地壳形变与地震, 1999, 19(1): 32-37 (Xiang Wen, Li Hui. Study on the Inner Relation of Gravity to Tectonic Stress Fields[J]. Crustal Deformation and Earthquake, 1999, 19(1): 32-37 DOI:10.3969/j.issn.1671-5942.1999.01.005) (0)
[6]
李忠亚, 胡敏章, 李辉. 华北地区重力势能差水平构造应力分布特征[J]. 大地测量与地球动力学, 2017, 37(9): 908-912 (Li Zhongya, Hu Minzhang, Li Hui. Research on Horizontal Tectonic Stress Distribution Contributed by Gravitational Potential Energy Differences in North China Area[J]. Journal of Geodesy and Geodynamics, 2017, 37(9): 908-912) (0)
[7]
祝意青, 闻学泽, 孙和平, 等. 2013年四川芦山MS7.0地震前的重力变化[J]. 地球物理学报, 2013, 56(6): 1 887-1 894 (Zhu Yiqing, Wen Xueze, Sun Heping, et al. Gravity Change before the Lushan, Sichuan, MS=7.0 Earthquake of 2013[J]. Chinese Journal of Geophysics, 2013, 56(6): 1 887-1 894) (0)
[8]
祝意青, 刘芳, 李铁明, 等. 川滇地区重力场动态变化及其强震危险含义[J]. 地球物理学报, 2015, 58(11): 4 187-4 196 (Zhu Yiqing, Liu Fang, Li Tieming, et al. Dynamics Variation of the Gravity Field in the Sichuan-Yunnan Region and Its Implication for Seismic Risk[J]. Chinese Journal of Geophysics, 2015, 58(11): 4 187-4 196) (0)
[9]
NGA(2010). Earth Gravitational Model 2008 (EGM2008)[EB/OL].http://earth-info.Nga.mil/GandG/wgs84/gravitymod/egm 2008/index.html (0)
[10]
Sandwell D T, Smith W H F. Global Marine Gravity from Retracked Geosat and ERS-1 Altimetry:Ridge Segmentation versus Spreading Rate[J]. J Geophy Res, 2009, 114: B01411 (0)
[11]
章传银, 郭春喜, 陈俊勇, 等. EGM2008地球重力场模型在中国大陆适用性分析[J]. 测绘学报, 2009, 38(4): 283-289 (Zhang Chuanyin, Guo Chunxi, Chen Junyong, et al. EGM2008 and Its Application Analysis in Chinese Mainland[J]. Acta Geodaetica et Cartographica Sinica, 2009, 38(4): 283-289 DOI:10.3321/j.issn:1001-1595.2009.04.001) (0)
[12]
邓起东, 高翔, 杨虎. 断块构造、活动断块构造与地震活动[J]. 地质科学, 2009, 44(4): 1 083-1 093 (Deng Qidong, Gao Xiang, Yang Hu. Fault-Block Tectonics, Active Fault-Block Tectonics and Earthquake Activity[J]. Chinese Journal of Geology, 2009, 44(4): 1 083-1 093) (0)
[13]
孟文, 陈群策, 吴满路, 等. 龙门山断裂带现今构造应力场特征及分段性研究[J]. 地球物理学进展, 2013, 28(3): 1 150-1 160 (Meng Wen, Chen Qunce, Wu Manlu, et al. Research on Segmentation and Characteristic of Tectonic Stress Field of Longmenshan Fault Zone[J]. Progress in Geophysics, 2013, 28(3): 1 150-1 160) (0)
[14]
丁开华.青藏高原中东部地区的现今地壳形变研究[D].武汉: 武汉大学, 2010 (Ding Kaihua. Research on Present-Day Crustal Deformation of Central and East Qinghai-Tibet Plateau[D]. Wuhan: Wuhan University, 2010) (0)
[15]
赵静, 江在森, 牛安福, 等. 川滇菱形块体东边界断层闭锁程度与滑动亏损动态特征研究[J]. 地球物理学报, 2015, 58(3): 872-885 (Zhao Jing, Jiang Zaisen, Niu Anfu, et al. Study on Dynamic Characteristics of Fault Locking and Fault Slip Deficit in the Eastern Boundary of the Sichuan-Yunnan Rhombic Block[J]. Chinese Journal of Geophysics, 2015, 58(3): 872-885) (0)
Research on Horizontal Tectonic Stress Field Based on EGM2008 in Western Sichuan
MAO Jinglun1,2     LI Hui1     ZHU Yiqing2     LI Zhongya1     
1. Key Laboratory of Earthquake Geodesy, Institute of Seismology, CEA, 40 Hongshance Road, Wuhan 430071, China;
2. The Second Monitoring and Application Center, CEA, 316 Xiying Road, Xi'an 710054, China
Abstract: We studied the method of computing the horizontal tectonic stress field by the earth gravity field, and then the horizontal tectonic stress field in western Sichuan is computed with the EGM2008. The results show that we can quickly obtain the tectonic stress field using gravity field data. The horizontal tectonic of the Longmenshan and Xianshuihe-Anninghe fault zones is greater, and the accumulation degree is higher than any other area, which reflects that tectonic movement is more active in these areas. The horizontal tectonic stress from the Bayanhar block to the southern China block to the Sichuan-Yunnan block shows clockwise rotation. In the Bayanhar block and southern China block boundary, the strong stop of Sichuan basin in the south China block results in the formation of Longmenshan foreland thrust zone. The tectonic stress field gradually decreases with the increase of depth, and the frequency of seismic activity decreases accordingly. The tectonic stress field is densely distributed with active tectonic activity and frequent seismic activity, which indicates that there is a corresponding relationship between seismic activity and tectonic stress accumulation.
Key words: EGM2008; gravity field; tectonic stress field; western Sichuan area