張 濤,佘 剛,李 苗,艾合買提江,葉雙江

(中國石油化工股份有限公司石油勘探開發(fā)研究院,北京100083)

A油田白堊系碳酸鹽巖裂縫型儲層綜合預測技術研究

張 濤,佘 剛,李 苗,艾合買提江,葉雙江

(中國石油化工股份有限公司石油勘探開發(fā)研究院,北京100083)

扎格羅斯盆地A油田白堊系裂縫-孔隙型碳酸鹽巖油藏厚度大于500m,裂縫分布復雜,孔隙受成巖控制,有效儲層預測難度大。針對上述難點,研究提出采用地質、測井、地震、油藏工程等多種信息對裂縫型碳酸鹽巖儲層進行綜合預測的思路。首先利用巖心的鑄體薄片及實測孔、滲資料從微觀上識別出裂縫的發(fā)育特征及其控制因素和發(fā)育模式,利用常規(guī)和成像測井定量解釋儲層的巖性、物性以及裂縫的產狀,為地震預測提供地質和巖石物理模型;然后優(yōu)選出預測裂縫的敏感疊后地震屬性(曲率屬性、對稱性)，利用縱波方位各向異性裂縫識別技術對儲層裂縫進行預測;最后借助與地震縱向分辨率相當?shù)脑嚲忉尦晒?、試油提供的流量和生產指數(shù)等動態(tài)資料,對各屬性及多屬性融合技術的預測結果進行了驗證。研究結果表明：①最大曲率和照明3D(Illuminator 3D,I3D)屬性在大尺度斷裂及裂縫發(fā)育帶識別方面具有優(yōu)勢；②縱波方位各向異性(AVAZ)裂縫檢測方法能夠提供關于裂縫密度及分布范圍方面較為詳實的信息；③地震多屬性融合技術可綜合各屬性的優(yōu)點,實現(xiàn)從小尺度裂縫到大尺度裂縫帶的預測,刻畫白堊系碳酸鹽巖裂縫型儲層的分布。

碳酸鹽巖;裂縫型儲層;多學科綜合預測;屬性融合;扎格羅斯盆地

碳酸鹽巖裂縫型油氣藏在世界范圍內廣泛分布,裂縫的發(fā)育分布情況與區(qū)域構造背景、褶皺強度、儲層巖性及厚度、層序組合等關系密切,儲集體非均質性強,分析、認識和預測裂縫發(fā)育帶分布規(guī)律是當今油氣勘探和開發(fā)中的重點和難點,相關的技術研究是業(yè)界的熱點。袁士義等[1]系統(tǒng)總結了裂縫的地質特征和裂縫空間分布預測技術;陳冬等[2]探討了裂縫在常規(guī)測井和成像測井曲線上的不同響應特征及響應機理,綜合利用雙側向電阻率和成像測井圖像識別裂縫發(fā)育段;劉振峰等[3]、曲壽利等[4]利用縱波資料開展方位各向異性檢測，并將其作為地震裂縫預測技術的重要組成部分;Ganguly等[5]提出的體曲率估算技術基于地震資料的傾角和方位角,避免了對連續(xù)解釋層位的需求,克服了二維層面曲率的局限性。齊晴等[6]、耿福蘭[7]和孔選林等[8]采用螞蟻追蹤技術、本征值相干技術和曲率屬性識別大尺度裂縫,采用基于方位各向異性的疊前地震屬性預測技術識別小尺度裂縫的方向和密度。本文以扎格羅斯盆地A油田白堊系巨厚、非均質性強的裂縫-孔隙型碳酸鹽巖儲層為目標,綜合地質控制因素、測井評價、疊前/疊后地震屬性、壓力恢復資料和試油成果等多種信息,開展碳酸鹽巖儲層預測技術研究,建立適合該油田的儲層綜合預測技術,為該油田開發(fā)技術的確定提供依據(jù)。

1 研究區(qū)概況

