2003年大姚6.2级、6.1级地震序列震源位置及震源区速度结构的联合反演

采用震源位置和速度结构的联合反演方法确定2003年大姚6.2级、6.1级地震序列的分布和震源区的速度结构.该方法首先确定研究区的速度结构,然后在该速度结构的基础上对地震重新定位.结果表明:① 大姚6.2级地震序列明显分为地震较密集的东南段和地震较稀疏的北西段,而且这两段存在一定的错动,速度结构也显示出高、低速交界带在错动处具有转折特点;② 大姚6.2级地震序列在深度上呈"V"字形分布,这与该地区呈倒三角形分布的低速体相对应,即余震主要分布在高、低速交界带附近;③ 大姚6.2级、6.1级地震的破裂区都位于倒三角形低速体的东南侧,6.2级地震沿高、低速交界带往下破裂,而6.1级地震沿高、低速交界带往上破裂;④ 6 km深度的速度分布显示震源区具有高、低速度体交替呈四象限分布特征,而且余震主要沿北西向的高低速交界带分布;⑤ 大姚6.2级主震的东南侧的高、低速交界带处的低速体速度值较北西侧低速值高,其对6.2级地震往东南方向破裂具有阻碍作用,即东南侧处于能量积累状态,有利于6.1级地震的发生.     

利用尾波包络线反演方法研究伽师强震群区地壳的非均匀结构

以单次散射模型、球形辐射源、一维线性速度模型为基础计算理论包络线,通过尾波归一化方法把理论与观测包络线结合起来形成非线性方程组,用迭代最小二乘方法反演伽师强震群区的散射系数空间分布.结果表明,伽师强震群区下方中上地壳中存在明显的非均匀结构.在10~25 km各深度图像上,研究区中部都存在一低散射系数区,这个低散射系数区在10 km、15 km深度被高散射系数区包围,在20 km深度仅其东侧存在高散射系数区.1997年的7次M_s≥6.0级地震发生在研究区中部的低散射系数区内或低散射系数区向高散射系数区过渡带的低散射系数一侧.     

2003年大姚6.2、6.1级地震构造应力环境及静态应力触发研究

利用2003年7月21日、10月16日大姚6.2、6.1级地震余震序列的震源机制解资料,采用滑动方向拟合法反演出大姚两次强震的构造应力,定量给出其构造应力变化量。以反演得到的2003年大姚6.2、6.1级地震构造应力为基础,分别计算两次强震后产生的静态库仑破裂应力变化,研究大姚6.2级地震引起的应力变化是否触发了6.1级地震。结果显示:两次地震发生后震源区构造应力场最大主压应力S1作用由北东向近水平作用为主转为北东东向水平作用,6.2级地震对6.1级地震有明显的触发作用。     

2003年云南大姚6.2、6.1级地震序列特征分析及地震触发研究

本文采用双差定位法对2003年7月21日、10月16日云南大姚先后发生的6.2、6.1级地震及其余震序列进行精定位。结果表明,两个主震的震中位置相当接近,余震序列的空间分布呈明显的线性分布,地震序列精定位后的结果同震源机制解吻合很好,两个地震的发震断层均是近乎直立的右旋走滑断层。此次大姚6.2、6.1级双震,和云南地区大多数其它双震有所不同,余震分布并不呈现共轭分布,而近乎一条直线,两者间有部分重叠,表明两者受控于同一发震断层。相对主震而言,余震都呈不均匀分布,近乎一种单侧分布,6.2级地震的余震集中于主震的西北方向,6.1级地震的余震集中于主震的东南方向,并且6.1级地震是6.2级地震在同一断层向东南端延伸的另一次破裂。库仑破裂静应力变化Δσf研究结果表明,7月21日6.2级地震对10月16日6.1级地震有明显的触发作用,两个主震对断层外余震的发生均有不同程度的应力触发作用,并且2000年姚安6.5级地震对大姚两次地震的发生也有触发作用。     

2014年川滇交界东段鲁甸M_S 6.5地震序列重新定位及发震构造分析

为了研究鲁甸M_S 6.5地震的余震空间分布特征和探讨震区深部发震构造,本文利用四川、云南省固定地震台网观测资料和震后流动监测资料,采用双差定位方法对2014年8月3日至11日期间发生的鲁甸M_S 6.5主震及506个M_l≥1.5余震序列进行重新定位,获得了359个事件的重定位结果.研究结果显示:鲁甸主震的震源深度为13.8 km,与震源破裂过程显示的初始破裂深度较为接近,余震序列分布纵剖面还显示出余震主要分布在主震上方,且余震序列呈现出近EW向和NW向的不对称共轭状分布,优势分布方向为NW向的地震序列延展范围约20 km左右,并与昭通—鲁甸断裂相交,近EW向地震序列的延展长度约16 km左右.其中,在垂直NW向密集条带的震源深度剖面上,余震序列呈近于直立的条带状聚集,且余震序列震源深度分布的优势区间为3~15 km,较浅的震源深度说明了地震事件大多发生在脆性上地壳内部.地震重定位结果还表明了发震断层呈现高倾角分布的特征,结合震源机制解和地面地震地质调查及科考成果,综合表明和论证了此次鲁甸MS6.5地震的发震断层与NW向的包谷垴—小河断裂密切相关.     

2013年8月香格里拉德钦—得荣M_S5.9地震序列震源机制与应力场特征

利用中国区域台网地震波形记录,采用CAP方法反演了香格里拉德钦(位于云南省)—得荣(属于四川省)2013年8月28日MS5.1、8月31日MS5.9地震及8次MS4余震的震源双力偶断层面解和震源质心深度.结合震区地质构造、余震分布、烈度分布、动力学背景等资料,分析了此次地震序列的震源机制和应力场特征.反演结果表明,此次地震序列为节面倾角倾斜的正断层型地震,发震断层为NWW向活动构造带.序列中最大地震MS5.9和次大地震MS5.1地震的破裂节面分别为走向299°、倾角53°、滑动角-73°;走向290°、倾角55°、滑动角-72°.震源区受到强烈的水平拉张力、垂直挤压力作用.MS5.9地震后续余震T、P轴方位角随时间变化强烈,表明MS5.9地震后震源区应力调整作用明显.震源区应力场反演结果显示,地震发生的构造带上最大主拉应力为NNE-SSW向,最大主压应力为NW-SE向,与GPS观测所反映的地表最大主应力分布方向基本一致,表明震源区的应力状态可能主要受到背景大尺度构造应力场的控制.此次地震序列填充了川滇地区震源机制及应力场的空间分布图像,1976年以来可靠的震源机制解资料表明香格里拉次级块体是川滇块体及周边区域显著的拉张作用区域.香格里拉次级块体和保山次级块体正断层地震的断层节面及震源应力轴分布的空间变化,与GPS观测反映的地表最大主拉应力分布较一致,其空间分布特征反映了在青藏高原物质挤出背景下,块体之间相互作用、地势差异等作用对构造活动的影响.     

1999年山西大同Ms 5.6地震的震源断层

大同震区先后在 1989、1991和 1999年发生MS >5地震 ,利用大同遥测地震台网的记录资料进行比较精确的地震序列震源定位 ,结合宏观烈度分布和震源机制解资料 ,详细地分析对比了 3次子序列的异同。结果显示 ,1999年MS5 .6地震的震源断层是走向NWW、长 16km、宽12km、埋深 5km以下、倾角近直立的左旋走滑断层。而前 2个子序列是NNE为主的右旋走滑断层活动所致 ,表明地震破裂方向发生了变化。这种 2个以上方向先后出现、并且强弱有别的地震破裂是普遍存在的 ,表明震源环境的复杂程度与地震序列的类型有关。虽然震区存在NE向的大王村断裂和NW向的团堡断裂 ,但目前没有证据说明震源断层和 2条构造断层连通。 3次子序列的震源断层都是走滑断层 ,也和 2条构造正断层有别。 1999年的子序列可能属于新破裂。     

