基于空間插值方法的2021年云南漾濞MS6.4地震儀器地震烈度影響場(chǎng)分析

于江 和仕芳 盧永坤 龐衛(wèi)東 羅偉東 鐘玉盛

摘要：基于2021年云南漾濞MS6.4主震預(yù)警臺(tái)網(wǎng)地震動(dòng)數(shù)據(jù)，使用多種空間插值方法探討獲得最優(yōu)儀器地震烈度影響場(chǎng)，并與宏觀地震烈度調(diào)查結(jié)果進(jìn)行對(duì)比分析。結(jié)果表明：基于指數(shù)模型變異函數(shù)的普通克里金插值法獲得的儀器地震烈度影響場(chǎng)與宏觀地震烈度具有較好的一致性，能夠反映實(shí)際地震烈度空間總體影響場(chǎng)范圍及特征;除了兩者因概念和屬性不同而造成的影響場(chǎng)差異之外，臺(tái)站分布、場(chǎng)地條件、空間插值方法、前震和余震震害疊加等因素也是導(dǎo)致本次地震影響場(chǎng)差異的主要原因。

關(guān)鍵詞：儀器烈度;地震烈度;空間插值;預(yù)警臺(tái)網(wǎng);漾濞MS6.4地震

中圖分類號(hào)：P315.9?? 文獻(xiàn)標(biāo)識(shí)碼：A?? 文章編號(hào)：1000-0666（2021）03-0414-08

0 引言

地震烈度是衡量地震對(duì)地表和工程建筑結(jié)構(gòu)影響和破壞程度的重要參數(shù)。在破壞性地震發(fā)生后快速評(píng)定地震烈度影響場(chǎng)可以為開展地震應(yīng)急救援、力量部署、災(zāi)害損失評(píng)估等工作提供重要的科學(xué)依據(jù)（周光全等，2006;盧永坤等，2019;張建國(guó)等，2012）。傳統(tǒng)的宏觀地震烈度評(píng)定是建立在人工調(diào)查的基礎(chǔ)上，根據(jù)多種地震宏觀影響后果來(lái)評(píng)定地震烈度，往往花費(fèi)較大的人力物力、耗時(shí)較長(zhǎng)（田秀豐等，2020）。為了滿足震后應(yīng)急快速?zèng)Q策的需要，利用強(qiáng)震動(dòng)臺(tái)網(wǎng)觀測(cè)記錄計(jì)算儀器地震烈度獲取地震烈度影響場(chǎng)的方法逐步發(fā)展起來(lái)。

儀器地震烈度是指根據(jù)儀器觀測(cè)記錄計(jì)算得到的地震動(dòng)強(qiáng)度，并能直接反映地震的影響程度，具有直觀、簡(jiǎn)便、快速等特點(diǎn)（金星等，2013;劉如山等，2021）。近年來(lái)，隨著國(guó)家地震烈度速報(bào)和預(yù)警項(xiàng)目的建設(shè)，我國(guó)大陸地區(qū)正在布設(shè)數(shù)千個(gè)速度計(jì)臺(tái)站和上萬(wàn)個(gè)加速度計(jì)臺(tái)站，其中部分已投入運(yùn)行，極大地提升了目前地震動(dòng)觀測(cè)臺(tái)站的數(shù)量和密度，使得震后利用高密度預(yù)警臺(tái)網(wǎng)儀器烈度進(jìn)行空間插值快速獲取有效的儀器地震烈度影響場(chǎng)成為了可能。2021年5月21日21時(shí)48分34秒，云南省漾濞縣（25.67°N，99.87°E）發(fā)生MS6.4地震。本次地震是云南地震烈度速報(bào)和預(yù)警項(xiàng)目建成后發(fā)生的首次強(qiáng)震，高密度的預(yù)警臺(tái)網(wǎng)獲得了豐富的地震動(dòng)數(shù)據(jù)。本文以漾濞MS6.4地震為例，使用預(yù)警臺(tái)網(wǎng)地震動(dòng)數(shù)據(jù)，嘗試運(yùn)用空間插值方法獲得與宏觀地震烈度有較好吻合度的儀器地震烈度影響場(chǎng)，然后將其插值結(jié)果與宏觀地震烈度調(diào)查結(jié)果進(jìn)行對(duì)比分析，并進(jìn)一步探討儀器地震烈度影響場(chǎng)空間插值修正方法。

1 地震動(dòng)觀測(cè)數(shù)據(jù)

云南省地震局作為國(guó)家地震烈度速報(bào)和預(yù)警項(xiàng)目的先行先試單位，于2020年底完成了全省預(yù)警臺(tái)網(wǎng)的建設(shè)工作，省內(nèi)共建設(shè)一千余個(gè)臺(tái)站，其中烈度計(jì)站占比80%，重點(diǎn)預(yù)警區(qū)內(nèi)臺(tái)均間距約為10 km，可以達(dá)到秒級(jí)的地震預(yù)警，目前正處于試運(yùn)行狀態(tài)。據(jù)云南地震臺(tái)網(wǎng)測(cè)定，截至2021年5月29日12時(shí)00分漾濞地震序列共發(fā)生M≥3.0地震46次，其中6.0～6.9級(jí)地震1次，5.0～5.9級(jí)地震3次，4.0～4.9級(jí)地震13次，3.0～3.9級(jí)地震29次。漾濞MS6.4主震350 km范圍內(nèi)共有241個(gè)預(yù)警臺(tái)站記錄到數(shù)據(jù)，其中一般站219個(gè)、基本站19個(gè)、基準(zhǔn)站3個(gè)。位于漾濞縣城的強(qiáng)震動(dòng)臺(tái)獲得主震記錄，最大加速度為719.6 gal。對(duì)數(shù)據(jù)進(jìn)行初步篩查，剔除了與震中距離較遠(yuǎn)且相對(duì)孤立的臺(tái)站13個(gè)，對(duì)兩個(gè)空間點(diǎn)位重復(fù)的臺(tái)站僅保留了一個(gè)與周邊站點(diǎn)數(shù)據(jù)變化范圍相近的樣本點(diǎn)，最終保留227個(gè)臺(tái)站的地震動(dòng)數(shù)據(jù)，最小臺(tái)間距為2 km，平均臺(tái)間距約為13 km。圖1為漾濞MS6.4地震宏觀地震烈度與預(yù)警臺(tái)站分布圖。

2 儀器地震烈度影響場(chǎng)空間插值方法

2.1 儀器地震烈度計(jì)算

目前，國(guó)內(nèi)外眾多學(xué)者根據(jù)不同地震動(dòng)特性的地震動(dòng)參數(shù)與地震烈度的關(guān)系確定了相應(yīng)的儀器地震烈度算法（關(guān)曙淵等，2020;金星等，2013;Wu et al，2003;崔建文等，2006）?！吨袊?guó)地震烈度表》（GB/T 17442—2020）給出了我國(guó)現(xiàn)行儀器地震烈度的計(jì)算流程及其定量化測(cè)定方法，在對(duì)數(shù)據(jù)進(jìn)行基線校正、記錄轉(zhuǎn)換、數(shù)字濾波后，三分向合成地震動(dòng)參數(shù)計(jì)算公式如下：

PGA=maxa2（ti）E-W+a2（ti）N-S+a2（ti）U-D（1）

PGV=maxv2（ti）E-W+v2（ti）N-S+v2（ti）U-D（2）

根據(jù)式（1）（2）分別計(jì)算IA、IV值：

IA=3.17lg（PGA）+6.59（3）

IV=3.00lg（PGV）+9.77（4）

最后根據(jù)IA和IV值計(jì)算儀器測(cè)定的地震烈度II為：

II=IV IA≥6.0且IV≥6.0IA+IV2，? （IA<6.0且IV<6.0）（5）

式中：PGA為合成地震動(dòng)加速度記錄的最大值，單位為m/s2; PGV為合成地震動(dòng)加速度記錄的最大值，單位為m/s; a（ti）為ti時(shí)刻點(diǎn)合成的三分向加速度值; v（ti）為ti時(shí)刻點(diǎn)合成的三分向速度值;下標(biāo)E-W、N-S、U-D分別表示東西方向、南北方向、垂直方向。

