1、分頻 處理工作流程1 本章主要講述分頻處理的主要工作流程,以及解釋和分析分頻數(shù)據(jù)。我們相信經(jīng)過本章工作流程的描述，會使你更清晰的了解儲層的精細模型。 在你獲得清楚的認識之前，你需要明白分頻處理中的一些參數(shù)問題。對于一個新手來說，我們推薦你從第二步（二維模塊分析）開始。第二步從模型數(shù)據(jù)開始，演示各種參數(shù)對輸出效果的影響。這將有助于我們剛接觸的新手對分頻處理理論進一步的理解。 注意：流程中使用的一些數(shù)據(jù)在CD中，請注意版本說明如何去加載一般的數(shù)據(jù)。“real-data” 是Landmark訓練類數(shù)據(jù)體的命名樣式。 通過模型的操作，初學者會看到疊加道和多樣化地形地層是如何干擾影響分頻中的離散F變的。

2、在所有的流程中，你會知道如何去拿捏這些參數(shù)，諸如一個單形體的時窗的選擇，都會大大的影響分頻處理的結(jié)果。 當你對模型數(shù)據(jù)的處理感到很愜意的時候，你就有能力進行你的實際勘探數(shù)據(jù)的參數(shù)化，并很快得跑完整個流程。 具體的流程如下： 數(shù)據(jù)準備初測調(diào)諧二維模塊三維流程時窗目標調(diào)諧圖離散頻率體儲層厚度估算 2數(shù)據(jù)準備（Data Preparation）3 輸入的地震數(shù)據(jù)體須經(jīng)精細的處理，去建立正確的子波相位和真實的道間振幅信息。一個定態(tài)的子波對分頻處理是至關(guān)重要的。地表一致性處理和二維的濾波去噪，會加強你的分頻效果。輸入數(shù)據(jù)的帶通濾波器也許會消除一些地質(zhì)信息，而這部分信息有的卻是你準備去分解的。如果可能盡量

3、用沒有未濾波的偏移數(shù)據(jù)作為輸入去處理。雖然分頻處理是在振幅和相位域操作的，但是最好選擇16道32位的數(shù)據(jù)，因為8位的數(shù)據(jù)不能做細節(jié)分析，雖然可以做初步的調(diào)諧工作。 在具體處理之前，你需作到所選擇的層必須經(jīng)過內(nèi)插而且要貫穿所有斷層。頻譜分析通常能使斷層更加明了，這使你能重新去認識這些斷層，在頻譜成像上。如果你的層位具有相對高的信噪比和穩(wěn)定的地震波同相軸， 還得考慮對數(shù)據(jù)作剩余靜校正。 這一步將清除許多由于工作站上自動追蹤的同相軸所引起的跳躍現(xiàn)象。你也可以通過簡單的圓滑地層來處理。頻譜分析前的這步會敏銳地提高成像質(zhì)量。這種提高能使你恰好得到地質(zhì)信息的一個重要部分。 輸入數(shù)據(jù)須適當?shù)钠?，分頻處理假

4、定已經(jīng)校正傾斜地層和菲涅耳帶的影響。簡單的入射模型能夠決定你的砂體是否發(fā)育在目的層段。假設某個砂體可以成像，然而，不能說她就是有形的。對于各種地震技術(shù)都顯示的砂體，她必然有一個突出的波阻抗差異。從聲波曲線的分析到合成地震記錄形成，你可以測定你的砂體是不是有形的。 這一步回答了這個問題：你的數(shù)據(jù)準備好了嗎？4初測調(diào)諧（Reconnaissance Tuning Cube）5 在單位厚度、不連續(xù)疊加的層和斷層上，最優(yōu)參數(shù)化的調(diào)諧體擁有最大的幾率顯示出露變化。細微的道間變化無論在厚度或者疊層上，都會使目的層段的數(shù)據(jù)元素調(diào)協(xié)成不同的分立的頻率成分。高頻部分顯示細微的、低分辨率的斷層；低頻部分則趨向于去

5、噪和展示宏觀結(jié)構(gòu)和斷層趨勢。 這步流程給予一個目的層的調(diào)諧的過程的瀏覽，她將回答如下問題： 你看到目的層段了嗎？ 流程主要分三部分： 一、Parameterization（參數(shù)化） 使用一個完全的長時窗（100200ms），中心在目的層段層上，有10到20%的余弦遞變。在確定可以獲得頻率從DC到Nyquists 后，使用一個寬泛的頻率范圍讓你去限定在解釋種有用的信號帶寬。 二、Animation（制作） 通過啟動制作調(diào)諧體，觀察并最終決定帶寬；也就決定這相干的地質(zhì)信息能否被觀測到。 三、Selection（選擇） 選擇一對相對獨立的頻率體成像顯示當前的目的層地質(zhì)地形。 如果在這個流程中沒有相干

