基于遙感初判的建筑物震害Fisher判別法研究

摘要： 破壞性地震發(fā)生后，迅速準確地預(yù)測建筑物破壞程度，對快速科學(xué)地開展地震應(yīng)急指揮、救援力量部署等工作具有重大意義。針對現(xiàn)有建筑震害預(yù)測模型存在的評估結(jié)果粗糙、數(shù)據(jù)獲取難度較高、計算工作量大、模型構(gòu)建難度大、普適性不強等問題，提出一種基于遙感初判的建筑物震害Fisher判別法。首先，從地震強度和建筑物抗震能力兩方面選取震級、震中距、場地條件及建筑物抗震能力4種震害作為判別因子；然后，基于判別分析理論，構(gòu)建建筑物震害Fisher判別模型；最后，以四川瀘定縣6.0級地震為例，對文章提出方法進行驗證。實驗結(jié)果表明：該方法的預(yù)測結(jié)果與實際震害基本一致，準確率高達80%以上，證明該方法具有較高的可靠性，能較準確地對建筑物震害程度進行科學(xué)預(yù)測。

關(guān)鍵詞： 遙感初判; 建筑物; 震害預(yù)測; Fisher判別法; 瀘定縣地震

中圖分類號： TU315.9 文獻標志碼：A 文章編號： 1000-0844（2024）06-1405-10

DOI：10.20000/j.1000-0844.20221120001

A Fisher discriminant method for seismic damage of buildings

based on preliminary judgment from remote sensing

ZHAO Zhen¹， GUO Hongmei¹， ZHANG Ying¹， YIN Wen'gang²，

LU Changjiang¹， FAN Kaihong¹， ZHANG Yi¹

（1.Sichuan Earthquake Agency， Chengdu 610041， Sichuan， China;

2. Officers College of PAP， Chengdu 610041， Sichuan， China）

Abstract： After a destructive earthquake， quickly and accurately predicting the damage degree of buildings is of great importance to quickly and scientifically carry out earthquake emergency command and rescue force deployment. Existing building seismic damage prediction models yield imprecise evaluation results， struggle with data acquisition， and require a large amount of computing power; therefore， they are difficult to construct and lack universality. To address these issues， this paper proposes a Fisher discriminant method for seismic damage of buildings based on preliminary judgment from remote sensing. First， four seismic damage factors， including magnitude， epicentral distance， site condition， and seismic resistance of the building， were selected as the discriminant factors. Then， based on the discriminant analysis theory， the Fisher discriminant model of building damage was constructed. Finally， the method proposed in this paper was verified by taking the Luxian M6.0 earthquake as an example. The experimental results show that the prediction results of the proposed method closely align with the actual earthquake damage， with an accuracy rate of 80%， which proves that this method is highly accurate， reliable， and can accurately predict the earthquake damage degree of buildings.

Keywords： preliminary judgment from remote sensing; buildings; seismic damage prediction; Fisher discriminant method; Luxian earthquake

0 引言

歷次地震災(zāi)害經(jīng)驗表明，導(dǎo)致人員傷亡和經(jīng)濟損失最重要的原因是建構(gòu)筑物破壞及倒塌^［1］。因此，科學(xué)合理地預(yù)測建筑物在地震作用下的破壞程度，并將其直觀地展示出來，對震后快速了解災(zāi)區(qū)建筑物詳細震害情況，準確定位災(zāi)區(qū)位置，進行救援力量科學(xué)部署等應(yīng)急處置行為提供可靠的數(shù)據(jù)支持，從而有效減輕地震造成的損失^［2-3］。

目前，國內(nèi)外常用的建筑物震害預(yù)測模型主要有經(jīng)驗統(tǒng)計法^［4-8］、基于力學(xué)計算的理論計算法^［9-13］、半經(jīng)驗半統(tǒng)計法^［14-17］。經(jīng)驗統(tǒng)計法是基于震害調(diào)查資料，運用統(tǒng)計和判別總結(jié)得出的震害影響因素與破壞程度之間的關(guān)系，從而進行建筑震害預(yù)測的方法。該方法直觀、簡單，便于理解和應(yīng)用，但評估的結(jié)果是各類建筑物群體總體上的破壞情況，并不能反映建筑物破壞的空間分布?；诹W(xué)計算的理論計算法是運用地震動參數(shù)理論對預(yù)先經(jīng)過結(jié)構(gòu)簡化的力學(xué)模型進行非線性分析的一種方法^［18］，該方法可獲得單體建筑物破壞詳細情況。但在計算過程中，力學(xué)模型的構(gòu)建及其參數(shù)的確定，需要單體建筑物詳盡的結(jié)構(gòu)性參數(shù)及大量實驗的累積，前期需耗費大量的人力、物力及財力進行現(xiàn)場調(diào)查量測獲取數(shù)據(jù)，后期數(shù)據(jù)提煉及計算工作量也很大。半經(jīng)驗半理論方法是在基于歷史震害資料對同類事物統(tǒng)計分析的基礎(chǔ)上，通過理論方法，用嚴格的數(shù)學(xué)方法處理各變量之間關(guān)系。該方法綜合了經(jīng)驗法和理論法各自的優(yōu)點，拓寬了其應(yīng)用范圍^［19］。但模型的構(gòu)建難度非常大，很難同時滿足不同結(jié)構(gòu)類型建筑物的要求，普適性差。因此，需探索一種既能充分利用已有的地震現(xiàn)場調(diào)查資料，又能夠簡潔高效地開展評估，且評估結(jié)果能準確反映建筑物破壞情況的方法。

