CN104570119A - 一种三维垂直地震剖面反射波拉伸校正方法 - Google Patents
一种三维垂直地震剖面反射波拉伸校正方法 Download PDFInfo
- Publication number
- CN104570119A CN104570119A CN201310520617.1A CN201310520617A CN104570119A CN 104570119 A CN104570119 A CN 104570119A CN 201310520617 A CN201310520617 A CN 201310520617A CN 104570119 A CN104570119 A CN 104570119A
- Authority
- CN
- China
- Prior art keywords
- offset
- road
- zero
- reflection wave
- vsp
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种三维垂直地震剖面反射波拉伸校正方法,根据零偏移距(或近偏移距)与远偏移距垂直地震剖面反射波速度确定零偏移距道(近偏移距)与远偏移距道反射点的对应关系,选取一定的时窗,采用最小平方法,进行最佳逼近,求取反算子,并将反算子应用到相应的远偏移距垂直地震剖面道上,消除拉伸畸变对反射波成像的影响。本发明可以使得浅层以及远偏移距拉伸较大的同相轴得到校正,有效改善成像的分辨率,提高三维垂直地震剖面反射成像的质量。
Description
技术领域
本发明涉及地球勘探技术,属于地震资料处理过程中反射波校正的技术领域,具体涉及一种三维垂直地震剖面反射波拉伸校正方法。
背景技术
垂直地震剖面(VSP)是一种井中地震观测技术。与地面地震相比,垂直地震剖面资料的信噪比高,分辨率高,波的运动学和动力学特征明显。垂直地震剖面技术提供了地下地层结构同地面测量参数之间最直接的对应关系,可以为地面地震资料处理解释提供精确的时深转换及速度模型,为零相位子波分析提供支持。
国内外石油工业界从20世纪80年代就开始了对VSP技术的大量研究开发工作,在井旁精细构造成像、井旁断层识别、井旁地层岩性描述、地震波衰减、速度各向异性以及孔隙压力预测和孔隙度估算等方面取得了许多实际的应用效果。进入90年代,VSP数据采集方式明显增多,为适应不同的地质目的,设计出了各种各样的观测系统,从零偏VSP到非零偏VSP,再到多方位VSP、多偏移距VSP和变偏VSP,直至三维垂直地震剖面(3DVSP)。这些新的观测方式,为储层预测和描述(精细构造特征、岩性、孔隙度、各向异性、地层弹性参数等研究)提供了有利的支持。
对于三维垂直地震剖面(3DVSP)观测资料,由于地震子波在不同偏移距、不同深度存在较大差异,导致地下任何点的反射在利用共检波点道集成像过程中在时间方向都存在着拉伸畸变,降低了三维垂直地震剖面反射波成像的分辨率。目前,由于缺少针对三维垂直地震剖面资料有效的子波校正手段,在进行三维垂直地震剖面资料处理时,对于拉伸严重的浅层远偏三维垂直地震剖面资料常常采用切除的方式进行处理。
发明内容
本发明的目的是,提供一种三维垂直地震剖面反射波拉伸校正方法,以消除地震子波对反射波的拉伸畸变影响,特别是解决拉伸畸变效应造成的远偏移距资料频率大大降低的问题。
本发明根据三维垂直地震剖面数据中零偏移距(或近偏移距)与远偏移距反射波的速度确定零偏移距道(或近偏移距道)与远偏移距道的反射点的对应关系,选取一定的时窗,采用最小平方法,进行最佳逼近,求取反算子,并将反算子应用到相应的远偏移距垂直地震剖面(VSP)道上,消除拉伸畸变对反射波成像的影响。
本发明具体技术方案是,一种三维垂直地震剖面反射波拉伸校正方法,其特征在于,包括如下步骤:
步骤1,数据选择
对三维垂直地震剖面观测资料中的共检波点道集资料进行分析,如果存在零偏移距道,选出零偏移距垂直地震剖面反射波的速度资料,如果不存在零偏移距道,选出近偏移距垂直地震剖面反射波的速度资料代替零偏移距垂直地震剖面反射波的速度资料;
步骤2,确定反射点对应关系
如果步骤1中选取的是零偏移距垂直地震剖面的反射波速度资料,利用选中的零偏移距垂直地震剖面的反射波速度资料与共检波点道集上的远偏移距垂直地震剖面的反射波速度资料,建立零偏移距道与远偏移距道反射点的对应关系,其反射波时距曲线方程为:
其中,t为反射波旅行时,t0为零偏移距反射波旅行时,v为反射波在地面到地下某一反射层H之间的平均传播速度,x为炮点位置;
求解公式(1)可得
即确立了零偏移距道与非零偏移距道反射点的对应关系;
如果步骤1中选取的是近偏移距垂直地震剖面的反射波速度资料,利用选中的近偏移距垂直地震剖面的反射波速度资料与共检波点道集上的远偏移距垂直地震剖面的反射波速度资料,建立近偏移距道与远偏移距道的反射点的对应关系,其反射波时距曲线方程为:
其中,t1为近偏移距反射波旅行时,x1为其相对应炮点位置,v1为近偏移距的平均速度;tn为第n道非零偏移距反射波旅行时,xn为其相对应炮点位置,vn为第n道非零偏移距的平均速度。
由上式,可得
由于垂直地震剖面(VSP)成像范围较小,在垂直地震剖面(VSP)成像范围内(即最大井源距范围的一半),针对地下某一反射界面,可忽略成像范围内的速度横向变化,则
v1=vn (5)
由公式(2)、(4)可得
即确立了近偏移距道与非零偏移距道反射点的对应关系;
步骤3,反算子计算
根据步骤2确立的零偏移距道或近偏移距道与非零偏移距道反射点的对应关系,选取一定的时窗范围,采用最小平方法,计算反算子,具体方法是:
设计反算子f(t),使实际和期望输出之间的最小平方误差最小,定义误差L为
其中dt为期望输出,yt为实际输出;t为每道记录的时间;
实际输出是反算子与输入的褶积,即有
yt=∑fτxt-τ (8)
将式(8)带入式(7),得
取其偏导数,令其为零,得到方程
其中,i为1,2,3,…,n;
利用Levinson递归算法解上述方程,可得到反算子f(t);
步骤4,远偏移距VSP道校正
按照公式
y′t=f-1(t)*yt (11)
计算出修正后的输出y′t,用修正后输出y′t代替实际输出yt,即完成了远偏移距VSP道的校正。
进一步,所述近偏移距是指,偏移距在100米范围之内的道集。
进一步,所述远偏移距是指,偏移距在100米到最大井源距的一半范围内的道集。
进一步,所述近浅层是指,根据拉伸率计算公式 当拉伸率大于1.3时对应的深度。
本发明的有益效果是:利用三维垂直地震剖面中的零偏移距(或近偏移距)资料对远偏移距的反射波进行校正,可以使得浅层以及远偏移距拉伸较大的同相轴得到校正,有效改善成像的分辨率,提高三维垂直地震剖面反射成像的质量。有效地消除地震子波对反射波的拉伸畸变影响,特别是解决拉伸畸变效应造成的远偏移距资料频率大大降低的问题。本发明有效利用零偏移距道(或近偏移距道)资料的高精度的特点,并用其校正远偏移距道的资料,相比现有技术中采用切除的方式进行处理,无论成像精度和资料信噪比都得到了大大提高。
附图说明
图1是本发明的技术流程图;
图2是垂直地震剖面反射波传播示意图;
图3a为原始共检波点道集数据图;
图3b为校正后的共检波点道集数据图;
图4a为图3a的A部分的局部放大图;
图4b为图3b的A’部分的局部放大图;
图5a为图3a的B部分的局部放大图;
图5b为图3b的B’部分的局部放大图;
图6a是反射子波校正前垂直地震剖面水平叠加(VSPCDP)的效果图;
图6b是反射子波校正后垂直地震剖面水平叠加(VSPCDP)的效果图。
其中:O-井口,X-炮点位置,h-检波点位置,H-地下某一反射层,R-反射点位置。
具体实施方式
以下,结合附图详细说明本发明。
实施例1。
图1出示了本发明的技术流程。
具体步骤是:
(一)确定零(近)偏移距与远偏移距的反射点关系
其方法具体是:
如图2所示VSP反射波传播示意图,O为井口,X为炮点位置,h为检波点位置,H为地下某一反射层,R为反射点位置。设反射波旅行时为t,零偏移距反射波旅行时为t0,从地面到地层H的平均速度为v,则VSP反射波时距曲线方程为
求解公式(1)可得
但是,在实际处理过程中,常常用近偏移距道近似代替零偏移距道。
对于近偏移距道与非零偏移距道,有
其中,t1为近偏移距反射波旅行时,x1为其相对应炮点位置,v1为近偏移距的平均速度;tn为第n道非零偏移距反射波旅行时,xn为其相对应炮点位置,vn为第n道非零偏移距的平均速度。
由上式,可得
由于VSP成像范围较小,在VSP成像范围内(即最大井源距范围的一半),针对地下某一反射界面,可忽略成像范围内的速度横向变化,则
v1=vn (5)
由(3)、(4)可得
即确立了近偏移距道与非零偏移距道反射点的对应关系。
(二)反算子计算
根据近偏移距道与非零偏移距道反射点的对应关系,选取一定的时窗范围,采用最小平方法,计算反算子。
设计反算子f(t),使实际和期望输出之间的最小平方误差最小,定义误差L为
其中dt为期望输出,yt为实际输出。
实际输出是反算子与输入的褶积,即有
yt=∑fτxt-τ (8)
带入式(6),得
取其偏导数,令其为零,得到
利用Levinson递归算法解上述方程,可得到反算子f(t)。
(三)远偏移距VSP道校正
将计算出的反算子f(t)作用于远偏移距VSP道,即
y′t=f-1(t)*yt (11)
即完成了远偏移距VSP道的校正。
试验验证
使用本方法对大庆地区B2-1-031井的3DVSP资料进行了反射子波校正处理,获得了较好的效果,检验了方法的正确性、有效性和稳定性。反射子波校正前后对比如图3a、图3b,图4a、图4b,图5a、图5b、图6a、图6b所示。从图中可以看出,反射子波校正处理后,浅层拉伸较大的同相轴得到校正。通过3DVSP反射波拉伸校正,可以使得浅层(根据拉伸率计算公式 当拉伸率大于1.3时对应的深度都可以作为浅层深度)以及远偏移距(最远范围不能超过最大井源距的一半)拉伸较大的同相轴得到校正,拉伸效应得到压制,道集质量和VSPCDP成像质量均得到改善,有效改善并提高了成像的分辨率。经过本方法处理后,拉伸率20%以内的反射点均能够得到较好的校正。
Claims (4)
1.一种三维垂直地震剖面反射波拉伸校正方法,其特征在于:包括如下步骤:
步骤1,数据选择
对三维垂直地震剖面观测资料中的共检波点道集资料进行分析,如果存在零偏移距道,选出零偏移距垂直地震剖面反射波的速度资料,如果不存在零偏移距道,选出近偏移距垂直地震剖面反射波的速度资料代替零偏移距垂直地震剖面反射波的速度资料;
步骤2,确定反射点对应关系
如果步骤1中选取的是零偏移距垂直地震剖面的反射波速度资料,利用选中的零偏移距垂直地震剖面的反射波速度资料与共检波点道集上的远偏移距垂直地震剖面的反射波速度资料,建立零偏移距道与远偏移距道反射点的对应关系,其反射波时距曲线方程为:
其中,t为反射波旅行时,t0为零偏移距反射波旅行时,v为反射波在地面到地下某一反射层H之间的平均传播速度,x为炮点位置;
求解公式(1)可得
即确立了零偏移距道与非零偏移距道反射点的对应关系;
如果步骤1中选取的是近偏移距垂直地震剖面的反射波速度资料,利用选中的近偏移距垂直地震剖面的反射波速度资料与共检波点道集上的远偏移距垂直地震剖面的反射波速度资料,建立近偏移距道与远偏移距道的反射点的对应关系,其反射波时距曲线方程为:
其中,t1为近偏移距反射波旅行时,x1为其相对应炮点位置,v1为近偏移距的平均速度;tn为第n道非零偏移距反射波旅行时,xn为其相对应炮点位置,vn为第n道非零偏移距的平均速度。
由上式,可得
由于垂直地震剖面成像范围较小,在垂直地震剖面成像范围内(即最大井源距范围的一半),针对地下某一反射界面,可忽略成像范围内的速度横向变化,则
v1=vn (5)
由公式(2)、(4)可得
即确立了近偏移距道与非零偏移距道反射点的对应关系;
步骤3,反算子计算
根据步骤2确立的零偏移距道或近偏移距道与非零偏移距道反射点的对应关系,选取一定的时窗范围,采用最小平方法,计算反算子,具体方法是:
设计反算子f(t),使实际和期望输出之间的最小平方误差最小,定义误差L为
其中dt为期望输出,yt为实际输出;t为每道记录的时间;
实际输出是反算子与输入的褶积,即有
yt=∑fτxt-τ (8)
将式(8)带入式(7),得
取其偏导数,令其为零,得到方程
其中,i为1,2,3,…,n;
利用Levinson递归算法解上述方程,可得到反算子f(t);
步骤4,远偏移距垂直地震剖面道校正
按照公式
y′t=f-1(t)*yt (11)
计算出修正后的输出y′t,用修正后输出y′t代替实际输出yt,即完成了远偏移距垂直地震剖面道的校正。
2.如权利要求1所述的一种三维垂直地震剖面反射波拉伸校正方法,其特征在于:所述近偏移距是指偏移距在100米范围之内的道集。
3.如权利要求1所述的一种三维垂直地震剖面反射波拉伸校正方法,其特征在于:所述远偏移距是指,偏移距在100米到最大井源距的一半范围内的道集。
4.如权利要求1所述的一种三维垂直地震剖面反射波拉伸校正方法,其特征在于:所述近浅层是指,根据拉伸率计算公式 当拉伸率大于1.3时对应的深度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310520617.1A CN104570119B (zh) | 2013-10-29 | 2013-10-29 | 一种三维垂直地震剖面反射波拉伸校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310520617.1A CN104570119B (zh) | 2013-10-29 | 2013-10-29 | 一种三维垂直地震剖面反射波拉伸校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104570119A true CN104570119A (zh) | 2015-04-29 |
CN104570119B CN104570119B (zh) | 2017-11-17 |
Family
ID=53086633
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310520617.1A Active CN104570119B (zh) | 2013-10-29 | 2013-10-29 | 一种三维垂直地震剖面反射波拉伸校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104570119B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106526668A (zh) * | 2016-11-14 | 2017-03-22 | 中国石油化工股份有限公司 | 原始波形提取与成像方法 |
CN107179552A (zh) * | 2016-03-11 | 2017-09-19 | 中国石油化工股份有限公司 | 一种基于波形动态匹配的子波拉伸校正处理方法 |
CN107992646A (zh) * | 2017-11-15 | 2018-05-04 | 中国石油集团东方地球物理勘探有限责任公司 | 垂直地震剖面观测系统炮检点分布范围动态设计方法 |
CN108508487A (zh) * | 2018-03-29 | 2018-09-07 | 中国石油大学(华东) | 一种基于多子波分解的地震道集子波拉伸校正方法和装置 |
CN109490962A (zh) * | 2018-11-30 | 2019-03-19 | 中国海洋大学 | 一种浅地层剖面非零偏移距的消除方法 |
CN111474584A (zh) * | 2020-05-29 | 2020-07-31 | 核工业北京地质研究院 | 一种基于相关型地震干涉的聚焦叠加成像方法及系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002008792A1 (en) * | 2000-07-22 | 2002-01-31 | Services Petroliers Schlumberger | A method of processing vertical seismic profile data using effective vti models |
CN101071175A (zh) * | 2006-05-11 | 2007-11-14 | 中国石油集团东方地球物理勘探有限责任公司 | 零井源距垂直地震剖面纵横波数据深度域走廊叠加剖面处理方法 |
US7508735B2 (en) * | 2006-09-21 | 2009-03-24 | Shell Oil Company | Method of analyzing vertical seismic profile data, method of producing a hydrocarbon fluid, and a computer readable medium |
CN101598811A (zh) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | 一种二维垂直地震剖面数据计算炮点静校正的方法 |
CN101630014A (zh) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | 一种利用垂直地震剖面数据对各向异性介质成像的方法 |
US20120185171A1 (en) * | 2008-12-03 | 2012-07-19 | Yonghe Sun | Multiple anisotropic parameter inversion for a tti earth model using well data |
CN103149592A (zh) * | 2013-03-07 | 2013-06-12 | 天津城市建设学院 | 一种变偏移距vsp波场分离方法 |
-
2013
- 2013-10-29 CN CN201310520617.1A patent/CN104570119B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002008792A1 (en) * | 2000-07-22 | 2002-01-31 | Services Petroliers Schlumberger | A method of processing vertical seismic profile data using effective vti models |
CN101071175A (zh) * | 2006-05-11 | 2007-11-14 | 中国石油集团东方地球物理勘探有限责任公司 | 零井源距垂直地震剖面纵横波数据深度域走廊叠加剖面处理方法 |
US7508735B2 (en) * | 2006-09-21 | 2009-03-24 | Shell Oil Company | Method of analyzing vertical seismic profile data, method of producing a hydrocarbon fluid, and a computer readable medium |
CN101598811A (zh) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | 一种二维垂直地震剖面数据计算炮点静校正的方法 |
CN101630014A (zh) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | 一种利用垂直地震剖面数据对各向异性介质成像的方法 |
US20120185171A1 (en) * | 2008-12-03 | 2012-07-19 | Yonghe Sun | Multiple anisotropic parameter inversion for a tti earth model using well data |
CN103149592A (zh) * | 2013-03-07 | 2013-06-12 | 天津城市建设学院 | 一种变偏移距vsp波场分离方法 |
Non-Patent Citations (3)
Title |
---|
伍敦仕 等: "一种基于共接收点道集的非零偏VSP速度分析方法", 《中国西部科技》 * |
苏媛媛 等: "VSP速度分析及动校正方法", 《中国石油勘探》 * |
魏国华 等: "井间地震广角反射波形校正技术研究与应用", 《石油物探》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107179552A (zh) * | 2016-03-11 | 2017-09-19 | 中国石油化工股份有限公司 | 一种基于波形动态匹配的子波拉伸校正处理方法 |
CN107179552B (zh) * | 2016-03-11 | 2019-06-18 | 中国石油化工股份有限公司 | 一种基于波形动态匹配的子波拉伸校正处理方法 |
CN106526668A (zh) * | 2016-11-14 | 2017-03-22 | 中国石油化工股份有限公司 | 原始波形提取与成像方法 |
CN106526668B (zh) * | 2016-11-14 | 2020-03-31 | 中国石油化工股份有限公司 | 原始波形提取与成像方法 |
CN107992646A (zh) * | 2017-11-15 | 2018-05-04 | 中国石油集团东方地球物理勘探有限责任公司 | 垂直地震剖面观测系统炮检点分布范围动态设计方法 |
CN107992646B (zh) * | 2017-11-15 | 2021-02-12 | 中国石油集团东方地球物理勘探有限责任公司 | 垂直地震剖面观测系统炮检点分布范围动态设计方法 |
CN108508487A (zh) * | 2018-03-29 | 2018-09-07 | 中国石油大学(华东) | 一种基于多子波分解的地震道集子波拉伸校正方法和装置 |
CN109490962A (zh) * | 2018-11-30 | 2019-03-19 | 中国海洋大学 | 一种浅地层剖面非零偏移距的消除方法 |
CN111474584A (zh) * | 2020-05-29 | 2020-07-31 | 核工业北京地质研究院 | 一种基于相关型地震干涉的聚焦叠加成像方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN104570119B (zh) | 2017-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104570119A (zh) | 一种三维垂直地震剖面反射波拉伸校正方法 | |
CN105974470B (zh) | 一种多分量地震资料最小二乘逆时偏移成像方法及系统 | |
CN106405651B (zh) | 一种基于测井匹配的全波形反演初始速度模型构建方法 | |
CN102890290B (zh) | 一种起伏地表条件下的叠前深度偏移方法 | |
CN102012521B (zh) | 一种地震储层预测中叠前裂缝的检测方法 | |
CN109738945B (zh) | 一种利用叠前深度偏移成果直接生成构造图的方法 | |
CN109669212B (zh) | 地震数据处理方法、地层品质因子估算方法与装置 | |
CN108196305B (zh) | 一种山地静校正方法 | |
CN104597490A (zh) | 基于精确Zoeppritz方程的多波AVO储层弹性参数反演方法 | |
CN105182408A (zh) | 一种合成地震记录的制作方法和装置 | |
CN102033242A (zh) | 一种深层倾斜裂缝储层地震振幅预测方法 | |
CN101201409B (zh) | 一种地震数据变相位校正方法 | |
CN104459794A (zh) | 共反射点道集时变时间差值的校正方法及装置 | |
CN104536042A (zh) | 一种二维叠后地震资料振幅补偿方法及装置 | |
CN104483704B (zh) | 基于avo异常类型约束的剩余相位校正方法 | |
CN106291693A (zh) | 一种基于广义s变换的叠前q值反演方法及系统 | |
CN103576200A (zh) | 低信噪比区浅层波阻抗界面静校正方法 | |
CN103984011A (zh) | 一种动态q补偿偏移方法 | |
CN103412324B (zh) | 一种估计介质品质因子的epifvo方法 | |
CN104297800B (zh) | 一种自相控叠前反演方法 | |
CN103616723A (zh) | 基于avo特征的crp道集真振幅恢复方法 | |
CN105093301A (zh) | 共成像点反射角角道集的生成方法及装置 | |
CN102565852B (zh) | 针对储层含油气性检测的角度域叠前偏移数据处理方法 | |
Vigh et al. | Essential steps for successful full-waveform inversion using land data | |
CN104459787A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |