一種針對(duì)疊前地震數(shù)據(jù)共成像點(diǎn)道集的能量校正方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明屬于地球物理成像領(lǐng)域，具體涉及一種針對(duì)疊前地震數(shù)據(jù)共成像點(diǎn)道集的能量校正方法，針對(duì)kirchhoff時(shí)間偏移得到的共成像點(diǎn)道集實(shí)現(xiàn)保持AVA特征。
【背景技術(shù)】
[0002]隨著油氣勘探的不斷深入，油氣藏類型由簡(jiǎn)單的構(gòu)造油氣藏轉(zhuǎn)變?yōu)閺?fù)雜的巖性油氣藏，單一的疊后數(shù)據(jù)已經(jīng)不能滿足儲(chǔ)層預(yù)測(cè)的需要，基于具有豐富信息的疊前地震數(shù)據(jù)的反演技術(shù)(即AV0/AVA反演技術(shù))越來越得到人們的重視。AV0/AVA反演技術(shù)利用地震反射振幅在不同偏移距或者角度的相對(duì)改變特征量化反射界面上下介質(zhì)彈性參數(shù)的變化，其理論基礎(chǔ)是Zoeppritz方程和它的各種簡(jiǎn)化式。
[0003]疊前時(shí)間偏移后的共成像點(diǎn)道集基本消除了界面彎曲和傾斜地層對(duì)振幅的影響，并且從理論上講，反射波的振幅更接近于AV0/AVA分析的假設(shè)條件(反射層是水平的)，因此共成像點(diǎn)道集被廣泛的用于疊前儲(chǔ)層預(yù)測(cè)中。但是在實(shí)際地震資料應(yīng)用當(dāng)中，地震成像仍以構(gòu)造成像為主，很少為地震儲(chǔ)層描述而進(jìn)行保持動(dòng)力學(xué)特征的處理尤其是疊前彈性參數(shù)反演需要的共成像點(diǎn)道集數(shù)據(jù)，而Zoeppritz方程認(rèn)為有效信號(hào)在寬角度條件下隨著偏移距變化應(yīng)表現(xiàn)為拋物線特征，在共成像點(diǎn)道集中常常由于采集等因素影響到共成像點(diǎn)道集的能量因此無法得到這種拋物線特征。
[0004]目前對(duì)共成像點(diǎn)道集進(jìn)行能量補(bǔ)償?shù)某Ｒ?guī)做法大體分為2種:基于模型AVO特征的背景趨勢(shì)能量補(bǔ)償，該方法為模型驅(qū)動(dòng)，對(duì)一定時(shí)窗內(nèi)道集偏移距振幅能量變化進(jìn)行統(tǒng)計(jì)，通過模型AVO趨勢(shì)與實(shí)際AVO趨勢(shì)進(jìn)行比較，求取隨偏移距變化的刻度因子，最后將刻度因子用于實(shí)際道集，實(shí)現(xiàn)目的層背景趨勢(shì)意義下的偏移距振幅能量補(bǔ)償。但應(yīng)用該方法的前提是區(qū)域地質(zhì)特征明確對(duì)有利儲(chǔ)層區(qū)域和非有利儲(chǔ)層區(qū)域有一定的了解，但該條件往往無法滿足。另外一種是大時(shí)窗多道能量均衡方法，即根據(jù)CMP道集中覆蓋次數(shù)信息對(duì)疊前時(shí)間偏移后共成像點(diǎn)道集進(jìn)行基于覆蓋次數(shù)的能量調(diào)整。覆蓋次數(shù)通常指的是采集的時(shí)候，所有的炮對(duì)某一個(gè)共中心點(diǎn)的滿覆蓋次數(shù)，而KirchhofT偏移的原理是將每一個(gè)地震道對(duì)偏移孔徑內(nèi)的地下成像點(diǎn)的貢獻(xiàn)累加起來，簡(jiǎn)單的用覆蓋次數(shù)進(jìn)行校正，明顯不能滿足偏移孔徑內(nèi)的地震數(shù)據(jù)量變化復(fù)雜的要求。

【發(fā)明內(nèi)容】

