一種基于地震相位微變分析的致密薄儲層地震預(yù)測方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及石油地球物理勘探領(lǐng)域，特別是基于地震相位微變分析的致密薄儲層地震預(yù)測方法。
【背景技術(shù)】
[0002]一般把厚度小于1/4波長的儲層稱為薄儲層，目前薄儲層地震預(yù)測方法有兩類，一類是只用地震信息，如有分頻成像(譜分解)、地層切片等，另一類就是井震聯(lián)合反演，如井約束波阻抗反演、地質(zhì)統(tǒng)計(jì)學(xué)反演等。
[0003]分頻技術(shù)主要通過譜分解，將地震資料中的高頻成份分離出來，通過不同頻率的掃描，觀察振幅能量變化，再依據(jù)不同頻率對應(yīng)的不同調(diào)諧厚度來識別薄儲層(孫學(xué)棟，趙建儒，白軍等.譜分解技術(shù)在營爾凹陷長沙嶺地區(qū)薄儲層預(yù)測中的應(yīng)用.石油地球物理勘探，2011，46(增刊):72-75;徐勝峰，劉春園、季玉新等.分頻技術(shù)在塔河油田石炭系薄儲層預(yù)測中的應(yīng)用.石油天然氣學(xué)報，2011，33 (2):65-69;徐麗英，徐鳴杰，陳振巖.利用譜分解技術(shù)進(jìn)行薄儲層預(yù)測.石油地球物理勘探，2006，41 (3):299-303 ;鄧和平.薄層預(yù)測技術(shù)概述.內(nèi)蒙古石油化工，2009，17:82-83)。國內(nèi)各主要盆地目的層的地震波速度為3500-5000m/s之間，地震波的頻寬10_75Hz，因此該方法理論上能識別的最薄儲層為
11.7m，而在實(shí)際應(yīng)用中，由于高頻資料信噪比低、能量弱，不可能用到75Hz頻率數(shù)據(jù)預(yù)測儲層，因此，該方法適合于15m以上薄儲層預(yù)測；地層切片它以解釋的兩個等時界面為頂?shù)祝诘貙拥捻數(shù)捉缑骈g按照厚度等比例內(nèi)插出一系列的層面，沿這些內(nèi)插出的層面逐一生成切片，這些切片的屬性就代表了同一沉積時期地質(zhì)體的特征(王學(xué)習(xí)，畢建軍，王利功等.地層切片穿時現(xiàn)象對地震屬性的影響.物探與化探，2012，36 (I):94-98) ? 一般認(rèn)為地層切片技術(shù)可以突破地震1/4波長分辨率的極限，達(dá)到1/8波長(張軍華，周振曉，譚明友等.地震切片解釋中的幾個理論問題.石油地球物理勘探，2007，42 (3) 348-352)，即使這樣，該方法分辨儲層的能力也與分頻成像相近。
[0004]第二類薄儲層預(yù)測的方法就是井震聯(lián)合反演，該方法借助于井上豐富的低頻及高頻信息，提高了地震的分辨率，它分辨薄儲層的能力可以達(dá)到測井水平，但缺點(diǎn)也是顯而易見的，模型占的比例大了，受井的影響較大，基本等同于井間內(nèi)插，模型占的比例小了又達(dá)不到分辨薄層的目的(劉占族，張雷，霍麗娜等.地質(zhì)統(tǒng)計(jì)學(xué)反演在煤層氣薄儲層識別中的應(yīng)用.石油地球物理勘探，2012，47 (I):30-34;張明振.對測井約束地震波阻抗反演的理解與應(yīng)用.油氣地球物理，2006，4 (3):1-5.) ο
[0005]鄂爾多斯盆地三疊系延長組是一套以湖相沉積為主的地層，自下而上分為十個油層組，分別為長10 —長I油層組。巖性油藏主要分布在長4+5-長8油層組，為三角洲前緣沉積，沉積微相為水下分流河道，砂體薄、致密，為典型的致密儲層。鄂爾多斯盆地油區(qū)地震資料主頻30-40HZ，頻寬為10-70HZ，大套油層組在地震剖面上為一個地震反射同相軸，地震可識別，油層組中一般發(fā)育1-5套砂體，單砂體厚度5-20m，利用目前的薄儲層預(yù)測技術(shù)難以預(yù)測油層組中的單砂體厚度。

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

[0006]本發(fā)明的目的是提供一種基于地震相位微變分析的致密薄儲層地震預(yù)測方法，解決儲層厚度小于1/8波長致密薄儲層的預(yù)測難題，為致密薄層勘探、油田開發(fā)單砂體預(yù)測、水平井軌跡設(shè)計(jì)及油藏精細(xì)描述提供依據(jù)，提高致密薄儲層勘探開發(fā)效率。
[0007]本發(fā)明解決其技術(shù)問題采用的技術(shù)方案是:一種基于地震相位微變分析的致密薄儲層地震預(yù)測方法，其特征是:所述方法包括以下步驟:
[0008]步驟1:正演分析，確定同一反射軸內(nèi)砂體厚度變化引起地震反射軸相位變化幅度及同相軸在縱向上的時差變化大小，選擇目的層不同砂體層厚度的井進(jìn)行正演分析；
[0009]步驟2:基于巖性變化分析的層位解釋追蹤，對目的層位解釋追蹤；
[0010]利用地震相位變化預(yù)測薄儲層，在地震相位穩(wěn)定時，層位解釋仍追蹤波峰或波谷，當(dāng)相位發(fā)生變化時，首先判斷是否是由區(qū)域構(gòu)造引起，如果相鄰地震反射軸發(fā)生同樣的變化，則是由區(qū)域構(gòu)造引起，層位解釋仍然與之前一致，如果相鄰地震反射軸無變化，則是由砂體變化引起，層位解釋不能再追蹤波峰或是波谷，要按照原構(gòu)造趨勢追蹤同相軸的任意位置；
[0011]步驟3:時窗選擇，依據(jù)步驟2的層位解釋結(jié)果，根據(jù)正演分析得到了不同砂體厚度變化同相軸的時差確定時窗；
[0012]步驟4:屬性分析，在步驟3時窗內(nèi)，利用地震均方根振幅、分頻成像、地震相分析這些屬性進(jìn)行分析，并利用實(shí)際的井資料進(jìn)行對比標(biāo)定，來預(yù)測致密薄儲層的砂體厚度。
[0013]所述的屬性分析的時窗選擇依據(jù)為做步驟I正演分析時將目的層放到同一時間位置，通過正演，比較不同砂體厚度地震反射同相軸相位的變化，找出同相軸在縱向上時差大小。
[0014]所述的正演分析選擇目的層砂體較發(fā)育井，通過改變砂體的阻抗信息，使它變?yōu)槟康膶由绑w厚度不同多口井，進(jìn)行對比分析。
[0015]所述的時窗選擇中，如果砂體厚度的變化引起相位增加，反射同相軸下移，沿層的時窗不能包含向下移動的同相軸的最大振幅，反之，如果砂體厚度的變化引起相位減小，反射同相軸上移，那邊沿層的時窗不能包含向上移動的同相軸的最大振幅。
[0016]本發(fā)明的有益效果是:本方法基于巖性變化分析的地震層位解釋追蹤技術(shù)，以往對層位的解釋是在井的標(biāo)定下，多數(shù)是追蹤波峰或是波谷，特別是某一套地層對應(yīng)于一個反射同相軸的情況下，更是嚴(yán)格追蹤波峰或是波谷，本方法與普通層位解釋的區(qū)別在于利用地震相位變化預(yù)測薄儲層，判斷是否是由區(qū)域構(gòu)造引起，可以追蹤同相軸的任意位置，依此來突出地震相位的微小變化，進(jìn)而對致密薄儲層進(jìn)行預(yù)測。
【附圖說明】
[0017]下面結(jié)合實(shí)施例附圖對本發(fā)明進(jìn)一步說明。
[0018]圖1是本發(fā)明實(shí)施例的流程圖；
[0019]圖2是Yl井長81段測井解釋結(jié)果示意圖；
[0020]圖3是利用25Hz雷克子波制作的合成地震記錄；
[0021]圖4a是把長81段上部砂巖的阻抗變?yōu)榕c泥巖相當(dāng)后的合成地震記錄；
[0022]圖4b是把長81下段部分砂巖的阻抗變?yōu)榕c泥巖相當(dāng)后的合成地震記錄；
[0023]圖5a是Xl井長31段測井解釋成果圖；
[0024]圖5b是x2井長31段測井解釋成果圖；
[0025]圖6是過xl、x2井地震剖面圖。
【具體實(shí)施方式】
[0026]地震反射同相軸的相位變化多用來進(jìn)行斷層描述，當(dāng)?shù)貙又写嬖跀鄬訒r，地震反射同相軸的相位錯斷明顯，依此來進(jìn)行斷層識別，但利用地震相位變化進(jìn)行薄儲層預(yù)測還未見報道。
[0027]地震的一個反射同相軸是同相軸所對應(yīng)的一套地層內(nèi)多個地震反射的疊合響應(yīng)，當(dāng)?shù)貙觾?nèi)砂泥巖結(jié)構(gòu)穩(wěn)定時，地震反射同相軸的相位比較穩(wěn)定，當(dāng)砂泥結(jié)構(gòu)發(fā)生變化時，地震反射軸的相位就會發(fā)生微小變化，通過這一微小的變化，結(jié)合地質(zhì)特點(diǎn)，就可以對單一砂體進(jìn)行預(yù)測。
[0028]如圖1的本發(fā)明基于地震相位微變分析的致密薄儲層地震預(yù)測方法流程圖所示，本發(fā)明通過以下步驟實(shí)現(xiàn):
[0029]步驟1:正演分析，確定同一反射軸內(nèi)砂體厚度變化引起地震反射軸相位變化幅度及同相軸在縱向上的時差變化大小，選擇目的層不同砂體層厚度的井進(jìn)行正演分析。
[0030]做正演時將目的層放到同一時間位置，通過正演，比較不同砂體厚度地震反射同相軸相位的變化，找出同相軸在縱向上時差大小，為后續(xù)屬性分析的時窗選擇提供依據(jù)；
[0031]在做正演分析時，也可以選擇目的層砂體較發(fā)育井，通過改變砂體的阻抗信息，使它變?yōu)槟康膶由绑w厚度不同多口井，這樣更有利于對比分析。
[0032]通過地震正演并結(jié)合實(shí)際地震資料情況，建立不同厚度砂體變化與地震相位變化之間的關(guān)系。正演分析最方便快捷的方法是在strata中作合成記錄，在正演過程中選擇richer子波即可，子波的主頻要依據(jù)實(shí)際工區(qū)地震