李煜偉, 呂宗剛, 彭海潤(rùn), 王 琳, 李 軍

(1.油氣藏地質(zhì)及開(kāi)發(fā)工程國(guó)家重點(diǎn)實(shí)驗(yàn)室(成都理工大學(xué))，成都 610059；2.中國(guó)石油西南油氣田公司 蜀南氣礦，四川 瀘州 642450)

雙相介質(zhì)理論檢測(cè)頁(yè)巖裂縫及其油氣

李煜偉1, 呂宗剛2, 彭海潤(rùn)2, 王 琳2, 李 軍2

(1.油氣藏地質(zhì)及開(kāi)發(fā)工程國(guó)家重點(diǎn)實(shí)驗(yàn)室(成都理工大學(xué))，成都 610059；2.中國(guó)石油西南油氣田公司 蜀南氣礦，四川 瀘州 642450)

嘗試用雙相介質(zhì)BISQ理論，在疊后地震數(shù)據(jù)上計(jì)算以頻率為自變量的瞬時(shí)頻散譜和瞬時(shí)吸收譜，用瞬時(shí)譜的頻散和衰減信息檢測(cè)頁(yè)巖裂縫及油氣。用迭代優(yōu)化方法設(shè)計(jì)頻寬為2 Hz的25個(gè)分頻算子，用分頻算子對(duì)地震數(shù)據(jù)分頻，用復(fù)數(shù)道分析技術(shù)對(duì)所有分頻道進(jìn)行“三瞬”計(jì)算，集合成瞬時(shí)振幅譜、瞬時(shí)頻率譜、瞬時(shí)相位譜，用這3個(gè)瞬時(shí)譜計(jì)算瞬時(shí)頻散譜和瞬時(shí)吸收譜。計(jì)算結(jié)果表明：在瞬時(shí)吸收譜的低頻段可見(jiàn)到強(qiáng)吸收被弱吸收封閉的異常特征，在瞬時(shí)頻散譜上伴有頻散突變異常特征，兩種異常共同構(gòu)成檢測(cè)頁(yè)巖裂縫及油氣的標(biāo)志。用瞬時(shí)頻散譜和瞬時(shí)吸收譜的衰減和頻散的異常特征對(duì)8口已知井頁(yè)巖儲(chǔ)層的檢測(cè)結(jié)果與鉆井吻合?；陔p相介質(zhì)BISQ理論的瞬時(shí)頻散譜和瞬時(shí)吸收譜可預(yù)測(cè)油氣分布范圍，提供鉆探目標(biāo)。

頁(yè)巖裂縫;雙相介質(zhì);瞬時(shí)頻散譜;瞬時(shí)吸收譜;頻率;油氣預(yù)測(cè)

影響頁(yè)巖氣分布并決定其是否具有工業(yè)勘探開(kāi)發(fā)價(jià)值的重要因素有：泥頁(yè)巖厚度與體積、有機(jī)質(zhì)類型與豐度、熱歷史與有機(jī)質(zhì)成熟度、孔隙度與滲透率、斷裂與裂縫以及構(gòu)造運(yùn)動(dòng)與現(xiàn)今埋藏深度等[1]。

裂縫在頁(yè)巖氣藏中的作用表現(xiàn)在：(1)有利于后期對(duì)頁(yè)巖層的改造工作；(2)如果裂縫規(guī)模過(guò)大，特別是穿層斷層發(fā)育時(shí)，可能導(dǎo)致頁(yè)巖氣藏中的游離氣和吸附氣逸散而減少；(3)裂縫是油氣運(yùn)移和儲(chǔ)集的主要空間[2]。裂縫性油氣藏在油氣生產(chǎn)中占有非常重要的地位。

裂縫性油氣藏由于孔隙度低、非均質(zhì)性強(qiáng)且裂縫分布復(fù)雜，使得裂縫預(yù)測(cè)成為公認(rèn)的難題。國(guó)內(nèi)外對(duì)這一難題進(jìn)行了深入研究，在20世紀(jì)70年代，國(guó)外有裂隙性油氣藏的勘探方法研究專著發(fā)表。20世紀(jì)90年代，地震三維三分量、地震三維多方位等資料的利用取得一些效果。近幾年來(lái)，由單相各向異性研究轉(zhuǎn)到了雙相各向異性的研究，雙相各向異性更符合客觀實(shí)際。而中國(guó)裂隙性油氣藏的研究主要是在20世紀(jì)80年代以來(lái)才有所進(jìn)展，對(duì)以上領(lǐng)域研究差距也逐漸縮小[2]。

