蔣永生，陸小明，趙林林，肖 鋒，薛 坤，盧小平

（1.河南理工大學(xué)測(cè)繪與國(guó)土信息工程學(xué)院，河南 焦作 454003；2.中國(guó)科學(xué)院南京地理與湖泊研究所，江蘇 南京 210008；3.江蘇省水文水資源勘測(cè)局，江蘇 南京 210029；4.江蘇省水文水資源勘測(cè)局揚(yáng)州分局，江蘇 揚(yáng)州 225000；5.河南省測(cè)繪工程院，河南 鄭州 450003）

1 概述

水文情勢(shì)反映了湖泊等自然水體水文要素在長(zhǎng)時(shí)間序列中的趨勢(shì)及周期變化情況，是合理對(duì)水域水資源進(jìn)行利用與管控的關(guān)鍵［1-2］。水位是湖泊水文情勢(shì)的主要特征指標(biāo)，直觀反映了湖泊總體水量的收支平衡，也是推算其他水文要素并掌握其變化過(guò)程的間接資料［3-4］。因此，研究水位的動(dòng)態(tài)變化以及水文要素的演變特征，對(duì)深入把控湖泊水文情勢(shì)和建立可持續(xù)的湖泊生態(tài)區(qū)有著重要意義。

高郵湖屬淮河流域的內(nèi)陸水域，北與洪澤湖水系相連接，南與長(zhǎng)江水系相連通，同時(shí)也是南水北調(diào)東線工程重要樞紐和蘇北地區(qū)重要的水源。高郵湖地跨江蘇省高郵市、寶應(yīng)縣、金湖縣和安徽省天長(zhǎng)市，湖水儲(chǔ)量約4.6 億m3，常年水位5～5.5 m。高郵湖屬于亞熱帶溫潤(rùn)季風(fēng)氣候，具有氣候溫和、四季分明等特點(diǎn)，年平均降水量約1010.5 mm，平均蒸發(fā)量890 mm，湖區(qū)生態(tài)環(huán)境優(yōu)美，物產(chǎn)資源豐富［5］。

近年來(lái)，高郵湖的水文情勢(shì)特別是水位發(fā)生了一定程度的改變，本文就1990—2020年高郵湖6 個(gè)水文站的水位數(shù)據(jù)綜合分析水情。主要研究高郵湖水位時(shí)間序列的年際、月際變化趨勢(shì)以及水位突變性檢驗(yàn)、周期分析，進(jìn)一步討論水文情勢(shì)影響因子等與水位變化的響應(yīng)關(guān)系，以期認(rèn)清高郵湖水文情勢(shì)不同階段變化的主要驅(qū)動(dòng)因素，從而為高郵湖水資源管理與利用、濕地生態(tài)環(huán)境保護(hù)以及南水北調(diào)東線工程通水狀態(tài)下高郵湖的水閘調(diào)度、水情預(yù)測(cè)提供參考。

2 數(shù)據(jù)與方法

2.1 相關(guān)數(shù)據(jù)

高郵湖6 個(gè)水文站包括南運(yùn)西閘（閘上游）站、南運(yùn)西閘（閘下游）站、廟溝站、崗板頭站、高郵（高）站、六閘（三）站的逐日水位數(shù)據(jù)，以及南運(yùn)西閘(閘上游)站的逐日入湖和出湖流量數(shù)據(jù)，時(shí)間序列均為1990—2020年。高郵湖氣象數(shù)據(jù)包括1990—2016年的逐日降水量、蒸散發(fā)、氣溫、風(fēng)速、濕度和水汽壓數(shù)據(jù)，以及2016—2020年的逐月平均降水量、蒸散發(fā)和氣溫資料。

2.2 研究方法

