波動方程疊前深度偏移直接產(chǎn)生角道集

波動方程疊前深度偏移直接產(chǎn)生角道集

1基于角道集的程疊前深度偏移成像技術隨著復雜結構成像和巖石勘探的深入，在地震資料的重疊前偏移圖像中，根據(jù)不同炮或不同位移的位移結果，得到不同角度下反射點的反射強度。此外，通過avo或ava響應，直接識別液體和氣體的含量已成為一個越來越受關注的環(huán)節(jié)。波動方程疊前深度偏移是用較準確的波動理論來描述地震波的傳播過程,能正確模擬復雜上覆地層導致的多值走時等復雜的波傳播現(xiàn)象和地震波傳播的幾何擴散現(xiàn)象,因此既可對復雜構造正確成像,又有望得到正確反映地下構造反射強弱的成像幅值.但服務于復雜構造成像的波動方程疊前深度偏移多是通過分炮偏移實現(xiàn)的,且每炮偏移實現(xiàn)時需要對原始炮集做鑲邊處理以滿足偏移孔徑的要求;這種數(shù)據(jù)集組織方式易于得到疊加剖面,但由這些炮集偏移結果抽取隨偏移距變化的共反射點(CRP)道集將產(chǎn)生巨大的硬盤存儲需求.此外,既然真正反映地下構造物性參數(shù)對比的是隨入射角變化的角度域成像道集,人們就希望能由分炮偏移結果直接得到角道集.已發(fā)展了一系列通過炮域波動方程疊前深度偏移產(chǎn)出角度域成像道集的方法.這些方法可歸納為如下兩類:一是在應用成像條件前對入射和反傳波場進行局部平面波分解,對每對平面波分別成像本文在文獻[7]所發(fā)展的單程波算子地震波入射角計算方法的基礎上,通過計算震源入射波場與地層界面夾角,發(fā)展了由炮域偏移的結果直接通過部分疊加得到角道集的方法;該方法通過應用“保幅”的反褶積成像條件和考慮角度區(qū)間內(nèi)累加的炮數(shù),解決了炮點覆蓋不均勻?qū)е碌某上穹嫡`差問題.與前述兩類方法相比,這一方法所需的計算量及內(nèi)存幾乎可以忽略,角道集的有效頻率范圍也與偏移剖面相同,因而更適于現(xiàn)行的計算能力和工業(yè)資料的實際情況.本文首先簡述基于單程波算子的地震波入射波之波前面法線、地層界面傾角估計以及以此為基礎的地震波入射角計算方法,然后建議了一種考慮炮點覆蓋不均勻的基于反褶積成像條件的直接累加生成角度域成像道集的方法,最后給出了實際應用的計算流程和基于角道集改善地震成像效果的實用方法.文中算例表明,這一炮域波動方程深度偏移直接產(chǎn)生角道集方法既能剔除復雜上覆地層的影響,又可補償觀測系統(tǒng)非均勻覆蓋對成像幅值的影響;而與常規(guī)炮域波動方程疊前深度偏移相比,并不增加過多內(nèi)存需求與計算量.2反射角計2.1基于規(guī)則波動的入射角方程單程波算子是通過對標準的聲波方程進行分解而得到的,一般情況人們總是沿垂直方向分解單程波算子.單程波算子可正確模擬復雜介質(zhì)中地震波的單向傳播,與標準的聲波(雙程波)方程相比,單程波算子沿其分解方向不發(fā)生地震波反射.由于這一特征,單程波算子所得到的簡諧波場就對應著地震波傳播過程中的波前面.地震波場的波數(shù)是與波前面的法線方向相聯(lián)系的,而波場的空間導數(shù)又與波數(shù)存在著對應關系,因此,我們可利用求取空間導數(shù)前后簡諧波場的變化,簡單求得全場的地震波入射角方向.由單程波方程可知,頻率-波數(shù)域下行的波場P式中k=ω/c(x,y,z)為波數(shù),c(x,y,z)是空間各位置處的介質(zhì)波速,角度θ式中符號Re代表求復波場值的實部.若令脈沖響應由此可得到(x2.2ensors的方向定義基于局部結構的張量分析(localstructuretensors)為簡單信號提供了一個較簡潔的描述其方向的方法.這里所謂的簡單信號,是指該信號在局部范圍內(nèi)只能沿一個方向變化.對于這類信號,可以把它的局部方向定義成一個張量T式中f3觀測系統(tǒng)非均勻覆蓋的消除方法—“保幅”角道集的形成就反射波偏移成像而言,震源入射波場與反射界面法線的夾角,即入射角,就是反映偏移道集AVA特征的角度.基于此角度的角道集完全反映了由界面波阻抗和泊松比對比決定的隨角度變化的界面反射特征基于第2節(jié)獲得的反射角,可以按反射角對炮域偏移結果直接累加來得到角道集.有學者從理論上證明本文因此建議了一種可同時處理幾何擴散效應和觀測系統(tǒng)非均勻覆蓋的“保幅”的角道集生成方法.角道集是按反射角對炮域偏移結果直接累加形成的;在偏移成像中用反褶積成像條件來補償幾何擴散,同時應用統(tǒng)計按角度累加過程中落入不同角度區(qū)間內(nèi)的炮數(shù)的方法來剔除觀測系統(tǒng)非均勻覆蓋的影響.炮域偏移的反褶積成像條件可表達為式中A疊加是壓制隨機噪音的有效手段,對地震資料整個處理流程,疊加起到了關鍵的作用;但疊加也改變了地震資料的幅值特征,特別是在觀測系統(tǒng)非均勻覆蓋的情況下.解決觀測系統(tǒng)非均勻覆蓋是地震資料處理中一直受到關注的問題,目前尚沒有非常有效的方法.在生成角度道集的過程中,為獲得較高信噪比的角道集,也需要采用疊加處理,即保證有足夠多的反射角(每個角度對應一個炮點)落在一個角度區(qū)間內(nèi).對等間距的角度區(qū)間而言,即使是規(guī)則的三維觀測系統(tǒng),由于(各向同性情況下)不考慮方位角,落在各個角度區(qū)間內(nèi)的炮點數(shù)也是不相同的.由于采用反褶積成像條件,可認為每個單炮偏移結果對應某一個角度的(有限帶寬)的反射系數(shù);因此,在累加過程中統(tǒng)計落入不同角度區(qū)間內(nèi)的炮數(shù),可籍此剔除非均勻覆蓋對角道集幅值的影響.具體實現(xiàn)如下:定義角度采樣間距Δγ,對每一成像點定義一個一維數(shù)組,數(shù)組長度小于90/Δγ;對全部炮循環(huán),對每炮的偏移結果,讀取由上節(jié)得到的該炮點震源對各成像點的地層界面法向量夾角γ,計算l=1+int{γ/Δγ},將各成像點的幅值累加到對應數(shù)組的l元素上,同時統(tǒng)計累加到元素上的炮數(shù);完成對全部炮循環(huán),然后用統(tǒng)計得到的各元素上的累加炮數(shù)除各成像點對應數(shù)組中的各元素值.這樣就得到各成像點隨角度變化的成像幅值.需要指出的是,這一操作并不能簡單推廣應用到常規(guī)的疊加處理;以炮域偏移的疊加過程為例,由于不能保證任一炮記錄的偏移結果對某一同相軸有實質(zhì)的貢獻,簡單的統(tǒng)計、剔除并不能消除非均勻覆蓋的影響.4基于高效可獲得性的角道集成像結合炮域波動方程偏移流程以及集群式并行機的軟硬件分布特點,提出了如圖1所示的一種直接由炮域偏移得到角道集的實現(xiàn)流程.這一流程具有如下幾個特點:首先,就各成像點震源入射波場與地層界面的夾角的計算環(huán)節(jié)而言,因其使用了與波動方程炮域疊前深度偏移處理完全相同的偏移速度模型和單程波算子,因此得到的角度更真實地反映了入射波與反射波之間的夾角;而其所需計算量僅為單個頻率的震源波場延拓加上求導及平滑、插值等處理,與疊前偏移的全頻帶波場延拓相比,幾乎可以忽略;其次,這一環(huán)節(jié)主要的計算(圖1中兩虛線箭頭中間部分所示)均是在集群并行機計算節(jié)點上并行進行的,因而可獲得較高的計算效率.就這一實現(xiàn)流程的內(nèi)存需求而言,各成像點震源入射波場與地層界面的夾角的計算環(huán)節(jié)所需內(nèi)存與單炮偏移基本相同,且這一操作與單炮偏移是先后進行的,不需要額外的內(nèi)存,這克服了空間移動或時間移動的成像條件產(chǎn)生角道集方法巨大內(nèi)存需求導致不能應用于大規(guī)模工業(yè)數(shù)據(jù)的障礙.本文流程是先對單炮偏移結果直接在各個計算節(jié)點部分疊加,然后在主節(jié)點對部分疊加結果進行疊加產(chǎn)生各成像點的角道集,這一實現(xiàn)方法規(guī)避了首先由炮集偏移結果抽取共成像點道集,然后再結合入射角信息將共成像點道集轉(zhuǎn)化為角道集時產(chǎn)生的巨大文件存儲及數(shù)據(jù)交換需求.與現(xiàn)行各類波動方程偏移生成角道集的方法相比,本文提出的方法幾乎不增加額外計算量和內(nèi)存需求,可容易與現(xiàn)行炮域偏移流程結合,獲得有較好幅值特征的角道集.本文方法獲得的角道集,除服務于疊前巖石物理參數(shù)反演,還可用于進一步改善偏移成像效果.就改進偏移成像效果而言,現(xiàn)行積分法偏移流程都包括一個后續(xù)的剩余動校及拉伸切除步驟;這一步驟將存在微小彎曲的同相軸拉平并切除無法拉平和存在拉伸畸變的部分.盡管剩余動校量較小,但它使得同相軸的波峰和波谷對得更齊,因而在疊加過程中更好地保持了高頻成分;而切除也提高了成像的信噪比.但對常規(guī)的炮域偏移而言,由于產(chǎn)出的是隨炮點變化的CIG道集,同相軸的彎曲沒有一定的規(guī)律,有效成像段并不一定在中心部分,將工業(yè)界經(jīng)常使用的剩余動校及拉伸切除應用于炮域波動方程偏移就產(chǎn)生了困難.而提取角道集后使得在炮域波動方程疊前深度偏移處理中應用剩余動校和拉伸切除變得與Kirchhoff積分偏移一樣容易.就本文提取的角道集而言,其剩余動校量可定量描述為式中z是零角度處同相軸的深度,γ是角度.在1.0附近對ρ進行掃描并將角道集沿角度疊加,可得到不同ρ值對應的能量譜,由能量譜中拾取最大能量對應的ρ,由式(12)計算不同角度處同相軸的移動量Δz,同時將大角度處的無效信息切除,這樣就簡單地實現(xiàn)了剩余動校及拉伸切除.5計算值的示例5.1理論炮集的構建為驗證本文發(fā)展的方法能夠得到相對保幅的角道集,文中設計了一個四層的水平層狀均勻速度模型,將層位按照從上到下排列,分別編號為1,2,3,4.各層位深度-厚度及速度密度信息見表1,其中第3層為一個低速薄層.先應用聲波方程正演方法生成理論炮集記錄圖3進一步給出在三個界面處,圖2c所示角道集反映的AVA特性與理論值的對比,可見兩者擬合得很好;大角度出現(xiàn)偏差是因為本文正演結果的長偏移距數(shù)據(jù)有限導致的.5.2角道集的疊加與同相軸的旋轉(zhuǎn)本節(jié)將文中提出的波動方程疊前深度偏移直接生成角道集流程應用于我國東部油田某地區(qū)滿覆蓋11.2km利用圖1所示流程,由該數(shù)據(jù)集3600個炮集偏移結果直接生成了滿覆蓋區(qū)各成像點的角道集.出于計算效率以及計算精度的考慮,流程中的炮集深度偏移及入射角計算中的波場延拓均采用了最優(yōu)分裂Fourier方法圖4a展示了某成像點初始的角道集,在1.0附近對ρ進行掃描并將角道集沿角度疊加,可得如圖4b所示的能量譜,由能量譜中拾取最大能量對應的ρ,由式(12)計算不同角度處同相軸的移動量Δz,可將圖4a所示的角道集轉(zhuǎn)換到圖4c,比較兩個圖可見同相軸變得更加平直.由圖4c可見,大角度處已無有效信息,因此可將它們?nèi)鐖D4d所示那樣切除.圖5展示了滿覆蓋區(qū)由剩余動校前后角道集疊加得到的數(shù)據(jù)體抽取的深度為750m的水平切片對比.比較可見在圖5a顯示的經(jīng)過剩余動校后的水平切片上,信噪比得到較大提高,斷層及斷點顯示更為清晰.我們將所生成的角道集與這一地區(qū)的17口測井數(shù)據(jù)進行了AVA特性對比.圖6給出了某一井點處所生成的角度域成像道集與典型層位的AVA特性曲線對比,其中理論曲線的計算參數(shù)中,速度與密度信息是由測井數(shù)據(jù)得到的,泊松比參數(shù)是由有關地層的經(jīng)驗參數(shù)給出的.圖6a為該點所有深度(自地面至深度為5000m)的角度域成像道集;圖6b與圖6c為圖中所示的深度處典型層位的AVA特性曲線與依據(jù)測井信息計算的理論曲線對比.考察可見,所生成的角度域道集與理論計算曲線變化趨勢吻合較好,從而表明所生成的角度域成像道集展現(xiàn)了良好的AVA特性.6實際應用效果本文給出了基于炮域波動方程疊前深度偏移結果直接生成角道集的方法.該方法可將炮域偏移結果直接轉(zhuǎn)換成反映界面AVA特征的角道集.文中也分別建議了校正因炮點覆蓋不均勻?qū)е碌某上穹嫡`差的方法和基于角道集改善波動方程深度偏移成像效果的方法.與現(xiàn)行各類提取角道集方法相比,本文方法計算效率高,內(nèi)存需求及存儲需求甚少,可靈活應用于二維和三維問題,更適用于工業(yè)規(guī)模資料.文中給出了該方法的實際應用流程和三維資料的應用實例.理論及實際數(shù)值算例表明,本文方法提取的角道集真實反映了界面的AVA特征,而改善偏移成像效果的方法提高了偏移剖面的信噪比.較大規(guī)模的實際三維資料應用表明,本