陳兆明

(中海石油(中國)有限公司深圳分公司 廣東廣州 510240)

惠州凹陷東沙隆起區(qū)三角洲前緣砂體受后期抬升隆起形成上傾尖滅，NE—SW向巖性尖滅線與NW—SE向構(gòu)造等值線之間形成的巖性圈閉是珠江口盆地進行地層巖性圈閉勘探最有利目標［1］。然而，受地震分辨率的限制，巖性尖滅線的位置很難準確落實，導(dǎo)致巖性圈閉封閉的范圍很難確定，迫切需要對本地區(qū)的有利砂體進行油的檢測，以證實巖性圈閉封閉性與成藏條件，減小勘探風險。

目前烴類檢測的方法有“亮點”、“平點”、“暗點”技術(shù)，AVO技術(shù)，吸收系數(shù)、阻抗、分形、波形、三瞬技術(shù)，多波多分量勘探技術(shù)，儲層反演烴類特征重構(gòu)，聚類分析，模式識別神經(jīng)網(wǎng)絡(luò)等技術(shù)［2-9］，這些技術(shù)可分為疊前、疊后屬性兩大類。隨著地震屬性分析從疊后走向疊前，20世紀90年代以來應(yīng)用P波反射特征隨偏移距變化(AVO)識別流體已成為研究熱點，并廣泛應(yīng)用于國內(nèi)外碎屑巖、碳酸鹽巖乃至致密砂巖儲層［2-8］，珠江口盆地通過地震剖面“亮點”或AVO異常識別也在番禺-流花和白云深水區(qū)取得了一系列天然氣的重大發(fā)現(xiàn)［9-10］。

疊前AVO特征分析以前多用于氣的檢測［10］，2013年汪瑞良 等［11］在珠江口盆地中淺層油藏中通過提取實際地震資料的甜點屬性、均方根振幅和反射強度交流分量等疊后屬性實現(xiàn)了對富油層段的識別。在此基礎(chǔ)上，本文結(jié)合惠州地區(qū)原油低密度特性，開展了疊前AVO特征分析在該地區(qū)原油檢測中的應(yīng)用工作，預(yù)測了目標井的含油性，與鉆井情況吻合較好。

1 疊前AVO特征分析方法

當儲層孔隙流體不同時，反射PP波的反射系數(shù)隨入射角度變化的特征不同，前人據(jù)此劃分出了4種AVO特征，并從AVO特征判斷儲層中的含油性。疊前AVO特征分析方法有以下2個關(guān)鍵步驟。

1)針對目的層巖石物理特征，利用式(1)對測井數(shù)據(jù)進行正演，分析已鉆井的AVO特征，初步判定研究區(qū)的疊前AVO特征。當入射角θ＜50°時，假定vP=2vS，Zoeppritz方程式中反射PP波的反射系數(shù)可寫成Shuey近似形式［12］

式(1)中：θ為入射角度；IP1、IP2分別為上下介質(zhì)的縱波阻抗，(IP1-IP2)/(IP1+IP2)反映縱波穿過地層界面時的相對變化；IS1、IS2分別為上下介質(zhì)的橫波阻抗，(IS1-IS2)/(IS1+IS2)反映橫波穿過地層界面時的相對變化；ρ1、ρ2分別為上下介質(zhì)的密度，(ρ1-ρ2)/(ρ1+ρ2)反映介質(zhì)密度在地層界面的相對變化。

2)對目的層進行流體替換，再利用式(1)進行正演，分析流體變化對AVO特征的影響，進一步界定研究區(qū)的疊前AVO特征。基于Gassmann方程［13-14］進行流體替換的步驟為：

①由測井數(shù)據(jù)得到流體1的縱波速度vP1、橫波速度vS1和體積密度ρ1。

②提取流體1的體積模量K1和剪切模量μ1，即

③進行流體替換，即

式(4)中：φ為孔隙度；Km為固體基質(zhì)體積模量(利用V-R-H模型計算)；Kfl1為孔隙中流體1的體積模量，Kfl2為孔隙中流體2的體積模量(均用Batzle-Wang方法計算)。

④保持剪切模量不變，即μ2=μ1。

⑤計算修正流體變化的體積密度，即

⑥得到流體2的速度，即

2 在惠州地區(qū)油檢測中的應(yīng)用

2.1 方法適用性分析

疊前AVO特征分析以前大多是針對氣的流體檢測，很少用于對油的流體檢測，這是由于通常情況下油與水的密度比較接近而很難區(qū)分。惠州地區(qū)M洼附近幾個油田地面原油密度大多在0.81～0.91 g/cm3，地層原油密度為 0.52～0.90 g/cm3，與水的密度差異較大(表1)。分析認為，這可能與M洼原油的氣油比較高有關(guān)，如靠近洼陷中心的H6油田的M32層地層原油密度低至0.52 g/cm3，原始氣油比高達345.4m3/m3，其他油田雖遠離洼陷中心，但DST測試時也產(chǎn)出不少氣量。另外，研究區(qū)主要目的層埋深為2 000～2 600 m，目的層段入射角在35°左右，小于50°，滿足Zoeppritz方程正演條件。

表1 惠州地區(qū)原油密度和氣油比統(tǒng)計結(jié)果Table 1 Statistics of crude oil density and gas-oil ratio in Huizhou area

2.2 疊前AVO特征分析

以A井為例，該井是研究區(qū)內(nèi)一口油井，其中L30up層是主要油層，L35、L45為水層。對該井進行疊前AVO正演模擬，結(jié)果顯示L30up油層為III類AVO特征，L35、L45水層呈IV類AVO特征(圖1)。通過對惠州地區(qū)多口井的正演模擬結(jié)果進行統(tǒng)計，初步判定，該區(qū)油層呈III類AVO特征，水層呈IV類AVO特征。

將A井L30up層從油層替換成水層，這時P波速度從3 570m/s增加到3 690m/s，增加約120 m/s，變化約3.3%；S波速度變化較?。幻芏葟?.25 g/cm3增加到 2.29 g/cm3，增加約 0.04 g/cm3，變化約1.8%(圖2)。沿L30up層底面向上提取AVO屬性，L30up層含油時為III類AVO特征，L30up層含水時為IV類AVO特征(圖2)。

將A井L35層從水層替換成油層，這時P波速度減小約200 m/s，變化約-5.7%；S波速度變化較?。幻芏葴p小約0.06 g/cm3，變化約-2.7%(圖3)；沿L35層底面向上提取AVO屬性，L35層含油時為III類AVO，L35層含水時為IV類AVO(圖3)。

流體替換結(jié)果表明，惠州地區(qū)油層一般呈III類AVO特征，水層呈IV類AVO特征，應(yīng)用地震數(shù)據(jù)體的疊前AVO特征分析能夠進行油檢測。

圖1 A井AVO正演結(jié)果Fig.1 AVO forward modeling of W ell A

圖2 A井L30up油層流體替換前后AVO正演結(jié)果Fig.2 AVO forward modeling before and after fluid substitution of L30up oil layer in W ell A

圖3 A井L35水層流體替換前后AVO正演結(jié)果Fig.3 AVO forward modeling before and after fluid substitution of L35 water layer in Well A

2.3 檢測效果

1)剖面檢測效果。B-1井L60層有明顯的III類AVO特征，L10low層有微弱的III類AVO特征，判斷為油層，實鉆結(jié)果顯示L60層、L10low層均為油層，說明檢測結(jié)果準確(圖4)。C井L60層、L10low層有明顯的IV類AVO特征，判斷為水層，實鉆結(jié)果與檢測結(jié)果相吻合(圖5)。

圖4 過B-1井疊前道集Fig.4 Pre-stack gathers acrossW ell B-1

圖5 過C井疊前道集Fig.5 Pre-stack gathers acrossW ell C

2)平面檢測效果。為了在平面上進行疊前AVO特征分析，沿L60層分別提取遠道數(shù)據(jù)體(10°～20°數(shù)據(jù)疊加)與近道數(shù)據(jù)體(1°～10°數(shù)據(jù)疊加)的振幅屬性(圖6a、b)，并用前者與后者相減(圖6c)，找到差值小于0的范圍，即遠角度數(shù)據(jù)體振幅加強的范圍(圖6d中綠色區(qū)域)，再找出近道數(shù)據(jù)體負振幅屬性的范圍(圖6d中紅色區(qū)域)，兩者重合的區(qū)域即為III類AVO區(qū)域?？梢钥闯觯珺-1井在III類AVO區(qū)域內(nèi)，判斷為油層；B-2井不在III類AVO區(qū)域內(nèi)，判斷為水層，與實鉆結(jié)果吻合。

圖6 HB構(gòu)造L60層振幅屬性Fig.6 Amp litude attribute of L60 layer in HB structure

3 結(jié)束語

針對惠州地區(qū)原油高氣油比、低密度的特性，利用疊前AVO特征分析方法界定了研究區(qū)的疊前AVO特征，該區(qū)油層一般呈III類AVO特征，水層呈IV類AVO特征。據(jù)此在該區(qū)無井區(qū)利用疊前地震數(shù)據(jù)進行AVO特征分析，有效區(qū)分了油層和水層，其檢測結(jié)果與鉆井結(jié)果相吻合，減小了勘探風險。

［1］龍更生，施和生，杜家元，等.珠江口盆地惠州地區(qū)中新統(tǒng)地層巖性圈閉形成條件分析［J］.中國海上油氣，2006，18(4)：229-235.Long Gengsheng，Shi Hesheng，Du Jiayuan，et al.An analysis of creation conditions for Miocene stratigraphic and lithologic traps in Huizhou area，Pearl River Mouth basin［J］.China Offshore Oil and Gas，2006，18(4)：229-235.

［2］秦秉讓.黔南碳酸鹽巖地震“亮點”、“平點”與“暗點”分析［J］.貴州地質(zhì)，1998，15(4)：338-344.Qin Bingrang.Analyses for bright，faded and darken spotson seismic profiles of carbonate rocks in Southern Guizhou［J］.Guizhou Geology，1998，15(4)：338-344.

［3］王達昌.陸相地層中亮點的檢測效果［J］.石油地球物理勘探，1982(4)：30-35.Wang Dachang.The detecting effect of bright spot from continental sedimentation［J］.Oil Geophysical Prospecting，1982(4)：30-35.

［4］柴緒兵.“亮點”、“暗點”技術(shù)在太平川南地區(qū)葡萄花油層油氣檢測中的應(yīng)用效果分析［J］.科學技術(shù)與工程，2012，12(9)：2132-2136.Chai Xubing.The effectiveness analysis of using“Bright Spots”，“Dim Spots”technique to testoiland gas in Putaohua reservoir of Taipingchuan region［J］.Science Technology and Engineering，2012，12(9)：2132-2136.

［5］生如巖，李振泉，孟陽.利用地震“亮點”技術(shù)實現(xiàn)淺層氣藏滾動勘探開發(fā)的一個實例［J］.中國海上油氣(地質(zhì))，1999，13(1)：41-47.Sheng Ruyan，Li Zhenquan，Meng Yang.A case study of progressive exploration and development of shallow gas reservoir utilizing technique of seismic bright spot［J］.China Offshore Oil and Gas(Geology)，1999，13(1)：41-47.

［6］張云鵬，李東，張新民，等.“平點反射”在濮城深層天然氣滾動勘探中的應(yīng)用［J］.江漢石油學院學報，2000，22(4)：21-23.Zhang Yunpeng，Li Dong，Zhang Xinmin，et al.Application of"Flat Spot Reflection"in deep gas rolling exploration in Pucheng oilfield［J］.Journal of Jianghan Petroleum Institute，2000，22(4)：21-23.

［7］吳超，王曉云，葉翠華，等.利用暗點技術(shù)確定氣藏氣水界面：以柯克亞E2K氣藏為例［J］.新疆石油地質(zhì)，2002，23(6)：528-530.Wu Chao，Wang Xiaoyun，Ye Cuihua，et al.Determination of gaswater boundary in gas reservoir using dim spot technique：an example of Kekeya E2k gas reservoir［J］.Xinjiang Petroleum Geology，2002，23(6)：528-530.

［8］巫芙蓉，李亞林，王聃，等.AVO屬性在致密砂巖儲層流體檢測中的應(yīng)用［J］.天然氣工業(yè)，2011，31(5)：55-57.Wu Furong，LiYalin，Wang Dan，etal.Application of AVO attribute to the detection of tight sandstone reservoirs［J］.Natural Gas Industry，2011，31(5)：55-57.

［9］軒義華，秦成崗，汪瑞良，等.分頻AVO技術(shù)在珠江口盆地番禺天然氣區(qū)含氣性分析中的應(yīng)用［J］.石油地球物理勘探，2010，45(1)：79-84.Xuan Yihua，Qin Chenggang，Wang Ruiliang，et al.Application of frequency division AVO technique in gas bearing analysis in Fanyu natural gas area，Pearl River Mouth Basin［J］.Oil Geophysical Prospecting，2010，45(1)：79-84.

［10］張衛(wèi)衛(wèi)，何敏，顏承志，等.AVO疊前反演技術(shù)在珠江口盆地白云深水區(qū)烴類檢測中的應(yīng)用［J］.中國海上油氣，2010，22(3)：160-167.Zhang Weiwei，He Min，Yan Chengzhi，et al.An application of AVO pre-stack inversion technique to DHI in Baiyun deep-water area，Pearl River Mouth basin［J］.China Offshore Oil and Gas，2010，22(3)：160-167.

［11］汪瑞良，劉軍，李麗，等.M凹陷富油層段地震響應(yīng)特征分析［J］.石油地球物理勘探，2013，48(5)：776-784.Wang Ruiliang，Liu Jun，Li Li，etal.Seismic response characteristics of rich reservoir in MDepression［J］.Oil Geophysical Prospecting，2013，48(5)：776-784.

［12］SHUEY R T.A simplification of the Zoeppritzequations［J］.Geophysics，1985，50(4)：609-614.

［13］GASSMANN F.Elastic waves through a Packing of spheres［J］.Geophysics，1951，16(4)：673-682.

［14］林凱，賀振華，熊曉軍，等.基于Gassmann方程的流體替換流程［J］.長江大學學報：自然科學版，2009，6(1)：180-182.Lin Kai，He Zhenhua，Xiong Xiaojun，et al.Based on the Gassmann equation fluid substitution Flow［J］.Journal of Yangtze University：Natural Science Edition，2009，6(1)：180-182.