龔麗文 李 霓 魏海泉 樊祺誠

1）中國地震局地質(zhì)研究所活動構(gòu)造與火山重點實驗室，北京 100029

2）重慶市地震局，重慶 401147

綜述與評述

淺談2014年全球火山動態(tài)及我國火山災(zāi)害預(yù)防對策＊

龔麗文1，2）李 霓1）魏海泉1）樊祺誠1）

1）中國地震局地質(zhì)研究所活動構(gòu)造與火山重點實驗室，北京 100029

2）重慶市地震局，重慶 401147

由于我國現(xiàn)代火山活動案例較少，主要通過搜集全球活動火山區(qū)的監(jiān)測數(shù)據(jù)，了解全球火山活動狀況及應(yīng)對策略。通過對2014年全球火山活動動態(tài)的總結(jié)，逐漸對全球現(xiàn)今活動火山分布位置、活動規(guī)律和災(zāi)害預(yù)防有了新的認(rèn)識；通過對冰島巴達(dá)本加（Bardarbunga）火山的監(jiān)測和追蹤，總結(jié)了冰島政府對火山災(zāi)害的監(jiān)測手段和應(yīng)對策略，希望彌補我國火山監(jiān)測及應(yīng)急等方面的經(jīng)驗空白；通過對最近發(fā)表的火山監(jiān)測技術(shù)及火山信息提取方法的總結(jié)，了解了各種先進(jìn)技術(shù)在預(yù)防和減輕火山災(zāi)害方面的應(yīng)用。最后，結(jié)合國內(nèi)火山監(jiān)測現(xiàn)狀及相關(guān)技術(shù)的發(fā)展趨勢，提出作者對火山監(jiān)測的思考，并提出幾點建議：一方面，我們要充分利用先進(jìn)的計算機(jī)技術(shù)，將實時的監(jiān)測數(shù)據(jù)進(jìn)行智能化處理，提取出必要的信息；另一方面，我們要利用網(wǎng)絡(luò)優(yōu)勢共享最新監(jiān)測成果和先進(jìn)監(jiān)測技術(shù)，為公眾提供一個交流平臺，使得災(zāi)害信息的發(fā)布更智能、更及時、更直觀和更人性化。

全球火山動態(tài)；火山監(jiān)測；災(zāi)害預(yù)防

引言

近些年來，全球構(gòu)造活動趨于活躍，地震和火山活動頻率不斷上升，造成大量地質(zhì)災(zāi)害和人員傷亡。國際上，大量的火山監(jiān)測研究正在進(jìn)行著［1－3］，許多新的技術(shù)手段也加入進(jìn)來，包括數(shù)字模擬、遙感圖像、火山灰建模、火山前兆地震及地表形變等方面［410］，同時大量與火山相關(guān)的網(wǎng)站和機(jī)構(gòu)（如，USGS、NOAA、WOVO dat、VHUB、VOGRIPA、GVP、IAVCEI等）也在不斷地完善和更新，主要用于觀測和研究全球火山活動規(guī)律。這些都為我們更好地了解火山噴發(fā)過程，評估火山危險區(qū)，預(yù)防火山災(zāi)害提供了理論上的支持。

與全球火山活動相比，我國現(xiàn)代火山噴發(fā)相對較少，但新生代以來噴發(fā)的火山數(shù)目較多，分布范圍廣。張大泉［11］根據(jù)前人研究成果，將我國分為7個火山區(qū)，共1 250座火山，其中活火山5座，休眠火山6座；劉若新等［12］分析了長白山天池、金龍頂子、鏡泊湖、五大連池、騰沖、新疆阿什山及臺灣大屯等火山區(qū)發(fā)生新噴發(fā)的可能性，認(rèn)為這些火山仍然處于活動狀態(tài)；吳建平［13］通過對長白山天池和騰沖兩個地區(qū)的地震、地表形變及水熱活動進(jìn)行研究，認(rèn)為這些地區(qū)存在巖漿系統(tǒng)，但是活動性相對較穩(wěn)定。此后，我國對火山監(jiān)測日漸重視，于2006年建成國家火山監(jiān)測臺網(wǎng)中心，到目前為止，已有6處火山受到監(jiān)測［14］，并主要在火山地質(zhì)學(xué)、年代學(xué)、巖石地球化學(xué)和火山活動性的研究上取得了新的進(jìn)展［15－16］。但對這些地區(qū)火山災(zāi)害監(jiān)測仍然處于起步階段，監(jiān)測成果匱乏，監(jiān)測經(jīng)驗不足，且這些地區(qū)大部分已經(jīng)開發(fā)成旅游城市，人口密集，人類活動較多，而相應(yīng)的火山監(jiān)測設(shè)備較少和災(zāi)害防范意識較弱。因此，有必要引進(jìn)國外火山監(jiān)測經(jīng)驗，學(xué)習(xí)其監(jiān)測流程和監(jiān)測手段，預(yù)防和減輕火山災(zāi)害。本文基于這種需求，通過搜集大量相關(guān)資料，總結(jié)火山監(jiān)測各個方面的進(jìn)展，提出對火山監(jiān)測的看法，希望為我國火山監(jiān)測貢獻(xiàn)一份力量。

1 火山活動概述

火山活動是指與火山噴發(fā)有關(guān)的巖漿活動，它包括巖漿噴出地表、產(chǎn)生爆炸、流出熔巖、噴射氣體、散發(fā)熱量、析離出氣體、水分和噴發(fā)碎屑物質(zhì)等活動。火山活動是一個長期的過程，其位置相對固定，具有復(fù)發(fā)的特征。其產(chǎn)生災(zāi)害的范圍受地表條件和氣候環(huán)境的影響較大，如火山的爆炸性受巖漿的揮發(fā)成分和地下水影響，熔巖流的分布受地勢的影響，火山灰分布主要受風(fēng)的影響。此外，火山噴發(fā)后由于巖漿的抽空和火山碎屑的局部堆積，容易引起崩塌或者泥石流等次生災(zāi)害?；鹕降膰姲l(fā)一般具有較多的火山活動前兆，比如大量的中小地震、地表溫度異常、地表形變異常及火山氣體增加，這一切都有利于火山災(zāi)害的預(yù)測。

1.1 火山分布

據(jù)Siebert等［17］統(tǒng)計，全球活火山主要分布在環(huán)太平洋俯沖邊界上，約占活火山總數(shù)的80%（圖1）。此外，在夏威夷和冰島等洋中脊、歐亞板塊與非洲板塊邊界及東非裂谷附近，也分布著一些活火山，但這些地區(qū)（除了夏威夷和冰島）火山活動的頻率和強(qiáng)度遠(yuǎn)不如環(huán)太平洋俯沖帶。而我國近代具有活動跡象的火山主要集中在長白山地區(qū)、五大連池地區(qū)、內(nèi)蒙東部、騰沖地區(qū)和臺灣地區(qū)，這些地區(qū)近些年來雖然沒有噴發(fā)活動，但是中小地震活動和熱泉活動仍然很明顯，且靠近我國的幾個火山活動國家（如，日本、菲律賓和印度尼西亞）仍在頻繁地發(fā)生火山活動，這就要求我們必須時刻做好應(yīng)對火山噴發(fā)的準(zhǔn)備。

1.2 火山災(zāi)害

由于不同類型火山的噴發(fā)方式、噴發(fā)產(chǎn)物和噴發(fā)規(guī)模都有很大的不同，其產(chǎn)生的災(zāi)害類型也有很大的不同。Myers等［18］通過簡化圖（圖2）直觀地展現(xiàn)出復(fù)合火山可能造成的火山災(zāi)害，經(jīng)總結(jié)可以分為：①熔巖流的破壞作用，占據(jù)農(nóng)田，破壞房屋和道路；②火山泥石流，對地表植被破壞作用較強(qiáng)；③崩塌，造成錐體損壞；④熔巖穹丘，可能是新的溢出口；⑤火山彈，爆炸噴射出大量火山彈；⑥火山灰柱，噴發(fā)產(chǎn)生火山灰云，影響航空運營，同時含大量火山灰和酸性氣體，形成酸雨，對地表建筑和植被有很大影響；⑦滑坡或者雪崩，大量火山碎屑或者積雪不穩(wěn)定向下滑動；⑧碎屑流，含大量熱的火山碎屑，掩埋附近地表物質(zhì)；⑨蒸汽噴發(fā)，巖漿沿裂隙上升產(chǎn)生的火山蒸汽。

2 火山監(jiān)測動態(tài)

由于我國缺少現(xiàn)代火山噴發(fā)，有關(guān)國內(nèi)火山監(jiān)測數(shù)據(jù)較少，只能通過搜集國際上一些網(wǎng)站和組織發(fā)表的數(shù)據(jù)來研究火山的噴發(fā)規(guī)律。利用全球火山活動動態(tài)數(shù)據(jù)，反映全球構(gòu)造背景下火山的活動特征和發(fā)展趨勢；利用火山頻發(fā)國家（冰島）對活動火山的監(jiān)測數(shù)據(jù)處理及災(zāi)害的應(yīng)對流程，獲得單個火山活動的前兆信息、噴發(fā)規(guī)律、演化過程及相應(yīng)的監(jiān)測手段和經(jīng)驗。2.1 全球火山監(jiān)測動態(tài)

圖1 全球活火山分布及其與構(gòu)造板塊的關(guān)系（改自Siebert等［17］）

圖2 典型復(fù)合火山可能造成的火山災(zāi)害示意圖［18］

本文通過實時收集史密森全球火山計劃（Smithsonian′s Global Volcanism Program）和美國地質(zhì)調(diào)查局火山災(zāi)害計劃（US Geological Survey′s Volcano Hazards Program）及其他相關(guān)火山網(wǎng)站監(jiān)測數(shù)據(jù)，通過數(shù)學(xué)統(tǒng)計、ArcGis投圖和Excel作圖等方法，綜合分析了2014年全球火山活動情況，包括活動火山位置分布、火山發(fā)生頻率和火山活動隨時間的波動情況，總結(jié)一年內(nèi)全球火山活動的整體規(guī)律，為今后火山監(jiān)測提供參考。2.1.1 2014年活動火山位置分布圖

據(jù)統(tǒng)計，2014年共報道有85座活動火山，通過GoogleEarth建立對應(yīng)火山地標(biāo)文件，用Arcgis進(jìn)行矢量化，然后將其疊加在經(jīng)過配準(zhǔn)的世界地形圖之上（圖3）。圖中清晰地顯示出該時間段內(nèi)活動火山分布的準(zhǔn)確位置，主要分布在太平洋板塊邊界處，集中分布在3個地區(qū)：① 北美洲與南美洲連接處；② 太平洋西岸的日本島和俄羅斯海岸；③印度尼西亞群島。其他地區(qū)的火山活動相對分散，和圖1對比，不難發(fā)現(xiàn)，歷史上火山分布較集中的東非裂谷、歐亞板塊與非洲板塊交界處和大西洋中脊的火山，目前都處于平靜期。

