本發(fā)明涉及地震資料處理及成像，特別涉及一種針對(duì)陸地資料窄方位、低覆蓋采集地震數(shù)據(jù)的裂縫預(yù)測方法。

背景技術(shù)：

1、裂縫型油氣藏勘探技術(shù)對(duì)于石油工業(yè)的發(fā)展具有重要意義。裂縫的地質(zhì)成因復(fù)雜，其發(fā)育受到多方面的因素控制，具有分布復(fù)雜，規(guī)律性差，預(yù)測和識(shí)別難度大的特點(diǎn)。目前，基于高密度寬方位采集地震資料的疊前、疊后裂縫預(yù)測理論與技術(shù)方法已經(jīng)趨于完善，得到廣泛應(yīng)用。但是受到實(shí)際采集環(huán)境與勘探精度的制約，有些地區(qū)的資料是窄方位角，低覆蓋次數(shù)采集的，而針對(duì)窄方位采集的地震數(shù)據(jù)，無法實(shí)現(xiàn)橢圓擬合，增加了裂縫預(yù)測的難度。因此，需要一套完整的技術(shù)體系來解決窄方位、低覆蓋次數(shù)采集資料的裂縫預(yù)測問題。

技術(shù)實(shí)現(xiàn)思路

1、本發(fā)明所要解決的技術(shù)問題是克服背景技術(shù)中存在的現(xiàn)有方法針對(duì)窄方位采集的地震數(shù)據(jù)無法實(shí)現(xiàn)橢圓擬合、增加了裂縫預(yù)測的難度與預(yù)測結(jié)果的不確定性的問題，而提供一種針對(duì)陸地資料窄方位、低覆蓋采集地震數(shù)據(jù)的裂縫預(yù)測方法。該針對(duì)陸地資料窄方位、低覆蓋采集地震數(shù)據(jù)的裂縫預(yù)測方法，在現(xiàn)有采集資料的基礎(chǔ)上，能夠從本質(zhì)上提高裂縫預(yù)測的精度，包含的地震信息更加全面，為裂縫預(yù)測提供了更準(zhǔn)確的基礎(chǔ)數(shù)據(jù)。

2、本發(fā)明解決其問題可通過如下技術(shù)方案來達(dá)到：該針對(duì)陸地資料窄方位、低覆蓋采集地震數(shù)據(jù)的裂縫預(yù)測方法，包括以下步驟：

3、s1：建立工區(qū)時(shí)間域的浮動(dòng)基準(zhǔn)面模型；

4、s2：將時(shí)間域cmp道集校正至小平滑面；

5、s3：全方位局部角度域疊前深度偏移成像的數(shù)據(jù)準(zhǔn)備及加載；

6、s4：將主要偏移參數(shù)填寫到偏移成像系統(tǒng)；

7、s5：全方位局部角度域疊前深度偏移成像，生成傾角道集和反射角道集；

8、s6：對(duì)s5得到的反射角道集進(jìn)行方位各向異性校正處理，生成經(jīng)過方位各向異性校正處理的反射角道集；

9、s7：對(duì)方位各向異性校正后反射角道集的方位角組進(jìn)行劃分，得到劃分好方位角組的反射角道集數(shù)據(jù)：

10、s8：疊前裂縫預(yù)測。

11、進(jìn)一步的,所述s1建立工區(qū)時(shí)間域的浮動(dòng)基準(zhǔn)面模型的方法，包括:

12、將近地表高程值進(jìn)行平滑處理，建立小平滑面；

13、再應(yīng)用小平滑面建立深度域浮動(dòng)基準(zhǔn)面模型；

14、然后應(yīng)用對(duì)應(yīng)區(qū)塊的替換速度，進(jìn)行時(shí)深轉(zhuǎn)換，將深度域的浮動(dòng)基準(zhǔn)面模型轉(zhuǎn)換到時(shí)間域，得到時(shí)間域的浮動(dòng)基準(zhǔn)面模型。

15、進(jìn)一步的,所述s2將時(shí)間域cmp道集校正至小平滑面的方法，包括：

16、將時(shí)間域的疊前cmp道集加載到偏移成像系統(tǒng)；

17、然后應(yīng)用步驟s1所建立的時(shí)間域浮動(dòng)基準(zhǔn)面模型更新時(shí)間域cmp道集中的浮動(dòng)基準(zhǔn)面道頭信息，將時(shí)間域cmp道集校正至小平滑面。

18、進(jìn)一步的,所述s3全方位局部角度域疊前深度偏移成像的數(shù)據(jù)準(zhǔn)備及加載，包括：

19、將s1建立的小平滑面，作為全方位局部角度域疊前深度偏移起始面；

20、加載疊前深度偏移速度體、tti各向異性參數(shù)到偏移成像系統(tǒng)；

21、所述tti各向異性參數(shù)包括傾角體、方位角體、δ和ε。

22、進(jìn)一步的,所述s4的主要偏移參數(shù)包括偏移孔徑、反假頻參數(shù)、偏移步長為、淺層半開角以及深層半開角。

23、進(jìn)一步的,所述s5全方位局部角度域疊前深度偏移成像的具體方法，包括：

24、對(duì)已校正至小平滑面的cmp道集的每個(gè)繞射點(diǎn)進(jìn)行一維射線追蹤，形成一個(gè)扇狀射線區(qū)域；

25、由這個(gè)扇狀射線區(qū)域構(gòu)成射線對(duì)；

26、由所述射線對(duì)模擬反射和散射的同相軸，將每個(gè)射線對(duì)在菲涅耳帶范圍內(nèi)構(gòu)成射線束；

27、然后將射線束直徑內(nèi)的所有炮點(diǎn)和接收點(diǎn)的信息合成一個(gè)地震道去偏移成像，生成傾角道集和反射角道集。

28、進(jìn)一步的,所述s6對(duì)得到的反射角道集進(jìn)行方位各向異性校正處理的方法，包括以下步驟：

29、61、將全方位局部角度域疊前深度偏移得到的深度域反射角道集轉(zhuǎn)換至?xí)r間域；

30、62、對(duì)時(shí)間域反射角道集進(jìn)行切除和動(dòng)平衡處理；對(duì)時(shí)間域反射角道集進(jìn)行疊加、去噪處理后作為外部模型道；

31、輸入經(jīng)過切除和動(dòng)平衡處理后的時(shí)間域反射角道集以及外部模型道到地震資料處理系統(tǒng)，反射角道集內(nèi)的各道與輸入的外部模型道在定義時(shí)窗內(nèi)做互相關(guān)處理，根據(jù)互相關(guān)系數(shù)最大原則，得到每個(gè)互相關(guān)時(shí)窗點(diǎn)中的相對(duì)時(shí)移，即時(shí)窗中點(diǎn)的剩余時(shí)差，然后進(jìn)行內(nèi)插，得到整道數(shù)據(jù)的剩余時(shí)差和相關(guān)系數(shù)數(shù)據(jù)；

32、63、應(yīng)用上一步產(chǎn)生的相關(guān)系數(shù)數(shù)據(jù)和拾取的剩余時(shí)差數(shù)據(jù)，輸入到疊前時(shí)間偏移速度場作為參考速度場，進(jìn)行方位各向異性參數(shù)估算，生成方位各向異性參數(shù)；所述方位各向異性參數(shù)包括快波速度、慢波速度和慢波速度方位；

33、64、輸入反射角道集、參考速度場以及前面產(chǎn)生的方位各向異性參數(shù)到地震資料處理系統(tǒng)，進(jìn)行方位各向異性剩余時(shí)差校正處理，生成經(jīng)過方位各向異性校正處理的反射角道集。

34、進(jìn)一步的,所述步驟s7中對(duì)方位各向異性校正后反射角道集的方位角組進(jìn)行劃分的方案為：

35、71、將經(jīng)過方位各向異性校正處理的反射角道集輸入到裂縫預(yù)測系統(tǒng)中，進(jìn)行方位角組劃分試驗(yàn)；

36、72、按照保證覆蓋次數(shù)盡量一致的原則，當(dāng)一個(gè)方位角組的覆蓋次數(shù)高于其他方位角組時(shí)，對(duì)其范圍進(jìn)行調(diào)整，縮小方位角組的范圍；當(dāng)一個(gè)方位角組的覆蓋次數(shù)低于其他方位角組時(shí)，擴(kuò)大此方位角組的范圍，整個(gè)方位角組的劃分呈現(xiàn)扇形分布，而不是等分的對(duì)稱分布，從而得到劃分好方位角組的反射角道集數(shù)據(jù)。

37、進(jìn)一步的,所述步驟s8疊前裂縫預(yù)測，包括以下步驟：

38、81、將hti介質(zhì)模型輸入裂縫預(yù)測系統(tǒng)；

