CN111553041B - 一种燃气管网参数校验方法、系统、终端及存储介质 - Google Patents
一种燃气管网参数校验方法、系统、终端及存储介质 Download PDFInfo
- Publication number
- CN111553041B CN111553041B CN202010326984.8A CN202010326984A CN111553041B CN 111553041 B CN111553041 B CN 111553041B CN 202010326984 A CN202010326984 A CN 202010326984A CN 111553041 B CN111553041 B CN 111553041B
- Authority
- CN
- China
- Prior art keywords
- parameters
- verified
- resistance loss
- theta
- measured data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000003860 storage Methods 0.000 title claims abstract description 14
- 238000004364 calculation method Methods 0.000 claims abstract description 34
- 238000002940 Newton-Raphson method Methods 0.000 claims abstract description 16
- 238000013461 design Methods 0.000 claims abstract description 8
- 239000007789 gas Substances 0.000 claims description 102
- 239000011159 matrix material Substances 0.000 claims description 30
- 239000002737 fuel gas Substances 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 3
- 238000004088 simulation Methods 0.000 abstract description 21
- 238000005259 measurement Methods 0.000 abstract description 10
- 238000005457 optimization Methods 0.000 description 10
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000018109 developmental process Effects 0.000 description 3
- 239000012530 fluid Substances 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000013480 data collection Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013524 data verification Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000006116 polymerization reaction Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Networks & Wireless Communication (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本申请提供一种燃气管网参数校验方法、系统、终端及存储介质,所述方法包括:根据经验公式确定燃气管道摩擦阻力损失计算公式;确定燃气管道摩擦阻力损失计算公式的待校验参数及实测数据;根据管道设计图纸和经验设定待校验参数的初始值;获取与待校验参数数量相等的实测数据对;利用所述实测数据对组成的实测数据建立燃气管道摩擦阻力损失计算公式的非线性方程组,并采用牛顿‑拉夫逊法求解非线性方程组得到多组参数校验值;获取多组参数校验值并计算参数校验值的均值,得到最终的参数的估计值;本申请基于燃气网管SCADA采集数据系统和燃气管道实际参数相结合对比校验的方法,并通过牛顿‑拉夫逊方法进行求解分析,实现真实工况参数实时模拟测量。
Description
技术领域
本申请涉及数据处理技术领域,尤其是涉及一种燃气管网参数校验方法、系统、终端及存储介质。
背景技术
随着我国城镇化水平的提升,城市规模不断扩大,为城市提供服务的基础设施也越来越多。天然气作为清洁能源,已成为城市运行和生产、生活的主体能源。随着城市建设发展,城市燃气管网的规模逐渐增大,拓扑结构趋向多样化、复杂化,同时一些老旧管网建设年代较长,数据采集难度大;城市建设的快速发展,部分管网更新改造变动多,仅根据人工经验和后期部署的SCADA系统数据来分析管网拓扑结构,很大程度上不能满足使用需要,同时也为燃气管网的实时调度和智能化系统升级带来一定障碍。
随着信息技术的发展和推广应用,国内天然气管网基础设施的自动化、信息化程度逐步提高,相关信息系统也逐渐建立与完善,使得借助实测数据和机理模型相结合,推断校验燃气管网拓扑结构准确性的相关方法应用成为可能,并且能够为燃气管网调度提供可靠的数据支撑。目前,国内关于燃气管网拓扑结构确定和参数校验问题的解决方案大多是基于宏观水力模型展开的,由于宏观水力模型是以数据驱动的,以此为基础的调度方案必定存在缺陷,例如无法得到非测点的运行状态,遇到改扩建、检修等引起管网拓扑结构改变的情况时,原调度方法会因模型失准而失效。
例如基于数据实时性采集而提出的应用于自来水管网调度的“基于在线构架的供水实时调度方法”(ZL201510078366.5),该发明专利是直接利用实时水量数据开展供水调度,虽在一定程度上可以提高响应速度,但由于缺乏供水管网实时模型的支撑,无法全面了解供水管网的实时工况,必然会影响调度决策效果。也有类似用于供热管网的“一种供热管网等效建模方法”(ZL201611164890.5),是按流体网络建模技术建立集中供热一次管网系统数学模型,以及建立的数学模型的等效聚合,得到模型参数和校验模型参数与实测数据是否存在偏差,但是对于模型数据的校验没有涉及,当数据出现问题时,容易产生误解偏差。
现有的燃气管网仿真与优化方法中,在进行仿真优化前,需要预先指定管网参数,包括却不仅限于管道长度、管道内径、管道摩擦因子等,但所使用的参数多为经验值或参考值、标称值。而实际工程中,一方面,工程实施误差和管道老化等因素会导致上述预设参数与实际有较大出入,由于缺少对于实际参数的采集和校正会导致仿真优化结果与实际工况不符合的情况,使得最终的仿真优化结果失真率高,不具有参考意义,对实际应用和策略决策造成影响;另一方面,现场实际工况数据和参数的采集和核验往往耗时耗力,同时城市燃气管网数量众多,数据采集点不具有可持续性,所以对于燃气管网参数的实时模拟测量和更新显得尤为重要。
有鉴于此,在燃气管网全面数字化背景下,亟需一种燃气管网参数校验方法、系统、终端及存储介质,以解决现有技术中燃气管网参数与实际工况偏差较大、实际工况数据和参数的采集和核验实时性差导致仿真优化结果失真率高的问题。
发明内容
针对现有技术的不足,本申请提供一种燃气管网参数校验方法、系统、终端及存储介质,解决了现有技术中燃气管网参数与实际工况偏差较大、实际工况数据和参数的采集和核验实时性差导致仿真优化结果失真率高的问题。
为解决上述技术问题,第一方面,本申请提供一种燃气管网参数校验方法,包括:
根据经验公式确定燃气管道摩擦阻力损失计算公式;
确定燃气管道摩擦阻力损失计算公式的待校验参数及实测数据;
根据管道设计图纸和经验设定待校验参数的初始值;
获取与待校验参数数量相等的实测数据对;
利用所述实测数据对组成的实测数据建立燃气管道摩擦阻力损失计算公式的非线性方程组,并采用牛顿-拉夫逊法求解所述非线性方程组得到多组参数校验值;
获取多组参数校验值并计算参数校验值的均值,得到最终的参数的估计值。
可选的,所述利用所述实测数据对组成的实测数据建立燃气管道摩擦阻力损失计算公式的非线性方程组,并采用牛顿-拉夫逊法求解所述非线性方程组,包括:
根据待校验参数的初始值初始化所述待校验参数;
计算燃气管道摩擦阻力损失数值矩阵F与雅各比矩阵J;
判断燃气管道摩擦阻力损失数值是否满足预设精度;若是,则返回得到当前待校验参数;若否,则更新待校验参数,直到燃气管道摩擦阻力损失数值满足预设精度。
可选的,所述计算燃气管道摩擦阻力损失数值矩阵F与雅各比矩阵J,包括:
设定燃气管道摩擦阻力损失数值公式为f(θ),待校验参数个数为i;
根据所述实测数据对组成的实测数据与待校验参数的当前值计算所述燃气管道摩擦阻力损失数值即fi(θ)数值;
根据各组燃气管道摩擦阻力损失数值即fi(θ)数值得到燃气管道摩擦阻力损失数值矩阵F;
根据公式f(θ)计算得到雅各比矩阵J。
可选的,所述所述判断燃气管道摩擦阻力损失数值是否满足预设精度;若是,则返回得到当前待校验参数;若否,则更新待校验参数,直到燃气管道摩擦阻力损失数值满足预设精度,包括:
判断max(|fi|)<ε是否成立;
若是,则返回当前θ;
若否,则根据公式θ′=θ-J-1F计算得到更新后的待校验参数,直到max(|fi|)<ε成立;
其中,ε为预设精度,例如可取1e-5。
第二方面,本申请还提供一种燃气管网参数校验系统,包括:
第一确定单元,配置用于根据经验公式确定燃气管道摩擦阻力损失计算公式;
第二确定单元,配置用于确定燃气管道摩擦阻力损失计算公式的待校验参数及实测数据;
初始设定单元,配置用于根据管道设计图纸和经验设定待校验参数的初始值;
数据获取单元,配置用于获取与待校验参数数量相等的实测数据对;
参数计算单元,配置用于利用所述实测数据对组成的实测数据建立燃气管道摩擦阻力损失计算公式的非线性方程组,并采用牛顿-拉夫逊法求解所述非线性方程组得到多组参数校验值;
参数估值单元,配置用于获取多组参数校验值并计算参数校验值的均值,得到最终的参数的估计值。
可选的,所述参数计算单元具体包括:
初始赋值单元,配置用于根据待校验参数的初始值初始化所述待校验参数;
矩阵计算单元,配置用于计算燃气管道摩擦阻力损失数值矩阵F与雅各比矩阵J;
参数更新单元,配置用于判断燃气管道摩擦阻力损失数值是否满足预设精度;若是,则返回得到当前待校验参数;若否,则更新待校验参数,直到燃气管道摩擦阻力损失数值满足预设精度。
可选的,所述矩阵计算单元具体用于:
设定燃气管道摩擦阻力损失数值公式为f(θ),待校验参数个数为i;
根据所述实测数据对组成的实测数据与待校验参数的当前值计算所述燃气管道摩擦阻力损失数值即fi(θ)数值;
根据各组燃气管道摩擦阻力损失数值即fi(θ)数值得到燃气管道摩擦阻力损失数值矩阵F;
根据公式f(θ)计算得到雅各比矩阵J。
可选的,所述参数更新单元具体用于:
判断max(|fi|)<ε是否成立;
若是,则返回当前θ;
若否,则根据公式θ′=θ-J-1F计算得到更新后的待校验参数,直到max(|fi|)<ε成立;
其中,ε为预设精度,例如可取1e-5。
第三方面,本申请提供一种终端,包括:
处理器、存储器,其中,
该存储器用于存储计算机程序,
该处理器用于从存储器中调用并运行该计算机程序,使得终端执行上述的终端的方法。
第四方面,本申请提供了一种计算机存储介质,所述计算机可读存储介质中存储有指令,当其在计算机上运行时,使得计算机执行上述各方面所述的方法。
与现有技术相比,本申请具有如下有益效果:
1、本申请基于燃气网管SCADA采集数据系统和燃气管道实际参数相结合对比校验的方法,将SCADA系统监测点收集的实时数据作为实际参数的体现,并结合机理模型通过牛顿-拉夫逊方法进行求解分析,最终实现真实工况参数的实时模拟测量,确保后续的管网仿真与优化是建立在较为准确的参数之上。
2、本申请通过将关键的压降方程转化为非线性方程组,从而可以使用牛顿-拉夫逊方法进行求解,求解速度快,其对于非线性方程组的求解而言,非常迅速。这对于大规模的燃气管网的实时模拟参数更新有重大作用。
3、本申请根据传热学、流体力学原理和相关经验公式等,建立符合实际燃气管网物理特性的相关计算方程;通过关键参数辨识和现场实测参数进行数据对比,可以实现大型复杂燃气管网的参数简化与校核,提高仿真和优化调度所需数据的可靠性和精确度。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1为本申请实施例所提供的一种燃气管网参数校验方法的流程图;
图2为本申请实施例所提供的另一种燃气管网参数校验方法的流程图;
图3为本申请另一实施例所提供的一种燃气管网参数校验系统的结构示意图;
图4为本申请实施例所提供的一种终端系统的结构示意图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
请参考图1,图1为本申请实施例所提供的一种燃气管网参数校验方法的流程图,该方法100包括:
S101:根据经验公式确定燃气管道摩擦阻力损失计算公式;
S102:确定燃气管道摩擦阻力损失计算公式的待校验参数及实测数据;
S103:根据管道设计图纸和经验设定待校验参数的初始值;
S104:获取与待校验参数数量相等的实测数据对;
S105:利用所述实测数据对组成的实测数据建立燃气管道摩擦阻力损失计算公式的非线性方程组,并采用牛顿-拉夫逊法求解所述非线性方程组得到多组参数校验值;
S106:获取多组参数校验值并计算参数校验值的均值,得到最终的参数的估计值。
基于上述实施例,作为可选的实施例,如图3所示,所述S105利用所述实测数据对组成的实测数据建立燃气管道摩擦阻力损失计算公式的非线性方程组,并采用牛顿-拉夫逊法求解所述非线性方程组,包括:
根据待校验参数的初始值初始化所述待校验参数;
计算燃气管道摩擦阻力损失数值矩阵F与雅各比矩阵J;
判断燃气管道摩擦阻力损失数值是否满足预设精度;若是,则返回得到当前待校验参数;若否,则更新待校验参数,直到燃气管道摩擦阻力损失数值满足预设精度。
基于上述实施例,作为可选的实施例,所述计算燃气管道摩擦阻力损失数值矩阵F与雅各比矩阵J,包括:
设定燃气管道摩擦阻力损失数值公式为f(θ),待校验参数个数为i;
根据所述实测数据对组成的实测数据与待校验参数的当前值计算所述燃气管道摩擦阻力损失数值即fi(θ)数值;
根据各组燃气管道摩擦阻力损失数值即fi(θ)数值得到燃气管道摩擦阻力损失数值矩阵F;
根据公式f(θ)计算得到雅各比矩阵J。
基于上述实施例,作为可选的实施例,所述判断燃气管道摩擦阻力损失数值是否满足预设精度;若是,则返回得到当前待校验参数;若否,则更新待校验参数,直到燃气管道摩擦阻力损失数值满足预设精度,包括:
判断max(|fi|)<ε是否成立;
若是,则返回当前θ;
若否,则根据公式θ′=θ-J-1F计算得到更新后的待校验参数,直到max(|fi|)<ε成立;
其中,ε为预设精度,例如可取1e-5。
具体的,本申请可通过SCADA系统采集到的燃气管网节点处的燃气压强、温度、流速数据,推导计算出管段的长度、内径、管道内壁的当量绝对粗糙度等实际工况参数。本申请将以上述参数为例进行说明,但本申请专利不仅限于计算这些参数。具体过程如下:
步骤1、根据经验公式确定燃气管道摩擦阻力损失计算公式(公式1),例如钢管:
步骤2、确定上述的燃气管道摩擦阻力损失计算公式的待校验参数L、△、d、ρ0四个参数,确定上述的燃气管道摩擦阻力损失计算公式的实测数据p1、p2、Q0、T;
步骤3、将步骤1的公式1设定为f(θ)(公式2),并计算f(θ)对θ的偏导数(公式3):
步骤4、根据管道设计图纸和经验,设定L、△、d、ρ0的初始值;
步骤5、获取与待校验参数数量相等的实测数据对,即获取4组实测数据对(p1,p2,Q0,T);
步骤6、利用所述步骤5获取的4组实测数据对,建立关于fi(θ)=0的非线性方程组,每个i对应一组实测数据对,本例中包含4个非线性方程,采用牛顿-拉夫逊法求解所述非线性方程组。
步骤6.1、将所述步骤4的待校验参数的初始值赋值到待校验参数的当前值;
步骤6.2、根据所述步骤5获取的实测数据与待校验参数的当前值计算所述步骤6中的各个fi(θ)数值,根据公式4得到燃气管道摩擦阻力损失数值矩阵F,并根据公式5计算当前的雅各比矩阵:
步骤6.3、若max(|fi|)<ε,即可终止步骤6,并返回当前θ,否则继续进行步骤6.4。其中ε为预设精度,例如可取1e-5;
步骤6.4、根据公式6计算一次更新后的待校验参数θ′,并重复进行步骤6.2。
θ′=θ-J-1F (公式6)
步骤7、重复多次步骤6,获取多组参数θ的校验值,取均值为最终的参数的估计值。
需要说明的是,针对已有的实际参数测量问题,创新性的将实测数据与参数校验相结合,通过节点的压强、流量和温度等数据,来辅助求解管网的工况参数,从而加强仿真优化的精度。将关键压降方程转化为可以用牛顿-拉夫逊方法求解的参数非线性方程组,并结合实际数据进行参数求解。牛顿-拉夫逊方法对于非线性方程组的求解非常迅速,从而能够实现实际工况参数的超实时模拟测量,也为超实时仿真和优化打下了基础。
请参考图3,图3为本申请实施例所提供的一种燃气管网参数校验系统的结构示意图,该系统300,包括:
第一确定单元301,配置用于根据经验公式确定燃气管道摩擦阻力损失计算公式;
第二确定单元302,配置用于确定燃气管道摩擦阻力损失计算公式的待校验参数及实测数据;
初始设定单元303,配置用于根据管道设计图纸和经验设定待校验参数的初始值;
数据获取单元303,配置用于获取与待校验参数数量相等的实测数据对;
参数计算单元305,配置用于利用所述实测数据对组成的实测数据建立燃气管道摩擦阻力损失计算公式的非线性方程组,并采用牛顿-拉夫逊法求解所述非线性方程组得到多组参数校验值;
参数估值单元306,配置用于获取多组参数校验值并计算参数校验值的均值,得到最终的参数的估计值。
基于上述实施例,作为可选的实施例,所述参数计算单元305具体包括:
初始赋值单元,配置用于根据待校验参数的初始值初始化所述待校验参数;;
矩阵计算单元,配置用于计算燃气管道摩擦阻力损失数值矩阵F与雅各比矩阵J;
参数更新单元,配置用于判断燃气管道摩擦阻力损失数值是否满足预设精度;若是,则返回得到当前待校验参数;若否,则更新待校验参数,直到燃气管道摩擦阻力损失数值满足预设精度。
基于上述实施例,作为可选的实施例,所述矩阵计算单元具体用于:
设定燃气管道摩擦阻力损失数值公式为f(θ),待校验参数个数为i;
根据所述实测数据对组成的实测数据与待校验参数的当前值计算所述燃气管道摩擦阻力损失数值即fi(θ)数值;
根据各组燃气管道摩擦阻力损失数值即fi(θ)数值得到燃气管道摩擦阻力损失数值矩阵F;
根据公式f(θ)计算得到雅各比矩阵J。
基于上述实施例,作为可选的实施例,所述参数更新单元具体用于:
判断max(|fi|)<ε是否成立;
若是,则返回当前θ;
若否,则根据公式θ′=θ-J-1F计算得到更新后的待校验参数,直到max(|fi|)<ε成立;
其中,ε为预设精度,例如可取1e-5。
请参考图4,图4为本申请实施例所提供的一种终端系统400的结构示意图,该终端系统400可以用于执行本发明实施例提供的软件多语言显示及输入同步切换方法。
其中,该终端系统400可以包括:处理器401、存储器402及通信单元403。这些组件通过一条或多条总线进行通信,本领域技术人员可以理解,图中示出的服务器的结构并不构成对本发明的限定,它既可以是总线形结构,也可以是星型结构,还可以包括比图示更多或更少的部件,或者组合某些部件,或者不同的部件布置。
其中,该存储器402可以用于存储处理器401的执行指令,存储器402可以由任何类型的易失性或非易失性存储终端或者它们的组合实现,如静态随机存取存储器(SRAM),电可擦除可编程只读存储器(EEPROM),可擦除可编程只读存储器(EPROM),可编程只读存储器(PROM),只读存储器(ROM),磁存储器,快闪存储器,磁盘或光盘。当存储器402中的执行指令由处理器401执行时,使得终端系统400能够执行以下上述方法实施例中的部分或全部步骤。
处理器401为存储终端的控制中心,利用各种接口和线路连接整个电子终端的各个部分,通过运行或执行存储在存储器402内的软件程序和/或模块,以及调用存储在存储器内的数据,以执行电子终端的各种功能和/或处理数据。所述处理器可以由集成电路(Integrated Circuit,简称IC)组成,例如可以由单颗封装的IC所组成,也可以由连接多颗相同功能或不同功能的封装IC而组成。举例来说,处理器401可以仅包括中央处理器(Central Processing Unit,简称CPU)。在本发明实施方式中,CPU可以是单运算核心,也可以包括多运算核心。
通信单元403,用于建立通信信道,从而使所述存储终端可以与其它终端进行通信。接收其他终端发送的用户数据或者向其他终端发送用户数据。
本申请还提供一种计算机存储介质,其中,该计算机存储介质可存储有程序,该程序执行时可包括本发明提供的各实施例中的部分或全部步骤。所述的存储介质可为磁碟、光盘、只读存储记忆体(英文:read-only memory,简称:ROM)或随机存储记忆体(英文:random access memory,简称:RAM)等。
本申请基于燃气网管SCADA采集数据系统和燃气管道实际参数相结合对比校验的方法,将SCADA系统监测点收集的实时数据作为实际参数的体现,并结合机理模型通过牛顿-拉夫逊方法进行求解分析,最终实现真实工况参数的实时模拟测量,确保后续的管网仿真与优化是建立在较为准确的参数之上。本申请通过将关键的压降方程转化为非线性方程组,从而可以使用牛顿-拉夫逊方法进行求解,求解速度快,其对于非线性方程组的求解而言,非常迅速。这对于大规模的燃气管网的实时模拟参数更新有重大作用。本申请根据传热学、流体力学原理和相关经验公式等,建立符合实际燃气管网物理特性的相关计算方程;通过关键参数辨识和现场实测参数进行数据对比,可以实现大型复杂燃气管网的参数简化与校核,提高仿真和优化调度所需数据的可靠性和精确度。
说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例提供的系统而言,由于其与实施例提供的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本申请的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本申请的方法及其核心思想。应当指出,对于本技术领域的普通技术人员来说,在不脱离本申请原理的前提下,还可以对本申请进行若干改进和修饰,这些改进和修饰也落入本申请权利要求的保护范围内。
还需要说明的是,在本说明书中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
Claims (4)
1.一种燃气管网参数校验方法,其特征在于,包括:
根据经验确定燃气管道摩擦阻力损失计算公式1;
确定燃气管道摩擦阻力损失计算公式的待校验参数L、Δ、d、ρ0及实测数据p1、p2、Q0、T;
根据管道设计图纸和经验设定待校验参数L、Δ、d、ρ0的初始值;
获取与待校验参数数量相等的实测数据对(p1,p2,Q0,T);
根据待校验参数的初始值初始化所述待校验参数;
设定燃气管道摩擦阻力损失数值公式1为f(θ),其中θ=[L Δ d ρ0]T;
建立关于fi(θ)=0的非线性方程组,待校验参数个数为i,每个i对应一组实测数据对,根据所述获取的实测数据与待校验参数的当前值计算每个i对应的fi(θ)数值;
根据所述实测数据对组成的实测数据与待校验参数的当前值计算所述燃气管道摩擦阻力损失数值即fi(θ)数值;
根据各组燃气管道摩擦阻力损失数值即fi(θ)数值得到燃气管道摩擦阻力损失数值矩阵F;
根据各组燃气管道摩擦阻力损失数值即fi(θ)数值计算得到雅各比矩阵J;
判断max(|fi(θ)|)<ε是否成立;
若是,则返回当前θ;
若否,则根据公式θ′=θ-J-1F计算得到更新后的待校验参数,直到max(|fi(θ)|)<ε成立;
其中,ε为预设精度,可取1e-5。
2.一种燃气管网参数校验系统,其特征在于,包括:
第一确定单元,配置用于根据经验公式确定燃气管道摩擦阻力损失计算公式1,
第二确定单元,配置用于确定燃气管道摩擦阻力损失计算公式的待校验参数L、Δ、d、ρ0及实测数据p1、p2、Q0、T;
初始设定单元,配置用于根据管道设计图纸和经验设定待校验参数L、Δ、d、ρ0的初始值;
数据获取单元,配置用于获取与待校验参数数量相等的实测数据对(p1,p2,Q0,T);
参数计算单元,配置用于利用所述实测数据对组成的实测数据建立燃气管道摩擦阻力损失计算公式的非线性方程组,并采用牛顿-拉夫逊法求解所述非线性方程组得到多组参数校验值;
参数估值单元,配置用于根据待校验参数的初始值初始化所述待校验参数;
设定燃气管道摩擦阻力损失数值公式1为f(θ),其中θ=[L Δ d ρ0]T;
建立关于fi(θ)=0的非线性方程组,待校验参数个数为i,每个i对应一组实测数据对,根据所述获取的实测数据与待校验参数的当前值计算每个i对应的fi(θ)数值;
根据所述实测数据对组成的实测数据与待校验参数的当前值计算所述燃气管道摩擦阻力损失数值即fi(θ)数值;
根据各组燃气管道摩擦阻力损失数值即fi(θ)数值得到燃气管道摩擦阻力损失数值矩阵F;
根据各组燃气管道摩擦阻力损失数值即fi(θ)数值计算得到雅各比矩阵J;
判断max(|fi(θ)|)<ε是否成立;
若是,则返回当前θ;
若否,则根据公式θ′=θ-J-1F计算得到更新后的待校验参数,直到max(|fi(θ)|)<ε成立;
其中,ε为预设精度,可取1e-5。
3.一种终端,其特征在于,包括:
处理器;
用于存储处理器的执行指令的存储器;
其中,所述处理器被配置为执行权利要求1所述的方法。
4.一种存储有计算机程序的计算机可读存储介质,其特征在于,该程序被处理器执行时实现如权利要求1所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010326984.8A CN111553041B (zh) | 2020-04-23 | 2020-04-23 | 一种燃气管网参数校验方法、系统、终端及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010326984.8A CN111553041B (zh) | 2020-04-23 | 2020-04-23 | 一种燃气管网参数校验方法、系统、终端及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111553041A CN111553041A (zh) | 2020-08-18 |
CN111553041B true CN111553041B (zh) | 2021-02-26 |
Family
ID=72003955
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010326984.8A Active CN111553041B (zh) | 2020-04-23 | 2020-04-23 | 一种燃气管网参数校验方法、系统、终端及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111553041B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112307909B (zh) * | 2020-10-15 | 2024-02-20 | 上海东方投资监理有限公司 | 一种钢结构计算方法、系统及其存储介质 |
CN117235950B (zh) * | 2023-11-15 | 2024-01-30 | 上海叁零肆零科技有限公司 | 基于牛顿迭代法的天然气管网稳态仿真方法、介质及设备 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107590336A (zh) * | 2017-09-13 | 2018-01-16 | 哈尔滨理工大学 | 燃气管道泄漏对内部流场影响的数值模拟方法 |
CN110532517A (zh) * | 2019-08-28 | 2019-12-03 | 北京理工大学 | 基于改进的arukf的燃气管道参数估计方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010115942A (ja) * | 2008-11-11 | 2010-05-27 | Yoshiaki Takahashi | 摩擦抵抗低減船およびその運転方法 |
CN104359520A (zh) * | 2014-11-28 | 2015-02-18 | 成都千嘉科技有限公司 | 一种用于流量计的参数存储系统 |
CN106777731A (zh) * | 2016-12-26 | 2017-05-31 | 广东技术师范学院 | 一种燃气管网水力计算方法 |
CN108062599B (zh) * | 2017-12-12 | 2020-08-18 | 西安交通大学 | 一种电力系统和天然气系统协同运行的建模方法 |
KR102005186B1 (ko) * | 2018-04-23 | 2019-07-29 | 목포해양대학교 산학협력단 | 배기가스를 이용한 마찰저항저감장치 및 이를 구비한 단차형 선박 |
CN109919452B (zh) * | 2019-02-15 | 2023-11-28 | 三峡大学 | 一种基于多主体博弈的电力-天然气综合能源系统联合规划方法 |
CN110083901B (zh) * | 2019-04-17 | 2020-09-11 | 北京讯腾智慧科技股份有限公司 | 燃气管线校验方法和装置 |
CN110881042A (zh) * | 2019-11-29 | 2020-03-13 | 株洲新奥燃气有限公司 | 燃气管网监测数据的采集终端及其输出电压调整方法 |
-
2020
- 2020-04-23 CN CN202010326984.8A patent/CN111553041B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107590336A (zh) * | 2017-09-13 | 2018-01-16 | 哈尔滨理工大学 | 燃气管道泄漏对内部流场影响的数值模拟方法 |
CN110532517A (zh) * | 2019-08-28 | 2019-12-03 | 北京理工大学 | 基于改进的arukf的燃气管道参数估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111553041A (zh) | 2020-08-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110108328B (zh) | 一种供水管网漏损区域漏水量的获取方法 | |
CN108984873B (zh) | 供水管网实时漏损检测方法、设备、系统及存储介质 | |
US10387590B2 (en) | Techniques for iterative reduction of uncertainty in water distribution networks | |
CN111553041B (zh) | 一种燃气管网参数校验方法、系统、终端及存储介质 | |
Preis et al. | Efficient hydraulic state estimation technique using reduced models of urban water networks | |
Zhang et al. | Efficient numerical approach for simultaneous calibration of pipe roughness coefficients and nodal demands for water distribution systems | |
Díaz et al. | Observability analysis in water transport networks: Algebraic approach | |
CN102968529B (zh) | 一种供水管网模型计算结果不确定性区间的量化方法 | |
CN105912844A (zh) | 一种天然气管网供气可靠性的评价方法及装置 | |
CN111626003B (zh) | 一种供暖系统热负荷分层预测方法、存储介质及预测设备 | |
CN108197725A (zh) | 基于需水量先验信息的供水管网节点需水量校核方法 | |
CN113837451B (zh) | 油气管道的数字孪生体构建方法、装置、设备及存储介质 | |
CN111563332B (zh) | 基于夜间调压实测的城镇供水管网漏失定位方法 | |
CN112113146B (zh) | 供水管网管道粗糙系数和节点需水量同步自适应校核方法 | |
Wang et al. | Eliminating the bridge modal variability induced by thermal effects using localized modeling method | |
CN114896891A (zh) | 一种基于核极限学习机误差修正的蒸汽仿真计算方法 | |
CN103778306B (zh) | 一种基于ei和逐步消减法的传感器优化布设方法 | |
Takahashi et al. | Water distribution system model calibration under uncertainty environments | |
CN113849943A (zh) | 一种耦合压力先验信息的供水管网节点需水量校核方法 | |
CN106575312B (zh) | 用于检测配送网络特别是水配送网络中的异常的方法 | |
CN111046458B (zh) | 一种基于马尔科夫链的金属屋面性能退化模型及其方法 | |
CN115982918B (zh) | 阻力特性辨识方法、装置、设备及存储介质 | |
CN110083933A (zh) | 一种考虑随机效应的腐蚀管道贝叶斯退化分析方法 | |
CN111241647B (zh) | 核电厂电缆长度预估方法及装置 | |
Zimoch et al. | Process of hydraulic models calibration |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |