劉 云,林登榮,王 斌,包中進(jìn)

(1.浙江省水利河口研究院,浙江 杭州 310020;2.溫州市水利投資開發(fā)有限公司,浙江 溫州 325000)

排澇工程規(guī)模計(jì)算及論證一般采用簡單對(duì)比法或全部試驗(yàn)法[1-6]。簡單對(duì)比法是在確定其他影響因素的基礎(chǔ)上,對(duì)單一因素進(jìn)行逐個(gè)比較,探索較為合理的排澇規(guī)模,這種方法的代表性很差;另外,排澇規(guī)模的論證牽涉因素眾多,各種因素之間關(guān)系錯(cuò)綜復(fù)雜,簡單對(duì)比法不能綜合考慮這些因素的影響,結(jié)果的可靠性也就無法驗(yàn)證。全面試驗(yàn)法雖然涵蓋了各種影響因素,但是試驗(yàn)次數(shù)太多,工作量較大。正交試驗(yàn)法[7]則能夠解決這些問題,它是利用“正交表”科學(xué)地安排與分析多因素試驗(yàn)的方法,其主要優(yōu)點(diǎn)是能在很多試驗(yàn)方案中挑選出代表性強(qiáng)的幾個(gè)試驗(yàn)方案,用統(tǒng)計(jì)計(jì)算的方法對(duì)試驗(yàn)結(jié)果進(jìn)行分析研究,分辨出主次因素,推斷出最優(yōu)方案。采用正交試驗(yàn)法既能縮短周期、降低成本,又能全面科學(xué)地分析問題,使得結(jié)論更加可靠。正交試驗(yàn)法在醫(yī)藥、食品、化工等領(lǐng)域應(yīng)用較多且比較成熟,也有研究者將正交試驗(yàn)法應(yīng)用到機(jī)械、巖土、水利等領(lǐng)域中[8-12]。

本文將正交試驗(yàn)法應(yīng)用到溫州甌飛一期圍墾工程排澇規(guī)模優(yōu)化計(jì)算中,充分考慮閘門控制、水利調(diào)度等因素的影響,確定了影響圍區(qū)排澇能力的主次因素,在此基礎(chǔ)上進(jìn)一步優(yōu)化排澇組合,最后確定了合理的排澇規(guī)模。

1 工程概況

溫州甌飛一期圍墾工程位于甌飛灘高灘區(qū)域,南、北邊界分別是飛云江及甌江河口邊界的外延線,圍墾面積約為116km2。圍區(qū)以中央橫河為界,西片設(shè)置南1號(hào)閘、北1號(hào)閘;圍區(qū)東片東面設(shè)置東1號(hào)閘、東2號(hào)閘,南北向設(shè)置南2號(hào)閘、北2號(hào)閘及兩個(gè)調(diào)蓄湖及兩個(gè)排澇泵站,東片主體為“二橫七縱兩湖”的排澇格局。設(shè)計(jì)方案南1號(hào)和北1號(hào)閘閘室凈寬均為80m,閘底板高程-2.0m;南2號(hào)和北2號(hào)閘閘室凈寬均為60m,閘底板高程-3.0m;東1號(hào)閘和東2號(hào)閘閘室凈寬均為40m,閘底板高程-3.0m。圍區(qū)近期為功能養(yǎng)殖期;中期擬采取“分區(qū)排澇、分片解決”的排澇總體格局,上游溫瑞平原客水作為一個(gè)排水系統(tǒng),高水高排,圍區(qū)東片澇水由圍區(qū)設(shè)置的排水閘及泵站等設(shè)施排出,與客水排澇體系相互獨(dú)立;遠(yuǎn)期則考慮相互貫通[13]。圍區(qū)排澇布局見圖1[14]。圍區(qū)排澇洪水標(biāo)準(zhǔn)除中期圍區(qū)東片按10年一遇控制,其他均按照50年一遇控制,各階段各閘調(diào)度控制水位不同,見表1,控制原則為高于起調(diào)水位開閘,低于最低控制水位關(guān)閘。

圖1 甌飛一期圍墾工程圍區(qū)排澇布局

表1 各閘調(diào)度控制水位

2 模型建立

數(shù)學(xué)模型主要模擬甌飛一期圍墾規(guī)劃水系,計(jì)算不同排澇規(guī)模組合下圍區(qū)最高洪水位。由于計(jì)算區(qū)域范圍較大、水動(dòng)力條件復(fù)雜、涉水建筑物較多,并需要考慮水閘調(diào)度運(yùn)行方式等因素,因此平面二維數(shù)學(xué)模型采用基于三角形網(wǎng)格的有限體積模式[15],并根據(jù)水閘調(diào)度運(yùn)行原則設(shè)置了閘門開關(guān),當(dāng)圍區(qū)水位達(dá)到排澇控制條件時(shí),閘門開啟,低于最低控制水位時(shí)閘門關(guān)閉,水閘僅允許單向排澇出流,不考慮納潮?；痉匠倘缦?