扎格羅斯盆地A油田位于伊拉克庫爾德地區(qū),構造上處于扎格羅斯盆地疊瓦狀褶皺帶與簡單褶皺帶之間,是上新世扎格羅斯造山運動形成的北西-南東向展布的傾伏背斜構造,主要斷層與背斜構造軸部平行,與背斜構造和斷裂相伴生的構造裂縫發(fā)育,且以高角度裂縫為主。生產層段為白堊系碳酸鹽巖裂縫-孔隙型儲層(圖1),厚約500m,自上而下分成3組:①Shiranish組,為淺灰色-深棕色含有孔蟲的灰?guī)r,含有一定量的泥質,巖性致密,裂縫發(fā)育;②Kometan組,是由厚約100m的灰白色-淺灰色泥晶灰?guī)r組成,巖性較純,泥質含量少;③Qamchuga組,由多孔白云巖段、致密白云巖與灰?guī)r互層段組成。微裂縫系統(tǒng)及大型裂縫通道在3個組都有發(fā)育,部分井在白堊系進行酸洗處理后原油測試產量(1～3)×104桶/d。但不同層系儲層類型、裂縫發(fā)育程度存在的差異造成油井產量差異較大,這是由于該油田裂縫-孔隙型儲層發(fā)育分布的非均質性和隨機性、裂縫發(fā)育控制因素的多樣性造成裂縫系統(tǒng)較為復雜,裂縫的識別、預測和定量表征存在較大難度。

圖1 白堊系地層巖性剖面

研究區(qū)各井均有試油產量資料和壓力數(shù)據(jù),測井資料齊全。常規(guī)測井為標準測井系列,包括自然伽馬、井徑、光電吸收截面指數(shù)、雙側向電阻率、微球形聚焦、聲波時差、巖性密度和中子密度,多數(shù)井有微電阻率成像測井。所使用的地震數(shù)據(jù)為全方位角采集,面元大小為25m×25m,滿覆蓋次數(shù)為80次以上,采集系統(tǒng)的排列橫縱比為0.9,2013年完成保幅疊前時間偏移處理，這為裂縫綜合預測及驗證提供了數(shù)據(jù)基礎。

2 巖石物理特征測試分析

2.1 巖相特征

經過薄片觀察和測試分析可知,白堊系Shiranish組為較深水環(huán)境下含浮游有孔蟲的灰?guī)r,上段(S1—S4)平均泥質含量約為10%,下段(S5—S7)為含泥灰?guī)r,平均泥質含量15%,主要為伊利石/云母,少量高嶺石和綠泥石。Kometan組(細分為K1和K2)為較純的具生物潛穴、擾動構造的含?；?guī)r和粒質灰?guī)r,在局部井區(qū)為中-粗晶白云巖,熱液改造產物,泥質含量較低,平均為2%。Qamchuqa組主要為白云巖組成,微晶白云石交代原始泥晶灰?guī)r，常見腹足類和雙殼類生物碎屑,缺少底棲有孔蟲,說明為相對咸水局限環(huán)境。頂部見厚層巖溶角礫巖,經歷長期大氣水暴露。基質孔隙從致密到孔隙發(fā)育,孔隙度2%～21%,次生鑄?？?、晶間溶孔。局部晶間孔發(fā)育段滲透率較高,高孔低滲段一般以孤立鑄?？缀图毦Щ|孔為主。

2.2 巖石儲集空間類型及特征

結合巖心鑄體薄片觀察和實驗測試分析,歸納出A油田白堊系碳酸鹽巖油藏主要存在兩類儲集空間類型(圖2)。

1) 裂縫。從形態(tài)上分為不連續(xù)微裂縫(裂縫寬度0.05～0.10mm)和連續(xù)性裂縫。不連續(xù)微裂縫只能在巖心和薄片上識別,多數(shù)被方解石膠結；連續(xù)性裂縫分為兩類，一類是早期張性縫,充填方解石和白云石,另一類是晚期構造擠壓形成的與背斜相關的張性裂縫,多為高角度裂縫。裂縫在白堊系Shiranish,Kometan和Qamchuqa組巖心上常有顯現(xiàn)。Shiranish層為裂縫型儲層,裂縫類型豐富,從小尺度的不連續(xù)閉合微縫到大尺度充滿油的開啟大縫均有分布;Kometan層裂縫類型與上部的Shiranish層有明顯差異,主要為中等尺度的裂縫,較少發(fā)現(xiàn)不連續(xù)的微裂縫;Qamchuqa儲層主要為中等規(guī)模的裂縫,并在大斷層或裂縫附近伴生形成塌陷洞穴或溶蝕孔洞，形成的裂縫網絡與較大的斷裂具有較好的連通性,整體屬于裂縫-孔隙型油藏。

