CN1332225C - 一种采用深度域滤波消除线性与非线性干扰波的方法 - Google Patents
一种采用深度域滤波消除线性与非线性干扰波的方法 Download PDFInfo
- Publication number
- CN1332225C CN1332225C CNB2004101026490A CN200410102649A CN1332225C CN 1332225 C CN1332225 C CN 1332225C CN B2004101026490 A CNB2004101026490 A CN B2004101026490A CN 200410102649 A CN200410102649 A CN 200410102649A CN 1332225 C CN1332225 C CN 1332225C
- Authority
- CN
- China
- Prior art keywords
- data
- wave
- continuation
- omega
- reference surface
- 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
- 238000001914 filtration Methods 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000013213 extrapolation Methods 0.000 claims abstract description 16
- 230000009466 transformation Effects 0.000 claims abstract description 6
- 230000010363 phase shift Effects 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 4
- 238000012423 maintenance Methods 0.000 claims description 3
- 230000005534 acoustic noise Effects 0.000 claims description 2
- 230000003068 static effect Effects 0.000 claims description 2
- 230000010355 oscillation Effects 0.000 abstract 1
- 208000027418 Wounds and injury Diseases 0.000 description 2
- 230000006378 damage Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 208000014674 injury Diseases 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
一种采用深度域滤波消除线性和非线性干扰波的方法,对于常规采集到的地震炮集数据,对采集到的炮集数据进行数据增益,对数据沿空间方向进行道均衡和沿时间方向进行振幅补偿,使数据适应波动方程处理;对炮集进行傅立叶正变换和向下波场延拓,从地表面延拓到参考面;对延拓到参考面上的炮集进行深度域滤波,对深度域滤波后的参考面的炮集进行向上波场延拓和傅立叶反变换,从参考面延拓到地表面;能有效识别复杂近地表产生的线性和非线性干扰波并采用深度域滤波消除线性和非线性干扰波。
Description
技术领域
本发明涉及石油地震勘探数据处理技术,具体是一种采用深度域滤波消除地震记录中线性与非线性干扰波的方法。
背景技术
在地震勘探领域,传统的信号加强或噪声压制技术一般是利用信号与噪声在频率(低通、高通或带通滤波)、视速度(FK滤波)、极化方向(极化滤波)和统计特性(随机噪声衰减)等方面的差异来分离或预测信号与噪声,然后再对信号或噪声进行相应的处理,这类方法称为数据驱动的去噪方法。
然而,随着勘探区域地表与近地表地震地质条件的日趋复杂,激发源附近各种地面障碍物及近地表岩性变化造成种类繁多的干扰波、沙漠的沙丘鸣震、海底产生的多次或微屈多次波,以及油田开发区各种地面与地下设施产生的干扰波,这些复杂近地表区产生的线性与非线性干扰波经常表现为在频率、视速度、传播方向、极化方向及统计特性方面很难与反射信号有效区别开来,因此利用数据中信号与噪声差异的传统数据驱动的去噪方法不能完全分离噪声与信号,达不到理想的去噪效果。
发明内容
本发明所要解决的技术问题是提供一种可有效识别复杂近地表产生的线性与非线性干扰波并采用深度域滤波消除线性与非线性干扰波的方法。
本发明是采用深度域滤波消除线性与非线性干扰波的方法,对于采集的地震炮集数据,其具体步骤包括,
(1)对炮集数据进行数据增益,对数据沿空间方向进行道均衡并沿时间方向进行振幅补偿,使数据适应波动方程处理;
(2)对炮集数据进行傅立叶正变换与向下波场延拓,从地表面延拓到参考面;
(3)对延拓到参考面上的炮集数据进行深度域滤波;
(4)对深度域滤波后的炮集数据进行向上波场延拓及傅立叶逆变换,从参考面延拓到地表面;
(5)进行数据逆增益,恢复反射信号的原始形状,得到消除地表与近地表干扰波的有效地震波并成图。
本发明当数据的空间采样间隔较大时,对炮集数据进行增益后,可以选择对数据进行空间插值以减少空间假频,从而消除对假频数据的不正确延拓。
本发明在进行向下与向上延拓过程中,选择信号保持的相位移算子,利用倏逝波数衰减噪声。
本发明向下延拓过程中,改变信号的相位,不改变信号的振幅并识别噪声。
本发明在向上延拓的过程中,恢复信号的相位。
本发明对于二维勘探所得到的地震记录,采用二维傅立叶正变换与逆变换。
本发明对于二维勘探所得到的地震记录,所采用的相位移算子在滤波前进行振幅补偿与野外静校正。滤波前对记录的处理不应破坏记录的地表一致性。
本发明对于三维勘探得到的地震记录,进行三维傅立叶正变换与向下波场延拓,相位移算子自动扩展到三维情况,进行深度域滤波后,进行三维傅立叶逆变换并向上波场延拓,即完成了三维深度域滤波。
本发明以二维地震数据为例,采用的具体原理为:
倏逝波数的公式为
振幅衰减因子为e-kz。
观测波场p(x,z,t)符合以下方程:
如近地表速度C(z)是只沿垂直方向变化的,近地表模型离散为厚度为Δz,速度为cn的N层均匀介质,对方程(1)的x和t坐标应用傅利叶变换,可以得出:
定义垂直波数为
,那么方程(2)的解就是:
式中的两个系数A,B由所应用的地表条件确定。
依据式(3)进行波场向下延拓,公式为:
依据式(3)进行波场向上延拓,公式为:
本发明在综合考虑计算效率、适应观测***能力及信号保持等方面因素后,采用基于波场延拓的模型驱动去噪方法,压制与倏逝波数相关联的噪声,消除近地表线性与非线性干扰,完全保持信号的振幅与相位信息,可以得到无干扰的地震数据。
附图说明
图1基于波场延拓的模型驱动滤波方法示意图,经过向下延拓与向上延拓之后,在振幅保持区域内,信号不变,而在相位保持区域,噪声振幅被压制。
图2为含有两层介质及一个散射体的近地表模型,其中散射体尺度为10m*10m、埋深为100m、速度为1500m/s,两层介质的速度cn分别为2500m/s和3500m/s。激发与接收面设置在深度15m处,参考面在深度195m处,波场向下延拓到参考面后,避免了散射体4的波场干扰。
其中1-炮点,2-参考面,3-有效信号波,4-散射体。
图3原始炮集(左)、沿深度轴延拓的深度域滤波结果(中)和差道集(右),从差道集剖面可以看出滤掉了干扰而保存了有效反射波。
其中5-有效波,6-干扰波,7-直达波。
图4是对二维地震勘探数据的一个叠前单炮地震记录滤波前后的对比图象。其中(a)是采用本发明滤波前的地震波记录剖面,所圈的区域近地表干扰很强;(b)是采用本发明方法处理后的炮集记录图象;(c)是采用本发明滤除的干扰波图象。
图5是对上述二维地震勘探数据的一个叠前单炮地震记录滤波前后所抽取的CMP道集对比图象。其中(a)是采用本发明滤波前的地震波记录剖面,有效波被干扰波掩盖,所圈的区域近地表干扰很强;(b)是采用本发明方法处理后的炮集记录图象,有效信号增强;(c)是采用本发明滤除的干扰波图象。
图6是本发明对同一数据滤波后的叠加剖面对比,(a)是滤波前的叠加结果,(b)是滤波后的叠加结果,(c)是噪声。
图7是三维地震数据深度域滤波对比图,(a)是滤波前原始炮集数据,(b)滤波后的炮集数据,(c)被滤除的噪声。
具体实施方式
实施例1
对于采集的地震炮集数据,选取一条二维测线的地震数据实施深度域滤波处理,数据的基本情况为:
数据道长为4000ms,采样间隔为4ms,
测线数据共包括351炮,每炮800道。
1)首先进行数据增益,对数据沿空间x方向进行道均衡并沿时间方向进行振幅补偿,使数据适应波动方程处理。
2)其次进行二维傅立叶正变换及向下波场延拓,对炮集按公式(4)和(5)从地表面延拓到参考面,所选的参考面深度为150m,分为三层,每层的厚度与速度分别为:80m、40m与30m,2000m/s,2400m/s,2800m/s。
式中,ω和kx分别表示圆频率及水平方向的圆波数;P(kx,z,ω)为ω-kx域内在深度z处的波场。P(kx,0,ω)表示在ω-kx域的地面炮集记录;cn表示在z与z+Δz之间的介质层速度。
在向下延拓的过程中,与倏逝波数相关联的噪声置为零,且被识别,有效信号延拓到参考面上。
3)然后进行深度域滤波,设从震源出发到达参考面的直达波旅行时间为DT(x,t),对延拓到参考面上的炮集,时间小于DT(x,t)的波场即为干扰波,将其置零。
4)进行向上波场延拓及二维傅立叶逆变换,对深度域滤波后的炮集按公式(6)和(7)从参考面延拓到地表面。
5)最后进行数据逆增益,恢复反射信号的原始形状,得到消除地表与近地表干扰波的有效地震波,由于沿时间方向的增益会改变信号的波形,逆增益的目的是恢复反射信号的原始形状,以保证在去噪过程中对反射信号没有任何伤害。
经过深度域滤波处理,与倏逝波数相关联的噪声被压制了,而反射有效信号得到了保留,如图4所示。
为了进一步说明滤波的效果,将滤波后的数据抽取为CMP道集,图5是滤波前后所抽取的CMP道集对比图象,从图中可以看出,在第一和第三个CMP道集中,原本不可见的反射波信息,在滤波以后清晰可见。第二个CMP道集,近炮检距处呈双曲线状的反射波信噪比有所提高,远炮检距处的反射从无到有。
图6给出了滤波前与滤波后的一段叠加剖面以及相应的差道集叠加剖面。叠加剖面的浅、中、深层反射波同相轴都变得更加清晰。而差道集叠加剖面中,没有反射波成像,进一步说明本发明信号保持的特点。
实施例2
对于采集的地震炮集数据,选取一炮三维测线的地震数据实施深度域滤波处理,数据的基本情况为:
数据道长为6000ms,采样间隔为2ms,
数据包括8条测线,每条线192道。
Dx=40m,Dy=240m
1)首先进行数据增益,对数据沿空间x、y方向进行道均衡并沿时间方向进行振幅补偿,使数据适应波动方程处理。
2)其次进行三维傅立叶正变换并向下波场延拓,对炮集按公式(5)和(6)从地表面延拓到参考面,所选的参考面深度为150m,分为三层,每层的厚度与速度分别为:80m、40m与20m,2000m/s,2400m/s,2700m/s。
式中,ω,kx和ky分别表示圆频率,水平方向与垂直方向的圆波数;P(kx,ky,z,ω)为ω-kx域内在深度z处的波场。P(kx,ky,0,ω)表示在ω-kx-ky域的地面炮集记录;cn表示在z与z+Δz之间的介质层速度。
3)然后进行深度域滤波,设从震源出发到达参考面的直达波的旅行时间为DT(x,t),对延拓到参考面上的炮集,把时间小于DT(x,t)的部分置零。
4)进行向上波场延拓及三维傅立叶逆变换,对深度域滤波后的炮集按公式(10)和(11)从参考面延拓到地表面。
5)最后进行数据逆增益,恢复反射信号的原始形状,得到消除地表与近地表干扰波的有效地震波,由于沿时间方向的增益会改变信号的波形,逆增益的目的是恢复反射信号的原始形状,从而保证在去噪过程中对反射信号没有任何伤害。
经过深度域滤波处理,与倏逝波数相关联的噪声被压制了,而有效反射信号得到了保留,如图7所示。
Claims (3)
1、一种采用深度域滤波消除线性与非线性干扰波的方法,其特征在于包括以下步骤:
(1)采集地震炮集数据,对采集的炮集数据进行数据增益,对数据沿空间方向进行道均衡并沿时间方向进行振幅补偿,使数据适应波动方程处理;
当数据的空间采样间隔较大时,对炮集数据进行增益后,选择对数据进行空间插值以减少空间假频,消除对假频数据的不正确延拓;
(2)对炮集数据进行傅立叶正变换并向下波场延拓,从地表面延拓到参考面,向下延拓过程中,改变信号的相位,不改变信号的振幅并识别噪声;
(3)对延拓到参考面上的炮集数据进行深度域滤波;
(4)对深度域滤波后的炮集数据进行向上波场延拓及傅立叶逆变换,从参考面延拓到地表面,在向上延拓的过程中,恢复信号的相位;
在进行向下与向上延拓过程中,选择信号保持的相位移算子,利用倏逝波数衰减噪声;
(5)进行数据逆增益,恢复反射信号的原始形状,得到消除地表与近地表干扰波的有效地震波并成图,
对于二维勘探所得到的地震记录,采用以下二维傅立叶正变换与逆变换:
二维傅立叶正变换及向下波场延拓,对炮集按以下公式从地表面延拓到参考面,参考面深度为150m,分为三层,每层的厚度与速度分别为80m、40m与30m,2000m/s,2400m/s,2800m/s;
P(kx,z+Δz,ω)=P(kx,z,ω)e-ikzΔz (4)
式中,ω和kx分别表示圆频率及水平方向的圆波数;P(kx,z,ω)为ω-kx域内在深度z处的波场。P(kx,0,ω)表示在ω-kx域的地面炮集记录;cn表示在z与z+Δz之间的介质层速度,
在向下延拓的过程中,与倏逝波数相关联的噪声置为零,且被识别,有效信号延拓到参考面上,
然后进行深度域滤波,从震源出发到达参考面的直达波旅行时间为DT(x,t),对延拓到参考面上的炮集,时间小于DT(x,t)的波场即为干扰波,将其置零,
进行向上波场延拓及二维傅立叶逆变换,对深度域滤波后的炮集按以下公式从参考面延拓到地表面:
P(kx,z-Δz,ω)=P(kx,z,ω)eikzΔz (6)
2、根据权利要求1所述的一种采用深度域滤波消除线性与非线性干扰波的方法,其特征在于对于二维勘探所得到的地震记录,所采用的相位移算子在滤波前进行振幅补偿和野外静校正。
3、根据权利要求1所述的一种采用深度域滤波消除线性与非线性干扰波的方法,其特征在于对于三维勘探得到的地震记录,对地震记录进行三维傅立叶正变换与逆变换,相位移算子自动扩展到三维情况。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2004101026490A CN1332225C (zh) | 2004-12-29 | 2004-12-29 | 一种采用深度域滤波消除线性与非线性干扰波的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2004101026490A CN1332225C (zh) | 2004-12-29 | 2004-12-29 | 一种采用深度域滤波消除线性与非线性干扰波的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1797041A CN1797041A (zh) | 2006-07-05 |
CN1332225C true CN1332225C (zh) | 2007-08-15 |
Family
ID=36818244
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2004101026490A Expired - Fee Related CN1332225C (zh) | 2004-12-29 | 2004-12-29 | 一种采用深度域滤波消除线性与非线性干扰波的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1332225C (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101592738B (zh) * | 2008-05-30 | 2011-04-20 | 中国石油集团东方地球物理勘探有限责任公司 | 二维起伏海底多次波的确定方法 |
CN101604028B (zh) * | 2008-06-11 | 2011-09-07 | 中国石油天然气集团公司 | 一种利用面炮合成记录确定海底高程的方法 |
CN102338890B (zh) * | 2010-10-22 | 2013-04-24 | 中国石油天然气股份有限公司 | 一种地球物理勘探中的圆形窗带通保幅滤波数据处理方法 |
CN102338884B (zh) * | 2010-10-22 | 2013-11-06 | 中国石油天然气股份有限公司 | 物探中的椭圆窗方向带通保幅滤波数据处理方法 |
CN103293553B (zh) * | 2013-04-17 | 2016-08-17 | 中国海洋石油总公司 | 一种复杂海底上下缆地震采集数据边界元延拓校正方法 |
CN104459770B (zh) * | 2013-09-24 | 2017-06-16 | 中国石油化工股份有限公司 | 一种高维地震数据规则化方法 |
CN104597501B (zh) * | 2013-11-01 | 2017-07-07 | 中国石油天然气集团公司 | 炮集记录处理方法及炮集记录处理装置 |
CN104181600B (zh) * | 2014-08-25 | 2017-06-06 | 中国石油天然气集团公司 | 地震数据线性噪声衰减方法及装置 |
CN104991268B (zh) * | 2015-07-03 | 2017-08-29 | 中国地质大学(北京) | 一种真振幅偏移成像方法 |
CN110895351B (zh) * | 2018-09-12 | 2022-03-08 | 中国石油化工股份有限公司 | 一种压制次生干扰的方法 |
CN113721297B (zh) * | 2020-05-26 | 2024-02-20 | 中国石油化工集团有限公司 | 利用速度特征压制沙丘鸣震的方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU817627A1 (ru) * | 1979-06-22 | 1981-03-30 | Уральское Отделение Всесоюзногогосударственного Ордена Ленинаи Ордена Октябрьской Революциипроектного Института "Тепло-Электропроект" | Способ сейсмической разведки |
US4611312A (en) * | 1983-02-09 | 1986-09-09 | Chevron Research Company | Method of seismic collection utilizing multicomponent receivers |
SU1513409A1 (ru) * | 1987-07-17 | 1989-10-07 | Комплексная Геофизическая Экспедиция Научно-Производственного Объединения "Нефтегеофизика" | Способ сейсморазведки отраженными волнами |
US5095466A (en) * | 1989-10-30 | 1992-03-10 | Total Compagnie Francaise Des Petroles | Marine seismic prospecting |
RU2172003C1 (ru) * | 2000-12-20 | 2001-08-10 | Лаврик Александр Сергеевич | Способ определения скоростных неоднородностей пород в верхней части геологического разреза |
US6418380B1 (en) * | 1997-10-10 | 2002-07-09 | Compagnie Generale De Geophysique | Method for seismic processing and in particular for three-dimensional seismic exploration using seismic data migration |
CN1472547A (zh) * | 2002-08-01 | 2004-02-04 | 张叔伦 | 一种利用合成震源记录相位编码的叠前深度偏移方法 |
WO2004090573A2 (en) * | 2003-04-07 | 2004-10-21 | Paradigm Geophysical (Luxembourg) S.A.R.L.. | Seimsic imaging by wave migration using a krylov space expansion of the square root exponent operator |
-
2004
- 2004-12-29 CN CNB2004101026490A patent/CN1332225C/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU817627A1 (ru) * | 1979-06-22 | 1981-03-30 | Уральское Отделение Всесоюзногогосударственного Ордена Ленинаи Ордена Октябрьской Революциипроектного Института "Тепло-Электропроект" | Способ сейсмической разведки |
US4611312A (en) * | 1983-02-09 | 1986-09-09 | Chevron Research Company | Method of seismic collection utilizing multicomponent receivers |
SU1513409A1 (ru) * | 1987-07-17 | 1989-10-07 | Комплексная Геофизическая Экспедиция Научно-Производственного Объединения "Нефтегеофизика" | Способ сейсморазведки отраженными волнами |
US5095466A (en) * | 1989-10-30 | 1992-03-10 | Total Compagnie Francaise Des Petroles | Marine seismic prospecting |
US6418380B1 (en) * | 1997-10-10 | 2002-07-09 | Compagnie Generale De Geophysique | Method for seismic processing and in particular for three-dimensional seismic exploration using seismic data migration |
RU2172003C1 (ru) * | 2000-12-20 | 2001-08-10 | Лаврик Александр Сергеевич | Способ определения скоростных неоднородностей пород в верхней части геологического разреза |
CN1472547A (zh) * | 2002-08-01 | 2004-02-04 | 张叔伦 | 一种利用合成震源记录相位编码的叠前深度偏移方法 |
WO2004090573A2 (en) * | 2003-04-07 | 2004-10-21 | Paradigm Geophysical (Luxembourg) S.A.R.L.. | Seimsic imaging by wave migration using a krylov space expansion of the square root exponent operator |
Also Published As
Publication number | Publication date |
---|---|
CN1797041A (zh) | 2006-07-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100349011C (zh) | 地震数据处理中压制与激发源无关的背景噪声的方法 | |
US9103935B2 (en) | Method of first arrival picking of seismic refraction data | |
CN107526101B (zh) | 一种获取地震反射波的采集和处理方法 | |
CN101652678B (zh) | 海洋地震采集方法和*** | |
US8116168B1 (en) | Hybrid one-way and full-way wave equation migration | |
CN102305941B (zh) | 由叠前时间偏移直接扫描确定地层叠加品质因子方法 | |
CN102176053B (zh) | 提升波动方程叠前深度偏移成像效果的方法 | |
Sun et al. | Recognizing surface scattering in ground-penetrating radar data | |
CN104570125A (zh) | 一种利用井数据提高成像速度模型精度的方法 | |
CN1332225C (zh) | 一种采用深度域滤波消除线性与非线性干扰波的方法 | |
EP2419761A1 (en) | Interferometric seismic data processing | |
CN106896409B (zh) | 一种基于波动方程边值反演的变深度缆鬼波压制方法 | |
CN102313900A (zh) | 三维地震采集观测***的激发位置确定方法 | |
CN111290017B (zh) | 一种震电波场联合提取瑞雷波频散特征的面波勘探方法 | |
Masoomzadeh et al. | Nonstretch moveout correction of long-offset multichannel seismic data for subbasalt imaging: Example from the North Atlantic | |
CN102590862A (zh) | 补偿吸收衰减的叠前时间偏移方法 | |
CN104620132A (zh) | 对于逆时迁移使用多样性拍摄堆积的信号增强(DeSSeRT) | |
CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
CN104330826A (zh) | 一种去除复杂地表条件下多种噪音的方法 | |
CN101609167B (zh) | 基于起伏地表的井间地震波动方程叠前深度偏移成像方法 | |
CN104977615B (zh) | 一种基于模型统计拾取的深水obc资料多次波压制方法 | |
CN104570114B (zh) | 一种基于波场分解的逆时偏移噪音压制方法 | |
CN108957551B (zh) | 基于重构地面力信号的可控震源谐波压制方法 | |
CN103135133A (zh) | 一种多分量地震资料的矢量降噪方法及设备 | |
Li et al. | Two robust imaging methodologies for challenging environments: Wave-equation dispersion inversion of surface waves and guided waves and supervirtual interferometry+ tomography for far-offset refractions |
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: 20070815 |