式中:h為水深,m;u,v分別為x,y方向上的垂向平均流速分量,m/s;t為時(shí)間,s;z0為河床高程,m;g為重力加速度;f為柯氏力參數(shù),f=2ωsinφ,其中φ為緯度,ω為地球自轉(zhuǎn)速度;εx,εy分別為x,y方向的渦動(dòng)擴(kuò)散系數(shù);Wx,Wy分別為x,y方向的風(fēng)應(yīng)力分量;Cz為謝才系數(shù),Cz=h1/6/n,其中 n為糙率系數(shù)。方程(1)為水流連續(xù)方程,方程(2)和方程(3)分別為x,y方向的動(dòng)量守恒方程。

模型上游客水邊界模擬至工程臨近圍區(qū)的排澇閘,下游則模擬至外排至甌江及東海的排澇閘附近。模型計(jì)算上游客水邊界采用流量過程線;下游飛云江、甌江以及東海附近6個(gè)閘的閘下邊界則采用對(duì)應(yīng)的外海潮位過程線;圍區(qū)降水直接按照降水邊界設(shè)置。

3 排澇規(guī)模優(yōu)化計(jì)算

排澇規(guī)模優(yōu)化計(jì)算主要考慮水閘配套規(guī)模問題,在水閘設(shè)計(jì)中,閘室寬度及閘底板高程對(duì)水閘的排澇能力影響較大,因此,本文在各閘閘室寬度及閘底板高程多個(gè)水平基礎(chǔ)上對(duì)圍區(qū)西片及東片排澇進(jìn)行正交組合計(jì)算,比選出對(duì)圍區(qū)排澇水位影響較大的因素并進(jìn)行優(yōu)化,結(jié)合工程投資提出較為合理的水閘布置規(guī)模,最后用遠(yuǎn)期貫通方案對(duì)優(yōu)化方案進(jìn)行復(fù)核。

3.1 中期及遠(yuǎn)期西片客水排澇工況

中期及遠(yuǎn)期西片排澇工況排澇控制標(biāo)準(zhǔn)相同,因此計(jì)算組合亦相同。圍區(qū)西片客水排澇影響因素主要考慮南1號(hào)閘、北1號(hào)閘的閘室寬度及閘底板高程,結(jié)合工程實(shí)際,對(duì)各因素分別選取了3個(gè)水平,即采用四因素三水平進(jìn)行計(jì)算,如表2所示。

表2 排澇規(guī)模優(yōu)化計(jì)算因素水平

按照常規(guī)計(jì)算方法需要計(jì)算81組次,為提高效率,加快進(jìn)度,將正交試驗(yàn)法用于排澇優(yōu)化計(jì)算中,采用L9(34)正交試驗(yàn)表,將試驗(yàn)組次減少為9組,詳見表3,試驗(yàn)結(jié)果見表4。

表3 西片客水閘室規(guī)模正交試驗(yàn)組次

表4 西片客水閘室規(guī)模正交試驗(yàn)結(jié)果 m

根據(jù)正交設(shè)計(jì)的特性,由Ki的大小可以判斷各水平對(duì)試驗(yàn)指標(biāo)的影響大小,本試驗(yàn)統(tǒng)計(jì)指標(biāo)為圍區(qū)最高洪水位,因此Ki越小,試驗(yàn)結(jié)果越好。由極差大小可以看出南1號(hào)閘閘寬、南1號(hào)閘閘底板高程對(duì)圍區(qū)最高洪水位影響最大;由ˉK i值大小得出各因素水平3為最優(yōu)排澇組合,即閘室寬度取最大,閘底板高程取最低,該結(jié)論也符合工程實(shí)際情況。

另外,從表3可以看出,除組次5、組次7、組次9外,其余組次圍區(qū)最高洪水位均高于設(shè)計(jì)最高洪水位。從投資來看,組次5為相對(duì)較優(yōu)方案,但考慮組次4方案圍區(qū)最高洪水位與設(shè)計(jì)最高洪水位較為接近,僅相差0.04m,且由極差分析得出,南1號(hào)閘閘底板高程對(duì)圍區(qū)最高洪水位的影響僅次于南1號(hào)閘閘寬,因此根據(jù)閘室附近實(shí)際地形高程嘗試將組次4的南1號(hào)閘閘底板降為-2.0m,即比較方案(組次10)。計(jì)算結(jié)果表明,組次10圍區(qū)最高洪水位為3.74m,等同于設(shè)計(jì)最高洪水位,滿足排澇要求。因此,根據(jù)圍區(qū)西片排澇計(jì)算,推薦組次5和組次10為較優(yōu)方案。

3.2 中期及遠(yuǎn)期東片圍區(qū)排澇工況

中期及遠(yuǎn)期東片圍區(qū)排澇工況仍采用以上方法進(jìn)行計(jì)算分析,計(jì)算組次都由81組減少為9組,計(jì)算所得優(yōu)化方案為南2號(hào)、北2號(hào)閘寬40m,東1號(hào)、東2號(hào)閘寬45m,閘底板高程均為-3.5m,圍區(qū)東片閘室總凈寬170m。

3.3 遠(yuǎn)期聯(lián)合排澇工況

該方案將圍區(qū)東片及西片通過節(jié)制閘進(jìn)行貫通組合,用此方案對(duì)以上計(jì)算結(jié)果進(jìn)行復(fù)核,見表5。計(jì)算結(jié)果表明,貫通后圍區(qū)優(yōu)化方案1及2的最高洪水位分別為3.60m和3.61m,低于設(shè)計(jì)最高洪水位3.65m,滿足設(shè)計(jì)排澇標(biāo)準(zhǔn)。兩方案閘室總凈寬分別為330m和320m,均比設(shè)計(jì)總凈寬360m小,減少了工程投資。

表5 遠(yuǎn)期貫通閘室規(guī)模復(fù)核

4 結(jié) 論

a.將正交試驗(yàn)法用于排澇規(guī)模優(yōu)化計(jì)算中,減少了計(jì)算組次,提高了效率,根據(jù)極差分析法找出了圍區(qū)最高洪水位影響因素的主次順序,從而進(jìn)一步優(yōu)化了排澇組合方案。

b.通過計(jì)算推薦兩個(gè)較優(yōu)方案,兩個(gè)方案的排澇寬度均比設(shè)計(jì)方案總凈寬小,減少了工程投資,為設(shè)計(jì)及工程運(yùn)行提供了依據(jù)。

[1]張惠平,吳嬌容.衢江新區(qū)3號(hào)排澇工程規(guī)模論證方法探討[J].浙江水利科技,2009(4):64-66.

[2]梁金艷.三鋪湖排澇泵站規(guī)模核算分析[J].安徽水利水電職業(yè)技術(shù)學(xué)報(bào),2011,11(1):24-25.

[3]柏菊.沿江圩區(qū)排澇水文計(jì)算實(shí)例[J].江淮水利科技,2007(2):45-46.

[4]林錫雄,陳映.澄饒聯(lián)圍排澇計(jì)算程序[J].廣東水利水電,2000(1):38-41.

[5]徐昌杰,劉俊,陶雷.混合圩區(qū)排澇布局優(yōu)化研究[J].中國農(nóng)村水利水電,2009(12):84-86.

[6]吳立新,韓愛華.對(duì)2005年淮河入海水道排澇調(diào)度的回顧與思考[J].水利水電科技進(jìn)展,2006,26(增刊1):34-35.

[7]蘇均和.試驗(yàn)設(shè)計(jì)[M].上海:上海財(cái)經(jīng)大學(xué)出版社,2005.

[8]趙惠菊.正交試驗(yàn)優(yōu)選烴含量分析條件[J].現(xiàn)代科學(xué)儀器,2011(1):108-111.

[9]彭安華,孫旭東,王智明.基于正交試驗(yàn)法的PID控制器參數(shù)整定[J].機(jī)械科學(xué)與技術(shù),2011,30(6):1028-1032.

[10]倪恒,劉佑榮,龍治國.正交設(shè)計(jì)在滑坡敏感性分析中的應(yīng)用[J].巖石力學(xué)與工程學(xué)報(bào),2002,21(7):989-992.

[11]翟遠(yuǎn)征,王金生,蘇小四.正交試驗(yàn)法在地下水?dāng)?shù)值模擬敏感性分析中的應(yīng)用[J].工程勘察,2011(1):46-50.

[12]陳惠玲.水工試驗(yàn)設(shè)計(jì)[M].南京:河海大學(xué)出版社,1995.

[13]王靈敏,呂中明,邵偉.溫州市甌飛一期圍墾工程圍區(qū)排澇控制規(guī)劃[R].杭州:浙江省水利水電勘測設(shè)計(jì)院,2011.

[14]包中進(jìn),王斌,徐崗.溫州市甌飛一期圍墾工程水閘水工模型試驗(yàn)研究階段報(bào)告[R].杭州:浙江省水利河口研究院,2011.

[15]王敏,楊國錄,史英標(biāo),等.一、二維數(shù)學(xué)模型在飛云江河口航道建設(shè)中的運(yùn)用[J].中國農(nóng)村水利水電,2006(8):109-112.