專利名稱：基于方位保真角度域成像的裂縫型油氣儲(chǔ)層地震探測方法
技術(shù)領(lǐng)域：
本發(fā)明涉及一種裂縫型油氣儲(chǔ)層地震探測方法，尤其是涉及基于方位保真角度域成像的裂縫型油氣儲(chǔ)層地震探測方法。
背景技術(shù)：
窄方位觀測不能完整地記錄三維地震波場及夾雜著的噪聲，降低了地震成像能力，也影響噪聲的壓制。為了提高復(fù)雜地質(zhì)體成像精度，改善裂縫系統(tǒng)的預(yù)測與描述，提高儲(chǔ)層及其含油氣性預(yù)測，寬方位三維地震采集技術(shù)越來越受到工業(yè)界的重視(凌云， 2007)。當(dāng)?shù)卣鸩▊鞑ソ?jīng)過定向排列的裂縫系統(tǒng)時(shí)，會(huì)誘導(dǎo)波場的方位各向異性，表現(xiàn)為速度或走時(shí)隨方位的變化，也可能引起振幅乃至衰減隨方位角的變化。在地震勘探中，利用正常時(shí)差(NMO)速度、走時(shí)、振幅等的方位各向異性分析，可以推斷裂縫系統(tǒng)的走向與發(fā)育密度(Lynn et al.，1996 ；Li, 1997 ；Ruger, 1998) 0當(dāng)前，無論是陸地還是海洋油氣勘探中，三維寬方位地震數(shù)據(jù)越來越多。然而傳統(tǒng)考慮方位的地震資料處理(包括偏移成像)存在諸多不足。比如，常規(guī)地震數(shù)據(jù)處理與解釋方法還不能充分地利用這些高代價(jià)獲得的方位信息豐富的波場數(shù)據(jù)(Manning et al.，2008)。通過地震偏移獲得的成像數(shù)據(jù)體被用于構(gòu)造解釋與層序分析、聲阻抗反演與屬性提取等過程中，但該數(shù)據(jù)體是不同傳播方向波場成像值的某種平均。雖然人們也意識(shí)到分方位成像有助于對(duì)復(fù)雜斷層、裂縫系統(tǒng)等的描述，但這些應(yīng)用通常都是以地面炮點(diǎn)-接收點(diǎn)的方位角(即采集方位)分選的疊前數(shù)據(jù)子集(俗稱限制方位數(shù)據(jù)體)為出發(fā)點(diǎn)，獨(dú)立進(jìn)行正常時(shí)差(NMO)校正或疊前偏移等，然后分析成像結(jié)果隨采集方位的變化特征。這種做法理論上假設(shè)各個(gè)限制方位子數(shù)據(jù)體在波場傳播過程中是獨(dú)立的，并認(rèn)為采集方位代表了波路徑的方位特征，不利于側(cè)向散射、平面外反射能量的歸位與聚焦，影響三維復(fù)雜地質(zhì)結(jié)構(gòu)、不規(guī)則界面、裂縫系統(tǒng)以及一些特殊散射體的成像與描述(Perez and Marfurt，2008)。在方位保真成像時(shí)，除了區(qū)分觀測地震數(shù)據(jù)采集方位的差異，還應(yīng)當(dāng)區(qū)分成像點(diǎn)處地震波局部傳播方位的差異。近年美國Paradigm地球物理公司根據(jù)國際上的角度域疊前深度偏移理論研究成果，正在研制與測試一種基于疊前深度偏移的全方位角度域成像軟件。不過，受該技術(shù)的復(fù)雜的理論與高昂的計(jì)算成本的制約，近期僅在其公司內(nèi)部試驗(yàn)性使用，尚未被業(yè)界認(rèn)可，也與我國油氣資源勘探需求不太吻合。

發(fā)明內(nèi)容
本發(fā)明的目的就是為了克服上述現(xiàn)有技術(shù)存在的缺陷而提供一種針對(duì)裂縫型油氣儲(chǔ)層地震探測的方位保真角度域成像數(shù)據(jù)處理方法。本發(fā)明的目的可以通過以下技術(shù)方案來實(shí)現(xiàn)一種針對(duì)裂縫型油氣儲(chǔ)層地震探測的方位保真角度域成像數(shù)據(jù)處理方法，其特征在于，包括以下步驟1)模擬地面三維地震，數(shù)據(jù)采集器采集地震波信號(hào)，按振幅相對(duì)保真處理流程進(jìn)行預(yù)處理；
2)輸入設(shè)備根據(jù)地質(zhì)任務(wù)確定局部目標(biāo)及其成像范圍，并對(duì)分方位成像數(shù)組和方位/角度域成像數(shù)組進(jìn)行初始化；3)數(shù)據(jù)采集器將步驟1)預(yù)處理后的地震波數(shù)據(jù)逐道輸入到處理器中，處理器根據(jù)中心點(diǎn)位置與偏移孔徑限制，確定該道地震波對(duì)應(yīng)的最大成像范圍，判斷該道地震波對(duì)應(yīng)的最大成像范圍是否在選定的局部目標(biāo)成像范圍內(nèi)，若為是，執(zhí)行步驟4)，若為否，輸入下一道；4)處理器對(duì)該道地震波進(jìn)行每個(gè)成像時(shí)間孔徑內(nèi)所有點(diǎn)的角度域成像處理，判斷該道地震波所有的成像時(shí)間的角度域成像處理是否完成，若為是，執(zhí)行步驟幻，若為否，繼續(xù)執(zhí)行步驟4)；5)處理器判斷所有的地震波數(shù)據(jù)是否都理處完成，若為是，執(zhí)行步驟6)，若為否， 返回步驟3)；6)輸出脈沖響應(yīng)疊加后的局部入射角度域成像結(jié)果，包括局部入射方位角域多方位成像數(shù)據(jù)體Imagein。_azimuth(X，y，Φ，τ)、入射方位與入射角域共成像道集ADCIGine (x，y, Φ, Y, τ )，同時(shí)也得到了局部照明角度域成像結(jié)果，包括局部照明方位角域多方位成像數(shù)據(jù)體Image,〃_— (x，φ, τ)、局部照明方位與傾角域共成像道集ADCIG, (χ, 乂 φ, τ)；7)輸出多方位成像數(shù)據(jù)體 Imageinc_azifflUth(x, y，Φ，τ )與Image^zi-(χ,少,外，
基于分方位成像數(shù)據(jù)體，多視角描述復(fù)雜地質(zhì)體，全面地分析沉積現(xiàn)象與沉積過程，并可提取隨方位變化的疊后屬性；8)輸出方位保真的入射方位與入射角域共成像道集ADCIGine (X，y，φ，Υ，τ )，可由ADCIGine (x，y，Φ，Y, τ)實(shí)施方位時(shí)差分析、振幅隨入射角與方位角變化AVAZ分析或反演，如果判斷不存方位振幅變化，通過進(jìn)一步疊加得到入射角域共成像點(diǎn)道集ADCIGin。(X， 1，Y, τ)，可由 ADCIGinc(X，y，Y, τ)實(shí)施三維 AVA 分析或反演，也可由 ADCIGinc (x，y，Φ, Y, τ)與ADCIGine (x，y，γ, τ)獲得部分角度疊加剖面，實(shí)施彈性阻抗EI反演，通過AVAZ 分析或反演可獲得裂縫走向與密度的參數(shù)信息，通過AVA或EI反演可獲得反映儲(chǔ)層及其含油氣性的屬性數(shù)據(jù)與圖表，通過顯示器顯示。所述的步驟4)中的對(duì)該道地震波進(jìn)行成像時(shí)間的角度域成像處理過程如下1)從角度域成像空間開始時(shí)刻，根據(jù)當(dāng)前成像時(shí)間(τ )的偏移孔徑范圍以及中心點(diǎn)坐標(biāo)，確定當(dāng)前成像時(shí)間成像點(diǎn)橫向分布區(qū)間；2)在局部一維各向同性或具有垂直對(duì)稱軸的橫向各向同性VTI介質(zhì)中，判定下行射線始終位于垂直平面內(nèi)，上行射線始終位于另一垂直平面，成像點(diǎn)橫向分布區(qū)間內(nèi)，針對(duì)每個(gè)成像點(diǎn)(X，1，τ)，根據(jù)非雙曲時(shí)差方程計(jì)算炮點(diǎn)-成像點(diǎn)單程射線走時(shí)％以及成像點(diǎn)-接收點(diǎn)單程射線走時(shí)、，總的雙程走時(shí)為t = ts+tr,再利用公式ρ = 計(jì)算入射與散射射線參數(shù)Ps與Py它們相當(dāng)于入射與散射慢度矢量的水平分量，根據(jù)上述共平面假設(shè)，計(jì)算入射射線X方向與y方向的水平慢度分量&與4，以及散射射線χ方向與y方向的水平慢度分量Av與Av，同時(shí)基于非雙曲時(shí)差方程計(jì)算入射和散射射線垂直慢度分量qs與l·，得到入射慢度矢與散射慢度矢量^ =OV^VA)‘3)根據(jù)矢量運(yùn)算規(guī)則，計(jì)算入射慢度矢量與散射慢度矢量的夾角θ，若假設(shè)y軸與正北方向平行，則按下式計(jì)算局部入射平面與正北方向所呈的方位角Φ ；4)入射慢度矢量與散射慢度矢量之和稱為照明慢度矢量，它可按下式計(jì)算二 (八+八，4+八，仏+1)，再根據(jù)矢量運(yùn)算規(guī)則，計(jì)算照明滿度矢量與垂向的夾角θ 以及照明矢量的方位角識(shí)；5)將輸入地震道對(duì)應(yīng)走時(shí)的振幅乘以幾何擴(kuò)散校正權(quán)系數(shù)后疊加到入射角度域成像空間(x，y，Φ，Y，1)，其中乂，又，τ為時(shí)間域成像坐標(biāo)，Φ為局部入射方位角，Y為局部入射角(Y = Θ/2)；將輸入地震道對(duì)應(yīng)走時(shí)的振幅乘以幾何擴(kuò)散校正權(quán)系數(shù)后疊加到照明角度域成像空間(X,L釣·9,Γ )，其中ρ為局部照明方位角，θ為局部照明傾角；6)為了實(shí)現(xiàn)多方位成像，將輸入地震道對(duì)應(yīng)走時(shí)的振幅乘以幾何擴(kuò)散校正權(quán)系數(shù)后疊加到局部入射方位角域成像空間(x，y，Φ，τ)；將輸入地震道對(duì)應(yīng)走時(shí)的振幅乘以幾何擴(kuò)散校正權(quán)系數(shù)后疊加到局部照明方位角域成像空間(Χ,γ,φ,τ )，完成當(dāng)前成像時(shí)間孔徑內(nèi)所有點(diǎn)的角度域成像后，進(jìn)入下一個(gè)成像時(shí)間，返回步驟1)。所述的步驟8)中的部分角度包括小角度、中角度與大角度。所述的步驟2)中的局部目標(biāo)及其成像范圍包括成像空間起止iLine號(hào)與xLine 號(hào)，起止自激自收時(shí)間范圍。所述的步驟1)中的預(yù)處理包括疊前去噪、靜校正、地表一致性校正。與現(xiàn)有技術(shù)相比，本發(fā)明具有以下優(yōu)點(diǎn)1)根據(jù)地震波局部方向特征實(shí)現(xiàn)的方位保真疊前時(shí)間偏移成像，考慮了目前針對(duì)裂縫系統(tǒng)的方位各向異性分析主要基于疊前時(shí)間偏移成果數(shù)據(jù)的現(xiàn)狀，以及疊前時(shí)間偏移相對(duì)于疊前深度偏移的成本優(yōu)勢，可以為復(fù)雜地質(zhì)體高精度成像與描述、面向裂縫儲(chǔ)層的方位各向異性分析與參數(shù)提取產(chǎn)生高保真的多方位成像數(shù)據(jù)體與方位共成像點(diǎn)道集；2)技術(shù)流程實(shí)用性強(qiáng)，國內(nèi)外尚無同類方法技術(shù)公開發(fā)表，它使得油氣勘探開發(fā)中能夠以可承受的成本投入提升地震波成像對(duì)碳?xì)浠衔飪?chǔ)層(尤其是裂縫系統(tǒng))的預(yù)測與描述能力，尤其是有助于生物礁、碳酸鹽巖溶洞、火成巖侵入體、鹽丘等地質(zhì)體的成像與描述，提高致密砂巖、碳酸鹽巖、頁巖、火成巖和變質(zhì)巖油氣藏中裂縫系統(tǒng)的預(yù)測與描述精度；3)遵循地震傳播規(guī)律實(shí)現(xiàn)方位保真成像，恢復(fù)儲(chǔ)層“原位”的波形(時(shí)差、振幅、頻率等)隨入射方位與照明方位的變化信息；4)獲得的分方位或多方位成像數(shù)據(jù)體支持復(fù)雜地質(zhì)體(如邊緣不規(guī)則地質(zhì)體、復(fù)雜斷裂與裂縫系統(tǒng)等)的多視角精細(xì)描述；5)為裂縫儲(chǔ)層預(yù)測與描述提供方位/角度域共成像點(diǎn)道集，有力地支撐方位時(shí)差分析、方位速度分析、方位振幅與衰減分析等；6)為三維AVA分析與反演提供入射角度域共成像點(diǎn)道集。


