基于回归分析的坝区流态的推求方法

1.本发明属于计算水力学领域，基于回归分析的坝区流态的推求方法。


背景技术：

2.大坝作为水利工程枢纽的重要组成部分，不仅对防洪发电等发挥着巨大的工程效益，同时对河道生态系统也有着重大影响，在水生生物生境和人类生产生活所依靠的社会、经济、文化、生态等方面扮演重要的角色。大坝作为承担挡水功能的水利枢纽，长期在河道中运行，其对河道水流条件产生了重大的影响，尤其是靠近坝体的河道部分。推求坝体附近区域流态分布情况是人们了解坝体所在区域河道水沙流动状态和保证水生生物生存环境及通航等亲水活动安全的有效手段，也是保证大坝运行安全的非工程措施。为及时掌握水流流态，及时了解存在的河道运行风险，通常采用二维水动力模型模拟坝区水流流态，以获得坝上及坝下流速分布。
3.研究发现，影响水流流态的因素是多方面的，而最直接的影响因素是河道断面形态，其主要参数包括泥沙淤积带来的河床变化、蓄水期水位带来的水位变化、水利枢纽的布设带来的湿周的变化，这些参数的变化直接导致了水流在河道不同位置处的垂向流速分布不同，如在宽谷河道与狭窄河道都呈现指数型的分布，且分布系数不同，又如坝区河道中间呈现指数型分布，靠岸的部分水流垂向流速分布不呈现指数型。
4.现有二维水流流态计算过程中，存在的主要问题是：河段中任意位置垂向流速分布概化为指数分布，且分布系数一致，而没有考虑因水流存在于河道中的位置不同，或者是水利枢纽等人类活动或非人类活动带来的河道条件的变化导致的垂向流速分布的变化，因而在水力计算中得出的坝体上下游河段水流流态与实际流态有一定的出入，准确率不高。


技术实现要素：