因波場(chǎng)差異小造成頁(yè)巖裂縫識(shí)別困難，需要用敏感參數(shù)。速度頻散和吸收衰減是雙相介質(zhì)BISQ理論的核心信息，也是識(shí)別頁(yè)巖裂縫及油氣的敏感信息。用敏感的核心信息可望檢測(cè)到波場(chǎng)特征差異小的頁(yè)巖裂縫及油氣。因此，有必要開(kāi)發(fā)基于雙相介質(zhì)理論檢測(cè)頁(yè)巖裂縫及油氣的新方法。

1 理論依據(jù)

地震波在每一類型的介質(zhì)中傳播都有其獨(dú)特的規(guī)律。雙相介質(zhì)BISQ理論指出流體的存在和運(yùn)動(dòng)是彈性波速度頻散和吸收衰減的主要原因，頻散和衰減揭示流體的存在和運(yùn)動(dòng)[3]，流體的存在和運(yùn)動(dòng)寓含裂縫的存在。受地震分辨率限制無(wú)法檢測(cè)到頁(yè)巖納米級(jí)儲(chǔ)集空間的游離態(tài)和吸附態(tài)頁(yè)巖氣，但可檢測(cè)含油氣的頁(yè)巖裂縫。雙相介質(zhì)BISQ理論是檢測(cè)頁(yè)巖裂縫及油氣的理論依據(jù)。

1.1 雙相介質(zhì)BISQ理論

固體和流體組成的雙相介質(zhì)能更全面地揭示地震波的衰減機(jī)制。Dvorikn和Nur于1993年建立了流體沿縱向上的Biot流動(dòng)和沿橫向上的噴射流動(dòng)引起的頻散和衰減參數(shù)與彈性參數(shù)、儲(chǔ)層參數(shù)、流體參數(shù)和頻率參數(shù)為函數(shù)的BISQ理論[3]：

(1)

式中：VP為速度頻散;αd為吸收衰減; J0為零階貝塞爾函數(shù); J1為一階貝塞爾函數(shù);ρs為固體密度;φ為孔隙度;ρf為流體密度;η為流體黏滯度;k為滲透率;S為飽和度;M為干噪骨架的單軸應(yīng)變模量;Fsq為固-流耦合系統(tǒng)的綜合壓縮性;α為有效應(yīng)力孔隙彈性系數(shù);R為特征噴射流動(dòng)長(zhǎng)度;F為固-流耦合系統(tǒng)的壓縮性;K為干噪骨架體積模量;ω為圓頻率[3,4]。

(1)式表明：速度頻散和吸收衰減與彈性參數(shù)、儲(chǔ)層參數(shù)、流體參數(shù)和頻率參數(shù)的函數(shù)關(guān)系，以參數(shù)形式表述其內(nèi)涵。然而，方程參數(shù)太多，又以不可測(cè)量者居多，目前無(wú)法實(shí)現(xiàn)用地震資料在反演中以參數(shù)形式進(jìn)行描述和應(yīng)用。根據(jù)用頻率作自變量表示信號(hào)為頻譜的定義，可以用地震偏移數(shù)據(jù)計(jì)算出以頻率為自變量的速度頻散譜和吸收衰減譜，以譜的形式表述其內(nèi)涵?？紤]到地震信號(hào)是時(shí)變信號(hào)和開(kāi)時(shí)窗計(jì)算引起的平均效應(yīng)，把速度頻散譜和吸收衰減譜定義在地震道的采樣點(diǎn)上，而具有瞬時(shí)性。以瞬時(shí)速度頻散譜和瞬時(shí)吸收衰減譜的形式表述其內(nèi)涵。瞬時(shí)頻散譜和瞬時(shí)吸收譜上的頻散和衰減信息揭示流體的存在和運(yùn)動(dòng)，流體的存在和運(yùn)動(dòng)寓含裂縫的存在。

瞬時(shí)頻散譜和瞬時(shí)吸收譜的計(jì)算，首先用小波多分辨分析導(dǎo)出的帶通濾波把地震信號(hào)分解成具有瞬時(shí)特征的單分量信號(hào)，用Hilbert變換計(jì)算單分量信號(hào)的瞬時(shí)振幅、瞬時(shí)頻率和瞬時(shí)相位。用“三瞬”計(jì)算瞬時(shí)頻散譜和瞬時(shí)吸收譜。

