陳占國(guó),陳 林,張衛(wèi)紅

(中國(guó)石油化工股份有限公司石油物探技術(shù)研究院,江蘇南京211103)

基于walkaround VSP的裂縫檢測(cè)方法及應(yīng)用

陳占國(guó),陳 林,張衛(wèi)紅

(中國(guó)石油化工股份有限公司石油物探技術(shù)研究院,江蘇南京211103)

研究油氣儲(chǔ)層裂縫發(fā)育情況是降低油氣勘探風(fēng)險(xiǎn),有效開(kāi)發(fā)油氣資源的重要任務(wù)之一。研究了利用walkaround VSP資料下行qP波的振幅、速度和極性方位各向異性特征分析裂縫參數(shù)的方法和流程。其中通過(guò)反褶積和偽透射系數(shù)( pseudo transmission coefficient,PTC)處理消除了上覆地層對(duì)振幅的影響;通過(guò)計(jì)算與各向異性有關(guān)的局部速度避開(kāi)炮點(diǎn)高程和偏移距及上覆地層對(duì)旅行時(shí)的影響;提出的相對(duì)極化角方法校正了偏振極性受炮點(diǎn)位置的影響;最后通過(guò)最小二乘橢圓擬合方法分析了qP波的振幅、速度和偏振極性隨方位的變化特征,統(tǒng)計(jì)了儲(chǔ)層的裂縫走向和發(fā)育情況。應(yīng)用結(jié)果表明,3種屬性分析的結(jié)果吻合度較高,且和橫波分裂方法的分析結(jié)果一致,從而證明該方法可靠性較高,是一種通過(guò)walkaround VSP預(yù)測(cè)裂縫的有效方法。

walkaround VSP;裂縫;各向異性;振幅;速度;極性

與裂縫方向和流體分布非均質(zhì)性有關(guān)的各向異性研究是地震勘探的重要內(nèi)容,已成為現(xiàn)今地球物理領(lǐng)域的研究熱點(diǎn)之一。縫洞型油氣藏是我國(guó)油氣勘探開(kāi)發(fā)中一種十分重要的油氣藏。在碳酸鹽巖和火山巖儲(chǔ)層中,油氣藏主要為縫洞型,即使在碎屑巖儲(chǔ)層中,也有相當(dāng)一部分油氣藏為裂縫型。因此識(shí)別裂縫型儲(chǔ)層對(duì)于油氣勘探開(kāi)發(fā),特別是碳酸鹽巖地區(qū)的油氣勘探開(kāi)發(fā)具有重要的現(xiàn)實(shí)意義[1]。

VSP裂縫儲(chǔ)層預(yù)測(cè)技術(shù)多以基于橫波分裂的檢測(cè)技術(shù)為主,但橫波受地表影響嚴(yán)重,信噪比和頻率較低,同時(shí)國(guó)內(nèi)外利用VSP檢測(cè)裂縫較多采用雙橫波震源激發(fā)的四分量技術(shù),這種橫波勘探方法激發(fā)條件苛刻,勘探成本較高,難以廣泛應(yīng)用[2-8]?；诜轿桓飨虍愋缘目v波方法在地面地震裂縫檢測(cè)技術(shù)中扮演著重要的角色。與其它技術(shù)比較,縱波更低的采集成本和相對(duì)成熟的處理技術(shù),使得利用縱波來(lái)進(jìn)行裂縫預(yù)測(cè)更具現(xiàn)實(shí)性[9-12]。walkaround VSP,其多方位采集的特點(diǎn),使地震波縱波速度、振幅等多種屬性攜帶了裂縫的方位各向異性信息;而井中接收的方式,使地震波受地面影響少,衰減小,信噪比較高;同時(shí)多分量多方位接收的方式能夠記錄與裂縫方位各向異性有關(guān)的縱波偏振信息。因此walkaround VSP縱波的這些特點(diǎn)對(duì)分析預(yù)測(cè)地層裂縫非常有利。但是國(guó)內(nèi)外通過(guò)walkaround VSP縱波信息研究裂縫各向異性的應(yīng)用實(shí)例較少,而且walkaround VSP觀測(cè)方式在國(guó)內(nèi)也少有采集。因此通過(guò)借鑒地面地震中縱波屬性分析裂縫的一些方法和手段,結(jié)合walkaround VSP自身特點(diǎn),研究利用其qP波信息分析井旁裂縫情況的方法和流程,對(duì)提高井旁裂縫預(yù)測(cè)的可靠性具有實(shí)際意義。

