Geolog-全波列聲波測(cè)井中文手冊(cè)-

1、Geolog 軟件技術(shù)手冊(cè)Full Sonic Wave Processing -SWB帕拉代姆公司北京代表處2006 年 12 月1 、綜述 - 1 -1.1 預(yù)備知識(shí) -. .1. .-1.2 數(shù)據(jù) -. .1. .-2、陣列聲波全波形 - 2 -2.1 數(shù)據(jù)準(zhǔn)備 -. .3. .-2.1.1 查看 /創(chuàng)建一個(gè)聲波列陣工具模版 -. .3. .- 2.1.2 練習(xí)指導(dǎo) 2- 創(chuàng)建其他波形屬性 -. .5. .- 2.1.3 波形分解 -. .6. .-2.1.4 深度轉(zhuǎn)換 -. .7. .-2.2 處理 -. .8. .-2.2.1 數(shù)據(jù)分析 -. .8. .-2.2.2 去噪 - 11

2、-2.2.3 設(shè)計(jì)濾波器 - 17 -2.2.4 振幅恢復(fù) - 19 -2.3 陣列聲波處理 - 20 -2.3.1 處理模塊簡(jiǎn)介 - 20 -2.3.2 偶極波形處理 - 21 -2.3.3 單極波形處理 - 23 -2.3.4 拾取標(biāo)志波至 - 26 -2.4 后期處理 322.4.1 綜述 322.4.2 頻散校正 332.4.3 傳播時(shí)間疊加 362.4.4 相關(guān)性顯示 382.4.5 陣列聲波重處理 393 、機(jī)械性質(zhì) 443.1 綜述 443.2 計(jì)算動(dòng)力學(xué)彈性性質(zhì) 44附錄 I- 快速運(yùn)行 46附錄 II- 頻散校正討論 471、綜述歡迎閱讀 Geolog軟件SWB指導(dǎo)教程。這是

3、帕拉代姆公司 Geolog軟件 SWB模塊的使用教程，通過典型工作流程和程序模塊的使 用，使用戶熟悉陣列聲波處理過程的基本工具的使用。熟悉陣列聲波波形 查看/創(chuàng)建一個(gè)陣列聲波工具模版 創(chuàng)建一個(gè)給定波形的其他屬性 解裝測(cè)井記錄 預(yù)處理及波形數(shù)據(jù)分析 使用平均深度、平均時(shí)間和頻率濾波器去噪 去除數(shù)據(jù)采集中振幅增益和標(biāo)準(zhǔn)化 處理偶極和單極聲波測(cè)井儀下波形 自動(dòng)和交互式的拾取波至頻散修正傳播時(shí)間覆蓋查看結(jié)果的相關(guān)性圖形1.1 預(yù)備知識(shí)Geolog 基本測(cè)井記錄中文檔，所有的輸入文件為 Bold Courier New ，輸出為 Bold Courier New，不加粗。1.2 數(shù)據(jù)教程中使用了下面的附

4、加文件（文件并非由軟件提供）數(shù)據(jù)： geophysics_master.unl曲線模板： (copy from layouts ) swp_array_sonic_comparison swp_dispersion swp_depth_average swp_frequency03 swp_frequency05 swp_other_attributes swp_projection swp_raw_sonic_waveforms swp_reprocessing swp semblance swp traveltime函數(shù)：圖像：井：polaris 03, polaris 042、陣列聲波全

5、波形典型的陣列聲波工具在接收剖面上有 8個(gè)接收器和 3個(gè)不同的發(fā)射器。也要注意的是，不 同服務(wù)公司出產(chǎn)的陣列工具也不相同。在數(shù)據(jù)采集上，基于不同的調(diào)查研究目的采用不同的工作模式。在Geolog 里，“WF”是慣用的波形記錄名字后跟兩位數(shù)字，首位數(shù)代表操作模式索引，第二位數(shù)表示接收器索引。例 如， WF21代表波形記錄為模式 2的第一記錄道（基于偶極模式）。每一個(gè)接收器的聲波波形以數(shù)組記錄存載。 用Geolog layout（ 模版）可以在圖像道和陣列道 顯示一個(gè)聲波波形。在這兩個(gè)道里，水平軸表示信號(hào)的旅行時(shí)，垂直軸為參考深度。當(dāng)選定 一個(gè)圖像道后，用戶可以打開一個(gè)聲波陣列查看窗口，對(duì)數(shù)據(jù)做進(jìn)一

6、步分析。圖1：原始聲波波形顯示圖1為原始聲波波形顯示圖像道（第一道和第三道），陣列道（第二道和最后一道）。最 初的兩道顯示的是第一接收器的偶極波形記錄，最后兩道是第一接收器的單極波形記錄2.1 數(shù)據(jù)準(zhǔn)備通過本節(jié)可以熟悉聲波全波形數(shù)據(jù)準(zhǔn)備 查看/創(chuàng)建一個(gè)列陣聲波工具模版 創(chuàng)建給定波形的其他屬性 手動(dòng)解裝波形 進(jìn)行深度轉(zhuǎn)換2.1.1 查看/ 創(chuàng)建一個(gè)聲波列陣工具模版在波形處理前先了解工具的相關(guān)信息， 它是以文件夾形式存儲(chǔ)。 斯倫貝謝公司的 DSI測(cè)得 的井資料在 Geolog中已經(jīng)存儲(chǔ)在 Site目錄下，作為四個(gè)缺省文件（ dsi_m1 to dsi_m4）。當(dāng)有 改動(dòng)或者初次運(yùn)行此工具，相關(guān)信息

7、會(huì)保存在本工具文件夾下的 specs路徑中。工具模版的數(shù)據(jù)信息可以在 Text中的 Constants 表格中看到1、啟動(dòng) Geolog軟件打開 STARS工區(qū)2、點(diǎn)開 Well 3、點(diǎn)擊 well>open ，打開 Polaris_03 4 、選擇 Petrophysics > Full Waveform Sonic > Data Preparation > Create Tool Specification模塊顯示斯倫貝謝公司的 spec文件dsi_m2 ，給定的為缺省值，如下圖 2圖 2： array_sonic_create_toolspec中 dsi_m2 的

8、缺省值Spec文件中的參數(shù)如下表所示：參數(shù)描述TOOL NAME聲波采集儀器的名稱如： DSI, XMACTOOL SPEC聲波工具模版，以模版文件命名TOOL MODE聲波工具的工作模式，通常在 Text> Constants 表格中定義TOOL R聲波儀器的直徑SAMPLE RATE聲波采樣率，通常以微次生法測(cè)量TR DISTANCE發(fā)射器與第一接收器的距離X INTERV兩個(gè)鄰近接收器之間的距離DEPTH_SHIFT工具深度的平差不同于工具的參考深度和井的參考深度。負(fù)的值表示工具參考深度在井 的參考深度之下5、如果需要的話可以對(duì)缺省值做調(diào)整6、點(diǎn)擊 Start 運(yùn)行模塊 一個(gè)新的模

9、塊文件在本工區(qū)的 SPECS目錄下產(chǎn)生了7、選擇 Launcher > Close 來關(guān)閉陣列聲波工具模塊8 、檢查新模塊 dsi_m2 的值 ，如圖 3 所示圖3： dsi_m2 的存儲(chǔ)位置以及值提示：不管Geolog 中使用的量度制，米制的或是特定的，聲波陣列工具模版總是會(huì)保存為米制的。當(dāng)要 運(yùn)行這個(gè)工具，首先請(qǐng)確認(rèn)工具模版的值，以確保得到正確的結(jié)果 。2.1.2 練習(xí)指導(dǎo) 2- 創(chuàng)建其他波形屬性聲波波形可以用反射振幅正常顯示。頻率譜、瞬時(shí)相位、和瞬時(shí)振幅屬性可以從解釋角 度很好的量化突出這些特性。原始聲波頻率譜顯示的是不同到達(dá)時(shí)的主要頻率圍，這樣可以 在頻率域查看各到達(dá)模式的屬性特

10、征。瞬時(shí)相使弱信號(hào)加強(qiáng)，但是也同樣使噪音增強(qiáng)。瞬時(shí) 振幅是包絡(luò)的振幅，包含更多側(cè)向變化，但是降低了垂向分辨率。 tp_array_sonic_attribute 模 塊用于計(jì)算一個(gè)給定波形的這些屬性。 三個(gè)測(cè)井記錄可以從這個(gè)模塊輸出 - 頻率譜、瞬時(shí)相位 以及瞬時(shí)振幅，可以在圖像道或者陣列道中查看這些屬性。1、選擇 Petrophysics > Full Waveform Sonic > Data Preparation > CreateOther Attributes 來顯示 tp_array_sonic_attribute 模塊。2、如圖 4所示，對(duì)窗口做如下值改變Inp

11、ut Set改為SONICOutput Set改為SONIC_ADJLOG_IN類型為 WF21FREQ_OUT類型為 WF21FREQINST_PHASE類型為 WF21PHASEINST_AMP類型為 WF21AMP圖 4: 改變屬性模型參數(shù)3、點(diǎn)擊 Start 運(yùn)行4、打開曲線模板 swp_other_attributes.layout 來查看結(jié)果 ，如圖 5所示。原始數(shù)據(jù) 原始數(shù)據(jù)頻率譜 原始數(shù)據(jù)瞬時(shí)相位 原始數(shù)據(jù)瞬時(shí)振幅 圖 5 ：波形屬性顯示2.1.3 波形分解來自每個(gè)接收器的聲波波形有時(shí)被疊加記錄到一起。這個(gè)模塊用于把波形分解為每個(gè)接 收器的單個(gè)記錄。通常，每個(gè)模塊有自有的測(cè)深記

12、錄，這個(gè)記錄要作為輸出的參考深度。在 這個(gè)模塊里，需要提供分解波形的名字前綴。要遵守 Geolog 的命名約定。為了做示，本例中 的未分解波形都刻意不進(jìn)行自動(dòng)分解。1、選擇 Petrophysics > Full Waveform Sonic > Data Preparation > LogUnpacking 來顯示 tp_array_sonic_unpack 模塊2、如圖 4所示，對(duì)窗口做如下值改變 Input Set改為SONIC_PACKED PREFIX類型為WG DEPTH選擇 PWD2 LOG_WF選擇 PWG2圖 6: 改變分解波形模型參數(shù)3、點(diǎn)擊 Start

