海岸泥沙运动的突出特征是受破碎波浪影响。海岸破碎波浪使水底泥沙悬移或层移而形成高浓度含沙量，悬移的泥沙在破波带内水流（水底回流或沿岸流、裂流等海岸环流）作用下，在垂直海岸方向和沿岸方向输移，从而形成沙坝、沙嘴和连岛沙洲等复杂海岸地貌形态（邹志利，2009；Dronkers，2009）。研究海岸破碎波作用下泥沙悬移的机理是研究上述海岸泥沙运动和地貌形态的关键问题。但由于波浪破碎和泥沙运动的复杂性，目前对于破波带泥沙悬移机理的研究还存在较多不清楚之处。有关研究中表达水底悬沙浓度一般的方法还是基于水底剪应力（其正比于波浪水质点速度或水流速度的平方） （Nielsen，1986；Van Rijn，1993，2007），但也有研究结果表明，若采用波浪破碎生成的湍流来表达水底处悬沙浓度则更能反映破波带内泥沙悬移的特征 （Steetzel，1990； Mocke et al，1993）。此外，研究也发现破波带内悬沙的水平对流作用也可能使得悬沙浓度不只与当地流态有关（Yoon et al，2012），但关于此贡献的定量结果还没有得到。

本研究通过物理模型实验研究破波带内水底悬沙浓度形成机理，即对于水底泥沙冲淤强度有关因素（波浪水质点速度、湍流能量）的影响进行了定量研究，给出了时间平均水底悬沙浓度与波浪水质点动能相关或者与湍流能量相关的条件。分析针对规则波、波群和不规则波3种波浪形态，并对破波带内爬坡区、内破波区和沙坝区的3个分区进行了讨论。根据分析结果讨论了适合于时间平均水底悬沙浓度计算公式的适用性，包括将其扩展到时间变化水底悬沙浓度的可行性。

1 实验内容

实验在大连理工大学海岸及近海工程国家重点实验室的大尺度波浪水槽进行。水槽长60 m、宽4.0 m、高2.5 m，实验平底处水深为2.0 m。水槽一端配置推板式造波机，另一端为天然沙构成的海岸模型，见图1（a），水槽中间加入两水泥隔离墙将海岸模型限制在宽度为1.0 m的范围，以保证实验的垂向二维特征。海岸模型由中值粒径D50=0.25 mm的天然沙组成，模型坡脚距离造波板18.5 m，总长26 m，沙层厚度0.4 m，初始剖面为1∶10的平面斜坡。实际剖面取决于实验过程的塑造。取坐标原点为海岸剖面与静水线的交点（不同波况原点位置不同），x坐标轴离岸方向为正，z坐标轴垂直向上为正。

Python是具有解释性、编译性、互动性和面向对象的高层次脚本语言。Python的设计风格具有很强的可读性，具有更有特色语法结构。支持面向对象技术，支持广泛的应用程序开发。Django框架是一个Python定制框架，它的核心组件包括用于创建模型的对象关系映射，为用户设计的完美管理界面，优秀的URL设计，友好的模板语言和缓存系统。

  

图1 实验水槽及实验布置（a）、实验海岸剖面及测点位置（b）

 

黑色竖线为浪高仪，红色竖线为流速和悬沙浓度测量断面。造波板位于x=41.5 m

为了考虑波浪形态的影响，实验波浪包括规则波、波群和不规则波 (规则波可以看作调制系数为零的波群)。波群调制系数δ取为0.5，其波面升高由下式确定

 

式中a和T分别为波群中短波波幅和周期；TG为波群周期，TG=nT，n为波群中短波的个数，实验中取n=10。不规则波生成谱采用JONSWAP谱，谱峰因子γ=3.30。表1给出了实验的波浪要素。

 

表1 实验波浪要素

  

波浪类型 Hrms/m T/s TG/s δ规则波 0.34 3.3波群 0.34 3.3 33 0.5不规则波 0.2 2.2

为了测量波面升高，由距造波板5.5 m到波浪破碎点布置10个浪高仪，从破碎点到岸线布置12～13个（依赖于测量波况）与流速和悬沙浓度测量剖面同位置的浪高仪。前者位置见图1（b）中破波带外测点标记，后者位置见图1（b）中破波带内测点标记。波面升高采用天津水运工程科学研究院生产的电容式浪高仪测量。

在校学生很多时候存在对本专业具体情况和发展动态了解不透彻、不全面的情况，而医学专业教师对本专业有更深入的了解和认识，一方面可以结合自身经验给予学生有效渡过大学阶段的建议；另一方面可利用自身对制度的了解与认识，如学校关于攻读研究生的条件及相关规定、医院对研究生的培养或人员招聘条件、导师对学生的要求等，从专业选择和就业等角度给予学生更准确、有用的信息，引导学生尽早做好职业生涯规划。

图4为规则波、波群和不规则波情况的（纵轴）与时间平均单位质量波浪水质点动能（=横轴）的相关性结果。从图中实验测点分布可以看出二者相关性：若实验点集中在一条斜线附近，则表明二者线性相关。图中实验点分别针对爬坡区、内破波区和沙坝区给出。由图可见，对3个波况与的相关性整体上是存在的，但规则波情况的沙坝区的实验点离散性很大，几乎不存在图中斜线所表达的变化趋势。为了表达图中实验点与图中斜线的接近程度，图中给出了3个波况拟合结果的决定系数r2。结果表明，不规则波最好（r2=0.83），波群与其接近（r2=0.8），规则波最差（r2=0.62）。后者结果最差的原因是由规则波情况沙坝区实验点的离散性较大所致，这是因为规则波情况波浪破碎更加集中，沙坝区卷破波产生的湍动能更为强烈，其可以由水面直达水底，所以对水底泥沙悬移（的产生）具有不可忽略的影响，而横轴中的波浪水质点动能没有包含该部分能量，所以造成其与相关性结果变差，因此实验点表现为具有很大离散性。与此相对比，不规则波情况沙坝区实验点离散性变得比规则波情况的要小，这是因为不规则波情况仅是大波在沙坝区破碎，而较小波在内破波区破碎，波浪破碎并不集中在沙坝区，所以在沙坝区与仍然具有较好相关性，但由于含有高强度的大波破碎，该相关性没有内破波区和爬坡区的好。

海岸地形采用武汉大学研制的URT-III三维动态底床地形测量仪测量，测量方式为超声式。地形测量仪测量范围为从静水岸线到破波点的地形，测量间距为0.01 m。对于静水线以上以及破波点外地形，利用直尺进行人工测量，测量间距为0.1 m。

Introduction of mechanism, application and evaluation methods of defoaming agents 2 40

  

图2 多普勒流速仪（ADV）和光学浊度仪（OTS）

每一波况的实验过程是，首先将波浪作用于海岸模型4小时形成平衡沙坝剖面并测量平衡剖面地形，然后测量平衡剖面的波浪以及各测点流速和悬沙浓度。测量分组进行，每一波况共作用17-19组波浪（波浪作用组次数目依赖于各波况流速和悬沙浓度测点数量，见附表1-3）。每组波浪作用40 min，共包括730个（规则波和波群）或1 090个（不规则波）短波。各组的测量在波浪作用海岸10 min之后开始，以使得泥沙起动并达到平衡状态。各测点测量时间为11min，并将其作为分析结果的样本长度。

  

图3 时间平均悬沙浓度空间分布

图3给出规则波、波群和不规则波3种波况各测量断面上的悬沙浓度时均值）的测量结果。为了分析（z）沿空间分布规律，根据水底地形特征将破波带分为3个不同形态的区域，即爬坡区、内破波区和沙坝区。爬坡区的特征是受波浪上爬和下落运动的影响，存在波浪爬坡回落所形成的靠近岸线的泥沙淤积地形，悬沙沿水深分布呈现幂函数型；沙坝区是波浪开始破碎所产生的沙坝和凹槽所在区域，该区域存在波浪以卷破波形式破碎产生的较大湍流，其对水底冲击作用强烈，因而对泥沙悬移和运动影响较大，悬沙沿水深分布也呈现幂函数型；内破波区则是位于爬坡区和沙坝区之间的区域，该区域波浪破碎处于饱和状态（破碎波波高受水深控制），悬沙沿水深分布较爬坡区和沙坝区更为不均匀，呈现出指数型分布。由于上述3个区域波浪运动状态不同，因此悬沙的分布和产生机理也存在差异，所以在下面的讨论中将分别考虑这3个区域的情况。

2 时均水底悬沙浓度生成机理

（φ=qb/，ψ = τbw/[ρ（s-1）gD]为希尔兹数，s=ρs/ρ，ρs为泥沙颗粒密度，g是重力加速度，D为泥沙粒径，ψc为临界希尔兹数），且考虑ψ＞＞ψc的情况，可知有qb正比于，代入式（7）可知正比于 τbw，即依赖于水底剪应力或者波浪水质点动能。

图5给出规则波、波群和不规则波情况的（纵轴）与时间平均水底湍动能（横轴）的相关性结果。由图可见，3种波浪情况的爬坡区、内破波区和沙坝区的实验点都可以分布在一条斜线附近，并且比图4中实验点对斜线的接近程度要高得多，这可由拟合结果的决定系数r2值看出：三个波况r2的值都在0.9以上，图4中规则波情况沙坝区实验点的较大离散在这里是不存在的。

 

式中C和U分别为水深平均悬沙浓度和水平流速，h为静水水深；Eb=为泥沙扬动函数（也称为冲刷函数），w′和C′分别为湍流周期波动垂向速度和悬沙浓度。该表达式中仅是将湍流脉动部分通过时间平均光滑掉（雷诺平均），波浪的波动部分还包含在其中。Db为泥沙沉降函数，由于沉降是由于泥沙在水底落淤引起，所以有表达式Db=ωCa（ω为沉速）。将Db表达式代入方程（2）可得Ca表达式：

2016年3月，广东省纪委发现邓强严重违纪，于同年3月22日对其立案并采取“两规”调查措施。其间，邓强如实交代了受贿事实。同年6月28日，邓强的亲属代为退出赃款260万元。

 

上式的左端中，第一项表示泥沙扬动函数Eb对Ca的贡献，即水底处的水流扰动会导致水底悬沙浓度Ca的形成；第二项（括弧内两项）表示泥沙沿水平方向输移（相邻的两沿水深垂直断面输沙通量存在差别U∂C/∂x≠0），加上水体含沙能力的改变（水深平均含沙量随时间变化∂C/∂t≠0）对水底悬沙浓度Ca的贡献，即二者的不平衡会导致泥沙向水底的沉降或水底泥沙的悬移。式中Ca分成两部分为随时间波动部分。时间平均部分是通常人们在一般的波浪平均模型所研究的对象，而时间波动部分在目前研究中关心的还不多，但其对于研究水底泥沙悬移机理是重要的，是需要考虑的对象。本节给出产生机理的分析，而C产生机理的分析将在下一节给出。

气象科技服务是公共气象的一个部分，但并不是完全无偿的。气象科技服务利用一定的科技信息，向人们提供一些知识，有一定的公益性。应从气象科技服务的公益性入手，为民众更好地提供公共气象服务，从而促进农村经济的发展。

为了分析时均水底悬沙浓度，对方程（3）做时间平均，有=0，对波浪垂直入射海岸情况又有时间平均流速=0（因为沿水深断面流量为零），在忽略水平扩散项（由波浪部分的乘积项关性不大，所以该项的值相对较小）后可得方程：

 

上式表明，时间平均悬沙浓度仅依赖于时间平均扬动函数，即水底悬沙浓度是由泥沙扬动函数所产生。在目前的研究结果中，一般是由水底剪应力τb来表达的（Partheniades，1965）：

 

式中 M 为常数 （单位：kg/（m2·s））；τe为临界冲刷剪应力；τbw为τb的幅值，其大小正比于水底边界层上方波浪水质点速度幅值um的平方：

 

式中，fw 为波浪摩阻系数，ρ为流体密度。由于式（6） 可表达为依赖于单位质量波浪水质点动能这样，根据式（4），也是依赖于波浪水质点动能。又由于水底剪应力所表征的是水底边界层湍流引起的作用力，所以以上结果的物理意义是：水底边界层产生的湍动能引起了对水底泥沙的扬动作用。这一作用也可从以下的表达式看出：

该公式中Ca和Eb是包含时间平均值和波动值的总的水底悬沙浓度和总的扬动函数。本节将通过论证与时变湍动能之间相关性来说明上式成立的条件于相关，即式（9）对时间变化也成立），即式 （3） 中 U∂C/∂x+ ∂C/∂t项的贡献可以忽略的条件。

 

上式中，qb是推移质输沙率，a是推移质厚度。若将 qb用Meyer-Peter公式（Meyer-Peter，1948）表达，

 

影响图3中水体含沙量的关键因素之一是水底悬沙浓度Ca（悬沙参考浓度）。波浪垂向速度以及湍流引起的垂向扩散作用使水底泥沙可以向上扩散，这样水底悬沙浓度的大小将决定整个水深含沙量的大小。本研究的目的是通过实验中得到的水底悬沙浓度Ca、波浪水质点速度和湍流速度（其形成湍动能ke的影响）的测量结果，来研究Ca与后两者之间的依赖关系，从而得到Ca的产生机理。分析中，将Ca取为最靠近水底测点的悬沙浓度，该距离为0.5 cm，对应的波浪水质点速度和湍流速度也取为同一位置处的ADV的测量结果。

以上依赖于水底剪应力这一结果的不足之处是其仅考虑了水底边界层湍动能对水底泥沙的扬动作用，事实上，在破波带内，波浪破碎产生的湍动能也可能作用到水底而引起泥沙扬动，若将这一作用考虑到的表达式中，则可采用下式表达（Kobayashi et al，2001）：

(3)重介质选和浮选效果差，实际生产数据与理论数据偏差大。主要表现在介耗高达4.5 kg/t，远高于同行业水平；循环水浓度高，为保证矿井生产，不得不强制洗煤生产，导致重介质选和浮选效果恶化；重介质操作困难，不稳定，导致矸石中小于1.80 kg/L密度级达到4%左右，洗混煤中小于1.40 kg/L密度级在10%上下。

“毫无疑问，对焦在这里是个难点，况且我用的还是一款为胶片机设计的古董镜头。但好在最后拍摄效果不错，这让我不得不对这台相机刮目相看。”

 

式中，WB=eBDB+efDf，eB和ef为效率系数，DB为波浪破碎产生的波能耗散，Df是水底边界层产生的波能耗散，二者分别可采用下式计算（Battjes et al，1978）：

 

式中α为大小接近1的常数，T为波浪周期，E=ρgH2/8为波能，Qb为波浪破碎概率，H为波高，k为波数，γb=H/h为波浪破碎指标。式（9）不仅考虑了水底边界层湍流产生的泥沙扬动作用（由Df表达），而且也考虑了波浪破碎产生的湍流对水底泥沙扬动作用（由DB表达），所以其比式（5）更能体现破波带内水流对泥沙扬动作用的特征。该式也说明事实上正比于波能耗散，因为DB和Df分别代表了波能耗散的两个方面，即波浪破碎产生的波能耗散和水底摩擦产生的波能耗散。而这两个波能耗散都可通过靠近水底边界层的湍动能来代表，因这一湍动能既来自于水底边界层耗散的波能（相当于式（5）所代表的波能耗散），又来自于波浪破碎耗散的波能（而不是如式（5）所示的仅来自于水底波浪动能耗散）。已有的研究结果表明，波浪水底边界层的上边界可取为泥沙层移运动的上边界，该边界一般距离水底有几个毫米的距离。Zanden等（2017）根据大水槽实验流速和浓度断面扫描仪测量结果分辨出该边界在距水底约0.005 m（5毫米）位置处。本文实验也在该位置进行悬沙浓度和流速的测量，所以该位置悬沙浓度结果可以作为水底处悬沙浓度而该位置流速测量结果给出的湍动能则可以作为反映产生的湍动能，或者波能耗散量。该湍动能可以通过下式计算

 

式中u′、v′和w′分别为水平、沿岸方向和垂向湍流流速。应用以上和的测量结果可以研究与的相关性，这样可以验证公式（9）的适用性，同理也可研究与时间平均波浪水质点动能的相关性来验证公式（5）的适用性。这两个结果分别由图4和图5给出。

流速、悬沙浓度和波面升高采取同点同步测量（图2）以便进行同步时域分析。为了测量流速和悬沙浓度的垂向分布，从岸线到沙坝区布置了12～13个垂向测量断面（见图1（b）红色竖线）。每个断面布置垂向间隔为0.5～10 cm的3～8个测点（根据当地静水水深确定，水深小则布点密度大，见图3，流速和悬沙浓度测点具体位置见附录附表1-3）。每个测点同时布置流速仪（ADV）和浊度仪（OTS）同点同步测量流速和悬沙浓度。流速仪为Nortek生产的Vectrino三维多普勒流速仪（探头为侧向和竖向型），测量精度为测量值的±0.5%～±1 mm/s。浊度仪为日本东京测量株式会社生产的PCT-2503光学浊度仪，测量精度为±0.02 g/L。

影响Ca的有关因素可从以下水深平均和湍流周期平均的悬沙输移方程（可由沿水深积分的泥沙浓度输移方程得到）看出

  

图4 时间平均水底悬沙浓度和水底波浪水质点动能拟合直线

以上两图的结果表明，与的相关性要明显高于与的相关性，所以由此可以推断采用式（9）计算比由式（5）计算精度更高，更适合整个破波带内（包括爬坡区、内破波区和沙坝区）的计算，因为式（9）中计算的量Wb是直接与湍动能成正比的，即WB既包含了水底边界层产生的能量耗散（Df，其产生了边界层的湍动能），又包含了波浪破碎产生的能量耗散（DB，其产生了由水面向水底扩散的湍动能），而式（5）仅包含了水底边界层产生的湍动能的贡献，即仅与水底剪应力有关（水底剪应力是边界层湍流流动对泥沙作用力的度量，也是水底边界层消耗水流能量的度量，即与Df有关），而不包含波浪破碎产生湍流的贡献。但是式（5）的不适用性主要体现在沙坝区，而在爬坡区和内破波区还能给出具有一定精度的结果，这表明仅考虑水底边界层湍动能对的贡献也有一定的适用性。其原因是爬坡区和内破波区，波浪破碎强度减弱，特别是破碎波湍动能直达水底的穿透能力较弱（波浪破碎不存在卷破），波浪破碎产生湍动能对贡献所占比例减小，所以式（5）仅考虑水底边界层湍流的贡献也是具有一定近似程度的。

  

图5 时间平均水底悬沙浓度和水底湍动能拟合直线

Zanden等（2017）仅针对规则波情况实验结果研究了相关性。他们的结果表明在内破波区水底悬沙浓度具有很好相关性（实验点离散性较小，他们也给出了破波带外的实验结果，离散性也较小）。但在沙坝区的实验点远远偏离内破波区的实验点，沙坝区不能用内破波区的的拟合直线来表达对（或者）的相关性。本文的沙坝区结果虽然也表现出较大离散性，但没有Zanden等（2017）的结果离散性大（无论是但从趋势上来说，本文结果与该研究结果是一致的，只是不同的实验条件使所得结果有所差别。

3 时变水底悬沙浓度生成机理

水底悬沙浓度Ca时变部分可由式（3）来分析，但与上节时均水底悬沙浓度仅由决定（式 （4）） 不同，因为∂C/∂t≠0 和 U≠0，还与式（3） 右端的 U∂C/∂x+ ∂C/∂t项有关，即不仅与当地水流作用有关，而且还与临近区域向本地的泥沙输移U∂C/∂x有关，也与本地泥沙对该输移的响应∂C/∂x有关。但后两者有抵消作用，二者之和的贡献可能并不占主要部分，若是这样，则主要取决于扬动函数，即仍然可以像一样采用式（4）计算，即有

 

随着互联网的高速运转，移动通讯及平板电脑的普遍使用，手机图书馆引进，OPAC系统的便利，微信的转发，使大数据在医学院校得以实现。

这里研究对依赖性的方法是：分析Ca时间历程的峰值出现时刻与湍动能ke时间历程峰值出现时刻的对应性，若二者出现时刻一一对应，则说明是由（或直接激发产生，而与U∂C/∂x+∂C/∂t项关系不大。与上节分析产生机理方法不同的是，这里不再给出和如图4和图5的相关性曲线，而仅讨论二者时间历程峰值出现时刻的相关性。

图6-图8分别给出了规则波、波群和不规则波情况波面升高η、水底湍动能ke和水底悬沙浓度Ca的时间历程。为了讨论这些结果在破波带内不同位置的差别，图中结果分别针对爬坡区、内破波区和沙坝区（凹槽和坝峰）给出。因为波群和不规则波情况和不明显随短波周期波动，所以图中利用短波平均水底悬沙浓度〈Ca〉和短波平均湍动能〈ke〉时间历程来研究波群和不规则波情况和特征。由图中结果可以看出，对规则波情况，峰值的出现是由η波峰破碎所引起，与η波峰出现时刻对应性在爬坡区和沙坝坝峰最好，在内破波区和沙坝凹槽较差；峰值与峰值出现时刻对应关系与此类似。对波群情况，短波平均湍动能〈ke〉峰值在内破波区与长波波峰对应，在爬坡区和沙坝区与长波波谷对应；悬沙浓度〈Ca〉较大峰值与湍动能〈ke〉较大峰值出现时刻基本具有一一对应关系。对不规则波情况，〈ke〉峰值出现时刻与波面、悬沙浓度〈Ca〉峰值没有以上简单对应关系，但悬沙浓度Ca与湍动能ke能量谱谱峰位置基本是对应的，后一结果在图9给出，所以在统计意义上（总能量分布上）也可以看作是和具有对应关系。

为了定量说明上述峰值与峰值对应关系，表2给出Ca与ke的相关系数rkC。表中规则波情况相关系数是将各物理量时间历程滤掉湍流周期波动项然后做相关性分析得到的，波群和不规则波情况相关系数为将各物理量时间历程短波平均然后做相关性分析得到的。从表中可以看出，对规则波情况，内破波区和沙坝凹槽处Ca和ke相关性较弱，对应的rkC值均较小（rkC=0.104，0.392），而爬坡区和沙坝坝峰处Ca和ke相关性较强，对应的rkC量值均较大（rkC=0.451，0.599）。这与上面所述图6中时间历程二者的峰值对应关系相符。对波群情况，在破波带内各区域rkC均量值较大（rkC=0.739，0.429，0.483，0.808），表明波群情况 〈Ca〉 对 〈ke〉相关性在整个破波带内都很好，这与上面图7中〈Ca〉与〈ke〉峰值出现时刻相同是一致的。对不规则波情况，爬坡区和内破波区的rkC也都有类似于规则波情况的较大值（rkC=0.501，0.428），沙坝坝峰处rkC量值有所降低（rkC=0.356），但沙坝凹槽处rkC接近零。所以，不规则波情况〈Ca〉对〈ke〉相关性总体来看要比图8中时间历程所表达的要好，且表现出rkC随水深减小而增加的趋势。

  

图6 规则波情况波面升高η、湍动能ke和水底悬沙浓度Ca时间历程（黑色实线为滤波后值）

  

图7 波群情况波面η（虚点线为长波，虚线为包络线）、平均湍动能〈ke〉和平均悬沙浓度〈Ca〉时间历程

  

图8 不规则波情况波面η（虚点线为长波，虚线为包络线）、平均湍动能〈ke〉和平均悬沙浓度〈Ca〉时间历程

  

图9 不规则波情况湍动能ke和水底悬沙浓度Ca能量谱

 

表2 水底悬沙浓度Ca和湍动能ke相关系数rkC

  

注：规则波情况rkC为ke与Ca相关系数；波群和不规则波情况rkC为〈ke〉与〈Ca〉相关系数。

 

波浪类型 爬坡区 内破波区 沙坝凹槽 沙坝坝峰规则波 0.451 0.104 0.392 0.599波群 0.739 0.429 0.483 0.808不规则波 0.501 0.428 -0.087 0.356

以上结果总的来看仍然可以认为与有近似对应关系，在工程计算中可以采用式（12）来近似计算规则波、波群和不规则波情况的水底悬沙浓度Ca（包括时间平均值和时间变化部分）。

综上所述，妊娠糖尿病高危孕妇极容易产生巨大儿、胎膜早破、胎儿窘迫和早产等不良妊娠结局。为了防范不良妊娠结局产生，医务人员应加强对产妇各个妊娠阶段的综合管理，为其提供科学的饮食、运动干预指导，促进高危孕妇合理膳食，形成母婴安全的良好前提。

4 结论

