Geolog全波列聲波測井中文手冊

1、Geolog軟件技術手冊 Full Sonic Wave Processing -SWB帕拉代姆公司北京代表處2006年12月1、綜述- 1 -1.1 預備知識- 1 -1.2數據- 1 -2、陣列聲波全波形- 2 -2.1數據準備- 3 -查看/創(chuàng)建一個聲波列陣工具模版- 3 -2.1.2 練習指導2-創(chuàng)建其他波形屬性- 5 -波形分解- 6 -2.1.4深度轉換- 7 -2.2 處理- 8 - 數據分析- 8 -去噪- 11 -2.2.3 設計濾波器- 17 -2.2.4 振幅恢復- 19 -2.3陣列聲波處理- 20 -處理模塊簡介- 21 -偶極波形處理- 21 -2.3.3 單極波形

2、處理- 23 -2.3.4 拾取標志波至- 26 -2.4后期處理34綜述34頻散校正352.4.3 傳播時間疊加372.4.4 相關性顯示382.4.5 陣列聲波重處理403、機械性質443.1綜述443.2 計算動力學彈性性質44附錄I-快速運行45附錄II-頻散校正討論461、綜述歡迎閱讀Geolog軟件SWB指導教程。這是帕拉代姆公司Geolog軟件SWB模塊的使用教程，通過典型工作流程和程序模塊的使用，使用戶熟悉陣列聲波處理過程的基本工具的使用。l 熟悉陣列聲波波形l 查看/創(chuàng)建一個陣列聲波工具模版l 創(chuàng)建一個給定波形的其他屬性l 解裝測井記錄l 預處理及波形數據分析l 使用平均深度

3、、平均時間和頻率濾波器去噪l 去除數據采集中振幅增益和標準化l 處理偶極和單極聲波測井儀下波形l 自動和交互式的拾取波至l 頻散修正l 傳播時間覆蓋l 查看結果的相關性圖形1.1 預備知識Geolog基本測井記錄中文檔，所有的輸入文件為Bold Courier New，輸出為Bold Courier New，不加粗。1.2數據教程中使用了下面的附加文件（文件并非由軟件提供）數據：geophysics_master.unl曲線模板：(copy from layouts) swp_array_sonic_comparison swp_dispersion swp_depth_average swp

4、_frequency03 swp_frequency05 swp_other_attributes swp_projection swp_raw_sonic_waveforms swp_reprocessing swp_semblance swp_traveltime函數：圖像：井：polaris_03, polaris_042、陣列聲波全波形典型的陣列聲波工具在接收剖面上有8個接收器和3個不同的發(fā)射器。也要注意的是，不同服務公司出產的陣列工具也不相同。在數據采集上，基于不同的調查研究目的采用不同的工作模式。在Geolog里，“WF”是慣用的波形記錄名字后跟兩位數字，首位數代表操作模式索引，第

5、二位數表示接收器索引。例如，WF21代表波形記錄為模式2的第一記錄道（基于偶極模式）。每一個接收器的聲波波形以數組記錄存載。用Geolog layout(模版)可以在圖像道和陣列道顯示一個聲波波形。在這兩個道里，水平軸表示信號的旅行時，垂直軸為參考深度。當選定一個圖像道后，用戶可以打開一個聲波陣列查看窗口，對數據做進一步分析。圖1：原始聲波波形顯示圖1為原始聲波波形顯示圖像道（第一道和第三道），陣列道（第二道和最后一道）。最初的兩道顯示的是第一接收器的偶極波形記錄，最后兩道是第一接收器的單極波形記錄。2.1數據準備通過本節(jié)可以熟悉聲波全波形數據準備l 查看/創(chuàng)建一個列陣聲波工具模版l 創(chuàng)建給定

6、波形的其他屬性l 手動解裝波形l 進行深度轉換2.1.1查看/創(chuàng)建一個聲波列陣工具模版在波形處理前先了解工具的相關信息，它是以文件夾形式存儲。斯倫貝謝公司的DSI測得的井資料在Geolog中已經存儲在Site目錄下，作為四個缺省文件（dsi_m1 to dsi_m4）。當有改動或者初次運行此工具，相關信息會保存在本工具文件夾下的specs路徑中。) 工具模版的數據信息可以在Text中的Constants 表格中看到。1、啟動Geolog軟件打開STARS工區(qū)2、點開 Well 3、點擊well>open，打開Polaris_034、選擇Petrophysics > Full Wav

7、eform Sonic > Data Preparation > Create Tool Specification.模塊顯示斯倫貝謝公司的spec文件dsi_m2，給定的為缺省值，如下圖2圖2：array_sonic_create_toolspec 中dsi_m2的缺省值Spec文件中的參數如下表所示：參數描述TOOL_NAME聲波采集儀器的名稱如：DSI, XMACTOOL_SPEC聲波工具模版，以模版文件命名TOOL_MODE聲波工具的工作模式，通常在Text>Constants表格中定義TOOL_R聲波儀器的直徑SAMPLE_RATE聲波采樣率，通常以微次生法測量TR

8、_DISTANCE發(fā)射器與第一接收器的距離X_INTERV兩個鄰近接收器之間的距離DEPTH_SHIFT工具深度的平差不同于工具的參考深度和井的參考深度。負的值表示工具參考深度在井的參考深度之下5、如果需要的話可以對缺省值做調整6、點擊Start運行模塊一個新的模塊文件在本工區(qū)的SPECS目錄下產生了。7、選擇Launcher > Close來關閉陣列聲波工具模塊8、檢查新模塊dsi_m2的值，如圖3所示圖3：dsi_m2的存儲位置以及值提示：不管Geolog中使用的量度制，米制的或是特定的，聲波陣列工具模版總是會保存為米制的。當要運行這個工具，首先請確認工具模版的值，以確保得到正確的結

9、果。2.1.2 練習指導2-創(chuàng)建其他波形屬性聲波波形可以用反射振幅正常顯示。頻率譜、瞬時相位、和瞬時振幅屬性可以從解釋角度很好的量化突出這些特性。原始聲波頻率譜顯示的是不同到達時的主要頻率范圍，這樣可以在頻率域查看各到達模式的屬性特征。瞬時相使弱信號加強，但是也同樣使噪音增強。瞬時振幅是包絡的振幅，包含更多側向變化，但是降低了垂向分辨率。tp_array_sonic_attribute模塊用于計算一個給定波形的這些屬性。三個測井記錄可以從這個模塊輸出-頻率譜、瞬時相位以及瞬時振幅，可以在圖像道或者陣列道中查看這些屬性。1、選擇Petrophysics > Full Waveform So

10、nic > Data Preparation > CreateOther Attributes來顯示tp_array_sonic_attribute模塊。2、如圖4所示，對窗口做如下值改變 Input Set改為SONIC Output Set改為SONIC_ADJ LOG_IN類型為WF21 FREQ_OUT類型為WF21FREQ INST_PHASE類型為WF21PHASE INST_AMP類型為WF21AMP圖4: 改變屬性模型參數3、點擊Start運行4、打開曲線模板swp_other_attributes.layout來查看結果，如圖5所示。原始數據 原始數據頻率譜 原始

11、數據瞬時相位 原始數據瞬時振幅圖5：波形屬性顯示2.1.3波形分解來自每個接收器的聲波波形有時被疊加記錄到一起。這個模塊用于把波形分解為每個接收器的單個記錄。通常，每個模塊有自有的測深記錄，這個記錄要作為輸出的參考深度。在這個模塊里，需要提供分解波形的名字前綴。要遵守Geolog的命名約定。為了做示范，本例中的未分解波形都刻意不進行自動分解。1、選擇Petrophysics > Full Waveform Sonic > Data Preparation > LogUnpacking來顯示tp_array_sonic_unpack模塊2、如圖4所示，對窗口做如下值改變 Inp

12、ut Set改為SONIC_PACKED PREFIX類型為WG DEPTH選擇PWD2 LOG_WF選擇PWG2圖6: 改變分解波形模型參數3、點擊Start運行4、將PREFIX值變?yōu)閃N，LOG_WF值變?yōu)镻WN2，再運行模型5、打開Text，檢查未分解記錄值。如下圖7所示圖7：分解波形輸出記錄2.1.4深度轉換DSI工具數據的參考深度與井記錄的參考深度不同，差值從測量工具的底部到接收器剖面的中段進行修正，如上節(jié)的depth-shift。1、選擇曲線道模版swp_other_attributes2、選擇Tools > Depth Shift > Apply Logs3、如圖8

13、所示，對窗口值做如下改變 Input Set 改為SONIC Output Set 改為SONIC_DEPTH_SHIFT DEPTH_OFFSET 輸入參數為-7.9248 LOG_IN 選擇波形WF21_1為WF28_1 WF41_1為WF48_1 CORRECTED 選擇DEPTH圖8： 深度轉換模型4、點擊Start運行5、關閉所有窗口。2.2 處理l 通過本節(jié)可以使用戶熟悉全波形聲波處理流程l 在聲波數組查看窗口操作和解釋數據l 使用平均時間和平均深度以及Butterworth頻率濾波器去除波形噪音l 設計一個濾波器l 去除數據采集中的振幅增益和標準化 數據分析為了保證處理結果質量，

14、在數據處理前要檢查數據確保無誤。波形數據通常包含噪聲。通常需要選擇適當的工具來盡可能多的去噪。1、打開模板swp_raw_sonic_waveforms.layout。2、選擇WF21的圖像道，選擇View > Array Sonic View。數組聲波查看窗口顯示波形WF21以及它當前深度頻率域的位置。（見圖9）) 要打開陣列聲波查看窗口，需要選擇打開帶數據的圖像道。陣列聲波查看窗口圖9：圖像道和聲波數組窗口陣列聲波窗口用來詳細的分析波形，它有兩個顯示區(qū)域。顯示區(qū)在當前深度位置顯示信號波形，可以估算出時間間隔（接收器的時間延時）。頻率域顯示區(qū)顯示每個波形相應的的頻率范圍。3、點擊聲波數

15、組記錄- Log Select 圖標打開可選擇多條曲線。提示：挑選出來的記錄必須來自同一個集合和模型工具（例如：有相同的模型索引-字符名字后的第一位數相同）在圖像道上顯示。4、選擇記錄SONIC.WF21_1 to SONIC.WF28_1，點擊OK。多重曲線以及其對應的頻率域在數組聲波的窗口顯示。5、在模版上第一圖像道(SONIC.WF21)點擊任何位置。聲波數組顯示窗口的信號在每次選擇深度點的時候會刷新。) 另一個方式是通過深度區(qū)域，鍵入范圍值或者使用上下按鍵。6、設置深度點和時間范圍如下所示：在顯示區(qū)放大或者改變時間范圍7、將光標定位在信號顯示區(qū)任意位置，按住鼠標左鍵拖曳定義一個矩形區(qū)域

16、，放開鼠標左鍵。通過顯示的時間范圍值來反映被選擇的區(qū)域。) 將時間范圍值的起始和終點值輸入與用鼠標改變信號顯示區(qū)域的效果是一樣的。8、再次點擊鼠標左鍵使區(qū)域變?yōu)樵瓉泶笮　Ｔ诓ㄐ紊瞎浪銜r間間隔：提示：這個功能只在多重波形在當前活動的時候起作用。)當時間間隔固定時，功能不可選時首先將信號顯示區(qū)放大。9、固定時間間隔，在慢度剖面顯示。10、顯示區(qū)，在第一個接收器（WF21）的信號，按住鼠標左鍵，拖至到最后一個接收器的信號處（WF28），放開鼠標左鍵。如圖10粗淡紅色的陡線顯示第一個和最后一個接收器初至時間間隔。淡紅色顯示估算時間間隔，波至慢度及接收間隔圖10：時間間隔估算8個接收器的信號以及頻率域在

17、當前選擇的深度顯示。11、不選Moveout復選欄，不顯示時間間隔數據。12、關閉陣列聲波顯示窗口。去噪波形圖總是會包含噪聲，太多的噪聲會不能識別可檢測的波至。預處理是提高數據質量保證初至被識別的重要步驟。Geolog 有一系列基本的預處理工具可以滿足一個波形的任意去噪要求。l 去除時間均值，在時間方向使用平均時濾波器l 去除深度均值，在深度方向使用平均深度濾波器l Butterworth 濾波器，是一個頻率濾波器根據具體的需求，用戶可決定針對某個波形使用對應的濾波器（1）去除時間均值噪聲這個模塊提供沿時間方向去除噪聲濾波器。先在每個窗口計算出平均值，然后從原始波形移除噪聲。) 這是用于移除低

18、頻噪聲背景的濾波器。1、選擇 Petrophysics > Full Waveform Sonic > Pre-Processing > RemoveTime-Average.2、如圖10，對打開窗口中值做如下修改。 Input Set 改為SONIC_DEPTH_SHIFT UTIME WINDOW 例如，設置濾波窗口值為520 in micro-seconds (usecs)忽略采樣場。（只是在Geolog不能識別采樣率時需要采樣場）改變窗口的取值，多次運行這個模塊以取得合適的結果。 LOG_IN 設為WF21（可以同時導入多個文件）圖11：時間平均濾波器參數設置3、點擊

