本發(fā)明屬于滑坡災(zāi)害預(yù)報(bào)預(yù)警領(lǐng)域，具體涉及一種基于均值-標(biāo)準(zhǔn)差的加權(quán)馬爾可夫滑坡預(yù)警方法。
背景技術(shù)：
：滑坡作為一種頻發(fā)的自然地質(zhì)災(zāi)害，會(huì)造成巨大的生命財(cái)產(chǎn)損失，其危害已成為僅次于地震的第二大自然災(zāi)害。據(jù)初步統(tǒng)計(jì)，我國(guó)每年由崩塌、滑坡、泥石流等災(zāi)害所造成的直接經(jīng)濟(jì)損失約200億元人民幣，間接損失更是無(wú)法估量。監(jiān)測(cè)預(yù)警作為規(guī)避滑坡等地質(zhì)災(zāi)害的重要手段之一，越來(lái)越受到人們的重視。盡管我國(guó)近10年來(lái)在地質(zhì)災(zāi)害監(jiān)測(cè)預(yù)警方面取得了一定的研究成果，積累了一定的經(jīng)驗(yàn)和教訓(xùn)。但是，在滑坡預(yù)警方面并不理想：一是在成功預(yù)警實(shí)例中專業(yè)預(yù)警比例過(guò)低；二是預(yù)警的滯后性，即預(yù)警不夠及時(shí)；三是錯(cuò)誤預(yù)警較多，即在并無(wú)滑坡危險(xiǎn)時(shí)給出預(yù)警提示，對(duì)生產(chǎn)及生活造成一定干擾；此外，滑坡預(yù)警包括滑坡時(shí)間預(yù)警和滑坡區(qū)域預(yù)警，以往研究中比較側(cè)重滑坡時(shí)間預(yù)警，而對(duì)滑坡區(qū)域預(yù)警的研究不足，而時(shí)間預(yù)警與區(qū)域預(yù)警應(yīng)該同時(shí)實(shí)現(xiàn)。技術(shù)實(shí)現(xiàn)要素：針對(duì)上述現(xiàn)有技術(shù)存在的不足，本發(fā)明提供一種基于均值-標(biāo)準(zhǔn)差的加權(quán)馬爾可夫滑坡預(yù)警方法。本發(fā)明的技術(shù)方案：一種基于均值-標(biāo)準(zhǔn)差的加權(quán)馬爾可夫滑坡預(yù)警方法，包括如下步驟：步驟1：布置滑坡監(jiān)測(cè)點(diǎn)，針對(duì)每一個(gè)滑坡監(jiān)測(cè)點(diǎn)，執(zhí)行步驟2到步驟7；步驟2：每日收集一次滑坡監(jiān)測(cè)點(diǎn)位移速度，構(gòu)成位移速度集合{x1，x2，...，xi，...，xn}，其中，xi為第i日收集到的位移速度，i＝{1,2，...，n}；步驟3：利用均值-標(biāo)準(zhǔn)差法確定位移速度集合中各位移速度所屬狀態(tài)，包括正常狀態(tài)或異常狀態(tài)，n個(gè)位移速度的狀態(tài)組成位移速度狀態(tài)序列；步驟3-1：計(jì)算位移速度的均值和標(biāo)準(zhǔn)差s：x‾=1nΣi=1i=nxis=1n-1Σi=1i=n(xi-x‾)2]]>步驟3-2：確定狀態(tài)劃分的分界點(diǎn)D：D＝+α*s其中，α為可變參數(shù)；步驟3-3：計(jì)算位移速度的最大值max和最小值min；步驟3-4：利用“D”“max”和“min”依次確定位移速度狀態(tài)：若位移速度在區(qū)間[min，D)內(nèi)，則位移速度狀態(tài)為正常狀態(tài)，用“1”表示；若位移速度在區(qū)間[D，max]內(nèi)，則位移速度狀態(tài)為異常狀態(tài)，用“2”表示；n個(gè)位移速度狀態(tài)組成位移速度狀態(tài)序列；步驟4：根據(jù)位移速度狀態(tài)序列，使用加權(quán)馬爾可夫鏈預(yù)測(cè)次日位移速度xn+1所屬狀態(tài)，即位移速度次日預(yù)測(cè)狀態(tài)：步驟4-1：使用馬爾可夫鏈計(jì)算位移速度次日預(yù)測(cè)狀態(tài)的分布矩陣；步驟4-1-1：根據(jù)位移速度狀態(tài)序列，統(tǒng)計(jì)不同步長(zhǎng)的狀態(tài)轉(zhuǎn)移頻數(shù)矩陣q(k)：q(k)=q11q12q21q22]]>其中，k為步長(zhǎng)，k＝{1,2，...，m}，m為階數(shù)，m≤5；q11為位移速度狀態(tài)序列中從正常狀態(tài)間隔k-1個(gè)位移速度狀態(tài)轉(zhuǎn)移到下一個(gè)正常狀態(tài)的轉(zhuǎn)移次數(shù)；q12為位移速度狀態(tài)序列中從正常狀態(tài)間隔k-1個(gè)位移速度狀態(tài)轉(zhuǎn)移到異常狀態(tài)的轉(zhuǎn)移次數(shù)；q21為位移速度狀態(tài)序列中從異常狀態(tài)間隔k-1個(gè)位移速度狀態(tài)轉(zhuǎn)移到正常狀態(tài)的轉(zhuǎn)移次數(shù)；q22為位移速度狀態(tài)序列中從異常狀態(tài)間隔k-1個(gè)位移速度狀態(tài)轉(zhuǎn)移到下一個(gè)異常狀態(tài)的轉(zhuǎn)移次數(shù)；步驟4-1-2：根據(jù)狀態(tài)轉(zhuǎn)移頻數(shù)矩陣計(jì)算狀態(tài)轉(zhuǎn)移概率矩陣P(k)：P(k)=p11p12p21p22]]>其中，步驟4-1-3：計(jì)算初始狀態(tài)概率向量AkT(n)=(p1,p2)]]>其中，步驟4-1-4：計(jì)算m階位移速度次日預(yù)測(cè)狀態(tài)的分布矩陣B：B=B1(n+1)B2(n+1)...Bk(n+1)...Bm(n+1)]]>其中，步驟4-2：使用加權(quán)馬爾可夫鏈預(yù)測(cè)位移速度次日預(yù)測(cè)狀態(tài)：步驟4-2-1：計(jì)算位移速度的自相關(guān)系數(shù)rk：rk=Σt=1t=n-k(xt-x‾)(xt+k-x‾)Σt=1t=n-k(xt-x‾)2·Σt=1t=n-k(xt+k-x‾)2]]>其中，為位移速度均值，xt+k為第t+k日的位移速度，t＝[1，n-k]；步驟4-2-2：計(jì)算不同步長(zhǎng)權(quán)重wk，構(gòu)成m階權(quán)重向量wm＝(w1,w2,…,wk,…,wm)，其中，步驟4-2-3：計(jì)算m階位移速度次日預(yù)測(cè)狀態(tài)的概率分布向量C(n+1)：C(n+1)＝(c1,c2,…,cr)＝wm*B步驟4-2-4：將向量C(n+1)中最大元素所在列對(duì)應(yīng)的狀態(tài)作為m階加權(quán)馬爾可夫鏈預(yù)測(cè)的次日位移速度xn+1所屬狀態(tài)，即位移速度次日預(yù)測(cè)狀態(tài)；步驟4-3：重復(fù)步驟4-1到步驟4-2，計(jì)算m取不同值時(shí)各階位移速度次日預(yù)測(cè)狀態(tài)，如果各階位移速度次日預(yù)測(cè)狀態(tài)均為“異?！?，位移速度次日預(yù)測(cè)狀態(tài)為“異常”；如果各階位移速度次日預(yù)測(cè)狀態(tài)不全為“異?！?，位移速度次日預(yù)測(cè)狀態(tài)為“正常”；步驟5：進(jìn)行滑坡預(yù)警判斷：當(dāng)滿足位移速度前日狀態(tài)為異常狀態(tài)、位移速度當(dāng)日狀態(tài)為異常狀態(tài)且位移速度次日預(yù)測(cè)狀態(tài)為異常狀態(tài)時(shí)，發(fā)出滑坡預(yù)警信號(hào)；步驟6：判斷是否有新收集到的位移速度數(shù)據(jù)，有，則將第一日位移速度從位移速度集合中刪除，將第n+1日收集到的位移速度加入位移速度集合中，對(duì)新的位移速度集合重復(fù)步驟3到步驟5；否則，執(zhí)行步驟7；步驟7：對(duì)滑坡預(yù)警結(jié)果進(jìn)行評(píng)價(jià)：步驟7-1：評(píng)價(jià)滑坡預(yù)警的及時(shí)性：滑坡預(yù)警及時(shí)性以滑坡發(fā)生前滑坡前應(yīng)急準(zhǔn)備時(shí)間內(nèi)滑坡預(yù)警密集度衡量；滑坡預(yù)警密集度越接近1，滑坡預(yù)警的及時(shí)性越強(qiáng)；滑坡預(yù)警密集度為：在滑坡發(fā)生之前R天進(jìn)行的R次滑坡預(yù)警判斷中，發(fā)出滑坡預(yù)警信號(hào)的次數(shù)為R*，則滑坡之前R天內(nèi)的滑坡預(yù)警密集度步驟7-2：評(píng)價(jià)滑坡預(yù)警的抗干擾性：錯(cuò)誤預(yù)警率越低，滑坡預(yù)警的抗干擾性越好；錯(cuò)誤預(yù)警率為：如果在滑坡應(yīng)急準(zhǔn)備時(shí)間之前共做出M次滑坡預(yù)警判斷，其中有M*次發(fā)出了滑坡預(yù)警信號(hào)，錯(cuò)誤預(yù)警率為步驟7-3：評(píng)價(jià)滑坡預(yù)警的可信度：滑坡預(yù)警的可信度由預(yù)測(cè)一致率決定，預(yù)測(cè)一致率越大，滑坡預(yù)警的可信度越好；預(yù)測(cè)一致率為：在滑坡發(fā)生前一共對(duì)位移速度狀態(tài)做Q次預(yù)測(cè)，其中Q*次與實(shí)際位移速度狀態(tài)相符合，則預(yù)測(cè)一致率為有益效果：一種基于均值-標(biāo)準(zhǔn)差的加權(quán)馬爾可夫滑坡預(yù)警方法與現(xiàn)有技術(shù)相比，具有如下優(yōu)勢(shì)：(1)從滑坡演變過(guò)程中獲得的位移參數(shù)的隨機(jī)特性出發(fā)，基于馬爾可夫鏈理論，綜合運(yùn)用滑坡演化三階段理論和均值-標(biāo)準(zhǔn)差分類法，超越了以往工程界對(duì)滑坡預(yù)警的認(rèn)識(shí)；(2)給出了具體預(yù)警過(guò)程，可用于從位移監(jiān)測(cè)數(shù)據(jù)角度出發(fā)的滑坡預(yù)警。附圖說(shuō)明圖1為本發(fā)明一種實(shí)施方式的一種基于均值-標(biāo)準(zhǔn)差的加權(quán)馬爾可夫滑坡預(yù)警方法流程圖；圖2為本發(fā)明一種實(shí)施方式的滑坡區(qū)及監(jiān)測(cè)點(diǎn)布置平面圖；圖3為本發(fā)明一種實(shí)施方式的滑坡前各監(jiān)測(cè)點(diǎn)預(yù)警情況示意圖，(a)為滑坡前第8天及之前各監(jiān)測(cè)點(diǎn)預(yù)警情況示意圖，(b)滑坡前第7天各監(jiān)測(cè)點(diǎn)預(yù)警情況示意圖，(c)滑坡前第6天各監(jiān)測(cè)點(diǎn)預(yù)警情況示意圖，(d)滑坡前第5天各監(jiān)測(cè)點(diǎn)預(yù)警情況示意圖，(e)滑坡前第4天各監(jiān)測(cè)點(diǎn)預(yù)警情況示意圖，(f)滑坡前第3天各監(jiān)測(cè)點(diǎn)預(yù)警情況示意圖，(g)滑坡前第2天各監(jiān)測(cè)點(diǎn)預(yù)警情況示意圖，(h)滑坡前第1天各監(jiān)測(cè)點(diǎn)預(yù)警情況示意圖。具體實(shí)施方式下面結(jié)合附圖對(duì)本發(fā)明的一種實(shí)施方式作詳細(xì)說(shuō)明。本實(shí)施方式中，以露天礦為例，該露天礦于2013年4月17日發(fā)生了侏羅系地層滑坡?；w高度近200m，走向長(zhǎng)700多m，傾向?qū)挾冉?00m，滑體平面投影面積約12萬(wàn)m2，滑面最大深度約50m，體積約300萬(wàn)m3。如圖1所示，基于均值-標(biāo)準(zhǔn)差的加權(quán)馬爾可夫滑坡預(yù)警方法，包括如下步驟：步驟1：布置滑坡監(jiān)測(cè)點(diǎn)，針對(duì)每一個(gè)滑坡監(jiān)測(cè)點(diǎn)，執(zhí)行步驟2到步驟7；如圖2所示，滑坡區(qū)附近分別在▽438、▽462、▽512及▽548共布置了28個(gè)滑坡監(jiān)測(cè)點(diǎn)，其中滑坡區(qū)外布置了7個(gè)滑坡監(jiān)測(cè)點(diǎn)；步驟2：每日收集一次滑坡監(jiān)測(cè)點(diǎn)位移速度，構(gòu)成位移速度集合{x1，x2，…，xi，…，xn}，其中，xi為第i日收集到的位移速度，i＝{1,2，…，n}；本實(shí)施方式中，選擇監(jiān)測(cè)時(shí)間區(qū)段為2013年1月1日至2014年4月17日，選擇樣本容量為n＝20，將前20天獲得的20個(gè)位移速度數(shù)據(jù)組成位移速度集合；鑒于該滑坡過(guò)程中的位移以水平方向?yàn)橹?，本?shí)施方式只對(duì)水平位移數(shù)據(jù)進(jìn)行研究；另一方面，由于數(shù)據(jù)量龐大，僅列出滑坡區(qū)內(nèi)WY2300-512監(jiān)測(cè)點(diǎn)數(shù)據(jù)。步驟3：利用均值-標(biāo)準(zhǔn)差法確定位移速度集合中各位移速度所屬狀態(tài)，包括正常狀態(tài)或異常狀態(tài)，n個(gè)位移速度的狀態(tài)組成位移速度狀態(tài)序列：步驟3-1：計(jì)算位移速度的均值和標(biāo)準(zhǔn)差s：x‾=1nΣi=1i=nxis=1n-1Σi=1i=n(xi-x‾)2]]>本實(shí)施方式中，均值標(biāo)準(zhǔn)差s＝0.6982；步驟3-2：確定狀態(tài)劃分的分界點(diǎn)D：D＝+α*s其中，α為可變參數(shù)；本實(shí)施方式中，α＝0.4，分界點(diǎn)D＝1.3008。步驟3-3：計(jì)算位移速度的最大值max和最小值min；步驟3-4：利用“D”“max”和“min”依次確定位移速度狀態(tài)：若位移速度在區(qū)間[min，D)內(nèi)，則位移速度狀態(tài)為正常狀態(tài)，用“1”表示；若位移速度在區(qū)間[D，max]內(nèi)，則位移速度狀態(tài)為異常狀態(tài)，用“2”表示；n個(gè)位移速度狀態(tài)組成位移速度狀態(tài)序列；本實(shí)施方式中，位移速度的最小值min為-0.12，最大值max為2.39，若位移速度在區(qū)間[-0.12，1.30)內(nèi)，則位移速度狀態(tài)為異常狀態(tài)，用“1”表示；若位移速度在區(qū)間[1.30，2.39]內(nèi)，則位移速度狀態(tài)為異常狀態(tài)，用“2”表示。位移速度狀態(tài)結(jié)果見表1：表1WY2300-512監(jiān)測(cè)點(diǎn)位移速度狀態(tài)表步驟4：根據(jù)位移速度狀態(tài)序列，使用加權(quán)馬爾可夫鏈預(yù)測(cè)次日位移速度xn+1所屬狀態(tài)，即位移速度次日預(yù)測(cè)狀態(tài)：步驟4-1：使用馬爾可夫鏈計(jì)算m階位移速度次日預(yù)測(cè)狀態(tài)的分布矩陣；步驟4-1-1：根據(jù)位移速度狀態(tài)序列，統(tǒng)計(jì)不同步長(zhǎng)的狀態(tài)轉(zhuǎn)移頻數(shù)矩陣q(k)：q(k)=q11q12q21q22]]>其中，k為步長(zhǎng)，k＝{1,2，…，m}，m為階數(shù)，m≤5；q11為位移速度狀態(tài)序列中從正常狀態(tài)間隔k-1個(gè)位移速度狀態(tài)轉(zhuǎn)移到下一個(gè)正常狀態(tài)的轉(zhuǎn)移次數(shù)；q12為位移速度狀態(tài)序列中從正常狀態(tài)間隔k-1個(gè)位移速度狀態(tài)轉(zhuǎn)移到異常狀態(tài)的轉(zhuǎn)移次數(shù)；q21為位移速度狀態(tài)序列中從異常狀態(tài)間隔k-1個(gè)位移速度狀態(tài)轉(zhuǎn)移到正常狀態(tài)的轉(zhuǎn)移次數(shù)；q22為位移速度狀態(tài)序列中從異常狀態(tài)間隔k-1個(gè)位移速度狀態(tài)轉(zhuǎn)移到下一個(gè)異常狀態(tài)的轉(zhuǎn)移次數(shù)；步驟4-1-2：根據(jù)狀態(tài)轉(zhuǎn)移頻數(shù)矩陣計(jì)算狀態(tài)轉(zhuǎn)移概率矩陣P(k)：P(k)=p11p12p21p22]]>其中，本實(shí)施方式中，m＝5，步長(zhǎng)k分別為1，2，3，4，5的狀態(tài)轉(zhuǎn)移頻數(shù)矩陣和狀態(tài)轉(zhuǎn)移概率矩陣如下：q1=7651P1=0.53850.46150.83330.1667q2=7542P2=0.58330.41670.66670.3333]]>q3=6551P3=0.54550.45450.83330.1667q4=8233P4=0.80000.20000.50000.5000]]>q5=5550P5=0.50000.500010]]>步驟4-1-3：計(jì)算初始狀態(tài)概率向量AkT(n)=(p1,p2)]]>其中，本實(shí)施方式中，第20個(gè)位移速度狀態(tài)為2，步長(zhǎng)為1的初始狀態(tài)概率向量同理，第19，18，17，16個(gè)位移速度狀態(tài)為分別為1，1，1，2，步長(zhǎng)2、3、4、5對(duì)應(yīng)的初始狀態(tài)概率向量分別為：步驟4-1-4：計(jì)算m階位移速度次日預(yù)測(cè)狀態(tài)的分布矩陣B：B=B1(n+1)B2(n+1)...Bk(n+1)...Bm(n+1)]]>其中，步驟4-2：使用加權(quán)馬爾可夫鏈預(yù)測(cè)m階位移速度次日預(yù)測(cè)狀態(tài)：步驟4-2-1：計(jì)算位移速度的自相關(guān)系數(shù)rk：rk=Σt=1t=n-k(xt-x‾)(xt+k-x‾)Σt=1t=n-k(xt-x‾)2·Σt=1t=n-k(xt+k-x‾)2]]>其中，為位移速度均值，xt+k為第t+k日的位移速度，t＝[1，n-k]；本實(shí)施方式中，各步長(zhǎng)自相關(guān)系數(shù)如表2：表2各步長(zhǎng)自相關(guān)系數(shù)表步驟4-2-2：計(jì)算不同步長(zhǎng)權(quán)重wk，構(gòu)成m階權(quán)重向量wm＝(w1,w2,…,wk,…,wm)，其中，本實(shí)施方式中，各步長(zhǎng)的權(quán)重向量如表3：表3各步長(zhǎng)權(quán)重向量表步驟4-2-3：計(jì)算m階位移速度次日預(yù)測(cè)狀態(tài)的概率分布向量C(n+1)：C(n+1)＝(c1,c2,…,cr)＝wm*B步驟4-2-4：將向量C(n+1)中最大元素所在列對(duì)應(yīng)的狀態(tài)作為m階加權(quán)馬爾可夫鏈預(yù)測(cè)的次日位移速度xn+1預(yù)測(cè)狀態(tài)；步驟4-3：重復(fù)步驟4-1到步驟4-2，計(jì)算m取不同值時(shí)各階位移速度次日預(yù)測(cè)狀態(tài)，如果各階位移速度次日預(yù)測(cè)狀態(tài)均為“異?！?，位移速度次日預(yù)測(cè)狀態(tài)為“異?！?；如果各階位移速度次日預(yù)測(cè)狀態(tài)不全為“異?！?，位移速度次日預(yù)測(cè)狀態(tài)為“正?！?；本實(shí)施方式中，各階位移速度次日預(yù)測(cè)狀態(tài)的概率分布如表4：表4位移速度數(shù)據(jù)狀態(tài)概率分布表步驟5：進(jìn)行滑坡預(yù)警判斷：當(dāng)滿足位移速度前日狀態(tài)為異常狀態(tài)、位移速度當(dāng)日狀態(tài)為異常狀態(tài)且位移速度次日預(yù)測(cè)狀態(tài)為異常狀態(tài)時(shí)，發(fā)出滑坡預(yù)警信號(hào)；未來(lái)位移速度的預(yù)測(cè)狀態(tài)信息以及當(dāng)前位移速度的真實(shí)狀態(tài)信息，概括了馬爾可夫鏈預(yù)測(cè)理論對(duì)邊坡穩(wěn)定性狀況的評(píng)價(jià)，基于這些信息即可對(duì)是否給出預(yù)警做出判斷。為確保預(yù)測(cè)的可信度，僅在各階加權(quán)馬爾可夫鏈的位移速度次日預(yù)測(cè)狀態(tài)均為異常狀態(tài)時(shí)，綜合預(yù)測(cè)結(jié)果才確定為異常狀態(tài)。由于監(jiān)測(cè)數(shù)據(jù)不斷更新，應(yīng)及時(shí)納入新數(shù)據(jù)，剔除歷史數(shù)據(jù)，并不斷以加權(quán)馬爾可夫鏈預(yù)測(cè)未來(lái)位移速度狀態(tài)。在給出每個(gè)次日位移速度預(yù)測(cè)狀態(tài)的同時(shí)，也要做出預(yù)警與否的決定。當(dāng)連續(xù)出現(xiàn)真實(shí)狀態(tài)和預(yù)測(cè)狀態(tài)均為異常狀態(tài)時(shí)，應(yīng)發(fā)出預(yù)警信號(hào)。這個(gè)連續(xù)出現(xiàn)異常狀態(tài)的階段對(duì)應(yīng)于滑坡演化三階段理論中位移加速度大于零的階段，根據(jù)加速度a≥0這一判據(jù)的滑坡預(yù)警判斷方法為：步驟6：判斷是否有新收集到的位移速度數(shù)據(jù)，有，則將第一日位移速度從位移速度集合中刪除，將第n+1日收集到的實(shí)際位移速度加入位移速度集合中，對(duì)新的位移速度集合重復(fù)步驟3到步驟5；否則，執(zhí)行步驟7；本實(shí)施方式中，監(jiān)測(cè)點(diǎn)位移速度預(yù)測(cè)狀態(tài)結(jié)果及預(yù)警情況如表5：表5WY2300-512監(jiān)測(cè)點(diǎn)位移速度預(yù)測(cè)狀態(tài)結(jié)果及預(yù)警情況表步驟7：對(duì)滑坡預(yù)警結(jié)果進(jìn)行評(píng)價(jià)：步驟7-1：計(jì)算滑坡預(yù)警的及時(shí)性：滑坡預(yù)警及時(shí)性以滑坡發(fā)生前滑坡前應(yīng)急準(zhǔn)備時(shí)間內(nèi)滑坡預(yù)警密集度衡量；本實(shí)施方式中，滑坡前應(yīng)急準(zhǔn)備時(shí)間為7天，故滑坡預(yù)警及時(shí)性以滑坡發(fā)生前7天內(nèi)滑坡預(yù)警密集度，即7日滑坡預(yù)警密集度來(lái)衡量；對(duì)于發(fā)生滑坡的監(jiān)測(cè)點(diǎn)，滑坡預(yù)警密集度應(yīng)是一個(gè)動(dòng)態(tài)變化逐漸接近1的變量；滑坡預(yù)警密集度越接近1，滑坡預(yù)警的及時(shí)性越強(qiáng)；而對(duì)于未發(fā)生滑坡的監(jiān)測(cè)點(diǎn)，滑坡預(yù)警密集度則應(yīng)是一個(gè)接近0的變量?；骂A(yù)警密集度為：在滑坡發(fā)生之前R天進(jìn)行的R次滑坡預(yù)警判斷中，發(fā)出滑坡預(yù)警信號(hào)的次數(shù)為R*，則滑坡之前R天內(nèi)的滑坡預(yù)警密集度步驟7-2：評(píng)價(jià)滑坡預(yù)警的抗干擾性：根據(jù)實(shí)際需要，錯(cuò)誤預(yù)警率越低，滑坡預(yù)警的抗干擾性越好；錯(cuò)誤預(yù)警率為：如果在滑坡應(yīng)急準(zhǔn)備時(shí)間之前共做出M次滑坡預(yù)警判斷，其中有M*次發(fā)出了滑坡預(yù)警信號(hào)，錯(cuò)誤預(yù)警率為步驟7-3：評(píng)價(jià)滑坡預(yù)警的可信度：滑坡預(yù)警的可信度由預(yù)測(cè)一致率決定，預(yù)測(cè)一致率越大，滑坡預(yù)警的可信度越好；預(yù)測(cè)一致率為：在滑坡發(fā)生前一共對(duì)位移速度狀態(tài)做Q次預(yù)測(cè)，其中Q*次與實(shí)際位移速度狀態(tài)相符合，則預(yù)測(cè)一致率為本實(shí)施方式中，由表5可知，對(duì)監(jiān)測(cè)點(diǎn)的7日預(yù)警密集度為57％，錯(cuò)誤預(yù)警率為0，預(yù)測(cè)一致率為58％。實(shí)際情況是2013/4/17發(fā)生滑坡，在83次預(yù)測(cè)中僅在滑坡前的第4、3、2、1日連續(xù)給出四次預(yù)警，滿足預(yù)警及時(shí)性標(biāo)準(zhǔn)，而錯(cuò)誤預(yù)警率為零則滿足了預(yù)警抗干擾性標(biāo)準(zhǔn)，預(yù)測(cè)一致率為58％，具有一定的可信度。對(duì)滑坡區(qū)域內(nèi)全部監(jiān)測(cè)點(diǎn)的預(yù)警結(jié)果評(píng)價(jià)見表6：表6基于均值-標(biāo)準(zhǔn)差的加權(quán)馬爾可夫滑坡預(yù)警方法預(yù)警評(píng)價(jià)表滑坡區(qū)域內(nèi)外各監(jiān)測(cè)點(diǎn)臨滑前15天的預(yù)警判斷結(jié)果，見表7。表中“--”表示當(dāng)日未給出預(yù)警信號(hào)，灰色框表示當(dāng)日給出預(yù)警信號(hào)。表7監(jiān)測(cè)點(diǎn)在滑坡前15天內(nèi)的預(yù)警結(jié)果表監(jiān)測(cè)區(qū)域內(nèi)全部監(jiān)測(cè)點(diǎn)在滑坡前的預(yù)警情況見圖3(a)-圖3(h)，從圖中監(jiān)測(cè)點(diǎn)的動(dòng)態(tài)預(yù)警情況可以看出，本方法具有較好的時(shí)間預(yù)警和區(qū)域預(yù)警效果，能在一定程度上揭示滑坡的啟動(dòng)及演變過(guò)程。當(dāng)前第1頁(yè)1 2 3