第2章-叠前时间偏移
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第二章叠前时间偏移
地震波成像在油气勘探中占据重要位置。
它的作用是使反射波或绕射波返回到产生它们的地下位置,从而得到地下地质构造的精确成像。
从二十世纪60年代偏移过程由计算机实现以来,已从常规偏移即叠后时间偏移发展到了目前的叠前深度偏移。
偏移方法的研究和应用是受油气勘探的实际需求驱动的,同时它又受到人们对偏移成像的认识程度和计算机处理能力的制约。
常规偏移(即叠后时间偏移)在以往的油气勘探过程中起到了重要作用,但随着勘探难度的提高,在构造较为复杂或/和强横向变速的地区,基于常规偏移的处理方法再也难见成效。
究其原因,一方面是由于常规处理是先叠加后偏移,水平叠加过程受水平层状介质假设制约,在复杂地质构造条件下,这种叠加过程很难实现同相叠加,这样会对波场产生破坏,所以用这种失真了的叠后数据去进行偏移处理难以取得好的成像效果就很自然了。
为了克服非同相叠加给后续偏移带来的麻烦,人们提出使用叠前偏移,即先偏移处理使波场归位,再把同一地下点的偏移波场相叠加。
这样,在横向速度中等变化的较为复杂构造成像中叠前时间偏移可以弥补常规偏移的不足。
另一方面是由于时间偏移是建立在均匀介质或水平层状介质的速度模型的基础上的,当速度存在横向变化,或速度分界面不是水平层状的情况下,常规偏移不能满足Snell定律,因此不能进行正确的反射波的偏移成像。
为了解决这个问题,出现了深度偏移。
这样,在强横向变速的一般构造成像中,叠后深度偏移可以弥补常规偏移的不足;而在强横向变速的复杂构造成像中,叠前深度偏移可以弥补常规偏移的不足。
迄今为止,人们已对叠前时间偏移进行了20多年的研究工作,而对叠前深度偏移也进行了十几年的研究和探索工作。
本章重点讨论叠前时间偏移。
叠前深度偏移将在第四章和第五章讨论。
近年来,随着叠前时间偏移方法和技术的不断成熟和与之配套技术的不断完善以及计算机性能的不断提高,实现叠前时间偏移已成为现实。
目前,国内外有多家地球物理处理公司和计算中心已进行叠前时间偏移处理,部分公司还把叠前时间偏移作为常规处理软件加入到常规处理流程中,使之成为常规处理的一个重要内容。
叠前时间偏移技术之所以受到如此重视和关注,主要是因为这种技术相对叠后时间偏移和叠前深度偏移技术有如下的几个特点:1)实现这种技术所需的软硬件成本合理,多家处理公司和计算中心都能接受和承受。
2)叠前时间偏移相对叠前深度偏移而言,对偏移速度场无过高的要求,假设条件少,经对常规法进行简单的改进或/和修正使之能够适应中等横向变速的介质,由此可以满足大多数探区的精度要求;相对叠后时间偏移来说,更适用于复杂构造,对目的层和储层的成像有较好的保幅性,所得结果能够更好地进行属性分析、A VO/A V A/A VP反演和其它参数反演。
3)实现叠前时间偏移的配套技术比较成熟和完善,如静校正和去噪等。
上述特点充分说明了我们应用叠前时间偏移技术的可行性、必要性和重要性。
下面就叠前时间偏移的基本情况、方法原理、方法技术、应用和与其它技术的比较以及应用该技术的可行性和必要性等做详细讨论和分析。
§2.1 概述
叠前时间偏移已进行了多年研究,上世纪九十年代初期开始初步应用,中后期在不少探区的地震勘探中发挥了重要作用,进入本世纪后开始了较为广泛的应用,目前部分处理公司和计算中心已把该技术作为常规软件加入到常规处理流程中,成为获取保幅信息实现属性分析、A VO/A V A/A VP反演和其它参数反演的重要步骤和依据。
自从上世纪九十年代以来,叠前时间偏移在国外取得了很大发展。
在理论研究方面,
Bleistein、Bortfeld和Hubral等进行了一系列有关真振幅叠前时间偏移理论的研究工作,Schneider给出了Kirchhoff型真振幅偏移权函数的一般公式。
Graham A. Winbow(1999)推出了控制振幅的三维叠前时间偏移的权函数的显式公式,并且利用真振幅权函数估计进行了振幅补偿。
另外,在权函数改进的基础上,提出了高精度的弯曲射线法绕射走时计算方法。
在应用方面,最近几年,在常规叠前时间偏移基础上,研究开发了多种保幅型叠前时间偏移软件,尤其是Kirchhoff保幅型叠前时间偏移软件取得了巨大成功。
它的特点是:计算效率高,目前的向量并行机和PC-Cluster机群使其计算效率又成倍提高,相对叠后时间偏移有较好的保幅性,更适合随后的属性分析、A VO/A V A/A VP反演和其它参数反演。
另外,保幅型的有限差分法和Fourier变换法叠前时间偏移软件也在应用中。
相信随着地震勘探技术的不断发展和计算机并行化的不断提高,Kirchhff型(甚至有限差分型、Fourier变换型和联合应用型)真振幅叠前时间偏移一定能应用于多波多分量地震资料及各向异性介质中。
国内在叠前时间偏移的研究上基本与国外保持同步,像吉林大学的孙建国博士在Kirchhoff型真振幅叠前时间偏移的理论研究方面取得了不少成果。
他定量表征了多种不同形式的真振幅权函数,研究了最佳有限偏移孔径的快速算法,分析了透射、衰减、薄互层、焦散和各向异性等的校正问题。
相对理论研究,国内的应用要落后于国外,应用成果还不是很多,尽管目前国内的许多处理公司和计算中心也都在做叠前时间偏移处理,但还没有真正把这种处理作为一种常规处理来看待。
随着计算机性能的不断提高和向量并行机(包括PC-Cluster机群)的广泛应用,加上叠前时间偏移相对叠后时间偏移较好的构造成像效果和保幅性以及属性分析、A VO/A V A/A VP反演和其它参数反演对地震资料保幅性的要求,叠前时间偏移处理在未来的几年一定会得到广泛应用,它解决构造和岩性问题的能力也会越来越强。
§2.2 方法原理
叠前时间偏移与叠后时间偏移和叠前深度偏移一样,都是基于三大数学工具,即Kirchhoff积分、有限差分和Fourier变换。
从原理和适用性上分析,叠后时间偏移是基于爆炸反射面的思想,做了水平层状介质的假设,不能对同一深度层具有不同叠加速度的不同倾斜层正确成像;而叠前时间偏移是基于绕射叠加或Claerbout的反射波成像原则,是一种成像射线成像(DMO是法向射线成像),能够解决叠后时间偏移存在的问题,适于V(z)介质和横向速度中等变化的V(x,y,z)介质,它对偏移速度场不是很敏感,具有较好的构造成像效果和保幅性,能满足大多数探区对地震资料的精度要求;至于叠前深度偏移,它也是基于绕射叠加或Claerbout的反射波成像原则,是一种没有横向偏差的曲射线成像,能适于横向速度急剧变化的V(x,y,z)介质,有好的构造成像效果和保幅性,但对偏移速度场非常敏感,成像效果直接依赖于偏移速度场的精度。
下面详细叙述有关叠前时间偏移的各种方法,其它方法从略。
一.Kirchhoff积分法叠前时间偏移
利用Kirchhff积分法作叠前时间偏移,一般在共炮点道集上进行,它比在其它道集上进行叠前偏移要方便些。
对二维和三维叠前偏移做法是一致的。
该方法的步骤是首先将共炮点记录从接收点上向地下外推。
外推时要先确定本道集可能产生反射波的地下空间范围,这个范围可以根据倾角、记录长度和道集的水平范围进行估算。
这个过程实际上是一个估算偏移孔径的反过程。
如果范围估计太大,一般会增加计算工作量,还会造成较多的偏移噪声背景。
如果把范围估计的太小,又会把反射界面丢失。
因此对向地下延拓的空间范围做一些模拟估算是必要的。
外推时使用一般Kirchhoff积分表达式:
⎰⎰∂+∂++-=A dxdy t
v R t y x u v R t y x u R v Rv t z y x u ]),0,,(),0,,([cos {21),,,(0000θπ (2-1) 2/122020])()/[(/cos z y y x x z R z +-+-==θ
式中R 为从地下(x,y,z)点到地面点(0,,000=z y x )的距离。
这样求出的结果,等于从地面某个炮点激发,在地下(x,y,z)点上接收的反射波记录。
在这个记录上有(x,y,z)点产生的反射波和z 深度以下的界面产生的反射波。
我们应当做的是把 (x,y,z) 点处的反射波放到该点上。
但是,在该点的记录还有很多其它深度点上的反射波。
因此,如何从这个点用积分公式延拓计算出地震道u(x,y,z,t),并从中取出用于在该点成像的波场值,这就是下一步的工作。
第二步,计算从炮点O 到地下R(x,z)点的地震波入射射线的走时d t 。
这可以用均方根速度rms v 去除炮点至地下R 点的距离近似求出。
或者用射线追踪法求取,就更准确。
用求出的下行波的走时d t 到u(x,y,z,t)的延拓记录的d t 时刻取出波场值做为该点的成像值。
第三步,将所有的深度点上的延拓波场都如第二步那样提取成像值,组成偏移剖面就完成了一个炮道集的Kirchhff 积分法偏移。
第四步,将所有的炮道集记录都做过上述三步处理后进行按地面点相重合的记录相叠加的原则进行叠加,即完成了叠前时间偏移。
二.有限差分法叠前时间偏移
在三维情况下,反射点轨迹变为一个旋转椭球面。
该椭球是绕炮检距方向由二维时的椭圆旋转而成。
如果取炮检距方向为x 方向,则椭球面的方程为:
222222224)/411(n n
t v z y x t v h =+++ (2-2) 通过波动方程的频散关系或波动方程的象征方程以及Fourier 变换,我们可以得到对应的三维波动方程
2222222222224)41(n
n t u v z u y u x u t v h ∂∂=∂∂+∂∂+∂∂+ (2-3) 如果炮检距方向与观测纵测线的方向成一定的α角度时需要进行坐标变换。
新坐标系下的方程为:
22222222224111n
t u v z u y x u C y u B x u A ∂∂=∂∂+∂∂∂+∂∂+∂∂ (2-4) 用有限差分法解(2-4)式有一定的难度,但它是可解的。
因此对三维面积观测的数据体用该方法进行叠前时间偏移在理论上是可实现的,目前尚未使用。
虽然各个方向的共炮检距道集也可以用(2-3)式进行偏移而且容易实现。
但是由于要在不同的α方向上抽取新的共炮检距道集,并要重新采样,同时剖面长度会长短不等,因此对处理效率会有影响。
(2-3)式虽然容易求解,但在炮检距方向有转角0≠α时,首先要将数据沿α方向和垂直α方向进行内插重排,这样内插重排后的三维数据体的水平切片将是某种菱形,造成纵横测线长短不一,
给处理带来不便。
如果仍按原坐标进行三维叠前偏移处理则必须用(2-4)式进行偏移。
三.Fourier 变换法叠前时间偏移
频率-波数(f-k )域叠前偏移是实现叠前时间偏移的一种有效方法。
Li(1991)用一组常速实现了叠前偏移。
用横向不变的速度偏移常炮检距数据可以在Fourier 域进行,与Kirchhoff 偏移相比,它具有成像速度快,能处理陡倾角且不会产生算子假频(是一宽带算子)的特点。
另外,该算子考虑了由于通过层状介质而发生折射弯曲所造成的相位和振幅变化。
另外,F-K 偏移算子可以分解为NMO+DMO+ZOM ,在常速偏移下,分解正确。
若速度随深度变化,这种分解对NMO+DMO 部分只是近似值。
二维情况下,F-K 域叠前时间偏移的向下延拓波场为:
),,,(ˆ),,,(ˆ)exp(),,,(ωωωz k k S z k k D h ik m ik dk dk z h m P h
m h m h m h m +=⎰⎰ (2-5) 对层状v(z)介质,传播算子),,,(ωz k k S h m 由下式给出:
⎥⎦
⎤⎢⎣⎡Φ=∑j h m j h m z k k i z k k S ),,,(exp ),,,(ωω (2-6)
其中, []
22
22
)()(2),,,(h m j h m j j h m j k k k k k k z z k k +-+--=Φω (2-7a ) ∑=j
j z z (2-7b )
j
j v k ω2= (2-7c ) v j 是层速度。
(2-7)式是常速频散关系的一扩展形式。
§2.3 方法技术
在上述有关叠前时间偏移的三大方法中,为保证储层成像的准确性和储层特征的保真性,达到为随后的各种属性分析、A VO/A V A/A VP 反演和其它参数反演提供高质量地震数据的目的,目前常用的方法技术包括:Kirchhoff 积分法中的真振幅权函数估计技术和Fourier 变换法中的稳相技术。
下面分别加以叙述。
一.Kirchhoff 积分法中的真振幅权函数估计技术
要正确进行属性分析和A VO/A V A/A VP 反演,必须利用真振幅地震资料。
真振幅的基本思想是在实现构造成像的同时给出反射系数和震源子波信息。
从80年代开始,Bortfeld 和Hubral 等进行了一系列的真振幅叠前时间偏移理论的研究工作,Schneider (1993)给出了Kirchhoff 型真振幅偏移权函数的一般公式,Graham A. Winbow (1999)推出了三维保幅型叠前时间偏移的权函数的显式公式,并且利用真振幅权函数估计进行了振幅补偿。
Schneider (1993)给出的Kirchhoff 积分型权函数的一般表达式为
2/112/1'11det det )det()cos (cos ),(GR
SR GR SR N N N N v x w +=θθξ (2-8) 其中,ξ 是炮-检中点,x 是成像点,1θ与'1θ分别是震源和接收点射线在地面的出射角,1
v
是地表速度,SR N 与GR N 是Schneider 定义的矩阵。
对上式进行分析、整理和简化,可以得到适用于不同介质的权函数。
Graham A. Winbow (1999)推出的三维保幅型叠前时间偏移的权函数的显式公式为:对常速介质
'
111'1212'11cos cos )
cos )(cos cos (cos θθθθθθv w += (2-9a )
对速度随深度线性变化的介质
⎥⎦
⎤⎢⎣⎡-+++=-αλδθθ22'11cos 2sin cos )1(cos )1(1
11X X v w (2-9b ) 对于水平层或近似水平层反射,简单的权函数就很好了。
而对于陡倾角地层的反射若不使用保幅型权函数则A VO 的误差就很大。
另外,走时计算可以使用射线追踪或求解程函方程。
在复杂的V(z)和VTI 介质中,为提高走时计算的精度,采用了弯曲射线法。
二.Fourier 变换法中的稳相技术
为完成常偏移距-常方位角数据的叠前时间偏移成像,Etgen (1998)和Rietveld 等人(1998)设计了一F-K 偏移算子。
它不需要使用NMO 和DMO ,对v(z)介质有正确的运动学特征。
在Etgen (1998)工作的基础上,Dai 等人(1999)提供了两种替代方法,他们相对Etgen (1998)设计的算子具有一定的优点。
这两种方法是:
(1)在炮检距-波数域中的稳相法,其偏移算子为
)],,,,(exp[),,,(ωψωz h k k i A z h k S h m m = (2-10)
(2)从空间算子开始,得到在中点域中的稳相的一显式解,其偏移算子为
(2-11) 该算子由空间上的相移和时间上的相移组成。
其相位中的第一项把能量沿炮检线移动一个DMO 距离;第二项代表DMO 校正;第三项是零炮检距相移偏移。
上述两种算子给出的总相位完全一致。
尽管第二种方法是一常速解,但由于它的高计算效率,并通过使用一RMS 模型或时间偏移速度模型使它能应用于v(z)函数和v(x,z)函数。
三.叠前时间偏移的保真度分析
1.Kirchhoff 型叠前时间偏移
它可以直接计算偏移响应。
因为该方法是非递归的,可以用如下的方程来计算地下不同成像点i i y x ,的走时,0t 采用4阶Taner Koehler 展开:
2/16443022
2
))((r c r c t v r t t rms ∆⋅+∆⋅+∆+= (2-12) 其中,2
/12244))()((i i y y x x r -+-=∆。
(2-12)式中的rms v 可以代表成像点的速度,也
可以表示成像点与地表中点的速度,前者更易于实现,而后者精度更高。
用t δ表示射线追踪绕射响应与实际偏移响应间的偏差,进而可用相对速度误差v v /δ来表示这种偏差
3228rms tv r v t ∆-=δδ, 224r t tv v v rms ∆-=δδ (2-13)
其中衡量保真度的速度场为空间坐标x,y,z 的线性梯度函数:
y x z z y x v 07.007.06.06000),,(-++= (2-14)
模型的空间范围为:x 方向15km ,y 方向7.5km 。
图2-1是用v v
δ测定的Kirchhoff 型叠前时
间偏移的成像保真度图,描述的是3000m 深度处成像点的保真度。
图2-1 Kirchhoff 型叠前时间偏移算法的成像保真度 图2-2 扩展的Stolt 法叠前时间偏移的成像保真度
2.F-K 域叠前时间偏移
为了让F-K 偏移适应速度的横向变化,常采用扩展的Stolt 偏移方法,这种扩展常通过非线性优化来实现,所使用的优化函数对第i 个响应表示如下:
2/16443022
2
))((),(),,(r c r c t v r t y x y x f rms i i ∆⋅+∆⋅+∆+=+ττ (2-15) 其中,),,(i y x f τ为优化扩展函数,),(y x i τ为第i 个响应的走时,对方程右边采用4阶Taner
Koehler 展开可用来近似偏移响应。
扩展的Stolt 偏移的成像保真度的分析方法同Kirchhoff 型偏移中的分析方法,速度模型也相同,所描述的3000m 深度处的成像保真度如图2-2所示,从中可以看出该方法的保真度要好于图2-1中所示的情况。
因为扩展的Stolt 偏移比Kirchhoff 方法更加有效,因此在该速度模型下可用扩展的Stolt 偏移法代替Kirchhoff 方法,但对其它的速度模型不能以此推广。
§2.4 应用与比较
基于上述方法原理和方法技术,通过国内外的一些实例,说明各种方法的应用效果,并与叠后时间偏移和叠前深度偏移进行对比分析。
一.应用实例
应用实例1
图2-3、2-4和2-5给出了一三维数据体的叠后叠前时间偏移结果。
图2-3是纵测线37和55的叠后时间偏移剖面,图2-4是基于一开发软件得到的相同纵测线的叠前时间偏移(PSTM )剖面,图2-5是通过另一开发软件得到的相同纵测线的叠前时间偏移剖面。
对比图2-3和图2-4,可以看出叠前偏移效果得到了改善(见同相轴2和3),但图2-4中的一些浅层同相轴(如同相轴1)的成像效果不如叠后偏移好;对比图2-4和图2-5,可以很明显地看到图2-5的成像效果不如图2-4好,特别是图2-5的浅层同相轴没有正确成像,例如同
相轴1下面存在一个很强的假同相轴。
可见对这种软件,浅层成像对成像控制相当敏感。
因此开发一套精确有效且实用的叠前时间偏移软件仍是一件极具挑战性的工作。
图2-3 叠后时间偏移剖面(a) 纵测线37;(b) 纵测线55
图2-4 基于一开发软件得到的叠前时间偏移剖面(a) 纵测线37;(b) 纵测线55
图2-5 通过另一开发软件得到的叠前时间偏移剖面(a) 纵测线37;(b) 纵测线55
图2-6给出了来自一三维数据体一纵测线时间偏移剖面的部分结果。
2-6a图是叠前时间偏移剖面;2-6b图是叠后时间偏移剖面。
比较两图可以看到叠前偏移的成像效果明显好于叠后偏移。
图2-7是共成像道集(CIG),结果表明了所用速度的正确性。
叠前时间偏移改善成像质量是不难理解的,但如果在叠加中失去了一信号,不管偏移方法如何准确,在叠后偏移中不可能恢复该信号。
应用实例3
图2-8是利用直射线PSTM(左边)和优化的6阶走时方程PSTM(右边)得到的共成像道集(CIG)。
数据体来源于墨西哥湾探区。
结果表明曲射线的效果明显好于直射线。
图2-9a和2-9b分别是利用直射线和曲射线叠前时间偏移得到的成像剖面。
观测中所用的最大偏移距是8km。
图2-9b显示出的断层要比图2-9a清晰得多。
在这种情况下,我们是把偏移RMS速度单独来调整以适合直射线和曲射线偏移算法。
图2-10给出了利用曲射线PSTM(左边的两道集)和曲射线+VTI校正(右边的两道集)得到的共成像道集(CIG)。
结果表明若没有VTI校正,仅靠调整速度很难使反射层同相轴拉平(左边的两道集)。
考虑了VTI效应后,CIG上的同相轴基本在同一水平线上。
图2-11a给出了没有考虑各向异性效应的PSTM 剖面的一部分;图2-11b给出了考虑各向异性效应的PSTM剖面的相同部分(曲射线+VTI 校正)。
注意到在两剖面的下部有明显的差异。
图2-6 (a) 纵测线35的叠前时间偏移剖面
图2-6 (b) 纵测线35的叠后时间偏移剖面
图2-12a 和图2-12b 分别给出了稳相法常炮检距-常方位角相移叠前时间偏移和常规Kirchhoff 积分法叠前时间偏移的成像结果。
二者均基于同一横向变RMS 速度的偏移速度场。
比较两图可以看到二者的成像效果基本相同,在这种情况下,可以充分利用相移法的高计算效率来解决3D 叠前时间偏移问题。
应用实例5
对叠前时间偏移后的地震资料进行多角度和多尺度的A V A 反演可以求取弹性和岩石物性参数,建立多尺度参数与地质相类型的关系,进一步估算储层参数。
图2-13给出了一构造成像图和一A V A 反射系数剖面的联合显示形式。
其中,d 是用于A V A 反演的数据向量。
图2-14是估算出的P 波速度反差剖面p p c c / 。
图2-15是对叠前偏移数据进行A V A 反演的
结果:(a) 没有做薄层校正;(b) 做
了薄层校正。
其中,兰色和绿色曲
线表示反演结果;黑色曲线是声波
曲线经滤波后的结果。
图2-16是
一构造成像图、一A V A 反射系数
剖面及其多尺度表示的联合显示
图。
其中p-σ面用于多尺度A V A
反演。
二.与其它技术的比较
通过上面的理论分析和应用
实例,对叠前时间偏移和它与叠后
时间偏移和叠前深度偏移的对比
有如下认识与结论:
1)水平层状介质条件下,叠
后时间偏移、叠前时间偏移、叠前
深度偏移都可以很好的成像。
2)横向速度变化不大的条件
下,叠前时间偏移与叠前深度偏移
成像效果相差不大,但对于陡倾角
反射层,叠前时间偏移的归位效果
优于常规的叠后时间偏移。
叠后时
间偏移对同一深度层具有不同叠
加速度的不同倾斜层不能正确成
像。
3)横向速度变化剧烈的情况
下,常规CMP 道集内的反射同相
轴具有复杂的非双曲线动态时差,
常规的叠后时间偏移和叠前时间
偏移都很难正确成像,而叠前深度
偏移却可以较好的成像。
图
图2-7 在CDP2895(纵测线36,横测线60)和CDP4313
(纵测线54,横测线20)处的叠前时间偏移的共成像道集。
图2-8 利用直射线(左边两图)和曲射线(右边图2-10 利用曲射线(左边两图)和曲射线+VTI 两图)算法得到的共成像道集校正(右边两图)算法得到的共成像道集
(a) 直射线算法(b) 曲射线算法
图2-9 叠前时间偏移剖面
(a) (b)
图2-11 叠前时间偏移剖面的一部分
(a) 没有考虑各向异性效应(仅用曲射线算法)(b) 考虑了各向异性效应(曲射线+VTI校正)
(a) (b)
图2-12 叠前时间偏移结果
(a) 基于稳相原理的常炮检距-常方位角相移法 (b) 一般Kirchhoff 积分法
图2-13 一构造成像图和一A V A 反射系数剖面的联合显示形式。
其中,d 是用于A V A 反演的数据向量
图2-14 估算出的P 波速度反差剖面p p c c /
图2-15 对叠前偏移数据进行A V A反演的结果:(a) 没有做薄层校正;(b) 做了薄层校正。
其中,兰色和绿色曲线表示反演结果;黑色曲线是声波速度曲线经滤波后的结果
图2-16 一构造成像图、一A V A反射系数剖面及其多尺度表示的联合显示图。
其中,p-σ面用于多尺度A V A反演
4)叠前深度偏移无疑是解决成像问题最好的方法,但由于依赖速度的准确性和计算量大的问题,使之很难大面积应用于生产,而叠前时间偏移运算速度快,求取速度较准确,它们之间联合应用,减少了深度偏移迭代的次数,提高了效率。
运用偏移的方法求取速度,不需要精确的解释,不需要层层剥皮,能很快准确地对新区建立速度模型,偏移的联合应用,势必为深度偏移的应用开辟一条新路。
因此,从某种意义上讲,叠前时间偏移是准确求取偏移速度场做好叠前深度偏移的一种主要方法。
近年来,向量并行巨型机和微机集群(尤其是微机集群)在地震叠前成像中的广泛应用也为提高二者联合应用的计算效率提供了硬件条件。
5)叠前时间偏移与叠前深度偏移相比,有明显高的计算效率;改进/修正后的保幅型叠前时间偏移能满足大多数探区对偏移成像精度的要求。
6)叠前时间偏移和叠前深度偏移在保幅性上明显好于叠后时间偏移。
两种叠前偏移资料更适合进行属性分析、A VO/A V A/A VP反演和其它参数反演。
§2.5 可行性与必要性
叠前时间偏移是在叠前进行,所用的各种道集(如共炮集、共炮检距道集等)均可独立偏移,这为并行处理提供了条件;另外,计算机技术发展很快,而软硬件成本和价格却在不断降低,再加上叠前时间偏移技术的不断改进、完善和成功的应用以及对储层成像的准确性和储层特征的保真性。
这些都说明了在地震成像中应用叠前时间偏移技术的可行性和必要性。
一.可行性分析
1.叠前叠后联合偏移能降低成本提高叠前时间偏移的成像精度
近年来,随着勘探目标的日益复杂化,地震资料处理的难度越来越大,对地震成像的精度要求也越来越高。
常规的叠后时间偏移方法由于受到其理论本身的制约,不能够解决这类复杂构造的成像问题,而叠前偏移技术由于其成像精度高,越来越受到人们的重视。
因为在叠加剖面上比较容易知道处理区的基本地质情况,建模相对精确,因此一般可以利用叠后深度偏移建立起初始模型,由此作为叠前偏移的输入模型。
另外,叠前偏移能够进行并行处理并具有以下四个适应性条件。
1)理论波场的适应性条件
从理论上讲,叠前深度偏移的波场理论更适应复杂地质体的情况,但是叠前深度偏移的成像效果必须依赖于准确的速度-深度模型,而模型的迭代和修改是一个十分复杂和费时的过程。
因此,叠前深度偏移的周期长,花费也是相当昂贵的,至今尚未大面积用于生产实践。
所以,对于横向速度变化并不是十分剧烈的地质构造,叠前时间偏移是一种经济有效的方法。
叠前时间偏移方法并不需要解释层位和建立速度-深度模型。
通过常规的速度分析及偏移后的速度修正,经过几次迭代就可以得到较准确的偏移速度场。
所以比较起来,叠前时间偏移的周期短,效果好,成本也是比较低的。
叠前时间偏移得到的叠前道集,可以消除倾角因素的影响和由于速度横向变化所引起的叠加成像问题,从而为储层研究、叠前A VO/A V A/A VP分析、裂缝分析及反演提供准确的依据,尤其是在A VO处理中的应用,成功解决了A VO资料的归位问题,特别是道集及属性剖面的归位,使在A VO道集及属性剖面上直接确定A VO异常的分布范围和位置成为可能,同时使绕射能量归位变好,减弱了随机噪音的干扰,使A VO响应更清晰,提高了A VO资料的分析质量;另外该叠前道集又能为深度偏移提供准确的层位及速度模型,减少深度偏移的迭代次数和处理时间,极大地提高迭代效率和成像精度。
从而降低叠前深度偏移的费用,提高其生产可行性。
2)硬件环境的适应性条件。