地震數(shù)據(jù)的斷裂識別方法和裝置的制造方法
【專利摘要】本發(fā)明提供了一種地震數(shù)據(jù)的斷裂識別的方法及裝置，其中，該方法包括：獲取工區(qū)的地震數(shù)據(jù)；對工區(qū)的方位角進(jìn)行分組后，對各組方位角的地震數(shù)據(jù)分別進(jìn)行疊加，得到多組疊加地震數(shù)據(jù)；壓縮各組疊加地震數(shù)據(jù)的斷裂識別屬性后，對所得到的多組斷裂識別屬性進(jìn)行融合，得到地震數(shù)據(jù)的斷裂識別屬性；根據(jù)該結(jié)果，對工區(qū)進(jìn)行斷裂識別。在本發(fā)明實施例中，通過對不同類型的地震數(shù)據(jù)的斷裂識別屬性進(jìn)行壓縮，保留各種屬性對斷裂識別的共性，去除各種斷裂識別屬性中的異常信息，大大提高了地震屬性信噪比；利用融合技術(shù)，將不同方位的斷裂信息融合到一起，突出了每個方位的斷裂識別差異，提高斷裂識別精度以及識別能力，最終提高油氣勘探的效率。
【專利說明】
地震數(shù)據(jù)的斷裂識別方法和裝置
技術(shù)領(lǐng)域
[0001]本發(fā)明涉及石油勘探技術(shù)領(lǐng)域，特別涉及一種地震數(shù)據(jù)的斷裂識別方法和裝置。
【背景技術(shù)】
[0002]斷裂是一種普遍存在的較為復(fù)雜的地質(zhì)現(xiàn)象，對于油氣的運移、聚集起著很重要的控制作用，與油藏的形成、分布、富集有著十分密切的關(guān)系。因此，正確解釋斷裂、明確其分布與形態(tài)對于識別和描述油氣藏至關(guān)重要。
[0003]斷裂的存在會引起地震反射波同向軸錯斷、同向軸數(shù)目變化或者振幅相位突變等異?，F(xiàn)象，這些現(xiàn)象會導(dǎo)致地震屬性產(chǎn)生相應(yīng)的變化。常規(guī)斷裂識別的方法主要是在精細(xì)構(gòu)造解釋的基礎(chǔ)上，通過計算和分析疊前偏移數(shù)據(jù)體的相干體屬性和曲率屬性來進(jìn)行識別的。
[0004]但是，相干體屬性一般只能反映較大尺度的斷裂，曲率屬性雖然包含更多較小尺度的斷裂信息，但信噪比低、多解性強，很難得到隱藏在三維地震數(shù)據(jù)體中的微小斷裂信息。
[0005]針對上述僅采用相干體屬性或者曲率屬性進(jìn)行地震數(shù)據(jù)的斷裂識別時，很難得到隱藏在三維地震數(shù)據(jù)體中的微小斷裂信息的問題，目前尚未提出有效的解決方案。

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

[0006]本發(fā)明實施例提供了一種斷裂識別的方法和裝置，以解決現(xiàn)有技術(shù)中僅采用相干體屬性或者曲率屬性進(jìn)行地震數(shù)據(jù)的斷裂識別時，很難得到隱藏在三維地震數(shù)據(jù)體中的微小斷裂信息的問題。
[0007]本發(fā)明實施例提供了一種地震數(shù)據(jù)的斷裂識別方法，包括:步驟1:獲取工區(qū)的地震數(shù)據(jù);步驟2:對所述工區(qū)的方位角進(jìn)行分組，得到多組方位角范圍；步驟3:對所述多組方位角范圍內(nèi)的各組方位角范圍的地震數(shù)據(jù)分別進(jìn)行疊加，得到多組疊加地震數(shù)據(jù);步驟4:分別計算所述多組疊加地震數(shù)據(jù)中各組疊加地震數(shù)據(jù)的斷裂識別屬性，并對各組疊加地震數(shù)據(jù)的斷裂識別屬性進(jìn)行壓縮，得到多組壓縮后的斷裂識別屬性;步驟5:對所述多組壓縮后的斷裂識別屬性進(jìn)行融合，得到所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性;步驟6:根據(jù)所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性，對所述工區(qū)進(jìn)行斷裂識別。
[0008]在一個實施例中，根據(jù)所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性，對所述工區(qū)進(jìn)行斷裂識別之后，所述方法還包括:確定對所述工區(qū)的斷裂識別結(jié)果是否滿足預(yù)設(shè)要求;如果滿足所述預(yù)設(shè)要求，將當(dāng)前得到的斷裂識別結(jié)果作為所述工區(qū)的斷裂識別結(jié)果;如果不滿足所述預(yù)設(shè)要求，則重新選擇斷裂識別屬性，并重復(fù)執(zhí)行步驟4至步驟6，直至得到的所述工區(qū)的斷裂識別結(jié)果滿足所述預(yù)設(shè)要求。
[0009]在一個實施例中，根據(jù)所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性，對所述工區(qū)進(jìn)行斷裂識別之后，所述方法還包括:確定對所述工區(qū)的斷裂識別結(jié)果是否滿足預(yù)設(shè)要求;如果滿足所述預(yù)設(shè)要求，將當(dāng)前得到的斷裂識別結(jié)果作為所述工區(qū)的斷裂識別結(jié)果;如果不滿足所述預(yù)設(shè)要求，則重新對方位角進(jìn)行分組，并重復(fù)執(zhí)行步驟3至步驟6，直至得到的所述工區(qū)的斷裂識別結(jié)果滿足所述預(yù)設(shè)要求。
[0010]在一個實施例中，所述預(yù)設(shè)要求包括:所述工區(qū)的斷裂識別結(jié)果和預(yù)定的斷裂識別結(jié)果的差值小于預(yù)設(shè)閾值。
[0011]在一個實施例中，對各組疊加地震數(shù)據(jù)的斷裂識別屬性進(jìn)行壓縮，得到多組壓縮后的斷裂識別屬性，包括:通過核主成分分析方法對各組疊加地震數(shù)據(jù)的斷裂識別屬性進(jìn)行壓縮，得到多組壓縮后的斷裂識別屬性。
[0012]在一個實施例中，對所述多組壓縮后的斷裂識別屬性進(jìn)行融合，得到所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性，包括:確定所述多組壓縮后的斷裂識別屬性中各組壓縮后的斷裂識別屬性中表征斷裂特性的數(shù)值區(qū)域;剔除各組壓縮后的斷裂識別屬性中除表征斷裂特性的數(shù)值區(qū)域之外的數(shù)值區(qū)域;對所述各組壓縮后的斷裂識別屬性中表征斷裂特性的數(shù)值區(qū)域進(jìn)行疊加覆蓋，得到所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性。
[0013]在一個實施例中，所述斷裂識別屬性包括以下至少之一:相干體屬性、曲率屬性和能量梯度屬性。
[0014]在一個實施例中，在選取的斷裂識別屬性為曲率屬性的情況下，對各組疊加地震數(shù)據(jù)的斷裂識別屬性進(jìn)行壓縮，包括:對各組疊加地震數(shù)據(jù)的構(gòu)造曲率和振幅曲率進(jìn)行壓縮。
[0015]本發(fā)明實施例提供了一種地震數(shù)據(jù)的斷裂識別裝置，包括:獲取模塊，用于獲取工區(qū)的地震數(shù)據(jù);分組模塊，用于對所述工區(qū)的方位角進(jìn)行分組，得到多組方位角范圍；疊加模塊，用于對所述多組方位角范圍內(nèi)的各組方位角范圍的地震數(shù)據(jù)分別進(jìn)行疊加，得到多組疊加地震數(shù)據(jù);壓縮模塊，用于分別計算所述多組疊加地震數(shù)據(jù)中各組疊加地震數(shù)據(jù)的斷裂識別屬性，并對各組疊加地震數(shù)據(jù)的斷裂識別屬性進(jìn)行壓縮，得到多組壓縮后的斷裂識別屬性;融合模塊，用于對所述多組壓縮后的斷裂識別屬性進(jìn)行融合，得到所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性;識別模塊，用于根據(jù)所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性，對所述工區(qū)進(jìn)行斷裂識別。
[0016]在一個實施例中，所述識別模塊還包括:第一確定單元，用于確定對所述工區(qū)的斷裂識別結(jié)果是否滿足預(yù)設(shè)要求;第一識別單元，用于在所述工區(qū)的斷裂識別結(jié)果滿足所述預(yù)設(shè)要求的情況下，將當(dāng)前得到的斷裂識別結(jié)果作為所述工區(qū)的斷裂識別結(jié)果;第一控制單元，用于在所述工區(qū)的斷裂識別結(jié)果不滿足所述預(yù)設(shè)要求的情況下，則重新選擇斷裂識別屬性，并控制所述壓縮模塊，分別計算所述多組疊加地震數(shù)據(jù)中各組疊加地震數(shù)據(jù)重新選擇的斷裂識別屬性，并對各組疊加地震數(shù)據(jù)的斷裂識別屬性進(jìn)行壓縮，得到多組壓縮后的斷裂識別屬性，控制所述融合模塊對所述多組壓縮后的斷裂識別屬性進(jìn)行融合，得到所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性，控制所述識別模塊根據(jù)所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性，對所述工區(qū)進(jìn)行斷裂識別，直至得到的所述工區(qū)的斷裂識別結(jié)果滿足所述預(yù)設(shè)要求。
[0017]在本發(fā)明實施例中，在進(jìn)行地震數(shù)據(jù)的斷裂識別的過程中，首先，獲取工區(qū)的地震數(shù)據(jù);對所述工區(qū)的方位角進(jìn)行分組，得到多組方位角范圍；并對所述多組方位角范圍內(nèi)的各組方位角范圍的地震數(shù)據(jù)分別進(jìn)行疊加，得到多組疊加地震數(shù)據(jù);其次，分別計算所述多組疊加地震數(shù)據(jù)中各組疊加地震數(shù)據(jù)的斷裂識別屬性，并對各組疊加地震數(shù)據(jù)的斷裂識別屬性進(jìn)行壓縮，得到多組壓縮后的斷裂識別屬性;對所述多組壓縮后的斷裂識別屬性進(jìn)行融合，得到所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性;最后，根據(jù)所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性，對所述工區(qū)進(jìn)行斷裂識別。通過對不同類型的地震數(shù)據(jù)的斷裂識別屬性進(jìn)行壓縮，保留了各種屬性對斷裂識別的共性，去除了各種斷裂識別屬性中的異常信息，大大提高了地震屬性信噪比；其次，利用融合技術(shù)，將不同方位的斷裂信息融合到一起，突出了每個方位的斷裂識別差異，提高了斷裂識別精度以及識別能力，進(jìn)一步的提高了油氣勘探的效率。
【附圖說明】
[0018]此處所說明的附圖用來提供對本發(fā)明的進(jìn)一步理解，構(gòu)成本申請的一部分，并不構(gòu)成對本發(fā)明的限定。在附圖中:
[0019]圖1是本發(fā)明實施例的一種地震數(shù)據(jù)的斷裂識別方法的流程圖；
[0020]圖2是本發(fā)明實施例的高密度寬方位地震數(shù)據(jù)的斷裂識別方法的流程圖；
[0021 ]圖3是本發(fā)明實施例的目的層段應(yīng)用該發(fā)明方法后斷裂識別效果示意圖；
[0022]圖4是本發(fā)明實施例的地震數(shù)據(jù)的斷裂識別裝置的一種結(jié)構(gòu)框圖。
【具體實施方式】
[0023]為使本發(fā)明的目的、技術(shù)方案和優(yōu)點更加清楚明白，下面結(jié)合實施方式和附圖，對本發(fā)明做進(jìn)一步詳細(xì)說明。在此，本發(fā)明的示意性實施方式及其說明用于解釋本發(fā)明，但并不作為對本發(fā)明的限定。
[0024]考慮到目前地質(zhì)勘探中，對于找油找氣的目標(biāo)越來越精細(xì)，對于微小斷裂識別的精度要求越來越高，僅僅采用一種斷裂識別屬性進(jìn)行地震數(shù)據(jù)的斷裂識別時，很難得到隱藏在三維地震數(shù)據(jù)體中的微小斷裂信息的問題，而高密度寬方位地震數(shù)據(jù)能夠從不同方位上接收到不同走向斷裂系統(tǒng)的地震信息，因此發(fā)明人通過綜合利用不同方位以及不同地震屬性上與斷裂有關(guān)的信息來提高微小斷裂的識別精度。具體地，在本例中，提供了一種地震數(shù)據(jù)的斷裂識別方法，如圖1所示，可以包括以下步驟:
[0025]步驟101:獲取工區(qū)的地震數(shù)據(jù)；
[0026]在具體實現(xiàn)的時候，上述地震數(shù)據(jù)可以是疊后地震數(shù)據(jù)。其中，疊后是針對地震數(shù)據(jù)處理過程中的偏移來說的。當(dāng)?shù)叵陆橘|(zhì)的產(chǎn)狀不是水平狀的時候，反射地震的同相軸會發(fā)生偏移，此時的地震數(shù)據(jù)不能夠反映該處地下介質(zhì)的真實產(chǎn)狀，這就要求在資料處理的過程中進(jìn)行偏移歸位。通常，疊后地震數(shù)據(jù)指的是利用疊前時間偏移方法處理后的疊加數(shù)據(jù)。
[0027]步驟102:對所述工區(qū)的方位角進(jìn)行分組，得到多組方位角范圍；
[0028]方位角又稱地平經(jīng)度，是在平面上量度物體之間的角度差的方法之一，具體指的是從某點的指北方向線起，依順時針方向到目標(biāo)方向線之間的水平夾角。
[0029]在進(jìn)行高密度寬方位地震數(shù)據(jù)的斷裂識別時，可以根據(jù)工區(qū)內(nèi)區(qū)域及研究區(qū)應(yīng)力、斷裂發(fā)育情況和斷裂特征，對工區(qū)的方位角進(jìn)行分組，從而可以得到多組方位角范圍。
[0030]步驟103:對所述多組方位角范圍內(nèi)的各組方位角范圍的地震數(shù)據(jù)分別進(jìn)行疊加，得到多組疊加地震數(shù)據(jù)；
[0031]例如，工區(qū)的方位角為0°至180°，可以根據(jù)工區(qū)內(nèi)區(qū)域及研究區(qū)應(yīng)力、斷裂發(fā)育情況和斷裂特征，將工區(qū)的方位角分為0°至30°、30°至60°、60°至90°、90°至120°、120°至150°、150°至180°，再對工區(qū)內(nèi)進(jìn)行方位角分組后的六組地震數(shù)據(jù)進(jìn)行疊加，可以得到六組疊加地震數(shù)據(jù)。
[0032]步驟104:分別計算所述多組疊加地震數(shù)據(jù)中各組疊加地震數(shù)據(jù)的斷裂識別屬性，并對各組疊加地震數(shù)據(jù)的斷裂識別屬性進(jìn)行壓縮，得到多組壓縮后的斷裂識別屬性；
[0033]隨著高密度寬方位采集處理技術(shù)的不斷發(fā)展和完善，地震資料品質(zhì)得到了明顯的提高，無論是地震分辨率、振幅保真度及斷裂識別的能力均有了明顯的改善。一般情況下，不同工區(qū)發(fā)育多組斷裂，并且斷裂的走向以及傾向也不盡相同。根據(jù)地震波傳播理論，受各向異性影響，在平行或者垂直斷裂方向的地震數(shù)據(jù)上，P波地震屬性也會隨之變化。高密度寬方位地震數(shù)據(jù)能夠從不同方位上接收到不同走向斷裂系統(tǒng)的地震信息，因此可以利用這些地震屬性信息來提高斷裂的識別精度。
[0034]具體地，在本實施例中，地震數(shù)據(jù)的斷裂識別屬性可以包括但不限于以下至少之一:相干體屬性、曲率屬性以及能量梯度屬性等，其中，曲率屬性具體可以包括:構(gòu)造曲率以及振幅曲率等。
[0035]具體地，當(dāng)選取的斷裂識別屬性為曲率屬性的時候，可以對各組疊加地震數(shù)據(jù)的構(gòu)造曲率和振幅曲率進(jìn)行壓縮。然而值得注意的是，該具體實施例僅僅是為了更好地說明本發(fā)明，上述選取的斷裂識別屬性也可以是相干體屬性或者能量梯度屬性，相應(yīng)的，所對應(yīng)的進(jìn)行壓縮的相干體屬性、曲率屬性以及能量梯度屬性也不同，本發(fā)明對此不做限定。
[0036]在本發(fā)明實施例中，在計算得到工區(qū)內(nèi)多組疊加地震數(shù)據(jù)中各組疊加地震數(shù)據(jù)的斷裂識別屬性之后，可以通過核主成分分析方法對各組疊加地震數(shù)據(jù)的斷裂識別屬性進(jìn)行壓縮，并得到多組壓縮后的斷裂識別屬性。
[0037]核主成分分析是一種非線性推廣的主成分分析。由于在非線性統(tǒng)計和高階統(tǒng)計特征提取具有較好的效果，相對于傳統(tǒng)的主成分分析觀測數(shù)據(jù)的特點，具有更好的分類能力。
[0038]在本實施例中，可以采用核主成分分析的方法，通過引入非線性函數(shù)，將一定數(shù)量的斷裂識別屬性通過坐標(biāo)轉(zhuǎn)換的方式映射到高維空間。在這個新的高維坐標(biāo)系統(tǒng)下，變換數(shù)據(jù)點的方差，使得該方差可以沿著新的高維坐標(biāo)系統(tǒng)的坐標(biāo)軸達(dá)到最大化，以最有利于分類為原則，將同一方位的多個斷裂識別屬性壓縮成一種屬性，經(jīng)過該變換后得到的斷裂識別屬性維數(shù)降低、數(shù)量減少。
[0039]利用核主成分分析方法對各組疊加地震數(shù)據(jù)的斷裂識別屬性進(jìn)行壓縮，從而既保留了各種斷裂識別屬性對斷裂識別的共性，即每種屬性共同反映出來的有效的斷裂信息，還去除了相干體屬性、曲率屬性以及能量梯度屬性中所包含的異常信息，大大提高了地震屬性信噪比，有效地提高了斷裂的識別能力。
[0040]步驟105:對所述多組壓縮后的斷裂識別屬性進(jìn)行融合，得到所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性；
[0041]在得到多組壓縮后的斷裂識別屬性后，可以按照以下步驟對所述多組壓縮后的斷裂識別屬性進(jìn)行融合，得到所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性:
[0042]S1:確定所述多組壓縮后的斷裂識別屬性中各組壓縮后的斷裂識別屬性中表征斷裂特性的數(shù)值區(qū)域；
[0043]S2:剔除各組壓縮后的斷裂識別屬性中除表征斷裂特性的數(shù)值區(qū)域之外的數(shù)值區(qū)域；
[0044]S3:對所述各組壓縮后的斷裂識別屬性中表征斷裂特性的數(shù)值區(qū)域進(jìn)行疊加覆蓋，得到所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性。
[0045]步驟106:根據(jù)所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性，對所述工區(qū)進(jìn)行斷裂識別。
[0046]結(jié)合工區(qū)內(nèi)地質(zhì)研究成果，對以上多信息融合的斷裂識別屬性進(jìn)行斷裂識別，可以確定斷裂識別結(jié)果的可靠性與精度。
[0047]如果上述識別結(jié)果滿足預(yù)設(shè)要求，則可以將該識別結(jié)果作為工區(qū)的斷裂識別結(jié)果;如果上述識別結(jié)果不滿足預(yù)設(shè)要求，可以對方位角進(jìn)行重新分組，或者是采用其它的斷裂識別屬性進(jìn)行斷裂識別，直至識別得到的結(jié)果可以滿足精度要求。其中，上述預(yù)設(shè)要求指的是:工區(qū)的斷裂識別結(jié)果和預(yù)定的斷裂識別結(jié)果相比小于預(yù)定的閾值。具體地，上述預(yù)設(shè)要求可以包括以下幾個方面:
[0048](I)工區(qū)地震剖面的縱向上出現(xiàn)同向軸錯斷以及同向軸撓曲的現(xiàn)象、或者，振幅橫向有突變；
[0049](2)斷裂剖面樣式以及平面組合形態(tài)滿足預(yù)定情況，具體地，可通過分析工區(qū)的地質(zhì)背景和應(yīng)力場方向，然后根據(jù)斷裂形成機理進(jìn)行判斷；
[0050](3)存在斷裂，具體地，可以利用錄井資料和測井資料等鉆井信息驗證斷裂是否存在。
[0051]進(jìn)一步地，針對上述基于重新分組和重新選取斷裂識別屬性以提高識別精度的方法說明如下:
[0052]I)重新對方位角進(jìn)行分組:
[0053]S卩，采用與上一次識別過程中不同的分組結(jié)果，例如，上一次的分組結(jié)果為0°至30。、30°至60°、60°至90°、90°至120°、120°至150°、150°至180°等。那么本次可以將重新將方位角分為:15° 至45°、45° 至75°、75° 至 105°、105° 至 135°、135° 至 165°、165° 至 15° 等，也可以在保證覆蓋次數(shù)的前提下，將數(shù)據(jù)分成更多的方位。
[0054]2)采用不同的斷裂識別屬性
[0055]例如，上次選擇的斷裂識別屬性為相干體屬性，那么本次可以選擇曲率屬性作為斷裂識別屬性。
[0056]現(xiàn)舉一具體實例來進(jìn)一步地說明本發(fā)明，具體可以描述如下:該工區(qū)的地震數(shù)據(jù)為高密度寬方位采集地震數(shù)據(jù)，且根據(jù)區(qū)域斷裂展布規(guī)律，以水平方向為基礎(chǔ)，將工區(qū)內(nèi)的方位角劃分為6組，從而可以得到6組不同方位角范圍的地震數(shù)據(jù)。有這樣的資料基礎(chǔ)，就可以開展地震數(shù)據(jù)的斷裂識別。然而，值得注意的是，該實例僅僅是為了更好的說明本發(fā)明，并不構(gòu)成對本發(fā)明的不當(dāng)限定。如圖2所示，可以包括以下步驟:
[0057]步驟201:獲取工區(qū)的高密度寬方位地震數(shù)據(jù)；
[0058]步驟202:對工區(qū)的方位角進(jìn)行分組，可以得到6組方位角范圍；
[0059]步驟203:對6組方位角范圍內(nèi)的各組方位角范圍的地震數(shù)據(jù)分別進(jìn)行疊加，從而可以得到6組疊加地震數(shù)據(jù)；
[0060]步驟204:分別計算6組疊加地震數(shù)據(jù)中各組疊加地震數(shù)據(jù)的斷裂識別屬性，并對各組疊加地震數(shù)據(jù)的斷裂識別屬性進(jìn)行壓縮，可以得到6組壓縮后的斷裂識別屬性；
[0061 ]具體地，可以分別計算該工區(qū)內(nèi)6個方位的構(gòu)造曲率、振幅曲率等體屬性共計12套，平面屬性12張;并將工區(qū)內(nèi)同一方位的所有12套構(gòu)造曲率、振幅曲率等體屬性進(jìn)行屬性壓縮。
[0062]步驟205:對6組壓縮后的斷裂識別屬性進(jìn)行融合，得到所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性；
[0063]具體地，將6個方位的屬性進(jìn)行壓縮后所得到的6張平面屬性進(jìn)行融合，該融合結(jié)果綜合了各種屬性的有效信息。
[0064]步驟206:根據(jù)所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性，對所述工區(qū)進(jìn)行斷裂識別。
[0065]具體地，采用本發(fā)明實施例所提出的地震數(shù)據(jù)的斷裂識別方法對高密度寬方位采集地震數(shù)據(jù)進(jìn)行識別后可以得到如圖3所示的目的層段應(yīng)用該發(fā)明方法后斷裂識別效果示意圖。
[0066]基于同一發(fā)明構(gòu)思，本發(fā)明實施例中還提供了一種地震數(shù)據(jù)的斷裂識別裝置，如下面的實施例所述。由于地震數(shù)據(jù)的斷裂識別裝置解決問題的原理與地震數(shù)據(jù)的斷裂識別方法相似，因此地震數(shù)據(jù)的斷裂識別裝置的實施可以參見地震數(shù)據(jù)的斷裂識別方法的實施，重復(fù)之處不再贅述。以下所使用的，術(shù)語“單元”或者“模塊”可以實現(xiàn)預(yù)定功能的軟件和/或硬件的組合。盡管以下實施例所描述的裝置較佳地以軟件來實現(xiàn)，但是硬件，或者軟件和硬件的組合的實現(xiàn)也是可能并被構(gòu)想的。圖4是本發(fā)明實施例的地震數(shù)據(jù)的斷裂識別裝置的一種結(jié)構(gòu)框圖，如圖4所示，包括:獲取模塊401、分組模塊402、疊加模塊403、壓縮模塊404、融合模塊405、識別模塊406，下面對該結(jié)構(gòu)進(jìn)行說明。
[0067]獲取模塊401，用于獲取工區(qū)的地震數(shù)據(jù)；
[0068]分組模塊402，用于對所述工區(qū)的方位角進(jìn)行分組，得到多組方位角范圍；
[0069]疊加模塊403，用于對所述多組方位角范圍內(nèi)的各組方位角范圍的地震數(shù)據(jù)分別進(jìn)行疊加，得到多組疊加地震數(shù)據(jù)；
[0070]壓縮模塊404，用于分別計算所述多組疊加地震數(shù)據(jù)中各組疊加地震數(shù)據(jù)的斷裂識別屬性，并對各組疊加地震數(shù)據(jù)的斷裂識別屬性進(jìn)行壓縮，得到多組壓縮后的斷裂識別屬性；
[0071]融合模塊405，用于對所述多組壓縮后的斷裂識別屬性進(jìn)行融合，得到所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性；
[0072]識別模塊406，用于根據(jù)所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性，對所述工區(qū)進(jìn)行斷裂識別。
[0073]在一個實施例中，所述識別模塊還包括:第一確定單元，用于確定對所述工區(qū)的斷裂識別結(jié)果是否滿足預(yù)設(shè)要求;第一識別單元，用于在所述工區(qū)的斷裂識別結(jié)果滿足所述預(yù)設(shè)要求的情況下，將當(dāng)前得到的斷裂識別結(jié)果作為所述工區(qū)的斷裂識別結(jié)果;第一控制單元，用于在所述工區(qū)的斷裂識別結(jié)果不滿足所述預(yù)設(shè)要求的情況下，則重新選擇斷裂識別屬性，并控制所述壓縮模塊，分別計算所述多組疊加地震數(shù)據(jù)中各組疊加地震數(shù)據(jù)重新選擇的斷裂識別屬性，并對各組疊加地震數(shù)據(jù)的斷裂識別屬性進(jìn)行壓縮，得到多組壓縮后的斷裂識別屬性，控制所述融合模塊對所述多組壓縮后的斷裂識別屬性進(jìn)行融合，得到所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性，控制所述識別模塊根據(jù)所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性，對所述工區(qū)進(jìn)行斷裂識別，直至得到的所述工區(qū)的斷裂識別結(jié)果滿足所述預(yù)設(shè)要求。
[0074]在一個實施例中，所述識別模塊還包括:第二確定單元，用于確定對所述工區(qū)的斷裂識別結(jié)果是否滿足預(yù)設(shè)要求;第三識別單元，用于在所述工區(qū)的斷裂識別結(jié)果滿足所述預(yù)設(shè)要求的情況下，將當(dāng)前得到的斷裂識別結(jié)果作為所述工區(qū)的斷裂識別結(jié)果;第二控制單元，用于在所述工區(qū)的斷裂識別結(jié)果不滿足所述預(yù)設(shè)要求的情況下，則重新對方位角進(jìn)行分組，并控制所述疊加模塊對重新分組后所得到的多組方位角范圍內(nèi)的各組方位角范圍的地震數(shù)據(jù)分別進(jìn)行疊加，得到多組疊加地震數(shù)據(jù)，控制所述壓縮模塊分別計算所述多組疊加地震數(shù)據(jù)中各組疊加地震數(shù)據(jù)的斷裂識別屬性，并對各組疊加地震數(shù)據(jù)的斷裂識別屬性進(jìn)行壓縮，得到多組壓縮后的斷裂識別屬性，控制所述融合模塊對所述多組壓縮后的斷裂識別屬性進(jìn)行融合，得到所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性，控制所述識別模塊根據(jù)所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性，對所述工區(qū)進(jìn)行斷裂識別，直至得到的所述工區(qū)的斷裂識別結(jié)果滿足所述預(yù)設(shè)要求。
[0075]在一個實施例中，所述預(yù)設(shè)要求包括:所述工區(qū)的斷裂識別結(jié)果和預(yù)定的斷裂識別結(jié)果的差值小于預(yù)設(shè)閾值。
[0076]在一個實施例中，所述壓縮模塊包括:數(shù)據(jù)壓縮單元，用于通過核主成分分析方法對各組疊加地震數(shù)據(jù)的斷裂識別屬性進(jìn)行壓縮，得到多組壓縮后的斷裂識別屬性。
[0077]在一個實施例中，所述融合模塊包括:數(shù)值區(qū)域確定單元，用于確定所述多組壓縮后的斷裂識別屬性中各組壓縮后的斷裂識別屬性中表征斷裂特性的數(shù)值區(qū)域;數(shù)值區(qū)域剔除單元，用于剔除各組壓縮后的斷裂識別屬性中除表征斷裂特性的數(shù)值區(qū)域之外的數(shù)值區(qū)域;數(shù)據(jù)識別單元，用于對所述各組壓縮后的斷裂識別屬性中表征斷裂特性的數(shù)值區(qū)域進(jìn)行疊加覆蓋，得到所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性。
[0078]在一個實施例中，所述斷裂識別屬性包括以下至少之一:相干體屬性、曲率屬性、和能量梯度屬性。
[0079]在一個實施例中，當(dāng)所述壓縮模塊用于壓縮所述曲率屬性的時候，所述裝置包括:曲率壓縮單元，用于對各組疊加地震數(shù)據(jù)的構(gòu)造曲率和振幅曲率進(jìn)行壓縮。
[0080]從以上的描述中，可以看出，本發(fā)明實施例實現(xiàn)了如下技術(shù)效果:在進(jìn)行地震數(shù)據(jù)的斷裂識別的過程中，首先，獲取工區(qū)的地震數(shù)據(jù);對所述工區(qū)的方位角進(jìn)行分組，得到多組方位角范圍；并對所述多組方位角范圍內(nèi)的各組方位角范圍的地震數(shù)據(jù)分別進(jìn)行疊加，得到多組疊加地震數(shù)據(jù);其次，分別計算所述多組疊加地震數(shù)據(jù)中各組疊加地震數(shù)據(jù)的斷裂識別屬性，并對各組疊加地震數(shù)據(jù)的斷裂識別屬性進(jìn)行壓縮，得到多組壓縮后的斷裂識別屬性;對所述多組壓縮后的斷裂識別屬性進(jìn)行融合，得到所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性;最后，根據(jù)所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性，對所述工區(qū)進(jìn)行斷裂識別。通過對不同類型的地震數(shù)據(jù)的斷裂識別屬性進(jìn)行壓縮，保留了各種屬性對斷裂識別的共性，去除了各種斷裂識別屬性中的異常信息，大大提高了地震屬性信噪比；其次，利用融合技術(shù)，將不同方位的斷裂信息融合到一起，突出了每個方位的斷裂識別差異，提高了斷裂識別精度以及識別能力，進(jìn)一步的，提尚了油氣勘探的效率。
[0081]顯然，本領(lǐng)域的技術(shù)人員應(yīng)該明白，上述的本發(fā)明實施例的各模塊或各步驟可以用通用的計算裝置來實現(xiàn)，它們可以集中在單個的計算裝置上，或者分布在多個計算裝置所組成的網(wǎng)絡(luò)上，可選地，它們可以用計算裝置可執(zhí)行的程序代碼來實現(xiàn)，從而，可以將它們存儲在存儲裝置中由計算裝置來執(zhí)行，并且在某些情況下，可以以不同于此處的順序執(zhí)行所示出或描述的步驟，或者將它們分別制作成各個集成電路模塊，或者將它們中的多個模塊或步驟制作成單個集成電路模塊來實現(xiàn)。這樣，本發(fā)明實施例不限制于任何特定的硬件和軟件結(jié)合。
[0082]以上所述僅為本發(fā)明的優(yōu)選實施例而已，并不用于限制本發(fā)明，對于本領(lǐng)域的技術(shù)人員來說，本發(fā)明實施例可以有各種更改和變化。凡在本發(fā)明的精神和原則之內(nèi)，所作的任何修改、等同替換、改進(jìn)等，均應(yīng)包含在本發(fā)明的保護范圍之內(nèi)。
【主權(quán)項】
1.一種地震數(shù)據(jù)的斷裂識別方法，其特征在于，包括: 步驟1:獲取工區(qū)的地震數(shù)據(jù)； 步驟2:對所述工區(qū)的方位角進(jìn)行分組，得到多組方位角范圍； 步驟3:對所述多組方位角范圍內(nèi)的各組方位角范圍的地震數(shù)據(jù)分別進(jìn)行疊加，得到多組疊加地震數(shù)據(jù)； 步驟4:分別計算所述多組疊加地震數(shù)據(jù)中各組疊加地震數(shù)據(jù)的斷裂識別屬性，并對各組疊加地震數(shù)據(jù)的斷裂識別屬性進(jìn)行壓縮，得到多組壓縮后的斷裂識別屬性； 步驟5:對所述多組壓縮后的斷裂識別屬性進(jìn)行融合，得到所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性； 步驟6:根據(jù)所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性，對所述工區(qū)進(jìn)行斷裂識別。2.如權(quán)利要求1所述的方法，其特征在于，根據(jù)所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性，對所述工區(qū)進(jìn)行斷裂識別之后，所述方法還包括: 確定對所述工區(qū)的斷裂識別結(jié)果是否滿足預(yù)設(shè)要求； 如果滿足所述預(yù)設(shè)要求，將當(dāng)前得到的斷裂識別結(jié)果作為所述工區(qū)的斷裂識別結(jié)果； 如果不滿足所述預(yù)設(shè)要求，則重新選擇斷裂識別屬性，并重復(fù)執(zhí)行步驟4至步驟6，直至得到的所述工區(qū)的斷裂識別結(jié)果滿足所述預(yù)設(shè)要求。3.如權(quán)利要求1所述的方法，其特征在于，根據(jù)所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性，對所述工區(qū)進(jìn)行斷裂識別之后，所述方法還包括: 確定對所述工區(qū)的斷裂識別結(jié)果是否滿足預(yù)設(shè)要求； 如果滿足所述預(yù)設(shè)要求，將當(dāng)前得到的斷裂識別結(jié)果作為所述工區(qū)的斷裂識別結(jié)果； 如果不滿足所述預(yù)設(shè)要求，則重新對方位角進(jìn)行分組，并重復(fù)執(zhí)行步驟3至步驟6，直至得到的所述工區(qū)的斷裂識別結(jié)果滿足所述預(yù)設(shè)要求。4.如權(quán)利要求2至3中任一項所述的方法，其特征在于，所述預(yù)設(shè)要求包括: 所述工區(qū)的斷裂識別結(jié)果和預(yù)定的斷裂識別結(jié)果的差值小于預(yù)設(shè)閾值。5.如權(quán)利要求1所述的方法，其特征在于，對各組疊加地震數(shù)據(jù)的斷裂識別屬性進(jìn)行壓縮，得到多組壓縮后的斷裂識別屬性，包括: 通過核主成分分析方法對各組疊加地震數(shù)據(jù)的斷裂識別屬性進(jìn)行壓縮，得到多組壓縮后的斷裂識別屬性。6.如權(quán)利要求1所述的方法，其特征在于，對所述多組壓縮后的斷裂識別屬性進(jìn)行融合，得到所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性，包括: 確定所述多組壓縮后的斷裂識別屬性中各組壓縮后的斷裂識別屬性中表征斷裂特性的數(shù)值區(qū)域； 剔除各組壓縮后的斷裂識別屬性中除表征斷裂特性的數(shù)值區(qū)域之外的數(shù)值區(qū)域； 對所述各組壓縮后的斷裂識別屬性中表征斷裂特性的數(shù)值區(qū)域進(jìn)行疊加覆蓋，得到所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性。7.如權(quán)利要求1至6中任一項所述的方法，其特征在于，所述斷裂識別屬性包括以下至少之一:相干體屬性、曲率屬性和能量梯度屬性。8.如權(quán)利要求7所述的方法，其特征在于，在選取的斷裂識別屬性為曲率屬性的情況下，對各組疊加地震數(shù)據(jù)的斷裂識別屬性進(jìn)行壓縮，包括: 對各組疊加地震數(shù)據(jù)的構(gòu)造曲率和振幅曲率進(jìn)行壓縮。9.一種地震數(shù)據(jù)的斷裂識別裝置，其特征在于，包括: 獲取模塊，用于獲取工區(qū)的地震數(shù)據(jù)； 分組模塊，用于對所述工區(qū)的方位角進(jìn)行分組，得到多組方位角范圍； 疊加模塊，用于對所述多組方位角范圍內(nèi)的各組方位角范圍的地震數(shù)據(jù)分別進(jìn)行疊加，得到多組疊加地震數(shù)據(jù)； 壓縮模塊，用于分別計算所述多組疊加地震數(shù)據(jù)中各組疊加地震數(shù)據(jù)的斷裂識別屬性，并對各組疊加地震數(shù)據(jù)的斷裂識別屬性進(jìn)行壓縮，得到多組壓縮后的斷裂識別屬性；融合模塊，用于對所述多組壓縮后的斷裂識別屬性進(jìn)行融合，得到所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性； 識別模塊，用于根據(jù)所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性，對所述工區(qū)進(jìn)行斷裂識別。10.如權(quán)利要求9所述的裝置，其特征在于，所述識別模塊還包括: 第一確定單元，用于確定對所述工區(qū)的斷裂識別結(jié)果是否滿足預(yù)設(shè)要求； 第一識別單元，用于在所述工區(qū)的斷裂識別結(jié)果滿足所述預(yù)設(shè)要求的情況下，將當(dāng)前得到的斷裂識別結(jié)果作為所述工區(qū)的斷裂識別結(jié)果； 第一控制單元，用于在所述工區(qū)的斷裂識別結(jié)果不滿足所述預(yù)設(shè)要求的情況下，則重新選擇斷裂識別屬性，并控制所述壓縮模塊，分別計算所述多組疊加地震數(shù)據(jù)中各組疊加地震數(shù)據(jù)重新選擇的斷裂識別屬性，并對各組疊加地震數(shù)據(jù)的斷裂識別屬性進(jìn)行壓縮，得到多組壓縮后的斷裂識別屬性，控制所述融合模塊對所述多組壓縮后的斷裂識別屬性進(jìn)行融合，得到所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性，控制所述識別模塊根據(jù)所述工區(qū)的地震數(shù)據(jù)的斷裂識別屬性，對所述工區(qū)進(jìn)行斷裂識別，直至得到的所述工區(qū)的斷裂識別結(jié)果滿足所述預(yù)設(shè)要求。
【文檔編號】G01V1/30GK105954797SQ201610245138
【公開日】2016年9月21日
【申請日】2016年4月19日
【發(fā)明人】趙小輝, 何開泉, 于寶利, 曹小璐
【申請人】中國石油天然氣集團公司, 中國石油集團東方地球物理勘探有限責(zé)任公司