度的裂縫進(jìn)行整體集合，圈定儲層最有效發(fā)育區(qū)，評價(jià)最有利油氣鉆 探的"甜點(diǎn)"。在傳統(tǒng)的疊后地震資料的基礎(chǔ)上，盡可能挖掘更多的方位、尺度信息，自動完 成尺度選擇與方位綜合，使其在多階屬性變化率屬性中得以體現(xiàn)，相比傳統(tǒng)的疊后地震數(shù) 據(jù)得到的屬性變化率體，可以得到更為精確的裂縫預(yù)測結(jié)果。
【附圖說明】
[0036] 圖1，4階屬性變化率計(jì)算過程中，5x5點(diǎn)陣選擇和中心點(diǎn)、相鄰點(diǎn)地震道位置的重 置過程，并按縱橫、斜交4個(gè)方向取5點(diǎn)屬性計(jì)算變化率。為各k階屬性變化率的計(jì)算提供 類似指導(dǎo)；
[0037] 圖2,本發(fā)明方法的步驟框圖；
[0038] 圖3-1，某油田某區(qū)2階方位最大振幅變化率屬性圖；
[0039] 圖3-2,某油田某區(qū)4階方位最大振幅變化率屬性圖；
[0040] 圖3-3,某油田某區(qū)12階方位最大振幅變化率屬性圖；
[0041] 圖3-4,某油田某區(qū)120階方位最大振幅變化率屬性圖；
[0042] 圖4是本發(fā)明實(shí)施例中某油田某區(qū)2~20階多階方位屬性變化率逐次晉級屬性 圖。
【具體實(shí)施方式】
[0043] 下面結(jié)合附圖對本發(fā)明作進(jìn)一步詳細(xì)描述：
[0044] 本發(fā)明針對裂縫型油氣田勘探開發(fā)中有效裂縫儲層描述問題，設(shè)計(jì)了一種利用地 震屬性變化率的逐次晉級來檢測裂縫發(fā)育帶的方法，可應(yīng)用于油田勘探開發(fā)中裂縫型油氣 減有效延伸范圍的圈定。
[0045] 本發(fā)明通過采用高精度變尺度的多階表達(dá)式、多方位、逐次晉級的屬性變化率計(jì) 算方法，可以發(fā)掘并累計(jì)疊加疊后地震資料中與方位變化、尺度變化、累計(jì)變化相關(guān)的與裂 縫發(fā)育的相關(guān)信息，獲得與現(xiàn)有傳統(tǒng)的相干算法、方位振幅變化率等方法更高的有效性信 息，使得裂縫發(fā)育有效帶的延展范圍愈加明確，為油藏范圍的圈定與儲量計(jì)算提供更好的 效果。
[0046] (1)技術(shù)原理
[0047] ①屬性變化率計(jì)算方法
[0048] 通常屬性變化率被定義為：
[0050] 其中A(x，y，t)為三維地震體中一點(diǎn)處的某種屬性值，VAR(x，y，t)表示該點(diǎn)上沿 x，y方向上的變化量之矢量模，表示該屬性變化的強(qiáng)度。由于地震觀測在平面、縱向的離散 采樣性，在實(shí)際應(yīng)用中只能用不同的差分網(wǎng)格及正交組合進(jìn)行計(jì)算和平差。按數(shù)學(xué)的計(jì)算 原理，可選擇求導(dǎo)差分網(wǎng)格包括：前差商、后差商、中心前后差商、正交和交叉前后差商等多 種形式，在每個(gè)方向還可考慮兩種或多種差商的平均，提高有效精度；此外還需要對數(shù)據(jù)進(jìn) 行縱向的平差，提高對噪音抑制與平抑縱向采樣點(diǎn)的影響。通常采用計(jì)算中心點(diǎn)兩側(cè)上下 一定小時(shí)窗內(nèi)屬性變化率的平均，逐次滑動來反映該點(diǎn)的計(jì)算結(jié)果。
[0051] ②相關(guān)屬性變化率計(jì)算方法
[0052] 當(dāng)?shù)貙悠鸱兓瘯r(shí)，在高精度的計(jì)算中需要引入任一點(diǎn)的傾角信息，以使計(jì)算反 映的是儲層的變化而非地層結(jié)構(gòu)的變化。在此引入計(jì)算道間的相關(guān)時(shí)移從而將完全沿縱 橫向平面計(jì)算的屬性變化率，沿展到三維空間考慮屬性的空變矢量，其模即屬性的變化率 為：
[0054] 其中Atx，Aty表示差商計(jì)算兩點(diǎn)間的波形道的相關(guān)時(shí)移量，取中心點(diǎn)道與被計(jì) 算另一道波形相關(guān)計(jì)算取得最大值的時(shí)移量。在計(jì)算過程中對相關(guān)性高的道才進(jìn)行部分時(shí) 移處理，表示基值反射確有起伏變化的影響；一般的相關(guān)系數(shù)則不予進(jìn)行時(shí)移上的處理，從 而考慮一些雜亂相下的部分屬性的變化。這樣的改進(jìn)使得屬性變化率的計(jì)算更加適合于地 層反射有一定起伏變化的裂縫區(qū)域。
[0055] ③4階屬性變化率計(jì)算方法
[0056] 數(shù)學(xué)上一個(gè)簡單的一次連續(xù)導(dǎo)數(shù)的計(jì)算，由于離散的關(guān)系，可采用不同網(wǎng)格的差 分格式進(jìn)行計(jì)算。為了提高計(jì)算精度并反映尺度的影響，作為一次導(dǎo)數(shù)的高階近似的一個(gè) 例證，以高次中4階屬性變化率計(jì)算方法來推導(dǎo)其計(jì)算公式。
[0057] 函數(shù)u(X)、u(x土h)和u(x±2h)(其中h為相鄰道采樣間隔，即⑶P間隔）為任一 個(gè)方向上以u(x)為中心點(diǎn)兩側(cè)的相鄰2個(gè)屬性，按Taylor展開可以得到一次導(dǎo)數(shù)的表達(dá) 式。按Taylor展開有：
[0062] 將上述四個(gè)展開式合并有：
[0065] 一次導(dǎo)數(shù)的4階表達(dá)式濾波算子為[1/6 -4/3 0 4/3 -1/6] = [0· 167-1. 33 0 1. 33 -0. 167]。
[0066] ④多階屬性變化率計(jì)算方法
[0067] 仿照上述（3)式的4項(xiàng)推導(dǎo)方式，分別按從2項(xiàng)、4項(xiàng)、6項(xiàng)……，逐項(xiàng)逐次遞推，就 可以獲得一次導(dǎo)數(shù)的k階表達(dá)式（k = 2,4,6,8,……）下不同階次的濾波算子（見表1)。 由于差分格式的奇對稱性，原點(diǎn)兩側(cè)各道系數(shù)值相同而符號相反，如4階[0.167 -1.33 0 1.33 -0.167]算子，可用[0 1.33 -0.167]標(biāo)識正分支系數(shù)，負(fù)分支系數(shù)取其相反符號即 可。這樣就獲得了多階高精度的屬性變化率的全體表達(dá)式，這為發(fā)展變尺度屬性變化率提 供了計(jì)算算子，表1為一次導(dǎo)數(shù)各階近視濾波算子表。
[0070] ⑵技術(shù)實(shí)現(xiàn)步驟
[0071] 本發(fā)明針對實(shí)際油田勘探開發(fā)中面臨的裂縫型儲層發(fā)育帶整體圈劃問題，采用的 方法是針對實(shí)際疊后地震資料開展具有方位的、多階屬性變化率計(jì)算、多階累計(jì)綜合描述 研究，對疊后中心點(diǎn)地震道展開變點(diǎn)陣的道集相關(guān)校正、對相鄰點(diǎn)的地震道開展變尺度的 高階方位振幅變化率計(jì)算、累計(jì)一定發(fā)育強(qiáng)度的不同尺度裂縫形成有效的裂縫發(fā)育區(qū)帶， 并最終將該步驟推廣到整個(gè)工區(qū)，得到整個(gè)地震數(shù)據(jù)的最大方位高階屬性變化率屬性、多 階逐次晉級的不同尺度裂縫發(fā)育綜合預(yù)測結(jié)果，進(jìn)而對該工區(qū)的有效裂縫發(fā)育區(qū)帶的圈定 提供手段。
[0072] 技術(shù)實(shí)現(xiàn)流程如圖2所示，包括第k階計(jì)算時(shí)選取道集進(jìn)行傾斜時(shí)差校正、計(jì)算不 同方位屬性變化率，取方位最大值與方位平均值，對k進(jìn)行循環(huán)，并取得多階累計(jì)逐次晉級 求和值，完成一個(gè)樣點(diǎn)的計(jì)算，在時(shí)間空間完成整個(gè)數(shù)據(jù)體后形成計(jì)算成果。具體步驟如 下：
[0073] ①地震三維疊后數(shù)據(jù)體
[0074] 按點(diǎn)距疊前采集的地震數(shù)據(jù)，對其進(jìn)行常規(guī)道集預(yù)處理、保幅處理與成像，可以獲 得反映地下地層反射特征的疊后成像地震數(shù)據(jù)體，形成縱向、橫向、深度三個(gè)方向的三維地 震數(shù)據(jù)體；通常該數(shù)據(jù)體覆蓋上百-上千平方公里的地面面積，由幾十萬-幾百萬個(gè)地震道 組成；
[0075] ②不同方位4階屬性變化率計(jì)算實(shí)現(xiàn)（4階例證）
[0076] 按k = 4階屬性變化率計(jì)算算子所涉及樣點(diǎn)，為一網(wǎng)格范圍為（5點(diǎn)X5點(diǎn)）的計(jì) 算方陣（如圖1所示），圖中被計(jì)算中心點(diǎn)屬性u坐標(biāo)被置為（0,0)，相鄰24點(diǎn)屬性坐標(biāo) 依次為 u (_2, _2)、u (_2, _1)、u (_2,0)、u (_2,1)、u (_2, 2) ;u (_1，_2)、u (_1，_1)、u (_1，0)、 u(-l，l)、u(-l，2) ;u(0, -2)、u(0, -l)、u(0，l)、u(0,2) ;u(l，-2)、u(l，-l)、u(l，0)、u(l， 1)、u (1，2) ;u (2, -2)、u (2, -1)、u (2,0)、u (2,1)、u (2, 2);-次導(dǎo)數(shù)的 4 階表達(dá)式濾波算子 為[01.33-0. 167]，也即[0. 167-1.3301.33-0. 167]，分別按 4 個(gè)方向：西-東向（w-e)、西 北-東南向（wn-es)、北-南向（n-s)、北東-南西向（ne-sw)計(jì)算不同方位的屬性變化率：
[0081] 形成不同方位4階屬性變化率計(jì)算結(jié)果。
[0082] ③不同方位多階屬性變化率計(jì)算實(shí)現(xiàn)
[0083] 4階計(jì)算用5x5點(diǎn)陣，6階計(jì)算使用7x7點(diǎn)陣，......，k階使用（k+1)x(k+Ι)點(diǎn) 陣，……；分別