根據(jù)裂縫介質(zhì)中地震波的傳播理論[13],包含一組垂直定向排列的裂縫系統(tǒng)影響地震波的方式取決于裂縫的尺寸與占主導(dǎo)地位的地震波波長(zhǎng)的相對(duì)關(guān)系。如果裂縫尺寸比信號(hào)中的主波長(zhǎng)小得多,則裂縫導(dǎo)致介質(zhì)像一個(gè)有著與裂縫走向垂直對(duì)稱(chēng)軸的等效各向異性地層,因此qP波的振幅、速度和偏振極性的方位各向異性與裂縫的走向有關(guān)[14]。本文介紹了利用qP波振幅、速度和偏振極性隨方位的變化特征來(lái)綜合分析地層的裂縫發(fā)育情況的理論;結(jié)合walkaround VSP的采集方式的特點(diǎn),研究了利用walkaround VSP多種屬性檢測(cè)裂縫的方法,總結(jié)出一套實(shí)用的處理流程;最后通過(guò)一口井實(shí)測(cè)的walkaround VSP資料驗(yàn)證了方法技術(shù)的實(shí)用性。

1 qP波地震屬性與裂縫關(guān)系

假設(shè)裂縫只有一個(gè)主走向,則多方位觀測(cè)的VSP的qP波的振幅、速度和偏振極性等屬性可以用一個(gè)最小二乘橢圓來(lái)擬合,橢圓擬合的長(zhǎng)、短軸的走向與裂縫方位有關(guān),而長(zhǎng)軸A和短軸B的比值可以表示各向異性百分比程度(圖1)。

(1)

研究表明,對(duì)于簡(jiǎn)單介質(zhì),qP波沿裂縫走向和垂直裂縫走向傳播時(shí)的速度不同,從而導(dǎo)致qP波振幅響應(yīng)的變化[13]。這種qP波通過(guò)裂縫介質(zhì)表現(xiàn)出來(lái)的方位各向異性特征是進(jìn)行裂縫檢測(cè)的直接參數(shù)[15-16]。對(duì)于固定炮檢距,直達(dá)qP波通過(guò)裂縫介質(zhì)時(shí)的振幅響應(yīng)R與炮檢方向和裂縫走向之間的夾角φ(圖1)有如下關(guān)系[17-21]:

(2)

式中:C為與方位振幅有關(guān)的偏振因子(均勻介質(zhì)下透射強(qiáng)度);D為與方位振幅有關(guān)的調(diào)制因子；θ表示裂縫走向。顯而易見(jiàn),當(dāng)φ=0時(shí),振幅R最大;隨著φ的增大,R變小;當(dāng)φ=90°時(shí),R最小。因此,振幅橢圓擬合的長(zhǎng)軸對(duì)應(yīng)著裂縫的走向[9]。

圖1 橢圓擬合裂縫走向

與振幅與裂縫方位的關(guān)系一樣,速度v與炮檢方向和裂縫走向之間夾角φ也滿足公式(2)的形式[22],只不過(guò)R用v代替,C表示為與速度有關(guān)的偏振因子(方位速度平均值),D表示為與方位速度有關(guān)的調(diào)制因子。因此速度橢圓擬合的長(zhǎng)軸也對(duì)應(yīng)著裂縫的走向[23]。

最近的研究表明,沿著垂直井方向測(cè)量的qP波的極化角ψ與介質(zhì)的弱各向異性有關(guān)[24]:

(3)

式中:Q為垂直慢度;Φ為方位角;vP為垂向速度;Sn是與局部各向異性有關(guān)的參數(shù)。方程(3)表明可以通過(guò)測(cè)量極化角來(lái)確定介質(zhì)的各向異性參數(shù)。由于極化角ψ的余弦與局部速度vP成比例(一階近似),所以極化角的橢圓擬合的短軸對(duì)應(yīng)了裂縫的走向[14]。

2 裂縫參數(shù)提取方法

2.1 qP波波場(chǎng)預(yù)處理

在裂縫分析之前要對(duì)原始數(shù)據(jù)做一些預(yù)處理,以分離下行qP波場(chǎng)。這些處理包括：通過(guò)偏振分析方法對(duì)水平分量檢波器方位的旋轉(zhuǎn);對(duì)垂直分量和徑向分量的旋轉(zhuǎn)以及中值濾波等。由于要通過(guò)振幅來(lái)分析方位各向異性特征,因此不需要采用幾何擴(kuò)散補(bǔ)償?shù)仁侄螌?duì)振幅進(jìn)行處理。

2.2 qP波振幅屬性提取

為了降低隨機(jī)干擾對(duì)初至波能量的影響,對(duì)下行qP波場(chǎng),在初至附近的一個(gè)時(shí)窗里提取初至的均方根(RMS)振幅。一般情況下,尤其對(duì)于陸上walkaround VSP采集,受地形條件限制,不能采用等高程、等偏移距的規(guī)則排列進(jìn)行觀測(cè),同時(shí)也受地表激發(fā)條件及上覆地層等因素影響,qP波的RMS振幅與炮點(diǎn)方位表現(xiàn)出很強(qiáng)的相關(guān)性。這種非裂縫因素對(duì)振幅的影響往往遠(yuǎn)強(qiáng)于地層裂縫對(duì)方位振幅的影響,因此必須予以消除。

首先消除近地表激發(fā)條件對(duì)qP波振幅的影響。如圖2所示,不同方位的炮點(diǎn)處地表地質(zhì)條件不同,同時(shí)近地表也存在非均勻體,這些因素導(dǎo)致震源激發(fā)能量不同,子波的形態(tài)和能量也存在差異。為此設(shè)計(jì)一個(gè)零相位反濾波算子,通過(guò)反褶積將下行qP波能量相對(duì)均衡，同時(shí)也對(duì)子波整形。具體做法是選擇每一炮集的最淺層檢波點(diǎn)下行qP波記錄作為反褶積參考道,對(duì)該炮集qP波記錄進(jìn)行反褶積處理。這相當(dāng)于將每一炮的能量用最淺層檢波點(diǎn)歸一化,用最淺層檢波點(diǎn)對(duì)波形整形,以此消除最淺層檢波點(diǎn)以上地表部分及激發(fā)條件對(duì)qP波振幅的影響。

接下來(lái)消除地形和上覆地層各向異性對(duì)目標(biāo)層振幅的影響。在陸上勘探中,受地形限制,各炮點(diǎn)的高程及偏移距往往不一致,使地震波傳播距離不同,衰減也不一樣，同時(shí)上覆地層的各向異性影響也會(huì)傳遞到目標(biāo)層。在井中,相鄰檢波點(diǎn)垂直深度差別不大,同一炮兩個(gè)檢波點(diǎn)接收的qP波傳播路徑可以認(rèn)為是相同的(圖2)。因此,相鄰兩個(gè)檢波點(diǎn)之間地層的偽透射系數(shù)(pseudo transmission coefficient,PTC)可用目標(biāo)地層頂、底部檢波點(diǎn)的圴方根(RMS)振幅AT和BT的比值表示:

(4)

經(jīng)過(guò)計(jì)算偽透射系數(shù),相當(dāng)于將同一炮相鄰兩個(gè)檢波點(diǎn)之間對(duì)振幅影響相同的因素(如傳播距離、子波形態(tài)、目標(biāo)層上覆地層等)抵消掉,此時(shí)的PTC值就可以認(rèn)為只受一個(gè)因素影響,即目標(biāo)層處的裂縫或各向異性引起的qP波振幅變化。

圖2 上覆地層對(duì)qP波振幅的影響示意

2.3 qP波局部速度計(jì)算

變化較大的炮點(diǎn)高程及偏移距是造成qP波初至旅行時(shí)和炮點(diǎn)位置表現(xiàn)出很強(qiáng)相關(guān)性的主要原因,其次才是上覆地層各向異性。這些影響造成的旅行時(shí)時(shí)差較難校正。高志凌等[25]提出了一種沿射線方向計(jì)算非零偏VSP局部層速度的方法,這種層速度比沿垂直方向計(jì)算的VSP速度更接近偏移速度。本文將這種速度計(jì)算方法引入到walkaround VSP層速度的計(jì)算,并通過(guò)該局部層速度分析walkaround VSP裂縫各向異性,這種層速度比沿垂直方向計(jì)算的層速度對(duì)局部方位各向異性更為敏感。

如圖3所示,一般情況下,有偏VSP的偏移距遠(yuǎn)大于井中兩個(gè)相鄰檢波器的間距,我們可認(rèn)為相鄰兩個(gè)接收點(diǎn)的直達(dá)qP波傳輸路徑相似,兩接收點(diǎn)直達(dá)qP波傳輸路徑之差(Δs)為兩接收點(diǎn)垂直間距(ΔH)與直達(dá)qP波相對(duì)垂直方向入射角α(通過(guò)偏振分析可以得到)余弦的乘積,則兩個(gè)檢波點(diǎn)之間的局部速度可近似的表示為:

(5)

式中:ti和ti+1表示相鄰兩個(gè)接收點(diǎn)的初至?xí)r間。通過(guò)偏振信息和走時(shí)信息計(jì)算的層速度與常規(guī)的沿垂直方向計(jì)算的層速度相比,更能反映地層各向異性分析特征變化。

圖3 qP波局部速度計(jì)算示意

2.4 qP波相對(duì)極性計(jì)算

通過(guò)偏振極性分析裂縫是VSP等多分量地震的一大優(yōu)勢(shì)。在各向同性介質(zhì)中,P波的偏振方向與波的傳播方向一致;在各向異性介質(zhì)中,qP波偏振方向不再平行或垂直于波傳播方向(圖4)。由公式(3)也表明,接收點(diǎn)處質(zhì)點(diǎn)的運(yùn)動(dòng)引起的qP波偏振現(xiàn)象與地層的局部各向異性有關(guān)。通過(guò)三分量記錄的偏振分析可以獲得每個(gè)檢波點(diǎn)處的qP波相對(duì)井軸(z軸)的偏振角,這個(gè)偏振角除了與局部各向異性有關(guān)外,也與炮點(diǎn)位置有關(guān)。為了消除炮點(diǎn)位置的影響,提出了一個(gè)相對(duì)偏振極性來(lái)表示各向異性引起的偏振方向變化。如圖4,當(dāng)?shù)卣鸩ń?jīng)過(guò)裂縫地層后,在裂縫地層中的檢波器接收到的qP波的垂向偏振方向?yàn)榧t色射線所指,假設(shè)檢波點(diǎn)以上地層各向同性,P波的垂向偏振方位為藍(lán)色射線所指,這兩個(gè)偏振方向之間有相對(duì)偏差,我們利用這個(gè)角的相對(duì)變化來(lái)表示各向異性:

(6)

式中,Pani是偏振分析的qP波相對(duì)垂直方向?qū)嶋H偏振角;Piso是在假設(shè)地層是各向同性介質(zhì)條件下,通過(guò)射線追蹤的方法計(jì)算接收點(diǎn)處的射線垂直入射角;ψr的大小表示各向異性的程度。如果通過(guò)橢圓擬合ψr值,那么和公式(3)有所區(qū)別的是,通過(guò)相對(duì)偏振處理后,橢圓的長(zhǎng)軸指示的方向?yàn)榱芽p的走向方向。

圖4 qP波偏振極性示意

2.5 橢圓擬合

經(jīng)過(guò)對(duì)振幅、速度和極性進(jìn)行處理,消除了各種觀測(cè)痕跡和上覆地層的影響后,這些屬性就能反映目標(biāo)地層的裂縫引起的各向異性情況。因此,分別就目標(biāo)層中同一深度、不同方位的屬性通過(guò)最小二乘橢圓進(jìn)行擬合,然后計(jì)算橢圓長(zhǎng)軸方向及長(zhǎng)短軸比值。橢圓擬合時(shí)可以對(duì)屬性中的一些異常值進(jìn)行剔除,有必要時(shí)做適當(dāng)?shù)钠交幚怼Ｗ詈髮?duì)于目標(biāo)地層,統(tǒng)計(jì)各個(gè)屬性計(jì)算的橢圓長(zhǎng)軸方向和長(zhǎng)短軸比,綜合分析后得到裂縫的主方位和各向異性程度。

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

實(shí)際應(yīng)用中，我們采集了山西某煤層一口氣井的walkaround VSP數(shù)據(jù)。該數(shù)據(jù)共有24炮,從正北順時(shí)針?lè)较蚱?～345°,每15°一炮(圖5);10m間距32級(jí)三分量檢波器接收;接收井段在距井口650～960m深度處。目標(biāo)層位于富含煤層氣的790～850m深度的山西組和860～960m深度的太原組。由于工區(qū)的山地地形條件復(fù)雜,從圖5中可以看到,炮點(diǎn)偏移距不是等距離的,同時(shí)炮點(diǎn)高程也不在同一水平面上,這樣的觀測(cè)條件為各向異性分析帶來(lái)了困難。

在裂縫分析之前對(duì)原始數(shù)據(jù)進(jìn)行水平分量旋轉(zhuǎn)、徑向和垂向分量旋轉(zhuǎn)、波場(chǎng)分離等預(yù)處理,最后得到下行qP波場(chǎng)(圖6)。

首先對(duì)振幅進(jìn)行分析。對(duì)下行qP波場(chǎng),我們選擇每一炮最淺層650m深度處檢波點(diǎn)記錄作為該炮記錄的參考道,對(duì)整炮記錄反褶積。圖7是最淺層檢波點(diǎn)反褶積前、后的記錄,從圖中看到,反褶積后記錄能量得到均衡,后續(xù)的多次波基本消失。然后我們?cè)诔踔粮浇?0ms的一個(gè)時(shí)窗里,提取目標(biāo)層不同方位反褶積后qP波的RMS振幅(圖8b)。對(duì)比反褶積前目標(biāo)層qP波的RMS振幅(圖8a),可以看到原始RMS振幅和炮點(diǎn)的方位表現(xiàn)出高度相關(guān)性,而且不同方位振幅級(jí)別相差很大;反褶積后的不同方位qP波RMS振幅還存在一定的相關(guān)性,但是不同方位之間能量的級(jí)別有所減小。如果仔細(xì)觀察圖8b和圖5就會(huì)發(fā)現(xiàn),偏移距小、高程低的方位上振幅較大;反之振幅較小。也就是說(shuō),炮點(diǎn)偏移距和高程對(duì)反褶積后RMS振幅影響還比較嚴(yán)重,必須通過(guò)計(jì)算偽透射系數(shù)(PTC)的辦法消除。從計(jì)算的反褶積后RMS振幅的PTC值(圖8c)可以看出,不同方位的PTC值與炮點(diǎn)方位相關(guān)性得以消除,此時(shí)可以認(rèn)為PTC是裂縫各向異性引起的qP波振幅的變化。最后通過(guò)橢圓擬合不同深度各個(gè)方位的PTC值(圖9),并統(tǒng)計(jì)的山西組地層的各向異性的平均值為7.5%,裂縫走向?yàn)?2.5°和127.0°,太原組地層各向異性平均值為6.4%,裂縫主走向97.0°。

圖5 walkaround VSP炮點(diǎn)偏移距分布(a)以及炮點(diǎn)的高程分布(b)

圖6 預(yù)處理分離出的下行qP波部分炮集

圖7 最淺層檢波點(diǎn)道集(深度650m)反褶積前(a)和反褶積后(b)記錄

圖8 qP波原始RMS振幅(a);反褶積后的RMS振幅(b)以及qP波偽透射系數(shù)(c)

圖9 振幅各向異性橢圓擬合結(jié)果

隨后對(duì)速度進(jìn)行分析。首先拾取qP波初至?xí)r間,通過(guò)圖10a的雷達(dá)圖看到,初至?xí)r間依然與炮點(diǎn)偏移距及高程(圖5)存在很強(qiáng)的相關(guān)性,這是不能利用走時(shí)分析方法分析復(fù)雜地形條件下walkaround VSP資料各向異性的主要原因。因此利用初至和偏振分析得到的偏振信息計(jì)算與各向異性有關(guān)的局部層速度,從局部層速度(圖10b)中看到,炮點(diǎn)偏移距及高程與局部層速度的相關(guān)性基本消除。結(jié)合橢圓擬合的結(jié)果(圖11),統(tǒng)計(jì)在山西組和太原組地層中,各向異性平均值都為5.6%左右,而且裂縫的主方位角都為127°。

最后對(duì)極性進(jìn)行分析。通過(guò)偏振分析得到qP波相對(duì)垂直方向的偏振極性。炮點(diǎn)偏移距和高程對(duì)偏振極性的影響雖然不像其對(duì)振幅和旅行時(shí)的影響那么明顯(圖12a),但是還是需要經(jīng)過(guò)相對(duì)偏振處理預(yù)以消除。消除后的偏振極性(圖12b)的變化就表示裂縫地層各向異性特征。圖13是橢圓擬合的偏振極性各向異性結(jié)果,從結(jié)果看橢圓長(zhǎng)軸基本都指向北西方向。最后通過(guò)統(tǒng)計(jì),山西組地層的各向異性程度平均值為7.1%,裂縫主方位角為112°,太原組地層各向異性程度平均值為5.6%,裂縫主方位為142°。

圖14是對(duì)3種屬性計(jì)算的裂縫參數(shù)綜合分析的總結(jié)。目標(biāo)層振幅、速度和偏振極性所表現(xiàn)的方位各向異性雖然存在具體差異,但是基本方位總體一致。綜合統(tǒng)計(jì)后,山西組地層的裂縫方位集中指向120°,與通過(guò)橫波分裂算法計(jì)算的該層的裂縫方位120°(圖15)完全一致。太原組地層的裂縫主方位127°也與橫波分裂算法計(jì)算結(jié)果一致。對(duì)于裂縫引起的各向異性程度,3種屬性的各向異性程度曲線趨勢(shì)很好的吻合(尤其對(duì)于太原組地層)。對(duì)3種屬性的各向異性程度值求平均，得到山西組地層的各向異性平均值為6.7%,太原組地層的各向異性程度平均值為5.8%。

圖10 qP波初至波旅行時(shí)(a)以及局部速度(b)

圖11 速度各向異性橢圓擬合結(jié)果

圖12 qP波原始偏振極性(a)及其相對(duì)偏振極性(b)

圖13 偏振極性各向異性橢圓擬合結(jié)果

圖14 通過(guò)振幅、速度和極性綜合分析的裂縫參數(shù)

圖15 橫波分裂方法計(jì)算的目標(biāo)層的裂縫方位

4 結(jié)論

文中結(jié)合實(shí)際walkaround VSP資料研究了通過(guò)下行qP波振幅、速度和偏振極性的方位各向異性特征分析儲(chǔ)層裂縫參數(shù)的方法和流程。應(yīng)用研究表明,炮點(diǎn)激發(fā)條件、地形條件和上覆地層中各種因素嚴(yán)重影響了walkaround VSP資料方位各向異性特征,進(jìn)而影響到對(duì)目標(biāo)層裂縫參數(shù)的預(yù)測(cè),因此消除這些觀測(cè)痕跡的影響是利用qP波的各種屬性正確分析地層裂縫參數(shù)的關(guān)鍵。文中通過(guò)最淺層檢波點(diǎn)反褶積和PTC處理有效地校正了目的層上覆地層對(duì)振幅的影響。通過(guò)引入計(jì)算與各向異性有關(guān)的局部速度方法分析各向異性能避開(kāi)炮點(diǎn)位置和上覆地層對(duì)旅行時(shí)時(shí)差影響較難消除的弊端。在偏振極性分析各向異性時(shí),提出的相對(duì)偏振極化方法消除了炮點(diǎn)位置對(duì)偏振極性的影響。最后通過(guò)最小二乘橢圓擬合3種屬性得到裂縫的主方位和各向異性程度參數(shù),并進(jìn)行綜合統(tǒng)計(jì),能避免單個(gè)屬性的分析誤差,提高了統(tǒng)計(jì)結(jié)果的可靠性。最終3種屬性的統(tǒng)計(jì)分析表明,統(tǒng)計(jì)后目標(biāo)層的裂縫主方位清晰集中,和橫波分裂算法計(jì)算的裂縫方位能很好吻合,而且3種屬性的各向異性曲線也有很好的一致性。這證明通過(guò)qP波的振幅、速度和偏振極性屬性分析裂縫參數(shù)的方法和流程可靠有效,在實(shí)際生產(chǎn)中具有應(yīng)用價(jià)值。

[1] 劉衛(wèi)華,雍學(xué)善,楊午陽(yáng),等.縱波地震疊前方位數(shù)據(jù)裂縫預(yù)測(cè)方法及應(yīng)用[C]∥中國(guó)地球物理學(xué)會(huì).中國(guó)地球物理(2011).合肥:中國(guó)科學(xué)技術(shù)大學(xué)出版社,2011:445 Liu W H,Yong X S,Yang W Y,et al.Method and application of fracture prediction with the P-wave azimuth prestack seismic data[C]∥Chinese Geophysical Society.The Chinese geophysical (2011).Hefei:University of Science and Technology of China Press,2011:445

[2] 曹立斌,李亞林,梁波,等.三分量VSP資料在裂縫檢測(cè)中的應(yīng)用[J].石油物探,2008,47(1):67-71 Cao L B,Li Y L,Liang B,et al.Application of three-component VSP data in fracture detection[J].Geophysical Prospecting for Petroleum,2008,47(1):67-71

[3] 解建建.基于橫波分裂分析的裂縫檢測(cè)方法研究及應(yīng)用[D].成都:成都理工大學(xué),2011 Xie J J.The research and application of fracture detection methods based on shear wave splitting[D].Chengdu:Chengdu University of Technology,2011

[4] 薛梅.橫波分裂方法研究[C]∥中國(guó)地球物理學(xué)會(huì).中國(guó)地球物理(2008).北京:中國(guó)大地出版社,2008:438 Xue M.Shear wave splitting method research [C]∥ Chinese Geophysical Society.The Chinese geophysical (2008).Beijing:China Land Press,2008:438

[5] 羅省賢,李錄明.基于橫波分裂的地層裂縫預(yù)測(cè)方法與應(yīng)用[J].成都理工大學(xué)學(xué)報(bào)(自然科學(xué)版),2003,30(1):52-59 LuoS X,Li L M.Method and application of fredicting medium fracture based on shear-wave selarating theory[J].Journal of Chengdu University of Technology(Science & Technology Edition),2003,30(1):52-59

[6] 黎書(shū)琴,李忠,張白林.利用橫波分裂預(yù)測(cè)裂縫[J].石油地球物理勘探,2009,44(S1):130-134 Li S Q,Li Z,Zhang B L.Application of shear wave splitting to predict fracture[J].Oil Geophysical Prospecting,2009,44( S1):130- 134

[7] 候安寧,馬在田.垂直地震剖面橫波分裂的理論研究[J].同濟(jì)大學(xué)學(xué)報(bào)(自然科學(xué)版),1996,24(5):574-579 Hou A N,Ma Z T.Shear wave splitting in VSP’s data[J].Journal of TongJi Universitity (Natural Science),1996,24(5):574-579

[8] 呂功訓(xùn),劉合年,鄧常念,等.多方位非零偏三分量垂直地震剖面方法在阿姆河右岸裂縫檢測(cè)中的應(yīng)用[J].石油天然氣學(xué)報(bào)(江漢石油學(xué)院學(xué)報(bào)),2014,36(7):47-64 Lv G X,Liu H N,Deng C N,et al.Application of multiple azimuth 3-component nonzero-offset VSP methods in fracture detection on the right bank of Amudarya Basin[J].Journal of Oil and Gas Technology(Journal of Jianghan Petroleum Institute),2014,36(7):47-64

[9] 蒲靜,秦啟榮.油氣儲(chǔ)層裂縫預(yù)測(cè)方法綜述[J].特種油氣藏,2008,15(3):9-13Pu J,Qin Q R.An overview of fracture prediction methods for oil and gas reservoirs[J].Special Oil & Gas Reservoirs,2008,15(3):9-13

[10] 尹志恒,狄?guī)妥?李向陽(yáng),等.國(guó)外應(yīng)用縱波各向異性技術(shù)檢測(cè)裂縫的研究進(jìn)展[J].科技導(dǎo)報(bào),2011,29(30):73-79 Yin Z H,Di B R,Li X Y,et al.Progress in P-wave anisotropy technology for fracture detection[J].Science & Technology Review,2011,29(30):73-79

[11] 劉振峰,曲壽利,孫建國(guó),等.地震裂縫預(yù)測(cè)技術(shù)研究進(jìn)展[J].石油物探,2012,51(2):191-198 Liu Z D,Qu S L,Sun J G,et al.Progress of seismic fracture characterization technology[J].Geophysical Prospecting for Petroleum,2012,51(2):191-198

[12] 王允清,侯亞彬,劉金平,等.縱波方位各向異性在河道砂體識(shí)別中的應(yīng)用[J].石油物探,2007,46 (3):267-271 Wang Y Q,Hou Y B,Liu J P,et al.Application of P-wave azimuthal anisotropy in channel-sand recognition[J].Geophysical Prospecting for Petroleum,2007,46(3):267-271

[13] Hudson J A.Wave speeds and attenuation of elastic waves in material containing cracks[J].Geophysical Jonrnal of the Royal Astronomical Society,1981,64(1):133-150

[14] Owusu J C,Mubarak S M,Liu Q L.High-fidelity walkaround VSP anisotropy analysis[J].The Leading Edge,20009,28(8):966-972

[15] Craft K L,Mallick S,Meister L J,et a1.Azimuthal anisotropy analysis from P wave seismic traveltime data[J].Expanded Abstracts of 67thAnnual Internat SEG Mtg,1997,1214-1217

[16] Mallick S,Frazer L N.Reflection/transmission coefficients and azimuthal anisotropy in marine seismic studies[J].Geophysical Jonrnal International,1991,105(1):241-252

[17] 甘其剛,楊振武,彭大鈞.振幅隨方位角變化裂縫檢測(cè)技術(shù)及其應(yīng)用[J].石油物探,2004,43(4):373-376 Gan Q G,Yang Z W,Peng D J.AVA fracture detection and its application in JM area[J].Geophysical Prospecting for Petroleum,2004,43(4):373-376

