亚星游戏官网-www.yaxin868.com

专利名称：利用非对称走时进行动校正速度分析的方法
技术领域：
本发明主要涉及石油勘探、天然气勘探、煤田勘探、矿产勘探等地震勘探领域，是一种提高地震速度分析精度的新方法，对于后续地震偏移成像、构造成图及岩性分析都具有重大影响。
背景技术：
地震勘探经过几十年的发展，浅层，大型构造，横向均勻介质背景下的油气藏陆续被发现，勘探领域逐渐向深层，隐伏构造，横向非均勻介质背景下的油气藏方向拓展，其中面临的一个最重要问题就是地震动校正速度分析问题。速度分析是地震资料处理中的一个重要环节，是实现动校正及叠加的基本保障， 是储层预测和反演的基础。Taner等1969年首次提出了常规叠加速度谱速度分析方法，之后速度分析技术得到了飞速发展。Davisl972年基于Taner叠加速度谱速度分析方法，探讨了速度谱的拾取问题，并第一次将速度谱绘制成能量团形式；May等1979年提出了高阶动校正速度谱，至今该方法仍在长排列资料处理中发挥作用；1989年，Biondi等提出了特征值速度分析方法；之后Key等1990年对该方法进行了改进；Kitdinl992年对多道相关法和特征值法进行了分析和比较；林小竹1993年实现了协方差速度谱；Morozov等1996年根据概率统计和假设检验理论，提出了统计性相位相关速度谱法，其数理基础不同于相似法和特征值法；王德利等2001年很好地实现了 Morozov等提出的统计性相位相关速度谱法。 Sarkar等2002年引入了一种较好的相干窗设置与计算方法，解决了振幅存在AVO现象时速度谱的计算问题；Larner等2007年将互相关作为控制参数，在最小和最大偏移距之间选取部分分析数据，提高了速度谱的准确度。但是以上速度分析方法大都是基于层状介质、横向均勻介质假设，而这种思想已经逐渐不能满足日益复杂的地下介质背景的油气勘探需要。

发明内容
为解决传统方法不能满足复杂介质背景下的油气勘探动校正速度分析的问题，本发明目的是试图提供一种更准确的非对称走时动校正速度分析方法，使得动校正速度分析的结果在横向变速条件下更加准确，成像更加清晰，非均质特性更加明显，从而为后期的储层预测及井位部署提供更高精度的基础数据。为实现上述发明目的，本发明采取的技术方案是一种利用非对称走时进行动校正速度分析的方法，其特征在于包括以下步骤1)对于一个给定的地质区域和相应的地震记录，首先利用共中心点道集进行常规速度估计，作为该区域的初始估计速度；2)根据公式(1)计算走时，即为非对称走时，Ti2 (X，y，ζ，χ-χ'，y-y')= t0j02+(c2j0(x-x' )2+Cljl(x-x' ) (y-y' )+c0j2(y-y' )2)
+(c3,0(x-x' )3+c2j1(x-x' )2(y-y' )+clj2(x-x' ) (y-y' )2+c0j3(y-y' )3)+…… (1)其中，Ti(x，y，ξ，χ-χ'，y-y')是炮点位于(χ' , y')时(x，y)处的走时，是关于(X-X'，y-y')的多项式，系数(^,」00是(x，y，ξ)的函数；ξ为射线坐标系的拟深度，即时间域走时，(χ，y)是CDP点的位置，(χ'，y')是炮点位置；3)根据所得的走时结果进行非对称走时叠前时间偏移，得到共成像点道集；4)利用非对称走时动校正重新估计速度场模型；5)如果更新后的速度与上一步的速度变化超过设定的阈值，重复进行步骤2 4， 继续进行速度场更新；6)如果更新的速度场与上一步的速度变化不超过设定的阈值或达到一定的迭代次数，则终止迭代，得到最终动校正速度。在步骤2)中，公式(1)是走时的精确表达式，在实际计算中，对公式(1)截取其前 4次项Ti2 (X，y，ζ，χ-χ'，y-y')= t0j02+(c2j0(x-x' )2+Cljl(x-x' ) (y-y' )+c0j2(y-y' )2)+(c3,0(x_x' )3+c2j1(x-x' )2(y-y' )+clj2(x-x' ) (y-y' )2+c0j3(y-y' )3)+(c4j0(x-x ' )4+c3a(x-x ' )3 (y-y ' )+c2,2(x_x ' )2(y-y ' )2+c1j3(x-x ')
(y-y' )3+c0,4(y-y' )4
(2)其中，
权利要求
1.一种利用非对称走时进行动校正速度分析的方法，其特征在于包括以下步骤1)对于一个给定的地质区域和相应的地震记录，首先利用共中心点道集进行常规速度估计，作为该区域的初始估计速度；2)根据公式(1)计算走时，即为非对称走时， Ti2 (X, y, ζ , X-X‘ , y-y')=t0j02+(c2j0(x-x' )2+cia(x-x' ) (y-y' )+c0j2(y-y' )2)+ (c3,0(x_x' )3+c2a (x-x' )2 (y-y' )+clj2(x-x' ) (y-y' )2+c0j3(y-y' )3)+......(1)其中，Ti (χ, y，ξ，x-x'，y-y')是炮点位于(χ'，y')时(x，y)处的走时，是关于 (x-x' ,1-1’ )的多项式，系数CiJx)是(x，y，ξ)的函数；ξ为射线坐标系的拟深度，即时间域走时，(X，y)是⑶P点的位置，(X'，1’ )是炮点位置；3)根据所得的走时结果进行非对称走时叠前时间偏移，得到共成像点道集；4)利用非对称走时动校正重新估计速度场模型；5)如果更新后的速度与上一步的速度变化超过设定的阈值，重复进行步骤2 4，继续进行速度场更新；6)如果更新的速度场与上一步的速度变化不超过设定的阈值或达到一定的迭代次数， 则终止迭代，得到最终动校正速度。
2.如权利要求1所述的利用非对称走时进行动校正速度分析的方法，其特征在于步骤2)中，公式(1)是走时的精确表达式，在实际计算中，对公式(1)截取其前4次项Ti2 (x, y, ζ,X-X1，y-y')=t 2+(c L0,0 乂。2,0(X-X')、,!(x-x')(y-y')+c0,2(y-y' )2)+ (c3,0(x_x')3+C2a(x-x')2(y--y' )+cL,2 (χ-χ')(y-y' )2+c0,3(y-y')3)+ (C4,0(x-x')4+C3.,!(x-x')3(y--y' )+c2,:2 (χ-χ')2(y-y' )2+Cl,3(χ-χ')(y-y4(y-y' )4 (2)其中，_^0,0^0,2 二 ^0,0C°'22<pl2 2φ01_ -^0,0^3,0 二 -^0,0^0,3-2^0,0^2,1 -^0,0^2,1 1 二 - 二 -22,1 8 ^02 ^22 0 4钓-2^0 0^1,2 -ψ^ψχ, Cj 2 二 - — ---''2 8 ^20 ^22 4φ20φ12
3.如权利要求1所述的利用非对称走时进行动校正速度分析的方法，其特征在于步骤4)的具体做法是，在上一步结果的基础上计算非对称走时动校正分析中的能量谱，通过速度扫描获取对应位置上的速度，从而进行非对称走时动校正速度分析。
4.如权利要求2所述的利用非对称走时进行动校正速度分析的方法，其特征在于步骤4)的具体做法是，在上一步结果的基础上计算非对称走时动校正分析中的能量谱，通过速度扫描获取对应位置上的速度，从而进行非对称走时动校正速度分析。
5.如权利要求1或2或3或4所述的利用非对称走时进行动校正速度分析的方法，其特征在于在步骤幻中通过不断的迭代进行非对称走时动校正速度分析，提高速度分析的精度。
全文摘要
本发明涉及一种利用非对称走时进行动校正速度分析的方法，包括以下步骤1)对于一个给定的地质区域和相应的地震记录，首先利用共中心点道集进行常规速度估计，作为该区域的初始估计速度；2)根据非对称走时公式计算走时；3)根据所得的走时结果进行非对称走时叠前时间偏移，得到共成像点道集；4)利用非对称走时动校正重新估计速度场模型；5)如果更新后的速度与上一步的速度变化超过设定的阈值，重复进行步骤2～4，继续进行速度场更新；6)如果更新的速度场与上一步的速度变化不超过设定的阈值或达到一定的迭代次数，则终止迭代，得到最终动校正速度。本方法更适用于复杂地下介质背景的油气勘探需要。
文档编号G01V1/36GK102162858SQ20101058238
公开日2011年8月24日 申请日期2010年12月6日 优先权日2010年12月6日
发明者刘国锋, 刘洪 , 夏忠谋, 宋亮, 张云鹏, 朱振宇, 李博, 王小六, 薛东川, 赵伟, 邹振, 郝振江 申请人:中国海洋石油总公司, 中国科学院地质与地球物理研究所, 中海石油研究中心