CN112069634A - 一种基于放宽约束的气网扩容规划系统和方法 - Google Patents

一种基于放宽约束的气网扩容规划系统和方法 Download PDF

Info

Publication number
CN112069634A
CN112069634A CN202010819675.4A CN202010819675A CN112069634A CN 112069634 A CN112069634 A CN 112069634A CN 202010819675 A CN202010819675 A CN 202010819675A CN 112069634 A CN112069634 A CN 112069634A
Authority
CN
China
Prior art keywords
network
air
module
gas
constraint
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
Application number
CN202010819675.4A
Other languages
English (en)
Other versions
CN112069634B (zh
Inventor
符政鑫
段意强
许斯滨
吕应龙
许方园
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangdong University of Technology
Original Assignee
Guangdong University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Guangdong University of Technology filed Critical Guangdong University of Technology
Priority to CN202010819675.4A priority Critical patent/CN112069634B/zh
Publication of CN112069634A publication Critical patent/CN112069634A/zh
Application granted granted Critical
Publication of CN112069634B publication Critical patent/CN112069634B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/18Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/14Pipes
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Business, Economics & Management (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • Water Supply & Treatment (AREA)
  • Tourism & Hospitality (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Public Health (AREA)
  • General Business, Economics & Management (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供一种基于放宽约束的气网扩容规划系统,包括气网数据模块、气网优化计算模块和气网成本模块;其中,气网优化计算模块包括气网约束模块、气网成本函数模块、数学建模模块和优化输出模块。本发明还提供一种基于放宽约束的气网扩容规划方法,包括以下步骤:S1:采集气网数据;S2:建立气网约束和气网成本函数;S3:构建扩容规划系统的数学模型;S4:对数学模型进行优化求解,得到优化结果;S5:计算得到规划总成本。本发明提供一种基于放宽约束的气网扩容规划系统和方法,不用使用复杂的算法,保证了求解的速度和最优解的存在,解决了现有的气网扩容规划系统建立的模型容易出现优化无解的情况的问题。

Description

一种基于放宽约束的气网扩容规划系统和方法
技术领域
本发明涉及气网规划技术领域,更具体的,涉及一种基于放宽约束的气网扩容规划系统和方法。
背景技术
气网规划是天然气网络规划的重要组成部分,其主要任务是根据规划期间的气负荷增长情况以及气源规划情况以确定相应的最佳气网网架结构。然而,现在的气网扩容规划系统在建立模型时往往采用了经典的气网暂稳态优化模型,并没有对约束条件进行修改,而是采用多目标函数并加以权重值以达到扩容规划的目的。这种扩容规划系统建立的模型比较复杂,权重值的选取非常重要,容易出现优化无解的情况。
现有技术中,如2017年5月17日公开的中国专利,一种燃气管网规划辅助决策系统,公开号为CN106682804A,通过建立现状管网模型与现状管网负荷数据之间的关联,选出合适的规划任务,但没有改进现有的气网扩容规划系统建立的模型。
发明内容
本发明为克服现有的气网扩容规划系统建立的模型容易出现优化无解的情况的技术缺陷,提供一种基于放宽约束的气网扩容规划系统和方法。
为解决上述技术问题,本发明的技术方案如下:
一种基于放宽约束的气网扩容规划系统,包括气网数据模块、气网优化计算模块和气网成本模块;其中,所述气网优化计算模块包括气网约束模块、气网成本函数模块、数学建模模块和优化输出模块;
所述气网数据模块的输出端分别与所述气网约束模块、气网成本函数模块的输入端连接,所述气网约束模块、气网成本函数模块的输出端均与所述数学建模模块的输入端连接,所述数学建模模块的输出端与所述优化输出模块的输入端连接,所述优化输出模块的输出端与所述气网成本模块连接;
所述气网数据模块用于提供气网网架结构数据、气网气源容量上下限数据、气网管道潮流上下限数据和气网负荷数据;
所述气网优化计算模块用于接收所述气网数据模块提供的数据,并结合气网约束条件和规划成本函数建立扩容规划系统模型,最后求解扩容规划系统模型获得优化结果并输出;
所述气网约束模块用于提供气网约束条件;
所述气网成本函数模块用于计算气网规划成本;
所述数学建模模块用于根据气网约束条件和气网规划成本函数构建扩容规划系统的数学模型;
所述优化输出模块用于优化求解扩容规划系统的数学模型并输出优化结果;
所述气网成本模块用于接收所述优化输出模块输出的优化结果,计算得到规划总成本。
优选的,所述气网约束模块包括气网负荷平衡约束模块、气网节点气压约束模块、气网潮流约束模块和气网气源容量约束模块;所述气网负荷平衡约束模块、气网节点气压约束模块、气网潮流约束模块、气网气源容量约束模块的输入端均与所述气网数据模块的输出端连接;所述气网负荷平衡约束模块、气网节点气压约束模块、气网潮流约束模块、气网气源容量约束模块的输出端均与所述数学建模模块的输入端连接;
所述气网负荷平衡约束模块用于提供气网负荷平衡约束条件;
所述气网节点气压约束模块用于提供气网节点气压约束条件;
所述气网潮流约束模块用于提供气网潮流约束条件;
所述气网气源容量约束模块用于提供气网气源容量约束条件。
优选的,所述气网负荷平衡约束模块提供的气网负荷平衡约束条件为:
Figure BDA0002634012050000021
所述气网节点气压约束模块提供的气网节点气压约束条件为:
Figure BDA0002634012050000022
所述气网潮流约束模块提供的气网潮流约束条件为:
Figure BDA0002634012050000023
所述气网气源容量约束模块提供的气网气源容量约束条件为:
GWwt≥GWw,min
其中,GWwt表示气网气源w在时刻t的供气量,W为气网气源数目,GDt表示气网在时刻t的气网负荷矩阵;wm,wn分别表示气网管道的首端节点和末端节点的气压,wm,min,wn,min分别表示该节点的最小气压值;GLl为气网管道l潮流,Cp为韦茅斯常数,通常大于0,wl,m,wl,n分别为气网管道l首端节点和末端节点的气压,GLl,min为气网对应管道l潮流下限值;GWw,min表示气网气源w供气量的下限值。
优选的,所述气网成本函数模块包括气网运营购气成本函数模块、气网气源容量扩容建设成本函数模块和气网管道扩容建设成本函数模块;所述气网运营购气成本函数模块、气网气源容量扩容建设成本函数模块和气网管道扩容建设成本函数模块的输入端均与所述气网数据模块的输出端连接;所述气网运营购气成本函数模块、气网气源容量扩容建设成本函数模块和气网管道扩容建设成本函数模块的输出端均与所述数学建模模块的输入端连接;
所述气网运营购气成本函数模块用于计算气网运营时为满足气网负荷的购气成本;
所述气网气源容量扩容建设成本函数模块用于计算当气网气源无法满足当前气网负荷需要对气源容量进行扩容所需成本;
所述气网管道扩容建设成本函数模块用于计算气网管道容量无法满足当前气网气源供气容量需求时对管道容量进行扩容所需成本。
优选的,所述气网运营购气成本函数模块中计算气网运营时为满足气网负荷的购气成本的函数为
Figure BDA0002634012050000031
所述气网气源容量扩容建设成本函数模块中计算气网气源容量扩容建设成本的函数为:
SWwt=u(GWwt-GWw,max)*(GWwt-GWw,max)*Φw
所述气网管道扩容建设成本函数模块中计算气网管道扩容建设成本的函数为:
SLlt=u(GLlt-GLl,max)*(GLlt-GLl,max)*ωl
其中,GWwt表示气网气源w在时刻t的供气量,W为气网气源数目,GW_pricew表示气网气源w的供气成本;SWwt为第w个气源在时刻t的容量扩容建设成本,函数u()表示替代可导函数,GWw,max为第w个气源的供气容量上限,Φw为第w个气源的扩容建设价格系数;SLlt为第l条气网管道在时刻t的扩容建设成本,GLlt为第l条气网管道在时刻t的潮流,GLl,max为线路l潮流上限,ωl为第l条气网管道的扩容建设价格系数。
优选的,所述数学建模模块构建的扩容规划系统的数学模型为:
目标函数:
Figure BDA0002634012050000041
约束条件:
Figure BDA0002634012050000042
Figure BDA0002634012050000043
Figure BDA0002634012050000044
GWwt≥GWw,min
其中,NW为气网气源数量,NL为气网管道数量,GLl为气网管道l潮流。
优选的,所述优化输出模块包括气网气源供气量数据模块和气网管道潮流数据模块,所述气网气源供气量数据模块和气网管道潮流数据模块的输入端均与所述数学建模模块的输出端连接,所述气网气源供气量数据模块和气网管道潮流数据模块的输出端均与所述气网成本模块的输入端连接;
所述气网气源供气量数据模块通过对数学建模模块中的数学模型进行优化计算得到气网气源供气量数据;
所述气网管道潮流数据模块通过对数学建模模块中的数学模型进行优化计算得到气网管道潮流数据。
优选的,所述气网成本模块包括气网购气成本矩阵模块、气源容量扩容建设成本矩阵模块、气网管道扩容建设成本矩阵模块和规划总成本模块;所述气网购气成本矩阵模块、气源容量扩容建设成本矩阵模块的输入端均与所述气网气源供气量数据模块的输出端连接,所述气网管道扩容建设成本矩阵模块的输入端与所述气网管道潮流数据模块的输出端连接;所述气网购气成本矩阵模块、气源容量扩容建设成本矩阵模块、气网管道扩容建设成本矩阵模块的输出端均与所述规划总成本模块的输入端连接;
所述气网购气成本矩阵模块用于提供气网购气成本矩阵;
所述气源容量扩容建设成本矩阵模块用于提供气源容量扩容建设成本矩阵;
所述气网管道扩容建设成本矩阵模块用于提供气网管道扩容建设成本矩阵;
所述规划总成本模块用于计算规划总成本。
优选的,气网购气成本矩阵Cost1为:
Cost1=GWNW×1⊙GW_priceNW×1
气源容量扩容建设成本矩阵Cost2为:
Figure BDA0002634012050000051
气网管道扩容建设成本矩阵Cost3为:
Figure BDA0002634012050000052
则规划总成本T_Cost的计算公式为:
T_Cost=∑Cost1+∑Cost2+∑Cost3
其中,⊙表示矩阵的点乘;GWNW×1为气网优化计算模块输出的NW×1气源供气量矩阵;
Figure BDA0002634012050000053
表示NW×1的气源供气量上限矩阵;GLNL×1为气网优化计算模块输出的NL×1气网潮流数据矩阵;
Figure BDA0002634012050000054
为NL×1的气网管道潮流上限矩阵;GW_priceNW×1为对应的NW×1气源供气价格系数矩阵;ΦNW×1为NW×1的气网气源容量扩容建设价格系数矩阵;ωNL×1为NL×1的气网管道扩容建设价格系数矩阵。
一种基于放宽约束的气网扩容规划方法,通过所述的一种基于放宽约束的气网扩容规划系统实现,包括以下步骤:
S1:通过气网数据模块采集气网数据,并将气网数据输入到气网约束模块和气网成本函数模块中;
S2:分别在气网约束模块和气网成本函数模块中建立气网约束和气网成本函数;其中,气网约束包括气网负荷平衡约束、气网节点气压约束、气网潮流约束和气网气源容量约束;气网成本函数包括气网运营购气成本函数、气网气源容量扩容建设成本函数和气网管道扩容建设成本函数;
S3:将气网约束和气网成本函数输入数学建模模块中构建扩容规划系统的数学模型;
S4:通过优化输出模块对数学模型进行优化求解,得到优化结果;
S5:将优化结果输入气网成本模块中,计算得到规划总成本。
与现有技术相比,本发明技术方案的有益效果是:
本发明提供了一种基于放宽约束的气网扩容规划系统和方法,提出了新的约束条件,并将新的约束条件结合规划成本构建数学模型,该数学模型线性化程度高,可以通过一些经典的优化算法进行快速求解,不用使用复杂的算法,保证了求解的速度和最优解的存在。
附图说明
图1为本发明中一实施例的模块连接示意图;
图2为本发明中一实施例的实施步骤流程示意图;
其中:1、气网数据模块;2、气网优化计算模块;21、气网约束模块;211、气网负荷平衡约束模块;212、气网节点气压约束模块;213、气网潮流约束模块;214、气网气源容量约束模块;22、气网成本函数模块;221、气网运营购气成本函数模块;222、气网气源容量扩容建设成本函数模块;223、气网管道扩容建设成本函数模块;23、数学建模模块;24、优化输出模块;241、气网气源供气量数据模块;242、气网管道潮流数据模块;3、气网成本模块;311、气网购气成本矩阵模块;312、气源容量扩容建设成本矩阵模块;313、气网管道扩容建设成本矩阵模块;314、规划总成本模块。
具体实施方式
附图仅用于示例性说明,不能理解为对本专利的限制;
为了更好说明本实施例,附图某些部件会有省略、放大或缩小,并不代表实际产品的尺寸;
对于本领域技术人员来说,附图中某些公知结构及其说明可能省略是可以理解的。
下面结合附图和实施例对本发明的技术方案做进一步的说明。
实施例1
如图1所示,一种基于放宽约束的气网扩容规划系统,包括气网数据模块1、气网优化计算模块2和气网成本模块3;其中,所述气网优化计算模块2包括气网约束模块21、气网成本函数模块22、数学建模模块23和优化输出模块24;
所述气网数据模块1的输出端分别与所述气网约束模块21、气网成本函数模块22的输入端连接,所述气网约束模块21、气网成本函数模块22的输出端均与所述数学建模模块23的输入端连接,所述数学建模模块23的输出端与所述优化输出模块24的输入端连接,所述优化输出模块24的输出端与所述气网成本模块3连接;
所述气网数据模块1用于提供气网网架结构数据、气网气源容量上下限数据、气网管道潮流上下限数据和气网负荷数据;
所述气网优化计算模块2用于接收所述气网数据模块1提供的数据,并结合气网约束条件和规划成本函数建立扩容规划系统模型,最后求解扩容规划系统模型获得优化结果并输出;
所述气网约束模块21用于提供气网约束条件;
所述气网成本函数模块22用于计算气网规划成本;
所述数学建模模块23用于根据气网约束条件和气网规划成本函数构建扩容规划系统的数学模型;
所述优化输出模块24用于优化求解扩容规划系统的数学模型并输出优化结果;
所述气网成本模块3用于接收所述优化输出模块24输出的优化结果,计算得到规划总成本。
在具体实施过程中,气网数据模块1收集气网网架数据、气网负荷数据、气源容量约束数据,气网管道潮流数据、气网运营购气价格数据、气源扩容价格数据、气网管道扩容价格数据等基础数据。
更具体的,所述气网约束模块21包括气网负荷平衡约束模块211、气网节点气压约束模块212、气网潮流约束模块213和气网气源容量约束模块214;所述气网负荷平衡约束模块211、气网节点气压约束模块212、气网潮流约束模块213、气网气源容量约束模块214的输入端均与所述气网数据模块1的输出端连接;所述气网负荷平衡约束模块211、气网节点气压约束模块212、气网潮流约束模块213、气网气源容量约束模块214的输出端均与所述数学建模模块23的输入端连接;
所述气网负荷平衡约束模块211用于提供气网负荷平衡约束条件;
所述气网节点气压约束模块212用于提供气网节点气压约束条件;
所述气网潮流约束模块213用于提供气网潮流约束条件;
所述气网气源容量约束模块214用于提供气网气源容量约束条件。
更具体的,所述气网负荷平衡约束模块211提供的气网负荷平衡约束条件为:
Figure BDA0002634012050000081
所述气网节点气压约束模块212提供的气网节点气压约束条件为:
Figure BDA0002634012050000082
所述气网潮流约束模块213提供的气网潮流约束条件为:
Figure BDA0002634012050000083
所述气网气源容量约束模块214提供的气网气源容量约束条件为:
GWwt≥GWw,min
其中,GWwt表示气网气源w在时刻t的供气量,W为气网气源数目,GDt表示气网在时刻t的气网负荷矩阵;wm,wn分别表示气网管道的首端节点和末端节点的气压,wm,min,wn,min分别表示该节点的最小气压值;GLl为气网管道l潮流,Cp为韦茅斯常数,通常大于0,wl,m,wl,n分别为气网管道l首端节点和末端节点的气压,GLl,min为气网对应管道l潮流下限值;GWw,min表示气网气源w供气量的下限值。
在具体实施过程中,气网负荷平衡约束为典型的气网负荷平衡约束条件。气网节点气压约束与经典气网优化模型的约束条件相比,其去掉了节点气压上限值,节点气压与流经该节点的气网潮流有关系,当气网供气量增加时,通过该约束条件可以保证气网扩容规划系统模型不会破约束。气网潮流约束模块213与典型的气网优化模型的约束条件相比,其去掉了上限值,气网管道潮流与气网节点气压相关,两者为正相关的关系,当节点气压增加时,气网管道潮流也会随之增加。由于在未来的规划中,气源的供气量有可能是一个非常大的数值,相应的气网节点的气压值也会非常大,所以气网管道的潮流上限必须是一个弹性值。气网气源容量约束与经典气网优化模型的约束条件相比,去掉了上限值,由于气网负荷具有不确定性,在规划中,负荷往往都是处于增加的情况。因此采用该气网气源容量约束就可以在气网负荷增加的情况下对气网气源进行扩容以满足未来气网运行需求。
更具体的,所述气网成本函数模块22包括气网运营购气成本函数模块221、气网气源容量扩容建设成本函数模块222和气网管道扩容建设成本函数模块223;所述气网运营购气成本函数模块221、气网气源容量扩容建设成本函数模块222和气网管道扩容建设成本函数模块223的输入端均与所述气网数据模块1的输出端连接;所述气网运营购气成本函数模块221、气网气源容量扩容建设成本函数模块222和气网管道扩容建设成本函数模块223的输出端均与所述数学建模模块23的输入端连接;
所述气网运营购气成本函数模块221用于计算气网运营时为满足气网负荷的购气成本;
所述气网气源容量扩容建设成本函数模块222用于计算当气网气源无法满足当前气网负荷需要对气源容量进行扩容所需成本;
所述气网管道扩容建设成本函数模块223用于计算气网管道容量无法满足当前气网气源供气容量需求时对管道容量进行扩容所需成本。
更具体的,所述气网运营购气成本函数模块221中计算气网运营时为满足气网负荷的购气成本的函数为
Figure BDA0002634012050000091
所述气网气源容量扩容建设成本函数模块222中计算气网气源容量扩容建设成本的函数为:
SWwt=u(GWwt-GWw,max)*(GWwt-GWw,max)*Φw
所述气网管道扩容建设成本函数模块223中计算气网管道扩容建设成本的函数为:
SLlt=u(GLlt-GLl,max)*(GLlt-GLl,max)*ωl
其中,GWwt表示气网气源w在时刻t的供气量,W为气网气源数目,GW_pricew表示气网气源w的供气成本;SWwt为第w个气源在时刻t的容量扩容建设成本,函数u()表示替代可导函数,GWw,max为第w个气源的供气容量上限,Φw为第w个气源的扩容建设价格系数;SLlt为第l条气网管道在时刻t的扩容建设成本,GLlt为第l条气网管道在时刻t的潮流,GLl,max为线路l潮流上限,ωl为第l条气网管道的扩容建设价格系数。
在具体实施过程中,气网气源容量扩容建设成本函数模块222接收气网数据,并用阶跃函数选择是否需要对气源进行扩容,以气网气源w在时刻t的供气容量与其最大供气容量的差值作为判断数据,当判断数据为正值时,判断数据大小不变,反之将判断数据修正为零;
气网管道扩容建设成本函数模块223接收管道潮流数据,并用阶跃函数选择是否需要对气网管道进行扩容,以气网管道l在时刻t的潮流与管道潮流最大值的差值作为判断数据,当判断数据为正值时,判断数据大小不变,反之将判断数据修正为零;
由于阶跃函数不可导,因此由以下替代可导函数函数u()近似获得判断数据的修正值,计算公式为:
Figure BDA0002634012050000101
其中,x为判断数据,α为一个大于零的数据,当α的值越大,该替代可导函数越接近阶跃函数。
更具体的,所述数学建模模块23构建的扩容规划系统的数学模型为:
目标函数:
Figure BDA0002634012050000102
约束条件:
Figure BDA0002634012050000103
Figure BDA0002634012050000104
Figure BDA0002634012050000111
GWwt≥GWw,min
其中,NW为气网气源数量,NL为气网管道数量,GLl为气网管道l潮流。
在具体实施过程中,以总成本最低作为目标函数构建数学模型。通过如分支定界法、拉格朗日松弛法、遗传算法等优化算法求解该数学模型,获得优化结果。
更具体的,所述优化输出模块24包括气网气源供气量数据模块241和气网管道潮流数据模块242,所述气网气源供气量数据模块241和气网管道潮流数据模块242的输入端均与所述数学建模模块23的输出端连接,所述气网气源供气量数据模块241和气网管道潮流数据模块242的输出端均与所述气网成本模块3的输入端连接;
所述气网气源供气量数据模块241通过对数学建模模块23中的数学模型进行优化计算得到气网气源供气量数据;
所述气网管道潮流数据模块242通过对数学建模模块23中的数学模型进行优化计算得到气网管道潮流数据。
在具体实施过程中,输出的优化结果为气网气源供气量数据和气网管道潮流数据。
更具体的,所述气网成本模块3包括气网购气成本矩阵模块311、气源容量扩容建设成本矩阵模块312、气网管道扩容建设成本矩阵模块313和规划总成本模块314;所述气网购气成本矩阵模块311、气源容量扩容建设成本矩阵模块312的输入端均与所述气网气源供气量数据模块241的输出端连接,所述气网管道扩容建设成本矩阵模块313的输入端与所述气网管道潮流数据模块242的输出端连接;所述气网购气成本矩阵模块311、气源容量扩容建设成本矩阵模块312、气网管道扩容建设成本矩阵模块313的输出端均与所述规划总成本模块314的输入端连接;
所述气网购气成本矩阵模块311用于提供气网购气成本矩阵;
所述气源容量扩容建设成本矩阵模块312用于提供气源容量扩容建设成本矩阵;
所述气网管道扩容建设成本矩阵模块313用于提供气网管道扩容建设成本矩阵;
所述规划总成本模块314用于计算规划总成本。
更具体的,气网购气成本矩阵Cost1为:
Cost1=GWNW×1⊙GW_priceNW×1
气源容量扩容建设成本矩阵Cost2为:
Figure BDA0002634012050000121
气网管道扩容建设成本矩阵Cost3为:
Figure BDA0002634012050000122
则规划总成本T_Cost的计算公式为:
T_Cost=∑Cost1+∑Cost2+∑Cost3
其中,⊙表示矩阵的点乘;GWNW×1为气网优化计算模块输出的NW×1气源供气量矩阵;
Figure BDA0002634012050000123
表示NW×1的气源供气量上限矩阵;GLNL×1为气网优化计算模块输出的NL×1气网潮流数据矩阵;
Figure BDA0002634012050000124
为NL×1的气网管道潮流上限矩阵;GW_priceNW×1为对应的NW×1气源供气价格系数矩阵;ΦNW×1为NW×1的气网气源容量扩容建设价格系数矩阵;ωNL×1为NL×1的气网管道扩容建设价格系数矩阵。
在具体实施过程中,利用优化结果GWNW×1与已知的
Figure BDA0002634012050000125
相减,得到一个NW×1的差数矩阵,利用阶跃函数
Figure BDA0002634012050000126
进行判断,当差数小于0为0,反之为1;然后,点乘该差数矩阵与气网气源容量扩容建设价格系数矩阵便可得到气网气源容量扩容建设成本矩阵Cost2。同理,可得气网管道扩容建设成本矩阵Cost3。该规划系统考虑了气网运营购气成本、气网气源扩容建设成本以及气网管道扩容建设成本,将各个成本加入到目标函数,构建了以总成本最低为目标函数的数学模型,该数学模型线性化程度高,可以通过一些经典的优化算法进行快速求解,不用使用复杂的算法,保证了求解的速度和最优解的存在,实现了气网扩容规划经济性最优。
实施例2
如图2所示,一种基于放宽约束的气网扩容规划方法,通过所述的一种基于放宽约束的气网扩容规划系统实现,包括以下步骤:
S1:通过气网数据模块1采集气网数据,并将气网数据输入到气网约束模块21和气网成本函数模块22中;
S2:分别在气网约束模块21和气网成本函数模块22中建立气网约束和气网成本函数;其中,气网约束包括气网负荷平衡约束、气网节点气压约束、气网潮流约束和气网气源容量约束;气网成本函数包括气网运营购气成本函数、气网气源容量扩容建设成本函数和气网管道扩容建设成本函数;
S3:将气网约束和气网成本函数输入数学建模模块23中构建扩容规划系统的数学模型;
S4:通过优化输出模块24对数学模型进行优化求解,得到优化结果;
S5:将优化结果输入气网成本模块3中,计算得到规划总成本。
在具体实施过程中,通过提出新的约束条件,并将新的约束条件结合规划成本构建数学模型,构建的数学模型以总成本最低为目标函数,且其线性化程度高,可以通过一些经典的优化算法进行快速求解,不用使用复杂的算法,保证了求解的速度和最优解的存在。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明权利要求的保护范围之内。

Claims (10)

1.一种基于放宽约束的气网扩容规划系统,其特征在于,包括气网数据模块、气网优化计算模块和气网成本模块;其中,所述气网优化计算模块包括气网约束模块、气网成本函数模块、数学建模模块和优化输出模块;
所述气网数据模块的输出端分别与所述气网约束模块、气网成本函数模块的输入端连接,所述气网约束模块、气网成本函数模块的输出端均与所述数学建模模块的输入端连接,所述数学建模模块的输出端与所述优化输出模块的输入端连接,所述优化输出模块的输出端与所述气网成本模块连接;
所述气网数据模块用于提供气网网架结构数据、气网气源容量上下限数据、气网管道潮流上下限数据和气网负荷数据;
所述气网优化计算模块用于接收所述气网数据模块提供的数据,并结合气网约束条件和规划成本函数建立扩容规划系统模型,最后求解扩容规划系统模型获得优化结果并输出;
所述气网约束模块用于提供气网约束条件;
所述气网成本函数模块用于计算气网规划成本;
所述数学建模模块用于根据气网约束条件和气网规划成本函数构建扩容规划系统的数学模型;
所述优化输出模块用于优化求解扩容规划系统的数学模型并输出优化结果;
所述气网成本模块用于接收所述优化输出模块输出的优化结果,计算得到规划总成本。
2.根据权利要求1所述的一种基于放宽约束的气网扩容规划系统,其特征在于,所述气网约束模块包括气网负荷平衡约束模块、气网节点气压约束模块、气网潮流约束模块和气网气源容量约束模块;所述气网负荷平衡约束模块、气网节点气压约束模块、气网潮流约束模块、气网气源容量约束模块的输入端均与所述气网数据模块的输出端连接;所述气网负荷平衡约束模块、气网节点气压约束模块、气网潮流约束模块、气网气源容量约束模块的输出端均与所述数学建模模块的输入端连接;
所述气网负荷平衡约束模块用于提供气网负荷平衡约束条件;
所述气网节点气压约束模块用于提供气网节点气压约束条件;
所述气网潮流约束模块用于提供气网潮流约束条件;
所述气网气源容量约束模块用于提供气网气源容量约束条件。
3.根据权利要求2所述的一种基于放宽约束的气网扩容规划系统,其特征在于,所述气网负荷平衡约束模块提供的气网负荷平衡约束条件为:
Figure FDA0002634012040000021
所述气网节点气压约束模块提供的气网节点气压约束条件为:
Figure FDA0002634012040000022
所述气网潮流约束模块提供的气网潮流约束条件为:
Figure FDA0002634012040000023
所述气网气源容量约束模块提供的气网气源容量约束条件为:
GWwt≥GWw,min
其中,GWwt表示气网气源w在时刻t的供气量,W为气网气源数目,GDt表示气网在时刻t的气网负荷矩阵;wm,wn分别表示气网管道的首端节点和末端节点的气压,wm,min,wn,min分别表示该节点的最小气压值;GLl为气网管道l潮流,Cp为韦茅斯常数,wl,m,wl,n分别为气网管道l首端节点和末端节点的气压,GLl,min为气网对应管道l潮流下限值;GWw,min表示气网气源w供气量的下限值。
4.根据权利要求3所述的一种基于放宽约束的气网扩容规划系统,其特征在于,所述气网成本函数模块包括气网运营购气成本函数模块、气网气源容量扩容建设成本函数模块和气网管道扩容建设成本函数模块;所述气网运营购气成本函数模块、气网气源容量扩容建设成本函数模块和气网管道扩容建设成本函数模块的输入端均与所述气网数据模块的输出端连接;所述气网运营购气成本函数模块、气网气源容量扩容建设成本函数模块和气网管道扩容建设成本函数模块的输出端均与所述数学建模模块的输入端连接;
所述气网运营购气成本函数模块用于计算气网运营时为满足气网负荷的购气成本;
所述气网气源容量扩容建设成本函数模块用于计算当气网气源无法满足当前气网负荷需要对气源容量进行扩容所需成本;
所述气网管道扩容建设成本函数模块用于计算气网管道容量无法满足当前气网气源供气容量需求时对管道容量进行扩容所需成本。
5.根据权利要求4所述的一种基于放宽约束的气网扩容规划系统,其特征在于,所述气网运营购气成本函数模块中计算气网运营时为满足气网负荷的购气成本的函数为
Figure FDA0002634012040000031
所述气网气源容量扩容建设成本函数模块中计算气网气源容量扩容建设成本的函数为:
SWwt=u(GWwt-GWw,max)*(GWwt-GWw,max)*Φw
所述气网管道扩容建设成本函数模块中计算气网管道扩容建设成本的函数为:
SLlt=u(GLlt-GLl,max)*(GLlt-GLl,max)*ωl
其中,GWwt表示气网气源w在时刻t的供气量,W为气网气源数目,GW_pricew表示气网气源w的供气成本;SWwt为第w个气源在时刻t的容量扩容建设成本,函数u()表示替代可导函数,GWw,max为第w个气源的供气容量上限,Φw为第w个气源的扩容建设价格系数;SLlt为第l条气网管道在时刻t的扩容建设成本,GLlt为第l条气网管道在时刻t的潮流,GLl,max为线路l潮流上限,ωl为第l条气网管道的扩容建设价格系数。
6.根据权利要求5所述的一种基于放宽约束的气网扩容规划系统,其特征在于,所述数学建模模块构建的扩容规划系统的数学模型为:
目标函数:
Figure FDA0002634012040000032
约束条件:
Figure FDA0002634012040000033
Figure FDA0002634012040000034
Figure FDA0002634012040000041
GWwt≥GWw,min
其中,NW为气网气源数量,NL为气网管道数量,GLl为气网管道l潮流。
7.根据权利要求6所述的一种基于放宽约束的气网扩容规划系统,其特征在于,所述优化输出模块包括气网气源供气量数据模块和气网管道潮流数据模块,所述气网气源供气量数据模块和气网管道潮流数据模块的输入端均与所述数学建模模块的输出端连接,所述气网气源供气量数据模块和气网管道潮流数据模块的输出端均与所述气网成本模块的输入端连接;
所述气网气源供气量数据模块通过对数学建模模块中的数学模型进行优化计算得到气网气源供气量数据;
所述气网管道潮流数据模块通过对数学建模模块中的数学模型进行优化计算得到气网管道潮流数据。
8.根据权利要求7所述的一种基于放宽约束的气网扩容规划系统,其特征在于,所述气网成本模块包括气网购气成本矩阵模块、气源容量扩容建设成本矩阵模块、气网管道扩容建设成本矩阵模块和规划总成本模块;所述气网购气成本矩阵模块、气源容量扩容建设成本矩阵模块的输入端均与所述气网气源供气量数据模块的输出端连接,所述气网管道扩容建设成本矩阵模块的输入端与所述气网管道潮流数据模块的输出端连接;所述气网购气成本矩阵模块、气源容量扩容建设成本矩阵模块、气网管道扩容建设成本矩阵模块的输出端均与所述规划总成本模块的输入端连接;
所述气网购气成本矩阵模块用于提供气网购气成本矩阵;
所述气源容量扩容建设成本矩阵模块用于提供气源容量扩容建设成本矩阵;
所述气网管道扩容建设成本矩阵模块用于提供气网管道扩容建设成本矩阵;
所述规划总成本模块用于计算规划总成本。
9.根据权利要求8所述的一种基于放宽约束的气网扩容规划系统,其特征在于,气网购气成本矩阵Cost1为:
Cost1=GWNW×1⊙GW_priceNW×1
气源容量扩容建设成本矩阵Cost2为:
Figure FDA0002634012040000042
气网管道扩容建设成本矩阵Cost3为:
Figure FDA0002634012040000051
则规划总成本T_Cost的计算公式为:
T_Cost=∑Cost1+∑Cost2+∑Cost3
其中,⊙表示矩阵的点乘;GWNW×1为气网优化计算模块输出的NW×1气源供气量矩阵;
Figure FDA0002634012040000052
表示NW×1的气源供气量上限矩阵;GLNL×1为气网优化计算模块输出的NL×1气网潮流数据矩阵;
Figure FDA0002634012040000053
为NL×1的气网管道潮流上限矩阵;GW_priceNW×1为对应的NW×1气源供气价格系数矩阵;ΦNW×1为NW×1的气网气源容量扩容建设价格系数矩阵;ωNL×1为NL×1的气网管道扩容建设价格系数矩阵。
10.一种基于放宽约束的气网扩容规划方法,其特征在于,包括以下步骤:
S1:通过气网数据模块采集气网数据,并将气网数据输入到气网约束模块和气网成本函数模块中;
S2:分别在气网约束模块和气网成本函数模块中建立气网约束和气网成本函数;其中,气网约束包括气网负荷平衡约束、气网节点气压约束、气网潮流约束和气网气源容量约束;气网成本函数包括气网运营购气成本函数、气网气源容量扩容建设成本函数和气网管道扩容建设成本函数;
S3:将气网约束和气网成本函数输入数学建模模块中构建扩容规划系统的数学模型;
S4:通过优化输出模块对数学模型进行优化求解,得到优化结果;
S5:将优化结果输入气网成本模块中,计算得到规划总成本。
CN202010819675.4A 2020-08-14 2020-08-14 一种基于放宽约束的气网扩容规划系统和方法 Expired - Fee Related CN112069634B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010819675.4A CN112069634B (zh) 2020-08-14 2020-08-14 一种基于放宽约束的气网扩容规划系统和方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010819675.4A CN112069634B (zh) 2020-08-14 2020-08-14 一种基于放宽约束的气网扩容规划系统和方法

Publications (2)

Publication Number Publication Date
CN112069634A true CN112069634A (zh) 2020-12-11
CN112069634B CN112069634B (zh) 2022-08-09

Family

ID=73661765

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010819675.4A Expired - Fee Related CN112069634B (zh) 2020-08-14 2020-08-14 一种基于放宽约束的气网扩容规划系统和方法

Country Status (1)

Country Link
CN (1) CN112069634B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050057211A1 (en) * 2003-09-11 2005-03-17 Stefan Roepke Power control device of an electric tool, a method of power control of an electric tool and an electric tool provided with a power control device
CN106529740A (zh) * 2016-12-08 2017-03-22 西安交通大学 一种天然气网络、电力网络和电源的联合规划方法
US9890908B1 (en) * 2017-04-18 2018-02-13 Air Products And Chemicals, Inc. Control system in a gas pipeline network to increase capacity factor
CN109598377A (zh) * 2018-11-28 2019-04-09 国网江苏省电力有限公司 一种基于故障约束的交直流混合配电网鲁棒规划方法
CN109975119A (zh) * 2019-04-15 2019-07-05 四川大学 一种岩石双轴压缩爆破设计方法
CN110070213A (zh) * 2019-03-28 2019-07-30 广东工业大学 一种电-气综合能源系统的日前调度方法
CN110263387A (zh) * 2019-05-29 2019-09-20 国网福建省电力有限公司 基于能源集线器的配电网、天然气网建立的综合能源系统优化运行模型及其线性处理方法
CN111371091A (zh) * 2020-04-10 2020-07-03 南京润北智能环境研究院有限公司 基于市场效益驱动的光伏发电智能预测系统的计算方法
CN111401628A (zh) * 2020-03-13 2020-07-10 南京润北智能环境研究院有限公司 一种基于市场效益驱动的光伏发电智能预测系统
US20200257971A1 (en) * 2019-01-11 2020-08-13 Chongqing University Full-linear model for optimal power flow of integrated power and natural-gas system based on deep learning methods

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050057211A1 (en) * 2003-09-11 2005-03-17 Stefan Roepke Power control device of an electric tool, a method of power control of an electric tool and an electric tool provided with a power control device
CN106529740A (zh) * 2016-12-08 2017-03-22 西安交通大学 一种天然气网络、电力网络和电源的联合规划方法
US9890908B1 (en) * 2017-04-18 2018-02-13 Air Products And Chemicals, Inc. Control system in a gas pipeline network to increase capacity factor
CN109598377A (zh) * 2018-11-28 2019-04-09 国网江苏省电力有限公司 一种基于故障约束的交直流混合配电网鲁棒规划方法
US20200257971A1 (en) * 2019-01-11 2020-08-13 Chongqing University Full-linear model for optimal power flow of integrated power and natural-gas system based on deep learning methods
CN110070213A (zh) * 2019-03-28 2019-07-30 广东工业大学 一种电-气综合能源系统的日前调度方法
CN109975119A (zh) * 2019-04-15 2019-07-05 四川大学 一种岩石双轴压缩爆破设计方法
CN110263387A (zh) * 2019-05-29 2019-09-20 国网福建省电力有限公司 基于能源集线器的配电网、天然气网建立的综合能源系统优化运行模型及其线性处理方法
CN111401628A (zh) * 2020-03-13 2020-07-10 南京润北智能环境研究院有限公司 一种基于市场效益驱动的光伏发电智能预测系统
CN111371091A (zh) * 2020-04-10 2020-07-03 南京润北智能环境研究院有限公司 基于市场效益驱动的光伏发电智能预测系统的计算方法

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
A TALEBI ET AL: "Co-optimization of Electricity and Natural Gas Networks Considering AC Constraints and Natural Gas Storage", 《2018 SMART GRID CONFERENCE (SGC)》 *
JIANG YUEWEN WT AL: "Research on Wind Power Accommodation for an Electricity-Heat-Gas Integrated Microgrid System With Power-to-Gas", 《IEEE ACCESS》 *
M.SEDGHI ET AL: "Distribution network expansion considering distributed generation and storage units using modified PSO algorithm", 《INTERNATIONAL JOURNAL OF ELECTRICAL POWER & ENERGY SYSTEMS》 *
ODETAYO B ET AL: "A chance constrained programming approach to the integrated planning of electric power generation, natural gas network and storage", 《IEEE TRANSACTIONS ON POWER SYSTEMS》 *
周贤正: "城市综合能源系统配电-气-热网联合规划研究", 《中国博士学位论文全文数据库》 *
杨经纬 等: "多能源网络的广义电路分析理论——(一)支路模型", 《电力系统自动化》 *
洪彬倬 等: "考虑交通信息与配电网结构的充电站规划", 《广东电力》 *
葛晓艳: "城市轨道项目运营成本预测模型研究", 《中国优秀硕士学位论文全文数据库 经济与管理科学辑》 *

Also Published As

Publication number Publication date
CN112069634B (zh) 2022-08-09

Similar Documents

Publication Publication Date Title
Rastgou et al. Improved harmony search algorithm for electrical distribution network expansion planning in the presence of distributed generators
Simab et al. Multi-objective programming of pumped-hydro-thermal scheduling problem using normal boundary intersection and VIKOR
CN109636052B (zh) 一种气电联合系统的协同规划方法
CN107968439B (zh) 基于混合整数线性规划的主动配电网联合优化算法
Hu et al. Transmission expansion planning considering the deployment of energy storage systems
Zhou et al. An optimal network constraint-based joint expansion planning model for modern distribution networks with multi-types intermittent RERs
Gu et al. Joint planning of electrical storage and gas storage in power-gas distribution network considering high-penetration electric vehicle and gas vehicle
Mitridati et al. Power systems flexibility from district heating networks
CN109558990B (zh) 一种基于斯坦纳树模型的配电网防灾骨干网架规划方法
Gan et al. A low-carbon planning method for joint regional-district multi-energy systems: From the perspective of privacy protection
CN107480913A (zh) 一种基于改进灰狼算法的分布式电源选址定容系统及方法
CN113901653A (zh) 一种电-气综合能源系统及其潮流求解方法
CN111724064A (zh) 一种基于改进免疫算法的含储能配电网规划方法
Domínguez et al. Analysing decarbonizing strategies in the European power system applying stochastic dominance constraints
Alayo et al. Economic externalities in transmission network expansion planning
CN112069634B (zh) 一种基于放宽约束的气网扩容规划系统和方法
CN109359830A (zh) 一种水电站电力交易计划电量分解方法
Wei et al. Transmission network planning with N-1 security criterion based on improved multi-objective genetic algorithm
CN117767407A (zh) 一种分布式电源联合可信容量评估方法
CN117114776A (zh) 一种面向省间日前现货交易的量价申报方法
Wang et al. Optimal double Q AC-DC hybrid distribution system planning with explicit topology-variable-based reliability assessment
Fan et al. A post-disaster restoration model for power-gas-transportation distribution networks considering spatial interdependency and energy hubs
Gendreau et al. Multi-period capacity expansion for a local access telecommunications network
Lou et al. A coordinated restoration method of three-phase AC unbalanced distribution network with DC connections and mobile energy storage systems
CN115545340A (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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20220809