劉禹良，郭劍峰，鄒 燕

(1.重慶電力高等?？茖W(xué)校，重慶400053;2.四川省電力公司內(nèi)江電業(yè)局，四川內(nèi)江641000)

0 引言

國內(nèi)外已在配電網(wǎng)可靠性評(píng)估方面取得較大進(jìn)展，常用的方法是故障模式后果分析法(FMEA，F(xiàn)ailure mode and effect analysis)［1］，此外還有故障遍歷算法［2］、最小路法［3］、最小割集法［4］和分塊算法［5］等。但上述算法均基于元件可靠性參數(shù)為確定值的假設(shè)。實(shí)際上，原始參數(shù)的數(shù)值及其概率分布可能會(huì)因?yàn)榻y(tǒng)計(jì)資料不足或統(tǒng)計(jì)誤差以及對(duì)電網(wǎng)未來運(yùn)行條件和環(huán)境預(yù)測(cè)不足而具有不確定性［6］。為使評(píng)估結(jié)果能更準(zhǔn)確地反映系統(tǒng)實(shí)際運(yùn)行的情況，應(yīng)充分考慮參數(shù)不確定性的影響［7］。近年來，區(qū)間分析方法將可靠性參數(shù)作為區(qū)間數(shù)處理，可通過一次運(yùn)算評(píng)估負(fù)荷點(diǎn)和系統(tǒng)的區(qū)間可靠性指標(biāo)，且從全局角度實(shí)現(xiàn)可靠性評(píng)估［8］。

貝葉斯網(wǎng)絡(luò)［9］是20世紀(jì)80年代后期新興起的一種人工智能方法，基于其推理模式能辨識(shí)系統(tǒng)薄弱環(huán)節(jié)，并提出有針對(duì)性的改進(jìn)措施［10］。貝葉斯網(wǎng)絡(luò)已應(yīng)用于配電網(wǎng)可靠性評(píng)估，提高了配電網(wǎng)的經(jīng)濟(jì)效益。文獻(xiàn)［11］基于故障樹分析法建立貝葉斯網(wǎng)絡(luò)，實(shí)現(xiàn)配電網(wǎng)的可靠性評(píng)估。文獻(xiàn)［12］以配電網(wǎng)各元件、分塊、負(fù)荷點(diǎn)及系統(tǒng)狀態(tài)為節(jié)點(diǎn)，以它們之間的影響關(guān)系為有向邊，建立配電網(wǎng)可靠性評(píng)估的4層貝葉斯網(wǎng)絡(luò)，實(shí)現(xiàn)配電網(wǎng)的可靠性計(jì)算和薄弱環(huán)節(jié)辨識(shí)。但上述算法均未計(jì)及元件可靠性參數(shù)的不確定性，且未計(jì)及開關(guān)設(shè)備的故障。

基于貝葉斯網(wǎng)絡(luò)的雙向推理模式，且計(jì)及參數(shù)的不確定性，本文提出計(jì)及開關(guān)故障的配電網(wǎng)可靠性評(píng)估的貝葉斯網(wǎng)絡(luò)區(qū)間算法，能更加高效地進(jìn)行參數(shù)為區(qū)間數(shù)的配電網(wǎng)可靠性計(jì)算和薄弱環(huán)節(jié)辨識(shí)。

1 區(qū)間數(shù)的基本運(yùn)算

本文主要研究非退化正區(qū)間數(shù)［13］，利用其進(jìn)行配電網(wǎng)可靠性評(píng)估時(shí)，可能涉及其四則運(yùn)算和排序相關(guān)知識(shí)。其加法、減法和乘法運(yùn)算均采用文獻(xiàn)［14］的定義，排序則采用文獻(xiàn)［15］定義的矩陣排序法。

配電網(wǎng)可靠性評(píng)估中，可靠性指標(biāo)計(jì)算的區(qū)間算法涉及正區(qū)間數(shù)的除法。由于區(qū)間除法運(yùn)算存在“相關(guān)性問題”［16］，其不再是區(qū)間乘法的逆運(yùn)算，且會(huì)引起運(yùn)算結(jié)果超寬度。為解決上述問題，須對(duì)參與除法運(yùn)算的區(qū)間數(shù)標(biāo)準(zhǔn)化［13］，并用重新定義的區(qū)間除法進(jìn)行計(jì)算［17］。