6、的地質(zhì)信息被觀測到，那么接著，SpecDecomp工具就不能被用來加深你對目的層層位的解釋/理解。如果這個流程KO了，你就能夠最優(yōu)化你的調(diào)諧參數(shù)，作進一步的詳細分析。 6快速瀏覽流程一 1.在啟動OpenWorks中選擇Application(模塊)SpecDecomp(分頻)Tuning Cube (調(diào)諧體) 2.選擇輸入的項目名稱如tns_train。 3.選擇輸入的數(shù)據(jù)體到seimic.cmp(CMP是由EJB容器自動完成對數(shù)據(jù)庫的操作，你所有做的，就是在實體bean重寫入SetXXX或getXXX方法，然后在ejb-jar.xml中定義cmp-field 4.設置分析時窗基準層到實際層

7、位。 5.改變cosine（余弦）梯度有20%左右的遞變。 6.給予數(shù)據(jù)體前綴名。 注意輸出的幅度體被命名位Workflow_1_0_DFT_AMP.bri. DFT(離散F變)表征著變換的方法；Amp 表示體屬性為振幅體。SpecDecomp輸出的地震數(shù)據(jù)體都被寫成程序塊格式。77.默認其他參數(shù)并選擇需計算數(shù)據(jù)體。計算會給予你5分多鐘的休息時間，具體在于你機器的性能。如果你的機器不是那么快，你可以去小憩一會。88.確定計算完成計算每道DFT過濾、標準、歸類頻率切片制作： 現(xiàn)在你預備去啟動地震波剖面，選擇一些分立的頻率對目的層段地質(zhì)地形成像。1.運行SeisWorks 3D,選擇tns_tra

8、in 項目。2.打開SD_wf_1-0-1.t.w3s 時域。3.選擇解釋人員（程序），不要井和斷層數(shù)據(jù)。4.從剖面上看，選擇 TaskSeismic 打開一個新地震剖面顯示窗口。95.改變輸入地震數(shù)據(jù)設置為Workflow_1_0_DFT_Amp.bri,并重新選擇時間頻率至0到250Hz 。頻率諧和體剖面圖實際工作中上圖是很難被確切解釋的，然而你能注意到在35Hz左右有很多大的幅度差異現(xiàn)象。注意：在（seismic view）地震波剖面圖上，35Hz顯示的是35ms同時可以看到，70Hz以下的數(shù)據(jù)都顯得像噪聲。106.從新的地震剖面窗口，選擇 SeismicSelect from List

9、。并選擇T0的時間切片。這實際上是0Hz的頻率切片。7.改變（color bar）色標到blkwht（自行調(diào)節(jié)）。8.選擇ViewFrame Control,在seismic view工具條中。9.單擊Frame Control forward 按紐圖標“”啟動制作調(diào)諧體。1110. 制作停止在30Hz你可以看到一個清楚的河道軌跡。當你在這種類型的顯示圖上看到了非地質(zhì)或者是很線性化的信息時，你的原始數(shù)據(jù)可能有點問題。那樣的話，你必須檢測你的輸入層位是否完全內(nèi)差并覆蓋目的層段；當然你還需要去確認你的地震數(shù)據(jù)是不是符合數(shù)據(jù)準備流程的要求。注意：Motion （引導線）是一個非常強有力的工具，解釋人

10、員可以用來在地震勘探中提取數(shù)值。你可以啟動制作整個數(shù)據(jù)體來探測細微的變化，這些變化在靜態(tài)地震剖面圖上是很難被檢測到的。在頻譜分析中，你可以發(fā)揮motion的作用去解釋儲層各向異性和厚度的復雜變化。1211. 向前并復查從0到折疊頻率的數(shù)據(jù)體制作。如果你想去控制制作速度，你需要在EarthCube 或者 OpenVision中去設置。如果數(shù)據(jù)的通頻帶寬在60Hz左右，你就基本上得不到60Hz以外的信息。當你操作你自己的未濾波的數(shù)據(jù)時，你會發(fā)現(xiàn)有用信息的頻率是從0到折疊頻率的。13選取從上面的分析中，你會意識到頻率從10到40Hz包含了有價值的地質(zhì)信息，也就是你想進一步作研究的這部分。14為了充分

11、理解分析這些調(diào)諧體的成像，你必須明白砂體越厚，振幅譜的V形切痕（凹口）就越多。注意： 振幅譜上，楔形模型的最厚部分有12個切口，最薄部分僅就一個。如果你通過楔形模型（從低頻到高頻）來制作調(diào)諧體，你會發(fā)現(xiàn)幅度等值線是從厚到薄。在下章中將具體分析楔形模型。你能注意到等值線從河道中心漂移嗎，Page34 當你在從10到40Hz制作的時候？ 15二維模塊（2D Scope-out Analysis）16 一旦你已經(jīng)通過初測調(diào)諧，你可以最優(yōu)化頻譜分析參數(shù)，這步解決這個問題： 時窗和（taper）梯度的關(guān)系？ 最優(yōu)參數(shù)化流程的目的是為你提供一些參數(shù)，這些參數(shù)是用作覆蓋整個區(qū)塊的三維時窗分析流程。這步流程還

12、用來作訓練流程，為了大家更好地了解頻譜分析的細枝末節(jié)。這一步一般只能作橫截面上一個獨立單元的分析。二維斷面測試方法效果不錯，因為時窗的移動范圍小，并且DFTs的迭代算法可以很快，還可以在原地震剖面上看到是否有對地質(zhì)現(xiàn)象突出的頻率響應。 流程分二步： .選擇一個橫截面作分析 .使用一個遞變范圍和分析時窗的長度來生成和觀測頻譜分析數(shù)據(jù) 盡管流程是分為二步的，我們將指導您通過一個楔形模型分形的詳細分析來加強你對頻譜分析的理解。17數(shù)據(jù)質(zhì)量控制 先提供一個模型數(shù)據(jù)的基本情況，楔形模型在下面相應地作出分析。 瞬時厚度1.從一個地震剖面窗口，改變輸入地震數(shù)據(jù)設置為wedge 3dv。 2.從 地震剖面窗口

13、，選擇SeismicSelect from List. 并選擇線號1131，重新選擇時窗 900到1150ms。 2.你現(xiàn)在可以看到這個楔形的模型經(jīng)過濾波整形 ，注意到瞬時厚度怎么沒有楔形模型上那么尖呢，你也許想改變一下顯示的比例。 183.高亮顯示the Top_Wedge_Interpreted 和 base_Wedge_Interpreted 地層。 注意為什么楔形解釋的精度不能和實際的模型相匹配呢？地下地形因薄層調(diào)諧而發(fā)生扭曲。如果你把沒有帶限的地震楔型變換到頻率域，你就能看到如下所示的圖象。注意振幅譜上切口會隨著厚度的增加而增多。這是由于地層厚度大余1/4時引起的干擾和破壞，高頻成分

14、短波長。也要看到?jīng)]有頻率能使0厚度的地層成像。（不能分辨尖滅現(xiàn)象）19上面圖中的白色條帶代表最主要的頻率。在下一節(jié)中你會建立一個相似的振幅譜分析圖，第三步在帶限的地震剖面上。 初始調(diào)諧分析本節(jié)指導我們在帶限楔型上去貫通一個調(diào)諧體的過程。 1.選擇Applications SpecDecompTuning Cube 從 Openworks 平臺上。 2.選擇tsp_train 項目。 3.設置輸入體為 wedge 3dv。 4.把開始結(jié)束線號調(diào)到1131。 5.改變分析時窗間斷型到Constant Window(連續(xù)型). 6.設置基準時間為1025ms。 7.選擇最大時窗長度，使用MB3 在時

15、窗長度輸入欄中。注意你可以使用MB3在輸入各項，可以得到可用輸入值的列表。 8.最大頻率改為250. 9.將余弦梯度給定為20%。 20 10.為Workflow_2_1設置數(shù)據(jù)體前綴. 11.選擇Output Horizons. 12.選擇Autocreate (自動生成) 如下層位。 13.完成后關(guān)閉輸出設置。2114.默認其他參數(shù)選擇Compute Volume。計算數(shù)據(jù)體15.關(guān)掉Computations Completed 信息，這步一分鐘就搞定了。16.顯示 1131線 在地震剖面圖上Workflow_2_1_DFT_Amp.bri.重新設置列表作新的brick文件。2217.你需

16、要再選擇一次時窗范圍在地震剖面圖上（從0到250Hz）注意：為什么通頻帶下面的所有數(shù)據(jù)都成了噪聲。也還要看到噪聲在72Hz左右每10Hz成一條帶，你能想到是什么引起的嗎？也許是歸一化引起這些條帶，因為這個問題出現(xiàn)在模型的1560到1570道之間。2318.地層演示（1025_ConstantTime_DFT_FisrstPeak_Freq）憑直覺你能看到已調(diào)諧地層的高振幅部分，得注意在此圖上色標是反轉(zhuǎn)的。是不是FirstPeak_Frep 層看來就像是沿著首峰的軌跡而來的？移動光標沿著并觀看顯示圖底部振幅值的輸出，你會看到一些振幅的峰值在計算出的FirstPeak_Frep 層之上。為什么呢？

17、首峰的計算是正確的，問題在于你的顯示是在歸一化數(shù)據(jù)上的。因此說明FirstPeak_Frep 層是在歸一化前計算出來的。歸一化后，所有的頻率切片被提高到同一個幅度水平。歸一化的細節(jié)具體討論在P68，首峰計算在P136。24為了分析數(shù)據(jù)，研究楔型體的瞬時厚度(t) 是如何決定振幅譜伴隨頻率的削減周期(Pf )就很重要。同樣的，你需要知道頻率成分（f）的數(shù)值是如何決定振幅譜伴隨地層厚度的削減周期(Pt )。在你作任何有關(guān)厚度估算流程前，必須充分理解這些關(guān)系。25調(diào)諧厚度的數(shù)學表達式如下：舉例來說，一個10Hz的子波可能有50ms的調(diào)諧厚度：20Hz的子波的調(diào)諧厚度為25ms：30Hz的子波的調(diào)諧厚

18、度為17ms：注意： Widess(1973)和Wood(1982)論述了一個求取調(diào)諧厚度的方法表達式：2610，20和30Hz的調(diào)諧厚度（前面計算的），其帶限振幅譜和頻率域圖示圖下。你需要確定你明白這2條振幅和頻率示意曲線的相互關(guān)系。如30Hz時，0振幅在2條曲線上的交叉點在哪兒？（34ms）27你也許會得到這樣的結(jié)論：如果我簡單的反演最主要的頻率，就會得到相對的調(diào)諧厚度。那就是說，主頻越高，地層越薄。對于初步的分析，你可以使用地震層位估算并反演xxx_DFT_FirstPeak_Freq 層，就可以獲得一個相對的調(diào)諧厚度圖。流程6.1指導你進行這個計算的細節(jié)。然而，你必須保持這樣的觀點：2

