Determination of the surface-wave magnitude of the Wenchuan earthquake and its seismic Doppler effect
-
摘要: 本文介绍了全球主要地震机构对2008年5月12日汶川地震参数的速报和修订情况，分析了美国地质调查局国家地震信息中心测定的面波震级。通过对比198个全球地震台站测定的面波震级和面波周期，得出如下结论：测定面波震级偏大的台站主要分布在震中的东北方向，测定面波震级偏小的台站主要分布在震中的西南和东南方向，面波周期偏小的台站主要分布在震中东北方向。由于此次地震破裂方向是以北东向单侧破裂为主，且地震多普勒效应导致震中东北方向振动加强，因此该方向上的面波震级偏大，地震烈度衰减慢；而震中西南方向的振动减弱，此方向面波震级偏小，地震烈度衰减快。从而造成地震烈度沿中央断裂带的北东方向衰减慢，而南西方向衰减快的特征分布。Abstract: The fast report and revision of the parameters of the Wenchuan earthquake on 12 May 2008 by some major international seismological institutions are introduced, and the surface-wave magnitude measured by National Earthquake Information Center of United States Geological Survey is analyzed in this paper. The comparison of the measured surface-wave magnitudes and periods of 198 global seismic stations reveals that, those stations with larger magnitudes are mainly in NE direction of the epicenter while those with smaller magnitudes are in SW and SE direction, and stations with smaller periods are mainly in NE direction. The seismic intensity map of the Wenchuan earthquake shows that the intensity attenuated slowly in the NE direction of the central fracture zone while the intensity attenuated fast in the SW direction. The reason is that Wenchuan earthquake is a unilateral rupture mainly in NE direction, due to the seismic Doppler effect, the stronger shaking leads to the larger surface-wave magnitudes and slower intensity attenuation in the NE direction of the epicenter, vice versa the weaker shaking leads to the smaller surface-wave magnitudes and faster intensity attenuation in the SW direction.
-
Keywords:
- Wenchuan earthquake /
- surface-wave magnitude /
- Doppler effect /
- seismic intensity
-
引言
地震Lg波对地壳的横向变化非常敏感，在稳定的地壳区域，如克拉通和地盾等，Lg波可以传播至3 000 km以外。但当Lg波穿过大陆地壳较薄区、地壳厚度突变区、巨厚的沉积盆地或存在部分熔融物质的地区时，Lg波会呈现出明显的衰减差异（Nuttli，1980；Kennett，1986；Mitchell，Hwang，1987；Bostock，Kennett，1990；Mitchell et al，1997；Rodgers et al，1997；Schurr et al，2003；Xie et al，2004）.例如在青藏高原北部的强烈衰减区，Lg波的传播不超过700 km （McNamara et al，1996；Xie，2002）。地震波振幅衰减的多少可以用地下介质品质因子Q值来衡量，它是了解地下裂隙的数量、孔隙密度与分布，孔隙中流体含量以及地壳温度的重要参数。例如青藏高原的Lg波Q值极低（约为100—200），Chen和Xie （2017）认为这可能是青藏高原地壳内的高温造成的。由于构造因素导致的Q值变化比速度更加敏感，因此借助Lg波对Q值进行研究可以为了解地下介质提供更多的约束和解释（Akinci et al，1995；Rietbrock，2001；Gung，Romanowicz，2004）。
Q值层析成像是一个非常复杂的问题，因为Q值反演所使用的地震波振幅数据不仅受有限的Q值的影响，而且还受三维速度结构中波前扩展所导致的未知的三维几何扩散的影响（Nolet，1987）。在反演Q值时，一般先用简化过的一维几何扩散修正观测到的地面运动振幅，然后再取对数，将其结果作为Q值反演的输入数据，所得到的输入数据称为“折衷振幅数据”（Menke et al，2006；Chen，Xie，2017）。即使通过一维几何扩散修正，目前未知的三维几何扩散影响依旧保留在输入数据中。Chen和Xie （2017）将这些残留的三维几何扩散的影响当作数据误差来处理并称其为建模误差，相较于振幅测量不准导致的误差，建模误差才是Q值数据的主要误差来源（Xie，1998）。本研究使用双台法测量Q值，并将两个台站相对地震的方位角差也引入了建模误差。
在许多地球物理反演研究中，直接运用中心极限定理隐含地假设输入数据的误差呈正态分布，反演问题就能简化至用最小二乘算法求解。中心极限定理指出：如果观测数据的误差来源很多，无论某一种来源的误差呈何种形态分布，所有来源误差的综合效应可以假设呈正态分布（盛骤等，2008）。如果一组数据在对数域呈正态分布，那么数据本身就可能不是正态分布，所以在使用振幅数据的Q值反演中，对于数据的正态性检验不可忽略。Chen和Xie （2017）首次开展了对“折衷振幅数据”建模误差的统计分析工作，使用中国东南地区的Lg波，得出“折衷振幅数据”的建模误差可以用正态分布描述的结论，此结论是仅适用于当前研究区或所使用的数据，还是一个具有普适性的结论，尚待商榷。本文将使用在川滇及其周边地区收集到的更为密集的地震数据，对Lg波的降幅数据误差进行更详尽的统计分析工作。
另一方面，地球物理层析成像方法主要包括同时迭代重构法（simultaneous iterative reconstruction technique）、反投影法（back-projection）和最小平方QR分解法（least-squares QR）等。当Q值层析成像所使用的振幅数据符合最小二乘算法的要求时，这些算法就能够被广泛使用在Q值层析成像中（Pei et al，2006；Xie et al，2006；Hearn et al，2008； Pasyanos et al，2009；Zhao et al，2010，2013；Zhao，Xie，2016）。这些方法都是利用一个简单的初始Q值模型，使用交互式的方式改进模型以使数据得到更好的拟合，从而使反演快速收敛，得到令数据残差最小化的Q值最终模型，然后再通过“棋盘格”测试（Gök et al，2003；Zhao et al，2010，2013；Liang et al，2014）、合成“点扩散函数”（Dai et al，2020）或者反映射线路径的采样密度的区域（Zhao ，Xie，2016）等各种分辨率测试方法来找出这些模型的近似分辨率和近似误差。需要指出的是，这些方法并不能定量地给出模型的分辨率和误差值，但是在实际应用这些模型时，定量的模型分辨率和误差是非常重要的。例如，对模型的分辨率和误差定量分析的结果有助于推断出所得到的Q模型的可靠部分，从而通过可靠的Q模型恢复特征来推断地壳属性。尽管理论上可以量化Q模型的分辨率和误差，并检查它们之间的影响（Backus，Gilbert，1970），但实际应用中它们会受到计算能力的限制。大规模矩阵奇异值（singular value decomposition，缩写为SVD）分解的数值计算需要使用迭代方法，但由于有限的计算精度，所得到的奇异矢量在迭代过程中会逐渐失去正交性，从而使最终的SVD分解是不完整的，即并非所有的奇异矢量都还保持正交。如果在迭代过程中保持所有奇异矢量的正交性，则需要庞大的计算量和内存。Chen和Xie （2017）开发了一种基于PROPACK算法的新的2D层析成像方法，该算法的优点在于它在迭代过程中只针对失去正交性的奇异矢量做调整，由于失去正交性的只是部分奇异矢量，因此计算量和内存需求大量减少.而正因为PROPACK算法得到的所有奇异矢量都保持了正交性，所以它们不仅可以用来计算模型参数，还可以计算模型的分辨率矩阵和协方差矩阵（Larsen，1998）。
本文拟利用双台法得到川滇及其周边地区地震数据的”折衷振幅数据”，并对其建模误差统计样本进行了提取和分析，采用 K-S （kolmogorov-Smirnov），Q-Q （Quartile-Quartile）图和正态分布图等三种统计检验方法得到研究区建模误差的统计特征，以期验证Q值反演数据中建模误差统计特征的普适性。在证明建模误差近似服从正态分布的基础上，使用最小二乘法求解Q值层析成像问题，并采用基于PROPACK 算法的 2D 层析成像方法得到模型的参数、分辨率矩阵与协方差矩阵，以期得到定量的模型分辨率和误差。最后，采用不同的层析成像方法来求解Q值模型，通过所得到模型结果的异同来检验新层析成像方法的有效性和可靠性。
1. 数据和方法
1.1 数据筛选
川滇地区位于华南地块、印度板块和青藏高原的交会处（95°E—108°E，22°N—35°N），区域内深大断裂纵横交错，地震活动强烈。川滇地区也是印度板块和欧亚板块碰撞带东南缘的一部分及特提斯—喜马拉雅造山系的转折点（邓起东等，2002），川滇地区的主要断裂带和地块分布如图1所示。本研究收集了来自国家数字测震台网数据备份中心（郑秀芬等，2009）提供的2007年7月至2013年5月期间262个固定台站（图1b）记录到的地震波形数据。台站分布在川滇及其周边地区（图1a），挑选了震中距范围在200—3 000 km、震级在M3.5—6.5之间且具有清晰的P波初至和Lg波波包发育的地震事件资料用于Lg波Q值的计算。通过手动拾取震相后，再进行去均值、去趋势、滤波等数据处理，选取133次地震的24 001条地震波垂直分量波形数据用来计算Lg波频谱。Wei和Zhao （2019）在川滇地区分别对频率为0.5，1.0和1.5 Hz的Lg波进行反演，得出在该区域内
$ {Q}_{\mathrm{L}\mathrm{g}} $ 没有很强的频率依赖性。所以在本研究中，我们只分析了1 Hz附近的Lg谱幅值数据，使用一个中心频率为1.0 Hz的窄带宽（0.5—1.5 Hz）来获得平均谱。![]()
图 1 研究所用地震事件震中（a）和研究区台站（b）分布图图中红色矩形框代表研究区域和台站所在区域；黑线为主要地块分界线，灰线为主要断裂带，下同Figure 1. Distribution of epicenters of events （a） and seismic stations （b） in this studyRed rectangular represents the extent of the study area；the black solid lines indicate the major tectonic block boundaries and the grey dashed lines indicate faults，the same below1.2 Q值的测量方法
在特定频率
$ f $ 处观测到的Lg波振幅可以通过正演建模得到，即：$$ A （ f ） ＝ S （ f ） G （ \varDelta ） \exp \left( { - \frac{{{\pi }f\varDelta }}{{vQ （ f ） }}} \right)R （ f ） \text{，} $$ (1) 式中
$ ，A （ f ） $ 表示振幅项，$ S （ f ） $ 为震源谱（Sereno et al，1988），$ G （ \mathrm{\varDelta } ） $ 是几何扩散项，$\varDelta$ 表示震中距，$ v $ 表示Lg波的群速度，$ Q （ f ） $ 为衰减值，$ R （ f ） $ 为台站的场地响应项（Street et al，1975；Xie，1993；Yang，2002；Chen，Xie，2017）。真实的几何扩散与3D地壳结构有关，但一维几何扩散可以简化为震中距$ \varDelta $ 的函数（Street et al，1975），表示为$$ G （ \varDelta ） ＝ \left\{ {\begin{array}{*{20}{l}} {\dfrac{1}{\varDelta }} \,\;\;\;\quad\qquad\qquad \varDelta ＜ {\varDelta _0}，\\ {\left({\dfrac{1}{{{\varDelta _0}}}} \right){{\left( {\dfrac{{{\varDelta _0}}}{\varDelta }} \right)}^m}} \qquad \varDelta {\text{≥}} {\varDelta _0}{\text{.}} \end{array}} \right. $$ (2) 对于Lg波来说，一维几何扩散的幂指数取0.5，临界距离
$ {\varDelta }_{0} $ 设为100 km。Q值是通过指数函数影响振幅的，经过1D几何扩散校正并取对数后得到$$ \ln {{A_i} （ f ） } － \ln {G （ {{\varDelta _i}} ） } ＝ \ln {{S_{ i}} （ f ） } － \frac{{{\pi }f{\varDelta _i}}}{{v{Q_i} （ f ） }} ＋ \ln {{R_i} （ f ） } \text{，} $$ (3) 式（3）等号左边是反演数据，称为“折衷振幅数据”。本研究使用双台法 （two station based method，缩写为TS），将式（1）应用于与震源位于同一大圆弧的两个台站，因为这两个台站记录了同一次地震，而路径衰减的差别主要体现在两个台站之间，所以两个台站的振幅比直接消除了震源谱，并求出两个台站之间的Q值（Xie，Mitchell，1990；Xie et al，2006；Zor et al，2007；Gallegos et al，2014；Chen，Xie，2017）。但传统双台法不能消除场地效应的影响，只是假设两个台站的场地效应相同。为了避免必须进行Q值和震源参数之间的权衡所导致的Q模型的非唯一性 （Menke et al，2006），研究中又使用了逆双台法（reversed two-station method，缩写为RTS）（Chun et al，1987）和逆双事件法（reversed two-event method，缩写为RTE）（Chen，Xie，2017）。逆双台法是在双台法的基础上，从双台路径的对拓方向再引进一次地震，两个台站记录的两个地震的交叉振幅比，可以将两个台站的场地效应也直接消除，从而得到两个台站之间的Q值。逆双事件法是逆双台法的拓展，即将逆双台法中台站和地震的位置互换就是逆双事件法，它可以测量两个地震事件之间的Q值。
2. 数据误差的统计分析
2.1 建模误差的统计方法
双台法的Q值计算表示为
$$ \frac{{\tilde \varDelta }}{{{Q_{\rm{e}}} （ f ） }} ＝ \frac{{v \ln \tilde A}}{{{\pi }f}} \text{，} $$ (4) 式中，
$ {Q}_{\mathrm{e}} （ f ） $ 为进行一维几何扩散校正后包含误差的$ Q （ f ） $ 的估计值，$\tilde{{\varDelta }}$ 为两个台站之间的有效距离，$ \tilde{A} $ 为双台法的振幅项与几何扩散项乘积的比值。本研究中将建模误差的定义范围进行了拓展，主要包括：① 残留的三维几何扩散导致的误差；② 不能完全去除台站响应部分造成的误差；③ 由Sn尾波引起的“噪声”和地震事件前无法模拟的背景噪声造成的误差；④ 两个台站和地震的方位角差别导致的误差。$ {Q}_{{\rm{e}}} （ f ） $ 中包含的误差可以表述为（Chen，Xie，2017）：$$ \varDelta \left[ {\frac{{\tilde \varDelta }}{{Q （ f ） }}} \right] ＝ \frac{{\tilde \varDelta }}{{{Q_e} （ f ） }} － \frac{{\tilde \varDelta }}{{Q （ f ） }} ＝ \frac{v}{{{\pi }f}}\varDelta （ {\ln {\tilde A} } ） ＝ \frac{v}{{{\pi }f}} （ 1 ＋ \varepsilon ） {\text{≈}} \frac{v}{{{\pi }f}}\varepsilon \text{，} $$ (5) 式中，
$\mathrm{ln}\tilde{A}$ 为“折衷振幅数据”，假设误差很小，对$\mathrm{ln}\tilde{A}$ 进行一阶近似，${{\tilde \varDelta }}/{{Q}_{{\rm{e}}} （ f ） }$ 的误差$\varDelta ［ {\tilde{\varDelta }}/{Q （ f ） } ］ $ ，可以用一个随机变量$\varepsilon$ 表示。由于岩石圈在学者们感兴趣的研究频率范围（＞0.1 Hz）是高度不均匀的，所以由$ \varepsilon $ 表示${\tilde{\varDelta }}/{Q （ f ） }$ 中的误差在不同的路径测量中是不相关的，即使是从同一个台站对测量收集的数据，它们的方位角（15°以内）也会略有变化.即使不同的震源位于台站相同的方位角上，仅仅只是震中距不同，也会因为局部波场的组成成分的差异导致产生不相关性（Mitchell，1995；Romanowicz，1998；Ji et al，2005；Gibbons et al，2010）。实验证明，使用不同地震与两个台站的组合，重复测量的两个台站之间Lg波的Q值的随机扰动平均值是没有任何规律的（Chen，Xie，2017），所以在研究过程中，对每一次${\tilde{{\varDelta }}}/{{Q}_{{\rm{e}}} （ f ） }$ 测量都提供一个随机变量$ \varepsilon $ 的样本。假设${\tilde{{\varDelta }}}/{{Q}_{{\rm{e}}} （ f ） }$ 的平均值接近真实的${\tilde{\varDelta }}/{Q （ f ） }$ ，可以用平均值减去单个样本值来近似的估计$\varDelta ［ {\tilde{\varDelta }}/{Q （ f ） } ］ $ ，从而根据式（4）得到单个$ \varepsilon $ 的值。2.2 建模误差样本采样
在计算1 Hz的频谱后，就可以求出“折衷振幅数据”
$ \mathrm{ln}\tilde{A} $ ，用此数据就可以估计两个台站间的$ {\tilde{\varDelta }}/{Q （ f ） } $ ，其中频率$f$ 设为1 Hz。在研究中分别得到了有效的9 633条双台法（TS）路径、3927条逆双台法（RTS）路径、2888条逆双事件（RTE）路径用于$ {\tilde{\varDelta }}/{Q （ f ） } $ 测量，其中有3 356条TS路径，722条RTS路径和302条RTE路径的唯一路径，在这些基于双台法的唯一路径中，分别有1 943条TS路径，547条RTS路径和253条RTE路径被多次采样，从而使得${\tilde{\varDelta }}/{Q （ f ） }$ 被重复测量.由于逆双台法和逆双事件法对台站和事件位置的要求比双台法更严格，使得可使用的数据大大减少，所以逆双台法和逆双事件法的有效路径相比于双台法要少很多。图2展示了三种方法的各个路径上
${\tilde{{\varDelta }}}/{{Q}_{{\rm{e}}} （ f ） }$ 的重复测量的次数分布图，最大重复测量次数达到100次以上.通过${\tilde{\varDelta }}/{Q （ f ） }$ 的重复测量来收集$ \varepsilon $ 样本，分别用$ {\varepsilon }_{\mathrm{T}\mathrm{S}} $ ，$ {\varepsilon }_{\mathrm{R}\mathrm{T}\mathrm{S}} $ 和$ {\varepsilon }_{\mathrm{R}\mathrm{T}\mathrm{E}} $ 来表示用TS，RTS和RTE方法估计的$ \varepsilon $ 值，如表1所示对于三种方法测量${\tilde{\varDelta }}/{Q （ f ） }$ 获得的$ {\varepsilon }_{\mathrm{T}\mathrm{S}} $ ，$ {\varepsilon }_{\mathrm{R}\mathrm{T}\mathrm{S}} $ 和$ {\varepsilon }_{\mathrm{R}\mathrm{T}\mathrm{E}} $ 的样本数量分别为8 220，3 652和2 839。![]()
图 2 基于双台法的各个路径上Q值重复测量次数分布（a） TS方法；（b） RTS方法；（C） RTE方法Figure 2. Distribution of the number of repeated two-stations based method over individual paths（a） TS method；（b） RTS method；（c） RTE method表 1 使用TS，RTS和RTE方法估计ε的详细信息Table 1. Estimated ε information using TS， RTS and RTE methods方法 $ {\varepsilon } $样本的
数量误差在两个标准差
内的$ \varepsilon $占比ECDF和NCDF的
相关系数K-S测试结果 标准差
std$ （ \varepsilon ） $方差Var$ （ \varepsilon ） $ P值 显著性水平$ \alpha $ TS 8 220 87.09% 99.99% 在置信度为95%下接受零假设 0.158 5 0.025 1 0.693 1 0.05 RTS 3 652 88.12% 99.98% 在置信度为95%下接受零假设 0.081 6 0.006 7 0.324 0 0.05 RTE 2 839 93.59% 99.64% 在置信度为95%下接受零假设 0.104 9 0.011 0 0.280 1 0.05 2.3 统计学方法检验结果
原始的
$ {\varepsilon }_{\mathrm{T}\mathrm{S}} $ ，$ {\varepsilon }_{\mathrm{R}\mathrm{T}\mathrm{S}} $ 和$ {\varepsilon }_{\mathrm{R}\mathrm{T}\mathrm{E}} $ 的样本分布具有中心峰值的特征，是显著的正态分布的特征，如图3左侧所示，然而远离峰值的分布在正负方向上都有较长的尾部，这些分布与Nolet （1987）讨论的远震P波走时残差分布非常相似；$ \varepsilon $ 的特征也符合Chen和Xie （2017）的统计特征结果。![]()
图 3 建模误差$\varepsilon$ 样本分布（左）和K-S检测图（右）（a） TS方法；（b） RTS方法；（c）RTE方法。左图中两条黑色直线标记了数据筛选准则|ε-median（ε）|＞2std（ε）的区间，红色线表示最佳的正态概率密度函数值；右图中灰色线代表着95%置信区间的CDF的上下界限，绿色小线段代表着K-S检验的经验CDF和理论CDF之间的最大垂向差值Figure 3. Distribution of samples of modeling error$\varepsilon$ （left）and results of K-S test （right）（a） TS method；（b） RTS method ；（c） RTS method. The black lines in the left figure mark the interval of |ε-median（ε）|＞2std（ε），which we use as criteria for data screening and the best-fit theoretical normal PDF of the screened samples is plotted over the distribution map with the red line. The grey lines in the right figure indicate lower and upper 95 percent confidence bounds for the CDF and the green segment indicates the maximum distance from the empirical CDF to the theoretical normal CDF which is used in the K-S test在对
$ \varepsilon $ 的样本进行了几次筛选试验后，确定了任何$ \varepsilon $ 偏离中位数超过两个$ \varepsilon $ 的标准差的样本都被剔除的筛选准则，即$\left| {\varepsilon － {\text{median}} （ \varepsilon ） } \right| ＞ 2 {{{\rm{std}} （ }}\varepsilon ） $ ，这也符合统计学上的“68-95-99.7”准则（Pukelsheim，1994），68.27% （n＝1），95.45% （n＝2）和99.73% （n＝3）的值是分别距离样本数据的平均值小于一个、两个和三个标准差内的百分比，超过三个标准差的数据可能是异常值，并且如果有许多数据超过三个标准差，则有理由质疑假定的正态性（Wheelan，2013）。两个标准差范围代表着95%的置信水平，可以同时平衡样本选择的精度和准确度。如表1所示，经过数据筛选后，对于所有的样本$ \varepsilon $ ，剔除的数据量不超过总体的12%。方差$ {\rm{Var}} （ {\varepsilon }_{\mathrm{T}\mathrm{S}} ） ，{\rm{Var}} （ {\varepsilon }_{\mathrm{R}\mathrm{T}\mathrm{S}} ） $ 和$ {\rm{Var}} （ {\varepsilon }_{\mathrm{R}\mathrm{T}\mathrm{E}} ） $ 分别为0.025 1，0.006 7和0.011 0，$ {\varepsilon }_{\mathrm{T}\mathrm{S}} $ 的方差远大于$ {\varepsilon }_{\mathrm{R}\mathrm{T}\mathrm{S}} $ 和$ {\varepsilon }_{\mathrm{R}\mathrm{T}\mathrm{E}} $ 的方差，说明$ {\varepsilon }_{\mathrm{T}\mathrm{S}} $ 的分布比$ {\varepsilon }_{\mathrm{R}\mathrm{T}\mathrm{S}} $ 和$ {\varepsilon }_{\mathrm{R}\mathrm{T}\mathrm{E}} $ 更广，$ {\rm{Var}} （ {\varepsilon }_{\mathrm{R}\mathrm{T}\mathrm{S}} ） $ 和$ {\rm{Var}} （ {\varepsilon }_{\mathrm{R}\mathrm{T}\mathrm{E}} ） $ 之间的微小差异很可能是由每个方法收集的重复测量样本数不同造成的，$ {\rm{Var}} （ {\varepsilon }_{\mathrm{T}\mathrm{S}} ） $ 的值最大可以归因于TS方法在测量中不能消除场地效应的影响。K-S （Kolmogorov-Smirnov）检验是一种非参数测试，测试分布的正态性时，将已知的正态概率分布与数据生成的分布（经验分布函数）进行量化比对（Massey，1951；Vrbik，2020），研究中的零假设为数据是正态分布。本研究采取单样本K-S检验， p值为否定零假设的最低显著性水平（Vrbik，2020）。K-S检验通过p值与显著性水平α进行比较，从而来决定是否拒绝零假设.如表1所示，p值都大于显著性水平，所以在95%的置信水平下接受零假设，建模误差是正态分布。剔除异常数据中新的ε样本的K-S检验图如图3，右侧图所示。使用筛选后的
$ {\varepsilon }_{\mathrm{T}\mathrm{S}} $ 样本计算累积分布函数（cumulation distribution function，缩写为CDF）如图3a右图所示，并与正态分布的CDF进行比较，观测和筛选的经验累积分布函数（empirical cumulation distribution function，缩写为ECDF）与理论最佳正态累积分布函数（normal cumulation distribution function，缩写为NCDF）之间的相关系数是99.99% （见表1），通过K-S检验不仅证明了筛选准则的可靠性，还证明了筛选后的$ {\varepsilon }_{{\rm{TS}}} $ 样本非常接近于正态分布函数。$ {\varepsilon }_{\mathrm{R}\mathrm{T}\mathrm{S}} $ 和$ {\varepsilon }_{\mathrm{R}\mathrm{T}\mathrm{E}} $ 的样本的检测结果如图3b，c右图所示，也得出了相同的结论。研究中还使用Q-Q （Quartile-Quartile）分位数图来检测
$ \varepsilon $ 是否服从正态分布。Q-Q分位数图是一种图形统计检验方法，用此方法可以找出两组数据是否来自同一分布。如果在Q-Q分位数图中纵轴表示理论正态分布分位数，横轴上输入验证数据，拟合的数据如果聚类在$y ＝ x$ 的线上，这表明样本数据呈正态分布（Marden，2004；Wheelan，2013）。如图4a所示，$ {\varepsilon }_{\mathrm{T}\mathrm{S}} $ 的Q-Q分位数图显示$ {\varepsilon }_{\mathrm{T}\mathrm{S}} $ 与正态分布基本一致，仅左下角和右上角存在一定偏差，代表着样本中剔除了长尾异常，即尾部未采样，但总体上可以证明$\varepsilon $ 服从正态分布。$ {\varepsilon }_{\mathrm{R}\mathrm{T}\mathrm{S}} $ 和$ {\varepsilon }_{\mathrm{R}\mathrm{T}\mathrm{E}} $ 的样本的检测结果如图4b，c所示，也得出了相同的结论。![]()
图 4 样本ε的Q-Q图（a） TS方法；（b） RTS方法；（c） RTE方法Figure 4. Q-Q plots of samples of$\varepsilon$ （a） TS method；（b） RTS method；（c） RTE method研究中还使用了其它图形检验方法如正态概率图以及数值检验方法如Shapiro-Wilk测试应用于检验筛选的
$\varepsilon $ 样本，所有的这些方法都是为了检验随机样本在一定置信水平下符合正态分布的假设而设计的，并互相验证测试结果。本研究中的检验结果与Chen和Xie （2017）开展的Q值模型建模误差的统计特征一致，但本次研究中使用K-S检验、Q-Q分位图并没有像Chen和Xie （2017）的结果一样对TS和RTS方法的正态性检测失效，因为这些方法是为相对较小（不超过几千）的样本量设计的，川滇研究区的Q值模型的建模误差样本都在几千的范围内，可以很好地呈现出正态性。所以通过多种统计检验方法，都证明了本文所采用的数据筛选准则是合适的，可以剔除建模误差样本的异常值，产生一个接近于完美正态分布的新的数据集，从而使用最小二乘算法求解Q值模型。3. 川滇地区Q值模型及评价
3.1 Q值层析成像方法
研究区Lg波Q值的横向变化可以用二维衰减层析成像解决，将研究区域划分为M个小单元来参数化Q值的横向变化，假设每个单元中的Q值恒定. 在收集N条唯一路径上的
${\tilde{\varDelta }}/{Q （ f ） }$ 测量值之后，对每条路径都可以构建线性方程（Chen，Xie，2017）$$ \frac{{{\varDelta _n}}}{{{Q_{{n}}} （ f ） }} ＝ \sum\limits_{m ＝ 1}^M {\frac{{{\varDelta _{nm}}}}{{{Q_m} （ f ） }}}， $$ (6) 式中，n＝1，2，···，N代表着第n条射线路径，m＝1，2，···，M代表着第m个模型的小单元。
$ {Q}_{n} （ f ） $ 是第n条射线路径上的Q（$ f $ ）的估计值，第n条射线路径长度为$ {{\varDelta }}_{n} $ ，当它穿过研究区域的第m个小单元的时候，相交的小段距离长度为$ {{\varDelta }}_{nm} $ 。式（6）的矩阵形式为$$ {{\boldsymbol{d}}_N} ＝ {{\boldsymbol{A}}_{N {\text{×}} M}}{{\boldsymbol{m}}_M} $$ (7) 式中，数据矢量d中的元素为
$ {{\varDelta }_{n}}/{{Q}_{n} （ f ） } $ ，转换矩阵中的元素为$ {\varDelta }_{nm} $ ，而模型参数矢量m中的元素为$ {1}/{{Q}_{m} （ f ） } $ 。转换矩阵A是一个稀疏矩阵，因为任意第n条射线路径穿过单元格数量比总的单元格数量M少得多。另外，实际反演的模型值是1/Q，只是按照习惯在图像中转换为Q。一旦证明了Q值反演的数据集受正态分布的误差影响，则可以使用最小二乘法来反演1/Q 值模型。数据的方差可以由
$ \varepsilon $ 的方差得出：$$ {\rm{Var}}\left( {\frac{{ \tilde{\varDelta} }}{{{Q_{\rm{e}}} （ f ） }}} \right) ＝ \frac{{{V^2}}}{{{{\pi }^2}{f^2}}}{\rm{Var}} （ \varepsilon ） {\text{.}}$$ (8) 为了更好地反映通过改变
${{\Delta }_{n}}/{{Q}_{n} （ f ） }$ 测量值之间的方差来量化模型的不确定性，通过式（8）可以计算基于双台方法的任何特定路径上的${\rm{Var}} ［ {\tilde{{\varDelta }}}/{{Q}_{{\rm{e}}} （ f ） } ］ $ 值。在研究中引入一个$ N{\text{×}} N $ 的数据协方差矩阵${\boldsymbol{C}}$ ，它的第n个对角线元素由${\rm{Var}} ［ {{\varDelta }_{n}}/{{Q}_{n} （ f ） } ］ $ 给出，然后归一化得到下式：$$ {{\boldsymbol{C}}^{ - \frac{1}{2}}}{\boldsymbol{A}}m ＝ {{\boldsymbol{C}}^{\frac{1}{2}}}{\boldsymbol{d}} {\text{.}}$$ (9) 最小二乘意义下求解模型
$ \mathit{m} $ 需要求解矩阵A的逆矩阵，但是在层析成像中这通常需要很大的计算量和内存，取而代之的是各种迭代技术，如反投影法（Xie et al，2004，2006）。反投影方法主要着力于迭代更新一个测试解$ \mathit{m} $ ，直到预测的数据不匹配而停止改善。研究中使用的新的2D层析成像方法不是迭代技术，而是通过PROPACK算法（Larsen，1998）对矩阵A进行精确的SVD分解，从而求得矩阵A的逆矩阵。3.2 川滇地区Q值模型和分辨率
可视化的棋盘格测试几乎是所有的层析成像方法均采用的一种评估反演可靠性的方法。为了测试研究使用的Lg波频谱振幅数据集能否很好地恢复Lg波衰减特征，使用棋盘格测试作为检测模型的可靠性方法之一，如图5所示。测试了不同网格尺寸的模型参数化范围，最终选择1°×1°作为网格反演大小。周连庆等（2008）研究表明，川滇地区最高Q值不超过800，最低Q值不小于100，因此设置了Q值交错为100和800的棋盘格来生成合成数据。合成数据的路径位置和数量与实际数据相同，然后通过合成数据来恢复原始的棋盘图案来评估反演的质量。基于SVD的新的层析成像方法可以完整地恢复出频率为1 Hz的
$ {Q}_{{\rm{Lg}}} $ 模型的棋盘格，并且恢复的棋盘格图案与射线覆盖范围（图6a）一致。![]()
图 5 （a） 1°×1°棋盘格结果；（b）研究区棋盘格恢复Figure 5. （a） The checkboard with 1°×1° grid；（b） The recovery of the original checkboard![]()
图 6$ {Q}_{{\rm{Lg}}} $ 模型的Backus-Gilbert扩散函数（a） 基于双台方法的射线路径覆盖图；（b） 用Backus-Gilbert扩散函数来表示研究区域的模型分辨率Figure 6. Backus-Gilbert spread functions of the tomographic Lg Q model（a） The path coverage of two-station based method over the studied region；（b） Spread function of the tomographic Q model尽管棋盘格测试说明反演过程是稳定可靠的，但它对模型分辨率大小的说明是相对模糊的，仅能说明射线覆盖的好坏，但是密集的射线覆盖并不一定意味着模型的高分辨率。一种极端情况是，如果模型地区仅仅被同一方向的射线覆盖，而缺乏交叉方向的射线，那么模型则不会被很好地解析。
新的层析成像方法同时给出了模型的分辨率矩阵，可以通过构建Backus-Gilbert扩散函数（spread function，缩写为SF）来定量评价模型的分辨率（Chen，Xie，2017），通过分辨率矩阵的元素所控制的均方根距离，可以对模型的有限分辨率进行近似度量，从而给出每个反演格点的可分辨范围。图6b显示出本研究中通过Backus-Gilbert扩散函数对川滇地区Q模型的分辨率进行评估，在地震射线覆盖良好的研究区域中，分辨精度可以达到100 km。通过棋盘格测试和Backus-Gilbert扩散函数都证明了研究中得到的川滇地区Lg波的1/Q模型是可靠的。
基于SVD的新的层析成像方法和反投影方法获得的川滇地区的地壳Q值模型如图7所示。由于两种方法原理不同，但都需要做光滑处理，为方便比较，两种方法采用了相同的光滑参数.在相同的光滑处理下，不同方法的
$ {Q}_{{\rm{Lg}}} $ 的区域变化似乎相当一致，最大值大于700，最小值小于80，横向非均匀性较强，模式具有与该地区速度结构相似的特征（如Tian 等（2021））。由图7可以看出，川滇地区西部及其邻近地区Q值普遍低于东部，四川盆地出现明显的东西向变化与Chen和Xie （2017）的结果相似，同时盆地东部的高Q值与Zhao等（2013）估算的整个四川盆地相似。扬子克拉通的Q值很高（低衰减），松潘甘孜褶皱带Q值很低（高衰减）。本研究中得到的$ {Q}_{{\rm{Lg}}} $ 结果与Wei和Zhao （2019）的川滇地区1 Hz的$ {Q}_{{\rm{Lg}}} $ 的区域分布基本一致。在滇中地块西侧和松潘—甘孜褶皱带西侧，反投影法得到的结果似乎有更精细的结果，但是从分辨率结果可以看出，滇中地块的Backus-Gilbert 扩散函数值（SF）约在100 km，而松潘—甘孜褶皱带西侧以及三江断层西侧的SF值在200至300 km。所以在相同的分辨率下进行构造讨论，新方法相比于反投影方法能够体现出更精细的地壳横向变化结果。不同的反演方法涉及不同的参数选择，在大多数情况下参数的设置是经验决定的，两种方法的比较过程中尽可能使用类似的参数，所以细节方面不尽相同。影响层析成像结果的因素很多，例如反投影方法需要设置一个初始模型，不同的初始条件，反演的结果也可能不同。新方法可以通过PROPACK算法直接对转换矩阵进行SVD分解，所以不需要初始模型，其优越性表现在模型分辨率的定量分析，从而决定了模型所反映的结构的可靠性。![]()
图 7 川滇地区及其邻近地区的Q值模型（a） 新的层析成像方法结果；（b） 反投影方法结果Figure 7. The Lg Q model of Sichuan-Yunnan area and its neighborhood（a） Result from the new tomographic method；（b） Result from the back projection method3.3 模型的不确定性
通过SVD奇异值分解，得到了Lg波1/Q 模型的协方差矩阵。因为证实了数据误差服从正态分布，误差的大小完全可以由误差的方差表示（式（8））。反演过程中，数据包含的误差会传递到模型中，反映在模型的协方差矩阵中，从而由数据误差导致的1/Q模型的误差可以由模型的协方差矩阵完全量化。Q值层析成像模型的不确定度由模型的分辨率和误差共同决定。一阶近似下，用δQ/Q来标定模型的相对误差，这个值可以使用协方差矩阵对角线元素的平方根来近似估计。在川滇地区，1 Hz的Q模型相对误差如图8所示，一个标准差下的模型的相对误差小于9%，在射线路径覆盖良好的区域，相对误差小于3%，这得益于研究中严格的数据筛选准则和基于SVD的新的层析成像方法。完整的Q值层析成像模型的协方差矩阵可以用来定量估计特定路径预测Q值的误差，当使用路径Q值校正振幅时，这些误差估计可用于建立地震事件识别和震源大小估计的置信水平。
4. 讨论与结论
使用川滇地区及其邻近地区2007年7月至2013年5月的地震数据，对Lg波的“折衷振幅数据”进行了严格的正态性检验工作，这项工作在应用最小二乘算法之前十分重要，但经常被忽略，尤其在应用类似“折衷振幅数据”这样被处理过的数据时。用一个随机变量
$ \varepsilon $ 来表示建模误差，将“折衷振幅数据”误差统计分析的复杂任务转化为更容易处理的对$ \varepsilon $ 样本的统计分析工作，获得了如下结论：1） Q层析成像使用“折衷振幅数据”作为输入数据导致的建模误差主要由几何扩散函数、地震事件前和相关相位前尾波的噪声以及场地效应的方位变化影响，通过多种统计检验方法证明，以前Q值成像研究中仅仅被隐含假设为正态分布的建模误差是真实服从正态分布的，尽管那些研究可能使用了不同的数据筛选方法。
$ \varepsilon $ 的分布以一个拟正态分布的峰为主，辅以向峰两侧扩散的弱长尾异常值。2） 通过恰当的数据筛选准则，可以剔除建模误差样本的异常值，从而产生一个接近于完美正态分布的新的数据集，用于在最小二乘意义下求解Q值模型。
3） 研究中所使用的川滇地区地震数据表明，方差
$ {\rm{Var}} （ {\varepsilon }_{\mathrm{T}\mathrm{S}} ） ，{\rm{Var}} （ {\varepsilon }_{\mathrm{R}\mathrm{T}\mathrm{S}} ） $ 和$ {\rm{Var}} （ {\varepsilon }_{\mathrm{R}\mathrm{T}\mathrm{E}} ） $ 分别为0.025 1，0.006 7和0.011 0，$ {\rm{Var}} （ {\varepsilon }_{\mathrm{T}\mathrm{S}} ） $ 值最大，这可以归因于TS方法在测量中不能消除场地效应的影响.4） 研究中使用了基于SVD分解的新的Q层析成像方法来反演Q值模型，相比于反投影等迭代方法，它最大优点是可以得到模型分辨率矩阵以及协方差矩阵，能够把模型的分辨率和误差定量化。
5） 新方法比反投影法可以获得更精细的地壳Q值横向结构，在射线覆盖良好的区域，新方法的QLg模型可分辨达到100 km左右，模型相对误差小于3%。研究区定量的不确定度结果表明，新的Q层析成像方法是稳定可靠的，可以用来解释地下介质的横向变化，同时通过使用模型的协方差矩阵，可以定量估计特定路径上Q预测的误差，这样的定量误差估计可以在未来继续开展在一定置信水平下的地震动预测、震源谱估计、震源类型识别和震源大小的估计。
国家测震台网数据备份中心为本文提供了数据支持，各位审稿专家为本文提出了宝贵的修改意见，作者在此一并表示感谢。
-
![]()
图 1 在震中距为R、方位角为φ的NEIC台上所测定震级大小的分布图
Figure 1. Magnitude distribution measured by NEIC stations where R is the epicentral distance and φ is the azimuth
![]()
表 1 国际地震机构速报的汶川地震的参数
Table 1 Parameters of the Wenchuan earthquake quickly reported by major international seismological institutions
序号 发震时刻 （UTC）
时：分：秒震中位置 震源
深度/kmMS 测定机构 北纬/° 东经/° 中文名称 代码 1 06：28：04.1 30.95 103.40 14 8.0 中国地震台网中心 CENC 2 06：28：00.9 31.10 103.30 10 7.8 美国国家地震信息中心 NEIC 3 06：27：59.3 31.10 103.30 10 8.0 俄罗斯科学院 RAS 4 06：27：58.9 31.10 103.20 10 7.5 欧洲地中海地震中心 EMSC 5 06：28：00.8 30.80 103.40 10 6.8 （mb） 罗马尼亚地球物理研究所 RNIEP 注：表中序号1数据来自中国地震台网中心 （2008a），其它数据来自Swiss Seismological Service （2008）。 
下载: 导出CSV
表 2 国际地震机构测定的汶川大地震的参数
Table 2 Parameters of the Wenchuan earthquake observed by major international seismological institutions
序号 发震时刻 （UTC）
时：分：秒震中位置 震源深
度/km震级类型 测定机构 数据来源 北纬/° 东经/° mb MS MW 名称 代码 1 06:27:59.5 31.01 103.42 14 6.4 8.2 中国地震台网中心 CENC 中国地震台网中心 （2008a） 2 06:28:01.8 31.00 103.32 19 6.9 8.1 7.9 美国国家地震信息中心 NEIC NEIC （2008b） 3 06:28:41.4 31.49 104.11 12 7.8 7.9 全球矩心矩张量项目
数据中心GCMT GCMT （2008） 4 06:27:59.0 31.10 103.20 10 7.9 欧洲地中海地震中心 EMSC EMSC （2008） 5 06:28:03.7 31.60 103.70 33 7.2 8.4 德国格拉芬堡地震
观测中心SZGRF SZGRF （2008） 注：全球矩心矩张量项目数据中心测定的是“矩心”（即所释放的地震矩的“时－空几何中心”）的位置，其物理意义与传统的震源位置及发震时刻（地震初始破裂的位置与时刻）不相同，不具有简单的可比性。
下载: 导出CSV
表 3 《地震数据报告》中列出的汶川地震的209个台站的参数和MSZ （NEIC，2008a）
Table 3 The related data of 209 stations and MSZ of Wenchuan earthquake in Earthquake Data Report（NEIC，2008a）
序号 台站代码 震中距/° 方位角/° 周期/s MSZ 序号 台站代码 震中距/° 方位角/° 周期/s MSZ 1 INCN 20.30 65.1 19.0 8.4 37 KONO 65.59 326.2 20.0 8.0 2 HIA 22.06 29.5 20.0 7.7 38 CTAO 65.42 135.3 21.0 7.7 3 MDJ 24.68 49.2 20.0 7.9 39 TIR 65.73 304.9 22.0 7.9 4 TKM2 25.03 306.1 22.0 7.8 40 AKUT 66.02 40.1 18.0 8.1 5 KULM 25.70 186.1 20.0 7.9 41 FOO 66.98 329.1 18.9 8.3 6 AAK 25.71 304.9 22.0 7.9 42 BER 67.32 327.8 18.0 8.3 7 AML 26.14 303.4 19.0 8.0 43 TRI 68.27 311.0 19.0 7.9 8 EKS2 26.22 304.6 22.0 7.8 44 MIDW 67.56 69.7 19.0 7.9 9 KURK 26.93 324.1 19.0 8.1 45 GRF 68.33 315.7 20.0 8.4 10 KKM 27.62 151.4 20.0 7.6 46 GRA1 68.33 315.7 20.0 8.4 11 MAJO 29.41 69.7 19.0 8.0 47 TIP 68.75 303.6 22.0 7.9 12 KSM 30.10 166.0 20.0 7.7 48 SDPT 68.43 37.6 20.0 8.1 13 BRVK 32.53 322.3 20.0 7.9 49 RER 69.10 227.5 22.0 7.7 14 ERM 33.58 59.7 22.0 7.8 50 COLA 69.81 25.5 20.0 8.1 15 YSS 34.14 50.8 19.0 8.3 51 CEL 69.73 303.0 21.0 7.7 16 YAK 35.44 21.2 19.0 8.0 52 AQU 69.82 307.8 20.0 8.0 17 ARU 40.10 322.6 22.0 7.9 53 OHAK 71.22 34.0 21.0 8.3 18 GUMO 41.88 105.0 20.0 7.6 54 VLC 71.09 310.6 22.0 8.1 19 COCO 43.39 189.3 19.0 7.9 55 KDAK 71.29 33.3 19.0 8.3 20 PET 45.29 44.4 21.0 8.3 56 ECH 71.38 315.4 22.0 7.9 21 GNI 47.76 297.9 20.0 7.7 57 WLF 71.34 317.1 21.0 8.0 22 DGAR 48.30 223.0 19.0 7.6 58 WDD 71.72 301.2 20.0 7.6 23 KIV 48.81 303.2 19.0 8.2 59 EGAK 72.27 24.0 19.0 8.3 24 BILL 51.55 25.3 19.0 8.2 60 TARA 71.87 99.6 21.0 7.6 25 MBWA 54.20 161.0 22.0 7.7 61 BNI 73.08 312.6 21.0 8.3 26 SMY 54.64 44.8 19.0 8.2 62 ABPO 73.31 235.1 19.0 7.8 27 BR13 56.21 299.5 21.0 7.6 63 MID 73.27 30.0 18.0 8.3 28 KEV 56.55 336.1 20.0 8.2 64 ESK 73.69 325.2 20.0 8.2 29 MSEY 57.80 240.8 21.0 7.8 65 LOR 73.79 315.6 18.5 8.4 30 WAKE 57.84 85.6 20.0 7.8 66 VSL 73.80 306.5 22.0 7.7 31 ISP 59.07 298.1 22.0 7.9 67 SSB 74.31 313.6 20.0 8.1 32 KBS 60.08 347.1 20.0 8.1 68 FLN 75.60 318.4 20.0 8.3 33 ADK 60.35 44.6 22.0 8.1 69 BORG 75.31 338.6 19.0 8.0 34 PSZ 63.55 311.7 21.0 7.8 70 RJF 76.16 314.6 19.0 8.2 35 NAO0 64.60 327.2 20.0 7.9 71 SKAG 77.73 26.5 18.0 8.3 36 NWAO 64.93 167.1 22.0 8.1 72 CAN 78.72 143.4 20.0 7.9 73 UCH 25.58 304.0 21.0 7.8 109 SNZO 97.50 133.9 20.0 7.8 74 SIT 79.33 28.3 21.0 8.2 110 LAO 98.19 20.1 20.0 8.2 75 SFJD 80.32 349.9 21.0 8.3 111 HOPS 98.07 34.7 21.0 8.1 76 CRAG 81.31 28.7 18.0 8.2 112 LKWY 98.75 23.7 22.0 8.4 77 PAB 82.91 312.1 19.0 8.1 113 RLMT 98.64 22.7 21.0 8.2 78 ESLA 82.59 312.0 19.0 8.8 114 AGMN 99.21 12.8 19.0 8.2 79 FUNA 82.43 104.6 20.0 7.6 115 BMN 99.61 30.2 22.0 8.1 80 TAU 83.91 149.1 22.0 7.8 116 EYMN 100.29 10.1 20.0 8.3 81 MTE 84.41 314.2 20.0 8.0 117 CMB 100.12 33.7 22.0 8.0 82 PAF 85.26 201.0 21.0 7.7 118 AHID 100.07 25.0 19.0 8.2 83 LSZ 85.50 249.3 22.0 7.8 119 ELK 100.19 28.8 21.0 8.2 84 SFS 85.88 310.4 20.0 8.5 120 BW06 100.62 24.0 20.0 8.2 85 KIP 86.24 67.4 22.0 7.6 121 HWUT 100.91 25.9 21.0 8.2 86 RTC 87.76 308.7 20.0 7.8 122 RSSD 101.17 19.7 19.0 7.6 87 POHA 89.09 67.5 20.0 7.8 123 DUG 101.64 27.5 19.0 7.9 88 NLWA 91.17 29.9 21.0 8.2 124 COWI 102.43 8.8 21.0 8.4 89 FFC 91.88 14.4 21.0 8.0 125 DBIC 101.94 285.3 20.0 8.0 90 LBTB 92.86 242.6 20.0 7.8 126 LIC 102.33 285.0 21.0 7.7 91 NEW 93.17 25.7 21.0 8.2 127 MAW 102.89 194.8 20.0 8.1 92 HAWA 93.77 28.1 20.0 8.1 128 ECSD 103.42 14.7 19.0 8.3 93 SCHQ 94.10 354.3 20.0 8.2 129 OGNE 104.67 19.6 20.0 8.2 94 MSO 95.54 24.7 22.0 8.0 130 ISCO 104.61 22.7 19.0 8.2 95 EGMT 95.92 21.6 21.0 8.1 131 LONY 104.72 358.5 20.0 8.1 96 HUMO 95.30 32.4 22.0 8.1 132 JFWS 105.40 10.3 19.0 8.2 97 BMO 95.93 27.8 20.0 8.0 133 NCB 105.36 358.2 20.0 8.3 98 YBH 96.10 32.8 19.0 7.9 134 SCIA 105.87 12.8 21.0 8.3 99 TSUM 96.16 251.4 21.0 7.9 135 MVCO 105.91 26.0 19.0 8.1 100 CASY 97.13 177.1 21.0 8.0 136 WUAZ 106.23 28.9 20.0 8.0 101 ULM 97.33 12.3 20.0 8.3 137 SDCO 106.50 23.5 21.0 8.1 102 DGMT 97.14 18.1 22.0 8.3 138 BINY 107.16 359.5 20.0 8.2 103 WDC 97.08 33.4 22.0 8.0 139 KSU1 107.83 16.2 20.0 8.3 104 BOZ 97.38 23.9 20.0 8.1 140 HDIL 107.85 10.1 20.0 8.2 105 MOD 97.07 31.3 19.0 8.1 141 ANMO 108.66 25.5 22.0 8.0 106 XMAS 96.86 83.5 19.0 7.7 142 ACSO 108.89 5.1 20.0 8.2 107 WVOR 97.36 29.9 20.0 8.0 143 CCM 109.97 12.2 19.0 7.8 108 HLID 98.15 26.7 20.0 8.1 144 TUC 109.23 30.2 21.0 8.1 145 WVT 112.41 9.8 19.0 8.1 178 LVZ 53.66 334.0 21.0 8.5 146 WCI 110.54 8.1 19.0 7.4 179 JOHN 79.19 77.0 21.0 7.8 147 AMTX 110.36 21.8 21.0 8.1 180 MCCM 98.82 35.2 20.0 8.0 148 MNTX 111.96 26.2 20.0 8.0 181 SAO 100.61 35.1 20.0 7.9 149 WMOK 111.38 19.4 20.0 8.4 182 TPH 101.55 31.6 20.0 7.9 150 SBA 114.81 168.0 19.0 8.0 183 DAC 102.85 32.9 20.0 8.0 151 MIAR 112.89 15.1 21.0 8.2 184 MVU 103.33 27.9 22.0 8.1 152 NATX 115.29 16.9 19.0 8.3 185 GLMI 104.19 5.8 20.0 8.2 153 GOGA 115.61 6.3 22.0 8.3 186 LBNH 105.00 356.5 19.0 8.1 154 VBMS 115.77 13.0 19.0 8.5 187 AAM 106.78 5.4 19.0 8.2 155 LRAL 115.57 9.6 19.0 8.3 188 CBKS 107.30 18.7 21.0 8.2 156 HKT 116.71 18.6 20.0 7.6 189 MCWV 109.64 2.6 20.0 8.4 157 BRAL 117.54 10.0 19.0 8.2 190 BLA 112.05 3.2 20.0 8.4 158 KVTX 118.36 21.4 19.0 8.3 191 CNNC 114.10 1.1 22.0 8.2 159 DWPF 121.04 4.9 20.0 7.9 192 SHEL 114.20 265.4 20.0 7.6 160 GRTK 127.52 353.5 20.0 8.0 193 BBSR 115.96 348.7 22.0 8.1 161 TGUH 134.06 14.4 19.0 8.0 194 NHSC 116.15 3.3 20.0 8.3 162 ANWB 129.54 341.5 20.0 8.0 195 MOL 65.53 329.6 23.1 8.2 163 SDDR 130.03 353.3 21.0 8.2 196 KMBO 70.23 256.2 20.0 7.5 164 SJG 130.13 346.9 20.0 8.4 197 COR 93.6 31.4 21.0 7.3 165 GTBY 129.34 358.1 22.0 8.2 198 PFO 105.57 12.8 21.0 7.5 166 MTDJ 131.05 1.1 20.0 8.1 199 BJO1 58.90 341.8 17.7 8.5 167 FDF 132.14 339.6 21.0 8.1 200 TRO 59.37 336.1 13.8 8.5 168 JTS 138.18 12.3 21.0 8.1 201 NSS 63.04 331.0 17.2 8.3 169 GRGR 134.75 339.1 22.0 8.1 202 PMR 70.95 28.8 17.0 8.2 170 RCBR 134.32 294.7 19.0 8.0 203 CLL 66.67 316.9 22.1 8.1 171 BCIP 139.95 4.9 22.0 7.9 204 CPUP 162.4 280.5 21.0 7.2 172 PMSA 145.30 189.5 21.0 7.9 205 NNA 161.08 0.5 21.0 7.8 173 PAYG 147.19 25.7 20.0 7.9 206 LPAZ 163.4 330.0 22.0 8.1 174 OTAV 148.89 3.4 20.0 7.9 207 TRQA 166.0 235.9 21.0 8.2 175 RPN 151.25 90.6 21.0 7.8 208 PLCA 169.1 205.4 19.0 7.9 176 SPB 153.04 278.7 20.0 8.2 209 LCO 174.5 289.4 20.0 8.2 177 EFI 155.20 208.3 21.0 8.0 平均 8.1 注：表中平均值计算未包括增加的11个台站资料。
下载: 导出CSV
-
陈运泰. 2008. 汶川地震的成因断层、破裂过程与成灾机理[C]//中国科学院第十四次院士大会学部学术报告汇编. 北京: 中国科学院: 38–39. Chen Y T. 2008. Generating faults, rupture process and disaster mechanism of the Wenchuan earthquake[C]//Collection of the Academician Reports on the 14th General Assembly, Chinese Academy of Sciences. Beijing: Chinese Academy of Sciences: 38–39 （in Chinese）.
陈运泰, 许力生, 张勇, 杜海林, 冯万鹏, 刘超, 李春来, 张红霞. 2009. 汶川特大地震震源特性分析报告[G]//汶川大地震工程震害调查分析与研究. 北京: 中国岩石力学与工程学会: 6–17. Chen Y T, Xu L S, Zhang Y, Du H L Feng W P, Liu C, Li C L , Zhang H X. 2009. Analysis report of the source characters of the great Wenchuan earthquake[G]//Investigation, Analysis and Research on the Engineering Earthquake Damages of the Great Wenchuan Earthquake. Beijing: Chinese Society for Rock Mechanics and Engineering: 6–17 （in Chinese）.
郭履灿, 庞明虎. 1981. 面波震级和它的台基校正值[J]. 地震学报, 3 （3）: 312-320. Guo L C, Pang M H. 1981. Surface-wave magnitude and its station correction [J]. Acta Seismologica Sinica, 3 （3）: 312-320 （in Chinese）.
胡聿贤. 2006. 地震工程学[M]. 北京: 地震出版社: 54–56. Hu Y X. 2006. Earthquake Engineering [M]. Beijing: Seismological Press: 54–56 （in Chinese）.
刘瑞丰, 陈运泰, 任枭, 徐志国, 王晓欣, 邹立晔, 张立文. 2015. 震级的测定[M]. 北京: 地震出版社: 1–10. Liu R F, Chen Y T, Ren X, Xu Z G, Wang X X, Zou L Y, Zhang L W. 2015. Determination of Magnitude [M]. Beijing: Seismological Press: 1–10 （in Chinese）.
刘瑞丰, 陈运泰, Peter Bormann, 任枭, 侯建民, 邹立晔. 2006. 中国地震台网与美国地震台网测定震级的对比(Ⅱ): 面波震级[J]. 地震学报, 28（1）: 1-7. Liu R F, Chen Y T, Bormann P, Ren X, Hou J M, Zou L Y. 2006. Comparison between earthquake magnitudes determined by China seismograph network and U.S. seismograph network （Ⅱ）: Surface wave magnitude [J]. Acta Seismologica Sinica, 28 （1）: 1-7（in Chinese）.
许绍燮. 1999. 中华人民共和国国家标准《地震震级的规定》（GB 17740—1999）宣贯教材[M]. 北京: 中国标准出版社: 1–10. Xu S X. 1999. Teaching Book for Propagandizing and Implementing the National Standard of the People’s Republic of China: General Ruler for Earthquake Magnitude （GB 17740–1999） [M]. Beijing: Standards Press of China: 1–10 （in Chinese）.
袁一凡. 2008. 四川汶川8.0级地震灾害损失评估[J]. 地震工程与工程振动, 28（5）: 10-19. Yuan Y F. 2008. The earthquake disaster loss evaluation[J]. Earthquake Engineering and Engineering Vibration, 28（5）: 10-19（in Chinese）.
张勇, 许力生, 陈运泰. 2009. 2008年汶川大地震震源机制的时空变化[J]. 地球物理学报, 52 （2）: 379-389. Zhang Y, Xu L S, Chen Y T. 2009. Spatio-temporal variation of the source mechanism of the 2008 great Wenchuan earthquake [J]. Chinese Journal of Geophysicis, 52 （2）: 379-389（in Chinese）.
中国地震台网中心. 2008a. 地震信息[DB/OL]. [2008–05–12]. http://www.cenc.ac.cn. China Earthquake Networks Center. 2008a. Earthquake information [DB/OL]. [2008–05–12]. http://www.cenc.ac.cn （in Chinese）.
中国地震台网中心. 2008b. 中国地震台站观测报告[R]. 北京: 中国地震台网中心: 1–10. China Earthquake Networks Center. 2008b. Seismological Report of the Chinese Seismic Stations[R]. Beijing: China Earthquake Networks Center: 1–10 （in Chinese）.
中国地震台网中心. 2008c. 中国数字地震台网观测报告[R]. 北京: 中国地震台网中心: 1–10. China Earthquake Networks Center. 2008c. Seismological Report of the Chinese Digital Seismic Network[R]. Beijing: China Earthquake Networks Center: 1–10 （in Chinese）.
周庆, 张春山, 陈献程. 2011. 汶川MS8.0地震灾害的非对称分布与成因分析[J]. 地震学报, 33（4）: 492-504. Zhou Q, Zhang C S , Chen X C. 2011. Asymmetric disaster distribution and its cause analysis of the MS8.0 Wenchuan earthquake. Acta Seismologica Sinica, 33 （4）: 492-504（in Chinese）.
EMSC. 2008. M7.9 − Eastern Sichuan, China − 2008-05-12 06:28:00 UTC[EB/OL]. [2017–10–22]. https://www.emsc-csem.org/Earthquake/earthquake.php?id=85969.
GCMT. 2008. Monthly CMT solutions[EB/OL]. [2017–10–22]. http://www.ldeo.columbia.edu/~gcmt/projects/CMT/catalog/NEW_MONTHLY/2008/mar08.ndk.
Gutenberg B. 1945.Amplitudes of surface waves and magnitudes of shallow earthquakes[J].Bull. Seism. Soc. Am. 35: 3-12.
NEIC. 2008a. The preliminary determination of epicenters (PDE) bulletin[EB/OL]. [2014–10–01]. ftp://hazards.cr.usgs.gov/NEICPDE/olderPDEdata/mchedr/.
NEIC. 2008b. M7.9−eastern Sichuan, China [EB/OL]. [2017–10–22]. https://earthquake.usgs.gov/earthquakes/eventpage/usp000g650
Somerville P G, Smith N F, Graves R W, Abrahamson N A. 1997.Modification of empirical strong ground motion attenuation relations to include the amplitude and duration effects of rupture direction[J]. Seis Res Lett, 68（1）: 199-222.
Swiss Seismological Service. 2008. Recent earthquakes Switzerland[DB/OL]. [2008–05–12]. http://www.seismo.ethz.ch.
SZGRF. 2008. Seismological center observatory[DB/OL]. [2008–07–01]. https://www.szgrf.bgr.de/bulletins.html.
计量
- 文章访问数: 1540
- HTML全文浏览量: 946
- PDF下载量: 62
