于紅梅 許建東 林傳勇

(中國地震局地質(zhì)研究所，北京 100029)

長白山天池火山千年大噴發(fā)空降浮巖碎屑的形貌特征和最終沉降速度

于紅梅 許建東 林傳勇

(中國地震局地質(zhì)研究所，北京 100029)

碎屑的形貌參數(shù)是計(jì)算碎屑最終沉降速度公式中的一個(gè)重要參數(shù)，而以往在計(jì)算中只是把碎屑假設(shè)為球形或橢球形，未對(duì)其進(jìn)行詳細(xì)研究。通過對(duì)長白山天池火山千年大噴發(fā)產(chǎn)生的空降碎屑的形貌分析得到一些新的認(rèn)識(shí):空降碎屑以浮巖為主，浮巖碎屑的形態(tài)不規(guī)則，從等軸狀到拉長狀，從次圓狀到棱角狀都有。隨著浮巖碎屑粒徑的減小，顆粒拉長現(xiàn)象明顯，拉長碎屑的比例從15.02%上升到47.5%。并且，浮巖碎屑的粒度越小，顆粒棱角狀越明顯。通過對(duì)樣品的形貌分析，得到其平均形態(tài)參數(shù)F=0.72。如果假設(shè)長白山天池火山千年大噴發(fā)浮巖碎屑為球形(F=1)時(shí)，計(jì)算的最終沉降速度是浮巖碎屑形態(tài)參數(shù)F=0.72時(shí)最終沉降速度的1.52倍，而當(dāng)假設(shè)顆粒為橢球(F=0.5)時(shí)，計(jì)算的最終沉降速度是浮巖碎屑形態(tài)參數(shù)F=0.72時(shí)最終沉降速度的0.89倍。風(fēng)速為10m/s時(shí)，不同形態(tài)參數(shù)的碎屑的沉積等厚線圖顯示，火山碎屑的擴(kuò)散范圍隨著形態(tài)參數(shù)F值的減小而擴(kuò)大。以上分析充分說明，碎屑的形貌特征對(duì)其最終沉降速度和擴(kuò)散范圍都具有較大的影響，今后在模擬火山碎屑擴(kuò)散和災(zāi)害預(yù)測時(shí)應(yīng)該充分加以考慮。

浮巖碎屑 形貌參數(shù) 最終沉降速度 擴(kuò)散模型 長白山天池火山

0 引言

爆炸式火山噴發(fā)產(chǎn)生的碎屑的擴(kuò)散和堆積與火山噴發(fā)柱(volcanic plume)的動(dòng)力學(xué)密切相關(guān)，而火山噴發(fā)柱的動(dòng)力學(xué)則強(qiáng)烈取決于風(fēng)和噴發(fā)柱的相互作用:當(dāng)噴發(fā)柱上升速率大于風(fēng)速時(shí)形成強(qiáng)噴發(fā)柱，而當(dāng)上升速率小于風(fēng)速時(shí)形成弱噴發(fā)柱(Sparks et al.，1997;Bonadonna et al.，2003)。碎屑的沉降過程主要受風(fēng)的平流和火山碎屑顆粒最終沉降速度的影響。Bonadonna等(1998)指出，碎屑顆粒的沉降狀態(tài)明顯影響火山碎屑堆積的減薄速率。因此，最終沉降速度的正確計(jì)算對(duì)火山碎屑堆積的模擬起著關(guān)鍵作用。Kunii等(1969)、Wilson等(1979)、Suzuki(1983)、Dellino等(2005)等均提出過計(jì)算火山碎屑顆粒最終沉降速度的公式。但是，目前火山碎屑擴(kuò)散范圍的模擬通常采用Suzuki(1983)的二維擴(kuò)散模型(趙誼等，2002;趙誼，2003;于紅梅等，2007)，而該模型應(yīng)用的碎屑最終沉降速度是Suzuki(1983)通過對(duì)實(shí)際測量數(shù)據(jù)的統(tǒng)計(jì)得到的經(jīng)驗(yàn)公式:

式(1)中，V0為碎屑顆粒最終沉降速度，d是近橢球狀碎屑顆粒主軸平均直徑，d=(a+b+c)/3，a、b和c分別為近橢球狀碎屑的長軸、中間軸和短軸，ηa、Ψa分別為空氣的黏度和密度，Ψp為碎屑顆粒的密度，g為重力加速度，F(xiàn)為碎屑顆粒的形態(tài)參數(shù)，F(xiàn)=(b+c)/2a。

由式(1)可見，碎屑顆粒的形態(tài)參數(shù)是計(jì)算碎屑顆粒最終沉降速度的一個(gè)重要參數(shù)。但是，在對(duì)火山碎屑的擴(kuò)散進(jìn)行模擬時(shí)，這些數(shù)據(jù)經(jīng)常缺失，火山碎屑通常被假設(shè)為球形或橢球形(Heffter et al.，1993;Searcy etal.，1998;Bonadonna etal.，2005;Macedonio etal.，2005;Pfeiffer et al.，2005)。因此，開展火山碎屑顆粒的形貌特征分析以獲得準(zhǔn)確的形貌參數(shù)和最終沉降速率，對(duì)火山碎屑擴(kuò)散和堆積的模擬具有十分重要的意義。于紅梅等(2007)也曾對(duì)長白山天池火山千年大噴發(fā)的火山碎屑擴(kuò)散范圍進(jìn)行過模擬，模擬時(shí)同樣將碎屑顆粒假定為橢球形，采用的碎屑形態(tài)參數(shù)F為0.5，未對(duì)碎屑形貌特征進(jìn)行詳細(xì)的研究。為此，本文以長白山天池火山千年大噴發(fā)空降浮巖碎屑為例，探討這些碎屑的形貌特征及其對(duì)最終沉降速度和擴(kuò)散范圍的影響。

1 采樣概況

長白山天池火山坐落于吉林省東部中朝邊境上，是一座具有潛在噴發(fā)危險(xiǎn)性的中央式復(fù)合火山，其全新世噴發(fā)初步分為5000年前、4000年前、公元1199—1200年、公元1668年及1702年等(劉若新，2000)。其中，1199—1200年爆發(fā)時(shí)噴出的火山灰擴(kuò)散到遠(yuǎn)至日本海及日本北部，造成的破壞相當(dāng)嚴(yán)重，被認(rèn)為是近2000年來全球最大規(guī)模的火山噴發(fā)之一(Machida et al.，1983;Machida et al.，1990;劉若新，2000)。由此，該次噴發(fā)也被稱為長白山天池火山千年大噴發(fā)(魏海泉等，2004;史蘭斌等，2005;樊祺誠等，2005;楊清福等，2007)。

圖1 長白山天池火山千年大噴發(fā)空降碎屑采樣位置圖Fig.1 Sampling site of the tephra fallout erupted during Millennium eruption of the Changbaishan Tianchi volcano.

由于受到噴發(fā)時(shí)風(fēng)向的影響，天池火山噴發(fā)形成的空降碎屑堆積多集中分布于火山口以東及東南方向?？紤]到本研究的主要目的是探討千年大噴發(fā)空降碎屑的形貌特征及其最終沉降速度，取樣位置選擇在距離火口東約35km的圓池至釣魚臺(tái)的公路旁(圖1)。此處公路兩側(cè)均覆蓋了約1m厚的灰白色空降浮巖碎屑，碎屑呈正粒序，分選好，較干凈，且出露面積較廣。同時(shí)，根據(jù)劉若新等(2000)和尹金輝等(2005)的14C測年數(shù)據(jù)，該處確實(shí)屬于千年大噴發(fā)的空降碎屑堆積物。

2 碎屑形貌的特征分析

根據(jù)對(duì)樣品的宏觀和鏡下觀測，長白山天池火山千年大噴發(fā)空降碎屑由浮巖、巖屑、晶屑及少量的玻屑等組成。其中，堿流質(zhì)浮巖是最為主要的成分，尤其是2mm以上的碎屑幾乎均為浮巖，而＜2mm的碎屑，浮巖也占20%～30%。因此，本次主要分析浮巖碎屑的形貌特征。首先利用振動(dòng)篩把碎屑篩分為不同的粒級(jí)，并以φ=-log2(d)來表示，d為顆粒的直徑，單位:mm。根據(jù)篩分結(jié)果，碎屑粒度在0.125～20mm之間(-4φ～3φ)。然后，以1φ為粒級(jí)間距對(duì)每一粒級(jí)的碎屑進(jìn)行稱重，得到每一粒級(jí)下的碎屑質(zhì)量分?jǐn)?shù)，并做出粒度分布直方圖 (圖2)。由圖2可見，本樣品碎屑粒徑較大，1mm以上的碎屑顆粒占99%。根據(jù)粒度分析得到該碎屑樣品的粒度標(biāo)準(zhǔn)偏差 σ =0.85，中間粒徑 Mφ=-2.55。

在粒徑＞0.25mm(2φ)的浮巖中選擇1,178個(gè)碎屑進(jìn)行形貌特征分析。由于＜0.25mm的浮巖碎屑不易挑選，而且處理時(shí)容易破碎，所以未對(duì)＜0.25mm粒徑的碎屑進(jìn)行研究。首先對(duì)碎屑顆粒進(jìn)行圖像拍攝，對(duì)于粒徑＞4mm的顆?？刹捎脭?shù)碼相機(jī)直接拍攝，而粒徑＜4mm的顆粒則需用安裝有數(shù)碼相機(jī)的體視顯微鏡來拍攝。這樣做是因?yàn)轶w視鏡的視域有限，對(duì)于較大的顆粒不能完全拍攝，而對(duì)于細(xì)小的顆粒直接利用數(shù)碼相機(jī)拍攝分辨率不足。拍攝時(shí)注意避免顆粒之間的接觸和重疊，并標(biāo)出比例尺。所獲圖像可利用Photoshop軟件進(jìn)行灰度、色階、亮度等處理，使圖像轉(zhuǎn)化為清晰的二元圖像，如圖3所示。注意，處理時(shí)不要把顆粒邊緣信息丟失或修改，盡量保持顆粒的原有形貌。另外，保存二元圖像時(shí)要保存為NIH Image軟件可識(shí)別的格式，例如tif、emf格式等。最后把二元圖像導(dǎo)入可在網(wǎng)上免費(fèi)下載的NIH Image軟件中，利用該軟件獲得二元圖像中顆粒的形態(tài)參數(shù)(Russ，1995;Dellino etal.，1996;Riley etal.，2003)。

圖2 長白山天池火山千年大噴發(fā)空降碎屑粒度分布直方圖Fig.2 Particle size distribution histogram of tephra particles erupted during Millennium eruption of Changbaishan Tianchi volcano.

利用NIH Image圖像分析軟件可以測定碎屑顆粒最小外切矩形的長度、寬度，顆粒最佳擬合橢圓的長軸和短軸，顆粒的面積和周長等。其中，面積為顆粒所包含的所有像素的總和，周長為圍繞顆粒邊界的長度。在此基礎(chǔ)上，便可計(jì)算表征顆粒形貌的幾個(gè)主要參數(shù)(Dellino et al.，1996;Coltelli et al.，2008):長寬比(aspect ratio)(AR)，定義為顆粒最佳擬合橢圓的短軸/最佳擬合橢圓的長軸，它反映顆粒的整體拉長程度，但不能區(qū)分方形和圓形;形狀因子(form factor)(FF)，定義為4×π×面積/周長2，它表征顆粒的不規(guī)則程度或粗糙度(Russ，1995);壓縮度(compactness)(CC)，定義為面積/(最小外切矩形的長度×最小外切矩形的寬度);矩形度(rectangularity)(RT)定義為周長/(2×最小外切矩形的長度+2×最小外切矩形的寬度)。

上述用來表征碎屑顆粒的形貌特征的參數(shù)，僅取決于碎屑顆粒寬度、長度、面積和周長之間的比值，因此它們是無量綱的參數(shù)，不取決于碎屑顆粒的大小。這些無量綱參數(shù)就避免了可能產(chǎn)生的測定誤差，這也正是該方法的優(yōu)點(diǎn)之所在。

圖3 浮巖碎屑形貌圖像Fig.3 Images of pumice pyroclasts.

3 形貌參數(shù)討論

計(jì)算所有碎屑形貌參數(shù)(AR、FF、CC和RT)的平均值、標(biāo)準(zhǔn)差和中位數(shù)，結(jié)果列于表1中，并制作不同粒級(jí)下浮巖碎屑的AR、FF、CC和RT頻率分布圖(圖4)。

由圖4和表1可見，AR值范圍較寬，介于0.29～1之間。隨著粒徑的減小AR的平均值呈下降趨勢，從0.760降為0.660，中位數(shù)從0.765降為0.663?？梢姼r碎屑的粒徑越小，顆粒拉長現(xiàn)象越明顯。根據(jù)Folk(1974)，火山碎屑顆粒可以分為極拉長狀(AR＜0.6)、拉長狀(0.6＜AR ＜0.63)、次拉長狀(0.63 ＜ AR ＜0.66)、中間形狀(0.66 ＜ AR ＜0.69)、次等軸狀(0.69＜AR＜0.72)、等軸狀(0.72＜AR ＜0.75)和極等軸狀(AR ＞0.75)。我們統(tǒng)計(jì)每個(gè)粒級(jí)范圍內(nèi)的AR值，結(jié)果見表2。由表2可見，隨著粒徑的減小，極等軸狀顆粒的比例呈減小趨勢，極拉長狀顆粒的比例呈增加趨勢。對(duì)于粒徑＞8mm(φ＜-3)的浮巖碎屑，極拉長狀碎屑占9.44%，極等軸狀碎屑占55.36%;而粒徑在0.25～0.5mm(1＜ φ ＜2)之間的浮巖碎屑，極拉長狀碎屑占35%，極等軸狀的碎屑只占25%。

圖4 不同浮巖碎屑粒徑下AR、FF、CC和RT頻率分布Fig.4 Frequency distribution of the AR，F(xiàn)F，CC and RT values calculated for different grain-size classes of the pumice pyroclasts.

CC值范圍在0.44～0.84之間。CC值接近1，代表近矩形;接近0.785，代表近橢圓形，并且CC值越小，顆粒表面越不光滑(Dellino et al.，1996)。隨著粒徑的減小CC平均值呈下降趨勢，從0.707降為0.664，中位數(shù)從0.709降為0.654。從圖4也可以看出，隨著粒度減小，直方圖向左偏移，說明碎屑顆粒表面越來越不光滑。

FF值范圍在0.42～0.86之間，隨著粒徑的減小，F(xiàn)F平均值減小。RT值范圍在0.79～1.05之間，隨著粒徑的減小，RT平均值增大。利用FF和RT可以確定顆粒表面的不規(guī)則程度(Coltelli etal.，2008)。Muller(1967)根據(jù)碎屑形狀把碎屑從極圓狀到棱角狀分為5個(gè)等級(jí)，并給出 Russel，Taylor，and Pettijohn對(duì)比圖(圖5)。利用NIH Image軟件，求得圖5中顆粒的形貌參數(shù)FF和RT值 (表3)，并制作FF和RT值與碎屑形狀關(guān)系圖(圖6)。把本次統(tǒng)計(jì)的長白山天池火山空降碎屑的FF和RT值投點(diǎn)到圖6中，可以確定浮巖碎屑的粒度越小，顆粒棱角狀越明顯。由圖6可見，粒徑＞1mm(φ＜0)的浮巖碎屑屬于次圓狀，0.5～1mm(0＜φ＜1)的碎屑屬于次棱角狀，0.25～0.5mm(1＜φ ＜2)的碎屑屬于棱角狀。

根據(jù)以上對(duì)浮巖碎屑形貌的分析可知，隨著粒徑的減小，拉長碎屑的比例呈增加趨勢，從15.02%上升到47.5%。碎屑表面不規(guī)則，尤其是0.5～1mm(0＜φ ＜1)碎屑屬于次棱角狀，＜0.5mm(φ ＞1)的碎屑屬于棱角狀。顯然，在計(jì)算這些拉長的、棱角狀的浮巖碎屑的最終沉降速度時(shí)，需要充分考慮碎屑形貌的影響。

圖 5 Russel，Taylor，and Pettijohn 對(duì)比圖(引自 Muller，1967)Fig.5 Comparative chart of Russel，Taylor，and Pettijohn(after Muller，1967).

4 顆粒形態(tài)對(duì)最終沉降速度和擴(kuò)散范圍的影響

由于我們只是進(jìn)行二維分析，根據(jù)以往文獻(xiàn)(Wilson et al.，1979;Coltelli et al.，2008)，一般假設(shè)碎屑第3個(gè)軸等于二維確定的短軸，即b=c，所以顆粒的形態(tài)參數(shù)F=(b+c)/2a=b/a，即AR值。因此，由本文所分析的浮巖碎屑的平均AR值可以獲得形態(tài)參數(shù)的平均值為F=0.72。然后，應(yīng)用Suzuki經(jīng)驗(yàn)公式(公式1)計(jì)算最終沉降速度。為了對(duì)比碎屑形態(tài)參數(shù)對(duì)其最終沉降速度的影響，同時(shí)還計(jì)算了顆粒為圓球(F=1)和橢球(F=0.5)時(shí)的最終沉降速度。計(jì)算時(shí)浮巖密度取為1.0g/cm3①陳文寄，樊祺誠，尹功明，等，2000，三個(gè)火山區(qū)的噴發(fā)歷史、規(guī)模、物理及化學(xué)過程和構(gòu)造地質(zhì)背景研究，課題(95-11-03)研究報(bào)告。，得到的結(jié)果見圖7。

