CN107255837A - 一种二维地形正演和改正方法 - Google Patents
一种二维地形正演和改正方法 Download PDFInfo
- Publication number
- CN107255837A CN107255837A CN201710575229.1A CN201710575229A CN107255837A CN 107255837 A CN107255837 A CN 107255837A CN 201710575229 A CN201710575229 A CN 201710575229A CN 107255837 A CN107255837 A CN 107255837A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- value
- coordinate
- forward modeling
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明提供了一种二维地形正演和改正方法,在绘制视电阻率等值断面图中,对均匀介质的地形和非均匀介质的地形进行坐标点数据采集并获取其电阻率值,通过正演和改正绘制纯地形剖面的视电阻率等值断面图。本发明能正确地计算出均匀起伏地形剖面视电阻率异常的位置、形态和数值的大小,客观显现地形起伏在剖面形成的假异常,进行直接的数据运算,免去了按一定比例制作导电纸、电阻网络、水槽和土槽等模型来模拟野外情况的工艺难、制作苦,同时也免除了人为带来的偶然误差。
Description
技术领域
本发明涉及一种二维地形正演和改正方法,属于地理信息系统技术领域。
背景技术
“常规直流电法任意对称四极装置(含高密度电法相同装置)”在矿产地质、水文地质、工程地质等领域应用越来越广泛,视电阻率资料的解释是在水平(或单斜)半空间进行的。然而实际地形存在各种各样的起伏情况,理论与实践表明,地形不但可以引起假异常,而且会掩盖地下由矿体或目标地质体引起的真异常。如果按照水平(或单斜)地形来处理数据、解释,势必会引起一定的误差,甚至导致完全错误的结果。对勘探和解释工作造成误导。因此有必要对起伏地形状态下二维空间人工电场分布的影响进行研究和评估,用地形改正方式削弱或消除地形对视电阻率观测值的影响,使视电阻率异常客观反映地质体的真实存在。
早期,日本、前苏联及我国,乃至西方国家的地球物理学者曾采用导电纸模拟、电阻网络模拟、水槽模拟、土槽模拟等物理模拟方法研究地形影响问题,如按一定比例制作导电纸、电阻网络,水槽、土槽等模型来模拟野外情况,得到地形影响的定性认识,并在野外实际测量中加以识别。有些学者也应用模拟结果,对野外观测数据作些简单的计算和改正,试验将地形影响去掉。这些做法基本属定性解释,不能从根本上解决问题。
其后,也有学者采用一些解析计算方法研究地形影响,一些反演软件也开发了地形改正功能,但基于原理是以地形对静态电流分布方式对异常形态按地形进行扭曲改正,不是真正意义上的地形改正,对地形影响的电阻率数值大小并未作任何改为变,仅改变异常存在形态,从本质上未达到真正地形改正的目的,所得视电阻率异常成果并非地质体的客观存在方式,而是不能定性、半定量、定量的对异常体作出判断和解释。
发明内容
为解决上述技术问题,本发明提供了一种二维地形正演和改正方法,该二维地形正演和改正方法省去了繁琐的模拟改正工作、减少了工作环节、节约了经费、缩短了工期,同时可以丰富专业教科书的内容,从本质上还原地下地质体异常的真实状状态。
本发明通过以下技术方案得以实现。
本发明提供的一种二维地形正演和改正方法,在绘制视电阻率等值断面图中,对均匀介质的地形和非均匀介质的地形进行坐标点数据采集并获取其电阻率值,通过正演和改正绘制纯地形剖面的视电阻率等值断面图。
基于一种二维地形正演和改正方法,包括以下步骤:
①获取电阻率值:以直流电法任意对称四极装置和高密度电法相同装置所测视电阻率值的剖面视电阻率平均值作为均匀介质的地形剖面视的电阻率值,以直流电法任意对称四极装置和高密度电法相同装置所测视电阻率值为非均匀介质的地形剖面视的电阻率值;
②采集地形点坐标:对均匀介质的地形和非均匀介质的地形进行剖面地形坐标点数据采集;
③数据地形正演:根据步骤①和步骤②中的参数,通过正演对称四极地形函数和正演公式,获取正演视电阻率值;
④绘制正演线图:根据正演视电阻率值制出正演的纯地形剖面的视电阻率等值断面图;
⑤数据地形改正:根据步骤①、步骤②中的参数和对正演视电阻率值的改正,通过改正对称四极地形函数和改正公式,获取改正后视电阻率值;
⑥绘制改正线图:根据改正后视电阻率值绘制出改正的纯地形剖面的视电阻率等值断面图。
所述获取电阻率值,具体包括:
通过常规直流电法任意对称四极装置,含高密度电法相同装置,获取均匀介质的地形和非均匀介质的地形剖面视电阻率二维假异常的位置、形态和数值的大小,获取均匀介质的地形和非均匀介质的地形剖面视的电阻率值。
所述采集地形点坐标,具体包括:
通过建立坐标系,获取均匀介质的地形和非均匀介质的地形的坐标值,对均匀介质的地形和非均匀介质的地形的坐标值进行剖面地形点坐标数据采集。
所述正演对称四极地形函数为:
其中,f(x,y)为地形变化函数,x为坐标点相对1号电极,x=0的剖面横坐标,y为坐标点剖面纵坐标,A、B、M、N为均匀介质的地形坐标上的点,xA、yA为A点的坐标值,xB、yB为B点的坐标值,xM、yM为M点的坐标值,xN、yN为N点的坐标值。
所述正演公式为:
ρs=ρ0×[1+f(x,y)],
其中,f(x,y)为地形变化函数,ρs为地形正演视电阻率,ρ0为均匀介质剖面视电阻率值。
所述改正对称四极地形函数为:
其中,f(xi,yi)为地形变化函数,xi为坐标点点相对1号电极,x=0的剖面横坐标,yi为坐标点剖面纵坐标,iA、iB、iM、iN为均匀介质的地形坐标上的点,xiA、yiA为iA点的坐标值,xiB、yiB为iB点的坐标值,xiM、yiM为iM点的坐标值,xiN、yiN为iN点的坐标值。
所述改正公式为:
ρsi=F[ρ0i,f(xi,yi)]=ρsi=[log([ρ0i)×[1-f(xi,yi)]13/5,
其中,ρsi为地形改正后视电阻率,ρ0i为视电阻率实测值,f(xi,yi)为地形变化函数,xi为坐标点点相对1号电极,x=0的剖面横坐标,yi为坐标点剖面纵坐标。
本发明的有益效果在于:能正确地计算出均匀起伏地形剖面视电阻率异常的位置、形态和数值的大小,客观显现地形起伏在剖面形成的假异常,进行直接的数据运算,免去了按一定比例制作导电纸、电阻网络、水槽和土槽等模型来模拟野外情况的工艺难、制作苦,同时也免除了人为带来的偶然误差。
附图说明
图1是本发明的流程图。
具体实施方式
下面进一步描述本发明的技术方案,但要求保护的范围并不局限于所述。
如图1所示,一种二维地形正演和改正方法,在绘制视电阻率等值断面图中,对均匀介质的地形和非均匀介质的地形进行坐标点数据采集并获取其电阻率值,通过正演和改正绘制纯地形剖面的视电阻率等值断面图。
基于一种二维地形正演和改正方法,包括以下步骤:
①获取电阻率值:以直流电法任意对称四极装置和高密度电法相同装置所测视电阻率值的剖面视电阻率平均值作为均匀介质的地形剖面视的电阻率值,以直流电法任意对称四极装置和高密度电法相同装置所测视电阻率值为非均匀介质的地形剖面视的电阻率值;
②采集地形点坐标:对均匀介质的地形和非均匀介质的地形进行剖面地形坐标点数据采集;
③数据地形正演:根据步骤①和步骤②中的参数,通过正演对称四极地形函数和正演公式,获取正演视电阻率值;
④绘制正演线图:根据正演视电阻率值制出正演的纯地形剖面的视电阻率等值断面图,确定假异常的规模、形态、产状等;
⑤数据地形改正:根据步骤①、步骤②中的参数和对正演视电阻率值的改正,通过改正对称四极地形函数和改正公式,获取改正后视电阻率值;
⑥绘制改正线图:根据改正后视电阻率值绘制出改正的纯地形剖面的视电阻率等值断面图,确定异常的规模、强度、形态、产状等。
进一步地,所述获取电阻率值,具体包括:通过常规直流电法任意对称四极装置,含高密度电法相同装置,获取均匀介质的地形和非均匀介质的地形剖面视电阻率二维假异常的位置、形态和数值的大小,客观显现地形起伏在二维剖面上形成的假异常,获取均匀介质的地形和非均匀介质的地形剖面视的电阻率值。
进一步地,所述采集地形点坐标,具体包括:通过建立坐标系,获取均匀介质的地形和非均匀介质的地形的坐标值,对均匀介质的地形和非均匀介质的地形的坐标值进行剖面地形点坐标数据采集。
进一步地,所述正演对称四极地形函数为:其中,f(x,y)为地形变化函数,x为坐标点相对1号电极,x=0的剖面横坐标,y为坐标点剖面纵坐标,A、B、M、N为均匀介质的地形坐标上的点,xA、yA为A点的坐标值,xB、yB为B点的坐标值,xM、yM为M点的坐标值,xN、yN为N点的坐标值。
进一步地,所述正演公式为:ρs=ρ0×[1+f(x,y)],其中,f(x,y)为地形变化函数,ρs为地形正演视电阻率,ρ0为均匀介质剖面视电阻率值。
进一步地,所述改正对称四极地形函数为:其中,f(xi,yi)为地形变化函数,xi为坐标点点相对1号电极,x=0的剖面横坐标,yi为坐标点剖面纵坐标,iA、iB、iM、iN为均匀介质的地形坐标上的点,xiA、yiA为iA点的坐标值,xiB、yiB为iB点的坐标值,xiM、yiM为iM点的坐标值,xiN、yiN为iN点的坐标值。
进一步地,所述改正公式为:ρsi=F[ρ0i,f(xi,yi)]=ρsi=[log([ρ0i)×[1-f(xi,yi)]13/5,其中,ρsi为地形改正后视电阻率,ρ0i为视电阻率实测值,f(xi,yi)为地形变化函数,xi为坐标点点相对1号电极,x=0的剖面横坐标,yi为坐标点剖面纵坐标。
综上所述,本发明能正确地计算出均匀起伏地形剖面视电阻率异常的位置、形态和数值的大小,客观显现地形起伏在剖面形成的假异常,进行直接的数据运算。
Claims (8)
1.一种二维地形正演和改正方法,其特征在于:在绘制视电阻率等值断面图中,对均匀介质的地形和非均匀介质的地形进行坐标点数据采集并获取其电阻率值,通过正演和改正绘制纯地形剖面的视电阻率等值断面图。
2.如权利要求1所述的二维地形正演和改正方法,其特征在于:包括以下步骤:
①获取电阻率值:以直流电法任意对称四极装置和高密度电法相同装置所测视电阻率值的剖面视电阻率平均值作为均匀介质的地形剖面视的电阻率值,以直流电法任意对称四极装置和高密度电法相同装置所测视电阻率值为非均匀介质的地形剖面视的电阻率值;
②采集地形点坐标:对均匀介质的地形和非均匀介质的地形进行剖面地形坐标点数据采集;
③数据地形正演:根据步骤①和步骤②中的参数,通过正演对称四极地形函数和正演公式,获取正演视电阻率值;
④绘制正演线图:根据正演视电阻率值制出正演的纯地形剖面的视电阻率等值断面图;
⑤数据地形改正:根据步骤①、步骤②中的参数和对正演视电阻率值的改正,通过改正对称四极地形函数和改正公式,获取改正后视电阻率值;
⑥绘制改正线图:根据改正后视电阻率值绘制出改正的纯地形剖面的视电阻率等值断面图。
3.如权利要求2所述的二维地形正演和改正方法,其特征在于:所述获取电阻率值,具体包括:
通过常规直流电法任意对称四极装置,含高密度电法相同装置,获取均匀介质的地形和非均匀介质的地形剖面视电阻率二维假异常的位置、形态和数值的大小,获取均匀介质的地形和非均匀介质的地形剖面视的电阻率值。
4.如权利要求2所述的二维地形正演和改正方法,其特征在于:所述采集地形点坐标,具体包括:
通过建立坐标系,获取均匀介质的地形和非均匀介质的地形的坐标值,对均匀介质的地形和非均匀介质的地形的坐标值进行剖面地形点坐标数据采集。
5.如权利要求2所述的二维地形正演和改正方法,其特征在于:所述正演对称四极地形函数为:
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>y</mi>
<mi>A</mi>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mi>M</mi>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mi>N</mi>
</msub>
<mo>+</mo>
<msub>
<mi>y</mi>
<mi>B</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>x</mi>
<mi>M</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>N</mi>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
其中,f(x,y)为地形变化函数,x为坐标点相对1号电极,x=0的剖面横坐标,y为坐标点剖面纵坐标,A、B、M、N为均匀介质的地形坐标上的点,xA、yA为A点的坐标值,xB、yB为B点的坐标值,xM、yM为M点的坐标值,xN、yN为N点的坐标值。
6.如权利要求2所述的二维地形正演和改正方法,其特征在于:所述正演公式为:
ρs=ρ0×[1+f(x,y)],
其中,f(x,y)为地形变化函数,ρs为地形正演视电阻率,ρ0为均匀介质剖面视电阻率值。
7.如权利要求2所述的二维地形正演和改正方法,其特征在于:所述改正对称四极地形函数为:
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mi>A</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mi>M</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mi>N</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mi>B</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>M</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>N</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
其中,f(xi,yi)为地形变化函数,xi为坐标点点相对1号电极,x=0的剖面横坐标,yi为坐标点剖面纵坐标,iA、iB、iM、iN为均匀介质的地形坐标上的点,xiA、yiA为iA点的坐标值,xiB、yiB为iB点的坐标值,xiM、yiM为iM点的坐标值,xiN、yiN为iN点的坐标值。
8.如权利要求2所述的二维地形正演和改正方法,其特征在于:所述改正公式为:
ρsi=F[ρ0i,f(xi,yi)]=ρsi=[log([ρ0i)×[1-f(xi,yi)]13/5,
其中,ρsi为地形改正后视电阻率,ρ0i为视电阻率实测值,f(xi,yi)为地形变化函数,xi为坐标点点相对1号电极,x=0的剖面横坐标,yi为坐标点剖面纵坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710575229.1A CN107255837B (zh) | 2017-07-14 | 2017-07-14 | 一种二维地形正演和改正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710575229.1A CN107255837B (zh) | 2017-07-14 | 2017-07-14 | 一种二维地形正演和改正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107255837A true CN107255837A (zh) | 2017-10-17 |
CN107255837B CN107255837B (zh) | 2019-12-13 |
Family
ID=60026342
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710575229.1A Expired - Fee Related CN107255837B (zh) | 2017-07-14 | 2017-07-14 | 一种二维地形正演和改正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107255837B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109870723A (zh) * | 2019-03-18 | 2019-06-11 | 云南航天工程物探检测股份有限公司 | 基于山区地形校正的大功率电测深方法及系统 |
CN115267919A (zh) * | 2022-09-27 | 2022-11-01 | 山东省鲁南地质工程勘察院(山东省地质矿产勘查开发局第二地质大队) | 一种基于分布式高密度电法的地球物理勘探系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005337746A (ja) * | 2004-05-24 | 2005-12-08 | National Institute For Rural Engineering | 電気探査方法 |
CN101614828A (zh) * | 2009-07-15 | 2009-12-30 | 上海大学 | 高密度电法堤坝探测空洞位置的自动确定方法 |
CN102778699A (zh) * | 2011-05-13 | 2012-11-14 | 中国石油天然气集团公司 | 一种电磁数据地形改正方法 |
CN103064124A (zh) * | 2013-01-31 | 2013-04-24 | 江苏大学 | 一种校正电磁勘探中地形影响的比值方法 |
CN104714254A (zh) * | 2015-02-02 | 2015-06-17 | 中国科学院地质与地球物理研究所 | 一种基于瞬变电磁法探测煤矿积水采空区的方法和装置 |
CN105277992A (zh) * | 2015-12-02 | 2016-01-27 | 山东科技大学 | 利用赤道偶极装置对采场底板水情动态监测方法 |
-
2017
- 2017-07-14 CN CN201710575229.1A patent/CN107255837B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005337746A (ja) * | 2004-05-24 | 2005-12-08 | National Institute For Rural Engineering | 電気探査方法 |
CN101614828A (zh) * | 2009-07-15 | 2009-12-30 | 上海大学 | 高密度电法堤坝探测空洞位置的自动确定方法 |
CN102778699A (zh) * | 2011-05-13 | 2012-11-14 | 中国石油天然气集团公司 | 一种电磁数据地形改正方法 |
CN103064124A (zh) * | 2013-01-31 | 2013-04-24 | 江苏大学 | 一种校正电磁勘探中地形影响的比值方法 |
CN104714254A (zh) * | 2015-02-02 | 2015-06-17 | 中国科学院地质与地球物理研究所 | 一种基于瞬变电磁法探测煤矿积水采空区的方法和装置 |
CN105277992A (zh) * | 2015-12-02 | 2016-01-27 | 山东科技大学 | 利用赤道偶极装置对采场底板水情动态监测方法 |
Non-Patent Citations (1)
Title |
---|
敖怀欢 等: "高密度电法的地形影响和地形改正探讨", 《贵州地质》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109870723A (zh) * | 2019-03-18 | 2019-06-11 | 云南航天工程物探检测股份有限公司 | 基于山区地形校正的大功率电测深方法及系统 |
CN115267919A (zh) * | 2022-09-27 | 2022-11-01 | 山东省鲁南地质工程勘察院(山东省地质矿产勘查开发局第二地质大队) | 一种基于分布式高密度电法的地球物理勘探系统 |
CN115267919B (zh) * | 2022-09-27 | 2022-12-30 | 山东省鲁南地质工程勘察院(山东省地质矿产勘查开发局第二地质大队) | 一种基于分布式高密度电法的地球物理勘探系统 |
Also Published As
Publication number | Publication date |
---|---|
CN107255837B (zh) | 2019-12-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110333543B (zh) | 基于反射系数分析的低阻体解释及成像方法与系统 | |
CN110007357B (zh) | 一种航空tem和航空mt联合反演方法 | |
CN103064124A (zh) | 一种校正电磁勘探中地形影响的比值方法 | |
CN108427741B (zh) | 一种基于大量高精度控制点的dem相对误差评价方法 | |
CN107255837A (zh) | 一种二维地形正演和改正方法 | |
CN103954842A (zh) | 一种确定大型接地系统接地阻抗值的方法 | |
CN108008456B (zh) | 一种圈定热液型铀矿深部三维重点铀成矿有利靶区的方法 | |
CN117092702A (zh) | 孔-隧激发极化探水结构的施工方法及反演探水方法 | |
CN104704396A (zh) | 处理重力梯度数据的方法 | |
CN113761457B (zh) | 基于测量的重力异常数据提取局部重力异常的方法 | |
CN114488314B (zh) | 一种基于陆地与水下直流电联合测量的地质反演方法 | |
CN106248740B (zh) | 一种土壤电阻率湿度校正方法 | |
CN113433595B (zh) | 基于自然电场隧道裂隙水的超前预报方法 | |
CN111597752B (zh) | 平衡孔间敏感性的跨孔电阻率ct深度学习反演方法及系统 | |
CN108562952A (zh) | 一种野外岩层产状数据采集装置及方法 | |
Kamah et al. | Enhancement of subsurface geologic structure model based on gravity, magnetotelluric, and well log data in Kamojang geothermal field | |
CN103926626A (zh) | 航磁数据高度矫正方法及系统 | |
Zhang et al. | Uncertainties in estimating regional methane emissions from rice paddies due to data scarcity in the modeling approach | |
CN113297526A (zh) | 一种基于Wenner四极和大地电磁数据的水平分层土壤结构联合反演方法 | |
CN113325482A (zh) | 一种时间域电磁数据反演成像方法 | |
CN107677237B (zh) | 一种三维观测系统方位角的计算方法 | |
Martinez et al. | Gravity gradiometry processing and interpretation at the Kauring Test Site | |
CN106324691A (zh) | 一种用于铀矿勘探的土壤氡浓度数据调平校正方法 | |
CN114152989B (zh) | 一种二维非对称四极电法的拟断面图坐标计算方法 | |
CN106383363B (zh) | 一种基于二维充电法的地质勘探方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for 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: 20191213 Termination date: 20200714 |