深部咸水層CO2地質(zhì)封存數(shù)值模擬參數(shù)的全局敏感性分析
——以蘇北盆地鹽城組為例

鄭 菲，施小清，吳吉春，趙 良，陳 旸

南京大學(xué)地球科學(xué)與工程學(xué)院，南京 210093

深部咸水層CO2地質(zhì)封存數(shù)值模擬參數(shù)的全局敏感性分析
——以蘇北盆地鹽城組為例

鄭 菲，施小清，吳吉春，趙 良，陳 旸

南京大學(xué)地球科學(xué)與工程學(xué)院，南京 210093

基于對蘇北盆地鹽城組下段砂巖儲層概化建立二維徑向剖面模型，運用TOUGH2/ECO2N程序模擬深部咸水層中CO2遷移分布過程。采用定性(Morris法)和定量(Sobol和EFAST法)全局敏感性分析方法，以儲層中注入井處的壓力、氣相CO2總量及CO2氣相羽擴散距離作為模型響應(yīng)變量，對儲層的水平滲透率、孔隙度、殘余液體飽和度、孔隙分布指數(shù)、壓縮系數(shù)、進氣壓力的倒數(shù)和鹽度7個參數(shù)進行全局敏感性分析，討論了咸水層的水文地質(zhì)參數(shù)、鹽度及其他模型參數(shù)對CO2封存運移泄漏過程的影響。Morris和Sobol法的敏感性分析結(jié)果都表明，對于不同的響應(yīng)變量，參數(shù)的敏感性排序不同：以注入井處的壓力為響應(yīng)變量時，孔隙度的敏感性最高；以氣相CO2總量和CO2羽擴散距離為響應(yīng)變量時，水平滲透率的敏感性最高。EFAST 1階及總敏感度分析結(jié)果與Sobol分析結(jié)果基本一致，但EFAST法相對Sobol法計算更高效穩(wěn)健，需要的樣本數(shù)較少。

CO2地質(zhì)封存；全局敏感性分析；Morris法；Sobol法；EFAST法;蘇北盆地

0 前言

隨著社會經(jīng)濟的發(fā)展、人類活動的增加及化石能源的大量利用，大氣中以CO2為主的溫室氣體的濃度持續(xù)增加，導(dǎo)致全球氣候變暖和海洋酸化，對生態(tài)系統(tǒng)和人類的生存環(huán)境產(chǎn)生嚴重影響[1]。由于咸水層廣泛分布于世界上的沉積盆地中且咸水不能作為飲用水源等原因，CO2深部咸水層封存被認為是目前減少CO2排放最具發(fā)展前景的方法之一。我國華北平原大部，四川盆地、準噶爾盆地東南部等許多盆地都是將來優(yōu)先考慮的CO2封存場地，封存潛力巨大[2]。其中蘇北盆地是近海陸相沉積盆地，面積約3.5×104km2，蘊含豐富的油氣資源，具有咸水層封存CO2的巨大潛力。由于人類生產(chǎn)活動對巖層完整性的破壞，不密封的斷層、天然裂縫等的存在[3-4]，注入到咸水層的CO2將在浮力作用下通過滲漏通道向上部逃逸，導(dǎo)致CO2封存存在一定的泄漏風(fēng)險。為確保注入的CO2能安全、有效、長久地封存在地下咸水層中，有必要對影響CO2泄露的因素進行敏感性分析，從而確定影響CO2封存泄露的主要因子。

目前國內(nèi)外已有多位學(xué)者應(yīng)用TOHGH2-ECO2N/TOUGHREACT[5-6]開展了深部咸水層CO2地質(zhì)封存影響因素的研究。Zhou等[7]以Illinois盆地西蒙山含水層為例，建立了盆地尺度和CO2羽尺度下的整合模型，表明低滲透率和高進氣壓力具有二次密封效應(yīng)，可阻滯CO2向上遷移，而注入井附近壓力增量會促使儲層中CO2-咸水向上運移進入上覆含水層；Brirkholzer等[8]選取蓋層滲透率及壓縮系數(shù)，進行了大尺度CO2封存區(qū)域壓力增量的敏感性分析，表明這2個參數(shù)的變化對儲層壓力增量的影響均較大；Wiese等[9]開展了影響CO2單井注入速率的參數(shù)局部敏感性分析，表明滲透率的影響最大，注入壓力的影響次之，垂向各向異性的影響較小；李義連等[10]分析了江漢盆地高鹽度鹵水對CO2封存過程中所發(fā)生的物理化學(xué)變化的影響，表明高鹽度鹵水層中單純地注入CO2會造成注入井附近巖鹽大量沉淀，不利于CO2安全儲存及持續(xù)注入；柯怡兵等[11]開展了江漢盆地咸水層CO2封存注入數(shù)值模擬研究，對鹽度的敏感性分析表明儲層鹽度越高壓力積累越嚴重；趙銳銳等[12]進行了松遼盆地深部咸水層CO2注入的流體遷移模擬研究及參數(shù)敏感性分析，表明殘余氣體飽和度、注入層的各向異性數(shù)值對CO2溶解捕獲的影響最大，鹽度次之，殘余液體飽和度和初始CO2濃度的影響最小。這些成果為研究咸水層CO2地質(zhì)封存影響因素提供了定性認識，但已有研究往往僅局限于單個參數(shù)因素進行情景或局部敏感性分析。

全局敏感性分析方法目前主要包括Morris篩選法[13]、傅里葉幅度敏感度檢驗法(FAST/EFAST)[14-15]以及基于方差分析的Sobol法[16]等。近年來的研究成果[17-19]表明，全局敏感性分析可同時分析多個參數(shù)變化及參數(shù)間相互耦合作用對模型輸出結(jié)果的直接和間接影響，較局部敏感性分析(只能分析單個參數(shù)對模型結(jié)果的直接影響)更接近實際，可全面地檢測輸入?yún)?shù)對模擬結(jié)果的影響，有助于模型使用者更好地理解識別模型結(jié)構(gòu)和解釋模擬結(jié)果。目前關(guān)于全局敏感性分析方法已廣泛應(yīng)用于地理空間分析模型[20]、氣候模型[21]、水質(zhì)模型[22]等模型敏感性分析中，但針對CO2地質(zhì)封存模型進行全局敏感性分析的研究甚少。

筆者以蘇北盆地鹽城組下段砂層為研究對象，建立簡單的二維徑向模型，應(yīng)用TOUGH2/ECO2N[5]模擬注入的CO2在地質(zhì)儲層中的溶解遷移以及CO2-咸水驅(qū)替過程，重點運用定性(Morris法)和定量(Sobol和EFAST法)全局敏感性分析方法針對不同的響應(yīng)變量，對選取的模型參數(shù)進行全局敏感性分析，探討和分析各參數(shù)的變化對深部咸水層CO2地質(zhì)封存泄露的直接和間接影響。

1 CO2地質(zhì)封存的數(shù)值模型

1.1 研究區(qū)地質(zhì)概況

蘇北盆地內(nèi)部發(fā)育一系列北東向的南陡北緩的斷陷和凹陷，是油氣聚集的良好場所[23-24]，也為CO2地質(zhì)封存提供了良好的儲存空間。鹽城組是泛濫平原河流相沉積，以砂礫巖-中粗砂巖、細砂巖為主，屬高孔隙高滲透性儲集巖層。鹽一段(下段)(Ny1)由3個沉積旋回組成，為粉砂質(zhì)泥巖與中粗砂巖、細礫巖互層，頂部埋深一般在1 000.0 m以上，均厚300.0～400.0 m，含水層中的鹽度約為1.5%，含水層孔隙度范圍為12%～25%，滲透率范圍為50.0～800.0 mD*mD為非法定計量單位， 1 mD=0.986 9×10-15 m2，下同。[25]；鹽二段(上段)(Ny2)中、上部是泥巖發(fā)育段，一般單層厚度為20.0～30.0 m，可作為有利的阻擋蓋層[26]。該組地層基本符合CO2地質(zhì)儲存的選址條件，從儲層埋深、巖層特征、區(qū)域蓋層封堵性、地質(zhì)圈閉條件等因素判斷，可作為一個潛在的CO2地質(zhì)封存場所。

1.2 概念模型建立

圖1 一維CO2注入概念模型示意圖(a)及網(wǎng)格剖分圖(b)Fig.1 Schematic diagram for the 1D injection model (a) and the mesh (b)

筆者選取鹽城組下段砂巖地層，建立二維徑向模型，模擬CO2注入到均質(zhì)各向異性咸水層中的運移演化過程(圖1)。模型水平半徑R為2.5 km，共劃分為52個單元，網(wǎng)格間距隨離井口距離增大而增大，最外側(cè)單元設(shè)為無流量邊界；垂直方向厚度H

為100.0 m，頂部及底部網(wǎng)格厚度均為2.0 m，設(shè)為無流量邊界，其余網(wǎng)格厚度均為6.0 m，共18層。注入井半徑設(shè)為0.5 m，CO2注入點位于距底部17.0 m處，恒定注入速率Qinj為0.1 kg/s，注入時間為10 a，總模擬時間為100 a?？紤]到后續(xù)敏感性分析需要大量的計算時間，本文模擬范圍較小但網(wǎng)格仍精細剖分，注入速率及模擬時間均設(shè)為較小值，以節(jié)省計算時耗同時避免影響到邊界。鹽城組下段砂層頂部埋深為1 000.0 m，模型頂部靜水壓強p約為10.0 MPa，能使注入的CO2以超臨界狀態(tài)存在。系統(tǒng)溫度T為45 ℃，巖石顆粒密度、熱傳導(dǎo)率及比熱分別為2 600 kg/m3、2.51 W/(m·℃)、920 J/(kg·℃)。相對滲透率及毛管壓力的概化分別采用Van Genuchten-Mualem模型和Van Genuchten模型，輸入?yún)?shù)詳見表1[25]。