[18] 張立勤,彭蘇萍,李國(guó)發(fā),等.方位AVO技術(shù)檢測(cè)儲(chǔ)層各向異性的方法和實(shí)踐[J].天然氣工業(yè),2005,25(10):38-39 Zhang L Q,Peng S P,Li G F,et al.The method and practice of detecting reservoir anisotropy with 3-D P-wave azimuth AVO[J].Natural Gas Industry,2005,25(10):38-39

[19] 楊勤勇,趙群,王世星,等.縱波方位各向異性及其在裂縫檢測(cè)中的應(yīng)用[J].石油物探,2006,45(2):177-181 Yang Q Y,Zhao Q,Wang S X,et al.P-wave azimuthal anisotropy its application in detection frcature[J].Geophysical Prospecting for Petroleum,2006,45(2):177-181

[20] 杜惠平.利用縱波方位AVO技術(shù)進(jìn)行裂縫檢測(cè)[J].新疆石油地質(zhì),2008,29(5):569-571 Du H P.Fracture detection with 3D P-wave azimuthal AVO[J].Xinjiang Petroleum Geology,2008,29(5):569-571

[21] 甘其剛,高志平.寬方位AVA裂縫檢測(cè)技術(shù)應(yīng)用研究[J].天然氣工業(yè),2005,25(5):42-43 Gan Q G,Gao Z P.Application of AVA frcatrue detection technique[J].Natural Gas Industry,2005,25(5):42-43

[22] 黃躍,馬昭軍,甘其剛,等.基于HTI介質(zhì)理論的縱波速度分析方法[J].成都理工大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,38(4):478-480 Huang Y,Ma Z J,Gan Q G,et al.P wave velocity analysis method based on the HTI media theory[J].Journal of Chengdu University of Technology(Science & Technology Edition),2011,38(4):478-480

[23] Qian Z P,Li X Y,Chapman M.Azimuthal variations of PP- and PS-wave attributes:a synthetic study[J].Expanded Abstracts of 77thAnnual Internat SEG Mtg,2007,184-188

[24] Grechka V,Mateeva A,Gentry C.Estimation of local seismic anisotropy from P-wave VSP data[J].The Leading Edge,2007,26(6):756-759

[25] 高志凌,王舉紅.非零偏VSP層速度的計(jì)算方法[C]∥中國(guó)地球物理學(xué)會(huì).中國(guó)地球物理(2006).成都:四川科學(xué)技術(shù)出版社,2006:72 Gao Z L,Wang J H.Non-zero offset VSP Interval velocity calculation method[C]∥ Chinese Geophysical Society.The Chinese geophysical (2006).Chengdu:Sichuan Publishing House of Science & Technology,2006:72

(編輯：朱文杰)

Fracture detection method based on walkaround VSP and its application

Chen Zhanguo,Chen Lin,Zhang Weihong

(SinopecGeophysicalResearchInstitute,Nanjing211103,China)

Research on the development of reservoir fractures is to reduce exploration risk and effectively develop hydrocarbon reservoirs.The analysis method and workflow on fracture parameters are studied by using the amplitude,velocity and polarity azimuthal anisotropy of down-going qP wave from walkaround VSP data.The impact of overlaying formations on the amplitude is removed by deconvolution and pseudo-transmission coefficient (PTC) processing.The impact of shot point elevation,offset and overlaying formations on travelling-time is avoided by calculating local velocity related with anisotropy.The impact of shot position on the polarization of polarity is corrected by the proposed relative polarization angle method.Finally,we analyze the variation characteristics of the amplitude,velocity and polarity azimuthal anisotropy of down-going qP wave via azimuth by least square ellipse fitting method.Moreover,the statistics of fracture strike and development in reservoirs is carried out.The application indicates that the analysis results of the attributes match with each other and are consistent with the results obtained from S-wave splitting method,which proves the higher reliability of the method and is effective for fracture prediction by walkaround VSP data.

walkaround VSP,fracture,anisotropy,amplitude,velocity,polarity

2015-04-21;改回日期：2015-07-16。

陳占國(guó)(1985—),男,碩士,助理工程師,現(xiàn)主要從事VSP相關(guān)方法研究。

國(guó)家科技重大專(zhuān)項(xiàng)項(xiàng)目“煤層氣井中地球物理勘探技術(shù)研究”(2011ZX05035003)資助。

P631

A

1000-1441(2015)06-0745-10

10.3969/j.issn.1000-1441.2015.06.013