李海明，劉穎華，侯伯剛

(1.承德石油高等?？茖W校 數(shù)理部，河北 承德 067000；2.中國石油勘探開發(fā)研究院,北京 100083)

測井資料在井震聯(lián)合標定中的應(yīng)用

李海明1，劉穎華1，侯伯剛2

(1.承德石油高等?？茖W校 數(shù)理部，河北 承德 067000；2.中國石油勘探開發(fā)研究院,北京 100083)

井震聯(lián)合標定地震反射目的層是構(gòu)造解釋、巖性反演中最基礎(chǔ)的工作，是連接地震、測井數(shù)據(jù)的橋梁。建立目標層位與地震反射之間對應(yīng)關(guān)系的工具是利用測井資料合成地震記錄。從合成地震記錄的原理及制作方法入手，指出影響井震聯(lián)合層位標定質(zhì)量的關(guān)鍵因素是正確處理測井數(shù)據(jù)和準確提取地震子波。利用杏樹崗油田杏六東Ⅰ區(qū)塊x5-4-135井資料和該區(qū)地震資料試驗研究了測井曲線與地震子波對合成地震記錄質(zhì)量的影響。證明利用測井曲線提取確定性子波標定效果優(yōu)于Ricker子波和統(tǒng)計性子波。

井震聯(lián)合標定；測井數(shù)據(jù)；反射系數(shù)；子波

利用合成地震記錄這個地震資料解釋基本工具將地震資料與測井資料有機聯(lián)系起來，井震聯(lián)合精確標定目的層是構(gòu)造解釋的基礎(chǔ)[1，2]。地震剖面上層位標定的目的是賦予地震同相軸特定的地質(zhì)含義，建立地震同相軸與地層、地層組合之間對應(yīng)關(guān)系。也是開發(fā)階段利用地震資料和高密度井網(wǎng)資料開展精細構(gòu)造解釋、巖性解釋和預測，重構(gòu)地下地質(zhì)認識體系基礎(chǔ)[3]，高精度的合成地震記錄將決定著構(gòu)造解釋和儲層精細描述研究工作的質(zhì)量與水平。

1 合成地震記錄原理與方法

在不考慮干擾的情況下，制作合成地震記錄的數(shù)學模型為地層反射系數(shù)與地震子波的褶積[4]，公式如下：

x(t)=w(t)*r(t)

(1)

式中，x(t)表示地震道，w(t)代表地震子波，r(t)表示地下反射系數(shù)序列。

制作方法為：

1)速度測井曲線和密度測井曲線相乘得到波阻抗曲線；2)對波阻抗曲線進行深時轉(zhuǎn)換，由深度域轉(zhuǎn)到時間域；3)由波阻抗曲線轉(zhuǎn)換為反射系數(shù)曲線；4)提取地震子波；5)反射系數(shù)與子波進行褶積，得到井旁合成地震記錄道。

顯然反射系數(shù)的準確性和地震子波的精度主要地決定了人工合成地震記錄的精度。

2 合成記錄制作的關(guān)鍵

制作合成地震記錄，涉及測井曲線的編輯校正、測井曲線的方波化、拉伸控制、子波的選擇、低信噪比資料的標定、高分辨率資料的標定等多個方面[5]。但最為關(guān)鍵的是測井曲線的處理和地震子波的選擇。

2.1 測井曲線處理

2.1.1 測井曲線的編輯和濾波

測井曲線編輯和濾波是一項重要的基礎(chǔ)工作，測井曲線上存在著由于井徑變化、油氣與泥漿侵入、儀器因素等引起的干擾、孤立野值等[6]，對測井曲線的編輯和濾波就是要剔除和壓制這些干擾數(shù)據(jù)?？梢赃x擇的濾波有：中值濾波、帶通濾波和滑動平均。

2.1.2 測井曲線的方波化處理

測井曲線的高頻變化使得很難從測井曲線上找出地質(zhì)意義上的地層，通常意義上，無論地層的厚度如何，同一地層的測井數(shù)據(jù)應(yīng)該是大致穩(wěn)定的。測井曲線方波化就是根據(jù)測井曲線的變化特征，對測井曲線自動分層，在每一個分層內(nèi)，測井數(shù)據(jù)是相同的，在自動分層后，可以根據(jù)其它信息，對自動分層的結(jié)果進行人為的手工調(diào)整。用測井波阻抗曲線合成的地震記錄與井旁地震道比較，若反射時間、波組關(guān)系等基本接近，就可以判斷標定是正確的；否則標定就是錯誤的。