設(shè)［x］ = ［xmin，xmax］、［y］= ［ymin，ymax］均為非退化正區(qū)間數(shù)，將其化為標(biāo)準(zhǔn)式為:

其中

令［x1］=［μ1，1］，［x2］=［μ2，1］，文獻(xiàn)［15］稱之為基本區(qū)間數(shù)，并定義新的除法公式(用■表示)為:

利用文［15］的方法進(jìn)行區(qū)間數(shù)除法運(yùn)算可以得到正確合理的結(jié)果，但步驟較繁。實(shí)際上，聯(lián)立式(1)～(4)，可得非退化正區(qū)間數(shù)除法的較為簡(jiǎn)單的公式:

2 配電網(wǎng)可靠性評(píng)估貝葉斯網(wǎng)絡(luò)區(qū)間算法

貝葉斯網(wǎng)絡(luò)是隨機(jī)事件間概率關(guān)系的圖形表示［12］，基于其前向推理模式，可計(jì)算出負(fù)荷點(diǎn)和系統(tǒng)的區(qū)間可靠性指標(biāo);基于其后向推理模式，可通過元件的可靠性參數(shù)和負(fù)荷點(diǎn)的節(jié)點(diǎn)類型辨識(shí)負(fù)荷點(diǎn)的薄弱環(huán)節(jié)，并可通過各元件分?jǐn)偟膮^(qū)間ENS辨識(shí)系統(tǒng)的薄弱環(huán)節(jié)。

圖1為一個(gè)簡(jiǎn)單中壓配電網(wǎng)示意圖。為簡(jiǎn)化起見，本文以該配電網(wǎng)為例，建立其可靠性評(píng)估的貝葉斯網(wǎng)絡(luò)，并在此基礎(chǔ)上建立其可靠性計(jì)算及薄弱環(huán)節(jié)辨識(shí)的區(qū)間模型。

圖1 簡(jiǎn)單配電網(wǎng)示意圖

2.1 貝葉斯網(wǎng)絡(luò)的建立

基于文獻(xiàn)［12］所述流程，建立圖1配電網(wǎng)可靠性評(píng)估的四層貝葉斯網(wǎng)絡(luò)，如圖2所示:

2.2 基于貝葉斯網(wǎng)絡(luò)的可靠性計(jì)算及薄弱環(huán)節(jié)分析

若用“1”表示元件或分塊的故障狀態(tài)及負(fù)荷點(diǎn)的停運(yùn)狀態(tài)、“0”表示正常狀態(tài)，則圖2中的貝葉斯網(wǎng)絡(luò)是一個(gè)典型的“或”節(jié)點(diǎn)模型。這是由于分塊中任意一個(gè)元件故障，均可造成其所在分塊的元件停運(yùn);同樣地，對(duì)某負(fù)荷點(diǎn)的供電有影響的任一分塊或開關(guān)元件故障，均可造成該負(fù)荷點(diǎn)停運(yùn)。

2.2.1 分塊的區(qū)間可靠性參數(shù)計(jì)算

圖2 圖1配電網(wǎng)可靠性評(píng)估的貝葉斯網(wǎng)絡(luò)

基于以上分析，利用貝葉斯網(wǎng)絡(luò)的前向推理功能，并引入第1節(jié)所述的區(qū)間數(shù)四則運(yùn)算法則，采用多元件串聯(lián)公式計(jì)算分塊的等效區(qū)間可靠性參數(shù)，如下:

式中，sub為待分析的分塊;［λi］、［γi］分別為元件 i的區(qū)間故障率、區(qū)間修復(fù)時(shí)間;［λsub］、［γsub］分別為分塊的區(qū)間等值故障率、區(qū)間等值修復(fù)時(shí)間;Ns為分塊包含的元件數(shù)。

分塊的區(qū)間等值修復(fù)時(shí)間的計(jì)算涉及區(qū)間除法運(yùn)算，為避免相關(guān)性問題導(dǎo)致運(yùn)算結(jié)果超寬度［16］，本文利用第1節(jié)的分析求解式(8)。

2.2.2 負(fù)荷點(diǎn)的區(qū)間可靠性指標(biāo)計(jì)算及薄弱環(huán)節(jié)分析

