径流和潮流共同作用下的福姜沙河段分流分沙不平衡性_刘高峰.pdf

收稿日期：；修回日期：基金项目：上海市科技创新行动计划（）作者简介：刘高峰（），男，江西省德兴市人，研究员，博士，主要从事河口海岸水动力科研工作。：，（）：径流和潮流共同作用下的福姜沙河段分流分沙不平衡性 刘高峰（上海河口海岸科学研究中心，上海）摘 要：分流分沙比是研究分汊河段的一个重要指标。通过收集和分析 年长江口福姜沙河段实测水文资料，并采用了三维潮流泥沙数学模型 对分流分沙进行计算研究。结果表明：上游径流增加或下游潮差减小均会导致福北水道的分流分沙不平衡度增加，航道的维护量和分流分沙不平衡度呈正相关性。福北水道仅有约 的流量难以输移 沙量，这种分流分沙比的不平衡性会导致福北水道航道严重淤积。根据研究成果，可在福姜沙目前的三座丁坝上实施加长或加高等工程措施，以降低福北水道的分流分沙不平衡性来减少福北航道的淤积量。关键词：长江口；福姜沙河段；福北水道；分流分沙比；不平衡度中图分类号：文献标志码：文章编号：（）（，）：，：；研究背景长江口福姜沙河段位于长江下游澄通河段，上起江阴长江大桥下至九龙港，全长约 ，该河段上距南京约 ，距离吴淞口约 ，处于长江口潮流界内。福姜沙河段是典型的弯曲分汊河型，在平面上呈现“两级分汊、三汊并存”的格局，形成了长江下游典型的多级分汊多级汇流的复杂河势格局。南京以下 深水航道二期工程于 年 月开工，年 月竣工验收。二期工程建设范 第 卷 第 期长 江 科 学 院 院 报 年 月 围为长江干线南通天生港区至南京新生圩港区 河段，主要建设项目为：在福姜沙水道、口岸直水道、和畅洲水道和仪征水道建设航道整治工程，将航道深度建设成为（理论深度基准面）。福姜沙水道二期工程实施后，现状条件下福姜沙河段 深水航道的维护疏浚段主要集中在福北水道，福北航道的维护疏浚量在 年为万，在 年为万。有很多学者研究了福北航道回淤量大的问题，目前这些研究主要是关注靖江边滩淤涨和周期性切割下移影响福北航道中上段航道水深条件。但是靖江边滩淤涨和周期性切割下移的时间周期是 ，所以这些研究无法解释福北航道的维护量在年内洪季和枯季也存在较大差异的问题，福北航道的维护量年内洪季和枯季的变化主要是受到上游流量变化和潮汐变化影响。对于福姜沙河段这样的多级分汊河段，各汊道分流分沙变化对汊道发育演变意义重大。图 研究河段和断面位置 河道分汊是一种广泛存在的自然现象，汊道的分流分沙变化对于主支汊冲淤兴衰的影响无疑起着关键的作用。汊道的冲刷引起河槽容积的增大与分流比的变化并非完全对应关系，主支汊河槽容积相对关系的变化才是引起分流比变化的根本原因。汊道会发生一定的兴衰与交替变化，一般而言，弯曲度较大的汊道因为分流比小、分沙比大而逐渐趋向衰弱，而弯曲度较小的汊道则相反。主汊和支汊兴衰交替是分汊河段河道演变的共同特点，大多数汊道主支汊交替周期长，是一种相对稳定的河床形态。随着河道的这种变形，各汊道分流道的分流分沙也随之改变，而分流分沙的改变，又反作用于分汊河道的演变。分汊口门处的环流强弱和悬移质垂线分布是影响汊道分沙比与分流比关系的主要因素。在径流与潮汐共同作用的感潮分汊河段，如果分汊河道的分流分沙能保持基本平衡，那么该河道容易保持长久稳定，比如长江口比较稳定的南北港河段分流分沙比基本保持在 左右波动。但是实际上受多种因素的影响，多级分汊、弯曲的河道分流分沙总是呈现出不平衡性特征。产生分流分沙不平衡性的原因很多，在径流与潮汐共同作用的分汊河段，汊道的分流、分沙比主要受地形、含沙量、上游径流和潮差等因素影响。本文以径流、潮流共同作用下的福姜沙河段分流分沙比为研究对象，并采用实测资料分析和河口、海岸和海洋 模型（，）的改进型 研究上游径流、潮汐变化以及两者相互作用对分流分沙不平衡性的影响。研究结果可为研究福姜沙河段减少航道淤积量的延长或者加高现有丁坝措施提供科学依据。研究区域和方法 研究区域福姜沙河段目前航道情况是福中双向航道（双行），航道宽度为 ，福北航道是单向航道（上行），航道宽度为（航道位置见图）。福姜沙河段整治工程研究开展以来的 年开始至 年期间，一共开展了 次大规模水文泥沙观测，观测资料中布置有福北水道（）和福中水道（）两条断面（断面位置见图）。这两个断面采用了 第 期刘高峰 径流和潮流共同作用下的福姜沙河段分流分沙不平衡性 测量断面流量，并采用了垂线法测量了代表垂线的含沙量，据此计算了分流、分沙比，这些丰富的分流分沙观测资料为本研究提供了良好的基础。研究方法 资料分析方法本文收集整理了 年的福姜沙河段的 次洪、枯季的水文测验资料，上游大通站流量最小为 ，最大为 。在内河分汊河道已经有比较规范的分流比和分沙比计算公式，实践中在潮汐河段采用潮量的方法进行统计，利用“面积包围法”计算出两涨两落闭合潮段内的总潮量，然后分别用同一河段主、支汊的潮流量与此河段的潮流量相比，从而计算得出主、支汊的分流比。分汊道河段的分流比计算式如下：；（）。（）式中：、分别为主汊和支汊分流比；、分别为主汊和支汊一个潮周期内某时刻 的流量。根据定义，实践中在潮汐分汊河段的分沙比计算式如下；（）。（）式中：、分别为主汊、支汊的分沙比；、分别为主汊和支汊在一个潮周期内某时刻 的含沙量。本文定义汊道分沙比与分流比的比值为分流分沙不平衡度，计算式为。（）该指标可以表示某一个汊道分流量和分沙量差异的相对大小，是可以表征一个汊道水沙关系的重要指标之一。分汊河道的分流分沙能保持基本平衡的时候，汊道的分流分沙不平衡度。分汊河道的分流分沙不平衡的时候，汊道的分流分沙不平衡度，如果，表明该汊道输沙量占比多于流量占比，反之表明该汊道的输沙量占比小于流量占比。潮流泥沙模型介绍三维潮流 模型是在普林斯顿大学的海洋模式（，）的基础上发展起来的一个适用于计算河口海岸区域的三维模式。年，等在此基础上发展了一个非正交坐标系下的 模式，这个模式垂向上采用 坐标系（可以确保不同水深有一样的分层数的一种坐标系），水平方向上采用非正交曲线网格，并采用“”网格差分格式。等对于该模型的物质输运方面作了一些改进，提高了物质平流输送的计算精度。刘高峰在 模型加入了泥沙输运模型，在泥沙模块的处理中采用总变差非增（，）格式来计算泥沙沉降过程，考虑了盐度和含沙量对絮凝沉降速度的影响，把模型发展成为 模型，并且使用该模型计算了径流对长江口泥沙变化的影响和净输移模式。本文即采用该三维潮流泥沙 模型建立了一个福姜沙河段局部数学模型。采用大小模型嵌套计算，大模型包含整个长江口，上边界是长江大通站，小模型的下边界取在徐六泾，由大模型计算徐六泾潮位过程，再插值到小模型下游边界上。小模型中垂线水深分 层，网格编号 方向节点数为，方向的节点数为，一共 个网格节点（见图）。福姜沙河段的进行局部加密，最小网格间距为 ，能够较好地贴合岸线，更好地模拟工程。5 2 05 3 05 4 05 5 05 6 0 x/k m3 5 3 03 5 3 53 5 4 03 5 4 53 5 5 0y/k m福姜沙民主沙长青沙图 福姜沙河段数学模型网格图 采用 年 月水文测验资料对模型进行了验证，潮位和流速的验证结果见图，图中有研究河段 个水文站点的潮位验证，以及实测水文资料中、和 三个站点流速流向验证，可见模型很好地验证了潮位和流速过程，可以用来开展该河段的研究工作。实测资料分析 福北水道的分流分沙比收集 年在福姜沙河段开展的 次水文测验资料，统计福北和福中水道的分流比和分沙比见表。根据表 数据可知，福北水道的分流比为 ，均值为，分沙比为 长江科学院院报 年 图 福姜沙河段数学模型验证 ，均值为，分流比小于分沙比。福中水道的分流比 ，均值为，分沙比为 ，均值为，分流比大于分沙比。在枯季期间，福北水道如皋水道用仅有的过流量承担了 输沙量，因此出现了分流比和分沙比的不协调、不匹配的现象。表 福姜沙水道分流比、分沙比 年月流量（）分沙比 分流比 不平衡度福北断面（）福中断面（）福北断面（）福中断面（）福北断面（）福中断面（）特别是在洪季期间，福北水道的平均分流比为，分沙比为，福中水道的平均分流比为，分沙比为，因此，在洪季福北水道仅有约 的过流量承担了 的输沙量输移任务，因此这种分流比和分沙比的不匹配和不平衡性在洪季更加突出。由于航道整治工程的缩窄，福北水道的宽度并不大，泥沙的输移主要集中在航道范围内向下输移。所以，进入福北水道的泥沙主要沿着航道向下输移，因此这种分流分沙的不平衡性是导致福北水道的输沙压力比较大、易造成航道的淤积的重要原因之一。径流和分流分沙比的关系在径流与潮汐共同作用的感潮分汊河段，汊道的分流分沙比受上游径流、地形、含沙量和潮汐等因素影响，本节主要研究上游径流和潮汐对于分流分沙的影响。对收集到的福姜沙河段 年期间的 次水文测验资料进行了整理分析，绘制了分流分沙比和上游流量的关系见图。1.01.52.02.53.03.54.04.55.00.2 50.3 00.3 50.4 00.4 50.5 50.6 00.6 50.7 00.7 50.8 0分流比(分沙比)(a)福北水道分沙比分流比分流比拟合线分沙比拟合线分流比(分沙比)(b)福中水道流量/(m3s-1)1.01.52.02.53.03.54.04.55.0流量/(m3s-1)图 福北和福中水道断面分流比、分沙比和流量关系 由图 及表 数据经过计算可知：（）分流比变化特征。福北断面和福中断面的分流比约为 ，大通流量从 增加到 时，福北水道分流比由增加到。福中水道分流比从减小为。随着上游大通流量增加，福北水道断面的落潮分流比呈现增加趋势，即呈现线性正相关关系，该结果和已有的研究结论一致。（）分沙比变化特征。总体上福北水道断面和福中水道断面的分沙比平均是 ，大通流量从 增加到 时，福北水道落潮分沙比由增加到。福中水道落潮分沙比为减小为。随着上游流量增加，福北水道断面的落潮分沙比呈现增加趋势，呈现线性正相关关系，而福中水道断面分沙比呈现逐渐减小的趋势，呈现线性负相关关系。第 期刘高峰 径流和潮流共同作用下的福姜沙河段分流分沙不平衡性 总之，随着上游流量的增加，福北水道落潮分流比和分沙比呈现线性增加趋势，而福中水道落潮分流比和分沙比呈现线性减小趋势。径流和分流分沙不平衡度及航道疏浚维护量关系 如果一个汊道分沙比约等于分流比，此时不平衡度，即来一定量的沙就有相应比例水量将其输移走，河段易趋于冲淤平衡；如果该汊道分沙比大于分流比，即，上游来一定量的沙却没有相应比例水量将其输移走，那么该汊道易趋于淤积；如果该汊道分沙比小于分流比，即，上游来的水量超过需要输移的沙量，则泥沙易于输移走，该汊道易趋于冲刷。根据式（）对实测资料进行计算得到福北水道的分流分沙不平衡度 ，点绘上游大通流量和分流分沙不平衡度 关系见图。由图 可见：福北水道的分流分沙不平衡度 和上游流量呈正相关性，即上游大通流量增加会导致不平衡度 增加。虽然目前实测数据尚不够多，但资料中有洪季和枯季数据，数据也有一定的代表性。1.01.52.02.53.03.54.04.55.01.01.21.4不平衡度趋势线不平衡度流量/(万m3s-1)1.31.10.9图 福北水道分流分沙不平衡度和上游流量的关系 年 月 年 月福北水道航道维护水深为；年 月 年 月福北水道航道维护水深为；年 月 年 月实施 航道疏浚工程（基建期），福北水道航道维护水深从原有的 逐步增深至；年 月到现在为止福北水道航道维护水深为。年福北水道航道维护量在（）万 之间，航道的维护压力比较大，而且现场航道维护时洪季疏浚量往往大于枯季疏浚量。收集福北水道的航道维护量和上游流量资料，绘制两者的关系见图。上游大通站流量影响到达长江口区域耗时是 ，一般上游流量变化并不会在 之内发生巨大变化。采用绘图的航道疏浚量数据是多日积累的疏浚量，而不是上游流量在时间上完全同步疏浚产生的疏浚量，因此采用航道进行疏浚时段的流量来表征流量的影响基本可行，这里面包含了泥沙输移和淤积的滞后性影响。由图 可见上游流量和维护量存在一定正相关关系，这个研究结果与王建军等关于上游洪季径流量增加会导致福北水道淤积量增加的研究结论相符。图 福北水道航道疏浚维护量和上游流量的关系 因为现场维护疏浚量受到疏浚船只、航道监测图频率和疏浚工艺等多方面因素影响，所以图 中散点比较分散，但是总体上流量越大维护量就越大这个基本规律比较清晰，其原因是上游大通站的流量越大，水流输移进入本河段的泥沙就越多，且分流分沙的不平衡性就越大，淤积增强导致航道淤积增加，从而导致航道的维护量增加。数学模型试验研究 径流对分流分沙不平衡性影响通过数学模型 计算了上游流量分别为 万、万、万 而下游分别为大潮和小潮情况下的福北断面（）分流、分沙比变化。计算统计福北水道断面在各种流量情况下大潮、小潮的分流分沙比，见图。由图 可知：福北水道断面分沙比大于分流比，且随流量的增加，福北水道分沙比增加的幅度要大于分流比增加幅度，大潮、小潮均遵从该规律。福北水道和福中水道是福姜沙河段的第二级分汊水道。福北水道虽是支汊，但其河道深泓线与上游主流衔接更加平顺，福中水道是主汊，但其河道深泓线和主流有一定夹角，因此和上游主流方向的弯曲率更大。与上游主流平顺衔接支汊的分流比和分沙比会随着流量增大而增加，而与主流有一定夹角、与上游主流衔接不够平顺的主汊分流比和分沙比会随着流量增大而减小，这符合水流具有“大水趋直，小水坐弯”的规律。长江科学院院报 年 1234560.2 40.2 60.2 80.3 00.3 20.3 41234560.2 20.2 60.3 00.3 40.3 8(a)大潮分流比分沙比拟合线(分流比)拟合线(分沙比)(b)小潮流量/(万m3s-1)流量/(万m3s-1)()()图 模型计算大潮、小潮情况下福北水道分流分沙比和流量关系 泥沙输移的直接动力是潮流动力，输沙会出现和输水类似的变化趋势，同时由于水流挟沙力一般是水流流速的 次方关系，所以分沙比变化并非等同于分流比的线性变化。由于福北水道分沙比大于分流比，据式（）计算分流分沙不平衡度 值，绘制模型计算的不平衡度 值与上游流量的关系见图。1234561.0 01.0 21.0 41.0 61.0 81.1 01.1 21.1 41.1 6不平衡度趋势线不平衡度流量/(万m3s-1)图 模型计算福北水道分流分沙不平衡度和上游流量的关系 由图 可以看出，福北水道的分流分沙比不平衡度 和上游流量呈正相关关系。随着上游流量的增加，福北水道的分流分沙不平衡度 值也增加，即表明上游流量越大分流比和分沙比的不匹配度就越大，所以实测资料和模型计算均反映出相同变化规律。由图 和图 可以看出，流量增加会导致更多泥沙进入福北水道，而福北水道又比较弯曲，分流比分沙比不平衡度增加，水流难以把泥沙输移下去，因此该航道维护量就越大。潮汐对分流分沙不平衡性的影响通过数学模型 计算了福北水道的分流分沙比和潮差的关系见图。2.52.62.72.82.93.03.13.23.30.2 40.2 60.2 80.3 00.3 20.3 40.80.91.01.11.21.31.41.51.60.2 20.2 60.3 00.3 40.3 8潮差/m分流比分沙比拟合线(分流比)拟合线(分沙比)(b)小潮(a)大潮潮差/m()()图 大潮和小潮期福北水道分流比、分沙比和潮差关系 由图 可见，福北水道的分沙比在小潮期大于大潮期，故小潮期间分流分沙比的不平衡性越大，越容易导致航道淤积。因此，潮差对于分流比和分沙比也存在影响，但是对分流比影响程度小于对分沙比的影响程度。径潮流相互作用对分流分沙不平衡性影响长江口福姜沙河段位于潮流界以下，是径流和潮汐共同作用明显的河段。潮流界不仅反映潮流动力对河流水体及其荷载的泥沙阻滞、截留作用，从而形成感潮河段与之适应的河床淤积形态，也反映了径潮动力作用的空间界限，是潮汐动力的重要参考指标。径流单向指向下游，主要起着增强落潮流、削弱涨潮流的作用。潮流一直往复循环，落潮时潮流和径流同方向，从而增强落潮动力，涨潮时潮流和径流反方向，从而弱化涨潮动力。通过数学模型 计算了上游流量分别为 万、万、万 的全潮过程，分析发现随着上游径流量的增加，低潮位和高潮位都不断抬高，而且低潮位抬高幅度大于高潮位抬高幅度，因此导致潮差减小，而且潮波的相位也略有不断被推后趋势。统计大潮和小潮情况下分流分沙不平衡度 和潮差变化，见图。在小潮时的分流分沙不平衡度 值普遍大于在大潮时的 值，潮差在 时 值在，潮差在 时 值在。总体上 值随着潮差的减小而增加，且 值在小潮增加速率大于大潮增加速率。由此推断，在大潮向中潮变化过程中，由于分流分沙不平衡度增加，福北航道淤积量趋于增加。在中潮向小潮变化的过程中，由于 第 期刘高峰 径流和潮流共同作用下的福姜沙河段分流分沙不平衡性 2.52.62.72.82.93.03.13.23.31.0 01.0 21.0 41.0 61.0 81.1 01.1 21.1 41.1 6不平衡度拟合线不平衡度0.80.91.01.11.21.31.41.51.61.1 01.2 01.3 0不平衡度潮差/m潮差/m(a)大潮(b)小潮1.3 51.2 51.1 51.0 5图 大小潮期间分流分沙不平衡度和潮差的关系 分流分沙不平衡度增加速率更大，福北航道淤积量增幅也越大。通过以上分析可知，洪季福北水道仅约 的过流量承担了约 的输沙量，分流分沙比的不平衡性出现显著不平衡。这种分流分沙比不平衡性在洪季也更加突出，现场航道维护量洪季大于枯季。从逻辑推理上即可知分流分沙比不平衡性增加会增加航道淤积，不平衡性减少会减少航道淤积。所以，解决福北航道淤积量大的一种思路是研究在福姜沙目前已建的三座丁坝上实施加长或加高等工程措施，使福北水道分流比增加、分沙比减少，以降低福北水道的分流分沙不平衡性。结 论（）洪季福姜沙河段的福北水道仅有的约 过流量承担了约 输沙量，分流、分沙比的不平衡性是导致福北水道航道易淤积、维护量大的原因之一。（）上游径流的增加会导致福北水道的分流分沙不平衡度越大，即上游流量越大，分流比和分沙比的不匹配度就越大，航道也就越容易淤积。大流量是易造成福北水道的航道淤积原因之一。（）福北水道的分流分沙不平衡度随着潮差的减小而增加，且其在小潮增加的速率大于在大潮增加的速率。（）在大潮向中潮变化的过程中，由于分流分沙不平衡度增加，福北航道维护量趋于增加。在中潮向小潮变化的过程中，由于分流分沙不平衡度增加速率更大，福北航道维护量增加的幅度也越大。参考文献：曹民雄，夏云峰，马启南长江福姜沙水道河床演变分析人民长江，（）：付 桂长江南京以下 深水航道治理工程福姜沙水道整治工程效果对比分析水运工程，（）：刘高峰，王统泽福姜沙河段输沙、边滩输移特征和动力机制研究水运工程，（）：张世钊，夏云峰，徐华，等 长江福姜沙水道深水航道整治工程分流分沙影响试验研究中国海洋工程学会第十八届中国海洋岸工程学术讨论会论文集北京：海洋出版社，：黄卫东，刘同宦，栾华龙水沙条件变化下长江下游典型江心洲汊道形态演变及趋势长江科学院院报，（）：胡向阳，潘庆燊长江中下游分汊河段的整治长江科学院院报，（）：余新明 河道冲淤情况下分流分沙模式初步研究武汉：武汉大学，：李琳琳，余锡平分汊河道分沙的三维数值模型清华大学学报（自然科学版），（）：刘传杰，李 保，李昌生，等 年长江口南北港水文 要 素 变 化 分 析 水 利 水 电 快 报，（）：，谢鉴衡 河床演变及整治北京：中国水利水电出版社，长江水利委员会水文局长江口水文水资源勘测局 年洪季长江口水域主要汊道分流分沙比测验成果报告上海：长江水利委员会水文局长江口水文水资源勘测局，韦立新，曹 双，张涛涛，等感潮分汊河道平均分流、分沙比确定方法探讨人民长江，（）：杨万伦，道付海，栾华龙，等长江口洪季南北槽落潮分流分沙比观测研究华东师范大学学报（自然科学版），（）：，：，：，（）：刘高峰 长江口水沙运动及三维泥沙模型研究 上海：华东师范大学，：刘高峰 长江口泥沙净输移模式初步研究浙江省水利河口研究院第九届全国泥沙基本理论研究学术讨论会论文集北京：中国水利水电出版社：姜宁林，张长宽，陈永平上游径流变化对长江澄通河段汊道分流比的影响人民长江，（）：王建军，杨云平，申 霞，等长江下游福姜沙河段边心滩演变及对航槽冲淤影响研究应用基础与工程科学学报，（）：（编辑：占学军）长江科学院院报 年