CN111724026A - 一种多能源网络与配水网耦合运行的优化方法 - Google Patents

一种多能源网络与配水网耦合运行的优化方法 Download PDF

Info

Publication number
CN111724026A
CN111724026A CN202010377682.3A CN202010377682A CN111724026A CN 111724026 A CN111724026 A CN 111724026A CN 202010377682 A CN202010377682 A CN 202010377682A CN 111724026 A CN111724026 A CN 111724026A
Authority
CN
China
Prior art keywords
network
water
representing
gas
energy
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
CN202010377682.3A
Other languages
English (en)
Other versions
CN111724026B (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.)
North China Electric Power University
Original Assignee
North China Electric Power University
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 North China Electric Power University filed Critical North China Electric Power University
Priority to CN202010377682.3A priority Critical patent/CN111724026B/zh
Publication of CN111724026A publication Critical patent/CN111724026A/zh
Application granted granted Critical
Publication of CN111724026B publication Critical patent/CN111724026B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
    • 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
    • G06Q30/00—Commerce
    • G06Q30/02—Marketing; Price estimation or determination; Fundraising
    • G06Q30/0201—Market modelling; Market analysis; Collecting market data
    • G06Q30/0206—Price or cost determination based on market factors
    • 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
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A20/00—Water conservation; Efficient water supply; Efficient water use
    • Y02A20/20—Controlling water pollution; Waste water treatment
    • Y—GENERAL 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
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS 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/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Strategic Management (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Development Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Business, Economics & Management (AREA)
  • Finance (AREA)
  • Marketing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Accounting & Taxation (AREA)
  • Tourism & Hospitality (AREA)
  • Game Theory and Decision Science (AREA)
  • Health & Medical Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Educational Administration (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种多能源网络与配水网耦合运行的优化方法，首先根据配水网络与多能源系统中的输电网络、输气网络和配热网络之间的耦合关系，搭建一个相互依赖和作用的水‑能源联系运行模型；所搭建的模型为能源和资源的双目标运行模型；采用双目标决策框架来实现能源和资源两个目标之间的公平权衡，将所述水‑能源联系运行模型表示为一个混合整数非线性模型；在利用双目标优化的纳什议价法来处理所述水‑能源联系运行模型中的双目标竞争关系，实现多能源网络与配水网耦合运行的优化。该方法可实现能源生产成本与水资源消耗之间的平衡，为能源和资源的综合运行下的协调优化提出了解决方案。

Description

一种多能源网络与配水网耦合运行的优化方法
技术领域
本发明涉及能源网络研究技术领域，尤其涉及一种多能源网络与配水网耦合运行的优化方法。
背景技术
在过去的几十年中，由于化石能源储量的减少以及人口和经济的增长，水资源和能源作为人类生存和社会发展的基础，引起了工业界和学术界的高度重视，逐渐形成水能联系的概念：water-energy nexus，二者的优化协调对能源系统的优化以及资源的节约都具有较高的研究价值。当下在对水资源与多能流关系的研究中，一些研究学者对配水网与电网之间的关系进行了一定的研究，其中较为普遍的是对最优泵调度问题的研究。但由于不同学科之间存在差异性，目前对配水网和电网之间相互作用的研究仍较浅，仅停留在对水泵运行能耗的优化层面，鲜有将水泵的灵活性应用于电力系统的供能优化研究中，更直观地来说，缺乏跨部门的协调来引导对水和能源的有效利用。
现有技术中研究水泵的调度对配水网络与供电网络的影响较为多见，利用水泵的柔性特征形成水-能综合优化框架；同时由于发电过程中存在一定的耗水量，部分学者从宏观角度分析水和电能的相关性政策对自然环境、社会经济的影响，并进行了水电依赖综合评价；在规划层面也有部分研究结合当地供水压力、发电资产扩张和长期气候变化等因素设计了可再生能源驱动的水能耦合系统框架，为未来水-能场景研究提供了研究支撑。长期以来，能源和资源两个系统大多数是相互独立运行的，将两个系统耦合起来进行优化运行的工作较少且不成熟，现有技术中对能源网络和配水网之间相互作用的研究仍较浅，对能源和资源在耦合运行过程中的成本竞争研究更少，因此该研究亟待进行。
发明内容
本发明的目的是提供一种多能源网络与配水网耦合运行的优化方法，该方法可实现能源生产成本与水资源消耗之间的平衡，为能源和资源的综合运行下的协调优化提出了解决方案。
本发明的目的是通过以下技术方案实现的：
一种多能源网络与配水网耦合运行的优化方法，所述方法包括：
步骤1、根据配水网络与多能源系统中的输电网络、输气网络和配热网络之间的耦合关系，搭建一个相互依赖和作用的水-能源联系运行模型；所搭建的模型为能源和资源的双目标运行模型；
步骤2、采用双目标决策框架来实现能源和资源两个目标之间的公平权衡，将所述水-能源联系运行模型表示为一个混合整数非线性模型；
步骤3、在利用双目标优化的纳什议价法来处理所述水-能源联系运行模型中的双目标竞争关系，实现多能源网络与配水网耦合运行的优化。
由上述本发明提供的技术方案可以看出，上述方法可实现能源生产成本与水资源消耗之间的平衡，为能源和资源的综合运行下的协调优化提出了解决方案，对将来构建智能一体化的社会综合能源与资源网络具有良好的社会意义。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域的普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他附图。
图1为本发明实施例提供的多能源网络与配水网耦合运行的优化方法流程示意图；
图2为本发明实施例所述纳什议价的几何原理示意图。
具体实施方式
下面结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明的保护范围。
下面将结合附图对本发明实施例作进一步地详细描述，如图1所示为本发明实施例提供的多能源网络与配水网耦合运行的优化方法流程示意图，所述方法包括：
步骤1、根据配水网络与多能源系统中的输电网络、输气网络和配热网络之间的耦合关系，搭建一个相互依赖和作用的水-能源联系运行模型；
在该步骤中，配水网络是本模型的核心，各个能源网络的耦合元件将配水网与其他能源系统紧密耦合在一起，例如：配水网的各项服务需要水泵维系，而水泵需要从电网索取电能；输电网中的燃煤机组、燃气机组、CHP机组以及输气网中的电转气设备的运作过程均需要配水网供给水资源；同时配热网中由于连接用户用水，也需要向配水网索取水资源，如此便形成了水-能联系网络。
在本实例中，配水网络与多能源系统中的输电网络、输气网络和配热网络之间的耦合关系具体为：
在输电网络中，将配水网络中的水泵和输气网络中的P2G装置接至电负荷节点处；
在输气网络中，将输电网络中的CHP接至气负荷节点处；
在配热网络中，将CHP设备放至加热站作为热源，并连接有储热罐用于补给水；
在配水网络中，将P2G耗水、火电、风电和CHP机组耗水接至配水网络的负荷节点处，以此构建具有紧密耦合关系的水-能源联系运行模型。
具体实现中，上述所搭建的模型为能源和资源的双目标运行模型，所建立的模型具体为输电网模型、输气网模型、配热网模型和配水网模型共同约束搭建的模型，具体包括：
在输电网模型中，电力系统全网平衡约束如式(1)所示；输电网络发电出力约束、弃风约束分别如式(2)、(3)所示；火电机组的爬坡约束如式(4)所示：
Figure BDA0002480614880000031
Figure BDA0002480614880000032
Figure BDA0002480614880000033
Figure BDA0002480614880000034
其中，Θ表示一系列接于节点e的设备；en/ew/es/ep2g分别表示火电机组节点、风电机组节点、燃气机组节点、电转气设备节点；
Figure BDA0002480614880000035
表示发电机出力；
Figure BDA0002480614880000036
表示管道上的输电功率；
Figure BDA0002480614880000037
表示流入和流出节点e的输电管道；
Figure BDA0002480614880000038
热网和水网水泵耗电量；
Figure BDA0002480614880000039
和
Figure BDA00024806148800000310
分别表示风力发电功率和弃风量；
Figure BDA00024806148800000311
发电机组en的出力上下限；
Figure BDA00024806148800000312
表示发电机组的爬坡约束值；
Figure BDA00024806148800000313
电负荷；pct表示CHP机组出力；
Figure BDA00024806148800000314
表示输电线路功率；θet表示电压相角；
Figure BDA00024806148800000315
表示机组启停状态；
Figure BDA00024806148800000316
表示线路传输功率极限；
加热站处CHP机组的产能关系式为式(5)-(8)所示，其中产出的电功率和热功率之间满足一定的关系，可以通过特定的多边形凸区域的顶点线性组合进行计算：
Figure BDA00024806148800000317
Figure BDA0002480614880000041
Figure BDA0002480614880000042
Figure BDA0002480614880000043
输电网中线路功率限制表示为式(9)：
Figure BDA0002480614880000044
其中，
Figure BDA0002480614880000045
分别为燃煤机组、燃气机组、CHP机组、风力机组、水网和热网水泵、电转气机组以及电力负荷的转移分布因子；
在输气网模型中，为了简便并未考虑输气管道内存气的情况，且压缩机模型的成本函数为线性，输气网节点平衡方程如式(10)所示：
Figure BDA0002480614880000046
其中，
Figure BDA0002480614880000047
λg分别表示电转气机组和燃气机组产气量和耗气量与电能之间的比例；Λ表示一系列接于节点g的设备；
Figure BDA0002480614880000048
表示CHP机组的耗气量；
Figure BDA0002480614880000049
表示气网管道流量；
Figure BDA00024806148800000410
表示储气装置储气量；
Figure BDA00024806148800000411
表示向储气装置注入和从储气装置提取的天然气量；
Figure BDA00024806148800000412
表示从气井取出的天然气量；
Figure BDA00024806148800000413
表示电转气设备的产气量；
Figure BDA00024806148800000414
表示气负荷；
Figure BDA00024806148800000415
表示气井的取气量上下限；pct/hct表示CHP机组的产电/产热功率；
压缩机管道耗气量计算式为式(11)，气井的供气量约束为式(12)，CHP设备的耗气量计算式为式(13)，该式体现了CHP机组在生产电功率和热功率时动态耗气的特性，
Figure BDA00024806148800000416
Figure BDA00024806148800000417
Figure BDA00024806148800000418
其中，
Figure BDA00024806148800000419
表示CHP机组的耗气系数；
输气网中管道压力和气流量的关系如Weymouth方程(14)所示，其中为了简便将气流量的平方值Πgt来代替二次项；
Figure BDA0002480614880000051
对于装有压缩机的管道，管道始端节点
Figure BDA0002480614880000052
的压力应比管道末端节点
Figure BDA0002480614880000053
的压力要小，节点压力约束以及含有压缩机管道的节点压力关系如式(15)和(16)所示：
Figure BDA0002480614880000054
Figure BDA0002480614880000055
天然气系统中的储气装置容量约束如式(17)所示,其中储气装置的天然气注入和提取约束如(18)所示：
Figure BDA0002480614880000056
Figure BDA0002480614880000057
其中，
Figure BDA0002480614880000058
表示接于节点g的储气装置；
Figure BDA0002480614880000059
表示储气装置的进气和出气效率；
Figure BDA00024806148800000510
表示气网管道压力平方值的上下限；
Figure BDA00024806148800000511
表示储气装置注入和提取限制；Πgt表示气网节点气压的平方值；
Figure BDA00024806148800000512
表示气网管道中始节点
Figure BDA00024806148800000513
和末节点
Figure BDA00024806148800000514
气压的平方值；
Figure BDA00024806148800000515
表示气网压缩机压力系数；
在配水网模型中，配水网中的节点平衡约束如式(19)所示，其中Ξ表示一系列接于节点w的设备，管道中的流量约束以及水库取水量约束分别如式(20)和(21)所示：
Figure BDA00024806148800000516
Figure BDA00024806148800000517
Figure BDA00024806148800000518
其中，
Figure BDA00024806148800000519
表示以节点w为管道首节点和末节点的水网管道集合；
Figure BDA00024806148800000520
表示水网负荷；
Figure BDA00024806148800000521
表示水网管道流量；
Figure BDA00024806148800000522
表示出水装置的提取和注入流量；
Figure BDA00024806148800000523
表示水库的取水流量；
Figure BDA00024806148800000524
表示管道流量的上限；
Figure BDA00024806148800000525
表示水库取水量上限；
进一步的，考虑配水网中的水泵能耗与管道流量的关系，普通管道和装有水泵管道的水力特性分别如式(22)和(23)所示，而且本文中采用恒速的水泵模型，水泵的水头压力表达式为式(24)，在式(23)中，假设该管道上的水泵处于打开的状态，那么它会在该管道上产生压头，如上所述；但若该水泵处于关断的状态，则该管道不可流通，整个配水网中节点压头约束为式(25)，水泵的耗能计算式为式(26)：
Figure BDA0002480614880000061
Figure BDA0002480614880000062
Figure BDA0002480614880000063
Figure BDA0002480614880000064
Figure BDA0002480614880000065
其中，
Figure BDA0002480614880000066
表示水泵压力系数；
Figure BDA0002480614880000067
表示管道压力损耗系数；
Figure BDA0002480614880000068
表示水网节点压力上下限；
Figure BDA0002480614880000069
表示水泵效率；
Figure BDA00024806148800000610
表示水泵能耗；
Figure BDA00024806148800000611
表示水泵所在管道流量；
Figure BDA00024806148800000612
表示水泵产生的压头；πwt表示水网节点压力；
在配热网模型中，为了构建与配水网紧密耦合的供暖网络，本文在开式供热系统的基础上，做出了以下合理假设：①该开式供热系统配有热水存储装置、配有淋浴装置的供热用户，以及供热和回热管道；②该热网模型是一个开放的供暖系统，这意味着给水管网和回水管网中的水流可能是不同的；③该供热系统的热源来自于CHP机组，同时在该过程中CHP机组产出的电功率用于输电系统中，而且本模型中的CHP机组的采用燃气供能的形式；④在本模型中，为了降低建模难度，热网模型中采用质调节的方式，固定供热网和回热网的温度，仅对水流变化进行分析。
上述配热网模型中供热管网和回热管网节点的平衡约束分别为式(26)和(27)：
Figure BDA00024806148800000613
Figure BDA00024806148800000614
其中，
Figure BDA00024806148800000615
表示接在换热站he处的热负荷序数和集合；
Figure BDA00024806148800000616
表示接在换热站he处的回热网节点
Figure BDA00024806148800000617
序数和集合；
Figure BDA00024806148800000618
表示接在加热站hs处的回热网节点序数和集合；
Figure BDA00024806148800000619
Figure BDA00024806148800000620
表示接在换热站he处的供热网节点
Figure BDA00024806148800000621
序数和集合；
Figure BDA00024806148800000622
表示接在加热站hs处的供热网节点
Figure BDA00024806148800000623
序数和集合；
Figure BDA00024806148800000624
表示加热站hs处水流提取和注入速率上限；
Figure BDA00024806148800000625
表示加热站hs处蓄热罐的容量上下限；
Figure BDA00024806148800000626
表示加热站处蓄水罐从配水网引入的流量；
Figure BDA0002480614880000071
表示换热站处的流量；
Figure BDA0002480614880000072
表示换热站连接回热网的管道的流量；
Figure BDA0002480614880000073
表示换热站连接供热网管道的流量；
Figure BDA0002480614880000074
表示加热站的流量；
Figure BDA0002480614880000075
表示加热站连接回热网管道的流量；
Figure BDA0002480614880000076
表示加热站连接供热网管道的流量；
Figure BDA0002480614880000077
表示热网负荷流量；
Figure BDA0002480614880000078
表示热网管道流量；
Figure BDA0002480614880000079
表示加热站处蓄水池对其注入的流量；
Figure BDA00024806148800000710
表示加热站流入蓄水池的流量；
加热站中注入和提取的流量约束如式(28)和(29)所示，同时根据式(30)计算出蓄水池的每个时刻的水量：
Figure BDA00024806148800000711
Figure BDA00024806148800000712
Figure BDA00024806148800000713
加热站处安装水泵的管道水力特性由式(31)计算得到，管道的水头增益计算如式(32)所示，加热站与回热网相连接的管道并无安装水泵，因此其水力特性由式(33)计算得到：
Figure BDA00024806148800000714
Figure BDA00024806148800000715
Figure BDA00024806148800000716
供热网中水泵能耗如式(34)所示：
Figure BDA00024806148800000717
本模型中供热网仅在加热站处配有一台水泵，其耗电量以及整个网络的管道流量约束如式(35)和(36)所示：
Figure BDA00024806148800000718
Figure BDA00024806148800000719
其中，
Figure BDA00024806148800000720
表示热网管道压力损耗系数；
Figure BDA00024806148800000721
表示加热站的水泵能耗系数；
Figure BDA00024806148800000722
表示加热站水泵能耗效率；
Figure BDA00024806148800000723
表示热网管道流量上下限；
Figure BDA00024806148800000724
表示加热站的水泵能耗限制；
Figure BDA00024806148800000725
表示加热站水泵产生的压力；
Figure BDA00024806148800000726
表示供热网管道始节点和末节点的压力；
加热站的供热平衡约束式为式(37)所示：
Figure BDA0002480614880000081
其中，Tin表示从配水网取出后注入至加热站的水流的温度；
同时换热站处流量以及换热平衡约束为式(38)-(39)所示：
Figure BDA0002480614880000082
Figure BDA0002480614880000083
其中，
Figure BDA0002480614880000084
表示CHP机组运行区域顶点序数和集合；
Figure BDA0002480614880000085
表示CHP机组运行区域顶点的线性组合常数；c表示水的比热容；
Figure BDA0002480614880000086
表示CHP机组运行区域顶点值；
Figure BDA0002480614880000087
表示CHP机组产热功率；
Figure BDA0002480614880000088
表示换热站的交换热功率；
Figure BDA0002480614880000089
表示表示CHP机组的产电功率；
Figure BDA00024806148800000810
表示蓄热罐的温度；
Figure BDA00024806148800000811
表示供热网温度；
Figure BDA00024806148800000812
表示回热网温度；
上述每个模型之间并非独立，他们之间通过前文所述的耦合元件如P2G、CHP等进行耦合。
步骤2、采用双目标决策框架(下称P-1)来实现能源和资源两个目标之间的公平权衡，将所述水-能源联系运行模型表示为一个混合整数非线性模型；
在该步骤中，水-能源联系运行模型的目标是尽量减少能源生产成本和水消耗成本，由于这两个目标之间可能存在冲突，而且它们的优先选择顺序也可能不同，因此本实施采用双目标决策框架(下称P-1)来实现能源和资源两个目标之间的公平权衡，该双目标决策框架P-1表示为：
Figure BDA00024806148800000813
其中，Coste,Costw分别表示能源生产成本和水资源成本；x和
Figure BDA00024806148800000819
分别是决策向量及其可行域；
Coste,Costw的详细表达式为下式(41)和(42)：
Figure BDA00024806148800000814
Figure BDA00024806148800000815
其中，
Figure BDA00024806148800000816
和
Figure BDA00024806148800000817
分别表示发电成本系数和气井气价；
Figure BDA00024806148800000818
表示了水资源的价格；
Figure BDA0002480614880000091
和
Figure BDA0002480614880000092
分别表示储气装置中充气/放气成本和储水装置中蓄水和抽水成本。
步骤3、在利用双目标优化的纳什议价法来处理所述水-能源联系运行模型中的双目标竞争关系，实现多能源网络与配水网耦合运行的优化。
在该步骤中，所述双目标优化的纳什议价模型表示为：
Figure BDA0002480614880000093
其中，不同成本的计算过程为：
Figure BDA0002480614880000094
其中，LP(0)和LP(1)可通过式(45)进行计算，分别表示以水成本、能源成本为单目标时的优化结果；
如图2所示为本发明实施例所述纳什议价的几何原理示意图，纳什议价过程可以被看作是在两个目标的帕累托前沿找到一个点，从而最大化“交易矩形”的面积；
其中，帕累托前沿通过两个目标函数加权而得，具体为：
Figure BDA0002480614880000095
通过仿真结果表明，与单目标成本最小化方案相比，纳什议价方案能使不同子系统之间的交互更加均衡，使水-能关系更不易受到设备停机或故障的影响。
下面以具体的实例对上述方法的实施过程进行详细说明，本实例所引入的算例仅仅是示例性的，而不是为了限制本发明的范围及应用，本实例所引入的算例模型由39节点电力网络、20节点天然气网络、13节点供热网络以及15节点配水网耦合而成，在MATLAB中运用IPOTP工具包进行仿真求解，这四个网络相互间通过耦合元件进行连接。
为了探索能源和资源两个目标之间在优化过程中的竞争关系，本实例探索了六个场景下的资源和能源的成本，其中Coste表示煤炭和天然气消费的成本，Costw表示水资源成本。在算例1和算例2中，仅有一个目标被最小化，即在此场景下求得的1.1294×107和3.3110×105分别为Coste和Costw的下界。
除了采用纳什议价法(算例6)外，加权方法也是处理多目标竞争情况的一种常用方法，探究在三组权重系数下的算例3、算例4和算例5的仿真结果的不同，如下表1所示：
表1目标成本对比
Figure BDA0002480614880000101
Costw的权系数增大时对应的Costw的成本相应减小，符合常理。但是权重系数的选择是主观的。在算例6中，两个目标的值与算例4非常接近，说明了在加权方法中选择权值的困难，特别是当目标的数量级不同时；尽管算例3提供了两个目标的最小和，但是并不是一个公平的解决方案。
在算例3中，Coste和Costw与它们的下界的相对偏差率分别为0％和$82.2％，因为Coste的数量级相对较高；然而，在本申请所提供的纳什讨价还价算例6中，Coste和Costw与它们的下界的相对偏差率更均匀，分别为6.85％和11.99％，这表明了所提出的方法讨价还价结果的有效性。
值得注意的是，本发明实施例中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明披露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求书的保护范围为准。

Claims (5)

1.一种多能源网络与配水网耦合运行的优化方法，其特征在于，所述方法包括：
步骤1、根据配水网络与多能源系统中的输电网络、输气网络和配热网络之间的耦合关系，搭建一个相互依赖和作用的水-能源联系运行模型；所搭建的模型为能源和资源的双目标运行模型；
步骤2、采用双目标决策框架来实现能源和资源两个目标之间的公平权衡，将所述水-能源联系运行模型表示为一个混合整数非线性模型；
步骤3、在利用双目标优化的纳什议价法来处理所述水-能源联系运行模型中的双目标竞争关系，实现多能源网络与配水网耦合运行的优化。
2.根据权利要求1所述多能源网络与配水网耦合运行的优化方法，其特征在于，在步骤1中，所述配水网络与多能源系统中的输电网络、输气网络和配热网络之间的耦合关系具体为：
在输电网络中，将配水网络中的水泵和输气网络中的P2G装置接至电负荷节点处；
在输气网络中，将输电网络中的CHP接至气负荷节点处；
在配热网络中，将CHP设备放至加热站作为热源，并连接有储热罐用于补给水；
在配水网络中，将P2G耗水、火电、风电和CHP机组耗水接至配水网络的负荷节点处，以此构建具有紧密耦合关系的水-能源联系运行模型。
3.根据权利要求1所述多能源网络与配水网耦合运行的优化方法，其特征在于，在步骤1中，所建立的水-能源联系运行模型具体为输电网模型、输气网模型、配热网模型和配水网模型共同约束搭建的模型，具体包括：
在输电网模型中，电力系统全网平衡约束如式(1)所示；输电网络发电出力约束、弃风约束分别如式(2)、(3)所示；火电机组的爬坡约束如式(4)所示：
Figure FDA0002480614870000011
Figure FDA0002480614870000012
Figure FDA0002480614870000013
Figure FDA0002480614870000014
其中，Θ表示一系列接于节点e的设备；en/ew/es/ep2g分别表示火电机组节点、风电机组节点、燃气机组节点、电转气设备节点；
Figure FDA0002480614870000015
表示发电机出力；
Figure FDA0002480614870000016
表示管道上的输电功率；
Figure FDA0002480614870000021
表示流入和流出节点e的输电管道；
Figure FDA0002480614870000022
热网和水网水泵耗电量；
Figure FDA0002480614870000023
和
Figure FDA0002480614870000024
分别表示风力发电功率和弃风量；
Figure FDA0002480614870000025
发电机组en的出力上下限；
Figure FDA0002480614870000026
表示发电机组的爬坡约束值；
Figure FDA0002480614870000027
电负荷；pct表示CHP机组出力；
Figure FDA0002480614870000028
表示输电线路功率；θet表示电压相角；
Figure FDA0002480614870000029
表示机组启停状态；
Figure FDA00024806148700000210
表示线路传输功率极限；
加热站处CHP机组的产能关系式为式(5)-(8)所示：
Figure FDA00024806148700000211
Figure FDA00024806148700000212
Figure FDA00024806148700000213
Figure FDA00024806148700000214
输电网中线路功率限制表示为式(9)：
Figure FDA00024806148700000215
其中，
Figure FDA00024806148700000216
分别为燃煤机组、燃气机组、CHP机组、风力机组、水网和热网水泵、电转气机组以及电力负荷的转移分布因子；
在输气网模型中，输气网节点平衡方程如式(10)所示：
Figure FDA00024806148700000217
其中，
Figure FDA00024806148700000218
λg分别表示电转气机组和燃气机组产气量和耗气量与电能之间的比例；Λ表示一系列接于节点g的设备；
Figure FDA00024806148700000219
表示CHP机组的耗气量；
Figure FDA00024806148700000220
表示气网管道流量；
Figure FDA00024806148700000221
表示储气装置储气量；
Figure FDA00024806148700000222
表示向储气装置注入和从储气装置提取的天然气量；
Figure FDA00024806148700000223
表示从气井取出的天然气量；
Figure FDA00024806148700000224
表示电转气设备的产气量；
Figure FDA00024806148700000225
表示气负荷；
Figure FDA00024806148700000226
表示气井的取气量上下限；pct/hct表示CHP机组的产电/产热功率；
压缩机管道耗气量计算式为式(11)，气井的供气量约束为式(12)，CHP设备的耗气量计算式为式(13)：
Figure FDA00024806148700000227
Figure FDA0002480614870000031
Figure FDA0002480614870000032
其中，
Figure FDA0002480614870000033
表示CHP机组的耗气系数；
输气网中管道压力和气流量的关系如Weymouth方程(14)所示，其中为了简便将气流量的平方值Πgt来代替二次项；
Figure FDA0002480614870000034
对于装有压缩机的管道，管道始端节点
Figure FDA0002480614870000035
的压力应比管道末端节点
Figure FDA0002480614870000036
的压力要小，节点压力约束以及含有压缩机管道的节点压力关系如式(15)和(16)所示：
Figure FDA0002480614870000037
Figure FDA0002480614870000038
天然气系统中的储气装置容量约束如式(17)所示,其中储气装置的天然气注入和提取约束如(18)所示：
Figure FDA0002480614870000039
Figure FDA00024806148700000310
其中，
Figure FDA00024806148700000311
表示接于节点g的储气装置；
Figure FDA00024806148700000312
表示储气装置的进气和出气效率；
Figure FDA00024806148700000313
表示气网管道压力平方值的上下限；
Figure FDA00024806148700000314
表示储气装置注入和提取限制；Πgt表示气网节点气压的平方值；
Figure FDA00024806148700000315
表示气网管道中始节点
Figure FDA00024806148700000316
和末节点
Figure FDA00024806148700000317
气压的平方值；
Figure FDA00024806148700000318
表示气网压缩机压力系数；
在配水网模型中，配水网中的节点平衡约束如式(19)所示，其中Ξ表示一系列接于节点w的设备，管道中的流量约束以及水库取水量约束分别如式(20)和(21)所示：
Figure FDA00024806148700000319
Figure FDA00024806148700000320
Figure FDA00024806148700000321
其中，
Figure FDA0002480614870000041
表示以节点w为管道首节点和末节点的水网管道集合；
Figure FDA0002480614870000042
表示水网负荷；
Figure FDA0002480614870000043
表示水网管道流量；
Figure FDA0002480614870000044
表示出水装置的提取和注入流量；
Figure FDA0002480614870000045
表示水库的取水流量；
Figure FDA0002480614870000046
表示管道流量的上限；
Figure FDA0002480614870000047
表示水库取水量上限；
进一步的，考虑配水网中的水泵能耗与管道流量的关系，普通管道和装有水泵管道的水力特性分别如式(22)和(23)所示，水泵的水头压力表达式为式(24)，在式(23)中，假设该管道上的水泵处于打开的状态，那么它会在该管道上产生压头；但若该水泵处于关断的状态，则管道不可流通，整个配水网中节点压头约束为式(25)，水泵的耗能计算式为式(26)：
Figure FDA0002480614870000048
Figure FDA0002480614870000049
Figure FDA00024806148700000410
Figure FDA00024806148700000411
Figure FDA00024806148700000412
其中，
Figure FDA00024806148700000413
表示水泵压力系数；
Figure FDA00024806148700000414
表示管道压力损耗系数；
Figure FDA00024806148700000415
表示水网节点压力上下限；
Figure FDA00024806148700000416
表示水泵效率；
Figure FDA00024806148700000417
表示水泵能耗；
Figure FDA00024806148700000418
表示水泵所在管道流量；
Figure FDA00024806148700000419
表示水泵产生的压头；πwt表示水网节点压力；
在配热网模型中，供热管网和回热管网节点的平衡约束分别为式(26)和(27)：
Figure FDA00024806148700000420
Figure FDA00024806148700000421
其中，
Figure FDA00024806148700000422
表示接在换热站he处的热负荷序数和集合；
Figure FDA00024806148700000423
表示接在换热站he处的回热网节点
Figure FDA00024806148700000424
序数和集合；
Figure FDA00024806148700000425
表示接在加热站hs处的回热网节点序数和集合；
Figure FDA00024806148700000426
表示接在换热站he处的供热网节点
Figure FDA00024806148700000427
序数和集合；
Figure FDA00024806148700000428
表示接在加热站hs处的供热网节点
Figure FDA00024806148700000429
序数和集合；
Figure FDA00024806148700000430
表示加热站hs处水流提取和注入速率上限；
Figure FDA0002480614870000051
表示加热站hs处蓄热罐的容量上下限；
Figure FDA0002480614870000052
表示加热站处蓄水罐从配水网引入的流量；
Figure FDA0002480614870000053
表示换热站处的流量；
Figure FDA0002480614870000054
表示换热站连接回热网的管道的流量；
Figure FDA0002480614870000055
表示换热站连接供热网管道的流量；
Figure FDA0002480614870000056
表示加热站的流量；
Figure FDA0002480614870000057
表示加热站连接回热网管道的流量；
Figure FDA0002480614870000058
表示加热站连接供热网管道的流量；
Figure FDA0002480614870000059
表示热网负荷流量；
Figure FDA00024806148700000510
表示热网管道流量；
Figure FDA00024806148700000511
表示加热站处蓄水池对其注入的流量；
Figure FDA00024806148700000512
表示加热站流入蓄水池的流量；
加热站中注入和提取的流量约束如式(28)和(29)所示，同时根据式(30)计算出蓄水池的每个时刻的水量：
Figure FDA00024806148700000513
Figure FDA00024806148700000514
Figure FDA00024806148700000515
加热站处安装水泵的管道水力特性由式(31)计算得到，管道的水头增益计算如式(32)所示，加热站与回热网相连接的管道并无安装水泵，因此其水力特性由式(33)计算得到：
Figure FDA00024806148700000516
Figure FDA00024806148700000517
Figure FDA00024806148700000518
供热网中水泵能耗如式(34)所示：
Figure FDA00024806148700000519
本模型中供热网仅在加热站处配有一台水泵，其耗电量以及整个网络的管道流量约束如式(35)和(36)所示：
Figure FDA00024806148700000520
Figure FDA00024806148700000521
其中，
Figure FDA00024806148700000522
表示热网管道压力损耗系数；
Figure FDA00024806148700000523
表示加热站的水泵能耗系数；
Figure FDA00024806148700000524
表示加热站水泵能耗效率；
Figure FDA00024806148700000525
表示热网管道流量上下限；
Figure FDA00024806148700000526
表示加热站的水泵能耗限制；
Figure FDA0002480614870000061
表示加热站水泵产生的压力；
Figure FDA0002480614870000062
表示供热网管道始节点和末节点的压力；
加热站的供热平衡约束式为式(37)所示：
Figure FDA0002480614870000063
其中，Tin表示从配水网取出后注入至加热站的水流的温度；
同时换热站处流量以及换热平衡约束为式(38)-(39)所示：
Figure FDA0002480614870000064
Figure FDA0002480614870000065
其中，
Figure FDA0002480614870000066
表示CHP机组运行区域顶点序数和集合；
Figure FDA0002480614870000067
表示CHP机组运行区域顶点的线性组合常数；c表示水的比热容；
Figure FDA0002480614870000068
表示CHP机组运行区域顶点值；
Figure FDA0002480614870000069
表示CHP机组产热功率；
Figure FDA00024806148700000610
表示换热站的交换热功率；
Figure FDA00024806148700000611
表示表示CHP机组的产电功率；
Figure FDA00024806148700000612
表示蓄热罐的温度；
Figure FDA00024806148700000613
表示供热网温度；
Figure FDA00024806148700000614
表示回热网温度；
上述每个模型之间并非独立，他们之间通过耦合元件进行耦合。
4.根据权利要求1所述多能源网络与配水网耦合运行的优化方法，其特征在于，在步骤2中，所述双目标决策框架P-1表示为：
Figure FDA00024806148700000615
其中，Coste,Costw分别表示能源生产成本和水资源成本；x和
Figure FDA00024806148700000616
分别是决策向量及其可行域；
Coste,Costw的详细表达式为下式(41)和(42)：
Figure FDA00024806148700000617
Figure FDA00024806148700000618
其中，
Figure FDA00024806148700000619
和
Figure FDA00024806148700000620
分别表示发电成本系数和气井气价；
Figure FDA00024806148700000621
表示了水资源的价格；
Figure FDA00024806148700000622
和
Figure FDA00024806148700000623
分别表示储气装置中充气/放气成本和储水装置中蓄水和抽水成本。
5.根据权利要求1所述多能源网络与配水网耦合运行的优化方法，其特征在于，在步骤3中，所述双目标优化的纳什议价模型表示为：
Figure FDA0002480614870000071
其中，不同成本的计算过程为：
Figure FDA0002480614870000072
其中LP(0)和LP(1)通过式(45)进行计算，分别表示以水成本、能源成本为单目标时的优化结果；
上述纳什议价过程是在两个目标的帕累托前沿找到一个点，从而最大化交易矩形的面积；其中，帕累托前沿通过两个目标函数加权而得，具体为：
Figure FDA0002480614870000073
CN202010377682.3A 2020-05-07 2020-05-07 一种多能源网络与配水网耦合运行的优化方法 Active CN111724026B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010377682.3A CN111724026B (zh) 2020-05-07 2020-05-07 一种多能源网络与配水网耦合运行的优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010377682.3A CN111724026B (zh) 2020-05-07 2020-05-07 一种多能源网络与配水网耦合运行的优化方法

Publications (2)

Publication Number Publication Date
CN111724026A true CN111724026A (zh) 2020-09-29
CN111724026B CN111724026B (zh) 2023-06-16

Family

ID=72564270

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010377682.3A Active CN111724026B (zh) 2020-05-07 2020-05-07 一种多能源网络与配水网耦合运行的优化方法

Country Status (1)

Country Link
CN (1) CN111724026B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112381267A (zh) * 2020-10-28 2021-02-19 武汉供电设计院有限公司 一种面向工业园区的供水-供能耦合规划方法
CN114065668A (zh) * 2021-11-25 2022-02-18 重庆大学 一种基于图论的配水系统沿线流量和水头压力的量化计算方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109978625A (zh) * 2019-03-28 2019-07-05 河海大学 一种计及电热气网络的综合能源系统多目标运行优化方法
CN110457745A (zh) * 2019-06-27 2019-11-15 华北电力大学 一种具有紧密耦合关系的水-能联系网络模型的建立方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109978625A (zh) * 2019-03-28 2019-07-05 河海大学 一种计及电热气网络的综合能源系统多目标运行优化方法
CN110457745A (zh) * 2019-06-27 2019-11-15 华北电力大学 一种具有紧密耦合关系的水-能联系网络模型的建立方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112381267A (zh) * 2020-10-28 2021-02-19 武汉供电设计院有限公司 一种面向工业园区的供水-供能耦合规划方法
CN114065668A (zh) * 2021-11-25 2022-02-18 重庆大学 一种基于图论的配水系统沿线流量和水头压力的量化计算方法
CN114065668B (zh) * 2021-11-25 2024-04-05 重庆大学 一种基于图论的配水系统沿线流量和水头压力的量化计算方法

Also Published As

Publication number Publication date
CN111724026B (zh) 2023-06-16

Similar Documents

Publication Publication Date Title
CN108154309B (zh) 计及冷热电多负荷动态响应的能源互联网经济调度方法
CN110163443B (zh) 考虑电-气综合需求响应的天然气调压站微能网优化调度方法
Cheng et al. Optimization of multi-carrier energy system based on new operation mechanism modelling of power-to-gas integrated with CO2-based electrothermal energy storage
Ju et al. Multi-objective electro-thermal coupling scheduling model for a hybrid energy system comprising wind power plant, conventional gas turbine, and regenerative electric boiler, considering uncertainty and demand response
CN111738502A (zh) 促进富余风电消纳的多能互补系统需求响应运行优化方法
CN108229025A (zh) 一种冷热电联供型多微网主动配电系统经济优化调度方法
Awad et al. Integrated optimal power flow for electric power and heat in a microgrid
CN106532709A (zh) 一种含配电网重构的区域综合能源系统最优潮流计算方法
CN109523065A (zh) 一种基于改进量子粒子群算法的微能源网优化调度方法
CN108053104A (zh) 区域电力-天然气-热力综合能源系统能量流优化方法
CN112035984B (zh) 一种电-气-储区域综合能源系统协同规划方法
CN111967659B (zh) 一种基于光伏消纳的区域综合能源系统配置优化方法
CN114595868A (zh) 一种综合能源系统源网荷储协同规划方法及系统
CN105955931A (zh) 面向高密度分布式光伏消纳的区域能源网络优化调度方法
CN113131513B (zh) 考虑碳排放的电、热、气转换优化运行配置方法和存储介质
CN113487188A (zh) 考虑电气联合价格引导机制的综合能源系统优化调度方法
CN114077934A (zh) 一种综合能源微网互联系统及其调度方法
CN111724026B (zh) 一种多能源网络与配水网耦合运行的优化方法
Qin et al. An adaptive piecewise linearized weighted directed graph for the modeling and operational optimization of integrated energy systems
CN114066204A (zh) 一种综合能源系统的一体化优化规划与运行方法及装置
Tan et al. Exploiting integrated flexibility from a local smart energy hub
CN114066056A (zh) 考虑热电协同综合能源系统灵活性的优化调度方法及系统
CN110119855A (zh) 一种气电耦合系统运行优化方法
CN112214883B (zh) 一种微能源网综合能源系统运行优化方法和系统
CN110457745B (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