圖6 碎屑形貌分析獲得的FF與RT圖Fig.6 FF vs RT plot of the values obtained by themorphological analysis.

由圖7可見，碎屑的形貌對(duì)其最終沉降速度影響較大。例如，當(dāng)顆粒直徑為8mm時(shí)，浮巖碎屑形態(tài)參數(shù)F=0.72時(shí)最終沉降速度為8.2m/s，而假設(shè)碎屑為球形(F=1)的最終沉降速度為12.3m/s，假設(shè)碎屑為橢球(F=0.5)的最終沉降速度為7.3m/s;當(dāng)顆粒直徑為0.25mm時(shí)，3個(gè)速度值分別為0.97m/s、1.27m/s和0.86m/s。擬合3組最終沉降速度的數(shù)據(jù)，得到它們的關(guān)系如下

式中，VF=1是浮巖碎屑為球形(F=1)時(shí)的最終沉降速度，VF=0.72是浮巖碎屑形態(tài)參數(shù)為F=0.72時(shí)的最終沉降速度，VF=0.5是浮巖碎屑為橢球形(F=0.5)時(shí)的最終沉降速度。

圖7 不同形態(tài)參數(shù)(F)得到的最終沉降速度Fig.7 Terminal settling velocities for different shape parameters(F).

為了研究碎屑形貌對(duì)最終空降碎屑擴(kuò)散范圍的影響，我們對(duì)碎屑形態(tài)參數(shù)F=1、F=0.72和F=0.5分別進(jìn)行了擴(kuò)散模擬。模擬時(shí)將風(fēng)速設(shè)為10m/s，其它參數(shù)參考劉若新等(1998)，見表4。根據(jù)模擬結(jié)果畫出碎屑沉積厚度等值線圖(圖8)。由圖8中可見，厚度分別在1000cm、100cm、20cm和10cm處，F(xiàn)=0.5時(shí)擴(kuò)散的范圍最廣，F(xiàn)=0.72次之，F(xiàn)=1時(shí)擴(kuò)散范圍最小，說明形態(tài)參數(shù)F值越小擴(kuò)散范圍越遠(yuǎn)。由此也證明了碎屑的形貌特征對(duì)擴(kuò)散范圍具有較大的影響。

應(yīng)當(dāng)指出，上述模擬結(jié)果與長白山天池千年大噴發(fā)空降碎屑的實(shí)測等厚圖有一定的差異，這可能和所采用的參數(shù)有關(guān)。但是，我們這里強(qiáng)調(diào)的是碎屑形貌特征對(duì)擴(kuò)散范圍的影響，盡管模擬結(jié)果有所偏差，仍可明顯看出碎屑顆粒形貌對(duì)擴(kuò)散范圍的較大影響。

圖8 不同碎屑形狀參數(shù)時(shí)沉積厚度等值線對(duì)比圖Fig.8 Isopach map of tephra for different shape parameters.點(diǎn)劃線為F=0.5;實(shí)線為F=0.72;虛線為F=1;單位:cm

5 初步討論

根據(jù)以上分析，本文獲得了以下初步結(jié)論:

(1)長白山天池火山千年大噴發(fā)空降碎屑以浮巖為主，浮巖碎屑的形態(tài)很不規(guī)則，從等軸狀到拉長狀，從次圓狀到棱角狀都有。隨著浮巖碎屑粒徑的減小，顆粒拉長現(xiàn)象明顯，拉長碎屑的比例從15.02%上升到47.5%。并且，浮巖碎屑的粒度越小，顆粒棱角狀越明顯，＞1mm的浮巖碎屑屬于次圓狀，0.5～1mm的碎屑屬于次棱角狀，0.25～0.5mm的碎屑屬于棱角狀。

(2)通過對(duì)樣品的形貌分析，得到其平均形態(tài)參數(shù)F=0.72。如果假設(shè)長白山天池火山千年大噴發(fā)浮巖碎屑為球形(F=1)時(shí)，計(jì)算的最終沉降速度是浮巖碎屑形態(tài)參數(shù)F=0.72時(shí)最終沉降速度的1.52倍左右，而當(dāng)假設(shè)顆粒為橢球(F=0.5)時(shí)，計(jì)算的最終沉降速度是浮巖碎屑形態(tài)參數(shù)F=0.72時(shí)最終沉降速度的0.89倍左右。

(3)風(fēng)速為10m/s時(shí)，不同的碎屑形態(tài)參數(shù)的沉積厚度等值線圖顯示:形態(tài)參數(shù)F=0.5時(shí)，擴(kuò)散的范圍最廣，F(xiàn)=0.72時(shí)次之，F(xiàn)=1時(shí)擴(kuò)散范圍最小，說明火山碎屑的擴(kuò)散范圍隨著形態(tài)參數(shù)F值的減小而擴(kuò)大。

以上的分析表明，碎屑的形貌特征對(duì)最終沉降速度和擴(kuò)散范圍都具有較大的影響。本次對(duì)＜0.25mm的碎屑未進(jìn)行形貌分析，而體視鏡下觀察它們的形態(tài)更為不規(guī)則。因此，如果按球形計(jì)算這些碎屑的最終沉降速度，產(chǎn)生的誤差必然會(huì)更大。顯然，為了更為準(zhǔn)確地模擬和預(yù)測火山噴發(fā)碎屑的擴(kuò)散范圍，需要對(duì)火山碎屑的形貌進(jìn)行更為詳盡的特征分析。此外，本次研究的主要目的是討論碎屑形貌對(duì)最終沉降速度和擴(kuò)散范圍的影響。在研究過程中我們發(fā)現(xiàn)，不同粒徑下碎屑顆粒的形貌參數(shù)有所不同，而以往在火山碎屑擴(kuò)散模型中浮巖碎屑形態(tài)參數(shù)值取的是一個(gè)定值，未詳細(xì)區(qū)分不同粒級(jí)下的形貌參數(shù)的變化，顯然也是會(huì)造成較大的誤差，更何況空降碎屑中還含有巖屑、晶屑、玻屑等其它成分，它們具有不同的形貌參數(shù)(Riley etal.，2003)。因此，在實(shí)際火山空降物災(zāi)害模擬時(shí)這些因素都應(yīng)該加以考慮。還應(yīng)指出的是，本研究中的最終沉降速度公式中只應(yīng)用了AR參數(shù)，相信今后隨著對(duì)顆粒最終沉降速度與其形貌參數(shù)之間關(guān)系的更深入了解，其它的參數(shù)(如CC、FF和RT等)也可以應(yīng)用到公式中。

樊祺誠，隋建立，孫謙，等.2005.天池火山千年大噴發(fā)的巖漿混合作用與噴發(fā)機(jī)制初步探討［J］.巖石學(xué)報(bào)，21(6):1703—1708.

FAN Qi-cheng，SUI Jian-li，SUN Qian，et al.2005.Preliminary research ofmagmamixing and explosivemechanism of themillennium eruption of Tianchi volcano［J］.Acta Petrologica Sinica，21(6):1703—1708(in Chinese).

劉若新，魏海泉，李繼泰，等.1998.長白山天池火山近代噴發(fā)［M］.北京:科學(xué)出版社.83—107.

LIU Ruo-xin，WEIHai-quan，LI Ji-tai，et al.1998.The Recent Eruption of Tianchi Volcano，Changbaishan Mountains［M］.Science Press，Beijing.83—107(in Chinese).

劉若新.2000.中國的活火山［M］.北京:地震出版社.17—31.

LIU Ruo-xin.2000.The Active Volcanoes in China［M］.Seismological Press，Beijing.17—31(in Chinese).

史蘭斌，陳孝德，楊清福，等.2005.長白山天池火山千年大噴發(fā)不同顏色浮巖的巖石化學(xué)特征［J］.地震地質(zhì)，27(1):73—82.

SHILan-bin，CHEN Xiao-de，YANG Qing-fu，et al.2005.Petrochemistry of pumices of various colors produced by the eruption of Changbaishan Tianchi volcano at 1000 years ago［J］.Seismology and Geology，27(1):73—82(in Chinese).

魏海泉，金伯祿，劉永順.2004.長白山天池火山地質(zhì)學(xué)研究的若干進(jìn)展與災(zāi)害分析［J］.巖石礦物學(xué)雜志，23(4):305—312.

WEIHai-quan，JIN Bo-lu，LIU Yong-shun.2004.Some advances in the study of volcanic geology and a hazards analysis of Tianchi volcano［J］.Acta Petrologica etMineralogica，23(4):305—312(in Chinese).

楊清福，史蘭斌，張羽，等.2007.長白山天池火山千年大噴發(fā)火山碎屑流堆積的粒度特征與地質(zhì)意義［J］.地震地質(zhì)，29(3):480—491.

YANG Qing-fu，SHI Lan-bin，ZHANG Yu，et al.2007.Grain-size characteristics of themillennium pyroclastic flow depoits of Tianchivolcano and their geological implications［J］.Seismology and Geology，29(3):480—491(in Chinese).

尹金輝，鄭勇剛，劉粵霞.2005.長白山天池火山炭化木的14C年代及其意義［J］.地震地質(zhì)，27(1):83—88.

YIN Jin-hui，ZHENG Yong-gang，LIU Yue-xia.2005.The radiocarbon age of carbinized wood in Tianchi volcano，Changbaishan Mountains and its implication［J］.Seismology and Geology，27(1):83—88(in Chinese).

于紅梅，許建東，趙誼.2007.長白山天池火山千年大噴發(fā)空降碎屑物的數(shù)值模擬［J］.地震地質(zhì)，29(3):522—534.

YU Hong-mei，XU Jian-dong，ZHAO Yi.2007.A numerical simulation of tephra transport and deposition formillennium eruption of Changbaishan Tianchi volcano［J］.Seismology and Geology，29(3):522—534(in Chinese).

趙誼.2003.火山危險(xiǎn)性和空降碎屑災(zāi)害的研究［D］:［學(xué)位論文］.安徽:中國科學(xué)技術(shù)大學(xué).

ZHAO Yi.2003.Study of volcanic hazard assessmentand the risk of airborne substances［D］.Dissertation for the Doctoral Degree.University of Science and Technology of China，Anhui(in Chinese).

趙誼，張程遠(yuǎn)，席道瑛.2002.火山碎屑空降沉積的二維數(shù)值模擬［J］.地震地質(zhì)，24(3):377—386.

ZHAO Yi，ZHANG Cheng-yuan，XIDao-ying.2002.2D-numerical simulation of dispersion of tephra fallout［J］.Seis-mology and Geology，24(3):377—386(in Chinese).

Bonadonna C，Ernst G G J，Sparks R S J.1998.Thickness variations and volume estimates of tephra fall deposits:The importance of particle Reynolds number［J］.JVolcanol Geotherm Res，81:173—187.

Bonadonna C，Connor CB，Houghton B F，etal.2005.Probabilisticmodeling of tephra dispersion:Hazard assessmentof amulti-phase rhyolitic eruption at Tarawera，New Zealand ［J］.JGeophys Res，110(B3):Art.No.B03203.doi:10.1029/2003JB002896.

Bonadonna C，Phillips JC.2003.Sedimentation from strong volcanic plumes［J］.JGeophys Res，108(B7):Art.No.2340.doi:10.1029/2002JB002034.

ColtelliM，Miraglia L，Scollo S.2008.Characterization of shape and terminal velocity of tephra particles erupted during the 2002 eruption of Etna volcano，Italy［J］.Bull Volcanol，70(9):1103—1112.doi:10.1007/s 00445-007-0192—8.

Dellino P，La Volpe L.1996.Image processing analysis in reconstructing fragmentation and transportationmechanisms of pyroclastic deposits:The case of Monte Pilato-Rocche Rosse eruptions，Lipari(Aeolian Islands，Italy)［J］.JVolcanol Geotherm Res，71:13—29.

Dellino P，Mele D，Bonaria R，et al.2005.The analysis of the influence of pumice shape on its terminal velocity ［J］.Geophys Res Lett32，Art.No.L21306.doi:10.1029/2005GL023954.

Folk R L.1974.Petrology of Sedimentary Rocks［M］.Hemphill，Austin，Texas.182.

Heffter JL，Stunder B JB.1993.Volcanic ash forecast transport and dispersion(Vaftad)model［J］.Weather Forecast，8:533—541.

