JP3033466B2 - 高炉の操業方法 - Google Patents
高炉の操業方法Info
- Publication number
- JP3033466B2 JP3033466B2 JP7106498A JP10649895A JP3033466B2 JP 3033466 B2 JP3033466 B2 JP 3033466B2 JP 7106498 A JP7106498 A JP 7106498A JP 10649895 A JP10649895 A JP 10649895A JP 3033466 B2 JP3033466 B2 JP 3033466B2
- Authority
- JP
- Japan
- Prior art keywords
- furnace
- distribution
- model
- state
- blast furnace
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
Landscapes
- Manufacture Of Iron (AREA)
Description
を含む3次元方向の非定常状態をモデル化した3次元非
定常モデルと炉頂からの装入物の分布モデルとを用いて
非定常状態における高炉の操業をシミュレートし、シミ
ュレート結果が実炉における計測値と一致するようにモ
デルの一部のプロセス定数を修正した後、修正したモデ
ルを用いて、溶銑温度、溶銑中Si濃度、炉内の圧力分
布、燃料比の時間的変化等、操業結果の良否に影響を及
ぼす因子の非定常状態が目標の操業結果が得られる非定
常状態となるような送風条件及び装入物分布を推定して
実炉の非定常状態における操業条件とする高炉の操業方
法に関する。
存し、互いに反応、相変化、及び熱交換を行いながら移
動する複雑な反応容器である。このような高炉の炉内現
象を、基礎実験、解体調査等の知見を基に、移動速度論
的に整理し、組み合わせることにより、高炉を近似的に
数学モデル化することが可能である。数学モデルとして
は1次元定常モデル(特公昭61−6122号、特公昭
61−60124号、特公昭63−24044号、特公
平1−12805号公報)、2次元定常モデル(特公平
6−2886号公報)等が開発されてきた。これらのモ
デルでは、高炉の3次元的な挙動を高さ方向のみ取り扱
い(1次元モデル)、又は円周方向の分布を均一とし、
さらに時間的にも状態の変化がないとした定常状態のみ
の解析に限定し(2次元モデル)、モデルを簡略化する
ことで炉内状態値を算出する計算量の減少を図って計算
機の応答性を上げている。
下部の円周方向に設けられている複数の羽口より熱風を
送風するため、炉壁,充填層等の固体温度は円周方向に
分布し、高温域における反応挙動(SiOガスのような
気体の固体からの蒸発反応等)に差を生じさせると考え
られるが、上述のような1次元モデルでは高さ方向のみ
の分布を考慮し、また2次元モデルでは周方向の分布を
均一とし、定常状態のみの解析を行って計算を簡略化し
ているので、炉内状態の解析精度が低い。従って、この
ようなモデルによる炉内状態の解析結果を用いた操業シ
ミュレーションで高炉の操業条件を決定した場合、高炉
のスタートアップ時等のように、炉内状態の時間経過に
ともなう変化が大きい非定常状態における高炉の的確な
操業条件を推定できないという問題がある。
になされたものであって、高炉の高さ方向、半径方向の
みならず、周方向を含む3次元方向の非定常状態を記述
する3次元非定常モデルを用いて高炉の操業をシミュレ
ートし、実炉での計測値と一致するように3次元非定常
モデルのプロセス定数及び装入物分布モデルによって算
出される炉頂での鉄鉱石及びコークスのそれぞれの堆積
層厚さ等の分布状態を修正し、修正後のモデルにより高
炉の操業を高精度にシミュレートすることにより、実炉
の操業条件が的確に決定できる高炉の操業方法の提供を
目的とする。
は、その壁面の周方向に複数個設けられている羽口から
それぞれ炉内に吹き込まれた熱風が羽口から炉芯にかけ
て堆積しているコークスを通過して炉頂に上昇し、炉頂
から鉄鉱石とともに装入されているコークスをガス化し
て鉄鉱石を還元するとともに、還元された鉄鉱石を溶融
し、溶銑が炉芯のコークス層を通って滴下し、炉床に湯
溜りを形成すべくなしてある高炉の、時間経過に伴って
変化する非定常状態を気相・液相・固相それぞれの物質
収支、運動量収支、及びエネルギー収支に関連付けて表
した状態モデルを用いて推測して非定常状態における実
炉の操業条件を決定する高炉の操業方法において、炉頂
における鉄鉱石及びコークスの分布状態を表した装入物
分布モデルにより鉄鉱石及びコークスの炉頂における粒
子径及び堆積層厚さの分布状態を算出し、炉内の高さ方
向、径方向、及び周方向における炉内の非定常状態を気
相・液相・固相それぞれの物質収支、運動量収支、及び
エネルギー収支に関する微分方程式群により表した3次
元状態モデルに、時間経過に伴ってその値が変化し得る
還元反応速度定数、炉心の空隙率、固気間熱伝達係数及
び気液接触面積を含むプロセス定数を与え、前記装入物
分布モデルにより算出された前記分布状態を炉頂におけ
る境界条件として、また、羽口への送風条件を羽口にお
ける境界条件として与えて高炉の操業をシミュレート
し、銑鉄生産速度、炉頂ガス組成分布、炉頂温度分布、
炉内圧力分布、溶銑温度及び、溶銑中Si濃度等を含
む、高炉の操業結果の良否に影響を及ぼす要因の将来の
時間変化状態を推定し、該要因の時間変化状態を実炉に
て測定し、推定値と実測値とが一致するように、3次元
状態モデルの前記プロセス定数と、前記装入物分布モデ
ルにより算出された前記分布状態とを修正した後、前記
要因の時間変化状態が、目的の操業結果が得られる時間
変化状態となるように、修正後の3次元状態モデルによ
る高炉の操業シミュレーションを基にして装入物分布と
羽口への送風条件とを推定し、推定した装入物分布及び
羽口への送風条件を実炉の操業条件とすることを特徴と
する。
方法という)の手順の概略を示すフローチャート、図2
は本発明方法を実施する高炉の模式図である。炉内の気
相・液相・固相の状態を記述する物質収支・運動量収支
・エネルギー収支の数式からなり、高炉の周方向を含む
3次元方向の非定常状態を推定する3次元非定常モデル
(後に詳述する)に、鉄鉱石の還元反応速度定数・炉芯
の空隙率・固気間熱伝達係数・気液接触面積等、時間経
過に伴ってその値が変化し得るプロセス定数と、鉱石/
コークスの落下高さ,落下速度,鉱石/コークスの装入
具である大ベルの角度,装入物と大ベル面との摩擦係数
等を導入した所定の計算式からなる装入物分布モデルに
より算出された、炉頂における鉄鉱石・コークスのそれ
ぞれの粒子径、堆積層厚さ等の分布状態と、解析対象の
実炉における羽口からの送風条件とを与えて解析対象の
高炉の操業をシミュレートし、高炉の操業結果の良否に
影響を及ぼす要因である、銑鉄生産速度・炉頂のガス組
成分布・炉頂温度分布・炉内の圧力分布・溶銑温度・溶
銑中Si濃度を推定する(S1〜S3)。一方、これら
の銑鉄生産速度・炉頂のガス組成分布・炉頂温度分布・
炉内の圧力分布・溶銑温度・溶銑中Si濃度を、ガス分
析計、垂直ゾンデ、圧力計、水平ゾンデ、温度計等の計
測端により実炉において計測する(S4)。
とともに変化する流体の流れ状況、固体の接触状況等に
よりその値が変化し得るプロセス定数、例えば、鉄鉱石
の還元反応速度定数、炉芯の空隙率、固気間熱伝達係
数、気液接触面積等と、装入物分布モデルにより算出し
た堆積層厚さとを、前述の両モデルによる推定値と実炉
での実測値とが一致するように修正する(S5)。その
後、修正した3次元非定常モデルを用いて、操業結果に
影響を及ぼす要因である溶銑温度、溶銑中Si濃度、炉
内の圧力分布、コークスの燃料比の時間的変化等の非定
常状態が目標の操業結果が得られる非定常状態となる装
入物分布及び羽口への送風条件を推測し(S6)、これ
らの条件を実炉の操業条件とする(S7)。
モデルの構成概要を示す模式図である。羽口から炉芯に
かけて逆V字状に堆積しているコークスのうち、羽口に
近いレースウェイのコークスが羽口からの送風によって
燃焼・運動し、羽口から送風された熱風とコークスとの
反応熱によりガス温度が上昇する。ガスはコークス層の
表面層を通って炉頂に上昇し、炉頂から交互に装入され
ている鉱石及びコークスに付着している水分を蒸発させ
て鉱石・コークスを乾燥、除々に加熱し、さらにコーク
スをH2 OとCO2 とによりガス化し、ガス化により発
生するH2 とCOとによって鉱石中の鉄と結合している
酸素の大半が除かれる(間接還元)。間接還元された鉱
石は熱と荷重により半溶融状になり、隣り合った粒子同
志が溶着し、逆V字状に堆積しているコークス層の表面
に融着帯が形成される。一旦溶着した鉱石はより高温の
環境下でさらに還元され(直接還元)、溶銑は逆V字状
に堆積しているコークスの表面層(滴下帯)を浸炭反応
を起こしながら、炉芯コークスを通過して炉底に滴下
し、湯溜りを形成する。
デルを記述する支配方程式である気液固3相それぞれの
物質・運動量・エネルギーの収支式について具体的に説
明する。なお、式中の添え字、gはガスを、sは固体
を、またlは液体を表す。 (a) 気相の物質収支 気相の成分kの質量収支は式(1) で表される。
気相に関する連続の式として式(2)を得る。
の圧力損失をモデル化したERGUN 式に従うと仮定する
と、気相における運動量収支は式(4) で表される。
及びコークス粒子内の成分kの物質収支は式(5) で表さ
れる。但し、nは焼結鉱及びコークスを表す。
得る。
固相に関する連続の式として式(7)を得る。
相のエネルギー収支は式(8) で表される。
ルを用いる。Kinematic モデルの支配方程式は、重力方
向をy方向、流速ベクトルの成分を(Us ,V s ,
Ws )とした場合、式(9) で表される。
である式(7) に代入して求める。 (g) 粒子体積に関する保存式 ガス化や、溶銑への浸炭反応等の固体粒子の消滅反応を
界面反応として記述すれば、粒子1個の体積をVpnとし
た場合、粒子体積変化はラグランシュ微分を用いて式(1
0)のように表される。
と式(12)が得られる。
支は式(13)で表される。なお、液相は分散相であるか
ら、熱拡散は考慮しないものとする。
支として滴下速度を与える。
る3次元非定常モデルを解析して炉内の3次元非定常状
態を推定する手順を示すフローチャートである。モデル
は完全自立型であり、装入物の炉頂での条件、羽口への
送風条件、及び炉壁での熱的条件を与えると、実炉運転
と基本的に同じ高炉の操業をシミュレートし、銑鉄生産
速度、及び3次元方向での炉頂のガス組成分布、炉頂温
度分布、炉内の圧力分布、溶銑温度、溶銑中Si濃度等
を算出する。
件と初期条件との下に連成し、時間進展することによ
り、各格子点の3次元方向の非定常状態が算出される。
このとき、格子系にはスタガード格子を採用し、任意形
状領域の処理は、境界適合格子による差分法を用いた。
解析アルゴリズムは流体解析で一般的な流速と圧力との
同時緩和法であるSOLA法を改良して用いた。図4に示し
た順に、ガス流れモデル、固体流れモデル、ガス相温度
Tg ・固相温度Ts ・液相温度Tl の連立解、ガス相の
質量分率xgk、液相の質量分率xlk、固相の質量分率x
sk、及び粒子径を求める。以上を格子点毎に繰り返し、
3次元方向の炉内状態が算出される。
の周方向の各断面におけるガス流れと温度分布との状態
を示す図であり、図5(b) は羽口近傍の液相及び固相に
おけるガス流れと温度分布との解析結果を2次元解析と
比較して示した図である。図から明らかなように、3次
元解析では羽口近傍の固相における周方向の温度分布が
解析されている。従って、高温域における反応挙動(S
iOガスのような気体の固体からの蒸発反応)の解析精
度が高いと考えられる。
定常モデルにより高炉のスタートアップ時における液相
及びガス相での非定常な流れ及び温度を解析した結果を
示す図である。このように、急激な変化を伴う非定常な
炉内挙動についても安定に計算することができ、実炉で
の計測値(温度の推移、炉頂ガス組成等)との一致も良
いことが分かった。なお、167 時間以上では炉内の状態
変数はほとんど変化しなくなるので、定常解が得られ
る。
さ方向、半径方向のみならず、周方向を含む3次元方向
の非定常状態を記述する3次元非定常モデルを用いて高
炉の操業をシミュレートし、実炉での計測値と一致する
ように3次元非定常モデルのプロセス定数及び装入物分
布モデルによって算出される炉頂での鉄鉱石及びコーク
スのそれぞれの堆積層厚さ等の分布状態を修正し、修正
後のモデルにより高炉の操業を高精度にシミュレートす
るので、実炉の操業条件が的確に決定できるという優れ
た効果を奏する。
る。
概要を示す模式図である。
手順を示すフローチャートである。
較を示す図である。
高炉のスタートアップ時の計算例(気体及び液体の流れ
と温度)を示す図である。
Claims (1)
- 【請求項1】 その壁面の周方向に複数個設けられてい
る羽口からそれぞれ炉内に吹き込まれた熱風が羽口から
炉芯にかけて堆積しているコークスを通過して炉頂に上
昇し、炉頂から鉄鉱石とともに装入されているコークス
をガス化して鉄鉱石を還元するとともに、還元された鉄
鉱石を溶融し、溶銑が炉芯のコークス層を通って滴下
し、炉床に湯溜りを形成すべくなしてある高炉の、 時間経過に伴って変化する非定常状態を気相・液相・固
相それぞれの物質収支、運動量収支、及びエネルギー収
支に関連付けて表した状態モデルを用いて推測して非定
常状態における実炉の操業条件を決定する高炉の操業方
法において、 炉頂における鉄鉱石及びコークスの分布状態を表した装
入物分布モデルにより鉄鉱石及びコークスの炉頂におけ
る粒子径及び堆積層厚さの分布状態を算出し、 炉内の高さ方向、径方向、及び周方向における炉内の非
定常状態を気相・液相・固相それぞれの物質収支、運動
量収支、及びエネルギー収支に関する微分方程式群によ
り表した3次元状態モデルに、 時間経過に伴ってその値が変化し得る還元反応速度定
数、炉心の空隙率、固気間熱伝達係数及び気液接触面積
を含むプロセス定数を与え、前記装入物分布モデルによ
り算出された前記分布状態を炉頂における境界条件とし
て、また、羽口への送風条件を羽口における境界条件と
して与えて高炉の操業をシミュレートし、銑鉄生産速
度、炉頂ガス組成分布、炉頂温度分布、炉内圧力分布、
溶銑温度及び、溶銑中Si濃度等を含む、高炉の操業結
果の良否に影響を及ぼす要因の将来の時間変化状態を推
定し、 該要因の時間変化状態を実炉にて測定し、 推定値と実測値とが一致するように、3次元状態モデル
の前記プロセス定数と、前記装入物分布モデルにより算
出された前記分布状態とを修正した後、 前記要因の時間変化状態が、目的の操業結果が得られる
時間変化状態となるように、修正後の3次元状態モデル
による高炉の操業シミュレーションを基にして装入物分
布と羽口への送風条件とを推定し、 推定した装入物分布及び羽口への送風条件を実炉の操業
条件とすることを特徴とする高炉の操業方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP7106498A JP3033466B2 (ja) | 1995-04-28 | 1995-04-28 | 高炉の操業方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP7106498A JP3033466B2 (ja) | 1995-04-28 | 1995-04-28 | 高炉の操業方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH08295910A JPH08295910A (ja) | 1996-11-12 |
JP3033466B2 true JP3033466B2 (ja) | 2000-04-17 |
Family
ID=14435107
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP7106498A Expired - Lifetime JP3033466B2 (ja) | 1995-04-28 | 1995-04-28 | 高炉の操業方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3033466B2 (ja) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AT407993B (de) | 1999-03-03 | 2001-07-25 | Voest Alpine Ind Anlagen | Verfahren zur optimierung von auslegung und betrieb eines reduktionsverfahrens |
JP4968388B2 (ja) * | 2010-02-05 | 2012-07-04 | Jfeスチール株式会社 | 流体系の温度推定方法、流体系における物質成分の濃度および温度推定方法、流体系の温度分布推定方法、流体系の温度モニタリング方法、溶融金属設備の溶融金属温度制御方法、溶融金属設備の濃度および温度推定方法、溶融亜鉛めっき鋼板ならびに流体系の温度推定装置 |
WO2011122610A1 (ja) | 2010-03-30 | 2011-10-06 | 新日鉄ソリューションズ株式会社 | 情報処理装置、情報処理方法及びプログラム |
JP4984000B1 (ja) * | 2010-09-30 | 2012-07-25 | Jfeスチール株式会社 | 流体系の温度推定方法、流体系の温度分布推定方法、流体系の温度分布モニタリング方法、温度推定装置、溶融亜鉛めっきポット内の溶融亜鉛温度制御方法、溶融亜鉛めっき鋼板、およびタンディッシュ内の溶鋼温度制御方法 |
JP5583614B2 (ja) * | 2011-02-23 | 2014-09-03 | 新日鐵住金株式会社 | 高炉操業状態制御システム、操業制御装置、高炉操業状態制御方法およびプログラム |
JP5583615B2 (ja) * | 2011-02-23 | 2014-09-03 | 新日鐵住金株式会社 | 高炉操業特徴量演算装置、高炉操業特徴量演算方法及びプログラム |
JP5862470B2 (ja) * | 2012-06-13 | 2016-02-16 | 新日鐵住金株式会社 | 高炉休風方法 |
CN102703627B (zh) * | 2012-06-19 | 2013-10-30 | 首钢总公司 | 一种预测高炉煤气利用率的方法 |
JP6248550B2 (ja) * | 2013-11-01 | 2017-12-20 | 新日鐵住金株式会社 | 高炉操業条件の決定方法 |
KR101482469B1 (ko) * | 2013-12-24 | 2015-01-13 | 주식회사 포스코 | 고로 공정 가상 현실 교육 시스템 |
JP6311659B2 (ja) * | 2015-06-22 | 2018-04-18 | Jfeスチール株式会社 | 高炉内の層厚分布の推定方法、高炉の操業方法、および高炉内の層厚分布の推定装置 |
JP6531782B2 (ja) * | 2016-08-02 | 2019-06-19 | Jfeスチール株式会社 | 溶銑温度予測方法、溶銑温度予測装置、高炉の操業方法、操業ガイダンス装置、溶銑温度制御方法、及び溶銑温度制御装置 |
JP6493447B2 (ja) * | 2016-08-02 | 2019-04-03 | Jfeスチール株式会社 | 溶銑温度予測方法、溶銑温度予測装置、高炉の操業方法、操業ガイダンス装置、溶銑温度制御方法、及び溶銑温度制御装置 |
JP7207118B2 (ja) * | 2019-04-11 | 2023-01-18 | 日本製鉄株式会社 | 装入方法決定方法、装入方法決定装置及び装入方法決定プログラム |
JP7335500B2 (ja) * | 2019-09-17 | 2023-08-30 | 日本製鉄株式会社 | パラメータ決定装置、パラメータ決定プログラム及びパラメータ決定方法 |
JP7464033B2 (ja) * | 2021-01-13 | 2024-04-09 | Jfeスチール株式会社 | 高炉融着層の通気性推定方法、高炉融着層の通気性推定装置、及び高炉の操業方法 |
CN114058749B (zh) * | 2021-11-29 | 2023-02-10 | 山西晋煤集团技术研究院有限责任公司 | 一种无烟块煤在高炉中的劣化研究方法 |
-
1995
- 1995-04-28 JP JP7106498A patent/JP3033466B2/ja not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
JPH08295910A (ja) | 1996-11-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3033466B2 (ja) | 高炉の操業方法 | |
Yang et al. | Numerical modeling of reaction and flow characteristics in a blast furnace with consideration of layered burden | |
Kitaev et al. | Heat exchange in shaft furnaces | |
Fu et al. | CFD modeling of multiphase reacting flow in blast furnace shaft with layered burden | |
Panjkovic et al. | Numerical modelling of iron flow and heat transfer in blast furnace hearth | |
Austin et al. | Analysis of actual blast furnace operations and evaluation of static liquid holdup effects by the four fluid model | |
Wen et al. | Fluidized‐bed heat transfer: A generalized dense‐phase correlation | |
Majumder et al. | Virtual indurator: A tool for simulation of induration of wet iron ore pellets on a moving grate | |
JP5583614B2 (ja) | 高炉操業状態制御システム、操業制御装置、高炉操業状態制御方法およびプログラム | |
Pintowantoro et al. | Numerical analysis of static holdup of fine particles in blast furnace | |
CN116306375A (zh) | 高炉生产运行状态的数字孪生方法 | |
Jindal et al. | A reduced order thermo-chemical model for blast furnace for real time simulation | |
JP2727563B2 (ja) | 高炉操業方法 | |
CN110427693B (zh) | 带有喷枪的工业反应器内流体仿真计算方法 | |
JP2012172222A (ja) | 高炉操業特徴量演算装置、高炉操業特徴量演算方法及びプログラム | |
Panjkovic et al. | Computational Fluid Dynamics Modelling of Iron Flow and Heat Transfer in the Iron Blast Furnace Hearth | |
Yang et al. | CFD prediction for the performance of a heat treatment furnace | |
Pomerleau et al. | A first principle simulator of an iron oxide pellet induration furnace–an application to optimal tuning | |
Nisioka et al. | Development of Mathematical Models for Blast Furnaces | |
JPS6055561B2 (ja) | 高炉の操業方法 | |
WO1997012064A9 (en) | A method for the determination of the gas flux distribution in a blast furnace | |
Dash et al. | Sinter Strand Control; A Simulation Investment | |
Yang et al. | Numerical modelling of temperature distribution and pressure drop through the layered burden loading in a blast furnace | |
Hekkala et al. | Mathematical model of heat and mass transfer in the steel belt sintering process | |
Makkiabadi et al. | Comparative Study on Induration of Wet Iron Ore Pellets in Different Straight Grate Furnaces in Iran Using Mathematical Models |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080218 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090218 Year of fee payment: 9 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100218 Year of fee payment: 10 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100218 Year of fee payment: 10 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110218 Year of fee payment: 11 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120218 Year of fee payment: 12 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120218 Year of fee payment: 12 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130218 Year of fee payment: 13 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130218 Year of fee payment: 13 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140218 Year of fee payment: 14 |
|
EXPY | Cancellation because of completion of term |