徐長江 徐高洪 戴明龍 張冬冬

摘要：洞庭湖湖區(qū)的水文情勢變化直接影響到區(qū)域洪水災(zāi)害防治、水資源利用、水環(huán)境保護(hù)和水生態(tài)安全維護(hù)，意義重大?；趯?shí)測資料分析了三峽水庫運(yùn)行前后洞庭湖區(qū)水文情勢變化，受三口分流量減少和來水偏枯等綜合因素影響，2003～2016年8～11月洞庭湖入、出湖水量較1981～2002年分別減少26%和23.7%，湖區(qū)水位下降0.76～1.27 m。建立了長江與洞庭湖一、二維耦合水動力模型，模擬計(jì)算了三峽水庫蓄水期初設(shè)調(diào)度方案、優(yōu)化調(diào)度方案和規(guī)程調(diào)度方案下長江干流及洞庭湖區(qū)的水文過程，3種方案對應(yīng)的城陵磯站蓄水期多年旬平均流量分別減少1 220，928，900 m?3/s，湖區(qū)鹿角站蓄水期多年旬平均水位分別降低1.47，1.23，1.20 m。實(shí)際調(diào)度流量減小812 m?3/s，水位降低1.14 m，表明實(shí)際調(diào)度最優(yōu)。從不同調(diào)度方案比較來看，實(shí)際調(diào)度提前了起蓄時間，減緩了對洞庭湖區(qū)水文情勢的不利影響。

關(guān)鍵詞：水文情勢; 蓄水期; 水動力模型; 三峽水庫; 洞庭湖

中圖法分類號：P333文獻(xiàn)標(biāo)志碼： ADOI：10.16232/j.cnki.1001-4179.2019.02.002

1研究背景

洞庭湖位于長江中游荊江河段南岸，是吞吐長江的通江湖泊，南納湘、資、沅、澧四水來水，北承長江松滋、太平、藕池三口分流，湖區(qū)水系錯綜紊亂[1]，入湖徑流經(jīng)湖區(qū)調(diào)蓄后，再由城陵磯注入長江，江湖關(guān)系復(fù)雜。洞庭湖水文情勢受長江干流與湖區(qū)水系支流來水的雙重影響，如2006，2009，2011年，受三口分流減少和四水來水偏少的影響，湖區(qū)發(fā)生了大面積干旱，致使湖區(qū)生態(tài)、生活用水出現(xiàn)一定程度的困難，引起了社會各界的廣泛關(guān)注[2]。

江湖關(guān)系變化是關(guān)系到長江及洞庭湖治理、開發(fā)與保護(hù)的關(guān)鍵性科學(xué)問題，特別是三峽水庫運(yùn)行以來江湖關(guān)系的變化受到學(xué)者的廣泛關(guān)注[3]。眾多學(xué)者基于實(shí)測資料，從水情水沙時空演變[4-8]、長江與洞庭湖水體交換能力[9-10]、槽蓄特性[11]等角度分析了新條件下的江湖關(guān)系變化。江湖關(guān)系受到氣候變化以及人類活動的多要素影響，三峽水庫運(yùn)行調(diào)度作為影響因素之一備受關(guān)注，孫思瑞等[12]基于BP神經(jīng)網(wǎng)絡(luò)，分析了三峽水庫不同調(diào)度方案對洞庭湖出口水位的影響;李景保等[13]通過運(yùn)用三峽水庫運(yùn)行前后洞庭湖區(qū)實(shí)測原型水文資料，分析了三峽水庫不同調(diào)度方式對洞庭湖區(qū)典型年的水情影響。本文在已有研究基礎(chǔ)上，重點(diǎn)研究蓄水期三峽水庫不同調(diào)度方案對湖區(qū)水位以及出湖水量的影響，通過構(gòu)建水動力學(xué)模型定量評估水庫蓄水的影響，分析不同蓄水方案情景下湖區(qū)水位及出湖水量的變化規(guī)律，為減輕水庫蓄水對湖區(qū)生產(chǎn)、生活、生態(tài)等用水的影響提供可靠的技術(shù)支撐，為長江經(jīng)濟(jì)帶城市規(guī)劃建設(shè)、水資源優(yōu)化配置決策提供科學(xué)可靠的參考依據(jù)。

2區(qū)域概況

洞庭湖位于東經(jīng)111°14′～113°10′，北緯28°30′～30°23′，為我國第二大淡水湖。洞庭湖水系主要由湘江、資水、沅江、澧水四大水系和長江松滋口、太平口、藕池口三口分流水系組成，另有汨羅江、新墻河等入?yún)R。

洞庭湖區(qū)是指荊江河段以南，湘、資、沅、澧四水尾閭控制站以下，跨湘、鄂兩省的廣大平原、湖泊水網(wǎng)區(qū)，總面積20 109 km?其中天然湖泊面積約2 625 km?洪道面積1 418 km?受堤防保護(hù)面積16 066 km?2。

洞庭湖區(qū)的地勢西高東低，分成東洞庭湖、南洞庭湖和西洞庭湖（由目平湖、七里湖組成），自西向東形成一個傾斜的水面。

3模型的構(gòu)建

3.1一、二維耦合水動力數(shù)學(xué)模型

長江干流、荊江三口分流河道以及四水尾閭的水流模擬采用MIKE11一維水動力數(shù)學(xué)模型進(jìn)行，洞庭湖湖區(qū)水流模擬采用MIKE21二維水動力數(shù)學(xué)模型，長江干流與洞庭湖區(qū)一、二維耦合通過MIKE FLOOD模型實(shí)現(xiàn)[14]。耦合模型涉及到的一維河道和二維區(qū)域模型采用標(biāo)準(zhǔn)連接形式。模型地形采用2011年長江干流河道和洞庭湖區(qū)地形資料。模型概化示意見圖1。

長江與洞庭湖一、二維耦合模型中，上邊界條件為宜昌站的流量，下邊界條件為螺山站的水位流量綜合線，洞庭四水、汨羅江及清江主要控制水文站的流量作為點(diǎn)源匯入模型。螺山站的水位流量關(guān)系非常復(fù)雜，影響因素較多，本次為了開展不同方案間的比較研究，采用近年的綜合線。二維模型的上游邊界條件由一維模型給出，而一維模型的下游邊界由二維模型提供，2個模型的數(shù)值求解過程交替進(jìn)行，在耦合邊界上傳遞計(jì)算結(jié)果，實(shí)現(xiàn)耦合。

湖區(qū)內(nèi)的降雨徑流過程采用MIKE 21中降雨產(chǎn)流模塊進(jìn)行模擬。模塊中的輸入為湖區(qū)鹿角、南咀、小河咀、營田以及自治局各站的實(shí)測逐日降水過程，采用泰森多邊形的方法插值形成湖面的逐日降水，蒸發(fā)采用湖區(qū)的逐月潛在蒸散發(fā)時間序列。

為準(zhǔn)確模擬洞庭湖湖區(qū)的水動力變化過程，模型采用三角形網(wǎng)格對湖區(qū)進(jìn)行網(wǎng)格劃分，網(wǎng)格總數(shù)為23 436個。對于不規(guī)則網(wǎng)格，最大網(wǎng)格面積一般不宜超過0.1 km?2。綜合考慮洞庭湖區(qū)面積、模擬精度、計(jì)算時間及軟件性能等因素的限制，本次計(jì)算采用的最大網(wǎng)格面積不超過0.08 km?一般網(wǎng)格面積控制在0.015～0.02 km?2左右，重要地區(qū)、河道及其它地形變化劇烈的區(qū)域，計(jì)算網(wǎng)格適當(dāng)加密。

3.2模型率定

模型參數(shù)的率定主要是河道和湖泊糙率的率定，以合理反映河道和湖泊的阻力。河道阻力不但與河道自身的物理特性有關(guān)，還與其上、下游的水位、流量、水面坡降等因素有關(guān)，就長江中游荊江和洞庭湖區(qū)而言，不同的洪水組成、量級及水位漲落率均會對水流阻力產(chǎn)生影響，因此進(jìn)行糙率率定時必須選擇具有一定代表性的典型洪水過程。綜合分析，選擇近期2008，2011，2012 年洪水過程作為率定的依據(jù)，3 a的洪水過程具有不同的洪水組成和特征，代表性較好。2008年洪水屬于典型的四水遭遇型洪水，2011年屬于該區(qū)域典型的枯水年份，2012年洪水為長江中游區(qū)域性較大洪水。

由于計(jì)算區(qū)域內(nèi)不同河段、洞庭湖不同湖區(qū)的河床邊界組成差異較大，對各區(qū)域進(jìn)行糙率率定時，在對河道及湖泊特性分析的基礎(chǔ)上，將長江中游干流河段分為枝城-沙市-新廠-監(jiān)利-城陵磯-螺山5個河段，松滋河系、虎渡河系、藕池河系、澧水洪道分為22個河段，洞庭湖區(qū)及四水尾閭地區(qū)分為目平湖、南洞庭湖西側(cè)、南洞庭湖東側(cè)、東洞庭湖南側(cè)、東洞庭湖北側(cè)、湘江尾閭、資水尾閭、沅江尾閭 8 塊分別按照地形等高線進(jìn)行糙率率定。

3.3模型驗(yàn)證與精度評定

為量化評估模型的模擬精度，采用Nash-Sutcliffe系數(shù)[15]和相對誤差RE[16]2個指標(biāo)評價流量過程的模擬精度。

（1） Nash-Sutcliffe系數(shù)ENS。用以衡量模擬流量過程和實(shí)測流量過程間的擬合度，ENS值越接近于1，模擬流量越接近實(shí)測流量，公式為

ENS=1-ni=1（Qpi-Qoi）?2ni=1（Qoi-Qo）?2（1）

（2） 相對誤差RE。用以評價模擬流量和實(shí)測流量的差值與實(shí)測流量之間的百分比，能反映模擬的可信程度，公式為

RE=ni=1（Qpi-Qoi）ni=1Qoi×100%（2）

式中，Q0為實(shí)測流量;Qp為模擬流量;Q0為實(shí)測平均流量;n為實(shí)測個數(shù)。

參考Moriasi等[17]研究中采用的模型精度評價等級標(biāo)準(zhǔn)（表1），定量評價本模型的模擬精度。

根據(jù)率定的參數(shù)，模擬2013～2014年長江干流、洞庭湖湖區(qū)及出口各站點(diǎn)的逐年流量和水位過程，并與實(shí)測流量、水位過程進(jìn)行對比驗(yàn)證，模擬效果見表2和圖2，其中水位是基于1985國家高級基準(zhǔn)。由表2可以看出，除管家鋪站為良好外，其他各個站點(diǎn)模擬流量過程與實(shí)測過程擬合程度均為優(yōu)秀，滿足模型模擬精度要求。由圖2可以看出，湖區(qū)各站點(diǎn)的水位模擬過程與實(shí)測過程擬合較好，高低水位及過程走勢比較一致，南咀、小河咀和鹿角站的平均水位絕對誤差分別為0.03，0.06，0.04 m。

4洞庭湖區(qū)水文情勢變化統(tǒng)計(jì)特征

以三峽水庫運(yùn)用為時間節(jié)點(diǎn)，分運(yùn)行前（1981～2002年）和運(yùn)行后（2003～2016年）2個階段，洞庭湖區(qū)代表水位站選擇南咀、小河咀和鹿角站，出湖代表水文站選擇城陵磯站，根據(jù)各站的實(shí)測資料，對比分析三峽水庫建庫前后湖區(qū)入湖水量、出湖水量和水位的變化特征。結(jié)合長江上中游梯級水庫蓄水調(diào)度運(yùn)行實(shí)踐，研究時段采用8～11月。