圖1給出的是一個測井曲線方波化處理的例子，在1 220 m之下的測井曲線進行了方波化處理，并利用方波化的測井曲線制作的合成記錄，對比合成地震記錄與井旁地震道波組對應(yīng)關(guān)系良好，說明處理是正確的。

2.1.3 時深轉(zhuǎn)換與反射系數(shù)計算

1)時深轉(zhuǎn)換：測井數(shù)據(jù)表示的是深度域的地層巖性信息(速度、密度、電性、孔隙等)；地震數(shù)據(jù)表示地下構(gòu)造和地下巖性在時間域的地震響應(yīng)。為了能使二者在時間域內(nèi)交互對比，就需要對測井資料進行時深轉(zhuǎn)換。直接利用聲波時差曲線積分進行時深轉(zhuǎn)換的公式為[7]：

(2)

式中：H0為起始深度，m；H為終止深度,m；T0起始深度對應(yīng)的時間,s；Th為深度H對應(yīng)的時間，s；Δth為深度對應(yīng)的聲波時差，ms。對不同深度的聲波時差測井曲線進行計算，就得到深時轉(zhuǎn)換關(guān)系。

2)反射系數(shù)：上下層介質(zhì)的波阻抗差異決定地層反射系數(shù)[8]，其表達式為：

(3)

式中，ρi，vi為測井曲線上第i層的密度和速度；ρi+1，vi+1為第i+1層的密度和速度；ρ，v來自于速度和密度測井曲線，因此測井資料對合成記錄有著很大的影響。

2.2 地震子波

地震子波是合成地震記錄最重要的影響因素之一，常用子波類型有：Ricker子波;由井口附近地震道數(shù)據(jù)提取的統(tǒng)計子波[4];確定性子波[6];時變子波。

制作合成地震記錄時，子波極性要與地震資料一致；子波頻率應(yīng)與井旁地震記錄主頻一致，需要對近井地震資料進行頻譜分析，確定主頻。

3 合成地震記錄過程中的拉伸問題及其原則

合成地震記錄過程中的拉伸就是將合成地震道與井旁地震道的標準層反射對齊。地震資料處理過程中的速度是通過速度分析獲得的。這種速度分析方法得到的速度與聲波測井得到的速度存在差異。一般而言，聲波測井的速度要高于地震資料分析得出的速度，因此對于同一套地層，合成記錄要比井旁實際地震記錄“短”一些。為了將兩者對應(yīng)起來，需要將合成記錄“拉長”一些。在標定過程中，拉伸從最淺的標志層開始，依次向下拉伸。

合成記錄制作過程中拉伸原則：1)用井徑曲線檢查聲波曲線可靠性，是否存在井眼垮塌，是否確實需要拉伸；2)拉伸壓縮時要遠距離調(diào)整，兼顧淺中深的各波組關(guān)系；3)避免局部小的調(diào)整；4)拉伸后，注意用慢度曲線監(jiān)控聲波曲線變形情況。

4 標定試驗與精度分析

依據(jù)上面的的分析，利用杏樹崗油田杏六東Ⅰ區(qū)塊x5-4-135井聲波測井資料和該區(qū)地震資料進行合成地震記錄和層位標定。密度使用Gardner公式ρ=0.31v0.25來求取[9]。

4.1 標定試驗

首先對井點附近地震資料的主頻進行了分析估計，目的層主頻在38～40 Hz左右，有效高頻能量大約為80 Hz，因此確定選用38 Hz的Ricker子波進行合成地震記錄，以最上面的主要標志層(T1)為標準，將合成地震記錄與井旁道對齊(見圖2)。

用井旁地震道提取統(tǒng)計性子波標定，試驗時地震道范圍選取為7×7，時窗長度250 ms。在仔細對比的基礎(chǔ)上，對合成地震記錄進行適度的拉伸對齊合成地震記錄與井旁地震道之間的波組關(guān)系，如圖3所示。

利用統(tǒng)計性子波對齊主要反射層之后，利用測井曲線提取確定性子波來合成的地震記錄(見圖4)。圖2、圖3、圖4給出的分別是利用Ricker子波、統(tǒng)計子波和利用測井曲線提取確定性子波合成記錄與井旁地震道的對比結(jié)果。

