大风影响下的莱州湾西岸单站底层泥沙输运特征分析

李海 ，万凯，宋新，宋德海 ，王楠，鲍献文

(1.中国海洋大学 物理海洋教育部重点实验室，山东 青岛 266100；2.青岛海洋科学与技术试点国家实验室 海洋动力过程与气候功能实验室，山东 青岛 266237；3.中国海洋大学 海洋与大气学院，山东 青岛 266100；4.国家海洋局北海海洋工程勘察研究院，山东 青岛 266061；5.潍坊市海洋发展研究院，山东 潍坊 261100)

1 引言

海洋泥沙在海洋动力学环境和生态环境中发挥着重要的作用，以悬浮泥沙为主的海域可通过输运、沉降和再悬浮过程对周围海底地形地貌和泥沙组成产生影响。另外，由于泥沙颗粒物对营养盐和污染物的吸附/解吸附作用，泥沙迁移亦会影响整个海域的生物化学循环过程。近岸泥沙输运的机制研究表明，泥沙输运与粒径大小、海流流速、海底切应力大小等因素相关，易受波浪、潮流和外部风场条件的影响，故不同海域的泥沙输运呈现不同的特点[1-2]。

在渤海近岸地区，悬浮泥沙受水动力影响很大，一般情况下受潮流作用控制，但泥沙长期输运与多种因素相关。有关于渤海泥沙输运的周期性变化和通量分解研究已有一定的成果，如庞重光和于炜[3]对2000-2004 年渤海表层悬沙的卫星反演数据进行经验正交分解，发现了悬沙浓度（Suspended Sediment Concentration,SSC）存在年际和半年变化的显著周期。王海龙等[4]基于ECOMSED 模型研究分析认为，潮流作用下渤海大部分海域表层SSC 较低，且呈现明显的半日潮周期特征，潮流决定了近岸泥沙的空间分布特征。刘波等[5]发现莱州湾龙口附近海域潮致余流对于悬沙输运的贡献较大，在时间上以12～16 h 的变化周期为主。

国内外的最新研究表明，风暴是渤海泥沙输运的主导因素之一，特别是在冬季，东亚季风较强，由强流作用导致的沉积物再悬浮和泥沙输运易受到强风浪的影响，强烈的垂直混合在弱层结的条件下使得沉积物再悬浮，SSC 升高[6]。同时，除了风的强弱，水深也是影响再悬浮和输运情况的重要因素。渤海海峡南部较浅区域的表层SSC 显著大于北部，在南北不对称的大风再悬浮情况和交换流作用下，流出的海峡沉积物通量大于流入通量[7]。在2017 年2 月秦皇岛海域的观测中记录了一次风速大于14 m/s、持续时间大于6 h 的风暴事件，沉积物再悬浮分为湍流控制的再悬浮过程和底层切应力控制的再悬浮过程，前者导致底层SSC 呈脉冲式增大，后者将其扩散到水体中上层[8]。

前人研究重点探讨了泥沙在渤海及渤、黄海间的再悬浮和输运[6-10]，而聚焦于莱州湾的研究较少，且主要为SSC 的季节性变化和年际变化，往往忽略了更短周期的影响因素。基于此，本文拟采用集合经验模态分解（Ensemble Empirical Mode Decomposition,EEMD）、希尔伯特-黄变换（Hilbert-Huang Transform,HHT）和小波分析法对莱州湾底层泥沙的小尺度特征进行分析，并探究大风过程的影响，从而进一步认识莱州湾海域泥沙输运量的变化情况，同时为研究区域海岸工程建设与防护、科学管理与利用等提供参考，为海洋生态环境保护提供理论依据。

2 研究区域和数据介绍

莱州湾是渤海三大海湾之一，位于山东半岛北部，总面积约6 967 km2，水深在10～20 m，以不规则半日潮为主。由于其海岸线较长，在东部形成若干个较小的港湾，有黄河、小清河等多条河流入海[11]（图1）。莱州湾营养盐和有机物含量较高，渔业资源物产丰富，同时河流携带泥沙入海使得莱州湾拥有较高的泥沙含量。近几十年来，由于海岸工程建设、海水养殖设施修建、沿岸挖沙以及河流中上游大中型水库修建等人为因素的影响，入海输沙量减少，海岸受到了不同程度的侵蚀[12]。

图1 研究区域水深、连续站和风场数据站位Fig.1 Bathymetry map of the study area and the location of the measurement site and wind data site

课题组于2018 年10 月21 日至11 月6 日期间在莱州湾西部的C1 站（37.48°N，119.10° E，水深约5.5 m）进行了为期半个月的底基三脚架观测（图1），获得了水位、潮流、波浪和悬浮泥沙浓度等数据。观测时两台ADCP（Nortek Aquadopp 2 MHz）探头盲区为0.3 m，垂直分层单元0.1 m，仰视ADCP 实际测量距海底2.1 m至海面的流速剖面，俯视ADCP 实际测量距离海底0.02 m 到0.62 m 的流速剖面；一台Nortek ADV 用于测量单点流速并用于湍动能计算；一台RBR TD-WAVE用于获取水位和波浪数据；在支架的不同位置固定两个Seapoint OBS 浊度计用于测量水体浊度（NTU），测量仪器的设置参数见表1。此外，本文中的风场数据采用了欧洲中期天气预报中心的ERA-Interim 数据，其每天提供4 个时刻（00：00，06：00，12：00，18：00）的风场数据，空间分辨率为0.125°×0.125°，取距离C1 站最近的位置37.50°N，119.13°E。

表1 测量仪器设置Table 1 Settings of observation instruments

3 研究方法和数据处理

3.1 单宽输沙率和净输沙量

同期在2018 年10 月21 日和10 月30 日采集了海底近OBS 处的水样，并进行了实验室抽滤、烘干、称重，获取了SSC 数据。通过与现场OBS 测量的46 组浊度数据进行线性拟合，得到了浊度转换为泥沙浓度的关系式为

两者拟合的决定系数R2=0.88，表明两者具有较好的线性拟合关系，将上式应用到连续观测的浊度数据中得到SSC。取两个深度（距海底0.28 m 和1.02 m）的SSC 平均值作为底层SSC，利用俯视ADCP 测量的底层有效流速剖面（距海底0.02～0.62 m）做垂线平均作为当前时刻底层水体的流速，底层单宽输沙率可表示为

式中，z为距海底高度；SSC(z,t)为t时刻距底z层的泥沙浓度；u(z,t)为t时刻距底z层的流速u、v；F(t)为t时刻的单宽输沙率，单位为kg/(m·s)；C(t)和U(t)分别为垂线平均的泥沙浓度和流速u、v；H为底层水体厚度。将单宽输沙率进行时间积分，得到站点在观测时间段的底层净输沙量：

式中，T为整个观测时间段。

3.2 Hilbert-Huang 变换和边际谱

Hilbert-Huang 变换是一种信号时频分析方法，将1998 年Huang 等[13]提出的经验模态分解（Empirical Mode Decomposition，EMD）方法和Hilbert 变换相结合。EEMD 是EMD 的改进方案，用于分析非线性、非平稳的时间序列。它把一个复杂的数据分解为有限个不同时间尺度的振荡分量，通过在原数据中加入适当大小的白噪声来模拟多次观测的情景，经多次计算后作集合平均[14]。具体过程如下：

（1）将正态分布的白噪声ω(t) 加到原始单宽输沙率信号F(t)中，则总信号为

（2）将加入白噪声的信号作为一个整体，然后进行EMD 分解，具体可以表示为

式中，imf(i)为单次分解得到的第i个分量，Res 为残差项；

（3）重复步骤（1）和（2），每次加入新的正态分布白噪声序列，重复次数为N次；

（4）在添加N次白噪声后，对每一个imf 分量作集成平均处理，使得白噪声相互抵消，得到最终的EEMD分解结果：

式中，IMF(i)和Res 分别为N次平均后的第i个分量和残差项。将EEMD 分解后得到的每个IMF 分量做Hilbert 变换，得到相应的Hilbert 谱H(ω,t)。对于离散信号来说，边际谱表示信号中某一频率在各个时刻的幅值之和，对H(ω,t)进行时间积分可得到Hilbert 边际谱：

式中，h(ω,t)代表边际谱，ω为瞬时频率，t为时间，T为整个观测时间段。

3.3 显著性检验

在引入白噪声对原始信号进行EEMD 分解后通常要对固有模态函数(Intrinsic Mode Function，IMF)分量进行显著性检验，目的是分析白噪声对分解结果的影响，也有助于对IMF 分量物理意义的判断。采用IMF能量谱密度周期的分布规律分析其显著性[15-16]，第k个IMF 分量的能量谱密度为

式中，m为IMF 序列的长度，IMFk(j)表示第k个分量的第j个数据。以蒙特卡洛法对白噪声进行实验得到与平均周期的关系如下：

白噪声能量谱分布的置信区间表达式为

式中，α代表显著性水平。

3.4 交叉小波分析法

交叉小波分析法是以小波变换为基础，用于处理变量之间周期和能量关系的方法，广泛应用于信号处理领域[17]。定义任意两个时间序列X(t)和Y(t)，则它们的交叉小波功率谱（XWT）为

式中，a为尺度伸缩参数，τ为时间平移参数，CX(a,τ)为序列X(t)的小波变换系数，为序列Y(t)小波变换系数的复共轭。

交叉小波谱（WTC）可用相干性大小来表征，反映两个时间序列在时频空间的相干程度，表达式为：

式中，S为平滑算子。

3.5 湍动能计算

湍动能（TKE）由距离海底0.34 m 的ADV 测量并计算所得，3 个方向的流速可表示为：

4 结果分析

4.1 流速和泥沙浓度变化

图2 为U分量、V分量的剖面和底层流速，变化范围为-0.5～0.5 m/s。总体上看，流速U分量大于V分量，观测期间最大水深在6 m 左右，平均水深为5.55 m，水位和流速变化呈现了不规则半日潮的特征。在每个潮周期内，涨落潮流呈现明显的往复流特征，主轴方向为东西向。图2c 和图2d 显示了ADCP 及ADV所测底层流速的对比，测量深度分别为距海底0.32 m和0.34 m，ADV 所测流速由于进行了脉冲平均而更加平滑，两者的结果具有一致性。

图2 ADCP 观测的流速剖面U 分量（a）和V 分量（b），ADV 观测的近底层流速U 分量（c）和V 分量（d）及临近层位的ADCP 观测数据对比Fig.2 The ADCP measured U-component (a) and V-component (b) of the velocity profile,and the ADV measured U-component (c) and V-component (d) of the near-bottom velocity and corresponding ADCP measurement

观测期间的风矢量如图3a 所示，总体上东西分量大于南北分量，而在其他时间段以偏南风为主。由于潮流流速相对较大，为探究风对于海流的作用，特别是表层流对于风的响应，将余流做33 h 低通滤波以消除潮流的影响。图3b 和图3c 中不同的高值区体现了风场对于海流的作用强度和范围，流速大小随着风速的增加而增加，流向取决于风向。高值区大多位于水体上层，影响深度大约为距海底5 m 左右，从时间上看，V分量的高值区更多，受风影响的时间大于U分量；从强度上看，在风的驱动下，两者的最大流速相近，U、V分量分别为0.15 m/s 和0.11 m/s。图3d显示了观测期间两层SSC 的变化，下层略大于上层，弱风时小于0.02 kg/m3，在10 月26-30 日，11 月4-6 日出现较强西北风和东南风的情况下，SSC 间歇式增长，最大达0.25 kg/m3。同时，上下层SSC 在峰值处的差值随之增大，表明在外界风场的作用下出现了沉积物再悬浮现象。

图3 风矢量（a），亚潮频率的流速剖面U 分量（b）、V 分量（c）和悬沙浓度（d）随时间的变化Fig.3 Temporal variation of wind vector (a),U-component (b) and V-component (c) of the sub-tidal velocity profile,and suspended sediment concentration (d)

4.2 单宽输沙率多尺度分析

图4a 为东西方向和南北方向单宽输沙率的时间序列，正值分别代表东向输沙和北向输沙，负值分别代表西向输沙和南向输沙，表征泥沙输运量的瞬时变化。如图所示，平静天气下具有周期性的特征，主要受潮流影响，东西分量略大于南北分量，输沙率均小于5×10-3kg/(m·s)，10 月26-30 日以及11 月5 日后受大风影响，东西分量输沙率出现间歇式增长，其中西向输沙率的最大值超过0.02 kg/(m·s)。南北分量输沙率变化较小，在原有的周期上产生一定的高频脉动。定义当前时刻的累积输沙率为之前所有时刻单宽输沙率之和，东西方向累积输沙率在10 月26-28 日显著增大，11 月5 日以后在波动中减小，南北方向输沙率则一直呈现缓慢增长的趋势，最终两个方向的累积输沙率分别为0.26 kg/(m·s)和-0.50 kg/(m·s)（图4b）。

图4 单宽输沙率（a）和累积输沙率（b）随时间的变化Fig.4 Temporal variation of the sediment transport rate (a) and cumulative sediment transport rate (b) within unit-width

为研究单宽输沙率的多尺度特征，将两个方向的时间序列进行EEMD 分解，白噪声标准差为0.3，集合次数为100，分别得到9 个IMF 分量和1 个趋势项,各分量的平均周期和方差贡献率如表2 所示。为了使得到的分量物理意义更加明确，基于潮周期和观测时间对它们进行尺度划分，将平均周期约为1～7 h 的IMF 分量（IMF1、IMF2、IMF3）合并，记为高频分量，代表风浪和湍流作用下输沙率的变化；将平均周期约为13 h 和25 h 的IMF 分量（IMF4、IMF5）合并，记为潮周期分量，表示主要由半日及全日潮流控制的输沙率变化，因为该地半日潮特征显著，其方差贡献率大于全日潮分量；将平均周期约为2～5 d 的IMF 分量（IMF6、IMF7）合并，记为低频分量，主要受到小清河等河流径流量变化和天气尺度下水位低频波动的影响；将平均周期约为8～16 d 的IMF 分量（IMF8、IMF9）合并，记为长周期分量，显示了由于海水密度差异引起的定常流动所导致的输沙率变化。Res 保持不变（图5）。

图5 单宽输沙率EEMD 分解重构结果Fig.5 The reconstructed intrinsic mode functions of the EEMD on unit-width sediment transport rate

表2 单宽输沙率EEMD 分解各IMF 分量平均周期和方差贡献率Table 2 The period and variance contribution rate of each IMF of the EEMD decomposition on unit-width sediment transport rate

上述重构分量的平均周期逐渐增大，平均幅值依次减小，在不同尺度的分量中仍可观察到大风对于单宽输沙率的影响。在10 月26-30 日，11 月3 日以后各分量的振幅明显增大，随着平均周期增加，引起的振幅依次减小，表明风输入的能量对于单宽输沙率不同频率的变化均有促进作用，效果上高频大于低频。另外，低频和长周期的南北分量中出现了明显的振幅不对称现象，这是亚潮尺度下的流速变化和沉积物再悬浮差异共同作用的结果。Res 显示在该时间段内输沙率的变化趋势表现为西向输沙和南向输沙，其中南北方向Res 的方差贡献率达到了23.11%，是原序列的重要组成部分。在4 个分量中，东西方向高频分量的方差贡献率最大，南北方向潮周期分量的方差贡献率最大，分别代表各自单宽输沙率的显著变化，低频和长周期分量方差贡献率均较小，周期尺度特征较弱。

分别取式（10）中α=0.01、0.05、0.2 作各分量的显著性检验（图6），横纵坐标分别为平均归一化能量的自然对数（lnE）和平均周期的自然对数（lnT），若分量在置信度曲线之上，则表明其通过了该置信水平下的显著性检验。东西方向高频分量、潮周期分量通过了99%的显著性检验，高频分量能量较高，与方差贡献率对应，并印证了其显著性，低频分量和长周期分量也通过了80%显著性检验，二者所含物理意义相对略低；南北方向4 个分量均通过了99%的显著性检验，在对原信号的分解过程中，白噪声经过多次平均后，其影响较小。

图6 东西方向（a）、南北方向（b）各分量显著性检验Fig.6 Significance test of each component in east-west direction (a) and north-south direction (b)

4.3 大风影响下的单宽输沙率时频变化

将两个方向原始的9 个IMF 分量分别进行Hilbert 变换，将实数信号变换成解析信号，把一维信号变成二维复平面上的信号，得到Hilbert 时频谱（图7b，图7c），用于分析输沙率在不同时间尺度上所具有的能量，图中颜色代表能量，颜色越深表示信号在该时刻该周期的能量越大。在观测期间，最大风速为12 m/s，平均风速为5.5 m/s，将风速大于8 m/s 称作大风过程（图7a）。图7d 和图7e 分别为有效波高和距海底0.34 m处的湍动能变化，有效波高最大达到0.73 m；湍动能在风的影响下超过20 cm2/s2。从图7 中可以看到，单宽输沙率的时频变化与有效波高、湍动能变化并不是一一对应的。相对来说，湍动能的相位变化与单宽输沙率较为相符。

图7 风速（a），东西方向（b）和南北方向（c）单宽输沙率的Hilbert 谱，有效波高（d）和湍动能（e）随时间的变化Fig.7 Temporal variation of the wind speed (a),Hilbert spectrum of east-west (b) and north-south (c) unit-width sediment transport rate,significant wave height (d),and turbulent kinetic energy (e)

着重研究3 个大风过程，每个过程分为增长期（图7 中红色矩形方框）和衰减期（图7 中绿色矩形方框）。增长期的最大风速均大于8 m/s：10 月24 日06:00 至26 日06:00，风速从2 m/s 增加至12 m/s；10 月28 日00:00 至06:00，风速从5.8 m/s 增加至9.6 m/s；11 月4 日00:00 到18:00，风速从3 m/s 增加至9 m/s。在增长期内，有效波高变化不明显，第一个增长期内湍动能有略微的上升趋势，其他两个增长期内无明显变化，而单宽输沙率时频分布几乎保持不变。在衰减期内，风速缓慢减小并趋于稳定，波高和湍动能达到极值，单宽输沙率0～0.5 h-1频段的能量增强，并且出现相对高频的变化，周期为1 h 左右。第一个大风过程尽管达到了观测期间的最大风速，但由于跨度最大、风速梯度最小，只有波高变化较为明显，而湍动能规律性较差，输沙率也仅仅表现为潮周期和高频能量加强。也可能是由于前两个大风过程时间间隔很短，第一次大风过程的影响会部分体现在第二次中。

为反映单宽输沙率的实际频率成分，探究其在大风影响下的显著周期。依据式（7）得到两个方向的单宽输沙率的边际谱如图8 所示，可以看到幅值较大的部分都集中在相对低频的频段，其中东西方向最大幅值的频率为0.075 h-1，即观测期间单宽输沙率变化的显著频率，对应的显著周期约为13.3 h，与半日潮周期相符。南北方向高幅值频段为0～0.09 h-1，体现了11 h以上的显著变化周期，包含了更多亚潮尺度的变化，流速U分量的潮信号强于V分量。