13、運(yùn)行4、將PREFIX值變?yōu)?WN ，LOG_WF值變?yōu)?PWN2 ，再運(yùn)行模型5、打開 Text ，檢查未分解記錄值。如下圖 7所示圖 7 ：分解波形輸出記錄2.1.4 深度轉(zhuǎn)換DSI工具數(shù)據(jù)的參考深度與井記錄的參考深度不同， 差值從測(cè)量工具的底部到接收器剖面 的中段進(jìn)行修正，如上節(jié)的 depth-shift 。1、選擇曲線道模版 swp_other_attributes2、選擇 Tools > Depth Shift > Apply Logs3、如圖 8所示，對(duì)窗口值做如下改變 Input Set 改為 SONICOutput Set改為SONIC_DEPTH_SHIFTDEP

14、TH_OFFSET輸入?yún)?shù)為 -7.9248LOG_IN選擇波形 WF21_1為WF28_1 WF41_1為WF48_1CORRECTED選擇DEPTH圖 8： 深度轉(zhuǎn)換模型4、點(diǎn)擊 Start 運(yùn)行5、關(guān)閉所有窗口2.2 處理通過本節(jié)可以使用戶熟悉全波形聲波處理流程 在聲波數(shù)組查看窗口操作和解釋數(shù)據(jù) 使用平均時(shí)間和平均深度以及 Butterworth 頻率濾波器去除波形噪音 設(shè)計(jì)一個(gè)濾波器去除數(shù)據(jù)采集中的振幅增益和標(biāo)準(zhǔn)化2.2.1 數(shù)據(jù)分析 為了保證處理結(jié)果質(zhì)量，在數(shù)據(jù)處理前要檢查數(shù)據(jù)確保無誤。波形數(shù)據(jù)通常包含噪聲。通 常需要選擇適當(dāng)?shù)墓ぞ邅肀M可能多的去噪。1、打開模板 swp_raw_so

15、nic_waveforms.layout 。2、選擇 WF21 的圖像道，選擇 View > Array Sonic View 。 數(shù)組聲波查看窗口顯示波形 WF21 以及它當(dāng)前深度頻率域的位置。（見圖 9）要打開陣列聲波查看窗口，需要選擇打開帶數(shù)據(jù)的圖像道陣列聲波查看窗口圖 9 ：圖像道和聲波數(shù)組窗口陣列聲波窗口用來詳細(xì)的分析波形，它有兩個(gè)顯示區(qū)域。顯示區(qū)在當(dāng)前深度位置顯示信 號(hào)波形，可以估算出時(shí)間間隔（接收器的時(shí)間延時(shí)）。頻率域顯示區(qū)顯示每個(gè)波形相應(yīng)的的 頻率圍。3、點(diǎn)擊聲波數(shù)組記錄 - Log Select 圖標(biāo)打開可選擇多條曲線。 提示：挑選出來的記錄必須來自同一個(gè)集合和模型工具

16、（例如：有相同的模型索引- 字符名字后的第一位數(shù)相同）在圖像道上顯示。4、選擇記錄 SONIC.WF21_1 to SONIC.WF28_1 ，點(diǎn)擊 OK。 多重曲線以及其對(duì)應(yīng)的頻率域在數(shù)組聲波的窗口顯示。5、在模版上第一圖像道 （SONIC.WF21） 點(diǎn)擊任何位置。聲波數(shù)組顯示窗口的信號(hào)在每次 選擇深度點(diǎn)的時(shí)候會(huì)刷新。另一個(gè)方式是通過深度區(qū)域，鍵入圍值或者使用上下按鍵。6、設(shè)置深度點(diǎn)和時(shí)間圍如下所示：在顯示區(qū)放大或者改變時(shí)間圍7、將光標(biāo)定位在信號(hào)顯示區(qū)任意位置， 按住鼠標(biāo)左鍵拖曳定義一個(gè)矩形區(qū)域，放開鼠標(biāo) 左鍵。通過顯示的時(shí)間圍值來反映被選擇的區(qū)域。將時(shí)間圍值的起始和終點(diǎn)值輸入與用鼠標(biāo)改變

17、信號(hào)顯示區(qū)域的效果是 一樣的。8、再次點(diǎn)擊鼠標(biāo)左鍵使區(qū)域變?yōu)樵瓉泶笮　?在波形上估算時(shí)間間隔： 提示： 這個(gè)功能只在多重波形在當(dāng)前活動(dòng)的時(shí)候起作用。當(dāng)時(shí)間間隔固定時(shí)，功能不可選時(shí)首先將信號(hào)顯示區(qū)放大9、固定時(shí)間間隔，在慢度剖面顯示。10、顯示區(qū)，在第一個(gè)接收器 （WF21）的信號(hào)，按住鼠標(biāo)左鍵，拖至到最后一個(gè)接收器 的信號(hào)處 （WF28 ），放開鼠標(biāo)左鍵。如圖 10淡紅色顯示 估算時(shí)間間 隔，波至慢度 及接收間隔粗淡紅色的 陡線顯示第 一 個(gè)和 最后 一個(gè)接收器 初至?xí)r間間圖 10 ：時(shí)間間隔估算 8個(gè)接收器的信號(hào)以及頻率域在當(dāng)前選擇的深度顯示。11、不選 Moveout 復(fù)選欄，不顯示時(shí)間間