1.2 小波多分辨分析

小波多分辨分析可導(dǎo)出帶通濾波概念。小波多分辨分析用時(shí)間-尺度聯(lián)合函數(shù)用“變焦距”的方法分析時(shí)變信號(hào)。地震道f(t)定義在信號(hào)空間v：f(t)∈v，尺度函數(shù)φ(t)對(duì)f(t)進(jìn)行分解[5]

(2)

信號(hào)空間被分解為鑲嵌空間〔Vi,(i=1,2,…)〕。a空間Vi+1包含在它的前一級(jí)a空間之中t

V0?V1?……Vi?Vi+1…

(3)

地震道在鑲嵌空間中被分解后，在Vi的部分fi(t)比在Vi+1中的fi+1(t)含有更多的信息Δfi+1

Δfi+1(t)=fi(t)-fi+1(t)=f(t)Δφi+1(t)

(4)

(2)式中φ(a-jt)相當(dāng)于低通濾波，(4)式中Δφi+1相當(dāng)于帶通濾波。小波多分辨分析導(dǎo)出的帶通濾波把地震信號(hào)分解成單分量信號(hào)而具有瞬時(shí)特征，是計(jì)算瞬時(shí)頻散譜和瞬時(shí)吸收譜的基礎(chǔ)。

2 瞬時(shí)頻散譜和瞬時(shí)吸收譜計(jì)算

地震道的有效頻帶為50 Hz，設(shè)計(jì)25個(gè)分頻帶通濾波器，其頻帶寬度為2 Hz。用李慶忠院士的迭代優(yōu)化方法[6]計(jì)算25個(gè)帶通濾波算子就可把一個(gè)地震道分成25個(gè)分頻道，合成后仍為50 Hz，分頻道近似為單分量信號(hào)具有瞬時(shí)特征。用復(fù)數(shù)道分析技術(shù)的Hilbert變換對(duì)25個(gè)分頻道分別進(jìn)行瞬時(shí)振幅、瞬時(shí)頻率、瞬時(shí)相位計(jì)算。用瞬時(shí)振幅作因變量，瞬時(shí)頻率作自變量用點(diǎn)斜式可計(jì)算出瞬時(shí)吸收α；當(dāng)速度是頻率的函數(shù)v(f)與相位函數(shù)θ的關(guān)系是[4]：v(f)=ω(f)/θ(f)，用瞬時(shí)相位譜作因變量，瞬時(shí)頻率作自變量可計(jì)算瞬時(shí)頻散。瞬時(shí)道集可合成頻帶寬度為50 Hz，頻率采樣率為2 Hz的5個(gè)瞬時(shí)譜：瞬時(shí)振幅譜A(f)，瞬時(shí)頻率譜ω(f)，瞬時(shí)相位譜θ(f)，瞬時(shí)吸收譜((f)，瞬時(shí)頻散譜v(f)。計(jì)算公式如下

(5)

(6)

(7)

(8)

(9)

式中：x(t)、y(t)分別為解析信號(hào)的實(shí)部、虛部；“′”表示求導(dǎo)運(yùn)算；i為分辨道序號(hào)[4, 6]。

瞬時(shí)頻散譜和瞬時(shí)吸收譜的分辨率高，信噪比高，精度高。地震道的有效頻帶為50 Hz，瞬時(shí)譜頻帶寬度也為50 Hz，保真度高。瞬時(shí)頻散譜和瞬時(shí)吸收譜整個(gè)計(jì)算過(guò)程屬于數(shù)值計(jì)算，有解存在且穩(wěn)定，是適定的。突出特點(diǎn)是：2個(gè)譜定義在地震道的釆樣點(diǎn)上，具有瞬時(shí)特征，信息量大，多解性小。

流體的存在和運(yùn)動(dòng)寓意著裂縫的存在，因此，用瞬時(shí)頻散譜和瞬時(shí)吸收譜的頻散和衰減信息可預(yù)測(cè)頁(yè)巖裂縫及油氣。

3 瞬時(shí)頻散譜和瞬時(shí)吸收譜檢測(cè)頁(yè)巖裂縫及油氣

瞬時(shí)頻散譜和瞬時(shí)吸收譜的吸收和頻散信息的解釋是本方法的靈魂。頻散和衰減受頻率影響很大；相速度(速度頻散)急劇變化帶，與滲透率相關(guān)，衰減峰值的位置依賴于滲透率的大小；流體黏滯系數(shù)增大，衰減峰值向低頻移動(dòng)；衰減具有各向異性，衰減和滲透率各向異性密切相關(guān)，P波的最強(qiáng)衰減發(fā)生在最大滲透率方向[2]。地層存在流體，在地層壓力下不逃逸，則封閉。因此，在瞬時(shí)吸收譜的低頻段(<30 Hz)強(qiáng)吸收封閉，并且低頻段伴有頻散強(qiáng)異?；蝾l散弱異常，為解釋裂縫及油氣模式。

3.1 用已知井建立預(yù)測(cè)裂縫及油氣模式

對(duì)川中、川西南、川南3個(gè)地區(qū)的7個(gè)構(gòu)造和1個(gè)區(qū)塊，5個(gè)層段(須四、須二、嘉四1-嘉三、嘉二1-嘉一、陽(yáng)新統(tǒng))。其中氣井15口，氣水同產(chǎn)井2口，水井4口，干井13口，氣枯竭井2口[4]，共36口已知井的常規(guī)儲(chǔ)層進(jìn)行檢驗(yàn)分析，得出瞬時(shí)吸收譜的低頻段強(qiáng)吸收封閉模式適用于不同的構(gòu)造和儲(chǔ)層，是固定異常模式。低頻段頻散強(qiáng)異常(頻散弱異常)模式適用于一個(gè)地區(qū)的同一儲(chǔ)層。對(duì)吸收和頻散2種異常共同擁有的CDP范圍作裂縫及油氣解釋。

對(duì)鉆遇頁(yè)巖非常規(guī)儲(chǔ)層裂縫的8口已知井進(jìn)行檢驗(yàn)分析，低頻段強(qiáng)吸收封閉模式適用于不同的構(gòu)造和儲(chǔ)層，是固定異常模式。川西南威遠(yuǎn)構(gòu)造寒武系筇竹寺組頁(yè)巖和川南地區(qū)3個(gè)局部構(gòu)造上的志留系龍馬溪組黑色頁(yè)巖的檢測(cè)，都表現(xiàn)為低頻段強(qiáng)吸收封閉模式。但速度頻散異常有差異，川西南威遠(yuǎn)構(gòu)造筇竹寺組頁(yè)巖為低頻段頻散強(qiáng)異常模式；川南地區(qū)龍馬溪組黑色頁(yè)巖為低頻段頻散弱異常模式。這種差異可能與所處的構(gòu)造部位、巖性、巖相、地層組合等因素有關(guān)。不同地區(qū)和不同的儲(chǔ)層，其速度有差異，速度頻散也有差異，這種多解性，屬于正常的地質(zhì)響應(yīng)。根據(jù)鉆井的有關(guān)信息，這種多解性是可排除的。

地震不能分辨頁(yè)巖非常規(guī)儲(chǔ)層納米級(jí)儲(chǔ)集空間的流體，可以分辨裂縫及油氣。因此，低頻段強(qiáng)吸收封閉和低頻段頻散強(qiáng)異常(頻散弱異常)模式構(gòu)成頁(yè)巖裂縫及油氣預(yù)測(cè)模式。對(duì)吸收和頻散2種異常共同擁有的CDP范圍作裂縫及油氣解釋。

吸收封閉：瞬時(shí)吸收譜上低頻段(<30 Hz)弱吸收封閉強(qiáng)吸收。封閉與被封閉的吸收差值，稱為封閉強(qiáng)度；封閉縱向“厚度”稱為封閉頻寬。吸收和頻散2種異常共同擁有的CDP個(gè)數(shù)稱為封閉范圍。頻散異常：瞬時(shí)頻散譜上相同頻段橫向上強(qiáng)弱的變化。圖1-圖8彩色圖由上而下分別為沿目的層的層位線在1/3周期范圍內(nèi)處理的瞬時(shí)吸收譜和瞬時(shí)頻散譜，縱坐標(biāo)是頻率坐標(biāo)分別為0～50 Hz，橫坐標(biāo)為CDP序號(hào)，彩色色標(biāo)表示頻散和衰減相對(duì)大小的百分?jǐn)?shù)，用黑色點(diǎn)線標(biāo)出異常，黑色直線為分析瞬時(shí)頻散譜異常的參考線。

