[发明专利]真地表起伏叠前深度域双程波成像方法

摘要

本发明公开了一种真地表起伏叠前深度域双程波成像方法，具体步骤包括：①以起伏地表基准面形态为基础，获取地震单炮数据所在面的真地表高程；②重新定义To值并修正速度模型；③构建震源子波采用优化差分方法实现震源波场正演并应用真地表边界条件；④记录数据规则化、沿时间方向逆时外推，应用起伏地表边界条件；地震波场相关以获取深度域地震成像；该方法通过记录数据在真地表进行深度域波场外推及成像，可克服常规技术下因地表高差较大导致的波场畸变，构造形态不真实的缺陷，可获得复杂地表条件下的精确地震成像及其更真实的地质构造形态。

主权项

1.一种真地表起伏叠前深度域双程波成像方法，其特征在于：该方法包括如下步骤：(1)、对起伏地表基准面进行处理，获取地震数据CMP面高程值，获取炮点、检波点真地表高程值及坐标信息，定义激发点空间位置；DATUM为固定基准面，CMP surface为CMP道集及速度所在面，Turce Surface为真地表高程面；(2)、构建震源子波采用差分方法实现震源波场正演，并应用真地表边界条件；w(t)=1w2-w1{w2exp[-(πw2t)2]-w1exp[-(πw1t)2]}---(1)]]>根据地震记录频率和公式(1)构建震源子波，利用高阶有限差分方法公式(2)正演震源波场到最大记录时间，在炮波场正推过程中把整个记录时间T分成若干段，每延拓一步应用起伏地表自由边界条件。构建滤波函数F_il(x，z)对起伏地表以上波场进行滤波；同时存储每一段起始时的两层波场数据到机器内存；Pi,j,kn+1=2Pi,j,kn-Pi,j,kn-1+c2Δt2Δx2Σm=1Mam(Pi+m,j,kn+Pi-m,j,kn-2Pi,j,kn)]]>+c2Δt2Δy2Σm=1Mam(Pi,j+m,kn+Pi,,j-m,kn-2Pi,j,kn)]]>+c2Δt2Δy2Σm=1Mam(Pi,j+m,kn+Pi,,j-m,kn-2Pi,j,kn)]]>+c2Δt2δi,j,k---(2)]]>式中δ_i，j为震源项，f(n)为离散的震源函数，i_s，j_s，k_s分别为震源点对应的x，y，z方向的网格点序号；(3)、以炮集记录为输入目标，对记录波场利用有限差分方法沿时间方向逆时外推，公式(4)，同时应用起伏地表自由边界条件，公式(5)；获取起伏地表的逆时外推记录波场；Pi,j,kn-1=2Pi,j,kn-Pi,j,kn+1+c2Δt2Δx2Σm=1Mam(Pi+m,j,kn+Pi-m,j,kn-2Pi,j,kn)]]>+c2Δt2Δy2Σm=1Mam(Pi,j+m,kn+Pi,,j-m,kn-2Pi,j,kn)]]>+c2Δt2Δz2Σm=1Mam(Pi,j,k+mn+Pi,j,k-mn-2Pi,j,kn)]]>+c2Δt2δi,j,k---(4)]]>起伏地表自由边界条件：P_Ω(x，d)＝0     (5)(4)、以应用起伏边界条件的震源波场和记录波场为输入，应用波场相关成像条件，以获取该时刻的成像结果。