第章

3

石炭ガス化炉の特性と炉内現象評価技術の開発

第３章

石炭ガス化炉の特性と炉内現象評価技術の開発 ● 目

横須賀研究所

エネルギー機械部

犬丸

淳

横須賀研究所

エネルギー機械部

原三郎

横須賀研究所

エネルギー機械部

芦沢正美

横須賀研究所

エネルギー機械部

大高

円

横須賀研究所

エネルギー機械部

渡辺裕章

横須賀研究所

エネルギー機械部

梶谷史朗

横須賀研究所

エネルギー機械部

市川和芳

横須賀研究所

エネルギー機械部

沖

裕壮

３−１

研究の背景

３−２

石炭ガス化炉の特性と炉内現象

３−３

石炭ガス化炉数値シミュレーション技術

３−４

石炭ガス化反応特性の解明

３−５

石炭ガス化炉における灰生成・付着挙動の解明

３−６

今後の展開

…………………………………………………………………………………………………………………39…………………………………………………………………………………………40………………………………………………………………………………42………………………………………………………………………………………………46………………………………………………………………………51…………………………………………………………………………………………………………………57

コラム：石炭中の鉱物について

犬丸……………………………………………………………………………………………………58

淳（10ページに掲載）

芦沢 正美（28ページに掲載）

38

次

原

三郎（10ページに掲載）

大高円（1994年入所）
これまで、石炭ガス化に関する研究に携わり、数値シミュレーション技術による石炭ガ
ス化炉内現象の解明に取り組んできました。
今後は、石炭ガス化複合発電の実用化に向け、石炭ガス化炉の設計および運転支援を目的とした炉内現象評価ツールの開発を進めて行く
予定です。

渡辺 裕章（1998年入所）
石炭ガス化炉内の伝熱流動状態の把握，ガス化炉性能評価等のための数値シミュレーシ
ョン技術の開発を行ってきた。現在、開発した数値シミュレーション技術を使った石炭ガ
ス化炉設計・運転支援ツールの開発，超重質
油ガス化炉の数値シミュレーション技術の開発等に従事している。

梶谷 史朗（1993年入所）
これまで、石炭ガス化の研究に携わり、高温高圧におけるガス化反応性を解明するとともに炭種による反応性の違いを明らかにしてきました。今後は、石炭ガス化複合発電プラントの実用化と燃料拡大技術に貢献するために、ガス化反応性の予測技術の確立に取り組んでゆきたいと思います。

市川 和芳（1991年入所）
入所後、２トン/日炉を用いた適合炭種拡大技術や灰付着性評価に関する研究に携わってきました。また、1995年より２年間、石炭ガス化複合発電技術研究組合に出向し、200
トン/日パイロットプラントによる運転試験
研究において、ガス化炉およびプラント全体
性能の評価を行いました。現在は、ガス化炉
内および熱交換器管への灰付着現象の解明およびその予測手法の開発を目指した研究に取り組んでいます。今後は、さらに、スラグ有効利用方策の拡大を目的にスラグ排出プロセスの高度化を図って行きたいと考えています。

沖裕壮（1989年入所）
これまで、微粉炭火力や石炭ガス化に関する研究に携わり、低負荷対応型微粉炭バーナの開発や、ボイラおよびガス化炉で生成する石炭灰の特性評価を行ってきました。今後は、ガス化炉生成灰の特性を予測する技術の開発
研究を通じ、石炭ガス化複合発電技術の実用化に貢献したいと思います。

３−１

研究の背景

石炭ガス化炉は、石炭ガス化複合発電（IGCC）プラ

まず、種々の実験設備を活用してガス化炉内の諸現象

ントの中核をなす技術であり、石炭ガス化炉の性能およを解明し、現象のモデル化を行う。それらを、数値シミ

び信頼性を確保することが、プロジェクトを成功に導くュレーションに導入することによって、実験結果に基づ

必要条件となる。

いた、確度の高いガス化炉内現象予測・評価技術の開発

当研究所では、これまでの２トン/日炉、200 トン/日を目標としている。また、石炭・石炭灰の性状、微細構

炉ならびに海外プラントでの運転経験から、石炭ガス化

造、形態等の基礎物性の解明およびデータベース化は、炉の実用化に向けた研究課題を整理し、当研究所が取り現象の理解や炭種による影響を解明する上で極めて重要

組む課題を設定した（第２章２節参照）。

である。図 3-1-1 中の石炭ガス化研究炉は、平成 15 年度

実証機成功の鍵を握るガス化炉の信頼性をより一層向上させるためには、スケールアップ技術、炭種による影に当研究所横須賀研究所に設置予定の研究設備である（3-6 節参照）。

響評価技術、灰に起因する諸障害の対策技術等を確立す

これにより、「IGCC 設計・運転評価ツール」の確立る必要があり、当研究所ではこれらの課題を解決するたを図り、大型ガス化炉設計に必須のスケールアップ効果

め、図 3-1-1 に示すような、)ガス化炉数値解析技術のの予測、炭種や運転条件が変化した場合の炉特性変化の開発、*実験による炉内現象の解明とモデル化、+石

予測・評価、トラブル発生時の現象解明と対策検討、コ

炭・灰基礎物性の解明、を３本柱とする支援研究を実施

ンパクト化等によるコストダウンの検討などに貢献してしている。

いきたいと考える。

溶融灰流
動モデル

灰付着
モデル

PDTF
反応性

熱交伝熱
モデル

数値シュミレーション
全体モデル：ガス反応モデル
灰付着モデル他

加圧 TG-DTA
反応性

石炭ガス化
研究炉

IGCC 設計・運転
評価ツール

レーザラマン
炭素構造分析

炉内現象の解明（実験）石炭・灰基礎物性解明

メカニズム解明、モデル化、
影響因子寄与度分析

炭素構造、微細組織
鉱物組成、灰溶融蒸発等

スラッキング装置
灰付着性実験

超高温灰融点測定装置
灰蒸発・凝縮
測定装置

CCSEM
鉱物分析

マセラル
分析装置

図3-1-1 当研究所における石炭ガス化研究の取り組み

電中研レビュー No.44 ● 39

３−２

石炭ガス化炉の特性と炉内現象

石炭ガス化炉には第１章で述べたよう様々な方式があとチャー（未燃炭素と灰分から成る微小粒子）になる。

るが、いずれの方式についても求められる機能は、極論

チャー中の未燃炭素は、さらに CO 2 や H 2O によりガスすれば、石炭を効率良くガス化し、所定の性状のガスを化されるが、完全にはガス化されずに生成ガスとともに得るとともに、石炭中の灰分を安定して系外に排出する

ガス化炉より排出される。この生成ガス中のチャーは、ことである。この機能を満足し、より高効率で信頼性のガス化炉の効率向上のため、回収・リサイクルされるが、高いガス化炉とするためには、炉内で生じている様々な

その量は炭種や運転条件により大きく変化する。従って、現象を解明することが重要である。

ガス化炉の性能予測、ガス化炉反応部の容積やチャーの

本節では、当研究所がこれまでに実施してきた石炭ガ

回収・リサイクル設備の容量の適正化等を行うためには、ス化研究で得られた知見、200 トン/日炉での経験や国ガス化炉から排出されるチャーの量を精度良く見積もる

内外で開発されているガス化炉に関する情報などに基づ必要があり、それを可能にするには、実機に相当する高き、特に重要と考えられる、石炭のガス化反応特性とガ

温・高圧下でのガス化反応特性を解明することが重要とス化炉内での灰の挙動について、その概要を述べる。ななる。

お、次項以降では、（酸素富化）空気吹き二段噴流床方式ガス化炉を念頭に置いて記すが、他方式のガス化炉に

3-2-2

ガス化炉内での灰の挙動（図 3-2-2）

ついても基本的には共通するものと考えられる。
ガス化炉内での灰の流れは基本的には以下の通りであ

石炭のガス化反応特性（図 3-2-1）

3-2-1

る。コンバスタ部に投入された石炭・リサイクルチャー
中灰分は溶融スラグとして壁面に捕捉され系外に排出さガス化炉内に投入された石炭は、熱分解により揮発分れるものと捕捉されずにリダクタ部へ運ばれるもの（キャリーオーバー灰）とに分けられる。このキャリーオー
バー灰とリダクタ部に投入された石炭中の灰は、チャー

生成ガス＋チャー粒子としてガス化炉から排出され、コンバスタ部にリサ
イクルされる。
このような灰の流れの中で、２トン/日炉や 200 トン/

リダクタ
チャー＋CO2→2CO
チャー＋H2O→CO＋H2
CO＋H2O⇔CO2＋H2
石炭→揮発分＋チャー
石炭

圧
力
容
器

リダクタ壁面
への灰付着

搬送ガス
チャー
石炭
空気

溶融スラグ

図3-2-140

溶融スラグの再飛散

コンパスタ
石炭→揮発分＋チャー
揮発分＋O2→CO2＋H2O
チャー＋O2→CO＋CO2

搬送ガス

溶融灰粒子のキャリーオーバー

搬送ガス

溶融スラグの
流動・排出スラグホール

炉内でのガス化反応

熱交換器部へのチャー付着・堆積

ガス精製設備へ

図3-2-2

灰に関する重要な炉内現象

日炉で経験したトラブル、あるいは、国内外のガス化プ差異が生じていると考えられる。さらに、リダクタ内をラントの情報から、解明すべき灰に関する重要な現象と通過する間に性状が変化すると考えられる。従って、こしては、)コンバスタ部での溶融スラグの流動・排出、

れらリダクタ内に存在する灰粒子の壁面への付着性を評*リダクタ壁面への灰付着、+ガス化炉後流の熱交換器価するためには、リダクタ内での灰粒子の性状とその変でのチャー付着・堆積、が考えられる。

化を CCSEM などの高度分析機器を活用して明らかにするとともに、壁面への付着性との定量的な関連付けを行¸

溶融スラグの流動・排出特性

う必要がある。

コンバスタ部では溶融スラグは、壁面に捕捉されコン
バスタ底部のスラグホールから排出される。ガス化炉をº

熱交換器部へのチャー付着・堆積特性

運転する上で、このスラグホールからの溶融スラグの安ガス化炉設備としては、ガス化炉本体の後流側に熱交定排出は必要不可欠なものである。スラグホールからの換器が設置され、熱回収を行う。200 トン/日炉では、溶融スラグの排出性を評価するためには、灰の溶融特性

伝熱管へのチャーの付着・堆積により、時間とともに伝の検討が必要であり、コンバスタ部を模擬した還元雰囲

熱性能が低下するため、除煤装置による伝熱管の清掃を

気下での各種石炭灰の軟化・溶融特性を明らかにするこ

定期的に行っていた。この伝熱性能の低下傾向は、炭種とが重要となる。また、スラグホール形状の影響に関すにより大きく異なっており、チャー量の多い炭種よりもる検討も必要となる。

チャー量の少ない炭種の方が、低下傾向が激しくなるケ

このようなスラグホール周りの現象に加え、200 トン

ースも確認された。このことから、伝熱性能への影響を/日炉では、一旦壁面に捕捉された溶融スラグの一部が、評価するためには、チャー量の違いだけではなく、チャコンバスタ出口部（スロート部）の強い旋回上昇流によーの性状（組成、粒径、安息角などの粉体特性等）が付り飛散し、リダクタ部壁面に付着・成長するという現象

着・堆積特性に及ぼす影響を明らかにする必要がある。

（スラギング現象）が生じ、大きな運転障害となった。
200 トン/日炉では、炉形状に変更を加え、スロート部

前項までに述べたよう、石炭ガス化炉の特性を評価すの流速を低減することにより、この現象は解決された。

るためには、各種現象を解明することが極めて重要であ

しかしながら、溶融スラグ飛散現象に関する定量的検討

る。さらに、実際のガス化炉内では、個々の現象は単独

は十分ではなく、溶融スラグの物性（粘性、表面張力等）に起こるものではなく、互いに影響を及ぼし合っており、を考慮した模擬流体によるコールド試験などから、各種一つの条件が変われば全てに影響を及ぼすため、現象解因子が飛散現象に及ぼす影響を定量的に評価し、発生条明研究で得られた結果を、実機の設計や運転に具体的に件を明らかにする必要がある。

反映していく手法として、数値シミュレーション技術の活用も必要になる。

¹

リダクタ壁面への灰付着特性

石炭ガス化炉内の各種現象の実験的解明とモデル化を

前述の溶融スラグの飛散現象が無くとも、リダクタ内行い、その結果を反映した確度の高い数値シミュレーシに存在する、コンバスタ部からのキャリーオーバー灰あョン技術を開発することにより、スケールアップの影響

るいはリダクタ投入石炭中の灰が壁面に付着し、長時間も含めたガス化炉の最適設計、安定運転条件の検討やト運転の妨げになる可能性がある。灰融点以上の高温雰囲

ラブルシューティング等の運転支援が可能となる。当研

気であるコンバスタ部を通過したキャリーオーバー灰と究所が進めている現象解明研究と数値シミュレーション

リダクタ投入石炭中の灰とでは粒子の温度履歴が大きく技術の開発について、次節以降に述べる。

異なるため、その性状（粒径・形状・含有鉱物質等）に

電中研レビュー No.44 ● 41

３−３

石炭ガス化炉数値シミュ
レーション技術ュレーションに取り込むことにより、さらに確度の高い

3-3-1

数値シミュレーション技術の必要性と位置付け

数値シミュレーション技術を開発することが可能となる。
当研究所では、この技術を応用することにより、炉スケ
ールアップ効果を含む炉設計評価および炭種・運転条件

石炭ガス化炉を設計するためには、炉内現象の予測・評価技術が不可欠である。例えば、ガス化炉の形状・寸変更、トラブルシューティング等の運転支援に向けたツール（図 3-3-1）構築を目指し、研究を進めている。

法を決定するためには、炉内に投入される微粉炭の滞留時間および反応速度を予測し、その反応性を評価する必要がある。また、石炭ガス化炉の運転においても、炉内

3-3-2

石炭ガス化炉数値シミュレーションの概要

現象の予測技術は重要であり、炉内現象を把握すること
により、炭種に応じた適切な運転条件を設定することが可能となる。
炉内現象予測・評価技術のひとつとして、計算機によ

ガス化炉内現象は、流動、伝熱、反応の３つに大別される。それぞれの現象に対し、本数値シミュレーションで採用した主な解析手法と解析モデルを以下に示す。

る数値シミュレーションが挙げられる。数値シミュレー

まず流動であるが、気相は弱圧縮性乱流とし、乱流モションの長所は、実験に比べ条件変更が容易な点にある。

デルには標準的な k -εモデル ºを採用した。粒子輸送は、

例えば、実際の炉において炉形状やバーナ配置等を変更

粒子群の移動を代表粒子によって模擬する Lagrangian

するためには、長い工事期間と多額の費用が必要となる

モデル»によって解いた。気相と粒子の相互作用は粒子が、数値ミュレーションでは、解析格子や入力条件を変に作用する流体力の反作用を気相流動において考慮し、更するだけである。石炭供給量や炭種等の条件変更も、粒子間の衝突は無視した。壁面に衝突した粒子は、粒子

運転条件や石炭分析値に基づいて入力条件を変更するこの温度履歴と液相率から壁面への付着を判定する灰付着とにより、容易に行うことができる。

モデルにより取り扱った。

石炭ガス化炉を対象とした数値シミュレーションでは、

伝熱については、乱流熱拡散を含む対流伝熱、粒子散

ガス流動、粒子輸送、伝熱、化学反応等の様々な炉内現

乱を伴う輻射伝熱および化学反応に伴う反応熱を考慮し象を考慮するため、多くの解析モデルが導入される。解た。輻射伝熱は Discrete Transfer 法¼によって解かれ、析モデルの導入には、目的とする炉内現象の実験に基づ気相および粒子相の輻射物性は、ガス組成分布や粒子濃くモデルパラメータの設定が必要となる。また、数値シ度分布によらず一定とした。粒子温度は、気相と粒子間

ミュレーション結果の妥当性を確認するには、比較のたの対流伝熱および後述のチャーガス化反応に伴う反応熱めの試験データが不可欠である。さらに、数値シミュレのバランスから求められる。

ーション結果から炉内現象を評価するためには、石炭ガス化炉に対する多くの知見が必要とされる。

気相反応は、シフト反応を含む５つの総括反応を考慮した。反応速度は渦消散モデル½による乱流混合速度と当研究所は、２トン/日試験炉および 200 トン/日パイ

化学反応速度¾¿Àから求めている。粒子反応は、熱分解

ロット炉の運転を通じ、多くの知見を取得し、膨大な試反応とチャーガス化反応を考慮した。後者のモデル化に

験データを保有している。また、高温高圧下における微不可欠な高温高圧下のチャーガス化反応速度には、当研

粉炭の反応性¸、石炭中灰分の壁面付着性¹といった個々

究所の PDTF と TG による実験データを用いた梶谷らのの炉内現象解明を目的とする実験装置を開発し、研究を提案する実験式¸を採用した。

行っている。これらの石炭ガス化炉に対する知見や試験以上を取りまとめ、図 3-3-2 に示す。なお、解析手法およ

データをベースに、実験に基づく解析モデルを数値シミ

び解析モデルの詳細については文献Áを参照して頂きたい。

42

入力データ・炉設計データ（炉寸法、バーナ配置）

石炭ガス化炉数値解析用
並列計算機

出力データ
・流速分布
・濃度分布

・温度分布
・粒子軌跡

・石炭分析値（石炭性状、粒径分布）
・石炭反応特性（チャーガス化反応速度）

ガス化性能予測・灰付着特性・炭素転換率・生成ガス発熱量
・生成チャー量 ・冷ガス効率

・運転条件（空気比、石炭供給量）
・etc.
炉設計・運転条件指標
・スケールアップの影響予測

・炭種・運転条件変更時の炉特性変化予測
・炉内現象再現によるトラブルシュウティング・コンパクト化・コストダウンの可能性検討

図3-3-1

炉設計・運転支援ツールの概念図

ガス化炉内現象数値解析技術の概要
流 れ 場：三次元弱圧縮性乱流
解 析 手 法：Finite Volume Method
Hybrid Upwind Differencing Method
SIMPLEC Algorithm
乱流モデル：k-ε２方程式モデル
固気二相流：Eulerian-Lagrangian Method
輻 射 伝 熱：Discrete Transfer Method
熱分解反応：CmHnOl →（CH4、H2、CO、CO2、H2O）
気 相 反 応：CH4 ＋ 1/2O2 → CO＋2H2
H2 ＋ 1/2O2 → H2O
CO ＋ 1/2O2 → CO2
CH4 ＋ H2O →← CO ＋ 3H2
CO ＋ H2O →← CO2 ＋ H2γチャーガス化反応：C ＋（1−―）O
2 →γCO ＋（1−γ）CO2
2
C ＋ H2O → CO ＋ H2
C ＋ CO2 → 2CO
解析格子：Multi Block、Body Fitted Coordinates

図3-3-2

炉内現象モデル様な傾向を示しており、良い一致を示している。図 3-3-

3-3-3

数値シミュレーション結果の例

5 に炉中心軸を含む縦断面におけるガス温度分布および
H 2、CO、CO 2、H 2O の質量濃度分布を示す。リダクタ

数値シミュレーションによって明らかとなった炉内現象の例を以下に紹介する。

バーナの後流でガス温度が低下し、CO2 と H2O が消費さ
れ H2 と CO の濃度が上昇しているのが分かる。これは熱最初に示すのは、２トン/日炉に対する数値シミュレ

分解反応を終えた微粉炭粒子（チャー粒子）のガス化反ーション結果Áである。参考のため、図 3-3-3 に２トン/

応（吸熱反応）によるものと考えられる。図 3-3-6 と図

日炉の解析格子を示す。図 3-3-4 は、炉中心軸上のガス

3-3-7 は、炉内空気比（注１）をパラメータとした場合の炉温度分布である。数値シミュレーション結果は実験値に比べ多少低めの温度分布となっているが、両者はほぼ同（注１）：ガス化炉投入空気量/石炭およびチャーの理論燃焼空気量

電中研レビュー No.44 ● 43

2000

0.015

500
ガス温度 K

図3-3-3

0.00
CO2濃度

90
λ0=0.469｝
λ0=0.495｝
λ0=0.525
｝

1500

1000

炉内炭素転換率 ηp ％

解析結果
実験結果
解析結果
実験結果
解析結果
実験結果

80
70
60
国内A炭（解析結果）
国内A炭（実験結果）
豪州B炭（解析結果）
豪州B炭（実験結果）

50
40

0.0

図3-3-41.0
2.0
ガス化炉高さ

0.00
H2O濃度

100

RB

2000

500

0.05

ガス化炉内ガス組成分布（質量分率）

2500
中心軸上ガス温度Tg K

0.20

0.00
CO濃度

図3-3-5

CB：コンバスタバーナ
RB：リダクタバーナ
CB

0.000
H2濃度

解析格子

0.35

3.0
Hm

4.0

5.0

ガス化炉中心軸上ガス温度分布

内炭素転換率（注２）と生成チャー量を示したものである。

30
0.3

0.4
炉内空気比 λp

図3-3-60.5

0.6

炉内炭素転換率

次に、200 トン/日クラスの石炭ガス化炉において、炉内炭素転換率は、炉内空気比の増加に伴い高くなるこ

炉形状およびバーナ配置変更を行った場合の数値シミュとが、実験から明らかにされている。数値シミュレーシ

レーション結果Âを示す。図 3-3-8 は、炉中心軸を含む

ョンは、その定性的な傾向を良く再現していることが分

縦断面の流速分布を比較したものである。炉形状およびかる。また、数値シミュレーションは、国内Ａ炭に比べバーナ配置の変更により、スロート部の中心に弱い下降

豪州Ｂ炭の方が低い炉内炭素転換率を示す実験結果の傾流が発生している。図 3-3-9 は、炉壁面の剪断応力（垂向も再現している。同様に、生成チャー量についても、

直成分）を比較したものである。剪断応力は壁面を流下

炉内空気比に対する傾向および炭種による傾向の違いが、する溶融スラグを飛散させる原因のひとつと考えられる。

数値シミュレーションにより良く捉えられている。

炉形状およびバーナ配置の変更により、溶融スラグの飛散が懸念されるスロート部壁面での剪断応力は、約 1/4（注２）：生成ガス中炭素量/石炭およびチャー中炭素量44に低減されているのがわかる。

0.0
壁面せん断応力 Pa

60
kg/h

コンバスタ
バーナ

国内A炭（解析結果）
国内A炭（実験結果）
豪州B炭（解析結果）
豪州B炭（実験結果）

70

生成チャー量炉形状変更前（スロート径小）

4.0

80

50
40
30
20

リダクタバーナ
コンバスタ スロート

リダクタ

−4.0
炉形状変更後（スロート径大）4.0 コンバスタ
バーナ
0.0

10

リダクタバーナ
00.30.4
炉内空気比 λp0.5

0.6

コンバスタ スロート
0.0

図3-3-72.0
4.0
コンバスタ底面からの無次元距離

生成チャー量

図3-3-9

リダクタバータ

6.0

壁面せん断応力分布

リダクタバータ

下降流
なし

下降流
あり

D2

D2

チャーバーナD1

コールバーナ

リダクタ

−4.0

D1

D2/D1＝0.4

図3-3-8

D2/D1＝0.8

チャーバーナコールバーナ

流速ベクトル分布

良等により、数値シミュレーションの予測精度が向上す

3-3-4 ま と め

れば、定量的な炉内現象の予測も可能と思われる。
今後は、平成 15 年度に設置予定の石炭ガス化研究炉

本数値シミュレーション技術を活用することにより、を活用し、数値シミュレーションの予測精度向上に向け様々な条件での炉内現象を定性的に比較することが可能

た解析モデルの改良を行う。さらに、石炭ガス化炉の長となる。現状において、炉内現象の定量的な予測は難し

期連続運転において、非常に重要な課題である溶融スラいと思われるが、定性的な傾向を把握することにより、グの炉外への安定排出に関し、数値シミュレーションに

石炭ガス化炉の設計支援、運転条件の最適化、トラブル

より、相変化を伴う溶融スラグの流動伝熱現象を予測し、発生時の原因解明が期待される。また、解析モデルの改溶融スラグ排出性やスラグホール形状の検討を行う。

電中研レビュー No.44 ● 45

３−４

石炭ガス化反応特性の解明
発分や一部のチャーは燃焼し、その燃焼熱を使ってチャ

3-4-1 石炭ガス化炉内における反応の概要

ーのガス化反応が起こる。チャーとは石炭が熱分解されて残った固定炭素分と灰分から成る固体である。チャー

石炭の反応性は炭種により大きく異なるため、ガス化

は酸素や二酸化炭素、水蒸気によってガス化され、これ炉の設計や運転条件の決定、性能予測や性能評価には石らのガスをガス化剤と呼ぶ。なお、空気吹き二段噴流床

炭の反応性に関する知見が必要不可欠である。例えば、ガス化炉内の各部での反応は図 3-2-1 で示した通りであ

空気吹き二段噴流床ガス化炉では、多くの未燃炭素（チる。チャーのガス化は燃焼よりも著しく反応速度が遅く、ャー）が炉出口から排出されるのでそのチャーをリサイ

また、炭種などによって反応性が大きく異なるため、ガクルする。そこで、チャー生成量を正確に予測すること

ス化炉の設計にはチャーのガス化反応速度が必要である。

がチャー系容量の最適設計につながり、また、負荷や空しかし、噴流床ガス化炉の場合では、流動層ガス化炉な気比などの運転条件や炭種を変更したときの生成チャー

どとは異なり炉内が高温になるため、ガス化反応は迅速量の変化を予測することがガス化炉の安定運転のためにに進むがガス化剤の拡散の影響が大きくなり、これまで必要である。

に豊富に行われている低温で大気圧下でのガス化反応性

噴流床ガス化炉内での主要な反応を表 3-4-1 にまとめの研究結果をそのまま使うことができない。噴流床ガスる。高温の炉内に投入された石炭は急速熱分解され、揮化炉はガスタービンと組み合わせるために炉内圧力が

表3-4-1

石炭急速熱分解

石炭ガス化炉内で起こる主な反応

揮発分の燃焼チャーの燃焼チャー

石炭

石炭熱分解反応
石炭→揮発成分＋チャー

気相平衡反応チャーガス化反応

灰

石炭を加熱するとガスやタール分が放出される。高温炉内での急速熱分解は数百ミリ秒で完了する迅速な反応である。

燃焼反応
揮発成分＋aO2→b CO2＋c H2O
C＋O2→CO2
2C＋O2→2CO（O2ガス化）

ガス化炉では理論量よりも少ない空気や酸素が投入される。燃焼反応で生じた熱と二酸化炭素や水蒸気がチャーのガス化反応に使われる。

2CO＋O2→2CO2
チャーガス化反応

チャーガス化反応は吸熱反応なので、ガス化の進行とともに炉内温度は低下する。空気吹き２

C＋CO2→2CO（O2ガス化）段噴流床炉ではコンバスタ出口の高温の燃焼ガスにリダクタ石炭が投入されるため、リダクタでC＋H2O→H2＋CO（水性ガス化）の急速熱分解とCO2ガス化が支配的となる。また、スラリーフィード炉などでは水性ガス化の効C＋2H2→CH4（水素化反応）

果が大きくなる。

気相反応→CO2＋H2（シフト反応）
CO＋H2O←
CH4＋H2O←→CO＋3H2（メタン改質）

46

高温炉内では迅速にシフト平衡へ達する。高温ほど一酸化炭素は増加し、メタンは減少する。

2.5 〜 3MPa 程度で運転され、炉内温度は 1800 ℃程度ま高温高圧下での実験が可能である。また、炉底から挿入で達することがある。そこで、当研究所では 1800 ℃、したサンプリングプローブを上下にトラバースすること

2.5MPa（25 気圧）までの実験が可能な世界的にもトッ

によって滞留時間を制御し、炉内の流動はプラグフロー

プクラスの仕様をもつ超高温・加圧型燃料反応実験設備と仮定して解析を容易にした。実験にあたっては、まず、（PDTF）を開発し、これまでにデータがほとんどない

常圧型 DTF（ADTF）を用いて窒素ガス中で微粉炭を高温高圧下での石炭ガス化反応性の研究を行っている。

熱分解してチャーを調整し、次に、PDTF を用いて窒素希釈したガス化剤中に再びチャーを投入してガス化実験

3-4-2

石炭の反応実験を行い、サンプリングプローブによって反応凍結して捕集したガスやチャーの組成分析や性状分析から反応解析

石炭の反応解析には熱天秤（TG）が一般によく用いを行った。

られている。TG は天秤の試料皿の部分に電気炉を備え、試料を加熱しながら重量変化を測る装置である。容易で3-4-3

ガス化反応速度解析

かつ高精度な測定が可能だが、高温ガス化のような速度の速い反応実験には適さない。一方、DTF（Drop

チャーのガス化は気固反応なので、反応の進行ととも

Tube Furnace）は縦型の管状電気炉で、炉頂から固体にチャー粒子の反応界面が変化して反応速度が変わるた

試料を落とし、炉内で反応した生成ガスや固体を分析しめ、粒子の反応をモデル化する必要がある。基本的なモて反応解析を行う装置である。石炭やチャーがガス中にデルとして、図 3-4-2 _の容積モデル、グレインモデル、分散される気流層反応炉なので、実炉と同等な反応場で細孔モデルなどが広く用いられている。容積モデルは、の高温実験が可能である。図 3-4-1 に示す当研究所の粒子内で均一に反応が起こるとしたもので、グレインモ

PDTF は加圧型の DTF で、実機のガス化炉に相当する

デルは一つの粒子はさらに小さい球形微粒子（グレイン）

種々の燃料
の供給

ADTF ADTF

反応ゾーン
0〜1200mm

電気ヒータ

PDTF
PDTF

広範囲な実験条件●温度：Max. 1800 ℃●圧力：Max. 2.5 MPa●各種反応ガス
サンプリングブローブ●任意の反応時間

反応過程の生成ガス・粒子の捕集と性状解析

図3-4-1

超高温・加圧型燃料反応実験設備（PDTF）の概要

電中研レビュー No.44 ● 47

灰分層

ガス境膜

容積モデル
外部
拡散律速

グレインモデル

CAb
CAs
CAm

粒子内
拡散律速

反応律速

log（反応速度）

ガス化剤濃度

反応領域

細孔モデル

r0

高温

0

半径位置

_

粒子反応モデル

