中国火山灾害区划研究历史回顾与未来展望

许建东，万 园，王新茹，潘 波，于红梅，赵 波，杨文健

1.中国地震局 地质研究所，北京 100029；2.北京大学 地球与空间科学学院，北京 100081

0 引言

火山灾害是人类主要的自然灾害之一，国际减灾十年委员会将火山灾害列为第六位主要自然灾害［1-2］.全球有大约10%的人口生活在火山危险地带及其附近.最近1 000 a 以来就有1 300 多座火山喷发，其中约有5%的喷发活动造成了灾难［3-4］.目前有记载的最为严重的火山灾害是1815 年的印度尼西亚坦博尔火山喷发［5］，直接死亡人数9 万多，但喷发造成全球气候变化导致粮食绝产又饿死数百万人.最近引起全球关注的火山喷发事件是2022 年汤加王国洪阿哈阿帕伊岛Hunga Tonga-Hunga Ha’apai（HTHH）火山发生喷发，喷发柱高出海平面以上18～20 km，并迅速向四周扩散至火口周围240 km 的范围，火山喷发所产生的能量引发了整个太平洋的海啸［6］.近年来，由于地球上人口密度不断增大，大量人口不可避免地居住在土壤肥沃的火山区域内，火山灾害就成为了必须面临的重大危险之一.

面对火山的潜在喷发危险和巨大灾害，如何减轻和防御火山灾害成为国际上研究的重点，而火山灾害区划就是其中一种有效的方法［2-3，7］.火山灾害区划是在对某一特定火山的历史火山灾害研究的基础上，通过数值模拟或经验方法预测未来火山喷发时可能造成的破坏或影响区域，为地方政府的土地使用规划提供参考，也为居民的防灾和避灾提供依据.火山灾害区划包括3 个因素：火山灾害的类型、火山灾害的范围、火山灾害的程度.不同的火山灾害种类构成不同类型的危害，不同规模的灾害所造成的影响范围和破坏程度也各不相同.这些通常应综合考虑，编制形成一张火山灾害区划图，对上述各类灾害覆盖的区域、活动类型和规模以不同颜色的方式综合表达.

长白山天池火山是中国境内规模最大的一座活火山，在公元946—947 年发生了巨大的喷发活动，影响到了北半球的大部分地区［8-10］.另外，长白山火山地下岩浆活跃，2002—2006 年间地下小震频发，地表形变异常，具有再次喷发活动的危险［11］.长白山火山不仅引起了中国政府的高度关注，而且东北亚其他各国也担心再次喷发形成巨大灾害，纷纷开展了针对本国可能遭受灾害的研究，如韩国开展了火山灰的灾害区划和应急方案的制定.本研究综合前期对长白山天池火山灾害区划的研究成果，并以此为基础回顾中国火山灾害区划研究历史，总结不同历史阶段各种不同类型火山灾害区划图的优缺点，并结合目前国际上火山灾害区划的研究现状和编图技术，对我国未来编制具有概率含义的火山灾害区划图的思路提出展望.

1 长白山火山概况

长白山天池火山位于中国吉林省东部的中朝边境上，是长白山火山区内3 座大型破火山口中最年轻的一座，其喷发活动约始于5 Ma，经历了造盾阶段（5～2 Ma）、造锥阶段（2～0.05 Ma）和爆炸式喷发阶段（0.05 Ma 以来）［12］.造盾阶段主要为玄武质岩浆的溢流式喷发，区域内形成了广泛而平坦的熔岩台地；造锥阶段由于岩浆的演化，从玄武质逐渐演化至粗面质，黏度增大，流动性变差，原地堆积逐渐形成了高耸的锥体；而晚期岩浆演化成富含挥发份且易爆炸的碱流质，发生了数次大规模的爆炸喷发活动，将原来高耸的锥体破坏并形成了如今的破火山口，且积水成湖形成了天池.天池目前湖面海拔2 194 m，最大水深373 m，蓄水量20×108m3，成为具有巨大潜在危险的天上悬河［13-15］.

长白山天池火山目前处于爆炸式喷发阶段，该阶段大约始于50 ka 前的天文峰期喷发，其后又发生了气象站期和千年大喷发期两次大规模的喷发活动［16-18］.除此之外，部分学者还认为存在其他喷发活动，如公元1018、1124、1265、1373、1401、1573、1668、1702 和1903年的喷发［14，16］，但这些喷发目前没有发现火山喷发物，可能仅是无喷发的岩浆扰动事件.其中最为著名的是千年大喷发，该次喷发发生于公元946—947 年［19-20］，是一次规模巨大的爆炸式喷发，形成了多种火山灾害，如火山灰灾害，喷发柱进入平流层，火山灰扩散至大半个北半球；火山碎屑流灾害横扫方圆50 km 内的一切，至今仍可见裹挟在碎屑流内的大量碳化木.另外，长白山天池火山还在大约距今8 ka 前发生了一次以溢流为主的气象站期熔岩喷发活动，熔岩流覆盖在锥体北坡上，是天池火山熔岩流灾害的典型代表［21］.对于天文峰期喷发目前研究较少，主要认识有其喷发时间大约距今50 ka，且是形成如今天池破火山口的关键喷发.该期喷发地表出露较少，但基本认同其应该也是类似于千年大喷发的一次大规模爆炸式喷发活动［18］.

在历经20 多年系统性研究与监测工作之后，对于长白山天池火山的喷发历史、喷发过程与深部岩浆机制已经取得了一系列成果和基本认识.长白山天池火山的地球物理探测资料表明，在浅地壳中存在低速异常体，分析认为应是地壳岩浆房.同时，多方资料证实天池火山地下深部仍不断向地壳岩浆房进行物质补充，且两者混合后不断地分异演化.当未来达到一定规模后，具有极高的再次触发喷发的危险性.且当地经过多年的发展，旅游业发达，人口相对密集，一旦发生火山喷发，将产生灾害性的巨大影响，需要给与足够的重视.

2 长白山火山灾害区划历史与现状

2.1 火山灾害的主要类型

火山喷发是岩浆等喷出物在短时间内从火山口向地表的释放.由于岩浆中含大量挥发份，加之上覆岩层的围压，使这些挥发份溶解在岩浆中无法逸出，当岩浆上升靠近地表时，压力减小，挥发份急剧被释放出来，于是形成火山喷发.火山喷发产生的灾害有两大类：一类是由于火山喷发本身造成直接灾害，另一类是由于火山喷发而引起的间接灾害.实际上，在火山喷发时，这两类灾害常常是兼而有之.直接灾害主要有火山空降碎屑物灾害、火山碎屑流灾害和火山熔岩流灾害；火山喷发引起间接灾害包括火山泥石流、山崩、海啸等.

2.1.1 火山空降碎屑物灾害

火山爆发式的喷发向空中释放大量火山碎屑物，因受重力作用而下降，在降落过程中产生冲击、掩埋和燃烧现象，在其覆盖范围内造成建筑物、耕地、农作物、森林等被破坏和人畜伤亡等，这类灾害在世界上平均每年发生一次.最著名的是公元79 年意大利维苏威火山喷发，最先喷出的是高黏度的岩浆，继而涌出灼热的火山碎屑流和碎屑堆积物，使当时繁华的庞贝古城被火山渣所掩埋.其他几个著名的海滨城市，如赫库兰尼姆、斯塔比亚等地也遭到严重破坏.1980 年5 月，美国圣海伦斯火山爆发，在距火山口10 km 范围内，无一林木得以保存.此外，降落的火山灰使很大范围的空气、水质受到污染，特别是硫、盐酸等化合物粘在火山灰上，对动植物极其有害.如1970 年冰岛海克拉火山爆发时，富含氯化物的火山灰降落在牧草上达1～6 mm 厚，致使7 000 头食草的羊群死去.

2.1.2 火山碎屑流灾害

在火山喷发时，由热气托着固体颗粒，以高速密集的射流喷出.山坡等地形对炽热火山碎屑流的行进几乎没有任何摩擦阻力，因此所到之处灾害严重.1902年培雷火山爆发时，充满气体且黏度很高的岩浆猛烈喷发和爆炸，炽热的气体形成飓风，伴随着山崩，以约100 km/h 速度从斜坡冲下来，构成一股由温度高达800 ℃的气体（主要是CO2）、玻璃质和尘土组成的火山碎屑流，顷刻吞没了山下的圣·皮尔市.具有破坏性的火山碎屑流，在世界上平均每百年约发生40 起.由于火山碎屑流的速度快，避离困难，横扫一切，破坏力大，加之温度高易发生重大火灾，所以它是火山喷发灾害严重的致灾因素.

2.1.3 火山熔岩流灾害

熔岩流是地下岩浆喷出地表后的高温流体，其温度一般是在几百至1 200 ℃左右.当岩浆喷溢流动时，因其黏度不同，其流速也不一样，一般为几千米每小时，如不及时躲开也很危险.熔岩流所经之处，能破坏建筑物、森林、耕地等，使之变成不毛之地.熔岩覆盖面积的大小，与熔岩黏度、地形条件等有关.熔岩流进入沼泽和水中，能产生第二次爆炸，在靠近冰川的地方，也有发生洪水泛流的情况.历史上有记载的最大熔岩流是1783 年冰岛拉基火山喷发，从16 km 长的裂隙中同时喷出玄武岩，溢满了附近的冰川河床，并蔓延到周围的村庄和海滨地区，覆盖面积达565 km2.

2.1.4 火山泥石流灾害

火山泥石流是由火山喷发碎屑物与水相混合的流动体，同火山碎屑流具有同样快速和猛烈的破坏作用，而且在火山灾害中发生的几率较高.世界上平均每年发生一次以上.因火山喷发死亡的人数中，至少有10%是由泥石流灾害造成的.如哥伦比亚鲁伊斯火山1985 年11 月13 日的喷发虽然发生于没有人类活动的山区，但火山喷发融化了冰雪，形成的泥石流厚达40 m，携带土壤和岩块，以50 km/h 的速度向山口快速流动.火山喷发2.5 h 后，到达了距喷发口74 km 阿梅罗镇，经过短短的4 h，阿梅罗镇被摧毁，造成了23 000 居民死亡、5 000 人受伤、5 000 幢房屋被毁的严重火山次生灾害.

2.2 长白山火山灾害区划研究历史

2.2.1 早期的火山灾害区划

中国早期的火山灾害区划研究开始于20 世纪90年代，1997 年“中国若干近代活动火山的监测与研究”项目（95211）正式启动，中国地震局首次对长白山、五大连池与腾冲火山进行了较系统的监测与研究.在项目实施过程中，本着“以古论今”基本原则，即假设大体相同的火山区很可能受同类喷发事件的影响，火山的未来行为在事件的类型、频率和量级上与过去的喷发行为是相似的.基于这种历史事件重演原则，以长白山天池火山历史上大规模喷发所产生的灾害范围为参考，对于未来同样规模的火山喷发的可能的灾害范围给出了半定量化的厘定，完成了第一张长白山、腾冲火山灾害区划图［22］.

根据降落堆积物的类型和厚度，把降落灾害程度分为重灾区、中等灾害区、轻灾区（图1）.空降重灾区是指遵循弹道轨迹抛射的火山碎屑降落物以及自火山喷发云中降落的火山碎屑降落物，该区的分布范围以天池火山口为中心半径13 km 的范围.空降中等灾害区是厚度大于10 cm 的碎屑物降落区.因为10 cm 的空降层相当于100～200 kg/m2的载荷，可引起许多平顶或尖顶建筑物坍塌.在火山口外30 km 处，浮岩温度仍可在燃点以上.该区空降碎屑物分布主要受风的影响.天池火山历史上的大喷发受西北风影响，未来的降落灾害在天池火山东南侧.空降轻灾区是厚度小于10 cm 的碎屑物降落区，该区主要是对人类的生活和生产造成影响，部分土质房屋会受到破坏，森林会受到影响.细粒火山尘对引擎、管道的堵塞、粘附作用影响范围会更远，虽然对人类生命不构成严重威胁，但会严重影响人的呼吸系统并会刺激灼伤眼睛.如果喷发时火山灰含有大量的有毒物质（如历史上曾经发生的富F、Cl、S 的喷发），还需要考虑火山灰对于水体、食物链的污染作用.

