劉明陽，張大海，劉 苗，周興海，李 倩

（中國石油化工股份有限公司上海海洋油氣分公司勘探開發(fā)研究院，上海 200120）

?

西湖凹陷W構造鉆井非均質儲層識別可行性分析

劉明陽，張大海，劉 苗，周興海，李 倩

（中國石油化工股份有限公司上海海洋油氣分公司勘探開發(fā)研究院，上海 200120）

摘 要：西湖凹陷W構造的主要目的層花港組砂體橫向變化快，埋藏深，屬于低孔低滲儲層。通過對構造區(qū)內已鉆井的巖石物理分析確定：低縱橫波速度比可以指示砂巖儲層；低縱橫波速度比、高泊松阻尼因子可以指示高含氣儲層。通過構造區(qū)已鉆井測井數(shù)據(jù)的AVO正反演研究可知：縱橫波速度比、泊松阻尼因子雙參數(shù)分析可以有效區(qū)分低孔低滲儲層和特低孔低滲儲層；且子波頻率越高，準確度越高。

關鍵詞：AVO正演；非均質性；低孔低滲儲層；儲層識別；泊松阻尼因子

1 研究背景

東海西湖凹陷是我國近海重要的油氣勘探區(qū)塊，根據(jù)構造、沉積等特征，西湖凹陷自西向東可分為西部斜坡帶、西部次凹帶、中央反轉構造帶、東部次凹帶和東部斷階帶［1］，目前已發(fā)現(xiàn)多個油氣田和含油氣構造。W構造位于西湖凹陷三潭區(qū)帶中北部，斜坡帶向凹陷中心的過渡帶，主要目的層為花港組?；ǜ劢M發(fā)育河流—湖泊—三角洲沉積體系，其中花港組上段為三角洲前緣沉積，下段為三角洲平原沉積，水下分流河道發(fā)育，砂體橫向變化快。

目前研究區(qū)已鉆遇五口井：X1井～X5井，完鉆井深最深達到5 300 m，儲層埋深在4 300 m左右，共解釋了H3、H5、H6、H7、H8、H9六套油氣層，主要目的層為H3、H5儲層。表1統(tǒng)計了五口已鉆井H3儲層砂體的厚度、孔隙度、滲透率，表2統(tǒng)計了X1井、X5井H5儲層砂體的厚度、孔隙度、滲透率。受壓實作用的影響，儲層孔隙度、滲透率都比較低，整體致密。

表1 W構造已鉆井H3儲層砂體厚度、孔隙度、滲透率統(tǒng)計

表2 W構造已鉆井H5儲層砂體厚度、孔隙度、滲透率統(tǒng)計

根據(jù)實際鉆井數(shù)據(jù)：在X1井H3儲層內部，隨著深度增加，鉆遇不同的沉積微相，巖性為細砂巖、含礫細砂巖、夾砂礫巖、中砂巖。

圖1為X1井巖心、壁心、測井數(shù)據(jù)的統(tǒng)計，由圖可知：儲層整體致密，隨深度增加，孔隙度、滲透率均呈減小趨勢，上下差別不大；全井段分析，H3儲層物性條件相對較好。

圖1 X1井巖心、壁心、測井孔隙度、滲透率特征

圖2為X1井H3儲層段巖心的實測孔隙度、滲透率分布。由圖可知：H3巖心實測孔隙度分布范圍為5% ～ 11.7%，平均值8%，垂向上差異不大；而實測滲透率在4 320 ～ 4 330 m之間均值為0.799×10-3μm2，在4 330 ～ 4 340 m之間均值為9.308×10-3μm2，上下差異較大，顯示H3單套儲層內部非均質性強。

圖2 X1井H3儲層段實測孔隙度、滲透率特征

結合測井解釋與勘探成果：確定了以下巖性劃分標準：泥巖：Vsh＞30%；特低孔特低滲儲層：0.2≤per≤0.5且φ≥7%；低孔低滲儲層：per ＞0.5且φ＞7.2%；干層：per＜0.2且φ＜7%。其中：Vsh為泥質質量分數(shù)；per為滲透率；φ為孔隙度。

通過前面分析可以發(fā)現(xiàn)：構造區(qū)內目的層埋深大，儲層整體致密，低孔低滲儲層與特低孔特低滲儲層差異小，識別難度大；垂向上，H3單套儲層內上部與下部之間物性特征差別比較大，非均質性強，地震屬性的敏感性區(qū)分難度大。

針對構造區(qū)目的層特征，主要研究內容為：確定巖性、物性、含氣性的敏感性參數(shù)；單套儲層內部低孔低滲儲層與特低孔低滲儲層的識別。

2 巖石物理分析

巖石物理分析技術是聯(lián)系地震彈性參數(shù)（速度、阻抗、振幅等）與儲層特性（孔隙度、滲透率、飽和度等）的橋梁，通過彈性參數(shù)的交會分析可以確定巖性敏感性參數(shù)，進一步可以確定孔隙度、滲透率高的有利儲層，是地震反演可行性分析中不可或缺的部分。

首先，研究彈性參數(shù)的巖性特征，圖3為X1井縱波阻抗與縱橫波速度比的交會圖，色標為泥質含量（紅色表示泥質含量高，為泥巖；藍色代表泥質含量低，為砂巖），從圖中可以看出，縱橫波速度比可以較好的區(qū)分砂巖、泥巖，識別儲層。

圖3 X1井彈性參數(shù)巖性分析

其次，對彈性參數(shù)的物性分析。此處主要用到了泊松阻尼因子［2］。泊松阻尼因子PDF（Poisson Dampening Factor）是縱波阻抗、橫波阻抗的函數(shù)：

式中：Ip為縱波阻抗，Is為橫波阻抗。

泊松阻尼因子越大，表明油氣豐度高或者儲層物性條件好［3］。

圖4為X1井泊松阻尼因子與縱橫波速度比的交會圖，色標分別代表孔隙度（藍色表示孔隙度低，紅色表示孔隙度高），交會分析可知泊松阻尼因子大值對應的孔隙度也比較大，反映儲層物性條件好。

圖4 X1井彈性參數(shù)物性分析

最后，在儲層識別、物性分析的基礎上，對儲層的含氣性進行分析。圖5為泊松阻尼因子與縱橫波速度比的交會圖，色標為含氣飽和度（紅色表示含氣飽和度高，藍色表示含氣飽和度低），分析可知，泊松阻尼因子越大、縱橫波速度比越低的位置，對應的含氣飽和度越高，劃定截止值，當PDF＞0.12和Vp/Vs＜1.65時含氣飽和度在40%左右，為有利儲層。

圖5 X1井彈性參數(shù)含油氣分析

3 疊前AVA同步反演原理及技術路線

疊前反演充分利用了疊前道集中的振幅信息，通過反演得到速度、密度、流體因子等彈性參數(shù)，提高了儲層預測的精度［4-5］。疊前反演的基本原理為：當平面波以非法向入射到介質分界面時，在界面上會產(chǎn)生四種波：反射PP波、反射SV波、透射PP波、透射SV波，可以用全Zoeppritz方程［6］精確地描述平面波入射介質分界面時反射系數(shù)、透射系數(shù)與界面上下介質參數(shù)之間的關系。但是由于完全Zoeppritz方程的非直觀性，在實際應用中更多的是其線性近似方程［7］，如公式（2）所示。這些近似方程形式簡潔，物理意義明確，是進行AVO反演、AVO交匯圖分析、烴類檢測的基礎。

式中：θ為入射角；Rpp（θ）當入射角為θ時的縱波反射系數(shù)；Vp為反射界面上下縱波速度和的平均數(shù)；Vs為反射界面上下橫波速度和的平均數(shù)；ρ為反射界面上下密度和的平均數(shù)；ΔVp、ΔVs、Δρ分別為界面上下地層縱波速度、橫波速度、密度之差。

通過對工區(qū)資料分析，結合疊前反演的基本原理，確定了如下技術路線，如圖6所示。主要的技術環(huán)節(jié)有：巖石物理分析，基于測井數(shù)據(jù)的道集正演，子波提取，低頻模型的建立，反演參數(shù)實驗，儲層識別。

圖6 疊前AVA同步反演技術路線

4 應用效果分析

以X1井為例，輸入的井曲線：縱波時差、橫波時差、密度。表3展示了目的層H3儲層的解釋結果，H3儲層上部（4 278.1 ～ 4 364.7 m）的孔隙度均值在7.8%左右，滲透率在0.75×10-3μm2左右，屬于低孔低滲儲層；H3儲層下部（4 364.7 ～ 4 416.4 m）的孔隙度均值在7.4%左右，滲透率在0.4×10-3μm2左右，屬于特低孔低滲儲層。

表3 X1井H3儲層物性及測井解釋結果

研究區(qū)目的層埋深在4 300 m左右，考慮到地震子波能量衰減，首先選取主頻為20 Hz的子波，正演角度范圍為0 ～ 45°（圖7）。在4 360 m有微弱反射軸，其對應的縱波阻抗上并沒有奇異值，屬于正常反射；提高子波頻率到30 Hz時（圖8），4 360 m反射軸加強，根據(jù)解釋成果，該軸是低孔低滲儲層與特低孔低滲儲層的分界面。

圖7 X1井AVO正演記錄（0 ～ 45°，子波主頻20 Hz）

依據(jù)巖石物理分析的結果，聯(lián)合泊松阻尼因子和縱橫波速度比對H3儲層非均質性進行識別。

圖9為子波頻率為20 Hz時的反演結果，從左到右依次為：孔隙度、滲透率、含氣飽和度、泊松阻尼因子、縱橫波速度比。由圖可知，在泊松阻尼因子相對較大的位置，對應的孔隙度、滲透率都比較大，根據(jù)巖石物理分析劃定的截止值PDF＞0.12且Vp/ Vs＜1.65，可以識別出兩套氣層，其含氣飽和度都在25%以上，厚度分別為90 m（4 280 ～ 4 370 m）、30 m（4 380 ～ 4 410 m），與解釋結果基本一致。圖10中提高子波頻率到30 Hz，同樣利用雙參數(shù)的截止值PDF＞0.12且Vp/ Vs＜1.65，可以識別的第一套氣層厚度為70 m （4 290 ～ 4 360 m），對應的含氣飽和度在40%以上，相較于子波頻率為20 Hz時，分辨率有所提高，同樣的對低孔低滲儲層和特低孔低滲儲層識別的準確度也提高。

圖8 X1井AVO正演記錄（0 ～ 45°，子波主頻30 Hz）

圖9 X1井反演結果（子波頻率20 Hz）

圖10 X1井反演結果（子波頻率30 Hz）

4.1 不確定性分析

地震反演是由地震數(shù)據(jù)求地下彈性參數(shù)的過程，由于環(huán)境、儀器、測井誤差等等因素的影響，其結果具有高度的不唯一性。本文主要研究了基于測井數(shù)據(jù)正反演的儲層非均質性識別方法，如果測井數(shù)據(jù)誤差較大，就會增加反演結果的多解性，進而增加儲層識別的不確定性。

4.2 結論

研究區(qū)目的層儲層埋深比較大，通過鉆井巖石物理分析可知：利用縱橫波速度比，可以區(qū)分砂泥巖，縱橫波速度比低值指示砂巖；聯(lián)合應用泊松阻尼因子與縱橫波速度比可以識別物性比較好、含氣飽和度比較高的儲層。在實際應用中，泊松阻尼因子與縱橫波速度比雙參數(shù)分析可以有效識別低孔低滲儲層與特低孔低滲儲層，且正反演子波頻率越高，識別越準確。

參考文獻：

［1］陳志勇，葛和平.西湖凹陷反轉構造與油氣聚集［J］.中國海上油氣（地質），2003，17（1）：20-24.

［2］MAZUMDARP P. Poisson Dampening Factor ［J］. The Leading Edge， 2007，26（7）： 850-852.

［3］高偉義，林貴康，李城堡，等.泊松阻尼因子在平湖地區(qū)儲層流體檢測中的應用——一種定量地震解釋的新方法［J］. 中國石油勘探，2013（2）：50-53.

［4］高偉義.平湖油氣田儲層及含油氣預測及應用［J］.海洋石油，2003，23（增刊）：7-12.

［5］姜勇，張磊，鄒偉，等.西湖凹陷B構造水下分流河道預測技術及應用［J］. 海洋石油，2015， 35（2）：35-39.

［6］ZOEPPRITZ K， ERDBEBENWELLEN V B. Uber Reflexion and durchgang seismischer wellen duch unstetigkeitsflachen［J］. Gottinger Nachr， 1919，17（1）：66-84.

［7］SHUEY R T. A simplification of the Zoeppritz equations［J］. Geophysics， 1985， 50（4）： 609-614.

中圖分類號：P631.4

文獻標識碼：A

DOI:10.3969/j.issn.1008-2336.2016.02.019

文章編號：1008-2336（2016）02-0019-06

收稿日期：2015-08-18；改回日期：2016-01-27

第一作者簡介：劉明陽，女，1987生，碩士，2014年畢業(yè)于中國石油大學（北京）地球探測與信息技術專業(yè)，從事地震反演研究工作。E-mail：liumy.shhy@sinopec.com。

Identifcation of the Heterogeneous Reservoirs in Structure W in Xihu Sag

LIU Mingyang， ZHANG Dahai， LIU Miao， ZHOU Xinghai， LI Qian
（Institute of Exploration and Development， SINOPEC Shanghai Offshore Oil & Gas Company， Shanghai 200120， China）

Abstract：The sandstone reservoir of Huagang Formation in Structure W in Xihu Sag is characterized by deep burial depth， lateral alteration in lithology， and classifed as the low porosity and low permeability reservoir. The analysis of the petrophysical properties in the wells drilled in this area shows that low P-S wave velocity ratio indicates the existence of sandstone reservoirs， low P-S wave velocity ratio and high Poisson dampening factor may indicate high possibility of gas bearing reservoir. Based on the study of AVO modeling and inversion on logging data， it can be concluded that the comprehensive analysis of P-S wave velocity ratio and Poisson damping factor can be used to effectively distinguish low porosity-low permeability reservoir and ultra-low porosity-low permeability reservoir. In addition， the higher accuracy can be expected with the increase of the wavelet frequency.

Keywords：AVO modeling； heterogeneity； low porosity and low permeability reservoir； reservoir identifcation； Poisson dampening factor