CN104240034A - 一种地震导致滑坡所造成的经济损失的快速评估方法 - Google Patents
一种地震导致滑坡所造成的经济损失的快速评估方法 Download PDFInfo
- Publication number
- CN104240034A CN104240034A CN201410490475.3A CN201410490475A CN104240034A CN 104240034 A CN104240034 A CN 104240034A CN 201410490475 A CN201410490475 A CN 201410490475A CN 104240034 A CN104240034 A CN 104240034A
- Authority
- CN
- China
- Prior art keywords
- county
- earthquake
- data
- caused
- landslide
- 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
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种地震导致滑坡所造成的经济损失的快速评估方法,其首先将研究区域按地质类别划分为不同的地质区域,然后获得研究区域的坡度数据;并计算各地质区域的安全系数,又计算出各地质区域的极限加速度及地震动峰值加速度数据,之后计算可能造成各地质区域经济损失的滑坡危险性参数并以县为单位求得每个县滑坡危险性参数平均值;然后计算每个县的经济脆弱性参数,对于每个县所辖区域中地震烈极限度大于等于VI或地震动峰值加速度大于等于0.05g的区域,计算折减系数,最后计算在所述研究区域由地震引起的次生滑坡可能造成的经济损失数量本发明所提供的一种地震导致滑坡所造成的经济损失的快速评估方法,其可快速评估次生滑坡灾害可能造成的经济损失情况。
Description
技术领域
本发明涉及一种评估地震造成的损失的方法,特别是一种对地震导致的次生滑坡灾害可能造成的经济损失的快速评估方法。
背景技术
滑坡是斜坡上的岩土体,在重力和地表水,地下水等作用下,沿一定的滑坡面整体向下滑动的一种现象。强烈地震时,地震诱发的滑坡次生地质灾害,特别是在山地、丘陵地区,其经济损失和人员伤亡甚至比地震直接造成的还要大。在中国大陆地区,尤其是地形相对复杂的山区,地震造成的滑坡是最为常见且破坏力最强的次生地震灾害。
然而,地震造成的滑坡具有发生突然、机理复杂、运动形态多样的特征,预测比较困难。各国政府为减轻地震造成的滑坡灾害,都采取了积极有效的措施。其中,在地震应急体系中,特别提到“防御次生灾害”,通过加强次生灾害监测预警,组织专家开展险情排查、评估和除险加固等方式,提高政府应对震后滑坡灾害的应急指挥和救援能力,减少人民的生命财产损失。
在针对地震造成的滑坡应急指挥***中,地震造成的滑坡灾害快速评估是地震应急决策***的重要基础和科学依据。地震导致滑坡造成的损失,也就是地震导致滑坡造成的直接经济损失,是指地震导致的滑坡灾害造成的建筑物和其他工程结构、设施、设备、财物等破坏而引起的经济损失,其折算价值以整修、恢复重建或重置所需费用来表示。它不包括非实物财产,如货币、有价证券等损失。场地和文物古迹破坏不折算为经济损失,其主要描述地震导致滑坡造成的破坏状态。
如果能够快速预估出地震导致滑坡灾害造成的灾害财产损失情况,就可以为合理启动应急响应级别、快速有效调配应急资源和科学制定灾后重建规划提供参考数据;从而能够帮助决策者针对性的制定资源调配部署。
我国在“九五”期间初步建立了地震灾害快速评估模型,其中包括了地震财产损失评估模型。然而,在现有评估方法中,关于地震造成的滑坡灾害经济损失快速评估方法研究却寥寥无几,更没有针对地震造成的滑坡次生灾害可能造成经济损失的有效评估方法,特别是适于震后应急期间的快速评估方法。
发明内容
本发明要解决的技术问题是提供一种地震导致滑坡所造成的经济损失的快速评估方法,以减少或避免前面所提到的问题。
为解决上述技术问题,本发明提供了一种地震导致滑坡所造成的经济损失的快速评估方法,其用于在地震灾害发生后,快速评估研究区域内由地震引起的次生滑坡可能造成的经济损失数量C,其包括如下步骤,
步骤A,在计算机上利用现有地质图将所述研究区域按地质类别划分为不同的地质区域,并获得各所述地质区域的内聚力c、土容重r及内摩擦角参数;
步骤B,在计算机上根据DEM的30米数据,利用Arcgis软件获得所述研究区域的坡度数据;
步骤C,在计算机上根据公式
获得各所述地质区域的安全系数Fs,上式中β为坡度;
步骤D,在计算机上根据公式
ac=(FS-1)gsinβ
获得各所述地质区域的极限加速度ac;
步骤E,将地震时监测台站采集到的南北向aEW和东西向aNS两个峰值加速度数据,利用公式求得平均值,然后进行插值计算,获得各所述地质区域的地震动峰值加速度amax数据,并利用Arcgis软件获得各所述地质区域amax-ac的数据;
步骤F,在计算机上根据下列公式计算可能造成各地质区域经济损失概率的滑坡危险性参数P1数据,
当amax-ac≥0时,求得P1=0.05-0.12(amax-ac-0.446)2,
当上式P1<0时,则P1取值为0;
当amax-ac<0时,则P1=0,
根据行政区划数据,以县为单位,利用Arcgis软件求得每个县滑坡危险性参数平均值P1i,i表示所述研究区域的第i县;
步骤G,按照下列公式计算每个县的经济脆弱性参数,
P'2i=(2.8-3.3K)1.9-0.2M,其中,
上式中M为地震震级,Swd为该县地震时未倒坍房屋面积,S为该县房屋总面积,i表示所述研究区域的第i县;
步骤H,对于每个县所辖区域中地震烈度大于等于VI或地震动峰值加速度大于等于0.05g的区域,按照下列公式计算折减系数S折:
上式中,S总为所述县总面积,SVI为所述县地震烈度大于等于VI的区域或地震动峰值加速度大于等于0.05g的面积;
步骤I,根据下列公式计算在所述研究区域由地震引起的次生滑坡可能造成的经济损失数量C,
C=∑S折i×P1i×P’2i×Di
本发明所提供的一种地震导致滑坡所造成的经济损失的快速评估方法,其可在地震发生后,快速评估由地震可能引发的次生滑坡灾害位置以及次生滑坡灾害可能造成的经济损失情况,对选择启动应急预案的应急响应级别具有较高参考价值,为决策者指定灾民的营救和安置计划提供依据。
附图说明
以下附图仅旨在于对本发明做示意性说明和解释,并不限定本发明的范 围。其中,
图1显示的是根据本发明的一个具体实施例的一种地震导致滑坡所造成的经济损失的快速评估方法的研究区域的地质区域划分图;
图2为图1基础上的研究区域的坡度图;
图3为图2基础上的安全系数数据图;
图4为图3基础上的临界加速度图;
图5为图1基础上的地震动峰值加速度图;
图6为图4、5基础上的amax-ac差值图;
图7为图6基础上的地震导致滑坡危险性参数图;
图8为图7基础上的县级地震导致滑坡危险性参数平均值数据图。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解,现说明本发明的具体实施方式。
本发明提供了一种地震导致滑坡所造成的经济损失的快速评估方法,其用于在地震灾害发生后,快速评估研究区域内由地震引起的次生滑坡可能造成的经济损失数量C,其包括如下步骤,
步骤A,在计算机上利用现有地质图将所述研究区域按地质类别划分为不同的地质区域,并获得各所述地质区域的内聚力c、土容重r及内摩擦角参数;不同的地质类别对应的参数如下表所示;
类别 | 内聚力(MPa) | 土容重(kN/m3) | 内摩擦角(°) |
火山岩 | 0.22 | 27.5 | 40 |
变质岩 | 0.17 | 26.5 | 33 |
碎屑岩 | 0.1 | 25.5 | 24 |
碳酸盐岩 | 0.07 | 23.5 | 16 |
第四系松散堆积物 | 0.05 | 21.5 | 13 |
步骤B,在计算机上根据DEM的30米数据,利用Arcgis软件获得所述研究区域的坡度数据;
基于30米的DEM数字高程模型(Digital Elevation Model)可通过“地理空间数据云”的ASTER GDEM获取,并利用Arcgis软件获得所述研究区 域的坡度数据,
步骤C,在计算机上根据公式
获得各所述地质区域的安全系数Fs,上式中β为坡度;
根据步骤A、B所获得的数据,可在Arcgis软件中根据上述公式计算出各所述地质区域的分辨率为30米的各基础地面数据单元的安全系数Fs,在计算之前,可先将内摩擦角和坡度β先换算为同一单位,例如,Arcgis软件在计算三角函数时需要以弧度为单位输入数据,则在计算前可先将内摩擦角和坡度β均换算为弧度进行输入。
步骤D,在计算机上根据公式
ac=(FS-1)gsinβ
获得各所述地质区域的极限加速度ac;在计算得到步骤C的安全系数Fs数据后,可在Arcgis软件中根据上述公式进一步计算出各所述地质区域的分辨率为30米的各基础地面数据单元的极限加速度ac。
步骤E,将地震时监测台站采集到的南北向aEW和东西向aNS两个峰值加速度数据,利用公式求得平均值,然后进行插值计算,获得各所述地质区域的地震动峰值加速度amax数据,并利用Arcgis软件获得各所述地质区域amax-ac的数据;也就是说,可在Arcgis软件中进一步计算出各所述地质区域的分辨率为30米的各基础地面数据单元的amax-ac的数据。
步骤F,在计算机上,根据下列公式计算可能造成各地质区域经济损失概率的滑坡危险性参数P1数据,
当amax-ac≥0时,求得P1=0.05-0.12(amax-ac-0.446)2,
当上式P1<0时,则P1取值为0;
当amax-ac<0时,则P1=0,
根据行政区划数据,以县为单位,利用Arcgis软件求得每个县滑坡危险性参数平均值P1i,i表示所述研究区域的第i县;
利用步骤E所获得的分辨率为30米的各基础地面数据单元的amax-ac的数据,可以根据上述公式计算出各基础地面数据单元的经济损失概率P1的数值,然后可将行政区划数据导入Arcgis软件,从而可获得每个县所包含的分辨率为30米的各基础地面数据单元的经济损失概率P1的数值,并进一步可每个县的滑坡危险性参数平均值P1i,
步骤G,按照下列公式计算每个县的经济脆弱性参数,
2其中,
上式中M为地震震级,Swd为该县地震时未倒坍房屋面积,S为该县房屋总面积,i表示所述研究区域的第i县;
步骤H,对于每个县所辖区域中地震烈度大于等于VI或地震动峰值加速度大于等于0.05g的区域,按照下列公式计算折减系数S折:
上式中,S总为所述县总面积,SVI为所述县地震烈度大于等于VI或地震动峰值加速度大于等于0.05g的区域的面积;
步骤I,根据下列公式计算在所述研究区域由地震引起的次生滑坡可能造成的经济损失数量C,
C=∑S折i×P1i×P2i×Ci,
也即是通过公式
即可计算出地震引起的次生滑坡可能造成的经济损失数量C。
图1显示的是根据本发明的一个具体实施例的一种地震导致滑坡所造成的经济损失的快速评估方法的研究区域的地质区域划分图;图2为图1基础上的研究区域的坡度图;图3为图2基础上的安全系数数据图;图4为图3基础上的极限加速度图;图5为图1基础上的地震动峰值加速度图;图6为图4、5基础上的amax-ac差值图;图7为图6基础上的地震导致滑坡危险 性参数图;图8为图7基础上的县级地震导致滑坡危险性参数平均值数据图。
参见图1-8所示,下面以2013年07月22日甘肃岷县、漳县Ms6.6地震为例具体说明本发明所提供的一种地震导致滑坡所造成的经济损失的快速评估方法。
首先,如图1所示,在计算机上利用现有地质图将遭受地震的研究区域按地质类别划分为不同的地质区域,并获得各所述地质区域的内聚力c、土容重r及内摩擦角参数;
然后,如图2所示,在计算机上根据所述区域的DEM的30米数据,利用Arcgis软件获得所述研究区域的坡度数据;
之后,如图3所示,在计算机上根据公式
获得各所述地质区域的安全系数Fs
之后,如图4所示,在计算机上根据公式
ac=(FS-1)gsinβ
获得各所述地质区域的极限加速度ac;
之后,如图5、6所示,将地震时监测台站采集到的南北向aEW和东西向aNS两个峰值加速度数据,利用公式求得平均值,然后进行插值计算,获得各所述地质区域的地震动峰值加速度amax数据,并利用Arcgis软件获得各所述地质区域amax-ac的数据;也就是说,可在Arcgis软件中进一步计算出各所述地质区域的分辨率为30米的各基础地面数据单元的amax-ac的数据。
之后,如图7、8所示,在计算机上,根据下列公式计算可能造成各地质区域经济损失概率的滑坡危险性参数P1数据,
当amax-ac0时,求得P1=0.05-0.12(amax-ac-0.446)2,
当上式P1<0时,则P1取值为0;
当amax-ac<0时,则P1=0,
根据行政区划数据,以县为单位,利用Arcgis软件求得每个县滑坡危险性参数平均值P1i,
之后,按照下列公式计算每个县的经济脆弱性参数,
P'2i=(2.8-3.3K)1.9-0.2M,其中,
上式中M为地震震级,Swd为该县地震时未倒坍房屋面积,S为该县房屋总面积,i表示所述研究区域的第i县;
之后,对于每个县所辖区域中地震烈度大于等于VI或地震动峰值加速度大于等于0.05g的区域,按照下列公式计算折减系数S折:
上式中,S总为所述县总面积,SVI为所述县地震烈度大于等于VI或地震动峰值加速度大于等于0.05g的区域的面积;
最后,根据下列公式计算在所述研究区域由地震引起的次生滑坡可能造成的经济损失数量C,
C=∑S折i×P1i×P’2i×Di
也即是通过公式
即可计算出地震引起的次生滑坡可能造成的经济损失数量C。
在本具体实施例中,岷县、漳县的相关数据计算后如下表所示:
也即是,通过本发明提供的一种地震导致滑坡所造成的经济损失的快速评估方法,可估算出由于本次地震所导致的次生滑坡灾害可能造成的经济损失数为660.3万,而本次地震所导致的次生滑坡灾害实际造成的经济损失为720.46万,也就是说,通过本发明提供的一种地震导致滑坡所造成的经济损 失的快速评估方法所评估结果的误差小于10%,因此,本发明所提供的一种地震导致滑坡所造成的经济损失的快速评估方法可以用于未来地震灾害发生时,对地震导致的次生滑坡灾害的发生位置以及其可能造成的经济损失进行评估。
本发明所提供的一种地震导致滑坡所造成的经济损失的快速评估方法,其可在地震发生后,快速评估由地震可能引发的次生滑坡灾害位置以及次生滑坡灾害可能造成的经济损失情况,对选择启动应急预案的应急响应级别具有较高参考价值,为决策者指定灾民的营救和安置计划提供依据。
本领域技术人员应当理解,虽然本发明是按照多个实施例的方式进行描述的,但是并非每个实施例仅包含一个独立的技术方案。说明书中如此叙述仅仅是为了清楚起见,本领域技术人员应当将说明书作为一个整体加以理解,并将各实施例中所涉及的技术方案看作是可以相互组合成不同实施例的方式来理解本发明的保护范围。
以上所述仅为本发明示意性的具体实施方式,并非用以限定本发明的范围。任何本领域的技术人员,在不脱离本发明的构思和原则的前提下所作的等同变化、修改与结合,均应属于本发明保护的范围。
Claims (1)
1.一种地震导致滑坡所造成的经济损失的快速评估方法,其特征在于,其用于在地震灾害发生后,快速评估研究区域内由地震引起的次生滑坡可能造成的经济损失数量C,其包括如下步骤,
步骤A,在计算机上利用现有地质图将所述研究区域按地质类别划分为不同的地质区域,并获得各所述地质区域的内聚力c、土容重r及内摩擦角参数;
步骤B,在计算机上根据DEM的30米数据,利用Arcgis软件获得所述研究区域的坡度数据;
步骤C,在计算机上根据公式
获得各所述地质区域的安全系数Fs,上式中β为坡度;
步骤D,在计算机上根据公式
ac=(FS-1)gsinβ
获得各所述地质区域的极限加速度ac;
步骤E,将地震时监测台站采集到的南北向aEW和东西向aNS两个峰值加速度数据,利用公式求得平均值,然后进行插值计算,获得各所述地质区域的地震动峰值加速度amax数据,并利用Arcgis软件获得各所述地质区域amax-ac的数据;
步骤F,在计算机上根据下列公式计算可能造成各地质区域经济损失概率的滑坡危险性参数P1数据,
当amax-ac0时,求得P1=0.05-0.12(amax-ac-0.446)2,
当上式P1<0时,则P1取值为0;
当amax-ac<0时,则P1=0,
根据行政区划数据,以县为单位,利用Arcgis软件求得每个县滑坡危险性参数平均值P1i,i表示所述研究区域的第i县;
步骤G,按照下列公式计算每个县的经济脆弱性参数,
其中,
上式中M为地震震级,Swd为该县地震时未倒坍房屋面积,S为该县房屋总面积,i表示所述研究区域的第i县;
步骤H,对于每个县所辖区域中地震烈度大于等于VI或地震动峰值加速度大于等于0.05g的区域,按照下列公式计算折减系数S折:
上式中,S总为所述县总面积,SVI为所述县地震烈度大于等于VI或地震动峰值加速度大于等于0.05g的区域的面积;
步骤I,根据下列公式计算在所述研究区域由地震引起的次生滑坡可能造成的经济损失数量C,
C=∑S折i×P1i×P’2i×Di。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410490475.3A CN104240034B (zh) | 2014-09-23 | 2014-09-23 | 一种地震导致滑坡所造成的经济损失的快速评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410490475.3A CN104240034B (zh) | 2014-09-23 | 2014-09-23 | 一种地震导致滑坡所造成的经济损失的快速评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104240034A true CN104240034A (zh) | 2014-12-24 |
CN104240034B CN104240034B (zh) | 2017-03-29 |
Family
ID=52228055
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410490475.3A Expired - Fee Related CN104240034B (zh) | 2014-09-23 | 2014-09-23 | 一种地震导致滑坡所造成的经济损失的快速评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104240034B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114239108A (zh) * | 2021-12-17 | 2022-03-25 | 四川大学 | 基于监控物联网的城市建筑群震后损失分布计算方法 |
CN114662950A (zh) * | 2022-03-31 | 2022-06-24 | 中南大学 | 一种适用于区域地震滑坡危险性评估的斜坡单元细化方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102005104A (zh) * | 2009-09-02 | 2011-04-06 | 吴立新 | 滑坡***移和姿态远距离快速监测报警装置及方法 |
CN103364830A (zh) * | 2013-07-24 | 2013-10-23 | 北京师范大学 | 一种基于多因子的地震后发生崩滑灾害位置的预测方法 |
CN103605898A (zh) * | 2013-11-29 | 2014-02-26 | 民政部国家减灾中心 | 一种区域泥石流灾害风险估计方法 |
-
2014
- 2014-09-23 CN CN201410490475.3A patent/CN104240034B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102005104A (zh) * | 2009-09-02 | 2011-04-06 | 吴立新 | 滑坡***移和姿态远距离快速监测报警装置及方法 |
CN103364830A (zh) * | 2013-07-24 | 2013-10-23 | 北京师范大学 | 一种基于多因子的地震后发生崩滑灾害位置的预测方法 |
CN103605898A (zh) * | 2013-11-29 | 2014-02-26 | 民政部国家减灾中心 | 一种区域泥石流灾害风险估计方法 |
Non-Patent Citations (4)
Title |
---|
GODT J. ET AL.: "Rapid assessment of earthquake-induced landsliding", 《ASSOCIATION OF DISASTER PREVENTION RESEARCH》 * |
王涛: "汶川地震重灾区地质灾害危险性评估研究", 《中国博士学位论文全文数据库》 * |
王涛等: "基于简化Newmark位移模型的区域地震滑坡危险性快速评估-以汶川Ms8.0级地震为例", 《工程地质学报》 * |
高晓路: "中国农村房屋震灾脆弱性评估及其成因分析", 《地理学报》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114239108A (zh) * | 2021-12-17 | 2022-03-25 | 四川大学 | 基于监控物联网的城市建筑群震后损失分布计算方法 |
CN114239108B (zh) * | 2021-12-17 | 2022-09-13 | 四川大学 | 基于监控物联网的城市建筑群震后损失分布计算方法 |
CN114662950A (zh) * | 2022-03-31 | 2022-06-24 | 中南大学 | 一种适用于区域地震滑坡危险性评估的斜坡单元细化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104240034B (zh) | 2017-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Karimzadeh et al. | A GIS-based seismic hazard, building vulnerability and human loss assessment for the earthquake scenario in Tabriz | |
Hashemi et al. | A GIS-based earthquake damage assessment and settlement methodology | |
Del Gaudio et al. | An approach to time-probabilistic evaluation of seismically induced landslide hazard | |
Putra et al. | Seismic hazard analysis for Indonesia | |
Sitharam et al. | Comprehensive seismic hazard assessment of Tripura and Mizoram states | |
Bhatti et al. | Probabilistic seismic hazard analysis of Islamabad, Pakistan | |
CN104133996A (zh) | 一种基于云模型和数据场的地面沉降风险等级评估方法 | |
Wang et al. | A critique of probabilistic versus deterministic seismic hazard analysis with special reference to the New Madrid seismic zone | |
Erdik et al. | Assessment of seismic risk in tashkent, uzbekistan and bishkek, kyrgyz republic | |
CN104615847A (zh) | 一种基于概率方法的南海海啸危险性预测方法 | |
Mazumder et al. | Seismic damage assessment using RADIUS and GIS: A case study of Sylhet City, Bangladesh | |
Wang et al. | Incorporating triggering and environmental factors in the analysis of earthquake-induced landslide hazards | |
Ashadi et al. | Probabilistic Seismic‐Hazard Analysis for Central Java Province, Indonesia | |
CN104200126A (zh) | 一种地震导致滑坡所造成的人员死亡的快速评估方法 | |
Niño et al. | Seismic-event-based methodology to obtain earthquake-induced translational landslide regional hazard maps | |
Al-Dogom et al. | Geostatistical seismic analysis and hazard assessment; United Arab Emirates | |
Rivas-Medina et al. | The role of GIS in urban seismic risk studies: application to the city of Almería (southern Spain) | |
CN104240034A (zh) | 一种地震导致滑坡所造成的经济损失的快速评估方法 | |
Hancilar et al. | Earthquake risk assessment for the building inventory of Muscat, Sultanate of Oman | |
Sekac et al. | GIS based evaluation in earthquake hazard micro-zonation-a case study of madang and morobe province, Papua New Guinea | |
Bellalem et al. | Seismic risk assessment for the downtown of the city of Blida, Algeria | |
Petersen et al. | New USGS Southeast Asia seismic hazard maps | |
James et al. | Seismic zonations at micro and macro-level for regions in the peninsular India | |
Jaimes et al. | An earthquake-event-based method for mapping tsunami hazards | |
Kontoes et al. | Mapping seismic vulnerability and risk of cities: the MASSIVE project |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170329 Termination date: 20210923 |