ラグランジュの未定乗数法を用いた等増分燃料費の法則

問題

【問題】ある電力系統に火力発電所が2基連携している。それぞれの発電電力[MW]をP1P2、発電費用[円]をF1(P1)F2(P2)とし、それぞれ下記で表すことができるとする。
(1)F1(P1)=100+8P1+6P12(2)F2(P2)=200+4P2+2P22
負荷で消費する電力が400[MW]であり、P1+P2=400が常に成立するとき、合計の発電費用が最小となるP1P2を求めよ。

はじめに

電力関係の問題の解説は、先駆者様のサイトが多々あります。本分野において、新規性のある記事がどこまで書けるか未知ですが、第1回目は等増分燃料費の法則について扱います。

作成した理由は下記です。

  1. 等増分燃料費の法則の式の結果のみを使用して解説していることが多いから
  2. 背景、イメージ論からの説明が少ないと感じたから

確かに、工学系の技術者としては結果だけ知っていれば与えられた問題に対して答えを出すことはできます。しかし、それが正しいことを他者に説明するうえでイメージ、理由について語ることはほぼ必須です。

こういった時に躓かないように、数学的な素養を身に着ける材料の一つになれば良いなと考え解説することにしました。

本記事は4章に分かれていますが、2章目の解説がネットに少ないと感じます。是非そこだけでも参考に下さると幸いです。

等増分燃料費の法則とは

(3)F1P1=F2P2==FiPi

の関係式になります。

それぞれの発電機のある負荷点における発電費用の微分係数が全て一致するとき、制約条件下で燃料費が最小になる。ことを示しています。

ラグランジュの未定乗数法という数学的法則から立式できます。(下記)

ラグランジュの未定乗数法

P1Piを変数に持つ関数L(P1,P2,Pi,λ)が制約条件g=0の下で極値を取るとき、その極値を与えるP1,P2Piは以下の条件で求められる。

(4)P1(Lλg)=P2(Lλg)=Pi(Lλg)=0

これだけ見るとよく分からないと思います。本問を例に説明していきます。

ラグランジュの未定乗数法を用いて例題を解いてみる

本問の目的は、総発電費用L=F1(P1)+F2(P2)を制約条件g=P1+P2=400の下で最小化することです。

これにより、偏微分対象のラグランジュ関数Lλgは下記で与えられます。

(5)Lλg=300+8P1+6P12+4P2+2P22λ(P1+P2300)

これを各変数P1P2λで偏微分すると以下になります。

(6)8+12P1λ=0(7)4+4P2λ=0(8)P1+P2400=0

(6),(7)式目を使用しλを削除すると、P1P2の関係式を求められます。

その結果を(8)式に代入すると、答えが出そうです。

しかし、これらの式はどこかで見たことがありませんか?

ラグランジュの未定乗数法で立式した式の意味

まず、(6),(7)式ですが、等増分燃料費の法則を表します。

前章の説明の通り、λを削除すると、以下の関係式が生まれます。

(9)8+12P1=4+4P2=F1P1=F2P2

前節では、制約条件を加味したラグランジュ関数Lλgを作成しました。しかし極値を与える変数を求める際に、結局λは消えます。

この結果から、極値を与える変数を求めるだけならば、ラグランジュ関数を作成せずとも各Fi(Pi)Piで偏微分して=で結ぶだけで良いことが分かります。

これが、等増分燃料費の法則導出の背景です。

(8)式は、問題文中に与えられた制約条件で最初から分かっています。未定乗数法を使わずとも分かる式です。

解答例とそのイメージ

大体答えは出たようなものですが、最後の値まで算出します。

解答例

(8)式を整理して、(10)12P14P2=4(11)P2=3P1+1

これを制約条件P1+P2=400に代入して

(12)4P1=399(13)P1=99.75

これより、P2=300.25であることが分かりました。

イメージ

「等増分燃料費の法則」の名前から察する方もいらっしゃるかもしれません。

解答例で求めたP1P2は、その点におけるF1(P1)F2(P2)の傾きが等しいことを意味しています。

P1P2をさらに発電費用を最小化できる点が無いか点を動かすが発電機1の発電費用の増加/減少分に対し、発電機2の費用の減少/増加分が不足する。

➡合計発電費用は高くなってしまうことが分かります。

これにより、極値を与える点が発電費用として最小であることが分かります。

補足

ラグランジュの未定乗数法は最小値を求めるものではない

ラグランジュの未定乗数法は、ある制約条件下で「極値」を与える変数を算出することができる方法です。

今回、例題で与えた関数は下に凸だったため、極小値しかありませんでした。ですので、最小値が求まります。

逆に、上に凸のグラフが与えられたときは、極大値が求まり、最大値が求められます。

得られたパラメータは現実的に取りうる値とは限らない

例えば、負荷で電力消費が無い状態P1+P2=0を考えます。

このとき、発電費用を最小にするP1,P2は、P1=P2=0より300円と計算できます。

しかし、ラグランジュの未定乗数法で計算してみるとどうなるのでしょうか。

(10)式より、4P1=1P1=0.25と負の値になります。問で考える対象は発電機のため、負に発電することは一般的に考えていないです。

このように、変数が取りうる値の範囲外で極値を求める可能性もあります。

また、元も子も無いですが、発電しないならば、そもそも発電機を停止するという考え方もあります。

\(F{P_{1},F{P_{2}\)の定数項(100,200)は、発電機を駆動させるための補機電力を指しています。発電機を停止するならば、補機を動かす必要もなくなり、\(F{P_{1},F{P_{2}=0)\になる算段です。

(ただ、再始動には大きな電源を必要とします。発電しないからと言って補機を停止することが最適解とは限りません。)

最後に

本問では、ラグランジュの未定乗数法を用いた最適発電負荷配分と留意事項を述べました。

資格試験で出題される可能性は低いですが、電気設備設計の実業務を行う際に使用するかもしれません。

タイトルとURLをコピーしました