5.根据现有技术的不足，本发明的目的是提供基于回归分析的坝区流态的推求方法,能较准确的计算出复杂条件下坝区内坝体上、下游流场分布情况，修正后的结果与实际情况更加符合。
6.为了解决上述技术问题，本发明采用的技术方案为：基于回归分析的坝区流态的推求方法，包括以下步骤：步骤1、确定包含大坝的上、下游目标河段的研究区域，从目标河段内提取各个断面及相应测点的实测垂向流速值，建立二维水流流态计算域地理信息数据库；步骤2、基于实测垂向流速值，对测点进行划分，对不同测点的垂向流速分布曲线进行分类整合，构建不同测点的垂线流速分布模型，得到不同测点表面流速和垂线流速平均分布的比值，作为计算水流表面流速的垂向流速分布边界条件；步骤3、采用深度学习，结合河道位置特征，确定校正域，在所述校正域基于支持向量机分类回归理论，构建非线性分类回归模型，得到垂向流速在研究区域不同位置处的分布特征，用作计算水流表面流速的垂向流速分布边界条件；
步骤4、对研究区域进行网格划分，得到网格节点的高程值，用于后面计算研究区域水流表面流速的地形条件及单元格划分；步骤5、进行水文资料分析，并分类列出边界条件，边界条件包括网格节点的地形条件、岸线走向、网格节点垂向流速分布，根据大坝的校核水位、设计洪水位或发电任务对应的流量以及水位作为计算条件，选取典型的计算条件；步骤6、根据网格节点垂向流速分布得到各自垂向平均流速与表面流速的比值，分别建立大坝的上、下游的平面二维水动力数学模型并进行验证，率定出研究区域不同河道断面的糙率；步骤7、考虑大坝水电站运行的影响，采用步骤5中的边界条件，坝上、下游区域模型进口断面给定流量边界条件，并根据水电的调度方式及计算区域特征给定出口的水位控制条件，为后面计算提供了进口流量边界条件和出口水位初始条件；步骤8、假定初始流速，采用步骤6中建立的平面二维水动力数学模型，按照步骤7中的流量边界条件，计算得出大坝的上、下游水流流态。
7.进一步地，所述步骤1具体包括：步骤11、收集所述目标河段地理信息，包括岸线、河道走向、河段滩涂分布，并计算入口断面及出口断面；步骤12、收集所述目标河段内二维水流测点信息，包括断面方向、断面起点到原点的岸线距离、测线到断面起点的距离以及测点在测线上的位置，建立二维水流流态测点信息数据集。
8.进一步地，所述步骤2中，对每一测点，以该测点不同的水深及对应的流速作为输入值，以垂线流速分布特征因子作为输出值构建垂线流速分布模型，垂线流速分布为指数形式、对数形式、抛物线形式以及椭圆形式。
9.进一步地，所述指数形式表示为：式中：表示测点水深，为系数，表示垂线水深，为该垂线最大流速，为测点流速。
10.进一步地，所述步骤3中，对于非线性分类回归模型，根据所述检验集地理位置自动筛选校正域及所述非线性分类回归模型参数，估测垂向流速分布参数估计误差场，并对所述检验集的背景场进行误差校正和精度评价，在精度评价过程中，有效性和符合预测精度的预测值的选取包括复相关系数和标准差。
11.进一步地，复相关系数的计算公式为：式中，为回归平方和，，表示垂线上测点序号，为垂线上测点数；表示第点处非线性分类回归模型预测的流速，为实测流速的平均值；
为总离差的平方和：式中，表示第点处实测垂向流速值，取值在0～1之间，越靠近1，拟合程度越好；复相关系数采用检验，对用下列统计量其中，和分别为统计量服从分布的第一、二自由度；根据以及、的自由度得到的临界值，、分别为统计量服从分布的第一、二自由度；的范围是1%～5%；当时，则判定非线性分类回归模型的回归曲线有效；当时，则判定非线性分类回归模型回归曲线无效；选取有效性和符合预测精度的垂线流速分布模型。
12.进一步地，所述步骤5中，典型的计算条件为一年一遇的洪水、五十年一遇的洪水或一百年一遇的洪水。
13.进一步地，所述步骤6中，根据曼宁公式率定出研究区域不同河道断面的糙率。
14.进一步地，所述步骤4中，使用反距离插值法对实测地形资料进行插值，得到网格节点的高程值。
15.进一步地，所述步骤3中，在深度学习过程中，确定训练集以及检验集，选取部分测点作为背景场，邻近的部分测点作为检验场，根据所述背景场合所述观测场，计算所述目标河段位置实测垂向流速值的初始误差校正场。
16.与现有技术相比，本发明具有以下优点和有益效果：本发明提供基于回归分析的坝区流态的推求方法，在特殊河道条件下，使用实测垂向流速值，对各个特征点进行了垂向水流分布分析，对于靠近坝前存在的垂向流速分布与宽谷河段流速分布线型不一致的情况，可以提供更加准确的参考；对于各处系数不同的情况，通过采用分类回归的算法，可以对河道不同位置进行分类，从而得出相对准确的各处垂向流速分布。
17.本发明能够将坝体附近河道不同位置垂向流速分布对水流流态的影响分别进行分析，以及能够将水利枢纽、泥沙冲淤等造成的河床和湿周变化以及水位变化等造成的影响通过非线性分类回归方法进行分区域对于坝区上下游水流流态进行预测。
附图说明
18.为了更清楚地说明本实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一个实施例，对于本领域普
通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
19.图1为本发明基于回归分析的坝区流态的推求方法的总体流程图。
20.图2为本发明实施过程的流程图。
21.图3为垂线分布示意图。
22.图4为上游流态计算示例。
23.图5为下游流态计算示例。
具体实施方式
24.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。
25.本发明提供的基于回归分析的坝区流态的推求方法，如图1和图2所示，包括以下步骤：步骤1、确定包含大坝的上、下游目标河段的研究区域，从目标河段内提取各个断面及相应测点的实测垂向流速值，建立二维水流流态计算域地理信息数据库；步骤2、基于实测垂向流速值，对测点进行划分，对不同测点的垂向流速分布曲线进行分类整合，构建不同测点的垂线流速分布模型，得到不同测点表面流速和垂线流速平均分布的比值；步骤3、采用深度学习，结合河道位置特征，确定校正域，在所述校正域基于支持向量机（svm，support vector machine）分类回归理论，构建非线性分类回归模型，得到垂向流速在研究区域不同位置处的分布特征；步骤4、对研究区域进行网格划分，得到网格节点的高程值；步骤5、进行水文资料分析，并分类列出边界条件，边界条件包括网格节点的地形条件、岸线走向、网格节点垂向流速分布，根据大坝的校核水位、设计洪水位或发电任务对应的流量以及水位作为计算条件，选取典型的计算条件；步骤6、根据网格节点垂向流速分布得到各自垂向平均流速与表面流速的比值，分别建立大坝的上、下游的平面二维水动力数学模型并进行验证，率定出研究区域不同河道断面的糙率；步骤7、考虑大坝水电站运行的影响，采用步骤5中的边界条件，坝上、下游区域模型进口断面给定流量边界条件，并根据水电的调度方式及计算区域特征给定出口的水位控制条件；步骤8、假定初始流速，采用步骤6中建立的平面二维水动力数学模型，按照步骤7中的流量边界条件，计算得出大坝的上、下游水流流态。
26.本发明提供的基于回归分析的坝区流态的推求方法，使用实测垂向流速值，对各个特征点进行了垂向水流分布分析，对于靠近坝前存在的垂向流速分布与宽谷河段流速分布线型不一致的情况，可以提供更加准确的参考；对于各处系数不同的情况，通过采用分类回归的算法，可以对河道不同位置进行分类，从而得出相对准确的各处垂向流速分布。
27.本发明能够将坝体附近河道不同位置垂向流速分布对水流流态的影响分别进行分析，以及能够将水利枢纽、泥沙冲淤等造成的河床和湿周变化以及水位变化等造成的影响通过非线性分类回归方法进行分区域对于坝区上下游水流流态进行预测。
28.本发明中，所述步骤1具体包括：步骤11、收集所述目标河段地理信息，包括岸线、河道走向、河段滩涂分布，并计算入口断面及出口断面；步骤12、收集所述目标河段内二维水流测点信息，包括断面方向、断面起点到原点的岸线距离、测线到断面起点的距离以及测点在测线上的位置，建立二维水流流态测点信息数据集。
29.本发明中，所述步骤2中，对每一测点，以该测点不同的水深及对应的流速作为输入值，以垂线流速分布特征因子作为输出值构建垂线流速分布模型，垂线流速分布为指数形式、对数形式、抛物线形式以及椭圆形式。
30.其中，最常用的垂线流速分布模型为指数形式，所述指数形式表示为：式中：表示测点水深，为系数，表示垂线水深，为该垂线最大流速，为测点流速；本发明中，所述步骤3中，对于非线性分类回归模型，根据所述检验集地理位置自动筛选校正域及所述非线性分类回归模型参数，估测垂向流速分布参数估计误差场，并对所述检验集的背景场进行误差校正和精度评价，在精度评价过程中，有效性和符合预测精度的预测值的选取包括复相关系数和标准差。
31.具体地，复相关系数的计算公式为：式中，为回归平方和，，表示垂线上测点序号，为垂线上测点数；表示第点处非线性分类回归模型预测的流速，为实测流速的平均值；为总离差的平方和：式中，表示第点处实测垂向流速值，取值在0～1之间，越靠近1，拟合程度越好；复相关系数采用检验，对用下列统计量其中，和分别为统计量服从分布的第一、二自由
度；根据以及、的自由度得到的临界值，、分别为统计量服从分布的第一、二自由度；的范围是1%～5%；当时，则判定非线性分类回归模型的回归曲线有效；当时，则判定非线性分类回归模型回归曲线无效；选取有效性和符合预测精度的垂线流速分布模型。
32.如图3所示，为有效性和符合预测精度的指数分布的垂线流速分布模型，其中方向表示流速方向，方向表示垂向，坐标原点表示河床底部任一点，为水深，为垂线平均流速，为测点在垂线上的位置。
33.本发明中，所述步骤5中，典型的计算条件为一年一遇的洪水、五十年一遇的洪水或一百年一遇的洪水，具体由计算河段水流情况决定。
34.所述步骤4中，使用反距离插值法对实测地形资料进行插值，得到网格节点的高程值。
35.所述步骤6中，根据曼宁公式率定出研究区域不同河道断面的糙率。
36.所述步骤3中，在深度学习过程中，确定训练集以及检验集，选取部分测点作为背景场，邻近的部分测点作为检验场，根据所述背景场合所述观测场，计算所述目标河段位置实测垂向流速值的初始误差校正场。
37.在本发明的一个实施例中，如图4所示，为本发明的上游流态计算示例，如图5所示为本发明的下游流态计算示例，各网格点加入了垂向流速分布边界条件之后，表面流速计算结果更加精确。
38.显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。