基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的方法及裝置的制造方法
【專利摘要】本發(fā)明提供了一種基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的方法及裝置，涉及地下裂縫空間展布預(yù)測(cè)技術(shù)領(lǐng)域。方法包括：獲取地震數(shù)據(jù)和測(cè)井?dāng)?shù)據(jù)；進(jìn)行井震標(biāo)定，并對(duì)地震數(shù)據(jù)進(jìn)行目的層構(gòu)造解釋，確定待預(yù)測(cè)的目的層的頂?shù)走吔?；?duì)地震數(shù)據(jù)進(jìn)行去噪處理，形成去噪后地震數(shù)據(jù)體；確定去噪后地震數(shù)據(jù)體的相干屬性和螞蟻?zhàn)粉檶傩裕粚⑽浵佔(zhàn)粉檶傩缘闹涤蜃兓秶Ｕ秊橄喔蓪傩缘闹涤蜃兓秶?；將值域變化范圍校正后的螞蟻?zhàn)粉檶傩耘c相干屬性進(jìn)行融合，形成裂縫預(yù)測(cè)敏感屬性體；根據(jù)裂縫預(yù)測(cè)敏感屬性體在空間范圍內(nèi)的數(shù)值大小預(yù)測(cè)頂?shù)走吔鐑?nèi)的地下裂縫空間展布形態(tài)。本發(fā)明可以實(shí)現(xiàn)利用地震數(shù)據(jù)同時(shí)準(zhǔn)確預(yù)測(cè)地下不同尺度裂縫的空間展布形態(tài)的效果。
【專利說(shuō)明】
基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的方法及裝置
技術(shù)領(lǐng)域
[0001] 本發(fā)明涉及地下裂縫空間展布預(yù)測(cè)技術(shù)領(lǐng)域，尤其涉及一種基于地震數(shù)據(jù)預(yù)測(cè)地 下裂縫空間展布的方法及裝置。
【背景技術(shù)】
[0002] 近年來(lái)，隨著油氣勘探程度的不斷提高，地下儲(chǔ)層類型也日趨多樣化。目前的地下 儲(chǔ)層類型有孔隙型、孔洞型，還出現(xiàn)了裂縫一孔隙型、裂縫一孔洞型和裂縫型儲(chǔ)層等，而裂 縫發(fā)育的儲(chǔ)層往往伴隨著油氣高產(chǎn)。當(dāng)前研究表明，地下裂縫不僅是滲流通道，也是重要的 儲(chǔ)集空間，是影響儲(chǔ)層發(fā)育的重要因素之一。因此，準(zhǔn)確預(yù)測(cè)地下裂縫空間展布是油氣勘探 的重要研究?jī)?nèi)容。
[0003] 由于裂縫空間展布規(guī)律性差，且地下構(gòu)造較為復(fù)雜，對(duì)其準(zhǔn)確預(yù)測(cè)難度較大。目前 常用的地震裂縫預(yù)測(cè)技術(shù)包括相干、曲率、螞蟻?zhàn)粉櫋Ⅲw積膨脹系數(shù)和疊前方位各向異性分 析等技術(shù)。其中，相干技術(shù)可以較好預(yù)測(cè)大中尺度的地下裂縫，但預(yù)測(cè)細(xì)小裂縫能力仍然較 差。而螞蟻?zhàn)粉櫦夹g(shù)可以較好刻畫細(xì)小裂縫的空間展布形態(tài)，但裂縫空間展布趨勢(shì)不清晰。 其余上述地震裂縫預(yù)測(cè)技術(shù)均各有優(yōu)缺點(diǎn)，然而均不能準(zhǔn)確預(yù)測(cè)地下裂縫的空間展布。當(dāng) 前，如何利用地震數(shù)據(jù)同時(shí)準(zhǔn)確預(yù)測(cè)地下不同尺度裂縫的空間展布形態(tài)成為一個(gè)亟待解決 的問(wèn)題。

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

[0004] 本發(fā)明的實(shí)施例提供一種基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的方法及裝置，以 實(shí)現(xiàn)利用地震數(shù)據(jù)同時(shí)準(zhǔn)確預(yù)測(cè)地下不同尺度裂縫的空間展布的形態(tài)的效果。
[0005] 為達(dá)到上述目的，本發(fā)明采用如下技術(shù)方案：
[0006] -種基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的方法，包括：
[0007] 獲取地震資料保幅處理后的偏移疊加地震數(shù)據(jù)，并獲取地震工區(qū)內(nèi)鉆穿目的層的 測(cè)井?dāng)?shù)據(jù)；
[0008] 根據(jù)所述地震數(shù)據(jù)和所述測(cè)井?dāng)?shù)據(jù)進(jìn)行井震標(biāo)定，并對(duì)所述地震數(shù)據(jù)進(jìn)行目的層 構(gòu)造解釋，確定待預(yù)測(cè)的目的層的頂?shù)走吔纾?br>[0009] 對(duì)所述地震數(shù)據(jù)進(jìn)行去噪處理，形成去噪后地震數(shù)據(jù)體；
[0010] 確定去噪后地震數(shù)據(jù)體的相干屬性；
[0011]確定去噪后地震數(shù)據(jù)體的螞蟻?zhàn)粉檶傩裕?br>[0012]將所述螞蟻?zhàn)粉檶傩缘闹涤蜃兓秶Ｕ秊樗鱿喔蓪傩缘闹涤蜃兓秶?br>[0013] 將值域變化范圍校正后的螞蟻?zhàn)粉檶傩耘c所述相干屬性進(jìn)行融合，形成裂縫預(yù)測(cè) 敏感屬性體；
[0014] 根據(jù)所述裂縫預(yù)測(cè)敏感屬性體在空間范圍內(nèi)的數(shù)值大小預(yù)測(cè)所述頂?shù)走吔鐑?nèi)的 地下裂縫空間展布形態(tài)。
[0015] 具體的，對(duì)所述地震數(shù)據(jù)進(jìn)行去噪處理，形成去噪后地震數(shù)據(jù)體，包括：
[0016] 根據(jù)中值濾波法對(duì)所述地震數(shù)據(jù)進(jìn)行去噪處理，形成去噪后地震數(shù)據(jù)體；
[0017] 所述根據(jù)中值濾波法對(duì)所述地震數(shù)據(jù)進(jìn)行去噪處理，形成去噪后地震數(shù)據(jù)體，包 括：
[0018] 根據(jù)預(yù)先設(shè)置的時(shí)窗在所述地震數(shù)據(jù)中選取2y+l個(gè)樣本點(diǎn)；
[0019] 根據(jù)公式：5了，，義，丨，…，St，…，5' , x=丨，2?….N:確定所述預(yù) 先設(shè)置的時(shí)窗中的2y+l個(gè)樣本點(diǎn)中第y+1個(gè)樣本點(diǎn)的值；
[0020] 將所述5T"確定為地震數(shù)據(jù)中的樣本點(diǎn)Sx的去噪處理結(jié)果。
[0021] 具體的，確定去噪后地震數(shù)據(jù)體的相干屬性，包括：
[0022] 獲取去噪后地震數(shù)據(jù)體在一預(yù)先設(shè)置的分析窗口內(nèi)的多道地震數(shù)據(jù)的多個(gè)采樣 點(diǎn)，形成地震數(shù)據(jù)矩陣D:
[0024] 其中，k為每道地震數(shù)據(jù)的采樣點(diǎn)個(gè)數(shù);n為所述預(yù)先設(shè)置的分析窗口內(nèi)的地震數(shù) 據(jù)的道數(shù)；
[0025] 確定地震數(shù)據(jù)矩陣D的協(xié)方差矩陣C: K
[0026] C = DrD = Y^kd[ k=[
[0027] 其中，<為地震數(shù)據(jù)矩陣D的第n行向量，劣=[0,;2，…，<4,]，丨 < 尺油<為< 的轉(zhuǎn)置，4 f ;
[0028] 確定地震數(shù)據(jù)矩陣D對(duì)應(yīng)的相干屬性E:
[0030] 其中，Ak為協(xié)方差矩陣C的第k個(gè)特征值;Amax為協(xié)方差矩陣C的最大特征值。
[0031] 具體的，確定去噪后地震數(shù)據(jù)體的螞蟻?zhàn)粉檶傩裕ǎ?br>[0032] 根據(jù)去噪后地震數(shù)據(jù)體確定多個(gè)地震數(shù)據(jù)節(jié)點(diǎn)；
? / e螞蟻g允許走的節(jié)點(diǎn)
[0033] 根據(jù)公式 螞蟻g允許走的節(jié)點(diǎn)
[0034] 確定螞蟻算法尋求路徑節(jié)點(diǎn)的轉(zhuǎn)移概率，并將所述轉(zhuǎn)移概率確定為所述螞蟻?zhàn)粉?屬性M;
[0035] 其中，i為螞蟻當(dāng)前所在地震數(shù)據(jù)節(jié)點(diǎn)位置；j為螞蟻行進(jìn)到下一地震數(shù)據(jù)節(jié)點(diǎn)位 置;T^U+n)為從t時(shí)刻開始經(jīng)過(guò)n次迀移后在地震數(shù)據(jù)節(jié)點(diǎn)i和地震數(shù)據(jù)節(jié)點(diǎn)j連線上殘留 的信息素;nu(t)為t時(shí)刻由地震數(shù)據(jù)節(jié)點(diǎn)i到地震數(shù)據(jù)節(jié)點(diǎn)j的期望程度;a、e分別表示信息 素和螞蟻行進(jìn)路線的相對(duì)重要程度；
[0036] 其中，
^為信息素的殘留程度；1-P為信息素的減少 程度；A 為一次循環(huán)中留在從地震數(shù)據(jù)節(jié)點(diǎn)i到地震數(shù)據(jù)節(jié)點(diǎn)j路徑上的信息素；Tlj(t) 為t時(shí)刻在地震數(shù)據(jù)節(jié)點(diǎn)i和地震數(shù)據(jù)節(jié)點(diǎn)j連線上殘留的信息素。
[0037] 具體的，將所述螞蟻?zhàn)粉檶傩缘闹涤蜃兓秶Ｕ秊樗鱿喔蓪傩缘闹涤蜃兓?圍，包括：
[0038]獲取螞蟻?zhàn)粉檶傩訫的各采樣點(diǎn)Mi;
[0039] 根據(jù)公式：
P1，2,…，況；確定值域變化范圍校正 后的螞蟻?zhàn)粉檶傩?，的各采樣點(diǎn)，i;其中，Mmin為螞蟻?zhàn)粉檶傩訫的各采樣點(diǎn)的最小值;M max 為螞蟻?zhàn)粉檶傩訫的各采樣點(diǎn)的最大值;Emin為相干屬性E的最小值。
[0040] 具體的，將值域變化范圍校正后的螞蟻?zhàn)粉檶傩耘c所述相干屬性進(jìn)行融合，形成 裂縫預(yù)測(cè)敏感屬性體，包括：
[0041 ]根據(jù)公式:Zii XEi，i = 1，2，…，N;確定裂縫預(yù)測(cè)敏感屬性體Zi;其中，Ei為相干 屬性E與各采樣點(diǎn)姐對(duì)應(yīng)位置的采樣點(diǎn)。
[0042] -種基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的裝置，包括：
[0043]數(shù)據(jù)獲取單元，用于獲取地震資料保幅處理后的偏移疊加地震數(shù)據(jù)，并獲取地震 工區(qū)內(nèi)鉆穿目的層的測(cè)井?dāng)?shù)據(jù)；
[0044]頂?shù)走吔绱_定單元，用于根據(jù)所述地震數(shù)據(jù)和所述測(cè)井?dāng)?shù)據(jù)進(jìn)行井震標(biāo)定，并對(duì) 所述地震數(shù)據(jù)進(jìn)行目的層構(gòu)造解釋，確定待預(yù)測(cè)的目的層的頂?shù)走吔纾?br>[0045]去噪處理單元，用于對(duì)所述地震數(shù)據(jù)進(jìn)行去噪處理，形成去噪后地震數(shù)據(jù)體；
[0046] 相干屬性確定單元，用于確定去噪后地震數(shù)據(jù)體的相干屬性；
[0047] 螞蟻?zhàn)粉檶傩源_定單元，用于確定去噪后地震數(shù)據(jù)體的螞蟻?zhàn)粉檶傩裕?br>[0048]值域變化范圍校正單元，用于將所述螞蟻?zhàn)粉檶傩缘闹涤蜃兓秶Ｕ秊樗鱿?干屬性的值域變化范圍；
[0049] 裂縫預(yù)測(cè)敏感屬性體生成單元，用于將值域變化范圍校正后的螞蟻?zhàn)粉檶傩耘c所 述相干屬性進(jìn)行融合，形成裂縫預(yù)測(cè)敏感屬性體；
[0050] 地下裂縫空間展布形態(tài)預(yù)測(cè)單元，用于根據(jù)所述裂縫預(yù)測(cè)敏感屬性體在空間范圍 內(nèi)的數(shù)值大小預(yù)測(cè)所述頂?shù)走吔鐑?nèi)的地下裂縫空間展布形態(tài)。
[0051]此外，所述去噪處理單元，具體用于：
[0052]根據(jù)中值濾波法對(duì)所述地震數(shù)據(jù)進(jìn)行去噪處理，形成去噪后地震數(shù)據(jù)體；
[0053] 所述去噪處理單元，具體還用于：
[0054]根據(jù)預(yù)先設(shè)置的時(shí)窗在所述地震數(shù)據(jù)中選取2y+l個(gè)樣本點(diǎn)；
[0055] 根據(jù)公式：5W=/船//"/<Sr ,,S、，丨，…,S、',…,丨,S、，)，x=l.2,…，N:確定所述預(yù) 先設(shè)置的時(shí)窗中的2y+l個(gè)樣本點(diǎn)中第y+1個(gè)樣本點(diǎn)的值；
[0056] 將所述確定為地震數(shù)據(jù)中的樣本點(diǎn)Sx的去噪處理結(jié)果。
[0057] 此外，所述相干屬性確定單元，具體用于：
[0058]獲取去噪后地震數(shù)據(jù)體在一預(yù)先設(shè)置的分析窗口內(nèi)的多道地震數(shù)據(jù)的多個(gè)采樣 點(diǎn)，形成地震數(shù)據(jù)矩陣D:
[0060] 其中，k為每道地震數(shù)據(jù)的采樣點(diǎn)個(gè)數(shù);n為所述預(yù)先設(shè)置的分析窗口內(nèi)的地震數(shù) 據(jù)的道數(shù)；
[0061] 確定地震數(shù)據(jù)矩陣D的協(xié)方差矩陣C:
[0063] 其中，< 為地震數(shù)據(jù)矩陣D的第n行向量，< …，丨 尺此為< 的轉(zhuǎn)置，4
[0064] 確定地震數(shù)據(jù)矩陣D對(duì)應(yīng)的相干屬性E:
[0066] 其中，Ak為協(xié)方差矩陣C的第k個(gè)特征值;Amax為協(xié)方差矩陣C的最大特征值。
[0067] 另外，所述螞蟻?zhàn)粉檶傩源_定單元，具體用于：
[0068] 根據(jù)去噪后地震數(shù)據(jù)體確定多個(gè)地震數(shù)據(jù)節(jié)點(diǎn)；
_/ e螞蟻g允許走的節(jié)點(diǎn)
[0069] 根據(jù)公式： /運(yùn)螞蟻g允許走的節(jié)點(diǎn)
[0070] 確定螞蟻算法尋求路徑節(jié)點(diǎn)的轉(zhuǎn)移概率，并將所述轉(zhuǎn)移概率確定為所述螞蟻?zhàn)粉?屬性M;
[0071 ]其中，i為螞蟻當(dāng)前所在地震數(shù)據(jù)節(jié)點(diǎn)位置；j為螞蟻行進(jìn)到下一地震數(shù)據(jù)節(jié)點(diǎn)位 置;T^U+n)為從t時(shí)刻開始經(jīng)過(guò)n次迀移后在地震數(shù)據(jù)節(jié)點(diǎn)i和地震數(shù)據(jù)節(jié)點(diǎn)j連線上殘留 的信息素;n^(t)為t時(shí)刻由地震數(shù)據(jù)節(jié)點(diǎn)i到地震數(shù)據(jù)節(jié)點(diǎn)j的期望程度;a、0分別表示信息 素和螞蟻行進(jìn)路線的相對(duì)重要程度； r,7 (t-rn)-= p zu (/) + (1 - /.V) A r..
[0072]其中，j ^ . k 4為信息素的殘留程度；1-P為信息素的減少 Ar.. = 2,Ar.. 、 .k =1 u 程度；A 為一次循環(huán)中留在從地震數(shù)據(jù)節(jié)點(diǎn)i到地震數(shù)據(jù)節(jié)點(diǎn)j路徑上的信息素；Tlj(t) 為t時(shí)刻在地震數(shù)據(jù)節(jié)點(diǎn)i和地震數(shù)據(jù)節(jié)點(diǎn)j連線上殘留的信息素。
[0073]另外，所述值域變化范圍校正單元，具體用于：
[0074]獲取螞蟻?zhàn)粉檶傩訫的各采樣點(diǎn)Mi;
[0075] 根據(jù)公式：
i=l，2，???，況;確定值域變化范圍校正 后的螞蟻?zhàn)粉檶傩?，的各采樣點(diǎn);其中，Mmin為螞蟻?zhàn)粉檶傩訫的各采樣點(diǎn)的最小值;Mmax 為螞蟻?zhàn)粉檶傩訫的各采樣點(diǎn)的最大值;Emin為相干屬性E的最小值。
[0076] 另外，所述裂縫預(yù)測(cè)敏感屬性體生成單元，具體用于：
[0077] 根據(jù)公式:Zii XEi，i = 1，2，…，N;確定裂縫預(yù)測(cè)敏感屬性體Zi;其中，Ei為相干 屬性E與各采樣點(diǎn)姐對(duì)應(yīng)位置的采樣點(diǎn)。
[0078] 本發(fā)明實(shí)施例提供的基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的方法及裝置，首先， 獲取了地震資料保幅處理后的偏移疊加地震數(shù)據(jù)，并獲取了地震工區(qū)內(nèi)鉆穿目的層的測(cè)井 數(shù)據(jù);之后，根據(jù)所述地震數(shù)據(jù)和所述測(cè)井?dāng)?shù)據(jù)進(jìn)行井震標(biāo)定，并對(duì)所述地震數(shù)據(jù)進(jìn)行目的 層構(gòu)造解釋，確定待預(yù)測(cè)的目的層的頂?shù)走吔?之后，對(duì)所述地震數(shù)據(jù)進(jìn)行去噪處理，形成 去噪后地震數(shù)據(jù)體;確定去噪后地震數(shù)據(jù)體的相干屬性;確定去噪后地震數(shù)據(jù)體的螞蟻?zhàn)?蹤屬性;將所述螞蟻?zhàn)粉檶傩缘闹涤蜃兓秶Ｕ秊樗鱿喔蓪傩缘闹涤蜃兓秶粚⒅?域變化范圍校正后的螞蟻?zhàn)粉檶傩耘c所述相干屬性進(jìn)行融合，形成裂縫預(yù)測(cè)敏感屬性體； 根據(jù)所述裂縫預(yù)測(cè)敏感屬性體在空間范圍內(nèi)的數(shù)值大小預(yù)測(cè)所述頂?shù)走吔鐑?nèi)的地下裂縫 空間展布形態(tài)?？梢姡景l(fā)明可以將地震相干屬性和地震螞蟻?zhàn)粉檶傩赃M(jìn)行融合，在保留裂 縫空間展布趨勢(shì)的同時(shí)提高了地下裂縫預(yù)測(cè)的精度，從而可以實(shí)現(xiàn)利用地震數(shù)據(jù)同時(shí)準(zhǔn)確 預(yù)測(cè)地下不同尺度裂縫的空間展布的形態(tài)的效果。
【附圖說(shuō)明】
[0079] 為了更清楚地說(shuō)明本發(fā)明實(shí)施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對(duì)實(shí)施例或現(xiàn) 有技術(shù)描述中所需要使用的附圖作簡(jiǎn)單地介紹，顯而易見地，下面描述中的附圖僅僅是本 發(fā)明的一些實(shí)施例，對(duì)于本領(lǐng)域普通技術(shù)人員來(lái)講，在不付出創(chuàng)造性勞動(dòng)性的前提下，還可 以根據(jù)這些附圖獲得其他的附圖。
[0080] 圖1為本發(fā)明實(shí)施例提供的一種基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的方法的流 程圖一；
[0081] 圖2為本發(fā)明實(shí)施例提供的一種基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的方法的流 程圖二；
[0082] 圖3為本發(fā)明實(shí)施例中的原始地震剖面結(jié)構(gòu)示意圖；
[0083] 圖4為基于相干算法得到的裂縫預(yù)測(cè)剖面示意圖；
[0084] 圖5為基于相干算法得到的Z段地層裂縫預(yù)測(cè)結(jié)果示意圖；
[0085] 圖6為基于螞蟻?zhàn)粉櫵惴ǖ玫降牧芽p預(yù)測(cè)剖面示意圖；
[0086] 圖7為基于螞蟻?zhàn)粉櫵惴ǖ玫降腪段地層裂縫預(yù)測(cè)結(jié)果示意圖；
[0087] 圖8為基于本發(fā)明實(shí)施例的方法得到的裂縫預(yù)測(cè)剖面示意圖；
[0088] 圖9為基于本發(fā)明方法得到的Z段地層裂縫預(yù)測(cè)結(jié)果示意圖；
[0089] 圖10為本發(fā)明實(shí)施例提供的一種基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的裝置的 結(jié)構(gòu)示意圖。
【具體實(shí)施方式】
[0090] 下面將結(jié)合本發(fā)明實(shí)施例中的附圖，對(duì)本發(fā)明實(shí)施例中的技術(shù)方案進(jìn)行清楚、完 整地描述，顯然，所描述的實(shí)施例僅僅是本發(fā)明一部分實(shí)施例，而不是全部的實(shí)施例。基于 本發(fā)明中的實(shí)施例，本領(lǐng)域普通技術(shù)人員在沒(méi)有做出創(chuàng)造性勞動(dòng)前提下所獲得的所有其他 實(shí)施例，都屬于本發(fā)明保護(hù)的范圍。
[0091] 如圖1所示，本發(fā)明實(shí)施例提供一種基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的方法， 包括：
[0092] 步驟101、獲取地震資料保幅處理后的偏移疊加地震數(shù)據(jù)，并獲取地震工區(qū)內(nèi)鉆穿 目的層的測(cè)井?dāng)?shù)據(jù)。
[0093] 步驟102、根據(jù)所述地震數(shù)據(jù)和所述測(cè)井?dāng)?shù)據(jù)進(jìn)行井震標(biāo)定，并對(duì)所述地震數(shù)據(jù)進(jìn) 行目的層構(gòu)造解釋，確定待預(yù)測(cè)的目的層的頂?shù)走吔纭?br>[0094]步驟103、對(duì)所述地震數(shù)據(jù)進(jìn)行去噪處理，形成去噪后地震數(shù)據(jù)體。
[0095] 步驟104、確定去噪后地震數(shù)據(jù)體的相干屬性。
[0096] 步驟105、確定去噪后地震數(shù)據(jù)體的螞蟻?zhàn)粉檶傩浴?br>[0097] 步驟106、將所述螞蟻?zhàn)粉檶傩缘闹涤蜃兓秶Ｕ秊樗鱿喔蓪傩缘闹涤蜃兓?范圍。
[0098] 步驟107、將值域變化范圍校正后的螞蟻?zhàn)粉檶傩耘c所述相干屬性進(jìn)行融合，形成 裂縫預(yù)測(cè)敏感屬性體。
[0099] 步驟108、根據(jù)所述裂縫預(yù)測(cè)敏感屬性體在空間范圍內(nèi)的數(shù)值大小預(yù)測(cè)所述頂?shù)?邊界內(nèi)的地下裂縫空間展布形態(tài)。
[0100] 本發(fā)明實(shí)施例提供的基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的方法，首先，獲取了 地震資料保幅處理后的偏移疊加地震數(shù)據(jù)，并獲取了地震工區(qū)內(nèi)鉆穿目的層的測(cè)井?dāng)?shù)據(jù)； 之后，根據(jù)所述地震數(shù)據(jù)和所述測(cè)井?dāng)?shù)據(jù)進(jìn)行井震標(biāo)定，并對(duì)所述地震數(shù)據(jù)進(jìn)行目的層構(gòu) 造解釋，確定待預(yù)測(cè)的目的層的頂?shù)走吔?之后，對(duì)所述地震數(shù)據(jù)進(jìn)行去噪處理，形成去噪 后地震數(shù)據(jù)體;確定去噪后地震數(shù)據(jù)體的相干屬性;確定去噪后地震數(shù)據(jù)體的螞蟻?zhàn)粉檶?性;將所述螞蟻?zhàn)粉檶傩缘闹涤蜃兓秶Ｕ秊樗鱿喔蓪傩缘闹涤蜃兓秶?；將值域?化范圍校正后的螞蟻?zhàn)粉檶傩耘c所述相干屬性進(jìn)行融合，形成裂縫預(yù)測(cè)敏感屬性體;根據(jù) 所述裂縫預(yù)測(cè)敏感屬性體在空間范圍內(nèi)的數(shù)值大小預(yù)測(cè)所述頂?shù)走吔鐑?nèi)的地下裂縫空間 展布形態(tài)?？梢姡景l(fā)明可以將地震相干屬性和地震螞蟻?zhàn)粉檶傩赃M(jìn)行融合，在保留裂縫空 間展布趨勢(shì)的同時(shí)提高了地下裂縫預(yù)測(cè)的精度，從而可以實(shí)現(xiàn)利用地震數(shù)據(jù)同時(shí)準(zhǔn)確預(yù)測(cè) 地下不同尺度裂縫的空間展布的形態(tài)的效果。
[0101] 為了使本領(lǐng)域的技術(shù)人員更好的了解本發(fā)明，下面列舉一個(gè)更為詳細(xì)的實(shí)施例， 如圖2所示，本發(fā)明實(shí)施例提供一種基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的方法，包括：
[0102] 步驟201、獲取地震資料保幅處理后的偏移疊加地震數(shù)據(jù)，并獲取地震工區(qū)內(nèi)鉆穿 目的層的測(cè)井?dāng)?shù)據(jù)。
[0103] 此處根據(jù)該地震數(shù)據(jù)形成的原始地震剖面結(jié)構(gòu)可以如圖3所示。
[0104] 步驟202、根據(jù)所述地震數(shù)據(jù)和所述測(cè)井?dāng)?shù)據(jù)進(jìn)行井震標(biāo)定，并對(duì)所述地震數(shù)據(jù)進(jìn) 行目的層構(gòu)造解釋，確定待預(yù)測(cè)的目的層的頂?shù)走吔纭?br>[0105]此處的井震標(biāo)定是為了建立測(cè)井?dāng)?shù)據(jù)與地震數(shù)據(jù)之間的聯(lián)系，落實(shí)測(cè)井測(cè)量深度 對(duì)應(yīng)于地震時(shí)間域深度的位置。而進(jìn)行地震數(shù)據(jù)目的層構(gòu)造解釋的作用，主要是落實(shí)目的 層段頂?shù)椎目臻g位置，用于最后對(duì)目的層段裂縫空間展布形態(tài)描述設(shè)定頂?shù)走吔纭?br>[0106] 步驟203、根據(jù)中值濾波法對(duì)所述地震數(shù)據(jù)進(jìn)行去噪處理，形成去噪后地震數(shù)據(jù) 體。
[0107] 此處根據(jù)中值濾波法對(duì)所述地震數(shù)據(jù)進(jìn)行去噪處理，可以采用如下方式：
[0108] 根據(jù)預(yù)先設(shè)置的時(shí)窗在所述地震數(shù)據(jù)中選取2y+l個(gè)樣本點(diǎn)。
[0109] 根據(jù)公式：57? 收//"《(、. 叫,…,S V，……S、，，j, x=丨上…凡 先設(shè)置的時(shí)窗中的2y+l個(gè)樣本點(diǎn)中第y+1個(gè)樣本點(diǎn)的值。
[0110] 將所述Sf"確定為地震數(shù)據(jù)中的樣本點(diǎn)sx的去噪處理結(jié)果。
[0111] 根據(jù)上述步驟203可以去除地震數(shù)據(jù)中包含的隨機(jī)噪音，為提高相干屬性和螞蟻 追蹤屬性的裂縫預(yù)測(cè)精度奠定了基礎(chǔ)。
[0112] 步驟204、獲取去噪后地震數(shù)據(jù)體在一預(yù)先設(shè)置的分析窗口內(nèi)的多道地震數(shù)據(jù)的 多個(gè)采樣點(diǎn)，形成地震數(shù)據(jù)矩陣D。 ^11 ^12 '* ^ln
[0113] D= d：1 "； k\ dk2 …d^ 一
[0114] 其中，k為每道地震數(shù)據(jù)的采樣點(diǎn)個(gè)數(shù);n為所述預(yù)先設(shè)置的分析窗口內(nèi)的地震數(shù) 據(jù)的道數(shù)。
[0115] 步驟205、確定地震數(shù)據(jù)矩陣D的協(xié)方差矩陣C。 K
[0116] C - D' D ~ Zdkd[ k.=l.
[0117] 其中，<為地震數(shù)據(jù)矩陣D的第n行向量，t// =|^H,…,丨5； A'S A： ;dk為< 的轉(zhuǎn)置，dk=[dki，dk2,…，dkn]T。
[0118] 步驟206、確定地震數(shù)據(jù)矩陣D對(duì)應(yīng)的相干屬性E。
[0120] 其中，Ak為協(xié)方差矩陣C的第k個(gè)特征值;Amax為協(xié)方差矩陣C的最大特征值。
[0121] 步驟207、根據(jù)去噪后地震數(shù)據(jù)體確定多個(gè)地震數(shù)據(jù)節(jié)點(diǎn)。
[0122] 步驟208、確定螞蟻算法尋求路徑節(jié)點(diǎn)的轉(zhuǎn)移概率，并將所述轉(zhuǎn)移概率確定為所述 螞蟻?zhàn)粉檶傩訫。
[0123] 在螞蟻算法中，路徑的選擇依賴于信息素濃度，信息素濃度越大路徑被選擇的概 率越大，螞蟻根據(jù)轉(zhuǎn)移概率尋求路徑節(jié)點(diǎn)。
[0124] 此處，具體可以根據(jù)公式：
螞蟻g允許走的節(jié)點(diǎn) /運(yùn)螞蟻g允許走的節(jié)點(diǎn)
[0126] 確定螞蟻算法尋求路徑節(jié)點(diǎn)的轉(zhuǎn)移概率，并將所述轉(zhuǎn)移概率確定為所述螞蟻?zhàn)粉?屬性M。
[0127] 其中，i為螞蟻當(dāng)前所在地震數(shù)據(jù)節(jié)點(diǎn)位置；j為螞蟻行進(jìn)到下一地震數(shù)據(jù)節(jié)點(diǎn)位 置;T^U+n)為從t時(shí)刻開始經(jīng)過(guò)n次迀移后在地震數(shù)據(jù)節(jié)點(diǎn)i和地震數(shù)據(jù)節(jié)點(diǎn)j連線上殘留 的信息素，該T^(t+n)說(shuō)明，螞蟻經(jīng)過(guò)后路徑上的信息素會(huì)隨之增加，但隨著時(shí)間增加，信息 素也會(huì)減少;n^(t)為t時(shí)刻由地震數(shù)據(jù)節(jié)點(diǎn)i到地震數(shù)據(jù)節(jié)點(diǎn)j的期望程度;a、0分別表示信 息素和螞蟻行進(jìn)路線的相對(duì)重要程度。
[0128] 其中，^ . k 4為信息素的殘留程度；1-P為信息素的減少 {_ k^x u 程度；A 為一次循環(huán)中留在從地震數(shù)據(jù)節(jié)點(diǎn)i到地震數(shù)據(jù)節(jié)點(diǎn)j路徑上的信息素；Tlj(t) 為t時(shí)刻在地震數(shù)據(jù)節(jié)點(diǎn)i和地震數(shù)據(jù)節(jié)點(diǎn)j連線上殘留的信息素。
[0129] 值得說(shuō)明的是，對(duì)于螞蟻算法，每只螞蟻經(jīng)過(guò)n次迀移后就可得到一條路徑，如果 該路徑滿足最短路徑要求即停止，否則進(jìn)行重新搜索直至收斂于最優(yōu)路徑。當(dāng)螞蟻從某個(gè) 節(jié)點(diǎn)出發(fā)將所有可能連續(xù)轉(zhuǎn)移的位置均嘗試后，通過(guò)轉(zhuǎn)移概率值的大小即可確定出最短路 徑。
[0130]步驟209、獲取螞蟻?zhàn)粉檶傩訫的各采樣點(diǎn)Mi。
[0131 ]步驟210、確定值域變化范圍校正后的螞蟻?zhàn)粉檶傩訵的各采樣點(diǎn)，i。
[0132] 此處，可以根據(jù)公式：
K，L…，TV;確定值域變 化范圍校正后的螞蟻?zhàn)粉檶傩?，的各采樣點(diǎn)，i;其中，Mmin*螞蟻?zhàn)粉檶傩訫的各采樣點(diǎn)的 最小值;MmaA螞蟻?zhàn)粉檶傩訫的各采樣點(diǎn)的最大值;Emin為相干屬性E的最小值。
[0133] 此處由于螞蟻?zhàn)粉檶傩耘c相干屬性的值域變化范圍不同，因此需要將螞蟻?zhàn)粉檶?性的值域變化范圍校正到相干屬性的值域變化范圍。一般情況下，相干屬性的最大值為1， 表示完全無(wú)相關(guān)性，最小值介于〇至1之間。
[0134] 步驟211、將值域變化范圍校正后的螞蟻?zhàn)粉檶傩耘c所述相干屬性進(jìn)行融合，形成 裂縫預(yù)測(cè)敏感屬性體。
[0135] 具體的，此處可以根據(jù)公式:21 = 11/以￡1，1 = 1，2，~，1確定裂縫預(yù)測(cè)敏感屬性體 Z1;其中，E:為相干屬性E與各采樣點(diǎn)姐對(duì)應(yīng)位置的采樣點(diǎn)。
[0136] 步驟212、根據(jù)所述裂縫預(yù)測(cè)敏感屬性體在空間范圍內(nèi)的數(shù)值大小預(yù)測(cè)所述頂?shù)?邊界內(nèi)的地下裂縫空間展布形態(tài)。
[0137] 此處，該裂縫預(yù)測(cè)敏感屬性體在空間范圍內(nèi)的數(shù)值范圍為0至1，若該數(shù)值越接近 1，則表示裂縫越不發(fā)育；而若該數(shù)值越接近0,則表示裂縫越發(fā)育完好。具體可以參見后續(xù) 圖8、9所示。純黑色表示數(shù)值為0，純白色表示數(shù)值為1。
[0138] 為了使本發(fā)明實(shí)施例提供的基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的方法的效果 更加明確，下面結(jié)合現(xiàn)有技術(shù)和本發(fā)明實(shí)施例的效果圖進(jìn)行說(shuō)明。此處以四川盆地Z段地層 裂縫預(yù)測(cè)為例進(jìn)行說(shuō)明，該Z段地層巖性為白云巖，脆性較大，容易產(chǎn)生裂縫。
[0139] 圖4為基于相干算法得到的裂縫預(yù)測(cè)剖面；圖5為基于相干算法得到的Z段地層裂 縫預(yù)測(cè)結(jié)果。圖6為基于螞蟻?zhàn)粉櫵惴ǖ玫降牧芽p預(yù)測(cè)剖面；圖7為基于螞蟻?zhàn)粉櫵惴ǖ玫?的Z段地層裂縫預(yù)測(cè)結(jié)果；圖8為基于本發(fā)明實(shí)施例的方法得到的裂縫預(yù)測(cè)剖面;圖9為基于 本發(fā)明方法得到的Z段地層裂縫預(yù)測(cè)結(jié)果。可見，圖4、5裂縫空間展布規(guī)律性較好，但是局部 細(xì)節(jié)連續(xù)性較差。圖6、7預(yù)測(cè)的裂縫空間展布散亂無(wú)規(guī)律。圖8、9與上述圖4、5相比不僅保持 裂縫原有的空間展布規(guī)律性，而且在局部細(xì)節(jié)更加連續(xù)清晰，描述裂縫更為準(zhǔn)確。
[0140] 本發(fā)明實(shí)施例提供的基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的方法，首先，獲取了 地震資料保幅處理后的偏移疊加地震數(shù)據(jù)，并獲取了地震工區(qū)內(nèi)鉆穿目的層的測(cè)井?dāng)?shù)據(jù)； 之后，根據(jù)所述地震數(shù)據(jù)和所述測(cè)井?dāng)?shù)據(jù)進(jìn)行井震標(biāo)定，并對(duì)所述地震數(shù)據(jù)進(jìn)行目的層構(gòu) 造解釋，確定待預(yù)測(cè)的目的層的頂?shù)走吔?之后，對(duì)所述地震數(shù)據(jù)進(jìn)行去噪處理，形成去噪 后地震數(shù)據(jù)體;確定去噪后地震數(shù)據(jù)體的相干屬性;確定去噪后地震數(shù)據(jù)體的螞蟻?zhàn)粉檶?性;將所述螞蟻?zhàn)粉檶傩缘闹涤蜃兓秶Ｕ秊樗鱿喔蓪傩缘闹涤蜃兓秶?；將值域?化范圍校正后的螞蟻?zhàn)粉檶傩耘c所述相干屬性進(jìn)行融合，形成裂縫預(yù)測(cè)敏感屬性體;根據(jù) 所述裂縫預(yù)測(cè)敏感屬性體在空間范圍內(nèi)的數(shù)值大小預(yù)測(cè)所述頂?shù)走吔鐑?nèi)的地下裂縫空間 展布形態(tài)?？梢?，本發(fā)明可以將地震相干屬性和地震螞蟻?zhàn)粉檶傩赃M(jìn)行融合，在保留裂縫空 間展布趨勢(shì)的同時(shí)提高了地下裂縫預(yù)測(cè)的精度，從而可以實(shí)現(xiàn)利用地震數(shù)據(jù)同時(shí)準(zhǔn)確預(yù)測(cè) 地下不同尺度裂縫的空間展布的形態(tài)的效果。
[0141] 對(duì)應(yīng)于上述圖1和圖2所示的方法實(shí)施例，如圖10所示，本發(fā)明實(shí)施例提供一種基 于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的裝置，包括：
[0142] 數(shù)據(jù)獲取單元31，用于獲取地震資料保幅處理后的偏移疊加地震數(shù)據(jù)，并獲取地 震工區(qū)內(nèi)鉆穿目的層的測(cè)井?dāng)?shù)據(jù)。
[0143] 頂?shù)走吔绱_定單元32,用于根據(jù)所述地震數(shù)據(jù)和所述測(cè)井?dāng)?shù)據(jù)進(jìn)行井震標(biāo)定，并 對(duì)所述地震數(shù)據(jù)進(jìn)行目的層構(gòu)造解釋，確定待預(yù)測(cè)的目的層的頂?shù)走吔纭?br>[0144] 去噪處理單元33,用于對(duì)所述地震數(shù)據(jù)進(jìn)行去噪處理，形成去噪后地震數(shù)據(jù)體。
[0145] 相干屬性確定單元34,用于確定去噪后地震數(shù)據(jù)體的相干屬性。
[0146] 螞蟻?zhàn)粉檶傩源_定單元35,用于確定去噪后地震數(shù)據(jù)體的螞蟻?zhàn)粉檶傩浴?br>[0147] 值域變化范圍校正單元36,用于將所述螞蟻?zhàn)粉檶傩缘闹涤蜃兓秶Ｕ秊樗?相干屬性的值域變化范圍。
[0148] 裂縫預(yù)測(cè)敏感屬性體生成單元37,用于將值域變化范圍校正后的螞蟻?zhàn)粉檶傩耘c 所述相干屬性進(jìn)行融合，形成裂縫預(yù)測(cè)敏感屬性體。
[0149] 地下裂縫空間展布形態(tài)預(yù)測(cè)單元38,用于根據(jù)所述裂縫預(yù)測(cè)敏感屬性體在空間范 圍內(nèi)的數(shù)值大小預(yù)測(cè)所述頂?shù)走吔鐑?nèi)的地下裂縫空間展布形態(tài)。
[0150] 此外，該去噪處理單元33,具體可以：
[0151] 根據(jù)中值濾波法對(duì)所述地震數(shù)據(jù)進(jìn)行去噪處理，形成去噪后地震數(shù)據(jù)體。
[0152] 所述去噪處理單元33,具體還可以：
[0153] 根據(jù)預(yù)先設(shè)置的時(shí)窗在所述地震數(shù)據(jù)中選取2y+l個(gè)樣本點(diǎn)。
[0154] 根據(jù)公式：=膨也川"&-糾，…，S,，….5，，，丨，先^ …,N;確定所述預(yù) 先設(shè)置的時(shí)窗中的2y+l個(gè)樣本點(diǎn)中第y+1個(gè)樣本點(diǎn)的值。
[0155] 將所述確定為地震數(shù)據(jù)中的樣本點(diǎn)Sx的去噪處理結(jié)果。
[0156] 此外，該相干屬性確定單元34,具體用于：
[0157] 獲取去噪后地震數(shù)據(jù)體在一預(yù)先設(shè)置的分析窗口內(nèi)的多道地震數(shù)據(jù)的多個(gè)采樣 點(diǎn)，形成地震數(shù)據(jù)矩陣D: dl2 …dln d，' d11 dry
[0158] D= ;l ;2 ? Xl 4:2 …
[0159] 其中，k為每道地震數(shù)據(jù)的采樣點(diǎn)個(gè)數(shù);n為所述預(yù)先設(shè)置的分析窗口內(nèi)的地震數(shù) 據(jù)的道數(shù)。
[0160] 確定地震數(shù)據(jù)矩陣D的協(xié)方差矩陣C: K
[0161 ] C = O1 〇 - ^ dk d[ k-\
[0162] 其中，<為地震數(shù)據(jù)矩陣D的第n行向量，(:// …,<4j，尺;dk為必 的轉(zhuǎn)置，dk=[dki，dk2,…，dkn]T。
[0163] 確定地震數(shù)據(jù)矩陣D對(duì)應(yīng)的相干屬性E:
[0165] 其中，Ak為協(xié)方差矩陣C的第k個(gè)特征值;Amax為協(xié)方差矩陣C的最大特征值。
[0166] 另外，所述螞蟻?zhàn)粉檶傩源_定單元35,具體可以：
[0167] 根據(jù)去噪后地震數(shù)據(jù)體確定多個(gè)地震數(shù)據(jù)節(jié)點(diǎn)。
/e媽蟻g允許走的節(jié)點(diǎn)
[0168] 根據(jù)公式： j運(yùn)螞蟻g允許走的節(jié)點(diǎn)
[0169] 確定螞蟻算法尋求路徑節(jié)點(diǎn)的轉(zhuǎn)移概率，并將所述轉(zhuǎn)移概率確定為所述螞蟻?zhàn)粉?屬性M。
[0170] 其中，i為螞蟻當(dāng)前所在地震數(shù)據(jù)節(jié)點(diǎn)位置；j為螞蟻行進(jìn)到下一地震數(shù)據(jù)節(jié)點(diǎn)位 置;T^U+n)為從t時(shí)刻開始經(jīng)過(guò)n次迀移后在地震數(shù)據(jù)節(jié)點(diǎn)i和地震數(shù)據(jù)節(jié)點(diǎn)j連線上殘留 的信息素;n^(t)為t時(shí)刻由地震數(shù)據(jù)節(jié)點(diǎn)i到地震數(shù)據(jù)節(jié)點(diǎn)j的期望程度;a、0分別表示信息 素和螞蟻行進(jìn)路線的相對(duì)重要程度。 ~卜十") = /)~(，) + (!-/))Ar"
[0171] 其中， a k ^為信息素的殘留程度；1-P為信息素的減少 at.. = 2^Av..
[ 'i=i y 程度；a 為一次循環(huán)中留在從地震數(shù)據(jù)節(jié)點(diǎn)i到地震數(shù)據(jù)節(jié)點(diǎn)j路徑上的信息素；Tlj(t) 為t時(shí)刻在地震數(shù)據(jù)節(jié)點(diǎn)i和地震數(shù)據(jù)節(jié)點(diǎn)j連線上殘留的信息素。
[0172]另外，所述值域變化范圍校正單元36,具體可以：
[0173] 獲取螞蟻?zhàn)粉檶傩訫的各采樣點(diǎn)Mi。
[0174] 根據(jù)公式
丨=1，2，…，7V;確定值域變化范圍校正 后的螞蟻?zhàn)粉檶傩?，的各采樣點(diǎn)，i;其中，Mmin為螞蟻?zhàn)粉檶傩訫的各采樣點(diǎn)的最小值;M max 為螞蟻?zhàn)粉檶傩訫的各采樣點(diǎn)的最大值;Emin為相干屬性E的最小值。
[0175] 另外，所述裂縫預(yù)測(cè)敏感屬性體生成單元37,具體可以：
[0176] 根據(jù)公式:Zii XEi，i = 1，2，…，N;確定裂縫預(yù)測(cè)敏感屬性體Zi;其中，Ei為相干 屬性E與各采樣點(diǎn)姐對(duì)應(yīng)位置的采樣點(diǎn)。
[0177] 值得說(shuō)明的是，本發(fā)明實(shí)施例提供的基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的裝置 的具體實(shí)現(xiàn)方式可以參見上述圖1和圖2對(duì)應(yīng)的方法實(shí)施例，此處不再贅述。
[0178] 本發(fā)明實(shí)施例提供的基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的裝置，首先，獲取了 地震資料保幅處理后的偏移疊加地震數(shù)據(jù)，并獲取了地震工區(qū)內(nèi)鉆穿目的層的測(cè)井?dāng)?shù)據(jù)； 之后，根據(jù)所述地震數(shù)據(jù)和所述測(cè)井?dāng)?shù)據(jù)進(jìn)行井震標(biāo)定，并對(duì)所述地震數(shù)據(jù)進(jìn)行目的層構(gòu) 造解釋，確定待預(yù)測(cè)的目的層的頂?shù)走吔?之后，對(duì)所述地震數(shù)據(jù)進(jìn)行去噪處理，形成去噪 后地震數(shù)據(jù)體;確定去噪后地震數(shù)據(jù)體的相干屬性;確定去噪后地震數(shù)據(jù)體的螞蟻?zhàn)粉檶?性;將所述螞蟻?zhàn)粉檶傩缘闹涤蜃兓秶Ｕ秊樗鱿喔蓪傩缘闹涤蜃兓秶?；將值域?化范圍校正后的螞蟻?zhàn)粉檶傩耘c所述相干屬性進(jìn)行融合，形成裂縫預(yù)測(cè)敏感屬性體;根據(jù) 所述裂縫預(yù)測(cè)敏感屬性體在空間范圍內(nèi)的數(shù)值大小預(yù)測(cè)所述頂?shù)走吔鐑?nèi)的地下裂縫空間 展布形態(tài)?？梢?，本發(fā)明可以將地震相干屬性和地震螞蟻?zhàn)粉檶傩赃M(jìn)行融合，在保留裂縫空 間展布趨勢(shì)的同時(shí)提高了地下裂縫預(yù)測(cè)的精度，從而可以實(shí)現(xiàn)利用地震數(shù)據(jù)同時(shí)準(zhǔn)確預(yù)測(cè) 地下不同尺度裂縫的空間展布的形態(tài)的效果。
[0179] 本領(lǐng)域內(nèi)的技術(shù)人員應(yīng)明白，本發(fā)明的實(shí)施例可提供為方法、系統(tǒng)、或計(jì)算機(jī)程序 產(chǎn)品。因此，本發(fā)明可采用完全硬件實(shí)施例、完全軟件實(shí)施例、或結(jié)合軟件和硬件方面的實(shí) 施例的形式。而且，本發(fā)明可采用在一個(gè)或多個(gè)其中包含有計(jì)算機(jī)可用程序代碼的計(jì)算機(jī) 可用存儲(chǔ)介質(zhì)（包括但不限于磁盤存儲(chǔ)器、CD-ROM、光學(xué)存儲(chǔ)器等）上實(shí)施的計(jì)算機(jī)程序產(chǎn) 品的形式。
[0180] 本發(fā)明是參照根據(jù)本發(fā)明實(shí)施例的方法、設(shè)備（系統(tǒng)）、和計(jì)算機(jī)程序產(chǎn)品的流程 圖和/或方框圖來(lái)描述的。應(yīng)理解可由計(jì)算機(jī)程序指令實(shí)現(xiàn)流程圖和/或方框圖中的每一流 程和/或方框、以及流程圖和/或方框圖中的流程和/或方框的結(jié)合。可提供這些計(jì)算機(jī)程序 指令到通用計(jì)算機(jī)、專用計(jì)算機(jī)、嵌入式處理機(jī)或其他可編程數(shù)據(jù)處理設(shè)備的處理器以產(chǎn) 生一個(gè)機(jī)器，使得通過(guò)計(jì)算機(jī)或其他可編程數(shù)據(jù)處理設(shè)備的處理器執(zhí)行的指令產(chǎn)生用于實(shí) 現(xiàn)在流程圖一個(gè)流程或多個(gè)流程和/或方框圖一個(gè)方框或多個(gè)方框中指定的功能的裝置。
[0181] 這些計(jì)算機(jī)程序指令也可存儲(chǔ)在能引導(dǎo)計(jì)算機(jī)或其他可編程數(shù)據(jù)處理設(shè)備以特 定方式工作的計(jì)算機(jī)可讀存儲(chǔ)器中，使得存儲(chǔ)在該計(jì)算機(jī)可讀存儲(chǔ)器中的指令產(chǎn)生包括指 令裝置的制造品，該指令裝置實(shí)現(xiàn)在流程圖一個(gè)流程或多個(gè)流程和/或方框圖一個(gè)方框或 多個(gè)方框中指定的功能。
[0182] 這些計(jì)算機(jī)程序指令也可裝載到計(jì)算機(jī)或其他可編程數(shù)據(jù)處理設(shè)備上，使得在計(jì) 算機(jī)或其他可編程設(shè)備上執(zhí)行一系列操作步驟以產(chǎn)生計(jì)算機(jī)實(shí)現(xiàn)的處理，從而在計(jì)算機(jī)或 其他可編程設(shè)備上執(zhí)行的指令提供用于實(shí)現(xiàn)在流程圖一個(gè)流程或多個(gè)流程和/或方框圖一 個(gè)方框或多個(gè)方框中指定的功能的步驟。
[0183] 本發(fā)明中應(yīng)用了具體實(shí)施例對(duì)本發(fā)明的原理及實(shí)施方式進(jìn)行了闡述，以上實(shí)施例 的說(shuō)明只是用于幫助理解本發(fā)明的方法及其核心思想；同時(shí)，對(duì)于本領(lǐng)域的一般技術(shù)人員， 依據(jù)本發(fā)明的思想，在【具體實(shí)施方式】及應(yīng)用范圍上均會(huì)有改變之處，綜上所述，本說(shuō)明書內(nèi) 容不應(yīng)理解為對(duì)本發(fā)明的限制。
【主權(quán)項(xiàng)】
1. 一種基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的方法，其特征在于，包括： 獲取地震資料保幅處理后的偏移疊加地震數(shù)據(jù)，并獲取地震工區(qū)內(nèi)鉆穿目的層的測(cè)井 數(shù)據(jù)； 根據(jù)所述地震數(shù)據(jù)和所述測(cè)井?dāng)?shù)據(jù)進(jìn)行井震標(biāo)定，并對(duì)所述地震數(shù)據(jù)進(jìn)行目的層構(gòu)造 解釋，確定待預(yù)測(cè)的目的層的頂?shù)走吔纾? 對(duì)所述地震數(shù)據(jù)進(jìn)行去噪處理，形成去噪后地震數(shù)據(jù)體； 確定去噪后地震數(shù)據(jù)體的相干屬性； 確定去噪后地震數(shù)據(jù)體的螞蟻?zhàn)粉檶傩裕? 將所述螞蟻?zhàn)粉檶傩缘闹涤蜃兓秶Ｕ秊樗鱿喔蓪傩缘闹涤蜃兓秶? 將值域變化范圍校正后的螞蟻?zhàn)粉檶傩耘c所述相干屬性進(jìn)行融合，形成裂縫預(yù)測(cè)敏感 屬性體； 根據(jù)所述裂縫預(yù)測(cè)敏感屬性體在空間范圍內(nèi)的數(shù)值大小預(yù)測(cè)所述頂?shù)走吔鐑?nèi)的地下 裂縫空間展布形態(tài)。2. 根據(jù)權(quán)利要求1所述的基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的方法，其特征在于，對(duì) 所述地震數(shù)據(jù)進(jìn)行去噪處理，形成去噪后地震數(shù)據(jù)體，包括： 根據(jù)中值濾波法對(duì)所述地震數(shù)據(jù)進(jìn)行去噪處理，形成去噪后地震數(shù)據(jù)體； 所述根據(jù)中值濾波法對(duì)所述地震數(shù)據(jù)進(jìn)行去噪處理，形成去噪后地震數(shù)據(jù)體，包括： 根據(jù)預(yù)先設(shè)置的時(shí)窗在所述地震數(shù)據(jù)中選取2y+l個(gè)樣本點(diǎn)； 根據(jù)公式：*5^=騰也"<兄，d 、，…,S、·.丨義，i，奸1，2,…， N;確定所述預(yù)先設(shè) 置的時(shí)窗中的2y+l個(gè)樣本點(diǎn)中第y+Ι個(gè)樣本點(diǎn)*的值； 將所述確定為地震數(shù)據(jù)中的樣本點(diǎn)Sx的去噪處理結(jié)果。3. 根據(jù)權(quán)利要求2所述的基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的方法，其特征在于，確 定去噪后地震數(shù)據(jù)體的相干屬性，包括： 獲取去噪后地震數(shù)據(jù)體在一預(yù)先設(shè)置的分析窗口內(nèi)的多道地震數(shù)據(jù)的多個(gè)采樣點(diǎn)，形 成地震數(shù)據(jù)矩陣D:其中，k為每道地震數(shù)據(jù)的采樣點(diǎn)個(gè)數(shù);η為所述預(yù)先設(shè)置的分析窗口內(nèi)的地震數(shù)據(jù)的 道數(shù)； 確定地震數(shù)據(jù)矩陣D的協(xié)方差矩陣C:其中，為地震數(shù)據(jù)矩陣D的第η行向量，<4? i，Κ Κ [ ;dk為4的轉(zhuǎn) 置，dk - [ dkl，dk2，· · ·，dkn ]; 確定地震數(shù)據(jù)矩陣D對(duì)應(yīng)的相干屬性E:其中，為協(xié)方差矩陣C的第k個(gè)特征值;Amax為協(xié)方差矩陣C的最大特征值。4. 根據(jù)權(quán)利要求3所述的基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的方法，其特征在于，確 定去噪后地震數(shù)據(jù)體的螞蟻?zhàn)粉檶傩?，包括?根據(jù)去噪后地震數(shù)據(jù)體確定多個(gè)地震數(shù)據(jù)節(jié)點(diǎn)；確定螞蟻算法尋求路徑節(jié)點(diǎn)的轉(zhuǎn)移概率，并將所述轉(zhuǎn)移概率確定為所述螞蟻?zhàn)粉檶傩?M； 其中，i為螞蟻當(dāng)前所在地震數(shù)據(jù)節(jié)點(diǎn)位置;j為螞蟻行進(jìn)到下一地震數(shù)據(jù)節(jié)點(diǎn)位置; (t+n)為從t時(shí)刻開始經(jīng)過(guò)η次迀移后在地震數(shù)據(jù)節(jié)點(diǎn)i和地震數(shù)據(jù)節(jié)點(diǎn)j連線上殘留的信息 素;n^(t)為t時(shí)刻由地震數(shù)據(jù)節(jié)點(diǎn)i到地震數(shù)據(jù)節(jié)點(diǎn)j的期望程度;α、β分別表示信息素和螞 蟻行進(jìn)路線的相對(duì)重要程度；，ρ為信息素的殘留程度；1 -ρ為信息素的減少程 度；A 為一次循環(huán)中留在從地震數(shù)據(jù)節(jié)點(diǎn)i到地震數(shù)據(jù)節(jié)點(diǎn)j路徑上的信息素；Tlj(t)為t 時(shí)刻在地震數(shù)據(jù)節(jié)點(diǎn)i和地震數(shù)據(jù)節(jié)點(diǎn)j連線上殘留的信息素。5. 根據(jù)權(quán)利要求4所述的基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的方法，其特征在于，將 所述螞蟻?zhàn)粉檶傩缘闹涤蜃兓秶Ｕ秊樗鱿喔蓪傩缘闹涤蜃兓秶?，包括?獲取螞蟻?zhàn)粉檶傩驭母鞑蓸狱c(diǎn)Mi;/=1，2,TV;確定值域變化范圍校正后的 螞蟻?zhàn)粉檶傩?，的各采樣點(diǎn)，i;其中，Mmin為螞蟻?zhàn)粉檶傩驭母鞑蓸狱c(diǎn)的最小值;M max為螞 蟻?zhàn)粉檶傩驭母鞑蓸狱c(diǎn)的最大值;Emin為相干屬性E的最小值。6. 根據(jù)權(quán)利要求5所述的基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的方法，其特征在于，將 值域變化范圍校正后的螞蟻?zhàn)粉檶傩耘c所述相干屬性進(jìn)行融合，形成裂縫預(yù)測(cè)敏感屬性 體，包括： 根據(jù)公式:Zii X Ei，i = 1，2，…，N;確定裂縫預(yù)測(cè)敏感屬性體Zi;其中，Ei為相干屬性 E與各采樣點(diǎn)Mi對(duì)應(yīng)位置的采樣點(diǎn)。7. -種基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的裝置，其特征在于，包括： 數(shù)據(jù)獲取單元，用于獲取地震資料保幅處理后的偏移疊加地震數(shù)據(jù)，并獲取地震工區(qū) 內(nèi)鉆穿目的層的測(cè)井?dāng)?shù)據(jù)； 頂?shù)走吔绱_定單元，用于根據(jù)所述地震數(shù)據(jù)和所述測(cè)井?dāng)?shù)據(jù)進(jìn)行井震標(biāo)定，并對(duì)所述 地震數(shù)據(jù)進(jìn)行目的層構(gòu)造解釋，確定待預(yù)測(cè)的目的層的頂?shù)走吔纾? 去噪處理單元，用于對(duì)所述地震數(shù)據(jù)進(jìn)行去噪處理，形成去噪后地震數(shù)據(jù)體； 相干屬性確定單元，用于確定去噪后地震數(shù)據(jù)體的相干屬性； 螞蟻?zhàn)粉檶傩源_定單元，用于確定去噪后地震數(shù)據(jù)體的螞蟻?zhàn)粉檶傩裕? 值域變化范圍校正單元，用于將所述螞蟻?zhàn)粉檶傩缘闹涤蜃兓秶Ｕ秊樗鱿喔蓪?性的值域變化范圍； 裂縫預(yù)測(cè)敏感屬性體生成單元，用于將值域變化范圍校正后的螞蟻?zhàn)粉檶傩耘c所述相 干屬性進(jìn)行融合，形成裂縫預(yù)測(cè)敏感屬性體； 地下裂縫空間展布形態(tài)預(yù)測(cè)單元，用于根據(jù)所述裂縫預(yù)測(cè)敏感屬性體在空間范圍內(nèi)的 數(shù)值大小預(yù)測(cè)所述頂?shù)走吔鐑?nèi)的地下裂縫空間展布形態(tài)。8. 根據(jù)權(quán)利要求7所述的基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的裝置，其特征在于，所 述去噪處理單元，具體用于： 根據(jù)中值濾波法對(duì)所述地震數(shù)據(jù)進(jìn)行去噪處理，形成去噪后地震數(shù)據(jù)體； 所述去噪處理單元，具體還用于： 根據(jù)預(yù)先設(shè)置的時(shí)窗在所述地震數(shù)據(jù)中選取2y+l個(gè)樣本點(diǎn)； 根據(jù)公式心仏川(Λ,丨，…，..丨，弋x=丨么_.·，Ν:確定所述預(yù)先設(shè) 置的時(shí)窗中的2y+l個(gè)樣本點(diǎn)中第y+Ι個(gè)樣本點(diǎn)的值； 將所述確定為地震數(shù)據(jù)中的樣本點(diǎn)Sx的去噪處理結(jié)果。9. 根據(jù)權(quán)利要求8所述的基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的裝置，其特征在于，所 述相干屬性確定單元，具體用于： 獲取去噪后地震數(shù)據(jù)體在一預(yù)先設(shè)置的分析窗口內(nèi)的多道地震數(shù)據(jù)的多個(gè)采樣點(diǎn)，形 成地震數(shù)據(jù)矩陣D:其中，k為每道地震數(shù)據(jù)的采樣點(diǎn)個(gè)數(shù);η為所述預(yù)先設(shè)置的分析窗口內(nèi)的地震數(shù)據(jù)的 道數(shù)； 確定地震數(shù)據(jù)矩陣D的協(xié)方差矩陣C:其中，4為地震數(shù)據(jù)矩陣D的第η行向量，ii/ lS/tS A-;dk為 <的轉(zhuǎn) 置，dk - [ dkl，dk2，· · ·，dkn ]; 確定地震數(shù)據(jù)矩陣D對(duì)應(yīng)的相干屬性E:其中，為協(xié)方差矩陣C的第k個(gè)特征值;Amax為協(xié)方差矩陣C的最大特征值。10. 根據(jù)權(quán)利要求9所述的基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的裝置，其特征在于， 所述螞蟻?zhàn)粉檶傩源_定單元，具體用于： 根據(jù)去噪后地震數(shù)據(jù)體確定多個(gè)地震數(shù)據(jù)節(jié)點(diǎn)；確定螞蟻算法尋求路徑節(jié)點(diǎn)的轉(zhuǎn)移概率，并將所述轉(zhuǎn)移概率確定為所述螞蟻?zhàn)粉檶傩?M； 其中，i為螞蟻當(dāng)前所在地震數(shù)據(jù)節(jié)點(diǎn)位置;j為螞蟻行進(jìn)到下一地震數(shù)據(jù)節(jié)點(diǎn)位置; (t+n)為從t時(shí)刻開始經(jīng)過(guò)η次迀移后在地震數(shù)據(jù)節(jié)點(diǎn)i和地震數(shù)據(jù)節(jié)點(diǎn)j連線上殘留的信息 素;n^(t)為t時(shí)刻由地震數(shù)據(jù)節(jié)點(diǎn)i到地震數(shù)據(jù)節(jié)點(diǎn)j的期望程度;α、β分別表示信息素和螞 蟻行進(jìn)路線的相對(duì)重要程度；^為信息素的殘留程度；1 -ρ為信息素的減少程 度；A 為一次循環(huán)中留在從地震數(shù)據(jù)節(jié)點(diǎn)i到地震數(shù)據(jù)節(jié)點(diǎn)j路徑上的信息素；Tlj(t)為t 時(shí)刻在地震數(shù)據(jù)節(jié)點(diǎn)i和地震數(shù)據(jù)節(jié)點(diǎn)j連線上殘留的信息素。11. 根據(jù)權(quán)利要求10所述的基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的裝置，其特征在于， 所述值域變化范圍校正單元，具體用于： 獲取螞蟻?zhàn)粉檶傩驭母鞑蓸狱c(diǎn)Mi;?.，:2, ·"，iV;確定值域變化范圍校正后的 螞蟻?zhàn)粉檶傩?，的各采樣點(diǎn)，i;其中，Mmin為螞蟻?zhàn)粉檶傩驭母鞑蓸狱c(diǎn)的最小值;M max為螞 蟻?zhàn)粉檶傩驭母鞑蓸狱c(diǎn)的最大值;Emin為相干屬性E的最小值。12. 根據(jù)權(quán)利要求11所述的基于地震數(shù)據(jù)預(yù)測(cè)地下裂縫空間展布的裝置，其特征在于， 所述裂縫預(yù)測(cè)敏感屬性體生成單元，具體用于： 根據(jù)公式:Zii X Ei，i = 1，2，…，N;確定裂縫預(yù)測(cè)敏感屬性體Zi;其中，Ei為相干屬性 E與各采樣點(diǎn)Mi對(duì)應(yīng)位置的采樣點(diǎn)。
【文檔編號(hào)】G01V1/50GK105929452SQ201610384344
【公開日】2016年9月7日
【申請(qǐng)日】2016年6月2日
【發(fā)明人】馬曉宇, 楊志祥, 李勁松, 于豪, 徐光成, 郭繼亮
【申請(qǐng)人】中國(guó)石油天然氣股份有限公司