黄河口滨海区冲淤演变与潮流不对称

2021-08-10 14:58杨洋陈沈良徐丛亮
海洋学报 2021年6期
关键词:输沙黄河口涨潮

杨洋,陈沈良*,徐丛亮

(1.华东师范大学 河口海岸学国家重点实验室,上海 200241;2.黄河河口海岸科学研究所,山东 东营 257091)

1 引言

在河口近岸浅水区域,潮波传播受岸线地形、底床摩擦和径流等的影响而普遍存在变形,导致浅水分潮的产生。浅水分潮与天文分潮相互作用或者与部分天文潮之间的相互作用会产生潮流不对称[1],另外余流以及余流和分潮的相互作用亦会影响潮流不对称[2]。三角洲滨海区的流场成分包括潮流、径流、风生流、波浪破碎产生的沿岸流、密度流等,其中,潮流具有周期往复性质,其他流在一定时间尺度内具有定常流性质,形成余流。单纯潮汐作用下的潮流不对称的影响因素主要包括不规则的地形岸线和不均匀的底床摩擦,而实际流场还包括余流对潮流不对称的影响。关于河口及近岸潮不对称的研究,Friedrichs和Aubrey[3]提出用M4分潮与M2分潮的振幅比来计算潮汐不对称的大小,用两者的相位差来定义其方向。Nidzieko[4]将统计学中的样本偏度计算引入潮汐不对称的量化中,并成功应用于加利福尼亚海岸的3个不同形态的河口。Song等[5]扩展了该方法,认为只有满足频率条件的分潮组合才会引起不对称。Gong等[2]基于Nidzieko的方法,定量了潮汐不对称性、潮流历时不对称性以及潮流流速不对称性,并研究了珠江三角洲黄茅海河口的这3种不对称性以及相互之间的关系。

涨落潮不对称现象是河口滨海区流场的重要特征,是导致区域涨落潮输沙不对称的重要动力因素[6],流速和历时的差异使得涨落潮输沙能力不同,从而产生潮周期净输沙[7]。由于泥沙输运效率与流速的高次幂成正比,因此即使较小的流速差异也会导致较大的泥沙输运差异。在河口滨海区,除了余流影响泥沙的净输运,涨落潮不对称也起到重要作用[3−4],泥沙的净输运又直接影响着区域的冲淤变化。因此,开展涨落潮不对称的研究有助于深入了解泥沙输运特征,深化对三角洲地貌演变机制的认识,对于河口治理、滨海湿地保护以及海洋工程建设与维护等都有重要的意义。

黄河三角洲普遍存在“波浪掀沙,潮流输沙”模式[8],且受区域气候条件的影响,黄河入海泥沙具有“夏储冬输”的显著特征:夏季入海泥沙在河口近岸区域沉积,在冬季较强水动力条件下,沉积物再悬浮并向外输运[9]。黄河口滨海区流场主要受地形、风场、径流、水体密度和地转偏向力等因素的影响,其中,潮流受渤海潮波系统和地形相互作用的影响,风生流主要受当地风场的控制,径流则主要在入海口附近。近年来黄河入海径流锐减,而且自2007年,汛期黄河自然出汊向北入海后,径流对莱州湾海域的影响变得更小,因此研究区域主要受潮流和风生流的影响,两者对泥沙的输运和再分配起到重要的作用。

本文基于Delft 3D模型建立黄河口滨海区数值模型进行水动力模拟,获取可靠的流场数据,计算涨落潮流速不对称的空间分布,既考虑了单纯潮汐作用下的潮流不对称,也考虑了流场中包含了风生流和波生流后的潮流不对称,并结合沉积物粒度特征和趋势分析、悬沙输运等,综合分析涨落潮流速不对称对黄河口滨海区冲淤变化的影响。

2 研究区域

黄河是我国第二大河,历史上流路摆动、尾闾出汊频繁,河口位置不断变迁。自1855年黄河回归渤海以来,共形成了多个三角洲叶瓣,其中最近形成的一个叶瓣体是从1976年发育至今的清水沟三角洲叶瓣。期间,口门亦发生多次变动,1996年汛前因胜利油田开发油区需要在清8断面附近实施人工改汊,使得口门向北偏移;2007年汛期在汊3断面以下300 m上下的位置自然出汊,主流向北入海;从2008年至今流路未发生显著摆动[10]。