18、隔數(shù)據(jù)。12、關(guān)閉陣列聲波顯示窗口。2.2.2 去噪 波形圖總是會(huì)包含噪聲，太多的噪聲會(huì)不能識(shí)別可檢測(cè)的波至。預(yù)處理是提高數(shù)據(jù)質(zhì)量 保證初至被識(shí)別的重要步驟。Geolog 有一系列基本的預(yù)處理工具可以滿足一個(gè)波形的任意去噪要求。 去除時(shí)間均值，在時(shí)間方向使用平均時(shí)濾波器 去除深度均值，在深度方向使用平均深度濾波器Butterworth 濾波器，是一個(gè)頻率濾波器 根據(jù)具體的需求，用戶可決定針對(duì)某個(gè)波形使用對(duì)應(yīng)的濾波器 （1）去除時(shí)間均值噪聲 這個(gè)模塊提供沿時(shí)間方向去除噪聲濾波器。先在每個(gè)窗口計(jì)算出平均值，然后從原始波形移除噪聲。) 這是用于移除低頻噪聲背景的濾波器 。1 、選擇 Petrophy

19、sics > Full Waveform Sonic > Pre-Processing > RemoveTime-Average2、如圖 10，對(duì)打開窗口中值做如下修改。改為SONIC_DEPTH_SHIFT Input Set UTIME WINDOW LOG_IN例如，設(shè)置濾波窗口值為 520 in micro-seconds (usecs)忽略采樣場(chǎng)。這個(gè)模塊以取得合適的結(jié)果。設(shè)為WF21 (可以同時(shí)導(dǎo)入多個(gè)文件)圖11 ：時(shí)間平均濾波器參數(shù)設(shè)置是在Geolog不能識(shí)別采樣率時(shí)需要采樣場(chǎng))改變窗口的取值，多次運(yùn)行3、點(diǎn)擊 Start 運(yùn)行模塊4、打開繪圖模版 swp_a

20、rray_sonic_comparison.layout 。圖13 ：設(shè)置平均深度濾波器模型值圖12：用時(shí)間平均濾波器濾波前后對(duì)比圖（2）去除深度均值噪聲 這個(gè)模塊提供沿深度方向去除噪聲濾波器。先在每個(gè)活動(dòng)窗口計(jì)算出平均值，然后從原 始波形移除噪聲。這個(gè)濾波器可去除典型的諸如套管信號(hào)或者儀器帶來的噪聲。1、打開井 Polaris_04 ，保存對(duì)先前井所作的改變。 2、關(guān)閉所有打開窗口。3、選擇 Petrophysics>Full Waveform Sonic>Pre-Processing> RemoveDepth-Average.4、在 tp_array_sonic_dept

21、h_average 窗口對(duì)值做如下改變（如圖 13 ） Input Set改為SONIC WINDOW設(shè)置濾波器長度為 21 LOG_INWF415、點(diǎn)擊 Start 運(yùn)行模型6、打開繪圖模版 swp_depth_average.layout 。套管和儀器形成的信號(hào)可以很強(qiáng)，使得重要的波至難以識(shí)別。圖 14 左邊道顯示的波形是一個(gè)這樣的例子，直直的垂線占主要部分。右邊的道顯示的是不需要信息被移除后的清晰的 波至。運(yùn)行濾波器后的結(jié)果圖圖14：用深度濾波器濾波前后對(duì)比圖3）使用頻率濾波器Butterworth 濾波器是一個(gè)頻率濾波器， 可以通過給定一個(gè)高的和低的截止值以及創(chuàng)建濾 過順序來實(shí)現(xiàn)。 B

22、utterworth 濾波器選擇的順序 （從 1到8）決定了從通頻帶到抑制帶的過渡陡度。單位設(shè)為千赫茲而不是赫茲，是為了便于進(jìn)行聲波波形的處理。如果數(shù)據(jù)是來自于標(biāo)準(zhǔn) 的DSI工具采集，可以忽略采樣率，因?yàn)?Geolog 可以自動(dòng)地讀取。有關(guān)濾波器的調(diào)整適應(yīng)更多信息參考菜單選項(xiàng)下的 Design Filter 。1、保存對(duì) Polaris_04的改動(dòng)，打開井 Polaris_03 。2、選擇 Petrophysics > Full Waveform Sonic > Pre-Processing > Frequency Filter3、對(duì)打開窗口做如下改變（如圖 15所示） In

23、put SetSONIC_DEPTH_SHIFT Output Set HIGH_END改為FILTERED03 SAMPLE_RATE 10 LOG_IN從輸入集合選擇8個(gè)波形： SONIC_DEPTH_SHIFT.WF21_1 to設(shè)為 3HZ（移除所有高于這個(gè)值之外的頻率信號(hào) ）WF28_1圖 15 ：改變 butterworth 濾波器的頻率圍和大小4、點(diǎn)擊 Start 運(yùn)行這個(gè)模型5、打開 swp_frequency03.layout6、選擇 View > Array Sonic View 來查看濾波后圖形7、設(shè)置當(dāng)前深度到如圖 16所示。在 3HZ以上的頻率已經(jīng)在當(dāng)前深度位置

24、被移除圖16：WF濾波后圖形8、關(guān)閉陣列聲波查看窗口。9、選擇 tp_array_sonic_filter 模塊，重復(fù)第二步和第三步改變頻率圍為 0-5 赫茲，輸出名 為FILTERED05。10、點(diǎn)擊 Start 運(yùn)行模塊11、打開 swp_frequency05.layout12、打開聲波數(shù)組查看窗口。大于5Hz 的頻率信號(hào)都被移除，如圖 17所示，在當(dāng)前深度位置顯示的頻率譜。圖17：濾除高于 5Hz波后的 WF21 波形圖13、關(guān)閉陣列聲波窗口2.2.3 設(shè)計(jì)濾波器Butterworth 濾波器將作為頻率濾波器用在下面的例子中。 Butterworth 是一個(gè)很好的帶通 濾波器， 在損失

25、了部分頻率域后有較為合理的時(shí)間衰減。 高頻和低頻截止值用于帶通和帶止。 對(duì)于低通濾波器， 高頻波被截止， 反之亦然。 Butterworth 濾波器的排列控制著在通頻帶之外 信號(hào)變?nèi)醯目炻?。設(shè)計(jì)一個(gè)合適的濾波器或是子波，必須要注意采樣率和頻率值。設(shè)計(jì)一個(gè)Butterworth 子波，采樣率( SR)要使用 0.001秒，這樣可以得到 Nyquist frequency = 1 / (2 * SR) = 500 Hz 。 輸入頻率必須要小于尼奎斯特頻率( Nyquist frequency )頻率。也就是說，如果提供的值大 于尼奎斯特頻率( Nyquist frequency )，采樣率也要更好

26、。1、選擇 Petrophysics > Full Waveform Sonic > Pre-Processing > DesignFilter2、在 filter_create 窗口對(duì)值做如下改變（如圖 18 ）FILTER_TYPE改為ButterworthFILTER_NAMEbw1040LOW_START1000HIGH_END3000PLOT_SR0.00004DISP_SCALEAMPLITUDE圖18：濾波器創(chuàng)建模塊3、點(diǎn)擊 Start 運(yùn)行這模塊。4、查看創(chuàng)建的濾波器（如圖 19）。 Input Set NORM_OPT GAIN_OPT REFERENCE選

27、擇UNPACK.DEPTH圖19 ：運(yùn)行創(chuàng)建濾波器模塊后的結(jié)果2.2.4 振幅恢復(fù)在數(shù)據(jù)采集中，振幅增益以及標(biāo)準(zhǔn)化是為了增強(qiáng)波形信號(hào)的質(zhì)量。可基于有效的增益和標(biāo)準(zhǔn)化因素（歸一化因素），使波形回到原來形狀。默認(rèn)的，去除標(biāo)準(zhǔn)化因素選擇inverse ；去除增益選擇為 multiply 。1、選擇Petrophysics > Full Waveform Sonic > Pre-Processing > Amplitude Recovery.2、在 tp_array_sonic_tar 窗口對(duì)值做改變（如圖 20所示）改為SONIC_DEPTH_SHIFT設(shè)為INVERSE設(shè)為MUL

28、TIPLY LOG_WF選擇 WF21_1至WF28_1波形 LOG_NORM 選擇井 UNPACK.WN21_1 到 UNPACK.WN28_1 LOG_GAIN選擇井 UNPACK.WG21_1 到UNPACK.WG28_1圖20 ：振幅恢復(fù)模塊3、點(diǎn)擊 Start 運(yùn)行模型。4、打開text 窗口在原來的集合中查看井?dāng)?shù)據(jù)(如圖 21所示)。圖21 ：Text 窗口查看去除振幅增益后的波形記錄2.3 陣列聲波處理通過本節(jié)，用戶可以：處理偶極波形處理單極波形使用自動(dòng)交互式拾取波至2.3.1 處理模塊簡(jiǎn)介聲波波形處理基于慢度時(shí)間相關(guān)法。這個(gè)模塊在活動(dòng)窗口過時(shí)間和慢度方向可以查看多個(gè)聲波。默認(rèn)的

29、，時(shí)間圍可以包括所有的記錄時(shí)間；慢度的圍至少包括目的區(qū)。例如，慢度圍為100800 us/m (1200-10000 m/s) 可以包含大多數(shù)區(qū)域的縱波波至。建議用較大的慢度圍這樣可以盡可能包含所有波至，但這也帶來更多的計(jì)算量首先選擇少量數(shù)據(jù)來決定時(shí)間和慢度方向的最優(yōu)選擇要注意慢度步長的選擇，過大的步長會(huì)降低最終聲波分辨率。在每個(gè)深度點(diǎn)計(jì)算出的最大值作為與深度相關(guān)的函數(shù)映射到慢度軸上；這是最終輸出Projection log 。這個(gè)模塊有很多控制參數(shù)，詳細(xì)參數(shù)解釋可以查看模塊幫助。2.3.2 偶極波形處理在本節(jié)，被處理的波形filtered03 和 filtered05 將用于下文的離散。1、