[0005]本發(fā)明的目的在于解決上述現(xiàn)有技術(shù)中存在的難題，提供一種針對(duì)疊前地震數(shù)據(jù)共成像點(diǎn)道集的能量校正方法，針對(duì)由于地震采集的觀測(cè)系統(tǒng)導(dǎo)致的，疊前時(shí)間偏移得到的共成像點(diǎn)道集存在中入射角能量強(qiáng)、近遠(yuǎn)入射角能量弱的問題，通過在偏移的過程中統(tǒng)計(jì)照明次數(shù)，而不是簡(jiǎn)單的以覆蓋次數(shù)作為校正標(biāo)準(zhǔn)，對(duì)疊前時(shí)間偏移得到的共成像點(diǎn)道集進(jìn)行AVA特征恢復(fù)，旨在得到不受觀測(cè)系統(tǒng)影響的、保留真實(shí)AVA特征的疊前道集，利用該道集進(jìn)行后續(xù)的疊前彈性參數(shù)反演和AVA分析可以得到更準(zhǔn)確的結(jié)果。
[0006]本發(fā)明是通過以下技術(shù)方案實(shí)現(xiàn)的:
[0007]—種針對(duì)疊前地震數(shù)據(jù)共成像點(diǎn)道集的能量校正方法，通過統(tǒng)計(jì)每個(gè)成像點(diǎn)對(duì)應(yīng)的照明次數(shù)，對(duì)疊前時(shí)間偏移得到的共成像點(diǎn)道集以照明次數(shù)為基礎(chǔ)進(jìn)行保持AVA特征的能量校正。
[0008]所述方法包括以下步驟:
[0009](I)利用時(shí)間偏移速度對(duì)疊前未動(dòng)較共中心點(diǎn)道集進(jìn)行疊前時(shí)間偏移，得到偏移后的共成像點(diǎn)道集，同時(shí)統(tǒng)計(jì)每一個(gè)成像點(diǎn)處的照明次數(shù)，得到照明次數(shù)體；
[0010](2)對(duì)照明次數(shù)體設(shè)定最大閾值和最小閥值；
[0011 ] (3)將偏移后的共成像點(diǎn)道集中的每一成像點(diǎn)處的振幅值除以該成像點(diǎn)對(duì)應(yīng)的照明次數(shù)值，得到校正后的疊前時(shí)間偏移道集；
[0012](4)判斷所述校正后的疊前時(shí)間偏移道集的效果是否保持了 AVA特征，如果否，則修改照明次數(shù)體的最小閾值，然后返回步驟(3)，如果是，則所述校正后的疊前時(shí)間偏移道集為最后的校正結(jié)果。
[0013]所述步驟(I)中的照明次數(shù)是指在疊前偏移成像過程中根據(jù)地震射線路徑信息統(tǒng)計(jì)照射到成像點(diǎn)處各個(gè)照明角度面元的次數(shù)，即偏移孔徑內(nèi)所有地震數(shù)據(jù)對(duì)成像點(diǎn)處有一個(gè)次的貢獻(xiàn)，照明次數(shù)就增加一次；所述照明次數(shù)體與疊前時(shí)間偏移得到的偏移后的共成像點(diǎn)道集的格式、大小一樣，其每一個(gè)值代表的是與之對(duì)應(yīng)的成像點(diǎn)的照明次數(shù)。
[0014]所述步驟⑵具體如下:
[0015]篩選出所述照明次數(shù)體中照明次數(shù)的最大值和最小值，設(shè)置照明次數(shù)的最小閥值，最大值保留不變作為最大閾值，對(duì)于照明次數(shù)體中照明次數(shù)小于最小閾值的，則使其等于最小閾值。
[0016]所述步驟(4)中，如果最終得到的校正后道集中存在過校正的情況，則對(duì)存在過校正問題的道集進(jìn)行切除，切除后的道集AVA特征與共成像點(diǎn)道集對(duì)應(yīng)井資料計(jì)算得到目的層的AVA特征類似，可以驗(yàn)證切除后的道集上的AVA特征是正確的。
[0017]與現(xiàn)有技術(shù)相比，本發(fā)明的有益效果是:
[0018](I)本發(fā)明基于照明次數(shù)實(shí)現(xiàn)了疊前共成像點(diǎn)道集保持AVA特征的方法，對(duì)于疊前時(shí)間偏移得到的共成像點(diǎn)道集，由于地震采集的觀測(cè)系統(tǒng)導(dǎo)致的疊前時(shí)間偏移角度域道集存在中入射角能量強(qiáng)、近遠(yuǎn)入射角能量弱的問題，傳統(tǒng)的解決方法是利用覆蓋次數(shù)(通常指的是地震采集的時(shí)候，所有的炮對(duì)某一個(gè)共中心點(diǎn)的滿覆蓋次數(shù))對(duì)道集進(jìn)行一個(gè)近似的校正，而基于照明次數(shù)的校正是在疊前偏移成像過程中根據(jù)地震射線路徑信息統(tǒng)計(jì)照射到各個(gè)照明角度面元的次數(shù)，從理論上比簡(jiǎn)單的覆蓋次數(shù)校正更加合理、更加精確。
[0019](2)照明次數(shù)校正方法對(duì)所有疊前時(shí)間偏移的結(jié)果無論是偏移距道集或者是角度域道集都適用，并且具有易于實(shí)現(xiàn)、算法穩(wěn)定和計(jì)算效率高的特點(diǎn)，對(duì)于疊前海量數(shù)據(jù)不增加計(jì)算成本。
【附圖說明】
[0020]圖1本發(fā)明方法的步驟框圖。
[0021]圖2疊前未動(dòng)較共中心點(diǎn)道集。
[0022]圖3疊前時(shí)間偏移得到的共成像點(diǎn)道集。
[0023]圖4疊前時(shí)間偏移過程中統(tǒng)計(jì)的照明次數(shù)體。
[0024]圖5共成像點(diǎn)道集沿目的層提取的AVA特征(振幅隨入射角變化特征)。
[0025]圖6對(duì)照明次數(shù)體沿目的層提取照明次數(shù)隨入射角變化的分布特征。
[0026]圖7共成像點(diǎn)道集對(duì)應(yīng)井資料計(jì)算得到目的層的AVA特征。
[0027]圖8設(shè)置了最大、最小值閥值的照明次數(shù)體。
[0028]圖9經(jīng)過照明次數(shù)校正后的共成像點(diǎn)道集。
[0029]圖10校正后的共成像點(diǎn)道集沿目的層提取的AVA特征。
【具體實(shí)施方式】
[0030]下面結(jié)合附圖對(duì)本發(fā)明作進(jìn)一步詳細(xì)描述:
[0031]本發(fā)明用于對(duì)疊前時(shí)間偏移得到的共成像點(diǎn)道集進(jìn)行保持AVA特征的能量校正。其特征在于在kirchhoff時(shí)間偏移中植入代碼，統(tǒng)計(jì)每個(gè)成像點(diǎn)對(duì)應(yīng)的照明次數(shù)，對(duì)疊前時(shí)間偏移得到的共成像點(diǎn)道集以照明次數(shù)為基礎(chǔ)進(jìn)行保持AVA特征的能量校正:所述方法包括以下步驟:
[0032](I)常規(guī)kirchhoff時(shí)間偏移植入統(tǒng)計(jì)照明次數(shù)的代碼；
[0033](2)利用時(shí)間偏移速度對(duì)疊前未動(dòng)較共中心點(diǎn)道集進(jìn)行疊前時(shí)間偏移，得到偏移后的共成像點(diǎn)道集并統(tǒng)計(jì)每一個(gè)成像點(diǎn)處的照明次數(shù)，得到照明次數(shù)體；
[0034](3)對(duì)照明次數(shù)體設(shè)定最大最小閥值；
[0035](4)疊前時(shí)間偏移得到的共成像點(diǎn)道集中每一成像點(diǎn)處的振幅值除以對(duì)應(yīng)照明次數(shù)值，即共成像點(diǎn)道集與照明次數(shù)體做除法；
[0036](5)校正后的疊前時(shí)間偏移道集如果效果不理想則重復(fù)(3) (4)步驟，直到得到可靠的校正結(jié)果。
[0037]所述步驟(I)照明次數(shù)是指在疊前偏移成像過程中根據(jù)地震射線路徑信息統(tǒng)計(jì)照射到成像點(diǎn)處各個(gè)照明角度面元的次數(shù)，相當(dāng)于偏移孔徑內(nèi)所有地震數(shù)據(jù)對(duì)成像點(diǎn)處有一個(gè)次貢獻(xiàn)，照明次數(shù)就增加一次。在偏移循環(huán)過程中加入count參數(shù)，對(duì)于第i個(gè)成像點(diǎn)，令count[i] = O,在偏移成像的時(shí)候每次輸入道在角道集上實(shí)現(xiàn)一次疊加時(shí),使count [i]=count[i]+l，最終形