1. 引言
世纪以来全球地震灾害频发，其中最为突出的是2008年汶川Ms8.0地震和2015年尼泊尔Ms8.1地震，均造成了重大的人员伤亡和财产损失，其灾情引起了国际社会的广泛关注。一次次的地震巨灾应急救援实际告诉我们，如何科学研判灾情对于合理部署救援力量、最大限度的抢救人民生命具有重要意义。
2008年5月12日的汶川Ms8.0地震，震源深度14 km，死亡69227人，最高烈度Ⅺ度，是中国发生的破坏性最为严重的地震之一；2015年4月25日的尼泊尔Ms8.1地震，震源深度20 km，死亡9084人，最高烈度IX度(图1)。同样，2010年2月27日在智利发生的Mw8.8地震，震源深度33 km，震中地表投影位于智利比奥大区，距离智利第二大城市康塞普西翁91 km，然而地震仅造成795人死亡。这三次8级以上强烈地震均属于潜源地震，震中投影位于人口相对密集区，然而造成的震害损失相差悬殊，主要原因是由于发震构造存在巨大差异，从而对强震震害分布产生重要影响。
地震发生后一段时间为黑箱期 [1]，外界信息非常少，如果忽视强震发生的构造背景而开展快速研判，有可能会造成一些误判。在地震快速评估工作中不仅需要重视震区建筑物抗震性能特征，还应该重视地质构造特征。基于地质构造的震害分布特征研究是地震宏观灾情空间分布研判的基础，从而为抗震救灾指挥提供支撑。因此紧密结合区域构造特征，对震害分布和灾害救援指挥决策具有重要意义。
前人对于地震灾害特征开展了大量研究工作，尤其是近年来愈加重视地震现场调查工作。过去两年发生了数次破坏性地震，学者们开展了大量现场调查工作，针对震害分布特征、机理开展了深入研究。2021年云南漾濞Ms6.4地震位于青藏高原东南缘，区域地质构造环境复杂，导致地震的发生与表现受多种因素共同控制 [2]。2022年门源Ms6.9地震造成的破坏形式以雁列式组合的张裂隙、张剪裂隙、挤压鼓包、断层陡坎等为主 [3]，沿区域内长约22 km的托来山–冷龙岭断裂带附近，道路形成多处裂缝与挤压鼓包，桥梁整体倾斜移位，房屋不同程度破坏 [4]。
本文针对地质构造与强震震害分布特征关系问题，通过对比汶川地震与尼泊尔地震震害分布的差异性，探讨地质构造对震害分布的影响作用，进而为灾情研判、应急救援对策分析提供技术支撑。
2. 构造背景
2.1. 汶川地震构造背景
龙门山推覆构造带由4条断层组成，主要表现为逆冲叠瓦状形式。自西向东依次为汶川–茂汶断裂、映秀–北川断)、灌县–江油断裂和广元–大邑断裂。对于龙门山推覆构造带东西分带的划分，综合不同学者认识 [5] [6] [7] [8]，本文将龙门山褶皱冲断带划分为5个构造带：松潘甘孜褶皱带(汶川–茂县断裂以西)、后山推覆变形带(汶川–茂县断裂与北川–映秀断裂之间)、前山推覆变形带(北川–映秀断裂与灌县–江油断裂之间)、前陆扩展变形带(灌县–江油断裂与广元–大邑断裂之间)和前陆坳陷带(广元–大邑断裂以东) (图2(a))。
龙门山推覆构造带沿南北方向可以分为3段：北段的分界线为北川–安县一线 [9] [10] [11] [12] [13]，南段分界线基本位于绵阳、都江堰地区，具体为卧龙–怀远一线 [10] [13]，中段的典型特征是出露了彭灌杂岩体，以及前缘发育飞来峰 [7] [8]。
Figure 1. Tectonic background map of Tibet Plateau
图1. 青藏高原构造背景图
(a) 为龙门山推覆构造带剖面图(图1中AB位置) (修改自 [16])；A带：龙门山后山推覆变形带，B带：龙门山前山推覆变形带，C带：前陆扩展变形带；(b) 为尼泊尔地区大陆俯冲带地质构造剖面图(图1中CD位置) (修改自 [17])
Figure 2. Geological structure profile
图2. 地质构造剖面图
2.2. 尼泊尔地震构造背景
尼泊尔为南亚内陆山国，位于喜马拉雅山南麓，北邻中国，其余三面与印度接壤。地势北高南低，境内大部分属丘陵地带，海拔1000 m以上的土地占总面积近一半。东、西、北三面多高山，北部为喜马拉雅地区，海拔高度在4877~8844 m之间；中部河谷区，多小山，中部山区占尼泊尔国土面积的68%，主要由的马哈布哈拉山系和相对较矮的丘日山系组成；南部是冲积平原，分布着森林和草原，其中特莱低地占尼泊尔领土面积的17%。
青藏高原南缘发育多条活动断裂，依次为喜马拉雅主前缘逆冲断裂(MFT)、主边界逆冲断裂(MBT)、主中央逆冲断裂(MCT)、主喜马拉雅推覆构造(MHT)和西藏南部拆离断裂系(STDS)。尼泊尔Ms8.1地震发生在主喜马拉雅逆冲断裂(MHT)之上(图2(b))。主喜马拉雅逆冲断裂为多层结构，由北往南逐渐加深 [14]，在主喜马拉雅逆冲断裂长期的挤压运动中，上地壳发生了明显的缩短 [15]。
3. 烈度分布特征
3.1. 汶川地震烈度特征
汶川地震发生后，中国地震局第一时间组织开展现场调查，调查面积达500,000 km²，调查点4150个，编绘了汶川8.0级地震烈度分布图 [18]，即图3中烈度圈范围。
XI度区面积约2419 km2，存在映秀和北川两个极震区；XI度区内的一般民用建筑物几无完者(图4(a)~(c)；Ⅹ度区面积约3144 km2，呈NE向展布；IX度区面积约7738 km2，呈NE向狭长展布；IX度区呈现明显的带状结构，与地表破裂带比较一致；VIII度区面积约27,786 km2，呈NE向不规则椭圆形状展布，SE方向受地形影响不规则衰减；VII度区面积约84,449 km2，与Ⅷ度区形态基本一致；VI度区面积约314,906 km2，呈NE向不均匀椭圆形展布 [18]。
Figure 3. Distribution map of Wenchuan earthquake damage (intensity circle information data [18])
图3. 汶川地震震害分布图(烈度圈信息据 [18])
Figure 4. Building structure type and seismic damage investigation photos
图4. 房屋结构类型及其震害现场调查照片
3.2. 尼泊尔地震烈度特征
尼泊尔地震发生后，中国地震局工作队赴尼泊尔开展了震害调查工作，针对房屋建筑、生命线工程以及滑坡等震害开展调查，修订了“尼泊尔8.1级地震烈度分布图” [19]。
Figure 5. Distribution map of Nepal earthquake damage (intensity circle information data [19])
图5. 尼泊尔地震震害分布图(烈度圈信息据 [19])
尼泊尔地区建筑结构类型主要为片石结构、砖木结构、砖混结构和自建钢筋混凝土框架结构等4种形式。
在8.1级与7.5级地震仪器震中附近出现两个宏观震中(图5)。烈度圈长轴走向北西西。VII度区以上面积约168,290 km2，主要涉及热索瓦、辛都巴尔乔克、廓尔喀以及我国聂拉木、吉隆等16个县。IX度区面积约900 km2。片石结构大部分倒塌(图4(d))；砖混结构部分倒塌，多数墙体剪裂；框架结构个别倒塌，少数框架柱开裂，多数填充墙体开裂(图4(e))。VIII度区面积约2330 km2。片石结构多数倒毁，少数局部墙体倒塌、酥裂；砖混结构少数倒塌，大部分墙体开裂、位错；个别框架结构的柱剪切破坏，少数填充墙局部倒塌，部分开裂(图4(f))。VII度区面积约165,060 km2。少数片石结构房屋倒塌或者存在房屋局部倒塌，部分墙体开裂；个别砖混结构墙体贯通裂缝，少数墙体开裂明显；少数框架结构填充墙开裂、抹灰层脱落。
4. 滑坡分布特征
4.1. 汶川地震滑坡
汶川地震发生在龙门山推覆构造带之上，触发了数以万计的地质灾害，这些地震滑坡危害极大，导致约2万人死亡 [20]。
许冲等 [21] 根据地震前后高分影的解译结果与野外验证工作，研究表明，汶川地震触发了超过197,000处滑坡，覆盖面积约1160 km2，这是目前最详细完整的汶川地震滑坡编录成果。本文采用滑坡数据源自许冲提供的滑坡编录成果。
滑坡体主要集中在龙门山后山推覆变形带(A带)中段，龙门山前山推覆变形带(B带)中段区也存在滑坡集中分布。在松潘–甘孜褶皱带的东南边界滑坡发育相对较多，远离边界地区，滑坡逐渐消失。整个前陆扩展变形区(C带)未发育滑坡。滑坡密度最高值分布在漩口镇、虹口乡和红白镇辖区。滑坡的高密度区域与彭灌杂岩体基本吻合。
滑坡的空间分布受断裂影响非常明显。滑坡沿中央断裂密集分布，随着与断裂距离的增大而急剧减少，沿着后山断裂与前山断裂滑坡的密集程度减弱，中央断裂附近的滑坡发育最为密集。断裂上盘的滑坡数量显著多于下盘。
4.2. 尼泊尔地震滑坡
通过对比尼泊尔震前震后不同时段的高清遥感影像，解译出尼泊尔地震形成的滑坡。图5中棕色区域为解译出地震滑坡数据，总共2072处。学者们在尼泊尔地震后开展了大量研究 [22] [23] [24] [25] [26]，依据不同判定标准解译出的滑坡数量存在一定差异，但是滑坡总体特征基本一致。
沿阿尼科公路与沙拉公路，选择热索瓦县与新都巴尔乔克县数十个典型滑坡地点开展野外调查，验证遥感解译的可靠性，也发现存在一些早期滑坡体在此次地震中重新活动。
尼泊尔地震滑坡密集成带，与地形、地貌、烈度以及余震展布基本一致，沿烈度圈长轴方向展布，发育在高喜马拉雅与低喜马拉雅过渡区域(图5)。滑坡体沿烈度圈长轴方向展布，在2个极震区分布最为密集，而震中连线以南区域地势快速降低，坡度减小，不利于滑坡产生，因此滑坡分布明显减少。
5. 讨论
5.1. 地质构造对震害分布的影响作用
1) 汶川地震震害分布受映秀–北川断裂控制。
汶川地震震害分布受龙门山推覆构造带横向分带特征的影响比较明显。地表破裂带主要沿着后山推覆变形带的东边界断裂(映秀–北川断裂)展布，在前山推覆变形带的东边界(灌县–江油断裂)也出现一段地表破裂。滑坡体主要分布在后山推覆变形带内，前山推覆变形带上滑坡分布较多。总体来说，震害主要集中在后山推覆变形带之上，前山推覆变形带上震害分布次之，松潘–甘孜褶皱带与前陆扩展变形区上震害较轻。
汶川地震震害分布与龙门山推覆构造带纵向分段特征具有一致性。震害分布整体上也存在分段特征：龙门山断裂带的中段发育2条地表破裂带，分别位于映秀–北川断裂与灌县–江油断裂之上；北段发育1条地表破裂带，而南段无地表破裂。同样，滑坡分布也具有明显的分段特征：滑坡体主要分布在中段区域，北段相对较少，主要沿地表破裂两侧分布，南段仅在分界线处发育，往西南方向逐渐消失。总体来说，震害分布与龙门山推覆构造带的南北分段特征一致，主要分布在中段区域，北段次之，南段震害最轻。
2) 尼泊尔地震震害分布受地形地貌影响较大。
震源机制解揭示发震断层产状比较平缓，图5中余震数据为USGS提供，由于发震断层为低倾角，断层面的地表投影面积较大，小震总体分布沿断裂带展布，基本位于Ⅶ度区域内，近似矩形展布(图5)，与烈度圈相近，推测断层破裂面积约为100 × 200 km²。地震滑坡密集成带，与地形地貌、地震烈度、余震分布特征一致，Ms8.1与Ms7.5地震震中连线位于高喜马拉雅与低喜马拉雅交界区域，连线以北区域地势落差较大，悬崖峭壁随处可见，均已在野外调查中验证。
5.2. 两次强震震害差异特征
1) 两次强震的烈度范围相差较大。
汶川地震Ⅶ度以上区域面积约125,536 km²，长轴约566 km，短轴约267 km，表现为NE向的不规则椭圆形状。尼泊尔地震烈度Ⅶ度以上区域面积约16,829 km²，长轴约224 km，短轴约105 km。
汶川地震Ⅶ度以上灾区面积明显大于尼泊尔地震，是其7.5倍。汶川地震极震区为Ⅺ度，而尼泊尔地震极震区仅为Ⅸ度，汶川地震极震区的破坏程度明显强于尼泊尔地震。同样，汶川地震死亡人数也是尼泊尔地震的7.6倍。
2) 地表破裂特征截然不同。
地表破裂带一般在极震区出现，是地震震害的重要特征之一。汶川地震的地表破裂带主要沿映秀–北川断裂、灌县–江油断裂展布。映秀–北川断裂上的地表破裂长度约220 km，南部以逆冲运动为主，北部以走滑运动为主；灌县–江油断裂的地表破裂长度达70 km，表现为纯逆冲运动的破裂特征 [27]。这两条地表破裂带之间还存在一条8 km长的小鱼洞破裂带，NW向，逆冲兼左旋走滑性质 [28]。
低倾角的主喜马拉雅逆冲断裂(MHT)被认为尼泊尔地震的发震断层。但是，遥感影像解译结合现场调查，均未发现尼泊尔Ms8.1地震形成地震地表破裂带 [29]。
3) 地震滑坡分布差异显著。
汶川地震中，滑坡沿中央断裂密集分布，随着与断裂距离的增大而急剧减少，沿着后山断裂与前山断裂滑坡的密集程度减弱，中央断裂附近的滑坡发育最为密集。断裂上盘的滑坡数量显著多于下盘。
与汶川推覆构造带上地震滑坡分布特征不同，喜马拉雅俯冲带中发震断层为低倾角，难以产生地表破裂。地震滑坡分布明显受深部近水平的断层破裂面控制，全部分布在断层上盘，呈相对较宽的矩形区域展布 [26]。
这两次地震诱发滑坡数量、面积存在显著差异。汶川地震解译出197,000处滑坡，覆盖面积约1163 km²，尼泊尔地震解译出2072处滑坡，覆盖面积约67 km²，二者面积相差16.3倍。
5.3. 震害差异成因讨论
映秀–北川断裂与地表破裂带走向NE35˚，倾向NW，近地表倾角达75˚ (徐锡伟等，2008)。余震资料揭示龙门山断裂带在地壳上部为高倾角，往下逐渐变缓，呈“犁形”或“铲形” [16]，在地壳深部15~20 km处并入地壳中部滑脱层 [28]。汶川地震断层科学钻探成果 [30] [31] 发现在地下500~1500 m度处的断层面倾角一般为65˚。震源机制反演认为汶川地震破裂面的倾角由地表向深部变缓，在起始破裂位置断层倾角减缓到30˚~40˚ [32] [33]。
Wang和Fialko [34] 利用地表测量数据反演了尼泊尔地震发震断层的几何形态和同震滑动分布，认为发震断层为一条倾向北、倾角7˚的逆冲断层。单新建等 [35] 基于InSAR和GPS观测数据联合反演获得尼泊尔地震的同震形变场主要集中在150 km × 100 km的范围内，断层面倾角在浅部约为7˚，倾角随着深度增加而逐渐变大，在20 km深度时接近12˚。刘刚等 [36] 利用GPS数据和地震波联合反演揭示显示尼泊尔地震发生在喜马拉雅主滑脱断层上，倾角11˚。
由此看来，两次地震发震断层的倾角相差较大，汶川地震表现为逆冲带有走滑分量，发震断层在震源深度倾角为30˚~40˚，往上逐渐增大，在地表破裂带出倾角高达75˚，局部地区近直立。而尼泊尔地震发生在MHT之上，该断层近水平，不同学者计算震源深处的断层倾角为7˚~12˚。高倾角断层易于将能量传导至地表，而低倾角尤其是近水平的断层，能量易于消耗在破裂面，传导至地表后衰减很大。最明显的例子就是汶川地震地表破裂带长度近300 km，而尼泊尔地震地表竟然未出现破裂带。
6. 结论
本文通过对比汶川地震与尼泊尔地震震害分布的差异性，获得了有关地质构造与强震震害分布关系的一些认识：
1) 两次强震震害差异悬殊，在受灾面积与地表破裂长度方面均相差较大。地质构造对震害分布的影响非常明显。汶川地震震害分布主要受映秀–北川断裂控制，呈线状展布，龙门山断裂带的分带性与分段性影响了震害的空间分布；尼泊尔地震震害分布受地形地貌影响较大，震害基本沿高喜马拉雅与低喜马拉雅交界区域展布，呈矩形状。
2) 汶川地震由震源深度至地表，发震断裂倾角由40˚增大到75˚，局部地区近直立。而尼泊尔地震发震断层为近水平的MHT，不同学者计算震源深处倾角为7˚~12˚。高角度断层易于能量传导至地表，而低角度尤其是近水平的断层，能量多消耗在破裂面。两次地震发震断层的倾角相差较大，可能影响震害差异的主要原因之一。
致谢
感谢应急管理部国家自然灾害防治研究院许冲研究员提供汶川地震滑坡数据。
基金项目
国家重点研发计划(2021YFC3002104)；河北省地震动力学重点实验室开放基金项目(FZ212101)。