張 可，崔敬浩，黃華愛，豐景春，左 媛

(1.河海大學(xué)商學(xué)院，江蘇 南京 211100； 2.國(guó)際河流研究中心，江蘇 南京 211100；3.江蘇省“世界水谷”與水生態(tài)文明協(xié)同創(chuàng)新中心，江蘇 南京 211100；4.河海大學(xué)項(xiàng)目管理研究所，江蘇 南京 211100；5.廣西壯族自治區(qū)防汛抗旱指揮部辦公室，廣西 南寧 530023)

近年來，我國(guó)由于突發(fā)性強(qiáng)降雨而致使泥石流、山洪災(zāi)害與日俱增，其中，由于山洪災(zāi)害而死亡的人數(shù)已經(jīng)超過所有洪澇災(zāi)害死亡率的50%[1-2]。發(fā)生洪災(zāi)幾率較高的小型水庫絕大部分均誕生于20世紀(jì)50年代，由于當(dāng)時(shí)技術(shù)、經(jīng)濟(jì)較弱，因而這些小型水庫存在在建設(shè)過程中建設(shè)標(biāo)準(zhǔn)較低、管理不全面，以及現(xiàn)在的年久失修等各種情況。這些小型水庫無法保證良好的安全性，大小事故頻發(fā)，不利于水庫流域下游安全。據(jù)1949—2009年史料記載，全國(guó)各類水庫垮壩失事3 487座，其中小型水庫為 3 361座，占垮壩失事總數(shù)的96.4%[3]。

國(guó)家高度重視災(zāi)后水利薄弱環(huán)節(jié)建設(shè)工作。水利部于2017年通報(bào)了《加快災(zāi)后水利薄弱環(huán)節(jié)實(shí)施方案》，其中，小型水庫除險(xiǎn)加固是4個(gè)重點(diǎn)內(nèi)容之一。小型水庫在灌溉、發(fā)電等諸多方面均產(chǎn)生了極其可觀的生產(chǎn)效益，然而，由于大部分小型水庫未聘請(qǐng)專業(yè)的機(jī)構(gòu)進(jìn)行監(jiān)管，國(guó)家也未下?lián)芫S養(yǎng)經(jīng)費(fèi)，加上年久失修等各種問題，這些水庫往往達(dá)不到較高的防洪標(biāo)準(zhǔn)。雖然當(dāng)前國(guó)家已經(jīng)逐步圍繞小型水庫進(jìn)行了一定的除險(xiǎn)加固工作，使得某些硬件設(shè)施日益完善，安全性較以往有了極大的提高，然而，依然存在大部分小型水庫由于缺乏資金而無法基于國(guó)家標(biāo)準(zhǔn)進(jìn)行加固的問題。

我國(guó)小型水庫眾多，如截至2013年，廣西壯族自治區(qū)水利系統(tǒng)注冊(cè)登記的水庫共4 362座，其中小型水庫為4 165座[4]?，F(xiàn)有小型水庫多建于20世紀(jì)，普遍存在設(shè)計(jì)標(biāo)準(zhǔn)不統(tǒng)一、建設(shè)資料不全、安全責(zé)任不落實(shí)、管理粗放等問題，且缺乏相應(yīng)的檢測(cè)儀器和設(shè)備，難以獲取反映大壩安全狀態(tài)的指標(biāo)，給大壩安全風(fēng)險(xiǎn)判斷和預(yù)警帶來困難，從而造成了水庫安全性評(píng)估的困難。如何對(duì)全國(guó)93 308座小型水庫的安全性進(jìn)行分級(jí)量化，精準(zhǔn)地開展高風(fēng)險(xiǎn)水庫的除險(xiǎn)加固工作是目前迫切需要研究的問題。

1 研究綜述

1.1 國(guó)外研究現(xiàn)狀