2.1.2 2014年活躍火山排名

圖3 2014年全球活動火山分布圖

文中火山的活躍程度以火山持續(xù)活動時間為標(biāo)準(zhǔn)，通過統(tǒng)計火山報道的次數(shù)確定火山持續(xù)時間，當(dāng)報道頻率（該火山的報道次數(shù)/報道總次數(shù)）大于50%，我們就認(rèn)為該火山較活躍，圖3為活躍火山分布位置圖。經(jīng)統(tǒng)計，2014年共有46次（46周）報道，其中報道超過20次的有15座火山，約占總報道火山的1/6，這些火山的報道次數(shù)約占總次數(shù)的1/2（465/808），活躍火山的具體報道次數(shù)如圖4所示，其中基拉韋厄（Kilauea）火山活動每周均有報道。

圖4 2014年全球活躍火山排名

圖5 周報火山數(shù)目隨時間的變化趨勢

圖6 地震頻次隨時間變化趨勢圖

2.1.3 火山活動時間分布特征

通過總結(jié)每月的活動火山數(shù)，可以發(fā)現(xiàn)2014年6月和9月活動的火山數(shù)目分別是35座和36座，其他時間在30座左右。統(tǒng)計每周報道的火山數(shù)目也可以發(fā)現(xiàn)這一規(guī)律（圖5），在6月19日至7月17日和8月14日至9月25日這段時間里，火山周報報道的火山數(shù)目最多，基本上都超過20座，呈波峰狀態(tài)。楊學(xué)祥等［19］認(rèn)為，這種現(xiàn)象有可能與太空中的外界因素有關(guān)，他們認(rèn)為，9月下旬秋分附近太陽又回到赤道上空，太陽潮使地球扁率變?yōu)樽畲?，自轉(zhuǎn)變?yōu)樽盥?，兩極流體向赤道流動，赤道突起增加，兩極收縮下降，而后發(fā)生反向變化，激發(fā)北極地區(qū)火山噴發(fā)。

2.2 活動火山追蹤——冰島巴達(dá)本加火山

巴達(dá)本加火山位于冰島最大的范特納冰原西北角。2014年6月火山開始膨脹，地表明顯抬升，自8月16日開始出現(xiàn)了強(qiáng)烈的火山不穩(wěn)定性：地震活動異常和GPS測量異常，暗示著巴達(dá)本加火山可能噴發(fā)?；鹕綄＜医⒍喾N模型，利用地震數(shù)據(jù)和GPS數(shù)據(jù)恢復(fù)地表以下巖漿活動狀況，結(jié)果顯示巖漿活動較明顯，且可能存在巖漿上隆，火山地質(zhì)學(xué)家利用GPS數(shù)據(jù)推測巖漿的位置、規(guī)模及走向。8月29日00時02分，噴發(fā)發(fā)生于Holuhraun熔巖區(qū)上一個老的火山裂隙，位于Dyngjujokull冰川邊緣以北大約5km。火山噴發(fā)裂隙長約600m，巖漿噴發(fā)率200～1000m3/s，此后就一直保持著裂隙式噴發(fā)，同時伴隨著一些地震活動，具體跟蹤總結(jié)如下。

表1 2014年8月16日至12月9日巴達(dá)本加火山附近各位置3級以上地震統(tǒng)計

2.2.1 地震活動監(jiān)測

在火山噴發(fā)前期，出現(xiàn)了大量火山噴發(fā)前兆，特別是地震發(fā)生的頻率和強(qiáng)度都呈直線上升，8月16日地震活動性明顯增強(qiáng)，每天有地震上千次，震級多在1～2級，最大超過5級，有兩個震群，大量證據(jù)表明，這與地下巖漿的水平向遷移有關(guān)。伴隨著火山的噴發(fā)和巖漿的溢出，地震活動也大幅度減弱，地震總數(shù)從8月29日的1 194次減少為9月9日的229次，此后都低于200次/日（圖6a）?；鹕絿姲l(fā)后，地震活動總體穩(wěn)定，但是，局部也呈周期性波動，每間隔2～3周，地震活動有微弱的增強(qiáng)，這可能與火山活動有關(guān)。在每天發(fā)生的地震總數(shù)中，0～2級地震約占80%，變化幅度較大，2～3級地震變化幅度相對較平緩（圖6b）。

通過總結(jié)3級以上地震發(fā)生位置（表1），不難看出，3級以上地震主要發(fā)生在破火山口附近，占地震總和的90%，特別是4級以上地震，幾乎只在破火山口附近發(fā)生，這說明破火山口附近更不穩(wěn)定，巖漿活動更強(qiáng)。

2.2.2 熔巖流活動跟蹤

圖7 巴達(dá)本加火山熔巖流面積隨時間動態(tài)變化

2014年8月29日00時38分，衛(wèi)星紅外圖像出現(xiàn)熱異常，這意味著熔巖流上升到地表，此后巖漿一直沿著裂隙溢出，流向東北方向的河流附近，大部分時間溢出速率相對穩(wěn)定，在11月底開始出現(xiàn)脈動式溢出，可以通過衛(wèi)星熱紅外圖像監(jiān)測熔巖流的動態(tài)，包括面積、范圍、溢出速率及產(chǎn)生的火山氣體等其他地質(zhì)災(zāi)害。圖7為熔巖流面積隨時間的動態(tài)變化，從圖中可以看出，熔巖流的面積整體在不斷遞增，遞增速率在不斷減小，這有可能是由于后期巖漿直接覆蓋在前期熔巖流之上造成的。結(jié)合該火山以往的活動規(guī)律（活動時間一般為6～8個月）和熔巖流面積的變化趨勢，可以推出最終熔巖流面積可能接近90km2。

2.2.3 火山灰跟蹤

冰島氣象局利用微分吸收光譜計（DOAS）和傅里葉變換紅外光譜計（FTIR）儀器，估算了巴達(dá)本加火山云中的氣體流量，氣體含量較高的是SO2、CO2、HCl、HF和H2O。噴發(fā)前一個半月平均氣體流量是400kg/s（35kt/d），其峰值為1 300kg/s （112kt/d）。假設(shè)氣體至今仍是連續(xù)性釋放，火山噴發(fā)釋放到大氣中的SO2總量在3.5 Mt（平均流量）～11.2Mt（考慮峰值流量）之間。這是2014年霍魯勞恩火山排放量的5 ～15倍。此外，火山灰云對冰島航空造成的影響較小，其高度大部分集中在1～3km的范圍內(nèi)。

2.2.4 破火山口沉降跟蹤

由于巖漿的抽空或冰原下部融化，造成破火山口附近地表不斷發(fā)生下陷，冰島氣象局主要通過雷達(dá)干涉圖像監(jiān)測火山口附近沉降動態(tài)，包括沉降深度、沉降體積及沉降速度，其沉降速度由一開始的80cm/d下降到后期的25cm/d，監(jiān)測時間范圍內(nèi)沉降深度達(dá)到56m，沉降體積達(dá)到1.7km3，目前沉降仍在繼續(xù)。

圖8 火山災(zāi)害應(yīng)對策略［20］

3 火山災(zāi)害監(jiān)測現(xiàn)狀

Huff等［20］總結(jié)了火山災(zāi)害的應(yīng)急與管理（圖8），認(rèn)為火山噴發(fā)位置一般比較固定，時間比較集中。對火山災(zāi)害的預(yù)防不僅要對火山噴發(fā)后應(yīng)急階段重視，還應(yīng)該提前對火山區(qū)進(jìn)行地質(zhì)調(diào)查，找出潛在噴發(fā)危險區(qū)域，根據(jù)之前噴發(fā)的產(chǎn)物推斷出火山噴發(fā)類型和噴發(fā)規(guī)模，并確定火山復(fù)發(fā)的頻率，通過現(xiàn)有的技術(shù)手段評估火山的危險性，圈定火山產(chǎn)物可能影響的區(qū)域，確定火山周圍土地使用計劃及人員撤離方案，給當(dāng)?shù)厝罕娖占盎鹕奖O(jiān)測知識，最后可生成一個火山災(zāi)害危險等級圖。

隨著科技的發(fā)展，將會有越來越多的監(jiān)測手段應(yīng)用于火山監(jiān)測中。以美國火山監(jiān)測為例，Guffanti等［21］總結(jié)了幾種火山監(jiān)測方法，包括實時地震處理系統(tǒng)、衛(wèi)星遙測技術(shù)、全球定位系統(tǒng)、地電方法和火山氣體監(jiān)測方法，各種地球物理方法和航空技術(shù)在火山災(zāi)害監(jiān)測領(lǐng)域內(nèi)都得到了普遍應(yīng)用。據(jù)美國2008年統(tǒng)計，主要的監(jiān)測手段包括地震監(jiān)測、形變監(jiān)測和水文監(jiān)測（圖9）。具體方法總結(jié)如下。

3.1 地球物理學(xué)方法

近年來，隨著地球物理學(xué)的快速發(fā)展，大量的地球物理方法用于監(jiān)測殼內(nèi)及深部構(gòu)造及巖漿活動，其中在火山監(jiān)測方面應(yīng)用較多的是地震監(jiān)測、重力監(jiān)測和電磁監(jiān)測。

圖9 2008年美國監(jiān)測火山噴發(fā)幾種方法的比例［21］

3.1.1 地震監(jiān)測

地震層析成像可以反演火山下面的三維結(jié)構(gòu)，Kuznetsov等［22］以墨西哥波波卡佩特（Popocatepetl）火山為例，用P波和S波的到時差及vP/vS值建立模型，進(jìn)行層析成像反演殼內(nèi)巖漿房的三維結(jié)構(gòu)。此外，Sicali等［23］用地震時間間隔分布（The space－time inter－event time（IET））分析意大利埃特納火山（Mount Etna）地區(qū)不同深度范圍地殼的地震屬性，推斷出不同構(gòu)造塊的動力環(huán)境。Agus等［24］分析了印尼默拉皮（Merapi）火山2010年大噴發(fā)與地震之間的關(guān)系，詳細(xì)描述了火山從孕育到噴發(fā)各階段不同的地震現(xiàn)象，為以后該地區(qū)預(yù)測大型火山爆炸提供了有用的前兆信息。Inza等［25］利用地震多分量接收數(shù)據(jù)分析了烏比拉斯火山活動的動力學(xué)特征。Dawson等［10］研究了夏威夷基拉韋厄火山噴發(fā)高峰期時長周期的地震特征。3.1.2 重力監(jiān)測

由于不同物質(zhì)成分的密度不同，根據(jù)物理學(xué)知識可以最終推出密度與重力之間的關(guān)系，而通過高精度重力儀測出的重力異常，可以簡單地分析出火山地區(qū)的密度異常，即不同物質(zhì)的分布范圍，Niklas等［26］綜合采用陸地和海洋重力數(shù)據(jù)反演了斯通博利火山殼內(nèi)三維密度，得出垂向剖面和125m深度水平剖面重力異常范圍，并通過野外采樣巖石標(biāo)本的密度特征分析推測出火山地質(zhì)演化信息。