同上述分析類似，得負(fù)荷點(diǎn)的區(qū)間可靠性指標(biāo)為:

式中，LP 為待分析的負(fù)荷點(diǎn);［λsubi］、［λj］分別為對(duì)負(fù)荷點(diǎn)的供電有影響的塊i和開關(guān)j的區(qū)間等值故障率;［γLP-subi］、［γLP-j］分別為塊 i和開關(guān) j故障時(shí)負(fù)荷點(diǎn)的區(qū)間缺電時(shí)間，其取值決定于其故障時(shí)負(fù)荷點(diǎn)的節(jié)點(diǎn)類型;［λLP］、［γLP］分別為負(fù)荷點(diǎn)的區(qū)間故障率和區(qū)間修復(fù)時(shí)間;Nb為對(duì)負(fù)荷點(diǎn)有影響的分塊數(shù)目;Nq為對(duì)負(fù)荷點(diǎn)有影響的開關(guān)數(shù)目。

利用貝葉斯網(wǎng)絡(luò)的后向推理模式辨識(shí)負(fù)荷點(diǎn)的薄弱環(huán)節(jié)，即:若觀測(cè)到系統(tǒng)中某負(fù)荷點(diǎn)缺電，反過

根據(jù)圖2貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)，知線路L1屬于塊1，而塊1故障會(huì)引起負(fù)荷點(diǎn) LP1停運(yùn)，因此上式中的P(LP1=1|L1=1)=1。式(12)可寫為:來分析系統(tǒng)中哪個(gè)元件對(duì)該事件影響最大。

以負(fù)荷點(diǎn)LP1為例進(jìn)行說明。已知負(fù)荷點(diǎn)LP1停運(yùn)，其由線路L1故障引起的后驗(yàn)概率［12］為:

當(dāng)元件可靠性參數(shù)為區(qū)間數(shù)時(shí)，根據(jù)第1節(jié)的分析，得區(qū)間后驗(yàn)概率為:

其中，

分別為L(zhǎng)1故障對(duì)LP1停運(yùn)的區(qū)間貢獻(xiàn)概率和LP1的區(qū)間停運(yùn)概率，［γLP1-L1］為L(zhǎng)1故障引起LP1的區(qū)間停運(yùn)時(shí)間。

對(duì)其他對(duì)LP1供電有影響的元件，都能得出類似(14)～(16)的式子。

由上可知:區(qū)間后驗(yàn)概率取決于元件的區(qū)間可靠性參數(shù)及負(fù)荷點(diǎn)的節(jié)點(diǎn)類型。計(jì)算每個(gè)元件故障引起某一負(fù)荷點(diǎn)停運(yùn)的區(qū)間后驗(yàn)概率，并按文獻(xiàn)［15］所述方法排序，其中最大者對(duì)應(yīng)的元件故障最可能造成該負(fù)荷點(diǎn)停運(yùn)，其為該負(fù)荷點(diǎn)的薄弱環(huán)節(jié)。

2.2.3 系統(tǒng)的區(qū)間可靠性指標(biāo)計(jì)算及薄弱環(huán)節(jié)分析

在計(jì)算負(fù)荷點(diǎn)的區(qū)間可靠性指標(biāo)的基礎(chǔ)上，可進(jìn)一步計(jì)算系統(tǒng)的區(qū)間可靠性指標(biāo)。常見系統(tǒng)指標(biāo)計(jì)算公式的區(qū)間形式如下:

①系統(tǒng)平均停電頻率指標(biāo)SAIFI(System average interruption frequency index)

其中［λi］和Ni分別為負(fù)荷點(diǎn)i的用戶區(qū)間停運(yùn)率和用戶數(shù)，NLP為系統(tǒng)負(fù)荷點(diǎn)總數(shù)。

②系統(tǒng)平均停電持續(xù)時(shí)間SAIDI(System average interruption duration index)

其中［Ui］= ［λi］［γi］，為負(fù)荷點(diǎn)i的區(qū)間年平均停運(yùn)時(shí)間，小時(shí)/年。

③用戶平均停電持續(xù)時(shí)間CAIDI(Customer average interruption duration index)

④平均供電可用率指標(biāo)ASAI(Average service availability index)

⑤系統(tǒng)缺供電量指標(biāo)ENS(Energy not supplied)

Lai為負(fù)荷點(diǎn)i的平均負(fù)荷(kW)。

由前述分析知:后向推理只能辨識(shí)負(fù)荷點(diǎn)的薄弱環(huán)節(jié)，而不能從全局的角度辨識(shí)整個(gè)配電網(wǎng)的薄弱環(huán)節(jié)。故選擇配電網(wǎng)的ENS指標(biāo)進(jìn)行分析，根據(jù)各元件分?jǐn)偟膮^(qū)間ENS確定配電網(wǎng)的薄弱環(huán)節(jié)。

配電網(wǎng)中第j個(gè)元件分?jǐn)偟膮^(qū)間ENS為:

式中，Lai為第i個(gè)負(fù)荷點(diǎn)的負(fù)荷，［λj］為第 j個(gè)元件的區(qū)間故障率，［γi-j］為僅第j個(gè)元件故障時(shí)第i個(gè)負(fù)荷點(diǎn)的區(qū)間缺電時(shí)間。

以線路L3為例進(jìn)行分析。根據(jù)圖2貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)，確定L3屬于塊2，其故障時(shí)負(fù)荷點(diǎn)LP2和LP3的缺電時(shí)間均為L(zhǎng)3的修復(fù)時(shí)間，LP1的缺電時(shí)間為塊2的前向分段開關(guān)s1的操作時(shí)間，LP4的缺電時(shí)間為切換開關(guān)c1的操作時(shí)間。L3分?jǐn)偟膮^(qū)間ENS為:

根據(jù)式(23)計(jì)算每個(gè)元件分?jǐn)偟膮^(qū)間ENS，并按文［15］所述方法排序，分?jǐn)偭孔畲笳邽橄到y(tǒng)中最薄弱的元件。

3 算例分析

利用以上算法對(duì)RBTS系統(tǒng)母線2的饋線F3進(jìn)行區(qū)間可靠性計(jì)算及薄弱環(huán)節(jié)分析。

假設(shè)各線路、變壓器的故障率和修復(fù)時(shí)間、分段開關(guān)和切換開關(guān)的操作時(shí)間在文獻(xiàn)［18］原始參數(shù)的±10%范圍內(nèi)變化，斷路器、分段開關(guān)和熔斷器的故障率為［0.004，0.008］次/年，修復(fù)時(shí)間為［3，5］小時(shí)/次。按文［18］的第5種接線方式進(jìn)行計(jì)算，即主饋線上配置分段開關(guān)、分支線保護(hù)、備用電源和備用變壓器。系統(tǒng)的接線圖如圖3所示。

圖3 RBTS母線2的系統(tǒng)接線圖

3.1 負(fù)荷點(diǎn)區(qū)間可靠性計(jì)算及薄弱環(huán)節(jié)分析

利用式(6)～(8)，對(duì)計(jì)及開關(guān)故障的情形計(jì)算各負(fù)荷點(diǎn)的區(qū)間可靠性指標(biāo)，結(jié)果如表1所示:

表1 RBTS-BUS2饋線3負(fù)荷點(diǎn)的區(qū)間可靠性指標(biāo)

從以上結(jié)果可知:區(qū)間算法考慮了元件可靠性參數(shù)的可能波動(dòng)范圍，可同時(shí)給出最好和最差兩種情況下的負(fù)荷點(diǎn)可靠性指標(biāo)。

應(yīng)用2.2.2節(jié)負(fù)荷點(diǎn)的薄弱環(huán)節(jié)辨識(shí)技術(shù)，分析該配電網(wǎng)的負(fù)荷點(diǎn)薄弱環(huán)節(jié)。

以負(fù)荷點(diǎn)LP10為例，利用式(14)計(jì)算各元件故障引起其停運(yùn)的區(qū)間后驗(yàn)概率，列于表2。

表2 RBTS-BUS2饋線3各元件故障引起負(fù)荷點(diǎn)10停運(yùn)的區(qū)間后驗(yàn)概率

利用文獻(xiàn)［15］所述的區(qū)間數(shù)矩陣排序法對(duì)上述計(jì)算結(jié)果進(jìn)行排序，可知:元件L16對(duì)應(yīng)的區(qū)間后驗(yàn)概率最大，其為負(fù)荷點(diǎn)LP10的薄弱環(huán)節(jié)。

