本發(fā)明涉及石油勘探開發(fā)技術(shù)領(lǐng)域，特別涉及一種粗面巖分布的確定方法及裝置。

背景技術(shù)：

20世紀(jì)70年代中國首次在遼河盆地發(fā)現(xiàn)了大規(guī)模粗面巖，并在粗面巖中發(fā)現(xiàn)了大量的油氣資源，引起了油氣勘探領(lǐng)域的廣泛重視。粗面巖形成于火山活動(dòng)的溢流相和爆發(fā)相，屬于中性火山噴出巖。具體的，粗面巖是指以粗面巖為主要成分(含量)的巖石，成分相當(dāng)于正長巖，基質(zhì)為隱晶質(zhì)，主要礦物是粗面結(jié)構(gòu)的長石。其中，長石中含有大量新生鉀長石和鈉長石系的斜長石；具有斑狀結(jié)構(gòu)，塊狀構(gòu)造、角礫狀構(gòu)造、氣孔構(gòu)造；多呈灰色、灰綠色，它的顯著特點(diǎn)是斑晶發(fā)育。

粗面巖作為中性火山巖的一個(gè)重要種屬，不僅可以作為儲(chǔ)層賦存油氣，而且影響著砂巖儲(chǔ)層的分布。準(zhǔn)確刻畫粗面巖的分布范圍對(duì)于油氣勘探與開發(fā)而言具有重要的實(shí)際意義。然而由于粗面巖具有巖性特殊，地球物理響應(yīng)特征復(fù)雜，橫向變化快，分布規(guī)律性差等特點(diǎn)，粗面巖分布預(yù)測依舊是世界性難題。

從國內(nèi)外研究現(xiàn)狀來看，火山巖預(yù)測常采用重磁勘探、微地震技術(shù)和地震剖面的直接解釋等方法，而對(duì)于粗面巖來講，這些方法不具有針對(duì)性。此外，國內(nèi)有些學(xué)者嘗試以巖相分析結(jié)合地震反演技術(shù)進(jìn)行預(yù)測，但是由于單一地震技術(shù)的局限性，預(yù)測效果并不十分理想。

綜上所述，目前有必要提出一種粗面巖分布的確定方法，以準(zhǔn)確預(yù)測粗面巖的分布，從而指導(dǎo)粗面巖油氣藏的勘探與開發(fā)。

技術(shù)實(shí)現(xiàn)要素：

本發(fā)明的目的是提供一種粗面巖分布的確定方法及裝置，以準(zhǔn)確預(yù)測粗面巖的分布，從而指導(dǎo)粗面巖油氣藏的勘探與開發(fā)。

本發(fā)明的上述目的可采用下列技術(shù)方案來實(shí)現(xiàn)：

一種粗面巖分布的確定方法，包括：

通過對(duì)目標(biāo)油藏的巖心測井分析，選取指定電性參數(shù)；其中，所述指定電性參數(shù)包括：自然伽馬、密度、聲波時(shí)差；

將所述聲波時(shí)差與密度轉(zhuǎn)換為縱波阻抗，再將所述縱波阻抗與所述自然伽馬交會(huì)處理，獲得粗面巖地球物理響應(yīng)的閾值；

基于所述粗面巖地球物理響應(yīng)的閾值，根據(jù)所述自然伽馬響應(yīng)特征和縱波阻抗特征，基于純波保幅地震數(shù)據(jù)，聯(lián)合進(jìn)行自然伽馬反演和縱波阻抗反演，獲得粗面巖分布的第一預(yù)測結(jié)果。

在優(yōu)選的實(shí)施方式中，所述方法還包括：

進(jìn)行地震屬性分析處理，獲得粗面巖分布的第二預(yù)測結(jié)果；

進(jìn)行頻譜分解分析處理，獲得粗面巖分布的第三預(yù)測結(jié)果；

對(duì)所述粗面巖分布的第一預(yù)測結(jié)果、第二預(yù)測結(jié)果、第三預(yù)測結(jié)果求交集，確定粗面巖分布的最終預(yù)測結(jié)果。

在優(yōu)選的實(shí)施方式中，所述地震屬性分析處理，獲得粗面巖分布的第二預(yù)測結(jié)果的步驟包括：

從具有平均算法的振幅統(tǒng)計(jì)類屬性和復(fù)地震道統(tǒng)計(jì)類地震屬性中選取指定地震屬性；

根據(jù)所述指定地震屬性的變化特征，獲得所述粗面巖分布的第二預(yù)測結(jié)果。