30、保存井文件 Polaris_03。2、關(guān)閉所有查看窗口。3、選擇 Petrophysics> Full Waveform Sonic > Processing > Waveform Processing來打開 tp_array_sonic_process 窗口。4、設(shè)置輸入?yún)?shù)（如圖 22所示） Input SetSONIC_DEPTH_SHIFT Output SetDIPOLE TOOL_SPECdsi_m2 START0 STOP10000SLOWNESS_FROM 35 SLOWNESS_TO500 SLOWNESS_STEP4 LOG_IN由SONIC_DEPTH_

31、SHIFT集合，選擇8個(gè)波形： WF21_1 to WF28_1圖22：偶極波形處理模型窗口5、點(diǎn)擊 Start 運(yùn)行模型。當(dāng)模型運(yùn)行結(jié)束，結(jié)果以相關(guān)性圖件顯示，如圖 23所示。數(shù)據(jù)在當(dāng)前深度位置顯示（當(dāng)深度未規(guī)定時(shí)是在井記錄的頂部）圖23 ：相關(guān)性圖 - WF21 - WF28 波形處理結(jié)果圖6、保存相關(guān)性圖為 semblance_plot_01 。將相關(guān)性圖以唯一文件名保存，便于以后查看，否則圖件會(huì)被覆蓋7 、改變輸入和輸出集合的名字為FILTERED03，重復(fù) 4-6 步，將相關(guān)性圖形保存為semblance_plot_03 。8、改變輸入和輸出集合的名字為FILTERED05，重復(fù)4-

32、5 步，將相關(guān)性圖形為semblance_plot_05 。9、關(guān)閉相關(guān)性圖。10、保存井文件。2.3.3 單極波形處理11、在 tp_array_sonic_process窗口，設(shè)置輸入?yún)?shù)，如圖 24所示： Input SetSONIC_DEPTH_SHIFT Output SetMONOPOLE TOOL_SPECdsi_m4 START0 STOP5120 UTIME_WINDOW400 UTIME_STEP50 SLOWNESS_FROM30 SLOWNESS_TO298.7 SLOWNESS_STEP1.98 LOG_IN從SONIC_DEPTH_SHIFT集合，選擇 8 個(gè)波形

33、WF41_1to WF48_1圖24：?jiǎn)螛O波形處理模型窗口2D 相關(guān)12、點(diǎn)擊 Start 運(yùn)行模塊。當(dāng)模型運(yùn)行結(jié)束輸出一個(gè)映射和一個(gè)相關(guān)性圖件，圖件在每個(gè)深度點(diǎn)存儲(chǔ)的是 性圖像， 可以通過陣列聲波數(shù)組查看窗口進(jìn)行質(zhì)量控制。如下圖所示圖25 ：陣列聲波數(shù)組窗口查看圖件13、關(guān)閉相關(guān)性圖件，不保存改變。14、打開模版 swp_projection.layout 。上面例子中的偶極和單極子映射記錄如圖 26所示偶極波形處理后圖形 單極子波形處理后圖形 水平軸顯示慢度 (START)到 (END) 圖26：深度相關(guān)函數(shù)映射15、在WELL下點(diǎn)擊 SAVE.2.3.4 拾取標(biāo)志波至映射記錄在慢度 -

34、深度平面代表了一個(gè)或多個(gè)相關(guān)性。它們與某些波至有較好的一致性。 模型在輸入的映射記錄中自動(dòng)或者交互式的尋找相關(guān)最大值位置，從慢度投影記錄拾取相應(yīng) 一組波至。搜索從一個(gè)起始值開始(初始慢度值在映射記錄里)來標(biāo)定某一初至，搜尋圍以 從映射記錄估計(jì)的起始值為中心。估算輸入起始值最簡(jiǎn)便方法是將投影記錄插入到圖像道， 打開陣列生波可視窗口， 移動(dòng)鼠標(biāo)標(biāo)定出信號(hào)顯示區(qū)的目的位置， 光標(biāo)位置 X的值會(huì)在位置對(duì) 話框的狀態(tài)欄顯示出來。這個(gè)值就是從圖上讀取的慢度值，它可以用來作為搜索起始值。對(duì)于交互式拾取，測(cè)井曲線首先是自動(dòng)被拾取(參見 Curve Insert 聯(lián)機(jī)幫助)，編輯，最 后作為備用初始記錄被保存。

35、主要有：1、自動(dòng)拾取一個(gè)偶極測(cè)井記錄2、自動(dòng)拾取一個(gè)單極測(cè)井記錄3、對(duì)于第二步中不合適點(diǎn)做交互式拾取自動(dòng)拾取一個(gè)偶極測(cè)井記錄1 、選擇 Petrophysics > Full Waveform Sonic > Processing > Picking 來打開 the tp_array_sonic_label 窗口。2、選取輸入起始值 (SLOWNESS_S1)。讓swp_projection 為活動(dòng)窗口，選擇偶極投影道來打開陣列聲波可視窗口。在信號(hào)顯示區(qū)移動(dòng)鼠標(biāo)到拾取位置， 讀取顯示出的 X值。這個(gè)值即是從圖上讀取的慢度 值，可作為起始值 (SLOWNESS_S1)，如例子中

