CN110412312B - 涌潮流速垂向分布的一种计算方法 - Google Patents
涌潮流速垂向分布的一种计算方法 Download PDFInfo
- Publication number
- CN110412312B CN110412312B CN201910480043.7A CN201910480043A CN110412312B CN 110412312 B CN110412312 B CN 110412312B CN 201910480043 A CN201910480043 A CN 201910480043A CN 110412312 B CN110412312 B CN 110412312B
- Authority
- CN
- China
- Prior art keywords
- tidal bore
- tidal
- flow velocity
- formula
- bore
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M10/00—Hydrodynamic testing; Arrangements in or on ship-testing tanks or water tunnels
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
- G01P5/24—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave
- G01P5/241—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave by using reflection of acoustical waves, i.e. Doppler-effect
Abstract
本发明公开了涌潮流速垂向分布的一种计算方法,该方法包括以下步骤:获取相关数据;判定涌潮类型;若判定为波状涌潮,则利用式(4)计算涌潮流速的垂向分布;若判定为漩滚涌潮,则利用式(8)计算涌潮流速的垂向分布。本发明以涌潮实测数据与水槽试验数据为依据,详细研究了涌潮流速在垂线上的分布特性,针对波状涌潮与漩滚涌潮在垂向分布特性的明显差异,基于数理统计方法分别提出了涌潮流速沿垂线分布的计算公式。经资料验证,该公式较好地反映了涌潮流速在垂线上的分布特征,计算误差较小。本发明可为工程设计或涌潮研究中计算垂线上各位置处的涌潮流速值提供有力工具。
Description
技术领域
本发明涉及流体力学领域,尤其涉及涌潮流速垂向分布的计算方法、组件及存储介质。
背景技术
世界上的许多强潮河口都存在涌潮现象,如巴西亚马逊河、英国的塞文河、印度的胡格里河等。涌潮存在两种典型的形态:一种是强度较小,涌潮是由一系列平行向前推进的波构成的波列,称为“波状涌潮”;另一种是强度大时,自由表面破碎,形成向前推进的水滚,称为“漩滚涌潮”或“强涌潮”。钱塘江涌潮在各河口中最为典型,是极具特色的旅游资源和自然奇观。钱塘江涌潮形成后,水流特性发生极大的改变。涌潮到达前后,水位骤然上涨2~3m,水流急速从落潮状态转为涨潮状态,并迅速达到极值,极值流速达6~10m/s,实测最大测点流速12m/s。
涌潮流速垂向分布是涌潮的重要特征,影响着河口水域的泥沙输移、河床冲淤以及涉水建筑物的安全运行,一直深受涌潮研究者的关注。随着声学剖面流速仪的使用,涌潮流速观测有了突破性的进展,先后获取了钱塘江河口、澳大利亚戴利河口、英国迪河河口等涌潮流速垂向分布的数据。但是上述成果大多是零星垂线的涌潮观测资料或室内试验数据,尚未形成具有普适意义的涌潮流速垂向分布的计算公式或方法。
发明内容
本发明的一个目的在于针对现有技术的不足,提供一种涌潮流速垂向分布的计算方法。该方法可根据判定的涌潮类型,利用已知点位的流速值计算涌潮在垂线上的流速分布。
为实现上述目的,本发明采用了以下技术方案:涌潮流速垂向分布的计算方法,包括以下步骤:
步骤一,获取相关数据。通过现场观测或试验测量的方法,获取某一涌潮的潮前水深h、潮前流速v、涌潮高度H、离河底距离为Z1处的涌潮流速V1。
步骤二,判定涌潮类型。通过涌潮Froude数(Fr)判定该涌潮是波状涌潮还是漩滚涌潮,具体方法如下:
a、计算涌潮的Froude数,即式(1),
式中,U为涌潮传播速度,v为潮前流速,h为潮前水深,H为涌潮高度。
b、根据公式(1)计算得到的涌潮Froude数,如果Fr<1.7,则判断该涌潮为波状涌潮,如果Fr≥1.7,则判断该涌潮为漩滚涌潮。
步骤三,当判定为波状涌潮时,计算涌潮流速垂向分布,具体方法如下:
a、根据已测得的离河底距离为Z1处的涌潮流速V1计算垂向最大流速Vmax,采用的是式(3),
b、根据公式(3)得到垂向最大流速Vmax,计算离河底任一距离Z(Z≤h+H)处的涌潮流速V,采用的是式(4),
步骤四,当判定为漩滚涌潮时,计算涌潮流速垂向分布,具体方法如下:
a、由已知的离河底距离Z1,计算该处的相对位置Z'1,即式(5),
b、根据上述公式(5)得到相对位置Z'1及涌潮流速V1,计算垂向最大流速Vmax,如果Z'1≤0.6,则采用的是式(6),
如果Z'1>0.6,则采用的是式(7),
c、根据式(6)或式(7)得到垂向最大流速Vmax,计算离河底任一距离Z(Z≤h+H)处的涌潮流速V,采用的是式(8),
本发明的原理是:本发明采用实测资料分析与涌潮水槽试验相结合的方法对涌潮流速在垂线上的分布特性进行了系统研究。实测资料来源于近期开展的钱塘江河口涌潮观测数据,其中,涌潮流速采用专用浮动平台加装多普勒流速剖面仪进行自动流速流向观测,记录了水平方向正北、正东方向和铅垂方向三个方向的流速,每1min记录1次,垂向上根据涨落潮期间的水深情况分5层~10层,每层间隔为0.5m。试验数据来源于涌潮玻璃水槽的试验结果,该水槽通过Bore2010涌潮测控系统生成不同潮前水深、潮前流速及涌潮高度的涌潮,采用三维声学多普勒流速仪(ADV)测量垂线上不同位置的涌潮流速。基于上述的数据资料,采用相对位置、相对流速等无量纲指标对涌潮流速沿垂线分布进行优化拟合,提出了计算公式,即式(4)和式(8)。采用该公式可利用已知点位的流速值计算出涌潮在垂线上任何点位的流速大小。
涌潮流速垂向分布的计算组件,其包括,
存储模块,用于存储实现上述方法的程序和/或参与运算的数据;
运算模块,用于实现上述方法的运算过程。
还包括,
输入输出组件,用于讲所述的参与运算的数据输入和结果的输出;
通讯模块,用于从远端接收所述的运算的数据和/或将结果传到远端。
一种存储介质,用以储存一计算机程序,该计算机程序包括多个程序代码,其用以加载至一电子装置,使得该电子装置执行本申请的方法。
本发明的有益效果是:本发明以涌潮实测数据与水槽试验数据为依据,通过数理统计方法,提出了涌潮流速沿垂线分布的计算公式,即式(4)和式(8)。该计算公式表明:波状涌潮到达时刻的垂向最大流速位于表层,其垂线分布特征来看,与一般潮汐河口的流速垂线分布特征相似;漩滚涌潮在垂向上最大流速出现在潮前静水位附近的中层,其垂线分布成抛物型曲线,与一般潮汐河口的流速垂线分布存在明显差异。经实测数据或试验数据验证,所提出的计算公式较好的反映了涌潮流速在垂线上的分布特征,计算误差较小。该方法能够作为工程设计或涌潮研究中用于计算垂线上各位置处的涌潮流速值的有力工具。
附图说明
图1为本发明的实施流程图;
图2为波状涌潮流速垂线分布计算值与实测或试验值对比图;
图3为漩滚涌潮流速垂线分布计算值与实测或试验值对比图;
具体实施方式
以下结合附图及实施例对本发明作进一步详细描述。
实施例1:涌潮流速垂向分布的计算方法,实施流程图如图1所示,包括以下步骤:
步骤一,获取相关数据。即获取计算所需的参数,采用现场观测或试验测量的方法,获取某一涌潮的潮前水深h、潮前流速v、涌潮高度H、离河底距离为Z1处的涌潮流速V1。通过历史数据的对比,进一步核实上述数据的合理性。
步骤二,判定涌潮类型。通过涌潮Froude数判定该涌潮是波状涌潮还是漩滚涌潮,具体方法如下:
a、计算涌潮的Froude数,即式(1),
式中,U为涌潮传播速度,v为潮前流速,h为潮前水深,H为涌潮高度。
b、根据公式(1)计算得到的涌潮Froude数,如果Fr<1.7,则判断该涌潮为波状涌潮,如果Fr≥1.7,则判断该涌潮为漩滚涌潮。
步骤三,当判定为波状涌潮时,计算涌潮流速垂向分布,具体方法如下:
a、根据已知的离河底距离为Z1处的涌潮流速V1计算垂向最大流速Vmax,采用的是式(3),
b、根据公式(3)得到垂向最大流速Vmax,计算离河底任一距离Z(Z≤h+H)处的涌潮流速V,采用的是式(4),
式(4)是根据波状涌潮的实测和试验所获取的垂向分布数据经拟合后得到的计算公式,其垂线分布特征是垂向最大流速位于表层,经与实测或试验数据对比,该式具有较高的计算精度,如图2所示。
步骤四,当判定为漩滚涌潮时,计算涌潮流速垂向分布,具体方法如下:
a、由已知的离河底距离Z1,计算该处的相对位置Z'1,即式(5),
b、根据上述公式(5)得到相对位置Z'1及涌潮流速V1,计算垂向最大流速Vmax,如果Z'1≤0.6,则采用的是式(6),
如果Z'1>0.6,则采用的是式(7),
c、根据式(6)或式(7)得到垂向最大流速Vmax,计算离河底任一距离Z(Z≤h+H)处的涌潮流速V,采用的是式(8),
式(8)是根据漩滚涌潮的实测和试验所获取的垂向分布数据经拟合后得到的计算公式,其垂线分布特征是垂向最大流速出现在潮前静水位附近的中层,成抛物型曲线,经与实测或试验数据对比,该公式具有较高的计算精度,如图3所示。
涌潮流速垂向分布的计算组件,其特征在于:包括,
存储模块,用于存储实现上述方法的程序和/或参与运算的数据;
运算模块,用于实现上述方法的运算过程。
还包括,
输入输出组件,用于讲所述的参与运算的数据输入和结果的输出;
通讯模块,用于从远端接收所述的运算的数据和/或将结果传到远端。
一种存储介质,用以储存一计算机程序,该计算机程序包括多个程序代码,其用以加载至一电子装置,使得该电子装置执行本实施例前述的方法。
由于方法的代码化,代码的封装及硬件的组合采用的本领域普通常规选择和实现方式,本领域普通技术人员所熟知,因此不再赘述。
实施例2:以下以涌潮实测数据进行进一步说明。
通过自记水位仪与三维声学多普勒流速仪(ADV)现场观测,获得某一涌潮的潮前水深h=1.8m、潮前流速v=-1.47m/s、涌潮高度H=2.6m、离河底距离为Z1=3.4m处的涌潮流速V1=4.15m/s。
首先,利用上述的潮前水深、涌潮高度、潮前流速,由式(1)、(2)计算得到Fr=2.02。由于Fr≥1.7,则判断该涌潮为漩滚涌潮。
接着,由已知的离河底距离Z1=3.4m及潮前水深、涌潮高度,代入式(5),计算其相对位置Z'1=0.77。由于Z'1>0.6,则将V1和Z'1代入式(7)计算得到垂向最大流速Vmax=4.18m/s。
最后,由式(8)计算得到的离河底距离Z=0.5、1.0、1.5、2.0、2.5、3.0、3.5、4.0及4.4m的流速分别为3.39、3.53、3.68、3.86、4.06、4.18、4.13、3.91及3.61m/s。
尽管已结合优选的实施例描述了本发明,然其并非用以限定本发明,任何本领域技术人员,在不脱离本发明的精神和范围的情况下,能够对在这里列出的主题实施各种改变、同等物的置换和修改,因此本发明的保护范围当视所提出的权利要求限定的范围为准。
Claims (4)
1.涌潮流速垂向分布的计算方法,适用于一电子装置,其特征在于:包括如下步骤:
步骤一:获取相关数据
通过现场观测或试验测量的方法,获取某一涌潮的潮前水深h、潮前流速v、涌潮高度H、离河底距离为Z1处的涌潮流速V1;
步骤二:判定涌潮类型
通过涌潮Froude数判定该涌潮是波状涌潮还是漩滚涌潮,具体方法如下:
a、计算涌潮的Froude数,公式如下:
式(1)和式(2)中,U为涌潮传播速度,v为潮前流速,h为潮前水深,H为涌潮高度;
b、计算得到涌潮Froude数Fr,如果Fr<1.7,则判断该涌潮为波状涌潮,如果Fr≥1.7,则判断该涌潮为漩滚涌潮;
步骤三:涌潮类型为波状涌潮的流速垂向分布,具体方法如下:
a、根据已知的离河底距离为Z1处的涌潮流速V1计算垂向最大流速Vmax,采用的是式(3),
b、根据公式(3)得到垂向最大流速Vmax,计算离河底任一距离Z处的涌潮流速V,采用的是式(4),
式(4)中,Z≤h+H,
步骤四:涌潮类型为漩滚涌潮的流速垂向分布,具体方法如下:
a、由已知的离河底距离Z1,计算该处的相对位置Z'1,
b、根据计算得到相对位置Z'1及涌潮流速V1,计算垂向最大流速Vmax,若Z'1≤0.6,则采用的是式(6),若Z'1>0.6,则采用式(7),
c、根据式(6)或式(7)得到垂向最大流速Vmax,计算离河底任一距离Z处的涌潮流速V,采用的是式(8),
式(8)中,Z≤h+H。
2.涌潮流速垂向分布的计算组件,其特征在于:包括,
存储模块,用于存储实现如权利要求1所述的方法的程序和/或参与运算的数据;
运算模块,用于实现如权利要求1所述的方法的运算过程。
3.如权利要求2所述的涌潮流速垂向分布的计算组件,其特征在于:还包括,
输入输出组件,用于将所述的参与运算的数据输入和结果的输出;
通讯模块,用于从远端接收所述的运算的数据和/或将结果传到远端。
4.一种存储介质,用以储存一计算机程序,其特征在于:所述的计算机程序包括多个程序代码,其用以加载至一电子装置,使得该电子装置执行如权利要求1所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910480043.7A CN110412312B (zh) | 2019-06-04 | 2019-06-04 | 涌潮流速垂向分布的一种计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910480043.7A CN110412312B (zh) | 2019-06-04 | 2019-06-04 | 涌潮流速垂向分布的一种计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110412312A CN110412312A (zh) | 2019-11-05 |
CN110412312B true CN110412312B (zh) | 2021-05-11 |
Family
ID=68358427
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910480043.7A Active CN110412312B (zh) | 2019-06-04 | 2019-06-04 | 涌潮流速垂向分布的一种计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110412312B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113344953B (zh) * | 2021-04-21 | 2023-07-04 | 中国计量大学 | 一种基于无人机的机器视觉涌潮流速测量方法 |
CN114218662B (zh) * | 2022-02-21 | 2022-06-07 | 浙江省水利河口研究院(浙江省海洋规划设计研究院) | 一种回头潮涌潮造景设计方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5515721A (en) * | 1993-08-25 | 1996-05-14 | Changmin Technology Co., Ltd. | Local river flow velocity measuring method and apparatus therefor |
CN202648661U (zh) * | 2012-06-21 | 2013-01-02 | 杭州市水文水资源监测总站 | 一种基于流速的江河涌潮检测装置 |
CN103090855A (zh) * | 2013-01-17 | 2013-05-08 | 杭州电子科技大学 | 一种基于水流流速判断涌潮是否到达的方法 |
CN103424106A (zh) * | 2013-08-27 | 2013-12-04 | 杭州市水文水资源监测总站 | 一种基于水位和流速的江河涌潮自动检测及实时预报方法 |
CN103981833A (zh) * | 2014-05-27 | 2014-08-13 | 浙江省水利河口研究院 | 高含沙量涌潮模拟试验方法 |
CN107632168A (zh) * | 2017-09-27 | 2018-01-26 | 杭州开闳环境科技有限公司 | 一种河道二维流速及流量测量系统及方法 |
CN109060056A (zh) * | 2018-08-20 | 2018-12-21 | 长江水利委员会长江科学院 | 一种非接触式雷达测流的河道断面流量计算方法 |
CN109582996A (zh) * | 2018-08-19 | 2019-04-05 | 珠江水利委员会珠江水利科学研究院 | 一种小尺度岸滩剖面与大尺度岸线变化的耦合模拟方法 |
CN109706881A (zh) * | 2018-12-11 | 2019-05-03 | 水利部交通运输部国家能源局南京水利科学研究院 | 河口海岸工程风暴潮灾害全动力实验室模拟系统及其方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201615909U (zh) * | 2010-02-02 | 2010-10-27 | 杭州电子科技大学 | 一种单杆潮汐流速测量装置 |
-
2019
- 2019-06-04 CN CN201910480043.7A patent/CN110412312B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5515721A (en) * | 1993-08-25 | 1996-05-14 | Changmin Technology Co., Ltd. | Local river flow velocity measuring method and apparatus therefor |
CN202648661U (zh) * | 2012-06-21 | 2013-01-02 | 杭州市水文水资源监测总站 | 一种基于流速的江河涌潮检测装置 |
CN103090855A (zh) * | 2013-01-17 | 2013-05-08 | 杭州电子科技大学 | 一种基于水流流速判断涌潮是否到达的方法 |
CN103424106A (zh) * | 2013-08-27 | 2013-12-04 | 杭州市水文水资源监测总站 | 一种基于水位和流速的江河涌潮自动检测及实时预报方法 |
CN103981833A (zh) * | 2014-05-27 | 2014-08-13 | 浙江省水利河口研究院 | 高含沙量涌潮模拟试验方法 |
CN107632168A (zh) * | 2017-09-27 | 2018-01-26 | 杭州开闳环境科技有限公司 | 一种河道二维流速及流量测量系统及方法 |
CN109582996A (zh) * | 2018-08-19 | 2019-04-05 | 珠江水利委员会珠江水利科学研究院 | 一种小尺度岸滩剖面与大尺度岸线变化的耦合模拟方法 |
CN109060056A (zh) * | 2018-08-20 | 2018-12-21 | 长江水利委员会长江科学院 | 一种非接触式雷达测流的河道断面流量计算方法 |
CN109706881A (zh) * | 2018-12-11 | 2019-05-03 | 水利部交通运输部国家能源局南京水利科学研究院 | 河口海岸工程风暴潮灾害全动力实验室模拟系统及其方法 |
Non-Patent Citations (5)
Title |
---|
Experimental hydrodynamic study of the Qiantang River tidal bore;HUANG Jing et al.;《Journal of Hydrodynamics》;20130330;第25卷(第3期);第481-490页 * |
基于定点连续观测的钱塘江涌潮特性研究;张巍等;《水动力学研究与进展》;20170330;第32卷(第2期);第253-259页 * |
基于实测资料的钱塘江涌潮水动力学特性研究;谢东风等;《水动力学研究与进展》;20120930;第27卷(第5期);第501-508页 * |
涌潮水力学特性试验研究;杨火其等;《水电能源科学》;20080830;第26卷(第4期);第136-138页 * |
钱塘江涌潮流速研究;潘存鸿;《海洋工程》;20170930;第35卷(第5期);第33-41页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110412312A (zh) | 2019-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110412312B (zh) | 涌潮流速垂向分布的一种计算方法 | |
Robertson et al. | OC5 Project Phase Ib: Validation of hydrodynamic loading on a fixed, flexible cylinder for offshore wind applications | |
CN109932711B (zh) | 一种雷达测量的大气折射修正方法 | |
CN111008465A (zh) | 考虑粘聚力与内摩擦角贡献度的边坡安全系数计算方法 | |
CN112085270A (zh) | 基于随机统计模型与水动力模型耦合的风暴潮极端水位预测方法 | |
US10579756B2 (en) | Simulation method of surface water flow movement process in surface irrigation | |
CN112818512A (zh) | 一种基于风暴模式的最大单波波高Hmax计算方法 | |
CN102426352B (zh) | 风廓线雷达风计算方法 | |
CN104200110A (zh) | 基于去耦合的二维波束扫描测角方法 | |
Kim et al. | Nonlinear time-domain NWT simulations for two types of a backward bent duct buoy (BBDB) compared with 2D wave-tank experiments | |
CN110388899B (zh) | 用于涌潮流速垂向分布计算的专用组件 | |
Lou et al. | Multiobjective equivalent static wind loads on complex tall buildings using non-gaussian peak factors | |
CN210293214U (zh) | 一种越浪水体上挑角度测量装置 | |
CN110398235B (zh) | 一种越浪水体上挑角度测量装置、方法及应用 | |
CN105160156A (zh) | 一种融合数据预处理的垂测电离图反演方法 | |
CN114166399B (zh) | 一种用于反演坝面水下爆炸压力的方法 | |
CN106092415A (zh) | 多坡度、多角度海浪冲击力现场测量装置及使用方法 | |
CN114045886B (zh) | 一种gil管廊基础沉降在线监测方法 | |
CN110057488B (zh) | 涌潮冲击圆柱桩的最大压强及其垂向分布的计算方法 | |
CN111487616B (zh) | 一种流量计量方法 | |
CN110093895A (zh) | 用于圆柱型构筑物涌潮作用力计算的专用组件 | |
Abualtayef et al. | Three-dimensional hydrostatic modeling of a bay coastal area | |
CN110134993B (zh) | 涌潮作用于圆柱桩的最大压强及其垂向分布计算的专用组件 | |
CN205981518U (zh) | 多坡度、多角度海浪冲击力现场测量装置 | |
ÖZMEN et al. | EXPERIMENTAL INVESTIGATION OF WIND PRESSURE FIELDS ON BUILDINGS WITH GABLED ROOFS HAVING DIFFERENT PITCH ANGLES. |
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 |