湖泊水動力模型外部輸入條件不確定性和敏感性分析

李一平,邱 利,唐春燕,布旻晟,田 威,余鐘波,Kumud Acharya(1.河海大學(xué)淺水湖泊綜合治理與資源開發(fā)教育部重點實驗室,江蘇 南京 210098；2.河海大學(xué)環(huán)境學(xué)院,江蘇 南京 210098；.河海大學(xué)水文水資源與水利工程科學(xué)國家重點實驗室,江蘇 南京 210098；.河海大學(xué)機電工程學(xué)院,江蘇 南京 210098；.江蘇省水文水資源勘測局,江蘇 南京 210029)

湖泊水動力模型外部輸入條件不確定性和敏感性分析

李一平1,2,3*,邱 利1,2,唐春燕1,2,布旻晟4,田 威5,余鐘波3,Kumud Acharya3(1.河海大學(xué)淺水湖泊綜合治理與資源開發(fā)教育部重點實驗室,江蘇 南京 210098；2.河海大學(xué)環(huán)境學(xué)院,江蘇 南京 210098；3.河海大學(xué)水文水資源與水利工程科學(xué)國家重點實驗室,江蘇 南京 210098；4.河海大學(xué)機電工程學(xué)院,江蘇 南京 210098；5.江蘇省水文水資源勘測局,江蘇 南京 210029)

以我國典型的大型淺水湖泊太湖為研究區(qū)域,采用國內(nèi)外常用的環(huán)境流體動力學(xué)模型(EFDC),結(jié)合拉丁超立方取樣(LHS)方法,研究湖泊水動力模型中4個重要的外部輸入條件,即3個邊界輸入條件(出入湖流量、風(fēng)速、風(fēng)向)和1個初始輸入條件(初始水位),對模型水動力模擬結(jié)果(水位、水齡以及流場)的影響與貢獻.結(jié)果表明,初始水位的設(shè)定對模擬全湖水位和水齡產(chǎn)生決定性影響,不確定性的貢獻率分別達到85.73%和66.125%,對垂向平均流速影響的貢獻率只有3%;風(fēng)速對表面流速模擬結(jié)果影響較大,貢獻率達到58.70%,而對水位和水齡的貢獻率分別為 5.25%和 3.00%.在垂向上,各層流速受外部輸入條件不確定性的影響規(guī)律相似,貢獻率排序為風(fēng)速(55%～60%)>風(fēng)向(10%～15%)>初始水位≈出入湖流量(1%～5%).因此在模擬大型淺水湖泊水動力過程時,可以根據(jù)不同的輸出目標(biāo)能夠有針對性地提高外部輸入條件的準(zhǔn)確度,為提高模型精確度提供有效信息.

輸入條件；不確定性分析；敏感性分析；environmental fluid dynamic code (EFDC)；大型淺水湖泊

