CN101419287B - 海雾厚度和低云云底高度的实时获取方法 - Google Patents

海雾厚度和低云云底高度的实时获取方法 Download PDF

Info

Publication number
CN101419287B
CN101419287B CN2008102380940A CN200810238094A CN101419287B CN 101419287 B CN101419287 B CN 101419287B CN 2008102380940 A CN2008102380940 A CN 2008102380940A CN 200810238094 A CN200810238094 A CN 200810238094A CN 101419287 B CN101419287 B CN 101419287B
Authority
CN
China
Prior art keywords
layer
thickness
height
sea fog
quasi saturation
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 - Fee Related
Application number
CN2008102380940A
Other languages
English (en)
Other versions
CN101419287A (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.)
Ocean University of China
Original Assignee
Ocean University of China
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 Ocean University of China filed Critical Ocean University of China
Priority to CN2008102380940A priority Critical patent/CN101419287B/zh
Publication of CN101419287A publication Critical patent/CN101419287A/zh
Application granted granted Critical
Publication of CN101419287B publication Critical patent/CN101419287B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Control Or Security For Electrophotography (AREA)
  • Laminated Bodies (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及一种海雾厚度和低云云底高度的实时获取方法。本发明充分利用气象台L波段二次测风雷达和数字式探空仪资料,依次分析了海雾和低云在垂直方向上温度、湿度的分布特征,计算了低层大气有海雾和有低云时的稳定度条件和湍流条件,统计分析归纳出将海雾和低云分离,然后再对海雾厚度和低云云底高度进行计算,获取海雾厚度和低云云底高度。形成的计算机软件将结果自动显示。本发明提供了实时监测的有效方法,即同时建立了一种客观可行的低云云底高度实时自动提取,和对海雾厚度进行实时监测的方法,其实施性强,可信度高,为海面上空飞行安全、海上交通运输、海港和沿海空港作业提供了有效气象资料。

Description

海雾厚度和低云云底高度的实时获取方法
技术领域
本发明属于海洋环境监测技术,具体涉及一种海雾厚度和低云云底高度的实时获取方法。
背景技术
有海雾时,水平能见度不足1000米,严重影响海上交通运输、海港和沿海空港作业。研究表明,海雾垂直厚度越厚,雾性就越浓,能见度就越差,就越不容易消散。所以海雾厚度是海雾监测的重要内容。但目前还没有对海雾厚度进行实时监测的有效方法。
另外,在中纬度的海面上空,经常出现低云。由于低云云底高度(距离海面的高度)一般在几百米,给飞机起飞、降落造成很大困难。因此,需要对云底高度做出足够准确的判断。虽然,云底高度可以通过云幕仪(ceilometer)测得,但是云幕仪价格昂贵,我国目前还不能生产,而气象人员一般通过目测对低云云底的高度进行估计。人工目测估计低云云底高度精度难以保证,而且主观成分较大,易造成估计值偏高或者偏低,或者贻误起飞降落时机,或者对起飞降落带来隐患。因此,建立一种客观可行的低云云底高度实时获取方法是十分必要的。
目前,我国沿海气象部门陆续配备了新一代探空设备——GFE(L)1型二次测风雷达(工作频率为1675MHz,简称L波段测风雷达),它与GTS1型数字式电子探空仪相配合,能连续自动测定地面至对流层上层的风向、风速、气温、气压、湿度等气象要素,其垂直分层最多可以每升高7米获取一套数据,从地面开始就可以获取资料,探测精度较上一代无线电探空设备有明显提高。但是,迄今还没有利用L波段二次测风雷达和数字式电子探空仪的资料(以下简称L波段雷达资料)对低云云底高度进行实时自动提取。
发明内容
本发明的目的是提供一种海雾厚度和低云云底高度的实时获取方法,以弥补现有技术的不足。
本发明充分利用了L波段二次测风雷达和数字式探空仪资料垂直分辨率高的特点,考虑了4种出现海雾的和3种出现低云的大气边界层垂直条件,可以有效地将海雾和低云区分开,用温度、湿度和风向风速在垂直方向上每隔30m的实测资料计算出海雾厚度和低云云底高度。
本发明的具体技术路线是利用L波段二次测风雷达和数字式探空仪资料,分析海雾和低云在垂直方向上温度、湿度的分布特征,计算低层大气有海雾和有低云时的稳定度和湍流条件,通过大样本统计分析,归纳出将海雾和低云分开的方法,将海雾与低云分离,然后再对海雾厚度和低云云底高度进行计算,获取海雾厚度和低云云底高度。该方法已经形成计算机软件,读入实时资料后,计算结果在计算机上自动显示。
本发明的海雾厚度和低云云底高度的获取方法或步骤如下:
1、资料预处理:读入数字探空仪实时资料,解码并生成该时刻ASCII码探空数据文件,该文件包含了从地面到对流层上层每30m间隔的高度z、温度t、相对湿度RH、露点温度td、气压P、风向ddd、风速ff等,再由生成的探空数据文件中,读入1000m以下各高度层的上述数据,用常规的方法将风向ddd、风速ff分解为经向风分量v和纬向风分量u;
2、依据准饱和层判据:RH≥93%,或者温度露点差t-td≤1.2℃的点(层),判断地面至1000m高度是否有准饱和层和准饱和层的厚度;如果第一层是准饱和层,再判断准饱和层厚度是否大于600m,如果大于600m,,采用判则1确定云底高度,判则1:从第一层向上寻找,找到t-td≤0.5℃的高度为云底高度;
3、根据不同的海雾发生条件,分别判断是否是海雾,如果是海雾,再给出海雾的厚度:如果第一层是准饱和层,且饱和层厚度小于600m,则计算准饱和层之上是否有干层,干层的判则是:从RH≥90%的最高高度向上,在200m以内突然降至RH≤60%;如果有干层,则为海雾,且海雾的厚度就是准饱和层的厚度;
如果第一层是准饱和层,但1000m以下没有干层,则计算准饱和层以下60m至以上300m是否有逆温层,其判则如下:在准饱和层高度确定以后,计算准饱和层附近的温度直减率 γ = - Δt Δz ,如果γ≤0,则为逆温。如果逆温层厚度≥100m,则为海雾,且海雾的厚度就是准饱和层的厚度;
如果第一层是准饱和层,但不符合干层和逆温层的判则条件,则计算准饱和层以下60m至以上400m是否为湿绝热稳定层;方法如下:首先计算湿绝热温度直减率γm
γ m = - dT dz = g C pd + L V C pd dw s dz
其中g=9.8m·s-2为重力加速度,Cpd=1004(J·kg-1·K-1)为干空气的定压比热,Lv=2.5×106(J/kg)为水汽凝结潜热, w s &ap; 0.622 e s P 为饱和混合比, e s = E 0 10 at b + t 为温度t下的饱和水汽压,E0=6.1hPa为t=0℃时的饱和水汽压,a=7.5,b=273.7;如果满足0<γ≤γm,则为湿绝热稳定层;如果湿绝热稳定层的厚度≥200m,则为海雾,且海雾的厚度就是准饱和层的厚度;
如果第一层是准饱和层,但不符合干层和逆温层判则,也不符合湿绝热稳定层条件,则计算准饱和层以下60m至以上200m的Richardson数Ri:
Ri = g &theta; v &OverBar; = &PartialD; &theta; v / &PartialD; z ( &PartialD; u / &PartialD; z ) 2 + ( &PartialD; v / &PartialD; z ) 2
其中 &theta; v = T v ( 1000 P ) R / C P 为虚位温,Tv=T(1+0.608q)为虚温(K), q = 0.622 e P - 0.378 e 为比湿(g/kg),Cp=Cpd(1+0.86q)为湿空气的定压比热, e = E 0 10 a t d b + t d 为水汽压,td为露点温度(℃)R=Rd(1+0.608q)为湿空气比气体常数,u,v分别为纬向风和经向风分量(m.s-1),g为重力加速度(m.s-2);如果至少有2层为Ri≤0,则为海雾,且海雾的厚度就是准饱和层的厚度;
如果第一层是准饱和层,但不符合以上干层、逆温层、湿绝热稳定层和Ri≤0的条件,则判断为低云,用判则1确定云底高度;
4、根据下面的湿度条件,得到低云云底高度:如果第一层没有准饱和层,则从地面向上至1000m寻找有无湿层,如果有湿层,再判断湿层厚度是否在300m以上,如果符合,则判断为低云,用判则2确定云底高度。低云云低高度判则2:从第一层向上至1000m,计算
Figure G2008102380940D00041
,最大
Figure G2008102380940D00042
出现的高度为云底高度;
5、如果从地面向上至1000m不符合湿层的条件,则在1000m以下没有低云,或者云底高度在1000m以上;将以上判识和计算结果直接在计算机上显示即得。
本发明的特点和优点,本发明同时建立了一种客观可行的低云云底高度实时自动提取,和对海雾厚度进行实时监测,其实施性强,可信度高,为海面上空飞行安全、海上交通运输、海港和沿海空港作业提供了有效气象资料。
附图说明
图1是本发明的海雾厚度和低云云底高度获取方法的流程示意图。
具体实施方式
如图1,本发明的海雾厚度和低云云底高度的获取方法或步骤如下:
1、资料预处理。读入数字探空仪实时资料,解码并生成该时刻ASCII码探空数据文件,该文件包含了从地面到对流层上层每30m间隔的高度z、温度t、相对湿度RH、露点温度td、气压P、风向ddd、风速ff。再由生成的探空数据文件中,读入1000m以下各高度层的上述数据,用常规的方法将风向ddd、风速ff分解为经向风分量v和纬向风分量u。
2、判断地面至1000m高度是否有准饱和层和准饱和层的厚度。依据准饱和层判据:RH≥93%,或者温度露点差t-td≤1.2℃的点(层)。如果第一层(距离地面30米高度处)是准饱和层,再判断准饱和层厚度是否大于600m,如果大于600m,则是低云。采用判则1确定云底高度。低云云低高度判则1:从第一层向上寻找,找到t-td≤0.5℃的高度为云底高度。
3、根据下面不同的海雾发生条件,分别判断是否是海雾,如果是海雾,再给出海雾的厚度:如果第一层是准饱和层,且饱和层厚度小于600m,则计算准饱和层之上是否有干层,干层的判则是:从RH≥90%的最高高度向上,在200m以内突然降至RH≤60%。如果有干层,则为海雾,且海雾的厚度就是准饱和层的厚度。
如果第一层是准饱和层,但1000m以下没有干层,则计算准饱和层附近(准饱和层以下60m至以上300m)是否有逆温层,其判则如下:在准饱和层高度确定以后,计算准饱和层附近的温度直减率 &gamma; = - &Delta;t &Delta;z ,如果γ≤0,则为逆温。如果逆温层厚度≥100m,则为海雾,且海雾的厚度就是准饱和层的厚度。
如果第一层是准饱和层,但不符合干层和逆温层的判则条件,则计算准饱和层附近(准饱和层以下60m至以上400m)是否为湿绝热稳定层。方法如下:首先计算湿绝热温度直减率γm
&gamma; m = - dT dz = g C pd + L V C pd dw s dz
其中g=9.8m·s-2为重力加速度,Cpd=1004(J·kg-1·K-1)为干空气的定压比热,Lv=2.5×106(J/kg)为水汽凝结潜热, w s &ap; 0.622 e s P 为饱和混合比, e s = E 0 10 at b + t 为温度t下的饱和水汽压,E0=6.1hPa为t=0℃时的饱和水汽压,a=7.5,b=273.7。如果满足0<γ≤γm,则为湿绝热稳定(含湿绝热中性)层。如果湿绝热稳定层的厚度≥200m,则为海雾,且海雾的厚度就是准饱和层的厚度。
如果第一层是准饱和层,但不符合干层和逆温层判则,也不符合湿绝热稳定层条件,则计算准饱和层附近(准饱和层以下60m至以上200m)的Richardson数Ri:
Ri = g &theta; v &OverBar; = &PartialD; &theta; v / &PartialD; z ( &PartialD; u / &PartialD; z ) 2 + ( &PartialD; v / &PartialD; z ) 2
其中 &theta; v = T v ( 1000 P ) R / C P 为虚位温,Tv=T(1+0.608q)为虚温(K), q = 0.622 e P - 0.378 e 为比湿(g/kg),Cp=Cpd(1+0.86q)为湿空气的定压比热, e = E 0 10 a t d b + t d 为水汽压,td为露点温度(℃)R=Rd(1+0.608q)为湿空气比气体常数,u,v分别为纬向风和经向风分量(m.s-1),g为重力加速度(m.s-2)。如果至少有2层为Ri≤0,则为海雾,且海雾的厚度就是准饱和层的厚度。
如果第一层是准饱和层,但不符合以上干层、逆温层、湿绝热稳定层和Ri≤0的条件,则判断为低云,用判则1确定云底高度。
4、根据下面的湿度条件,得到低云云底高度:如果第一层没有准饱和层,则从地面向上至1000m寻找有无湿层(湿层的条件是RH≥85%)。如果有湿层,再判断湿层厚度是否在300m以上,如果符合,则判断为低云。用判则2确定云底高度。低云云低高度判则2:从第一层向上至1000m,计算
Figure G2008102380940D00065
,最大
Figure G2008102380940D00066
出现的高度为云底高度。
5、如果从地面向上至1000m不符合湿层的条件,则在1000m以下没有低云,或者云底高度在1000m以上。
6、将以上判识和计算结果直接在计算机上显示。
不难将以上方法编制成计算机软件,探空完成后,启动该程序,就可以在计算机终端上显示出海雾厚度和云底高度,或者显示出云底高度在1000米以上。

Claims (2)

1.一种海雾厚度和低云云底高度的实时获取方法,其步骤如下:
(1)、资料预处理:读入数字探空仪实时资料,解码并生成该时刻ASCII码探空数据文件,该文件包含了从地面到对流层上层每30m间隔的高度z、温度t、相对湿度RH、露点温度td、气压P、风向ddd、风速ff,再由生成的探空数据文件中,读入1000m以下各高度层的上述数据,用常规的方法将风向ddd、风速ff分解为经向风分量v和纬向风分量u;
(2)、依据准饱和层判据:相对湿度RH≥93%的点,判断地面至1000m高度是否有准饱和层和准饱和层的厚度;如果第一层是准饱和层,再判断准饱和层厚度是否大于600m,如果大于600m,采用判则1确定云底高度,判则1:从第一层向上寻找,找到t-td≤0.5℃的高度为云底高度;
(3)、根据不同的海雾发生条件,分别判断是否是海雾,如果是海雾,再给出海雾的厚度:如果第一层是准饱和层,且饱和层厚度小于600m,则计算准饱和层之上是否有干层,干层的判则是:从RH≥90%的最高高度向上,在200m以内突然降至RH≤60%;如果有干层,则为海雾,且海雾的厚度就是准饱和层的厚度;
如果第一层是准饱和层,但1000m以下没有干层,则计算准饱和层以下60m至以上300m是否有逆温层,其判则如下:在准饱和层高度确定以后,计算准饱和层附近的温度直减率 
Figure FSB00000289967400011
如果γ≤0,则为逆温;如果逆温 层厚度≥100m,则为海雾,且海雾的厚度就是准饱和层的厚度;
如果第一层是准饱和层,但不符合干层和逆温层的判则条件,则计算准饱和层以下60m至以上400m是否为湿绝热稳定层;具体方法如下:首先计算湿绝热温度直减率γm
其中g=9.8m·s-2为重力加速度,Cpd=1004(J·kg-1·K-1)为干空气的定压比热,Lv=2.5×106(J/kg)为水汽凝结潜热, 
Figure FSB00000289967400022
为饱和混合比, 
Figure FSB00000289967400023
为温度t下的饱和水汽压,E0=6.1hPa为t=0℃时的饱和水汽压,a=7.5,b=273.7;如果满足0<γ≤γm,则为湿绝热稳定层;如果湿绝热稳定层的厚度≥200m,则为海雾,且海雾的厚度就是准饱和层的厚度;
如果第一层是准饱和层,但不符合干层和逆温层判则,也不符合湿绝热稳定层条件,则计算准饱和层以下60m至以上200m的Richardson数Ri:
Figure FSB00000289967400024
其中 
Figure FSB00000289967400025
为虚位温,Tv=T(1+0.608q)为虚温,其单位为K, 
Figure FSB00000289967400026
为比湿,其单位为g/kg,Cp=Cpd(1+0.86q)为湿空气的定压比热, 
Figure FSB00000289967400027
为水汽压,其中E0=6.1hPa;td为露点温度,其单位为℃,R=Rd(1+0.608q)为湿空气比气体常数,u,v分别为纬向风和经向风分量,其单位为m.s-1,g为重力加速度,其单位为m.s-2;如果至少有2层为Ri≤0, 则为海雾,且海雾的厚度就是准饱和层的厚度;
如果第一层是准饱和层,但不符合以上干层、逆温层、湿绝热稳定层和Ri≤0的条件,则判断为低云,用判则1确定云底高度;
(4)、根据湿度条件,得到低云云底高度:如果第一层没有准饱和层,则从地面向上至1000m寻找有无湿层,如果有湿层,再判断湿层厚度是否在300m以上,如果符合,则判断为低云,用判则2确定云底高度。低云云底高度判则2:从第一层向上至1000m,计算 最大 
Figure FSB00000289967400032
出现的高度为云底高度;
(5)、如果从地面向上至1000m不符合湿层的条件,则云底高度在1000m以上;将以上判识和计算结果直接在计算机上显示即得。
2.如权利要求1所述的海雾厚度和低云云底高度的实时获取方法,其特征是步骤(2)中判据准饱和层的依据是温度露点差t-td≤1.2℃的点。 
CN2008102380940A 2008-12-08 2008-12-08 海雾厚度和低云云底高度的实时获取方法 Expired - Fee Related CN101419287B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2008102380940A CN101419287B (zh) 2008-12-08 2008-12-08 海雾厚度和低云云底高度的实时获取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2008102380940A CN101419287B (zh) 2008-12-08 2008-12-08 海雾厚度和低云云底高度的实时获取方法

Publications (2)

Publication Number Publication Date
CN101419287A CN101419287A (zh) 2009-04-29
CN101419287B true CN101419287B (zh) 2010-12-22

Family

ID=40630172

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2008102380940A Expired - Fee Related CN101419287B (zh) 2008-12-08 2008-12-08 海雾厚度和低云云底高度的实时获取方法

Country Status (1)

Country Link
CN (1) CN101419287B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102230796B (zh) * 2011-04-06 2013-03-27 中国人民解放军理工大学气象学院 基于无线电掩星探测资料的云底高反演方法
CN105091851B (zh) * 2015-04-24 2018-07-17 广东小天才科技有限公司 一种测量高度的方法及装置
CN106772385B (zh) * 2016-11-28 2019-03-15 北京无线电测量研究所 一种基于毫米波雷达的港口航道海雾探测方法及系统
CN109657188B (zh) * 2018-11-23 2020-09-01 中国气象科学研究院 一种基于单目标变量判断云区的方法
CN112213703B (zh) * 2020-08-18 2023-07-18 成都信息工程大学 一种利用毫米波云雷达反演云参数的方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0578544B1 (fr) * 1992-07-09 1997-09-10 Rhea Procédé et dispositif pour l'évaluation des précipitations sur une zone de terrain
CN101008673A (zh) * 2007-01-24 2007-08-01 南京大桥机器有限公司 二次测风雷达接收控制、数据处理方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0578544B1 (fr) * 1992-07-09 1997-09-10 Rhea Procédé et dispositif pour l'évaluation des précipitations sur une zone de terrain
CN101008673A (zh) * 2007-01-24 2007-08-01 南京大桥机器有限公司 二次测风雷达接收控制、数据处理方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JP昭59-204786A 1984.11.20
姚雯.L波段电子探空仪相对湿度误差研究及其应用.《应用气象学报》.2008,第19卷(第3期),第356~361页. *

Also Published As

Publication number Publication date
CN101419287A (zh) 2009-04-29

Similar Documents

Publication Publication Date Title
Vickery et al. A hurricane boundary layer and wind field model for use in engineering applications
Chen et al. A modeling case study of heavy rainstorms along the Mei-Yu front
CN105607063B (zh) 一种机场低空风切变探测方法及系统
Emeis et al. Measurement and simulation of the 16/17 April 2010 Eyjafjallajókull volcanic ash layer dispersion in the northern Alpine region
Frey et al. The impact of overshooting deep convection on local transport and mixing in the tropical upper troposphere/lower stratosphere (UTLS)
Solari Thunderstorm downbursts and wind loading of structures: Progress and prospect
Johnson et al. An overview of the Lagrangian experiments undertaken during the North Atlantic regional Aerosol Characterisation Experiment (ACE-2)
CN101419287B (zh) 海雾厚度和低云云底高度的实时获取方法
De Wekker et al. A comparison of aerosol-layer and convective boundary-layer structure over a mountain range during STAAARTE'97
Hon et al. Application of LIDAR‐derived eddy dissipation rate profiles in low‐level wind shear and turbulence alerts at H ong K ong I nternational A irport
Viana et al. Propagation and effects of a mesoscale gravity wave over a weakly-stratified nocturnal boundary layer during the SABLES2006 field campaign
Chubb et al. Observations of high droplet number concentrations in Southern Ocean boundary layer clouds
Martínez-Alvarado et al. Sting jets in simulations of a real cyclone by two mesoscale models
Davenport et al. Creating a more realistic idealized supercell thunderstorm evolution via incorporation of base-state environmental variability
Soler et al. Modelling local sea-breeze flow and associated dispersion patterns over a coastal area in north-east Spain: a case study
Scully et al. Characterization and modulation of Langmuir circulation in Chesapeake Bay
Stith et al. NO signatures from lightning flashes
Eckermann et al. Mountain wave–induced polar stratospheric cloud forecasts for aircraft science flights during SOLVE/THESEO 2000
Juliano et al. An atmospheric hydraulic jump in the Santa Barbara Channel
Hallgren et al. Classification and properties of non-idealized coastal wind profiles-an observational study
Pinsky et al. Convective and turbulent motions in nonprecipitating Cu. Part II: LES simulated cloud represented by a starting plume
Frech et al. High-resolution weather database for the terminal area of Frankfurt Airport
Fesquet et al. Impact of terrain heterogeneity on coherent structure properties: numerical approach
Finn et al. Plume dispersion in low-wind-speed conditions during project Sagebrush phase 2, with emphasis on concentration variability
Ramana et al. Thermodynamic structure of the Atmospheric Boundary Layer over the Arabian Sea and the Indian Ocean during pre-INDOEX and INDOEX-FFP campaigns

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: 20101222

Termination date: 20111208