[發明專利]真地表起伏疊前深度域雙程波成像方法

無效

• 申請號: 201110459375.0
• 申請日: 2011-12-31
• 公開（公告）號: CN102590857A
• 公開（公告）日: 2012-07-18
• 發明（設計）人: 劉文卿;雍學善;胡自多;王宇超;劉紅偉;曾華會;徐新榮
• 申請（專利權）人: 中國石油集團西北地質研究所
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 蘭州振華專利代理有限責任公司 62102
• 代理人: 陳寶玲

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)對震源波場和記錄波場進行濾波，滿足起伏地表邊界應力為零的邊界條件，以獲取基于起伏地表的地震波場。