隨著計算、監(jiān)測和通訊技術(shù)的進步,河湖水環(huán)境模擬預(yù)測技術(shù)也不斷提高,目前模擬河湖水動力過程的模型有很多,比如QUAL2K(River and Stream Water Quality Model)、MIKE11、WASP (Water Quality Analysis Simulation Program)和EFDC(Environmental Fluid Dynamic Code)[1].為了使水動力模型預(yù)測的結(jié)果能夠更好地服務(wù)于水環(huán)境預(yù)測及管理,提高模型輸出結(jié)果的準(zhǔn)確度顯得十分重要.然而很多因素會影響模型結(jié)果的準(zhǔn)確度,總的來說不確定性來源于三大方面:參數(shù)估計、外部輸入條件以及模型結(jié)構(gòu)[2].外部輸入條件包括邊界輸入條件和初始輸入條件,然而監(jiān)測技術(shù)、監(jiān)測設(shè)備以及操作人員水平等的局限性導(dǎo)致外部輸入條件值本身存在不確定性[3],而且在不同文獻中的取值也存在差別. Ahsan等[3]研究表明Mississippi海峽的水動力模型對水深和淡水入流量十分敏感,Mississippi海峽的淺水區(qū)域?qū)λ顥l件十分敏感,水深對該淺水區(qū)域溫度和鹽度輸出結(jié)果的不確定性貢獻率達到 88%.大型淺水湖泊(以太湖為例)因為“面大、水淺”的特性和周邊河網(wǎng)的密集性,其水動力過程復(fù)雜多變,增加了水動力過程的模擬難度,因此研究外部輸入條件對大型淺水湖泊水動力模擬結(jié)果的不確定性和敏感性影響具有重要的意義[4-6].目前用于不確定性分析方法很多,如敏感性分析、一階誤差分析、蒙特卡羅方法、Bootstrap 方法、最大似然方法、貝葉斯法、離散Bayes 法、區(qū)域敏感性分析方法、神經(jīng)網(wǎng)絡(luò)法、遺傳算法、模糊數(shù)學(xué)法、傅立葉敏感性檢驗法等. 國內(nèi)外學(xué)者在這方面做了很多研究工作,其中敏感性分析、FOEA、蒙特卡羅方法和Bootstrap方法是不確定性分析最常用的方法[7].而拉丁超立方抽樣法(Latin hypercube sampling, LHS)是一種分層的蒙特卡羅(Monte Carlo)抽樣技術(shù),與 Monte Carlo隨機抽樣法相比,LHS的樣本能夠更加精確地反映外部輸入條件的概率函數(shù)分布,能夠利用最少的樣本數(shù)反映外部輸入條件最大的信息量[8]. LHS最早由McKay提出,由Stein給出比較數(shù)學(xué)化的表述.LHS是一種可以替代Monte Carlo方法的效果好的方差縮減技術(shù),在仿真模擬、優(yōu)化計算和可靠性計算方面得到較為廣泛的應(yīng)用[9].本研究以我國典型的大型淺水湖泊太湖為例,選取湖泊水動力模型中4個重要的外部輸入條件,即3個邊界輸入條件(出入湖流量、風(fēng)速、風(fēng)向)和1個初始輸入條件(初始水位),利用LHS方法對外部輸入條件取值隨機取樣并進行不確定分析,并且利用標(biāo)準(zhǔn)秩逐步回歸法對外部輸入條件進行敏感性分析,量化每個外部輸入條件對模擬結(jié)果不確定性的貢獻率.以期通過定量分析每個外部輸入數(shù)據(jù)對模型結(jié)果不確定性的貢獻率,找出對模型結(jié)果不確定性影響較大的外部輸入條件,從而為有效地提高模型模擬精度提供可靠的信息.

1 研究區(qū)域

圖1 研究區(qū)域Fig.1 Study area

太湖(119°08′E～122°55′E,30°05′N～32°08′)是我國第三大淺水湖泊,湖泊總面積 2338km2.太湖平均水深 1.9m,最大水深為 2.6m,太湖底部呈現(xiàn)坡度為0°0′19.66′的平底地形,水深小于1.5m的淺水區(qū)域主要位于東部湖區(qū),占總表面積的19.3%.水深最大(>2.5m)的區(qū)域主要位于西北湖區(qū),占總表面積的 8.4%.太湖周邊水系復(fù)雜,本研究將環(huán)湖河道概化為30條河,編號如圖1所示.太湖流域夏季主導(dǎo)風(fēng)向是東南風(fēng),冬季主導(dǎo)風(fēng)向是西北風(fēng),平均風(fēng)速 3.5～5.0m/s[10-11].根據(jù)太湖各個湖區(qū)的特點,本研究將太湖劃分為 8個子區(qū)域(圖1),竺山湖、梅梁灣、貢湖、西北湖區(qū)、西南湖區(qū)、湖心區(qū)、東部湖區(qū)和東太湖灣[12].

2 研究方法

2.1 模型建立