19、0ms的瞬時厚度，相對于更新世沉積的砂體，在8000f/s圖示為80ft，當20ms相對泥盆紀石灰?guī)r時，在18000f/s圖示為180ft。就像我們所知道的，地質(zhì)學會變得很容易，假使我們所打的井是時間意義上的而不是深度的時候。另一個分析得出的結(jié)論是：你花在主頻分析上的時間越長，你的解就越向惟一靠攏。28二維時窗分析在區(qū)塊輸出參數(shù)測試中，我們推薦你從觀測二維線的DFT剖面開始。你首先要通過使用一個窄的窗口（或者，你可以使用梯度減小有當前的時窗）聚焦一個小的時窗。縮短時窗或拉長梯度對DFTs 的輸出效果是相似的。其效果就是DFTs 在垂向上（在頻率的方向上）明顯得圓滑了。圓滑后，能很容易地清除缺切

20、并高度了聯(lián)合振幅譜分析，因而，一些目的層段的地質(zhì)地形在35Hz下的，可以顯示在100ms窗口中，當在160ms時窗中觀測時。（P2-50）因此，我們建議“window carpentry” 窗口工序化，在選擇任何分立的頻率作進一步研究之前。1.在地震工作中，顯示W(wǎng)orkflow_2_1_DFT_Top, Workflow_2_1_DFT_Middle, Workflow_2_1_DFT_Bottom分析窗口在wedge01.3dv seismic view 上。29302.在SpecDecomp-Runing Cube中，改變窗口大小從218到100ms。 這是小時窗的最好輸出。如果小時窗不能

21、顯示任何地質(zhì)信息，則逐漸增大其大小直到你看到相干的同相軸為止。接著，采準這個時窗直到很多冗余信息的出現(xiàn)。 如果你有時間，你可以從30ms時窗開始，按每次增加10ms一直增加到200ms。3.關(guān)閉并歸一化切片。4.默認其他參數(shù)開始計算，按Compute Volume.如果你選擇輸出的數(shù)據(jù)體是已經(jīng)存在的，你需準備是否刪除或者覆蓋原來的體。5.選擇覆蓋/合成文件。當你測試/分析你的調(diào)諧體的時候,這個覆蓋和刪除按紐使你能夠使用規(guī)范化易于管理的輸出數(shù)據(jù)設置。警告： 雖然覆蓋和刪除按紐是很有用途的,你使用的時候還是必須很小心.尤其是,當你不能夠再生成特殊的切片(你正在地震工作中觀測)的時候，刪除足以引起問

22、題。請參考日志文件去解決已有的問題。316.演示/刷新 加載在地震工作中的數(shù)據(jù)體。你會發(fā)現(xiàn)默認的顯示比例對沒有歸一化的數(shù)據(jù)是不能正確顯示的。7.從地震演示上，選擇Seismic Parameters 變化Dynamic Scale&Clip （動態(tài)比例按紐）從No Scaling 到 Sample*15/1。這些與218ms分析的相比會怎么樣呢？它會不會更圓滑？首峰地層看起來象不象確切地計算出來的？地層噪聲條帶時夠已經(jīng)被清楚？你能想到在左邊的演示中當楔型mod變得非常薄的時候會發(fā)生什么變化嗎？328.重復操作生成一個50ms的分析窗口。這些與218和100ms的時窗分析比較又會怎么樣呢？請把注

23、意力集中到第一條譜帶密集區(qū)的凹口處，當凹口特別明顯的時候，你已經(jīng)對你的時窗分析參數(shù)做到了最優(yōu)化。9.使用這個程序去測試不同的遞變項。10.繼續(xù)測試直到你對所設置的參數(shù)滿意為止。33二維教學模型 下面的實例在SpecDecomp的安裝盤上。這些實物是用來讓你測試不同的參數(shù)如何反應不同的地質(zhì)條件。請查看發(fā)布的日志文件來加載這些模型。 wedge01.3dv 是在前面練習中使用的楔形模型。這個模型涵蓋了40條線（1111到1150）和460道（1110到1561）。時窗范圍是900到1150ms，而頻率范圍是從10-20Hz（1111線）到10-100Hz（1150線）。wedge01.3dv對于w

24、edge01.3dv，嘗試生成并闡明下面的展示圖。34 chanv201.3dv 是從1101到1229線的一個河道模型，時窗2800到3200ms。這個模型從薄到厚隨著線號的增加。chanv201.3dv chanv401.3dv 是一個河道疊加模型，就一條線1195。是chanv201.3dv的3個拷貝的混成，時窗從2900到3100ms。所有道都有不同的振幅。chanv401.3dv35 thinch01.3dv 是一個抽道集河道模型從1101道1229線，時窗是2800到3200ms。它和chanv201.3dv的道集模型厚度遞減是一樣，所以一般看來都是薄的。thinch01.3dv

