基于Monte Carlo方法数值反演区域初始构造应力场——以巴颜喀拉块体为例
详细信息    查看全文 | 推荐本文 |
  • 英文篇名:Numerical inversion of regional initial tectonic stress based on Monte Carlo method——A case study of Bayan Har block
  • 作者:董培育 ; 程惠红 ; 石耀霖 ; 柳畅 ; 乔学军
  • 英文作者:DONG PeiYu;CHENG HuiHong;SHI YaoLin;LIU Chang;QIAO XueJun;Key Laboratory of Earthquake Geodesy,Institute of Seismology,China Earthquake Administration;Key Laboratory of Computational Geodynamics,University of Chinese Academy of Sciences;State Key Laboratory of Marine Geology,Tongji University;
  • 关键词:巴颜喀拉块体 ; 有限元数值模拟 ; 初始构造应力场 ; Monte ; Carlo方法
  • 英文关键词:Bayan Har block;;Finite element modeling;;Initial tectonic stress;;Monte Carlo method
  • 中文刊名:DQWX
  • 英文刊名:Chinese Journal of Geophysics
  • 机构:中国地震局地震研究所地震大地测量重点实验室;中国科学院计算地球动力学重点实验室中国科学院大学;同济大学海洋地质国家重点实验室;
  • 出版日期:2019-08-12
  • 出版单位:地球物理学报
  • 年:2019
  • 期:v.62
  • 基金:国家自然科学青年科学基金(41704101);; 国家重点研发计划(2016YFC0600504);; 中国地质科学院基本科研业务专项(JKY201901-017)共同资助
  • 语种:中文;
  • 页:DQWX201908010
  • 页数:13
  • CN:08
  • ISSN:11-2074/P
  • 分类号:80-92
摘要
构造应力场往往对地震活动性具有控制作用,应力快速集中的地方常常是地震频繁发生的地方.本文以巴颜喀拉块体及其边界断裂带近20年来的7次中强震为例,结合区域历史地震震源信息、地质背景及GPS等观测数据,利用Monte Carlo方法和库仑-摩尔破裂准则为计算依据,反演该块体的震前初始构造应力场.通过将初始应力场反演中不确定部分限定在一个合理的上下限范围内进行独立的重复性随机试验,并运用统计学方法得到了巴颜喀拉块体1997年玛尼M_W7.5地震震前区域初始应力场.计算结果显示:(1)巴颜喀拉块体10km深度处最大水平主应力方向自西向东呈顺时针旋转趋势,由NS向转变为近EW向,与浅部实测地应力数据、历史地震类型和板块运动方向吻合较好.(2)最大/最小水平主应力和二者差值自西向东均逐渐增加,最大水平主压应力值~400 MPa,最小水平主压应力值~250MPa.差应力在昆仑山断裂带与阿尔金断裂带交汇处及甘孜—玉树断裂带西段较低(~150MPa);在昆仑山断裂带东端和甘孜—玉树断裂带的东南段局部地区较高(~220MPa).
        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 M_W7.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).