図3-4-2`

1／T（粒子温度）

高温における反応速度の変化

石炭（チャー）の反応の概念図の集合体と考え、ぞれぞれのグレインの表面で反応が起タ示す ¸Ä。ガス化反応は全圧よりもガス化剤分圧の影こるとする。細孔モデルでは粒子内部には多くの細孔が響が強いことが分かり、その圧力影響（反応次数 n）や存在し、細孔の内表面で反応が起こるとするÃ。これら

温度の影響（活性化エネルギー E）などを明らかにした。

を検討したところ、チャーのガス化反応には細孔モデル

なお、図 3-4-3 _のように酸素による燃焼・ガス化は迅が適しており、細孔モデルに使われる細孔構造係数によ速な反応であるので、反応速度の遅い二酸化炭素と水蒸って様々な炭種の構造の違いを表現することが可能であ気によるガス化が律速となり、ガス化炉設計には酸素よる ¸Ä。しかし、反応速度の変化と比表面積の変化が必

りも二酸化炭素や水蒸気によるガス化反応速度が重要でずしも一致せず、また、細孔構造係数を物理的な分析かある。以上の速度解析の結果、二酸化炭素と水蒸気につら求めることができないため、細孔モデルにも改良の余

いて同じ温度、ガス化剤分圧において比較すると、水蒸地があると考えらえる。

気によるガス化反応速度は二酸化炭素よりも５倍以上早高温炉内での反応は、上記モデルによる活性点での化いことが分かった。また、炭種による反応性を比較する

学反応だけでなく、ガス化剤の拡散速度が影響を与えると、炭種によりチャーのガス化反応速度が 10 倍近く違ようになる。図 3-4-2 `のように、高温になると化学反

うこと、高温において拡散律速領域へ遷りやすいものと応速度は加速するが、チャー粒子の中をガス化剤が拡散遷りにくいもののあることが分かった。反応律速領域でする速度が化学反応速度と同程度かむしろ遅くなるため、

E や n の値に対して、粒子内拡散律速領域では見かけの化学反応速度が粒子の反応速度を決める領域から粒子内活性化エネルギーが約 E/2、見かけの反応次数が（1+n）

拡散速度が律速となる領域に遷り、反応速度の温度依存/2 となることが解析的に示されており、本結果はほぼ性が小さくなる。さらなる高温では、粒子の周囲（バルこれに従っている。前節で述べたように、炉内数値解析

ク）から粒子外表面へのガス化剤の供給が最も遅い外部において本反応速度データをそのまま導入することで精拡散律速と呼ばれる領域になり、反応速度の温度依存性

度のよい計算結果を得ることができた。

がほとんどなくなる。
PDTF と TG によるチャーのガス化実験結果を解析し、

3-4-4ガス化反応性に及ぼす影響因子

低温から高温までの反応速度を求めた。図 3-4-3 _には
各ガス化剤に対するガス化反応速度を、図 3-4-3 `には

上述のようにチャーの反応速度は炭種により大きく異二酸化炭素ガス化反応速度の炭種による違いを、表 3-なる。このようにチャーの反応性に影響を与える主な因

4-2 には細孔モデルに基づく反応速度式と速度パラメー子として、炭素構造、細孔構造、ミネラル成分の触媒効481600℃

102

1200℃

1000℃

800℃

600℃
A炭チャーのガス化

101dx/dt｜x＝0/PAn（s−1）

400℃

PDTF 0.5MPa
H2Oガス化100

TG 0.1MPa

10−1
10−2
CO2ガス化
10−3

O2ガス化

10−4
10−5

5

6

7

_
1600℃

102

8

9
10
11
1/T（K−1×10−4）

13

14

15

各種ガス化剤による反応速度の違い

1400℃

1200℃

1000℃

800℃

石炭チャーの二酸化炭素ガス化

101dx/dt｜x＝0/PAn（s−1）

12

PDTF 0.5MPa
CO2分圧PA＝0.2MPa
TG 0.1MPa
CO2分圧PA＝0.1MPa

100
10−1
粒子内拡散律速10−2

反応律速A炭チャー（高温時n＝0.73）（低温時n＝0.54）
B炭チャー（n＝0.49）

10−3
10−4

C炭チャー（n＝0.44）

10−5
5

6

7

8

9

10

1/T（K−1×10−4）`

図3-4-3炭種による反応速度の違い

チャーの初期ガス化反応速度のアレニウスプロット

表3-4-2

チャーのガス化反応速度式

細孔モデルに基づく反応速度式

x：反応率［−］Edx−＝A 0・PAn・e RT ・（1−x）・√１−Ψ・ln（1−x）
dt

炭

種

気体定数 R＝8.314×10−3kJ/mol K
A 炭チャー

ガス化剤

CO2

H2 O

高温
（1200℃以上）
粒子内拡散律速領域

低温（1200℃以下）
反応律速領域

高温

細孔構造係数 Ψ33反応次数 n

0.73

0.54

適用範囲

活性化エネルギー E
頻度因子 A0t：時間［s］

PA：ガス化剤分圧［MPa］、T：温度［K］

O2

B 炭チャー

C 炭チャー

CO2

CO2

低温

高温〜低温

高温〜低温

反応律速領域

反応律速領域

反応律速領域

3

14

0.1

3

0.86

0.68

0.49

0.44

163 kJ/mol

283 kJ/mol

214 kJ/mol

130 kJ/mol

261 kJ/mol

4

9

7

6

9

6.78×10

1.09×102.45×10

1.36×10

1.23×10266 kJ/mol
2.64×108

供試チャーは常圧DTFにより窒素気流中1400℃で微粉炭で急速熱分解して調整した。
平均粒径はA炭チャーが43μm、B炭チャーが44μm、C炭チャーが20μmであった。

電中研レビュー No.44 ● 49

果などが考えられる。まず、炭素構造に関してはチャー

度は高温での粒子内拡散律速領域の現れ方に影響するとの熱履歴との関係が深いÅ。図 3-4-4 _のように、微粉思われる。図 3-4-3（b）の A 炭チャーは C 炭チャーより

炭を異なる温度で熱分解して生成したチャーに対して、も複雑な細孔構造をもつことが分かっている。また、粒二酸化炭素によるガス化反応速度を同一条件で測定した子内拡散は粒径の影響も受けると考えられる。

ところ、高温で熱分解されたチャーほど反応性が悪く、さらに、触媒効果をもつミネラル成分としてはイオン

また、PDTF により急速加熱されたチャーよりも低速加交換性のカルシウムとナトリウムがよく知られている。

熱で調整したチャーほど反応性が悪かった。この原因は、

図 3-4-3 `の B 炭にはカルシウムの含有量が多く、その図 3-4-4 `に示すように炭素の結晶構造の成長度合（黒触媒効果が現れていると考えられる。
以上のような因子がガス化反応へ影響をおよぼすこと

鉛化度）に起因することがわかった。つまり、黒鉛化のが分かったが、未だその定量的な把握がなされていない。

進んだチャーほどガス化しにくい。
また、チャー粒子表面には、顕微鏡では見ることのでそこで、われわれはそれぞれの因子ごとに定量的に検討きないマイクロ孔（孔径２ nm 以下）と呼ばれる微小なし、高温高圧下での反応速度を予測する手法を確立して細孔構造がよく発達している。この細孔構造の発達の程

ゆきたい。

d
d

Lc

大

黒鉛化度

小

0.002
0.002

熱天秤によるチャーのCO2ガス化速度0.001
ガス化条件：1100℃、常圧
チャー生成条件
PDTF生成チャー（0.15〜2.5MPa）
TG生成チャー
0
1000

1400

_

熱分解温度とチャーの反応性

図3-4-40.001

ガス化条件：1100℃、常圧
チャー生成条件
PDFT生成チャー（1400〜1700℃、0.15〜2.5MPa）
TG生成チャー
（1100〜1400℃、常圧〜0.6MPa）（常圧〜0.9MPa）

チャー生成温度（℃）

50

ガス化初期反応速度dX/dt（s−1）

ガス化初期反応速度dX/dt（s−1）

熱天秤によるチャーのCO2ガス化速度

1800

0
0.2

0.4

0.6

ラマンスペクトルV/G（ー）`炭素構造と反応性の相関

石炭の熱分解温度および炭素構造がチャーのガス化反応性へ及ぼす影響

0.8

Lc

３−５

石炭ガス化炉における灰生成・付着挙動の解明
解明およびモデル化を行い、それらとガス化炉内数値解

3-5-1

ガス化炉内灰付着・成長現象

析コードを組み合わせた灰付着性予測ツールの開発を進めている。図 3-5-2 にガス化炉内灰付着性予測ツールのガス化炉を運転する上で特に問題となるのは、チャーフローを示す。本ツールでは、まず、CCSEM

および灰粒子がガス化炉内壁や熱交換器に付着・堆積す（Computer Controlled SEM）により微粉炭中鉱物の分析ることによる諸障害である。これらは、伝熱阻害を引きを行い、この結果からチャー生成モデルを用いて、チャ起こすばかりでなく、運転継続をも困難に至らしめるこー中鉱物の粒径、組成および含有形態（可燃成分に内包Æとがある 。微粉炭ボイラ内の灰挙動については国内外された Included mineral か反応に伴い可燃成分の外部にで多くの研究例があるが、ガス化炉内については未解明放出された Excluded mineral かの区別、図 3-5-1）を予測な点が多く、灰付着性評価手法が十分確立されているとする。次にこれら生成粒子性状を初期条件として数値解は言い難い。ガス化炉内の灰付着・成長現象は、図 3-5-

析を実施し、粒子の壁面衝突時に灰付着モデルにて付着/

1 に示すような３つの現象に分類できる。炉内に投入さ反発を判別し、炉内位置毎の付着量を求める。さらに、れた石炭中鉱物質（石炭から独立に存在する鉱物質も含焼結モデルにて付着層の構造、有効熱伝導率を予測するむ）は、まず、熱分解および燃焼・ガス化反応の進行過程で、石炭から放出され、溶融・未溶融鉱物粒子、未反応炭素が生成される（灰生成現象）。次にこれら粒子は、ガス流れによる慣性力あるいは熱泳動力によりリダクタ
壁面へ衝突し、付着あるいは反発する（灰付着現象）。さらに付着粒子同士が時間の経過とともに焼結し、構造が緻密化する（灰焼結現象）。当研究所では、各基礎特性の

反発粒子

付着層成長
衝突

チャー
溶融鉱物

微粉炭中の鉱物分析（CCSEM）
・粒径分布
・鉱物組成
・含有形態
（Included/Excluded）

未溶融鉱物

粒子間の焼結緻密化

チャー生成モデル・鉱物粒径分布
・鉱物組成
・含有形態

灰中液相率の推算（平衡計算、DTAデータ）③灰焼結現象②灰付着現象INPUT
・石炭性状・ガス化炉運転条件（空気比、給炭量比等）

ガス化炉内現象予測
基本数値解析コード
・炉内ガス温度分布・粒子温度・粒子衝突位置

衝突

反発

灰付着モデル
液相率による
灰付着判別
付着

再凝縮粒子
アルカリ等の蒸発
Na K
S

灰焼結モデル
・構造強度
・有効熱伝導率

未反応炭素
燃焼・ガス化反応①灰生成現象OUTPUT◎灰付着位置◎灰付着量
・付着層構造・付着層有効熱伝導率

揮発分放出

石炭
内部鉱物質
（Included Mineral）

図3-5-1

外部鉱物質（Excluded Mineral）

石炭ガス化炉内灰挙動

図3-5-2

灰付着性予測ツールのフロー図

電中研レビュー No.44 ● 51

ものである。本節では、これらの内、CCSEM を用いた

Si

灰生成現象解明に向けた研究と灰付着モデルならびにそれをガス化炉内数値解析コードに導入した灰付着解析結果について紹介する。

3-5-2
¸

試料中の灰平均組織（推定融点1500℃以上）

ガス化生成灰特性の予測Ç

リダクタ部で生成する灰の特性把握
前述したように石炭ガス化炉では、リダクタ壁面へのCa灰の付着やガス化炉後流へのチャー（灰を含む）の堆積¸

Al

従来の分析法による三相面が問題となる。コンバスタ生成灰の大半は壁面で捕集さSiれスラグとして排出されるため、リダクタから後流における灰の付着、堆積特性を評価するには、リダクタ生成

融点1200℃以下の低融点鉱物粒子

灰の諸特性（特に溶融状況）に対する炭種や運転条件の寄与度を把握する必要がある。しかし、リダクタ生成灰の特性を検討した例はなく、ガス化炉捕集サンプルの詳細な分析が求められている。

¹

石炭灰の特性を正確に予測するために石炭灰の炉内への付着特性は、従来も微粉炭ボイラの

Ca

Al

¹

伝熱性能確保等のために検討されており、種々の付着特

CCSEM による三相面性指標が提案されてきた。しかし、実炉における灰の付
融点1200℃以下の低融点鉱物粒子着性と必ずしも一致した傾向になく、高い信頼性を持っ30た指標でないことが指摘されている ÈÉ。これは従来の該当粒子数（−）

指標が灰分の平均組成から算出されるためと考えられる。
もともと石炭は、植物が長時間かけて炭化した、いわば植物の化石であり、灰分は付近の土砂の混入などに起因する鉱物粒子である。従って、産炭地の土質の影響を強

20

Al
10

く受け、生成灰の特性も炭種によって大きく異なる。

Si

様々な鉱物粒子の混合体の特性を、平均組成だけから評0価しようとしても限界があり、鉱物的観点からの解析がCa有効である。付着特性を例に取ると、平均的な融点の高º

CCSEM で作成できる立体三相図

低よりも、リダクタ部の温度域で溶融する低融点鉱物の
図3-5-3粒子量の方が指標として適するものと考えられる。

CCSEM 分析と従来法の比較

CCSEM は、走査電子顕微鏡をベースとした自動分析装置で、数千点の灰粒子を１つ１つ分析し、その粒径、組成１点のデータしかえられなかった従来の分析（図 3-

形状(溶融状況)、鉱物組成(融点)、含有形態（図 3-5-1、5-3 ¸）と比べ、低融点の鉱物粒子がどの程度含まれてIncluded Mineral か Excluded Mineral か）等のデータをいるかを詳細に把握できることがわかる（図 3-5-3 ¹、º）。

収集し、統計的解析によりサンプル中の灰の特徴を評価することができる。図 3-5-3 は灰の溶融性などの検討に用いられる三相図であるが、CCSEM を用いれば、平均52º

CCSEM による解析事例
図 3-5-4 は、粒子の軌跡などに影響する生成灰粒径のなかった FS は、チャー中には 2.2 〜 10 μｍ程度の比較

ふるい下（wt％）

100

的小さな粒子として検出された。このことから、チャー80

中灰粒子微細化の一因は、石炭中の鉱物粒子が他の鉱物

60と接触、反応して溶融することによる微細なシリケート

40
微粉炭
チャー
（空気比0.50）
チャー（空気比0.55）

20
0

1−2.2

2.2−4.6 4.6−10

10−22

22−46

46−100粒径（μm）

全灰粒子中FE比率（wt％）

当研究所 PDTF チャーの XRD 解析結果などから、
Excluded Mineral は高温ガスに直接さらされるが、可燃成分（炭素分）に内包された Included Mineral は、吸¸ 微粉炭、チャー中の鉱物粒子の粒径

熱反応が支配的なガス化反応場では高温ガスの熱的影響

5
ワークワースが比較的小さくなると考えられる。従って低融点鉱物で

4
微粉炭
チャー
（空気比0.50）
チャー（空気比0.55）

3
2

も Included Mineral の形態であれば高温下でも溶融しない可能性があり、灰の含有形態を把握することは、灰付

1着特性を評価する上で重要といえる。当研究所は

0
1−2.2

2.2−4.6 4.6−10

10−22

22−46

46−100CCSEM により、空気比が高いほどガス化炉出口における全灰粒子中の Excluded 率（wt ％）が増大することを粒径（μm）¹
全灰粒子中FS比率（wt％）

粒子の形成であると考えられる。

FE鉱物の粒径別含有量明らかとした（図 3-5-5）。

5
ワークワース
4
微粉炭
チャー（空気比0.50）
チャー（空気比0.55）

32

また、CCSEM では粒子の形状認識も可能である。前述したチャーサンプル中に含まれる Excluded 粒子の一部は球状となっており、これはコンバスタで高温に曝され溶融し、リダクタへキャリオーバしたものと推察され

1

る。CCSEM 分析結果とチャー生成量・性状から求めた、

0
1−2.2

2.2−4.6 4.6−10

10−22

22−46

46−100粒径（μm）º

図3-5-4

FS鉱物の粒径別含有量

微粉炭、チャー中の鉱物粒径分布

ガス化炉出口における球状 Excluded 粒子量は、炉内温度が高くキャリオーバ灰量が増える高空気比条件ほど増大することがわかった（図 3-5-6）。このような溶融粒
子は高い付着性を持つと考えられ、上記の情報はガス化
炉内の灰付着モデルを検討する上で極めて有効である。

変化を考察した例である。図 3-5-4 ¸はガス化炉出口捕今後は、炭種や運転条件を様々に変化させたデータを集チャーとリダクタ投入微粉炭中の灰粒子の粒径分布を蓄積し、ガス化によって生成する灰粒子の特性を予測す比較したもので、ガス化反応により灰粒子の粒径が微細化すること、空気比の影響はあまり大きくないことがわ
80

次に灰粒径が微細化した原因について、鉱物的観点から考察した。CCSEM 分析では、一般に石炭中の鉱物質を三十数種の鉱物に分類するが、このうち Pyrite（黄鉄鉱 FeS2）などもともと石炭中に含まれる鉄系鉱物（FE）と、これら鉄系鉱物が高温雰囲気下で反応して生成すると考えられる鉄シリケート系鉱物（FS）に注目し、その粒径分布を整理した（図 3-5-4 ¹､º）。
その結果、微粉炭中の FE は 10 〜 46 μｍ程度の比較

全灰粒子中Excluded率（wt％）

かる。

60

40
モーラ
太平洋
神木
イリノイ＃5
ワークワース

20

0
0.4

0.5

0.6

空気比（−）

的大きな粒子であったが、、チャー中にはほとんどみられなかった。一方、反応前の微粉炭にはほとんどみられ

図3-5-5空気比と全灰粒子中Excluded率の関係

電中研レビュー No.44 ● 53

1500 ℃（1773K）以上の高温であるが、リダクタ部に投入

球状Excluded粒子量（kg/h）

7された石炭の熱分解・ガス化反応により、急速に低下し、

6

ガス化炉出口では約 1000 ℃（1273K）程度となる。この
5

ため、ガス化炉内には、明確に温度履歴が異なる、多様4

な溶融状態の灰粒子が存在する。従来、灰付着性評価に
モーラ

3
2

太平洋

は、灰融点や塩基度などから推測する手法がとられてい神木るが、これらの断片的な情報では、急激に変化し、かつイリノイ
1

連続した温度分布を持つガス化炉内全域への灰付着性評

ワークワース

価は困難である。また、粘性による推算も提案されてい

0
0.4

0.5
空気比（−）

0.6

るがÌ〜Î、ガス化炉のリダクタ領域は灰粘性の測定が困難
な灰軟化点以下であるため適用できない。本研究では、

図3-5-6

空気比と球状Excluded粒子量の関係

温度変化に伴う灰溶融性変化に着目し、灰軟化点以下の領域を対象に灰付着性と灰溶融性の相関について調べ、るモデルの構築を目指す。

炭種間で生じる付着性の相違について検討を行った。

3-5-3

¹

灰付着判別モデル

灰溶融性指標灰の溶融開始から流動までの連続した溶融状態変化を¸

研究のアプローチ

定量的に示す指標として、灰中に生成した溶融相の割合、

２トン/日炉および 200 トン/日炉での経験から、炉内のすなわち灰中液相率を算出した。算出手法としては、平灰付着状況は、運転条件や炭種により様々に変化するこ

衡計算、示差走査熱量分析（DSC）、画像解析などによとがわかっているÊË。すなわち、灰付着現象には、灰のる試みもあるが、当研究所では、比較的簡便であること

溶融特性と炉内の温度場および流れ場が強く影響を及ぼから、超高温示差熱分析（DTA）Ï より算出する手法をすと言える。図 3-5-7 にガス化炉内の温度分布を示す。炉採用した。図 3-5-8 に示差熱分析結果の一例を示す。こ内ガス温度は、コンバスタ部では灰を溶融排出するための手法は、DTA 曲線で囲まれた全吸熱量と各温度毎の２
２段噴流床炉内ガス温度分布

生成ガス＋チャー

サイクロン

リダクタ部
軟化点温度以下領域

リダクタ

圧
力
容
器

コンバスタ〜リダクタ
軟化点〜流動点領域（遷移領域）

石炭
搬送ガス
チャー
石炭
搬送ガス

搬送ガス

コンバスタ部
流動点領域

空気

溶融スラグ

1300 1400 1500 1600 1700 1800 1900
ガス温度（Ｋ）

図3-5-754

二段噴流床ガス化炉内温度分布例

溶融開始温度

溶融終了温度

ガス供給装置

TJ炭例

全吸熱量

ガス制御
装置

粉体供給
装置

アルミナ管

示差熱曲線

−20

コンプレッサー電気ヒーター

−40

１段
N2

−60

示差熱曲線と溶融終了温度で囲まれる面積
を灰の融解に要した全吸熱量（融解潜熱）と見なし、各温度毎の吸熱面積比から灰中の液相割合を算出した

−80

−100
1000

1200

1400

２段

CO2

1600

1800

高速ビデオカメラ

2000

1000mmDTAヒートフロー、 μV

0

３段
灰付着プ
ローブ

温度、 K

冷却媒体

図3-5-8 示差熱分析結果の一例（TJ炭）
４段

φ50mm

吸熱量の比から算出するものである。灰溶融性の異なる
国内 TJ 炭、豪州 MC 炭など５炭種の灰中液相率を図 3-

排気5-9 に示す。図からも明らかなように、灰中液相率は、

図3-5-10

軟化点温度以下から急激に増加し、その傾きは炭種によ

灰付着特性基礎試験装置

り大きく異なる。
各温度条件において、テストプローブ上に付着した灰量
º

灰付着実験

と衝突した灰量の比から付着率を算出し、前述の DTA

灰付着実験には、図 3-5-10 に示す粉体落下型電気炉より算出した灰中液相率との相関をとると、図 3-5-11を使用した。電気炉内には直径 50mm、長さ１ m のアの通り、良い相関が見られ、炭種によらず灰中液相率かルミナ管を設置し、管上部より約 650mm の位置にガスら灰付着性を推定可能であることがわかったÐ。

化炉内壁を模擬した灰付着テストプローブを挿入した。
本プローブは、SUS304 製円盤の上に、２トン/日炉リ»

灰付着判別モデルの構築

ダクタ壁と同材質の耐火材を施工し、付着面を水平面に以上の結果から、灰中液相率を基に灰付着率を判断す

対して約 60 度の角度（安息角）に傾けて設置した。ア

る基本概念（図 3-5-12）をたて、灰付着判別モデルを

ルミナ管内のガス温度を 1473K 〜 1873K 間で 100K 毎に構築した。図 3-5-13 に石炭ガス化シミュレーションに変えて灰の溶融状態を変化させ、表面温度を 873K のテ組み込むための付着判別アルゴリズムを示す。コンバス

ストプローブに衝突させて４炭種の灰の付着性を調べた。

SB

ULTJ MCNL

灰中液相率φL；wt％

80

灰付着率 f dep；wt％

100

100
：IDT

TJ
MC
NL
SB
UL

60

40

60
40

灰付着率 f dep；27％
20
0

20

灰付着特性基礎試験装置にて計測した4炭種の灰付着率と灰中液相率との相関曲線

TJ
MC
NL
UL

80

灰中液相率φL＝55％
の場合
0

20

40

60

80

100

灰中液相率φL；wt％
0
1200

1300

1400

1500

1600
1700
温度T、K

1800

図3-5-9 灰中液相率の比較

1900

2000

灰付着率 f dep（wt ％）＝

付着灰量
衝突灰量

m dep（kg）
×100
m dep（kg）

図3-5-11 灰付着率と灰中液相率の関係

電中研レビュー No.44 ● 55

炉壁

炉内壁に衝突した灰粒子の内、付着する確率（灰付着率；wt％）
は最高履歴温度T Pmaxで液相となった粒子の割合（灰中液相率；wt％）と相関

粒子が壁面に衝突

付着粒子

YES

コンバスタ壁面？
反発粒子
溶融粒子
未溶融粒子

衝突

粒子最高履歴温度TPmax
粒子最高履歴温度TPmax
から付着率fdepを計算から付着率fを計算

最高履歴温度TPmaxの灰粒子群

図3-5-12

灰付着判別の基本概念

一様乱数Ｒを発生

タ壁は、溶融スラグが流下しているため、衝突した灰粒

YES

Ｒ≦f
dep
Ｒ≦f

子は 100%付着するとした。リダクタ壁に衝突した灰粒子は、最大履歴温度に基づく灰中液相率に従い灰付着率が算出され、付着判別される。

3-5-4

粒子追跡を継続

灰付着数値解析

図3-5-13

粒子付着

灰付着判別アルゴリズム

構築した灰付着判別アルゴリズムを当研究所で開発しの２条件について解析を行い、運転条件が灰付着性に及た石炭ガス化炉数値解析コード Áに導入し、２トン/日

ぼす影響をみた。解析結果として、図 3-5-14 に示すガ

ガス化炉を対象に TJ 炭の空気比λ =0.46 およびλ =0.59

ス化炉壁面への単位面積当たりの灰付着量分布が得られ

5.0
A

A

約2〜3mmの付着4.0

約10〜15mmの付着

A

リダクタ

B

B

3.0
約15mmの付着

約2〜3mmの付着

B
2.0
C
リダクタ
バーナー

角状スラグは形成されず

1.0

バーナーに角状スラグ形成

コンバスタ

C
コンバスタ
バーナー
0

1
2
3
4
壁面灰付着量、kg/m2・s（×E-03）_

図3-5-1456C

低空気比（λ=0.46）運転時

0

1
2
3
4
壁面灰付着量、kg/m2・s（×E-03）`

高空気比（λ=0.59）運転時

２トン/日炉を対象とした灰付着解析結果

た。低空気比のλ =0.46 の場合、リダクタバーナ近傍へ

以上、灰粒子の温度変化に伴う連続的溶融過程を灰中の著しい付着が見られ、バーナ先端部に角状スラグの形

液相率で与え、付着/反発挙動をこれに基づき確率的に成が予想される。リダクタ壁面においても、コンバスタ

判別するモデルを構築した。

を通過した灰およびリダクタ供給石炭中の灰による付着

当所で開発した石炭ガス化炉数値解析コードに同モデが予測される。一方、高空気比のλ =0.59 の場合、0.46

ルを導入することにより、従来困難であった灰軟化点以の場合に比べて、リダクタバーナ周辺への灰付着量は少下の領域を含むリダクタ部全域にわたり、灰付着性予測

なく、むしろリダクタ中段から上段にかけての付着量がが可能となった。

増加した。これらの解析結果は、ガス化試験後の炉内灰付着結果とも定性的に良く一致したÑ。

今後はさらに、付着物の焼結性を考慮し、伝熱阻害量を評価するモデルの構築へと研究をすすめる予定である。

３−６

今後の展開

引き続き、実証機および商用機の設計・運転支援に向IGCC の導入が予想されるため、これらの炉内現象を解けた評価技術の確立をめざす。数値解析手法の改良・高析・評価できるツールの開発が望まれる。当研究所では、度化はもとより、基礎実験装置や高度分析装置を駆使しこれらの課題に対応するため、様々なガス化炉方式に対てガス化反応特性の解明、灰・チャーの生成・付着機構

応でき、炉内現象の詳細な計測・解析が可能な２〜３トの解明・モデル化を進め、シミュレーションへの導入を

ン/日規模の研究用ガス化炉を平成 15 年度に設置し、実図る。また、炭種の違いによりガス化炉特性が大きく変験技術と計算技術を両輪とするガス化炉の総合評価技術

化することが予想されるため、炭種評価技術の開発に重を開発することによって、IGCC の実用化・本格導入に点を置き検討を進める予定である。

貢献していきたいと考えている。

IGCC さらに、将来は海外機を含む様々な方式の電中研レビュー No.44 ● 57

コラム：石炭中の鉱物について
石炭中の鉱物は、¸石炭の元となった植物に含晶構造が変化し別の鉱物となる。図１¹）や高温でまれていた無機成分、¹混入した周囲の土砂、ºの平衡組成（図２¹）について研究が進められてい

石炭層形成後に、地下水や熱水から移行した無機る。

成分から形成されているが、量的には、¹が大部分を占めるといわれており¸、石炭中に含まれる

参考文献

主 な 鉱 物 と し て は 、 石 英 （ S i O 2） や 粘 土 鉱 物¸

二宮、「火炉内の石炭灰の溶融・付着現象の解（Kaolinite（Al2Si2O72H2O）、Illite （KAl2（Si3Al）

明 」、 日 本 エ ネ ル ギ ー 学 会 誌 、 V o l . 7 7 、 N o . 3 、O10（OH）2）など）、Calcite（CaCO3）や Dolomite（1998）、pp. 177-186、（Ca,Mg（CO3）2）などの炭酸塩鉱物、Pyrite（FeS2）¹

Gpysum（（CaSO42H2O）などがあげられる¹。

Godon Couch、「Understanding slagging andfouling during pf combustion」、IEA Coal

高温の炉内で生成する石炭灰の溶融特性（付着research、（1994）

特性）などを精度よく予測するためには、化合物
組成ではなく、鉱物組成に基づいて検討することが重要と考えられており、石炭中に含まれる鉱物

SiO2
1698℃

が高温にさらされたときに発生する形態変化（結〜1590℃
Cristoballite

Kaolinite

Metakaolinite

IlliteMullite+Amorphous Quartz

Lime

CalcitePytiteHematite/Magnetite

Gysum

Anhydrite

0

1698℃

Amorphous alumioslicates

500

1083℃
1178℃
Favalite58

Mullite

1205℃
1380℃

1177℃

1000
Temperature（℃）

2000〜1840℃
1148℃

Wustite

1329℃
FeO

図１

Tridymite¹

石炭含有鉱物の高温での反応

Hercynite

Corundum〜1750℃
Al2O3

図２

FeO-Al2O3-SiO2系の高温平衡組成