在優(yōu)選的實(shí)施方式中，所述指定地震屬性包括：平均波峰振幅和平均反射強(qiáng)度。

在優(yōu)選的實(shí)施方式中，所述頻譜分解分析處理，獲得粗面巖分布的第三預(yù)測結(jié)果的步驟包括：

對(duì)地震數(shù)據(jù)進(jìn)行頻譜分解；

根據(jù)所述頻譜分解結(jié)果，獲取預(yù)定頻率的振幅調(diào)諧體；

根據(jù)所述振幅調(diào)諧體頻率切片特征，獲得所述粗面巖分布的第三預(yù)測結(jié)果。

在優(yōu)選的實(shí)施方式中，所述對(duì)地震數(shù)據(jù)進(jìn)行頻譜分解包括：

通過短時(shí)窗離散傅里葉變換將時(shí)間域的地震數(shù)據(jù)轉(zhuǎn)換到頻率域，獲得轉(zhuǎn)換后的數(shù)據(jù)體；其中，所述轉(zhuǎn)換后的數(shù)據(jù)體垂向上是連續(xù)變化的頻率，平面上是歸一化后的調(diào)諧振幅。

在優(yōu)選的實(shí)施方式中，所述預(yù)定頻率為5赫茲至125赫茲。

在優(yōu)選的實(shí)施方式中，所述預(yù)定頻率為21赫茲。

在優(yōu)選的實(shí)施方式中，閾值為自然伽馬大于95API，縱波阻抗大于10000克/立方厘米×米/秒。

一種粗面巖分布的確定裝置，包括：

電性參數(shù)選取模塊，用于通過對(duì)目標(biāo)油藏的巖心測井分析，選取指定電性參數(shù)；其中，所述指定電性參數(shù)包括：自然伽馬、密度、聲波時(shí)差；

閾值獲取模塊，用于將所述聲波時(shí)差與密度轉(zhuǎn)換為縱波阻抗，再將所述縱波阻抗與所述自然伽馬交會(huì)處理，獲得粗面巖地球物理響應(yīng)的閾值；

第一預(yù)測結(jié)果獲取模塊，用于基于所述粗面巖地球物理響應(yīng)的閾值，根據(jù)所述自然伽馬響應(yīng)特征和縱波阻抗特征，基于純波保幅地震數(shù)據(jù)，聯(lián)合進(jìn)行自然伽馬反演和縱波阻抗反演，獲得粗面巖分布的第一預(yù)測結(jié)果。

在優(yōu)選的實(shí)施方式中，所述裝置還包括：

第二預(yù)測結(jié)果獲取模塊，用于進(jìn)行地震屬性分析處理，獲得粗面巖分布的第二預(yù)測結(jié)果；

第三預(yù)測結(jié)果獲取模塊，用于進(jìn)行頻譜分解分析處理，獲得粗面巖分布的第三預(yù)測結(jié)果；

最終預(yù)測結(jié)果獲取模塊，用于對(duì)所述粗面巖分布的第一預(yù)測結(jié)果、第二預(yù)測結(jié)果、第三預(yù)測結(jié)果求交集，確定粗面巖分布的最終預(yù)測結(jié)果。

本發(fā)明的特點(diǎn)和優(yōu)點(diǎn)是：通過對(duì)目標(biāo)油藏的巖心測井分析，選取指定電性參數(shù)；然后根據(jù)指定電性參數(shù)獲得粗面巖地球物理響應(yīng)的閾值；進(jìn)一步基于所述粗面巖地球物理響應(yīng)的閾值，根據(jù)所述自然伽馬響應(yīng)特征和縱波阻抗特征，基于純波保幅地震數(shù)據(jù)，聯(lián)合進(jìn)行自然伽馬反演和縱波阻抗反演，獲得粗面巖分布的第一預(yù)測結(jié)果。由于所述第一預(yù)測結(jié)果為根據(jù)所述粗面巖地球物理響應(yīng)的閾值，并聯(lián)合進(jìn)行自然伽馬反演和縱波阻抗反演聯(lián)合獲得，能夠?qū)⒋置鎺r和圍巖中的干擾巖石種類進(jìn)行明顯區(qū)分，因此，能夠較為準(zhǔn)確地預(yù)測粗面巖的分布，從而指導(dǎo)粗面巖油氣藏的勘探與開發(fā)。

附圖說明

圖1是本申請實(shí)施方式中一種粗面巖分布的確定方法的流程圖；

