基于水力波动理论的长输供热管道流量参数预测方法

技术领域

本申请涉及数学建模技术及优化算法技术领域，特别是涉及基于水力波动理论的长输供热管道流量参数预测方法、系统、电子设备和存储介质。

背景技术

随着中国经济的稳步发展，居民生活水平显著提升，对于供热的需求大幅度增加，我国的集中供热系统正处于一个前所未有的飞速发展阶段。

在相关技术中，长输供热管网普遍采用质调节、量调节的运行调节方式。但是，在长输供热管道中，流量的变化会引起管道中的压力发生高频震荡，甚至远远超过管网的正常运行压力，因此，可能会造成水锤事故，轻者会造成管道位移、阀门损坏、爆管等，严重时会造成管网瘫痪。所以，对于长输供热管道来说，通过建立水力特性模型进行流量参数分析是十分必要的。

目前针对相关技术中，如何对长输供热管道流量参数进行预测的问题，尚未提出有效的解决方案。

发明内容

本申请实施例提供了一种基于水力波动理论的长输供热管道流量参数预测方法、系统、电子设备和存储介质，以至少解决相关技术中长输供热管道稳定性差的问题。

第一方面，本申请实施例提供了一种基于水力波动理论的长输供热管道流量参数预测方法，所述方法包括：

按照预设管道步长对所述长输供热管道进行节点划分，得到边界节点和内部节点；

基于预先获取的水力常数构建流量参数预测模型，所述水力常数包括管道参数和流体参数；

获取所述边界节点的初始流量参数，通过管网稳态水力-热力耦合模型，根据所述边界节点的初始流量参数以及所述水力常数，获取各个所述内部节点的初始流量参数；

通过所述流量参数预测模型，根据所述初始流量参数和预设时间步长，获取水力波动过程中的不同时刻下，所述边界节点和所述内部节点的流量参数。

在其中一些实施例中，通过所述流量参数预测模型，根据所述初始流量参数和预设时间步长，获取水力波动过程中的不同时刻下，所述边界节点和所述内部节点的流量参数包括：

以所述边界节点和所述内部节点的初始流量参数为迭代初始值，根据所述预设时间步长通过所述流量参数预测模型进行迭代更新，以获取所述内部节点的流量参数；

以所述边界节点的初始流量参数为迭代初始值，根据所述预设时间步长通过边界条件方程和所述流量参数预测模型，以获取所述边界节点的流量参数。

在其中一些实施例中，所述基于预先获取的水力常数构建流量参数预测模型包括：

通过预设有限差分方程，根据所述水力常数，得到流量参数预测系数；

根据所述流量参数预测系数，确定所述流量参数预测模型，其中，所述预设有限差分方程表征了相邻时间步长所述流量参数的映射关系。

在其中一些实施例中，所述方法还包括：

根据所述水力常数，确定所述长输供热管道的水力波动运动方程和水力波动连续性方程；

通过特征线法，对所述水力波动运动方程和所述水力波动连续性方程进行转换，得到所述长输供热管道内流体流速的全微分方程；

对所述全微分方程进行一阶近似变换，得到所述预设有限差分方程。

在其中一些实施例中，所述水力常数包括：流体密度、热扩散系数、管道长度、管道直径、管道摩擦因数、管道压力和管道输热量，基于预先获取的水力常数构建流量参数预测模型包括：

根据预设有限差分方程和所述管道直径，得到第一系数；

根据所述管道摩擦因数和所述管道直径，得到第二系数；

根据所述第一系数和所述第二系数，确定所述流量参数预测模型。

在其中一些实施例中，所述水力常数包括：流体密度、热扩散系数、管道长度、管道直径、管道摩擦因数、管道压力和管道输热量，根据所述水力常数，确定所述长输供热管道的水力波动运动方程和水力波动连续性方程包括：

根据所述管道直径、所述管道摩擦因数和所述管道压力，得到所述水力波动运动方程；

根据所述管道直径、液体体积的弹性模量和应力分布面积，确定压力波波速；

通过所述压力波波速，得到所述水力波动连续性方程。

在其中一些实施例中，所述水力波动运动方程如下：

其中，H是所述管道压力，V是流体流速，f是所述管道摩擦因素，D是所述管道直径，g是重力加速度，x是正方向。

