• 検索結果がありません。

デュアルレートサンプリング適応形負荷周波数制御: University of the Ryukyus Repository

N/A
N/A
Protected

Academic year: 2021

シェア "デュアルレートサンプリング適応形負荷周波数制御: University of the Ryukyus Repository"

Copied!
7
0
0

読み込み中.... (全文を見る)

全文

(1)

Title

デュアルレートサンプリング適応形負荷周波数制御

Author(s)

山下, 勝己; 平安, 充; 宮城, 隼夫

Citation

琉球大学工学部紀要(42): 85-90

Issue Date

1991-09

URL

http://hdl.handle.net/20.500.12000/5493

Rights

(2)

85

A Load Frequency Self-Tuning Regulator

Based on Dual-Rate Sampling

Katsumi YAMASHITA*.

Mituru HIRAYASU··

Hayao MIYAGI*

Summary

This paper presents a new method of designing a load frequency

self - tuning regulator for power system based on dual- rate sampling.

The method is applied to a two - area power system provided with

reheat

type

turbines and

then

the control

effect of

the

proposed

regulator is compared with that of standard sampling self - tuning

regulator using digital simulasions.

Key

Words : Self-Tuning Regulator, Load Frequency Control. Dual- Rate

Sampling, Power Systems

77 '/

H7)ib4'<1M:~;t, m~Q)m([!~H:cJ:~TlI{U-C~tet~' ~ft.~T~:6~. ~Q)~Itb7.)-:6~1tt2l'RJJ\~~

'C

~ ~:~;t)f/I

'litQ)

mJ~

?'

1 '/

QY7 -{ -

t:

I~'7~

1tfJt1V:n;t-c

-t

1M.Ii!

~

ftlJfJ1IM:W~~Q':: c:6~-C~

Q. l-:6'

l-t~:6~;.

77 './

~Q)~~M:Q)M.~~+~-C~~~~, 77,/~Q)~ t.1M::6~J.QJftQ)~fl::':cJ: ~ ~Jltrt Q1J}-f;~:~;t, fim~?'

1

'/Q)-;J -{ -

t:

I~'7~i1iUfl1J-C,;t~-r

L--t

llGtetftllJfao4'<1M:

~MIDET ~

.:

c~;t

-C

~ t~ ~

\,

JlI~, '::Q)rD'RlH:~~T~~~'il.1J~mQ)ft~.m!~ ft~~Mm~~l-~~, 77,/~Q)D~M:Q)~.~~ tftjlJ~.mt~*

'/71

'/tA~: ~IbPJ]ll,~Q;a~~A OO.mliltlliltlUfao~Q)i}f~:6~!l~

<

.i3'~n-C~'

Q(1)-(8).

L-~L-~~;. ~mmQ)~<~.~cftjlJ~Q)~'/7~

'./~mWJ:6~(UJ--c~Q~~. ~tA~'./71) './~PI1Wj

(1.5~ Ta~ 2.5)011 J:cls~\-C

,;t

(UJ~':~

l- -C

~fl::~

7)v~~tQflIaH:~Q.

-:no

m~'~'/7I) '/~/li3

Wj,:1-t

l--C~;t,

• •

t~ftllJ_fl:J:cJ:~ftlJ~mQ)~fl::~

~ tQMIaH:~Q.

*Mre~;t.

;:

Q)rD'IlH:~~TQt=~.mc 1IiIJti1lQ)~

'./7~'./~/li3M~~~Q7~7N~-~~'/7~'/~

~#

:

1991~5f.113B

I~IB.T·

fWftI!!f'::W

Dept. of Electronics and Information Eng., Fac. of Eng.

(3)

デュアルレートサンプリング適応形負荷周波数制御:山下・平山・宮城 86 適応形負荷周波数制御器を文献⑪を基盤に構築する とともに,本制御器を再熱式火力系からなる2地域電 力系統モデルに適用し,本制御方式の有効性を周波数 偏差および連系線潮流偏差などの時間応答波形より検 証する. j毎に次のようにまとめる. Andu(t-(k-1)j-1)=…=Andu(t+kj)(4) (k=1,2]…Ⅲs) 但し1sは(m+j-1)/jを切り上げした整数値 である.このとき,(2)式に対するディアルレートサ ンプリングモデルが次式のように得られる. 2.制御則の決定 再熱式火力系統からなる2地域電力系統の第i地 域の動特性式は,地域制御誤差βlAfi(t)+ APtioi(t)をyi(t)に,スピードチェンジャ位置 の基準値からの偏差をu,(t)にとることにより次式 のように与えることができる.

X(q-1)Andy(t+j)=F(q-j)Andu(t)(5)

但し,

X(q-1)=1+五jq-j+颪zq-j-I+…+anq-n-j+’

へ-