3.1.3 電磁監(jiān)測

Caracciolo等［27］利用高分辨率航磁對埃特納火山進(jìn)行分析，并結(jié)合地貌圖分析該火山的磁異常，認(rèn)為磁異常與該地區(qū)構(gòu)造應(yīng)力場有關(guān)，同時與火山產(chǎn)物分布及輪廓有很大關(guān)系。在國內(nèi)，白登海等［28］利用電磁測深方法探測了熱海地區(qū)下部可能存在低阻體。

3.2 大地測量學(xué)方法

巖漿的上升乃至噴發(fā)伴隨著大幅度的地殼變形，我們可以利用GPS和雷達(dá)干涉圖像對火山進(jìn)行實時形變監(jiān)測，獲得許多火山噴發(fā)的前兆信息，為火山災(zāi)害預(yù)防提供數(shù)據(jù)支持；而火山噴發(fā)后，我們可以通過遙感衛(wèi)星圖像實時監(jiān)測熔巖流和火山灰動態(tài)，預(yù)防次生地質(zhì)災(zāi)害。

3.2.1 形變監(jiān)測

Gerald等［29］利用1998—2005年記錄的地表GPS形變數(shù)據(jù)（基于Mogi模型）反演出菲律賓塔爾火山巖漿房的位置、形狀和深度，并對比了粘彈性殼層模型與等效彈性半空間模型的優(yōu)缺點。除了GPS形變數(shù)據(jù)外，用得比較多的另一種方法是雷達(dá)干涉圖像。巴達(dá)本加火山噴發(fā)后，由于巖漿抽空或冰川底部融化使得地表下沉，通過合成孔徑雷達(dá)圖像實時監(jiān)測了地表沉降幅度、沉降范圍和沉降體積。

3.2.2 遙感監(jiān)測

隨著遙感技術(shù)的不斷發(fā)展，空間分辨率和波譜分辨率不斷提高，遙感技術(shù)逐漸成為實時監(jiān)控地表動態(tài)的重要技術(shù)手段。在災(zāi)害預(yù)測領(lǐng)域，它的應(yīng)用范圍也在逐年增大，主要包括：①火山灰云監(jiān)測，如火山灰柱的高度、漂流方向、體積、覆蓋范圍、成分及是否對航空有影響；②熔巖流監(jiān)測，如熔巖流的溫度、流向、覆蓋面積、體積、溢出速率及未來可能到達(dá)的區(qū)域；③地表形變監(jiān)測，如巖漿上隆造成的形變異常、巖漿房抽空造成地表沉降、沉降速率和沉降體積，這都有利于建模預(yù)測巖漿房的動態(tài)；④地質(zhì)災(zāi)害預(yù)測，如崩塌、雪崩、火山灰泥石流及圈定火山危險區(qū)域范圍等，具體可詳見Oppenheimer［30］的總結(jié)（圖10）。

此外，Christopher［31］通過數(shù)字?jǐn)z影測量和基于GIS分析了日本鹿兒島火山1947—2006年的地貌變化情況，對火山構(gòu)造、地表植被、熔巖流和火山碎屑的分布情況及演化規(guī)律進(jìn)行了研究，為進(jìn)一步判斷火山演化趨勢提供了強(qiáng)有力的證據(jù)。

3.3 地球化學(xué)方法

地球化學(xué)手段應(yīng)用較廣，如利用巖石地化數(shù)據(jù)反演火山形成的巖漿來源、火山成因、巖漿作用、噴發(fā)特征、形成年代及噴發(fā)期次等。通過以上火山形成及噴發(fā)歷史規(guī)律預(yù)測未來火山可能的噴出口位置、噴發(fā)規(guī)模、噴發(fā)方式、噴發(fā)周期性及巖漿演化規(guī)律。

3.3.1 水文監(jiān)測

圖10 遙感在火山監(jiān)測方面的應(yīng)用［30］

水文監(jiān)測包括巖漿揮發(fā)分的監(jiān)測和地下水環(huán)境的監(jiān)測，由于火山噴發(fā)前會有大量火山氣體逸出，地下水的成分和溫度也有很大的異常，因此，水文監(jiān)測可以作為一個很好的火山噴發(fā)前兆監(jiān)測手段。趙慈平等［32］通過對溫泉的基本要素和溫泉水化學(xué)分析數(shù)據(jù)分析，用克里金插值方法獲得了騰沖火山區(qū)的相對地?zé)崽荻鹊钠矫娣植?；趙慈平等［33］還通過對溫泉逸出氣體CO2和CH4碳同位素樣品的采集、分析測試，利用Horita通過實驗矯正的Richet平衡分餾系數(shù)的理論計算數(shù)據(jù)，再通過擬合得到碳同位素平衡分餾方程，反演了巖漿溫度。

3.3.2 噴發(fā)周期性研究

Avellán等［34］通過研究尼加拉瓜尼加帕火山的地層特征、地貌特征及地球化學(xué)特征對火山災(zāi)害的指示意義，推斷該地區(qū)未來有60%的可能性為蒸汽巖漿噴發(fā)，噴發(fā)位置在斷層附近，兩次火山噴發(fā)的間隔周期在400～700年之間。

4 淺談我國火山災(zāi)害預(yù)防

火山活動是由于地下巖漿的運移與噴發(fā)造成的，過去，人們對地下物質(zhì)看不見，摸不著，火山爆發(fā)造成一些毀滅性災(zāi)害和古文明的掩埋。近些年來，隨著科技的高速發(fā)展，監(jiān)測方法和技術(shù)手段越來越多，數(shù)據(jù)精度和監(jiān)測質(zhì)量有很大提高。但是，目前我國在現(xiàn)代火山噴發(fā)，火山監(jiān)測方面經(jīng)驗仍然不足。

4.1 所處現(xiàn)狀

在國家“九五”計劃實施之前，由于對我國火山活動狀況缺乏認(rèn)識，對火山監(jiān)測幾乎是空白，1997年隨著“中國若干近代活動火山的監(jiān)測與研究”項目正式啟動，首次對長白山、五大連池和騰沖火山進(jìn)行了較系統(tǒng)的監(jiān)測與研究。但由于缺乏現(xiàn)代火山噴發(fā)案例，我國火山監(jiān)測研究至今仍處于初級階段［14］。試以同屬地學(xué)的氣象預(yù)報為參照系進(jìn)行比較，氣象預(yù)報已經(jīng)歷了從 “經(jīng)驗性預(yù)報”、“天氣圖預(yù)報”到 “數(shù)值預(yù)報”3個階段。依靠氣象衛(wèi)星、各種高空探測及密集的地面站網(wǎng)等工具，獲得覆蓋廣闊空間的隨時間演進(jìn)的天氣圖和 “鋒面”圖等反映天氣變化基本特征的預(yù)測學(xué)科理論，使氣象預(yù)報由經(jīng)驗性預(yù)報躍進(jìn)到天氣圖預(yù)報，并為進(jìn)入數(shù)值預(yù)報提供了基礎(chǔ)事實和動力學(xué)模型約束。與氣象預(yù)報相對比，可看出火山災(zāi)害預(yù)報當(dāng)前仍處于依靠對時空域殘缺信息的不完全歸納的初級階段。目前，我國火山研究現(xiàn)狀主要是關(guān)于新生代火山錐和火山巖的巖石地化特征及噴發(fā)動力學(xué)模型反演的基礎(chǔ)性研究。

4.2 應(yīng)對策略

由于火山的噴發(fā)具有一定的復(fù)發(fā)性，火山的分布具有很強(qiáng)的地域性，其影響范圍相對集中，其噴發(fā)的前兆較明顯，前兆時間相對較長。加之我們已進(jìn)入地球空間信息學(xué)和數(shù)字地球的時代，多種空間測地、數(shù)字和信息技術(shù)正高速發(fā)展，實時為我們提供過去無法獲得的大空間尺度的、整體動態(tài)的地球變化的多種四維定量信息。因此，我們應(yīng)抓住這一機(jī)遇，充分發(fā)揮好這一優(yōu)勢。經(jīng)總結(jié)，認(rèn)為應(yīng)從以下幾方面著手。

4.2.1 火山災(zāi)害性評估

通過火山基礎(chǔ)性研究，確定出幾個復(fù)發(fā)概率較高的火山和地區(qū)，對其附近的地貌和構(gòu)造進(jìn)行分析，圈定出火山口可能形成的位置及巖漿和火山碎屑可能到達(dá)的危險區(qū)域，人類的活動盡量避開這些危險區(qū)域，從根本上降低火山的災(zāi)害性。

4.2.2 土地利用計劃和建筑選址

對于新城市的規(guī)劃和新建筑的建設(shè)選址一定要避開火山危險區(qū)域和可能造成次生災(zāi)害的區(qū)域。對于已經(jīng)建立在火山附近的城市，遷移難度和成本較大，一定要建設(shè)相應(yīng)的避險裝置，規(guī)劃好避險路線。

4.2.3 建立空、地、深相結(jié)合的立體動態(tài)監(jiān)測系統(tǒng)

建立各種火山噴發(fā)前兆的監(jiān)測臺站，精確地監(jiān)測火山噴發(fā)前兆，如火山地震、地表形變及火山流體等，并把地面觀測臺站逐步與空間測地系統(tǒng)統(tǒng)一在共同的地學(xué)參考框架內(nèi)，這種框架可以是幾何的、物理的，甚至是信息的，但必須是統(tǒng)一的，改變過去分散孤立互無聯(lián)系的狀態(tài)，以便以互補的方式共同為巖石圈中不同空間尺度區(qū)域現(xiàn)今時空演化過程的研究做出貢獻(xiàn)。

4.2.4 前兆信息參量的積累和災(zāi)害識別系統(tǒng)的建立

不同的前兆指示著不同火山噴發(fā)信息，充分利用前兆反演巖漿活動狀況。積累大量前兆信息參量，并建立現(xiàn)場應(yīng)急監(jiān)測前兆臺陣和災(zāi)害識別系統(tǒng)，使得災(zāi)害監(jiān)測信息更全面，災(zāi)害識別更智能。

4.2.5 開發(fā)集成數(shù)據(jù)處理軟件

由于監(jiān)測技術(shù)變得越來越先進(jìn)和復(fù)雜，數(shù)據(jù)處理方法就越來越專業(yè)，所以有必要開發(fā)一個集成數(shù)據(jù)處理軟件，簡化處理步驟，明確處理結(jié)果，對抽象混合的原始數(shù)據(jù)精心處理和轉(zhuǎn)化，經(jīng)若干中心處理，去粗取精、去偽存真，做必要預(yù)處理和規(guī)范化之后，獲得更直觀、更形象的結(jié)果。