Fisher判別法是根據(jù)已有典型樣本資料，基于選取的判別影響因素提取指標統(tǒng)計值，對待測樣本屬性判別預(yù)測的一種統(tǒng)計分析方法^［20-21］。該方法現(xiàn)已廣泛用于人臉識別^［22-23］、故障診斷^［24-25］、煤礦工程^{［20-21，26］}、水文地質(zhì)^［27］等領(lǐng)域，這一方法的優(yōu)勢在于對原始樣本整體分布無特殊要求，且能全面考慮影響判別的各種因素，構(gòu)造簡單，預(yù)測精度也較高?；诖?，本文將Fisher判別方法應(yīng)用于建筑物震害預(yù)測，提出了一種建筑物快速震害預(yù)測方法。該方法通過歸納總結(jié)影響建筑物震害的主要因素，并選取了震級、震中距、場地條件及建筑物抗震能力4個震害因子作為判別因子，在此基礎(chǔ)上構(gòu)建了Fisher判別函數(shù)，建立了適用于不同結(jié)構(gòu)類型建筑的快速便捷的建筑物震害預(yù)測方法。

1 建筑物震害影響因素分析

建筑物震害是指某一地區(qū)建筑物在遭受某一強度地震作用下造成的破壞情況，其震害等級是由多種因素共同決定的，科學(xué)地選取震害影響因素是準確評估建筑物震害等級的前提。從建筑物破壞的成因來看，可將其分成兩方面：一是建筑物所在場地遭受地震作用的強度，與地震能量大小、工程場地條件等相關(guān);二是建筑物抗震能力的強弱，受建筑結(jié)構(gòu)類型、設(shè)防標準、是否加固、建造年代等諸多因素影響。因此，本文主要從這兩方面來選擇并確定建筑物震害影響因素。

1.1 地震強度

地震強度大小既可用震級來表示，其反映地震釋放能量的多少，也可用烈度表示，反映地震影響和破壞的程度。它不僅與本次地震的釋放能量（即震級）有關(guān)，而且和震源深度、距離震中的距離（震中距）、工程地質(zhì)條件等因素有關(guān)。

1.1.1 震級

地震震級是造成建筑物破壞的一個重要因素，其大小取決于地震時所釋放出的能量多少。釋放的能量越大，其對應(yīng)的震級越大，反之，則對應(yīng)的震級越小。一般而言，地震震級越大，在同一位置的建筑物遭受的破壞程度也就越大。

1.1.2 震中距

距離是影響建筑物震害的另一個重要因素，震中距直接影響建筑物遭受地震作用的強度大小，距離震中近的和距離震中遠的建筑物在地震中造成的破壞程度具有顯著差異。

1.1.3 場地條件

通過對歷史震害資料進行分析歸納可以發(fā)現(xiàn)：在不同的場地條件上，建筑物在遭受相同地震作用時的破壞程度是不同的^［28］。通常來說，在軟弱的地基上，弱性結(jié)構(gòu)易遭到破壞而剛性結(jié)構(gòu)不易受損;而對于堅硬基礎(chǔ)而言，則恰恰相反。目前，場地分類主要是根據(jù)覆蓋土層厚度和場地剪切波速（或場地土類型）兩個指標綜合確定的，采用的方法主要包括地形坡度法、地質(zhì)地貌、譜比法等。本文選擇Wald等^［29］提出的地形坡度法，利用DEM數(shù)據(jù)提取坡度，然后根據(jù)地形坡度、v_S30和場地類型的對應(yīng)關(guān)系（表1），對場地類型進行劃分，從而得到場地數(shù)據(jù)。

1.2 建筑物抗震能力

歷史震害統(tǒng)計數(shù)據(jù)表明，同一地區(qū)建筑物在承受相同地震作用強度時，致使其出現(xiàn)不同破壞程度的原因是房屋本身抗震能力的差異。房屋抗震能力強弱受諸多因素影響，主要包括結(jié)構(gòu)類型、設(shè)防標準、建造年代、建筑物用途、結(jié)構(gòu)層數(shù)、材料強度、砌筑類型等等。這些因素對數(shù)據(jù)的詳細程度要求高，實地調(diào)研工作量大，獲取難度高，適用于小范圍區(qū)域，且影響因素是否考慮齊全，數(shù)據(jù)間是否存在冗余等問題仍需考慮^［30］。因此，為適用于各類建筑物，并考慮到數(shù)據(jù)的可獲取性及完整性，盡可能選用最少的因子來解釋變量，綜合分析下選擇建筑物抗震能力這個單因素綜合指標。

本文根據(jù)建筑物在遙感影像上的光譜特征、紋理特征、幾何特征、空間位置特征等對建筑物抗震能力進行描述，并結(jié)合實地調(diào)研獲取的區(qū)域房屋結(jié)構(gòu)特點，建立區(qū)域房屋遙感影像解譯標志，按照建筑物抗震能力良好、建筑物抗震能力一般及建筑物抗震能力較差三個等級，判別房屋的抗震性能。（1）建筑物抗震能力好指經(jīng)過正規(guī)抗震設(shè)計的設(shè)防磚混、框架、框架-剪力墻及木結(jié)構(gòu)房屋。在影像上表現(xiàn)為屋頂為平頂，色調(diào)均勻偏亮，平面呈方形、T形、Y形、十字形等，形狀多規(guī)則且對稱，單體面積大且伴隨有較長陰影，單體呈點狀分布，多個單體組合呈簇狀或團狀分布，多分布在市區(qū)、各城鎮(zhèn)中心臨街、臨路位置。（2）建筑物抗震能力一般指未經(jīng)過正規(guī)抗震設(shè)防設(shè)計的磚混結(jié)構(gòu)房屋。在影像上表現(xiàn)為屋頂一般為平頂，色調(diào)不均勻且亮度偏暗，平面多為長方形、偶有L形或井字形，主要分布在老舊的城區(qū)、城鄉(xiāng)接合部及農(nóng)村部分道路旁邊。（3）建筑物抗震能力較差指老舊的磚混（預(yù)制板）、磚木及土木結(jié)構(gòu)房屋等。表現(xiàn)為屋頂多為人字坡頂，屋頂顏色多為深灰色或灰色，單體面積小且無陰影或陰影特征不明顯，獨立分布或零星分布在農(nóng)村偏遠地區(qū)，或房屋緊密相連無法區(qū)分邊界，呈塊狀分布。

2 基于Fisher判別的建筑物震害預(yù)測方法

2.1 Fisher判別法基本原理

Fisher判別法由費歇爾于1936年提出，屬于多元統(tǒng)計方法中判別分析方法的一種，根據(jù)訓(xùn)練樣例的多個特征確定待測樣本歸屬類別的一種監(jiān)督分類方法^［31-32］。其基本思想是通過尋求一個適當?shù)耐队胺较?，將高維特征空間數(shù)據(jù)投影到低維空間上^［33-34］，根據(jù)同一類樣本間離散度最小、不同類樣本間離散度最大的原則構(gòu)建判別分析函數(shù)，根據(jù)構(gòu)建的判別函數(shù)，以確定待測樣本歸屬類別的一種方法。

若假定樣本類別為G_i（i=1，2，…，n），分別從各類別中隨機抽取樣品數(shù)為k_i，判別函數(shù)為u（x）=a′x，對于任意給定的一個樣本x∈G_i的情況下，樣本均值E_i、樣本方差D_i、組內(nèi)距B₀、組間距C₀計算如下所示：

E_i=E（y|G_i）=E（a′x|G_i）=a′u_i， i=1，2，…，n（1）

D_i=D（y|G_i）=D（a′x|G_i）=a′v_ia， i=1，2，…，n（2）

B₀=∑n/i=1∑k_i/α=1（a′xⁱ_α-a′u_i）²=

a′∑n/i=1∑k_i/α=1（xⁱ_α-u_i）（xⁱ_α-u_i）′a=a′Ba（3）

C₀=∑n/i=1k_i（a′uⁱ-a′u）=

a′∑n/i=1k_i（u_i-u）（u_i-u）′a=a′Ca（4）

式中：u_i表示類別G_i的樣本均值，記u_i=1/k_i∑k_i/α=1xⁱ_α;v_i表示類別G_i的樣本協(xié)方差；u為總體均值，記u=1/5∑n/i=1∑k_i/α=1xⁱ_α；B為組內(nèi)離差矩陣，記B=∑n/i=1∑k_i/α=1（xⁱ_α-u_i）（xⁱ_α-u_i）′；C為組間離差矩陣，記C=∑n/i=1k_i（u_i-u）（u_i-u）′。