由于无潮点的存在使黄河三角洲海域潮汐特征较为复杂,潮流多为半日潮型往复流,潮汐涨落和潮流涨落之间存在时间差。由于M2分潮的潮流椭圆长轴近似平行于浅海等深线,故各区域潮流的往复流流向基本平行于等深线走向[11]。从黄河海港无潮点向两侧,潮型依次为正规全日潮−不正规全日潮−不正规半日潮,其中,现行黄河口海域的潮汐类型为不规则全日潮。无潮点两侧涨落潮具有摆动特征:一侧涨潮,则另一侧落潮。

渤海属于相对封闭的内海,外海的波浪不易传到黄河口滨海区,因而该海区的波浪以本地风生浪为主,波浪与风场具有高度的相关性,波浪随风场的变换具有明显的季节性特征。黄河三角洲海域波浪具有风区短、波高且周期短的特点[12]。

黄河三角洲是由入海泥沙快速堆积而成的,沉积物较为松散,在不断的落淤−再悬浮−输运−沉积的过程中发生筛选和再分配。受搬运和沉积过程的动力条件控制,沉积物的量和粒度分布在不同区域存在明显差异。

3 资料与方法

3.1 数据资料

本文使用的数据主要包括水深数据、沉积物粒度和水文泥沙定点观测数据。水深数据为2007年和2015年黄河三角洲滨海区水下地形130条断面实测数据,水深基准面均为黄海基准面。沉积物粒度数据来源于水深测量时同期获得的表层沉积物样品测试资料。水文泥沙观测资料来源于2006年11月(SA1、SA2站)、2009年 7月(A1、A2站)现行河口附近区域的同步观测(图1)。

图1 研究区水文泥沙观测站位和底质沉积物位置Fig.1 Locations of hydrological observation stations and sediment sampling in the study area

3.2 悬沙通量机制分解

悬沙通量机制分解在世界大河河口和海湾均有广泛应用[13–16],包括黄河三角洲海域[17]及莱州湾[18]。

悬沙通量机制分解是将悬沙通量分解成多个动力项,以探究各个动力因子对悬沙输移贡献的大小。潮周期平均单宽悬沙输运通量计算公式为[19]

式中,h表示水深;u表示瞬时流速;c表示瞬时含沙量;T 为潮周期;上划线“—”表示垂向平均;“′”表示垂向偏差;下脚0表示潮周期平均;下脚t表示潮变化。其中,T1为欧拉余流输沙项;T2为潮汐与潮流相关项,即斯托克斯漂移输沙项;T3为潮位与含沙量的潮变化相关项;T4为潮流与潮变化含沙量相关项,通常被称为潮泵输沙项;T5为潮汐、潮流和潮变化含沙量三者相关项;T6为垂向环流输沙项;T7为潮流和潮平均含沙量的垂向切变与潮位的相关项;T8为余流和潮变化含沙量的垂向切变与潮位的相关项;T9为潮流和潮变化含沙量两者的垂向切变相关项;T10为潮流和潮变化含沙量两者的垂向切变与潮汐相关项。根据输沙的驱动力,可以把输沙项分为两类:平流输沙项和潮流输沙项,其中T1、T2、T6和T7属于平流输沙项,其余6项属于潮流输沙项。

3.3 潮流不对称的计算

Nidzikeo[4]提出,用水位时间导数的偏度来定量潮汐不对称性的大小。之后,他又将类似的方法推广到了潮流不对称性的量化[20]。计算公式如下:

式中,当x为水位的时间导数时,γ表示涨落潮历时不对称;当x为流速时,γ表示涨落潮流速不对称。γ大于0表示涨潮主导,γ小于0表示落潮主导。

3.4 数值模型

3.4.1 计算区域与网格

本文采用Delft 3D数值模型,模型计算区域包含整个渤海海域,为了保证研究区域的计算精度并提高计算效率,网格分辨率从黄河三角洲近岸的约200 m逐渐过渡至外海边界附近的约3 km,共有511×191个格点。模型网格与黄河三角洲岸线拟合较好,同时也保证了网格的正交性和平滑性,提高了模式计算的精度和稳定性。垂向采用σ坐标,等分为5层。基于CFL数(Courant-Friedrichs-Lewy Number),模型选取时间步长为 1 min。

