CN107748395A - 一种多频电阻率勘探方法 - Google Patents

一种多频电阻率勘探方法 Download PDF

Info

Publication number
CN107748395A
CN107748395A CN201710974753.6A CN201710974753A CN107748395A CN 107748395 A CN107748395 A CN 107748395A CN 201710974753 A CN201710974753 A CN 201710974753A CN 107748395 A CN107748395 A CN 107748395A
Authority
CN
China
Prior art keywords
frequency
current
mrow
msub
electrode
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
CN201710974753.6A
Other languages
English (en)
Other versions
CN107748395B (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.)
Central South University
Original Assignee
Central South University
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 Central South University filed Critical Central South University
Priority to CN201710974753.6A priority Critical patent/CN107748395B/zh
Publication of CN107748395A publication Critical patent/CN107748395A/zh
Application granted granted Critical
Publication of CN107748395B publication Critical patent/CN107748395B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Electromagnetism (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供一种多频电阻率勘探方法。该方法步骤包括：1)采用能发射n个相互独立、不同频率值的电流信号的电法发送机通过接地的n个供电电极和1个无穷远供电电极分别发送不同频率的电流信号；2)在待测目标区通过电法接收机采集接地的两个测量电极之间的n个不同频率的电场总信号，并对不同频率的电场总信号进行分离，获得n个频率的电位差数据；3)采集供电电极和测量电极位置的三维坐标；4)基于电法勘探的视电阻率算法开展不同频率的电位差数据与视电阻率数据之间的换算；6)对所有视电阻率数据开展反演处理，获取待测目标区的真实电阻率分布特征；7)基于真实电阻率分布特征，分析待测目标区的地质特征。该方法有助于明显提高电法勘探的野外工作效率和勘探效果、降低电法勘探成本。

Description

