最新地震处理教程——第七章 倾斜叠加及其应用
- 格式:doc
- 大小:2.42 MB
- 文档页数:24
地震作用最大方向附加斜交构件的问题STAWE算出的最大地震作用方向是根据地震作用和坐标角度进行刚度矩阵的算法得到一个函数的。
这个函数我查阅过相关资料碍于自己的数学水平有限说句实话实在是看不懂。
详见附件的一篇论文《地震最大作用方向的确定》不过在实际设计工作中,经过自己的实践我发现最大地震作用角度和自己刚度布置有关。
如果是结构本身是比较规则的情况下。
如果刚度布置分部不是按建筑形状的分部布置的。
比较极端的情况如下图,刚度布置是由45度对称布置的。
左下角和右上角是布置的剪力墙截图可能看不得不是特别清楚。
振型号周期转角平动系数 (X+Y)扭转系数1 1.1296 45.00 1.00 ( 0.50+0.50 ) 0.002 0.7428 135.00 1.00 ( 0.50+0.50 ) 0.003 0.490 10.00 0.00 ( 0.00+0.00 ) 1.00地震作用最大的方向 =45.001 (度)由以上一个简单的例子可以发现,地震作用最大方向是通常是由刚度对称分部的主轴方向来的。
通过看X向和Y向转角可以判断这个结构的主轴是45度的坐标系了,而不是我们默认90度的。
(当然会有例外的情况,这个还和体量大小,刚心质心多方面因素相关。
只不过这个主轴影响是最大的,如果有例外的情况不符合这个经验,在从其他地方着手,不要和我抬杠)。
我只要改变一下结构布置,这个结构布置按建筑形状来布置刚度左下角右下角布置了剪力墙。
地震最不利方向就会正常。
主轴方向也正常了变成0 度和90度。
考虑扭转耦联时的振动周期(秒)、X,Y 方向的平动系数、扭转系数振型号周期转角平动系数 (X+Y) 扭转系数1 1.3212 0.00 0.79 ( 0.79+0.00 ) 0.212 0.8733 90.00 1.00 ( 0.00+1.00 ) 0.003 0.4374 180.00 0.21 ( 0.21+0.00 ) 0.79地震作用最大的方向 = 0.000 (度)单个构件的刚度没有变,只是改变了布置的位置就会改变最大地震作用角度。
§7.2地震映像及其应用地震映像(又称高密度地震勘探和地震多波勘探),是基于反射波法中的最佳偏移距技术发展起来的。
这种方法可以利用多种波作为有效波来进行探测,也可以根据探测目的要求仅采用一种特定的波作为有效波。
除常见的折射波、反射波、绕射波外,还可以利用有一定规律的面波、横波和转换波。
在这种方法中,每一个地震记录都采用相同的偏移距激发和接收。
在该偏移距处接收到的有效波具有较好的信噪比和分辩率,能够反映出地质体沿垂直方向和水平方向的变化。
可以用波形图或彩色振幅图显示结果,同时进行运动学和动力学方面的解释分析,数据处理的基本依据为信号在空间和时间域的对比,图示直观。
目前一些地震仪器(如北京水电物探研究所的SWS 系列)已采用了特殊的数据采集技术,可以方便、快速地获得地震映像记录。
一、 野外工作方法1 测量方法在测量过程中,每次激发,在接收点采用单个或多个检波器接收。
仪器记录后,激发点和接收点同时向前移动一定的距离(或称为点距),重复上述过程可获得测线上的一条或多条地震映像时间剖面。
2 记录点的位置这种装置的记录点位于激发和接收距离的中点,反映中点两侧射线传播范围内地下的岩层、岩性的变化。
3 最佳偏移距(窗口)在地震映像数据采集中,最佳偏移距(窗口)已不仅局限于纵波反射,而是扩展为对全波列而言。
为了获得具有高信噪比和分辩率的地震映像记录,需要做试验剖面,进行干扰波调查,,分析各种波的传播规律,确定能够最好地反映探测目标的有效波,以及该有效波在时间域和空间域的最佳时空段。
在最佳偏移距(窗口)内有效波在空间距离和时间上与其它干扰波分离,信号清晰。
二、 各种波在地震映像波形图上的反映1.折射波设水平二层大地模型如图1所示,设偏移距L 大于临界距离,则时距方程为: 21cos 2V L V i z T += (7.2.1) 式中各参数如图7.2.1地质模型所示,其中i 为临界角;z 为V 1介质的厚度。
地震资料叠前偏移处理技术应用2001-05-18收到2001-06-30改回地震资料叠前偏移处理技术应用陈宝书李松康张丽焕(中海石油研究中心勘探研究院河北高碑店074010)摘要叠前偏移处理技术是解决精细速度分析和复杂构造成像的有效手段之一。
叠前时间偏移处理是近年来国内外地震资料常规处理的发展趋势,它可获得偏移归位后的速度场,适用于陡倾角构造和深部构造的准确成像。
在渤南海区垦利11-2构造上,叠前时间偏移处理的剖面较叠后偏移处理的剖面有较大改善,构造成像质量得到提高,断面及地层不整合关系较清晰,据此得到了新的更准确的解释方案。
在南黄海盆地北部凹陷应用叠前深度偏移处理技术很好地解决了逆断层的归位和下盘的准确成像。
在琼东南盆地松涛36-1构造上,叠前深度偏移处理消除了由于海底崎岖造成横向速度变化所引起的下伏地层变形和不成像问题,证实了构造北倾存在,恢复了构造的真实面貌,获得了很好的地质效果,使得该构造的解释工作有了突破性的进展。
采用先进的二维多测线深度域叠前偏移处理和成图技术成功地对乐东22-1构造中深层“模糊带”进行了精细速度反演,实现了深度构造图的成图,落实了乐东22-1构造中深层存在圈闭,为发现乐东22-1气田提供了准确的构造图,为勘探项目提供了坚实有力的技术支持。
关键词:叠前时间偏移叠前深度偏移法向速度垂直速度图偏移崎岖海底逆断层二维多测线深度域处理地震资料处理是油气勘探过程中最重要的基础工作,其处理成果的质量直接影响着勘探的全过程。
海上油气勘探,由于钻井数量很少,地震资料的应用尤为重要。
因此,搞好地震资料处理,提高处理技术水平,不断改进和提高处理质量,是我们始终追求的目标。
1 问题与思路地震资料处理工作长期存在两个难以解决的问题,即速度分析精度和复杂构造成像。
这两个问题又是相互依赖的,复杂构造的成像建立在精确的速度分析基础之上。
根据地震资料处理技术的发展现状,解决上述两个问题的最佳方法是实施叠前精细速度反演和叠前时间(或深度)偏移处理。
第七章倾斜叠加及其应用7.1 引置在1.6节中,我们已知道二维傅氏变换是将波场分解成具有各自不同的频率,并以各自不同的与垂直方向夹角传播的各个平面波分量的一种方法。
本章要讨论射线参数域并介绍把波场分离成平面波分量的其他办法。
在炮检距轴线上面应用线性时差(LMO)和振幅求和能够实现波场[例如一个共炮点道集(CSG)]的平面波分解。
这种方法叫做倾斜叠加。
倾斜叠加的基本假设是水平的层状地层。
处理是在中点-炮检距座标系中做的倾斜叠加,并用射线参量P轴代替炮检距轴(这个射线参量是水平相速度的倒数)。
一个P值范围内的一簇地震道称作倾斜叠加道集。
在射线参量中点坐标中已设计出几种处理技术。
其实例包括道内插(7.2节),倾角滤波(7.4节),多次波压制(7.5节),折射反演(附录E)及偏移和速度分析。
Taner(1977)首先提出了中点-射线参数座标系。
他论述了平面波叠加的解释应用,其中在限定的P值范围内,将若干个常数P的剖面叠加起来以便加强倾斜同相轴。
随后还研究了其他的处理方法,例如偏移(Ottolini,1982)和速度分析(Schultz和Claerbout,1978,Diebold和Stoffa 1981以及Gonzalez-serrano,1982)。
Alam和Lasocki(1981a)以及Alam和Austin(1981b)分别讨论了应用于道内插和多次波压制的可能性。
Clayton和McMechan(1981)设计了一个折射波场反演方法,它包括倾斜叠加域中的向下延拓。
McMechan和Yedlin(1981)设计了一种应用倾斜叠加变换获得散频波的相速度曲线的方法。
根据倾斜叠加道集的向下延拓,Schlttz(1982)研究出了估算层速度的技术。
现在我们研究构成一个倾斜叠加道集,通常称做τ-P道集或就称为p道集的物理特点。
这种道集中的每一道代表一个在与铅垂方向呈某一角度传播的平面波。
实际上,用炸药震源时,能量是全方位传播的(图7-1)。
由于震源和接收器之间的炮检距不同,所以反射能量以不同的角度到达不同检波器组。
炮检距愈大或反射界面愈浅上行波的角度也就愈大。
为了帮助确定建立倾斜叠加道集的方案,首先考虑怎样才能产生平面波。
图7-2显示了一条点震源线。
假设这条线上震源的激发,是使所有点源同时激励,并且每个点产生一个球面波场。
在距地面一定距离的地方,球面波前重叠并形成一个垂直向下传播的平面波。
该平面波在界面上反射并由地面上的检波器记录下来(有一些震源类型如Geoflex和Primacord,近似短的线震源)。
使用同样的点震源线(图7-3)也可产生与垂向呈所希望角度传播的平面波。
为了做到这一点,必须从这条测线的一端开始、以一个等时间延迟依次激发这些点震源。
当一个单独点震源激发时,前一个震源位置所产生的波前早已向地下传播了一定距离。
当由不同震源产生的所有球面波前重叠时,其结果是一个倾斜了的平面波前(图7-3)。
然后这个平面波就传播,由界面反射且被地面上的检波器记录。
波前的倾斜度(或平面波的传播角度)是可以控制的。
我们研究图7-4的几何图形。
在震源位置S1产生的波前到达地下点A时,在S2位置的点震源应被激发,这样就可以获得所希望的角度。
规定凡S1S2间的距离为Δx,波在介质中传播的速度为v,如果波前从S1到A点所用的时间为Δt,那么这个平面波的倾角θ为, Sinθ=vΔt/Δx (7.1) 激发震源的位置经以面Δx/Δt=v/sinθ的速度沿水平方向传播,S2位置的点震源应该按时激发,这样,我们才能在Sl震源点激发所产生的波前到达地下反射点A时赶上它。
震源位置移动所必须用的速度v/Sinθ,叫做水平相速度。
根据图7-2和图7-3的试验可以看出用以下步骤能产生与垂向呈某一角度传播的平面波。
1.在地面上将点震源排成一条直线。
2.以同一个时间延迟依次激发这些点震源。
3.将球面波前形状的响应叠加。
单个检波器(图7-3)记录这个叠加了的响应,这种响应是由界面反射后的平面波的形式。
对于一个给定的检波器位置,叠加就意味着在炮点轴线上求和。
应用互换原理,求和也可以用对给定炮点位置在检波器轴线上求和来完成。
上述讨论说明,一个CSG做为单独一个波场,怎样能分解成其平面波分量。
用检波器的排列轴代替图7-4的炮点排列轴,得到图7-5的射线路径几何图形。
与垂向呈θ角的平面波的时间延迟由方程(7.2)给出:Δt=(sinθ/v)Δx (7.2) 斯奈尔定律说明,sinθ/v即水平相速度的倒数的值在层状介质中(图7-6)沿射线路径是常数,这个常数叫做射线参量P,那么方程(7.2)可改写成:Δt=pΔx (7.3)平面波传播的角度是通过调整p值来控制的。
置p=O时,相当于平面波垂向传播。
给定层状地层的P和速度模型,就可追踪与具体的P值有关的一簇射线路径,如图7-7所示。
在层状地层中传播的平面波,叫做斯奈尔波(Claerbout,1978)。
这种平面波按照斯奈尔定律(图7-6)在每一个层的分界面上改变它的传播方向。
对于单个的p值,注意信号是在许多炮检距上记录的。
一般情况下,所有炮检距上的检波器都记录许多p值的平面波。
为了将炮检距道集分解成平面波分量,道集中的所有道的振幅必须沿着若干个倾斜路径用方程(7.3)确定的各自独有的时间延迟进行求和。
我们已讨论了把CSG波场分解成它的平面波分量。
只要没有倾角,CSG和CDP 道集的旅行时曲线是没有区别的(图7-8)。
由于一个CDP道集不是单独一个波场,所以似乎平面波分解不能应用于CDP道集。
然而水平层状地层中CDP道集和CSG的等效性为在两种类型的道集中应用平面波分解提供了理论依据。
7.2 倾斜叠加的建立通常,在炮检距域中通过倾斜路径的振幅求和人工合成平面波有两个步骤。
第一步,通过座标变换将线性时差校正值加到数据中τ= t-px (7.4)式中p为射线参量,x为炮检距,t为双程旅行时,τ为线性时差时间,LMO以后,输入端斜率为p的同相轴是平直的。
第二步是,对炮检距轴线上的数据求和得到:p(x,τ+px) (7.5)s(p,τ)=x式中 S(P、τ)表示具有射线参量p=Sinθ/v的平面波。
对各种p值重复LMO并进行求和[方程(7.5)],由原来的炮检距数据中的所有倾斜分量构成的完整的倾斜叠加道集(或P道集),就建立起来了。
倾斜叠加和波场精确地平面波分解之间是有差别的。
Treitel等人(1982)从数学上分析了平面波的分解过程,并区分了这里所述的普通倾斜叠加和正确的倾斜叠加之间的差异。
当我们论及线性震源时,常规叠加提供精确的平面波分解;当我们涉及到点震源时,正确的倾斜叠加提供精确的平面波分解。
除求和前用一个滤波算子与线性时差的波场进行褶积外,使用上述相同的步骤就产生正确的倾斜叠加。
这个算子对由点震源获得的波场,变换成线震源获得的波场,所产生的三维效应进行校正。
就所涉及的运动学来说,两种类型的倾斜叠加是等效的。
只有在讨论振幅时,它们才有所差别(Treitel,个人通信)。
图7-9为方程(7.4)和(7.5)所描述的平面波映射的示意图。
由沿p=0的水平路径在炮检距域中振幅求和开始。
这条路径在反射双曲线的顶点A附近与其相交。
因此,A点映射到(p、τ)平面上的A′点,把求和线倾斜,在B点处与双曲线相交,B点映射到B′点上。
注意,沿倾斜路径求和的主要贡献来自切点B的区域。
这个相切区叫做菲涅尔带。
对于较高的速度和较深的同相轴,菲涅尔带是较宽的。
实际上,方程(7.5)的求和可限定在菲涅尔带中。
求和必须的最陡路径是沿p=1/v 它是双曲线的渐近线,这路径相当于与船垂方向呈90°的射线。
能量没渐近线映射到p轴的C′点上。
利用上述映射,(x,t)域中的双曲线轨迹就映射成(p,τ)域中的椭园轨迹(Schultz和Claerbout,1978;参考练习7.1)。
实际上,我们不可能记录一条无限长的双曲线,也不可能记录零炮检道。
因此,在倾斜叠加域的椭园路径从A′到C′不可能是完全的。
图7-10显示了一个更为复杂的情况。
当超临界的反射C(广角反射)映射到较高的p值地方时,临界点以下的反射A和D(这些反射的入射角比临界角小)映射到较低p值的地方。
理论上,炮检距域内的直线同相轴如折射波至B,映射成倾斜叠加域中的一个点。
相反,倾斜叠加域中直线同相轴映射为炮检距域中的一个点(练习7.2)。
图7-11显示了一个野外数据的实例,主要包括水底和短程层间多次反射。
除了水底反射W外,有两个明显的一次反射P1和P2。
多次反射映射沿着收敛在P=(1/1500)s/m(水速度的倒数)的椭园轨迹。
图7-12显示的是包含线性同相轴的。
注意倾斜叠加道集上的强振幅对应于炮检距数据中观测到的导波。
在两个野外数据的实例中,都是只用正p值建立倾斜叠加道集。
因此,图7-12中炮检距域内的反向散射的能量在倾斜叠加道集中没有表现出来。
现在研究用于地震资料处理的各种定义域之间的相互关系。
研究如图7-13所示的炮检距(x,t)域内有限带宽的倾斜同相轴。
炮检距的范围是从250—5000m,道间距是50m。
这个同相轴是在f-k域(ω,kx)中沿着清楚的径向射线映射的,径向射线的斜率与水平相速度有关,用下列关系表示:ω/kx=v/Sinθ (7.6)用p=sinθ/v置换找出变换域中各个变量之间的关系。
kx =Pω (7.7)图7-13也表示倾斜同相轴映射到倾斜叠加域。
在时间方向上倾斜叠加道的一维傅氏变换给出振幅谱(p,ω).它也显示在图7-10中。
这个面描述水平相速度的)面上,沿着径向方向AA′的频率函数关系和用于导波分析(7.3节)。
在(ω,kx能量,等于(P,ω)面上沿垂直方向BB′的能量。
图7-14显示空间假频了的倾斜分量。
在(ω,k)面中观测到的重影是由于同x相轴不恰当的空间采样引起的。
注意,未假频和假频了的两个分量(分别为1和2)都映射成单独一个单个P迹线。
我们希望把空间假频了的部分映射到一些负P值上。
如果是这种情况的话,那么假频范围(21Hz到42Hz)在只包含正p值的(p,ω)平面上不存在(图7-14是一个单一倾斜;在图7-17中讨论倾斜范围的重建)。
一旦在倾斜叠加域中实行了特殊处理,利用逆映射重建炮检距域的数据:Thorson(1978)提供了详细的重建步骤。
为了严格地进行振幅恢复,在逆映射前应用P滤波。
通过用频率绝对值乘以每个倾斜叠加道的振幅谱来完成这项工作。
这点与包括偏移积分表达式(方程4.5)中的求和之前的波场微分有点相似。
倾斜叠加处理流程表1炮检距数据P(x,t)↓从(x,t) 到(P,τ)步骤1:在指定的p值应用LMO[方程7.4]步骤2:在炮检距上求和[方程7.5]步骤3:对一个p值范围重复步骤1和步骤2输出是倾斜叠加,S(p,τ)↓在倾斜叠加域中应用所要求的处理如,反褶积,时变切除↓从(p,τ)到(x,t)步骤1:应用ρ滤波步骤2:对指定炮检距值应用反LMO步骤3:在ρ范围内求和步骤4:对于炮检距范围重复步骤2和步骤3,输出是倾斜叠加处理后的炮检距数据图7-15显示人工合成的炮检距道集,相应的倾斜叠加道集和除ρ滤波以外未作任何其他处理的重建的炮检距道集。