由于悬沙垂向扩散作用，水底悬沙浓度Ca是决定水体含沙量的关键因素。本文通过大尺度水槽模型实验测量了水底悬沙浓度Ca和流速的空间和时间分布，以此为基础研究了Ca的产生机理。分析针对时间平均和时间变化水底悬沙浓度分别进行，对前者给出了其与波浪水质点动能（表征水底剪应力大小）、湍动能（表征破碎波能量耗散和边界层能量耗散）的相关性曲线，对后者分析了其峰值与湍动能峰值出现时刻对应性。以这两个结果为基础，讨论了合适计算Ca的有关公式。主要结论如下：

（1）给出了时间平均水底悬沙浓度与时间平均波浪水质点动能、时间平均湍动能的相关性结果。结果表明，与的相关性要大于与相关性，特别在沙坝区。这说明采用波浪耗散计算式（9））比采用水底剪应力计算（式（5））更加合理。这一结果的原因是沙坝区内波浪卷破产生的湍流可以作用于水底，其对水底悬沙浓度的贡献是主要的。

（2）讨论了时间变化水底悬沙浓度峰值与湍动能峰值出现时刻的对应性，发现规则波情况峰值的出现是由波峰破碎所引起，二者对应性在爬坡区、内破波区和沙坝区有所差别，峰值与峰值出现时刻对应关系与此类似。短波平均的波群的湍动能〈ke〉峰值或与长波波峰对应（内破波区），或与长波波谷对应（爬坡区和沙坝区），悬沙浓度〈Ca〉较大峰值与湍动能〈ke〉较大峰值出现时刻基本具有一一对应关系。不规则波的湍动能〈ke〉峰值出现时刻与波面、悬沙浓度〈Ca〉峰值没有以上简单对应关系，但湍动能ke和悬沙浓度Ca能量谱谱峰位置基本是对应的。鉴于以上结果，可以将与正比关系推广到a与之间的关系：即以下关系既可以用于之间，也可以用于与之间

2019年，《青春期健康》杂志社将继续与北京走进崇高研究院合作，围绕传播崇高理念，报道先进典型人物与事迹，促进崇高文化建设，开展丰富校园活动，旨在为青少年树立良好的世界观、价值观、人生观起到指导作用。

 

式中Eb可通过WB=eBDB+efDf来计算，见式（9）。

犍为县母猪存栏较少，自产苗猪难以满足繁殖的需求，为满足肉猪市场需求，每年需要从外地调入大量的苗猪。部分养猪户存在侥幸心理，在自产苗猪生产后未能及时注射猪瘟疫苗，增加了疫病发生的几率；部分养猪户从外地采购苗猪，这些苗猪未经过检验、隔离观察，就与原猪舍的生猪混群饲养，一旦疫病发生，将对养殖场造成较大的危害，严重影响生猪养殖事业的发展。另外，生猪病死后，养殖户未能采取无害化处理，只将尸体掩埋在土壤内，加剧疫病的传播速度。

（3）应用以上公式时要注意到该公式不同情况下的精度是不同的。对时间平均水底悬沙浓度和时间平均水底扬动函数该公式具有较高精度，但对时变部分对不同波浪形态具有不同精度：规则波和不规则波情况下精度不如波群情况下精度高。因为波群情况整个破波带内〈Ca〉与〈ke〉有明显对应关系。而规则波情况仅在爬坡区和沙坝坝峰与有明显对应关系（二者相关系数rkC分别为0.451和0.599），但内破波区和沙坝凹槽情况对应关系较弱（rkC分别为0.104和0.392）；不规则波情况仅在爬坡区和内破波区〈Ca〉与〈ke〉有明显对应关系（rkC分别为0.501和0.428），在沙坝区〈Ca〉与〈ke〉没有明显对应关系（凹槽和坝峰rkC分别为-0.087和0.356）。虽然不规则波情况〈Ca〉与〈ke〉峰值没有明显对应关系，但Ca和ke能量谱谱峰几乎一一对应。所以，利用上述公式计算〈Ca〉也可以达到一定精度。

附录：

 

附表1 规则波情况流速和悬沙浓度测点位置

  