2008年攀枝花6.1级地震序列精定位

利用四川和云南区域数字地震台网震相到时资料,并结合Hypo2000+Hypo DD对2008年8月30日攀枝花MS6.1地震序列进行了定位。定位结果表明,地震序列的震中在空间上呈近SN向展布,余震密集区长度约为30km,主震震源深度约为14km。序列深度的分布范围主要为0~4、5~20km,而4~5km范围显示为明显的少震层。沿序列长轴的深度剖面显示,余震区中段存在1个不规则的少震的"空区",为1955年634级地震的破裂区,该地震发生在空区的南端。分析认为,2008年8月30日攀枝花MS6.1地震是因汶川8.0级地震后的应力调整造成未破裂的小凹凸体发生破裂所致。余震密集分布区沿垂直于破裂长轴的两个剖面则显示了在其北端地震震源分布更深,且断层面向NW倾斜,与已知的红格断裂的倾向一致。     

2003年大姚6.2和6.1级地震前三维波速结构的演化

尝试在年度尺度上对云南地区地壳速度结构进行成像,以确定2003年大姚地震前三维地壳速度结构的演化,并对其发震成因的介质物性变化进行探讨.检测板分辨试验显示,在年度尺度上,大姚附近区域15km深度上节点解的分辨率在0.6左右;误差分析显示,在震源附近的误差(约0.02km/s)远小于速度变化的幅度(约0.15km/s).研究结果表明,大姚地震前震源区附近形成一条北北西向的高、低波速交界带,其走向与两次大姚地震震源机制解的走向以及通过余震精定位确定的断层走向基本吻合,且震源位于交界带的高波速一侧.此外,大姚地震震后在震源下方形成一低速体.本文得到的大姚地震前的波速结构演化可为探讨其孕震、发震条件提供约束.      

2021年5月21日漾濞MS6.4地震震源区三维P和S波速度结构与地震重定位研究

2021年5月21日漾濞MS6.4地震震源区位于川滇块体西边界的维西—乔后—巍山断裂西侧,该地区近些年来发生了多次中强地震,地震活动较为活跃.对漾濞地震序列重定位和漾濞震源区及邻区的地壳精细结构研究,有助于深入理解漾濞地震的孕震环境、发震机理和破裂过程.本文基于2008年 1月1日到2021 年6 月3日区域固定台站接收到的 36938条Pg和 32111 条 Sg波到时数据,采用新发展的三重差地震层析成像算法(tomoTD)开展了漾濞MS 6.4地震震源区三维速度结构成像与地震重定位研究.结果显示:(1)余震活动主要集中在维西—乔后—巍山断裂的西侧,整体呈现沿北北西向的条带状分布,结合已有走滑型震源机制解特征,揭示了北西向隐伏断裂是发震断层,其北西段表现为倾角较陡、结构相对简单的走滑断裂,南东段由两条分支断裂组成.(2)主震的发生及地震序列分布与地壳速度结构不均匀性有着密切的关系.主震及4级以上的地震发生在高速边界上或高低速过渡区域,余震主要发生在低速、高VP/VS 区,主震上方与下方均显示高VP/VS 异常,推测在区域构造应力场的作用下,应力在孕震区的刚性介质中积累,中下地壳流体(或者部分熔融地壳物质)侵入发震断层区,弱化了漾濞 6.4 级地震的主震区.另外,余震东南侧的低VP/VS 区可能代表介质刚性强,可能阻碍了余震向南东方向继续扩展.(3)结合2013 洱源 5.5 级地震研究结果,推测维西—乔后—巍山断裂西侧可能存在着较大的北北西向隐伏断层.     

云南地区散射衰减、吸收衰减及尾波衰减的综合研究

