CN115346427A - 一种波浪作用下沙质海岸岸滩演变物理模型试验方法 - Google Patents

一种波浪作用下沙质海岸岸滩演变物理模型试验方法 Download PDF

Info

Publication number
CN115346427A
CN115346427A CN202211006619.4A CN202211006619A CN115346427A CN 115346427 A CN115346427 A CN 115346427A CN 202211006619 A CN202211006619 A CN 202211006619A CN 115346427 A CN115346427 A CN 115346427A
Authority
CN
China
Prior art keywords
wave
sand
model
test
water
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202211006619.4A
Other languages
English (en)
Other versions
CN115346427B (zh
Inventor
孙波
王宁舸
唐磊
厉泽逸
郑金海
戴鹏
陈大可
刘均卫
张弛
孙林云
肖凯
钱龙
郭天润
吴炳良
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Hydraulic Research Institute of National Energy Administration Ministry of Transport Ministry of Water Resources
Original Assignee
Nanjing Hydraulic Research Institute of National Energy Administration Ministry of Transport Ministry of Water Resources
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing Hydraulic Research Institute of National Energy Administration Ministry of Transport Ministry of Water Resources filed Critical Nanjing Hydraulic Research Institute of National Energy Administration Ministry of Transport Ministry of Water Resources
Priority to CN202211006619.4A priority Critical patent/CN115346427B/zh
Publication of CN115346427A publication Critical patent/CN115346427A/zh
Application granted granted Critical
Publication of CN115346427B publication Critical patent/CN115346427B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G09—EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
    • G09B—EDUCATIONAL OR DEMONSTRATION APPLIANCES; APPLIANCES FOR TEACHING, OR COMMUNICATING WITH, THE BLIND, DEAF OR MUTE; MODELS; PLANETARIA; GLOBES; MAPS; DIAGRAMS
    • G09B23/00—Models for scientific, medical, or mathematical purposes, e.g. full-sized devices for demonstration purposes
    • G09B23/40—Models for scientific, medical, or mathematical purposes, e.g. full-sized devices for demonstration purposes for geology

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Geology (AREA)
  • Mathematical Optimization (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Paleontology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Algebra (AREA)
  • Computational Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Mathematical Analysis (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Educational Administration (AREA)
  • Educational Technology (AREA)
  • Theoretical Computer Science (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明公开一种波浪作用下沙质海岸岸滩演变物理模型试验方法，确定模型的试验范围、定床边界位置和动床边界位置；选择模型沙并确定模型沙的中值粒径、模型沙的级配、模型沙的颗粒密度和模型沙的表观干密度；在试验水池中制作定床地形和动床地形；开展波浪率定工作：将试验水池加水到试验水位；待试验水池的水面稳定后，根据代表波要素设置波浪参数，开启造波机造波，调试造波机使试验水池中波浪要素满足试验要求；开展模型验证试验；根据模型验证试验中最后实测的模型沿岸输沙率，计算岸滩冲淤时间比尺；开展试验要求的岸滩演变试验组次。本发明方法能够更准确地模拟沙质海岸沿岸输沙过程，更合理地反映水沙动力与涉水建筑物的相互作用。

Description

一种波浪作用下沙质海岸岸滩演变物理模型试验方法
技术领域
本发明涉及一种波浪作用下沙质海岸岸滩演变物理模型试验方法，属于海岸工程技术领域。
背景技术
沙质海岸具有丰富的自然资源和优美的风光，是人类生产生活和休闲娱乐的主要聚居地之一。沙质海岸泥沙运动可以分为垂直海岸方向的横向(向离岸)运动和平行海岸的沿岸运动。前者控制岸滩的短期(小时-月)演变，后者控制岸滩的长期(年-数百年)演变。海平面上升、波浪作用、风暴、入海泥沙减小、人工采砂以及人类活动则是泥沙运动和岸滩演变的主要驱动因素或作用。沙质海岸岸滩演变，对于保护和合理开发沙质海岸资源、防灾减灾、规划和管理滨海生存空间具有重要影响，因此合理地预测沙质海岸岸滩演变一直是海岸工程界关心的核心问题之一。
由于近岸波浪浅化、破碎及其引起的波生流和强烈的水体紊动，沙质海岸近岸泥沙运动十分复杂，岸滩演变具有显著的多元性和非线性，准确预测沙质海岸岸滩演变存在极大的困难。目前，预测沙质海岸岸滩演变的方法主要包括经验法则和经验公式、数学模型试验以及物理模型试验。
经验法则往往仅能提供一些趋势性的判断，无法给出较为详细的演变过程和具体的泥沙冲淤量预测。经验公式往往针对某一具体的案例或某一特定的动力条件和海岸，缺乏通用性和可移植性。数学模型方便易用，近年来发展迅速，深受研究者的欢迎。但受制于目前对近岸泥沙运动机理认识有限，数学模型在反演破波带和冲流带泥沙运动过程等方面仍有较大欠缺，一定程度上影响数学模型预测海岸演变的准确性。此外，数学模型在精细化求解泥沙对流扩散过程时耗时较长，也限制了数学模型在模拟海岸长期演变方面的应用。整体物理模型试验能够比较直观且完整地重现原型海岸岸滩演变过程和变化，在岸滩演变研究中具有得天独厚的优势，因而得到广泛关注和应用。
受制于比尺效应、参数概化以及近海水动力和泥沙本身具有的复杂性，沙质海岸岸滩演变物理模型试验技术还远不成熟。已有沙质海岸岸滩演变物理模型试验在代表波确定、模型沙选取、沿岸输沙率计算、涉水建筑物处理、验证方法、加沙方法等诸多方面存在不足。这些不足影响物理模型模拟原型海岸水沙物理过程的准确性，进而影响模型预测岸滩演变和泥沙冲淤的可靠程度。
发明内容
本发明所要解决的技术问题是克服现有沙质海岸岸滩演变物理模型试验技术的缺陷，提供一种波浪作用下沙质海岸岸滩演变物理模型试验方法。
为达到上述目的，本发明提供一种波浪作用下沙质海岸岸滩演变物理模型试验方法，包括：
(1)确定模型的试验范围、定床边界位置和动床边界位置；
(2)确定原型海岸的代表波向、代表波高、代表波周期、泥沙中值粒径、粒径级配、泥沙颗粒密度、泥沙表观干密度以及总质量沿岸输沙率；
(3)确定模型的平面比尺、垂直比尺、波长比尺、波高比尺、波周期比尺、泥沙起动波高比尺和沉速比尺；
(4)根据泥沙的起动波高比尺和沉速比尺，选择模型沙并确定模型沙的中值粒径、模型沙的级配、模型沙的颗粒密度和模型沙的表观干密度；
(5)计算模型预估总质量沿岸输沙率；
(6)在试验水池中，制作定床地形和动床地形，并在动床地形和定床地形交接处设置过渡带；
若定床地形或动床地形上有涉水建筑物和人为地形，则在试验水池中相应的制作涉水建筑物地形和人为地形；
在动床地形的上游边界处设置加沙装置，在动床地形的下游边界处设置收沙装置，并按试验要求的波浪方向在外海边界处设置造波机；
在造波机两侧设置与造波机垂直的导波板；
(7)开展波浪率定工作：首先将试验水池加水到试验水位；
待试验水池的水面稳定后，根据代表波要素设置波浪参数，开启造波机造波，调试造波机使试验水池中波浪要素满足试验要求；
(8)开展模型验证试验；
(9)根据模型验证试验中最后实测的模型沿岸输沙率，计算岸滩冲淤时间比尺；
(10)开展试验要求的岸滩演变试验组次。
优先地，确定模型的试验范围，包括以下内容：
模型的试验范围包括岸线、岸滩和岸滩外海边界；
根据试验要求确定模型所需容纳的岸线和岸滩，岸滩外海边界应超出泥沙起动水深；
确定定床边界位置和动床边界位置，通过以下步骤实现：
动床范围以预设的研究区域为中心，动床边界设置在研究区域的岸滩演变范围之外，动床外海边界应满足水深大于泥沙起动水深的条件；
设定物理模型试验范围中除去动床范围的部分为定床范围，定床沿岸上下游边界与动床沿岸上下游边界间设置过渡段，过渡段长度大于5倍模型中试验波长。
优先地，确定原型海岸的代表波周期，通过以下步骤实现：
根据当地波高与波周期关系，确定代表波周期；
若无法获得当地波高与波周期关系，则按下式确定代表波周期：
Figure BDA0003809120830000031
式中：Tr为代表波周期；Ti为波高大于起动波高的第i级有效波高对应的波周期，Hi为波高大于起动波高的第i级有效波高，Pi为为波高大于起动波高的第i级有效波高对应的频率；
确定原型海岸的代表波向和原型海岸的代表波高，通过以下步骤实现：
若有原型海岸波浪测站的现场波浪观测资料，则代表波向和代表波高按下式确定：
Figure BDA0003809120830000032
Figure BDA0003809120830000033
式中：αr为代表波向，Hr为代表波高；Hi为波高大于起动波高的第i级有效波高，Pi为波高大于起动波高的第i级有效波高对应的频率，αi为波高大于起动波高的第i级有效波高对应的波向角；
若没有原型海岸波浪测站的现场波浪观测资料，则根据原型海岸波浪测站的深水波浪观测资料或后报的再分析资料，通过公式(2)和(3)分别计算代表波向和代表波高，并根据波浪浅水变形将深水代表波推算至原型海岸近岸浅水区，获得研究区域的代表波向和代表波高；
确定原型海岸总质量沿岸输沙率，通过以下步骤实现：
若有原型沿岸输沙率实测资料，则根据原型沿岸输沙率实测资料，确定原型海岸的总质量沿岸输沙率；
若无原型沿岸输沙率实测资料，则根据下式计算原型海岸的总质量沿岸输沙率：
Figure BDA0003809120830000034
式中，Qm为总质量沿岸输沙率；ρs为泥沙颗粒密度，ρ为水的密度；m为岸滩坡度；Ir为波浪破碎因子，
Figure BDA0003809120830000041
其中Hb为破波波高，L0为深水波长；ω为泥沙沉速；Eb为破碎处的波能，Eb＝0.125ρgHb 2，其中g为重力加速度；Cgb为破碎处的波群速度；αb为破波波向角；umb为破波底质点最大速度，
Figure BDA0003809120830000042
其中，hb为破波水深；Hb采用均方根波高。
优先地，确定平面比尺、垂直比尺、波长比尺、波高比尺和波周期比尺，通过以下步骤实现：
基于试验场地尺寸、研究区域范围和造波机造波能力，确定平面比尺λl和垂直比尺λh；当若试验要求的波浪绕射相似强于波浪折射相似，则平面比尺λl等于垂直比尺λh，波长比尺λL和波高比尺λH均等于垂直比尺λh；
按下式计算波周期比尺λT：
λT＝λh 1/2 (17)
泥沙的起动波高比尺等于模型的垂直比尺；
确定泥沙的沉速比尺，通过以下步骤实现：
基于悬沙对流输运相似和悬浮相似，确定泥沙的沉速比尺；
悬沙对流输运相似对应的沉速比尺为：
λω＝λh 3/2/λl (18)
式中：λω为沉速比尺；λh为垂直比尺；λl为平面比尺；
悬浮相似对应的沉速比尺为：
λω＝λh/λl 1/2 (19)
若模型为正态模型，则同时满足悬沙对流输运相似和悬浮相似；
若模型为变态物理模型且试验要求的悬浮相似强于淤积部位相似，则采用公式(5)确定泥沙的沉速比尺；
若模型为变态物理模型且要求淤积部位相似，则采用公式(6)确定泥沙的沉速比尺。
优先地，根据泥沙的起动波高比尺和沉速比尺，选择模型沙并确定模型沙的中值粒径、模型沙的级配、模型沙的颗粒密度和模型沙的表观干密度，通过以下步骤实现：
步骤4.1：选取刘家驹提出的波浪作用下的泥沙起动公式作为泥沙起动波高计算依据，得到模型沙的中值粒径、相对浮密度与垂直比尺的关系：
Figure BDA0003809120830000051
式中：Δ为相对浮密度，Δ＝(ρs/ρ)-1，其中ρs为泥沙颗粒密度，ρ为水的密度；λd为泥沙粒径比尺，d为中值粒径；λh为垂直比尺；
Figure BDA0003809120830000052
为在原型沙中Δgd+0.486×10-6/d的值和在模型沙中Δgd+0.486×10-6/d的值之比；
步骤4.2：选取武汉水利电力学院提出的泥沙沉速公式作为泥沙沉速计算依据，得到模型沙的中值粒径、相对浮密度与垂直比尺的关系：
Figure BDA0003809120830000053
式中：
Figure BDA0003809120830000054
为在原型海岸中
Figure BDA0003809120830000055
的值和模型海岸中
Figure BDA0003809120830000056
的值之比，其中υ为水的运动粘滞系数，λω为沉速比尺；
步骤4.3：若要求淤积部位相似，则分别根据公式(6)、(8)和(9)计算粒径比尺，确定模型沙的中值粒径和模型沙的颗粒密度；
若要求悬浮相似，则分别根据公式(7)、(8)和(9)计算粒径比尺，确定模型沙的中值粒径和模型沙的颗粒密度；
步骤4.4：根据下式确定模型沙的级配：
Figure BDA0003809120830000057
式中：di表示粒径累计频率等于i时对应的粒径；
Figure BDA0003809120830000058
为在原型沙di/d的值和模型沙中di/d的值之比；
模型沙的表观干密度由模型沙干燥后使用环刀称重法测得。
优先地，根据公式(4)，计算模型预估总质量沿岸输沙率；
在试验水池中，制作定床地形和动床地形，通过以下步骤实现：
根据原型海岸实测地形，按照平面比尺和垂直比尺制作定床地形和动床地形；
定床地形采用米沙、砖块或石子填充后进行水泥抹面，动床地形采用在以米沙、砖块或石子衬底上铺模型沙的形式制作，模型沙厚度超过原型海岸历史最大冲刷深度经垂直比尺换算后的模型冲刷深度；
在动床地形和定床地形交接处设置过渡带，通过以下步骤实现：
过渡带设置在动床地形和定床地形接触处，设置过渡带的方法为使用胶水将一层模型沙黏贴在定床地形表面；
若原型海岸上有涉水建筑物，则在试验水池中相应的制作涉水建筑物地形，通过以下步骤实现：
涉水建筑物根据垂直比尺和正态模型，采用不易变形材料制作涉水建筑物并进行水泥抹面，建筑物地形边坡坡度大于1:1的部分按1:1制作；
动床上游边界设置加沙装置，加沙装置包括PVC圆管、阀门、电磁流量计、水泵、水槽、搅拌器以及连接水管，水槽中存放模型沙与水混合形成的水沙混合液，PVC圆管、阀门、电磁流量计、水泵、水槽、搅拌器以及连接水管水平放置，PVC圆管、阀门、电磁流量计和水泵由连接水管按顺序连通，PVC圆管连通水槽，PVC圆管两端封闭，PVC圆管面向模型中心区域的一面均匀开设若干个圆孔，PVC圆管垂直输沙方向布置；
优先地，加沙装置的工作方法，通过以下步骤实现：
首先将模型沙和水以一定比例放置在水槽中；
通过搅拌器搅拌形成一定浓度的水沙混合液；
通过水泵将水槽中的水沙混合液抽至PVC圆管中；
水沙混合液通过PVC圆管上的圆孔流进试验水池的动床地形的上游；
加沙装置工作时，水槽中水沙混合液中悬沙浓度与水泵抽水流速的关系满足下式：
CSUadd＝Qadd (23)
式中，C为水槽中水沙混合液的悬沙浓度；S为连接水泵与PVC圆管的连接水管横断面过水面积；Uadd为连接水泵与PVC圆管的连接水管中水流平均流速；Qadd为试验加沙率；
动床下游边界设置收沙装置，收沙装置为一封底不封口且具有梯形横断面的槽体，收沙装置埋置在动床输沙方向的末端，收沙装置的封口与动床地形表面齐平。
优先地，开展模型验证试验，通过以下步骤实现：
步骤8.1：将试验水池加水到设定的试验水位；
步骤8.2：待试验水池的水面稳定后，设置造波机造出不使动床泥沙起动的波浪，造波机运行一段时间，直到模型沙密实后关闭造波机；
步骤8.3：待试验水池的水面稳定后，按代表波要素设置造波机参数，造波机开始造波；
步骤8.4：待造波机造出的波浪稳定后，加沙装置开始加沙，其中初次加沙率取值为模型预估总质量沿岸输沙率；
步骤8.5：运行时间tv后，测量模型沿岸输沙率和模型沿岸收沙率，关闭造波机和加沙装置；
待试验水池的水面稳定后，排干试验水池中的水，测量动床地形，获得试验结果；
试验结果包括模型沿岸输沙率和岸滩冲淤情况，岸滩冲淤情况由开展模型验证试验后的动床地形与开展模型验证试验前动床地形相减得到；
检验试验结果是否满足设定的验证要求，恢复初始的涉水建筑物地形、人为地形、动床地形；
步骤8.6：若步骤8.5中试验结果未达到验证要求，则通过调节加沙率和代表波要素，重复步骤8.1至步骤8.5，直至试验结果满足验证要求；
步骤8.3中按代表波要素设置造波机时，若采用不规则波，则采用适合当地波浪条件的波浪谱，并设置均方根波高等于代表波高，谱峰周期等于代表波周期；
步骤8.5中时间tv由公式tv＝tpv/λt确定，其中tpv为模型验证时间tv对应的原型岸滩演变时间；λt为岸滩冲淤时间比尺；
步骤8.5中检验试验结果是否满足验证要求，通过以下步骤实现：
(a)对于原型海岸地形稳定的岸滩，验证模型验证试验获得的模型沿岸输沙率是否达到平衡以及岸滩冲淤情况是否稳定，若模型验证试验获得的模型沿岸输沙率达到平衡且岸滩冲淤情况稳定，则认为试验结果满足验证要求；
(b)对于原型海岸地形不稳定的岸滩，则验证对应的实际岸滩演变过程，若模型验证试验中岸滩演变与实际岸滩演变过程一致，则认为试验结果满足验证要求。
优先地，开展试验要求的岸滩演变试验组次，通过以下步骤实现：
步骤10.1，将动床地形设置为试验要求的初始地形；
步骤10.2，将试验水池加水到试验要求的水位；
待试验水池的水面稳定后，设置造波机造出不使动床泥沙起动的波浪，造波机运行一段时间，直到模型沙密实后关闭造波机；
步骤10.3，待试验水池的水面稳定后，按试验波浪要求设置造波机参数，造波机开始造波；
待造波机造出的波浪稳定后，加沙装置开始加沙，加沙率取值为模型验证试验中最后实测的模型沿岸输沙率；
步骤10.4，运行时间tt后，关闭造波机和加沙装置，排干试验水池中的水，测量动床地形；
步骤10.5，如试验尚未完成，加水至试验要求水位，重复步骤10.3和步骤10.4直至该组试验完成；
步骤10.4中时间tt由试验要求的原型岸滩演变典型时间除以岸滩冲淤时间比尺得到。
优先地，按如下公式计算岸滩冲淤时间比尺：
Figure BDA0003809120830000081
式中，λt为岸滩冲淤时间比尺，
Figure BDA0003809120830000082
为即原型沙表观干密度与模型沙表观干密度之比；
Figure BDA0003809120830000083
为实际的总质量沿岸输沙率比尺。
本发明所达到的有益效果：
本发明改进了已有沙质海岸岸滩演变物理模型试验中有关代表波确定、模型沙选取、沿岸输沙率计算、涉水建筑物处理、验证方法、加沙装置及加沙方法方面的不足，本发明方法能够更准确地模拟沙质海岸沿岸输沙过程，更合理地反映水沙动力与涉水建筑物的相互作用。本发明方法不仅适用于沙质海岸岸滩在自然状态下的演变过程预测，也适用于人类活动影响下的岸滩演变预测；本发明方法不仅适用于平常浪动力，也适用于极端大浪动力；本发明方法不仅适用于规则波，也适用于不规则波；且具有准确性高、步骤简单、实施方便、节约成本等优点，具备广阔的应用前景。
本发明方法与已有沙质海岸岸滩演变物理模型试验方法相比，本发明提出的方法具有更好的适用性和更高的准确性。
附图说明
图1为本发明实施例二中物理模型试验范围、平面布置及初始地形示意图；
图2为本发明中加沙装置的示意图；
图3为本发明实施例二的原型海岸中港口、航道、防浪挡沙堤以及T型防护丁坝建设3年后岸滩地形图；
图4为本发明实施例二的原型海岸中港口、航道、防浪挡沙堤以及T型防护丁坝建设9年后岸滩地形图。
图中标识：1-拟建港口，2-拟建航道，3-拟建防浪挡沙堤，4-拟建T型防护丁坝，5-动床上边界，6-动床下边界，7-动床外海边界，8-加沙装置，9-收沙装置，10-造波机，11-导波板，12-PVC圆管，13-阀门，14-电磁流量计，15-水泵，16-水槽，17-连接水管，18-搅拌器。
具体实施方式
以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。
需要说明，若本发明实施例中有方向性指示(诸如上、下、左、右、前、后......)，则其仅用于解释在某一特定姿态下各部件之间的相对位置关系和运动情况等，如果该特定姿态发生改变时，则该方向性指示也相应地随之改变。
实施例一
本发明提供一种波浪作用下沙质海岸岸滩演变物理模型试验方法，包括：
波浪作用下沙质海岸岸滩演变物理模型，简称模型；
(1)确定模型的试验范围、定床边界位置和动床边界位置；
(2)确定原型海岸的代表波向、代表波高、代表波周期、泥沙中值粒径、粒径级配、泥沙颗粒密度、泥沙表观干密度以及总质量沿岸输沙率；
(3)确定模型的平面比尺、垂直比尺、波长比尺、波高比尺、波周期比尺、泥沙起动波高比尺和沉速比尺；
(4)根据泥沙的起动波高比尺和沉速比尺，选择模型沙并确定模型沙的中值粒径、模型沙的级配、模型沙的颗粒密度和模型沙的表观干密度；
(5)计算模型预估总质量沿岸输沙率；
(6)在试验水池中，制作定床地形和动床地形，并在动床地形和定床地形交接处设置过渡带；
若定床地形或动床地形上有涉水建筑物和人为地形，则在试验水池中相应的制作涉水建筑物地形和人为地形；
在动床地形的上游边界处设置加沙装置，在动床地形的下游边界处设置收沙装置，并按试验要求的波浪方向在外海边界处设置造波机；
在造波机两侧设置与造波机垂直的导波板；
(7)开展波浪率定工作：首先将试验水池加水到试验水位；
待试验水池的水面稳定后，根据代表波要素设置波浪参数，开启造波机造波，调试造波机使试验水池中波浪要素满足试验要求；
(8)开展模型验证试验；
(9)根据模型验证试验中最后实测的模型沿岸输沙率，计算岸滩冲淤时间比尺；
(10)开展试验要求的岸滩演变试验组次。
进一步地，本实施例中确定模型的试验范围，包括以下内容：
模型的试验范围包括岸线、岸滩和岸滩外海边界；
根据试验要求确定模型所需容纳的岸线和岸滩，岸滩外海边界应超出泥沙起动水深；
确定定床边界位置和动床边界位置，通过以下步骤实现：
动床范围以预设的研究区域为中心，动床边界设置在研究区域的岸滩演变范围之外，动床外海边界应满足水深大于泥沙起动水深的条件；
设定物理模型试验范围中除去动床范围的部分为定床范围，定床沿岸上下游边界与动床沿岸上下游边界间设置过渡段，过渡段长度大于5倍模型中试验波长。
进一步地，本实施例中确定原型海岸的代表波周期，通过以下步骤实现：
根据当地波高与波周期关系，确定代表波周期；
若无法获得当地波高与波周期关系，则按下式确定代表波周期：
Figure BDA0003809120830000101
式中：Tr为代表波周期；Ti为波高大于起动波高的第i级有效波高对应的波周期，Hi为波高大于起动波高的第i级有效波高，Pi为为波高大于起动波高的第i级有效波高对应的频率；
确定原型海岸的代表波向和原型海岸的代表波高，通过以下步骤实现：
当若有原型海岸波浪测站的现场波浪观测资料，则代表波向和代表波高按下式确定：
Figure BDA0003809120830000102
Figure BDA0003809120830000103
式中：αr为代表波向，Hr为代表波高；Hi为波高大于起动波高的第i级有效波高，Pi为波高大于起动波高的第i级有效波高对应的频率，αi为波高大于起动波高的第i级有效波高对应的波向角；
若没有原型海岸波浪测站的现场波浪观测资料，则根据原型海岸波浪测站的深水波浪观测资料或后报的再分析资料，通过公式(2)和(3)分别计算代表波向和代表波高，并根据波浪浅水变形将深水代表波推算至原型海岸近岸浅水区，获得研究区域的代表波向和代表波高；
确定原型海岸总质量沿岸输沙率，通过以下步骤实现：
若有原型沿岸输沙率实测资料，则根据原型沿岸输沙率实测资料，确定原型海岸的总质量沿岸输沙率；
若无原型沿岸输沙率实测资料，则根据下式计算原型海岸的总质量沿岸输沙率：
Figure BDA0003809120830000111
式中，Qm为总质量沿岸输沙率；ρs为泥沙颗粒密度，ρ为水的密度；m为岸滩坡度；Ir为波浪破碎因子，
Figure BDA0003809120830000112
其中Hb为破波波高，L0为深水波长；ω为泥沙沉速；Eb为破碎处的波能，Eb＝0.125ρgHb 2，其中g为重力加速度；Cgb为破碎处的波群速度；αb为破波波向角；umb为破波底质点最大速度，
Figure BDA0003809120830000113
其中，hb为破波水深；Hb采用均方根波高。
进一步地，本实施例中确定平面比尺、垂直比尺、波长比尺、波高比尺和波周期比尺，通过以下步骤实现：
基于试验场地尺寸、研究区域范围和造波机造波能力，确定平面比尺λl和垂直比尺λh；
当若试验要求的波浪绕射相似强于波浪折射相似，则平面比尺λl等于垂直比尺λh，波长比尺λL和波高比尺λH均等于垂直比尺λh；
按下式计算波周期比尺λT：
λT＝λh 1/2 (29)
泥沙的起动波高比尺等于模型的垂直比尺；
确定泥沙的沉速比尺，通过以下步骤实现：
基于悬沙对流输运相似和悬浮相似，确定泥沙的沉速比尺；
悬沙对流输运相似对应的沉速比尺为：
λω＝λh 3/2/λl (30)
式中：λω为沉速比尺；λh为垂直比尺；λl为平面比尺；
悬浮相似对应的沉速比尺为：
λω＝λh/λl 1/2 (31)
若模型为正态模型，则同时满足悬沙对流输运相似和悬浮相似；
若模型为变态物理模型且试验要求的悬浮相似强于淤积部位相似，则采用公式(5)确定泥沙的沉速比尺；
若模型为变态物理模型且要求淤积部位相似，则采用公式(6)确定泥沙的沉速比尺。
进一步地，本实施例中根据泥沙的起动波高比尺和沉速比尺，选择模型沙并确定模型沙的中值粒径、模型沙的级配、模型沙的颗粒密度和模型沙的表观干密度，通过以下步骤实现：
步骤4.1：选取刘家驹提出的波浪作用下的泥沙起动公式作为泥沙起动波高计算依据，得到模型沙的中值粒径、相对浮密度与垂直比尺的关系：
Figure BDA0003809120830000121
式中：Δ为相对浮密度，Δ＝(ρs/ρ)-1，其中ρs为泥沙颗粒密度，ρ为水的密度；λd为泥沙粒径比尺，d为中值粒径；λh为垂直比尺；
Figure BDA0003809120830000122
为在原型沙中Δgd+0.486×10-6/d的值和在模型沙中Δgd+0.486×10-6/d的值之比；
步骤4.2：选取武汉水利电力学院提出的泥沙沉速公式作为泥沙沉速计算依据，得到模型沙的中值粒径、相对浮密度与垂直比尺的关系：
Figure BDA0003809120830000123
式中：
Figure BDA0003809120830000124
为在原型海岸中
Figure BDA0003809120830000125
的值和模型海岸中
Figure BDA0003809120830000126
的值之比，其中υ为水的运动粘滞系数，λω为沉速比尺；
步骤4.3：若要求淤积部位相似，则分别根据公式(6)、(8)和(9)计算粒径比尺，确定模型沙的中值粒径和模型沙的颗粒密度；
若要求悬浮相似，则分别根据公式(7)、(8)和(9)计算粒径比尺，确定模型沙的中值粒径和模型沙的颗粒密度；
步骤4.4：根据下式确定模型沙的级配：
Figure BDA0003809120830000127
式中：di表示粒径累计频率等于i时对应的粒径；
Figure BDA0003809120830000128
为在原型沙di/d的值和模型沙中di/d的值之比；
模型沙的表观干密度由模型沙干燥后使用环刀称重法测得。
进一步地，本实施例中根据公式(4)，计算模型预估总质量沿岸输沙率；
在试验水池中，制作定床地形和动床地形，通过以下步骤实现：
根据原型海岸实测地形，按照平面比尺和垂直比尺制作定床地形和动床地形；
定床地形采用米沙、砖块或石子填充后进行水泥抹面，动床地形采用在以米沙、砖块或石子衬底上铺模型沙的形式制作，模型沙厚度超过原型海岸历史最大冲刷深度经垂直比尺换算后的模型冲刷深度；
在动床地形和定床地形交接处设置过渡带，通过以下步骤实现：
过渡带设置在动床地形和定床地形接触处，设置过渡带的方法为使用胶水将一层模型沙黏贴在定床地形表面；
若原型海岸上有涉水建筑物，则在试验水池中相应的制作涉水建筑物地形，通过以下步骤实现：
涉水建筑物根据垂直比尺和正态模型，采用不易变形材料制作涉水建筑物并进行水泥抹面，建筑物地形边坡坡度大于1:1的部分按1:1制作；
动床上游边界设置加沙装置，加沙装置包括PVC圆管、阀门、电磁流量计、水泵、水槽、搅拌器以及连接水管，水槽中存放模型沙与水混合形成的水沙混合液，PVC圆管、阀门、电磁流量计、水泵、水槽、搅拌器以及连接水管水平放置，PVC圆管、阀门、电磁流量计和水泵由连接水管按顺序连通，PVC圆管连通水槽，PVC圆管两端封闭，PVC圆管面向模型中心区域的一面均匀开设若干个圆孔，PVC圆管垂直输沙方向布置；
进一步地，本实施例中加沙装置的工作方法，通过以下步骤实现：
首先将模型沙和水以一定比例放置在水槽中；
通过搅拌器搅拌形成一定浓度的水沙混合液；
通过水泵将水槽中的水沙混合液抽至PVC圆管中；
水沙混合液通过PVC圆管上的圆孔流进试验水池的动床地形的上游；
加沙装置工作时，水槽中水沙混合液中悬沙浓度与水泵抽水流速的关系满足下式：
CSUadd＝Qadd (35)
式中，C为水槽中水沙混合液的悬沙浓度；S为连接水泵与PVC圆管的连接水管横断面过水面积；Uadd为连接水泵与PVC圆管的连接水管中水流平均流速；Qadd为试验加沙率；
动床下游边界设置收沙装置，收沙装置为一封底不封口且具有梯形横断面的槽体，收沙装置埋置在动床输沙方向的末端，收沙装置的封口与动床地形表面齐平。
进一步地，本实施例中开展模型验证试验，通过以下步骤实现：
步骤8.1：将试验水池加水到设定的试验水位；
步骤8.2：待试验水池的水面稳定后，设置造波机造出不使动床泥沙起动的波浪，造波机运行一段时间，直到模型沙密实后关闭造波机；
步骤8.3：待试验水池的水面稳定后，按代表波要素设置造波机参数，造波机开始造波；
步骤8.4：待造波机造出的波浪稳定后，加沙装置开始加沙，其中初次加沙率取值为模型预估总质量沿岸输沙率；
步骤8.5：运行时间tv后，测量模型沿岸输沙率和模型沿岸收沙率，关闭造波机和加沙装置；
待试验水池的水面稳定后，排干试验水池中的水，测量动床地形，获得试验结果；
试验结果包括模型沿岸输沙率和岸滩冲淤情况，岸滩冲淤情况由开展模型验证试验后的动床地形与开展模型验证试验前动床地形相减得到；
检验试验结果是否满足设定的验证要求，恢复初始的涉水建筑物地形、人为地形、动床地形；
步骤8.6：若步骤8.5中试验结果未达到验证要求，则通过调节加沙率和代表波要素，重复步骤8.1至步骤8.5，直至试验结果满足验证要求；
步骤8.3中按代表波要素设置造波机时，若采用不规则波，则采用适合当地波浪条件的波浪谱，并设置均方根波高等于代表波高，谱峰周期等于代表波周期；
步骤8.5中时间tv由公式tv＝tpv/λt确定，其中tpv为模型验证时间tv对应的原型岸滩演变时间；λt为岸滩冲淤时间比尺；
步骤8.5中检验试验结果是否满足验证要求，通过以下步骤实现：
(a)对于原型海岸地形稳定的岸滩，验证模型验证试验获得的模型沿岸输沙率是否达到平衡以及岸滩冲淤情况是否稳定，若模型验证试验获得的模型沿岸输沙率达到平衡且岸滩冲淤情况稳定，则认为试验结果满足验证要求；
(b)对于原型海岸地形不稳定的岸滩，则验证对应的实际岸滩演变过程，若模型验证试验中岸滩演变与实际岸滩演变过程一致，则认为试验结果满足验证要求。
进一步地，本实施例中开展试验要求的岸滩演变试验组次，通过以下步骤实现：
步骤10.1，将动床地形设置为试验要求的初始地形；
步骤10.2，将试验水池加水到试验要求的水位；
待试验水池的水面稳定后，设置造波机造出不使动床泥沙起动的波浪，造波机运行一段时间，直到模型沙密实后关闭造波机；
步骤10.3，待试验水池的水面稳定后，按试验波浪要求设置造波机参数，造波机开始造波；
待造波机造出的波浪稳定后，加沙装置开始加沙，加沙率取值为模型验证试验中最后实测的模型沿岸输沙率；
步骤10.4，运行时间tt后，关闭造波机和加沙装置，排干试验水池中的水，测量动床地形；
步骤10.5，如试验尚未完成，加水至试验要求水位，重复步骤10.3和步骤10.4直至该组试验完成；
步骤10.4中时间tt由试验要求的原型岸滩演变典型时间除以岸滩冲淤时间比尺得到。
进一步地，本实施例中按如下公式计算岸滩冲淤时间比尺：
Figure BDA0003809120830000151
式中，λt为岸滩冲淤时间比尺，
Figure BDA0003809120830000152
为即原型沙表观干密度与模型沙表观干密度之比；
Figure BDA0003809120830000153
为实际的总质量沿岸输沙率比尺。
试验要求包括预先设定的波浪要素、波浪方向、模型所需容纳的岸线和岸滩、悬浮相似、初始地形、水位和原型岸滩演变典型时间；人为地形包括航道和挖槽。
加沙装置8、收沙装置9、造波机10、导波板11、PVC圆管12、阀门13、电磁流量计14、水泵15、水槽16、连接水管17和搅拌器18上述部件在现有技术中可采用的型号很多，本领域技术人员可根据实际需求选用合适的型号，本实施例不再一一举例。
实施例二
某沙质海岸面向开敞大西洋，海域波浪强、潮流弱。近岸波浪破碎及其引起的沿岸流较强，是沿岸泥沙运动的主要动力。沿岸泥沙运动以悬移质输运为主，主要发生在在斜向波作用下的波浪破碎带内。在季节性波浪作用下，海滩剖面存在明显的季节性变化。但多年实测资料表明，长期来看，岸滩比较稳定。岸线多年累计进退幅度不足15m。现拟在该海岸段建设挖入式港池，并采用防浪挡沙堤和T型防护丁坝掩护进出航道。港口、航道、防浪挡沙堤和T型防护丁坝的建设将打破该海域动力泥沙与水下地形的平衡，引起泥沙的冲刷与淤积，致使岸滩地形在年计的时间尺度内发生较大变化，进而可能对港口的营运和周边地区的人民亲水生活及生命安全造成威胁。因此，采用科学的手段预测港口、航道、防浪挡沙堤和T型防护丁坝建设后岸滩的稳定性及泥沙冲淤情况极为必要。
下面结合该实施例说明如何利用本发明方法来开展沙质海岸岸滩演变物理模型试验，并预测拟建港口、航道、防浪挡沙堤和T型防护丁坝建设3年和9年后岸滩的稳定性和泥沙冲淤情况。基本实施步骤如下。
(1)该实施例旨在研究港口1、航道2、防浪挡沙堤3和T型防护丁坝4建设后岸滩的稳定性及泥沙冲刷和淤积情况。根据邻近工程经验，类似规模的工程引起的泥沙冲淤影响大致在工程上下游各5km范围内。水下地形资料显示，11m水深以浅区域坡度较陡，平均坡度约为1/54，水深11米以外的深水区水下地形坡度大幅减缓，平均坡度仅约1/391。根据当地水文泥沙气象条件分析，原型海岸泥沙起动水深约10.5m～12m，工程海域泥沙运动和岸滩演变主要限于12m水深以浅区域。
根据上述分析，确定模型的试验范围为：以拟建港口1为中心，包括港口上下游岸线5.6公里，模型外海边界设置在水深16m附近，平均离岸距离为5.75公里。动床地形范围为：以拟建港口1为中心，动床地形的上游边界5和动床地形的下游边界6分别定在港口上下游5.0公里(对应模型上下游20m)处，动床外海边界7设置在水深12m处。模型试验范围内除去动床的部分即为定床模型范围。定床沿岸上下游边界与动床上下游边界间设置有过渡段，长度大于5倍模型中试验波长。模型试验范围和动床边界位置如图1所示。
(2)该研究海域无波浪测站，亦无当地波高与波周期关系。根据研究区域附近的深水波浪后报资料获得海域附近水深100米处的波浪资料。进一步根据公式(1)、(2)和(3)计算深水代表波向、深水代表波高和代表波周期，深水代表波向、深水代表波高和代表波周期分别等于303°、2.05m和8.2s。
若没有原型海岸波浪测站的现场波浪观测资料，则根据原型海岸波浪测站的深水波浪观测资料或后报的再分析资料，通过公式(2)和(3)分别计算代表波向和代表波高，并根据波浪浅水变形将深水代表波推算至原型海岸近岸浅水区，获得研究区域的代表波向和代表波高；研究海域的代表波向为301°、代表波高为1.99m、代表波周期为8.2s。
根据原型沿岸实测资料，直接获取泥沙中值粒径、粒径级配、泥沙颗粒密度和泥沙表观干密度；
具体地，岸滩泥沙中值粒径为0.23mm(即0.0023m)，粒径级配如表1。泥沙颗粒密度为2650kg/m3，泥沙表观干密度为1590kg/m3。总体积沿岸输沙率约为100万方每年，即0.032m3/s，对应的总质量沿岸输沙率为159万t/a，即50.4kg/s。
表1原型海岸泥沙粒径级配表
粒径累计频率(％) 95 75 50 25 5
粒径(m) 0.00031 0.00026 0.00023 0.00019 0.00013
(3)根据试验场地尺寸(试验港池长45m、宽35m)、模型的试验范围及造波机能力，确定模型的平面比尺λl＝250。试验未要求波浪绕射相似强于波浪折射相似，考虑模型水深不宜过小，取模型变率ξ＝2，垂直比尺λh＝125。
波长比尺λL和波高比尺λH均等于垂直比尺λh，即λL＝λH＝125，确定波周期比尺λT＝λh 1/2＝11.18。
(4)取泥沙起动波高比尺λHcr等于垂直比尺λh，即λHcr＝λh＝125。试验未要求悬浮相似强于淤积部位相似，泥沙沉速比尺根据公式λω＝λh 3/2/λl计算为5.59。
拟选择颗粒密度为1270～1400kg/m3的煤粉作为模型沙，若要求淤积部位相似，根据公式(6)、(8)和(9)，计算得到粒径比尺在0.56和1.74之间。若要求悬浮相似，则分别根据根据公式(6)和(9)，计算得到粒径比尺在1.04和1.27之间。进一步考虑模型沙资源，选择粒径比尺为1.17。
综上，采用某精制煤粉作为本实施例物理模型试验的模型沙，煤粉颗粒密度为1350kg/m3，中值粒径为0.20mm，对应的泥沙粒径比尺为1.17。进一步根据公式(10)，获得模型沙级配，见表2。
表2模型沙粒径级配表
粒径累计频率(％) 95 75 50 25 5
粒径(m) 0.00026 0.00022 0.0002 0.00016 0.00011
煤粉表观干密度由干燥后的煤粉使用环刀称重法测得，为870kg/m3。
(5)根据公式(4)，预估模型总质量沿岸输沙率Qmm为0.0061kg/s。
(6)在试验水池中，根据实测地形资料按照平面比尺和垂直比尺制作定床地形和动床地形，其中定床地形部分采用砖块或石子填充并进行水泥抹面，动床地形部分采用以砖块衬底上铺模型沙的形式制作。通过胶水将一薄层的模型沙黏贴在连接动床部分的定床表面上作为定动床过渡带，宽度1.5m，大于5倍模型试验中的波长。
模型的试验范围和初始动床地形示意图如图1所示。在动床上下游分别设置加沙装置8和收沙装置9，按代表波方向(301°)在外海边界处设置不规则造波机10，同时在造波机两侧设置与其垂直的导波板11，以防止波浪能量耗散。
收沙装置9为一封底不封口具有梯形横断面的水槽，埋置在动床地形输沙方向的末端，垂直输沙方向布置，上开口与动床表面齐平。
加沙装置8如图2所示，由水平放置的PVC圆管12、阀门13、电磁流量计14、水泵15、存放模型沙与水混合浊液的水槽16、连接水管17及搅拌器18组成。其中PVC圆管12、阀门13、电磁流量计14、水泵15由一个连接水管17按顺序连通；另一个连接水管17一端连接PVC圆管12，连接水管17另一端伸入水槽16。PVC圆管12两端封闭，面向模型中心区域的一面均匀布置一系列圆孔，PVC圆管垂直输沙方向布置。
(7)开展试验前的波浪率定工作，首先将试验水池加水到试验水位，待水面稳定后，设置不规则造波机参数，采用JONSWAP谱作为不规则波波浪谱，谱峰升高因子取3.3；当圆频率小于等于谱峰圆频率时，谱峰参数取0.07，当圆频率大于谱峰圆频率时，谱峰参数取0.09；均方根波高取代表波高对应的模型值(0.0164m)；谱峰周期取代表波周期对应的模型值(0.73s)，开启造波机，调试设备使水池中波浪要素满足试验要求。
(8)尽管海滩剖面存在明显的季节性变化，但实测资料表明，长期来看，岸滩比较稳定。因此，该实施例中将验证沿岸输沙率。连续2年实测资料显示年平均总体积沿岸输沙率为100万方每年，即0.032m3/s，对应的总质量沿岸输沙率为159万t/a，即50.4kg/s。根据原型海岸总质量沿岸输沙率和模型预估总质量沿岸输沙率以及公式
Figure BDA0003809120830000181
计算得到模型对应的泥沙冲淤时间比尺为λt＝1728，即模型运行10.14h相当于原型海岸演变2年。实施例中验证试验时长tv取10.14个小时，验证试验按如下程序开展：
(a)首先将试验水池加水到试验水位；
(b)待水面稳定后，设置造波机造某一不使动床泥沙起动的波浪，运行一段时间使模型沙充分密实后，关闭造波机；
(c)待试验水池的水面稳定后，按代表波要素设置造波机参数：
采用JONSWAP谱作为不规则波波浪谱，谱峰升高因子取3.3；当圆频率小于等于谱峰圆频率时，谱峰参数取0.07，当圆频率大于谱峰圆频率时，谱峰参数取0.09；均方根波高取代表波高对应的模型值(0.0164m)；谱峰周期取代表波周期对应的模型值(0.73s)，开启造波机开始造波；
(d)待波浪稳定后，启动加沙装置开始加沙，初次加沙率取模型预估总质量沿岸输沙率，即0.0061kg/s；
(e)运行时间10.14h后，测量模型沿岸输沙率和模型沿岸收沙率，关闭造波机和加沙装置，待水面平静后，缓慢排干试验水池中的水，测量动床地形，检验岸滩冲淤情况是否稳定，检验模型沿岸输沙率是否平衡，然后恢复初始涉水建筑物地形、人为地形和动床地形；
(f)通过调节加沙率和代表波要素，重复过程(a)-(e)，直至当调节均方根波高为0.02m、谱峰周期为0.73s、加沙率为0.0083kg/s时，发现岸滩冲淤情况稳定、模型沿岸输沙率达到平衡，此时验证试验符合要求。
试验中加沙装置的工作方法为首先将模型沙和水以一定比例放置在水槽中，通过搅拌器不停搅拌形成一定浓度的水沙混合液，然后通过水泵将水槽中的水沙混合液抽至PVC圆管中，并通过圆管上的小孔流进试验水池的动床上游。水槽中水沙混合液中悬沙浓度与水泵抽水流速的关系应满足下式：
CSUadd＝Qadd
式中：C为水槽中水沙混合液的悬沙浓度；S为连接水泵与PVC圆管的水管横断面过水面积；Uadd为连接水管中水流平均流速；Qadd为试验加沙率。
(9)验证试验中最后实测的模型沿岸输沙率为0.0083kg/s，根据公式
Figure BDA0003809120830000191
计算岸滩冲淤时间比尺为2351。
(10)在完成验证的基础上，进一步开展试验以预测在与验证资料相同动力条件下拟建港口、航道及防浪挡沙堤建设后3和9年后该段岸滩的稳定性及泥沙冲刷和淤积情况。本实施例中代表波参数同验证试验中最后取的代表波参数，即波向为301°、均方根波高为2.0cm、谱峰周期为0.73s。根据岸滩冲淤时间比尺2351，3年和9年对应的模型试验时间分别为11.18h和33.53h。试验程序如下。
(A)将动床地形设置为试验需要的初始地形，同时制作港池、航道及防浪挡沙堤，其中防浪挡沙堤按垂直比尺依照正态模型采用不易变形材料制作并进行表面糙率处理，且坡度大于1:1的地方按1:1制作，按波浪方向(301°)在外海边界处设置造波机；
(B)将试验水池加水到试验要求水位，待水面稳定后，设置造波机造某一小波浪(不使动床泥沙起动)，运行一段时间待模型沙密实后，关闭造波机；
(C)待水面稳定后，设置不规则造波机参数，采用JONSWAP谱作为不规则波波浪谱，谱峰升高因子取3.3；当圆频率小于等于谱峰圆频率时，谱峰参数取0.07，当圆频率大于谱峰圆频率时，谱峰参数取0.09；均方根波高取0.02m；谱峰周期0.73s；开启造波机开始造波，待波浪稳定后，启动加沙装置开始加沙，加沙率取0.0083kg/s，加沙装置的工作方法同验证试验中的加沙方法；
(D)运行时间11.18h后，关闭造波机和加沙装置，排干水池中的水，测量动床地形，该地形即为拟建港口、航道及防浪挡沙堤影响下岸滩演变3年后的预测地形，见图3。运行后的动床地形减去初始地形，即可得到3年后与3年前相比的泥沙冲淤分布。
重复(C)和(D，其中(D)中运行时间为22.36h(33.53-11.18＝22.36)。(D)中所测地形即为拟建港口、航道及防浪挡沙堤影响下岸滩演变9年后的预测地形，见图4。该动床地形减去初始地形，即可得到9年后与9年前相比的泥沙冲淤分布。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。

Claims (10)

1.一种波浪作用下沙质海岸岸滩演变物理模型试验方法，其特征在于，包括：
(1)确定模型的试验范围、定床边界位置和动床边界位置；
(2)确定原型海岸的代表波向、代表波高、代表波周期、泥沙中值粒径、粒径级配、泥沙颗粒密度、泥沙表观干密度以及总质量沿岸输沙率；
(3)确定模型的平面比尺、垂直比尺、波长比尺、波高比尺、波周期比尺、泥沙起动波高比尺和沉速比尺；
(4)根据泥沙的起动波高比尺和沉速比尺，选择模型沙并确定模型沙的中值粒径、模型沙的级配、模型沙的颗粒密度和模型沙的表观干密度；
(5)计算模型预估总质量沿岸输沙率；
(6)在试验水池中，制作定床地形和动床地形，并在动床地形和定床地形交接处设置过渡带；
若定床地形或动床地形上有涉水建筑物和人为地形，则在试验水池中相应地制作涉水建筑物地形和人为地形；
在动床地形的上游边界处设置加沙装置，在动床地形的下游边界处设置收沙装置，并按试验要求的波浪方向在外海边界处设置造波机；
在造波机两侧设置与造波机垂直的导波板；
(7)开展波浪率定工作：首先将试验水池加水到试验水位；
待试验水池的水面稳定后，根据代表波要素设置波浪参数，开启造波机造波，调试造波机使试验水池中波浪要素满足试验要求；
(8)开展模型验证试验；
(9)根据模型验证试验中最后实测的模型沿岸输沙率，计算岸滩冲淤时间比尺；
(10)开展试验要求的岸滩演变试验组次。
2.根据权利要求1所述的一种波浪作用下沙质海岸岸滩演变物理模型试验方法，其特征在于，确定模型的试验范围，包括以下内容：
模型的试验范围包括岸线、岸滩和岸滩外海边界；
根据试验要求确定模型所需容纳的岸线和岸滩，岸滩外海边界应超出泥沙起动水深；
确定定床边界位置和动床边界位置，通过以下步骤实现：
动床范围以预设的研究区域为中心，动床边界设置在研究区域的岸滩演变范围之外，动床外海边界应满足水深大于泥沙起动水深的条件；
设定物理模型试验范围中除去动床范围的部分为定床范围，定床沿岸上下游边界与动床沿岸上下游边界间设置过渡段，过渡段长度大于5倍模型中试验波长。
3.根据权利要求1所述的一种波浪作用下沙质海岸岸滩演变物理模型试验方法，其特征在于，确定原型海岸的代表波周期，通过以下步骤实现：
根据当地波高与波周期关系，确定代表波周期；
若无法获得当地波高与波周期关系，则按下式确定代表波周期：
Figure FDA0003809120820000021
式中：Tr为代表波周期；Ti为波高大于起动波高的第i级有效波高对应的波周期，Hi为波高大于起动波高的第i级有效波高，Pi为为波高大于起动波高的第i级有效波高对应的频率；
确定原型海岸的代表波向和原型海岸的代表波高，通过以下步骤实现：
若有原型海岸波浪测站的现场波浪观测资料，则代表波向和代表波高按下式确定：
Figure FDA0003809120820000022
Figure FDA0003809120820000023
式中：αr为代表波向，Hr为代表波高；Hi为波高大于起动波高的第i级有效波高，Pi为波高大于起动波高的第i级有效波高对应的频率，αi为波高大于起动波高的第i级有效波高对应的波向角；
若没有原型海岸波浪测站的现场波浪观测资料，则根据原型海岸波浪测站的深水波浪观测资料或后报的再分析资料，通过公式(2)和(3)分别计算代表波向和代表波高，并根据波浪浅水变形将深水代表波推算至原型海岸近岸浅水区，获得研究区域的代表波向和代表波高；
确定原型海岸总质量沿岸输沙率，通过以下步骤实现：
若有原型沿岸输沙率实测资料，则根据原型沿岸输沙率实测资料，确定原型海岸的总质量沿岸输沙率；
若无原型沿岸输沙率实测资料，则根据下式计算原型海岸的总质量沿岸输沙率：
Figure FDA0003809120820000024
式中，Qm为总质量沿岸输沙率；ρs为泥沙颗粒密度，ρ为水的密度；m为岸滩坡度；Ir为波浪破碎因子，
Figure FDA0003809120820000031
其中Hb为破波波高，L0为深水波长；ω为泥沙沉速；Eb为破碎处的波能，Eb＝0.125ρgHb 2，其中g为重力加速度；Cgb为破碎处的波群速度；αb为破波波向角；umb为破波底质点最大速度，
Figure FDA0003809120820000032
其中，hb为破波水深；Hb采用均方根波高。
4.根据权利要求1所述的一种波浪作用下沙质海岸岸滩演变物理模型试验方法，其特征在于，确定平面比尺、垂直比尺、波长比尺、波高比尺和波周期比尺，通过以下步骤实现：
基于试验场地尺寸、研究区域范围和造波机造波能力，确定平面比尺λl和垂直比尺λh；
若试验要求的波浪绕射相似强于波浪折射相似，则平面比尺λl等于垂直比尺λh，波长比尺λL和波高比尺λH均等于垂直比尺λh；
按下式计算波周期比尺λT：
λT＝λh 1/2 (5)
泥沙的起动波高比尺等于模型的垂直比尺；
确定泥沙的沉速比尺，通过以下步骤实现：
基于悬沙对流输运相似和悬浮相似，确定泥沙的沉速比尺；
悬沙对流输运相似对应的沉速比尺为：
λω＝λh 3/2/λl (6)
式中：λω为沉速比尺；λh为垂直比尺；λl为平面比尺；
悬浮相似对应的沉速比尺为：
λω＝λh/λl 1/2 (7)
若模型为正态模型，则同时满足悬沙对流输运相似和悬浮相似；
若模型为变态物理模型且试验要求的悬浮相似强于淤积部位相似，则采用公式(5)确定泥沙的沉速比尺；
若模型为变态物理模型且要求淤积部位相似，则采用公式(6)确定泥沙的沉速比尺。
5.根据权利要求1所述的一种波浪作用下沙质海岸岸滩演变物理模型试验方法，其特征在于，根据泥沙的起动波高比尺和沉速比尺，选择模型沙并确定模型沙的中值粒径、模型沙的级配、模型沙的颗粒密度和模型沙的表观干密度，通过以下步骤实现：
步骤5.1：选取刘家驹提出的波浪作用下的泥沙起动公式作为泥沙起动波高计算依据，得到模型沙的中值粒径、相对浮密度与垂直比尺的关系：
Figure FDA0003809120820000041
式中：Δ为相对浮密度，Δ＝(ρs/ρ)-1，其中ρs为泥沙颗粒密度，ρ为水的密度；λd为泥沙粒径比尺，d为中值粒径；λh为垂直比尺；
Figure FDA0003809120820000042
为在原型沙中Δgd+0.486×10-6/d的值和在模型沙中Δgd+0.486×10-6/d的值之比；
步骤5.2：选取武汉水利电力学院提出的泥沙沉速公式作为泥沙沉速计算依据，得到模型的中值粒径、相对浮密度与垂直比尺的关系：
Figure FDA0003809120820000043
式中：
Figure FDA0003809120820000044
为在原型海岸中
Figure FDA0003809120820000045
的值和模型海岸中
Figure FDA0003809120820000046
的值之比，其中υ为水的运动粘滞系数，λω为沉速比尺；
步骤5.3：若要求淤积部位相似，则分别根据公式(6)、(8)和(9)计算粒径比尺，确定模型沙的中值粒径和模型沙的颗粒密度；
若要求悬浮相似，则分别根据公式(7)、(8)和(9)计算粒径比尺，确定模型沙的中值粒径和模型沙的颗粒密度；
步骤5.4：根据下式确定模型沙的级配：
Figure FDA0003809120820000047
式中：di表示粒径累计频率等于i时对应的粒径；
Figure FDA0003809120820000048
为在原型沙di/d的值和模型沙中di/d的值之比；
模型沙的表观干密度由模型沙干燥后使用环刀称重法测得。
6.根据权利要求3所述的一种波浪作用下沙质海岸岸滩演变物理模型试验方法，其特征在于，根据公式(4)，计算模型预估总质量沿岸输沙率；
在试验水池中，制作定床地形和动床地形，通过以下步骤实现：
根据原型海岸实测地形，按照平面比尺和垂直比尺制作定床地形和动床地形；
定床地形采用米沙、砖块或石子填充后进行水泥抹面，动床地形采用在以米沙、砖块或石子衬底上铺模型沙的形式制作，模型沙厚度超过原型海岸历史最大冲刷深度经垂直比尺换算后的模型冲刷深度；
在动床地形和定床地形交接处设置过渡带，通过以下步骤实现：
过渡带设置在动床地形和定床地形接触处，设置过渡带的方法为使用胶水将一层模型沙黏贴在定床地形表面；
若原型海岸上有涉水建筑物，则在试验水池中相应的制作涉水建筑物地形，通过以下步骤实现：
涉水建筑物根据垂直比尺和正态模型，采用不易变形材料制作涉水建筑物并进行水泥抹面，建筑物地形边坡坡度大于1:1的部分按1:1制作；
动床上游边界设置加沙装置，加沙装置包括PVC圆管、阀门、电磁流量计、水泵、水槽、搅拌器以及连接水管，水槽中存放模型沙与水混合形成的水沙混合液，PVC圆管、阀门、电磁流量计、水泵、水槽、搅拌器以及连接水管水平放置，PVC圆管、阀门、电磁流量计和水泵由连接水管按顺序连通，PVC圆管连通水槽，PVC圆管两端封闭，PVC圆管面向模型中心区域的一面均匀开设若干个圆孔，PVC圆管垂直输沙方向布置。
7.根据权利要求1所述的一种波浪作用下沙质海岸岸滩演变物理模型试验方法，其特征在于，加沙装置的工作方法，通过以下步骤实现：
首先将模型沙和水以一定比例放置在水槽中；
通过搅拌器搅拌形成一定浓度的水沙混合液；
通过水泵将水槽中的水沙混合液抽至PVC圆管中；
水沙混合液通过PVC圆管上的圆孔流进试验水池的动床地形的上游；
加沙装置工作时，水槽中水沙混合液中悬沙浓度与水泵抽水流速的关系满足下式：
CSUadd＝Qadd (11)
式中，C为水槽中水沙混合液的悬沙浓度；S为连接水泵与PVC圆管的连接水管横断面过水面积；Uadd为连接水泵与PVC圆管的连接水管中水流平均流速；Qadd为试验加沙率；
动床下游边界设置收沙装置，收沙装置为一封底不封口且具有梯形横断面的槽体，收沙装置埋置在动床输沙方向的末端，收沙装置的封口与动床地形表面齐平。
8.根据权利要求1所述的一种波浪作用下沙质海岸岸滩演变物理模型试验方法，其特征在于，开展模型验证试验，通过以下步骤实现：
步骤8.1：将试验水池加水到设定的试验水位；
步骤8.2：待试验水池的水面稳定后，设置造波机造出不使动床泥沙起动的波浪，造波机运行一段时间，直到模型沙密实后关闭造波机；
步骤8.3：待试验水池的水面稳定后，按代表波要素设置造波机参数，造波机开始造波；
步骤8.4：待造波机造出的波浪稳定后，加沙装置开始加沙，其中初次加沙率取值为模型预估总质量沿岸输沙率；
步骤8.5：运行时间tv后，测量模型沿岸输沙率和模型沿岸收沙率，关闭造波机和加沙装置；
待试验水池的水面稳定后，排干试验水池中的水，测量动床地形，获得试验结果；
试验结果包括模型沿岸输沙率和岸滩冲淤情况，岸滩冲淤情况由开展模型验证试验后的动床地形与开展模型验证试验前动床地形相减得到；
检验试验结果是否满足设定的验证要求，恢复初始的涉水建筑物地形、人为地形、动床地形；
步骤8.6：若步骤8.5中试验结果未达到验证要求，则通过调节加沙率和代表波要素，重复步骤8.1至步骤8.5，直至试验结果满足验证要求；
步骤8.3中按代表波要素设置造波机时，若采用不规则波，则采用适合当地波浪条件的波浪谱，并设置均方根波高等于代表波高，谱峰周期等于代表波周期；
步骤8.5中时间tv由公式tv＝tpv/λt确定，其中tpv为模型验证时间tv对应的原型岸滩演变时间；λt为岸滩冲淤时间比尺；
步骤8.5中检验试验结果是否满足验证要求，通过以下步骤实现：
(a)对于原型海岸地形稳定的岸滩，验证模型验证试验获得的模型沿岸输沙率是否达到平衡以及岸滩冲淤情况是否稳定，若模型验证试验获得的模型沿岸输沙率达到平衡且岸滩冲淤情况稳定，则认为试验结果满足验证要求；
(b)对于原型海岸地形不稳定的岸滩，则验证对应的实际岸滩演变过程，若模型验证试验中岸滩演变与实际岸滩演变过程一致，则认为试验结果满足验证要求。
9.根据权利要求1所述的一种波浪作用下沙质海岸岸滩演变物理模型试验方法，其特征在于，开展试验要求的岸滩演变试验组次，通过以下步骤实现：
步骤9.1，将动床地形设置为试验要求的初始地形；
步骤9.2，将试验水池加水到试验要求的水位；
待试验水池的水面稳定后，设置造波机造出不使动床泥沙起动的波浪，造波机运行一段时间，直到模型沙密实后关闭造波机；
步骤9.3，待试验水池的水面稳定后，按试验波浪要求设置造波机参数，造波机开始造波；
待造波机造出的波浪稳定后，加沙装置开始加沙，加沙率取值为模型验证试验中最后实测的模型沿岸输沙率；
步骤9.4，运行时间tt后，关闭造波机和加沙装置，排干试验水池中的水，测量动床地形；
步骤9.5，如试验尚未完成，加水至试验要求水位，重复步骤10.3和步骤10.4直至该组试验完成；
步骤9.4中时间tt由试验要求的原型岸滩演变典型时间除以岸滩冲淤时间比尺得到。
10.根据权利要求1所述的一种波浪作用下沙质海岸岸滩演变物理模型试验方法，其特征在于，按如下公式计算岸滩冲淤时间比尺：
Figure FDA0003809120820000071
式中，λt为岸滩冲淤时间比尺，
Figure FDA0003809120820000072
为即原型沙表观干密度与模型沙表观干密度之比；
Figure FDA0003809120820000073
为实际的总质量沿岸输沙率比尺。
CN202211006619.4A 2022-08-22 2022-08-22 一种波浪作用下沙质海岸岸滩演变物理模型试验方法 Active CN115346427B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211006619.4A CN115346427B (zh) 2022-08-22 2022-08-22 一种波浪作用下沙质海岸岸滩演变物理模型试验方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211006619.4A CN115346427B (zh) 2022-08-22 2022-08-22 一种波浪作用下沙质海岸岸滩演变物理模型试验方法

Publications (2)

Publication Number Publication Date
CN115346427A true CN115346427A (zh) 2022-11-15
CN115346427B CN115346427B (zh) 2023-06-13

Family

ID=83953498

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211006619.4A Active CN115346427B (zh) 2022-08-22 2022-08-22 一种波浪作用下沙质海岸岸滩演变物理模型试验方法

Country Status (1)

Country Link
CN (1) CN115346427B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116451413A (zh) * 2023-01-12 2023-07-18 珠江水利委员会珠江水利科学研究院 基于水沙动力自适应特性的河口治理模拟方法
CN117852147A (zh) * 2024-01-15 2024-04-09 河海大学 一种考虑实际潮位波浪变化过程的海滩剖面冲淤方法
CN118010310A (zh) * 2024-04-10 2024-05-10 交通运输部天津水运工程科学研究所 一种模拟原型海况的波流-结构物-地基耦合作用试验方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2594969A1 (fr) * 2011-11-17 2013-05-22 IFP Energies Nouvelles Procédé d'exploration d'un bassin sédimentaire à partir d'une simulation stratigraphique
CN107421716A (zh) * 2017-09-07 2017-12-01 交通运输部天津水运工程科学研究所 水位变化下的波浪物理模型试验无间断造波装置及方法
CN112033641A (zh) * 2020-08-26 2020-12-04 大连理工大学 一种海岸地貌演变物理模型试验方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2594969A1 (fr) * 2011-11-17 2013-05-22 IFP Energies Nouvelles Procédé d'exploration d'un bassin sédimentaire à partir d'une simulation stratigraphique
CN107421716A (zh) * 2017-09-07 2017-12-01 交通运输部天津水运工程科学研究所 水位变化下的波浪物理模型试验无间断造波装置及方法
CN112033641A (zh) * 2020-08-26 2020-12-04 大连理工大学 一种海岸地貌演变物理模型试验方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
吴舸;杨振勇;: "评价绥中发电厂取水口泥沙淤积趋势的三维海岸动床模型的试验研究" *
夏益民: "沙质海岸波浪动床模型设计──毛里塔尼亚友谊港下游冲刷试验模型" *
孙林云;吴炳良;郭天润;: "波流共同作用下细沙粉沙质海岸复合沿岸输沙率计算" *
陈勇;刘家驹;: "沙质海岸突堤下游岸滩演变数值计算" *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116451413A (zh) * 2023-01-12 2023-07-18 珠江水利委员会珠江水利科学研究院 基于水沙动力自适应特性的河口治理模拟方法
CN116451413B (zh) * 2023-01-12 2023-11-14 珠江水利委员会珠江水利科学研究院 基于水沙动力自适应特性的河口治理模拟方法
CN117852147A (zh) * 2024-01-15 2024-04-09 河海大学 一种考虑实际潮位波浪变化过程的海滩剖面冲淤方法
CN118010310A (zh) * 2024-04-10 2024-05-10 交通运输部天津水运工程科学研究所 一种模拟原型海况的波流-结构物-地基耦合作用试验方法
CN118010310B (zh) * 2024-04-10 2024-06-07 交通运输部天津水运工程科学研究所 一种模拟原型海况的波流-结构物-地基耦合作用试验方法

Also Published As

Publication number Publication date
CN115346427B (zh) 2023-06-13

Similar Documents

Publication Publication Date Title
CN115346427A (zh) 一种波浪作用下沙质海岸岸滩演变物理模型试验方法
CN108256137B (zh) 一种强潮河口湾人工岛作业区港池航道回淤模拟方法
CN106096203B (zh) 一种枢纽大坝的基底施工方法
Van Rijn Estuarine and coastal sedimentation problems
Liu et al. Prediction of concerted sediment flushing
CN111289219B (zh) 一种超临界流条件下基础局部冲刷的水槽试验方法
Ghomeshi Reservoir sedimentation modelling
Caballero et al. Advances in the understanding of the hydraulic behavior of wedge-shape block spillways
CN115630507A (zh) 侵蚀基面下降后复式河道河床下切突变响应临界条件预测方法
Haghjouei et al. Experimental investigation of the diameter and length effects of the dendritic, bottomless, extended structure on reservoir sediment removal efficiency by flushing
Dam et al. Long term process-based morphological model of the Western Scheldt Estuary
Esmaeili et al. Numerical study on flushing channel evolution, case study of Dashidaira reservoir, Kurobe river
Schleiss et al. Physical model experiments on reservoir sedimentation
Bahgat Optimum use of dredged materials for sustainable shoreline management in Nile Delta
Luijendijk et al. Integrated modelling of the morphological evolution of the sand engine mega-nourishment
Karagiannis et al. Numerical simulation of scour in front of a breakwater using OpenFoam
Gonzalez et al. Mixing depth experiments on an estuarine dissipative beach; St. Georges Beach, Gironde (France)
Ding et al. Three-dimensional numerical simulation of tidal flows in the Yangtze River Estuary
CN111501656B (zh) 一种基于水泵射流的水动力提升装置及其性能测试方法和在水质提升的应用
Walstra et al. Modelling of sedimentation of dredged trenches and channels under the combined action of tidal currents and waves
Trouw et al. Design of storm return walls for the masterplan for coastal safety: from conceptual to detailed design
de Anda et al. The use of digital elevation models (dems) for bathymetry development in large tropical reservoirs
Zhu et al. Analysis on the influence of dike engineering on the discharge of artificial lake dam at the upstream end
Dong et al. A study on the riprap movement based on marine hydrodynamic model and field observation
Villaret et al. Large scale morphodynamic modeling of the Gironde estuary

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant