李昱春， 張文兵， 任 杰， 李舸航， 甘 磊， 徐力群

(1.蘇州市吳中區(qū)水務(wù)局， 江蘇 蘇州 215128； 2.上海海事大學(xué) 海洋科學(xué)與工程學(xué)院， 上海 201306；3.河海大學(xué) 水利水電學(xué)院， 江蘇 南京 210098； 4.西安理工大學(xué) 水利水電學(xué)院， 陜西 西安 710048；5.上海市堤防泵閘建設(shè)運(yùn)行中心， 上海 200080)

1 研究背景

近年來(lái)，隨著河長(zhǎng)制、湖長(zhǎng)制以及長(zhǎng)江大保護(hù)等一系列重大戰(zhàn)略舉措的推行，我國(guó)對(duì)于河流的保護(hù)、治理和管理進(jìn)入了新的發(fā)展階段，對(duì)江河湖泊水域岸線的管理保護(hù)提出了更高的要求。河岸帶作為河流生態(tài)系統(tǒng)與陸地生態(tài)系統(tǒng)的過(guò)渡帶，對(duì)削減河流污染、調(diào)蓄洪水和保護(hù)水土環(huán)境具有重要作用，同時(shí)它也是動(dòng)植物生長(zhǎng)和棲息的重要場(chǎng)所[1-3]。因此，針對(duì)河岸帶相關(guān)問(wèn)題的研究受到越來(lái)越多學(xué)者的關(guān)注，并成為生態(tài)環(huán)境水力學(xué)研究中的熱點(diǎn)[4-6]。

河岸帶地下水與河流地表水時(shí)刻發(fā)生著能量的交換，溫度作為能量的直觀載體，是理想的天然示蹤劑[7]。目前，已有眾多學(xué)者將溫度示蹤法應(yīng)用于河流地表水與地下水交換速率及過(guò)程模式的研究中[8-13]。有關(guān)溫度示蹤法應(yīng)用于河岸帶相關(guān)問(wèn)題的研究已有大量報(bào)道，但多側(cè)重于試驗(yàn)手段。隨著河岸帶水熱耦合理論研究的深入，通過(guò)建立水熱耦合模型的數(shù)值方法來(lái)研究河岸帶內(nèi)部水熱動(dòng)態(tài)變化過(guò)程成為一種重要手段，例如，張文兵等[3]構(gòu)建了考慮土體非均質(zhì)傳熱的河岸帶水熱耦合模型，對(duì)比了不同土體有效導(dǎo)熱系數(shù)模型下河岸帶水熱耦合模型的模擬效果；Ren等[14]利用洞庭湖河岸帶某典型斷面的溫度和水位監(jiān)測(cè)資料，驗(yàn)證了所建立的水熱耦合模型的合理性。然而，河岸帶水熱耦合模型涉及諸多參數(shù)，參數(shù)反演分析比較困難，同時(shí)也會(huì)給模型的應(yīng)用增添較多的校正工作。因此，開(kāi)展河岸帶水熱耦合模型參數(shù)敏感性研究有利于快速確定相關(guān)參數(shù)的取值，為模型參數(shù)的校準(zhǔn)提供參考，對(duì)模型在實(shí)際工程中的廣泛應(yīng)用具有重要意義。

現(xiàn)有的敏感參數(shù)識(shí)別方法大致可分為單因素分析法和多因素分析法兩種。單因素分析法多適用于模型參數(shù)較少的情況，并且該方法假設(shè)模型參數(shù)之間無(wú)相互作用，而這一假設(shè)與實(shí)際情況不符[15]。相較之下，多因素分析法則彌補(bǔ)了單因素分析法的不足，能夠更加準(zhǔn)確、全面地反映敏感參數(shù)識(shí)別結(jié)果。在利用多因素分析法進(jìn)行敏感參數(shù)識(shí)別時(shí)，若對(duì)參數(shù)進(jìn)行全面組合，會(huì)帶來(lái)較大的工作量，正交試驗(yàn)設(shè)計(jì)則可以避免多因素試驗(yàn)的全面組合，并且能夠?qū)で笞顑?yōu)的參數(shù)組合來(lái)反映試驗(yàn)結(jié)果。為此，本文采用正交試驗(yàn)方法來(lái)識(shí)別河岸帶水熱耦合模型的敏感參數(shù)，以期為河岸帶水熱耦合模型參數(shù)選取、校正及應(yīng)用提供參考。