2.2 空間插值方法選取

空間插值是一種通過離散空間數(shù)據(jù)計(jì)算未知空間數(shù)據(jù)的方法，一般分為確定性方法和地質(zhì)統(tǒng)計(jì)學(xué)方法兩類（李海濤，邵澤東，2019）。

確定性插值方法是基于已知樣本點(diǎn)之間的相似程度或者整個(gè)曲面的光滑性來(lái)創(chuàng)建一個(gè)擬合曲面。主要有反距離權(quán)重法趨勢(shì)面法、樣條函數(shù)法等;其中反距離權(quán)重法作為一種精確插值方法，算法簡(jiǎn)便易行，當(dāng)樣本數(shù)據(jù)分布均勻、密度較高足以反映局部差異時(shí)具有較好的插值精度（朱求安等，2004;張海平等，2017）。

地質(zhì)統(tǒng)計(jì)學(xué)插值方法是利用樣本點(diǎn)的統(tǒng)計(jì)規(guī)律，使樣本點(diǎn)之間的空間自相關(guān)性定量化，從而在待預(yù)測(cè)的點(diǎn)周圍構(gòu)建樣本點(diǎn)的空間結(jié)構(gòu)模型，比如克立金插值法?？死锝鸩逯捣ㄗ鳛榈刭|(zhì)統(tǒng)計(jì)學(xué)的主要內(nèi)容是一種以變異函數(shù)理論和結(jié)構(gòu)分析為基礎(chǔ)，在有限區(qū)域內(nèi)對(duì)區(qū)域化變量進(jìn)行無(wú)偏最優(yōu)估計(jì)的一種方法，一般分為普通克里金、簡(jiǎn)單克里金、協(xié)克里金、泛克里金、指示克里金等（牟乃夏等，2012;鄭向向，帥向華，2013）。

克里金插值法與反距離權(quán)重法類似，兩者均通過已知樣本點(diǎn)權(quán)重來(lái)求得未知樣本點(diǎn)的值，可統(tǒng)一表示為：

Z（x0）=∑ni=1λiZ（xi）（6）

式中：Z（x0）為未知樣本點(diǎn)的值;Z（xi）為未知樣本點(diǎn)周圍的已知樣本點(diǎn)的值;λi為第i個(gè)已知樣本點(diǎn)對(duì)未知樣本點(diǎn)的權(quán)重;n為已知樣本點(diǎn)的個(gè)數(shù)。

相對(duì)于反距離權(quán)重法只考慮已知樣本點(diǎn)與未知樣本點(diǎn)的距離遠(yuǎn)近，克里金插值法還通過變異函數(shù)和結(jié)構(gòu)分析，考慮了已知樣本點(diǎn)的空間分布及與未知樣本點(diǎn)的空間方位關(guān)系，能較好反映客觀地質(zhì)規(guī)律。其次，地質(zhì)統(tǒng)計(jì)學(xué)方法的插值精度與選擇的變異函數(shù)模型密切相關(guān)，常用的變異函數(shù)理論模型有圓形模型、球形模型、指數(shù)模型、高斯模型等（湯國(guó)安，楊昕，2012）。

本文基于以上兩類插值方法基本原理，考慮到此次漾濞地震具有較密集的預(yù)警臺(tái)網(wǎng)數(shù)據(jù)，為尋找最優(yōu)空間插值方法，分別選擇反距離權(quán)重法和普通克里金法對(duì)離散儀器烈度值進(jìn)行空間插值。

