Strata軟件操作手冊-8.10_第1頁
Strata軟件操作手冊-8.10_第2頁
Strata軟件操作手冊-8.10_第3頁
Strata軟件操作手冊-8.10_第4頁
Strata軟件操作手冊-8.10_第5頁
已閱讀5頁,還剩75頁未讀, 繼續(xù)免費(fèi)閱讀

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡介

1、中國石油STRATA地震反演軟件常規(guī)操作手冊編寫人:付建元川慶物探公司物探研究中心二0一0年七月目 錄一、地震反演原理及方法1(一)反演簡介1(二)STRATA軟件概述2(三)反演特點(diǎn)3(四)反演方法4二、地震反演的準(zhǔn)備工作9(一)項(xiàng)目工區(qū)建立9(二)數(shù)據(jù)準(zhǔn)備10三、疊后地震反演操作流程12(一)地震反演步驟12(二)數(shù)據(jù)輸入13(三)子波提取31(四)井震相關(guān)40(五)模型建立45(六)地震反演51(七)結(jié)果顯示輸出63四、疊后地震反演小結(jié)70(一)地震反演質(zhì)控70(二)反演方法推薦70應(yīng)用STRATA軟件進(jìn)行常規(guī)地震反演Hampson-Russell公司自1987年起從事地球物理勘探方面應(yīng)

2、用軟件的設(shè)計(jì)開發(fā)和咨詢服務(wù),為地球物理工程師提供了一系列卓越的技術(shù)平臺,涉及領(lǐng)域包括測井?dāng)?shù)據(jù)分析管理和合成記錄制作、地震反演、油藏屬性參數(shù)預(yù)測、AVO技術(shù)巖性與流體分析檢測、地質(zhì)統(tǒng)計(jì)制圖等。各軟件包功能獨(dú)特強(qiáng)大,可在PC和UNIX兩種平臺上快速、穩(wěn)定運(yùn)行,兩種平臺上的工區(qū)可相互轉(zhuǎn)換。Hampson-Russell卓越的技術(shù)平臺和強(qiáng)勁的咨詢支持為石油勘探工程師們提供了一系列優(yōu)秀的研究工具和手段。STRATA是一套交互的2D/3D模型和反演的軟件,它將疊后的地震道轉(zhuǎn)換成波阻抗地震道,然后可利用相互的關(guān)系計(jì)算速度道。除了地震反演外,STRATA軟件還包括其它一些功能如地震預(yù)處理,子波提取和井編輯等。

3、所有的任務(wù)可在多個連接的窗口中完成和顯示,用戶可以建立地質(zhì)模型,分析地震數(shù)據(jù)和交互地研究反演結(jié)果。在軟件中提供了多種反演方法如基于模型的反演、稀疏脈沖反演和人工智能神經(jīng)網(wǎng)絡(luò)等反演。一、地震反演原理及方法(一)反演簡介 地震反演是正演模型處理的反問題(圖1-1)。正演模型相對簡單,目前對正演技術(shù)是無可非議的。對某一給定的地質(zhì)模型,其地震響應(yīng)是相對穩(wěn)定的;而地震反演相對復(fù)雜,對某一給定的地震響應(yīng),可能對應(yīng)有多個地質(zhì)模型。圖1-1反演是從地震數(shù)據(jù)中提取它所包含的潛在地質(zhì)信息的過程。傳統(tǒng)上,反演是在疊后地震數(shù)據(jù)上進(jìn)行,目的是提取聲波速度和阻抗體。目前反演已經(jīng)被擴(kuò)展到疊前數(shù)據(jù)體,目的是既提取聲波阻抗又提

4、取橫波阻抗體,這樣就允許計(jì)算孔隙流體性質(zhì)。目前也可以利用反演結(jié)果直接預(yù)測巖性參數(shù)如孔隙度體。從而綜合疊前和疊后地震反演對儲層的巖性和流體分布進(jìn)行定量預(yù)測,最終實(shí)現(xiàn)精細(xì)預(yù)測有利儲層的目的。地震反演的常用軟件有HRS、EPS、Seislog和Jason等,本次將通過STRATA軟件介紹疊后地震反演的常規(guī)操作,熟悉最常見的地震反演的思路和流程。(二)STRATA軟件概述STRATA軟件如何為我們解決地震反演問題可容易地從圖1-2中看出,由疊加地震數(shù)據(jù)和速度/密度信息(以測井曲線或均方根速度形式)組成的兩個數(shù)據(jù)流導(dǎo)出一個基本的地震速度模型,然后用這個模型指導(dǎo)地震數(shù)據(jù)的整個反演。盡管這個目標(biāo)簡單,但在實(shí)

5、現(xiàn)時可能會遇見許多麻煩。該軟件主要包含以下工作。1)生成合成記錄2)交互壓縮/拉伸測井曲線3)提取子波4)疊后地震處理5)地震解釋6)通過插值建立模型7)用反演算法進(jìn)行地震反演它們之間的關(guān)系如圖1-3所示,具體在軟件使用過程中進(jìn)一步理解。圖1-2圖1-3(三)反演特點(diǎn)STRATA軟件疊后地震反演主要特點(diǎn)在于與反演相關(guān)的子波影響以及反演結(jié)果非唯一性的問題。疊后地震反演是我們分析地震道并試圖重構(gòu)地下的速度和阻抗結(jié)構(gòu)。反演基于的基礎(chǔ)模型是一維褶積模型: (1-1)其中: 以時間序列表示的零偏移距地下反射系數(shù) 地震子波,假設(shè)為不變的 隨機(jī)噪聲,與地震道不相關(guān)注意:1)在這個模型中,假設(shè)不存在多次波;2

6、)沒有考慮傳輸損失和幾何擴(kuò)散;3)沒有考慮頻率吸收;4)子波為常數(shù)。圖1-4反演可以看作是已知地震道,估算反射系數(shù)的過程。因?yàn)樵诘叵?,反射系?shù)是與地層的聲阻抗相關(guān)的: (1-2)其中:, 為密度,為縱波速度。反演可以認(rèn)為是相當(dāng)于估算下伏地層的聲阻抗。從圖1-4中(a)、(b)、(c)分別指示的子波、反射系數(shù)和地震道頻譜的變化可以看出,反射系數(shù)與子波褶積的結(jié)果是將地震頻譜中高頻和低頻消除了。由此理論上講,地震反演就是試圖將這些失去的頻率區(qū)域進(jìn)行恢復(fù)。所有地震反演算法也都存在非唯一性問題。可能存在多于一種地質(zhì)模型可以與地震數(shù)據(jù)相一致。要在這些可能的模型中決定一個模型,就需要去選擇地震數(shù)據(jù)以外的的一

7、些信息。通常用以下兩種方法來使用地震以外的信息:1)初始猜測模型的建立;2)最終結(jié)果與初始猜測模型背離的幅度約束。最后的反演結(jié)果既依賴于以上“其它信息”,也取決于地震數(shù)據(jù)(如圖1-2、1-3)。(四)反演方法STRATA軟件中疊后地震反演為我們提供了遞歸法(Bandlimited/Colored Inversion)、稀疏脈沖法(Sparse Spike)和基于模型法(Model Based)等反演方法。1、遞歸反演遞歸反演也稱為帶限反演,是最簡單和最早的一種疊后反演方法。在原理上,這種技術(shù)很明確。如果我們假設(shè)地震道描繪了地下反射系數(shù)的一個近似值,可以反演出反射系數(shù),最終得到聲阻抗。不過,因?yàn)?/p>

8、地震道是帶限的,反射頻譜的低頻(010Hz)和高頻(80250Hz)丟失,所以,這種處理不如想象的那樣容易。使它復(fù)雜化的可能還有地震道的相位誤差和噪聲污染。在STRATA程序中,帶限反演結(jié)果的產(chǎn)生涉及到三個步驟:1)用聲波曲線或均方根速度,或者二者組合導(dǎo)出一個低頻速度模型,約束高切頻率默認(rèn)為10Hz;2)用遞歸反演程序反演地震道,它是方程1-2的一個反問題(這幾乎等同于相位的-90 度轉(zhuǎn)換),其結(jié)果為聲阻抗或速度的中頻帶(1060Hz)數(shù)據(jù);3)將低頻和中頻信息結(jié)合,計(jì)算帶限反演的結(jié)果。在地震數(shù)據(jù)上拾取主要的反射層,并在測線上的一些有選擇的點(diǎn)提供詳細(xì)的速度信息,可以得到模型。對速度函數(shù)或控制點(diǎn)

9、進(jìn)行插值,得到要被反演的每個地震道的每個樣點(diǎn)速度值。拾取的同相軸作為插值的一個指導(dǎo)。假如沒有拾取同相軸,每個樣點(diǎn)均按線性進(jìn)行插值。因?yàn)樗俣饶P鸵呀?jīng)被定義,它描述了一個全帶寬的、但也許不均勻的速度模型。帶限反演的唯一參數(shù)是約束高切頻率(Constraint High-Cut Frequency),這個參數(shù)的出現(xiàn)是因?yàn)?,在遞歸反演中,我們希望把速度模型中的低頻成分加到地震推出的中頻模型中。缺省值為10Hz,是因?yàn)榈卣饠?shù)據(jù)中通常情況下不包括10Hz以下的有用能量。反演的第二部分涉及到反演地震道,生成中頻范圍數(shù)據(jù)。這一步按下述公式完成: (1-3)在實(shí)際應(yīng)用此公式之前,需將地震振幅歸一化到真反射系數(shù)-

10、0.1,0.1范圍內(nèi)。在帶限反演中,它按下述方法完成,計(jì)算地震道的均方根或最大樣點(diǎn)值,確定一個比例因子,用這個因子將數(shù)值歸一化為用戶指定的值。反演地震道遇到的另一個問題是,反演的值是聲阻抗而不是速度,因此,我們必須按照某種關(guān)系將聲阻抗轉(zhuǎn)換為速度。選擇的關(guān)系式為下述的Gardner公式: (1-4)其中,和是常數(shù)。假如沒有密度曲線,則采用=109和=4,與Gardner方程一致。假如有密度曲線,要做一個回歸分析,擬合出最佳的和值。將和代入方程1-4,即可將波阻抗轉(zhuǎn)換為速度。我們現(xiàn)在一般在PC機(jī)上做反演,處理時間可能是一個主要的考慮因素。因此,建議在一個窗口上完成反演。假如每次遞歸反演的起始點(diǎn)為不

11、同的采樣值,這樣做會使問題復(fù)雜化,加上道與道之間速度值在空間上的變化,會導(dǎo)致最終反演結(jié)果出現(xiàn)“條紋圖案”現(xiàn)象。STRATA程序中,反演前采用在反演道前加一個余弦斜坡的方法,使這種現(xiàn)象最小化。最后,將低頻模型值加到遞歸反演的中頻值上,就完成了帶限反演。顏色反演是對遞歸反演的修改,最早是由英國石油公司BP的Lancaster and Whitcombe在2000年SEG年會上提出來的。這種反演過程,尋找一個簡單操作因子O,對地震道進(jìn)行直接轉(zhuǎn)換為反演結(jié)果。通過比較地震數(shù)據(jù)和實(shí)際反演結(jié)果,他們得出結(jié)論,在頻率域確定的簡單操作因子相位是-90度。顏色反演結(jié)果非常類似于遞歸反演。其中一個區(qū)別是結(jié)果是相對聲

12、阻抗,既有正值也有負(fù)值(圖1-5)。2、稀疏脈沖反演稀疏脈沖反演是假定實(shí)際反射可以認(rèn)為是由一系列大脈沖里夾雜有小脈沖的背景,并假定只有大脈沖有意義,該方法通過檢查地震道來尋找大脈沖的位置。稀疏脈沖反演每次建立反射序列為一個脈沖,增加脈沖直到地震道被足夠準(zhǔn)確地進(jìn)行反演。圖1-5稀疏脈沖反演的要點(diǎn)是:1)只有當(dāng)?shù)卣饠?shù)據(jù)有脈沖才得到同相軸;2)試圖利用地震數(shù)據(jù)獲得符合地震數(shù)據(jù)的最簡單的可能模型;3)通常得到的結(jié)果比地質(zhì)本身實(shí)際產(chǎn)生的同相軸要少;4)它比基于模型反演比較少的依賴于初始猜測模型。3、基于模型反演STRATA中的基于模型反演又被稱為“塊狀”反演算法,是因?yàn)樗M梢粋€“塊狀”的偽速度曲線

13、。每個方塊的平均大小由用戶自定義,通常大于輸入數(shù)據(jù)的采樣率。反演中“塊狀”的偽速度曲線結(jié)果較從井信息計(jì)算得到的聲波或密度曲線的分辨率要低。圖1-3展示的就是適合塊狀反演的工作流程。對每一個要反演的道,STRATA在程序的建模部分給出一個最初始的猜測。根據(jù)一條測線上數(shù)個點(diǎn)的測井曲線或速度/時間對推出它們,并將其插值,最終使每個地震道都有一個對應(yīng)的阻抗道。在STRATA中有兩種使用約束的方法。一種是將其它信息做為一種“軟”約束,意思是,將最初猜測的聲阻抗作為一個單獨(dú)的信息部分,并將其與地震道以一定的加權(quán)的形式加在地震道上,這種方法稱為隨機(jī)法。第二種方法是,將其他信息做為一種“硬”約束,它設(shè)置一些最

