海洋学报

	振幅/m	频率/（（°）·h⁻¹）	相位差/（°）
仅M₂作用	工况1	M₂	2.0	28.985 5	0
工况2	M₂	2.5	28.985 5	0
工况3	M₂	3.0	28.985 5	0
M₂与S₂叠加作用	工况4	M₂	2	28.985 5	0
S₂	0.4	30	0
工况5	M₂	2	28.985 5	0
S₂	0.8	30	0
工况6	M₂	2	28.985 5	0
S₂	1	30	0

	振幅/m	频率/（（°）·h⁻¹）	相位差/（°）
仅M₂作用	工况1	M₂	2.0	28.985 5	0
工况2	M₂	2.5	28.985 5	0
工况3	M₂	3.0	28.985 5	0
M₂与S₂叠加作用	工况4	M₂	2	28.985 5	0
S₂	0.4	30	0
工况5	M₂	2	28.985 5	0
S₂	0.8	30	0
工况6	M₂	2	28.985 5	0
S₂	1	30	0

	黏土含沙量/ (kg·m⁻³)	粉砂含沙量/ (kg·m⁻³)	细砂含沙量/ (kg·m⁻³)
工况一	0.006	0.004	0.001
工况二	0.004	0.016	0.001
工况三	0.040	0.160	0.001

	黏土含沙量/ (kg·m⁻³)	粉砂含沙量/ (kg·m⁻³)	细砂含沙量/ (kg·m⁻³)
工况一	0.006	0.004	0.001
工况二	0.004	0.016	0.001
工况三	0.040	0.160	0.001

大小潮作用对潮滩沉积物层理影响的数值模拟研究

PDF下载

徐孟飘 ¹^,² , 东培华 ² , 马骏 ³ , 罗锋 ¹ , 张长宽 ¹ , 范代读 ⁴ , 周曾 ¹^,⁴^,⁵^,^*

海洋学报 | 论文 2021,43(10): 70-80

收起

海洋学报 | 论文 2021, 43(10): 70-80

大小潮作用对潮滩沉积物层理影响的数值模拟研究

全屏

徐孟飘¹^,², 东培华², 马骏³, 罗锋¹, 张长宽¹, 范代读⁴, 周曾¹^,⁴^,⁵^,^*

作者信息

¹河海大学水文水资源与水利工程科学国家重点实验室，江苏南京 210098

²华设设计集团股份有限公司，江苏南京 210014

³杭州市港航管理服务中心航道管理处，浙江杭州 310016

⁴同济大学海洋与地球科学学院，上海 200092

⁵河海大学江苏省海岸海洋资源开发与环境安全重点试验室，江苏南京 210098

徐孟飘（1993-），女，江苏省扬州市人，主要从事河口海岸水动力泥沙模拟研究。E-mail：mengpiao.xu@outlook.com

通讯作者:

周曾（1986-），教授，主要从事河口海岸地貌学研究。E-mail：zeng.zhou@hhu.edu.cn

The effects of spring-neap tide on sediment bedding on tidal flats: A numerical study

Mengpiao Xu¹^,², Peihua Dong², Jun Ma³, Feng Luo¹, Changkuan Zhang¹, Daidu Fan⁴, Zeng Zhou¹^,⁴^,⁵^,^*

Affiliations

¹State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering, Hohai University, Nanjing 210098, China

²China Design Group, Nanjing 210014, China

³Channel Management Office of Hangzhou Port and Waterway Administration, Hangzhou 310016, China

⁴School of Ocean and Earth Sciences, Tongji University, Shanghai 200092, China

⁵Jiangsu Key Laboratory of Coast Ocean Resources Development and Environment Security, Hohai University, Nanjing 210098, China

出版时间: 2021-10-25 doi: 10.12284/hyxb2021125

文章导航

摘要

收起

潮滩垂向沉积韵律层的形成主要取决于周期性的潮汐条件，包括涨落潮、大小潮、季节性及更长时间尺度的潮汐特征，为探究大小潮周期对潮滩沉积物垂向层理形成机制的影响，应用一维潮流泥沙与底床分层数学模型，对周期性潮汐条件作用下潮滩垂向沉积韵律层形成机制进行了数值模拟研究。结果表明，大小潮的周期性是模型中沉积层理表现韵律性的主要原因之一，韵律层中单个层理结构对应于1个大小潮周期过程，层理结构由形成于小潮期间的泥质层及形成于大潮期间的砂质层组成，层理的厚度也呈旋回性变化，大潮时层理较厚而小潮时层理较薄。水体边界含沙量是影响潮汐层理结构的重要因子，边界含沙量中粉砂占比增大会使潮汐韵律层整体粗化且砂质层厚度增大，当边界含沙量整体显著增大时，潮滩上的垂向潮汐韵律层会更加完整且厚度明显增大。潮汐层理的形成与特征是多种因子共同作用的结果，后续需进一步探究包括波浪、风暴潮、潮滩生物等其他因子的作用。

关键词

潮滩 / 泥沙分选 / 潮汐韵律层 / 数值模拟

