CN112054515A - 一种基于多目标优化的受端电网接纳直流能力检测方法 - Google Patents

一种基于多目标优化的受端电网接纳直流能力检测方法 Download PDF

Info

Publication number
CN112054515A
CN112054515A CN202010882430.6A CN202010882430A CN112054515A CN 112054515 A CN112054515 A CN 112054515A CN 202010882430 A CN202010882430 A CN 202010882430A CN 112054515 A CN112054515 A CN 112054515A
Authority
CN
China
Prior art keywords
direct current
index
receiving
total
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
Application number
CN202010882430.6A
Other languages
English (en)
Other versions
CN112054515B (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.)
State Grid Corp of China SGCC
Wuhan University WHU
China Electric Power Research Institute Co Ltd CEPRI
State Grid Tianjin Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Wuhan University WHU
China Electric Power Research Institute Co Ltd CEPRI
State Grid Tianjin Electric Power Co Ltd
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 State Grid Corp of China SGCC, Wuhan University WHU, China Electric Power Research Institute Co Ltd CEPRI, State Grid Tianjin Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202010882430.6A priority Critical patent/CN112054515B/zh
Publication of CN112054515A publication Critical patent/CN112054515A/zh
Application granted granted Critical
Publication of CN112054515B publication Critical patent/CN112054515B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • 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
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • 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
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393Score-carding, benchmarking or key performance indicator [KPI] analysis
    • 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
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/36Arrangements for transfer of electric power between ac networks via a high-tension dc link
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J5/00Circuit arrangements for transfer of electric power between ac networks and dc networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/60Arrangements for transfer of electric power between AC networks or generators via a high voltage DC link [HVCD]

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Development Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Power Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Game Theory and Decision Science (AREA)
  • Health & Medical Sciences (AREA)
  • Educational Administration (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明涉及一种基于多目标优化的受端电网接纳直流能力检测方法,本发明首先进行直流落点的初步筛选,缩减场景数;接着采用广义等效短路比(GESCR)指标、暂态电压支撑强度(TVSI)指标以及网损指标作为直流落点和机组出力的优化目标,以及考虑多约束条件建立直流馈入受端电网的多目标优化模型;然后通过直流馈入规模的等步长递增并重复进行直流落点和机组优化,直至不满足安全运行约束条件停止;最后对计算分析结果逆序进行N‑1安全稳定校验,最先通过校验结果作为受端电网的最大直流承载规模。本发明提出的计算分析方法逻辑清晰、考虑因素全面、易于工程实现、可有效得出受端电网的最大直流接纳能力。

Description

一种基于多目标优化的受端电网接纳直流能力检测方法
技术领域
本发明涉及一种受端电网接纳直流能力计算方法,尤其是涉及考虑直流落点和机组出力分配对受端电网承载直流能力影响,基于多目标优化方法和构建多约束条件,建立受端电网最大接纳直流能力分析检测方法。
背景技术
高压直流输电在大容量、远距离输电方面具有独特优势,以及在我国能源优化配置中起着关键作用。对于受端电网而言,多直流的同时馈入易引发电压崩溃、频率失稳、换向失败等一系列安全稳定问题,特别是随着馈入直流规模的不断增大,受端电网面临的安全稳定问题日益突出,因此研究建立科学有效的受端电网最大直流接纳能力分析方法,得出受端电网最大直流承载规模,对交直流混联电网的规划调度具有重要参考意义。
目前研究关于受端电网接纳直流能力的评估方法主要为基于多馈入短路比(MSCR)指标,通过不断增大直流馈入规模,直至MSCR指标等于临界值,并进行一些安全稳定校验,能够通过各项安全约束的承载规模即为受端电网的最大直流接纳能力。由此可见,现有研究主要利用MSCR指标量化评估受端电网直流接纳能力,并进行一些安全稳定的辅助验证,但由于多馈入短路比的临界值依赖于工程经验,实际应用中不能较准确的量化受端电网的电压稳定裕度,得出的受端电网直流接纳能力可能有偏差,同时现有研究也没有形成完善的基于多目标优化的受端电网接纳直流能力检测方法体系。此外,现有研究在评估受端电网的直流接纳能力时,没有考虑和涉及发电机组出力的优化分配对直流承载规模的影响。
本发明从直流落点和机组出力优化两个层面,通过循环迭代方式和考虑多约束条件建立评估受端电网最大直流承载规模的指标体系和计算方法。首先进行直流落点的初步筛选,缩减场景数。接着采用广义等效短路比(GESCR)指标、暂态电压支撑强度(TVSI)指标以及网损指标建立直流落点和机组出力为优化目标,以及考虑多约束条件建立直流馈入受端电网的多目标优化模型。通过直流馈入规模的循环递增,重复进行落点和机组出力优化,直至不满足所建立的安全运行约束条件停止,并对计算分析结果逆序进行N-1安全稳定校验,最先通过校验结果作为受端电网的最大直流接纳能力。最后利用IEEE 39节点***作为多直流馈入的受端电网,计算分析得出了受端的最大直流接纳能力,验证了所建立受端电网最大直流接纳能力分析方法的可行性和有效性。
发明内容
本发明的上述技术问题主要是通过下述技术方案得以解决的:
一种基于多目标优化的受端电网接纳直流能力检测方法,其特征在于,所提出方法包括以下步骤:
步骤1、基于确定的受端电网结构和直流馈入的典型运行方式,对受端电网所有可能的直流落点进行初步筛选得到备选落点集D;
步骤2、定义初始受端电网的直流馈入数目Ndc=N0,并设定每条直流输电线路都以最大输送功率Pdmax运行;
步骤3、利用所建立直流落点和机组出力组合优化模型,通过直流落点和机组出力优化确定最优落点和机组的有功、无功出力方案;
步骤4、在优化得出的最优运行方式下,判断是否满足各项安全运行约束条件,若满足条件,进入步骤5,若不满足条件,进入步骤7;
步骤5、第Ndc条直流以步长ΔP增加输送功率,即Pd(Ndc)=Pd(Ndc)+ΔP,同时判断Pd(Ndc)与Pdmax大小,若Pd(Ndc)≤Pdmax,进入步骤3,否则,进入步骤6;
步骤6、执行Ndc=Ndc+1,且设置第Ndc条直流输送功率为最小限定值,即以Pd(Ndc)=Pdmin运行,进入步骤3;
步骤7、对得到的迭代优化结果,按照直流馈入总功率的大小按倒序方式进行N-1安全稳定校验,直至得出在N-1校验中无电压越限、支路潮流越限、潮流不收敛等安全稳定问题的直流馈入方案,以此作为受端最大直流受电规模,进入步骤8;
步骤8、根据步骤7中结果确定最终馈入直流数目Ndc,并统计得出的受端电网最大接纳直流规模PD_total,即
Figure BDA0002654484050000021
在上述方法,为缩减直流落点场景数和提高优化准确率,在步骤1中,定义负荷集中度指标和节点的度数指标对落点进行初步筛选,如式一所示:
Figure BDA0002654484050000031
式一中,F(i)、M(i)分别为节点i的负荷集中度和节点度数指标值;J为节点i的相邻节点集;L为***负荷节点集;PLj、PLk分别为节点j和节点k的有功负荷;lbij为i与j节点之间支路;count为计数函数。
在上述方法,在步骤3中建立直流落点和机组出力组合优化模型时,利用广义等效短路比(GESCR)指标、暂态电压支撑强度指标(ITVSI)以及网损指标构建优化目标,分别为:
①GESCR指标量化评估受端电网的电压支撑能力时,因其具有理论的临界值为1,相较于多馈入短路比(MSCR)指标具有更好准确性和实用性,选用GESCR指标作为评估受端电网电压支撑能力的指标,并分为均衡性指标和安全裕度指标分别为:
Figure BDA0002654484050000032
式二中,Ibal(i)、Imar(i)分别为第i个多直流馈入方案下的均衡性和安全裕度指标值;Ndc为直流馈入总数;GESCR(j)表示第j条直流馈入节点对应的GESCR;GESCRave为第i个方案下各直流馈入对应的GESCR平均值;GESCRmin为节点电压临界稳定对应的GESCR,大小等于1,其中,GESCR指标的定义为:对于LCC型直流,当采用送端定功率、受端定熄弧角控制时,第j条馈入直流对应的GESCR指标为:
Figure BDA0002654484050000033
式三中,ULj为负荷节点j的节点电压;Zeqj为将其它直流馈入对本直流的影响通过阻抗形式折算到换流母线j后,得到的馈入点等值阻抗大小;Sdj为直流馈入节点j的视在功率;QVCF为无功电压补偿系数,表示直流控制特性对GESCR的影响,表达式为:
Figure BDA0002654484050000034
式四中,Qdj为第j条直流换流母线吸收的无功功率;θeqj为等值阻抗Zeqj的阻抗角;
②为反映各换流母线间的相互影响程度,采用暂态电压支撑强度指标(ITVSI)作为优化目标之一,并定义第i种直流馈入方案下的平均暂态电压支撑强度指标为:
Figure BDA0002654484050000041
式五中,Ndc为直流馈入总数;其中,ITVSI指标定义为:
Figure BDA0002654484050000042
式六中,Xekj为只保留换流母线k和j的等值电网中k、j节点之间的转移阻抗;Xek为等值网中节点k与等值电源之间的支路阻抗;
③优化过程中应保证电网运行经济性,定义第i种运行方案下网损指标为:
Iloss(i)=PG_total(i)+PD_total(i)-PL_total(i) 式七
式七中,PG_total为***常规发电机组总的有功出力;PD_total为受端直流馈入总功率;PL_total为***总负荷。
在上述方法,在步骤一和步骤三中,为处理指标数量级和方向上存在的不一致,用极差变换法对第i个方案下的指标进行处理,方法如式八所示,并采用层次分析法和反熵权法对指标进行主客观组合赋权;
Figure BDA0002654484050000043
式八中,I1、I2分别表示负向指标和正向指标;
Figure BDA0002654484050000044
分别为标准化处理后的指标;Imax、Imin表示指标归一化后取值范围,为避免归一化到闭区间[0,1]时指标取值出现为0情况,可取Imax=0.996,Imin=0.002。
在上述方法,在步骤4中构建了***潮流约束、发电机有功和无功约束、节点电压约束、短路电流约束、支路最大传输功率约束、GESCR指标不小于临界值、有效直流惯性时间常数约束、最大单馈入约束等多约束条件,如式九所示:
s.t.
Figure BDA0002654484050000051
式九中,Pi和Qi是节点i注入有功和无功;δij是节点i和j之间的电压相位差;Gij和Bij分别是节点i和j之间的互电导和互导纳;
Figure BDA0002654484050000052
Figure BDA0002654484050000053
分别表示发电机有功和无功出力的上下限;N为***节点总数;NG为***发电机节点总数;Nl为***支路总数;
Figure BDA0002654484050000054
Figure BDA0002654484050000055
分别表示节点电压的上下限;
Figure BDA0002654484050000056
表示节点短路电流控制上限;
Figure BDA0002654484050000057
表示支路传输功率上限;Hdc为有效直流惯性时间常数,取大于2~3;γdcmax,i为因受端***故障使第i条直流闭锁,***在恢复稳态后,为保证***的稳态频率能够维持在允许偏差范围,第下i条直流输送功率占***所有发电机组输出功率的最大比值,其中Hdc和γdcmax,i的计算表达式为:
Figure BDA0002654484050000058
式十中,JA_total为受端***的总转动惯量;PD_total为受端馈入直流总功率;Req为***发电机组等效调差系数;DL为***负荷有功频率调节系数;γ为馈入直流总功率占其与***所有发电机有功出力之和的比值;Δfmax为***最大允许稳态频率偏差;fN为***额定频率。
在上述方法,在步骤3中,利用MOEA/D优化算法对直流落点和机组出力组合优化模型进行求解,为增强MEOA/D优化算法的局部搜索能力,提高优化效果,进化操作后采用高斯变异增加种群多样性;对进化个体中的第i个元素vi变异操作为:
Figure BDA0002654484050000061
式十一中,rand为产生随机数的函数;C为随机数矩阵;Nv表示自变量个数;Gussiant(t)表示自变量为t的高斯函数;δi 2为给定高斯变异方差;Pi为变异率。
因此,本发明建立的基于多目标优化的受端电网接纳直流能力分析方法,相较于目前的技术具有如下优点:基于GESCR指标,***性的建立了受端电网最大直流接纳能力计算分析方法,克服传统方法以MSCR指标评估受端电压支撑能力的模糊性和不准确性;考虑了不同机组出力优化分配对受端承载直流规模的影响,并随着直流馈入规模的迭代增加,考虑最优落点和机组出力分配的变化,进行直流落点和机组出力的重复优化,更具有合理性,相较于先固定最优落点在通过迭代增加直流规模方法得到的受电规模不偏向于保守;另外,考虑了直流馈入后的有效直流惯性时间常数约束、最大单馈入约束等多约束条件,以及对结果按逆序进行安全稳定校核,该方法计算得到的受端最大接纳直流能力更具可靠性和准确性。
附图说明
图1为基于多目标优化的受端电网最大接纳直流能力计算方法流程图。
图2为IEEE 39节点***拓扑图。
具体实施方式
下面通过实施例,并结合附图,对本发明的技术方案作进一步具体的说明。
1.基于多目标优化的受端电网最大接纳直流能力计算方法流程。
本发明基于多目标优化方法提出了受端电网最大接纳直流能力的分析方法,其流程如图1所示,具体方法步骤为:
步骤1、基于确定的受端电网结构和直流馈入的典型运行方式,对受端电网所有可能的直流落点进行初步筛选得到备选落点集D;
步骤2、设定初始受端电网的直流馈入数目Ndc=N0,并设定各条直流线路都以最大输送功率Pdmax运行;
步骤3、利用所建立直流落点和机组出力组合优化模型,通过落点和机组出力优化,确定最优落点及机组的有功、无功出力方案;
步骤4、在得出的最优运行方式下,判断是否满足各项安全运行约束条件,若满足条件,进入步骤5,若不满足条件,进入步骤7;
步骤5、第Ndc条直流以步长ΔP增加输送功率,即Pd(Ndc)=Pd(Ndc)+ΔP,同时判断Pd(Ndc)与Pdmax大小,若Pd(Ndc)≤Pdmax,进入步骤3,否则,进入步骤6;
步骤6、执行Ndc=Ndc+1,且设置第Ndc条直流输送功率为最小限定值,即以Pd(Ndc)=Pdmin运行,进入步骤3;
步骤7、对得到的迭代优化结果,按照直流馈入总功率的大小按倒序方式进行N-1安全稳定校验,直至得出在N-1校验中无电压越限、支路潮流越限、潮流不收敛等安全稳定问题的直流馈入方案,以此作为受端最大直流受电规模,进入步骤8;
步骤8、根据步骤7中结果确定最终馈入直流数目Ndc,并统计得出的受端电网最大接纳直流规模PD_total,即
Figure BDA0002654484050000071
2.直流落点初步筛选。
根据1中所述方法,为缩减场景数和提高优化准确率,在步骤1中,定义负荷集中度指标和节点的度数指标对落点进行初步筛选,如式1所示:
Figure BDA0002654484050000072
式中,F(i)、M(i)分别为节点i的负荷集中度和节点度数指标值;J为节点i的相邻节点集;L为***负荷节点集;PLj、PLk分别为节点j和节点k的有功负荷;lbij为i与j节点之间支路;count为计数函数。
3.直流落点和机组出力优化目标构建。
根据1中所述方法,在步骤3的直流落点和机组出力组合优化模型中,利用广义等效短路比(GESCR)指标、暂态电压支撑强度指标(ITVSI)以及网损指标构建优化目标,分别为:
①GESCR指标量化评估受端电网的电压支撑能力时,因其具有理论的临界值为1,相较于多馈入短路比(MSCR)指标具有更好准确性和实用性,选用GESCR指标作为评估受端电网电压支撑能力的指标,并分为均衡性指标和安全裕度指标分别为:
Figure BDA0002654484050000081
式中,Ibal(i)、Imar(i)分别为第i个多直流馈入方案下的均衡性和安全裕度指标值;Ndc为直流馈入总数;GESCR(j)表示第j条直流馈入节点对应的GESCR;GESCRave为第i个方案下各直流馈入对应的GESCR平均值;GESCRmin为节点电压临界稳定对应的GESCR,大小等于1,其中,GESCR指标的定义为:对于LCC型直流,当采用送端定功率、受端定熄弧角控制时,第j条馈入直流对应的GESCR指标为:
Figure BDA0002654484050000082
式中,ULj为负荷节点j的节点电压;Zeqj为将其它直流馈入对本直流的影响通过阻抗形式折算到换流母线j后,得到的馈入点等值阻抗大小;Sdj为直流馈入节点j的视在功率;QVCF为无功电压补偿系数,表示直流控制特性对GESCR的影响,表达式为:
Figure BDA0002654484050000083
式中,Qdi为第i条直流换流母线吸收的无功功率;θeqi为等值阻抗Zeqj的阻抗角。
②为反映各换流母线间的相互影响程度,采用暂态电压支撑强度指标(ITVSI)作为优化目标之一,并定义第i种直流馈入方案下的平均暂态电压支撑强度指标为:
Figure BDA0002654484050000084
式中,Ndc为直流馈入总数;其中,ITVSI指标定义为:
Figure BDA0002654484050000085
式中,Xekj为只保留换流母线k和j的等值电网中k、j节点之间的转移阻抗;Xek为等值网中节点k与等值电源之间的支路阻抗。
③直流落点和机组出力组合后应保证电网运行经济性,定义第i种运行方案下网损指标为:
Iloss(i)=PG_total(i)+PD_total(i)-PL_total(i) (7)
式中,PG_total为***常规发电机组总的有功出力;PD_total为受端直流馈入总功率;PL_total为***总负荷。
4.指标处理和指标权重确定方法。
根据1中所述方法,在步骤一和步骤三中,为处理指标数量级和方向上的不一致,用极差变换法对第i个方案下的指标进行处理,方法如式(8)所示,并采用层次分析法和反熵权法对指标进行主客观赋权。
Figure BDA0002654484050000091
式中,I1、I2分别表示负向指标和正向指标;
Figure BDA0002654484050000092
分别为标准化处理后的指标;Imax、Imin表示指标归一化后取值范围,为避免归一化到闭区间[0,1]时指标取值出现为0情况,可取Imax=0.996,Imin=0.002。
综上可得,直流落点初步筛选总目标函数和直流落点和机组出力优化总目标函数分别为:
Figure BDA0002654484050000093
式中,fD,i、fi分别为第i个方案下直流落点初步筛选总目标值和落点与机组出力组合优化总目标值;λ1、λ2、α1、α2、α3和α4都为权重因子,其中λ12=1,α1234=1。
5.多约束条件构建。
根据1中所述方法,在步骤4中构建了***潮流约束、发电机有功和无功约束、节点电压约束、短路电流约束、支路最大传输功率约束、GESCR指标不小于临界值、有效直流惯性时间常数约束、最大单馈入约束等多约束条件,如式(10)所示:
s.t.
Figure BDA0002654484050000101
式中,Pi和Qi是节点i注入有功和无功;δij是节点i和j之间的电压相位差;Gij和Bij分别是节点i和j之间的互电导和互导纳;
Figure BDA0002654484050000102
Figure BDA0002654484050000103
分别表示发电机有功和无功出力的上下限;N为***节点总数;NG为***发电机节点总数;Nl为***支路总数;
Figure BDA0002654484050000104
Figure BDA0002654484050000105
分别表示节点电压的上下限;
Figure BDA0002654484050000106
表示节点短路电流控制上限;
Figure BDA0002654484050000107
表示支路传输功率上限;Hdc为有效直流惯性时间常数,一般大于2~3;γdcmax,i为因***故障使第i条直流闭锁时,为保证***在恢复稳态时的稳态频率能够维持在允许偏差范围,第下i条直流输送功率占***所有发电机组输出功率的最大比值,其中Hdc和γdcmax,i的计算表达式为:
Figure BDA0002654484050000108
式中,JA_total为受端***的总转动惯量;Req为***发电机组等效调差系数;DL为***负荷有功频率调节系数;γ为馈入直流总功率占其与***所有发电机有功出力之和的比值;Δfmax为***最大允许稳态频率偏差;fN为***额定频率。
6.优化模型求解。
根据1中所述方法,在步骤3中,利用MOEA/D优化算法对直流落点和机组出力组合优化模型进行求解,为增强MEOA/D优化算法的局部搜索能力,提高优化效果,进化操作后采用高斯变异增加种群多样性。对进化个体中的元素vi变异操作为:
Figure BDA0002654484050000111
式中,rand为产生随机数的函数;C为随机数矩阵;Nv表示自变量个数;Gussiant(t)表示自变量为t的高斯函数;
Figure BDA0002654484050000112
为给定高斯变异方差;Pi为变异率。
7、基于IEEE 39受端网架的实施例。
7.1、参数设置。
采用IEEE 39节点***作为多直流馈入的受端电网,进行最大直流承载规模分析。***拓扑结构如图2所示,***包含29个负荷节点、10个发电机节点,其中节点31为发电机平衡节点,***总负荷为6254.23MW,基准值为100MW。
设馈入的直流都为LCC型直流,且都采用送端定功率、受端定熄弧角控制方式。不失一般性,设馈入直流具有相同的参数,且逆变侧都由12脉动逆变器构成,其中主要直流参数见表1。
表1馈入直流控制参数
Figure BDA0002654484050000113
另外,在对受端电网的最大直流承载规模进行评估和计算过程中,设定***节点电压Ui的允许范围为0.94pu~1.06pu,额定值为1.0pu;节点短路电流控制上限
Figure BDA0002654484050000114
为63kA;支路最大传输功率
Figure BDA0002654484050000115
为1000MW;***总转动惯量15000MW·S;设定***等效调差系数Req为0.036,负荷有功频率调节系数DL为1.5,***额定频率fN为50Hz,***允许最大稳态频率偏差Δfmax为0.5Hz;发电机出力允许范围参见下表3。
7.2、最大接纳能力计算。
为避免换流站与发电厂之间的相互影响,当直流馈入时,暂不考虑直流落入发电机节点,所以29个负荷节点为可能的直流落点,利用负荷集中度指标和节点的度数标准化处理后对所有落点进行筛选排序,选取排名前15个节点作为备选落点集,排序结果及对应的指标值大小如表2所示。
表2直流落点备选集
Figure BDA0002654484050000121
设置直流增加步长为ΔP=0.5pu,初始直流馈入数目设为N0=1,Pdmax、Pdmin分别为1000MW和500MW。MOEA/D算法的最大迭代次数设置为100,高斯变异率为0.1。对于循环迭代结果,从最大馈入直流规模开始按倒序进行N-1安全稳定校验并作适当调整,直至选出能通过安全稳定校验的方案,其最终结果为:直流输送功率Pd1、Pd2和Pd3分别为1000MW、1000MW和650MW,且最优直流落点为节点16、2和9,发电机有功、无功出力允许范围和机端电压如表3所示。
表3机组出力范围及优化结果
Figure BDA0002654484050000131
此直流馈入规模下,在***允许频率偏差Δfmax=0.5Hz时,单直流允许的最大输送功率为1095MW,大于Pdmax=1000MW,所以此时受端电网的最大直流承载规模为3回,共2650MW。
本申请中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。

Claims (6)

1.一种基于多目标优化的受端电网接纳直流能力检测方法,其特征在于,所提出方法包括以下步骤:
步骤1、基于确定的受端电网结构和直流馈入的典型运行方式,对受端电网所有可能的直流落点进行初步筛选得到备选落点集D;
步骤2、定义初始受端电网的直流馈入数目Ndc=N0,并设定每条直流输电线路都以最大输送功率Pdmax运行;
步骤3、利用所建立直流落点和机组出力组合优化模型,通过直流落点和机组出力优化确定最优落点和机组的有功、无功出力方案;
步骤4、在优化得出的最优运行方式下,判断是否满足各项安全运行约束条件,若满足条件,进入步骤5,若不满足条件,进入步骤7;
步骤5、第Ndc条直流以步长ΔP增加输送功率,即Pd(Ndc)=Pd(Ndc)+ΔP,同时判断Pd(Ndc)与Pdmax大小,若Pd(Ndc)≤Pdmax,进入步骤3,否则,进入步骤6;
步骤6、执行Ndc=Ndc+1,且设置第Ndc条直流输送功率为最小限定值,即以Pd(Ndc)=Pdmin运行,进入步骤3;
步骤7、对得到的迭代优化结果,按照直流馈入总功率的大小按倒序方式进行N-1安全稳定校验,直至得出在N-1校验中无电压越限、支路潮流越限、潮流不收敛等安全稳定问题的直流馈入方案,以此作为受端最大直流受电规模,进入步骤8;
步骤8、根据步骤7中结果确定最终馈入直流数目Ndc,并统计得出的受端电网最大接纳直流规模PD_total,即
Figure FDA0002654484040000011
2.根据权利要求1中所述方法,其特征在于,为缩减直流落点场景数和提高优化准确率,在步骤1中,定义负荷集中度指标和节点的度数指标对落点进行初步筛选,如式一所示:
Figure FDA0002654484040000012
式一中,F(i)、M(i)分别为节点i的负荷集中度和节点度数指标值;J为节点i的相邻节点集;L为***负荷节点集;PLj、PLk分别为节点j和节点k的有功负荷;lbij为i与j节点之间支路;count为计数函数。
3.根据权利要求1中所述方法,其特征在于,在步骤3中建立直流落点和机组出力组合优化模型时,利用广义等效短路比(GESCR)指标、暂态电压支撑强度指标(ITVSI)以及网损指标构建优化目标,分别为:
①GESCR指标量化评估受端电网的电压支撑能力时,因其具有理论的临界值为1,相较于多馈入短路比(MSCR)指标具有更好准确性和实用性,选用GESCR指标作为评估受端电网电压支撑能力的指标,并分为均衡性指标和安全裕度指标分别为:
Figure FDA0002654484040000021
式二中,Ibal(i)、Imar(i)分别为第i个多直流馈入方案下的均衡性和安全裕度指标值;Ndc为直流馈入总数;GESCR(j)表示第j条直流馈入节点对应的GESCR;GESCRave为第i个方案下各直流馈入对应的GESCR平均值;GESCRmin为节点电压临界稳定对应的GESCR,大小等于1,其中,GESCR指标的定义为:对于LCC型直流,当采用送端定功率、受端定熄弧角控制时,第j条馈入直流对应的GESCR指标为:
Figure FDA0002654484040000022
式三中,ULj为负荷节点j的节点电压;Zeqj为将其它直流馈入对本直流的影响通过阻抗形式折算到换流母线j后,得到的馈入点等值阻抗大小;Sdj为直流馈入节点j的视在功率;QVCF为无功电压补偿系数,表示直流控制特性对GESCR的影响,表达式为:
Figure FDA0002654484040000023
式四中,Qdj为第j条直流换流母线吸收的无功功率;θeqj为等值阻抗Zeqj的阻抗角;
②为反映各换流母线间的相互影响程度,采用暂态电压支撑强度指标ITVSI作为优化目标之一,并定义第i种直流馈入方案下的平均暂态电压支撑强度指标为:
Figure FDA0002654484040000024
式五中,Ndc为直流馈入总数;其中,ITVSI指标定义为:
Figure FDA0002654484040000031
式六中,Xekj为只保留换流母线k和j的等值电网中k、j节点之间的转移阻抗;Xek为等值网中节点k与等值电源之间的支路阻抗;
③优化过程中应保证电网运行经济性,定义第i种运行方案下网损指标为:
Iloss(i)=PG_total(i)+PD_total(i)-PL_total(i) 式七
式七中,PG_total为***常规发电机组总的有功出力;PD_total为受端直流馈入总功率;PL_total为***总负荷。
4.根据权利要求1中所述方法,其特征在于,在步骤一和步骤三中,为处理指标数量级和方向上存在的不一致,用极差变换法对第i个方案下的指标进行处理,方法如式八所示,并采用层次分析法和反熵权法对指标进行主客观组合赋权;
Figure FDA0002654484040000032
式八中,I1、I2分别表示负向指标和正向指标;
Figure FDA0002654484040000033
分别为标准化处理后的指标;Imax、Imin表示指标归一化后取值范围。
5.根据权利要求1中所述方法,其特征在于,在步骤4中构建了***潮流约束、发电机有功和无功约束、节点电压约束、短路电流约束、支路最大传输功率约束、GESCR指标不小于临界值、有效直流惯性时间常数约束、最大单馈入约束等多约束条件,如式九所示:
Figure FDA0002654484040000041
式九中,Pi和Qi是节点i注入有功和无功;δij是节点i和j之间的电压相位差;Gij和Bij分别是节点i和j之间的互电导和互导纳;
Figure FDA0002654484040000042
Figure FDA0002654484040000043
分别表示发电机有功和无功出力的上下限;N为***节点总数;NG为***发电机节点总数;Nl为***支路总数;
Figure FDA0002654484040000044
Figure FDA0002654484040000045
分别表示节点电压的上下限;
Figure FDA0002654484040000046
表示节点短路电流控制上限;
Figure FDA0002654484040000047
表示支路传输功率上限;Hdc为有效直流惯性时间常数,取大于2~3;γdcmax,i为因受端***故障使第i条直流闭锁,***在恢复稳态后,为保证***的稳态频率能够维持在允许偏差范围,第下i条直流输送功率占***所有发电机组输出功率的最大比值,其中Hdc和γdcmax,i的计算表达式为:
Figure FDA0002654484040000048
式十中,JA_total为受端***的总转动惯量;PD_total为受端馈入直流总功率;Req为***发电机组等效调差系数;DL为***负荷有功频率调节系数;γ为馈入直流总功率占其与***所有发电机有功出力之和的比值;Δfmax为***最大允许稳态频率偏差;fN为***额定频率。
6.根据权利要求1中所述方法,其特征在于,在步骤3中,利用MOEA/D优化算法对直流落点和机组出力组合优化模型进行求解,为增强MEOA/D优化算法的局部搜索能力,提高优化效果,进化操作后采用高斯变异增加种群多样性;对进化个体中的第i个元素vi变异操作为:
Figure FDA0002654484040000051
式十一中,rand为产生随机数的函数;C为随机数矩阵;Nv表示自变量个数;Gussiant(t)表示自变量为t的高斯函数;
Figure FDA0002654484040000052
为给定高斯变异方差;Pi为变异率。
CN202010882430.6A 2020-08-28 2020-08-28 一种基于多目标优化的受端电网接纳直流能力评估方法 Active CN112054515B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010882430.6A CN112054515B (zh) 2020-08-28 2020-08-28 一种基于多目标优化的受端电网接纳直流能力评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010882430.6A CN112054515B (zh) 2020-08-28 2020-08-28 一种基于多目标优化的受端电网接纳直流能力评估方法

Publications (2)

Publication Number Publication Date
CN112054515A true CN112054515A (zh) 2020-12-08
CN112054515B CN112054515B (zh) 2022-12-06

Family

ID=73606473

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010882430.6A Active CN112054515B (zh) 2020-08-28 2020-08-28 一种基于多目标优化的受端电网接纳直流能力评估方法

Country Status (1)

Country Link
CN (1) CN112054515B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102723726A (zh) * 2012-06-01 2012-10-10 中国电力科学研究院 一种确定交流受端***最大直流受入量的方法
CN105048521A (zh) * 2015-04-08 2015-11-11 国电南瑞科技股份有限公司 一种考虑多约束的受端电网直流接入能力计算方法
CN105281325A (zh) * 2015-10-10 2016-01-27 国网天津市电力公司 一种提升受端电网外受电能力的网架优化方法
CN105356481A (zh) * 2015-11-18 2016-02-24 中国电力科学研究院 一种基于多馈入短路比的动态无功补偿选点方法
CN106845752A (zh) * 2016-11-08 2017-06-13 国家电网公司 一种大规模特高压互联电网受电规模评价体系
CN107134799A (zh) * 2017-04-19 2017-09-05 国家电网公司 一种基于短路比的特高压电网换流器无功配置方法
JP2017158265A (ja) * 2016-02-29 2017-09-07 パナソニックIpマネジメント株式会社 電力供給システム、及び電力変換システム
CN107681682A (zh) * 2017-10-25 2018-02-09 国家电网公司 一种基于ward等值的交直流***等值方法
CN109873442A (zh) * 2019-04-03 2019-06-11 广东电网有限责任公司电力调度控制中心 一种多馈入直流***新增直流落点方案综合评估方法
CN110048446A (zh) * 2019-03-21 2019-07-23 中国电力科学研究院有限公司 一种确定分层直流接入***受端最佳落点的方法及***

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102723726A (zh) * 2012-06-01 2012-10-10 中国电力科学研究院 一种确定交流受端***最大直流受入量的方法
CN105048521A (zh) * 2015-04-08 2015-11-11 国电南瑞科技股份有限公司 一种考虑多约束的受端电网直流接入能力计算方法
CN105281325A (zh) * 2015-10-10 2016-01-27 国网天津市电力公司 一种提升受端电网外受电能力的网架优化方法
CN105356481A (zh) * 2015-11-18 2016-02-24 中国电力科学研究院 一种基于多馈入短路比的动态无功补偿选点方法
JP2017158265A (ja) * 2016-02-29 2017-09-07 パナソニックIpマネジメント株式会社 電力供給システム、及び電力変換システム
CN106845752A (zh) * 2016-11-08 2017-06-13 国家电网公司 一种大规模特高压互联电网受电规模评价体系
CN107134799A (zh) * 2017-04-19 2017-09-05 国家电网公司 一种基于短路比的特高压电网换流器无功配置方法
CN107681682A (zh) * 2017-10-25 2018-02-09 国家电网公司 一种基于ward等值的交直流***等值方法
CN110048446A (zh) * 2019-03-21 2019-07-23 中国电力科学研究院有限公司 一种确定分层直流接入***受端最佳落点的方法及***
CN109873442A (zh) * 2019-04-03 2019-06-11 广东电网有限责任公司电力调度控制中心 一种多馈入直流***新增直流落点方案综合评估方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
BINJIE CHENG等: "Optimal DC-Segmentation for Multi-Infeed HVDC Systems Based on Stability Performance", 《IEEE TRANSACTIONS ON POWER SYSTEMS》 *
HAMMAM ABDELWAHAB等: "Maximizing Solar Energy Input for Cubesat Using Sun Tracking System and A Maximum Power Point Tracking", 《2017 INTERNATIONAL CONFERENCE ON COMMUNICATION, CONTROL, COMPUTING AND ELECTRONICS ENGINEERING》 *
汤奕等: "特高压直流分层接入方式下受端交流***接纳能力分析", 《中国电机工程学报》 *

Also Published As

Publication number Publication date
CN112054515B (zh) 2022-12-06

Similar Documents

Publication Publication Date Title
CN108418255B (zh) 一种适用于含高渗透率新能源的特高压直流外送电网规划方法及***
CN113595158B (zh) 配售电竞争态势下区域配电网的供电能力评估方法
CN107679658B (zh) 一种高比例清洁能源接入下的输电网规划方法
CN110896231B (zh) 一种扶贫区配电网接纳分布式光伏能力计算方法及***
CN111446721B (zh) 一种基于暂态电压灵敏度的配电网调压控制方法
CN113241803B (zh) 一种基于新能源消纳的储能调度方法及计算机介质
CN110970905B (zh) 一种优化风电场电压控制能力的多电源无功配合方法及***
CN111614110B (zh) 一种基于改进多目标粒子群算法的受端电网储能优化配置方法
CN113962159A (zh) 一种基于合理弃光的配电网光伏最大接入容量的评估方法
CN109390971B (zh) 一种基于门当户对遗传算法的配电网多目标主动重构方法
CN115021308A (zh) 一种考虑负荷脱网的配电网中分布式光伏承载能力计算方法
CN104732302A (zh) 一种基于免疫算法的多级电磁环网优化解环方法
CN108512219B (zh) 一种电压稳定约束下的多通道受端***受电能力评估方法
CN112054515B (zh) 一种基于多目标优化的受端电网接纳直流能力评估方法
CN107732972A (zh) 一种地区电网风电接纳能力评估方法
CN114996908B (zh) 考虑智能软开关接入的主动配电网扩展规划方法及***
CN113363991B (zh) 一种基于综合灵敏度的潮流控制方法
CN111697607A (zh) 一种多端柔性直流输电受端电网接入方法及***
Elende et al. Distribution network optimization by optimal sizing and placement of d-statcom using teaching and learning based optimization algorithm
CN115051392A (zh) 一种含新能源、储能及柔性负荷的无功协调优化控制方法
Reddy et al. Optimal placement of interline power flow controller (IPFC) to enhance voltage stability
CN113241793A (zh) 一种计及风电场景的含ipfc电力***预防控制方法
CN111884220A (zh) 一种考虑区域源-荷均衡的电磁环网开环分区方案生成方法
CN112446521B (zh) 考虑经济性和安全性的风电场接入***多目标规划方法
Yang et al. A fuzzy inference based method for evaluating coordinated generation-grid-load-storage control ability in receiving-end power system

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