Strata軟件操作手冊(cè)-810資料

1、中國(guó)石油STRATA地震反演軟件常規(guī)操作手冊(cè)編寫人：付建元川慶物探公司物探研究中心二0一0年七月1目 錄一、地震反演原理及方法1（一）反演簡(jiǎn)介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（二）反演方法推薦702應(yīng)用STRATA軟件進(jìn)行常規(guī)地震反演Hampson-Russell公司自1987年起從事地球物理勘探方

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

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

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

5、在實(shí)現(xiàn)時(shí)可能會(huì)遇見許多麻煩。該軟件主要包含以下工作。1）生成合成記錄2）交互壓縮/拉伸測(cè)井曲線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ù) 地震子波，假設(shè)為不變的 隨機(jī)噪聲，與地震道不相關(guān)注意：1）在這個(gè)模型中，假設(shè)不存在多次波

6、；2）沒有考慮傳輸損失和幾何擴(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ù)。所有地震反演算法也都存在非唯一性問題?？赡艽嬖诙嘤谝环N地質(zhì)模型可以與地震數(shù)據(jù)相一致。要在這些可能的模型中決定一個(gè)模型，就需要去選擇地震數(shù)據(jù)以外的

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

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

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

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

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

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

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

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

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

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

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

18、礎(chǔ)上進(jìn)一步精細(xì)解釋，以滿足地震反演的層位約束的需要。對(duì)于井?dāng)?shù)據(jù)，除了最基本的井位信息，我們還需要常用的測(cè)井?dāng)?shù)據(jù)（包括聲波時(shí)差、密度、自然伽瑪、孔隙度等），測(cè)井分層以及測(cè)井解釋成果。對(duì)于已有鉆井的油氣水顯示、測(cè)試、試采等資料也需收集。對(duì)于斜井，其井斜（井軌跡）數(shù)據(jù)也是必需的。同樣以上數(shù)據(jù)均需通過一定的編輯形成符合軟件錄入的分列式ASCII文本數(shù)據(jù)或常用的標(biāo)準(zhǔn)LAS格式文件。當(dāng)然，前人對(duì)該地區(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)生的子波、地震體等。可以對(duì)以上數(shù)據(jù)進(jìn)行輸入、輸出、顯示、更改以及刪除等。對(duì)于地震數(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、

20、Analysis和Inversion模塊完成；反演結(jié)果的顯示及輸出主要在菜單File和View里進(jìn)行調(diào)整和操作。（二）數(shù)據(jù)輸入數(shù)據(jù)輸入主要包括與地震和測(cè)井相關(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ù)的幾何形狀（觀測(cè)系統(tǒng)）和文件格式的一系列界面。首先進(jìn)入圖3-7所示的選擇觀測(cè)系統(tǒng)類型，按實(shí)際工區(qū)類型選擇即可；點(diǎn)擊Next進(jìn)入圖3-8所示的選擇道頭中有無(wú)Inline&Xline號(hào)和X&Y坐標(biāo)，對(duì)于三維工區(qū)一般Inline&Xline號(hào)和X&Y坐標(biāo)均有，對(duì)于二維工區(qū)一般無(wú)Inline&Xline號(hào)，僅有X&Y坐標(biāo)；點(diǎn)擊Next進(jìn)入圖3-9的數(shù)據(jù)格式和測(cè)線道頭信息選擇的界面，一般數(shù)據(jù)格式多采用IBM浮點(diǎn)，少數(shù)采用IEEE浮點(diǎn)格式。下面的Inline

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

23、角有一個(gè)Display in new seismic window，建議勾選，地震剖面將在一個(gè)新的窗口顯示（圖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擇一個(gè)層位文件加載，也可以一次選擇所有的層位文件一次性加入，還可以選擇一個(gè)包含多個(gè)層位的文件實(shí)現(xiàn)一次性加載所有的層位。點(diǎn)擊OK進(jìn)圖3-14圖3-15如圖3-15所示的層位

24、文件格式的指定窗口，若層位是從Landmark按默認(rèn)格式輸出可直接選擇文件類型Default Landmark，一般的文本文件默認(rèn)選擇Free Format，文件里是否有多個(gè)層位選擇No，對(duì)于三維工區(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)確的填寫線道號(hào)及拾取層位的列數(shù)以及輸入層位數(shù)據(jù)跳過的行數(shù)，單位和層位值的范圍。最下面的為層位名稱及顯

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

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

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ù)）、測(cè)井分層、測(cè)井?dāng)?shù)據(jù)（包括聲波時(shí)差、密度、自然伽瑪、孔隙度等）信息，完成地震儲(chǔ)層反演還需要加載測(cè)井的解釋成果。如前所述，這些數(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)即可加載地震反演需要的測(cè)井曲線、分層及井斜數(shù)據(jù)。圖3-27圖3-28如圖3-28所示，Import Logs界面中選擇要加入的測(cè)井?dāng)?shù)據(jù)，下面有選擇測(cè)井文件格式的選項(xiàng)，一般收集的資料多為非標(biāo)準(zhǔn)的LAS文本文件，此時(shí)選擇General ASCII格式；數(shù)據(jù)類型分為三類（測(cè)井、分層、井斜），按選擇的文件選上對(duì)應(yīng)的數(shù)據(jù)類型。針對(duì)一個(gè)地震工區(qū)的整理好的所有測(cè)井?dāng)?shù)據(jù)，文件格式相同，對(duì)于井?dāng)?shù)不多且測(cè)井文件中測(cè)井曲線類型較少的情況可以一次性全選，完成同一

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

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

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

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

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

