李悅+張紹良

doi：10.15889/j.issn.1002-1302.2017.13.063[HT9.]

摘要：以某CO2驅(qū)油技術(shù)（CO2-EOR）示范工程為研究對(duì)象，通過(guò)對(duì)CO2-EOR示范區(qū)內(nèi)玉米長(zhǎng)勢(shì)的監(jiān)測(cè)，采用對(duì)比分析、顯著性檢驗(yàn)等方法，定量分析了CO2-EOR對(duì)玉米生長(zhǎng)的影響。結(jié)果表明，和對(duì)照區(qū)相比，監(jiān)測(cè)區(qū)內(nèi)玉米的指標(biāo)均值波動(dòng)為0.23%～5.39%，但長(zhǎng)勢(shì)沒(méi)有顯著差異；油井附近的玉米長(zhǎng)勢(shì)總體優(yōu)于其他區(qū)域，指標(biāo)均值的波動(dòng)幅度在0.33%～3.00%，但是也沒(méi)有顯著差異。究其原因，CO2-EOR過(guò)程中從抽油井中攜帶出來(lái)的CO2輕微提高了油井近地表的CO2濃度，產(chǎn)生施肥效應(yīng)。這一結(jié)果顯示，該CO2-EOR示范區(qū)內(nèi)玉米沒(méi)有受到CO2溢出的負(fù)面影響，反而在一定程度上促進(jìn)了農(nóng)作物的生長(zhǎng)。

關(guān)鍵詞：CO2地質(zhì)封存；CO2驅(qū)油技術(shù)（CO2-EOR）；玉米；影響分析

中圖分類號(hào)： S162.5文獻(xiàn)標(biāo)志碼： A[HK]

文章編號(hào)：1002-1302（2017）13-0238-04[HS）][HT9.SS]

收稿日期：2016-12-29

基金項(xiàng)目：國(guó)家科技支撐計(jì)劃（編號(hào)：2012BAC24B05）

作者簡(jiǎn)介：李悅（1993—），女，四川廣元人，碩士研究生，從事生態(tài)環(huán)境監(jiān)測(cè)研究。E-mail：1658826806@qq.com。

通信作者：張紹良，博士，教授，從事生態(tài)監(jiān)測(cè)、CO2監(jiān)測(cè)研究。Tel：（0516）83591301；E-mail：cumtzslly@126.com。

[ZK）]

近年來(lái)，人們特別關(guān)注CO2人為排放引起的全球氣候變暖問(wèn)題，將CO2埋存于地下是解決該問(wèn)題的有效措施[1]。CO2驅(qū)油技術(shù)（CO2-EOR，CO2-enhanced oil recovery）是CO2捕獲與存儲(chǔ)優(yōu)先發(fā)展的技術(shù)，也最具商業(yè)推廣前景[2-3]。但是CO2在運(yùn)輸、注入地下、驅(qū)油以及封存過(guò)程中，可能會(huì)發(fā)生泄漏，而不同途徑的泄漏對(duì)植被、生態(tài)系統(tǒng)的影響至今還不清楚[4-5]。在中國(guó)，很多油井都分布在平原農(nóng)業(yè)生產(chǎn)地區(qū)，因此人們擔(dān)心CO2-EOR工程會(huì)影響農(nóng)作物的生長(zhǎng)。

國(guó)內(nèi)外學(xué)者研究發(fā)現(xiàn)，地質(zhì)封存CO2泄漏導(dǎo)致的大氣CO2濃度升高，對(duì)農(nóng)作物的影響主要表現(xiàn)在促進(jìn)植物光合作用、增加其生物量累積；顯著提高C3作物產(chǎn)量，但對(duì)C4作物產(chǎn)量的影響很小；降低了C3和C4作物氣孔導(dǎo)度，非常顯著地提高了所有作物的水分利用率等方面[6-9]。但也有研究指出，在淹水、地面板結(jié)導(dǎo)致的土壤通氣不暢情景下，水稻、豆類植物和柑橘等作物并未受根部CO2濃度升高、O2缺乏的影響，玉米受高濃度土壤CO2短期影響后能夠得到恢復(fù)[10-11]。因此，地下CO2泄漏對(duì)不同農(nóng)作物的影響機(jī)理仍存在爭(zhēng)論，未能得出明確的結(jié)論。此外，CO2泄露對(duì)生態(tài)環(huán)境影響的研究方法主要停留在實(shí)驗(yàn)室模擬、開放系統(tǒng)中的模擬以及數(shù)值模型模擬等方面[4，6，10-11]，都是在假設(shè)一定濃度CO2泄漏的基礎(chǔ)上，分析其對(duì)植被、環(huán)境的影響，而缺少對(duì)實(shí)際CO2封存工程的監(jiān)測(cè)與研究。

為此，本試驗(yàn)以一個(gè)CO2-EOR工程項(xiàng)目為例，通過(guò)對(duì)監(jiān)測(cè)區(qū)、對(duì)照區(qū)、油井附近區(qū)和油井遠(yuǎn)離區(qū)的玉米作物長(zhǎng)勢(shì)進(jìn)行監(jiān)測(cè)，調(diào)查CO2-EOR工程實(shí)施對(duì)玉米生長(zhǎng)的影響，以此探討CO2-EOR技術(shù)的安全性

1材料與方法

1.1研究區(qū)概況

[JP2]CO2-EOR示范區(qū)位于山東省某油田范圍內(nèi)（圖1），屬黃河下游的魯北平原，北溫帶季風(fēng)大陸性氣候，春秋短而冬夏長(zhǎng)，[JP]光照充足，年平均氣溫13.10 ℃，年平均降水量 539.4 mm，[JP]年均無(wú)霜期198 d，夏季多雨，冬春多旱。油田位于青城凸起東緣和東營(yíng)凹陷的西南緣，主要油氣層埋藏集中在900～1 500 m，埋深大于1 500 m的中深層油藏主要分布于東區(qū)斷裂帶，探明的含油面積7.9 km2，石油地質(zhì)儲(chǔ)量0.196億t。區(qū)域內(nèi)主要作物是玉米和小麥，農(nóng)地水利灌溉設(shè)施良好。

1.2研究方法

1.2.1對(duì)比區(qū)域設(shè)置

本研究設(shè)置了監(jiān)測(cè)區(qū)與對(duì)照區(qū)（圖2），監(jiān)測(cè)區(qū)位于CO2-EOR工程實(shí)施區(qū)域內(nèi)，面積約 50 km2。工程采取一注多采的方式，核心注采區(qū)位于監(jiān)測(cè)區(qū)中心位置，50余口采出井則分布在監(jiān)測(cè)區(qū)各處。監(jiān)測(cè)區(qū)主要分布著耕地和村莊居民點(diǎn)；對(duì)照區(qū)則位于工程實(shí)施區(qū)2 km以外的西南角，氣候條件、耕作方式、灌溉系統(tǒng)及作物種植類型與監(jiān)測(cè)區(qū)相同。

為了監(jiān)測(cè)油井附近玉米作物的生長(zhǎng)狀況，劃分油井附近區(qū)和遠(yuǎn)離油井區(qū)。根據(jù)實(shí)地CO2濃度監(jiān)測(cè)結(jié)果，采出井口有較高濃度（6 000 μL/L）（圖3）的CO2逸散，但CO2擴(kuò)散很快，距離采出井50 m處沒(méi)有影響。因此，以50 m作為油井附近區(qū)和遠(yuǎn)離油井區(qū)劃分界限。

1.2.2樣本調(diào)查與數(shù)據(jù)采集

采用隨機(jī)抽樣方法，在對(duì)照區(qū)范圍內(nèi)隨機(jī)設(shè)置5個(gè)固定樣地，在監(jiān)測(cè)區(qū)范圍內(nèi)隨機(jī)設(shè)置18個(gè)[CM（25]固定樣地。在每一個(gè)樣地，布設(shè)3～5個(gè)樣方，每個(gè)樣方內(nèi)

采集玉米作物數(shù)據(jù)。連續(xù)2年監(jiān)測(cè)，分別于2013、2014年8月20日實(shí)地調(diào)查采樣。

葉面積指數(shù)、株高是對(duì)CO2濃度變化敏感的玉米作物長(zhǎng)勢(shì)指標(biāo)[12-13]。葉面積指數(shù)是反映作物群體大小的較好的宏觀動(dòng)態(tài)指標(biāo)，在一定的范圍內(nèi)，作物的產(chǎn)量隨葉面積指數(shù)的增大而提高，CO2濃度增加使植被葉面積指數(shù)增加，且促進(jìn)作用比較明顯[12-13]。株高是一項(xiàng)可反映植被形態(tài)、健康狀況植被生理指標(biāo)，大氣CO2濃度升高能顯著增加作物株高，促使作物生長(zhǎng)，提高產(chǎn)量[13]。采用冠層分析儀測(cè)量樣方內(nèi)玉米葉面積指數(shù)，軟尺測(cè)量株高。