14、終結(jié)果的絕對邊界,使得推出的最終模型與最初模型之間的誤差不能超出此邊界,這種方法叫做約束法。隨機(jī)反演認(rèn)為地震道和最初猜測的阻抗是數(shù)據(jù)的兩部分(可能會出現(xiàn)沖突),必須結(jié)合它們以提供最終的反演結(jié)果。約束反演把初始猜測約束作為反演的起始點(diǎn),設(shè)置任意參數(shù)與初始猜測偏差的絕對值或 “硬”邊界。遞歸反演的結(jié)果其頻帶與地震數(shù)據(jù)頻帶一樣?;谀P头囱莴@得了相對寬頻高頻的結(jié)果。問題是高頻成分可能來自于初始猜測模型而不是來自地震數(shù)據(jù)。這是遞歸反演和基于模型反演的比較(圖1-6)。通?;谀P头囱菘梢垣@得更細(xì)致信息,但這個結(jié)果實(shí)際上很類似。圖1-6稀疏脈沖反演生成了寬帶高頻反演結(jié)果,其結(jié)果與基于模型反演結(jié)果類似(圖

15、1-7),主要區(qū)別是缺少非常薄層的細(xì)節(jié)。圖1-7基于模型反演的要點(diǎn)是:1)地震有效分辨率得到提高; 2)子波估算的正確與否,直接影響反演的結(jié)果;3)反演結(jié)果可能在很大程度上依賴原始猜測初始模型,解決的方法是對初始模型進(jìn)行濾波;4)與其它反演方法一樣,存在多解性問題,在井少的工區(qū)通過增加地震的作用可以降低反演結(jié)果的多解性;5)在陸相盆地應(yīng)用該反演方法的效果相對更加實(shí)用。二、地震反演的準(zhǔn)備工作(一)項(xiàng)目工區(qū)建立STRATA可以在GEOVIEW主窗口下啟動,在WINDOWS(UNIX)系統(tǒng)的程序運(yùn)行處鍵入geoview或從開始/所有程序/HRS applications (Windows系統(tǒng))啟動G

16、EOVIEW 軟件,打開GEOVIEW的主窗口(圖2-1),出現(xiàn)STRATA的啟動按鈕。圖2-1圖2-2圖2-3點(diǎn)擊STRATA出現(xiàn)圖2-2的工區(qū)選擇模式窗口,可以啟動一個新的工區(qū)、打開已經(jīng)存在的工區(qū)以及打開先前正進(jìn)行的工區(qū)(均可能有多個)。這樣已經(jīng)運(yùn)行的地震項(xiàng)目工區(qū)均可直接啟動進(jìn)入。如圖2-3所示,新建的項(xiàng)目工區(qū)需選擇工區(qū)存放的路徑及工區(qū)名。在STRATA中完成一個項(xiàng)目的地震反演通常需要結(jié)合地震和井?dāng)?shù)據(jù),由此進(jìn)行地震反演需建立一個地震項(xiàng)目工區(qū)(*.prj),同時需建立一個井?dāng)?shù)據(jù)庫(*.wdb)。在圖2-1的GEOVIEW的主菜單里我們看到數(shù)據(jù)庫Database的選項(xiàng),點(diǎn)擊其中的New就會出現(xiàn)

17、圖2-4的建立一個新的井?dāng)?shù)據(jù)庫。其實(shí)在啟動GEOVIEW主程序時就會自動提示用戶打開已有的丼工區(qū)或新建一個井工區(qū)(圖2-5)。注意:1)地震反演項(xiàng)目工區(qū)或井工區(qū)的名稱以及存放其相關(guān)數(shù)據(jù)的路徑均不能含有中文,否則無法讀取相關(guān)信息;2)考慮到地震數(shù)據(jù)及反演過程產(chǎn)生的數(shù)據(jù)會越來越大,建議將工區(qū)及數(shù)據(jù)存放于剩余存儲空間較大的磁盤分區(qū)。圖2-4圖2-5(二)數(shù)據(jù)準(zhǔn)備如前所述,完成一個項(xiàng)目的地震反演通常需要結(jié)合地震數(shù)據(jù)和井?dāng)?shù)據(jù),那地震數(shù)據(jù)和井?dāng)?shù)據(jù)具體需要哪些呢?對于地震數(shù)據(jù),除了最基本的地震偏移數(shù)據(jù)以外,我們經(jīng)常已經(jīng)在Landmark等常規(guī)地震解釋軟件中完成的地震時間層位對比解釋成果,在該時間層位的基礎(chǔ)上

18、進(jìn)一步精細(xì)解釋,以滿足地震反演的層位約束的需要。對于井?dāng)?shù)據(jù),除了最基本的井位信息,我們還需要常用的測井?dāng)?shù)據(jù)(包括聲波時差、密度、自然伽瑪、孔隙度等),測井分層以及測井解釋成果。對于已有鉆井的油氣水顯示、測試、試采等資料也需收集。對于斜井,其井斜(井軌跡)數(shù)據(jù)也是必需的。同樣以上數(shù)據(jù)均需通過一定的編輯形成符合軟件錄入的分列式ASCII文本數(shù)據(jù)或常用的標(biāo)準(zhǔn)LAS格式文件。當(dāng)然,前人對該地區(qū)或鄰近區(qū)域所做的地質(zhì)研究、地震解釋成果圖件及相關(guān)成果報(bào)告等也進(jìn)行查詢和收集,通過查閱前人的工作,可以指導(dǎo)最新的地震反演工作。三、疊后地震反演操作流程(一)地震反演步驟在STRATA軟件概述部分就說明了疊后地震反演

19、通常需要進(jìn)行的操作步驟,如圖3-1所示,其主要的順序可以從STRATA軟件模塊界面左側(cè)的菜單看出,如圖3-2所示。圖3-1圖3-2第一項(xiàng)Data Manager模塊就是數(shù)據(jù)的管理器,包括地震、層位、井以及在反演過程中產(chǎn)生的子波、地震體等??梢詫σ陨蠑?shù)據(jù)進(jìn)行輸入、輸出、顯示、更改以及刪除等。對于地震數(shù)據(jù)和層位數(shù)據(jù)的輸入也可以通過Data Manager菜單下面的Seismic和Horizon模塊中進(jìn)行。井?dāng)?shù)據(jù)輸入在Database菜單里打開相應(yīng)井工區(qū)進(jìn)入Well Explorer里完成(圖2-1)。子波的設(shè)置和井震相關(guān)分析在Well和Wavelet模塊里完成;地震反演的相關(guān)操作在Model、An

20、alysis和Inversion模塊完成;反演結(jié)果的顯示及輸出主要在菜單File和View里進(jìn)行調(diào)整和操作。(二)數(shù)據(jù)輸入數(shù)據(jù)輸入主要包括與地震和測井相關(guān)的數(shù)據(jù)。在第二章數(shù)據(jù)準(zhǔn)備一節(jié)中已有較詳細(xì)的說明。這是進(jìn)行地震反演的基礎(chǔ)。1、加載地震數(shù)據(jù)進(jìn)入STRATA主窗口,單擊Data ManagerImport DataOpen SeismicFrom SEG-Y File(圖3-3)或單擊SeismicOpen Seismic from File(圖3-4),均可進(jìn)入圖3-3圖3-4地震數(shù)據(jù)的選擇界面(圖3-5),選擇SEG-Y文件所在的路徑,選中SEG-Y文件名,默認(rèn)文件后綴為sgy,若文件后綴

