一种水平界面vti介质转换点ccp道集抽取方法
【专利摘要】本发明公开一种水平界面VTI介质转换点CCP道集抽取方法,首先需推导xp、xs和x三者之间的关系,得到迭代表达式,然后选取计算精度与迭代次数的上限,做循环运算找到符合条件的入射P波相角,从而得到转换点位置的坐标,在满足给定的精度范围内,能快速找到转换点的位置,同时也给出了理论时距曲线,在石油勘探PS转换波数据处理中具有重要的意义。
【专利说明】
一种水平界面VTI介质转换点CCP道集抽取方法
技术领域
[0001] 本发明具体涉及一种水平界面VTI介质转换点CCP道集抽取方法。
【背景技术】
[0002] 近年来随着陆地三维地震和海底地震(0BS)采集和处理技术的提高,多波多分量 研究已经进入生产领域中来。传统的P-P波资料在地下构造成像方面取得了巨大成功, 但在含气储层(Gas-bearing layer)构造成像、液体识别、岩性识别以及裂缝各向异性检 测等方面存在着不足。而P-S转换波在穿越含气储层时衰减相比纵波要小、对各向异性 也相对较纵波敏感,它对于含气储层下成像、提取地下裂缝的信息和岩性识别等方面具有 一定的优势。一般而言,P-S转换波资料的处理需要解决两个问题,其一是时差校正分析 (moveout);其二是共转换点(common converted point)位置提取,其涉及到共转换点道集 CCP抽取的问题。就现有的CCP道集抽取技术而言,要想抽取转换点道集,必须要求取转换 点位置,目前有三种方法:第一种方法就是严格解析方法;第二种是基于泰勒级数展开的 近似解析方法;第三是迭代类方法。在上述三种方法中,方法一的优点是能得到精确解,缺 点是要解四次方程,过程比较烦琐,最后得到的解析表达式复杂以致难以判断哪个解是真 实解,且难以扩展到各向异性介质,因此在运用范围上受到限制;方法二是基于泰勒级数展 开的近似解析方法,Thomsen做了先导性的工作,他将x p(转换点到炮点的水平距离)展开 成x的级数,其中的系数跟介质参数有关,数值实验表明该表达式子在x/z < 1的情况下 比较精确,做些修改也适用于VTI介质的计算。在此思路的基础上,Li和Yuan把这个限制 推广到x/z < 3的情形,扩大了该方法的使用范围,并且扩展到层状VTI介质,为各向异性 介质转换波资料的处理奠定了理论基础。此方法的优点是给出了近似解析解、物理意义明 确,在x/z < 3范围之内,计算结果是较为精确的,但不易扩展到其它较为复杂的TI介质。 方法三是属于射线追踪类型的迭代算法如王维佳等人利用牛顿迭代算法求取各向同性倾 斜界面情形下转换点位置的工作,此类方法的缺点是不能从解析的角度来分析参数对转换 点位置的影响。
【发明内容】
[0003] 本发明针对上述问题,我们提出了用弦截迭代算法计算了单层水平界面VTI介质 P-S转换波转换点的位置和旅行时我们技术效果显示:当介质退化到各向同性介质时,用 弦截迭代算法计算的结果与苑春方等的转换点位置解析结果完全一致;在计算VTI介质情 形时,本文将弦截迭代算法与Thomsen以及Li和Yuan的基于泰勒级数展开的近似解析方 法进行了对比,表明在短排列(x/h< 1)时,三种方法结果基本一致;在长排列时Thomsen 的近似解析结果误差较大,而Li和Yuan的近似解析结果在x/h < 3内与弦截迭代算法基 本一致,在3倍偏移距/深度比之外,Li和Yuan的结果比弦截迭代方法的结果要大些。
[0004] 本发明的目的是通过以下技术方案实现的:
[0005] 1. -种水平界面VTI介质转换点CCP道集抽取方法,其步骤如下:
[0006] 第一步:已知入射P波相角,求入射P波群角
[0007] 假设知道入射P波的相角0 p,现在来求其群角Wp,根据Tsvankin给出的VTI介 质P波相速度严格解析表达式为:
(1)
[0009] Vn代表入射P波相速度、0 p代表入射P波相角、e和S为Thomsen参数,其中:
[0010]贝U入射P波的群角 '的正切函数可表达为;
(2)
[0012] 由(2)式,可得转换点坐标为:
[0013] xp = htan(^p) (3)
[0014] 第二步:求取反射SV波的相角
[0015] 在反射界面处,P波和SV波的相速度和相角满足Snell定理,于是有
(4)
[0017] 方程⑷中,Vsv和0 sv分别是反射SV波的相速度和相角,px为水平慢度,同样根 据Tsvankin的工作,可以得到反射SV波的相速度:
[0019] 若已知入射P波相角0 p,就可以知道入射P波相速度Vp,就能得到水平慢度px,根 据(5)式可可得关于sin 0 sv的一元四次方程,解这个方程并且保证-1 < sin 0 sv < 1成 立,则可得到反射SV波的相角:
[0021] 若已知入射P波相角0 p,就可以知道入射P波相速度Vp,就能得到水平慢度px,根 据(4)式可可得关于sin 0 sv的一元四次方程,解这个方程并且保证-1 < sin 0 sv < 1成 立,则可得到反射SV波的相角0SV,
[0022] 第三步:已知反射SV波的群角,求偏移距与共转换点的关系
[0023] 求出反射SV波的相角后,利用下面的公式(7)可以求取反射SV波的群角, (7)
[0025] 方程(7)中,Wsv为反射SV波的群角。这样可以得到xs
[0026] xs = htan(Wsv) (8)
[0027] 于是就得到了单一水平界面VTI介质偏移距x与xp的对应关系
[0028] x = xp+xs = h[tan(^p)+tan(^sv)] (9)
[0029] 第四步:弦截法得到迭代关系式并最后求解:
[0030] 对于(9)式,有下面的目标函数f(ep)
[0031] f ( 0 p) = xp+xs-x = httanC^^+tanC^^J-x (10)
[0032] 其中x为偏移距(固定值)、xp和xs均与0 p相关,那么已知偏移距x求取转换点 的问题变成求取方程f (e p) = 〇的根e p的问题;
[0033] 对于f( 0 p) = 0,利用弦截迭代方法来求取0 p,其迭代表达式为:
[0035] 先选取两个初始值0。和,求取f(0。)和f(0i),再代入到(11)式,作循
[0036] 环运算并判断绝对值| 0k+1_0k| < 为满足的精度要求,这里取€ = 1. Oe-9),如果满足该精度要求则输出结果0 p,即找到入射P波的相角,如果不满足精度,继 续做循环运算,直到满足精度要求为止。获得了 ep后,就可以根据(1)、(2)和(3)式而得 到转换点位置坐标xp,
[0037] 得到0 p和xp后,就可以得到入射P波群速度和反射SV波群速度F;〖,于是P-S 转换波走时T可表示为
(12)
[0039] 另外还可以得到共中心点与转换点CP之间的距离D :
(13)
[0041] 本发明具有如下有益的效果:
[0042] 本发明给出了单一水平界面VTI介质PS转换波转换点计算的一种方法-弦截法, 在满足给定的精度范围内,能快速找到转换点的位置。同时也给出了理论时距曲线。这个方 法有下面几个优点:第一,适用性很强,方法可以推广到更为复杂的介质,比如平层VTI介 质,任意取向的TTI介质(对称轴在入射面内),以及弱裂隙介质。而其它的方法比如基于 泰勒级数展开的近似方法和严格的解析方法很难作更进一步的扩展;第二,方法的计算公 式简单;第三,计算精度可以保证,由于在计算的过程中,没有引入任何近似条件,计算的结 果在满足精度的情形下是可靠的,在石油勘探PS转换波数据处理中具有重要的意义。
【附图说明】
[0043] 图1为本发明的VTI介质P-SV波反射示意图;
[0044] 图2为本发明的VTI介质x,x/h与转换点之间的关系;其中,a偏移距与转换点位 置的关系,b偏移距,深度比与转换点位置的关系;
[0045] 图3为本发明的单一水平界面VTI介质PS转换波时距曲线对比,其中,a来自于Li 和Yuan (2003)的结果;b弦截迭代方法计算的结果,其中,1、弦截法,2、Li和Yuan (2003), 3、Thomsen (1999)〇
【具体实施方式】
[0046] 下面结合附图对本发明作进一步的说明:
[0047] 实施例1
[0048] 考虑单一水平界面VTI介质模型,如图1所示,选取直角坐标系,其中VTI介质对 称轴n与水平反射界面和地表垂直。源点到接收点之间的距离为X,转换点CP (Converted point)到源点的水平距离为xp,转换点到接收点的水平距离为xs,CMP (Common mid point) 为共中心点。问题归结为已知x、反射界面的深度h,描述VTI介质的3个Thomsen参数[21] 分别为e、S、Y和沿着对称轴方向的传播速度Vp。和Vs。,然后来求取x p和P-S转换波理 论走时T的问题。为了解决上述问题,首先需推导xp、xs和x三者之间的关系,得到迭代表 达式,然后选取计算精度与迭代次数的上限,做循环运算找到符合条件的入射P波相角,从 而得到转换点位置的坐标。下面分几个步骤:
[0049] 第一步:已知入射P波相角,求入射P波群角
[0050] 假设知道入射P波的相角0 p,现在来求其群角Wp,根据Tsvankin给出的VTI介 质P波相速度严格解析表达式为:
[0052] Vp代表入射P波相速度、0 p代表入射P波相角、e和S为Thomsen参数,其中:
[0053] 则入射P波的群角的正切函数可表达为;
(2)
[0055] 由(2)式,可得转换点坐标为:
[0056] xp = htan(^p) (3)
[0057] 源点到接收点之间的距离为x,转换点CP(Converted point)到源点的水平距离 为xp,转换点到接收点的水平距离为xs ;
[0058] 第二步:求取反射SV波的相角
[0059] 在反射界面处,P波和SV波的相速度和相角满足Snell定理,于是有 (4)
[0061] 方程⑷中,Vsv和0 sv分别是反射SV波的相速度和相角,px为水平慢度,同样根 据Tsvankin的工作,可以得到反射SV波的相速度:
[0063] 若已知入射P波相角0 p,就可以知道入射P波相速度Vp,就能得到水平慢度px,根 据(5)式可可得关于sin 0 sv的一元四次方程,解这个方程并且保证-1 < sin 0 sv < 1成 立,则可得到反射SV波的相角:
[0065] 若已知入射P波相角0 p,就可以知道入射P波相速度Vp,就能得到水平慢度px,根 据(4)式可可得关于sin 0 sv的一元四次方程,解这个方程并且保证-1 < sin 0 sv < 1成 立,则可得到反射SV波的相角0SV,
[0066] 第三步:已知反射SV波的群角,求偏移距与共转换点的关系
[0067] 求出反射SV波的相角后,利用下面的公式(7)可以求取反射SV波的群角,
[0069] 方程(7)中,Wsv为反射SV波的群角。这样可以得到xs
[0070] xs = htan(Wsv) (8)
[0071] 于是就得到了单一水平界面VTI介质偏移距x与xp的对应关系
[0072] x = xp+xs = h[tan(^p)+tan(^sv)] (9)
[0073] 第四步:弦截法得到迭代关系式并最后求解:
[0074] 对于(9)式,有下面的目标函数f(0p)
[0075] f ( 0 p) = xp+xs-x = httanC^^+tanC^^J-x (10)
[0076] 其中x为偏移距(固定值)、\和\均与0P相关,那么已知偏移距x求取转换点 的问题变成求取方程f (e p) = 〇的根e p的问题;
[0077] 对于f( 0 p) = 0,利用弦截迭代方法来求取0 p,其迭代表达式为:
[0079] 先选取两个初始值0。和,求取f(0。)和f(0i),再代入到(11)式,作循环运 算并判断绝对值I ek+1-ek| < 为满足的精度要求,这里取€ = i.0e-9),如果满足 该精度要求则输出结果e P,即找到入射p波的相角,如果不满足精度,继续做循环运算,直 到满足精度要求为止。获得了 ep后,就可以根据(1)、(2)和(3)式而得到转换点位置坐 丰不Xp,
[0080] 得到0 p和Xp后,就可以得到入射P波群速度和反射sv波群速度U,于是P-S 转换波走时T可表示为 (12)
[0082] 另外还可以得到共中心点与转换点CP之间的距离D :
[0084] (13)。
【主权项】
1. 一种水平界面VTI介质转换点CCP道集抽取方法,首先需推导Xp、Xs和X三者之间 的关系,得到迭代表达式,然后选取计算精度与迭代次数的上限,做循环运算找到符合条件 的入射P波相角,从而得到转换点位置的坐标,其具体步骤如下: 第一步:已知入射P波相角,求入射P波群角 假设知道入射P波的相角θ p,现在来求其群角Ψρ,根据Tsvankin给出的VTI介质P 波相速度严格解析表达式为: 由(2)式,可得转换点坐标为:xp = htan(Wp) (3) 源点到接收点之间的距离为X,转换点CP(Converted point)到源点的水平距离为xp, 转换点到接收点的水平距离为Xs ; 第二步:求取反射SV波的相角 在反射界面处,P波和SV波的相速度和相角满足Snell定理,于是有(4) 方程⑷中,Vsv和Θ sv分别是反射SV波的相速度和相角,Px为水平慢度,同样根据 Tsvankin的工作,可以得到反射SV波的相速度:(5) 若已知入射P波相角θ p,就可以知道入射P波相速度vp,就能得到水平慢度Px,根据 (5)式可可得关于sin Θ sv的一元四次方程,解这个方程并且保证-I < Sin Θ sv < 1成立, 则可得到反射SV波的相角:(6) 若已知入射P波相角θ p,就可以知道入射P波相速度vp,就能得到水平慢度Px,根据 (4)式可可得关于sin Θ sv的一元四次方程,解这个方程并且保证-I < Sin Θ sv < 1成立, 则可得到反射SV波的相角0sv, 第三步:已知反射SV波的群角,求偏移距与共转换点的关系 求出反射SV波的相角后,利用下面的公式(7)可以求取反射SV波的群角,(7) 方程(7)中,Ψ3ν为反射SV波的群角。这样可以得到Xs xs = htanWsv) (8) 于是就得到了单一水平界面VTI介质偏移距X与xp的对应关系 X = xp+xs = h[tan(Wp)+tan(Wsv)] (9) 第四步:弦截法得到迭代关系式并最后求解: 对于(9)式,有下面的目标函数f(0p) f ( θP) = xP+xs_x = h[tan(Wp)+tan(Wsv)]_x (10) 其中x为偏移距(固定值)、\和xs均与θρ相关,那么已知偏移距χ求取转换点的问 题变成求取方程f ( θ ρ) = 0的根θ ρ的问题; 对于f( Θ J = 〇,利用弦截迭代方法来求取θ η,其迭代表达式为:免=1,2,3…… (11) 先选取两个初始值9。和Q1,求取?·(θ。)和He1),再代入到(11)式,作循环运算并 判断绝对值I 0k+1-0k| < ξ (ξ为满足的精度要求,这里取ξ = l.〇e-9),如果满足该精 度要求则输出结果Θ P,即找到入射P波的相角,如果不满足精度,继续做循环运算,直到满 足精度要求为止。获得了 θρ后,就可以根据(1)、(2)和(3)式而得到转换点位置坐标xp, 得到Θ p和Xp后,就可以得到入射P波群速度K/和反射SV波群速度,于是P-S转换 波走时T可表示为(12) 另外还可以得到共中心点与转换点CP之间的距离D :(13)。
【文档编号】G01V1/28GK105891881SQ201410204442
【公开日】2016年8月24日
【申请日】2014年5月15日
【发明人】蔡晓刚
【申请人】蔡晓刚