ri(q-j)=Ii,+F2q-j+Br3q-2j+…+Bsq-(s-1》’

A(q-L)y(t)=B(q-1)u(t)+C(q~')の(t)(1) 但し, A(q-`)=1+alq-1+a2q-2+…+anq-n B(q~')=b1q-1+b2q-‘+…+bmq-m c(q-1)=CO+clq-'十czq-z+…+cpq-P 次に,(5)式に対する評価関数を次式で定義する. J=y(t+j)2+r(AndU(t))2 (6) 但し,rは制御量に対する重みを示す. なお,上記の評価関数は制御入力の代りに制御入力の 変化率を抑えることにより,間接的に大きな制御入力 の抑制をはかっている このときⅢ評価関数Jを最小にする最適制御入力は. 評価関数Jを制御入力AndU(t)で偏微分した. 上式において!u(t)およびy(t)は時刻tでのプラ ントの入出力を表し,⑩(t)はn.次の時間多項式の 負荷外乱を表す.なおり記述の煩雑さをさけるため添 字iを省略している. 今,(1)式にt=t+jを代入して得られる関係式と(1) 式との差のn.階差分をとることにより,次式のよう な負荷外乱項が除去された入出力関係式が得られる. 0J

=25(l-ZI(q-1)And)y(t+j)

OAndU(t)

+2B,(百(q-l)AndU(t))+2rムndU(t)(7)

A(q~')Andy(t+j)=B(q~j)Andu(t+j)(2) 但し, And=(1-q-j)、。 を零と極き]AndU(t)について解くことにより,次 式のように求められる

一石,(1-X(q-1)And)y(t+j)

AndU(t)= (8) r+blB(q-」) 更に,文献⑪の考え方に基づき次式で示す なお、上記のパラメータaj,a,!…,an1bl,b2,…, bsは各サンプリング周期毎にオンラインで同定され,

同定されたパラメータa1,a2,…,且、,Bi,B2,…,

bsを用いて(8)式の制御則を計算することになる. それ故,パラメータ推定計算を実施するため(5)式を P(q-j)=1+plq-l+p2q-`+…+p」-1q-q-l) (3) を(2)式の両辺に乗じAndy(t+j ̄1),Andy(t+j-2), …,Andy(t+1)の項を除去すると共に,制御入力を

(4)

琉球大学工学部紀要第42号,1991年 87 次式のように変形する. Andy(t)=OTG(t-j)(9) 但し,

0?=[且,,且幽…Ⅲanbhb2b…,6s]

竺悸H

ど(t-j)T=

[_Andy(t-j)…-A、。y(t-n-j+1)

AndU(t-j)…Andyq-Sj)] 【】 このとき.未知パラメータOは次式の繰返し形パラ メータ推定アルゴリズムにより同定される. Fig.1Blockdiagramo『atwo-areapower system. r(t-j)ど(t-j)

o(t)=o(t-j)+,+E(t-j)でTq-j)ど(t-j)

xundy(t)-0(t-j)6(t-j)]⑩

但し, TablelSystemparameters. 「(t-j)e(t-j)S(t-jyr(t-j)

「(t)=「(t-j)-,+G(t-j)で「(t-j)G(t-j)

O<八(t)≦1,0≦ん(t)<2.「⑪>0 3.シミュレーション結果 本制御方式の有効性を検証するために,図1に示す 再熱式火力系からなる2地域電力系統を対象モデルと する.なお,各パラメータ値については表1に示す. また,シミュレーションに用いる負荷外乱を次式で示 すステップ負荷外乱とする.

△u(1)一一:÷[伽畠且,q-"M)

r+b,

+三,`,q-…Au(')]鰯

但し,

(0≦t<40)⑪

(40≦t) 0.01puMW qO2puMW

AP(t)

△y(t)=y(t)一yq-j)IAML)=u(1)-u(t-j)

なお,推定パラメータの初期値は6-1しと,それ

以外は零とした.また,行列「の初期値は「(0)=I ×103,町(t)およびル(t)の値をそれぞれ1と このとき,最適制御則は(8)式に、。=1を代入した次 式となる

(5)

デュアルレートサンプリング適応形負荷周波数制御:山下・平山・宮城 88 した. まず,観測と制御のサンプリング周期が同一である シングルレートサンプリング[No.1]およびサンプリ ング周期が異なるディアルレートサンプリング[No.2] の適応形負荷周波数制御器を、=m=3および「=1 に対して櫛簗する.このとき観測サンプリング周期 Tsに対する制御器No.1およびNo.2の制御効果を 調べるため,次式の2乗積分を定義する 0000000000 000000000 ■P■■凸●●●■ 432101234 『|』『 一 刺‐三×閏ニロコ

liWpやG W鞭1

'1

00 01 ]08.0m