39、82、基于hti介質(zhì)模型，進(jìn)行方位屬性計(jì)算：通過應(yīng)用已經(jīng)劃分完方位角組的疊前方位角道集，進(jìn)行五個(gè)疊前方位屬性的計(jì)算，所述五個(gè)疊前方位屬性分別為：地震能量屬性，地震頻率屬性，指定頻率能量占總能量百分比，頻率衰減梯度及地震波阻抗屬性；

40、83、方位屬性分析及優(yōu)選：通過應(yīng)用區(qū)塊測井資料和巖石物理數(shù)據(jù)，建立關(guān)鍵井位置裂縫儲(chǔ)層的地質(zhì)模型和巖石物理模型，應(yīng)用全波動(dòng)模擬技術(shù)，計(jì)算地震波的地震響應(yīng)，分析由裂縫引起的地震屬性隨偏移距和方位角的變化特征，從而確定可用于裂縫預(yù)測的地震屬性體；

41、84、方位屬性橢圓擬合求取橢圓扁率：

42、進(jìn)一步的,所述方位屬性橢圓擬合求取橢圓扁率的方法，包括：

43、101、在hti介質(zhì)條件下的裂縫為高角度裂縫，裂縫傾角接近90度，當(dāng)縱波入射角較小時(shí)，各向異性介質(zhì)中縱波反射振幅與裂縫走向之間的關(guān)系可以表述為通過佐普里茲方程簡化成的aki-richards方程：

44、

45、式中：rpp為縱波反射振幅，g振幅隨偏移距的變化率，p為縱波垂直入射時(shí)的反射振幅(常量)；giso為振幅隨偏移距的變化率，即各向同性梯度；gani為振幅隨方位角的變化率，即各向異性梯度；為裂縫走向；θ為地震波的入射角；為炮點(diǎn)到檢波點(diǎn)的方位角；

46、102、通過對(duì)上述方程式(1)分析，在入射角θ為0°時(shí)，裂縫型地層的反射系數(shù)與方位角無關(guān)，說明近偏移距不能反映裂縫介質(zhì)的各向異性；隨著入射角的增大，各個(gè)方位角上反射系數(shù)相對(duì)變化增大，就越能觀測到地下介質(zhì)的各向異性特征；在方位角為0°，即垂直裂縫走向時(shí)，p波反射系數(shù)隨入射角變化最顯著，在方位角為90°，即平行裂縫走向時(shí)，p波反射系數(shù)隨入射角變化最??；

47、由此可以看出，不同入射角的地震波對(duì)裂縫的敏感度是不同的，選擇對(duì)裂縫敏感度較大的某個(gè)入射角范圍，通常選擇一個(gè)較小的入射角范圍，此范圍內(nèi)sin2θ變化不大，可視為常量c，上面的公式(1)可以簡化為：

48、

49、式中：a＝p+gisoc，b＝ganic；

50、103、在裂縫預(yù)測系統(tǒng)中，通過應(yīng)用不同方位角組的波阻抗屬性體，進(jìn)行平面擬合，可以得到一個(gè)橢圓，應(yīng)用公式(2)可以解出a和b；a表示橢圓短軸，b表示橢圓長軸，橢圓扁率b/a作為裂縫密度的相對(duì)度量，可以預(yù)測得到裂縫密度，當(dāng)b＝0時(shí)，介質(zhì)為各向同性介質(zhì)，即認(rèn)為不存在裂縫時(shí)，p波反射系數(shù)擬合為一半徑為a的圓形；b的絕對(duì)值越大，橢圓的形狀越扁，表示裂縫越發(fā)育，裂縫密度越大。

51、本發(fā)明與上述背景技術(shù)相比較可具有如下有益效果：該針對(duì)陸地資料窄方位、低覆蓋采集地震數(shù)據(jù)的裂縫預(yù)測方法，其核心在于在裂縫預(yù)測之前，對(duì)窄方位、低覆蓋次數(shù)采集的地震資料進(jìn)行全方位局部角度域疊前深度偏移處理，此種偏移方法從地下真實(shí)反射點(diǎn)出發(fā)，進(jìn)行射線的掃描和積分，是真正的全方位處理方法，相對(duì)于傳統(tǒng)的ovt域疊前時(shí)間偏移方法，成像的精度更高，包含的地震信息更加全面，為裂縫預(yù)測提供了更準(zhǔn)確的基礎(chǔ)數(shù)據(jù)。結(jié)合非等分方位角組劃分方式，能夠消除傳統(tǒng)等分方位角組劃分方式導(dǎo)致的覆蓋次數(shù)對(duì)能量的干擾，最終得到更加精確的裂縫密度預(yù)測結(jié)果。