第二方面，本申请实施例提供了一种基于水力波动理论的长输供热管道流量参数预测系统，所述系统包括：节点划分模块、模型构建模块、初始参数确定模块和流量参数确定模块：

所述节点划分模块，用于按照预设管道步长对所述长输供热管道进行节点划分，得到边界节点和内部节点；

所述模型构建模块，用于基于预先获取的水力常数构建流量参数预测模型，所述水力常数包括管道参数和流体参数；

所述初始参数确定模块，用于获取所述边界节点的初始流量参数，通过管网稳态水力-热力耦合模型，根据所述边界节点的初始流量参数以及所述水力常数，获取各个所述内部节点的初始流量参数；

所述流量参数确定模块，用于通过所述流量参数预测模型，根据所述初始流量参数和预设时间步长，获取水力波动过程中的不同时刻下，所述边界节点和所述内部节点的流量参数。

第三方面，本申请实施例提供了一种计算机设备，包括存储器、处理器以及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如上述第一方面所述的一种基于水力波动理论的长输供热管道流量参数预测方法。

第四方面，本申请实施例提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现如上述第一方面所述的一种基于水力波动理论的长输供热管道流量参数预测方法。

相比于相关技术，本申请实施例提供的一种基于水力波动理论的长输供热管道流量参数预测方法，通过预设管道步长对长输供热管道进行节点划分，根据预先获取的水力常数，确定各个管道节点在水力波动前的初始流量参数，通过流量参数预测模型和初始流量参数，得到不同时刻下管道节点的流量参数。解决了如何对长输供热管道流量参数进行预测的问题，通过建立长输供热管道非稳态水力模型，获取长输供热管道各节点的流量参数，精确预测供热管道水力波动过程，为实际工程提供可靠的技术支持。

附图说明

此处所说明的附图用来提供对本申请的进一步理解，构成本申请的一部分，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的不当限定。在附图中：

图1是根据本申请实施例的基于水力波动理论的长输供热管道流量参数预测方法的应用环境示意图；

图2是根据本申请实施例的基于水力波动理论的长输供热管道流量参数预测方法的流程图；

图3是根据本申请实施例的一种长输供热管道微小控制体受力分情况意图；

图4是根据本申请实施例的一种基于特征线法的管道节点水力波动时层推进计算示意图；

图5是根据本申请实施例的一种基于特征线法的描述长输供热管道节点水力波动的有限差分方程计算方法的矩形网格图；

图6是根据本申请实施例的一种长输供热管道水力波动特性分析方法流程图；

图7是根据本申请实施例的基于水力波动理论的长输供热管道流量参数预测系统的结构框图；

图8是根据本申请实施例的电子设备的内部结构示意图。

具体实施方式

为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行描述和说明。应当理解，此处所描述的具体实施例仅仅用以解释本申请，并不用于限定本申请。基于本申请提供的实施例，本领域普通技术人员在没有作出创造性劳动的前提下所获得的所有其他实施例，都属于本申请保护的范围。

显而易见地，下面描述中的附图仅仅是本申请的一些示例或实施例，对于本领域的普通技术人员而言，在不付出创造性劳动的前提下，还可以根据这些附图将本申请应用于其他类似情景。此外，还可以理解的是，虽然这种开发过程中所作出的努力可能是复杂并且冗长的，然而对于与本申请公开的内容相关的本领域的普通技术人员而言，在本申请揭露的技术内容的基础上进行的一些设计，制造或者生产等变更只是常规的技术手段，不应当理解为本申请公开的内容不充分。

在本申请中提及“实施例”意味着，结合实施例描述的特定特征、结构或特性可以包含在本申请的至少一个实施例中。在说明书中的各个位置出现该短语并不一定均是指相同的实施例，也不是与其它实施例互斥的独立的或备选的实施例。本领域普通技术人员显式地和隐式地理解的是，本申请所描述的实施例在不冲突的情况下，可以与其它实施例相结合。

