導航：安全管理網>> 安全論文>> 電力安全論文>>正文

基于均勻設計的供熱機組安全運行區的計算方法

文檔作者：王桐田亮文檔來源：華北電力大學控制與計算機工程學院		點擊數：	更新時間： 2021年06月06日
下載地址：點擊這里	文件大小： 375.32 KB 共7頁	文檔格式： PDF	下載點數： 1 點（VIP免費）

全屏查看

內容預覽 [文件共7頁]

本文件共7頁，如需編輯使用，請下載。

注：預覽效果可能會出現部分文字亂碼（如口口口）、內容顯示不全等問題，下載是正常的。

	文件大小：375.32 KB 共7頁文件格式：PDF 下載點數：1 點（VIP會員免費）

下一篇：基于機組運行安全的水布埡水電站機組負荷…

上一篇：計及機組不確定損耗成本的魯棒安全經濟調度

文本預覽

僅提取頁面文字內容，供快速閱讀使用。

第43卷第5期
2016年9月
華北電力大學學報
Journal of North China Electric Power University
Vo1．43．No．5
Sep．，2016
doi：10．3969／j．ISSN．1007—2691．2016．05．14
基于均勻設計的供熱機組安全運行區的計算方法
王 桐， 田 亮
(華北電力大學控制與計算機工程學院，河北保定071003)
摘要：在線計算汽輪機安全運行區對供熱機組參與電網“調峰調頻” 至關重要。傳統熱力性質工況圖以供熱
抽汽流量、汽輪機低壓缸排汽流量確定汽輪機工作點，但現場缺少這兩個信號的測點或測量不準確。分析機
組發電功率、主蒸汽流量、供熱抽汽壓力、供熱抽汽流量、汽輪機低壓缸排汽流量之間的約束關系，建立發
電功率、主蒸汽流量、供熱抽汽壓力對供熱抽汽流量、汽輪機低壓缸排汽流量的多元擬合模型， 并結合均勻
設計確定決定機組安全運行區的典型工況點。誤差分析表明，該方法擬合精度接近熱平衡方法。
關鍵詞：供熱機組；安全工作區；均勻設計； 多元擬合
中圖分類號：TK272 文獻標識碼：A 文章編號：1007—2691(2016)05—0088—06
Calculation M ethod for Safe Operation Zone of Heat
Supply Units Based on Uniform Design
W ANG Tong，TIAN Liang
(School of Control and Computer Engineering，North China Electric Power University，
Baoding 071003，China)
Abstract： Online calculation of the safe operation zone of turbine plays a significant role in the participation of heat
supply units in the frequency control and peak load regulation in power grid．Working points of turbine is usually deter—
mined by heat extraction steam flow and low-·pressure cylinder exhaust steam flow in traditional working conditions dia--
gram of thermodynamic properties；however，either there is no measuring point of those steam flow or the measuring re—
suhs are inaccurate． By analyzing the constrained relationship of generated power of the units，main steam flow，heat
extraction steam pressure and flow and low-pressure cylinder exhaust steam flow，multiple fitting model of generated
power，main steam flow and heat extraction steam pressure to heat extraction steam flow and low—pressure cylinder ex—
haust steam flow is established．Typical operating points which determine the safe operation zone of the units are found
by combining the uniform design． Error analysis shows that fitting accuracy of this method is similar to heat balance
method．
Key words：heat supply units；safe operation zone；uniform design；multiple fitting
0 引 言
“十一五”期間我國熱電聯產機組的裝機容量
收稿日期：2015—12-22．
基金項目： 國家重點基礎研究發展計劃項目(973 計劃)
(2Ol2cB215200)； 中央高校基本科研業務費專項資
金資助項目(11QG73)．
進一步擴大，供暖季節，熱電聯產機組一般工作在
“以熱定電”方式下，基本不參與調頻調峰電網，但
是，由于新能源的不穩定性，電網側需要供熱機組
承擔一定的調峰調頻任務¨ 。風電的大規模并
網，使電網波動加劇，尤其在供暖季節，風電高發，
“棄風”現象嚴重，更加需要利用供熱機組的儲能，
深度參與電網的調峰調頻，補償風電造成的擾
動 。除發電負荷外，供熱抽汽流量和汽輪機低
第5期 王 桐，等：基于均勻設計的供熱機組安全運行區的計算方法 89
壓缸排汽流量也直接影響供熱機組汽輪機運行狀
態，在線計算供熱機組汽輪機安全運行區間，是實
現這一目標的必要條件 。
現場依據汽輪機廠家提供的汽輪機熱力性質
工況圖判斷其安全運行區 ，但在工況圖中以供
熱抽汽流量和汽輪機低壓缸排汽流量反向確定主
汽流量與發電功率合理范圍，供熱抽汽流量與汽
輪機低壓缸排汽流量測量困難、精度低，造成這一
方法應用困難。
在線計算供熱機組運行安全區的研究主要集
中于機理建模方面：文獻[7]建立了供熱機組的發
電負荷與抽汽流量之間的非線性動態模型，文獻
[8—9]分別用熱平衡法和循環函數法迭代計算出
供熱機組的熱、電負荷關系以及不同抽汽流量下
電負荷的范圍，文獻[10]利用分段弗留蓋爾公式
分析了變工況下供熱機組負荷特性，討論了供熱
機組發電負荷與供熱抽汽壓力、流量等關系。但
都沒有解決需要直接或間接利用供熱抽汽流量或
汽輪機低壓缸排汽流量確定汽輪機安全運行區間
的問題。
采用數據擬合方法，建立供熱抽汽流量與汽
輪機低壓缸排汽流量簡化計算公式，并通過多元
線性擬合的方法確定各流量與發電功率的關系。
在此基礎上，采用均勻設計? 的方法篩選出能夠
決定機組安全運行區的典型工況點，提高擬合的
精確性，推算出供熱機組的安全運行區，為機組深
度參與電網調峰、調頻提供支持。
1 供熱機組特性分析
1．1 熱力系統結構
本文的研究對象為抽汽式供熱機組，汽輪
機與熱網部分如圖1所示。與純凝式機組相
較，供熱機組在中壓缸與低壓缸之間增加了
LV，汽輪機中壓缸的排汽分為兩部分，一部分
經LV進入低壓缸繼續做功，另一部分經過EV
進人熱網加熱器提供熱負荷，換熱后返回熱網
加熱器。熱網加熱器分為管側和殼側兩部分，
熱網循環水在管側流動，通過熱交換吸收殼側
內供熱抽汽的熱量，并在二級換熱站釋放熱量；
供熱抽汽在殼側內釋放能量并凝結成當前壓力
下的飽和水進入除氧器，加入凝結水系統進行
新一輪的鍋爐水循環。
再
1一汽輪機高壓缸；2一汽輪機中壓缸；3一汽輪機低壓缸；4一
LV(汽輪機供熱調節蝶閥)；5一EV(熱網供熱調節蝶閥)；6一
熱網加熱器；7一熱網疏水泵；8一除氧器；9一給水泵
圖1 抽汽式汽輪機熱力系統簡圖
Fig．1 Thermal system diagram of extraction turbine
1．2 運行約束條件
供熱機組在供暖季工作在“以熱定電”方式
下，同時提供熱、電兩種負荷，因此，供熱機組與純
凝機組的運行范圍不盡相同。純凝式機組只受最
大最小負荷限制¨ ，理論上而言，最大負荷的約束
條件為汽輪機四閥全開工況下(VWO)最大進汽
量，但因環境溫度等客觀條件的影響，實際的約束
量為汽輪機最大連續工況(TNL)下的最大進汽
量；最小負荷的約束條件為鍋爐穩燃且汽水均勻
條件下的最小蒸發量。
供熱機組因受供熱負荷的限制，其運行工況
相較于純凝機組，約束條件更多：最大負荷的約束
條件還包括汽輪機額定熱耗考核工況(THA)下，
鍋爐負荷達到最大連續出力狀態(BMCR)時的飽
和蒸汽流量；最小負荷的約束條件還包括汽輪機
低壓缸的最小進汽量。同時，中壓缸排汽壓力過
高會出現中壓缸末級“鼓風”現象；低壓缸進氣壓
力過低會出現低壓缸進汽量不足的現象，汽輪機
反動，振動加劇、軸向推力變化，即所謂“悶缸”¨
現象，嚴重時會危及汽輪機運行安全，導致汽輪機
“跳閘”停機事故。
1．3 熱力系統工況圖
供熱機組熱力系統工況圖是計算安全運行區
的主要手段之一，其簡圖如圖2所示。圖中，虛線
為非供熱季機組的安全運行區域，實線為供熱季
機組的安全運行區域，其圍成的面積即為供熱工
況下機組安全運行區，從圖中可以明顯看出，供暖
季機組的運行區間明顯減小。
90 華北電力大學學報 2016矩
圖2 簡化機組工況圖
Fig．2 Simple operating mode chart of turbine
從上圖中可以看出決定機組運行范圍的參量
主要為：主汽流量、供熱壓力、供熱流量與汽輪機
低壓缸排汽流量。當機組在最小電力負荷時，低
壓缸排汽流量為最小冷卻蒸汽流量，據此可以推
算出機組的最小出力與主蒸汽、供熱抽汽流量的
關系；在最大負荷時，主蒸汽流量為TMCR工況下
的主蒸汽流量，根據額定工況下的各段流量焓值
可以推算出最大出力與供熱抽汽流量、低壓缸排
氣流量之間的關系。
圖3為某供熱汽輪機額定抽汽壓力下的工況
圖，在測得機組供熱抽汽流量與低壓缸排汽流量
之后，即可圈出機組的安全運行區域。但這種方
法在應用上有兩個難點：
(1)現場沒有供熱抽汽流量與低壓缸排汽流
量的測點，由于EV、LV的非線性，難以根據閥的
開度來判斷當前流量。
(2)供熱壓力變化時，工況線發生變化，需引
入壓力修正，因此，在圈算安全區時需要找到對應
的工況圖，難以編程，實現在線計算。
因此，構造供熱抽汽流量與汽輪機低壓缸排
GEX／t·h-I：供熱抽汽量
l GE 00
2GE _300
3 G 00
4 G臣 100
5 GE
GK／t·h-ll低壓缸排汽量
llGI 0
22GJ O0
33 GK=400
44 GK-500
0 4O 80120160 200 240 280 320 360400 發電功率／MW
圖3 供熱汽輪機熱力系統工況圖
Fig．3 Operating mode chart of thermal system
汽流量的簡化計算公式并引入供熱壓力修正是計
算安全運行區的必要前提。
2 均勻設計
2．1 均勻設計的基本理論
設 “’∈R 表示m個因素(待考察變量)的
取值向量，m表示因素的個數，因素在取值范圍內
的標幺值 E R：，n稱之為水平，試驗結果稱之
為響應n引。均勻設計的目的就是建立一個易于表
達與計算的近似模型式(1)來逼近實際模型。
Y = “’
，
’
， ?
， )=Σ盧 “ )(1)
式中：JB 為插值系數。
若存在一個最大可容忍誤差 ，使得
I Y—f( ¨’，?， )l≤6 (2)
則為通過閾值檢驗的模型。
記試驗響應Y在R 的平均值E(Y)為
E(Y)=』 u ，?， )dx¨ ，?，dx ’(3)
若在試驗區域內做了It次試驗，其試驗點集
為D = { =( ¨，? ¨， ，k=1，? ，n}，并
用響應Y在這n個點的均值 來估計實際模型，
則這個估計的偏差由kolsma．Hlawaka不等式¨ 給
出：
E(y)一Y。≤ ( D(D ) (4)
其中：
一 1
Y。 n =÷n Σ，( “ )
上式中：
( 是函數，在R 上的全變差，D(D )為點
集D 在實驗區域R 上的星偏差(star disctepancy)，
它是度量D 在 上散布均勻的測量度。
由不等式(4)可知：
(1)D(D )偏差越小，則估計的誤差越小，試
驗點集均勻分布時偏差達D(D )到極小。
(2)kolsma—Hlawaka不等式右邊對函數的依
賴只存在于v(f)一項，若現有另外一個模型Y=
g( ‘¨，? ， ‘ ’)+8，且V(g)= ( 。如果利用
同一組試驗點集D 來估計，與g，這時此時兩者有
著相同的誤差上界：
E(g)一言(D )≤V(g)D(O )= ∽D(D )(5)
均勻設計對模型的變化有穩健性，它不僅能
第5期 王桐，等：基于均勻設計的供熱機組安全運行區的計算方法 91
夠對某個特定的，( ‘¨，? ， )+ 做出估計，并
且可以對任意多個 ( 有界的g( ‘¨，? ， ‘“ )+
占做出估計。
2．2 均勻性的度量
現定義：在一個／'t的試次驗中，令Ⅳ(D ，[0，
])為D 中的點落在落在[0， ]中的個數，當D
中的點在R 中均勻分布時，B(D ，[0， ])應與
[0， ]的體積相似 ” 。
記兩者的偏差為D( )，則
D( )：l ． 二 一v。l([o， ])l(6)
記兩者的 一星偏差為D(A)，則
。(A)=(【 l 一V。·([。， ])I 】Ti)
(7)
式(6)定義的偏差衡量均勻性不夠靈敏，在均
勻設計中，通常采用數論中的偽蒙特卡洛方法的
p一星偏差來描述均勻度。由于 p一星偏差忽略
了低于s維投影的偏差，對體現點集均勻性十分重
要的一維和二維沒有體現，有時會出現不合理的
結果，因此，Hickerenell將 p一星偏差的修正為中
心化 p一星偏差，記做CD ：
c。 c。 ，=c( ) 一 毫 c 2+ —
+ Σ11(1+ + 一[ )]寺
(8)
其中：
一丟】
CD2滿足kolsma—Hlawaka不等式，并且將D
中點的次序改變，或將每個試驗點的m個坐標任
意調換，相應的偏差不變，適合作為衡量一個設計
均勻度的指標，因此本文將使用CD 的大小作為
衡量均勻度的標準。
3 計算實例
3．1 機組簡介
抽汽式機組與純凝機組的主要區別在于供熱
部分，本文以JQ電廠單抽凝汽式機組為例計算其
供熱安全區，其熱力特性如圖4所示。供熱抽汽
管道布置在中壓缸與低壓缸之間，中壓缸排汽分
為兩路，一路通過蝶閥進入低壓缸，另一路則進入
供熱抽汽管道，經熱網換熱后由熱網疏水泵輸送
進入除氧器。其額定供熱壓力為0．35 MPa，額定
供熱抽汽流量為400 t／h。
1一再熱蒸汽(出)；2一再熱蒸汽(人)；3一供熱抽汽；
圖4 單抽凝汽機組的熱力特性
Fig．4 Thermed odynamic characteristics of extraction tur—
bine
3．2 均勻設計與工況點選擇
模型因素數的選擇與工況點數量密切相關，
最少工況點與因素數、待擬合的模型有關，關系如
表1所示。
表1 最少工況點與模型類型關系
Tab．1 Relationship of minimum operating points and mode
types
理論上，水平數較少時，無法擬合高階多項式
模型，故對模型的變化缺乏“穩健性”，水平數較多
就可以有更多的模型選擇，但實際應用中，由于歷
史運行數據庫中的數據的隨機性較差，分布不均
勻等問題，隨著水平數的增加，模型的均勻程度先
升高后下降¨ ，如圖5所示。
圖5 不同因素數下模型的￡p一星偏差
Fig．5 Lp— star deviation in different models of different
factors
92 華北電力大學學報
在選擇模型的因素與水平時，可按照表1所
示的最少工況點，從線性模型開始計算，如若計算
結果滿足精度，則停止計算；若精度無法滿足需
求，則計算高一級別的模型，以此類推，直至模型
精度滿足要求為止。
現已知發電功率Q范圍與4個變量(主蒸汽
流量q 、供熱蒸汽流量q 低壓缸排氣流量q 供
熱壓力P )有關，故選取因素數m =4，最少典型
工況點的數量為5組；供熱抽汽流量的取值與機
故選取因素數m，=3，由機理分析可知關系為線性
或者二次關系，最少典型工況點的數量為4組；汽
輪機低壓缸排汽壓力的取值與機組當前的發電功
率、主蒸汽流量、供熱壓力有關，故選取因素數m，
=3，由機理分析可知關系為線性或者二次關系，
最少典型工況點的數量為4。
均勻設計的數據選擇可以借助均勻設計表或
采用好格子法、拉丁方法計算。本文中借助文獻
[17]出列出的均勻設計表選擇數據，所使用的均
組當前的發電功率、主蒸汽流量、供熱壓力有關， 勻數據表如表2所示。
襲2 基于均勻設計的典型工況點數據
Tab．2 Typical operating point data of uniform design experimentation
供熱抽汽流量擬合模型為
q =一122．45p +2．46q 一6．62Q一194．03 (9)
汽輪機低壓缸排汽流量擬合模型為
q。=105．75p 一1．Olq +4．43Q+78．04 (10)
式(9)～(10)中：g 為供熱蒸汽流量，t／h；g 為
低壓缸排汽流量，t／h；p 為供熱抽汽壓力，MPa；q
為當前工況下的主蒸汽流量，t／h；Q為當前機組發
電功率，Mw。
得到機組的供熱抽汽流量與低壓缸排汽流量
之后即可圈出安全運行區，但是圈出的安全運行
區難以編程實現在線計算，因此，需要擬合出發電
功率模型。
可調度發電功率擬合模型為
Q =一26．5lp +0．25q 一
0．047q +0．11qp+9．43 (11)
得到當前機組的發電功率模型后，即可通過
最大最小條件的約束獲得當前機組的安全運行
區，具體為：
最小發電功率：低壓缸排汽流量為最小冷卻
蒸汽流量。
最大發電功率：主蒸汽流量為汽輪機最大連
續出力(TMCR)工況下的主蒸汽流量。
將上述約束條件帶入上述模型即可得到機組
運行區間，表3為8組不同工況下的測量數據與
計算數據，其中抽汽流量、汽輪機低壓缸排汽流
量、發電功率區間均為本文模型計算得出，工況圖
區間為機組相應工況圖中讀出。從表中可以看
出，模型的計算結果與工況圖十分接近，所需數據
較少計算簡單，易于在機組DCS(分散控制系統)
中組態實現。
表3 模型計算值
Tab．3 Caculation results of models
第5期 王 桐，等：基于均勻設計的供熱機組安全運行區的計算方法 93
若不考慮尋找具有決定性的典型工況點數
據，而是采用最d~-<-乘法擬合全部的工況點，得到
的擬合模型為
供熱抽汽流量擬合模型為
g =一169．93p +1．71q 一5．10Q +83．62(12)
汽輪機低壓缸排汽流量擬合模型為
q。 =195．67p 一0．83q +4．12Q 一46．95(13)
發電功率擬合模型為
Q ：一46．14p +0．219 3q +
0．198 3q +0．001 5qP +8．83 (14)
式(12)一(14)中：g 為供熱抽汽流量，t／h；p，
為供熱抽汽壓力，MPa；g 為主蒸汽流量，t／h；Q 為
發電功率，MW。
3．3 模型的評價標準
本文采用交叉檢驗的方式評價模型，相應的
指標為關聯度、絕對差值率、均方根誤差率。
關聯度：
R ： (15)
√Σ(y —y) ( 一 )
絕對差值率
= Σ 1 Y 一y l (16)
均方根誤差率
— — — — — 一
R眥 √亡 (y —Yi)2 (17)
式中：，， 為實際的功率輸出，Y 為模型計算的功率
輸出，Y 為實際功率輸出的平均值，Y 為模型計算
功率輸出的平均值。
其數學意義為：關聯度表示擬合模型與真實
對象特性之間的關聯函數，關聯度越大，模型的擬
合效果越好；絕對差值率用于評價模型的差異度，
體現著樣本個體的差異對模型的影響，其值越小，
模型的預測能力越好；均方根誤差率反映著樣本
整體的差異對模型的影響，其值越小，模型越精
確。
由表4可以看出，基于均勻設計的模型的相
關度較高，已接近機理計算，另外，均勻設計得到
的模型的精確度要優于最小二乘擬合。
表4 模型的評價指標
Tab．4 Evaluation indexes of models
本文模型 最dx-乘擬合模型
q qp Q gr qp Q
／t．h一 ／t．h一 ／MW ／t·h一 ／t．h一 ／MW
R 0．973 0．968 0．955 0．955 0．942 0．938
M^E／％ 2．13 2．75 3．83 5．12 4．14 7．35
RMSE／％ 3．11 3．37 4．92 5．73 4．72 8．84
4 結 論
本文研究了一種基于均勻設計的供熱機組運
行安全區的在線計算方法，得到如下結論：
(1)建立了供熱機組的主汽流量、發電負荷、
供熱抽汽壓力的擬合模型。
(2)機組運行工況眾多，找出能決定機組安全
區的典型工況點，是提高擬合模型精度的關鍵。
通過均勻設計的方法尋找典型工況點，得到的擬
合模型的精度滿足要求，且數據量少，計算過程簡
單，易于實現在線計算。
(3)通過JQ電廠300 MW 機組的運行數，驗
證了算法的可行性，模型的計算精度與工況圖相
當。
參考文獻：
[1]徐則林，姜燕妮，徐磊，等．熱電聯產機組熱能深
度梯級利用的研究[J]． 華北電力大學學報，
2014。41(5)：107一l12．
[2]付殿崢，黃國和．基于因子交互分析的熱電廠經濟
性優化研究[J]．華北電力大學學報，2014，41
(4)：101—106．
[3]劉鑫屏，田亮，王琪．補償風電擾動的供熱機組快
速變負荷控制方法[J]． 電力系統自動化，2014，
38 (6)：26—32．
[4]鄧拓宇，田亮，劉吉臻．利用熱網儲能提高供熱機
組調頻調峰能力的控制方法[J]． 中國電機工程學
報，2015，35 (14)：3626—3633．
[5]畢慶生，呂項羽，李德鑫，等．基于熱網及建筑物
蓄熱特性的大型供熱機組深度調峰能力研究[J]．
汽輪機技術，2014，56 (2)：141—144．
[6]王宏偉，馮矗飛，馮林魁，等．供熱機組調峰性能
試驗[J]．熱力發電，2014，43(4)：106—109．
[7]王漪，薛永鋒，張敏，等．基于能量平衡法的供熱
機組電量分析數學模型[J]． 電力系統自動化，
2014，38 (8)：108一l12．
(下轉第104頁)
104 華北電力大學學報 2016年
(上接第93頁)
[8]李建剛，李麗萍，楊小琨．供熱汽輪機新汽耦合做功
計算的研究[J]．汽輪機技術，2010，52(6)：431—434
+442．
[9]呂瀝鋒，王培紅，金旭英，等．供熱機組熱化發電率算
法研究[J]．汽輪機技術，2002，44(4)：210—212．
[10]吳龍，袁奇，劉昕．供熱機組熱電負荷最佳分配方法分
析[J]．中國電機工程學報，2012，32(35)：6—12．
[11]黃永青，梁昌勇，張祥德．基于均勻設計的蟻群算法參
數設定[J]．控制與決策，2006，21(1)：93—96．
[12]林振嫻，楊勇平，何堅忍．熱網加熱器在熱電聯產系統
中的全工況分析[J]．中國電機工程學報，2010，30
(23)：l4—18．
[13]王琪．風電規模化并網條件下供熱機組優化控制研
究[D]．北京：華北電力大學，2013．
[14]方開泰，劉民千，周永道．試驗設計與建模[M]．北京：
高等教育出版社，2011．
[15]方開泰，馬長興．正交與均勻試驗設計[M]．北京：科
學出版社，2001．
[16]王學棟，魏東，孫書耀，等．不同類型供熱機組的電熱
負荷優化分配和調峰性能[J]．汽輪機技術，2010，52
(5)：327—390．
[17]趙曙光，焦李成，王宇平，等．基于均勻設計的多目標
自適應遺傳算法及應用[J]．電子學報，2004，32
(1O)：1723—1725+1729．
作者簡介：王桐(1991一)，女，碩士研究生，研究方向為模式識
別與數據融合。

網友評論 more

電力安全論文最新內容

05-01