王永生 （中石油青海油田分公司勘探事業(yè)部，甘肅 敦煌736202）

蔡涵鵬 （中國石油川慶鉆探工程有限公司地球物理勘探公司，四川 成都610213

成都理工大學(xué)油氣藏地質(zhì)及開發(fā)工程國家重點(diǎn)實(shí)驗(yàn)室，四川 成都610059）

何佳 （中石油青海油田分公司鉆采工藝研究院，甘肅 敦煌736202）

朱波 （中石油青海油田分公司勘探事業(yè)部，甘肅 敦煌736202）

劉維俊 （中石油青海油田分公司采油二廠，青海 茫崖816500）

程濤 （中石油青海油田分公司鉆采工藝研究院，甘肅 敦煌736202）

王文卓 （中石油青海油田分公司勘探開發(fā)研究院，甘肅 敦煌736202）

1 基本理論

當(dāng)縱波在孔隙介質(zhì)中傳播時(shí)，巖石的部分孔隙收縮，而部分孔隙擴(kuò)張。如果巖石孔隙飽含流體，則壓縮部分的孔隙壓力要比擴(kuò)張部分的孔隙壓力大。當(dāng)孔隙彼此連通時(shí)，高孔壓孔隙中的流體將向低孔壓孔隙中流動(dòng)。如果流體是黏滯性的，那么孔隙中流體的相對流動(dòng)將促使地震波能量衰減，即地震波能量被吸收，同時(shí)出現(xiàn)速度頻散。

2 模型分析

為了分析頻散對反射系數(shù)和相位角影響的特征，根據(jù)上覆巖層和下伏儲(chǔ)層砂巖波阻抗差異大小可以將油氣勘探中遇到的砂巖儲(chǔ)層分為3類：①下伏儲(chǔ)層砂巖波阻抗大于上覆頁巖波阻抗，該類儲(chǔ)層用于描述固結(jié)良好的砂巖儲(chǔ)層；②下伏砂巖儲(chǔ)層波阻抗與上覆頁巖波阻抗接近，用于描述中等壓實(shí)和固結(jié)的砂巖儲(chǔ)層；③下伏砂巖儲(chǔ)層波阻抗小于上覆頁巖波阻抗，該類儲(chǔ)層代表未固結(jié)或者固結(jié)差的砂巖儲(chǔ)層。對于每一個(gè)試驗(yàn)?zāi)Ｐ?，?chǔ)層由相同巖石骨架的薄層構(gòu)成，且飽含鹽水層和飽含氣層交替排列。盡管層狀模型可能不符合實(shí)際情況，但是層狀模型較好地描述了與White周期層狀非飽和狀態(tài)模型[1]有密切聯(lián)系的吸收。

砂巖顆粒的體積模量為38GPa，密度為2.65g/cm3。根據(jù)3類砂巖儲(chǔ)層的特征，設(shè)計(jì)3個(gè)對應(yīng)模型，其3類砂巖儲(chǔ)層巖石骨架和孔隙流體的物性參數(shù)見表1[2]。表1中，Kdry和Gdry為干燥巖石的體積模量和剪切模量； 為孔隙度；K為滲透率；cm為曲度參數(shù)；a為孔隙尺度因子；K、ρ、η分別為孔隙流體的體積模量、密度和黏滯系數(shù)；下標(biāo)g和w代表孔隙中流體成分為氣和水。由于假設(shè)上覆頁巖無吸收和頻散，因此上覆頁巖參數(shù)相對簡單，其參數(shù)見表2。表2中，vp和vs分別為縱波速度和橫波速度。

表1 巖石骨架及孔隙流體的物性參數(shù)

當(dāng)飽含氣層厚度d1和飽含水層厚度d2均等于1m時(shí)，利用Biot理論模型預(yù)測多孔彈性介質(zhì)吸收和頻散的運(yùn)動(dòng)方程[3～7]可計(jì)算獲得孔隙介質(zhì)的相速度和品質(zhì)因子倒數(shù) （1/Q）隨頻率變化的曲線（如圖1、2）。圖1、2表明，模型3未固結(jié)砂巖具有最強(qiáng)的速度頻散和吸收，其次是中等壓實(shí)程度的砂巖 （模型2），模型1砂巖的速度頻散和吸收最弱。頻率影響相速度和吸收的頻段主要集中在低頻段。3類模型砂巖的相速度均表現(xiàn)為隨著頻率的增加而增加，頻率越低相速度的變化越大，高頻段相速度隨頻率變化緩慢，逐漸趨于一個(gè)極大值。3類模型砂巖均表現(xiàn)為低頻對吸收影響較大，并在某個(gè)低頻處出現(xiàn)最強(qiáng)吸收 （圖2中箭頭標(biāo)記處），但隨著頻率逐漸增加，吸收逐漸減弱且趨于穩(wěn)定。

圖1 相速度隨頻率變化的曲線

圖2 1/Q隨頻率變化的曲線

表2 上覆頁巖物性參數(shù)

計(jì)算獲得的縱波法向入射反射系數(shù)的振幅絕對值和相位角隨頻率變化曲線見圖3、4所示。模型1中，儲(chǔ)層砂巖具有壓實(shí)程度高、孔隙度和滲透率低的特點(diǎn)，巖石的頻散和吸收相對于中等壓實(shí)和未壓實(shí)砂巖的頻散和吸收小。在零頻率時(shí)，砂巖波阻抗值高于上覆蓋層波阻抗值，而相速度隨著頻率逐漸增加而增加，使得波阻抗差異隨著頻率增加而增加，因此反射系數(shù)振幅隨著頻率增加而增加。頻率從1Hz至100Hz，反射系數(shù)振幅增加約0.008。由于頻散微弱，致使相位角隨頻率變化相應(yīng)微弱，且在所有的頻率上均為較小的正值。

圖3 反射系數(shù)振幅絕對值隨頻率變化的曲線

圖4 相位角隨頻率變化的曲線

對于模型2，砂巖具有中等壓實(shí)、中等孔隙度和較好滲透率的特征，其巖石的吸收和速度頻散比模型1中的大。頻率從1Hz至100Hz，相速度增加約145m/s。在零頻率時(shí)，砂巖波阻抗略低于上覆巖層波阻抗。由于速度頻散，相速度隨著頻率增加而增加，致使砂巖波阻抗隨頻率的增加而逐漸增加。因此從低頻到高頻相位發(fā)生反轉(zhuǎn)，從負(fù)轉(zhuǎn)向正 （圖4中箭頭標(biāo)記處）。

對于模型3，在零頻率時(shí)，砂巖波阻抗比上覆巖層的波阻抗小很多，相速度頻散強(qiáng)，頻率從1Hz至100Hz，相速度增加約210m/s。隨著頻率的增加相速度增加，致使砂巖的波阻抗與上覆泥巖波阻抗差異減小，因此反射系數(shù)振幅減少。頻率從1Hz至100Hz，反射系數(shù)振幅絕對值減少約0.08；所有與頻率相關(guān)的相位角均為較小的負(fù)值。

厚度d1＝d2交替組成的儲(chǔ)層模型相當(dāng)于儲(chǔ)層含水飽和度為50%的情況。為了分析不同飽和狀態(tài)時(shí)介質(zhì)的吸收對反射系數(shù)和相位角影響的效果，可以利用飽含氣和飽含水層厚度的不同來定性地分析吸收對反射系數(shù)和相位角影響的效果[23]。當(dāng)d1＝0.5m，d2＝1m時(shí)，計(jì)算獲得的3類儲(chǔ)層相速度、1/Q、反射系數(shù)振幅絕對值和相位角隨頻率的變化曲線如圖5所示。當(dāng)d1＝1.5m，d2＝0.8m時(shí)，計(jì)算獲得的3類儲(chǔ)層相速度、吸收、法向反射系數(shù)振幅絕對值和相位角隨頻率變化的曲線如圖6所示。圖5、6表明，d1和d2差異（即含水飽和度的不同）對3個(gè)模型的相速度、1/Q和反射系數(shù)振幅絕對值隨頻率變化趨勢的影響較小，僅在各個(gè)參數(shù)趨于穩(wěn)定的頻率發(fā)生變化，隨著含水飽和度的增加頻率趨于穩(wěn)定增大，但在低頻段對各個(gè)參數(shù)的變化率影響較大；d1和d2的變化對模型1和模型3的相位角影響不大，而對模型2的相位角有較大影響，即對于模型2代表的儲(chǔ)層類型，由于含水飽和度的不同，可能導(dǎo)致在地震波的有效頻帶內(nèi)，觀察不到相位的反轉(zhuǎn)。

圖5 相速度、1/Q、反射系數(shù)振幅絕對值和相位角隨頻率變化的曲線 （d1＝0.5m，d2＝1m）

圖6 相速度、1/Q、反射系數(shù)振幅絕對值和相位角隨頻率變化的曲線 （d1＝1.5m，d2＝0.8m）

由于3個(gè)模型中水的黏滯系數(shù)均是氣的黏滯系數(shù)的100倍以上，致使d2增加時(shí)（含水飽和度增加），相速度、1/Q和反射系數(shù)振幅絕對值隨頻率變化的變化率增加，即儲(chǔ)層孔隙中包含黏滯系數(shù)越大的流體，其相速度、1/Q和反射系數(shù)振幅絕對值在低頻段形成的異常越強(qiáng)。模型2中，飽含流體的黏滯系數(shù)越大，發(fā)生相位反轉(zhuǎn)的頻率向低頻段移動(dòng)，即在地震波有效頻帶內(nèi)，相位發(fā)生反轉(zhuǎn)的可能性越大。

由于孔隙中流體的相對流動(dòng)才能促使地震波能量衰減，導(dǎo)致速度頻散。Michael等[8]利用巖石物理試驗(yàn)證實(shí)了，甚至在地震頻帶內(nèi)，低的流體流動(dòng)可能產(chǎn)生強(qiáng)烈的速度頻散?？刂屏黧w流動(dòng)特征的參數(shù)主要是巖石的滲透率和孔隙中流體的黏滯性。Michael等[8]利用巖石物理試驗(yàn)證實(shí)了流體的黏滯系數(shù)隨著溫度的上升急劇降低，呈指數(shù)降低的趨勢。相對于巖石的滲透率對流體流動(dòng)的影響，流體的黏滯性對流體的影響較小。根據(jù)Biot理論模型運(yùn)動(dòng)方程[3～7]和Trapeznikova建立的縱波法向入射到兩種吸收介質(zhì)之間界面的反射系數(shù)和相位角的公式[9]及表1的參數(shù)，筆者探討了巖石滲透率的變化對相速度、1/Q、反射系數(shù)振幅和相位角隨頻率的影響。圖7（a）、（b）和 （c）分別為3個(gè)模型的各個(gè)參數(shù)隨滲透率變化時(shí)，以頻率為變量的特征曲線 （箭頭指示了滲透率增大的方向）。圖7表明，隨著巖石滲透率減小，即孔隙內(nèi)流體的流動(dòng)能力降低，孔隙內(nèi)壓力達(dá)到平衡所需的時(shí)間增加，致使發(fā)生頻散和高吸收的頻率更低，可能超出地震波的有效頻帶。因此巖石的滲透率控制著速度頻散和強(qiáng)吸收是否在地震波的有效頻帶內(nèi)發(fā)生，其分析結(jié)果與Michael等[8]利用巖石物理試驗(yàn)獲得結(jié)果是一致的。

圖7 滲透率對各個(gè)參數(shù)的影響

3 儲(chǔ)層流體識別準(zhǔn)則和步驟

根據(jù)縱波法向入射時(shí)，反射系數(shù)振幅和相位角隨頻率變化的特征，可以將振幅和相位角隨頻率變化曲線劃分為3類：低頻暗點(diǎn)消失儲(chǔ)層 （模型1），相位可能反轉(zhuǎn)儲(chǔ)層 （模型2），低頻亮點(diǎn)儲(chǔ)層 （模型3）。

3.1 識別準(zhǔn)則

