《共聚焦点(CFP)偏移速度估计和AVO分析方法》第3章方法技术和算法.ppt

对速度-深度模型的平滑我们采用了三次样条和样条函数。上述迭代方程是基于参数化速度函数推出的。若给定不同参数化速度函数，迭代方程的系数矩阵会有所不同。但它们的方法原理都是共同的，都是基于最小二乘原理的广义线性反演。为此，我们可以根据不同的情况，选用不同的参数化速度函数，实现最优约束参数反演。(6-14)2.实现步骤采用层剥离法，通过走时分析（DTS）实现速度估计。从顶层开始，对每一层应用如下步骤。(1)计算合成聚焦算子和逆时聚焦算子，计算方法包括射线追踪和波动方程法等，在此，我们采用了有限差分法走时计算；(2)合成聚焦算子与炮记录做褶积生成CFP道集（聚焦第一步）；(3)CFP道集与逆时聚焦算子做互相关计算DTS（聚焦第二步），判断DTS是否为零；(4)DTS0时，应用迭代反演更新速度，用更新后的速度重复(1)-(3);(5)实现CFP偏移。七.基于拉东域的正确聚焦算子的生成和基于正确聚焦算子的偏移速度建模的思路1.相位误差的对称性2.走时分析-------速度估计1）生成CFP道集（聚焦第一步）2）估计正确的聚焦算子（CFP偏移速度分析）3）计算CFP叠加（聚焦第二步）DTS=04）重复1）------3）得到共聚型偏移结果七.基于拉东域的正确聚焦算子的生成和基于正确聚焦算子的偏移速度建模的思路图7-1.误差对逆时聚焦算子和对应聚焦点响应的影响.a)逆时聚焦算子（无误差）；b)聚焦点响应（无误差）；c)–10%速度误差；d)+10%速度误差；e)–20%深度误差；f)+20%深度误差；g)–10%速度误差++20%深度误差；h)–10%速度误差+-20%深度误差图7-2:单层的算子更新过程图7-3:倾斜反射层的算子更新步骤图7-4正确速度模型下的逆时算子和CFP道集图7-5+10%速度误差时的逆时算子和CFP道集图7-6正确速度模型下的逆时算子和CFP道集在RADON域中的结果图7-7+10%速度误差时的逆时算子和聚焦点响应在RADON域中的结果图7-8在RADON域求出的新算子和其RADON反变换后的结果图7-9更新算子边缘效应切除后的结果和正确的逆时聚焦算子图7-10更新后的算子与其CFP道集做零延迟互相关的结果正确速度模型下逆时聚焦算子和CFP道集做零延迟互相关的结果3.基于正确聚焦算子进行速度估计的思路1)速度扫描法a.给定速度扫描范围，利用常速或变速射线追踪、有限差分等方法生成估算的聚焦算子；b.通过振幅、能量或相似系数准则，比较和拟合估计的和正确的聚焦算子，由此实现偏移速度建模。2)迭代反演法a.给定一初始速度模型，基于迭代反演将估算和正确的聚焦算子之间的走时误差转换为速度和深度误差，实现速-深模型的更新；b.反复迭代，直到满足精度要求为止。2模型试算计算中选取了一均匀模型。其速度为2000m/s，网格大小为100*100，当间距与之比小于3时，计算的相对走时误差小于0.05%；当网格大小为80*100时，网格间距之比越小，计算的相对走时误差也越小。图3是一层状介质模型。上、下层速度分别为1000m/s和2500m/s，网格大小为100*80，纵横向网格间距分别为2m和4m，震源位于（40,160）点处的初至波等时线图重叠在该层状介质上。图4为均匀介质中含一高速体，背景速度为1000m/s，高速体速度为2500m/s,网格大小为80*100，纵横向网格间距分别为2m和4m，震源位于（160,0）点处的初至波等时线重叠在该高速体介质上。图4含高速体的均匀介质和重叠在其上的震源位于(100,40)点的等时线图。背景速度为1000m/s，高速体速度为2500m/s图3层状介质和重叠在其上的震源位于(10,100)点的等时线图。注意，上、下层速度分别为1000m/s和2500m/s。2500m/s1000m/s1000m/s2500m/s横向距离(x/m)深度(z/m)深度(z/m)横向距离(x/m)3Kirchhoff叠前深度偏移基于Marmousi模型进行了初至波走时的计算。Marmousi模型是一比较复杂的2-D地质模型，如图5(a)所示。模型上部主要有三大断裂组成；中部有一盐丘；