度確定裂縫發(fā)育規(guī)律(S106)dS103和S104沒(méi)有次序上的關(guān) 系,它們是兩個(gè)分支。S103是利用正演的關(guān)系圖版得到,S104是來(lái)自實(shí)際地震數(shù)據(jù)計(jì)算得 至IJ。即圖1體現(xiàn)出本發(fā)明提供的海上窄方位角資料的方位各向異性屬性差裂縫預(yù)測(cè)方法主 要包括六個(gè)步驟:步驟一,利用測(cè)井?dāng)?shù)據(jù)設(shè)計(jì)基于巖石物理模型的方位各向異性正演模型, 獲得不同裂縫密度下正演的井旁道疊加道集。測(cè)井?dāng)?shù)據(jù)包括:縱波時(shí)差、橫波時(shí)差、密度測(cè) 井?dāng)?shù)據(jù)。
[0030] 步驟二,利用正演的井旁道疊加道集計(jì)算地震屬性,獲取屬性異常度與裂縫走向 的關(guān)系圖版。屬性異常度定義為:不同走向的裂縫屬性值相對(duì)于固定方位夾角的屬性值的 變化與平行于裂縫走向的方位屬性值的百分比絕對(duì)值。步驟Ξ,通過(guò)所述關(guān)系圖版獲取可 檢測(cè)的裂縫方位范圍;
[0031 ]步驟四,對(duì)疊前地震數(shù)據(jù)體進(jìn)行方位角劃分,獲取屬性差絕對(duì)值數(shù)據(jù)體。在本發(fā)明 中,根據(jù)實(shí)際地震數(shù)據(jù)的方位角范圍,把實(shí)際地震數(shù)據(jù)按照方位角度和覆蓋次數(shù)均衡的劃 分為兩個(gè)方位角道集數(shù)據(jù),然后對(duì)每一個(gè)位置點(diǎn)上,相應(yīng)的兩個(gè)分屬于不同方位角道集的 地震屬性值進(jìn)行相減后取絕對(duì)值,得到各個(gè)位置點(diǎn)上的屬性差絕對(duì)值,從而形成了一個(gè)包 含每個(gè)位置點(diǎn)上的屬性差絕對(duì)值的數(shù)據(jù)體。
[0032] 運(yùn)些屬性包括:均方根振幅、起始衰減頻率、能量比、85 %能量對(duì)應(yīng)的頻率、65 %能 量對(duì)應(yīng)的頻率、瞬時(shí)頻率、最大能量、低頻能量、總能量。
[0033] 疊前指的是地震數(shù)據(jù)進(jìn)行疊加W前。
[0034] 步驟五,通過(guò)所述屬性差絕對(duì)值數(shù)據(jù)體與井上儲(chǔ)層進(jìn)行對(duì)比,獲取裂縫密度;
[0035] 步驟六,通過(guò)已獲取的可檢測(cè)的所述裂縫方位范圍和裂縫密度分析裂縫發(fā)育特 征。
[0036] 其中,對(duì)應(yīng)圖1的步驟S101中,首先要利用聲波時(shí)差、密度測(cè)井?dāng)?shù)據(jù),進(jìn)行層位標(biāo) 定,在單井上確定裂縫儲(chǔ)層的位置,此處的"單井"指的是某一 口井。
[0037] 標(biāo)識(shí)出裂縫儲(chǔ)集層組的數(shù)量和深度間隔。本發(fā)明中使用填充流體類(lèi)型、裂縫密度、 巖石密度和縱橫波速度建立方位各向異性正演模型。通常,裂縫能夠引起嚴(yán)重的各向異性, 運(yùn)里的"各向異性"指的是:在地層中,巖石物理性質(zhì)(如:彈性模量)隨著方位的變化而變化 的性質(zhì)。
[0038]巖石中裂縫的密度、寬度和所含流體的成份都是影響地震各向異性的因素,因此 若要分析裂縫對(duì)各向異性的影響,疊前正演模擬是非常必要的。在疊前正演模擬中,需要建 立裂縫儲(chǔ)層的巖石物理模型。要建立巖石物理模型,需要知道巖石的縱橫波速度和密度,運(yùn) 些參數(shù)可W從井資料中獲取,用于計(jì)算巖石的彈性參數(shù)。運(yùn)里的"疊前"指的是疊加之前,即 地震數(shù)據(jù)還沒(méi)有疊加在一起。
[0039] 在疊前的正演模擬中,通過(guò)建立儲(chǔ)層的巖石物理模型,即:首先通過(guò)填充流體類(lèi) 型、裂縫密度、巖石密度和縱橫波速度參數(shù),創(chuàng)建巖石物理模型,然后通過(guò)對(duì)運(yùn)些參數(shù)的計(jì) 算,獲得巖石的各向異性彈性模量,最后通過(guò)射線(xiàn)追蹤,獲得地震射線(xiàn)在空間上的傳播位 置、走時(shí)及振幅,從而模擬出在運(yùn)個(gè)巖石物理模型中的地震道。
[0040]運(yùn)里的"縱橫波速度"指的是縱波速度和橫波速度。
[0041]運(yùn)里的"井資料"指的是聲波時(shí)差和密度測(cè)井?dāng)?shù)據(jù)。
[0042] 在建立巖石物理模型中,本發(fā)明中采用的是化dson楠圓縫等效介質(zhì)模型,運(yùn)里的 "Hudson楠圓縫等效介質(zhì)模型"指的是:為了從地震或測(cè)井?dāng)?shù)據(jù)中獲得裂縫的物理屬性,需 要用到某種理論把巖石裂縫的動(dòng)力學(xué)屬性與它的物理屬性聯(lián)系起來(lái),現(xiàn)今有Ξ種理論比較 流行,化dson楠圓縫等效介質(zhì)理論就是其中一種,化dson楠圓縫等效介質(zhì)模型就是在運(yùn)種 理論下的介質(zhì)模型。
[0043]假定整個(gè)介質(zhì)由未包含裂縫的彈性介質(zhì)和它內(nèi)部分布的狹長(zhǎng)楠圓形裂隙組成,使 用裂縫密度Θ和裂縫的扁率α來(lái)描述裂縫系統(tǒng),等效的剛度系數(shù)為:
[0044]
…
[004引在等式(1)中,楚I表示各向同性背景的剛度系數(shù),C|為一階校正量。由于多數(shù)裂 縫型儲(chǔ)層的裂縫發(fā)育角度通常較大且成組出現(xiàn),因此可簡(jiǎn)化為HTI(水平橫向各向同性)介 質(zhì)模型,在運(yùn)種情況下,化dson模型的各向同性背景的剛度系數(shù)Cj為:
[0051]在等式(2a)-(2e)中,λ,μ表示沒(méi)有裂隙巖石的拉梅系數(shù),可通過(guò)所述的巖石縱橫 波速度和密度求得,
[0054]在等式(3a)和(3b)中,p、Vp、Vs分別為所述從井資料中獲取的巖石密度和縱橫波速 度。
[00巧]此時(shí),化dson模型的剛度系數(shù)一階校正量為:
[0061] 在等式(4a)-(4e)中,化和化由裂隙來(lái)確定,當(dāng)包含流體時(shí),化和化為:
[0067] 在等式(6a)和(6b)中,霧,分別為包含流體的體積模量和剪切模量。
[0068] 通過(guò)確定填充流體類(lèi)型,計(jì)算出所包含流體的體積模量和剪切模量。進(jìn)一步把從 井資料中獲取的巖石密度和縱橫波速度及裂縫密度作為參數(shù),計(jì)算出井中的彈性參數(shù),建 立井中的彈性參數(shù)模型。使用建立的井中彈性參數(shù)模型,確定方位角、入射角,通過(guò)正演模 擬計(jì)算疊前地震反射在各個(gè)方位角的響應(yīng),獲得不同裂縫密度、不同方位角、不同入射角的 正演井旁道道集。在裂縫密度的選擇上,當(dāng)裂縫體密度小于1%時(shí),裂縫的地震各向異性強(qiáng) 度通常較小,當(dāng)裂縫體密度大于1%后,裂縫的地震各向異性強(qiáng)度迅速增大。因此,裂縫密度 增量不宜過(guò)大,參數(shù)可選擇為0 %、1 %、2%、3%。在方位角的選擇上,通常方位角在0° -18〇D 之間,因此方位角的增量可選擇為10°。在入射角的選擇上,可采用最小值0°、最大值30°、增 量為3°的辦法確定入射角。此處的方位角的增量是在步驟1中,建立正演模型時(shí)設(shè)置的不同 方位角,即:在使用填充流體類(lèi)型、裂縫密度、巖石密度和縱、橫波速度建立正演模型時(shí),設(shè) 置的方位角,其方位角分別為:0°、10°、20°-180°,增量為10°。通過(guò)正演模擬獲得不同裂縫 密度、不同方位角、不同入射角的正演井旁道道集,固定裂縫密度和方位角,把不同入射角 的正演井旁道道集進(jìn)行疊加,得到不同裂縫密度、不同方位角條件下的正演井旁道疊加道 集。在實(shí)施步驟S102之前,需要對(duì)疊前地震數(shù)據(jù)體進(jìn)行分析獲得實(shí)際的采集資料方位角。在 步驟S102中,利用所述正演的井旁道疊加道集計(jì)算各種地震屬性,獲取各種屬性異常度與 裂縫走向的關(guān)系圖版。利用所述不同裂縫密度、不同方位角條件下的正演井旁道疊加道集 計(jì)算各種地震屬性,地震屬性包括:均方根振幅、起始衰減頻率、能量比、85%能量對(duì)應(yīng)的頻 率、65%能量對(duì)應(yīng)的頻率、瞬時(shí)頻率、最大能量、低頻能量、總能量。分別計(jì)算所述各種地震 屬性的屬性異常度,屬性異常度定義為:不同走向的裂縫屬性值相對(duì)于固定方位夾角的屬 性值的變化與平行于裂縫走向的方位屬性值的百分比絕對(duì)值。
[0069]圖2為固定采集方位與變化裂縫走向模型。"固定"即"設(shè)定"的意思,設(shè)定實(shí)際采集 方位固定不變。在此模型中,實(shí)際采集的方位角在84°-134°之間,方位角間隔為50%裂縫的 起始走向與實(shí)際采集方位的84°相對(duì)應(yīng),此模型的屬性異常度計(jì)算公式為:
[0070]
在等式(7)中,Θ為裂縫走向,Α(θ)為裂 縫走向上的屬性值,F(xiàn)(0)為裂縫走向上的屬性異常度,Α(I50-ΘI)為固定方位夾角的屬性 值,Α(θ)為平行于裂縫走向的方位屬性值。在本發(fā)明中,Θ取值為0°、1〇°、20°、···180°。
[0071]利用各種屬性異常度,分別建立各種屬性異常度與裂縫走向的關(guān)系圖版,橫坐標(biāo) 為裂縫走向,縱坐標(biāo)為某