JPH11259450A - Optimal output deciding method and device therefor - Google Patents
Optimal output deciding method and device thereforInfo
- Publication number
- JPH11259450A JPH11259450A JP5643098A JP5643098A JPH11259450A JP H11259450 A JPH11259450 A JP H11259450A JP 5643098 A JP5643098 A JP 5643098A JP 5643098 A JP5643098 A JP 5643098A JP H11259450 A JPH11259450 A JP H11259450A
- Authority
- JP
- Japan
- Prior art keywords
- constraint
- value
- input
- output
- area
- 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.)
- Pending
Links
Classifications
-
- Y04S10/54—
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
- Feedback Control In General (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
Description
【０００１】[0001]
【発明の属する技術分野】本発明は、設備の運転計画，
制御システムにおいてその変数である生産量，流通量を
決定するための方法及び装置に関する。TECHNICAL FIELD The present invention relates to an operation plan of equipment,
The present invention relates to a method and an apparatus for determining production variables and distribution volumes which are variables in a control system.
【０００２】[0002]
【従来の技術】発電機の出力を決定する方法として、平
成７年電気学会電力系統技術研究会ＰＥ−９５−３７
「ＭＡＲＫＡＬ型エネルギーモデルへの主双対アフィン
変換法の適用検討」や電気学会電力・エネルギー部門論
Ｖol.１１６−Ｂ,Ｎo.２，'９６「電力系統最適化問題
への内点法二次計画アルゴリズムの適用」で提案されて
いるように、内点法を適用して解くとき、変数の線形和
の制約式を新たな変数とし、この新たな変数に対する上
下限制約を加えることで、変数の上下限制約以外の制約
条件を等式制約で表現している。この定式化により内点
法を用いて最適な発電機出力計算している。ここで、従
来の主双対内点法を用いた発電機出力を決定する方法に
ついて説明しておく。2. Description of the Related Art As a method of determining the output of a generator, the Institute of Electrical Engineers of Japan, 1995, PE-95-37.
"Examination of application of primary dual affine transformation method to MARKAL type energy model" and IEEJ Power & Energy Division, Vol. When applying the interior point method to solve, as shown in `` Applying the algorithm '', the constraint equation of the linear sum of the variables is set as a new variable, and by adding upper and lower bound constraints on this new variable, the The constraint conditions other than the upper and lower limit constraints are expressed by equation constraints. With this formulation, the optimum generator output is calculated using the interior point method. Here, a method of determining the generator output using the conventional primary dual interior point method will be described.
【０００３】目的関数が線形、又は二次式で表され、制
約条件が線形となる問題で目的関数を最小化するときの
標準形は一般に、[0003] In a problem where the objective function is represented by a linear or quadratic expression and the constraint condition is linear, the standard form for minimizing the objective function is generally
【０００４】[0004]
【数１】 (Equation 1)
【０００５】で表現される。ここで、ｘは変数のベクト
ルであり、例えば発電機の出力である。ｃは変数ｘによ
り生じる線形のコスト係数ベクトルである。Ｑは変数ｘ
により生じる二次のコスト係数マトリックスである。発
電機出力によるコストは、一般に発電機出力ｘの二乗に
比例するので、この係数ベクトルをｑとすると、マトリ
ックスＱは、[0005] Here, x is a vector of variables, for example, the output of a generator. c is a linear cost coefficient vector caused by the variable x. Q is the variable x
Is a quadratic cost coefficient matrix generated by Since the cost of the generator output is generally proportional to the square of the generator output x, if this coefficient vector is q, the matrix Q is
【０００６】[0006]
【数２】 (Equation 2)
【０００７】で表される。マトリックスＡは変数ｘの出
力上下限を除いた制約の数をｍ，変数の数をｎとする
と、ｍ×ｎ(ｍ＜ｎ)のマトリックスとなる。変数ｘの下
限をｘs，上限をｘu とする。但し、ｘs は０以上とす
る。この不等式制約はスラック変数ｓ，ｕを導入するこ
とで、[0007] The matrix A is an m × n (m <n) matrix, where m is the number of constraints excluding the upper and lower limits of the output of the variable x, and n is the number of variables. The lower limit of the variable x is x s and the upper limit is x u . However, xs is set to 0 or more. This inequality constraint introduces slack variables s and u,
【０００８】[0008]
【数３】 (Equation 3)
【０００９】と表現できる。変数ｓ，ｕは共に０以上と
する。キューン・カッター条件，相補性定理等により、
数１の最適解（ｘ*，ｙ*，ｚ*，ｗ*）は、Can be expressed as The variables s and u are both 0 or more. By the Kuhn-cutter condition, the complementarity theorem, etc.,
The optimal solution (x *, y *, z *, w *) of Equation 1 is
【００１０】[0010]
【数４】 (Equation 4)
【００１１】でμ＝０を満たすことが必要十分条件であ
る。この数４にニュートン法を適用した、It is a necessary and sufficient condition that μ = 0 is satisfied. Applying Newton's method to this equation 4,
【００１２】[0012]
【数５】 (Equation 5)
【００１３】において、μを０に近づけながら、（ｄ
ｘ，ｄｙ，ｄｚ，ｄｗ）だけ（ｘ，ｙ，ｚ，ｗ）を更新
することを繰り返すことにより、最適解（ｘ*，ｙ*，ｚ
*，ｗ*）に到達することができる。数５でマトリックス
Ａはｍ×ｎ(ｍ＜ｎ)のマトリックスであり、正方行列で
はないので逆行列は存在しないため、この式の２番目の
式からｄｘを求めることはできない。よって、この式
で、In (d), while μ approaches 0, (d)
By repeating updating (x, y, z, w) by x, dy, dz, dw), the optimal solution (x *, y *, z) is obtained.
*, W *). In Equation 5, the matrix A is an m × n (m <n) matrix, and since it is not a square matrix, there is no inverse matrix. Therefore, dx cannot be obtained from the second expression of this expression. Thus, in this equation:
【００１４】[0014]
【数６】 (Equation 6)
【００１５】を代入して変数ｄｘ，ｄｙ以外を除くと、And excluding variables other than variables dx and dy,
【００１６】[0016]
【数７】 (Equation 7)
【００１７】となる。この式から、## EQU1 ## From this equation,
【００１８】[0018]
【数８】 (Equation 8)
【００１９】により、最初にｄｙを計算することができ
る。このｄｙを用いてｄｘを、ｄｘを用いてｄｚ，ｄｗ
を求めることができる。Thus, dy can be calculated first. Dx is obtained by using this dy, and dz and dw are obtained by using dx.
Can be requested.
【００２０】以上の主双対内点法のアルゴリズムの処理
ステップを図２に示す。処理1001では最大反復回数ｋma
x を設定する。処理１００２では、ｋ＝０として、初期
内点となるようにｙ(ｋ)以外は正の値となる任意の初期
点(ｘ(ｋ)，ｓ(ｋ)，ｕ(ｋ),ｙ(ｋ)，ｚ(ｋ)，ｗ(ｋ)）
を設定する。処理１００３では、最適解と現在得られて
いる解との最大偏差の指標である双対ギャップｇ(ｋ)を
計算する。処理1004では、ｇ(ｋ)と双対ギャップの収束
定数εｇと比較し、εｇより小さいならば終了する。処
理１００５では、バリアパラメータμ(ｋ)を計算する。
処理１００６では、現在の解と、μ(ｋ)を元に、数６に
より、Ｄ(ｋ)，ｒ(ｋ)，ｈ(ｋ)を計算する。処理１００
７では、数７，数８から最初にｄｙが求まり、ｄｙから
ｄｘを計算する。さらに、数５，数６にｄｘを代入し
て、ｄｓ，ｄｕ，ｄｚおよびｄｗを計算する。処理１０
０８では、ｓ，ｕ，ｚおよびｗが正でαが最大となるス
テップサイズαｐとαｄを計算する。λは０から１の間
の定数である。処理１００９では、現在の解とステップ
サイズを元に変数を更新する。処理１０１０では最大反
復回数ｋmax よりｋが大きいとき終了する。そうでなけ
れば、処理１００３へ移る。以上を繰り返すことによ
り、最適解を得ることができる。FIG. 2 shows the processing steps of the above algorithm of the primary dual interior point method. In processing 1001, the maximum number of iterations kma
Set x. In the process 1002, k = 0, and any initial point (x (k), s (k), u (k), y (k) having a positive value except y (k) so as to become an initial interior point. ), Z (k), w (k))
Set. In processing 1003, a dual gap g (k), which is an index of the maximum deviation between the optimal solution and the currently obtained solution, is calculated. In the process 1004, g (k) is compared with the convergence constant εg of the dual gap. In processing 1005, a barrier parameter μ (k) is calculated.
In the process 1006, D (k), r (k), and h (k) are calculated according to Equation 6 based on the current solution and μ (k). Processing 100
In 7, dy is obtained first from Equations 7 and 8, and dx is calculated from dy. Furthermore, ds, du, dz, and dw are calculated by substituting dx into Equations 5 and 6. Processing 10
In step 08, the step sizes αp and αd at which s, u, z and w are positive and α is maximum are calculated. λ is a constant between 0 and 1. In the process 1009, a variable is updated based on the current solution and the step size. The process 1010 ends when k is larger than the maximum number of iterations kmax. Otherwise, the process moves to processing 1003. By repeating the above, an optimal solution can be obtained.
【００２１】また、内点法を用いた例として特開平9−2
52542 号公報では、発電機の出力変化速度の制約を考慮
して、１つの時間断面の発電機出力を決定し、平成７年
電気学会電力系統技術研究会ＰＥ−９５−３７「ＭＡＲ
ＫＡＬ型エネルギーモデルへの主双対アフィン変換法の
適用検討」では期と期を関係づける複数期に跨る変数の
制約を考慮してエネルギー源の開発を決定していた。An example using the interior point method is disclosed in
In Japanese Patent No. 52542, the generator output in one time section is determined in consideration of the constraint on the output change speed of the generator, and the Institute of Electrical Engineers of Japan, 1995, PE-95-37, "MAR
In "Examination of application of primary dual affine transformation method to KAL-type energy model", the development of energy sources was determined in consideration of the constraints of variables over multiple periods that relate periods.
【００２２】[0022]
【発明が解決しようとする課題】平成７年電気学会電力
系統技術研究会ＰＥ−９５−３７「ＭＡＲＫＡＬ型エネ
ルギーモデルへの主双対アフィン変換法の適用検討」や
電気学会電力・エネルギー部門論Ｖol.１１６−Ｂ，Ｎ
o.２，'９６「電力系統最適化問題への内点法二次計画
アルゴリズムの適用」の従来技術では内点法を適用する
に当たって、以下のように制約式を変換している。[Problems to be Solved by the Invention] In 1995, the Institute of Electrical Engineers of Japan, Power System Technical Research Group PE-95-37, "Examination of application of primary dual affine transformation method to MARKAL type energy model" and IEEJ Power and Energy Division, Vol. 116-B, N
o.2, '96 In the prior art of "application of the interior point method quadratic programming algorithm to the power system optimization problem", in applying the interior point method, the constraint equation is converted as follows.
【００２３】[0023]
【数９】 (Equation 9)
【００２４】の等式制約は発電機出力を表す変数ｘによ
りそのまま表現する。The equation constraint is directly expressed by a variable x representing the generator output.
【００２５】[0025]
【数１０】 (Equation 10)
【００２６】の下限のみの制約は線形の制約式を新たな
変数ｆ1 とすることにより、新たな変数ｆを含めた等式
となる。内点法では、変数は０以上であることから、こ
の式の右辺の下限値は０以上であるとする。上限制約の
みのときも同様に表現できる。The constraint of only the lower limit becomes an equation including the new variable f by using a linear constraint equation as the new variable f 1 . In the interior point method, since the variable is 0 or more, it is assumed that the lower limit of the right side of this equation is 0 or more. The same can be said for only the upper limit constraint.
【００２７】[0027]
【数１１】 [Equation 11]
【００２８】の上下限制約は変数の線形和を新たな１つ
の変数ｆ2 として等式制約に変換し、この新たな変数は
下限値以上，上限値以下とする。変数ｘがｎ１次元ベク
トル，等式制約の数がｎ２，下限のみまたは上限のみの
制約の数がｎ３，上下限制約の数がｎ４とすると、制約
マトリックスＡは、The upper / lower limit constraint converts a linear sum of variables into an equality constraint as a new variable f 2 , and the new variable has a value equal to or greater than a lower limit and equal to or less than an upper limit. If the variable x is an n1-dimensional vector, the number of equality constraints is n2, the number of constraints with only lower and upper bounds is n3, and the number of upper and lower bounds is n4, the constraint matrix A is
【００２９】[0029]
【数１２】 (Equation 12)
【００３０】のように、（ｎ２＋ｎ３＋ｎ４)×(ｎ１＋
ｎ３＋ｎ４）の行列となる。数８の左辺のマトリックス
の積はＡとその転置行列の積があるため、この左辺のマ
トリックスは（ｎ１＋ｎ３＋ｎ４)×(ｎ１＋ｎ３＋ｎ
４）の正方行列になる。制約の数が増加とともに、行列
が大きくなる。この結果、記憶する要素の数が多くなる
こと、要素の数が多いため全ての要素を計算する時間が
増大することになる。As described above, (n2 + n3 + n4) × (n1 +
n3 + n4). Since the product of the matrix on the left side of Equation 8 is the product of A and its transposed matrix, the matrix on the left side is (n1 + n3 + n4) × (n1 + n3 + n
4) is a square matrix. As the number of constraints increases, the matrix grows. As a result, the number of elements to be stored increases, and the time required to calculate all elements increases because of the large number of elements.
【００３１】そこで、本発明では、等式制約を上限また
は下限制約に置換する、または等式制約および上限また
は下限のみの制約式を上下限制約とすることにより、制
約マトリックスＡをブロック化した行列の各行に単位行
列を配置することができ、従来より解を更新するための
行列の計算時間を短縮することができる出力決定方法と
該方法を用いた作成装置および記憶媒体とを提供するこ
とを目的とする。Therefore, in the present invention, the equality constraint is replaced with an upper limit or a lower limit, or the equality constraint and a constraint expression having only the upper or lower limit are set as the upper and lower limits, thereby forming a block of the constraint matrix A. It is an object of the present invention to provide an output determining method capable of arranging a unit matrix in each row of the matrix and shortening the calculation time of a matrix for updating a solution, and a creating apparatus and a storage medium using the method. Aim.
【００３２】[0032]
【課題を解決するための手段】等式制約の例として、需
給バランス制約がある。目的関数を燃料費最小化とする
ならば、発電機の出力を必要最小限とすることになる。
よって、この特性を利用すると、等式制約は不等式制約
に置き換えることができる。この結果、マトリックスＡ
は、An example of an equation constraint is a supply-demand balance constraint. If the objective function is to minimize fuel cost, the output of the generator will be minimized.
Therefore, using this property, the equality constraint can be replaced with the inequality constraint. As a result, matrix A
Is
【００３３】[0033]
【数１３】 (Equation 13)
【００３４】の様に表現することができる。行列Ａをブ
ロック化すると、It can be expressed as follows. When the matrix A is blocked,
【００３５】[0035]
【数１４】 [Equation 14]
【００３６】のように表現できる。各ブロックの行方向
に単位行列Ｉがあることを利用して、この行列を数７に
代入すると、ｄｘは、Can be expressed as follows. Substituting this matrix into Equation 7 using the fact that there is a unit matrix I in the row direction of each block, dx becomes
【００３７】[0037]
【数１５】 （Ｑ＋Ａ１TＤ１Ａ１＋Ａ２TＤ２Ａ２＋Ａ３TＤ３Ａ３)ｄｘ＝ｇ …（数１５） の形式により求めることができる。ｄｘにかかる行列は
行列Ｑ，Ａ１とその転値行列の積，Ａ２とその転値行列
の積およびＡ３とその転値行列の積の和となる。各ブロ
ックの行列は単位行列部分がないので、ｄｘにかかる行
列を求めるために必要な計算回数を減少できる。さら
に、各正方行列の大きさが小さくなることから、行列の
要素の値を保持するために必要なメモリを従来より低減
することができる。(Q + A1 T D1A1 + A2 T D2A2 + A3 T D3A3) dx = g (Expression 15) The matrix related to dx is the sum of the products of the matrices Q and A1 and its transposition matrix, the product of A2 and its transposition matrix, and the product of A3 and its transposition matrix. Since the matrix of each block does not have a unit matrix portion, the number of calculations required to find a matrix related to dx can be reduced. Further, since the size of each square matrix is reduced, the memory required to hold the values of the elements of the matrix can be reduced as compared with the conventional case.
【００３８】同様に、制約マトリックスＡの各ブロック
の行方向に単位行列を配置する方法として、Similarly, as a method of arranging the unit matrix in the row direction of each block of the constraint matrix A,
【００３９】[0039]
【数１６】 Ａ＝［Ａ０ １］ …（数１６） のようにブロックＡ０と単位行列で表す方法、全てを上
下限制約とした、A = [A0 1] (Equation 16) A method expressed by a block A0 and a unit matrix as shown in Expression 16
【００４０】[0040]
【数１７】 [Equation 17]
【００４１】の方法がある。前者ではｄｘにかかる行列
は行列ＱとＡ０とその転置行列の積の和となり、後者で
は行列ＱとＡ４とその転置行列の積の和となる。いずれ
の方法においても、従来より、計算時間および行列の要
素を保持するメモリを低減することができる。The method is as follows. In the former case, the matrix related to dx is the sum of the products of the matrices Q and A0 and its transposed matrix, and in the latter case it is the sum of the products of the matrices Q and A4 and its transposed matrix. In any of the methods, the calculation time and the memory for holding the elements of the matrix can be reduced.
【００４２】特開平9−252542 号公報では、需給バラン
ス制約，発電機出力の出力上下限制約，発電機の出力変
化速度の制約を考慮して、現在の発電機出力から次の１
時点の発電機出力を決定していた。電力系統をより最適
に運用する上では、複数の時点を同時に考慮して、最適
な発電機出力を決定する必要がある。In Japanese Patent Application Laid-Open No. 9-252542, the following 1
The generator output at the time was determined. In order to operate the electric power system more optimally, it is necessary to determine an optimal generator output by simultaneously considering a plurality of time points.
【００４３】[0043]
【発明の実施の形態】以下、本発明を発電機出力対象と
した出力決定方法に適用した場合の実施例について、図
面を用いて説明する。DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS An embodiment in which the present invention is applied to an output determination method for a generator output will be described below with reference to the drawings.
【００４４】Ａ．装置構成 （１）ハードウェア構成 本実施例の出力決定装置１は、図３に示すように、中央
演算処理装置(ＣＰＵ)２と、主記憶装置３と、入出力装
置４と、外部記憶装置５とを備える。A. 3. Device Configuration (1) Hardware Configuration As shown in FIG. 3, the output determination device 1 of the present embodiment includes a central processing unit (CPU) 2, a main storage device 3, an input / output device 4, an external storage device 5 is provided.
【００４５】入出力装置４は、入力装置６（図４に図
示）として、キーボードおよびマウスを備え、出力装置
として、ＣＲＴ（cathode ray tubu）を備える表示装置
７（図４に図示）を備える。なお、入出力装置４とし
て、これらの代わりに、あるいは、これらと併用して、
ポインティングディバイス，タッチセンサ等の入力手段
や、液晶表示装置，プリンタ，スピーカ等の出力手段を
備えてもよい。The input / output device 4 includes a keyboard and a mouse as an input device 6 (shown in FIG. 4), and a display device 7 (shown in FIG. 4) having a CRT (cathode ray tubu) as an output device. In addition, as the input / output device 4, instead of these or in combination with these,
Input means such as a pointing device and a touch sensor and output means such as a liquid crystal display device, a printer and a speaker may be provided.
【００４６】外部記憶装置５としては、ハードディスク
装置，フロッピィディスク装置，ＣＤ−ＲＯＭ(compact
disc−read only memory）装置，ＤＡＴ（digital vid
eotape）装置，ＲＡＭ(randam access memory)装置，Ｄ
ＶＤ(digital video disc)装置，不揮発性メモリ等を用
いることができる。なお、本実施例では、外部記憶装置
５として、データベース８（図４に図示）などを保持す
るための大容量記憶装置と、処理プログラムなどを保持
する記憶媒体と、該記憶媒体に保持された情報を読み取
るための読取装置９（図４に図示）とを用いるが、本発
明はこれらには限られず、一つの外部記憶装置にデータ
ベースと処理プログラムとが両方保持されているように
してもよい。記憶媒体としては、フロッピィディスク，
ＣＤ−ＲＯＭ，磁気テープ，光ディスク，光磁気ディス
ク，ＤＡＴ，ＲＡＭ，ＤＶＤ，不揮発性メモリ等を用い
ることができる。As the external storage device 5, a hard disk device, a floppy disk device, a CD-ROM (compact
Disc-read only memory (DAT), digital vid
eotape) device, RAM (randam access memory) device, D
A VD (digital video disc) device, a nonvolatile memory, or the like can be used. In this embodiment, as the external storage device 5, a large-capacity storage device for holding a database 8 (shown in FIG. 4) and the like, a storage medium for holding processing programs and the like, and a storage medium for holding the processing program and the like. A reading device 9 (shown in FIG. 4) for reading information is used, but the present invention is not limited thereto, and one external storage device may hold both a database and a processing program. . As storage media, floppy disks,
A CD-ROM, a magnetic tape, an optical disk, a magneto-optical disk, a DAT, a RAM, a DVD, a nonvolatile memory, or the like can be used.
【００４７】（２）機能構成 本発明の出力決定装置１の機能構成を図４に示す。本実
施例の出力決定装置１は、上述の入力装置６，表示装置
７，データベース８および読取装置９と、データ処理部
１７とを備える。なお、本実施例では、データ処理部１
７は、あらかじめ外部記憶装置５の記憶媒体に保持さ
れ、読取装置９を介して主記憶装置３に読み込まれたプ
ログラムを、ＣＰＵ２が実行することにより実現される
が、本発明はこのようなプログラムされた汎用プロセッ
サによるものに限られるわけではない。例えば、本発明
の各ステップを実行するハードワイヤードロジックを含
む特定のハードウェア装置との組合せによって、データ
処理部１７を実現しても良い。(2) Functional Configuration FIG. 4 shows the functional configuration of the output determination device 1 of the present invention. The output determination device 1 of the present embodiment includes the above-described input device 6, display device 7, database 8, and reading device 9, and a data processing unit 17. In this embodiment, the data processing unit 1
7 is realized by the CPU 2 executing a program stored in the storage medium of the external storage device 5 in advance and read into the main storage device 3 via the reading device 9. However, the present invention is not limited to the general-purpose processor. For example, the data processing unit 17 may be realized by a combination with a specific hardware device including hard wired logic that executes each step of the present invention.
【００４８】入力装置６は表示装置７に表示された選択
肢の選択、データの入力などを受け付け、データ処理部
１７に送る手段である。表示装置７は、入力装置６から
送られたデータを表示する手段である。データ処理部１
７は入力装置６通知されたデータと、データベース８か
ら読み込んだデータと、読取装置９から読み込んだ処理
プログラムと、電力システム１６から通知されたデータ
とを基に、最適な出力を決定する手段である。The input device 6 is a means for receiving selection of an option displayed on the display device 7, input of data, and the like, and sending the data to the data processing section 17. The display device 7 is a means for displaying data sent from the input device 6. Data processing unit 1
Numeral 7 is a means for determining an optimum output based on the data notified from the input device 6, the data read from the database 8, the processing program read from the reading device 9, and the data notified from the power system 16. is there.
【００４９】なお、本実施例では、データ処理部１７に
よる処理結果は、表示装置７に送られて表示されるとと
もに、データベース１４に格納される。また、電力シス
テム１０の運用・計画システム，制御システムが設定し
た条件での出力決定が要求された場合には、データ処理
部１７の作成した出力および評価結果は電力システム１
０へも通知される。この通知を受けた電力システム１０
は、通知された出力をもとに電力系統１１に発電機出力
等の制御信号を出して、発電機を制御するとともに、電
力系統１１から、該系統の情報（発電機出力，連系線潮
流，送電線潮流，ダムの池水位，ＬＮＧタンクレベル
等）を取り込み、その取り込んだデータと、計画システ
ム（図示せず）において作成したデータとを、内部デー
タベース（図示せず）に格納する。In this embodiment, the processing result by the data processing unit 17 is sent to the display device 7 for display and stored in the database 14. When an output determination is required under the conditions set by the operation / planning system and the control system of the power system 10, the output created by the data processing unit 17 and the evaluation result are stored in the power system 1.
0 is also notified. Power system 10 receiving this notification
Sends a control signal, such as a generator output, to the power system 11 based on the notified output to control the generator, and from the power system 11, information on the system (generator output, interconnection line power flow) , Power line power flow, dam pond water level, LNG tank level, etc.) and store the captured data and data created by the planning system (not shown) in an internal database (not shown).
【００５０】本実施例のデータ処理部１７は、コントロ
ール部１２，データ設定部１３，データ読込部１４，出
力計算部１５および表示データ作成部１６を備える。ま
た、データ処理部１７は、通信手段１８を介して外部の
電力システム１０に接続されており、この電力システム
は、実際の送発変電設備である電力系統１１に接続され
ている。なお、電力システム１０はこの電力系統１１の
計画・運用のための支援システム、電力系統１１、およ
び電力系統１の状態などを示す情報を保持するためのデ
ータベースを備える。電力系統１１の状態は、リレーや
センサなどにより検出され、通信手段１９を介して電力
システム１０に通知され、該システム１０のデータベー
ス（図示せず）に格納される。The data processing section 17 of this embodiment includes a control section 12, a data setting section 13, a data reading section 14, an output calculation section 15, and a display data creation section 16. The data processing unit 17 is connected to an external power system 10 via a communication unit 18, and the power system is connected to a power system 11 which is an actual transmission / reception substation. The power system 10 includes a support system for planning and operation of the power system 11, a database for holding information indicating the state of the power system 11, the state of the power system 1, and the like. The state of the power system 11 is detected by a relay, a sensor, or the like, notified to the power system 10 via the communication unit 19, and stored in a database (not shown) of the system 10.
【００５１】コントロール部１２は、電力システム１
０、および、各処理部１３〜１６の間のデータや処理プ
ログラムなどの授受を円滑に行うためのデータの加工・
処理を行い、その授受をコントロールして、全体の処理
を正常に動作させるための手段である。The control unit 12 controls the power system 1
0, and data processing / smoothing for smoothly transferring data and processing programs between the processing units 13 to 16.
This is a means for performing processing, controlling the transfer thereof, and operating the entire processing normally.
【００５２】データ設定部１３は、コントロール部１２
を介して、既にデータベース８および／または読取装置
９の記憶媒体に保持されている、発電機出力を決定し使
用する電力系統のデータ，発電機特性，送変電設備特
性，目的関数等の条件や、処理プログラムなどを読み込
む手段である。また、データ設定部１４は、読み込んだ
条件を変更し、新規の条件を作成するとともに、コント
ロール部１２を介して、作成した条件をデータベース１
４に格納する。The data setting unit 13 includes a control unit 12
And the conditions of the power system data, generator characteristics, transmission / transformation facility characteristics, objective function, etc., which are already stored in the database 8 and / or the storage medium of the reading device 9 and which determine and use the generator output. , And a processing program. Further, the data setting unit 14 changes the read condition, creates a new condition, and stores the created condition into the database 1 via the control unit 12.
4 is stored.
【００５３】データ読込部１４は、コントロール部１２
を介して、データベース８に保持されているデータを読
み込み、コントロール部１２を介して、データ設定部１
４または出力計算部１５に通知する。なお、本実施例で
は、特に説明しない限り、データベース８からのデータ
の読み込みは、このデータ読込部１４を介して行われ
る。The data reading section 14 includes the control section 12
The data stored in the database 8 is read via the control unit 12, and the data setting unit 1 is read via the control unit 12.
4 or the output calculation unit 15 is notified. In the present embodiment, reading of data from the database 8 is performed via the data reading unit 14 unless otherwise specified.
【００５４】出力計算部１５は、主送対内点法を用いて
目的関数が最も良い発電機出力，電力融通量，託送量，
ＩＰＰ（独立電気事業者）からの電力購入量等をを自動
的に計算する。The output calculation unit 15 calculates the generator output, the power interchange amount, the contracted amount,
Automatically calculate the amount of power purchased from an IPP (Independent Electricity Utility).
【００５５】表示データ作成部１６は、出力計算部１５
が計算した出力や目的関数の推移，制約充足度，入力装
置６からのデータ入力の支援のための情報等を、表示装
置７表示する。なお、本実施例では、特に説明しない限
り、表示装置７への出力は、この表示データ作成部１６
を介して行われる。The display data creation unit 16 includes the output calculation unit 15
The display device 7 displays the calculated output, the transition of the objective function, the degree of constraint satisfaction, information for supporting data input from the input device 6, and the like. In this embodiment, the output to the display device 7 is output to the display data creation unit 16 unless otherwise specified.
Done through.
【００５６】（３）データベース８のデータ構造 データベース８は図５に示すように、計画対象期間の記
憶領域１００と、目的関数の記憶領域１０１と、制約条
件の記憶領域１０２と、主双対内点法のパラメータ記憶
領域１０３と、計算結果記憶領域１０４とを備える。(3) Data Structure of Database 8 As shown in FIG. 5, the database 8 has a storage area 100 for the planning target period, a storage area 101 for the objective function, a storage area 102 for the constraint condition, and a primary dual interior point. It has a parameter storage area 103 for the modulus and a calculation result storage area 104.
【００５７】計画対象期間記憶領域は、出力を決定する
年月日と開始時刻、および終了時刻を保持するための記
憶領域である。The planning period storage area is a storage area for holding the date, start time, and end time for determining the output.
【００５８】目的関数記憶領域１０１は時間毎の目的関
数の情報を保持するための記憶領域であり、発電機等の
設備特性から出力に応じて増減するコスト関数の係数を
保持するための発電機等設備特性記憶領域１０１１と、
制約違反が生じたとき、制約と設備とその違反量に応じ
てペナルティコストを加算するときの重み係数を保持す
るためのペナルティ重み係数記憶領域１０１２とを備え
る。The objective function storage area 101 is a storage area for storing information on the objective function for each time. An equal facility characteristic storage area 1011;
A penalty weight coefficient storage area 1012 is provided for holding a weight coefficient when a penalty cost is added according to a constraint, a facility, and a violation amount when a constraint violation occurs.
【００５９】制約条件記憶領域１０２は、出力を決定す
るときに考慮する時間毎の制約条件を保持するための記
憶領域であり、想定された需要を保持するための需要記
憶領域１０２１，制約条件を満たす解がないときの出力
の初期値，需給バランスを維持できるチェックするため
の最大供給力，最小供給力を保持するための初期供給力
記憶領域１０２２と、送電ロスを考慮したとき需給バラ
ンスを維持できるかどうかチェックするための送電ロス
率を保持するための初期送電ロス率１０２３と、発電機
等の起動停止状態を表す数値を保持する発電機等起動停
止状態記憶領域１０２４，計画対象期間全体で共通であ
る計画対象期間で発電量と水力，揚水の出発水位を保持
する供給力積分制約記憶領域１０２５と、潮流計算で使
用するために送電線と変圧器をブランチで表現し、母線
をノードで表現した系統構成でノード間を接続するブラ
ンチの使用／未使用（停止）状態を数字で表した値と、
各発電機等の供給力が接続するノード番号と各ノードに
かかる負荷の大きさを保持するための系統接続状態記憶
領域１０２６と、系統設備のブランチ毎のインピーダン
ス、潮流値に応じた送電ロスを計算するための送電ロス
係数，潮流制約の上下限値を保持するための系統設備運
用特性記憶領域１０２７と、発電機の出力上下限，発電
か負荷かを番号で示す供給力特性，発電機の種類を番号
で示す設備特性を保持するための発電機運用特性記憶領
域１０２８とを備える。計算結果記憶領域１０４は、内
点法を用いて問題を解く途中ででてきた計算結果である
変数の値と最終結果である変数の値とを保持するための
記憶領域である。The constraint storage area 102 is a storage area for retaining constraint conditions for each time to be considered when determining an output, and stores a demand storage area 1021 for retaining an assumed demand. The initial value of output when there is no solution to be satisfied, the maximum supply power for checking that the supply and demand balance can be maintained, the initial supply power storage area 1022 for keeping the minimum supply power, and the supply and demand balance when the power transmission loss is considered. The initial power transmission loss rate 1023 for holding the power transmission loss rate for checking whether it can be performed, and the generator etc. start / stop state storage area 1024 holding the numerical value indicating the start / stop state of the generator / etc. A power supply integral constraint storage area 1025 that holds the power generation amount, the hydropower, and the starting water level of the pumped water in the common planning period, and the power transmission for use in the tidal flow calculation A value representing transformer represented by branch, used / unused branches connecting the nodes in the system configuration representing the bus at the node (stop) state numbers and,
A system connection state storage area 1026 for holding a node number to which the supply power of each generator or the like is connected and a magnitude of a load applied to each node, and an impedance for each branch of the system equipment and a power transmission loss according to a power flow value. A power transmission loss coefficient for calculation, a system facility operation characteristic storage area 1027 for holding upper and lower limits of a power flow constraint, an upper and lower limit of a generator output, a supply power characteristic indicating a power generation or a load by a number, and a A generator operating characteristic storage area 1028 for holding equipment characteristics indicating the type by number. The calculation result storage area 104 is a storage area for holding the value of a variable that is a calculation result and the value of a variable that is a final result obtained while solving a problem using the interior point method.
【００６０】Ｂ．処理の流れ 本実施例の出力決定装置１は、電力系統の発電機出力等
の供給力の値を計算するための、各種データを設定・編
集する機能と、設定・編集されたデータを用いて電力系
統の発電機出力等の供給力を自動的に計算する機能とを
備える。B. Process Flow The output determining apparatus 1 of the present embodiment uses a function of setting and editing various data for calculating a value of a supply power such as a generator output of a power system, and the set and edited data. A function of automatically calculating a supply power such as a generator output of the power system.
【００６１】本実施例の出力決定装置１は、まず表示装
置７に、図６に示す処理選択メニュー画面２０を表示
し、入力装置６を介して選択およびデータの入力・設定
を受け付ける。処理選択メニュー画面２０は、電力系統
の計画作成条件設定(項目１)，出力計算（項目２）およ
び結果の表示（項目３）からなり、さらに計画作成条件
（項目１）は計画対象期間の入力，設定（項目ａ），目
的関数と考慮する制約条件の選択（項目ｂ），需要デー
タの入力，設定（項目ｃ），発電機等起動停止計画の入
力，設定（項目ｄ），初期出力と供給力チェックの入
力，設定(項目ｅ)，発電機等設備特性の入力，設定（項
目ｆ），供給力積分特性の入力，設定（項目ｇ），系統
接続状態の入力，設定（項目ｈ），系統設備の運用特性
の入力，設定（項目ｉ），収束判定条件の入力，設定
（項目ｊ），内点法制御変数の入力，設定（項目ｋ），
内点法の初期値の入力，設定（項目ｌ），等式制約の許
容偏差の入力，設定（項目ｍ）からなり、総計１５項目
のメニューを備える。ここで、入力装置６は、タッチ入
力，ポインティングデバイス，マウスなどによる項目毎
の選択領域の選択入力，キーボードによる項目番号また
は記号入力等による、表示された１５項目のいずれかの
選択を受け付け、選択された項目の処理を実行する。以
下に、各項目の処理について説明する。The output determining device 1 of this embodiment first displays the process selection menu screen 20 shown in FIG. 6 on the display device 7 and accepts selection and input / setting of data via the input device 6. The process selection menu screen 20 includes a power system plan creation condition setting (item 1), an output calculation (item 2), and a result display (item 3). , Setting (item a), selection of constraint conditions to be considered as the objective function (item b), input of demand data, setting (item c), input of start / stop plan of generators, etc., setting (item d), initial output Input and setting of supply power check (item e), input and setting of equipment such as generator (item f), input and setting of supply power integration characteristics (item g), input and setting of system connection status (item h) , Input and setting of operation characteristics of system equipment (item i), input and setting of convergence judgment condition (item j), input and setting of interior point method control variables (item k),
The menu consists of inputting and setting of initial values of the interior point method (item 1), inputting and setting of allowable deviation of equation constraints (item m), and has a total of 15 items. Here, the input device 6 receives and selects any one of the displayed 15 items by touch input, selection input of a selection area for each item by a pointing device, a mouse, or the like, and item number or symbol input by a keyboard. Execute the processing of the selected item. The processing of each item will be described below.
【００６２】（１）計画作成条件の設定 入力装置６を介してメニュー項目１のなか ａ．計画対象期間の入力，設定 処理選択メニュー画面２０において、計画対象期間の入
力，設定（項目ａ）が入力装置６を介して選択される
と、データ設定部１３は処理対象の電力系統のデータベ
ース８の計画対象期間記憶領域１００から年，月，日，
開始時刻および終了時刻を読み込む。計画対象期間記憶
領域１００にデータが格納されているときはその値を格
納されていないときは、デフォルト値である０が入力さ
れた設定画面３０を表示される。ここで、表示される計
算対象期間の入力，設定画面３０を図７に示す。この画
面３０は、計画対象期間の年，月，日，開始時刻および
終了時刻の入力領域３１〜３５を備え、さらにデータ格
納３６と設定終了３７の領域がある。必要に応じて、入
力領域３１〜３５に入力される。(1) Setting of plan creation conditions In menu item 1 via input device 6 a. Input and setting of target period When input and setting (item a) of the target period are selected via the input device 6 on the process selection menu screen 20, the data setting unit 13 sets the database 8 of the power system to be processed. Year, month, day,
Read the start time and end time. When data is stored in the planning target period storage area 100, and when the value is not stored, the setting screen 30 in which a default value of 0 is input is displayed. Here, the input and setting screen 30 for the calculation target period displayed is shown in FIG. The screen 30 includes input areas 31 to 35 for the year, month, day, start time, and end time of the plan target period, and further includes a data storage 36 and a setting end 37 area. Input is made to the input areas 31 to 35 as necessary.
【００６３】なお、表示された年，月，日，開始時刻お
よび終了時刻と異なる値を入力するには、入力装置６を
介して入力領域３１〜３５へ数値を入力する。入力領域
３１〜３５の値は「データ格納」３６の領域を選択する
ことにより、計画対象期間領域１００に格納される。
「設定終了」３７の領域が選択されると、データ設定部
１３は、計画対象期間の入力，設定を終了し、処理選択
メニュー画面２０に戻る。In order to input a value different from the displayed year, month, day, start time and end time, a numerical value is input to the input areas 31 to 35 via the input device 6. The values of the input areas 31 to 35 are stored in the planning target period area 100 by selecting an area of “data storage” 36.
When the area of “setting end” 37 is selected, the data setting unit 13 ends the input and setting of the planning target period, and returns to the processing selection menu screen 20.
【００６４】ｂ．目的関数と考慮する制約条件の選択 処理選択メニュー画面２０において、目的関数と考慮す
る制約条件の選択（項目ｂ）が入力装置６を介して選択
されると、データ設定部１３は処理対象のデータベース
８の目的関数記憶領域１０１から目的関数の重み係数，
項目の考慮，制約条件の重み係数，制約の考慮の値を読
み込む。目的関数記憶領域１０１にデータが格納されて
いるときは、格納されている値を、格納されていないと
きはデフォルト値である０が入力された設定画面４０が
表示される。ここで、表示される目的関数と考慮する制
約条件の画面４０を図８に示す。この画面４０は、目的
関数の項目４１を表示し、目的関数の重み係数４２及び
その目的関数の項目を考慮するかどうかのフラグ４３の
入力領域と、目的関数の項目をスクロールさせるための
スクロール４７と、制約条件の項目４４を表示し、制約
条件の重み係数と制約を考慮するかどうかのフラグ４６
の入力領域と、制約条件の項目をスクロールさせるため
のスクロール４８とを備え、さらにデータ格納４９と設
定終了５０の領域がある。スクロール４７，４８を選択
することにより、その選択した方向に応じて、目的関
数，制約条件の項目が上，下へスクロールして、全体の
項目及び値を確認することができる。B. Selection of Constraint Condition to be Considered as Objective Function When selection of a constraint condition to be considered as an objective function (item b) is selected via the input device 6 on the process selection menu screen 20, the data setting unit 13 sets the database to be processed. 8 from the objective function storage area 101,
The values of item consideration, constraint weighting factor, and constraint consideration are read. When data is stored in the objective function storage area 101, the stored screen is displayed, and when not stored, a setting screen 40 in which a default value of 0 is input is displayed. Here, a screen 40 of the objective function and the constraint condition to be considered is shown in FIG. The screen 40 displays an objective function item 41, an input area of a weight coefficient 42 of the objective function and a flag 43 for determining whether to consider the objective function item, and a scroll 47 for scrolling the objective function item. And a constraint condition item 44, and a flag 46 indicating whether to consider the constraint weighting factor and the constraint.
, And a scroll 48 for scrolling the items of the constraint condition. Further, there are areas for data storage 49 and setting end 50. By selecting the scrolls 47 and 48, the items of the objective function and the constraints are scrolled up and down in accordance with the selected direction, so that all items and values can be confirmed.
【００６５】目的関数の項目として、電気を調達するの
に必要なコストを含めた燃料費，送電ロス，ダムからの
水の放流パターン，ダムの溢水量等がある。使用する目
的関数は、各項目に重み係数を掛けた総和である。但
し、「重み係数」４２として有効なのは、「項目の考
慮」４３が１のときのみである。０のときは考慮しな
い。目的関数が電気を調達するのに必要なコストを含め
た燃料費とするときを、Items of the objective function include a fuel cost including a cost necessary to procure electricity, a power transmission loss, a water discharge pattern from a dam, an overflow amount of a dam, and the like. The objective function to be used is a sum total obtained by multiplying each item by a weight coefficient. However, the “weight coefficient” 42 is effective only when the “item consideration” 43 is 1. When it is 0, it is not considered. When the objective function is the fuel cost including the cost required to procure electricity,
【００６６】[0066]
【数１８】 (Equation 18)
【００６７】に示す。Is shown in FIG.
【００６８】制約条件として、需給バランス，出力上下
限，潮流制約，計画対象期間で発電する電力量（燃料消
費制約，ＩＰＰ利用率等は電力量に換算），池水位制
約，出力変化速度，発電所単位での最低出力，計画対象
期間内で放流する水量の変動制約等がある。これらの制
約は３つに分類できる。１つは時間断面毎の制約であ
り、需給バランス，潮流制約がこれに当たる。需給バラ
ンスは、The constraints include supply and demand balance, output upper and lower limits, power flow constraints, the amount of power generated during the planning period (fuel consumption constraints, IPP utilization, etc. are converted to power), pond water level constraints, output change speed, power generation There is a minimum output per unit, a restriction on the amount of water discharged during the target period, and other factors. These constraints can be divided into three categories. One is a constraint for each time section, which corresponds to a supply-demand balance and a tidal current constraint. The supply and demand balance is
【００６９】[0069]
【数１９】 [Equation 19]
【００７０】となる。図２１の系統での潮流制約は、Is obtained. The power flow constraint in the system of FIG.
【００７１】[0071]
【数２０】 (Equation 20)
【００７２】となる。スラック変数を用いると、Is obtained. Using slack variables,
【００７３】[0073]
【数２１】 (Equation 21)
【００７４】となる。スラック変数を導入することによ
り、等式制約にすることができることがわかる。２つめ
の制約は差分制約で、この制約として出力変化速度制約
がある。出力変化速度制約は、Is obtained. It can be seen that by introducing a slack variable, it can be made an equality constraint. The second constraint is a difference constraint, which is an output change speed constraint. The output change rate constraint is
【００７５】[0075]
【数２２】 (Equation 22)
【００７６】で表される。スラック変数を導入すると、Is represented by Introducing the slack variable,
【００７７】[0077]
【数２３】 (Equation 23)
【００７８】のように等式制約変換できる。３つめは積
分制約であり、Equation conversion can be performed as follows. The third is the integral constraint,
【００７９】[0079]
【数２４】 (Equation 24)
【００８０】に示す水位制約がこれにあたる。スラック
変数を導入することにより、The water level constraint shown in FIG. By introducing slack variables,
【００８１】[0081]
【数２５】 (Equation 25)
【００８２】のように、等式制約に変換できる。変数を
考慮する制約条件は、「制約の考慮」４６が１又は２の
ものである。０は考慮しない。「制約の考慮」４６のフ
ラグが１および２の制約を全て満たす解がないことが、
出力計算部１５の処理で分かったときは、「制約の考
慮」４６の値が２の制約のみを考慮して、「制約の考
慮」４６の値が１の制約値は出力計算部１５の処理によ
り緩和した制約値に変更するとともに、制約式に「重み
係数」４５を加算した総和をペナルティコストとして目
的関数に加算する。制約違反のペナルティ加算方法とし
て、違反量の２次式とすることもできる。この重み係数
が大きいほど、ペナルティコストが大きくなるので、目
的関数を最小化するときは、その制約を満たすように出
力を決定することになり、制約の重要度，優先順位を表
す指標となる。As described above, it can be converted into an equality constraint. The constraint condition in which the variable is considered is one in which “consider constraint” 46 is 1 or 2. 0 is not considered. The fact that there is no solution in which the flag of "consider constraint" 46 satisfies all of the constraints of 1 and 2,
When it is found by the processing of the output calculation unit 15, only the constraint having the value of “consider constraint” 46 of 2 is considered, and the constraint value of 1 of “consider constraint” 46 is 1 , And the total sum of the constraint formula and the “weight coefficient” 45 is added to the objective function as a penalty cost. As a penalty addition method for the constraint violation, a quadratic expression of the violation amount may be used. The penalty cost increases as the weight coefficient increases. Therefore, when minimizing the objective function, the output is determined so as to satisfy the constraint, and is an index indicating the importance and priority of the constraint.
【００８３】なお、表示された重み係数４２と、項目の
考慮４３と、重み係数４５と制約の考慮４６と異なる値
を入力するには、入力装置６を介して入力領域４２，４
３，４５，４６へ数値を入力する。入力領域４２，４
３，４５，４６の値は「データ格納」４９の領域を選択
することにより、目的関数領域１０１のペナルティ重み
係数記憶領域１０１２に格納される。「設定終了」５０
の領域が選択されると、データ設定部１３は、目的関数
と考慮する制約条件の選択を終了し、処理選択メニュー
画面２０に戻る。In order to input a value different from the displayed weighting coefficient 42, item consideration 43, and weighting coefficient 45 and constraint consideration 46, the input areas 42, 4 are input via the input device 6.
Enter numerical values in 3, 45, 46. Input areas 42, 4
The values of 3, 45, and 46 are stored in the penalty weight coefficient storage area 1012 of the objective function area 101 by selecting the area of "data storage" 49. “End setting” 50
Is selected, the data setting unit 13 ends the selection of the constraint condition to be considered as the objective function, and returns to the process selection menu screen 20.
【００８４】ｃ．需要データの入力，設定 処理選択メニュー画面２０において、需要データの入
力，設定（項目ｃ）が入力装置６を介して選択される
と、データ設定部１３は処理対象のデータベース８の計
画対象期間記憶領域１００の計画対象年月日と、制約条
件記憶領域１０２の需要記憶領域１０２１から需要の値
を読み込む。需要記憶領域１０２１にデータが格納され
ているときは、格納されている値を、格納されていない
ときはデフォルト値である０が入力された設定画面５１
が表示される。ここで、表示される需要データの画面５
１を図９に示す。この画面５１は、各時刻の需要値の入
力領域５６とその需要値のグラフ表示領域５９と、需要
データを読み込む先である「サーバからのデータ入力」
選択ボタン５３と「ディスクからのデータ入力」選択ボ
タン５４を表示し、各時刻の需要値をスクロールさせる
ためのスクロール５７を備え、さらに「データ反映」５
８と「データ格納」６０と「設定終了」６１の領域があ
る。「サーバからのデータ入力」選択ボタン５３を選択
すると電力システム１０の対象年月日の１日分の需要デ
ータを読み込んで、需要の入力領域５６に表示する。ま
た、「ディスクからのデータ入力」選択ボタン５４を選
択するとデータベース８の需要記憶領域１０２１から需
要データを読み込んで、需要の入力領域５６に表示す
る。また、表示されている需要データを修正するために
は入力装置６を介して需要データの入力領域５６に需要
データを入力する。スクロール５７を選択することによ
り、その選択した方向に応じて、需要データを表示する
時刻が右，左へスクロールして、全体の需要値を確認す
ることができる。「データ反映」ボタン５８を選択する
と、需要データ入力領域５７に表示されている値が需要
値のグラフ表示領域５９に表示される。入力領域５６の
需要値は「データ格納」ボタン６０を選択することによ
り、需要記憶領域１０２１に格納される。「設定終了」
６１の領域が選択されると、データ設定部１３は、目的
関数と考慮する制約条件の選択を終了し、処理選択メニ
ュー画面２０に戻る。C. Input and setting of demand data When input and setting of demand data (item c) are selected via the input device 6 on the process selection menu screen 20, the data setting unit 13 stores the plan target period in the processing target database 8. The demand value is read from the target date in the area 100 and the demand storage area 1021 in the constraint condition storage area 102. When the data is stored in the demand storage area 1021, the stored value is set, and when the data is not stored, the default value 0 is input to the setting screen 51.
Is displayed. Here, the screen 5 of the demand data to be displayed
1 is shown in FIG. The screen 51 includes a demand value input area 56 at each time, a demand value graph display area 59, and a "data input from server" from which demand data is read.
A selection button 53 and a "data input from disk" selection button 54 are displayed, and a scroll 57 for scrolling the demand value at each time is provided.
8, “data storage” 60, and “setting end” 61. When the “data input from server” selection button 53 is selected, the demand data for the target date of the power system 10 for one day is read and displayed in the demand input area 56. When the “input data from disk” selection button 54 is selected, demand data is read from the demand storage area 1021 of the database 8 and displayed in the demand input area 56. In order to correct the displayed demand data, the demand data is input to the demand data input area 56 via the input device 6. By selecting the scroll 57, the time at which the demand data is displayed is scrolled right and left in accordance with the selected direction, so that the entire demand value can be confirmed. When the “data reflection” button 58 is selected, the value displayed in the demand data input area 57 is displayed in the demand value graph display area 59. The demand value in the input area 56 is stored in the demand storage area 1021 by selecting the “data storage” button 60. "End of setting"
When the region 61 is selected, the data setting unit 13 ends the selection of the constraint condition to be considered as the objective function, and returns to the process selection menu screen 20.
【００８５】ｄ．発電機等起動停止計画の入力，設定 処理選択メニュー画面２０において、発電機等起動停止
計画の入力，設定（項目ｄ）が入力装置６を介して選択
されると、データ設定部１３は処理対象のデータベース
８の計画対象期間記憶領域１００の計画対象年月日と、
制約条件記憶領域１０２の発電機等起動停止状態記憶領
域１０２３から発電機等の起動停止状態を表す値を読み
込む。発電機等起動停止状態記憶領域１０２３にデータ
が格納されているときは、格納されている値を、格納さ
れていないときはデフォルト値である１が入力された設
定画面６２が表示される。ここで、表示される需要デー
タの画面６２を図１０に示す。この画面６２は、各時刻
の発電機等の起動停止状態の入力領域６６と、起動停止
状態を読み込む先である「サーバからのデータ入力」選択
ボタン６４と「ディスクからのデータ入力」選択ボタン
６５を表示し、発電機名称方向の起動停止状態をスクロ
ールさせるためのスクロールボタン６７と時刻方向の起
動停止状態をスクロールさせるためのスクロールボタン
６８を備え、さらに「データ格納」選択ボタン６９と「設
定終了」選択ボタン７０の領域がある。「サーバからの
データ入力」選択ボタン６４を選択すると電力システム
１０の対象年月日の１日分の起動停止計画を読み込ん
で、起動停止状態の入力領域６６に表示する。また、
「ディスクからのデータ入力」選択ボタン６５を選択す
るとデータベース８の発電機等起動停止計画記憶領域１
０２４から起動停止状態を読み込んで、起動停止状態の
入力領域６６に表示する。また、表示されている起動停
止状態を修正するためには入力装置６を介して起動停止
状態の入力領域６６に起動停止状態を表す数値を入力す
る。起動停止状態を表す数値は、０と１と２である。０
は停止の状態を、１は運転可能状態を、２は必ず運転す
ることを表す。起動停止状態の０と２は出力計算部１５
により発電機の状態を変更できないことを示し、起動停
止状態の１は制約条件を満たす解がないときに出力計算
部１５による０に変更できる状態であることを示す。
「サーバからのデータ入力」選択ボタン６４を選択して
読み込んだ起動停止状態は０または２である。スクロー
ル６７を選択することにより、その選択した方向に応じ
て、発電機名称を上，下へスクロールさせて、表示して
いる時間帯の全体の起動停止状態を確認することができ
る。同様に、スクロール６８を選択することにより、そ
の選択した方向に応じて、時間軸方向を右，左へスクロ
ールさせて、表示している発電機の全体の起動停止状態
を確認することができる。「データ反映」ボタン５８を
選択すると、需要データ入力領域５７に表示されている
値が需要値のグラフ表示領域５９に表示される。入力領
域５６の需要値は「データ格納」ボタン６０を選択する
ことにより、需要記憶領域１０２１に格納される。「設
定終了」６１の領域が選択されると、データ設定部１３
は、目的関数と考慮する制約条件の選択を終了し、処理
選択メニュー画面２０に戻る。D. Input / setting of generator / startup / stop plan When input / setting of generator / startup / stop plan (item d) is selected via the input device 6 on the process selection menu screen 20, the data setting unit 13 sets the processing target Target date in the target period storage area 100 of the database 8
A value indicating the start / stop state of the generator or the like is read from the generator / startup / stop state storage area 1023 of the constraint condition storage area 102. When data is stored in the generator or other start / stop state storage area 1023, the stored screen is displayed, and when not stored, a setting screen 62 in which a default value of 1 is input is displayed. Here, the displayed demand data screen 62 is shown in FIG. This screen 62 includes an input area 66 for the start / stop state of the generator or the like at each time, a “data input from server” selection button 64 and a “data input from disk” selection button 65 to which the start / stop state is read. And a scroll button 67 for scrolling the start / stop state in the generator name direction and a scroll button 68 for scrolling the start / stop state in the time direction. "Selection button 70 area. When the “input data from server” selection button 64 is selected, the start / stop plan for one day on the target date of the power system 10 is read and displayed in the start / stop state input area 66. Also,
When the “data input from disk” selection button 65 is selected, the generator and other start / stop plan storage area 1
The start / stop state is read from 024 and displayed in the start / stop state input area 66. In order to correct the displayed start / stop state, a numerical value representing the start / stop state is input to the start / stop state input area 66 via the input device 6. Numerical values representing the start / stop state are 0, 1, and 2. 0
Indicates a stopped state, 1 indicates an operable state, and 2 indicates that the apparatus is always operated. 0 and 2 in the start / stop state are output calculation unit 15
Indicates that the state of the generator cannot be changed, and 1 in the start / stop state indicates that the output calculation unit 15 can change the state to 0 when there is no solution satisfying the constraint condition.
The start / stop state read by selecting the “input data from server” selection button 64 is 0 or 2. By selecting the scroll 67, the generator name can be scrolled up or down in accordance with the selected direction, and the start / stop state of the entire displayed time zone can be confirmed. Similarly, by selecting the scroll 68, the time axis direction can be scrolled right and left in accordance with the selected direction, and the start / stop state of the displayed generator as a whole can be confirmed. When the “data reflection” button 58 is selected, the value displayed in the demand data input area 57 is displayed in the demand value graph display area 59. The demand value in the input area 56 is stored in the demand storage area 1021 by selecting the “data storage” button 60. When the area of “setting end” 61 is selected, the data setting unit 13
Terminates the selection of the constraint condition to be considered as the objective function, and returns to the processing selection menu screen 20.
【００８６】ｅ．初期出力と供給力チェックの入力，設
定 処理選択メニュー画面２０において、初期出力と供給力
チェックの入力，設定（項目ｅ）が入力装置６を介して
選択されると、データ設定部１３は処理対象のデータベ
ース８の計画対象期間記憶領域１００の計画対象年月日
と、制約条件記憶領域１０２の需要記憶領域１０２１の
需要と発電機等起動停止状態記憶領域１０２３から発電
機等の起動停止状態を表す値と、初期供給力記憶領域１
０２２から発電機等の初期出力と、初期送電ロス率記憶
領域１０２３から送電ロス率と、発電機運用特性記憶領
域１０２８から発電か負荷（揚水ポンプと対象系統から
外部へでる供給力）かを示すフラグである供給力特性と
発電機の種類を番号で示す設備特性と出力上下限とを読
み込む。需要記憶領域１０２１，発電機等起動停止状態
記憶領域１０２３，初期供給力記憶領域１０２２，初期
送電ロス率記憶領域１０２３，発電機運用特性記憶領域
１０２８にデータが格納されているときは、格納されて
いる値を、格納されていないときはデフォルト値である
１が入力された設定画面７１が表示される。ここで、表
示される画面７１を図１１に示す。この画面７１は、計
画対象期間を示す表示領域７２と、供給力初期値を読み
込む元である「サーバからのデータ入力」選択ボタン７
３と「ディスクからのデータ入力」選択ボタン７４を表
示し、供給力初期値設定領域２０２と、サーバまたはデ
ィスクから読み込んだデータの表示をさせるための「表
示」ボタン７６と表示する時刻を入力する「表示時刻」
入力領域７６とコピー元の時刻入力領域７７とコピー先
の開始時刻入力領域７８とコピー先の終了時刻入力領域
７９とコピーすることを指示する「コピー開始」ボタン
８０と供給力特性を入力および表示する領域８１と運転
状態を表示する領域８２と設備特性を入力および表示す
る領域８３と出力下限を入力および表示する領域８４と
出力上限を入力および表示する領域８５と制約条件を満
たす解がないときに使用する初期値を入力および表示す
る領域８４と発電機名称方向の初期値等をスクロールさ
せるためのスクロールボタン８７と計画対象時刻の供給
力のチェックを指示する「供給力チェック」ボタン８８
と供給力チェック結果を表示する領域８９と時刻方向の
供給力チェック結果をスクロールさせるためのスクロー
ルボタン９９と供給力のチェックによりエラーとなる
「エラー断面数」を表示する領域９７と最大供給力の計
算結果を表示する領域９０と最小供給力の計算結果を表
示する領域９１と需要値を表示する領域９２とロスの計
算結果を表示する領域９３と需要とロスを加えた必要供
給力を表示する領域９４と最大供給力から必要供給力を
引いた上げ余裕量を表示する領域９５と必要供給力から
最小供給力を引いた下げ余裕量を表示する領域９６と初
期ロス率を入力および表示する領域９８を備え、さらに
「データ格納」選択ボタン２００と「設定終了」選択ボ
タン２０１の領域がある。「サーバからのデータ入力」
選択ボタン７３を選択すると電力システム１０の対象年
月日の１日分の供給力特性，運転状態，設備特性，出力
上限，出力下限，初期値および初期ロス率を読み込む。
また、「ディスクからのデータ入力」選択ボタン７４を
選択すると需要記憶領域１０２１の需要と発電機等起動
停止状態記憶領域１０２３から発電機等の起動停止状態
を表す値と、初期供給力記憶領域１０２２から発電機等
の初期出力と、初期送電ロス率記憶領域１０２３から送
電ロス率と、発電機運用特性記憶領域１０２８から供給
力特性と設備特性と出力上下限とを読み込む。これらの
読み込んだ値は、「表示時刻」入力領域７５に時刻を入
力して、「表示」ボタン７６を選択することにより、指
定した時刻の数値が領域２０２に表示される。表示され
ている数値を修正するためには入力装置６を介して入力
領域の特性を持つ「供給力特性」領域８１，「設備特
性」領域８３，「出力下限」領域８４，「出力上限」領
域８５，「初期値」領域８６および「初期ロス率」領域
９８に数値を入力する。また、「コピー元の時刻」入力
領域７７に時刻を入力し、「コピー先開始時刻」入力領
域７８と「コピー先終了時刻」入力領域７９に時刻を入
力して、「コピー開始」ボタン８０を選択すると、コピ
ー元の時刻の入力領域の特性を持つ「供給力特性」領域
８１，「設備特性」領域８３，「出力下限」領域８４，
「出力上限」領域８５，「初期値」領域８６および「初
期ロス率」領域９８のデータがコピー先の時刻のデータ
領域にコピーされる。「供給特性」８１の値が１のとき
対象系統に対して供給力が発電であることを表し、「供
給特性」８１の値が−１のとき対象系統に対して供給力
が負荷となることを表す。「供給特性」８１の値が−１
となる供給力には、揚水ポンプ、連系線からでていく融
通，託送による供給力である。「設備特性」８３は供給
力の種類を表し、火力発電機を１，揚水発電を２，揚水
ポンプを３，水力を４，ＩＰＰを５，融通受電を６，融
通送電を７，原子力を８としている。火力を燃料によ
り、重油，ＬＮＧ，石炭等に分類して番号を付けても良
い。E. Input and setting of initial output and supply power check When input and setting (item e) of initial output and supply power check are selected via the input device 6 on the process selection menu screen 20, the data setting unit 13 Of the planning target period storage area 100 of the database 8, the demand of the demand storage area 1021 of the constraint condition storage area 102, and the start / stop state of the generator / etc. From the generator / start / stop state storage area 1023. Value and initial supply capacity storage area 1
022 indicates the initial output of the generator or the like, the initial transmission loss rate storage area 1023 indicates the power transmission loss rate, and the generator operation characteristic storage area 1028 indicates the power generation or the load (supplying power from the pump and the target system to the outside). The power supply characteristics and the equipment characteristics indicating the type of the generator, which are flags, and the output upper and lower limits are read. When data is stored in the demand storage area 1021, the generator and the like start / stop state storage area 1023, the initial supply power storage area 1022, the initial power transmission loss rate storage area 1023, and the generator operation characteristic storage area 1028, the data is stored. When the stored value is not stored, the setting screen 71 in which the default value of 1 is input is displayed. Here, the displayed screen 71 is shown in FIG. This screen 71 includes a display area 72 indicating a planning target period, and a “input data from server” selection button 7 from which a supply initial value is read.
3 and a "input data from disk" selection button 74 are displayed, and a supply power initial value setting area 202, a "display" button 76 for displaying data read from the server or the disk, and a display time are input. "Display time"
An input area 76, a copy source time input area 77, a copy destination start time input area 78, and a copy destination end time input area 79, a "copy start" button 80 for instructing to copy, and a supply power characteristic are input and displayed. When there is no solution that satisfies the constraint conditions, a region 81 for displaying the operating state, a region 83 for inputting and displaying the facility characteristics, a region 84 for inputting and displaying the lower output limit, a region 85 for inputting and displaying the upper output limit, and a region 85 for inputting and displaying the upper output limit. 84 for inputting and displaying an initial value used for the operation, a scroll button 87 for scrolling the initial value in the generator name direction and the like, and a "supply power check" button 88 for instructing a check on the supply power at the planned time.
89, a scroll button 99 for scrolling the supply power check result in the time direction, an area 97 for displaying the “number of error cross-sections” resulting in an error due to the supply power check, and a maximum supply power An area 90 for displaying the calculation result, an area 91 for displaying the calculation result of the minimum supply capacity, an area 92 for displaying the demand value, an area 93 for displaying the calculation result of the loss, and a required supply capacity including the demand and the loss are displayed. An area 94, an area 95 for displaying a surplus amount obtained by subtracting the necessary supply power from the maximum supply power, an area 96 for displaying a reduction allowance amount obtained by subtracting the minimum supply power from the necessary supply power, and an area for inputting and displaying the initial loss rate. 98, and an area for a “data storage” selection button 200 and a “setting end” selection button 201. "Data input from server"
When the selection button 73 is selected, the supply power characteristics, the operating state, the equipment characteristics, the output upper limit, the output lower limit, the initial value, and the initial loss rate of the power system 10 for one day are read.
When the “input data from disk” selection button 74 is selected, the demand in the demand storage area 1021 and the value indicating the start / stop state of the generator and the like from the start / stop state storage area 1023 of the generator and the like, and the initial supply capacity storage area 1022 From the initial power transmission loss rate storage area 1023, and the power supply characteristics, facility characteristics, and output upper and lower limits from the generator operation characteristic storage area 1028. By inputting a time in the “display time” input area 75 and selecting a “display” button 76, a numerical value of the designated time is displayed in the area 202. In order to correct the displayed numerical value, the "supply force characteristic" region 81, the "equipment characteristic" region 83, the "output lower limit" region 84, the "output upper limit" region having the characteristics of the input region via the input device 6. Numerical values are input to an “initial value” area 86 and an “initial loss rate” area 98. The user inputs a time in the “copy source time” input area 77, inputs a time in the “copy destination start time” input area 78 and a “copy destination end time” input area 79, and presses the “copy start” button 80. When selected, the “supply power characteristic” area 81, the “equipment characteristic” area 83, the “output lower limit” area 84,
The data in the “output upper limit” area 85, the “initial value” area 86, and the “initial loss rate” area 98 are copied to the data area at the copy destination time. When the value of the “supply characteristic” 81 is 1, the supply power is power generation for the target system, and when the value of the “supply characteristic” 81 is −1, the supply power is a load for the target system. Represents The value of the “supply characteristic” 81 is −1
The supply power is the supply power from pumping pumps, interchanges and consignments that go out of the interconnection line. The “equipment characteristics” 83 indicates the type of supply power, and indicates a thermal power generator 1, a pumped-storage power generator 2, a pumped pump 3, a hydraulic power supply 4, an IPP 5, a flexible power reception 6, a flexible power transmission 7, and a nuclear power supply 8 And The thermal power may be classified by fuel into heavy oil, LNG, coal, etc., and may be numbered.
【００８７】「供給力チェック」ボタン８８を選択する
と、「表示時刻」７６の需要を表示領域９２に、初期ロ
ス率を表示し、需要と初期ロス率および領域２０２に表
示されているデータから領域２０３の値を計算して、表
示する。最大供給力９０は運転状態が１または２かつ供
給特性は１の発電機等の出力上限値の和である。最小供
給力９１は運転状態が１または２でさらに設備特性が
１，５，６または８の出力下限と運転状態が１または２
でさらに設備特性が３または７の出力上限と運転状態が
２でさらに設備特性が２または４の出力下限との和であ
る。「送電ロス」９３は需要９２と初期ロス率９８の積
である。必要供給力９５は需要９２と送電ロス９３の和
である。上げ余裕量９５は最大供給力９０と必要供給力
９４の差である。下げ余裕量９６は必要供給力９４と最
小供給力９０の差である。需給バランス制約を満たす解
が存在するかどうか事前にチェックするための機能であ
り、上げ余裕量９５および下げ余裕量９６が正の値で有
れば良い。供給力チェック判定領域８９は各時刻の上げ
余裕量および下げ余裕量の両方ともが正の値のとき、フ
ラグとして１を表示し、そうでないときは０を表示す
る。エラー断面数９７は供給力チェック判定領域８９に
ある０の数を表示する。需給バランスを制約条件として
考慮し、かつエラー断面が存在するときは、出力上限，
出力下限または初期ロス率または運転状態８２を変更す
る必要がある。但し、運転状態８２の変更は発電機等起
動停止計画の入力，設定（項目ｄ）にて変更する必要が
ある。When the “supply capacity check” button 88 is selected, the demand at the “display time” 76 is displayed in the display area 92 and the initial loss rate is displayed. The value of 203 is calculated and displayed. The maximum supply force 90 is the sum of the output upper limit values of a generator or the like whose operation state is 1 or 2 and whose supply characteristic is 1. The minimum supply power 91 is 1 or 2 when the operating state is 1 or 2, and the output lower limit and the operating state are 1 or 2 when the equipment characteristic is 1, 5, 6 or 8.
Is the sum of the output upper limit when the equipment characteristic is 3 or 7 and the output lower limit when the operation state is 2 and the equipment characteristic is 2 or 4. The “power transmission loss” 93 is the product of the demand 92 and the initial loss rate 98. The necessary supply power 95 is the sum of the demand 92 and the power transmission loss 93. The raising margin 95 is the difference between the maximum supply force 90 and the required supply force 94. The lowering allowance 96 is the difference between the required supply force 94 and the minimum supply force 90. This is a function for checking in advance whether there is a solution that satisfies the supply-demand balance constraint, and it is sufficient that the allowance 95 and the allowance 96 have positive values. The supply capacity check determination area 89 displays 1 as a flag when both the allowance for increase and the allowance for decrease at each time are positive values, and displays 0 otherwise. The number of error cross sections 97 indicates the number of 0 in the supply capacity check determination area 89. When the supply and demand balance is considered as a constraint and an error cross section exists,
It is necessary to change the lower output limit or the initial loss rate or the operating state 82. However, it is necessary to change the operation state 82 by inputting and setting (item d) the start / stop plan of the generator or the like.
【００８８】スクロール９７を選択することにより、そ
の選択した方向に応じて、発電機名称を上，下へスクロ
ールさせて、表示対象の発電機を変更して供給力の初期
値等を確認することができる。同様に、スクロール９９
を選択することにより、その選択した方向に応じて、時
間軸方向を右，左へスクロールさせて、表示している供
給力チェックフラグを確認することができる。表示，入
力領域８１，８３，８４，８５，８６，９８の値は、
「データ格納」ボタン２００を選択することにより、発
電機等運用特性記憶領域１０２１に格納される。「設定
終了」２０１の領域が選択されると、データ設定部１３
は、初期出力と供給力チェックの入力，設定を終了し、
処理選択メニュー画面２０に戻る。By selecting the scroll 97, the generator name is scrolled up or down in accordance with the selected direction, and the generator to be displayed is changed to confirm the initial value of the supply power and the like. Can be. Similarly, scroll 99
Is selected, the time axis direction can be scrolled right and left in accordance with the selected direction, and the displayed supply force check flag can be confirmed. The values of the display and input areas 81, 83, 84, 85, 86, 98 are
By selecting the “data storage” button 200, the data is stored in the operation characteristic storage area 1021 such as a generator. When the area of “setting end” 201 is selected, the data setting unit 13
Finishes input and setting of initial output and supply power check,
It returns to the process selection menu screen 20.
【００８９】ｆ．発電機等設備特性の入力，設定 処理選択メニュー画面２０において、発電機等設備特性
の入力，設定(項目ｆ)が入力装置６を介して選択される
と、データ設定部１３は処理対象のデータベース８の計
画対象期間記憶領域１００の計画対象年月日と、目的関
数記憶領域101の発電機等設備特性記憶領域１０１１の
発電機等のコスト特性と揚水効率（１ＭＷｈで汲み上げ
た水により発電できるＭＷｈ）と電水比（１ton ／秒で
ダムから水を放流したときの発電量）を表す値と、供給
力積分制約記憶領域１０２５から揚水，揚発および水系
の各時刻の水位下限と上限（ＭＷｈ換算値）とを読み込
む。表示時刻はデフォルトの１が最初に入力され、この
時刻の目的関数記憶領域１０１１，供給力積分制約記憶
領域１０２５にデータが格納されているときは、格納さ
れている値を、格納されていないときはデフォルト値で
ある１が入力された設定画面２２０が表示される。ここ
で、表示される画面２２０を図１２に示す。この画面２
２０は、計画対象期間を示す表示領域２２１と、発電機
等設備特性と揚水，水系の水位上下限制約を読み込む元
である「サーバからのデータ入力」選択ボタン２２２と
「ディスクからのデータ入力」選択ボタン２２３を表示
し、コスト特性設定領域２２６と、サーバまたはディス
クから読み込んだデータの表示をさせるための「表示」
ボタン２２５と表示する時刻を入力する「表示時刻」入
力領域２２４とコピー元の時刻入力領域２３７とコピー
先の開始時刻入力領域２３８とコピー先の終了時刻入力
領域２３９とコピーすることを指示する「コピー開始」
ボタン２４０とコスト特性の２次係数を入力および表示
する領域２２７とコスト特性の１次係数を入力および表
示する領域２２８とコスト特性の定数を入力および表示
する領域２２９と揚水，水系の水位上下限を表示，入力
する領域２３０と水位下限の入力，表示領域２３１と水
位上限の入力，表示領域２３２と効率電水比の入力，表
示領域２３３と渓流から水系のダムに流れ込む流入量の
入力，表示２３４と発電機名称方向のコスト特性等をス
クロールさせるためのスクロールボタン２３６と水位上
下限値を表示する揚水，水系の名称をスクロールさせる
ためのスクロールボタン２３６とを備え、さらに「デー
タ格納」選択ボタン２４１と「設定終了」選択ボタン２
４２の領域がある。「サーバからのデータ入力」選択ボ
タン２２２を選択すると電力システム１０の対象年月日
の１日分のコスト特性，揚水，水系の水位上下限，効
率，電水比および流入量を読み込む。また、「ディスク
からのデータ入力」選択ボタン２２３を選択すると発電
機等設備特性記憶領域１０１１の発電機等のコスト特性
と揚水効率と電水比を表す値と、供給力積分制約記憶領
域１０２５から揚水，揚発および水系の各時刻の水位下
限と上限を表す値を読み込む。これらの読み込んだ値
は、「表示時刻」入力領域２２４に時刻を入力して、
「表示」ボタン２２５を選択することにより、指定した
時刻の数値が領域２２６に表示される。表示されている
数値を修正するためには入力装置６を介して入力領域の
特性を持つ「コスト特性の２次係数」領域２２７，「コス
ト特性の１次係数」領域２２８，「コスト特性の定数」
領域２２９，「水位下限」領域２３１，「水位上限」領
域２３２，「効率電水比」領域２３３，「流入量」領域
２３４に数値を入力する。また、「コピー元の時刻」入
力領域２３７に時刻を入力し、「コピー先開始時刻」入
力領域２３８と「コピー先終了時刻」入力領域２３９に
時刻を入力して、「コピー開始」ボタン２４０を選択す
ると、コピー元の時刻の入力領域の特性を持つ「コスト
特性の２次係数」領域２２７，「コスト特性の１次係
数」領域２２８，「コスト特性の定数」領域２２９，
「水位下限」領域２３１，「水位上限」領域２３２，「効
率電水比」領域２３３，「流入量」領域２３４のデータ
がコピー先の時刻のデータ領域にコピーされる。例え
ば、火力１のコスト特性は、出力の二乗の２次の係数
０.０２ を掛けた値と出力に１次の係数２.０１ を掛け
た値と定数１００の和である。揚水，揚発はコスト係数
が全て０となっているが、揚水ポンプにより水を汲み上
げるには、火力等による出力増によるコストアップが生
じるため、揚水をすると結果的にコストが増加する。こ
れらのコスト係数に０以外の値を入力しても良い、例え
ば揚水，揚発により設備の寿命が減少する場合に値を入
力する等である。原子力についてもコスト係数が０とな
っているのは、現在は固定供給力として考えられている
ためであり、出力変動が許容されるならば、コスト係数
を設定して需要に応じて出力を変化させることもでき
る。ＩＰＰのコスト係数は１次係数以外は０であるが、
これも２次，定数に０以外の値を入力してもよい。コス
ト係数の１次，定数の値は負でもよいが、２次の係数は
０または正とする。F. Input and setting of generator and other equipment characteristics When the input and setting of generator and other equipment characteristics (item f) is selected via the input device 6 on the process selection menu screen 20, the data setting unit 13 sets the database to be processed. 8, the planning target date in the planning target period storage area 100, the cost characteristics of the generator in the facility characteristic storage area 1011 in the objective function storage area 101, and the pumping efficiency (MWh that can be generated by water pumped at 1 MWh). ) And the electric water ratio (the amount of power generated when water is discharged from the dam at 1 ton / sec), and the water level lower limit and upper limit (MWh) at each time of pumping, pumping, and water systems from the supply power integration constraint storage area 1025. And the converted value). As the display time, a default of 1 is input first, and when the data is stored in the objective function storage area 1011 and the supply force integration restriction storage area 1025 at this time, the stored value is not stored. Is displayed on the setting screen 220 in which 1 as a default value is input. Here, the screen 220 displayed is shown in FIG. This screen 2
Reference numeral 20 denotes a display area 221 indicating a planning target period, a "data input from server" selection button 222, and a "data input from disk" from which the characteristics of generators and other facilities and the water level upper and lower limits of pumping and water systems are read. A selection button 223 is displayed, and a “display” for displaying the cost characteristic setting area 226 and data read from the server or the disk.
A button 225 and a “display time” input area 224 for inputting a display time, a copy source time input area 237, a copy destination start time input area 238, and a copy destination end time input area 239 are instructed to copy. Start Copy "
A button 240, an area 227 for inputting and displaying the second order coefficient of the cost characteristic, an area 228 for inputting and displaying the first order coefficient of the cost characteristic, an area 229 for inputting and displaying the constant of the cost characteristic, and upper and lower limits of the water level of the pumping and water system. Display and input area 230 and input of water level lower limit, display area 231 and input of water level upper limit, display area 232 and input of efficient water-to-water ratio, display area 233 and input and inflow of inflow from mountain stream to water system dam 234, a scroll button 236 for scrolling the cost characteristics and the like in the direction of the generator name, a scroll button 236 for scrolling the pumping and water system names for displaying the upper and lower limits of the water level, and a "data storage" selection button. 241 and "Setting end" selection button 2
There are 42 areas. When the "input data from server" selection button 222 is selected, the cost characteristics, pumping, water level upper and lower limits, efficiency, power-to-water ratio, and inflow of the power system 10 for one day are read. When the “input data from disk” selection button 223 is selected, the cost characteristics of the generator and the like, the pumping efficiency and the electric water ratio in the generator and the like characteristic storage area 1011, and the supply power integration constraint storage area 1025 are displayed. The values representing the lower and upper water levels at each time of pumping, pumping, and water systems are read. These read values are input into the “display time” input area 224 by inputting the time,
By selecting the “display” button 225, the numerical value of the designated time is displayed in the area 226. In order to correct the displayed numerical value, the “second-order coefficient of cost characteristic” region 227, the “first-order coefficient of cost characteristic” region 228, and the “constant of cost characteristic” having the characteristics of the input region via the input device 6 "
Numerical values are input to the area 229, the “lower water level” area 231, the “higher water level” area 232, the “efficiency electric water ratio” area 233, and the “inflow amount” area 234. Further, the user inputs a time in the “copy source time” input area 237, inputs a time in the “copy destination start time” input area 238 and a “copy destination end time” input area 239, and presses the “copy start” button 240. When selected, a “secondary coefficient of cost characteristic” region 227, a “first coefficient of cost characteristic” region 228, a “constant of cost characteristic” region 229 having the characteristics of the input region of the copy source time,
The data in the “lower water level” area 231, the “higher water level” area 232, the “efficient water-to-water ratio” area 233, and the “inflow amount” area 234 are copied to the data area at the copy destination time. For example, the cost characteristic of the thermal power 1 is a sum of a value obtained by multiplying the output by a quadratic coefficient of 0.02, a value obtained by multiplying the output by a primary coefficient of 2.01, and a constant of 100. Pumping and pumping all have a cost coefficient of 0, but pumping up water with a pump raises costs due to an increase in output due to thermal power or the like. Values other than 0 may be input to these cost coefficients, for example, when the life of the equipment is reduced due to pumping and pumping. The reason why the cost coefficient for nuclear power is 0 is that it is currently considered as a fixed supply capacity. If output fluctuations are allowed, the cost coefficient is set and the output changes according to demand. It can also be done. The IPP cost coefficient is 0 except for the primary coefficient,
Again, a value other than 0 may be input to the secondary and constants. The primary and constant values of the cost coefficient may be negative, but the secondary coefficient is 0 or positive.
【００９０】スクロール２３５，２３６を選択すること
により、その選択した方向に応じて、発電機名称を上，
下へスクロールさせて、表示対象の発電機を変更してコ
スト特性，水位上下限等を確認することができる。「デ
ータ格納」ボタン２４１を選択することにより、表示，
入力領域２２７，２２８，２２９，２３３，２３４の値
は、発電機等設備特性記憶領域１０１１に、表示，入力
領域２３１，２３２の値は、供給力積分制約記憶領域１
０１１にそれぞれ格納される。「設定終了」２４２の領
域が選択されると、データ設定部１３は、発電機等設備
特性の入力，設定を終了し、処理選択メニュー画面２０
に戻る。By selecting the scrolls 235 and 236, the generator name is changed to “up” or “down” according to the selected direction.
By scrolling down, the generator to be displayed can be changed to check cost characteristics, water level upper and lower limits, and the like. By selecting the “store data” button 241, the display,
The values of the input areas 227, 228, 229, 233, and 234 are stored in the equipment characteristic storage area 1011. The values of the display and input areas 231, 232 are stored in the supply force integration constraint storage area 1
011. When the area of “setting end” 242 is selected, the data setting unit 13 ends the input and setting of the equipment characteristics such as the generator, and the processing selection menu screen 20
Return to
【００９１】ｇ．供給力積分特性の入力，設定 処理選択メニュー画面２０において、供給力積分特性の
入力，設定(項目ｇ)が入力装置６を介して選択される
と、データ設定部１３は処理対象のデータベース８の計
画対象期間記憶領域１００の計画対象年月日と、制約条
件記憶領域１０２の発電機運用特性記憶領域１０２８の
設備特性と、供給力積分制約記憶領域1025の積分制約で
同一グループであることを示す制約グループ番号と定格
出力と計画対象期間で発電する日量下限と日量上限と揚
水池，貯水池等の初期水位と日量制約の重み係数とを読
み込む。発電機運用特性記憶領域１０２８，供給力積分
制約記憶領域１０２５にデータが格納されているとき
は、格納されている値を、格納されていないときはデフ
ォルト値である１が入力された設定画面２０２が表示さ
れる。ここで、表示される画面２０２を図１３に示す。
この画面２０２は、計画対象期間を示す表示領域２０６
と、供給力積分制約を読み込む元である「サーバからの
データ入力」選択ボタン２０４と、「ディスクからのデ
ータ入力」選択ボタン２０５を表示し、設備特性表示領
域２０７と、制約グループ表示入力表示領域２０８と、
定格出力入力表示領域２０９と、日量下限入力表示領域
２１０と、日量上限入力表示領域２１１と、利用率下限
表示領域２１２と、利用率上限表示領域２１３と、初期
水位入力表示領域２１４と、重み係数入力表示領域２１
５と、発電機名称方向に制約領域２０６をスクロールさ
せるためのスクロールボタン２１７と日量から利用率の
計算を指示する「利用率計算」ボタン２０６とを備え、
さらに「データ格納」選択ボタン２００と「設定終了」
選択ボタン２０１の領域がある。「サーバからのデータ
入力」選択ボタン２０４を選択すると電力システム１０
の対象年月日の設備特性，制約グループ，定格出力，日
量下限，日量上限，初期水位および重み係数を読み込ん
で、領域２０６に表示する。また、「ディスクからのデ
ータ入力」選択ボタン２０５を選択すると発電機運用特
性記憶領域１０２８の設備特性と、供給力積分制約記憶
領域１０２５の積分制約で同一グループであることを示
す制約グループ番号と定格出力と計画対象期間で発電す
る日量下限と日量上限と揚水池，貯水池等の初期水位と
日量制約の重み係数とを読み込んで、領域２０６に表示
する。表示されている数値を修正するためには入力装置
６を介して入力領域の特性を持つ「制約グループ」領域
２０８，「定格出力」領域２０９，「日量下限」領域２１
０，「日量上限」領域２１１，「初期水位」領域２１４
および「重み係数」領域２１５に数値を入力する。制約
グループは積分制約を考慮するときに、同時に積分（加
算）すべき発電機等であることを示す。たとえば、図１
２の例では、火力３と火力４が同一のグループ番号１と
なっている。このとき、計画対象期間の日量下限値は火
力３と火力４の日量下限を加算したものとなる。すなわ
ち、火力３の日量が６０００であっても、火力４と火力
４の和が１８０００以上あれば、下限制約を満たすこと
になる。同様に、初期水位も、同一グループ番号の初期
水位を加算したものがグループでの初期水位になる。同
一グループ番号が１つしかないものは、単独の制約であ
ることを表している。但し、設備特性が２，３の揚水，
揚発，設備特性が４の水力の日量制約はないので、対象
外である。揚水，揚発，水系の日量に相当するものは、
池水位であり、これらについては池水位の上下限として
扱うためである。重み係数２１５は出力計算部１５で出
力を計算するとき、制約条件を満たす解がないときに、
日量制約に対する違反にペナルティコストを目的関数に
加算するときの重み係数として使用する。この大きさ
は、日量制約の中での制約の重要なものほど大きくする
ことになる。最終的な重み係数は図８の重き係数４５の
日量制約の値との積で決まる。「利用率計算」ボタン２
０６を選択すると、定格出力２０９と計画対象期間（こ
の例では２４時間）の積で「日量下限」２１０を割った
値を利用率下限として表示する。同様に定格出力２０９
と計画対象期間（この例では２４時間）の積で「日量上
限」２１１を割った値を利用率上限として表示する。G. Input and Setting of Supply Force Integration Characteristics When the input and setting (item g) of the supply force integration characteristics are selected via the input device 6 on the process selection menu screen 20, the data setting section 13 stores the data of the database 8 to be processed. It shows that they belong to the same group based on the planning target date in the planning target period storage area 100, the facility characteristics in the generator operation characteristic storage area 1028 in the constraint condition storage area 102, and the integration constraint in the supply power integration restriction storage area 1025. The constraint group number, the rated output, the lower limit of daily power and the upper limit of daily power to be generated during the planning period, the initial water level of the pumping basin, the reservoir, etc., and the weighting factor of the daily constraint are read. When the data is stored in the generator operation characteristic storage area 1028 and the supply power integration restriction storage area 1025, the stored value is set, and when the data is not stored, the default value 1 is input to the setting screen 202. Is displayed. Here, a screen 202 to be displayed is shown in FIG.
This screen 202 includes a display area 206 indicating a plan target period.
And a "data input from server" selection button 204 and a "data input from disk" selection button 205 from which the supply force integration constraint is read, are displayed, and a facility characteristic display area 207 and a constraint group display input display area are displayed. 208,
A rated output input display area 209, a daily amount lower limit input display area 210, a daily amount upper limit input display area 211, a utilization rate lower limit display area 212, a utilization rate upper limit display area 213, an initial water level input display area 214, Weight coefficient input display area 21
5, a scroll button 217 for scrolling the restricted area 206 in the generator name direction, and a "utilization rate calculation" button 206 for instructing the utilization rate to be calculated from the daily amount.
Furthermore, a "data storage" selection button 200 and a "setting end"
There is an area for the selection button 201. When the “input data from server” selection button 204 is selected, the power system 10
, The facility characteristics, the constraint group, the rated output, the daily amount lower limit, the daily amount upper limit, the initial water level, and the weighting factor are read and displayed in the area 206. When the “input data from disk” selection button 205 is selected, the equipment characteristics in the generator operation characteristic storage area 1028 and the restriction group number and rating indicating that they belong to the same group by the integration restriction in the supply power integration restriction storage area 1025. The output, the daily amount lower limit, the daily amount upper limit, the initial water level of the pumping basin, the reservoir, etc., and the weighting factor of the daily amount constraint in the planned period are read and displayed in the area 206. In order to correct the displayed numerical value, the “restriction group” area 208, the “rated output” area 209, and the “daily lower limit” area 21 having the characteristics of the input area via the input device 6.
0, “daily amount upper limit” area 211, “initial water level” area 214
Then, a numerical value is input to the “weight coefficient” area 215. The constraint group indicates a generator or the like to be integrated (added) at the same time when considering the integration constraint. For example, FIG.
In the example of 2, the thermal power 3 and the thermal power 4 have the same group number 1. At this time, the daily amount lower limit of the planning target period is the sum of the daily amount lower limits of the thermal power 3 and the thermal power 4. That is, even if the daily amount of the thermal power 3 is 6000, if the sum of the thermal power 4 and the thermal power 4 is 18000 or more, the lower limit constraint is satisfied. Similarly, the initial water level of the group is obtained by adding the initial water levels of the same group number. If there is only one same group number, it indicates that it is a single constraint. However, pumping with equipment characteristics of a few
It is out of scope because there is no restriction on the daily amount of the hydropower with the emergence and equipment characteristics of 4. The daily equivalent of pumping, pumping, and water systems is
Pond water level, because these are treated as upper and lower limits of pond water level. The weight coefficient 215 is calculated by the output calculation unit 15 when there is no solution satisfying the constraint condition.
It is used as a weighting factor when the penalty cost is added to the objective function for the violation of the daily amount constraint. This magnitude will be greater for the more important of the daily constraints. The final weighting factor is determined by the product of the weighting factor 45 in FIG. 8 and the daily amount constraint value. "Calculation of utilization rate" button 2
When 06 is selected, the value obtained by dividing the “daily lower limit” 210 by the product of the rated output 209 and the planned period (24 hours in this example) is displayed as the utilization rate lower limit. Similarly, rated output 209
And the value obtained by dividing the “daily amount upper limit” 211 by the product of the schedule target period (24 hours in this example) is displayed as the utilization rate upper limit.
【００９２】スクロール２１７を選択することにより、
その選択した方向に応じて、発電機名称を上，下へスク
ロールさせて、表示対象の発電機を変更して日量制約等
を確認することができる。表示，入力領域２０６，２０
９，２１０，２１１，２１４，２１５の値は、「データ
格納」ボタン２１８を選択することにより、供給力積分
制約記憶領域１０２５に格納される。「設定終了」２１
９の領域が選択されると、データ設定部１３は、供給力
積分制約の入力，設定を終了し、処理選択メニュー画面
２０に戻る。By selecting the scroll 217,
In accordance with the selected direction, the generator name can be scrolled up and down to change the generator to be displayed and confirm the daily amount constraint and the like. Display, input area 206, 20
The values of 9, 210, 211, 214, and 215 are stored in the supply force integration constraint storage area 1025 by selecting a “data storage” button 218. "End setting" 21
When the area 9 is selected, the data setting unit 13 ends the input and setting of the supply force integration constraint, and returns to the process selection menu screen 20.
【００９３】ｈ．系統接続状態の入力，設定 処理選択メニュー画面２０において、系統接続状態の入
力，設定(項目ｈ)が入力装置６を介して選択されると、
データ設定部１３は処理対象のデータベース８の計画対
象期間記憶領域１００の計画対象年月日と、制約条件記
憶領域１０２の系統接続状態記憶領域１０２６のブラン
チの運用状態のフラグとブランチの潮流制約の上下限と
発電機の接続先ノード番号と負荷の固定分の負荷率と変
動分の負荷率の値と発電機等起動停止状態記憶領域１０
２４から発電機の運転状態を読み込む。表示時刻はデフ
ォルトの１が入力され、この時刻の系統接続状態記憶領
域１０２６と発電機等起動停止状態記憶領域１０２４に
データが格納されているときは、格納されている値を、
格納されていないときはデフォルト値である１が入力さ
れた設定画面２４３が表示される。ここで、表示される
画面２４３を図１４に示す。この画面２４３は、計画対
象期間を示す表示領域２４４と、系統接続状態を読み込
む元である「サーバからのデータ入力」選択ボタン２４
５と「ディスクからのデータ入力」選択ボタン２４６を
表示し、系統接続図表示領域２４９と、サーバまたはデ
ィスクから読み込んだデータの表示をさせるための「表
示」ボタン２４７と表示する時刻を入力する「表示時
刻」入力領域２４８とコピー元の時刻入力領域２６２と
コピー先の開始時刻入力領域２６３とコピー先の終了時
刻入力領域２６４とコピーすることを指示する「コピー
開始」ボタン２６５と「データ反映」ボタン２６６とブ
ランチの運用状態と潮流制約値の領域２５０とブランチ
の使用／停止の運用状態をフラグで示す入力，表示領域
２５１と潮流制約の下限値を入力，表示する領域２５２
と潮流制約の上限値を入力，表示する領域２５３と発電
機の運転状態と接続先ノードの領域２５５と発電機の運
転状態を示す領域２５６と発電機の接続先ノード番号を
入力，表示する領域２５７とノードの負荷比率の領域２
５９と各ノードの負荷比率の固定分を入力，表示する領
域２６０と各ノードの負荷比率の変動分を入力，表示す
る領域２６１と運用状態，潮流制約をブランチ番号方向
にスクロールさせるためのスクロールボタン２５４と運
転状態，接続先ノード番号を発電機名称方向にスクロー
ルさせるためのスクロールボタン２５８を備え、さらに
「データ格納」選択ボタン２６７と「設定終了」選択ボ
タン２６８の領域がある。「サーバからのデータ入力」
選択ボタン２４５を選択すると電力システム１０の対象
年月日の１日分のブランチの運用状態，ブランチの潮流
制約値，発電機の運転状態とその接続先ノード番号，負
荷の固定分比率と変動分比率を読み込む。また、「ディ
スクからのデータ入力」選択ボタン２４６を選択すると
系統接続状態記憶領域１０２６のブランチの運用状態の
フラグとブランチの潮流制約の上下限と発電機の接続先
ノード番号と負荷の固定分の負荷率と変動分の負荷率の
値と発電機等起動停止状態記憶領域１０２４から発電機
の運転状態を表す値を読み込む。これらの読み込んだ値
は、「表示時刻」入力領域２４７に時刻を入力して、「デ
ータ反映」ボタン２６６を選択することにより、指定し
た時刻のブランチの接続状態が領域２４９に表示され
る。ここで、ブランチとは送変電設備のことであり、送
電線，変圧器等のことである。使用されているブランチ
は実線，停止しているブランチは点線で表示される。黒
いまるはノードを示し、母線等を表す。領域２４９に表
示されている接続状態を修正するためには入力装置６を
介して入力領域の特性を持つ「運用状態」領域２５１の
値を変更し、「データ反映」ボタン２６６を選択する。
また、表示されている数値は入力装置６を介して入力領
域の特性を持つ領域２５１，２５２，２５３，２５７,
２６０,２６１の値を修正することができる。本実施例
では、領域２４９の表示で、潮流の正の方向を左から
右、中心から外側へとしている。このため、潮流制約の
下限２５１は負の値となっている。ブランチ番号は領域
２４９では数字の最初にＬを付けて表している。ノード
番号は領域２４９では数字の最初にＮを付けて表してい
る。負荷とその負荷が接続するノードの番号は同じもの
としている。各ノードの負荷は、需要と送電ロスの和か
ら負荷比率の固定負荷の和を引いた値が変動負荷であ
る。負荷比率の和に対する各ノードの負荷比率にこの変
動負荷を掛けた値が各ノードの変動負荷である。この変
動負荷と固定負荷を加算した値が各時刻，各ノードの負
荷となる。送電ロスは各ブランチで発生するが、ここで
は送電ロス分も負荷に組み込んでいる。H. Input and setting of system connection state When input and setting (item h) of the system connection state are selected via the input device 6 on the process selection menu screen 20,
The data setting unit 13 sets a planning target date in the planning target period storage area 100 of the processing target database 8, a flag of a branch operation state in the system connection state storage area 1026 of the constraint condition storage area 102, and a branch flow restriction. The upper and lower limits, the connection node number of the generator, the load factor for the fixed load and the load factor for the variation, and the generator start / stop state storage area 10
From 24, the operation state of the generator is read. As the display time, a default of 1 is input, and when data is stored in the system connection state storage area 1026 and the generator / startup / stop state storage area 1024 at this time, the stored value is
If it is not stored, a setting screen 243 in which a default value of 1 has been input is displayed. Here, the screen 243 displayed is shown in FIG. This screen 243 includes a display area 244 indicating a planning target period, and a “data input from server” selection button 24 from which the system connection state is read.
5 and a "data input from disk" selection button 246 are displayed, and a system connection diagram display area 249, a "display" button 247 for displaying data read from the server or the disk, and a display time are input. A “copy start” button 265 for instructing to copy the “display time” input area 248, the copy source time input area 262, the copy destination start time input area 263, and the copy destination end time input area 264, and a “data reflection” A button 266, an operation state of the branch, an area 250 of the power flow constraint value, and an input indicating the operation state of use / stop of the branch by a flag, a display area 251 and an area 252 for inputting and displaying the lower limit value of the power flow restriction
253, an area for inputting and displaying the upper limit value of the power flow restriction, an area 255 for displaying the operating state of the generator and the connection destination node, an area 256 for displaying the operating state of the generator, and a node number for the connection of the generator 257 and the area 2 of the node load ratio
59 and a fixed portion of the load ratio of each node, and a display area 260 and a change amount of the load ratio of each node are input and displayed. A display button 261 and a scroll button for scrolling the operation state and the power flow restriction in the branch number direction. 254, a scroll button 258 for scrolling the operation state and the connection destination node number in the direction of the generator name, and an area for a "data storage" selection button 267 and a "setting end" selection button 268. "Data input from server"
When the selection button 245 is selected, the operation state of the branch for one day on the target date of the power system 10, the power flow restriction value of the branch, the operation state of the generator and its connection destination node number, the fixed ratio of the load and the fluctuation amount Read the ratio. When the “input data from disk” selection button 246 is selected, a flag indicating the operation status of the branch in the system connection status storage area 1026, the upper and lower limits of the power flow restriction of the branch, the connection destination node number of the generator, and the fixed load. The value indicating the operating state of the generator is read from the load factor, the value of the load factor for the variation, and the start / stop state storage area 1024 for the generator. With these read values, the time is input to the “display time” input area 247 and the “data reflection” button 266 is selected, so that the connection state of the branch at the specified time is displayed in the area 249. Here, the branch is a transmission and transformation facility, such as a transmission line and a transformer. The used branch is indicated by a solid line, and the stopped branch is indicated by a dotted line. A black circle indicates a node and indicates a bus or the like. To correct the connection state displayed in the area 249, the value of the “operation state” area 251 having the characteristics of the input area is changed via the input device 6, and the “data reflection” button 266 is selected.
The displayed numerical values are input to the areas 251, 252, 253, 257, and 257 having the characteristics of the input area via the input device 6.
The values of 260 and 261 can be modified. In the present embodiment, in the display of the area 249, the positive direction of the tidal current is from left to right and from the center to the outside. For this reason, the lower limit 251 of the power flow constraint is a negative value. In the area 249, the branch number is indicated by prefixing L to the number. In the area 249, the node number is indicated by adding N to the beginning of the number. The load and the number of the node to which the load connects are the same. The load of each node is a variable load obtained by subtracting the sum of the fixed load of the load ratio from the sum of the demand and the transmission loss. The value obtained by multiplying the load ratio of each node with respect to the sum of the load ratios by this variable load is the variable load of each node. The value obtained by adding the variable load and the fixed load is the load of each node at each time. Although power transmission loss occurs in each branch, the power transmission loss is incorporated in the load here.
【００９４】また、「コピー元の時刻」入力領域２６２
に時刻を入力し、「コピー先開始時刻」入力領域２６３
と「コピー先終了時刻」入力領域２６４に時刻を入力し
て、「コピー開始」ボタン２６５を選択すると、コピー
元の時刻の入力領域の特性を持つ「運用状態」領域２５
１，「潮流下限」領域２５２，「潮流上限」領域２５３，
「接続先ノード」領域２５７，「負荷比率固定」領域２
６０および「負荷比率変動」領域２６１のデータがコピ
ー先の時刻のデータ領域にコピーされる。「データ反
映」ボタン２６６を選択すると、「表示時刻」２４７の
ブランチの運用状態２５１を表示領域２４９に表示す
る。Also, a "copy source time" input area 262
In the “copy destination start time” input area 263
When the user inputs the time in the “copy destination end time” input area 264 and selects the “copy start” button 265, the “operation state” area 25 having the characteristics of the input time of the copy source time is displayed.
1, the “lower tide” area 252, the “upper tide” area 253,
"Connection destination node" area 257, "Load ratio fixed" area 2
60 and the data in the “load ratio fluctuation” area 261 are copied to the data area at the copy destination time. When the “data reflection” button 266 is selected, the operation state 251 of the “display time” 247 branch is displayed in the display area 249.
【００９５】スクロール２５４，２５８を選択すること
により、その選択した方向に応じて、ブランチ番号，発
電機名称を上，下へスクロールさせて、表示対象の領域
とすることができる。表示，入力領域２５１，２５２，
２５３，２５７，２６０，２６１の値は、「データ格
納」ボタン２６７を選択することにより、系統接続状態
記憶領域１０２６に格納される。「設定終了」２６８の
領域が選択されると、データ設定部１３は、系統接続状
態の入力，設定を終了し、処理選択メニュー画面２０に
戻る。By selecting the scrolls 254 and 258, the branch number and the generator name can be scrolled up and down in accordance with the selected direction to make the display area. Display, input areas 251, 252,
The values of 253, 257, 260, and 261 are stored in the system connection state storage area 1026 by selecting a "data storage" button 267. When the “setting end” area 268 is selected, the data setting unit 13 ends the input and setting of the system connection state, and returns to the process selection menu screen 20.
【００９６】ｉ．系統設備の運用特性の入力，設定 処理選択メニュー画面２０において、系統設備の運用特
性の入力，設定（項目ｉ）が入力装置６を介して選択さ
れると、データ設定部１３は処理対象のデータベース８
の計画対象期間記憶領域１００の計画対象年月日と、制
約条件記憶領域１０２の系統設備運用特性記憶領域１０
２７のブランチのインピーダンス特性,潮流値に対する
送電ロス係数およびブランチの潮流制約違反に対する重
み係数と系統接続状態記憶領域１０２６からブランチの
運用状態のフラグを読み込む。表示時刻はデフォルトの
１が最初に入力され、この時刻の系統接続状態記憶領域
１０２６と系統設備運用特性記憶領域１０２７にデータ
が格納されているときは、格納されている値を、格納さ
れていないときはデフォルト値である１が入力された設
定画面２７６が表示される。ここで、表示される画面２
７６を図１５に示す。この画面２７６は、計画対象期間
を示す表示領域２７０と、系統設備運用特性状態を読み
込む元である「サーバからのデータ入力」選択ボタン２
７１と「ディスクからのデータ入力」選択ボタン２７２
を表示し、系統接続図表示領域２７５と、サーバまたは
ディスクから読み込んだデータの表示をさせるための
「表示」ボタン２７４と表示する時刻を入力する「表示
時刻」入力領域２７３とコピー元の時刻入力領域２８５
とコピー先の開始時刻入力領域２８６とコピー先の終了
時刻入力領域２８７とコピーすることを指示する「コピ
ー開始」ボタン２８８と「データ反映」ボタン２８４と
ブランチの特性領域２７６とブランチのインピーダンス
特性の実部の入力，表示領域２７７とブランチのインピ
ーダンス特性の虚部の入力，表示領域２７８とブランチ
のロス係数の２次係数の入力，表示領域２７９とブラン
チのロス係数の１次係数の入力，表示領域２８０とブラ
ンチのロス係数の定数の入力，表示領域２８１と潮流制
約違反時のペナルティコストのための潮流制約違反の重
み係数の入力，表示領域２８２とブランチの特性をブラ
ンチ番号方向にスクロールさせるためのスクロールボタ
ン２８３を備え、さらに「データ格納」選択ボタン２８
９と「設定終了」選択ボタン２９０の領域がある。「サ
ーバからのデータ入力」選択ボタン２７１を選択すると
電力システム１０の対象年月日の１日分のブランチのイ
ンピーダンス特性，送電ロス特性および潮流制約違反時
のペナルティの重み係数とブランチの運用状態のフラグ
を読み込む。また、「ディスクからのデータ入力」選択
ボタン２７２を選択すると系統接続状態記憶領域１０２
６のブランチの運用状態のフラグと系統設備運用特性記
憶領域１０２７からブランチのインピーダンス特性の実
部と虚部、ブランチのロス係数の２次係数，１次係数，
定数と潮流制約違反時のペナルティコストのための潮流
制約違反の重み係数を読み込む。これらの読み込んだ値
は、「表示時刻」入力領域２７３に時刻を入力して、
「データ反映」ボタン２８４を選択することにより、表
示領域２７５に対象時刻のブランチの接続状態を表示す
る。また、「表示時刻」入力領域２７３に時刻を入力し
て、「表示」ボタン２７４を選択することにより、入
力，表示領域２７６に対象時刻のブランチの特性等を表
示する。また、表示されている数値は入力装置６を介し
て入力領域の特性を持つ領域２７７，２７８，２７９，
２８０，２８１，２８２の値を修正することができる。
重み係数２８２は出力計算部１５で制約違反を満たす解
が存在しないときに、制約違反の量に応じて目的関数に
ペナルティコストとして、加算するときに使用する。ペ
ナルティコストは違反量に重み係数２８２と図８の潮流
制約の重み係数４５を掛けた値となる。重み係数２８２
はブランチの中でも、電圧階級が高いブランチの潮流制
約が重要であるとき、その重み係数を大きくすることに
より、他の潮流制約違反より優先して違反を解消しやす
くするためのものであり、この値が大きいブランチほど
制約違反が解消しやすくなる。送電ロス係数は、出力計
算部１５で出力を計算し、各ブランチの潮流値が計算で
きたとき、この潮流に基づいてロス係数を用いた２次式
により各ブランチの送電ロスを計算するために使用す
る。この計算した送電ロスの和が全体の送電ロスとな
る。各時間の全体の送電ロスの初期値は需要と図１３の
初期ロス率９８の積により計算する。また、「コピー元
の時刻」入力領域２８５に時刻を入力し、「コピー先開
始時刻」入力領域２８６と「コピー先終了時刻」入力領
域２８７に時刻を入力して、「コピー開始」ボタン２８
８を選択すると、コピー元の時刻の入力領域の特性を持
つ「インピーダンス特性の実部」領域２７７，「インピ
ーダンス特性の虚部」領域２７８，「ブランチの２次の
ロス係数」領域２７９，「ブランチの１次のロス係数」
領域２８０，「ブランチの定数のロス係数」領域２８
１，「潮流制約違反時のペナルティの重み係数」領域２
８２のデータがコピー先の時刻のデータ領域にコピーさ
れる。「データ反映」ボタン２８４を選択すると、「表
示時刻」２７４のブランチの運用状態を表示領域２７５
に表示する。I. Input and setting of operation characteristics of system equipment When input and setting (item i) of operation characteristics of system equipment are selected via the input device 6 on the process selection menu screen 20, the data setting unit 13 sets the database to be processed. 8
Of the planning target period storage area 100, and the system equipment operation characteristic storage area 10 of the constraint condition storage area 102.
27, the branch characteristic is read from the system connection state storage area 1026, and the impedance characteristics of the 27 branches, the power transmission loss coefficient for the power flow value, the weight coefficient for violation of the power flow restriction of the branch, and the system connection state storage area 1026. As the display time, a default of 1 is input first, and when data is stored in the system connection state storage area 1026 and the system equipment operation characteristic storage area 1027 at this time, the stored value is not stored. At this time, a setting screen 276 in which 1 as a default value is input is displayed. Here, the displayed screen 2
76 is shown in FIG. This screen 276 includes a display area 270 indicating a planning target period and a “data input from server” selection button 2 from which the system equipment operation characteristic state is read.
71 and "Data input from disk" selection button 272
Is displayed, a system connection diagram display area 275, a "display" button 274 for displaying data read from a server or a disk, a "display time" input area 273 for inputting a display time, and a copy source time input Region 285
A "copy start" button 288, a "data reflection" button 284 for instructing to copy, a branch characteristic area 276, and the impedance characteristic of the branch are input to the start time input area 286 of the copy destination and the end time input area 287 of the copy destination. Input of real part, input of imaginary part of impedance characteristic of display area 277 and branch, input of second order coefficient of display area 278 and loss coefficient of branch, input and display of first order coefficient of display area 279 and loss coefficient of branch To input the constant of the loss coefficient of the area 280 and the branch, input the display area 281 and the weight coefficient of the power flow constraint violation for penalty cost when the power flow restriction is violated, and scroll the display area 282 and the characteristic of the branch in the branch number direction. Scroll button 283, and furthermore, a "data storage" selection button 28
9 and an “end setting” selection button 290. When the “input data from server” selection button 271 is selected, the impedance characteristic of the branch for one day on the target date of the power system 10, the transmission loss characteristic, the weighting factor of the penalty when the power flow constraint is violated, and the operation state of the branch are displayed. Read the flags. When the “input data from disk” selection button 272 is selected, the system connection status storage area 102
6, the real part and the imaginary part of the impedance characteristic of the branch, the second-order coefficient, the first-order coefficient of the loss coefficient of the branch, and the like.
Load constant and weight factor of power constraint violation for penalty cost when power constraint is violated. These read values are input into the “display time” input area 273 by inputting the time,
By selecting the “data reflection” button 284, the connection state of the branch at the target time is displayed in the display area 275. By inputting a time in the “display time” input area 273 and selecting a “display” button 274, the characteristics of the branch at the target time are displayed in the input / display area 276. Also, the displayed numerical values are input to the areas 277, 278, 279,
The values of 280, 281, 282 can be modified.
The weight coefficient 282 is used when the output calculation unit 15 does not find a solution that satisfies the constraint violation, and adds it as a penalty cost to the objective function according to the amount of the constraint violation. The penalty cost is a value obtained by multiplying the violation amount by the weight coefficient 282 and the weight coefficient 45 of the power flow constraint in FIG. Weight coefficient 282
When the power flow restriction of a branch with a high voltage class is important among branches, the weight coefficient is increased to make it easier to resolve the violation in preference to other power flow constraint violations. The larger the value of the branch, the easier it is to resolve the constraint violation. The power transmission loss coefficient is calculated by calculating the output by the output calculation unit 15 and calculating the power flow value of each branch by a quadratic expression using a loss coefficient based on the power flow when the power flow value of each branch can be calculated. use. The sum of the calculated power transmission losses is the total power transmission loss. The initial value of the entire power transmission loss at each time is calculated by the product of the demand and the initial loss rate 98 in FIG. Further, the user inputs a time in the “copy source time” input area 285, inputs a time in the “copy destination start time” input area 286 and a “copy destination end time” input area 287, and enters a “copy start” button 28.
When “8” is selected, the “real part of impedance characteristic” region 277, the “imaginary part of impedance characteristic” region 278, the “second-order loss coefficient of branch” region 279, and the “branch” having the characteristics of the input region of the copy source time are selected. First-order loss coefficient "
Region 280, “Branch constant loss coefficient” region 28
1, "Weight coefficient of penalty when power flow constraint is violated" area 2
The data 82 is copied to the data area at the copy destination time. When the “reflect data” button 284 is selected, the operation status of the branch of “display time” 274 is displayed in the display area 275.
To be displayed.
【００９７】スクロール２８３を選択することにより、
その選択した方向に応じて、ブランチ番号を上，下へス
クロールさせて、表示対象の領域とすることができる。
表示，入力領域２７７，２７８，２７９，２８０，２８
１，２８２の値は、「データ格納」ボタン２８９を選択
することにより、系統設備運用特性記憶領域１０２７に
格納される。「設定終了」２９０の領域が選択される
と、データ設定部１３は、系統設備運用特性の入力，設
定を終了し、処理選択メニュー画面２０に戻る。 ｊ．収束判定条件の入力，設定 処理選択メニュー画面２０において、収束判定条件の入
力，設定（項目ｊ）が入力装置６を介して選択される
と、データ設定部１３は処理対象のデータベース８の主
双対内点法のパラメータ記憶領域１０３の収束判定条件
記憶領域１０３１の最大反復回数，主問題の許容制約違
反量，双対問題の許容制約違反量，許容双対ギャップ
量，許容コストギャップ率を読み込む。主双対内点法の
パラメータ記憶領域１０３に格納されているデフォルト
値が設定画面２９１に表示される。ここで、表示される
画面２９１を図１６に示す。この画面２９１は、最大反
復回数の入力，表示領域２９２，主問題の許容制約違反
量の入力，表示領域２９３，双対問題の許容制約違反量
の入力，表示領域２９４，許容双対ギャップ量の入力，
表示領域２９５，許容コストギャップ率の入力，表示領
域２９６を備える。また、表示されている数値は入力装
置６を介して入力領域の特性を持つ領域２９２，２９
３，２９４，２９５，２９６の値を修正することができ
る。これらの値は、出力計算部１５により主双対内点法
で問題を解くときに使用する変数である。最大反復回数
２９２は主双対内点法で解を更新する最大回数である。
主問題の制約違反量２９３は、主問題の制約違反量であ
り、（数４）の第２番目の式の左辺のベクトルの各要素
の和である。双対問題の制約違反量２９４は、双対問題
の制約違反量であり、この式の第１番目の式の左辺のベ
クトルの各要素の和である。双対ギャップ２９５は、最
適解との最大偏差を示す指標であり、主問題の制約違反
量であり、この式の第５番目と第６番目の式のＳz とＵ
w のベクトルの各要素の和である。コストギャップ率２
９６は、双対ギャップを双対ギャップと目的関数の和で
割った値の比率（％）である。表示，入力領域２９２，
２９３，２９４，２９５，２９６の値は、「データ格
納」ボタン２９７を選択することにより、収束判定条件
記憶領域１０３１に格納される。「設定終了」２９８の
領域が選択されると、データ設定部１３は、収束判定条
件の入力，設定を終了し、処理選択メニュー画面２０に
戻る。By selecting the scroll 283,
According to the selected direction, the branch number can be scrolled up and down to be a display target area.
Display and input areas 277, 278, 279, 280, 28
By selecting the “data storage” button 289, the values of 1,282 are stored in the system facility operation characteristic storage area 1027. When the “setting end” area 290 is selected, the data setting unit 13 ends the input and setting of the system facility operation characteristics, and returns to the process selection menu screen 20. j. Input and Setting of Convergence Judgment Condition When input and setting (item j) of the convergence judgment condition are selected via the input device 6 on the process selection menu screen 20, the data setting unit 13 sets the main dual of the database 8 to be processed. The maximum number of iterations of the convergence determination condition storage area 1031 of the parameter storage area 103 of the interior point method, the allowable constraint violation amount of the main problem, the allowable constraint violation amount of the dual problem, the allowable dual gap amount, and the allowable cost gap ratio are read. The default value stored in the parameter storage area 103 of the primary dual interior point method is displayed on the setting screen 291. Here, the screen 291 displayed is shown in FIG. This screen 291 is used to enter the maximum number of iterations, display area 292, input the allowable constraint violation amount of the main problem, display area 293, input the allowable restriction violation amount of the dual problem, display area 294, input the allowable dual gap amount,
A display area 295 is provided for inputting an allowable cost gap rate and a display area 296. In addition, the displayed numerical values are input to the areas 292 and 29 having the characteristics of the input area via the input device 6.
3,294,295,296 values can be modified. These values are variables used when the output calculation unit 15 solves the problem by the primary dual interior point method. The maximum number of iterations 292 is the maximum number of times the solution is updated by the prime dual interior point method.
The constraint violation amount 293 of the main problem is the constraint violation amount of the main problem, and is the sum of the elements of the vector on the left side of the second expression of (Equation 4). The constraint violation amount 294 of the dual problem is the constraint violation amount of the dual problem, and is the sum of the respective elements of the vector on the left side of the first expression of this expression. The dual gap 295 is an index indicating the maximum deviation from the optimal solution, is a constraint violation amount of the main problem, and is Sz and U in the fifth and sixth equations of this equation.
The sum of each element of the vector w. Cost gap ratio 2
Reference numeral 96 denotes a ratio (%) of a value obtained by dividing the dual gap by the sum of the dual gap and the objective function. Display, input area 292,
The values of 293, 294, 295, and 296 are stored in the convergence determination condition storage area 1031 by selecting a “data storage” button 297. When the area of “setting end” 298 is selected, the data setting unit 13 ends the input and setting of the convergence determination condition, and returns to the process selection menu screen 20.
【００９８】ｋ．内点法制御変数の入力，設定 処理選択メニュー画面２０において、内点法制御変数の
入力，設定（項目ｋ）が入力装置６を介して選択される
と、データ設定部１３は処理対象のデータベース８の主
双対内点法のパラメータ記憶領域１０３のパラメータ記
憶領域１０３３の最大ステップサイズ上限，ステップサ
イズ係数，δ係数，μ分母を読み込む。主双対内点法の
パラメータ記憶領域１０３３に格納されているデフォル
ト値が設定画面２９９に表示される。ここで、表示され
る画面２９９を図１７に示す。この画面２９９は、最大
ステップサイズ上限の入力表示領域３００，ステップサ
イズ係数の入力表示領域３２６，δ係数の入力表示領域
３０１，μ分母の入力表示領域３０２を備える。また、
表示されている数値は入力装置６を介して入力領域の特
性を持つ領域３００，３２６，３０１，３０２の値を修
正することができる。これらの値は、出力計算部１５に
より主双対内点法で問題を解くときに使用する変数であ
る。最大ステップサイズ上限３００は更新するニュート
ン方向のベクトルによる現在の解が内点となる条件を満
たすステップ幅の上限値であり、ステップ幅がこれ以上
のときは、最大ステップ幅に置き換えて減少させる。ス
テップサイズ係数３２６は、ステップサイズに対して掛
ける係数である。δ係数３０１とμ分母（ｋの値）はμ
を計算するために使用する係数である。図１７のｎは
（数４）の第５番目の式のベクトルＳz の要素の数と第
６番目の式のベクトルＵｗのベクトルの要素の数の和で
ある。μはδ係数と双対ギャップの積をμ分母で割った
値である。表示，入力領域３００，３２６，３０１，３
０２の値は、「データ格納」ボタン３０３を選択すること
により、パラメータ記憶領域１０３３に格納される。
「設定終了」３０４の領域が選択されると、データ設定
部１３は、内点法制御変数の入力，設定を終了し、処理
選択メニュー画面２０に戻る。K. Input and Setting of Interior Point Method Control Variable When the input and setting (item k) of the interior point method control variable is selected via the input device 6 on the process selection menu screen 20, the data setting section 13 sets the database to be processed. 8, the maximum step size upper limit, the step size coefficient, the δ coefficient, and the μ denominator of the parameter storage area 1033 of the parameter storage area 103 of the 8 principal dual interior point method are read. The default value stored in the parameter storage area 1033 of the primary dual interior point method is displayed on the setting screen 299. Here, the displayed screen 299 is shown in FIG. This screen 299 includes an input display area 300 for the upper limit of the maximum step size, an input display area 326 for the step size coefficient, an input display area 301 for the δ coefficient, and an input display area 302 for the μ denominator. Also,
The displayed numerical values can modify the values of the regions 300, 326, 301, and 302 having the characteristics of the input region via the input device 6. These values are variables used when the output calculation unit 15 solves the problem by the primary dual interior point method. The maximum step size upper limit 300 is the upper limit value of the step width that satisfies the condition that the current solution based on the vector in the Newton direction to be updated is an inner point. If the step width is larger than this, the maximum step size is reduced. The step size coefficient 326 is a coefficient by which the step size is multiplied. δ coefficient 301 and μ denominator (value of k) are μ
Is the coefficient used to calculate In FIG. 17, n is the sum of the number of elements of the vector Sz of the fifth equation of (Equation 4) and the number of elements of the vector Uw of the sixth equation. μ is a value obtained by dividing the product of the δ coefficient and the dual gap by the μ denominator. Display, input area 300, 326, 301, 3
The value of 02 is stored in the parameter storage area 1033 by selecting a “data storage” button 303.
When the “setting end” 304 area is selected, the data setting unit 13 ends the input and setting of the interior point method control variables, and returns to the process selection menu screen 20.
【００９９】ｌ．内点法の初期値の入力，設定 処理選択メニュー画面２０において、内点法の初期値の
入力，設定（項目ｌ）が入力装置６を介して選択される
と、データ設定部１３は処理対象のデータベース８の主
双対内点法のパラメータ記憶領域１０３の初期値記憶領
域１０３４の供給力変数の決定方法の番号，供給力の実
変数初期値，スラック変数初期値，ラグランジェ変数の
決定方法の番号，ラグランジェ変数の初期値，双対変数
の決定方法の番号，双対変数初期値を読み込む。主双対
内点法の初期値記憶領域１０３４に格納されているデフ
ォルト値が設定画面３０５に表示される。ここで、表示
される画面３０５を図１８に示す。この画面３０５は、
供給力実変数の決定方法の選択ボタン３０６〜３１０
と、実変数初期値の入力，出力領域３１２とスラック変
数初期値の入力，表示領域３１１とラグランジェ変数の
決定方法の選択ボタン３１３，３１４とラグランジェ変
数初期値の入力，表示領域３１５と、双対変数の決定方
法の選択ボタン３１６，３１７と双対変数初期値の入
力，表示領域316とデータ格納選択領域３１９と設定終
了選択領域３２０とを備える。また、表示されている数
値，選択ボタンは入力装置６を介して入力領域の特性を
持つ領域３１２，３１１，３１５，３１８の値や選択ボ
タンを修正することができる。供給力実変数の決定方法
の選択ボタン３０６〜３１０はトルグスイッチになって
おり、選択されて黒丸となるとその他は白丸になる。等
λ配分は、等λ法により需給バランス制約を満たすよう
に発電機の出力を配分した値を初期値として使用するこ
とを表している。また、出力上限とは、発電機の出力上
限値を初期値として使用することを表している。出力中
間とは、発電機の出力上限と出力下限の平均値を初期値
として使用することを表してる。出力下限とは、発電機
の出力上限値を初期値として使用することを表してい
る。但し、出力下限が０のときは、内点とするために、
図１９の上側許容誤差３２２の値を使用する。実数変数
初期値を指定すると、入力，表示領域３１２に入力され
た値が初期値として使用される。スラック変数は初期値
３１１に入力された値を初期値として使用される。ラグ
ランジェ変数の決定方法の選択ボタン３１３，３１４は
トルグスイッチになっており、選択されて黒丸となると
その他は白丸になる。ラグランジェ変数の初期値３１３
が選択されると、入力，表示領域３１５に表示されてい
る値が初期値として使用される。双対問題の制約充足３
１４が選択されると、(数４)の第１式を満たす変数を計
算し、この値を使用する。双対変数の決定方法の選択ボ
タン３１６，３１７はトルグスイッチになっており、選
択されて黒丸となるとその他は白丸になる。双対変数の
初期値３１６が選択されると、入力，表示領域３１８に
表示されている値が初期値として使用される。双対問題
の制約充足３１７が選択されると、この式の第１式を満
たす変数を計算し、この値を使用する。供給力変数の決
定方法選択ボタン３０６〜３１０で選択された番号，供
給力の実変数初期値，スラック変数初期値，ラグランジ
ェ変数の決定方法選択ボタン３１３，３１４で選択され
た番号，ラグランジェ変数の初期値，双対変数の決定方
法選択ボタン３１６，３１７で選択された番号，双対変
数初期値は、「データ格納」ボタン３１９を選択するこ
とにより、初期値記憶領域１０３４に格納される。「設
定終了」３２０の領域が選択されると、データ設定部１
３は、内点法の初期値の入力，設定を終了し、処理選択
メニュー画面２０に戻る。L. Input and Setting of Initial Value of Interior Point Method When the input and setting of the initial value of the interior point method (item 1) are selected via the input device 6 on the processing selection menu screen 20, the data setting unit 13 Of the method of determining the supply force variable in the initial value storage area 1034 of the parameter storage area 103 of the primary dual interior point method of the database 8 of the above, the method of determining the actual variable initial value of the supply force, the initial value of the slack variable, and the method of determining the Lagrange variable Read the number, the initial value of the Lagrange variable, the number of the method for determining the dual variable, and the initial value of the dual variable. The default value stored in the initial value storage area 1034 of the primary dual interior point method is displayed on the setting screen 305. Here, the displayed screen 305 is shown in FIG. This screen 305
Selection buttons 306 to 310 for determining the actual supply power variable
Input of real variable initial value, output area 312 and input of slack variable initial value, display area 311 and selection buttons 313 and 314 of method for determining Lagrange variable, input of Lagrange variable initial value, display area 315, The system includes selection buttons 316 and 317 for determining a dual variable, input of a dual variable initial value, a display area 316, a data storage selection area 319, and a setting end selection area 320. The displayed numerical values and selection buttons can modify the values and selection buttons of the regions 312, 311, 315, and 318 having the characteristics of the input region via the input device 6. The selection buttons 306 to 310 for the method of determining the actual supply power variable are toggle switches, and when selected and turned black, the others become white circles. The equal λ distribution indicates that a value obtained by allocating the output of the generator so as to satisfy the supply-demand balance constraint by the equal λ method is used as an initial value. The output upper limit indicates that the output upper limit value of the generator is used as an initial value. The term “intermediate output” indicates that the average value of the upper and lower limits of the output of the generator is used as an initial value. The output lower limit indicates that the output upper limit value of the generator is used as an initial value. However, when the output lower limit is 0, in order to set the inside point,
The value of the upper tolerance 322 in FIG. 19 is used. When a real variable initial value is specified, the value input to the input / display area 312 is used as the initial value. As the slack variable, the value input to the initial value 311 is used as an initial value. The selection buttons 313 and 314 of the Lagrange variable determination method are toggle switches, and when selected and turned to a black circle, the others become white circles. Initial value of Lagrange variable 313
Is selected, the value displayed in the input / display area 315 is used as an initial value. Constraint satisfaction of dual problem 3
When 14 is selected, a variable that satisfies the first expression of (Equation 4) is calculated, and this value is used. The selection buttons 316 and 317 for the method of determining the dual variable are toggle switches, and when they are selected and become black circles, others become white circles. When the initial value 316 of the dual variable is selected, the value displayed in the input / display area 318 is used as the initial value. When the constraint satisfaction 317 of the dual problem is selected, a variable that satisfies the first expression of this expression is calculated, and this value is used. Number selected with supply power variable determination method selection buttons 306 to 310, actual variable initial value of supply power, initial slack variable, number selected with Lagrange variable determination method selection buttons 313, 314, Lagrange variable , The number selected by the dual variable determination method selection buttons 316 and 317, and the dual variable initial value are stored in the initial value storage area 1034 by selecting the "data storage" button 319. When the area of “setting end” 320 is selected, the data setting unit 1
Step 3 terminates the input and setting of the initial value of the interior point method, and returns to the process selection menu screen 20.
【０１００】ｍ．等式制約の許容偏差の入力，設定 処理選択メニュー画面２０において、等式制約の許容偏
差の入力，設定（項目ｍ）が入力装置６を介して選択さ
れると、データ設定部１３は処理対象のデータベース８
の主双対内点法のパラメータ記憶領域１０３の上限，下
限許容範囲記憶領域１０３２の上限側許容偏差と下限側
許容偏差の値を読み込む。主双対内点法の上限，下限許
容範囲記憶領域１０３２に格納されているデフォルト値
が設定画面３２１に表示される。ここで、表示される画
面３２１を図１９に示す。この画面３２１は、上限側許
容偏差の入力，表示領域３２２と下限側許容偏差の入
力，表示領域３２３とデータ格納選択領域３２４と設定
終了選択領域３２５を備える。また、表示されている数
値は入力装置６を介して入力領域の特性を持つ領域３２
２，３２３の値を修正することができる。上側許容偏差
は等式制約の定数項に上側許容偏差を加算して、この加
算値以下として不等式制約とするために使用する値であ
る。また、下側許容偏差は等式制約の定数項に下側許容
偏差を減算して、この減算値以上として不等式制約とす
るために使用する値である。値は、０以上の値とする。
上限側許容偏差３２２と下限側許容偏差３２３の値は、
「データ格納」ボタン３２４を選択することにより、上
限，下限許容範囲記憶領域1032に格納される。「設定終
了」３２５の領域が選択されると、データ設定部１３
は、等式制約の許容偏差の入力，設定を終了し、処理選
択メニュー画面２０に戻る。M. Input and Setting of Permissible Deviation of Equation Constraint When the input and setting (item m) of the permissible deviation of equation constraint are selected via the input device 6 on the process selection menu screen 20, the data setting unit 13 Database 8
The values of the upper and lower allowable tolerances of the upper and lower allowable range storage areas 1032 of the parameter storage area 103 of the primary dual interior point method are read. The default value stored in the upper / lower allowable range storage area 1032 of the primary dual interior point method is displayed on the setting screen 321. Here, the screen 321 to be displayed is shown in FIG. The screen 321 includes an upper limit side allowable deviation input, a display area 322 and a lower limit side allowable deviation input / display area 323, a data storage selection area 324, and a setting end selection area 325. Further, the displayed numerical value is input to the area 32 having the characteristics of the input area via the input device 6.
The value of 2,323 can be modified. The upper allowable deviation is a value used to add the upper allowable deviation to the constant term of the equality constraint and set the inequality constraint to be equal to or less than the added value. Further, the lower allowable deviation is a value used for subtracting the lower allowable deviation from the constant term of the equality constraint to make the inequality constraint equal to or greater than the subtracted value. The value is a value of 0 or more.
The values of the upper limit tolerance 322 and the lower limit tolerance 323 are
By selecting the “store data” button 324, the data is stored in the upper limit / lower limit allowable range storage area 1032. When the area of “setting end” 325 is selected, the data setting unit 13
Ends the input and setting of the allowable deviation of the equation constraint, and returns to the process selection menu screen 20.
【０１０１】ｎ．制約マトリックス分割方法の入力，設
定 処理選択メニュー画面２０において、制約マトリックス
分割方法の入力，設定（項目ｍ）が入力装置６を介して
選択されると、データ設定部１３は処理対象のデータベ
ース８の主双対内点法のパラメータ記憶領域１０３の制
約マトリックス分割方法記憶領域１０３５の分割の組合
せ番号，時間軸方向の分割数、上下限制約への変更方法
の番号を読み込む。制約マトリックス分割方法記憶領域
１０３５に格納されているデフォルト値が設定画面３２
７に表示される。ここで、表示される画面３２７を図２
０に示す。この画面３２７は、分割の組合せ方法の選択
ボタン３２８〜３３２と、時間軸方向の分割数の入力，
出力領域３３３と時間で分類した制約マトリックスの構
成表示領域３３４と上下限制約への変更方法選択ボタン
３３５〜３４２と等式，不等式による制約の分類表示領
域３４３とデータ格納選択領域３４５と設定終了選択領
域３４６とを備える。また、表示されている数値，選択
ボタンは入力装置６を介して入力領域の特性を持つ領域
３３３の値や選択ボタン３２８〜３３２の選択番号と選
択ボタン３３５〜３４２の選択番号を修正することがで
きる。分割方法の選択ボタン３２８〜３３２はトルグス
イッチになっており、選択されて黒丸となるとその他は
白丸になる。１番を選択すると、時間断面制約，差分制
約，積分制約の３つの行方向に制約マトリックスを分割
することを表す。２番を選択すると、時間断面制約と差
分制約を合わせて１つの行方向とし、積分制約を１つの
行方向として、２つの行方向に制約マトリックスを分割
することを表す。３番を選択すると、時間断面制約を１
つの行方向とし、差分制約と積分制約を合わせて１つの
行方向として、２つの行方向に制約マトリックスを分割
することを表す。４番を選択すると、時間断面制約と積
分制約を合わせて１つの行方向とし、差分制約を１つの
行方向として、２つの行方向に制約マトリックスを分割
することを表す。５番を選択すると、時間断面制約と積
分制約と差分制約を合わせて１つの行方向、すなわち、
行方向には分割しないことを表す。分割数入力表示領域
３３３の数値は、時間断面方向（列方向）に制約マトリ
ックスを分割するときの数を表す。時間で分類した制約
マトリックスの構成表示領域３３４は制約マトリックス
の非ゼロ要素ブロック行列とゼロ行列の配置構造を表し
ている。本実施例では、制約マトリックスと掛ける変数
ベクトルの要素を並べる順序を１時の発電機１の出力，
１時の発電機２の出力，１時の発電機３の出力，１時の
発電機４の出力，２時の発電機１の出力，２時の発電機
２の出力，２時の発電機３の出力，２時の発電機４の出
力，３時の発電機１の出力，３時の発電機２の出力，３
時の発電機３の出力，３時の発電機４の出力とうい順序
としている。このため、潮流制約，需給バランス制約等
の時間断面制約は時間毎に固まった行列（ハッチング部
分）が非ゼロ要素がある部分となる。時間断面制約を示
す領域で、空白の行列部分はゼロ行列である。時間断面
制約の１つのブロック行列Ａ１，Ａ２，Ａ３は、N. Input and Setting of Constraint Matrix Division Method When input and setting (item m) of the restriction matrix division method are selected via the input device 6 on the processing selection menu screen 20, the data setting unit 13 sets the data of the database 8 to be processed. The combination number of the division in the constraint matrix division method storage area 1035 of the parameter storage area 103 of the primary dual interior point method, the number of divisions in the time axis direction, and the number of the change method to the upper and lower limit restrictions are read. The default value stored in the constraint matrix division method storage area 1035 is the setting screen 32
7 is displayed. Here, the displayed screen 327 is shown in FIG.
0 is shown. This screen 327 includes selection buttons 328 to 332 for a combination method of division, input of the number of divisions in the time axis direction,
The output area 333, the configuration display area 334 of the restriction matrix classified by time, the change method selection buttons 335 to 342 to the upper and lower restriction, the classification display area 343 of the restriction by the equation and the inequality, the data storage selection area 345, and the setting end selection And an area 346. The displayed numerical values and selection buttons can be used to correct the values of the area 333 having the characteristics of the input area, the selection numbers of the selection buttons 328 to 332, and the selection numbers of the selection buttons 335 to 342 via the input device 6. it can. The division method selection buttons 328 to 332 are toggle switches, and when selected and turned into a black circle, the others become white circles. Selecting No. 1 indicates that the constraint matrix is divided in the three row directions of the time section constraint, the difference constraint, and the integral constraint. When No. 2 is selected, the time section constraint and the difference constraint are combined into one row direction, and the integration constraint is set as one row direction, and the constraint matrix is divided into two row directions. Select No. 3 to set the time section constraint to 1
This means that the constraint matrix is divided into two row directions, and the difference matrix and the integration constraint are combined into one row direction. When No. 4 is selected, the time section constraint and the integration constraint are combined into one row direction, and the difference constraint is set as one row direction, and the constraint matrix is divided into two row directions. When No. 5 is selected, the time section constraint, the integration constraint, and the difference constraint are combined in one row direction, that is,
Represents no division in the row direction. The numerical value in the division number input display area 333 represents the number when dividing the constraint matrix in the time section direction (column direction). The configuration display area 334 of the constraint matrix classified by time represents the arrangement structure of the non-zero element block matrix and the zero matrix of the constraint matrix. In the present embodiment, the order of arranging the elements of the constraint matrix and the variable vector to be multiplied is the output of the generator 1 at 1 o'clock,
1 o'clock generator 2 output, 1 o'clock generator 3 output, 1 o'clock generator 4 output, 2 o'clock generator 1 output, 2 o'clock generator 2 output, 2 o'clock generator Output at 3 o'clock, output of generator 4 at 2 o'clock, output of generator 1 at 3 o'clock, output of generator 2 at 3 o'clock, 3 o'clock
The output of the generator 3 at the time and the output of the generator 4 at 3 o'clock are in order. For this reason, in the time cross section constraint such as the tidal flow constraint and the supply / demand balance constraint, a matrix (hatched portion) that is fixed for each time is a portion having a non-zero element. In the region indicating the time section constraint, a blank matrix portion is a zero matrix. One block matrix A1, A2, A3 of the time section constraint is
【０１０２】[0102]
【数２６】 (Equation 26)
【０１０３】のようになる。この式の第１行は需給バラ
ンス制約で、全ての発電機の出力の和であることを示し
ている。揚水ポンプ、融通送電は要素の値が−１とな
る。２行目以降は潮流制約の例を示している。系統構成
が同一であれば、時間断面制約の各時刻のブロック行列
は同一となる。出力変化制約は差分型の制約で、前後の
時間との出力の差なので、非ゼロ要素は表示領域３３４
の差分制約の構造になる。全ての発電機で出力変化速度
制約がある場合はブロック行列Ｂ１，Ｂ３は、Is as follows. The first line of this equation is a supply-demand balance constraint, and indicates that it is the sum of the outputs of all the generators. The element value of the pump and the power transmission is -1. The second and subsequent lines show examples of power flow restrictions. If the system configuration is the same, the block matrix at each time of the time section constraint is the same. The output change constraint is a difference-type constraint, and the difference between the output before and after the difference is the non-zero element.
Of the difference constraint. If there is an output change speed constraint in all generators, the block matrices B1 and B3 are
【０１０４】[0104]
【数２７】 [Equation 27]
【０１０５】のように単位行列になり、行列Ｂ２，Ｂ４
は単位行列に−１を掛けたものとなる。出力変化速度制
約がない発電機の出力を変数ベクトルの要素で並べる順
序を最後に持ってくると、ブロック行列Ｂ１，Ｂ３はこ
の式のように単位行列とゼロ行列になる。As shown in the following, a unit matrix is obtained, and matrices B2 and B4
Is the unit matrix multiplied by -1. When the order of arranging the outputs of the generator without the output change speed constraint by the elements of the variable vector is finally brought, the block matrices B1 and B3 become a unit matrix and a zero matrix as in this equation.
【０１０６】[0106]
【数２８】 [Equation 28]
【０１０７】は２つの発電機の出力変化速度制約がない
ときの例を示している。池水位制約は積分型の制約で、
各時刻までの発電機の出力の和が池水位となる。よっ
て、ブロック行列Ｃ１〜Ｃ６は、FIG. 11 shows an example in which there is no restriction on the output change speed of the two generators. The pond water level constraint is an integral type constraint,
The sum of the generator output up to each time is the pond water level. Therefore, the block matrices C1 to C6 are
【０１０８】[0108]
【数２９】 (Equation 29)
【０１０９】とゼロ行列と単位行列となる。（数２９）
は水力発電機が３機の場合の例である。但し、揚水発電
により低下する水位は出力を揚水効率で割った値なの
で、この揚水発電の行の値は１ではなくなる。しかし、
行列の構造は対角要素のみが非ゼロとなる対角行列のま
まである。同一揚水発電機において揚水効率が時間によ
らず同一であれば、ブロック行列Ｃ１〜Ｃ６は全て同じ
になる。また、積分制約として、最終時刻でのみ制約が
あるＬＮＧ燃料消費制約やＩＰＰ利用率がある。これら
の制約は、Ｃ１〜Ｃ３がなく、Ｃ４〜Ｃ６が、（数２
３）の構造の行列となる。以上のように制約マトリック
スをブロック化すると、同一の行列や単位行列が現れる
ので、制約マトリックスの非ゼロ要素のみを記憶するよ
りも、さらに必要なメモリを減少できる。また、出力計
算部１５のマトリックスの計算でブロックの行列の転置
行列とブロックの行列の積がでてくるが、上記の差分制
約や積分制約のブロック行列はゼロ行列を除くと単位行
列や対角行列となるケースがあるので、これらのブロッ
クの行列の転置行列とブロックの行列の積は対角行列と
なることにより、行列の積の計算に要する時間を短縮で
きる。Then, a zero matrix and a unit matrix are obtained. (Equation 29)
Is an example in the case of three hydraulic power generators. However, since the water level lowered by the pumped-storage power generation is a value obtained by dividing the output by the pumping efficiency, the value of the line of the pumped-storage power generation is not 1. But,
The structure of the matrix remains a diagonal matrix in which only the diagonal elements are non-zero. If the pumping efficiencies are the same regardless of time in the same pumping generator, all the block matrices C1 to C6 are the same. Further, as integration constraints, there are LNG fuel consumption constraints and IPP utilization rates that are limited only at the last time. These restrictions are such that C1 to C3 do not exist and C4 to C6 are (Equation 2)
The matrix has the structure of 3). As described above, when the constraint matrix is divided into blocks, the same matrix and unit matrix appear, so that the required memory can be further reduced as compared with storing only the non-zero elements of the constraint matrix. The product of the transposed matrix of the block matrix and the matrix of the block is obtained in the calculation of the matrix by the output calculation unit 15. However, the block matrix of the difference constraint and the integration constraint is a unit matrix and a diagonal matrix except for the zero matrix. Since there is a case where the matrix is used, the product of the transposed matrix of these blocks and the matrix of the blocks is a diagonal matrix, so that the time required for calculating the matrix product can be reduced.
【０１１０】上記の例では、制約マトリックスと掛ける
変数ベクトルの要素を並べる順序を１時の発電機１の出
力，１時の発電機２の出力，１時の発電機３の出力，１
時の発電機４の出力，２時の発電機１の出力，２時の発
電機２の出力，２時の発電機３の出力，２時の発電機４
の出力，３時の発電機１の出力，３時の発電機２の出
力，３時の発電機３の出力，３時の発電機４の出力とい
う順序としたが、変数ベクトルの要素を並べる順序を１
時の発電機１の出力，２時の発電機１の出力，３時の発
電機１の出力，１時の発電機２の出力，２時の発電機２
の出力，３時の発電機２の出力，１時の発電機３の出
力，２時の発電機３の出力，３時の発電機３の出力，１
時の発電機４の出力，２時の発電機４の出力，３時の発
電機４の出力となる順序でならべて、ブロック行列を、
時間断面制約，差分制約，積分制約に分割する方法もあ
る。また、行列ブロックＡｉ，Ｂｉ，Ｃｉをさらに分割
することもできる。In the above example, the order of arranging the elements of the variable vector to be multiplied by the constraint matrix is the output of the generator 1 at 1 o'clock, the output of the generator 2 at 1 o'clock, the output of the generator 3 at 1 o'clock, and 1
Output of generator 4 at time, output of generator 1 at 2 o'clock, output of generator 2 at 2 o'clock, output of generator 3 at 2 o'clock, generator 4 at 2 o'clock
, The output of the generator 1 at 3 o'clock, the output of the generator 2 at 3 o'clock, the output of the generator 3 at 3 o'clock, and the output of the generator 4 at 3 o'clock, but the elements of the variable vector are arranged. Order 1
Output of generator 1 at time, output of generator 1 at 2 o'clock, output of generator 1 at 3 o'clock, output of generator 2 at 1 o'clock, generator 2 at 2 o'clock
Output, output of generator 2 at 3 o'clock, output of generator 3 at 1 o'clock, output of generator 3 at 2 o'clock, output of generator 3 at 3 o'clock, 1
The output of the generator 4 at time, the output of the generator 4 at 2 o'clock, and the output of the generator 4 at 3 o'clock are arranged in this order.
There is also a method of dividing into time section constraints, difference constraints, and integration constraints. Further, the matrix blocks Ai, Bi, and Ci can be further divided.
【０１１１】上下限制約への変更選択ボタン３３５〜３
４２はトルグスイッチになっており、選択されて黒丸と
なるとその他は白丸になる。１番を選択すると、全ての
制約を上下限制約とすることを表す。２番を選択する
と、等式制約と下限制約を上下限制約に変更することを
表している。３番を選択すると、下限制約と上限制約を
上下限制約に変更することを表している。４番を選択す
ると、等式制約と上限制約を上下限制約に変更すること
を表す。５番を選択すると、等式制約を上下限制約に変
更することを表す。６番を選択すると、下限制約を上下
限制約に変更することを表す。７番を選択すると、上限
制約を上下限制約に変更することを表す。８番を選択す
ると、新たに上下限制約に変更しないことを表す。上下
限制約は、非負のスラック変数を導入することにより、
スラック変数が０以上、上限と下限の差以下となるの
で、制約マトリックスに必要なのは、スラック変数の上
限制約のみとなり、制約マトリックスの１行のみを必要
とすることがわかる。元の上限，下限制約をそれぞれ、
１行ずつで数式化して、制約マトリックスに入れると、
制約マトリックスの行数が増加するが計算上は問題は生
じない。Select button 335-3 for change to upper / lower limit constraint
Reference numeral 42 denotes a toggle switch, and when a black switch is selected and the other switches are turned, the other switches become open circles. Selecting No. 1 indicates that all constraints are upper and lower limit constraints. Selecting No. 2 indicates that the equality constraint and the lower limit constraint are changed to upper and lower limit constraints. Selecting No. 3 indicates that the lower and upper constraints are changed to upper and lower constraints. Selecting No. 4 indicates that the equality constraint and the upper limit constraint are changed to upper and lower limit constraints. Selecting No. 5 indicates that the equation constraint is changed to upper and lower limit constraints. Selecting No. 6 indicates that the lower limit constraint is changed to the upper and lower limit constraint. Selecting No. 7 indicates that the upper limit constraint is changed to the upper and lower limit constraint. When No. 8 is selected, it means that the upper / lower limit is not newly changed. Upper and lower bound constraints are introduced by introducing a non-negative slack variable,
Since the slack variable is equal to or greater than 0 and equal to or less than the difference between the upper limit and the lower limit, it is understood that only the upper limit constraint of the slack variable is required for the constraint matrix, and only one row of the constraint matrix is required. The original upper and lower constraints, respectively,
Formulating each line one by one and putting it in the constraint matrix,
Although the number of rows in the constraint matrix increases, there is no problem in calculation.
【０１１２】分割の組合せ選択ボタン３２８〜３３２で
選択された番号，分割数３３３，上下限制約への変更方
法選択ボタン３３５〜３４２で選択された番号は、「デ
ータ格納」ボタン３４５を選択することにより、制約マ
トリックス分割方法記憶領域１０３５に格納される。
「設定終了」３４６の領域が選択されると、データ設定
部１３は、制約マトリックス分割方法の入力，設定を終
了し、処理選択メニュー画面２０に戻る。The "data storage" button 345 is selected for the number selected by the division combination selection buttons 328 to 332, the number of divisions 333, and the number selected by the upper / lower limit change method selection buttons 335 to 342. Is stored in the constraint matrix division method storage area 1035.
When the “setting end” area 346 is selected, the data setting unit 13 ends the input and setting of the constraint matrix division method, and returns to the process selection menu screen 20.
【０１１３】（２）出力の計算 計画作成条件が全て設定されると、発電機等の出力を計
算する。処理選択メニュー画面２０において、出力計算
が入力装置６を介して選択されると、データ設定部１３
は処理対象のデータベース８の制約条件記憶領域１０２
と目的関数記憶領域１０１と主双対内点法のパラメータ
記憶領域１０３の全てのデータを読み込む。次に、図１
の処理フローに従って、出力を計算する。処理２０００
は計画作成条件の設定であり、既に設定済みである。処
理２００１では、全ての制約条件を満たすかどうかのフ
ラグｆlgを０にセットする。最終結果まで、ｆlgが０な
らば、全ての制約条件を満たす出力が得られたことにな
り、ｆlgが１ならば、全ての制約条件を満たす出力がな
かったことを示すことになる。処理２００２は処理２０
００で設定した条件から２次計画問題を定式化し、内点
法により出力を計算する。目的関数の式，制約条件のマ
トリックスおよび内点法のパラメータを設定する。計画
対象期間全てで、停止している発電機は、出力決定の対
象外とする。停止している時刻がある発電機の出力上限
は０であるが、内点法を適用するためには、出力が正で
ある必要があるので、出力上限は、図１９の上側許容偏
差322の値を設定する。ここでは、目的関数は電力を調
達する費用を含めた燃料費として説明する。このとき、
目的関数は（数１９）となる。制約条件マトリックスの
分割方法は、図２０の分割の組合せボタン３３２が選択
され、分割数３３３が１のときを説明する。すなわち、
制約マトリックスは（数１６）の形になる。揚水，揚発
および水力の起動停止状態が１の場合は出力下限をここ
では０とする。揚水，揚発および水力の場合、運転可能
であっても必ずしも運転するとは限らないためである。
揚水，揚発および水力に限らず起動停止が決められない
発電機に対しては、一旦最低出力をゼロにしておいて、
出力配分結果により、最低出力以下であれば停止に変更
することもできる。等式制約の需給バランス（数１９）
は、供給力が増加すると、コストが増加するので、右辺
の需要と送電ロスの和以上に供給力が必要とする不等式
にできる。また、制約条件等により不必要に供給力が増
加しないように、需要と送電ロスに図１９で設定した上
限側許容値３２２を加算した値以下とする不等式に変え
るこれらを、スラック変数を用いて、等式制約に置換す
るものとする。（数５）のニュートン方程式は行列で表
すと、(2) Calculation of output When all the conditions for creating a plan are set, the output of the generator and the like is calculated. When the output calculation is selected via the input device 6 on the process selection menu screen 20, the data setting unit 13
Is the constraint storage area 102 of the processing target database 8
Then, all the data in the objective function storage area 101 and the parameter storage area 103 of the primary dual interior point method are read. Next, FIG.
The output is calculated according to the processing flow of (1). Processing 2000
Is the setting of the plan creation condition, which has already been set. In the process 2001, the flag flg of whether or not all the constraint conditions are satisfied is set to 0. Until the final result, if flg is 0, an output satisfying all constraints is obtained, and if flg is 1, it indicates that no output satisfies all constraints. Processing 2002 is processing 20
A quadratic programming problem is formulated from the conditions set at 00, and the output is calculated by the interior point method. Set the objective function formula, constraint matrix, and interior point method parameters. Generators that are stopped during the entire planning period are excluded from output determination. The upper limit of the output of the generator with the time at which it is stopped is 0, but the output must be positive in order to apply the interior point method. Set the value. Here, the objective function will be described as a fuel cost including the cost of procuring power. At this time,
The objective function is (Equation 19). The method of dividing the constraint matrix will be described when the combination button 332 in FIG. 20 is selected and the number of divisions 333 is one. That is,
The constraint matrix has the form of (Equation 16). When the pumping, pumping, and start / stop states of hydraulic power are 1, the output lower limit is set to 0 here. This is because, in the case of pumping, pumping, and hydraulic power, even if operable, operation is not necessarily performed.
For generators that cannot be started or stopped regardless of pumping, pumping and hydraulic power, set the minimum output to zero once,
According to the output distribution result, if the output is equal to or less than the minimum output, it can be changed to stop. Supply-demand balance of equality constraint (Equation 19)
Since the cost increases when the supply capacity increases, the inequality that the supply capacity requires more than the sum of the demand on the right side and the transmission loss can be obtained. Further, in order to prevent the supply power from unnecessarily increasing due to the constraint condition or the like, the demand and the transmission loss are changed to an inequality of not more than a value obtained by adding the upper limit allowable value 322 set in FIG. , Equality constraints. The Newton equation of (Equation 5) can be expressed by a matrix,
【０１１４】[0114]
【数３０】 [Equation 30]
【０１１５】となる。この式に（数１６）を代入する
と、Is obtained. By substituting (Equation 16) into this equation,
【０１１６】[0116]
【数３１】 (Equation 31)
【０１１７】となる。この式は図２のフローで解くこと
ができる。但し、ステップ１００５のバリアパラメータ
は図１７のδ係数３０１とμ分母３０２の値を用いるこ
ととする。処理１００４の収束条件として、図１６で設
定したコストギャップ率２９６以下となったときも終了
条件に追加する。処理１００７で制約マトリックスＡを
ブロックＡ０と単位行列に分けたので、ブロック化した
行列形式でニュートン方向ベクトルを計算すると、ｄｘ
が最初に計算できる。この値を代入することにより、ｄ
ｙ，ｄｓ，ｄｕ，ｄｚ，ｄｗが計算できる。処理１００
８でλは図１７のステップサイズ係数３２６の値を使用
し、αｐ，αｄの上限は、図１７の最大ステップサイズ
上限３００を越えない様にする。すなわち、上限以上と
なるときは、αｐ、αｄは上限値で置き換えるものとす
る。処理１０１１でｋが図１６の最大反回数２９２(＝
ｋmax）より大きいときは内点法による解法を終了とす
る。初期値は、図１８で設定した方法又は値を用いるも
のとする。供給力変数の実数値は選択した方法で計算ま
たは設定した値を使用する。スラック変数は指定した値
を使用する。ラグランジェ変数，双対変数は初期値で設
定した値か双対問題の制約（数４）の第１式である、Is obtained. This equation can be solved by the flow of FIG. However, the values of the δ coefficient 301 and the μ denominator 302 in FIG. 17 are used as the barrier parameters in step 1005. As the convergence condition of the process 1004, when the cost gap rate becomes equal to or less than the cost gap rate 296 set in FIG. Since the constraint matrix A is divided into the block A0 and the unit matrix in the processing 1007, when the Newton direction vector is calculated in the form of a block matrix, dx
Can be calculated first. By substituting this value, d
y, ds, du, dz, dw can be calculated. Processing 100
In FIG. 8, λ uses the value of the step size coefficient 326 in FIG. 17, and the upper limits of αp and αd do not exceed the maximum step size upper limit 300 of FIG. That is, when the value is equal to or more than the upper limit, αp and αd are replaced with the upper limit values. In the process 1011, k is the maximum number 292 (=
If it is larger than (kmax), the solution by the interior point method is terminated. As the initial value, the method or value set in FIG. 18 is used. The real value of the supply variable uses the value calculated or set by the selected method. The slack variable uses the specified value. The Lagrange variable and the dual variable are the values set by the initial values or the first expression of the constraint (equation 4) of the dual problem.
【０１１８】[0118]
【数３２】 (Equation 32)
【０１１９】が０となるように設定する。Is set to 0.
【０１２０】[0120]
【数３３】 [Equation 33]
【０１２１】を満たすように、ラグランジェ変数ｙ，双
対変数ｚ１，ｚ２，ｗ１，ｗ２を決めることにより、双
対問題の制約条件を満たす初期点から出発することがで
き、収束までの反復回数を低減することができる。By determining the Lagrange variable y and the dual variables z1, z2, w1, and w2 so as to satisfy the conditions, it is possible to start from the initial point satisfying the constraint condition of the dual problem, and reduce the number of iterations until convergence. can do.
【０１２２】制約マトリックス分割方法に応じて、（数
３０）の制約マトリックスを小ブロックに分けて計算を
実施するものとする。According to the constraint matrix division method, the calculation is performed by dividing the constraint matrix of (Equation 30) into small blocks.
【０１２３】処理２００４で、全ての制約条件を満たす
出力がなかったかどうかｆlgで判定して、制約緩和をし
ていないならば、処理２００５の最適解があったかどう
かチェックする。この判定は、処理２００２を抜けたの
が、処理１００４の収束条件を満たすかどうかチェック
することで最適解かどうか判定する。最適解であるなら
ば、処理２００６で、結果を図５の計算結果記憶領域１
０４に格納する。格納する項目は、反復毎の双対ギャッ
プ，バリアパラメータ，（数４）の第１式の双対制約違
反量、この式の第２式の主問題制約違反量および最適解
である。処理２００７では、処理２０００で設定した制
約条件に対して、制約違反があればその制約条件の名
称，違反が発生した時刻、および違反の超過量を図５の
計算結果記憶領域１０４に格納する。処理２００３で
は、運転可能である揚水，揚発，水力等に発電機の出力
が出力下限より小さいものがあるかどうかチェックす
る。なければ、出力計算の処理は終了して、図６の処理
選択メニュー画面２０に戻る。処理２００４でｆlgが１
のときは、制約を緩和した制約で、制約の上限または下
限となったものがあるかチェックする。これは、制約違
反のペナルティを目的関数に加算することにより、制約
違反を解消できた制約が存在するかどうか判定するため
のものである。制約違反を解消できた制約がなければ、
これ以上処理を続けても、良い解を得ることはできない
ので制約を緩和して、現状の制約を満たした解を結果と
して格納することになる。処理２００９のチェックによ
り、制約違反を解消できたものがあれば、処理２０１１
で、処理２０００で設定した制約条件でチェックして、
制約違反となっているものを抽出する。処理２０１２で
は、違反している制約の項目を格納する。さらに違反し
ている制約に対して制約条件を緩和するとともに、制約
違反のペナルティを設定する。ペナルティ係数は図８の
制約条件に対する重み係数４５と図１３の重み係数２１
５や図１５の重み係数２８２の設備に対する重み係数の
積を使用する。ペナルティにより、目的関数の線形項が
変化する。すなわち、線形の制約式のうち、発電機出力
の線形式のペナルティ係数を掛けた式を目的関数に加算
する。これにより、目的関数の線形式が変化する。ま
た、次の計算で、可能解を見つけるために、制約違反が
生じている制約に対しては制約値を緩和する。すなわ
ち、処理２０００で設定した制約値に対して、上限制約
を超過しているときは、制約値の下限値を上限値とし、
上限値をもとの上限値に違反量を加算したものにする。
また、処理２０００で設定した制約値に対して、下限制
約を超過しているときは、制約値の上限値を下限値と
し、下限値をもとの下限値に違反量を減算したものにす
る。これにより、緩和した制約に対しては最適解を求め
ることができる。次に処理２０１３では、制約の緩和を
実施したことを判定できるようにｆlgを１にする。再
度、処理２００２にて、緩和した制約に対する最適解を
求めることになる。In the process 2004, it is determined by flg whether or not there is an output satisfying all the constraint conditions. If the constraint is not relaxed, it is checked whether or not the process 2005 has an optimal solution. In this determination, it is determined whether or not the process exited from the process 2002 is an optimal solution by checking whether or not the convergence condition of the process 1004 is satisfied. If it is the optimal solution, the result is stored in the calculation result storage area 1 in FIG.
04. The items to be stored are the dual gap, the barrier parameter, the dual constraint violation amount of the first expression of (Equation 4), the main problem constraint violation amount of the second expression of this expression, and the optimal solution for each iteration. In the process 2007, if there is a constraint violation with respect to the constraint set in the process 2000, the name of the constraint, the time when the violation occurred, and the excess amount of the violation are stored in the calculation result storage area 104 in FIG. In the process 2003, it is checked whether or not there is any operable pumping, pumping, hydraulic power or the like whose output of the generator is smaller than the lower output limit. If not, the output calculation process ends, and the process returns to the process selection menu screen 20 of FIG. Flg is 1 in process 2004
In the case of, it is checked whether there is a constraint whose constraint has been relaxed and has reached the upper limit or the lower limit of the constraint. This is to determine whether there is a constraint that has been able to resolve the constraint violation by adding the constraint violation penalty to the objective function. If there are no constraints that can resolve the constraint violation,
Even if the processing is further continued, a good solution cannot be obtained, so that the constraint is relaxed and a solution satisfying the current constraint is stored as a result. If any of the constraint violations can be resolved by the check in the process 2009, the process 2011
Then, check with the constraint conditions set in the process 2000,
Extract those that violate constraints. In the process 2012, the item of the constraint in violation is stored. Furthermore, the constraint condition is relaxed for the violating constraint, and a penalty for the constraint violation is set. The penalty coefficient is a weight coefficient 45 for the constraint condition of FIG. 8 and a weight coefficient 21 of FIG.
5 or the product of the weighting factors for the equipment of the weighting factor 282 in FIG. The penalty changes the linear term of the objective function. That is, of the linear constraint expressions, an expression multiplied by a linear penalty coefficient of the generator output is added to the objective function. This changes the linear form of the objective function. In the next calculation, in order to find a possible solution, the constraint value for the constraint in which the constraint violation has occurred is relaxed. That is, when the upper limit is exceeded for the constraint value set in the process 2000, the lower limit of the constraint value is set as the upper limit,
The upper limit value is obtained by adding the violation amount to the original upper limit value.
If the constraint value set in the process 2000 exceeds the lower limit constraint, the upper limit value of the constraint value is set as the lower limit value, and the lower limit value is obtained by subtracting the violation amount from the original lower limit value. . As a result, an optimal solution can be obtained for the relaxed constraint. Next, in process 2013, flg is set to 1 so that it can be determined that the restriction has been relaxed. Again, in process 2002, an optimal solution to the relaxed constraint is determined.
【０１２４】処理２００５で最適解がないとき、最初に
発電機出力を計算したときに生じるケースであるが、処
理２０１０に進む。処理２０１０では、出力上下限制約
と需給バランスの制約を満たすように出力を決定し、こ
のときの出力を基準にして、違反している制約を最初に
同定するための処理である。コスト関数がない発電機に
対しては、図１１の初期値８６の値を用いて決定し、そ
の他の発電機は等λ法により発電機出力を決定する。If there is no optimal solution in the process 2005, this case occurs when the generator output is calculated first, but the process proceeds to the process 2010. In the process 2010, the output is determined so as to satisfy the output upper / lower limit constraint and the supply / demand balance constraint, and based on the output at this time, the violating constraint is first identified. The generator having no cost function is determined using the value of the initial value 86 in FIG. 11, and the other generators determine the generator output by the equal λ method.
【０１２５】処理２００３で最低出力より小さい出力に
配分されている発電機があるときは、最低出力未満の発
電機の１つを運転可能から停止に変更して、需給バラン
スを満たす供給力があるものがあるかチェックする。停
止により供給力不足になるならば、処理２００８に進
む。停止できる発電機があるならば、処理２０１５でそ
の発電機を停止にする。When there is a generator that is allocated to an output smaller than the minimum output in the process 2003, one of the generators having the minimum output is changed from operable to stopped to have a supply capacity that satisfies the supply and demand balance. Check for things. If the supply becomes insufficient due to the stop, the process proceeds to processing 2008. If there is a generator that can be stopped, the generator is stopped in step 2015.
【０１２６】以上により、制約条件を満たす最適な発電
機出力か緩和した制約条件での最適な電機出力を決定す
ることができる。As described above, it is possible to determine the optimum generator output that satisfies the constraint condition or the optimum electric machine output under the relaxed constraint condition.
【０１２７】（３）結果の表示 計画作成条件が全て設定されると、図１の処理２００８
により発電機等の出力を計算する。処理選択メニュー画
面２０において、出力計算が入力装置６を介して選択さ
れると、データ設定部１３は処理対象のデータベース８
の計算結果記憶領域１０４の全てのデータを読み込む。
次に、これらのデータを表示データ作成部１６に転送す
る。表示データ作成部では、計算結果を表示装置７に出
力する。出力画面の例を図２２〜図２６に示す。図２２
は目的関数，収束変数等を表示する画面３７０であり、
最終の計算結果である目的関数（燃料費），内点法の収
束状態の評価指標である双対ギャップ，主問題の制約違
反量，双対問題の制約違反量，満足できなかった制約違
反名称とその時刻，設備名称および違反量を含む。水力
の出力において、出力下限以下となったために、運転可
能から停止に変更した発電機名称とその時刻も含む。図
２３は各時刻の発電機の出力を表にした画面３５０であ
り、発電機名称と各時刻の出力を含む。図２４は１つの
発電機に対して計画対象時間（２４時間）の出力をグラ
フで表示した画面３６０であり、表示する発電機はいず
れでも選択できるものとする。図２５は等式制約の偏差
のグラフを表す画面３８０で等式制約を不等式制約とし
たとき下限又は上限側許容偏差のみを考慮すると、一方
の制限以上に変化する可能性があるので、等式制約から
の各時刻の偏差および最大偏差を含む。図２６は反復回
数毎の双対ギャップ，主問題の制約違反量，双対問題の
制約違反量および目的関数の値を含み、それらの変化を
グラフで表示する画面である。これらを表示することに
より、内点法で求めた解の妥当性を判断することができ
る。一連の表示を終了すると、処理選択メニュー画面に
戻る。(3) Display of Results When all the plan creation conditions are set, the processing 2008 in FIG.
To calculate the output of the generator and the like. When an output calculation is selected via the input device 6 on the process selection menu screen 20, the data setting unit 13 sets the database 8 to be processed.
All the data in the calculation result storage area 104 is read.
Next, these data are transferred to the display data creation unit 16. The display data creation unit outputs the calculation result to the display device 7. Examples of the output screen are shown in FIGS. FIG.
Is a screen 370 displaying an objective function, a convergence variable, and the like.
Objective function (fuel cost) that is the final calculation result, dual gap that is an evaluation index of the convergence state of the interior point method, constraint violation amount of the main problem, constraint violation amount of the dual problem, name of the constraint violation that was not satisfied and Includes time, equipment name, and violation amount. The output of the hydropower also includes the name of the generator whose operation was changed from operable to stopped because the output became lower than the output lower limit, and the time. FIG. 23 is a screen 350 showing the output of the generator at each time in a table, and includes the generator name and the output at each time. FIG. 24 is a screen 360 that displays the output of one generator in the planned time (24 hours) in a graph, and any of the generators to be displayed can be selected. FIG. 25 is a screen 380 showing a graph of the deviation of the equality constraint. When the equality constraint is set as the inequality constraint, if only the lower limit or the upper limit allowable deviation is considered, there is a possibility that the value may change more than one limit. Includes the deviation of each time from the constraint and the maximum deviation. FIG. 26 is a screen that displays the dual gap, the amount of constraint violation of the main problem, the amount of constraint violation of the dual problem, and the value of the objective function for each number of iterations, and displays their changes in a graph. By displaying these, the validity of the solution obtained by the interior point method can be determined. When a series of displays is completed, the process returns to the process selection menu screen.
【０１２８】以上では、１日２４時間の発電機出力を決
定する例で説明したが、期間および時間間隔はそれぞ
れ、１日，１時間とことなる場合も同じ処理で発電機出
力を決定することができる。In the above description, an example in which the generator output for 24 hours a day is determined has been described. Can be.
【０１２９】Ｃ．本実施例の効果 本実施例によれば、等式制約を不等式制約に変換するこ
とで、制約マトリックスの行方向に少なくとも１つの単
位行列を含ませることができ、解を更新するために制約
マトリックスの転置行列と元の行列との積のマトリック
スを小さくできるので、行列要素を記憶するメモリを少
なくできると同時に計算量を少なくできるので、処理を
高速化できる。C. According to the present embodiment, by converting the equality constraint to the inequality constraint, it is possible to include at least one unit matrix in the row direction of the constraint matrix. Since the matrix of the product of the transposed matrix and the original matrix can be reduced, the memory for storing the matrix elements can be reduced, and the amount of calculation can be reduced, so that the processing can be speeded up.
【０１３０】[0130]
【発明の効果】以上述べたように、本発明によれば、等
式制約を不等式制約に変換することで、制約マトリック
スの行方向に少なくとも１つの単位行列を含ませること
ができ、解を更新するために制約マトリックスの転置行
列と元の行列との積を計算する計算量を少なくできるの
で、制約条件を満たし目的関数を最小化する発電機等の
出力を高速化に計算することができる。As described above, according to the present invention, by converting an equality constraint into an inequality constraint, at least one unit matrix can be included in the row direction of the constraint matrix, and the solution is updated. Therefore, the amount of calculation for calculating the product of the transposed matrix of the constraint matrix and the original matrix can be reduced, so that the output of a generator or the like that satisfies the constraint conditions and minimizes the objective function can be calculated at high speed.
【図１】本発明の処理フローを示す図である。FIG. 1 is a diagram showing a processing flow of the present invention.
【図２】内点法の処理フローを示す図である。FIG. 2 is a diagram showing a processing flow of the interior point method.
【図３】出力決定装置のハードウェア構成を示す図であ
る。FIG. 3 is a diagram illustrating a hardware configuration of an output determination device.
【図４】機能構成を示す図である。FIG. 4 is a diagram showing a functional configuration.
【図５】データベースの構成を示す図である。FIG. 5 is a diagram showing a configuration of a database.
【図６】処理選択メニュー画面を示す図である。FIG. 6 is a diagram showing a process selection menu screen.
【図７】計画対象期間の入力，設定画面を示す図であ
る。FIG. 7 is a diagram showing a screen for inputting and setting a planning target period.
【図８】目的関数と考慮する制約条件を選択する画面を
示す図である。FIG. 8 is a diagram illustrating a screen for selecting a constraint condition to be considered as an objective function.
【図９】需要データの入力，設定画面を示す図である。FIG. 9 is a diagram showing a demand data input / setting screen.
【図１０】発電機等起動停止計画の入力，設定画面を示
す図である。FIG. 10 is a diagram showing a screen for inputting and setting a start / stop plan for a generator or the like.
【図１１】初期出力と供給力チェックの入力，設定画面
を示す図である。FIG. 11 is a diagram showing input and setting screens for initial output and supply power check.
【図１２】発電機等設備特性の入力，設定画面を示す図
である。FIG. 12 is a diagram showing a screen for inputting and setting equipment characteristics such as a generator.
【図１３】供給力積分制約の入力，設定画面を示す図で
ある。FIG. 13 is a diagram showing a screen for inputting and setting a supply force integration constraint.
【図１４】系統接続状態の入力，設定画面を示す図であ
る。FIG. 14 is a diagram showing a screen for inputting and setting a system connection state.
【図１５】系統設備の運用特性の入力，設定画面を示す
図である。FIG. 15 is a diagram showing a screen for inputting and setting operation characteristics of system equipment.
【図１６】収束判定条件の入力，設定画面を示す図であ
る。FIG. 16 is a diagram showing a screen for inputting and setting a convergence determination condition.
【図１７】内点法制御変数の入力，設定画面を示す図で
ある。FIG. 17 is a diagram showing a screen for inputting and setting interior point method control variables.
【図１８】内点法の初期値入力，設定画面を示す図であ
る。FIG. 18 is a diagram showing an initial value input / setting screen of the interior point method.
【図１９】等式制約の許容偏差の入力，設定画面を示す
図である。FIG. 19 is a diagram showing a screen for inputting and setting an allowable deviation of an equation constraint.
【図２０】制約マトリックス分割方法の入力，設定画面
を示す図である。FIG. 20 is a diagram showing an input / setting screen for a constraint matrix division method.
【図２１】簡易系統モデルを示す図である。FIG. 21 is a diagram showing a simplified system model.
【図２２】目的関数等の計算結果を示す図である。FIG. 22 is a diagram illustrating calculation results of an objective function and the like.
【図２３】発電機等の出力値を表す図である。FIG. 23 is a diagram illustrating output values of a generator and the like.
【図２４】発電機等の出力値をグラフで表す図である。FIG. 24 is a graph showing output values of a generator and the like.
【図２５】等式制約の等式からの偏差を示すグラフであ
る。FIG. 25 is a graph showing the deviation of the equation constraints from the equations.
【図２６】反復毎の収束指標の変化のグラフを示す図で
ある。FIG. 26 is a graph showing a change in a convergence index for each iteration.
１…出力決定装置、２…ＣＰＵ、３…主記憶装置、４…
入出力装置、５…外部記憶装置、６…入力装置、７…表
示装置、８…データベース、９…読取装置、１０…電力
システム、１２…コントロール部、１３…データ設定
部、１４…データ読込部、１５…出力計算部、１６…表
示データ作成部。DESCRIPTION OF SYMBOLS 1 ... Output determination device, 2 ... CPU, 3 ... Main storage device, 4 ...
Input / output device, 5 external storage device, 6 input device, 7 display device, 8 database, 9 reading device, 10 power system, 12 control unit, 13 data setting unit, 14 data reading unit , 15 ... Output calculation unit, 16 ... Display data creation unit.
───────────────────────────────────────────────────── フロントページの続き (72)発明者 田村 滋 茨城県日立市大みか町五丁目２番１号 株 式会社日立製作所大みか工場内 ────────────────────────────────────────────────── ─── Continuing on the front page (72) Inventor Shigeru Tamura 5-2-1 Omika-cho, Hitachi City, Ibaraki Prefecture Inside the Hitachi, Ltd. Omika Plant
Claims (17)
数の最適値を決定する問題を数理計画法を用いて解く上
で、出力変数が制約値と一致する等式制約を出力変数が
制約値から非負の定数Δ１を減算した値以上または／お
よび出力変数が制約値と非負の定数Δ２を加算した値以
下として定式化することを特徴とする最適な出力決定方
法。In solving a problem of deciding an optimum value of an output variable for minimizing or maximizing an objective function by using a mathematical programming, an equality constraint in which the output variable coincides with a constraint value is used for a constraint in the output variable. An optimal output determination method characterized by formulating a value equal to or greater than a value obtained by subtracting a non-negative constant Δ1 from a value and / or a value equal to or less than a value obtained by adding a constraint value and a non-negative constant Δ2 to an output variable.
て、電力系統の供給力を決定する上で、複数時間断面に
またがる積分制約を考慮することを特徴とする最適な出
力決定方法。2. An optimum output determination method according to claim 1, wherein an integration constraint over a plurality of time sections is considered when determining the supply power of the power system.
て、複数時間断面にまたがる積分制約がＬＮＧ消費制
約，ＩＰＰ利用率，揚水・水力発電所の池の水位制約で
あることを特徴とする最適な出力決定方法。3. An optimum output determination method according to claim 2, wherein the integral constraints over a plurality of time sections are LNG consumption constraints, IPP utilization rates, and water level constraints on ponds of pumped / hydroelectric power plants. Optimal output determination method.
て、変数の線形和で表される不等式制約式に下限制約値
を加算した式をスラック変数とし、該スラック変数の上
限制約値を該不等式の上限制約値と下限制約値の差と
し、下限制約値を０として定式化することを特徴とする
最適な出力決定方法。4. An optimum output determination method according to claim 1, wherein an expression obtained by adding a lower limit constraint value to an inequality constraint expression expressed by a linear sum of variables is set as a slack variable, and the upper limit constraint value of the slack variable is set as the slack variable. An optimal output determination method characterized by formulating a difference between an upper limit value and a lower limit value of an inequality and setting the lower limit value to 0.
て、該変数が生産量，流通量，電力、水道，ネットワー
クトラフィック，鉄鋼生産量または物流であり、制約値
が消費量，配達量または需要であることを特徴とする最
適な出力決定方法。5. The optimum output determination method according to claim 1, wherein said variables are production amount, distribution amount, electric power, water supply, network traffic, steel production amount or distribution, and said constraint value is consumption amount, delivery amount or distribution amount. An optimal output determination method characterized by demand.
て、定数Δ１または定数Δ２を入力する手段を有するこ
とを特徴とする最適な出力決定装置。6. An optimum output determining apparatus according to claim 1, further comprising means for inputting a constant Δ1 or a constant Δ2.
て、等式制約式の値と制約値との差を数値またはグラフ
として表示する手段を有することを特徴とする最適な出
力決定装置。7. An optimum output determining apparatus according to claim 1, further comprising means for displaying a difference between the value of the equality constraint equation and the constraint value as a numerical value or a graph.
て、該変数の線形和と該制約値の差がゼロ以上の定数Δ
３より大きいおよび／または該制約値と該変数の線形和
の差がゼロ以上の定数Δ４より大きいときに、その差を
数値またはグラフ，アラームメッセージとして表示する
手段を有することを特徴とする最適な出力決定装置。8. A method according to claim 1, wherein a difference between a linear sum of said variables and said constraint value is a constant Δ not less than zero.
When the difference between the constraint value and the linear sum of the variable and the constraint value is greater than zero or more than a constant Δ4, the means for displaying the difference as a numerical value, a graph, or an alarm message is provided. Output determination device.
て、変数の値または／および制約値を数値またはグラフ
として表示する手段を有することを特徴とする最適な出
力決定装置。9. An optimal output determining apparatus according to claim 1, further comprising means for displaying a value of a variable and / or a constraint value as a numerical value or a graph.
いて、決定した変数の値を用いて運転計画または運転制
御する手段を有することを特徴とする最適な設備の運転
計画制御システム。10. An optimal facility operation plan control system according to claim 1, further comprising means for performing an operation plan or operation control using the determined variable values.
現するためのプログラムを保持することを特徴とする機
械読み取り可能な記憶媒体。11. A machine-readable storage medium storing a program for realizing the optimum output determination method according to claim 1.
いて、適用する数理計画法が内点法であることを特徴と
する最適な出力決定方法。12. The optimum output determination method according to claim 1, wherein the applied mathematical programming is an interior point method.
いて、制約条件の変数にかかる係数をマトリックスで表
現したとき、該マトリックスをブロック化した小マトリ
ックスの列方向のブロックに少なくとも１つの単位行列
が含まれるように制約条件マトリックスを変換すること
を特徴とする最適な変数決定方法。13. The optimal variable determination method according to claim 1, wherein when a coefficient relating to a variable of the constraint condition is represented by a matrix, at least one unit matrix is added to a block in the column direction of the small matrix obtained by blocking the matrix. Optimum variable determination method characterized by transforming a constraint condition matrix so that is included.
いて、制約マトリックスを１つの時間以上毎に分割し、
制約条件を時間断面制約，差分制約および積分制約に分
割して、小ブロックのマトリックスにすることを特徴と
する最適な出力決定方法。14. The optimal output determination method according to claim 1, wherein the constraint matrix is divided every one time or more,
An optimal output determination method characterized by dividing a constraint condition into a time section constraint, a difference constraint, and an integral constraint to form a matrix of small blocks.
の最適値を決定する線形または二次計画問題を主双対内
点法を用いて解く上で、制約マトリックスをブロック化
し、主変数ｘ，ラグランジェ変数ｙと双対変数ｚの更新
にニュートン法を適用して、ニュートン方向の変化ベク
トルを計算するとき、主変数ｘの変化ベクトルを最初に
計算し、該主変数ｘの値を代入することにより他の変数
を決定することを特徴とする最適な出力決定方法。15. A linear or quadratic programming problem for determining an optimal value of a variable for minimizing or maximizing an objective function is solved using a principal dual interior point method. When calculating the change vector in the Newton direction by applying the Newton method to the update of the Lagrange variable y and the dual variable z, first calculate the change vector of the main variable x and substitute the value of the main variable x. An optimal output determination method characterized by determining other variables by
て、解候補に対して、満たせない制約の違反量を計算す
る手段と制約値の上限制約値を超過した制約に対して
は、下限制約値を上限制約値に置き換え、該上限制約値
に違反量を加算した値を上限制約値とする処理と制約値
の下限を下回った制約に対しては、該下限値から違反量
を減算した値を下限制約値をとし上限制約値を下限制約
値に置き換える手段とを有することを特徴とする最適な
出力決定装置。16. An optimal output determination method according to claim 1, wherein means for calculating a violation amount of a constraint that cannot be satisfied for a solution candidate and a lower limit constraint for a constraint that exceeds an upper limit constraint value. The value is replaced with the upper limit value, and the value obtained by adding the violation amount to the upper limit value is used as the upper limit value. And a means for replacing the upper limit value with the lower limit value.
て、元の制約値に対して違反した量に応じたペナルティ
を目的関数に加算する手段を有することを特徴とする最
適な出力決定装置。17. An optimal output determining apparatus according to claim 16, further comprising means for adding a penalty according to an amount violating an original constraint value to an objective function.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP5643098A JPH11259450A (en) | 1998-03-09 | 1998-03-09 | Optimal output deciding method and device therefor |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP5643098A JPH11259450A (en) | 1998-03-09 | 1998-03-09 | Optimal output deciding method and device therefor |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH11259450A true JPH11259450A (en) | 1999-09-24 |
Family
ID=13026885
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP5643098A Pending JPH11259450A (en) | 1998-03-09 | 1998-03-09 | Optimal output deciding method and device therefor |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH11259450A (en) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001184334A (en) * | 1999-12-24 | 2001-07-06 | Nec Corp | Method and device for finding solution of linear scheduling question |
| JP2002373002A (en) * | 2001-06-14 | 2002-12-26 | Toshiba Corp | Device and method for process simulater application and control |
| JP2004265381A (en) * | 2002-09-11 | 2004-09-24 | Fisher Rosemount Syst Inc | Feasibility handling of constraint and limit in process control system optimizer |
| JP2006034004A (en) * | 2004-07-16 | 2006-02-02 | Chugoku Electric Power Co Inc:The | Apparatus and method for determining quantity of generator output and program |
| US7330767B2 (en) | 2002-09-11 | 2008-02-12 | Fisher-Rosemount Systems, Inc. | Configuration and viewing display for an integrated model predictive control and optimizer function block |
| JP2010079323A (en) * | 2008-09-23 | 2010-04-08 | Mitsubishi Electric Corp | Secondary plan problem calculation device, program for secondary plan problem calculation device, periodic current plan calculation device, generator output value calculation device, and portfolio optimization device |
| JP2010226879A (en) * | 2009-03-24 | 2010-10-07 | Yamatake Corp | Method and device for generating operation plan of energy supply system |
| JP2018191482A (en) * | 2017-05-11 | 2018-11-29 | 三菱電機株式会社 | Optimum tidal current calculation device, optimum tidal current calculation method, and optimum tidal current calculation program |
| CN113498523A (en) * | 2019-03-06 | 2021-10-12 | 三菱电机株式会社 | Apparatus and method for controlling operation of machine object and storage medium |
| JP2023100054A (en) * | 2022-01-05 | 2023-07-18 | 東芝エネルギーシステムズ株式会社 | Hydroelectric power station operation schedule generation device and hydroelectric power station operation schedule generation method |
-
1998
- 1998-03-09 JP JP5643098A patent/JPH11259450A/en active Pending
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001184334A (en) * | 1999-12-24 | 2001-07-06 | Nec Corp | Method and device for finding solution of linear scheduling question |
| JP2002373002A (en) * | 2001-06-14 | 2002-12-26 | Toshiba Corp | Device and method for process simulater application and control |
| JP2004265381A (en) * | 2002-09-11 | 2004-09-24 | Fisher Rosemount Syst Inc | Feasibility handling of constraint and limit in process control system optimizer |
| US7330767B2 (en) | 2002-09-11 | 2008-02-12 | Fisher-Rosemount Systems, Inc. | Configuration and viewing display for an integrated model predictive control and optimizer function block |
| US7376472B2 (en) | 2002-09-11 | 2008-05-20 | Fisher-Rosemount Systems, Inc. | Integrated model predictive control and optimization within a process control system |
| JP2006034004A (en) * | 2004-07-16 | 2006-02-02 | Chugoku Electric Power Co Inc:The | Apparatus and method for determining quantity of generator output and program |
| JP2010079323A (en) * | 2008-09-23 | 2010-04-08 | Mitsubishi Electric Corp | Secondary plan problem calculation device, program for secondary plan problem calculation device, periodic current plan calculation device, generator output value calculation device, and portfolio optimization device |
| JP2010226879A (en) * | 2009-03-24 | 2010-10-07 | Yamatake Corp | Method and device for generating operation plan of energy supply system |
| JP2018191482A (en) * | 2017-05-11 | 2018-11-29 | 三菱電機株式会社 | Optimum tidal current calculation device, optimum tidal current calculation method, and optimum tidal current calculation program |
| CN113498523A (en) * | 2019-03-06 | 2021-10-12 | 三菱电机株式会社 | Apparatus and method for controlling operation of machine object and storage medium |
| JP2022513416A (en) * | 2019-03-06 | 2022-02-07 | 三菱電機株式会社 | Devices and methods for controlling the operation of machines, as well as storage media |
| CN113498523B (en) * | 2019-03-06 | 2024-04-30 | 三菱电机株式会社 | Apparatus and method for controlling operation of machine object and storage medium |
| JP2023100054A (en) * | 2022-01-05 | 2023-07-18 | 東芝エネルギーシステムズ株式会社 | Hydroelectric power station operation schedule generation device and hydroelectric power station operation schedule generation method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kong et al. | An overview on formulations and optimization methods for the unit-based short-term hydro scheduling problem | |
| Nguyen et al. | A bidding strategy for virtual power plants with the intraday demand response exchange market using the stochastic programming | |
| Jenkins et al. | Enhanced decision support for a changing electricity landscape: the GenX configurable electricity resource capacity expansion model | |
| Bustos et al. | Energy storage and transmission expansion planning: substitutes or complements? | |
| Catalao et al. | Scheduling of head-sensitive cascaded hydro systems: A nonlinear approach | |
| CN112016747B (en) | Optimization method suitable for source-load-storage flexible resource overall planning and operation | |
| JP4658089B2 (en) | Generator operation plan formulation device, operation plan formulation method, and program for executing the device and method | |
| JP2001211548A (en) | Method for power-generation planning and device | |
| CN114498639A (en) | Day-ahead scheduling method for multi-microgrid combined mutual aid considering demand response | |
| JP2001258157A (en) | Device and method of developing operation plan for generator and recording medium storing operation plan program for generator | |
| JPH11259450A (en) | Optimal output deciding method and device therefor | |
| Sadek et al. | Adaptive robust energy management for isolated microgrids considering reactive power capabilities of distributed energy resources and reactive power costs | |
| Bhuvanesh et al. | Application of differential evolution algorithm and its variants for solving energy storage technologies integrated generation expansion planning | |
| Nikpour et al. | Stochastic-risk based approach for Microgrid participation in joint active, reactive, and ancillary services markets considering demand response | |
| Nikpour et al. | Hybrid stochastic risk-based approach for a microgrid participating in coupled active and reactive power market | |
| CN110247392A (en) | More standby resources robust Optimal methods of meter and wind-powered electricity generation marginal capacity and Demand Side Response | |
| Bjorndal et al. | Flow-based market coupling in the european electricity market–a comparison of efficiency and feasibility | |
| Nalini et al. | Towards prosumer flexibility markets: A photovoltaic and battery storage model | |
| JP2004274851A (en) | Energy accommodation system | |
| Gacitúa et al. | The role of fast-acting energy storage for contingency grid support in the transmission planning | |
| Gebrekiros | Analysis of integrated balancing markets in northern Europe under different market design options | |
| Afraz et al. | Active distribution network planning considering shared demand management | |
| Gea-Bermúdez et al. | North sea offshore grid-effects of integration towards 2050 | |
| Laws et al. | Valuing distributed energy resources for non-wires alternatives | |
| Sadighmanesh et al. | Probabilistic dispatch in hybrid-microgrid system with considering energy arbitrage |