CN1797028A - 一种消除起伏地表成像误差的方法 - Google Patents
一种消除起伏地表成像误差的方法 Download PDFInfo
- Publication number
- CN1797028A CN1797028A CN 200410102632 CN200410102632A CN1797028A CN 1797028 A CN1797028 A CN 1797028A CN 200410102632 CN200410102632 CN 200410102632 CN 200410102632 A CN200410102632 A CN 200410102632A CN 1797028 A CN1797028 A CN 1797028A
- Authority
- CN
- China
- Prior art keywords
- extrapolation
- wave field
- wave
- simple harmonic
- curved 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000003384 imaging method Methods 0.000 title claims abstract description 8
- 238000013213 extrapolation Methods 0.000 claims abstract description 57
- 230000008030 elimination Effects 0.000 claims description 10
- 238000003379 elimination reaction Methods 0.000 claims description 10
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 3
- 230000005012 migration Effects 0.000 claims description 2
- 238000013508 migration Methods 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract 1
- 230000010363 phase shift Effects 0.000 description 4
- 230000003068 static effect Effects 0.000 description 4
- 238000005553 drilling Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 241000220317 Rosa Species 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
一种涉及地球物理勘探地震数据的利用简谐波拟合技术一种消除起伏地表成像误差的方法,采集观测系统地震数据;对需要进行波场外推的地震数据抽成共炮点道集或共检波点道集整理,获得相应的已知波场;对起伏地表至外推目的曲面之间的区域进行外推网格的剖分;将原始波场数据从时间域变换到频率域;用简谐波拟合的方法计算出对应于剖分网格的外推算子;利用外推算子对地震数据进行波场的外推计算;将外推到目的曲面的波场数据从频率域变换到时间域并成像,在地震数据处理中误差小,干扰小,精度高的消除起伏地表成像误差。
Description
技术领域
本发明涉及地球物理勘探地震数据的处理方法,特别是涉及一种利用简谐波拟合技术消除起伏地表成像误差的方法。
背景技术
起伏地表是油气地震勘探所面临的一个难题,当地表剧烈起伏时,常规的静校正方法(静态时移)难以完全消除成像存在的误差,使速度分析、叠加成像都受到不利影响。同样,当(静态时移)校正量较大时,基准面校正叠后数据仍存在类似情况。解决该问题的方法通常是用波场外推来代替静态时移,使不同倾角的波都能得到正确归位,从而使其后的处理步骤可达到良好成像效果。波场外推是地震数据处理中的一项非常重要的技术,长期以来已形成了许多适于平面波场外推的技术方法,如差分法和相移法等。这些算法在进行非平面波场外推时存在两个主要缺陷:其一在进行非平面波场外推时不得不对波场进行非波动的近似处理,这种处理实际上是对原波场的一个扰动,有时会造成较大成像误差;其二在许多地区由于地表条件的影响造成测线上观测点分布的非规则性,而差分法和相移法要求样点必须是等间隔的,为此,有时会在测线内插大量零数据道,从而在波场外推过程中产生干扰。此外,差分法不能随步长变化而保持较高精度,相移法对速度变化的适应能力差。
发明内容
本发明的目的在于提供一种在地震数据处理中误差小、干扰小、精度高的消除起伏地表成像误差的方法。
本发明的消除起伏地表成像误差的方法,包括以下步骤:
(1)采集观测系统地震数据;
(2)对需要进行波场外推的地震数据抽成共炮点道集或共检波点道集整理,获得相应的已知波场;
(3)将原始波场数据从时间域变换到频率域;
(4)对起伏地表至外推目的曲面之间的区域进行外推网格剖分;
(5)用简谐波拟合的方法计算出对应于剖分网格的外推算子;
(6)利用外推算子对地震数据进行波场外推计算;
(7)将外推到目的曲面的波场从频率域变换到时间域并成像。
本发明所述的消除起伏地表成像误差的方法,所使用的偏移算法公式为:
Mp=0
其中,M为一线性数值矩阵形式的外推算子,p为任一波场,简谐波拟合的计算公式为:
Mei=0
式中ei为简谐单波,i=1,2,......。
本发明所述的消除起伏地表成像误差的方法,求解一局部范围内外推算子时使用简谐波来拟合实际波场,公式为:
ETMT=0
其中,MT为外推算子的转置矩阵,E为简谐单波集合,ET为波场的转置矩阵,ET=(e1,e2,……)T。
本发明所述的消除起伏地表成像误差的方法,实现方法为:设在一局部周围的A、B、C、D、E、F六点,对应p,ei和M分别为
p=(pA,pB,pC,pD,pE,pF)T
M=(MA,MB,MC,MD,ME,MF)
由已知波场曲面l0(......A,B,C,......)向未知波场曲面l1(......D,E,F,......)外推是通过求解一个复(频率—空间域)三角方程组
本发明所述的消除起伏地表成像误差的方法,波场外推从地表曲面逐步外推,直到目的曲面,实现非平面波场外推过程。
本发明采用简谐波拟合实现自由曲面的波场外推,具有明显的效果:
(1)本发明可以很自然地在自由曲面上实现外推,克服了现有平面波场外推算法存在的缺陷;
(2)本发明可以很自然地适应外推步长的变化而保持较高的成像精度;
(3)本发明可以很自然地适应速度的变化;
(4)本发明可以很自然地实现吸收边界;
(5)本发明可以很自然地在外推过程中适应采样间隔的不均匀。
附图说明
图1为自由曲面波场外推网格示意图,波场由已知波场曲面l0(......A,B,C,......)外推向未知波场曲面l1(......D,E,F......);
图2为一起伏地表的速度模型,此模型地下介质速度的变化范围从1500m/s至3000m/s之间;
图3为位于速度模型底部用于正演的初始波场,波场由相交于一点的三组深度角分别为-45°、0°、45°的同相轴所组成;
图4为图3所示模型底部初始波场经正演(用相位移法实现)在起伏地表处观测得到的波场;
图5为把图4所示在起伏地表处观测的波场用简谐波拟合法(自由曲面外推法)由模型起伏地表外推至模型底部的结果。
图6为一个实际资料的地表高程,波场外推目的面高程为1000m,从地表至波场外推目的面之间的外推速度为2500m/s;
图7为波场外推前的实际资料;
图8为波场外推后的结果。
具体实施方式
利用本方法进行外推的实施步骤为:
①首先对需要进行波场外推的地震数据,如图7所示的数据进行输入整理,利用快速傅里叶变换(FFT)将原始波场数据从时间域变换到频率域;
②如图6所示对起伏地表曲面至外推目的面之间的区域进行外推网格剖分,L1为起始的起伏地表曲面,L2为第一次外推曲面,L3为第二次的外推曲面,以此类推,Ln为第n-1次的外推曲面,直至到波场外推目的面;
③根据剖分后的外推网格,利用如下的公式:
ETMT=0
分别计算出从第n个外推曲面到第n+1个外推曲面所相对应的外推算子;
④在每一个频率,对于一个网格的上下六点,具体如图6所示外推曲面L1上的A、B、C三点和外推曲面L2上的D、E、F三点,利用如下的公式:
通过求解频率—空间域复三角方程组,就可以由已知的波场(PA、PB、PC)得到下一个外推面上的波场(PD、PE、PF),依次计算外推曲面L1和外推曲面L2之间的所有网格,就可以将波场从外推曲面L1外推至外推曲面L2,循环进行外推曲面之间的外推,就可以将波场从起始的起伏地表曲面外推至目的面,从而实现了频率域波场数据外推计算,循环以上过程直到完成所有的频率计算;
⑤如图8所示将外推到目的曲面的波场数据利用逆快速傅里叶变换(FFT)从频率域变换到时间域,得到外推波场。
Claims (5)
1.一种消除起伏地表成像误差的方法,其特征在于:
(1)采集观测系统地震数据;
(2)对需要进行波场外推的地震数据抽成共炮点道集或共检波点道集整理,获得相应的已知波场;
(3)对起伏地表至外推目的曲面之间的区域进行外推网格剖分;
(4)将原始波场数据从时间域变换到频率域;
(5)用简谐波拟合的方法计算出对应于剖分网格的外推算子;
(6)利用外推算子对地震数据进行波场外推计算;
(7)将外推到目的曲面的波场数据从频率域变换到时间域并成像。
2.如权利要求1所述的消除起伏地表成像误差的方法,其特征在于:所使用的偏移算法的公式为:
Mp=0
其中,M为一线性数值矩阵形式的外推算子,p为任一波场,简谐波拟合的计算公式为:
Mei=0
式中ei为简谐单波,i=1,2,......。
3.如权利要求1所述的消除起伏地表成像误差的方法,其特征在于:求解一局部范围内外推算子时使用简谐波拟合实际波场,公式为:
ETMT=0
其中,MT为外推算子的转置矩阵,E为简谐单波集合,ET为波场的转置矩阵。
4.如权利要求1所述的消除起伏地表成像误差的方法,其特征在于:实现波场外推的公式为:
5.如权利要求1所述的消除起伏地表成像误差的方法,其特征在于:波场外推从地表曲面逐步外推,直到目的曲面。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2004101026325A CN100349003C (zh) | 2004-12-29 | 2004-12-29 | 一种消除起伏地表成像误差的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2004101026325A CN100349003C (zh) | 2004-12-29 | 2004-12-29 | 一种消除起伏地表成像误差的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1797028A true CN1797028A (zh) | 2006-07-05 |
CN100349003C CN100349003C (zh) | 2007-11-14 |
Family
ID=36818231
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2004101026325A Expired - Fee Related CN100349003C (zh) | 2004-12-29 | 2004-12-29 | 一种消除起伏地表成像误差的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100349003C (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102338887A (zh) * | 2010-07-14 | 2012-02-01 | 中国石油天然气集团公司 | 不规则尺寸空变网格层析成像静校正方法 |
CN102520438A (zh) * | 2011-12-15 | 2012-06-27 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 分析干扰影响范围的方法 |
CN102859392A (zh) * | 2010-06-29 | 2013-01-02 | 雪佛龙美国公司 | 炮点道集数据的波束法和去波束法 |
CN104570069A (zh) * | 2013-10-12 | 2015-04-29 | 中国石油集团东方地球物理勘探有限责任公司 | 一种从起伏地表直接成像的地震成像方法和装置 |
CN112799126A (zh) * | 2020-12-30 | 2021-05-14 | 中国地质大学(北京) | 沿起伏地表的地震数据重构方法、装置、介质及设备 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB9717409D0 (en) * | 1997-08-15 | 1997-10-22 | Geco Prakla Uk Ltd | A method of processing seismic data |
CA2428743C (en) * | 2001-01-31 | 2009-05-12 | Jack Lees | System and method for analyzing and imaging an enhanced three-dimensional volume data set using one or more attributes |
GB2387225B (en) * | 2001-12-22 | 2005-06-15 | Westerngeco As | A method of seismic surveying and a seismic surveying arrangement |
-
2004
- 2004-12-29 CN CNB2004101026325A patent/CN100349003C/zh not_active Expired - Fee Related
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102859392A (zh) * | 2010-06-29 | 2013-01-02 | 雪佛龙美国公司 | 炮点道集数据的波束法和去波束法 |
CN102338887A (zh) * | 2010-07-14 | 2012-02-01 | 中国石油天然气集团公司 | 不规则尺寸空变网格层析成像静校正方法 |
CN102338887B (zh) * | 2010-07-14 | 2013-09-04 | 中国石油天然气集团公司 | 不规则尺寸空变网格层析成像静校正方法 |
CN102520438A (zh) * | 2011-12-15 | 2012-06-27 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 分析干扰影响范围的方法 |
CN104570069A (zh) * | 2013-10-12 | 2015-04-29 | 中国石油集团东方地球物理勘探有限责任公司 | 一种从起伏地表直接成像的地震成像方法和装置 |
CN104570069B (zh) * | 2013-10-12 | 2017-10-27 | 中国石油集团东方地球物理勘探有限责任公司 | 一种从起伏地表直接成像的地震成像方法和装置 |
CN112799126A (zh) * | 2020-12-30 | 2021-05-14 | 中国地质大学(北京) | 沿起伏地表的地震数据重构方法、装置、介质及设备 |
CN112799126B (zh) * | 2020-12-30 | 2022-04-08 | 中国地质大学(北京) | 沿起伏地表的地震数据重构方法、装置、介质及设备 |
Also Published As
Publication number | Publication date |
---|---|
CN100349003C (zh) | 2007-11-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101614826B (zh) | 三维地震数据处理中实现面元均化的方法和装置 | |
CN111158049B (zh) | 一种基于散射积分法的地震逆时偏移成像方法 | |
CN108845317B (zh) | 一种基于分层介质格林函数的频域逆时偏移方法 | |
CN109633752B (zh) | 基于三维快速Radon变换的海上拖缆资料自适应鬼波压制方法 | |
CN113625337A (zh) | 一种极浅水高精度地震资料快速成像方法 | |
CN112099082B (zh) | 一种共面元共方位角道集的地震回折波走时反演方法 | |
CN111665556B (zh) | 地层声波传播速度模型构建方法 | |
CN1797028A (zh) | 一种消除起伏地表成像误差的方法 | |
CN106338766A (zh) | 基于分步傅里叶算法的叠前时间偏移方法 | |
CN104898166B (zh) | 一种处理地震数据的方法及装置 | |
CN109557588B (zh) | 一种煤矿井下二维矿震波速反演降维方法 | |
CN102998702B (zh) | 保幅平面波叠前深度偏移方法 | |
CN108303741A (zh) | 基于水检积分的地震双检合并鸣震压制方法 | |
CN110888158B (zh) | 一种基于rtm约束的全波形反演方法 | |
Mao et al. | Fast image decomposition in dip angle domain and its application for illumination compensation | |
CN110161565A (zh) | 一种地震数据重建方法 | |
CN111538084B (zh) | Ovt域数据转换成方位角度域成像道集的方法及系统 | |
Mao et al. | Target oriented 3D acquisition aperture correction in local wavenumber domain | |
CN111929731B (zh) | 地表一致性和非一致性联合静校正方法 | |
CN112327361B (zh) | 基于线性同相轴迭代追踪衰减的倾斜干扰剔除方法 | |
Yang et al. | 3D surface-related multiples elimination based on improved apex-shifted Radon transform | |
CN111983685B (zh) | τ-p域地表非一致性静校正方法 | |
CN115542395A (zh) | 一种基于正负倾角偏移体乘积的绕射波成像方法 | |
CN115903042A (zh) | 一种基于构造整形正则化的波形反演方法及其设备 | |
CN116125535A (zh) | 三维vsp成像的方法及装置 |
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: 20071114 |
|
CF01 | Termination of patent right due to non-payment of annual fee |