圖2 儲集空間類型圖版

2) 鑄?？住⒕чg孔、晶間溶孔。局部晶間孔發(fā)育段滲透率較高,高孔低滲段一般以孤立鑄?？缀图毦Щ|孔為主,孔隙度15%～20%,整體上滲透率不好,范圍(0.01～10.00)×10-3μm2,主要在Qamchuqa組上段發(fā)育。

巖心資料統(tǒng)計表明,白云巖層段的裂縫密度是灰?guī)r裂縫密度的3～10倍,泥質灰?guī)r最差。裂縫以高角度裂縫為主，泥質含量越高,巖石脆性越低,在相同應力條件下發(fā)育的裂縫密度就越低。

3 裂縫的測井信息識別分析

根據(jù)白堊系巖相特征,采用巖心刻度測井資料的方法對白堊系目的層段的礦物成分、泥質含量、基質孔隙度、裂縫產狀、裂縫密度等信息進行分析,總結出各參數(shù)之間的關系,為有利儲集體預測評價提供依據(jù)，主要是通過總結裂縫發(fā)育段在雙側向-微球形聚焦、聲波時差、密度和雙井徑等曲線上的響應特征來識別裂縫發(fā)育段。一般高角度裂縫發(fā)育段呈現(xiàn)出高電阻背景下的中低電阻率異常特征,且深、淺雙側向電阻率為正異常,微球形聚焦曲線顯著降低。密度測井值異常減小也指示裂縫的發(fā)育[2]。當?shù)貙又辛芽p十分發(fā)育時,聲波時差值異常增大。當裂縫寬度較大時,聲波時差曲線會出現(xiàn)周期性跳躍，經巖心檢驗的情況下,聲波時差曲線的小幅度增大也可以作為裂縫的識別標志(圖3)。圖3a為某井裂縫型儲層的測井響應特征,1875～2100m裂縫發(fā)育段在雙側向曲線(RS和RD)上正異?，F(xiàn)象明顯,微球形聚焦曲線(R0)幅值降低,聲波時差(Δt)增大,密度(ρ)降低。

圖3 裂縫的測井信息識別(1ft=0.3048m;1in=2.54cm)

A油田地層微電阻率掃描成像測井資料較多,解釋中剔除誘導縫、紋層、充填縫,張開縫的資料顯示本區(qū)裂縫產狀以NE—SW向為主(平行于區(qū)域主應力方向),其張開縫施密特等面積投影圖顯示多為高角度裂縫,NW—SE方向(平行背斜軸部)的裂縫組不發(fā)育。但不同層系的裂縫密度差異大,Shiranish組裂縫從小尺度的不連續(xù)閉合微縫到大尺度充滿油的開啟縫均有分布,裂縫密度較低;Kometan組發(fā)育中等尺度裂縫,較少發(fā)現(xiàn)不連續(xù)微裂縫,裂縫發(fā)育密度中等;Qamchuqa組構造裂縫發(fā)育(圖4),裂縫網絡具有較好的連通性。成像測井解釋的裂縫產狀近似直立,為高角度縫,裂縫密度2～8條/m,在常規(guī)測井曲線上有高角度縫對應的響應特征。

圖4 巖心統(tǒng)計的裂縫密度

4 地震屬性裂縫預測

研究區(qū)主要采用疊后地震屬性(體曲率、照明3D屬性)和疊前縱波方位各向異性(AVAZ)進行裂縫預測,同時利用已知的地質、測井、試井等先驗裂縫信息對不同地震屬性進行標定,以確定合適的輸入參數(shù),優(yōu)選適合的屬性,使得屬性預測結果與裂縫特征吻合。

4.1 曲率屬性

研究區(qū)裂縫的發(fā)育與上新世以來扎格羅斯造山運動形成的背斜構造密切相關,而地震曲率屬性是利用脆性地層的彎曲程度進行構造解釋和裂縫分析的方法[8],因此可以利用曲率屬性預測研究區(qū)裂縫。為了抑制噪聲對曲率計算的影響,使曲率屬性能更加準確地反映斷裂-微斷裂的空間分布特征,需要首先沿著構造傾角方向對地震資料進行去噪處理,并進行邊緣保存的構造增強濾波[9]。濾波后的地震資料能更好地滿足提取地震幾何屬性的要求。