美國(guó)是最早將風(fēng)險(xiǎn)理論與風(fēng)險(xiǎn)評(píng)價(jià)應(yīng)用于大壩安全管理評(píng)價(jià)的國(guó)家。20世紀(jì)70年代，美國(guó)在圍繞核電站而展開的風(fēng)險(xiǎn)評(píng)估工作中，首次引入相關(guān)的風(fēng)險(xiǎn)方法，在此之后，風(fēng)險(xiǎn)評(píng)價(jià)日益受到各國(guó)學(xué)者的關(guān)注。學(xué)者卡薩格蘭特提出，需要基于風(fēng)險(xiǎn)方面，對(duì)大壩進(jìn)行科學(xué)合理的安全評(píng)估。這一提議受到美國(guó)學(xué)者的認(rèn)可，20世紀(jì)80年代，美國(guó)正式在大壩安全管理方面引入完善的風(fēng)險(xiǎn)評(píng)估方法[5]。

1973年，ASCE圍繞溢洪道具備的泄洪能力進(jìn)行科學(xué)有效的評(píng)估，并在出具評(píng)估報(bào)告時(shí)，引入某特定的應(yīng)用風(fēng)險(xiǎn)分析方法[6]。1982年，Hagen首次提及風(fēng)險(xiǎn)的基本定義，并借助相對(duì)風(fēng)險(xiǎn)指數(shù)對(duì)某大壩的風(fēng)險(xiǎn)進(jìn)行科學(xué)評(píng)價(jià)[7]。此外，美國(guó)相關(guān)部門在應(yīng)用風(fēng)險(xiǎn)分析等方法綜合管理350個(gè)大壩的過程中，基于現(xiàn)場(chǎng)評(píng)分法，對(duì)大壩本身的安全風(fēng)險(xiǎn)進(jìn)行科學(xué)衡量[8]。同時(shí)，美國(guó)氣象局為深入探究大壩潰壩可能會(huì)對(duì)下游產(chǎn)生的實(shí)際影響，認(rèn)為有必要細(xì)致剖析潰壩洪水，并制定科學(xué)豐富的潰壩模型，自最初的DAMBRK至后期的FLDWAV，均有助于進(jìn)行精準(zhǔn)的潰壩洪水計(jì)算[9]。1991年，BC Hydro首次在針對(duì)大壩安全風(fēng)險(xiǎn)進(jìn)行科學(xué)評(píng)估的過程中引入應(yīng)用概率的評(píng)估方法，其中風(fēng)險(xiǎn)評(píng)估程序不僅涉及風(fēng)險(xiǎn)分析，而且關(guān)聯(lián)風(fēng)險(xiǎn)評(píng)估[10]。

在大壩管理過程中，澳大利亞為首批引入風(fēng)險(xiǎn)分析的國(guó)家，在風(fēng)險(xiǎn)分析方面成果頗豐[11-14]，尤其是風(fēng)險(xiǎn)排序技術(shù)等。澳大利亞1994年出臺(tái)了《澳大利亞風(fēng)險(xiǎn)評(píng)價(jià)指南》，1998年出臺(tái)了《大壩環(huán)境管理指南》，2000年出臺(tái)了《大壩潰決后果評(píng)價(jià)指南》，正是由于這些法規(guī)的相繼出臺(tái)，才使得現(xiàn)有的大壩管理工作極為完善。

2013年，在英國(guó)出臺(tái)的大壩管理指南中，將土石壩需要進(jìn)行的風(fēng)險(xiǎn)評(píng)估工作細(xì)致劃分為：①層次1：可應(yīng)用基本的若干定性方法和各項(xiàng)技術(shù)，同時(shí)可細(xì)致描述風(fēng)險(xiǎn)；②層次2：可應(yīng)用基本的若干定量方法，并可判斷出險(xiǎn)概率；③層次3：可應(yīng)用復(fù)雜細(xì)致的相關(guān)定量方法，并可進(jìn)行不確定分析[15]。在該指南中，以若干估算方法為基礎(chǔ)，并綜合考量風(fēng)險(xiǎn)后果，根據(jù)民眾對(duì)出險(xiǎn)模式的領(lǐng)會(huì)程度逐步深入。曾有學(xué)者對(duì)該指南的實(shí)踐應(yīng)用效果進(jìn)行評(píng)價(jià)，并歸納英國(guó)HSE公布的若干研究成果，確認(rèn)該指南能基于個(gè)體及社會(huì)方面有效評(píng)估大壩風(fēng)險(xiǎn)。

