CN117371154B - 含m3c变频站选址定容的工分频互联系统电网规划系统 - Google Patents
含m3c变频站选址定容的工分频互联系统电网规划系统 Download PDFInfo
- Publication number
- CN117371154B CN117371154B CN202311659037.0A CN202311659037A CN117371154B CN 117371154 B CN117371154 B CN 117371154B CN 202311659037 A CN202311659037 A CN 202311659037A CN 117371154 B CN117371154 B CN 117371154B
- Authority
- CN
- China
- Prior art keywords
- power
- frequency conversion
- conversion station
- node
- frequency division
- 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
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 95
- 238000013439 planning Methods 0.000 title claims abstract description 57
- 230000005540 biological transmission Effects 0.000 claims abstract description 43
- 238000000034 method Methods 0.000 claims description 49
- 230000008569 process Effects 0.000 claims description 16
- 238000012423 maintenance Methods 0.000 claims description 8
- 238000011084 recovery Methods 0.000 claims description 4
- 238000004513 sizing Methods 0.000 claims 6
- 235000009508 confectionery Nutrition 0.000 claims 1
- 238000011161 development Methods 0.000 abstract description 9
- 238000010276 construction Methods 0.000 abstract description 8
- 238000005516 engineering process Methods 0.000 abstract description 7
- 238000004364 calculation method Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 239000003990 capacitor Substances 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
- 238000012384 transportation and delivery Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 238000010792 warming Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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
- 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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/36—Arrangements for transfer of electric power between ac networks via a high-tension dc link
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/04—Power grid distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Business, Economics & Management (AREA)
- General Engineering & Computer Science (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Mathematical Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Pure & Applied Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明提供了一种含M3C变频站选址定容的工分频互联系统电网规划系统,属于交直流电网规划技术领域。该模型适用于目前分频输电技术的发展需求,其目标函数针对M3C变频站的并网点和建设容量进行规划,得到考虑M3C变频站选址定容的工分频互联系统电网规划模型。
Description
技术领域
本发明属于交直流电网规划技术领域,具体涉及一种含M3C变频站选址定容的工分频互联系统电网规划系统。
背景技术
在全球变暖问题和节能环保倡议早已被纳入各国发展计划的今天,用可再生能源取代化石能源一直是各国的努力目标。可再生能源包括水能、风能和太阳能等。其中风能被人类利用的历史非常悠久,最早甚至可以追溯到公元前。在二十一世纪,风能在我国和全世界已经有了非常大规模的利用。从2015年至2021年,我国海上风电装机容量保持着年均近60%的增长率,2021年中国海上风电新增并网容量达169万千瓦,累计装机容量达278万千瓦,在风电总装机中占比达到35.5%,成为世界第一大海上风电装机量国。相对于陆上风电,海上风电具有发电利用效率高、不占用土地资源、适宜大规模开发、风机水路运输方便、靠近沿海电力负荷中心等优势。我国海上风能资源丰富,大部分近海海域90米高度年平均风速在7~8.5米/秒之间,具备较好的风能资源条件,适合大规模开发建设海上风电场。
目前海上风电并网主要有三种技术:高压交流输电(high voltage alternatingcurrent,HVAC)、高压直流输电(high voltage direct current,HVDC)以及分频输电(fractional frequency transmission system,FFTS)技术。分频输电方式又称为低频输电方式,最早是由西安交通大学电气工程学院王锡凡教授于1994年针对水电开发提出的一种新型输电方式。FFTS早期研究利用倍频变压器的磁饱和原理实现分频侧与工频电网的互联,低频电网频率为工频电网的三分之一。相比于高压交流输电,FFTS技术由于输电频率的降低,海底电缆的容性充电电流会大大降低,线路传输容量相比工频线路会成倍提高;相比高压直流输电,分频输电不需要建立海上换流站,建设投资和运行维护投资会大大减少,而且其作为交流输电类型比直流具有更强的组网能力。综上,分频风力发/输电系统已成为海上风力研究的一个新浪潮。
电网规划是在机组组合的基础上,对候选线路进行优化选择的问题,优化目标是线路投资成本和系统功率失衡成本综合最小。目前暂无针对工分频互联系统的电网规划模型,但针对新能源接入的交直流混合电网规划问题,现有技术提出一种针对新能源通过高压直流线路送出的场景,对直流送出线路和原交流网架进行规划的方案。
然而,该方案所提出的考虑新能源高压直流接入的交直流电网规划模型针对的场景为新能源高压直流接入,由于直流输电和分频输电在规划、运行和控制上存在较大差异,并不能适用于目前分频输电技术的发展需求,且其在直流送出线路与主网并网时,采用的是固定并网点的方式,没有对并网点位置进行优化。
发明内容
有鉴于此,本发明旨在针对目前分频输电技术的发展需求,提出一套工分频互联系统电网规划模型,该模型针对M3C变频站的并网点和建设容量进行规划,提供一种考虑M3C变频站选址定容的工分频互联系统电网规划模型。
为了解决上述技术问题,本发明提供以下技术方案:
含M3C变频站选址定容的工分频互联系统电网规划系统,工分频互联系统电网规划模型的目标函数,具体如下:
;
式中,是候选M3C变频站编号下标,是候选M3C变频站集合,是M3C变频站规格下标,是M3C变频站规格集合,是M3C变频站第种规格等年值,是变频站第种规格新建状态二进制变量,是线路下标,是候选线路集合,是线路投资等年值,是线路新建状态二进制变量,是场景下标,是场景集合,是场景年等效小时数,是时间下标,是时间集合,是场景第小时总惩罚费用;
该目标函数用于在工分频互联系统电网规划中,在满足负荷需求的前提下,线路一次投资、维护和损耗费用最小。
进一步的,在目标函数中,等年值的折算方式按下式计算:
;
式中,是总费用的等年值,是一次投资费用,是贴现率,是回收年限,是年损耗费用,是年维护费用。
进一步的,在目标函数中,总惩罚费用按下式计算:
;
其中,
;
式中,是场景第小时功率不平衡惩罚费用,是节点下标,是节点集合,分别为有功功率和无功功率不平衡惩罚单位费用,分别为场景第小时节点的有功功率超负荷、欠负荷松弛因子和无功功率超负荷、欠负荷松弛因子;
;
式中,是场景第小时新能源弃用惩罚费用,是新能源机组下标,是新能源机组集合,是新能源弃用惩罚单位费用,是场景第小时新能源机组有功功率弃用因子;
;
式中,是场景第小时电压偏移惩罚费用,是电压偏移惩罚单位费用,分别是场景第小时节点电压上偏移因子和下偏移因子。
进一步的,工分频互联系统电网规划模型还包括约束条件,约束条件具体包括:
机组出力约束、节点功率平衡约束、线路传输容量约束、线路潮流约束、M3C变频站稳态潮流约束以及惩罚项约束。
进一步的,线路潮流约束基于线性化的支路功率方程确定,线性化的支路功率方程的线性化过程如下:
获取节点流向节点的功率方程,如下:
;
式中,是支路传输有功、无功功率,是节点和节点的电压幅值,是节点和节点的电压相角,是支路阻抗参数;
将支路阻抗转换为支路导纳,得到支路功率方程,如下:
;
;
其中是支路导纳参数;
将节点电压幅值和相角作为独立变量在处对上式进行泰勒展开,得到线性化支路功率方程,如下:
。
进一步的,线路潮流约束,具体如下:
已建线路的潮流方程约束为:
;
式中,和是场景第小时线路传输有功功率变量和无功功率变量,分别为线路的电导值和电纳值,下标分别表示线路的起始节点和终止节点,为已建线路集合;
候选线路由二进制变量决定是否建设,其潮流方程约束如下所示:
;
;
式中,为引入的大值。
进一步的,M3C变频站稳态潮流约束基于M3C变频站稳态简化模型得到,在M3C变频站稳态简化模型中,M3C变频站提供的无功支撑受到工/分频侧换流变压器容量的限制,如下式:
;
式中,分别为M3C工频侧网侧有功功率和无功功率,分别为M3C分频侧有功功率和无功功率,分别为工/分频侧换流变压器的容量。
进一步的,M3C变频站稳态潮流约束,具体如下:
M3C变频站传输的有功功率约束:
;
式中,为M3C变频站的有功功率,为M3C变频站第种规格容量,下标中分别表示变频站的工频并网点和分频并网点,为M3C变频站合集;
M3C变频站提供的无功功率约束:
;
式中,为M3C变频站的无功功率,为M3C变频站第种规格换流变压器容量。
进一步的,M3C容量型号的二进制变量满足如下约束:
。
进一步的,在惩罚项约束中,节点电压偏移惩罚约束,如下:
;
式中,U为节点电压幅值,下标i为节点编号,和为节点在场景第小时的上偏移因子和下偏移因子,N为节点集合。
综上,本发明提供了一种含M3C变频站选址定容的工分频互联系统电网规划系统,该模型适用于目前分频输电技术的发展需求,其目标函数针对M3C变频站的并网点和建设容量进行规划,得到考虑M3C变频站选址定容的工分频互联系统电网规划模型。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1为本发明实施例提供的节点流向节点的功率示意图;
图2为本发明实施例提供的背靠背VSC换流器的结构示意图;
图3为本发明实施例提供的一算例下的电网示意图;
图4为本发明实施例提供的一算例下拓展电网示意图。
具体实施方式
为使得本发明的目的、特征、优点能够更加的明显和易懂,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,下面所描述的实施例仅仅是本发明一部分实施例,而非全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
以下首先对现有技术进行详细介绍。
前述提到现有技术提供一种针对新能源通过高压直流线路送出的场景,对直流送出线路和原交流网架进行规划的方案。在该方案中,首先由于交流潮流的非凸问题,现有求解算法并不总是能保证模型收敛和结果的最优,因此该方案使用二阶圆锥松弛对其进行凸化处理;针对交直流混合网络的耦合问题,该方案对换流站中的换流变压器、滤波器和换流器等核心元件利用凸化后的网络参数进行了重新建模;最后与常规电网规划约束条件一起,构成了考虑新能源高压直流接入的交直流混合电网规划模型。
然而,上述现有技术所提出的考虑新能源高压直流接入的交直流电网规划模型存在以下问题:
1)现有技术针对的场景为新能源高压直流接入,由于直流输电和分频输电在规划、运行和控制上存在较大差异,且现阶段分频输电的核心设备——M3C变频站和柔性直流输电的核心设备——VSC换流器在电路结构、数学模型和运行控制策略上完全不同,因此针对目前分频输电技术的发展需求,需要提出一套工分频互联系统电网规划模型。
2)现有技术在处理交流潮流方程非凸问题时,引入了中间变量来进行凸化处理,此方法虽然比利用直流潮流法计算精确,却也是得模型更加复杂,且后续还需要针对提出的中间变量来对VSC换流器进行重新建模。本发明利用线性化思路,对支路功率方程形式的潮流方程进行线性化处理,在不引入中间变量的基础上对非凸的交流潮流方程进行线性化处理,既比直流潮流法更加精确,又没有引入文章所述的中间变量。
3)现有技术在直流送出线路与主网并网时,采用的是固定并网点的方式,没有对并网点位置进行优化。因此本发明针对M3C变频站的并网点和建设容量进行规划,建立考虑M3C变频站选址定容的工分频互联系统电网规划模型。
基于此,以下对本发明的一种含M3C变频站选址定容的工分频互联系统电网规划系统的实施例进行详细介绍。
本实施例提供一种含M3C变频站选址定容的工分频互联系统电网规划系统,其目标函数,具体如下:
;
式中,是候选M3C变频站编号下标,是候选M3C变频站集合,是M3C变频站规格下标,是M3C变频站规格集合,是M3C变频站第种规格等年值,是变频站第种规格新建状态二进制变量,是线路下标,是候选线路集合,是线路投资等年值,是线路新建状态二进制变量,是场景下标,是场景集合,是场景年等效小时数,是时间下标,是时间集合,是场景第小时总惩罚费用;
该目标函数用于在工分频互联系统电网规划中,在满足负荷需求的前提下,线路一次投资、维护和损耗费用最小。
以下结合潮流方程线性化过程、适用规划模型的M3C变频站稳态简化模型、规划模型目标函数和规划模型约束条件4个方面对上述工分频互联系统电网规划模型进行介绍。
1)潮流方程线性化
根据基尔霍夫电流定律和基尔霍夫电压定律,可以得到任意电路各个节点电压和各支路电流的相互关系,考虑支路功率形式的潮流方程。具体分析如下:
考虑如图1所示线路,可以得到从节点流向节点的功率为:
;
其中是支路传输有功、无功功率,是节点和节点的电压幅值,是节点和节点的电压相角,。将支路阻抗转换为支路导纳,可以得到支路功率方程如下所示:
;
;
其中是支路导纳参数。以为独立变量在处对上式进行泰勒展开,得到如下所示的线性化支路功率方程:
;
由此可以得到线性化后的支路功率方程,作为规划模型的潮流方程约束。
2)适用规划模型的M3C变频站稳态简化模型
模块化多电平矩阵式换流器电路稳态运行时可等效为如图2所示背靠背VSC换流器,图中引入了一虚拟直流母线电压。
图2中分别为工频侧、分频侧换流变压器等效阻抗,分别为工频侧、分频侧换流变压器变比,分别为工频侧、分频侧桥臂等效阻抗,与分别为M3C工频侧网侧、等效阀侧有功功率和无功功率,与为分频侧对应功率,分别为工/分频侧换流变压器一次侧电压,分别为工/分频侧换流变压器二次侧电压,分别为变频器工/分频侧等效阀电压。
将桥臂阻抗与换流变压器阻抗合并,即,可由图2可以推导得到M3C变频站稳态运行时注入工频、分频侧电网的功率如下:
;
实验表明,在M3C触发策略确定后,其阀损耗可以等效为下式:
;
其中是由电力电子器件参数确定的损耗率,由制造厂家提供。由图2和上式所述阀损耗表达式可以得到M3C工分频两侧功率耦合关系:
;
M3C变频站提供的无功支撑受到工/分频侧换流变压器容量的限制,如下所示:
;
3)规划模型目标函数
考虑经济性最优,即在满足负荷需求的前提下,线路一次投资、维护和损耗费用最小。模型目标函数如下:
;
上式中是候选M3C变频站编号下标,是候选M3C变频站集合,是M3C变频站规格下标,是M3C变频站规格集合,是M3C变频站第种规格等年值,是变频站第种规格新建状态二进制变量,是线路下标,是候选线路集合,是线路投资等年值,是线路新建状态二进制变量,是场景下标,是场景集合,是场景年等效小时数,是时间下标,是时间集合,是场景第小时总惩罚费用。
等年值的折算方式如下式所示:
;
其中是总费用的等年值,是一次投资费用,是贴现率,是回收年限,是年损耗费用,是年维护费用。
惩罚项根据其定义进行计算,惩罚项包括节点功率不平衡惩罚、新能源弃用惩罚和节点电压偏移惩罚,其中节点功率不平衡惩罚如下所示:
;
式中是场景第小时功率不平衡惩罚费用,是节点下标,是节点集合,分别为有功功率和无功功率不平衡惩罚单位费用,分别为场景第小时节点的有功功率超负荷、欠负荷松弛因子和无功功率超负荷、欠负荷松弛因子。
新能源弃用惩罚计算方式如下所示:
;
其中是场景第小时新能源弃用惩罚费用,是新能源机组下标,是新能源机组集合,是新能源弃用惩罚单位费用,是场景第小时新能源机组有功功率弃用因子。电压偏移惩罚计算方式如下所示:
;
其中是场景第小时电压偏移惩罚费用,是电压偏移惩罚单位费用,分别是场景第小时节点电压上偏移因子和下偏移因子。综上,目标函数中惩罚项如下:
;
4)规划模型约束条件
本发明所述规划模型的约束条件包括机组出力约束、节点功率平衡约束、线路传输容量约束、线路潮流约束、M3C变频站稳态潮流约束以及惩罚项约束。
1.机组出力约束
火电机组出力约束如下所示:
;
式中为火电机组下标,为火电机组集合,分别为火电机组有功出力上限和下限,分别为火电机组无功出力上限和下限,分别是场景第小时火电机组有功出力和无功出力变量,是第小时机组运行状态二进制变量,
新能源机组出力约束如下:
;
式中是场景第小时新能源机组有功出力变量,是新能源机组装机容量,是场景第小时新能源机组出力系数。
2.节点功率平衡约束
有功功率平衡约束如下:
;
式中代表节点所连的所有火电机组集合,代表节点所连的所有新能源机组集合,是场景第小时线路传输有功功率变量,代表从节点流出功率的所有线路,代表向节点流入功率的所有线路,为联络线下标,表示节点所连联络线,是场景第小时对外联络线传输有功功率,是场景第小时节点变频站有功功率变量,是场景第小时节点有功负荷。且满足在非M3C变频站并网点处为0,即如下式所示:
;
式中表示在节点集合内不含变频站并网点的子节点集合。
无功功率平衡约束如下:
;
式中是场景第小时线路传输无功功率变量,是场景第小时节点变频站无功功率变量,是场景第小时节点无功负荷。且满足在非M3C变频站并网点处为0,即如下式所示:
;
3.线路传输容量约束
线路传输视在功率限制为:
;
式中是线路最大输送视在功率,是线路额定电压,是线路长期运行电流。线路静稳极限约束为:
;
式中是线路静稳极限,是线路电抗值。考虑线路电容充电功率后线路有功传输能力如下:
;
式中是线路有功最大传输能力,是工作频率,是线路单位长度电容,是线路的长度。因此已有线路传输容量约束如下:
;
式中为已建线路集合。待建线路传输容量约束如下:
;
4.线路潮流方程约束
已建线路的潮流方程约束为:
;
其中前者为有功功率方程,后者为无功功率方程。式中分别为线路的电导值和电纳值,分别为节点的电压幅值和相角,其中下标分别表示线路的起始节点和终止节点。候选线路由二进制变量决定是否建设,其潮流方程约束如下所示:
;
;
式中,为引入的大值,用以解耦未建线路的传输功率和端电压。
5.M3C稳态潮流约束
M3C有功功率变量满足功率耦合方程:
;
下标中分别表示变频站的工频并网点和分频并网点,分别为M3C变频器损耗率和换流变压器损耗率。M3C变频站传输的有功功率需要在变频器容量范围内:
;
式中为M3C变频站第种规格容量。M3C变频站提供的无功功率需要满足换流变压器容量的限制:
;
式中为M3C变频站第种规格换流变压器容量。同时每一座变频站只能选择一个容量规格进行建造,因此M3C容量型号二进制变量需要满足以下约束:
;
6.惩罚项约束
新能源弃用惩罚由弃用因子引入,它的定义为当前风速条件下新能源机组理想最大出力与实际出力之差,如下所示:
;
节点电压偏移惩罚同功率不平衡惩罚,由上偏移因子和下偏移因子构成,其与节点电压的关系可以表示为:
;
除上述约束外,所有惩罚因子还都应满足非负这一条件,如下所示:
;
以下提供一个5节点算例对上述模型进行验证。
为验证本发明所述规划模型的可行性,现利用5节点算例进行验证。对该5节点算例进行新增海上风电拓展,针对该场景下海上风电设置分频和工频送出两套方案进行优化。考虑如下表所示的四种典型场景,设置M3C变频器功率损耗率,换流变压器功率损耗率。算例共5个节点,有火电机组5台,总装机为3060MW,有输电线路6条,电压等级为220kV,额定负荷为2000MW。现规划新建两个海上风电场,装机容量分别为800MW和100MW,至入网点敷设距离分别设置为100km和40km,负荷增长为原额定负荷1.3倍,海上风电岸上候选接入点为A1、A5节点,设置回收年限为20年,贴现率设置为0.05,针对以上场景进行电网规划。
表1 四种场景参数
算例的电网示意图如3所示。
本算例主网电压等级为220kV,根据相关市场数据设置主网架空线路为112万元/km。海底电缆根据不同电压等级和不同规格,价格也不相同,本算例海上风电规模接近百万千瓦级别,因此考虑220kV和330kV两种输电等级,敷设费用都为120万元/km。具体投资费用如下表所示:
表2架空线与海底电缆一次投资价格清单
M3C变频器投资费用根据MMC换流器造价进行估算,将M3C变频器、工分频换流变压器等费用相加,得到M3C变频站的一次投资如下:
;
模型计算结果如下:100MW海上风电通过工频220kV电压等级3×500mm2规格的海底电缆1回线送出,接入A1节点;800MW海上风电通过分频330kV等级下3×630mm2规格海底电缆2回线送出,接入A1节点,在A1节点建造800MW变频站,两侧换流变压器容量为1200MVA,分频侧电压等级为330kV,工频侧为220kV。同时工频主网需要进行拓展,新建四条220kV电压等级架空线路,其中三条为由风电接入点A1节点送出到A2、A3、A5三个节点,另一条为由A4节点至A5节点。拓展后电网示意图如4所示,A6和A7节点为800MW和100MW海上风电节点。
规划结果线路及变频站一次投资费用如下表所示:
表3分频组各元件一次投资
表4一次投资费用对比
综上可以得到总一次投资费用为26.623亿元,其中架空线投资1.19504亿元,海底电缆投资16.4760亿元,M3C变频站投资8.952亿元。与之对应对照组一次投资为35.581亿元,高于分频组,由此可以看出海上风电经过分频送出相比工频能够有效降低一次投资。
由算例计算结果可以看出,本发明所提出的电网规划模型能够有效的对海上风电分频并网场景进行网架规划。
以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (9)
1.含M3C变频站选址定容的工分频互联系统电网规划系统,其特征在于,所述工分频互联系统电网规划系统的目标函数,具体如下:
式中,m是候选M3C变频站编号下标,是候选M3C变频站集合,j是M3C变频站规格下标,CSpec是M3C变频站规格集合,Cj是M3C变频站第j种规格的等年值费用,wm,j是变频站m第j种规格新建状态二进制变量,l是线路下标,Lcandi是候选线路集合,Cl是线路l投资等年值,vl是线路l新建状态二进制变量,b是场景下标,S是场景集合,Hb是场景b年等效小时数,t是时间下标,T是时间集合,Pect,b是场景b第t小时总惩罚费用;
所述目标函数用于在工分频互联系统电网规划中,在满足负荷需求的前提下,使线路一次投资、维护和损耗费用最小;
在所述目标函数中,总惩罚费用按下式计算:
其中,
式中,Pipt,b是场景b第t小时功率不平衡惩罚费用,i是节点下标,N是节点集合,Pnsp,Pnsq分别为有功功率和无功功率不平衡惩罚单位费用,分别为场景b第t小时节点i的有功功率超负荷、欠负荷松弛因子和无功功率超负荷、欠负荷松弛因子;
式中,Rdpt,b是场景b第t小时新能源弃用惩罚费用,ri是新能源机组下标,CR是新能源机组集合,Pnsres是新能源弃用惩罚单位费用,是场景b第t小时新能源机组ri有功功率弃用因子;
式中,Vopt,b是场景b第t小时电压偏移惩罚费用,Pnsu是电压偏移惩罚单位费用,分别是场景b第t小时节点i电压上偏移因子和下偏移因子。
2.根据权利要求1所述的含M3C变频站选址定容的工分频互联系统电网规划系统,其特征在于,在所述目标函数中,等年值的折算方式按下式计算:
式中,C是总费用的等年值,Inv是一次投资费用,σ是贴现率,n是回收年限,L是年损耗费用,M是年维护费用。
3.根据权利要求1所述的含M3C变频站选址定容的工分频互联系统电网规划系统,其特征在于,所述工分频互联系统电网规划系统的约束条件,具体包括:
机组出力约束、节点功率平衡约束、线路传输容量约束、线路潮流约束、M3C变频站稳态潮流约束以及惩罚项约束。
4.根据权利要求3所述的含M3C变频站选址定容的工分频互联系统电网规划系统,其特征在于,所述线路潮流约束基于线性化的支路功率方程确定,所述线性化的支路功率方程的线性化过程如下:
获取节点i流向节点j的功率方程,如下:
式中,Pij,Qij是支路传输有功、无功功率,Ui,Uj是节点i和节点j的电压幅值,θi,θj是节点i和节点j的电压相角,rij,xij是支路阻抗参数;
将支路阻抗转换为支路导纳,得到支路功率方程,如下:
其中gij,bij是支路导纳参数;
将节点电压幅值和相角U,θ作为独立变量在U=1,θ=0处对上式进行泰勒展开,得到线性化支路功率方程,如下:
5.根据权利要求4所述的含M3C变频站选址定容的工分频互联系统电网规划系统,其特征在于,所述线路潮流约束,具体如下:
已建线路的潮流方程约束为:
式中,pfl,t,b和qfl,t,b是场景b第t小时线路l传输有功功率变量和无功功率变量,gl,bl分别为线路l的电导值和电纳值,下标fr(l),to(l)分别表示线路l的起始节点和终止节点,Lexist为已建线路集合;
候选线路由二进制变量vl决定是否建设,其潮流方程约束如下所示:
式中,M为引入的大值。
6.根据权利要求3所述的含M3C变频站选址定容的工分频互联系统电网规划系统,其特征在于,所述M3C变频站稳态潮流约束基于M3C变频站稳态简化模型得到,在所述M3C变频站稳态简化模型中,M3C变频站提供的无功支撑受到工/分频侧换流变压器容量的限制,如下式:
式中,Ps M、分别为M3C工频侧网侧有功功率和无功功率,分别为M3C分频侧有功功率和无功功率,ST,s,ST,f分别为工/分频侧换流变压器的容量。
7.根据权利要求6所述的含M3C变频站选址定容的工分频互联系统电网规划系统,其特征在于,所述M3C变频站稳态潮流约束,具体如下:
M3C变频站传输的有功功率约束:
式中,pM为M3C变频站的有功功率,PNj为M3C变频站第j种规格容量,下标中sys(m),ff(m)分别表示变频站m的工频并网点和分频并网点,CM为M3C变频站合集;
M3C变频站提供的无功功率约束:
式中,qM为M3C变频站的无功功率,STj为M3C变频站第j种规格换流变压器容量。
8.根据权利要求7所述的含M3C变频站选址定容的工分频互联系统电网规划系统,其特征在于,M3C容量型号的二进制变量满足如下约束:
9.根据权利要求3所述的含M3C变频站选址定容的工分频互联系统电网规划系统,其特征在于,在所述惩罚项约束中,节点电压偏移惩罚约束,如下:
式中,U为节点电压幅值,下标i为节点编号,和为节点i在场景b第t小时的上偏移因子和下偏移因子,N为节点集合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311659037.0A CN117371154B (zh) | 2023-12-06 | 2023-12-06 | 含m3c变频站选址定容的工分频互联系统电网规划系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311659037.0A CN117371154B (zh) | 2023-12-06 | 2023-12-06 | 含m3c变频站选址定容的工分频互联系统电网规划系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117371154A CN117371154A (zh) | 2024-01-09 |
CN117371154B true CN117371154B (zh) | 2024-05-03 |
Family
ID=89394847
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311659037.0A Active CN117371154B (zh) | 2023-12-06 | 2023-12-06 | 含m3c变频站选址定容的工分频互联系统电网规划系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117371154B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112968441A (zh) * | 2021-03-12 | 2021-06-15 | 内蒙古科技大学 | 一种应用于大规模风电基地的电网规划方法 |
CN114202186A (zh) * | 2021-12-06 | 2022-03-18 | 国网江苏省电力有限公司 | 一种海上风电场并网输送系统经济评估性方法 |
CN117035335A (zh) * | 2023-08-17 | 2023-11-10 | 华北电力大学 | 一种多阶段储能与输电网协同规划方法及系统 |
CN117117870A (zh) * | 2023-08-24 | 2023-11-24 | 华北电力大学(保定) | 一种考虑网络优化重构降损的旋转潮流控制器选址定容方法 |
-
2023
- 2023-12-06 CN CN202311659037.0A patent/CN117371154B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112968441A (zh) * | 2021-03-12 | 2021-06-15 | 内蒙古科技大学 | 一种应用于大规模风电基地的电网规划方法 |
CN114202186A (zh) * | 2021-12-06 | 2022-03-18 | 国网江苏省电力有限公司 | 一种海上风电场并网输送系统经济评估性方法 |
CN117035335A (zh) * | 2023-08-17 | 2023-11-10 | 华北电力大学 | 一种多阶段储能与输电网协同规划方法及系统 |
CN117117870A (zh) * | 2023-08-24 | 2023-11-24 | 华北电力大学(保定) | 一种考虑网络优化重构降损的旋转潮流控制器选址定容方法 |
Non-Patent Citations (2)
Title |
---|
Comprehensive optimization framework for transmission expansion planning with fractional frequency transmission system;Cai, Zhongqi et al.,;《IET GENERATION TRANSMISSION & DISTRIBUTION》;第5355-5365页 * |
考虑配电网韧性的储能系统选址定容优化方法;杨火明 等;《电力建设》(第01期);第30-31页 * |
Also Published As
Publication number | Publication date |
---|---|
CN117371154A (zh) | 2024-01-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108039726B (zh) | 一种基于多智能体系统的能源局域网分布式协同控制方法 | |
CN106503839B (zh) | 一种海上风电场环形集电网络分层规划方法 | |
CN102541621B (zh) | 风光储联合发电系统的仿真等值方法 | |
CN112487659A (zh) | 一种海上风电场集电系统优化设计方法和系统 | |
CN104600713A (zh) | 含风力/光伏发电配电网日前无功调度的生成装置及方法 | |
Buticchi et al. | The role of renewable energy system in reshaping the electrical grid scenario | |
CN104135030A (zh) | 一种用于智能电网的柔性孤岛-并网控制装置及方法 | |
CN105375513A (zh) | 一种基于实时在线等值的110千伏风电场自动电压控制方法 | |
CN111952958A (zh) | 一种考虑控制模式转换的配电网柔性软开关优化配置方法 | |
CN104362650A (zh) | 一种考虑成本因素的电力系统无功优化方法 | |
CN114996908B (zh) | 考虑智能软开关接入的主动配电网扩展规划方法及系统 | |
CN108233408B (zh) | 一种mmc-mtdc系统自适应下垂控制方法 | |
CN108594658A (zh) | 一种电-气耦合系统概率最大负荷裕度多目标优化模型及其求解方法 | |
CN109768543B (zh) | 一种基于混合整数线性规划的弹性保底网架搜索建模方法 | |
Wu et al. | Equivalent modeling method for regional decentralized photovoltaic clusters based on cluster analysis | |
Liao et al. | Model predictive control based coordinated voltage control for offshore radial DC-connected wind farms | |
CN106780119A (zh) | 基于多主动管理的配电网可再生能源消纳方法 | |
CN118117564A (zh) | 一种考虑直流配电形态下配电网多层级光伏消纳模式的协同运行优化方法 | |
CN117371154B (zh) | 含m3c变频站选址定容的工分频互联系统电网规划系统 | |
CN112052543A (zh) | 一种基于混合整数二阶锥规划的保底网架搜索建模方法 | |
Rios et al. | Planning multi-terminal direct current grids based graphs theory. | |
CN109002938B (zh) | 考虑n-1安全准则的交直流混合配电网双层规划方法 | |
CN105811435A (zh) | 一种智能储能发电系统的无功补偿方法 | |
CN116681149A (zh) | 基于全生命周期的海上风电输电方式经济技术分析方法 | |
CN112736913B (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 |