程超
哈爾濱電機(jī)廠有限責(zé)任公司 黑龍江哈爾濱 150000
水電是電網(wǎng)中具有靈活調(diào)節(jié)能力的可再生能源,可以有效補(bǔ)償風(fēng),光等間歇性電源對(duì)電網(wǎng)穩(wěn)定性的有害影響,并保持電網(wǎng)的工頻穩(wěn)定。至于水輪發(fā)電機(jī)組,其運(yùn)行模式可能會(huì)發(fā)生變化,可以根據(jù)需要在隔離電網(wǎng)模式下運(yùn)行或集成到大型電網(wǎng)中。在將一個(gè)單元集成到大型電網(wǎng)中的運(yùn)行模式下,單元側(cè)的穩(wěn)定性對(duì)電網(wǎng)影響很小;在隔離電網(wǎng)模式下,單元的調(diào)節(jié)性能直接決定著電網(wǎng)的穩(wěn)定性和電能質(zhì)量。因此,研究與水電機(jī)組隔網(wǎng)布置有關(guān)的穩(wěn)定性和動(dòng)力特性尤為重要。
孤網(wǎng)模式下的水輪機(jī)調(diào)節(jié)系統(tǒng)主要包含調(diào)速器、水輪機(jī)及引水系統(tǒng)、發(fā)電機(jī)和負(fù)載等。機(jī)組在孤網(wǎng)模式下控制器主要采用PID控制,隨動(dòng)系統(tǒng)主要包含主接力器和延遲環(huán)節(jié)??刂破骱碗S動(dòng)系統(tǒng)的數(shù)學(xué)模型分別為:u 為控制器輸出;e 為跟蹤誤差;KP、KI、KD 分別為控制器比例、積分和微分增益;Ty 為接力器反應(yīng)時(shí)間常數(shù),s;s 為拉普拉斯變量;Td 為隨動(dòng)系統(tǒng)的延遲時(shí)間,s;y 為導(dǎo)葉開度;xr、x 分別為機(jī)組轉(zhuǎn)速的給定值、測(cè)量值;bp 為永態(tài)轉(zhuǎn)差系數(shù);yr 為導(dǎo)葉開度的給定值[1]。
水輪機(jī)及引水系統(tǒng)模型。水輪機(jī)及引水系統(tǒng)模型包含水輪機(jī)模型和引水系統(tǒng)模型兩個(gè)部分。本文采用線性水輪機(jī)模型,利用模型綜合特性曲線計(jì)算某一工況點(diǎn)的傳遞系數(shù)來表征水輪機(jī)在該工況點(diǎn)處的流量和力矩特性:q為水輪機(jī)流量;mt 為水輪機(jī)力矩;h 為水輪機(jī)水頭;eqx、eqy、eqh 分別為流量對(duì)轉(zhuǎn)速、開度和水頭的傳遞系數(shù);eqy、eyx、eqh 分別為力矩對(duì)轉(zhuǎn)速、開度和水頭的傳遞系數(shù)。
引水模型包含剛性水擊模型、彈性水擊模型及特征線模型。當(dāng)引水管道長(zhǎng)度較短時(shí),一般采用剛性水擊模型即可滿足工程要求發(fā)電機(jī)和負(fù)載模型。在孤網(wǎng)模式下,發(fā)電機(jī)模型通常采用不計(jì)電磁暫態(tài)的一階動(dòng)力學(xué)模型。該模型通過發(fā)電機(jī)轉(zhuǎn)子運(yùn)動(dòng)方程導(dǎo)出。相應(yīng)的負(fù)載可用常數(shù)表示,由此得到發(fā)電機(jī)的傳遞函數(shù)??傻玫焦戮W(wǎng)模式下水輪機(jī)調(diào)節(jié)系統(tǒng)的整體數(shù)學(xué)模型。
水輪機(jī)調(diào)節(jié)系統(tǒng)的動(dòng)態(tài)行為一般可表示為代數(shù)微分方程形式:
z為系統(tǒng)的狀態(tài)變量,如發(fā)電機(jī)轉(zhuǎn)速;a 為系統(tǒng)的代數(shù)變量,如發(fā)電機(jī)機(jī)端電壓;r 為系統(tǒng)的輸入或外部擾動(dòng),如頻率給定;p為系統(tǒng)的參數(shù),如水流慣性時(shí)間常數(shù)。軌跡靈敏度是變量軌跡對(duì)系統(tǒng)參數(shù)的導(dǎo)數(shù),反映了系統(tǒng)參數(shù)微小變化時(shí)對(duì)變量動(dòng)態(tài)軌跡變化的影響程度,因此在小擾動(dòng)條件下將狀態(tài)變量或代數(shù)變量的解軌跡在參數(shù)p 的當(dāng)前值處進(jìn)行泰勒展開,并忽略二階及以上高階微量,可得到參數(shù)對(duì)變量的軌跡靈敏度指標(biāo):t 為仿真的當(dāng)期時(shí)刻,s;Δp 為參數(shù)的攝動(dòng)量。為便于比較和分析,采用相對(duì)軌跡靈敏度,并取其絕對(duì)平均值作為指標(biāo):分別為初始狀態(tài)和參數(shù);N 為仿真的數(shù)據(jù)長(zhǎng)度;Sp 為相對(duì)軌跡靈敏度;Mmax為將系統(tǒng)某一變量在過渡過程中相對(duì)變化量的絕對(duì)平均值作為因變量,將某一參數(shù)值作為自變量時(shí),所構(gòu)成的函數(shù)在該參數(shù)值處的斜率[2]。
為便于計(jì)算敏感性指標(biāo),仿真的條件設(shè)計(jì)為水電機(jī)組受到負(fù)載擾動(dòng)作用,負(fù)載擾動(dòng)設(shè)置為一個(gè)脈寬為0.1s、方向向下的脈沖信號(hào),此時(shí)可近似認(rèn)為機(jī)組所帶的負(fù)荷瞬間減少至某一定值,且在極短的時(shí)間后負(fù)荷恢復(fù)至初始狀態(tài)。在該條件下,系統(tǒng)狀態(tài)將在控制器的調(diào)節(jié)作用下,從初始平衡點(diǎn)歷經(jīng)短時(shí)的波動(dòng)并再次回到原始平衡點(diǎn)。負(fù)荷擾動(dòng)的形式之所以未選擇階躍信號(hào),一方面是為了便于敏感性指標(biāo)的計(jì)算,根據(jù)軌跡靈敏度計(jì)算公式,僅需計(jì)算系統(tǒng)變量動(dòng)態(tài)變化部分的面積即可;另一方面,在階躍擾動(dòng)后,系統(tǒng)的平衡點(diǎn)會(huì)發(fā)生變化,意味著參數(shù)敏感性既受動(dòng)態(tài)的影響,也受穩(wěn)態(tài)的影響,而實(shí)際上研究關(guān)注的是參數(shù)變化對(duì)系統(tǒng)動(dòng)態(tài)產(chǎn)生的影響。通過比較參數(shù)攝動(dòng)前后系統(tǒng)變量變化的差值即可得到系統(tǒng)變量對(duì)參數(shù)的敏感程度[3]。
水輪機(jī)調(diào)節(jié)系統(tǒng)基本參數(shù),根據(jù)電站運(yùn)行參數(shù)設(shè)置,該模式下的永態(tài)轉(zhuǎn)差系數(shù)bp=0.01。待分析的系統(tǒng)參數(shù)共有14 個(gè),待分析的系統(tǒng)變量共有6 個(gè),分別為控制器輸出u、隨動(dòng)系統(tǒng)輸出y、引水系統(tǒng)輸出h、水輪機(jī)輸出q、mt 及發(fā)電機(jī)輸出x。參數(shù)敏感性分析流程,對(duì)各參數(shù)在其初始值的基礎(chǔ)上分別施加±5%的攝動(dòng),并進(jìn)行負(fù)載擾動(dòng)的仿真試驗(yàn)。仿真總時(shí)長(zhǎng)為50s,其中負(fù)載擾動(dòng)的開始時(shí)刻為1s,持續(xù)時(shí)間為0.1s。由于參數(shù)較多,僅以對(duì)系統(tǒng)狀態(tài)影響最大的參數(shù)(即Ta)為例,分析參數(shù)攝動(dòng)對(duì)系統(tǒng)變量的影響。
綜上所述,目前,水電機(jī)組的穩(wěn)定性和動(dòng)力特性研究主要包括穩(wěn)定性分析,非線性動(dòng)力分析和參數(shù)敏感性分析。其中,參數(shù)靈敏度分析是系統(tǒng)穩(wěn)定性研究的重要內(nèi)容,參數(shù)靈敏度分析結(jié)果對(duì)系統(tǒng)辨識(shí),穩(wěn)定性分析和最優(yōu)控制具有重要意義。在分析水輪機(jī)調(diào)節(jié)系統(tǒng)的參數(shù)敏感性方面,當(dāng)前的研究主要集中在單邊參數(shù)(例如流道參數(shù))對(duì)單邊系統(tǒng)(例如單元振動(dòng))的影響。這些研究集中在一些主要參數(shù)對(duì)特定變量的影響機(jī)理上,但是他們?nèi)狈?duì)主要和次要作用以及系統(tǒng)參數(shù)之間關(guān)系的系統(tǒng)研究。