波场合成叠前偏移层速度反演.pdfVIP

波场合成叠前偏移层速度反演 张文生1张关泉1朱海斌2 (I中国科学院数学与系统科学研究院，计算数学与科学工程计算研究所，科学与工程计算国家重点实验室，北京 100080；2中国科学院地质与地球物理研究所，北京100029) 常规的正常时差速度分析，由于采用简单的层状速度模型，当模型复杂时会有较大误差，叠前速度分析可用 来改善这一情况。叠前速度分析的主要方法有：深度聚焦分析(DFA)、共反射点道集(CRG)或共成像点道集(CIG) 剩余时差分析和共聚焦点分析(CFP)“。31。在各种叠前偏移速度分析方法中，CRG道集提供一个重要准则，即在CRG 道集上的时差曲线向上或向下弯曲可提供偏移速度比真正介质速度高还是低的信息。通常CRG道集由炮记录和共 偏移距剖面的叠前偏移得到。这里，我们用平面波合成叠前深度偏移来产生用于剩余速度分析CRG道集。与其他 叠前偏移速度方法所得到的CRG道集相比，通过波场合成成像所得的CRG道集能提供分辨率较好的剩余速度面板。 当地下的一个位置被正确成像时，其成像结果应当在同一个深度而不论震源的位置和形态，因此，我们可以 将平面波合成震源作为震源，并作为入射角而不是偏移距来确定一个CRG道集。假定介质的速度为v，反射点深 度为z，偏移速度为V。，由于速度不准所导致的新的反射点的深度为Z卅，经过推导，可得如下不确定性的关系 式(剩余时差方程) sine (1) 口2z：+bz。=(y2—1)x2+4zjy2+bzo 其中y=v／v。，是介质速度与偏移速度之比；≯是反射面倾角：‰是源点到反射面的垂直距离；x、a和b是 与平面波入射角、反射面倾角≯有关的系数，有精确表达式。方程(1)就是深度一y方程，它给出了视深度Z。与 真实深度z之间的关系，它们通过参数y联系在一起。当偏移速度与介质的速度相等即y=1时，参数y与偏移 距和平面波的入射角无关。当y不等于1时，将出现深度上的误差，但在零偏移距处，关系式(1)变为z。=yz， 即波场用不正确的速度传播到正确的深度等同于用正确的速度传播到错误的速度，这就解释了为什么偏移速度分 析不能通过合成一个平面波来完成，正如不能由零偏移距偏移来完成一样。 ClIP道集中的曲率随深度而变化，因此参数，，不是一个常数而是深度的函数。ClIP道集中的曲率可以用于估 计速度模型。我们采用搜索的方法来度量曲率。即在每～个深度，在CRP道集中将数据沿由关系式(1)确定的轨 迹求和。求和对某一范围的y进行，最大的求和值所对应的y值即与曲率相匹配。由CRP道集可以量化剩余时差， 剩余时差再转化成关于y的速度图。由此确定的速度，称为剩余速度，它是一种平均速度。由于在y面板上所拾 取的是最大能量团所对应的，，值，从而在各个深度上的各个能量团之和应最大，所以根据均方根速度(其中包含 层速度的信息)与平均速度的关系，可建立一个关于优化层速度的目标函数。通过优化目标函数，可导出一个反 演层速度的线性代数方程组，并用预处理共轭梯度方法来反演层速度。在数值计算中，均采用了MPI并行编程，进 一步提高了计算效率。 本研究由国家重点项目和国家自然科学基金项目资助。 参考文献 S [1]Mackayet．a1．Geophysics，57(也)，1992．． CF [2]Lafondet．a1．Geophysics，58(1)，1993． [3]A1一Yahyaet．a1．Geophysics，54(6)，1989． ·536·