通常在插值前需要根據(jù)所選的空間插值方法，對(duì)樣本數(shù)據(jù)分布特征、內(nèi)在規(guī)律性以及趨勢(shì)效應(yīng)等進(jìn)行分析，驗(yàn)證確定權(quán)重系數(shù)及相關(guān)參數(shù)或選擇合適的變異函數(shù)模型，從而提高空間插值的準(zhǔn)確率。在使用地質(zhì)統(tǒng)計(jì)學(xué)方法分析時(shí)，首先需要對(duì)數(shù)據(jù)是否符合正態(tài)分布進(jìn)行驗(yàn)證，若不符合正態(tài)分布，則需對(duì)數(shù)據(jù)進(jìn)行變換。如圖2a為本次地震儀器烈度數(shù)據(jù)直方圖，柱狀圖整體呈左偏的鐘形曲線，顯示數(shù)據(jù)接近于正態(tài)分布。正態(tài)QQ分布圖通常用來(lái)評(píng)估兩個(gè)數(shù)據(jù)集分布的相似程度，圖2b顯示了該組數(shù)據(jù)整體延45°方向近似呈直線展布，同樣顯示數(shù)據(jù)具有正態(tài)分布特征。因而該組數(shù)據(jù)滿足地質(zhì)統(tǒng)計(jì)學(xué)要求，可以使用克里金插值方法。

在最優(yōu)權(quán)重或變異函數(shù)理論模型選擇方面，使用反距離權(quán)重法插值時(shí)，通過預(yù)覽輸出和對(duì)比驗(yàn)證確定插值最優(yōu)權(quán)重為3。使用普通克里金插值方法時(shí)，根據(jù)步長(zhǎng)及步長(zhǎng)數(shù)選擇標(biāo)準(zhǔn)，一般最大步長(zhǎng)取值應(yīng)該小于研究區(qū)域采樣點(diǎn)之間最大尺度的一半，并通過確定點(diǎn)與最近的相鄰要素之間的平均距離來(lái)確定步長(zhǎng)大小，取步長(zhǎng)值為0.169，同時(shí)為了滿足步長(zhǎng)數(shù)乘以步長(zhǎng)大小的值是所有點(diǎn)之間最大距離的一半左右的經(jīng)驗(yàn)關(guān)系，確定步長(zhǎng)組數(shù)為12;最后，根據(jù)數(shù)據(jù)空間變量的相關(guān)性，普通克里金變異函數(shù)理論模型分別選擇球面模型、指數(shù)模型、高斯模型進(jìn)行空間插值。

2.3 空間插值結(jié)果分析

為了評(píng)價(jià)不同方法獲得的空間插值結(jié)果，本文將空間插值結(jié)果按照地震烈度值進(jìn)行分類顯示，以空間插值結(jié)果與宏觀地震烈度的吻合程度作為評(píng)價(jià)依據(jù)，從而判斷插值結(jié)果的優(yōu)劣。本次漾濞地震現(xiàn)場(chǎng)共調(diào)查了218個(gè)居民點(diǎn)的震害情況，確定了本次地震最高烈度為Ⅷ度，Ⅵ度區(qū)及以上面積約為6 600 km2，共涉及大理州6個(gè)縣（市）。圖3為反距離權(quán)重法和普通克里金法對(duì)應(yīng)的4種儀器地震烈度空間插值結(jié)果與宏觀地震烈度的對(duì)比圖。

如圖3a所示，使用反距離權(quán)重法獲得的儀器烈度影響場(chǎng)與實(shí)際宏觀烈度邊界在Ⅷ、Ⅶ度區(qū)的東邊界擬合度較好，Ⅵ度區(qū)擬合度較差。擬合得到的儀器烈度影響場(chǎng)長(zhǎng)軸走向不明顯，尤其是Ⅶ度區(qū)的長(zhǎng)軸走向整體呈北東向，與宏觀地震烈度差別較大。在巍山盆地東緣存在一個(gè)孤立Ⅶ度異常區(qū)，該區(qū)域儀器地震烈度值明顯高于周圍儀器烈度值，與宏觀地震烈度有所差別。

對(duì)于地質(zhì)統(tǒng)計(jì)學(xué)方法，不同變異函數(shù)理論模型所得到的插值結(jié)果差距較大。如球面模型與高斯模型（圖3b、c）對(duì)高烈度區(qū)的插值效果較差，基本沒有擬合出Ⅷ度區(qū)的范圍，儀器烈度影響場(chǎng)邊界與實(shí)際宏觀烈度邊界擬合吻合度較差。但其影響場(chǎng)范圍的整體形態(tài)與實(shí)際宏觀調(diào)查所獲得的長(zhǎng)軸走向擬合較好，可以明顯看出儀器烈度影響場(chǎng)長(zhǎng)軸走向?yàn)楸蔽飨颉?/p>

