CN105223608A - 一种含煤强屏蔽缝洞型储层的地震预测与描述方法 - Google Patents
一种含煤强屏蔽缝洞型储层的地震预测与描述方法 Download PDFInfo
- Publication number
- CN105223608A CN105223608A CN201510478904.XA CN201510478904A CN105223608A CN 105223608 A CN105223608 A CN 105223608A CN 201510478904 A CN201510478904 A CN 201510478904A CN 105223608 A CN105223608 A CN 105223608A
- Authority
- CN
- China
- Prior art keywords
- reservoir
- force
- fractured
- analysis
- curve
- 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
- 238000000034 method Methods 0.000 title claims abstract description 50
- 239000003245 coal Substances 0.000 title claims abstract description 36
- 238000010521 absorption reaction Methods 0.000 claims abstract description 12
- 238000001228 spectrum Methods 0.000 claims abstract description 9
- 230000002159 abnormal effect Effects 0.000 claims description 5
- 238000000605 extraction Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 abstract description 4
- 238000004064 recycling Methods 0.000 abstract description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000005553 drilling Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 208000010392 Bone Fractures Diseases 0.000 description 1
- 206010017076 Fracture Diseases 0.000 description 1
- 208000013201 Stress fracture Diseases 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种含煤强屏蔽缝洞型储层的地震预测与描述方法，主要用于该类储层地震资料目标处理以及精细描述。该方法首先通过正演模拟分析煤层特征，为实际资料煤层定位提供基础；再利用广义S变换时变子波谱模拟，提高目的层分辨率；接着利用层位控制多道动态匹配追踪方法，进行强屏蔽剥离，突出弱信号；最后利用裂缝检测方法、吸收衰减分析以及波阻抗反演进行对比分析，获得储层预测结果。本发明可以较好的解决裂缝难以识别问题、缝洞型储层横向非均质性强以及储层低孔、低渗的问题，提高储层预测的精度。
Description
技术领域
本发明涉及一种含煤强屏蔽缝洞型储层的地震预测与描述方法，属于地震资料解释领域。
背景技术
受多套煤层或煤线屏蔽的影响，碳酸盐岩地层的地震反射普遍具有能量弱、多次波干扰严重等缺陷；地震资料主频低，分辨率低，这无疑增大了精细解释的难度；并且该类储层储集空间类型多样，孔隙和裂缝叠合发育，裂缝基本上为微裂缝，裂缝密度小，充填程度高，测井响应很弱，构成了复杂的配置关系。碳酸盐岩缝洞型储层一般为低孔、低渗致密储层，储层横向非均质性强，常规属性存在方法敏感性问题，反演方法在横向预测上也存在较大不确定性问题。因此，利用合理的储层预测方法对含煤强屏蔽缝洞型储层进行精细描述至关重要。
发明内容
本发明的目的在于提供一种含煤强屏蔽缝洞型储层的地震预测与描述方法，该方法主要用于地震资料目标处理以及储层精细刻画。
本发明所采用的技术解决方案是：
一种含煤强屏蔽缝洞型储层的地震预测与描述方法，包括以下步骤：
S1通过对地震资料进行频谱和波形特征分析得到资料基本信息；通过对井曲线中声波时差曲线的分析，获取煤层以及目标层段速度和密度值，根据速度以及密度值进行煤层正演，分析煤层对目标层段的影响；
S2采用广义S变换时变子波谱模拟提高分辨率技术进行目标处理，同时采用层位控制多道动态匹配追踪方法进行强屏蔽剥离，得到新的地震数据体；
S3对步骤S2目标处理后的地震数据体利用方差体对异常区域进行检测，进一步采用蚂蚁体算法识别裂缝，在前面分析的基础上，采用单频体检测低频伴影，利用吸收衰减属性分析目的层衰减，圈定储层范围，并利用井点信息进行印证；
S4通过拟声波曲线约束稀疏脉冲反演进行阻抗分析，并结合吸收衰减分析对储层进行精细刻画。
优选的，所述拟声波曲线约束稀疏脉冲反演包括以下步骤：测井曲线重构步骤，地震子波提取步骤，低频模型建立步骤，反演参数选取步骤和稀疏脉冲反演步骤。
本发明的有益技术效果是：
本发明针对煤层造成的强屏蔽进行正演分析，在有效识别煤层的基础上进行强屏蔽剥离，同时进行提高分辨率处理，能突出弱信号，有效提高薄层识别能力，提高主频，展宽频带范围；裂缝检测方法、吸收衰减方法以及拟声波曲线约束稀疏脉冲反演的联合使用可以较好的解决裂缝难以识别问题、缝洞型储层横向非均质性强以及储层低孔、低渗的问题，提高储层预测的精度。
附图说明
图1为本发明的处理流程图。
图2为本发明的实际数据剖面以及频谱分析图。图中：(a)为地震剖面；(b)为时间1000-2000ms多地震道频谱图；(c)为时间1400-1600ms多地震道频谱图；(d)为时间1400-1550ms单地震道频谱图。
图3为某井点处煤层替换前后正演分析。图中：(a)为井曲线含煤层正演分析图；(b)为井曲线1470ms以上煤层替换后正演分析图；(c)为井曲线1480ms以上煤层替换后正演分析图。
图4为含煤系地层正演模拟图。图中：(a)含煤层地层速度和密度统计表；(b)含煤层模型(对应表a)；(c)25Hz子波正演结果；(d)50Hz子波正演结果。
图5为广义S变换时变子波谱模拟提高分辨率处理前后剖面对比图。图中：(a)为原始地震剖面；(b)为提高分辨率后剖面。
图6为去强屏蔽前后连井剖面对比图。图中：(a)为去强屏蔽前原始地震剖面；(b)为去强屏蔽后剖面。
图7为去强屏蔽前后沿层(T9b层位下25-40ms)属性对比图。(a)为去强屏蔽前能量半时属性；(b)为去强屏蔽后能量半时属性；(c)为去强屏蔽前均方根振幅属性；(d)为去强屏蔽后均方根振幅属性。
图8为原始地震数据以及方差体三维立体图。图中：(a)为原始地震数据体；(b)为对原始地震数据处理后得到的方差体。
图9为蚂蚁体数据沿层切片图。(a)为蚂蚁体T9b层位下18-28ms切片；(b)为蚂蚁体T9b层位下23-33ms切片；(c)为蚂蚁体T9b层位下28-38ms切片。
图10为过D1-502、D1-501、D4、D93、D92、D3、D52井单频剖面(去强屏蔽后)图。(a)为48Hz单频剖面；(b)为28Hz单频剖面。
图11为本发明中面积差值法吸收衰减原理图。
图12为面积差值法吸收衰减方法处理后数据连井剖面以及沿层切片(T9b层位下45-50ms)图。图中：(a)为高频衰减异常剖面图；(b)为高频衰减异常沿层切片图。
图13为拟声波曲线约束稀疏脉冲反演中曲线重构、反演阻抗连井剖面以及沿层切片图。图中：(a)为利用GR曲线重构声波时差曲线；(b)为反演阻抗连井剖面；(c)为反演阻抗沿层切片(T9b层位下25-30ms)。
具体实施方式
本发明提供一种含煤强屏蔽缝洞型储层的地震预测与描述方法，主要用于该类储层地震资料目标处理以及精细描述。该方法首先通过正演模拟分析煤层特征，为实际资料煤层定位提供基础；再利用广义S变换时变子波谱模拟，提高目的层分辨率；接着利用层位控制多道动态匹配追踪方法，进行强屏蔽剥离，突出弱信号；最后利用裂缝检测方法、吸收衰减分析以及波阻抗反演进行对比分析，获得储层预测结果。
下面结合附图对本发明作进一步说明。
如图1所示，本发明含煤强屏蔽缝洞型储层的地震预测与描述方法，主要包括以下步骤：
S1通过对地震资料进行频谱和波形特征分析得到资料主频等基本信息；通过对测井资料井曲线中声波时差曲线的分析，获取煤层以及目标层段速度以及密度值，根据速度以及密度值进行煤层正演，分析煤层对目标层段的影响，为后续实际地震资料煤层识别以及煤层引起的强屏蔽剥离打下基础。
S2针对地震资料纵向分辨率低的问题，采用广义S变换时变子波谱模拟提高分辨率技术进行提高分辨率处理，同时针对煤层强屏蔽掩盖有效信号问题，采用层位控制多道动态匹配追踪方法进行强屏蔽剥离。经过基于GST提高分辨率以及匹配追踪强屏蔽剥离目标处理后，得到新的地震数据体，为下一步储层预测提供条件。
S3对步骤S2目标处理后的地震数据体进行缝洞型储层精细描述，先利用方差体对异常区域进行检测，进一步采用蚂蚁体算法识别裂缝。考虑到缝洞型含气储层频率衰减较快，在前面分析的基础上，采用单频体检测低频伴影现象，以及利用吸收衰减属性分析目的层衰减，圈定储层范围，并利用井点进行印证。
S4通过拟声波曲线约束稀疏脉冲反演进行阻抗分析，并结合吸收衰减分析对储层进行精细刻画，获得储层预测结果。
上述拟声波曲线约束稀疏脉冲反演包括以下步骤：测井曲线重构步骤，地震子波提取步骤，低频模型建立步骤，反演参数选取步骤和稀疏脉冲反演步骤。
本发明获取的最终储层预测结果可参见图12和图13。经对比研究可以发现，本发明能够有效圈定碳酸盐岩缝洞型储层发育范围，提高了地震储层预测的精度。
下面是本发明的具体应用实例：
将本发明应用于某含煤层碳酸盐岩储层工区，利用目标处理完的数据，进行裂缝检测以及缝洞型储层预测。最终得到如图12(b)所示沿层吸收衰减切片以及图13(c)所示沿层波阻抗切片，通过综合分析图12(b)和图13(c)，可以较好地划分储层范围，提高储层预测的精度。
上述方式中未述及的有关技术内容采取或借鉴已有技术即可实现。
需要说明的是，在本说明书的教导下本领域技术人员还可以做出这样或那样的容易变化方式，诸如等同方式，或明显变形方式。上述的变化方式均应在本发明的保护范围之内。
Claims (2)
1.一种含煤强屏蔽缝洞型储层的地震预测与描述方法，其特征在于包括以下步骤：
S1通过对地震资料进行频谱和波形特征分析得到资料基本信息；通过对井曲线中声波时差曲线的分析，获取煤层以及目标层段速度和密度值，根据速度以及密度值进行煤层正演，分析煤层对目标层段的影响；
S2采用广义S变换时变子波谱模拟提高分辨率技术进行目标处理，同时采用层位控制多道动态匹配追踪方法进行强屏蔽剥离，得到新的地震数据体；
S3对步骤S2目标处理后的地震数据体利用方差体对异常区域进行检测，进一步采用蚂蚁体算法识别裂缝，在前面分析的基础上，采用单频体检测低频伴影，利用吸收衰减属性分析目的层衰减，圈定储层范围，并利用井点信息进行印证；
S4通过拟声波曲线约束稀疏脉冲反演进行阻抗分析，并结合吸收衰减分析对储层进行精细刻画。
2.根据权利要求1所述的一种含煤强屏蔽缝洞型储层的地震预测与描述方法，其特征在于，所述拟声波曲线约束稀疏脉冲反演包括以下步骤：测井曲线重构步骤，地震子波提取步骤，低频模型建立步骤，反演参数选取步骤和稀疏脉冲反演步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510478904.XA CN105223608B (zh) | 2015-08-07 | 2015-08-07 | 一种含煤强屏蔽缝洞型储层的地震预测与描述方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510478904.XA CN105223608B (zh) | 2015-08-07 | 2015-08-07 | 一种含煤强屏蔽缝洞型储层的地震预测与描述方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105223608A true CN105223608A (zh) | 2016-01-06 |
| CN105223608B CN105223608B (zh) | 2017-08-04 |
Family
ID=54992663
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510478904.XA Expired - Fee Related CN105223608B (zh) | 2015-08-07 | 2015-08-07 | 一种含煤强屏蔽缝洞型储层的地震预测与描述方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105223608B (zh) |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106353803A (zh) * | 2016-10-17 | 2017-01-25 | 易宗富 | 基于谱模拟的地震资料宽频处理方法 |
| CN106526671A (zh) * | 2016-09-29 | 2017-03-22 | 中国石油天然气股份有限公司 | 一种碳酸盐岩储层地震波响应特征分析方法 |
| CN107340536A (zh) * | 2016-05-03 | 2017-11-10 | 中国石油化工股份有限公司 | 地震反演剥离煤层强反射的方法和装置 |
| CN107355219A (zh) * | 2016-05-09 | 2017-11-17 | 中国石油化工股份有限公司 | 裂缝性地层模型和其使用方法 |
| CN107632319A (zh) * | 2017-07-31 | 2018-01-26 | 成都理工大学 | 基于gst的碳酸盐岩非均质储层中溶洞识别标定方法 |
| CN107643539A (zh) * | 2016-07-21 | 2018-01-30 | 中国石油化工股份有限公司 | 一种基于煤层地震响应特征分析剥离强屏蔽层的方法 |
| CN107797145A (zh) * | 2016-08-31 | 2018-03-13 | 中国石油化工股份有限公司 | 消除煤系地层影响恢复下伏地层地震反射振幅的方法 |
| CN107831544A (zh) * | 2017-01-16 | 2018-03-23 | 中国石油化工股份有限公司 | 针对沙一段生物灰岩薄储层地震混叠强信号处理方法 |
| CN108318921A (zh) * | 2018-01-11 | 2018-07-24 | 中国石油大学(华东) | 一种基于横向约束的快速地震随机反演方法 |
| CN108756867A (zh) * | 2018-05-11 | 2018-11-06 | 中国地质调查局油气资源调查中心 | 基于声波测井曲线和电阻率测井曲线进行压裂选层的方法 |
| CN109143326A (zh) * | 2018-07-18 | 2019-01-04 | 西安恒泰艾普能源发展有限公司 | 识别煤层顶底界面的方法和装置 |
| CN110907990A (zh) * | 2018-09-18 | 2020-03-24 | 中国石油化工股份有限公司 | 一种叠后地震裂缝定量预测方法及系统 |
| CN111208562A (zh) * | 2020-01-17 | 2020-05-29 | 中国海洋石油集团有限公司 | 一种碳酸盐岩储层中火成岩侵入范围的识别方法 |
| CN112965109A (zh) * | 2021-02-05 | 2021-06-15 | 中国石油化工股份有限公司 | 一种强反射屏蔽下的储层的预测方法及装置 |
| CN113325468A (zh) * | 2020-02-28 | 2021-08-31 | 中国石油天然气集团有限公司 | 储层缝洞分布范围预测方法及装置 |
| CN114428361A (zh) * | 2020-09-21 | 2022-05-03 | 中国石油化工股份有限公司 | 基于低频模型优化的表层岩溶缝洞体预测方法及电子设备 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103675913A (zh) * | 2012-09-04 | 2014-03-26 | 中国石油天然气集团公司 | 一种裂缝型储层检测方法 |
| CN103728654A (zh) * | 2012-10-16 | 2014-04-16 | 中国石油化工股份有限公司 | 一种碳酸盐岩油气储层预测的方法 |
-
2015
- 2015-08-07 CN CN201510478904.XA patent/CN105223608B/zh not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103675913A (zh) * | 2012-09-04 | 2014-03-26 | 中国石油天然气集团公司 | 一种裂缝型储层检测方法 |
| CN103728654A (zh) * | 2012-10-16 | 2014-04-16 | 中国石油化工股份有限公司 | 一种碳酸盐岩油气储层预测的方法 |
Non-Patent Citations (3)
| Title |
|---|
| ZHANG JUN-HUA,ET AL: "Low-frequency data analysis and expansion", 《APPLIED GEOPHYSICS》 * |
| 石战战，等: "基于匹配追踪算法的碳酸盐岩储层低频伴影识别方法研究", 《岩性油气藏》 * |
| 黄捍东，等: "广义Ｓ变换地震高分辨率处理方法研究", 《石油地球物理勘探》 * |
Cited By (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107340536B (zh) * | 2016-05-03 | 2019-02-01 | 中国石油化工股份有限公司 | 地震反演剥离煤层强反射的方法和装置 |
| CN107340536A (zh) * | 2016-05-03 | 2017-11-10 | 中国石油化工股份有限公司 | 地震反演剥离煤层强反射的方法和装置 |
| CN107355219A (zh) * | 2016-05-09 | 2017-11-17 | 中国石油化工股份有限公司 | 裂缝性地层模型和其使用方法 |
| CN107643539A (zh) * | 2016-07-21 | 2018-01-30 | 中国石油化工股份有限公司 | 一种基于煤层地震响应特征分析剥离强屏蔽层的方法 |
| CN107797145A (zh) * | 2016-08-31 | 2018-03-13 | 中国石油化工股份有限公司 | 消除煤系地层影响恢复下伏地层地震反射振幅的方法 |
| CN106526671A (zh) * | 2016-09-29 | 2017-03-22 | 中国石油天然气股份有限公司 | 一种碳酸盐岩储层地震波响应特征分析方法 |
| CN106526671B (zh) * | 2016-09-29 | 2018-11-16 | 中国石油天然气股份有限公司 | 一种碳酸盐岩储层地震波响应特征分析方法 |
| CN106353803B (zh) * | 2016-10-17 | 2018-08-10 | 易宗富 | 基于谱模拟的地震资料宽频处理方法 |
| CN106353803A (zh) * | 2016-10-17 | 2017-01-25 | 易宗富 | 基于谱模拟的地震资料宽频处理方法 |
| CN107831544A (zh) * | 2017-01-16 | 2018-03-23 | 中国石油化工股份有限公司 | 针对沙一段生物灰岩薄储层地震混叠强信号处理方法 |
| CN107632319A (zh) * | 2017-07-31 | 2018-01-26 | 成都理工大学 | 基于gst的碳酸盐岩非均质储层中溶洞识别标定方法 |
| CN108318921A (zh) * | 2018-01-11 | 2018-07-24 | 中国石油大学(华东) | 一种基于横向约束的快速地震随机反演方法 |
| CN108756867A (zh) * | 2018-05-11 | 2018-11-06 | 中国地质调查局油气资源调查中心 | 基于声波测井曲线和电阻率测井曲线进行压裂选层的方法 |
| CN108756867B (zh) * | 2018-05-11 | 2021-11-19 | 中国地质调查局油气资源调查中心 | 基于声波测井曲线和电阻率测井曲线进行压裂选层的方法 |
| CN109143326A (zh) * | 2018-07-18 | 2019-01-04 | 西安恒泰艾普能源发展有限公司 | 识别煤层顶底界面的方法和装置 |
| CN110907990A (zh) * | 2018-09-18 | 2020-03-24 | 中国石油化工股份有限公司 | 一种叠后地震裂缝定量预测方法及系统 |
| CN110907990B (zh) * | 2018-09-18 | 2021-11-19 | 中国石油化工股份有限公司 | 一种叠后地震裂缝定量预测方法及系统 |
| CN111208562A (zh) * | 2020-01-17 | 2020-05-29 | 中国海洋石油集团有限公司 | 一种碳酸盐岩储层中火成岩侵入范围的识别方法 |
| CN113325468A (zh) * | 2020-02-28 | 2021-08-31 | 中国石油天然气集团有限公司 | 储层缝洞分布范围预测方法及装置 |
| CN114428361A (zh) * | 2020-09-21 | 2022-05-03 | 中国石油化工股份有限公司 | 基于低频模型优化的表层岩溶缝洞体预测方法及电子设备 |
| CN112965109A (zh) * | 2021-02-05 | 2021-06-15 | 中国石油化工股份有限公司 | 一种强反射屏蔽下的储层的预测方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105223608B (zh) | 2017-08-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105223608A (zh) | 一种含煤强屏蔽缝洞型储层的地震预测与描述方法 | |
| US8427902B2 (en) | Method for monitoring a subsoil zone, particularly during simulated fracturing operations | |
| CN104142519B (zh) | 一种泥岩裂缝油藏预测方法 | |
| CN105652323B (zh) | 一种储层预测方法 | |
| CN105353411A (zh) | 超深层煤上储层砂体刻画与描述方法 | |
| CN110529087B (zh) | 一种评价地层水力压裂效果的方法及装置 | |
| CN111158044B (zh) | 潜山裂缝储集体油藏预测方法 | |
| CN102914793B (zh) | 一种非均质储层的非等效体建模方法及装置 | |
| CA2640176A1 (en) | Analysis of time-series data using singularities | |
| CN104458774A (zh) | 一种利用生物内部纳米微粒寻找隐伏矿床的方法 | |
| CN103616390A (zh) | 一种胶结充填体顶板裂隙状态无损探测方法 | |
| Li et al. | Automatic fracture–vug identification and extraction from electric imaging logging data based on path morphology | |
| CN105301647A (zh) | 区分灰质泥岩和砂岩的方法 | |
| Zhao et al. | Fracture propagation induced by hydraulic fracturing using microseismic monitoring technology: Field test in CBM wells in Zhengzhuang region, Southern Qinshui Basin, China | |
| CN113970785B (zh) | 暗河裂缝发育预测方法、系统、存储介质以及电子设备 | |
| Hei et al. | Evaluation of the downhole hydraulic fracturing using single-well shear wave scattering imaging-modeling and field applications | |
| CN112711068B (zh) | 一种砂岩中油气有效储层预测方法及装置 | |
| CN105204066A (zh) | 一种基于谱分解的煤层火成岩侵入位置直接指示方法 | |
| Cheng et al. | A Method of Reconstructing Ground Penetrating Radar Bscan for Advanced Detection Based on High-Order Synchrosqueezing Transform | |
| CN114167515B (zh) | 一种岩性圈闭有效性识别的方法 | |
| CN114325820A (zh) | 一种煤层下伏岩层富水性的分析方法及其系统 | |
| Clifford et al. | Gas detection from absorption attributes and amplitude versus offset with artificial neural networks in Grand Bay Field | |
| Lu et al. | Exploring the Influence of Seismic Source and Improvement Methods on Tunnel Seismic Prediction | |
| Li et al. | Processing of random roadway source signals based on a cross-correlation algorithm in the deconvolution domain | |
| He et al. | The research on the recognition technique of rock AE and noise signal |
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 |
Granted publication date: 20170804 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |