摘要
在复杂区域应力场的作用下,大范围扩容各向异性(extensive-dilatancy anisotropy,EDA)造成的裂隙定向排列形成了地壳内部介质的不均匀性。因此,地壳介质各向异性的描述不仅限于横向各向同性(transversely isotropic,TI)的理论。我们首先利用自相关函数随机扰动得到水平向小尺度速度不均匀,然后利用这种不均匀性描述地壳介质中EDA裂隙定向排列所造成的介质各向异性,利用数值模拟弹性波在该介质中的传播过程验证了利用介质的不均匀性构建介质的各向异性是一个有效、可靠的方法。可以将此方法应用到复杂地壳介质结构中,了解地震波在复杂地壳各向异性介质中的传播特征。
Under the effect of complicated local stress,the aligned micro-cracks from extensive-dilatancy anisotropy(EDA)are the main anisotropic mechanism in the upper crust.Therefore,the seismic anisotropy study is not only based on transversely isotropic(TI)medium.In this study,the EDA micro-cracks are regarded as the directional arrangement of molecular scale heterogeneity and be described by heterogeneous medium exercising an appropriate smooth disturbances autocorrelation function.According to the simulation of seismic waves by pseudospectral time-domain(PSTD)method in this equivalent anisotropic model,the directional heterogeneity is an effective and reliable method to describe anisotropy essentially,and can be extensively used in geological structure research and analysis of the characteristics of seismic wave propagation.
引文
[1]Crampin S,Evans R,Atkinson B K.Earthquake prediction:A new physical basis[J].Geophys J R Astron Soc,1984a,76:147-15.
[2]Crampin S,Gao Y.Two species of microcracks[J].Applied Geophysics.2014,11(1):1-8.
[3]高原,郑斯华,孙勇.唐山地区地壳裂隙各向异性[J].地震学报,1995,17(3):283-293.
[4]Gao Y,Wu J,Fukao Y,et al.Shear-wave splitting in the crust in North China:stress,faults and tectonic implications[J].Geophys J Int,2011,187(2):642-654.DOI:10.1111/j.1365-246X.2011.05200.x.
[5]Peng Z G,Ben-Zion Y.Systematic analysis of crustal anisotropy along the Karadere-Duzce branch of the North Anatolian fault[J].Geophys J Int,2004,159:253-272.
[6]Peng Z G,Ben-Zion Y.Spatiotemporal variations of crustal anisotropy from similar events in aftershocks of the 1999 M7.4Izmit and M7.1Düzce,Turkey,earthquake sequences[J].Geophys J Int,2005,160:1 027-1 043.
[7]Cochran E S,Li Y G,Vidale J E.Anisotropy in the shallow crust observed around the San Andress Fault before and after the 2004 M6.0Parkfield earthquake[J].Bull Seism Soc Am,2006,96:S364-S375.
[8]石玉涛,高原,张永久,等.松潘—甘孜地块东部、川滇地块北部与四川盆地西部的地壳剪切波分裂[J].地球物理学报,2013,56(2):481-494.
[9]石玉涛,高原,赵翠萍,等.汶川地震余震序列的地震各向异性[J].地球物理学报,2009,52(2):398-407.
[10]Wu J,Gao Y,Chen,Y T.Shear-wave splitting in the crust beneath the southeast Capital area of North China[J].J Seismol,2009,13:277-286.
[11]Maupin V,Park J.Chapter 1.09Theory and Observations-Wave Propagation in Anisotropic Media.In:Treatise on Geophysics,Vol1:Seismology and Structure of the Earth[M].(Eds.A.M.Dziewonski,B.Romanowicz),Elsevier,2007,289-321.
[12]Crampin S.Effective anisotropic elastic-constants for wave propagation through cracked solids[J].Geophys J R astron Soc,1984b,76:135-145.
[13]Orszag S A.Comparison of pseudospectral and spectral approximation[J].Stu App1 Math,1972,51:253-259.
[14]Tessmer E,Kosloff D.3-D elastic modeling with surface topography by a Chebyshev spectral method[J].Geophysics,1994,59:464-473.
[15]Wang Y,Takenaka H,Furumura T.Modelling seismic wave propagation in a two-dimensional cylindrical whole earth model using the pseudospectral method[J].Geophys J Int,2001,145:689-708.
[16]Huang Q H,Li Z H,Wang Y B.A parallel 3-D staggered grid pseudospectral time domain method for ground-penetrating radar wave simulation[J].J Geophys Res,2010,115:B12101,doi:10.1029/2010JB007711.
[17]Liu Q H.The PSTD algorithm:a time-domain method requiring only two cells per wave-length,Microwave Opt[J].Tech Lett,1997,15(3):158-165.
[18]Liu L B,Xiao L,Liu H,et al.Numerical simulation of the effect of a DC electric field on seismic wave propagation with the pseudospectral time domain method[J].Pure and Applied Geophysics,2006,163(9):1 893-1 913.
[19]侯安宁,等.各向异性介质中弹性波传播特征的伪谱法模拟研究[J].石油物探,1995,34(2):36-45.
[20]Ergintav S,Canitez N.Modeling of multi-scale media in discrete form[J].J Seis Expl,1997,6:77-96.
[21]Furumura T,Kennett B L N,Furumura M.Seismic wave field calculation for laterally heterogeneous whole earth models using the pseudospectral method[J].Geophys J Int,1998,135:845-860.
[22]Thomsen L.Weak elastic anisotropy[J].Geophysics,1986,51(10):1 954-1 996.
[23]Andreas Rüger.P-wave reflection coefficients for transversely isotropic models with vertical and horizontal axis of symmetry[J].Geophysicists.1997,62(3);713-722.