CN1963453A - 研究煤热解过程中半焦收缩特性的方法 - Google Patents

研究煤热解过程中半焦收缩特性的方法 Download PDF

Info

Publication number
CN1963453A
CN1963453A CNA2006101022094A CN200610102209A CN1963453A CN 1963453 A CN1963453 A CN 1963453A CN A2006101022094 A CNA2006101022094 A CN A2006101022094A CN 200610102209 A CN200610102209 A CN 200610102209A CN 1963453 A CN1963453 A CN 1963453A
Authority
CN
China
Prior art keywords
semicoke
coal
sample
shrinkage
temperature
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.)
Granted
Application number
CNA2006101022094A
Other languages
English (en)
Other versions
CN100395536C (zh
Inventor
谢克昌
常海洲
曾凡桂
李文英
王传格
李军
张代林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Taiyuan University of Technology
Original Assignee
Taiyuan University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Taiyuan University of Technology filed Critical Taiyuan University of Technology
Priority to CNB2006101022094A priority Critical patent/CN100395536C/zh
Publication of CN1963453A publication Critical patent/CN1963453A/zh
Application granted granted Critical
Publication of CN100395536C publication Critical patent/CN100395536C/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

一种研究煤热解过程中半焦收缩特性的方法，属于煤科学和煤化工领域，具体来讲涉及一种用图象分析的方法准确地得到在煤的热解过程中半焦收缩特性的方法的技术方案，其特征在于是一种以显微镜和热台为主要实验手段，以煤及其显微组分颗粒为研究对象，通过显微镜下拍摄到的显微图片获得样品在热解过程中的动态变化特点，用图像分析方法准确得到煤热解过程中半焦收缩特性的方法。本发明克服了用热重分析研究煤热分解的局限性，采用图像分析的方法，既直观又定量地获得了半焦收缩的分段特点、影响结束温度的因素及其它动力学参数。避免了延时排焦，有利于节能并提高生产效率和产品质量。

Description

研究煤热解过程中半焦收缩特性的方法
技术领域
本发明研究煤热解过程中半焦收缩特性的方法，属于煤科学和煤化工领域，具体来讲涉及一种用图象分析的方法准确地得到在煤的热解过程中半焦收缩特性的方法的技术方案。
背景技术
因为煤的热解是煤的各种热转化过程中重要和基本的过程，所以使其成为煤化工研究所关注的一个重要问题。煤在热解中表现出的半焦收缩特性是炼焦技术中需要掌握的重要技术特征。热重分析是研究煤热分解的一种常用方法，以此能够获得煤热分解中脱挥发分的动力学特征，但却难以准确地得到半焦收缩过程的动力学特征，原因是半焦收缩过程虽然伴随着一定程度的热失重变化，但其更为重要的特点是固相体积收缩，而热重分析不能反映体积的变化。该发明所使用的显微热台装置能够观察到体积的变化。显微观察也是研究煤热转化过程的一种有效手段。Ma等人使用CCD数字摄像机拍摄到煤热解初始阶段挥发分析出的显微图像，观察了挥发分释放的微观现象，Alonso等人通过煤显微组分的光学显微图片研究了煤显微组分在热处理中的物理结构和反射率的变化，他们还利用光学显微图片阐明了半焦在不同温度下具有不同的形态形貌及其横剖面不同的光学特性；D.Diniyati等人利用显微图片考察了单一颗粒煤在燃烧中是如何形成煤灰的。我国学者康西栋等和叶道敏等利用显微热台将煤样制成块煤光片和煤砖光片考察了显微组分的粘结性、结焦性和软化流动温度等。图像分析不仅在煤岩组成的研究中是一种很有效的手段，而且也应用于其他方面，如煤的浮选、煤焦油中间相的形成、煤对CO2的吸附等。但是在显微图片基础上运用图像分析的方法准确得到煤热解过程中半焦收缩特性的问题一直尚未解决，其有关的技术方案还未见报道和公开。所以该发明研究半焦收缩特性所得到的结果将对炼焦中涉及到的排焦工艺提供依据和指导。
发明内容
本发明研究煤热解过程中半焦收缩特性的方法，目的在于解决上述现有技术中一直尚未解决的问题，从而公开一种以显微镜和热台为主要实验手段，以煤及其显微组分颗粒为研究对象，在显微图片基础上应用图像分析的方法准确得到煤热解过程中半焦收缩特性的技术方案。
本发明研究煤热解过程中半焦收缩特性的方法，其特征在于是一种以显微镜和热台为主要实验手段，以煤及其显微组分颗粒为研究对象，通过显微镜下拍摄到的显微图片获得样品在热解过程中的动态变化特点，用图像分析方法准确得到煤热解过程中半焦收缩特性的方法，其具体的方法步骤为：
I、取样及装置的调试
用毛细管取样，以保证对各样品的研究具有可比性将煤样横放在样品池中，各样品具有相同的厚度，将装有小于0.1mg样品颗粒的铂坩埚置于显微热台中，调节显微镜光强度为10，将热台通入流量为120mL/min的N2，接通循环冷却水，开启升温控温程序，在升温速率为5～40K/min时自室温加热至终温；
II、使用图像分析软件绘制半焦收缩曲线
利用徕卡显微镜自身配置的Qwin图像处理软件对热分解过程不同温度时采集的显微图片进行分析；由该软件控制的数字摄像机的拍摄条件为：曝光时间51.5ms，增益2.1，白平衡49、2.90、16；数码摄像机采集不同温度下样品颗粒的显微图片，用Qwin图像处理软件将起始温度时显微图片上煤样的灰度作为基准，灰度域值为1～10，求得热解过程中某一温度时样品颗粒具有此灰度域值的面积数值，从而算出该温度下的样品颗粒面积与起始时煤粒面积的比值，称为“剩余面积分数”，以剩余面积分数对温度做图，所得曲线称为半焦收缩曲线(见图1)；
III、由半焦收缩曲线考察半焦收缩特性
通过观察半焦收缩曲线，发现该曲线能够反映热解过程中半焦面积随温度的变化情况；在升温过程中半焦收缩的程度随温度而变化，半焦收缩过程分为三个阶段，半焦开始收缩时收缩程度较小，称为缓慢收缩阶段；随温度升高收缩程度变大，称为过渡阶段；最后收缩程度急剧增大，称为快速收缩阶段；根据半焦收缩曲线确定在显微镜下观察到的半焦收缩程度最大时的温度即结束温度，结束温度越高，表明煤样半焦收缩的性能越弱，半焦收缩的程度越小，反之，煤样的半焦收缩的程度越大；考察升温速率分别为5，10，20和40K/min的影响，发现结束温度皆随升温速率的增大而变高；变质程度较高的煤及其显微组分的结束温度高于变质程度较低的煤及对应显微组分；对于原煤和显微组分的结束温度高低顺序为：惰质组＞原煤＞镜质组；
IV、依据半焦收缩曲线求得半焦收缩的动力学参数
依据半焦收缩曲线，利用一级反应动力学方程：
1 n ( ( - 1 n ( 1 - x ) ) / T 2 ) = 1 n [ AR φE ( 1 - 2 RT E ) ] - E RT
式中，1-x：样品的剩余面积分数；φ：升温速率，K·s-1；E：活化能，J·mol-1；A：频率因子，s-1，用上式左边对 作图，得一直线，斜率为
Figure A20061010220900073
求得活化能E；截距为 据此求得频率因子A，再依据k＝Aexp(-E/RT)求得反应速率常数k，s-1，各样品的活化能、频率因子和反应速率常数皆随缓慢收缩阶段、过渡阶段和快速收缩阶段依次增大。
本发明研究煤热解过程中半焦收缩特性的方法的优点在于：克服了用热重分析研究煤热分解的局限性，采用图像分析的方法，通过半焦收缩曲线既直观又定量地获得了半焦收缩的分段特点、影响结束温度的因素(升温速率，煤种和显微组分)及动力学参数等。煤的收缩度是炼焦中重要指标之一，半焦收缩结束温度的高低反映了煤的收缩度的大小，同时对准确把握各种配煤的排焦时间提供参考，避免延时排焦，有利于节能并提高生产效率和产品质量，所以可对整个排焦工艺提供依据和指导；
附图说明
图1是升温速率为5K/min时神东惰质组的剩余面积分数随温度的变化曲线。FA代表剩余面积分数，T代表温度，S代表半焦收缩曲线，a、b和c分别代表缓慢收缩阶段、过渡阶段和快速收缩阶段，Te代表半焦收缩的结束温度。
具体实施方式
具体实施方式I
用毛细管取样，以保证对各样品的研究具有可比性将煤样横放在样品池中，各样品具有相同的厚度，将装有神东惰质组颗粒的铂坩埚置于显微热台中，调节显微镜光强度为10，将热台通入N2(流量为120mL/min)，接通循环冷却水，开启升温控温程序在升温速率为5K/min时自室温加热至终温；利用显微镜系统配置的图像分析软件，选择曝光时间51.5ms，增益2.1，白平衡49、2.90、16使图像质量最佳，每隔一定温度由数码摄像机同步采集显微图片并保存以用于图像分析。煤粒的热分解在横向和纵向上都有所收缩，表现在显微图片上便是煤粒面积的缩小。利用图像处理软件对热分解过程不同温度时采集的显微图片进行分析。以起始温度时显微图片上煤样的灰度作为基准，域值为1～10，求得该灰度域值时煤样的面积。表1是升温速率为5K/min时神东惰质组颗粒在不同温度时的面积及剩余面积分数：
表1：神东惰质组的面积及剩余面积分数随温度的变化
温度/℃ 面积/像素   剩余面积分数
20 111406   1
30 111642   1.002178649
40 111655   1.002178649
    50     111712     1.002904866
    60     111517     1.000726216
    70     111564     1.001452433
    80     111491     1.000726216
    90     111537     1.001452433
    100     111359     0.999273784
    110     111265     0.998547567
    120     111140     0.997821351
    130     110792     0.994190269
    140     110501     0.992011619
    150     109881     0.986201888
    160     109245     0.980392157
    170     108478     0.973856209
    180     106607     0.957153232
    200     105471     0.946986202
    210     104210     0.935366739
    220     102853     0.92302106
    230     101176     0.907770516
    240     98663     0.885257807
    250     95379     0.85620915
    260     91469     0.821 350763
    270     87275     0.783587509
    280     84584     0.758896151
    290     83207     0.746550472
    300     82799     0.74291939
    310     83065     0.745824256
    320     84159     0.755265069
    330     85386     0.766158315
    340     86686     0.777777778
    350     88484     0.794480755
    360     89907     0.806826434
    370     91392     0.820624546
    380     92586     0.830791576
    390     93590     0.840232389
    400     93986     0.843863471
    410     94173     0.845315904
    420     94321     0.846768337
    430     94140     0.845315904
    440     93754     0.841684822
    450     93338     0.83805374
    460     92896     0.833696442
    470     92268     0.82788671
    480     91567     0.822076979
    490     90702     0.814088598
    500     89665     0.804647785
    510     88391     0.793028322
    520     86779     0.779230211
    530     84174     0.755265069
    532     83949     0.753812636
    534     83452     0.748729121
    536     83084     0.745824256
    538     82437     0.740014524
    540     82011     0.736383442
    542     81259     0.729121278
    544     80426     0.721859114
    546     79741     0.716049383
    548     78800     0.707334786
    550     78075     0.700798838
    552     76359     0.685548293
    554     75114     0.673928831
    556     71967     0.645606391
    558     70401     0.631808279
    560     69170     0.620915033
    562     68007     0.610748003
    564     65843     0.59114016
    566     63859     0.572984749
    568     61660     0.553376906
    570     59691     0.535947712
    572     57019     0.511982571
    574     54251     0.487291213
    576     51269     0.460421206
    578     48067     0.431372549
    580     45392     0.407407407
    582     42098     0.377632534
    584     39508     0.354393609
    586     36663     0.328976035
    588     34504     0.309368192
    590     30657     0.27523602
    592     27816     0.249818446
    594     24209     0.217138707
    596     21099     0.189542484
    598     14048     0.126361656
    600     8058     0.072621641
    602     3537     0.031953522
    604     256     0.002178649
    606     3     0
根据表1中的数据，以剩余面积分数对温度做图，所得曲线称为半焦收缩曲线(如图1所示)；半焦收缩过程分为三个阶段，a：缓慢收缩阶段，b：过渡阶段，c：快速收缩阶段；根据半焦收缩曲线确定半焦收缩结束温度为606℃。依据半焦收缩曲线，利用一级反应动力学方程：
1 n ( ( - 1 n ( 1 - x ) ) / T 2 ) = 1 n [ AR φE ( 1 - 2 RT E ) ] - E RT
式中，1-x：样品的剩余面积分数；φ：升温速率，K·s-1；E：活化能，J·mol-1；A：频率因子，s-1，用上式左边对 作图，得一直线，斜率为
Figure A20061010220900113
求得活化能E；截距为 据此求得频率因子A，再依据k＝Aexp(-E/RT)求得反应速率常数k，s-1。神东惰质组颗粒在半焦收缩过程中三个阶段动力学参数的计算结果见表2：
表2 5K/min时神东惰质组半焦收缩过程的动力学参数
样品 E/kJ·mol-1 A/s-1 k/(10-3s-1)
a b c a b c a b c
神东惰质组 42.2 69.4 184.0 0.587  45.95 1.8E+9 0.98 1.69 4.11
具体实施方式II
将装有神东原煤颗粒的铂坩埚置于显微热台中，升温速率为10K/min，其它同具体实施方式I。该样品的半焦收缩结束温度为650℃，在半焦收缩过程中三个阶段动力学参数的计算结果见表3：
表3 10K/min时神东原煤半焦收缩过程的动力学参数
样品     E/kJ·mol-1   A/s-1  k/(10-3s-1)
    a   b   c   a  b  c  a  b  c
    神东原煤     20.4   35.2   155.4   0.013  0.136  1.1E+7  0.46  0.70  2.56
具体实施方式III
将装有平朔镜质组颗粒的铂坩埚置于显微热台中，升温速率为20K/min，其它同具体实施方式I。该样品的半焦收缩结束温度为790℃，在半焦收缩过程中三个阶段动力学参数的计算结果见表4：
表4 20K/min时平朔镜质组半焦收缩过程的动力学参数
样品 E/kJ·mol-1  A/s-1 k/(10-3s-1)
a     b    c  a  b  c  a  b  c
平朔镜质组 29.7     66.0    242.1  0.045  9.86  2.6E+10  0.44  1.11  14.64
具体实施方式IV
将装有平朔原煤颗粒的铂坩埚置于显微热台中，升温速率为40K/min，其它同具体实施方式I。该样品的半焦收缩结束温度为870℃，在半焦收缩过程中三个阶段动力学参数的计算结果见表5：
表5 40K/min时平朔原煤半焦收缩过程的动力学参数
样品   E/kJ·mol-1   A/s-1   k/(10-3s-1)
  a   b   c   a   b   c   a   b  c
  平朔原煤   40.9   83.9   580.1   0.124   24.81   5.8E+21   0.67   2.23  6.0E-3

Claims (1)

1、一种研究煤热解过程中半焦收缩特性的方法，其特征在于是一种以显微镜和热台为主要实验手段，以煤及其显微组分颗粒为研究对象，通过显微镜下拍摄到的显微图片获得样品在热解过程中的动态变化特点，用图像分析方法准确得到煤热解过程中半焦收缩特性的方法，其具体的方法步骤为：
I、取样及装置的调试
用毛细管取样，以保证对各样品的研究具有可比性将煤样横放在样品池中，各样品具有相同的厚度，将装有小于0.1mg样品颗粒的铂坩埚置于显微热台中，调节显微镜光强度为10，将热台通入流量为120mL/min的N2，接通循环冷却水，开启升温控温程序，在升温速率为5～40K/min时自室温加热至终温；
II、使用图像分析软件绘制半焦收缩曲线
利用徕卡显微镜自身配置的Qwin图像处理软件对热分解过程不同温度时采集的显微图片进行分析；由该软件控制的数字摄像机的拍摄条件为：曝光时间51.5ms，增益2.1，白平衡49、2.90、16；数码摄像机采集不同温度下样品颗粒的显微图片，用Qwin图像处理软件将起始温度时显微图片上煤样的灰度作为基准，灰度域值为1～10，求得热解过程中某一温度时样品颗粒具有此灰度域值的面积数值，从而算出该温度下的样品颗粒面积与起始时煤粒面积的比值，称为“剩余面积分数”，以剩余面积分数对温度做图，所得曲线称为半焦收缩曲线；
III、由半焦收缩曲线考察半焦收缩特性
通过观察半焦收缩曲线，发现该曲线能够反映热解过程中半焦面积随温度的变化情况；在升温过程中半焦收缩的程度随温度而变化，半焦收缩过程分为三个阶段，半焦开始收缩时收缩程度较小，称为缓慢收缩阶段；随温度升高收缩程度变大，称为过渡阶段；最后收缩程度急剧增大，称为快速收缩阶段，根据半焦收缩曲线确定在显微镜下观察到的半焦收缩程度最大时的温度即结束温度，结束温度越高，表明煤样半焦收缩的性能越弱，半焦收缩的程度越小，反之，煤样的半焦收缩的程度越大；考察升温速率分别为5，10，20和40K/min的影响，发现结束温度皆随升温速率的增大而变高；变质程度较高的煤及其显微组分的结束温度高于变质程度较低的煤及对应显微组分；对于原煤和显微组分的结束温度高低顺序为：惰质组＞原煤＞镜质组；
IV、依据半焦收缩曲线求得半焦收缩的动力学参数
依据半焦收缩曲线，利用一级反应动力学方程：
ln ( ( - ln ( 1 - x ) ) / T 2 ) = ln [ AR Φ E ( 1 - 2 RT E ) ] - E RT
式中，1-x：样品的剩余面积分数；φ：升温速率，K·s-1；E：活化能，J·mol-1；A：频率因子，s-1，用上式左边对 作图，得一直线，斜率为
Figure A2006101022090003C3
求得活化能E；截距为 ln [ AR ΦE ( 1 - 2 RT E ) ] , 据此求得频率因子A，再依据k＝Aexp(-E/RT)求得反应速率常数k，s-1，各样品的活化能、频率因子和反应速率常数皆随缓慢收缩阶段、过渡阶段和快速收缩阶段依次增大。
CNB2006101022094A 2006-12-01 2006-12-01 研究煤热解过程中半焦收缩特性的方法 Expired - Fee Related CN100395536C (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2006101022094A CN100395536C (zh) 2006-12-01 2006-12-01 研究煤热解过程中半焦收缩特性的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2006101022094A CN100395536C (zh) 2006-12-01 2006-12-01 研究煤热解过程中半焦收缩特性的方法

Publications (2)

Publication Number Publication Date
CN1963453A true CN1963453A (zh) 2007-05-16
CN100395536C CN100395536C (zh) 2008-06-18

Family

ID=38082594

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2006101022094A Expired - Fee Related CN100395536C (zh) 2006-12-01 2006-12-01 研究煤热解过程中半焦收缩特性的方法

Country Status (1)

Country Link
CN (1) CN100395536C (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102507493A (zh) * 2011-10-24 2012-06-20 上海理工大学 一种测定煤中的镜质组与惰质组在热解过程中相互作用的方法
CN105842386A (zh) * 2016-05-31 2016-08-10 扬州大学 消除燃煤污染物排放的两段分离式试验系统
CN108519301A (zh) * 2018-03-12 2018-09-11 沈阳环境科学研究院 一种利用热重分析仪评价煤及生物质焦反应性的方法
CN109580421A (zh) * 2019-01-03 2019-04-05 东北大学 难热解材料的挥发分释放特性指数计算方法
CN111693558A (zh) * 2019-03-11 2020-09-22 上海梅山钢铁股份有限公司 一种炼焦过程中半焦收缩率的检测方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2553070B2 (ja) * 1987-03-20 1996-11-13 興亜石油株式会社 か焼コ−クスの製造方法
JP2001324460A (ja) * 2000-05-12 2001-11-22 Nippon Steel Corp 石炭の品質評価方法
JP3892666B2 (ja) * 2000-12-25 2007-03-14 新日本製鐵株式会社 石炭の品質評価方法

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102507493A (zh) * 2011-10-24 2012-06-20 上海理工大学 一种测定煤中的镜质组与惰质组在热解过程中相互作用的方法
CN102507493B (zh) * 2011-10-24 2013-12-25 上海理工大学 一种测定煤中的镜质组与惰质组在热解过程中相互作用的方法
CN105842386A (zh) * 2016-05-31 2016-08-10 扬州大学 消除燃煤污染物排放的两段分离式试验系统
CN105842386B (zh) * 2016-05-31 2018-03-27 扬州大学 消除燃煤污染物排放的两段分离式试验系统
CN108519301A (zh) * 2018-03-12 2018-09-11 沈阳环境科学研究院 一种利用热重分析仪评价煤及生物质焦反应性的方法
CN109580421A (zh) * 2019-01-03 2019-04-05 东北大学 难热解材料的挥发分释放特性指数计算方法
CN111693558A (zh) * 2019-03-11 2020-09-22 上海梅山钢铁股份有限公司 一种炼焦过程中半焦收缩率的检测方法

Also Published As

Publication number Publication date
CN100395536C (zh) 2008-06-18

Similar Documents

Publication Publication Date Title
CN100395536C (zh) 研究煤热解过程中半焦收缩特性的方法
Soria-Verdugo et al. Comparison of wood pyrolysis kinetic data derived from thermogravimetric experiments by model-fitting and model-free methods
CN102156147B (zh) 用热分析方法检测聚丙烯在双向拉伸中的成膜稳定性
CN110455803A (zh) 一种文物鉴定系统及鉴定方法
CN103811701B (zh) 锂离子电池隔膜的材料一致性快速评价方法
Boddapati et al. Gel time prediction of multifunctional acrylates using a kinetics model
CN106383934B (zh) 一种焦炭质量预测模型的构建方法和炼焦配煤方法
CN105606641A (zh) 一种在线实时监测锂电池隔膜热收缩率系统及监测方法
CN101603819B (zh) 一种木材形变微观结构特征实时测量方法
CN105838407A (zh) 一种控制配合煤收缩度的炼焦煤配用方法
CN109102035B (zh) 一种基于聚类分析的炼焦煤多维指标相似性细化分类方法
Chen et al. Chemical structure transformation during the later stage of plastic layers during coking using Synchrotron infrared microspectroscopy technique
CN101750279B (zh) 氟碳彩涂板涂层中pvdf树脂的含量红外峰面积比值测定法
WO2023102893A1 (zh) 一种硅橡胶复合绝缘子老化状态评估方法及装置
CN107807222A (zh) 一种确定生排烃实验热成熟和生烃进程的方法和装置
Pozzobon et al. Biomass gasification under high solar heat flux: Experiments on thermally thick samples
CN105158177B (zh) 高光谱图像技术定量检测固态发酵水分分布均匀性的方法
CN109060878A (zh) 基于燃烧动力学参数的沥青组分热稳定性评价方法
Thompson et al. Estimating the temperature experienced by biomass particles during fast pyrolysis using microscopic analysis of biochars
CN105954285A (zh) 一种快速测定催化裂化油浆固含量的方法
CN106248612A (zh) 乙丙丁三元共聚聚丙烯中共聚单体含量的测定方法
CN105784545B (zh) 一种煤的膨胀与收缩性能评价方法
CN105866065B (zh) 一种乌洛托品-醋酸溶液中乌洛托品含量分析方法
CN113008931A (zh) 一种材料热分析动力学参数的获取方法
Pin et al. Exploration of Polymer Calorimetric Glass Transition Phenomenology by Two-Dimensional Correlation Analysis

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20080618

Termination date: 20100104