知財求人 - 知財ポータルサイト「IP Force」
特表2025-502588ＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正する方法
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公表特許公報(A)
(11)【公表番号】
(43)【公表日】2025-01-28
(54)【発明の名称】ＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正する方法
(51)【国際特許分類】
F01D 17/24 20060101AFI20250121BHJP
F01D 17/10 20060101ALI20250121BHJP
【ＦＩ】
F01D17/24 N
F01D17/10 C
F01D17/24 F
F01D17/24 H
【審査請求】有
【予備審査請求】未請求
(21)【出願番号】P 2024519686
(86)(22)【出願日】2023-10-31
(85)【翻訳文提出日】2024-03-28
(86)【国際出願番号】 CN2023128612
(87)【国際公開番号】W WO2024120074
(87)【国際公開日】2024-06-13
(31)【優先権主張番号】202211573195.X
(32)【優先日】2022-12-08
(33)【優先権主張国・地域又は機関】CN
(81)【指定国・地域】
【公序良俗違反の表示】
（特許庁注：以下のものは登録商標）
１．ＭＡＴＬＡＢ
(71)【出願人】
【識別番号】524120376
【氏名又は名称】云南電力試験研究院（集団）有限公司
(74)【代理人】
【識別番号】110002077
【氏名又は名称】園田・小林弁理士法人
(72)【発明者】
【氏名】宋小竜
(72)【発明者】
【氏名】劉友寛
(72)【発明者】
【氏名】伍陽陽
(72)【発明者】
【氏名】文天舒
(72)【発明者】
【氏名】袁華宇
【テーマコード（参考）】
3G071
【Ｆターム（参考）】
3G071BA09
3G071BA33
3G071DA05
3G071EA03
3G071EA06
3G071EA07
3G071FA03
3G071FA05
3G071FA06
3G071HA04
3G071JA04
(57)【要約】
本発明は、ＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正する方法に関し、火力発電所の熱自動制御の技術分野に属する。該システムは、データ収集モジュール、データ処理モジュール、単一バルブ制御モードの各分岐高圧制御バルブのパラメータ補正モジュール、シーケンスバルブ制御モードの各分岐高圧制御バルブのパラメータ補正モジュールから構成されｒ、本発明は、蒸気タービンの制御バルブ非線形補償データを提供し、ＤＥＨ制御線形補正を実施した後、ユニット負荷後蒸気タービンの総合バルブ位置指令変化は、発生した蒸気流量およびユニット負荷に線形に対応して変化するのが保証され、ユニット負荷制御の迅速性および安定性が向上し、ユニット負荷制御の振動の発生を防止し、ユニットＡＧＣと一次周波数変調の制御品質の向上に良好な促進作用があり、火力発電ユニットの安全、安定、経済的、信頼性の高い運転に重要な意義がある。
【選択図】図２
【選択図】図２
【特許請求の範囲】
【請求項1】
データ収集：ユニットをＤＥＨバルブ制御モードに切り替え、ユニットの一次周波数変調制御を切り替え、ユニットのＡＧＣ制御を切り替え、機械と炉の協調制御を手動に切り替え、ユニット燃料自動を切り替え、吸気酸素自動を切り替え、ユニット主蒸気温度自動、給水自動、炉圧自動を入れ、単一高圧制御バルブ１００～０％および０～１００％の単一ステップ５％のステップ外乱を行い、ある制御バルブが動作しているとき、他の制御バルブはすべて開いており、単一高圧制御バルブの動作が完了した後すべての制御バルブを全開に保ち、その後５％のステップでシーケンスバルブ投入の逆順序に従って最終段階で制御バルブを全開に保つまで制御バルブを徐々に閉じ、調整過程全体中のユニット負荷、機械前主蒸気圧力、調圧器圧力、主蒸気温度および調圧器後温度を記録するステップ（１）と、流量特性曲線計算：ステップ（１）で得られたデータを用いて、フリューゲル式に従って各高圧制御バルブ流量特性曲線を計算するステップ（２）と、
単一バルブ制御モードの各分岐高圧制御バルブのパラメータ補正：単一バルブ制御モードで、実際の単一バルブ正規化相対流量を計算し、実際の単一バルブ正規化相対流量をｙ、対応のバルブ開度をｘとしてフィッティングし、高次連続関数ｙ＝ｆ(ｘ)、ＤＣＳで現在ユニット単一バルブで使用されるポリゴン関数を見つけ、高次連続関数ｙ＝ｆ(ｘ)を利用して総合バルブ位置指令またはバルブ開度を補正するステップ（３）と、
シーケンスバルブ制御モードの各分岐高圧制御バルブのパラメータ補正：シーケンスバルブ制御モードで、すべての制御バルブの全開から全閉までａ個のステップがあり、バルブ全開の場合、最終段階全開まで閉じ順序に従ってバルブを順次閉じる過程中、発生した変曲点数がａ－１であり、変曲点に対応する相対流量ｙｍを計算するステップ（４）と、を含み、
ＤＣＳでシーケンスバルブ制御モードで、各支制御バルブで使用されるポリゴン関数を見つけ、該関数において、Ｙｑｉは離散統合されたバルブ位置指令であり、Ｘｑｉは対応のバルブ開度であり、
実際の各高圧制御バルブ正規化相対流量をｙ、対応のバルブ開度をｘとしてフィッティングし、各高圧制御バルブ高次連続関数ｙｉ＝ｆｉ(ｘ)を得、ｉはｉ番目の高圧制御バルブであり、
変曲点に対応する相対流量ｙｍおよび各高圧制御バルブ高次連続関数ｙ＝ｆ(ｘ)、ポリゴン関数に基づいて、総合バルブ位置と相対流量の対照関係を計算して得、該対照関係に基づいて、総合バルブ位置指令またはバルブ開度を補正する、ことを特徴とするＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正する方法。
【請求項2】
ステップ（１）において、記録時、記録サンプリング間隔は１ｓである、ことを特徴とする請求項１に記載のＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正する方法。
【請求項3】
ステップ（２）は、具体的に、ユニットがｉ番目の高圧調整バルブ開度ｊ下で、ｊ∈［０、１００％］、該高圧調整バルブ前後の
Ｐ１ｉｊはｉ番目の高圧制御バルブ開度ｊの時の調圧器圧力であり、単位ＭＰａ、
ＴＳｉｊはｉ番目の高圧制御バルブ開度ｊの時の機械前主蒸気温度であり、単位℃、
Ｔ１ｉｊはｉ番目の高圧制御バルブ開度ｊの時の調圧器後温度であり、単位℃である。
試験範囲内の相対流量を０～１００％正規化し：
ＦＲｉ０はｉ番目の高圧制御バルブのバルブ開度０％時の補正後蒸気圧力比であり、
ＦＲｉ１００はｉ番目の高圧制御バルブのバルブ開度１００％時の補正後蒸気圧力比であり、
得られた相対流量に基づいて、ｉ番目の高圧制御バルブ流量特性曲線を描く、ことを特徴とする請求項１に記載のＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正する方法。
【請求項4】
ステップ（３）において、実際の単一バルブ正規化相対流量を計算することは、具体的に、ＤＣＳで現在ユニット単一バルブで使用されるポリゴン関数を見つけ、該関数において、Ｙｎは離散統合されたバルブ位置指令であり、ＸｎはＹｎに対応するバルブ開度であり、実際の単一バルブ正規化相対流量
ポリゴン関数中のＸｎを上記フィッティング後の高次連続関数に代入し、補正後の総合バルブ位置指令Ｙｎ’を得、または逆関数法ｘ＝ｆ－１(ｙ)により、ポリゴン関数中のＹｎを高次連続関数を解くことに代入し、補正後のバルブ開度Ｘｎ’を得る、ことを特徴とする請求項３に記載のＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正する方法。
【請求項5】
ステップ（４）の方法は具体的に、シーケンスバルブ制御モードで、すべての制御バルブの全開から全閉までａ個のステップがあり、バルブ全開の場合、最終段階全開まで閉じ順序に従ってバルブを順次閉じる過程中、発生した変曲点数がａ－１であり、最初に閉じたバルブに対応する変曲点番号がｍ＝１であり、２番目に閉じたバルブに対応する変曲点番号がｍ＝２であり、以降同様であり、最後の変曲点番号がｍ＝ａ－１であり、
変曲点に対応する相対流量ｙｍを計算し、
ＤＣＳでシーケンスバルブ制御モードで、各支制御バルブで使用されるポリゴン関数を見つけ、該関数において、Ｙｑｉは離散統合されたバルブ位置指令であり、Ｘｑｉは対応のバルブ開度であり、
実際の各高圧制御バルブ正規化相対流量ＦＲｉｊ％をｙ、対応のバルブ開度をｘとしてフィッティングし、各高圧制御バルブ高次連続関数ｙｉ＝ｆｉ(ｘ)を得、ｉはｉ番目の高圧制御バルブであり、
変曲点に対応する相対流量ｙｍおよび各高圧制御バルブ高次連続関数ｙ＝ｆ(ｘ)、ポリゴン関数に基づいて、総合バルブ位置と相対流量の対照関係を計算して得ることは、具体的に、
１、各変曲点の総合バルブ位置指令および相対流量を決定し、すなわち、ポリゴン関数からすべての変曲点の総合バルブ位置指令を見つけ、計算された変曲点に対応する相対流量に対応させ、
２、総合バルブ位置が［０,ａ－１番目の変曲点に対応する総合バルブ位置指令］の範囲内にある場合、該総合バルブ位置をｘとしてこの範囲内で開く必要のある各高圧制御バルブの高次連続関数に代入し、対応のｙｉ値を求め、その後、該総合バルブ位置下でこの範囲内で開く必要のあるすべての高圧制御バルブのｙｉ値の平均値を求め、ａ－１番目の変曲点に対応する相対流量を乗じることで、該総合バルブ位置に対応する実際の相対流量を得、
３、総合バルブ位置が(ａ－１番目の変曲点に対応する総合バルブ位置指令,ａ－２番目の変曲点に対応する総合バルブ位置指令］の範囲内にある場合、該総合バルブ位置をｘとしてこの範囲内で開く必要のある各高圧制御バルブの高次連続関数に代入し、対応のｙｉ値を求め、ここで、全開した高圧制御バルブを除外し、その後、該総合バルブ位置下でこの範囲内で開く必要のあるすべての高圧制御バルブのｙｉ値の平均値を求め、（ａ－２番目の変曲点に対応する相対流量－ａ－１番目の変曲点に対応する相対流量）＋ａ－１番目の変曲点に対応する相対流量を乗じ、該総合バルブ位置に対応する実際の相対流量を得、以降同様であり、
４、総合バルブ位置が(第１変曲点に対応する総合バルブ位置指令,１］の範囲内にある場合、該総合バルブ位置をｘとしてこの範囲内で開く必要のある各高圧制御バルブの高次連続関数に代入し、対応のｙｉ値を求め、ここで、全開した高圧制御バルブを除外し、その後、該総合バルブ位置下でこの範囲内で開く必要のあるすべての高圧制御バルブのｙｉ値の平均値を求め、（第１変曲点に対応する相対流量－第２変曲点に対応する相対流量）＋第２変曲点に対応する相対流量を乗じ、該総合バルブ位置に対応する実際の相対流量を得、
その後、該対照関係に基づいて、総合バルブ位置指令またはバルブ開度を補正する、ことを特徴とする請求項３に記載のＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正する方法。
【請求項6】
対照関係において理論曲線からの偏差が最も大きい点を補正する、ことを特徴とする請求項５に記載のＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正する方法。
【請求項7】
請求項１～４のいずれか１項に記載のＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正する方法を採用する、ＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正するシステムであって、試験過程中、ユニット負荷、機械前主蒸気圧力、調圧器圧力、主蒸気流量、給水流量、主蒸気脱熱水流量、主蒸気温度、再熱蒸気温度、蒸気ドラム圧力および蒸気タービン全エネルギー流量を収集するように構成されたデータ収集モジュールと、
データ収集モジュールによって収集されたデータに基づいて、フリューゲル式を利用して各高圧制御バルブ流量特性曲線を計算するように構成されたデータ処理モジュールと、
実際の単一バルブ正規化相対流量を計算し、実際の単一バルブ正規化相対流量をｙ、対応のバルブ開度をｘとしてフィッティングし、高次連続関数ｙ＝ｆ(ｘ)を得、ＤＣＳで現在ユニット単一バルブで使用されるポリゴン関数を見つけ、高次連続関数ｙ＝ｆ(ｘ)を利用して総合バルブ位置指令またはバルブ開度を補正するように構成された単一バルブ制御モードの各分岐高圧制御バルブのパラメータ補正モジュールと、
変曲点に対応する相対流量ｙｍを計算し、ＤＣＳでシーケンスバルブ制御モードで、各支制御バルブで使用されるポリゴン関数を見つけ、該関数において、Ｙｑｉは離散統合されたバルブ位置指令であり、Ｘｑｉは対応のバルブ開度であり、実際の各高圧制御バルブ正規化相対流量をｙ、対応のバルブ開度をｘとしてフィッティングし、各高圧制御バルブ高次連続関数ｙｉ＝ｆｉ(ｘ)を得、ｉはｉ番目の高圧制御バルブであり、変曲点に対応する相対流量ｙｍおよび各高圧制御バルブ高次連続関数ｙ＝ｆ(ｘ)、ポリゴン関数に基づいて、総合バルブ位置と相対流量の対照関係を計算して得、該対照関係に基づいて、総合バルブ位置指令またはバルブ開度を補正するように構成されたシーケンスバルブ制御モードの各分岐高圧制御バルブのパラメータ補正モジュールと、を備える、ことを特徴とするＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正するシステム。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、火力発電所の熱自動制御の技術分野に関し、具体的に、ＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正する方法に関する。
【背景技術】
【0002】
「カーボンピーク・カーボンニュートラル」の目標の提案により、風力発電と太陽光発電を主体とする新しいタイプの電力システムが構築されつつある。風力発電と太陽光発電の変動性と不確実性の影響を受け、大規模な新エネルギー系統連系は多くの問題をもたらし、特に顕著なのは周波数の問題である。風力発電と太陽光発電の系統連系がもたらす周波数問題は、一方ではそれぞれの出力特性、変動性、不確実性によって引き起こされ、他方では回転慣性の欠如と調整後のサグ調整の問題によって、周波数のさらなる悪化につながる。高品質のＦＭ資源である火力発電ユニットは、基本的に約３０秒でその場に反応することができ、より大きな慣性モーメントを頼りに、より高いＦＭ効果を確保することができる。
【0003】
火力発電設備は、新電力システムの「バラスト」として良い役割を果たすだけでなく、現在の新電力システムは建設に不可欠な部分でもある。大多数のユニットのバルブ特性曲線は、タービンが工場から出荷されるときに策定された特性曲線に拡張されており、各発電企業がタービンのバルブをオーバーホールするユニットのオーバーホールを考慮せず、バルブのストロークが多かれ少なかれ調整され、バルブの流量特性が変化し、特性曲線の実験的補正が行われず、元の特性曲線が使用され続けているため、ＤＣＳで使用されるバルブ流量特性曲線は実際のバルブ流量特性と一致せず、ユニットの統合バルブ位置指令は、発生蒸気流量とユニット負荷に対して大きな非線形性を示し、火力発電ユニットの一次ＦＭの応答速度と調整効果に大きな影響を与える。ＤＥＨバルブの流量特性を再発見し、シングルバルブとシーケンスバルブの制御パラメータを修正することは、火力発電ユニットの一次ＦＭ能力を向上させ、新しい電力システムの安全で安定した運転を確保するために大きな意義がある。従って、火力発電所の熱自動化制御技術の分野では、既存技術の欠点をどのように克服するかが解決すべき課題である。
【発明の概要】
【0004】
ＤＣＳで使用されるバルブ流量特性曲線は実際のバルブ流量特性と一致せず、ユニットの総合バルブ位置指令変化は発生した蒸気流量およびユニット負荷に対して大きな非線形を示し、火力発電ユニットの一次周波数変調応答速度および調整効果に影響を与え、この問題は解決されていない。本発明の目的は、先行技術の欠点を解決するために、ＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正する方法を提供し、試験によりユニットＤＥＨバルブ流量特性を分析し、バルブ流量特性曲線を確立し、単一バルブとシーケンスバルブの制御パラメータを補正することにより、ユニットの総合バルブ位置指令変化は発生した蒸気流量およびユニット負荷に対してほぼ線形関係を示す。
【0005】
上記目的を達成するために、本発明で採用した技術的解決策は以下の通りである。
【0006】
ＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正する方法は、
データ収集：ユニットをＤＥＨバルブ制御モードに切り替え、ユニットの一次周波数変調制御を切り替え、ユニットのＡＧＣ制御を切り替え、機械と炉の協調制御を手動に切り替え、ユニット燃料自動を切り替え、吸気酸素自動を切り替え、ユニット主蒸気温度自動、給水自動、炉圧自動を入れ、単一高圧制御バルブ１００～０％および０～１００％の単一ステップ５％のステップ外乱を行い、ある制御バルブが動作しているとき、他の制御バルブはすべて開いており、単一高圧制御バルブの動作が完了した後すべての制御バルブを全開に保ち、その後５％のステップでシーケンスバルブ投入の逆順序に従って最終段階で制御バルブを全開に保つまで制御バルブを徐々に閉じ、調整過程全体中のユニット負荷、機械前主蒸気圧力、調圧器圧力、主蒸気温度および調圧器後温度を記録するステップ（１）と、
流量特性曲線計算：ステップ（１）で得られたデータを用いて、フリューゲル式に従って各高圧制御バルブ流量特性曲線を計算するステップ（２）と、
単一バルブ制御モードの各分岐高圧制御バルブのパラメータ補正：単一バルブ制御モードで、実際の単一バルブ正規化相対流量を計算し、実際の単一バルブ正規化相対流量をｙ、対応のバルブ開度をｘとしてフィッティングし、高次連続関数ｙ＝ｆ(ｘ)、ＤＣＳで現在ユニット単一バルブで使用されるポリゴン関数を見つけ、高次連続関数ｙ＝ｆ(ｘ)を利用して総合バルブ位置指令またはバルブ開度を補正するステップ（３）と、
シーケンスバルブ制御モードの各分岐高圧制御バルブのパラメータ補正：シーケンスバルブ制御モードで、すべての制御バルブの全開から全閉までａ個のステップがあり、バルブ全開の場合、最終段階全開まで閉じ順序に従ってバルブを順次閉じる過程中、発生した変曲点数がａ－１であり、変曲点に対応する相対流量ｙｍを計算するステップ（４）と、を含み、
ＤＣＳでシーケンスバルブ制御モードで、各支制御バルブで使用されるポリゴン関数を見つけ、該関数において、Ｙｑｉは離散統合されたバルブ位置指令であり、Ｘｑｉは対応のバルブ開度であり、
実際の各高圧制御バルブ正規化相対流量をｙ、対応のバルブ開度をｘとしてフィッティングし、各高圧制御バルブ高次連続関数ｙｉ＝ｆｉ(ｘ)を得、ｉはｉ番目の高圧制御バルブであり、
変曲点に対応する相対流量ｙｍおよび各高圧制御バルブ高次連続関数ｙ＝ｆ(ｘ)、ポリゴン関数に基づいて、総合バルブ位置と相対流量の対照関係を計算して得、該対照関係に基づいて、総合バルブ位置指令またはバルブ開度を補正する。
データ収集：ユニットをＤＥＨバルブ制御モードに切り替え、ユニットの一次周波数変調制御を切り替え、ユニットのＡＧＣ制御を切り替え、機械と炉の協調制御を手動に切り替え、ユニット燃料自動を切り替え、吸気酸素自動を切り替え、ユニット主蒸気温度自動、給水自動、炉圧自動を入れ、単一高圧制御バルブ１００～０％および０～１００％の単一ステップ５％のステップ外乱を行い、ある制御バルブが動作しているとき、他の制御バルブはすべて開いており、単一高圧制御バルブの動作が完了した後すべての制御バルブを全開に保ち、その後５％のステップでシーケンスバルブ投入の逆順序に従って最終段階で制御バルブを全開に保つまで制御バルブを徐々に閉じ、調整過程全体中のユニット負荷、機械前主蒸気圧力、調圧器圧力、主蒸気温度および調圧器後温度を記録するステップ（１）と、
流量特性曲線計算：ステップ（１）で得られたデータを用いて、フリューゲル式に従って各高圧制御バルブ流量特性曲線を計算するステップ（２）と、
単一バルブ制御モードの各分岐高圧制御バルブのパラメータ補正：単一バルブ制御モードで、実際の単一バルブ正規化相対流量を計算し、実際の単一バルブ正規化相対流量をｙ、対応のバルブ開度をｘとしてフィッティングし、高次連続関数ｙ＝ｆ(ｘ)、ＤＣＳで現在ユニット単一バルブで使用されるポリゴン関数を見つけ、高次連続関数ｙ＝ｆ(ｘ)を利用して総合バルブ位置指令またはバルブ開度を補正するステップ（３）と、
シーケンスバルブ制御モードの各分岐高圧制御バルブのパラメータ補正：シーケンスバルブ制御モードで、すべての制御バルブの全開から全閉までａ個のステップがあり、バルブ全開の場合、最終段階全開まで閉じ順序に従ってバルブを順次閉じる過程中、発生した変曲点数がａ－１であり、変曲点に対応する相対流量ｙｍを計算するステップ（４）と、を含み、
ＤＣＳでシーケンスバルブ制御モードで、各支制御バルブで使用されるポリゴン関数を見つけ、該関数において、Ｙｑｉは離散統合されたバルブ位置指令であり、Ｘｑｉは対応のバルブ開度であり、
実際の各高圧制御バルブ正規化相対流量をｙ、対応のバルブ開度をｘとしてフィッティングし、各高圧制御バルブ高次連続関数ｙｉ＝ｆｉ(ｘ)を得、ｉはｉ番目の高圧制御バルブであり、
変曲点に対応する相対流量ｙｍおよび各高圧制御バルブ高次連続関数ｙ＝ｆ(ｘ)、ポリゴン関数に基づいて、総合バルブ位置と相対流量の対照関係を計算して得、該対照関係に基づいて、総合バルブ位置指令またはバルブ開度を補正する。
【0007】
さらに、好ましくは、ステップ（１）において、記録時、記録サンプリング間隔は１ｓである。
【0008】
さらに、好ましくは、ステップ（２）は、具体的に、
ユニットがｉ番目の高圧調整バルブ開度ｊ下で、ｊ∈［０、１００％］、該高圧調整バルブ前後の
高圧調整バルブ前後の蒸気温度に基づいて補正し、
ここで、ＰＴｉｊはｉ番目の高圧制御バルブ開度ｊの時の機械前主蒸気圧力であり、単位ＭＰａ、
Ｐ１ｉｊはｉ番目の高圧制御バルブ開度ｊの時の調圧器圧力であり、単位ＭＰａ、
ＴＳｉｊはｉ番目の高圧制御バルブ開度ｊの時の機械前主蒸気温度であり、単位℃、
Ｔ１ｉｊはｉ番目の高圧制御バルブ開度ｊの時の調圧器後温度であり、単位℃である。
ユニットがｉ番目の高圧調整バルブ開度ｊ下で、ｊ∈［０、１００％］、該高圧調整バルブ前後の
Ｐ１ｉｊはｉ番目の高圧制御バルブ開度ｊの時の調圧器圧力であり、単位ＭＰａ、
ＴＳｉｊはｉ番目の高圧制御バルブ開度ｊの時の機械前主蒸気温度であり、単位℃、
Ｔ１ｉｊはｉ番目の高圧制御バルブ開度ｊの時の調圧器後温度であり、単位℃である。
【0009】
試験範囲内の相対流量を０～１００％正規化し：
ここで、ＦＲｉｊ％はｉ番目の高圧制御バルブ開度ｊの時の相対流量であり、
ＦＲｉ０はｉ番目の高圧制御バルブのバルブ開度０％時の補正後蒸気圧力比であり、
ＦＲｉ１００はｉ番目の高圧制御バルブのバルブ開度１００％時の補正後蒸気圧力比であり、
得られた相対流量に基づいて、ｉ番目の高圧制御バルブ流量特性曲線を描く。
ＦＲｉ０はｉ番目の高圧制御バルブのバルブ開度０％時の補正後蒸気圧力比であり、
ＦＲｉ１００はｉ番目の高圧制御バルブのバルブ開度１００％時の補正後蒸気圧力比であり、
得られた相対流量に基づいて、ｉ番目の高圧制御バルブ流量特性曲線を描く。
【0010】
さらに、好ましくは、ステップ（３）において、実際の単一バルブ正規化相対流量を計算することは、具体的に、
ここで、
は実際の単一バルブ正規化相対流量であり、ｉは高圧調整バルブ数であり、ｊは高圧調整バルブ開度であり、ｊ∈［０、１００％］、
ＤＣＳで現在ユニット単一バルブで使用されるポリゴン関数を見つけ、該関数において、Ｙｎは離散統合されたバルブ位置指令であり、ＸｎはＹｎに対応するバルブ開度であり、実際の単一バルブ正規化相対流量
をｙ、対応のバルブ開度をｘとしてフィッティングし、高次連続関数ｙ＝ｆ(ｘ)を得、
ポリゴン関数中のＸｎを上記フィッティング後の高次連続関数に代入し、補正後の総合バルブ位置指令Ｙｎ’を得、または逆関数法ｘ＝ｆ－１(ｙ)により、ポリゴン関数中のＹｎを高次連続関数を解くことに代入し、補正後のバルブ開度Ｘｎ’を得る。
ＤＣＳで現在ユニット単一バルブで使用されるポリゴン関数を見つけ、該関数において、Ｙｎは離散統合されたバルブ位置指令であり、ＸｎはＹｎに対応するバルブ開度であり、実際の単一バルブ正規化相対流量
ポリゴン関数中のＸｎを上記フィッティング後の高次連続関数に代入し、補正後の総合バルブ位置指令Ｙｎ’を得、または逆関数法ｘ＝ｆ－１(ｙ)により、ポリゴン関数中のＹｎを高次連続関数を解くことに代入し、補正後のバルブ開度Ｘｎ’を得る。
【0011】
さらに、好ましくは、ステップ（４）の方法は具体的に、
シーケンスバルブ制御モードで、すべての制御バルブの全開から全閉までａ個のステップがあり、バルブ全開の場合、最終段階全開まで閉じ順序に従ってバルブを順次閉じる過程中、発生した変曲点数がａ－１であり、最初に閉じたバルブに対応する変曲点番号がｍ＝１であり、２番目に閉じたバルブに対応する変曲点番号がｍ＝２であり、以降同様であり、最後の変曲点番号がｍ＝ａ－１であり、
変曲点に対応する相対流量ｙｍを計算し、
式において、Ｐｍはｍ番目の変曲点に対応するユニット負荷であり、Ｐｅはバルブ全開時のユニット負荷であり、
ＤＣＳでシーケンスバルブ制御モードで、各支制御バルブで使用されるポリゴン関数を見つけ、該関数において、Ｙｑｉは離散統合されたバルブ位置指令であり、Ｘｑｉは対応のバルブ開度であり、
実際の各高圧制御バルブ正規化相対流量ＦＲｉｊ％をｙ、対応のバルブ開度をｘとしてフィッティングし、各高圧制御バルブ高次連続関数ｙｉ＝ｆｉ(ｘ)を得、ｉはｉ番目の高圧制御バルブであり、
変曲点に対応する相対流量ｙｍおよび各高圧制御バルブ高次連続関数ｙ＝ｆ(ｘ)、ポリゴン関数に基づいて、総合バルブ位置と相対流量の対照関係を計算して得ることは、具体的に、
１、各変曲点の総合バルブ位置指令および相対流量を決定し、すなわち、ポリゴン関数からすべての変曲点の総合バルブ位置指令を見つけ、計算された変曲点に対応する相対流量に対応させ、
２、総合バルブ位置が［０,ａ－１番目の変曲点に対応する総合バルブ位置指令］の範囲内にある場合、該総合バルブ位置をｘとしてこの範囲内で開く必要のある各高圧制御バルブの高次連続関数に代入し、対応のｙｉ値を求め、その後、該総合バルブ位置下でこの範囲内で開く必要のあるすべての高圧制御バルブのｙｉ値の平均値を求め、ａ－１番目の変曲点に対応する相対流量を乗じることで、該総合バルブ位置に対応する実際の相対流量を得、
３、総合バルブ位置が(ａ－１番目の変曲点に対応する総合バルブ位置指令,ａ－２番目の変曲点に対応する総合バルブ位置指令］の範囲内にある場合、該総合バルブ位置をｘとしてこの範囲内で開く必要のある各高圧制御バルブの高次連続関数に代入し、対応のｙｉ値を求め、ここで、全開した高圧制御バルブを除外し、その後、該総合バルブ位置下でこの範囲内で開く必要のあるすべての高圧制御バルブのｙｉ値の平均値を求め、（ａ－２番目の変曲点に対応する相対流量－ａ－１番目の変曲点に対応する相対流量）＋ａ－１番目の変曲点に対応する相対流量を乗じ、該総合バルブ位置に対応する実際の相対流量を得、以降同様であり、
４、総合バルブ位置が(第１変曲点に対応する総合バルブ位置指令,１］の範囲内にある場合、該総合バルブ位置をｘとしてこの範囲内で開く必要のある各高圧制御バルブの高次連続関数に代入し、対応のｙｉ値を求め、ここで、全開した高圧制御バルブを除外し、その後、該総合バルブ位置下でこの範囲内で開く必要のあるすべての高圧制御バルブのｙｉ値の平均値を求め、（第１変曲点に対応する相対流量－第２変曲点に対応する相対流量）＋第２変曲点に対応する相対流量を乗じ、該総合バルブ位置に対応する実際の相対流量を得、
その後、該対照関係に基づいて、総合バルブ位置指令またはバルブ開度を補正する。
シーケンスバルブ制御モードで、すべての制御バルブの全開から全閉までａ個のステップがあり、バルブ全開の場合、最終段階全開まで閉じ順序に従ってバルブを順次閉じる過程中、発生した変曲点数がａ－１であり、最初に閉じたバルブに対応する変曲点番号がｍ＝１であり、２番目に閉じたバルブに対応する変曲点番号がｍ＝２であり、以降同様であり、最後の変曲点番号がｍ＝ａ－１であり、
変曲点に対応する相対流量ｙｍを計算し、
ＤＣＳでシーケンスバルブ制御モードで、各支制御バルブで使用されるポリゴン関数を見つけ、該関数において、Ｙｑｉは離散統合されたバルブ位置指令であり、Ｘｑｉは対応のバルブ開度であり、
実際の各高圧制御バルブ正規化相対流量ＦＲｉｊ％をｙ、対応のバルブ開度をｘとしてフィッティングし、各高圧制御バルブ高次連続関数ｙｉ＝ｆｉ(ｘ)を得、ｉはｉ番目の高圧制御バルブであり、
変曲点に対応する相対流量ｙｍおよび各高圧制御バルブ高次連続関数ｙ＝ｆ(ｘ)、ポリゴン関数に基づいて、総合バルブ位置と相対流量の対照関係を計算して得ることは、具体的に、
１、各変曲点の総合バルブ位置指令および相対流量を決定し、すなわち、ポリゴン関数からすべての変曲点の総合バルブ位置指令を見つけ、計算された変曲点に対応する相対流量に対応させ、
２、総合バルブ位置が［０,ａ－１番目の変曲点に対応する総合バルブ位置指令］の範囲内にある場合、該総合バルブ位置をｘとしてこの範囲内で開く必要のある各高圧制御バルブの高次連続関数に代入し、対応のｙｉ値を求め、その後、該総合バルブ位置下でこの範囲内で開く必要のあるすべての高圧制御バルブのｙｉ値の平均値を求め、ａ－１番目の変曲点に対応する相対流量を乗じることで、該総合バルブ位置に対応する実際の相対流量を得、
３、総合バルブ位置が(ａ－１番目の変曲点に対応する総合バルブ位置指令,ａ－２番目の変曲点に対応する総合バルブ位置指令］の範囲内にある場合、該総合バルブ位置をｘとしてこの範囲内で開く必要のある各高圧制御バルブの高次連続関数に代入し、対応のｙｉ値を求め、ここで、全開した高圧制御バルブを除外し、その後、該総合バルブ位置下でこの範囲内で開く必要のあるすべての高圧制御バルブのｙｉ値の平均値を求め、（ａ－２番目の変曲点に対応する相対流量－ａ－１番目の変曲点に対応する相対流量）＋ａ－１番目の変曲点に対応する相対流量を乗じ、該総合バルブ位置に対応する実際の相対流量を得、以降同様であり、
４、総合バルブ位置が(第１変曲点に対応する総合バルブ位置指令,１］の範囲内にある場合、該総合バルブ位置をｘとしてこの範囲内で開く必要のある各高圧制御バルブの高次連続関数に代入し、対応のｙｉ値を求め、ここで、全開した高圧制御バルブを除外し、その後、該総合バルブ位置下でこの範囲内で開く必要のあるすべての高圧制御バルブのｙｉ値の平均値を求め、（第１変曲点に対応する相対流量－第２変曲点に対応する相対流量）＋第２変曲点に対応する相対流量を乗じ、該総合バルブ位置に対応する実際の相対流量を得、
その後、該対照関係に基づいて、総合バルブ位置指令またはバルブ開度を補正する。
【0012】
さらに、好ましくは、対照関係において理論曲線からの偏差が最も大きい点を補正する。
【0013】
本発明は同時にＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正するシステムを提供し、上記ＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正する方法を採用し、
試験過程中、ユニット負荷、機械前主蒸気圧力、調圧器圧力、主蒸気流量、給水流量、主蒸気脱熱水流量、主蒸気温度、再熱蒸気温度、蒸気ドラム圧力および蒸気タービン全エネルギー流量を収集するように構成されたデータ収集モジュールと、
データ収集モジュールによって収集されたデータに基づいて、フリューゲル式を利用して各高圧制御バルブ流量特性曲線を計算するように構成されたデータ処理モジュールと、
実際の単一バルブ正規化相対流量を計算し、実際の単一バルブ正規化相対流量をｙ、対応のバルブ開度をｘとしてフィッティングし、高次連続関数ｙ＝ｆ(ｘ)を得、ＤＣＳで現在ユニット単一バルブで使用されるポリゴン関数を見つけ、高次連続関数ｙ＝ｆ(ｘ)を利用して総合バルブ位置指令またはバルブ開度を補正するように構成された単一バルブ制御モードの各分岐高圧制御バルブのパラメータ補正モジュールと、
変曲点に対応する相対流量ｙｍを計算し、ＤＣＳでシーケンスバルブ制御モードで、各支制御バルブで使用されるポリゴン関数を見つけ、該関数において、Ｙｑｉは離散統合されたバルブ位置指令であり、Ｘｑｉは対応のバルブ開度であり、実際の各高圧制御バルブ正規化相対流量をｙ、対応のバルブ開度をｘとしてフィッティングし、各高圧制御バルブ高次連続関数ｙｉ＝ｆｉ(ｘ)を得、ｉはｉ番目の高圧制御バルブであり、変曲点に対応する相対流量ｙｍおよび各高圧制御バルブ高次連続関数ｙ＝ｆ(ｘ)、ポリゴン関数に基づいて、総合バルブ位置と相対流量の対照関係を計算して得、該対照関係に基づいて、総合バルブ位置指令またはバルブ開度を補正するように構成されたシーケンスバルブ制御モードの各分岐高圧制御バルブのパラメータ補正モジュールと、を備える。
試験過程中、ユニット負荷、機械前主蒸気圧力、調圧器圧力、主蒸気流量、給水流量、主蒸気脱熱水流量、主蒸気温度、再熱蒸気温度、蒸気ドラム圧力および蒸気タービン全エネルギー流量を収集するように構成されたデータ収集モジュールと、
データ収集モジュールによって収集されたデータに基づいて、フリューゲル式を利用して各高圧制御バルブ流量特性曲線を計算するように構成されたデータ処理モジュールと、
実際の単一バルブ正規化相対流量を計算し、実際の単一バルブ正規化相対流量をｙ、対応のバルブ開度をｘとしてフィッティングし、高次連続関数ｙ＝ｆ(ｘ)を得、ＤＣＳで現在ユニット単一バルブで使用されるポリゴン関数を見つけ、高次連続関数ｙ＝ｆ(ｘ)を利用して総合バルブ位置指令またはバルブ開度を補正するように構成された単一バルブ制御モードの各分岐高圧制御バルブのパラメータ補正モジュールと、
変曲点に対応する相対流量ｙｍを計算し、ＤＣＳでシーケンスバルブ制御モードで、各支制御バルブで使用されるポリゴン関数を見つけ、該関数において、Ｙｑｉは離散統合されたバルブ位置指令であり、Ｘｑｉは対応のバルブ開度であり、実際の各高圧制御バルブ正規化相対流量をｙ、対応のバルブ開度をｘとしてフィッティングし、各高圧制御バルブ高次連続関数ｙｉ＝ｆｉ(ｘ)を得、ｉはｉ番目の高圧制御バルブであり、変曲点に対応する相対流量ｙｍおよび各高圧制御バルブ高次連続関数ｙ＝ｆ(ｘ)、ポリゴン関数に基づいて、総合バルブ位置と相対流量の対照関係を計算して得、該対照関係に基づいて、総合バルブ位置指令またはバルブ開度を補正するように構成されたシーケンスバルブ制御モードの各分岐高圧制御バルブのパラメータ補正モジュールと、を備える。
【0014】
先行技術と比較すると、本発明は以下の有益な効果を有する。
【0015】
現在ＤＣＳ（Ｄｉｓｔｒｉｂｕｔｅｄ Ｃｏｎｔｒｏｌ Ｓｙｓｔｅｍ、分散制御システム）で使用されるバルブ流量特性曲線は実際のバルブ流量特性と一致せず、ユニットの総合バルブ位置指令変化は発生した蒸気流量およびユニット負荷に対して大きな非線形を有し、一次周波数変調の過度調整または調整不足が深刻になり、協調制御システムは逆調整効果を生じ、ひいては系統低周波発振を誘発し、本発明は、ＤＥＨ（Ｄｉｇｉｔａｌ Ｅｌｅｃｔｒｏ－Ｈｙｄｒａｕｌｉｃ、蒸気タービン用デジタル電気油圧制御システム）バルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正する方法を提案する。
【0016】
本発明の方法は実地試験データに基づいており、得られた結果はより的を射た実用的なものであり、ユニットの実際の状況によく適応することができ、本発明の方法の試験時間は約２時間であり、試験工程は簡単で操作しやすく、相対流量計アルゴリズムとカーブフィッティングの方法は、他の方法と比べて、より明確で、直感的で、時間を節約することができ、異なる火力発電ユニットに適用でき、普遍的に適用できる。
【図面の簡単な説明】
【0017】
【図1】蒸気タービンバルブノズルの配置を示す図である。
【図2】正規化後総合バルブ位置指令と相対流量の理論および実際の対応関係を示す図である。
【図3】本発明のＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正するシステムの構造を示す概略図である。矢印方向はデータまたは信号方向である。
【発明を実施するための形態】
【0018】
以下で、実施例に関連して本発明をさらに詳細に説明する。
【0019】
当業者であれば、以下の実施例は本発明を説明する目的でのみ使用され、本発明の範囲を限定するものではないことを理解されたい。実施例では、具体的な技術または条件が示されていない場合、本分野内の文献に記載されている技術または条件または製品仕様書に従うものである。使用される材料または装置にメーカの記載がない場合、購入により入手可能な従来品である。
【0020】
実施例１
ＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正する方法は、
データ収集：ユニットをＤＥＨバルブ制御モードに切り替え、ユニットの一次周波数変調制御を切り替え、ユニットのＡＧＣ制御を切り替え、機械と炉の協調制御を手動に切り替え、ユニット燃料自動を切り替え、吸気酸素自動を切り替え、ユニット主蒸気温度自動、給水自動、炉圧自動を入れ、単一高圧制御バルブ１００～０％および０～１００％の単一ステップ５％のステップ外乱を行い、ある制御バルブが動作しているとき、他の制御バルブはすべて開いており、単一高圧制御バルブの動作が完了した後すべての制御バルブを全開に保ち、その後５％のステップでシーケンスバルブ投入の逆順序に従って最終段階で制御バルブを全開に保つまで制御バルブを徐々に閉じ、調整過程全体中のユニット負荷、機械前主蒸気圧力、調圧器圧力、主蒸気温度および調圧器後温度を記録するステップ（１）と、
流量特性曲線計算：ステップ（１）で得られたデータを用いて、フリューゲル式に従って各高圧制御バルブ流量特性曲線を計算するステップ（２）と、
単一バルブ制御モードの各分岐高圧制御バルブのパラメータ補正：単一バルブ制御モードで、実際の単一バルブ正規化相対流量を計算し、実際の単一バルブ正規化相対流量をｙ、対応のバルブ開度をｘとしてフィッティングし、高次連続関数ｙ＝ｆ(ｘ)、ＤＣＳで現在ユニット単一バルブで使用されるポリゴン関数を見つけ、高次連続関数ｙ＝ｆ(ｘ)を利用して総合バルブ位置指令またはバルブ開度を補正するステップ（３）と、
シーケンスバルブ制御モードの各分岐高圧制御バルブのパラメータ補正：シーケンスバルブ制御モードで、すべての制御バルブの全開から全閉までａ個のステップがあり、バルブ全開の場合、最終段階全開まで閉じ順序に従ってバルブを順次閉じる過程中、発生した変曲点数がａ－１であり、変曲点に対応する相対流量ｙｍを計算するステップ（４）と、を含み、
ＤＣＳでシーケンスバルブ制御モードで、各支制御バルブで使用されるポリゴン関数を見つけ、該関数において、Ｙｑｉは離散統合されたバルブ位置指令であり、Ｘｑｉは対応のバルブ開度であり、
実際の各高圧制御バルブ正規化相対流量をｙ、対応のバルブ開度をｘとしてフィッティングし、各高圧制御バルブ高次連続関数ｙｉ＝ｆｉ(ｘ)を得、ｉはｉ番目の高圧制御バルブであり、
変曲点に対応する相対流量ｙｍおよび各高圧制御バルブ高次連続関数ｙ＝ｆ(ｘ)、ポリゴン関数に基づいて、総合バルブ位置と相対流量の対照関係を計算して得、該対照関係に基づいて、総合バルブ位置指令またはバルブ開度を補正する。
ＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正する方法は、
データ収集：ユニットをＤＥＨバルブ制御モードに切り替え、ユニットの一次周波数変調制御を切り替え、ユニットのＡＧＣ制御を切り替え、機械と炉の協調制御を手動に切り替え、ユニット燃料自動を切り替え、吸気酸素自動を切り替え、ユニット主蒸気温度自動、給水自動、炉圧自動を入れ、単一高圧制御バルブ１００～０％および０～１００％の単一ステップ５％のステップ外乱を行い、ある制御バルブが動作しているとき、他の制御バルブはすべて開いており、単一高圧制御バルブの動作が完了した後すべての制御バルブを全開に保ち、その後５％のステップでシーケンスバルブ投入の逆順序に従って最終段階で制御バルブを全開に保つまで制御バルブを徐々に閉じ、調整過程全体中のユニット負荷、機械前主蒸気圧力、調圧器圧力、主蒸気温度および調圧器後温度を記録するステップ（１）と、
流量特性曲線計算：ステップ（１）で得られたデータを用いて、フリューゲル式に従って各高圧制御バルブ流量特性曲線を計算するステップ（２）と、
単一バルブ制御モードの各分岐高圧制御バルブのパラメータ補正：単一バルブ制御モードで、実際の単一バルブ正規化相対流量を計算し、実際の単一バルブ正規化相対流量をｙ、対応のバルブ開度をｘとしてフィッティングし、高次連続関数ｙ＝ｆ(ｘ)、ＤＣＳで現在ユニット単一バルブで使用されるポリゴン関数を見つけ、高次連続関数ｙ＝ｆ(ｘ)を利用して総合バルブ位置指令またはバルブ開度を補正するステップ（３）と、
シーケンスバルブ制御モードの各分岐高圧制御バルブのパラメータ補正：シーケンスバルブ制御モードで、すべての制御バルブの全開から全閉までａ個のステップがあり、バルブ全開の場合、最終段階全開まで閉じ順序に従ってバルブを順次閉じる過程中、発生した変曲点数がａ－１であり、変曲点に対応する相対流量ｙｍを計算するステップ（４）と、を含み、
ＤＣＳでシーケンスバルブ制御モードで、各支制御バルブで使用されるポリゴン関数を見つけ、該関数において、Ｙｑｉは離散統合されたバルブ位置指令であり、Ｘｑｉは対応のバルブ開度であり、
実際の各高圧制御バルブ正規化相対流量をｙ、対応のバルブ開度をｘとしてフィッティングし、各高圧制御バルブ高次連続関数ｙｉ＝ｆｉ(ｘ)を得、ｉはｉ番目の高圧制御バルブであり、
変曲点に対応する相対流量ｙｍおよび各高圧制御バルブ高次連続関数ｙ＝ｆ(ｘ)、ポリゴン関数に基づいて、総合バルブ位置と相対流量の対照関係を計算して得、該対照関係に基づいて、総合バルブ位置指令またはバルブ開度を補正する。
【0021】
実施例２
ＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正する方法は、
データ収集：ユニットをＤＥＨバルブ制御モードに切り替え、ユニットの一次周波数変調制御を切り替え、ユニットのＡＧＣ制御を切り替え、機械と炉の協調制御を手動に切り替え、ユニット燃料自動を切り替え、吸気酸素自動を切り替え、ユニット主蒸気温度自動、給水自動、炉圧自動を入れ、単一高圧制御バルブ１００～０％および０～１００％の単一ステップ５％のステップ外乱を行い、ある制御バルブが動作しているとき、他の制御バルブはすべて開いており、単一高圧制御バルブの動作が完了した後すべての制御バルブを全開に保ち、その後５％のステップでシーケンスバルブ投入の逆順序に従って最終段階で制御バルブを全開に保つまで制御バルブを徐々に閉じ、調整過程全体中のユニット負荷、機械前主蒸気圧力、調圧器圧力、主蒸気温度および調圧器後温度を記録するステップ（１）と、
流量特性曲線計算：ステップ（１）で得られたデータを用いて、フリューゲル式に従って各高圧制御バルブ流量特性曲線を計算するステップ（２）と、
単一バルブ制御モードの各分岐高圧制御バルブのパラメータ補正：単一バルブ制御モードで、実際の単一バルブ正規化相対流量を計算し、実際の単一バルブ正規化相対流量をｙ、対応のバルブ開度をｘとしてフィッティングし、高次連続関数ｙ＝ｆ(ｘ)、ＤＣＳで現在ユニット単一バルブで使用されるポリゴン関数を見つけ、高次連続関数ｙ＝ｆ(ｘ)を利用して総合バルブ位置指令またはバルブ開度を補正するステップ（３）と、
シーケンスバルブ制御モードの各分岐高圧制御バルブのパラメータ補正：シーケンスバルブ制御モードで、すべての制御バルブの全開から全閉までａ個のステップがあり、バルブ全開の場合、最終段階全開まで閉じ順序に従ってバルブを順次閉じる過程中、発生した変曲点数がａ－１であり、変曲点に対応する相対流量ｙｍを計算するステップ（４）と、を含み、
ＤＣＳでシーケンスバルブ制御モードで、各支制御バルブで使用されるポリゴン関数を見つけ、該関数において、Ｙｑｉは離散統合されたバルブ位置指令であり、Ｘｑｉは対応のバルブ開度であり、
実際の各高圧制御バルブ正規化相対流量をｙ、対応のバルブ開度をｘとしてフィッティングし、各高圧制御バルブ高次連続関数ｙｉ＝ｆｉ(ｘ)を得、ｉはｉ番目の高圧制御バルブであり、
変曲点に対応する相対流量ｙｍおよび各高圧制御バルブ高次連続関数ｙ＝ｆ(ｘ)、ポリゴン関数に基づいて、総合バルブ位置と相対流量の対照関係を計算して得、該対照関係に基づいて、総合バルブ位置指令またはバルブ開度を補正する。
ＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正する方法は、
データ収集：ユニットをＤＥＨバルブ制御モードに切り替え、ユニットの一次周波数変調制御を切り替え、ユニットのＡＧＣ制御を切り替え、機械と炉の協調制御を手動に切り替え、ユニット燃料自動を切り替え、吸気酸素自動を切り替え、ユニット主蒸気温度自動、給水自動、炉圧自動を入れ、単一高圧制御バルブ１００～０％および０～１００％の単一ステップ５％のステップ外乱を行い、ある制御バルブが動作しているとき、他の制御バルブはすべて開いており、単一高圧制御バルブの動作が完了した後すべての制御バルブを全開に保ち、その後５％のステップでシーケンスバルブ投入の逆順序に従って最終段階で制御バルブを全開に保つまで制御バルブを徐々に閉じ、調整過程全体中のユニット負荷、機械前主蒸気圧力、調圧器圧力、主蒸気温度および調圧器後温度を記録するステップ（１）と、
流量特性曲線計算：ステップ（１）で得られたデータを用いて、フリューゲル式に従って各高圧制御バルブ流量特性曲線を計算するステップ（２）と、
単一バルブ制御モードの各分岐高圧制御バルブのパラメータ補正：単一バルブ制御モードで、実際の単一バルブ正規化相対流量を計算し、実際の単一バルブ正規化相対流量をｙ、対応のバルブ開度をｘとしてフィッティングし、高次連続関数ｙ＝ｆ(ｘ)、ＤＣＳで現在ユニット単一バルブで使用されるポリゴン関数を見つけ、高次連続関数ｙ＝ｆ(ｘ)を利用して総合バルブ位置指令またはバルブ開度を補正するステップ（３）と、
シーケンスバルブ制御モードの各分岐高圧制御バルブのパラメータ補正：シーケンスバルブ制御モードで、すべての制御バルブの全開から全閉までａ個のステップがあり、バルブ全開の場合、最終段階全開まで閉じ順序に従ってバルブを順次閉じる過程中、発生した変曲点数がａ－１であり、変曲点に対応する相対流量ｙｍを計算するステップ（４）と、を含み、
ＤＣＳでシーケンスバルブ制御モードで、各支制御バルブで使用されるポリゴン関数を見つけ、該関数において、Ｙｑｉは離散統合されたバルブ位置指令であり、Ｘｑｉは対応のバルブ開度であり、
実際の各高圧制御バルブ正規化相対流量をｙ、対応のバルブ開度をｘとしてフィッティングし、各高圧制御バルブ高次連続関数ｙｉ＝ｆｉ(ｘ)を得、ｉはｉ番目の高圧制御バルブであり、
変曲点に対応する相対流量ｙｍおよび各高圧制御バルブ高次連続関数ｙ＝ｆ(ｘ)、ポリゴン関数に基づいて、総合バルブ位置と相対流量の対照関係を計算して得、該対照関係に基づいて、総合バルブ位置指令またはバルブ開度を補正する。
【0022】
ステップ（１）において、記録時、記録サンプリング間隔は１ｓであるユニットがｉ番目の高圧調整バルブ開度ｊ下で、ｊ∈［０、１００％］、該高圧調整バルブ前後の
高圧調整バルブ前後の蒸気温度に基づいて補正し、補正後の
ここで、ＰＴｉｊはｉ番目の高圧制御バルブ開度ｊの時の機械前主蒸気圧力であり、単位ＭＰａ、
Ｐ１ｉｊはｉ番目の高圧制御バルブ開度ｊの時の調圧器圧力であり、単位ＭＰａ、
ＴＳｉｊはｉ番目の高圧制御バルブ開度ｊの時の機械前主蒸気温度であり、単位℃、
Ｔ１ｉｊはｉ番目の高圧制御バルブ開度ｊの時の調圧器後温度であり、単位℃である。
Ｐ１ｉｊはｉ番目の高圧制御バルブ開度ｊの時の調圧器圧力であり、単位ＭＰａ、
ＴＳｉｊはｉ番目の高圧制御バルブ開度ｊの時の機械前主蒸気温度であり、単位℃、
Ｔ１ｉｊはｉ番目の高圧制御バルブ開度ｊの時の調圧器後温度であり、単位℃である。
【0023】
試験範囲内の相対流量を０～１００％正規化し：
ここで、ＦＲｉｊ％はｉ番目の高圧制御バルブ開度ｊの時の相対流量であり、
ＦＲｉ０はｉ番目の高圧制御バルブのバルブ開度０％時の補正後蒸気圧力比であり、
ＦＲｉ１００はｉ番目の高圧制御バルブのバルブ開度１００％時の補正後蒸気圧力比であり、
得られた相対流量に基づいて、ｉ番目の高圧制御バルブ流量特性曲線を描く。
ＦＲｉ０はｉ番目の高圧制御バルブのバルブ開度０％時の補正後蒸気圧力比であり、
ＦＲｉ１００はｉ番目の高圧制御バルブのバルブ開度１００％時の補正後蒸気圧力比であり、
得られた相対流量に基づいて、ｉ番目の高圧制御バルブ流量特性曲線を描く。
【0024】
ステップ（３）において、実際の単一バルブ正規化相対流量を計算することは、具体的に、
ここで、
は実際の単一バルブ正規化相対流量であり、ｉは高圧調整バルブ数であり、ｊは高圧調整バルブ開度であり、ｊ∈［０、１００％］、
ＤＣＳで現在ユニット単一バルブで使用されるポリゴン関数を見つけ、該関数において、Ｙｎは離散統合されたバルブ位置指令であり、ＸｎはＹｎに対応するバルブ開度であり、実際の単一バルブ正規化相対流量
をｙ、対応のバルブ開度をｘとしてフィッティングし、高次連続関数ｙ＝ｆ(ｘ)を得、
ポリゴン関数中のＸｎを上記フィッティング後の高次連続関数に代入し、補正後の総合バルブ位置指令Ｙｎ’を得、または逆関数法ｘ＝ｆ－１(ｙ)により、ポリゴン関数中のＹｎを高次連続関数を解くことに代入し、補正後のバルブ開度Ｘｎ’を得る。
ＤＣＳで現在ユニット単一バルブで使用されるポリゴン関数を見つけ、該関数において、Ｙｎは離散統合されたバルブ位置指令であり、ＸｎはＹｎに対応するバルブ開度であり、実際の単一バルブ正規化相対流量
ポリゴン関数中のＸｎを上記フィッティング後の高次連続関数に代入し、補正後の総合バルブ位置指令Ｙｎ’を得、または逆関数法ｘ＝ｆ－１(ｙ)により、ポリゴン関数中のＹｎを高次連続関数を解くことに代入し、補正後のバルブ開度Ｘｎ’を得る。
【0025】
ステップ（４）の具体方法：
シーケンスバルブ制御モードで、すべての制御バルブの全開から全閉までａ個のステップがあり、バルブ全開の場合、最終段階全開まで閉じ順序に従ってバルブを順次閉じる過程中、発生した変曲点数がａ－１であり、最初に閉じたバルブに対応する変曲点番号がｍ＝１であり、２番目に閉じたバルブに対応する変曲点番号がｍ＝２であり、以降同様であり、最後の変曲点番号がｍ＝ａ－１であり、
変曲点に対応する相対流量ｙｍを計算し、
式において、Ｐｍはｍ番目の変曲点に対応するユニット負荷であり、Ｐｅはバルブ全開時のユニット負荷であり、
ＤＣＳでシーケンスバルブ制御モードで、各支制御バルブで使用されるポリゴン関数を見つけ、該関数において、Ｙｑｉは離散統合されたバルブ位置指令であり、Ｘｑｉは対応のバルブ開度であり、
実際の各高圧制御バルブ正規化相対流量ＦＲｉｊ％をｙ、対応のバルブ開度をｘとしてフィッティングし、各高圧制御バルブ高次連続関数ｙｉ＝ｆｉ(ｘ)を得、ｉはｉ番目の高圧制御バルブであり、
変曲点に対応する相対流量ｙｍおよび各高圧制御バルブ高次連続関数ｙ＝ｆ(ｘ)、ポリゴン関数に基づいて、総合バルブ位置と相対流量の対照関係を計算して得ることは、具体的に、
１、各変曲点の総合バルブ位置指令および相対流量を決定し、すなわち、ポリゴン関数からすべての変曲点の総合バルブ位置指令を見つけ、計算された変曲点に対応する相対流量に対応させ、
２、総合バルブ位置が［０,ａ－１番目の変曲点に対応する総合バルブ位置指令］の範囲内にある場合、該総合バルブ位置をｘとしてこの範囲内で開く必要のある各高圧制御バルブの高次連続関数に代入し、対応のｙｉ値を求め、その後、該総合バルブ位置下でこの範囲内で開く必要のあるすべての高圧制御バルブのｙｉ値の平均値を求め、ａ－１番目の変曲点に対応する相対流量を乗じることで、該総合バルブ位置に対応する実際の相対流量を得、
３、総合バルブ位置が(ａ－１番目の変曲点に対応する総合バルブ位置指令,ａ－２番目の変曲点に対応する総合バルブ位置指令］の範囲内にある場合、該総合バルブ位置をｘとしてこの範囲内で開く必要のある各高圧制御バルブの高次連続関数に代入し、対応のｙｉ値を求め、ここで、全開した高圧制御バルブを除外し、その後、該総合バルブ位置下でこの範囲内で開く必要のあるすべての高圧制御バルブのｙｉ値の平均値を求め、（ａ－２番目の変曲点に対応する相対流量－ａ－１番目の変曲点に対応する相対流量）＋ａ－１番目の変曲点に対応する相対流量を乗じ、該総合バルブ位置に対応する実際の相対流量を得、以降同様であり、
４、総合バルブ位置が(第１変曲点に対応する総合バルブ位置指令,１］の範囲内にある場合、該総合バルブ位置をｘとしてこの範囲内で開く必要のある各高圧制御バルブの高次連続関数に代入し、対応のｙｉ値を求め、ここで、全開した高圧制御バルブを除外し、その後、該総合バルブ位置下でこの範囲内で開く必要のあるすべての高圧制御バルブのｙｉ値の平均値を求め、（第１変曲点に対応する相対流量－第２変曲点に対応する相対流量）＋第２変曲点に対応する相対流量を乗じ、該総合バルブ位置に対応する実際の相対流量を得、
その後、該対照関係に基づいて、総合バルブ位置指令またはバルブ開度を補正する。
シーケンスバルブ制御モードで、すべての制御バルブの全開から全閉までａ個のステップがあり、バルブ全開の場合、最終段階全開まで閉じ順序に従ってバルブを順次閉じる過程中、発生した変曲点数がａ－１であり、最初に閉じたバルブに対応する変曲点番号がｍ＝１であり、２番目に閉じたバルブに対応する変曲点番号がｍ＝２であり、以降同様であり、最後の変曲点番号がｍ＝ａ－１であり、
変曲点に対応する相対流量ｙｍを計算し、
ＤＣＳでシーケンスバルブ制御モードで、各支制御バルブで使用されるポリゴン関数を見つけ、該関数において、Ｙｑｉは離散統合されたバルブ位置指令であり、Ｘｑｉは対応のバルブ開度であり、
実際の各高圧制御バルブ正規化相対流量ＦＲｉｊ％をｙ、対応のバルブ開度をｘとしてフィッティングし、各高圧制御バルブ高次連続関数ｙｉ＝ｆｉ(ｘ)を得、ｉはｉ番目の高圧制御バルブであり、
変曲点に対応する相対流量ｙｍおよび各高圧制御バルブ高次連続関数ｙ＝ｆ(ｘ)、ポリゴン関数に基づいて、総合バルブ位置と相対流量の対照関係を計算して得ることは、具体的に、
１、各変曲点の総合バルブ位置指令および相対流量を決定し、すなわち、ポリゴン関数からすべての変曲点の総合バルブ位置指令を見つけ、計算された変曲点に対応する相対流量に対応させ、
２、総合バルブ位置が［０,ａ－１番目の変曲点に対応する総合バルブ位置指令］の範囲内にある場合、該総合バルブ位置をｘとしてこの範囲内で開く必要のある各高圧制御バルブの高次連続関数に代入し、対応のｙｉ値を求め、その後、該総合バルブ位置下でこの範囲内で開く必要のあるすべての高圧制御バルブのｙｉ値の平均値を求め、ａ－１番目の変曲点に対応する相対流量を乗じることで、該総合バルブ位置に対応する実際の相対流量を得、
３、総合バルブ位置が(ａ－１番目の変曲点に対応する総合バルブ位置指令,ａ－２番目の変曲点に対応する総合バルブ位置指令］の範囲内にある場合、該総合バルブ位置をｘとしてこの範囲内で開く必要のある各高圧制御バルブの高次連続関数に代入し、対応のｙｉ値を求め、ここで、全開した高圧制御バルブを除外し、その後、該総合バルブ位置下でこの範囲内で開く必要のあるすべての高圧制御バルブのｙｉ値の平均値を求め、（ａ－２番目の変曲点に対応する相対流量－ａ－１番目の変曲点に対応する相対流量）＋ａ－１番目の変曲点に対応する相対流量を乗じ、該総合バルブ位置に対応する実際の相対流量を得、以降同様であり、
４、総合バルブ位置が(第１変曲点に対応する総合バルブ位置指令,１］の範囲内にある場合、該総合バルブ位置をｘとしてこの範囲内で開く必要のある各高圧制御バルブの高次連続関数に代入し、対応のｙｉ値を求め、ここで、全開した高圧制御バルブを除外し、その後、該総合バルブ位置下でこの範囲内で開く必要のあるすべての高圧制御バルブのｙｉ値の平均値を求め、（第１変曲点に対応する相対流量－第２変曲点に対応する相対流量）＋第２変曲点に対応する相対流量を乗じ、該総合バルブ位置に対応する実際の相対流量を得、
その後、該対照関係に基づいて、総合バルブ位置指令またはバルブ開度を補正する。
【0026】
対照関係において理論曲線からの偏差が最も大きい点を補正する。
【0027】
図３に示すように、ＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正するシステムは、上記ＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正する方法を採用し、
試験過程中、ユニット負荷、機械前主蒸気圧力、調圧器圧力、主蒸気流量、給水流量、主蒸気脱熱水流量、主蒸気温度、再熱蒸気温度、蒸気ドラム圧力および蒸気タービン全エネルギー流量を収集するように構成されたデータ収集モジュール１０１と、
データ収集モジュールによって収集されたデータに基づいて、フリューゲル式を利用して各高圧制御バルブ流量特性曲線を計算するように構成されたデータ処理モジュール１０２と、
実際の単一バルブ正規化相対流量を計算し、実際の単一バルブ正規化相対流量をｙ、対応のバルブ開度をｘとしてフィッティングし、高次連続関数ｙ＝ｆ(ｘ)を得、ＤＣＳで現在ユニット単一バルブで使用されるポリゴン関数を見つけ、高次連続関数ｙ＝ｆ(ｘ)を利用して総合バルブ位置指令またはバルブ開度を補正するように構成された単一バルブ制御モードの各分岐高圧制御バルブのパラメータ補正モジュール１０３と、
変曲点に対応する相対流量ｙｍを計算し、ＤＣＳでシーケンスバルブ制御モードで、各支制御バルブで使用されるポリゴン関数を見つけ、該関数において、Ｙｑｉは離散統合されたバルブ位置指令であり、Ｘｑｉは対応のバルブ開度であり、実際の各高圧制御バルブ正規化相対流量をｙ、対応のバルブ開度をｘとしてフィッティングし、各高圧制御バルブ高次連続関数ｙｉ＝ｆｉ(ｘ)を得、ｉはｉ番目の高圧制御バルブであり、変曲点に対応する相対流量ｙｍおよび各高圧制御バルブ高次連続関数ｙ＝ｆ(ｘ)、ポリゴン関数に基づいて、総合バルブ位置と相対流量の対照関係を計算して得、該対照関係に基づいて、総合バルブ位置指令またはバルブ開度を補正するように構成されたシーケンスバルブ制御モードの各分岐高圧制御バルブのパラメータ補正モジュール１０４と、を備える。
試験過程中、ユニット負荷、機械前主蒸気圧力、調圧器圧力、主蒸気流量、給水流量、主蒸気脱熱水流量、主蒸気温度、再熱蒸気温度、蒸気ドラム圧力および蒸気タービン全エネルギー流量を収集するように構成されたデータ収集モジュール１０１と、
データ収集モジュールによって収集されたデータに基づいて、フリューゲル式を利用して各高圧制御バルブ流量特性曲線を計算するように構成されたデータ処理モジュール１０２と、
実際の単一バルブ正規化相対流量を計算し、実際の単一バルブ正規化相対流量をｙ、対応のバルブ開度をｘとしてフィッティングし、高次連続関数ｙ＝ｆ(ｘ)を得、ＤＣＳで現在ユニット単一バルブで使用されるポリゴン関数を見つけ、高次連続関数ｙ＝ｆ(ｘ)を利用して総合バルブ位置指令またはバルブ開度を補正するように構成された単一バルブ制御モードの各分岐高圧制御バルブのパラメータ補正モジュール１０３と、
変曲点に対応する相対流量ｙｍを計算し、ＤＣＳでシーケンスバルブ制御モードで、各支制御バルブで使用されるポリゴン関数を見つけ、該関数において、Ｙｑｉは離散統合されたバルブ位置指令であり、Ｘｑｉは対応のバルブ開度であり、実際の各高圧制御バルブ正規化相対流量をｙ、対応のバルブ開度をｘとしてフィッティングし、各高圧制御バルブ高次連続関数ｙｉ＝ｆｉ(ｘ)を得、ｉはｉ番目の高圧制御バルブであり、変曲点に対応する相対流量ｙｍおよび各高圧制御バルブ高次連続関数ｙ＝ｆ(ｘ)、ポリゴン関数に基づいて、総合バルブ位置と相対流量の対照関係を計算して得、該対照関係に基づいて、総合バルブ位置指令またはバルブ開度を補正するように構成されたシーケンスバルブ制御モードの各分岐高圧制御バルブのパラメータ補正モジュール１０４と、を備える。
【0028】
実施例３
ＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正する方法は、以下のステップを含み：
（一）現場試験。試験ユニットの有効出力をバルブ全開時の負荷Ｐｅにし、ユニットをＤＥＨバルブ制御モードに切り替え、ユニットの一次周波数変調制御を切り替え、切りユニットのＡＧＣ制御、機械と炉の協調制御を手動に切り替え、ユニット燃料自動を切り替え、吸気酸素自動を切り替え、ユニット主蒸気温度自動、給水自動、炉圧自動を入れ、単一高圧制御バルブ１００～０％および０～１００％の単一ステップ５％のステップ外乱を行い、ある制御バルブが動作しているとき、他の制御バルブはすべて開いており、単一高圧制御バルブの動作が完了した後すべての制御バルブを全開に保ち、その後５％のステップでシーケンスバルブ投入の逆順序に従って最終段階で制御バルブを全開に保つまで制御バルブを徐々に閉じ、調整過程中のユニット負荷Ｐ、機械前主蒸気圧力ＰＴ、調圧器圧力Ｐ１、主蒸気温度ＴＳ、調圧器後温度Ｔ１などの重要なパラメータを記録する。
ＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正する方法は、以下のステップを含み：
（一）現場試験。試験ユニットの有効出力をバルブ全開時の負荷Ｐｅにし、ユニットをＤＥＨバルブ制御モードに切り替え、ユニットの一次周波数変調制御を切り替え、切りユニットのＡＧＣ制御、機械と炉の協調制御を手動に切り替え、ユニット燃料自動を切り替え、吸気酸素自動を切り替え、ユニット主蒸気温度自動、給水自動、炉圧自動を入れ、単一高圧制御バルブ１００～０％および０～１００％の単一ステップ５％のステップ外乱を行い、ある制御バルブが動作しているとき、他の制御バルブはすべて開いており、単一高圧制御バルブの動作が完了した後すべての制御バルブを全開に保ち、その後５％のステップでシーケンスバルブ投入の逆順序に従って最終段階で制御バルブを全開に保つまで制御バルブを徐々に閉じ、調整過程中のユニット負荷Ｐ、機械前主蒸気圧力ＰＴ、調圧器圧力Ｐ１、主蒸気温度ＴＳ、調圧器後温度Ｔ１などの重要なパラメータを記録する。
【0029】
（二）適切なＤＥＨバルブ流量特性の数値計算方法を選択し、ユニットがｉ番目の高圧調整バルブ開度ｊ（ｊ∈［０、１００％］）下で、該高圧調整バルブ前後の
異なる境界条件は高圧調整バルブ前後の蒸気温度に基づいて補正した後の
よりも大きいはずである。
ここで、ＰＴｉｊはｉ番目の高圧制御バルブ開度ｊの時の機械前主蒸気圧力であり、単位ＭＰａ、
Ｐ１ｉｊはｉ番目の高圧制御バルブ開度ｊの時の調圧器圧力であり、単位ＭＰａ、
ＴＳｉｊはｉ番目の高圧制御バルブ開度ｊの時の機械前主蒸気温度であり、単位℃、
Ｔ１ｉｊはｉ番目の高圧制御バルブ開度ｊの時の調圧器後温度であり、単位℃である。
Ｐ１ｉｊはｉ番目の高圧制御バルブ開度ｊの時の調圧器圧力であり、単位ＭＰａ、
ＴＳｉｊはｉ番目の高圧制御バルブ開度ｊの時の機械前主蒸気温度であり、単位℃、
Ｔ１ｉｊはｉ番目の高圧制御バルブ開度ｊの時の調圧器後温度であり、単位℃である。
【0030】
試験範囲内の相対流量を０～１００％正規化し：
ここで、ＦＲｉｊ％はｉ番目の高圧制御バルブ開度ｊの時の相対流量であり、
ＦＲｉ０はｉ番目の高圧制御バルブのバルブ開度０％時の補正後蒸気圧力比であり、
ＦＲｉ１００はｉ番目の高圧制御バルブのバルブ開度１００％時の補正後蒸気圧力比である。
ＦＲｉ０はｉ番目の高圧制御バルブのバルブ開度０％時の補正後蒸気圧力比であり、
ＦＲｉ１００はｉ番目の高圧制御バルブのバルブ開度１００％時の補正後蒸気圧力比である。
【0031】
（三）単一バルブ制御モードの各分岐高圧制御バルブのパラメータを補正する。単一バルブ運転モード下で、蒸気は高圧調整バルブおよびノズル室を通過し、３６０℃全周の調圧段ベーンに入り、調圧段ベーンが均一に加熱され、調圧段ベーンの応力分布が効果的に改善され、ユニットは比較的早く負荷を変えることができ、その時、すべての調整バルブは同じ開度で部分的に開く（例えば、単一バルブモードにおける各バルブの開度は５０％で同じように開く）。この場合、実際の単一バルブ正規化相対流量は、
であり、ここで、
は実際の単一バルブ正規化相対流量であり、ｉは高圧調整バルブ数であり、ｊは高圧調整バルブ開度であり、（ｊ∈［０、１００％］）、
ＤＣＳで現在ユニット単一バルブで使用されるポリゴン関数を見つけ、該関数において、重要なパラメータは離散統合されたバルブ位置指令Ｙｎおよび対応のバルブ開度Ｘｎであり、ここでのＹｎとＸｎはすべて正規化後の値である。
ＤＣＳで現在ユニット単一バルブで使用されるポリゴン関数を見つけ、該関数において、重要なパラメータは離散統合されたバルブ位置指令Ｙｎおよび対応のバルブ開度Ｘｎであり、ここでのＹｎとＸｎはすべて正規化後の値である。
【0032】
ｍａｔｌａｂでｐｏｌｙｆｉｔ関数およびｐｏｌｙｖａｌ関数を呼び出し、離散された
を高次連続関数ｙ＝ｆ(ｘ)にフィッティングする。正規化後のＸｎを上記フィッティング後の関数に代入して補正後の総合バルブ位置指令Ｙｎ’を求めることができる。または、逆関数法ｘ＝ｆ－１(ｙ)を通じて、Ｙｎを高次連続関数の解くことに代入して補正後のバルブ開度Ｘｎ’を得ることができる。
【0033】
（四）シーケンスバルブ制御モードの各分岐高圧制御バルブのパラメータを補正する。
【0034】
シーケンスバルブ制御モードで、調整バルブを予め設定された各バルブのシーケンスバルブ流量特性曲線に従って順次開き、例えば６つの高圧調整バルブを有するハルビン蒸気タービン工場が製造した３００ＭＷユニット：開き順序は、ＧＶ４とＧＶ５同時に開き→ＧＶ６→ＧＶ３→ＧＶ２→ＧＶ１順次開き、閉じ順序は逆である。定義：シーケンスバルブ制御モードで、すべての制御バルブの全開から全閉までａ個のステップがあり、例えば上記例ではａは５である。現場試験と併せて、バルブ全開の場合、閉じ順序に従って、最後段階全開までバルブを閉じ、この過程中、発生した変曲点数がａ－１である。変曲点番号をｍとすると、最初に閉じたバルブに対応する変曲点番号がｍ＝１であり、最後の変曲点番号がｍ＝ａ－１である。計算式は以下のとおりであり：
式において、Ｐｍはｍ番目の変曲点に対応するユニット負荷であり、Ｐｅはバルブ全開時のユニット負荷であり、
ＤＣＳにおいて、シーケンスバルブ制御モードで、各支制御バルブで使用されるポリゴン関数を見つけ、該関数において、重要なパラメータは離散統合されたバルブ位置指令Ｙｑｉおよび対応のバルブ開度Ｘｑｉであり、ここでのＹｑｉとＸｑｉはすべて正規化後の値である。
ＤＣＳにおいて、シーケンスバルブ制御モードで、各支制御バルブで使用されるポリゴン関数を見つけ、該関数において、重要なパラメータは離散統合されたバルブ位置指令Ｙｑｉおよび対応のバルブ開度Ｘｑｉであり、ここでのＹｑｉとＸｑｉはすべて正規化後の値である。
【0035】
ｍａｔｌａｂでｐｏｌｙｆｉｔ関数およびｐｏｌｙｖａｌ関数を呼び出し、実際の各高圧制御バルブ正規化相対流量ＦＲｉｊ％をｙ、対応のバルブ開度をｘとしてフィッティングし、各高圧制御バルブ高次連続関数ｙｉ＝ｆｉ(ｘ)を得、ｉはｉ番目の高圧制御バルブであり、
変曲点に対応する相対流量ｙｍおよび各高圧制御バルブ高次連続関数ｙ＝ｆ(ｘ)、ポリゴン関数に基づいて、総合バルブ位置と相対流量の対照関係を計算して得ることは、具体的に、
１、各変曲点の総合バルブ位置指令および相対流量を決定し、すなわち、ポリゴン関数からすべての変曲点の総合バルブ位置指令を見つけ、計算された変曲点に対応する相対流量に対応させ、
２、総合バルブ位置が［０,ａ－１番目の変曲点に対応する総合バルブ位置指令］の範囲内にある場合、該総合バルブ位置をｘとしてこの範囲内で開く必要のある各高圧制御バルブの高次連続関数に代入し、対応のｙｉ値を求め、その後、該総合バルブ位置下でこの範囲内で開く必要のあるすべての高圧制御バルブのｙｉ値の平均値を求め、ａ－１番目の変曲点に対応する相対流量を乗じることで、該総合バルブ位置に対応する実際の相対流量を得、
３、総合バルブ位置が(ａ－１番目の変曲点に対応する総合バルブ位置指令,ａ－２番目の変曲点に対応する総合バルブ位置指令］の範囲内にある場合、該総合バルブ位置をｘとしてこの範囲内で開く必要のある各高圧制御バルブの高次連続関数に代入し、対応のｙｉ値を求め、ここで、全開した高圧制御バルブを除外し、その後、該総合バルブ位置下でこの範囲内で開く必要のあるすべての高圧制御バルブのｙｉ値の平均値を求め、（ａ－２番目の変曲点に対応する相対流量－ａ－１番目の変曲点に対応する相対流量）＋ａ－１番目の変曲点に対応する相対流量を乗じ、該総合バルブ位置に対応する実際の相対流量を得、以降同様であり、
４、総合バルブ位置が(第１変曲点に対応する総合バルブ位置指令,１］の範囲内にある場合、該総合バルブ位置をｘとしてこの範囲内で開く必要のある各高圧制御バルブの高次連続関数に代入し、対応のｙｉ値を求め、ここで、全開した高圧制御バルブを除外し、その後、該総合バルブ位置下でこの範囲内で開く必要のあるすべての高圧制御バルブのｙｉ値の平均値を求め、（第１変曲点に対応する相対流量－第２変曲点に対応する相対流量）＋第２変曲点に対応する相対流量を乗じ、該総合バルブ位置に対応する実際の相対流量を得、
その後、該対照関係に基づいて、総合バルブ位置指令またはバルブ開度を補正する。
変曲点に対応する相対流量ｙｍおよび各高圧制御バルブ高次連続関数ｙ＝ｆ(ｘ)、ポリゴン関数に基づいて、総合バルブ位置と相対流量の対照関係を計算して得ることは、具体的に、
１、各変曲点の総合バルブ位置指令および相対流量を決定し、すなわち、ポリゴン関数からすべての変曲点の総合バルブ位置指令を見つけ、計算された変曲点に対応する相対流量に対応させ、
２、総合バルブ位置が［０,ａ－１番目の変曲点に対応する総合バルブ位置指令］の範囲内にある場合、該総合バルブ位置をｘとしてこの範囲内で開く必要のある各高圧制御バルブの高次連続関数に代入し、対応のｙｉ値を求め、その後、該総合バルブ位置下でこの範囲内で開く必要のあるすべての高圧制御バルブのｙｉ値の平均値を求め、ａ－１番目の変曲点に対応する相対流量を乗じることで、該総合バルブ位置に対応する実際の相対流量を得、
３、総合バルブ位置が(ａ－１番目の変曲点に対応する総合バルブ位置指令,ａ－２番目の変曲点に対応する総合バルブ位置指令］の範囲内にある場合、該総合バルブ位置をｘとしてこの範囲内で開く必要のある各高圧制御バルブの高次連続関数に代入し、対応のｙｉ値を求め、ここで、全開した高圧制御バルブを除外し、その後、該総合バルブ位置下でこの範囲内で開く必要のあるすべての高圧制御バルブのｙｉ値の平均値を求め、（ａ－２番目の変曲点に対応する相対流量－ａ－１番目の変曲点に対応する相対流量）＋ａ－１番目の変曲点に対応する相対流量を乗じ、該総合バルブ位置に対応する実際の相対流量を得、以降同様であり、
４、総合バルブ位置が(第１変曲点に対応する総合バルブ位置指令,１］の範囲内にある場合、該総合バルブ位置をｘとしてこの範囲内で開く必要のある各高圧制御バルブの高次連続関数に代入し、対応のｙｉ値を求め、ここで、全開した高圧制御バルブを除外し、その後、該総合バルブ位置下でこの範囲内で開く必要のあるすべての高圧制御バルブのｙｉ値の平均値を求め、（第１変曲点に対応する相対流量－第２変曲点に対応する相対流量）＋第２変曲点に対応する相対流量を乗じ、該総合バルブ位置に対応する実際の相対流量を得、
その後、該対照関係に基づいて、総合バルブ位置指令またはバルブ開度を補正する。
【0036】
応用実施例
ＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正する方法は、以下のステップを含み：
（一）現場試験を行い、ある３００ＭＷ石炭火力発電ユニットを例にして、蒸気タービン製造工場およびモデルは東方蒸気タービン工場／Ｎ３００－１６．７／５３７／５３７であり、蒸気タービンバルブノズル配置は図１に示す通りで、合計４つの高圧制御バルブがあり、シーケンスバルブモードで開き順序はＧＶ１、ＧＶ２が同時に開き→ＧＶ３→ＧＶ４である。
ＤＥＨバルブ流量特性に基づいて単一バルブおよびシーケンスバルブパラメータを補正する方法は、以下のステップを含み：
（一）現場試験を行い、ある３００ＭＷ石炭火力発電ユニットを例にして、蒸気タービン製造工場およびモデルは東方蒸気タービン工場／Ｎ３００－１６．７／５３７／５３７であり、蒸気タービンバルブノズル配置は図１に示す通りで、合計４つの高圧制御バルブがあり、シーケンスバルブモードで開き順序はＧＶ１、ＧＶ２が同時に開き→ＧＶ３→ＧＶ４である。
【0037】
試験条件下でバルブ全開時負荷は基本的に約２５０ＭＷで安定し、ユニットをＤＥＨバルブ制御モードに切り替え、ユニットの一次周波数変調制御に切り替え、ユニットのＡＧＣ制御に切り替え、機械と炉の協調制御を手動に切り替え、ユニット燃料自動を切り替え、吸気酸素自動を切り替え、ユニット主蒸気温度自動、給水自動、炉圧自動を入れ、単一高圧制御バルブ１００～０％および０～１００％の単一ステップ５％のステップ外乱を行う。
【0038】
収集した結果、ＰＴ１＝［１１．８９９ １１．８３４３ １１．８０２５ １１．７８１６ １１．７３９ １１．７２８９ １１．７０８３ １１．６８７７ １１．６８７７ １１．６７０１ １１．６７０１ １１．６７０１ １１．６７０１ １１．６８１４ １１．７０３３ １１．７５５８ １１．８０１８ １１．８８０５ １１．９９９７ １２．０７９９ １２．２５４３、１２．３１６５ １２．３９８５ １２．４２８９ １２．４６１７ １２．４９２９ １２．５０３４ １２．５５１８ １２．５５１８ １２．５５１８ １２．５５１８ １２．５５１８ １２．５５１８ １２．５５１８ １２．５３１７ １２．５２１６ １２．５０１ １２．４８ １２．４６８４ １２．４５７８ １２．３６９９ １２．２５４３］。
【0039】
Ｐ１１＝［９．９０８６ ９．８５４４ ９．８２０７ ９．７９８７ ９．７７７５ ９．７５６９ ９．７３２ ９．７１９６ ９．７１９６ ９．７０１３ ９．６９１ ９．６９１ ９．６８ ９．６６６１ ９．６４３４ ９．６０７５ ９．５２５５ ９．３９６６ ９．２２０１ ９．２０４６ ９．３２２５、１０．２５１２ １０．３１４９ １０．３５３ １０．３７５ １０．４１３１ １０．４１３１ １０．４４６８ １０．４４６８ １０．４４６８ １０．４４６８ １０．４４６８ １０．４３２１ １０．４１６８ １０．３８０９ １０．３２８９ １０．２２４９ １０．０８０６ ９．８５７２ ９．５７２３ ９．４２７３ ９．３２２５］。
【0040】
ＴＳ１＝［５３３．３６５３ ５３３．１０７８ ５３３．３２２４ ５３３．３２２４ ５３３．９２３ ５３４．１８０５ ５３４．９０９９ ５３５．６３９３ ５３６．４１１６ ５３７．２２６９ ５３７．５２７２ ５３８．１７０９ ５３８．５１４２ ５３８．９８６１ ５３９．２４３６ ５３９．２８６５ ５３９．２８６５ ５３９．２８６５ ５３９．５０１ ５３９．７１５６５３９．９３０２、５２９．８４７ ５２６．８４３８ ５２６．５８６３ ５２６．３７１８ ５２６．３７１８ ５２６．８４３８ ５２７．６１６ ５２８．３０２５ ５２９．５０３８ ５３１．１３４２ ５３２．４２１３ ５３３．９６５９ ５３５．５９６４ ５３６．８８３６ ５３８．２１３８ ５３９．０２９１ ５３９．５０１ ５３９．９７３１ ５３９．９７３１ ５３９．９３０２ ５３９．９３０２］。
【0041】
Ｔ１１＝［５１０．０２４６ ５１０．０２４６ ５１０．２８２ ５１０．２８２ ５１０．６２５３ ５１１．２２６ ５１２．０８４２ ５１２．７２７８ ５１３．３２８６ ５１３．８８６４ ５１４．４４４２ ５１４．７８７５ ５１５．１３０７ ５１５．１３０７ ５１５．１３０７ ５１４．６１５８ ５１３．４１４４ ５１１．３１１９ ５０７．８７９ ５０６．５０５７ ５０５．８１９１、５０６．８０６１ ５０３．７５９１ ５０３．４５８７ ５０２．９８６６ ５０２．９８６６ ５０３．２４４１ ５０３．７５９１ ５０４．３６ ５０５．２６１２ ５０６．３７７ ５０７．８３６１ ５０８．９５１８ ５１０．２３９１ ５１１．５６９３ ５１２．１７ ５１２．２１３ ５１１．３５４８ ５０９．３８０９ ５０６．８０６１ ５０５．８１９１ ５０５．８１９１］。
【0042】
以上は、１番目の高圧制御バルブを例にして収集されたデータであり、他の３つの高圧制御バルブが収集したデータと同様であり、繰り返さない。
【0043】
シーケンスバルブ開きの反対方向に従って閉じるとき、ユニット負荷変化は、２５３．８６ＭＷ→２３７．３５ＭＷ→１８７．８５ＭＷである。
【0044】
（二）式から計算した結果、ＦＲ１ｊ％＝［１．０００ １．０００ １．０００ １．０００ １．０００ １．０００ １．０００ １．０００ １．０００ １．０００ １．０００ ０．９８６ ０．９６８ ０．９４２ ０．８９０ ０．７９４ ０．６５０ ０．４１６ ０．１００ ０．０１５ ０．０００］、ここで、ｊは（１００％、５％、０）、すなわち、１００から始まり、５間隔で０まで減少すると、以下のような２１個の数字がある。ＦＲ２ｊ％＝［１．０００ １．０００ １．０００ １．０００ １．０００ １．０００ １．０００ １．０００ １．０００ １．０００ １．０００ ０．９９５ ０．９８４ ０．９４５ ０．８８８ ０．７９４ ０．６３９ ０．４１３ ０．１０５ ０．０２３ ０．０００］。ＦＲ３ｊ％＝［１．０００ １．０００ １．０００ １．０００ １．０００ １．０００ １．０００ １．０００ １．０００ １．０００ １．０００ ０．９９３ ０．９８３ ０．９５９ ０．９１８ ０．８４９ ０．６９７ ０．４８８ ０．２１３ ０．０３２ ０．０００］。ＦＲ４ｊ％＝［１．０００ １．０００ １．０００ １．０００ １．０００ １．０００ １．０００ １．０００ １．０００ １．０００ ０．９９３ ０．９８４ ０．９６２ ０．９２６ ０．８６０ ０．７５６ ０．５８６ ０．３６０ ０．０９４ ０．０３２ ０．０００］。
【0045】
（三）単一バルブ制御パラメータを補正する。
【0046】
式から計算した結果、ＦＲ単一バルブｊ％＝［１．０００ １．０００ １．０００ １．０００ １．０００ １．０００ １．０００ １．０００ １．０００ １．０００ ０．９９８ ０．９９０ ０．９７４ ０．９４３ ０．８８９ ０．７９８ ０．６４３ ０．４１９ ０．１２８ ０．０２６ ０．０００］。
【0047】
ＤＣＳから見つけた元の単一バルブ正規化制御パラメータは表１に示される。
ｍａｔｌａｂフィッティング関数により、例えば５次の多項式関数を選択してフィッティングし、得られた関数関係は、
ｙ＝－２８５．１２８２ｘ５＋４８８．５３１５ｘ４－３０３．２１９１ｘ３＋７５．８５７９ｘ２－３．３８６２ｘ＋０．００７９７２である。
ｙ＝－２８５．１２８２ｘ５＋４８８．５３１５ｘ４－３０３．２１９１ｘ３＋７５．８５７９ｘ２－３．３８６２ｘ＋０．００７９７２である。
【0048】
例えば、バルブ開度が０．２２８である場合、補正後の総合バルブ位置指令は０．７３であるべきであり、または総合バルブ位置指令が０．８２９である場合、補正後のバルブ開度は０．３７８であるべきであり、または機能的な関係を見つけた後、カスタマイズに従ってポリゴン関数全体を変更してもよい。
【0049】
上記の原則に従って、必要に応じて単一バルブ制御パラメータを補正することができる。
【0050】
（四）シーケンスバルブ制御パラメータを補正する。
【0051】
試験データによると、２つの変曲点があり、１番目の変曲点はｙ１＝２３７．３５／２５３．８６＝０．９３４９に対応し、２番目の変曲点はｙ２＝１８７．３５／２５３．８６＝０．７３８に対応する。
【0052】
ＤＣＳで見つけた現在で使用されるシーケンスバルブ制御正規化パラメータは、表２に示される。
表２において、実際用の変曲点は、１番目は０．８９９、２番目は０．７１２であり、現場試験において測定された１番目の変曲点が０．９３４９、２番目が０．７３８である状況に近く、オーバーラップの要件も満たしているため、発電所のパラメータの大半の使用に合わせて引き続き使用することができ、本明細書で計算したものを直接使用してもよく、すなわち、表２中のＧＶ３総合バルブ位置指令欄の０．７１２を０．７３８に置き換え、ＧＶ４中の０．８９９を０．９３４９に置き換えることもできる。
【0053】
ステップ（二）で得られた離散データをｍａｔｌａｂでフィッティングした結果、
ｙ１＝－３２３．０７６９ｘ１ ５＋５５０．７９２５ｘ１ ４－３３９．４２３１ｘ１ ３＋８４．６５０１ｘ１ ２－４．１４１２ｘ１＋０．００９５５５９
ｙ２＝－３０１．０２５６ｘ２ ５＋５１０．１７４８ｘ２ ４－３１５．１２８８ｘ２ ３＋７９．１２３６ｘ２ ２－３．７５３ｘ２＋０．００９５３３２
ｙ３＝－３９３．８４６２ｘ３ ５＋６２１．９５８ｘ３ ４－３５４．９３２４ｘ３ ３＋８０．９２１１ｘ３ ２－２．８６２ｘ３＋０．００２５２４５
ｙ４＝－１２３．５８９７ｘ４ ５＋２７３．８３４５ｘ４ ４－２０５．２５５８ｘ４ ３＋５９．２１０５ｘ４ ２－２．８２２７ｘ４＋０．００９９７７３
上記の表および変曲点データを参照して、各バルブ開度を上記公式に代入し、対応の相対流量を計算して得、さらに現在で使用されるシーケンスバルブパラメータ下で相対流量と総合バルブ位置指令間の対照表は表３に示される。
ｙ１＝－３２３．０７６９ｘ１ ５＋５５０．７９２５ｘ１ ４－３３９．４２３１ｘ１ ３＋８４．６５０１ｘ１ ２－４．１４１２ｘ１＋０．００９５５５９
ｙ２＝－３０１．０２５６ｘ２ ５＋５１０．１７４８ｘ２ ４－３１５．１２８８ｘ２ ３＋７９．１２３６ｘ２ ２－３．７５３ｘ２＋０．００９５３３２
ｙ３＝－３９３．８４６２ｘ３ ５＋６２１．９５８ｘ３ ４－３５４．９３２４ｘ３ ３＋８０．９２１１ｘ３ ２－２．８６２ｘ３＋０．００２５２４５
ｙ４＝－１２３．５８９７ｘ４ ５＋２７３．８３４５ｘ４ ４－２０５．２５５８ｘ４ ３＋５９．２１０５ｘ４ ２－２．８２２７ｘ４＋０．００９９７７３
上記の表および変曲点データを参照して、各バルブ開度を上記公式に代入し、対応の相対流量を計算して得、さらに現在で使用されるシーケンスバルブパラメータ下で相対流量と総合バルブ位置指令間の対照表は表３に示される。
【0054】
具体的な計算過程は以下のとおりであり：
１、２つの変曲点の総合バルブ位置指令および相対流量を決定し、すなわち総合バルブ位置指令が０．７１２である場合、相対流量が０．７３８であり、総合バルブ位置指令が０．７１２である場合、相対流量が０．９３４９である。
１、２つの変曲点の総合バルブ位置指令および相対流量を決定し、すなわち総合バルブ位置指令が０．７１２である場合、相対流量が０．７３８であり、総合バルブ位置指令が０．７１２である場合、相対流量が０．９３４９である。
【0055】
２、総合バルブ位置が０．７１２以下である場合、総合バルブ位置指令が０．５８５である場合を例にすると、元のパラメータ設定においてＧＶ１、ＧＶ２の開度がすべて０．３１９であり、それぞれ上記関数関係に代入して求めた２つのｙ値はそれぞれ０．９２０７および０．９２３であり、２つの値の平均値に０．７３８を乗じて計算した結果、このパラメータ下での実際の相対流量が０．６８であり、他の値は計算値を参照すればよい。
【0056】
３、総合バルブ位置が０．７１２を超え、０．８９９以下である場合、総合バルブ位置が０．８３５である場合を例にすると、元のパラメータ設定においてＧＶ３の開度が０．２０６であり、上記関数関係に代入して求めたｙ値は０．７１８であり、この場合、ＧＶ１、ＧＶ２が全開であるため、該値０．７１８に（０．９３４９－０．７３８）を乗じた後０．７３８を加算して計算した結果、このパラメータ下での実際の相対流量が０．８７９であり、他の値は計算値を参照すればよい。
【0057】
４、総合バルブ位置が０．８９９を超える場合、総合バルブ位置が０．９６０である場合を例にすると、元のパラメータ設定においてＧＶ４の開度が０．２０６であり、上記関数関係に代入して求めたｙ値は０．５９４であり、この場合、ＧＶ１、ＧＶ２、ＧＶ３が全開であるため、該値０．５９４に（１－０．９３４９）を乗じた後０．９３４９を加算して計算した結果、このパラメータ下での実際の相対流量は０．９７４であり、他の値は計算値を参照すればよい。
理論的に正規化後の総合バルブ位置指令と相対流量は
の関係を有し、理論的な関数と実際の曲線を同じグラフに描き、図２を参照する。
【0058】
例えば、比較した結果、総合バルブ位置が０．５８５である場合に対応のバルブ開度を補正する必要があることが判明され、表３から分かるように、現在で使用されるパラメータ計算で得られた相対流量は０．６８であり、理論的に相対流量と総合バルブ位置指令は正比例関係を有するべきであるため、理論値は０．５８５であり、補正する必要があり、表２から分かるように、総合バルブ位置指令が０．５８５である場合、動作するのが１、２番目の高圧制御バルブだけであるため、１、２番目の高圧制御バルブを補正する必要があり、逆関数法を利用して、１、２番目の高圧制御バルブ開度は０．３１９から０．２５に補正する必要がある。他の総合バルブ位置指令に対応するバルブ開度の補正方法は同様である。
【0059】
厳密に言えば、すべてのポイントが理論曲線上、すなわち正比例関係曲線にない限り補正することができるが、現場で一般的に最小範囲にしたがって補正し、すなわち理論曲線から顕著に偏差しているポイントのみ補正しているため、実施例では偏差が最も大いいポイントを例にして説明する。
【0060】
これまで、シーケンスバルブ制御パラメータ補正の方法を説明した。
【0061】
以上、本発明の基本原理、主要特徴および本発明の利点を示し、説明した。当業者であれば、本発明が上記実施例によって制限されなく、上記実施例および明細書の説明は本発明の原理を例示するものに過ぎず、本発明の精神および範囲から逸脱することなく、本発明に対して様々な変更や改良を加えることが可能であり、これらの変更および改良はすべて保護しようとする本発明の範囲内に含まれることを理解されたい。本発明が保護しようとする範囲は、添付の特許請求の範囲およびその等価物によって定義されるものとする。
【図1】
【図2】
【図3】
【国際調査報告】