~
ただし，密度関数
g ( c )
概形の具体的表示を要する際は， 4次までのキユ ムラントについて次のエッジワース展開を用いる。『 マ
rさらに，キュムラントによる粒度分布記述は粉砕機構 (Fig.502中ではロー ラ近傍)模擬の計算量低減に効果的である。すなわち，粉砕機構における粒子 径変化の確率は通常対数軸で整理して粉砕特性関数市
‑ x )
と称され，その 概形の一例を Fig.504に示す。(/)
﹀パ 干一 ω 心七一心の﹄CKAHO
﹂ 門 ]
‑4 ‑2
(1/2980) (1/54. 59)
ratio in logarithmic ξ‑x
F i g . 5 0 4 Example
01P u l v e r i z i n g C h a r a c t e r i s t i c s
このとき，粉砕機構通過による粉体粒度分布の変化は次式の重畳積分で模擬 され，これを従来の区間分割法で数値計算すると多大な計算量となる。
g
o p ( c )
=f s ( c ‑
x ) g山
(503)ーコ士
ここに，
g o p ( c )
，g ; p ( c )
はそれぞれ粉砕機構出口，入口の粒度分布である。87
~
一方，キュムラントを用いると上式の重畳積分は次の加算式に帰着し，著 しい計算量低減となるc
μop μ/j μ5 σ op σ 戸 σ S
一 +
K30p K3ip K35 I (504)
K40P K4ip κ45
ここに，添字 op，/pヲ S は順に粉砕機構出口，入口，粉砕特性関数を示す。 Figure 505はFig.504の特性の粉砕機構について出入口の粒度分布の関係の
(504 )式による算出例であり，簡単な演算にもかかわらず，良好な模擬結果を 示している。
一心何﹄
O
﹂丘
¥A
H一ω
むC
℃
﹀ リ
F
particle size in logarithmic
F i g . 5 0 5 E x a m p l e
01P a r t i c l e S i z e D i s t r i b u t i o n a t P u l v e r i z i n g P r o c e s s
粉砕機高速版動特性モデルは以上のとおり，計算量の低減に特徴を有し，さらに 既報[3]に述べた解析精度の検証を終えている。
88
実機試験 5.2.3
~
本サブシステムはlOOOMWe発電用ボイラにおける全7台の国内最大級の石 炭粉砕機の制御に適用され， Fig.506には状態推定ステップの動作履歴の一例 を粉砕機の負荷変化から停止操作を通じて示している。本例により粉砕機差圧，
の推定， お モータ動力の予測値と計測値の偏差に着目した石炭粉砕性(HGI)
よび粉砕性推定値に基づく粉砕機出口微粉炭流量，粒度分布の推定の状況が把 握できる。このうち粉砕性および粉砕機出口粒度分布については実機よリサン
プリングした石炭の分析値を追記してあり，本図より状態推定機能の精度を確 認できる。
F i e l d Testing 0 1 the P u l v e r i z e r C o n t r o l l i n g
89
Ishut‑downl
75 time [minJ
pulverized coal flow(estimated) (samp I ed)
Subsystem
本under 76μ(200mesh pass)
￨ l
叫 ‑swingi叩 ￨
。
60
[工
¥ド]玄O
一半
一句
00
。
﹁1 11
寸﹂1lll
ハUハUハ
U R J
[諒
]c o一
ち円
三﹂
μ ω 一 ℃ ωN
一ω
︒ 一
ο一
H L E
z; 70 寸
コ
ニ 1
〉、 1
+‑' I
4コ i
何 I
E50 J
b且
F i g . 5 0 6
，・ー
さらに Fig.507は本サブシステムによる石炭粉砕性同定精度の検証にあたり，
種々の炭種について手分析値との関係の一例を示す。本例において各炭種の HGI同定値は手分析値と 15程度の偏差で一致しており，実用上十分な精度と 評価できる。
F i g . 5 0 7 V a l i d a t i o n f o r t h e G r i n d a b i l i t y I d e n t i f i c a t i o n
90
また， Fig.508は本サブシステムによる適応状態機能の例証として，石塊扮 砕性 (HGI)の相違による粉砕機起動時の出炭量応答の相違を示す。Figure501
に示すとおり石炭粉砕性の同定値に応じた出炭量推定値が変化する状況が把
C o a l E ( H G 1 = 5 5 ) c o a l D
( H G 1 = 7 0 )
握できる。
‑・ー
Outlet Coal Flow CSimulated)
I n l e t
2 0 1 ‑ 60 40
[ 工 ¥ μ ]
﹀﹀
O 一 比
一 c
o O
1 5
E s t i m a t e d c o a l f l o w i n d i f f e r e n t c o a l s
91