基于水沙动力自适应特性的河口治理模拟方法与流程

1.本发明涉及河口治理分析决策技术领域，尤其涉及基于水沙动力自适应特性的河口治理模拟方法。


背景技术：

2.随着国际经贸的兴盛，促使我国沿海经济带不断壮大，珠三角、长三角都依托内陆河及出海口而发展，河口的开发力度进一步加大。1)河岸两侧港口、码头、围填海等工程的建设，使得岸线不断向海(水域)推进，导致近岸的浅滩水域被侵占，河床形态发生不可逆的改变；2)同时，由于河道通航要求，需要从港口、码头位置向出海口水域开挖航道，会改变河床的自然形态；3)为了提高上游防洪能力，开展河口治理，进行清淤，开挖主槽，扩大行洪断面，保证上游河道和网河区洪水安全宣泄。
3.如何开展河口治理，采用什么样的开挖方案，是设计的重点。常规的研究方式：设计一条线路，开挖深槽，定期进行地形监测，如果淤积严重，影响到通航、行洪，再开展清淤，维护主槽稳定。


技术实现要素：

4.本发明为解决上述技术问题，提出了一种基于水沙动力自适应特性的河口治理模拟方法，以解决至少一个上述技术问题。
5.一种基于水沙动力自适应特性的河口治理模拟方法，包括以下步骤：
6.步骤s1：获取河口基础数据，并根据河口基础数据构建河口物理模型或数学模型；
7.步骤s2：开展物理模型清水动床冲刷试验或数学模型清水动床冲刷模拟计算，获取冲刷后的河床冲刷形态数据；
8.步骤s3：根据河床冲刷后的形态数据生成河口治理方案；
9.步骤s4：根据河口治理方案构建河口治理模型，并根据河口治理模型进行物理模型浑水定床淤积试验或数学模型悬沙淤积模拟计算，根据主槽淤积特性，生成主槽淤积特征数据；
10.步骤s5：对主槽淤积特征数据进行稳定性计算，生成主槽稳定性参数，以供河口治理。
11.本实施例不进行人为设计路线，采用物理模型动床试验的方式，在径流、潮汐动力的作用下，基于水沙相互作用、自我演变调节为主导，分析河口水域在自然动力的作用下，哪个区域会发展出主槽，根据自然演变的规律，按照自然深槽的走向，设计河口治理主槽的分布位置，再按照冲深，设计河底的开挖高程。
12.在本说明书的一个实施例中，步骤s1包括以下步骤：
13.步骤s11：获取原始河口基础数据，其中原始河口基础数据包括原始潮流流向数据、原始涨潮潮流数据、原始落潮潮流数据以及原始潮流间拓扑关系数据；
14.步骤s12：获取原始河口图像信息序列以及修建堤防后获取的堤防河口图像信息
序列，并对原始河口图像信息序列以及堤防河口图像信息序列进行特征提取，生成原始河口图像特征运动数据集以及堤防河口图像特征运动数据集；
15.步骤s13：根据原始河口图像特征运动数据集以及堤防河口图像特征运动数据集通过潮流系数加权计算公式进行计算，生成潮流加权系数，并将原始河口基础数据通过潮流加权系数进行加权计算并进行率定，生成河口基础数据；
16.步骤s14：根据河口基础数据生成河口水动力特性数据；
17.步骤s15：根据河口水动力特性数据以及河口基础数据构建河口物理模型或数学模型。
18.本实施例通过图像数据分析，通过原始河口图像信息以及设立水利工程后的堤防河口图像信息进行比对计算，生成加权系数，以根据原始河口图像信息的记录参数进行加权计算，从而生成稳定可靠的河口基础数据，以建立准确的河口物理模型。
19.在本说明书的一个实施例中，其中潮流系数加权计算公式具体为：
[0020][0021]
p
t
为潮流加权系数，ai为原始河口图像特征运动数据集中的第i个原始河口图像特征运动数据，αi为原始河口图像特征运动数据集中的第i个原始河口图像特征运动数据相应的权重信息，bi为堤防河口图像特征运动数据集中的第i个堤防河口图像特征运动数据，βi为堤防河口图像特征运动数据集中的第i个堤防河口图像特征运动数据的参数信息的权重信息，ui为原始河口图像特征运动数据集中的第i个原始河口图像特征运动数据对应的流量速度信息，θ为调整项的权重信息，为调整项的误差系数，x为根据标准河口图像特征生成的调整项，i的取值为1、2、3...n，n为潮流加权系数的修正项。
[0022]
本实施例提供一种潮流系数加权计算公式，该公式充分考虑了原始河口图像特征运动数据集中的第i个原始河口图像特征运动数据ai、原始河口图像特征运动数据集中的第i个原始河口图像特征运动数据相应的权重信息αi、堤防河口图像特征运动数据集中的第i个堤防河口图像特征运动数据bi、堤防河口图像特征运动数据集中的第i个堤防河口图像特征运动数据的参数信息的权重信息βi、原始河口图像特征运动数据集中的第i个原始河口图像特征运动数据对应的流量速度信息ui、调整项的权重信息θ、调整项的误差系数根据标准河口图像特征生成的调整项x以及之间的相互关系，以形成函数关系并通过潮流加权系数的修正项n进行修正，从而提供准确可靠的数据支撑。
[0023]
在本说明书的一个实施例中，步骤s12包括以下步骤：
[0024]
步骤s121：获取原始河口图像信息序列以及堤防河口图像信息序列；
[0025]
步骤s122：对原始河口图像信息序列以及堤防河口图像信息序列进行关键像素点提取，生成原始河口图像信息序列的原始河口潮流分布特征点集序列以及堤防河口图像信息序列的堤防河口潮流分布特征点集序列；
[0026]
步骤s123：将原始河口图像信息序列中的当前原始河口图像信息的原始河口潮流分布特征点集与原始河口图像信息序列中的下一张原始河口图像信息的原始河口潮流分布特征点集进行匹配，生成原始河口潮流匹配特征点集，并将堤防河口图像信息序列中的
当前堤防河口图像信息的堤防河口潮流分布特征点集与堤防河口图像信息序列中的下一张堤防河口图像信息的堤防河口潮流分布特征点集进行匹配，生成堤防河口潮流匹配特征点集；
[0027]
步骤s124：根据步骤s123的方式依次将原始河口图像信息序列中的原始河口潮流分布特征点集与下一张原始河口图像信息的原始河口潮流分布特征点集进行匹配以及堤防河口图像信息序列中的堤防河口潮流分布特征点集与下一张堤防河口图像信息序列中的堤防河口潮流分布特征点集进行匹配，生成原始河口图像全匹配特征点集以及堤防河口图像全匹配特征点集，且根据原始河口图像全匹配集在原始河口图像信息中的相对位置以及堤防河口图像全匹配特征点集在堤防河口图像信息中的相对位置分布进行位置空间计算，生成原始河口图像特征运动数据集以及堤防河口图像特征运动数据集。
[0028]
本实施例通过对高清运动摄像头拍摄的图像信息序列进行依次比对，依次对相近两张的图像信息的特征进行比对，从而生成图像运动特征信息，并根据在图像的相对位置以及预设的参数分析方式，从而生成图像运动数据信息，提供准确的水动力分布特征信息。
[0029]
在本说明书的一个实施例中，步骤s2包括以下步骤：
[0030]
步骤s21：获取清水流量数据；
[0031]
步骤s22：基于河口物理模型或数学模型，通过清水动床试验的河口床沙计算公式对清水流量数据进行河口流量与河床推移质的运动关系的模拟计算，生成河口断面流量流速推移质运动数据；
[0032]
步骤s23：根据河口断面流量流速推移质运动数据对河口物理模型或数学模型进行修正，获取冲刷后的河床冲刷形态数据。
[0033]
本实施例通过清水动床试验计算，实现对断面形态变化的模拟，从而分析水沙自适应特性，研究深槽在自然动力作用下的发育、演变和分布。
[0034]
在本说明书的一个实施例中，其中河口床沙计算公式具体为：
[0035][0036]
q为河口断面流量流速推移质运动数据，v为清水流速信息，f为清水流量水动力含沙平均加权量，p为清水流速信息对应的平均含沙量，z为清水流量水动力的调整项，r为清水流速信息的平均含沙量，为清水流速信息的平均含沙量的实际误差调整项，为河口断面流量流速推移质运动数据的修正项。
[0037]
本实施例提供一种河口床沙计算公式，该公式充分考虑清水流速信息v、清水流量水动力含沙平均加权量f、清水流速信息对应的平均含沙量p、清水流量水动力的调整项z、清水流速信息的平均含沙量清水流速信息的平均含沙量的实际误差调整项r以及相互之间的作用关系，以形成函数关系并通过河口断面流量流速推移质运动数据的修正项进行修正，从而提供准确可靠的数据支撑。
[0038]
在本说明书的一个实施例中，步骤s3包括以下步骤：
[0039]
根据河床冲刷后的形态数据中的河床主槽数据以及河床深泓线分布数据，生成河口预治理方案；
[0040]
获取河道整治历史经验数据，并根据河道整治历史经验数据对河口预治理方案进
行调整，以生成河口治理方案。
[0041]
本实施例通过根据河床冲刷后的形态数据中的河床主槽数据以及河床深泓线分布数据，生成河口预治理方案，并获取河道整治历史经验数据，并根据河道整治历史经验数据对河口预治理方案进行调整，以生成河口治理方案，从而为下一步做好前提准备工作。
[0042]
在本说明书的一个实施例中，步骤s4包括以下步骤：
[0043]
根据河口治理方案构建河口治理模型；
[0044]
根据河口治理模型进行物理模型浑水定床淤积试验或数学模型悬沙淤积模拟计算，根据主槽淤积特性，生成主槽淤积特征数据，其中浑水定床试验计算的步骤为通过浑水淤积计算公式进行计算；
[0045]
其中浑水淤积计算公式具体为：
[0046][0047]
s为主槽淤积特性数据中的浑水淤积数据，s为浑水流速信息，γ为浑水流速水动力的当前含沙量的调整项，m为浑水流量水动力的调整项，h为浑水流量水动力含沙平均加权量，t为浑水流速信息的平均含沙量的实际误差调整项,c为浑水流速信息的平均含沙量，为浑水流速信息的平均含沙量的当前预估加权信息，为主槽淤积特性数据中的浑水淤积数据的修正项。
[0048]
本实施例通过根据河口治理方案构建河口治理模型，并根据河口治理模型进行物理模型浑水定床淤积试验或数学模型悬沙淤积模拟计算，根据主槽淤积特性，生成主槽淤积特征数据，从而为下一步做好前提准备。
[0049]
本实施例提供一种浑水淤积计算公式，该公式充分考虑浑水流速信息s、浑水流速水动力的当前含沙量的调整项γ、浑水流量水动力的调整项m、浑水流量水动力含沙平均加权量h、浑水流速信息的平均含沙量的实际误差调整项t、浑水流速信息的平均含沙量c、浑水流速信息的平均含沙量的当前预估加权信息以及相互之间的作用关系，以形成函数关系并通过主槽淤积特性数据中的浑水淤积数据的修正项进行修正，以提供准确的数据支持。
[0050]
在本说明书的一个实施例中，步骤s5包括以下步骤：
[0051]
通过主槽稳定性计算公式对主槽淤积特征数据进行计算，生成主槽稳定性参数；
[0052]
根据主槽稳定性参数生成主槽稳定性报告，以供河口治理；
[0053]
其中主槽稳定性计算公式具体为：
[0054][0055]
p为主槽稳定性参数，li为第i个主槽横截面的平均长度信息，wi为第i个主槽横截面的平均宽度信息，gi为第i个主槽横截面的预估平均淤积量，u为调整项，a为第i个主槽横截面的面积收敛系数，o为预估加权淤积量，ω为根据当前主槽流量信息的预估淤积信息，t为根据当前主槽流量信息的预估淤积信息的权重调整信息，ε为主槽稳定性参数的修正项。
[0056]
本实施例通过主槽稳定性计算，以生成主槽稳定性参数，从而根据主槽稳定性参
数生成主槽稳定性报告，以实现科学可靠的河口主槽治理。
[0057]
本实施例提供一种主槽稳定性计算公式，该公式充分考虑了第i个主槽横截面的平均长度信息li、第i个主槽横截面的平均宽度信息wi、第i个主槽横截面的预估平均淤积量gi、调整项u、第i个主槽横截面的面积收敛系数a、预估加权淤积量o、根据当前主槽流量信息的预估淤积信息ω、根据当前主槽流量信息的预估淤积信息的权重调整信息t，以及相互之间的作用关系，以形成函数关系并通过主槽稳定性参数的修正项ε进行修正，从而提供准确可靠的数据支持。
[0058]
在本说明书的一个实施例中，其中根据主槽稳定性参数，生成主槽稳定性报告的步骤包括以下步骤：
[0059]
步骤s51：判断主槽稳定性参数是否大于预设的第一主槽稳定性阈值；
[0060]
步骤s52：确定主槽稳定性参数大于预设的第一主槽稳定性阈值时，生成高等级稳定性主槽报告；
[0061]
步骤s53：确定主槽稳定性参数小于或等于预设的第一主槽稳定性阈值时，判断主槽稳定性参数是否大于第二主槽稳定性阈值；
[0062]
步骤s54：确定主槽稳定性参数大于预设的第二主槽稳定性阈值时，生成中等级稳定性主槽报告；
[0063]
步骤s55：确定主槽稳定性参数小于或等于预设的第二主槽稳定性阈值时，生成低等级稳定性主槽报告。
[0064]
本实施例通过不同的主槽稳定性跟预设的不同的主槽稳定性阈值的比对，以生成不同的稳定性报告，从而为实施人员提供直观的河口治理方案报告。
[0065]
本发明不进行人为设计路线，采用物理模型动床试验的方式，在径流、潮汐动力的作用下，基于水沙相互作用、自我演变调节为主导，分析河口水域在自然动力的作用下，哪个区域会发展出主槽，根据自然演变的规律，按照自然深槽的走向，设计河口治理主槽的分布位置，再按照冲深，设计河底的开挖高程，以实现主槽治理方案与河口的水沙动力分布、河势格局相适应。
附图说明
[0066]
通过阅读参照以下附图所作的对非限制性实施所作的详细描述，本技术的其它特征、目的和优点将会变得更明显：
[0067]
图1示出了一实施例的一种基于水沙动力自适应特性的河口治理模拟方法的步骤流程图；
[0068]
图2示出了一实施例的河口物理模型构建方法的步骤流程图；
[0069]
图3a-3b示出了一实施例的河口图像特征运动数据集生成方法的步骤流程图；
[0070]
图4示出了一实施例的河床冲刷形态数据生成方法的步骤流程图；
[0071]
图5示出了一实施例的稳定性主槽报告生成方法的步骤流程图；
[0072]
图6示出了一实施例的试验流态示踪的效果展示图；
[0073]
图7示出了一实施例的1号断面冲刷前后泥面线的效果展示图；
[0074]
图8示出了一实施例的2号断面冲刷前后泥面线的效果展示图；
[0075]
图9示出了一实施例的3号断面冲刷前后泥面线的效果展示图；
[0076]
图10示出了一实施例的4号断面冲刷前后泥面线的效果展示图。
具体实施方式
[0077]
下面结合附图对本发明专利的技术方法进行清楚、完整的描述，显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域所属的技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0078]
此外，附图仅为本发明的示意性图解，并非一定是按比例绘制。图中相同的附图标记表示相同或类似的部分，因而将省略对它们的重复描述。附图中所示的一些方框图是功能实体，不一定必须与物理或逻辑上独立的实体相对应。可以采用软件形式来实现功能实体，或在一个或多个硬件模块或集成电路中实现这些功能实体，或在不同网络和/或处理器方法和/或微控制器方法中实现这些功能实体。
[0079]
应当理解的是，虽然在这里可能使用了术语“第一”、“第二”等等来描述各个单元，但是这些单元不应当受这些术语限制。使用这些术语仅仅是为了将一个单元与另一个单元进行区分。举例来说，在不背离示例性实施例的范围的情况下，第一单元可以被称为第二单元，并且类似地第二单元可以被称为第一单元。这里所使用的术语“和/或”包括其中一个或更多所列出的相关联项目的任意和所有组合。
[0080]
请参阅图1至图10，一种基于水沙动力自适应特性的河口治理模拟方法，包括以下步骤：
[0081]
步骤s1：获取河口基础数据，并根据河口基础数据构建河口物理模型或数学模型；
[0082]
步骤s2：开展物理模型清水动床冲刷试验或数学模型清水动床冲刷模拟计算，获取冲刷后的河床冲刷形态数据；
[0083]
步骤s3：根据河床冲刷后的形态数据生成河口治理方案；
[0084]
步骤s4：根据河口治理方案构建河口治理模型，并根据河口治理模型进行物理模型浑水定床淤积试验或数学模型悬沙淤积模拟计算，根据主槽淤积特性，生成主槽淤积特征数据；
[0085]
步骤s5：对主槽淤积特征数据进行稳定性计算，生成主槽稳定性参数，以供河口治理。
[0086]
具体地，例如
①
围填工程实施后，主槽被切断，开展清水动床试验，研究河床冲刷形态。
[0087]
②
根据冲刷后主槽的分布位置，制定河口治理方案。
[0088]
③
将河口治理方案建模，开展浑水定床试验，分析主槽的淤积特性。
[0089]
④
研究该方法确定的主槽是否稳定，是否符合区域的水动力特性。
[0090]
本实施例不进行人为设计路线，采用物理模型动床试验的方式，在径流、潮汐动力的作用下，基于水沙相互作用、自我演变调节为主导，分析河口水域在自然动力的作用下，哪个区域会发展出主槽，根据自然演变的规律，按照自然深槽的走向，设计河口治理主槽的分布位置，再按照冲深，设计河底的开挖高程。
[0091]
在本说明书的一个实施例中，请参阅图2，步骤s1包括以下步骤：
[0092]
步骤s11：获取原始河口基础数据，其中原始河口基础数据包括原始潮流流向数据、原始涨潮潮流数据、原始落潮潮流数据以及原始潮流间拓扑关系数据；
[0093]
具体地，例如生成可视化用户输入界面，以获取原始河口基础数据。
[0094]
步骤s12：获取原始河口图像信息序列以及修建堤防后获取的堤防河口图像信息序列，并对原始河口图像信息序列以及堤防河口图像信息序列进行特征提取，生成原始河口图像特征运动数据集以及堤防河口图像特征运动数据集；
[0095]
具体地，例如通过如下图的历史测试数据进行测试，获取测试图像信息；
[0096]
模拟试验径潮水文组合
[0097][0098]
具体地，例如通过cfog算法生成图像特征运动数据集。
[0099]
步骤s13：根据原始河口图像特征运动数据集以及堤防河口图像特征运动数据集通过潮流系数加权计算公式进行计算，生成潮流加权系数，并将原始河口基础数据通过潮流加权系数进行加权计算并进行率定，生成河口基础数据；
[0100]
步骤s14：根据河口基础数据生成河口水动力特性数据；
[0101]
步骤s15：根据河口水动力特性数据以及河口基础数据构建河口物理模型或数学模型。
[0102]
具体地，例如所有上、下游边界的过渡段都按实测地形模拟，保证上下游可较好的模拟径流、潮汐过程，尽可能的保证模型与原型相似。
[0103]
根据试验场地面积、供水能力，选定模型的平面比尺为1:700。垂直比尺根据模型水流处于紊流阻力平方区的水深比尺条件来确定，并参考国内大型河口模型设计，确定模型的垂直比尺为1:100，相应的变率为7；水流时间比尺为1:70，冲刷时间比尺为1:559。
[0104]
本实施例通过图像数据分析，通过原始河口图像信息以及设立水利工程后的堤防河口图像信息进行比对计算，生成加权系数，以根据原始河口图像信息的记录参数进行加权计算，从而生成稳定可靠的河口基础数据，以建立准确的河口物理模型。
[0105]
在本说明书的一个实施例中，其中潮流系数加权计算公式具体为：
[0106]
[0107]
p
t
为潮流加权系数，αi为原始河口图像特征运动数据集中的第i个原始河口图像特征运动数据，αi为原始河口图像特征运动数据集中的第i个原始河口图像特征运动数据相应的权重信息，bi为堤防河口图像特征运动数据集中的第i个堤防河口图像特征运动数据，βi为堤防河口图像特征运动数据集中的第i个堤防河口图像特征运动数据的参数信息的权重信息，ui为原始河口图像特征运动数据集中的第i个原始河口图像特征运动数据对应的流量速度信息，θ为调整项的权重信息，为调整项的误差系数，x为根据标准河口图像特征生成的调整项，i的取值为1、2、3...n，n为潮流加权系数的修正项。
[0108]
本实施例提供一种潮流系数加权计算公式，该公式充分考虑了原始河口图像特征运动数据集中的第i个原始河口图像特征运动数据ai、原始河口图像特征运动数据集中的第i个原始河口图像特征运动数据相应的权重信息αi、堤防河口图像特征运动数据集中的第i个堤防河口图像特征运动数据bi、堤防河口图像特征运动数据集中的第i个堤防河口图像特征运动数据的参数信息的权重信息βi、原始河口图像特征运动数据集中的第i个原始河口图像特征运动数据对应的流量速度信息ui、调整项的权重信息θ、调整项的误差系数根据标准河口图像特征生成的调整项x以及之间的相互关系，以形成函数关系并通过潮流加权系数的修正项n进行修正，从而提供准确可靠的数据支撑。
[0109]
在本说明书的一个实施例中，请参阅图3a-3b，步骤s12包括以下步骤：
[0110]
步骤s121：获取原始河口图像信息序列以及堤防河口图像信息序列；
[0111]
具体地，例如通过高清运动摄像头获取图像信息。
[0112]
步骤s122：对原始河口图像信息序列以及堤防河口图像信息序列进行关键像素点提取，生成原始河口图像信息序列的原始河口潮流分布特征点集序列以及堤防河口图像信息序列的堤防河口潮流分布特征点集序列；
[0113]
具体地，例如加入罗丹明染色剂，观测落潮时径流输移路线，如图6所示，对超过预设的rbg值的图像像素点进行识别并提取。
[0114]
步骤s123：将原始河口图像信息序列中的当前原始河口图像信息的原始河口潮流分布特征点集与原始河口图像信息序列中的下一张原始河口图像信息的原始河口潮流分布特征点集进行匹配，生成原始河口潮流匹配特征点集，并将堤防河口图像信息序列中的当前堤防河口图像信息的堤防河口潮流分布特征点集与堤防河口图像信息序列中的下一张堤防河口图像信息的堤防河口潮流分布特征点集进行匹配，生成堤防河口潮流匹配特征点集；
[0115]
具体地，例如通过多模态图像融合匹配算法进行计算，从而生成匹配的特征点集。
[0116]
步骤s124：根据步骤s123的方式依次将原始河口图像信息序列中的原始河口潮流分布特征点集与下一张原始河口图像信息的原始河口潮流分布特征点集进行匹配以及堤防河口图像信息序列中的堤防河口潮流分布特征点集与下一张堤防河口图像信息序列中的堤防河口潮流分布特征点集进行匹配，生成原始河口图像全匹配特征点集以及堤防河口图像全匹配特征点集，且根据原始河口图像全匹配集在原始河口图像信息中的相对位置以及堤防河口图像全匹配特征点集在堤防河口图像信息中的相对位置分布进行位置空间计算，生成原始河口图像特征运动数据集以及堤防河口图像特征运动数据集。
[0117]
具体地，例如其中位置空间计算通过将原图像分布特征点进行匹配后进行图像相对位置、帧数时间差以及实际等比例转化，生成图像特征运动数据集。
[0118]
本实施例通过对高清运动摄像头拍摄的图像信息序列进行依次比对，依次对相近两张的图像信息的特征进行比对，从而生成图像运动特征信息，并根据在图像的相对位置以及预设的参数分析方式，从而生成图像运动数据信息，提供准确的水动力分布特征信息。
[0119]
在本说明书的一个实施例中，请参阅图4，步骤s2包括以下步骤：
[0120]
步骤s21：获取清水流量数据；
[0121]
具体地，例如生成用户输入界面，并通过用户输入界面获取清水流量数据。
[0122]
步骤s22：基于河口物理模型或数学模型，通过清水动床试验的河口床沙计算公式对清水流量数据进行河口流量与河床推移质的运动关系的模拟计算，生成河口断面流量流速推移质运动数据；
[0123]
步骤s23：根据河口断面流量流速推移质运动数据对河口物理模型或数学模型进行修正，获取冲刷后的河床冲刷形态数据。
[0124]
具体地，例如在无开挖(现状地形)情况下，a河口现状主槽被围填，在常年洪水作用下，会形成新的主槽，a河出口处主槽靠岸，至下游出口处略微向西偏转；如图7至图10所示，根据河口淤积的起点距、高度绘制断面冲刷前后泥面线，其中实线为现状河床，虚线为新岸线在自然冲刷下的结果，1号、2号、3号断面位于现状河出口深槽被占用段，4号断面位于工程下游段。1号断面深泓线靠左岸，断面最深位置-3.8m，约40％的深槽被侵占；岸线外推后，在自然冲刷演变条件下，深泓线西侧滩面冲刷变深，深槽向西侧发展、扩宽，但深槽仍然靠近左侧堤岸，冲刷后河床深槽底高程-3.3m。2号断面位置，现状深槽最深位置为-3.3m；围填后-3m深槽全部被侵占，冲刷试验后，近岸20-190m范围内冲刷形成-2m深槽，深泓线最低高程-2.6m，深泓线距离堤脚位置约100m，深槽基本贴岸分布。3号断面位置，现状-3m深槽距离岸线约490-680m，最深位置-3.6m，近岸被侵占后，深槽大部分得以保留；冲刷后，深泓线向西摆动约70m，近岸发生淤积，西侧浅滩贴近深槽段形成冲刷，断面最低点高程-3.4m，-3m深槽缩窄约65m。4号断面位置，从岸线向西侧水域逐渐变深，现状滩面无明显深槽分布；冲刷试验后，河床整体表现为冲刷，近岸滩面平均下切约0.09m/a。
[0125]
本实施例通过清水动床试验计算，实现对断面形态变化的模拟，从而分析水沙自适应特性，研究深槽在自然动力作用下的发育、演变和分布。
[0126]
在本说明书的一个实施例中，其中河口床沙计算公式具体为：
[0127][0128]
q为河口断面流量流速推移质运动数据，v为清水流速信息，f为清水流量水动力含沙平均加权量，p为清水流速信息对应的平均含沙量，z为清水流量水动力的调整项，r为清水流速信息的平均含沙量，为清水流速信息的平均含沙量的实际误差调整项，为河口断面流量流速推移质运动数据的修正项。
[0129]
本实施例提供一种河口床沙计算公式，该公式充分考虑清水流速信息v、清水流量水动力含沙平均加权量f、清水流速信息对应的平均含沙量p、清水流量水动力的调整项z、清水流速信息的平均含沙量清水流速信息的平均含沙量的实际误差调整项r以及相互之
间的作用关系，以形成函数关系并通过河口断面流量流速推移质运动数据的修正项进行修正，从而提供准确可靠的数据支撑。
[0130]
在本说明书的一个实施例中，步骤s3包括以下步骤：
[0131]
根据河床冲刷后的形态数据中的河床主槽数据以及河床深泓线分布数据，生成河口预治理方案；
[0132]
具体地，例如通过试验，与现状相比，新冲刷的深槽束窄，-2m深槽宽度150-280m，-3m深槽宽度70-130m；现状深泓线位于-3.3至-4.2m，新岸线条件下，深槽变浅，介于-2.6至-3.4m。
[0133]
获取河道整治历史经验数据，并根据河道整治历史经验数据对河口预治理方案进行调整，以生成河口治理方案。
[0134]
具体地，例如结合新冲刷后的主槽、深泓线分布，制订河口治理方案。根据河道整治经验，深槽开挖时预留回淤空间，根据现状及冲刷试验后深泓线分布，可将挖深初定在-4m。
[0135]
整治起点为河出口，与上游深槽对接，下游接现状-4m等高线，全长约4.0km，推荐深槽开挖至-4.0m，基本沿着本次研究成果中新冲刷出现的深槽布置，上游底宽150m，下游底宽200m，边坡1:7。
[0136]
本实施例通过根据河床冲刷后的形态数据中的河床主槽数据以及河床深泓线分布数据，生成河口预治理方案，并获取河道整治历史经验数据，并根据河道整治历史经验数据对河口预治理方案进行调整，以生成河口治理方案，从而为下一步做好前提准备工作。
[0137]
在本说明书的一个实施例中，步骤s4包括以下步骤：
[0138]
根据河口治理方案构建河口治理模型；
[0139]
根据河口治理模型进行物理模型浑水定床淤积试验或数学模型悬沙淤积模拟计算，根据主槽淤积特性，生成主槽淤积特征数据，其中浑水定床试验计算的步骤为通过浑水淤积计算公式进行计算；
[0140]
具体地，例如1号断面，-4m深槽底宽159m，开挖主槽距离岸线50至209m；淤积试验后，平均淤积厚度0.75m，-3m深槽宽146m，距离岸线37至183m，主槽最深位置高程-3.40m，距离岸线103m。2号断面，-4m深槽底宽175m，开挖主槽距离岸线50至225m；淤积试验后，平均淤积厚度0.68m，-3m深槽宽155m，距离岸线50至205m，主槽最深位置高程-3.56m，距离岸线88m。3号断面，-4m深槽底宽190m，开挖主槽距离岸线50至240m；淤积试验后，平均淤积厚度0.53m，-3m深槽宽185m，距离岸线29至214m，主槽最深位置高程-3.62m，距离岸线119m。4号断面位于治理河段下游900m，断面呈微淤态势。
[0141]
主槽平均淤积厚度单位：m
[0142][0143]
其中浑水淤积计算公式具体为：
[0144][0145]
s为主槽淤积特性数据中的浑水淤积数据，s为浑水流速信息，γ为浑水流速水动力的当前含沙量的调整项，m为浑水流量水动力的调整项，h为浑水流量水动力含沙平均加权量，t为浑水流速信息的平均含沙量的实际误差调整项,c为浑水流速信息的平均含沙量，为浑水流速信息的平均含沙量的当前预估加权信息，为主槽淤积特性数据中的浑水淤积数据的修正项。
[0146]
本实施例通过根据河口治理方案构建河口治理模型，并根据河口治理模型进行物理模型浑水定床淤积试验或数学模型悬沙淤积模拟计算，根据主槽淤积特性，生成主槽淤积特征数据，从而为下一步做好前提准备。
[0147]
本实施例提供一种浑水淤积计算公式，该公式充分考虑浑水流速信息s、浑水流速水动力的当前含沙量的调整项γ、浑水流量水动力的调整项m、浑水流量水动力含沙平均加权量h、浑水流速信息的平均含沙量的实际误差调整项t、浑水流速信息的平均含沙量c、浑水流速信息的平均含沙量的当前预估加权信息以及相互之间的作用关系，以形成函数关系并通过主槽淤积特性数据中的浑水淤积数据的修正项进行修正，以提供准确的数据支持。
[0148]
在本说明书的一个实施例中，步骤s5包括以下步骤：
[0149]
通过主槽稳定性计算公式对主槽淤积特征数据进行计算，生成主槽稳定性参数；
[0150]
具体地，例如通过本实施例提供的主槽稳定性计算公式进行计算。
[0151]
根据主槽稳定性参数生成主槽稳定性报告，以供河口治理；
[0152]
具体地，例如根据历史稳定性参数数据同主槽稳定性参数进行比对，生成主槽稳定性报告。
[0153]
其中主槽稳定性计算公式具体为：
[0154][0155]
p为主槽稳定性参数，li为第i个主槽横截面的平均长度信息，wi为第i个主槽横截面的平均宽度信息，gi为第i个主槽横截面的预估平均淤积量，u为调整项，a为第i个主槽横截面的面积收敛系数，o为预估加权淤积量，ω为根据当前主槽流量信息的预估淤积信息，t
为根据当前主槽流量信息的预估淤积信息的权重调整信息，ε为主槽稳定性参数的修正项。
[0156]
本实施例通过主槽稳定性计算，以生成主槽稳定性参数，从而根据主槽稳定性参数生成主槽稳定性报告，以实现科学可靠的河口主槽治理。
[0157]
本实施例提供一种主槽稳定性计算公式，该公式充分考虑了第i个主槽横截面的平均长度信息li、第i个主槽横截面的平均宽度信息wi、第i个主槽横截面的预估平均淤积量gi、调整项u、第i个主槽横截面的面积收敛系数a、预估加权淤积量o、根据当前主槽流量信息的预估淤积信息ω、根据当前主槽流量信息的预估淤积信息的权重调整信息t，以及相互之间的作用关系，以形成函数关系并通过主槽稳定性参数的修正项ε进行修正，从而提供准确可靠的数据支持。
[0158]
在本说明书的一个实施例中，请参阅图5，其中根据主槽稳定性参数，生成主槽稳定性报告的步骤包括以下步骤：
[0159]
步骤s51：判断主槽稳定性参数是否大于预设的第一主槽稳定性阈值；
[0160]
具体地，例如判断主槽稳定性参数如23.5是否大于预设的第一主槽稳定性阈值如60。
[0161]
步骤s52：确定主槽稳定性参数大于预设的第一主槽稳定性阈值时，生成高等级稳定性主槽报告；
[0162]
具体地，例如确定主槽稳定性参数如70大于预设的第一主槽稳定性阈值如60时，生成高等级稳定性主槽报告。
[0163]
步骤s53：确定主槽稳定性参数小于或等于预设的第一主槽稳定性阈值时，判断主槽稳定性参数是否大于第二主槽稳定性阈值；
[0164]
具体地，例如确定主槽稳定性参数如56小于或等于预设的第一主槽稳定性阈值如60时，判断主槽稳定性参数是否大于第二主槽稳定性阈值。
[0165]
步骤s54：确定主槽稳定性参数大于预设的第二主槽稳定性阈值时，生成中等级稳定性主槽报告；
[0166]
具体地，例如确定主槽稳定性参数如38.3大于预设的第二主槽稳定性阈值如30时，生成中等级稳定性主槽报告。
[0167]
步骤s55：确定主槽稳定性参数小于或等于预设的第二主槽稳定性阈值时，生成低等级稳定性主槽报告。
[0168]
具体地，例如确定主槽稳定性参数如28.1小于或等于预设的第二主槽稳定性阈值如30时，生成低等级稳定性主槽报告。
[0169]
本实施例通过不同的主槽稳定性跟预设的不同的主槽稳定性阈值的比对，以生成不同的稳定性报告，从而为实施人员提供直观的河口治理方案报告。
[0170]
本发明不进行人为设计路线，采用物理模型动床试验的方式，在径流、潮汐动力的作用下，基于水沙相互作用、自我演变调节为主导，分析河口水域在自然动力的作用下，哪个区域会发展出主槽，根据自然演变的规律，按照自然深槽的走向，设计河口治理主槽的分布位置，再按照冲深，设计河底的开挖高程，以实现主槽治理方案与河口的水沙动力分布、河势格局相适应。
[0171]
因此，无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本
发明的范围由所附权利要求而不是上述说明限定，因此旨在将落在权利要求的等同要件的含义和范围内的所有变化涵括在本发明内。不应将权利要求中的任何附关联图标记视为限制所涉及的权利要求。
[0172]
以上所述仅是本发明的具体实施方式，使本领域技术人员能够理解或实现本发明。对这些实施例的多种修改对本领域的技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所发明的原理和新颖特点相一致的最宽的范围。

技术特征：
1.一种基于水沙动力自适应特性的河口治理模拟方法，其特征在于，包括以下步骤：步骤s1：获取河口基础数据，并根据河口基础数据构建河口物理模型或数学模型；步骤s2：开展物理模型清水动床冲刷试验或数学模型清水动床冲刷模拟计算，获取冲刷后的河床冲刷形态数据；步骤s3：根据河床冲刷后的形态数据生成河口治理方案；步骤s4：根据河口治理方案构建河口治理模型，并根据河口治理模型进行浑水定床淤积试验或数学模型模拟计算，根据主槽淤积特性，生成主槽淤积特征数据；步骤s5：对主槽淤积特征数据进行稳定性计算，生成主槽稳定性参数，以供河口治理。2.根据权利要求1所述的方法，其特征在于，步骤s1包括以下步骤：步骤s11：获取原始河口基础数据，其中原始河口基础数据包括原始潮流流向数据、原始涨潮潮流数据、原始落潮潮流数据以及原始潮流间拓扑关系数据；步骤s12：获取原始河口图像信息序列以及修建堤防后获取的堤防河口图像信息序列，并对原始河口图像信息序列以及堤防河口图像信息序列进行特征提取，生成原始河口图像特征运动数据集以及堤防河口图像特征运动数据集；步骤s13：根据原始河口图像特征运动数据集以及堤防河口图像特征运动数据集通过潮流系数加权计算公式进行计算，生成潮流加权系数，并将原始河口基础数据通过潮流加权系数进行加权计算并进行率定，生成河口基础数据；步骤s14：根据河口基础数据生成河口水动力特性数据；步骤s15：根据河口水动力特性数据以及河口基础数据构建河口物理模型或数学模型。3.根据权利要求2所述的方法，其特征在于，其中潮流系数加权计算公式具体为：p
t
为潮流加权系数，a
i
为原始河口图像特征运动数据集中的第i个原始河口图像特征运动数据，α
i
为原始河口图像特征运动数据集中的第i个原始河口图像特征运动数据相应的权重信息，b
i
为堤防河口图像特征运动数据集中的第i个堤防河口图像特征运动数据，β
i
为堤防河口图像特征运动数据集中的第i个堤防河口图像特征运动数据的参数信息的权重信息，u
i
为原始河口图像特征运动数据集中的第i个原始河口图像特征运动数据对应的流量速度信息，θ为调整项的权重信息，为调整项的误差系数，x为根据标准河口图像特征生成的调整项，i的取值为1、2、3...n，n为潮流加权系数的修正项。4.根据权利要求2所述的方法，其特征在于，步骤s12包括以下步骤：步骤s121：获取原始河口图像信息序列以及堤防河口图像信息序列；步骤s122：对原始河口图像信息序列以及堤防河口图像信息序列进行关键像素点提取，生成原始河口图像信息序列的原始河口潮流分布特征点集序列以及堤防河口图像信息序列的堤防河口潮流分布特征点集序列；步骤s123：将原始河口图像信息序列中的当前原始河口图像信息的原始河口潮流分布特征点集与原始河口图像信息序列中的下一张原始河口图像信息的原始河口潮流分布特征点集进行匹配，生成原始河口潮流匹配特征点集，并将堤防河口图像信息序列中的当前堤防河口图像信息的堤防河口潮流分布特征点集与堤防河口图像信息序列中的下一张堤防河口图像信息的堤防河口潮流分布特征点集进行匹配，生成堤防河口潮流匹配特征点
集；步骤s124：根据步骤s123的方式依次将原始河口图像信息序列中的原始河口潮流分布特征点集与下一张原始河口图像信息的原始河口潮流分布特征点集进行匹配以及堤防河口图像信息序列中的堤防河口潮流分布特征点集与下一张堤防河口图像信息序列中的堤防河口潮流分布特征点集进行匹配，生成原始河口图像全匹配特征点集以及堤防河口图像全匹配特征点集，且根据原始河口图像全匹配集在原始河口图像信息中的相对位置以及堤防河口图像全匹配特征点集在堤防河口图像信息中的相对位置分布进行位置空间计算，生成原始河口图像特征运动数据集以及堤防河口图像特征运动数据集。5.根据权利要求1所述的方法，其特征在于，步骤s2包括以下步骤：获取清水流量数据；基于河口物理模型或数学模型，通过清水动床试验的河口床沙计算公式对清水流量数据进行河口流量与河床推移质的运动关系的模拟计算，生成河口断面流量流速推移质运动数据；根据河口断面流量流速推移质运动数据对河口物理模型或数学模型进行修正，获取冲刷后的河床冲刷形态数据。6.根据权利要求5所述的方法，其特征在于，其中河口床沙计算公式具体为：q为河口断面流量流速推移质运动数据，v为清水流速信息，f为清水流量水动力含沙平均加权量，p为清水流速信息对应的平均含沙量，z为清水流量水动力的调整项，r为清水流速信息的平均含沙量，为清水流速信息的平均含沙量的实际误差调整项，为河口断面流量流速推移质运动数据的修正项。7.根据权利要求1所述的方法，其特征在于，步骤s3包括以下步骤：根据河床冲刷后的形态数据中的河床主槽数据以及河床深泓线分布数据，生成河口预治理方案；获取河道整治历史经验数据，并根据河道整治历史经验数据对河口预治理方案进行调整，以生成河口治理方案。8.根据权利要求1所述的方法，其特征在于，步骤s4包括以下步骤：根据河口治理方案构建河口治理模型；根据河口治理模型进行物理模型浑水定床淤积试验或数学模型悬沙淤积模拟计算，根据主槽淤积特性，生成主槽淤积特征数据，其中浑水定床试验计算的步骤为通过浑水淤积计算公式进行计算；其中浑水淤积计算公式具体为：s为主槽淤积特性数据中的浑水淤积数据，s为浑水流速信息，γ为浑水流速水动力的当前含沙量的调整项，m为浑水流量水动力的调整项，h为浑水流量水动力含沙平均加权量，t为浑水流速信息的平均含沙量的实际误差调整项，c为浑水流速信息的平均含沙量，为浑水流速信息的平均含沙量的当前预估加权信息，为主槽淤积特性数据中的浑水淤积数
据的修正项。9.根据权利要求1所述的方法，其特征在于，步骤s5包括以下步骤：通过主槽稳定性计算公式对主槽淤积特征数据进行计算，生成主槽稳定性参数；根据主槽稳定性参数生成主槽稳定性报告，以供河口治理；其中主槽稳定性计算公式具体为：p为主槽稳定性参数，l
i
为第i个主槽横截面的平均长度信息，w
i
为第i个主槽横截面的平均宽度信息，g
i
为第i个主槽横截面的预估平均淤积量，u为调整项，a为第i个主槽横截面的面积收敛系数，o为预估加权淤积量，ω为根据当前主槽流量信息的预估淤积信息，t为根据当前主槽流量信息的预估淤积信息的权重调整信息，ε为主槽稳定性参数的修正项。10.根据权利要求9所述的方法，其特征在于，其中根据主槽稳定性参数，生成主槽稳定性报告的步骤包括以下步骤：判断主槽稳定性参数是否大于预设的第一主槽稳定性阈值；确定主槽稳定性参数大于预设的第一主槽稳定性阈值时，生成高等级稳定性主槽报告；确定主槽稳定性参数小于或等于预设的第一主槽稳定性阈值时，判断主槽稳定性参数是否大于第二主槽稳定性阈值；确定主槽稳定性参数大于预设的第二主槽稳定性阈值时，生成中等级稳定性主槽报告；确定主槽稳定性参数小于或等于预设的第二主槽稳定性阈值时，生成低等级稳定性主槽报告。

技术总结
本发明涉及河口治理分析决策技术领域，尤其涉及一种基于水沙动力自适应特性的河口治理模拟方法。该方法包括以下步骤：获取河口基础数据，并根据河口基础数据构建河口物理模型或数学模型；开展物理模型清水动床试验或数学模型清水动床冲刷模拟计算，获取冲刷后的河床形态数据；根据形态数据生成河口治理方案；根据河口治理方案构建河口治理模型，并根据河口治理模型进行浑水定床淤积试验或数学模型模拟计算，根据主槽淤积特性，获取主槽淤积特性数据；对主槽淤积特征数据进行稳定性计算，生成主槽稳定性参数，以供河口治理。本发明通过河口水沙动力分布特性模拟，实现河床冲、淤强度的自适应性，弱化人工干预，以提供稳定可靠的河口治理方案。的河口治理方案。的河口治理方案。


技术研发人员：刘国珍 佟晓蕾 温舒茵 袁菲 吴尧 杨裕桂 李慧婧 龙晓飞 陈娟 林平 王海俊 黄淞宣 刘佳琪 罗陈威 陈浩南
受保护的技术使用者：珠江水利委员会珠江水利科学研究院
技术研发日：2023.01.12
技术公布日：2023/7/22