CN102128028A - 用于评估水泥与地层界面胶结质量的超声信号反演方法 - Google Patents
用于评估水泥与地层界面胶结质量的超声信号反演方法 Download PDFInfo
- Publication number
- CN102128028A CN102128028A CN2010100227252A CN201010022725A CN102128028A CN 102128028 A CN102128028 A CN 102128028A CN 2010100227252 A CN2010100227252 A CN 2010100227252A CN 201010022725 A CN201010022725 A CN 201010022725A CN 102128028 A CN102128028 A CN 102128028A
- Authority
- CN
- China
- Prior art keywords
- signal
- cement
- sleeve pipe
- echo
- interface
- 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
Images
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明涉及用于评估水泥与地层界面胶结质量的超声信号反演方法，根据声波在套管井结构中的传播模型，仿真垂直入射套管内壁的换能器回波信号，与实测信号相比较，调整模型参数，使二者均方根误差最小，分离回波信号中的套管内壁反射波、套管内共振回波、二界面反射回波，最后用归一化的二界面回波信号幅度作为二界面水泥胶结的质量参数。与现有技术相比，本发明能将二界面回波信号从相互叠加的反射信号中单独提取出来，并消除了泥浆衰减、换能器发射能量、仪器偏心、椭圆井眼等变化因素的影响。
Description
技术领域
本发明涉及一种套管井水泥胶结质量评价技术，尤其是涉及一种用于评估水泥与地层界面胶结质量的反演方法。
背景技术
石油、天然气井采用套管/水泥/地层结构来保证井眼的稳定，其中水泥一方面起到固定套管的作用，另一方面阻止不同地层中的流体混合，确保产油(气)区不被污染。通过下井仪器可对水泥的胶结质量进行检测和评估。最早采用的CBL、VDL声波方法在井周上没有方向性，所以不能确定胶结缺陷所在的方位，且不能检测出竖直方向的沟槽缺陷。近年来，国外公司采用宽频带换能器组或绕轴旋转的单个超声换能器，垂直于套管内壁激发和接收波形，对井周进行全覆盖检测，目前采用此检测方法的仪器(斯伦贝谢USI，哈利伯顿CAST-F)已经可以检测出套管/水泥界面(以下称一界面)胶结缺陷，然而由于套管的高声阻抗，窄脉冲声波在套管内发生多次反射，持续时间长，淹没了后续的水泥/地层界面(以下称二界面)，因此无法对二界面胶结质量进行评价。同济大学胡文祥等人针对这个问题，根据声半波透射理论，在原有的检测方式不变的基础上，选用多周期(周期等于套管厚度共振频率倒数)猝发信号激励换能器，产生的声波可有效地透过钢套管，获得增强的二界面回波信号，为评价二界面胶结质量奠定了基础。但是，在许多情况下，二界面回波脉冲仍然与套管内壁反射信号、套管内共振反射信号在时间上部分叠加，仍无法直接用来评价套管井二界面水泥胶结质量。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种消除了不利因素影响的用于评估水泥与地层界面胶结质量的超声信号反演方法。
本发明的目的可以通过以下技术方案来实现：
一种用于评估水泥与地层界面胶结质量的超声信号反演方法，其特征在于，该方法包括以下步骤：
(1)利用换能器发射超声脉冲对实际套管内壁进行检测，得到实测信号，同时选取换能器垂直套管内壁发射超声波的声传播模型，仿真出换能器接收的信号，获得仿真信号；
(2)对实测信号和仿真信号进行归一化处理和波形时间对齐处理，并从实测信号和仿真信号中截取一段时间窗口中的信号；
(3)对所截取的实测信号和仿真信号进行信号变换处理，并计算二者均方根误差，调整仿真模型中传播时间和水泥表面等效阻抗，反演参数使得上述误差达到最小；
(4)通过反演参数和选取的声传播模型计算套管内壁反射回波、套管内共振回波以及水泥/地层界面回波，最后用归一化的二界面回波信号评价界面的胶结质量。
所述的超声波的声传播模型包括一维声传输线模型或三维声波动模型。
所述的实测信号和仿真信号进行波形时间对齐处理时通过计算二者相关系数方法得到时间差，然后对实测信号进行延时。
所述的截取一段时间窗口中的信号包括套管内壁反射回波、套管内共振反射回波。
所述的信号变换处理方法包括Hilbert变换求包络或小波分解变换。
与现有技术相比，本发明能将二界面回波信号从各波成分相互叠加的反射信号中单独提取出来，评价套管井二界面水泥胶结质量，并消除了泥浆衰减、换能器发射能量、仪器偏心、椭圆井眼等变化因素的影响。
附图说明
图1为超声换能器垂直入射套管井检测方式的剖面示意图；
图2为图1所示检测方式的超声换能器输出电信号；
图3为反演时间t1和水泥表面阻抗Z2流程；
图4为图2中电信号通过反演结果分离之后的独立脉冲。
图中11为多周期猝发电信号、12为超声换能器、13为辐射声波、14为泥浆、15为套管、16为水泥、17为地层、18为反射声波、19为仿真电信号、W0(t)为套管内壁的反射回波、W1(t)为套管内的共振反射回波、W2(t)为二界面的回波。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
实施例
图1是超声换能器垂直入射套管井的剖面图，多周期猝发电信号11激励超声换能器12振动，向泥浆14中辐射声波13，在泥浆14、套管15、水泥16、地层17中传播，遇到界面后产生反射声波18，由超声换能器12转换为仿真电信号19。
为得到仿真电信号19，首先要确定声传播模型，这里所考虑的声传播模型，可以是近似的一维声传输线模型，或者更精确的三维声波动模型，但并不限于此。
采用一维声传输线模型，则“泥浆/套管/水泥”的系统响应函数为：
其中，f为频率；c0，d0为泥浆的纵波速度、换能器与套管内壁的距离；Z0为泥浆声阻抗；Z1’为等效的套管表面阻抗；
c1，d1，t1为套管的纵波速度、厚度，声波单次传播时间；Z1为套管声阻抗
Z2为水泥表面声阻抗；最终，电信号19的频域表达式为
W(f)＝E(f)T(f)H(f)R(f)＝E(f)T(f)R(f)H(f)＝ETR(f)H(f)，
其中，E(f)，激励电信号；T(f)，换能器发射响应函数；R(f)，换能器接收响应函数；ETR(f)，换能器自身响应函数；H(f)，泥浆/套管/水泥结构系统响应函数。
ETR(f)事先使用厚钢块(厚度＞150mm)作为反射体，用相同的激励电信号激发超声换能器，记录反射信号ETR(t)后进行傅立叶变换得到，换能器与厚钢块之间的距离等于实际测量中换能器与套管内壁距离。
对频谱进行反傅立叶变换得到时域波形：
应用以上所述方法，给定初始参数c0，d0，Z0，Z1，t1，Z2就可以仿真“泥浆/套管/水泥”中的声波传播反射信号，信号中包含了套管内壁反射回波、套管内共振反射回波，但不包含二界面反射回波。将仿真信号与测量信号进行比较，调整水泥表面阻抗Z2、时间t1，其他参数不变，均方根误差最小时，模型参数最接近实际值。这是一个反演的过程。
由于实际测量环境中c0，d0这些值不可能精确知道，反演之前预先进行处理，可以消除所带来的影响，首先对信号进行归一化处理，消除泥浆衰减所带来的影响。其次，将实测信号与仿真信号的时间调整一致，消除由于仪器偏心、椭圆井眼、以及泥浆声速不同所造成c0，d0变化的影响。采用相关的方法计算实测波形与仿真波形的时间差，由于仪器电路采样间隔限制，采用抛物线拟合方法更精确估计相关曲线的极值和信号时差，将测量信号调整至与仿真信号对齐。
进一步，确定二者比较的时间窗口，该时间窗口内包含套管内壁的反射回波、套管共振回波，但不包含二界面反射回波。如图2所示，窗口的起始时间Tstart为套管内壁反射波起始位置，而时间窗口的宽度为Tw1。
Tw1＝(N+3)*Tres
其中，N为猝发电信号的周期数；Tres套管厚度共振频率的倒数；Tw2时间窗口紧接着窗口Tw1，该窗口截至到信号结束，二界面回波的幅度采用该窗口内的信号计算。
另外在实际情况中，还可根据波形分离的效果，手动修改时间窗口的宽度。
对所选择的窗口内，可以设定不同的权重，使得误差对部分信号段较敏感，当权重均为1时，表示各部分相同，而如果想增大窗口内后半段信号的权重，可使用逐渐增大的权重函数(如线性函数)，对于其他的要求，可选相应的权重函数。
接下来，计算所选窗口内的仿真信号与测量信号的均方根差。可以直接计算信号的均方根差，也可对信号进行变换之后计算均方根差，例如Hilbert变换后的包络的均方根差等。
接着，初始化声传播模型中泥浆声阻抗Z0、声套管Z1，并设定反演的参数t1和Z2初始值，使用单纯形方法调整参数使得仿真波形Wsim(t)逐渐逼近实测波形Wexp(t)，其流程为图3所示。
最后由反演参数值和声传播模型计算各分波，如图4所示：
套管内壁反射回波：
套管内共振反射回波：
二界面反射回波：W2(t)＝Wexp(t)-W0(t)-W1(t)
在反演之前，已经对测量信号进行了归一化处理，直接采用二界面反射回波幅度W2(t)的Tw2窗口内的最大幅度或均方根幅度作为二界面水泥胶结质量的评价参数。
因此这里所述的本发明，完全能够实现上述所提到的目的，虽然由于公开的目的，给出了当前本发明的一个实例，为了完成所要求的结果，在过程的细节中存在大量可以变化的地方，本领域的技术人员对这些变化是清楚的，这些可能的变化也包括在本发明所附的权利要求中。
Claims (5)
1.一种用于评估水泥与地层界面胶结质量的超声信号反演方法，其特征在于，该方法包括以下步骤：
(1)利用换能器发射脉冲对实际套管内壁进行检测，得到实测信号，同时选取换能器垂直套管内壁发射超声波的声传播模型，仿真出换能器接收的信号，获得仿真信号；
(2)对实测信号和仿真信号进行归一化处理和波形时间对齐处理，并从实测信号和仿真信号中截取一段时间窗口中的信号；
(3)对所截取的实测信号和仿真信号进行信号变换处理，并计算二者均方根误差，调整仿真模型中传播时间和水泥表面等效阻抗，反演参数使得上述误差达到最小；
(4)通过反演参数和选取的声传播模型计算套管内壁反射回波、套管内共振回波以及水泥/地层界面回波，最后用归一化的二界面回波信号评价界面的胶结质量。
2.根据权利要求1所述的一种用于评估水泥与地层界面胶结质量的超声信号反演方法，其特征在于，所述的超声波的声传播模型包括一维声传输线模型或三维声波动模型。
3.根据权利要求1所述的一种用于评估水泥与地层界面胶结质量的超声信号反演方法，其特征在于，所述的实测信号和仿真信号进行波形时间对齐处理时通过计算二者相关系数方法得到时间差，然后对实测信号进行延时。
4.根据权利要求1所述的一种用于评估水泥与地层界面胶结质量的超声信号反演方法，其特征在于，所述的截取一段时间窗口中的信号包括套管内壁反射回波、套管内共振反射回波。
5.根据权利要求1所述的一种用于评估水泥与地层界面胶结质量的超声信号反演方法，其特征在于，所述的信号变换处理方法包括Hilbert变换求包络或小波分解变换。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201010022725 CN102128028B (zh) | 2010-01-12 | 2010-01-12 | 用于评估水泥与地层界面胶结质量的超声信号反演方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201010022725 CN102128028B (zh) | 2010-01-12 | 2010-01-12 | 用于评估水泥与地层界面胶结质量的超声信号反演方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102128028A true CN102128028A (zh) | 2011-07-20 |
| CN102128028B CN102128028B (zh) | 2013-06-26 |
Family
ID=44266277
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 201010022725 Expired - Fee Related CN102128028B (zh) | 2010-01-12 | 2010-01-12 | 用于评估水泥与地层界面胶结质量的超声信号反演方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102128028B (zh) |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103033154A (zh) * | 2012-12-21 | 2013-04-10 | 北京工业大学 | 一种用于提高超声检测时间分辨率的自回归谱分析方法 |
| CN104295285A (zh) * | 2014-08-14 | 2015-01-21 | 中国海洋石油总公司 | 水泥环第二界面胶结状况的判定方法及系统 |
| CN105545292A (zh) * | 2015-12-30 | 2016-05-04 | 山东大学 | 一种泥浆液连续波信号的处理方法 |
| CN106383173A (zh) * | 2016-09-13 | 2017-02-08 | 中国海洋石油总公司 | 一种水泥声阻抗计算方法和装置 |
| CN106761679A (zh) * | 2017-02-28 | 2017-05-31 | 中国地质大学（北京） | 一种可偏心测试固井第一界面胶结质量的评价装置及方法 |
| EP2778673A3 (en) * | 2013-03-13 | 2018-01-10 | Rolls-Royce plc | Ultrasonic inspection method for diffusion bonded articles |
| CN110376287A (zh) * | 2019-06-27 | 2019-10-25 | 湖南科技大学 | 一种地应力测量用空芯包体安装质量检测装置及检测方法 |
| CN111335888A (zh) * | 2020-03-28 | 2020-06-26 | 西安石油大学 | 一种确定地质层中井筒性质的方法 |
| CN114274366A (zh) * | 2022-01-06 | 2022-04-05 | 南京航空航天大学 | 一种基于数字孪生的超声辅助切削表面微织构创建方法 |
| US11326445B2 (en) | 2015-10-09 | 2022-05-10 | Darkvision Technologies Inc. | Devices and methods for imaging wells using phased array ultrasound |
| CN115078479A (zh) * | 2022-06-24 | 2022-09-20 | 西南石油大学 | 一种油基泥浆介电常数测量方法 |
| CN115822564A (zh) * | 2023-02-20 | 2023-03-21 | 中海油田服务股份有限公司 | 套管厚度测量方法、装置、电子设备及计算机存储介质 |
| US11619125B2 (en) | 2015-03-16 | 2023-04-04 | Darkvision Technologies Inc | Device and method to image flow in oil and gas wells using phased array doppler ultrasound |
| CN116660997A (zh) * | 2023-08-02 | 2023-08-29 | 中海油田服务股份有限公司 | 套管内外介质声阻抗反演方法、装置及电子设备 |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104747178B (zh) * | 2015-02-10 | 2018-02-13 | 中国海洋石油总公司 | 一种获取声阻抗的方法和装置 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1038489A (zh) * | 1988-06-09 | 1990-01-03 | 中国科学院声学研究所 | 声波全面评价水泥胶结固井质量的方法 |
| US5216638A (en) * | 1989-04-26 | 1993-06-01 | Schlumberger Technology Corporation | Method and apparatus for the acoustic investigation of a casing cemented in a borehole |
| CN1098473A (zh) * | 1993-08-06 | 1995-02-08 | 西安石油勘探仪器总厂 | 二次界面固井质量判定方法 |
| CN1098472A (zh) * | 1993-08-06 | 1995-02-08 | 西安石油勘探仪器总厂 | 声波水泥胶结固井质量测井方法 |
| US5907131A (en) * | 1997-08-27 | 1999-05-25 | Computalog U.S.A., Inc. | Method and system for cement bond evaluation high acoustic velocity formations |
| US6041861A (en) * | 1997-12-17 | 2000-03-28 | Halliburton Energy Services, Inc. | Method to determine self-calibrated circumferential cased bond impedance |
| CN101042046A (zh) * | 2006-11-30 | 2007-09-26 | 中国石油天然气集团公司 | 一种套管井方位声波测井方法 |
-
2010
- 2010-01-12 CN CN 201010022725 patent/CN102128028B/zh not_active Expired - Fee Related
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1038489A (zh) * | 1988-06-09 | 1990-01-03 | 中国科学院声学研究所 | 声波全面评价水泥胶结固井质量的方法 |
| US5216638A (en) * | 1989-04-26 | 1993-06-01 | Schlumberger Technology Corporation | Method and apparatus for the acoustic investigation of a casing cemented in a borehole |
| CN1098473A (zh) * | 1993-08-06 | 1995-02-08 | 西安石油勘探仪器总厂 | 二次界面固井质量判定方法 |
| CN1098472A (zh) * | 1993-08-06 | 1995-02-08 | 西安石油勘探仪器总厂 | 声波水泥胶结固井质量测井方法 |
| US5907131A (en) * | 1997-08-27 | 1999-05-25 | Computalog U.S.A., Inc. | Method and system for cement bond evaluation high acoustic velocity formations |
| US6041861A (en) * | 1997-12-17 | 2000-03-28 | Halliburton Energy Services, Inc. | Method to determine self-calibrated circumferential cased bond impedance |
| CN101042046A (zh) * | 2006-11-30 | 2007-09-26 | 中国石油天然气集团公司 | 一种套管井方位声波测井方法 |
Cited By (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103033154B (zh) * | 2012-12-21 | 2015-06-03 | 北京工业大学 | 一种用于提高超声检测时间分辨率的自回归谱分析方法 |
| CN103033154A (zh) * | 2012-12-21 | 2013-04-10 | 北京工业大学 | 一种用于提高超声检测时间分辨率的自回归谱分析方法 |
| EP2778673A3 (en) * | 2013-03-13 | 2018-01-10 | Rolls-Royce plc | Ultrasonic inspection method for diffusion bonded articles |
| CN104295285A (zh) * | 2014-08-14 | 2015-01-21 | 中国海洋石油总公司 | 水泥环第二界面胶结状况的判定方法及系统 |
| CN104295285B (zh) * | 2014-08-14 | 2017-03-08 | 中国海洋石油总公司 | 水泥环第二界面胶结状况的判定方法及系统 |
| US11619125B2 (en) | 2015-03-16 | 2023-04-04 | Darkvision Technologies Inc | Device and method to image flow in oil and gas wells using phased array doppler ultrasound |
| US11326445B2 (en) | 2015-10-09 | 2022-05-10 | Darkvision Technologies Inc. | Devices and methods for imaging wells using phased array ultrasound |
| CN105545292A (zh) * | 2015-12-30 | 2016-05-04 | 山东大学 | 一种泥浆液连续波信号的处理方法 |
| CN106383173B (zh) * | 2016-09-13 | 2019-12-06 | 中国海洋石油集团有限公司 | 一种水泥声阻抗计算方法和装置 |
| CN106383173A (zh) * | 2016-09-13 | 2017-02-08 | 中国海洋石油总公司 | 一种水泥声阻抗计算方法和装置 |
| CN106761679A (zh) * | 2017-02-28 | 2017-05-31 | 中国地质大学（北京） | 一种可偏心测试固井第一界面胶结质量的评价装置及方法 |
| CN106761679B (zh) * | 2017-02-28 | 2020-08-28 | 中国地质大学（北京） | 一种可偏心测试固井第一界面胶结质量的评价装置及方法 |
| CN110376287A (zh) * | 2019-06-27 | 2019-10-25 | 湖南科技大学 | 一种地应力测量用空芯包体安装质量检测装置及检测方法 |
| CN110376287B (zh) * | 2019-06-27 | 2024-01-23 | 湖南科技大学 | 一种地应力测量用空芯包体安装质量检测装置及检测方法 |
| CN111335888A (zh) * | 2020-03-28 | 2020-06-26 | 西安石油大学 | 一种确定地质层中井筒性质的方法 |
| CN114274366A (zh) * | 2022-01-06 | 2022-04-05 | 南京航空航天大学 | 一种基于数字孪生的超声辅助切削表面微织构创建方法 |
| CN115078479A (zh) * | 2022-06-24 | 2022-09-20 | 西南石油大学 | 一种油基泥浆介电常数测量方法 |
| CN115078479B (zh) * | 2022-06-24 | 2023-12-12 | 西南石油大学 | 一种油基泥浆介电常数测量方法 |
| CN115822564A (zh) * | 2023-02-20 | 2023-03-21 | 中海油田服务股份有限公司 | 套管厚度测量方法、装置、电子设备及计算机存储介质 |
| CN116660997A (zh) * | 2023-08-02 | 2023-08-29 | 中海油田服务股份有限公司 | 套管内外介质声阻抗反演方法、装置及电子设备 |
| CN116660997B (zh) * | 2023-08-02 | 2023-09-29 | 中海油田服务股份有限公司 | 套管内外介质声阻抗反演方法、装置及电子设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102128028B (zh) | 2013-06-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102128028B (zh) | 用于评估水泥与地层界面胶结质量的超声信号反演方法 | |
| US10012749B2 (en) | Fast model based inversion of acoustic impedance of annulus behind casing | |
| US9829597B2 (en) | Model based inversion of acoustic impedance of annulus behind casing | |
| US10191173B2 (en) | Systems and methods for evaluating annular material using beamforming from acoustic arrays | |
| US10042069B2 (en) | Systems, methods, and computer-readable media for determining shear-wave to compressional-wave velocity ratios in well casings | |
| MXPA06001468A (es) | Formacion de imagenes acusticas en modos multiples en pozos revestidos. | |
| WO2013138956A1 (zh) | 一种确定可控震源最佳低频扫描信号的方法 | |
| CN102128029B (zh) | 一种用于套管井二界面的超声检测成像方法 | |
| CN103244104A (zh) | 利用频散校正提取套管波衰减率的方法 | |
| CA2579953C (en) | Determining impedance of material behind a casing in a borehole | |
| Viggen et al. | Simulation and modeling of ultrasonic pitch-catch through-tubing logging | |
| CN105223271A (zh) | 一种有限空间条件下水声无源材料隔声量时反聚焦测量方法 | |
| EP3879311B1 (en) | Method for determining the integrity of a solid bonding between a casing and a wellbore | |
| CN109958432B (zh) | 利用超声回波测井评价固井ⅱ界面胶结质量方法和装置 | |
| CN111980676A (zh) | 一种阵列声波测井评价固井质量的方法及处理装置 | |
| CN115992691A (zh) | 基于超声兰姆波的固井质量检测方法及装置 | |
| CN106383173B (zh) | 一种水泥声阻抗计算方法和装置 | |
| CN101270659A (zh) | 对套管井外地层声阻抗不连续界面检测的系统及方法 | |
| US11905816B2 (en) | Casing thickness estimation by frequency correlation | |
| Tang et al. | An acoustic-wave theory for casing bond evaluation using viscoelastic slip boundary modeling | |
| CN104295285A (zh) | 水泥环第二界面胶结状况的判定方法及系统 | |
| Robb et al. | Measurement of the in situ compressional wave properties of marine sediments | |
| Mosland et al. | Using spectrum-of-spectrum (SoS) filtering to extract direct and multipath arrivals from a frequency domain simulation. Comparison with cepstrum and time-gating methods | |
| More et al. | Multipath Modeling under Tank Environment for Underwater Acoustic Communication | |
| Hill et al. | Reinterpretation of the reciprocity theorem for the calibration of acoustic emission transducers operating on a solid |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130626 Termination date: 20220112 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |