CN100349009C - 一种起伏地表地震数据处理的叠前深度偏移方法 - Google Patents
一种起伏地表地震数据处理的叠前深度偏移方法 Download PDFInfo
- Publication number
- CN100349009C CN100349009C CNB2004101026448A CN200410102644A CN100349009C CN 100349009 C CN100349009 C CN 100349009C CN B2004101026448 A CNB2004101026448 A CN B2004101026448A CN 200410102644 A CN200410102644 A CN 200410102644A CN 100349009 C CN100349009 C CN 100349009C
- Authority
- CN
- China
- Prior art keywords
- wave field
- field
- imaging
- wave
- partiald
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 230000008569 process Effects 0.000 title claims description 17
- 238000003384 imaging method Methods 0.000 claims abstract description 50
- 238000013213 extrapolation Methods 0.000 claims abstract description 13
- 230000021615 conjugation Effects 0.000 claims abstract description 8
- 238000012545 processing Methods 0.000 claims abstract description 5
- 238000013508 migration Methods 0.000 claims description 38
- 230000005012 migration Effects 0.000 claims description 38
- 230000001419 dependent effect Effects 0.000 claims description 6
- 238000012423 maintenance Methods 0.000 claims description 5
- 230000000694 effects Effects 0.000 abstract description 3
- 238000005096 rolling process Methods 0.000 abstract 5
- 239000000523 sample Substances 0.000 abstract 2
- 238000013459 approach Methods 0.000 abstract 1
- 230000003068 static effect Effects 0.000 description 5
- 238000005452 bending Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 230000001788 irregular Effects 0.000 description 3
- 238000012360 testing method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
Images
Abstract
一种起伏地表地震数据处理的叠前深度偏移方法,根据采集区域起伏地表的最高点确定成像的基准面,将起伏地表投影到规则网格上,使起伏地表之上填充一个接近近地表的常数速度;利用单程波动方程的傅立叶有限差分延拓算子分别将炮点的震源波场和检波点的接收波场沿着深度轴进行波场延拓,且保持输入基准面和输出基准面之间的波场为零;将延拓后的炮点和检波点的波场记录应用共轭相关成像条件进行成像,并且将每一个频率互相关的结果叠加,通过常规的显示软件将成像结果处理为剖面图像,对于起伏地表条件下的复杂构造的成像具有明显的成像效果。
Description
技术领域
本发明涉及反射波地震数据处理过程中的叠前偏移成像技术范畴,是一种起伏地表地震数据处理的叠前深度偏移方法。
背景技术
在地表条件比较复杂地区进行地震勘探的过程中,数据的采集是在起伏的地表上进行的,获得数据所在的大地基准面是一个曲面,而不是一个平面。这种地震数据在进行处理时,为提高成像质量需要采用波动方程叠前偏移成像技术,这时会面临两个问题:第一,常规的波动方程叠前深度偏移是基于水平地表的前提条件,与实际起伏地表的条件不符;第二,如果用基于水平地表的波动方程叠前深度偏移方法进行成像,它将扭曲地下构造的实际形态。如果用基于水平地表的方法进行波动方程叠前偏移成像,获得的将是一个弯曲的界面,其弯曲的规律是起伏地表的镜像形态。如果地下界面本来就是弯曲变化的,那么得到就是更为复杂的脱离了实际构造形态的图像,会给后续的地下岩层构造解释成图和钻井造成误导。
对于起伏地表条件下的构造成像,以往的办法之一是先把数据通过静校正移动到一个水平的基准面上,在水平的基准面和起伏的地表之间采用一个替换速度,然后用基于水平地表的波动方程叠前深度偏移来实现构造成像。这种方法的问题在于如果将水平基准面放在起伏地表最高点之上,那么静校正的移量过大;如果放在起伏地表的中线,那么替换速度和实际地下介质之间的速度差别较大。同时静校正时移只是考虑了射线垂直方向的移动,没有考虑射线水平分量的移动,和波的实际传播路径存在误差。以往的另一个办法是应用克希霍夫积分法叠前深度偏移技术,优点是可以灵活地处理起伏地表的条件和不规则的观测点分布,但是由于该方法是基于高频射线近似的,有些波动现象难以实现,因此对于复杂构造的成像精度较低。
发明内容
本发明是要提供一种基于起伏地表的炮域波动方程叠前深度偏移方法,它不但能够适应起伏的地表条件,解决垂向时移静校正误差,而且能够有效地提高复杂构造的成像精度。
本发明的基于起伏地表的炮域波动方程叠前深度偏移技术,具体步骤包括:
(1)输入叠前炮集地震数据。
(2)根据起伏地表的最高点确定成像输出的基准面,将起伏地表投影到规则网格上,起伏地表之上填充一个与近地表的常数速度接近的速度值。
(3)利用单程波动方程的傅立叶有限差分延拓算子分别将炮点的震源波场和检波点的接收波场沿着深度轴进行波场延拓,且保持输入基准面和输出基准面之间的波场为零。
(4)将延拓后的炮点和检波点的波场记录应用共轭相关成像条件进行成像。
(5)通过常规的显示软件将成像结果处理为地层剖面图像。
本发明的基于起伏地表的炮域波动方程叠前深度偏移方法,在输入基准面和输出基准面之间保持一个数值为零的波场层来解决波动方程叠前深度偏移中的起伏地表问题,吸取了克希霍夫积分叠前深度偏移中灵活处理起伏地表的优点。
本发明的基于起伏地表的炮域波动方程叠前深度偏移方法,将叠前数据设定在炮域,它能适应不规则的观测系统分布,尤其是陆地观测系统和海底电缆观测系统,吸取了克希霍夫积分叠前深度偏移中灵活处理不规则观测系统的优点。
本发明的基于起伏地表的炮域波动方程叠前深度偏移方法,利用傅立叶有限差分延拓算子来延拓波场,以提高复杂构造成像精度的,改善了克希霍夫积分叠前深度偏移对于复杂构造成像精度低的缺点。
本发明的具体实现原理如下:
选定叠前深度偏移成像输出的基准面为水平,且位于起伏地表的最高点之上,对水平基准面之下的空间进行规则的网格剖分,于是起伏的地表形态被离散化后成为规则化网格的一部分。
对于任何一个在起伏地表上采集的叠前炮集记录R(x,y,z(x,y),t),采用傅立叶变换将它变换到频率域后为R(x,y,z(x,y),ω),用上行波方程对其进行向下延拓,延拓后的记录为U(x,y,z,ω),该过程可以表示为:U(x,y,z,ω)=Wz[R(x,y,z(x,y),ω)],其中x,y分别为空间直角坐标系中两个水平方向的坐标,z是深度方向的坐标,ω是圆频率,Wz是延拓算子。对于炮点的震源波场为S(x0,y0,z(x0,y0),t),将它变换到频率域后为S(x0,y0,z(x0,y0),ω),可以采用下行波方程对它进行向下传播,传播后的记录为D(x,y,z,ω),该过程可以写为:
其中Wz *是共轭延拓算子。延拓过程中下一个深度的波场为上一个深度波场的延拓结果与原有的波场在该位置的叠加,同时输出基准面和输入基准面之间通过设置一个零波场层参与正常的波场传播和延拓过程。
对延拓后的波场进行共轭相关成像,并且对所有的频率求和,就得到起伏地表条件下叠前深度偏移的结果:
然后通过常规的显示软件将成像结果处理为地层剖面图像。
上述延拓算子是根据单程波方程通过逼近得到的,采用稳定的傅立叶有限差分算子实现,该算子可以分为:W=W1+W2+W3,其中:
本发明采用适应于起伏地表的炮域波动方程叠前深度偏移方法,具有如下特点,主要表现为:
(1)可以正确处理偏移成像中的起伏地表条件。
(2)将起伏地表、炮域和傅立叶有限差分算子波动方程延拓有机结合。
(3)对于起伏地表条件下的复杂构造具有明显的成像效果。
附图说明
图1是一个只有两地层的二维起伏地表模型,地表的最大起伏高差为250米。
图2是合成地震记录的激发接收射线路径示意图,只激发了8炮。
图3是前4炮的合成地震记录,可以看出,尽管地下反射界面为水平,但是炮记录的反射同相轴形态发生了严重弯曲。
图4是经过垂向静态时移将炮记录由起伏地表移动到地形中线上的炮记录,可以看出炮记录变得比较接近双曲线形态。
图5是对时移后的记录用基于水平地表的方法进行叠前深度偏移成像的结果,可以看出水平反射界面同相轴不连续并且存在变形。
图6是用本技术进行叠前深度偏移成像后的结果,成像的反射界面平直并且连续。
图7是对一个起伏地表的复杂构造模型用积分法叠前深度偏移的结果,可以看出成像剖面上存在着显著的成像噪声,同时,诸多小的倾斜界面没有成像。
图8是对起伏地表复杂构造模型用本方法进行叠前深度偏移的结果,可以看出,噪声水平显著降低,成像背景清楚,诸多小的成像界面成像清晰,成像效果比积分法具有明显的改善。
图9是实际资料的积分法叠前深度偏移结果,构造形态不明显。
图10是该实际资料用本方法进行叠前深度偏移的结果,一个向斜构造明显呈出来。
具体实施方式
一种起伏地表地震数据处理的叠前深度偏移方法,具体为以下步骤:
1)采集读取叠前炮集地震数据;
2)根据采集区域起伏地表的最高点确定成像基准面,将起伏地表投影到规则网格上,使起伏地表之上填充一个与近地表的常数速度接近的速度值;
3)利用单程波动方程的傅立叶有限差分延拓算子分别将炮点的震源波场和检波点的接收波场沿着深度轴进行波场延拓,且保持输入基准面和输出基准面之间的波场为零;
具体是:将激发点震源子波设置在炮点相应的规则化网格位置上,并将它变换到频率域;将接收点接收的地震数据波场设置在相应规则化网格位置上并且变换到频率域;对地震数据进行傅立叶分析;用初始波和终止波确定波场延拓的频率范围f1和f2,f1为5赫兹,取f2为60赫兹,确定延拓算子W=W1+W2+W3,其中:
参数a,b的取值和偏移的最大倾角为:
对f1和f2之间每一个频率的震源波场和接收点波场,都分别采用前面叙述的稳定的傅立叶有限差分算子对震源波场向下延拓;
延拓过程中下一个深度上的波场为上一个深度延拓后的波场和原有的波场在该位置的上的叠加,使用的延拓步长和速度场的网格相同;
同时输出基准面和输入基准面之间通过一个零数据层参与正常的波场传播和延拓过程。
4)将延拓后的炮点和检波点的波场记录应用共轭相关成像条件进行成像,并且将每一个频率互相关的结果叠加,这就完成了成像过程,输出叠前深度偏移的结果即完成了整个叠前深度偏移。
图7和图8是一个起伏地表条件下复杂构造合成数据的测试例子。
1)输入叠前炮集地震数据。
2)根据起伏地表的最高点确定成像输出的基准面,将起伏地表投影到规则网格上,根据道距确定网格间距为20m,起伏地表之上填充一个接近近地表的常数速度4000m/s。
3)利用单程波动方程的傅立叶有限差分延拓算子分别将炮点的震源波场和检波点的接收波场沿着深度轴进行波场延拓,且保持输入基准面和输出基准面之间的波场为零。具体地说,将激发点震源子波设置在炮点相应的规则化网格位置上,并将它变换到频率域,将接收点接收的地震数据波场设置在相应规则化网格位置上并且变换到频率域,通过对地震数据进行傅立叶分析,设定波场延拓的频率范围f1和f2,一般将它设定为地震数据的有效频带,一般地,取f1为5赫兹,取f2为60赫兹,可满足绝大多数据成像的要求。参数a,b,取65°倾角即可满足实际要求,即a=0.4761,b=0.3767。对f1和f2之间每一个频率的震源波场和接收点波场,都分别采用前面叙述的稳定的傅立叶有限差分算子对震源波场向下延拓;延拓过程中下一个深度上的波场为上一个深度延拓后的波场和原有的波场在该位置的上的叠加,使用的延拓步长和速度场的网格相同;同时输出基准面和输入基准面之间通过一个零数据层参与正常的波场传播和延拓过程。
4)将延拓后的炮点和检波点的波场记录应用共轭相关成像条件进行成像,并且将每一个频率互相关的结果叠加,这就完成了成像过程,输出叠前深度偏移的结果即完成了整个叠前深度偏移。
5)通过常规的显示软件将成像结果处理为地层剖面图像如图8所示,图7是积分法偏移的结果用来比较。
图9和图10是一个起伏地表条件下复杂构造的实际数据的测试例子。
1)输入叠前炮集地震数据。
2)根据起伏地表的最高点确定成像输出的基准面,将起伏地表投影到规则网格上,根据道距确定网格间距为25m,起伏地表之上填充一个接近近地表的常数速度2500m/s。
3)利用单程波动方程的傅立叶有限差分延拓算子分别将炮点的震源波场和检波点的接收波场沿着深度轴进行波场延拓,且保持输入基准面和输出基准面之间的波场为零。具体地说,将激发点震源子波设置在炮点相应的规则化网格位置上,并将它变换到频率域,将接收点接收的地震数据波场设置在相应规则化网格位置上并且变换到频率域,通过对地震数据进行傅立叶分析,设定波场延拓的频率范围f1和f2,一般将它设定为地震数据的有效频带,一般地,取f1为5赫兹,取f2为50赫兹,可满足绝大多数据成像的要求。参数a,b,取65°倾角即可满足实际要求,即a=0.4761,b=0.3767。对f1和f2之间每一个频率的震源波场和接收点波场,都分别采用前面叙述的稳定的傅立叶有限差分算子对震源波场向下延拓;延拓过程中下一个深度上的波场为上一个深度延拓后的波场和原有的波场在该位置的上的叠加,使用的延拓步长和速度场的网格相同;同时输出基准面和输入基准面之间通过一个零数据层参与正常的波场传播和延拓过程。
4)将延拓后的炮点和检波点的波场记录应用共轭相关成像条件进行成像,并且将每一个频率互相关的结果叠加,这就完成了成像过程,输出叠前深度偏移的结果即完成了整个叠前深度偏移。
5)通过常规的显示软件将成像结果处理为地层剖面图像如图10所示,图9是积分法偏移的结果。
Claims (3)
1、一种起伏地表地震数据处理的叠前深度偏移方法,其特征在于采用以下步骤:
1)用常规的手段采集并读取叠前炮集地震数据;
2)根据采集区域起伏地表的最高点确定成像的基准面,将起伏地表投影到规则网格上,使起伏地表之上填充一个与近地表的常数速度接近的速度值;
3)利用单程波动方程的傅立叶有限差分延拓算子分别将炮点的震源波场和检波点的接收波场沿着深度轴进行波场延拓,且保持输入基准面和输出基准面之间的波场为零;
4)将延拓后的炮点和检波点的波场记录应用共轭相关成像条件进行成像,并且将每一个频率互相关的结果叠加,通过常规的显示软件将成像结果处理为剖面图像。
所述3)的波场延拓采用以下步骤:
将激发点震源子波设置在炮点相应的规则化网格位置上,并将它变换到频率域;
将接收到的地震波场数据设置在相应规则化网格位置上并且变换到频率域;
对地震数据进行傅立叶分析;
用初始波和终止波确定波场延拓的频率范围;
对频率范围内每一个频率的震源波场和接收点波场,分别采用稳定的傅立叶有限差分算子对震源波场向下延拓;
延拓过程中下一个深度的波场为上一个深度波场延拓后的结果与原有的波场在该位置的叠加,使用的延拓步长和速度场的网格相同;
同时在基准面之间通过一个零数据层参与正常的波场传播和延拓过程。
2、根据权利要求1所述的一种起伏地表地震数据处理的叠前深度偏移方法,其特征在于:初始频率为5赫兹,终止频率波为60赫兹。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2004101026448A CN100349009C (zh) | 2004-12-29 | 2004-12-29 | 一种起伏地表地震数据处理的叠前深度偏移方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2004101026448A CN100349009C (zh) | 2004-12-29 | 2004-12-29 | 一种起伏地表地震数据处理的叠前深度偏移方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1797038A CN1797038A (zh) | 2006-07-05 |
CN100349009C true CN100349009C (zh) | 2007-11-14 |
Family
ID=36818241
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2004101026448A Expired - Fee Related CN100349009C (zh) | 2004-12-29 | 2004-12-29 | 一种起伏地表地震数据处理的叠前深度偏移方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100349009C (zh) |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101021567B (zh) * | 2007-02-07 | 2010-05-19 | 徐兆涛 | 地震资料处理方法 |
CN101285894B (zh) * | 2008-05-30 | 2011-02-09 | 中国科学院地质与地球物理研究所 | 起伏地表下采集的地震资料的直接叠前时间偏移方法 |
CN101840000B (zh) * | 2009-03-20 | 2011-12-07 | 中国石油天然气集团公司 | 一种平面波叠前深度偏移方法 |
CN102033244B (zh) * | 2010-10-22 | 2012-09-26 | 中国石油化工股份有限公司 | 一种适于浅层高精度的曲地表叠加成像方法 |
CN102176053B (zh) * | 2011-01-27 | 2012-10-17 | 中国科学院地质与地球物理研究所 | 提升波动方程叠前深度偏移成像效果的方法 |
CN102183786A (zh) * | 2011-02-12 | 2011-09-14 | 中国石油大学(华东) | 双复杂条件下保真振幅高斯束叠前深度偏移方法 |
CN102778690B (zh) * | 2011-05-13 | 2015-10-07 | 中国石油化工股份有限公司 | 一种基于混合基离散傅氏变换的波动方程叠前偏移性能优化方法 |
CN103064110B (zh) * | 2011-10-18 | 2015-11-18 | 中国石油化工股份有限公司 | 一种波动方程叠前偏移中的分层延拓成像方法 |
CN102636809B (zh) * | 2012-03-27 | 2014-05-28 | 中国科学院地质与地球物理研究所 | 一种传播角度域共成像点道集的生成方法 |
CN103293553B (zh) * | 2013-04-17 | 2016-08-17 | 中国海洋石油总公司 | 一种复杂海底上下缆地震采集数据边界元延拓校正方法 |
CN105093318B (zh) * | 2014-05-14 | 2018-05-25 | 中国石油化工股份有限公司 | 一种自适应波动方程波场延拓静校正方法 |
CN104808242B (zh) * | 2014-10-20 | 2017-12-19 | 重庆丰银包装材料有限公司 | 一种地震资料成像处理方法 |
CN104749623B (zh) * | 2014-10-20 | 2017-12-08 | 重庆硕奥科技有限公司 | 一种地震资料成像处理方法 |
CN105842731B (zh) * | 2016-05-11 | 2018-02-23 | 吉林大学 | 基于波场延拓原理的起伏地表组合震源波场定向方法 |
CN111610559B (zh) * | 2019-02-22 | 2022-05-06 | 中国石油天然气股份有限公司 | 目的层叠前深度偏移成像方法及装置 |
CN113126153A (zh) * | 2019-12-30 | 2021-07-16 | 中国石油天然气集团有限公司 | 基于数据组合的叠前深度偏移方法及装置 |
CN111999767B (zh) * | 2020-07-21 | 2023-09-26 | 中国石油天然气集团有限公司 | 起伏地表的偏移成像方法及装置 |
CN112799126B (zh) * | 2020-12-30 | 2022-04-08 | 中国地质大学(北京) | 沿起伏地表的地震数据重构方法、装置、介质及设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1118441A (zh) * | 1994-09-02 | 1996-03-13 | 埃克森生产研究公司 | 处理具有多次反射噪声的地震数据的方法 |
CN1266195A (zh) * | 2000-04-14 | 2000-09-13 | 西安石油学院 | 一种提高测量地震检波器精度的方法 |
WO2003022141A1 (en) * | 2001-09-13 | 2003-03-20 | Imagyn Medical Technologies, Inc. | A signal processing method and device for signal-to-noise improvement |
-
2004
- 2004-12-29 CN CNB2004101026448A patent/CN100349009C/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1118441A (zh) * | 1994-09-02 | 1996-03-13 | 埃克森生产研究公司 | 处理具有多次反射噪声的地震数据的方法 |
CN1266195A (zh) * | 2000-04-14 | 2000-09-13 | 西安石油学院 | 一种提高测量地震检波器精度的方法 |
WO2003022141A1 (en) * | 2001-09-13 | 2003-03-20 | Imagyn Medical Technologies, Inc. | A signal processing method and device for signal-to-noise improvement |
Also Published As
Publication number | Publication date |
---|---|
CN1797038A (zh) | 2006-07-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100349009C (zh) | 一种起伏地表地震数据处理的叠前深度偏移方法 | |
US5671136A (en) | Process for seismic imaging measurement and evaluation of three-dimensional subterranean common-impedance objects | |
CN107526101B (zh) | 一种获取地震反射波的采集和处理方法 | |
CN104570125B (zh) | 一种利用井数据提高成像速度模型精度的方法 | |
CN102193109B (zh) | 起伏地表采集的三维地震资料的直接叠前时间偏移方法 | |
CN102879819B (zh) | 保持地震波场运动学特征的地震数据处理方法及装置 | |
CN102841379B (zh) | 一种基于共散射点道集的叠前时间偏移与速度分析方法 | |
CN101957455A (zh) | 三维保幅叠前时间偏移方法 | |
CN102313900A (zh) | 三维地震采集观测系统的激发位置确定方法 | |
CN105093292A (zh) | 一种地震成像的数据处理方法和装置 | |
CN101598811A (zh) | 一种二维垂直地震剖面数据计算炮点静校正的方法 | |
CN102305941A (zh) | 由叠前时间偏移直接扫描确定地层叠加品质因子方法 | |
CN106896409A (zh) | 一种基于波动方程边值反演的变深度缆鬼波压制方法 | |
CN107817526A (zh) | 叠前地震道集分段式振幅能量补偿方法及系统 | |
CN103576200A (zh) | 低信噪比区浅层波阻抗界面静校正方法 | |
Wilson et al. | Single‐chamber silicic magma system inferred from shear wave discontinuities of the crust and uppermost mantle, Coso geothermal area, California | |
Farrell et al. | Refraction statics | |
CN104570073A (zh) | 一种适用于复杂高陡构造的双反射地震波成像方法 | |
CN105137479A (zh) | 一种面元覆盖次数的计算方法及装置 | |
CN101609163B (zh) | 基于波动理论的多尺度地震资料联合成像方法 | |
CN104977615A (zh) | 一种基于模型统计拾取的深水obc资料多次波压制方法 | |
CN106125139A (zh) | 一种三维地震数据处理方法及系统 | |
CN101551464B (zh) | 一种复杂地表三维勘探初至波时间确定方法 | |
US20200333485A1 (en) | Randomizing sweeps in a marine survey | |
CN102798888B (zh) | 一种利用非零井源距数据计算纵横波速度比的方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20071114 |