JP4577232B2 - Power plant operation plan creation device and method - Google Patents

Power plant operation plan creation device and method Download PDF

Info

Publication number
JP4577232B2
JP4577232B2 JP2006044682A JP2006044682A JP4577232B2 JP 4577232 B2 JP4577232 B2 JP 4577232B2 JP 2006044682 A JP2006044682 A JP 2006044682A JP 2006044682 A JP2006044682 A JP 2006044682A JP 4577232 B2 JP4577232 B2 JP 4577232B2
Authority
JP
Japan
Prior art keywords
demand
stop
power plant
objective function
plan
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP2006044682A
Other languages
Japanese (ja)
Other versions
JP2007228676A (en
Inventor
敏之 澤
満男 鶴貝
司 大西
康生 佐藤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Ltd
Original Assignee
Hitachi Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hitachi Ltd filed Critical Hitachi Ltd
Priority to JP2006044682A priority Critical patent/JP4577232B2/en
Publication of JP2007228676A publication Critical patent/JP2007228676A/en
Application granted granted Critical
Publication of JP4577232B2 publication Critical patent/JP4577232B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • Y04S10/54

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Description

本発明は火力発電所,揚水発電所等の発電所運用計画を作成する上で、組合せ問題となる火力発電所の起動停止計画を作成することにより、需給バランス,予備力,燃料消費量制約,ダム水位制約,潮流制約等の制約条件を満たし、発電コストを最小化する発電所の運用計画を作成する方法および装置に関する。   The present invention creates a power plant operation plan such as a thermal power plant, a pumped storage power plant, etc., and creates a start / stop plan for a thermal power plant that becomes a combination problem, thereby providing a balance between supply and demand, reserve power, fuel consumption constraints, The present invention relates to a method and apparatus for creating a power plant operation plan that satisfies constraints such as dam water level constraints and tidal current constraints and minimizes power generation costs.

電力系統の信頼度を維持するためには、電力需要と発電機による供給力を一致させながら、予備力を確保した発電機の起動停止計画を作成する必要がある。しかしながら、計画段階では予測した電力需要を使用するため、予測誤差を伴う。信頼性と同時に経済的な運用を行うためには、この予測誤差を考慮して、発電コストが最小となるように、発電機の起動停止計画を作成する必要がある。発電所運用計画問題は、火力発電所の起動停止計画は組合せ問題となっており、また発電機出力は非線形計画問題となっている。このように、非線形混合整数計画問題であると同時に、発電機台数および計画時点数を考慮すると問題の規模が非常に大きいため、事実上最適化することは難しい問題となっている。このため、特開2004−274956号公報のようにいろいろな需要パターンシナリオのそれぞれに対応する計画を作成して、シナリオ毎のコストの分布を計算している。また、特開2005−12912号公報では同一特性の発電機に対して、需要誤差を考慮して運転台数を決定している。   In order to maintain the reliability of the power system, it is necessary to create a start-up / shutdown plan for the generator that ensures reserve power while matching the power demand with the supply capacity of the generator. However, since the predicted power demand is used in the planning stage, it involves a prediction error. In order to perform economic operation at the same time as reliability, it is necessary to create a start-up / stop plan for the generator so that the power generation cost is minimized in consideration of this prediction error. As for the power plant operation plan problem, the thermal power plant start / stop plan is a combination problem, and the generator output is a nonlinear programming problem. As described above, it is a nonlinear mixed integer programming problem, and at the same time, considering the number of generators and the number of planning points, the problem is very large, and therefore it is difficult to optimize in practice. For this reason, a plan corresponding to each of various demand pattern scenarios is created as in JP-A-2004-274956, and the cost distribution for each scenario is calculated. In Japanese Patent Laid-Open No. 2005-12912, the number of operating units is determined in consideration of a demand error for generators having the same characteristics.

特開2004−274956号公報Japanese Patent Laid-Open No. 2004-274958 特開2005−12912号公報JP 2005-12912 A

火力発電所の起動停止計画は組合せ最適化問題であるため、例えば、火力発電機10台で24時間の起動停止の全組合せ数は10の72乗にもなる。1つの計画を1マイクロ秒で評価できたとしても全ての計画を評価するためには10の59乗の年数を要することとなり、事実上全組合せをチェックして、最適解を求めることは不可能である。このため、この火力発電所の起動停止計画を作成すると同時に、火力発電所,揚水発電所の最適な発電機出力を決定することも事実上不可能である。そこで、火力発電機の起動停止計画問題を組合せ問題から実数変数の問題へと緩和して、二次計画法,線形計画法を用いて解くことにより、組合せ爆発の問題を回避できる。このように実数問題に緩和することにより、燃料消費制約,電力潮流制約、あるいはダム水位制約等の線形制約で表現できる制約を同時に考慮して、最適解を求めることができる。但し、火力発電機の起動停止変数が起動
(1)または停止(0)に確定できない未確定の状態(例えば、0.3 )のものも発生するという課題が発生する。このため、繰り返し収束計算の途中で得られる起動停止変数の値や発電機出力が大きいものほど起動しやすく、小さいものほど停止しやすくなるように、目的関数,制約条件を修正することにより、起動停止変数の値を起動、または停止に確定していくことができる。更に、使用する電力需要は予測値を用いるため、実際の需要とは誤差が存在する。この誤差は過去の需要予測結果と需要実績からがどれくらいの誤差分布かを推定することができる。この需要誤差は必ず発生するものであるから、これを考慮して、発電所の運用計画を作成することが必要となる。複数の需要に対する問題に定式化して、二次計画法あるいは線形計画法を用いることにより、火力発電所の起動停止変数の値を確定すると同時に最適な火力発電機,揚水発電機に出力が決定できる。
Since the thermal power plant startup / shutdown plan is a combination optimization problem, for example, the total number of startup / shutdowns for 24 hours with 10 thermal power generators is 10 to the 72nd power. Even if one plan can be evaluated in 1 microsecond, it takes 10 59 years to evaluate all plans, and it is impossible to find the optimal solution by checking all combinations in practice. It is. For this reason, it is virtually impossible to determine the optimum generator output of the thermal power plant and the pumped-storage power plant at the same time as creating the startup / shutdown plan of this thermal power plant. Therefore, the problem of combined explosion can be avoided by relaxing the start / stop planning problem of a thermal power generator from a combination problem to a problem of real variables and solving it using a quadratic programming method and a linear programming method. By relaxing to the real number problem in this way, it is possible to obtain an optimal solution while simultaneously taking into consideration constraints that can be expressed by linear constraints such as fuel consumption constraints, power flow constraints, or dam water level constraints. However, there arises a problem that an uncertain state (for example, 0.3) in which the start / stop variable of the thermal power generator cannot be determined to start (1) or stop (0) occurs. For this reason, it is possible to start by correcting the objective function and constraint conditions so that the larger the start / stop variable value and generator output obtained during the repeated convergence calculation, the easier it is to start and the smaller it is, the easier it is to stop. The value of the stop variable can be confirmed to start or stop. Furthermore, since the power demand to be used uses a predicted value, there is an error from the actual demand. This error can estimate how much error distribution is based on past demand prediction results and actual demand. Since this demand error always occurs, it is necessary to prepare an operation plan for the power plant taking this into account. By formulating a problem for multiple demands and using quadratic programming or linear programming, it is possible to determine the values of startup and shutdown variables of a thermal power plant and at the same time determine the output to the optimal thermal generator or pumped-storage generator .

特開2004−274956号公報では、需要ごとに計画を作成しているため、平均的な需要に対する平均的なコストを求めていることになる。しかし、需要予測誤差を考慮すると、需要が大きいほうが発電コストも大きくなり、平均的な発電コストよりも発電コストの期待値は大きい方にずれるため、需要予測誤差を適切に評価して計画を作成していない。また、特開2005−12912号公報では同一特性の発電機を対象としているが、電力系統で運用されている火力発電機の特性は需要変動に適切に対応できるように、例えばベース火力,ミドル火力,ピーク火力のように異なると特性の発電機が導入されているため、同一特性の発電機を台数制御する方法では発電コストを最小化することはできない。   In Japanese Patent Application Laid-Open No. 2004-274956, since a plan is created for each demand, an average cost for the average demand is obtained. However, considering the demand forecast error, the larger the demand, the larger the power generation cost, and the expected value of the power generation cost shifts to the larger one than the average power generation cost. Not done. Japanese Patent Application Laid-Open No. 2005-12912 targets generators having the same characteristics, but the characteristics of thermal power generators operated in the power system can be appropriately adapted to demand fluctuations, for example, base thermal power, middle thermal power , Since generators with different characteristics such as peak thermal power are introduced, the method of controlling the number of generators with the same characteristics cannot minimize power generation costs.

本発明は上記点に対処して成されたもので、その目的とするところは発電機の運用制約,発電機の特性,燃料消費制約,揚水ダム水位制約のみならず、電力需要予測誤差を考慮して発電費用を最小化できる発電所運用計画作成方法および装置を提供することにある。   The present invention has been made in response to the above points, and its objective is to consider not only generator operation constraints, generator characteristics, fuel consumption constraints, pumping dam water level constraints, but also power demand prediction errors. It is an object of the present invention to provide a power plant operation plan creation method and apparatus that can minimize power generation costs.

本発明の主な特徴によれば、予測誤差に基づいた複数の需要カーブを入力して、火力発電機の起動停止変数を実数変数に緩和することにより、1つの起動停止計画問題と需要毎の負荷配分計画問題を実数変数問題として二次計画法あるいは線形計画法を適用し、燃料消費制約,揚水ダム水位制約等を考慮して最適な発電所運用計画作成方法および装置により発電コスト最小となる起動停止計画と負荷配分計画を作成できる。   According to the main feature of the present invention, by inputting a plurality of demand curves based on the prediction error and relaxing the start / stop variable of the thermal power generator to a real variable, one start / stop plan problem and each demand Applying quadratic programming or linear programming to load distribution planning problems as real variable problems, and taking into account fuel consumption constraints, pumping dam water level constraints, etc., and optimal power plant operation plan creation method and equipment minimize power generation costs Start / stop plans and load distribution plans can be created.

本発明によれば、電力需要誤差が発生しても、燃料消費制約,揚水ダム水位制約および発電機の運用制約等を考慮して、発電費用を少なくすることができる起動停止計画と負荷配分計画を作成することができる。   According to the present invention, even if a power demand error occurs, a start / stop plan and a load distribution plan that can reduce power generation costs in consideration of fuel consumption constraints, pumping dam water level constraints, generator operation constraints, and the like. Can be created.

以下、本発明の実施例について図面を参照しながら詳細に説明する。   Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings.

図1は、本発明の第一の実施例の適用対象である発電所運用計画方法の計算処理フロー、図2は発電所需給運用計画作成装置1の概略全体構成図、図3は運用計画装置の詳細構成図を示す。   FIG. 1 is a calculation processing flow of a power plant operation planning method to which the first embodiment of the present invention is applied, FIG. 2 is a schematic overall configuration diagram of a power plant supply and demand operation plan creation device 1, and FIG. The detailed block diagram of an apparatus is shown.

発電所需給運用計画作成装置1は中央演算処理装置(CPU)2,主記憶装置3,入出力装置4および外部記憶装置5とから構成される。   The power plant supply and demand operation plan creation device 1 includes a central processing unit (CPU) 2, a main storage device 3, an input / output device 4 and an external storage device 5.

図3は本発明における発電所需給運用計画作成装置の詳細な機能構成図である。発電所需給運用計画作成装置1は、上述の入力装置6,表示装置7,データベース9および読取装置10と、運用計画作成処理部20とを備える。なお、運用計画作成処理部20は、予め外部記憶装置5の記憶媒体に保持され、読取装置10を介して主記憶装置3に読み込まれたプログラムを、CPU2が実行することにより実現されるが、本発明はこのようなプログラムされた汎用プロセッサによるものに限られるわけではない。例えば、本発明の各処理を実行するワイヤードロジックを含む特定のハードウェア装置との組合せによって運用計画作成処理部20を構成することもできる。   FIG. 3 is a detailed functional configuration diagram of the power plant supply and demand operation plan creation apparatus according to the present invention. The power plant supply and demand operation plan creation device 1 includes the input device 6, the display device 7, the database 9, the reading device 10, and the operation plan creation processing unit 20. The operation plan creation processing unit 20 is realized by the CPU 2 executing a program that is stored in advance in the storage medium of the external storage device 5 and read into the main storage device 3 via the reading device 10. The present invention is not limited to such a programmed general purpose processor. For example, the operation plan creation processing unit 20 can be configured by a combination with a specific hardware device including a wired logic that executes each processing of the present invention.

入出力装置4は、図3に示すキーボードやマウスを備えた入力装置6と、出力装置としての表示装置7を有する。なお、入出力装置4としては、これらの代わりに、あるいは、これらと併用して、ポインティングデバイス,タッチセンサ等の入力装置や、液晶表示装置,プリンタ,スピーカ等の出力装置を設けることもできる。外部記憶装置5としては、ハードディスク装置,FD装置,CD−ROM(compact disc-read onlymemory)装置,DAT(digital videotape)装置,RAM(randam access memory) 装置,DVD
(digital video disc)装置,不揮発性メモリ等を用いることができる。外部記憶装置5は、図3に示すデータベース9を保持するための大容量記憶装置と、処理プログラムなどを保持する記憶媒体と、該記憶媒体に保持された情報を読み取るための読取装置10とを用いているが、一つの外部記憶装置にデータベースと処理プログラムとが両方保持されているようにすることもできる。また、記憶媒体としては、FD,CD−ROM,磁気テープ,光ディスク,光磁気ディスク,DAT,RAM,DVD,不揮発性メモリ等を用いることができる。
The input / output device 4 includes an input device 6 having a keyboard and a mouse shown in FIG. 3 and a display device 7 as an output device. The input / output device 4 may be provided with an input device such as a pointing device or a touch sensor, or an output device such as a liquid crystal display device, a printer, or a speaker instead of or in combination with these. As the external storage device 5, a hard disk device, an FD device, a CD-ROM (compact disc-read only memory) device, a DAT (digital videotape) device, a RAM (randam access memory) device, a DVD
A (digital video disc) device, a non-volatile memory, or the like can be used. The external storage device 5 includes a mass storage device for holding the database 9 shown in FIG. 3, a storage medium for holding a processing program, and a reading device 10 for reading information held in the storage medium. Although used, both the database and the processing program can be held in one external storage device. As the storage medium, FD, CD-ROM, magnetic tape, optical disk, magneto-optical disk, DAT, RAM, DVD, nonvolatile memory, and the like can be used.

入力装置6は表示装置7に表示された選択肢の選択,データの入力などを受付け、運用計画作成処理部20に伝送する。表示装置7は、入力装置6から送られたデータを表示する。運用計画作成処理部20は、入力装置6から伝送されるデータと、データベース9から読み込んだデータと、読取装置10から読み込んだ処理プログラムと、電力システム
11から伝送されるデータとに基づいて発電所需給運転計画を作成する。電力システム
11は図示していないデータベースを備えている。電力システム11は電力系統12の監視,制御を行うものである。
The input device 6 accepts selection of options displayed on the display device 7, input of data, and the like, and transmits them to the operation plan creation processing unit 20. The display device 7 displays the data sent from the input device 6. The operation plan creation processing unit 20 generates the power plant based on the data transmitted from the input device 6, the data read from the database 9, the processing program read from the reading device 10, and the data transmitted from the power system 11. Create a supply and demand plan. The power system 11 includes a database (not shown). The power system 11 monitors and controls the power system 12.

運用計画作成処理部20の処理結果は、表示装置7に送られて表示されると共にデータベース9に格納される。また、電力システム11の運用・計画システム,制御システムが設定した条件での運用計画作成が要求された場合には、運用計画作成処理部20の作成した計画および評価結果は電力システム11へも通知される。この通知を受けた電力システム11は、通知された計画をもとに電力系統12に供給力のひとつである発電機の出力制御信号を出して発電機を制御するとともに、電力系統から発電機出力データを取り込み、その取り込んだデータと図示していない計画システムにおいて作成したデータとを内部データベース(図示せず)に格納する。   The processing result of the operation plan creation processing unit 20 is sent to the display device 7 for display and stored in the database 9. In addition, when an operation plan creation under the conditions set by the operation / planning system and control system of the power system 11 is requested, the plan and evaluation result created by the operation plan creation processing unit 20 are also notified to the power system 11. Is done. Upon receiving this notification, the power system 11 outputs an output control signal of a generator, which is one of the supply powers, to the power system 12 based on the notified plan, and controls the generator. Data is fetched, and the fetched data and data created in a planning system (not shown) are stored in an internal database (not shown).

運用計画作成処理部20は、コントロール部21,計画作成条件設定部22,データ読込部23,定式化部24,計画計算部25,起動停止変数上下限修正部26,目的関数修正部27および計算結果処理部28を備える。また、運用計画作成処理部20は、通信線31を介して外部の電力システム11に接続されており、電力システム11は、実際の送発変電設備である電力系統12に接続されている。電力システム11は電力系統12の計画・運用及び制御のためのシステム、電力系統12および電力系統12の状態などを示す情報を保持するための図示していないデータベースを備えている。電力系統12の状態はリレーやセンサなどにより検出され、通信線32を介して電力システム11に通知され、電力システム11のデータベース(図示せず)に格納される。   The operation plan creation processing unit 20 includes a control unit 21, a plan creation condition setting unit 22, a data reading unit 23, a formulation unit 24, a plan calculation unit 25, a start / stop variable upper / lower limit modification unit 26, an objective function modification unit 27, and a calculation. A result processing unit 28 is provided. In addition, the operation plan creation processing unit 20 is connected to an external power system 11 via a communication line 31, and the power system 11 is connected to a power system 12 that is an actual transmission and substation facility. The power system 11 includes a system for planning / operation and control of the power system 12, a power system 12 and a database (not shown) for holding information indicating the state of the power system 12, and the like. The state of the power system 12 is detected by a relay, a sensor, or the like, notified to the power system 11 via the communication line 32, and stored in a database (not shown) of the power system 11.

コントロール部21は、電力システム11、および、各処理部22〜28の間のデータや処理プログラムなどの授受を円滑に行うためのデータの加工・処理を行い、その授受をコントロールして、全体の処理を正常に動作させる。   The control unit 21 processes and processes data for smoothly transmitting and receiving data and processing programs between the power system 11 and each of the processing units 22 to 28, and controls the transmission and reception of the entire system. Make the process work properly.

計画作成条件設定部22はコントロール部21を介してデータベース9および/または読取装置10の記憶媒体に保持されている、運用計画を作成するためのデフォルトの条件等を読み込んで表示するとともに、これらの条件を表示装置7に表示して、条件を変更する。また、この変更した新規の条件をコントロール部21を介してデータベース9に格納する。   The plan creation condition setting unit 22 reads and displays the default conditions for creating the operation plan, which are held in the storage medium of the database 9 and / or the reading device 10 via the control unit 21, and displays these conditions. The condition is displayed on the display device 7 to change the condition. Further, the changed new condition is stored in the database 9 via the control unit 21.

データ読込部23は、コントロール部21を介してデータベース9に保持されているデータを読み込み、これらのデータをコントロール部21を介して各処理部24〜28に伝送する。なお、本実施例では、特に説明しない限り、データベース9からのデータの読み込みはデータ読込部23を介して行われる。   The data reading unit 23 reads data held in the database 9 via the control unit 21 and transmits these data to the processing units 24 to 28 via the control unit 21. In this embodiment, data is read from the database 9 through the data reading unit 23 unless otherwise specified.

定式化部24は、コントロール部21を介して計画作成条件設定部22で作成された条件やデータ読込部23からのデータをもとに計画計算部25で計算するための前処理および定式化を実施する。   The formulation unit 24 performs preprocessing and formulation for calculation by the plan calculation unit 25 based on the conditions created by the plan creation condition setting unit 22 and the data from the data reading unit 23 via the control unit 21. carry out.

計画計算部25はコントロール部21を介して計画作成条件設定部22で作成された条件あるいは起動停止変数上下限修正部26,目的関数修正部27で修正された制約条件,目的関数やデータ読込部23からのデータあるいは定式化部24で作成された前処理および定式化により発電所運用計画の計算処理を実施する。   The plan calculation unit 25 is a condition created by the plan creation condition setting unit 22 via the control unit 21, or the constraint condition, objective function and data reading unit modified by the start / stop variable upper / lower limit modification unit 26 and the objective function modification unit 27. The calculation processing of the power plant operation plan is executed by the data from 23 or the pre-processing and formulation created by the formulation unit 24.

起動停止変数上下限修正部26はコントロール部21を介して計画計算部25で作成した運用計画をもとに繰り返し計算して、新たな計画を作成させるときの起動停止変数の上下限制約値を修正する。   The start / stop variable upper / lower limit correction unit 26 repeatedly calculates the upper / lower limit constraint values for the start / stop variable when creating a new plan by repeatedly calculating based on the operation plan created by the plan calculation unit 25 via the control unit 21. Correct it.

目的関数修正部27はコントロール部21を介して計画計算部25で作成した運用計画をもとに繰り返し計算して、新たな計画を作成させるときの目的関数を修正,設定する。   The objective function correction unit 27 repeatedly calculates based on the operation plan created by the plan calculation unit 25 via the control unit 21, and corrects and sets the objective function when creating a new plan.

計算結果処理部28は、各処理部22〜27で設定,計算した計画作成条件,データ,作成した運用計画や目的関数の推移等や入力装置6からのデータ入力の支援のための情報等を、表示装置7に表示する。   The calculation result processing unit 28 stores the plan creation conditions and data set and calculated by each of the processing units 22 to 27, the created operation plan, the transition of the objective function, information for supporting data input from the input device 6, and the like. Is displayed on the display device 7.

次に、動作について図1に示す処理フロー図を参照して説明する。   Next, the operation will be described with reference to the processing flowchart shown in FIG.

処理S101の計画作成条件の設定では、計画作成条件設定部22によりデフォルトの計画再生条件をデータベース9から読み込んで、表示装置7に表示するとともに、入力装置6を介して条件を変更する。ここで設定する条件には、複数の電力需要予測カーブとそのカーブの重み係数,需要増加時間帯,需要減少時間帯,繰り返し計算回数がある。図4にはこれらの値を設定する画面の例を示している。また、設定した値をグラフで表示するものである。他に設定するものとして、予備力,LNGなど消費制約,揚水ダムの水位制約,系統構成,送電線の潮流制約値等がある。計画作成期間はここでは1時間毎24時間の電力需要予想のデータを設定している。需要予測データは、通常の方法で予想したものをカーブ0とし、過去の予測誤差の分布をもとにカーブ1から4を設定し、その確率を重み係数として設定している。予備力は需要のカーブ0をもとに設定する。燃料消費制約については、計画期間の最終時刻までに消費する量や必要に応じて各時刻における燃料消費量の上下限を設定する。揚水ダム水位についても、初期水位と最終時刻の水位や各時刻の水位上下限を設定する。系統構成についても各時刻の送電線と発電設備,負荷との接続状態を設定し、各時刻の送電線,予備力を設定する。   In setting the plan creation conditions in step S101, the plan creation condition setting unit 22 reads default plan regeneration conditions from the database 9, displays them on the display device 7, and changes the conditions via the input device 6. The conditions set here include a plurality of power demand prediction curves and their weighting factors, a demand increase time zone, a demand decrease time zone, and the number of repeated calculations. FIG. 4 shows an example of a screen for setting these values. The set value is displayed in a graph. Other items to be set include reserve capacity, consumption restrictions such as LNG, water level restrictions for pumping dams, system configuration, power flow restriction values for transmission lines, and the like. In the plan creation period, power demand forecast data for 24 hours per hour is set here. The demand forecast data is a curve 0 that is predicted by a normal method, curves 1 to 4 are set based on the distribution of past prediction errors, and the probability is set as a weighting factor. The reserve capacity is set based on the demand curve 0. Regarding the fuel consumption restriction, the amount consumed by the final time of the planning period and the upper and lower limits of the fuel consumption amount at each time are set as necessary. For the pumping dam water level, the initial water level, the final water level, and the upper and lower water level at each time are set. Regarding the system configuration, the connection state between the transmission line at each time, the power generation equipment, and the load is set, and the transmission line and reserve capacity at each time are set.

処理S102は発電機等の供給力の特性や運用制約をデータベース9から読み込む。供給力のなかで火力発電機の特性として出力の上下限値,燃料費特性(出力に対する二次関数,線形関数),起動費や同時起動ができない発電機の組みがある。各発電機の使用する燃料種別もデータベースから読み出す。また、揚水発電所については、発電時とポンプ負荷の上下限値および揚水効率がある。   In step S102, the characteristics of supply power of the generator and the operation restrictions are read from the database 9. Among the supply power, there are thermal generator characteristics such as output upper and lower limits, fuel cost characteristics (secondary function and linear function for output), start-up costs and generator sets that cannot be started simultaneously. The fuel type used by each generator is also read from the database. For pumped storage power plants, there are upper and lower limits and pumping efficiency during power generation and pump load.

処理S103では目的関数,制約条件を用いて定式化を行う。ここでは、まず一般的な定式化、次に複数の需要を考慮したときの定式化、最後に本発明で使用する定式化の順で説明する。   In process S103, formulation is performed using an objective function and constraint conditions. Here, the general formulation will be described first, then the formulation when considering a plurality of demands, and finally the formulation used in the present invention.

一般的な定式化における目的関数TCは式(1)である。この式は、需要カーブが1つのときのものである。起動停止変数はuitで発電機i,時刻tの状態を示しており、0のとき停止、1のとき運転であることを示す。また、燃料費関数Fi(Pit) は発電機出力
itの二次関数となっている。起動費SUCi は最小需要時刻nで停止、最大需要時刻pで運転しているときに発生する。このため、最小需要時刻nと最大需要時刻pの起動停止変数のuinとuipを使用している。燃料費関数のうち定数項を除いた関数をfi(Pit) としている。
The objective function TC in the general formulation is Equation (1). This equation is for a single demand curve. The start / stop variable indicates the state of the generator i and the time t with u it. When the variable is 0, it is stopped and when it is 1, the operation is stopped. The fuel cost function F i (P it ) is a quadratic function of the generator output P it . The start-up cost SUC i is incurred when stopping at the minimum demand time n and operating at the maximum demand time p. Therefore, the start / stop variables u in and u ip of the minimum demand time n and the maximum demand time p are used. A function excluding a constant term in the fuel cost function is defined as f i (P it ).

Figure 0004577232
Figure 0004577232

式(2)は火力発電機の出力上下限,揚水発電所の揚水発電,揚水ポンプの出力上下限を示している。火力発電機の起動停止変数uitは0または1であるため、0のときは火力発電機の出力上下限はともに0となり、1のときはPi,minからPi,maxとなる。 Equation (2) shows the output upper and lower limits of the thermal power generator, the pumped power generation of the pumped storage power plant, and the output of the pump. Since the start-stop variable u it thermal generator is 0 or 1, becomes 0 both the output on the lower limit of the thermal power machine when 0, P i is set to 1, the min P i, and max.

Figure 0004577232
Figure 0004577232

式(3)は需給バランスの式である。時刻tの需要がDtである。 Formula (3) is a formula for supply and demand balance. Demand of the time t is D t.

Figure 0004577232
Figure 0004577232

式(4)は予備力制約である。時刻tで必要な予備力はRt である。ここでは、揚水発電所の予備力を考慮しているが、火力のみとすることもある。 Equation (4) is a reserve force constraint. The reserve force required at time t is R t . Here, the reserve capacity of the pumped storage power plant is considered, but it may be only thermal power.

Figure 0004577232
Figure 0004577232

式(5)は同一中央の同時起動禁止制約を示している。古い型の制御装置は、2台の発電機を起動するときに、1日に1台しか起動できないことを考慮するものである。式(5)では火力発電機iと火力発電機(i+1)が同一中央の制約があるときの例を示している。   Equation (5) shows the same central simultaneous prohibition constraint. Older control devices take into account that when starting two generators, only one can be started per day. Equation (5) shows an example when the thermal power generator i and the thermal power generator (i + 1) have the same central constraint.

Figure 0004577232
Figure 0004577232

式(6)は火力発電機でも燃料がLNGの発電機の発電量合計LNGが上下限範囲内となるようにする制約である。   Expression (6) is a constraint that the total power generation amount LNG of the generator with LNG as the fuel is within the upper and lower limit range even in the thermal power generator.

Figure 0004577232
Figure 0004577232

式(7)は揚水ダムの水位制約を示している。ηは揚水効率を示している。発電量に換算したダム水位WLtはダム水位の上下限範囲内となるようにする制約である。 Equation (7) shows the water level constraint of the pumping dam. η indicates the pumping efficiency. The dam water level WL t converted to the power generation amount is a constraint that the dam water level is within the upper and lower limits of the dam water level.

Figure 0004577232
Figure 0004577232

以上の目的関数および制約条件をもとに発電所の起動停止計画を作成し、次にこの起動停止計画をもとに負荷配分計画を作成するのが一般的である。   It is common to create a power plant start / stop plan based on the above objective function and constraints, and then create a load distribution plan based on this start / stop plan.

次に、複数の需要予想カーブDt kがあるときについて説明する。式(8)が目的関数である。各需要カーブでの発電コスト期待値を加算したものである。式(1)にその需要発生確率に相当する重み係数αk をかけたものである。このとき、重み係数αk の合計は式(9)に示すように1となる。式(8)で火力の発電機出力はPit k となっており、同一発電機、同一時刻でも需要カーブkごと出力が異なることがわかる。一方、起動停止変数はuitと需要カーブkに依存しないことがわかる。すなわち、1つの起動停止計画に対して、需要予測誤差を考慮したカーブを使用すると発電機の出力がそれぞれの需要カーブで異なり、発電コストが変化する。このため、需要カーブの発生確率を重み係数として、期待値を計算し、これが最小となる計画を作成している。 Next, a case where there are a plurality of demand forecast curves D t k will be described. Equation (8) is the objective function. This is the sum of the expected power generation costs for each demand curve. Equation (1) is multiplied by a weighting factor α k corresponding to the demand occurrence probability. At this time, the sum of the weighting factors α k is 1 as shown in equation (9). In equation (8), the thermal power generator output is P it k, and it can be seen that the output varies with the demand curve k even at the same generator and the same time. On the other hand, it can be seen that the start / stop variable does not depend on u it and the demand curve k. That is, if a curve that takes into account a demand prediction error is used for one start / stop plan, the output of the generator differs in each demand curve, and the power generation cost changes. For this reason, the expected value is calculated using the probability of occurrence of the demand curve as a weighting factor, and a plan that minimizes this is created.

Figure 0004577232
Figure 0004577232

Figure 0004577232
Figure 0004577232

式(10)に示すように、各需要カーブDt kごとに需給バランス制約が成り立つ必要がある。 As shown in equation (10), the supply-demand balance constraints need to hold each demand curve D t k.

Figure 0004577232
Figure 0004577232

式(11)は予備力制約を示している。ここで、D0 tは基準需要カーブのときのみでの予備力を考えている。ここでは、需要カーブ0を基準需要カーブとしている。 Equation (11) shows the reserve power constraint. Here, D 0 t contemplates reserve only when the reference demand curve. Here, the demand curve 0 is used as a reference demand curve.

Figure 0004577232
Figure 0004577232

式(5),式(6)および式(7)の制約は基準カーブのときの発電機出力を用いるものとする。   For the constraints of Equation (5), Equation (6), and Equation (7), the generator output at the time of the reference curve is used.

最後に、本発明の定式化を説明する。本発明では、式(12)の目的関数を採用している。元の式(8)に対して、起動停止変数を収束させるために、一例としてβおよびq,rを導入している。起動停止変数を収束させるためには、起動停止変数が0に近いものは0へ、1に近いものは1へ近づくように、目的関数あるいは制約条件を修正すればよい。ここで、燃料費関数のなかで、定数項の部分のみを修正する理由は、起動停止変数に直接影響を与えていることおよび発電機出力の線形、二次の項は最適な負荷配分を維持するために等λ配分となる条件を維持させるためである。また、後者のq,rはそれぞれ、1に近づけるあるいは0に近づけるためのものもである。ここで、使用されている定数は図5の画面で設定する。表形式の項目は直接数値を設定する。また、「目的関数の係数修正」あるいは「制約条件の修正」はどちらか一方を選択するものとする。この選択は処理S101のときに実施されているものとする。前者の「目的関数の係数修正」を選択すると、繰り返し計算の途中で得られている起動停止変数uitの値により、式(12)の目的関数の係数を修正するものとする。また、後者の「制約条件の修正」を選択すると、後で出てくる式(14)の制約条件を修正するとともに、式(12)の目的関数の係数も修正する。 Finally, the formulation of the present invention will be described. In the present invention, the objective function of Expression (12) is adopted. For the original equation (8), β and q, r are introduced as an example in order to converge the start / stop variables. In order to converge the start / stop variable, the objective function or the constraint condition may be modified so that the start / stop variable close to 0 is close to 0 and that close to 1 is close to 1. Here, the reason for correcting only the constant term part in the fuel cost function is that it has a direct influence on the start / stop variable and that the linearity of the generator output and the second order term maintain the optimal load distribution. This is to maintain the condition of equal λ distribution. Further, the latter q and r are each for approaching 1 or approaching 0. Here, the constants used are set on the screen of FIG. The tabular items are set directly. Also, either “Correction of objective function coefficient” or “Correction of constraint” is selected. It is assumed that this selection is performed at the time of processing S101. When the former “Correction of coefficient of objective function” is selected, the coefficient of the objective function of Expression (12) is corrected by the value of the start / stop variable u it obtained during the repeated calculation. When the latter “correction of constraint” is selected, the constraint of equation (14) that appears later is corrected and the coefficient of the objective function of equation (12) is also corrected.

Figure 0004577232
Figure 0004577232

図5に出てくる係数βは式(13)で使用する。図5のa,b,c,d,e,fの値は、起動停止変数を収束させるために設定するものであり、起動停止変数が0に近いものは0へ、1に近いものは1へ近づくように、設定するものとする。   The coefficient β shown in FIG. 5 is used in equation (13). The values of a, b, c, d, e, and f in FIG. 5 are set in order to converge the start / stop variables. When the start / stop variable is close to 0, the value is 0. It shall be set so that it approaches.

Figure 0004577232
Figure 0004577232

式(14)は起動停止変数は0または1であるが、二次計画法で解ける定式化とするために、0から1の実数変数に緩和している。   In equation (14), the start / stop variable is 0 or 1, but it is relaxed to a real variable from 0 to 1 in order to formulate it by a quadratic programming method.

Figure 0004577232
Figure 0004577232

最後に、式(15)は図4の需要増加時間帯,需要減少時間帯を利用して実数変数に緩和したuitの時間的な連続性を維持するために設定している。つまり、1時間おきに起動停止を繰り返したりしないようにするためである。例えば、図4のように設定すると、1日の間で起動,停止はそれぞれ高々1回となり、実際の運用と同じくすることができる。但し、稀ではあるが小規模のガスタービン等では、1日に2回起動,停止することもある。 Finally, equation (15) is set in order to maintain increased demand time slot of FIG. 4, the temporal continuity of u it was relaxed to real variables by using the lower demand times. That is, it is to prevent repeated start and stop every hour. For example, when setting is made as shown in FIG. 4, the start and stop are each performed once a day at most, and can be the same as the actual operation. However, although it is rare, a small-scale gas turbine or the like may start and stop twice a day.

Figure 0004577232
Figure 0004577232

例えば、図4の需要のカーブ0に対して、式(15)の制約条件を設定すると、需要増加時間帯では起動停止変数は減少せず、需要減少時間帯では起動停止変数は増加しない。また、発電機には一旦停止すると、最小連続停止時間は連続して停止する必要がある。同様に、一旦起動すると、最小連続運転時間は連続して運転する必要がある。通常は、最大需要の時刻を中心として、需要が大きい時刻から最小連続運転時間の間は起動停止変数の変化をゼロにし、最小需要の時刻を中心として、需要が小さい時刻から最小連続停止時間の間は起動停止変数の変化をゼロにすることになる。この結果、最終的な起動停止変数は0または1になるが、この起動停止変数を実数変数に緩和していることから、発電機によっては、起動停止変数は図5−1に示すようになっている。   For example, if the constraint condition of Expression (15) is set for the demand curve 0 in FIG. 4, the start / stop variable does not decrease during the demand increase time period, and the start / stop variable does not increase during the demand decrease time period. Further, once the generator is stopped, the minimum continuous stop time needs to be continuously stopped. Similarly, once activated, the minimum continuous operation time must be continuously operated. Normally, the change in the start / stop variable is set to zero between the time when the demand is high and the minimum continuous operation time, centering on the time of the maximum demand, and the time of the minimum continuous stop time from the time when the demand is low, centering on the time of the minimum demand. During this time, the change in the start / stop variable is set to zero. As a result, the final start / stop variable becomes 0 or 1, but since this start / stop variable is relaxed to a real variable, the start / stop variable is as shown in FIG. 5-1, depending on the generator. ing.

これら図4,図5の設定は入力装置6により運用者等が個々に設定することができる。   4 and 5 can be individually set by an operator or the like by the input device 6.

以上の制約条件を追加して、定式化する。第1回目の繰り返し計算が処理S103としている。   The above constraints are added and formulated. The first iterative calculation is the processing S103.

処理S104では、計画計算部25により二次計画法,線形計画法を用いて、処理S103あるいは処理S107,処理S108で設定した目的関数,制約条件のもので、火力発電機の起動停止変数と火力,揚水発電機出力の最適解を同時に求める。ただし、火力発電機の起動停止変数は0または1に収束しているとは限らない。このため、処理S107,処理S108を繰り返していくことにより、全ての起動停止変数を停止0または運転1に確定させる。また、制約条件を満たす計画が作成できないときには、エラーメッセージを表示するものとする。   In the process S104, the plan calculation unit 25 uses the quadratic programming method and the linear programming method, and the objective function and the constraint conditions set in the process S103, the process S107, and the process S108. , The optimum solution of the pumped generator output is obtained simultaneously. However, the start / stop variable of the thermal power generator does not necessarily converge to 0 or 1. For this reason, by repeating the processing S107 and the processing S108, all the start / stop variables are fixed to the stop 0 or the operation 1. An error message is displayed when a plan that satisfies the constraint conditions cannot be created.

処理S105では計算結果処理部28により、処理S104での計算結果を表示装置7に表示する。また、計算結果をデータベース9に格納する。格納するものには、火力発電機の起動停止変数,火力および揚水発電機の出力,発電コスト,燃料費,起動費および使用した制約条件,目的関数等である。画面表示例を図6,図7に示す。図6は起動停止変数の値を示す。横軸方向が時間軸である。縦軸方向が火力発電機のユニット番号である。表の中の数値が対象時刻の火力発電機の起動停止変数を示している。停止0または運転1の他に0.6 などの起動停止が確定していないものがある。また、同一ユニットに対して、起動停止確定している確定率を表示している。例えば、1時の確定率は80%となっているが、ユニット10台のうち、8台の起動停止変数が0または1になっていることを示している。同様に各時間での起動停止の確定率も表示している。表の上部には全体の起動停止の確定率78.4% を表示している。起動停止が確定してくると、確定率が100%に近くなる。通常は同一時間においては、1つのユニットの起動停止が確定できないことが多い。また、図7は起動停止変数が図6のときの図火力発電機の出力を示している。図6と同様に横軸が時間で、縦軸が火力ユニットの番号である。この出力も本来の火力ユニットの出力下限制約より小さい値になっていることもある。   In process S105, the calculation result processing unit 28 displays the calculation result in process S104 on the display device 7. The calculation result is stored in the database 9. The stored items include thermal generator start / stop variables, thermal power and pumped-storage generator output, power generation cost, fuel cost, startup cost, constraints used, objective function, and the like. Examples of screen display are shown in FIGS. FIG. 6 shows the values of the start / stop variables. The horizontal axis direction is the time axis. The vertical axis direction is the unit number of the thermal power generator. The numerical values in the table indicate the start / stop variables of the thermal power generator at the target time. In addition to stop 0 or operation 1, there are some that have not been confirmed to start or stop, such as 0.6. In addition, the fixed rate at which start / stop is fixed for the same unit is displayed. For example, although the determination rate at 1 o'clock is 80%, the start / stop variable of 8 units among the 10 units is 0 or 1. Similarly, the fixed rate of starting and stopping at each time is also displayed. At the top of the table, the overall start / stop decision rate of 78.4% is displayed. When start / stop is confirmed, the confirmation rate becomes close to 100%. Normally, it is often impossible to determine the start / stop of one unit at the same time. FIG. 7 shows the output of the thermal power generator when the start / stop variable is that shown in FIG. Similar to FIG. 6, the horizontal axis represents time, and the vertical axis represents the number of the thermal power unit. This output may also be smaller than the output lower limit constraint of the original thermal power unit.

処理S106では計画計算部25にてすべての火力発電機の起動停止変数が停止0または起動1に収束しているかどうかをチェックする。すべての起動停止変数が0または1に収束していれば、すべての制約条件を満たす起動停止計画を作成できたことになり、処理S109に進む。そうでなければ、処理S113に進む。まず処理S113に進んだときについて説明する。   In process S106, the plan calculation unit 25 checks whether the start / stop variables of all the thermal power generators have converged to stop 0 or start 1. If all the start / stop variables have converged to 0 or 1, it means that a start / stop plan that satisfies all the constraint conditions has been created, and the process proceeds to step S109. Otherwise, the process proceeds to process S113. First, the case where the processing proceeds to step S113 will be described.

処理S113では、数理計画法での繰り返し計算をカウントアップしている。繰り返し回数が図4で設定した回数を超過しているときには、処理S107,処理S108では図5で設定した条件を利用して、強制的に起動停止を確定させる。すなわち、図5の(4)の条件を満たすuitはすべて停止(0)、それ以外のuitは運転(1)にする。これにより、次に判定処理S106に来たときには、終了方向の処理S109へ進むことになる。繰り返し回数に達していないときには、図5で設定した条件と式(12)の目的関数、式(14),式(15)および式(2)から式(7)の制約条件を用いて、制約条件,目的関数を修正する。すなわち、処理S107では、起動停止変数上下限修正部26により図5の設定にしたがって、起動停止変数の上下限制約を修正する。次に、処理S108では目的関数修正部27により、図5の設定にしたがって、目的関数を修正する。 In the process S113, the repeated calculation by the mathematical programming method is counted up. When the number of repetitions exceeds the number set in FIG. 4, in step S107 and step S108, the start / stop is forcibly confirmed using the conditions set in FIG. That is, all u it satisfying the condition of (4) in FIG. 5 is stopped (0), and other u it is set to operation (1). Thereby, when it comes to determination process S106 next, it will progress to process S109 of an end direction. When the number of repetitions has not been reached, the constraints set using the conditions set in FIG. 5 and the objective function of equation (12), the constraints of equations (14), (15), and (2) to (7) are used. Correct the condition and objective function. That is, in the process S107, the start / stop variable upper / lower limit correction unit 26 corrects the upper / lower limit constraint of the start / stop variable according to the setting of FIG. In step S108, the objective function correcting unit 27 corrects the objective function according to the setting shown in FIG.

以上の処理S107,処理S108により制約条件,目的関数を修正することになる。この修正した結果をもとに、処理S104により再度計画を作成する。このループに沿って、繰り返し計画を作成することにより、起動停止変数を停止0または起動1に収束させていく。   The constraint conditions and the objective function are corrected by the above processing S107 and processing S108. Based on this corrected result, a plan is created again in step S104. By creating a repetitive plan along this loop, the start / stop variable is converged to stop 0 or start 1.

以上の繰り返し計算により、処理S106ですべての起動停止変数0または1に確定した場合は処理S109に進む。   If it is determined in the process S106 that all the start / stop variables are 0 or 1, the process proceeds to the process S109.

以上により起動停止変数は0または1に確定できたことになる。そこで、起動停止変数を0または1に収束させるために目的関数の式(12)でβ,q,rを導入したが、本来の目的を達成できたので元の目的関数にするために、処理S109では目的関数式(12)でβ=1,q=0,r=0に設定する。また制約条件の式(14)は上下限が0または1に変更となる。あるいは、変数としての扱いをやめて、定数として扱うこともできる。   Thus, the start / stop variable can be determined to be 0 or 1. Therefore, β, q, and r were introduced in the objective function equation (12) in order to converge the start / stop variable to 0 or 1, but since the original purpose was achieved, the processing was performed to obtain the original objective function. In S109, β = 1, q = 0, and r = 0 are set in the objective function equation (12). Further, the upper and lower limits of the constraint condition expression (14) are changed to 0 or 1. Alternatively, you can stop treating it as a variable and treat it as a constant.

処理S109で設定した目的関数,制約条件をもとに、計画計算部25を用いて、処理S110により二次計画法等の数理計画法により制約条件を満たし、最適化は発電機出力を決定する。   Based on the objective function and constraints set in step S109, the plan calculation unit 25 is used to satisfy the constraints by mathematical programming such as quadratic programming in step S110, and optimization determines the generator output. .

最後に計算結果処理部28を用いて処理S111の計算結果を表示装置7に表示する、あるいは電力システム11に出力する。図8から図14に計算結果の画面表示の例を示す。図8,図9は図6,図7と同じ形式のもので、起動停止変数が収束したときの値を表示している。図8では、0,1確定率は100%になっており、すべての火力発電機の起動停止が0または1に確定していることを確認できる。図9では、発電機を運転しているときには、最低出力から最大出力の範囲内になっている。図10は繰り返し計算の経過を確認するための画面である。時刻のボタンをクリックして増加あるいは減少あるいは直接見たい時刻を入力して、「表示」をクリックすると、その下に選択した時刻の全火力ユニットの起動停止変数の繰り返し計算によりどう変化したかを棒グラフで表示する。黒くなっているユニットは最終的に運転(1)になったもので、棒グラフの高さが、起動停止変数の値を示している。白い棒グラフがあるものは、最終的に停止(0)になったことを示している。図11は各繰り返し計算での各需要カーブの発電コストの推移を示している。太い線が基準需要カーブのときの発電コストを示している。図12は各繰り返し計算での全火力ユニット、全時刻を対象としたときの起動停止変数が0または1に確定した確定率の推移を示している。この例では5回目で0,1に確定していることがわかる。図13は最終結果の時間毎の発電コストを示している。図10と同様に、時刻を設定して、「表示」をクリックすることで、設定した時刻の発電コストを表示させることができる。発電コストは需要カーブ毎の発電コストの棒グラフと表を表示している。また、需要カーブの確率を考慮したときの重み付きコスト分布の棒グラフを示す。表は各需要カーブの重み付き発電コストを示しており、合計がその時刻のコスト期待値になる。図14は図13の時間毎の発電コストを計画対象期間の全時刻分を加算したものである。右下の表の発電コストが全体のコスト期待値で、その内訳として燃料費および起動費を表示している。   Finally, the calculation result processing unit 28 is used to display the calculation result of the process S111 on the display device 7 or output it to the power system 11. FIG. 8 to FIG. 14 show examples of calculation result screen display. FIGS. 8 and 9 are of the same format as FIGS. 6 and 7, and display values when the start / stop variables converge. In FIG. 8, the 0, 1 determination rate is 100%, and it can be confirmed that the start / stop of all the thermal power generators is fixed to 0 or 1. In FIG. 9, when the generator is operating, the range is from the minimum output to the maximum output. FIG. 10 is a screen for confirming the progress of repeated calculations. Click the time button and enter the time you want to increase or decrease, or click on “Display”, and click “Display” to see how it has changed due to repeated calculation of the start / stop variables of all thermal power units at the time selected below. Display as a bar graph. The unit that is black is finally in operation (1), and the height of the bar graph indicates the value of the start / stop variable. The one with the white bar graph indicates that it has finally stopped (0). FIG. 11 shows the transition of power generation cost for each demand curve in each repetitive calculation. The thick line shows the power generation cost when the standard demand curve is shown. FIG. 12 shows the transition of the fixed rate in which the start / stop variable is fixed to 0 or 1 when all the thermal power units and all times are subjected to each repetition calculation. In this example, it can be seen that the values are fixed to 0 and 1 at the fifth time. FIG. 13 shows the power generation cost for each hour of the final result. Similarly to FIG. 10, the power generation cost at the set time can be displayed by setting the time and clicking “Display”. The power generation cost displays a bar graph and a table of power generation cost for each demand curve. Moreover, the bar graph of weighted cost distribution when the probability of a demand curve is considered is shown. The table shows the weighted power generation costs for each demand curve, and the total is the expected cost at that time. FIG. 14 is obtained by adding the power generation costs for each hour in FIG. 13 to all the times in the planning target period. The power generation cost in the table on the lower right is the expected cost for the entire project. The breakdown is fuel cost and startup cost.

最後に、処理S110ではこれらの計算結果をデータベース9に格納する。格納するものは、起動停止変数,発電機出力,発電コスト,燃料費,起動費等である。   Finally, in process S110, these calculation results are stored in the database 9. The stored items are start / stop variables, generator output, power generation cost, fuel cost, start-up cost, and the like.

上記本発明の実施例によれば、需要予測誤差がある場合においても、各需要に対する発電コストにその需要の発生確率を重み係数として、両者の積の総和を最小化する目的関数に設定し、また、0,1変数である火力発電機の起動停止変数を実数変数に緩和して、基準となる電力需要予測の需要変化を考慮した起動停止変数の時間変化制約を導入し、LNG消費制約および揚水ダム水位制約を考慮した発電費用最小となる全時刻の発電所計画を作成することができる。   According to the embodiment of the present invention, even when there is a demand prediction error, the generation probability of the demand is set as a weighting factor in the power generation cost for each demand, and the objective function for minimizing the sum of the products of both is set. In addition, the start and stop variables of thermal power generators, which are 0 and 1 variables, are relaxed to real variables, and time change constraints of start and stop variables that take into account changes in the demand for power demand prediction as a reference are introduced. It is possible to create a full-time power plant plan that minimizes the power generation cost considering the pumping dam water level constraint.

第二の実施例は、発電機の燃料費関数が気温等により効率が変化することを考慮するための方法である。例えば、火力発電機がガスタービンの場合、気温が高くなると効率が低下し、気温が低いと効率が向上する特性がある。また、夏期に気温が高くなると冷房需要の増加により電力需要は増加し、逆に冬期は気温が高くなると暖房需要が減少して、電力需要が減少する。このように、気温により、火力発電機の燃料消費特性および電力需要が変化する。このため、予想気温の確率分布をもとに、各予想気温に対する火力発電機の燃料消費特性および予想需要が変化することを考慮して、発電費用を最小化することが必要になる。第一の実施例では、需要が複数あって、発電機特性は固定であったが、第二の実施例は、需要および発電機特性がともに気温により確率的に変化する例である。この場合も、第一の実施例と同様の方法で、処理することにより、最適な起動停止計画および負荷配分計画を作成することができる。目的関数,制約条件を設定すればよい。但し、目的関数の燃料費関数の係数が異なることおよび各燃料費関数に対して重み係数を設定して、同様の処理で最適な起動停止計画および負荷配分計画を作成することができる。第一の実施例と異なる点について以下に述べる。処理S101では、新たに次の計画作成条件の設定を追加する。予想気温の確率分布情報を追加して入力する。次に処理S102では、気温と発電機の燃料費特性の関係を読み込む。更に、予想気温に対する燃料費特性を決定する。更に処理S103では、予想気温の確率分布を用いて、予想気温の燃料費特性を用いる確率を予想気温の確率分布を用いて、目的関数を設定することになる。第一の実施例の目的関数は式(8)になるのに対して、第二の実施例での目的関数は式(16)のようになる。第二の実施例の目的関数では予想気温の違いを変数,関数の上付きのkで表しており、この上付きのkが燃料費関数fi k(P)に含まれていることが第一の実施例と異なる点である。 The second embodiment is a method for taking into consideration that the efficiency of the fuel cost function of the generator varies with temperature or the like. For example, when the thermal power generator is a gas turbine, the efficiency decreases when the temperature rises, and the efficiency improves when the temperature is low. In addition, when the temperature rises in summer, the demand for electric power increases due to an increase in cooling demand. Conversely, in winter, when the temperature rises, the demand for heating decreases and the demand for electric power decreases. Thus, the fuel consumption characteristics and power demand of the thermal power generator change depending on the temperature. For this reason, it is necessary to minimize the power generation cost in consideration of changes in the fuel consumption characteristics and the predicted demand of the thermal power generator for each predicted temperature based on the probability distribution of the predicted temperature. In the first embodiment, there are a plurality of demands and the generator characteristics are fixed, but the second embodiment is an example in which both the demand and the generator characteristics change stochastically with the temperature. Also in this case, an optimum start / stop plan and load distribution plan can be created by processing in the same manner as in the first embodiment. An objective function and constraint conditions can be set. However, it is possible to create optimum start / stop plans and load distribution plans by the same process by setting different coefficients of the fuel cost function of the objective function and setting a weight coefficient for each fuel cost function. Differences from the first embodiment will be described below. In process S101, the setting of the next plan creation condition is newly added. Enter the probability distribution information of the expected temperature. Next, in process S102, the relationship between the temperature and the fuel cost characteristic of the generator is read. Furthermore, the fuel cost characteristic with respect to the expected temperature is determined. Further, in step S103, an objective function is set using the probability distribution of the predicted temperature and the probability distribution of the fuel temperature characteristic of the predicted temperature using the probability distribution of the predicted temperature. The objective function in the first embodiment is expressed by equation (8), whereas the objective function in the second embodiment is expressed by equation (16). In the objective function of the second embodiment, the difference in the predicted temperature is represented by a superscript k of a variable or function, and the superscript k is included in the fuel cost function f i k (P). This is different from the first embodiment.

Figure 0004577232
Figure 0004577232

よって、燃料費関数が含まれない数式である、式(9),式(10),式(11),式(13),式(14)および式(15)は第二の実施例でも第一の実施例と同じでよいことになる。また、収束させるために使用する目的関数式(12)は燃料費関数が含まれるので、第二の実施例では、式(17)のようになる。   Therefore, the formulas (9), (10), (11), (13), (14), and (15) that do not include the fuel cost function are the same in the second embodiment. It may be the same as one embodiment. Further, since the objective function formula (12) used for convergence includes the fuel cost function, the formula (17) is obtained in the second embodiment.

Figure 0004577232
Figure 0004577232

上記のように目的関数を修正して、第一の実施例と同様に繰り返し計算により、設定した制約条件を満たし、各発電機の起動停止変数uitを0または1に収束させることができる。 Modify the objective function as described above, by repeatedly calculated as in the first embodiment, satisfy the set constraints, a start-stop variable u it for each generator can be converged to 0 or 1.

以上の方法により、予想気温の確率分布等を考慮して、発電機の燃料費関数の変化確率を考慮して発電費用を最小化できる発電所運用計画を作成することができる。   By the above method, it is possible to create a power plant operation plan that can minimize the power generation cost in consideration of the probability distribution of the predicted temperature and the change probability of the fuel cost function of the generator.

発電所運用計画作成の処理フロー。Process flow for creating a power plant operation plan. 発電所運用計画作成装置の全体構成。Overall configuration of power plant operation plan creation device. 発電所運用計画作成装置の詳細構成。Detailed configuration of the power plant operation plan creation device. 電力需要カーブ等の設定画面の例。An example of a setting screen such as a power demand curve. 目的関数,制約条件を修正する画面の例。Example of a screen for modifying objective functions and constraints. 起動停止変数の時間変化のグラフの例。The example of a graph of the time change of a start stop variable. 起動停止変数を表示する画面の例。An example of a screen displaying start / stop variables. 発電機出力を表示する画面の例。An example of a screen displaying the generator output. 起動停止変数を表示する画面の例。An example of a screen displaying start / stop variables. 発電機出力を表示する画面の例。An example of a screen displaying the generator output. 起動停止変数の推移を表示する画面の例。An example of a screen that displays the transition of start / stop variables. 各需要カーブの発電コストの推移を表示する画面の例。The example of the screen which displays transition of the power generation cost of each demand curve. 起動停止変数の確定率の推移を表示する画面の例。The example of the screen which displays the transition of the fixed rate of the start / stop variable. 作成した計画の時間毎の発電コストを表示する画面の例。The example of the screen which displays the power generation cost for every hour of the created plan. 作成した計画の全時間の発電コストを表示する画面の例。An example of a screen that displays the total cost of power generation for the created plan.

符号の説明Explanation of symbols

1…発電所需給運用計画作成装置、2…中央演算処理装置CPU、3…主記憶装置、4…入出力装置、5…外部記憶装置、6…入力装置、7…表示装置、8…印字装置、9…データベース、10…読取装置、11…電力システム、12…電力系統、20…運用計画作成処理部、21…コントロール部、22…計画作成条件設定部、23…データ読込部、
24…定式化部、25…計画計算部、26…起動停止変数上下限修正部、27…目的関数修正部、28…計算結果処理部。
DESCRIPTION OF SYMBOLS 1 ... Power station supply and demand operation plan preparation device, 2 ... Central processing unit CPU, 3 ... Main storage device, 4 ... Input / output device, 5 ... External storage device, 6 ... Input device, 7 ... Display device, 8 ... Printing Device: 9 ... Database 10 ... Reading device 11 ... Power system 12 ... Power system 20 ... Operation plan creation processing unit 21 ... Control unit 22 ... Plan creation condition setting unit 23 ... Data reading unit
24 ... Formulation unit, 25 ... Plan calculation unit, 26 ... Start / stop variable upper / lower limit correction unit, 27 ... Objective function correction unit, 28 ... Calculation result processing unit.

Claims (10)

2つ以上の同一時刻における需要予測と、前記需要予測発生確率に相当する重み係数を受付け記憶装置に記憶させる手段と、
前記記憶装置に記憶された前記2つ以上の同一時刻における需要予測と、前記需要予測発生確率に相当する重み係数により目的関数を設定する手段を有する発電所運計画の作成装置において、
各需要に対して、0から1の実数変数として緩和された所定発電機を含んで供給バランスの制約を満たすような発電機起動の組合せ問題に係る演算を介して定まる起動停止計画および負荷配分したときの発電コストの重み係数を受付け記憶装置に記憶する手段と、
前記需要毎の発電コストに需要に対する発電コストの重み係数をかけたものを前記目的関数に設定し、前記目的関数に基づいて発電所運用計画を行う手段とを有することを特徴とする発電所運用計画の作成装置。
Means for storing demand prediction at two or more same times and a weighting factor corresponding to the demand forecast occurrence probability in an accepting storage device;
In the power plant operation plan creation device having means for setting an objective function with a demand coefficient corresponding to the demand prediction occurrence probability and the demand prediction at the two or more same times stored in the storage device,
For each demand, the start / stop plan and load distribution determined through the calculation related to the generator start-up combination problem including the predetermined generator relaxed as a real variable from 0 to 1 and satisfying the supply balance constraint Means for storing the weighting factor of the power generation cost at the time in the receiving storage device;
A power plant operation comprising: means for multiplying the power generation cost for each demand by a weighting factor of the power generation cost with respect to the demand as the objective function, and performing a power plant operation plan based on the objective function Plan creation device.
請求項1において、
前記目的関数を数理計画法を用いて起動停止計画および負荷配分計画を最適化する手段を有することを特徴とする発電所運用計画の作成装置。
In claim 1,
An apparatus for creating a power plant operation plan, comprising means for optimizing a start / stop plan and a load distribution plan using a mathematical programming method for the objective function.
請求項1において、
需要増加時間帯または需要減少時間帯を受付け記憶装置に記憶する手段とを有することを特徴とする発電所運用計画の作成装置。
In claim 1,
An apparatus for creating a power plant operation plan, comprising means for receiving a demand increase time zone or a demand decrease time zone in a storage device.
請求項1において
需要増加時間帯と需要減少時間帯の起動停止変数の時間制約を受付け記憶装置に記憶する手段を有することを特徴とする発電所運用計画の作成装置。
The power plant operation plan creation device according to claim 1, further comprising means for receiving and storing in a storage device time constraints of start / stop variables in a demand increase time zone and a demand decrease time zone.
請求項1において、
起動停止変数の値あるいは各発電機の発電機出力をもとに、目的関数あるいは起動停止変数の上下限制約値を修正する手段を有することを特徴とする発電所運用計画の作成装置。
In claim 1,
An apparatus for creating a power plant operation plan, comprising means for correcting an upper limit / lower limit constraint value of an objective function or a start / stop variable based on a value of a start / stop variable or a generator output of each generator.
請求項1において、
予備力,同一中央同時起動禁止制約または燃料消費制約の少なくとも1つの制約を受付け記憶装置に記憶すことを特徴とする発電所運用計画の作成装置。
In claim 1,
Reserve, the same central simultaneous start prohibition constraints or power plant operation plan creation device, characterized in that you stored in a storage device accepting at least one constraint of fuel consumption restrictions.
請求項1において、
目的関数の係数,目的関数修正関数,起動停止変数の値,発電機出力の値,起動停止確定率,目的関数の係数の繰返回数に対する推移,起動停止変数の繰返回数に対する推移,発電機出力の繰返回数に対する推移,起動停止確定率の繰返回数に対する推移,燃料消費量の繰返回数に対する推移,ダム水位の繰返回数に対する推移のうち少なくとも1つを出力し表示する手段を有することを特徴とする発電所運用計画の作成装置。
In claim 1,
Objective function coefficient, objective function correction function, start / stop variable value, generator output value, start / stop decision rate, transition of objective function coefficient with respect to the number of repetitions, transition with respect to the number of repetitions of the start / stop variable, generator A means for outputting and displaying at least one of a transition with respect to the number of repetitions of output, a transition with respect to the number of repetitions of the start / stop fixed rate, a transition with respect to the number of repetitions of fuel consumption, and a transition with respect to the number of repetitions of the dam water level; This is a power plant operation plan creation device.
請求項1において、
予想気温の確率分布を受付け記憶装置に記憶する手段と、
気温と発電機の燃料特性の関係を受付け記憶装置に記憶する手段と、
前記予想気温に対する燃料費特性を作成し記憶装置に記憶する手段と、
前記予想気温に対する燃料費特性を反映させて目的関数を作成する手段を有することを特徴とする発電所運用計画の作成装置。
In claim 1,
Means for accepting and storing the probability distribution of the expected temperature in the receiving storage device;
Means for receiving and storing the relationship between the temperature and the fuel characteristics of the generator in a storage device;
Means for creating a fuel cost characteristic for the expected temperature and storing it in a storage device;
An apparatus for creating a power plant operation plan, comprising means for creating an objective function reflecting the fuel cost characteristic with respect to the predicted temperature.
2つ以上の同一時刻における需要予測と、前記需要予測発生確率に相当する重み係数を受付け記憶装置に記憶させるステップと、
前記記憶装置に記憶された前記2つ以上の同一時刻における需要予測と、前記需要予測発生確率に相当する重み係数により目的関数を設定するステップを有する発電所運計画の作成方法において、
各需要に対して、0から1の実数変数として緩和された所定発電機を含んで供給バランスの制約を満たすような発電機起動の組合せ問題に係る演算を介して定まる起動停止計画および負荷配分したときの発電コストの重み係数を受付け記憶するステップと、
前記需要毎の発電コストに需要に対する発電コストの重み係数をかけたものを前記目的関数に設定し、前記目的関数に基づいて発電所運用計画を行うステップとを有することを特徴とする発電所運用計画の作成方法。
Two or more demand forecasts at the same time, and a weighting coefficient corresponding to the demand forecast occurrence probability is stored in an accepting storage device;
In the method for creating a power plant operation plan, comprising the step of setting an objective function by means of a weighting coefficient corresponding to the demand prediction occurrence probability and the demand prediction at the two or more same times stored in the storage device,
For each demand, the start / stop plan and load distribution determined through the calculation related to the generator start-up combination problem including the predetermined generator relaxed as a real variable from 0 to 1 and satisfying the supply balance constraint Accepting and storing the weighting factor of the power generation cost when
Power plant operation, characterized in that a step of performing a power plant operation plan on the basis of a multiplied by the weighting factor of the power generation cost for demand generation cost for each of the demand to the set objective function, the objective function How to create a plan.
請求項9において、
前記目的関数を数理計画法を用いて起動停止計画および負荷配分計画を最適化するステップを有することを特徴とする発電所運用計画の作成方法。
In claim 9,
A method for creating a power plant operation plan, comprising: a step of optimizing a start / stop plan and a load distribution plan using a mathematical programming method for the objective function.
JP2006044682A 2006-02-22 2006-02-22 Power plant operation plan creation device and method Expired - Fee Related JP4577232B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2006044682A JP4577232B2 (en) 2006-02-22 2006-02-22 Power plant operation plan creation device and method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2006044682A JP4577232B2 (en) 2006-02-22 2006-02-22 Power plant operation plan creation device and method

Publications (2)

Publication Number Publication Date
JP2007228676A JP2007228676A (en) 2007-09-06
JP4577232B2 true JP4577232B2 (en) 2010-11-10

Family

ID=38549935

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2006044682A Expired - Fee Related JP4577232B2 (en) 2006-02-22 2006-02-22 Power plant operation plan creation device and method

Country Status (1)

Country Link
JP (1) JP4577232B2 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20220057768A1 (en) * 2018-12-06 2022-02-24 Hitachi, Ltd. Power generation planning apparatus
US12001181B2 (en) * 2018-12-06 2024-06-04 Hitachi, Ltd. Power generation planning apparatus

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4980303B2 (en) * 2008-06-26 2012-07-18 三菱電機株式会社 Power generation cost evaluation method, power generation cost evaluation apparatus, power generation cost evaluation program, power generation plan creation system, and power generation facility start / stop control system
JP4896084B2 (en) * 2008-06-30 2012-03-14 株式会社日立製作所 Generator operation plan creation method
JP5248372B2 (en) * 2009-03-11 2013-07-31 株式会社東芝 Power generation plan creation method, device, program, and storage device
JP5412332B2 (en) * 2010-03-05 2014-02-12 株式会社日立製作所 Generator operation plan determination device, generator operation plan determination method, and generator operation plan determination program
JP5452714B2 (en) * 2010-05-10 2014-03-26 三菱電機株式会社 Power generation plan creation device
JP5494279B2 (en) * 2010-06-23 2014-05-14 東京電力株式会社 Thermal power generator operation planning system
JP5494282B2 (en) * 2010-06-24 2014-05-14 東京電力株式会社 Thermal power generator operation planning system
US8565929B2 (en) * 2010-11-17 2013-10-22 General Electric Company Enhanced plant level support of grid reconstruction
CN102611106B (en) * 2012-03-31 2014-04-16 凯里供电局 Maximum load supply capability evaluation method of medium-voltage power distribution network for loop power supply
JP6401553B2 (en) * 2014-09-17 2018-10-10 一般財団法人電力中央研究所 Operation plan formulation device, operation plan formulation method, and operation plan formulation program
JP2016126453A (en) * 2014-12-26 2016-07-11 川崎重工業株式会社 Device and method for making optimum energy system operation plan, program, and recording medium
JP7018287B2 (en) * 2017-10-06 2022-02-10 株式会社日立製作所 Market-ready reliability management device and market-ready reliability management system
CN107832979B (en) * 2017-12-06 2021-05-11 浙江大学 Factory comprehensive energy system economic optimization scheduling method considering energy cascade utilization
JP7023766B2 (en) * 2018-03-26 2022-02-22 東芝エネルギーシステムズ株式会社 Hydropower plant output distribution device and hydropower system
JP7008580B2 (en) * 2018-06-06 2022-01-25 三菱電機株式会社 Power generation planning equipment
CN108984875B (en) * 2018-07-02 2022-11-08 中国电力工程顾问集团西北电力设计院有限公司 Method for calculating capacity of photo-thermal unit participating in power balance
JP7339747B2 (en) * 2019-03-13 2023-09-06 株式会社日立製作所 Control device for energy storage equipment
CN111415041B (en) * 2020-03-20 2024-01-05 海南电网有限责任公司 Method for evaluating economy of power grid planning scheme
CN112686467A (en) * 2021-01-12 2021-04-20 南方电网科学研究院有限责任公司 Large-user production operation optimization method and device containing self-contained power plant and terminal equipment

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001258157A (en) * 2000-03-09 2001-09-21 Hitachi Ltd Device and method of developing operation plan for generator and recording medium storing operation plan program for generator
JP2003016374A (en) * 2001-06-28 2003-01-17 Toshiba Corp Evaluating method and drawing-up method for power generating facility plan, and program
JP2006050730A (en) * 2004-08-03 2006-02-16 Hitachi Ltd Operation plan creating method and apparatus of supply capability

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001258157A (en) * 2000-03-09 2001-09-21 Hitachi Ltd Device and method of developing operation plan for generator and recording medium storing operation plan program for generator
JP2003016374A (en) * 2001-06-28 2003-01-17 Toshiba Corp Evaluating method and drawing-up method for power generating facility plan, and program
JP2006050730A (en) * 2004-08-03 2006-02-16 Hitachi Ltd Operation plan creating method and apparatus of supply capability

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20220057768A1 (en) * 2018-12-06 2022-02-24 Hitachi, Ltd. Power generation planning apparatus
US12001181B2 (en) * 2018-12-06 2024-06-04 Hitachi, Ltd. Power generation planning apparatus

Also Published As

Publication number Publication date
JP2007228676A (en) 2007-09-06

Similar Documents

Publication Publication Date Title
JP4577232B2 (en) Power plant operation plan creation device and method
Moretti et al. An efficient robust optimization model for the unit commitment and dispatch of multi-energy systems and microgrids
Cobos et al. Robust energy and reserve scheduling considering bulk energy storage units and wind uncertainty
Nikoobakht et al. Flexible power system operation accommodating uncertain wind power generation using transmission topology control: an improved linearised AC SCUC model
Hu et al. Robust SCUC considering continuous/discrete uncertainties and quick-start units: A two-stage robust optimization with mixed-integer recourse
US8554382B2 (en) Multivariable control of regulation and fast demand response in electrical grids
US20210381120A1 (en) Hydrogen-system control device and hydrogen-system control method
US10452041B2 (en) Gas turbine dispatch optimizer real-time command and operations
Zhang et al. Two-stage robust security-constrained unit commitment model considering time autocorrelation of wind/load prediction error and outage contingency probability of units
Khoshjahan et al. Advanced bidding strategy for participation of energy storage systems in joint energy and flexible ramping product market
US20180284706A1 (en) Gas turbine dispatch optimizer
JP6197689B2 (en) Operation plan support program, operation plan support method, and operation plan support apparatus
JP4179293B2 (en) Generator operation plan creation device and generator operation plan creation program
JP2006352933A (en) Power supply control method and system
US10410301B2 (en) Operation plan planning system for planning an operation plan of a power plant and operation plan planning method for planning an operation plan of a power plant
JP2022045625A (en) Power generation plan creation system, power generation plan creation method, and program
JP5412332B2 (en) Generator operation plan determination device, generator operation plan determination method, and generator operation plan determination program
JP4075872B2 (en) Supply plan operation plan creation method and operation plan creation device
Rabbanifar et al. Stochastic multi‐objective tie‐line power flow and frequency control in market clearing of multi‐area electricity markets considering power system security
CN101330210B (en) Device and method for formulating employing plan of power station
JP4296140B2 (en) Plant optimum operation support system and method, program
US11921482B2 (en) Method of controlling a microgrid, power management system, and energy management system
JP2019213374A (en) Power generation plan formulation device
Andrianesis et al. Medium-Term Unit Commitment in a pool market
JP2023523979A (en) Power grid resource allocation

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20080328

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20091112

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20091124

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20100121

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20100727

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20100809

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130903

Year of fee payment: 3

R151 Written notification of patent or utility model registration

Ref document number: 4577232

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130903

Year of fee payment: 3

LAPS Cancellation because of no payment of annual fees