方前程，商 麗，商擁輝，2，宋 譯

爆破振動(dòng)誘發(fā)民房結(jié)構(gòu)損傷識(shí)別的隨機(jī)森林模型＊

方前程1，商 麗1，商擁輝1，2，宋 譯3

（1.黃淮學(xué)院建筑工程學(xué)院，河南 駐馬店463000；2.中南大學(xué)土木工程學(xué)院，湖南 長(zhǎng)沙410075；3.湖南科技大學(xué)能源與安全工程學(xué)院，湖南 湘潭411201）

為快速、準(zhǔn)確地評(píng)價(jià)爆破振動(dòng)誘發(fā)民房結(jié)構(gòu)損傷效應(yīng)，借鑒隨機(jī)森林理論并結(jié)合工程實(shí)際，建立露采爆破振動(dòng)誘發(fā)民房結(jié)構(gòu)損傷識(shí)別的隨機(jī)森林模型；綜合考慮爆破參數(shù)、爆破振動(dòng)特征參量及房屋結(jié)構(gòu)動(dòng)力特性等因素，選取質(zhì)點(diǎn)峰值振動(dòng)速度、主頻率、主頻率持續(xù)時(shí)間、段藥量、爆心距、施工質(zhì)量參數(shù)、場(chǎng)地條件參數(shù)、屋蓋形式參數(shù)、磚墻面積率、民房高度、灰縫強(qiáng)度和圈梁構(gòu)造柱參數(shù)等12個(gè)影響因素作為模型輸入，將磚混結(jié)構(gòu)建筑物的損害等級(jí)作為模型輸出；基于多分類器集成的思想，以108組爆破振動(dòng)實(shí)測(cè)數(shù)據(jù)作為學(xué)習(xí)樣本進(jìn)行訓(xùn)練，建模過(guò)程中由多個(gè)決策樹(shù)集成隨機(jī)森林、用投票的方式實(shí)現(xiàn)對(duì)民房結(jié)構(gòu)損傷有效識(shí)別；用12組現(xiàn)場(chǎng)數(shù)據(jù)驗(yàn)證模型的有效性；在對(duì)樣本分類的同時(shí)，計(jì)算預(yù)測(cè)變量的重要性值，發(fā)現(xiàn)質(zhì)點(diǎn)峰值振動(dòng)速度為最重要的評(píng)價(jià)指標(biāo)，其后依次為爆心距，主頻率持續(xù)時(shí)間，主頻率，圈梁構(gòu)造柱參數(shù)，灰縫強(qiáng)度，屋蓋形式參數(shù)，民房高度，段藥量，施工質(zhì)量參數(shù)，磚墻面積率和場(chǎng)地條件參數(shù)。研究結(jié)果表明：隨機(jī)森林模型預(yù)測(cè)結(jié)果學(xué)習(xí)樣本準(zhǔn)確度是87.97%，而測(cè)試集準(zhǔn)確度是91.67%，與實(shí)際情況吻合較好，預(yù)測(cè)精度較高。

爆破振動(dòng)；民房結(jié)構(gòu)損傷；隨機(jī)森林；質(zhì)點(diǎn)峰值振動(dòng)速度；預(yù)測(cè)

爆破振動(dòng)產(chǎn)生的負(fù)面效應(yīng)是礦山巖土開(kāi)挖爆破難以回避的難題，尤其對(duì)近鄰建構(gòu)筑物的影響。而隨著經(jīng)濟(jì)發(fā)展和我國(guó)人均生活水平的提高，人口逐漸從農(nóng)村向城市轉(zhuǎn)移，因此，人們對(duì)民用住房和商用建筑的需求越來(lái)越大，尤其是在繁華地段，往往建筑物的密度比較高，所以這就要求相鄰的工程不能對(duì)已有的建筑物構(gòu)成安全影響，因此需要建立更高的安全量化標(biāo)準(zhǔn)。

在露天開(kāi)采的過(guò)程中，通常會(huì)采用爆破的方法，除了會(huì)對(duì)臺(tái)階的穩(wěn)定性造成一定的影響之外；爆破產(chǎn)生的震動(dòng)還會(huì)影響到附近的建筑物，為此，學(xué)者們力求對(duì)爆破振動(dòng)誘發(fā)建筑物損害機(jī)制進(jìn)行研究、對(duì)建筑物破壞機(jī)理進(jìn)行分析。近些年來(lái)，國(guó)內(nèi)外許多專家對(duì)爆破振動(dòng)誘發(fā)民房破壞的機(jī)理進(jìn)行了深入研究，從源頭和傳播途徑上對(duì)爆破震動(dòng)等級(jí)進(jìn)行了預(yù)報(bào)，防止爆破振動(dòng)引起建筑破壞，并劃分了爆破振動(dòng)等級(jí)，以上工作均取得了不錯(cuò)的成效。此外，國(guó)內(nèi)外相繼出臺(tái)了一系列爆破振動(dòng)安全判據(jù)、標(biāo)準(zhǔn)和規(guī)程等法律法規(guī)［1－3］，但是還沒(méi)有形成一套完整的關(guān)于爆破震動(dòng)誘發(fā)建筑物損害的論述體系，在實(shí)踐方面缺少相關(guān)指導(dǎo)，項(xiàng)目工程方面也缺少相關(guān)數(shù)據(jù)分析。因此理論與實(shí)際往往會(huì)存在一定程度的偏差，尤其是通過(guò)觀察大量的爆破工程實(shí)例后，發(fā)現(xiàn)某些爆破震動(dòng)已經(jīng)超出安全標(biāo)準(zhǔn)，卻并未對(duì)周圍的民房構(gòu)成損害，分析其原因，一方面是因?yàn)槟壳暗陌踩袚?jù)尚不完善，對(duì)于不同環(huán)境下的爆破振動(dòng)誘發(fā)民房破壞可能會(huì)產(chǎn)生差異，而安全判據(jù)也應(yīng)隨之變化；另一方面是，爆破振動(dòng)誘發(fā)建筑物損害程度受各種因素影響，而且誘因往往比較復(fù)雜多樣，誘因和損害程度之間也呈非線性關(guān)系；因此，判別準(zhǔn)則不具有普遍適用性。近些年來(lái)國(guó)內(nèi)外一些專家通過(guò)反應(yīng)譜法［4］、時(shí)程分析法［5］和波動(dòng)理論分析法［6］、判別分析法［7－9］、梯度提升機(jī)［10］、模糊神經(jīng)網(wǎng)絡(luò)［11－12］和支持向量機(jī)［13］等軟科學(xué)方法對(duì)其進(jìn)行評(píng)估，雖然這些方法對(duì)爆破振動(dòng)誘發(fā)風(fēng)險(xiǎn)評(píng)估產(chǎn)生了積極的影響，但是這些方法都不具有普遍適用性。具體表現(xiàn)在：（1）輸入?yún)⒘康碾S機(jī)性與模糊性，各參量間相互關(guān)系（包括定性與定量）尚不明晰，以及統(tǒng)計(jì)資料不完整性；（2）個(gè)模型型本身固有不足與局限性，如神經(jīng)網(wǎng)絡(luò)方法存在隱含層確定具有主觀性，且易陷入局部最優(yōu)和收斂速度慢等不足；（3）支持向量機(jī)核函數(shù)較多且難以確定合適的核函數(shù)參數(shù)等。

隨機(jī)森林（random forest，RF）的算法是由 Leo Breiman和 Adele Cutler發(fā)展并推論出的［14］。作為一種包含多個(gè)決策樹(shù)的分類器，隨機(jī)森林算法簡(jiǎn)潔且具有較好的解釋能力，近些年來(lái)被應(yīng)用于自然科學(xué)和工程技術(shù)的各個(gè)領(lǐng)域。鑒于誘發(fā)民房破壞因素具有多樣不一的特點(diǎn)，部分誘因難以去清晰的量化，且評(píng)價(jià)對(duì)象不是單一固定的，本文中從開(kāi)挖區(qū)爆破振動(dòng)對(duì)附近房屋影響的角度出發(fā)，構(gòu)建露采爆破震動(dòng)誘發(fā)建筑物損害的RF判別模型，以分析民房破壞等級(jí)和誘發(fā)民房破壞因素之間的關(guān)系，以便對(duì)可能產(chǎn)生的爆破振動(dòng)進(jìn)行預(yù)報(bào)，為工程操作人員提供一個(gè)可參考的依據(jù)。