利用云南省地震局遥测台网记录的地震波形资料,采用多流逝时间窗分析方法,得到了云南地区不同台站的散射衰减和吸收衰减,并通过单次散射模型得到了不同台站的尾波衰减.研究结果显示,云南地区的散射衰减、吸收衰减、总衰减和尾波衰减空间分布以金沙江—红河断裂带为界,均表现出西部低,东部高的特征.云南地区地震波衰减以吸收为主,2~4 Hz、4~8 Hz、8~16 Hz三个频带各台站的吸收衰减Q^－1_i均大于散射衰减Q^－1_c,1~2 Hz除个别台站的吸收衰减明显小于散射衰减,其它台站均大于或接近散射衰减.尾波衰减Q^－1_c在1~2 Hz稍大于吸收衰减而小于总衰减,在2~4 Hz、4~8 Hz、8~16 Hz与总衰减接近,个别台站大于总衰减Q^－1_t.云南地区散射衰减与频率的关系为Q^－1_s=0.0138f^－1.61,吸收衰减与频率的关系为Q^－1_i=0.0119f^－0.86,总衰减与频率的关系为Q^－1_t=0.0247f^－1.1,尾波衰减与频率的关系为Q^－1_c=0.0129f^－0.71.与世界其它地区相比,云南地区的散射衰减处于中等偏低的水平,吸收衰减处于中等偏高的水平.     

Acoustic scattering in the coastal ocean at Monterey Bay,CA, USA: Fine-scale vertical structures

Fine-structure (centimeters to meters) in vertical profiles of acoustic volume scattering strength is a common and ecologically significant characteristic of the coastal marine water column. The processes that create these structures modify the availability of food and exposure to predation for secondary producers at and below spatial and temporal scales that characterize their daily ambits. Thin acoustic scattering layers may persist for weeks at a particular coastal location, but they may also appear and disappear in only a few hours. These layers are usually evidence of mesozooplankton and micronekton having aggregated at peaks, gradients, or boundaries in the vertical distribution of various water column properties that characterize a marine ecosystem. The behaviors of both predators and prey are implicated in the generation of complex, time-dependent patterns of fine-structure. Specifically, diel vertical migration to layers of prey, isolumes, isotherms, isopycnals, or chemoclines can be responsible for the nighttime formation of thin acoustic layers. Physical processes such as horizontal shear, internal waves, water mass intrusions, tidal forcing, wind mixing of the upper water column, and horizontal advection modify acoustic scattering patterns by changing the vertical distributions of organisms that scatter sound.     

Backus averaging,scattering and drift1

Differences between traveltimes from sonic to seismic frequencies, commonly known as drift, can be attributed to a combination of multiple scattering and absorption. The portion due to scattering can be estimated directly by calculating synthetic seismograms from sonic logs. A simple alternative approach is suggested by the long-wave equivalent averaging formulae for the effective elastic properties of a stack of thin layers, which gives the same traveltime delays as the low-frequency limit of the scattering dispersion. We consider the application of these averaging formulae over a frequency-dependent window with the hope of extending their use to frequencies higher than those allowed by the original validity conditions. However, comparison of the time delay due to window-averaging with the scattering dispersion predicted by the O'Doherty-Anstey formula reveals that it is not possible to specify a form of window that will fit the dispersion across the spectrum for arbitrary log statistics. A window with a width proportional to the wavelength squared matches the behaviour at the low-frequency end of the dispersive range for most logs, and allows an almost exact match of the drift across the entire spectrum for exponential correlation functions. We examine a real log, taken from a hole in nearly plane-layered geology, which displays strong quasi-cyclical variations on one scale as well as more random, smaller-scale fluctuations. The details of its drift behaviour are studied using simple models of the gross features. The form of window which gave a good theoretical fit to the dispersion for an exponential log correlation function can only fit the computed drift at high or low frequencies, confirming that there are at least two significant scale-lengths of fluctuation. A better overall fit is obtained for a window whose width is proportional to the wavelength. The calculated scattering drift is significantly less than that observed from a vertical seismic profile, but the difference cannot be wholly ascribed to absorption. This is because the source frequency of the sonic tool is not appropriate for its resolution (receiver spacing) so that the scattering drift from sonic to seismic frequencies cannot be fully estimated from the layer model derived from the log.     