4.2.6 開發(fā)數(shù)字模擬實驗室

通過長時間的基礎(chǔ)研究，我們對火山的形成與演化都有了直觀的認(rèn)識。我們可以利用之前的研究成果，結(jié)合現(xiàn)在高速發(fā)展的計算機(jī)模擬技術(shù)，開發(fā)一個數(shù)字模擬實驗室，根據(jù)前兆數(shù)據(jù)進(jìn)行火山的宏觀模擬實驗，彌補傳統(tǒng)小型實驗室無法完成的大型實驗，為更科學(xué)的研究火山提供一個實驗平臺和理論依據(jù)。

4.2.7 建立實時監(jiān)測數(shù)據(jù)庫

數(shù)據(jù)種類繁多，數(shù)據(jù)量大，實時性強(qiáng)，為了減少冗余，且讓各種數(shù)據(jù)得到充分利用，我們可以通過Arcgis或SQL sever建立數(shù)據(jù)庫，共享數(shù)據(jù)，提高數(shù)據(jù)的利用率。

4.2.8 信息資源共享

隨著各種遙感數(shù)據(jù)和航空數(shù)據(jù)精度越來越高，各領(lǐng)域新的研究成果不斷涌現(xiàn)。為了高效利用豐富的資源數(shù)據(jù)，我們應(yīng)該充分應(yīng)用好當(dāng)前的網(wǎng)絡(luò)優(yōu)勢，實現(xiàn)數(shù)據(jù)共享，提高數(shù)據(jù)利用率，加大各學(xué)科之間的聯(lián)系，充分發(fā)揮各學(xué)科優(yōu)勢。

4.2.9 開發(fā)用戶信息交流平臺

監(jiān)測結(jié)果最終面對的是廣大非專業(yè)人士，因此，開發(fā)對應(yīng)的網(wǎng)頁及手機(jī)APP等應(yīng)用軟件，實時為用戶提供共享、交流和查詢平臺，提供短信提醒服務(wù)。用戶可以更清晰地了解災(zāi)害，咨詢更全面的信息，研究人員也可以了解更多的用戶需求，使得產(chǎn)品更人性化。

（作者電子信箱，李霓：lini67＠sina.com）

［1］Cashman K V，Sparks R S J.How volcanoes work：A 25year perspective.Geol.Soc.Am.Bull.，2013，125（5－6）：664－690.doi：10.1130/B30720.1

［2］Hickson C，Spurgeon T，Tilling R，et al.Factors influencing volcanic hazards and the morphology of volcanic landforms.Treatise on Geomorphology，2013，13：219－242

［3］Thouret J C，Németh K.Special issue on volcano geomorphology‘Landforms，processes andhazards’：Introduction.Geomorphology，2012，136（1）：1－5

［4］Alessandro F，Marina B，Patrizia L，et al.A LiDAR survey of Stromboli volcano（Italy）：Digital elevation model－based geomorphology and intensity analysis.Int.J.Remote Sens.，2010，31（12）：3177－3194.doi：10.1080/01431160903154416

［5］Kervyn M，Ernst G G J，Goossens R，et al.Mapping volcano topography with remote sensing：ASTER vs.SRTM.Int.J.Remote Sens.，2008，29（22）：6515－6538.doi：10.1080/01431160802167949

［6］Parrot J F.Tri－dimensional parameterization：An automated treatment to study the evolution of volcanic cones.Géomorphologie，2007，3（3）：247－258

［7］Yosuke A，Teguh P S.Ground deformation associated with the eruption of Lumpur Sidoarjo mud volcano，east Java，Indonesia.J.Volcanol.Geotherm.Res.，2014，278－279：96－102

［8］Moiseenko K B，Malik N A.Estimates of total ash content from 2006and 2009explosion events at Bezymianny volcano with use of a regional atmospheric modeling system.J.Volcanol.Geotherm.Res.，2014，270（2）：53－75

［9］Tony H，Arthur D J，Steven S.Precursory characteristics of the seismicity before the 6August 2012 eruption of Tongariro volcano，North Island，New Zealand.J.Volcanol.Geotherm.Res.，2014，286：294－302

［10］Dawson P，Chouet B.Characterization of very－long－period seismicity accompanying summit activity atKīlauea Volcano，Hawai′i：2007—2013.J.Volcanol.Geotherm.Res.，2014，278－279：59－85

［11］張大泉.中國第四紀(jì)火山.西南師范大學(xué)學(xué)報：自然科學(xué)版，1986（3）：111－119

［12］劉若新，樊祺誠，魏海泉，等.中國活火山研究.地質(zhì)論評，1999，45（增刊）：3－15

［13］吳建平.中國大陸火山監(jiān)測研究進(jìn)展.地震地磁觀測與研究，2005，26（5）：1－10

［14］許建東.中國活動火山監(jiān)測進(jìn)展回顧.礦物巖石地球化學(xué)通報，2011，30（4）：390－392

［15］隋建立，樊祺誠，趙勇偉，等.中國活動火山研究進(jìn)展.礦物巖石地球化學(xué)通報，2011，30（4）：361－364

［16］李霓.淺談我國西部的火山活動及可能的火山災(zāi)害問題.礦物巖石地球化學(xué)通報，2002，21（1）：35－38

［17］Siebert L，Simkin T，Kimberly P.Volcanoes of the World.California：University of California Press，2010，3：551－559

［18］Myers B，Brantley S R，Stauffer P，et al.What are volcano hazards？U.S.Geological Survey Fact Sheet 002－97.http：∥pubs.usgs.gov/fs/fs002－97