4.1入湖徑流變化

洞庭湖入湖水量主要來自湘、資、沅、澧四水和荊江三口分流。采用湘江湘潭站、資水桃江站、沅江桃源站、澧水石門站合成流量代表洞庭湖四水入湖水量，松滋河新江口站和沙道觀站、虎渡河彌陀寺站、藕池河康家港站和管家鋪站合成流量代表荊江三口入湖水量。洞庭湖入湖徑流（不含未控區(qū)間入流，下同）不同時段徑流量及其變化統(tǒng)計(jì)見表3。

與1981～2002年相比，2003～2016年8～11月多年平均入湖徑流量偏少210.6億m?3，減幅26.6%，其中三口分流來水偏少132.3億m?3，洞庭四水來水偏少78.3億m?3，說明總體上入湖徑流偏少是三口分流量減小和湖區(qū)水系來水偏枯共同影響而致。

與1981～2002年相比，2003～2016年多年平均情況下，荊江三口分流來水、洞庭四水來水和入湖徑流在8，9，10月均呈偏少現(xiàn)象，尤以10月最為突出，多年平均入湖徑流偏少71.8億m?3，減幅45.6%，其中荊江三口減少39.7億m?3，洞庭湖四水減少32.1億m?3。

4.2出湖徑流變化

采用城陵磯站流量代表洞庭湖出湖入江水量。城陵磯站不同時段徑流量及其變化見表4。較1981～2002年，2003～2016年8～11月平均出湖流量偏少236億m?3，減幅23.7%，8～11月各月均減少，以10月減少百分比最大。結(jié)合表3分析，10月入、出湖水量的差值在三峽工程運(yùn)行前后基本相當(dāng)。

4.3湖區(qū)水位變化

洞庭湖湖區(qū)由東、南、西洞庭組成，出口在東洞庭的城陵磯，本次以南咀、小河咀、鹿角水位站為湖區(qū)代表站，城陵磯水文站為洞庭湖出口代表站，根據(jù)各站實(shí)測水位資料，分析洞庭湖區(qū)的水位變化規(guī)律。

以2003年為時間分界點(diǎn)，分別統(tǒng)計(jì)1981～2002年以及2003～2016年各代表站8～11月各個時段月平均水位變化，統(tǒng)計(jì)結(jié)果見表5?？梢钥闯?，與1981～2002年相比，各站月平均水位在不同時期均有一定的下降，南咀、小河咀、鹿角和城陵磯站2003～2016年8～11月多年平均水位分別下降0.82，0.76，1.24，1.27 m;從湖區(qū)各站點(diǎn)變化幅度來看，靠近出湖點(diǎn)處的城陵磯站變化幅度最大，小河咀站和南咀站變化幅度相當(dāng);各月水位變化中，各站均以10月變化幅度最大，城陵磯站2003～2016年10月平均水位較之1981～2002年下降2.24 m。

5蓄水期洞庭湖區(qū)水文情勢變化研究

5.1不同蓄水方案下宜昌站流量過程

三峽水庫以上流域已建成若干座具有較強(qiáng)調(diào)節(jié)能力的大型水庫，鑒于水庫群的影響較為復(fù)雜[18]，且洞庭湖入湖徑流的重要組成之荊江三口分流就位于三峽壩下游河段，故本次重點(diǎn)研究三峽水庫蓄水對洞庭湖水文情勢的影響，宜昌站流量還原時主要考慮三峽水庫的調(diào)蓄影響。

2008年和2009年三峽工程曾進(jìn)行過2次試驗(yàn)性蓄水，由于蓄水目標(biāo)及當(dāng)年蓄水期水情等因素，最高蓄水位分別為172.74 m和171.41 m?？紤]到2008年和2009年三峽水庫并沒有蓄滿，實(shí)際調(diào)度的蓄水量會較其他蓄水方案偏小，影響不同調(diào)度方案間的比較分析，因此，本次選取2010～2016年作為分析時段。根據(jù)三峽水庫2010～2016年的實(shí)際調(diào)度運(yùn)行資料，由三峽水庫壩前水位、庫容曲線以及出庫流量，采用水量平衡法反推入庫流量，再采用馬斯京根法將入庫流量過程由清溪場演算到宜昌站，得宜昌站天然流量過程。

由于不同調(diào)度方案主要從9月以后才開始蓄水，本次研究時段選為9～11月。在實(shí)際蓄水調(diào)度中，批復(fù)的三峽試驗(yàn)性蓄水實(shí)施計(jì)劃中的蓄水方式為：三峽水庫9月上旬可在承接8月下旬防洪調(diào)度運(yùn)用水位基礎(chǔ)上逐漸上浮水位，上浮期間控制下泄流量滿足中下游各方用水需求。依據(jù)三峽水庫初設(shè)調(diào)度方案、優(yōu)化調(diào)度方案及規(guī)程調(diào)度方案對三峽水庫進(jìn)行2010～2016年蓄水模擬調(diào)度，得到不同情境下的宜昌站9～11月逐日流量過程，與宜昌站還原后的流量過程相比，統(tǒng)計(jì)影響特征如表6所示。

由于初設(shè)調(diào)度方案、優(yōu)化調(diào)度方案和規(guī)程調(diào)度方案的蓄水起始時間分別為10月1日、9月15日和9月10日，因此以上方案對宜昌站9月份流量無影響或有部分旬時段的影響。

總體而言，實(shí)際調(diào)度對宜昌站各旬流量的影響要小于其它幾種調(diào)度方案。據(jù)統(tǒng)計(jì)，實(shí)際調(diào)度、初設(shè)調(diào)度方案、優(yōu)化調(diào)度方案和規(guī)程調(diào)度方案的宜昌站多年旬平均流量較天然情況分別減少3 160，4 740，3 580，3 420 m?3/s。從9～11月，實(shí)際調(diào)度反映出較其它調(diào)度方案相對更為均勻的蓄水過程，對各旬的流量影響相對更緩和，多年平均最大影響出現(xiàn)在9月中旬，旬流量減少4 300 m?3/s，10月底已基本蓄滿。初設(shè)調(diào)度方案、優(yōu)化調(diào)度方案和規(guī)程調(diào)度方案對旬平均流量的減少超過5 000 m?3/s的情形均有2個旬，旬平均流量的最大減小值分別為9 920，8 670，9 090 m?3/s，且多在11月中旬蓄滿。長江中下游在11月已進(jìn)入枯水期，為了避免水庫蓄水對下游水資源開發(fā)利用的不利影響，上游水庫宜在汛后盡早蓄滿。

5.2不同蓄水方案對湖區(qū)水位的影響

在固化其他邊界條件的情況下，將宜昌站2010～2016年還原后的天然流量過程和三峽水庫初設(shè)調(diào)度方案、優(yōu)化調(diào)度方案以及規(guī)程調(diào)度方案下的宜昌站流量過程作為上游邊界條件，分別模擬洞庭湖湖區(qū)典型站點(diǎn)水位以及出湖站點(diǎn)流量過程。

根據(jù)模擬結(jié)果，統(tǒng)計(jì)出9～11月不同情境下的湖區(qū)各站點(diǎn)旬平均水位的特征值，分析三峽水庫蓄水對洞庭湖湖區(qū)的水位影響。

不同情境下洞庭湖湖區(qū)各站的旬平均水位影響見圖3。由圖可以看出，總體而言，實(shí)際調(diào)度對水位的影響要小于其它幾種調(diào)度方案。據(jù)統(tǒng)計(jì)，實(shí)際調(diào)度使洞庭湖區(qū)各站蓄水期旬平均水位下降0.24～1.14 m，其次是規(guī)程調(diào)度方案，各站蓄水期旬平均水位下降0.33～1.20 m，優(yōu)化調(diào)度方案下各站蓄水期旬平均水位下降0.35～1.23 m，初設(shè)調(diào)度方案下各站蓄水期旬平均水位下降0.36～1.47 m。

在相同調(diào)度方案下，分析湖區(qū)不同位置處的水位變化幅度，處于東洞庭湖的鹿角站變化幅度最大，處于西洞庭湖的南咀站和處于南洞庭湖的小河咀站變化幅度相當(dāng)，均小于鹿角站。

以鹿角站為典型，統(tǒng)計(jì)分析不同調(diào)度情境下各旬的水位特征變化，見表7。

由表7可以看出，實(shí)際調(diào)度、初設(shè)調(diào)度方案、優(yōu)化調(diào)度方案和規(guī)程調(diào)度方案下的鹿角站蓄水期多年旬平均水位分別較天然情況下降1.14，1.47，1.23，1.20 m，實(shí)際調(diào)度對水位的影響相對最小?？傮w上，各個調(diào)度方案對鹿角旬水位的影響主要集中在9～10月下旬。從9月至10月，實(shí)際調(diào)度使鹿角站的水位出現(xiàn)一定的下降，下降過程較其它調(diào)度方案更趨緩和，旬平均水位最大下降值出現(xiàn)在10月上旬，降低1.51 m，9月上旬影響相對較小，下降0.43 m。初設(shè)調(diào)度方案、優(yōu)化調(diào)度方案和規(guī)程調(diào)度方案對旬平均水位的影響相對較大，旬平均水位最大下降值均超過了2 m，分別為2.83，2.49，2.07 m，并有多個旬平均水位下降超過1.5 m;且多年平均情況下對水位的下降影響一直持續(xù)到了11月中旬，對于湖區(qū)枯水期的水資源開發(fā)利用和水環(huán)境保護(hù)較為不利。

5.3不同蓄水方案對出湖流量的影響

根據(jù)5.2節(jié)模擬結(jié)果，統(tǒng)計(jì)出2010～2016年9～11月不同情境下的城陵磯站出湖流量變化特征見表8，分析三峽水庫蓄水對洞庭湖出湖流量的影響。

由表6和表8可以看出，優(yōu)化調(diào)度方案自9月15日開始蓄水，規(guī)程調(diào)度方案自9月10日開始蓄水，宜昌站9月中旬平均流量在相應(yīng)調(diào)度方案下較工程蓄水前的天然情況已呈明顯下降現(xiàn)象，但城陵磯站在相應(yīng)調(diào)度方案下9月下旬平均流量才出現(xiàn)顯著下降，說明三峽水庫的蓄水調(diào)度對洞庭湖的出流影響具有一定滯后性。

考慮三峽水庫蓄水對洞庭湖出湖流量影響的滯后性，統(tǒng)計(jì)分析城陵磯站的多年旬平均流量影響值可知，總體而言，實(shí)際調(diào)度對城陵磯站的流量影響要小于其他幾種調(diào)度方案：實(shí)際調(diào)度、初設(shè)調(diào)度方案、優(yōu)化調(diào)度方案和規(guī)程調(diào)度方案的多年旬平均流量減少值分別為812，1 220，928，900 m?3/s。

6結(jié) 論

本文以洞庭湖區(qū)為重點(diǎn)研究對象，根據(jù)湖區(qū)代表水位站南咀、小河咀和鹿角站以及出湖代表站城陵磯站的實(shí)測資料，對比分析了三峽水庫建庫前后湖區(qū)入湖水量、出湖水量和水位的變化特征。并通過構(gòu)建長江干流與洞庭湖區(qū)一、二維耦合的水動力數(shù)學(xué)模型，模擬計(jì)算了不同蓄水調(diào)度方案下洞庭湖湖區(qū)的水文過程，分析了不同蓄水方案對湖區(qū)水位和出湖流量的影響，主要結(jié)論如下。

（1） 與1981～2002年相比，2003～2016年8～11月多年平均入湖徑流量偏少210.6億m?3，減幅26.6%，城陵磯站8～11月平均出湖流量偏少236億m?3，減幅23.7%，湖區(qū)各代表站月平均水位在相應(yīng)時期均有一定的下降。其中，三口分流來水偏少132.3億m?3，洞庭四水來水偏少78.3億m?3，說明總體上入湖徑流偏少是三口分流量減小和湖區(qū)水系來水偏枯共同影響而致。

（2） 在固化其他邊界條件的情況下，利用模型模擬計(jì)算和統(tǒng)計(jì)分析三峽水庫不同蓄水調(diào)度方案對宜昌站流量和湖區(qū)的水位影響。結(jié)果表明，從9月上旬至11月下旬，實(shí)際調(diào)度較其它調(diào)度方案表現(xiàn)出相對更為均勻的蓄水過程，對宜昌站各旬的流量影響相對更為平緩，繼而對洞庭湖區(qū)水位的影響也要小于其他幾種調(diào)度方案。在相同調(diào)度方案下，處于東洞庭湖的鹿角站水位變幅最大，南咀和小河咀站水位變幅相當(dāng)。

（3） 在其他條件不變情況下，城陵磯出湖水量變化主要是由荊江三口進(jìn)入洞庭湖的水量變化、湖泊水體自身調(diào)蓄量變化和干流水位變化共同影響所致。三峽水庫蓄水期不同調(diào)度方案下洞庭湖出湖流量均有所減少，實(shí)際調(diào)度對城陵磯站的流量變化影響要小于其他幾種調(diào)度方案。

（4） 三峽水庫運(yùn)行以來，發(fā)揮了巨大的防洪、發(fā)電、通航、供水等綜合效益，三峽水庫汛末蓄水是實(shí)現(xiàn)枯期補(bǔ)水的必然過程，蓄水期間對長江中游干流和洞庭湖區(qū)的水文情勢產(chǎn)生了一定影響，為了盡量減小蓄水期的影響，優(yōu)化蓄水調(diào)度非常重要。從本研究可以看出，三峽水庫提前蓄水有利于減緩對洞庭湖區(qū)水文情勢的不利影響，考慮到三峽水庫以上流域若干大型水庫也都在汛末蓄水，影響更為復(fù)雜，今后需要在確保防洪安全的前提下，進(jìn)一步優(yōu)化三峽水庫預(yù)報(bào)預(yù)蓄過程，深入開展長江上游水庫群聯(lián)合蓄水調(diào)度研究。

參考文獻(xiàn)：

[1]竇鴻身，姜加虎.洞庭湖[M].北京：中國科學(xué)技術(shù)大學(xué)出版社，2000.

[2]孫占東，黃群，姜加虎，等.洞庭湖近年干旱與三峽蓄水影響分析[J].長江流域資源與環(huán)境，2015，24（2）：251-256.

[3]于亞文.洞庭湖水沙變化過程及其對人類活動的響應(yīng)研究[D].上海：華東師范大學(xué)，2018.

[4]孫思瑞，謝平，趙江艷，等.洞庭湖三口洪峰流量和水位變異特性分析[J].湖泊科學(xué)，2018，30（3）：812-824.

[5]代穩(wěn)，王金鳳，仝雙梅，等.近60年長江荊江三口水沙變化過程及對洞庭湖的影響[J].水土保持研究，2017，24（6）：255-261.

[6]劉曉群，戴斌祥.三峽水庫運(yùn)行以來洞庭湖水文條件變化與對策[J].水利水電科技進(jìn)展，2017，37（6）：25-31.

[7]謝文君，岳翠瑩，張文.洞庭湖1996-2016年時空特征變化研究[J].水利信息化，2017（5）：32-38.

[8]渠庚，唐峰，劉小斌.荊江與洞庭湖水沙變化及影響[J].水資源與水工程學(xué)報(bào)，2007，18（3）：94-100.

[9]李景保，周永強(qiáng)，歐朝敏，等.洞庭湖與長江水體交換能力演變及對三峽水庫運(yùn)行的響應(yīng)[J].地理學(xué)報(bào)，2013，68（1）：108-117.

[10]戴雪，楊桂山，萬榮榮，等.長江通過江湖水文水動力交互對其支流水域東洞庭湖的影響[J].地理學(xué)報(bào)，2018，28（8）：1072-1084.

[11]廖小紅，朱楓，黎昔春，等.典型年洪水的洞庭湖槽蓄特性研究[J].中國農(nóng)村水利水電，2018（5）：134-137.

[12]孫思瑞，謝平，陳柯兵，等.三峽水庫蓄水期不同調(diào)度方案對洞庭湖出口水位的影響[J].長江流域資源與環(huán)境，2018，27（8）：1819-1826.

[13]李景保，張照慶，歐朝敏，等.三峽水庫不同調(diào)度方式運(yùn)行期洞庭湖區(qū)的水情響應(yīng)[J].地理學(xué)報(bào)，2011，66（9）：1251-1260.

[14]劉氚，張可能，戴明龍，等.荊江—洞庭湖河網(wǎng)一二維嵌套水動力學(xué)模型研究[J].沈陽農(nóng)業(yè)大學(xué)學(xué)報(bào)，2017，48（5）：576-583.

[15]Nash J E，Sutcliffe J V.River flow forecasting through conceptual models part I — A discussion of principles[J].Journal of Hydrology，1980，10（3）：282-290.

[16]Dole W P，Jackson D L.Relative error and variability in blood flow measurements with radiolabeled microspheres[J].American Journal of Physiology，1982，243（3）：371-378.

[17]Moriasi D N，Arnold J G，Van Liew M W，et al.Model evaluation guidelines for systematic quantification of accuracy in watershed simulations[J].Transactions of the ASABE，2007，50（3）：885-900.

[18]朱玲玲，陳劍池，袁晶.洞庭湖和鄱陽湖泥沙沖淤特征及三峽水庫對其影響[J].水科學(xué)進(jìn)展，2014，25（3）：348-357.

引用本文：徐長江，徐高洪，戴明龍，張冬冬.三峽水庫蓄水期洞庭湖區(qū)水文情勢變化研究[J].人民長江，2019，50（2）：6-12.

Study on hydrological regime variation in Dongting Lake during water storage period of Three Gorges Reservoir

XU Changjiang， XU Gaohong， DAI Minglong， ZHANG Dongdong

（Bureau of Hydrology， Changjiang Water Resources Commission， Wuhan 430010， China）

Abstract： The variation of hydrological regime in Dongting Lake would directly affect regional flood prevention， water resources utilization， water environment protection and water ecology protection， which is of great significance. Based on the measured data， the variation of the hydrological situation in Dongting Lake area before and after the Three Gorges Reservoir operation was analyzed. The results showed that the inflowing and outflowing water volume of Dongting Lake from 2003 to 2016 had decreased by 26% and 23.7% respectively and the water level in lake area dropped 0.76～1.27m from August to November， compared with the period from 1981 to 2002， due to the combined factors， such as the diversion decrease of three distributaries from Jinjiang River and the decrease of incoming water， etc. A 1D and 2D coupled hydrodynamic model of the Yangtze River and Dongting Lake was established to simulate the hydrological processes of the study area during the water storage period of Three Gorges Reservoir under the initial-designed scheduling scheme， optimal scheduling scheme and protocol scheduling scheme. The corresponding average annual flow of Chenglingji station decreased 1 220，928 m?3/s and 900 m?3/s respectively， and the average annual water level of the Lujiao Station in the lake area decreased 1.47，1.23 m and 1.20 m respectively; however， the average flow decreased 812 m?3/s and the water level decreased by 1.14 m under the actual scheduling practice， indicating that the actual scheduling is optimal. From the comparison of different scheduling schemes， the actual dispatching advanced the storage starting time and slowed down the adverse impact on the hydrological situation in Dongting Lake area.

Key words：hydrological regime; water storage period; hydrodynamic model; Three Gorges Reservoir;Dongting Lake