本區(qū)目的層段為碳酸鹽巖儲層,巖性變化小,裂縫發(fā)育,基本不存在其它地質異常體,小斷層在曲率屬性上表現(xiàn)為線性構造,故曲率屬性能更好地描述垂向上的非連續(xù)性以及與斷裂相關的裂縫通道。這些裂縫通道往往是近于垂直的裂縫面,在橫向上切割儲層,長度一般為幾十米。通過對比用不同Inline/Xline相鄰步長計算出的曲率過井剖面與鉆井試井解釋結果和試油結果,選取合適的步長,利用最大正曲率屬性研究了該區(qū)斷層和裂縫帶的發(fā)育情況。圖5為Kometan向下50m曲率屬性,顯示較為連續(xù)且延伸較長的曲率異常清楚地指示NW—SE向的斷層,這與該區(qū)解釋的斷層在剖面和平面上一致；延伸較短的NE—SW向曲率異常指示裂縫發(fā)育帶。圖6a為NE—SW向裂縫發(fā)育帶上TT-8井白堊系儲層的裂縫發(fā)育程度曲率預測結果；圖6b為TT-8井成像測井解釋結果對比，顯示S1—S4段裂縫密度較小,S5以下至Q6裂縫發(fā)育,且裂縫走向為NE—SW向。試井解釋Shiranish組滲透率947×10-3μm2,Kometan組滲透率2680×10-3μm2,Q1段滲透率30242×10-3μm2,Q2—Q6段滲透率5250×10-3μm2。曲率剖面上顯示的裂縫發(fā)育情況與成像測井解釋結果及試井解釋結果完全吻合。

4.2 縱波方位各向異性裂縫檢測

針對A油田白堊系Shiranish,Kometan和Qamchuqa三套層系巖性橫向變化小、分布均勻、厚度大、裂縫發(fā)育的特點,選用縱波方位各向異性裂縫檢測技術(AVAZ)預測裂縫。首先依據(jù)本區(qū)成像測井解釋的裂縫發(fā)育優(yōu)勢方位NE-SW,平行于背斜軸部的斷層走向NW-SE,在對比不同個數(shù)方位角疊加數(shù)據(jù)的信噪比后,選取有足夠信噪比、方位角盡可能多的8個方位角(5.0°,27.5°,50.0°,72.5°,95.0°,117.5°,140.0°,162.0°)進行疊加,得到用于方位特征分析的方位角疊加數(shù)據(jù),其振幅隨方位角的變化明顯;然后分別將利用3,4,6,8個方位數(shù)據(jù)體的預測結果與地震波形疊合,顯示利用8個方位數(shù)據(jù)體預測的裂縫信息豐富且地震資料品質可靠。預測結果分析表明,縱波方位各向異性預測方法能夠有效識別剖面上同相軸大的錯斷及微弱的不連續(xù)性。試井解釋成果和試油結果很好地驗證了該方法識別裂縫的有效性。圖7a為Kometan以下50m AVAZ預測的裂縫密度平面圖,可以看到TT-13井在裂縫發(fā)育區(qū)邊緣。過TT-13井白堊系AVAZ預測的裂縫密度剖面上Kometan組裂縫不發(fā)育,Qamchuqa組裂縫非常發(fā)育，而TT-13井成像測井解釋結果顯示S3—S4段裂縫不發(fā)育,與AVAZ預測結果略有差異,其它層段對應較好(圖7b,圖7c)。試井解釋Shiranish下段滲透率為13.8×10-3μm2,Kometan組滲透率為495.0×10-3μm2,Qamchuqa組上段滲透率為2940×10-3μm2,與AVAZ和成像測井所反映的裂縫發(fā)育情況基本吻合。

圖5 Kometan組向下50m曲率屬性

圖6 TT-8井白堊系裂縫發(fā)育程度曲率預測結果(a)與成像測井解釋結果對比(b)

圖7 縱波方位各向異性裂縫檢測

4.3 照明3D屬性

利用IHS Kingdom (2015版)提供的基于照明技術的對稱性(Symmetry)和照明3D屬性(Illuminator 3D,I3D),可以進行裂縫帶識別。對稱性用來度量地震數(shù)據(jù)的雜亂程度,使用高階統(tǒng)計量方法分析地震信號樣點分布概率的不對稱程度可以作為識別斷層、裂縫、河道及其它地質特征的工具。對稱性、相似性、曲率等地震屬性可以用來表征斷層,然而由于地震資料品質和自身算法的局限,計算的結果不具有很好的連續(xù)性,與地質規(guī)律有差異。此外,采集腳印的存在會進一步降低這些地表屬性的分辨能力。I3D屬性是在對稱性基礎上,將原始屬性中不連續(xù)、不規(guī)則的斷層面進行組合和連接,使斷層面更符合地質規(guī)律,同時可以消除采集腳印的影響,提高分辨能力。圖8a顯示了I3D處理結果與原始地震剖面解釋的斷層,可見兩者吻合較好。圖8b為平面上沿Kometan層以下0～50m I3D屬性,與解釋斷層一致。本區(qū)Shiranish,Kometan和Qamchuqa三套層系為海相碳酸鹽巖沉積,巖性橫向變化小,排除了河道砂體等其它地質體的可能,因此I3D屬性指示裂縫帶的發(fā)育。

4.4 多屬性融合技術

用單一地震屬性預測裂縫儲集體往往具有多解性,該問題可以通過多屬性融合技術加以改善。融合技術是將多種屬性在一定的數(shù)學運算的基礎上,考慮每一種屬性對于儲層的響應特征及相關程度,在地震屬性歸一化的情況下得出最優(yōu)的裂縫預測結果[10]。通過4.1—4.3節(jié)分析可知,I3D屬性對于識別大的斷層較為可靠,最大曲率屬性可以識別出大—中級別的斷裂,疊前方位各向異性在識別中等級別的裂縫帶時有優(yōu)勢。本文采用線性加權多屬性融合技術將這3種屬性進行融合,如圖9所示。首先計算各屬性的正態(tài)分布范圍,對各屬性做歸一化處理,使每個屬性值范圍處于一個量綱內;然后取各井點周圍一定范圍內(不能太大或太小)屬性的平均值,計算所有井點處屬性的總平均值,將二者之差除以總平均值作為能量偏差(取絕對值)；根據(jù)能量偏差與加權系數(shù)成反比,得到加權系數(shù)k1,k2…,kn(k1+k2+…+kn=1);最后將加權后的能量偏差之和為最小的一組加權系數(shù)作為多屬性融合的加權系數(shù)。由于3種屬性相互取長補短,融合結果與實際情況吻合較好(圖10)。圖10a為過TT-14井的融合體剖面,顯示S上段(S1—S4)裂縫不發(fā)育,成像測井圖顯示該段裂縫密度低,對應的試油產量也不高,1000～2000桶/d。TT-14井融合體剖面還顯示S4以下至Kometan組位于裂縫帶上,成像測井解釋該段裂縫密度高,試油產量可達7000～16000桶/d(圖10)。可見三種信息之間吻合較好。

圖8 照明3D屬性分析

圖9 多屬性融合過程

圖10 過TT-14井多屬性融合結果與實際情況對比

5 動態(tài)驗證及綜合預測

5.1 預測效果動態(tài)驗證

試井分析是通過生產動態(tài)測試來研究油藏參數(shù)、生產能力的一種測試方法,反映較厚地層的宏觀生產能力。A油田采用封隔器逐層試油,試油井段一般為30～100m,因此,用試井解釋成果和試油結果來驗證地震裂縫預測的結果是可行的。從Kometan組以下0～50m范圍內的地震屬性融合圖上(圖11)看,多數(shù)鉆井位于預測的NW—SE和NE—SW向裂縫發(fā)育帶上,試井解釋的滲透率很高(>3000×10-3μm2),在Kometan層產量高(超過10000桶/d)。經過多年生產,近井地帶地層壓力有小幅度下降,說明油井連通的裂縫網絡在平面和縱向分布范圍廣。TT-06井位于相對弱的裂縫發(fā)育帶上,試井解釋的滲透率和產量低一些,但該井仍靠近NW—SE向斷層,產量在2000桶/d以上。NW—SE斷裂平行于背斜軸部,是開啟性裂縫發(fā)育帶,規(guī)模大,地震屬性上響應特征明顯(圖11);NE—SW是平行于現(xiàn)今主應力方向的開啟性裂縫帶,規(guī)模小,但發(fā)育密度大,成像測井解釋出的裂縫主要為該組裂縫,地震屬性上也有響應(圖12)。這兩組裂縫彼此交織,形成裂縫網絡。TT-06井干擾測試結果表明,盡管TT-06井與TT-04井相距4.2km,但TT-06井與TT-04井及TT-05井之間存在橫向連通,白堊系Shiranish,Kometan,Qamchuqa三個地層之間存在垂向連通(干擾由Shiranish地層至Qamchuqa地層)。TT-24井是針對Shiranish組裂縫部署的水平井,已完鉆,在Shiranish組的水平段為1000m,獲得高產。

5.2 有利目標預測及井位建議

通過對預測成果的驗證分析,結合3套目的層系的儲集體特征及本區(qū)油水界面的情況,指出研究區(qū)裂縫發(fā)育的有利目標和建議井位。Kometan和Shiranish組基質致密,裂縫發(fā)育,裂縫產狀以NE-SW向為主,且油水界面已進入Kometan組,構造低部位的油水界面已進入Shiranish組,因此建議在Kometan和Shiranish組的融合地震屬性圖所指示的裂縫發(fā)育區(qū)中部署新的開發(fā)井位及井軌跡,對應圖12上YT30和YT31井。在Shiranish組中段(對應S4段)采用NW—SE向的水平井軌跡,可以鉆遇更多的裂縫帶。

圖11 Kometan組以下0～50m融合屬性

圖12 S4段以下60m融合屬性(TT-23,TT-24為已鉆水平井；YT-30，YT-31為建議水平井井位及軌跡)

6 結論

針對扎格羅斯盆地A油田白堊系巨厚裂縫-孔隙型碳酸鹽巖油藏的地質特點及開發(fā)需求,綜合應用曲率屬性、I3D屬性、方位各向異性及屬性融合等預測方法開展了儲集體綜合預測研究,得出以下結論:

1) 基于不連續(xù)性檢測的疊后屬性對大尺度斷裂系統(tǒng)的刻畫效果顯著,對斷裂輪廓之外的裂縫信息識別能力相對弱,適用于斷層的自動解釋及斷裂帶的識別;

2) 疊前方位各向異性包含地震資料的方位、偏移距信息,對裂縫帶等小尺度斷裂刻畫效果顯著,預測結果分辨率更高,適用于有效刻畫小尺度斷裂及裂縫帶;

3) 多屬性融合技術可綜合各屬性的優(yōu)點,實現(xiàn)從微觀尺度到宏觀尺度的裂縫檢測,提高裂縫型儲集體預測的精度。

[1] 袁士義,宋新民,冉啟全.裂縫性油藏開發(fā)技術[M].北京:石油工業(yè)出版社,2004:3-5 Yuan S Y,Song X M,Ran Q Q.Development technique on fractured reservoir[M].Beijing:Petroleum Industry Press,2004:3-5

[2] 陳冬,魏修成.塔河地區(qū)碳酸鹽巖裂縫型儲層的測井評價技術[J].石油物探,2010,49(2):147-152 Chen D,Wei X C.The evaluation technique on carbonate fractured reservoir from well logging in Tahe Oilfield[J].Geophysical Prospecting for Petroleum,2010,49(2):147-152

[3] 劉振峰,曲壽利,孫建國,等.地震裂縫預測技術研究進展[J].石油物探,2012,51(2):191-198 Liu Z F,Qu S L,Sun J G,et al.The research progress on seismic prediction of fractured reservoir[J].Geophysical Prospecting for Petroleum,2012,51(2):191-198

[4] 曲壽利,季玉新,王鑫.全方位P波屬性裂縫檢測方法[J].石油地球物理勘探,2001,36(4):390-397 Qu S L,Ji Y X,Wang X.Seismic method for using full-azimuth P-wave attribution to detect fracture[J].Oil Geophysical Prospecting,2001,36(4):390-397

[5] Ganguly N,Dearborn D,Moore M,et al.Application of seismic curvature attributes in the appraisal of the Tishrine-West field,North-East Syria[EB/OL].[2015-04-01].http:∥csegrecorder.com/articles/view/application-of-seismic-curvature-attribute-in-appraisal-of-tishrine-west

[6] 齊晴,曲壽利,李振春,等.也門S2區(qū)塊疊前、疊后裂縫綜合預測方法應用研究[J].石油物探,2014,53(5):603-608 Qi Q,Qu S L,Li Z C,et al.Fracture prediction by using pre-stack and post-stack seismic techniques in S2 block,Yemen[J].Geophysical Prospecting for Petroleum,2014,53(5):603-608

[7] 耿福蘭.南堡2號潛山碳酸鹽巖儲層疊前和疊后裂縫綜合預測[J].石油物探,2013,52(1):89-96 Geng F L.Integrated prediction of fracture in carbonate reservoir using pre- and post-stack seismic attributes in Nanpu No.2 buried hill[J].Geophysical Prospecting for Petroleum,2013,52(1):89-96

[8] 孔選林,唐建明,徐天吉.曲率屬性在川西新場地區(qū)裂縫檢測中的應用[J].石油物探,2011,50(5):517-520 Kong X L,Tang J M,Xu T J.Application of seismic curvature attributes in fracture detection in Xinchang area,Chuanxi[J].Geophysical Prospecting for Petroleum,2011,50(5):517-520

[9] 楊威,賀振華,陳學華.結構方位濾波在體曲率屬性中的應用[J].石油物探,2011,50(1):27-32 Yang W,He Z H,Chen X H.Application of structural orientation filtering in volumetric curvature attributes extraction[J].Geophysical Prospecting for Petroleum,2011,50(1):27-32

[10] 趙虎,尹成,朱仕軍.多屬性融合技術研究[J].勘探地球物理進展,2013,32(2):119-121 Zhao H,Yin C,Zhu S J.Study on multiple seismic attributes fusion technique[J].Progress in Exploration Geophysics,2013,32(2):119-121

(編輯：戴春秋)

Multi-disciplinary integration prediction of fractured carbonate reservoir in A Oilfield,Zagros Basin

Zhang Tao,She Gang,Li Miao,AhmatjanAbudrahman,Ye Shuangjiang

(SinopecPetroleumExploration&ProductionResearchInstitute,Beijing100083,China)

A Oilfield located in Zagros Basin is porous-fractured reservoirs with oil column of more than 500m.Complex fracture distribution and diagenesis result in the difficulty in effective reservoir prediction.In order to solve the problem,an integrated study of geological,well logging,seismic attributes and reservoir engineering information on predicting the naturally fractured carbonate reservoir was carried out.The development characteristics of fractures and their controlling factors & development pattern can be identified from the core,casting thin sections,the measured porosity and permeability data in microscopicially.Meanwhile,the conventional and imaging log are applied to quantitatively interpret the lithology,physical properties and the orientation of fractures to provide geologic and petrophysical model for seismic prediction.According to the results of geological and petrophysical interpretation,seismic attributes sensitive to fractures are selected and verified by results of imaging log interpretation,pressure build-up well test interpretation and oil production rates.Traditional 3D post-stack attribute show an advantage of predicting the faults and large-scale fractured belts,such as maximum curvature and illuminator 3D attribute.P-wave azimuthal anisotropy (AVAZ) method can provide more details about the intensity and distribution of the fractured reservoir.Multiple seismic attributes fusion techniques can be used to integrate the advantages of various attributes and reduce the uncertainty of results.Multi-disciplinary research result shows that fractured reservoir can be accurately predicted and evaluated by different scale.

carbonate,fractured reservoir,multi-disciplinary integration identification,seismic attributes fusion,Zagros Basin

2015-04-2;改回日期：2015-08-13。

張濤(1973—),男,博士,高級工程師,從事油藏地質評價、碳酸鹽巖儲層建模等方面的研究工作。

國家科技重大專項項目(2011ZX05031-002)資助。

P631

A

1000-1441(2015)06-0770-10

10.3969/j.issn.1000-1441.2015.06.016