均匀弹性分层介质模型中的同震位移、应变以及应力
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
地震发生时产生的同震位移、应变以及应力变化,特别是同震应力的变化,在地震触发等问题的研究中有着重要的意义.本文进一步发展了基于均匀弹性水平层状介质,利用广义反射透射系数矩阵和离散波数计算同震位移的方法,使之可以计算相应的应变、应力以及同震库仑应力变化.可适用于多种情况,接收点可以位于地表以及地表以下,震源类型可以是剪切位错源以及拉张位错源.通过与半无限介质的解析解相比较,结果一致,验证了方法的可靠性以及计算精度,可以用于计算地震之后库仑应力变化,为判断余震分布提供参考.在计算同震位移时,使用了梯形积分与Filon积分相结合的积分方式,即提高了同震位移计算的速度,又保证了计算精度,有利于反演问题研究.
The coseismic displacement,strain and stress,especially the change of the coseismic stress,are significant in the earthquake triggering investigation.Based on the layered elastic model,we develop a method to calculate the coseismic dislocation by the generalized reflection-transmission coefficient matrix and discrete wavenumber method.We can use it to calculate the coseismic strain,stress,and the change of Coulomb stress.It can be applied in the general situation.The receiver can be located on the surface or underground.The earthquake source can be pure shear dislocation source or tensile source.For the half-space model,the result of our method coincides with the analytical solution,so the precision of our method is reliable.The method can be used to calculate the Coulomb stress after an earthquake,and the change of Coulomb stress can give a reference to estimate the distribution of the aftershocks.We combine the trapezoidal integration method with the Filon integration method to calculate the coseismic displacement.It improves the speed of the calculation and there is no accuracy loss.It is useful for the inversion research based on the coseismic displacement.
引文
[1]Stein R S,King G C P,Lin J.Change in failure stress on thesouthern San Andreas fault system caused by the 1992magnitude=7.4 Landers earthquake.Science,1992,258(5086):1328-1332.
    [2]Bouchon M.The state of stress on some faults of the SanAndreas system as inferred from near-field strong motiondata.J.Geophys.Res.,1997,102(B6):11731-11744.
    [3]沈正康,万永革,甘卫军等.东昆仑活动断裂带大地震之间的黏弹性应力触发研究.地球物理学报,2003,46(6):786-795.Shen Z K,Wan Y G,Gan W J,et al.Viscoelastic triggeringamong large earthquakes along the east Kunlun fault system.Chinese J.Geophys.(in Chinese),2003,46(6):786-795.
    [4]King G C P,Stein R S,Lin J.Static stress changes and thetriggering of earthquakes.Bull.Seism.Soc.Am.,1994,84(3):935-953.
    [5]石耀霖,曹建玲.库仑应力计算及应用过程中若干问题的讨论—以汉川地震为例.地球物理学报,2010,53(1):102-110.Shi Y L,Cao J L.Some aspects in static stress changecalculation—case study on Wenchuan earthquake.Chinese J.Geophys.(in Chinese),2010,53(1):102-110.
    [6]Toda S,Stein R S,Reasenberg P A,et al.Stress transferredby the 1995 Mw=6.9Kobe,Japan,shock:effect on aftershocksand future earthquake probabilities.J.Geophys.Res.,1998,103:24543-24565.
    [7]万永革.“地震静态应力触发”问题的研究.北京:中国地震局地球物理研究所,2001.Wan Y G.The study about“Static stress triggering”(inChinese).Beijing:Institute of Geophysics China EarthquakeAdministration,2001.
    [8]Barka A.The 17 August 1999Izmit earthquake.Science,1999,258(5435):1858-1859.
    [9]Steketee J A.On volterra′s dislocations in a semi-infiniteelastic medium.Can.J.Phys.,1958,36(2):192-205.
    [10]Maruyama T.Static elastic dislocation in an infinite and semi-infinite medium.Bull.Earthquake Res.Inst.Tokyo Univ.,1964,42:289-368.
    [11]Press F.Displacements,strains,and tilts at teleseismicdistances.J.Geophys.Res.,1965,70(10):2395-2412.
    [12]Okada Y.Surface deformation due to shear and tensile faultsin a half-space.Bull.Seism.Soc.Am.,1985,75(4):1135-1154.
    [13]Okada Y.Internal deformation due to shear and tensile faultsin a half-space.Bull.Seism.Soc.Am.,1992,82:1018-1040.
    [14]Sato R,Matsu′ura M.Static deformations due to the faultspreading over several layers in a multilayered medium,I,displacement.J.Phys.Earth,1973,21(3):227-249.
    [15]Ben-Menahem A,Singh S J.Multipolar elastic fields in alayered half space.Bull.Seism.Soc.Am.,1968,58(5):1519-1572.
    [16]Singh S J.Static deformation of a multilayered half-space byinternal sources.J.Geophys.Res.,1970,75(17):3257-3263.
    [17]Wang R,Martin F L,Roth F.Computation of deformationinduced by earthquakes in a multi-layered elastic crust:FORTRAN programs EDGRN/EDCMP.Comp.Geosci.,2003,29(2):195-207.
    [18]Israel M,Ben-Menahem A,Singh S J.Residual deformationof real Earth models with application to the Chandler wobble.Geophysical Journal of the Royal Astronomical Society,1973,32(2):219-247.
    [19]Pollitz F F.Coseismic deformation from earthquake faultingon a layered spherical earth.Geophys.J.Int.,1996,125(1):1-14.
    [20]Nur A,Mavko G.Postseismic viscoelastic rebound.Science,1974,183(4121):204-206.
    [21]Wang R,Martin F L,Roth F.PSGRN/PSCMP—a new codefor calculating co-and post-seismic deformation,geoid andgravity changes based on the viscoelastic-gravitationaldislocation theory.Comp.Geosci.,2006,32(4):527-541.
    [22]Pollitz F F.Postseismic relaxation theory on the sphericalearth.Bull.Seism.Soc.Am.,1992,82(1):422-453.
    [23]Masterlark T,Wang H F.Transient stress coupling betweenthe 1992Landers and 1999Hector Mine,California,earthquakes.Bull.Seism.Soc.Am.,2002,92(2):1470-1480.
    [24]张怀,吴忠良,张东宁等.虚拟川滇—基于千万网格并行有限元计算的区域强震演化过程数值模型设计和构建.中国科学D辑,2009,39(3):260-270.Zhang H,Wu Z L,Zhang D N,et al.Quasi Sichuan-Yunnan-Based on thousands of grid parallel finite elementcomputations strong earthquake evolution process numericalmodel design and construction.Science Chian(D)(inChinese),2009,39(3):260-270.
    [25]Kennett B L N.Reflections,rays,and reverberations.Bull.Seism.Soc.Am.,1974,64(6):1685-1696.
    [26]Kennett B L N.Seismic wave propagation in stratified media.New York:Cambridge U.Press,1983.
    [27]Bouchon M.Discrete wave number representation of elasticwave fields in three-space dimensions.J.Geophys.Res.,1979,84(B7):3609-3614.
    [28]Yao Z X,Harkrider D G.A generalized reflectiontransmission coefficient matrix and discrete wavenumbermethod for synthetic seismograms.Bull.Seism.Soc.Am.,1983,73(6A):1685-1699.
    [29]谢小碧,姚振兴.计算分层介质中位错点源静态位移场的广义反射、透射系数矩阵和离散波数方法.地球物理学报,1989,32(3):270-280.Xie X B,Yao Z X.A generalized reflection-transmitioncoefficient matrix method to calculate static displacement fieldof a stratified half-space by dislocation source.Chinese J.Geophys.(in Chinese),1989,32(3):270-280.
    [30]He Y M,Wang W M,Yao Z X.Static deformation due toshear and tensile faults in a layered half-space.Bull.Seism.Soc.Am.,2003,93(5):2253-2263.
    [31]Mallick S,Frazer L N.Practical aspects of reflectivitymodeling.Geophysics,1988,55(10):1355-1364.
    [32]纪晨,姚振兴.区域地震范围的宽频带理论地震图算法研究.地球物理学报,1995,38(4):460-468.Ji C,Yao Z X.The study of the method for broadbandregional synthetic seismogram.Chinese J.Geophys.(in Chinese),1995,38(4):460-468.
    [33]Chen X F,Zhang H M.An efficient method for computingGreen′s functions for a layered half-space at large epicentraldistances.Bull.Seism.Soc.Am.,2001,91(4):858-869.
    [34]Cotton F,Coutant O.Dynamic stress variations due to shearfaults in a plane-layered medium.Geophys.J.Int.,1997,128(3):676-688.
    [35]Aki K,Richards P G.Quantitative Seismology:Theory andMethods.San Francisco:W.H.Freeman and Co.,1980.
    [36]汪建军.同震、震后和震间应力触发.武汉:武汉大学,2010.Wang J J.Coseismic,Postseismic and Interseismic StressTriggerings(in Chinese).Wuhan:Wuhan University,2010.
    [37]周宇明,单斌,熊熊.静态应力触发中影响库仑应力变化的参数敏感性分析.大地测量与地球动力学,2008,28(5):21-26.Zhou Y M,Shan B,Xiong X.Parameters sensitivity analysisof Coulomb stress change in static stress triggering.Journalof Geodesy and Geodynamics(in Chinese),2008,28(5):21-26.

版权所有:© 2023 中国地质图书馆 中国地质调查局地学文献中心