射线层析与波动方程联合成像的速度和深度构造同时成像.docxVIP

射线层析与波动方程联合成像的速度和深度构造同时成像 一般来说，波动方程偏移是用于获得结构成像（深度成像）的方法。 1 修正前的速度解决方案 当偏移中的速度模型正确时, 各炮在各反射点正确成像, 得到清晰的反射层 (图1a) 。当速度错误时, 则各炮在各反射点成像不正确, 就得不到清晰的反射层 (图1b) 当出现图1b所示的错误成像时, 根据各共成像点 (a、b、c、d) 道集的各道成像深度偏差, 算出各道的速度校正量, 作层析成像速度反演, 即按其射线反向路径均匀修正速度。对于图1b中的绿色区域, 有多条射线交叉, 要作多次修正, 取各次修正的加权平均或算术平均作为修正结果。把修正后的速度作偏移的速度输入, 再作偏移处理, 偏移后再作层析速度反演。以此反复迭代, 最终求得正确的速度和深度。 针对单斜模型 (图2) 来说明实现原理的方法步骤。 (1) 给定一个初始速度模型, 可以是叠加速度或其他常规处理得出的速度。 (2) 对各炮集作波动方程叠前深度偏移。 (3) 根据具体情况确定共成像点CIP的地面位置 (图2a) , 从各炮集的偏移结果中抽取各点的共成像道集。如图2b所示, 在3种速度情况下, 在44点位置抽取共成像道集, 道集中的各道依据其所在炮集的炮点与共成像点的距离大小顺序排列。 (4) 对某一共成像道集m, 由浅入深分析共成像道集同相轴。 (1) 首先设定一个小的方差值σ. (2) 由浅入深分析, 如果是第一次迭代, 就从第一条同相轴开始, 否则按选定的同相轴做。e表示同相轴由浅入深的序号, 计算各共成像道集中各道的深度D (5) 改变共成像点, 重复 (4) 中的各步, 对于同一条深度同相轴在不同共成像点所修正的速度v (6) 用新的速度, 重复 (2) 到 (5) 步的过程。 为突出射线层析与波动方程联合主题, 在射线层析计算中始终以直线代替弯曲射线, 这种替代在地层速度变化较大的情况下将产生误差, 此时可直接以弯曲射线作反向射线追踪 2 模型试验 实验设计的模型如图3a所示, 首先用单程声波方程多次波可控制的共炮记录模拟法 (定位原理 3 速度和构造同时成像的可行性之路 把波动方程叠前深度偏移可正确进行深度构造成像和射线追踪层析成像, 并与反演细微速度的优点结合起来, 由深度偏移作构造成像和计算深度偏差, 进而由深度偏差计算速度修正量, 用这个修正量分层去修正速度。由此反复迭代, 达到速度和构造的同时正确成像。这条思路避开了射线法层析成像迭代反演要给定深度 (构造) 的假设, 充分发挥了波动方程偏移对深度成像精确的优势, 把焦点集中在解决射线法反演速度问题上来, 减少了相应的近似。同时, 波动方程和射线层析理论有机的结合, 本身就是一个相互印证的过程, 是速度和构造同时成像的一种可行之路。 本文只讨论了直射线追踪反演, 弯曲射线问题有待继续深入研究。