根据火山碎屑流的搬运机制和堆积物的特征，将其分为碎屑流重灾区和碎屑流中等灾害区（图1）.火山碎屑流重灾区是火山碎屑流的主体部分，以层流搬运，堆积物不显层理.它是一种炽热的气固体混合高度流体化的流动体，具有很高的流速，可达160 m/s，可翻越很高的障碍.这个区是以天池火山口为中心，半径40～60 km 的范围.在这个区带将发生流动冲撞、掩埋、高温灼烧等破坏作用，造成火灾、建筑物、桥梁、植物、水电站、微波站、道路等毁灭性破坏以及人员伤亡等.火山碎屑流中等灾害区是火山碎屑流派生的灰云浪部分，以湍流搬运，堆积物具有交错层理.灰云浪影响范围将会超过60 km 半径范围，可能会达到75～90 km 范围.与远缘相空降型火山灰灾害类型相似，灰云浪对空气与水源的污染作用、引擎堵塞作用是主要的灾害效应.

图1 长白山天池火山灾害区划图（据文献［22］修改）Fig.1 Volcanic hazard map of Changbaishan Tianchi volcano（Modified from reference［22］）

根据火山泥石流的特点、河流的密集程度、原生火山碎屑的分布将火山泥石流灾害分为重灾区和中等灾害区带.重灾区分布岩屑流，即岩屑比例高、岩屑颗粒大，河流密布，是原生火山碎屑堆积如火山碎屑流堆积物的主要分布区.以火山口为中心分布的半径为75 km左右（图1）.在这个区内形成高速冲撞、掩埋作用，造成居民点、建筑物、桥梁、铁路、水电站等毁灭性破坏和人员伤亡.中等灾害区分布着火山泥石流及火山泥石流消退流，即岩屑含量变低、浮岩及泥质含量增高，主要沿三江流域分布，其分布距火山口240 km 以外，在这个区内以掩埋、冲击、堵塞作用为主.造成水库淤积、村庄淹没等.

从总体上看，早期“以古论今”为原则的火山灾害区划图编制，通过研究火山的过去喷发行为，以过去火山喷发行为或最具破坏性的喷发事件为参照系，评估预测未来火山喷发的可能影响范围和灾害程度.虽然以定性、半定量评价方法为主，仅仅给出了不同灾害类型影响的大致宏观范围，缺乏详细和具体的评价结果，但初步提出了火山灾害区划的原则、工作程序、灾害区划工作的内容和方法，为进一步的火山灾害区划研究打下了基础.

2.2.2 晚期的火山灾害区划

火山灾害区划的进一步研究开始于中国地震局地质研究所在实施地震行业科研专项“中国主要活动火山喷发序列研究与灾害预测”（200708027）项目中的实践.这项研究使用了从地质调查和地质年代学中获得的过去喷发数据，来模拟火山空降物、火山碎屑流、熔岩流和火山泥流分布范围，采用了几种被广泛接受的数值模型来模拟不同喷发规模的每种灾害.随后按照风险程度进行了灾害区域划分，绘制了长白山-天池火山的综合灾害风险区图［23］（图2）.

图2 长白山天池火山的综合灾害风险区图Fig.2 Comprehensive volcanic hazard zoning map of Changbaishan Tianchi volcano