3.4.2 边界条件设置

模型底边界采用的近岸地形由2007年黄河三角洲滨海区水下地形实测数据插值所得。

模型共设有两个开边界,其中外海开边界采用潮汐调和常数以确定水位,由 M2、S2、K1、O1、N2、K2、P1和Q1共8个主要分潮驱动,各分潮的调和常数由美国俄勒冈州立大学开发的 OTIS(OSU Tidal Prediction Software)潮汐预报模型获得,所使用的数据为TPXO全球大洋潮波模式结果。另一个开边界为河流开边界,边界设在现行河口,边界水动力条件以通量形式给定,采用的是利津水文站的月平均径流量,数据来源于黄河水利委员会《黄河泥沙公报》。

模型中使用的风场数据为欧洲中期天气预报中心的ERA-interim再分析资料,包含了每6 h的气象、水文格点数据。为了覆盖模型的计算区域,选取范围为 37°~ 41°N,117°~ 123°E,空间分辨率为0.125°×0.125°的10 m高度U、V方向风速和气压数据。

3.4.3 重要参数设置

模型中的底床摩擦力根据曼宁系数来确定,而曼宁系数M则由水深计算所得[21]:

式中,h为水深。当h≤1 m时,M统一为40。

3.4.4 模型验证

模型验证采用的数据通过黄河口附近定点水文观测(站位如图1所示)获得。为了使外海边界传入的潮波和研究区域内的波浪场充分发展,使流场更加稳定,模型模拟开始的时间比实测数据观测时间提前 30 d。

首先对渤海潮波系统进行率定分析。对模型后30 d的模拟结果进行调和分析,得到研究海区的潮汐调和常数。图2为模拟得到的渤海M2、S2、K1和O1分潮的同潮图,各分潮振幅和相位线与Ji等[22]的研究吻合。神仙沟外和秦皇岛附近各有1个半日潮无潮点,附近为正规全日潮区;渤海海峡处有1个全日潮无潮点,附近为正规半日潮区。渤海主要的半日潮由M2分潮主导,而主要全日分潮由K1分潮主导。

图2 M2、S2、K1和O1分潮的同潮图Fig.2 Simulated co-tidal charts of M2,S2,K1,and O1constituents

利用2006年11月SA1、SA2站同步观测数据,对模型进行流速流向验证。图3显示,各测站表底层的流速、流向的计算值与实测值吻合较好。

图3 2006年11月黄河口近岸海域定点表层和底层流速、流向模拟值与实测值对比Fig.3 Comparison of measured flow velocities and directions with the computed results in the coastal area of the Huanghe River Estuary in November 2006

采用 SS(Skill Score)对模型精度进行统计分析,SS的计算公式为

式中,Xmod和Xobs分别为统计的模拟值和观测值;n为统计量个数。SS大于0.65表示精度非常好,SS为0.5~0.65表示精度很好,SS为0.2~0.5表示精度好,SS小于0.2则表示精度差[23]。

统计结果显示,SA1站表层流速和流向的SS值分别为0.822 1和0.704 9,底层流速和流向的SS值分别为 0.775 3 和 0.649 2;SA2站对应的 SS 值依次为0.768 8、0.546 6、0.604 7 和 0.478 6,表明该模型模拟黄河口邻近海域流场具有很好的精度。

4 结果与讨论

4.1 冲淤变化

基于2007年和2015年黄河三角洲滨海区水下地形实测数据点,运用反距离加权法(IDW)插值得出空间分辨率为0.001°×0.001°网格点的水深值,得到黄河口滨海区冲淤变化的空间分布(图4)。

图4 2007−2015 年黄河口滨海区冲淤变化空间分布Fig.4 Spatial distribution of erosion and deposition in the coastal area of the Huanghe River Estuary from 2007 to 2015

2007−2015年,从现行河口至老清水沟河口南部的滨海区存在数个淤积中心和侵蚀中心,分别为:现行河口外泥沙大量淤积,成为最大的淤积中心,最大淤积厚度达到9 m以上;老清8口门外显著侵蚀,最大侵蚀深度大于4 m;老清8河口与老清水沟河口之间为淤积区,最大淤积厚度超过2 m;老清水沟河口外向南呈侵蚀−淤积−侵蚀马鞍状,侵蚀深度均大于2 m;最南端的淤积中心位于37.5°N,119.25°E附近,最大淤积厚度也超过2 m。此外,老清水沟沙嘴以南近岸海域轻微侵蚀,莱州湾西北部浅水海域有轻微的淤高。