除美國(guó)和澳大利亞之外，瑞典、葡萄牙以及日本等國(guó)家[16-20]紛紛在本國(guó)的水利工程中引入風(fēng)險(xiǎn)評(píng)價(jià)體系，正因如此，風(fēng)險(xiǎn)分析才得以在當(dāng)前的水利領(lǐng)域中迅猛發(fā)展。

1.2 國(guó)內(nèi)研究現(xiàn)狀

20世紀(jì)50年代，中國(guó)已經(jīng)圍繞大壩的安全性能進(jìn)行各項(xiàng)研究，最初主要面向三峽工程，然而，由于技術(shù)等條件的局限，當(dāng)時(shí)尚未清晰大壩的潰壩機(jī)理。受國(guó)外水庫大壩風(fēng)險(xiǎn)管理理念的影響，20世紀(jì)80年代，中國(guó)首次引入風(fēng)險(xiǎn)分析對(duì)大壩進(jìn)行深入探究。當(dāng)前中國(guó)逐步提高對(duì)大壩風(fēng)險(xiǎn)問題的研究力度，諸多學(xué)者圍繞大壩風(fēng)險(xiǎn)評(píng)價(jià)進(jìn)行各項(xiàng)研究，并在理論與方法上均獲得突出的成效。徐祖信等[21]應(yīng)用JC法并引入可靠度公式對(duì)開敞式溢洪道所具備的泄洪風(fēng)險(xiǎn)進(jìn)行深入剖析；姜樹海[22-23]引入可靠性分析對(duì)面流消能進(jìn)行相應(yīng)的工優(yōu)化設(shè)計(jì)；李君純等[24]提出了評(píng)估大壩安全度的合理方法，并擬定了科學(xué)的判別標(biāo)準(zhǔn)，基于可靠度將判別元素設(shè)定為源自經(jīng)濟(jì)等條件的影響，由此獲得工程的匹配級(jí)別；尉維斌等[25]提出可借助于模糊數(shù)學(xué)法，獲得大壩實(shí)測(cè)性態(tài)評(píng)價(jià)體系所涉及的多樣化元素權(quán)重；謝崇寶等[26]針對(duì)水庫自身的防洪風(fēng)險(xiǎn)進(jìn)行科學(xué)細(xì)致的分析，并深入探究不確定風(fēng)險(xiǎn)的實(shí)際分布；李雷等[27]基于可靠度的定義，針對(duì)某特定土石壩含有的安全度進(jìn)行科學(xué)合理的綜合評(píng)價(jià)。1999年，李君純等[28]基于西方發(fā)達(dá)國(guó)家已擬定的先進(jìn)體系，提出安全度的具體判別標(biāo)準(zhǔn)，根據(jù)該標(biāo)準(zhǔn)，基于大壩表現(xiàn)的總體安全度SD，即可獲得其風(fēng)險(xiǎn)程度。該方法實(shí)則已將安全度量化，基于可靠度并綜合考量社會(huì)影響等因素，極具系統(tǒng)化，也使得中國(guó)自工程安全管理逐步轉(zhuǎn)向風(fēng)險(xiǎn)管理進(jìn)行探究。