19、Start運行模塊4、打開繪圖模版swp_array_sonic_comparison.layout。有低頻噪聲的原始圖像運行濾波器后的結果圖圖12：用時間平均濾波器濾波前后對比圖（2）去除深度均值噪聲這個模塊提供沿深度方向去除噪聲濾波器。先在每個活動窗口計算出平均值，然后從原始波形移除噪聲。) 這個濾波器可去除典型的諸如套管信號或者儀器帶來的噪聲。1、打開井Polaris_04，保存對先前井所作的改變。2、關閉所有打開窗口。3、選擇Petrophysics>Full Waveform Sonic>Pre-Processing> RemoveDepth-Average.4、在

20、tp_array_sonic_depth_average窗口對值做如下改變（如圖13） Input Set改為SONIC WINDOW設置濾波器長度為21 LOG_INWF41圖13：設置平均深度濾波器模型值5、點擊Start運行模型6、打開繪圖模版swp_depth_average.layout。套管和儀器形成的信號可以很強，使得重要的波至難以識別。圖14左邊道顯示的波形是一個這樣的例子，直直的垂線占主要部分。右邊的道顯示的是不需要信息被移除后的清晰的波至。運行濾波器后的結果圖有低頻噪聲的原始圖像圖14：用深度濾波器濾波前后對比圖（3）使用頻率濾波器Butterworth 濾波器是一個頻率濾

21、波器，可以通過給定一個高的和低的截止值以及創(chuàng)建濾過順序來實現。Butterworth 濾波器選擇的順序（從1到8）決定了從通頻帶到抑制帶的過渡陡度。單位設為千赫茲而不是赫茲，是為了便于進行聲波波形的處理。如果數據是來自于標準的DSI工具采集，可以忽略采樣率，因為Geolog可以自動地讀取。有關濾波器的調整適應更多信息參考菜單選項下的Design Filter。1、保存對Polaris_04的改動，打開井Polaris_03。2、選擇Petrophysics > Full Waveform Sonic > Pre-Processing > Frequency Filter.3、

22、對打開窗口做如下改變（如圖15所示） Input Set SONIC_DEPTH_SHIFT Output Set 改為FILTERED03 HIGH_END 設為3HZ(移除所有高于這個值之外的頻率信號) SAMPLE_RATE10 LOG_IN 從輸入集合選擇8個波形：SONIC_DEPTH_SHIFT.WF21_1 to WF28_1圖15 ：改變butterworth濾波器的頻率范圍和大小4、點擊Start運行這個模型5、打開swp_frequency03.layout6、選擇View > Array Sonic View來查看濾波后圖形7、設置當前深度到如圖16所示。在3HZ以

23、上的頻率已經在當前深度位置被移除。圖16：WF濾波后圖形8、關閉陣列聲波查看窗口。9、選擇tp_array_sonic_filter模塊，重復第二步和第三步改變頻率范圍為0-5赫茲，輸出名為FILTERED05。10、點擊Start運行模塊11、打開swp_frequency05.layout12、打開聲波數組查看窗口。大于5Hz的頻率信號都被移除，如圖17所示，在當前深度位置顯示的頻率譜。圖17：濾除高于5Hz波后的WF21波形圖13、關閉陣列聲波窗口。 設計濾波器Butterworth濾波器將作為頻率濾波器用在下面的例子中。Butterworth 是一個很好的帶通濾波器，在損失了部分頻率域

24、后有較為合理的時間衰減。高頻和低頻截止值用于帶通和帶止。對于低通濾波器，高頻波被截止，反之亦然。Butterworth 濾波器的排列控制著在通頻帶之外信號變弱的快慢。設計一個合適的濾波器或是子波，必須要注意采樣率和頻率值。設計一個Butterworth子波，采樣率（SR）要使用0.001秒，這樣可以得到Nyquist frequency = 1 / (2 * SR) = 500 Hz。輸入頻率必須要小于尼奎斯特頻率（Nyquist frequency）頻率。也就是說，如果提供的值大于尼奎斯特頻率（Nyquist frequency），采樣率也要更好。1、選擇Petrophysics >

25、Full Waveform Sonic > Pre-Processing > DesignFilter.2、在filter_create窗口對值做如下改變（如圖18） FILTER_TYPE 改為Butterworth FILTER_NAME bw1040 LOW_START 1000 HIGH_END 3000 PLOT_SR 0.00004 DISP_SCALE AMPLITUDE圖18：濾波器創(chuàng)建模塊3、點擊Start運行這模塊。4、查看創(chuàng)建的濾波器（如圖19）。圖19：運行創(chuàng)建濾波器模塊后的結果 振幅恢復在數據采集中，振幅增益以及標準化是為了增強波形信號的質量。可基于有效的

26、增益和標準化因素（歸一化因素），使波形回到原來形狀。默認的，去除標準化因素選擇inverse；去除增益選擇為multiply。1、選擇Petrophysics > Full Waveform Sonic > Pre-Processing > Amplitude Recovery.2、在tp_array_sonic_tar窗口對值做改變（如圖20所示） Input Set 改為SONIC_DEPTH_SHIFT NORM_OPT 設為INVERSE GAIN_OPT 設為MULTIPLY REFERENCE 選擇UNPACK.DEPTH LOG_WF 選擇WF21_1至WF28

27、_1波形 LOG_NORM 選擇井UNPACK.WN21_1 到 UNPACK.WN28_1 LOG_GAIN 選擇井UNPACK.WG21_1 到UNPACK.WG28_1圖20：振幅恢復模塊3、點擊Start運行模型。4、打開text窗口在原來的集合中查看井數據（如圖21所示）。圖21：Text窗口查看去除振幅增益后的波形記錄2.3陣列聲波處理通過本節(jié)，用戶可以：l 處理偶極波形l 處理單極波形l 使用自動交互式拾取波至處理模塊簡介聲波波形處理基于慢度時間相關法。這個模塊在活動窗口中通過時間和慢度方向可以查看多個聲波。默認的，時間范圍可以包括所有的記錄時間；慢度的范圍至少包括目的區(qū)。例如，

28、慢度范圍為100800 us/m (1200-10000 m/s)可以包含大多數區(qū)域的縱波波至。建議用較大的慢度范圍這樣可以盡可能包含所有波至，但這也帶來更多的計算量。) 首先選擇少量數據來決定時間和慢度方向的最優(yōu)選擇。要注意慢度步長的選擇，過大的步長會降低最終聲波分辨率。在每個深度點計算出的最大值作為與深度相關的函數映射到慢度軸上；這是最終輸出Projection log。這個模塊有很多控制參數，詳細參數解釋可以查看模塊幫助。偶極波形處理在本節(jié)，被處理的波形filtered03 和filtered05將用于下文的離散。1、保存井文件Polaris_03。2、關閉所有查看窗口。3、選擇Petr

29、ophysics> Full Waveform Sonic > Processing > Waveform Processing來打開tp_array_sonic_process窗口。4、設置輸入參數（如圖22所示） Input Set SONIC_DEPTH_SHIFT Output Set DIPOLE TOOL_SPEC dsi_m2 START 0 STOP 10000 SLOWNESS_FROM 35 SLOWNESS_TO 500 SLOWNESS_STEP 4 LOG_IN 由SONIC_DEPTH_SHIFT集合，選擇8個波形： WF21_1 to WF28_

30、1圖22：偶極波形處理模型窗口5、點擊Start運行模型。當模型運行結束，結果以相關性圖件顯示，如圖23所示。數據在當前深度位置顯示（當深度未規(guī)定時是在井記錄的頂部）圖23：相關性圖- WF21 - WF28波形處理結果圖6、保存相關性圖為semblance_plot_01。)將相關性圖以唯一文件名保存，便于以后查看，否則圖件會被覆蓋。7、改變輸入和輸出集合的名字為FILTERED03，重復4-6步，將相關性圖形保存為semblance_plot_03。8、改變輸入和輸出集合的名字為FILTERED05，重復4-5步，將相關性圖形為semblance_plot_05。9、關閉相關性圖。10、保

31、存井文件。 單極波形處理11、在tp_array_sonic_process窗口，設置輸入參數，如圖24所示： Input Set SONIC_DEPTH_SHIFT Output Set MONOPOLE TOOL_SPEC dsi_m4 START 0 STOP 5120 UTIME_WINDOW 400 UTIME_STEP 50 SLOWNESS_FROM 30 SLOWNESS_TO 298.7 SLOWNESS_STEP 1.98 LOG_IN 從SONIC_DEPTH_SHIFT集合，選擇8個波形WF41_1to WF48_1圖24：單極波形處理模型窗口12、點擊Start運行模