PI-/WI[(“(`)ド+(APM。((〃]。'鰯

. 000000000 543210987 ●●● 句●□□●● I11111000 、Iご【×『H{ Fig3ResponsesofAf,. irate 0.200.601.001.401.002.20 SumplingperiodTs[secl 0000000000 000000000 ●■ ■●●。■●● 4327’01.234 一一 一『 Fig.2Ts-PIcharacteristic ㈹’三×昌一N]q UWbP P匹旧P 図2は,制御サンプリング周期Tjを2.1秒に固定し, 観測サンプリング周期Tsをパラメータに上記の2乗 積分値を描いたものである.同図より,No.2の実用 的サンプリング周期における最適値がNo.1の最適 値と同程度の値となり,No.2の本制御方式がNo.1 制御方式に十分代替えし得ることが明らからになった. また,図3~図6は図2のNo.1およびNo.2に対す るそれぞれの最適値Tsを用いたときのAfI,Af2, APM。およびAPCを示したものであり,同図から も明らかなように本制御方式の有効性が見られる. I DOB.O、10.DC 7Z.U【 Fig.4ResponsesofAf2.

(6)

琉球大学工学部紀要第42号,1991年 89 4.むすび 本論では,観測と制御のサンプリング周期が異なる ディアルレートサンプリング適応形負荷周波数制御器 を構築すると共に,本制御器を再熱式火力系からなる 2地域電力系統モデルに適用し,木制御方式の有効性 を検証した. 0.20 5050050 11000l ●●●の■● 000000 一’一 一【‐三xエエコロ一・【抄』勺 参考文献

〕DID.n[ [1]COMMITTEEonLOADDISPATCHINq lEE叩JAPAN:mLoadfrequencycontrol ofpowersystemUlTechnicalReport,IEE, Japan,DivⅡ’No.40(1976) [2]ELGERD,0.1.,andFOSHA,CE.:'’ 0primummegawatt-frequencycontrol ofmultiareaelectricenergySystems,,, IEEETreas.,1970,PAS-89,pp556-562 [3]CABIN,R、K,BUDGE,M、0,and RASMUSSEN,P.:',Anoptimallinear systemsapproachtoload-frequency control'0,ibid.,1971,PAS-9 [4]SHEIRAH,M,A、,MALIK,0.P.,and HOPE1G.S,:,,Minimumvariance stratBgyforload-frequencycontrolIIInt. 』・ofE1ectr・Power&EnergySsyst., 1986,8,(2),ppl20-126. [5]CHENG,SJ,,MALIK,0.P.,and HOPE,GS.:,'Self-tuningstabiliserfor amulti-machinepowersystemI,,Proc・ IEE,1986,133,(4)]pp,176-185. [6]CHANDRA,A、ⅡMALIK,OP.,and HOPE,G,S:,'Aself-tuningcontroller forthecontrolofmultimachinepower systems叩,IEEETrans.,1988,3,pp、1065 -1071.

[7]ROSTAMKOLA1,N,,PHADKE,AG.,

LONG,W,E,andTHORP,』.S:”An adaptiveoptimalcontrolstrategyfor dynamicstabilityenhancementofac/dc powersystems,',ibid.,1988,3,pp、1139-1145. [8]YAMASInTA,K,andMIYAGI,H、:', Multivariableself-tuningregulatorfor 0.15 0.20 ofAPtie. 5Responses Fig 6543210 {㈱1つ「xエエコユlU』。 ヘP )

-1i;,に

【Ⅱ |に9群一一戸Ⅵu予告…一戸一戸 0,2.004.006.0080010.0[ 2 Fig.6ResponsesofムPC.

(7)

90 デュアルレートサンプリング適応形負荷周波数制御:山下・平山・宮城 [10]M、L、Kothari,P,SSatsangi,J・Nanda :',Sampled-dataautomaticgeneration controlofinterconnectedreheatthermal systemsconsideringgenerationrate constants,IEEETransiPAS-100,2334 (1981) [11]J、Kanniah,0.P、Malik1GSHope:w Self-tuningregulatorbasedoludual-rate sampling,IEEETrans.,AC-29,755(1984) loadfrequencycontrolwithintcraction ofvoltageonloaddemand”,IEEProc., 1991,138,ppl77-l83. [9]YAMASHITA,K、,andMIYAGLH.:', Loadfrequencyself-tuningregulator forinterconnectedpowersystemswith unknowndeterministicloaddisturbanc‐ es",InLJ・ofControL1989,49,pp、1555 -1568

参照

関連したドキュメント

If condition (2) holds then no line intersects all the segments AB, BC, DE, EA (if such line exists then it also intersects the segment CD by condition (2) which is impossible due

The inclusion of the cell shedding mechanism leads to modification of the boundary conditions employed in the model of Ward and King (199910) and it will be

We proposed an additive Schwarz method based on an overlapping domain decomposition for total variation minimization.. Contrary to the existing work [10], we showed that our method

Keywords: continuous time random walk, Brownian motion, collision time, skew Young tableaux, tandem queue.. AMS 2000 Subject Classification: Primary:

Kilbas; Conditions of the existence of a classical solution of a Cauchy type problem for the diffusion equation with the Riemann-Liouville partial derivative, Differential Equations,

This paper develops a recursion formula for the conditional moments of the area under the absolute value of Brownian bridge given the local time at 0.. The method of power series

To derive a weak formulation of (1.1)–(1.8), we first assume that the functions v, p, θ and c are a classical solution of our problem. 33]) and substitute the Neumann boundary

In order to be able to apply the Cartan–K¨ ahler theorem to prove existence of solutions in the real-analytic category, one needs a stronger result than Proposition 2.3; one needs