3.2 川西南威遠(yuǎn)構(gòu)造檢測(cè)頁(yè)巖裂縫及油氣

川西南威遠(yuǎn)構(gòu)造上的W76、W112、W3井鉆遇筇竹寺組頁(yè)巖見(jiàn)到井涌、氣浸，W201井在筇竹寺組頁(yè)巖日產(chǎn)氣1.08×104m3，表明4口井鉆遇頁(yè)巖裂縫。低頻段強(qiáng)吸收封閉是固定裂縫模式，是由低頻段頻散強(qiáng)異常伴隨它。按照模式對(duì)2種異常共同擁有的CDP范圍作裂縫解釋與鉆井吻合。頻散強(qiáng)異常模式在這個(gè)地區(qū)的筇竹寺組頁(yè)巖是穩(wěn)定的(圖1-圖4，表1)。

首先按照裂縫解釋模式作裂縫解釋，然后作油氣及其產(chǎn)能解釋。相速度(速度頻散)急劇變化帶，與滲透率相關(guān)；衰減峰值的位置依賴于滲透率的大??；流體黏滯系數(shù)增大，衰減峰值向低頻移動(dòng)[3]。由油礦地質(zhì)理論可知，滲透率反映儲(chǔ)層巖石的滲流能力，與儲(chǔ)層產(chǎn)能直接相關(guān)[7]。表1中的封閉強(qiáng)度、封閉頻寬越大，向低頻移動(dòng)量越大，則產(chǎn)量越大，對(duì)威遠(yuǎn)構(gòu)造的頁(yè)巖裂縫及油氣檢測(cè)獲得令人滿意的效果。

圖1 W201井裂縫及油氣檢測(cè)

圖2 W76井裂縫及油氣檢測(cè)

圖3 W112井裂縫及油氣檢測(cè)

圖4 W3井裂縫及油氣檢測(cè)

圖5 Y63井裂縫及油氣檢測(cè)

圖6 L32井裂縫及油氣檢測(cè)

圖7 Y深2井裂縫及油氣檢測(cè)

圖8 P1井裂縫及油氣檢測(cè)

表1 威遠(yuǎn)構(gòu)造鉆井情況與吸收特征

3.3 川南地區(qū)檢測(cè)頁(yè)巖裂縫

川南地區(qū)鉆遇龍馬溪組黑色頁(yè)巖層段時(shí)，普遍見(jiàn)氣測(cè)異常、井涌、氣浸和井噴。其中Y63井下志留統(tǒng)龍馬溪組黑色頁(yè)巖段，酸化后獲日產(chǎn)氣3 500 m3；L32井龍馬溪組黑色碳質(zhì)頁(yè)巖段日產(chǎn)氣1 948 m3；Y深2井井涌、氣浸；P1井在龍馬溪組為氣測(cè)異常。按照低頻段強(qiáng)吸收封閉和頻散弱異常模式作裂縫解釋與鉆井吻合。頻散弱異常模式在川南地區(qū)的3個(gè)局部構(gòu)造上對(duì)龍馬溪組頁(yè)巖的應(yīng)用是穩(wěn)定的，說(shuō)眀檢測(cè)方法是穩(wěn)定的。表2中的封閉強(qiáng)度、封閉頻寬越大，向低頻移動(dòng)量越大，產(chǎn)量越大，也與鉆井吻合。對(duì)于P1井(圖8)，距離區(qū)域性大斷層很近，大斷層斷開(kāi)陽(yáng)頂，可能使頁(yè)巖中的游離氣和吸附氣逃逸而剩余較少，因此吸收不封閉，無(wú)異常存在，與鉆井吻合。對(duì)川南地區(qū)的頁(yè)巖裂縫及油氣檢測(cè)獲得令人滿意的效果(圖5～圖8，表2)。

4 瞬時(shí)頻散譜和瞬時(shí)吸收譜潛在能力討論

瞬時(shí)頻散譜和瞬時(shí)吸收譜的吸收衰減和速度頻散揭示流體的存在和運(yùn)動(dòng)，影響流體流動(dòng)的因素主要有流動(dòng)方向、流體黏度、裂縫連通性等。這些流動(dòng)因素又影響吸收衰減和速度頻散的變化，具體表現(xiàn)在：裂縫滲透率與流動(dòng)方向有關(guān)，滲透率取決于流體流動(dòng)方向，在不同方向上，具有不同的滲透率值。吸收衰減和滲透率各向異性密切相關(guān)，P波的最強(qiáng)衰減發(fā)生在最大滲透率方向，吸收衰減峰值的位置依賴于滲透率的大小。速度頻散急劇變化帶，與滲透率相關(guān)[3,7]。石油、天然氣、地層水的黏度差別很大，大部分原油的黏度都高于5 mPa·s，水在20℃時(shí)的黏度為1.005×10-3mPa·s，天然氣的黏度在20℃時(shí)為12×10-3mPa·s。流體黏度差異影響流動(dòng)能力，表現(xiàn)在流體黏滯系數(shù)增大，衰減峰值向低頻移動(dòng)[3,7]；裂縫具有孔滲性，裂縫的孔隙度與裂縫的密度成正比。裂縫的連通性對(duì)滲透率的主導(dǎo)控制作用很強(qiáng)，如果裂縫的孔隙度很小，裂縫連通性好則滲透率很高。若基巖孔隙度比裂縫的孔隙度大，當(dāng)連通性相對(duì)較差，基巖滲透率較低。通常裂縫滲透率很高，基巖滲透率相對(duì)較低[3,7]。流體的流動(dòng)方向、流體黏度、裂縫連通性影響流體流動(dòng)，流體流動(dòng)又影響吸收衰減和速度頻散的變化，因此，瞬時(shí)頻散譜和瞬時(shí)吸收譜吸收衰減和速度頻散的變化作為檢測(cè)裂縫及油氣的依據(jù)。

由油礦地質(zhì)學(xué)可知：在一定壓差下，巖石允許流體通過(guò)的能力稱之為滲透性，孔喉大小的宏觀表現(xiàn)為滲透率，滲透率反映儲(chǔ)層巖石的滲流能力，與儲(chǔ)層產(chǎn)能直接相關(guān)[7]。在檢測(cè)到裂縫的基礎(chǔ)上，用封閉強(qiáng)吸收的封閉強(qiáng)度、封閉頻寬，低頻方向移動(dòng)的頻率量，可進(jìn)一步作油氣產(chǎn)能的解釋。

a.有產(chǎn)能的W201、Y63、L32井封閉強(qiáng)度大、封閉頻寬大，向低頻移動(dòng)量越大，產(chǎn)量大。

b.無(wú)產(chǎn)能的W76、W112、W3、Y深2井封閉強(qiáng)度小、封閉頻寬小。

表2 川南地區(qū)鉆井情況與吸收特征

c.用三維地震資料圈出油氣分布范圍為油藏描述提供參考信息。

依靠水平井多級(jí)壓裂、微地震等技術(shù)“人造滲透率”增大儲(chǔ)集體中油氣流動(dòng)能力，才能產(chǎn)出具經(jīng)濟(jì)價(jià)值的天然氣。雙相介質(zhì)理論研究成果指出，最強(qiáng)吸收衰減發(fā)生在最大滲透率方向[3]，瞬時(shí)吸收譜和瞬時(shí)頻散譜預(yù)測(cè)最大滲透率方向有潛在的積極意義。

5 結(jié) 論

雙相介質(zhì)BISQ波動(dòng)方程參數(shù)太多，又以不可測(cè)量者居多，目前無(wú)法實(shí)現(xiàn)其地震反演的應(yīng)用。用地震數(shù)據(jù)計(jì)算的瞬時(shí)吸收譜和瞬時(shí)頻散譜實(shí)現(xiàn)了雙相介質(zhì)BISQ理論在檢測(cè)頁(yè)巖裂縫及油氣中的應(yīng)用，獲得了令人滿意的效果。

a.瞬時(shí)吸收譜和瞬時(shí)頻散譜信息量擴(kuò)大了25倍。信息類型多，信息量大，多解性小，獲得地質(zhì)信息更客觀。

b.瞬時(shí)吸收譜和瞬時(shí)頻散譜可建立裂縫預(yù)測(cè)模式，檢測(cè)裂縫密度相對(duì)小、宏觀規(guī)模小的頁(yè)巖裂縫及其油氣，也可檢測(cè)常規(guī)儲(chǔ)層裂縫及油氣。