根據(jù)吸收和速度頻散引起的3類儲(chǔ)層頂部反射系數(shù)振幅和相位角的變化規(guī)律，建立基于自激自收地震數(shù)據(jù)識別3類儲(chǔ)層流體的識別準(zhǔn)則。為了避免大的濾波作用、擴(kuò)散作用、儲(chǔ)層上覆地層吸收等因素對激發(fā)地震子波影響的非一致性，假設(shè)到達(dá)儲(chǔ)層頂部反射界面的子波特征相對穩(wěn)定。在地震頻帶范圍內(nèi)，3類儲(chǔ)層流體識別準(zhǔn)則如表3。

表3 3類儲(chǔ)層流體識別準(zhǔn)則

3.2 識別步驟

基于吸收和速度頻散特征的儲(chǔ)層流體識別的主要步驟為：

1）利用研究區(qū)內(nèi)鉆井和測井資料分析儲(chǔ)層類型。

2）利用測井或者巖石物理測試評價(jià)儲(chǔ)層的滲透率。

3）地震資料相對振幅保持處理，盡可能保持地震波信號的頻寬和頻譜中的低頻信息成分。

4）根據(jù)合成記錄標(biāo)定儲(chǔ)層在地震剖面上的位置，追蹤相應(yīng)層位，并提取儲(chǔ)層頂部反射界面的地震子波。

5）利用瞬時(shí)譜分解技術(shù)提取儲(chǔ)層頂部反射層位處與角頻率相關(guān)的瞬時(shí)振幅譜和瞬時(shí)相位譜。

6）利用子波振幅譜對提取的瞬時(shí)振幅譜進(jìn)行加權(quán)處理，獲得反射系數(shù)瞬時(shí)振幅譜。

7）利用識別準(zhǔn)則分析加權(quán)處理后的瞬時(shí)振幅譜或瞬時(shí)相位譜，或者兩者結(jié)合確定儲(chǔ)層流體的橫向展布。

4 實(shí)例分析

為了檢驗(yàn)利用吸收和速度頻散特征識別儲(chǔ)層流體的有效性，利用上述流體識別準(zhǔn)則及識別步驟對塔河油田S區(qū)塊的三維地震數(shù)據(jù)進(jìn)行了處理和分析。地震數(shù)據(jù)的采樣率為2ms，主頻率約35Hz。圖8是從三維數(shù)據(jù)體中抽取的過THN1井的垂直地震剖面，圖8中標(biāo)注了T64的位置，T64為中三疊統(tǒng)阿克庫勒組三段 （T2a3）對應(yīng)的地震反射波組，黑色曲線為THN1井的合成記錄。儲(chǔ)層位于T64的下方 （箭頭指示段），即T64層位為上覆泥巖與下伏含氣砂巖的反射界面。圖9為THN1鉆井的測井曲線圖，曲線包括自然電位（Usp）、密度（ρ）、縱波速度（vp）和計(jì)算的縱波阻抗（Zp）。該井中含油氣砂巖段儲(chǔ)層厚度約11m，波阻抗約為679kg/（s·cm2）。上覆泥巖層厚度約為45m，波阻抗約為810kg/（s·cm2）。根據(jù)波阻抗差異，該地區(qū)的儲(chǔ)層屬于前述的第3類儲(chǔ)層。根據(jù)測井和巖心薄片統(tǒng)計(jì)分析，該儲(chǔ)層的滲透率約為0.35mD。

圖8 過THN1井的垂直地震剖面

圖9 THN1井測井曲線

測井分析表明，該地區(qū)縱波阻抗能較好地區(qū)分砂巖和泥巖，圖10是沿T64層提取的縱波阻抗切片，相變明顯 （黑色實(shí)線所示），有利相帶內(nèi)砂巖分布范圍廣，但其流體分布特征雜亂，難以直接識別和解釋，而從利用上述原理生成的反射系數(shù)瞬時(shí)振幅譜數(shù)據(jù)體抽取的12Hz和8Hz沿T64層的單頻切片中（圖11），THN1井鉆遇的含氣砂體在沿T64層的反射系數(shù)低頻瞬時(shí)振幅切片中表現(xiàn)為強(qiáng)能量 （圖11中虛線框內(nèi)），且8Hz切片中的能量強(qiáng)于12Hz切片中的能量，與第3類儲(chǔ)層流體判別準(zhǔn)則相符，即反射系數(shù)振幅隨頻率的增大而減小。而在8Hz和12Hz切片中能量中西北方向上的大片高能量區(qū)位于相變線以外的區(qū)域，利用沉積相分析可以判斷該高異常區(qū)可能不是油氣的因素造成的?？梢娫诜瓷湎禂?shù)低頻瞬時(shí)振幅單頻切片中目的層的流體特征突出地展現(xiàn)出來了，清晰地展布了砂巖體中含氣分布范圍，其分布形態(tài)完整、邊緣清楚、橫向非均質(zhì)性亦清晰可見，但這些特征信息在原始切片和波阻抗切片中并不明顯。圖11中虛線框內(nèi)展布的氣藏已被多口鉆井的地質(zhì)資料所證實(shí)，并已確定是一個(gè)高產(chǎn)氣藏。