2 資料來(lái)源與研究方法

2.1 研究區(qū)域概況

沃克湖(Walker Lake)是美國(guó)內(nèi)華達(dá)州西部的一個(gè)重要水源，為當(dāng)?shù)氐霓r(nóng)業(yè)提供用水，并為當(dāng)?shù)囟鄠€(gè)社區(qū)使用的地下水含水層補(bǔ)水。沃克湖是大盆地中為數(shù)不多的多年生天然終端湖泊之一，終端湖泊是地形封閉盆地中地表水排水的終點(diǎn)。在自然條件下，湖面的蒸發(fā)量通常是湖水流出的全部組成部分。由于大盆地的高蒸發(fā)率，終端湖泊的水位和鹽度對(duì)河流流入的變化極為敏感。19世紀(jì)60年代開(kāi)始，為了支持日益增長(zhǎng)的農(nóng)業(yè)發(fā)展，當(dāng)?shù)剞r(nóng)耕者開(kāi)始從沃克湖引水。隨著時(shí)間的推移，從上游水庫(kù)和引水口流入沃克湖的水量減少，導(dǎo)致湖面下降了51.816 m，因此，沃克河流域的大部分水流只能來(lái)自內(nèi)華達(dá)山脈的融雪，但這也導(dǎo)致水體中總?cè)芙夤腆w(total dissolved solids, TDS)濃度大大增加，嚴(yán)重威脅水生生態(tài)系統(tǒng)的平衡，包括瀕危物種拉洪頓鱒魚(yú)的生存。由于消耗性用水，沃克湖和大盆地的其他終端湖泊的生態(tài)系統(tǒng)已處于危險(xiǎn)之中。為此，美國(guó)地質(zhì)調(diào)查局自2012年3月至2013年10月在沃克湖流域進(jìn)行了廣泛的實(shí)地研究，以量化史密斯谷、梅森谷和沃克湖谷灌溉區(qū)相關(guān)河流的滲漏損失。圖1為沃克湖流域圖。

圖1 沃克湖流域概況

2.2 數(shù)據(jù)來(lái)源

本文數(shù)據(jù)主要來(lái)源于美國(guó)沃克河流域2012年4月至2012年6月相關(guān)河流的水位及溫度變化數(shù)據(jù)。該部分?jǐn)?shù)據(jù)公開(kāi)于美國(guó)地質(zhì)調(diào)查局官網(wǎng)，允許研究人員下載獲取(網(wǎng)址：https://pubs.er.usgs.gov/publication/sir20165133)。

2.3 研究方法

2.3.1 河岸帶水熱耦合數(shù)學(xué)模型 河岸帶飽和-非飽和瞬態(tài)滲流場(chǎng)采用Richards方程描述[16]：

(1)

水力特性對(duì)飽和度和壓力的依賴性導(dǎo)致Richards方程的高度非線性，因此，學(xué)者們提出了各種數(shù)值和分析方法用于評(píng)估Richards方程，而最常用的是通過(guò)土壤水分特征曲線(soil water characteristic curve，SWCC)來(lái)研究水力特性。迄今，在用以描述土壤水分特征曲線的經(jīng)驗(yàn)公式中，具有代表性的有van Genuchten模型[17]、Gardner模型[18]和Brooks-Corey模型[19]。在眾多的模型中，van Genuchten模型具有精度高、參數(shù)物理意義明確、適用性較強(qiáng)[20-21]的優(yōu)點(diǎn)，其對(duì)負(fù)壓Hp<0(即非飽和條件)下水力特性的解析定義為：

θ=θr+Sr(n-θr)

(2)

(3)

(4)

kr(θ)=Sr1/2[1-(1-Sr1/m)m]2

(5)

式中：θr和θs分別為殘余與飽和體積含水率，m3/m3；n為孔隙率；hp為壓力水頭，m；α和β為與SWCC曲線有關(guān)的van Genuchten模型參數(shù)；m為與β相關(guān)的參數(shù)，m=1-1/β。

值得注意的是，當(dāng)Hp≥0時(shí)，假定θ=n，則有Sr=1，Cm=0和kr(θ)=1，方程式(1)將自動(dòng)轉(zhuǎn)化為完全飽和的達(dá)西方程。因而，對(duì)于河岸帶整體模型而言，該方程既能考慮河岸帶中飽和土體的水流運(yùn)動(dòng)，又能考慮河岸帶上覆非飽和土體的水流運(yùn)動(dòng)。

采用對(duì)流-傳熱方程描述河岸帶水體與土體間的熱量交換[22]：

(6)

式中：Cw、Cs分別為水與土體的體積熱容，J/(m3·℃)；λe為土體有效導(dǎo)熱系數(shù)，W/(m·℃)；DH為水動(dòng)力彌散系數(shù)，m2/s；u為平均流速，m/s；Qs為熱量源匯項(xiàng)，W/m3。

方程式(6)中的水動(dòng)力彌散系數(shù)DH可表示為：

(7)

式中：αT和αL分別為橫向和縱向彌散度，m； |v|為流速矢量的大小，m/s；δij為克羅內(nèi)克爾二階單位張量；vi、vj分別為速度矢量的第i個(gè)和第j個(gè)分量，m/s。

同樣地，當(dāng)孔隙處于飽和狀態(tài)，則θ=n，方程式(6)退化為飽和含水層對(duì)流-傳熱方程。

由于土體有效導(dǎo)熱系數(shù)在干燥狀態(tài)和飽和狀態(tài)下差異較大，僅設(shè)定為固定常數(shù)是不準(zhǔn)確的，為此引入土體有效導(dǎo)熱系數(shù)經(jīng)驗(yàn)?zāi)Ｐ陀糜诜从骋蚝什煌瑢?dǎo)致的非均質(zhì)參數(shù)空間分布問(wèn)題。目前，有關(guān)土體有效導(dǎo)熱系數(shù)預(yù)測(cè)的經(jīng)驗(yàn)?zāi)Ｐ捅姸?，其中Johansen土體有效導(dǎo)熱模型能夠較好地反映河岸帶土體非均質(zhì)傳熱[3]，其表達(dá)式為：

λe(θ)=(λsat-λdry)Ke+λdry

(8)

式中：λsat和λdry分別為飽和與干燥狀態(tài)的土體導(dǎo)熱系數(shù)，W/(m·℃)；Ke為歸一化導(dǎo)熱系數(shù)。這些系數(shù)的計(jì)算如公式(9)～(11)所示。

(9)

(10)

(11)

式中：λw為水的導(dǎo)熱系數(shù)，W/(m·℃)；λs為土體的導(dǎo)熱系數(shù)，W/(m·℃)；λs由土體中的石英含量(q)及其導(dǎo)熱系數(shù)(λq=7.7 W/(m·℃))和其他礦物的導(dǎo)熱系數(shù)(λ0)計(jì)算得到，即λs=λqqλ01-q(當(dāng)q>0.2時(shí)，λ0=2.0 W/(m·℃)；當(dāng)q≤0.2時(shí)，λ0=3.0 W/(m·℃))；ρb為土體的堆積密度，kg/m3。

2.3.2 敏感性分析 正交試驗(yàn)法是從大量的試驗(yàn)點(diǎn)中選取部分具有代表性的點(diǎn)進(jìn)行多因素試驗(yàn)，可以大大減少試驗(yàn)次數(shù)，減少工作量，目前已廣泛應(yīng)用于科學(xué)試驗(yàn)研究[20,23-26]。正交試驗(yàn)表是正交試驗(yàn)法的核心，以“七因素三水平”情況為例，其正交試驗(yàn)表可參考文獻(xiàn)[15]進(jìn)行設(shè)計(jì)。

極差分析法和方差分析法是分析正交試驗(yàn)結(jié)果的兩種常用方法。相較于方差分析法，極差分析法難以區(qū)分影響試驗(yàn)結(jié)果變化的因素，并且只能給出各因素的影響次序，而無(wú)法定量判斷該因素的敏感程度[24]。方差分析法則能較好地估計(jì)試驗(yàn)誤差的大小，并且能夠反映模型結(jié)果輸出對(duì)各參數(shù)敏感程度的量化結(jié)果。方差分析法的基本原理如下：

以Ln(rc)試驗(yàn)為例，將第k次試驗(yàn)的結(jié)果記為Yk(k=1, 2, ···,n)，令Tij為第j列因素第i個(gè)水平試驗(yàn)結(jié)果Yk之和，T為總試驗(yàn)結(jié)果之和，pij為因素j在i水平下的試驗(yàn)次數(shù)，則：

(12)

(13)

pij=n/r

(14)

(15)