图8 东西方向（a）和南北方向（b）单宽输沙率的边际谱Fig.8 The marginal spectrum of the east-west (a) and north-south (b) sediment transport within unit-width

4.4 风速与单宽输沙率的交叉小波分析

时频分析的结果显示，单宽输沙率在大风作用下会产生更强的周期性变化，且与风速的变化在时间上并不是同步的。为定量化两者的相位关系，接下来通过交叉小波分析研究风速和单宽输沙率在时频空间内的相干性。图9 为风场的东西和南北分量分别与该方向单宽输沙率的交叉小波功率谱和小波相干谱，图中黑色的锥形细实线为小波影响锥（COI）区域，内部黑色粗实线区域通过了95%的显著性检验，这里采用红噪声作为背景谱。交叉小波功率谱颜色代表时域内的频率强度，小波相干谱颜色为在该频率上二者之间的相干性强弱，箭头方向表示二者之间的相位关系，向右为正相位关系，向左为负相位关系。

图9 东西方向（a，b）和南北方向（c，d）单宽输沙率与风速的交叉小波功率谱（a，c）和小波相干谱（b，d）Fig.9 The cross-wavelet power spectrum (a,c) and wavelet coherence spectrum (b,d) between east-west (a,b) and north-south (c,d) unitwidth sediment transport rate and wind speed

交叉小波功率谱显示两者共同的高能量区集中在大于潮周期的低频部分，在该尺度下两者具有较强的共振周期。关注通过显著性校验的高值区域，主要位于风速较强的时间段，东西方向为28-30 日期间，南北方向在26-30 日期间，而11 月3 日之后的大风阶段由于边缘效应不予考虑。高值区内箭头指向表明东西方向上大风提前输沙率变化约1/4 个周期，南北方向上约为1/2 个周期。小波相干谱显示相干性较强的周期较为散乱，大多位于高能量区上，通过显著性检验的相干时频区域为低频区，即在潮周期尺度和亚潮尺度上两者相干性较强，箭头所示的输沙率滞后周期也与功率谱吻合。

5 讨论

5.1 渤海泥沙输运的小尺度特征

海洋环境中，泥沙的垂向运动与水平输运受潮流、海浪、环流及温盐结构等诸多因素影响[4]，不同周期尺度的热动力学条件综合决定了当地的泥沙运动情况。渤海悬沙浓度的潮周期尺度研究显示潮流输沙难以进行远距离输送，往往沉积在近岸[4]，故从图10 中可看到余流输沙远大于潮输沙，在所有分量中占主导地位，表明在比半月更长的时间尺度下，泥沙输运的方向和大小主要取决于渤海的余流特征。从EEMD 分解的结果出发，发现在半月时间尺度下，1～7 h 的高频振荡在单宽输沙率的变化中尤为突出，经过显著性校验后，可以排除白噪声的干扰，可能是在落潮期间垂向湍流混合不均匀，垂向速度梯度比涨潮期间大，垂向切变不稳定增强，泥沙浓度的高频扰动更加频繁和剧烈所致[18]。另外，注意到方差贡献率大小和边际谱显著周期不匹配（表2，图8），表明在进行时间积分后，落潮所导致的高频振荡信号和大风的能量输入影响减弱，因此输沙率呈现明显的潮周期尺度及亚潮尺度特征。

图10 各分量及余流分量的净输沙量Fig.10 Net sediment transport of each component and Res

5.2 大风过程对于渤海泥沙再悬浮及输运的影响

观测站点底部表层沉积物为粉砂，中值粒径为38.48 μm。通过LISST（C 型）激光粒度仪测得距海底0.9 m 处悬沙的中值粒径为38.65 μm，该地的沉积物粒径特征和悬沙粒径特征类似，水体中的泥沙基本为局地再悬浮泥沙。渤海悬沙浓度的季节性变化主要取决于风场等气候因素，且悬沙浓度与风速之间存在一定时间段的滞后相关[19]，而对于水深很浅的海域，风浪和湍流会加剧沉积物再悬浮的发生，经过垂向混合和扩散过程后悬沙浓度增大。图7 显示，湍动能在大风成长和衰弱过程中增强，而波浪则是在风速达到峰值后迅速攀升，进而共同增大了单宽输沙率并影响其时频分布特征，且输沙率相对于风速存在一定的滞后周期。根据崔廷伟等[20]基于遥感数据的研究，大风过程可能在几天内改变渤海表层悬沙浓度的分布结构，从HHT 时频分析的结果来看，风场对于输沙率的周期性振荡有显著的促进作用，若冬季发生范围广、强度大、持续时间长的大风事件，底层泥沙在强动力作用下大面积再悬浮至表层，渤海表层悬沙浓度的分布结构可能会发生变化。

5.3 泥沙输运主控机制

如图10 所示，底层净泥沙输运以东向和南向为主，总输沙量分别为305.77 kg/m 和597.25 kg/m。经过EEMD 分解得到的残差项Res 已经不具有明显的周期性特征，故可以将其看作余流分量，则东西方向输沙量从大到小依次为：余流分量、潮周期分量、低频分量、高频分量和长周期分量；南北方向输沙量从大到小依次为：余流分量、长周期分量、潮周期分量、低频分量和高频分量。在观测时间段内还是以余流输沙为主，低频分量和高频分量输沙量较小，另外，由于流速U分量的潮信号较强，其潮输沙量仅次于余流输沙量。尽管各潮周期分量在不同潮周期对泥沙输运的贡献值有差异（图11），但在平静天气下，涨落潮泥沙输运量大体是对称的，涨潮为西向和北向输运，落潮为东向和南向输运，导致整个全潮的净输沙量较低。在大风影响下，第14、15、32、33 和34 潮周期内输沙量较高，是由于风浪周期小于潮周期，风浪掀起的泥沙浓度高值存续期往往短于潮周期，如风浪掀沙出现在涨潮期间，则落潮期间风浪产生的再悬浮泥沙大部分已经沉降，这就会导致涨潮期间净输沙量增大，而落潮期间净输沙量基本不变；反之亦然。风浪掀沙的相位与涨落潮相位的不一致破坏了涨落潮期间的原本输沙量平衡，加剧了涨落潮输沙不对称，进而增加了潮周期内的净输沙量。

图11 潮流分量涨潮、落潮和全潮净输沙量Fig.11 The net sediment transport of tide current component during each flood tide,ebb tide and full tide

6 结论

本文根据2018 年10 月21 日至11 月6 日莱州湾连续观测站的泥沙数据，结合波浪、水位、流速，对底层单宽输沙率进行了EEMD 分解和HHT 变换，并探究了风对其影响，具体结果如下。

（1）两个方向EEMD 分解的重构结果分别得到9 个IMF 分量和1 个趋势项，根据其显著周期，可进一步分为高频分量、潮周期分量、低频分量和长周期分量，趋势项可看作余流分量。研究显示单宽输沙率的能量主要集中在高频和潮周期分量，表明了高频扰动和潮流对于泥沙输运速率的影响较大。

（2）以HHT 变换对单宽输沙率进行时频分析，风的影响主要发生在风速衰减期，风浪和湍流使得输沙率的0～0.5 h-1低频特征的显著性大大增强，并出现一定程度的弱高频变化。小波分析显示风速和单宽输沙率共同的高能量区及高相干区在低频频段，单宽输沙率的变化滞后约1/4 至1/2 个周期。

（3）边际谱通过对时间的积分弱化了高频的扰动和外界风场的影响，受流速分量潮信号的强弱影响，单宽输沙率的显著周期为半日潮周期和亚潮周期。

（4）底层总泥沙输运表现为东向和南向输沙，输运量分别为305.77 kg/m 和597.25 kg/m，浪潮频率的差异导致泥沙输运量在大风过程中出现显著的涨落潮不对称情况。