灰色關(guān)聯(lián)分析是根據(jù)因素之間發(fā)展趨勢(shì)的相似或相異程度，作為衡量因素間關(guān)聯(lián)程度的一種方法，首先確定所研究的系統(tǒng)特征的參考數(shù)列和比較數(shù)列，再對(duì)原始數(shù)列進(jìn)行無(wú)量綱化處理，并求參考數(shù)列與比較數(shù)列的灰色關(guān)聯(lián)系數(shù)，最后求解關(guān)聯(lián)度并對(duì)關(guān)聯(lián)度排序?；瑒?dòng)T檢驗(yàn)［6］通過(guò)考察兩段子樣本序列的均值是否有顯著差異來(lái)判斷突變的發(fā)生。Mann-Kendall 檢驗(yàn)［7-8］是被廣泛應(yīng)用于水文和氣象時(shí)間序列的趨勢(shì)和突變點(diǎn)分析的非參數(shù)檢驗(yàn)方法。Morlet 小波分析［9-10］通過(guò)小波變換將時(shí)間系列分解到時(shí)間頻率域內(nèi)，得到周期變化動(dòng)態(tài)和時(shí)間格局，小波方差反映波動(dòng)的能量隨尺度的分布，可確定1個(gè)時(shí)間序列中存在的主要周期，方差越大，信號(hào)在對(duì)應(yīng)尺度上的周期性特征就越顯著。累計(jì)距平法中的距平值主要用來(lái)確定時(shí)間序列內(nèi)的各水文或氣象數(shù)據(jù)，相對(duì)于該數(shù)據(jù)總體均值的高低趨勢(shì)，累計(jì)距平［11-12］即距平值的累加，其曲線呈上升趨勢(shì)，表明距平值增加，下降則表示距平值減小，從曲線的明顯上下起伏還可以診斷出發(fā)生突變的時(shí)間。

3 研究與分析

3.1 水文站相關(guān)水情

由歷年統(tǒng)計(jì)數(shù)據(jù)可知，廟溝、崗板頭、高郵（高）、六閘（三）水文站年內(nèi)水位變幅相似，均存在明顯漲落變化，其每年7月、8月是豐水時(shí)段。六閘（三）站年最高水位、年最低水位和年均水位都最小，很大程度是由于其處在高郵湖附屬連通的邵伯湖之中。南運(yùn)西閘（閘上游）、南運(yùn)西閘（閘下游）2個(gè)水文站存在直接水文聯(lián)系，月際水位都保持著基本穩(wěn)定的水文情勢(shì)，且年最低水位始終高于其他水文站，標(biāo)準(zhǔn)差較小，這主要是緣于上接以寶應(yīng)湖為主的入湖流量的緩沖。高郵湖豐水期時(shí)，為防止水位越過(guò)警戒線，向下游長(zhǎng)江渠快速排水，南運(yùn)西閘（閘上游）通常被關(guān)閉，但此時(shí)在風(fēng)速等外力因素作用下仍存在著水流倒灌現(xiàn)象。枯水期時(shí)高郵湖水位較低，為使水位不低于生態(tài)水位，開(kāi)啟南運(yùn)西閘（閘上游）從上游引進(jìn)水源，南運(yùn)西閘（閘下游）以存量形式定期為湖區(qū)水源進(jìn)行補(bǔ)水?？傊?，南運(yùn)西閘為高郵湖生態(tài)水位可持續(xù)發(fā)展發(fā)揮至關(guān)重要的作用。在南運(yùn)西閘的及時(shí)、合理調(diào)控下，位于高郵湖區(qū)中央的廟溝、崗板頭和高郵（高）站水位變化過(guò)程一致性較高，年最高水位、最低水位、平均水位及水位變幅都相差不大。

3.2 代表水文站選取

結(jié)合實(shí)地考察和水文數(shù)據(jù)分析結(jié)果，發(fā)現(xiàn)6 個(gè)水文站的水位變化在一定程度上都能對(duì)高郵湖的水情有所表征，但在對(duì)高郵湖總體水情進(jìn)行觀測(cè)研究時(shí)，仍需選取最具代表性的水文站。本研究通過(guò)灰色關(guān)聯(lián)分析研究了1990—2020年6 個(gè)水文站的平均水位序列相關(guān)程度，得出關(guān)聯(lián)系數(shù)矩陣結(jié)果：

式中，i，j=1，2，…，6，分別代表南運(yùn)西閘（閘上游）、南運(yùn)西閘（閘下游）、廟溝、崗板頭、高郵（高）、六閘（三）水文站。

由關(guān)聯(lián)系數(shù)矩陣可知，南運(yùn)西閘（閘上游）、南運(yùn)西閘（閘下游）和六閘（三）站對(duì)其余水文站的關(guān)聯(lián)性普遍較弱。廟溝站能94.84%代表崗板頭水情，96.65%代表高郵（高）水情，78.68%代表六閘（三）水情，對(duì)南運(yùn)西閘上、下游站也分別有67.47% 和62.65%的水情解釋。將其與崗板頭、高郵（高）站對(duì)比發(fā)現(xiàn)，廟溝站總體上能夠較全面反映高郵湖的水文情勢(shì)，可選取廟溝站作為高郵湖代表水文站。