25、chx2v101.3dv 是一個由chanv201.3dv和thinch01.3dv疊加的河道模型。從1101道1229線，時窗是2800到3200ms。chx2v101.3dv36 chx2v201.3dv 是兩個淺河道模型的集合，線號為1196，時窗2900到3100ms。 chx2v201.3dv thinv201.3dv thinv201.3dv是淺河道的振幅變化模型，但沒有薄層！線道號是1175和1420，時窗2800到3200ms。37 thinv301.3dv 淺河道模型，從thinch01.3dv，上面有一個強同相軸。線道號是1175和1420，時窗2800到3200ms。 t

26、hinmd01.3dv 是尖峰模型用來形成淺河道模型，線道號是1175和1420，時窗2800到3200ms。38 fubar401.3dv 是.wedge01.3dv 1130線的靜態(tài)頻移，相位旋轉(zhuǎn)，振幅衰弱的變型。這個模型線號1130，時窗為900到1150ms。 除fubar401.3dv 和wedge01.3dv 是1ms采樣，其他的.3dv都是2ms采樣。39三維時窗（3D Running-window Analysis）40 對實際數(shù)據(jù)來說，與模型數(shù)據(jù)是有很多不同的，大多數(shù)勘探需要在可變分析時窗上進行3D分析。這步流程主要解決以下問題： “對多重疊加的產(chǎn)層怎么處理？”“如何與地質(zhì)現(xiàn)

27、象同步？”“能觀測到泛濫平原的分層面嗎？” “信道被時間約束的地方在哪兒？”“能否看到其他的地震相界面？”“時窗分析的最佳位置點在什么？” 我們在開眼界之前，就可以看到流程帶來的地形地貌。此流程有3個主要步驟。 Selecting(選取) 選擇一個區(qū)帶并在研究區(qū)上下適度地延伸一段距離。如果你的工區(qū)位于一個傾斜的地層，在生成運行時窗分析前需要把這個層拉平。 Generating(生成) 生成運行時窗分析，頻率成分從初測調(diào)諧體中選取，覆蓋全區(qū)。與初測調(diào)諧體相比，這里對每個時間點的局部分析是很重要的。一個80ms的時窗使用高斯分級可以提供所需的局部化。 Animating(演化) 通過得到的數(shù)據(jù)體演

28、化測定： 區(qū)帶間的信號滲濾（例如在已知現(xiàn)象作用前，你移動窗口的時距不再顯示在圖象中。）41工區(qū)是由實際的層組成的，層是需要拉平的，為了消除分析中傾斜干擾的可能性。在這個流程中，你僅分析那些顯現(xiàn)最高的信噪比的頻率。視具體情況使用10，20和30Hz，和初測調(diào)諧體定義的一樣（流程一）。選擇ToolsVolume Recon 從調(diào)諧體菜單欄。這個項目會被自動設置為當前的項目。如果必要的話，設置輸入數(shù)據(jù)體為real_data.cmp。在時窗范圍選項中，改變時窗范圍為Flatten To Horizon。對于層位，點擊List. 按紐，選擇實際層位。默認輸出采樣頻率為2ms。改變中心窗口定義層的時間長度

29、為200ms。生成初測調(diào)諧體在程序選項中，設置時窗為100ms。設置變換函數(shù)為Gaussian。高斯變換集中于窗口中心，P120 有各種變換的對比。設置體前綴為WF_3_1_Gaussian_100。體前綴中“100”標志著時窗長度?！癎aussian”標志變換類型。42設置頻率列表：10，20，30。默認其他參數(shù)并選擇Compute Volume。選擇OK到Computations Create Warning. 計算將花費幾個小時，就看你工作站的CPU，內(nèi)存和硬盤傳輸速度。本工區(qū)在800MHz Sun-Blade-1000 工作站上花費了40分鐘 。清除完成信息。43演化初測調(diào)諧體加載 W

30、F_3_1_Gaussian_100_Amp_030_Hz_DFT.bri 數(shù)據(jù)體到地震時間切片觀測圖上。你需要重新選擇地震時間，調(diào)諧體才會被顯示出來。 點擊“”演化整個分立的能量體。 分立的能量體在30Hz時的窗口。就像你在地質(zhì)時期上下旅行的演化一樣，0ms表示輸出體實際層的新位置點。嘗試用各種色標；blkwht_color_max,Sopacity2, 和Sspecdecomp。你看到有趣的地方了嗎？請注意邊灘（1260線，1170道附近，它什么時候出現(xiàn)和消失的呢？你又怎么解釋這種現(xiàn)象呢？44演化10到20Hz的調(diào)諧體。當你演化整個“運行時窗分析能量體”時，你須要能夠分選出在什么樣的級數(shù)上

31、最能涵蓋你的目的層段。你變更的時窗在什么時候不能使你看到主要目標的效果？在任何變更范圍內(nèi)底圖和切片的調(diào)節(jié)會降低分辨率嗎？你演化的時候，你可以去闡明其地質(zhì)時期發(fā)生的變化。例如，你可以看到由主河道發(fā)展的分叉河道的現(xiàn)象嗎？你可以看到邊灘在地質(zhì)年代上的發(fā)展嗎？當你通過地質(zhì)時期去演化時，你也能看到河道充填在其上的變化？當你演化一個河道的時候，你首先看到一個暗黑色條帶在河道下面的中心上。然后，暗黑色條帶相稱性地向河道兩邊轉(zhuǎn)動。河道沉積充填簡單模式檢索Geology 101會發(fā)現(xiàn)大多數(shù)的河道種砂體/沉積物水流較緩的邊緣。因此，砂體一般都是中間薄由于河流的快速流動。能量體的演化于調(diào)諧體不一樣，你看到同樣厚度的

32、河道轉(zhuǎn)化為不同的頻率切片。要緊記調(diào)諧體的演化在時間上是凍結(jié)的（窗口一直不動）。運行窗口種分立的能量體演化凍結(jié)的只是在單一的頻率上，而且它的諧和圖可以在截面上上下移動。更重要的是，當觀測復雜儲層的成像時，地質(zhì)學家，地球物理學家和儲量工作者以及現(xiàn)場工程師可以看到在各自獨立的學科上的直接效應。地震勘探隊伍就能夠盡快而協(xié)同地影響到全局。45當你有時間的時候，可以重復整個流程，測試不同的時窗的效果。當你有時間的時候，可以重復整個流程，測試不同的變換梯度。在本流程的最后，你也許對高保真度三維頻譜成像（通過分頻處理）已經(jīng)有了一個粗略的概念，三維高保真度來自于直接對反射系數(shù)的成像（如子波的影響被清除），而且抽

33、樣的反射系數(shù)包是獨立于傾角的。這個流程輸出的一些參數(shù)是為時窗分析和Taper設置的。下面列舉一些合理的參數(shù)：.Analysis Window Reference Horizon real_data horizon.Reference Horizon Time Shift 0ms.Window Length 100 ms.Taper 20% Cosine Taper 時窗分析 參考層位 實際層位參考層位時窗變更0ms時窗長度100ms變換20%余弦遞變未拉平當進行分析的時候，你沒有必要一定拉平層數(shù)據(jù)。然而，拉平有兩個主要的優(yōu)點： 可以快速生成小數(shù)據(jù)體。 可以在拉平的地質(zhì)年代上演化切片。 不拉平時

34、跑完整個流程，你首先要運行Volume Recon 從100ms到730ms。因為輸入層的范圍是200到730ms而你還需要一個上下100ms的延伸。然后你可以在Seismic View種選擇Horizons Flatten.交互式得拉平你輸出的橫截面。不幸的是，地震工作不允許在時間切片上層拉平。因此，你需要在演化你的時間切片前，使用疊后資料拉平輸出的數(shù)據(jù)體。在這個流程中歸一化是不一樣的，相鄰未拉平道集會影響歸一化的數(shù)量。（P65）46交互式拉平使用拉平數(shù)據(jù)的一個問題是：它使拉平數(shù)據(jù)和未拉平數(shù)據(jù)對比起來非常困難。如何去克服，一個辦法是在地震轉(zhuǎn)換中使用持續(xù)拉平，當然這種轉(zhuǎn)換要求能同時顯示拉平和未

35、拉平的地震數(shù)據(jù)體。這個流程發(fā)表在July/August 2002 UserNet article title:Spectral Decomposition Persistent Flattening Workflow for Volume Recon 你可以從兩方面獲得。這篇文章的另一個Adpbe Acrobat版本包含在$SPECDHOME/docs 文件夾里叫SpecD_PersistentFlattening_Workflow.pdf.你可以打開Adpbe Acrobat查看這個文件。47目標調(diào)諧圖（Targeted Tuning Map Analysis ）48目標調(diào)諧圖分析 這步也叫

36、做目標調(diào)諧體。此流程假定頻率范圍，時窗分析和變化都已經(jīng)定義好了。這個流程能回答如下問題：“哪兒才是鉆井的亮點？”流程輸出包括：. 調(diào)諧體和子波疊加分離.用疊加補償性頻譜化的調(diào)諧體制作全帶寬的振幅圖（在子波被分離之后）信號帶寬綜合部分經(jīng)?？梢陨筛玫母爬ㄐ缘臄鄬蛹捌溥吔鐖D。不同的斷層和邊界用不同的頻率成分（一般高頻時比較清晰）來成像?？偤停ǒB加）調(diào)諧體（在清除疊加的子波后）允許每個頻率成分給予一個權(quán)重，因而當強化邊界和斷層的時候，由于厚度的可變性，會導致疊加出的調(diào)諧效應。然而，這存在一個冒險，這種疊加中你可能會抹掉一些頻譜異常。如果你要成圖的地形恰好是這樣的，調(diào)入的大于這一半，而調(diào)出的大于另一