從指數(shù)模型的普通克里金插值結(jié)果可以看出，在高烈度區(qū)該方法依舊可以根據(jù)僅有的一個(gè)高烈度值擬合出較小的Ⅷ度區(qū)，在Ⅶ、Ⅵ度區(qū)的東南邊界插值結(jié)果和實(shí)際宏觀烈度也比較吻合，在Ⅵ度區(qū)西南邊界可以看出受預(yù)警臺(tái)站間距的影響，僅在有樣本數(shù)據(jù)的位置上與宏觀Ⅵ度區(qū)邊界擬合度較好，其邊界呈花瓣?duì)?。其影響?chǎng)Ⅵ度區(qū)與Ⅶ度區(qū)長(zhǎng)軸方向也較明顯，與宏觀地震烈度相符（圖3d）。

結(jié)合上述討論的4種插值結(jié)果可知，在極震區(qū)內(nèi)臺(tái)站數(shù)量較少的情況下，產(chǎn)出的儀器地震烈度空間分布圖更加依賴于插值算法的適用性，不同的插值方法、參數(shù)、變異函數(shù)理論模型均會(huì)對(duì)插值結(jié)果造成較大影響。本次漾濞地震儀器地震烈度空間插值方法中基于指數(shù)模型變異函數(shù)的普通克里金方法的插值結(jié)果與宏觀地震烈度整體范圍及特征吻合度較好，反距離權(quán)重法插值結(jié)果次之。

3 與宏觀地震烈度對(duì)比分析

根據(jù)基于指數(shù)模型變異函數(shù)的普通克里金插值方法獲得的最優(yōu)儀器烈度影響場(chǎng)與實(shí)際宏觀地震烈度進(jìn)行對(duì)比，分析導(dǎo)致其影響場(chǎng)差異的主要因素。

3.1 儀器烈度分布特征

對(duì)儀器烈度空間插值結(jié)果與宏觀地震烈度空間范圍進(jìn)行疊加分析，如圖4所示，在宏觀地震烈度Ⅵ度區(qū)范圍內(nèi)，存在部分儀器烈度值小于宏觀地震烈度的情況。Ⅵ度區(qū)內(nèi)共有預(yù)警臺(tái)14個(gè)，其中8個(gè)預(yù)警臺(tái)儀器烈度值小于Ⅵ度，共有7個(gè)預(yù)警臺(tái)位于宏觀地震烈度Ⅵ度區(qū)的西北部，僅有1個(gè)預(yù)警臺(tái)位于東南部。儀器烈度影響場(chǎng)存在沿長(zhǎng)軸方向在Ⅷ度、Ⅶ度、Ⅵ度區(qū)東南方向面積普遍較大的情況，其原因可能是在震中東南方向地震動(dòng)衰

減較慢，導(dǎo)致在主震東南側(cè)記錄到的儀器地震烈度普遍比西北側(cè)高，且距震中更遠(yuǎn)。而實(shí)際調(diào)查的宏觀地震烈度則沿長(zhǎng)軸方向在微觀震中兩側(cè)面積相差不大，其原因可能是儀器烈度數(shù)據(jù)僅為漾濞MS6.4主震產(chǎn)生的地震動(dòng)強(qiáng)度所計(jì)算的烈度值，而漾濞地震除主震以外，前震和余震頻發(fā)，其中5級(jí)以上地震就有3次，宏觀地震烈度為主震及多次前震、余震疊加所造成的震害，因而在漾濞地震中同一地區(qū)宏觀烈度值普遍高于儀器烈度值。

