殷 文

（1.新疆克拉瑪依職業(yè)技術(shù)學(xué)院，新疆 克拉瑪依834000；2.中國石化 勝利油田有限公司 博士后工作站，山東 東營257002；3.山東省東營職業(yè)學(xué)院，山東 東營257091）

目前，隱蔽油氣藏勘探已占油田勘探的80%以上，勘探的難度也越來越大，如何提高儲層預(yù)測的精度、尋找剩余儲量，是今后相當(dāng)長的時(shí)間內(nèi)勘探開發(fā)所需要解決的重要問題。提高地震儲層預(yù)測精度的主要途徑之一就是開展地震地層的精細(xì)解釋［1－3］，因此在油氣勘探開發(fā)中，精細(xì)地層旋回劃分與對比研究必將成為地震地質(zhì)解釋人員必不可少的研究項(xiàng)目。

合成記錄將地球物理信息和鉆井地質(zhì)信息聯(lián)系了起來，做好合成記錄標(biāo)定工作是地質(zhì)評價(jià)中非常重要的環(huán)節(jié)。但在實(shí)際工作中，往往會遇到合成記錄標(biāo)定效果不理想、時(shí)深關(guān)系不夠準(zhǔn)確的問題，除了與合成記錄制作和層位標(biāo)定過程簡單化以及實(shí)際資料因素有關(guān)外，也與傳統(tǒng)的標(biāo)定手段單一有關(guān)［4－7］。

母巖風(fēng)化產(chǎn)物以及其他來源的沉積物在運(yùn)移、沉積過程中會按照顆粒的礦物成分、相對密度、形狀和大小在地表依次沉積下來，稱之為地表沉積分異作用。普斯托瓦洛夫認(rèn)為沉積物會首先按照顆粒大小和相對密度發(fā)生分異，因此在沉積地層中表現(xiàn)出層序特征。相應(yīng)地在測井與地震資料上也表現(xiàn)出相應(yīng)的特征。時(shí)頻譜在分析沉積旋回特征研究中體現(xiàn)出一定的優(yōu)越性，因此本文利用測井資料與地震資料的時(shí)頻譜研究層位對應(yīng)關(guān)系，從而校正時(shí)深關(guān)系。

1 井震層序分析理論依據(jù)

1.1 基本原理

層序體內(nèi)的物質(zhì)成分與其構(gòu)造特征是密切相關(guān)的，根據(jù)層序體內(nèi)沉積相、巖性變化和沉積環(huán)境的特點(diǎn)，其物性特征的分布通常會引起與構(gòu)造相關(guān)聯(lián)的橫向和垂向的變化［8，9］。結(jié)合層序地層學(xué)、地質(zhì)和地震勘探方法，建立地震、測井層序分析方法。井震資料是各級層序體的綜合體，各層序體由相應(yīng)各級別沉積間歇面所分割，沉積層序的層理性及其沉積旋回體現(xiàn)了它們的普遍埋藏形式。

作為地震和測井旋回體，層序體的頻率化與相對應(yīng)反射系數(shù)的變化是相關(guān)聯(lián)的。時(shí)頻分析技術(shù)是把實(shí)際地質(zhì)體與井震資料的頻率特性相結(jié)合，利用井震資料的頻率特征在空間上的變化規(guī)律可以反映實(shí)際地質(zhì)體的結(jié)構(gòu)關(guān)系［10－12］。根據(jù)時(shí)頻分析結(jié)果，可以獲得有關(guān)層序體的沉積旋回性，劃分出沉積地層的旋回及內(nèi)部精細(xì)結(jié)構(gòu)和物質(zhì)成分，能夠?qū)崿F(xiàn)在同一級別研究目標(biāo)上的井震資料綜合解釋（圖1）。因此，利用井震資料的時(shí)頻譜研究層位對應(yīng)關(guān)系，基于一個(gè)參考時(shí)深關(guān)系，把測井資料轉(zhuǎn)換到時(shí)間域，對多種曲線（井旁地震道、測井曲線及合成地震記錄等）進(jìn)行時(shí)頻分析，利用時(shí)頻譜以尋找可對比特征，用于約束精細(xì)時(shí)深標(biāo)定，從而校正時(shí)深關(guān)系（圖2），為合成記錄層位標(biāo)定和時(shí)深關(guān)系校正提供有效的技術(shù)支持。

圖1 井震旋回對比分析解釋Fig.1 The interpretation，comparison and analysis of seismic and well cycles

1.2 時(shí)頻技術(shù)

1.2.1 三角型濾波器

通過傅里葉變換可以將時(shí)域的地震記錄轉(zhuǎn)換成對應(yīng)的頻譜。但常見頻譜計(jì)算方法有兩點(diǎn)不足之處：一是計(jì)算量較大，二是計(jì)算精度受采樣長度的影響。計(jì)算精度越高，計(jì)算量也就隨之越大；因此，可以采用快速傅里葉變換法來提高計(jì)算效率。時(shí)頻分析的效果取決于濾波器的特征參數(shù)。

圖2 井震時(shí)頻譜約束時(shí)深關(guān)系校正流程Fig.2 Flow of time-depth relationship constrained by seismic and well time-frequency sepctrum

濾波器的特征參數(shù)的不同會影響到時(shí)頻分析結(jié)果的好壞。其特征包括：①濾波器的頻率響應(yīng)旁極值低平、極大值突出；②為了能夠確保濾波器的輸出穩(wěn)定，需要保證頻帶寬度足夠；③采用零相位濾波器，確保地震記錄的時(shí)間特征，即使得時(shí)移量為零；④要求濾波器的左、右截頻對數(shù)陡度一定，也就是說濾波器具有相似的濾波掃描響應(yīng)。

綜合以上特征要求，可以采用2倍頻程零相位的三角型濾波器。掃描濾波的步長Δf可由薄層反射系數(shù)的周期來決定。濾波器的遞進(jìn)優(yōu)勢頻率公式為

其中：k為頻率遞歸梯度，由濾波器的個(gè)數(shù)和主頻決定；fn為第n個(gè)濾波器的優(yōu)勢頻率。

通常情況下，k取1.05或1.2，假設(shè)第一個(gè)濾波器的主頻是20Hz，那么第二個(gè)濾波器的主頻則是21Hz或24Hz（圖3）。由此可見，每次改變?yōu)V波器后，相鄰濾波器之間有2/3的部分是覆蓋的，擴(kuò)充了1/3的新信息。

圖3 濾波器遞推算法示意圖Fig.3 The diagram of recursion algorithm of filters

1.2.2 小波變換

小波變換具有良好的頻域局部性和時(shí)域局部性，是理想的時(shí)頻分析工具。小波函數(shù)還具備多樣性的特點(diǎn)，在對地震和測井資料進(jìn)行小波計(jì)算分析時(shí)，分析小波可以采用與子波波形相近的［13－18］。

Morlet小波是常見的復(fù)值小波，能提取被分析的信號的幅值和相位或時(shí)間信息，小波系數(shù)的選取一般應(yīng)滿足歸一化條件‖φ‖2＝1。尺度為σ的 Morlet小波φσ，0（t）的傅里葉變換是

該尺度小波φσ，0（ω）的支撐區(qū)域基本是ω＞0的整個(gè)頻率域，因?yàn)楫?dāng)ω＜0而ω0＞5時(shí)，，因此φ（ω）值趨于零，可忽略不計(jì)。φσ，0（ω）的支撐區(qū)中心在波形寬度δ＝1/σ，并隨σ的減小而向外擴(kuò)展；而小波φσ，τ（t）的中心在τ處，以的方式隨著σ增大向外擴(kuò)展。

在使用 Morlet小波f（t）時(shí)，需明確小波變換的相位和模的平方曲線，即和曲線，之后方可分析變換的實(shí)部表示的過程。

小波變換應(yīng)當(dāng)是在頻率域內(nèi)而不是在時(shí)間域內(nèi)進(jìn)行分析?？赏ㄟ^以下方法來確定尺度范圍：如果f（t）信號的采樣時(shí)間是Δt，則由香農(nóng)采樣定理可得，采樣頻率fs應(yīng)當(dāng)是信號f（t）所包含的最高頻率的2倍，因?yàn)閒s＝1/Δt，fs≥2（f0/σmin），σmin是分析中的最小尺度，由此可得

通過確定小波φσ，τ（t）的擴(kuò)展寬度δ可以得出最大尺度σmax。鑒于｜φσ，τ｜的幅值在3σφσ，τ處會下降到自身的99.9%，故可增加3σφσ，τ≤（tmax－tmin）/2條件，即Morlet小波支撐區(qū)限定的取值范圍如下式

1.2.3 S型變換

S型變換是一種可逆的局部時(shí)頻變換方法。S型變換引入了頻率與寬度成反比例變化的高斯窗，分辨率與頻率相關(guān)［19－22］。信號x（t）的S型變換定義為

其中，f是尺度因子，τ參數(shù)是用來控制高斯窗在時(shí)間軸的位置，w（τ－t，f）為高斯窗。在S型變換中，高斯窗的寬度和高度是隨著頻率而變化的，進(jìn)而克服了短時(shí)傅里葉變換窗口寬度和高度固定的缺陷。由于信號x（t）的S型變換與其傅里葉變換X（f）之間存在以下關(guān)系

由上式可以得到S型變換的離散形式

當(dāng)n＝0時(shí)，

當(dāng)n≠0時(shí)，

S型變換可以由快速傅里葉變換實(shí)現(xiàn)快速計(jì)算，經(jīng)變換后得到二維復(fù)時(shí)頻矩陣，其中列對應(yīng)于采樣點(diǎn)，行對應(yīng)于頻率，對矩陣中的各元素求?？傻玫较鄳?yīng)的模時(shí)頻矩陣［23－25］。

2 時(shí)深關(guān)系校正

由于測井資料中存在漂移，在井震標(biāo)定過程中地震層位是固定不變的。測井層位相對于地震層位進(jìn)行測井曲線的拉伸、壓縮與平移，因此需要根據(jù)圖像映射原理對測井曲線層位進(jìn)行校正。圖4為圖像映射示意圖，在圖中采用插值方法對原始圖像進(jìn)行放大。

圖5為井震標(biāo)定層位標(biāo)定示意圖，其中T為地震資料中層位對應(yīng)的時(shí)間，t為測井曲線層位對應(yīng)的時(shí)間，t′為經(jīng)過校正之后的測井曲線上層位對應(yīng)的時(shí)間。測井曲線校正公式如下

經(jīng)過校正之后測井上的層位與地震層位相對應(yīng)，從而為時(shí)深關(guān)系校正提供了理論依據(jù)。相應(yīng)的校正前后的時(shí)深關(guān)系如圖6所示。

3 實(shí)際資料應(yīng)用

本文采用某地區(qū)的實(shí)際資料對該方法和軟件應(yīng)用效果進(jìn)行測試驗(yàn)證。

地震數(shù)據(jù)：疊前時(shí)間偏移數(shù)據(jù)；測井?dāng)?shù)據(jù)：53井、25 井 數(shù) 據(jù)；測 線 范 圍 InLine：2949-4101 CrossLine 1658-2040。Syn Rec列為合成記錄，Seimic列為地震剖面，Maker列為井標(biāo)記，Spectrum（Well）列為井曲線時(shí)頻譜，Spectrum（Seismic）列為井旁道時(shí)頻譜。

圖4 圖像映射示意圖Fig.4 The diagram of projected image

圖5 井震層位標(biāo)定示意圖Fig.5 The calibration diagram of seismic and well layers

圖6 時(shí)深關(guān)系校正前后示意圖Fig.6 The diagram before and after the time-depth relationship correction

有些測井與地震層位的關(guān)系無法直接從原始曲線上識別，這就可以借助于時(shí)頻分析方法將地震數(shù)據(jù)和測井?dāng)?shù)據(jù)轉(zhuǎn)化為多尺度時(shí)頻譜剖面。圖7為25井處時(shí)頻譜約束的精細(xì)時(shí)深關(guān)系校正前后對比圖，圖8為53井處時(shí)頻譜約束的精細(xì)時(shí)深關(guān)系校正前后對比圖。經(jīng)過標(biāo)定測井與地震層位以及對井震時(shí)頻譜沉積旋回性的對比分析，進(jìn)而實(shí)現(xiàn)對時(shí)深關(guān)系的標(biāo)定和校正，最終為提高時(shí)深關(guān)系準(zhǔn)確度提供了更為有效的技術(shù)支持。

圖7 采用小波變換的井震聯(lián)合時(shí)頻譜時(shí)深關(guān)系校正實(shí)例Fig.7 The instances of time-depth relationship correction united seismic and well time-frequency spectrum using wavelet transform

圖8 采用S變換的井震聯(lián)合時(shí)頻譜時(shí)深關(guān)系校正實(shí)例Fig.8 The instances of time-depth relationship correction united seismic and well time-frequency spectrum using S transform

4 結(jié)論

常規(guī)合成記錄標(biāo)定往往效果不理想，進(jìn)而影響到時(shí)深關(guān)系的準(zhǔn)確性。地層的沉積模式在測井和地震資料中都表現(xiàn)出沉積旋回特征，測井資料與地震資料也存在對應(yīng)的層序關(guān)系。而時(shí)頻譜在分析沉積旋回特征研究中表現(xiàn)出優(yōu)越性，選擇合適的時(shí)頻分析工具，根據(jù)井震時(shí)頻譜層位對應(yīng)關(guān)系，以地震資料為標(biāo)準(zhǔn)道進(jìn)行井曲線的調(diào)整，從而確立了新的時(shí)深關(guān)系，為時(shí)深關(guān)系校正提供了新的思路和技術(shù)支持。

［1］吳國忱，康仁華，印興耀.三維時(shí)頻分析方法在地震層序分析中的應(yīng)用［J］.石油大學(xué)學(xué)報(bào)：自然科學(xué)版，2000，24（1）：81－84.Wu G C，Kang R H，Yin X Y.Application of 3D time-frequency analysis method to seismic sequence analysis［J］.Journal of the University of Petroleum，China，2000，24（1）：81－84.（In Chinese）

［2］武建新.地震層序時(shí)頻分析系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)［D］.成都：電子科技大學(xué)檔案館，2011.Wu J X.The Designer and Application of Time-Frequency Analysis System on Seismic Sequence［D］.Chengdu：The Archive of University of Electronic Science and Technology of China，2011.（In Chinese）

［3］周靜，殷文.井震聯(lián)合高分辨率層序地層分析方法及軟件實(shí)現(xiàn)［J］.山東科技大學(xué)學(xué)報(bào)：自然科學(xué)版，2012，31（3）：40－46.Zhou J，Yin W.The method and application of high resolution sequence stratigraphy analysis with wellseismic unite［J］.Journal of Shandong University of Science and Technology（Natural Science），2012，31（3）：40－46.（In Chinese）

［4］Lu S，McMechan G A.Elastic impedance inversion of multichannel seismic data from unconsolidated sediments containing gas hydrate and free gas［J］.Geophysics，2004，69（1）：164–179.

［5］Quakenbush M，Shang B，Tuttle C.Poisson impedance［J］.The Leading Edge，2006，25（2）：128－126.

［6］劉傳虎，劉福貴，李衛(wèi)忠.時(shí)頻分析方法及在儲層預(yù)測中的應(yīng)用［J］.石油地球物理勘探，1996，31（1）：11－20.Liu C H，Liu F G，Li W Z.The method of timefrequency analysis and the application in reservoir prediction［J］.Oil Geophysical Prospecting，1996，31（1）：11－20.（In Chinese）

［7］Hart B，Chen M A.Understanding seismic attributes through forwarding modeling［J］.The Leading Edge，2004，23（9）：834－841

［8］侯樹杰.時(shí)頻分析方法對比及S變換在地震數(shù)據(jù)處理中的應(yīng)用［J］.油氣地質(zhì)與采收率，2011，18（3）：44－46.Hou S J.Method comparison of time-frequency analysis and application of S transform in seismic data processing［J］.Petroleum Geology and Recovery Efficiency，2011，18（3）：44－46.（In Chinese）

［9］張春賀，劉雪軍，何蘭芳，等.基于時(shí)頻電磁法的富有機(jī)質(zhì)頁巖層系勘探研究［J］.地球物理學(xué)報(bào)，2013，56（9）：3173－3183.Zhang C H，Liu X J，He L F，etal.A study of exploration organic rich shales using time-frequency electromagnetic method［J］.Chinese Journal of Geophysics，2013，56（9）：3173－3183.（In Chinese）

［10］龐銳，劉百紅，孫成龍.時(shí)頻分析技術(shù)在地震勘探中的應(yīng)用綜述［J］.巖性油氣藏，2013，25（3）：92－96.Pang R，Liu B H，Sun C L.Review on timefrequency analysis technique and its application in seismic exploration［J］.Lithologic Reservoirs，2013，25（3）：92－96.（In Chinese）

［11］李小梅，余娟麗.時(shí)頻分析技術(shù)在層序旋回劃分中的應(yīng)用［J］.石油與天然氣地質(zhì)，2008，29（6）：273－278.Li X M，Yu J L.Application of time-frequency analysis to division of cyclical sequence［J］.Oil ＆Gas Geology，2008，29（6）：273－278.（In Chinese）

［12］Philippe S，Guy D.Seismic sequence analysis and attribute extraction using quadratic time-frequency representations［J］.Geophysics，2001，66（6）：1947－1959.

［13］唐曉初.小波分析及其應(yīng)用［M］.重慶：重慶大學(xué)出版社，2006：115－135.Tang X C.The Analysis and Application of Wavelet［M］.Chongqing：Chongqing University Press，2006：115－135.（In Chinese）

［14］王艷忠，操應(yīng)長，遠(yuǎn)光輝.小波分析在深水砂礫巖和泥頁巖地層層序劃分中的應(yīng)用［J］.天然氣地球科學(xué)，2012，23（2）：251－258.Wang Y Z，Cao Y C，Yuan G H.Application of wavelet analysis in sequence division of deepwater glutenite and shale layer［J］.Natural Gas Geoscience，2012，23（2）：251－258.（In Chinese）

