近100年来长江与洞庭湖汇流河段水文特征演变*

肖 潇 ，毛北平，吴时强

(1:河海大学水文水资源与水利工程科学国家重点实验室，南京 210098) (2:长江水利委员会水文局长江中游水文水资源勘测局，武汉 430012) (3:南京水利科学研究院，南京 210029)

长江在荆江三口分流(调弦口1958年建闸控制)进入洞庭湖后，与湘、资、沅、澧来水经洞庭湖调蓄作用后在城陵矶与长江再度汇流，构成江湖分合、相互影响、相互制约的错综复杂的江湖关系[1]. 在长期的演变过程中，尤其在近百年来江湖整治等人类活动和气候变化影响下，江湖关系持续发生改变，长江洞庭湖的水文特征不断发生变化，影响着整个长江中下游的防洪、水资源利用等问题. 因而，长江洞庭湖的江湖关系一直是学者们研究的热点.

周柏林、付湘等[2-3]通过建立数学模型研究在三峡工程影响下城陵矶水域水位变化情况. 唐金武、丛振涛、王鸿翔等[4-6]认为城陵矶水位是洞庭湖水位的重要标志，城陵矶水位发生变化主要是三峡工程引起的. 而李景保等[7]通过对长江洞庭湖的水量交换研究认为江湖关系是互补过程，丰水期荆江三口对洞庭湖补给强，枯水期洞庭湖对长江的补给强；何征等[8]发现洞庭湖在枯水期湖泊出流减缓，丰水期三口分流减少，江湖关系复杂多变. 柯文莉等[9]试图通过研究洞庭湖水面面积与城陵矶水位的相关性，建立城陵矶水位在丰水期与洞庭湖水面面积的相关关系. Cheng等[10]分析三峡工程对洞庭湖水文变化的影响，表明洞庭湖秋季水位下降0.52～2.26 m，最低年水位日提前5～35 d. Dai等[11]建立2009-2013年三峡工程运行对洞庭湖与长江水体交换影响的水动力模型，得出三峡工程运行明显降低了洞庭湖和干流的洪峰流量和洪涝量的结论. 以上针对长江、洞庭湖江湖关系的研究虽然方法多种多样，但结合监利、城陵矶、螺山3个水文站历史实测水文资料，分析长江洞庭湖汇流河段水文特征演变及相互机理的较少.

鉴于此，本文主要以城陵矶1904年建站以来100多年的长系列水文资料，结合监利站、螺山站的长系列水文资料，系统分析汇流河段水文要素变化特征，进一步梳理河段年内分配、年际变化趋势、突变特征及变化周期等，深化对江湖关系的认识. 长江洞庭湖汇流河段经历的水利历史背景特殊，是人类积极改造自然、利用自然的发展史，系统地对长江洞庭湖汇流河段水文要素变化特征进行深入分析，对今后进一步治理和改善江湖关系意义重大.

1 研究区和资料情况

1.1 研究区概况

长江洞庭湖汇流河段[12]，跨越湖北和湖南两省交界区域，位于长江中游，属于亚热带湿润季风气候，降水充足，但降水集中，降水量年际变化较大. 监利水文站位于长江下荆江河段，上距洞庭湖汇合口约80 km，是下荆江河段重要水文控制站. 城陵矶水文站位于洞庭湖出口段，上距长江汇合口约4 km[13]，是洞庭湖出口河段水文控制站. 而螺山水文站位于长江干流，下距洞庭湖与长江汇合口约30 km[14]. 研究区具体位置见图1.

图1 长江洞庭湖汇流河段示意图Fig.1 Diagram of the Changjiang River and Lake Dongting confluence reach

1.2 资料来源与处理

本文依据的水文数据资料为长江水利委员会水文局提供的长江洞庭湖汇流河段3个水文站的实测资料，包括逐日水位数据：监利(1950-1969、1974-2017年)、城陵矶(1904-1938、1945-2017年)和螺山(1953-2017年)；以及逐日流量数据：监利(1950-1969、1975-2017年)、城陵矶(1946-2017年)、螺山(1953-2017年). 资料的主要处理方法为：1)根据三站水位、流量逐日数据分析汇流河段年内水文分布特征及年际水文变化特征，下文中水位的高程基面均采用冻结基面；2)利用M-K检验[15]分析水文数据变化趋势及突变点；3)利用Morlet小波[16]函数分析水文数据的周期性变化；4)分析时段主要是考虑1966年下荆江裁弯取直、1981年葛洲坝运行和2003年三峡水库运行进行划分.