圖10 沿T64層縱波阻抗切片

圖11 沿T64層反射系數(shù)瞬時(shí)振幅譜單頻切片

5 結(jié) 語

White的周期層狀非飽和模型理論和Biot孔隙彈性理論表明，當(dāng)縱波從非頻散介質(zhì)法向入射至頻散介質(zhì)時(shí)，反射系數(shù)振幅和相位角會(huì)隨頻率發(fā)生規(guī)律性的變化。根據(jù)反射系數(shù)振幅和相位角隨頻率的變化特征，可以將含流體儲(chǔ)層分為3類。儲(chǔ)層頂部反射系數(shù)振幅和相位角隨頻率變化特征表明，在某種程度上，利用疊后地震數(shù)據(jù)識別流體是可行的，可以作為儲(chǔ)層中流體識別的一個(gè)依據(jù)，并證實(shí)地震資料處理中，合理地保留較低頻率成分信息的重要性。實(shí)例分析表明，巖性、儲(chǔ)層和地震資料處理技術(shù)應(yīng)用的不合理均可導(dǎo)致反射系數(shù)低頻強(qiáng)異常，因此在一定的條件下，利用儲(chǔ)層頂界面反射系數(shù)的異常特征識別儲(chǔ)層中的流體時(shí)，應(yīng)與沉積相、地震相、阻抗反演等技術(shù)相結(jié)合，多項(xiàng)技術(shù)研究結(jié)果相互印證，降低其多解性。此外，基于周期層狀非飽和模型理論計(jì)算反射系數(shù)振幅和相位角隨頻率的變化時(shí)，考慮了薄層、薄互層流體對其的影響，因此為薄層或者薄互層儲(chǔ)層流體識別提供了新的理論參考，具有重要的實(shí)際應(yīng)用價(jià)值。

[1]White J E.Computed seismic speeds and attenuation in rocks with partial gas saturation [J].Geophysics，1975，40 （2）：224～232.

[2]Ren H R，Goloshbin G，Hilterman F J.Poroelastic analysis of amplitude－versus－frequency variations [J].Geophysics，2009，74（6）：41～48.

[3]Carcione J M，Picotti S.P－wave seismic attenuation by slowwave diffusion：Effects of inhomogeneous rock properties [J].Geophysics，2006，71 （3）：1～8.

[4]Biot M A.Theory of propagation of elastic waves in a fluid－saturated porous solid：Part2—Higher frequency range [J].Journal of the Acoustical Society of America，1956，28 （2）：179～191.

[5]Schmidt H，Tango G.Efficient global matrix approach to the computation of synthetic seismograms[J].Geophysical Journal of the Royal Astronomical Society，1986，84 （2）：331～359.

[6]Stern M，Bedford A，Millwater H R.Wave reflection from a sediment layer with depth－dependent properties[J].Journal of the A－coustical Society of America，1985，77 （1）：1043～1063.

[7]Abramowitz M，Stegun A.Handbook of mathematical functions with formulas，graphs，and mathematical tables [M].Dover：Dover Publications，1972.

[8]Michael L B，De H H，Rronny H.Fluid mobility and frequency－dependent seismic velocity—Direct measurements [J].Geophysics，2006，71 （1）：1～9.

[9]Trapeznikova N A.Prognosis and interpretation of seismic wave dynamics[D].Russian：Nauka，1985.