［15］任金鋒，廖遠(yuǎn)濤，孫鳴，等.基于小波變換的高精度層序地層定量劃分研究及其應(yīng)用［J］.地球物理學(xué)進(jìn)展，2013，28（5）：2651－2657.Ren J F，Liao Y，Sun M，etal.A method for quantitative division of sequence stratigraphy with highresolution based on wavelet transform and its application［J］.Process in Geophysics，2013，28（5）：2651－2657.（In Chinese）

［16］段春節(jié)，艾敬旭，葉增爐.基于小波變換的相干技術(shù)在斷層解釋中的應(yīng)用［J］.地球物理學(xué)進(jìn)展，2013，28（1）：434－437.Duan C J，Ai J X，Ye Z L.The application of wavelet transform coherency algorithm in seismic data interpretation［J］.Preogress in Geophysics，2013，28（1）：434－437.（In Chinese）

［17］周靜毅，印興耀，吳國忱.小波變換與地震奇異性屬性［J］.物探化探計(jì)算技術(shù)，2006，28（2）：117－120.Zhou J Y，Yin X Y，Wu G C.Wavelet transform and seismic singularity attribute［J］.Computing Techniques for Geophysical and Geochemical Exploration，2006，28（2）：117－120.（In Chinese）

［18］李偉，馬欽忠，宋志平.小波變換在地電場數(shù)據(jù)分析中的應(yīng)用［J］.地震學(xué)報(bào)，2013，35（1）：26－35.Li W，Ma Q Z，Song Z P.Application of wavelet transform to the analysis of geo-electric field data［J］.Acta Seismologica Sinica，2013，35（1）：26－35.（In Chinese）

［19］張春寧，陳紅坤，黃絹，等.基于S變換和決策樹算法的電能質(zhì)量擾動識別［J］.武漢大學(xué)學(xué)報(bào)：工學(xué)版，2010，43（6）：800－804.Zhang C N，Chen H K，Huang J，etal.Electricity quality disturbance identification method based on stock well-transform and decision tree［J］.Engineering Journal of Wuhan University，2010，43（6）：800－804.（In Chinese）

［20］武國寧，曹思遠(yuǎn)，馬寧，等.S變換的時(shí)頻分析特性及其改進(jìn)［J］.地球物理學(xué)進(jìn)展，2011，26（5）：1661－1667.Wu G N，Cao S Y，Ma N，etal.The time frequency property of S transform and its generalization［J］.Progress in Geophysics，2011，26（5）：1661－1667.（In Chinese）

［21］王長江，楊培杰，羅紅梅，等.基于廣義S變換的時(shí)變分頻技術(shù)［J］.石油物探，2013，52（5）：489－494.Wand C J，Yang P J，Luo H M，etal.The variable frequency technology based on generalized S transform［J］.Geophysical Prospecting for Petroleum，2013，52（5）：489－494.（In Chinese）

［22］魏學(xué)強(qiáng)，李輝峰，楊超.基于廣義S變換的沉積旋回分析方法研究［J］.西安石油大學(xué)學(xué)報(bào)：自然科學(xué)版，2013，28（4）：35－40.Wei X Q，Li H F，Yang C.Study on sedimentary cycle analysis method based on generalized S transform［J］.Journal of Xi’an Shiyou University（Natural Science Edition），2013，28（4）：35－40.（In Chinese）

［23］高靜懷，陳文超，李幼銘，等.廣義S變換與薄互層地震響應(yīng)分析［J］.地球物理學(xué)報(bào)，2003，46（4）：526－536.Gao J H，Chen W C，Li Y M，etal.Generalized S transform and seismic response analysis of thin interbeds［J］.Chinese Journal of Geophysics，2003，46（4）：526－536.（In Chinese）

［24］趙軍龍，蒙靈飛，張斌成，等.S變換在常規(guī)測井信息處理中的應(yīng)用綜述［J］.地球物理學(xué)進(jìn)展，2011，24（6）：1554－1559.Zhao J L，Meng L F，Zhang B C，etal.A review of S transform applications in conventional information processing of logging［J］.Progress in Geophysics，2011，24（6）：1554－1559.（In Chinese）

［25］王寶江，張永常，王大興，等.廣義S變換用于含煤地層薄砂體預(yù)測研究［J］.地球物理學(xué)進(jìn)展，2012，27（3）：1219－1225.Wang B J，Zhang Y C，Wang D X，etal.Predicting and detecting of the thin sand bodies by general transform in coal-bearing strata［J］.Progress in Geophysics，2012，27（3）：1219－1225.（In Chinese）