［19］楊學(xué)祥，韓延本，陳震，等.強(qiáng)潮汐激發(fā)地震火山活動的新證據(jù).地球物理學(xué)報，2004，47（4）：616－621

［20］Huff W D，Owen L A.Volcanic landforms and hazards.Treatise on Geomorphology，2013，5：148－192

［21］Guffanti M，Diefenbach A K，Ewert J W，et al.Volcano－monitoring instrumentation in the United States.U.S.Geological Survey Open－File Report 2009－1165，32－p.http：//pubs.usgs.gov/of/2009/1165.USGS Open－File Report，2008，1－32

［22］Kuznetsov P K，Koulakov I Y.The three－dimensional structure beneath the Popocatépetl volcano （Mexico）based on local earthquake seismic tomography.J.Volcanol.Geotherm.Res.，2014，276（4）：10－21

［23］Sicali S，Barbano M S，Salvatore D A，et al.Characterization of seismicity at Mt.Etna volcano（Italy）by inter－event time distribution.J.Volcanol.Geotherm.Res.，2014，270（2）：1－9

［24］Agus B S，Philippe L，Sapari D，et al.Analysis of the seismic activity associated with the 2010eruption of Merapi Volcano，Java.J.Volcanol.Geotherm.Res.，2013，261（3）：153－170

［25］Inza L A，Métaxian J P，Mars J I，et al.Analysis of dynamics of vulcanian activity of Ubinas volcano，using multicomponent seismic antennas.J.Volcanol.Geotherm.Res.，2014，270（2）：35－52

［26］Niklas L，Ludovic B，Tullio R，et al.3－D density structure and geological evolution of Stromboli volcano （Aeolian Islands，Italy）inferred from land－based and sea－surface gravity data.J.Volcanol.Geotherm.Res.，2014，273（3）：58－69

［27］Caracciolo F D，Nicolosi I，Carluccio R，et al.High resolution aeromagnetic anomaly map of Mount Etna volcano，Southern Italy.J.Volcanol.Geotherm.Res.，2014，277（3）：36－40

［28］白登海，廖志杰.從MT探測結(jié)果推論騰沖熱海熱田的巖漿熱源.科學(xué)通報，1994，39（4）：344，347

［29］Gerald A G，Andrew V N，Michael W H，et al.Geodetic observations and modeling of time－varying deformation at Taal Volcano，Philippines.J.Volcanol.Geotherm.Res.，2014，271（2）：11－23

［30］Oppenheimer C.Surveillance and mapping of volcanoes and their emissions by satellite remote sensing.Geography，1997，82（4）：317－333

［31］Christopher G.Digital photogrammetry and GIS－based analysis of the bio－geomorphological evolution of Sakurajima Volcano，diachronic analysis from 1947to 2006.J.Volcanol.Geotherm.Res.，2014，280 （2）：1－13

［32］趙慈平，冉華，陳坤華.由相對地?zé)崽荻韧茢嗟尿v沖火山區(qū)現(xiàn)存巖漿囊.巖石學(xué)報，2006，22（6）：1517－1528

［33］趙慈平，冉華，陳坤華.騰沖火山區(qū)殼內(nèi)巖漿囊現(xiàn)今溫度：來自溫泉逸出氣體CO2、CH4間碳同位素分餾的估計.巖石學(xué)報，2011，27（10）：2883－2897

［34］Avellán D R，Macías J L，Pardo N，et al.Stratigraphy，geomorphology，geochemistry and hazard impli－cations of the Nejapa Volcanic Field，western Managua，Nicaragua.J.Volcanol.Geotherm.Res.，2012，213－214（2）：51－71

An elementary introduction to 2014 global volcanic activities and the strategy of volcanic disaster prevention in China

Gong Liwen1，2），Li Ni1），Wei Haiquan1），F(xiàn)an Qicheng1）
1）Key Laboratory of Active Tectonics and Volcano，Institute of Geology，China Earthquake Administration，Beijing 100029，China
2）Earthquake Administration of Chongqing Municipality，Chongqing 401147，China

As modern volcanic activities are rare in China，this paper mainly introduces the global volcanic activities and relevant strategies about disaster prevention by collecting the global volcano monitoring data.Through summarizing the global volcanic activities in 2014，we may have a new understanding about current active volcanoes distribution，patterns of their eruptions and volcanic disaster prevention.We can also get more information about monitoring methods and volcanic disaster mitigation by tracking the monitoring of Bardarbunga taken by Iceland government，hope to fill the blank of our experiences in volcano monitoring and emergency.We will learn some about advanced technologies in volcano monitoring and methods applied in volcanic information extraction based on the latest records.Considering the situation of volcano monitoring in China and the trend of relevant technology application，the authors give some suggestions as follows：on the one hand，we should make full use of advanced computer technology，do some intelligent processing for real time monitoring data and extract necessary information；on the other hand，we should share the latest monitoring records and advanced monitoring technology through the network，and provide a communication platform for the public in order to release the disaster information more intelligent，timely，visualized and humanized.

global volcanic dynamic；volcano monitoring；disaster prevention

P317；

A；

10.3969/j.issn.0235－4975.2015.11.004

2015－03－27；采用日期：2015－05－29。

國家自然科學(xué)基金項目 “云南騰沖大六沖火山機(jī)構(gòu)及火山碎屑巖研究”（No.41472305）和中國地震局監(jiān)測預(yù)報司項目 “全球典型火山事件跟蹤”資助。