4.2 沉积物特征与输运

4.2.1 粒度特征

粒度参数是泥沙输运、沉积和再分配等过程的集中体现,包含了沉积动力条件和运移趋势等信息[24]。泥沙搬运和沉积不仅受泥沙本身性状特征的影响,还受水动力条件等控制,因而粒度参数和沉积结构可指示水动力环境[25]。

从2007年和2015年表层沉积物中值粒径和分选系数的空间分布来看(图5),黄河三角洲海域表层沉积物整体上是近岸中值粒径大,离岸中值粒径小,而且普遍中值粒径大的区域分选性好。其中,在老清水沟沙嘴南侧存在一个中值粒径大、分选性好的表层沉积物中心,并且2015年相较于2007年该中心的位置向西移动。而在莱州湾中部,沉积物粒径较小,分选性较差。

图5 2007年和2015年黄河三角洲海域表层沉积物中值粒径、分选系数空间分布Fig.5 Spatial distribution of median grain size and sorting coefficient of surface sediment in the Huanghe River Delta in 2007 and 2015

2007−2015年中值粒径和组分变化显示(图6),在老清水沟沙嘴和老清8沙嘴外表层沉积物的中值粒径减小最显著。老清水沟河口外的侵蚀区与南部淤积区的中值粒径和粒径大于0.016 mm组分质量百分比均减小,整体上粒径大于0.016 mm组分质量百分比减小的中心比中值粒径粗化中心更偏西南。沙嘴岬角间和老清水沟南部海域则出现沉积物粗化,粒径大于0.016 mm组分也增多。

图6 2007−2015 年黄河三角洲海域底质变化Fig.6 Changes of submarine sediment composition of Huanghe River Delta between 2007 and 2015

4.2.2 粒度趋势

粒度参数在沉积物搬运过程中常由于物理、化学、生物等作用而发生沿程变化,通过对沉积物粒度的空间分布信息进行分析和挖掘,可以反演出沉积物的净输运趋势。McLaren和Bowles[25]提出,由于选择性起动、搬运和堆积,沿沉积物搬运方向,粒径将变小,分选变好,偏态更负偏;或者粒径将变大,分选变好,偏态更正偏。

本文选取2007年和2015年研究区部分样品,采用三次样条插值法[26],获得空间分辨率为 2 km×2 km的粒度参数,通过粒度趋势分析得到表层沉积物净输运矢量(图7),其中特征距离选取4 km。实测点所在区域的潮流属往复流性质,泥沙输运位于同一沉积体系中,符合粒径趋势分析的条件。结果显示,粒径趋势矢量在老清水沟沙嘴前缘指向南北两侧,其中向北指向了老清8河口和老清水沟河口之间的淤积区,向南较为复杂,远端的侵蚀区的粒径趋势矢量直接指向其西南侧的淤积区域,而近岸侵蚀区域的粒径趋势矢量先指向西南方向,后转向东南方向,同样指向了南部的淤积区。2015年与2007年相比,粒径趋势矢量汇聚的区域略微向西南方向偏移。在研究区的西南区域整体表现出向东北方向的运移趋势。另外,图中还显示老清水沟的东南区域有向东南方向的净输运趋势,主要是由于在强北风的作用下,切变锋减弱,泥沙更多地向莱州湾中部输运。

图7 2007年(a) 和 2015年(b) 表层沉积物粒度趋势分析净输运矢量分布Fig.7 Distribution of surface sediment transport vectors in 2007 (a) and 2015 (b)

4.3 悬沙输运机制

受涨落潮流的作用,黄河入海泥沙可分别向南、北输运。向南的泥沙沿着清水沟叶瓣体输运,对研究区的冲淤变化影响巨大。对河口入海泥沙输运机制的探究有助于更好地了解泥沙向南输运的机理。

根据2009年7月两次现行河口A1、A2站同步水文泥沙观测资料(站位如图1所示),进行悬沙通量机制分析。第1次在调水调沙期间,观测期间的利津水文站平均径流量为 3 445 m3/s,平均含沙量为 9.59 kg/m3;第2次在调水调沙之后,利津水文站的平均径流量为430 m3/s,平均含沙量为 3.14 kg/m3。起到主要作用的几个输移项列于表1,其他输移项由于量值很小,未予列出。A2站的第1次观测由于部分数据缺失未能计算各悬沙输移项。

