CN101630016B - 一种提高垂直地震剖面成像质量的方法 - Google Patents
一种提高垂直地震剖面成像质量的方法 Download PDFInfo
- Publication number
- CN101630016B CN101630016B CN2008101167092A CN200810116709A CN101630016B CN 101630016 B CN101630016 B CN 101630016B CN 2008101167092 A CN2008101167092 A CN 2008101167092A CN 200810116709 A CN200810116709 A CN 200810116709A CN 101630016 B CN101630016 B CN 101630016B
- Authority
- CN
- China
- Prior art keywords
- wave
- vsp
- descending
- interface
- wave field
- 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
Abstract
本发明是石油物探中提高垂直地震剖面成像质量的方法,在零偏及非零偏垂直地震剖面波场数据上拾取不同接收点处的下行纵波初至时间,建立无倾角的初始层速度模型,纵波初至时间外推至井旁点,获得下行波的时间场,求取下、上行波射线角度利用下行波射线角度,采用射线追踪方法形成理论上行波场,获取一系列的对应不同界面倾角值,确定最大对应的界面倾角为所要求取的井旁反射界面倾角。本发明可有效降低上覆地层对反射界面倾角估计结果的影响和一次反射波旅行时拾取误差对反射界面倾角估计结果的影响,提高利用垂直地震剖面成像质量。
Description
技术领域
本发明涉及石油物探技术,是垂直地震剖面(VSP)数据处理方法,具体为提高垂直地震剖面成像质量的方法。
背景技术
常规VSP资料处理通常包括初至拾取、速度反演、波场分离、成像等。VSP资料的成像处理一般基于水平速度模型。然而,在反射界面存在倾角的情况下,基于水平层状速度模型进行VSP成像会使归位不准确。对这一问题通常的考虑是通过倾角时差校正或迭代成像来解决。VSP数据的倾角时差校正方法是基于地下介质为单一均匀倾斜层的假设,并且需要先由倾角测井或VSP资料独立地估算出反射界面的倾角来。VSP数据的迭代成像实际上隐含着对反射界面的倾角估计(通过不断修改速度模型来降低反射界面倾角对成像结果的影响)。如果能在成像前获得反射界面的倾角信息,则可以形成一个合乎逻辑的VSP数据处理思路,即基于速度模型(反射层的几何形态和速度属性)的VSP数据处理。
VSP的波场结构与反射界面的倾角紧密相关,因此,可通过VSP数据来估算反射界面的倾角,用以修正速度模型,来提高VSP成像质量。
Kennett、Wyatt等曾提出利用VSP资料基于解析公式估算反射界面倾角的方法。这一方法实际上是通过测量两个接收深度位置之间的倾角时差来计算反射界面倾角。随后,Malloy、Noponen等将以一道或两道倾斜时差进行反射界面倾角估计变为以多道倾斜时差进行反射界面的倾角估计,提高了倾角估算的准确性。但这类倾斜时差法的理论基础实际上是基于单层倾斜平反射界面,这意味着实际估计反射界面倾角时要使用平均速度,没有考虑到射线偏折的影响,使得估算的反射界面倾角存在较大误差。
为了避免使用平均速度对反射界面倾角估算结果的影响,Lins等提出了基于层速度模型进行反射界面倾角估算的模型比较方法。以零倾角作为层速度模型的初始估计值,计算射线从激发点传播到接收点的旅行时,再比较根据模型计算的旅行时与实际观测的旅行时之差,不断调整模型的倾角估计值,直到旅行时差值达到给出的精度范围内。
为了减少反演计算的变量个数和反演过程中的多解性,Lins采用所谓“剥层”法,从浅层到深层逐步进行反射界面的倾角估算。但用“剥层”法反演深层反射界面的倾角时,浅层误差会累积进去,使得对于深层反射界面倾角的估算结果并不理想。另一方面,相比于初至波旅行时,一次反射波旅行时对于界面倾角的变化敏感,使用一次反射波旅行时进行反射界面倾角反演结果会更精确。但在Lins的模型比较方法中,由于进行比较的是射线旅行时间,这样,多次波的干扰及在界面附近下行直达波的干涉就会影响到一次反射波旅行时的准确拾取,从而影响到反射界面倾角的反演成像结果。
利用VSP数据估算反射界面倾角还有一类偏振测量法(Spencer等),通过对三分量VSP数据进行偏振分析估算出射线的传播方向,再由射线传播方向和反射界面倾角的关系获得反射界面的倾角。这一方法估算反射层倾角的准确性受偏振测量误差的影响很大。但由于波场干涉的影响,要想准确确定波场的偏振方向并非易事。
发明目的
本发明目的是提供一种用修正速度模型,提高利用垂直地震剖面成像质量的方法。
发明内容
本发明的具体实施步骤为:
1、在地表人工激发地震波,在井中放置地震波接收器串接收地震波场信号,再由电缆传到地面记录仪器记录VSP地震波波场;
2、在零偏及非零偏VSP地震波波场数据上拾取不同接收点处的下行纵波初至时间;
3、对VSP地震波波场进行波场分离,得到实际观测的反射波场;
4、根据已知的声波测井资料及零偏VSP下行纵波初至时间建立无倾角的初始层速度模型;
步骤4)所述的建立无倾角的初始层速度模型为通过声波测井资料划分层位,利用零偏VSP下行纵波初至时间反演各层速度。
5、根据下式将非零偏VSP纵波初至时间外推至井旁点,获得下行波的时间场t(x,z):
式中:τ(h;x,z)是用射线追踪方法确定的下行波由深度h处接收点到井旁点(x,z)的旅行时间,t(h)指非零偏VSP深度h处接收点下行纵波初至时间,x指井旁点横坐标,z指井旁点纵坐标,h代表接收点深度,t和τ用来表示时间,max表示找到沿不同射线路径所得一系列t(h)与τ(h;x,z)时间差中最大者;
其中:arcsin表示反正弦函数;
7、利用下行波射线角度α(x,z)求取出上行波射线角度β(x,z):
步骤7)所述的入射速度,反射速度(Vα、Vβ)是纵波(P)波速度并且相等。
8、利用下行波的时间场t(x,z)、上行波射线角度β(x,z),采用射线追踪方法形成理论上行波场
12、用井旁反射界面倾角值对步骤4中所得初始层速度模型进行修正,利用修改后的带有界面倾角信息的层速度模型进行垂直地震剖面成像。
本发明通过外推下行纵波(P)波初至时间以获得接收点附近上行P波时间场的方法,可有效降低上覆地层对反射界面倾角估计结果的影响和一次反射波旅行时拾取误差对反射界面倾角估计结果的影响,提高了利用垂直地震剖面成像质量。
附图说明
图1为用来产生理论VSP数据的二维构造模型。
图2为根据图1给出的模型所产生的理论三分量VSP上行波场。
图3为倾角扫描结果。其中a为倾角扫描图谱,b为根据振幅最大准则在图谱中所拾取的各反射层倾角值。
具体实施方式
以下结合附图进一步说明本发明。
在地表人工激发地震波,在井中放置地震波接收器串接收地震波场信号,再由电缆传到地面记录仪器记录VSP地震波波场。激发位置包括零偏和非零偏两种方式。
本发明实施例中使用的数据来自于一个二维构造模型。图1给出了该构造模型,图2是根据这一模型所产生的理论三分量VSP上行波场数据,炮点到井口的距离为1000米。
分别在所记录的零偏和非零偏两种数据中拾取下行纵波初至。然后,通过声波测井资料划分层位,利用零偏VSP下行纵波初至时间反演各层速度,建立无倾角信息的初始层速度模型,进行倾角扫描时所采用的层速度模型数据如下表1所示;
序号 | X坐标(m) | Y坐标(m) | Z坐标(m) | P波速度(m/s) | S波速度(m/s) | 界面倾角(°) |
1 | 0.0 | 0.0 | 0.0 | 1750.00 | 712.730 | 0.0 |
2 | 0.0 | 0.0 | 300.0 | 1800.00 | 717.42 | 0.0 |
3 | 0.0 | 0.0 | 700.0 | 2200.00 | 1044.80 | 0.0 |
4 | 0.0 | 0.0 | 1200.0 | 1900.00 | 882.85 | 0.0 |
5 | 0.0 | 0.0 | 1800.0 | 2600.00 | 914.15 | 0.0 |
6 | 0.0 | 0.0 | 2100.0 | 3000.00 | 953.88 | 0.0 |
7 | 0.0 | 0.0 | 2400.0 | 3400.00 | 1024.56 | 0.0 |
表中的坐标代表层界面与井轨迹交点处的坐标值。
对非零偏VSP数据进行波场分离,以得到实际观测的反射波场,如图2所示,根据图1给出的模型所产生的理论三分量VSP上行波场中的X、Z分量。纵坐标为深度(米),横坐标为时间(毫秒)。
式中x指井旁点横坐标,z指井旁点纵坐标,h代表接收点深度,t和τ用来表示时间。τ(h;x,z)是用射线追踪方法确定的下行波由深度h处接收点到井旁点(x,z)的旅行时间,射线追踪过程采用的是对初始层速度模型用界面倾角假设值修正后的层速度模型;t(h)指非零偏VSP深度h处接收点下行波初至时间;max表示找到沿不同射线路径所得一系列t(h)与τ(h;x,z)时间差中最大者。
其中arcsin表示反正弦函数。
然后,依据Snell定律,利用下行波射线角度α(x,z)求取出上行波射线角度β(x,z):
改变界面倾角的假设值重复上面的步骤,可以获取一系列的对应不同界面倾角假设值的值,这里称之为倾角扫描谱图,如图3中a所示,图中纵坐标为深度(米),横坐标为角度(度)。在倾角扫描谱图中最大振幅值位置处对应的倾角即为井旁反射界面的倾角。
最后,用依据振幅最大准则在倾角扫描谱图中所拾取的各反射界面倾角值(如图3中b所示)对初始层速度模型进行修正,利用修改后的带有界面倾角信息的层速度模型进行垂直地震剖面成像。
本发明实施例中使用的数据来自于一个二维构造模型。图1给出了该构造模型,图2是根据这一模型所产生的理论三分量VSP上行波场数据,炮点到井口的距离为1000米。本发明使用图2中的理论上行波场数据及表1的层速度模型数据作为输入,对反射界面的倾角进行扫描。倾角扫描结果如图3所示,其中a为扫描图谱,b为根据振幅最大准则在图谱中所拾取的各反射层倾角值。将图3中b所示的倾角拾取结果和图1给出的理论模型中各层倾角进行对比,可知拾取的倾角值是准确的,表明本发明是行之有效的。
Claims (3)
1.一种提高垂直地震剖面成像质量的方法,其特征在于具体实施步骤为:
1)在地表人工激发地震波,在井中放置地震波接收器串接收地震波场信号,再由电缆传到地面记录仪器记录VSP地震波波场;
2)在零偏及非零偏VSP地震波波场数据上拾取不同接收点处的下行纵波初至时间;
3)对VSP地震波波场进行波场分离,得到实际观测的反射波场;
4)根据已知的声波测井资料及零偏VSP下行纵波初至时间建立无倾角的初始层速度模型;
所述的建立无倾角的初始层速度模型为通过声波测井资料划分层位,利用零偏VSP下行纵波初至时间反演各层速度;
5)根据下式将非零偏VSP纵波初至时间外推至井旁点,获得下行波的时间场t(x,z):
式中:τ(h;x,z)是用射线追踪方法确定的下行波由深度h处接收点到井旁点(x,z)的旅行时间,t(h)指非零偏VSP深度h处接收点下行纵波初至时间,x指井旁点横坐标,z指井旁点纵坐标,h代表接收点深度,t和τ用来表示时间,max表示找到沿不同射线路径所得一系列t(h)与τ(h;x,z)时间差中最大者;
其中:arcsin表示反正弦函数;
7)利用下行波射线角度α(x,z)求取出上行波射线角度β(x,z):
式中:sin表示正弦函数;Vα为入射速度,Vβ为反射速度;φ为界面倾角假设值,x、z分别为计算点处的水平坐标和垂直坐标;
12)用井旁反射界面倾角值对步骤4)中所得初始层速度模型进行修正,利用修改后的带有界面倾角信息的层速度模型进行垂直地震剖面成像。
2.根据权利要求1所述的提高垂直地震剖面成像质量的方法,其特征在于步骤5)所述射线追踪过程采用的是对初始层速度模型用界面倾角假设值修正后的层速度模型。
3.根据权利要求1所述的提高垂直地震剖面成像质量的方法,其特征在于步骤7)所述的入射速度Vα、反射速度Vβ是纵波速度并且相等。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008101167092A CN101630016B (zh) | 2008-07-16 | 2008-07-16 | 一种提高垂直地震剖面成像质量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008101167092A CN101630016B (zh) | 2008-07-16 | 2008-07-16 | 一种提高垂直地震剖面成像质量的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101630016A CN101630016A (zh) | 2010-01-20 |
CN101630016B true CN101630016B (zh) | 2011-10-05 |
Family
ID=41575188
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008101167092A Active CN101630016B (zh) | 2008-07-16 | 2008-07-16 | 一种提高垂直地震剖面成像质量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101630016B (zh) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102313901B (zh) * | 2010-06-29 | 2014-01-15 | 中国石油天然气集团公司 | 一种初至波迭代拾取的方法 |
CN102053267B (zh) * | 2010-10-22 | 2012-11-07 | 中国石油化工股份有限公司 | 一种地震剖面资料处理中基于参数反演的vsp波场分离方法 |
CN102323620B (zh) * | 2011-07-29 | 2014-08-06 | 中国石油化工股份有限公司 | 一种利用钻井轨迹校正vsp上行波的方法 |
CN102393533A (zh) * | 2011-07-29 | 2012-03-28 | 中国石油化工股份有限公司 | 一种利用钻井轨迹校正vsp上行转换波的处理方法 |
CN102411155A (zh) * | 2011-07-29 | 2012-04-11 | 中国石油化工股份有限公司 | 一种利用钻井轨迹校正vsp下行透射转换波的处理方法 |
CN102914792B (zh) * | 2011-08-05 | 2015-10-07 | 中国石油化工股份有限公司 | 一种提高非零偏移距vsp三分量资料的成像效果的方法 |
CN102788993B (zh) * | 2012-08-13 | 2017-02-08 | 中国石油大学(华东) | 最优偏移孔径与免子波畸变有机融合的3d vsp叠前成像方法 |
CN104375176B (zh) * | 2013-08-15 | 2017-05-10 | 中国石油天然气集团公司 | 一种利用初至时间误差确定同相叠加最高截至频率的方法 |
CN103615239B (zh) * | 2013-11-20 | 2016-05-11 | 中国石油天然气集团公司 | 一种测井地层成像系统及其方法 |
CN105301639B (zh) * | 2015-10-21 | 2017-09-01 | 中国石油天然气集团公司 | 基于vsp旅行时双加权层析反演速度场的方法及其装置 |
CN105911586B (zh) * | 2016-04-13 | 2018-10-16 | 中国石油天然气集团公司 | 一种垂直地震数据成像方法及装置 |
CN107340541A (zh) * | 2017-07-10 | 2017-11-10 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种叠前深度偏移速度建模方法及其反射点优选方法 |
CN108845350B (zh) * | 2018-03-15 | 2020-04-10 | 中国石油天然气集团有限公司 | 反演二维速度模型的方法及装置 |
CN108693562B (zh) * | 2018-05-16 | 2019-08-06 | 中国石油大学(华东) | 一种用电成像测井测定古流向的方法 |
CN108828664B (zh) * | 2018-06-07 | 2019-12-10 | 中国石油天然气股份有限公司 | 一种多次波识别方法及装置 |
CN109100792B (zh) * | 2018-10-31 | 2021-06-22 | 中国石油化工股份有限公司 | 基于台站与三维地震联合采集资料的速度反演方法 |
CN113703040B (zh) * | 2021-08-30 | 2023-12-01 | 中油奥博(成都)科技有限公司 | 零偏垂直地震剖面多次波属性确定方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6766253B2 (en) * | 2001-03-14 | 2004-07-20 | Witten Technologies Inc. | Method for merging position information with measurements and filtering to obtain high-quality images that are positioned accurately with respect to global coordinates |
CN101109821A (zh) * | 2007-08-16 | 2008-01-23 | 中国石化集团胜利石油管理局 | 基于系统辨识提高地震资料分辨率的方法 |
-
2008
- 2008-07-16 CN CN2008101167092A patent/CN101630016B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6766253B2 (en) * | 2001-03-14 | 2004-07-20 | Witten Technologies Inc. | Method for merging position information with measurements and filtering to obtain high-quality images that are positioned accurately with respect to global coordinates |
CN101109821A (zh) * | 2007-08-16 | 2008-01-23 | 中国石化集团胜利石油管理局 | 基于系统辨识提高地震资料分辨率的方法 |
Non-Patent Citations (1)
Title |
---|
王宝彬等.利用VSP 资料估算反射界面倾角方法综述.《勘探地球物理进展》.2008,第31卷(第1期),16-24. * |
Also Published As
Publication number | Publication date |
---|---|
CN101630016A (zh) | 2010-01-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101630016B (zh) | 一种提高垂直地震剖面成像质量的方法 | |
CN104656142B (zh) | 一种利用垂直地震剖面与测井联合的地震层位标定方法 | |
CN102967882B (zh) | 地层的层速度模型的建模方法 | |
CN109738945B (zh) | 一种利用叠前深度偏移成果直接生成构造图的方法 | |
CN102841379B (zh) | 一种基于共散射点道集的叠前时间偏移与速度分析方法 | |
CN109669212B (zh) | 地震数据处理方法、地层品质因子估算方法与装置 | |
CN102879819B (zh) | 保持地震波场运动学特征的地震数据处理方法及装置 | |
CN106226818A (zh) | 地震数据处理方法和装置 | |
CN105093301B (zh) | 共成像点反射角角道集的生成方法及装置 | |
CN105182408A (zh) | 一种合成地震记录的制作方法和装置 | |
CN102053261A (zh) | 一种地震数据处理方法 | |
CN109839660A (zh) | 一种利用叠前道集数据建立速度深度模型的方法 | |
CN105093281A (zh) | 一种反演框架下的地震多波建模方法 | |
CN104122582B (zh) | 利用叠加速度求取高精度地震波速度的方法 | |
CN105093319A (zh) | 基于三维地震数据的地面微地震静校正方法 | |
CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
CN107656308B (zh) | 一种基于时间深度扫描的共散射点叠前时间偏移成像方法 | |
CN107515420A (zh) | 一种用于局部相关同相轴的走时与梯度精确拾取方法 | |
CN103149588A (zh) | 一种利用井震标定计算vti各向异性参数的方法及系统 | |
CN111123359B (zh) | 随钻测井与地层格架约束的井周地震成像探测方法及装置 | |
CN106324682A (zh) | 一种应用于冻土层地区的表层结构调查方法 | |
CN102565852B (zh) | 针对储层含油气性检测的角度域叠前偏移数据处理方法 | |
CN104570073A (zh) | 一种适用于复杂高陡构造的双反射地震波成像方法 | |
CN104977615A (zh) | 一种基于模型统计拾取的深水obc资料多次波压制方法 | |
CN111158050B (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 |