近断层应力变化快速计算方法及对强余震空间分布的指示意义:以M_w7.9汶川地震为例
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
构造地震一般由断层摩擦失稳所致.断层内部及周边所累积的剪切形变则通过同震滑动位移得到局部释放.因此,震后断层内部及近断层周边的静态剪切应力变化量的空间分布可通过断层面上的滑动位移分布计算得到.本文采用傅氏变换方法(FTM)计算单一有限断层同震滑移场所形成的静态剪切应力变化分布,近断层两侧的应力变化计算可由波数域内应力近似算法获得.结果表明,FTM快速有效、易于实现,有效地避免了常规应力计算中奇异值的出现.以2008年Mw7.9中国汶川大地震为例,采用前人所得有限断层滑动位移模型,得到了断层面和近断层周边准3D剪切应力分布解,并同主震后中强余震的空间分布特征作了比较.结果表明,大部分的中强余震震源位置处于剪切应力变化值为正的区域,由同震位移所产生的局部应力降峰值和均值大小同应力变化的正值大致相当,从而表明了快速且有效地计算断层内部及近断层附近的应力变化分布可以为主震后强余震发生的潜在区域提供指示意义.需要强调的是,应力变化空间分布特征的获取强烈地依赖于有限断层滑移模型解.有关滑动位移反演解的多解性对应力变化计算结果的影响,本文作了必要的讨论.
It is widely accepted that earthquake rupture occurs when the shear stress accumulated on the fault exceeds the fault frictional strength.The fault slip motion is accompanied with a partial shear strain release,and the shear stress change can be resolved from final slip distribution or co-seismic slip on the causative fault.We use the Fourier transform method(FTM)to fast calculate the shear stress changes on and near faults after the 2008 Mw7.9 Wenchuan,China,earthquake due to a heterogeneous slip distribution on the main fault and try to understand the relationship between stress change and aftershock distribution.For a single fault plane,the shear stress change on the main fault can be calculated efficiently by using the FTM,and the calculation of off-fault stress change is based on an approximate decay relation proposed by Helmstetter and Shaw.For a 2Danti-plane strain problem,we compare our numerical model with the analytical model given by Burridge and Halliday.The results show that the FTM algorithm is easy to be implemented and effectively avoid the singular values appearing in conventional stress calculation.Based on the finite fault slip model inverted from tele-seismic data for the Wenchuan event,we obtain the stress changes on and off the fault and create a 3D view of the aftershock distribution around the main fault together with the spatial stress change pattern resulted from heterogeneous fault slip.The result of 3D view of the aftershock distribution shows that the triggered events or aftershocks are mostly concentrated in the regions with a positive stress change,and the value of the maximum positive stress change caused by the co-seismic slip is roughly equal to the localized largest stress drop(stress change is negative).We calculate the seismicity rate against time following the Wenchuan main shock based on the temporal distribution of stress changes.The results my help us better understanding the aftershock duration and seismicity reduce caused by the stress drop.We emphasize that the fast and effective computations of stress changes on and off the fault immediately following a main shock can give us an insight into the source physics behind earthquake nucleation which provides us the plausible indicators of future potential areas of the moderate to strong aftershocks.In addition,the spatial distribution of stress change is strongly dependent on the finite fault slip model.So we discuss the influence of non-unique inversion of slip distribution on the stress change calculation.
引文
Aki K,Richards P G.1980.Quantitative Seismology.San Francisco:Freeman.
    Andrews D J.1974.Evaluation of static stress on a fault plane froma Green′s function.Bull.Seismol.Soc.Am.,64(6):1629-1633.
    Andrews D J.1976.Rupture velocity of plane strain shear cracks.J.Geophys.Res.,81(32):5679-5687.
    Andrews D J.1978.Coupling of energy between tectonic processesand earthquakes.J.Geophys.Res.,83(B5):2259-2264.
    Andrews D J.1980.A stochastic fault model:1.Static case.J.Geophys.Res.,85(B7):3867-3877.
    Bouchon M.1997.The state of stress on some faults of the SanAndreas System as inferred from near-field strong motion data.J.Geophys.Res.,102(B6):11731-11744.
    Burridge R,Halliday G S.1971.Dynamic shear cracks with frictionas models for shallow focus earthquakes.Geophysical JournalInternational,25(1-3):261-283.
    Chen J,Hayes G.2008.Preliminary result of the May 12,2008Mw7.9Eastern Sichuan,China earthquake(http://earthquake.usgs.gov/eqcenter/eqinthenews/2008/us2008ryan/finite_fault.php).
    Chen J H,Liu Q Y,Li S C,et al.2009.Seismotectonic study byrelocation of the Wenchuan Ms8.0earthquake sequence.ChineseJournal of Geophysics(in Chinese),52(2):390-397.
    Day S M,Yu G,Wald D J.1998.Dynamic stress changes duringearthquake rupture.Bull.Seismol.Soc.Am.,88(2):512-522.
    Dieterich J H,Linker M F.1992.Fault stability under conditions ofvariable normal stress.Geophys.Res.Lett.,19(16):1691-1694.
    Dieterich J H.1994.A constitutive law for rate of earthquakeproduction and its application to earthquake clustering.J.Geophys.Res.,99(B2):2601-2618.
    Dieterich J H,Kilgore B.1996.Implications of fault constitutiveproperties for earthquake prediction.Proceedings of the NationalAcademy of Sciences of the United States of America,93(9):3787-3794.
    Feng G C,Hetland E A,Ding X L,et al.2010.Coseismic fault slipof the 2008 Mw7.9Wenchuan earthquake estimated from InSARand GPS measurements.Geophys.Res.Lett.,37(1):L01302,doi:10.1029/2009GL041213.
    Freed A M.2005.Earthquake triggering by static,dynamic,andpostseismic stress transfer.Annual Review of Earth andPlanetary Sciences,33(1):335-367.
    Hainzl S,Marsan D.2008.Dependence of the Omori-Utsu lawparameters on main shock magnitude:Observations andmodeling.J.Geophys.Res.,113(B10):B10309,10.1029/2007JB005492.
    Harris R A,Simpson R W.1998.Suppression of large earthquakesby stress shadows:A comparison of Coulomb and rate-and-statefailure.J.Geophys.Res.,103(B10):24439-24451.
    Helmstetter A,Shaw B E.2006.Relation between stress heterogeneityand aftershock rate in the rate-and-state model.J.Geophys.Res.,111(B7):B07304,doi:10.1029/2005JB004077.
    Helmstetter A.Shaw B E.2009.Afterslip and aftershocks in therate-and-state friction law.J.Geophys.Res.,114(B1):B01308,doi:10.1029/2007JB005077.
    Hill D P,Prejean S G.2007.Dynamic triggering.//Schubert G.Treatise on Geophysics,Vol.4.Oxford:Elsevier:257-292.
    Huang Y,Wu J P,Zhang T Z,et al.2008.Relocation of the M8.0Wenchuan earthquake and its aftershock sequence.ScienceChina(Earth Sciences),51:1703-1711.
    Ide S,Takeo M.1997.Determination of constitutive relations offault slip based on seismic wave analysis.J.Geophys.Res.,102(B12):27379-27391.
    Kanamori H,Brodsky E E.2004.The physics of earthquakes.Rep.Prog.Phys.,67(8):1429-1496.
    King G C P,Stein R S,Lin J.1994.Static stress changes and thetriggering of earthquakes.Bull.Seismol.Soc.Am.,84(3):935-953.
    Lay T,Wallace T C.1995.Modern Global Seismology.New York:Academic Press.
    Liu B Y,Shi B P,Lei J S.2013.Effect of Wenchuan earthquake onprobabilities of earthquake occurrence of Lushan and surroundingfaults.Acta Seismologica Sinica(in Chinese),35(5):642-651.
    Mai P M,Beroza G C.2002.A spatial random field model tocharacterize complexity in earthquake slip.J.Geophys.Res.,107(B11):ESE10-1-ESE10-21.
    Okada Y.1985.Surface deformation due to shear and tensile faultsin a half-space.Bull.Seismol.Soc.Am.,75(4):1135-1154.
    Okada Y.1992.Internal deformation due to shear and tensile faultsin a half-space.Bull.Seismol.Soc.Am.,82(2):1018-1040.
    Parsons T.2002.Global Omori law decay of triggered earthquakes:Large aftershocks outside the classical aftershock zone.J.Geophys.Res.,107(B9):ESE9-1-ESE9-20.
    Parsons T,Velasco A A.2009.On near-source earthquaketriggering.J.Geophys.Res.,114(4):B10307,doi:10.1029/2008JB006277.
    Perfettini H,Avouac J P.2004.Postseismic relaxation driven bybrittle creep:A possible mechanism to reconcile geodeticmeasurements and the decay rate of aftershocks,application tothe Chi-Chi earthquake,Taiwan.J.Geophys.Res.,109(B2):B02304,doi:10.1029/2003JB002488.
    Perfettini H,Avouac J P.2007.Modeling afterslip and aftershocksfollowing the 1992Landers earthquake.J.Geophys.Res.,112(B7):B07409,doi:10.1029/2006JB004399.
    Peyrat S,Olsen K B.2004.Nonlinear dynamic rupture inversion ofthe 2000 Western Tottori,Japan,earthquake.Geophys.Res.Lett.,31(5):L05604,doi:10.1029/2003GL019058.
    Ripperger J,Mai P M.2004.Fast computation of static stresschanges on 2Dfaults from final slip distributions.Geophys.Res.Lett.,31(18):L18610,doi:10.1029/2004GL020594.
    Schaff D P,Beroza G C,Shaw B E.1998.Postseismic response ofrepeating aftershocks.Geophys.Res.Lett.,25(24):4549-4552.
    Scholz C H.1998.Earthquakes and friction laws.Nature,391(6662):37-42.
    Segall P.2010.Earthquake and Volcano Deformation.Princeton:Princeton University Press.
    Sladen A.2008.Preliminary result 05/12/2008(Mw7.9),eastSichuan(http://www.tectonics.caltech.edu/slip_history/2008_e_sichuan/e_sichuan.html).
    Starr A T.1928.Slip in a crystal and rupture in a solid due toshear.//Green B J.Mathematical Proceedings of the CambridgePhilosophical Society.Cambridge:Cambridge University Press,489-500.
    Stein R S.1999.The role of stress transfer in earthquakeoccurrence.Nature,402(6762):605-609.
    Stein R S,Barka A A,Dieterich J H.1997.Progressive failure onthe North Anatolian fault since 1939 by earthquake stresstriggering.Geophys.J.Int.,128(3):594-604.
    Stein S,Wysession M.2009.An Introduction to Seismology,Earthquakes,and Earth Structure.Chichester:John Wiley&Sons.
    Tong X P,Sandwell D T,Fialko Y.2010.Coseismic slip model ofthe 2008 Wenchuan earthquake derived from joint inversion ofinterferometric synthetic aperture radar,GPS,and field data.J.Geophys.Res.,115(B4):B04314,doi:10.1029/2009JB006625.
    Wang X,Yu X,Zhang W.2011.High-resolution 3D seismicimaging of the Longmenshan fault zone structure using doubledifference seismic tomography.AGU Fall Meeting Abstracts,7.
    Wu J P,Huang Y,Zhang T Z,et al.2009.Aftershock distributionof the Ms8.0 Wenchuan earthquake and three dimensional Pwave velocity structure in and around source region.ChineseJournal of Geophysics(in Chinese),52(2):320-328.
    Xu X W,Wen X Z,Ye J Q.2008.The Ms8.0 Wenchuanearthquake surface ruptures and its seismogenic structure.Seismology and Geology(in Chinese),30(3):597-629.
    Zhu A L,Xu X W,Diao G L,et al.2008.Relocation of the Ms8.0Wenchuan earthquake sequence in part:preliminary seismotectonicanalysis.Seismology and Geology(in Chinese),30(3):759-767.
    Ziv A.2012.Inference of coseismic slip via joint inversion of GPSand aftershock data:The 2004Parkfield example.J.Geophys.Res.,117(B3):B03307,doi:10.1029/2011JB008400.
    陈九辉,刘启元,李顺成等.2009.汶川Ms8.0地震余震序列重新定位及其地震构造研究.地球物理学报,52(2):390-397.
    刘博研,史保平,雷建设.2013.汶川地震对芦山地震及周边断层发震概率的影响.地震学报,35(5):642-651.
    吴建平,黄媛,张天中等.2009.汶川Ms8.0级地震余震分布及周边区域P波三维速度结构研究.地球物理学报,52(2):320-328.
    徐锡伟,闻学泽,叶建青.2008.汶川Ms8.0地震地表破裂带及其发震构造.地震地质,30(3):597-629.
    朱艾斓,徐锡伟,刁桂苓等.2008.汶川Ms8.0地震部分余震重新定位及地震构造初步分析.地震地质,30(3):759-767.

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