CN101545981B - 可控震源地震数据零相位子波最小相位化方法 - Google Patents
可控震源地震数据零相位子波最小相位化方法 Download PDFInfo
- Publication number
- CN101545981B CN101545981B CN2008101029622A CN200810102962A CN101545981B CN 101545981 B CN101545981 B CN 101545981B CN 2008101029622 A CN2008101029622 A CN 2008101029622A CN 200810102962 A CN200810102962 A CN 200810102962A CN 101545981 B CN101545981 B CN 101545981B
- Authority
- CN
- China
- Prior art keywords
- vibroseis
- phase
- wavelet
- zero
- factor
- 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.)
- Active
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及油田的勘探、开发、开采技术,是提供高分辨率的地震图形的可控震源地震数据零相位子波最小相位化方法。实现步骤为:对地震数据处理,转换为可控震源零相位子波地震数据,计算可控震源零相位子波和自相关函数:求解反滤波因子,计算纯相位因子,标定因子,相位转化因子:用相位转化因子把零相位子波转化为最小相位子波,用校正后的地震数据绘制相位校正后的地震数据剖面。本发明采用通常的地震数据处理系统可有效的把可控震源零相位子波地震数据转化为最小相位地震数据,为反褶积等地震数据处理提供了必要的输入地震数据。
Description
技术领域
本发明涉及油田的勘探、开发、开采技术,具体是为反映地下地层层位、油藏提供高分辨率的地震图形的一种可控震源地震数据零相位子波最小相位化方法。
背景技术
地震勘探的过程,就是在地面上的一系列点上,利用人工激发地震波,地震波向地下传播,当遇到波阻抗(地震波在地层介质中向地下传播的速度与介质密度的乘积)界面(即上下地层波阻抗不相等面)时,在波阻抗界面上地震波产生反射现象,地震波传播方向发生改变,地震波开始向上传播,在地面上的一系列接收点上安置着接收器,接收向上传播的地震波数据,完成野外勘探。当人工激发地震波采用可控震源激发时,可控震源激发仪器激发产生可控震源扫描信号,它是一个振幅确定的线性升频正弦扫描信号。地面接收点上接收器接收到的信号,是这个线性升频正弦扫描信号从可控震源激发仪器位置处向下传播,当遇到波阻抗界面时,在波阻抗界面上产生反射,传播方向发生改变,开始向上传播,到达地面的信号。地面接收点上接收器接收到的信号就是地震数据记录,它表示地下地层介质的反射系数序列与可控震源扫描信号的褶积。然而,实际接收到的地震数据还包含着激发点和接收点空间位置和排列位置的信息和各种噪声干扰等。地震数据处理就是对野外勘探过程中向上传播的地震数据记录进行处理,保留反映地下地层波阻抗界面的信息,而消除其它的信息,这种信息就是叠后地震数据。这种叠后地震数据仅反映地下地层的结构和构造。
对于可控震源激发的可控震源地震数据记录,第一步就是消除记录上可控震源扫描信号的影响。把可控震源地震数据记录与可控震源扫描信号进行互相关,就可以得到去除可控震源扫描信号影响的地震数据,它是可控震源零相位子波与地下地层介质的反射系数序列的褶积。而炸药震源激发的地震数据是炸药震源最小相位子波与地下地层介质的反射系数序列的褶积。为了把可控震源零相位子波地震数据与炸药震源最小相位子波地震数据进行对比,可以把可控震源零相位子波地震数据的零相位子波转换为最小相位子波,这样一来两种数据就可以对比,且一起处理。另一方面,为了提高地震数据的信噪比和分辨率,地震数据处理流程中都包含着地震反褶积处理。而反褶积处理理论的最根本前提和要求就是地震数据子波的相位是最小相位。这样一来,对于可控震源零相位子波地震数据,必须把可控震源零相位子波地震数据的零相位子波转换为最小相位子波,然后才可以做地震反褶积处理以提高地震数据的信噪比和分辨率。因此对于可控震源地震数据,可控震源地震数据最小相位化处理是必不可少的处理过程。
在地震数据处理中,包括反褶积在内的许多处理都要求地震数据中子波是最小相位,但在很多情况下难以满足。随着地震勘探的发展,地震数据采集过程中对环境污染和破坏要求越来越严格。因此与炸药震源激发相比,陆上越来越多地使用可控震源激发采集地震数据。可控震源激发采集的地震数据经过自相关后其子波是零相位,因此在反褶积之前首先需要把零相位震源子波的地震数据转化为最小相位震源子波的地震数据。通常地震数据处理系统并没有把可控震源零相位震源子波的地震数据转化为最小相位震源子波的功能。对于实际的可控震源零相位震源子波的地震数据,用户要么使用子波整型来实现这一功能处理,要么就忽视零相位震源子波这一现实直接进行反褶积处理。这样虽然可以得出反褶积处理结果,但反褶积处理的效果会受到很大影响。由于可控震源激发采集的地震数据越来越多,因此把零相位震源子波的地震数据转化为最小相位震源子波的问题就越来越突出。
发明内容
本发明目的在于提供一种满足反褶积等地震数据处理对地震数据必须是最小相位震源子波要求,得到反映地下地层层位、油藏提供高分辨率的地震图形的可控震源地震数据零相位子波最小相位化方法。
本发明采用如下技术方案实现,包括以下步骤:
(1)用通常的可控震源激发手段激发和采集地震数据,对地震数据置标签和定义观测系统以及对地震数据进行互相关的预处理;
(2)采用下式表达可控震源的扫描信号为线性升频正弦波信号:
s(t)=A(t)sin2πfit (1)
式中:
式中:
s(t) ----可控震源扫描信号
fl ----可控震源扫描信号的起始频率(HZ)
fu ----可控震源扫描信号的终止频率(HZ)
fi ----可控震源扫描信号的瞬时频率(HZ)
T ----可控震源扫描信号的扫描长度(MS)
T1 ----可控震源扫描信号起始段镶边斜坡长度(MS)
T2 ----可控震源扫描信号终止段镶边斜坡长度(MS)
A ----可控震源扫描信号的扫描振幅
(3)对接收到的地震数据记录用可控震源扫描信号进行自相关处理,把地震数据记录转换为可控震源零相位子波地震数据,采用下式计算可控震源零相位子波,
式中w(t)表示可控震源子波,是一个零相位子波;
(4)采用下式计算可控震源零相位子波的自相关函数:
它也是对应的最小相位子波的自相关函数rvv(t);
(5)由最小平方滤波方程求解反滤波因子,
(6)由可控震源零相位子波和反滤波因子计算纯相位因子,
式中“*”表示褶积;
(7)由纯相位因子计算相位因子标定因子:
(8)计算可控震源地震数据零相位子波最小相位化处理的相位转化因子:
(9)用相位转化因子把零相位子波转化为最小相位子波,可控震源子波和相位转化因子的互相关函数,即最小相位子波序列v(t),
(10)用相位转化因子把零相位子波的可控震源地震数据转化为最小相位子波可控震源地震数据;对可控震源子波地震数据x(t),其对应的最小相位子波地震数据y(t)为,
(11)采用通常的方法根据得出相位校正后的地震数据绘制相位校正后的地震数据剖面。
本发明采用通常的地震数据处理系统可有效的把可控震源零相位子波地震数据转化为最小相位地震数据,并应用于实际地震数据处理中。
本发明仅对可控震源地震数据零相位子波的相位进行处理,而不改变数据的振幅谱,并且处理后的数据子波一定是最小相位,为反褶积等地震数据处理提供了必要的输入地震数据。
本发明还可以实现把零相位子波转化为最小相位子波。既可以实现叠前地震数据的最小相位化,也可以实现叠后地震数据的最小相位化。
附图说明
图1本发明可控震源扫描信号图;
图2子波对比图,(a)可控震源零相位子波,(b)对应可控震源零相位子波的最小相位子波,(c)相位转换因子。
图3零相位子波和最小相位子波振幅谱的对比图,(a)是零相位子波振幅谱,(b)是最小相位子波振幅谱。
图4为合成数据对比图,(a)是反射系数序列,(b)是应用图2(a)的零相位子波所作的合成数据,(c)是对(b)进行最小相位化处理后的最小相位子波数据。
图5为实际炮集数据对比图,(a)是一个可控震源激发且经过互相关处理的零相位子波炮集数据,(b)是对(a)进行最小相位化处理后的最小相位子波数据。
图6与OMEGA系统处理叠加数据对比图,(a)是一个可控震源激发且经过互相关处理的零相位子波叠加数据,(b)是OMEGA系统对叠前数据进行最小相位化处理后的最小相位子波叠加数据,(c)是GEOEAST/GRISYS系统对叠前数据进行最小相位化处理后的最小相位子波叠加数据。
具体实施方式
可控震源地震记录就是由可控震源激发,地面检波器接收的地震数据,它是可控震源扫描信号和反射系数序列的褶积。可控震源扫描信号的自相关就是可控震源零相位地震子波。把可控震源地震记录与可控震源扫描信号进行自相关就得到可控震源零相位子波地震数据,它是可控震源零相位子波和反射系数序列的褶积,与炸药震源激发,地面检波器接收的地震数据意义相同,具有可比性。
在可控震源地震数据零相位子波最小相位化方法处理前,对可控震源地震记录已经进行了必要的自相关处理,地震数据已经是零相位子波的可控震源地震数据。本发明对可控震源地震数据进行最小相位化处理,根据可控震源扫描信号的参数,计算可控震源扫描信号、可控震源零相位地震子波以及对应的最小相位地震子波和它们之间的最小相位化转换因子,以实现可控震源地震数据零相位子波最小相位化方法处理。
本发明采用如下技术方案实现,包括以下步骤:
(1)用通常的可控震源激发手段激发和采集地震数据,对地震数据置标签和定义观测系统以及对地震数据进行互相关的预处理;
(2)采用下式表达可控震源的扫描信号为线性升频正弦波信号:
s(t)=A(t)sin2πfit (1)
式中:
式中:
s(t) ----可控震源扫描信号
fl ----可控震源扫描信号的起始频率(HZ)
fu ----可控震源扫描信号的终止频率(HZ)
fi ----可控震源扫描信号的瞬时频率(HZ)
T ----可控震源扫描信号的扫描长度(MS)
T1 ----可控震源扫描信号起始段镶边斜坡长度(MS)
T2 ----可控震源扫描信号终止段镶边斜坡长度(MS)
A ----可控震源扫描信号的扫描振幅。
(3)对接收到的地震数据记录用可控震源扫描信号进行自相关处理,把地震数据记录转换为可控震源零相位子波地震数据,采用下式计算可控震源零相位子波,
式中w(t)表示可控震源子波,是一个零相位子波;
(4)采用下式计算可控震源零相位子波的自相关函数:
它也是对应的最小相位子波的自相关函数rvv(t);
(5)由最小平方滤波方程求解反滤波因子,
(6)由可控震源零相位子波和反滤波因子计算纯相位因子,
式中“*”表示褶积;
(7)由纯相位因子计算相位因子标定因子:
(8)计算可控震源地震数据零相位子波最小相位化处理的相位转化因子:
(9)用相位转化因子把零相位子波转化为最小相位子波,可控震源子波和相位转化因子的互相关函数,即最小相位子波序列v(t),
(10)用相位转化因子把零相位子波的可控震源地震数据转化为最小相位子波可控震源地震数据;对可控震源子波地震数据x(t),其对应的最小相位子波地震数据y(t)为,
(11)采用通常的方法根据得出相位校正后的地震数据绘制相位校正后的地震数据剖面。
本发明实施情况如下:
首先模拟可控震源扫描信号的扫描长度为12000ms,可控震源扫描信号起始段和终止段镶边斜波长度均为500ms,可控震源扫描信号的起始和终止频率分别12Hz和50Hz,子波长度为300ms,可控震源扫描信号的扫描振幅为2000,地震数据采样间隔是2ms,生成的可控震源扫描信号如图1所示。对叠前地震数据记录用模拟的可控震源扫描信号进行互相关处理,得到可控震源零相位子波地震数据。对可控震源扫描信号进行自相关函数计算,得到了可控震源零相位子波,如图2(a)所示,利用零相位子波,通过求解矩阵方程,得出纯相位因子,如图2(c)所示,纯相位因子与零相位子波互相关,得出最小相位子波,如图2(b)所示。零相位子波和最小相位子波振幅谱的对比如图3所示,(a)是零相位子波振幅谱,(b)是最小相位子波振幅谱。图4为合成数据对比,(a)是反射系数序列,(b)是应用图2(a)的零相位子波所作的合成数据,(c)是对(b)进行最小相位化处理后的最小相位子波数据。图5为实际数据对比,(a)是一个可控震源激发且经过互相关处理的零相位子波炮集数据,(b)是对(a)进行最小相位化处理后的最小相位子波数据。图6与OMEGA系统处理叠加数据对比,(a)是一个可控震源激发且经过互相关处理的零相位子波叠加数据,(b)是OMEGA系统对叠前数据进行最小相位化处理后的最小相位子波叠加数据,(c)是GEOEAST/GRISYS系统对叠前数据进行最小相位化处理后的最小相位子波叠加数据。OMEGA系统对最小相位化处理算子进行滤波、平滑处理,本方法对转换算子并未做任何处理。从处理的叠加剖面效果上看,总体上本方法与OMEGA系统处理结果相当,但在局部细节上优于OMEGA系统。
Claims (1)
1.一种可控震源地震数据零相位子波最小相位化方法,其特征在于包括以下步骤:
1)用通常的可控震源激发手段激发和采集地震数据,对地震数据置标签和定义观测系统以及对地震数据进行互相关的预处理;
2)采用下式表达可控震源的扫描信号为线性升频正弦波信号:
s(t)=A(t)sin2πfit (1)
式中:
式中:
s(t) ----可控震源扫描信号
fi ----可控震源扫描信号的起始频率;
fu ----可控震源扫描信号的终止频率;
fi ----可控震源扫描信号的瞬时频率;
T ----可控震源扫描信号的扫描长度;
T1 ----可控震源扫描信号起始段镶边斜坡长度;
T2 ----可控震源扫描信号终止段镶边斜坡长度;
A ----可控震源扫描信号的扫描振幅;
3)对接收到的地震数据记录用可控震源扫描信号进行自相关处理,把地震数据记录转换为可控震源零相位子波地震数据,采用下式计算可控震源零相位子波,
式中w(t)表示可控震源子波,是一个零相位子波;
4)采用下式计算可控震源零相位子波的自相关函数:
它也是对应的最小相位子波的自相关函数rvv(t);
5)由最小平方滤波方程求解反滤波因子,
求解出反滤波因子为
6)由可控震源零相位子波和反滤波因子计算纯相位因子,
式中“*”表示褶积;
7)由纯相位因子计算相位因子标定因子:
8)计算可控震源地震数据零相位子波最小相位化处理的相位转化因子:
9)用相位转化因子把零相位子波转化为最小相位子波,可控震源子波和相位转化因子的互相关函数,即最小相位子波序列v(t),
式中表示互相关;
10)用相位转化因子把零相位子波的可控震源地震数据转化为最小相位子波可控震源地震数据;对可控震源子波地震数据x(t),其对应的最小相位子波地震数据y(t)为,
11)采用通常的方法根据得出相位校正后的地震数据绘制相位校正后的地震数据剖面。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008101029622A CN101545981B (zh) | 2008-03-28 | 2008-03-28 | 可控震源地震数据零相位子波最小相位化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008101029622A CN101545981B (zh) | 2008-03-28 | 2008-03-28 | 可控震源地震数据零相位子波最小相位化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101545981A CN101545981A (zh) | 2009-09-30 |
CN101545981B true CN101545981B (zh) | 2012-07-18 |
Family
ID=41193217
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008101029622A Active CN101545981B (zh) | 2008-03-28 | 2008-03-28 | 可控震源地震数据零相位子波最小相位化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101545981B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102269824B (zh) * | 2010-06-02 | 2014-02-05 | 中国石油天然气集团公司 | 一种地震数据子波相位转换处理的方法 |
CN102798887B (zh) * | 2011-05-27 | 2016-03-30 | 中国石油天然气集团公司 | 一种提高可控震源记录初至质量的滤波方法 |
CN102426386B (zh) * | 2011-10-27 | 2013-06-26 | 吉林大学 | 脉冲编码可控震源的多维匹配冲击方法 |
CN103217716B (zh) * | 2012-01-18 | 2015-07-22 | 中国石油天然气集团公司 | 一种地震数据快速相位校正处理的方法 |
CN103605154B (zh) * | 2013-11-06 | 2016-04-06 | 中国石油集团东方地球物理勘探有限责任公司 | 一种可控震源分频同时激发方法 |
CN105572723B (zh) * | 2014-10-14 | 2017-09-12 | 中石化石油工程地球物理有限公司胜利分公司 | 基于自相关子波的可控震源扫描信号的设计方法 |
CN104991273B (zh) * | 2015-07-03 | 2017-04-26 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 提取叠前校正道集地震子波的方法 |
CN106569275B (zh) * | 2015-10-10 | 2018-10-02 | 中国石油化工股份有限公司 | 子波零相位化处理方法和装置 |
CN105954802B (zh) * | 2016-06-20 | 2018-02-02 | 中国石油天然气股份有限公司 | 一种岩性数据体的转换方法及装置 |
CN106199713A (zh) * | 2016-07-07 | 2016-12-07 | 崔宏良 | 一种地震勘探数值模拟的分析方法 |
CN113156511B (zh) * | 2020-01-07 | 2024-06-25 | 中国石油天然气集团有限公司 | 井炮和可控震源混合采集的地震资料的匹配方法及装置 |
CN113109865A (zh) * | 2020-01-09 | 2021-07-13 | 中国石油天然气集团有限公司 | 混源激发地震波的处理方法及装置 |
CN111856580B (zh) * | 2020-07-27 | 2022-02-15 | 广州海洋地质调查局 | 一种超远偏移距obs数据初至波能量增强方法及处理终端 |
CN112379436A (zh) * | 2020-11-02 | 2021-02-19 | 中国石油天然气集团有限公司 | 混合震源激发地震资料静校正量计算方法及装置 |
CN114384580B (zh) * | 2021-12-31 | 2023-05-02 | 同济大学 | 一种基于可控震源的理想子波定制方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1237351C (zh) * | 2004-03-03 | 2006-01-18 | 吉林大学 | 相控可控震源系统 |
-
2008
- 2008-03-28 CN CN2008101029622A patent/CN101545981B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1237351C (zh) * | 2004-03-03 | 2006-01-18 | 吉林大学 | 相控可控震源系统 |
Non-Patent Citations (1)
Title |
---|
刘金俊等.子波零相位化、反褶积与地震记录分辨率的关系.《海洋地质动态》.2000,第16卷(第4期),5-8. * |
Also Published As
Publication number | Publication date |
---|---|
CN101545981A (zh) | 2009-09-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101545981B (zh) | 可控震源地震数据零相位子波最小相位化方法 | |
Li et al. | Wave-equation dispersion inversion | |
CN102112894B (zh) | 用地震表面波的波形评估土壤性质 | |
CN1065631C (zh) | 地震振动器特征反褶积 | |
EP2689273B1 (en) | System and method for seismic data modeling and migration | |
EP2356492B1 (en) | Method for separating independent simultaneous sources | |
US8953411B2 (en) | Apparatus and method for imaging a subsurface using frequency-domain elastic reverse-time migration | |
CN103487835B (zh) | 一种基于模型约束的多分辨率波阻抗反演方法 | |
EP2375268B1 (en) | Method for separating up and down propagating pressure and vertical velocity fields from pressure and three-axial motion sensors in towed streamers | |
CN105911587B (zh) | 一种利用单程波算子的双程波叠前深度偏移方法 | |
US20090292476A1 (en) | Method of seismic data interpolation by projection on convex sets | |
Tauzin et al. | Receiver functions from seismic interferometry: a practical guide | |
WO2007143355A2 (en) | Diplet-based seismic processing | |
CN102262243B (zh) | 一种滤波法可控震源地震数据谐波干扰压制方法 | |
CN104024890A (zh) | 通过加权叠加同时采集的波场来压制4d噪音 | |
EP2113792A1 (en) | Spectral shaping inversion and migration of seismic data | |
CN107884829A (zh) | 一种联合压制浅海obc地震资料多次波的方法 | |
CN102269824B (zh) | 一种地震数据子波相位转换处理的方法 | |
CN110967734B (zh) | 基于快速傅立叶变换的虚源重构方法及系统 | |
Novais et al. | GPR velocity determination by image-wave remigration | |
CN106257309A (zh) | 叠后地震数据体处理方法及装置 | |
CN117031539A (zh) | 一种自监督深度学习地震数据低频重建方法及系统 | |
EP3387466B1 (en) | Efficient internal multiple prediction methods | |
Etienne et al. | Three-dimensional acoustic full waveform inversion-algorithm and application to Valhall | |
WO2013116550A1 (en) | Geological medium exploration |
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 |