32、塊。當模型運行結束輸出一個映射和一個相關性圖件，圖件在每個深度點存儲的是2D 相關性圖像， 可以通過陣列聲波數組查看窗口進行質量控制。如下圖所示。圖25：陣列聲波數組窗口查看圖件13、關閉相關性圖件，不保存改變。14、打開模版swp_projection.layout。上面例子中的偶極和單極子映射記錄如圖26所示。偶極波形處理后圖形 單極子波形處理后圖形 水平軸顯示慢度(START)到(END)圖26：深度相關函數映射15、在WELL下點擊SAVE. 拾取標志波至映射記錄在慢度-深度平面內代表了一個或多個相關性。它們與某些波至有較好的一致性。模型在輸入的映射記錄中自動或者交互式的尋找相關最大值

33、位置，從慢度投影記錄拾取相應一組波至。搜索從一個起始值開始（初始慢度值在映射記錄里）來標定某一初至，搜尋范圍以從映射記錄估計的起始值為中心。估算輸入起始值最簡便方法是將投影記錄插入到圖像道，打開陣列生波可視窗口，移動鼠標標定出信號顯示區(qū)的目的位置，光標位置X的值會在位置對話框的狀態(tài)欄顯示出來。這個值就是從圖上讀取的慢度值，它可以用來作為搜索起始值。對于交互式拾取，測井曲線首先是自動被拾?。▍⒁奀urve Insert聯(lián)機幫助），編輯，最后作為備用初始記錄被保存。主要有：1、自動拾取一個偶極測井記錄2、自動拾取一個單極測井記錄3、對于第二步中不合適點做交互式拾取。自動拾取一個偶極測井記錄1、選擇

34、Petrophysics > Full Waveform Sonic > Processing > Picking來打開the tp_array_sonic_label窗口。2、選取輸入起始值(SLOWNESS_S1)。讓swp_projection為活動窗口，選擇偶極投影道來打開陣列聲波可視窗口。在信號顯示區(qū)移動鼠標到拾取位置，讀取顯示出的X值。這個值即是從圖上讀取的慢度值，可作為起始值（SLOWNESS_S1），如例子中圖27示。關閉陣列聲波可視窗口。圖27：確定搜索起始值3、選擇模型窗口，設置值如下圖28所示 Input/Output Set選擇DIPOLE，這是處理后

35、的輸出集合 SLOWNESS_WINDOW 本例中設為121 us/f（因為僅有一個波至，對于偶極使用一個大點的范圍是可靠的）窗口長度必須包括波至的最大波動幅度，但是也不能過大要注意排除鄰區(qū)波至。 SLOWNESS_S1在前一步操作中，位置對話框里確定的X值，必須要接近值121。圖28偶極自動拾取模型設置4、運行模型5、重復地3-4步二次，改變輸入輸出集合名字分別為FILTERED03，FILTERED05。6、設定swp_projection為活動窗口，查看曲線道（如圖29示）圖29：自動拾取的起始記錄覆蓋在偶極映射圖像上7、打開屬性對話框，注意DIPOLE.DT_PICK1井的左右邊界為處

36、理的慢度起始和終止值（35 and 500 US/F），輸入模型中。8、關閉對話框。單極投影記錄自動拾取9、確定輸入tp_array_sonic_label模型所需要的的起始值(SLOWNESS_S1 and SLOWNESS_S2)選擇單極圖像道打開陣列聲波可視窗口。在信號顯示區(qū)使用鼠標定位，在定位對話框確定兩個初始值，如圖30所示為確定第二波至的例子。圖30：確定SLOWNESS_S2的輸入初始值10、使tp_array_sonic_label為活動窗口，對模型值做如下設置（如圖31所示）： Input / Output Set 選擇MONOPOLE集合，它是處理中的輸出集合 Slowne

37、ss_Window 設為15.24 us/f為搜索窗口大小 PICK_OPTION BOTH SLOWNESS_FROM30 SLOWNESS_TO 298.7 SLOWNESS_STEP 1.98SLOWNESS_S1 60SLOWNESS_S2 85圖31：交互式拾取兩個井記錄模型界面11、運行模型。12、選擇模板swp_projection作為查看模板，如圖32所示。圖32：自動拾取井記錄覆蓋在單極圖像道上13、菜單欄WELL下點擊SAVE Interactive Picking交互式拾取14、用swp_projection模板，在單極投影道圖像道上選擇拾取記錄DT_PICK2。15、選