在长白山天池火山综合灾害风险图的编制过程中，综合考虑了长白山天池火山的历史灾害情况和相关研究成果，采用数值方法先分别对各个灾害类型进行数值模拟，编制分灾种的灾害区划图，最后汇总各灾种的灾害情况综合编制长白山天池火山灾害区划图.在综合对比各种灾害数值模型的基础上，选择了当前应用相对广泛和相对成熟的模拟方法，具体模拟方法与模拟结果如下.

火山空降灾害的模拟计算分为抛射部分和扩散部分.抛射部分的计算使用Charles（2001）编写的Eject模型，设定抛射初始速度为300 m/s，在不同角度下计算获得抛射的最远距离为7.2 km，即火口周围7.2 km应为极危险区域.扩散模型模拟参数的选择主要依据：1）野外调查千年大喷发期和天文峰期两次喷发空降浮岩的分布、粒径、颗粒密度和形状参数等.2）长白山区域内气象台1980—2009 年的风向资料与不同高度风速变化资料［24］.在此基础上，将喷发柱高度作为衡量喷发规模的一个参数，分别设定了喷发柱高度10、20、30 km 进行空降灾害的模拟，发现随着喷发规模的增大，影响范围也不断扩大.喷发柱为10 km 时，火山灰1 cm 厚度线顺风方向最远到达约200 km；而喷发柱为30 km 时，火山灰1 cm 厚度线顺风方向最远到达约450 km.

火山碎屑流模拟计算使用了Flow3D 模型.模拟中设定的物理参数主要来源于千年大喷发火山碎屑流的研究数据，如碎屑流密度790 kg/m3，单碎屑体积0.001 26 m3.另外，碎屑流是受地形强控制的流体，因此模拟中根据长白山火山区DEM 资料选择了12 条沟谷作为碎屑流流动的通道，并且测量了这12 条沟谷的地形高程和坡度.模拟以喷发柱高度来控制喷发的规模，设定了10、20、30 km 喷发规模，且假设喷发柱自天池火口上方坍塌，仅沿设定的11 条沟谷流动扩散.模拟计算结果显示，10 km 喷发柱形成的碎屑流的流动长度介于10.1～13.7 km 之间；20 km 喷发柱形成的碎屑流的流动长度介于31.8～35.4 km 之间；30 km喷发柱形成的碎屑流的流动长度介于50.3～57.8 km之间.研究发现3 个喷发规模下沿二道白河沟谷里流动的距离最长［25］.

火山熔岩流的模拟使用了热流变运动学FLOWGO 模型.该模型的参数设置包括熔岩参数、环境参数和地形参数3 个部分.熔岩参数的选取参照了气象站期熔岩流的研究数据，如熔岩密度2 380 kg/m3，溢出口温度1 328 K.环境参数的选取则根据区域气象台的历史数据，如大气温度使用长白山天池年平均气温-5 ℃.地形参数的设定与碎屑流模拟相似，选取设定了15 条熔岩流流动路径，并测绘了路径的高程和坡度.模拟计算中设定以火口周缘为熔岩溢出口或寄生火口发育点，利用溢出率来控制熔岩的喷发规模，分别设定了100、200、350、500 m3/s 四个喷发等级，计算出了不同规模喷发沿不同路径熔岩流流动所能达到的最远距离［26］.

火山泥石流的模拟使用了基于半经验模型的LAHARZ 软件.该模型的参数包括能量锥阈值、爆发总体积和地形参数3 个部分.能量锥阈值依据千年大喷发喷发柱高度和区域平均地形坡度计算得0.07.地形参数使用区域的DEM 高程资料，在LAHARZ 软件自动计算出的河流中选取4 条火山泥石流易发河道作为模拟路径.爆发总体积成为控制火山泥石流规模的主要参数，分别设定了108、109、1 010、1 011 m3四个喷发规模.依据上述参数开展了模拟计算，得出了不同喷发规模条件下泥石流沿各河道流动的最远距离［27］.