37、半，疊加就會嚴重的消除這些異常。在這種情況下，單獨的窄帶，調(diào)諧圖可能會提供優(yōu)秀的圖象。.分析時窗層位的上、中和下.（First）頻譜最大層.由工區(qū)中波長和多邊形形成的瞬時厚度圖.主要切片“調(diào)諧圖” 保存為地震屬性層位 49基礎：清除子波疊加 調(diào)諧體由3部分組成：薄層干涉，子波疊加和噪聲。既然地質(zhì)響應是解釋人員最感興趣的部分，所以人們都很謹慎在沒有降低地質(zhì)信息的情況下來抵消子波的幅度。這一步將調(diào)諧體縮減為薄層干涉，噪音。50通用頻譜均衡技術(shù)依靠稀疏不變量穩(wěn)態(tài)統(tǒng)計學。那就是，先為所有點假定一個連續(xù)的子波。如果在拉平的層上呈現(xiàn)出地質(zhì)調(diào)諧多樣化，然后你可以通過每個頻率切片的中間振幅值作均衡補償來抵消子

38、波波譜。歸一譜白化以后，子波疊加效應縮到最小程度，調(diào)諧體保留兩個部分：薄層干涉，噪音。51 在調(diào)諧圖上（頻率切片）地層中，薄層干涉是作為相干振幅的變量。隨機噪聲斑點化干涉的樣式就像低質(zhì)量的雪花點電視帶來的感覺一樣。在主頻上，相對高信噪比（SNR）帶來了相對清晰的薄層調(diào)諧的圖片。在有用帶寬外的頻率上，低信噪比則帶來了雪花電視。生成目標性調(diào)諧體 (流程4.0) 這步和初測調(diào)諧體的生成(流程一)是相似的，那是確定你的研究工區(qū)是否具有有利的地質(zhì)信息。然而，本流程目標調(diào)諧體在運行時，是有最優(yōu)的參數(shù)輸入的。如果你對3D時窗分析進程的分析流程三所依助的參數(shù)比流程一中的快速瀏覽差異較多的話，你就必須要再運行流

39、程1.0，使用更新過的參數(shù)。 生成目標性調(diào)諧圖 (流程4.1) 一旦生成了一個正常的調(diào)諧體，就伴隨著一些細節(jié)分析，應用于各個頻率切片上。儲層厚度估計就是這么一個流程，可以在頻率切片上顯示。然而，這個流程，需要把精力集中在調(diào)諧圖對地震解釋層位的轉(zhuǎn)換以及精確的頻率圖計算上。 選擇Application SpecDecomp Tuning Mapper 從藍瑪上啟動。 選擇spt_train 項目。 設置輸入的數(shù)據(jù)體為Workflow_1_0_DFT_Amp.bri。 設置切片范圍在Range List上為10，20，30，40。 在Merge Slices 下，打開Sum并生成一個輸出層位叫作WF

40、_1.0_Merged_Frequencies. 設置輸出層位前綴未“WF_1.0” “WF_1.0”表示輸入體的來源。52點擊Compute Horizons from Slices。點擊OK,對Horizon Create Warning。清除 the computations completed 信息。53在地震解釋地圖上顯示W(wǎng)F_1.0_030Hz層。 注意看清晰的河道特征。 54調(diào)整色標，只顯示河道。 55調(diào)諧成圖解釋 (流程4.2) 作為一個解釋人員，一個主要的任務就是生成圖件并計算圈閉面積。藍瑪?shù)腁uto-Polygon Generator 工具就是用來成圖和面積計算的，而且是更

41、加有效更為精確的。 在SeisWorks Map View上，選擇Horizons Auto-Polygon Generator. 這個工具將：計算所有polygon的面積；生成適合在EarthCube 或 OpenVision 中觀測的3D polygons；生成適合在Z-MAPPlus中成圖的區(qū)域?qū)φ眨刂疲﹫D。 參數(shù)化的Auto-Polygon Generator 顯示如右圖：56點擊Process Horizon 后，你會等到如下（左圖）的Completed信息。 觀察生成的Polygons的數(shù)量，以及最大的20個的面積。 再運行一次Auto-Polygon Generator 作業(yè)，

42、將 Minimum Polygon Area (acres) to 100 。新的作業(yè)將只計算出7個Polygons。57在SeisWorks Map view 中顯示這些Polygons及其面積。確認并關(guān)閉指示層。 58向前，關(guān)閉Auto-Polygon Generator 對話框，在地震橫截面，你可以開始審閱河道部分如何變化的情況。在Seismic View 顯示real_data.bri的1175線。在Seismic View 只打開Polygons層。這將顯示河道的全局范圍。如果觀看Polygons層對你來說有困難，請選擇View Contents 并選上 Horizon Inters

43、ection Circles?，F(xiàn)在顯示1475線。注意：在3D中查看Polygons ，你需要進入EarthCube 和OpenVision。流程6.1 在EarthCube2003.3 Help Workflows 下描述了一個功能不大的3D分析選項。59在Map View 和 Seismic View 之間使用光標追蹤你所看到河道的地方。60如果你還有時間，重復流程4.2的步驟，作出10，20，30和40Hz的Polygons 圖。10，20，30和40Hz的Polygons 圖能顯示河道隨頻率的變化嗎？這或許是一個暗示，暗示著河道厚度的不協(xié)調(diào)。調(diào)諧圖合成 (流程4.3) WF_1.0_M