测量断面位置x/m当地静水水深h/cm各测点距当地水底距离 z′/cm 0.2 2 0.5，1，2 0.5 5 0.5，2，3.5，5 1.2 10 0.5，3.5，6.5，10 2.1 20 0.5，3.5，6.5，10，15，20 3 28 0.5，5.5，10.5，16，22，28 3.9 30 0.5，5.5，10.5，17，23，30 4.8 35 0.5，5.5，10.5，16，22，28，35 5.7 26 0.5，5.5，10.5，15.5，20.5，26 6.6 42 0.5，5.5，10.5，18，26，34，42 7.05 34 0.5，5.5，10.5，18，26，34 7.5 30 0.5，5.5，10.5，16，23，30 8.4 35 0.5，5.5，10.5，15，20，25，30，35 9.3 60 10，20，30，40，50，60

 

附表2 波群情况流速和悬沙浓度测点位置

  

测量断面位置x/m当地静水水深h/cm各测点距当地水底距离 z′/cm 0.2 2 0.5，1，2 0.5 5 0.5，2，3.5，5 1.2 10 0.5，3.5，6.5，10 2.1 18 0.5，3.5，6.5，12，18 3 21 0.5，5.5，10.5，15.5，21 3.9 21 0.5，5.5，10.5，15.5，21 4.8 22 0.5，5.5，10.5，16，22 5.7 26 0.5，5.5，10.5，15.5，20.5，26 6.6 45 0.5，5.5，10.5，18，27，36，45 7.05 40 0.5，5.5，10.5，20，30，40 7.5 34 0.5，5.5，10.5，18，26，34 8.4 40 0.5，5.5，10.5，16，22，28，34，40 9.3 60 10，20，30，40，50，60

 

附表3 不规则波情况流速和悬沙浓度测点位置

  

测量断面位置x/m当地静水水深h/cm各测点距当地水底距离 z′/cm 0.2 2 0.5，1，2 0.5 5 0.5，2，3.5，5 1.2 10 0.5，3.5，6.5，10 2.4 18 0.5，3.5，6.5，12，18 3.7 18 0.5，3.5，6.5，12，18 4.9 35 0.5，5.5，10.5，18，26，35 6.1 43 0.5，5.5，10.5，21，32，43 7.4 50 0.5，5.5，10.5，20，30，40，50 8.6 52 0.5，5.5，10.5，20，30，41，52 9.2 40 0.5，5.5，10.5，20，30，40 9.8 32 0.5，5.5，10.5，17，24，32 11 50 0.5，5.5，10.5，18，26，34，42，50

参考文献

Battjes J A,Janssen J,1978.Energy loss and set-up due to breaking of random waves.Coastal Engineering,569-587.

Dronkers J,2009.Dynamics of coastal systems.World Scientific.

Kobayashi N,Johnson B D,2001.Sand suspension,storage,advection,and settling in surf and swash zones.Journal of Geophysical Research:Oceans,106(C5):9363-9376.

Meyer-Peter E,1948.Formulas for bed-load transport.Proc of Congress Iahr.

Mocke G P,Smith G G,1993.Wave breaker turbulence as a mechanism for sediment suspension.Coastal Engineering,2279-2292.

Nielsen P,1986.Suspended sediment concentrations under waves.Coastal Engineering,10(1):23-31.

Partheniades E A,1965.Erosion and Deposition of Cohesive Soils.World Journal of Biological Psychiatry the Official Journal of the World Federation of Societies of Biological Psychiatry,2(4):190-192.

Steetzel H J,1990.Cross-shore transport during storm surges.Coastal Engineering Proceedings,1(22).

Zanden J V D,Hurther D,Cáceres I,et al,2017.Suspended sediment transport around a large-scale laboratory breaker bar.Coastal Engineering,125:51-69.

Van Rijn L C,1993.Principles of sediment transport in rivers,estuaries and coastal seas.Amsterdam:Aqua publications.

Van Rijn L C,2007.Unified view of sediment transport by currents and waves,II:Suspended transport.Journal of Hydraulic Engineering,133(6):668-689.

Yoon H D,Cox D T,2012.Cross-shore variation of intermittent sediment suspension and turbulence induced by depth-limited wave breaking.Continental Shelf Research,47(10):93-106.

邹志利,2009.海岸动力学.北京:人民交通出版社.