1.3 模擬結(jié)果分析

注入階段，由于氣相CO2(超臨界狀態(tài)，以下統(tǒng)稱為氣相)密度比水小，注入的CO2在浮力作用下向上遷移，孔隙中氣相CO2飽和度逐漸增大，運移至蓋層底部時由于受到阻擋，開始在蓋層底部聚集并沿底部側(cè)向擴散遷移(圖2a)；整個注入階段，注入井的壓力增量最大，10 a末注入井壓力為1.06×107Pa(圖2b)；停止注入后，CO2繼續(xù)向上并沿蓋層底部側(cè)向遷移(圖2c)，100 a后氣相羽的遷移距離為335.4 m。CO2飽和度逐漸降低，最大飽和度由0.52減小為0.21，表明越來越多的CO2溶解在儲層咸水中，儲層壓力增量得以緩慢釋放。10 a后CO2

表1 模型中的水文地質(zhì)及其他參數(shù)

圖中θliq為氣相CO2飽和度。圖2 10 a后氣相CO2飽和度(a)、儲層壓力(b)以及100 a后氣相CO2飽和度(c)的分布情況Fig.2 CO2 gas saturation (a),Spatial distribution of pressure (b) after 10 years and CO2 gas saturation after 100 years (c)

注入總量為3.2×107kg，氣相及溶解相CO2各占總注入量的60%和40%；100 a后氣相CO2溶解于咸水中使得溶解相CO2含量增加，各占總注入量的38%和62%。

CO2注入到儲層后，由于壓力積累可能誘發(fā)巖石破裂，破壞地層的密封性及完整性，導(dǎo)致儲層中的氣相CO2在浮力的作用下持續(xù)向上部逃逸，發(fā)生CO2泄露，其可能的泄漏范圍與CO2驅(qū)替區(qū)域大小密切相關(guān)。為了評估模型參數(shù)對CO2封存泄漏的影響，筆者選取注入結(jié)束時(10 a后)注入井處的壓力(此時儲層壓力增量最大)、氣相CO2總量及模擬結(jié)束時(100 a后)CO2氣相羽沿蓋層底部的擴散距離作為下文敏感性分析的響應(yīng)變量。

2 全局敏感性分析方法

2.1 Morris法

Morris法[13]是定性全局敏感性分析方法，用于識別和篩選較敏感的輸入?yún)?shù)。假設(shè)模型有n個參數(shù)，參數(shù)的取樣點數(shù)為q，每個參數(shù)都從對應(yīng)的q個取樣點上隨機取值，獲得向量X=[x1,x2,…,xn]。構(gòu)造m×n(m=n+1)的矩陣，矩陣中相鄰2行只有1個參數(shù)不同。取構(gòu)造矩陣中相鄰兩行的組合參數(shù)作為模型參數(shù)值，輸入模型計算相應(yīng)的輸出結(jié)果y1和y2，獲得相應(yīng)參數(shù)的敏感度|Δy|/Δ。其中：Δy為輸出結(jié)果的變化量，Δ為輸入?yún)?shù)的變化量。由于Morris法的隨機性，很容易在一次隨機取樣及隨機化過程中出現(xiàn)誤差，所以需進行多次重復(fù)取樣，獲得單個參數(shù)敏感度的均值和方差。Morris[12]提出2個指標：1)均值，表征單個參數(shù)的敏感度，確定參數(shù)的敏感性排序；2)標準差，表征參數(shù)間相互作用的影響。該法的詳細原理和計算步驟見文獻[13,25]。

2.2 Sobol法

Sobol法的核心是把模型參數(shù)分解為單個參數(shù)及參數(shù)之間相互組合的函數(shù)[16]，當(dāng)參數(shù)正交時，模型具有唯一的分解形式。將模型總方差分解為單個參數(shù)單獨作用的方差和多個參數(shù)同時作用的方差：

(1)

其中：D為模型總方差；Di為參數(shù)xi單獨作用的方差；Di,j為參數(shù)xi和xj同時作用的方差；D1,2,…,n為n個參數(shù)同時作用的方差。

對上式歸一化：

(2)

可計算模型單個參數(shù)或參數(shù)之間相互作用的一階或高階敏感度，即Si=Di/D，Si,j=Di,j/D。 方程(1)可改寫為

(3)

其中：Si為1階敏感度；Si,j為2階敏感度；S1,2,…,n為n階敏感度。第i個參數(shù)的總敏感度定義為STi=∑S(i)，其中，S(i)為所有包含參數(shù)i的敏感度。

2.3 Extend FAST法

和

(4)

其中，Z0為非零整數(shù)。

該方法中參數(shù)的取樣數(shù)為

(5)

3 敏感性分析結(jié)果與討論

對第2節(jié)中已建立的數(shù)學(xué)模型采用Morris、Sobol及EFAST法探討各參數(shù)及參數(shù)間相互作用對模型輸出結(jié)果的影響。選取水平滲透率、孔隙度、殘余液體飽和度、孔隙分布指數(shù)、壓縮系數(shù)、進氣壓力的倒數(shù)和鹽度7個參數(shù)，各參數(shù)不確定性取值范圍見表1。

3.1 Morris分析結(jié)果

采用Morris法對7個參數(shù)重復(fù)抽樣10次獲得80組樣本，各響應(yīng)變量的敏感性分析結(jié)果如圖3所示。

圖3a為以注入井處的壓力為響應(yīng)變量的分析結(jié)果。由圖3a可見：孔隙度的敏感性最高，其與其他參數(shù)相互作用對模型輸出的影響也最為顯著；其次為壓縮系數(shù)和水平滲透率。無論是從均值還是標準差的角度來看，孔隙分布指數(shù)和鹽度的敏感性均較小?？紫抖燃皦嚎s系數(shù)影響CO2存儲空間，水平滲透率影響CO2的遷移速率，從宏觀上直接影響儲層中積累壓力的傳遞釋放，參數(shù)取值較大，可有效地減小儲層壓力積累。

圖3 以注入井處的壓力(a)、氣相CO2總量(b)及CO2氣相羽的擴散距離(c)為響應(yīng)變量的Morris敏感性分析結(jié)果Fig.3 Sensitivity analysis results of Morris based on the injection well pressure (a), the total mass of gas-phase CO2 (b) and the spread distance of CO2 gas plume (c) as response variables

圖3b為以氣相CO2總量為響應(yīng)變量的分析結(jié)果。由圖3b可見：水平滲透率的敏感性最大，其次為孔隙度和進氣壓力的倒數(shù)，這3個參數(shù)與其他參數(shù)相互作用對模型輸出的影響也最顯著；壓縮系數(shù)的敏感性最小。水平滲透率影響遷移速率，孔隙度影響儲存空間，進氣壓力表征毛管阻力，CO2橫向運移范圍越大，毛管阻力越小，孔隙可存儲的氣相CO2總量越多。

圖3c為以CO2氣相羽的擴散距離為響應(yīng)變量的分析結(jié)果。其中：水平滲透率的敏感性最高；其次為孔隙度和進氣壓力的倒數(shù)。滲透率與其他參數(shù)相互作用的影響最大，殘余液體飽和度與其他參數(shù)相互作用的影響次之。水平滲透率越大，CO2羽沿蓋層底部的橫向遷移速率及擴散距離越大；氣相CO2充滿空隙后會繼續(xù)向前運移，因此孔隙度越大，意味著驅(qū)替咸水后形成的CO2儲存空間越大，則CO2羽的擴散范圍越??；進氣壓力與毛管壓力有關(guān)，毛管壓力越大，CO2運移受到的阻力越大，遷移距離越小。

圖4 以注入井壓力(a)、氣相CO2總量(b)及CO2氣相羽擴散距離(c)為響應(yīng)變量的Sobol 1階敏感度及總敏感度分析結(jié)果Fig.4 First and total order sensitivity indices of Sobol based on the injection well pressure (a), the total mass of gas-phase CO2 (b) and the spread distance of CO2 gas plume (c) as response variables

總體而言，對于不同的響應(yīng)變量，水平滲透率與孔隙度對模型輸出的影響較大。考慮到Morris法不能定量說明參數(shù)間相互作用的不確定貢獻率，為此進一步采用Sobol和EFAST法進行參數(shù)的定量敏感性分析。

3.2 Sobol法分析結(jié)果

由于Sobol法中取樣數(shù)較少時，敏感度指數(shù)易受數(shù)值計算影響出現(xiàn)負值，所以本文經(jīng)過多次取樣計算測試穩(wěn)定性后，選取取樣個數(shù)為8 192。Sobol法的定量全局敏感性分析結(jié)果如圖4所示。對于不同的響應(yīng)變量，敏感參數(shù)的1階敏感性排序結(jié)果(圖4)與Morris分析結(jié)果(圖3)相近，參數(shù)間相互作用的影響程度有所不同。對比參數(shù)的1階敏感度及總敏感度可以看出，單個參數(shù)對模型輸出的影響貢獻是主要的，參數(shù)間相互作用對模型輸出也具有一定的影響，但計算結(jié)果顯示交互敏感度數(shù)值不大(以氣相CO2總量為響應(yīng)變量除外)，因此本文未進行2階敏感度分析。

Sobol的定量分析結(jié)果進一步表明，對于不同的響應(yīng)變量，參數(shù)的敏感性排序不同，說明參數(shù)敏感性也具有不確定性。敏感性分析結(jié)果揭示了各參數(shù)在CO2封存過程中的控制影響作用不同。如圖4a所示：以注入井處的壓力為響應(yīng)變量時，φ為最敏感參數(shù)，C為次敏感參數(shù)，其控制氣相CO2存儲空間，存儲空間大則不易產(chǎn)生壓力增量。圖4b中，以氣相CO2總量為響應(yīng)變量時，kx為最敏感參數(shù)，1/p0為次敏感參數(shù)。1/p0影響毛管壓力進而影響氣相CO2總量，但kx表現(xiàn)為正影響，而1/p0表現(xiàn)為負影響。需要指出的是，此時φ的1階敏感性較小，交互敏感性較大，與前述Morris的定性分析結(jié)果(圖3b)中φ排序第二差異顯著。這可能是因為Morris取樣較少導(dǎo)致個別參數(shù)的敏感性分析存在誤差，特別是當(dāng)該參數(shù)的交互敏感度較大時，數(shù)值計算的微小誤差使得Morris對孔隙度的分析結(jié)果產(chǎn)生了較大的偏差，而Sobol計算取樣次數(shù)較Morris多，結(jié)果更為可靠。圖4c中，以CO2氣相羽的擴散距離為響應(yīng)變量時，kx仍為最敏感參數(shù)，各參數(shù)間的交互敏感度較小，kx控制氣相CO2的遷移速率，遷移速率越大，CO2羽橫向擴散范圍越大。綜上所述，對于不同的響應(yīng)變量，各參數(shù)對決定CO2封存運移的主要因素的影響程度不同，在模型中的影響效果的表現(xiàn)也不同，如kx，對于注入井處的壓力表現(xiàn)為正影響，對于CO2羽擴散距離表現(xiàn)為負影響，易引發(fā)泄漏。對比分析Morris與Sobol的分析結(jié)果表明，當(dāng)模型參數(shù)較多、參數(shù)間的關(guān)聯(lián)性較小、相互作用較弱時，可先采用Morris篩選分析識別模型潛在的重要參數(shù)，以有效地減少參數(shù)維數(shù)降低計算時耗。

3.3 EFAST法分析結(jié)果

為了與Sobol法比較，依據(jù)EFAST法中參數(shù)的取樣算法(式(5))，設(shè)置取樣個數(shù)分別為2 023、3 031、4 039、5 047及7 007等，經(jīng)計算測試后發(fā)現(xiàn)EFAST法在取樣較少(3 031個)時計算結(jié)果即趨于穩(wěn)定。以注入井處的壓力和CO2氣相羽擴散距離為響應(yīng)變量的EFAST和Sobol的敏感度計算結(jié)果如表2所示(限于篇幅，其他結(jié)果未列出)，對于不同的響應(yīng)變量，EFAST計算所得出的各參數(shù)的主要敏感度及總敏感度分析結(jié)果與Sobol分析結(jié)果基本一致，由于多個參數(shù)(至少2個)的敏感度值在同一數(shù)量級上，數(shù)值誤差的影響較小，使得EFAST和Sobol的計算結(jié)果相近。Sobol和EFAST均定量給出了各參數(shù)的主要敏感度和總敏感度值，它們的差值即為參數(shù)的交互敏感度，即參數(shù)間相互耦合作用的間接影響，實質(zhì)上反映了模型“異參同效”的現(xiàn)象[18]。從表2中可見，當(dāng)以注入井處的壓力或以CO2氣相羽擴散距離為響應(yīng)變量時，參數(shù)的交互敏感度不大，說明參數(shù)的獨立性較強，參數(shù)間的相關(guān)性較弱，表現(xiàn)為“異參同效”性較小。

值得一提的是，從敏感度計算結(jié)果(表2)上來看，盡管Sobol和EFAST 2種方法的分析結(jié)果基本一致，但是EFAST較Sobol的取樣數(shù)要少得多，相應(yīng)的計算過程中前者較后者要節(jié)省大量的計算時耗。以Windows XP系統(tǒng)、Intel core2 E7400 CPU、主頻2.8 GHz、內(nèi)存2 GB為例，運行1次正問題需耗時2.0～3.5 min。本次研究中，對于以注入井壓力及CO2氣相羽擴散距離為響應(yīng)變量時，EFAST法較Sobol法計算分別節(jié)省時間約160.0和250.0 h。綜上所述，EFAST法需要的樣本數(shù)較少且計算高效、穩(wěn)健，對于多參數(shù)模型，可較大程度地降低計算時耗，因此建議當(dāng)模型參數(shù)較多時，可首選EFAST法進行定量全局敏感性分析。

表2 EFAST與Sobol敏感性分析結(jié)果比較

4 結(jié)論與展望

