Characteristics of source parameters of the Shanxi Reservoir earthquakes in Zhejiang region
-
摘要: 利用2014年1月-2015年12月浙江省区域数字地震台网记录的1 360次珊溪水库ML0.5-4.4地震资料，采用Brune模式，将台站获取的速度记录进行几何扩散校正、介质衰减校正及仪器校正后得到速度谱，对速度谱在频率域进行傅里叶积分得到震源位移谱，之后使用遗传算法计算得到拐角频率和震源位移谱零频极限值，进而计算出该区的地震矩为1010-1014 N·m，震源破裂半径为33-550 m，拐角频率为2.4-39.7 Hz，地震应力降为0.04-6.74 MPa，视应力为0-2.75 MPa.在此基础上，对各参数特征及其关系进行系统性分析，结果表明：各震源参数之间表现出一定的对数线性或半对数线性关系；空间上，应力降和拐角频率在库区西北段较高，东南段较低；应力降与地震矩在空间上呈反向关系，拐角频率与破裂半径在空间上呈负相关.Abstract: According to the seismic data of 1 360 earthquakes with ML0.5-4.4 occurred in Shanxi Reservoir from January of 2014 to December of 2015, which were recorded by the Regional Digital Seismic Network of Zhejiang Province, this paper obtains the source displacement spectra by using Fourier integration of velocity spectrum in the frequency field based on the Brune's model, after the geometrical diffusion correction, medium attenuation correction and instrument correction of the velocity records from the stations, and then the corner frequency and the source displacement spectra at the zero frequency are acquired with the genetic algorithm. And then the radius of the epicenter, the seismic moments and other small local earthquake source parameters are calculated. Furthermore, this paper systematically analyzes the characteristics of the source parameters and the relationship among the parameters. The results show that: seismic moment ranges 1010-1014 N·m, the range of source radius is 33-550 m, the corner frequency ranges 2.4-39.7 Hz, seismic stress drop is in the range of 0.04-6.74 MPa, and apparent stress is in the range of 0-2.75 MPa. The result also show that there is some log-linear or semi-logarithmic linear relationships among the seismic moment, corner frequency, Richter magnitude, moment magnitude and other source parameters. In the space domain, the stress drop and the corner frequency are higher at northwest section of Shanxi Reservoir area and lower at southeast section. In addition, there is opposite relationship between the stress drop and the seismic moment in space, meanwhile the corner frequency is negatively correlated to the source radius in space.
-
引言
库区地震是浙江省地区比较显著的地震事件，近年来，珊溪水库库区发生多次地震，引起了各方的关注.珊溪水库库区的断裂构造较为复杂，库区范围内约有14条断裂，断裂一般长达10余千米，断裂带宽约20-30 m，位于5 km深度以上.出露的断裂主要有NE向和NW向两组：NE向断裂的走向为NE40°-60°，倾向以NW为主，倾角为60°-80°，多为逆走滑断层；NW向断裂的走向为NW40°-50°，倾向以NE为主，倾角为60°-70°，多为逆断层和逆走滑断层，未发现NW向断裂切错水系的现象.
珊溪水库于2000年5月开始蓄水，2002年7月库区开始发生地震，截至目前一直有强弱不等的地震发生，其中最显著的高密度震群活动有3组：2002年7-10月为第一组，最大为ML3.7地震；2006年2-11月为第二组，最大为ML4.6，期间2-5月发生密集的震群活动，之后出现长达8年的ML≥3.0地震平静；2014年9月-2015年4月为第三组，期间2014年9-12月发生密集的震群活动，对比朱新运等(2010)定位第一、二组地震震中位置结果可知，本组地震填充了之前该库区地震活动的弱活动段，共发生ML≥2.0地震258次，最大地震为ML4.4，震中为(119.949°E，27.704°N)，震源深度为3 km.
浙江省地震局在珊溪水库地震发生期间建成了子台密度较大的数字地震台网，包括新浦、联云、包垟、云湖、泰顺、珊溪、黄坛等台站；这些台站近距离记录了大量地震资料，为计算震源参数以及分析各参数之间的拟合关系提供了详实的数据支持.针对本区域水库地震的前期地震已进行了大量研究工作，包括拐角频率对地震序列的依赖性研究(朱新运等，2004)、库区介质衰减研究(张帆等，2013)、震源参数研究(钟羽云等，2004；于俊谊，朱新运，2008)等.这些研究为该区域地震的发生发展提供了合理完整的理论解释. 2014年9月震群再次活动，其强度和频度基本接近于该区域前期的地震活动水平.本文拟利用这些地震的数字地震波资料研究其震源参数，对该区域前期地震的孕育过程进行必要补充，以期结合震源参数空间分布规律和震例，深刻认识水库区域地震的发生规律.
1. 原理与方法
台站记录中包含了地震震源、地震波的传播途径及场地响应等信息，通过有效简化或扣除台站记录中的非震源信息，基于震源谱理论模型获得震源参数，这是震源参数提取的有效手段(Aaron et al，2000；Bindi et al，2001；Giampiccolo et al，2007).
台站记录到的观测谱可以表示为
(1) 式中，Sobs(f)为台站观测谱，Oobs(f)为地震震源谱，G(f)为场地响应，I(f)为仪器响应，P(f)=R0-1exp[-πR0f/(Q(f)vS)]为路径效应，R0为震源距，f为频率，Q(f)为衰减系数，vS为横波速度.
取EW向和NS向均具有较高信噪比的S波记录数据计算傅里叶谱.计算中使用体波数据，几何衰减使用r-1校正，设定S波的非弹性、非均匀性衰减与其本身的尾波衰减相一致(Rautian，Khalturin，1978；Aki，1980；Herrmann，1980；Matsunmni，Nakamura，2004)，介质非均匀性衰减采用朱新运和张帆(2006)通过尾波衰减所获取的参数进行校正，场地效应采用朱新运和陈运泰(2007)利用反演获得的结果进行校正.
为了保证地震波谱数据的稳定性，在计算傅里叶谱时采用Chael (1987)提出的延时窗技术，即将S波段分割为包含n个点的小数据段计算傅里叶谱，之后以1/2数据段前移，前次计算数据段的一半被重叠，对S波重新分割后的n个小数据段进行第二次计算，依次类推，这样就获得了整个数据段分割、重叠的n段谱数据，用下式将n小段数据归算为全部S波段的谱数据：
(2) 式中：O(f)为观测位移谱；n为数据段总数，这里取n=256；T为整段S波窗长；t为数据分段窗长.使用速度记录进行谱计算时乘以系数1/(2πf)将其转换为位移谱.使用式(2)分别计算EW分量和NS分量的震源谱，并根据下式合成最终震源观测谱(Aaron et al，2000)：
(3) 式中，OEW(f)和ONS(f)分别为EW向和NS向的震源谱，Oobs(f)为单个台站对一个地震记录的震源谱.计算S波震源谱结果，舍弃低信噪比的地震波记录，将多台对同一地震记录的振幅谱均值作为该地震的有效振幅谱.
根据Brune (1970)的ω2模型(Shi et al，1998)，理论震源位移谱可表示为
(4) 式中，Otheo(f)为理论震源位移谱，Ω0为震源谱低频极限值，f0为拐角频率.以Ω0和f0为独立变量，利用遗传算法或最小二乘法，使Oobs(f)与Otheo(f)具有最小残差，从而确定参数Ω0和f0.
本文采用遗传算法，确定下式为目标函数(朱新运，于俊谊，2008).根据仪器性能和采样率，在1-18 Hz时，仪器响应的振幅稳定且相位较低，所以限定在1-18 Hz范围内搜索频率，使目标函数最小，进而获得Ω0和f0.
(5) 式中，m为频率点数，ε为目标值，Oobs为同一地震的波谱转化为震源位移谱后的多台平均值，Otheo为根据ω2模型计算所得的震源谱，k为频率点序号.
以2014年10月25日库区发生的ML4.4地震为例，计算泰顺台的单台位移谱、多台站位移谱平均值以及拟合得到的理论震源谱，如图 1所示.
![]()
图 1 2014年10月25日珊溪水库ML4.4地震的位移谱计算过程实例(a)单台(泰顺台)震源谱与噪声谱比较；(b)多台站震源谱及其平均值和多台数据拟合结果Figure 1. An example for calculating the seismic displacement spectra of the ML4.4 earthquake occurred in Shanxi Reservoir on 25 October 2014(a) Comparison of the source spectrum recorded by Taishun station with noise spectrum. Blue line represents north-south and east-west observation spectra in NS and EW direction; red line represents the average observation spectra, black dotted line indicates the noise spectrum. (b) Source spectra (blue lines) and their average spectrum (green line) as well as data fitting of spectra of a series of stations by using the ω2 model (red line)确定低频极限值Ω0和拐角频率f0后，即可由下式获取地震矩M0，破裂半径r0，应力降Δσ和视应力σapp(Brune，1970)：
(6) 式中，RS(φ，θ)为辐射方向性因子，Re为自由表面反射系数. 2014年10月25日珊溪水库ML4.4地震震源参数的计算结果为：拐角频率为3.5，地震矩为5.73×1014 N·m，震源破裂半径为386 m，应力降为4.34 MPa，矩震级为3.8，视应力1.77 MPa.
2. 数据与计算结果
本文所用资料来源于55个台站记录到的2014年1月-2015年12月珊溪水库库区(119.758-120.393°E，27.452-29.218°N)发生的1 360次ML0.5-4.4地震，其中1.0≤ML≤1.9地震846次，2.0≤ML≤2.9地震207次，3.0≤ML≤3.9地震43次，4.0≤ML≤4.9地震8次；最大地震为ML4.4，震中为(119.949°E，27.704°N)，震源深度为3.0 km.所用仪器为BBVS-60型、FBS-3B型宽频带地震计和FSS-3B短周期地震计，仪器对频率响应的平坦段在1 Hz以下，采样频率为100 Hz的仪器对频率响应的高段可达40 Hz，系统动态范围优于120 dB.震中与台站分布如图 2所示.
![]()
图 2 珊溪水库库区地震震中与台站的分布，图(b)为图(a)中方框的放大图Figure 2. Distribution of earthquake epicenters (dots) and stations (triangles) in Shanxi Reservoir Fig.(b) is amplification of the square in Fig.(a)计算中取
3. 结果分析
研究结果表明地震矩M0，震源破裂半径r0，应力降Δσ，MW以及ML之间呈一定的相关关系，这种地震震源参数之间的关系被称为定标关系.
3.1 ML与MW及ML与M0的定标关系
震级是用来衡量地震相对大小的物理量，本文基于地震事件和震级数据，使用最小二乘法拟合出ML与MW的关系为
(7) 地震矩M0是对断层滑动引起的地震强度的直接测量，反映了震源处破裂的大小，所计算地震的震级范围为ML0.5-4.4，相对应的地震矩范围为1010-1014 N·m.作为另一个描述地震绝对大小的力学量，M0与ML的关系如图 3所示，可以看出M0的对数与ML之间呈很好的线性关系，其表达式为
(8) ![]()
图 3 (a) MW与ML拟合的定标关系；(b) M0与ML拟合的定标关系Figure 3. (a) Scaling relation between MW and ML; (b) Scaling relation between M0 and ML结合式(7)和式(8)，可得M0与MW的关系为
(9) 这与Kanamori (1977)给出的M0与MW之间的关系
(10) 极为相近.
3.2 M0与r0及M0与f0的定标关系
库区内的震源破裂半径r0=33-550 m，其中r0=33-75 m范围内的地震有744次，r0=75-135 m范围内的有401次，大部分地震的破裂半径尺度较小. 图 4a给出了破裂半径的分布直方图，可以看出，地震次数N随着r0的增大而迅速减小；图 4b给出了M0与r0的定标关系，总体上看M0越大，r0越大，最小二乘拟合后M0与r0的定标关系为
(11) ![]()
图 4 破裂半径r0，拐角频率f0的分布直方图以及地震矩M0与r0和f0的定标关系(a) r0分布直方图；(b) M0与r0拟合的定标关系；(c) f0分布直方图；(d) f0与M0拟合的定标关系Figure 4. Histograms of source radius r0 and corner frequency f0 as well as scaling relation between M0, r0 and f0(a) Histogram of r0; (b) Scaling relation between M0 and r0; (c) Histogram of f0; (d) Scaling relation between f0 and M0图 4c给出了拐角频率f0的分布直方图，可以看出，f0=2.4-39.7 Hz，其中：f0=2.4-10 Hz范围内的地震有221次；f0=10-25 Hz范围内有881次，占总数的64.8%；f0=25-40 Hz的有258次.一般而言，震级愈大，地震震源破裂半径愈大，拐角频率越小.本文对拐角频率和震源破裂半径进行统计的结果显示，破裂半径小的地震次数多，这与序列中大小地震次数关系呈正相关，但拐角频率却出现了有意义的变化.通常地震震级越小，拐角频率愈大，按照地震频度G-R关系，拐角频率与频度的关系应该与震源破裂半径与频度的关系相反，但本文所计算的拐角频率却表现为类正态分布(图 4c)，主要集中在10-25 Hz之间.按照下式(Brune, 1970, 1971)
(12) 可知，在横波速度稳定的情况下，f0与r成反比关系，这说明区域内vS在地震时有过较大扰动.由于地震时vS降低，使得大部分地震的f0并未集中在其高段. 图 4d给出了地震矩M0与拐角频率f0的关系，可以看出，M0越大，f0越小，其拟合关系式为
(13) 3.3 Δσ与M0的定标关系
图 5a给出了Δσ与M0的定标关系，图 5b给出了珊溪水库库区Δσ的分布直方图，可以看出，库区内Δσ=0.04-6.74 MPa，其中：Δσ=0-0.5 MPa的地震有802次；Δσ=0.5-1.0 MPa的有299次，约占地震总数的81%；其余258次分布在Δσ=1.0-7.0 MPa范围内，整组地震的应力降均值约为0.87 MPa. 图 5c给出了珊溪水库库区视应力σapp的分布直方图，可以看出，σapp=0.000 2-2.75 MPa，其中：σapp=0-0.5 MPa的地震有1 184次；σapp=0.5-1.0 MPa的有124次，约占地震总数的95%；其余分布在1.0-2.75 MPa之间；整组地震的视应力均值约为0.37 MPa.水库库区低Δσ和低σapp说明震群属于低应力背景下的地震活动，其f0较低，这与陈培善和Duda (1993)的研究结果“若地震发生在低背景应力条件下，地震震源脉冲低而宽，相应的高频成分少，拐角频率小”相一致.
![]()
图 5 应力降Δσ与地震矩M0的定标关系及Δσ和σapp分布直方图(a) Δσ与M0拟合的定标关系；(b) Δσ分布直方图；(c) σapp分布直方图Figure 5. Distribution of stress drop Δσ, seismic moment M0 and apparent stress σapp(a) Scaling relation between Δσ and M0; (b) Histogram of Δσ; (c) Histogram of σapp水库区地震的应力降的对数与地震矩的对数具有较显著的线性正相关关系，拟合得到的定标关系式为
(14) 此外，珊溪水库库区重复发生的小地震，其应力降Δσ随时间的推移呈系统性的变化.如图 6所示，Δσ随时间的变化与M0随时间变化呈现出较为一致的变化，这与陈运泰等(2000)持有的“对于不同地震矩的小地震，应力降的大小随地震的增大而增大”的观点相一致.
![]()
图 6 应力降Δσ(a)和地震矩M0 (b)随时间的变化Figure 6. Variation of stress drop Δσ (a) and seismic moment M0 (b) with time3.4 Δσ与r0及Δσ与f0之间的关系
如图 7所示，Δσ与r0及Δσ与f0的对数拟合后均无明显的线性关系.
![]()
图 7 应力降Δσ与破裂半径r0 (a)和拐角频率f0 (b)的关系Figure 7. Distribution of data fitting between stress drop Δσ and source radius r0 (a) and Δσ and corner frequency f0 (b)3.5 震源参数空间分布特征
将2014-2015年共计1 360次地震的震源参数按照400 m×400 m进行网格化，绘出震源参数的空间分布，如图 8和图 9所示.由图 8a可以看出：研究区域应力降空间分布的范围为0.033-2.61 MPa，大部分区域应力降集中在0.5-1.0 MPa范围内变化，应力降最大区域位于库区东南部，并且其附近应力降均高于均值；应力降沿双溪-焦溪垟断裂分布不均匀，总体上靠近断裂且位于水库淹没区西北段的应力降较小，靠近断裂且位于水库淹没区东南段的较大.另一方面，应力降的变化与震中位置有关，在研究区域西北段发生地震较多，震级较大，该段地区应力降处于均值范围内；东南段发生地震较少，应力降较大.例如，珊溪水库2014年9-12月地震主要发生在双溪-焦溪垟断裂带淹没区西北段.如图 8b所示，水库库区的地震矩范围为3×1010-2.5×1013 N·m，库区西北段的地震矩数值较大，这与该区多次发生较大地震有关，而在库区东南段地震矩数值较小，大部分处于1×1010-5×1012 N·m之间.
![]()
图 8 珊溪水库库区应力降Δσ (a)和地震矩M0 (b)空间分布图Figure 8. Distribution of stress drop Δσ (a) and seismic moment M0 (b) in the Shanxi Reservoir region在空间上，珊溪水库库区应力降与地震矩存在负相关关系，即地震矩小的区域(库区东南段)应力降大，地震矩大的区域(库区西北段)应力降小；应力降与地震矩在空间上的这种负相关关系，可能与该地区的地震活动情况有关.库区西北段发生地震、释放应力后，应力水平较低，相应的应力降较低，而地震活动较低的西南段应力水平较高，应力降也相对较高.应力降与地震矩在数值上表现为正相关，而在空间上呈负相关关系，因此应力降较高的库区东南段也是未来可能发生较大地震的重点监测区域.
如图 9所示，震源区的地震拐角频率f0为2.4-39.7 Hz，表现为西北段较大，东南段较小，基本呈北低南高的特点；破裂半径r0为33-550 m.可以看出，f0与r0在空间上呈明显的负相关.
![]()
图 9 珊溪水库库区拐角频率f0(a)和破裂半径r0(b)的空间分布Figure 9. Distribution of corner frequency f0 (a) and source radius r0 (b) in the region of Shanxi Reservoir4. 讨论与结论
本文基于Brune (1970)提出的小震震源模型, 使用经验公式计算浙江省珊溪水库1 360次ML0.5-4.4地震的震源参数.由于这些经验公式本身具有一定的偏差，计算处理单次地震时利用多个近台结果平均得到该地震的震源参数，这样使得单次地震的震源参数误差更小，此外还通过统计大量的中、小地震震源参数减小计算结果误差所带来的影响.在此基础上，进一步研究各震源参数之间的标定关系，得到的主要结论如下：
1)珊溪水库库区地震的地震矩对数与破裂半径对数呈较显著的线性正相关，而与拐角频率的对数呈显著的线性负相关；地震应力降与地震矩表现出正相关关系，而与震源破裂半径未呈明显的相关.此外，应力降与地震矩在时间上呈正相关.
2) 于俊谊和朱新运(2008)使用本区域第一、二组地震中ML≥2.0地震资料进行震源参数研究.本研究与其标定关系(表 1)相差不大，这说明在同一区域不同时间尺度上、不同震级范围的地震震源参数不随地震序列的发展而变化.
表 1 本文与于俊谊和朱新运(2008)关于珊溪水库库区震源参数拟合定标关系结果对比Table 1. Comparison of fitting scaling relations of source parameters in the Shanxi Reservoir between this paper and Yu and Zhu (2008)关系式 k1 k2 本文 于俊谊和朱新运(2008) 本文 于俊谊和朱新运(2008) MW=k1ML+k2 0.708 0.74 0.732 0.70 lgM0=k1ML+k2 1.052 1.11 10.16 10.11 lgM0=k1lgr0+k2* 2.195 7.471 lgM0=k1r0* 0.0061 11.71 lgf0=k1lgM0+k2 -0.221 -0.26 3.779 3.97 lgΔσ=k1lgM0+k2 0.332 0.281 1.667 2.158 *关于M0与r0的定标关系，于俊谊和朱新运(2008)采用M0对数与r0拟合，本文采用M0对数与r0对数拟合；若采用M0对数与r0拟合，本文所得结果为k1=0.007 1和k2=10.967，与于俊谊和朱新运(2008)结果近似. 3)珊溪水库库区的地震震源破裂半径r0=33-550 m，大部分地震的r0较小，地震次数随着r0的增大而迅速减小；拐角频率f0=2.4-39.7 Hz，地震次数与f0并未呈现出正相关关系，而是大部分处于拐角频率中段10-25 Hz；f0与r0未成反比关系，说明区域内地震波横波速度在地震时有较大扰动.
4)珊溪水库库区的应力降Δσ=0.04-6.74 MPa，均值为0.87 MPa，小震的应力降集中分布在0-0.5 MPa；视应力σapp=0.000 2-2.75 MPa, 均值约为0.37 MPa，说明该震群属于低应力背景下的地震活动.
5)应力降和拐角频率在水库库区东南段较高，西北段较低，而地震矩和破裂半径在水库库区东南段较低，西北段较高；应力降与地震矩在空间上呈反向关系，这与二者在数值大小上的关系并不一致；拐角频率和破裂半径在空间上呈负相关.从地震分布的角度看，本组地震填补了朱新运等(2010)一文中所提到的空段，而本组地震应力释放较大的地方也正是第一、二组地震应力释放不充分的区段.
-
![]()
图 1 2014年10月25日珊溪水库ML4.4地震的位移谱计算过程实例
(a)单台(泰顺台)震源谱与噪声谱比较；(b)多台站震源谱及其平均值和多台数据拟合结果
Figure 1. An example for calculating the seismic displacement spectra of the ML4.4 earthquake occurred in Shanxi Reservoir on 25 October 2014
(a) Comparison of the source spectrum recorded by Taishun station with noise spectrum. Blue line represents north-south and east-west observation spectra in NS and EW direction; red line represents the average observation spectra, black dotted line indicates the noise spectrum. (b) Source spectra (blue lines) and their average spectrum (green line) as well as data fitting of spectra of a series of stations by using the ω2 model (red line)
![]()
图 2 珊溪水库库区地震震中与台站的分布，图(b)为图(a)中方框的放大图
Figure 2. Distribution of earthquake epicenters (dots) and stations (triangles) in Shanxi Reservoir Fig.(b) is amplification of the square in Fig.(a)
![]()
图 3 (a) MW与ML拟合的定标关系；(b) M0与ML拟合的定标关系
Figure 3. (a) Scaling relation between MW and ML; (b) Scaling relation between M0 and ML
![]()
图 4 破裂半径r0，拐角频率f0的分布直方图以及地震矩M0与r0和f0的定标关系
(a) r0分布直方图；(b) M0与r0拟合的定标关系；(c) f0分布直方图；(d) f0与M0拟合的定标关系
Figure 4. Histograms of source radius r0 and corner frequency f0 as well as scaling relation between M0, r0 and f0
(a) Histogram of r0; (b) Scaling relation between M0 and r0; (c) Histogram of f0; (d) Scaling relation between f0 and M0
![]()
图 5 应力降Δσ与地震矩M0的定标关系及Δσ和σapp分布直方图
(a) Δσ与M0拟合的定标关系；(b) Δσ分布直方图；(c) σapp分布直方图
Figure 5. Distribution of stress drop Δσ, seismic moment M0 and apparent stress σapp
(a) Scaling relation between Δσ and M0; (b) Histogram of Δσ; (c) Histogram of σapp
![]()
图 6 应力降Δσ(a)和地震矩M0 (b)随时间的变化
Figure 6. Variation of stress drop Δσ (a) and seismic moment M0 (b) with time
![]()
图 7 应力降Δσ与破裂半径r0 (a)和拐角频率f0 (b)的关系
Figure 7. Distribution of data fitting between stress drop Δσ and source radius r0 (a) and Δσ and corner frequency f0 (b)
![]()
图 8 珊溪水库库区应力降Δσ (a)和地震矩M0 (b)空间分布图
Figure 8. Distribution of stress drop Δσ (a) and seismic moment M0 (b) in the Shanxi Reservoir region
![]()
图 9 珊溪水库库区拐角频率f0(a)和破裂半径r0(b)的空间分布
Figure 9. Distribution of corner frequency f0 (a) and source radius r0 (b) in the region of Shanxi Reservoir
表 1 本文与于俊谊和朱新运(2008)关于珊溪水库库区震源参数拟合定标关系结果对比
Table 1 Comparison of fitting scaling relations of source parameters in the Shanxi Reservoir between this paper and Yu and Zhu (2008)
关系式 k1 k2 本文 于俊谊和朱新运(2008) 本文 于俊谊和朱新运(2008) MW=k1ML+k2 0.708 0.74 0.732 0.70 lgM0=k1ML+k2 1.052 1.11 10.16 10.11 lgM0=k1lgr0+k2* 2.195 7.471 lgM0=k1r0* 0.0061 11.71 lgf0=k1lgM0+k2 -0.221 -0.26 3.779 3.97 lgΔσ=k1lgM0+k2 0.332 0.281 1.667 2.158 *关于M0与r0的定标关系，于俊谊和朱新运(2008)采用M0对数与r0拟合，本文采用M0对数与r0对数拟合；若采用M0对数与r0拟合，本文所得结果为k1=0.007 1和k2=10.967，与于俊谊和朱新运(2008)结果近似. 
下载: 导出CSV
-
陈培善, Duda S J. 1993.峰值速度和加速度对环境剪应力的依赖性[J].地震学报, 36 (2): 185-194. http://www.dzxb.org/Magazine/Show?id=28525 Chen P S, Duda S J. 1993. Dependence of peak acceleration and peak velocity on the environment shear stress[J]. Acta Seismologica Sinica, 36 (2): 185-194. http://www.dzxb.org/Magazine/Show?id=28525
陈运泰, 吴忠良, 王培德, 许力生, 李鸿吉, 牟其铎. 2000.数字地震学[M].北京:地震出版社: 1-172. Chen Y T, Wu Z L, Wang P D, Xu L S, Li H J, Mu Q D. 2000. Digital Seismology[M]. Beijing: Seismological Press: 1-172 (in Chinese).
于俊谊, 朱新运. 2008.浙江珊溪水库地震震源参数研究[J].中国地震, 24 (4): 379-387. http://www.cnki.com.cn/Article/CJFDTOTAL-ZGZD200804008.htm Yu J Y, Zhu X Y. 2008. Hypocentral parameters research of the Shanxi reservoir area in Zhejiang[J]. Earthquake Research in China, 24 (4): 379-387 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-ZGZD200804008.htm
张帆, 朱新运, 钟羽云. 2013.珊溪水库地震小震震源机制解特征研究[J].地震工程学报, 35 (3): 684-691. http://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ201303045.htm Zhang F, Zhu X Y, Zhong Y Y. 2013. Study on the characteristic of focal mechanism solutions of the Shanxi reservoir earthquake[J]. China Earthquake Engineering Journal, 35 (3): 684-691 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ201303045.htm
钟羽云, 朱新运, 张震峰. 2004.温州珊溪水库ML3.9震群震源参数特征[J].地震, 24 (3): 107-114. http://www.cnki.com.cn/Article/CJFDTotal-DIZN200403015.htm Zhong Y Y, Zhu X Y, Zhang Z F. 2004. Source parameter characteristics of ML3.9 earthquake swarm occurred in Shanxi reservoir, Wenzhou, Zhejiang Province[J]. Earthquake, 24 (3): 107-114 (in Chinese). http://www.cnki.com.cn/Article/CJFDTotal-DIZN200403015.htm
朱新运, 钟羽云, 张震峰. 2004.地震波拐角频率对地震序列的依赖性研究[J].西北地震学报, 26 (2): 131-136. http://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ200402006.htm Zhu X Y, Zhong Y Y, Zhang Z F. 2004. Study on seismic wave corner frequency depending on earthquake sequences[J]. Northwestern Seismological Journal, 26 (2): 131-136 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ200402006.htm
朱新运, 张帆. 2006.浙江及邻区尾波Q值研究[J].西北地震学报, 28 (2): 108-113. http://www.cnki.com.cn/Article/CJFDTotal-ZBDZ200602002.htm Zhu X Y, Zhang F. 2006. Research on coda Q value in Zhejiang Province and neighbouring areas[J]. Northwestern Seismological Journal, 28 (2): 108-113 (in Chinese). http://www.cnki.com.cn/Article/CJFDTotal-ZBDZ200602002.htm
朱新运, 陈运泰. 2007.用Lg波资料反演场地效应与地震波衰减参数[J].地震学报, 29 (6): 569-580. http://www.dzxb.org/Magazine/Show?id=26451 Zhu X Y, Chen Y T. 2007. An inversion of site response and Lg attenuation using Lg waveform[J]. Acta Seismologica Sinica, 29 (6): 569-580 (in Chinese). http://www.dzxb.org/Magazine/Show?id=26451
朱新运, 于俊谊. 2008.基于MATLAB的小震震源参数计算软件研制[J].西北地震学报, 30 (4): 380-384. http://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ200804012.htm Zhu X Y, Yu J Y. 2008. The data processing software of source parameters for middle or small earthquakes based on MATLAB[J]. Northwestern Seismological Journal, 30 (4): 380-384 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ200804012.htm
朱新运, 张帆, 于俊谊. 2010.浙江珊溪水库地震精细定位及构造研究[J].中国地震, 26 (4): 380-390. http://www.cnki.com.cn/Article/CJFDTOTAL-ZGZD201004005.htm Zhu X Y, Zhang F, Yu J Y. 2010. Study on precise positioning and structure of the Shanxi reservoir earthquake, Zhejiang Province[J]. Earthquake Research in China, 26 (4): 380-390 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-ZGZD201004005.htm
Aaron M, Jorge A, Kojiro I. 2000. Inversion of source parameters and site effects from strong ground motion records using genetic algorithms[J]. Bull Seismol Soc Am, 90 (4): 977-992. doi: 10.1785/0119990007
Aki K. 1980. Attenuation of shear-waves in the lithosphere for frequencies from 0.05 to 25 Hz[J]. Phys Earth Planet Inter, 21 (1): 50-60. doi: 10.1016/0031-9201(80)90019-9
Bindi D, Spallarossa D, Augliera P, Cattaneo M. 2001. Source parameters estimated from the aftershocks of the 1997 Umbria Marche (Italy) seismic sequence[J]. Bull Seismol Soc Am, 91 (3): 448-455. doi: 10.1785/0120000084
Brune J N. 1970. Tectonic stress and the spectra of seismic shear waves from earthquakes[J]. J Geophys Res, 75 (26): 4997-5009. doi: 10.1029/JB075i026p04997
Brune J N. 1971. Correction[J]. J Geophys Res, 76 (20): 5002. doi: 10.1029/JB076i020p05002
Chael E P. 1987. Spectral scaling of earthquakes in the Miramichi region of New Brunswick[J]. Bull Seismol Soc Am, 77 (2): 347-365.
Giampiccolo E, D′Amico S, Patanè D, Gresta S. 2007. Attenuation and source parameters of shallow microearthquakes at Mt. Etna Volcano, Italy[J]. Bull Seismol Soc Am, 97 (1B): 184-197. doi: 10.1785/0120050252
Herrmann R B. 1980. Q estimates using the coda of local earthquakes[J]. Bull Seismol Soc Am, 70 (2): 447-468. http://cn.bing.com/academic/profile?id=2054846959&encoded=0&v=paper_preview&mkt=zh-cn
Kanamori H. 1977. The energy release in great earthquakes[J]. J Geophys Res, 82 : 2981-2987. doi: 10.1029/JB082i020p02981
Matsunmni K, Nakamura M. 2004. Seismic attenuation in a nonvolcanic swarm region beneath Wakayama, southwest Japan[J]. J Geophys Res, 109 (B9): B09302. http://cn.bing.com/academic/profile?id=2063020719&encoded=0&v=paper_preview&mkt=zh-cn
Rautian T G, Khalturin V I. 1978. The use of the coda for determination of earthquake source spectrum[J]. Bull Seismol Soc Am, 68 (4): 923-948.
Shi J H, Kim W Y, Paul G R. 1998. The corner frequencies and stress drops of intraplate earthquakes in the northeastern United States[J]. Bull Seismol Soc Am, 88 (2): 531-542. http://cn.bing.com/academic/profile?id=2277660255&encoded=0&v=paper_preview&mkt=zh-cn
Zhu X Y. 2014. An inversion of Lg-wave attenuation and site response from seismic spectral ratios in the eastern China region[J]. Bull Seismol Soc Am, 104 (3): 1389-1399. doi: 10.1785/0120120359.
-
期刊类型引用(6)
1. 汪建，朱新运，魏红梅，汤兰荣，杨亚运，向涯，朱祥. 重庆区域地震震源参数特征研究. 四川地震. 2024(03): 8-15+23 . 
百度学术
2. 龚俊，汪贞杰，戴煜暄. 联合测震和GPS资料反演双溪-焦溪垟断裂运动特征. 测绘科学. 2021(12): 1-6+15 . 百度学术
3. 郁建芳，谢石文，张炳，隆爱军，韩成成，杨波. 安徽地区震源参数的相关性研究. 华南地震. 2019(03): 32-41 . 百度学术
4. 汪贞杰，孙侃，朱新运. 单台sPL震相测定珊溪水库地震震源深度. 地震学报. 2019(06): 735-742+808 . 本站查看
5. 熊凡. 遗传算法优化支持向量机的网络流量混沌预测. 现代电子技术. 2018(18): 166-169 . 百度学术
6. 张锦玲，朱新运，马起杨，金春华. 宁夏区域震源参数特征研究. 大地测量与地球动力学. 2017(08): 808-812 . 百度学术
其他类型引用(1)
计量
- 文章访问数: 665
- HTML全文浏览量: 277
- PDF下载量: 39
- 被引次数: 7