圖2是本申請實(shí)施方式中一種自然伽馬與縱波阻抗(GR-IMP)交會(huì)圖；

圖3a是本申請實(shí)施方式中一種縱波阻抗反演平面圖，虛線范圍為預(yù)測粗面巖分布區(qū)；

圖3b是本申請實(shí)施方式中一種自然伽馬反演平面圖，虛線范圍為預(yù)測粗面巖分布區(qū)；

圖4是本申請實(shí)施方式中一種粗面巖分布的確定方法的流程圖；

圖5是本申請實(shí)施方式中一種粗面巖分布的確定方法的流程圖；

圖6a是本申請實(shí)施方式中一種平均反射強(qiáng)度平面圖，虛線范圍為預(yù)測粗面巖分布區(qū)；

圖6b是本申請實(shí)施方式中一種平均波峰振幅平面圖，虛線范圍為預(yù)測粗面巖分布區(qū)；

圖7是本申請實(shí)施方式中一種粗面巖分布的確定方法的流程圖；

圖8是本申請實(shí)施方式中一種振幅調(diào)諧體頻率切片，虛線范圍為預(yù)測粗面巖分布區(qū)；

圖9是本申請實(shí)施方式中一種最終預(yù)測粗面巖分布圖；

圖10是本申請實(shí)施方式中一種粗面巖分布的確定裝置的模塊示意圖；

圖11是本申請實(shí)施方式中一種粗面巖分布的確定裝置的模塊示意圖。

具體實(shí)施方式

下面將結(jié)合附圖和具體實(shí)施方式，對(duì)本發(fā)明的技術(shù)方案作詳細(xì)說明，應(yīng)理解這些實(shí)施方式僅用于說明本發(fā)明而不用于限制本發(fā)明的范圍，在閱讀了本發(fā)明之后，本領(lǐng)域技術(shù)人員對(duì)本發(fā)明的各種等價(jià)形式的修改均落入本申請所附權(quán)利要求所限定的范圍內(nèi)。

除非另有定義，本文所使用的所有的技術(shù)和科學(xué)術(shù)語與屬于本申請的技術(shù)領(lǐng)域的技術(shù)人員通常理解的含義相同。本文中在本申請的說明書中所使用的術(shù)語只是為了描述具體的實(shí)施方式的目的，不是旨在于限制本申請。

本發(fā)明的目的是提供一種粗面巖分布的確定方法及裝置，以準(zhǔn)確預(yù)測粗面巖的分布，從而指導(dǎo)粗面巖油氣藏的勘探與開發(fā)。

請參閱圖1，本申請實(shí)施方式中提供的粗面巖分布的確定方法可以包括如下步驟。

步驟S10：通過對(duì)目標(biāo)油藏的巖心測井分析，選取指定電性參數(shù)；其中，所述指定電性參數(shù)包括：自然伽馬、密度、聲波時(shí)差。

在本實(shí)施方式中，可以對(duì)目標(biāo)油藏進(jìn)行地球物理響應(yīng)特征分析。具體的，可以根據(jù)目標(biāo)油藏井內(nèi)獲取的巖心和測井資料，通過巖心刻度測井分析巖心的電性響應(yīng)特征，從而選取對(duì)粗面巖比較敏感的電性參數(shù)。例如，可以選取對(duì)粗面巖比較敏感的自然伽馬、密度和聲波時(shí)差三個(gè)電性參數(shù)作為指定電性參數(shù)。

當(dāng)然，指定電性參數(shù)均不限于上述描述，所屬領(lǐng)域技術(shù)人員在本申請的技術(shù)精髓啟示下，還可能做出其他的變更，但只要其實(shí)現(xiàn)的功能和效果與本申請相同或相似，均應(yīng)涵蓋于本申請保護(hù)范圍內(nèi)。

步驟S12：將所述聲波時(shí)差與密度轉(zhuǎn)換為縱波阻抗，再將所述縱波阻抗與所述自然伽馬交會(huì)處理，獲得粗面巖地球物理響應(yīng)的閾值。

在本實(shí)施方式中，可以將獲取的指定電性參數(shù)中的聲波時(shí)差和密度參數(shù)轉(zhuǎn)換為縱波阻抗，然后將所述縱波阻抗與所述自然伽馬交會(huì)處理，獲取粗面巖地球物理響應(yīng)的閾值。

其中，所述縱波阻抗是巖石中的縱波速度與巖石密度的乘積。具體的，可以將聲波時(shí)差取倒數(shù)，之后與密度相乘，再乘以10⁶從而轉(zhuǎn)換為縱波阻抗。

所述縱波阻抗可以表明應(yīng)力波在巖體中傳播時(shí)，運(yùn)動(dòng)著的巖石質(zhì)點(diǎn)產(chǎn)生單位速度所需要的擾動(dòng)力，其能反映巖石對(duì)能量傳遞的抵抗能力。

當(dāng)獲取所述縱波阻抗之后，可以將所述縱波阻抗與所述自然伽馬交會(huì)處理，獲得粗面巖地球物理響應(yīng)的閾值。所述閾值為所述粗面巖地球物理響應(yīng)的門檻值。其中，所述將所述縱波阻抗與所述自然伽馬交會(huì)處理具體可以為獲取所述縱波阻抗與所述自然伽馬交會(huì)圖。

具體的，請參閱圖2，本申請實(shí)施方式中提供了一種自然伽馬與縱波阻抗(GR-IMP)交會(huì)圖。其中，橫坐標(biāo)IMP(impedance)表示縱波阻抗，單位為克/立方厘米×米/秒(g/cm³×m/s)；縱坐標(biāo)GR(gamma-ray)表示自然伽馬，單位為(API)。從所述自然伽馬與縱波阻抗的交會(huì)圖中，可以獲取粗面巖地球物理響應(yīng)的閾值。其中，自然伽馬的閾值為95API，縱波阻抗的閾值為10000g/cm³×m/s。

步驟S14：基于所述粗面巖地球物理響應(yīng)的閾值，根據(jù)所述自然伽馬響應(yīng)特征和縱波阻抗特征，基于純波保幅地震數(shù)據(jù)，聯(lián)合進(jìn)行自然伽馬反演和縱波阻抗反演，獲得粗面巖分布的第一預(yù)測結(jié)果。

在本實(shí)施方式中，根據(jù)自然伽馬響應(yīng)特征及縱波阻抗特征，基于純波保幅地震數(shù)據(jù)聯(lián)合開展自然伽馬反演以及縱波阻抗反演，再結(jié)合粗面巖地球物理響應(yīng)特征的閾值可以確定粗面巖的分布，并生成初步預(yù)測結(jié)果。

其中，所述純波保幅地震數(shù)據(jù)是指未做振幅剪切、濾波以及動(dòng)平衡等疊后修飾性處理的地震數(shù)據(jù)。由于純波保幅地震數(shù)據(jù)未做上述處理，從而保留了更多關(guān)于巖性、孔隙度、流體成分等方面的信息，進(jìn)而有利于獲取準(zhǔn)確的預(yù)測結(jié)果。

其中，所述自然伽馬反演是指在分析儲(chǔ)層的巖性、電性特征的基礎(chǔ)上，選取多個(gè)對(duì)巖性區(qū)分敏感的儲(chǔ)層參數(shù)，與地震建立聯(lián)系，通過非線性函數(shù)映射和神經(jīng)網(wǎng)絡(luò)技術(shù)計(jì)算巖性參數(shù)數(shù)據(jù)體，以達(dá)到巖性識(shí)別和儲(chǔ)層預(yù)測的目的一種反演方法。具體的，所述自然伽馬反演主要包含合成記錄制作和層位標(biāo)定、主組分分析和模型估算等幾個(gè)步驟。

其中，所述縱波阻抗反演是指應(yīng)用地震數(shù)據(jù)基于褶積模型，根據(jù)反射系數(shù)推算地層波阻抗的地震反演方法。具體的，所述縱波阻抗反演過程可以包括：層位標(biāo)定，低頻模型建立，反演運(yùn)算等。

通過對(duì)目標(biāo)油藏進(jìn)行自然伽馬反演和縱波阻抗反演聯(lián)合分析，獲得的粗面巖分布的第一預(yù)測結(jié)果，具體的，請參閱圖3a至圖3b。其中，圖3a是本申請實(shí)施方式中一種縱波阻抗反演平面圖，虛線范圍為預(yù)測粗面巖分布區(qū)；圖3b是本申請實(shí)施方式中一種自然伽馬反演平面圖，虛線范圍為預(yù)測粗面巖分布區(qū)。

由于圍設(shè)在粗面巖周圍的圍巖包括縱波阻抗大于10000g/cm³×m/s的玄武巖。如若簡單地通過縱波阻抗對(duì)應(yīng)的閾值來判斷，可能將玄武巖包括在了預(yù)測結(jié)果中。在本實(shí)施方式中，通過所述縱波阻抗反演可以將縱波阻抗中閾值中大于10000g/cm³×m/s的玄武巖從粗面巖中排除。

由于圍設(shè)在粗面巖周圍的圍巖包括自然伽馬大于95API的泥巖。如若簡單地通過自然伽馬對(duì)應(yīng)的閾值來判斷，可能將泥巖包括在了預(yù)測結(jié)果中。在本實(shí)施方式中，通過所述自然伽馬反演可以將自然伽馬閾值中大于95API的泥巖從粗面巖中排除。

也就是說，通過聯(lián)合進(jìn)行自然伽馬反演和縱波阻抗反演，可以根據(jù)門檻值將反演結(jié)果中縱波阻抗大于10000g/cm³×m/s和自然伽馬大于95API的提取出來，較為準(zhǔn)確地獲得所述第一預(yù)測結(jié)果。

在本申請實(shí)施方式中，通過對(duì)目標(biāo)油藏的巖心測井分析，選取指定電性參數(shù)；然后根據(jù)指定電性參數(shù)獲得粗面巖地球物理響應(yīng)的閾值；進(jìn)一步基于所述粗面巖地球物理響應(yīng)的閾值，根據(jù)所述自然伽馬響應(yīng)特征和縱波阻抗特征，基于純波保幅地震數(shù)據(jù)，聯(lián)合進(jìn)行自然伽馬反演和縱波阻抗反演，獲得粗面巖分布的第一預(yù)測結(jié)果。由于所述第一預(yù)測結(jié)果為根據(jù)所述粗面巖地球物理響應(yīng)的閾值，并聯(lián)合進(jìn)行自然伽馬反演和縱波阻抗反演聯(lián)合獲得，能夠?qū)⒋置鎺r和圍巖中的干擾巖石種類進(jìn)行明顯區(qū)分，因此，能夠較為準(zhǔn)確地預(yù)測粗面巖的分布，從而指導(dǎo)粗面巖油氣藏的勘探與開發(fā)。

請參閱圖4，在一個(gè)實(shí)施方式中，所述粗面巖分布的確定方法還可以包括如下步驟。

步驟S16：進(jìn)行地震屬性分析處理，獲得粗面巖分布的第二預(yù)測結(jié)果；

步驟S18：進(jìn)行頻譜分解分析處理，獲得粗面巖分布的第三預(yù)測結(jié)果；

步驟S20：對(duì)所述粗面巖分布的第一預(yù)測結(jié)果、第二預(yù)測結(jié)果、第三預(yù)測結(jié)果求交集，確定粗面巖分布的最終預(yù)測結(jié)果。

在本實(shí)施方式中，可以基于獲取的地震數(shù)據(jù)，進(jìn)行地震屬性分析處理和頻譜分解分析處理，獲得粗面巖分布第二預(yù)測結(jié)果和第三預(yù)測結(jié)果，然后將所述第一預(yù)測結(jié)果、第二預(yù)測結(jié)果、第三預(yù)測結(jié)果求交集，從而確定粗面巖分布的最終預(yù)測結(jié)果。上述將所述第一預(yù)測結(jié)果、第二預(yù)測結(jié)果、第三預(yù)測結(jié)果求交集也可以理解為通過所述第二預(yù)測結(jié)果和第三預(yù)測結(jié)果對(duì)所述第一預(yù)測結(jié)果進(jìn)行修正，從而有利于獲得更為精確的預(yù)測結(jié)果。

當(dāng)然，對(duì)所述第一預(yù)測結(jié)果的修正方式均不限于上述描述，所屬領(lǐng)域技術(shù)人員在本申請的技術(shù)精髓啟示下，還可能做出其他的變更，但只要其實(shí)現(xiàn)的功能和效果與本申請相同或相似，均應(yīng)涵蓋于本申請保護(hù)范圍內(nèi)。

請參閱圖5，在一個(gè)實(shí)施方式中，所述地震屬性分析處理，獲得粗面巖分布的第二預(yù)測結(jié)果的步驟可以包括下述子步驟：

步驟S161：從具有平均算法的振幅統(tǒng)計(jì)類屬性和復(fù)地震道統(tǒng)計(jì)類地震屬性中選取指定地震屬性；

步驟S162：根據(jù)所述指定地震屬性的變化特征，獲得所述粗面巖分布的第二預(yù)測結(jié)果。

