径流和潮流共同作用下的福姜沙河段分流分沙不平衡性_刘高峰.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图 福北水道分流分沙不平衡度和上游流量的关系 年 月 年 月福北水道航道维护水深为；年 月 年 月福北水道航道维护水深为；年 月 年 月实施 航道疏浚工程（基建期），福北水道航道维护水深从原有的 逐步增深至；年 月到现在为止福北水道航道维护水深为。年福北水道航道维护量在（）万 之间，航道的维护压力比较大，而且现场航道维护时洪季疏浚量往往大于枯季疏浚量。收集福北水道的航道维护量和上游流量资料，绘制两者的关系见图。上游大通站流量影响到达长江口区域耗时是 ，一般上游流量变化并不会在 之内发生巨大变化。采用绘图的航道疏浚量数据是多日积累的疏浚量，而不是上游流量在时间上完全同步疏浚产生的疏浚量，因此采用航道进行疏浚时段的流量来表征流量的影响基本可行，这里面包含了泥沙输移和淤积的滞后性影响。由图 可见上游流量和维护量存在一定正相关关系，这个研究结果与王建军等关于上游洪季径流量增加会导致福北水道淤积量增加的研究结