1.2.3數(shù)據(jù)處理與計(jì)算

采用單因素方差分析方法，檢驗(yàn)監(jiān)測(cè)區(qū)和對(duì)照區(qū)、油井附近區(qū)和遠(yuǎn)離區(qū)的玉米葉面積指數(shù)、株高的差異顯著性，檢驗(yàn)量F值為：endprint

[JZ（]F=[SX（]SSA/（k-1）SSE/（n-1）[SX）]；[JZ）][JY]（1）

[JZ（]SSA=∑[DD（]ki=1[DD）]ni（x[TX-*5]i-x[TX-*5]）2；[JZ）][JY]（2）

[JZ（]SSE=∑[DD（]ki=1[DD）] ∑[DD（]nii=1[DD）]（xij-x[TX-*5]i）2。[JZ）][JY]（3）

式中：i表示控制變量的第i水平，即第i總體；k表示總體數(shù)；n表示樣本總數(shù)；xij表示第i總體的第j觀測(cè)值；x[TX-*5]i表示第i總體的樣本均值；x[TX-*5]表示樣本總均值；SSA為組間平方和，SSE為組內(nèi)平方和。

采用Excel 2010軟件統(tǒng)計(jì)數(shù)據(jù)，用SPSS 19.0軟件完成單因素方差分析。

2結(jié)果與分析

2.1監(jiān)測(cè)區(qū)和對(duì)照區(qū)內(nèi)玉米長(zhǎng)勢(shì)的對(duì)比

2.1.1描述統(tǒng)計(jì)量分析

從表1可知，2013年對(duì)照區(qū)玉米葉面積指數(shù)、株高均值稍大于監(jiān)測(cè)區(qū)均值；2014年則相反，監(jiān)測(cè)區(qū)玉米葉面積指數(shù)和株高均值略大于對(duì)照區(qū)均值；2年的綜合統(tǒng)計(jì)結(jié)果則是對(duì)照區(qū)葉面積指數(shù)均值稍大于監(jiān)測(cè)區(qū)均值，而監(jiān)測(cè)區(qū)株高均值稍大于對(duì)照區(qū)均值。以對(duì)照區(qū)為基準(zhǔn)，監(jiān)測(cè)區(qū)內(nèi)玉米長(zhǎng)勢(shì)指標(biāo)均值的波動(dòng)幅度在0.23%～5.39%。2年玉米長(zhǎng)勢(shì)指標(biāo)中值的比較結(jié)果和均值比較結(jié)果相同?？梢姡瑢?duì)照區(qū)和監(jiān)測(cè)區(qū)的玉米長(zhǎng)勢(shì)不存在明顯差異。

2.1.2單因素方差分析

由檢驗(yàn)結(jié)果（表2）可以看出，葉面積指數(shù)的F值為0.125，株高的F值為0.000，對(duì)照區(qū)域內(nèi)玉米的葉面積指數(shù)方差沒(méi)有顯著差異（P>0.05），株高方差也沒(méi)有顯著差異（P>0.05）。進(jìn)一步對(duì)2013年和2014年的玉米長(zhǎng)勢(shì)指標(biāo)進(jìn)行檢驗(yàn)（表3）可知，2年間監(jiān)測(cè)區(qū)和對(duì)照區(qū)內(nèi)玉米葉面積指數(shù)、株高的方差也不存在顯著差異。

方差檢驗(yàn)結(jié)果表明，葉面積指數(shù)和株高作為對(duì)CO2濃度變化敏感的玉米作物長(zhǎng)勢(shì)指標(biāo)，在CO2-EOR區(qū)域和對(duì)照區(qū)域中都沒(méi)有顯著差異，表明玉米的生長(zhǎng)沒(méi)有受到顯著影響。

2.2油井附近區(qū)與遠(yuǎn)離油井區(qū)內(nèi)玉米長(zhǎng)勢(shì)對(duì)比

2.2.1描述統(tǒng)計(jì)分析

由表4可以發(fā)現(xiàn)，2013年油井附近區(qū)玉米葉面積指數(shù)、株高的均值稍大于遠(yuǎn)離油井區(qū)的均值；2014年玉米2項(xiàng)指標(biāo)的均值也是油井附近區(qū)均值大于遠(yuǎn)離油井區(qū)均值。與遠(yuǎn)離油井區(qū)相比，油井附近區(qū)內(nèi)玉米長(zhǎng)勢(shì)指標(biāo)均值的波動(dòng)幅度在0.32%～2.92%。中值的比較結(jié)果和均值的比較結(jié)果相同，即油井附近區(qū)葉面積指數(shù)、株高中值略大于遠(yuǎn)離油井區(qū)中值。

2.2.2單因素方差分析

從表5可知，葉面積指數(shù)的F值為0.038，株高的F值為1.144，對(duì)應(yīng)的P值分別是0.846和0292，均大于0.05的顯著性水平，表明沒(méi)有顯著差異。同樣地，分別對(duì)2013、2014年對(duì)比區(qū)域內(nèi)玉米葉面積指數(shù)和株高進(jìn)行F檢驗(yàn)，結(jié)果如表6所示。由表6可以看出，在油井附近區(qū)和遠(yuǎn)離油井區(qū)內(nèi)，玉米2項(xiàng)長(zhǎng)勢(shì)指標(biāo)沒(méi)有顯著差異。

綜上所述，盡管表4的統(tǒng)計(jì)結(jié)果顯示，油井附近的玉米葉面積指數(shù)和株高的均值稍大于遠(yuǎn)離油井區(qū)的均值，但是表5和表6的檢驗(yàn)結(jié)果表明這種差異并不顯著（P>0.05），可見靠近油井的玉米作物長(zhǎng)勢(shì)比遠(yuǎn)離油井的玉米長(zhǎng)勢(shì)稍好，但不存在顯著差異。

3討論

現(xiàn)有研究表明，植被短期暴露于高濃度大氣CO2條件下，對(duì)其生長(zhǎng)具有肥效作用，長(zhǎng)期暴露情景下，植被生理、形態(tài)特征和基因均發(fā)生適應(yīng)性地變化[10，14-15]。CO2濃度升高對(duì)農(nóng)作物的影響作用主要表現(xiàn)在促進(jìn)植物光合作用、增加其生物量累積，顯著提高C3作物產(chǎn)量，但對(duì)玉米等C4作物產(chǎn)量的影響很小[6]。CO2驅(qū)油過(guò)程中，從油井逃逸出來(lái)的CO2導(dǎo)致油井附近地表CO2濃度增加，對(duì)農(nóng)作物產(chǎn)量稍有促進(jìn)作用，且影響范圍與風(fēng)向有關(guān)[16]。據(jù)研究區(qū)內(nèi)CO2濃度監(jiān)測(cè)，距離油井越近，近地表的CO2濃度越高，最大可達(dá)6 000 μL/L，但隨著距離的增加，地表CO2濃度迅速減小直至恢復(fù)一般水平（圖3）。從表4可以發(fā)現(xiàn)，油井附近區(qū)內(nèi)玉米的葉面積指數(shù)和株高均值都比遠(yuǎn)離油井區(qū)的均值稍高，而且連續(xù)2年監(jiān)測(cè)的結(jié)果均相同，這可能就是因?yàn)橛途椴蓵r(shí)CO2逸散出來(lái)，增高了油井附近的地表CO2濃度，促進(jìn)了這些區(qū)域內(nèi)的玉米生長(zhǎng)，但是表5和表6結(jié)果顯示這種間歇性的逸散對(duì)農(nóng)作物長(zhǎng)勢(shì)促進(jìn)作用不大，不存在顯著差異。