2000年，國(guó)家出臺(tái)的SL 258—2000《水庫大壩安全導(dǎo)則》，將大壩風(fēng)險(xiǎn)評(píng)價(jià)細(xì)分為若干子系統(tǒng)，其中包含但并不僅限于工程質(zhì)量評(píng)價(jià)、防洪安全評(píng)價(jià)、結(jié)構(gòu)安全評(píng)價(jià)以及相應(yīng)的金屬結(jié)構(gòu)安全評(píng)價(jià)等，還針對(duì)這些子系統(tǒng)依次擬定科學(xué)合理的評(píng)價(jià)標(biāo)準(zhǔn)，并將大壩安全細(xì)分為若干等級(jí)。何金平等[29]借助于模糊積分評(píng)價(jià)模型的作用，針對(duì)某特定大壩所表現(xiàn)的結(jié)構(gòu)實(shí)測(cè)性態(tài)進(jìn)行科學(xué)評(píng)價(jià)；姜樹海等[30]基于數(shù)值模擬等技術(shù)，自洪水災(zāi)害風(fēng)險(xiǎn)方面入手，針對(duì)洪水調(diào)度以及滯洪等環(huán)節(jié)進(jìn)行精準(zhǔn)的數(shù)值模擬，并對(duì)可能出現(xiàn)的防洪失控風(fēng)險(xiǎn)進(jìn)行分析評(píng)估；周元春等[31]論述了大壩所關(guān)聯(lián)的安全風(fēng)險(xiǎn)評(píng)估工作，需要秉持的基本原理和操作環(huán)節(jié)，并對(duì)可能影響大壩風(fēng)險(xiǎn)的一系列因素進(jìn)行探討，闡述事故樹分析法等多樣化方法中存在的問題；李清富等[32]針對(duì)大壩在風(fēng)險(xiǎn)分析方面涉及的各項(xiàng)理論及諸多方法進(jìn)行歸納，包含風(fēng)險(xiǎn)估計(jì)等方法；楊禮祥等[33]借助于多因素評(píng)判指標(biāo)法，對(duì)水庫所關(guān)聯(lián)的防洪調(diào)度風(fēng)險(xiǎn)展開科學(xué)的定量分析，并對(duì)風(fēng)險(xiǎn)主體彼此間的數(shù)量特征關(guān)系進(jìn)行深入探究，逐一獲得這些主體的風(fēng)險(xiǎn)率，以此來擬定科學(xué)合理的評(píng)判方式；劉亞蓮等[34]借助于模糊層次綜合評(píng)價(jià)法對(duì)某特定土石壩含有的安全性進(jìn)行科學(xué)合理的綜合評(píng)價(jià)；喻蔚然[35]將灰色系統(tǒng)設(shè)定為不確定性，并將其與隨機(jī)概率緊密銜接，以此來構(gòu)建相應(yīng)的功能函數(shù)，并通過持續(xù)優(yōu)化，僅計(jì)算一般隨機(jī)風(fēng)險(xiǎn)概率，簡(jiǎn)潔性極佳。

1.3 現(xiàn)有研究評(píng)述

在已建大壩的各國(guó)中，中國(guó)的大壩數(shù)量最多，面臨的問題也較多。由于受到經(jīng)濟(jì)等因素的局限，很難在短期內(nèi)對(duì)全部存在問題的大壩進(jìn)行全面加固，故而必須對(duì)潰壩型以及危險(xiǎn)型等諸多問題大壩所面臨的運(yùn)行風(fēng)險(xiǎn)進(jìn)行科學(xué)合理的分析，對(duì)更具風(fēng)險(xiǎn)度的少數(shù)大壩進(jìn)行退役評(píng)估，或?qū)踩L(fēng)險(xiǎn)最大的少數(shù)大壩作為除險(xiǎn)加固對(duì)象，才能有效解決當(dāng)前大壩面臨的一系列安全問題。