為達到較好的分類效果，應(yīng)盡可能使組內(nèi)樣本投影點分布集中、間距小，組間樣本投影點分布稀疏、間隔距離遠，即λ=C₀/B₀=∑n/i=1k_i（a′u_i-a′u）²/∑n/i=1∑k_i/α=1（a′xⁱ_α-a′u_i）²=a′Ca/a′Ba達到最大。根據(jù)極值存在的必要條件，即對其求導(dǎo)?λ/?α=0，從而轉(zhuǎn)換成求解方程|B-λC|=0的最大特征值，而a則對應(yīng)最大特征值的特征向量，從而得到判別函數(shù)。建立判別函數(shù)后，計算待測樣本的判別函數(shù)值，并與判別界值進行比較，進而判斷所屬類別。

2.2 建筑物震害Fisher判別模型

常用的建筑物震害預(yù)測方法是在選取歷史震害案例的基礎(chǔ)上，基于建筑物的震害影響因子，采用多元統(tǒng)計分析方法對未知建筑物破壞等級做出準確判斷。其基本原理與Fisher判別法一致，因此，本文將Fisher判別法應(yīng)用到建筑物震害預(yù)測中。首先，將選取的建筑物震害影響因子量化后作為震害特征;然后，根據(jù)建筑結(jié)構(gòu)的5種破壞狀態(tài)將樣本歸屬類別劃分成5類，即基本完好、輕微破壞、中等破壞、嚴重破壞及毀壞^［35］；最后，基于獲取的具有典型破壞特點的建筑物樣本進行訓(xùn)練，從而構(gòu)建建筑物震害Fisher判別模型（圖1）。

根據(jù)建筑物震害類別設(shè)5個母體類別為：G₁、G₂、G₃、G₄、G₅，并分別從各類別中隨機抽取樣品數(shù)為k₁，k₂，k₃，k₄，k₅，令k=k₁+k₂+k₃+k₄+k₅，則抽樣的樣品總數(shù)為k。將上一節(jié)確定的建筑物震害影響因素，即建筑物抗震能力x₁、場地類型x₂、震中距x₃和震級x₄作為判別因子，則xⁱ_α=（xⁱ₁，xⁱ₂，xⁱ₃，xⁱ₄）表示第i個總體的第α個樣本的觀測向量。根據(jù)式（1）～（4）計算各類別樣本均值E_i、樣本方差D_i、組內(nèi)距B₀、組間距C₀，進而求得判別函數(shù)，依據(jù)判別規(guī)則來判斷待測樣本的類別。

2.3 判別效果檢驗

在判別函數(shù)構(gòu)建完成后，還需要對其判別效果進行檢驗，從而驗證建立的判別函數(shù)是否具有合理性且符合實際。通常采用回代的方式進行檢驗，并以誤判率進行衡量^［36］，即誤判樣本數(shù)量與樣本總數(shù)之間的比值，具體計算如式（5）。如果計算得到的誤判率越小，表明該方法判別效果越好，反之，則說明該方法判別效果差。

β=t/k×100% （5）

3 以四川瀘定縣6.0級地震為例

本文以2021年9月16日四川瀘定縣6.0級地震為例，此次地震的最高烈度為Ⅷ度，等震線長軸呈NWW走向，震害波及四川省5個縣（市、區(qū)），51個鄉(xiāng)鎮(zhèn)（街道）;重慶市1個區(qū)，10個鄉(xiāng)鎮(zhèn)（街道）。瀘定縣地震發(fā)生后，筆者通過整理地震現(xiàn)場調(diào)查報告收集了80組建筑物震害數(shù)據(jù)作為樣本，在此基礎(chǔ)上對本文提出的方法進行驗證。

3.1 數(shù)據(jù)預(yù)處理

根據(jù)前述選取的建筑物震害影響因素，分析本文需要的數(shù)據(jù)，主要包括：震中位置、樣本點、航天雷達地形測繪（SRTM）數(shù)據(jù)、2020年四川省地理國情普查數(shù)據(jù)及高分遙感影像。