Abstract

收起

The formation of vertical sedimentary rhythmic layers of tidal flat mainly depends on periodic tidal conditions, including flood and ebb tide, spring and neap tide, seasonal and longer term scale tidal characteristics. In order to investigate the distribution and mechanism of sediment bedding on tidal flats, a one-dimensional numerical model was used to simulate the rhythmic layers of long-term tidal flat bedding layers under spring and neap tidal cycles. Results indicate that the periodicity of spring-neap tide is the main reason for the rhythmicity of sedimentary bedding. One couplet in the rhythm layer corresponds to the spring-neap tidal period, which is formed by the mud-dominated layer during the neap tide and the sand-dominated layer during the spring tide. The thicknesses of layers also show a cyclical change: bedding layers are thicker during the spring tide and thinner during the neap tide. The boundary sediment concentration is also an important factor affecting the structure of tidal couplets. An increasing boundary concentration of silt makes the tidal rhythm layer coarser and increases the overall thickness of the sand-dominated layer. When the boundary sediment concentration significantly increases, the vertical tidal rhythm layers on the tidal flat are more intact with an evident increase in layer thickness. The formation and characteristics of tidal bedding layers are the result of the joint action of many factors (e.g., waves, storms, biological factors and etc.), which await further research effort in the future.

Key words

tidal flat / sediment sorting / tidal rhythm layers / numerical simulation

引用本文

徐孟飘, 东培华, 马骏, 罗锋, 张长宽, 范代读, 周曾. 大小潮作用对潮滩沉积物层理影响的数值模拟研究. 海洋学报, 2021 , 43 (10) : 70 -80 . DOI: 10.12284/hyxb2021125

Mengpiao Xu, Peihua Dong, Jun Ma, Feng Luo, Changkuan Zhang, Daidu Fan, Zeng Zhou. The effects of spring-neap tide on sediment bedding on tidal flats: A numerical study[J]. Haiyang Xuebao, 2021 , 43 (10) : 70 -80 . DOI: 10.12284/hyxb2021125

正文

收起

1 引言

收起

潮滩是细颗粒沉积物在潮流为主导的水动力条件下形成的宽阔而平坦的浅滩^[1]。潮滩上的水动力过程及沉积物输移过程呈双向反馈，集中表现为潮滩剖面形态的演变及潮滩上多组沉积物的分选现象。不同类型、属性的沉积物颗粒在潮滩上呈现不同的空间分布特征，该现象通常被称为沉积物分选。沉积物分选现象主要包括水平向分布和垂向分层两个方面，常见的水平向分布现象例如潮滩沉积物粒径的陆向细化，垂向分层现象一般通过潮滩底床沉积物层理来描述，沉积物在水平向和垂向的分选过程一般是紧密联系的。

层理是描述沉积特征常用的术语^[2]，潮汐沉积物中的潮汐层理展现出的周期韵律性可用于估算风暴潮事件、推算古地月轨道参数，有助于推算地史时期地−月系的演化历史^[3]，同时有助于分析潮滩沉积历史及提高对海平面变化规律的认知等^[3]，因此开展粉砂淤泥质潮滩上多组分泥沙的垂向分选机制的研究具有重要的实际意义。潮滩上潮汐的周期性作用（涨落潮周期、大小潮周期、季节性及更长时间尺度的周期）在沉积地貌上表现出的冲刷与淤积的交替促使了互层层理（砂质层和泥质层构成）的形成^[4]。对于互层层理的形成机制及内在机理的认知仍未形成统一认识，许多学者从现场实际观测入手分析其机制并取得一些进展。Deloffre等^[5]在对法国、英国等河口潮间带潮滩的研究中发现沉积物表现出大小潮周期性特征，王建等^[6]对江苏中部淤泥质潮滩进行实地观测，认为毫米级薄水平互层层理由半日潮产生，而厘米级厚砂泥互层层理则为半月天文潮（大小潮）的产物。范代读等^[7]和Fan等^[8]对长江口潮滩进行观测，发现潮汐层偶的厚度与大小潮周期潮差旋回性有关，层偶厚度是潮差的函数，随潮差由小潮−大潮−小潮而发生由薄−厚−薄的规律性变化。相较于现场实际观测方面，数值模拟层面的研究相对较少，是有待进一步深入的方向。本文以潮流为主导作用的潮滩为例，建立了潮滩一维概化动力地貌模型，分析沉积物垂向分布规律，探究周期性潮汐条件及不同边界含沙量对沉积韵律层结构（层理组成、层理厚度等）的影响。

2 研究方法

收起

本文使用Delft3D模型系统中的水动力模块、沉积物输移模块以及动力地貌演变模块来进行一维水动力及地貌演变的数值模拟计算。Delft3D模型系统是一套可进行水流、波浪、水质、生态、泥沙输移、床底地貌及各个过程之间相互作用数值模拟的强大开源软件包，在内河、河口、海岸、沉积物、地貌演变、水质等各方面被广泛应用。

2.1 模型介绍

2.1.1 水动力模块

假设潮滩沿岸均匀，可使用一维模式对潮滩剖面进行模拟，在静压假定和Boussinesq近似两个假定的基础上，对不可压缩流体的Navier-Stokes方程进行求解。一维浅水方程包括连续性方程和动量守恒方程。

（1）基本方程

① 连续性方程：

(1) $\frac{{\partial \eta }}{{\partial t}} + \frac{{\partial (hu)}}{{\partial x}} = 0,$

② 动量守恒方程：

(2) $\frac{{\partial u}}{{\partial t}} + u\frac{{\partial u}}{{\partial x}} = - g\frac{{\partial \eta }}{{\partial x}} + \nu \frac{{{\partial ^2}u}}{{\partial {x^2}}} - g\frac{{u\left| u \right|}}{{{C^2}h}},$

式中，x为空间水平坐标；u为垂线平均流速在x方向的分量（单位：m/s）；$h$为全水深，即水面到水底的距离（单位：m）；$\eta $为相对于某标准（例如：平均海平面，MSL）的水位值（单位：m）；v为平面紊动黏性系数（单位：m²/s）；g为重力加速度（单位：m/s²）；C为谢才系数（单位：m^1/2/s）。

模型中设置某一时间步长中当计算单元的水深小于0.1 m（临界值）时，处于“干”状态，该计算单元不参与计算。当该计算单元水深大于0.2 m（临界值）时，处于“湿”状态，该计算单元重新被激活参与模型计算。

（2）由潮流引起的底部切应力

计算公式为：${\tau _c} = \rho g{u^2}/{C^2}$，式中，$\rho $为水的密度（单位：kg/m³）。

（3）波浪模块

Delft3D默认的波浪模块为SWAN（Simulating Waves Nearshore）波浪模型，由于长周期模拟需要的计算时长较长因此本文未使用SWAN模块，而是采用了基于Young 和Verhagen^[9]以及van Rijn^[10]提出的公式开发的简单的风浪模型，该公式于2009年被证实适用于模拟浅滩上风浪产生的作用^[11]，该公式将无量纲波能量（$\varepsilon = {g^2}E/{U_{\rm{W}}}^4$）及无量纲谱峰波频率（$\varsigma = {U_{\rm{W}}}{T_{\rm{P}}}^{ - 1}/g$）与无量纲风区（$\chi = gF/{U_{\rm{W}}}^2$）及无量纲水深（$\delta = gh/{U_{\rm{W}}}^2$）联系起来，对已知风场（例如：风区长度、风速及水深）的波高及波周期进行估算。

(3) $\varepsilon = 3.64 \times {10^{ - 3}}{\left[ {\tanh {A_1}\tanh \left( {\frac{{{B_1}}}{{\tanh {A_1}}}} \right)} \right]^{1.74}},$

(4) $\varsigma = 0.133{\left[ {\tanh {A_2}\tanh \left( {\frac{{{B_2}}}{{\tanh {A_2}}}} \right)} \right]^{ - 0.37}},$

式中，${A_1} = 0.493{\delta ^{0.75}}$，${B_1} = 3.13 \times {10^{ - 3}}{\chi ^{0.57}}$，${A_2} = 0.331{\delta ^{1.01}}$，${B_2} = 5.215 \times {10^{ - 4}}{\chi ^{0.73}}$；${U_{\rm{W}}}$为海拔10 m处的风速；F为风区长度；${T_{\rm{P}}}$为峰波周期；E为波能，其公式为$E = \rho g{H_{\rm{S}}}^2/8$（${H_{\rm{S}}}$为有效波高）。当风条件给定时，解以上方程可估算有效波高${H_{\rm{S}}}$和峰波周期${T_{\rm{P}}}$。由于波浪在传播过程中波陡增大，波浪在浅滩上会破碎，衡量波浪破碎的典型破碎指标$\gamma $为有效波高和水深的比值（$\gamma = {H_{\rm{S}}}/h$），在本文中，参考Roberts等^[12]参数设置，设定$\gamma = 0.5$。在波浪破碎后，波高变为受水深限制，即${H_{\rm{S}}} = \gamma h$。当波浪破碎时，波浪引起的底部切应力达到最大值且随着水深的向陆减小而减小^[13]。根据线性波理论，底部最大轨道速度为

(5) ${u_{{\rm{w,b}}}} = \frac{{{\text π}{H_{\rm{S}}}}}{{{T_{\rm{p}}}\sinh (kh)}}.$

且由波浪引起的底床切应力为

(6) ${\tau _{{\rm{w,b}}}} = \frac{1}{2}\rho {f_{\rm{w}}}u_{{\rm{w,b}}}^2,$

(7) ${f_{\rm{w}}} = 1.39{\left[ {\frac{{{U_{{\rm{w,b}}}}{T_{\rm{p}}}}}{{2{\text π} (D_{50}^{{\rm{mix}}}/12)}}} \right]^{ - 0.52}},$

式中，$k$为波数；${f_{\rm{w}}}$为由Soulsby ^[14]定义的波浪摩擦因子；$D_{50}^{{\rm{mix}}}$为非黏性沉积物的中值粒径。综合底床最大切应力是引起沉积物输移的主要因素，由潮流和风浪引起的底床切应力的叠加通常大于两者的线性总和，在本模型中参考文献[14]，综合底床最大切应力的公式为

(8) ${\tau _{\max }} = {\tau _{\rm{w}}} + {\tau _{\rm{c}}}\left[ {1 + 1.2{{\left( {\frac{{{\tau _{\rm{w}}}}}{{{\tau _{\rm{c}}} + {\tau _{\rm{w}}}}}} \right)}^{3.2}}} \right].$

2.1.2 沉积物输运模块

模块中的沉积物分为黏性和非黏性两种，在计算沉积物的输运时，黏性部分与非黏性部分使用不同方法进行计算，即该模块忽略了沉积物混合物的临界再悬浮条件与单个种类沉积物的临界条件的差别。黏性沙和非黏性沙的物理过程（例如侵蚀、淤积和输移）均分为不同部分考虑。

黏性沉积物的输运使用带源汇项的平流方程来描述黏性沙的侵蚀与淤积过程

(9) $\frac{{\partial (ch)}}{{\partial t}} + \frac{{\partial (uch)}}{{\partial x}} = {Q_{{\rm{mud,e}}}} - {Q_{{\rm{mud,d}}}},$

式中，${Q_{{\rm{mud,e}}}}$与${Q_{{\rm{mud,d}}}}$分别代表侵蚀量与淤积量，该冲淤量使用被广泛运用的经典的Patheniades-Krone公式来计算，公式如下：