44、erged_Frequencies 演示圖注意：為什么頻率合成在一起后，移除演示圖上一些河道信息?；氐?Tuning Mapper，不選Extract Slices as Horizons 并改變切片范圍從Range 0 到 80 。確認其他參數(shù)的匹配性（參考流程4.1）61保持其他的參數(shù)不變并點擊Compute Horizons from Slices。當完成調(diào)諧成圖工作后，刷新Map View ，演示修改后的 WF_1.0_Merged_Frequencies 圖。這個演示顯示地震切片到分立的頻率切片的強大的分解能力。通過分解后，你可以避免來自未調(diào)諧頻率成分的破壞性的干擾。62RGB超覽看

45、圖器 (流程4.4) 超覽允許你推算你的河道的RGB演示。在一張圖上，你可以把RGB演示當成3個不同的估值。 如果須要打開 PowerView，Tools PowerView 從Tuning Mapper上。 通過選擇 File Load Data 或者點擊Data Loader 釋放 Data Loader。如下所示：點擊層位欄，選擇WF_1.0_020Hz， WF_1.0_030Hz， WF_1.0_040Hz 層。點擊“ transfer arrow” 移動所選擇的層位到當前列表中。63從主窗口，選擇Tools Image Composer 彈開如下的對話框。64點擊RGB Blendi

46、ng 欄在下面的邊框中，再點Preview 欄。你的混合圖將顯示在PowerView 主窗口中。演示圖中20，30和40Hz分別對應紅，綠和藍。在圖上顯示為白色的地方，所有的頻率都有一個相似的較高的地質(zhì)響應；顯示為單色的地方（紅，綠和藍），僅其對應的頻率在區(qū)域內(nèi)已調(diào)諧。這樣的話， RGB混合的演示圖把效果調(diào)諧到單一的圖上。你可以疊加RGB演示，形成一個和分立頻率能量體相似的動畫。你可以通過更改最大最小值或者他們的滑動條，修改各種顏色演示的振幅范圍。當你對演示圖滿意的時候，輸入圖名并點擊Create Composite 按紐。65離散頻率體（Discrete Frequency Volume I

47、nterpretation ）66 離散數(shù)據(jù)體解釋、檢測，在目的層分析之后進行。調(diào)諧圖把問題限定在了一個局部的目的層段的范圍上。而表征大型地震數(shù)據(jù)體時，則需要另外不同的方法。對于超出單反射數(shù)據(jù)包或目的層段的分頻處理，我們推薦你使用“Discrete frequency energy cubes”（離散頻率能量調(diào)諧體）。這些能量的數(shù)據(jù)體，本質(zhì)來說是整個地震數(shù)據(jù)體的一個三維時窗分析，用連續(xù)的時窗。 離散頻率體是由一個單一的地震數(shù)據(jù)體計算到多元的離散頻率、振幅和相位的數(shù)據(jù)體。 通過一個滾動譜分析時窗,為每個模塊計算出其振幅或者相位。這些譜分量會被歸類成相同頻率成分的模塊。這種方法通常是再分析及參數(shù)化

48、目的層后進行。 67 本流程解決如下問題: 還有其他新的用來鉆井的井位嗎? 如何提高斷層的解釋精度? 流程將產(chǎn)生在單一頻率下的振幅信息的數(shù)據(jù)體。對體的解釋流程5.0 此流程與三維時窗分析一樣，不同的是一個頻率需要計算一個數(shù)據(jù)體。舉個例子，你可以圍繞產(chǎn)層（產(chǎn)油層）通過休整輸入的初始數(shù)據(jù)，進行分析，并運行“Volume Recon ”僅選30Hz。68儲層厚度估算（Thickness Estimation）69 頻譜分析提供一個獨立相位的穩(wěn)定的方法來估算地震上的儲層厚度。借助分頻處理，即使是一個相對低的頻率成分，如10Hz，都可以定量分析薄層的變化情況。 在油田上，任何定量分析的關(guān)鍵都在于有井控制

49、。正象地震反演技術(shù)一樣，分頻處理厚度分析需要井控來保證最終結(jié)果的合理性。 從塊狀模型薄層分析的理論來看，振幅譜的周期于頻率Pf 存在V型切口的關(guān)系,而頻率Pf與薄層的時間厚度t有關(guān),如下:t = (2 * z) / Vt = 1 / Pf 其中z等于薄層厚度, V等于速度.可求得時間厚度z:z = V / (2 * Pf) 直達波的頻率,f0等于1/2 Pf ,將這個關(guān)系式代入厚度等式,你會得到:z (厚度) = V / (4 * fo) t (雙程旅行時) = 1 / (2 * f0) 很多時候,基于頻譜分析的薄層厚度分析,可以確保薄層厚度與有效砂巖厚差在1米的范圍內(nèi).而當傳統(tǒng)的厚度分析產(chǎn)生歧義或者不能解決問題的時候,分頻處理的厚度分析就突顯了它的作用.顯然地,伴隨著反射系數(shù)模型復雜性和綜合性的不斷增加,它必將可以完善振幅譜的解釋,并結(jié)合其相互的聯(lián)系.這項研究正在進行中. 本流程解決如下問題: 儲層的相對厚度是多少?70儲層相對厚度圖 (流程6.1) SpecDecomp的(R2003.0)版, 主要集中于定性分析.分析結(jié)果將顯示地層是厚是薄,期望的有利井區(qū)等.但它不能使你獲得精確的儲集巖厚度. 再次運