首先，對獲取的SRTM數(shù)據(jù)進行拼接、裁剪、投影轉(zhuǎn)換等操作獲取研究區(qū)（本文特指烈度為Ⅵ范圍內(nèi)）DEM數(shù)據(jù)。運用ArcGIS中的Slope命令，坡度單位選擇度，獲取區(qū)域地形坡度，依據(jù)中國構(gòu)造活躍區(qū)及穩(wěn)定大陸劃分結(jié)果及1.1.3小節(jié)中的地形坡度與v_S30的對應(yīng)關(guān)系，得到研究區(qū)范圍內(nèi)的場地分類結(jié)果，如圖2所示。從圖2中可以發(fā)現(xiàn)，場地類型共有四類，其中D類場地的面積占比最大，B類場地次之，C類和E類場地面積較少。從地貌上看，震中位于四川盆地的川中丘陵與川東平行嶺谷的過渡地帶，處在華鎣山脈向西南方，向延伸的分支及余脈之間的寬谷中，震中西北側(cè)為螺觀山，東南側(cè)是古佛山（薄刀嶺），地勢相對陡峭，其余地區(qū)均較平坦。表現(xiàn)在場地分類圖上，螺觀山和古佛山附近的場地類型與其余地區(qū)不相同，驗證了本文場地分類結(jié)果的正準確性。其次，基于2020年四川省地理國情普查地表覆蓋分類代碼提取建筑物數(shù)據(jù)，結(jié)合高分遙感影像，根據(jù)1.2節(jié)構(gòu)建的建筑物抗震能力解譯標志，通過解譯得到建筑物抗震能力遙感初判圖（圖3）。

3.2 Fisher判別模型建立

基于樣本點空間位置，結(jié)合經(jīng)數(shù)據(jù)預(yù)處理后的建筑物震害影響因素，運用ArcGIS空間疊加分析工具，分析得到各樣本點各影響因素值。通過整理歸納后得到樣本統(tǒng)計數(shù)據(jù)（表2）。利用ArcGIS的統(tǒng)計分析工具，從樣本中隨機抽取80%的樣本（64組）作為訓(xùn)練樣本，剩余20%（16組）作為驗證樣本（圖4）。

為了構(gòu)建判別模型，需要將建筑物震害影響因素定量化，根據(jù)建筑物破壞程度劃分為5類：1代表基本完好，2代表輕微破壞，3代表中等破壞，4代表嚴重破壞，5代表毀壞。建筑物抗震能力分為3類：1代表抗震能力好，2代表抗震能力一般，3代表抗震能力差;場地類型分為4類：4表示B類場地，3表示C類場地，2表示D類場地，1表示E類場地。根據(jù)上節(jié)構(gòu)建的建筑物震害Fisher判別模型，將處理完后的64組建筑物震害實例作為訓(xùn)練樣本建立Fisher判別模型，通過訓(xùn)練、學(xué)習(xí)及計算后得到各Fisher判別函數(shù)的系數(shù)，如表3所列。在此基礎(chǔ)上，建立的判別函數(shù)分別為：

y₁=1.858x₁+0.355x₂-1.194x₃-2.510 ;

y₂=1.209x₁-0.093x₂+0.958x₃-3.569;

y₃=-0.232x₁+1.349x₂-0.064x₃-3.109。

從表4中可看出：第一個判別函數(shù)方差貢獻率為75.5%，表明利用該函數(shù)可以對大部分樣本進行所屬類別判定。當利用第一個判別函數(shù)無法進行判斷時，可以利用第二個判別函數(shù)，聯(lián)合利用前2個判別通過計算可得到各判別函數(shù)特征值（表4）。前兩個函數(shù)判別時，累積方差貢獻率已達96%以上，已基本能對所有樣本類別進行判別。當加入第3個判別函數(shù)進行判別時，累積方差貢獻率已達100%，表明運用這3個函數(shù)能對所有樣本類別進行判別。

3.3 判別效果檢驗

為檢查所求判別函數(shù)的準確度，以回代的方式進行檢驗，把參與訓(xùn)練的64組樣本代入已構(gòu)建的判別函數(shù)中，依據(jù)判別準則對其所屬類別進行判斷，并將預(yù)測結(jié)果與實際震害進行對比，再利用式（5）計算其誤判率，結(jié)果如表5所列。從表中可以看出，基本完好的誤判率為11.8%，輕微破壞的為25%，中等破壞的為36.4%，嚴重破壞的為7.7%，其中輕微破壞與中等破壞誤判率偏高。從整體來看，有12組訓(xùn)練樣本出現(xiàn)預(yù)判結(jié)果與實際震害情況不符，誤判率為18.75%，表明模型整體預(yù)測精度較高，證明該方法可以較好地對建筑物震害情況進行預(yù)測。另外，從表中可以發(fā)現(xiàn)判斷結(jié)果與實際不符的樣本基本都是預(yù)判等級偏高的，其原因可能在于建筑抗震能力分類比較粗，同一類抗震能力的建筑物間還是存在較大差異，如均為抗震能力基本達標的設(shè)防磚混結(jié)構(gòu)與框架-剪力墻結(jié)構(gòu)。設(shè)防磚混結(jié)構(gòu)的抗震設(shè)防標準為本地設(shè)防烈度，抗震能力良好，而框架-剪力墻結(jié)構(gòu)其抗震設(shè)防標準多高于本地設(shè)防烈度，抗震性能明顯高于設(shè)防磚混，在面對相同地震強度時遭受的破壞程度更輕。

