崔炯成

(中海油研究總院)

高壓地層條件下AVO異常陷阱研究

崔炯成

(中海油研究總院)

根據(jù)南海鶯瓊盆地某井高壓地層的典型砂巖樣品實(shí)驗(yàn)室?guī)r心測(cè)試數(shù)據(jù)，分析了砂巖速度和密度等參數(shù)隨壓力變化的特征，在此基礎(chǔ)上開(kāi)展了隨機(jī)AVO正演模擬研究，獲得了在不同壓力情況下不同巖性、物性和流體組合的AVO響應(yīng)特征。在圍壓不變情況下，縱、橫波速度隨孔隙壓力的增大而降低，而密度受壓力的影響較小;當(dāng)砂巖地層壓力大于上覆泥巖地層壓力時(shí)，飽和水砂巖可能出現(xiàn)第Ⅳ類AVO異常假象，因此在高壓地層中進(jìn)行儲(chǔ)層預(yù)測(cè)或烴類檢測(cè)時(shí)應(yīng)結(jié)合地層壓力預(yù)測(cè)和隨機(jī)AVO正演模擬研究來(lái)成功規(guī)避AVO異常陷阱。

高壓地層AVO異常陷阱隨機(jī)正演模擬鶯瓊盆地

勘探實(shí)踐表明，在常壓地層條件下利用AVO分析方法可以識(shí)別孔隙流體［1-3］，而在高(超)壓盆地中高壓對(duì)AVO特征的影響仍然是一個(gè)棘手的問(wèn)題［4］，例如我國(guó)南海鶯瓊盆地中大量存在非油氣引起的亮點(diǎn)和AVO異常。Verm等人［5］通過(guò)對(duì)墨西哥灣超過(guò)2 700口井?dāng)?shù)據(jù)的研究發(fā)現(xiàn)，超壓對(duì)AVO特征有重要影響，超壓地層頂界面之下的AVO響應(yīng)特征與頂界面之上的AVO響應(yīng)特征存在很大差異。因此，巖性及流體識(shí)別必須結(jié)合壓力預(yù)測(cè)研究成果。Rick Lindsay等人［6］發(fā)現(xiàn)，異常地層壓力會(huì)導(dǎo)致AVO異常由第Ⅲ類變?yōu)榈冖纛?Lloyd Weathers等人［7］通過(guò)對(duì)墨西哥灣38口深井研究發(fā)現(xiàn)，當(dāng)砂巖位于壓力過(guò)渡帶時(shí)，干砂巖的AVO特征可能出現(xiàn)第Ⅱ類異常假象，而產(chǎn)氣砂巖往往無(wú)AVO異常;Antonio等人［8］通過(guò)模型研究表明，相鄰砂泥巖地層間壓力的突變會(huì)導(dǎo)致AVO異常假象出現(xiàn)。筆者通過(guò)分析我國(guó)南海某高壓盆地砂巖樣本實(shí)驗(yàn)室測(cè)試數(shù)據(jù)，得到了縱、橫波速度和密度隨壓力變化的特征;在此基礎(chǔ)上，利用隨機(jī)AVO正演模型方法研究不同壓力、巖性和物性組合條件下的AVO響應(yīng)特征，探討了高壓盆地可能出現(xiàn)的與高壓相關(guān)的AVO異常陷阱。

1 縱、橫波速度和密度隨壓力變化特征

通過(guò)對(duì)南海鶯瓊盆地某井高壓地層的典型砂巖樣品進(jìn)行實(shí)驗(yàn)室?guī)r石物理測(cè)試分析，研究縱、橫波速度和密度隨壓力變化的特征。測(cè)試中采用固定圍壓(pc)，通過(guò)改變巖樣孔隙壓力(pp)測(cè)量得到不同差壓(pd=pc－pp)時(shí)聲波速度的變化特征。圖1為該井飽和水砂巖聲波速度隨差壓變化圖，可見(jiàn)縱、橫波速度均隨孔隙壓力增大(差壓降低)而呈現(xiàn)降低的趨勢(shì)，但在孔隙壓力較低時(shí)速度的變化不大。研究認(rèn)為，這種速度-壓力關(guān)系主要因巖石內(nèi)部微縫隙的開(kāi)閉狀態(tài)及巖石骨架顆粒接觸的致密程度隨壓力變化造成［9-12］。

圖1 南海鶯瓊盆地某井高壓地層典型飽和水砂巖樣品縱、橫波速度隨差壓變化特征

圖2為該井飽和水砂巖樣品密度隨孔隙壓力變化圖，可見(jiàn)測(cè)試樣品的密度隨孔隙壓力升高而略微增高;在不同圍壓與孔隙壓力的組合下，密度的變化均不超過(guò)0.5%，即測(cè)試樣品的密度受壓力的影響很小。另外，其他學(xué)者的研究成果也有類似結(jié)論，但只適用于已固結(jié)的巖石［13］。

圖2 不同圍壓條件下密度隨孔隙壓力的變化特征

2 高壓地層條件下AVO特征分析

為了進(jìn)一步分析壓力對(duì)AVO特征的影響，筆者對(duì)低壓泥巖-高壓砂巖及高壓泥巖-低壓砂巖組合(圖3)條件下地層壓力對(duì)砂巖頂界面反射AVO特征的影響進(jìn)行了正演模擬研究。

圖3 相鄰砂泥巖地層的2種壓力組合

2.1 低壓泥巖與高壓砂巖組合

選取鶯瓊盆地某井三亞組泥巖作為“低壓”泥巖代表，陵水組砂巖作為“高壓”砂巖代表，保持泥巖參數(shù)不變，隨機(jī)均勻抽樣選取砂巖的孔隙度(10%～25%)、泥質(zhì)含量(0～40%)、含水飽和度(0～100%)，代入R-H-G(Raymer-Hunt-Gardner)公式求得飽和水狀態(tài)下對(duì)應(yīng)的縱波速度

式(1)中:vp為巖石的縱波速度，km/s;φ為總孔隙度，%;vps為巖石礦物的縱波速度，km/s;vpf為巖石孔隙流體的縱波速度，km/s。

利用Greenberg-Castagna公式可以求得飽和水狀態(tài)下對(duì)應(yīng)的橫波速度

式(2)中:vs_wet為飽和水狀態(tài)下的橫波速度，km/s;L為巖石中的礦物數(shù)量;Xi為每個(gè)巖石礦物的體積分量;Ni為成分i的多項(xiàng)指數(shù);aij為實(shí)驗(yàn)回歸系數(shù); vp_wet為飽和水狀態(tài)下的縱波速度，km/s。

再利用Gassmann方程計(jì)算出含水飽和度變化時(shí)對(duì)應(yīng)的速度，同時(shí)計(jì)算出密度對(duì)應(yīng)的變化

式(3)中:Ksat為飽和巖石的有效體積模量，GPa;Kdry為干巖石有效體積模量，GPa;K0為巖石礦物的有效體積模量，GPa;φ為巖石總孔隙度，%;Kfl為孔隙流體的有效體積模量，GPa;μsat為飽和巖石的有效剪切模量，GPa;μdry為干巖石的有效剪切模量，GPa。

最后將速度和密度代入Aki-Richards方程中求得一系列反射系數(shù)隨角度(偏移距)的變化特征

式(4)中:α=(vp1+vp2)/2;β=(vs1+vs2)/2;ρ=(ρ1+ρ2)/2;Δα=vp2－vp1;Δβ=vs2－vs1;Δρ=ρ2－ρ1; Rpp為縱波入射縱波反射的反射系數(shù);下標(biāo)1、2分別代表上、下地層;θ為入射角，(°)。