基于各灾害类型的特征和数值模拟方法，分别完成了各灾种的灾害区划，并绘制了相应的灾害区划图，划定了不同规模灾害的影响范围，这为单一灾种的防御提供了参考.然而，火山的喷发活动是复杂的，通常一次喷发可能产生多种灾害类型，而同一区域也可能遭受多种灾害类型的破坏和影响，因此建立综合性的火山灾害区划图对于火山区地方政府和民众防御火山灾害会更加有效和实用.

最终完成的长白山天池火山灾害区划图，综合不同灾害类型的具体情况，围绕长白山天池火山口划分了3 个等级的风险灾害区域（图2）.高风险区位于近火口的锥体之上，除火山泥石流外，其余灾种在该区均可造成严重的破坏，是重点防御区域.中等风险区位于距火口15～60 km 的区域，主要为火山碎屑流的影响区域，小规模喷发时此区相对安全，但中大规模喷发时仍将遭受严重破坏.低风险区域位于距火口60～100 km 的河谷内，主要为火山碎屑流和火山泥石流的影响范围.

高风险灾害区：整体呈面状分布在距火口约15 km的范围内，主要为火山熔岩流、抛射碎屑、顺风向空降浮岩和火山碎屑流所形成的灾害区域.由于该区距离火口近，喷发物温度高，密度大，且速度快，对人类生命财产可造成巨大影响，属于高风险区域.但此区域范围整体有限，加之良好的火山监测可提前预警火山的喷发活动，使人员尽早撤离避免伤亡.

中等风险灾害区：分布在距火口15～60 km 的相对低洼区域，主要受火山碎屑流、火山泥石流和空降浮岩的破坏和影响，属于灾害泛滥区域.该区随着搬运距离的增大，火山喷发物的温度降低、体积变小，势能减弱，风险性有所降低，但对于中等规模以上的喷发活动，破坏性还是巨大，因此仍需做好财产的规避和人员的疏散.

低风险灾害区：分布在距火口60～100 km 的山麓地区和低洼河谷内，主要受火山碎屑流和火山泥石流的破坏和影响.该区远离火口，火山熔岩流和碎屑流物质难以到达该区，但在湖水或降雨的作用下产生火山泥石流，沿河道及两侧造成破坏，该区人员需提高警惕，当发出火山泥石流预警之时尽快撤离，避免不必要的人员伤亡.

长白山天池火山的综合灾害风险区图，充分考虑了长白山火山不同历史灾害类型的具体情况，通过利用新的火山灾害模拟技术，给出了天池火山爆发时在火口周缘所形成的灾害情况和危险程度.这些工作将为地方政府建立火山灾害应急方案和规划地方建设提供科学依据，也为提高火山区周围民众的防灾意识和疏散避险能力提供帮助.

3 结论与讨论

中国的火山监测与灾害评估工作起步较晚，由于缺少现代火山喷发，我国火山研究者没有直接观测火山喷发机会和实施火山减灾的经验，除了研究历史喷发所遗留的喷发物外，我国早期火山灾害区划的相关研究主要集中在野外调查和各种喷发物的鉴定、分期、灾害类型识别等方面，火山灾害区划图以定性和半定量化的表述为主.

近些年来开展的火山灾害评估与灾害区划图编制等尝试性的研究，通过开展中国大陆以长白山火山为代表的不同类型火山未来的火山作用特征、可能的灾害类型和发育程度的研究，参考国际流行的火山灾害模拟技术，建立了反映中国火山活动特点的空降物扩散模型、碎屑流流动模型、火山泥石流模型和火山熔岩流模型，并以此为基础编制了长白山火山灾害区划图，实现了中国火山灾害区划研究与国际火山学潮流的同步.

目前在国际上有关火山灾害区划图还没有形成完全统一的标准，其技术模型、分析方法和火山灾害风险管理方面也有所差异，但基本的方法与思路已经成形，并逐步走向成熟，比如说对灾害进行量化并利用概率论的方法预测火山灾害的“事件树”的发展与应用［28］.我们要继续借鉴国际上的先进经验，发挥我们的特点与优势，继续开展中国大陆不同活动火山地区的火山灾害评价与编图技术研究，提升我国火山灾害区划的研究水平和国际影响力.