38、擇Tools > Curve Insert。16、手動的重新拾取不合適的橫波波至。17、打開文本查看重命名被編輯的井記錄DT_PICK2 version 2為DT_SEED。18、菜單欄WELL下點擊SAVE。19、選擇tp_array_sonic_label模型，對值做如下改變（如圖33）： PICK OPTION SINGLE SLOWNESS_S1 使用DT_SEED作為初始記錄 DT_PICK1 使用DTS_PICK圖33：單井交互式拾取模型界面20、運行模塊21、在曲線模板中查看結果，使用圖34作為參考。如果效果不夠好，調整搜索范圍（例如：SLOWNESS_WINDOW和SLO

39、WNESS_S1值），重復先前的步驟。圖34：偶極和單極的投影記錄的自動式和交互式拾取顯示在單極投影區(qū)，綠色的曲線表示作為初始記錄的自動拾取曲線(DT_PICK2)。藍色曲線是基于初始子記錄的交互式拾取(DTS_PICK)。提示：處理的慢度記錄需要進行深度轉換，這涉及到井的參考深度。在Geolog中不會自動轉換，因為對于參考深度工具，沒有足夠的信息可供選擇。2.4后期處理通過本節(jié)，用戶將熟悉后期處理流程l 計算并應用頻散校正，確保每個井記錄在同一參考系下l 應用頻散校正l 導出旅行時l 查看相關性圖形顯示l 陣列聲波數組重處理，計算平均慢度綜述后期處理以及后期質量控制是聲波處理的重要部分。從處

40、理過的記錄中導出的旅行時可以作為質量控制來描述聲波波形波至。頻散以及井眼情況等將影響最終處理結果。頻散校正模型計算出頻散曲線，在每個深度位置使用校正值來修正頻散影響。波形將按照一定間隔被重處理（在后期處理中）來增強垂直分辨率或者補償井眼較大的直徑差值。頻散校正這個模塊輸出修正后的慢度記錄。它可以通過覆蓋陣列道隨意輸出頻散和頻譜記錄，也可以輸出壞道來標記有問題區(qū)域。1、菜單欄WELL下點擊SAVE。2、選擇Petrophysics > Full Waveform Sonic > Post-Processing > Dispersion Correction。提示：離散修正是一個

41、很費時的處理過程，本教程中為了實現這個操作只是用了少量資料來處理。3、如圖35所示，設置處理參數。圖35：設置處理范圍4、輸入參數如下：如圖36所示 Input SetSONIC_DEPTH_SHIFT Output SetRAW TOOL_SPECdsi_m2 START0 STOP10000 CALIWIRE.HCAL BSWIRE.BS RHOB_F1.24 G/C3(1240 K/M3) RHOB_SWIRE.RHOZ DT_F245.8 US/F (75 US/M) DT_PMONOPOLE.DT_PICK1 DT_SDIPOLE.DT_PICK1LOG_IN從 SONIC_DEPT

42、H_SHIFT 集合, 選擇8個波形WF21_1 to WF28_1圖36：頻散校正窗口5、運行模塊) 處理速度取決于計算機的性能和數據量的大小。6、重復第三和第四步，改變值。 Input /Output setFILTERED03 DT_SFILTERED03.DT_PICK17、重復第三和第四步，改變值 Input /Output setFILTERED05 DT_SFILTERED05.DT_PICK18、打開模板swp_dispersion.layout。提示：離散修正結果如圖37所示，處理范圍為13150 - 13450 ft.圖37：離散結果對比實例（數據為13150 - 1345

43、0 ft）1,2,3道分別表示覆蓋在疊加頻率域的波形離散曲線，4道表示低頻濾波后的波形，5道波形大于5 Khz的頻率被濾過，6道大于3 Khz頻率被濾過。對原始波形做處理時，發(fā)現在做頻散校正后，有不穩(wěn)定的波動現象。如上例，范圍如箭頭所示，通過0-5 KHz帶通濾波處理，修正的剪切波慢度穩(wěn)定，同時1-3 KHz濾波器給出的離散都很清晰。9、菜單欄WELL下點擊SAVE。 傳播時間疊加在前面步驟中得到的慢度記錄，可以從中導出波形的旅行時。導出的旅行時可以作為質量控制參數，稍后可以將這個旅行時覆蓋在波形上來形象化的檢查處理結果的質量。導出旅行時的模塊需要處理出來的慢度記錄以及鉆孔直徑記錄圖，這些資料

44、都必須與參考深度匹配。1、關閉所有打開的窗口，保存所做的改動。2、選擇Petrophysics > Full Waveform Sonic > Post-Processing > Traveltime打開旅行時模塊。3、如圖37所示，輸入以下參數（如圖38所示） Input/Output SetDIPOLE DTDT_PICK1_1 CALIWIRE.HCAL_1 DT_F203圖38：建立旅行時模塊4、運行模塊5、打開曲線模板swp_traveltime.layout。圖39：第一個接收器計算的旅行時（從處理得到的慢度中導出） 相關性顯示相關性圖形是由慢度時間相關分析（ST

45、C）得到，如果STC處理有差錯，可通過這個圖件做檢測分析。1、打開swp_semblance.layout.2、打開陣列聲波窗口改變下列值：從FILTERED03集合增加聲波陣列記錄(WF21 - WF28)。在顯示的Display Waveform(s) Plus里，選擇Semblance，并且選擇記錄為FILTERED03.SEMBLANCE。選擇the FILTERED03.SEMBLANCE log結果如圖40所示。圖40：在陣列聲波數組中顯示的相關性圖3、菜單欄WELL下點擊SAVE。4、關閉所有打開的窗口。 陣列聲波重處理波形將按照一定間隔被重處理（在后期處理中）來增強垂直分辨率或

46、者補償井眼的直徑差值。重處理是基于multiple-shot processing 技術，它包含了由于與工具位置點重合而帶來的冗余信息。選擇若干子陣列，這樣子陣列的大小可以覆蓋住調查區(qū)的地層厚度。有兩個方法(configurations)來進行處理，RECEIV和TRANSM，RECEIV使用的是共源極配置，可以用來提高垂向分辨率；TRANSM使用的是共接收配置。由這兩個配置得來的平均慢度用于補償因井眼大小改變而帶來的影響?？墒褂媚KCalculate Mean完成。提示：當信號波長非常長時，這個方法是無效的，厚度需要大于波長的1/4。這個技術也易于帶來噪聲，較少的接收器刪除的噪聲也少。RE_

47、PROC可以開啟和關閉，一旦打開重處理，子波列接收器NUM_R在subarray，即定義的地層厚度也可以輸入，但是圖像和相關性輸出是不可用的。1、選擇Petrophysics > Full Waveform Sonic > Processing > WaveformProcessing打開tp_array_sonic_process窗口。2、選擇Launcher > Datum設置處理參數Depth 13123.35 - Depth 13189.00 ft (4000 - 4020m)3、鍵入輸入參數如下圖41所示： Input Set SONIC_DEPTH_SHIF

48、T Output Set DIPOLE_RECEIV4 TOOL_SPECdsi_m2 START 0 STOP 8000 RE_PROC Yes1.Hsu, K. and Cheng, S., 1987; X.M.Tang and A.Cheng, 2004 METHOD RECEIV NUM_R 4 LOG_IN 選擇波形SONIC_DEPTH_SHIFT WF21_1 to WF28_1圖41：使用RECEIV方法重處理聲波波形4、運行模型。5、改變輸出集合為DIPOLE_TRANSM4，方法改為TRANSM，重復以上步驟。6、關閉相關性圖，不保存。7、對創(chuàng)建的投影記錄執(zhí)行自動拾取。選擇

49、Petrophysics > Full Waveform Sonic > Processing > Picking來打開tp_array_sonic_label窗口。改變輸入/輸出集合為DIPOLE_RECEIV4使用上步的慢度值(SLOWNESS_S1) 121運行模型改變輸入/輸出集合為DIPOLE_TRANSM4運行模型8、打開layout swp_reprocessing，如下圖42圖42：RECEIV和 TRANSM方法得到的慢度記錄對比圖9、選擇Petrophysics > Full Waveform Sonic > Post-Processing > Calculate Mean來打開tp_array_sonic mean模塊。10、鍵入輸入參數，如圖43所示 Input Set DIPOLE_RECEIV4 Output Set DIPOLE_MEAN LOGIN1 DT_PICK1 LOGIN2 D