(10) ${Q_{{\rm{mud,e}}}} = \left\{ {\begin{aligned}& {{M_{\rm{e}}}\left(\frac{\tau }{{{\tau _{{\rm{cr,e}}}}}} - 1\right)\;\;\;\;\;\;{\rm{ }}\tau > {\tau _{{\rm{cr,e}}}}} \\ & {0{\rm{ }}\;\;\;\;\;\;{\rm{ }}\tau \leqslant {\tau _{{\rm{cr,e}}}}} \end{aligned}} \right.,$

(11) ${Q_{{\rm{mud,d}}}} = \left\{ {\begin{aligned}& {{\omega _{\rm{s}}}c\left(1 - \frac{\tau }{{{\tau _{{\rm{cr,d}}}}}}\right)\;\;\;\;\;\;{\rm{ }}\tau < {\tau _{{\rm{cr,d}}}}} \\ & {0{\rm{ }}\;\;\;\;\;\;{\rm{ }}\tau \geqslant {\tau _{{\rm{cr,d}}}}} \end{aligned}} \right.,$

式中，${M_{\rm{e}}}$为冲刷系数（单位：kg/（m²·s））；${\omega _{\rm{s}}}$为沉降速度（单位：m/s）；$c$是垂向平均浓度（单位：kg/m³）；${\tau _{{\rm{cr,e}}}}$和${\tau _{{\rm{cr,d}}}}$（单位：Pa）分别表示黏性沉积物的临界起动切应力及临界沉降切应力。该模型中使用Winterwerp^[15]提出的“连续沉积”的概念，因此黏性沉积物的临界沉降切应力${\tau _{{\rm{cr,d}}}}$取较大值（默认为1 000 Pa），则上公式可近似等于${Q_{{\rm{mud,d}}}} = {\omega _{\rm{s}}}c$。

对于非黏性沉积物，使用Sousby-van Rijn公式来计算潮流与波浪共同作用下的输移量。总沉积物输移量分为推移质和悬移质两类来计算

(12) ${Q_{{\rm{sand}}}} = \left( {{A_{{\rm{sb}}}} + {A_{{\rm{ss}}}}} \right)u{\left[ {{{\left( {{u^2} + \frac{{0.018}}{{{C_{\rm{D}}}}}U_{{\rm{rms}}}^2} \right)}^{1/2}} - {u_{{\rm{cr}}}}} \right]^{2.4}},$

其中，

(13) ${A_{{\rm{sb}}}} = \frac{{0.005{A_{{\rm{cal}}}}h{{\left( {D_{50}^s/h} \right)}^{1.2}}}}{{{{\left[ {\left( {s - 1} \right)gD_{50}^s} \right]}^{1.2}}}},$

(14) ${A_{{\rm{ss}}}} = \frac{{0.012{A_{{\rm{cal}}}}D_{50}^sD_*^{ - 0.6}}}{{{{\left[ {\left( {s - 1} \right)gD_{50}^s} \right]}^{1.2}}}},$

式中，${A_{{\rm{cal}}}}$是校准系数；${C_{\rm{D}}}$是仅有潮流产生的阻力系数；${U_{{\rm{rms}}}}$为波质点轨道流速的均方根（${U_{{\rm{rms}}}} = \sqrt 2 {U_{{\rm{w,b}}}}$）；${u_{{\rm{cr}}}}$为泥沙颗粒起动的临界流速；$s$为沉积物的相对密度；${D_*}$无量纲粒径大小（${D_*} = D_{50}^s{\left[ {g\left( {s - 1} \right)/{\nu ^2}} \right]^{1/3}}$）。

2.1.3 地貌演变模块

床面变形方程为

(15) $\frac{{\partial {z_{\rm{b}}}}}{{\partial t}} = \frac{{D - E}}{{{\rho _{\rm{s}}}}},$

式中，$D$是沉积物的沉积通量；$E$是沉积物的侵蚀通量；${\rho _{\rm{s}}}$为沉积物的密度。地貌演变的计算是基于黏性沙与非黏性沙这两种沉积物在每个水动力计算步长上的侵蚀或淤积量的更新。

2.1.4 底床分层模块

模型中将初始底床分为若干薄层，分为3大类：冲淤层、交换层以及基准层。冲淤层的厚度${\delta _{\rm{a}}}$是一个用户根据模拟需要选取的参数，在模型计算时冲淤层中的沉积物发生侵蚀或者淤积与水体中的沉积物进行交换，但厚度始终保持不变来保证模型的稳定性，交换层在计算过程中厚度可变化，但存在一个最大厚度${\delta _{\rm{u}}}$。

模型开始计算之前，不同种类沉积物（黏性沙及非黏性沙）由初始设定的比例进行均匀混合且被分为若干层，当模型开始计算时，最顶层的冲淤层会发生侵蚀或淤积，且冲淤层下的交换层数量及厚度也会随之发生变化。底床分层模型及模型中的侵蚀或淤积的计算过程见图1。当床面发生侵蚀时，冲淤层会发生冲刷且部分沉积物会被带走，当侵蚀量大于${\delta _{\rm{a}}}$厚的冲淤层时，其下的交换层会向冲淤层进行沉积物补充使侵蚀过程继续执行，在下一个计算步长开始之前交换层会对被侵蚀过的冲淤层执行补充过程，使冲淤层厚度再次回到定值${\delta _{\rm{a}}}$，此时交换层厚度减小，以此反映床面的侵蚀过程。当床面发生淤积时，冲淤层会变厚，在下一个计算步长开始之前会执行补充过程，冲淤层厚度将保持定值${\delta _{\rm{a}}}$将多余的沉积物向下传递并入交换层，若交换层达到所设定的厚度值${\delta _{\rm{u}}}$，则底床中将会产生一层新的交换层，以此反映床面的淤积过程。

2.2 模型建立

本文基于我国沿海典型粉砂淤泥质潮滩概化模型对潮滩沉积物垂向分层机制开展研究，研究使用的底床分层模型由于模型设置垂向层数较多，长期沉积过程计算时间很长，且本研究主要关注的是垂向的沉积物分布特征，因此选择一维模型更为合适，在不影响研究需求的前提下大大提高了模型计算效率。

粉砂淤泥质潮滩坡度较缓，约为1‰～3‰^[16]，本文基于Zhou等^[17]的研究进行模型设定，采用一维模型，初始床面为高程从–14 m至–1 m（平均海平面以下）的缓坡，且垂向分层层厚1 mm。模型使用50 m的等距矩形网格，共计280个网格，时间步长$\Delta t$=0.25 min，满足模型计算稳定性及精度要求。图2为初始地形剖面示意图，假设初始底床由50%的黏土、25%粉砂和25%的细砂组成，且这3种沉积物沿深度方向均匀混合。

模型的东侧为开边界其余均为固边界，开边界水动力条件主要使用M₂及S₂两个分潮，M₂为太阴主要半日分潮，S₂为太阳主要半日分潮。为研究垂向分层结构对不同潮流条件的响应，设置了6组工况进行计算，工况设置如表1所示。

由于粉砂淤泥质潮滩主要由黏土、粉砂及细砂构成，因此模型同时考虑黏性沙与非黏性沙，设定3种沉积物：黏土、粉砂和细砂。粉砂介于黏性和非黏性之间，处理时有时被当做黏性砂，但也有一些学者发现由于粉砂的成分与砂相似，透水性强、黏性差，一些国内外同行在研究中也有将之处理为非黏性砂。本研究中的粉砂粒径设置为30 μm，在粉砂界定的3.9～62.5 μm范围内，颗粒相对较大，将其考虑为非黏性沙。其中，黏土设置为黏性沙且以悬移质的方式输移，而粉砂与细砂均设置为非黏性沙，既可以悬移质的方式输移也以推移质的方式输运。初始底床由2 m的黏土、1 m的粉砂和1 m的细砂组成，3种沉积物均匀混合。经过调试，黏土沉速取为0.5 mm/s，临界起动切应力为0.2 N/m²，冲刷率系数取为5×10⁻⁴ kg/(m²·s)，为了使黏性沙始终处于沉降状态，因此黏性沙的沉降临界切应力取较大值设定为1 000 N/m²（默认值），黏性沉积物干密度为500 kg/m³。属于非黏性沙的细砂的中值粒径D₅₀取为90 μm，粉砂的中值粒径D₅₀取为30 μm，由于Delft3D软件内默认的非黏性沙计算公式不能设置粒径小于 63 μm的颗粒，因此本文调用Van Rijn^[10]公式进行非黏性沙的计算，非黏性沉积物干密度为1 600 kg/m³。在东侧开边界给予泥沙含沙量边界条件：黏土为0.006 kg/m³，粉砂为0.004 kg/m³，细砂为0.001 kg/m³。由于本研究是概化模型，侧重于揭示一些机制，若边界含沙量设置过大，易形成边界的剧烈淤积造成模型无法长期运行，因此选取的工况边界含沙量相对较小，但不会影响所得结论。

3 结果与分析

收起

3.1 垂向分层结构对潮流的响应

为研究潮流对潮滩底床沉积物垂向结构的影响，模型采用了单M₂分潮作用及M₂、S₂共同作用两种工况，其模拟的沉积物分选、分层特征结果如图3所示，计算结果是10 a接近平衡态。由于边界含沙量相对较低，潮间带发育相对较慢，且表现出不完全，但这对本文侧重说明的泥沙垂向层理现象影响较小。

在仅M₂分潮作用下，选取潮间带上的3点A、B、C点来进行对比分析（A点距离岸线0.5 km，B点距离岸线1 km，C点距离岸线1.5 km），图4为3个代表点仅在M₂分潮作用下垂向沉积物占比图。

由图4可见，潮间带上A、B、C 3点，由于底床高程的不断抬高，同一位置的水动力条件发生改变，各类沉积物垂向含沙量也发生改变。A、B、C 3点细砂垂向含沙量占比几乎为0。A点位于高程最大处，潮间带上部黏性沉积物占主导，因此垂向沉积物含沙量中黏土占比均大于粉砂占比。B点位于黏土与粉砂分布的过渡区域，垂向上在高程较小处粉砂占比较多，但随着沉积物淤积底床抬升，动力条件减弱，在高程较小处沉积物的垂向含沙量发生变化（黏土占比大于粉砂占比），而C点处于粉砂主导区域，因此粉砂占比整体大于黏土占比。尽管3个代表点在垂向占比上有所不同，但在垂向上含沙量曲线较为平滑没有明显波动，且均未见到明显的层理现象，可能原因包括：（1）模型中底床分层层厚设置较大（1 mm）限制了单个潮周期内潮汐层偶现象的形成；（2）模型中将粉砂当作一个组分，未考虑粗、细粉砂的不同行为过程；（3）模型未考虑絮凝沉积作用和细粒沉积物在滩上沉积后暴露空气发生脱水固结作用对再侵蚀的影响等。

图5为M₂和S₂分潮叠加的大小潮作用下潮间带上离岸0.5 km、1 km、1.5 km的A、B、C 3点处的垂向沉积物占比分布，由于潮间带上的沉积物主要由粉砂与黏土组成，因此图中只对比了粉砂与黏土的垂向占比分布情况。由图可见，相较于图4，垂向上黏土与粉砂的含量分布线呈锯齿状表现出明显的波动且具有一定的规律性，这种周期性的规律主要是由大小潮的周期旋回造成的。A点处高程较小处的锯齿状较高程大处更明显，是由于随着沉积物的不断淤积地貌不断变化使得水深在减小，潮动力相对的越来越弱，相对较粗颗粒的沉积物（粉砂）更难以被搬运至近岸处。从水平向沉积物的分布规律来看，B点位于黏土区域与粉砂区域的过渡区域，在此区域中黏土与粉砂对于潮汐能量的变化响应较好沉积物占比变化剧烈，因此锯齿状明显。而C点位于粉砂与细砂的过渡带，该点垂向上高程较小处细砂含量相对较大，因此由黏土与细砂形成的层理结构表现较差。

为研究沉积物垂向分布与潮汐周期的对应关系，选取离岸0.6 km、高程约为–1.4 m处的P1点以及离岸0.2 km、高程约为–1.1 m处的P2点进行分析（即模型的一个网格50 m），图6为黏土与粉砂的垂向分布情况，图中的沉积物于29 d内形成，包含两个大小潮周期，该段时间内的潮位过程、P2点累计层数及含沙量变化如图7所示。

图6a代表黏土的占比分布，颜色越红代表黏土占比越多，越蓝代表占比越少，图6b代表粉砂的占比分布，同理越红表示粉砂越多，越蓝代表粉砂越少，由图可见，P1点与P2点均包含完整的潮汐韵律层，韵律层由位于上部的砂质层及位于下部的泥质层构成，且韵律层的个数与大小潮周期数量相一致。P2点垂向层理结构中黏土占比均大于P1点，这是由于P2点较P1点高程更大，水动力条件相对弱，粗颗粒沉积物更难输移至较高处。由图7可见，该潮位过程呈现“大潮–小潮–大潮–小潮”的变化趋势，与潮位过程相对应，垂向分层结构中沉积物颗粒均呈现“较粗–较细–较粗–较细”的变化规律。同时潮汐大小潮周期性的变化规律也体现于层理厚度的旋回性变化之中，由于大潮期间潮汐能量大，沉积物在较强的动力条件下更易起动并发生输运，潮差大时层理较厚而潮差小时能量衰减，因此层理较薄。

3.2 边界含沙量对垂向分层结构的影响

现场观测表明，水体中悬浮沉积物的浓度也会影响大小潮周期层偶结构的组成。因此，为研究水体含沙量对于沉积层偶结构特征的影响，在相同水动力条件的基础上设置3种边界含沙量工况，如表2所示。3种工况中，工况三的含沙量量级与通常认为的相对接近，工况一及工况二的设置实际上是对层理结构中“粗–细”形成机制的对比试验及研究，工况二及工况三的设置实际上是对层理结构中“厚–薄”形成机制的对比试验及研究。此外，由于本研究是概化模型，侧重于揭示一些机制，在边界含沙量设置过大时，容易形成边界的剧烈淤积造成模型无法长期运行，因此选取的工况边界含沙量相对较小，但不会影响所得结论。

图8为不同工况下潮间带区域（0～4 km）黏土与粉砂占比分布对比，黏土占比图中颜色越红代表黏土占比越多，颜色越蓝代表黏土占比越少，粉砂占比图中颜色含义同理。在工况一情况下，黏土边界含沙量略大于粉砂含沙量，潮汐韵律层结构于离岸0～2 km处较为明显，主要是由黏土与粉砂组成，而2～4 km处主要是由粉砂与细砂组成的层理。在粉砂边界含沙量大于黏土边界含沙量的工况二的情况中，可见在0.5～4 km范围内粉砂占比明显增大，在离岸0～2 km处虽然仍是黏土与粉砂组成的层理结构，但不论是大潮期形成的砂质层还是小潮期形成的泥质层，层理中的沉积物颗粒一定程度上均发生粗化，而在离岸2～4 km范围内由粉砂与细砂构成的韵律层中的沉积物颗粒在一定程度上均发生了细化。工况三中由于黏土与细砂的含沙量是工况二中的10倍，由图8明显可见潮间带区域内整体淤积且最高处达到1 m左右，且离岸2～4 km范围内黏土与粉砂占比均有所增大，其中粉砂占比增大更明显，该范围内变为黏土、粉砂与细砂3种沉积物构成的层理结构，韵律层内的沉积物颗粒发生了明显的细化。并且，当边界含沙量浓度增大时，整个潮间带区域内大小潮对应的层理更加完整。

图9为代表点（离岸1 km处）不同工况下两个大小潮周期内黏土与粉砂的垂向占比分布，由图可见，当粉砂边界含沙量大于黏土边界含沙量时，潮汐韵律层中的沉积物均发生粗化。而当黏土与粉砂边界含沙量整体增大时（工况三）该代表点处两个潮周期内产生的潮汐韵律层明显增厚，但整体沉积物颗粒粒径变小，这主要是由于粉砂的边界含沙量虽然较大，但水体的挟沙能力有限，很大一部分粉砂较难随潮流输移至高滩处而是停留在高程较小处，而较大部分黏土却能输移至高滩处并沉积，因此粉砂占比反而减小。

4 结论

收起

潮滩垂向沉积韵律层的形成是多种动力因子共同作用的结果，其中大小潮周期性变化是韵律层形成的重要原因之一，本文着重从数值模拟角度探讨了大小潮周期的影响，分析了潮滩垂向沉积物的分选与分层结构形成的机理。具体结论如下：

（1）大小潮作用下，潮间带区域可见明显的分层结构（主要由黏土与粉砂形成），黏土与粉砂含沙量曲线沿垂向表现出明显的波动呈锯齿状。

（2）结果表明一个层理结构由泥质层与砂质层构成，其中泥质层在小潮期间形成，而砂质层于大潮期间形成，同时层理的厚度也与大小潮周期相关，大潮时层理较厚而小潮时层理较薄。

（3）除了周期性的潮汐条件这个主要因素，悬浮泥沙浓度也是影响大小潮周期性层理结构的重要因子。边界含沙量中粉砂占比增大会使得大小潮周期性层理结构中泥质层与砂质层均粗化且砂质层厚度增大，当边界含沙量整体显著增大时，潮滩上的垂向潮汐韵律层会更加完整且厚度明显增大。

在本文的研究中，在计算沉积物的输运时，黏性与非黏性沙使用不同方法进行计算，忽略了泥沙混合物由于黏性沙的存在呈现出的一定的黏性，可能会使得泥沙输运通量略有增大。同时，在潮滩底床垂向沉积物观测中，除了大小潮表现出的韵律性层理结构，涨落潮期间形成的潮汐层偶（基本沉积单元）以及波浪、风暴（风暴来临的水位、风暴历时）等动力因子作用下的层理重塑或破坏也常被观测到的现象，是未来数值模拟研究的重点，进而可以更加深入地掌握潮汐层理的形成与演化机制。

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

Friedrichs C T. Tidal flat morphodynamics: A synthesis[M]//Wolanski E, McLusky D. Treatise on Estuarine and Coastal Science. Waltham: Academic Press, 2011: 137−170.

Choi K S, Park Y A. Late pleistocene silty tidal rhythmites in the macrotidal flat between Youngjong and Yongyou islands, west coast of Korea[J]. Marine Geology, 2000, 167(3/4): 231−241.

Fan Daidu, Li Congxian, Archer A W, et al. Temporal distribution of diastems in deposits of an open-coast tidal flat with high suspended sediment concentrations[J]. Sedimentary Geology, 2002, 152(3/4): 173−181.

龚小辉, 柏春广, 王建. 淤泥质潮滩沉积周期性研究综述[J]. 南京师大学报(自然科学版), 2012, 35(1): 117−121.

Gong Xiaohui, Bai Chunguang, Wang Jian. Review of research on sedimentary periodicity of tidal mud flat[J]. Journal of Nanjing Normal University (Natural Science Edition), 2012, 35(1): 117−121.

Deloffre J, Verney R, Lafite R, et al. Sedimentation on intertidal mudflats in the lower part of macrotidal estuaries: Sedimentation rhythms and their preservation[J]. Marine Geology, 2007, 241(1/4): 19−32.

王建, 柏春广, 徐永辉. 江苏中部淤泥质潮滩潮汐层理成因机理和风暴沉积判别标志[J]. 沉积学报, 2006, 24(4): 562−569.

Wang Jian, Bai Chunguang, Xu Yonghui. Mechanism of silt-mud couplet of mud tidal flat and discrimination criteria of storm surge sedimentation in the middle Jiangsu Province[J]. Acta Sedimentologica Sinica, 2006, 24(4): 562−569.

范代读, 李从先, 邓兵, 等. 潮汐周期在潮坪沉积中的记录[J]. 同济大学学报, 2002, 30(3): 281−285.

Fan Daidu, Li Congxian, Deng Bing, et al. Tidal cycles recorded in tidal-flat deposits[J]. Journal of Tongji University, 2002, 30(3): 281−285.

Fan Daidu, Guo Yanxia, Wang Ping, et al. Cross-shore variations in morphodynamic processes of an open-coast mudflat in the Changjiang Delta, China: With an emphasis on storm impacts[J]. Continental Shelf Research, 2006, 26(4): 517−538.

Young I R, Verhagen L A. The growth of fetch limited waves in water of finite depth. Part 1. Total energy and peak frequency[J]. Coastal Engineering, 1996, 29(1/2): 47−78.

Van Rijn L C. Unified view of sediment transport by currents and waves. I: Initiation of motion, bed roughness, and bed-load transport[J]. Journal of Hydraulic Engineering, 2007, 133(6): 649−667.

Fagherazzi S, Wiberg P L. Importance of wind conditions, fetch, and water levels on wave-generated shear stresses in shallow intertidal basins[J]. Journal of Geophysical Research Earth Surface, 2009, 114(F3): F03022.

Roberts W, Le Hir P, Whitehouse R J S. Investigation using simple mathematical models of the effect of tidal currents and waves on the profile shape of intertidal mudflats[J]. Continental Shelf Research, 2000, 20(10/11): 1079−1097.

Green M O, Coco G. Review of wave-driven sediment resuspension and transport in estuaries[J]. Reviews of Geophysics, 2014, 52(1): 77−117.

Soulsby R. Dynamics of Marine Sands: A Manual for Practical Applications[M]. Thomas: Telford, 1997.

Winterwerp J C. On the sedimentation rate of cohesive sediment[J]. Proceedings in Marine Science, 2007, 8: 209−226.

龚政, 靳闯, 张长宽, 等. 江苏淤泥质潮滩剖面演变现场观测[J]. 水科学进展, 2014, 25(6): 880−887.

Gong Zheng, Jin Chuang, Zhang Changkuan, et al. Surface elevation variation of the Jiangsu mudflats: Field observation[J]. Advances in Water Science, 2014, 25(6): 880−887.

Zhou Zeng, Coco G, van der Wegen M, et al. Modeling sorting dynamics of cohesive and non-cohesive sediments on intertidal flats under the effect of tides and wind waves[J]. Continental Shelf Research, 2015, 104: 76−91.

2021年第43卷第10期

PDF下载

129

引用本文

BibTeX

文章信息

doi: 10.12284/hyxb2021125

接收时间：2020-04-08
首发时间：2026-02-26
出版时间：2021-10-25

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2020-04-08
修回日期：2020-06-30

基金

作者信息

¹河海大学水文水资源与水利工程科学国家重点实验室，江苏南京 210098

²华设设计集团股份有限公司，江苏南京 210014

³杭州市港航管理服务中心航道管理处，浙江杭州 310016

⁴同济大学海洋与地球科学学院，上海 200092

⁵河海大学江苏省海岸海洋资源开发与环境安全重点试验室，江苏南京 210098

通讯作者:

周曾（1986-），教授，主要从事河口海岸地貌学研究。E-mail：zeng.zhou@hhu.edu.cn

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/hyxb/CN/10.12284/hyxb2021125

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT