孫 煒，李玉鳳，何巍巍，黑繼虎，孔凱波

(1.中國(guó)地質(zhì)大學(xué)地球物理信息技術(shù)學(xué)院，北京100083; 2.中國(guó)石化勝利油田分公司物探研究院，山東東營(yíng)257022; 3.中國(guó)石油勘探開(kāi)發(fā)研究院西北分院，甘肅蘭州730020; 4.中國(guó)石化中原油田分公司采油二廠，河南濮陽(yáng)457532; 5.中國(guó)石化中原油田普光分公司天然氣凈化廠，四川宣漢636150)

在碳酸鹽巖裂縫型油氣藏中，裂縫可以溝通溶孔和溶洞、有效提高儲(chǔ)層的滲流能力，因此，裂縫發(fā)育帶往往指示碳酸鹽巖儲(chǔ)層中的高產(chǎn)區(qū)。以方位各向異性理論為基礎(chǔ)的縱波方位各向異性裂縫預(yù)測(cè)技術(shù)成為近年來(lái)裂縫預(yù)測(cè)問(wèn)題的研究熱點(diǎn)［1-2］。

方位各向異性，是指地震波動(dòng)力學(xué)屬性隨著方位角的改變而發(fā)生變化。關(guān)于各向異性裂縫預(yù)測(cè)的相關(guān)研究，國(guó)外最早開(kāi)始于20 世紀(jì)70年代，經(jīng)過(guò)四十余年的不斷發(fā)展，已成為一套成熟的理論知識(shí)體系［3-12］。國(guó)內(nèi)的方位各向異性裂縫預(yù)測(cè)起步于20 世紀(jì)末，隨著裂縫型儲(chǔ)層的不斷發(fā)現(xiàn)，國(guó)內(nèi)的石油勘探工作者們?cè)絹?lái)越多的將其應(yīng)用于裂縫預(yù)測(cè)研究工作中［13-25］。本文以ZY 區(qū)塊碳酸鹽巖裂縫型儲(chǔ)層為例，進(jìn)行P 波各向異性裂縫預(yù)測(cè)方法的應(yīng)用實(shí)踐，通過(guò)測(cè)井裂縫及試油資料驗(yàn)證該方法在碳酸鹽巖儲(chǔ)層裂縫預(yù)測(cè)研究中的可靠性和有效性。

1 各向異性裂縫預(yù)測(cè)原理及思路

目前研究的介質(zhì)，通常是具有對(duì)稱(chēng)軸的介質(zhì)，即橫向各向同性介質(zhì)(transversely isotropy，簡(jiǎn)稱(chēng)TI)。TI 介質(zhì)，根據(jù)其對(duì)稱(chēng)軸在空間定向是垂直還是水平又分為VTI(transversely isotropy with a vertical axis of symmetry)介質(zhì)和HTI(transversely isotropy with a horizontal axis of symmetry)介質(zhì)［26］。VTI 介質(zhì)一般認(rèn)為是由周期性薄互層產(chǎn)生;HTI 介質(zhì)一般由平行排列的垂直裂隙、裂縫產(chǎn)生，本文的相關(guān)研究均基于HTI 介質(zhì)展開(kāi)。

1.1 方位各向異性裂縫預(yù)測(cè)原理

Ruger［10-11］基于弱各向異性的概念，并結(jié)合Thomsen 的各向異性系數(shù)，得到各向異性介質(zhì)中反射系數(shù)隨方位角和入射角變化的公式:

在小入射角的前提下，Ruger 對(duì)(1)式進(jìn)一步簡(jiǎn)化，并引入方位角變化的梯度項(xiàng)B(φk)，使反射系數(shù)與方位角變化的梯度項(xiàng)B(φk)建立關(guān)系，B(φk)由各向同性項(xiàng)系數(shù)Biso以及各向異性項(xiàng)系數(shù)Bani組成，它們的表達(dá)式分別為:

式中:φk為第k 次的觀測(cè)方位角，(°);φsym為沿裂縫帶對(duì)稱(chēng)軸方向的方位角，(°)，在沒(méi)有橫波數(shù)據(jù)的情況下一般為未知的。

在進(jìn)行寬方位地震采集時(shí)，可以得到多個(gè)φk，對(duì)3 個(gè)未知數(shù)Biso，Bani和φsym進(jìn)行非線性方程組的求解，并根據(jù)各向同性時(shí)(無(wú)裂縫)的Biso和各向異性時(shí)(裂縫發(fā)育)的Biso+Bani，來(lái)擬合方位各向異性橢圓，用橢圓長(zhǎng)短軸來(lái)指示裂縫方向，用橢圓長(zhǎng)短軸之比來(lái)確定表征裂縫發(fā)育情況的各向異性強(qiáng)度。Mallick 和Craft［9］等人研究認(rèn)為，裂縫強(qiáng)度越大，由各向異性擬合出的方位橢圓的扁率越扁，其長(zhǎng)軸或短軸方向代表裂縫走向，從而實(shí)現(xiàn)對(duì)裂縫強(qiáng)度和方向的預(yù)測(cè)。這樣，在單界面反射波的假設(shè)前提下，具有3 個(gè)方位或3 個(gè)方位以上的地震反射振幅數(shù)據(jù)就能夠進(jìn)行地層中某位置裂縫發(fā)育程度和方向的求取。

根據(jù)賀振華的研究成果［27］，對(duì)于裂縫，尤其是高角度裂縫，地震波動(dòng)力學(xué)特征(頻率、振幅、衰減和相對(duì)波阻抗等)比運(yùn)動(dòng)學(xué)特征更為敏感。研究表明，在這些動(dòng)力學(xué)屬性中，衰減屬性與裂縫的關(guān)系更為密切［12］。因此，研究中選用表征衰減屬性的85%能量比對(duì)應(yīng)頻率值來(lái)進(jìn)行各向異性裂縫預(yù)測(cè)，該屬性計(jì)算公式:

式中:Etotal為有效頻帶內(nèi)總能量;E85%為有效頻帶內(nèi)總能量的85%;F 算子代表在頻率能量坐標(biāo)系中，求取某能量值對(duì)應(yīng)的頻率值;f85%即為所求85%能量比對(duì)應(yīng)頻率值。

1.2 技術(shù)流程

首先，討論P(yáng) 波方位各向異性裂縫預(yù)測(cè)方法使用的地震數(shù)據(jù)。純波和成果數(shù)據(jù)由于經(jīng)過(guò)了全方位角疊加處理，提高地震數(shù)據(jù)信噪比的同時(shí)丟失了方位角的信息，因此純波和成果數(shù)據(jù)均不適合該方法的應(yīng)用;與之類(lèi)似，經(jīng)過(guò)疊前時(shí)間偏移的CRP 道集數(shù)據(jù)，由于丟失了方位角的信息，同樣不能用于方位各向異性的研究。在實(shí)際應(yīng)用中，選擇未做過(guò)疊前偏移，僅做了預(yù)處理、靜校正和動(dòng)校正的疊前NMO 道集數(shù)據(jù)來(lái)進(jìn)行該方法的裂縫預(yù)測(cè)研究。

其次，為了驗(yàn)證該方法的準(zhǔn)確性，需要有可靠的測(cè)井裂縫先驗(yàn)信息。在所有測(cè)井資料中，取心資料和FMI(微電阻率掃描成像測(cè)井)資料提供的裂縫信息最為可靠。因此，參考研究區(qū)的巖心資料以及FMI 測(cè)井資料，分析裂縫的井上分布特征。

該研究的主要發(fā)現(xiàn)有三個(gè)：第一，《新視野》整體AWL覆蓋率為4.63%，介于通用英語(yǔ)和學(xué)術(shù)英語(yǔ)之間，更傾向于通用英語(yǔ)。第二，《新視野》四冊(cè)書(shū)的AWL覆蓋率分別為4.30%、5.41%、3.67%和5.12%。變化范圍為3.67%到5.41%，介于1.4%至10%之間，各冊(cè)之間波動(dòng)幅度較大，沒(méi)有遵循由低到高、循序漸進(jìn)的規(guī)律。第三，每?jī)?cè)內(nèi)部單元間AWL覆蓋率數(shù)值變化較大，甚至由部分篇章語(yǔ)言接近學(xué)術(shù)英語(yǔ)。

圖1 基于HTI 介質(zhì)各向異性正演的裂縫預(yù)測(cè)方法的技術(shù)思路Fig.1 Technical process of fracture prediction based on anisotropy forwarding modeling of HTI media

在確定了關(guān)鍵技術(shù)及使用數(shù)據(jù)后，制定了如下的技術(shù)思路(圖1):①首先進(jìn)行測(cè)井資料的分析，確定井上裂縫的分布特征;②對(duì)NMO 地震道集數(shù)據(jù)進(jìn)行分方位角疊加，得到幾個(gè)含有方位角信息的分方位角疊加數(shù)據(jù)體，然后分別計(jì)算各分方位角疊加數(shù)據(jù)體的衰減屬性，并進(jìn)行方位各向異性橢圓的擬合，預(yù)測(cè)研究區(qū)裂縫的分布特征;③利用巖心和FMI 測(cè)井資料從剖面和平面來(lái)檢測(cè)裂縫預(yù)測(cè)的準(zhǔn)確性，根據(jù)驗(yàn)證情況確定最終的裂縫分布特征。

2 實(shí)例應(yīng)用

2.1 地質(zhì)背景

研究區(qū)構(gòu)造上位于T 油田ZY 區(qū)塊F 盆地古隆起的斜坡中央?yún)^(qū)，目的層為石炭系M 組碳酸鹽巖儲(chǔ)層，地層厚度約280 m。研究區(qū)目的層儲(chǔ)層以孔縫洞復(fù)合型為主，且裂縫以高角度裂縫為主，裂縫的溝通作用對(duì)儲(chǔ)層的貢獻(xiàn)很大。在單井上表現(xiàn)為高產(chǎn)井孔隙度值并不突出，而滲透率較高;且孔隙度高，滲透率偏低的井試油結(jié)果普遍不理想。

準(zhǔn)確的井上裂縫先驗(yàn)信息是進(jìn)行裂縫預(yù)測(cè)研究的基礎(chǔ)。區(qū)內(nèi)有五口有FMI 測(cè)井資料的井，分別為A1，A2，A3，A4，A5，并對(duì)A1 井進(jìn)行了取心觀察。根據(jù)巖心資料和FMI 測(cè)井資料的分析，研究區(qū)碳酸鹽巖儲(chǔ)層裂縫以構(gòu)造縫為主，在不同層段也見(jiàn)溶蝕縫發(fā)育。

2.2 方位各向異性裂縫預(yù)測(cè)

首先對(duì)疊前地震道集數(shù)據(jù)進(jìn)行偏移距與方位角的分析，按照滿覆蓋和寬方位角的要求，選取合適的偏移距和方位角范圍，依據(jù)能量均衡(由不同方位角范圍疊加得到的若干疊加數(shù)據(jù)體在相同側(cè)線或聯(lián)絡(luò)線剖面上無(wú)明顯的能量差異)的原則合理進(jìn)行方位角的劃分，通過(guò)分方位角疊加得到含有方位角信息的5 個(gè)方位角疊加數(shù)據(jù)體(圖2)，并對(duì)各方位角疊加數(shù)據(jù)體進(jìn)行頻譜分析(圖3)。根據(jù)頻譜分析結(jié)果，在有效頻帶內(nèi)計(jì)算各分方位角疊加數(shù)據(jù)體的衰減屬性，通過(guò)各方位角的衰減屬性擬合，得到研究區(qū)裂縫的分布特征。結(jié)合FMI成像測(cè)井及巖心的裂縫先驗(yàn)信息，從剖面和平面探討該方法的準(zhǔn)確性。

