北京周围地区地震的分布特点与地壳屈曲
SOME FEATURES OF EARTHQUAKES DISTRIBUTION AROUND PEKING AND CRUSTAL BUCKLING
-
摘要: 地震活动的规模常具有大尺度的特点,结合全球地震分布具有相当的似对称性,我们认为,地震的成因必然蕴含着地球的整体运动.地震活动具有块体性活动的特点.块体的边界线常与经纬线相交成45°.强震发生前也有这种情况.文中给出了与经纬线相交成45°的一个条块,近年华北的强震基本上都落在这个条块的两条边线上,且这些强震的分布具有较明显的等间隔分布特点,类似于地层中的褶皱所显示的一定的波长.我们分析了有历史记载以来我国全部八级地震,发现每次八级地震前后,都有另一次八级或七级半地震在它的一定距离上发生.这个间距有两个优势分布,一个短距的为480—550公里,另一个长距的为1100公里左右(后者约为前者的一倍).引用了柱壳轴向受力屈曲作为初步探讨.从分析可知,柱壳屈曲的半波长(1/2λ≈1.72$ \sqrt{Rh} $)主要决定于圆柱半径与柱壳厚度:以 R=6000公里、h=15公里估算,得半波长为530公里,与我国强震分布间距较为符合.故我们的初步认识是地震的分布可能受着地壳屈曲的控制.Abstract: The large spatial distribution of the activity of strong earthquakes is a characteristic feature. And it seems that the World-wide seismic zones are somehow arranged in several regularly directions on the spherical surface. We suppose that the cause of earthquakes might be related to the movements of the whole earth.A feature relating to earthquakes activity in blocks, the boundaries of which usually make 45° with the latitude or longitude was found. The same situation happened before strong earthquakes occurred in recent years. A stripe making an angle of 45° with the latitude or longitude has been delineated. Strong earthquakes occurred in recent years are located essentially along the two boundary lines of this stripe, and obviously distribute themselves at equal distance span, which means that the distance span between the strong earthquakes resembles equal wave length during the folding of geological strata. We have studied all of the M≥8 earthquakes occurred in China during history. We found that before or after each earthquake (M≥8), another earthquake (M≥7.5) generally occurred at a definite distance away. The distance span has two predominant distributions, one at a shorter distance of 480—550km, the other at a longer distance of about 1100km (the latter distance is about twice the former).We take for a preliminary discussion the buckling process of a simple cylindrical shell under compression in the axial direction. From analysis we know that, whenthe cylindrical shell buckles, the buckling wave length (1/2λ≈1.72 $ \sqrt{Rh} $) was de-termined only by the radius of the cylinder and the thickness of the shell. Now if we substitute R=6000 km; h=15 km, the halfwave length will be 530km. It agrees with the distance span for strong earthquake distribution.So we have the tentative opinion that the character of earthquakes distribution is controlled by crustal buckling.
-
引言
2011年3月11日日本东北地区太平洋沿岸的日本海沟发生MW9.0地震，该地震引发巨大海啸，造成日本重大的人员伤亡与经济损失.由于海啸设防不足，还导致了福岛核电站损坏、放射性物质外泄的严重次生灾害.这次巨大地震是日本海沟俯冲带上发生的前所未有的超级地震，十分罕见，震后启示日本地震学界重新评价他们一直所关注的西南日本太平洋海域的地震海啸危险性.该海域存在着所谓的“南海海槽”俯冲构造带，历史上8级左右地震反复发生，此前日本一直预报但迟迟未发生的震级可达8级的“东海地震”就位于该构造带东北段.经过对发震能力的重新评价，日本地震调查研究推进本部认为南海海槽如果全部发生联动性破裂，则具备发生罕遇超9级地震的能力，而不仅仅是此前认为的单次地震只会发生部分破裂，即只有8级左右的发震强度，从而修订了其原先的地震海啸危险性评价结果(地震調査研究推進本部事務局，2013).2012年日本内阁府首次向社会公开发布新的研究结果，日本沿海地区政府和企业表示将以此为依据提高防灾等级，重新进行地震海啸设防，制定相应防御对策.日本的研究结果经媒体、网络传播后，其对我国沿海地区可能产生的影响引起我国民众强烈关注，甚至造成恐慌(人民网，2012).据我国古籍记载，历史上日本南海海槽发生的巨大地震海啸对我国都有不同程度的影响，因而我们需要运用现代科学知识评价日本南海罕遇地震海啸对我国沿海造成的危险性.
本文在历史地震海啸影响评估和基于地震构造的罕遇海啸模拟法评估的基础上，综合评估潜在的日本南海海槽地震海啸对我国华东沿海地区的影响程度，并初步进行了其总体影响规模评级.一旦发生这种情况，可以为政府部门应急提供科学依据，为民众提供科学的地震海啸咨询服务.
1. 基于历史记载评估日本南海海槽地震海啸对我国的影响
地震导致海底变形，从海底到海面的海水整体受到扰动以波的形式传播出去即形成海啸.海啸的大小多以海啸规模表示.最常用的表示方法是今村-饭田海啸规模分级表(Iida，1963)(表 1).该表根据海啸到岸波高将海啸大小分为6级，规模级别以“m”表示. 2011年3月11日日本海沟发生MW9.0地震引发的巨大海啸，传到我国大陆沿岸时各海洋站监测到的海啸波幅多数在50 cm以下.其中江苏、上海沿海海洋站监测到约30 cm的海啸波，最大分别为浙江沈家门和石浦海洋站监测到的55 cm和52 cm的海啸波(国家海洋环境预报中心，2011).其对我国无灾害影响，在我国的海啸规模m为-1级(表 1).
表 1 今村-饭田海啸规模分级表Table 1. Imamura-Iida intensity scale of tsunami
据统计，日本南海海槽每隔100—150年左右就会发生8级以上巨大地震，每次地震都引发大海啸，而且影响到我国沿海.我国地方古籍中对此多有记载.1498年7月9日，日本发生大地震(震中推测在日本南海，M≥8.0)，广大地区有感，京都、三河、熊野最为强烈，同日我国江浙多处发生水溢.嘉定六月十一日申刻，邑中河渠池泽以及井泉悉皆震荡，涌高数丈(疑应为尺)，良久乃定(正德《练川备记》卷十)；金山六月十一日江海泖湖水溢(乾隆《金山县志》卷十八)；松江六月十一日海水溢(乾隆《娄县志》卷十五)；上海六月十一日海溢(同治《上海县志》卷三十)；崇明河水高涌数尺，凡城乡池沼、沟浍皆然，良久乃定(光绪《崇明县志》卷5)；常熟水溢池沼，涌起三四尺(崇祯《海虞别乘》册二)；桐乡水溢时无风雨，河水忽涌三尺，池沼亦然(康熙《桐乡县志》卷二).同时记载的尚有苏州、太仓、震泽、宝山、川沙、青浦、嘉兴、石门等十余处(卢振恒，1991；王锋等，2005；刘昌森等，2010).是时各处无风雨，排除了气象原因的可能性，均未构成海啸危害，总体上评估其影响江浙沪的海啸规模m≤0.
1707年10月28日日本南海发生8.4级地震(震后49天富士山火山大规模喷发)，自伊豆半岛至九州太平洋沿岸，包括大波湾、播摩、伊予、防长、八丈岛等地均遭海啸袭击，仅土佐一处就漂没房屋11 170栋，溺毙1 844人，失踪926人，破损、漂失船舶768艘，最大浪高25.7 m，海啸规模m为3.5级(卢振恒，1991).同日浙江吴兴县双林地震水涌；乌青镇河水暴涨；海盐县(治今盐官镇)地震水沸(误记十月十四日)；安徽巢县无风河塘忽然水起大浪，水面宽处波高丈余，窄处亦有三尺(雍正《巢县志》卷二十一)；无为河水斗，无风波浪自起(乾隆《无为州志》卷二)；贵池水沸逾时(康熙《池州府志》卷二十九).上海各县虽无记录亦当受影响之列.不仅地震波及浙江，钱塘江口亦有海啸反映，强度甚低，海啸规模m为-1级(谢毓寿，蔡美彪，1987).
1854年12月23日和24日，日本东海、南海先后发生两次8.4级地震，引发东京以东的房总半岛至九州太平洋沿岸海啸，久礼波高16.1 m，串本15 m，种崎11 m，古座、牟歧、阿波9 m，伊予西海岸3—4 m，室户3.3 m，共漂没房屋15 000余栋，损坏船舶800余艘，溺毙3 000余人，经济损失惨重，海啸规模m为3—4级.海啸甚至波及北美沿岸，对我国也造成一些灾害.例如，江苏如皋海啸，淹毙多人；上海黄浦水沸，有高二三尺者，南至嘉定，北至苏州皆同(同治《上海县志》卷三十)；松江十一月初四日午后，水忽涌起二尺者三次，船之泊者知之，木牌震响有声，禾郡更甚(1962年重印《葵丑纪闻录》页60)；吴兴双林河水忽涨一二尺，动荡不定，沟池皆然，周围数百里皆同(民国《双林镇志》卷九)；海宁河水无风自涌，如潮涨落，池沼皆然；慈溪河水骤腾三四尺；鄞县河水骤涨三四尺，城中湖无风而起巨浪，观者如堵(咸丰《鄞县志》卷九)；镇海河水骤涨三四尺，状如沸汤(光绪《镇海县志》卷三十七)；临海海潮泛滥，城乡沟池积潦，同时俱沸，历二时止，沿海庐舍多被淹没(咸丰《临海县志稿》卷三十八)(王锋等，2005；刘昌森等，2010).这次日本地震海啸对我国华东区域的影响最为明显，陆上河湖水涌振荡，沿海海溢，应是地震与海啸先后作用的结果，总体上看海啸规模m可达1级.从有详细记载的上述3次历史海啸事件的分析，初步认为发生在日本南海海槽的巨大地震海啸影响我国华东沿海区域的最大规模m不超过1级.
2. 基于构造环境与地震活动性的罕遇地震震级的判定
日本是环太平洋火山地震带上典型的“沟-弧-盆”构造体系，列岛东侧有很深的日本海沟，南侧有南海海槽、相模海槽，北、西侧有边缘海(日本海、中国东海).日本列岛被太平洋板块、菲律宾海板块和欧亚板块包围，处于三大板块衔接处，太平洋板块和菲律宾海板块正向日本列岛下面俯冲.日本列岛存在两条长的构造线，即大地沟带和中央构造线.大地沟带自日本海(138°E处)以南北方向横穿本州岛至太平洋，大地沟带以东地区称东北日本，以西地区称西南日本.中央构造线(33°—35°N)以南西西—北东东方向穿过西南日本.因此日本列岛及其周围构造复杂，地震多发，地震数目约占世界地震总数的1/10.根据历史资料统计，该地区8级以上地震平均15年发生1次，7—8级地震平均每年发生1—2次，6—7级地震每年发生15次左右.日本历史地震资料表明，很多大震发生在日本列岛靠太平洋的海底构造——日本海沟和南海海槽上(冈本舜三，1984).
日本南海海槽是从日本静冈县骏河湾至九州以东海面约700 km2、水深4 km以下的海槽，位于菲律宾海板块与亚欧板块交接处，处于该海槽沿线的日本东海—南海地区为地震多发区.南海海槽发生的地震深度大多在30 km以内，大地震每30—50年一次，特征地震为8级，重现期约为90—150年.从地震活动图像看，目前日本南海海槽地区东北段(即Tokai——“东海地区”)是个地震空区(图 1).从特征地震复发间隔来看，南海海槽地区已进入了危险期.目前日本海洋研究机构已在南海海槽布设海底观测网并开展综合大洋钻探进行“南海海槽发震带试验”研究(Nakano et al，2013).
![]()
图 1 (a) 日本西南地区M≥7.0地震震中分布图(宋治平等, 2011); (b)日本南海海槽分段和各段上历史大震年代及破裂范围示意图Koizumi N, Matsumoto N, Itaba S.2007.Integrated observation well network of Geological Survey of Japan, AIST for study of Tonankai and Nankai earthquakes (ppt).Tokyo:Geological Survey of Japan, AIST..历史记载的地震震级为M7.9-8.4, 自1605年后以90-150年间隔复发Figure 1. (a) Distribution of M≥7.0 earthquakes in southwestern of Japan (after Song et al, 2011); (b) Sketch of segments of Nankai Trough and corresponding rupture ranges of historical events①.The historical earthquakes are of M7.9-8.4 with the recurrence interval of 90-150 years after the year 1605日本地震调查研究推进本部对南海海槽的活动性进行了分段，将其分为东海、东南海和南海等3段，西南日本海域历史上大地震都是由该海槽部分段落破裂产生的.日本研究人员根据历史记载和海啸沉积物推断了历次地震事件的时间、强度和位置.结果表明，1498年发生的M 8.6地震破裂在东南海—东海段、1605年发生的M 7.9地震在南海—东南海段，1707年发生的M 8.4地震使得南海海槽3段大部分破裂，1854年东海—东南海段和南海段相隔一日破裂先后发生两次8.4级地震，1944年和1946年在东南海段和南海段上分别发生8.0级和8.1级地震①(卢振恒，1991).如果仅根据历史最大地震震级去判定南海海槽地震的潜在最大震级，将会低估其发震强度.依据地震构造规模和活动性类比原则，南海海槽具备与日本海沟和印尼巽他海沟相似的发震能力.若在一次地震中3段全部联合而发生一次性破裂，其强度将可达到罕遇的9.1级规模(地震調査研究推進本部事務局，2013).
①Koizumi N，Matsumoto N，Itaba S. 2007. Integrated observation well network of Geological Survey of Japan，AIST for study of Tonankai and Nankai earthquakes(ppt). Tokyo: Geological Survey of Japan，AIST.
3. 设定日本南海海槽发生罕遇地震时的海啸模拟
沿日本南海海槽的特大地震的震源机制为低角度逆断层，均由菲律宾海板块向日本岛弧下俯冲引起的.这些反复发生的地震中，对于1944年东南海地震和1946年南海道地震，Kato和Ando(1997)已求出其震源机制解.这两次地震的发震断层长度分别为100 km以上，为低角度逆断层.又根据海沟俯冲板片模型①，设定海啸模拟时发震断层倾角为25°.根据Yao等(2013)的研究，认为浅部地壳(深度15 km以上)的地震变形才能产生海啸，中深部的滑动位错不会产生海啸.因此设定震源深度平均为10 km.考虑2011年日本海沟MW9.0地震反演的断层面滑动量大小分布不均匀，不同研究人员得到的最大滑动量从30—60 m不等(Grilli et al，2012；Kozdon，Dunham，2013；Maeda et al，2013).为简化计算，我们假定计算时地震破裂面上各处滑动量大小相等，均取为30 m.为便于计算罕遇情况下MW9.1地震同震破裂引起的海底变形，将整个南海海槽沿走向分成6个接续的破裂子段(从东北到西南依次为破裂段1—6，图 2a).具体设定的震源破裂参数见表 2.
①Slab models for subduction zones. http://earthquake.usgs.gov/data/slab/.查询日期：2014-03-10.
![]()
图 2 模拟的日本南海海啸传播过程不同时刻示意图Figure 2. Snapshots of simulated tsunami propagation in Nankai Trough, Japan表 2 设定地震的发震断层破裂参数Table 2. Rupture parameters of seismogenic fault for the scenario earthquake
海啸是一种海洋浅水波.本文采用康奈尔大学Wang和Liu(2006)基于非线性浅水波方程开发的COMCOT数值模式进行海啸模拟.国内外研究人员用该模式多次模拟的实际海啸的结果与实测数据都较吻合，证明该模式是一个有效的海啸计算模式. 刘双庆等(2012)采用COMCOT模拟2011年日本东北太平洋海域地震海啸传播情况，其计算的海啸传到我国沿海的波高结果与我国东南沿海7个验潮站监测到的波高数据相符.
本文研究区域为25°—36°N，119°—141°E，空间网格尺度采用规则的1分弧度网格，网格数为1 321×661，时间步长为1 s.计算区域水深及陆地地形数据引自美国国家地球物理数据中心的ETOPO1数据库.该数值模式在陆地边界采用反射边界条件计算，其它在水区的边界采用辐射开边界条件以“透浪”消波.本文模拟了罕遇情况下产生的MW9.1地震海啸在我国海域及邻区传播的过程，海啸波场如图 2所示.
4. 华东沿海可能遭受的海啸危险性评价
我们在江浙沪沿海设置10个虚拟海啸监测站(编号及分布见图 2a中三角形所示). 图 3给出了10个虚拟监测站记录的本次模拟的MW9.1地震海啸时程曲线，各台站时程曲线记录了海啸到达各地的时间和波高变化信息.
![]()
图 3 10个虚拟监测站记录的模拟海啸时程曲线Figure 3. Curves of simulated travel time of tsunami in ten virtual tidal stations由海啸传播过程图和虚拟监测站数据可以看出，MW9.1地震可产生初始波幅约10 m的海啸，向日本和太平洋两个方向传播，西南日本沿岸各地在10分钟后将陆续受到海啸冲击.据日本内阁府最新公布的南海巨大地震海啸模拟结果，日本沿岸到岸波幅最大可达30多米(内閣府政策統括官(防災担当)，2012)，海啸规模m达到4级.我们计算的海啸传到本州、四国、九州沿岸最大波幅数据与上述结果相符.1小时后海啸全面穿越琉球群岛，2小时后抵达东海大陆架边缘.由于海啸传播速度只与海水深度有关，到大陆架后随着水深变浅，传播速度变缓.在宽广平缓的东海大陆架上，海啸传播速度比在太平洋深海中明显变缓，传播6小时后海啸才抵达浙江近海；7小时左右浙江沿岸普遍遭到海啸袭击，波幅大约为2—3 m；8小时后海啸逼近上海沿岸；9小时后海啸进入杭州湾和长江口传播，波幅达2 m左右，受地形影响局部地方波高会爬高至近3 m；11小时后海啸邻近苏北黄海沿岸，波幅已普遍在1 m左右.由于海啸的爬高与海岸附近的水深密切相关，且近岸传播时受浅水作用、港湾效应的强烈影响，因此海啸进入华东沿海各河口海湾内部后波幅会增大.考虑到近岸水深数据的分辨率不够及其计算误差，预计罕遇地震情况下实际传到华东海岸的海啸局地最大波幅可达3 m以上.因此综合考虑近岸各影响因素后认为，华东沿海区域的海啸危险性规模m为1—2级.
5. 结论
日本南海海槽每隔100—150年左右就会发生8级以上巨大地震，每次地震都会引发大海啸，而且影响到我国沿海.根据历史记载评估，日本南海海槽的地震海啸事件对我国江浙沪地区的影响规模m最大可达1级.数值模拟研究表明，日本南海海槽发生罕遇MW9.1地震情况下，波源区长度可达500 km以上，西南日本沿岸将遭受大规模海啸袭击.我国华东沿海即使有东海浅水大陆架摩擦衰减，大部分沿岸地区可能遭受的海啸仍接近2 m.若考虑部分地区的地形影响，则局地最大海啸可达3 m以上.综合评估认为，罕遇地震情况下我国华东沿海区域的海啸危险性规模m可达1—2级.更精确的评价结果有待详细考虑近岸传播的非线性效应，以及使用更精细的水深数据进行嵌套网格精细化模拟研究.此外，如果遇上天文高潮位与风暴潮叠加，可能会造成一定程度的海啸灾害.需要指出的是，本文模拟的南海海槽达到发生9级以上地震的可能性很小，但是今后一旦日本南海海域发生8级左右大地震，有关部门则可参考本文结果对海啸影响程度作出迅速估算，为应急提供科学依据.
感谢审稿专家对本文提出的宝贵意见.
-
[1] K. Mogi, Active periods in the world's chief seismic belts, Tectonophysics, 22, 3/4. 256—282, 1974.
[2] G. Dubourdieu, Seismic map of the world, Collége de France Paris, 1973
[3] 许忠淮等, 京津唐张地区地震应力场的方向特征, 地震学报, 1, 2, 1979.
[4] 中国科学院地球物理研究所, 中国强地震简目, 地图出版社出版, 1976.
[5] M. A. Biot, Theory of folding of stratified viscoelastic media and its implications in tectonics and orogenesis, GSAB, 72, 11, 1595—1620, 1961.
[6] S. Timoshenko, Theory of Elastic Stability, Megraw—Hill book conpany, Inc., 324—491.
[7] Brusk, Don O, & Almroth, Bo O, Buckling .of bars, Plates and shells, 1975. -
期刊类型引用(3)
1. 任智源, 侯京明, 王培涛, 李涛, 原野, 赵联大. 精细化地震海啸波流实时预警系统研究与应用. 海洋学报. 2019(09): 145-155 . 
百度学术
2. 闪迪, 王培涛, 任智源, 孙立宁, 赵联大, 范婷婷. 日本南海海槽地震海啸对我国东海沿岸影响. 海洋与湖沼. 2018(06): 1151-1158 . 百度学术
3. 侯京明, 高义, 王君成, 范婷婷, 闪迪, 王宗辰. 基于遥感的台州市海啸脆弱性评估. 海洋学报. 2016(08): 12-23 . 百度学术
其他类型引用(1)
下载:
计量
- 文章访问数: 1020
- HTML全文浏览量: 23
- PDF下载量: 73
- 被引次数: 4