21、不是sgy的可改后綴或在Filter里選擇*.*即可(圖3-6)。圖3-5圖3-6圖3-7圖3-8圖3-9點(diǎn)擊圖3-5右下角的Next進(jìn)入定義數(shù)據(jù)的幾何形狀(觀測系統(tǒng))和文件格式的一系列界面。首先進(jìn)入圖3-7所示的選擇觀測系統(tǒng)類型,按實(shí)際工區(qū)類型選擇即可;點(diǎn)擊Next進(jìn)入圖3-8所示的選擇道頭中有無Inline&Xline號和X&Y坐標(biāo),對于三維工區(qū)一般Inline&Xline號和X&Y坐標(biāo)均有,對于二維工區(qū)一般無Inline&Xline號,僅有X&Y坐標(biāo);點(diǎn)擊Next進(jìn)入圖3-9的數(shù)據(jù)格式和測線道頭信息選擇的界面,一般數(shù)據(jù)格式多采用IBM浮點(diǎn),少數(shù)采用IEEE浮點(diǎn)格式。下面的Inline&X

22、line號和X&Y坐標(biāo)在道頭中的起始位置是準(zhǔn)確讀取地震工區(qū)位置的重要信息,點(diǎn)擊Header Dump就進(jìn)入了圖3-10所示的頭文件的顯示窗口,一般我們在Display里選中Trace Headers查看道頭信息。SEG-Y記錄一般是按照線號方向記錄的,線號在一定范圍內(nèi)不變,而CDP按步長1的規(guī)律遞增。線道號、坐標(biāo)位置確定后,在圖3-9中的Start分別填入Inline&Xline號和X&Y坐標(biāo)的起始位置,點(diǎn)擊OK或Next進(jìn)入地震數(shù)據(jù)的掃描提示窗口點(diǎn)擊Yes。圖3-10圖3-11掃描結(jié)果如圖3-11所示,通過線道號、道間距、坐標(biāo)、單位、方位以及工區(qū)視圖綜合判斷地震數(shù)據(jù)信息讀取的正確性。左下角有

23、一個Display in new seismic window,建議勾選,地震剖面將在一個新的窗口顯示(圖3-12),便于后期方便的調(diào)用該數(shù)據(jù),有利于下一步工作的開展。 圖3-122、加載地震解釋數(shù)據(jù)在圖3-12中左側(cè)的導(dǎo)航模塊可以看到Seismic下面就是Horizon模塊,點(diǎn)擊進(jìn)入圖3-13所示的界面,點(diǎn)擊Import HorizonsFrom File打開地震解釋層位選擇的界面,如圖3-13圖3-14所示??梢砸淮芜x擇一個層位文件加載,也可以一次選擇所有的層位文件一次性加入,還可以選擇一個包含多個層位的文件實(shí)現(xiàn)一次性加載所有的層位。點(diǎn)擊OK進(jìn)圖3-14圖3-15如圖3-15所示的層位文件

24、格式的指定窗口,若層位是從Landmark按默認(rèn)格式輸出可直接選擇文件類型Default Landmark,一般的文本文件默認(rèn)選擇Free Format,文件里是否有多個層位選擇No,對于三維工區(qū)一般按Inline and Xline Number導(dǎo)入,可以點(diǎn)擊Display the first selected file查看文件的具體內(nèi)容。點(diǎn)擊Nextt進(jìn)入層位名及各數(shù)據(jù)所在的列數(shù)指定界面,如圖3-16所示。此處也可以點(diǎn)擊Display selected file查看文件的具體內(nèi)容,便于準(zhǔn)確的填寫線道號及拾取層位的列數(shù)以及輸入層位數(shù)據(jù)跳過的行數(shù),單位和層位值的范圍。最下面的為層位名稱及顯示顏

25、色的指定。對于一個文件包含多個層位的情形,在圖3-15中的文件里有多個層位處選擇Yes,點(diǎn)擊下一步指定文件包含的層位個數(shù),通過點(diǎn)擊Display selected file查看文件內(nèi)容確定層位個數(shù),如圖3-17、18所示。圖3-16圖3-17圖3-18點(diǎn)擊下一步進(jìn)入層位名及各數(shù)據(jù)所在的列數(shù)指定界面,如圖3-19所示,與圖3-16類似,只是每個層位所對應(yīng)的列數(shù)不同而已,需分別指定。點(diǎn)擊OK就完成層位數(shù)據(jù)的加載。圖3-19層位加載完成在圖3-13中點(diǎn)擊地震顯示參數(shù)啟動按鈕,點(diǎn)擊Horizons,系統(tǒng)默認(rèn)選擇All Horizons in Project,如果用戶需自定義選擇User Defined

26、,在選擇需要顯示的層位,如圖3-20所示,這樣就可以在地震剖面上看到相應(yīng)的層位。圖3-20如圖3-21所示,在地震剖面上顯示相應(yīng)的層位名稱及相應(yīng)的顏色。對于層位加載是否正確,可以點(diǎn)擊Horizon菜單或在左側(cè)的Horizon導(dǎo)航模塊進(jìn)入均可選擇Display Horizons,出現(xiàn)層位選擇的窗口,如圖3-22所示,選擇導(dǎo)入的時間層位T3x5,按OK就立即顯示彈出該層位的平面展布圖。圖3-21圖3-22如圖3-23所示,該層位的平面展布過渡自然,說明加載的層位基本正確,層位拾取沒有明顯的問題。如圖3-24中黑圈所示,如果層位局部發(fā)生突變,可能拾取的層位有問題,鼠標(biāo)移至相應(yīng)位置即可顯示其線道號,可

27、以在其地震剖面上進(jìn)行查看手工修改。圖3-23圖3-24如圖3-25所示,拾取的層位局部存在問題,在圖3-21中的Horizon菜單中打開Pick Horizons,選擇已經(jīng)導(dǎo)入的層位后通過圖3-25中的方框所示的層位拾取控制方式進(jìn)行編輯,結(jié)合圖3-24的層位平面顯示,最終完成該層位的修正。圖3-253、加載井?dāng)?shù)據(jù)井?dāng)?shù)據(jù),在前一章已提到,包括井位(斜井需要井軌跡數(shù)據(jù))、測井分層、測井?dāng)?shù)據(jù)(包括聲波時差、密度、自然伽瑪、孔隙度等)信息,完成地震儲層反演還需要加載測井的解釋成果。如前所述,這些數(shù)據(jù)均可以在圖3-26所示的Well Explorer里完成。井的所有相關(guān)數(shù)據(jù)均在Import Data模塊

28、里進(jìn)行,加入的井?dāng)?shù)據(jù)可以在Curve View和Base Map里進(jìn)行圖形顯示查看。圖3-26如圖3-27所示,點(diǎn)擊Import Data打開其中的第一項(xiàng)即可加載地震反演需要的測井曲線、分層及井斜數(shù)據(jù)。圖3-27圖3-28如圖3-28所示,Import Logs界面中選擇要加入的測井?dāng)?shù)據(jù),下面有選擇測井文件格式的選項(xiàng),一般收集的資料多為非標(biāo)準(zhǔn)的LAS文本文件,此時選擇General ASCII格式;數(shù)據(jù)類型分為三類(測井、分層、井斜),按選擇的文件選上對應(yīng)的數(shù)據(jù)類型。針對一個地震工區(qū)的整理好的所有測井?dāng)?shù)據(jù),文件格式相同,對于井?dāng)?shù)不多且測井文件中測井曲線類型較少的情況可以一次性全選,完成同一類型

29、的井?dāng)?shù)據(jù)加載;當(dāng)測井文件中的測井曲線類型較多時建議單個文件加載,多個文件同時加載讀取文件響應(yīng)時間較長且容易導(dǎo)致加載錯誤。點(diǎn)擊Next進(jìn)入ASCII文件格式的指定界面,如圖3-29所示。Display Log File中可以選擇每一個文本文件,并通過View File查看具體內(nèi)容,如圖3-30所示,從而在下面的文件格式中進(jìn)行準(zhǔn)確的定義。一般測井?dāng)?shù)據(jù)深度位于第1列,單位為m(米)相對固定,文件頭所跳過的行數(shù)一般需要指定,忽略的空值多為-9999或-999.25,需查看文件明確指定。用于地震反演的常用測井曲線數(shù)據(jù)多為正值,我們建議可在值的范圍最小一項(xiàng)改為0,自動過濾小于0的樣點(diǎn)值,即可在加載時去掉文

30、件中所包含的野值。點(diǎn)擊Next進(jìn)入圖3-31井名的選擇或輸入,再點(diǎn)擊Next進(jìn)入圖3-32井坐標(biāo)的輸入,圖3-29圖3-30圖3-31圖3-32同時可以選擇井的類型(油井、氣井、干井等),點(diǎn)擊Next進(jìn)入圖3-33所示的測井曲線類型、名稱及單位的詳細(xì)指定,有部分不需要的測井曲線類型可以在最后一欄里不選擇輸入。點(diǎn)擊OK完成測井曲線數(shù)據(jù)的加載。圖3-33收集的測井?dāng)?shù)據(jù)會出現(xiàn)一口井有多井段(多個文件)組成,這樣可以分別加入每一段的測井?dāng)?shù)據(jù)。如圖3-33所示的井就由兩段組成,按照以上步驟或開始就可以同時選擇其兩個測井?dāng)?shù)據(jù)文件一起加載,如圖3-34所示,在左側(cè)的井?dāng)?shù)據(jù)列表中雙擊兩個聲波測井名,其對應(yīng)的曲

31、線就可以在Curve View里顯示。我們可以將加載的同一測井曲線的兩段數(shù)據(jù)合并在一起,便于后面反演工作的開展。圖3-34關(guān)于同一測井曲線的兩段數(shù)據(jù)的合并,需要在圖2-1所示的eLog模塊中完成。點(diǎn)擊進(jìn)入eLog窗口,如圖3-35所示。左側(cè)的模塊欄的下部有一個Math選項(xiàng),點(diǎn)擊進(jìn)入圖3-36所示的計(jì)算操作及描述界面,選擇Merge就可以將兩個或兩個以上的測井曲線合并到一起。如圖3-37,選擇好需合并的文件,在下面的Output Log Type中選擇測井曲線類型及單位,點(diǎn)擊Next系統(tǒng)自動指定輸出合并的測井曲線名,也可以自己指定。如圖3-38,根據(jù)測井實(shí)際情況,一般的在合并時的數(shù)據(jù)重疊區(qū)選擇深

32、度較淺的井段。點(diǎn)擊OK完成測井曲線合并,形成一條完整的測井曲線,如圖3-39所示。圖3-35圖3-36圖3-37圖3-38圖3-39對于井的分層數(shù)據(jù),如圖3-40,在井?dāng)?shù)據(jù)列表中點(diǎn)擊某一口井的Tops,打開對應(yīng)表格,可以直接從收集的分層資料復(fù)制粘貼到對應(yīng)位置即可,建議盡量使用規(guī)范的地質(zhì)層位名稱,其實(shí)在地震反演過程中常需要的測井解釋層段數(shù)據(jù),其信息可在分層數(shù)據(jù)中一并加入。有相應(yīng)的井分層文本文件,從圖3-28中,選擇數(shù)據(jù)類型Tops即可加載分層數(shù)據(jù),在圖3-41中指定分層數(shù)據(jù)各列所代表的類型,操作時可通過View File Contents查看具體內(nèi)容(圖3-42)。圖3-40圖3-41圖3-42

33、如果測井?dāng)?shù)據(jù)有使用Geoview的Well Explorer輸出的LAS格式的測井文件,在輸出時可以選擇測井Logs和分層Tops數(shù)據(jù)同時導(dǎo)出,可以輸出兩個單獨(dú)的文件,也可以將分層數(shù)據(jù)置入測井LAS文件的道頭Header中,如圖3.43所示,這樣輸出的LAS文件里包括了分層數(shù)據(jù),在輸入LAS文件時,測井的分層數(shù)據(jù)就一起加入到井?dāng)?shù)據(jù)庫里。圖3-43對于井斜數(shù)據(jù),圖3-43所示可與測井及分層一起輸出,在圖3-28中選擇Deviated Geometry,如圖3-44通過View File Contents查看井斜數(shù)據(jù),然后進(jìn)行輸入信息的詳細(xì)指定,一般來說,井斜數(shù)據(jù)有兩種類型:一種以X、Y的坐標(biāo)表示

34、井的軌跡,一種是以井斜角Dip和方位角Azimuth來表示。在圖3-45中根據(jù)3-44的文件內(nèi)容在Input deviation amplitude中選擇Dip and Azimuth in Degree,如果是坐標(biāo)表示的井軌跡這里就選擇X&Y in Meters,下面還有指定井斜角Dip和方位角Azimuth所在的列及數(shù)據(jù)范圍,深度一般為測量深度。點(diǎn)擊Next加入井名的指定就完成了該井井斜數(shù)據(jù)的加載。圖3-44圖3-45圖3-46如圖3-43所示,井?dāng)?shù)據(jù)加載完成可在Base Map中查看井的平面分布情況,通過Table View查看各井加載的相關(guān)數(shù)據(jù),通過Curve View可查看測井曲線、

35、分層及測井解釋成果,通過Summary可對測井?dāng)?shù)據(jù)進(jìn)行簡要的統(tǒng)計(jì)分析工作。(三)子波提取在井震相關(guān)前需要計(jì)算提取準(zhǔn)確的子波,子波的設(shè)置和提取在實(shí)際操作中是與井震相關(guān)交互完成的。如圖3-47所示,子波的設(shè)置和提取均在STRATA的Wavelet模塊中進(jìn)行。井震相關(guān)在eLog模塊中進(jìn)行,所以在eLog中同樣有Wavelet模塊,功能也一樣。1、子波提取參數(shù)用于計(jì)算合成記錄的子波,通過點(diǎn)擊WaveletDisplay Current Wavelet可以查看,如圖3-48所示,在軟件里默認(rèn)的子波是帶通子波(濾波頻率:5/10-50/60Hz),子波可用振幅譜和它的相位譜來定義,點(diǎn)擊Frequency標(biāo)

36、簽可在限定的頻率域顯示子波(圖3-49),點(diǎn)擊History即可看到該子波的詳細(xì)信息。實(shí)際工作中常用雷克子波來制作合成記錄,且使用不同的主頻,通過點(diǎn)擊WaveletCreate Ricker Wavelet進(jìn)入創(chuàng)建設(shè)置,建議雷克子波名稱里寫入雷克子波主頻值(圖3-50),便于后期的調(diào)用和查看。子波相位的極性在軟件里默認(rèn)協(xié)定是聲阻抗增加代表零相位地震數(shù)據(jù)的波峰,另一種協(xié)定是聲阻抗增加代表零相位地震數(shù)據(jù)的波谷,在WaveletSynthetic Polarity中定義,如圖3-51所示。實(shí)際中的子波極性通常使用的就是該軟件默認(rèn)的極性協(xié)定。圖3-47圖3-48圖3-49圖3-50圖3-51實(shí)際操作中

37、,除用ASCII文件讀取外部子波(通過圖3-47中的Data Manager導(dǎo)入Wavelet),子波提取有2種方法:只利用地震數(shù)據(jù)估算統(tǒng)計(jì)子波(Statistical);利用地震數(shù)據(jù)和測井曲線估算的全子波或常相位子波(Use Well),最終獲得滿足地震反演所需的最佳子波。圖3-52圖3-53圖3-54統(tǒng)計(jì)子波的提取,不使用井,只對地震數(shù)據(jù)進(jìn)行自相關(guān),假定相位是已知的。主要參數(shù):道范圍(默認(rèn)是整個數(shù)據(jù)體,通??稍O(shè)的大一些以增加統(tǒng)計(jì)運(yùn)算),如圖3-52所示。時間窗口(至少要為兩倍子波長度),如圖3-53所示,可直接給定起始時間,也可以通過目的層頂?shù)咨舷卵由靵碇付ㄌ崛〉臅r間窗口。子波長度(默認(rèn)值

38、為200ms,一般可選100ms),如圖3-54所示,斜坡長度一般用默認(rèn)值25ms。點(diǎn)擊OK開始運(yùn)算,完成即顯示所提取的統(tǒng)計(jì)子波(圖3-55),該子波為零相位是在提取子波時設(shè)置好的??梢栽O(shè)定該子波用于井震相關(guān)分析。決定是否利用井來確定子波的相位的關(guān)鍵問題是井曲線是否與地震數(shù)據(jù)匹配。井的相關(guān)壓縮/拉伸會引起相位的問題,利用測井曲線估算子波的振幅譜和相位譜,估算的結(jié)果與井和地震之間是否匹配緊密聯(lián)系,因此通常在利用井提取子波之前需要人工進(jìn)行井震相關(guān)使井曲線與地震數(shù)據(jù)匹配。圖3-55圖3-56如圖3-47,點(diǎn)擊WaveletExtract WaveletUse Well,進(jìn)入圖3-56所示的利用井提取

39、子波的設(shè)置窗口,可提取單井的子波,亦可提取多井的子波。在左下角的Set Active Logs For Each Well選項(xiàng)就是選擇每一口井的通過相關(guān)分析后的測井曲線。圖3-57所示,選擇地震體及離井的范圍,可選相鄰的道數(shù)或距離。圖3-58所示,選擇時間窗口。圖3-59所示,選擇子波長度和提取方式(全子波、常相位和Roy White算法)。圖3-57圖3-58圖3-59在井曲線與地震的相關(guān)比較好時,全子波提取的振幅譜和相位譜在時間域的算子可以將測井曲線反射系數(shù)成形于地震復(fù)合道上,如圖3-60。圖3-60圖3-59中的第二種子波提取方式是常相位,這個選項(xiàng)是利用地震道自相關(guān)求取子波的振幅譜,與地

40、震統(tǒng)計(jì)法完全一樣,利用井曲線估算一個近似的常相位值(圖3-61),這種方式比全子波法更穩(wěn)健,特別是當(dāng)井曲線與地震的相關(guān)比較差。由此該軟件中將常相位設(shè)為默認(rèn)的井子波提取方式。圖3-61使用Roy White先生的算法,提取過程與全子波提取法類似,另外可形成一個診斷圖允許用戶估計(jì)子波提取的誤差(圖3-62)。圖3-62如圖3-60、61、62所示,如果井匹配比較好,以上3種方式求取的子波很類似。子波提取的一般問題是:利用井提取子波,首先要對井作適宜的相關(guān);要合適地進(jìn)行相關(guān),需要有已知子波。實(shí)際子波提取的過程是:1)利用統(tǒng)計(jì)法求子波以確定初步的子波,假定子波的相位是已知的。2)通過拉伸/壓縮井曲線使

41、之與地震匹配。3)利用井曲線提取新的子波。4)通常要重復(fù)步驟2)和3),提取最佳子波。2、子波顯示編輯在圖3-47中的Display Wavelet Tool是子波顯示工具。在子波顯示工具里雙擊子波圖標(biāo)即可顯示子波。如圖3-63所示,其中編輯菜單中有一工具是Edit Phase.,可以顯示子波所有特征并允許進(jìn)行改變。其它編輯工具包括拷貝、刪除子波以及更改子波名稱。在File菜單里也提供了設(shè)置當(dāng)前子波和輸出、輸入子波等功能。圖3-633、子波分析診斷子波提取的診斷方法之一是井的相關(guān)系數(shù)值,如果不進(jìn)行拉伸/壓縮,通過time shift后達(dá)到最大的相關(guān)系數(shù),點(diǎn)擊ApplyShift,如圖3-64,

42、點(diǎn)擊Parameters打開互相關(guān)圖Cross Correlation顯示剩余相位誤差(圖3-65),當(dāng)旁瓣對稱時,顯示無剩余相位誤差。圖3-64圖3-65對于多口井時,可以用多井分析診斷,點(diǎn)擊WaveletMulti-well Analysis即可進(jìn)入選擇用于分析的井、子波的窗口,最后出現(xiàn)圖3-66的相關(guān)分析圖,圖中顯示所用的子波以及所有井的總的相關(guān)系數(shù),也可以查看每一口井的相關(guān)系數(shù)。點(diǎn)擊Error AnalysisCorrelation profile進(jìn)入圖3-67所示的相關(guān)系數(shù)繪圖,整體顯示每口井的相關(guān)系數(shù)圖,這可以用于標(biāo)記劣質(zhì)井,在以后的子波提取中可以排除掉,從而提高最終提取的綜合子波

43、的質(zhì)量。圖3-66圖3-67(四)井震相關(guān)井震相關(guān)就是應(yīng)用手動(交互圖形化)法將井曲線(深度域)與地震數(shù)據(jù)(時間域)進(jìn)行相關(guān),對深-時曲線關(guān)系進(jìn)行人工校正,從而優(yōu)化初始模型和地震數(shù)據(jù)間的相關(guān)性。若有井的Check Shot信息,一個Check Shot表格是由一組深度點(diǎn)計(jì)算的實(shí)際雙程旅行時組成的。井震相關(guān)應(yīng)在Check Shot校正之后進(jìn)行,理想情況下井相關(guān)對井的改變是小的變化。井震相關(guān)改變深-時曲線的方式與Check Shot校正一樣。井震相關(guān)是在合成記錄上選擇幾個同相軸使之對應(yīng)于實(shí)際地震道上。由于使用合成記錄,子波的選擇是至關(guān)重要的,由此井震相關(guān)和子波提取是交互進(jìn)行的。如圖3-68,在ST

44、RATA的左側(cè)導(dǎo)航模塊里選擇Well,點(diǎn)擊Well出現(xiàn)Edit/Correlate Well,即可選擇要進(jìn)行井震相關(guān)的井(圖3-69),選中其中的一口井后點(diǎn)擊Edit就進(jìn)入該井的相關(guān)分析界面eLog,如圖3-70所示。也可以直接點(diǎn)擊Geoview的主界面的eLog模塊進(jìn)入井震相關(guān)的分析界面。圖3-68圖3-69如圖3-70,在eLog的左側(cè)導(dǎo)航模塊里選擇Open Well可選擇工區(qū)內(nèi)其它井進(jìn)行相關(guān)分析。子波的設(shè)置提取在Wavelet模塊中進(jìn)行,具體操作參照上一節(jié)相關(guān)內(nèi)容。Seismic用于選擇井曲線相關(guān)的地震數(shù)據(jù)體。具體的相關(guān)操作在Correlate模塊中進(jìn)行。圖3-70圖3-71點(diǎn)擊Corr

45、elate,進(jìn)入圖3-71所示的提取地震復(fù)合道的菜單,一般使用默認(rèn)的Neighbours,即使用Inline和Xline旁邊(+/-1)道提取復(fù)合道。若為斜井,可選擇Distance,填寫目的層偏離井口的水平距離。點(diǎn)擊Ok進(jìn)入圖3-72所示的相關(guān)分析界面。左側(cè)縱向軸為時間,右側(cè)縱向軸為垂深。內(nèi)框一般從左到右分別顯示測井分層、聲波曲線、合成記錄道、地震復(fù)合道(井旁地震道)、過井地震剖面及解釋層位。圖3-72在圖3-72中所示的數(shù)據(jù)信息準(zhǔn)備齊全就可以開始相關(guān)分析工作,下面簡述相關(guān)分析的基本過程。下面的Snap to:選項(xiàng)默認(rèn)為Snap Peak/Trough,也可以選擇Snap Peak、Snap

46、 Trough或Do Not Snap,一般來說使用默認(rèn)選項(xiàng),在合成記錄道與井旁地震道標(biāo)定時會自動拾取相應(yīng)的波峰或波谷,若不滿意,可左鍵點(diǎn)擊Clear Picks按鈕,消除以上的拾取線。若滿意,可點(diǎn)擊Stretch按鈕,將合成記錄道和井旁地震道的對應(yīng)波形拉伸到相同的位置,同時井曲線得到了自動校正,即可看到此時的相關(guān)系數(shù)大小。也可以點(diǎn)擊圖3-64所示的Parameters打開互相關(guān)圖Cross Correlation,對相關(guān)的子波進(jìn)行質(zhì)控。首先,通過設(shè)置不同主頻的雷克子波進(jìn)行合成記錄道與井旁地震道的標(biāo)定。在此基礎(chǔ)上,提取一個新子波,子波提取過程參照前一節(jié)內(nèi)容,先利用統(tǒng)計(jì)法Statistical

47、求零相位子波而使振幅譜與地震數(shù)據(jù)的一樣,進(jìn)行拉伸操作,最后利用井曲線提取子波,進(jìn)一步提高合成記錄道與井旁地震道的相關(guān)程度。如圖3-73,在進(jìn)行Stretch操作時,默認(rèn)選項(xiàng)為Spline 漂移曲線間兩點(diǎn)進(jìn)行樣條內(nèi)差,可以改變內(nèi)差方式為線性 Linear 或多項(xiàng)式法Polynomial。聲波曲線改變項(xiàng)常選擇僅改變深-時曲線,以保持原始曲線值。圖3-73圖3-74如圖3-74,漂移曲線兩點(diǎn)間的內(nèi)差方法有三種:線性擬合(用直線段內(nèi)差);樣條擬合(用平滑曲線內(nèi)差);多項(xiàng)式擬合(利用最小平方擬合一條光滑曲線)。不同的內(nèi)差方式可能使聲波曲線的變化是相當(dāng)大的,一般選擇樣條擬合效果最佳。同時通過圖3-75的C

48、heck Shot Analysis判斷通過相關(guān)改進(jìn)的深-時曲線是否合理,若滿意,點(diǎn)擊Ok。然后點(diǎn)擊圖3-72中左下側(cè)的Ok按鈕,出現(xiàn)輸出校正的聲波曲線名稱(圖3-76),默認(rèn)為P-wave_corr,建議加上做相關(guān)操作的時間,便于后期的編輯或調(diào)用。其它的每一口井也按照該步驟實(shí)現(xiàn)最佳的井震相關(guān),從而完成所有井的相關(guān)工作,這是地震道初始模型建立以及地震反演的基礎(chǔ)。圖3-75圖3-76(五)模型建立每個地震道初始猜測模型是由阻抗曲線組成的,阻抗曲線通常是將實(shí)際的聲波曲線乘以密度曲線。阻抗曲線模型是以雙程旅行時來計(jì)算的,原始的測井曲線是按深度測量的,關(guān)鍵的一步是作深-時轉(zhuǎn)換,這在前面的子波提取和井震

49、相關(guān)中實(shí)現(xiàn)。初始猜測模型中使用井而無層位時只能產(chǎn)生均衡的水平模型;如果有一個拾取的層位約束,只能在每道進(jìn)行簡單井曲線的塊狀移動;拾取多個層位進(jìn)行約束,相當(dāng)于在每道都運(yùn)用可變的Check Shot校正,在兩個層位之間的阻抗進(jìn)行拉伸/壓縮。由此對于導(dǎo)入的地震解釋層位,通過剖面和平面圖形質(zhì)控,對層位進(jìn)行一定的平滑和編輯,使層位的拾取精度達(dá)到儲層反演的要求。圖3-77如圖3-77所示,點(diǎn)擊Model模塊,可以進(jìn)行模型的建立,模型參數(shù)的設(shè)定以及模型的拷貝、刪除等操作。點(diǎn)擊Build/Rebuild a Model進(jìn)入模型建立的設(shè)置窗口(圖3-78)。默認(rèn)的模型名稱為Model,可以自己修改以便后期的編輯

50、和調(diào)用。STRATA可以建立聲阻抗、橫波阻抗、彈性阻抗、疊前反演等模型。默認(rèn)為聲阻抗模型建立,點(diǎn)擊Next,選擇振幅數(shù)據(jù)源,默認(rèn)選擇Use Well Logs from Database,點(diǎn)擊Next,選擇工區(qū)內(nèi)要用的做好相關(guān)的井,點(diǎn)擊Next,進(jìn)入每口井P-wave和Density激活的選項(xiàng)(圖3-79),一般P-wave選用井震相關(guān)最佳的數(shù)據(jù),Density一般用原始數(shù)據(jù),對于沒有密度測井?dāng)?shù)據(jù)的井而言,系統(tǒng)默認(rèn)使用Gardner公式從速度進(jìn)行換算,如圖3-80,軟件中第二個選項(xiàng)提供了Lindseth公式從速度換算密度的方法,當(dāng)然還可以設(shè)為常數(shù)。圖3-78圖3-79圖3-80點(diǎn)擊Next,進(jìn)

51、入選擇解釋層位的界面,選中需要的約束層位即可。點(diǎn)擊Next,進(jìn)入模型道過濾選項(xiàng),如圖3-81,默認(rèn)自動在內(nèi)差后對模型通過高切濾波進(jìn)行平滑,使模型只含有低頻成分,這樣可以避免模型中的高頻成分干擾反演結(jié)果,反演結(jié)果更可靠??梢詸z查通過內(nèi)差模型道得到的高頻模型(選擇不進(jìn)行任何平滑處理)。圖3-81點(diǎn)擊Ok,初步建立了聲阻抗的初始猜測模型(圖3-77)。通過調(diào)節(jié)色標(biāo)使模型便于觀察,查看模型的縱橫向分布是否與測井以及區(qū)域地質(zhì)巖性的總體趨勢一致??梢酝ㄟ^點(diǎn)擊圖3-77中的Set Model Parameters進(jìn)行模型相關(guān)參數(shù)的調(diào)整。圖3-82如圖3-82所示,前面3個選項(xiàng)主要是工區(qū)、時間域信息以及測井曲

52、線的選擇。Advanced選項(xiàng)和圖3-81是一樣的。井間的內(nèi)插選項(xiàng)有3種:反距離加權(quán)(Inverse Distance Power)、三角測量(Triangulation)和克里金法(Kriging)。系統(tǒng)默認(rèn)選用反距離加權(quán)法進(jìn)行內(nèi)插,三角測量內(nèi)插是在相鄰井間擬合一系列的平面單元(圖3-83),克里金法需要有變量圖的輸入(圖3-84),一般都較少使用。圖3-83圖3-84圖3-85圖3-86如圖3-85所示,設(shè)置地層之間的接觸關(guān)系,默認(rèn)除第一層和最后一層外均為整合面,第一層為頂超面,最后一層為底超面,可以自定義不同的接觸關(guān)系,其對應(yīng)的內(nèi)插方式如圖3-86所示。對于地層尖滅的情形,如圖3-87,

53、可以通過使兩個層位聚合到一起來處理,如果層位交叉,每個層位被賦予優(yōu)先權(quán)值來處理。H1優(yōu)先級高,故H2層位超出部分被刪除;相反H2優(yōu)先級高,H1層位超出部分被刪除(圖3-88)。層位的優(yōu)先權(quán)值在圖3-89所示的Priority Value中進(jìn)行設(shè)定,值越大,優(yōu)先級越高。在Use Flag里還可以選擇層位是否在模型建立中使用。目前STRATA不支持逆斷層的建模,通過斷層面兩側(cè)的層位相連來模擬有逆斷層的阻抗模型。圖3-87圖3-88圖3-89圖3-90如圖3-90所示,最后一個Info項(xiàng)主要查看模型建立使用的子波,通過子波設(shè)置模塊指定后,選擇Force a rebuild of the model

54、on Apply后點(diǎn)擊Apply即可應(yīng)用新設(shè)定的子波,建立新的模型。通過各參數(shù)的設(shè)置、測試和更新,建立較合理的阻抗模型。(六)地震反演1、反演分析建立了阻抗模型,進(jìn)入地震反演環(huán)節(jié)。在做反演前最好進(jìn)行兩步反演分析,以便在后面的反演參數(shù)選擇中能選擇合適的參數(shù),從而得到更好的反演結(jié)果。反演分析的兩個內(nèi)容就是檢測井點(diǎn)位置實(shí)際阻抗和預(yù)測阻抗的誤差和確定最佳的縮放調(diào)整比例因子。反演分析操作在圖3-91所示的AnalysisPost-stack Analysis中進(jìn)行。點(diǎn)擊Model Based.進(jìn)入地震數(shù)據(jù)體的選擇,點(diǎn)擊Next,進(jìn)入阻抗模型的選擇、子波的設(shè)定以及處理時窗和采樣率的設(shè)置(圖3-92),點(diǎn)擊

55、Next選擇工區(qū)內(nèi)使用的井,點(diǎn)擊Next進(jìn)入復(fù)合道的提取方法,選擇默認(rèn)選項(xiàng)(圖3-93),點(diǎn)擊OK進(jìn)入圖3-94所示的反演分析窗口,點(diǎn)擊井名右側(cè)的箭頭可以前后瀏覽每一口井的實(shí)際阻抗和反演阻抗間的誤差以及合成記錄道和井旁地震道的相關(guān)系數(shù)。為更好的顯示測井的真實(shí)阻抗,我們將檢測輸入的地震數(shù)據(jù)體的頻率范圍,如圖3-95所示,在STRATA的Process菜單里選擇Utility中的第一項(xiàng)Amplitude Spectrum即為振幅譜的分析,點(diǎn)擊進(jìn)入圖3-96所示的選擇分析的地震道的范圍,點(diǎn)擊OK計(jì)算的振幅譜,如圖3-97所示,可以看到地震數(shù)據(jù)的相對高的頻率范圍。點(diǎn)擊圖3-94中的,進(jìn)入圖3-98打開

56、Curves在High pass和High cut頻率中分別填入75Hz和100Hz,并選中在誤差分析前對原始測井?dāng)?shù)據(jù)進(jìn)行帶通濾波,點(diǎn)擊Ok,井點(diǎn)實(shí)際阻抗和反演阻抗間的誤差如圖3-99所示,可以看到局部反演結(jié)果稍稍低于真實(shí)的阻抗值,如果這恰好是目的層中所關(guān)注的區(qū)域,而且其它的井也確認(rèn)出現(xiàn)這種情形,可以通過調(diào)整比例因子來進(jìn)行修正(圖3-100),進(jìn)一步提高反演的精度。當(dāng)然還可以在反演分析界面點(diǎn)擊Error AnalysisLogsError profile(圖3-101),顯示每口井的阻抗誤差值(圖3-102),從中可以很直觀的看到誤差相對較大的井,并檢查其相關(guān)的原因,從而控制反演的質(zhì)量。圖3-91圖3-92圖3-93圖3-94圖3-95圖3-96圖3-97圖3-98圖3-99圖3-100圖3-101圖3-102圖3-103反演分析結(jié)束,退出分析窗口。如圖3-103所示,在STRATA主界面點(diǎn)擊InversionPost-stack Inversion進(jìn)入可以啟動其中的各種疊后反演程序。2

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論