3.4 模型精度驗證

為進一步驗證建筑物震害Fisher判別模型的準確性，利用未參與訓(xùn)練的16組樣本進行驗證，并以專題圖的形式展示預(yù)測結(jié)果與實際震害的對比情況（圖5）。從圖中可以看出：驗證樣本的預(yù)測結(jié)果與實際震害情況基本一致，只有3組樣本預(yù)判結(jié)果有誤，且存在預(yù)判等級均偏高的情況，誤判率為18.75%，準確率達到81.25%。以上驗證表明該模型具有較高的可信度，能對建筑物震害等級進行科學(xué)預(yù)測。

此外，筆者還收集到2022年9月5日四川瀘定6.8級地震的81組建筑物震害數(shù)據(jù)，進一步驗證本文模型的科學(xué)性，具體結(jié)果如表6所列。從整體來看，有61組訓(xùn)練樣本出現(xiàn)預(yù)判結(jié)果與實際震害情況一致，準確率為75.3%，表明模型整體預(yù)測精度較高，證明該方法可以較好地對建筑物震害情況進行預(yù)測。

4 結(jié)論

為震后快速了解災(zāi)區(qū)建筑物詳細震害情況，精準定位受災(zāi)位置，進行救援力量科學(xué)部署等應(yīng)急處置提供可靠的數(shù)據(jù)支撐，本文提出了一種基于遙感初判的建筑物震害Fisher判別法，并結(jié)合四川瀘定縣6.0級地震案例進行分析，得到結(jié)論如下：

（1） 為適用于各類建筑物震害預(yù)測，從地震強度和建筑物抗震能力兩方面分析影響建筑物震害影響因素，并結(jié)合數(shù)據(jù)的可獲取性及完整性，綜合選擇震級、震中距、場地條件及建筑物抗震能力4個指標作為判別因子。

（2） 基于判別分析理論，構(gòu)建了基于Fisher判別的建筑物震害預(yù)測模型。利用建立的判別模型對四川瀘定縣6.0級地震收集的64組建筑物震害實例進行訓(xùn)練學(xué)習(xí)得到相應(yīng)的判別函數(shù)，并以回代的方式進行判別效果檢驗，計算其誤判率為18.75%，說明該方法預(yù)測精度較高。利用未參與訓(xùn)練的16組樣本進行驗證，預(yù)測結(jié)果與實際震害基本一致，準確率高達80%以上，證明該方法可靠度較高，能較準確地對建筑物震害程度進行科學(xué)預(yù)測。

（3） 本文提出的模型造成誤差的原因主要在于建筑抗震能力分類比較粗獷，同一類抗震能力的建筑物間存在較大差異，且采用地形坡度法進行場地分區(qū)，分區(qū)結(jié)果比較粗糙等。此外，考慮的影響因素也不全面，如破裂方向、破裂深度等尚未考慮。未來可從上述角度進一步優(yōu)化模型，以期提高建筑物震害判別的準確性。

參考文獻（References）

［1］劉晶晶，寧寶坤，呂瑞瑞，等.震后典型建筑物倒塌分類及救援特點分析［J］.震災(zāi)防御技術(shù)，2017，12（1）：220-229.

LIU Jingjing，NING Baokun，Lü Ruirui，et al.Classification of typical building collapse and analysis of rescue characteristics after the earthquake［J］.Technology for Earthquake Disaster Prevention，2017，12（1）：220-229.

［2］ZHAI W，SHEN H F，HUANG C L，et al.Building earthquake damage information extraction from a single post-earthquakePolSAR image［J］.Remote Sensing，2016，8（3）：171.

［3］YANG F，WEN X T，WANG X S，et al.A model study of building seismic damage information extraction and analysis on ground-based LiDAR data［J］.Advances in Civil Engineering，2021，2021：5542012.

［4］ALGERMISSEN S T，STEINBRUGGE K，LAGORIO H.Estimation of earthquake losses to buildings （except single family dwellings）［J］.Insurance Services Offices，1978.

［5］ROSSETTO T，ELNASHAI A.Derivation of vulnerability functions for European-type RC structures based on observational data［J］.Engineering Structures，2003，25（10）：1241-1263.

［6］GAUTAM D，F(xiàn)ABBROCINO G，SANTUCCI DE MAGISTRIS F.Derive empirical fragility functions for Nepali residential buildings［J］.Engineering Structures，2018，171：617-628.

［7］ALLALI S A，ABED M，MEBARKI A.Post-earthquake assessment of buildings damage using fuzzy logic［J］.Engineering Structures，2018，166：117-127.

［8］FOYTONG P，ORNTHAMMARATH T.Empirical seismic fragility functions based on field survey data after the 5 May 2014 Mae Lao （Northern Thailand） earthquake［J］.International Journal of Disaster Risk Reduction，2020，42：101344.

［9］HOSSEINI HASHEMI B，HASSANZADEH M.Study of a semi-rigid steel braced building damaged in the Bam earthquake［J］.Journal of Constructional Steel Research，2008，64（6）：704-721.

［10］NAGASHREE B K，RAVI K C M，VENKAT R D.A parametric study on seismic fragility analysis of RC buildings［J］.Earthquakes and Structures，2016，10（3）：629-643.

［11］于曉輝，呂大剛，范峰.基于易損性指數(shù)的鋼筋混凝土框架結(jié)構(gòu)地震損傷評估［J］.工程力學(xué)，2017，34（1）：69-75，100.

YU Xiaohui，Lü Dagang，F(xiàn)AN Feng.Seismic damage assessment of RC frame structures based on vulnerabilty index［J］.Engineering Mechanics，2017，34（1）：69-75，100.

［12］PALAGALA V Y，SINGHAL V.Structural score to quantify the vulnerability for quick seismic assessment of RC framed buildings in India［J］.Engineering Structures，2021，243：112659.

［13］呂碩碩，尹訓(xùn)強，王桂萱.基于OpenSees的某超高層建筑結(jié)構(gòu)地震響應(yīng)分析研究［J］.地震工程學(xué)報，2023，45（1）：34-42.

Lü Shuoshuo，YIN Xunqiang，WANG Guixuan.Seismic response analysis of a super high-rise building structure based on OpenSees［J］.China Earthquake Engineering Journal，2023，45（1）：34-42.

［14 ］PAMPANIN S，NIGEL PRIESTLEY M J，SRITHARAN S.Analytical modelling of the seismic behaviour of precast concrete frames designed with ductile connections［J］.Journal of Earthquake Engineering，2001，5（3）：329-367.

［15］劉章軍，葉燎原，潘文.基于模糊貝葉斯的震害預(yù)測［J］.昆明理工大學(xué)學(xué)報（理工版），2002，27（6）：108-111.

LIU Zhangjun，YE Liaoyuan，PAN Wen.Prediction of seismic disasters based on fuzzy-Bayes theory［J］.Journal of Kunming University of Science and Technology，2002，27（6）：108-111.

［16］韓博，熊琛，許鎮(zhèn)，等.城市區(qū)域建筑物震害預(yù)測剪切層模型及其參數(shù)確定方法［J］.工程力學(xué)，2014，31（增刊1）：73-78，85.

HAN Bo，XIONG Chen，XU Zhen，et al.A multi-story concentrated-mass shear model for urban regional seismic damage prediction and associated parameter determination［J］.Engineering Mechanics，2014，31（Suppl01）：73-78，85.

［17］MOYA L，MAS E，KOSHIMURA S，et al.Synthetic building damage scenarios using empirical fragility functions：a case study of the 2016 Kumamoto earthquake［J］.International Journal of Disaster Risk Reduction，2018，31：76-84.

［18］崔玉紅，邱虎，聶永安，等.國內(nèi)外單體建筑物震害預(yù)測方法研究述評［J］.地震研究，2001，24（2）：175-182.

CUI Yuhong，QIU Hu，NIE Yongan，et al.Research review of predictive method about single building in the world［J］.Journal of Seismological Research，2001，24（2）：175-182.

［19］張桂欣，孫柏濤.基于模糊層次分析的建筑物單體震害預(yù)測方法研究［J］.工程力學(xué)，2018，35（12）：185-193，202.

ZHANG Guixin，SUN Baitao.Seismic damage prediction for a single building based on a fuzzy analytical hierarchy approach［J］.Engineering Mechanics，2018，35（12）：185-193，202.

［20］宮鳳強，李夕兵.礦區(qū)立井井筒非采動破裂的距離判別分析預(yù)測法［J］.煤炭學(xué)報，2007，32（7）：700-704.

GONG Fengqiang，LI Xibing.A distance discriminant analysis method of forecast for shaft-lining non-mining fracture of mine［J］.Journal of China Coal Society，2007，32（7）：700-704.

［21］董隴軍，李夕兵，白云飛.急傾斜煤層頂煤可放性分類預(yù)測的Fisher判別分析模型及應(yīng)用［J］.煤炭學(xué)報，2009，34（1）：58-63.

DONG Longjun，LI Xibing，BAI Yunfei.A Fisher discriminant analysis model for classifying top coal cavability of the steep seam［J］.Journal of China Coal Society，2009，34（1）：58-63.

［22］HU W P，HU H F.Heterogeneous face recognition based on modality-independent kernel Fisher discriminant analysis joint sparse auto-encoder［J］.Electronics Letters，2016，52（21）：1753-1755.

［23］GU X Q，NI T G，WANG W B，et al.Cross-domain transfer person re-identification via topology properties preserved local fisher discriminant analysis［J］.Journal of Ambient Intelligence and Humanized Computing，2020（4）：1-11.

［24］ZHU Z B，SONG Z H.A novel fault diagnosis system using pattern classification on kernel FDA subspace［J］.Expert Systems with Applications，2011，38（6）：6895-6905.

［25］MD NOR N，HUSSAIN M A，CHE HASSAN C R.Multi-scale kernel Fisher discriminant analysis with adaptive neuro-fuzzy inference system （ANFIS） in fault detection and diagnosis framework for chemical process systems［J］.Neural Computing and Applications，2020，32（13）：9283-9297.

［26］CHEN L，WANG E Y，F(xiàn)ENG J J，et al.Hazard prediction of coal and gas outburst based on Fisher discriminant analysis［J］.Geomechanics and Engineering，2017，13（5）：861-879.

［27］BI Y S，WU J W，ZHAI X R，et al.Discriminant analysis of mine water inrush sources with multi-aquifer based on multivariate statistical analysis［J］.Environmental Earth Sciences，2021，80（4）：1176-1188.

［28］趙密，賈智富，昝子卉，等.不同場地條件下地連墻對地鐵車站結(jié)構(gòu)地震動力響應(yīng)的影響［J］.地震工程學(xué)報，2023，45（2）：270-278，295.

ZHAO Mi，JIA Zhifu，ZAN Zhihui，et al.Influence of the diaphragm wall on the dynamic responses of subway station structures under different site conditions［J］.China Earthquake Engineering Journal，2023，45（2）：270-278，295.

［29］WALD D J，ALLEN T I.Topographic slope as a proxy for seismic site conditions and amplification［J］.Bulletin of the Seismological Society of America，2007，97（5）：1379-1395.

［30］張朝陽，李雪，劉珠妹.遙感自動提取技術(shù)在房屋抗震調(diào)查中的應(yīng)用［J］.地震工程學(xué)報，2023，45（6）：1478-1484.

ZHANG Zhaoyang，LI Xue，LIU Zhunmei.Application of the automatic extraction method of remote sensing images to the seismic capacity investigation of buildings［J］.China Earthquake Engineering Journal，2023，45（6）：1478-1484.

［31］錢兆明，任高峰，褚夫蛟，等.基于PCA法和Fisher判別分析法的巖體質(zhì)量等級分類［J］.巖土力學(xué)，2016，37（增刊2）：427-432，441.

QIAN Zhaoming，REN Gaofeng，CHU Fujiao，et al.Classification of rock mass quality grades based on PCA and Fisher discriminant analysis［J］.Rock and Soil Mechanics，2016，37（Suppl02）：427-432，441.

［32］ZHANG W Q，WANG Z Y，ZHU X X，et al.A risk assessment of a water-sand inrush during coal mining under a loose aquifer based on a factor analysis and the Fisher model［J］.Journal of Hydrologic Engineering，2020，25（8）：04020033.

［33］胡漢華，劉征，李孜軍，等.硫化礦石自燃傾向性等級分類的Fisher判別分析法［J］.煤炭學(xué)報，2010，35（10）：1674-1679.

HU Hanhua，LIU Zheng，LI Zijun，et al.Fisher discriminant analysis to the classification of spontaneous combustion tendency grade of sulphide ores［J］.Journal of China Coal Society，2010，35（10）：1674-1679.

［34］趙麗娜.引入權(quán)重因子的Fisher判別法的研究［J］.哈爾濱師范大學(xué)自然科學(xué)學(xué)報，2012，28（5）：24-26.

ZHAO Li'na.The study on introduction weighting factor of Fisher discriminant［J］.Natural Science Journal of Harbin Normal University，2012，28（5）：24-26.

［35］尹之潛.地震災(zāi)害及損失預(yù)測方法［M］.北京：地震出版社，1995.

YIN Zhiqian.Earthquake disaster and loss prediction method［M］.Beijing：Seismological Press，1995.

［36］DONG L J，WESSELOO J，POTVIN Y，et al.Discrimination of mine seismic events and blasts using the Fisher classifier，naive Bayesian classifier and logistic regression［J］.Rock Mechanics and Rock Engineering，2016，49（1）：183-211.

（本文編輯：任 棟）