36、圖 27示。關(guān)閉陣列聲波可視窗口。圖27 ：確定搜索起始值3、選擇模型窗口，設(shè)置值如下圖 28所示 Input/Output Set選擇 DIPOLE，這是處理后的輸出集合圖29：自動(dòng)拾取的起始記錄覆蓋在偶極映射圖像上 SLOWNESS_WINDOW 本例中設(shè)為 121 us/f （因?yàn)閮H有一個(gè)波至，對(duì)于 偶極使用一個(gè)大點(diǎn)的圍是可靠的）窗口長度必須包括波至的最大波動(dòng)幅 度，但是也不能過大要注意排除鄰區(qū)波至。 SLOWNESS_S1 在前一步操作中， 位置對(duì)話框里確定的 X 值，必須要接近值121圖28偶極自動(dòng)拾取模型設(shè)置4、運(yùn)行模型 5、重復(fù)地3-4 步二次，改變輸入輸出集合名字分別為FILT

37、ERED03，F(xiàn)ILTERED05。6、設(shè)定 swp_projection 為活動(dòng)窗口，查看曲線道（如圖 29示） PICK_OPTIONBOTH7、打開屬性對(duì)話框，注意 DIPOLE.DT_PICK1井的左右邊界為處理的慢度起 始和終止值(35 and 500 US/F )，輸入模型中。8、關(guān)閉對(duì)話框。單極投影記錄自動(dòng)拾取9、確定輸入 tp_array_sonic_label 模型所需要的的起始值 (SLOWNESS_S1 and SLOWNESS_S2)選擇單極圖像道打開陣列聲波可視窗口。在信號(hào)顯示區(qū)使用鼠標(biāo)定位， 在定位對(duì)話框確定兩個(gè)初始值， 如圖 30所示 為確定第二波至的例子。圖30

38、 ：確定 SLOWNESS_S2的輸入初始值10、使 tp_array_sonic_label 為活動(dòng)窗口，對(duì)模型值做如下設(shè)置(如圖 31所示)： Input / Output Set選擇MONOPOLE集合，它是處理中的輸出集合 Slowness_Window設(shè)為15.24 us/f 為搜索窗口大小 SLOWNESS_FROM 30 SLOWNESS_TO 298.7 SLOWNESS_STEP 1.98 SLOWNESS_S1 60 SLOWNESS_S2 85圖31 ：交互式拾取兩個(gè)井記錄模型界面11、運(yùn)行模型 12、選擇模板 swp_projection 作為查看模板，如圖 32所示。

39、圖32：自動(dòng)拾取井記錄覆蓋在單極圖像道上13、菜單欄 WELL下點(diǎn)擊 SAVEInteractive Picking 交互式拾取14 、 用 swp_projection 模 板 ， 在 單 極 投 影 道 圖 像 道 上 選 擇 拾 取 記 錄DT_PICK2。15、選擇 Tools > Curve Insert 。16、手動(dòng)的重新拾取不合適的橫波波至。17、打開文本查看重命名被編輯的井記錄 DT_PICK2 version 2 為DT_SEED。18、菜單欄 WELL下點(diǎn)擊 SAVE。19、選擇 tp_array_sonic_label 模型，對(duì)值做如下改變（如圖 33）： PICK

40、 OPTION SINGLE SLOWNESS_S1使用DT_SEED作為初始記錄 DT_PICK1使用DTS_PICK圖33：?jiǎn)尉换ナ绞叭∧Ｐ徒缑?0、運(yùn)行模塊21、在曲線模板中查看結(jié)果，使用圖 34作為參考。 如果效果不夠好， 調(diào)整搜 索圍（例如： SLOWNESS_WINDOW 和SLOWNESS_S1值），重復(fù)先前的步驟。圖34：偶極和單極的投影記錄的自動(dòng)式和交互式拾取顯示在單極投影區(qū)，綠色的曲線表示作為初始記錄的自動(dòng)拾取曲線 (DT_PICK2) 藍(lán)色曲線是基于初始子記錄的交互式拾取 (DTS PICK)。提示：處理的慢度記錄需要進(jìn)行深度轉(zhuǎn)換，這涉及到井的參考深度。在Geolog

41、中不會(huì)自動(dòng)轉(zhuǎn)換，因?yàn)閷?duì)于參考深度工具，沒有足夠的信息可供選擇。2.4 后期處理通過本節(jié)，用戶將熟悉后期處理流程 計(jì)算并應(yīng)用頻散校正，確保每個(gè)井記錄在同一參考系下 應(yīng)用頻散校正導(dǎo)出旅行時(shí)查看相關(guān)性圖形顯示 陣列聲波數(shù)組重處理，計(jì)算平均慢度2.4.1 綜述后期處理以及后期質(zhì)量控制是聲波處理的重要部分。 從處理過的記錄中導(dǎo)出 的旅行時(shí)可以作為質(zhì)量控制來描述聲波波形波至。 頻散以及井眼情況等將影響最終處理結(jié)果頻散校正模型計(jì)算出頻散曲線，在每個(gè)深度位置使用校正值來修正頻散影 響。波形將按照一定間隔被重處理 （在后期處理中） 來增強(qiáng)垂直分辨率或者補(bǔ)償 井眼較大的直徑差值。2.4.2 頻散校正這個(gè)模塊輸出修

