Research on the earthquakes in early Eastern Han Dynasty
-
摘要: 东汉早期(公元25—127年)的地震史料十分宝贵， 它反映了从远古零散、 不系统的地震记载转入连续记载并逐渐认识地震的过程. 这时的地震事件较多， 但研究不充分， 地震参数基本呈空白状. 史料的记载方式以“郡国地震”和“京师地震”为主； 时间上， 基本在汉和帝以后； 地点上， 以大华北为主. 本文具体分析了东汉早期的33次地震事件, 利用新的地震烈度衰减关系， 采用郡国平均地理面积和有感面积的概念， 比照地震范例的标准， 按《中国地震目录》的烈度震级表进行参数估算. 研究期间最主要的地震事件是4次震级≥6 1/2 的强震， 即公元46年10月和119年3月南阳地震、 121年10月冀南-鲁西地震和123年5月汉阳地震， 震后都有持续1—2年的余震活动. 公元118年前的地震主要分布于大华北， 震中至京师洛阳的距离多在200 km以外. 除两次5 1/2 —6级地震外， 均属5—5 1/2 级中等强度地震， 少数小于4 3/4 级. 京师洛阳虽经历过多次地震， 但其受影响程度并不强， 基本为有感或强有感的水平， 没有遭受过破坏性或中等强度的震害.Abstract: The historical literatures of earthquakes in early Eastern Han Dynasty (AD 25—127) are valuable. They reflect the transfer-process of earthquake documentary from scattered and nonsystematic records in early years to continuous records and being understood gradually. The documented earthquake events in this period are rather plentiful. However, the relevant information about these earthquakes has not been studied deeply and the parameters of few earthquakes in this period have been obtained so far. The existing investigations show that the reporting manners in this period are mainly expressed as “Junguo earthquake” and “Capital area earthquake”. The documented events occurred mainly in the region of the Great North Chinaafter He Emperor of Han Dynasty. In this paper, a new relationship of seismic intensity attenuation is used. By means of the intensity-magnitude table in Earthquake Catalogue of China, the earthquake parameters in this period are evaluated following the conceptions of Junguo-average geographic area and seismic felt area. In the light of earthquake example standard, 33 earthquake events in early Eastern Han Dynasty are analyzed particularly. The most important events in the studied period are four strong earthquakes with magnitude ≥6 1/2 , i.e. , the Nanyang earthquakes in October of AD 46 and March of AD 119, the south Hebei-West-Shandong earthquake in October of AD 121 and the Hanyang earthquake in May of AD 123; their aftershocks all occurred in 1—2 years. On the other hand, the earthquake events before AD 118 distribute principally in the Great North China, and their epicentral distances from the capital Luoyang are all larger than 200 km. Except for the two events with magnitude of 5 1/2 —6, all the others are of moderate magnitude of 5—5 1/2 . Though the capital Luoyang underwent many earthquakes, it has not subjected to destructive or medium catastrophes. The effects caused by these earthquakes in this period are not strong and they are basically seismically felt or strong felt events.
-
Keywords:
- historical earthquake /
- Junguo earthquake /
- magnitude estimation /
- capital Luoyang
-
引言
地壳中的应力是地震和其它构造运动过程的直接动力来源． 区域地壳应力场及其时空变化的研究为深入了解地震机理、 构造负载及地震应力的相互作用提供了可靠的信息．
岩石圈内真实构造应力的计算和测量都是非常困难的，但可以利用其它的一些地球物理观测资料反演出构造应力的某些特征，如应力方向、 应力相对大小等. 这一方面国内外早已开展了大量的研究，并取得了丰硕的成果(Zoback，1992； 谢富仁等，2004)． 20世纪80—90年代，许忠淮等利用小震初动符号采用区域综合震源机制解的方法，对华北和中国大陆其它地区的构造应力方向进行了大量的研究，取得了很多具有开创意义的成果(许忠淮， 1983，1989； 许忠淮，1985)； 利用野外观测的第四纪断层活动擦痕数据反演构造应力张量，也可以得到区域构造应力场的基本特征(Nemcok，Lishle，1995； 谢富仁， 1993，2001； 荆振杰等，2008； 张红艳等，2009). 此外，还可以使用跨断层形变观测资料(谢富仁等，2008)以及钻井的井孔崩落特征或水压致裂的情况对构造应力场进行研究(许忠淮，吴少武，1997； 许忠淮等，1999). 这些研究所得到的构造应力场的方向已经为近年来GPS观测的结果所证实(袁金荣等，1999； 刘峡等，2006)． 值得一提的是，1986—1990年，在国际岩石圈计划下开展了“世界应力图”项目的合作研究，综合各种资料完成了全球范围内应力图的绘制，取得了具有里程碑意义的成果(Zoback，1992； 许忠淮，1994)．
使用大量的地震震源机制解资料可以推断区域应力场的特征，对此地震学家已经发展了许多经典的方法. 例如，Michael(1984)，1987)的LSIB(linear stress inversion with bootstrapping)方法，Gephart和Forsyth(1984)的FMSI(focal mechanisms stress inversion)方法，Angelier(2002)基于滑动剪切应力分量的反演方法，以及崔效锋和谢富仁(1999)逐次收敛法等． 这些方法在实际研究中都有着广泛的应用(Petit et al，1996; Plenefisch，Bonjer，1997; Lund，Slunga，1999; Hardebeck，Hauksson，2001;Bressan et al，2003;Kastrup et al，2004;Boness，Zoback，2006;Carena，Moder，2009)．
中小地震的震源机制参数通常有很大的不确定性，但是大量的离散分布的数据可以准确地约束应力张量的方向(Hardebeck，Michael，2006)． 近年来，使用震源机制解推断构造应力特征的方法又有了新的发展(Loohuis，van Eck，1996;Abers，Gephart，2001)，一方面基于大量的资料进行小区域精细构造应力场的研究(Hardebeck，2006；Fay et al，2008)，另一方面通过增加阻尼约束、 应用概率信息准则等数学手段来更好地约束反演结果(Hardebeck，Michael，2006；Arnold，Townend，2007；Terakawa，Matsu ’ura，2008).
自20世纪80年代的“七五”、 “八五”攻关中集中计算小震震源机制以来，利用初动符号或振幅比资料求解区域中小地震的双力偶模型参数的工作逐渐成为国内分析预报部门的日常工作，积累了大量的小震震源机制解的资料． 尤其是近年来，随着数字观测台网的加密，区域小震震源机制的结果日见丰富，其质量也更加可靠(胡新亮等，2004). 也由此对不同区域开展了构造应力场的研究(吴建平等，2004； 康英等，2008； 李瑞莎等，2008； 刘泽民等，2011； 李莹甄等，2011)． 我们利用山东省地震预报研究中心所掌握的资料，使用Hardebeck和Michael(2006)的区域应力张量阻尼反演方法对山东地区的背景应力场进行研究，并分析探讨区域应力场与区域内的深大断裂 郯庐断裂带的关系．
1. 数据与资料
1.1 震源机制解数据
前人对山东地区中小地震的单震震源机制已经开展了大量的工作(魏光兴，李秉锋，1980； 魏光兴等，1982； 周翠英，魏光兴，1987； 魏光兴等，1993； 董旭光等，1999； 周翠英， 2003，2005； 山长仑等，2007； 郑建常等，2010； 郑建常，陈运泰，2012)． 自2006年以来，随着数字化测震台网的投入使用，山东省地震预报研究中心使用中国地震局数字地震实验室推广的Snoke的FOCMEC程序(Snoke et al，1984; 刘杰等，2004)以及MapSIS软件中的初动符号格点尝试法程序进行例行的震源机制求解，基本上每次省内ML≥2.0地震后都会给出相应的震源机制解． 通过搜集汇总这些资料，我们总计得到了山东地区1970—2012年5月共387次中小地震的震源机制解(图 1a)． 采用Frohlich(2001)的震源机制量化分类和评估方法(图 1b)，可以看出这些小震的震源机制类型主要是走滑和正断的断层作用，合计约占总数的76%．
![]()
图 1 本文使用的小震资料空间分布(a)及震源机制解分类(b)图Figure 1. Locations of 387ML≥2.0 earthquakes used in this paper (a) and classification of focal mechanisms on the octant equal-area projection according to Frohlich (2001) (b)1.2 山东地区整体应力场特征
图 2给出了本文所使用的震源机制解的P轴和T轴在震源球上的投影分布. 可以看出其P轴方位的优势分布为ENE-WSW向，T轴方位的优势分布为NNW-SSE向． 将P，T轴方位归一化到[0，180°]之间，统计研究区内所有震源机制解的P，T轴方位的优势方向(图 3)，结果得到山东地区的平均P，T轴取向的均值及标准差分别为P=86.9°±38.6°，T=173.2°±38.5°．
![]()
图 2 本文使用的小震震源机制解的P轴(a)、T轴(b)下半球极射投影图Figure 2. Lower hemisphere projection ofPaxes (a) andTaxes (b) used in this study![]()
图 3 本文使用的小震震源机制解的归一化P轴方位(a)与T轴方位(b)统计直方图Figure 3. Histogram of normalized azimuths ofPaxes (a) andTaxes (b) of the focal mechanism solutions used in this study使用ZMAP软件包中的FMSI程序反演山东地区整体的应力场，得到3个应力主轴的方向(图 4). 结果表明，山东地区的最大主应力方位为82°，最小主应力方位为172°，与前面统计结果相吻合，应力相对大小的平均值R=(σ1-σ2)/(σ1-σ3)=0.67．图 4给出了3个应力主轴方向的95%置信区间． 六省市震源机制小组(1981)对山东及周边应力场的研究表明，苏鲁皖豫地区处于北东东向水平压应力和北北西向水平张应力的现代构造应力场作用之中，其地震应力场的优势分布为主压应力轴约76°±28°，主张应力轴约344°±27°. Wan(2010)基于中国地壳应力数据库和美国哈佛大学矩心矩张量目录给出了中国大陆现代构造应力场. 其中山东地区的最大主应力轴方位为70°—90°，最小主应力轴方位为160°—194°. 本文结果与上述结果基本一致．
![]()
图 4 山东地区主应力轴及其置信区间的下半球等面积投影Figure 4. Lower-hemisphere equal-area projections of principle stress axesσ1， σ2andσ3， and their confidence regions2. 理论与方法
地球内部任一点的构造应力状态可以由6个独立的参数完全定义，表示为6个分量的3×3对称矩阵或者3个特征向量与地理坐标的夹角以及对应的3个特征值． 在理想情况下，钻孔应力测量可以推断6个应力张量的参数，但在通常的地质学和地震学研究中，一般只能给出其中的4个参数，即3个主应力的方向和中等主应力相对大小的参数．
受地下介质结构、 地质构造条件及断层展布情况等因素的影响，一定区域的应力空间分布是不均匀的(如王勤彩等(2009)对汶川地震余震区应力状态分段特征的研究)． 当需要考察应力状态的空间分布特征时，通常的方法是将区域分成很多小区，对每个小区的震源机制拟合应力张量(Petit et al，1996 Hardebeck，Hauksson，2001;Bressan et al，2003)． 分区方式不同和分区大小的变化都可能对反演的应力张量结果产生影响(Maury，Cornet，2011)． 并且由于每个小区都是独立反演，在某些情况下，计算结果显示出来的应力偏转等特征有可能相当的大． 相邻两个区内主应力方向的不同或许反映了应力场真实的变化，但也可能是震源机制数据的误差或反演的约束较差造成的假象． 受实际地震资料的限制，某些小区的震源机制可能对一组应力张量都拟合得很好，因而对这样的一个小区就有多个应力模型可以满足数据，如何增加约束来选择更合适的应力模型就成为必要的问题(Hardebeck，Michael，2006)．
由于地壳中的应力分布存在连续性，因此，为了解决反演的应力模式依赖于分区的问题，并且为了更准确地得到区域范围内的应力变化的特点，Hardebeck和Michael(2006)参考其它地球物理领域中的阻尼最小二乘反演方法，提出了区域尺度的应力张量阻尼反演方法． 他们构建一组可调整的阻尼参数模型，引入平滑约束来抑制相邻单元之间应力模式的差异，然后使用最小二乘法可以得到稳定解．
2.1 阻尼区域尺度应力反演
构造应力张量有6个独立分量，对于一组地震的震源机制，设应力张量不存在各向同性分量时，σ33=-(σ11+σ22)，则应力张量只有5个分量
Michael(1987)的反演方法通过最小化每个断层面上滑动向量和解得的剪切应力向量的差异，求解方程
式中，d 为由震源机制解求得的滑动矢量
式中，sij是第i次地震的单位滑动矢量的第j个分量(j=1，2，3)． 矩阵 G 由每个震源机制的断层法线矢量组成．
对于区域应力场，设研究范围分为I×J个二维空间网格． 第i行第j列的格点上的应力张量由向量 mij表示，各点上的观测数据由向量 dij表示，数据核矩阵记为 Gij，则新的模型向量
新的数据向量
新的数据核矩阵
则对于区域应力场求解问题，方程(2)可写为
对此方程增加“平滑”约束，最小化相邻各点的应力张量的差异，引入阻尼矩阵
式中，I是5×5的单位矩阵．
阻尼反演最小化两个值的权重和为数据错配、 模型长度． 模型长度表示为向量的二阶范数，包含了每对相邻格点的每个应力张量的差异
数据与模型理论值的错配写为
阻尼最小二乘问题的解，需要同时最小化等式(9)和式(10)，即
式中，e是标量阻尼参数，控制最小化问题中数据错配和模型长度的相对权重． 由于 G，D 都非常稀疏，大括号内的矩阵也是稀疏矩阵，故使用共轭梯度方法求解等式(11) ．
2.2 “真实”最大水平主压应力方向的确定
记最大水平主应力为σHmax. 识别最大水平应力的常用方法是简单地采用最大近水平主应力的水平投影方向，但该情况只有当应力张量中的其中一个主应力轴严格垂直时是适用的． 应力的Anderson模型可以很好地解释不同构造条件下观测到的应力状态. 该情况下，σHmax与最大主应力轴平行，除非该轴相对于地表水平面是垂直的(此时，σHmax平行于介质应力轴)． 但是，一般情况下由地震数据估计得到的应力张量的结果中，几乎没有应力轴是严格垂直于地表水平面的． 这时，σHmax不能简单地采用最大近水平应力的水平投影，二者之间在方向上可能存在最大数十度的偏差(Lund，Townend，2007)． Lund和Townend(2007)从数学上推导给出了由应力张量严格计算σHmax方向的公式和方法，可以得到“真实的”最大水平主应力方向．
已知构造应力张量 S，由此可以得到3个主应力σ1≥σ2≥σ3及其对应的主应力矢量 s1，s2，s3，在Aki和Richards(1980)定义的地理坐标系下(北、 东、 下为正)
作用在我们感兴趣的垂直平面上的法向正应力Sn与该法向应力的方位角α之间存在如下关系(Lund，Townend， 2007，式5)：
式中，α为法向正应力Sn的方位角． 一般情况下，Sn在区间[0，π)具有一个最大值和一个最小值，使得上式为零，分别对应了水平最大主应力和最小主应力．
3. 计算结果与分析
3.1 阻尼系数的确定
在对公式(11)进行实际反演计算过程中，阻尼参数e的取值对反演结果起着重要的作用．由于该参数控制着理论值与观测数据之间错配值和模型长度(即模型复杂程度)的相对权重，简化模型，则错配值升高，反演误差增大； 反之提高错配的相对权重，则反演误差减小，而模型逐渐变得复杂化，甚至失去了阻尼约束的意义． 因此通常是根据如图 5所示的折中曲线来选择阻尼参数． 折中曲线的拐点在e≈0.5附近，意味着低于该相对权重值，提高模型复杂程度对反演误差的改善几乎没有作用; 而提高相对权重时，随着模型的简化，反演误差会急剧增加． 因此，我们在反演过程中选择相对权重系数e≈0.5.
![]()
图 5 模型长度与数据拟合误差之间的折中曲线图虚曲线旁所标数字是阻尼参数e的取值(见式(11))Figure 5. Trade-off between model length and misfit Numbers along the dashed line are values of damping parameterse3.2 山东地区背景应力场
对山东地区(33.5°—38.5°N，114°—124.5°E)划分0.5°×0.5°的网格，使用前面引述的Hardebeck和Michael(2006)的区域应力场阻尼反演方法，计算每个网格的最佳拟合构造应力张量; 然后使用Lund和Townend(2007)给出的公式计算最大水平主应力的方向． 为了平滑地震空白的区域，没有数据的格点也包括在反演中，但这些点上的结果没有显示． 最终反演结果见图 6．
![]()
图 6 山东地区背景应力场反演结果(0.5°网格， 阻尼系数e=0.5)(a) 最大水平主压应力方向； (b) 最大水平主张应力方向； (c) 最大水平主压应力方向误差Figure 6. Inverted stress results for Shandong region with damping parametere=0.5, on a 0.5°×0.5° grid (a) The maximum horizontal principal compressive stress orientation; (b) The maximum horizontal principal extensive stress orientation; (c) Errors in estimated principal stress orientations图 6a给出的从鲁西到黄海海域地区最大水平主应力的方向变化，直观地反映了山东地区背景应力场的作用方式． 从图 6中可以看出，山东地区的最大主应力方向整体上为近东西方向，与前面使用ZMAP得到的结果一致，也与前人(许忠淮， 1983，1989； 王辉等，2004； 张国民等，2004)对华北地区研究得到的结果吻合； 并且最大主应力方向在空间上存在平稳的过渡变化，存在自西向东逐渐的偏转，由ENE向逐渐转化为近EW向，这一现象与汪素云等(1996)以及王绳祖和张流(2002)的研究结果一致．
对比前人对山东地区不同区域的研究结果(魏光兴等，1982； 周翠英，魏光兴，1987； 魏光兴等，1993； 周翠英等，2005)，本文结果与其也非常一致． 周翠英等(2005)根据区域综合机制解研究认为，以鲁中隆起为主的泰沂蒙山区及其周边地区主压应力轴约为60°—80°，主张应力轴约为330°—350°，本文统计结果显示该区的最大主压应力方向变化范围为67.2°—93.5°； 沂沭带以东的胶东隆起和南黄海北部凹陷地区主压应力轴优势方位约为80°—100°，主张应力轴约为350°—360°，本文的统计结果显示该区的最大主压应力方向变化范围为80.0°—112.8°，平均最大主应力方位为89.1°±7.3°． 沂沭带地区最大主压应力轴方位为72°±21°，主张应力轴方位为342°±20°(周翠英，魏光兴，1987)，我们得到这一地区的最大主应力方向为74.0°±5.5°，也与之非常一致．
图 6a显示，鲁豫交界地区主压应力方向差异变化较大，聊考断裂带中段的主压应力方向为WNW，而在其南侧的鲁西南及河南地区主压应力方向显示为ENE. 机制小组(1981)研究也曾发现，豫北地区地震应力主轴方向离散，并且认为造成应力主轴取向摆动较大、 错动类型多变的原因，可能是该区断裂性质、 运动方式复杂． 我们结合地质构造背景分析认为，除此以外，这一现象可能还与该地区的受力状态有关． 鲁西地区在地质构造属于渤海盆地的凹陷区包含了内黄隆起、 东濮凹陷、 开封凹陷、 鲁西隆起等构造单元，地质条件复杂，构造活动明显； 作为整体，鲁西块体一方面受到秦岭—大别构造带的北东向挤压作用，另一方面由于鄂尔多斯地块的旋转平移运动，对华北平原地块也产生了ESE方向的推挤作用．
南黄海北部海域的应力方向偏转较为明显. 其原因为一方面该区震源机制解的个数相对较少且处于研究区的边缘，反演误差较大(图 6c)； 另一方面可能反映了该区真实的复杂应力场(见Wan(2010)文中图 1)．
3.3 山东地区区域应力场与郯庐断裂带的关系探讨
由图 6可以发现，山东地区的构造应力场以郯庐断裂带为界，其东侧胶东半岛及南黄海北部海域地区最大主压应力方向整体呈近EW或ESE方向，而西侧沂沭断裂带及鲁西地区最大主压应力方向整体呈ENE方向． 郯庐断裂带山东段两侧最大主压应力方向的不同，反映了鲁东-黄海地块与鲁西块体受力状态的差异以及运动方式的不同．
作为超壳断裂，郯庐断裂带的断层运动性质直接反映了区域的构造应力状态． 研究发现，现今郯庐带的运动以右旋走滑为主，具有压扭性特征(魏光兴等，1993)． 本文得到的郯庐断裂带山东段两侧最大主压应力方向的不同，恰好与上述研究揭示的郯庐带的运动性质非常吻合． 应力场方向主要受两个因素的影响： 一个是主要的地质单元的构造性质； 另一是与断层的复杂性(Hardebeck，Hauksson，2001)，应力场与断层几何形态之间彼此影响、 相互作用(Townend，Zoback，2004)． 因此我们认为，山东地区郯庐带两侧应力场状态的差异，一方面反映了断裂带两侧地质单元的构造性质的不同，另一方面也与郯庐带的断层形态及其发展演化有关．
4. 讨论与结论
根据1970—2012年5月共387次中小地震的震源机制解资料，使用Hardebeck和Michael(2006)提出的区域应力场反演方法，结合共轭梯度法求解了应力张量的阻尼最小二乘最优化问题，并使用Lund和Townend(2007)推导的利用构造应力张量计算最大水平主应力方向的公式，研究了山东地区的背景构造应力场状态．
计算结果显示，山东地区整体呈现出相对统一的区域应力场. 其最大主应力方向为82°，最小主应力方向为172°； 最大水平主压应力方向主要是ENE—近EW方向，并且在空间上存在平稳的过渡变化，从西往东由ENE方向逐渐偏转为近EW方向. 这一结果与前人(许忠淮， 1983，1989； 汪素云等，1996； 王绳祖，张流，2002； 王辉等，2004； 张国民等，2004)研究结果一致; 一些小的构造单元的结果也与其他研究者(魏光兴等，1982； 周翠英，魏光兴，1987； 魏光兴等，1993； 周翠英等，2005)的结论有较好的吻合． 此外，郯庐断裂带两侧地区最大主压应力方向存在明显不同，反映了不同地质构造块体的不同受力状态．
由于山东地区地震活动水平不高且分布不均，除近几年由初动符号和振幅比得到的机制解资料相对丰富外，早期的机制解个数非常稀疏． 受原始资料的限制，我们对山东地区划分了0.5°×0.5°的网格，部分网格没有小震震源机制资料，个别子区仅有一两个震源机制，在相邻网格的约束下，这些网格也可以得到结果，但结果的可靠性肯定较差． 因此，要了解山东地区更加详细的应力场方面的信息，我们现有的数据质量是远远不够的，这有待未来进一步地研究．
审稿专家对本文的修改提出了中肯的意见和建议； 山东省地震局周翠英研究员提供了大量的原始数据，中国地震局地震预测研究所邵志刚博士提出了富有建设性的建议并给予了相关指导，张浪平博士给予了重要的帮助； 部分图件使用Stefan Wiemer和Max Wyss等开发的ZMAP 6.0生成； 本研究还得益于全国大形势分析预报青年小组成员之间的友好而融洽的讨论. 作者谨在此一并致谢．
-
![]()
图 1 东汉13个州的行政分布图(局部). 图中数字为该州的郡国平均地理面积S0(表1)
Figure 1. Distribution of 13 administrative divisions in Eastern Han Dynasty. The number indicates the geographical mean area S0 of a “Junguo” in this division
![]()
图 2 京师洛阳附近的郡国和大震位置分布概图
Figure 2. Simplified map of “Junguo” and epicenters of strong earthquakes nearby the capital Luoyang
表 1 东汉郡国的有关参数
Table 1 Some parameters for administrative divisions in Eastern Han Dynasty
地区 州名 郡国地理面积 S0/103km2 郡国总数 东汉时期有过地震记载的郡国 中部 司隶 31.418 7 河内,河南,雍州,京兆尹,右扶风,洛阳,长安 冀州 13.159 9 魏郡 东部 豫州 11.418 6 颍川 兖州 11.418 8 陈留,东郡,山阳,东平 青州 11.418 6 东莱,北海，琅琊 徐州 11.418 5 周边 凉州 53.014 12 陇西,汉阳,金城,北地,武威,张掖,武都 荆州 83.395 7 南阳,汉寿 扬州 79.168 6 并州 26.578 9 太原,雁门 益州 61.533 11 益州 注： 郡国平均地理面积的数据引自高建国(2001)，边远地区的交州和幽州未计入． 
下载: 导出CSV
表 2 3个典型震例
Table 2 Three earthquake events as reference
震例 科目 内容分析 震 例 1 引自《目录》 公元前70年6月1日，"河南以东四十九郡国皆震"． 震中山东诸城、 昌乐一带，烈度Ⅸ，震级7 计算 郡国数目N=49,震级M=7． 由式(1)求得地震有感半径R*=400 km，代入式(2)后可得郡国平均有感面积S*=10 254 km2． 因震中位置隶属青州，参照表1的参数,求出郡国平均有感面积与平均地理面积之比S*/S0=0.89 震 例 2 引自《目录》和《汇编》 公元46年10月23日,"郡国四十二地震"． 震中河南南阳，烈度Ⅷ，震级6 1/2 级． 史料原文为"郡国四十二地震，南阳尤甚，坏垣毁屋，压死官民，地震裂" 计算 步骤同震例1． N=42，M=6.5，R*=300 km，S*= 9 243 km2，S*/S0=0.81 震 例 3 引自《汇编》 公元123年，"京都、 郡国三十二地震"，震中位置未定 分析 《后汉书·左雄传》： "(左)雄复上疏谏曰，先帝封野王君，汉阳地震． " 按《资治通鉴》和《后汉纪校注》(周天游，1987)，安帝的乳母王圣被封野王君的时间为"延光 二年四月"． 故地震时间和地点可作为公元123年5月凉州汉阳(今甘肃天水甘谷一带) 计算 郡国数目N=32． 因为汉阳地震京都有感，其震中距可作为有感半径R*=600 km，由式(2)得S*=34 254 km2． 震中为凉州，参照表1的参数，有S*/S0=0.65． 将上述参数代入式(1)后，得M=6.8． 再按《目录》所给烈度震级表(见附录)，可由震级6.8查得对应的震中烈度I0为Ⅷ
下载: 导出CSV
表 3 东汉早期地震汇总
Table 3 Summary of earthquake information in early Eastern Han Dynasty
序号 发震日期和历史事件 史料位置 史料细节 估算震级 震中或震中至 京师的距离 /km 帝 年-月-日 罪己诏/地震诏 郡国 京师 其它 1 光武 46-10-21 罪己诏 郡国42 南阳尤甚， 地裂压杀人 [6 1/2] 河南南阳 2 章 76-05-02 罪己诏 山阳 东平 ≤5 山东山阳 3 和 92-06-09 京师 ≤4 3/4 京师 4 92-08-08 郡国13 5 1/2 >240 5 93-04-07 陇西 陇西 6 95-11-08 京师 ≤4 3/4 京师 7 97-04-08 陇西 陇西 8 105-06-19 雍县 雍县 9 安︱邓太后执政 107 郡国18 5 1/2 —6 >280 10 108 郡国12 5—5 1/2 >230 11 110-01-12 郡国 9 5—5 1/2 >200 12 110-04-14 郡国 4 ≤5 >130 13 110-10-02 云南益州 云南益州 14 111-02-01 郡国10 5—5 1/2 >200 15 113-02-06 郡国18 5 1/2 —6 >280 16 114 郡国15 5 1/2 >250 17 115-12-11 郡国10 5—5 1/2 >200 18 116-03-01— 03-29 郡国10 5—5 1/2 >200 19 117-01-17 郡国 9 5—5 1/2 >200 20 117 郡国13 5 1/2 >250 21 118 郡国14 5 1/2 >250 22 119-03-10 郡国42 京师 地坼裂涌水,坏败 城郭民屋,压人 6 1/2 ~300 河南南阳一带 23 120-01-17 郡国 8 5—5 1/2 南阳余震 24 120 郡国23 5 1/2 —6 南阳余震 25 120—121 灾谴频数，地坼天崩，高岸为谷 南阳余震 26 安 121-10-10 地震诏 郡国35 地坼裂，坏城郭，压杀人 6 1/2 ~300冀南-鲁西一带 27 122-05-23—06-20 京师 ≤4 3/4 冀南-鲁西余震 28 122-08-19 郡国13 京师 5 1/2 冀南—鲁西余震 29 122-10-23 郡国27 5 1/2 —6 冀南—鲁西余震 30 123-05 郡国32 京师 汉阳地震 6 1/2 —7 汉阳 31 124-01-06 京师 ≤4 3/4 京师 32 124 郡国23 京师 5 1/2 —6 ～320 33 125-12-15 郡国16 京师 5 1/2 ～280 34 顺 128-02-22 地震诏 京师 汉阳 汉阳屋坏杀人，地坼水涌出． 汉阳地陷裂 [6 1/2] 汉阳 35 133-06-18 罪己诏/撤高官 京师 地动山崩火灾，二高官以地震免 ≤5 京师 36 134-12-13 撤高官 陇西 二高官以地震免 (6 3/4 —7) 陇西 37 136-02-18 罪己诏/撤高官 京师 地摇京师 ≤5 京师 38 138-02-28 地震诏 京师 金城陇西 地震裂，城郭室屋多坏，坏杀人 [6 3/4] 金城陇西 注： 地震事件根据《汇编》，震级加[]引自《目录》，震级加( )引自冯锐和俞言祥(2006)． 史料中凡"郡国十八地震" 字样，表中简记为"郡国18"，其余类同．
下载: 导出CSV
表 4 震例4和震例5
Table 4 Earthquake events Nos. 4 and 5
震例 科目 内 容 分 析 震例4 引自《汇编》 公元119年3月10日"京都、 郡国四十二地震，或地坼裂，涌水，坏败城郭、 民室屋，压杀人" 计算 按烈度-震级表(见附录)，可得I0为Ⅷ，M=6.5． 此例N=42，由式(1)算出有感半径R* =300 km,有感面积S*=8 997 km2 震例 5 引自《汇编》 公元121年"郡国三十五地震"，"或地坼裂，坏城郭室屋，压杀人" 计算 按烈度-震级表，I0为Ⅷ，M=6.5． 此例N=35，得地震有感半径R*=300 km,有感面积S*= 11 116 km2,推出S0约为13×103km2 分析 参照2.3节的范例参数,该地震在京师洛阳东北300 km一带为合理
下载: 导出CSV
附录 烈度-震级简表
附录 Simplified list of seismic intensity and magnitude
烈度 震级 极震区破坏情况 灾情 最远破坏 建筑物 房屋 山崩 地裂 Ⅶ 5—5 3/4 坏城垛、 城楼、 墙垣多圮 民居多坏、 坏官民庐舍(约<25%) 黄土崖崩， 陡坎有滑坡 河滩等软湿地有裂缝， 间有出水者 有死伤 <30 km Ⅷ 6—6 3/4 墙垣边墙部分崩坏， 多倒塌， 坏沟渠， 桥梁倾等 庙堂仓库等损坏或部分倒塌， 公廨民 房多倾圮(约<50%)， 树木折倒 土岗山脚崩滑， 山石裂缝 平地多裂缝， 涌沙水， 山坡道路间有开裂， 出现新泉， 干枯老泉 人畜多死伤 <100 km 注： 表中内容摘自《中国地震目录》(顾功叙， 1983)的部分烈度项目.
下载: 导出CSV
-
陈达生, 刘汉兴. 1989. 地震烈度椭圆衰减关系[J]. 华北地震科学, 7 (3): 31-42. 冯锐, 武玉霞. 2013. 张衡地动仪的诞生[J]. 中国地震, 29 (2): 178-196. 冯锐, 俞言祥. 2006. 张衡地动仪与公元134年陇西地震[J]. 地震学报, 28 (6): 654-668. 范晔. 1935. 后汉书(顺帝纪)[M]. 上海: 世界书局: 1085-1108. 高建国. 2001. 汉代地震考(上)[J]. 城市与减灾, (5): 20-23. 顾功叙(主编). 1983. 中国地震目录[M]. 北京: 科学出版社: 1-100 钱国祥. 2002. 汉魏洛阳故城沿革与形制演变初探[C]//中国社会科学院考古研究所编. 纪念中国社会科学院考古研究所成立50周年大会暨21世纪中国考古学与世界考古学国际学术研讨会论文集. 北京: 中国社会科学出版社: 46-51. 万年历编写组. 1994. 中华两千年历书(1-2060)[M]. 北京: 气象出版社: 1-80. 汪素云, 时振梁. 1993. 有感半径与震级的关系及其应用[G]//国家地震局震害防御司编. 中国地震区划文集. 北京: 地震出版社: 179-184. 汪素云, 俞言祥, 高阿甲, 阎秀杰. 2000. 中国分区地震动衰减关系的确定[J]. 中国地震, 16 (2): 99-106. 王振铎. 1963a. 张衡候风地动仪的复原研究[J]. 文物, (2): 1-8. 王振铎. 1963b. 张衡候风地动仪的复原研究[J]. 文物, (4): 1-20. 王鹏飞. 2005. 张衡候风地动仪功能测试和感震原理的探讨[J]. 自然科学史研究, 24 (4): 291-318. 谢毓寿, 蔡美彪(主编). 1983. 中国地震历史资料汇编(第一卷)[M]. 北京: 科学出版社: 14-32. 袁道阳, 雷中生, 刘小凤, 刘百篪, 才树华, 郑文俊. 2004. 138年金城-陇西6 3/4 级地震的史料考证与构造背景探讨[J]. 地震地质, 26 (1): 52-60. 周天游(校注). 1987. 后汉纪校注[M]. 袁宏(撰). 天津: 天津古籍出版社: 466-467. 周振鹤. 1998. 中国历代行政区划的变迁[M]. 北京: 商务印书馆: 1-188. -
期刊类型引用(40)
1. 戴宗辉，郑建常，周连庆，曲均浩，王鹏，尹迪. 基于深度学习的山东平原M_S5.5地震余震检测及发震构造研究. 地球物理学报. 2025(03): 925-937 . 
百度学术
2. 彭永明，鞠玮，赵显勇，王洪池，周波，于文浩，于国栋，梁艳，赵习龙. 胶东半岛现今地应力特征与分布预测. 矿产勘查. 2024(04): 499-508 . 百度学术
3. 柴光斌，张辉，崔华伟，李翠芹，赵银刚，池国民，刘海林，李世莹. 潍坊青州ML4.1地震震源区发震构造分析. 震灾防御技术. 2024(01): 61-67 . 百度学术
4. 万永革. 震源机制水平应变花面应变的地震震源机制分类方法及序列震源机制总体特征分析. 地球科学. 2024(07): 2675-2684 . 百度学术
5. 龚萱，杨顺，任长风，Chit Thet Mon，何玉梅，姜明明，艾印双，陈棋福. 缅甸中部印缅山脉-中央盆地的孕震构造环境. 地球物理学报. 2024(09): 3412-3423 . 百度学术
6. 樊文杰. 2021年5月21日漾濞M_S6.4地震及周边的构造应力场特征和动力学意义. 地震地质. 2023(01): 208-230 . 百度学术
7. 仇瀚之，张学涛. 利用HASH方法反演美国三次中等地震的震源机制解. 山东石油化工学院学报. 2023(01): 19-23+41 . 百度学术
8. 李霞，王鹏，张志慧，董腾超. 山东长清4.1级地震序列的遗漏地震检测与发震断层. 大地测量与地球动力学. 2023(12): 1211-1217 . 百度学术
9. YANG Chengwei，WANG Chenghu. Tensile Fractures and in situ Stress Measurement Data Constraints on Cretaceous–Present Tectonic Stress Field Evolution of the Tanlu Fault Zone in Shandong Province, North China Craton. Acta Geologica Sinica(English Edition). 2023(06): 1616-1624 . 必应学术
10. 崔华伟，郑建常，柴光斌，赵金花，王晓山，闫睿. 2020年2月18日济南长清M4.1地震震源区发震构造分析. 地球物理学进展. 2022(01): 1-10 . 百度学术
11. 孙业君，黄耘，刘泽民，郑建常，江昊琳，李婷婷，杨浩，王俊菲. 日本M_W9.0地震前后郯庐断裂带中南段构造应力场动态变化特征. 地球物理学报. 2022(06): 2124-2136 . 百度学术
12. 戴宗辉，李冬梅，王鹏，郑建常，王志才，李霞. 2020年2月18日长清M_S4.1地震发震构造研究. 地震. 2022(01): 111-121 . 百度学术
13. 张正帅，周晨，郑建常，刘承雨. 长岛震群震源谱参数反演和震源参数特征. 地震. 2022(01): 85-98 . 百度学术
14. 孙强，郑建常，王鹏. 利用综合震源机制解法分析鲁中地区构造应力场特征. 华北地震科学. 2021(01): 56-63 . 百度学术
15. 刘同文，王慧敏，吴啸龙，张志进. 华北平原及其临区现今地壳形变特征. 测绘科学技术学报. 2021(02): 117-123 . 百度学术
16. 林眉，郑建常. 长清M4.1地震前区域小震调制比异常分析. 大地测量与地球动力学. 2021(11): 1207-1210 . 百度学术
17. 孙业君，黄耘，刘泽民，郑建常，江昊琳，李婷婷，叶青，方韬. 郯庐断裂带鲁苏皖段及邻区构造应力场特征及其动力学意义. 地震地质. 2021(05): 1188-1207 . 百度学术
18. 王少坡，孙云强，解孟雨，魏东平. 基于震源机制解的鄂霍次克微板块东部俯冲带地区构造应力场反演. 中国地震. 2020(02): 181-199 . 百度学术
19. 王莹，赵韬，刘春，冯希杰. 鄂尔多斯块体西南缘地区震源机制解与构造应力场特征分析. 地震地磁观测与研究. 2020(04): 18-25 . 百度学术
20. 崔华伟，郑建常，张正帅，李冬梅，柴光斌. 长岛地区小地震断层面参数拟合及应力场特征. 地震地质. 2020(06): 1432-1445 . 百度学术
21. 李振月，万永革，盛书中. 米亚罗断裂活动与汶川地震序列活动的关系. 地震地质. 2019(01): 72-83 . 百度学术
22. 崔华伟，万永革，黄骥超，盛书中，靳志同. 帕米尔—兴都库什地区构造应力场反演及拆离板片应力形因子特征研究. 地球物理学报. 2019(05): 1633-1649 . 百度学术
23. 贺振鹏，王泽利，宋矿银，朱华，牛庆明，张书凯. 烟台地下水封洞库工程区地应力场数值模拟. 科学技术与工程. 2019(21): 63-70 . 百度学术
24. 谢张迪，韩竹军. 2014年云南景谷M_S6.6地震发震断层及其动力学参数. 地震地质. 2019(04): 887-912 . 百度学术
25. 杨佳佳，张永庆，谢富仁. 日本海沟俯冲带M_W9.0地震震源区应力场演化分析. 地球物理学报. 2018(04): 1307-1324 . 百度学术
26. 朱成林，甘卫军，贾媛，李杰，殷海涛，孔向阳. 沂沭断裂带两侧地区地震能量释放与块体相对运动的关系. 地震地质. 2018(02): 299-309 . 百度学术
27. 郑建常，徐长朋，李冬梅，周翠英. 利用台站综合节面解研究华北地区应力场的变化. 国际地震动态. 2018(04): 48-57 . 百度学术
28. 刘泽民，张炳，周冬瑞，谢石文，倪红玉，赵朋，朱虹. 东大别地区现代构造应力场空间分布特征研究. 华北地震科学. 2018(02): 23-27 . 百度学术
29. 黄骥超，万永革，梁倩，靳志同. 2017年9月墨西哥M_S8.2地震震源区构造应力场特征. 华北地震科学. 2018(04): 10-15+62 . 百度学术
30. 苗庆杰，高原，刘希强，曲均浩，郑建常. 基于地壳介质各向异性分析乳山地震近场区构造应力特征. 大地测量与地球动力学. 2017(03): 234-239 . 百度学术
31. 丰成君，张鹏，孟静，戚帮申，秦向辉，牛琳琳，孙东生，孙炜峰，谭成轩，陈群策. 郯庐断裂带及邻区深孔地应力测量与地震地质意义. 地球物理学进展. 2017(03): 946-967 . 百度学术
32. 孙业君，赵小艳，黄耘，杨浩，李锋. 云南地区震源机制及应力场特征. 地震地质. 2017(02): 390-407 . 百度学术
33. 李兵，车金虎，马裕凤，张志国. 鲁南某抽水蓄能电站工程闪长岩体内的地应力测量. 地壳构造与地壳应力文集. 2017(00): 185-194 . 百度学术
34. 丰成君，张鹏，戚帮申，孟静，谭成轩，胡道功. 郯庐断裂带附近地壳浅层现今构造应力场. 现代地质. 2017(01): 46-70 . 百度学术
35. 范建柯，丁志峰，徐小明，董冬冬，苏道磊. 基于自适应量子遗传算法的地震重定位及其在山东地区的应用. 地球物理学报. 2016(11): 4075-4088 . 百度学术
36. 苗庆杰，刘希强. 山东乳山震群剪切波分裂参数的时间演化特征分析. 地震学报. 2016(02): 220-231 . 本站查看
37. 张致伟，阮祥，王晓山，王宇航，祁玉萍. 汶川、芦山地震前后四川地区应力场时空演化. 地震. 2015(04): 136-146 . 百度学术
38. 张致伟，周龙泉，程万正，阮祥，梁明剑. 芦山M_w6.6地震序列的震源机制及震源区应力场. 地球科学(中国地质大学学报). 2015(10): 1710-1722 . 百度学术
39. 郑建常，李冬梅，王鹏，赵金花，徐长朋. 2013年莱州M4.6地震序列震源机制与发震构造初探. 地震地质. 2015(02): 384-399 . 百度学术
40. 孙业君，刘红桂，江昊琳，詹小艳，王俊菲，丁烨，叶碧文. 江苏南部地区现今震源机制和应力场特征. 地震研究. 2015(02): 203-210 . 百度学术
其他类型引用(18)
计量
- 文章访问数: 754
- HTML全文浏览量: 296
- PDF下载量: 25
- 被引次数: 58