34、表示井的軌跡，一種是以井斜角Dip和方位角Azimuth來(lái)表示。在圖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ù)范圍，深度一般為測(cè)量深度。點(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可查看測(cè)井曲

35、線、分層及測(cè)井解釋成果，通過Summary可對(duì)測(cè)井?dāng)?shù)據(jù)進(jìn)行簡(jiǎ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)，子波可用振幅譜和它的相位譜來(lái)定義，點(diǎn)擊Frequenc

36、y標(biāo)簽可在限定的頻率域顯示子波（圖3-49），點(diǎn)擊History即可看到該子波的詳細(xì)信息。實(shí)際工作中常用雷克子波來(lái)制作合成記錄，且使用不同的主頻，通過點(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ù)和測(cè)井曲線估算的全子波或常相位子波（Use Well），最終獲得滿足地震反演所需的最佳子波。圖3-52圖3-53圖3-54統(tǒng)計(jì)子波的提取，不使用井，只對(duì)地震數(shù)據(jù)進(jìn)行自相關(guān)，假定相位是已知的。主要參數(shù)：道范圍（默認(rèn)是整個(gè)數(shù)據(jù)體，通?？稍O(shè)的大一些以增加統(tǒng)計(jì)運(yùn)算），如圖3-52所示。時(shí)間窗口（至少要為兩倍子波長(zhǎng)度），如圖3-53所示，可直接給定起始時(shí)間，也可以通過目的層頂?shù)咨舷卵由靵?lái)指定提取的時(shí)間窗口。子波長(zhǎng)度（默

38、認(rèn)值為200ms，一般可選100ms），如圖3-54所示，斜坡長(zhǎng)度一般用默認(rèn)值25ms。點(diǎn)擊OK開始運(yùn)算，完成即顯示所提取的統(tǒng)計(jì)子波（圖3-55），該子波為零相位是在提取子波時(shí)設(shè)置好的?？梢栽O(shè)定該子波用于井震相關(guān)分析。決定是否利用井來(lái)確定子波的相位的關(guān)鍵問題是井曲線是否與地震數(shù)據(jù)匹配。井的相關(guān)壓縮/拉伸會(huì)引起相位的問題，利用測(cè)井曲線估算子波的振幅譜和相位譜，估算的結(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)分析后的測(cè)井曲線。圖3-57所示，選擇地震體及離井的范圍，可選相鄰的道數(shù)或距離。圖3-58所示，選擇時(shí)間窗口。圖3-59所示，選擇子波長(zhǎng)度和提取方式（全子波、常相位和Roy White算法）。圖3-57圖3-58圖3-59在井曲線與地震的相關(guān)比較好時(shí)，全子波提取的振幅譜和相位譜在時(shí)間域的算子可以將測(cè)井曲線反射系數(shù)成形于地震復(fù)合道上，如圖3-60。圖3-60圖3-59中的第二種子波提取方式是常相位，這個(gè)選項(xiàng)是利用地震道自相關(guān)求取子波的振幅譜，