一种多频电阻率勘探方法
技术领域
本发明涉及一种电法勘探领域的电阻率勘探方法。
背景技术
常规视电阻率法是通过A、B两个供电电极向大地供电，用M、N两个测量电极接收因供电电极A、B供电产生的电场及地下不同电阻率介质对电场的影响而在MN电极间形成的电位差，勘探深度与AB间距离或AO间距离(O点为MN中点)、AB向大地的供电强度、接收仪器的灵敏度和精度等有关。
上述方法是采用两个供电电极供电，为了测量不同深度的电阻率，每次测量完测量电极M和N之间的电位差后，需要重新调整供电电极A和B的位置，即为了获取不同深度的电阻率数据，需要多次改变供电电极的位置，从而该方法工作效率比较低，勘探成本高。
为了提高电阻率勘探工作效率，《一种电法勘探的多电极多个供电单元系统及其勘探方法》(专利号：CN201310473315)专利提出一种勘探系统和方法，该专利的摘要内容为：“一种电法勘探的多电极多个供电单元系统及其勘探方法，其中系统中，供电电极A和供电电极B的数目均等于电流发送单元的个数；发射控制器包括相互之间连接的发射电流调制控制器和发射延时控制器；接收机与发射电流调制控制器之间连接，每个电流发送单元均与发射延时控制器之间连接；每个电流发送单元均与接收机之间通讯连接；接收机上设有分别用于连接测量电极M和测量电极N，并且用以测量测量电极M和测量电极N之间的电位差的端子；每个电流发送单元上均分别连接供电电极A和对应的供电电极B。本发明电流发送单元的供电电流不超过1A，供电电压不超过200V，各发射单元功率小，易实现，在隧道中使用和安装在全金属的隧道全断面掘进机上安全”。该专利采用多套发射单元通过不同供电电极发射电流信号，以图同时获得不同勘探深度的电位差数据，进而获得不同勘探深度的电阻率数据。然而该专利存在一些问题：1)当多个不同供电电极同时供电时，无法判断在测量电极上获取的电位差数据是哪个或哪两个供电电极的电流信号所贡献，因为测量电极上获取的电位差数据是所有供电电极所提供的电位响应，也就无法换算不同勘探深度的电阻率数据，从而获取的电位差数据无法获得反映勘探区不同深度的电阻率特征，也就无法获取勘探区的地质特征；2)该专利需要采用多套发送单元，从而增加设备的复杂程度、设备数量，不利于提高勘探工作效率，也增加勘探成本和工作难度。
发明内容：
基于现有技术的有关问题，本发明提供一种多频电阻率勘探方法。本方法的步骤包括：
a)采用能同时发射n个相互独立、不同频率的电流信号的电法发送机通过接地的n个供电电极(A1,A2,A3...An)和1个无穷远供电电极C分别发送不同频率的电流信号，且每个供电电极(A1,A2,A3...An)通过的电流信号频率不同，n个相互独立、不同频率的电流信号均通过无穷远供电电极C；其中n为大于1的正整数；
b)在待测目标区通过至少一台电法接收机采集接地的两个测量电极M、N之间的n个不同频率的电场信号，并对n个不同频率的电场信号进行分离，获取每个频率的电位差数据；
c)利用电法接收机完成整个待测目标区的所有测量点电场信号采集，并分离出所有测量点、所有频率的电位差数据；
d)采集所有供电电极和测量电极位置的三维坐标；并按照电法勘探的勘探装置系数公式计算与供电电极和测量电极位置有关的勘探装置系数Ki(i＝1,2,3...n)；其中 其中XAi,YAi,ZAi为供电电极Ai的三维大地坐标值，XN,YN,ZN为测量电极N的三维大地坐标值，XM,YM,ZM为测量电极M的三维大地坐标值，i＝1,2,3...n；
e)基于视电阻率算法开展所有测量点的不同频率的电位差数据与视电阻率数据之间的转换；其中i为某个频率电流信号的编号，取值为闭区间[1，n]的正整数；ρSi为第i个频率的电位差数据换算的视电阻率，Ki为第i个频率的电流信号通过的供电电极Ai与两个测量电极M、N的三维坐标计算的勘探装置系数，ΔUi为第i个频率的电流信号在两个测量电极M、N贡献的电位差，Ii第i个频率的电流强度；
f)对上述所有视电阻率数据开展反演处理，获取待测目标区的真实电阻率分布特征；
g)基于真实电阻率分布特征，分析待测目标区的地质特征。
上述方法中提到的n个相互独立、不同频率的电流信号优选地相邻频率的电流信号频差为即其中fi和fi-1分别表示第i个和第i-1个电流信号的频率，其中i为闭区间[1，n]的正整数。
采用上述方法能实现多个频率的电流信号通过不同供电电极同时发射，同一台电法接收机则同时接收所有频率的电位差信号，通过电法接收机的频率分离功能，从而获取到不同频率的电位差信号，实现同时获取不同勘探深度的视电阻率数据，从而实现明显提高勘探效率，降低勘探成本的目的。
附图说明：
图1为本发明的方法流程图；
图2为本发明的野外工作布置示意图；
图中(A1,A2,A3,A4,A5,A6)代表6个供电电极，C代表无穷远供电电极，O代表测点O，M代表测量电极M，N代表测量电极N，1代表能同时发射6个相互独立的、频率值不同的频率电流信号的电法发送机，2代表能接收电法发送机发出的电流信号的电法接收机，3代表连接无穷远供电电极C和电法发送机的导线，4代表连接电法发送机和供电电极的六芯导线，5代表连接测量电极M和电法接收机的导线，6代表连接测量电极N和电法接收机的导线。
具体实施方式：
为了使本发明所解决的技术问题、技术方案及有益效果更加清楚明白，以下结合附图1、附图2及具体实施方式，对本发明作进一步详细说明。
若获取某个测点O从浅往深的电阻率分布特征，需要在六个位置布置供电电极，则本发明所提出的方法可以按如下步骤进行实现：
a)在一条勘探线上按照常规三极电测深的方式布置六个供电电极(A1,A2,A3,A4,A5,A6)和一个无穷远供电电极C，以及连接无穷远供电电极C和能同时发射六个频率(0.1、0.2、0.4、赫兹)的电法发送机1的导线3，连接六个供电电极(A1,A2,A3,A4,A5,A6)和电法发送机1的六芯导线4，六芯导线4的每根芯与电法发送机的不同的频率电流信号输出口相连；
b)在O点两边布设测量电极M和N，以及能同时接收六个频率(0.1、0.2、0.4、赫兹)的电法接收机2，并布置分别连接两个测量电极M、N和电法接收机2的导线5和6；
c)采用能同时发射六个相互独立、不同频率(0.1、0.2、0.4、赫兹)的电流信号的电法发送机1通过接地的六个供电电极(A1,A2,A3,A4,A5,A6)和一个无穷远供电电极C发送不同频率(0.1、0.2、0.4、赫兹)的电流信号，且六个供电电极(A1,A2,A3,A4,A5,A6)通过的电流信号频率分别为：0.1、0.2、0.4、赫兹；六个相互独立、不同频率(0.1、0.2、0.4、赫兹)的电流信号均通过无穷远供电电极C，并记录每个频率的电流强度Ii(i＝1,2,3,4,5,6)；
d)通过电法接收机2采集接地的两个测量电极M、N之间的六个不同频率(0.1、0.2、0.4、赫兹)的电场总信号，并对六个不同频率(0.1、0.2、0.4、赫兹)的电场总信号进行分离，获取每个频率(0.1、0.2、0.4、赫兹)的电位差数据ΔUi(i＝1,2,3,4,5,6)；
e)采集所有供电电极(A1,A2,A3,A4,A5,A6)和测量电极MN位置的三维坐标，并采用电法勘探的勘探装置系数公式计算与供电电极和测量电极位置有关的所有勘探装置系数Ki(i＝1,2,3,4,5,6)；其中 其中XAi,YAi,ZAi为供电电极Ai的三维大地坐标值，XN,YN,ZN为测量电极N的三维大地坐标值，XM,YM,ZM为测量电极M的三维大地坐标值，i＝1,2,3...n；
f)基于视电阻率算法开展不同频率(0.1、0.2、0.4、赫兹)的电位差数据与视电阻率数据之间的换算，获得视电阻率数据；其中i为某个频率电流信号的编号，取值为闭区间[1,6]的正整数；ρSi为第i个频率的电位差数据换算的视电阻率，Ki为第i个频率的电流信号通过的供电电极Ai与两个测量电极M、N的三维坐标计算的勘探装置系数，ΔUi为第i个频率的电流信号在两个测量电极M、N贡献的电位差，Ii第i个频率的电流强度；
g)对上述所有视电阻率数据开展反演处理，获取待测目标区测量电极M、N之间的真实电阻率分布特征；
h)基于真实电阻率分布特征，分析待测目标区测量电极M、N之间的地质特征。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明内。

