首页 / 专利库 / 化学元素和化合物 / / 树状河网梯级水库水动力水质泥沙耦合模拟方法及系统

树状河网梯级库水动水质泥沙耦合模拟方法及系统

阅读:861发布:2021-09-19

专利汇可以提供树状河网梯级库水动水质泥沙耦合模拟方法及系统专利检索,专利查询,专利分析的服务。并且本 发明 涉及 水 利研究技术领域,具体涉及一种树状河网 梯级 水库水动 力 水质泥沙耦合模拟方法及系统,本发明解决河道地形断面和控制线水文 站点 资料的完整性和可靠性;加入梯级水库后如何解决河网与梯级水库的联解,依托最新实测干支流河道水深地形资料,建立了将梯级水库调度模式和一维河网水动力模型融合为一体的西江中游河网整体 水动力学 模拟计算平台。模型计算以一维圣维南方程为 基础 ,提出了从河网、子河网、河道、 节点 四个层面建立模型的计算顺序,对河网及梯级水库进行概化,通过水利枢纽调度方式、外部边界条件的传递关系连接成一个整体,将各个河段连接成一个计算整体,实现同步实时模拟和计算,丰富了河网梯级水库调度的研究手段。,下面是树状河网梯级库水动水质泥沙耦合模拟方法及系统专利的具体信息内容。

1.一种树状河网梯级库水动水质泥沙耦合模拟方法,其特征在于,所述方法包括以下步骤:
S1 根据梯级水库在树状河网中的分布,分割得到河网的控制节点
S2 对河网的控制节点进行分类;
S3 编制程序计算控制节点,实现河网与梯级水库的联合求解;
S4 通过河网整体水动力学模拟计算平台,实现同步实时模拟和计算。
2.根据权利要求1所述的树状河网梯级水库水动力水质泥沙耦合模拟方法,其特征在于,所述S1中,根据梯级水库在树状河网中的分布,将每个水库大坝所在断面作为分割整体河网的控制要素,分隔出若干子河网,按照从上游向下游的顺序来编子河网的编号,每个子河网的控制节点也按从上游向下游编排,大坝上游端和下游端各布置一个控制节点且编号连续,作为每个子河网的上边界或下边界控制节点。
3.根据权利要求1所述的树状河网梯级水库水动力水质泥沙耦合模拟方法,其特征在于,所述S2中,将河网的控制节点进行分类,一类是需要外部给定边界条件的节点,称为外边界节点;另一类为内边界节点,无需给定边界条件,是河道交汇口,外边界节点是河道的上、下游端和大坝前端的节点。
4.根据权利要求1所述的树状河网梯级水库水动力水质泥沙耦合模拟方法,其特征在于,所述S3中,编制程序计算的过程中,按照子河网的编号顺序计算,每个子河网按照树状河网的算法独立运算,运算结束后,根据大坝上游端和下游端控制节点间的关系,将该子河网的下游端计算结果通过大坝上、下游端的控制节点传入给下一个子河网,作为下一个子河网的上游端边界条件,依次类推,最后实现河网与梯级水库的联合求解。
5.根据权利要求1所述的树状河网梯级水库水动力水质泥沙耦合模拟方法,其特征在于,所述S4中,模型计算以一维圣维南方程为基础,从河网、子河网、河道、节点四个层面建立模型的计算顺序,对河网及梯级水库进行概化,通过水利枢纽调度方式、外部边界条件的传递关系连接成一个整体,从而将各个河段连接成一个计算整体。
6.根据权利要求5所述的树状河网梯级水库水动力水质泥沙耦合模拟方法,其特征在于,所述一维圣维南方程包括水流连续方程和水流运动方程;
所述水流连续方程为:
所述水流运动方程为:
式中s为流程坐标,单位m;Q为流量,单位m3/s;Z为水位,单位m;g为重力加速度,单位m/s2;t为时间,单位s;A为过水断面面积,单位m2;R为水力半径,单位m;n为河道糙率系数。
7.根据权利要求6所述的树状河网梯级水库水动力水质泥沙耦合模拟方法,其特征在于,采用Preissmann四点隐式差分格式对一维圣维南方程组进行离散,最后化简为:
离散方程(3)和(4)为编号为i和i+1相邻两断面间的水位和流量关系式,为内断面方程,a1i、c1i、e1i、a2i、c2i、d2i和e2i为隐含有z和Q的系数,需要内部迭代求解;
针对某一内部划分了n个断面的顺直河段,式(3)和(4)第i个矩形网格中建立的2个非线性代数方程式,带有4个未知数 和 n个断面有2n个未知量,可以列2×
(n-1)个方程,加上下游边界条件各一个方程做到闭合,获得的大型稀疏三对方程组采用追赶法求解。
8.根据权利要求1所述的树状河网梯级水库水动力水质及泥沙耦合模拟方法,其特征在于,所述S4中,同步实时模拟和计算包括以下步骤:
S41 读入河网拓扑结构数据文件;
S42 读入河道断面基本信息;
S43 读入初始水位和流量文件;
S44 读入计算时间步长等基本参数;
S45 在每个时间开始迭代前调用边界文件,开始计算水位和流量。
9.一种树状河网梯级水库水动力水质泥沙耦合模拟系统,所述系统用于实现权利要求
1-8所述的任意一项的树状河网梯级水库水动力水质泥沙耦合模拟方法,其特征在于,所述系统的计算软件采用fortran语言编写,包含主程序模、边界文件模块、计算模块、泥沙模块和温度模块。
10.根据权利要求9所述的树状河网梯级水库水动力水质泥沙耦合模拟系统,其特征在于,所述主程序模块运算中,首先读入河网拓扑信息数据和每个河段的断面几何信息数据,并给各断面和水位赋初始值;此时读入每个断面水位与过流面积、水面宽和湿周的对应关系以及每个流量与糙率的对应关系,对此求各断面过流面积、水面、湿周和流速初始值,并给予汇流点的水位赋初始值,最后读入计算需要的其他参数;
所述边界文件模块运算中,首先读入该时刻各节点和分、进水口的流量或水位值,并确定边界该时刻的边界条件;
所述计算模块运算中,首先从总数据库文件中找出上一时刻的数据,并读入该时刻的边界条件,求出各河段的系数,根据系数求出各河段的水位和流量,通过误差判断控制泥沙程序条用,并控制结果输出;
所述泥沙模块运算中,首先从总数据库文件中找出上一时刻的数据,并读入该时刻的边界条件,计算含沙量和河床冲淤高度,对节点赋值并误差判断,读入该时刻各节点的流量或水位值,将节点边界值赋值给各河段上下游断面,最后读入该时刻和分水口含沙量;
所述温度模块运算中,首先从总数据库文件中找出上一时刻的数据,并读入该时刻的边界条件,计算温度场并节点赋值和误差判断,读入该时刻各节点的流量或水位值,将节点边界值赋值给各河段上下游断面,最后读入该时刻和分水口含沙量。

说明书全文

树状河网梯级库水动水质泥沙耦合模拟方法及系统

技术领域

[0001] 本发明涉及水利研究技术领域,具体涉及一种树状河网梯级水库水动力水质泥沙耦合模拟方法及系统。

背景技术

[0002] 国内中上游内陆河网主要为树状形态,干支流河道上建设了大量梯级水库用于防洪、发电、供水、灌溉、航运等目的,科学合理的开展梯级水库的调度运行极为重要。目前已有研究基于实测洪枯季水文资料的洪水特性成果及斯京根法为基础对梯级水库进行调度计算,而传统的一维数学模型的计算由于水库大坝截断河道导致大坝上下游失去水力联系,圣维南方程组不能适用,导致河网不能进行整体连贯的计算。

发明内容

[0003] 针对现有技术的不足,本发明公开了一种树状河网梯级水库水动力水质泥沙耦合模拟方法及系统,用于开发河网水动力及水质、泥沙模型时需要将梯级水库的调度同步考虑进来,实现同步实时模拟和计算,丰富了河网梯级水库调度的研究手段。
[0004] 本发明通过以下技术方案予以实现:
[0005] 一种树状河网梯级水库水动力水质泥沙耦合模拟方法,其特征在于,所述方法包括以下步骤:
[0006] S1根据梯级水库在树状河网中的分布,分割得到河网的控制节点
[0007] S2对河网的控制节点进行分类;
[0008] S3编制程序计算控制节点,实现河网与梯级水库的联合求解;
[0009] S4通过河网整体水动力学模拟计算平台,实现同步实时模拟和计算。
[0010] 优选的,所述S1中,根据梯级水库在树状河网中的分布,将每个水库大坝所在断面作为分割整体河网的控制要素,分隔出若子河网,按照从上游向下游的顺序来编子河网的编号,每个子河网的控制节点也按从上游向下游编排,大坝上游端和下游端各布置一个控制节点且编号连续,作为每个子河网的上边界或下边界控制节点。
[0011] 优选的,所述S2中,将河网的控制节点进行分类,一类是需要外部给定边界条件的节点,称为外边界节点;另一类为内边界节点,无需给定边界条件,是河道交汇口,外边界节点是河道的上、下游端和大坝前端的节点。
[0012] 优选的,所述S3中,编制程序计算的过程中,按照子河网的编号顺序计算,每个子河网按照树状河网的算法独立运算,运算结束后,根据大坝上游端和下游端控制节点间的关系,将该子河网的下游端计算结果通过大坝上、下游端的控制节点传入给下一个子河网,作为下一个子河网的上游端边界条件,依次类推,最后实现河网与梯级水库的联合求解。
[0013] 优选的,所述S4中,模型计算以一维圣维南方程为基础,从河网、子河网、河道、节点四个层面建立模型的计算顺序,对河网及梯级水库进行概化,通过水利枢纽调度方式、外部边界条件的传递关系连接成一个整体,从而将各个河段连接成一个计算整体。
[0014] 优选的,所述一维圣维南方程包括水流连续方程和水流运动方程;
[0015] 所述水流连续方程为:
[0016] 所述水流运动方程为:
[0017]
[0018] 式中s为流程坐标,单位m;Q为流量,单位m3/s;Z为水位,单位m;g为重力加速度,单位m/s2;t为时间,单位s;A为过水断面面积,单位m2;R为水力半径,单位m;n为河道糙率系数。
[0019] 优选的,采用Preissmann四点隐式差分格式对一维圣维南方程组进行离散,最后化简为:
[0020]
[0021]
[0022] 离散方程(3)和(4)为编号为i和i+1相邻两断面间的水位和流量关系式,为内断面方程,a1i、c1i、e1i、a2i、c2i、d2i和e2i为隐含有z和Q的系数,需要内部迭代求解;
[0023] 针对某一内部划分了n个断面的顺直河段,式(3)和(4)第i个矩形网格中建立的2个非线性代数方程式,带有4个未知数 和 n个断面有2n个未知量,可以列2×(n-1)个方程,加上下游边界条件各一个方程便可以做到闭合,获得的大型稀疏三对方程组可以采用追赶法求解。
[0024] 优选的,所述S4中,同步实时模拟和计算包括以下步骤:
[0025] S41读入河网拓扑结构数据文件;
[0026] S42读入河道断面基本信息;
[0027] S43读入初始水位和流量文件;
[0028] S44读入计算时间步长等基本参数;
[0029] S45在每个时间开始迭代前调用边界文件,开始计算水位和流量。
[0030] 本发明还提出一种树状河网梯级水库水动力水质泥沙耦合模拟系统,所述系统用于实现上述的树状河网梯级水库水动力水质泥沙耦合模拟方法,其特征在于,所述系统的计算软件采用fortran语言编写,包含主程序模、边界文件模块、计算模块、泥沙模块和温度模块。
[0031] 优选的,所述主程序模块运算中,首先读入河网拓扑信息数据和每个河段的断面几何信息数据,并给各断面和水位赋初始值;此时读入每个断面水位与过流面积、水面宽和湿周的对应关系以及每个流量与糙率的对应关系,对此求各断面过流面积、水面、湿周和流速初始值,并给予汇流点的水位赋初始值,最后读入计算需要的其他参数;
[0032] 所述边界文件模块运算中,首先读入该时刻各节点和分、进水口的流量或水位值,并确定边界该时刻的边界条件;
[0033] 所述计算模块运算中,首先从总数据库文件中找出上一时刻的数据,并读入该时刻的边界条件,求出各河段的系数,根据系数求出各河段的水位和流量,通过误差判断控制泥沙程序条用,并控制结果输出;
[0034] 所述泥沙模块运算中,首先从总数据库文件中找出上一时刻的数据,并读入该时刻的边界条件,计算含沙量和河床冲淤高度,对节点赋值并误差判断,读入该时刻各节点的流量或水位值,将节点边界值赋值给各河段上下游断面,最后读入该时刻和分水口含沙量;
[0035] 所述温度模块运算中,首先从总数据库文件中找出上一时刻的数据,并读入该时刻的边界条件,计算温度场并节点赋值和误差判断,读入该时刻各节点的流量或水位值,将节点边界值赋值给各河段上下游断面,最后读入该时刻和分水口含沙量。
[0036] 本发明的有益效果为:
[0037] 本发明解决河道地形断面和控制线水文站点资料的完整性和可靠性;加入梯级水库后如何解决河网与梯级水库的联解,依托最新实测干支流河道水深地形资料,建立了将梯级水库调度模式和一维河网水动力模型融合为一体的西江中游河网整体水动力学模拟计算平台。模型计算以一维圣维南方程为基础,提出了从河网、子河网、河道、节点四个层面建立模型的计算顺序,对江中游河网及梯级水库进行概化,通过水利枢纽调度方式、外部边界条件的传递关系连接成一个整体,从而将各个河段连接成一个计算整体,实现同步实时模拟和计算,丰富了江中游洪枯季调度的研究手段。附图说明
[0038] 为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
[0039] 图1本发明树状河网梯级水库水动力水质泥沙耦合模拟方法步骤图;
[0040] 图2是本发明同步实时模拟和计算原理图;
[0041] 图3是本发明主程序模块运算原理图;
[0042] 图4是本发明边界文件模块运算原理图;
[0043] 图5是本发明计算模块运算原理图;
[0044] 图6是本发明泥沙模块运算原理图;
[0045] 图7是本发明温度模块运算原理图;
[0046] 图8是本发明实施例中河网基本计算单元示意图。

具体实施方式

[0047] 为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
[0048] 实施例1
[0049] 本实施例公开了如图1所示的一种树状河网梯级水库水动力水质泥沙耦合模拟方法,包括以下步骤:
[0050] S1根据梯级水库在树状河网中的分布,分割得到河网的控制节点;
[0051] S2对河网的控制节点进行分类;
[0052] S3编制程序计算控制节点,实现河网与梯级水库的联合求解;
[0053] S4通过河网整体水动力学模拟计算平台,实现同步实时模拟和计算。
[0054] S1中,根据梯级水库在树状河网中的分布,将每个水库大坝所在断面作为分割整体河网的控制要素,分隔出若子河网,按照从上游向下游的顺序来编子河网的编号,每个子河网的控制节点也按从上游向下游编排,大坝上游端和下游端各布置一个控制节点且编号连续,作为每个子河网的上边界或下边界控制节点。
[0055] S2中,将河网的控制节点进行分类,一类是需要外部给定边界条件的节点,称为外边界节点;另一类为内边界节点,无需给定边界条件,是河道交汇口,外边界节点是河道的上、下游端和大坝前端的节点。
[0056] S3中,编制程序计算的过程中,按照子河网的编号顺序计算,每个子河网按照树状河网的算法独立运算,运算结束后,根据大坝上游端和下游端控制节点间的关系,将该子河网的下游端计算结果通过大坝上、下游端的控制节点传入给下一个子河网,作为下一个子河网的上游端边界条件,依次类推,最后实现河网与梯级水库的联合求解。
[0057] S4中,模型计算以一维圣维南方程为基础,从河网、子河网、河道、节点四个层面建立模型的计算顺序,对河网及梯级水库进行概化,通过水利枢纽调度方式、外部边界条件的传递关系连接成一个整体,从而将各个河段连接成一个计算整体。
[0058] 本实施例模型计算以一维圣维南方程为基础,提出了从河网、子河网、河道、节点四个层面建立模型的计算顺序,对西江中游河网及梯级水库进行概化,通过水利枢纽调度方式、外部边界条件的传递关系连接成一个整体,从而将各个河段连接成一个计算整体。
[0059] 实施例2
[0060] 本实施例公开一维圣维南方程包括水流连续方程和水流运动方程;
[0061] 所述水流连续方程为:
[0062] 所述水流运动方程为:
[0063]
[0064] 式中s为流程坐标,单位m;Q为流量,单位m3/s;Z为水位,单位m;g为重力加速度,单位m/s2;t为时间,单位s;A为过水断面面积,单位m2;R为水力半径,单位m;n为河道糙率系数。
[0065] 优选的,采用Preissmann四点隐式差分格式对一维圣维南方程组进行离散,最后化简为:
[0066]
[0067]
[0068] 离散方程(3)和(4)为编号为i和i+1相邻两断面间的水位和流量关系式,为内断面方程,a1i、c1i、e1i、a2i、c2i、d2i和e2i为隐含有z和Q的系数,需要内部迭代求解;
[0069] 针对某一内部划分了n个断面的顺直河段,式(3)和(4)第i个矩形网格中建立的2个非线性代数方程式,带有4个未知数 和 n个断面有2n个未知量,可以列2×(n-1)个方程,加上下游边界条件各一个方程便可以做到闭合,获得的大型稀疏三对角方程组可以采用追赶法求解。
[0070] 天然状态下,流域的河系一般呈网状,根据其特征,可分为树状河网和环状河网。树状河网计算遵循的原则是:从支流到干流,从上游到下游。这样可以将河网依次分解成一系列的单一河道,用单一河道的方法求解。树状河网计算原则主要体现在河道编号和断面编号上。
[0071] 如图8所示的最基本的河网计算单元,总共有三条河道,4个河道节点①~④,编号为(一)和(二)的上游河道在节点③汇入河道(三),因此河道的计算顺序是从(一)至(二)到(三);河道(一)内部断面编号从1到N,河道(二)内部断面编号从N+1到K,河道(三)内部断面编号从K+1到M,因此河道断面计算顺序也是依次从1到M,即先计算河道(一)中从断面1到N的流量系数,再计算河道(二)从N+1到K的流量系数,在节点③应用水位和流量平衡原理,获取河道(三)初始断面K+1的系数,以此计算从断面K+2到断面M的系数,最后利用节点④的外部边界条件从断面M依次反推到断面1的水位和流量。
[0072] 实施例3
[0073] 本实施例公开了一种树状河网梯级水库水动力水质泥沙耦合模拟系统,所述系统的计算软件采用fortran语言编写,包含主程序模块、边界文件模块、计算模块、泥沙模块和温度模块。
[0074] 如图3所示的主程序模块运算中,首先读入河网拓扑信息数据和每个河段的断面几何信息数据,并给各断面和水位赋初始值;此时读入每个断面水位与过流面积、水面宽和湿周的对应关系以及每个流量与糙率的对应关系,对此求各断面过流面积、水面、湿周和流速初始值,并给予汇流点的水位赋初始值,最后读入计算需要的其他参数;
[0075] 如图4所示的边界文件模块运算中,首先读入该时刻各节点和分、进水口的流量或水位值,并确定边界该时刻的边界条件;
[0076] 如图5所示的计算模块运算中,首先从总数据库文件中找出上一时刻的数据,并读入该时刻的边界条件,求出各河段的系数,根据系数求出各河段的水位和流量,通过误差判断控制泥沙程序条用,并控制结果输出;
[0077] 如图6所示的泥沙模块运算中,首先从总数据库文件中找出上一时刻的数据,并读入该时刻的边界条件,计算含沙量和河床冲淤高度,对节点赋值并误差判断,读入该时刻各节点的流量或水位值,将节点边界值赋值给各河段上下游断面,最后读入该时刻和分水口含沙量;
[0078] 如图7所示的温度模块运算中,首先从总数据库文件中找出上一时刻的数据,并读入该时刻的边界条件,计算温度场并节点赋值和误差判断,读入该时刻各节点的流量或水位值,将节点边界值赋值给各河段上下游断面,最后读入该时刻和分水口含沙量。
[0079] 本实施例软件开发过程中引入子河网的概念,将水利枢纽作为特殊节点,河网水动力模型可分为河网、子河网、河道、节点四个层面,四者通过水利枢纽调度方式、外部边界条件的传递关系连接成一个整体,从而将各个河段连接成一个计算整体。
[0080] 目前该软件已用于全流域型、中上游型、中下游型洪水下的梯级水库调度、梯级水库调度对梧州断面压咸和生态流量的调度、洪水预报中的洪峰演进研究、突发性水质污染预警预报、西江中游各河段现状洪水水面线研究。
[0081] 以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