1)基于對蘇北盆地鹽城組下段砂巖儲層概化的數(shù)值模擬研究表明，CO2注入儲層后會產(chǎn)生壓力增量，氣相CO2在浮力作用下向上運移沿蓋層底部橫向擴散遷移。本文選擇儲層注入井處的壓力、氣相CO2總量及CO2氣相羽的擴散距離3個響應(yīng)變量進行定性和定量全局敏感性分析，篩選評估影響CO2封存和運移的決定因子。

2)Morris定性全局敏感性分析可有效篩選識別模型重要參數(shù)。Morris和Sobol全局敏感性分析結(jié)果都表明，對于不同的響應(yīng)變量，參數(shù)的敏感性排序不同，說明各參數(shù)對決定CO2封存過程主要因素的影響程度不同：以最大壓力為響應(yīng)變量時，孔隙度最敏感，其次為壓縮系數(shù)和水平滲透率；以氣相CO2總量為響應(yīng)變量時，水平滲透率最敏感；以CO2氣相羽的擴散距離為響應(yīng)變量時，水平滲透率最敏感，其次為孔隙度。參數(shù)的交互敏感度，實質(zhì)上體現(xiàn)了模型“異參同效”的現(xiàn)象，參數(shù)的交互敏感度小，表征參數(shù)的“異參同效”性較小。對于本例，以氣相CO2總量為響應(yīng)變量時，由于孔隙度的交互敏感度較高，取樣次數(shù)較少的Morris敏感性排序產(chǎn)生了較大的偏差，而Sobol計算取樣次數(shù)較多，結(jié)果更為可靠。因此對于Morris定性全局敏感性分析時需增加取樣次數(shù)以避免交互敏感度較高參數(shù)的影響。

3)與Morris法相比，Sobol和EFAST全局敏感性分析方法可進一步定量參數(shù)不確定性影響的貢獻率。EFAST需要的樣本數(shù)較少且計算高效、穩(wěn)健，建議當(dāng)模型參數(shù)較多時，可首選EFAST法進行定量全局敏感性分析。但需要說明的是，無論Sobol還是EFAST法都無法充分利用樣本的所有信息且計算量較大，對于多參數(shù)的模型評價計算費時。筆者擬下一步采用替代模型技術(shù)評價CO2封存過程中參數(shù)的敏感性以提高計算效率。

[1] Holloway S. Storage of Fossil Fuel-Derived Carbon Dioxide Beneath the Surface of the Earth[J/OL]. Annu Rev Energy Environ, 2001, 26:145-166,doi: 10.1146/annurev.energy.26.1.145.

[2] 李曉春, 劉延鋒, 白冰, 等. 中國深部咸水含水層CO2儲存優(yōu)先區(qū)選擇[J]. 巖石力學(xué)與工程學(xué)報, 2006, 25(5): 963-968. Li Xiaochun, Liu Yanfeng, Bai Bing, et al. Ranking and Screening of CO2Saline Aquifer Storage Zones in China[J]. Chinese Journal of Rock Mechanics and Engineering, 2006, 25(5): 963-968.

[3] 張旭輝, 魯曉兵, 劉慶杰. 蓋層特性對CO2埋存逃逸速度的影響[J]. 土工基礎(chǔ), 2009, 23(3): 67-70. Zhang Xuhui, Lu Xiaobing, Liu Qingjie. The Effect of the Characteristics of Cap on the Escaping Velocity of CO2[J]. Soil Eng and Foundation, 2009, 23(3): 67-70.

[4] 郭建強, 張森琦, 刁玉杰, 等. 深部咸水層CO2地質(zhì)儲存工程場地選址技術(shù)方法[J]. 吉林大學(xué)學(xué)報: 地球科學(xué)版, 2011, 41(4): 1084-1090. Guo Jianqiang, Zhang Senqi, Diao Yujie, et al. Site Selection Method of CO2Geological Storage in Deep Saline Aquifers[J]. Journal of Jilin University: Earth Science Edition, 2011, 41(4): 1084-1090.

[5] Pruess K.ECO2N:A TOUGH2 Fluid Property Mo-dule for Mixtures of Water, NaCl and CO2[R]. Berkeley: Lawrence Berkeley Laboratory，2005.

[6] Xu T F, Sonnenthal E, Nicolas S, et al. TOUGHREACT User’s Guide: A Simulation Program for Non-Isothermal Multiphase Reactive Geochemical Transport in Variably Saturated Geologic Media[R]. Berkeley: Lawrence Berkeley Laboratory，2008.

[7] Zhou Q L,Birkholzer J T,Mehnert E,et al.Modeling Basin and Plume-Scale Processes of CO2Storage for Full-Scale Deployment[J]. Groundwater, 2010, 48(4): 494-514.

[8] Birkholzer J T, Zhou Q, Cortis A, et al. A Sensitivity Study on Regional Pressure Buildup from Large-Scale CO2Storage Projects[C]//Gale J, Hendriks C, Turkenberg W, et al. Energy Procedia：Vol 4. Ne-therlands: Elsevier Ltd, 2011:4371-4378.

[9] Wiese B, Nimtz M, Klatt M, et al. Sensitivities of Injection Rates for Single Well CO2Injection into Saline Aquifers[J]. Chemie der Erde, 2010, 70(53): 165-172.

[10] 李義連, 房琦, 柯怡兵, 等. 高鹽度鹵水對CO2地質(zhì)封存的影響: 以江漢盆地潛江凹陷為例[J]. 地球科學(xué): 中國地質(zhì)大學(xué)學(xué)報, 2012, 37(2): 283-288. Li Yilian, Fang Qi, Ke Yibing, et al. Effect of High Salinity on CO2Geological Storage: A Case Study of Qianjiang Depression in Jianghan Basin[J]. Earth Science: Journal of China University of Geosciences, 2012, 37(2): 283-288.

[11] 柯怡兵, 李義連, 張煒, 等. 巖鹽沉淀對咸水層二氧化碳地質(zhì)封存注入過程的影響研究: 以江漢盆地為例[J]. 地質(zhì)科技情報, 2012, 31(3): 109-115. Ke Yibing, Li Yilian, Zhang Wei, et al. Impact of Halite Precipitation on CO2Injection into Saline Aquifers: A Case Study of Jianghan Basin[J]. Geological Science and Technology Information, 2012, 31(3): 109-115.

[12] 趙銳銳, 孟慶輝, 成建梅. 深部咸含水層CO2注入的流體遷移模擬研究: 以松遼盆地三肇凹陷為例[J]. 巖土力學(xué), 2012, 33(4): 1247-1252. Zhao Ruirui, Meng Qinghui, Cheng Jianmei. Fluid Migration Modeling of CO2Injection in Deep Saline Aquifers: A Case Study of the Sanzhao Depression, Songliao Basin[J]. Rock and Soil Mechanics, 2012, 33(4): 1247-1252.

[13] Morris M D. Factorial Sampling Plans for Preliminary Computational Experiments[J]. Technometrics, 1991, 33(2):161-174.

[14] Saltelli A, Anderes T H. Sensitivity Analysis of Model Output: An Investigation of New Techniques[J]. Comput Statist Data Anal, 1993, 15(2): 211-238.

[15] Saltelli A, Tarantola S, Chan K P S. A Quantitative Model-Independent Method for Global Sensitivity Analysis of Model Output[J]. Technometrics, 1999, 41(1): 39-56.

[16] Sobol I M.Sensitivity Estimates for Nonlinear Mathematical Model[J]. Mathematical Modeling and Computational Experiment, 1993, 1: 407-414.

[17] 宋曉猛, 孔凡哲, 占車生, 等. 基于統(tǒng)計理論方法的水文模型參數(shù)敏感性分析[J]. 水科學(xué)進展, 2012, 23(4): 503-510. Song Xiaomeng, Kong Fanzhe, Zhan Chesheng, et al. Sensitivity Analysis of Hydrological Model Parameters Using a Statistical Theory Approach[J]. Advances in Water Science, 2012, 23(4): 503-510.

[18] 任啟偉, 陳洋波, 舒曉娟. 基于Extend FAST方法的新安江模型參數(shù)全局敏感性分析[J]. 中山大學(xué)學(xué)報: 自然科學(xué)版, 2010, 49(3): 127-134. Ren Qiwei, Chen Yangbo, Shu Xiaojuan. Global Sensitivity Analysis of Xin’anjiang Model Parameters Based on Extend FAST Method[J]. Acta Scientiarum Naturalium Universitatis Sunyatseni, 2010, 49(3): 127-134.

[19] 任啟偉, 陳洋波, 周浩瀾, 等. 基于Sobol法的TOPMODEL模型參數(shù)全局敏感性分析[J]. 人民長江, 2010, 41(19): 91-107. Ren Qiwei, Chen Yangbo, Zhou Haolan, et al. Global Sensitivity Analysis of TOPMODEL Model Parameters Based on Sobol Method[J].Yangtze River, 2010, 41(19): 91-107.

[20] Crosetto M, Tarantola S. Uncertainty and Sensitivity Analysis: Tools for GIS-Based Model Implementation[J]. International Journal of Geographical Information Science, 2001, 15(5): 415-437.

[21] Meszaros R, Zsely I G, Szinyei D, et al. Sensitivity Analysis of an Ozone Deposition Model[J]. Atmospheric Environment, 2009, 43(3): 663-672.

[22] Pastres R, Chan K, Solidoro C, et al. Global Sensitivity Analysis of a Shallow Water 3-D Eutrophication Model[J]. Computer Physics Communications, 1999, 117(1): 62-74.

[23] 呂炳全, 王耀庭, 馮國清. 蘇北溱潼高郵拉分盆地的基本特征和油氣富集規(guī)律[J]. 西南石油學(xué)院學(xué)報, 1989, 11(2): 17-22. Lü Bingquan, Wang Yaoting, Feng Guoqing. Basic Characteristics of Qintong and Gaoyou Pull-Apart Basins in North Jiangsu and Their Laws of Oil and Gas Accumulation[J]. Journal of Southwestern Petroleum Institute, 1989, 11(2): 17-22.

[24] 朱靜昌, 張國棟, 王益友. 蘇北盆地阜寧群時期的古氣候和水介質(zhì)的物理化學(xué)條件分析[J]. 巖相古地理, 1992, 12(6): 8-16. Zhu Jingchang, Zhang Guodong, Wang Yiyou. Palaeoclimatic Conditions and Physicochemical Pro-perties of Water Bodies in the North Jiangsu Basin During the Funing Grounp Deposition[J]. Sedimentary Facies and Palaeogeography, 1992, 12(6): 8-16.

[25] 鄭菲, 施小清, 吳吉春, 等. 蘇北盆地鹽城組咸水層CO2地質(zhì)封存泄漏風(fēng)險的全局敏感性分析[J]. 高校地質(zhì)學(xué)報, 2012, 18(2): 232-238. Zheng Fei, Shi Xiaoqing, Wu Jichun, et al. Global Sensitivity Analysis of Leakage Risk for CO2Geolo-gical Sequestration in the Saline Aquifer of Yancheng Formation in Subei Basin[J]. Geological Journal of China Universities, 2012, 18(2): 232-238.

[26] 朱厚勤, 朱煜, 鄭開富. 蘇北盆地鹽城組天然氣藏成藏控制因素探討[J]. 海洋地質(zhì)動態(tài), 2003, 19(9): 22-26. Zhu Houqin, Zhu Yu, Zheng Kaifu. Pooling Conditions and Controlling Factors of Gas Pools in Yancheng Formation of Subei Basin[J]. Marine Geology Letters, 2003, 19(9): 22-26.

Global Parametric Sensitivity Analysis of Numerical Simulation for CO2Geological Sequestration in Saline Aquifers:A Case Study of Yancheng Formation in Subei Basin

Zheng Fei, Shi Xiaoqing, Wu Jichun, Zhao Liang, Chen Yang

SchoolofEarthSciencesandEngineering,NanjingUniversity,Nanjing210093,China

TOUGH2/ECO2N code was used to simulate the complex coupled processes taking place during and after CO2injection in the deep saline sandstone aquifer of Lower Yancheng Formation in Subei basin. Morris, Sobol and EFAST methods were used to analyze the global sensitivity of parameters, i.e.,kx,n,Slr,λ,C, 1/p0andS, on three response variables: 1)the pressure located in the injection well; 2)the amount of gas-phase CO2and 3)the spread distance of CO2gas plume. Results from the Morris and Sobol methods show thatnis the most sensitive parameter for the injection well pressure,kxis the most sensitive parameter both for the amount of gas-phase CO2and the spread distance of CO2gas plume. In summary, the sensitivity orders of these parameters are totally different for different corresponding variables. Comparative results from EFAST and Sobol show that the main and total effect obtained from the two quantitative sensitivity analysis methods are basically identical, however, EFAST is computationally more efficient than Sobol in numerical experiment.

CO2geological sequestration; global sensitivity analysis; Morris; Sobol; EFAST；Subei basin

10.13278/j.cnki.jjuese.201401209.

2013-05-04

國家自然科學(xué)基金項目(41172206)；江蘇省自然科學(xué)基金項目(BK2012313)；江蘇省產(chǎn)學(xué)研聯(lián)合創(chuàng)新資金計劃項目(BY2010136)

鄭菲(1987-)，女，博士研究生，主要從事多相流數(shù)值模擬方面的研究，E-mail:xiaomi2008happy@163.com

施小清(1979-)，男，副教授，博士，主要從事反應(yīng)溶質(zhì)運移模擬研究，E-mail:shixq@nju.edu.cn。

10.13278/j.cnki.jjuese.201401209

P641

A

鄭菲，施小清，吳吉春，等.深部咸水層CO2地質(zhì)封存數(shù)值模擬參數(shù)的全局敏感性分析:以蘇北盆地鹽城組為例.吉林大學(xué)學(xué)報:地球科學(xué)版,2014,44(1):310-318.

Zheng Fei,Shi Xiaoqing,Wu Jichun,et al.Global Parametric Sensitivity Analysis of Numerical Simulation for CO2Geological Sequestration in Saline Aquifers: A Case Study of Yancheng Formation in Subei Basin.Journal of Jilin University:Earth Science Edition,2014,44(1):310-318.doi:10.13278/j.cnki.jjuese.201401209.