40、與地震統(tǒng)計(jì)法完全一樣，利用井曲線估算一個(gè)近似的常相位值（圖3-61），這種方式比全子波法更穩(wěn)健，特別是當(dāng)井曲線與地震的相關(guān)比較差。由此該軟件中將常相位設(shè)為默認(rèn)的井子波提取方式。圖3-61使用Roy White先生的算法，提取過程與全子波提取法類似，另外可形成一個(gè)診斷圖允許用戶估計(jì)子波提取的誤差（圖3-62）。圖3-62如圖3-60、61、62所示，如果井匹配比較好，以上3種方式求取的子波很類似。子波提取的一般問題是：利用井提取子波，首先要對(duì)井作適宜的相關(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-6

42、4，點(diǎn)擊Parameters打開互相關(guān)圖Cross Correlation顯示剩余相位誤差（圖3-65），當(dāng)旁瓣對(duì)稱時(shí)，顯示無(wú)剩余相位誤差。圖3-64圖3-65對(duì)于多口井時(shí)，可以用多井分析診斷，點(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)用手動(dòng)（交互圖形化）法將井曲線（深度域）與地震數(shù)據(jù)（時(shí)間域）進(jìn)行相關(guān)，對(duì)深-時(shí)曲線關(guān)系進(jìn)行人工校正，從而優(yōu)化初始模型和地震數(shù)據(jù)間的相關(guān)性。若有井的Check Shot信息，一個(gè)Check Shot表格是由一組深度點(diǎn)計(jì)算的實(shí)際雙程旅行時(shí)組成的。井震相關(guān)應(yīng)在Check Shot校正之后進(jìn)行，理想情況下井相關(guān)對(duì)井的改變是小的變化。井震相關(guān)改變深-時(shí)曲線的方式與Check Shot校正一樣。井震相關(guān)是在合成記錄上選擇幾個(gè)同相軸使之對(duì)應(yīng)于實(shí)際地震道上。由于使用合成記錄，子波的選擇是至關(guān)重要的，由此井震相關(guān)和子波提取是交互進(jìn)行的。如圖3-68，在

44、STRATA的左側(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)擊Co

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

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

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

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

49、井震相關(guān)中實(shí)現(xiàn)。初始猜測(cè)模型中使用井而無(wú)層位時(shí)只能產(chǎn)生均衡的水平模型；如果有一個(gè)拾取的層位約束，只能在每道進(jìn)行簡(jiǎn)單井曲線的塊狀移動(dòng)；拾取多個(gè)層位進(jìn)行約束，相當(dāng)于在每道都運(yùn)用可變的Check Shot校正，在兩個(gè)層位之間的阻抗進(jìn)行拉伸/壓縮。由此對(duì)于導(dǎo)入的地震解釋層位，通過剖面和平面圖形質(zhì)控，對(duì)層位進(jìn)行一定的平滑和編輯，使層位的拾取精度達(dá)到儲(chǔ)層反演的要求。圖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ù)，對(duì)于沒有密度測(cè)井?dāng)?shù)據(jù)的井而言，系統(tǒng)默認(rèn)使用Gardner公式從速度進(jìn)行換算，如圖3-80，軟件中第二個(gè)選項(xiàng)提供了Lindseth公式從速度換算密度的方法，當(dāng)然還可以設(shè)為常數(shù)。圖3-78圖3-79圖3-80點(diǎn)擊Next

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

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

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

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

55、點(diǎn)擊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ù)。為更好的顯示測(cè)井的真實(shí)阻抗，我們將檢測(cè)輸入的地震數(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ù)的相對(duì)高的頻率范圍。點(diǎn)擊圖3-94中的，進(jìn)入圖3-98

56、打開Curves在High pass和High cut頻率中分別填入75Hz和100Hz，并選中在誤差分析前對(duì)原始測(cè)井?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)整比例因子來(lái)進(jìn)行修正（圖3-100），進(jìn)一步提高反演的精度。當(dāng)然還可以在反演分析界面點(diǎn)擊Error AnalysisLogsError profile（圖3-101），顯示每口井的阻抗誤差值（圖3-102），從中可以很直觀的看到誤差相對(duì)較大的井，并檢查其相關(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)入可以啟