利用上述方法可辨識(shí)該配電網(wǎng)所有負(fù)荷點(diǎn)的薄弱環(huán)節(jié)，列于表3。

在實(shí)際配電網(wǎng)中，可利用上述方法識(shí)別重要負(fù)荷點(diǎn)的薄弱環(huán)節(jié)，采取相應(yīng)改善措施以提高其可靠性。

3.2 系統(tǒng)區(qū)間可靠性計(jì)算及薄弱環(huán)節(jié)分析

對(duì)計(jì)及開關(guān)故障的情形進(jìn)行系統(tǒng)區(qū)間可靠性計(jì)算，其結(jié)果如表4所示。

表3 RBTS-BUS2饋線3各負(fù)荷點(diǎn)的薄弱環(huán)節(jié)

表4 計(jì)及開關(guān)故障的RBTS-BUS2饋線3的區(qū)間可靠性

從表4可知:區(qū)間算法考慮了元件可靠性參數(shù)的可能波動(dòng)范圍，可同時(shí)給出最好和最差兩種情況下的系統(tǒng)可靠性指標(biāo)。

根據(jù)2.2.3節(jié)分析，將計(jì)及開關(guān)故障時(shí)每個(gè)元件分?jǐn)偟膮^(qū)間ENS列于表5。

表5 RBTS-BUS2饋線3各元件分?jǐn)偟膮^(qū)間ENS

利用文獻(xiàn)［15］所述的區(qū)間數(shù)矩陣排序法對(duì)上述計(jì)算結(jié)果進(jìn)行排序，可知:干線L16、L18、L21分?jǐn)偟腅NS較多。其中，L18分?jǐn)偟腅NS最多，其為配電網(wǎng)的最薄弱環(huán)節(jié)。

表5的計(jì)算結(jié)果表明，L16是比L25更加薄弱的環(huán)節(jié)。為驗(yàn)證上述算法的正確性，將L16和L25的故障率分別置0并計(jì)算系統(tǒng)的區(qū)間可靠性指標(biāo)，結(jié)果如表6所示。

表6 RBTS-BUS2饋線3將L16和L25的故障率分別置0時(shí)的區(qū)間可靠性指標(biāo)

將表6與表4的結(jié)果對(duì)比可知:改善L16的可靠性比改善L25的可靠性對(duì)提高整個(gè)配電網(wǎng)的可靠性水平效果更為顯著，驗(yàn)證了2.2.3節(jié)薄弱環(huán)節(jié)分析方法的正確性和有效性。

在實(shí)際配電網(wǎng)中，可利用上述方法識(shí)別系統(tǒng)的薄弱環(huán)節(jié)，并采取相應(yīng)改善措施以提高其可靠性。

4 結(jié)論

本文提出一種配電網(wǎng)可靠性評(píng)估貝葉斯網(wǎng)絡(luò)區(qū)間算法。該算法以各元件、分塊、負(fù)荷點(diǎn)及系統(tǒng)的狀態(tài)為隨機(jī)變量，分析它們之間的關(guān)聯(lián)關(guān)系和影響關(guān)系，建立配電網(wǎng)可靠性評(píng)估的4層貝葉斯網(wǎng)絡(luò)。

為充分計(jì)及元件可靠性參數(shù)的不確定性，引入?yún)^(qū)間數(shù)運(yùn)算法則，基于貝葉斯網(wǎng)絡(luò)的前向推理模式計(jì)算負(fù)荷點(diǎn)和系統(tǒng)的區(qū)間可靠性指標(biāo)。

基于貝葉斯網(wǎng)絡(luò)的后向推理模式，提出通過元件的區(qū)間可靠性參數(shù)和負(fù)荷點(diǎn)節(jié)點(diǎn)類型辨識(shí)負(fù)荷點(diǎn)的薄弱環(huán)節(jié)的方法，及通過各元件分?jǐn)偟膮^(qū)間ENS辨識(shí)系統(tǒng)的薄弱環(huán)節(jié)的方法。對(duì)實(shí)際配電網(wǎng)，可根據(jù)前述方法辨識(shí)配電網(wǎng)薄弱環(huán)節(jié)，并采取相應(yīng)改善措施以提高系統(tǒng)可靠性。