由于粗面巖的波阻抗遠(yuǎn)遠(yuǎn)大于圍設(shè)在其周圍的圍巖，在粗面巖與圍巖界面處反射系數(shù)會(huì)增加，進(jìn)而會(huì)導(dǎo)致產(chǎn)生較大的屬性異常。其中，反射系數(shù)是指反射波與入射波的振幅比。反射系數(shù)的作用可以用于描述地下反射界面的物理性質(zhì)，反射系數(shù)越大，說明相鄰地層的波阻抗差異越大，反射波的反射作用越強(qiáng)。具體的，所述反射系數(shù)的計(jì)算公式如下：

$Ri=\frac{{\mathrm{\ρ}}_{i+1}{v}_{i+1}-{\mathrm{\ρ}}_i{v}_i}{{\mathrm{\ρ}}_{i+1}{v}_{i+1}+{\mathrm{\ρ}}_i{v}_i}$

其中，Ri是反射系數(shù)，ρ_i是各層的密度，v_i是各層的速度。

利用上述原理，在進(jìn)行地震屬性分析處理時(shí)，可以首先提取具有平均算法的振幅統(tǒng)計(jì)類屬性和復(fù)地震道統(tǒng)計(jì)類地震屬性。然后，從上述屬性中優(yōu)選對(duì)粗面巖反映敏感和能模糊小異常(如砂泥巖變化引起)的指定地震屬性，如平均波峰振幅和平均反射強(qiáng)度等。再根據(jù)粗面巖與圍巖阻抗差異大，反射系數(shù)突然增大的特征帶判斷產(chǎn)生大的屬性異常區(qū)域。該異常區(qū)域可能為粗面巖與圍巖的交界面，據(jù)此預(yù)測粗面巖的分布，生成第二預(yù)測結(jié)果。

在一個(gè)實(shí)施方式中，對(duì)粗面巖反映敏感和能模糊小異常的指定地震屬性可以包括：平均波峰振幅和平均反射強(qiáng)度。

所述振幅是指振動(dòng)的物理量可能達(dá)到的最大值，通常以A表示，它能表示振動(dòng)的范圍和強(qiáng)度的物理量。

所述反射強(qiáng)度：其中h(t)是虛地震道，f(t)是實(shí)地震道。

其中，所述平均波峰振幅是指時(shí)窗內(nèi)波峰振幅的平均數(shù)。

其中，所述平均反射強(qiáng)度是指對(duì)反射強(qiáng)度進(jìn)行時(shí)窗內(nèi)平均。

上述平均波峰振幅和平均反射強(qiáng)度不僅對(duì)粗面巖反映敏感，而且能模糊小異常(如砂泥巖變化引起的異常波動(dòng))。

具體的，請參閱圖6a和圖6b。其中，圖6a是本申請實(shí)施方式中一種平均反射強(qiáng)度平面圖，虛線范圍為預(yù)測粗面巖分布區(qū)；圖6b是本申請實(shí)施方式中一種平均波峰振幅平面圖，虛線范圍為預(yù)測粗面巖分布區(qū)。當(dāng)圖6a和圖6b中預(yù)測的區(qū)域相疊加可以獲得所述粗面巖分布的第二預(yù)測結(jié)果。

請參閱圖7，在一個(gè)實(shí)施方式中，所述頻譜分解分析處理，獲得粗面巖分布的第三預(yù)測結(jié)果的步驟可以包括如下子步驟：

步驟S181：對(duì)地震數(shù)據(jù)進(jìn)行頻譜分解；

步驟S182：根據(jù)所述頻譜分解結(jié)果，獲取預(yù)定頻率的振幅調(diào)諧體；

步驟S183：根據(jù)所述振幅調(diào)諧體頻率切片特征，獲得所述粗面巖分布的第三預(yù)測結(jié)果。

對(duì)于粗面巖而言，其相對(duì)于圍巖，巖層厚度大，密度大，地震波在巖層中傳播速度大，反射系數(shù)也較大，對(duì)應(yīng)的地震波的頻率較低。其中，所述頻率是指單位時(shí)間內(nèi)完成周期性變化的次數(shù)，是描述周期運(yùn)動(dòng)頻繁程度的量，常用符號(hào)f或ν表示，單位為秒分之一，符號(hào)為s^-1。其中，所述頻譜是指一個(gè)復(fù)雜振幅信號(hào)，可以看成是由許多簡諧分量疊加而成，各簡諧分量各自的振幅、頻率和初相就是該復(fù)雜信號(hào)的頻譜。