我國(guó)目前剛剛推進(jìn)風(fēng)險(xiǎn)管理的步伐，雖然已引入全球前沿的風(fēng)險(xiǎn)分析技術(shù)，且從事水利行業(yè)的工作者也逐步意識(shí)到大壩風(fēng)險(xiǎn)的重要性，然而，受監(jiān)測(cè)數(shù)據(jù)缺乏和相應(yīng)方法的限制，國(guó)內(nèi)對(duì)小型水庫大壩安全風(fēng)險(xiǎn)的評(píng)估主觀性較強(qiáng)。雖然在大壩的安全評(píng)估方面已研究出多樣化方法，但不夠系統(tǒng)化和統(tǒng)一化，已擬定的安全分類標(biāo)準(zhǔn)也未設(shè)定合理的定量界限，故而各學(xué)者的評(píng)判結(jié)果也各不相同，可能發(fā)生兩種錯(cuò)誤情況：其一，將安全水庫錯(cuò)誤劃分為病險(xiǎn)水庫，阻止其蓄水并投入使用，這樣不僅減少了水庫效益，還要耗費(fèi)額外的資金對(duì)其加固；其二，將病險(xiǎn)水庫錯(cuò)誤劃分為安全水庫，造成其蓄水直至潰壩，使得下游蒙受巨額的經(jīng)濟(jì)損失。為此，本文針對(duì)目前小型水庫安全風(fēng)險(xiǎn)評(píng)價(jià)存在的不足，研究并構(gòu)建基于支持向量的小型水庫安全風(fēng)險(xiǎn)等級(jí)評(píng)價(jià)模型。

2 研究方法

2.1 支持向量機(jī)原理

支持向量機(jī)涉及的一系列學(xué)習(xí)算法，屬于當(dāng)前典型的監(jiān)督學(xué)習(xí)模型之一，能分析各項(xiàng)數(shù)據(jù)，還可辨別諸多模式，適用于回歸分析，其核心思想如下：建立某特定的最優(yōu)決策超平面，將此平面兩側(cè)和其自身最近的兩個(gè)類進(jìn)行最大化，以泛化各分類問題。對(duì)表現(xiàn)為多維的某特定樣本集，系統(tǒng)將會(huì)隨機(jī)形成相應(yīng)超平面并持續(xù)移動(dòng)，從而實(shí)現(xiàn)對(duì)樣本的分類，一直持續(xù)到某些其余種類的樣本點(diǎn)也處在此平面兩側(cè)為止。實(shí)際上，可滿足上述條件的超平面或許存在若干個(gè)，SVM基于有效保障分類精準(zhǔn)性的前提下，需要找到滿足以上條件的超平面，并使其兩側(cè)最大化，以達(dá)到最優(yōu)分類的效果。

假定線性可分樣本集為：T={(X1，Y1)，…，(Xn，Yn)}，Xi∈Rm，Yi∈{-1，1}，i=1，2，…，n，其中Rm為m維向量空間。分類面方程為(wXi)+b=0，其中w為權(quán)重，b為函數(shù)的偏置量。分類判別如下：

Yi[(wXi)+b]-1≥0

(1)

滿足式(1)的向量被稱之為支持向量。當(dāng)處于樣本線性可分類的情況下，若想求解以最優(yōu)超平面為基礎(chǔ)而構(gòu)建的決策函數(shù)，僅需將其看作二次規(guī)劃問題即可?；诮馕鰩缀螌⒛艿弥?，類間間隔表現(xiàn)為D=2/‖ω‖，此時(shí)引入Lagrange函數(shù)，即可完成求解：

(2)

式中ai>0主要代表Lagrange乘子，基于Kuhn-Tucker條件，此求解過程需要滿足：

ai{[(wXi)+b]Yi-1}≥0

(3)

通常情況下，絕大部分樣本ai均等于0，若其不等于0，則匹配的樣本即為支持向量。經(jīng)由一系列求解，可獲得判別函數(shù)如下：

(4)

由于超平面無法對(duì)某些樣本進(jìn)行精準(zhǔn)分類，故而需引入懲罰參數(shù)C和松弛變量ξi等對(duì)現(xiàn)有函數(shù)進(jìn)行修正，使得模型保持極佳的容錯(cuò)性：

(5)