除非另作定义，本申请所涉及的技术术语或者科学术语应当为本申请所属技术领域内具有一般技能的人士所理解的通常意义。本申请所涉及的“一”、“一个”、“一种”、“该”等类似词语并不表示数量限制，可表示单数或复数。本申请所涉及的术语“包括”、“包含”、“具有”以及它们任何变形，意图在于覆盖不排他的包含；例如包含了一系列步骤或模块(单元)的过程、方法、系统、产品或设备没有限定于已列出的步骤或单元，而是可以还包括没有列出的步骤或单元，或可以还包括对于这些过程、方法、产品或设备固有的其它步骤或单元。本申请所涉及的“连接”、“相连”、“耦接”等类似的词语并非限定于物理的或者机械的连接，而是可以包括电气的连接，不管是直接的还是间接的。本申请所涉及的“多个”是指两个或两个以上。“和/或”描述关联对象的关联关系，表示可以存在三种关系，例如，“A和/或B”可以表示：单独存在A，同时存在A和B，单独存在B这三种情况。字符“/”一般表示前后关联对象是一种“或”的关系。本申请所涉及的术语“第一”、“第二”、“第三”等仅仅是区别类似的对象，不代表针对对象的特定排序。

本实施例提供的基于水力波动理论的长输供热管道流量参数预测方法，可以应用于如图1所示的应用环境，图1是根据本申请实施例的基于水力波动理论的长输供热管道流量参数预测方法的应用环境示意图。如图1所示，服务器10按照预设管道步长将长输供热管道划分为边界节点和内部节点，并基于预先获取的水力常数构建流量参数预测模型，服务器10获取边界节点的初始流量参数，根据边界节点的初始流量参数以及水力常数，得到各个内部节点的初始流量参数，通过流量参数预测模型，根据初始流量参数和预设时间步长，获取水力波动过程中的不同时刻下，边界节点和内部节点的流量参数。

本实施例提供了一种基于水力波动理论的长输供热管道流量参数预测方法。图2是根据本申请实施例的基于水力波动理论的长输供热管道流量参数预测方法的流程图，如图2所示，该流程包括如下步骤：

步骤S201，按照预设管道步长对长输供热管道进行节点划分，得到边界节点和内部节点。

本实施例中，将长输供热管道划分为若干个空间步长，得到对应的管道节点，其中，管道节点包括边界节点和内部节点。例如，将长为L的管道划分为N个空间步长，则第1与第N+1个节点为边界节点(也即，管道两端的节点为边界节点)，第2～N个节点为内部节点。

步骤S202，基于预先获取的水力常数构建流量参数预测模型，水力常数包括管道参数和流体参数。

其中，管道参数包括但不限于管道直径、管道长度和摩擦因数，流体参数包括但不限于是流体密度、热扩散系数、管道压力和管道输热量。由于上述常数均可通过本领域常规手段获取，因此，本实施例中对具体手段不再赘述。

步骤S203，获取边界节点的初始流量参数，通过管网稳态水力-热力耦合模型，根据边界节点的初始流量参数以及水力常数，获取各个内部节点的初始流量参数。

初始流量参数即水力波动发生前管道各个节点上的流量参数。边界节点的初始流量参数可以通过测试或远传等手段得到，流量参数包括管道各节点的流量和压力。可选的，通过管网稳态水力-热力耦合模型、起始端压力和流量，以及水力常数，确定水力波动发生前各管道节点的流量和压力。

步骤S204，通过流量参数预测模型，根据初始流量参数和预设时间步长，获取水力波动过程中的不同时刻下，边界节点和内部节点的流量参数。

水力波动主要指的是流体在管道传输过程中所出现的压力波动，这种波动通常是由于供热系统中的流体流动速度、流量变化而产生的。在长输供热管道中，造成水力波动的原因一般有：管道中阀门正常或不正常的开启、关闭和调节以及阀瓣的损坏脱落；水泵正常或事故状态下的开启、停止及调频，以及水泵叶轮在运转时不稳定的振动；管道中含气或被异物堵塞。

本实施例中，将初始流量参数带入流量参数预测模型，获取经过第一个时间步长的时刻下各管道节点的流量参数。进一步的，根据当前时刻的流量参数和流量参数预测模型，得到当前时刻经过一个时间步长后，各个管道节点的流量参数，并不断循环计算，直到管道内部的水力状态重新达到稳定，即水力波动结束。其中，该时间步长可以按照水力波动的传递特征时间与应用场景指定，例如，10s、20s等。

本实施例可以将分析得到的各管道节点的流量参数作为数据支持，对长输供热管道进行调节。可选的，通过调节水泵频率、阀门开度等来改变管道的流量与压力，使长输供热管道在相对平稳的工况下安全运行，避免水力波动对管道的破坏。

通过上述步骤，按照预设管道步长对长输供热管道进行节点划分，根据预先获取的水力常数，确定各个管道节点在水力波动前的初始流量参数，通过流量参数预测模型和初始流量参数，得到不同时刻下管道节点的流量参数。解决了如何对长输供热管道流量参数进行预测的问题，通过建立长输供热管道非稳态水力模型，获取长输供热管道各节点的流量参数，精确预测供热管道水力波动过程，为实际工程提供可靠的技术支持。

在其中一些实施例中，步骤S204具体包括：以边界节点和内部节点的初始流量参数为迭代初始值，根据预设时间步长通过流量参数预测模型进行迭代更新，获取内部节点的流量参数。以边界节点的初始流量参数为迭代初始值，根据预设时间步长通过边界条件方程和流量参数预测模型，以获取边界节点的流量参数。

本实施例中，在经过第一个时间步长的第一时刻，根据水力波动前的初始流量参数和流量参数预测模型，得到各内部节点在水力波动过程中第一时刻的流量参数。并根据第一时刻的水力常数，利用边界条件方程和对应的流量参数预测模型，确定第一时刻边界节点的流量参数。

将第一时刻各管道节点的流量参数作为第二时刻的初始参数，第二时刻为第一时刻经过第二个时间步长的时刻。获取第二时刻的流量参数，具体获取方法与第一时刻的流量参数获取方法相同，此处不再赘述。

后续经过每一个时间步长时刻下的流量参数，都根据前一个时刻的流量参数得到。随着时间步长逐步向前推进，直到水力波动结束，管道内部的水力状态重新达到稳态。

在其中一些实施例中，基于预先获取的水力常数构建流量参数预测模型包括：

通过预设有限差分方程，根据水力常数，得到流量参数预测系数，根据流量参数预测系数，确定流量参数预测模型，其中，预设有限差分方程表征了相邻时间步长流量参数的映射关系。

预设有限差分方程可以是特征线法的有限差分方程，特征线法是一种数值解偏微分方程的方法，有限差分方程是指将偏微分方程中的导数用差分形式进行逼近的方程。本实施例中，应用特征线法在空间上与时间上分别采用特征线进行离散化，采用该方式得到的流量参数预估模型的模型系数具有较好的稳定性和精度。进一步的，流量参数预测模型可以是相容性方程，该方程和其他条件结合使用，通过数值模拟方法来预测和估计后续的流量参数。通过本实施例得到的流量参数预测模型，可以精确预估不同时刻下长输供热管道各节点的流量参数。

在其中一些实施例中，预设有限差分方程的获取方法包括：根据水力常数，确定长输供热管道的水力波动运动方程和水力波动连续性方程；通过特征线法，对水力波动运动方程和水力波动连续性方程进行转换，得到长输供热管道内流体流速的全微分方程；对全微分方程进行一阶近似变换，得到预设有限差分方程。

在其中一些实施例中，步骤S202具体包括：根据预设有限差分方程和管道直径，得到第一系数，根据管道摩擦因数和管道直径，得到第二系数，根据第一系数和第二系数，确定流量参数预测模型。

本实施例中，当流量参数预测模型是相容性方程时，利用特征线法的有限差分方程，结合重力加速度、压力波波速和管道横截面积，计算相容性方程的第一常数系数，结合摩擦因数、计算步长、重力加速度、管道直径和管道横截面积，计算相容性方程的第二系数。

在管道发生水力波动时，压力通常是以秒为单位进行往复周期波动的。为了进行长输供热管道水力波动分析，需要根据弹性水柱理论建立运动方程与连续性方程。选取控制体进行受力分析，图3是根据本申请实施例的一种长输供热管道微小控制体受力情况示意图，如图3所示，管道流体受到前后表面流体压力、流体可压缩性产生的管壁压力、管壁摩擦力和重力的作用。依据牛顿第二定律得到公式一，

公式一：

其中，x表示正方向，A表示管道面积，P

非稳态流动中，与控制体压力变化相比，密度变化很小，可以忽略。将公式一展开，忽略高阶微量，得到公式二，

公式二：

非稳态流动时的切应力可采用相同流速时稳态流动的切应力表达式，

切应力表达式：

其中，V是流体流速，f是管道摩擦因素，在上述切应力表达式中引入绝对值符号，可以确保流体所受的摩擦力总与流速反向。

根据达西公式

公式三：

控制体的流速V是时间t和空间x的函数，故得到公式四，

公式四：

将切应力表达式和公式四带入公式二中得到公式五，

公式五：

需要说明的是，公式五同样适用于分流或合流的流体

因为与压力H或控制体标高z相比，密度ρ的变化非常小，故可当做常量。对公式五进行化简，得到长输供热管道的水力波动运动方程，

水力波动运动方程：

其中，H是管道压力，D是管道直径，V是流体流速，f是管道摩擦因素，g是重力加速度，x是正方向。p和H均表示压力，但单位不同，p的单位是帕或千帕，H的单位是米水柱。

根据质量守恒定律，单位时间内流入控制体的净质量，应等于单位时间内控制体由于升压引起管道膨胀和流体被压缩而增加的流体质量，则有公式六，

公式六：

其中，u表示管壁在x处的速度。

将公式六展开化简，得到公式七，

公式七：

AVρ

ρ和A的全微分如公式八，

公式八：

它们是沿着流体运动的全微分，将公式八带入式公式七，得到公式九，

公式九：

在公式九中，第一项代表流体的可压缩性，第二项代表管壁材料的弹性，第三项代表单位时间内流入控制体的质量与流出控制体的质量之差。

本实施例中，压力波波速a的计算公式如公式十，

公式十：

其中，K为液体体积的弹性模量，

将公式十沿流体运动方向进行全微分展开，忽略轴向变形，引入液体体积的弹性模量，并带入公式九，即可得到长输供热管道的水力波动连续性方程，

水力波动连续性方程：

在实际工程中，管道多采用金属材料，其压力波波速a一般为1000m/s，远远大于管道内的流速，因此，可以忽略影响较小的流速项，同时忽略因高度差产生的压强变化，则上述水力波动运动方程和水力波动连续性方程可以简化为：

为了便于模型求解，采用特征线法将上述偏微分方程转变为全微分方程，并将这两个方程式进行线性组合，得到公式十一，

公式十一：

其中，λ表示待定实数。

当取±g/a时，公式十一可转化为全微分方程组：

上述全微分方程组的不同形式，分别代表了当前时刻下相邻节点在正向压力波和反向压力波的作用下，沿斜率为1/a与-1/a的特征线经过一定时间步长到达同一流体状态点，如图4所示，图4是根据本申请实施例的一种基于特征线法的管道节点水力波动时层推进计算示意图。通过该方法，对所有管道节点进行分析，如图5所示，图5是根据本申请实施例的一种基于特征线法的描述长输供热管道节点水力波动的有限差分方程计算方法的矩形网格图。

将全微分方程组沿特征线进行积分，且用流量Q表示流速V，并进行一阶近似，认为在一个空间步长内流量不变，可得到特征线法的第一有限差分方程和第二有限差分方程，

第一有限差分方程：

第二有限差分方程：

其中，Q表示管道流量，单位m

在此基础上，可根据x-t平面坐标系中的特征网格描述有限差分的计算过程。从初始状态开始，只要初始时刻相邻节点的流量与压力已知，则可计算下一时刻特征网格节点的流量与压力，结合管道内阀门、水泵及关键节点的边界条件，可以得到每一时刻所有节点上的流量和压力。

图6是根据本申请实施例的一种长输供热管道水力波动特性分析方法流程图，如图6所示，该流程图包括如下步骤：

步骤S61，输入重力加速度g、流体密度σ、热扩散系数α、管道长度L、管道直径D、管道摩擦因数f、管道节点数N、水力波动空间长度t

步骤S62，根据水力常数，得到相容性方程的第一常数系数B和第二常数系数R；

步骤S63，根据水力常数，计算初始时刻t＝0、计数值Z＝0时各管道节点的流量Q

步骤S64，时间步长加一，判断当前时间t是否大于水力波动时间他们，若是，结束水力波动分析，若否，计数值Z加一，并利用相容性方程确定该时刻内部节点i＝2至i＝N处的流量和压力；