引文
Anderson E M.1951.The Dynamics of Faulting and Dyke Formation with Applications to Britain.Edinburgh,Oliver and Boyd,U.K.:Hafner Pub.Co.
    Cao J L,Shi Y L,Zhang H,et al.2009.Numerical simulation of GPS observed clockwise rotation around the eastern Himalayan syntax in the Tibetan Plateau.Chinese Science Bulletin,54(8):1398-1410,doi:10.1007/s11434-008-0588-7.
    Chen Q C,Feng C J,Meng W,et al.2012.Analysis of in situ stress measurements at the northeastern section of the Longmenshan fault zone after the 5.12Wenchuan earthquake.Chinese Journal of Geophysics(in Chinese),55(12):3923-3932,doi:10.6038/j.issn.0001-5733.2012.12.005.
    Chen Y T,Yang Z X,Zhang Y,et al.2013.From 2008Wenchuan earthquake to 2013Lushan earthquake.Scientia Sinica Terrae(in Chinese),43(6):1064-1072,doi:10.1360/zd-2013-43-6-1064.
    Cheng H H,Pang Y J,Dong P Y,et al.2014.Analysis of the stress environment of the 2008and 2014Yutian MS7.3earthquakes.Chinese Journal of Geophysics(in Chinese),57(10):3238-3246,doi:10.6038/cjg20141012.
    Cui X F,Xie F R,Zhao J T.2005.The regional characteristics of focal mechanism solutions in China and its adjacent areas.Seismology and Geology(in Chinese),27(2):298-307.
    Deng Q D,Cheng S P,Ma J,et al.2014.Seismic activities and earthquake potential in the Tibetan Plateau.Chinese Journal of Geophysics(in Chinese),57(7):2025-2042,doi:10.6038/cjg20140701.
    Dong P Y,Hu C B,Shi Y L.2016.Numerical simulation of longterm deformation of Tibetan Plateau and surrounding area.Seismology and Geology(in Chinese),38(2):410-422,doi:10.3969/j.issn.0253-4967.2016.02.014.
    Dong P Y,Shi Y L.2013.A discussion on“The mechanism of longdistance jumping and the migration of main active areas for strong earthquakes occurred in the Chinese continent”-Transverse isotropic“wounded element”is a better method.Chinese Journal of Geophysics(in Chinese),56(6):2133-2139,doi:10.6038/cjg20130633.
    Du J J,Chen Q C,Ma Y S,et al.2013.Faults activity and stress state in the northeast segment of Longmenshan faults zone.Progress in Geophysics(in Chinese),28(3):1161-1170,doi:10.6038/pg20130307.
    Funning G J,Parsons B,Wright T J.2007.Fault slip in the 1997Manyi,Tibet earthquake from linear elastic modelling of InSAR displacements.Geophysical Journal International,169(3):988-1008,doi:10.1111/j.1365-246X.2006.03318.x.
    Gan W J,Zhang P Z,Shen Z K,et al.2007.Present-day crustal motion within the Tibetan Plateau inferred from GPS measurements.Journal of Geophysical Research:Solid Earth,112(B8):B08416,doi:10.1029/2005JB004120.
    Goodman R E.1989.Introduction to Rock Mechanics.2nd ed.New York:Wiley.
    Guo Q L,Wang C H,Ma H S,et al.2009.In-situ hydro-fracture stress measurement before and after the Wenchuan MS8.0earthquake of China.Chinese Journal of Geophysics(in Chinese),52(5):1395-1401,doi:10.3969/j.issn.0001-5733.2009.05.029.
    Hu C B,Zhou Y J,Cai Y E.2009.A new finite element model in studying earthquake triggering and continuous evolution of stress field.Science in China Series D:Earth Sciences,52(7):994-1004,doi:10.1007/s11430-009-0082-3.
    Li Z F,Zhou B G,Ran H L.2012.Strong earthquake risk assessment of eastern segment on the East Kunlun fault in the next 100years based on paleo-earthquake data.Chinese Journal of Geophysics(in Chinese),55(9):3051-3065,doi:10.6038/j.issn.0001-5733.2012.09.023.
    Liao C T,Zhang C S,Wu M L,et al.2003.Stress change near the Kunlun fault before and after the MS8.1Kunlun earthquake.Geophysical Research Letters,30(20):2027,doi:10.1029/2003gl018106.
    Liu C,Shi Y L,Zhu B J,et al.2014.Crustal rheology control on the mechanism of the earthquake generation at the Longmen Shan fault.Chinese Journal of Geophysics(in Chinese),57(2):404-418,doi:10.6038/cjg20140207.
    Liu C,Zhu B J,Yang X L,et al.2016.Geodynamic background of the 2008 Wenchuan earthquake based on 3D visco-elastic numerical modelling.Physics of the Earth and Planetary Interiors,252:23-36,doi:10.1016/j.pepi.2016.01.003.
    Liu C L,Zheng Y,Ge C,et al.2013.Rupture process of the MS7.0Lushan earthquake,2013.Science China Earth Sciences,56(7):1187-1192,doi:10.1007/s11430-013-4639-9.
    Liu Q Y,Wu J C.2003.On numerical forecast of earthquakesthinking about the strategy for promoting earthquake prediction.Earth Science Frontiers(in Chinese),10(S1):217-224,doi:10.3321/j.issn:1005-2321.2003.z1.030.
    Meng W,Guo C B,Zhang C Y,et al.2017.In situ stress measurements and implications in the Lhasa Terrane,Tibetan Plateau.Chinese Journal of Geophysics(in Chinese),60(6):2159-2171,doi:10.6038/cjg20170611.
    Meng W,Chen Q C,Zhao Z,et al.2015.Characteristics and implications of the stress state in the Longmen Shan fault zone,eastern margin of the Tibetan Plateau.Tectonophysics,656:1-19,doi:10.1016/j.tecto.2015.04.010.
    Pang Y J,Cheng H H,Zhang H,et al.2017.Numerical modeling of crustal deformation in the eastern margin of the Bayan Har block and analysis of seismogenic environment of the 2017Jiuzhaigou earthquake.Chinese Journal of Geophysics(in Chinese),60(10):4046-4055,doi:10.6038/cjg20171030.
    Pang Y J,Zhang H,Gerya T V,et al.2018.The Mechanism and dynamics of N-S rifting in southern Tibet:insight from 3-Dthermomechanical modeling.Journal of Geophysical Research:Solid Earth,123(1):859-877,doi:10.1002/2017JB014011.
    Prieto G A,Beroza G C,Barrett S A,et al.2012.Earthquake nests as natural laboratories for the study of intermediate-depth earthquake mechanics.Tectonophysics,570-571:42-56,doi:10.1016/j.tecto.2012.07.019.
    Royden L H,Burchfiel B C,King R W,et al.1997.Surface deformation and lower crustal flow in eastern Tibet.Science,276(5313):788-790,doi:10.1126/science.276.5313.788.
    Scholz C H.1998.Earthquakes and friction laws.Nature,391(391):37-42,doi:10.1038/34097.
    Shi F,Li A,Yang X P,et al.2013.Research on late Quaternary activity of the southeastern segment of Ganzi-Yushu Fault zone.Seismology and Geology(in Chinese),35(1):50-63,doi:10.3969/j.issn.0253-4967.2013.01.004.
    Shi Y L,Zhang B,Zhang S Q,et al.2013.Numerical earthquake prediction.Physics(in Chinese),42(4):237-255,doi:10.7693/wl20130402.
    Shi Y L,Sun Y Q,Luo G,et al.2018.Roadmap for earthquake numerical forecasting in China-Reflection on the Tenth Anniversary of Wenchuan Earthquake.Chinese Science Bulletin(in Chinese),63(19):1865-1881,doi:10.1360/N972018-00335.
    Wan Y G,Shen Z K,Zeng Y H,et al.2007.Evolution of cumulative Coulomb failure stress in northeastern Qinghai-Xizang(Tibetan)Plateau and its effect on large earthquake occurrence.Acta Seismologica Sinica(in Chinese),29(2):115-129,doi:10.3321/j.issn:0253-3782.2007.02.001.
    Wang C H,Song C K,Guo Q L,et al.2014.Stress build-up in the shallow crust before the Lushan Earthquake based on the insitu stress measurements.Chinese Journal of Geophysics(in Chinese),57(1):102-114,doi:10.6038/cjg20140110.
    Wang H,Zhang G M,Shi Y L,et al.2006.Numerical simulation of movement and deformation of Qinghai-Tibet Plateau.Journal of Geodesy and Geodynamics(in Chinese),26(2):15-23,doi:10.3969/j.issn.1671-5942.2006.02.004.
    Wang Q,Qiao X J,Lan Q G,et al.2011.Rupture of deep faults in the 2008 Wenchuan earthquake and uplift of the Longmen Shan.Nature Geoscience,4(9):634-640,doi:10.1038/ngeo1210.
    Wen X Z,Du F,Zhang P Z,et al.2011.Correlation of major earthquake sequences on the northern and eastern boundaries of the Bayan Har block,and its relation to the 2008 Wenchuan earthquake.Chinese Journal of Geophysics(in Chinese),54(3):706-716,doi:10.3969/j.issn.0001-5733.2011.03.010.
    Wu M L,Zhang C S,Liao C T,et al.2005.The recent state of stress in the central Qinghai-tibet plateau according to in-situ stress measurements.Chinese Journal of Geophysics(in Chinese),48(2):327-332.
    Wu Z H,Zhou C J,Feng H,et al.2014.Active faults and earthquake around Yushu in eastern Tibetan Plateau.Geological Bulletin of China(in Chinese),33(4):419-469,doi:10.3969/j.issn.1671-2552.2014.04.003.
    Xie F R,Cui X F,Zhao J T,et al.2004.Regional division of the recent tectonic stress field in China and adjacent areas.Chinese Journal of Geophysics(in Chinese),47(4):654-662.
    Xu L S,Chen Y T.2004.The inversion of temporal and spatial rupture process of the 14Nov.2001Kunlun earthquake based on the global long-period waveform data.Science in China Series D:Earth Sciences(in Chinese),34(3):256-264.
    Yang S X,Lu Y Z,Chen L W,et al.2012a.The mechanism of long-distance jumping and the migration of main active areas for strong earthquakes occurred in the Chinese continent.Chinese Journal of Geophysics(in Chinese),55(1):105-116,doi:10.6038/j.issn.0001-5733.2012.01.010.
    Yang S X,Yao R,Cui X F,et al.2012b.Analysis of the characteristics of measured stress in Chinese mainland and its active blocks and North-South seismic belt.Chinese Journal of Geophysics(in Chinese),55(12):4207-4217,doi:10.6038/j.issn.0001-5733.2012.12.032.
    Zhang G H,Qu C Y,Shan X J,et al.2011.The coseismic InSARmeasurements of 2008Yutian earthquake and its inversion for source parameters.Chinese Journal of Geophysics(in Chinese),54(11):2753-2760,doi:10.3969/j.issn.0001-5733.2011.11.005.
    Zhang Y,Xu L,Chen Y T.2010.Fast inversion of rupture process for 14April 2010Yushu,Qinghai,earthquake.Acta Seismologica Sinica(in Chinese),32(3):361-365,doi:10.3969/j.issn.0253-3782.2010.03.011.
    Zhou Y,Wang W M,Xiong L,et al.2015.Rupture process of 12February 2014,Yutian MW6.9earthquake and stress change on nearby faults.Chinese Journal of Geophysics(in Chinese),58(1):184-193,doi:10.6038/cjg20150116.
    Zhu S B,Shi Y L.2005.Genetic algorithm-finite element inversion of topographic spreading forces and drag forces of lower crust to upper crust in Tibetan plateau.Acta Scicentiarum Naturalum Universitis Pekinesis(in Chinese),41(2):225-234,doi:10.3321/j.issn:0479-8023.2005.02.008.
    Zhu S B,Shi Y L.2007.Error analysis of strain rates from GPSmeasurements based on Monte Carlo method.Chinese Journal of Geophysics(in Chinese),50(3):806-811.
    Zhu S B,Yuan J,Miao M.2017.Dynamic mechanisms for supershear rupture processes of the Yushu earthquake(MS=7.1).Chinese Journal of Geophysics(in Chinese),60(10):3832-3843,doi:10.6038/cjg20171013.
    Zoback M D,Harjes H P.1997.Injection-induced earthquakes and crustal stress at 9km depth at the KTB deep drilling site,Germany.Journal of Geophysical Research:Solid Earth,102(B8):18477-18491,doi:10.1029/96jb02814.
    曹建玲,石耀霖,张怀等.2009.青藏高原GPS位移绕喜马拉雅东构造结顺时针旋转成因的数值模拟.科学通报,54(2):224-234.
    陈群策,丰成君,孟文等.2012.5.12汶川地震后龙门山断裂带东北段现今地应力测量结果分析.地球物理学报,55(12):3923-3932,doi:10.6038/j.issn.0001-5733.2012.12.005.
    陈运泰,杨智娴,张勇等.2013.从汶川地震到芦山地震.中国科学:地球科学,43(6):1064-1072,doi:10.1360/zd-2013-43-6-1064.
    程惠红,庞亚瑾,董培育等.2014.于田2008年和2014年两次MS7.3地震孕育的应力环境.地球物理学报,57(10):3238-3246,doi:10.6038/cjg20141012.
    崔效峰,谢富仁,赵建涛.2005.中国及邻区震源机制解的分区特征.地震地质,27(2):298-307.
    邓起东,程绍平,马冀等.2014.青藏高原地震活动特征及当前地震活动形势.地球物理学报,57(7):2025-2042,doi:10.6038/cjg20140701.
    董培育,石耀霖.2013.关于“用单元降刚法探索中国大陆强震远距离跳迁及主体活动区域转移”的讨论---横向各向同性“杀伤单元”才是更好的途径.地球物理学报,56(6):2133-2139,doi:10.6038/cjg20130633.
    董培育,胡才博,石耀霖.2016.青藏高原及周边区域地表长期变形数值模拟.地震地质,38(2):410-422,doi:10.3969/j.issn.0253-4967.2016.02.014.
    杜建军,陈群策,马寅生等.2013.龙门山断裂带东北段地应力状态与断裂活动性研究.地球物理学进展,28(3):1161-1170,doi:10.6038/pg20130307.
    郭啟良,王成虎,马洪生等.2009.汶川MS8.0级大震前后的水压致裂原地应力测量.地球物理学报,52(5):1395-1401,doi:10.3969/j.issn.0001-5733.2009.05.029.
    胡才博,周一杰,蔡永恩.2009.如何用有限元新模型研究地震触发和应力场连续演化.中国科学D辑:地球科学,39(5):546-555.
    李正芳,周本刚,冉洪流.2012.运用古地震数据评价东昆仑断裂带东段未来百年的强震危险性.地球物理学报,55(9):3051-3065,doi:10.6038/j.issn.0001-5733.2012.09.023.
    柳畅,石耀霖,朱伯靖等.2014.地壳流变结构控制作用下的龙门山断裂带地震发生机理.地球物理学报,57(2):404-418,doi:10.6038/cjg20140207.
    刘启元,吴建春.2003.论地震数值预报---关于我国地震预报研究发展战略的思考.地学前缘,10(S1):217-224,doi:10.3321/j.issn:1005-2321.2003.z1.030.
    孟文,郭长宝,张重远等.2017.青藏高原拉萨块体地应力测量及其意义.地球物理学报,60(6):2159-2171,doi:10.6038/cjg20170611.
    庞亚瑾,程惠红,张怀等.2017.巴颜喀拉块体东缘形变及九寨沟地震孕震环境数值分析.地球物理学报,60(10):4046-4055,doi:10.6038/cjg20171030.
    石峰,李安,杨晓平等.2013.甘孜-玉树断裂带东南段晚第四纪活动性研究.地震地质,35(1):50-63,doi:10.3969/j.issn.0253-4967.2013.01.004.
    石耀霖,孙云强,罗纲等.2018.关于我国地震数值预报路线图的设想---汶川地震十周年反思.科学通报,63(19):1865-1881,doi:10.1360/N972018-00335.
    石耀霖,张贝,张斯奇等.2013.地震数值预报.物理,42(4):237-255,doi:10.7693/wl20130402.
    万永革,沈正康,曾跃华等.2007.青藏高原东北部的库仑应力积累演化对大地震发生的影响.地震学报,29(2):115-129,doi:10.3321/j.issn:0253-3782.2007.02.001.
    王成虎,宋成科,郭启良等.2014.利用原地应力实测资料分析芦山地震震前浅部地壳应力积累.地球物理学报,57(1):102-114,doi:10.6038/cjg20140110.
    王辉,张国民,石耀霖等.2006.青藏活动地块区运动与变形特征的数值模拟.大地测量与地球动力学,26(2):15-23,doi:10.3969/j.issn.1671-5942.2006.02.004.
    闻学泽,杜方,张培震.2011.巴颜喀拉块体北和东边界大地震序列的关联性与2008年汶川地震.地球物理学报,54(3):706-716,doi:10.3969/j.issn.0001-5733.2011.03.010.
    吴满路,张春山,廖椿庭等.2005.青藏高原腹地现今地应力测量与应力状态研究.地球物理学报,48(2):327-332.
    吴中海,周春景,冯卉等.2014.青海玉树地区活动断裂与地震.地质通报,33(4):419-469,doi:10.3969/j.issn.1671-2552.2014.04.003.
    谢富仁,崔效锋,赵建涛等.2004.中国大陆及邻区现代构造应力场分区.地球物理学报,47(4):654-662.
    许力生,陈运泰.2004.从全球长周期波形资料反演2001年11月14日昆仑山口地震时空破裂过程.中国科学D辑:地球科学,34(3):256-264.
    杨树新,陆远忠,陈连旺等.2012a.用单元降刚法探索中国大陆强震远距离跳迁及主体活动区域转移.地球物理学报,55(1):105-116,doi:10.6038/j.issn.0001-5733.2012.01.010.
    杨树新,姚瑞,崔效锋等.2012b.中国大陆与各活动地块、南北地震带实测应力特征分析.地球物理学报,55(12):4207-4217,doi:10.6038/j.issn.0001-5733.2012.12.032.
    张国宏,屈春燕,单新建等.2011.2008年MS7.1于田地震InSAR同震形变场及其震源滑动反演.地球物理学报,54(11):2753-2760,doi:10.3969/j.issn.0001-5733.2011.11.005.
    张勇,许力生,陈运泰.2010.2010年4月14日青海玉树地震破裂过程快速反演.地震学报,32(3):361-365,doi:10.3969/j.issn.0253-3782.2010.03.011.
    周云,王卫民,熊林等.2015.2014年2月12日MW6.9于田地震震源破裂过程及对周围断层的应力影响.地球物理学报,58(1):184-193,doi:10.6038/cjg20150116.
    朱守彪,石耀霖.2005.青藏高原地形扩展力以及下地壳对上地壳的拖曳力的遗传有限单元法反演.北京大学学报(自然科学版),41(2):225-234,doi:10.3321/j.issn:0479-8023.2005.02.008.
    朱守彪,石耀霖.2007.基于Monte Carlo方法的由GPS观测计算地应变率的误差分析.地球物理学报,50(3):806-811.
    朱守彪,袁杰,缪淼.2017.青海玉树地震(MS=7.1)产生超剪切破裂过程的动力学机制研究.地球物理学报,60(10):3832-3843,doi:10.6038/cjg20171013.