鄭斯水,顧圣平,唐鳳珍,張佳丹
(河海大學(xué)水利水電學(xué)院,江蘇 南京 210098)
水電是我國(guó)能源組成中的重要部分。合理進(jìn)行水能開(kāi)發(fā)有利于我國(guó)優(yōu)化能源結(jié)構(gòu)、保障能源安全,以適應(yīng)我國(guó)當(dāng)前可持續(xù)發(fā)展的能源戰(zhàn)略方向。水電站運(yùn)行管理評(píng)價(jià)是對(duì)水電站工作人員的運(yùn)行調(diào)度和管理工作進(jìn)行等級(jí)評(píng)定,并作為水電站運(yùn)行管理考核的依據(jù),可以充分調(diào)動(dòng)水電站工作人員的積極性,有利于挖掘水電站發(fā)電潛力,并保障水電站安全穩(wěn)定運(yùn)行。
近年來(lái),有關(guān)學(xué)者對(duì)水電站運(yùn)行管理評(píng)價(jià)的研究主要針對(duì)水電站經(jīng)濟(jì)運(yùn)行方面,評(píng)價(jià)對(duì)象由單一水電站向梯級(jí)水電站群轉(zhuǎn)變,評(píng)價(jià)指標(biāo)逐漸從絕對(duì)指標(biāo)轉(zhuǎn)向相對(duì)指標(biāo)[1]。水能利用提高率[2]、理想發(fā)電量[3]、可用水能利用率[4]、發(fā)電完成率[5]等指標(biāo)先后被用于水電站經(jīng)濟(jì)運(yùn)行評(píng)價(jià)。由于各水電站庫(kù)容和機(jī)組性能的差異性、入庫(kù)徑流的隨機(jī)性、電網(wǎng)要求的不確定性等眾多因素的影響,水電站運(yùn)行調(diào)度呈現(xiàn)非結(jié)構(gòu)化的特點(diǎn)[6],因此對(duì)于水電站運(yùn)行管理評(píng)價(jià)須考慮其相對(duì)性、隨機(jī)性和模糊性?,F(xiàn)有研究對(duì)于水電站運(yùn)行的隨機(jī)性和模糊性的考慮主要體現(xiàn)在水庫(kù)優(yōu)化調(diào)度[7]和水電開(kāi)發(fā)風(fēng)險(xiǎn)評(píng)價(jià)[8]等方面,而將兩者結(jié)合考慮對(duì)水電站運(yùn)行管理行為進(jìn)行系統(tǒng)評(píng)價(jià)的研究目前還比較少見(jiàn)?;诖耍疚慕⑺娬具\(yùn)行管理綜合評(píng)價(jià)指標(biāo)體系,并利用云模型方法建立評(píng)價(jià)模型,將隨機(jī)性和模糊性有機(jī)結(jié)合,以實(shí)現(xiàn)對(duì)水電站運(yùn)行管理評(píng)價(jià)工作的改進(jìn)。
李德毅等[9]于1995年在隸屬函數(shù)概念的基礎(chǔ)上提出云模型,并被廣泛應(yīng)用于系統(tǒng)評(píng)價(jià)、算法改進(jìn)、數(shù)據(jù)挖掘等[10-12]。云是分析論域X在模糊集合A上的隸屬度μ的模型,在評(píng)價(jià)模型中,X即為評(píng)價(jià)對(duì)象集合,A為評(píng)語(yǔ)集,μ為X中元素隸屬于A的概率。云模型的特征參數(shù)主要有期望Ex、熵En、超熵He,分別表征概念的最典型點(diǎn)、概念的不確定性、熵的不確定性,分別反映各指標(biāo)的均值、各指標(biāo)值變化的可接受范圍、云滴的離散程度。
正態(tài)概率分布模型和正態(tài)隸屬函數(shù)的普適性,共同奠定了正云模型的普適性[13]。正態(tài)云有別于正態(tài)分布之處在于其特征參數(shù)超熵。正態(tài)云屬于泛正態(tài)分布,當(dāng)He=0時(shí),正態(tài)云就退化為正態(tài)分布。根據(jù)正態(tài)分布理論,若所有云滴表示的事件發(fā)生概率和為1,則有99.7%的云滴落于區(qū)間[Ex-3En,Ex+3En]上,相應(yīng)區(qū)間外的云滴視為小概率事件,一般不考慮,依此進(jìn)行云滴檢驗(yàn)。
云發(fā)生器是定性概念和定量概念之間相互轉(zhuǎn)化的算法,分為正向云發(fā)生器和逆向云發(fā)生器。正向云發(fā)生器已知云模型Ex、En、He3個(gè)參數(shù),生成一定數(shù)量的云滴,是定性到定量的過(guò)程;逆向云發(fā)生器為根據(jù)已知云滴分布確定3個(gè)參數(shù),是定量到定性的過(guò)程。
云相似度算法[14]用于計(jì)算相同論域內(nèi)兩個(gè)云的重合度,本文應(yīng)用于衡量指標(biāo)云對(duì)于評(píng)價(jià)等級(jí)云的歸屬度,以確定該指標(biāo)的評(píng)價(jià)等級(jí)。指標(biāo)云和評(píng)價(jià)等級(jí)云分別是指標(biāo)云化和評(píng)語(yǔ)集云化的結(jié)果,采用正向云發(fā)生器。
為更加系統(tǒng)、客觀地反映不同水電站在不同來(lái)水條件下進(jìn)行水電運(yùn)行的管理水平,評(píng)價(jià)指標(biāo)體系的構(gòu)建遵循以下原則[15]:①系統(tǒng)性,從多方面考慮,盡可能全面描述運(yùn)行管理效果水平;②客觀性,盡可能選取定量的可精確計(jì)算的指標(biāo);③相對(duì)性,盡量過(guò)濾水電規(guī)模和可用水量等客觀因素的影響;④可行性,指標(biāo)值計(jì)算方便、可行。
依據(jù)以上原則構(gòu)建水電站運(yùn)行管理綜合評(píng)價(jià)體系,一級(jí)指標(biāo)主要考慮經(jīng)濟(jì)和安全兩方面。①經(jīng)濟(jì)指標(biāo)主要針對(duì)水能的利用效率,包括水能利用提高率、棄水日均發(fā)電負(fù)荷率、相對(duì)發(fā)電耗水率、可用水量利用率、計(jì)劃報(bào)表合格率;②安全指標(biāo)包括機(jī)組運(yùn)行安全、壩體安全、下游防洪安全以及閘門(mén)控制安全等方面,包括洪水預(yù)報(bào)準(zhǔn)確率、下游防洪安全性、壩體安全性、機(jī)組工作安全性、閘門(mén)控制可靠性等,見(jiàn)表1。其中定性指標(biāo)的評(píng)定標(biāo)準(zhǔn)為:0.85~1為完全安全;0.70~0.85為基本安全;0.50~0.70為存在隱患但未造成損失;0~0.5為出現(xiàn)問(wèn)題并造成損失。
表1 水電站運(yùn)行管理綜合評(píng)價(jià)指標(biāo)體系
指標(biāo)云化是計(jì)算每個(gè)評(píng)價(jià)指標(biāo)的Ex、En、He等3個(gè)云模型特征參數(shù),從而將指標(biāo)特征通過(guò)云圖進(jìn)行具體描繪。收集待評(píng)單位評(píng)價(jià)期內(nèi)各指標(biāo)的有效數(shù)據(jù),假設(shè)指標(biāo)xj在評(píng)價(jià)期內(nèi)共有n個(gè)相互獨(dú)立的數(shù)據(jù)xji,其云模型參數(shù)Cj(Exj,Enj,Hej)滿足
(1)
式中:Exj為指標(biāo)xj的期望;Enj為指標(biāo)xj的熵;Hej為指標(biāo)xj的超熵。
實(shí)際計(jì)算中,如果缺乏某些指標(biāo)的具體數(shù)據(jù)資料,但可以獲得其指標(biāo)期望數(shù)據(jù)時(shí),熵和超熵依據(jù)經(jīng)驗(yàn)公式計(jì)算或根據(jù)經(jīng)驗(yàn)擬定[16]。
對(duì)所有指標(biāo)均給予優(yōu)秀、良好、合格、較差、極差5個(gè)評(píng)價(jià)等級(jí),根據(jù)數(shù)據(jù)統(tǒng)計(jì)和專家分析,確定指標(biāo)xj評(píng)價(jià)等級(jí)k(k=1,2,…,5)對(duì)應(yīng)的Exjk,則指標(biāo)xj評(píng)價(jià)等級(jí)k對(duì)應(yīng)的評(píng)價(jià)等級(jí)云CRjk(Exjk,Enjk,Hejk)滿足
(2)
式中:Enjk為指標(biāo)xj評(píng)語(yǔ)等級(jí)k對(duì)應(yīng)的熵值;Exj(k-1)、Exjk、Exj(k+1)分別為指標(biāo)xj比等級(jí)k低一等級(jí)評(píng)語(yǔ)、評(píng)語(yǔ)等級(jí)k、比等級(jí)k高一等級(jí)評(píng)語(yǔ)對(duì)應(yīng)的期望;Hejk為指標(biāo)xj評(píng)語(yǔ)等級(jí)k對(duì)應(yīng)的超熵。
根據(jù)對(duì)若干水電站各指標(biāo)的統(tǒng)計(jì)結(jié)果和專家估計(jì),確定各指標(biāo)各個(gè)評(píng)價(jià)等級(jí)的期望值,進(jìn)而根據(jù)式(2)確定熵和超熵的值。所得各指標(biāo)各評(píng)價(jià)等級(jí)云對(duì)應(yīng)的云模型特征參數(shù)如表2所示。
采用云相似度算法計(jì)算xj指標(biāo)云Cj(Exj,Enj,Hej)對(duì)于該指標(biāo)評(píng)價(jià)等級(jí)k對(duì)應(yīng)的評(píng)價(jià)等級(jí)云CRjk(Exjk,Enjk,Hejk)的相似度,記為δjk,步驟如下:
表2 水電站運(yùn)行管理各指標(biāo)評(píng)價(jià)等級(jí)云模型特征參數(shù)
步驟4:計(jì)算
(3)
步驟5:重復(fù)步驟1~4,直至產(chǎn)生足夠大數(shù)量的ηjk,其個(gè)數(shù)設(shè)為M個(gè)。
步驟6:計(jì)算
(4)
根據(jù)Cj對(duì)于各CRjk的δjk,可計(jì)算該指標(biāo)隸屬于該等級(jí)的概率,指標(biāo)xj隸屬于等級(jí)k的概率pjk滿足
(5)
若評(píng)價(jià)指標(biāo)總共有m個(gè),可采用層次分析法(AHP)和熵權(quán)分析法(EWM)分別計(jì)算指標(biāo)的主客觀權(quán)重,在此基礎(chǔ)上計(jì)算組合權(quán)重。該方法既考慮了評(píng)價(jià)者對(duì)水電站運(yùn)行管理評(píng)價(jià)指標(biāo)的偏好,也兼顧了對(duì)數(shù)據(jù)本身內(nèi)在屬性的挖掘。xj的組合權(quán)重值滿足
(6)
式中:wj為xj的組合權(quán)重值;wAj為xj通過(guò)AHP法得到的主觀權(quán)重值;wEj為xj通過(guò)EWM法得到的客觀權(quán)重值;m為評(píng)價(jià)指標(biāo)個(gè)數(shù)。
若m個(gè)評(píng)價(jià)指標(biāo)的云模型特征參數(shù)和權(quán)重分別為Cj(Exj,Enj,Hej)和wj(j=1,2,…,m),則綜合指標(biāo)云Cz(Exz,Enz,Hez)滿足[17]
(7)
式中:Exz、Enz和Hez分別為綜合指標(biāo)云的期望、熵和超熵。
綜合等級(jí)評(píng)價(jià)分為優(yōu)秀、良好、合格、較差、極差5個(gè)等級(jí)。在已知m個(gè)指標(biāo)k等級(jí)的評(píng)價(jià)等級(jí)云和權(quán)重分別為CRjk(Exjk,Enjk,Hejk)和wj的基礎(chǔ)上,采用與式(5)相同的方式計(jì)算綜合評(píng)價(jià)等級(jí)為k的評(píng)價(jià)等級(jí)云CRzk(Exzk,Enzk,Hezk),并計(jì)算Cz對(duì)于各個(gè)CRzk的云相似度δzk,從而得到水電站運(yùn)行管理綜合評(píng)價(jià)為k等級(jí)的概率pzk。選用其中概率最大值pzk_max=pzL對(duì)應(yīng)的評(píng)價(jià)等級(jí)L作為水電站運(yùn)行管理綜合評(píng)價(jià)等級(jí)。
以我國(guó)西南某省7個(gè)水電站(以水電站Ⅰ~Ⅶ表示)為例,應(yīng)用云模型方法對(duì)其運(yùn)行管理進(jìn)行綜合評(píng)價(jià),各水電站部分特征參數(shù)如表3所示。
通過(guò)對(duì)定量指標(biāo)的數(shù)據(jù)統(tǒng)計(jì)以及專家對(duì)定性指標(biāo)的獨(dú)立打分,根據(jù)式(1)以及相關(guān)經(jīng)驗(yàn)公式,得到7個(gè)水電站10個(gè)指標(biāo)的云模型特征參數(shù),結(jié)果見(jiàn)表4。定量指標(biāo)數(shù)據(jù)主要來(lái)自文獻(xiàn)[18-23];當(dāng)部分?jǐn)?shù)據(jù)資料缺乏,但可以獲得其指標(biāo)期望數(shù)據(jù)時(shí),熵和超熵依據(jù)經(jīng)驗(yàn)公式計(jì)算。定性指標(biāo)值由專家獨(dú)立評(píng)價(jià)得到。
表3 水電站Ⅰ~Ⅶ特征參數(shù)
表4 水電站Ⅰ~Ⅶ運(yùn)行管理評(píng)價(jià)指標(biāo)云參數(shù)
表5 水電站Ⅰ運(yùn)行管理指標(biāo)評(píng)價(jià)云相似度及評(píng)價(jià)等級(jí)概率 單位:%
根據(jù)本文方法,分別計(jì)算各水電站各個(gè)指標(biāo)云對(duì)評(píng)價(jià)等級(jí)云的云相似度,并計(jì)算各指標(biāo)隸屬于各等級(jí)的概率。以水電站Ⅰ各指標(biāo)為例,計(jì)算結(jié)果如表5所示。水電站ⅠA2、A3、A4、B2、B35個(gè)指標(biāo)均以95%以上概率隸屬于優(yōu)秀等級(jí);A1指標(biāo)分別以68.69%和31.06%的概率隸屬于優(yōu)秀和良好,A5指標(biāo)分別以43.18%、43.37%和11.32%的概率隸屬于合格、良好和優(yōu)秀;B1指標(biāo)分別以35.11%和64.89%的概率隸屬于良好和優(yōu)秀;B4指標(biāo)以98.35%的概率隸屬于合格;B5指標(biāo)分別以81.49%和18.51%的概率隸屬于較差和合格。鑒于水電站ⅠB4和B5兩個(gè)指標(biāo)評(píng)價(jià)等級(jí)較低,建議其在閘門(mén)安全控制和機(jī)組安全工作等方面加強(qiáng)管理。
根據(jù)AHP和EWM方法分別計(jì)算各指標(biāo)主觀權(quán)重wA和客觀權(quán)重wE,并根據(jù)式(4)計(jì)算組合權(quán)重w,結(jié)果見(jiàn)表6。由AHP方法一致性分析得到CI=0.028 8,CR=0.019 3<0.1,故認(rèn)為專家判斷矩陣一致性較好。且經(jīng)濟(jì)性指標(biāo)占47.4%,安全性指標(biāo)占52.6%,較為符合實(shí)際的決策情況。
表6 水電站運(yùn)行管理評(píng)價(jià)指標(biāo)權(quán)重分析
根據(jù)式(5)分別計(jì)算水電站Ⅰ~Ⅶ運(yùn)行管理綜合評(píng)價(jià)指標(biāo)云模型特征參數(shù)和評(píng)價(jià)等級(jí)云參數(shù),結(jié)果見(jiàn)表7和表8。
表7 水電站Ⅰ~Ⅶ運(yùn)行管理評(píng)價(jià)綜合指標(biāo)云參數(shù)
表8 水電站Ⅰ~Ⅶ運(yùn)行管理評(píng)價(jià)評(píng)價(jià)等級(jí)云模型特征參數(shù)
表9 水電站Ⅰ~Ⅶ運(yùn)行管理評(píng)價(jià)云相似度及評(píng)價(jià)等級(jí)概率 單位:%
分別計(jì)算各水電站綜合指標(biāo)云對(duì)各個(gè)等級(jí)的評(píng)價(jià)等級(jí)云的云相似度以及隸屬于各評(píng)價(jià)等級(jí)概率,結(jié)果見(jiàn)表9。以概率最大的隸屬等級(jí)作為評(píng)價(jià)結(jié)果,水電站Ⅰ、Ⅲ、Ⅳ、Ⅴ、Ⅵ、Ⅶ分別以86.78%、62.54%、96.77%、90.35%、58.08%、84.92%的概率隸屬于良好等級(jí),其中水電站Ⅰ隸屬于優(yōu)秀的概率最大,為13.13%,水電站Ⅵ隸屬于合格的概率最大,為41.55%,水電站Ⅱ以86.13%的概率隸屬于合格,其評(píng)價(jià)等級(jí)明顯低于其他水電站。
為更加直觀地體現(xiàn)各水電站間的差異,運(yùn)用Matlab繪圖工具繪制其綜合指標(biāo)云和評(píng)價(jià)等級(jí)云,每一簇云均由各自的2 000云滴組成,由于部分指標(biāo)優(yōu)秀等級(jí)的期望達(dá)到可能最大值,故優(yōu)秀等級(jí)云采用半升云形式,其余等級(jí)均采用完整的正態(tài)云形式。圖1為水電站Ⅰ~Ⅶ的運(yùn)行管理綜合評(píng)價(jià)結(jié)果云圖。從圖1可以看出,水電站Ⅰ期望值略高于良好等級(jí),電站Ⅳ的云圖與良好等級(jí)云圖基本吻合;水電站Ⅱ、Ⅲ、Ⅴ、Ⅵ、Ⅶ的期望值均介于良好與合格之間,水電站Ⅴ、Ⅶ的期望更靠近良好等級(jí),電站Ⅲ、Ⅵ的期望值處于良好和合格的中間位置,電站Ⅱ的期望更偏于合格等級(jí)。
云模型是在概率論和模糊集合理論交叉滲透的基礎(chǔ)上發(fā)展而來(lái),將隨機(jī)性和模糊性有機(jī)結(jié)合,從而實(shí)現(xiàn)了定量數(shù)據(jù)和定性概念之間的有效轉(zhuǎn)化。本文應(yīng)用云模型方法對(duì)7個(gè)水電站運(yùn)行管理進(jìn)行評(píng)價(jià),得到7個(gè)水電站各個(gè)評(píng)價(jià)指標(biāo)隸屬于各個(gè)等級(jí)的概率,并且繪制各水電站運(yùn)行管理綜合評(píng)價(jià)結(jié)果云圖,對(duì)水電站運(yùn)行管理效果有相對(duì)直觀的展示。若采用典型的綜合評(píng)價(jià)方法如模糊綜合評(píng)價(jià)法,則不能有效處理評(píng)價(jià)因素的隨機(jī)性,具體表現(xiàn)為其結(jié)果僅能得到各水電站在各個(gè)指標(biāo)上所隸屬的等級(jí)。而水電站運(yùn)行管理效果隨著各種因素的變化,理應(yīng)是波動(dòng)
(a) 水電站Ⅰ
(b) 水電站Ⅱ
(d) 水電站Ⅳ
(f) 水電站Ⅵ
不定的,這在云模型的結(jié)果上得到了更好地反映。兩者對(duì)比之下,對(duì)于具有模糊性和隨機(jī)性特點(diǎn)的水電站運(yùn)行管理綜合評(píng)價(jià),云模型具有更好的適用性。
本文建立水電站運(yùn)行管理綜合評(píng)價(jià)指標(biāo)體系,并將云模型方法引入水電站運(yùn)行管理評(píng)價(jià)研究中。應(yīng)用云模型方法進(jìn)行評(píng)價(jià),有效反映了水電站運(yùn)行過(guò)程中的模糊性和隨機(jī)性。并結(jié)合云相似度算法和評(píng)價(jià)等級(jí)概率公式,避免了直觀判定可能存在的誤差。將該方法應(yīng)用到7個(gè)水電站的運(yùn)行管理評(píng)價(jià),得到水電站運(yùn)行管理各指標(biāo)評(píng)價(jià)等級(jí)及概率、綜合評(píng)價(jià)等級(jí)及概率。結(jié)果表明了該模型方法的適用性。
對(duì)于有生態(tài)需求或其他社會(huì)效益顯著的水電站,應(yīng)用云模型方法進(jìn)行水電站運(yùn)行管理綜合評(píng)價(jià)時(shí),可考慮增加如庫(kù)區(qū)水質(zhì)合格率、下游生態(tài)流量保證率等生態(tài)指標(biāo),增加如調(diào)峰任務(wù)完成率、調(diào)頻任務(wù)完成率、間接防洪效益等社會(huì)指標(biāo)。