在本實(shí)施方式中，對(duì)地震數(shù)據(jù)進(jìn)行頻譜分解，可以通過短時(shí)窗離散傅里葉變換將時(shí)間域的地震數(shù)據(jù)轉(zhuǎn)換到頻率域，獲得轉(zhuǎn)換后的數(shù)據(jù)體。其中，所述轉(zhuǎn)換后的數(shù)據(jù)體垂向上是連續(xù)變化的頻率，平面上是歸一化后的調(diào)諧振幅。

其中，所述短時(shí)窗離散傅里葉變換是指以時(shí)間為自變量的信號(hào)和以頻率為自變量的頻譜函數(shù)之間的一種變換關(guān)系，其公式為：

$X\left(k\right)=\underset{n=0}{\overset{N-1}{\Σ}}x\left(n\right)e^{-j\frac{2\mathrm{\π}}{N}kn}$

其中，0≤k≤N-1,X(k)為變換后的頻率域函數(shù)，x(n)原時(shí)間域函數(shù)；

其中，所述數(shù)據(jù)體具體為振幅調(diào)諧體，其是指通過對(duì)地震數(shù)據(jù)進(jìn)行離散傅里葉變換得到一種新的三維數(shù)據(jù)體，其垂向是連續(xù)變化的頻率，平面是歸一化之后的調(diào)諧振幅。所述振幅調(diào)諧體是頻譜分解技術(shù)的成果數(shù)據(jù)之一，能用于確定粗面巖頻率異常，預(yù)測粗面巖分布。其中，所述振幅調(diào)諧體頻率切片為沿著某一頻率提取的振幅調(diào)諧體的切片。

基于上述獲得的地震數(shù)據(jù)的頻譜分解結(jié)果，根據(jù)粗面巖低頻率域特征，可以制作預(yù)定頻率范圍的振幅調(diào)諧體進(jìn)行預(yù)測，生成第三預(yù)測結(jié)果。具體的，如圖8所示，為本申請實(shí)施方式中一種振幅調(diào)諧體頻率切片，虛線范圍為預(yù)測粗面巖分布區(qū)。

在一個(gè)實(shí)施方式中，所述預(yù)定頻率可以為5赫茲(Hz)至125赫茲。

根據(jù)地震資料采樣間隔△t，可以確定可恢復(fù)的最高頻率是1/2△t，稱為尼奎斯特頻率。本地震數(shù)據(jù)的采樣間隔是4ms，最高頻率就是125Hz。而模擬檢波器無法接受到9Hz以下的低頻成分，所以選擇5-125Hz的頻率范圍。由于粗面巖具有低頻率域特征，因此在制作振幅調(diào)諧體時(shí)，可以選擇5Hz至125Hz范圍制作，以獲得能敏感反映粗面巖的振幅調(diào)諧體頻率切片。

進(jìn)一步的，當(dāng)頻率為21Hz時(shí)，振幅調(diào)諧體頻率切片能最為敏感的反映粗面巖。

所述頻率21Hz為優(yōu)選的頻率，其可以根據(jù)振幅調(diào)諧體頻率切片上的異常分布確定。當(dāng)異常規(guī)律性較強(qiáng)，與測井資料符合較好，與反演結(jié)果相對(duì)接近的時(shí)候可以確定為優(yōu)選的頻率。當(dāng)然，所述優(yōu)選的頻率并不限定為21Hz，其也可以包括在21Hz左右波動(dòng)的數(shù)值，例如20Hz等，本申請?jiān)诖瞬⒉蛔骶唧w的限定。

在本實(shí)施方式中，通過對(duì)所述粗面巖分布的第一預(yù)測結(jié)果、第二預(yù)測結(jié)果、第三預(yù)測結(jié)果求交集，確定粗面巖分布的最終預(yù)測結(jié)果。

具體的，請參閱圖9，將所述第一預(yù)測結(jié)果、第二預(yù)測結(jié)果、第三預(yù)測的預(yù)測區(qū)域進(jìn)行疊加，獲得所述最終的粗面巖的分布，并生成最終預(yù)測結(jié)果。