圖2 分方位角疊加數(shù)據(jù)剖面Fig.2 Section of five azimuthal-stack data volume

圖3 各方位角疊加數(shù)據(jù)體的頻譜分布Fig.3 Spectrum of each azimuthal-stack data volume

首先觀察裂縫強(qiáng)度過(guò)井剖面。研究區(qū)的測(cè)井?dāng)?shù)據(jù)中，只有A1 井同時(shí)具有FMI 測(cè)井及巖心資料，因此重點(diǎn)考察過(guò)A1 井裂縫強(qiáng)度剖面(圖4上)，并通過(guò)局部放大的方式，將預(yù)測(cè)的裂縫結(jié)果與井點(diǎn)處FMI 測(cè)井以及巖心裂縫進(jìn)行對(duì)比分析(圖4 下)。在過(guò)井剖面中，裂縫發(fā)育分布較為散亂，規(guī)律性較差;從井點(diǎn)局部放大對(duì)比圖中可以看到，預(yù)測(cè)裂縫強(qiáng)度大的位置對(duì)應(yīng)FMI 測(cè)井及巖心上的裂縫發(fā)育區(qū)(局部放大圖中的FMI 和巖心長(zhǎng)度不到1 m，是由于鉆井單塊巖心樣品一般不超過(guò)1 m，放大圖中的選取厚度包含在裂縫強(qiáng)度剖面紅色區(qū)域的15 ms 中)，預(yù)測(cè)的裂縫強(qiáng)度與FMI 測(cè)井及巖心裂縫信息吻合較好，產(chǎn)油井段與裂縫發(fā)育段的位置基本對(duì)應(yīng)。

圖4 過(guò)A1 井衰減各向異性預(yù)測(cè)裂縫密度剖面(上)及預(yù)測(cè)結(jié)果與測(cè)井裂縫對(duì)比(下)Fig.4 Attenuation anisotropy-based fracture density section crossing Well A1 (shown above)and comparison between the prediction results and log-based fracture detection(bottom)

然后，分析研究區(qū)的裂縫平面分布特征。對(duì)裂縫強(qiáng)度數(shù)據(jù)體進(jìn)行平面成圖，得到裂縫的平面分布規(guī)律(圖5)。圖中左側(cè)色標(biāo)表示裂縫發(fā)育情況，代表各向異性強(qiáng)度，即裂縫密度，無(wú)量綱，暖色代表裂縫較為發(fā)育;預(yù)測(cè)的裂縫方向通過(guò)玫瑰圖來(lái)指示(玫瑰圖中色標(biāo)代表百分比，即某一方向的裂縫占該位置目的層段所有裂縫方向的百分比，玫瑰圖中的紅色指示預(yù)測(cè)出的裂縫優(yōu)勢(shì)方向)。

研究區(qū)內(nèi)發(fā)育斷層若干條，裂縫的發(fā)育區(qū)與斷層的展布特征關(guān)系密切，且在主要斷層之間(圖5 中黑線范圍內(nèi))，裂縫尤為發(fā)育。同時(shí)根據(jù)FMI測(cè)井資料的裂縫顯示，A2 井全井段無(wú)裂縫;A3 井在不同深度發(fā)育裂縫兩條，裂縫不發(fā)育;A1，A4 和A5 井裂縫均較為發(fā)育，從圖中可以看出預(yù)測(cè)的裂縫平面分布規(guī)律與FMI 成像測(cè)井的裂縫發(fā)育情況吻合較好。研究區(qū)裂縫的方向以北東向?yàn)橹?，北北東為輔，裂縫方向受斷層走向控制。

綜合來(lái)看，研究區(qū)裂縫以構(gòu)造縫為主，與井上的認(rèn)識(shí)相符合。根據(jù)預(yù)測(cè)的裂縫平面分布特征，在裂縫發(fā)育區(qū)的邊界位置設(shè)定了新的井位SJ1 井(圖5)，鉆井完成后，在目的層段試油獲日產(chǎn)43 t工業(yè)油流，從而進(jìn)一步驗(yàn)證了預(yù)測(cè)方法的有效性。

3 結(jié)論

圖5 衰減各向異性預(yù)測(cè)研究區(qū)裂縫方向平面分布特征Fig.5 Plane view of azimuth distribution of fractures predicted based on attenuation anisotropy in the study area

圖6 衰減各向異性預(yù)測(cè)研究區(qū)裂縫方向與FMI 井上裂縫方向?qū)Ρ菷ig.6 Comparison of fracture azimuth obtained from FMI with those from attenuation anisotropy prediction

1)在利用該方法進(jìn)行裂縫預(yù)測(cè)時(shí)，方位角的劃分非常關(guān)鍵。對(duì)疊前縱波地震道集進(jìn)行方位角的劃分時(shí)，要注意以下兩點(diǎn):①至少要?jiǎng)澐殖? 個(gè)方位，以保證各向異性橢圓的擬合;②進(jìn)行劃分時(shí)要保證各分方位角疊加數(shù)據(jù)體的能量較為均衡，覆蓋次數(shù)大致相同，減少由地震數(shù)據(jù)方位覆蓋不均勻而帶來(lái)的各向異性差異，導(dǎo)致預(yù)測(cè)結(jié)果不準(zhǔn)。

2)利用衰減各向異性來(lái)預(yù)測(cè)裂縫，要考慮到AVO 現(xiàn)象對(duì)裂縫預(yù)測(cè)的影響。如果隨著方位角的增加能量增大或減小，而同時(shí)，隨著入射角的增加能量變化情況恰恰相反，即減小或增大，這時(shí)的裂縫預(yù)測(cè)問(wèn)題將變得更為復(fù)雜。

3)在利用縱波方位各向異性進(jìn)行裂縫預(yù)測(cè)時(shí)，由于Ruger 的公式(1)描述的是單界面HTI 介質(zhì)下的反射系數(shù)變化情況，因此，在進(jìn)行裂縫預(yù)測(cè)研究時(shí)，該方法主要適用于高角度裂縫的分布規(guī)律預(yù)測(cè)。

致謝:感謝中國(guó)地質(zhì)大學(xué)(北京)王彥春教授、中國(guó)科學(xué)院地質(zhì)與地球物理研究所白志明副研究員、史小東博士對(duì)文章提出的寶貴意見(jiàn)和建議。

［1］劉振峰，曲壽利，孫建國(guó).地震裂縫預(yù)測(cè)技術(shù)研究進(jìn)展［J］.石油物探，2012，51(2):191-198.Liu Zhenfeng，Qu Shouli，Sun Jianguo.Progress of seismic fracture characterization technology ［J］. Geophysical Prospecting for Petroleum，2012，51(2):191-198.

［2］尹志恒，狄?guī)妥專(zhuān)钕蜿?yáng)，等.國(guó)外應(yīng)用縱波各向異性技術(shù)檢測(cè)裂縫的研究進(jìn)展［J］.科技導(dǎo)報(bào)，2011，29(30):73-78.Yin Zhiheng，Di Bangrang，Li Xiangyang，et al. Progress in Pwave anisotropy technology for fracture detection［J］. Science and Technology Review，2011，29(30):73-78.

［3］Crampin S.A review of the effects of anisotropic layering on the propagation of seismic waves［J］.Geophysics，1977，49:9-27.

［4］Crampin S.Effective anisotropic elastic constants for wave propagation through cracked solids［J］. Geophysics，1984，76:135-145.

［5］Hudson J A. Wave speeds and attenuation of elastic waves in material containing cracks［J］. Geophysics，1981，64:133-150.

［6］Hudson J A.A higher-order approximation to the wave propagation constants for a cracked solid［J］. Geophysics，1986，87:265-274.

［7］Thomsen L.Weak elastic anisotropy［J］.Geophysics，1986，51:1954-1966.

［8］Thomsen L. Elastic anisotropy due to aligned cracks in porous rock［J］.Geophys Prosp，1995，43:805-829.

［9］Mallick S，Craft K，Meister L，et al.Determination of the principle direction of azimuthal anisotropy from P-wave seismic data［J］.Geophysics，1998，(63):692-706.

［10］Rüger A. Variation of P-wave reflectivity with offset and azimuth in anisotropic media［J］. Geophysics，1998，63(3):935-947.

［11］Rüger A. Using AVO for fracture detection:analytic basis and practical solutions［J］. The Leading Edge，1997，16(10):1429-1430，1432-1434.

［12］Shen F，Zhu X and Toksoz N. Anisotropy of aligned fractures and P-wave azimuthal AVO response［J］. SEG，1997:2001-2004.

［13］曲壽利，季玉新，王鑫，等. 全方位P 波屬性裂縫檢測(cè)方法研究［J］.石油地球物理勘探，2001，36(4):390-397.Qu Shouli，Ji Yuxin，Wang Xin，et al.Seismic method for using full-azimuth P-wave attributes to detect fracture［J］. Oil Geophysical Prospecting，2001，36(4):390-397.

［14］范國(guó)章，牟永光，金之鈞.裂縫介質(zhì)中地震波方位AVO 特征分析［J］.石油學(xué)報(bào)，2002，23(4):42-45.Fan Guozhang，Mou Yongguang，Jin Zhijun. Amplitude versus offset characteristics of azimuth of seismic wave in fractured reservoir［J］.Acta petrolei sinica，2002，23(4):42-45.

［15］樂(lè)紹東.AVA 裂縫檢測(cè)技術(shù)在川西JM 構(gòu)造的應(yīng)用［J］.天然氣工業(yè)，2004，24(4):22-24 Yue Shaodong.Application of AVA fracture detection technique in JM structure in west Sichuan［J］. Natural gas industry，2004，24(4):22-24.

［16］高云峰，李緒宣，陳桂華，等. 裂縫儲(chǔ)層地震預(yù)測(cè)技術(shù)在錦州25-1 南潛山的應(yīng)用研究［J］. 中國(guó)海上油氣，2008，20(1):22-26.Gao Yunfeng，Li Xuxuan，Chen Guihua，et al.An application of seismic techniques to predict fractured reservoirs in JZ25-1S buried hill［J］. China Offshore Oil and Gas，2008，20(1):22-26

［17］喻岳鈺，楊長(zhǎng)春，王彥飛，等. 瞬時(shí)頻域衰減屬性及其在碳酸鹽巖裂縫檢測(cè)中的應(yīng)用［J］. 地球物理學(xué)進(jìn)展，2009，24(5):1717-1722.Yu Yueyu，Yang Changchun，Wang Yanfei，et al. P-wave azimuthal attenuation attributes in wavelet-scale domain and its application to fracture detection in carbonate［J］. Progress in Geophysics，2009，24(5):1717-1722.

［18］黃偉傳，楊長(zhǎng)春，王彥飛.利用疊前地震數(shù)據(jù)預(yù)測(cè)裂縫儲(chǔ)層的應(yīng)用研究［J］.地球物理學(xué)報(bào)，2007，22(5):1602-1606.Huang Weichuan，Yang Changchun，Wang Yanfei.The application of pre-stack seismic data in predicting the fractured reservoir［J］.Progress in Geophysics，2007，22(5):1602-1606.

［19］孫煒，王彥春，李梅，等. 利用疊前地震數(shù)據(jù)預(yù)測(cè)火山巖儲(chǔ)層裂縫［J］.物探與化探，2010，34(2):229-232.Sun Wei，Wang Yanchun，Li Mei，et al. The detection of fractures in volcanic reservoir with pre-stack seismic data ［J］.Geophysical and Geochemical Exploration，2010，34(2):229-232.

［20］姜傳金，鞠林波，張廣穎，等. 利用地震疊前數(shù)據(jù)預(yù)測(cè)火山巖裂縫的方法和效果分析——以松遼盆地北部徐家圍子斷陷營(yíng)城組火山巖為例［J］. 地球物理學(xué)報(bào)，2011，54(2):515-523 Jiang Chuanjin，Ju Linbo，Zhang Guangying，et al. The method and effect analysis of volcanic fracture prediction with prestack seismic data——An example from the volcanic rocks of Yingcheng formation in Xujiaweizi fault depression，north of Songliao Basin［J］. Chinese Journal of Geophysics，2011，54(2):515-523.

［21］白斌，鄒才能，朱如凱，等. 四川盆地九龍山構(gòu)造須二段致密砂巖儲(chǔ)層裂縫特征、形成時(shí)期與主控因素［J］.石油與天然氣地質(zhì)，2012，33(4):526-535.Bai Bin，Zou Caineng，Zhu Rukai，et al. Characteristics，timing and controlling factors of structural fractures in tight sandstones of the 2ndmember of Xujiahe Formation in Jiulong Mountain structure，Sichuan Basin［J］.Oil＆Gas Geology.2012，33(4):526-535.

［22］陳冬，陳力群，魏修成，等. 火成巖裂縫性儲(chǔ)層測(cè)井評(píng)價(jià)——以準(zhǔn)噶爾盆地石炭系火成巖油藏為例［J］. 石油與天然氣地質(zhì)，2011，32(1):83-90.Chen Dong，Chen Liqun，Wei Xiucheng，et al.Log evaluation of fractured igneous reservoirs:a case study of the Carboniferous igneous reservoirs in the Junggar Basin［J］. Oil ＆ Gas Geology，2011，32(1):83-90.

［23］黃保綱，汪利兵，趙春明，等.JZS 油田潛山裂縫儲(chǔ)層形成機(jī)制及分布預(yù)測(cè)［J］.石油與天然氣地質(zhì)，2011，32(5):710-717.Huang Baogang，Wang Libing，Zhao Chunming，et al.Formation mechanism and distribution of fractured buried hill reservoir in JZS oilfiled［J］.Oil ＆ Gas Geology，2011，32(5):710-717

［24］涂彬，丁祖鵬，劉月田.火燒山油田H3 層裂縫發(fā)育特征與剩余油分布關(guān)系［J］. 石油與天然氣地質(zhì)，2011，32(2):229-235.Tu Bin，Ding Zupeng and Liu Yuetian.Correlation of remaining oil distribution and fracture development in H-3 reservoir of Huoshaoshan oilfield［J］. Oil ＆ Gas Geology，2009，30(2):229-235.

［25］葛瑞·馬沃可(美).巖石物理手冊(cè):孔隙介質(zhì)中地震分析工具［M］. 合肥:中國(guó)科學(xué)技術(shù)大學(xué)出版社，2008，113-118.Gary Mavko. The rock physics handbook:Tools for seismic analysis in porous media［M］. Hefei:University of Science and Technology of China Press，2008，113-118.

［26］吳國(guó)忱.各向異性介質(zhì)地震波傳播與成像［M］.東營(yíng):中國(guó)石油大學(xué)出版社，2006，24～25.Wu Guochen.Propagation and imaging of seismic wave in anisotropy media［M］. Dongying:China University of Petroleum Press，2006，24-25.

［27］賀振華，黃德濟(jì)，文曉濤. 裂縫油氣藏地球物理預(yù)測(cè)［M］.成都:四川科學(xué)技術(shù)出版社，2007.He Zhenhua，Huang Deji，Wen Xiaotao.Geophysical predictive theory and technique on fractured reservoir［M］.Chengdou:Sichuan Publishing House of Science and Technology，2007.