CN105610192A - 考虑大规模风电接入的在线风险评估方法 - Google Patents
考虑大规模风电接入的在线风险评估方法 Download PDFInfo
- Publication number
- CN105610192A CN105610192A CN201610052437.9A CN201610052437A CN105610192A CN 105610192 A CN105610192 A CN 105610192A CN 201610052437 A CN201610052437 A CN 201610052437A CN 105610192 A CN105610192 A CN 105610192A
- Authority
- CN
- China
- Prior art keywords
- wind
- electricity generation
- powered electricity
- wind power
- risk
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000012502 risk assessment Methods 0.000 title claims abstract description 27
- 230000010354 integration Effects 0.000 title claims abstract description 5
- 230000003068 static effect Effects 0.000 claims abstract description 16
- 230000005611 electricity Effects 0.000 claims description 85
- 239000011159 matrix material Substances 0.000 claims description 24
- 238000012545 processing Methods 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 7
- 230000009466 transformation Effects 0.000 claims description 4
- 230000001186 cumulative effect Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 14
- 230000008901 benefit Effects 0.000 abstract description 5
- 238000004422 calculation algorithm Methods 0.000 abstract description 5
- 238000003672 processing method Methods 0.000 abstract description 2
- 238000007519 figuring Methods 0.000 abstract 2
- 230000035699 permeability Effects 0.000 description 20
- 238000011156 evaluation Methods 0.000 description 6
- 238000004458 analytical method 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
- 238000011161 development Methods 0.000 description 3
- 230000004069 differentiation Effects 0.000 description 3
- 230000002265 prevention Effects 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000005431 greenhouse gas Substances 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 102000030592 phosphoserine aminotransferase Human genes 0.000 description 1
- 108010088694 phosphoserine aminotransferase Proteins 0.000 description 1
- 230000003449 preventive effect Effects 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
-
- 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/0635—Risk analysis of enterprise or organisation activities
-
- 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
- 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]
-
- 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
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Educational Administration (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Power Engineering (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Wind Motors (AREA)
Abstract
本发明公开了考虑大规模风电接入的在线风险评估方法,包括,利用基于风电波动的马尔科夫链模型预测未来时刻th+1时风电出力区间及相应的概率;对于每个预测的风电出力区间,利用局部线性化处理的方法来预估电网运行状态;利用所求结果以及能够反应静态电压稳定裕度的严重度函数求得相应风电出力区间的严重度;利用得到的结果求得一个计算周期内的运行风险指标,根据求得的运行风险指标进行在线风险评估。提出了新的风险指标在线评估短期风电波动对静态电压稳定的影响,所提风险指标融合了短期风电预测精度高和戴维南等值算法算量小的优点。
Description
技术领域
本发明涉及电网安全评估技术领域,尤其涉及一种考虑大规模风电接入的在线风险评估方法。
背景技术
风能作为新型的清洁能源,成为世界上很多国家能源战略的发展重心。中国风电装机容量在2014年达到114GW,排名世界第一位,并且政府计划直至2020年,风电每年的装机容量至少增加20GW。风电的大力发展可以减少温室气体的排放,节约化石能源的消耗,但由于风电出力的不确定性和间歇性,其大规模并网给电网安全运行带了了巨大的挑战。
风险评估方法可考虑电力系统的各种概率不确定因素,是近期用来研究电力系统静态电压稳定的热点。文献“RiskAssessmentofGenerationandTransmissionSystemsConsideringWindPowerPenetration”将风电出力分布及各元件故障作为状态量,运用蒙特卡洛模拟方法评估系统运行风险,此方法精度与模拟次数成正比,计算量较大,主要用于日前运行计划的确定。由于风电出力分钟级的波动已较为明显,为应对风电出力短期波动带来的影响,亟需提出在线安全风险评估方法。文献“RiskAssessmentModelforWindGeneratorTrippingOff”针对风场脱网场景,提出了在线评估系统安全风险的分析方法,风场在未脱网时的随机出力也会给系统安全运行带来风险,因此并网风场对系统安全运行的影响也需关注。文献“Onlinerisk-basedsecurityassessment”提出了在线风险评估的基础思路,采用预估的运行状态及相应严重度确定系统短期运行风险,其中运用了连续潮流方法获得负荷裕度,然而在针对节点数较多的大电力系统时,连续潮流计算量较大,对在线计算造成困难。
发明内容
本发明的目的就是为了解决上述问题,提供一种考虑大规模风电接入的在线风险评估方法,提出了新的风险指标在线评估短期风电波动对静态电压稳定的影响,所提风险指标融合了短期风电预测精度高和戴维南等值算法算量小的优点,且在风险指标中给出了新的严重度函数,对安全域内不同运行状态进行差异化处理,提高了评估的准确性。
为了实现上述目的,本发明采用如下技术方案:
考虑大规模风电接入的在线风险评估方法,包括以下步骤:
步骤一,当前时刻为th,利用基于风电波动的马尔科夫链模型预测未来时刻th+1时风电出力区间及相应的概率;
步骤二,对于每个预测的风电出力区间,利用局部线性化处理的方法来预估电网运行状态;
步骤三,利用所求结果以及能够反应静态电压稳定裕度的严重度函数求得相应风电出力区间的严重度;
步骤四,利用步骤一和步骤三得到的结果求得一个计算周期(th时刻到th+1时刻)内的运行风险指标,根据求得的运行风险指标进行在线风险评估。
所述步骤一中,风电出力变化为随机过程,形成一个状态空间数目为K的状态空间I,I中任何一个状态都有上下限;利用马尔科夫链将来的状态独立于过去的状态,且只依赖于现在的状态得到元素为Pi,j的状态转移矩阵P;每个元素Pi,j都能用转移频数矩阵N中的元素Ni,j表示出;预测的风场在th+1时刻的出力都为一个K维向量中的任一元素为对应的区间的中值,而对应的区间的上下限都能由当前状态及状态转移矩阵计算得到,由此得到th+1时刻风电出力的区间分布。
风电出力为对应的概率为Pj,i,将其作为风险评估中的不确定场景概率Pr(Xi),其中Xi为th+1时刻风电出力的第i个不确定区间。
所述步骤二的具体方法为,风场按恒功率因数控制,将系统的有功和无功的变化量用常规潮流计算迭代公式表示,对该公式雅可比矩阵求逆,得到节点电压向量的变化矩阵,由节点电压向量的变化矩阵得到预估的节点电压幅值和相角向量。
所述步骤四的具体方法为,根据戴维南等值理论得到负荷节点的戴维南等值参数,应用PV曲线中电压在鼻尖点处只有一个取值的原理计算负荷节点的负荷裕度λcr,将风电并网节点处理成恒功率因数控制的PQ节点;假设电力系统共包含N个PQ节点,将第m个PQ节点作为观察节点,计算出其负荷裕度为λcr,m,m=1,2,3…N,计算所有PQ节点的负荷裕度,选择其中负荷裕度最小值作为系统的负荷裕度λmin;结合th+1时刻第i个区间的风电出力,得到对应风电出力第i个区间时系统的最小负荷裕度λmin,i,将λmin,i平方后取倒数求得th+1时刻第i风电出力预测区间的严重度。
风险指标的求解方法为,假设在th+1时刻风电出力共有K个不确定区间,每个不确定区间的概率与该区间对应的严重度相乘后累加得到th时刻到th+1时刻周期内的风险指标。
本发明的有益效果:
提出了新的风险指标在线评估短期风电波动对静态电压稳定的影响,所提风险指标融合了短期风电预测精度高和戴维南等值算法算量小的优点,且在风险指标中给出了新的严重度函数,对安全域内不同运行状态进行差异化处理,提高了评估的准确性。
附图说明
图1为简单两节点系统;
图2为多端口网络模型;
图3为在线风险评估系统流程图;
图4为IEEE9节点系统单线图;
图5为30%渗透率下风险指标;
图6为50%渗透率下风险指标;
图7为IEEE39节点系统单线图;
图8为25%渗透率下风险指标;
图9为44.1%渗透率下风险指标。
具体实施方式
下面结合附图与实施例对本发明作进一步说明。
本发明提出新的风险指标在线评估短期风电波动对静态电压稳定的影响,所提风险指标融合了短期风电预测精度高和戴维南等值算法算量小的优点,且在风险指标中给出了新的严重度函数对安全域内不同运行状态进行差异化处理。最后采用IEEE9节点系统和IEEE39节点系统仿真验证了所提风险指标的有效性。
1大规模风电并网的风险评估
1.1风险指标定义
风险指标是衡量事件发生概率及后果的综合指数,用于评估对目标的不确定影响。在电力系统中,学者将风险的概念引入,用于评估元件在一定的强迫停运率下系统的运行风险,考虑不确定性的风险评估方法可以提供更准确的信息,发挥电网的运行潜力。随着风电渗透率的提高,风电出力的波动对电网静态安全的影响也日渐显著,甚至会引起电压失稳,因此需要将电网运行风险的概念延拓到风电出力波动对电网所造成的影响,而不仅仅研究风机脱网或电网元件出现不确定故障的情形。因此,本发明将风险指标用于评估风电出力波动对电网静态安全稳定的影响,风险指标的定义如下:
式中:K为不确定场景的总数目;Xi为第i个不确定场景,i∈K;Pr(Xi)为第i个场景发生的概率;Sev(Xi)为第i个场景发生时对应的严重度。式(1)是风险指标的通用公式,应用到电力系统中,场景概率和严重度可根据实际需要具备不同的含义。比如场景概率可以为某元件的故障率,也可以为某事故的发生率;严重度可以是低压电程度、频率偏移量、支路潮流越限量,也可以是负荷裕度。本发明将关注点放在风电出力波动对系统静态电压稳定的影响上,具体分析上述变量应用到此场景时的定义。
假定当前运行时刻为th,由于电力系统调度指令下达的周期为15min,15min后的时刻为th+1。将15min作为一个计算周期,运用短期风电出力预测方法分析在未来15min内风电出力变化的所有可能状态,获得的状态数目为K。其中第i个可能状态表示为Xi,i∈K,其对应的概率为Pr(Xi),严重度为Sev(Xi)。将各区间个概率与严重度相乘可得到每个风电出力状态时的风险值,累加可以获得电力系统的综合风险指标。调度人员不仅可以参考综合风险指标来衡量系统短期运行风险,也可以单独分析其中小概率却大严重度的场景,从源头处寻找解决方案,进一步降低此类小概率事件的发生甚至杜绝其发生。
1.2Pr(Xi)的计算
由于风电出力的波动性和间歇性,风电大规模并网造成电网状态的频繁波动,因此需要运用一些方法对风电出力进行预测,以便及时做好预防措施应对风电出力随机波动带来的可能后果。
本发明风险评估的时间框架为未来15分钟,因此采用短期风电出力预测方法。马尔可夫链模型在短期风电预测中具有良好的效果,其基于离散风电出力统计建立模型进行预测,由于本发明主要考虑的是风电波动所带来的影响,需要囊括风电出力在一个计算周期内的所有可能状态,基于风电波动的马尔科夫链模型统计了各时段风电出力的所有变化,然后细分状态区间,能够得到反映风电出力在一个时段波动的区间分布情况,将应用此模型计算Pr(Xi)的过程简述如下。
假定两个相邻的计算时刻为th-1和th,时间间隔为Δt,则有th=th-1+Δt,两个时刻对应的风电出力分别为和两个时刻风电出力变化为有
令为风电出力变化的随机过程,其状态空间为状态空间数目为K,其中任一状态的上下限为由于马尔科夫链将来的状态独立于过去的状态,且只依赖于现在的状态,因此可表达为下式:
式中P为状态转移矩阵,它的元素为Pi,j,要满足下式:
令矩阵N为转移频数矩阵,矩阵中的元素Ni,j表示从状态转移到状态的数目,由于状态空间数目为K,则矩阵P和N都为K×K维矩阵,矩阵P的元素Pi,j可由下式计算:
令为预测的风场在th+1时刻的出力,为一个K维向量,其中任一元素为对应的区间的中值,区间上下限可由当前状态及状态转移矩阵计算,假设th+1和th两个时刻风电出力的变化为公式如下:
得到th+1时刻风电出力的区间分布,取各区间中值作为该区间的风电出力计算,对风电出力若th和th-1两个时刻风电出力变化对应的状态为则风电出力为对应的概率为Pj,i,将其作为风险评估中的不确定场景概率,即:
Pr(Xi)=Pj,ii,j=1,2,…,K(6)
式中,Xi为风险评估中第i个不确定场景,在本发明的风险评估中,其详细含义为th+1时刻风电出力的第i个不确定区间,其概率为Pr(Xi),此概率由基于风电波动的马尔可夫链模型统计得到,包含了计算周期内风电波动转移的所有情况,从而可得到较为全面的系统短期运行风险信息。
1.3Sev(Xi)的计算
文献“Onlinerisk-basedsecurityassessment”构造了低电压严重度的函数,当电压满足上下限约束时,严重度为0,超过上下限时,严重度与越限程度呈线性关系。如果仅设定电压的上下界约束并不能真实反映当前运行状态距离电压失稳边界的距离,且在安全域内,不同的运行状态具备不同的电压稳定裕度,将严重度都定义为0无法反应出不同运行状态时电压稳定裕度的差异,因此本发明提出能够反应静态电压稳定裕度严重度函数,更具有实际意义。连续潮流(CPF)是一种比较成熟的计算静态电压稳定裕度的方法,但是CPF算量巨大,难以用于在线计算。随着PMU的推广,基于广域量测的戴维南等值方法可以计算静态电压稳定裕度,且计算量较连续潮流法小,因此本发明将戴维南等值方法扩展到含风电电力系统中,得到戴维南等值参数后再量化计算系统静态电压稳定的严重度。
在戴维南等值理论中,任意复杂的电力系统从某一负荷节点侧观察,都可以等值为图1中两节点系统。
在图1中,S∠Φ=P1+jQ1为观察节点的视在功率和功率因数,其中P1为有功负荷,Q1为无功负荷;Z∠α=R+jX为从观察负荷处的戴维南等值阻抗,R为戴维南等值电阻,X为戴维南等值电抗;E∠δ1和U∠δ2为相应的戴维南等值电势和负荷节点电压赋值和相角。
利用耦合单端口的概念来追踪每一个负荷节点的戴维南等值参数,多端口网络拓扑见图2。
图2中,VG1,VG2,…,VGm代表电力系统中m个发电厂,SL1,SL2,…SLn代表n个负荷,ZLL为发电厂和负荷之间的阻抗关联矩阵,它们之间的关系可用下式表示:
式中,Y矩阵代表电网导纳矩阵,V和I分别代表电压和电流向量,下标L,T和G分别表示负荷节点,联结节点和发电节点,联结节点即既没有发电厂也没有负荷的节点。式(7)经过变形,可推导得到下式:
VL=KVG-ZLLIL
其他负荷节点对观察节点的影响部分作为耦合因素,此耦合因素可以作为戴维南等值电势或戴维南等值阻抗的一部分,经验证作为耦合阻抗所得结果比较准确,因此对任一负荷节点j,有下式成立:
式中,Zcj为耦合阻抗,Ecoupled-j为耦合电势,Zeq为阻抗关联矩阵的对角元素,Zth和Eth即观察节点的戴维南等值阻抗和等值电势。
得到负荷节点的戴维南等值参数,应用PV曲线中电压在鼻尖点处只有一个取值的原理可以计算该负荷节点的负荷裕度λcr,公式如下:
在含风场电力系统中,通常可将风电并网节点处理成恒功率因数控制的PQ节点,上述公式同样适用于风场并网节点。假设电力系统共包含N个PQ节点,将第m个PQ节点作为观察节点,可以计算出其负荷裕度为λcr,m(m=1,2,3…N),计算所有PQ节点的负荷裕度,选择其中负荷裕度最小值作为系统的负荷裕度λmin。结合th+1时刻第i个区间的风电出力,可得到对应风电出力第i个区间时系统的最小负荷裕度λmin,i,有
λmin,i=min(λcr,m)m=1,2,...,Ni=1,2,...,K(10)
负荷裕度越小,当前运行状态距离静态电压失稳边界越近,后果越严重,严重度为式(11):
对于th+1时刻风电出力预测的区间i,其严重度与最小负荷裕度成反比,不仅能够说明负荷裕度越小,严重度越大,同时也能说明运行状态距离边界越近,严重度呈指数增长,相比线性函数更能引起调度人员的注意。一般运行时期望负荷裕度在10%以上,对应严重度应在10以下,因此严重度若在10附近时需密切关注。
1.4基于灵敏度的电网状态预估
利用1.2部分获得的th+1时刻风电出力分布区间,可以计算各区间对应的严重度,电网在th+1时刻对应各风电出力区间的潮流状态是严重度计算的前提,一般通过常规潮流计算方法获得。但对复杂大电力系统,用潮流迭代方法获得系统运行状态计算量太大,不适用于在线计算,因此本节给出了局部线性化处理方法来预估电网运行状态。
常规潮流计算迭代公式可简写如下:
式中,[J]代表th时刻电网雅可比矩阵,其本质为电网功率变化对节点电压状态变化的灵敏度矩阵;为系统有功和无功的变化向量;为系统节点电压幅值和相角的变化向量。
对th+1时刻预测第i个风电出力区间的风电出力波动为风场按恒功率因数控制,无功出力变化设为令将式(12)雅可比矩阵求逆,可得到节点电压向量的变化矩阵,即
针对预测的风电出力第i个区间,可以得到预估的节点电压向量为
式中,为th时刻系统节点电压幅值和相角向量,为对应预测的风电出力变化第i个区间的th+1时刻系统节点电压幅值和相角向量。按此方法预估系统运行状态可以显著减少风险评估过程的计算量。
2风险评估流程
本发明旨在构建一个能够在线评估由风电波动造成静态电压稳定风险的系统,提供风险指标供调度运行人员参考,可以作为辅助决策的一部分。整个风险评估流程图见图3。
系统当前运行时刻为th,需要评估15分钟后即th+1时刻的运行风险,其流程如下:
1.应用1.2部分介绍的马尔可夫链模型预估th+1时刻风电出力的区间及相应概率;
2.对每个预测的风电出力区间,用1.4部分介绍方法预估对应的系统运行状态;
3.基于上述结果,用1.3部分介绍算法计算相应系统严重度;
4.用(1)式计算获得系统在下一个计算周期内的运行风险指标;
3仿真算例
3.1IEEE9节点系统仿真
本发明采用PSAT工具进行仿真分析。首先采用IEEE9节点系统验证所提风险指标计算的有效性,IEEE9节点系统单线图如图4所示,将母线2由PV节点更改为风场并入系统,风场的历史运行数据由位于山东省济南市的一个风场提供,时间跨度为2012年1月1日至2014年12月31日,采样间隔为15分钟。
将该系统整体负荷水平提升到初始值的1.3倍,风电的渗透率设为30%,风场按PQ节点处理,采用恒功率因数0.96控制。将风场的历史数据按照1.2介绍方法转化为马尔可夫链模型,经过统计计算,可以获得时间间隔为15分钟的各时刻风电出力预测值、概率、严重度和风险指标。取2015年1月1日前16个时刻进行分析,执行前述风险指标计算流程,可获得30%渗透率下2015年1月1日前4个小时的风险指标见图5。
图中蓝色柱状体代表风险指标,由图5,可看出30%渗透率下,这16个时刻风险指标的波动较为平缓,以2015年1月1日00:30分为例,结果见表1。
表1风电30%渗透率下出力预测区间及风险
从表1可看出,各区间的严重度都较小,且该渗透率下系统运行状态具备较高的负荷裕度,风险指标良好,无需进行预防控制。而随着风电的发展,风电渗透率有所提高,风电渗透率达到50%时,可得到不同于图5结果的风险指标,见图6。
图6的风险指标均方差为21.8549,而图5均方差为0.0451,因此高渗透率下风电随机波动将导致风险指标的变化更剧烈,其中第二个时刻的风险指标最高,其结果见表2。
表2风电50%渗透率下出力预测区间及风险
表2中,风电出力各区间严重度数值相比表1高,风险指标也较表2高,表2数据表明50%渗透率下系统运行状态距离静态电压稳定边界比较接近,系统运行裕度较小,再受到扰动有可能导致发生电压崩溃现象,此时的运行状况比较恶劣,需采取预防控制措施进行调控。
3.2IEEE39节点系统仿真
前面使用IEEE9节点系统验证了所提风险指标和预防控制模型,但由于算例简单,计算量小,无法突显出所使用的基于灵敏度的方法在计算速度上的优势,本节使用IEEE39节点系统再次进行验证。IEEE39节点系统单线图见图7,将原PQ节点12、18、26更改为风场,这三个风场运行数据仍用前述的山东风场数据,即三个风场运行状态的变化一样,仅出力大小比例不同。三个风场出力所占比例分别为18%、42%、40%,将系统负荷水平设置为初始值的1.1倍,带来的不平衡出力由所有发电机按当前有功出力的比例共同承担。设置三个风场出力总渗透率分别为25%和44.1%,计算2015年月1日前4个小时共16个时刻的风险指标,比较高渗透率和低渗透率下的结果,见图8和图9。
图8为低渗透率25%时系统4个小时的风险指标图,图9为高渗透率44.1%时系统4个小时的风险指标图。图8中16个风险指标的均方差为0.0229,图9中16个风险指标的均方差为11.7095,因此高渗透率下风电出力随机波动范围变大,所引起的风险指标波动也更为剧烈。通过分析图9各风险指标,第6个时刻的风险值为3.3517,该风险指标及以下时系统运行状态较好,有较高负荷裕度,选择此风险作为风险指标门槛值,高于该风险指标的场景需采取预防控制措施。
4结语
风电作为清洁能源代表得到了重大关注,其发展十分迅速,风电在电网中的渗透率也越来越高。风电的波动和间歇性导致电网运行状态频繁波动,对电网运行提出了新的挑战。本发明延拓风险传统定义,提出新的风险指标在线评估由风电波动引起静态电压风险,分别用IEEE9节点系统和IEEE39节点系统验证了所提风险指标的有效性。本发明所提在线风险评估系统可同目前存在日前调度系统互为补充,为电网安全运行做出贡献。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (6)
1.考虑大规模风电接入的在线风险评估方法,其特征是,包括以下步骤:
步骤一,当前时刻为th,利用基于风电波动的马尔科夫链模型预测未来时刻th+1时风电出力区间及相应的概率;
步骤二,对于每个预测的风电出力区间,利用局部线性化处理的方法来预估电网运行状态;
步骤三,利用所求结果以及能够反应静态电压稳定裕度的严重度函数求得相应风电出力区间的严重度;
步骤四,利用步骤一和步骤三得到的结果求得一个计算周期内的运行风险指标,根据求得的运行风险指标进行在线风险评估。
2.如权利要求1所述考虑大规模风电接入的在线风险评估方法,其特征是,所述步骤一中,风电出力变化为随机过程,形成一个状态空间数目为K的状态空间I,I中任何一个状态都有上下限;利用马尔科夫链将来的状态独立于过去的状态,且只依赖于现在的状态得到元素为Pi,j的状态转移矩阵P;每个元素Pi,j都能用转移频数矩阵N中的元素Ni,j表示出;预测的风场在th+1时刻的出力都为一个K维向量中的任一元素为对应的区间的中值,而对应的区间的上下限都能由当前状态及状态转移矩阵计算得到,由此得到th+1时刻风电出力的区间分布。
3.如权利要求2所述考虑大规模风电接入的在线风险评估方法,其特征是,风电出力为对应的概率为Pj,i,将其作为风险评估中的不确定场景概率Pr(Xi),其中Xi为th+1时刻风电出力的第i个不确定区间。
4.如权利要求1所述考虑大规模风电接入的在线风险评估方法,其特征是,所述步骤二的具体方法为,风场按恒功率因数控制,将系统的有功和无功的变化量用常规潮流计算迭代公式表示,对该公式雅可比矩阵求逆,得到节点电压向量的变化矩阵,由节点电压向量的变化矩阵得到预估的节点电压幅值和相角向量。
5.如权利要求1所述考虑大规模风电接入的在线风险评估方法,其特征是,所述步骤四的具体方法为,根据戴维南等值理论得到负荷节点的戴维南等值参数,应用PV曲线中电压在鼻尖点处只有一个取值的原理计算负荷节点的负荷裕度λcr,将风电并网节点处理成恒功率因数控制的PQ节点;假设电力系统共包含N个PQ节点,将第m个PQ节点作为观察节点,计算出其负荷裕度为λcr,m,m=1,2,3…N,计算所有PQ节点的负荷裕度,选择其中负荷裕度最小值作为系统的负荷裕度λmin;结合th+1时刻第i个区间的风电出力,得到对应风电出力第i个区间时系统的最小负荷裕度λmin,i,将λmin,i平方后取倒数求得th+1时刻第i风电出力预测区间的严重度。
6.如权利要求1所述考虑大规模风电接入的在线风险评估方法,其特征是,风险指标的求解方法为,假设在th+1时刻风电出力共有K个不确定区间,每个不确定区间的概率与该区间对应的严重度相乘后累加得到th时刻到th+1时刻周期内的风险指标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610052437.9A CN105610192B (zh) | 2016-01-26 | 2016-01-26 | 考虑大规模风电接入的在线风险评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610052437.9A CN105610192B (zh) | 2016-01-26 | 2016-01-26 | 考虑大规模风电接入的在线风险评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105610192A true CN105610192A (zh) | 2016-05-25 |
CN105610192B CN105610192B (zh) | 2019-04-09 |
Family
ID=55989821
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610052437.9A Expired - Fee Related CN105610192B (zh) | 2016-01-26 | 2016-01-26 | 考虑大规模风电接入的在线风险评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105610192B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106682407A (zh) * | 2016-12-19 | 2017-05-17 | 国网上海市电力公司 | 基于戴维南等值和支路传输功率极限的电压稳定评估方法 |
CN107332251A (zh) * | 2017-06-15 | 2017-11-07 | 清华大学 | 一种双馈风机风电场参与电力系统动态电压控制的方法 |
CN107464007A (zh) * | 2016-06-02 | 2017-12-12 | 南京理工大学 | 基于马尔科夫理论和比例分配原理的连续时段概率潮流预测方法 |
CN108092262A (zh) * | 2017-12-07 | 2018-05-29 | 广东电网有限责任公司惠州供电局 | 一种计及风电机组脱网时序的故障电流双层迭代算法 |
CN108306287A (zh) * | 2018-02-09 | 2018-07-20 | 清华大学 | 一种含风电电力系统的小干扰失稳风险识别方法及系统 |
CN109274124A (zh) * | 2018-11-22 | 2019-01-25 | 国网黑龙江省电力有限公司电力科学研究院 | 基于场景马尔科夫法的风电就地消纳能力预测方法 |
CN109934486A (zh) * | 2019-03-11 | 2019-06-25 | 长沙理工大学 | 一种提取故障耦合传播过程中多能流负荷裕度演变及其交互特征的方法 |
CN110148935A (zh) * | 2019-05-10 | 2019-08-20 | 东北电力大学 | 一种基于长短期记忆神经网络的电网静态安全裕度评估方法 |
CN110460085A (zh) * | 2019-08-28 | 2019-11-15 | 华北电力大学(保定) | 一种考虑风电和负荷特性对电力系统影响的方法 |
CN110543693A (zh) * | 2019-08-13 | 2019-12-06 | 国网山西省电力公司电力科学研究院 | 一种计及运行风险的鲁棒机组组合方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103246806A (zh) * | 2013-04-25 | 2013-08-14 | 浙江大学 | 一种含风电场电力系统的运行风险评估方法 |
CN104659782A (zh) * | 2015-03-20 | 2015-05-27 | 太原理工大学 | 考虑负荷波动极限的电力系统电压稳定性风险评估方法 |
CN104899798A (zh) * | 2015-06-30 | 2015-09-09 | 天津大学 | 一种考虑旋转备用的风电接入系统的暂态风险控制方法 |
-
2016
- 2016-01-26 CN CN201610052437.9A patent/CN105610192B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103246806A (zh) * | 2013-04-25 | 2013-08-14 | 浙江大学 | 一种含风电场电力系统的运行风险评估方法 |
CN104659782A (zh) * | 2015-03-20 | 2015-05-27 | 太原理工大学 | 考虑负荷波动极限的电力系统电压稳定性风险评估方法 |
CN104899798A (zh) * | 2015-06-30 | 2015-09-09 | 天津大学 | 一种考虑旋转备用的风电接入系统的暂态风险控制方法 |
Non-Patent Citations (4)
Title |
---|
JIAN-HONG LIU等: "Wide-Area Measurement-Based Voltage Stability Indicators by Modified Coupled Single-Port Models", 《IEEE TRANSACTIONS ON POWER SYSTEMS》 * |
JINGWEN SUN等: "Wind Power Forecasting Based on a Markov Chain Model of Variation", 《2015 12TH INTERNATIONAL CONFERENCE ON FUZZY SYSTEMS AND KNOWLEDGE DISCOVERY》 * |
张元等: "风电场等值建模研究综述", 《电力系统保护与控制》 * |
户秀琼等: "静态电压稳定预防控制的风险评估", 《电力自动化设备》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107464007A (zh) * | 2016-06-02 | 2017-12-12 | 南京理工大学 | 基于马尔科夫理论和比例分配原理的连续时段概率潮流预测方法 |
CN106682407B (zh) * | 2016-12-19 | 2019-05-31 | 国网上海市电力公司 | 基于戴维南等值和支路传输功率极限的电压稳定评估方法 |
CN106682407A (zh) * | 2016-12-19 | 2017-05-17 | 国网上海市电力公司 | 基于戴维南等值和支路传输功率极限的电压稳定评估方法 |
CN107332251A (zh) * | 2017-06-15 | 2017-11-07 | 清华大学 | 一种双馈风机风电场参与电力系统动态电压控制的方法 |
CN107332251B (zh) * | 2017-06-15 | 2019-09-27 | 清华大学 | 一种双馈风机风电场参与电力系统动态电压控制的方法 |
CN108092262A (zh) * | 2017-12-07 | 2018-05-29 | 广东电网有限责任公司惠州供电局 | 一种计及风电机组脱网时序的故障电流双层迭代算法 |
CN108092262B (zh) * | 2017-12-07 | 2019-10-18 | 广东电网有限责任公司惠州供电局 | 一种计及风电机组脱网时序的故障电流双层迭代算法 |
CN108306287A (zh) * | 2018-02-09 | 2018-07-20 | 清华大学 | 一种含风电电力系统的小干扰失稳风险识别方法及系统 |
CN109274124B (zh) * | 2018-11-22 | 2022-08-12 | 国网黑龙江省电力有限公司电力科学研究院 | 基于场景马尔科夫法的风电就地消纳能力预测方法 |
CN109274124A (zh) * | 2018-11-22 | 2019-01-25 | 国网黑龙江省电力有限公司电力科学研究院 | 基于场景马尔科夫法的风电就地消纳能力预测方法 |
CN109934486A (zh) * | 2019-03-11 | 2019-06-25 | 长沙理工大学 | 一种提取故障耦合传播过程中多能流负荷裕度演变及其交互特征的方法 |
CN110148935A (zh) * | 2019-05-10 | 2019-08-20 | 东北电力大学 | 一种基于长短期记忆神经网络的电网静态安全裕度评估方法 |
CN110148935B (zh) * | 2019-05-10 | 2021-12-28 | 东北电力大学 | 基于长短期记忆神经网络的电网静态安全裕度评估方法 |
CN110543693A (zh) * | 2019-08-13 | 2019-12-06 | 国网山西省电力公司电力科学研究院 | 一种计及运行风险的鲁棒机组组合方法 |
CN110543693B (zh) * | 2019-08-13 | 2022-08-23 | 国网山西省电力公司电力科学研究院 | 一种计及运行风险的鲁棒机组组合方法 |
CN110460085A (zh) * | 2019-08-28 | 2019-11-15 | 华北电力大学(保定) | 一种考虑风电和负荷特性对电力系统影响的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105610192B (zh) | 2019-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105610192A (zh) | 考虑大规模风电接入的在线风险评估方法 | |
CN105741025B (zh) | 基于风电波动在线风险评估的预防控制方法 | |
CN103036230B (zh) | 一种基于工程应用的交直流混联大电网的动态等值方法 | |
CN103279804B (zh) | 超短期风电功率的预测方法 | |
EP2537222B1 (en) | Energy generating system and control thereof | |
CN101741086B (zh) | 基于pmu数据的稳定计算模型的综合评估方法 | |
CN108183512A (zh) | 一种接入新能源的电力系统的可靠性评估方法 | |
CN104599189A (zh) | 一种计入电力系统运行方式的电网规划方案风险评估方法 | |
CN110417011A (zh) | 一种基于互信息与迭代随机森林的在线动态安全评估方法 | |
Li | Framework of probabilistic power system planning | |
CN105098763A (zh) | 一种基于wams和ems的暂态电压稳定在线评估方法 | |
CN103258103A (zh) | 基于偏最小二乘回归的戴维南等值参数辨识方法 | |
Mo et al. | Impact of aging and performance degradation on the operational costs of distributed generation systems | |
CN104052058A (zh) | 一种基于马尔可夫链蒙特卡罗法的系统谐波概率评估方法 | |
CN101923685A (zh) | 一种基于断线故障概率预测的电力切负荷决策系统及方法 | |
Kaplan et al. | A novel method based on Weibull distribution for short-term wind speed prediction | |
CN103034787A (zh) | 一种微电网状态估计的方法 | |
CN110224392A (zh) | 一种用于分析含风电系统电压稳定概率的无迹变换方法 | |
CN105447620A (zh) | 一种电能量缺失值自动处理的方法 | |
Kim et al. | A probabilistic approach to potential estimation of renewable energy resources based on augmented spatial interpolation | |
CN103632314A (zh) | 基于概率统计的广义节点特征建模方法 | |
Kou et al. | PQ curve based voltage stability analysis considering wind power | |
Rahimi et al. | Quasi-Steady-State computation of voltage flicker with cloud motion simulator | |
Zeraati et al. | Meter placement algorithms to enhance distribution systems state estimation: Review, challenges and future research directions | |
Xiao et al. | Comprehensive Evaluation Index System of Distribution Network for Distributed Photovoltaic Access |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into 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: 20190409 |