Simulink 環境における小規模マイクログリッドのシミュレーション第 2 報 - フェーザー法による 24 時間のシミュレーション - 三田宇洋, テクニカルコンサルティング部, MathWorks Japan 1. 導入文献 [1] では 実用的な時間でシミュレーションを行う小規模のマイクログリッドのシミュレーションモデルを紹介した モデルはコンバータやインバータとそのスイッチング制御を省略し 電源 負荷を電流源扱いとしたことで シミュレーションを高速化し 数十秒規模のシミュレーションを行った さて 太陽光発電は 太陽を含めた外的環境により発電量が変動する 太陽光発電がマイクログリッドに及ぼす影響を太陽の挙動のサイクルで見るためには 評価する時間を少なくとも 24 時間としてシミュレーションを実行する必要がある 24 時間規模のシミュレーションを現実的な時間で行うためには 上記のモデルとシミュレーション手法では難しい Simulink の電気系ブロックライブラリの SimPowerSystems TM は シミュレーションの方法として システムの特定周波数 ( 例えば系統周波数 5[Hz]) の振幅と位相に注目するフェーザー法を提供する フェーザー法は SimPowerSystems のブロックにおいて振幅と位相以外の情報はシミュレーションから除外されるので 大幅なシミュレーション時間の改善が期待できる 本報告では 先報 [1] のモデルをフェーザー法に変更し 24 時間のシミュレーションを実行する シミュレーションでは マイクログリッドの各要素の供給 / 消費電力 蓄電池の使用量及び系統に売電できる金額を定量化する 以下にモデルとシミュレーションについて報告する 2. モデル化するマイクログリッド図 1 に想定するモデル化用のマイクログリッドを示す マイクログリッドは文献 [1] と同様である 単相交流 (2[V]) の電力ネットワークをマイクログリッドとする 再生可能エネルギーとして 太陽光発電を想定する 電力源は 系統電力 太陽光発電及び蓄電池である 蓄電池はマイクログリッド内の余剰電力の吸収 / 不足電力の供給を行う 電力を消費する負荷は一般家庭 3 軒である マイクログリッドは系統と 1 点 ( 柱上変圧器 ) で接続する マイクログリッドの定格を表 1 に示す 系統電力は 三相交流の電圧源 (66[kV]) が変圧器 (66[kV]/6.6[kV]) に接続する 高圧 (66[kV]) の電力は変圧器で降圧 (6.6[kV]) される 伝送路を介して 柱上変圧器 (6.6[kV]/2[V]) により 単相交流 (2[V]) に変換される 交流周波数は 5[Hz] とする 単相交流の電力ネットワークに 電力源の太陽光発電 ( 直流 ) 及び蓄電池 ( 直流 ) が接続される 1
三相 AC 66k[V] 変圧器系統電力三相電圧源柱上変圧器伝送路 : 電力供給網 ( 単相 AC) : 分散配置をつなぐネットワーク : 系統電源 ( 三相 AC) 三相 AC 6.6k[V] 一般家庭 1 太陽光発電 マイクログリッド 一般家庭 2 単相 AC 2[V] IT 制御システム 一般家庭 3 蓄電池 図 1 小規模マイクログリッド [1] 表 1 定格 施設太陽光発電 1 台蓄電池 1 台一般家庭 3 軒 定格最大 5[W] 容量 1[Ah] 最大 25[W]/1 軒 2
3. シミュレーションモデルシミュレーションモデルを図 2 に示す シミュレーションモデルの骨格は文献 [1] とほぼ同様である 電流の向きは以下とする 太陽光発電マイクログリッドへの流入 : 正 蓄電池マイクログリッドへの流入 : 正 一般家庭 1,2,3 電力を消費する方向 : 正 柱上変圧器 2 次側 : 系統へ流出する方向 : 正 文献 [1] と異なる変更点について説明する 変更点 (1) シミュレーションモードをフェーザー法に変更 (2) シミュレーション時間を 24 時間に変更 (3) (2) に伴いシミュレーションのシナリオを変更 (4) 蓄電池のモデルを変更 [2] 柱上変圧器 [3] 太陽光発電の電流注入 [1] 系統 [4] 蓄電池の放電 / 充電 [6] 一般家庭 1 [7] 一般家庭 2 [8] 一般家庭 3 [9] 蓄電池制御器 [5] 蓄電池の特性 [1] 注入電流量計算 (Phasor 法 ) 図 2 モデルの全体図 3
変更点 (1) シミュレーションモードをフェーザー法に変更 フェーザー法に対応するよう 電圧と電流の信号を振幅と位相から複素数表現に変換する ( 図 3 参照 ) (a) 太陽光発電電流値流入 (b) 蓄電池電流値流入 (c) 負荷 1 電流値流入 (d) 負荷 2 電流値流入 (e) 負荷 3 電流値流入 (b) 蓄電池電流値 (a) 太陽光発電電流値 (c) 負荷 1 電流値 (d) 負荷 2 電流値 (e) 負荷 3 電流値 フェーザー法で扱うため交流電流信号を複素数に変換 図 3 フェーザー法への信号変換 (2) シミュレーション時間を 24 時間に変更 (3) (2) に伴いシミュレーションのシナリオを変更想定する太陽光発電の時系列データを図 4 に示す 一般家庭の負荷の時系列データを図 5 に示す 3 軒の一般家庭の負荷は 3 軒とも同様に図 5 の負荷変動をとるものとする これらの時系列データがシミュレーションモデルに与えられる 図 4, 図 5 の形状は文献 [2] を参考にする 4
Load Power[W] Solar Power[W] 55 5 45 4 35 3 25 2 15 1 5 5 1 15 2 25 Time[Hour] 図 4 一日の太陽光発電の時系列データ 25 2 15 1 5 5 1 15 2 25 Time[Hour] 図 5 一日の一般家庭の負荷変動 (4) 蓄電池のモデルを変更 SimPowerSystems で提供される Battery Block は フェーザー法に対応しない そこで 蓄電池の充電状態を示す SOC(State Of Charge) と蓄電池の放電容量 [Ah] を示す簡単なモデルを Simulink で作成し使用する モデルを図 6 に示す 5
容量 Ah Ah SOC[%] 1 t i t 36 t i dt t dt q Q *36 q i t dt t Q: 最大蓄電池容量 [Ah] 注 : 下記モデルでは SOC の初期値 =8[%] にするよう q()=q*36*.2 とした 注 : 図 6 蓄電池の容量と SOC を表現する簡易モデル 4. シミュレーションのシナリオについて図 7 に一般家庭の総負荷と太陽光発電量を示す 区間 () では 太陽光発電量が一般家庭の総負荷を上回る 区間 () 近傍の 12 時から 18 時まで ( 区間 (1)) と区間 (1) 以外で分け 以下の異なる制御則を適用する 制御則 [1] 12 時から 18 時まで ( 区間 (1)) 以下 [2] の蓄電池制御を実行しない 蓄電池が余剰電力の吸収 / 不足電力の放出をしないので 太陽光発電と系統電力で負荷の需要をまかなう つまり 系統との売電 / 買電を行う マイクログリッドに余剰電力がある場合 系統に電力が流出し マイクログリッドの電力が不足する場合 系統から電力を吸収する [2] 区間 (1) 以外蓄電池制御を実行する 蓄電池制御はマイクログリッドと電力系統間で流入 / 流出する有効電力を に追従させる制御である 太陽光発電と蓄電池の電力で 負荷の需要をまかなう つまり 系統との売電 / 買電を行わない 太陽光発電の電力 > 負荷の電力の場合 余剰電力が発生する 蓄電池が余剰電力を吸収する 太陽光発電の電力 < 負荷の電力の場合 負荷の電力が不足する 蓄電池が不足電力を放出する 6
Power [W] 8 7 6 余剰電力 Total Power of Load Solar 不足電力 5 4 3 2 1 区間 (1) 区間 () -1 5 1 15 2 25 Time [h] 図 7 一般家庭分の総負荷と太陽光発電量 5. シミュレーション結果シミュレーション結果を図 8 に示す このときの蓄電池のステータスを図 9 に示す 図 8 図 9 の各データについて 表 2 に説明を示す なお シミュレーション時間は以下の環境で次を得た アクセラレータモード :39 分 環境 プロセッサ :CoreTMi7-264M CPU @ 2.8GHz 2.8GHz RAM: 8.GB System: 64 bit operating system OS: Windows 7 Enterprise 7
SOC [%] Ah 積算 [Ah] Power Price[JPYEN] Sum(Utility Power)[kWh] Utility Power [W] Power [W] Power [W] 1 5 一般家庭総負荷太陽光発電 -5 5 1 15 2 25 Time [h] 1 蓄電池 5-5 5 1 15 2 25 5 系統 買電 (3) 放電 充電 売電 (2) -5 5 1 15 2 25 15 1 5-5 5 1 15 2 25 6 4 2 系統電力積算 金額 -2 5 1 15 2 25 Time [h] 3 25 2 15 1 5 区間 (1) 図 8 シミュレーション結果 Ah 積算 12.3[kWh] 517[ 円 ] 5 1 15 2 25 8 75 7 65 6 55 区間 (1) SOC 5 5 1 15 2 25 Time [h] 図 9 蓄電池のステータス 8
表 2 図 8 図 9 の説明 名称蓄電池系統系統電力積算金額 Ah 積算 SOC 説明区間 (1) では蓄電池制御が実行されず 蓄電池の電力は となる 区間 (1) 以外では 一般家庭の総負荷と太陽光発電の大小に応じて マイクログリッドの余剰 / 不足電力を充電 / 放電する 区間 (1) では蓄電池制御が実行されず 系統が余剰 / 不足電力を補償する 売電 (2) では 太陽光発電 > 一般家庭の総負荷となり 余剰電力が系統に流出する 買電 (3) では 太陽光発電 < 一般家庭の総負荷となり 不足電力が系統から流入する 系統電力は 24 時間で 12.3[kWh] となる 売電金額は 24 時間で 517 円となる 但し売電単価 =42[ 円 /kwh] とする 蓄電池がどれだけの電流を放出したかの指標である 24 時間で積算量は 25[Ah] となる SOC は時刻 の初期値 8[%] から 24 時間で 55[%] 程度まで低下する 6. 使用ツール本モデル作成に使用した MATLAB プロダクトを表 3 に示す 表 3 ツールの役割ツール (Version:R212a) 役割 Simulink シミュレーション基本環境 制御器モデリング MATLAB 基本環境 SimPowerSystems 系統システム 電流源で表現した太陽光発電 / 蓄電池 電気系機器のモデリング Simscape TM 物理モデリングツールの基本環境 7. 結語太陽光発電 蓄電池を系に含む小規模マイクログリッドのモデルの 24 時間のシミュレーションを紹介した モデルは SimPowerSystems により作成し フェーザー法を適用する 太陽光発電量が大な 12 時から 18 時は売電し 太陽光発電量が小な時間 ( 時から 12 時 18 時から 24 時 ) は蓄電池が系の不足電力を供給する簡単な制御則により 蓄電池の充電 / 放電を制御する 蓄電池はフェーザー法で使える簡易なモデルを作り 放出電流の積算値 (Ah) や SOC を計算する シミュレーションでは マイクログリッドに接続する電源 負荷の供給 / 消費電力 蓄電池の使用量及び系統に売電できる金額が定量化された これらの値は物理的に妥当な値が得られ シミュレーション実行時間は実用に耐えうることが確認された 9
8. 付録フェーザー法について [3] 線形回路において 電圧や電流の振幅や位相の変化だけに注目したい場合がある この場合 R( レジスタンス ) L( インダクタンス ) C( キャパシタンス ) の素子間での相互作用から得られる微分方程式 ( 状態空間モデル ) をすべて解く必要はない 代わりに 電圧と電流の複素表現を使った簡単な代数方程式系を解くだけでよい これを フェーザー法と呼ぶ フェーザー法では電圧と電流を複素表現で計算する ある周波数における一定の周期性 ( 正弦波 余弦波 ) をもつ電圧と電流を複素数で表現する これらは また直交座標系 ( 実数部と虚数部 ) でも 極座標系 ( 大きさと位相 ) でも表現することが可能である フェーザー法では 電気系の状態量は無視されるので システムの電気部品個々の動作をシミュレーションする必要はない そのため このシミュレーションの実行は高速になる ただし フェーザー法は ある 1 つの周波数についてだけの解を得るものである また回路のスイッチング等の非線形要素には対応しない 9. 参考文献 [1] 三田,Simulink 環境における小規模マイクログリッドのシミュレーション, 技術資料,MathWorks Japan,212 [2] 太陽生活ドットコム, http://taiyoseikatsu.com/, 212 [3] SimPowerSystems User s guide,r211b, pp.1-42, MathWorks,211 1