JP7068557B2 - Power controller and power distribution plan setting method - Google Patents
Power controller and power distribution plan setting method Download PDFInfo
- Publication number
- JP7068557B2 JP7068557B2 JP2021534845A JP2021534845A JP7068557B2 JP 7068557 B2 JP7068557 B2 JP 7068557B2 JP 2021534845 A JP2021534845 A JP 2021534845A JP 2021534845 A JP2021534845 A JP 2021534845A JP 7068557 B2 JP7068557 B2 JP 7068557B2
- Authority
- JP
- Japan
- Prior art keywords
- distribution
- power
- power distribution
- plan
- unit
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims description 37
- 239000013598 vector Substances 0.000 claims description 67
- 238000006243 chemical reaction Methods 0.000 claims description 34
- 238000004364 calculation method Methods 0.000 claims description 22
- 238000011156 evaluation Methods 0.000 claims description 11
- 238000013507 mapping Methods 0.000 claims description 7
- 230000001131 transforming effect Effects 0.000 claims 1
- 230000006870 function Effects 0.000 description 38
- 238000010248 power generation Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000036962 time dependent Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Power Engineering (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Public Health (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Description
本願は、電力制御装置および電力配分計画設定方法に関するものである。 The present application relates to a power control device and a power distribution plan setting method.
再生可能エネルギー等の普及に伴い、配電側にも電源を持つようになっている。そこで、再生可能エネルギーの発電量も含めて電力需要を予測し、コストが最小になるように、配電系統の能力に応じて需要予測値と供給計画値の差分を分配する電力系統の制御システムが開示されている(例えば、特許文献1参照。)。しかし、再生可能エネルギーの発電可能量、または電力需要の予測値は、必ずしも実際の発電可能量、または電力需要と一致せず、通常は誤差を持つ。したがって、事前に計算した配電量に基づき配電系統を制御したとしても、実際の配電系統に関して最適な制御とはならない。 With the spread of renewable energy, etc., the distribution side also has a power source. Therefore, a power system control system that predicts power demand including the amount of power generated by renewable energy and distributes the difference between the demand forecast value and the supply plan value according to the capacity of the distribution system so that the cost is minimized. It is disclosed (see, for example, Patent Document 1). However, the power generation capacity of renewable energy or the predicted value of power demand does not always match the actual power generation capacity or power demand, and usually has an error. Therefore, even if the distribution system is controlled based on the distribution amount calculated in advance, the control is not optimal for the actual distribution system.
この問題に対し、あらかじめ発電可能量、または電力需要に関する変動幅を確率分布として表示し、電力供給コストをより現実的に見積もる手法が開示されている(例えば、特許文献2参照。)。 To solve this problem, a method of displaying the amount of power generation or the fluctuation range related to the power demand as a probability distribution in advance and estimating the power supply cost more realistically is disclosed (see, for example, Patent Document 2).
しかしながら、見積もりが現実的であったとしても、確率分布の状態では、実際の配電量の制御に活かすことは困難であり、実態に応じた見積もりに基づく制御を行うことは困難であった。 However, even if the estimation is realistic, it is difficult to utilize it for the control of the actual power distribution amount in the state of the probability distribution, and it is difficult to perform the control based on the estimation according to the actual situation.
本願は、上記のような課題を解決するための技術を開示するものであり、実態に応じた見積もりに基づく制御を可能とする電力制御装置および電力配分計画設定方法を得ることを目的とする。 The present application discloses a technique for solving the above-mentioned problems, and an object thereof is to obtain a power control device and a power distribution plan setting method capable of control based on an estimate according to an actual situation.
本願に開示される電力制御装置は、蓄電池を含む配電系統をリーマン多様体とみなし、前記配電系統において予測される電力供給量と電力需要から、前記蓄電池の残量を並べたベクトルと前記配電系統での制御可能な量を並べたベクトルとを変数とする前記配電系統での電力損失を評価する関数と、前記蓄電池の目標充電量からのずれを示す関数の和を示す評価関数を最小化する電力配分計画を生成する配分計画部、前記配電系統における実際の配分状態を取得する配分状態取得部、前記電力配分計画に基づいて前記配電系統での電力配分を制御する電力配分制御部、および前記電力配分計画が、前記配分状態に対してずれを有する際、フィードバック制御によって、前記ずれを有する電力配分計画と前記配分状態との距離が最小となるよう、接ベクトルと点を相互に滑らかに変換する写像を用いて前記ずれを有する電力配分計画の調整値を生成する配分調整部、を備えたことを特徴とする。 The power control device disclosed in the present application regards a distribution system including a storage battery as a Lehman variety, and has a vector in which the remaining amount of the storage battery is arranged from the power supply amount and power demand predicted in the distribution system and the distribution system. Minimize the evaluation function that indicates the sum of the function that evaluates the power loss in the distribution system and the function that indicates the deviation from the target charge amount of the storage battery, with the vector that arranges the controllable amounts in A distribution planning unit that generates a power distribution plan, a distribution state acquisition unit that acquires an actual distribution state in the distribution system, a power distribution control unit that controls power distribution in the distribution system based on the power distribution plan, and the above. When the power distribution plan has a deviation from the distribution state , the tangent vector and the point are smoothly converted from each other so that the distance between the power distribution plan having the deviation and the distribution state is minimized by feedback control. It is characterized by comprising a distribution adjustment unit for generating an adjustment value of a power distribution plan having the deviation by using a mapping.
本願に開示される電力配分計画設定方法は、蓄電池を含む配電系統をリーマン多様体とみなし、前記配電系統において、電力供給量と電力需要を予測する予測ステップ、前記予測された電力供給量と電力需要から、前記蓄電池の残量を並べたベクトルと前記配電系統での制御可能な量を並べたベクトルとを変数とする前記配電系統での電力損失を評価する関数と、前記蓄電池の目標充電量からのずれを示す関数の和を示す評価関数を最小化する電力配分計画を生成する電力配分計画生成ステップ、前記電力配分計画に基づいて前記配電系統での電力配分を制御する電力配分制御ステップ、前記配電系統における実際の配分状態を取得する配分状態取得ステップ、および前記電力配分計画が、前記配分状態に対してずれを有する際、フィードバック制御によって、前記ずれを有する電力配分計画と前記配分状態との距離が最小となるよう、接ベクトルと点を相互に滑らかに変換する写像を用いて前記ずれを有する電力配分計画の調整値を生成する調整値生成ステップ、を含むことを特徴とする。
The power distribution plan setting method disclosed in the present application regards a distribution system including a storage battery as a Lehman variety, and in the distribution system, a prediction step for predicting power supply amount and power demand, and the predicted power supply amount and power. From the demand, a function for evaluating the power loss in the distribution system with a vector in which the remaining amount of the storage battery is arranged and a vector in which the controllable amount in the distribution system is arranged are used as variables, and a target charge amount of the storage battery. A power distribution plan generation step that generates a power distribution plan that minimizes an evaluation function that indicates the sum of functions indicating deviations from, and a power distribution control step that controls power distribution in the distribution system based on the power distribution plan. When the distribution state acquisition step for acquiring the actual distribution state in the distribution system and the power distribution plan have a deviation from the distribution state, the power distribution plan having the deviation and the distribution state are subjected to feedback control. It is characterized by including an adjustment value generation step of generating an adjustment value of a power distribution plan having the deviation using a mapping that smoothly transforms a tangent vector and a point into each other so that the distance between the two is minimized.
本願に開示される電力制御装置あるいは電力配分計画設定方法によれば、予測と実態が一致しなかった場合でも、フィードバックが可能で、実態に応じた見積もりに基づく制御を可能とする。 According to the power control device or the power distribution plan setting method disclosed in the present application, even if the prediction and the actual situation do not match, feedback is possible and control based on the estimation according to the actual situation is possible.
実施の形態1.
図1~図7は、実施の形態1にかかる電力制御装置の構成と動作、および電力配分計画設定方法について説明するためのものであり、図1は電力制御装置の構成を説明するためのブロック図、図2は電力制御装置の対数変換部の動作、あるいは電力配分計画設定方法における第1の接ベクトルを算出する手順を示す模式図、図3は配分計画部の動作、あるいは電力配分計画設定方法における電力配分計画生成ステップを示すフローチャート、図4は配分調整部を構成する対数変換部の動作、あるいは電力配分計画設定方法における対数変換ステップを示すフローチャートである。
1 to 7 are for explaining the configuration and operation of the power control device according to the first embodiment, and the power distribution plan setting method, and FIG. 1 is a block for explaining the configuration of the power control device. FIG. 2 is a schematic diagram showing the operation of the logarithmic conversion unit of the power control device or the procedure for calculating the first tangent vector in the power distribution plan setting method, and FIG. 3 is the operation of the distribution planning unit or the power distribution plan setting. A flowchart showing the power distribution plan generation step in the method, and FIG. 4 is a flowchart showing the operation of the logarithmic conversion unit constituting the distribution adjustment unit or the logarithmic conversion step in the power distribution plan setting method.
また、図5は電力制御装置の全体動作、あるいは電力配分計画設定方法の全体動作を示すフローチャート、図6は電力制御装置、あるいは電力配分計画設定方法の効果を説明するための、状況変化が発生したときの、計画上の充電量の経時変化と、本実施の形態によるフィードバックを行った際の充電量の経時変化を示すグラフ形式の図である。さらに、図7は電力制御装置の演算処理を実行する部分の構成例を示すブロック図である。 Further, FIG. 5 is a flowchart showing the overall operation of the power control device or the overall operation of the power distribution plan setting method, and FIG. 6 is a situation change for explaining the effect of the power control device or the power distribution plan setting method. It is a graph-type diagram which shows the time-dependent change of the charge amount in the planned time, and the time-dependent change of the charge amount at the time of giving feedback by this embodiment. Further, FIG. 7 is a block diagram showing a configuration example of a portion that executes arithmetic processing of the power control device.
実施の形態1にかかる電力制御装置1は、図1に示すように、電力供給量を予測する供給予測部13、電力需要を予測する需要予測部14、供給予測部13と需要予測部14の予測に基づき、電力配分計画を生成する配分計画部15を備えている。そして、配電系統800の配分状態を取得する配分状態取得部16と、取得した配分状態と計算された電力配分計画とのずれを評価して、配分を調整する配分調整部12と、調整した配分により、配電系統800の電力配分を制御する電力配分制御部11とを備える。なお、配電系統800には、図示しない一つ以上の電力源と、負荷、および充放電を行う蓄電池が含まれているものとする。
As shown in FIG. 1, the
供給予測部13は、配電系統800に含まれる太陽光発電、風力発電、燃料電池等の電力源による電力供給量を予測し、予測供給量Qpを配分計画部15に出力する。また、需要予測部14は、配電系統800に含まれる負荷等の電力需要を予測し、予測需要量Dpを配分計画部15に出力する。また、配分状態取得部16は、配電系統800における実際の電力供給と電力需要の状態ハットztを取得し、配分調整部12に出力する。
The
<配分計画部>
配分計画部15は、供給予測部13から受信した予測供給量Qpと需要予測部14から受信した予測需要量Dpに基づき、最適な電力配分計画ut、ztとフィードバック係数Ht,νν、Ht,νzを算出する。電力配分計画ut、ztの算出には、蓄電池の残量を並べたベクトルxt、配電系統800内での再生可能エネルギーでの発電量、融通される電力量等を並べたベクトルyt、および蓄電池の充放電電力を含む電力制御装置1で制御可能な量を並べたベクトルである入力utを用いる。そして、ベクトルxt、ytをまとめたベクトルztを用い、目標充電量からのずれを表す関数と、電力ロス、コストロス等の経済的損失を表す関数の和として定義される、後述する評価関数Eを最小化する電力配分計画ut、ztを求める。そして、フィードバック係数Ht,νν、Ht,νzとともに、配分調整部12に出力する。<Distribution Planning Department>
The
<配分調整部>
配分調整部12は、第1の接ベクトルdzを算出する対数変換部121と、第2の接ベクトルduを算出する制御入力計算部122と、フィードバック入力uを算出して、電力配分制御部11に出力する指数変換部123とを備えている。そして、配電系統800がなす多様体において、電力配分計画から配電系統800の状態に向けた接ベクトルdzを用いて接空間におけるフィードバック入力duを算出し、それを実際のフィードバック入力uに変換して、電力配分制御部11に出力する。<Distribution adjustment department>
The distribution adjustment unit 12 calculates the
<対数変換部>
対数変換部121は、配分計画部15から出力された現在の電力配分計画ut、ztと、配分状態取得部16から出力された現在の状態ハットztから、対数写像の値である第1の接ベクトルdzを算出する。より具体的には、配電系統800がなす多様体において、現在の電力配分計画ut、ztから、現在の配電系統800の現在の状態であるハットztに向けた、第1の接ベクトルdzを算出し、制御入力計算部122に出力する。<Logarithmic conversion unit>
The
<制御入力計算部>
制御入力計算部122は、対数変換部121から出力された第1の接ベクトルdzと、配分計画部15から出力されたフィードバック係数Ht,νν、Ht,νzから、接空間におけるフィードバック入力である第2の接ベクトルduを算出する。そして、算出した第2の接ベクトルduを指数変換部123に出力する。<Control input calculation unit>
The control
<指数変換部>
指数変換部123は、第2の接ベクトルduを指数変換して、実際のフィードバック入力uを算出し、電力配分制御部11に出力する。以下、詳細に説明する。<Exponential conversion unit>
The
本実施の形態1にかかる電力制御装置1あるいは電力配分計画設定方法における電力配分問題は、評価関数E(式1)を最小化する入力ut:tfを求める問題として定式化される。The power distribution problem in the
ここで、xtは1つまたは複数の蓄電池の残量を並べたベクトルである。また、ytは太陽光、風力等の発電量、配電系統内で融通される電力等を並べたベクトルであってよい。入力utは蓄電池の充放電電力、発電機、あるいは上位電力系統からの受電量、買電量等、電力制御装置によって制御可能な量を並べたベクトルであってよい。関数Lは蓄電池の充放電損失、変換器の変換損失、電力あたりの買電価格等を表す関数であってよい。関数Sは蓄電池の目標充電量からのずれを表す関数であってよい。つまり、配電系統800内での電力損失あるいは買電価格等を表す関数Lの積分と、蓄電池の目標充電量からのずれを示す関数Sとの和を示す評価関数Eを最小化する入力u0:tfを求めることになる。Here, x t is a vector in which the remaining amount of one or a plurality of storage batteries is arranged. Further, yt may be a vector in which the amount of power generation such as solar power and wind power, the electric power interchanged in the distribution system, and the like are arranged. The input ut may be a vector in which the amounts that can be controlled by the power control device, such as the charge / discharge power of the storage battery, the generator, the amount of power received from the upper power system, and the amount of power purchased, are arranged. The function L may be a function representing the charge / discharge loss of the storage battery, the conversion loss of the converter, the purchase price per electric power, and the like. The function S may be a function representing a deviation from the target charge amount of the storage battery. That is, the input u 0 that minimizes the evaluation function E indicating the sum of the integral of the function L representing the power loss or the purchase price in the
このとき、任意のxに対し、h(x,y)=0を満たすyは少なくとも局所的に一意に定まる。ここでxt、ytをまとめたベクトルをztとして、ztとutからzt+1を与える写像をFとおく。評価関数Eの、時刻tからtfまでの部分和をJt(zt,ut:tf)とおくと、部分和は式(2)を満たす。
Jt(zt,ut:tf)=L(zt,ut)+Jt+1(F(zt,ut),ut+1:tf) ・・・(2)
ただし、Jtf(ztf)=S(ztf)である。そして、動的計画法を適用することで、式(3)が成り立つ。At this time, y satisfying h (x, y) = 0 is uniquely determined at least locally for any x. Here, let z t be a vector that summarizes x t and y t , and let F be a map that gives z t + 1 from z t and ut. If the partial sum of the evaluation function E from the time t to t f is set to J t (z t , ut : tf ), the partial sum satisfies the equation (2).
J t (z t , u t: tf ) = L (z t , u t ) + J t + 1 (F (z t , u t ), u t + 1: tf ) ・ ・ ・ (2)
However, J tf (z tf ) = S (z tf ). Then, by applying the dynamic programming method, the equation (3) is established.
関数Vt(zt)は、Ht(z,ν)のνに関する最小値であるから、式(5)を式(4)に代入して、式(6)を得る。
Since the function V t (z t ) is the minimum value of H t (z, ν) with respect to ν, the equation (5) is substituted into the equation (4) to obtain the equation (6).
∇i´V=∇j´Ht-(∇iHt)(∇ijHt)(∇ij´Ht)
∇i´j´V=∇i´j´Ht-(∇i´jHt)(∇ijHt)(∇ij´Ht) ・・・(7)
一方、式(8)の定義から、式群(9)および式群(10)を得る。
∇ i´ V = ∇ j´ H t- (∇ i H t ) (∇ ij H t ) (∇ ij´ H t )
∇ i´j´ V = ∇ i´j´ H t- (∇ i´j H t ) (∇ ij H t ) (∇ ij´ H t ) ・ ・ ・ (7)
On the other hand, from the definition of the formula (8), the formula group (9) and the formula group (10) are obtained.
さて、本願の実施の形態にかかる電力制御装置1あるいは電力配分計画設定方法における特徴は、ある入力u0:tfに対し、更新後の入力ハットu0:tfを式群(11)によって求める点にある。The feature of the
まず、等式制約h(z)=0のなす多様体において、zの測地線方程式は、式(12)で与えられる。 First, in the manifold with the equation constraint h (z) = 0, the geodesic equation of z is given by equation (12).
ここで、h(z)=0のなすリーマン多様体上の非負かつ非退化な関数D(z,z´)で、少なくともz´の近傍でzに関するヘッシアン(Hessian)が正定値となるようなものを考える。本実施形態において、D(z,z´)はユークリッド距離の2乗としてもよい。このD(z,z´)を用いて対数写像Logztハットztを効率的に求めることができる。Here, in the non-negative and non-degenerate function D (z, z ′ ) on the Riemannian manifold with h (z) = 0, the Hessian with respect to z becomes a positive-definite value at least in the vicinity of z ′ . Think about things. In this embodiment, D (z, z ′ ) may be the square of the Euclidean distance. Using this D (z, z ′ ), the logarithmic map Log zt hat zt can be efficiently obtained.
今、D(z,ハットzt)をzの関数と考えると、D(z,ハットzt)は、ハットztにおいて最小値0をとる凸関数である。従って、点ztにおけるテイラー(Taylor)展開(式(13))に基づくD(z,ハットzt)の最小化方向dzi=-(∇ijD)∇jDは、Expzt[dz]がハットztとなるようなdzを近似する。Now, considering D (z, hat z t ) as a function of z, D (z, hat z t ) is a convex function having a minimum value of 0 in hat z t . Therefore, the minimization direction dz i =-(∇ ij D) ∇ j D of D (z, hat z t ) based on the Taylor expansion at the point z t (Equation (13)) is Exp zz [dz]. Approximate dz such that is hat z t .
点ztから距離関数D(z,ハットzt)を最小化する方向へのベクトルがdzであり、dz方向への測地線の先がrである。rからD(z,ハットzt)を最小化する方向へのベクトルがdチェックz(r)であり、これを測地線にそって点ztに引き戻したベクトルがdチェックz(zt)である。新たな最小化方向はdz+dチェックz(zt)であり、その方向への測地線はrよりもハットztに近づくことがわかる。The vector in the direction of minimizing the distance function D (z, hat z t ) from the point z t is dz, and the tip of the geodesic line in the dz direction is r. The vector in the direction of minimizing D (z, hat z t ) from r is d check z (r), and the vector pulled back to the point z t along the geodesic line is d check z (z t ). Is. It can be seen that the new minimization direction is dz + d check z (z t ), and the geodesic in that direction is closer to the hat z t than to r.
等式制約h(z)=0、およびl(z,u)=0のなすリーマン多様体の、点utにおける指数写像Exputは、zとuを新たにwとおくことで、式(12)と同様の測地線方程式が得られる。得られた測地線方程式のうち、ドットzに関する部分を、対数写像Logztを求めたときの測地線とすることで得られる。The exponential map Expon at the point ut of the Riemannian manifold with the equation constraints h (z) = 0 and l (z, u) = 0 is expressed by setting z and u as w. A geodesic equation similar to 12) can be obtained. It is obtained by setting the part related to the dot z in the obtained geodesic equation as the geodesic line when the logarithmic map Log zt is obtained.
つぎに、図3のフローチャートを用いて、配分計画部15における電力配分計画u0:tf、z0:tfの生成手順について説明する。はじめに、初期状態x0、y0、および初期入力u0:tfを定める(ステップS201)。つづいて、zt+1=F(zt,ut)に従い、ztとutからzt+1を計算する(ステップS202)。なお、ステップS202は、t=0からtf-1まで繰り返される。Next, the procedure for generating the power distribution plans u 0: tf and z 0: tf in the
次に、式(7)、(9)、(10)を用い、Vt+1の微分からHtとVtの微分を求める(ステップS203)。ステップS203は、t=tf-1からt=0まで繰り返される。そして、式(15)に示す勾配のノルムを算出し、これが所定のしきい値αより小さいかどうかを判定する(ステップS204)。Next, using equations (7), (9), and (10), the derivative of H t and V t is obtained from the derivative of V t + 1 (step S203). Step S203 is repeated from t = t f -1 to t = 0. Then, the norm of the gradient shown in the equation (15) is calculated, and it is determined whether or not this is smaller than the predetermined threshold value α (step S204).
つぎに、図4のフローチャートを用いて、対数変換部121での第1の接ベクトルdzの計算手順(対数変換ステップ)を表す。配分計画部15が出力した時刻tでの電力配分計画ut、ztと配分状態取得部16から取得した現在の状態ハットztから、点(ut,zt)におけるD(z,ハットzt)の最小化方向dzを計算する(ステップS301)。つづいて、点(ut,zt)から、dz方向への測地線を計算し、r=Expzt[dz]を求める(ステップS302)。Next, using the flowchart of FIG. 4, the calculation procedure (logarithmic transformation step) of the first tangent vector dz in the
そして、D(r,ハットzt)のノルムが、所定のしきい値βより小さいかどうかを判定する(ステップS303)。D(r,ハットzt)のノルムが、しきい値βよりも小さいならば(「Yes」)、dzを出力する(ステップS307)。一方、しきい値β以上であれば(「No」)、点rにおける最小化方向dチェックzを計算し(ステップS304)、平行移動方程式(式(14))にしたがい、dチェックzを点ztまで平行移動する(ステップS305)。続いて、ステップS305によって得られた方向dチェックzをdzに足し合わせ、新たなdzとする(ステップS306)。Then, it is determined whether or not the norm of D (r, hat z t ) is smaller than the predetermined threshold value β (step S303). If the norm of D (r, hat z t ) is smaller than the threshold value β (“Yes”), dz is output (step S307). On the other hand, if it is equal to or higher than the threshold value β (“No”), the minimization direction d check z at the point r is calculated (step S304), and the d check z is set as a point according to the translation equation (Equation (14)). Translate to zt (step S305). Subsequently, the direction d-check z obtained in step S305 is added to dz to obtain a new dz (step S306).
制御入力計算部122は、配分計画部15が計算したフィードバック係数Ht,νν、Ht,νzから、式(16)により、接空間におけるフィードバック入力であるduを計算する。
du=-[Ht,νν
-1Ht,νz]udz ・・・(16)
ここで、dzは、第1の接ベクトルに相当し、duは、第2の接ベクトルに相当する。 The control
du =-[H t, ν ν -1 H t, ν z ] u dz ・ ・ ・ (16)
Here, dz corresponds to the first tangent vector, and du corresponds to the second tangent vector.
指数変換部123は、u=Exputduにより、接空間におけるフィードバック入力duから、実際のフィードバック入力uを計算し、電力配分制御部11は、計算されたフィードバック入力uに基づき、配電系統800における電力配分を制御する。The
上述した動作の全体、つまり電力配分計画設定方法について図5のフローチャートに基づき説明する。はじめに、供給予測部13による予測供給量Qpの生成、需要予測部14による予測需要量Dpの生成が行われ、配分計画部15に出力される(ステップS100)。
The entire operation described above, that is, the power distribution plan setting method will be described with reference to the flowchart of FIG. First, the
配分計画部15は、予測供給量Qpと予測需要量Dpをもとに、上述したベクトルxt、yt、utを用いて、蓄電池の目標充電量からのずれを示す関数と、配電系統800内での損失を示す関数の和として定義した評価関数Eを最小化する電力配分計画を求める。そして、写像Fによりu0:tf、z0:tfを、式(7)~(10)により、Ht、Vtの微分を計算し、条件(||∇H||<α)を満たすことを前提として、それぞれ電力配分計画、フィードバック係数として出力する(ステップS110(ステップS201~S206))。Based on the predicted supply amount Qp and the predicted demand amount Dp, the
このとき、電力配分計画は、配分調整部12の対数変換部121に出力され、フィードバック係数は、制御入力計算部122に出力される。なお、条件(||∇H||<α)を満たさない場合は、式群(11)を用いて、新たな入力と状態を算出し、条件を満たすまで再計算を行う。
At this time, the power distribution plan is output to the
電力配分制御部11は、少なくとも電力配分計画が生成された直後は、生成された電力配分計画u0:tf、z0:tfに基づいて、配電系統800内の電力配分制御を実行する(ステップS120)。The power
一方、配分状態取得部16は、配電系統800から現実の配分状態ハットztを取得し(ステップS130)、配分調整部12の対数変換部121に出力する。対数変換部121は、配分状態ハットztと電力配分計画u0:tf、z0:tfを受信すると、計画と実態(電力配分計画u0:tf、z0:tfと配分状態ハットzt)が一致するか否かを判定(ステップS140)し、一致している場合(Yes)は、ステップS120へ戻る。一方、一致しなかった場合(No)は、電力配分計画の修正値であるフィードバック入力uを演算するため、現在の状態ハットztから、対数写像の値である第1の接ベクトルdzを算出し(ステップS150(ステップS301~S307))、制御入力計算部122に出力する。On the other hand, the distribution
制御入力計算部122は、第1の接ベクトルdzと、配分計画部15からのフィードバック係数Ht,νν、Ht,νzを用い、接空間におけるフィードバック入力である第2の接ベクトルduを算出し(ステップS160)、指数変換部123に出力する。指数変換部123は、第2の接ベクトルduを指数変換して、実際のフィードバック入力uを算出し(ステップS170)、電力配分制御部11に出力し、ステップS120に戻る。The control
このような電力配分制御における効果について、図6を用いて説明する。図において、横軸は時刻であり、縦軸は蓄電池の充電量である。実線は配電計画時点における予測供給量と予測需要量に基づく最適な充放電曲線Pを表す。 The effect of such power distribution control will be described with reference to FIG. In the figure, the horizontal axis is the time and the vertical axis is the charge amount of the storage battery. The solid line represents the optimum charge / discharge curve P based on the predicted supply amount and the predicted demand amount at the time of distribution planning.
実際の電力供給量が計画を上回る等の理由により、早期に蓄電池を充電したほうが目的関数を減少させられる状況が発生したとする。例えば、これまでの手法によれば、そのような状況が発生した時点(時刻t)で、最終時刻までの蓄電池の充放電計画を再計算しなければならない。しかし、本願の実施の形態1および以降の実施の形態2にかかる電力制御装置1または電力配分計画設定方法では、オンラインでのフィードバック制御によって、充放電計画を再計算することなく、目的関数(評価関数E)を減少させるような充放電曲線M(破線)を実現できる。またその際、蓄電池に関する制約(充電量は0から1の間)が考慮される。
It is assumed that the objective function can be reduced by charging the storage battery at an early stage because the actual power supply exceeds the plan. For example, according to the conventional method, the charge / discharge plan of the storage battery up to the final time must be recalculated at the time when such a situation occurs (time t). However, in the
ここで、配電系統800が蓄電池を含む場合、最適な配電量は蓄電池の将来にわたる充電状態を考慮して決定されるため、誤差の補正もまた蓄電池の将来にわたる充電状態を考慮したものでなければならない。それに対して、上述した配分調整部12を備えたことで、蓄電池を含む配電系統800において、発電量または電力需要の予測が実際の発電量または需要と一致しない場合においても、蓄電池の将来にわたる充電状態を考慮に入れたフィードバック制御が実現できる。つまり、本実施の形態にかかる電力制御装置1あるいは電力配分計画設定方法では、蓄電池の将来にわたる充電状態を考慮して、最適な配電量を決定するので、実態への適応性と制御性を両立させることができる。
Here, when the
なお、本実施の形態に示す電力制御装置1は、オンラインの制御則が対数変換部121、制御入力計算部122、指数変換部123で構成されており、これらは微分方程式の初期値問題、およびその少数回の繰り返しと、単純な代数演算である。従って、最適化問題をオンラインで解く必要がなく、オンライン制御に適したものとなっている。
In the
そして、本実施の形態で例示する配分計画部15と、配分調整部12(対数変換部121、制御入力計算部122、指数変換部123)は、例えば計算機内のプログラムとして実現されてもよい。その際、各ステップが必要とする情報がそれ以前のステップにて計算される限りにおいて、計算順序の入れ替えは許容される。また、指数変換部123、および対数変換部121における計算は、かならずしも指数写像と対数写像に限らなくてもよく、接ベクトルと多様体上の点を相互に滑らかに変換する写像であれば許容される。
The
その場合、配分計画部15と配分調整部12は、図7に示すように、プロセッサ101と記憶装置102を備えた一つのハードウェア10によって構成することが考えられる。記憶装置102は、図示していないが、ランダムアクセスメモリ等の揮発性記憶装置と、フラッシュメモリ等の不揮発性の補助記憶装置とを具備する。また、フラッシュメモリの代わりにハードディスクの補助記憶装置を具備してもよい。プロセッサ101は、記憶装置102から入力されたプログラムを実行する。この場合、補助記憶装置から揮発性記憶装置を介してプロセッサ101にプログラムが入力される。また、プロセッサ101は、演算結果等のデータを記憶装置102の揮発性記憶装置に出力してもよいし、揮発性記憶装置を介して補助記憶装置にデータを保存してもよい。
In that case, as shown in FIG. 7, it is conceivable that the
実施の形態2.
本実施の形態2においては、電力配分問題が入力に関する不等式制約k(xt,yt,ut)>0を含む場合の電力配分における計算方法について説明する。なお、本実施の形態2にかかる電力制御装置あるいは電力配分計画設定方法では、不等式制約k(xt,yt,ut)>0を含むことにより、配分計画部と配分調整部で用いる計算式が異なる以外は、実施の形態1と同様である。そのため、実施の形態1で用いた図、あるいは同様部分の式については、援用し、同様部分の説明は繰り返さない。Embodiment 2.
In the second embodiment, the calculation method in the power distribution when the power distribution problem includes the inequality constraint k ( xt , yt , ut )> 0 regarding the input will be described. In the power control device or the power distribution plan setting method according to the second embodiment, the calculation used in the distribution planning unit and the distribution adjustment unit by including the inequality constraint k ( xt , yt , ut )> 0. It is the same as the first embodiment except that the formula is different. Therefore, the figure used in the first embodiment or the formula of the same part is incorporated, and the explanation of the same part is not repeated.
蓄電池の残量に関する不等式制約等は、蓄電池の充放電電力、または充放電電流が変数utに含まれる場合、容易に入力に関する不等式制約に変換できる。本実施の形態2において、不等式制約は等式制約k(xt,yt,ut)-st=0、および不等式制約st>0に置き換えられる。ただしst>0とは、stの成分がすべて正であることを意味する。また、等式制約k(xt,yt,ut)-st=0は、l(z,u,st)=0に含めて取り扱う。The inequality constraint on the remaining amount of the storage battery can be easily converted into the inequality constraint on the input when the charge / discharge power or the charge / discharge current of the storage battery is included in the variable ut . In the second embodiment, the inequality constraint is replaced by the equality constraint k (x t , y t , ut ) −st = 0, and the inequality constraint s t > 0. However, s t > 0 means that all the components of s t are positive. Further, the equation constraint k (x t , y t , ut) -s t = 0 is included in l (z, u, st ) = 0.
さて、s>0なるsのなす多様体は、k≧0とすると、計量gij(s)=k+(si)-2を導入することで、リーマン多様体と考えることができる。組(zt,νt,st)のなす多様体は、stのなす多様体と、ユークリッド空間の直積であり、これもまたリーマン多様体である。このリーマン多様体における、関数Htの微分は、νt=(ut,λt,st)を入力だとみなして、式(17)と変わることを除き、実施の形態1と同様である。
∇ijHt=∂ijHt-Γij
k∇kHt ・・・(17)Now, assuming that k ≧ 0, the manifold formed by s of s> 0 can be considered as a Riemannian manifold by introducing the metric gij (s) = k + ( si ) -2 . The manifold formed by the set (z t , ν t , st ) is the direct product of the manifold formed by st and the Euclidean space, which is also a Riemannian manifold. The derivative of the function H t in this Riemannian manifold is the same as that of the first embodiment except that ν t = ( ut, λ t , st ) is regarded as an input and is different from the equation (17). be.
∇ ij H t = ∂ ij H t --Γ ij k ∇ k H t・ ・ ・ (17)
ただし、δijHtは、Htのνi、νjに関する2階の偏微分であり、Γij kは、第2種クリストッフェル(Christoffel)記号である。ここでは、νiがsjである場合、Γii i=-1/sjであり、それ以外の添字の組み合わせに関してΓij k=0である。本実施の形態2では、実施の形態1においては∇ijHt=δijHtであった点に対して異なる。また、w=(z,u,s)に関する測地線方程式については、ユークリッド計量δijが計量gijに置き換わることを除き、実施の形態1と同様である。However, δ ij H t is a second-order partial derivative with respect to ν i and ν j of H t , and Γ ij k is a second-class Christoffel symbol. Here, when ν i is s j , Γ ii i = -1 / s j , and Γ ij k = 0 for other combinations of subscripts. The second embodiment is different from the point that ∇ ij H t = δ ij H t in the first embodiment. The geodesic equation for w = (z, u, s) is the same as that of the first embodiment except that the Euclidean metric δ ij is replaced with the metric g ij .
なお、本願は、様々な例示的な実施の形態および実施例が記載されているが、1つ、または複数の実施の形態に記載された様々な特徴、態様、および機能は特定の実施の形態の適用に限られるのではなく、単独で、または様々な組み合わせで実施の形態に適用可能である。従って、例示されていない無数の変形例が、本願明細書に開示される技術の範囲内において想定される。例えば、少なくとも1つの構成要素を変形する場合、追加する場合または省略する場合、さらには、少なくとも1つの構成要素を抽出し、他の実施の形態の構成要素と組み合わせる場合が含まれるものとする。 Although various exemplary embodiments and examples are described in the present application, the various features, embodiments, and functions described in one or more embodiments are specific embodiments. It is not limited to the application of, but can be applied to the embodiment alone or in various combinations. Therefore, innumerable variations not exemplified are envisioned within the scope of the techniques disclosed herein. For example, it is assumed that at least one component is modified, added or omitted, and further, at least one component is extracted and combined with the components of other embodiments.
以上のように、各実施の形態にかかる電力制御装置1によれば、蓄電池を含む配電系統800をリーマン多様体とみなし、配電系統800内において予測される電力供給量(予測供給量Qp)と電力需要(予測需要量Dp)から、蓄電池の残量を並べたベクトルxt(配電系統800内で融通される電力等を並べたベクトルytとまとめてベクトルzt)と配電系統800内での制御可能な量を並べたベクトルutとを変数とする配電系統800内での電力損失を評価する関数Lと、蓄電池の目標充電量からのずれを示す関数Sの和を示す評価関数E(式(1))を最小化する電力配分計画u0:tf、z0:tfを生成する配分計画部15、配電系統800内における実際の配分状態ハットztを取得する配分状態取得部16、電力配分計画u0:tf、z0:tfに基づいて配電系統800内での電力配分を制御する電力配分制御部11、および電力配分計画u0:tf、z0:tfと配分状態ハットztとの間にずれが生じた際、電力配分計画と配分状態との距離(例えば、D(z,ハットzt))が最小となるよう、接ベクトルと点を相互に滑らかに変換する写像を用いて電力配分計画の調整値(フィードバック入力u)を生成する配分調整部12、を備えるように構成したので、実態に応じた見積もりに基づく制御を可能とする。As described above, according to the
とくに、配分調整部12には、配分状態ハットztの対数写像の値である第1の接ベクトルdzを計算する対数変換部121、対数変換部121によって計算された第1の接ベクトルdzに係数(フィードバック係数Ht,νν、Ht,νz)を乗じ、接空間におけるフィードバック入力である第2の接ベクトルduを計算する制御入力計算部122、および制御入力計算部122によって計算された第2の接ベクトルduの指数写像の値である調整値(フィードバック入力u)を計算する指数変換部123、が設けられているように構成すれば、微分方程式の初期値問題、の小数回の繰り返しと単純な代数演算で調整値(フィードバック入力u)を算出し、容易にオンライン制御できる。In particular, in the distribution adjustment unit 12, the
また、各実施の形態にかかる電力配分計画設定方法によれば、蓄電池を含む配電系統800をリーマン多様体とみなし、配電系統800内において、電力供給量(予測供給量Qp)と電力需要(予測需要量Dp)を予測する予測ステップ(ステップS100)、予測供給量Qpと予測需要量Dpから、蓄電池の残量を並べたベクトルxt(配電系統800内で融通される電力等を並べたベクトルytとまとめてベクトルzt)と配電系統800内での制御可能な量を並べたベクトルutとを変数とする配電系統800内での電力損失を評価する関数Lと、蓄電池の目標充電量からのずれを示す関数Sの和を示す評価関数E(式(1))を最小化する電力配分計画u0:tf、z0:tfを生成する電力配分計画生成ステップ(ステップS110(ステップS201~S206))、電力配分計画u0:tf、z0:tfに基づいて配電系統800内での電力配分を制御する電力配分制御ステップ(ステップS120)、配電系統800内における実際の配分状態ハットztを取得する配分状態取得ステップ(ステップS130)、および電力配分計画u0:tf、z0:tfと配分状態ハットztとの間にずれが生じた際、電力配分計画と配分状態との距離例えば、D(z,ハットzt)が最小となるよう、接ベクトルと点を相互に滑らかに変換する写像を用いて電力配分計画の調整値(フィードバック入力u)を生成する調整値生成ステップ(ステップS140~S170)、を含むように構成したので、実態に応じた見積もりに基づく制御を可能とする。Further, according to the power distribution plan setting method according to each embodiment, the
1:電力制御装置、 11:電力配分制御部、 12:配分調整部、 13:供給予測部、 14:需要予測部、 15:配分計画部、 16:配分状態取得部、 121:対数変換部、 122:制御入力計算部、 123:指数変換部、 800:配電系統、 D(z,ハットzt):距離、 Dp:予測需要量、 du:第2の接ベクトル、 dz:第1の接ベクトル、 Ht,νν,Ht,νz:フィードバック係数、 Qp:予測供給量、 u:フィードバック入力(調整値)、 xt,yt,zt:ベクトル、 u0:tf,z0:tf:電力配分計画、 ハットzt:配分状態。1: Power control device, 11: Power distribution control unit, 12: Distribution adjustment unit, 13: Supply forecast unit, 14: Demand forecast unit, 15: Distribution planning unit, 16: Distribution status acquisition unit, 121: Coefficient conversion unit, 122: Control input calculation unit, 123: Exponential conversion unit, 800: Distribution system, D (z, hat z t ): Distance, Dp: Predicted demand, du: Second tangent vector, dz: First tangent vector , H t, νν , H t, ν z : feedback coefficient, Qp: predicted supply amount, u: feedback input (adjustment value), x t , y t , z t : vector, u 0: tf , z 0: tf : Power allocation plan, hat zt : Allocation status.
Claims (3)
前記配電系統における実際の配分状態を取得する配分状態取得部、
前記電力配分計画に基づいて前記配電系統での電力配分を制御する電力配分制御部、および
前記電力配分計画が、前記配分状態に対してずれを有する際、フィードバック制御によって、前記ずれを有する電力配分計画と前記配分状態との距離が最小となるよう、接ベクトルと点を相互に滑らかに変換する写像を用いて前記ずれを有する電力配分計画の調整値を生成する配分調整部、
を備えたことを特徴とする電力制御装置。 The distribution system including the storage battery is regarded as a Lehman variety, and the vector in which the remaining amount of the storage battery is arranged and the controllable amount in the distribution system are arranged from the power supply amount and the power demand predicted in the distribution system. A distribution planning unit that generates a power distribution plan that minimizes the evaluation function that indicates the sum of the function that evaluates the power loss in the distribution system and the function that indicates the deviation from the target charge amount of the storage battery.
Distribution status acquisition unit that acquires the actual distribution status in the distribution system,
When the power distribution control unit that controls the power distribution in the distribution system based on the power distribution plan and the power distribution plan have a deviation from the distribution state , the power distribution having the deviation is performed by feedback control. A distribution adjustment unit that generates an adjustment value for a power distribution plan with the deviation using a mapping that smoothly transforms a tangent vector and a point to each other so that the distance between the plan and the distribution state is minimized.
A power control device characterized by being equipped with.
前記配分状態の対数写像の値である第1の接ベクトルを計算する対数変換部、
前記対数変換部によって計算された第1の接ベクトルに係数を乗じ、接空間におけるフィードバック入力である第2の接ベクトルを計算する制御入力計算部、および
前記制御入力計算部によって計算された第2の接ベクトルの指数写像の値である前記調整値を計算する指数変換部、
が設けられていることを特徴とする請求項1に記載の電力制御装置。 In the distribution adjustment unit,
A logarithmic transforming unit that calculates a first tangent vector that is the value of the logarithmic map in the distributed state.
A control input calculation unit that calculates a second tangent vector that is a feedback input in tangent space by multiplying the first tangent vector calculated by the logarithmic conversion unit by a coefficient, and a second tangent vector calculated by the control input calculation unit. Exponential conversion unit that calculates the adjustment value, which is the value of the exponential map of the tangent vector of
The power control device according to claim 1, wherein the power control device is provided.
前記予測された電力供給量と電力需要から、前記蓄電池の残量を並べたベクトルと前記配電系統での制御可能な量を並べたベクトルとを変数とする前記配電系統での電力損失を評価する関数と、前記蓄電池の目標充電量からのずれを示す関数の和を示す評価関数を最小化する電力配分計画を生成する電力配分計画生成ステップ、
前記電力配分計画に基づいて前記配電系統での電力配分を制御する電力配分制御ステップ、
前記配電系統における実際の配分状態を取得する配分状態取得ステップ、および
前記電力配分計画が、前記配分状態に対してずれを有する際、フィードバック制御によって、前記ずれを有する電力配分計画と前記配分状態との距離が最小となるよう、接ベクトルと点を相互に滑らかに変換する写像を用いて前記ずれを有する電力配分計画の調整値を生成する調整値生成ステップ、
を含むことを特徴とする電力配分計画設定方法。 A prediction step in which a distribution system including a storage battery is regarded as a Riemannian manifold and power supply amount and power demand are predicted in the distribution system.
From the predicted power supply amount and power demand, the power loss in the distribution system is evaluated with the vector in which the remaining amount of the storage battery is arranged and the vector in which the controllable amount in the distribution system is arranged as variables. A power distribution plan generation step that generates a power distribution plan that minimizes an evaluation function that indicates the sum of the function and the function that indicates the deviation from the target charge of the storage battery.
A power distribution control step that controls power distribution in the distribution system based on the power distribution plan,
When the distribution state acquisition step for acquiring the actual distribution state in the distribution system and the power distribution plan have a deviation from the distribution state, the power distribution plan having the deviation and the distribution state are subjected to feedback control. An adjustment value generation step that generates an adjustment value for a power distribution plan with the deviation using a mapping that smoothly transforms the tangent vector and the point to each other so that the distance between them is minimized.
A power distribution plan setting method characterized by including.
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2019/028395 WO2021014483A1 (en) | 2019-07-19 | 2019-07-19 | Power control device and power distribution plan setting method |
Publications (2)
Publication Number | Publication Date |
---|---|
JPWO2021014483A1 JPWO2021014483A1 (en) | 2021-11-18 |
JP7068557B2 true JP7068557B2 (en) | 2022-05-16 |
Family
ID=74193734
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2021534845A Active JP7068557B2 (en) | 2019-07-19 | 2019-07-19 | Power controller and power distribution plan setting method |
Country Status (2)
Country | Link |
---|---|
JP (1) | JP7068557B2 (en) |
WO (1) | WO2021014483A1 (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114362262B (en) * | 2021-12-02 | 2023-09-29 | 国网电力科学研究院有限公司 | Voltage regulating system and regulating method based on active self-adaptive reduction |
JP2024008336A (en) * | 2022-07-08 | 2024-01-19 | 株式会社日立製作所 | Information processing apparatus |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003030172A (en) | 2001-07-11 | 2003-01-31 | Canon Inc | Optimization method with constraint condition and program |
JP2014143835A (en) | 2013-01-24 | 2014-08-07 | Toshiba Corp | Control system for power system |
JP2016163488A (en) | 2015-03-04 | 2016-09-05 | 株式会社東芝 | Device, method and program for power control |
CN108306297A (en) | 2018-01-11 | 2018-07-20 | 河海大学 | A kind of optimal load flow computational methods based on Riemann integral |
-
2019
- 2019-07-19 WO PCT/JP2019/028395 patent/WO2021014483A1/en active Application Filing
- 2019-07-19 JP JP2021534845A patent/JP7068557B2/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003030172A (en) | 2001-07-11 | 2003-01-31 | Canon Inc | Optimization method with constraint condition and program |
JP2014143835A (en) | 2013-01-24 | 2014-08-07 | Toshiba Corp | Control system for power system |
JP2016163488A (en) | 2015-03-04 | 2016-09-05 | 株式会社東芝 | Device, method and program for power control |
CN108306297A (en) | 2018-01-11 | 2018-07-20 | 河海大学 | A kind of optimal load flow computational methods based on Riemann integral |
Also Published As
Publication number | Publication date |
---|---|
WO2021014483A1 (en) | 2021-01-28 |
JPWO2021014483A1 (en) | 2021-11-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Montoya et al. | Sequential quadratic programming models for solving the OPF problem in DC grids | |
JP5178242B2 (en) | Energy storage device operation plan creation method and operation plan creation device | |
JP7068557B2 (en) | Power controller and power distribution plan setting method | |
CN105631528A (en) | NSGA-II and approximate dynamic programming-based multi-objective dynamic optimal power flow solving method | |
Adetunji et al. | An optimization planning framework for allocating multiple distributed energy resources and electric vehicle charging stations in distribution networks | |
Kouvaritakis et al. | Recent developments in stochastic MPC and sustainable development | |
CN109980636A (en) | Based on the geomantic omen fire coordination optimization dispatching method for improving Benders decomposition method | |
Chen et al. | An iterative and hierarchical approach to co-optimizing the velocity profile and power-split of plug-in hybrid electric vehicles | |
Jang et al. | Optimal HVAC system operation using online learning of interconnected neural networks | |
CN108154266B (en) | Economic dispatching method for wind power-containing power system | |
Behtash et al. | A Reliability-Based Formulation for Simulation-Based Control Co-Design Using Generalized Polynomial Chaos Expansion | |
Pinto et al. | Reliability‐constrained robust expansion planning of active distribution networks | |
Liu et al. | On decentralized optimization for a class of multisubsystem codesign problems | |
Jiménez et al. | A closed-loop data-driven optimization framework for the unit commitment problem: A Q-learning approach under real-time operation | |
Clever et al. | Combination of an adaptive multilevel SQP method and a space-time adaptive PDAE solver for optimal control problems | |
An et al. | Go beyond black-box policies: Rethinking the design of learning agent for interpretable and verifiable HVAC control | |
JP2001255903A (en) | Process control system | |
CN115136438A (en) | Distributed resource management device and distributed resource management method | |
CN115313422B (en) | Operation control method of electricity-gas comprehensive energy system for actively supporting main network frequency | |
Li et al. | Real-time dispatching performance improvement of multiple multi-energy supply microgrids using neural network based approximate dynamic programming | |
Kippelt et al. | Assessment of capacity management approaches in distribution grids | |
Dey et al. | Reinforcement Learning Building Control: An Online Approach With Guided Exploration Using Surrogate Models | |
Roy et al. | Trust-region based algorithms with low-budget for multi-objective optimization | |
CN114039346A (en) | Voltage regulation and current sharing method based on neural network | |
JP2022077459A (en) | Electric power distribution device, electric power distribution method, and program |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20210705 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20210705 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20220405 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20220428 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 7068557 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |