CN102162858A - 利用非对称走时进行动校正速度分析的方法 - Google Patents
利用非对称走时进行动校正速度分析的方法 Download PDFInfo
- Publication number
- CN102162858A CN102162858A CN2010105823810A CN201010582381A CN102162858A CN 102162858 A CN102162858 A CN 102162858A CN 2010105823810 A CN2010105823810 A CN 2010105823810A CN 201010582381 A CN201010582381 A CN 201010582381A CN 102162858 A CN102162858 A CN 102162858A
- Authority
- CN
- China
- Prior art keywords
- velocity
- walk
- asymmetric
- nmo
- travel time
- 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.)
- Pending
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种利用非对称走时进行动校正速度分析的方法,包括以下步骤:1)对于一个给定的地质区域和相应的地震记录,首先利用共中心点道集进行常规速度估计,作为该区域的初始估计速度;2)根据非对称走时公式计算走时;3)根据所得的走时结果进行非对称走时叠前时间偏移,得到共成像点道集;4)利用非对称走时动校正重新估计速度场模型;5)如果更新后的速度与上一步的速度变化超过设定的阈值,重复进行步骤2~4,继续进行速度场更新;6)如果更新的速度场与上一步的速度变化不超过设定的阈值或达到一定的迭代次数,则终止迭代,得到最终动校正速度。本方法更适用于复杂地下介质背景的油气勘探需要。
Description
技术领域
本发明主要涉及石油勘探、天然气勘探、煤田勘探、矿产勘探等地震勘探领域,是一种提高地震速度分析精度的新方法,对于后续地震偏移成像、构造成图及岩性分析都具有重大影响。
背景技术
地震勘探经过几十年的发展,浅层,大型构造,横向均匀介质背景下的油气藏陆续被发现,勘探领域逐渐向深层,隐伏构造,横向非均匀介质背景下的油气藏方向拓展,其中面临的一个最重要问题就是地震动校正速度分析问题。
速度分析是地震资料处理中的一个重要环节,是实现动校正及叠加的基本保障,是储层预测和反演的基础。Taner等1969年首次提出了常规叠加速度谱速度分析方法,之后速度分析技术得到了飞速发展。Davis1972年基于Taner叠加速度谱速度分析方法,探讨了速度谱的拾取问题,并第一次将速度谱绘制成能量团形式;May等1979年提出了高阶动校正速度谱,至今该方法仍在长排列资料处理中发挥作用;1989年,Biondi等提出了特征值速度分析方法;之后Key等1990年对该方法进行了改进;Kitdin1992年对多道相关法和特征值法进行了分析和比较;林小竹1993年实现了协方差速度谱;Morozov等1996年根据概率统计和假设检验理论,提出了统计性相位相关速度谱法,其数理基础不同于相似法和特征值法;王德利等2001年很好地实现了Morozov等提出的统计性相位相关速度谱法。Sarkar等2002年引入了一种较好的相干窗设置与计算方法,解决了振幅存在AVO现象时速度谱的计算问题;Larner等2007年将互相关作为控制参数,在最小和最大偏移距之间选取部分分析数据,提高了速度谱的准确度。
但是以上速度分析方法大都是基于层状介质、横向均匀介质假设,而这种思想已经逐渐不能满足日益复杂的地下介质背景的油气勘探需要。
发明内容
为解决传统方法不能满足复杂介质背景下的油气勘探动校正速度分析的问题,本发明目的是试图提供一种更准确的非对称走时动校正速度分析方法,使得动校正速度分析的结果在横向变速条件下更加准确,成像更加清晰,非均质特性更加明显,从而为后期的储层预测及井位部署提供更高精度的基础数据。
为实现上述发明目的,本发明采取的技术方案是:一种利用非对称走时进行动校正速度分析的方法,其特征在于包括以下步骤:
1)对于一个给定的地质区域和相应的地震记录,首先利用共中心点道集进行常规速度估计,作为该区域的初始估计速度;
2)根据公式(1)计算走时,即为非对称走时,
Ti 2(x,y,ζ,x-x′,y-y′)
=t0,0 2+(c2,0(x-x′)2+c1,1(x-x′)(y-y′)+c0,2(y-y′)2)
+(c3,0(x-x′)3+c2,1(x-x′)2(y-y′)+c1,2(x-x′)(y-y′)2+c0,3(y-y′)3)+……
(1)
其中,Ti(x,y,ξ,x-x′,y-y′)是炮点位于(x′,y′)时(x,y)处的走时,是关于(x-x′,y-y′)的多项式,系数ci,j(x)是(x,y,ξ)的函数;ξ为射线坐标系的拟深度,即时间域走时,(x,y)是CDP点的位置,(x′,y′)是炮点位置;
3)根据所得的走时结果进行非对称走时叠前时间偏移,得到共成像点道集;
4)利用非对称走时动校正重新估计速度场模型;
5)如果更新后的速度与上一步的速度变化超过设定的阈值,重复进行步骤2~4,继续进行速度场更新;
6)如果更新的速度场与上一步的速度变化不超过设定的阈值或达到一定的迭代次数,则终止迭代,得到最终动校正速度。
在步骤2)中,公式(1)是走时的精确表达式,在实际计算中,对公式(1)截取其前4次项:
Ti 2(x,y,ζ,x-x′,y-y′)
=t0,0 2+(c2,0(x-x′)2+c1,1(x-x′)(y-y′)+c0,2(y-y′)2)
+(c3,0(x-x′)3+c2,1(x-x′)2(y-y′)+c1,2(x-x′)(y-y′)2+c0,3(y-y′)3)
+(c4,0(x-x′)4+c3,1(x-x′)3(y-y′)+c2,2(x-x′)2(y-y′)2+c1,3(x-x′)(y-y′)3+c0,4(y-y′)4
(2)
其中,
在步骤4)的具体做法是,在上一步结果的基础上计算非对称走时动校正分析中的能量谱,通过速度扫描获取对应位置上的速度,从而进行非对称走时动校正速度分析。
在步骤5)中通过不断的迭代进行非对称走时动校正速度分析,提高速度分析的精度。
本发明由于采取以上技术方案,所获得的有益效果是:1、地震勘探中常规动校正方法是基于水平层状假设,利用单参数速度进行动校正速度分析;非对称走时动校正方法依据象征理论,获取了速度、速度的一阶导和二阶导三个参数进行非对称走时动校正速度分析,适用于横向变速介质情况。2、利用非对称走时动校正速度分析方法进行速度分析,使得速度扫描时,速度谱的能量团更集中,求取的速度场更适于后续的时间偏移成像、储层预测及井位部署提供更高精度的基础数据。3、利用非对称走时动校正速度分析得到的速度模型进行地震数据的偏移成像处理,其偏移剖面陡倾角构造分辨率和信噪比更高,同相轴连续性好、断点清楚、断层清晰、地层归位正确、接触关系清楚;中深层信号共像点道集同相轴更长,相位更一致。
附图说明
图1是非对称走时动校正方法原理图,图中曲线a为常规动校正速度分析方法的走时曲线,曲线b曲线为非对称走时方法求取的对应关系;
图2是非对称走时动校正速度分析速度谱;
图3是非对称走时动校正速度分析速度场;
图4是利用非对称走时动校正第一次迭代后的速度偏移得到的共成像点道集;
图5是利用非对称走时动校正第二次迭代后的速度偏移得到的共成像点道集;
图6是利用非对称走时动校正第三次迭代后的速度偏移得到的共成像点道集;
图7是利用非对称走时动校正第一次迭代后的速度进行偏移的结果;
图8是利用非对称走时动校正第二次迭代后的速度进行偏移的结果;
图9是利用非对称走时动校正第三次迭代后的速度进行偏移的结果。
具体实施方式
下面结合附图和实施例,具体说明本发明。
非对称走时是利用速度横向导数计算地震波传播的走时,算法原理包括格林函数走时计算(非均匀网格速度求导、交换算子系数计算、走时系数计算),基于此对未经动校正和叠加的地震数据进行偏移成像。非对称走时计算是从波场延拓的单平方根方程出发,利用两个变换:时间空间域到频率波数域的正反变换和向量场到指数流形上的正反变换,求得时间空间的格林函数解。第一个变换的具体形式是算子的频率波数域表达和傅里叶积分算子,第二个变换的具体形式是向量场相关和指数映射;计算的基础是利用速度的函数和横向导数,通过速度之间相乘和垂向积分,产生新的函数,这个过程是通过交换运算完成的;交换运算的基本运算是第一个函数乘以第二个函数的横向导数,减去第二个函数乘以第一个函数的横向导数;整个运算还包括Bernoulli数作为系数进行求和运算;这些计算将求出传统走时多项式所不含有的3次项系数和5次项系数;根据这些系数计算克希霍夫积分格林函数的走时,其具体计算公式为:
Ti 2(x,y,ζ,x-x′,y-y′)
=t0,0 2+(c2,0(x-x′)2+c1,1(x-x′)(y-y′)+c0,2(y-y′)2)
+(c3,0(x-x′)3+c2,1(x-x′)2(y-y′)+c1,2(x-x′)(y-y′)2+c0,3(y-y′)3)+……
(1)
Ti(x,y,ξ,x-x′,y-y′)是炮点位于(x′,y′)时(x,y)处的走时,是关于(x-x′,y-y′)的多项式,其相位系数ci,j(x)是(x,y,ξ)的函数;ξ为射线坐标系的拟深度,即时间域走时,(x,y)是CDP点的位置,(x′,y′)是炮点位置。
公式(1)是走时的精确表达式,有无穷多项,在实际计算中,对公式(1)截取其前4次项:
Ti 2(x,y,ζ,x-x′,y-y′)
=t0,0 2+(c2,0(x-x′)2+c1,1(x-x′)(y-y′)+c0,2(y-y′)2)
+(c3,0(x-x′)3+c2,1(x-x′)2(y-y′)+c1,2(x-x′)(y-y′)2+c0,3(y-y′)3)
+(c4,0(x-x′)4+c3,1(x-x′)3(y-y′)+c2,2(x-x′)2(y-y′)2+c1,3(x-x′)(y-y′)3+c0,4(y-y′)4)
(2)
其中,
在克希霍夫积分叠前偏移方法的实施中,走时多项式的系数在处理地震数据之前算好存储起来。在处理地震数据时,利用这些系数“现算”走时,用于非对称走时动校正。利用非对称走时来进行动校正速度分析,首先利用计算的走时通过克希霍夫积分对未经动校正和叠加的地震数据进行偏移,然后根据偏移结果反复迭代求取动校正速度。
非对称走时计算得到了3次、5次等奇数次项,传统的走时公式只有偶数次项,公式(1)、(2)中奇数次项的系数与速度的横向导数有关,能够反应速度的横向变化;在横向不变速情况下,非对称走时公式等价于4阶动校正公式。如图1所示,是非对称走时动校正方法原理图。图中曲线a为常规动校正速度分析方法的走时曲线,由图可知走时平方与偏移距平方间呈线性关系。曲线b为非对称走时方法求取的对应关系,由图可知,上述关系呈非线性,这更符合横向非均匀介质走时特征。
非对称走时动校正速度分析的具体步骤为:
1、对于一个给定的地质区域和相应的地震记录,首先利用共中心点(CMP)道集进行常规速度估计,作为该区域的初始估计速度;
2、根据公式(1)尤其是公式(2)计算走时,即为非对称走时;
3、根据所得的走时结果进行非对称走时叠前时间偏移,得到共成像点道集(如图4~6);
4、利用非对称走时动校正重新估计速度场模型,具体做法是在上一步结果的基础上计算非对称走时动校正分析中的能量谱(如图2),通过速度扫描获取对应位置上的速度,从而进行非对称走时动校正速度分析(如图3);
5、如果更新后的速度与上一步的速度变化超过设定的阈值,重复进行步骤2~4,继续进行速度场更新;
图4是利用非对称走时动校正第一次迭代后的速度偏移得到的共成像点道集,道集向下弯曲,说明仍需要进一步迭代来调整速度;图5是利用非对称走时动校正第二次迭代后的速度偏移得到的共成像点道集。道集向下上弯曲,说明仍需要进一步迭代来调整速度;图6是利用非对称走时动校正第三次迭代后的速度偏移得到的共成像点道集。可见道集被拉直,说明迭代速度已经基本趋于真实值。
每次速度场的更新都对偏移剖面的成像效果有较大改善(如图7~9),图7是利用非对称走时动校正第一次迭代后的速度进行偏移的结果;图8是利用非对称走时动校正第二次迭代后的速度进行偏移的结果,与第一次速度迭代后的偏移剖面相比,椭圆框内偏移结果得到了改进,同相轴更连续,能量更集中。图9是用第三次迭代后的速度进行偏移的结果,其成像质量更好。通过椭圆框区域的结果比较,可见利用非对称走时动校正分析得到的速度进行偏移,具有更好的成像效果。
6、如果更新的速度场与上一步的速度变化不超过设定的阈值或达到一定的迭代次数,则终止迭代,得到最终较为精确的动校正速度。
上述方法中,地震记录的振幅变化记录与时间有关,振幅变化范围也很大,几何扩散、传播过程中的能量损失、层间多次反射及非弹性损耗,都是振幅衰减的原因。在进行振幅处理中我们主要从球面扩散振幅补偿、Q补偿、地表一致性振幅补偿、透射损失补偿等几个方面进行处理。
利用非对称走时动校正速度分析方法进行速度分析,使得速度扫描时,能量团更集中,求取的速度更接近于真实速度。从而为后期的时间偏移成像、储层预测及井位部署提供更高精度的基础数据。
地震勘探中常规动校正方法是基于水平层状假设,利用单参数速度进行动校正分析,求取动校正量。非对称走时动校正方法依据象征理论,获取了速度、速度的一阶导和二阶导三个参数进行速度分析,适用于横向变速介质情况。
通常速度模型都是运用水平层状假设下的对称走时叠前时间偏移作为工具来迭代建立的,当速度发生横向变化时,往往速度谱存在焦散现象。运用非对称走时叠前时间偏移,通过不断修改完善程序,反复处理叠前地震资料,利用含速度横向导数的保幅叠前时间偏移,提高了聚焦分析的精度。
Claims (5)
1.一种利用非对称走时进行动校正速度分析的方法,其特征在于包括以下步骤:
1)对于一个给定的地质区域和相应的地震记录,首先利用共中心点道集进行常规速度估计,作为该区域的初始估计速度;
2)根据公式(1)计算走时,即为非对称走时,
Ti 2(x,y,ζ,x-x′,y-y′)
=t0,0 2+(c2,0(x-x′)2+c1,1(x-x′)(y-y′)+c0,2(y-y′)2)
+(c3,0(x-x′)3+c2,1(x-x′)2(y-y′)+c1,2(x-x′)(y-y′)2+c0,3(y-y′)3)+……
(1)
其中,Ti(x,y,ξ,x-x′,y-y′)是炮点位于(x′,y′)时(x,y)处的走时,是关于(x-x′,y-y′)的多项式,系数ci,j(x)是(x,y,ξ)的函数;ξ为射线坐标系的拟深度,即时间域走时,(x,y)是CDP点的位置,(x′,y′)是炮点位置;
3)根据所得的走时结果进行非对称走时叠前时间偏移,得到共成像点道集;
4)利用非对称走时动校正重新估计速度场模型;
5)如果更新后的速度与上一步的速度变化超过设定的阈值,重复进行步骤2~4,继续进行速度场更新;
6)如果更新的速度场与上一步的速度变化不超过设定的阈值或达到一定的迭代次数,则终止迭代,得到最终动校正速度。
2.如权利要求1所述的利用非对称走时进行动校正速度分析的方法,其特征在于:步骤2)中,公式(1)是走时的精确表达式,在实际计算中,对公式(1)截取其前4次项:
Ti 2(x,y,ζ,x-x′,y-y′)
=t0,0 2+(c2,0(x-x′)2+c1,1(x-x′)(y-y′)+c0,2(y-y′)2)
+(c3,0(x-x′)3+c2,1(x-x′)2(y-y′)+c1,2(x-x′)(y-y′)2+c0,3(y-y′)3)
+(c4,0(x-x′)4+c3,1(x-x′)3(y-y′)+c2,2(x-x′)2(y-y′)2+c1,3(x-x′)(y-y′)3+c0,4(y-y′)4
(2)
其中,
3.如权利要求1所述的利用非对称走时进行动校正速度分析的方法,其特征在于:步骤4)的具体做法是,在上一步结果的基础上计算非对称走时动校正分析中的能量谱,通过速度扫描获取对应位置上的速度,从而进行非对称走时动校正速度分析。
4.如权利要求2所述的利用非对称走时进行动校正速度分析的方法,其特征在于:步骤4)的具体做法是,在上一步结果的基础上计算非对称走时动校正分析中的能量谱,通过速度扫描获取对应位置上的速度,从而进行非对称走时动校正速度分析。
5.如权利要求1或2或3或4所述的利用非对称走时进行动校正速度分析的方法,其特征在于:在步骤5)中通过不断的迭代进行非对称走时动校正速度分析,提高速度分析的精度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010105823810A CN102162858A (zh) | 2010-12-06 | 2010-12-06 | 利用非对称走时进行动校正速度分析的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010105823810A CN102162858A (zh) | 2010-12-06 | 2010-12-06 | 利用非对称走时进行动校正速度分析的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102162858A true CN102162858A (zh) | 2011-08-24 |
Family
ID=44464188
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010105823810A Pending CN102162858A (zh) | 2010-12-06 | 2010-12-06 | 利用非对称走时进行动校正速度分析的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102162858A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104932016A (zh) * | 2015-06-23 | 2015-09-23 | 成都理工大学 | 一种Viterbi-BMC算法的自动速度分析方法 |
CN105652317A (zh) * | 2016-02-02 | 2016-06-08 | 中国矿业大学(北京) | 一种高分辨率地震波成像方法和装置 |
CN106569278A (zh) * | 2015-10-12 | 2017-04-19 | 中国石油化工股份有限公司 | 一种多道相似相干速度谱计算方法 |
CN107942385A (zh) * | 2017-11-15 | 2018-04-20 | 中国石油集团东方地球物理勘探有限责任公司 | 一种精细旅行时计算的动校正方法 |
CN108363102A (zh) * | 2018-01-16 | 2018-08-03 | 山东科技大学 | 一种生成正演地震记录随机噪声的方法及系统 |
CN108983289A (zh) * | 2018-05-02 | 2018-12-11 | 中国石油天然气集团有限公司 | 一种确定地震波走时的方法及装置 |
-
2010
- 2010-12-06 CN CN2010105823810A patent/CN102162858A/zh active Pending
Non-Patent Citations (1)
Title |
---|
刘洪等: "基于横向导数的走时计算方法及其在叠前时间偏移中的应用", 《石油物探》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104932016A (zh) * | 2015-06-23 | 2015-09-23 | 成都理工大学 | 一种Viterbi-BMC算法的自动速度分析方法 |
CN106569278A (zh) * | 2015-10-12 | 2017-04-19 | 中国石油化工股份有限公司 | 一种多道相似相干速度谱计算方法 |
CN106569278B (zh) * | 2015-10-12 | 2018-08-07 | 中国石油化工股份有限公司 | 一种多道相似相干速度谱计算方法 |
CN105652317A (zh) * | 2016-02-02 | 2016-06-08 | 中国矿业大学(北京) | 一种高分辨率地震波成像方法和装置 |
CN107942385A (zh) * | 2017-11-15 | 2018-04-20 | 中国石油集团东方地球物理勘探有限责任公司 | 一种精细旅行时计算的动校正方法 |
CN108363102A (zh) * | 2018-01-16 | 2018-08-03 | 山东科技大学 | 一种生成正演地震记录随机噪声的方法及系统 |
CN108983289A (zh) * | 2018-05-02 | 2018-12-11 | 中国石油天然气集团有限公司 | 一种确定地震波走时的方法及装置 |
CN108983289B (zh) * | 2018-05-02 | 2020-06-09 | 中国石油天然气集团有限公司 | 一种确定地震波走时的方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102707316B (zh) | 复杂地质目标的2.5次三维地震勘探方法 | |
CN102176053B (zh) | 提升波动方程叠前深度偏移成像效果的方法 | |
CN102841379B (zh) | 一种基于共散射点道集的叠前时间偏移与速度分析方法 | |
CN109407151B (zh) | 基于波场局部相关时移的时间域全波形反演方法 | |
CN109669212B (zh) | 地震数据处理方法、地层品质因子估算方法与装置 | |
CN102841375A (zh) | 一种复杂条件下基于角度域共成像点道集的层析速度反演方法 | |
CN102162858A (zh) | 利用非对称走时进行动校正速度分析的方法 | |
CN107272062B (zh) | 一种数据驱动的地下介质q场估计方法 | |
CN104237945B (zh) | 一种地震资料自适应高分辨处理方法 | |
CN102841376A (zh) | 一种基于起伏地表的层析速度反演方法 | |
CN105974479A (zh) | Gpu空间网格体的层析2d/3d各向异性深度域速度建模方法 | |
CN103984011A (zh) | 一种动态q补偿偏移方法 | |
CN108445538A (zh) | 基于反射地震资料建立深度域层q模型的方法和系统 | |
CN110471113A (zh) | 基于非稳态地震资料的反演动校正方法、装置及存储介质 | |
CN103048690B (zh) | 基于最优地震子波提取的快速匹配投影分解的地层反射拾取技术 | |
WO2005026776A1 (en) | Wide-offset-range pre-stack depth migration method for seismic exploration | |
CN103499836A (zh) | 空变-多时窗融合高精度剩余静校正方法 | |
US9798028B2 (en) | System and method for seismic adaptive optics | |
CN102053269A (zh) | 一种对地震资料中速度分析方法 | |
CN107942405A (zh) | 预测砂泥岩薄互层砂体累积厚度的方法 | |
CN105092343A (zh) | 一种基于叠前道集去除薄层调谐效应的方法 | |
CN106257309A (zh) | 叠后地震数据体处理方法及装置 | |
CN111474576B (zh) | 一种保持地层结构的叠前地震道集反演初始模型构建方法 | |
CN106908838A (zh) | 构建目标区域地层非弹性衰减品质因子立体模型的方法 | |
CN105093324B (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 | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20110824 |