EFDC(Environmental Fluid Dynamic Code)是國內(nèi)外常用的環(huán)境流體動力學(xué)模型,本研究采用EFDC模型來模擬太湖的水動力過程.模型采用笛卡爾直角坐標(biāo)網(wǎng)格,共劃分 4464個網(wǎng)格,每個網(wǎng)格單元邊長為 750m.為了較好地模擬湖底地形,垂直方向采用σ坐標(biāo),平均分為3層.根據(jù)流體靜力學(xué)連續(xù)性和避免產(chǎn)生 σ坐標(biāo)帶來的壓力梯度錯誤,應(yīng)使湖底坡度小于0.33[12].模型計算時間步長為10s,模擬時間為365d.

2.2 不確定性結(jié)果統(tǒng)計分析方法

在 LHS方法計算結(jié)果統(tǒng)計中,外部輸入條件個數(shù)記為 k,抽樣組數(shù)記為 n,每一個外部輸入條件X1、X2,…,Xk在各自的取值范圍內(nèi)等分成N組,每組隨機抽取一個值,n個k維變量組值共產(chǎn)生n個預(yù)測值,然后將n個預(yù)測值按大小排列,并分配給最小的預(yù)測值的累積概率為 1/n,分配給次最小的預(yù)測值的累積概率為 2/n,依此類推得預(yù)測值的經(jīng)驗分布函數(shù),此經(jīng)驗分布函數(shù)提供了子樣的分位值,即第m個預(yù)測值是m/n×100%的子樣分位值[13-14].其中5%、95%百分位值代表了由邊界條件引起的不確定性的邊界,5%代表下邊界, 95%代表上邊界[15].

2.3 敏感性分析方法

敏感性分析采用標(biāo)準(zhǔn)秩逐步回歸法,該方法將原有數(shù)據(jù)轉(zhuǎn)化為相關(guān)秩系數(shù),再進行標(biāo)準(zhǔn)化回歸分析,解決了輸入和輸出之間的非線性關(guān)系,得到標(biāo)準(zhǔn)回歸系數(shù)(SRRC)和決定系數(shù)(R2).SRRC2用來表示每個外部輸入條件對整個輸出結(jié)果方差的貢獻率,其值越大則表明該輸入條件越敏感,對模型結(jié)果不確定性的貢獻率也越大.R2是表示回歸模型分析的可行性程度,其值越大回歸分析越可行,一般R2值大于0.7,即可認(rèn)為回歸分析合理[16].

2.4 外部輸入條件的統(tǒng)計特征及LHS抽樣方法

本研究選取出3個邊界條件(出入湖流量、風(fēng)速、風(fēng)向)和1個初始條件(初始水位)作為外部輸入條件.表 1給出了各外部輸入條件取值的最大值和最小值,假設(shè)各外部輸入條件之間是相互獨立的,且在取值范圍內(nèi)服從均勻分布.環(huán)湖河道概化為30條河,其值以2005年實測數(shù)據(jù)為基礎(chǔ),模擬抽樣在流量 90%～110%范圍內(nèi)進行.太湖夏季主導(dǎo)風(fēng)向為東南風(fēng),本研究采用東風(fēng)到南風(fēng)(包括東南風(fēng))之間的風(fēng)向作為風(fēng)向取值范圍.

表1 LHS抽樣分析模型外部輸入條件的范圍Table 1 Ranges of boundary and initial conditions used for the Latin Hypercube Sampling

本研究采用 LHS抽樣方法隨機抽樣生成200組外部輸入條件[14-15],每組外部輸入條件數(shù)據(jù)都用EFDC模型運行一次,最后得到200組特定目標(biāo)的輸出結(jié)果,利用這200組結(jié)果定量分析4個外部輸入條件對模型模擬結(jié)果的不確定性影響.考慮到太湖是典型的風(fēng)生流,湖流垂向可能存在切變[17].因此本研究選取水位、水齡和垂向3層流速(表層、中層和底層)作為模型特定的模擬目標(biāo).最后利用 surfer 軟件分析外部輸入條件對整個湖區(qū)不同模擬目標(biāo)的不確定性影響[17].

3 結(jié)果與討論

3.1 外部輸入條件不確定性分析