按照上述方法計(jì)算得到低壓泥巖-高壓砂巖組合在砂巖飽和水與含氣條件下的隨機(jī)AVO特征，結(jié)果如圖4、5所示，圖中顏色代表不同孔隙度范圍內(nèi)計(jì)算結(jié)果，其中黑色代表孔隙度為10%～15%，藍(lán)色為15%～20%，紅色為20%～25%。由圖4可見(jiàn)，在低壓泥巖-高壓砂巖組合下，砂巖飽和水時(shí)可出現(xiàn)第Ⅳ類AVO異常陷阱(紅色部分);而在層間壓力平衡時(shí)，相同巖性物性條件下卻沒(méi)有出現(xiàn)明顯的第Ⅳ類AVO異常。由圖5可見(jiàn)，與層間壓力平衡時(shí)相比，低壓泥巖-高壓砂巖組合在砂巖含氣后AVO梯度絕對(duì)值有所減小(需要注意的是，層間壓力平衡時(shí)的第Ⅲ類AVO異常受到壓力的影響后將變?yōu)榈冖纛怉VO異常)。

2.2 高壓泥巖與低壓砂巖組合

選取該井梅山組泥巖作為“高壓”泥巖代表，三亞組砂巖作為“低壓”砂巖代表，采用相同方法計(jì)算得到高壓泥巖-低壓砂巖組合在砂巖飽和水與含氣條件下的隨機(jī)AVO特征，結(jié)果如圖6、7所示。

由圖6、7可見(jiàn)，與層間壓力平衡時(shí)相比，此種壓力組合在飽和水條件下沒(méi)有導(dǎo)致AVO異常出現(xiàn)，只是截距整體有所提高;但在砂巖含氣后，這種層間壓力差導(dǎo)致AVO截距有所提高，AVO梯度絕對(duì)值增大，第Ⅱ類和第Ⅲ類AVO異常更加明顯。

3 分析與討論

上述研究結(jié)果表明，在低壓泥巖-高壓砂巖組合下，水層可能出現(xiàn)第Ⅳ類AVO異常，在剖面上表現(xiàn)為“亮點(diǎn)”特征，這種現(xiàn)象主要由高壓砂巖的低速特征造成，此時(shí)若砂巖含氣，由于上覆泥巖速度相對(duì)較大，導(dǎo)致砂泥巖泊松比差異降低，更容易出現(xiàn)第Ⅳ類AVO異常(由Shuey方程可知)。在高壓盆地中，若砂巖高壓為它源機(jī)制，則其地層壓力一般要大于周圍泥巖，此時(shí)要注意這種地層間的壓力差異導(dǎo)致的AVO異常陷阱，這是因?yàn)楦邏旱貙泳哂小胺糯蟆钡卣甬惓５墓δ埽畬踊蚝瑲馑畬釉诟邏簵l件下常常表現(xiàn)出“亮點(diǎn)”和AVO異常，因此，在高壓地層進(jìn)行油氣勘探時(shí)需要結(jié)合地層壓力預(yù)測(cè)結(jié)果和地震正演模擬研究來(lái)排除這種假象。

4 結(jié)束語(yǔ)

巖心測(cè)試結(jié)果表明，異常高孔隙壓力會(huì)使縱、橫波速度降低，但對(duì)于已固結(jié)巖石密度的影響不大。AVO隨機(jī)正演模擬結(jié)果表明，地層間壓力的突變會(huì)明顯改變AVO反射特征，特別是在低壓泥巖-高壓砂巖組合下，含水砂巖可能會(huì)表現(xiàn)為第Ⅳ類AVO異常，高壓使得水層頂部反射表現(xiàn)為強(qiáng)反射，造成“亮點(diǎn)”假象。因此，在高壓地層中進(jìn)行儲(chǔ)層預(yù)測(cè)或烴類檢測(cè)時(shí)要充分考慮到壓力的影響，并結(jié)合壓力預(yù)測(cè)成果規(guī)避可能出現(xiàn)的假異常陷阱。

［1］謝濤，楊凱，金明霞，等．利用多方位地震采集資料預(yù)測(cè)BZ28-1油田裂縫儲(chǔ)層［J］．中國(guó)海上油氣，2011，23(3):158-162．

［2］張衛(wèi)衛(wèi)，何敏，顏承志，等．AVO疊前反演技術(shù)在珠江口盆地白云深水區(qū)烴類檢測(cè)中的應(yīng)用［J］．中國(guó)海上油氣，2010，22(3):160-163．

［3］田立新，周東紅，明君，等．渤海Q32-6北地區(qū)稠油層預(yù)測(cè)［J］．中國(guó)海上油氣，2010，22(1):12-16．

［4］BOX R，DOSS E．Typical AVO response as a function of depth and fluid pressure gradient:Gulf of Mexico shelf［J］．The Leading Edge，2008，27(10):1252-1262．

［5］VERM R，LIANG L，HILTERMAN F J．Significance of geopressure in predicting lithology［J］．The Leading Edge，1998，17(2): 227-234．

［6］LINDSAY R，TOWNER B．Pore pressure influence on rock property and reflectivity modeling［J］．The Leading Edge，2001，20 (2):184-187．

［7］WEATHERS L，LEROY S，LUNEN E V．AVO models in geopressure setting-deep Gulf of Mexico shelf［C］．SEG 2003，22:611-614．

［8］ANTONIO C B R，MARCO A S T．Induced AVO anomalies from pore pressure effects［C］．SEG 2004，23:252-255．

［9］NUR A，SIMMONS G．Stress-induced velocity anisotropy in rocks:an experimental study［J］．Journal of Geophysical Research，1969，74(27):6667．

［10］NUR A．Effects of stress on velocity anisotropy in rocks with rocks［J］．Journal of Geophysical Research，1971，76(8):2022-2034．

［11］SAYER C M．Stress-induced ultrasonic wave velocity anisotropy in fractured rock［J］．Ultrasonics，1988，26(6):311-317．

［12］MARKO G，MUKERJI T，GODFREY N．Predicting stress-induced velocity anisotropy in rocks［J］．Geophysics，1995，60(4):1081-1087．

［13］MACBETH C．A classification for the pressure-sensitivity properties of a sandstone rock frame［J］．Geophysics，2004，69(2)，497-510．

A study on AVO anomaly pitfalls in overpressured formations

Cui Jiongcheng
(CNOOC Research Institute，Beijing，100027)

Based on the laboratory testing data of cores exacted from a typical overpressured sandstone formation in Ying-Qiong basin，South China Sea，the changes of the sandstone parameters，such as velocity and density，with pressure were analyzed．Then，a stochastic forward modeling of AVO was carried out，by which the characteristic AVO responses were recognized for different combinations of lithology，petrophysics and fluids at different pressures．If the confining pressure is fixed，the velocity of P-wave and S-wave will decrease with increasing pore pressure，and the rock density will be more insensitive to pore pressure．When the pressure in a sandstone formation is higher than that in its overlying shale formation，a falseⅣtype of AVO anomaly may occur in the sandstone if it is saturated with water．Therefore，the formation pressure prediction and the stochastic forward modeling of AVO should be integrated into conducting reservoir prediction or hydrocarbon detection for overpressured formations，in order to successfully avoid the AVO anomaly pitfalls．

overpressured formation;AVO anomaly pitfall;stochastic forward modeling;Ying-Qiong basin

2011-11-14改回日期:2012-02-05

(編輯:周雯雯)

崔炯成，男，高級(jí)工程師，1984年畢業(yè)于同濟(jì)大學(xué)海洋地球物理勘探專業(yè)，2008年獲中國(guó)科學(xué)院地質(zhì)與地球物理研究所博士學(xué)位，現(xiàn)從事油氣田勘探開(kāi)發(fā)項(xiàng)目管理及研究工作。地址:北京市東城區(qū)東直門(mén)外小街6號(hào)海油大廈(郵編:100027)。