3.3 水位變化特征

3.3.1 水位年際變化

從枯期、汛期和年均水位3 個(gè)角度分析高郵湖1990—2020年間水位年際變化趨勢(shì)，高郵湖枯期水位序列波動(dòng)較小，枯期水位多年平均為5.22 m，最高值為2016年的5.71 m，最低值為1994年的4.59 m。汛期水位序列波動(dòng)幅度較大，汛期多年平均水位為7.22 m，最高值為2003年的9.26 m，最低值為1994年的5.54 m。高郵湖年均水位序列多年平均為5.73 m，最高水位值為2003年的6.37 m，最低水位值為2001年的5.09 m。從高郵湖3個(gè)特征水位趨勢(shì)線得知，高郵湖1990—2020年際水位序列均呈上升趨勢(shì)，且汛期水位上升速率為每10 a 0.07 m，枯期水位上升速率為每10 a 0.18 m，平均水位上升速率為每10 a 0.14 m。

為進(jìn)一步研究高郵湖年際水位的突變性，采用滑動(dòng)T檢驗(yàn)法與Mann-Kendall 法結(jié)合，同時(shí)對(duì)高郵湖1990—2020年平均水位序列進(jìn)行突變檢驗(yàn)，顯著性水平都選取0.05。在滑動(dòng)T檢驗(yàn)法中，當(dāng)子序列長(zhǎng)度為3 時(shí)，由滑動(dòng)T檢驗(yàn)統(tǒng)計(jì)量曲線可大致判斷突變發(fā)生在2013—2015年時(shí)段；當(dāng)子序列長(zhǎng)度為5 時(shí)，可大致判斷突變發(fā)生在2013年時(shí)段；當(dāng)子序列長(zhǎng)度為10 時(shí)，可大致判斷突變發(fā)生在2001—2003年時(shí)段。

1990—2020年高郵湖年均水位序列Mann-Kendall 法突變檢驗(yàn)結(jié)果可將年均水位序列分為4 個(gè)階段。

（1）第一階段為1990—1992年，其水位平均值為5.59 m，低于1990—2020 多年平均水位（5.73 m）。1990—1992年間UF值大于0，說(shuō)明高郵湖水位上升，但未超過(guò)0.05 顯著水平，即上升趨勢(shì)不顯著。這一階段內(nèi)最高水位出現(xiàn)在1991年（6.01 m），最低水位出現(xiàn)在1992年（5.15 m），除1991年水位外，其他年份水位均低于這一階段平均水位，整體處于較低的水平。

（2）第二階段為1992—1998年，其水位平均值為5.49 m，遠(yuǎn)低于1990—2020多年平均水位。1992年UF值由大于0 變?yōu)樾∮?，表明1992—1998年期間水位轉(zhuǎn)為下降，在置信度95%的水平上并不顯著，即處于減緩的趨勢(shì)。這一階段最低水位出現(xiàn)于1994年（5.13 m）。這一時(shí)期除1996、1998年外，水位均低于此階段的平均水平，水位整體比第一階段略低。

（3）第三階段為1998—2002年，其水位平均值為5.60 m，與第1階段相當(dāng)，低于1990—2020多年平均水位。2000年UF和UB曲線在0.05 顯著性兩臨界水平直線之間出現(xiàn)交點(diǎn)，但卻不在滑動(dòng)T檢驗(yàn)法的突變時(shí)域內(nèi)，故不是突變點(diǎn)。1998年UF值從小于0 變成大于0，得知從1998年起湖泊水位由下降轉(zhuǎn)為上升，也未通過(guò)顯著性水平檢驗(yàn)，表明上升不顯著。2001年UF值從大于0 變成小于0，2001—2002年期間水位變化有逐漸減緩的趨勢(shì)。這一階段內(nèi)最低水位出現(xiàn)于2001年（5.09 m），最高水位出現(xiàn)在1998年（6.07 m），水位整體略高于第二階段。

（4）第四階段為2002—2020年，其水位平均值為5.85 m，遠(yuǎn)高于1990—2020年多年平均水位。2002年UF和UB曲線在顯著性水平內(nèi)之間出現(xiàn)交點(diǎn)，符合滑動(dòng)T檢驗(yàn)法所判斷突變時(shí)域，因此2002年即為突變點(diǎn)。在2002年UF值從小于0變成大于0，表明2002年之后，湖泊水位呈現(xiàn)上升的趨勢(shì)。2002—2015年水位上升趨勢(shì)并不顯著，2015年之后UF統(tǒng)計(jì)量超過(guò)0.05 顯著性水平，說(shuō)明水位有明顯升高的趨勢(shì)。在第4 階段內(nèi)的最高水位6.37 m 出現(xiàn)于2003年，這一時(shí)期除個(gè)別年份外，水位均高于1990—2020年多年平均水位。

綜上所述，高郵湖水位序列年際變化總體呈“漲-退-回升-漲”趨勢(shì)。另外，高郵湖的年均水位大多高于多年平均水位（5.73 m），2002年前年均水位基本低于多年平均水位，進(jìn)一步說(shuō)明2002年即為突變開(kāi)始的年份。

3.3.2 水位年內(nèi)變化

高郵湖年內(nèi)水位序列波動(dòng)變化趨勢(shì)見(jiàn)圖1。從月最高水位變化趨勢(shì)看，高郵湖水位年內(nèi)變化呈單峰形，1—5月水位相當(dāng)，從5月開(kāi)始水位持續(xù)升高，直到7月達(dá)到最高水位8.78 m，此后7—12月最高水位持續(xù)下降。從月最低水位變化趨勢(shì)看，高郵湖水位年內(nèi)變化呈平穩(wěn)態(tài)勢(shì)，月最低水位先稍有所上升，后略有下降，再稍有回升。

圖1 1990—2020年高郵湖月均水位變化

從月平均水位變化趨勢(shì)看，其中7月的水位最高（6.43 m），12月最低（5.74 m）。7—9月是年中高水位時(shí)期，12月至次年4月是年中低水位時(shí)期。春夏季氣溫回升、降水緩增，促使5—6月水位稍有增加，徑流量增加；夏季放閘，導(dǎo)致湖泊水位急速抬升，至7月水位達(dá)最高點(diǎn)。8月以后湖泊水位有所下降，冬季氣溫下降，湖泊處于低水位狀態(tài)。因此，高郵湖年內(nèi)平均水位呈枯（12、1、2、3、4月）、漲（5、6月）、豐（7、8、9月）、退（10、11月）波動(dòng)。

3.3.3 水位周期性分析

在Matlab 軟件上運(yùn)用Morlet 復(fù)小波系數(shù)方法，可對(duì)高郵湖1990—2020年標(biāo)準(zhǔn)化水位序列進(jìn)行連續(xù)小波變換，進(jìn)而分析高郵湖水位周期性變化規(guī)律。其中，小波系數(shù)實(shí)部等值線可用于反映水位序列在不同時(shí)間尺度下的周期變化及其在時(shí)間域中的分布，以此判斷高郵湖水位波動(dòng)變化的特征。水位序列總體在3～5 a、8～14 a、18～24 a 時(shí)間尺度范圍內(nèi)有波動(dòng)跡象，且水位最明顯的振蕩周期在11 a 的時(shí)間尺度附近。

小波方差在11 a 處存在最大峰值，表明11a 為第一主周期，還有2 個(gè)波峰分別對(duì)應(yīng)著23 a 和3 a 尺度處，則23 a 為第二主周期，3 a 為第三主周期。根據(jù)3 個(gè)主周期尺度對(duì)應(yīng)的小波系數(shù)實(shí)部變化過(guò)程，表明在11 a 時(shí)間尺度下，高郵湖水位序列的周期為6～7 a；在23 a 時(shí)間尺度下，水位序列的周期為14 a 左右；在3a 時(shí)間尺度下，水位序列的周期為2 a左右。

3.3.4 水位變化分析

根據(jù)相關(guān)氣象數(shù)據(jù)，通過(guò)整理分析得到降水量和蒸散發(fā)量變化趨勢(shì)。其中降水量多年平均為1067.4 mm，年降水量總體呈減緩的趨勢(shì)（P＞0.05），減緩速率為每10 a 16.6 mm；蒸散發(fā)量多年平均為891.3 mm，年蒸散發(fā)量總體呈顯著上升的趨勢(shì)（P＜0.01），上升速率為每10 a 19.62 mm。另外，根據(jù)南運(yùn)西閘（閘上游）的流量數(shù)據(jù)，統(tǒng)計(jì)得到入湖、出湖流量變化趨勢(shì)。其中年入湖流量總體呈不顯著的上升趨勢(shì)（P＞0.05），上升速率為每10 a 522.6 m3；年出湖量總體也呈不顯著的上升趨勢(shì)（P＞0.05），上升速率為每10 a 175.6 m3。

通過(guò)以上統(tǒng)計(jì)發(fā)現(xiàn)，降水量作為氣候因素之一，其變化趨勢(shì)和高郵湖水位的波動(dòng)走勢(shì)具有一定相似性，說(shuō)明降水可能是高郵湖水位變化的主要驅(qū)動(dòng)因素。為此，建立了高郵湖累積水位變化和高郵湖面降水累積距平變化曲線。通過(guò)對(duì)比分析發(fā)現(xiàn)，1990—2002年，代表水文站實(shí)測(cè)的水位累積距平與降水累積距平的變化趨勢(shì)基本一致；2003—2013年，高郵湖水位累積距平變化和降水累積距平變化幅度不盡相同；2013—2020年間，兩者甚至有較大偏差，說(shuō)明自2003年以來(lái)降水不再是高郵湖水位變化的主要驅(qū)動(dòng)因素。

氣候變化、徑流量與其他因素（人類活動(dòng)、水利工程等）綜合影響著高郵湖水文情勢(shì)。結(jié)合實(shí)地勘察研究發(fā)現(xiàn)，氣候變化是影響高郵湖水位的長(zhǎng)期因素，且除降水外，還有風(fēng)速和氣溫等通過(guò)影響蒸散發(fā)來(lái)間接影響水位變化。高郵湖徑流量包括入湖流量和出湖流量2 種，徑流量也是高郵湖水位變化的重要驅(qū)動(dòng)因素，而人類活動(dòng)和水利工程等其他因素具有突發(fā)性（不具有規(guī)律），可使湖泊水位短期內(nèi)發(fā)生劇烈變化。因此，需要進(jìn)一步將實(shí)測(cè)水位累積距平變化與降水累積距平、蒸散發(fā)累積距平、（降水-蒸發(fā)）累積距平、入湖流量累積距平、出湖流量累積距平和入湖-出湖累積距平分階段進(jìn)行相關(guān)性分析（表1）。

表1 累積水位變化與影響因子的相關(guān)系數(shù)

由表1可知，1990—2002年，實(shí)測(cè)水位累積變化與降水累積距平和降水-蒸發(fā)累積距平變化分別呈0.01 顯著性相關(guān)，與入湖流量累積距平變化呈0.01顯著性相關(guān)，說(shuō)明此時(shí)段高郵湖水位變化可能主要受氣象降水的影響。2002—2013年，實(shí)測(cè)水位累積變化與入湖流量累積距平變化呈0.01 顯著性相關(guān)，與入湖-出湖累積距平變化呈0.05 顯著性相關(guān)，與降水累積距平和降水-蒸發(fā)累積距平變化分別呈0.1顯著性相關(guān)，說(shuō)明這期間的高郵湖水位變化可能主要受入湖流量的影響。2013—2020年，實(shí)測(cè)水位累積變化與蒸發(fā)累積距平呈0.05 顯著性水平相關(guān)，與降水累積距平、降水-蒸發(fā)累積距平和出湖流量累積距平變化分別呈0.01 顯著性相關(guān)，說(shuō)明該期間水位變化受氣溫、風(fēng)速導(dǎo)致蒸散發(fā)的影響較大，也同時(shí)與其他因素（包括人類活動(dòng)、水利工程等）有關(guān)。綜上所述，1990—2020年，實(shí)測(cè)水位累積變化與降水量累積距平和降水-蒸發(fā)累積距平變化分別呈0.05顯著性相關(guān)，與入湖流量累積距平變化呈0.1顯著性相關(guān)，說(shuō)明總時(shí)段的高郵湖水位變化受降水為主導(dǎo)的多重因素影響。

為探究1990—2020年不同時(shí)段高郵湖水位變化各驅(qū)動(dòng)因素的貢獻(xiàn)率，建立高郵湖水位累積距平與降水累積距平、蒸發(fā)累積距平、入湖流量累積距平和出湖流量累積距平之間的回歸方程：

式中，ΔH為水位累積距平，m；ΔP為降水累積距平，mm；ΔE為蒸發(fā)累積距平，mm；ΔRH為入湖流量累積距平，m3/a；ΔCH為出湖流量累積距平，m3/a。

通過(guò)上述回歸方程，運(yùn)算出預(yù)測(cè)水位序列累積變化，發(fā)現(xiàn)與實(shí)測(cè)水位序列累積變化具有很好的相關(guān)性（R2=0.9787），可用此公式計(jì)算出的預(yù)測(cè)水位作為實(shí)際自然水位累積。1990—2002年，高郵湖實(shí)測(cè)水位累積差值為-0.314 m，預(yù)測(cè)累積水位差值為-0.296 m，因此認(rèn)為氣象和徑流量因素可有94.3%解釋水位的變化，而其他因素可有5.7%解釋水位的變化，以該方法類推并作出其他時(shí)段的累積水位變化各因素貢獻(xiàn)百分比（表2）。

表2 不同影響因素解釋累積水位變化比例 單位：%

由上述分析可知，1990—1992年為高郵湖總體水位波動(dòng)的第一階段，降水、蒸散發(fā)、入湖流量、出湖流量和其他因素可分別有53.6%、10.5%、20.6%、9.6%和5.7%解釋水位變化，說(shuō)明此階段水位上升的主要驅(qū)動(dòng)因素是降水量；1992—1998年間的第二階段，相比第一階段總體的蒸散發(fā)和出湖流量解釋水位變化比例有所增加，導(dǎo)致總體水位稍有減退；1998—2002年間的第三階段，相比第二階段總體的降水和入湖流量解釋水位變化比例增加，而且出湖流量解釋水位變化比例減少，導(dǎo)致總體水位再次回升；第四階段中的2002—2013年和2013—2020年水位都在總體增加，但是可以看出2002—2013年水位變動(dòng)主要原因是入湖流量引起水位變化的因素比例的增加，而2013—2020年水位變動(dòng)主要原因是其他因素（人類活動(dòng)、水利工程等）引起水位變化因素比例的增加。

4 結(jié)語(yǔ)

高郵湖水文情勢(shì)的規(guī)律對(duì)高郵湖流域綜合開(kāi)發(fā)、高郵周邊城市的規(guī)劃、工農(nóng)業(yè)生產(chǎn)及建設(shè)等都有著重要的影響。本文通過(guò)1990—2020年6 個(gè)水文站的實(shí)測(cè)日水位、南運(yùn)西閘（閘上游）的日流量數(shù)據(jù)和相關(guān)氣象數(shù)據(jù)，結(jié)合灰色關(guān)聯(lián)分析、突變檢測(cè)、Morlet 小波分析、累積距平和多元回歸統(tǒng)計(jì)分析，研究蘇北南水北調(diào)東線工程主要調(diào)蓄湖泊高郵湖的水文要素的演變特征，得出如下主要結(jié)論：

（1）高郵湖不同水文站的月際、年際水位變化情勢(shì)各有不同波動(dòng)特點(diǎn)，其中南運(yùn)西閘對(duì)湖泊水位的控制作用是高郵湖生態(tài)水位可持續(xù)的重要因素，最能代表高郵湖總體水文情勢(shì)的是廟溝站。

（2）高郵湖1990—2020年水位年際序列趨勢(shì)變化可分為4 個(gè)階段，即“漲-退-回升-漲”，利用滑動(dòng)T檢驗(yàn)結(jié)合Mann-Kendall 突變檢測(cè)得知在2002年發(fā)生突變。1990—2020年水位年內(nèi)序列趨勢(shì)變化，總體呈枯（12、1、2、3、4月）、漲（5、6月）、豐（7、8、9月）、退（10、11月）波動(dòng)。

（3）高郵湖標(biāo)準(zhǔn)化水位存在3 個(gè)不同時(shí)間尺度內(nèi)的周期變化，在11 a時(shí)間尺度下，高郵湖水位序列的周期為6～7 a；在23 a 尺度下，水位序列的周期為14 a左右；在3 a尺度下，周期為2 a左右。

（4）通過(guò)建立累積水位變化與影響因子的多元回歸統(tǒng)計(jì)關(guān)系分析發(fā)現(xiàn)，1990—2002年降水量是高郵湖水位變化的主要驅(qū)動(dòng)因素，2002—2013年入湖流量是主要驅(qū)動(dòng)因素，2013—2020年則是其他因素。