42、正后的慢度記錄。 它可以通過覆蓋陣列道隨意輸出頻散和頻 譜記錄，也可以輸出壞道來標(biāo)記有問題區(qū)域。1、菜單欄 WELL下點(diǎn)擊 SAVE。2、選 擇Petrophysics > Full Waveform Sonic > Post-Processing >Dispersion Correction 。提示：離散修正是一個(gè)很費(fèi)時(shí)的處理過程，本教程中為了實(shí)現(xiàn)這個(gè)操作只是用了少量 資料來處理。3、如圖 35所示，設(shè)置處理參數(shù) 。圖35：設(shè)置處理圍4、輸入?yún)?shù)如下：如圖 36所示 Input SetSONIC_DEPTH_SHIFT Output SetRAWTOOL_SPECdsi_m

43、2START0STOP10000CALIWIRE.HCALBSWIRE.BSRHOB_F1.24 G/C3(1240 K/M3)RHOB_SWIRE.RHOZDT_F245.8 US/F (75 US/M )DT_PMONOPOLE.DT_PICK1DT_SDIPOLE.DT_PICK1LOG_IN從 SONIC_DEPTH_SHIFT集合, 選擇8個(gè)波形WF21_1 to WF28_ 1圖36：頻散校正窗口5、運(yùn)行模塊處理速度取決于計(jì)算機(jī)的性能和數(shù)據(jù)量的大小6、重復(fù)第三和第四步，改變值 Input /Output setFILTERED03 DT_S7、重復(fù)第三和第四步，改變值FILTERE

44、D03.DT_PICK1 Input /Output setFILTERED05FILTERED05.DT_PICK1 DT_S8、打開模板 swp_dispersion.layout 。提示： 離散修正結(jié)果如圖 37 所示，處理圍為 13150 - 13450 ft.圖37 ：離散結(jié)果對(duì)比實(shí)例（數(shù)據(jù)為 13150 - 13450 ft ）1,2,3道分別表示覆蓋在疊加頻率域的波形離散曲線， 4道表示低頻濾波后的 波形， 5道波形大于5 Khz的頻率被濾過， 6道大于3 Khz頻率被濾過。對(duì)原始波形做處理時(shí)， 發(fā)現(xiàn)在做頻散校正后， 有不穩(wěn)定的波動(dòng)現(xiàn)象。 如上例， 圍如箭頭所示，通過0-5 KH

45、z帶通濾波處理，修正的剪切波慢度穩(wěn)定， 同時(shí)1-3 KHz 濾波器給出的離散都很清晰。9、菜單欄 WELL下點(diǎn)擊 SAVE。2.4.3 傳播時(shí)間疊加在前面步驟中得到的慢度記錄，可以從中導(dǎo)出波形的旅行時(shí)。導(dǎo)出的旅行 時(shí)可以作為質(zhì)量控制參數(shù)， 稍后可以將這個(gè)旅行時(shí)覆蓋在波形上來形象化的檢查 處理結(jié)果的質(zhì)量。導(dǎo)出旅行時(shí)的模塊需要處理出來的慢度記錄以及鉆孔直徑記錄圖， 這些資料 都必須與參考深度匹配。1、關(guān)閉所有打開的窗口，保存所做的改動(dòng)。2、選擇 Petrophysics > Full Waveform Sonic > Post-Processing > Traveltime打開旅

46、行時(shí)模塊3、如圖 37所示，輸入以下參數(shù)（如圖 38所示） Input/Output SetDIPOLE5、打開曲線模板 swp_traveltime.layout DT CALIDT_PICK1_1203WIRE.HCAL_1 DT_F圖38 ：建立旅行時(shí)模塊4、運(yùn)行模塊圖39：第一個(gè)接收器計(jì)算的旅行時(shí)(從處理得到的慢度中導(dǎo)出)2.4.4 相關(guān)性顯示相關(guān)性圖形是由慢度時(shí)間相關(guān)分析( STC)得到，如果 STC處理有差錯(cuò)，可 通過這個(gè)圖件做檢測(cè)分析。1、打開 swp_semblance.layout .2、打開陣列聲波窗口改變下列值：從FILTERED03集合增加聲波陣列記錄 (WF21 -

47、WF28)在顯示的 Display Waveform(s) Plus 里，選擇 Semblance ，并且選擇記錄 為FILTERED03.SEMBLANCE。選擇the FILTERED03.SEMBLANCE log結(jié)果如圖 40所示。圖40 ：在陣列聲波數(shù)組中顯示的相關(guān)性圖3、菜單欄 WELL下點(diǎn)擊 SAVE。4、關(guān)閉所有打開的窗口2.4.5 陣列聲波重處理TRANSM波形將按照一定間隔被重處理 (在后期處理中) 來增強(qiáng)垂直分辨率或者補(bǔ)償 井眼的直徑差值。 重處理是基于 multiple-shot processing 技術(shù)，它包含了由于 與工具位置點(diǎn)重合而帶來的冗余信息。 選擇若干子陣列， 這樣子陣列的大小可以 覆蓋住調(diào)查區(qū)的地層厚度。有兩個(gè)方法 (configurations) 來進(jìn)行處理， RECEIV和TRANSM ，RECEIV使用的是共源極配置，可以用來提高垂向分辨率；使用的是共接收配置。 由這兩個(gè)配置得來的平均慢度用于補(bǔ)償因井眼大小改變而 帶來的影響?？墒褂媚K Calculate Mean 完成。提示： 當(dāng)信號(hào)波長非常長時(shí)，這個(gè)方法是無效的