李 軍，張軍華，譚明友，曲志鵬，崔世凌

(1.中國石油大學(xué)(華東)，山東 青島 266580；2.海洋國家實驗室礦產(chǎn)資源評價與探測技術(shù)功能實驗室，山東 青島 266071；3.中國石化勝利油田分公司，山東 東營 257015)

?

高89-4井區(qū)二氧化碳驅(qū)波及范圍地震預(yù)測

李 軍1，2，張軍華1，2，譚明友3，曲志鵬3，崔世凌3

(1.中國石油大學(xué)(華東)，山東 青島 266580；2.海洋國家實驗室礦產(chǎn)資源評價與探測技術(shù)功能實驗室，山東 青島 266071；3.中國石化勝利油田分公司，山東 東營 257015)

通過分析勝利油田高89-4井區(qū)CO2驅(qū)油前后2期地震資料發(fā)現(xiàn)，地震資料在高頻段的頻譜吸收衰減差異明顯。利用該特點提取單頻、有效帶寬能量及頻譜衰減指數(shù)等屬性，對CO2驅(qū)的波及范圍進行了預(yù)測。研究結(jié)果表明，預(yù)測的波及范圍與實際結(jié)果吻合率達到86.7%。預(yù)測結(jié)果對高89-4井區(qū)的CO2驅(qū)有一定指導(dǎo)作用，相關(guān)處理解釋一體技術(shù)對國內(nèi)外CO2驅(qū)油地震監(jiān)測的理論與實踐有一定的借鑒作用。

高89-4井區(qū)；CO2驅(qū)油；地震監(jiān)測；頻率類屬性；吸收衰減；波及面預(yù)測

0 引 言

自Whorton等[1]于1952年首次取得CO2驅(qū)油專利以來，注CO2提高油氣采收率的工作一直在進行。加拿大Pann西部公司、英國BP公司、法國道達爾公司等有許多成功的應(yīng)用實例[2-4]。中國的大慶油田、吉林油田、中原油田、江蘇油田、勝利油田也進行了工業(yè)試驗[5-8]。CO2驅(qū)油及其地震監(jiān)測方法在發(fā)達國家比較成熟，但國內(nèi)針對該技術(shù)的研究還存在一些問題[9-21]，CO2驅(qū)的地震動態(tài)監(jiān)測，特別是CO2驅(qū)波及范圍的預(yù)測方法研究還比較少。勝利油田高89-4井區(qū)CO2驅(qū)先導(dǎo)試驗是中國目前最大的CO2驅(qū)工業(yè)應(yīng)用實例，該區(qū)位于勝利油田正理莊油田的北部，區(qū)域構(gòu)造上處于博興洼陷金家—正理莊—樊家鼻狀構(gòu)造帶的中部。該區(qū)于1992年進行第1次資料采集，由于油田滾動勘探的需要，2012年進行了重新采集。由于先后2次地震資料的采集條件不同，導(dǎo)致國外常規(guī)使用的地震方法無法在該井區(qū)取得較好的研究效果。對2期資料進行一致性匹配處理后發(fā)現(xiàn)，注氣后高頻段存在比較明顯的衰減現(xiàn)象。通過提取單頻、有效帶寬能量及頻譜衰減指數(shù)等屬性，可較為準(zhǔn)確地預(yù)測出CO2驅(qū)的波及范圍。該研究方法對國內(nèi)外CO2驅(qū)油氣監(jiān)測研究有一定的參考價值。

1 CO2驅(qū)2期資料地震特征差異分析

1.1 采集參數(shù)與觀測系統(tǒng)比較

高89-4井區(qū)三維地震資料面積為100.8 km2，從2007年開始，油田選取沙四下段薄互層儲集層進行CO2驅(qū)油先導(dǎo)性試驗。注采層頂面是沙三段不整合面蓋層，兩側(cè)有大斷層夾持，在采油井兩側(cè)布設(shè)規(guī)則的注氣井網(wǎng)(圖1)。

高89-4井區(qū)2期地震資料采集時間相隔近20 a，采集參數(shù)和觀測系統(tǒng)差別很大(表1)。由表1可知，2期資料在覆蓋次數(shù)、偏移距、方位角、頻率及能量等方面均存在較大差異。需先對2期資料進行互約束反褶積、互約束剩余靜校正、覆蓋次數(shù)退化、道集內(nèi)插等互均化處理，才能獲得一致性較好的地震資料。

1.2 新老資料振幅差異比較