Kunii D K，Levenspiel O.1969.Fluidization Engineering［M］.Wiley，New York.97.

Macedonio G，Costa A，Longo V.2005.A computermodel for volcanic ash fallout and assessment of subsequent hazard［J］.Comput Geosci，31:837—845.

Machida H，Arai F.1983.Extensive ash falls in and around the Sea of Japan from large Quaternary eruption ［J］.Jour Volcan Geother Res，18:151—164.

Machida H，Sugiyama S.1990.The recentmajor eruption of Changbai Volcano and its environmental effects［J］.Geographical Report of Tokyo Metropolotan University，25:1—20.

Muller G.1967.Methods in Sedimentary Petrology［M］.Hafner，New York.

Pfeiffer T，Costa A，Macedonio G.2005.A model for the numerical simulation of tephra fall deposits［J］.JVolcanol Geotherm Res，140:273—294.

Riley CM，Rose W Ⅰ，Bluth G JS.2003.Quantitative shapemeasurements of distal volcanic ash［J］.JGeophys Res，108(B10)，Art.No.2504.doi:10.1029/2001JB000818.

Russ J.1995.The Image Processing Handbook ［M］.CRC，Boca Raton，F(xiàn)L.674.

Sparks R S J，Bursik M Ⅰ，Carey SN，et al.1997.Volcanic Plumes［M］.Wiley，Chichester.

Searcy C，Dean K，Stringer W.1998.PUFF:A high resolution volcanic ash trackingmodel［J］.JVolcanol Geotherm Res，80:1—16.

Suzuki T.1983.A theoretical model for dispersion of tephra ［A］.In:Shimozuru D，Yokoyama I eds.，Volcanism:Physics and Tectonics.Arc，Tokyo:95—113.

Wilson L，Huang TC.1979.The influence of shape on the atmospheric settling velocity of volcanic ash particles［J］.Earth Planet Sci Lett，44:311—324.

MORPHOLOGICAL CHARACTERIZATION AND TERM INAL VELOCITY OF PUM ICE PARTICLES ERUPTED DURING THE M ILLENNIUM ERUPTION OF CHANGBAISHAN TIANCHIVOLCANO，CHINA

YU Hong-mei XU Jian-dong LIN Chuan-yong
(Institute of Geology，China Earthquake Administration，Beijing 100029，China)

The shape parameter of particles is an important parameter in the formula of the terminal settling velocity，but itwas usually assumed to be spherical or ellipsoidal without detailed study.According to themorphological characterization of the tephra particles erupted during the Millennium eruption of Changbaishan Tianchi volcano，we got some new recognition:The tephra particles consistmainly of pumice，which have irregular shapes，varying from very equant to elongate and from subrounded to angular.With decreasing of the particle size，the elongation of particles becomes significant，and the proportion of elongated particles increases from 15.02%to 47.5%.Moreover，the smaller the size of particles is，themore angular the shape of particles.An average shape parameter F=0.72 has been obtained through themorphological characterization of the pumice particles.The terminal velocity of the pumice particles，which were assumed to be spherical(F=1)，erupted during the Millennium eruption of the Changbaishan Tianchi volcano is about1.52 times the terminal velocity of the pumice particles of F=0.72，and the terminal velocity of the pumice particles，which were assumed to be ellipsoidal(F=0.5)，is about0.89 times the terminal velocity of the pumice particles of F=0.72.When the wind speed is taken to be 10m/s，the isopach map of tephra particleswith different shape parameters shows that the dispersal range expands with the decrease of the shape parameter F.The aforementioned results show clearly that the morphology of particles influences strongly both the terminal settling velocity and dispersal range，and hence should be fully considered when simulating tephra dispersion and forecasting disasters in the future.

pumice particles，shape parameters，terminal settling velocity，tephra dispersal model，Changbaishan Tianchi volcano

P317.3

A

0253-4967(2011)02-0440-12

10.3969/j.issn.0253-4967.2011.02.017

2010-09-12收稿，2010-11-22改回。

國家自然科學(xué)基金(40972209)和地震行業(yè)科研專項(xiàng)(20070827)共同資助。

于紅梅，女，1981年生，中國地震局地質(zhì)研究所博士研究生，主要從事火山學(xué)研究，電話:010-62009135，E-mail:yuhongmei188@163.com。