Claims (2)

1.一种多频电阻率勘探方法，其特征在于，该方法包括如下步骤：
a)在开展电阻率勘探剖面上，采用能同时发射n个相互独立、不同频率值的电流信号的电法发送机通过接地的n个供电电极(A1,A2,A3...An)和1个无穷远供电电极C分别发送不同频率的电流信号，且每个供电电极(A1,A2,A3...An)通过的电流信号频率相互不同，n个相互独立、不同频率的电流信号均通过无穷远供电电极C；并记录每个频率的电流信号强度Ii(i＝1,2,3...n)；其中n为大于1的正整数；
b)在待测目标区通过至少一台电法接收机采集接地的两个测量电极M、N之间的n个相互独立、不同频率值的电场总信号，并对n个相互独立、不同频率值的电场总信号进行分离，获取每个频率的电位差数据ΔUi(i＝1,2,3...n)；
c)采集所有供电电极(A1,A2,A3...An)和测量电极位置的三维大地坐标，并按照电法勘探的勘探装置系数公式计算与供电电极和测量电极位置有关的勘探装置系数Ki(i＝1,2,3...n)；其中
<mrow> <msub> <mi>A</mi> <mi>i</mi> </msub> <mi>N</mi> <mo>=</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>X</mi> <mrow> <mi>A</mi> <mi>i</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>X</mi> <mi>N</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>Y</mi> <mrow> <mi>A</mi> <mi>i</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>Y</mi> <mi>N</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>Z</mi> <mrow> <mi>A</mi> <mi>i</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>Z</mi> <mi>N</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>,</mo> </mrow>
其中XAi,YAi,ZAi为供电电极Ai的三维大地坐标值，XN,YN,ZN为测量电极N的三维大地坐标值，XM,YM,ZM为测量电极M的三维大地坐标值，i＝1,2,3...n；
d)基于视电阻率算法开展测量点的不同频率的电位差数据与视电阻率数据之间的换算，获得视电阻率数据；其中i为某个频率电流信号的编号，取值为闭区间[1，n]的正整数；ρSi为第i个频率的电位差数据换算的视电阻率，Ki为第i个频率的电流信号通过的供电电极Ai与两个测量电极M、N的三维坐标计算的勘探装置系数，ΔUi为第i个频率的电流信号在两个测量电极M、N贡献的电位差，Ii第i个频率的电流强度；
e)对上述所有视电阻率数据开展反演处理，获取待测目标区的真实电阻率分布特征；
f)基于真实电阻率分布特征，分析待测目标区的地质特征。
2.一种多频电阻率勘探方法，其特征在于：所述的n个相互独立、不同频率的电流信号优选地相邻频率的电流信号频差为即其中fi和fi-1分别表示第i个和第i-1个电流信号的频率，其中i为闭区间[2，n]的正整数。
CN201710974753.6A 2017-10-19 2017-10-19 一种多频电阻率勘探方法 Expired - Fee Related CN107748395B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710974753.6A CN107748395B (zh) 2017-10-19 2017-10-19 一种多频电阻率勘探方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710974753.6A CN107748395B (zh) 2017-10-19 2017-10-19 一种多频电阻率勘探方法

