時(shí)移地震技術(shù)在CO2地質(zhì)儲存監(jiān)測中的適用性

由榮軍， 閆照濤

(1.中國礦業(yè)大學(xué)　銀川學(xué)院 礦業(yè)工程系， 銀川　750011；2. 寧夏地球物理地球化學(xué)勘查院， 銀川　750001)

?

時(shí)移地震技術(shù)在CO2地質(zhì)儲存監(jiān)測中的適用性

由榮軍1， 閆照濤2

(1.中國礦業(yè)大學(xué)銀川學(xué)院 礦業(yè)工程系， 銀川750011；2. 寧夏地球物理地球化學(xué)勘查院， 銀川750001)

摘要：目前，國內(nèi)、外正在興起對CO2地中隔離技術(shù)的研究和試驗(yàn)工作。為了解決地質(zhì)儲存之后，CO2氣體的流動(dòng)性及安全性監(jiān)測問題，結(jié)合我國實(shí)際地質(zhì)地層情況，同時(shí)借鑒了國內(nèi)外對CO2深部地下儲層隔離的方法和研究思路，這里采用正演模擬的方法，模擬CO2在地下的運(yùn)移擴(kuò)散情況，并通過編程模擬時(shí)移地震技術(shù)的處理過程對獲得的常規(guī)地震剖面進(jìn)行時(shí)移地震處理。從獲得的常規(guī)地震剖面和時(shí)移地震剖面的結(jié)果顯示，這里提出的時(shí)移地震技術(shù)在CO2地下儲存檢測中具有一定的適用性。

關(guān)鍵詞：CO2； 地質(zhì)儲存； 時(shí)移地震； 巖石物理學(xué)； 正演模擬

0引言

目前，雖然對于全球氣候異常的原因還沒有定論，世界各國對于節(jié)能減排，減少大氣中CO2含量的努力卻在不斷加強(qiáng)。就目前各國進(jìn)行的實(shí)際CO2儲集項(xiàng)目來看，砂巖咸水層和深部煤層成為儲集CO2的首選場所。對于CO2儲存到地下之后如何流動(dòng)，引起地層的巖性性質(zhì)怎么變化，以及是否存在泄漏風(fēng)險(xiǎn)，這些問題目前仍在實(shí)驗(yàn)探索中。要解決這些問題需要進(jìn)行實(shí)際的監(jiān)測實(shí)驗(yàn)來確定。這里的重點(diǎn)就是結(jié)合前人研究結(jié)果，綜合分析地層巖石物理學(xué)屬性，運(yùn)用正演模擬的方法探討時(shí)移地震技術(shù)在CO2地質(zhì)儲存監(jiān)測中的適用性。

1巖石物理學(xué)和流體替代理論

彈性波在固體中的傳播速度依賴于介質(zhì)的彈性變形。在地球中，孔隙巖石的彈性變形和礦物成分、壓實(shí)過程、孔隙度、膠結(jié)程度、壓力、溫度和孔隙中流體占有率有關(guān)。運(yùn)用巖石物理學(xué)從已知的巖石特性中預(yù)測多孔巖石的彈性特性。

1.1Gassman理論

對巖石物理學(xué)的研究理論有很多，其中應(yīng)用最普遍也最適用的理論為Gassman理論。Gassman理論的基本假設(shè)為觀測巖石(骨架和基質(zhì))宏觀上是均勻分布的，所有的孔隙是連通的，孔隙中充滿無摩擦流體(液體、氣體或者混合物)，研究的巖石-流體系統(tǒng)是封閉的(不排出液體)，孔隙流體不會(huì)對巖石骨架產(chǎn)生軟化或硬化作用。

經(jīng)過理論研究，Gassman理論的主體方程表示如式(1)所示。

(1)

其中：φ為孔隙度；Ksat是飽和巖石的體積模量；Kdry是干燥巖石的體積模量；Kgrain是礦物顆粒的體積模量；Kfluid是飽和流體的體積模量。

同時(shí)，Gassman理論假設(shè)在彈性波傳播期間孔隙中飽和的流體不與骨架產(chǎn)生任何相互作用。由于牛頓流體的剪切模量為“0”，則飽和巖石的剪切模量等于干燥巖石的剪切模量(在低頻情況下)如式(2)所示。

μsat= μdry

(2)

Gassman方程中涉及到的各種體積模量和剪切模量都可以經(jīng)過理論實(shí)驗(yàn)估計(jì)出來。

1.2煤層中流體替代理論

煤層中的孔隙可分為原生孔隙和次生孔隙，此外煤層中還含有大量的裂隙，這使得煤層形成雙重孔隙系統(tǒng)類型。裂隙在煤中的作用相當(dāng)重要，它不僅可以儲存注入到煤層的氣體，同時(shí)又可以將基質(zhì)間孔隙連通起來，增強(qiáng)煤層的滲透性。正是由于這種雙重孔隙系統(tǒng)類型，使得煤層具有很大的內(nèi)表面積，這成為氣體在煤層中聚集的有利條件。

煤層中原有的氣體主要有CH4、CO2和N2，根據(jù)Langmuir等溫吸附模型(圖1)，當(dāng)CO2注入到煤層中后，由于競爭吸附，CO2會(huì)替代原有的CH4和N2儲存在煤層中。

圖1　55℃時(shí)單組分等溫吸附線Fig．1　Single component adsorption 　　isotherm at 55℃

當(dāng)CO2注入地下后，根據(jù)壓力和溫度條件的不同，它可以以氣體、液體、固體或超臨界流體的狀態(tài)存在(圖2)。

圖2　CO2相態(tài)圖Fig．2　CO2 phase diagram

1.3砂巖層中流體圈閉理論

當(dāng)CO2被注入到砂巖層中后，由于砂巖層的巖性特征及其所含咸水的一些特性，使得CO2在砂巖層中發(fā)生各種圈閉作用，包括構(gòu)造圈閉作用、溶解度圈閉作用和礦物圈閉作用。這些圈閉作用使得砂巖層能儲存更多的CO2，同時(shí)由于CO2的注入也會(huì)引起砂巖層巖性性質(zhì)的變化，包括孔隙度變化、孔隙壓變化以及流體替代等引起的變化，這一變化最終會(huì)對地震波勘探速度產(chǎn)生影響。

2CO2儲集模擬

根據(jù)前面的理論分析，這里分別對CO2在砂巖層和煤層中的儲集情況進(jìn)行了模擬，并獲得理論地震波數(shù)據(jù)，用于隨后的分析研究。

根據(jù)Gassman理論和儲層流體替代理論，CO2的注入會(huì)引起儲層孔隙度、孔隙壓及氣體飽和度的變化，從而引起地震波速度的變化，經(jīng)理論分析并編程計(jì)算，得出速度變化分別與儲層孔隙度、孔隙壓、氣體飽和度變化之間的關(guān)系(圖3)

圖3　速度變化曲線Fig．3　The carve of velocity variation(a)速度隨孔隙度變化曲線;(b)速度隨孔隙壓變化曲線;(c)速度隨二氧化碳飽和度變化曲線

2.1砂巖層CO2儲集模擬

模型1為CO2注入前初始模型,模型中，砂泥巖互層總厚度為300m，其中砂巖和泥巖各自厚度分別為20 m，以互層形式交替沉疊，孔隙介質(zhì)為水，沒有氣體注入，模型參數(shù)如表3所示，地質(zhì)結(jié)構(gòu)模型示意圖如圖4所示，礦巖層基本地震參數(shù)如表1所示。

表1　砂巖層基本地震參數(shù)

表2　CO2儲集數(shù)值模擬觀測系統(tǒng)參數(shù)表

表3　砂巖層初始模型參數(shù)

表4　煤層基本地震參數(shù)

表5　煤層初始模型參數(shù)

圖4　CO2注入前初始砂巖層模型地質(zhì)結(jié)構(gòu)Fig．4　Geologic framework for initial model before injecting CO2 into sandstone

根據(jù)表2給出的觀測系統(tǒng)參數(shù)，采用聲波方程進(jìn)行數(shù)值模擬，得到了36炮的單炮記錄，運(yùn)用處理解釋軟件對單炮記錄進(jìn)行常規(guī)處理：直達(dá)波切除、抽道集、速度分析、動(dòng)校正、疊加和偏移等，得到了這一初始模型的地震剖面，如圖5所示。

從圖5以看出，當(dāng)?shù)貙铀剑瑳]有注入CO2氣體時(shí)，我們通過地震記錄得到的剖面同相軸是水平的，能較好地映地下巖層基本形態(tài)。因此，這一記錄剖面可以作為監(jiān)測不同時(shí)期注入氣體流動(dòng)狀態(tài)的基線剖面，用于與注入氣體剖面進(jìn)行對比。

模型2為CO2流動(dòng)范圍200 m模型。在該模型中，砂泥巖互層總厚度為300 m，其中砂巖和泥巖各自厚度分別為20 m，以互層形式交替沉疊，孔隙介質(zhì)為水，CO2氣體流動(dòng)范圍為200 m，地質(zhì)結(jié)構(gòu)模型示意圖如圖6所示。

圖6　CO2注入模型地震結(jié)構(gòu)圖Fig．6　Geologic framework for CO2injection mode

經(jīng)過與上述模型相同的處理，得到了這一模型的地震剖面，如圖7所示。與圖5相比，當(dāng)?shù)貙铀?，注入一定范圍的CO2氣體時(shí)，我們通過地震記錄得到的剖面同相軸在氣體注入?yún)^(qū)域是彎曲的，時(shí)間軸出現(xiàn)下拉現(xiàn)象，而且在氣體飽和度較大的中間部位出現(xiàn)煙囪狀的形狀。這些異常形態(tài)較真實(shí)地反映了圖6正演模型中的氣體注入情況。通過對模擬地震資料進(jìn)行簡單的處理分析，得出對比結(jié)果，直觀地說明運(yùn)用地震勘探技術(shù)能較真實(shí)的監(jiān)測CO2流體在砂巖層中的儲存、流動(dòng)情況。

圖7　CO2注入模型地震剖面Fig．7　Seismic profile for CO2 injection model

2.2煤層CO2儲集模擬

作者根據(jù)CO2不同階段的擴(kuò)散范圍，建立了4個(gè)模型。

在模型中，煤層厚度為10 m，頂板厚度為220 m，底板厚度為370 m，孔隙介質(zhì)為水，氣體注入范圍分別為0 m、20 m、100 m、200 m，具體的模型參數(shù)如表5所示，地質(zhì)結(jié)構(gòu)模型示意圖如圖8所示。

圖8　CO2不同注入范圍煤層模型地質(zhì)結(jié)構(gòu)及地震剖面Fig．8　Geologic framework and seismic profile for model of CO2 into coal(a)氣體流動(dòng)范圍0 m；(b)氣體流動(dòng)范圍20 m；(c)氣體流動(dòng)范圍100 m；(d)氣體流動(dòng)范圍200 m

從圖8可以看出：當(dāng)煤層水平， CO2氣體流動(dòng)范圍逐漸從0 m 大到200 m時(shí)，得到的地震記錄剖面上其異常范圍擴(kuò)大、同相軸下拉、能量增強(qiáng)，在注入點(diǎn)兩邊對稱位置還出現(xiàn)類似煙囪狀的同相軸上凹現(xiàn)象。

3屬性分析與時(shí)移地震處理

3.1屬性分析

前面對氣體不同流動(dòng)范圍進(jìn)行模擬，通過分析各個(gè)模型的地震波同相軸的差異，定性地描述了由于氣體注入引起的地震波異常。為了進(jìn)一步說明異常情況，作者將通過屬性分析的方法進(jìn)行論述，并重點(diǎn)分析了不同范圍模型振幅屬性的差異。通過振幅屬性分析，我們得到的各個(gè)模型的能量對比如圖9所示。

從圖9可以看出，沒有氣體注入時(shí)，地震波在煤層中的能量是均勻分布的；當(dāng)煤層中存在氣體時(shí)，能量大多集中在氣體聚集區(qū)域，并且其較強(qiáng)能量集中范圍能很好的與氣體流動(dòng)范圍相匹配。

圖9　煤層中不同氣體流動(dòng)范圍振幅屬性圖Fig．9　Amplitude attribute map for different range of gas flow in coal(a)為初始水平模型振幅屬性圖；(b)為氣體流動(dòng)范圍20 m的振幅屬性圖；(c)為氣體流動(dòng)范圍100 m的振幅屬性圖；(d)為氣體流動(dòng)范圍200 m的振幅屬性圖

3.2移地震處理

到目前為止，地震解釋已從提供精確的構(gòu)造圖發(fā)展到進(jìn)行精細(xì)儲層描述的階段。將地震數(shù)據(jù)與測井資料、巖心及其巖層物理數(shù)據(jù)相結(jié)合，可以實(shí)現(xiàn)儲層靜態(tài)精細(xì)成像。實(shí)際上，儲集層除了包含靜態(tài)信息外，還包括動(dòng)態(tài)信息，如流體飽和度、壓力和溫度等，這些對了解儲層的動(dòng)態(tài)變化大有益處。但是要在一次測量中同時(shí)得到這兩方面的信息是很難的，這就必須借助于多次地震觀測——時(shí)間推移地震。

時(shí)間推移地震的發(fā)展歷史，可以追溯到20多年前開始的蒸汽驅(qū)和注氣過程，主要用于研究油藏動(dòng)態(tài)變化，這里論述的氣體流動(dòng)動(dòng)態(tài)變化與油藏動(dòng)態(tài)變化有類似之處，因此可用此技術(shù)來監(jiān)測所研究的CO2動(dòng)態(tài)流動(dòng)過程。

這里正演模擬中均采用一致的采集、處理參數(shù)，為保證數(shù)據(jù)匹配，仍然進(jìn)行了一些互均化處理。這里互均化按時(shí)間校正、均方根振幅校正、帶寬互均化和相位匹配來描述[8-12]，其數(shù)學(xué)模型如式(3)所示。

SXEQ= S*f(tcorr,rmscorr,mcorr,pcorr)

(3)

其中：S為地震道；*表示褶積；f代表脈沖響應(yīng);影響因素有tcorr、rmscorr、mcorr、pcorr分別代表時(shí)間校正、RMS振幅校正、帶寬互均化和相位匹配項(xiàng)?；ゾ院?，在沒有儲層流體的地方，基礎(chǔ)勘探和監(jiān)測勘探應(yīng)該是一致的。作者編程實(shí)現(xiàn)互均化處理，對模擬取得的地震資料處理后得到的差異剖面見圖10。

從圖10中可以看出，沒有氣體注入時(shí)，差異為0；隨著氣體流動(dòng)范圍的不斷擴(kuò)大，差異范圍也相應(yīng)變大。該結(jié)果初步證明這里均化處理程序的準(zhǔn)確性，也一定程度上說明采用時(shí)移地震方法監(jiān)測氣體在儲層中動(dòng)態(tài)變化的可行性。

圖10　煤層中不同氣體流動(dòng)范圍時(shí)移剖面圖Fig．10　Time-lapse cross-sectional view for different range of gas flow in coal(a) 氣體流動(dòng)范圍0 m ；(b) 氣體流動(dòng)范圍20 m；(c) 氣體流動(dòng)范圍100 m；(d) 氣體流動(dòng)范圍200 m

4結(jié)論

作者以CO2地質(zhì)儲存為背景，通過正演模擬描述儲層注入CO2后其巖性變化以及引起的地震波差異，分析處理合成地震記錄，動(dòng)態(tài)監(jiān)測氣體流動(dòng)情況。對正演得到的地震資料進(jìn)行簡單的處理分析，通過同相軸對比、屬性分析對比及時(shí)移地震差異剖面對比，證明地震勘探技術(shù)，尤其是時(shí)移地震技術(shù)在監(jiān)測CO2儲存動(dòng)態(tài)變化方面的可行性及準(zhǔn)確性。

參考文獻(xiàn)：

[1]劉小燕.二氧化碳深含水層隔離的二相滲流模擬與巖石物理學(xué)研究[D].北京：中國科學(xué)技術(shù)大學(xué)，2005.

LIU X Y.Two-phase flow simulation and petrophysical study of carbon dioxide deep aquifers isolated[D].Beijing:China University of Science and Technology,2005.(In Chinese)

[2]MARCIA L.Couёslan.Processing and interpretation of time-lapse vertical seismic profile data from the penn west CO2monitoring project[D]. University of Calgary ,2007.

[3]林凱. 基于Gassman方程的彈性參數(shù)數(shù)值模擬研究[D].成都:成都理工大學(xué)，2010.

LIN K.Numerical Simulation of equations based on the elastic parameters Gassman[D].Chengdu:Chengdu University of Technology,2010.(In Chinese)

[4]張璐.基于巖石物理的地震儲層預(yù)測方法應(yīng)用研究[D].青島：中國石油大學(xué)，2009

ZHANG L.Application seismic reservoir prediction method based on petrophysical[D].Qingdao:China University of Petroleum,2009. (In Chinese)

[5]JERRY M. HARRIS .Reservoir modeling, rock physics, and seismic simulation of CO2sequestration in coal beds[R] .Stanford University，2003.

[6]BOB HARDAGE, DIANA SAVA .Seismic Model for Monitoring CO2Sequestration[J] .Search and Discovery Article,2009(05):130-136.

[7]崔永君.煤對CH4,N2, CO2及多組分氣體吸附的研究[D].西安:煤炭科學(xué)研究總院西安分院，2003.

CUI Y J.Coal of CH4, N2, CO2 and multicomponent gas adsorption[D].Xi’an:China Coal Research Institute Xi’an Branch,2003.(In Chinese)

[8]史亮.數(shù)值模擬和四維地震在巖性油藏開發(fā)中的協(xié)同應(yīng)用[D].成都：西南石油大學(xué)，2009.

SHI L.Collaborative applications about numerical simulation and 4D seismic in the olithologic reservoir development[D].Chengdu:Southwest Petroleum University, 2009. (In Chinese)

[9]李德棟.氣-水-鹽-礦體系相平衡耦合化學(xué)平衡及其在二氧化碳地質(zhì)儲存數(shù)值模擬中的應(yīng)用 [D].北京：中科院研究生院，2008.

LI D D.Gas-water-salt-mineral equilibrium system and its application in chemical equilibrium coupling numerical simulation of geological storage of carbon dioxide[D].Beijing:Graduate School of Chinese Academy of Sciences, 2008.(In Chinese)

[10]YINBUIN LIU,DON C. LAWTON .Well-seismic ban-dwidth and time-lapse seismic characterization related to CO2injection and fluid substitution: physical considerations[J].Well-seismic bandwidth and time-lapse seismic,2006(04):100-110.

[11]ANDY CHADWICK,GARETH WILLIAMS.Quantitative analysis of time-lapse seismic monitoring at the Sleipner CO2storage operation[J].SEG,2010,29(2):170-177.(In Chinese)

[12]陳小宏，牟永光.四維地震油藏監(jiān)測技術(shù)及其應(yīng)用[J].石油地球物理勘探，1998,(06):707-715.

CHEN X,MOU Y G.4D seismic reservoir monitoring technology and its application[J].Oil geophysical prospecting,1998(06):707-715.(In Chinese)

收稿日期：2015-04-08改回日期：2015-07-23

作者簡介：由榮軍(1986-)，女，碩士，主要從事礦業(yè)工程教學(xué)與研究工作及CO2地質(zhì)儲存方面的地震勘探方法研究，E-mail：745620803@qq.com。

文章編號：1001-1749(2016)03-0375-07

中圖分類號：P 631.4

文獻(xiàn)標(biāo)志碼：A

DOI：10.3969/j.issn.1001-1749.2016.03.13

The applicability of time-lapse seismic monitoring in the geological storage of CO2

YOU Rong-jun1, YAN Zhao-tao2

(1.Department of Mineral Engineering,China University of Mining and Technology YinChuan college,Yinchuan750011,China;2.Ningxia Academy of Geophysical and Geochemical Exploration,Yinchuan750001,China)

Abstract:Recently, a lot about sequestration of CO2 in subsurface have been studied in both china and other countries. More and more people are concerned about security of CO2 flowing and diffusing in subsurface. Combined with our actual geologic and stratigraphic situation, and also benefit from our and abroad studying in sequestration of CO2 by injection into deep reservoir, this paper mainly study the applicability of time-lapse seismic prospecting in CO2 flowing and diffusing in subsurface.

Key words:CO2; geology storing; time-lapse seismic prospecting; rock physics; forward modeling