[發(fā)明專利]一種屬性變化率逐次晉級的裂縫發(fā)育帶檢測方法

在審

• 申請?zhí)枺?/td>: 201410524363.5
• 申請日: 2014-10-08
• 公開（公告）號: CN105572732A
• 公開（公告）日: 2016-05-11
• 發(fā)明（設(shè)計）人: 王世星;梁志強(qiáng)
• 申請（專利權(quán)）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術(shù)研究院
• 主分類號: G01V1/30
• 分類號: G01V1/30;G06F19/00
• 代理公司: 北京思創(chuàng)畢升專利事務(wù)所 11218
• 代理人: 郭韞

技術(shù)領(lǐng)域

本發(fā)明屬于油田勘探開發(fā)領(lǐng)域，具體涉及一種屬性變化率逐次晉級的裂縫 發(fā)育帶檢測方法。

背景技術(shù)

三維地震振幅屬性的變化率作為一種非均勻儲層描述的手段，已成為碳酸 鹽巖地區(qū)縫洞型儲層密集發(fā)育帶及巖溶儲層發(fā)育有利目標(biāo)預(yù)測的重要技術(shù)手 段。但數(shù)學(xué)上連續(xù)的變化率計算由于被地質(zhì)觀測采樣中不得不采用的離散抽樣 的網(wǎng)格點(diǎn)所取代，不同的采樣數(shù)據(jù)與采用不同的步長，使得導(dǎo)數(shù)的計算與橫向 步長相關(guān)，以往采用的計算不能涵蓋不同尺度分量的全體，而使很多有效區(qū)帶 內(nèi)圈定的振幅變化率區(qū)不能成帶成片分布，不能準(zhǔn)確描述不同尺度的計算結(jié)果， 而使裂縫有效發(fā)育帶的圈定邊界模糊、有效性不明確。

基于屬性(或振幅)變化率的技術(shù)都是在地震疊后地震資料的基礎(chǔ)上展開 的，但是傳統(tǒng)的方法通常只是在相鄰的幾個方向、選用一套公式、在固定的網(wǎng) 格區(qū)間的地震道之間開展計算，這對于不同尺度、不同發(fā)育級別的變尺度裂縫 預(yù)測而言通常是不全面的、結(jié)果也是局限的，只代表大規(guī)模裂縫而不能獲取從 小到大不同尺度的裂縫的全體。現(xiàn)有技術(shù)存在的缺陷是：1)要么只能在相鄰的 一個方向的地震道之間開展計算，這對于具有一定方位的裂縫預(yù)測而言通常是 不準(zhǔn)確的；2)要么只能選擇一種高階計算式(如4階)，只能使用一種陣列(如 5x5點(diǎn))、一列濾波算子計算裂縫發(fā)育，不能預(yù)測不同尺度的裂縫，尤其是小尺 度的裂縫，還要需要預(yù)知裂縫的發(fā)育大致尺度。

綜上所述，屬性變化率的計算都是在疊后地震資料上展開的，原先連續(xù)的 計算被有限的觀測點(diǎn)限制而被離散化，傳統(tǒng)的方法通常只能靠人工選取與觀測 步長相關(guān)的一種固定的差分或高階差分格式、在固定點(diǎn)陣的相鄰地震道之間開 展計算，數(shù)學(xué)上這種處理技術(shù)對小尺度(高波數(shù))的裂縫有壓制作用，檢測不 了不同尺度尤其是小尺度裂縫的相對變化，因此檢測的成果體現(xiàn)不了裂縫的多 尺度性，在很多有效區(qū)帶內(nèi)只能斷斷續(xù)續(xù)檢測結(jié)果，而不能檢測成片成帶的不 同裂縫的整體發(fā)育區(qū)域，不能精確描述平面空間裂縫發(fā)育的整體分布區(qū)域。

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

本發(fā)明的目的在于解決上述現(xiàn)有技術(shù)中存在的難題，提供一種屬性變化率 逐次晉級的裂縫發(fā)育帶檢測方法，在繼承傳統(tǒng)的疊后地震資料的包括方位信息、 高階變化率等信息的基礎(chǔ)上，引入不同尺度、變高階各種級次計算，并使這種 信息疊置起來，挖掘出從大到小不同級次裂縫發(fā)育信息，提供空間有效范圍的 圈定，可以得到更為精確全面的裂縫發(fā)育有效帶的預(yù)測結(jié)果。

本發(fā)明是通過以下技術(shù)方案實(shí)現(xiàn)的：

一種屬性變化率逐次晉級的裂縫發(fā)育帶檢測方法，對疊后中心點(diǎn)地震道展 開變點(diǎn)陣的道集相關(guān)校正、對相鄰點(diǎn)的地震道開展變尺度的高階方位振幅變化 率計算、累計一定發(fā)育強(qiáng)度的不同尺度裂縫形成有效的裂縫發(fā)育區(qū)帶，并最終 推廣到整個工區(qū)，得到整個地震數(shù)據(jù)的最大方位高階屬性變化率屬性、多階逐 次晉級的不同尺度裂縫發(fā)育綜合預(yù)測結(jié)果。

所述方法包括以下步驟：

(1)獲取疊后三維地震數(shù)據(jù)體；

(2)第k階計算所需道集方陣選擇與坐標(biāo)重置：對步驟(1)得到的疊后 三維地震數(shù)據(jù)體中的每一個地震道的每一個樣點(diǎn)，取其為中心點(diǎn)，選擇相鄰縱 橫向(k+1)*(k+1)相鄰道，構(gòu)建第k階道集方陣，并將該中心點(diǎn)置為原點(diǎn)(0，0)， 然后對其相鄰各點(diǎn)的坐標(biāo)進(jìn)行重置；

(3)計算中心道與相鄰道相關(guān)時差：計算中心道與相鄰每一道的互相關(guān)， 取其相關(guān)最大值時的時移量進(jìn)行道間時差校正，消除地層起伏的影響；

(4)利用第k階濾波算子計算不同方位屬性變化率；

(5)提取第k階不同方位屬性變化率最大值；

(6)提取第k階不同方位屬性變化率平均值；

(7)提取各階不同方位屬性最大值與平均值：在每個時間樣點(diǎn)重復(fù)步驟(2) -(6)，即可獲得k＝2，4，6，…N各階不同方位屬性最大值與平均值；其中， N為用戶選定的任意最大整數(shù)，可以選擇任意高的階數(shù)，但到達(dá)一定階后，基本 就不變了，具體實(shí)施時，用N＝30效果較好；