本研究利用EFDC運行第365d的水位、水齡和垂向 3層流速作為輸出目標(biāo),通過分析 200組模型輸出結(jié)果,定量分析各外部輸入條件對模擬目標(biāo)的影響,并進行不確定性分析.外部輸入條件不確定性以模擬目標(biāo)值的均值、5%、95%分位值以及方差來進行量化描述.

圖2 模擬水位的均值, 5%分位值, 95%分位值和方差的空間分布Fig.2 Mean, 5%, 95% and variance of the simulated water level

3.1.1 以水位為輸出目標(biāo) 200組 4個外部輸入條件的 LHS抽樣組合的水位模擬結(jié)果的均值、5%和 95%百分位值如圖 2所示,結(jié)果表明±10%的出入湖流量變化以及初始水位、風(fēng)速和風(fēng)向在一定范圍內(nèi)的變化會造成模擬水位的分異性.均值、5%和95%的水位值都呈現(xiàn)出明顯的空間分異性,幾個典型湖區(qū)的水位分別如下:西北湖區(qū)分別為 3.4,2.8,3.9m,湖心區(qū)分別為 3.3,2.7, 3.8m,東太湖分別為 3.1,2.2,3.7m.不同頻率下的水位大致呈現(xiàn)出由東南至西北方向逐漸雍高的趨勢,其水位變化梯度也不相同, 95%水位空間分布變化梯度最大,水位均值空間變化梯度次之,5%水位空間變化梯度最小.圖2(b)給出了200組模擬水位方差的空間分布,量化了 4個外部輸入條件對水位結(jié)果的不確定性影響.方差值越大表明外部輸入條件對水位計算結(jié)果產(chǎn)生的不確定性也越大.湖心區(qū)和貢湖灣的水位方差相對較小約為 0.08,沿著東南和西北方向水位方差分別逐漸增加,特別是東南方向,方差變化范圍為[0.08,0.18].整個太湖水位都受到4個外部輸入條件不確定性的影響,特別是東太湖水位受不確定性的影響最大,原因可能是在風(fēng)的作用下,太湖東南部水位最低,可推測湖泊水位較低的區(qū)域受外部輸入條件的不確定性影響較大.綜上所述,當(dāng)以水位為輸出目標(biāo)時:整個湖區(qū)的水位都受到了 4個外部輸入條件的不確定性影響;東太湖受不確定性影響明顯大于湖心區(qū)和其他湖區(qū).

3.1.2 以水齡為輸出目標(biāo) 200組4個外部輸入條件的 LHS抽樣組合對輸出水齡值產(chǎn)生的不確定性結(jié)果如圖 3所示,結(jié)果同樣表明外部輸入條件對水齡產(chǎn)生了不確定性影響.模擬水齡的均值、5%分位和 95%分位的水齡值都呈現(xiàn)出明顯的空間分異性,幾個典型湖區(qū)的水齡分別如下:竺山灣分別為70,50,100d,西南湖區(qū)分別為180,150,200d,東太湖分別為190,200,240d.竺山灣和西北湖區(qū)的模擬水齡值最小,沿著西南方向水齡逐漸增加,在東太湖區(qū)域模擬水齡值最大.對于同一湖區(qū),不同頻率下的水齡表現(xiàn)出的規(guī)律是:95%分位的水齡值>平均水齡值>5%分位的水齡值.200組抽樣組合的水齡模擬結(jié)果方差空間分布[圖 3(b)]量化了不確定性影響.竺山灣和西北湖區(qū)的水齡方差最大約為800,而其他湖區(qū)的水齡方差相對較小約為300,整個太湖的水齡模擬結(jié)果方差變化范圍為[100,900],遠大于水位模擬結(jié)果的方差變化范圍.這表明水齡受外部輸入條件的不確定性影響很大,尤其是竺山灣和西北湖區(qū)受不確定性的影響最大.原因可能是太湖入湖河流主要集中在竺山灣和西北湖區(qū)附近,該區(qū)域水動力擾動最強,水體交換速率較快,因此可以推測在湖泊入湖河道附近水動力擾動較強的區(qū)域,水齡受外部輸入條件不確定性的影響較大.綜上所述,當(dāng)以水齡為輸出目標(biāo)時:外部輸入條件對水齡產(chǎn)生的不確定性影響很顯著(遠遠大于水位);竺山灣和西北湖區(qū)受外部輸入條件的不確定性影響明顯大于其他湖區(qū). 3.1.3 以流速為輸出目標(biāo) 200組4個外部輸入條件的LHS抽樣組合對表層流速產(chǎn)生的不確定性結(jié)果如圖4所示.模擬表層流速的均值、5%分位和95%分位的水齡值都呈現(xiàn)出明顯的空間分異性,在近岸湖區(qū)分別為 3,1,5m/s,在遠離湖岸區(qū)域分別為2,0.5,3m/s,在湖心區(qū)分別為1,0.3, 2m/s.表層流速在沿岸湖區(qū)相對較大,由湖岸到湖心方向流速逐漸減小,并且不同分位下表層流速的空間變化梯度也不相同,5%分位的流速值<均值<95%分位的流速值,這表明外部輸入條件的變化對表層流速模擬結(jié)果產(chǎn)生較大的不確定性.200組抽樣組合的表層流速模擬結(jié)果方差的空間分布如圖 4(b)所示,沿岸區(qū)域的表層流速方差最大約為10,而其他區(qū)域的方差相對較小,不確定性方差范圍[0,15].表明沿岸區(qū)域表層流速受到外部輸入條件的不確定性影響較大,而湖心區(qū)受到的不確定性的影響較小.

圖3 模擬水齡的均值, 5%分位值, 95%分位值和方差的空間分布Fig.3 Mean, 5%, 95% and variance of the simulated water age

表層流速、水位和水齡的不確定性空間分布形式各不相同,外部輸入條件對水齡產(chǎn)生的不確定性遠大于對表層流速和水位產(chǎn)生的不確定性,而對表層流速的不確定性大于對水位產(chǎn)生的不確定性.表明針對不同的模擬目標(biāo),邊界條件所產(chǎn)生的不確定性空間分布和影響程度也不同.當(dāng)以表層流速為輸出目標(biāo)時:表層流速受到外部輸入條件的不確定性影響介于以水位與水齡為輸出目標(biāo)時的影響之間;沿岸區(qū)域的表層流速受到外部輸入條件的不確定性影響最大.

圖4 模擬表層流速的均值, 5%分位值, 95%分位值和方差的空間分布Fig.4 Mean, 5%, 95% and variance of the simulated flow velocity in the surface layer

3.2 邊界條件敏感性分析

3.2.1 以水位為輸出目標(biāo) 對 200組外部輸入條件組合和輸出目標(biāo)水位值進行逐步秩回歸分析,得出決定系數(shù)R2為0.92(表2),表明該回歸分析是可行的.太湖初始水位 SRRC2值最大(85.73%),表明太湖初始水位對水位模擬結(jié)果最敏感,不確定性貢獻率占絕對主要的地位.其次是風(fēng)速的SRRC2值(5.25%),表明相對于其他3個邊界條件,風(fēng)速對水位的模擬結(jié)果比較敏感.出入湖流量和風(fēng)向的SRRC2值均小于1%,表明這兩個邊界條件對水位模擬結(jié)果不確定性的貢獻率較小,敏感性較差.4個外部輸入條件對模擬水位結(jié)果不確定性的貢獻率排序為:太湖初始水位(85.73%)>風(fēng)速(5.25%)>風(fēng)向(0.86%)>出入湖流量(0.10%).綜上所述,以太湖水位為模擬目標(biāo)時,初始水位對模擬結(jié)果的不確定性貢獻率最顯著,是最敏感的外部輸入條件,而風(fēng)速、風(fēng)向和出入湖流量對水位敏感性較差.

表2 4個外部輸入條件對水位和水齡不確定性的SRRC2均值和標(biāo)準(zhǔn)方差Table 2 Mean and standard deviation of SRRC2values of four input conditions on the uncertainty of the simulated water level and water age

3.2.2 以水齡為輸出目標(biāo) 對200組外部輸入條件組合和輸出目標(biāo)水齡值進行逐步秩回歸分析,得出決定系數(shù)R2為0.92(表2),表明該回歸分析是可行的.太湖初始水位的 SRRC2值最大為66.125%,表明初始水位對水齡模擬結(jié)果最敏感,不確定性的貢獻顯著.風(fēng)向和出入湖流量的SRRC2值分別為15.00%和8.37%,表明3個邊界條件中風(fēng)向?qū)λg模擬結(jié)果比較敏感,其次是出入湖流量,而風(fēng)速的 SRRC2值只有 3.00%,表明風(fēng)速對水齡模擬結(jié)果不確定性的貢獻較小,敏感性較差.4個外部輸入條件對模擬水齡結(jié)果不確定性的貢獻率排序為:太湖初始水位(66.125%)>風(fēng) 向 (15.00%)>出 入 湖 流 量(8.37%%)>風(fēng)速(3.00%).綜上所述,以水齡為模擬目標(biāo)時,初始水位對結(jié)果不確定性的貢獻率最顯著,是最敏感的外部輸入條件,其次為風(fēng)向和出入湖流量,而風(fēng)速對模擬結(jié)果的貢獻率較小,即敏感性較差.

圖5 4個外部輸入條件對各層流速的不確定性的SRRC2均值和標(biāo)準(zhǔn)方差Fig.5 Mean and stand deviation of SRRC2of four input conditions on the uncertainty of the simulated flow velocity in different layers

3.2.3 以流速為輸出目標(biāo) 對 200組外部輸入條件組合及輸出目標(biāo)表層流速進行逐步秩回歸分析,決定系數(shù)R2為0.754,表明逐步回歸分析是可行的(圖 5).風(fēng)速貢獻率為 58.70%,表明風(fēng)速對表層流速模擬結(jié)果很敏感,風(fēng)速的不確定性影響顯著.其次風(fēng)向的SRRC2值為13.52%,表明風(fēng)向也是對表層流速變化較敏感的一個邊界條件.風(fēng)速和風(fēng)向?qū)μ韺恿魉俨淮_定性的貢獻率占絕對主導(dǎo)地位.而出入湖流量和初始水位對太湖表層流速的不確定性貢獻率均小于 2%,表明其敏感性較差.4個外部輸入條件對表層流速不確定性的貢獻率排序為:風(fēng)速(58.70%)>風(fēng)向(13.52%)>太湖初始水位(1.93%)≈出入湖流量(1.27%).綜上所述:以表層流速為輸出目標(biāo)時,風(fēng)速對模擬結(jié)果不確定性的貢獻率最顯著,是最敏感的邊界條件,其次是風(fēng)向,而初始水位和出入湖流量的敏感性均較差.

4個外部輸入條件對垂向各層流速模擬結(jié)果的不確定性貢獻率的平均值和標(biāo)準(zhǔn)方差如圖5所示,各層的決定系數(shù)R2均大于0.7,表明外部輸入條件與各層流速之間的回歸分析是可行的.從圖中可知,風(fēng)速對各層流速不確定性的貢獻率都最大,均在55%～60%之間.其次是風(fēng)向,不確定性貢獻率均在 10%～15%之間.而初始水位和出入湖流量對各層流速不確定性的貢獻率都較小,均在 1%～5%之間.表明外部輸入條件對各層流速的不確定性影響規(guī)律相似,貢獻率排序均為風(fēng)速(55%～60%)>風(fēng)向(10%～15%)>初始水位≈出入湖流量(1%～5%).由圖5可知,外部輸入條件對流速不確定性的貢獻率在垂向上差距很小.原因可能是太湖屬于典型的淺水湖泊,其特點是風(fēng)生流,由于水淺的特性使得風(fēng)場對流速的影響在垂向上差別很小.由圖5(b)可知風(fēng)速和風(fēng)向?qū)γ繉恿魉俨淮_定性貢獻率的方差很大,表明這兩個邊界條件對各層流速的不確定性影響在空間上的變化較大.

綜上所述:對于大型淺水湖泊,外部輸入條件對各層流速的不確定性影響規(guī)律相似,風(fēng)速對各層流速不確定性的貢獻率最顯著,是最敏感的外部輸入條件,其次是風(fēng)向,而初始水位和出入湖流量的敏感性較差.并且風(fēng)速和風(fēng)向?qū)α魉俨淮_定性的貢獻率空間差異性很大.

4 結(jié)論

4.1 對于大型淺水湖泊的水動力模擬過程,外部輸入條件的不確定性對模擬結(jié)果(以水位、水齡以及流速表征)產(chǎn)生影響.對于不同的輸出目標(biāo),水齡受外部輸入條件的不確定性影響最大,其次是垂向平均流速,水位受到的不確定性影響最小.以水位為輸出目標(biāo)時,湖泊水位較低區(qū)域受外部輸入條件的不確定性影響較大;以水齡為輸出目標(biāo)時,在湖泊入湖河道附近水動力擾動較強的區(qū)域受到的不確定性最大;而以流速為輸出目標(biāo)時,湖泊沿岸區(qū)域受到的不確定性影響最大,且不確定性影響在垂向上差別很小.

4.2 針對不同的輸出目標(biāo),4個重要外部輸入條件對模擬結(jié)果產(chǎn)生的不確定性貢獻率也不相同.初始水位對全湖的水位和水齡模擬結(jié)果產(chǎn)生決定性影響,貢獻率分別達到85.73%和66.125%,而對垂向平均流速不確定性影響的貢獻率只有3%;風(fēng)速對垂向平均流速模擬結(jié)果產(chǎn)生決定性影響,貢獻率達到 55%～60%,而對水位和水齡的不確定性貢獻率分別只有5.25%和3.00%.

[1] Zhao X, Shen Z Y, Xiong M, et al. Key uncertainty sources analysis of

water quality model using the first order error method [J]. International

Journal of Environmental Science and Technology, 2011,8(1):137-148. [2] Lindenschmidt K E, Fleischbein K, Baborowski M. Structural uncertainty in a river water quality modelling system [J]. Ecological Modelling, 2007,204(3/4):289-300.

[3] Blumberg A F, Georgas N. Quantifying uncertainty in estuarine and coastal ocean circulation modeling [J]. Journal of Hydraulic Engineering, 2008,134(4):403-415.

[4] 秦伯強.太湖生態(tài)與環(huán)境若干問題的研究進展及其展望 [J].湖泊科學(xué), 2009,21(4):445-455.

[5] 王 鵬,王勝艷,郝少盼,等.模擬擾動條件下太湖沉積物的再懸浮特征 [J]. 水科學(xué)進展, 2010,21(3):309-404.

[6] 丁艷青,朱廣偉,秦伯強,等.波浪擾動對太湖底泥磷釋放影響模擬 [J]. 水科學(xué)進展, 2011,22(2):273-278.

[7] 余 紅,沈珍瑤.三峽水庫大寧河流域非點源污染參數(shù)的不確定性分析 [J]. 中國環(huán)境科學(xué), 2007,27(4):554-558.

[8] Mckay M D, Beckman R J, Conover W J. A comparison of three

methods for selecting values of input variables in the analysis of

output from a computer code [J]. Technometrics, 2000,42(1):55-61. [9] 劉紀(jì)濤,劉 飛,張為華.基于拉丁超立方抽樣及響應(yīng)面的結(jié)構(gòu)模糊分析 [J]. 機械強度, 2011,33(1):73-76.

[10] 李一平,逄 勇,羅瀲蔥.波流作用下太湖水體懸浮物輸運實驗及模擬 [J]. 水科學(xué)進展, 2009,20(5):701-709.

[11] Li Y P, Acharya K, Mark C S, et al. Spationtemporal patterns in nutrient loads, nutrient concentrations, and algal biomass in Lake Taihu, China [J]. Lake and Reservoir Management, 2011,27(4):298-309.

[12] Li Y P, Acharya K, Yu Z B. Modeling impacts of Yangtze River water transfer on water ages in Lake Taihu, China [J]. Ecological Engineering, 2011,37(2):325-334.

[13] 施小青,吳吉春,姜蓓蕾,等.基于LHS方法的地下水流模型不確定性分析 [J]. 水文地質(zhì)工程地質(zhì), 2009,2:1-6.

[14] Manache G, Melching C S. Sensitivity analysis of a water-quality modeling using Latin Hypercube Sampling [J]. Journal of Water Resources Planning and Management, 2004, 130(3):232-242.

[15] Pan F, Ye M, Zhu J T, et al. Incorporating layer- and local-scale heterogeneities in numerical simulation of unsaturated flow and tracer transport [J]. Journal of Contaminant Hydrology, 2009, 103(3/4):194-205.

[16] Pan F, Zhu J T, Ye M, et al. Sensitivity analysis of unsaturated flow and contaminant transport with correlated parameters [J]. Journal of Hydrology, 2011,397(3/4):238-249.

[17] 李一平,唐春燕,余鐘波,等.大型淺水湖泊水動力模型不確定性

和敏感性分析 [J]. 水科學(xué)進展, 2012,23:271- 277.

Uncertainty and sensitivity analysis of input conditions in large shallow lake hydrodynamic model.

LI Yi-ping1,2,3*,

QIU Li1,2, TANG Chun-yan1,2, BU Min-sheng4, TIAN Wei5, YU Zhong-bo3, Kumud Acharya3(1.Key Laboratory of Integrated Regulation and Resources Development of Shallow Lakes, Ministry of Education, Hohai University, Nanjing 210098, China；2.College of Environment, Hohai University, Nanjing 210098, China；3.State Key Laboratory of Hydrology Water Resources and Hydraulic Engineering, Nanjing 210098, China；4.College of Mechanical and Electrical Engineering, Hohai University, Nanjing 210098, China；5.Hydrology and Water Resources Investigation Bureau of Jiangsu Province, Nanjing 210029, China). China Environmental Science, 2014,34(2)：410～416

Uncertainty and sensitivity analysis of four important input conditions on the Environmental Hydrodynamic Fluid Code (EFDC) model results (i.e., water level, water age and currents) was investigated for a large shallow lake, Lake Taihu, China. The four input conditions included three boundary conditions (i.e., inflow/outflow, wind speed, wind direction) and an initial condition (i.e., initial water level). The Latin Hypercube sampling (LHS) as a global sensitivity method was used to estimate the uncertainty and sensitivity from the four input conditions. The results showed that uncertainties in the hydrodynamic process existed due to the uncertainties of model input conditions. Among the four input conditions, the initial water level was the most sensitive factor for the simulated water level and water age with the uncertainty contributions of 85.73% and 66.125% respectively, while it had barely 3% contributions to vertical averaged velocity. Wind speed played a significant role in the uncertainty of the velocity in the surface layer with a sensitivity coefficient of 58.70%, while it only had 5.25% and 3.00% contributions to the simulated water level and water age, respectively. Additionally, there was a similar impact of the four input conditions on the uncertainty of velocities in different layers. The four input conditions’ contributions to the velocities were as follows: wind speed (55%～60%) > wind direction (10%～15%) > initial water level ≈ inflow/outflow (1%～5%). Thus, the results provided reliable information forthe model prediction of large shallow lakes like Lake Taihu. For different output targets, improving the precision of the input conditions with priority can efficiently enhance the precision of the hydrodynamic model.

input conditions；uncertainty analysis；sensitive analysis；environmental fluid dynamic code (EFDC)；large shallow lake

X143

：A

：1000-6923(2014)02-0410-07

李一平(1978-),男,湖北荊州人,副教授,博士,主要從事河湖富營養(yǎng)化機理及水環(huán)境數(shù)學(xué)模型研究.發(fā)表論文50余篇.

2013-06-08

江蘇省高校“青藍工程”;河海大學(xué)創(chuàng)新人才計劃;國家“973”項目(2010CB951101);國家自然科學(xué)基金(51379061);江蘇省自然科學(xué)基金(BK20131370)

* 責(zé)任作者, 副教授, liyiping@hhu.edu.cn