經(jīng)過一致性處理后，根據(jù)圖1所示的注氣井網(wǎng)部署，選取從南到北過G89-9、G89-S3、G89-4及G891-7井儲層頂面的振幅值作為代表來研究2期資料幅值差異性(圖2)。由圖2可知：①新資料頂面能量整體較強，但在高89-4井以南到G891-7井之間，局部區(qū)域能力還有變小的現(xiàn)象；②離南北大斷層較近的兩側(cè)(G891-7井以南近2號斷層與G89-9井以北近1號斷層附近區(qū)域，見圖1)幅值差異很大，雖然波動范圍基本一致，但已不能滿足差異性比較的條件，故不能用其來研究CO2驅(qū)替。

圖1 高89-4井區(qū)CO2驅(qū)油先導(dǎo)性試驗井網(wǎng)

年份觀測系統(tǒng)放炮方式道距/m線距/m炮點/m線距/m最大偏移距/m道數(shù)覆蓋次數(shù)網(wǎng)格尺寸/m×m19924L6S單邊5020015020031502402025×100201218L12S雙邊2520025255500360022525×25

圖2 新老資料儲層頂面振幅值比較

1.3 新老資料頻譜差異比較

分析對比儲層處的頻譜特征(圖3)可知：新資料頻帶范圍稍窄為8.0～43.0 Hz，老資料頻帶范圍為6.0～50.0 Hz；新老資料主頻基本一致，大約在23.4 Hz左右；新資料相對于老資料在25.0～50.0 Hz頻段衰減明顯(圖3綠色虛線標(biāo)注范圍)。

圖3 儲層頻譜特征對比

分析總結(jié)2期資料的地震特征可知：2期資料的振幅差異沒有規(guī)律可循，但2期資料頻譜在25.0～50.0 Hz頻段存在明顯的差異。充分利用該特點，使用單頻屬性、頻率吸收衰減屬性等頻率類屬性提取技術(shù)，從2期地震資料中提取相應(yīng)的屬性數(shù)據(jù)體，而后沿T7層抽取屬性切片進行注氣前后的差異性分析，以期檢測出CO2驅(qū)的波及范圍。

2 頻率類屬性檢測方法

頻率類屬性可以有效利用地震數(shù)據(jù)中蘊含的頻率信息，表征需要的地震特征與差異。單頻屬性展示的是地震數(shù)據(jù)某單一頻率信息，不同的單一頻率屬性可以表征該頻率所包含的特有信息。另外，理論研究和實際應(yīng)用表明，在地質(zhì)體中，如果孔隙發(fā)育，充填油、氣、水(尤其對于含氣的情況)時，地震反射吸收加大，高頻吸收衰減加劇，含油氣地層吸收系數(shù)可比相同巖性不含油氣地層高幾倍甚至一個數(shù)量級。地層含有非飽和的油氣時，能量衰減更明顯地表現(xiàn)出異常衰減，而吸收衰減類屬性利用的就是這種特性，主要包括衰減梯度因子、有效頻帶寬度、有效帶寬能量及頻譜衰減指數(shù)等。根據(jù)上文的分析，在特定頻率段，注氣后的地震資料存在著頻率衰減現(xiàn)象。因此，采用單頻屬性和吸收衰減屬性來進行分析。

2.1 單頻屬性

根據(jù)前文注氣后新資料在25.0～50.0 Hz頻率段存在明顯衰減的認識，運用廣義S變換[22]分別提取2期資料的25.0、45.0 Hz單頻屬性數(shù)據(jù)體，而后沿T7層下延10～75 ms(時窗范圍為沙四下段注氣位置)抽取得到單頻屬性切片(圖4)。由于1、2號大斷層所起的封堵作用，只分析2個斷層之間的區(qū)域部分，根據(jù)頻率衰減效果，結(jié)合圖1所示注氣井網(wǎng)分布，標(biāo)注出利用25.0、45.0 Hz單頻屬性檢測得到的CO2驅(qū)替波及范圍，如圖4中藍色虛線圈定區(qū)域所示。

對比2種單頻屬性切片預(yù)測結(jié)果可知：45.0 Hz單頻屬性切片圖中藍色虛線圈定區(qū)域頻率衰減更為明顯，效果更好，即45.0 Hz單頻屬性的預(yù)測結(jié)果更為準(zhǔn)確。

圖4 單頻屬性沿層切片比較(左圖：老資料，右圖：新資料)

2.2 吸收衰減屬性

編程實現(xiàn)吸收衰減類屬性[23-24]提取算法，對2期地震數(shù)據(jù)體進行計算得到相應(yīng)的注氣前后頻率衰減屬性數(shù)據(jù)體，沿T7層下延10～75 ms提取有效帶寬能量和頻譜衰減指數(shù)屬性(圖5)。分析2個斷層之間的區(qū)域，根據(jù)頻率衰減效果，結(jié)合圖1所示注氣井網(wǎng)分布，標(biāo)注出利用有效帶寬能量、頻譜衰減指數(shù)屬性檢測得到的CO2驅(qū)替波及范圍，如圖5中藍色虛線圈定區(qū)域所示。分析對比2種屬性的預(yù)測結(jié)果可知：2種屬性藍色虛線圈定區(qū)域的衰減效果均很明顯，預(yù)測范圍也基本一致；相比于單頻屬性，頻率衰減屬性預(yù)測結(jié)果的波及范圍更大，方法一致性較好，預(yù)測可靠性更高。

圖5 頻譜衰減類屬性比較(左圖為老資料，右圖為新資料)

根據(jù)工區(qū)構(gòu)造特點及井網(wǎng)分布，實地考察了圖1井網(wǎng)中的34口井。分析以上幾種頻率類屬性的波及范圍預(yù)測結(jié)果可知：單頻屬性的預(yù)測結(jié)果中，只有11口井在預(yù)測的波及范圍內(nèi)，而實地考察結(jié)果顯示，共有15口井出氣，而且單頻屬性預(yù)測結(jié)果與注氣井點的灘壩砂分布規(guī)律并不吻合，效果相對較差；有效帶寬能量屬性的預(yù)測結(jié)果顯示：共有14口井在預(yù)測的波及范圍內(nèi)，其中13口井的結(jié)果與實地考察一致，但在儲層兩端與井點認識存在差異；頻譜衰減指數(shù)屬性預(yù)測準(zhǔn)確率更高，有13口井的預(yù)測結(jié)果與實地考察相吻合，準(zhǔn)確率達到了86.7%。

3 結(jié)論與認識

(1) 注氣前后采集的2期地震資料存在著采集條件不同的問題，經(jīng)過一致性處理后，2期資料時差基本一致且振幅值波動范圍基本一致，但在儲層兩端區(qū)域還存在較大差異，不能利用振幅或時差來分析注氣前后的差異性；2期資料的主頻基本一致，但注氣后發(fā)生了高頻衰減現(xiàn)象，使得2期資料在高頻段存在一些差異，可以用來對CO2驅(qū)波及范圍進行預(yù)測。

(2) 運用單頻、頻譜衰減指數(shù)及有效帶寬能量等多種頻率類屬性進行注氣后波及范圍的預(yù)測，并與實地考察結(jié)果進行了對比。單頻屬性預(yù)測結(jié)果準(zhǔn)確性較低，而有效帶寬能量及頻譜衰減指數(shù)等頻率衰減類屬性的檢測效果較高，特別是頻譜衰減指數(shù)屬性的波及范圍預(yù)測結(jié)果，相比于實地考察得到的波及范圍，其井點吻合率達到了86.7%。預(yù)測結(jié)果對高89-4井區(qū)的CO2驅(qū)有實際指導(dǎo)作用，對CO2驅(qū)油地震監(jiān)測的理論與實踐有一定的借鑒作用。

[1] WHORTON L P， BROWNSCOMBE E R， DYES A B. Method for producing oil by means of carbon dioxide： US， 2623596[P]. 1952.

[2] KHATIWADA M， VAN W K， ADAM L， et al. A numerical sensitivity analysis to monitor CO2sequestration in layered basalt with coda waves [C]. Houston， 2009 SEG Annual Meeting： Society of Exploration Geophysicists， 2009：3865-3869.

[3] GUTIERREZ M， KATSUKI D， ALMRABAT A. Effects of CO2injection on the seismic velocity of sandstone saturated with saline water[J]. International Journal of Geosciences， 2012， 3(5)： 908.

[4] WANG Z， CATES M E， LABGAN R T. Seismic monitoring of a CO2flood in a carbonate reservoir： A rock physics study [J]. Geophysics， 1998， 63(5)： 1604-1617.

[5] 劉炳官，朱平，雍志強，等.江蘇油田CO2混相驅(qū)現(xiàn)場實驗研究[J].石油學(xué)報，2002，23(4)： 56- 60.

[6] 劉必心，侯吉瑞，李本亮，等.二氧化碳驅(qū)特低滲油藏的封竄體系性能評價[J].特種油氣藏，2014，21(3)： 128- 132.

[7] 羅二輝，胡永樂，李保柱，等.中國油氣田注CO2提高采收率實踐[J].特種油氣藏，2013，20(2)： 1-7， 42.

[8] 李峰，楊永超.濮城油田沙一段油藏濮1-1井組CO2驅(qū)研究及先導(dǎo)實驗效果分析[J].海洋石油，2009，29(4)： 56- 60.

[9] 白素. 榆樹林油田CO2驅(qū)非達西滲流數(shù)值模擬研究[D]. 大慶：東北石油大學(xué)，2011.

[10] 張永剛，羅懿，劉岳龍，等.注CO2井筒溫度和壓力分布模型研究及現(xiàn)場應(yīng)用[J].巖性油氣藏，2014，26(2)：108-113.

[11] 李向良.溫度和注入壓力對二氧化碳驅(qū)油效果的影響規(guī)律實驗[J].油氣地質(zhì)與采收率，2015，22(1)：84-87,92.

[12] 陳德春，周淑娟，孟紅霞，等.陳家莊油田陳373塊蒸汽吞吐后轉(zhuǎn)CO2—化學(xué)劑復(fù)合吞吐研究[J].油氣地質(zhì)與采收率，2014，21(6)：76-78.

[13] 尚寶兵，廖新維，趙曉亮，等.雜質(zhì)氣體對二氧化碳驅(qū)最小混相壓力和原油物性的影響[J].油氣地質(zhì)與采收率，2014，21(6)：92-94.

[14] 尚寶兵，廖新維，盧寧，等.CO2驅(qū)水氣交替注采參數(shù)優(yōu)化——以安塞油田王窯區(qū)塊長6油藏為例[J].油氣地質(zhì)與采收率，2014，21(3)：70-72,77.

[15] 李菊，李實.基于多次接觸特性曲線的二氧化碳驅(qū)最小混相壓力預(yù)測方法[J].油氣地質(zhì)與采收率，2014，21(1)：51-54.

[16] 殷丹丹，李宜強，趙東鋒. CO2驅(qū)產(chǎn)出氣回注的驅(qū)油效果[J].大慶石油地質(zhì)與開發(fā)，2014，33 (2)： 137-140.

[17] 汪艷勇. 大慶榆樹林油田扶楊油層CO2驅(qū)油試驗[J]. 大慶石油地質(zhì)與開發(fā)，2015，34 (1)：136-139.

[18] 鐘立國，馬帥，魯淵，等. 稠油與CO2、CH4或N2體系高溫高壓界面張力測定分析[J]. 大慶石油地質(zhì)與開發(fā)，2015，34 (1)：140-145.

[19] 謝宗林. VLEFlash相態(tài)模型在昌德氣田CO2氣藏的應(yīng)用[J]. 大慶石油地質(zhì)與開發(fā)，2015，34 (3)：76-80.

[20] 李紹杰. 低滲透灘壩砂油藏CO2近混相驅(qū)生產(chǎn)特征及氣竄規(guī)律[J]. 大慶石油地質(zhì)與開發(fā)，2016，35 (2)：110-115.

[21] 王建波，高云叢，宗暢，等. 特低滲油藏CO2非混相驅(qū)水氣交替注入見效特征[J]. 大慶石油地質(zhì)與開發(fā)，2016，35 (2)：116-120.

[22] STOCKWELL R G， et al. Localization of the complex spectrum： the S-transform [J]. IEEE Transactions on Signal Processing， 1996， 44(4)：998-1001.

[23] MITCHELL J T， DERZHI N， LICHMA E. Energy absorption analysis： a case study [C]. Denver， Expanded Abstracts of 66th SEG Mtg， 1996： 1785-1788.

[24] 刁瑞，李振春，韓文功，等.基于廣義S變換的吸收衰減分析技術(shù)在油氣識別中的應(yīng)用[J].石油物探，2011，50(3)：260-265.

編輯 朱雅楠

20160405；改回日期：20160820

中國石油化工股份有限公司先導(dǎo)課題“CO2驅(qū)油地震監(jiān)測評價方法探索研究”(P14085)

李軍(1988-)，男，2010年畢業(yè)于中國石油大學(xué)(華東)信息與計算科學(xué)專業(yè)，現(xiàn)為該校地質(zhì)資源與地質(zhì)工程專業(yè)在讀博士研究生，主要從事地震解釋與儲層預(yù)測工作。

10.3969/j.issn.1006-6535.2016.06.009

TE122

A

1006-6535(2016)06-0040-05