4.2 質(zhì)量分析

合成地震記錄與井旁地震道的相關(guān)性是評價合成地震記錄質(zhì)量的重要標準[5]，觀察圖2可發(fā)現(xiàn)，除最上面T1標志層對齊外，其他波組關(guān)系均沒有對齊，四個標定層位效果很差。圖3、圖4目視差別不明顯，合成記錄與井旁地震道波組對應(yīng)關(guān)系明顯好于圖2，四個標定層位中最上與最下兩個層位對齊，中間兩個層位均有可比較的明確含義。通過計算，x5-4-135井Ricker理論子波(38 Hz)合成地震記錄與井旁地震道相關(guān)系數(shù)約為69%，統(tǒng)計子波合成地震記錄與井旁地震道相關(guān)系數(shù)約為84%，而用測井資料與地震資料提取的確定性子波合成記錄與井旁地震道相關(guān)系數(shù)達88%，合成記錄質(zhì)量最高。

5 結(jié)論

1)提高井震聯(lián)合層位標定質(zhì)量的關(guān)鍵是測井曲線的正確處理和子波選用。2)試驗分析表明，利用測井曲線和井旁地震道提取的確定性子波進行標定精度最高。在開發(fā)階段利用地震資料和高密度井網(wǎng)資料開展精細構(gòu)造解釋時應(yīng)盡量選用確定性子波進行層位標定，提高相關(guān)研究工作的質(zhì)量與水平。

[1] 周永波，秦向紅，袁燕，等.淺談合成地震記錄制作中的制約因素及解決方法[J].工程地球物理學報，2013,10(4)：476-479.

[2] 吳奎，周東紅，吳俊剛. 合成地震記錄制作在油氣勘探中的應(yīng)用[J].海洋石油，2012,32(1)：28-32.

[3] 孫振濤，孟憲軍，慎國強，等.高精度合成地震記錄制作技術(shù)研究[J].石油地球物理勘探， 2002, 37(6)：640-643.

[4] 高國民，胡望水，黃玉欣.合成地震記錄子波提取方法[J].內(nèi)蒙古石油化工，2008，34(1)：130-132.

[5] 屈紹忠.合成地震記錄標定反射層位方法簡析[J].中國煤炭地質(zhì)，2013，25(7)：48-52.

[6] 劉明洋，郭韜，沈銘成，等.多道統(tǒng)計法提取子波及其應(yīng)用[J].石油天然氣學報，2010，32(5)：238-240.

[7] 李錄明，李正文.地震勘探原理、方法和解釋[M].北京：北京地質(zhì)出版社，2007.

[8] 劉中奇，邱廣筠.一種新的聲波測井曲線環(huán)境校正方法[J].測井技術(shù),2000,24(4):283-286.

[9] 董震,潘和平.聲波測井資料與地震屬性關(guān)系研究綜述[J].工程地球物理學報，2007,4(5):488-493.

Application of Logging Data in Well-Seismic Integration Calibration

LI Hai-ming1, LIU Ying-hua1, HOU Bo-gang2

(1.Department of Mathematics and Physics, Chengde Petroleum College, Chengde 067000, Hebei, China; 2.Research Institute of Petroleum Exploration & Development， PetroChina, Beijing 100083,China)

The well-seismic integration calibration of seismic reflection layer is the basis for structure explanation and lithological inversion, and it is also the bridge between seisme and logging data. Drawing seismic record from logging data is the only tool to set the relation between target strata and seismic reflection. This paper starts from the theory and method of Synthetic Seismograms, then points out that the crucial factor influencing well-seismic integration calibration is the correct handling of logging data and the exact extraction of seismic wavelet. This paper takes the logging data of Xingliudong x5-4-135 and the seismic data test to study the influence of logging and seismic wavelet on Synthetic Seismograms quality. It is proved that extracting the determined wavelet calibration by well log is better than Ricker and statistical wavelet.

well-seismic integration calibration; logging data; reflection coefficient; wavelet

2016-05-17

李海明(1967-)，男，河北灤縣人，副教授，主要從事石油工程專業(yè)教學和非常規(guī)油氣資源研究,郵箱lihaiming@cdpc.edu.cn。

P618.13

B

1008-9446(2016)06-0001-04