一种基于神经网络的台风浪对防波堤破坏预警方法及系统与流程

1.本技术属于港口与海岸工程领域，具体涉及一种基于神经网络的台风浪对防波堤破坏预警方法及系统。

背景技术：

2.防波堤是港口、海岸工程建设中的重要水工建筑物之一，其作用是保护港口的水域不受波浪的侵袭或减小波浪的作用。由于台风引起的极端海况形成冲击能量，往往使得防波堤遭受严重的破坏，破坏的主因为强风引起的波浪和水位叠加，在破坏中的防波堤，91%以上是护面块体或者坡脚块石等散粒体滚落失稳造成的。另外破坏中的防波堤大部分都建设在水深大浪处，因此给破坏后的防波堤修复难度加大，以及散粒体修复费用成倍增加。因此面临台风浪极端海况破坏力强，破坏造成经济损失巨大难题，提供一种精确、实用的防波堤破坏智能预警技术具有极大的需求，也是工程界面临的一项重要而迫切的任务。
3.目前揭示在波浪作用下防波堤散粒体失稳破坏的技术手段，常用的有经验公式计算、物理模型试验和数值模拟3种，由于经验公式只能用来近似的估计，通用性不强；物理模型的研究，受测量仪器和试验场地、缩尺效应和模型制作时间等限制，存在模拟范围小、时间滞后；数值模拟较物理模型试验具有方便、快捷和经济等优点，但散粒体失稳受各种复杂因素的影响表现为非线性动力学过程，且因素之间的变化及相互影响关系也很难确定，目前仍无法建立起在各种情况下普遍采用的数学模型，同时上述3种研究技术手段，均无法快速反映未来时间段台风登陆时，实现快速预报和预警的功能，难以为后方港口管理者开展应急保障提供足够的时长。
4.基于上述急需问题和目标实现，鉴于人工神经网络具有很强大的处理复杂非线性问题，以及快速准确预报预警的能力，因此发明了一种基于神经网络的台风浪对防波堤破坏预警方法及系统。利用人工神经网络自学习功能，采用编制新程序的算法，以模型试验成果作为神经网络的训练和测试样本，将测模型试验结果、经验公式计算和现场实测数据进行检验，经多次迭代方式，最终达到期望目标差；解决了实际问题，实现了快速预警的目标。

技术实现要素：

5.本技术提出了一种基于神经网络的台风浪对防波堤破坏预警方法及系统。采用由一个输入层、一个隐含层和一个输出层构成的三层改进的bp（back-propagation）神经网络，网络训练方法采用l-m（非线性最小二乘法）算法，神经网络输出值与目标值之间的吻合程度用均方根误差进行衡量，隐含层和输出层的激活函数为分别为sigmoid型函数和线性函数。网络运行流程为：首先确定台风浪极端海况下防波堤散粒体破坏的主要影响因素，选取足够多的训练样本数据；其次设置参数：包括输入层及输出层神经元、隐含层的数目、神经元数、初始权值和激活函数、初始学习速率、误差速率和期望目标误差等；最后开始训练网络，获得计算结果后，将网络模拟结果与模型试验结果、经验公式计算进行对比检验，若两者相关系数较差，则进一步对预报模型设置参数优化后再次迭代计算，直到满足期望目
标误差允许范围。
6.为实现上述目的，本技术提供了如下方案：一种基于神经网络的台风浪对防波堤破坏预警方法，包括如下步骤：s1.获取防波堤位置波浪要素、水位组合条件和防波堤横断面参数；s2.基于所述防波堤位置波浪要素、所述水位组合条件和所述防波堤横断面参数，开展水动力环境-防波堤散粒体受力分析，获得散粒体稳定性影响的关键因素；其中，所述散粒体受力分析方法包括：对防波堤横断面斜坡上散粒体受力分析，得到散粒体受波浪作用下的水流浮托力，以及水下重量，也即散粒体重量与水对散粒体产生浮力之差；所述水流浮托力包括速度力fd和惯性力fm，
7.其中，分别为速度力和惯性力系数；分别为面积和体积系数；l为散粒体长度，g为重力加速度，u为波浪产生的水流速度，为水的重度；将惯性力fm的影响合并入速度力fd中，简化后水流浮托力公式为：
8.其中，为联合影响系数，是与和有关的函数；波浪破碎时计算式为：
9.其中，为波浪破碎系数，与防波堤断面坡度m、水深d、波高h有关，则
10.所述防波堤斜坡上散粒体水下重量计算式为：
11.其中，为材料重度，为水的重度；根据散粒体在防波堤斜坡上受力分析，得到极限平衡条件下的力的平衡方程为：
12.其中，为防波堤坡度夹角；散粒体重量，即：
13.将式（9）代入式（8）中，因此力的平衡方程式可变换为：
14.其中，上述公式右侧项为一综合系数，与防波堤断面坡度m、波浪参数、破波水流流态和散粒体特性有关；将上述力的平衡方程式，再次进行变换，最终得到波浪作用下，与波高、水深和防波堤尺度、材料性质参数有关的散粒体重量计算公式为：
15.其中，为散粒体重量，为波浪陡度，表示波高/波长，h表示波高，l表示波长；为相对水深，表示水深/波长，为水深，为防波堤断面坡度，为散粒体的稳定系数，为波浪传播过程中的波浪方向；由上述波浪作用下散粒体稳定重量关系式得到：防波堤断面坡度、波浪陡度、相对水深、稳定系数、材料重度、稳定系数和波浪方向为关键因素；s3.基于所述散粒体稳定性影响的关键因素，建立水动力环境—散粒体失稳重量数据库；s4.基于改进后的bp网络，并采用lm算法，结合水动力环境—散粒体失稳重量数据库，提取训练样本，并设置网络参数，利用针对预警预报精度的误差分析方法，获取隐含层数目，隐含层节点数；基于所述散粒体稳定性影响的关键因素，确定输入层、输出层的神经元类型和数量，开展网络样本训练、样本测试和样本检验，获得台风浪极端海况对防波堤破坏预警模型。
16.优选的，所述预警方法还包括：s5.对所述预警模型进行优化，包括精度优化和计算速度优化；所述优化的方法包括：通过对比所述预警模型的计算结果与试验结果的相关性，建立波浪与散粒体重量响应机制，若所述波浪与所述散粒体重量之间相关系数低于预设值的，则进行误差收敛，进一步对所述预警模型设置参数进行优化，直到满足预定目标误差允许范围。
17.优选的，所述s1中，获取所述防波堤位置波浪要素和所述水位组合条件的方法包括：利用台风预报数据文件，以及造成破坏的波浪和水位因素，采用台风浪swan模型进行数据处理，获得防波堤位置波浪、水位要素，并波浪、水位参数的模拟区间的范围。
18.优选的，所述s1中，获取所述防波堤横断面参数的方法包括：将整个防波堤结构划分为堤头、堤身和堤根三段，获得每段断面参数和散粒体比重相对值，所述断面参数包括坡度和散粒体护面的稳定系数。
19.优选的，所述s5中，计算速度优化的方法包括：基于改进的bp网络，采用附加动量法和自适应学习速率法，进行计算速度优化。
20.另一方面，为实现上述目的，本技术还提出了一种基于神经网络的台风浪对防波堤破坏预警系统，包括：参数获取单元、影响因素单元、数据库单元和预警模型单元；所述参数获取单元用于获取防波堤位置波浪要素、水位组合条件和防波堤横断面参数；所述影响因素单元用于基于所述防波堤位置波浪要素、所述水位组合条件和所述防波堤横断面参数，开展水动力环境-防波堤散粒体受力分析，获得散粒体稳定性影响的关键因素；其中，所述散粒体受力分析方法包括：对防波堤横断面斜坡上散粒体受力分析，得到散粒体受波浪作用下的水流浮托力，以及水下重量，也即散粒体重量与水对散粒体产生浮力之差；所述水流浮托力包括速度力fd和惯性力fm，
21.其中，分别为速度力和惯性力系数；分别为面积和体积系数；l为散粒体长度，g为重力加速度，u为波浪产生的水流速度，为水的重度；将惯性力fm的影响合并入速度力fd中，简化后水流浮托力公式为：
22.其中，为联合影响系数，是与和有关的函数；波浪破碎时计算式为：
23.其中，为波浪破碎系数，与防波堤断面坡度m、水深d、波高h有关，则
24.所述防波堤斜坡上散粒体的水下重量计算式为：
25.其中，为材料重度，为水的重度；根据散粒体在防波堤斜坡上受力分析，得到极限平衡条件下的力的平衡方程为：
26.其中，为防波堤坡度夹角；散粒体重量，即：
27.将式（9）代入式（8）中，因此力的平衡方程式可变换为：
28.其中，上述公式右侧项为一综合系数，与防波堤断面坡度m、波浪参数、破波水流流态和散粒体特性有关；将上述力的平衡方程式，再次进行变换，最终得到波浪作用下，与波高、水深和防波堤尺度、材料性质等参数有关的散粒体重量计算公式为：
29.其中，为波浪陡度，表示波高/波长，h表示波高，l表示波长；为相对水深，表示水深/波长，为水深，为防波堤断面坡度，为散粒体的稳定系数，为波浪传播过程中的波浪方向；由上述波浪作用下散粒体稳定重量关系式得到：防波堤断面坡度、波浪陡度、相对水深、稳定系数、材料重度、稳定系数和波浪方向为关键因素；所述数据库单元用于基于所述散粒体稳定性影响的关键因素，建立水动力环境—散粒体失稳重量数据库；所述预警模型单元用于基于改进后的bp网络，并采用lm算法，结合水动力环境—散粒体失稳重量数据库，提取训练样本，并设置网络参数，利用针对预警预报精度的误差分
析方法，获取隐含层数目，隐含层节点数；基于所述散粒体稳定性影响的关键因素，确定输入层、输出层的神经元类型和数量，开展网络样本训练、样本测试和样本检验，获得台风浪极端海况对防波堤破坏预警模型，所述预警模型在台风浪极端海况时对防波堤破坏做出预警。
30.优选的，所述预警系统还包括优化单元；所述优化单元用于对所述预警模型单元进行优化，包括精度优化和计算速度优化。
31.本技术的有益效果为：本技术公开了一种基于神经网络的台风浪对防波堤破坏预警方法及系统，通过开展系列组合模型试验，获得水动力环境—散粒体失稳破坏重量数据库，运用改进的bp网络模型，经网络训练和检验，可快速对台风浪极端海况的防波堤破坏进行预警，模拟结果一方面可为港口设计者在防波堤前期提供基础数据，另一方面为港口管理者在防波堤运营期提供更长的应急准备时间，大大降低了因台风浪导致防波堤破坏后对掩护区内设施次级灾害。
附图说明
32.为了更清楚地说明本技术的技术方案，下面对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅是本技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
33.图1为本技术实施例一的预警方法流程示意图；图2为本技术实施例一的开展系列模型试验工况组合示意图；图3为本技术实施例一的防波堤断面斜坡上散粒体受力分析示意图；图4为本技术实施例一的神经网络参数设置示意图；图5为本技术实施例一的神经网络模拟的流程示意图；图6为本技术实施例一的网络训练结果与模型结果相关性对比示意图；图7为本技术实施例一的模拟预警结果与模型试验结果、公式结果数据和设计值对比示意图；图8为本技术实施例二的整体系统结构示意图。
具体实施方式
34.下面将结合本技术实施例中的附图，对本技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本技术一部分实施例，而不是全部的实施例。基于本技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本技术保护的范围。
35.为使本技术的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本技术作进一步详细的说明。
36.实施例一：如图1所示，为本技术实施例的预警方法流程示意图；包括以下步骤：s1.获取防波堤位置波浪要素、水位组合条件和防波堤横断面参数；利用国家海洋预报台风网公开台风预报数据文件，以及造成破坏的波浪和水位等主要因素，采用成熟的
台风浪swan模型进行计算再分析，获得了防波堤位置波浪、水位要素等水动力条件，确定波浪、水位参数的模拟区间的范围，其中包括波高h、周期t、波浪方向α和一场台风浪持续作用的时间t,以及极端高水位h1、设计高水位h2、平均水位h3、极端低水位h4等。图2所示，具体包括：
①
、基于防波堤轴线沿程波浪要素分布统计，确定n1个需要计算的代表波浪方向=0
°
~315
°
，间隔为45
°
，共计8个组次；确定n2个需要计算的代表波高h=4.0~16.0m，间隔2.0m，共计7个组次；确定n3个需要计算的代表周期t=8.0~20.0s，间隔2.0s，共计7个组次；
②
、基于防波堤位置水位分析统计，确定n4个需要计算的代表水位h=2.0~10.0m，间隔为2.0m，共计5个组次；
③
、由
①
和
②
计算条件叠加，确定所需模拟的总组次为：n1
×
n2
×
n3
×
n4=8
×7×7×5×
5=9800，以下类似。
37.确定防波堤断面参数：受防波堤轴沿程水深不同，则波浪对散粒体冲击作用的程度和失稳形式不同，因此将整个防波堤按不同结构划分为堤头、堤身和接岸的堤根3段，从而分别获得每段横断面参数。受防波堤轴沿程水深不同，基于《防波堤与护岸设计规范》（jts154-2018）规定，因此将整个防波堤结构划分为堤头、堤身和堤根三段，从而获得每段断面参数：防波堤断面坡度m，稳定系数kd和散粒体比重相对值为材料重度，为水的重度，具体包括：
①
、基于《防波堤与护岸设计规范》（jts154-2018）中规定，给出所有结构断面型式坡度选择，确定n5个需要计算坡度m=1:1.25、1:1.5、1:2、1:3，共计4个组次；
②
、基于《防波堤与护岸设计规范》（jts154-2018）中规定，给出工程上常用的散粒体护面包括块石、四脚锥体、四脚空心方块、扭工字和扭王字块稳定系数kd和比重相对值选择，确定n6个需要计算稳定系kd=4.0、8.0、14.0、16.0，共计4个组次；确定n7个需要计算比重相对值=2.65、2.3，共计2个组次；
③
、由
①
和
②
计算条件叠加，确定所需模拟的总组次为：n5
×
n6
×
n7=4
×4×
2=32。
38.s2.基于防波堤位置波浪要素、水位组合条件和防波堤横断面参数，开展水动力环境-防波堤散粒体受力分析，获得散粒体稳定性影响的关键因素。
39.其中，对防波堤断面斜坡上散粒体受力分析的方法，如图3所示，散粒体主要受波浪作用下的水流浮托力，以及水下重量，也即散粒体重量与水对散粒体产生浮力之差。当波浪冲击作用下，坡面上散粒体失稳的瞬间，在忽略此时散粒体之间的摩擦力，利用力的平衡原理，获得了平衡力的等量方程式为：。
40.其中，水流浮托力包括两部分，分别为：
速度力：
41.惯性力：
42.其中，分别为速度力和惯性力系数；分别为面积体积系数；为散粒体长度，g为重力加速度，u为水流速度，为水的重度。
43.由于实际应用过程中和水流加速度的取值存在一定困难，因此将惯性力的影响合并入速度力fd中，j简化后水流浮托力公式为：
44.其中，为联合影响系数，是与和有关的函数；其它字母含义同上。
45.在波浪破碎时计算式为：
46.其中，为波浪破碎系数，与防波堤断面坡度m，也即垂直高度与水平距离的比值，水深d，也即工程位置处海底高程与水位之间的和，和波高h有关；因此水流浮托力可变换为：
47.散粒体水下重量计算式为：
48.其中，为材料重度，为水的重度；根据散粒体在防波堤斜坡上受力分析，得到极限平衡条件下的力的平衡方程为：
49.其中，为防波堤坡度夹角；散粒体重量，也即可以写成：
50.将式（9）代入式（8）中，因此力的平衡方程式可继续变换为：
51.其中，上述公式右侧项为一综合系数，与防波堤断面坡度m、波浪参数、破波水流流态和散粒体特性有关；将上述力的平衡方程式，再次进行变换，最终得到波浪作用下，与波高、水深和防波堤尺度、材料性质等参数有关的散粒体重量计算公式为：
52.其中，为散粒体重量，为波浪陡度，表示波高/波长，h表示波高，l表示波长；为相对水深，表示水深/波长，为水深，也即工程位置处海底高程与水位之间的和，为防波堤断面坡度，也即垂直高度与水平距离的比值，为散粒体的稳定系数，为波浪传播过程中的波浪方向；由式（11）得到：防波堤断面坡度、波浪陡度、相对水深、稳定系数、材料重度、稳定系数和波浪方向为关键因素。
53.s3.基于所述散粒体稳定性影响的关键因素，建立水动力环境—散粒体失稳重量数据库；基于s1中系列工况组合，开展系列模型（物理模型、数学模型）试验，具体如下：以防波堤堤头段为例进一步说明，首先查询并确定防波堤堤头段断面参数信息，选择适当的模型比尺或数学模型计算范围，针对防波堤堤头段断面形式和堤头区域的水下地形，开展试验模型制作或数值化建模；其次设置不同波浪(波高h、周期t、波浪方向α、波浪持续作用时间t)和水位（极端高水位h1、设计高水位h2、设计低水位h3、极端低水位h4）组合工况，进一步开展正式试验或计算；最后对研究成果进行分析处理，获得系列工况组合条件下，水动力环境—散粒体失稳重量的数据文件；同样针对防波堤其它段（堤身和堤根），利用上述同样试验方法，获得重量数据文件，一并存入数据库中。因此模型试验模拟的总组次：n=n1
×
n2
×
n3
×
n4
×
n5
×
n6
×
n7。
54.s4.基于改进后的bp网络，并采用lm算法，结合所述水动力环境—散粒体失稳重量数据库，构建并训练改进后的bp神经学习网络，得到台风浪极端海况对防波堤破坏预警模型，所述预警模型用于在台风浪极端海况时对防波堤破坏做出预警。
55.基于改进后的bp网络，采用lm算法，基于s2获得的水动力环境—散粒体失稳重量数据库，提取部分为训练样本，设置网络参数，利用针对预警预报精度的误差分析方法，获取隐含层数目，隐含层节点数，基于s3散粒体失稳重量因素，确定输入层、输出层的神经元类型和数量，开展样本网络训练、样本测试和样本检验，形成本技术基于人工神经网络的台
风浪极端海况对防波堤破坏预警方法。如图4、图5所示，具体包括：
①
将数据库分成两部分样本，训练样本和测试样本。各输入神经元信息依据如下函数表达式在神经网络中传递：
56.其中上标表示第层神经网络，表示所在层内第个神经元，为激活函数，为层上第个神经元与层上第个神经元之间的权重系数，为第层的偏置系数。
57.②
设置网络参数，利用针对预警预报精度的误差分析，获取隐含层数目为1个，隐含层节点数，依据经验公式，其中，为隐含层输入神经元数，为隐含层输入单元数，为[1,10]之间的常数，利用模型误差期望分析法，最终得到最优节点数为10个，训练过程见表1。
[0058]
表1
[0059]
③
初始权值及激活函数的选取：由于系统是非线性的，初始权值对于学习是否达到局部最小，是否能够收敛以及训练时间长短的关系很大，一般取初始权值在(0，1)之间的随机值。隐含层和输出层的激活函数为分别为型函数和线性函数。
[0060]
在本实施例中，预警方法还包括：s5.对所述预警模型进行优化，包括精度优化和计算速度优化；即通过对比计算结果与试验结果的相关性，建立波浪与散粒体重量响应机制，若两者相关系数低于预设值的，则认为两者相关系数较差，编写误差收敛程序，进一步对预警模型设置参数进行优化，直到满足期望目标误差允许范围。
[0061]
①
进一步加快训练速度方法，基于改进型bp网络，采用附加动量法和自适应学习速率法，计算表达式分别如下：附加动量因子的权值调整计算公式：
[0062]
其中，为权值变化量；为偏差变化量；为误差项；为输入矢量；
为学习速率；为动量因子，一般取0.95。
[0063]
自适应学习速率法的调整计算公式：
[0064]
其中，为误差平方和，初始学习速率的选取范围可以有很大的随意性。
[0065]
其它参数设置，初始学习速率=0.05，自适应学习速率递增乘因子，自适应学习速率递减乘因子；误差速率，为动量因子取为0.95；期望误差。
[0066]
②
模型误差评估和优化，进一步评估进行训练网络模型结果，对比数据随机在数据库选择部分数据，也可以用其它完整试验结果进行验证，按照预设期望误差的要求，由图6可以得到，神经网络计算结果与训练样本结果对比发现，结果表明拟合程度很好，基本都达到0.95以上，均分布在y＝x（1
±
10％）包络图内。
[0067]
检验网络模型预警效果，利用受台风浪作用下，防波堤发生破坏例子，以及所开展模型试验结果，进一步对网络模型预警效果进行检验。具体如下：利用以沿海某电厂防波堤在“山竹”台风浪作用下破坏为例子，在2018年受“山竹”台风浪的作用下，防波堤在堤头衔接处的位置遭到严重的破坏，为便于检验，选择台风浪近24个小时过程模拟（2018-9-160:00至2018-9-170:00），每个小时对应的堤头位置处的波高h以及同时刻的水位h值如表2所示,复演台风浪过境时防波堤破坏过程。
[0068]
防波堤堤头断面参数：散粒体为扭王字块，坡度m=1:2.0,稳定系数kd=16.0，散粒体相当比重=2.3，散粒体设计重量为8t。神经网络模型模拟结果、经验公式计算、物理模型试验和防波堤设计值对比见图7。
[0069]
由图7可知，神经网络模型模拟结果能很好的反映波浪强度与散粒体重量w变化关系；与物理模型试验结果和经验公式计算值对比，也具有很好的相关性，变化趋势基本一致；在防波堤设计的护面扭王字块散粒体为8t时，当台风浪波高达到5.5m时，则坡面上散粒体受波浪力大于其重量，最终导致扭王字块破坏失稳，网络模拟结果与现场监测结果基本吻合；对比图中扭王字块重量w随波高h变化的4条可知，神经网络模型模拟的重量趋势线，最先与设计扭王字块8t重线相交，说明本次所发明的神经网络台风浪极端海况预警方法，更能快反映防波堤散粒体发生破坏的风险。
[0070]
表2
[0071]
实施例二如图8所示，为本技术实施例的整体系统结构示意图；包括：参数获取单元、影响因素单元、数据库单元和预警模型单元；具体的，参数获取单元用于获取防波堤位置波浪要素、水位组合条件和防波堤横断面参数；影响因素单元用于基于防波堤位置波浪要素、水位组合条件和防波堤横断面参数，开展水动力环境-防波堤散粒体受力分析，获得散粒体稳定性影响的关键因素；其中，所述散粒体受力分析方法包括：对防波堤横断面斜坡上散粒体受力分析，得到散粒体受波浪作用下的水流浮托力，以及水下重量，也即散粒体重量与水对散粒体产生浮力之差；所述水流浮托力包括速度力fd和惯性力fm，
[0072]
其中，分别为速度力和惯性力系数；分别为面积和体积系数；l为散粒体长度，g为重力加速度，u为波浪产生的水流速度，为水的重度；将惯性力fm的影响合并入速度力fd中，简化后水流浮托力公式为：
[0073]
其中，为联合影响系数，是与和有关的函数；波浪破碎时计算式为：
[0074]
其中，为波浪破碎系数，与防波堤断面坡度m，也即垂直高度与水平距离的比值，水深d，也即工程位置处海底高程与水位之间的和，波高h有关，则
[0075]
所述防波堤斜坡上散粒体的水下重量计算式为：
[0076]
其中，为材料重度，为水的重度；根据散粒体在防波堤斜坡上受力分析，得到极限平衡条件下的力的平衡方程为：
[0077]
其中，为防波堤坡度夹角；散粒体重量，也即可以写成：
[0078]
将式（9）代入式（8）中，因此力的平衡方程式可变换为：
[0079]
其中，上述公式右侧项为一综合系数，与防波堤断面坡度m、波浪参数、破波水流流态和散粒体特性有关；将上述力的平衡方程式，再次进行变换，最终得到波浪作用下，与波高、水深和防波堤尺度、材料性质等参数有关的散粒体重量计算公式为：
[0080]
其中，为散粒体重量，为波浪陡度，表示波高/波长，h表示波高，l表示波长；为相对水深，表示水深/波长，为水深，为防波堤断面坡度，为散粒体的稳定系数，为波浪传播过程中的波浪方向；由上述波浪作用下散粒体稳定重量关系式得到：防波堤断面坡度、波浪陡度、相对水深、稳定系数、材料重度、稳定系数和波浪方向为关键因素；数据库单元用于基于散粒体稳定性影响的关键因素，建立水动力环境—散粒体失稳重量数据库；预警模型单元用于基于改进后的bp网络，并采用lm算法，结合水动力环境—散粒体失稳重量数据库，提取训练样本，并设置网络参数，利用针对预警预报精度的误差分析方法，获取隐含层数目，隐含层节点数；基于所述散粒体稳定性影响的关键因素，确定输入层、输出层的神经元类型和数量，开展网络样本训练、样本测试和样本检验，获得台风浪极端海况对防波堤破坏预警模型，该预警模型在台风浪极端海况时对防波堤破坏做出预警。
[0081]
进一步的，在本实施例中，预警系统还包括优化单元；优化单元用于对预警模型单元进行优化，包括精度优化和计算速度优化。
[0082]
以上所述的实施例仅是对本技术优选方式进行的描述，并非对本技术的范围进行限定，在不脱离本技术设计精神的前提下，本领域普通技术人员对本技术的技术方案做出的各种变形和改进，均应落入本技术权利要求书确定的保护范围内。