將試驗(yàn)中n個(gè)試驗(yàn)結(jié)果的總變差記為ST，表示全部試驗(yàn)結(jié)果之間的差異程度；第j列變差平方和記為Sj，表示第j列因素不同水平之間的差異程度；總和記為Se，表示試驗(yàn)條件的差異程度。則：

(16)

(17)

(18)

假定ST、Sj和Se的自由度分別為fT、fj和fe，則：

fT=n-1

(19)

fj=r-1

(20)

(21)

計(jì)算過(guò)程中，各試驗(yàn)結(jié)果(Y1，Y2，…,Yn)相互不影響，并且服從同方差σ2的正態(tài)分布，則可構(gòu)造F檢驗(yàn)的統(tǒng)計(jì)量：

(22)

通過(guò)將Fj值與F分布表中查找的檢驗(yàn)臨界值Fα(fj,fe)進(jìn)行比較，則可以判斷模型輸出結(jié)果對(duì)各因素變化的敏感程度。

2.3.3 計(jì)算模型 以文獻(xiàn)[3]中的工程實(shí)例即沃克湖福克斯1號(hào)灌渠作為本文的計(jì)算模型。

(1)幾何模型。?？怂构嗲敲绹?guó)內(nèi)達(dá)華州沃克湖的支流之一，美國(guó)內(nèi)華達(dá)水科學(xué)中心的研究人員在此布設(shè)有溫度和水位監(jiān)測(cè)裝置，用于監(jiān)測(cè)該河流地表水與地下水的溫度和水位[26]。本文選取?？怂?號(hào)灌渠的某河岸帶典型剖面進(jìn)行水熱耦合模型敏感參數(shù)識(shí)別，模型計(jì)算范圍如圖2所示。圖2中BP1和BP2分別為距離河道中心2.8和5.8 m的溫度測(cè)桿，每個(gè)測(cè)桿內(nèi)懸掛有3個(gè)溫度傳感器，用于測(cè)量河岸帶沉積層的溫度變化。另外，河道中還布設(shè)有溫度和水位記錄儀用于記錄河水溫度和水位。環(huán)境溫度則是由埋設(shè)在河岸帶表層土體中的溫度傳感器進(jìn)行記錄。圖3給出了?？怂?號(hào)灌渠在2012年4-5月的河水溫度、水位及環(huán)境溫度實(shí)測(cè)變化曲線。

(2)計(jì)算參數(shù)。本文河岸帶水熱耦合模型計(jì)算參數(shù)取自于文獻(xiàn)[3]、[27]～[29]中給出的相關(guān)參數(shù)，具體數(shù)值見(jiàn)表1。

圖2 模型計(jì)算范圍示意圖[3](A～F為邊界編號(hào))

圖3 2012年4-5月福克斯1號(hào)灌渠河水溫度、水位及環(huán)境溫度實(shí)測(cè)變化曲線

表1 模型計(jì)算參數(shù)

(3)邊界條件及初始條件。對(duì)于河岸帶滲流場(chǎng)，模型左(AF)、右(BC)及表層土體邊界(CD和EF)設(shè)為無(wú)流動(dòng)邊界；土-水接觸邊界(DE)設(shè)為水頭邊界；底邊界(AB)設(shè)為透水層邊界。對(duì)于河岸帶溫度場(chǎng)，模型左(AF)、右(BC)及底邊界(AB)設(shè)為絕熱邊界；表層土體邊界(CD和EF)和土-水接觸邊界(DE)均設(shè)為溫度邊界，將變化的環(huán)境溫度和河水溫度分別施加于對(duì)應(yīng)的邊界上。滲流場(chǎng)的初始條件是將前一時(shí)刻實(shí)測(cè)的壓力水頭通過(guò)線性插值施加于計(jì)算域中；溫度場(chǎng)的初始條件為各區(qū)域模擬前一時(shí)刻的實(shí)測(cè)平均溫度。

3 結(jié)果與分析

3.1 模型驗(yàn)證

為了驗(yàn)證河岸帶水熱耦合模型的準(zhǔn)確性及定量評(píng)價(jià)模型的模擬效果，引入均方根誤差(RMSE)、決定系數(shù)(R2)和相對(duì)誤差(Re)作為模型評(píng)價(jià)指標(biāo)[29]。圖4給出了河岸帶各監(jiān)測(cè)點(diǎn)溫度模擬值與實(shí)測(cè)值的對(duì)比。

圖4 2012年4-5月河岸帶各監(jiān)測(cè)點(diǎn)溫度模擬值與實(shí)測(cè)值對(duì)比

由圖4可以看出，基于水熱耦合模型模擬得到的河岸帶各測(cè)點(diǎn)的溫度模擬值與現(xiàn)場(chǎng)實(shí)測(cè)值結(jié)果較為吻合，并且具有較好的一致性，表明該模型能夠較為準(zhǔn)確地反映河岸帶水熱動(dòng)態(tài)變化過(guò)程。為了定量描述模型的準(zhǔn)確性，表2給出了河岸帶各監(jiān)測(cè)點(diǎn)溫度模擬值與實(shí)測(cè)值對(duì)比的評(píng)價(jià)指標(biāo)統(tǒng)計(jì)結(jié)果。

由表2可以看出，河岸帶水熱耦合模型模擬結(jié)果的RMSE值介于0.16～1.67 ℃，R2值介于0.65～0.99，Re值介于1.87%～12.65%，各項(xiàng)評(píng)價(jià)指標(biāo)均表現(xiàn)出較好的結(jié)果，僅在測(cè)點(diǎn)BP2-1.50 m處的結(jié)果表現(xiàn)稍差，但亦在指標(biāo)控制的合理范圍內(nèi)。因此，河岸帶水熱耦合模型的合理性得到了驗(yàn)證。

3.2 河岸帶水熱耦合模型敏感參數(shù)識(shí)別

基于上述河岸帶水熱耦合模型，將敏感性分析對(duì)象定為影響河岸帶溫度(T)和側(cè)向潛流交換速率(|vx|)所涉及的7個(gè)模型參數(shù)，具體包括水力傳導(dǎo)率(ks)、飽和含水率(θs)、殘余含水率(θr)、孔隙率(n)、van Genuchten模型參數(shù)(α、β)、土體的體積熱容(Cs)。計(jì)算過(guò)程中，每個(gè)參數(shù)選取3個(gè)水平，各水平分別取自基準(zhǔn)值以及基準(zhǔn)值分別增減10%以后的數(shù)值?？紤]到河岸帶的范圍[30]，將敏感性研究區(qū)域限定在圖2，x∈[5.13, 6.14]且y∈[-2.90,1.10]的模型尺寸范圍中，將因參數(shù)變動(dòng)的溫度(側(cè)向潛流交換速率)與額定溫度(額定潛流交換速率)在對(duì)象區(qū)域每個(gè)有限元單元節(jié)點(diǎn)上做差值，然后取平均值來(lái)表示溫度(側(cè)向潛流交換速率)的變化。

本文選擇L18(2×37)正交表設(shè)計(jì)正交試驗(yàn)，由于只有7個(gè)試驗(yàn)因素，故將第1列設(shè)為空列，并將各因素的不同水平對(duì)應(yīng)填入，得到識(shí)別模型敏感參數(shù)的正交表，如表3所示，表3中的每1行即為1組試驗(yàn)工況。

表2 模型評(píng)價(jià)指標(biāo)統(tǒng)計(jì)

表3 敏感參數(shù)識(shí)別正交試驗(yàn)方案及試驗(yàn)結(jié)果

按照表3設(shè)計(jì)的試驗(yàn)方案分別計(jì)算每種方案下的試驗(yàn)指標(biāo)T和|vx|，并將它們列在表3的最后2列。根據(jù)得到的各組工況下的試驗(yàn)結(jié)果，通過(guò)方差分析可得到各參數(shù)對(duì)試驗(yàn)指標(biāo)的敏感性。

根據(jù)方差分析原理，本文選用α=0.10和α=0.05作為顯著性檢驗(yàn)水平，并將計(jì)算所得的統(tǒng)計(jì)量Fj與選取的兩種顯著性檢驗(yàn)水平的F值進(jìn)行比較，從而確定各參數(shù)的顯著性。這里將顯著性水平劃分為3個(gè)等級(jí)：(1) 當(dāng)Fj>F0.05(2,3)時(shí)，表明因素敏感性非常高，影響高度顯著；(2)當(dāng)F0.05(2,3)>Fj>F0.1(2,3)時(shí)，表明因素敏感性中等，影響一般顯著；(3)Fj

表4 各因素對(duì)指標(biāo)T和|vx|的敏感性分析結(jié)果

由表4可以看出，對(duì)河岸帶溫度而言，參數(shù)敏感性由高到低依次為n、θs、θr、Cs、ks、β、α，其中孔隙率n對(duì)河岸帶溫度的影響高度顯著，即敏感程度高；參數(shù)θs、θr、Cs、ks、β、α對(duì)河岸帶溫度的影響不顯著，即敏感程度低。對(duì)于側(cè)向潛流交換速率而言，參數(shù)敏感性由高到低依次為ks、θs、β、Cs、α、n、θr，其中水力傳導(dǎo)率ks對(duì)河岸帶側(cè)向潛流交換速率的影響高度顯著，即敏感程度高；參數(shù)θs、β、Cs、α、n、θr對(duì)側(cè)向潛流交換速率的影響不顯著，即敏感程度低。

4 討 論

雖然已有部分學(xué)者評(píng)估了河岸帶水熱耦合模型參數(shù)對(duì)計(jì)算輸出結(jié)果的敏感性[9,26,31]，但不同的建模方法和計(jì)算方法均會(huì)對(duì)敏感性分析結(jié)果產(chǎn)生影響。Noranjo等[27]使用VS2DH軟件對(duì)構(gòu)建的水熱耦合模型參數(shù)敏感性進(jìn)行了分析，結(jié)果表明水力傳導(dǎo)率ks和VG模型參數(shù)α是影響模型輸出的敏感參數(shù)，其中ks是最重要和最敏感的參數(shù)。然而，該結(jié)果是基于單因素敏感性分析獲得，未考慮多參數(shù)間的相互作用。Ren等[9]基于HYDRUS-2D構(gòu)建了室內(nèi)水槽試驗(yàn)的河岸帶水熱耦合模型，分析了河岸帶水熱動(dòng)態(tài)變化過(guò)程，并對(duì)模型參數(shù)進(jìn)行了單因素和多因素敏感性分析。該研究結(jié)果不僅揭示了對(duì)模型溫度輸出敏感的內(nèi)部因素(水力傳導(dǎo)率ks和VG模型參數(shù)α)，而且研究了水頭和溫度等外部因素對(duì)模型輸出結(jié)果的影響。在最新的研究中，Ren等[31]采用Morris法對(duì)不同土體導(dǎo)熱系數(shù)模型下的水熱耦合模型參數(shù)進(jìn)行了敏感性分析，認(rèn)為水力傳導(dǎo)率ks和孔隙率n對(duì)模型輸出結(jié)果較為敏感。這與他們之前的敏感參數(shù)識(shí)別結(jié)果存在差異，而這種差異可能是由數(shù)學(xué)模型和導(dǎo)熱系數(shù)模型之間的差異引起的。因此，不同的建模方法或計(jì)算方法可能會(huì)對(duì)敏感性分析結(jié)果產(chǎn)生影響。但是，結(jié)合本文的結(jié)果，一般來(lái)說(shuō)，水力傳導(dǎo)率ks和孔隙率n仍然是最重要和最敏感的參數(shù)，這在一定程度上反映了本文結(jié)果的準(zhǔn)確性。

5 結(jié) 論

(1)對(duì)于河岸帶水熱耦合模型而言，孔隙率n是影響模型溫度輸出的主要參數(shù)，水力傳導(dǎo)率ks是影響模型側(cè)向潛流交換速率輸出的主要參數(shù)。其他參數(shù)相對(duì)這兩個(gè)指標(biāo)在模型輸出結(jié)果中的影響較小。

(2)在對(duì)河岸帶水熱耦合模型參數(shù)進(jìn)行反演分析時(shí)，可著重對(duì)敏感程度高的參數(shù)進(jìn)行反演，而對(duì)于敏感程度較低的參數(shù)，則可通過(guò)工程類(lèi)比的方法進(jìn)行推求，進(jìn)而提高模型參數(shù)的反演效率，并且減少模型校正的工作量。

(3)河岸帶水熱耦合模型處于不斷發(fā)展中，未來(lái)的模型將更符合實(shí)際工程，而模型敏感參數(shù)的識(shí)別受模型構(gòu)建方法等因素的影響，需對(duì)不同構(gòu)建方法下的模型參數(shù)進(jìn)行敏感識(shí)別，以期為模型在實(shí)際工程中的應(yīng)用提供指導(dǎo)。另外，現(xiàn)有的敏感性分析方法種類(lèi)繁多，并且模型所涉及的參數(shù)越來(lái)越多，需考慮發(fā)展適用于多參數(shù)耦合作用下的敏感性參數(shù)分析方法。