1 隨機(jī)森林分類原理

1.1 隨機(jī)森林方法

隨機(jī)森林（random forest，RF）方法本質(zhì)上將隸屬于集成學(xué)習(xí)范疇，其算法步驟為［14－15］：（a）在原始樣本集L通過(guò)自助法（bootstrap）重采樣等核心技術(shù)形成ntree個(gè)自助樣本集L1，L2，…Lntree，其大小約為2L／3；（b）為自助樣本訓(xùn)練子集分別建立分類回歸決策樹(shù)（CART）［15－16］，共產(chǎn)生ntree棵決策樹(shù)形成隨機(jī)森林，這些決策樹(shù)生長(zhǎng)過(guò)程中無(wú)剪枝環(huán)節(jié)，但從M 個(gè)屬性中隨機(jī)選擇mtry（mtry≤M）個(gè)最優(yōu)屬性作為內(nèi)部節(jié)點(diǎn)進(jìn)行分支操作；（c）集合步驟（a）和（b）產(chǎn)生的ntree棵決策樹(shù)預(yù)測(cè)結(jié)果，采用投票方法決定新樣本歸屬，見(jiàn)圖1。

1.2 RF變量重要性

RF方法相對(duì)于其他分類器的明顯優(yōu)勢(shì)是能夠識(shí)別出每個(gè)變量的重要性值，它提供了2種基本的屬性重要度計(jì)算依據(jù)［14－15］：Gini importance參數(shù)和 Permutation importance參數(shù)。在隨機(jī)森林節(jié)點(diǎn)分裂過(guò)程中，本文中采用Gini importance參數(shù)（G）來(lái)衡量各節(jié)點(diǎn)的樣本純度：爆破振動(dòng)誘發(fā)民房結(jié)構(gòu)損傷數(shù)據(jù)樣本集L共有3類樣本，并且每類的樣本數(shù)分別是pj（j＝1，2，3），則有：

1.3 RF預(yù)測(cè)模型的評(píng)價(jià)指標(biāo)

常用的精度評(píng)價(jià)方法有混淆矩陣方法，分類準(zhǔn)確率和κ指數(shù)［16－17］。混淆矩陣是一個(gè)M×M 矩陣（M 為分類數(shù)），用于比較參照點(diǎn)和分類點(diǎn)，如表1所示。而κ指數(shù)多用于分類預(yù)測(cè)模型的精度評(píng)價(jià)，即從被評(píng)價(jià)分類比完全隨機(jī)分類產(chǎn)生錯(cuò)誤減少的比例角度表征模型精度，即：

式中：xii是i行i列上的值，xi＋是第i行之和，x＋i是第i列之和，N 是樣本總數(shù)，r是誤差矩陣行數(shù)。

根據(jù)J.Landis等的［18］建議，κ統(tǒng)計(jì)值與分類精度對(duì)應(yīng)關(guān)系如下：κ最低允許判別精度為0.4，否則模型精度低；當(dāng)κ＝0.60～0.80時(shí)，分類精度較好；當(dāng)κ＝0.80～0.10時(shí)，分類效果非常好。

2 爆破振動(dòng)誘發(fā)民房結(jié)構(gòu)損傷預(yù)測(cè)的RF模型及其應(yīng)用

2.1 確定RF模型輸入輸出參量

爆破振動(dòng)誘發(fā)民房結(jié)構(gòu)損傷主要受兩方面影響一方面來(lái)自于振動(dòng)源；包括段藥量、裝藥不耦合系數(shù)、炸藥的性質(zhì)、起爆的方式等；另一方面取決于民房的各種參數(shù)，包括構(gòu)建民房材料的性質(zhì)、民房的結(jié)構(gòu)、民房附近的地形條件等。除此之外，振動(dòng)源與民房的距離也會(huì)影響到振動(dòng)傳感器的輸出結(jié)果。