在漾濞地震震中附近僅有漾濞臺(tái)儀器地震烈度大于Ⅷ度，所以導(dǎo)致僅通過儀器烈度插值獲得的Ⅷ度區(qū)范圍較小，與實(shí)際宏觀烈度Ⅷ度區(qū)范圍相差較大。其次，在宏觀地震烈度Ⅵ度區(qū)南部，受預(yù)警臺(tái)站分布不均勻且密度較低的影響，導(dǎo)致在該處空間插值形成的儀器地震烈度影響場(chǎng)Ⅵ度邊界呈花瓣?duì)?，與實(shí)際調(diào)查烈度存在一定差異?？梢娕_(tái)站密度和分布均勻程度是造成儀器地震烈度與宏觀地震烈度差異的另一個(gè)影響因素。

由于宏觀地震烈度反映的是震害水平的綜合指標(biāo)，影響因素較多，且在綜合制圖中有平均化的趨勢(shì)，因而其影響場(chǎng)邊界更為均勻光滑，而儀器烈度值則直接由觀測(cè)的地震動(dòng)記錄計(jì)算得到，兩者在概念和屬性上存在一定差異，因而兩者的影響場(chǎng)分布存在差異是不可避免的。

3.2 觀測(cè)臺(tái)站場(chǎng)地局部特征

震區(qū)屬橫斷山滇西高山峽谷區(qū)，地形起伏較大，現(xiàn)場(chǎng)調(diào)查發(fā)現(xiàn)，震害情況受地形、場(chǎng)地影響較大。通?，F(xiàn)場(chǎng)調(diào)查是以居民點(diǎn)為單位對(duì)該區(qū)域的房屋震害整體情況進(jìn)行評(píng)估，而預(yù)警臺(tái)記錄的儀器烈度值則反映的是臺(tái)站所處區(qū)域的場(chǎng)地地震動(dòng)特征。因目前大部分烈度計(jì)與通信基站合作架設(shè)，其周邊地形及場(chǎng)地條件可能往往與居民點(diǎn)所在的環(huán)境有所不同，從而造成儀器烈度值與實(shí)際調(diào)查結(jié)果有一定的差異。

在地形起伏緩和、場(chǎng)地覆蓋層較厚的地區(qū)，場(chǎng)地的放大作用對(duì)儀器烈度具有較明顯的影響。圖4中在Ⅵ度區(qū)東南部，巍山盆地東緣大倉(cāng)臺(tái)（震中距47 km）儀器烈度值為7.1，而位于其西南側(cè)山區(qū)馬鞍山臺(tái)（震中距43 km）儀器烈度值僅為5.9。在地形起伏較大的山區(qū)，臺(tái)站周邊場(chǎng)地通常覆蓋層較薄，而居民點(diǎn)通常位于平緩的坡地、谷地等覆蓋層相對(duì)較厚的地方，從而導(dǎo)致儀器烈度值往往低于實(shí)際調(diào)查烈度值。關(guān)坪臺(tái)位于Ⅵ度區(qū)西北部，臺(tái)站架設(shè)在覆蓋層較薄的坡地上（圖5a），本次記錄的儀器烈度值僅為4.1，而在其北側(cè)、東側(cè)約2 km的關(guān)坪村舊街小組與自新村實(shí)際宏觀烈度值均為Ⅵ度（圖5b），兩者烈度值相差較大。另外，臺(tái)站與居民點(diǎn)實(shí)際所處的地形地貌對(duì)震害情況也有著明顯的控制作用，往往造成兩者烈度值的差異。

4 儀器烈度影響場(chǎng)修正

本次漾濞地震在僅使用預(yù)警臺(tái)網(wǎng)儀器地震烈度進(jìn)行空間插值所獲得的Ⅵ度至Ⅷ度影響場(chǎng)與宏觀地震烈度仍有一定出入，尤其是在Ⅷ度高烈度區(qū)臺(tái)站數(shù)量很少的情況下，兩者差異更為明顯。為了獲得更加符合實(shí)際災(zāi)情的地震烈度影響場(chǎng)，可動(dòng)態(tài)利用有效數(shù)據(jù)對(duì)插值參數(shù)或樣本數(shù)據(jù)進(jìn)行修正。

儀器烈度影響場(chǎng)快速評(píng)估通常在震后5 min以內(nèi)進(jìn)行，此時(shí)可以根據(jù)主震震級(jí)及其烈度衰減關(guān)系，考慮斷層走向，初步確定影響場(chǎng)長(zhǎng)軸方向，調(diào)整空間插值參數(shù)，針對(duì)極震區(qū)觀測(cè)臺(tái)站較少的情況，可以增加微觀震中點(diǎn)作為高烈度區(qū)樣本數(shù)據(jù)，修正儀器烈度影響場(chǎng)。震后24 h，可以根據(jù)現(xiàn)場(chǎng)實(shí)際調(diào)查結(jié)果，補(bǔ)充震中高烈度區(qū)樣本點(diǎn)、臺(tái)站分布空區(qū)與邊界控制樣本點(diǎn)，結(jié)合余震分布方向或余震精定位數(shù)據(jù)，優(yōu)化儀器烈度影響場(chǎng)長(zhǎng)軸方向，修正儀器烈度影響場(chǎng)。根據(jù)漾濞現(xiàn)場(chǎng)第一天實(shí)際調(diào)查結(jié)果，在儀器烈度樣本數(shù)據(jù)基礎(chǔ)上，補(bǔ)充Ⅷ度區(qū)樣本點(diǎn)，以及在Ⅵ度區(qū)、Ⅶ度區(qū)臺(tái)站空區(qū)及控制邊界補(bǔ)充樣本點(diǎn)，結(jié)合余震分布優(yōu)勢(shì)方向，使用基于指數(shù)變異函數(shù)的克里金插值方法獲得修正后的地震烈度影響場(chǎng)。從圖6可見，經(jīng)過修正優(yōu)化后的儀器地震烈度影響場(chǎng)與宏觀地震烈度更為接近，顯示出較好的修正效果，可以為現(xiàn)場(chǎng)宏觀地震烈度調(diào)查工作部署提供依據(jù)，如第二天可以在Ⅵ度區(qū)西北和東南空區(qū)方向部署力量，從而促進(jìn)現(xiàn)場(chǎng)工作開展。

5 結(jié)論

本文利用多種空間插值方法探討獲得漾濞地震最優(yōu)儀器地震烈度影響場(chǎng)，并與宏觀地震烈度調(diào)查結(jié)果進(jìn)行對(duì)比分析，得到以下結(jié)論：

（1）在本文討論的4種方法中，最優(yōu)空間插值法是基于指數(shù)模型變異函數(shù)的普通克里金插值法，其儀器地震烈度插值結(jié)果與宏觀地震烈度整體吻合度較高。

（2）儀器地震烈度與宏觀地震烈度影響場(chǎng)的差異主要受臺(tái)站分布、場(chǎng)地條件、空間插值方法、前震和余震震害疊加等因素的影響。同時(shí)，由于儀器地震烈度與宏觀地震烈度的概念、屬性和物理意義不同，因而兩者影響場(chǎng)的分布存在一定差異也是不可避免的。

（3）基于預(yù)警臺(tái)網(wǎng)的儀器地震烈度影響場(chǎng)較好地反映了本次漾濞地震宏觀地震烈度長(zhǎng)軸方向與地震動(dòng)衰減方向，真實(shí)地反映了地震動(dòng)強(qiáng)弱的空間范圍，與傳統(tǒng)地震烈度調(diào)查相比，具有速度快，范圍廣等優(yōu)點(diǎn)。

（4）針對(duì)本次地震儀器烈度與宏觀烈度分布差異，可以通過動(dòng)態(tài)利用有效數(shù)據(jù)對(duì)插值參數(shù)或樣本數(shù)據(jù)進(jìn)行修正的方法來(lái)進(jìn)一步改進(jìn)儀器地震烈度影響場(chǎng)分布圖的繪制。

本文在撰寫過程中得到了中國(guó)地震局昆明地震預(yù)報(bào)研究所崔建文研究員的指導(dǎo)和幫助，使用了云南地震臺(tái)提供的漾濞地震預(yù)警臺(tái)網(wǎng)數(shù)據(jù)，在此表示衷心感謝。

參考文獻(xiàn)：