表1 各动力项潮周期输沙量和方向Table 1 Tidal-averaged suspended sediment flux and direction of each dynamic term

计算结果显示,A1站在低径流量时潮流输沙项的方向与涨潮流方向一致,量值是平流输沙项的2.5倍左右,而高径流量则增强了落潮流及落潮流阶段的含沙量,使得潮流输沙项转为落潮流方向。A1站的潮泵输沙项均指向东南方向,与涨潮流方向一致,尤其是在低径流量期间(A1_第2次),潮泵输沙项达到20 kg/(m·s)。而A2站的潮泵输沙方向和落潮流一致,不过量值很小。潮泵输沙的机理为涨落潮的流速差异导致的挟沙能力不同,泥沙的再悬浮与沉降及其滞后效应致使悬沙浓度的涨落潮不对称变化,并与流速时间过程线存在一定的相位差,潮周期内涨落潮的悬沙输移不能相互抵消,从而产生净输沙[27]。A1站于6 m水深处,A2站于11 m水深处,由此可见,涨落潮的不对称使得在口门落淤的入海泥沙再悬浮向南净输运,尤其是在非调水调沙期间,而远岸泥沙少量向北输运。另外,两站点的斯托克斯漂移输沙项均指向东南方向,量值较为接近,对欧拉余流输沙项起到一定的平衡作用。黄河口海域的泥沙输运主要为平流项和潮泵项,前者主要受余流控制,后者主要由涨落潮的不对称引起,而净输沙则是直接决定了区域的冲淤变化。因此,对余流和潮流不对称的研究对于解释研究区的泥沙输运和冲淤演变具有重要意义。

4.4 流场特征

为了揭示研究区余流和潮不对称情况,基于建立的Delft 3D模型,分别模拟单纯由天文潮驱动和不同风场条件(冬半年、夏半年及强风天气)下波流联合作用的流场。

黄河三角洲地区受东亚季风影响,夏季以偏南风为主,冬季以偏北风为主,且风力较强。根据统计的2007−2015年夏半年(4−9月)和冬半年(10月至翌年3月)地面以上10 m的风场状况(图8)可知,黄河口夏半年常风向为SSW和S,强风向为NE;冬半年常风向为NE、NW和SSW,强风向为NE和NW。

图8 2007−2015 年黄河口风场玫瑰图Fig.8 Wind speed rose chart of the Huanghe River Estuary from 2007 to 2015

4.4.1 潮周期内流场分布及变化

根据模拟的黄河口滨海区潮周期的流速流向情况(图9),从现行河口向南至莱州湾,近岸潮流呈往复流,尤其是5~15 m等深线之间的区域,且流向与等深线基本平行。而在莱州湾西北部、清水沟以南区域,尤其是5 m等深线以浅海域旋转流特征明显。

图9 黄河口近岸海域潮周期深度平均流速矢量Fig.9 Velocity vectors of depth-averaged currents in a tidal cycle in the coastal area of the Huanghe River Estuary

潮周期内涨急、转流和落急时的表层流场显示(图10),在突出的沙嘴外流速最大。其中,老清水沟沙嘴外的高流速区在涨潮过程中的中心位置大致在37.5°N,119.3°E(图 10b),而落潮时的高流速区比涨潮时的高流速区离岸更近。

图10b和图10c显示,转流时该区域存在低流速带,在近岸水域产生并逐渐向东南深水区域移动,主要是因为内外转流的不同步,相反流向水流的剪切作用形成潮流切变锋[28]。潮流切变锋在一个潮周期内存在约1/3时间,使得泥沙输移大部分被限制在切变锋向陆一侧,主要在大约10 m等深线以内。

图10 潮周期内不同时刻的流场平面图Fig.10 Flow fields at different moments in a tidal cycle

4.4.2 余流

余流是对海水流动进行潮平均,去除周期性流动后剩余的非周期性流动。模拟的余流场显示(图11),从现行河口至老清水沟河口,5~10 m等深线之间的潮余流方向大体指向南,近岸余流方向较为复杂,总体上有离岸的趋势,而莱州湾的大部分区域余流都很小。在老清水沟河口以南有顺时针余环流,以东有逆时针余环流。近岸两个余流流速高值区分别位于老清8河口外和老清水沟河口外大致5~10 m等深线之间,但两者之间的区域余流较小。

图11 黄河口近岸海域垂向平均余流分布Fig.11 Distribution of vertically-averaged residual currents in the coastal area of the Huanghe River Estuary

4.5 涨落潮流速不对称

涨落潮流速不对称对区域物质输运有直接的作用。通过纯天文潮作用和冬半年、夏半年及强风天气下波流联合作用的水动力模拟,获取模型稳定后29.5 d(1个朔望月)的流场数据,利用偏度计算公式(2)计算涨落潮流速不对称的大小,绘制了现行黄河口滨海区涨落潮流速不对称的空间分布(图12),并根据相同方法模拟并计算纯天文潮作用下2015年的涨落潮流速不对称的空间分布(图13),其中图上箭头表示涨落潮急流时的垂向平均流向,并以涨潮流流向为正。

图12 黄河口滨海区垂向平均涨落潮流速不对称分布Fig.12 Distribution of depth-areraged flood-ebb velocity asymmetry in the coastal area of the Huanghe River Estuary

图13 2015 年纯天文潮作用下黄河口滨海区垂向平均涨落潮流速不对称分布Fig.13 Distribution of depth-areraged flood-ebb velocity asymmetry in the coastal area of the Huanghe River Estuary only tidal action in 2015

结果显示,在纯天文潮作用下(图12a),黄河口滨海区存在显著的涨落潮流速不对称现象,现行河口外表现为涨潮优势流,并呈舌状向南部莱州湾方向伸展,而近岸区域和莱州湾海域则普遍为落潮主导;与单纯的潮汐作用下的流速不对称分布相比,冬半年涨潮优势流分布区显著扩大和增强(图12b),夏半年近岸的落潮主导有所增强(图12c)。强风作用对涨落潮流速不对称影响显著,其中强北风作用下(图12d),涨潮优势区向近岸扩展并向莱州湾西部海域延伸;而在强南风天气条件下(图12e),近岸的落潮优势显著增强,外侧涨潮优势减弱,沙嘴岬角间外侧区域甚至转为落潮优势。

潮流不对称和地貌演变彼此影响,随着冲淤变化,研究区的潮流不对称不断调整。对比2007年和2015年纯天文潮作用下黄河口滨海区涨落潮流速不对称分布,2015年涨潮优势区域的舌状部分在莱州湾海域进一步延伸,主要是由于该区域泥沙淤积,水深变浅,涨潮优势得到发展,而黄河入海口转向北,使得口门的泥沙堆积中心也向北移动,导致部分孤东海域由落潮优势转为涨潮优势。另外,由于老清水沟南岸的西侧区域潮滩淤高,使得其近岸水域涨潮优势显著增强。整体上而言,研究区域的涨落潮流速不对称的空间分布未发生根本性改变,现行河口外的涨潮优势流都呈舌状向莱州湾深入,而近岸区域和莱州湾海域均普遍为落潮主导,由此可见,潮流不对称对泥沙的输运沉积和冲淤演变而言是一个长期存在的机制。

4.6 涨落潮流速不对称对冲淤演变的影响

涨落潮不对称及其优势的转变往往导致泥沙的输运和沉积,10 m水深以浅海域涨落潮优势转换的区域成为泥沙显著淤积区。现行黄河口外为涨潮流主导的水域,涨潮优势使得泥沙沿着涨潮流方向向南输运。近岸潮流属往复流,主流流向与等深线基本平行,而潮流切变锋的存在使得泥沙输移大部分被限制在切变锋的向陆一侧。而当泥沙从涨潮优势区进入落潮优势区后,落潮优势阻碍了泥沙继续顺着涨潮流的方向输运,使得泥沙产生落淤。最南端莱州湾中部的淤积中心与涨潮优势向落潮优势的转变处相吻合。黄河口床沙质和冲泻质的临界粒径为0.016 mm[29],图6显示,侵蚀区域的中值粒径和粒径大于0.016 mm组分质量百分比均显著减小,但与老清8河口外两变化中心位置重叠不同的是,老清水沟河口外粒径大于0.016 mm组分质量百分比减小的中心比中值粒径减小的中心更偏西南,淤积区粒径小于等于0.016 mm组分质量百分比增大最多,说明向南输运的黄河入海细颗粒悬沙更多的在南部淤积区落淤,使得底质中细颗粒组分增多。涨落潮优势的空间分布同样会影响区域的沉积物再搬运。沉积物粒度趋势显示,优势流发生转换的区域分别存在来自西南方向和东北方向的沉积物净输入,其中老清水沟外侵蚀区的粒径趋势矢量指向西南方向的淤积区,主要是因为该区域流速大,为余流高值区,且长期受涨潮主导控制,沉积物再悬浮后向西南方向输运。可见,南部淤积区泥沙主要是黄河入海泥沙和老清水沟沙嘴前缘的侵蚀泥沙。另外,沙嘴岬角间外侧5~10 m等深线之间区域由北向南呈强侵蚀−弱侵蚀−弱淤积分布,是因为涨潮优势逐渐减弱,强南风天气条件下,弱淤积区域甚至转为落潮主导(图 12e)。

强北风作用增强涨潮优势区,是莱州湾泥沙淤积和西北部海域底质粒度粗化的主要动力机制。老清水沟沙嘴南部海域主要由落潮主导,其中,37.6°N,119.25°E附近区域落潮优势相对较弱,在强北风作用下,该区域由落潮主导转变为涨潮主导,涨潮流输入泥沙堆积,形成近东西向的淤积区,老清水沟近岸侵蚀区被侵蚀泥沙可能是该淤积区沉积物的主要来源。莱州湾西北部沉积物较粗,且大部分时间水动力比较弱,不易发生底部泥沙的再悬浮。不过在强北风天气条件下,涨潮流优势区域会深入该海域,有利于泥沙的输入,且强风天气下水动力较强,输入泥沙粒径相对较粗,而平时落潮主导时水动力相对较弱,粗颗粒泥沙不易被带走,细颗粒泥沙顺着落潮流方向输运。图6也显示,该海域沉积物的中值粒径、粒径大于0.016 mm组分质量百分比都增大。整体上该海域底质冲刷较少,粗颗粒泥沙在强北风天气条件下输入后不易被落潮流带走,从而使得该海域发生微淤。

沙嘴岬角间的局部落潮优势造成泥沙淤积。近岸区域落潮流速大于涨潮流速,为落潮优势区,并且老清水沟与老清8两沙嘴岬角间水动力较弱,余流很小,受沙嘴岬角地形影响,两沙嘴的近岸区依次出现顺时针环流,即所谓的岬角余流效应,泥沙易在岬间外落淤沉积,形成淤积中心。沉积物粒度趋势也显示,岬间的粒径趋势矢量由老清水沟沙嘴前缘指向沙嘴岬角间的淤积区。

5 结论

(1)2007−2015年,从现行入海口至老清水沟沙嘴南部的滨海区依次存在数个淤积和侵蚀中心,其中除了现行河口外的强淤积区,在老清8河口和老清水沟河口之间存在一个淤积区,而最南端的淤积中心位于莱州湾中部。

(2)黄河口滨海区存在显著的涨落潮流速不对称现象,突出的沙嘴外缘呈涨潮优势流,并呈舌状向南部莱州湾方向伸展,近岸以落潮优势为主。在强北风天气作用下,涨潮优势区显著增强扩展,舌状涨潮主导区伸向莱州湾西岸。

(3)黄河口滨海区冲淤变化很大程度上受涨落潮流速不对称空间分布的控制,老清水沟东南海域由涨潮主导转为落潮主导,使得大量向南输运的泥沙在莱州湾中部沉积;近岸的落潮优势促使老清水沟沙嘴侵蚀泥沙被搬运至老清8河口与老清水沟河口之间的淤积区;而强北风天气条件下涨潮优势的增强使得莱州湾近岸区域淤积,沉积物粗化。

猜你喜欢
输沙黄河口涨潮
跟您走向黄河口
雅鲁藏布江中游悬移质输沙率的周期性与不均匀性分析
黄河口的小康图景
悠悠黄河口,浓浓黄河情
涨潮的时候
美哉黄河口
“水城”被淹
黄河下游河道洪水期输沙规律研究
黄河下游高效输沙洪水调控指标研究
神奇的潮汐