1.一種地震勘探中的非水平地表直接延拓束偏移方法,其特征在于包括:獲取非水平真地表高程和建模速度場;讀入炮記錄、高程和速度場文件,確定頻帶、初始束寬和射線采樣;利用地表高程和速度場對單個束進行成像;所有束成像結(jié)果疊加得到單炮成像結(jié)果;所有單炮成像結(jié)果疊加得到數(shù)據(jù)成像結(jié)果;輸出最終的非水平地表束偏移結(jié)果。
2.根據(jù)權(quán)利要求1所述的地震勘探中的非水平地表直接延拓束偏移方法,其特征在于包括:
(1)對地震資料進行速度建模,獲得非水平真地表高程及疊前深度域速度場;
(2)讀入實際炮集記錄、非水平真地表高程以及疊前深度域速度場三個數(shù)據(jù)文件,確定成像頻帶范圍為10Hz~50Hz,然后根據(jù)速度場的幾何平均值和參考頻率確定束初始寬度,最后利用頻帶范圍和初始寬度確定束間隔以及初始射線參數(shù)間隔;
(3)根據(jù)束初始寬度和束間隔劃分成一系列的“束”,然后針對每個束進行如下操作:利用起伏地表高程和近地表速度對地震炮集記錄按照公式(1)所示進行加窗局部傾斜疊加,分別利用運動學(xué)和動力學(xué)射線追蹤計算震源和檢波器位置處出射的“束”的走時信息和振幅信息,得到震源到檢波點的延拓波場D(ω,x,xs)和檢波點到震源的延拓波場U(ω,x,xs),然后將對應(yīng)的延拓波場傾斜疊加道應(yīng)用公式(2)所示的反褶積成像條件進行成像,完成一個束的偏移;
其中,G(L,ω,pLx)為中心位于L處的“束”內(nèi)地震記錄的加窗局部傾斜疊加,ω和 ωr分別為頻率和參考頻率,Ω為非水平地表面,β為“束”出射角,α為地表傾角,U(ω,xr,xs)為檢波點xr到震源點xs的延拓波場,w0為束初始寬度,pLx和pLz為束射線的初始慢度的x和z分量,h為檢波點xr與束中心L的高程差;R(x,xs)為成像值,D(ω,x,xs)和U(ω,x,xs)分別為震源到檢波點的延拓波場和檢波點到震源的延拓波場,π為圓周率,*表示共軛;
(4)對每個“束”進行循環(huán),得到所有束的偏移結(jié)果,疊加后得到單個炮集記錄的成像結(jié)果;
(5)對每個單炮記錄進行循環(huán),得到所有單炮的偏移結(jié)果,疊加后得到整個數(shù)據(jù)的成像結(jié)果。