地震偏移成像基本原理.ppt
《地震偏移成像基本原理.ppt》由会员分享,可在线阅读,更多相关《地震偏移成像基本原理.ppt(114页珍藏版)》请在三一文库上搜索。
1、地震成像原理与方法,第一章 偏移成像 第二章 三维叠前深度偏移成像理论与方法 第三章 共聚集点偏移 第四章 共反射面叠加 第五章 偏移速度建模,第一章 偏移成像,1.1 偏移成像的基本原理 1.2 波动方程偏移 1.3 叠前偏移 1.4 偏移速度分析 1.5 深度偏移 1.6 三维偏移 1.7 二维和三维叠前深度偏移,地震技术的发展趋势:,三维叠前深度偏移(3DPSDM)-地震成像 (波动方程法3DPSDM, CRS叠加, CFP偏移) 四维地震-开发地震 (VSP技术, P-S技术, 井间地震, 3D_AVO技术, 4D地震, 弹性波阻抗反演, 裂缝分析, 岩石物理, 地震相与地震属性分析,
2、 油藏描述等),Reflection point smearing,NMO_DMO_PostMig Correction,三大处理技术: 反褶积、叠加、和偏移成像 反褶积和叠加引自其它相关学科 偏移成像基于古典技术 偏移成像: 1.具有地震勘探本身的特征。 2.计算机使其研究由地震波运动学特征 过渡到地震波动力学特征 3.提高地震空间分辨率和保真度,一偏移成像的概念,偏移,反偏移,反射地震方法: 1.激发弹性波,2.记录反射波, 3.研究地质岩层结构和物性特征。是一种反散射问题。 反射地震成像分做两步: 1. 记录反射波,2. 处理反射波。地震偏移技术是使反射界面最佳成像的一种技术。,1.1
3、偏移成像的基本原理,地震偏移: 叠前或/和叠后偏移 叠前偏移: 使CSP道集记录或COF道集记录中的反射波归位, 绕射波收敛 叠后偏移: 基于水平叠加剖面,采用爆炸反射面概念实现倾斜反射层归位和绕射波收敛 偏移原理和偏移效果见下图,1. 偏移成像的基本概念,偏移原理图,偏移过程的定量分析图,2. 发展史,1).古典的偏移技术(60年代前) -反射点的空间位置成像; 2).早期的计算机偏移技术(6070年代)-定性和概念性地对反射波运动学特征成像; 3).波动方程偏移技术(70年代后)-定性或/和定量地对反射波运动学或/和动力学特征成像.,1).有限差分法波动方程偏移:70年代初期,J.Clae
4、rbout教授首先提出了用有限差分法解单程波动方程的近似式,用地面观测的地震数据重建地震波在地下传播过程中的波场,从这些传播过程的波场中提取使地震界面成像的那些数据,组成地震偏移剖面。由于这种偏移方法在计算过程中要解波动方程或其近似式,所以被称为波动方程法偏移技术。 2). Kirchhoff积分法波动方程法偏移:70年代中期,French和Schneider等在绕射偏移法的基础上使用了波动方程解的Kirchhoff积分公式,发展为地震偏移的波动方程积分法。使绕射偏移建立在可靠的波的基本原理上。因而改善了偏移剖面,取得了良好的效果。 3).富里叶变换法波动方程法偏移:70年代后期,Stolt和
5、Gazdag等又先后提出了在频率-波数域解波动方程,外推地震波场的方法。这种方法被称为F-K域偏移方法。由于该方法计算简单,效率高,因而很快得到了推广。,波动方程偏移技术的发展,3. 偏移方法分类,二基于射线理论的叠后偏移与叠前偏移,经典的偏移方法和早期的计算机偏移方法 都是基于射线理论 经典的偏移方法只研究到达时间。叠后偏移有圆弧切线法和线段移动法;叠前偏移包括椭圆切线法和交会法等 早期的计算机偏移方法利用了波前、绕射等地震波传播的惠更斯原理,尽管只是定性的、概念性的,但与手工操作法相比偏移剖面除了归位精度提高外,还考虑了波形特征。叠后偏移有波前模糊法、绕射曲线叠加法; 叠前偏移有Rockw
6、ell偏移叠加法和Paturet-Tariel偏移叠加法等。,1. 叠后偏移,叠后偏移: 即叠加偏移,是对叠加后的地震记录做偏移。下面介绍圆弧切线法、波前模糊法和绕射曲线(面)叠加法。 1). 圆弧切线法 一次反射波NMO后, 得到时间叠加剖面 由 (1.1.1) 得到视深度剖面 如果界面的倾角 =0或者很小,例如只有 或更小,则视深度界面就是真深度界面。如果界面倾角不可忽略,则应当进行倾角校正,以求出反射界面的真实位置。,校正的做法是以地面各点为圆心,以各点下至视界面的垂直距离为半径做圆弧,其圆弧族的切线即为校正后的反射界面(v=cont)。,当速度是深度的函数时,例如 为常数时,则圆弧的圆
7、心不位于地面上,而位于地面点的正下方某深度上。这时,圆心的深度和圆弧的半径由下式求出:,(1.1.2),2). 波前模糊法,波前模糊法也可以称为波前切线法,它是对叠加后的地震剖面进行偏移的方法。这个方法是反推反射界面上的波场。,以地面接收点为中心,把相当于反射到达时间上的值送到以 的深度为半径的圆弧上去。如果我们把深度z仍以双程时间表示,就把反射数值送到以t为半径的圆弧上去(图1-4)。把各道上的所有反射波值都按这个原则去做,并把送到同一点的值叠加起来,就可以组成偏移剖面。把某道 上某时间t上的振幅值送到相邻各道上的时间 由下式算出:,(1.1.3),其中,用波前振幅叠加来求反射界面发出的波前
8、实际上就是用这种方法做切线。,要求:较密的地震道和较高的信噪比,以得到满意的偏移剖面。,3). 绕射曲线(面)叠加法,绕射曲线或绕射曲面叠加法是把地震剖面上的波场振幅值按绕射波时距曲线进行相加。因为绕射波时距曲线与所有反射波的时距曲线形状相比较,其凸率最大,故亦可称它为最大凸率法。,具体做法是,当要得到地震剖面上某个 点的偏移后的数据时,我们要计算一条以这点为顶点的绕射双曲线。它在各道上的时间t由下式算出:,(1.1.4),式中,在进行偏移时我们把各道上等于上式时间t的波场值取出来叠加在 点的波场值上,这就算完成了 点的偏移处理,如图1-5所示。,无论是波前模糊法还是 绕射叠加法,其基本原理都
9、 是根据惠更斯原理提出来的。 波前模糊法是把一个道 上的波场值送到各个道上去 叠加输出道法;而绕射叠加法是把各道上的相应值取来在一道上叠加输入道法。两者都符合反射波归 位和绕射波收敛的要求,而且它们的叠加值也相等。,波前弧或绕射曲线在x方向上的范围L称为偏移孔径。L的范围是由最大实际倾角来决定的。倾角越大,L越大;有效波越深,L也越大。L的大小可用下式来估算:,(1.1.5),孔径的中心,原则上 应当位于 处,但也 可以是不对称的。,图1-7是用绕射叠加偏移法处理前后的地震剖面。从对比中可以看出,偏移后剖面上的地层层位关系得到了正确的反映。有利于地质解释。,2叠前偏移 叠前偏移: 即偏移叠加,
10、是对叠加前的多次覆盖的地震记录先偏移,再叠加。下面介绍椭圆切线法、Rockwell偏移叠加法和Paturet-Tariel偏移叠加法。,1)椭圆切线法 当给定CSP记录时,可用椭圆切线法(图1-8)。反射点(2D)位于以炮点和接收点为焦点的椭圆上,这个椭圆的方程可表示为:,(1.1.6),对每个炮检距的记录上的反射波画好椭圆弧。做椭圆弧族的切线即为偏移后的剖面。,2)Rockwell偏移叠加法 Rockwell偏移叠加法实际上是叠后偏移所使用的波前模糊法的一个扩展。,具体做法:把每个记录道上任一t时刻的采样值,在以炮检距中点的地面点为原点的直角坐标系中送到以 为长轴, 为短轴的椭圆与各个地震记
11、录道垂直线相交的各个点上去,并且与其它地震道送至该交点上的采样振幅值相加,即得偏移叠加剖面。 偏移叠加实质上是用振幅叠加来做切线的。,3)Paturet-Tariel偏移叠加法 1971年Paturet-Tariel用相同炮检距的剖面进行叠前偏移,把所有相同炮检距的偏移后的剖面叠加得到偏移叠加剖面。叠前偏移的原理如图1-9所示。,绕射点M所产生的绕射波到达时曲线为:,(1.1.7),当炮检距 时,上式表现为:,(1.1.8),式中 为从M点到A点的双倍旅行时间。 和 的曲线表示在图1-9的右图中。 为了进行偏移,我们应当把 的曲线上的地震能量(即采样点振幅)送到零炮检距绕射双曲线的顶点M上去叠
12、加。这样, 把各个相同炮检距的剖面偏移后叠加在一起即得偏移叠加剖面。,图1-10. 偏移叠加剖面与叠后偏移剖面对比图 (a). 水平叠加剖面,(b). 叠后偏移剖面; (c). 偏移叠加剖面,三基于波动方程的波场外推与地震成像原理 使用波动方程进行偏移,首先就是要重建反射波的原来波场。反射界面上刚刚产生的反射波,就认为是该反射面的像。 为进行波场外推,把波动方程分解为上行波方程和下行波方程。 1上行波和下行波,波动方程有两个解, 一般表示为 。在地震勘探中一般取深度方向向下为正z的方向。向正z方向传播的地震波称为下行波,即用 所表示的波。向负z方向传播的波为上行波,即用 代表的波。下行波即入射
13、波,上行波为反射波。,只有在均匀各向同性完全弹性介质的情况下上行波和下行波才是分离的。分离过程如下: 二维波动方程为:,(1.1.9),对(1.1.9)式相对x和t做二维付里叶正变换,并进行算子分解得到:,(1.1.10),其中利用了波散关系:,(1.1.11),由(1.1.10)式得出:,(1.1.12),其中,正号代表上行波方程,负号代表下行波方程。,2波场外推 正向外推就是根据波在当前位置上的振动情况向波的自然传播方向用计算手段预测出波场。反向外推是向波的自然传播方向的反方向上重建原来的波场。对一个波场应是进行正向外推还是反向外推均有物理问题决定。,1)上行波的外推,(1.1.13),积
14、分结果为:,(1.1.14),由此得出上行波的正、反向外推式。,(1)上行波正向外推公式 上行波的正向外推式就是向负z方向的外推公式。从(1.1.14)式可求出为:,(1.1.15),根据这个公式可以计算模拟反射波的地震记录(地震图)。 (2)上行波反向外推公式 上行波的反向外推式就是向正z方向的外推公式。从(1.1.14)式可得出为:,(1.1.16),根据这个公式可以进行地震记录的向下半空间延拓,求出地下任何一点的波场,实现地震波偏移的目的。,2)下行波的外推,(1.1.17),积分结果为:,(1.1.18),据此可以得出下行波的正、反向外推公式。,(1)下行波正向外推公式 下行波的正向外
15、推式是指沿正z方向的外推。其外推式为:,(1.1.19),这个方程可用来模拟下行波的地震记录。 (2)下行波反向外推公式 下行波的反向外推是指沿负z方向的外推。其外推式为:,(1.1.20),上式可用来从下行波场进行反向求源的计算工作。,下面分析波场本身的条件对外推结果的影响,(1.1.21),当 时, 为正或负的实数,这时所有外推公式中存在虚指数。说明在外推过程中波场发生相位变化。一般都能得出正确的结果。,当 时, 值为虚数:,(1.1.22),(1.1.23),波场外推时只有振幅变化,而无相位变化。当指数项取负号时,外推的波场迅速衰减,称这种波为倏逝波。当指数项取正号时,外推波场迅速增大,
16、这是一种实际不存在的波,只是进行波场计算时发生,我们称它为耗损波。在计算中要避免发生这种情况。,当 时,上行波的外推式可写为:,(1.1.24),此时反向外推遇到倏逝波,正向外推发生耗损波。分别表示为:,(1.1.25),(1.1.26),由此可见,用上行波方程进行向下波场外推永远是计算稳定的。而用上行波方程进行正向外推就可能遇到耗损波,因此有可能是不稳定的。除非在计算中不断地把 的波场滤除掉。 同理可求出 时下行波的外推式为:,(1.1.27),此时也是反向外推遇到倏逝波,正向外推遇到耗损波。,3)波场外推的Kirchhoff积分法,Kirchhoff积分法并不直接解波动方程,而是用数学方法
17、来描述关于波的传播的惠更斯原理,从而求出空间上任一点波场值的。,Kirchhoff早在1883年就证明了,从扰动区向外某点 传播的波的t时刻的波场 可以从扰动区封闭表面上的 波场 以及该波场对时间和表面法线方向的导数通 过积分式求出来。因此要假定 在封闭面上和封闭面内有直至二阶导数的连续性。 Kirchhoff利用了格林定理:,(1.1.28),u取为波场函数,,当把观测点用包含有波前面在内的封闭曲面包围起来,如图1-11(a),(b)那样的封闭时,这样的封闭面S和它所包围的体积V作为(1.1.28)式的积分限,经过一定的推导后得出 点的正向外推波场为:,(1.1.29),这里 的方向取封闭表
18、面的外法线方向。 如果把观测点M移至封闭面外,则有:,(1.1.30),( 1.1.29)式中,(1.1.29)式就是著名的Kirchhoff积分。它描述了物理波场传播的过程,也满足奇次波动方程,是它的积分形式解。对我们来说,也可以称它为正向外推公式。 注意: Kirchhoff积分只满足均匀介质的情况。,是推迟场,下面讨论用Kirchhoff积分进行波场反向外推问题(地震偏移)。这时,所取的封闭体积V应在波前传播方向的反方向,计算点 就在这个封闭体内。根据格林定理同样可求出形式上相同的反向外推的Kirchhoff积分式:,(1.1.31),式中的u不再是推迟场,而是超前场 。,(1.1.31
19、)式为用于波场反向外推的Kirchhoff积分式。它可用于上行波的反向外推,也可用于下行波的反向外推。当然,这种外推与正向外推不同,它不代表一个物理过程,而只是一种重建波场的计算过程。,3地震反射波场成像,从波动场的观点叙述反射波成像的一般原理。 地震成像-地震偏移 反射系数值-反映该反射点反射系数相对值的反射波振幅 反射成像实际上就是把地面上观测到的反射波归位到产生它的反射点上去。地震偏移与地震成像在现阶段可以视为同一概念。,地震偏移成像: 一是上行波场的反向外推; 二是在外推波场中提取成像值。 Claerbout提出下述反射波成像原则: 反射面位于这些点上,其入射波 的初至与反射波的产生时
20、间相同。 如图1-12所示,反射波成像的基本公式可写为:,(1.1.32),(1.1.32)式没有考虑反射系数随着入射角变化的情况,它实质上是相位信息的公式。或者说,它对接近法线入射的情况时基本是正确的,能够反映反射系数在各点上的变化情况。 应用(1.1.32)式涉及到要选择下行波的初始时间。这是一个困难问题。我们通过假设下行波是最小相位而避开这个问题。我们把 作为初始时间,可推出如下的反射图象公式:,(1.1.38),当下行波是脉冲波时, (1.1.38)式很精确。但是,如果 是 一个短延续长度的子波时, 它只是一个很好的近似成像公式。,爆炸反射面的概念:,水平叠加剖面后的自激自收剖面等价于
21、在反射界面上同时爆炸产生地震波,并以半速度向外传播,在地面上观测到的上行波剖面。这就是Loewenthal等人首先 提出来的爆炸反射面的概念(图1-14)。 这个概念对于理解水平叠加剖面的偏移 成像是很重要的。因为它比较直观地说 明了这种剖面的成像原理。比前面所述 的反射波成像的一般原理要容易理解得 多。在这里我们引入了两种等价。一个 是水平叠加剖面等价于自激自收剖面。 另一个是自激自收剖面等价于界面上同 时激发在地面记录的上行波剖面。这种 等价只是概念上的,实际上只有一种水 平叠加剖面,并没有三种剖面。,1.2 波动方程偏移,地震偏移成像技术发展到今天已经产生了各种形式的在各种域实现的方法。
22、历史上曾经起过作用的根据几何光学原理的成像方法已经被淘汰。现在正在流行的是建立在波动方程基础上的三种方法,即Kirchhoff积分法,有限差分法和F-K法及其各种变形。这三种方法由于有相同的数理基础,因此它们的原理相同。同时,因计算方法不同,它们之间又有许多不同之处。下面讨论三种方法对水平叠加地震剖面的偏移。,一频率-波数域波动方程偏移,采用爆炸反射面的理论。为了成像,要求向地面以下反向外推地震波场。假定z轴垂直向下为正,测线沿x轴,则u(x,z,0)表示偏移后的真实剖面,而u(x,0,t)是未偏移的叠加剖面。,在均匀各向同性完全弹性介质中,用半速度代替地震波传播速度,则标量波动方程变为:,(
23、1.2.1),(1.2.2),对(1.2.1)式进行傅里叶变换并利用(1.2.2)式有,(1.2.3),其中正号代表上行波,负号是下行波。,1Stolt偏移法,设 为 的二维傅里叶变换,对(1.2.1) 式进行上述变换得到:,将(1.2.3)式代入上式有,按上行波求解,即取正值得,其中A与t无关。令t=0,上式变为:,从而, 是待求的偏移剖面 的傅里叶变换。,下面讨论用水平叠加剖面 如何求出 。对 做傅里叶逆变换得:,令z=0,上式变为:,(1.2.4),设水平叠加剖面 的二维傅里叶变换为 ,则,(1.2.5),其逆变换为:,(1.2.6),比较(1.2.4)与(1.2.6)有,这样,按上行波
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 地震 偏移 成像 基本原理
链接地址:https://www.31doc.com/p-2978053.html