根據(jù)前人研究成果，并結(jié)合爆破載荷下民房結(jié)構(gòu)損傷機(jī)制［7－11，13］，本模型選取爆破振動(dòng)三大特征參量：質(zhì)點(diǎn)峰值振動(dòng)速度（vppv）、主頻率（ν）及其持續(xù)時(shí)間（Δt），爆源因素參量：最大段藥量（Qmax）和距離（R）；民房結(jié)構(gòu)特征參量：如施工質(zhì)量參數(shù)（Qc）、場(chǎng)地條件參數(shù)（Sc）、屋蓋形式參數(shù)（Rs）、磚墻面積率（K）、房屋高度（H）、灰縫強(qiáng)度（S）和圈梁構(gòu)造柱參數(shù)（Bcf）等12個(gè)主要影響民房結(jié)構(gòu)損傷輸入?yún)⒘浚渲衯ppv、ν、Δt、Qmax、R、S、K 和H 等參量可用實(shí)測(cè)值直接輸入；對(duì)于Qc、Sc、Rs、Bcf等4個(gè)自變量和破壞等級(jí)（V）因變量均為狀態(tài)參量，為了清晰度量，按照表2標(biāo)準(zhǔn)對(duì)其進(jìn)行量化［7－8，10］。

表1 M 元分類問(wèn)題混淆矩陣Table 1 M－ary classification confusion matrix

表2 狀態(tài)參量數(shù)據(jù)量化建議值Table 2Recommended value for quantified input and output parameters

文獻(xiàn)［7－11］中從民房結(jié)構(gòu)破壞特征入手，將民房結(jié)構(gòu)損傷程度分為3個(gè)類別：（1）基本完好，記為V1；（2）輕微損傷，記為V2；（3）嚴(yán)重破壞，記為V3，其損傷類別描述見(jiàn)圖2，隨機(jī)森林模型輸出向量為V＝［V1，V2，V3］。

2.2 仿真數(shù)據(jù)采集及數(shù)據(jù)歸一化處理

為驗(yàn)證本文中提出RF模型的有效性和可靠性，以J.Zhou等［10］和史秀志等［11］現(xiàn)場(chǎng)實(shí)測(cè)的120組建筑物采動(dòng)損害實(shí)例數(shù)據(jù)作為總體樣本，其中隸屬于基本完好，輕微損傷和嚴(yán)重破壞3類民房危害類別的樣本數(shù)分別為65（V1），23（V2）和32（V3），并將原始數(shù)據(jù)集隨機(jī)分為2組：其中108組（編號(hào) X1，X2，…，X108）作為RF模型訓(xùn)練樣本集進(jìn)行訓(xùn)練（表3），其余部分（編號(hào)C1，C2，…，C6）作為檢驗(yàn)集（見(jiàn)表4）。

表3 RF模型學(xué)習(xí)樣本及識(shí)別結(jié)果Table 3Training samples and identification results of RF model

表4 RF模型測(cè)試樣本及識(shí)別結(jié)果對(duì)比Table 4RF model test samples and recognition results

2.3 建立民房危害等級(jí)評(píng)估的RF模型

RF模型輸入向量為（vppv，ν，Δt，Qmax，R，S，K，H，Rs，Bcf，Qc，Sc），輸出層為3種損傷等級(jí)類別V：V1，V2和V3。建立映射V→F（vppv，ν，Δt，Qmax，R，S，K，H ，Rs，Bcf，Qc，Sc）。應(yīng)用 RF理論，借助 R軟件編制相關(guān)數(shù)據(jù)分析程序建立模型，對(duì)各個(gè)樣本進(jìn)行歸類判別（見(jiàn)表4）。本文中用Random Jungle對(duì)模擬數(shù)據(jù)集生長(zhǎng)出1 000棵樹(shù)，用10折交叉驗(yàn)證對(duì)訓(xùn)練集進(jìn)行遍歷尋優(yōu)計(jì)算，并根據(jù)最低誤差選取RF模型最優(yōu)參數(shù)mtry，如圖3所示，在分類的過(guò)程中可生成一個(gè)泛化誤差的內(nèi)部無(wú)偏估計(jì)值，當(dāng)mtry＝12時(shí)，該值為15.74%，10折交叉驗(yàn)證準(zhǔn)確率達(dá)87.12%，κ＝0.772。RF混淆矩陣顯示訓(xùn)練集分類誤差見(jiàn)表5。在對(duì)樣本分類的同時(shí)，利用RF模型對(duì)解釋變量重要性排序，從圖4可以觀察到vppv為最重要的評(píng)價(jià)指標(biāo)，其后依次為Bcf、R、ν、Δt、Qmax、S、Rs、H 、Qc、K 和Sc。

根據(jù)學(xué)習(xí)好的民房結(jié)構(gòu)損傷識(shí)別RF模型對(duì)12個(gè)測(cè)試樣本進(jìn)行識(shí)別，其混淆矩陣見(jiàn)表6，可見(jiàn)識(shí)別結(jié)果與實(shí)際狀態(tài)較為相符，準(zhǔn)確率達(dá)91.67%（11／12），κ高達(dá)0.875，并與神經(jīng)網(wǎng)絡(luò)方法［11］預(yù)測(cè)結(jié)果較為一致（見(jiàn)表4）。從而說(shuō)明構(gòu)建的爆破振動(dòng)誘發(fā)民房結(jié)構(gòu)損傷隨機(jī)森林模型泛化能力強(qiáng)，魯棒性好。

與神經(jīng)網(wǎng)絡(luò)方法相比［11］，RF方法優(yōu)勢(shì)在于：（1）無(wú)需過(guò)多的數(shù)據(jù)預(yù)處理，RF預(yù)測(cè)精度高且運(yùn)算速度高效；（2）RF可自動(dòng)辨識(shí)各輸入變量重要度；（3）RF受異常值和噪聲的影響較小；（4）能同時(shí)處理連續(xù)型變量和分類變量。與此同時(shí)，RF缺點(diǎn)在于：（1）對(duì)若干數(shù)據(jù)集易呈現(xiàn)過(guò)擬合狀態(tài)因缺少剪枝步驟；（2）對(duì)于不均衡數(shù)據(jù)集，RF算法傾向于觀測(cè)值較多的類別（如V1的記錄較多，而且V2和V3間的差距不大，預(yù)測(cè)值會(huì)傾向于V1）。在實(shí)際工程中，結(jié)合本文模型確定民房的損害程度，可及時(shí)對(duì)民房進(jìn)行補(bǔ)償或者維護(hù)，降低爆破振動(dòng)帶來(lái)的后續(xù)危害。另一方面為了使該模型更有說(shuō)服力，可更廣泛的收集爆破震動(dòng)資料，建立更詳細(xì)的數(shù)據(jù)庫(kù)，使該方法的可信度更高。

表6 隨機(jī)森林的混淆矩陣顯示測(cè)試集分類誤差Table 6Confusion matrix drawn from Random forest showing the classification error of test set

3 結(jié) 論

（1）綜合考慮振動(dòng)源的各種參數(shù)和民房的各種材料、強(qiáng)度性質(zhì)，選取12大影響民房采動(dòng)損害程度的因素（vppv，ν，Δt，Qmax，R，S，K，H，Rs，Bcf，Qc，Sc）作為 RF模型輸入，將民房損害等級(jí)作為 RF模型輸出，提出并建立了露采爆破振動(dòng)誘發(fā)民房結(jié)構(gòu)損傷的RF識(shí)別模型。（2）利用108組實(shí)測(cè)數(shù)據(jù)計(jì)算表明，利用上述指標(biāo)建立的RF分析模型對(duì)爆破振動(dòng)誘發(fā)民房結(jié)構(gòu)損傷精度高，且方法簡(jiǎn)單實(shí)用、結(jié)果準(zhǔn)確可靠。（3）計(jì)算預(yù)測(cè)變量的重要性值發(fā)現(xiàn)vppv為最重要的評(píng)價(jià)指標(biāo)，其后依次為Bcf、R、ν、Δt、Qmax、S、Rs、H 、Qc、K 和Sc。

［1］ Khandelwal M，Singh T N.Evaluation of blast－induced ground vibration predictors［J］.Soil Dynamics and Earth－quake Engineering，2007，27（2）：116－125.

［2］ Siskind D E，Stagg M S，Kopp J W，et al.Structure response and damage produced by ground vibration from surface mine blasting［R］.USBM Report of Investigation 8507，1980：73－74.

i［3］ 中華人民共和國(guó)國(guó)家標(biāo)準(zhǔn)編寫(xiě)組.爆破安全規(guī)程：GB6722－2014［S］.北京：中國(guó)標(biāo)準(zhǔn)出版社，2015.

［4］ 史秀志，周健，崔松，等.露天采礦爆破振動(dòng)對(duì)民房危害預(yù)測(cè)的DDA模型及應(yīng)用［J］.中南大學(xué)學(xué)報(bào)（自然科學(xué)版），2011，42（2）：441－448.Shi Xiuzhi，Zhou Jian，Cui Song，et al.Distance discriminant analysis model and its application for prediction residential house＇s damage against blasting vibration of open pit mining［J］.Journal of Central South University（Science and Technology），2011，42（2）：441－448.

［5］ 史秀志，周健，杜坤，等.爆破振動(dòng)對(duì)民房破壞效應(yīng)預(yù)測(cè)的 BDA 模型及應(yīng)用［J］.振動(dòng)與沖擊，2010，29（7）：60－65.Shi Xiuzhi，Zhou Jian，Du Kun，et al.BDA model for predicting destructive effect of blast vibration on housing［J］.Journal of Vibration and Shock，2010，29（7）：60－65.

［6］ 中國(guó)生，徐國(guó)元，熊正明.基于小波變換的爆破地震信號(hào)能量分析法的應(yīng)用研究［J］.爆炸與沖擊，2006，26（3）：222－227.Zhong Guosheng，Xu Guoyuan，Xiong Zhengming.Application research of the energy analysis method for blasting seismic signals based on wavelet transform［J］.Explosion and Shock Waves，2006，26（3）：222－227.

［7］ 史秀志，周健，崔松，等.露天采礦爆破振動(dòng)對(duì)民房危害預(yù)測(cè)的DDA模型及應(yīng)用［J］.中南大學(xué)學(xué)報(bào)（自然科學(xué)版），2011，42（2）：441－448.Shi Xiuzhi，Zhou Jian，Cui Song，et al.Distance discriminant analysis model and its application for prediction residential house’s damage against blasting vibration of open pit mining［J］.Journal of Central South University（Science and Technology），2011，42（2）：441－448.

［8］ 史秀志，周健，杜坤，等.爆破振動(dòng)對(duì)民房破壞效應(yīng)預(yù)測(cè)的 BDA 模型及應(yīng)用［J］.振動(dòng)與沖擊，2010，29（7）：60－65.Shi Xiuzhi，Zhou Jian，Du Kun，et al.BDA model for predicting destructive effect of blast vibration on housing［J］.Journal of Vibration and Shock，2010，29（7）：60－65.

［9］ 董隴軍，李夕兵，趙國(guó)彥，等.露天采礦爆破振動(dòng)對(duì)砌體結(jié)構(gòu)破壞效應(yīng)預(yù)測(cè)的Fisher判別模型及應(yīng)用［J］.巖石力學(xué)與工程學(xué)報(bào)，2009，28（4）：750－756.Dong Longjun，Li Xibing，Zhao Guoyan，et al.Fisher discriminant analysis model and its application to predicting destructive effect of masonry structure under blasting vibration of open－pit mine［J］.Chinese Journal of Rock Mechanics and Engineering，2009，28（4）：750－756.

［10］ Zhou J，Shi X，Li X.Utilizing gradient boosted machine for the prediction of damage to residential structures owing to blasting vibrations of open pit mining［J］.Journal of Vibration and Control，2015，22（19）.DOI：10.1177／1077546314568172.

［11］ 史秀志.爆破振動(dòng)信號(hào)時(shí)頻分析與爆破振動(dòng)預(yù)測(cè)研究［D］.長(zhǎng)沙：中南大學(xué)，2007.

［12］ Melih I，Mahmut Y，Hakan A.Prediction of ground vibrations resulting from the blasting operations in an openpit mine by adaptive neuro－fuzzy inference system［J］.Environmental Geology，2008，56（11）：97－107.

［13］ Shi X Z，Zhou J.Prediction residential house’s damage effect near openpit against blasting vibration based on svm with grid searching method／genetic algorithm［J］.Advanced Science Letters，2012，11（1）：238－243.

［14］ Breiman L.Random forests［J］.Machine Learning，2001，45：5－32.

［15］ 方匡南，吳見(jiàn)彬，朱建平等.隨機(jī)森林方法研究綜述［J］.統(tǒng)計(jì)與信息論壇，2011，26（3）：32－38.Fang Kuangnan，Wu Jianbin，Zhu Jianping，et al.A review of technologies on random forests［J］.Statistics ＆Information Forum，2011，26（3）：32－38.

［16］ 秦鋒，楊波，程澤凱.分類器性能評(píng)價(jià)標(biāo)準(zhǔn)研究［J］.計(jì)算機(jī)技術(shù)與發(fā)展，2006，16（10）：85－88.Qin Feng，Yang Bo，Chen Zekai.Research on measure criteria in evaluating classification performance［J］.Computer Technology and Development，2006，16（10）：85－88.

［17］ Congalton R G.A review of assessing the accuracy of classifications of remotely sensed data［J］.Working Pages，1991，119（2）：270－279.

［18］ Landis J，Koch G.The measurement of observer agreement for categorical data［J］.Biometrics，1977，33（1）：159－174.

Random forest model for identification of residential structure damage induced by blast vibration

Fang Qiancheng1，Shang Li1，Shang Yonghui1，2，Song Yi3
（1.Institute of Architecture and Engineering，Huanghuai University，Zhumadian 463000，Henan，China；2.School of Civil Engineering，Central South University，Changsha 410075，Hunan，China；3.School of Energy and Safty Enginerring，Hunan University of Science and Technology，Xiangtan 411201，Hunan，China）

In this work，aiming to the prediction speed and accuracy，we established a random forest model for residential structure damage induced by blast vibration identification on the basis of the random forest（RF）theory.Twelve indexes，i.e.peak particle velocity，dominant frequency，dominant frequency duration，maximum charge per delay，distance，gray joints intensity，rate of brick walls，height of housing，roof structures parameter，beam－column frames parameter，quality parameter of construction and site conditions parameters，were considered as the criterion indices for this kind of damage in the proposed model based on the of analysis of the characteristic parameters of blasting vibration and dynamic characteristics of the housing structure.108sets of vibration measured data were investigated to create an RF classifier.RF was a combination of tree predictors，and variable importance was measured by gini importance parameter when the forest grows.A random tree was a combination of decision trees，and each tree is generated depending on the values of random vectors sampled independently，with the same distribution for all trees in the forest.The Gini importance value shows that the peak particle velocity is the most important discrimination indicator，followed by the distance，the dominant frequency duration，the dominant frequency，the beam－column frames parameter，the gray joints intensity，the roof structures parameter，the height of housing，the maximum charge per delay，the quality parameter of construction，the rate of brick walls and the site conditions parameters.Another twelve groups of residential structure damage instances were tested as forecast samples，and the predicted results were identical with the actual situation.Engineering practices indicate that the accuracy of the RF method of learning samples is 87.97%，and the accuracy of the test samples is 91.7%，effectively verifying and supplementing the existing methods for evaluating residential structure damage induced by blast vibration.

blasting vibration；residential structure damage；randomforest；peak particle velocity；prediction

O381 國(guó)標(biāo)學(xué)科代碼：13035

A

10.11883／1001－1455（2017）06－0939－07

2015－11－07；

2016－04－22

國(guó)家自然科學(xué)基金項(xiàng)目（11072072）

方前程（1982－ ），男，博士，fangqiancheng314＠126.com。

（責(zé)任編輯 王小飛）