2 长江与洞庭湖汇流河段水文要素的分布特征

2.1 年内水文要素的分布特征

2.1.1 水位 根据长江洞庭湖汇流河段监利、城陵矶、螺山三站的水位年内变化规律发现，3个水文站的水位波动特征相似(图2). 洞庭湖水位变化大致可分为4个阶段：12月至次年3月为枯水期，4-6月为涨水期，7-9月为丰水期，10-11月为退水期[16]. 监利、城陵矶、螺山三站多年月均最大值均出现在7月，分别为33.18、29.95、29.24 m；监利站多年月均最小值出现在2月，为24.47 m，而城陵矶、螺山二站多年月均最小值出现在1月，分别为19.85、18.89 m. 监利、城陵矶、螺山三站月均最大值与最小值的差值分别为8.71、10.10、10.35 m. 汇流河段的年内差值变化大，而洞庭湖水位比下荆江河段年内变化大，造成长江洞庭湖汇流后下游水位变化更明显. 监利、城陵矶、螺山三站涨水期的上升速率分别为2.00、2.24、2.22 m/月，退水期的下降速率分别为2.14、2.68、2.56 m/月，退水期下降速率快于涨水期的上升速率，涨水期和落水期的变化机理存在差异，季节性变化对洞庭湖水位的影响大于长江干流.

图2 水位的年内变化Fig.2 Annual variation of water level

长江洞庭湖汇流河段的水面比降季节性变化明显(图3). 城陵矶站与螺山站的月均水面比降在0.035‰～0.046‰之间，最大比降发生在2月，最小比降发生在8月；监利站与螺山站的月均水面比降在0.035‰～0.057‰之间，最大比降发生在1月，最小比降发生在7月；监利站与城陵矶站的月均水面比降在0.034‰～0.053‰之间，最大比降发生在1月，最小比降发生在5月. 仅涨水期城陵矶-螺山河段的水面比降会大于监利-城陵矶河段，丰枯水期及退水期均低于监利-城陵矶河段水面比降. 年均相关系数显示，城陵矶与螺山的年均相关系数R2为0.97，监利与螺山的年均相关系数R2为0.96，城陵矶与监利的年均相关系数R2为0.96. 3站相关性在涨水和退水期较稳定，此期间相关性强，在枯水期和丰水期，相关性较弱. 枯水期时，城陵矶站与监利站、城陵矶站与螺山站的水位关联性较强，但监利站与螺山站的水位关联性较弱，表明长江干流水位在枯水期受洞庭湖影响明显；丰水期时，螺山站与监利站的水位关联性较强，而螺山站与城陵矶站的水位关联性弱，表明丰水期长江城陵矶下游干流水位受洞庭湖来流的影响减弱，而受荆江来流影响增强.

图3 三站之间的水面比降及相关系数的年内变化Fig.3 Annual variation of the water surface slope and correlation coefficient

2.1.2 流量 根据长江洞庭湖汇流河段三站年内流量变化规律(图4)，将城陵矶出流与荆江来流即城陵矶站流量与监利站流量定义为汇流比，以下简称汇流比. 从图4可以看出，3-6月城陵矶站流量大于监利站，期间洞庭湖出流对下游螺山河段流量的贡献率较大，占48.9%～54.6%，汇流比在1.04～1.28之间. 该时段处于涨水期，洞庭湖出流对下游水位上涨驱动作用明显，相关性强. 从7月至次年2月监利站流量大于城陵矶站，即荆江来流对下游螺山河段流量的贡献率较大，占58.2%～69.3%，汇流比在0.48～0.84之间. 尤其是枯水期荆江来流占的比重最大，即汇流比最小.

图4 流量及其占比的年内变化Fig.4 Annual variation of flow and its percentage

汇流河段年径流量主要集中在5-9月(图5)，监利站、城陵矶站、螺山站各占其年径流总量的65.59%、73.44%、63.64%. 三站最大月径流量均出现在7月，各占年径流总量的16.50%、16.34%、16.45%. 城陵矶站最小月径流量一般出现在1月，占年径流总量的2.57%；而螺山站、监利站一般出现在2月，各占年径流总量的3.02%、3.06%. 可见枯水期经城陵矶流入长江的流量较小，螺山站流量来自长江干流.

图5径流量及其占比的年内变化Fig.5 Annual variation of runoff and its percentage

2.2 年际水文变化特征

2.2.1 水位 从长江洞庭湖汇流河段监利站、螺山站、城陵矶站历年水位时空分布变化(图6)可知，城陵矶站和螺山站水位在高洪期间关联密切，而监利站高洪时由于洞庭湖调蓄使得下游螺山站洪水明显缓和. 由于城陵矶站与螺山站之间距离较近，两者水位图谱关联性较强. 随着年份的变化，城陵矶站与螺山站低水变化明显，分界点是1980s，其整体上升显著. 监利站的低水位也在上升，不过明显的分界点是2003年后. 那么城陵矶站与螺山站低水位上升主要是由于地形的改变[17-18]，监利站的低水位上升则是由于三峡工程的调蓄作用[11,19]. 高洪期，城陵矶站发生时期相比1950s明显提前，监利站发生时期则在逐年提前.

城陵矶站和螺山站的年平均水位、最高水位、最低水位整体呈上升趋势，而监利站水位变化整体稳定. 1951-1966年相比1904-1950年，城陵矶站高、中、低水位均有所下降，变幅为0.02～0.36 m. 1967-1980年相比1951-1966年，监利站高水位上升，中、低水位下降，而城陵矶站、螺山站水位均呈上升趋势，变幅为0.17～0.58 m. 1966-1969年实施过下荆江裁弯取直工程，该工程应该是引起江湖汇流河段水位上升的主要因子. 1981-2002年相比1951-1966年，监利站、城陵矶站、螺山站水位均上升明显，变幅在0.23～0.97 m之间. 很多学者都探讨过汇流河段水位上升的原因，大多数认为是河床形态改变引起的；2003-2017年相比1981-2002年，监利站、城陵矶站、螺山站表现为低水位上升且波动幅度小，中、高水位下降，主要是三峡工程调蓄作用造成的. 从百年来看，1968-2017年相比1904-1967年，城陵矶站水位整体抬升，高、中、低水位分别抬升0.98、0.56、1.46 m，即汇流河段百年来水位特征变化明显.

图6 水位的年际变化Fig.6 Interannual variation of water level

2.2.2 流量 从长江洞庭湖江湖汇流河段监利站、城陵矶站、螺山站的年径流量及特征流量年际变化(图7)可以看出，监利站多年平均径流量、最小流量均逐渐增长，最大流量变化不明显. 城陵矶站多年平均径流量及最大流量一直呈减小趋势，即洞庭湖出流整体减小；而其最小流量在1951-1980年间整体呈减小趋势，1981-2017年间逐渐增大. 螺山站多年平均径流量、最大流量基本稳定，而最小流量呈比较稳定的增长趋势. 城陵矶站与监利站年径流量之比一直在减小，由1951-1966年的1.26已降低到2003-2017年的0.66.

图7 径流及特征流量的年际变化Fig.7 Interannual variation of runoff and characteristic flow

长江洞庭湖江湖汇流河段监利、城陵矶、螺山三站径流年际极值比(为使数据具有可比性，径流年际极值比同乘0.1)介于1.29～2.80之间，径流变差系数介于0.076～0.366之间(图8). 城陵矶站的径流年际极值比和变差系数变化较大，说明洞庭湖出流径流丰枯变化较剧烈，差异大；螺山站相对监利站而言变化较大.

图8 径流年际极值比和变差系数Fig.8 Runoff variation extreme value ratio and coefficient

2.3 水文特征参数变化的趋势特征及突变点

对长江洞庭湖江湖汇流段3个水文站的年径流量变化进行M-K突变点检验(图9)，其中UK统计量为正序列统计量，UB为逆序列统计量，0.05显著水平为95%置信区间. 在1954-2017年时间段内，监利河段年径流量变化呈先减小后增加的趋势，1970s是分界点，平均水位变化基本呈先减小后增加的趋势，1960年也出现突变点. 城陵矶站年径流量变化一直呈减小趋势，平均水位变化则呈先减小后增加的趋势，两者在1960和1973年均出现突变，年径流变化也在2008年出现明显突变. 螺山年径流量变化一直在波动，明显的突变点是2008年，平均水位变化与城陵矶站平均水位变化一致，两者也在1960年有一个突变点. 从M-K突变检验分析看出，城陵矶站与监利站的年径流量变化是相反的，也证实洞庭湖出流与下荆江来流相互影响、相互制约、此消彼涨. 两者相互改变的交界点是1970s，即荆江裁弯取直后，证明裁弯取直是造成两者改变的主要原因. 另外，城陵矶站和螺山站在2008年也出现了一个突变点，其与三峡蓄水密切相关，三峡水库开始进行175 m实验性蓄水. 城陵矶站与螺山站的平均水位变化趋势基本一致，即两者之间关联密切. 三站年径流量、平均水位变化均在1960年出现突变点，考虑到当时大规模的围湖造田及兴建排涝工程，应该是湖水库容减小及排涝引起的干流径流、水位增加.

图9 M-K趋势检验(a～c:年径流量；d～f：年平均水位)Fig.9 M-K mutation test (a-c: annual runoff; d-f: annual water level)

2.4 水文要素变化的周期特征

图10为采用Morlet小波函数对长江洞庭湖汇流河段水文序列数据进行周期性变化特征分析的成果，从数据波动周期的等高线图可知：小波变换系数实部为正值，对应径流量、水位增大；小波变化系数实部为负值，对应年径流量、水位减小. 城陵矶站和螺山站径流、水位变化周期比较一致，尤其是水位变化周期较为吻合. 监利站与城陵矶站、螺山站的周期波动差别较大. 城陵矶站、螺山站年径流量序列16、32 a尺度上的波动比较明显，分别出现了4、2个高低径流交替震荡中心；城陵矶站、螺山站水位序列16、64 a尺度上波动较为明显，分别出现了4、1个高低水位交替震荡中心. 监利站年径流量序列8～16、32 a尺度上波动明显，水位序列8、32 a尺度上波动明显.

监利站年径流量、水位平均周期为8 a，城陵矶站和螺山站的年径流量、水位平均周期为16 a. 监利站在32 a尺度上水位呈“减小-增加-减小”的趋势；城陵矶站与螺山站水位关联密切，在64 a尺度上水位均呈“减小-增加”的趋势.

图10 Morlet小波变换系数实部时频分布Fig.10 Real part time-frequently distribution from Morlet wavelet transform coefficients

3 江湖汇流河段水文相互影响机理

3.1 长江洞庭湖汇流河段水文关系变化原因探讨

1954-2017年，长江洞庭湖汇流河段受到几次较大的工程影响. 在1960s以前，洞庭湖区内涝灾害严重，造成湖区经济损失严重，为了改善内涝清空，1960s-1970s，修建了很多排涝工程，再加上此期间大规模的围湖造田，造成河道流量增加，引起汇流河段径流、水位的上涨. 1967-1969年，为了加大长江干流下泄能力，下荆江河段实施人工裁弯取直工程，该工程对江湖关系影响较大，造成长江入荆江三口分流减小，江湖汇流河段洞庭湖出流与下荆江来流汇流比发生较大改变. 1970s-1980s，撇洪渠工程的建设以及葛洲坝工程竣工通航，此期间汇流河段枯水期水位逐渐上升，众多学者[17]认为这是河段淤积引起的. 1986-1996年，实施洞庭湖一期治理，11个重点堤垸堤防巩固达标，湖区防洪标准得到提高，但也是引起河道内水位上涨的原因之一. 2003-2008年，三峡工程2003年开始蓄水，2008年完成175 m试验性蓄水，众多学者[19]均以2003年三峡工程开始蓄水为分界点来研究三峡工程的影响，本文M-K检验表明2008年才是径流的突变点，也就是2008年三峡完成175 m试验性蓄水对长江洞庭湖汇流河段的影响才更加明显.

3.2 长江洞庭湖汇流河段汇流比变化

汇流比为城陵矶站流量与监利站流量之比，图11a显示了1954-2017年4个时期汇流比此消彼涨的剧烈变化. 为研究汇流比与3个测站之间的复杂关系，以螺山站流量为横坐标，以汇流比为纵坐标，对4个不同时期的汇流比进行点绘(图11b). 结果表明：汇流比在年内整体上呈先增大再减小的趋势，1954-2017年间，汇流比最大值一直在降低，其发生时间在逐渐提前. 1967-1980年相比1954-1966年，汇流比最大值降低0.5左右. 1981-2002年相比1967-1980年，汇流比最大值发生时间提前约2个月. 即洞庭湖出流一直在减小，年内长江干流受洞庭湖出流影响明显的时间在不断缩短. 根据时间节点来看，引起汇流比降低提前主要是下荆江裁弯取直、三峡蓄水等工程引起的，导致三口入洞庭湖的流量减小[20]，洞庭湖出流相应减小.

汇流比与螺山站流量的关系则越来越集中，1954-2017年来，螺山站高洪时期汇流比值逐渐变小，表明近年来高洪时期洞庭湖出流对长江干流的影响减小；最大比值越来越集中在螺山流量为15000～25000 m3/s的区间，表明该流量下洞庭湖出流影响最大. 汇流比的变化表明洞庭湖出流对长江干流的影响越来越弱，尤其是高洪时期影响明显减小，年内影响时间也在不断缩短，仅在涨水期影响较大；相反，长江干流下荆江河段来流在不断增加，影响在加大，高洪时期影响明显加强，影响时间在不断增加. 洞庭湖出流与下荆江来流一直是此消彼涨、相互影响、相互制约的.

图11 城陵矶与监利站汇流比的变化Fig.11 The confluence ratio between Chenglingji and Jianli Stations

图12 监利(a)、城陵矶(b)、螺山(c)各站水位与螺山流量的关系Fig.12 Relationship between the water levels of Jianli (a)， Chenglingji (b), Luoshan (c) Stations and discharges of Luoshan

3.3 长江洞庭湖汇流河段水位流量关系变化

长江洞庭湖江湖汇流处水位流量关系的改变对长江中下游水文情势影响重大，研究水位流量关系的变化极为重要. 长江洞庭湖汇流河段监利、城陵矶、螺山3站水位与螺山站流量关系如图12所示，看出螺山站流量与3站的水位关联性是密切相关的，可以用螺山站流量来推算长江洞庭湖汇流河段的水位. 1954-2017年间，3站的水位流量关系整体变化基本一致，1967-1980年相比1954-1966年， 40000 m3/s以下流量级，对应城陵矶站、螺山站水位抬升，40000 m3/s以上流量级，对应监利站水位抬升、螺山站水位下降；1981-2002年相比1967-1980年，20000 m3/s以上流量级，对应监利站、城陵矶站、螺山站水位均有所下降；2003-2017年相比1981-2002年，3站水位流量关系变化明显，且水位流量关系规律性不强、不集中，线性上看20000 m3/s以下流量级，对应水位抬升，20000 m3/s以上流量级，对应水位下降. 相关数据表明：三峡工程对长江洞庭湖汇流河段水位流量关系有一定影响，改变了河道的自然流态，水位流量关系受调蓄影响较大.

4 结论

1)长江洞庭湖汇流河段水文特征年内分配有较大差异性，季节性变化明显. 枯水期汇流段下游水位与洞庭湖水位关联性强；丰水期则与荆江水位关联性强. 3-6月涨水期，洞庭湖出流较荆江来流大，其余时段下荆江来流比重大. 长江洞庭湖汇流河段水文特征年际变化表现为：城陵矶站和螺山站的年平均水位、最高水位、最低水位整体呈上升趋势，而监利站水位变化整体稳定. 监利站多年平均径流量、最小流量均逐渐增长，最大流量无明显变化；城陵矶站多年平均径流量及最大流量一直呈减小趋势，而最小流量在1951-1980年整体呈减小趋势，1981-2017年逐渐增大；螺山站多年平均径流量、最大流量基本稳定，最小流量一直稳定增长.

2)1954-2017年，长江洞庭湖汇流河段受到几次较大工程的影响. 结合M-K突变检验表明下荆江裁弯取直和三峡蓄水175 m对江湖水文特征的影响较为显著. Morlet小波分析长江洞庭湖汇流河段水文序列周期性变化特征表明：监利站年径流量、水位平均周期约为8 a，城陵矶站和螺山站的年径流量、水位平均周期约为16 a.

3)汇流比在年内整体上呈先增大再减小的趋势，1954-2017年间，汇流比的最大值一直在降低，其发生时间逐渐提前. 这些变化主要是下荆江裁弯取直、三峡蓄水等工程引起的，导致三口入洞庭湖的流量减小，洞庭湖出流相应减小.

本文从汇流河段3个水文站的相互关系结合汇流比等要素变化讨论了长江洞庭湖汇流河段的演变特性，深入研究了汇流河段水文演变特征. 但是，汇流河段水动力特性、泥沙输移等演变特征并未考虑，后期将引入相关研究分析，为江湖关系的研究提供更好的技术支撑.