通過對(duì)IEEE RBTS母線2等算例的分析，驗(yàn)證了本文貝葉斯網(wǎng)絡(luò)可靠性評(píng)估區(qū)間算法及薄弱環(huán)節(jié)辨識(shí)技術(shù)的正確性和有效性。

［1］ 陳文高.配電系統(tǒng)可靠性實(shí)用基礎(chǔ)［M］.北京:中國電力出版社，1998:117-124.

［2］ 李志民，李衛(wèi)星，劉迎春.輻射狀配電系統(tǒng)可靠性評(píng)估的故障遍歷算法［J］.電力系統(tǒng)自動(dòng)化，2002，(2):53-56.

［3］ 戴雯霞，吳捷.基于最小路的配電網(wǎng)可靠性快速評(píng)估法［J］.電力自動(dòng)化設(shè)備，2002，(7):29-31.

［4］ 楊文宇，余健明.基于最小割集的配電系統(tǒng)可靠性評(píng)估算法［J］. 西安理工大學(xué)學(xué)報(bào)，2001，(4):387-391.

［5］ 謝開貴，王小波.計(jì)及開關(guān)故障的復(fù)雜配電系統(tǒng)可靠性評(píng)估［J］. 電網(wǎng)技術(shù)，2008，(19):16-21.

［6］ 任震，萬官泉，黃雯瑩.參數(shù)不確定的配電系統(tǒng)可靠性區(qū)間評(píng)估［J］.中國電機(jī)工程學(xué)報(bào)，2003，(12):68-73.

［7］ 鄭彥芹，婁北，張鴻雁，等.基于區(qū)間數(shù)和價(jià)值理論的配電網(wǎng)可靠性評(píng)估［J］.電力系統(tǒng)保護(hù)與控制，2008，(18):19-23.

［8］ 張鵬，王守相.大規(guī)模配電系統(tǒng)可靠性評(píng)估的區(qū)間算法［J］.中國電機(jī)工程學(xué)報(bào)，2004，(3):77-84.

［9］ 霍利民，朱永利，張立國，等.用于電力系統(tǒng)可靠性評(píng)估的貝葉斯網(wǎng)絡(luò)時(shí)序模擬推理算法［J］.電工技術(shù)學(xué)報(bào)，2008，(6):89-95.

［10］ Yu D C，Nguyen T C，Haddawy P.Bayesian Network Model for Reliability Assessment of Power Systems［J］.IEEE Trans on Power Systems，1999，(2):426-432.

［11］ 盧錦玲，楊曉東，粟然，等.一種基于貝葉斯網(wǎng)絡(luò)的配電網(wǎng)可靠性評(píng)估方法［J］.華北電力大學(xué)學(xué)報(bào)，2004，(2):16-19.

［12］ 劉禹良，謝開貴，曹侃.基于貝葉斯網(wǎng)絡(luò)和分塊技術(shù)的配電網(wǎng)可靠性評(píng)估［J］.華東電力，2010，(5):90-95.

［13］ 王清印，呂瑞華.區(qū)間數(shù)的標(biāo)準(zhǔn)表示及其四則運(yùn)算法則與泛灰數(shù)的內(nèi)在聯(lián)系［J］.數(shù)學(xué)的實(shí)踐與認(rèn)識(shí)，2005，(6):216-222.

［14］ Alefeld G，Mayer.Interval Analysis:Theory and Applications［J］.Journal of Computational and Applied Mathematics，2000，(1-2):421-464.

［15］ 李志林.區(qū)間數(shù)的一種改進(jìn)的排序方法［J］.數(shù)學(xué)的實(shí)踐與認(rèn)識(shí)，2004，(6):124-127.

［16］ 張鵬.區(qū)間分析在電力系統(tǒng)可靠性評(píng)估中的應(yīng)用研究［D］.北京:清華大學(xué)，2002.

［17］ 周志仁，周俊建，張志海.區(qū)間數(shù)逆運(yùn)算和簡(jiǎn)單的區(qū)間數(shù)方程［J］.河北煤炭建筑工程學(xué)院學(xué)報(bào)，1994，(2):47-51.

［18］ R.N.Allan，R.Billinton，I.Sjarief，et al.A Reliability Test System for Educational Purposes:Basic Distribution System Data and Result［J］.IEEE Trans on Power Systems，1991，(2):813-820.