若支持向量機(jī)遇到的樣本為線性不可分樣本，則需要基于核函數(shù)K對(duì)其進(jìn)行相應(yīng)的非線性變換，將其轉(zhuǎn)化為線性問題進(jìn)行求解，最后再映射為非線性。此外，核函數(shù)有助于支持向量機(jī)在空間中進(jìn)行一定的積轉(zhuǎn)換運(yùn)算，不僅可有效解決線性不可分問題，而且還能切實(shí)解決維數(shù)災(zāi)難問題，簡(jiǎn)化計(jì)算量。一般情況下，核函數(shù)主要表現(xiàn)為多項(xiàng)式函數(shù)等。

2.2 研究思路

小型水庫存在建設(shè)年代久遠(yuǎn)、建設(shè)資料不全、監(jiān)測(cè)設(shè)備缺少等問題，造成了小型水庫普遍存在相關(guān)數(shù)據(jù)資料零散不全的現(xiàn)狀，而其他的機(jī)器學(xué)習(xí)方法往往對(duì)數(shù)據(jù)樣本的數(shù)量要求較高，少量零散多指標(biāo)的樣本數(shù)據(jù)無法建立精確度高且實(shí)用的模型。經(jīng)過對(duì)一系列算法進(jìn)行選擇和評(píng)估，選用機(jī)器學(xué)習(xí)中的支持向量機(jī)作為建模的基礎(chǔ)算法。

在小型水庫安全風(fēng)險(xiǎn)評(píng)估中，使用支持向量機(jī)對(duì)數(shù)據(jù)進(jìn)行分類使用。先將分類屬性設(shè)定為安全風(fēng)險(xiǎn)評(píng)估指標(biāo)體系，再將其安全風(fēng)險(xiǎn)等級(jí)作為分類類別進(jìn)行建模，通過模型參數(shù)的訓(xùn)練和測(cè)試，得到評(píng)估精確度最佳的模型。

支持向量機(jī)的模型設(shè)計(jì)主要包括兩個(gè)部分：①指標(biāo)體系的確立；②模型的訓(xùn)練和測(cè)試，通過訓(xùn)練確定參數(shù)。模型構(gòu)建完成，即可對(duì)小型水庫進(jìn)行風(fēng)險(xiǎn)評(píng)估。

3 模型構(gòu)建及結(jié)果

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

數(shù)據(jù)來源于廣西壯族自治區(qū)山洪災(zāi)害數(shù)據(jù)庫，數(shù)據(jù)庫表包括水庫基本信息表、水庫水文特征表、大壩表、下游影響表、水庫管理體制表、水庫運(yùn)行管理表。數(shù)據(jù)庫中各表根據(jù)水利部批準(zhǔn)的SL 700—2015《水利工程建設(shè)與管理數(shù)據(jù)庫表結(jié)構(gòu)及標(biāo)識(shí)符》建立，均以水庫代碼，即水庫工程的識(shí)別代碼作為主鍵。

根據(jù)之前的研究成果，在原有指標(biāo)體系的基礎(chǔ)上選取小型水庫安全風(fēng)險(xiǎn)指標(biāo)，剔除了地震基本烈度、設(shè)防地震烈度、注冊(cè)登記與否、壩型、壩基防滲措施和庫區(qū)重大滑坡體等6個(gè)不區(qū)分類別的指標(biāo)，最后共選取指標(biāo)14個(gè)。選取的指標(biāo)包括歷年最大泄洪量、多年平均最低氣溫、地震基本烈度、已淤積庫容、壩前淤積面高程等包含小型水庫所處的極端事件、工程地質(zhì)、氣象水文、自然地理，以及包含技術(shù)、組織管理、人員等多方面因素。選取的字段見表1。

表1 數(shù)據(jù)字段及說明

選定字段后，對(duì)各表進(jìn)行表連接操作，共得到946條數(shù)據(jù)，再對(duì)得到的數(shù)據(jù)進(jìn)行預(yù)處理。具體步驟如下。

a.數(shù)據(jù)清洗。水庫數(shù)據(jù)庫中的水庫類型包括大(1)型水庫、大(2)型水庫、中型水庫、小(1)型水庫和小(2)型水庫。先篩選出小(1)型水庫和小(2)型水庫數(shù)據(jù)。

b.抽取非空數(shù)據(jù)。根據(jù)各個(gè)指標(biāo)字段，提取非空數(shù)據(jù)。由于小型水庫數(shù)據(jù)較為匱乏，在完成數(shù)據(jù)提取和清洗操作后，共得到23個(gè)較為完整的樣本數(shù)據(jù)，包括16個(gè)小(1)型水庫，7個(gè)小(2)型水庫。對(duì)于部分空字段，采取數(shù)據(jù)修補(bǔ)方法對(duì)數(shù)據(jù)進(jìn)行補(bǔ)齊。若該字段指標(biāo)為分類數(shù)據(jù)，對(duì)空數(shù)據(jù)用其他類別表示；若為數(shù)值型數(shù)據(jù)，采用該指標(biāo)字段的平均值進(jìn)行補(bǔ)齊。在得到的小型水庫數(shù)據(jù)中，80%以上建于20世紀(jì)50—70年代，壩型均為均質(zhì)土壩，地震基本烈度均為5級(jí)，設(shè)防地震烈度為0，全都登記在冊(cè)，壩基防滲措施，且無庫區(qū)重大滑坡體事件發(fā)生，剔除這6個(gè)指標(biāo)，最終的數(shù)據(jù)總體情況見圖1。

圖1 數(shù)據(jù)展示

3.2 結(jié)果分析

3.2.1基于災(zāi)害通告的水庫等級(jí)劃分

依據(jù)廣西壯族自治區(qū)各水庫洪水災(zāi)害通告劃分水庫等級(jí)。等級(jí)劃分步驟如下。

a.對(duì)通告中提及的發(fā)生洪水災(zāi)害的水庫進(jìn)行文字整理。

b.依據(jù)水庫的洪水災(zāi)害抵御情況，將災(zāi)害后果的定性描述轉(zhuǎn)化為定量描述。

c.根據(jù)定量描述對(duì)水庫等級(jí)進(jìn)行劃分，共分為兩類，一類是近期沒有發(fā)生過潰壩、漫頂?shù)仁录陌踩畮?；另一類是遭受洪水?zāi)害且對(duì)下游有危害的病險(xiǎn)水庫。

由于實(shí)際能得到的有價(jià)值的小型水庫數(shù)據(jù)樣本量較少，最終只得到23個(gè)符合要求的樣本，從而極大地限制了在分類過程中對(duì)訓(xùn)練樣本的需求。為此，采用支持向量機(jī)這種適合小樣本識(shí)別的方法對(duì)小型水庫安全風(fēng)險(xiǎn)等級(jí)進(jìn)行評(píng)估。

3.2.2基于支持向量機(jī)的水庫等級(jí)劃分

考慮使用支持向量機(jī)的方法對(duì)小型水庫等級(jí)進(jìn)行分類，其步驟如下。

a.將數(shù)據(jù)分為訓(xùn)練集和測(cè)試集。其中，70%的數(shù)據(jù)作為訓(xùn)練集，30%的數(shù)據(jù)作為測(cè)試集。

b.實(shí)驗(yàn)分別采用徑向基函數(shù)、線性核函數(shù)和多項(xiàng)式函數(shù)作為核函數(shù)，并在不同的核函數(shù)條件下，選擇最優(yōu)參數(shù)對(duì)小型水庫進(jìn)行評(píng)估。實(shí)驗(yàn)的最優(yōu)參數(shù)、核函數(shù)選擇以及分類的準(zhǔn)確率見表2。

表2 參數(shù)選擇及分類正確率

從表2可以看出，在對(duì)3類已知核函數(shù)的試驗(yàn)中，linear線性核函數(shù)的結(jié)果稍優(yōu)于其他兩類，模型正確率為73.3%，其中懲罰常數(shù)cost為1，gamma為0.08。

為了進(jìn)一步對(duì)模型結(jié)果進(jìn)行分析，可以通過模型維度作圖對(duì)模型中一些變量的分類結(jié)果進(jìn)行描述。

當(dāng)管理單位的工程師人數(shù)在24人以上時(shí)，水庫建成年份和下游影響人口對(duì)小型水庫的風(fēng)險(xiǎn)等級(jí)劃分結(jié)果見圖2。從圖2可以看出，有失事風(fēng)險(xiǎn)的水庫大多建于1960年以前，下游影響的人口數(shù)目多樣；而1970年以后建成的小型水庫失事風(fēng)險(xiǎn)明顯減少。

圖2 建成年份-影響人口分類結(jié)果

多年平均降水量-影響人口分類結(jié)果見圖3。從圖3可以看出當(dāng)?shù)啬甓嗄昶骄邓看笥? 000 mm時(shí)，對(duì)水庫造成的影響較大。結(jié)合樣本所在地實(shí)際情況，水庫所在地區(qū)多年平均降水量多數(shù)分布在 1 500～2 000 mm，因此隨著降水量的增多小型水庫失事風(fēng)險(xiǎn)明顯增多。

圖3 多年平均降水量-影響人口分類結(jié)果

水情監(jiān)測(cè)方式-影響人口分類結(jié)果見圖4(圖中橫坐標(biāo)數(shù)字分別代表：1為無監(jiān)測(cè)設(shè)施；2為人工監(jiān)測(cè)；3為自動(dòng)+人工；4為自動(dòng)監(jiān)測(cè)；5為其他)。通過模型維度作圖得出圖4水情監(jiān)測(cè)方式-影響人口分類結(jié)果，樣本實(shí)際情況基本上以無監(jiān)測(cè)設(shè)施、人工監(jiān)測(cè)、自動(dòng)+人工、自動(dòng)監(jiān)測(cè)這4種監(jiān)測(cè)方式為主，從圖4可看出有失事風(fēng)險(xiǎn)的水庫大多集中于無監(jiān)測(cè)設(shè)施的水庫，該類水庫對(duì)下游影響的人口數(shù)目多樣，部分采取人工監(jiān)測(cè)方式的水庫也存在一定失事風(fēng)險(xiǎn)，采取人工+自動(dòng)以及自動(dòng)監(jiān)測(cè)的水庫基本全部屬于安全水庫。

圖4 水情監(jiān)測(cè)方式-影響人口分類結(jié)果

通過上述不同維度的數(shù)據(jù)分析方式，可以從不同的角度對(duì)小型水庫等級(jí)評(píng)估的影響因素進(jìn)行查看和分析，更加科學(xué)地對(duì)小型水庫安全風(fēng)險(xiǎn)等級(jí)進(jìn)行評(píng)判。

4 結(jié) 論

本文基于小型水庫現(xiàn)狀，綜合考慮小型水庫安全風(fēng)險(xiǎn)的主要影響因子，進(jìn)一步從機(jī)器學(xué)習(xí)角度對(duì)小型水庫安全風(fēng)險(xiǎn)等級(jí)進(jìn)行評(píng)估，探索了支持向量機(jī)在小樣本識(shí)別中的應(yīng)用，構(gòu)建了基于機(jī)器學(xué)習(xí)的小型水庫安全風(fēng)險(xiǎn)等級(jí)評(píng)價(jià)模型。

研究表明，基于支持向量機(jī)的小型水庫風(fēng)險(xiǎn)等級(jí)評(píng)估在小型水庫安全風(fēng)險(xiǎn)識(shí)別上具有一定的適用性，能夠有效解決小型水庫安全等級(jí)評(píng)價(jià)中存在的數(shù)據(jù)樣本較少的問題；與傳統(tǒng)的人工經(jīng)驗(yàn)計(jì)算方法相比，支持向量機(jī)具有簡(jiǎn)便、快捷、經(jīng)濟(jì)、適于推廣的特點(diǎn)。研究成果極大地提升了小型水庫下游山洪預(yù)警的準(zhǔn)確性，實(shí)現(xiàn)了小型水庫重要信息與山洪預(yù)警工作的緊密結(jié)合，為開展小型水庫除險(xiǎn)加固工作提供參考。