圖1為本發(fā)明的流程圖；圖2為本發(fā)明的硬件結(jié)構(gòu)示意圖；圖3為常規(guī)疊前偏移提供的成像切片圖；圖4為本發(fā)明提供的四個(gè)地面采集方位的成像切片圖；圖5為本發(fā)明提供的四個(gè)地下照明方位的成像切片圖；圖6為本發(fā)明提供的第1個(gè)照明方位角面元、入射角6-10°與16-20°成像剖面
6圖；圖7為本發(fā)明提供的第2個(gè)照明方位角面元、入射角6-10°與16_20°成像剖面圖；圖8為本發(fā)明提供的第3個(gè)照明方位角面元、入射角6-10°與16_20°成像剖面圖；圖9為本發(fā)明提供的第4個(gè)照明方位角面元、入射角6-10°與16_20°成像剖面圖·
具體實(shí)施例方式下面結(jié)合附圖和具體實(shí)施例對(duì)本發(fā)明進(jìn)行詳細(xì)說明。實(shí)施例針對(duì)復(fù)雜地質(zhì)體(如斷裂與裂縫系統(tǒng)、碳酸鹽巖生物礁與溶蝕孔洞、火成巖侵入體等特殊地質(zhì)體)的精細(xì)成像與描述、方位各向異性分析與裂縫儲(chǔ)層預(yù)測和描述的需要， 在方位與振幅相對(duì)保真前處理基礎(chǔ)上實(shí)施方位保真局部角度域疊前時(shí)間偏移成像，獲得多種常規(guī)方法無法獲得的成果數(shù)據(jù)體，有助于在解釋階段充分挖掘?qū)挿轿坏卣鹳Y料的潛力。實(shí)施角度域疊前時(shí)間偏移處理之前，在常規(guī)處理流程中對(duì)疊前地震資料做到振幅相對(duì)保真的噪音壓制、靜校正以及地表一致性校正等準(zhǔn)備工作。先由常規(guī)疊前時(shí)間偏移處理獲得合理的速度模型供偏移使用，并轉(zhuǎn)換生成相對(duì)準(zhǔn)確的層速度供角度計(jì)算使用。如圖 1、圖2所示，具體實(shí)施方式
如下1)模擬地面三維地震，數(shù)據(jù)采集器1采集地震波信號(hào)，按振幅相對(duì)保真處理流程進(jìn)行預(yù)處理；2)輸入設(shè)備4根據(jù)地質(zhì)任務(wù)確定局部目標(biāo)及其成像范圍，并對(duì)分方位成像數(shù)組和方位/角度域成像數(shù)組進(jìn)行初始化；3)數(shù)據(jù)采集器1將步驟1)預(yù)處理后的地震波數(shù)據(jù)逐道輸入到處理器2中，處理器 2根據(jù)中心點(diǎn)位置與偏移孔徑限制，確定該道地震波對(duì)應(yīng)的最大成像范圍，判斷該道地震波對(duì)應(yīng)的最大成像范圍是否在選定的局部目標(biāo)成像范圍內(nèi)，若為是，執(zhí)行步驟4)，若為否，輸入下一道；4)處理器2對(duì)該道地震波進(jìn)行每個(gè)成像時(shí)間孔徑內(nèi)所有點(diǎn)的角度域成像處理，判斷該道地震波所有的成像時(shí)間的角度域成像處理是否完成，若為是，執(zhí)行步驟幻，若為否， 繼續(xù)執(zhí)行步驟4)；5)處理器2判斷所有的地震波數(shù)據(jù)是否都理處完成，若為是，執(zhí)行步驟6)，若為否，返回步驟3)；6)輸出脈沖響應(yīng)疊加后的局部入射角度域成像結(jié)果，包括局部入射方位角域多方位成像數(shù)據(jù)體Imagein。_azimuth(X，y，Φ，τ)、入射方位與入射角域共成像道集ADCIGine (x，y, Φ, Y, τ )，同時(shí)也得到了局部照明角度域成像結(jié)果，包括局部照明方位角域多方位成像數(shù)據(jù)體Image,(χ,乂 φ, τ)、局部照明方位與傾角域共成像道集ADCIGw (χ,乂 φ, 3, τ)；7)輸出多方位成像數(shù)據(jù)體Image
inc-azimuth
(X，y，φ，τ )與Image,"
-azimuth
(χ,γ, φ, τ),
基于分方位成像數(shù)據(jù)體，多視角描述復(fù)雜地質(zhì)體，通過顯示器3顯示，全面地分析沉積現(xiàn)象與沉積過程，并可提取隨方位變化的疊后屬性(如圖3、圖4與圖5)；
8)輸出方位保真的入射方位與入射角域共成像道集ADCIGinc (X，y，φ，Y，τ )，可由ADCIGine (x，y，Φ，Y, τ)實(shí)施方位時(shí)差分析、振幅隨入射角與方位角變化AVAZ分析或反演，如果判斷不存方位振幅變化，通過進(jìn)一步疊加得到入射角域共成像點(diǎn)道集ADCIGin。(X， 1，Y, τ)，可由 ADCIGinc(X，y，Y, τ)實(shí)施三維 AVA 分析或反演，也可由 ADCIGinc (x，y，Φ, Y, τ)與ADCIGine (x，y，γ, τ)獲得部分角度疊加剖面，實(shí)施彈性阻抗EI反演，通過AVAZ 分析或反演可獲得裂縫走向與密度的參數(shù)信息，通過AVA或EI反演可獲得反映儲(chǔ)層及其含油氣性的屬性數(shù)據(jù)與圖表，通過顯示器3顯示(如圖6、圖7、圖8與圖9)。所述的步驟4)中的對(duì)該道地震波進(jìn)行成像時(shí)間的角度域成像處理過程如下1)從角度域成像空間開始時(shí)刻，根據(jù)當(dāng)前成像時(shí)間(τ )的偏移孔徑范圍以及中心點(diǎn)坐標(biāo)，確定當(dāng)前成像時(shí)間成像點(diǎn)橫向分布區(qū)間；2)在局部一維各向同性或具有垂直對(duì)稱軸的橫向各向同性VTI介質(zhì)中，判定下行射線始終位于垂直平面內(nèi)，上行射線始終位于另一垂直平面，成像點(diǎn)橫向分布區(qū)間內(nèi)，針對(duì)每個(gè)成像點(diǎn)(X，1，τ)，根據(jù)非雙曲時(shí)差方程計(jì)算炮點(diǎn)-成像點(diǎn)單程射線走時(shí)％以及成像點(diǎn)-接收點(diǎn)單程射線走時(shí)、，總的雙程走時(shí)為t = ts+tr,再利用公式P = 計(jì)算入射與散射射線參數(shù)Ps與Py它們相當(dāng)于入射與散射慢度矢量的水平分量，根據(jù)上述共平面假設(shè)，計(jì)算入射射線X方向與y方向的水平慢度分量與，以及散射射線X方向與y方向的水平慢度分，同時(shí)基于非雙曲時(shí)差方程計(jì)算入射和散射射線垂直慢度分量qs與l·，得到入射慢度矢量^ =(ZV^Vl)與散射慢度矢量、=(/VA，；3)根據(jù)矢量運(yùn)算規(guī)則，計(jì)算入射慢度矢量與散射慢度矢量的夾角θ，若假設(shè)y軸與正北方向平行，則按下式計(jì)算局部入射平面與正北方向所呈的方位角φ ；4)入射慢度矢量與散射慢度矢量之和稱為照明慢度矢量，它可按下式計(jì)算 λ, =(i\+/v;\+iVi + ^·)，再根據(jù)矢量運(yùn)算規(guī)則，計(jì)算照明滿度矢量與垂向的夾角θ 以及照明矢量的方位角識(shí)；5)將輸入地震道對(duì)應(yīng)走時(shí)的振幅乘以幾何擴(kuò)散校正權(quán)系數(shù)后疊加到入射角度域成像空間(x，y，Φ，Y，1)，其中乂，又，τ為時(shí)間域成像坐標(biāo)，Φ為局部入射方位角，Y為局部入射角(Y = Θ/2)；將輸入地震道對(duì)應(yīng)走時(shí)的振幅乘以幾何擴(kuò)散校正權(quán)系數(shù)后疊加到照明角度域成像空間外·9,r )，其中ρ為局部照明方位角，θ為局部照明傾角；6)為了實(shí)現(xiàn)多方位成像，將輸入地震道對(duì)應(yīng)走時(shí)的振幅乘以幾何擴(kuò)散校正權(quán)系數(shù)后疊加到局部入射方位角域成像空間(x，y，Φ，τ)；將輸入地震道對(duì)應(yīng)走時(shí)的振幅乘以幾何擴(kuò)散校正權(quán)系數(shù)后疊加到局部照明方位角域成像空間(Χ,ν,φ,τ )，完成當(dāng)前成像時(shí)間孔徑內(nèi)所有點(diǎn)的角度域成像后，進(jìn)入下一個(gè)成像時(shí)間，返回步驟1)。
權(quán)利要求
1.一種針對(duì)裂縫型油氣儲(chǔ)層地震探測的方位保真角度域成像數(shù)據(jù)處理方法，其特征在于，包括以下步驟1)模擬地面三維地震，數(shù)據(jù)采集器采集地震波信號(hào)，按振幅相對(duì)保真處理流程進(jìn)行預(yù)處理；2)輸入設(shè)備根據(jù)地質(zhì)任務(wù)確定局部目標(biāo)及其成像范圍，并對(duì)分方位成像數(shù)組和方位/ 角度域成像數(shù)組進(jìn)行初始化；3)數(shù)據(jù)采集器將步驟1)預(yù)處理后的地震波數(shù)據(jù)逐道輸入到處理器中，處理器根據(jù)中心點(diǎn)位置與偏移孔徑限制，確定該道地震波對(duì)應(yīng)的最大成像范圍，判斷該道地震波對(duì)應(yīng)的最大成像范圍是否在選定的局部目標(biāo)成像范圍內(nèi)，若為是，執(zhí)行步驟4)，若為否，輸入下一道；4)處理器對(duì)該道地震波進(jìn)行每個(gè)成像時(shí)間孔徑內(nèi)所有點(diǎn)的角度域成像處理，判斷該道地震波所有的成像時(shí)間的角度域成像處理是否完成，若為是，執(zhí)行步驟幻，若為否，繼續(xù)執(zhí)行步驟4)；5)處理器判斷所有的地震波數(shù)據(jù)是否都理處完成，若為是，執(zhí)行步驟6)，若為否，返回步驟3)；6)輸出脈沖響應(yīng)疊加后的局部入射角度域成像結(jié)果，包括局部入射方位角域多方位成像數(shù)據(jù)體Imagein?！甿uth(x，y，Φ, τ )、入射方位與入射角域共成像道集ADCIGinc (x，y，Φ, Y, τ),同時(shí)也得到了局部照明角度域成像結(jié)果，包括局部照明方位角域多方位成像數(shù)據(jù)體Image,,,,— (χ, φ, τ)、局部照明方位與傾角域共成像道集ADCIG,,, (x, φ, 3, τ)；7)輸出多方位成像數(shù)據(jù)體Imagein。_azimuth(X，y，Φ, τ)與 Image,〃_■_(>,}系 r)，基于分方位成像數(shù)據(jù)體，多視角描述復(fù)雜地質(zhì)體，全面地分析沉積現(xiàn)象與沉積過程，并可提取隨方位變化的疊后屬性；8)輸出方位保真的入射方位與入射角域共成像道集ADCIGin。(x，y，Φ,γ, τ)，可由 ADCIGinc (χ, y, Φ, γ, τ )實(shí)施方位時(shí)差分析、振幅隨入射角與方位角變化AVAZ分析或反演，如果判斷不存方位振幅變化，通過進(jìn)一步疊加得到入射角域共成像點(diǎn)道集ADCIGin。(X， 1，Y, τ)，可由 ADCIGinc(X，y，Y, τ)實(shí)施三維 AVA 分析或反演，也可由 ADCIGinc (x，y，Φ, Y, τ)與ADCIGine (x，y，γ, τ)獲得部分角度疊加剖面，實(shí)施彈性阻抗EI反演，通過AVAZ 分析或反演可獲得裂縫走向與密度的參數(shù)信息，通過AVA或EI反演可獲得反映儲(chǔ)層及其含油氣性的屬性數(shù)據(jù)與圖表，通過顯示器顯示。
2.根據(jù)權(quán)利要求1所述的一種針對(duì)裂縫型油氣儲(chǔ)層地震探測的方位保真角度域成像數(shù)據(jù)處理方法，其特征在于，所述的步驟4)中的對(duì)該道地震波進(jìn)行成像時(shí)間的角度域成像處理過程如下1)從角度域成像空間開始時(shí)刻，根據(jù)當(dāng)前成像時(shí)間(τ)的偏移孔徑范圍以及中心點(diǎn)坐標(biāo)，確定當(dāng)前成像時(shí)間成像點(diǎn)橫向分布區(qū)間；2)在局部一維各向同性或具有垂直對(duì)稱軸的橫向各向同性VTI介質(zhì)中，判定下行射線始終位于垂直平面內(nèi)，上行射線始終位于另一垂直平面，成像點(diǎn)橫向分布區(qū)間內(nèi)，針對(duì)每個(gè)成像點(diǎn)(x，y，τ )，根據(jù)非雙曲時(shí)差方程計(jì)算炮點(diǎn)-成像點(diǎn)單程射線走時(shí)％以及成像點(diǎn)-接收點(diǎn)單程射線走時(shí)、，總的雙程走時(shí)為t = ts+tr,再利用公式ρ = 計(jì)算入射與散射射線參數(shù)Ps與Py它們相當(dāng)于入射與散射慢度矢量的水平分量，根據(jù)上述共平面假設(shè)，計(jì)算入射射線χ方向與y方向的水平慢度分量&與八，以及散射射線χ方向與y方向的水平慢度分量，同時(shí)基于非雙曲時(shí)差方程計(jì)算入射和散射射線垂直慢度分量qs與l·，得到入射慢度矢量^ =(凡,，圪，仏)與散射慢度矢量、=OVjPvI)；3)根據(jù)矢量運(yùn)算規(guī)則，計(jì)算入射慢度矢量與散射慢度矢量的夾角θ，若假設(shè)y軸與正北方向平行，則按下式計(jì)算局部入射平面與正北方向所呈的方位角Φ ；4)入射慢度矢量與散射慢度矢量之和稱為照明慢度矢量，它可按下式計(jì)算 λ , =(/\+iVA+iV&+l)，再根據(jù)矢量運(yùn)算規(guī)則，計(jì)算照明滿度矢量與垂向的夾角θ 以及照明矢量的方位角爐；5)將輸入地震道對(duì)應(yīng)走時(shí)的振幅乘以幾何擴(kuò)散校正權(quán)系數(shù)后疊加到入射角度域成像空間(X，y，φ，Υ，1)，其中乂，7，τ為時(shí)間域成像坐標(biāo)，Φ為局部入射方位角，Y為局部入射角(Y = Θ/2)；將輸入地震道對(duì)應(yīng)走時(shí)的振幅乘以幾何擴(kuò)散校正權(quán)系數(shù)后疊加到照明角度域成像空間(X，乂約<9,r)，其中爐為局部照明方位角，θ為局部照明傾角；6)為了實(shí)現(xiàn)多方位成像，將輸入地震道對(duì)應(yīng)走時(shí)的振幅乘以幾何擴(kuò)散校正權(quán)系數(shù)后疊加到局部入射方位角域成像空間(x，y，Φ，τ)；將輸入地震道對(duì)應(yīng)走時(shí)的振幅乘以幾何擴(kuò)散校正權(quán)系數(shù)后疊加到局部照明方位角域成像空間(外r )，完成當(dāng)前成像時(shí)間孔徑內(nèi)所有點(diǎn)的角度域成像后，進(jìn)入下一個(gè)成像時(shí)間，返回步驟1)。
3.根據(jù)權(quán)利要求1所述的一種針對(duì)裂縫型油氣儲(chǔ)層地震探測的方位保真角度域成像數(shù)據(jù)處理方法，其特征在于，所述的步驟8)中的部分角度包括小角度、中角度與大角度。
4.根據(jù)權(quán)利要求1所述的一種針對(duì)裂縫型油氣儲(chǔ)層地震探測的方位保真角度域成像數(shù)據(jù)處理方法，其特征在于，所述的步驟2、中的局部目標(biāo)及其成像范圍包括成像空間起止 iLine號(hào)與xLine號(hào)，起止自激自收時(shí)間范圍。
5.根據(jù)權(quán)利要求1所述的一種針對(duì)裂縫型油氣儲(chǔ)層地震探測的方位保真角度域成像數(shù)據(jù)處理方法，其特征在于，所述的步驟1)中的預(yù)處理包括疊前去噪、靜校正、地表一致性校正。
全文摘要
本發(fā)明涉及一種針對(duì)裂縫型油氣儲(chǔ)層地震探測的方位保真角度域成像數(shù)據(jù)處理方法，包括以下步驟1)模擬地面三維地震，數(shù)據(jù)采集器采集地震波信號(hào)，按振幅相對(duì)保真處理流程進(jìn)行預(yù)處理；2)輸入設(shè)備根據(jù)地質(zhì)任務(wù)確定局部目標(biāo)及其成像范圍，并對(duì)分方位成像數(shù)組和方位/角度域成像數(shù)組進(jìn)行初始化；3)數(shù)據(jù)采集器將步驟1)預(yù)處理后的地震波數(shù)據(jù)逐道輸入到處理器中，處理器根據(jù)中心點(diǎn)位置與偏移孔徑限制，確定該道地震波對(duì)應(yīng)的最大成像范圍等。與現(xiàn)有技術(shù)相比，本發(fā)明具有復(fù)雜地質(zhì)體高精度成像與描述、面向裂縫儲(chǔ)層的方位各向異性分析與參數(shù)提取產(chǎn)生高保真的多方位成像數(shù)據(jù)體與方位共成像點(diǎn)道集等優(yōu)點(diǎn)。
文檔編號(hào)G01V1/36GK102540250SQ20101057949
公開日2012年7月4日 申請(qǐng)日期2010年12月8日 優(yōu)先權(quán)日2010年12月8日
發(fā)明者程玖兵, 馬在田 申請(qǐng)人:同濟(jì)大學(xué)