波动方程有限元叠前逆时偏移
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
为研究起伏地表、复杂构造以及复杂速度分布条件下地震数据的精确成像方法,本文采用有限元法完成观测波场的声波方程逆时外推,并提出了以节点波场最大振幅出现时间作为成像时间的最大振幅成像条件。该成像条件不同于激发时间成像条件,它基于波动理论而不是射线理论,可使能量更好地聚焦。有限元法采用六节点三角形单元,对单元内的场和速度均作二次插值以提高模拟精度;用对角化的集中质量矩阵代替协调质量矩阵,避免了矩阵求逆,提高了计算效率。通过对带地形的凹陷模型、Mormousi模型及带地形的Mormousi模型等三个理论模型的计算,结果表明该方法能够有效处理复杂地质条件下的地震成像问题,对于复杂地质体中的断层、断块、背斜构造、高速体和目标体成像都比较清晰,位置准确。尤其对于基于起伏地表直接进行偏移的模型数据也得到了较好的成像结果。
In order to study the preciously imaging method of seismic data in conditions of relief surface,complex structure and velocity distribution,the paper used finite-element algorithm to finish the wave-equation reverse-time extrapolation of surveyed wavefield,and proposed taking the time the max amplitude appeared in node wavefield for the imaging condition of max amplitude,which is different from the imaging condition of shooting time.This imaging condition is based on the wave theory rather than on ray theory,which can focus energy better.The finite-element method used 6-node triangular element,quadratic interpolation of field and velocity in element is carried out to improve the precision of simulation; Diagonal lumped mass matrix was used to replace consistent mass matrix to avoid matrix inverse and improve computation efficiency.It's shown by computing three theoretical models: sag model with topography,Marmousi model and Marmousi model with topography that the method can effectively handle the issue of seismic imaging in complex geologic condition.The imaging in such complex geologic bodies as fault,faulted block,anticline structure,high-velocity body and targets is characterized by identification and accurate position.Especially,the imaging results of model data based on ragged surface by direct migration are also better.
引文
[1]Schneider W A.Integral formulation for migration in two and three di mensions.Geophysics,1978,43(1):49~76
    [2]Gazdag J.Wave equation migration with the phase-shift method.Geophysics,1978,43(7):1342~1351
    [3]Stolt R H.Migration by Fourier transform.Geophys-ics,1978,43(1):23~48
    [4]Claerbout J F.I maging the Earth’s Interior.Blackwell Scientific Publications Inc,1985
    [5]马在田.地震成像技术——有限差分法偏移.北京:石油工业出版社,1989
    [6]Stoffa P L,Fokkema J T,de Luna Freire R Met al.Split-step Fourier migration.Geophysics,1990,55(4):410~421
    [7]Ristow D,R櫣hl T.Fourier finite-difference migration.Geophysics,1994,59(12):1882~1893
    [8]Chang W F,Mc Mechan G A.Reverse-ti me migration of offset vertical seismic profiling data using the exci-tation-ti me i maging condition.Geophysics,1986,51(1):67~84
    [9]Teng Y C,Dai T F.Finite-element prestack reverse-ti me migration for elastic waves.Geophysics,1989,54(9):1204~1208
    [10]Mc Mechan G A,Chen H W.I mplicit static corrections in prestack migration of common-source data.Geo-physics,1990,55(6):757~760
    [11]Raiaskaran S,Mc Mechan G A.Prestack processing of land data with complex topography.Geophysics,1995,60(6):1875~1886
    [12]Zhu J,Lines L R.Comparison of Kirchhoff and re-verse-ti me migration methods with applications to prestack depth i maging of complex structures.Geo-physics,1998,63(4):1166~1176
    [13]Mc Mechan G A.Migration by extrapolation of ti me-dependent boundary values.Geophysical Prospec-ting,1983,31:413~420
    [14]Vidale J E.Finite-difference travel ti me calculation.Bull Seis Soc Am,1988,78:2062~2076
    [15]Nichols DE.Maxi mumenergytravelti mescalculatedin the seismic frequency band.Geophysics,1996,61(1):253~263
    [16]Shin C et al.Wave equation calculation of most ener-getic travelti mes and amplitudes for Kirchhoff pres-tack migration.Geophysics,2003,68(6):2040~2042
    [17]Zienkiewicz O C,Taylor K L.The Finite Element Method,Fifth Edition.Butterworth Heinemann,2000
    [18]王勖成.有限单元法.北京:清华大学出版社,2003

版权所有:© 2023 中国地质图书馆 中国地质调查局地学文献中心