基于MonteCarlo方法数值反演区域初始构造应力-中国电机工程学报.PDF

基于基于Monte Carlo 方法数值方法数值反演区域初始构造应力场反演区域初始构造应力场 基于基于 方法方法数值数值反演区域初始构造应力场反演区域初始构造应力场 ——以巴颜喀拉块体为例——以巴颜喀拉块体为例 ————以巴颜喀拉块体为例以巴颜喀拉块体为例 1 2* 2 3 1 董培育 ，程惠红 ，石耀霖 ，柳畅 ，乔学军 1 中国地震局地震研究所，地震大地测量重点实验室，武汉430071 2 中国科学院计算地球动力学重点实验室，中国科学院大学，北京100049 3 同济大学，海洋地质国家重点实验室，上海 200092 摘要摘要 构造应力场往往对地震活动性具有控制作用，应力快速集中的地方常常是地震频繁发 摘要摘要 生的地方。本文以巴颜喀拉块体及其边界断裂带近二十年来的7 次中强震为例，结合区域历 史地震震源信息、地质背景及GPS 等观测数据，利用Monte Carlo 方法和库仑破裂准则为计 算依据，反演该块体的震前初始构造应力场。通过将初始应力场反演中不确定部分限定在一 个合理的上下限范围内进行独立的重复性随机试验，并运用统计学方法得到了巴颜喀拉块体 1997 年玛尼M 7.5 地震震前区域初始应力场。计算结果显示：（1）巴颜喀拉块体10 km 深 w 度处最大水平主应力方向自西向东呈顺时针旋转趋势，由NS 向转变为近EW 向，与浅部实 测地应力数据、历史地震类型和板块运动方向吻合较好。(2) 最大/最小水平主应力和二者差 值自西向东均逐渐增加，最大水平主压应力值～400MPa ，最小水平主压应力～250MPa 。差 应力在昆仑山断裂带与阿尔金断裂带交汇处及甘孜玉树断裂带西段较低(～150MPa)；在昆 仑山断裂带东端和甘孜-玉树断裂带的东南段局部地区较高(最大值～220MPa) 。 关键词：巴颜喀拉块体；有限元数值模拟；初始构造应力场；关键词 Monte Carlo 方法 关键词关键词 Numerical inversion of regional initial tectonic stress based on Monte Carlo Method --A case study of Bayan Har Block 1 2* 2 3 1 DONG Pei-Yu , CHENG Hui-Hong , SHI Yao-Lin , LIU Chang , QIAOXue-Jun 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 Abstract The regional tectonic stress controls the seismic activity, so earthquakes usually occur where stress is rapidly concentrated. We take 7 moderate-strong earthquakes occurring