消除海底起伏影响的海洋地震波场正反向延拓(英文)

英文篇名：Downward and upward continuation of 2-D seismic data to eliminate ocean bottom topography's effect
中文刊名：Applied Geophysics
英文刊名：应用地球物理(英文版)
作者：王祥春 ; 夏常亮 ; 刘学伟
英文作者：Wang Xiang-Chun1 ; Xia Chang-Liang2 ; and Liu Xue-Wei11.Key Laboratory of Geo-detection(China University of Geosciences ; Beijing) ; Ministry of Education ; Beijing 100083 ; China.2.Overseas Business Department of Geophysical Research Institute ; BGP of CNPC ; Zhuozhou 072750 ; China.
中文关键词：波动方程 ; 延拓 ; 地形起伏 ; 变换
英文关键词：wave equation ; continuation ; ocean bottom topography ; transformation
出版日期：2010-06-15
机构：地下信息探测技术与仪器教育部重点实验室(中国地质大学);中石油集团东方地球物理勘探有限责任公司研究院海外业务部;
年：2010
期：02
出版单位：Applied Geophysics

摘要

为了解决海底起伏变化对地震波场的影响问题,本文提出将(x-z)域中的曲网格映射成(ξ,η)域中的矩形网格,推导出(ξ,η)域中的二维标量声波方程,根据推导出来的波动方程采用逆时有限差分法将海面上采集到的地震波场在(ξ,η)域中向下延拓至海底面,延拓时采用海水的速度,然后采用顺时有限差分法将延拓后的地震波场再反延拓到海面上,延拓时采用海底面以下地层的速度,从而消除了海底起伏带来的负面影响。模型及实际地震资料的计算分析表明该方法不但能够校正由于海底起伏所引起的海底面下地层反射波场的不连续性还能够校正由于海底起伏所引起的地震波的动力学特征的变化。对延拓前后的地震波场进行速度反演,延拓后反演的地层速度比延拓前反演的地层速度的精度提高很多,延拓前后地震波场的叠加剖面对比表明该延拓方法能够明显提高地震波场的成像质量。
In order to eliminate the effect of ocean bottom topography on seismic wave field,we transformed curved(x,z)coordinate system grids into rectangular(ξ,η)coordinate system grids and derived a 2-D scalar acoustic wave equation in theξ,ηdomain.The seismic wave field collected at the sea surface was downward continued to the ocean bottom by the inverse finite difference method with the water velocity and then was reversely continued to the ocean surface by the finite difference method using the layer velocity from just below the ocean bottom in the(ξ,η)domain.Simulation calculations and practical application show that this method can not only remove the reflection travel time distortion but also correct the dynamic parameter changes caused by the ocean bottom topography.The inverted velocity after wave field continuation is much more accurate than before continuation and the image section was greatly improved compared to the original wave field.

引文

Ai,Y.S.,and Zheng,T.Y.,1998,Wavefield continuation depth filtering:Chinese Journal of Geophysics,41(3),400-407.
    Berryhill,J.R.,1979,Wave-equation datuming:Geophysics,44(8),1329-1344.
    Bevc,D.,1997,Flooding the topography;wave equation datuming of land data with rugged acquisition topography:Geophysics,62(5),1558-1569.
    Fornberg,B.,1998,The pseudospectral method:Accurate representation of interfaces in elastic wave calculations:Geophysics,53(5),625-637.
    Hestholm,S.,and Ruud,B.,1994,2D finite-difference elastic wave modeling including surface topography:Geophysical Prospecting,42,371-390.
    Hestholm,S.,and Ruud,B.,1998,3-D finite-difference elastic wave modeling including surface topography:Geophysics,63(2),613-622.
    Liu,Y.,Li,C.C.,and Mou Y.G.,1998,Finite-difference numerical modeling of any even-order accuracy,Oil Geophysics Prospecting,33(1),1-10.
    Nielsen,P.,and Skovgaard,O.,1994,Using the pseudospectral technique on curved grids for2D acoustic forward modeling:Geophysical Prospecting,42,321-341.
    Tessmer,E.,Kosloff,D.,and Behle,A.,1992,Elastic wave propagation simulation in the presence of surface topography:Geophysics.J.Int.,108,621-632.