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

权利要求书

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)、以应用起伏边界条件的震源波场和记录波场为输入，应用波场相关成像条件，以获取该时刻的成像结果。

2.根据权利要求1所述的真地表起伏叠前深度域双程波成像方法，其特征在于，真地表起伏叠前深度域双程波成像方法，处理技术包括：

a、起伏地表获取：对地震记录的基准面形态进行分解，获取地震记录对应的起伏地表面；

b、速度模型修正：建立地震速度模型，定义时间域、深度域基准面形态图时空关系、修正速度模型，获取起伏地表面对应的速度模型；

c、震源子波构建与波场正演：构建合适的震源子波，在震源子波基础上，进行波场正演，获取震源波场；

d、数据规则化：地震记录数据处理，以获取所述时间域规则的地震数据做为记录波场逆时外推的基础；

e、起伏地表自由边界处理：处理起伏地表边界问题，构建滤波函数F_il(x，z)对震源波场和记录波场进行滤波，满足起伏地表边界应力为零的边界条件，以获取基于起伏地表的地震波场。