崔建文，李世成，高東，等.2006.云南分區(qū)地震動(dòng)衰減關(guān)系[J].地震研究，29（4）：386-391.

關(guān)曙淵，單振東，馬榮.2020.日本儀器烈度計(jì)算3種方法對(duì)比[J].地震科學(xué)進(jìn)展，50（8）：20-24.

金星，張紅才，李軍，等.2013.地震儀器烈度標(biāo)準(zhǔn)初步研究[J].地球物理學(xué)進(jìn)展，28（5）：2336-2351.

李海濤，邵澤東.2019.空間插值分析算法綜述[J].計(jì)算機(jī)系統(tǒng)應(yīng)用，28（7）：1-8.

劉如山，熊明攀，馬強(qiáng)，等.2021.基于儀器地震烈度的變電站高壓電氣設(shè)備易損性研究[J].自然災(zāi)害學(xué)報(bào)，30（2）：14-23.

盧永坤，周洋，代博洋，等.2019.2018年云南墨江5.9級(jí)地震房屋震害特征與烈度評(píng)定[J].地震研究，42（2）：172-178.

牟乃夏，劉文寶，王海銀，等.2012.ArcGIS 10地理信息系統(tǒng)教程[M].北京：測(cè)繪出版社.

湯國(guó)安，楊昕.2012.ArcGIS地理信息系統(tǒng)空間分析實(shí)驗(yàn)教程[M].北京：科學(xué)出版社.

田秀豐，張衛(wèi)東，袁潔，等.2020.汶川8.0級(jí)地震儀器烈度與宏觀烈度對(duì)比分析[J].地震工程學(xué)報(bào)，42（5）：1225-1231.

張海平，周星星，代文.2017.空間插值方法的適用性分析初探[J].地理與地理信息科學(xué)，33（6）：14-18.

張建國(guó)，周光全，崔建文，等.2012.寧洱MS6.4地震考察研究綜述[J].地震研究，35（1）：151-155.

鄭向向，帥向華.2013.基于地質(zhì)統(tǒng)計(jì)方法與DEM的地震災(zāi)情空間插值研究[J].地震學(xué)報(bào)，35（4）：573-583.

周光全，非明倫，施偉華.2006.1992—2005年云南地震災(zāi)害損失與主要經(jīng)濟(jì)指標(biāo)研究[J].地震研究，29（2）：198-202.

朱求安，張萬(wàn)昌，余鈞輝.2004.基于GIS的空間插值方法研究[J].江西師范大學(xué)學(xué)報(bào)（自然科學(xué)版），25（2）：183-188.

Wu Y M，Teng T L，Shin T C，et al.2003.Relationship between peak ground acceleration，peak ground velocity，and intensity in Taiwan[J].Bulletin of the Seismological Society of America，93（1）：386-396.

GB/T 17442—2020，中國(guó)地震烈度表[S].

Analysis of the Influence Field of the Instrumental Seismic Intensity of the 2021 Yangbi，Yunnan MS6.4 Earthquake Basedon Spatial Interpolation Method

YU Jiang，HE Shifang，LU Yongkun，PANG Weidong，LUO Weidong，ZHONG Yusheng

（Yunnan Earthquake Agency，Kunming 650224，Yunnan，China）

Abstract

Based on the ground motion data of the Yangbi MS6.4 earthquake recorded by the Yunnan Early Warning Network，a variety of spatial interpolation methods are used to obtain the optimal instrumental intensity influence field，and compared with the results from the macroscopic intensity obtained from field survey.The results show that the influence field of instrumental intensity obtained by the ordinary Kriging interpolation method based on the variogram of the exponential model has a good agreement with the macroscopic seismic intensity，which can reflect the overall range of the influence field and characteristics of the actual seismic intensity.The main reasons for the differences between the instrumental intensity influence field and the macroscopic intensity field are the distribution of stations，site conditions，spatial interpolation methods，and the superposition of foreshocks damage and aftershocks damage.

Keywords：instrumental intensity;earthquake intensity;spatial interpolation;the Yunnan Early Warning Network;the Yangbi MS6.4 earthquake