程惠红2,,,
石耀霖2,
柳畅3,
乔学军1
1. 中国地震局地震研究所, 地震大地测量重点实验室, 武汉 430071
2. 中国科学院计算地球动力学重点实验室, 中国科学院大学, 北京 100049
3. 同济大学, 海洋地质国家重点实验室, 上海 200092
基金项目: 国家自然科学青年科学基金(41704101)、国家重点研发计划(2016YFC0600504)和中国地质科学院基本科研业务专项(JKY201901-017)共同资助
详细信息
作者简介: 董培育, 女, 助理研究员, 主要从事大地测量与地球动力学数值模拟研究.E-mail:dongpeiyu97@163.com
通讯作者: 程惠红, 女, 副教授, 主要从事地球动力学方面的研究.E-mail:chenghuihong@163.com
中图分类号: P313收稿日期:2018-07-05
修回日期:2018-09-03
上线日期:2019-08-05
Numerical inversion of regional initial tectonic stress based on Monte Carlo method-A case study of Bayan Har block
DONG PeiYu1,,CHENG HuiHong2,,,
SHI YaoLin2,
LIU Chang3,
QIAO XueJun1
1. Key Laboratory of Earthquake Geodesy, Institute of Seismology, China Earthquake Administration, Wuhan 430071, China
2. Key Laboratory of Computational Geodynamics, University of Chinese Academy of Sciences, Beijing 100049, China
3. State Key Laboratory of Marine Geology, Tongji University, Shanghai 200092, China
More Information
Corresponding author: CHENG HuiHong,E-mail:chenghuihong@163.com
MSC: P313--> Received Date: 05 July 2018
Revised Date: 03 September 2018
Available Online: 05 August 2019
摘要
摘要:构造应力场往往对地震活动性具有控制作用,应力快速集中的地方常常是地震频繁发生的地方.本文以巴颜喀拉块体及其边界断裂带近20年来的7次中强震为例,结合区域历史地震震源信息、地质背景及GPS等观测数据,利用Monte Carlo方法和库仑-摩尔破裂准则为计算依据,反演该块体的震前初始构造应力场.通过将初始应力场反演中不确定部分限定在一个合理的上下限范围内进行独立的重复性随机试验,并运用统计学方法得到了巴颜喀拉块体1997年玛尼MW7.5地震震前区域初始应力场.计算结果显示:(1)巴颜喀拉块体10 km深度处最大水平主应力方向自西向东呈顺时针旋转趋势,由NS向转变为近EW向,与浅部实测地应力数据、历史地震类型和板块运动方向吻合较好.(2)最大/最小水平主应力和二者差值自西向东均逐渐增加,最大水平主压应力值~400 MPa,最小水平主压应力值~250 MPa.差应力在昆仑山断裂带与阿尔金断裂带交汇处及甘孜—玉树断裂带西段较低(~150 MPa);在昆仑山断裂带东端和甘孜—玉树断裂带的东南段局部地区较高(~220 MPa).
关键词: 巴颜喀拉块体/
有限元数值模拟/
初始构造应力场/
Monte Carlo方法
Abstract:The regional tectonic stress controls the seismic activity, so earthquakes usually occurred where stress is rapidly concentrated. We took 7 moderate-strong earthquakes occurring in the Bayan Har block and the boundary fault zone as a case study to calculate the regional initial tectonic stress of the block. Combined with the focal mechanism of historical earthquakes, geological settings and GPS velocities, we used the Monte Carlo method and the Coulomb-Mohr fracture criterion to invert the initial tectonic stress field. We took abundant repetitive independent calculations within the reasonable range, and used statistical method to get an initial tectonic stress before the 1997 Manyi MW7.5 earthquake occurred. Our results show that at the depth of 10 km, from west to east, the direction of the maximum horizontal principal stress gradually changes from NS to EW direction in a clockwise rotation trend, which is consistent with the data from the observed shallow crustal stress, historical earthquakes and plate movement. The maximum/minimum horizontal principal stress and the their difference values gradually increase from west to east, with the maximum horizontal principal compressive stress value is about 400 MPa, and the minimum value is about 250 MPa. The differential stress value is lower (about 150 MPa) at the junction of the Kunlun fault and the Altyn fault zone and the western section of the Garzê-Yushu fault zone. However, the eastern end of the Kunlun fault zone and the southeastern segment of the Garzê-Yushu fault zone is distributed with the higher differential stress values (about 220 MPa).
Key words:Bayan Har block/
Finite element modeling/
Initial tectonic stress/
Monte Carlo method
PDF全文下载地址:
http://www.geophy.cn/data/article/export-pdf?id=dqwlxb_15111