步骤S65，重复步骤S64，直到当前时间t大于水力波动时间。

需要说明的是，在上述流程中或者附图的流程图中示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。

本实施例还提供了一种基于水力波动理论的长输供热管道调节系统，该装置用于实现上述实施例及优选实施方式，已经进行过说明的不再赘述。如以下所使用的，术语“模块”、“单元”、“子单元”等可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。

图7是根据本申请实施例的基于水力波动理论的长输供热管道流量参数预测系统，如图7所示，该系统包括：节点划分模块71、模型构建模块72、初始参数确定模块73和流量参数确定模块74：

节点划分模块71，用于按照预设管道步长对长输供热管道进行节点划分，得到边界节点和内部节点；

模型构建模块72，用于基于预先获取的水力常数构建流量参数预测模型，水力常数包括管道参数和流体参数；

初始参数确定模块73，用于获取边界节点的初始流量参数，通过管网稳态水力-热力耦合模型，根据边界节点的初始流量参数以及水力常数，获取各个内部节点的初始流量参数；

流量参数确定模块74，用于通过流量参数预测模型，根据初始流量参数和预设时间步长，获取水力波动过程中的不同时刻下，边界节点和内部节点的流量参数。

本实施例通过节点划分模块71按照预设管道步长对长输供热管道进行节点划分，得到边界节点和内部节点，模型构建模块72基于预先获取的水力常数构建流量参数预测模型，初始参数确定模块73获取边界节点的初始流量参数，通过管网稳态水力-热力耦合模型，根据边界节点的初始流量参数以及水力常数，获取各个内部节点的初始流量参数，流量参数确定模块74通过流量参数预测模型，根据初始流量参数和预设时间步长，获取水力波动过程中的不同时刻下边界节点和内部节点的流量参数。解决了如何对长输供热管道流量参数进行预测的问题，通过建立长输供热管道非稳态水力模型，获取长输供热管道各节点的流量参数，精确预测供热管道水力波动过程，为实际工程提供可靠的技术支持。

需要说明的是，上述各个模块可以是功能模块也可以是程序模块，既可以通过软件来实现，也可以通过硬件来实现。对于通过硬件来实现的模块而言，上述各个模块可以位于同一处理器中；或者上述各个模块还可以按照任意组合的形式分别位于不同的处理器中。

本实施例还提供了一种电子设备，包括存储器和处理器，该存储器中存储有计算机程序，该处理器被设置为运行计算机程序以执行上述任一项方法实施例中的步骤。

可选地，上述电子设备还可以包括传输设备以及输入输出设备，其中，该传输设备和上述处理器连接，该输入输出设备和上述处理器连接。

可选地，在本实施例中，上述处理器可以被设置为通过计算机程序执行以下步骤：

S1，按照预设管道步长对长输供热管道进行节点划分，得到边界节点和内部节点。

S2，基于预先获取的水力常数构建流量参数预测模型，水力常数包括管道参数和流体参数。

S3，获取边界节点的初始流量参数，通过管网稳态水力-热力耦合模型，根据边界节点的初始流量参数以及水力常数，获取各个内部节点的初始流量参数。

S4，通过流量参数预测模型，根据初始流量参数和预设时间步长，获取水力波动过程中的不同时刻下，边界节点和内部节点的流量参数。

需要说明的是，本实施例中的具体示例可以参考上述实施例及可选实施方式中所描述的示例，本实施例在此不再赘述。

在一个实施例中，图8是根据本申请实施例的电子设备的内部结构示意图，如图8所示，提供了一种电子设备，该电子设备可以是服务器，其内部结构图可以如图8所示。该电子设备包括通过系统总线连接的处理器、存储器、网络接口和数据库。其中，该电子设备的处理器用于提供计算和控制能力。该电子设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统、计算机程序和数据库。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该电子设备的数据库用于存储数据。该电子设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种基于水力波动理论的长输供热管道流量参数预测方法。

本领域技术人员可以理解，图8中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的电子设备的限定，具体的电子设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，该计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限，RAM以多种形式可得，诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。

本领域的技术人员应该明白，以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。

以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请专利的保护范围应以所附权利要求为准。