技术特征：

1.一种基于神经网络的台风浪对防波堤破坏预警方法，其特征在于，包括如下步骤：s1.获取防波堤位置波浪要素、水位组合条件和防波堤横断面参数；s2.基于所述防波堤位置波浪要素、所述水位组合条件和所述防波堤横断面参数，开展水动力环境-防波堤散粒体受力分析，获得散粒体稳定性影响的关键因素；其中，所述散粒体受力分析方法包括：对防波堤横断面斜坡上散粒体受力分析，得到散粒体受波浪作用下的水流浮托力，以及水下重量，也即散粒体重量与水对散粒体产生浮力之差；所述水流浮托力包括速度力f
d
和惯性力f
m
，（2）（3）其中，分别为速度力和惯性力系数；分别为面积和体积系数；l为散粒体长度，g为重力加速度，u为波浪产生的水流速度，为水的重度；将惯性力f
m
的影响合并入速度力f
d
中，简化后水流浮托力公式为：（4）其中，为联合影响系数，是与和有关的函数；波浪破碎时计算式为：（5）其中，为波浪破碎系数，与防波堤断面坡度m、水深d、波高h有关，则（6）所述防波堤斜坡上散粒体水下重量计算式为：（7）其中，为材料重度，为水的重度；根据散粒体在防波堤斜坡上受力分析，得到极限平衡条件下的力的平衡方程为：（8）其中，为防波堤坡度夹角；
散粒体重量，即：（9）将式（9）代入式（8）中，因此力的平衡方程式可变换为：（10）其中，上述公式右侧项为一综合系数，与防波堤断面坡度m、波浪参数、破波水流流态和散粒体特性有关；将上述力的平衡方程式，再次进行变换，最终得到波浪作用下，与波高、水深和防波堤尺度、材料性质参数有关的散粒体重量计算公式为：（11）其中，为散粒体重量，为波浪陡度，表示波高/波长，h表示波高，l表示波长；为相对水深，表示水深/波长，为水深，为防波堤断面坡度，为散粒体的稳定系数，为波浪传播过程中的波浪方向；由上述波浪作用下散粒体稳定重量关系式得到：防波堤断面坡度、波浪陡度、相对水深、稳定系数、材料重度、稳定系数和波浪方向为关键因素；s3.基于所述散粒体稳定性影响的关键因素，建立水动力环境—散粒体失稳重量数据库；s4.基于改进后的bp网络，并采用lm算法，结合水动力环境—散粒体失稳重量数据库，提取训练样本，并设置网络参数，利用针对预警预报精度的误差分析方法，获取隐含层数目，隐含层节点数；基于所述散粒体稳定性影响的关键因素，确定输入层、输出层的神经元类型和数量，开展网络样本训练、样本测试和样本检验，获得台风浪极端海况对防波堤破坏预警模型。2.根据权利要求1所述的基于神经网络的台风浪对防波堤破坏预警方法，其特征在于，所述预警方法还包括：s5.对所述预警模型进行优化，包括精度优化和计算速度优化；所述优化的方法包括：通过对比所述预警模型的计算结果与试验结果的相关性，建立波浪与散粒体重量响应机制，若所述波浪与所述散粒体重量之间相关系数低于预设值的，则进行误差收敛，进一步对所述预警模型设置参数进行优化，直到满足预定目标误差允许范围。3.根据权利要求1所述的基于神经网络的台风浪对防波堤破坏预警方法，其特征在于，所述s1中，获取所述防波堤位置波浪要素和所述水位组合条件的方法包括：利用台风预报数据文件，以及造成破坏的波浪和水位因素，采用台风浪swan模型进行数据处理，获得防波堤位置波浪、水位要素，并波浪、水位参数的模拟区间的范围。
4.根据权利要求1所述的基于神经网络的台风浪对防波堤破坏预警方法，其特征在于，所述s1中，获取所述防波堤横断面参数的方法包括：将整个防波堤结构划分为堤头、堤身和堤根三段，获得每段断面参数和散粒体比重相对值，所述断面参数包括坡度和散粒体护面的稳定系数。5.根据权利要求2所述的基于神经网络的台风浪对防波堤破坏预警方法，其特征在于，所述s5中，计算速度优化的方法包括：基于改进的bp网络，采用附加动量法和自适应学习速率法，进行计算速度优化。6.一种基于神经网络的台风浪对防波堤破坏预警系统，其特征在于，包括：参数获取单元、影响因素单元、数据库单元和预警模型单元；所述参数获取单元用于获取防波堤位置波浪要素、水位组合条件和防波堤横断面参数；所述影响因素单元用于基于所述防波堤位置波浪要素、所述水位组合条件和所述防波堤横断面参数，开展水动力环境-防波堤散粒体受力分析，获得散粒体稳定性影响的关键因素；其中，所述散粒体受力分析方法包括：对防波堤横断面斜坡上散粒体受力分析，得到散粒体受波浪作用下的水流浮托力，以及水下重量，也即散粒体重量与水对散粒体产生浮力之差；所述水流浮托力包括速度力f
d
和惯性力f
m
，（2）（3）其中，分别为速度力和惯性力系数；分别为面积和体积系数；l为散粒体长度，g为重力加速度，u为波浪产生的水流速度，为水的重度；将惯性力f
m
的影响合并入速度力f
d
中，简化后水流浮托力公式为：（4）其中，为联合影响系数，是与和有关的函数；波浪破碎时计算式为：（5）其中，为波浪破碎系数，与防波堤断面坡度m、水深d、波高h有关，则
（6）所述防波堤斜坡上散粒体的水下重量计算式为：（7）其中，为材料重度，为水的重度；根据散粒体在防波堤斜坡上受力分析，得到极限平衡条件下的力的平衡方程为：（8）其中，为防波堤坡度夹角；散粒体重量，即：（9）将式（9）代入式（8）中，因此力的平衡方程式可变换为：（10）其中，上述公式右侧项为一综合系数，与防波堤断面坡度m、波浪参数、破波水流流态和散粒体特性有关；将上述力的平衡方程式，再次进行变换，最终得到波浪作用下，与波高、水深和防波堤尺度、材料性质等参数有关的散粒体重量计算公式为：（11）其中，为波浪陡度，表示波高/波长，h表示波高，l表示波长；为相对水深，表示水深/波长，为水深，为防波堤断面坡度，为散粒体的稳定系数，为波浪传播过程中的波浪方向；由上述波浪作用下散粒体稳定重量关系式得到：防波堤断面坡度、波浪陡度、相对水深、稳定系数、材料重度、稳定系数和波浪方向为关键因素；所述数据库单元用于基于所述散粒体稳定性影响的关键因素，建立水动力环境—散粒体失稳重量数据库；所述预警模型单元用于基于改进后的bp网络，并采用lm算法，结合水动力环境—散粒体失稳重量数据库，提取训练样本，并设置网络参数，利用针对预警预报精度的误差分析方法，获取隐含层数目，隐含层节点数；基于所述散粒体稳定性影响的关键因素，确定输入层、输出层的神经元类型和数量，开展网络样本训练、样本测试和样本检验，获得台风浪极端海况对防波堤破坏预警模型，所述预警模型在台风浪极端海况时对防波堤破坏做出预警。
7.根据权利要求6所述的基于神经网络的台风浪对防波堤破坏预警系统，其特征在于，所述预警系统还包括优化单元；所述优化单元用于对所述预警模型单元进行优化，包括精度优化和计算速度优化。

技术总结

本申请公开一种基于神经网络的台风浪对防波堤破坏预警方法及系统，包括以下步骤：获取防波堤位置波浪要素、水位组合条件和防波堤横断面参数；基于防波堤位置波浪要素、水位组合条件和防波堤横断面参数，开展水动力环境-防波堤散粒体受力分析，获得散粒体稳定性影响的关键因素；基于散粒体稳定性影响的关键因素，建立水动力环境—散粒体失稳重量数据库；构建并训练改进后的BP神经学习网络，得到台风浪极端海况对防波堤破坏预警模型，在台风浪极端海况时对防波堤破坏做出预警。本申请开展系列组合模型试验，获得水动力环境—散粒体失稳破坏重量数据库，运用改进后的BP网络模型，经网络训练和检验，可对台风浪极端海况的防波堤破坏进行预警。破坏进行预警。破坏进行预警。