地質(zhì)封存的CO2泄漏路徑主要有火山等劇烈地質(zhì)活動(dòng)、構(gòu)造斷層或破碎帶、高滲透性的滲流帶、井口等[10，17-19]。研究區(qū)內(nèi)可能存在的泄漏途徑有2種：一是CO2在驅(qū)油過(guò)程中隨著石油到達(dá)井口從而逃逸。從油井附近的CO2濃度監(jiān)測(cè)結(jié)果來(lái)看，抽油的同時(shí)確實(shí)存在CO2的逸散情況，但是逸散范圍小，并對(duì)作物生長(zhǎng)有輕微的促進(jìn)作用。另一種途徑是從地質(zhì)斷層泄漏至地表。但本研究區(qū)內(nèi)沒(méi)有發(fā)現(xiàn)此類泄漏，原因可能是，該區(qū)域油藏深度達(dá)900～1 500 m，上面蓋層厚且完整性好，而且CO2從2009年開始注入，注入量?jī)H4萬(wàn)t/年，注入量也不大。從監(jiān)測(cè)區(qū)和對(duì)照區(qū)玉米長(zhǎng)勢(shì)的比較結(jié)果可知，CO2-EOR區(qū)域玉米長(zhǎng)勢(shì)沒(méi)有顯著變化，這說(shuō)明示范區(qū)實(shí)施CO2-EOR工程對(duì)農(nóng)作物生長(zhǎng)沒(méi)有威脅。

本研究是在連續(xù)注入5年后開展的監(jiān)測(cè)，且連續(xù)監(jiān)測(cè)了2年。統(tǒng)計(jì)分析表明，示范工程區(qū)域內(nèi)農(nóng)作物玉米沒(méi)有受到影響，不過(guò)本研究的的結(jié)論是建立在幾個(gè)主要生理指標(biāo)以及實(shí)地抽樣調(diào)查的數(shù)據(jù)的基礎(chǔ)上的，因此，只能表明在本研究區(qū)農(nóng)作物玉米種植區(qū)內(nèi)開展CO2-EOR工程是安全的。不過(guò)，這一結(jié)論在一定程度上消除了人們對(duì)CO2-EOR工程實(shí)施對(duì)地表農(nóng)作物生產(chǎn)的擔(dān)憂。

4結(jié)論

本試驗(yàn)以CO2-EOR示范區(qū)域玉米作物為對(duì)象，以葉面積指數(shù)、株高為指標(biāo)進(jìn)行研究，對(duì)比分析了監(jiān)測(cè)區(qū)和對(duì)照區(qū)、油井附近區(qū)和遠(yuǎn)離油井區(qū)兩組區(qū)域的玉米長(zhǎng)勢(shì)情況。結(jié)果表明，和對(duì)照區(qū)相比，監(jiān)測(cè)區(qū)內(nèi)玉米指標(biāo)均值波動(dòng)為0.23%～[JP2]539%，但方差分析表明玉米長(zhǎng)勢(shì)沒(méi)有顯著差異；油井附近的[JP2]玉米長(zhǎng)勢(shì)總體優(yōu)于其他區(qū)域，[JP3]指標(biāo)均值的波動(dòng)幅度在0.33%～[JP]3.00%，但是F檢驗(yàn)結(jié)果發(fā)現(xiàn)也沒(méi)有顯著差異；分析其原因，CO2-EOR過(guò)程中從抽油井中攜帶出來(lái)的CO2輕微提高了油井近地表的CO2大氣濃度，產(chǎn)生了施肥效應(yīng)。這一結(jié)果顯示，該CO2-EOR示范區(qū)玉米作物沒(méi)有受到CO2溢出的負(fù)面影響，反而在一定程度上促進(jìn)了農(nóng)作物的生長(zhǎng)。endprint

由于受到研究區(qū)本身?xiàng)l件、監(jiān)測(cè)指標(biāo)、調(diào)查方法等限制，只能證明本研究區(qū)內(nèi)開展CO2-EOR工程是安全的，在一定程度上消除了人們對(duì)CO2地質(zhì)封存的擔(dān)憂。

[HTK]致謝：感謝項(xiàng)目組楊永均、劉莘、徐雅晴、黃安平、郝紹金、任雪峰等對(duì)資料的收集和整理，為本研究提供了大量分析數(shù)據(jù)，感謝油田公司和工作人員對(duì)我們?cè)谡{(diào)研過(guò)程中提供的幫助。[HT]

[HS2]參考文獻(xiàn)：[HJ1.65mm]

[1][ZK（#]張亮，王舒，張莉，等. 勝利油田老油區(qū)CO2提高原油采收率及其地質(zhì)埋存潛力評(píng)估[J]. 石油勘探與開發(fā)，2009，36（6）：737-742.

[2]劉瑜. 二氧化碳地下封存與強(qiáng)化采油利用基礎(chǔ)研究[D]. 大連：大連理工大學(xué)，2011.

[3]江懷友，沈平平，盧穎，等. CO2提高世界油氣資源采收率現(xiàn)狀研究[J]. 特種油氣藏，2010，17（2）：5-10.

[4][JP2]聶莉娟，馬俊杰，趙雪峰，等. 模擬CCS技術(shù)CO2泄漏對(duì)C3、[JP3]C4作物土壤化學(xué)性質(zhì)的影響[J]. 水土保持學(xué)報(bào)，2015，29（5）：200-205.[JP]

[5]張森琦，刁玉杰，程旭學(xué)，等. 二氧化碳地質(zhì)儲(chǔ)存逃逸通道及環(huán)境監(jiān)測(cè)研究[J]. 冰川凍土，2010（6）：1251-1261.

[6]Kimball B A，朱建國(guó)，程磊，等. 開放系統(tǒng)中農(nóng)作物對(duì)空氣CO2濃度增加的響應(yīng)[J]. 應(yīng)用生態(tài)學(xué)報(bào)，2002，13（10）：1323-1338.

[7]龔偉，宮淵波，胡庭興，等. CO2濃度升高對(duì)濕地松針葉蒸騰特性和水分利用效率的影響[J]. 水土保持學(xué)報(bào)，2005，19（5）：178-182.

[8]李中陽(yáng)，樊向陽(yáng)，宋正國(guó)，等. CO2濃度升高對(duì)不同品種水稻Cu、Cd吸收及根際土Cu，Cd形態(tài)的影響[J]. 水土保持學(xué)報(bào)，2013，27（6）：78-82.

[9]郭建平，高素華. 植物對(duì)高CO2和土壤干旱的響應(yīng)趨勢(shì)及對(duì)策[J]. 水土保持學(xué)報(bào)，2004，18（2）：170-173.[ZK）]

[10][ZK（#]田地，馬欣，查良松，等. 地質(zhì)封存CO2泄漏對(duì)近地表陸地生態(tài)系統(tǒng)的影響綜述[J]. 生態(tài)與農(nóng)村環(huán)境學(xué)報(bào)，2013，29（2）：137-145.

[11]聶莉娟，馬俊杰，趙雪峰，等. 模擬CCS技術(shù)CO2泄漏對(duì)C3、C4作物土壤化學(xué)性質(zhì)的影響[J]. 水土保持學(xué)報(bào)，2015，29（5）：200-205.

[12]袁東敏，郭建平. CO2濃度增加對(duì)東北玉米生長(zhǎng)影響的數(shù)值模擬[J]. 自然資源學(xué)報(bào)，2010，25（5）：822-829.

[13]王修蘭，徐師華. CO2濃度增加對(duì)C3、C4作物生育和產(chǎn)量影響的實(shí)驗(yàn)研究[J]. 中國(guó)農(nóng)業(yè)科學(xué)，1998，31（1）：55-61.

[14]Anderson T H，Heinemeyer O，Weigel H J. Changes in the fungal-to-bacterial respiratory ratio and microbial biomass in agriculturally managed soils under free-air CO2 enrichment （FACE）—a six-year survey of a field study[J]. Soil Biology and Biochemistry，2011，43（5）：895-904.

[15]Nakamura I，Onoda Y，Matsushima N，et al. Phenotypic and genetic differences in a perennial herb across a natural gradient of CO2 concentration[J]. Oecologia，2011，165（3）：809-818.

[16]Xu Y Q，Zhang S L，Hou H P，et al. Influence of CO2 leakage from oil-producing wells on crop growth based on improved CASA model[JP3][J]. International Journal of Remote Sensing，2016，37（2）：290-308.[JP]

[17]Hill P M. Possible asphyxiation from Carbon dioxide of a cross-country skier in eastern California：a deadly volcanic hazard[J]. Wilderness & Environmental Medicine，2000，11（3）：192-195.

[18]Giammanco S，Gurrieri S，Valenza M. Fault-controlled soil CO2 degassing and shallow magma bodies：Summit and lower East Rift of Kilauea Volcano （Hawaii），1997[J]. Pure and Applied Geophysics，2006，163（4）：853-867.

[19]吳江莉，馬俊杰. 淺議CO2地質(zhì)封存的潛在風(fēng)險(xiǎn)[J]. 環(huán)境科學(xué)導(dǎo)刊，2012，31（6）：89-93.endprint