CN114936522A - 一种天然气***动态建模方法 - Google Patents
一种天然气***动态建模方法 Download PDFInfo
- Publication number
- CN114936522A CN114936522A CN202210534573.7A CN202210534573A CN114936522A CN 114936522 A CN114936522 A CN 114936522A CN 202210534573 A CN202210534573 A CN 202210534573A CN 114936522 A CN114936522 A CN 114936522A
- Authority
- CN
- China
- Prior art keywords
- model
- network
- dynamic
- air
- node
- 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.)
- Granted
Links
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 title claims abstract description 70
- 238000000034 method Methods 0.000 title claims abstract description 35
- 239000003345 natural gas Substances 0.000 title claims abstract description 35
- 239000011159 matrix material Substances 0.000 claims abstract description 54
- 239000007789 gas Substances 0.000 claims abstract description 38
- 230000008569 process Effects 0.000 claims abstract description 10
- 239000003795 chemical substances by application Substances 0.000 claims abstract description 9
- 238000013528 artificial neural network Methods 0.000 claims abstract description 6
- 238000005096 rolling process Methods 0.000 claims abstract description 6
- 238000013135 deep learning Methods 0.000 claims abstract description 5
- 238000002347 injection Methods 0.000 claims description 12
- 239000007924 injection Substances 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 2
- 238000005094 computer simulation Methods 0.000 claims 1
- 238000012549 training Methods 0.000 abstract description 16
- 238000000265 homogenisation Methods 0.000 abstract description 3
- 230000007246 mechanism Effects 0.000 abstract description 2
- 238000005260 corrosion Methods 0.000 description 3
- 230000007797 corrosion Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical group C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 239000011541 reaction mixture Substances 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000006467 substitution reaction Methods 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/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Artificial Intelligence (AREA)
- Software Systems (AREA)
- General Health & Medical Sciences (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Molecular Biology (AREA)
- Computational Linguistics (AREA)
- Mathematical Physics (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Economics (AREA)
- Computing Systems (AREA)
- Human Resources & Organizations (AREA)
- Computer Hardware Design (AREA)
- Medical Informatics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Geometry (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- Strategic Management (AREA)
- Primary Health Care (AREA)
- General Business, Economics & Management (AREA)
- Feedback Control In General (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了属于能源***技术领域的一种天然气***动态建模方法。该方法包括以下步骤:步骤1:获取气网实际的气源流量、节点压强和负荷流量的运行数据;步骤2:建立气网动态二端口时域模型,获得气网关联参数矩阵;步骤3:确定气网最大时间常数;步骤4:根据气网最大时间常数对气网实际的运行数据进行切分;步骤5:运用深度学习方法拟合天然气***动态,提取神经网络参数来构建气网动态代理模型;步骤6:将气网动态代理模型滚动使用,达到使用小规模模型刻画长时间气网动态过程。本发明相较于使用均一化参数的气网动态物理模型精度更高;引入滚动机制降低了代理模型规模与模型训练代价,减少了模型参数,提升了模型拟合精度。
Description
技术领域
本发明涉及能源***技术领域,尤其涉及一种天然气***动态建模方法。
背景技术
我国作为全球最大的能源消费国,能源发展面临环境污染严重、用能效率较低、能源安全形势严峻等诸多挑战,因此,我国能源结构转型势在必行。天然气作为最清洁的化石能源,可以通过多种耦合装置变换为其他能源形式,如:热能、电能、冷能等,而随着国家对于氢能源的关注,天然气掺氢运输成为新的研究方向,可见天然气***在能源***中具有十分重要的地位。
天然气沿管道传播动态可由一组时域内偏微分方程描述,计算代价较高,其在典型参数下千米级管道的时间常数为分钟级。若在最优能流问题中假定天然气***时刻处于稳态,虽然可避免求解时域偏微分方程,降低模型解算代价,但将导致较大的建模误差。为此,如何合理建模天然气动态成为了最优能流研究的重点与难点。
一般而言,天然气动态指天然气流量、气压沿管道的时域动态过程,表征了管道任意一端状态量(流量、气压)的“变化”向管道另一端的传导,直观表现为状态量变化时空传播的“非即时性”。为规避时域有限差分模型精度与模型解算代价难以平衡的问题,国内外学者采用信号变换思路,将天然气时域偏微分动态方程变换到频域、复频域、伯恩斯坦空间。也有工作对天然气时域偏微分动态方程做不同形式的时域有限差分格式近似,推导得到了天然气动态时域二端口模型。
以上天然气动态模型在推导过程有两点共同之处:其一,同一管道的偏微分方程系数(管道倾角、气体摩擦系数、管径)视作不变;其二,忽略了偏微分方程的某些项(如对流项和高程差)。然而,实际工程中天然气管路由大量长度较短的管道拼接而成,在施工过程中为适应地形、已有公共设施,管道局部可能存在倾角,同时,管道的腐蚀以内腐蚀为主,受内腐蚀的情况不同导致不同位置摩擦系数和管径也不相同。因此,前述模型关于管道参数均一化的假设以及模型简化条件可能造成较大的建模误差。由于天然气管道深埋于地下,每段管道的实际物理参数不易获得,增加了精细化建模天然气动态的难度。因此,需要一种全新的天然气***动态建模方法解决上述问题。
发明内容
本发明的目的是提出一种天然气***动态建模方法,其特征在于,包括以下步骤:
步骤1:获取气网实际的气源流量、节点压强和负荷流量的运行数据;
步骤2:基于显式差分法的管道动态二端口物理模型建立气网动态二端口时域模型,获得气网关联参数矩阵;
步骤3:根据步骤2所得的气网关联参数矩阵,确定气网最大时间常数;
步骤4:根据气网最大时间常数对气网实际的运行数据进行切分;
步骤5:运用深度学习方法拟合天然气***动态,提取神经网络参数来构建气网动态代理模型;
步骤6:将步骤5中的气网动态代理模型滚动使用,达到使用小规模模型刻画长时间气网动态过程。
所述步骤2具体如下:
管道p首末端压强和流量满足:
对式(1)做变形处理,有:
进一步建立节点-流出管道关联矩阵Aout、节点-流入管道关联矩阵Ain、管道首端-节点关联矩阵Apn1、管道末端-节点关联矩阵Apn2;将气网节点注入流量表示为:
Mn=AinMout-AoutMin (3)
同时,节点-管道压强满足:
Γin=Apn1Γn Γout=Apn2Γn (4)
将式(3)与式(4)代入式(2),得到:
式中:Msr为气源节点的注入流量;Mint为中间节点的注入流量,由KCL知此向量为0;Mld为负荷节点的注入流量;Γsr为气源节点的压强;Γint为中间节点的压强;Γld为负荷节点的压强,整理式(5)得:
至此,构建了气网动态二端口时域模型:
其中,矩阵Ymm,Ymp,Ypm,Ypp为Mld、Γsr与Msr、Γld的气网关联参数矩阵。
所述步骤3中的气网最大时间常数Πnet为Rnet×ΔT,其中Rnet为矩阵Ymm,Ymp,Ypm,Ypp元素坐标边界,ΔT为时间微元长度。
本发明的有益效果在于:
1、本发明提供的气网动态代理模型相较于使用均一化参数的气网动态物理模型精度更高;
2、本发明引入滚动机制降低了代理模型规模与模型训练代价,减少了模型参数,提升了模型拟合精度。
附图说明
图1为本发明天然气***动态建模方法的流程图;
图2为入口流量与出口流量的关系示意图;
图3(a)为系数矩阵A第n行列数曲线图;
图3(b)为系数矩阵A第n列行数曲线图;
图4为验证气网动态代理模型天然气***拓扑图;
图5(a)为气网动态代理模型与全时段代理模型在样本集内气源质量流量比较图;
图5(b)为气网动态代理模型与全时段代理模型在样本集内3节点压强比较图;
图5(c)为气网动态代理模型与全时段代理模型在样本集内4节点压强比较图;
图6(a)为气网动态代理模型与全时段代理模型在样本集外气源质量流量比较图;
图6(b)为气网动态代理模型与全时段代理模型在样本集外3节点压强比较图;
图6(c)为气网动态代理模型与全时段代理模型在样本集外4节点压强比较图。
具体实施方式
本发明提出一种天然气***动态建模方法,下面结合附图和具体实施例对本发明做进一步说明。
如图1所示,为本发明的实现流程,具体步骤如下:
步骤1:获取气网实际的运行数据,包括气源流量、节点压强和负荷流量等。
步骤2:基于显式差分法的管道动态二端口物理模型推导至气网二端口时域模型,获得气网关联参数矩阵。
管道p首末端压强和流量满足:
根据天然气管道参数生成A,B,C,D矩阵,上述矩阵均为严格下三角矩阵且靠近主对角线占优。其中A矩阵为Min和Mout关联系数矩阵,由于A矩阵为靠近主对角线占优的下三角矩阵,考虑压缩其维度,从而降低模型训练代价。如图2所示,等于Mout乘A矩阵的第n行与BΓin中对应元素之和,而从图3(a)不难发现,A矩阵第n行第l到n列数值远大于其他列(数值为0或接近0),即主要受到的影响,忽略Mout其他维度的影响,其中l为第n行、第1到n列中数值为ξ%最大值所对应列的序数。A矩阵的第n列表征Min对的影响,由图3(b)可知A矩阵的第n列在n到m行数值远大于其他行(数值为0或接近0),即主要受到的影响,忽略Min其他行对的影响,其中m为第n列,第n到N行中数值为ξ%最大值所对应行的序数。第n个时间微元对应的(n-l,m-n)维压缩后关联系数矩阵,表征了第n个时间微元状态量与气网其他维度状态量之间的近似关系。按照时间微元序数从低到高,依次计算各微元压缩后关联系数矩阵维度,通式为RA,n=max(n-l,m-n)。进一步,取各微元压缩后关联系数矩阵的最大维度即在满足精度条件下得到Min和Mout压缩后的关联系数矩阵维度,同理可得RB、RC、RD。按照上述对最大时间常数定义,可以推算出影响管道动态关联系数矩阵元素坐标边界R=max(RA,RB,RC,RD),单管道最大时间常数Π即为R×ΔT。
对式(9)做变形处理,有:
进一步建立节点-流出管道关联矩阵Aout、节点-流入管道关联矩阵Ain、管道首端-节点关联矩阵Apn1、管道末端-节点关联矩阵Apn2。将气网节点注入流量表示为:
Mn=AinMout-AoutMin (11)
同时,节点-管道压强满足:
Γin=Apn1Γn Γout=Apn2Γn (12)
将式(11)与式(12)代入式(10),得到:
式中:Msr为气源节点的注入流量;Mint为中间节点的注入流量,由KCL可知此向量为0;Mld为负荷节点的注入流量;Γsr为气源节点的压强;Γint为中间节点的压强,可由其他已知量表示;Γld为负荷节点的压强。整理式(13)可得:
至此,构建了气网动态时域二端口模型:
步骤3:根据步骤2所得的气网关联参数矩阵,确定气网最大时间常数。
矩阵Ymm,Ymp,Ypm,Ypp为Mld、Γsr与Msr、Γld的关联系数矩阵,关联系数矩阵Y由一个或多个Y1,Y2,…,Yn方阵组成,各方阵亦为下三角矩阵且靠近主对角线占优,与单管道动态系数矩阵类似,需要说明因此,可以根据气网最大时间常数定义,推算出影响气网动态的Ymm,Ymp,Ypm,Ypp矩阵元素坐标边界气网最大时间常数Πnet即为Rnet×ΔT。
步骤4:根据气网最大时间常数对气网实际的运行数据进行切分,以增加样本数量。
步骤5:应用深度学习方法拟合天然气***动态,提取神经网络参数构建气网动态代理模型。
步骤6:将步骤5中的代理模型滚动使用,以达到使用小规模模型刻画长时间气网动态过程。
以上步骤中,关键环节为基于显式差分法形成单管道关联参数矩阵,从气网拓扑图获取管道之间的连接关系,推导出气网关联参数矩阵,从关联参数矩阵中估算气网最大时间常数。
关键参数为气网最大时间常数,气网最大时间常数决定了如何对气网实际的运行数据进行切分、神经网络输入层与输出层的维数。
下面通过一个具体实例来对本发明方法进行说明。如图4所示,1#-3#管道长度均为2000m,节点1为气源节点,其气压维持0.4MPa不变,节点3与节点4为负荷节点。在训练样本集外选择10组工况生成10组测试样本集(与训练样本集中场景不同),在模型I(全时段代理模型)的训练样本集内选择10组数据共同组成测试样本集,比较模型I与模型II(本方法所提模型)的精度、模型参数和求解时间。
将模型I与模型II的模型精度比较结果汇总于表1。不难发现,模型I对训练样本拟合精度较高,质量流量和节点压强平均误差不超过0.4%且最大误差不超过2%,而其对训练样本集以外新样本的拟合精度较低,节点压强平均拟合误差接近5%、最大拟合误差接近10%,质量流量最大拟合误差超过1%;而模型II对训练样本与新样本的气网状态量拟合精度均较高,最大误差不超过万分之五。一方面源于模型I需要拟合的参数较多,另一方面在于模型II通过引入气网最大时间常数,实现了对训练样本的切分,间接增加了模型训练样本数量,提高了拟合精度。
表1气网动态代理模型与全时段代理模型对比
同时,从训练样本集和新样本中各选择一个工况,将模型I与模型II对气网状态量的拟合动态曲线绘制于图5(a)、图5(b)、图5(c)、图6(a)、图6(b)和图6(c)中。可以看出,模型I拟合的质量流量和压强曲线有毛刺,不平滑,对气网动态的刻画较为粗糙。相比之下,模型II在训练样本集内、外平均误差和最大误差均小于0.05%,模型精度高,泛化能力强,且其质量流量和压强曲线更平滑,贴近实际工况,准确地刻画了气网动态过程。
在模型拟合负担方面,由于模型I的输出量为样本总时长的气网未知状态量,对应代理模型的输入量和输出量维数也较大,导致参数矩阵W与B维数增加,模型I的参数数量为15,151,920,而模型II的参数数量仅为346,860,为模型I待拟合参数数量的2.29%。若样本时长增加,如由120分钟增加到天级,模型I的输出维度将持续增加,待拟合参数数量随之增长,模型训练负担大;由于模型II可对样本按气网最大时间常数切分,因此其待拟合参数数量与样本时长无关,届时优势将更加明显。由此可见,本发明方法基于天然气***历史运行数据与深度学习技术,拟合天然气动态代理模型,取消了管道参数均一化假设,提高了模型精度;既可以有效降低代理模型规模,又提高了天然气***历史运行数据的复用率,间接提升了模型训练阶段有效样本的数量。
此实施例仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (3)
1.一种天然气***动态建模方法,其特征在于,包括以下步骤:
步骤1:获取气网实际的气源流量、节点压强和负荷流量的运行数据;
步骤2:基于显式差分法的管道动态二端口物理模型建立气网动态二端口时域模型,获得气网关联参数矩阵;
步骤3:根据步骤2所得的气网关联参数矩阵,确定气网最大时间常数;
步骤4:根据气网最大时间常数对气网实际的运行数据进行切分;
步骤5:运用深度学习方法拟合天然气***动态,提取神经网络参数来构建气网动态代理模型;
步骤6:将步骤5中的气网动态代理模型滚动使用,达到使用小规模模型刻画长时间气网动态过程。
2.根据权利要求1所述的天然气***动态建模方法,其特征在于,所述步骤2具体如下:
管道p首末端压强和流量满足:
对式(1)做变形处理,有:
进一步建立节点-流出管道关联矩阵Aout、节点-流入管道关联矩阵Ain、管道首端-节点关联矩阵Apn1、管道末端-节点关联矩阵Apn2;将气网节点注入流量表示为:
Mn=AinMout-AoutMin (3)
同时,节点-管道压强满足:
Γin=Apn1Γn Γout=Apn2Γn (4)
将式(3)与式(4)代入式(2),得到:
式中:Msr为气源节点的注入流量;Mint为中间节点的注入流量,由KCL知此向量为0;Mld为负荷节点的注入流量;Γsr为气源节点的压强;Γint为中间节点的压强;Γld为负荷节点的压强,整理式(5)得:
至此,构建了气网动态二端口时域模型:
其中,矩阵Ymm,Ymp,Ypm,Ypp为Mld、Γsr与Msr、Γld的气网关联参数矩阵。
3.根据权利要求1所述的天然气***动态建模方法,其特征在于,所述步骤3中的气网最大时间常数Πnet为Rnet×ΔT,其中Rnet为矩阵Ymm,Ymp,Ypm,Ypp元素坐标边界,ΔT为时间微元长度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210534573.7A CN114936522B (zh) | 2022-05-17 | 2022-05-17 | 一种天然气***动态建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210534573.7A CN114936522B (zh) | 2022-05-17 | 2022-05-17 | 一种天然气***动态建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114936522A true CN114936522A (zh) | 2022-08-23 |
CN114936522B CN114936522B (zh) | 2023-04-07 |
Family
ID=82865595
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210534573.7A Active CN114936522B (zh) | 2022-05-17 | 2022-05-17 | 一种天然气***动态建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114936522B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150261893A1 (en) * | 2014-04-22 | 2015-09-17 | China University Of Petroleum - Beijing | Method and apparatus for determining pipeline flow status parameter of natural gas pipeline network |
CN111259547A (zh) * | 2020-01-16 | 2020-06-09 | 清华大学 | 一种用于综合能源***运行控制的天然气气路建模方法 |
CN112883658A (zh) * | 2021-01-14 | 2021-06-01 | 华南理工大学 | 一种复杂天然气潮流加速凸化调度方法 |
-
2022
- 2022-05-17 CN CN202210534573.7A patent/CN114936522B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150261893A1 (en) * | 2014-04-22 | 2015-09-17 | China University Of Petroleum - Beijing | Method and apparatus for determining pipeline flow status parameter of natural gas pipeline network |
CN111259547A (zh) * | 2020-01-16 | 2020-06-09 | 清华大学 | 一种用于综合能源***运行控制的天然气气路建模方法 |
CN112883658A (zh) * | 2021-01-14 | 2021-06-01 | 华南理工大学 | 一种复杂天然气潮流加速凸化调度方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114936522B (zh) | 2023-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Davidson et al. | The size distortion of bootstrap tests | |
CN110046408B (zh) | 基于bp神经网络微分与区间分析的区间不确定性优化方法 | |
CN109255173B (zh) | 考虑区间不确定性的结构失效概率区间计算方法 | |
Zhou et al. | Dynamic simulation of natural gas pipeline network based on interpretable machine learning model | |
CN110443407B (zh) | 一种提高供水串联管道瞬态水力模拟精度的节点流量优化分配方法 | |
CN110119853B (zh) | 基于时间序列监测数据分析的供水管网漏损报警阈值选取方法 | |
CN111475913A (zh) | 一种蒸汽动力***的运行操作优化方法及*** | |
CN114963029B (zh) | 一种基于负压波监测的多径管网泄漏检测方法 | |
CN112113146B (zh) | 供水管网管道粗糙系数和节点需水量同步自适应校核方法 | |
CN105260806A (zh) | 一种管路***的流固耦合动力学特性预测方法 | |
CN111290274B (zh) | 一种具有数据包丢失的网络控制***h∞控制方法 | |
CN103116698A (zh) | 一种基于三次样条的gm(1,1)模型预测方法 | |
CN113609721B (zh) | 多类型极端灾害电气互联***韧性计算方法及装置 | |
CN114936522B (zh) | 一种天然气***动态建模方法 | |
Zhang et al. | Efficient reliability-based design optimization for hydraulic pipeline with adaptive sampling region | |
CN112380606B (zh) | 土木结构体系中构件重要性和体系易损性评估方法 | |
CN112257355A (zh) | 一种掺氢天然气的中低压配气管网建模方法 | |
CN107808021B (zh) | 基于cfd的流体设备的阻力计算方法 | |
CN106202694A (zh) | 基于组合预测方法的组合Kriging模型构建方法 | |
CN102567606A (zh) | 给水管网等水压面实时绘制方法 | |
Ying et al. | Simulation modeling and interface parameter design of the semi‐physical braking intelligent test system | |
CN107133418B (zh) | 基于交替tvd差分算法的地球流体物质平流输运模拟方法 | |
CN110334425A (zh) | 一种用于工程结构不确定性的改进区间算法 | |
CN115688340B (zh) | 一种天然气输气管网***动态仿真的求解方法及*** | |
CN113962131B (zh) | 一种高效模拟大型天然气管网流动传热的方法 |
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 |