在本申請實(shí)施方式中，通過對(duì)目標(biāo)油藏的巖心測井分析，選取指定電性參數(shù)；然后根據(jù)指定電性參數(shù)獲得粗面巖地球物理響應(yīng)的閾值；進(jìn)一步基于所述粗面巖地球物理響應(yīng)的閾值，根據(jù)所述自然伽馬響應(yīng)特征和縱波阻抗特征，基于純波保幅地震數(shù)據(jù)，聯(lián)合進(jìn)行自然伽馬反演和縱波阻抗反演，獲得粗面巖分布的第一預(yù)測結(jié)果。由于所述第一預(yù)測結(jié)果為根據(jù)所述粗面巖地球物理響應(yīng)的閾值，并聯(lián)合進(jìn)行自然伽馬反演和縱波阻抗反演聯(lián)合獲得，能夠?qū)⒋置鎺r和圍巖中的干擾巖石種類進(jìn)行明顯區(qū)分，因此，能夠較為準(zhǔn)確地預(yù)測粗面巖的分布，從而指導(dǎo)粗面巖油氣藏的勘探與開發(fā)。

進(jìn)一步的，通過分別進(jìn)行地震屬性分析處理和頻譜分解分析處理，獲得粗面巖分布的第二預(yù)測結(jié)果和第三預(yù)測結(jié)果；然后對(duì)所述粗面巖分布的第一預(yù)測結(jié)果、第二預(yù)測結(jié)果、第三預(yù)測結(jié)果求交集，可以獲得更為準(zhǔn)確的粗面巖分布預(yù)測結(jié)果，從而確定粗面巖的分布位置，以指導(dǎo)粗面巖油氣藏的勘探與開發(fā)。

請參閱圖10，本申請實(shí)施方式中還提供了一種粗面巖分布的確定裝置，其可以包括：

電性參數(shù)選取模塊10，用于通過對(duì)目標(biāo)油藏的巖心測井分析，選取指定電性參數(shù)；其中，所述指定電性參數(shù)包括：自然伽馬、密度、聲波時(shí)差；

閾值獲取模塊12，用于將所述聲波時(shí)差與密度轉(zhuǎn)換為縱波阻抗，再將所述縱波阻抗與所述自然伽馬交會(huì)處理，獲得粗面巖地球物理響應(yīng)的閾值；

第一預(yù)測結(jié)果獲取模塊14，用于基于所述粗面巖地球物理響應(yīng)的閾值，根據(jù)所述自然伽馬響應(yīng)特征和縱波阻抗特征，基于純波保幅地震數(shù)據(jù)，聯(lián)合進(jìn)行自然伽馬反演和縱波阻抗反演，獲得粗面巖分布的第一預(yù)測結(jié)果。

請參閱圖11，所述粗面巖分布的確定裝置的另一種實(shí)施方式中，所述裝置還可以包括：

第二預(yù)測結(jié)果獲取模塊16，用于進(jìn)行地震屬性分析處理，獲得粗面巖分布的第二預(yù)測結(jié)果；

第三預(yù)測結(jié)果獲取模塊18，用于進(jìn)行頻譜分解分析處理，獲得粗面巖分布的第三預(yù)測結(jié)果；

最終預(yù)測結(jié)果獲取模塊20，用于對(duì)所述粗面巖分布的第一預(yù)測結(jié)果、第二預(yù)測結(jié)果、第三預(yù)測結(jié)果求交集，確定粗面巖分布的最終預(yù)測結(jié)果。

上述實(shí)施方式公開的粗面巖分布的確定裝置與本申請粗面巖分布的確定方法實(shí)施方式相對(duì)應(yīng)，可以實(shí)現(xiàn)本申請的粗面巖分布的確定方法實(shí)施方式并達(dá)到方法實(shí)施方式的技術(shù)效果。

本說明書中的上述各個(gè)實(shí)施方式均采用遞進(jìn)的方式描述，各個(gè)實(shí)施方式之間相同相似部分相互參照即可，每個(gè)實(shí)施方式重點(diǎn)說明的都是與其他實(shí)施方式不同之處。

以上所述僅為本發(fā)明的幾個(gè)實(shí)施方式，雖然本發(fā)明所揭露的實(shí)施方式如上，但所述內(nèi)容只是為了便于理解本發(fā)明而采用的實(shí)施方式，并非用于限定本發(fā)明。任何本發(fā)明所屬技術(shù)領(lǐng)域的技術(shù)人員，在不脫離本發(fā)明所揭露的精神和范圍的前提下，可以在實(shí)施方式的形式上及細(xì)節(jié)上作任何的修改與變化，但本發(fā)明的專利保護(hù)范圍，仍須以所附權(quán)利要求書所界定的范圍為準(zhǔn)。