2.5D modelling, inversion and angle migration in anisotropic elastic media

2.5D modelling approximates 3D wave propagation in the dip‐direction of a 2D geological model. Attention is restricted to raypaths for waves propagating in a plane. In this way, fast inversion or migration can be performed. For velocity analysis, this reduction of the problem is particularly useful. We review 2.5D modelling for Born volume scattering and Born–Helmholtz surface scattering. The amplitudes are corrected for 3D wave propagation, taking into account both in‐plane and out‐of‐plane geometrical spreading. We also derive some new inversion/migration results. An AVA‐compensated migration routine is presented that is simplified compared with earlier results. This formula can be used to create common‐image gathers for use in velocity analysis by studying the residual moveout. We also give a migration formula for the energy‐flux‐normalized plane‐wave reflection coefficient that models large contrast in the medium parameters not treated by the Born and the Born–Helmholtz equation results. All results are derived using the generalized Radon transform (GRT) directly in the natural coordinate system characterized by scattering angle and migration dip. Consequently, no Jacobians are needed in their calculation. Inversion and migration in an orthorhombic medium or a transversely isotropic (TI) medium with tilted symmetry axis are the lowest symmetries for practical purposes (symmetry axis is in the plane). We give an analysis, using derived methods, of the parameters for these two types of media used in velocity analysis, inversion and migration. The kinematics of the two media involve the same parameters, hence there is no distinction when carrying out velocity analysis. The in‐plane scattering coefficient, used in the inversion and migration, also depends on the same parameters for both media. The out‐of‐plane geometrical spreading, necessary for amplitude‐preserving computations, for the TI medium is dependent on the same parameters that govern in‐plane kinematics. For orthorhombic media, information on additional parameters is required that is not needed for in‐plane kinematics and the scattering coefficients. Resolution analysis of the scattering coefficient suggests that direct inversion by GRT yields unreliable parameter estimates. A more practical approach to inversion is amplitude‐preserving migration followed by AVA analysis. SYMBOLS AND NOTATION A list of symbols and notation is given in Appendix D .     

Strongly Frequency-Dependent Intrinsic and Scattering Q for the Western Margin of the Sichuan Basin, China, as Revealed by a Phenomenological Model

Through applying the energy-flux model to teleseismic P-wave data, seismological scattering Q and intrinsic Q for the western margin of the Sichuan basin in southwestern China are estimated separately. The results exemplify the frequency-dependent Q values in the seismic frequency band. The scattering Q is heavily frequency-dependent and reaches a minimum at ~1.5 Hz. The intrinsic Q is nearly invariant for frequencies above 2 Hz and rapidly increases at lower frequencies. A major advantage of the energy-flux model is that it is a phenomenological model based on the conservation of energy and requires no geometrical spreading correction. The results are good counterexamples against the claim of a frequency-independent seismological Q.     

罗布泊湖相沉积物极化散射特性及其环境意义

通过对罗布泊古湖盆区域("大耳朵")湖相沉积物含水量与含盐量、介电性质、地表微地貌等理化参量的分析,本文提出了罗布泊古湖盆区域的散射过程与次地表散射结构;利用全极化测量技术(PolSAR, Polarimetric Synthetic Aperture Radar),本文对不同纹理区域进行了物理散射机制的定量化剖析,同极化相关分析与极化度推算结果表明整个"大耳朵"区域后向散射强烈,亮条带比暗条带的散射机制更为复杂,地表微地貌状态是造成"大耳朵"特征的表观原因.在此基础上,本文利用极化分解技术进一步分离不同散射机制分量,通过定量比较各部分贡献对于"大耳朵"特征的表征能力,提出次地表介质性质(如含盐量)是"大耳朵"特征的根本性成因,相关性可达0.784,同时分析了次地表高含水含盐湖相沉积物在SAR响应中的重要性及其与地表形态之间的动力学形成机制.最后,本文论证了合成孔径雷达(SAR, Synthetic Aperture Radar)在罗布泊地区的穿透能力,预示着SAR遥感技术在干旱区环境研究中的应用潜力.     

Constant Q and a Fractal, Stratified Earth

— Frequency-dependent measurements of the quality factor Q typically show a constant behaviour for low frequencies and a positive power-law dependence for higher frequencies. In particular, the constant Q pattern is usually explained using intrinsic attenuation models due to anelasticity with either a single or multiple superposed relaxation mechanisms — each with a particular resonance peak.¶However, in this study, I show using wave localisation theory that a constant Q may also be due to apparent attenuation due to scattering losses. Namely, this phenomenon occurs if the earth displays fractal characteristics. Moreover, if fractal characteristics exist over a limited range of scales only, even an absorption band can be created—in accordance with observations. This indicates that it may be very difficult to distinguish between intrinsic and scattering attenuation on the basis of frequency-dependent measurements of the quality factor only.     

Solar almucantar radiance and aerosol scattering: Springtime data at Kensington,Australia

Solar irradiance, skylight, and aureole intensities were measured early in spring 1981 in Kensington, Australia. The aerosol optical thickness, the aureole intensities, and the scattering phase function at 5° suggest that scattering was due to aerosols of very small sizes. The estimates of the weighted aerosol mean square radius,p ², seem to lead to the same conclusion.     

High frequency S wave envelope synthesis using a multiple nonisotropic scattering model: application to aftershocks from the 2008 Wenchuan earthquake

Based on the formulation of a multiple non-isotropic scattering process, a characteristic source time is introduced to define the initial impulse width of energy density at the source. An analytical expression of the initial intensity spectral density of a seismic wave is incorporated into the integral equation of seismic wave energy density. And, a recursive formula of Green’s function is derived to obtain the higher order Green’s function, which is included to describe the stronger non-isotropic scattering process. Then, the effect of the scattering pattern on the energy density envelope is investigated by the modified scattering theory. Significant differences are found in the decay of the energy density envelopes with distances using different scattering patterns. The envelope synthesized by the forward dominated scattering pattern is larger than the results obtained by the isotropic and backward dominated scattering pattern. Different scattering patterns are also used to fit the observation data from the aftershocks of the 2008 Wenchuan earthquake. It is concluded that the envelopes synthesized by the forward scattering pattern can match the data better than the isotropic and backward dominated scattering cases, and a new interpretation of the coda wave is given. Finally, using the forward dominated scattering pattern, the envelope broadening of the observed data is reproduced.     

Modeling Full Seismogram Envelopes Using Radiative Transfer Theory with Born Scattering Coefficients

The equation of radiative transfer is used to model the transport of seismic energy in 2-D and 3-D acoustic random media. Monte-Carlo solutions of this equation using non-isotropic Born scattering coefficients are compared to three analytical solutions: Markov approximation, radiative transfer theory with isotropic scattering coefficients, and diffusion approximation. Additionally, we compare to finite differences solutions of the full wave equation in 2-D. We find a good correspondence of radiative transfer theory to Markov approximation for the case of multiple forward scattering. The comparison to radiative transfer theory with isotropic scattering coefficients, a model frequently used in data analysis, demonstrates that in the case of forward scattering the isotropic scattering model is not better than a diffusion approach. To compare radiative transfer theory with non-isotropic scattering coefficients to finite differences solutions of the full wave equation, the finite source duration and the bandpass filter process as well as the normalization of absolute amplitudes are explicitely taken into account. We find a good coincidence of both theories for scattering parameters, which are realistic for usual Earth crust. The theory correctly describes the unscattered direct wavefront, the envelope broadening caused by multiple forward scattering, as well as the late coda caused by multiple wide angle scattering. For strong scattering, which can be expected for very heterogeneous media such as strato volcanoes, the solutions of radiative transfer differ from the more complete solutions of the full wave equation.