c.吸收封閉強(qiáng)度、封閉頻寬信息有預(yù)測(cè)油氣及其產(chǎn)能的潛在能力。在應(yīng)用中既單一，又直接，異常特征突出，異常邊界清晰，效果明顯。

感謝我的導(dǎo)師劉樹根教授，感謝李智武、孫瑋副教授與李硯、岳宏高級(jí)工程師，他們對(duì)作者的指導(dǎo)和幫助及對(duì)本文提出寶貴的修改意見(jiàn)。

[1] 劉樹根,曾祥亮,黃文明,等.四川盆地頁(yè)巖氣藏和連續(xù)型-非連續(xù)型氣藏基本特征[J].成都理工大學(xué)學(xué)報(bào)：自然科學(xué)版,2009,36(6):578-590.

Liu S G, Zeng X L, Huang W M,etal. Basic characteristics of shale and continuous-discontinuous transition gas reservoirs in Sichuan Basin, China[J]. Journal of Chengdu University of Technology (Science & Technology Edition), 2009, 36(6): 578-590. (In Chinese)

[2] 劉成齋.泥巖裂縫油藏地球物理特征與檢測(cè)方法[D].成都:成都理工大學(xué)檔案館,2004:74-85.

Liu C Z. Shale Oil-gas Reservoir Geophysical Characteristics and Detection Methods[D]. The Archive of Chengdu University of Technology, 2004: 74-85. (In Chinese)

[3] 楊頂輝,陳小宏.含流體多孔介質(zhì)的BISQ模型[J].石油地球理勘探,2001,36(2):146-159.

Yang D H, Chen X H. BISQ model for fluid-filled, porous medium[J]. Oil Geophysical Prospecting, 2001, 36(2): 146-159. (In Chinese)

[4] 李煜偉,李硯,呂宗剛,等.雙相介質(zhì)BISQ地震油氣預(yù)測(cè)技術(shù)[J].天然氣工業(yè),2009,29(12):29-31.

Li Y W, Li Y, Lyu Z G,etal. Seismic hydrocarbon detection based on the unified Biot and squirt-flow mechanism (BISQ) model of two-phase medium[J]. Natural Gas Industry, 2009, 29(12): 29-31. (In Chinese)

[5] 曹聿明.地震道的全息特征[J].石油地球物理勘探,1998,33(6):716-722.

Cao L M. Holographic characteristic of seismic trace[J]. Oil Geophysical Prospecting, 1998, 33(6): 716-722. (In Chinese)

[6] 李慶忠.地震信號(hào)的內(nèi)插與噪聲剔除[J].地球物理學(xué)報(bào)，1987,30(5):517-519.

Li Q Z. Interpolating seismic signal and noise removed[J]. Chinese Journal of Geophysics, 1987, 30(5): 517-519.

[7] 吳元燕,吳勝和,蔡正旗.油礦地質(zhì)學(xué)[M].北京:石油工業(yè)出版社,2009:237-242.

Wu Y Y, Wu S H, Cai Z Q. Oil Field Geology[M]. Beijing: Petroleum Industry Press, 2009: 237-242. (In Chinese)

[8] 鄧?yán)^新,王尚旭.儲(chǔ)層砂巖縱波衰減的實(shí)驗(yàn)結(jié)果分析[J].石油物探,2008,47(4):346-351.

Deng J X, Wang S X. Experimental analysis of P-wave attenuation in sandstone of reservoir with ultrasonic frequency[J]. Geophysical Prospecting for Petroleum, 2008, 47(4): 346-351. (In Chinese)

[9] 席道瑛,張程遠(yuǎn),劉小燕.儲(chǔ)層砂巖的各向異性[J].石油地球物理勘探,2001,36(2):187-192.

Xi D Y, Zhang C Y, Liu X Y. Anisotropy in reservoir sandstone[J]. Oil Geophysical Prospecting, 2001, 36(2): 187-192.(In Chinese)

[10] 崔志文,王克協(xié)，曹正良,等.多孔介質(zhì)BISQ模型中的慢縱波[J].地球物理學(xué)報(bào)，2004,53(9): 3083-3087.

Cui Z W, Wang K X, Cao Z L,etal. Slow waves propagation in BISQ poroelastic model[J]. Acta Physica Sinica, 2004, 53(9): 3083-3087.(In Chinese).

[11] Biot M A. Theory of propagation of elastic waves in a fluid saturated porous solid, low frequency range and higher-frequency range[J]. Journal of the Acoustical Society of America, 1956, 28(2): 168-191.

[12] Biot M A. Mechanics of deformation and acoustic propagation in porous media[J]. Journal of Applied Physics, 1962, 33(4): 1482-1498.

[13] Dvorkin J, Nur A. Dynamic poroelasticity: a unified model with the Squirt flow and the Biot Mechanism[J]. Geophysics, 1993, 58(3): 524-533.

[14] Dvorkin J, Nolen-Hoeksema R, Nur A. The Squirt flow mechanism: macroscopic description[J]. Geo-physics, 1994, 59(3): 428-438.

[15] Johnston D H, Toksoz M N. Ultrasonic P and S wave attenuation in dry and saturated rocks underpressure[J]. Journal of Geophysical Research, 1980, 85(B4): 925-936.

[16] Mavko G, Nur A. Wave attenuation in partially saturated rocks[J]. Geophysics, 1979, 44:161-178.

[17] Pahner I D, Traviolia M L. Attenuation by squirt flow in undersaturated gas sands[J]. Geophysics, 1980, 45: 1781-1792.

[18] Murphy W F. Acoustic relaxation in sedimentary rocks: Dependence on grain contacts and fluid saturation[J]. Geophysics, 1986, 51: 757-766.

[19] Miksis M J. Effect of contact line movement on the dissipation of wave in partially saturated rocks[J]. J Geophys Res, 1986, 93: 6624-6634.

[20] Wang Z, Nur A. Dispersion analysis of acoustic velocities in rocks[J]. J Acoust Soc Am, 1980, 87: 2384-2395.

[21] Akbar N. Relating P-wave attenuation to permeability[J]. Geophysics, 1993, 58: 20-29.

[22] Dvorkin J, Nolen-Hoeksema R, Nur A. The squirt-flow mechanism: Macroscopic description[J]. Geophysics, 1994, 59: 428-438.

Detectionofshalefracturesandoil/gasbybiphasicmediumtheory

LI Yu-wei1, LYU Zong-gang2, PENG Hai-run2, WANG Lin2, LI Jun2

1.StateKeyLaboratoryofOilandGasReservoirGeologyandExploitation,ChengduUniversityofTechnology,Chengdu610059,China;2.SouthSichuanGasField,PetroChinaSouthwestOil-GasFieldCompany,Luzhou642450,China

Using the biphasic medium BISQ theory, this paper tries to calculate the instantaneous frequency dispersion spectrum and the absorption spectrum taking frequency as the independent variable based on the post-stacked seismic data and to detect the shale fractures and the oil and gas via spectrum frequency dispersion and attenuation information on the instantaneous spectrum. By using the iterative optimization method, 25 frequency division operators (whose bandwidth is 2 Hz) are designed, the frequency division is done to the seismic record using the frequency division operators, and tri-instant processing is conducted for all frequency division traces using the complex trace analytic method, which are then combined into instantaneous amplitude spectrum, instantaneous frequency spectrum and instantaneous phase spectrum. Then, calculation is done with these three instantaneous spectra to get the instantaneous spectral dispersion and absorption spectra. The calculation result has the following indications. First of all, in the low frequency of the transient absorption spectra, abnormal characteristics can be seen when the strong absorption is sealed by the weak absorption. Second, the instantaneous frequency dispersion spectrum is associated with the dispersion mutation. The above two anomalies jointly become the key signs in terms of detecting the shale fractures and the oil and gas. The test result of the shale reservoirs of 8 known wells by using the dispersion and attenuation anomalies of the instantaneous frequency dispersion spectrum and the absorption spectrum accords with the exploratory drilling result. So the instantaneous frequency dispersion spectrum and the absorption spectrum based on the biphasic medium BISQ theory can predict gas and oil distribution and locate drilling targets.

shale fracture; biphasic medium; instantaneous frequency dispersion spectrum; instantaneous frequency absorption spectrum; frequency; petroleum prediction

10.3969/j.issn.1671-9727.2013.06.04

1671-9727(2013)06-0648-10

P631.4; TE132

A

2013-06-03

李煜偉(1981-),男,碩士研究生,研究方向:石油地質(zhì)學(xué), E-mail:95554409@qq.com。