Publications (2)

Publication Number Publication Date
CN107748395A true CN107748395A (zh) 2018-03-02
CN107748395B CN107748395B (zh) 2019-07-02

Family

ID=61252684

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710974753.6A Expired - Fee Related CN107748395B (zh) 2017-10-19 2017-10-19 一种多频电阻率勘探方法

Country Status (1)

Country Link
CN (1) CN107748395B (zh)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108732630A (zh) * 2018-06-06 2018-11-02 中南大学 一种4个供电点同时供电的激电扫面方法
CN108845361A (zh) * 2018-06-06 2018-11-20 中南大学 一种多个供电点同时供电的激电扫面方法
CN108845362A (zh) * 2018-06-06 2018-11-20 中南大学 一种多个供电点同时供电的激电扫面方法
CN108873083A (zh) * 2018-05-06 2018-11-23 东华理工大学 一种人工场源频率域电磁视电阻率测量方法
CN108957560A (zh) * 2018-06-06 2018-12-07 中南大学 一种适用于城镇地下空间探测的电法勘探方法
CN108957561A (zh) * 2018-06-06 2018-12-07 中南大学 一种适用于城镇地下空间探测的电法扫面方法
CN109001816A (zh) * 2018-06-06 2018-12-14 中南大学 一种多对供电点同时供电的激电扫面方法
CN109597130A (zh) * 2018-12-21 2019-04-09 中国地质科学院地球物理地球化学勘查研究所 一种电法测量方法及系统
CN110907999A (zh) * 2019-12-26 2020-03-24 中南大学 一种基于被动源电场信号测量的主动源电法勘探方法
CN111551997A (zh) * 2020-03-12 2020-08-18 上海环联生态科技有限公司 一种隐伏断裂层的勘察系统及勘察方法
CN112433257A (zh) * 2020-12-01 2021-03-02 中南大学 一种基于频分电法的矢量电法扫面方法
CN113031078A (zh) * 2021-04-06 2021-06-25 中南大学 一种用于渗漏检测的频分电法校正方法
CN113608271A (zh) * 2021-08-02 2021-11-05 中国地质科学院 一种地质勘探数据处理方法及系统
CN116931102A (zh) * 2023-09-18 2023-10-24 山东大学 一种多水平错频激励充电探测方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101216566A (zh) * 2008-01-15 2008-07-09 湖南继善高科技有限公司 一种适合于准远区的电偶源频率域电测深方法及装置
CN105089656A (zh) * 2015-05-26 2015-11-25 广西大学 测井装置和测井方法
CN103487836B (zh) * 2013-10-11 2016-06-29 钟世航 一种电法勘探的多电极多个供电单元系统及其勘探方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101216566A (zh) * 2008-01-15 2008-07-09 湖南继善高科技有限公司 一种适合于准远区的电偶源频率域电测深方法及装置
CN103487836B (zh) * 2013-10-11 2016-06-29 钟世航 一种电法勘探的多电极多个供电单元系统及其勘探方法
CN105089656A (zh) * 2015-05-26 2015-11-25 广西大学 测井装置和测井方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
中国矿业学院物探教研室: "《中国煤田地球物理勘探》", 30 June 1981, 煤炭工业出版社 *

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108873083A (zh) * 2018-05-06 2018-11-23 东华理工大学 一种人工场源频率域电磁视电阻率测量方法
CN108732630A (zh) * 2018-06-06 2018-11-02 中南大学 一种4个供电点同时供电的激电扫面方法
CN108845361A (zh) * 2018-06-06 2018-11-20 中南大学 一种多个供电点同时供电的激电扫面方法
CN108845362A (zh) * 2018-06-06 2018-11-20 中南大学 一种多个供电点同时供电的激电扫面方法
CN108957560A (zh) * 2018-06-06 2018-12-07 中南大学 一种适用于城镇地下空间探测的电法勘探方法
CN108957561A (zh) * 2018-06-06 2018-12-07 中南大学 一种适用于城镇地下空间探测的电法扫面方法
CN109001816A (zh) * 2018-06-06 2018-12-14 中南大学 一种多对供电点同时供电的激电扫面方法
CN109597130A (zh) * 2018-12-21 2019-04-09 中国地质科学院地球物理地球化学勘查研究所 一种电法测量方法及系统
CN110907999A (zh) * 2019-12-26 2020-03-24 中南大学 一种基于被动源电场信号测量的主动源电法勘探方法
CN110907999B (zh) * 2019-12-26 2021-04-02 中南大学 一种基于被动源电场信号测量的主动源电法勘探方法
CN111551997A (zh) * 2020-03-12 2020-08-18 上海环联生态科技有限公司 一种隐伏断裂层的勘察系统及勘察方法
CN112433257A (zh) * 2020-12-01 2021-03-02 中南大学 一种基于频分电法的矢量电法扫面方法
CN113031078A (zh) * 2021-04-06 2021-06-25 中南大学 一种用于渗漏检测的频分电法校正方法
CN113608271A (zh) * 2021-08-02 2021-11-05 中国地质科学院 一种地质勘探数据处理方法及系统
CN113608271B (zh) * 2021-08-02 2024-01-30 中国地质科学院 一种地质勘探数据处理方法及系统
CN116931102A (zh) * 2023-09-18 2023-10-24 山东大学 一种多水平错频激励充电探测方法及系统
CN116931102B (zh) * 2023-09-18 2023-12-19 山东大学 一种多水平错频激励充电探测方法及系统

Also Published As

Publication number Publication date
CN107748395B (zh) 2019-07-02

Similar Documents

Publication Publication Date Title
CN107748395A (zh) 一种多频电阻率勘探方法
CN102183792B (zh) 人工源频率域电磁测深装置及测量方法
CN103064120B (zh) 煤矿井下磁电综合探测方法
CN104360399B (zh) 长导线接地源地下瞬变电磁探测方法及装置
CN102062876A (zh) 一种全区电偶源频率域电测深方法
CN103809204B (zh) 一种野外音频大地电磁的数据采集方法
CN105785475B (zh) 水力压裂震电联合探测系统、探测方法及野外工作方法
CN112415615B (zh) 基于分布式光纤传感的时频电磁压裂监测系统及监测方法
CN105301663A (zh) 时频电磁勘探数据空中采集装置及系统
CN110031901A (zh) 一种变深电法勘测系统和一种勘测方法
CN115097531B (zh) 全区观测交替覆盖积分差分混合激励全信息电磁勘探方法
CN103955000A (zh) 一种三维电测深方法
CN108957560A (zh) 一种适用于城镇地下空间探测的电法勘探方法
CN110488365A (zh) 一种多极化大地电磁测深方法
CN106646624A (zh) 非对称电法测深方法
CN203849425U (zh) 地裂缝三维电法探测装置
CN105891892A (zh) 全区测量电流源频率域磁场垂直分量的装置及方法
CN108732630B (zh) 一种4个供电点同时供电的激电扫面方法
CN105911594B (zh) Zh-8多功能电法工作站
CN108761540B (zh) 一种频率域天然电场三维勘探方法
CN213398937U (zh) 基于分布式光纤传感的时频电磁压裂监测系统
CN205139386U (zh) 时频电磁勘探数据空中采集装置及系统
CN217981876U (zh) 一种基于折线形电极结构的电磁测深勘探系统
CN108845361B (zh) 一种多个供电点同时供电的激电扫面方法
CN107024722A (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: 20190702