CN105743093B - 一种考虑风光出力不确定性的配电网潮流计算方法 - Google Patents

一种考虑风光出力不确定性的配电网潮流计算方法 Download PDF

Info

Publication number
CN105743093B
CN105743093B CN201610230133.7A CN201610230133A CN105743093B CN 105743093 B CN105743093 B CN 105743093B CN 201610230133 A CN201610230133 A CN 201610230133A CN 105743093 B CN105743093 B CN 105743093B
Authority
CN
China
Prior art keywords
mtd
msub
output
minizone
mrow
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
Application number
CN201610230133.7A
Other languages
English (en)
Other versions
CN105743093A (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.)
Hefei University of Technology
Original Assignee
Hefei University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hefei University of Technology filed Critical Hefei University of Technology
Priority to CN201610230133.7A priority Critical patent/CN105743093B/zh
Publication of CN105743093A publication Critical patent/CN105743093A/zh
Application granted granted Critical
Publication of CN105743093B publication Critical patent/CN105743093B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • H02J3/04Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
    • H02J3/06Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • 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
    • 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]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Power Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computing Systems (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Wind Motors (AREA)

Abstract

本发明公开了一种考虑风光出力不确定性的配电网潮流计算方法,其特征是按如下步骤进行:利用矩阵判断法得到盲数各区间的可信度;将风力发电机的出力和光伏电站的出力用盲数进行表示,由小区间和可信度组成,为简化计算,将小区间内的抽样频率代替小区间的可信度;将抽样得到的风光出力加入到配电网并结合重复潮流计算,输出潮流计算结果。本发明能够描述电网规划中不确定因素的多种不确定性,能更加准确具体的反映电网的状态。

Description

一种考虑风光出力不确定性的配电网潮流计算方法
技术领域
本发明涉及电力系统分析技术领域,具体涉及一种考虑风光出力不确定性的配电网潮流计算的方法。
背景技术
潮流计算是对电力系统运行状态进行量化分析的基础。微网的大力发展对传统的电力系统规划提出了挑战。由于风速、光照强度的随机性和波动性,风电场和光伏电站的出力存在着间歇性和随机性等特点。在电力系统规划和运行分析中,传统的确定性潮流计算已经不再适用,需要进行不确定性的潮流计算。
目前对于不确定信息的处理主要有两种方法:利用多场景技术和数学方法建立不确定信息的模型。多场景技术用来处理多种不确定性信息取得很好的效果,但其缺陷是难以准确地量化场景发生的可能性或由于场景过多而造成求解上的困难。利用数学来描述不确定信息主要有四种方法:概率统计、灰色数学、模糊数学和未确知数学。虽然概率统计、区间数和模糊数能够较好的反映不确定信息,但是它们只能反映单一的不确定信息。而在实际工程应用中,信息的不确定性往往不是单一的,而是具有多种不确定性,通常包含灰性、模糊性、随机性等两种或两种以上的不确定性。
发明内容
本发明为克服上述现有技术中存在的不足之处,提供一种考虑风光出力不确定性的配电网潮流计算方法,以期能描述电网中同时具有随机性、灰性、模糊性等多种不确定性的信息,从而能在电网的运行规划中给出更加准确全面的信息,并对电网的薄弱环节做出量化分析。
为了达到上述目的,本发明所采用的技术方案为:
本发明一种考虑风光出力不确定性的配电网潮流计算方法,所述配电网接入有风力发电机和光伏电站;其特点是,所述配电网潮流计算方法是按如下步骤进行:
步骤一、利用矩阵判断法得到盲数各区间的可信度:
步骤1.1、将所述风力发电机的出力分为m个小区间,记为{p1,p2,…,pi,…,pm},pi表示第i个小区间;1≤i≤m;
步骤1.2、对所述风力发电机的出力进行统计,获得第i个小区间pi对第j个小区间pj的可能性程度以及第j个小区间pj对第i个小区间pi的可能性程度1≤j≤m;
步骤1.3、将可能性程度进行比较,若则令判断值则令判断值a∈(1,9];
步骤1.4、利用式(1)构造判断矩阵S,其中第i行第j列的元素Sij为:
步骤1.5、求取所述判断矩阵S的特征值,并求出最大特征值所对应的特征向量,对最大特征值所对应的特征向量进行归一化处理,获得的归一化处理结果作为m个小区间的可信度,记为{α12,…,αi,…,αm};αi表示第i个小区间pi所对应的可信度;
步骤二、将风力发电机的出力用盲数表示;
利用式(2)将所述风力发电机的出力PWTG用盲数表示为:
式(2)中,m为盲数阶数,当i≠j时,有αi≠αj,且
步骤三、按照步骤一和步骤二获得所述光伏电站的出力PPVG用盲数表示为:
式(2)中,n为盲数阶数,{x1,x2,…,xk,…,xn}表示光伏电站的出力所划分出的n个小区间,xk表示第k个小区间,βk表示第k个小区间所对应的可信度;1≤k≤n;
步骤四、求解所述配电网的不确定性潮流:
步骤4.1、设定在所述风力发电机的第i个小区间pi内的抽样次数为li,且li=N×αi;m个小区间内的总的抽样次数为N;并在每个小区间内进行li次随机抽样;定义抽样次数为q,并初始化q=1;
步骤4.2、设定在所述光伏电站的第k个小区间xk内的抽样次数为Lk,且Lk=N×βk;n个小区间内的总的抽样次数为N;并在每个小区间内进行Lk次随机抽样;
步骤4.3、将第q次抽样得到的风力发电机的出力代入配电网进行潮流计算,获得第q次抽样的风力发电机接入时配电网的电压、电流和功率;
步骤4.4、将第q次抽样得到的光伏电站的出力代入配电网进行潮流计算,获得第q次抽样的光伏电站接入时配电网的电压、电流和功率;
步骤4.5、将第q次抽样得到的风力发电机的出力和光伏电站的出力代入配电网则获得第q次抽样的配电网的电压、电流和功率;
步骤4.6、将q+1赋值给q,并返回步骤4.3执行,直到q=N为止,从而获得风力发电机接入时配电网的电压、电流和功率、光伏电站接入时配电网的电压、电流和功率以及风力发电机和光伏电站同时接入时配电网的电压、电流和功率,即为潮流计算结果;
步骤五、利用盲数的求均值计算将盲数潮流的小区间用相应的数值代替,从而对潮流计算结果进行简化处理。
本发明所述的配电网潮流计算方法的特点也在于,
所述步骤4.5中,风力发电机和光伏电站同时接入配电网时,风力发电机的出力和光伏电站的出力由两个不同的盲数表示,同时对风力发电机出力和光伏电站出力进行区间进行抽样时,利用多个盲数的运算法则对风力发电机的出力和光伏电站的出力进行运算,得到新盲数{y1,y2,…,yh,…,yt}用于表示风力发电机和光伏电站总的出力,其中,yh表示第h个小区间;令{δ12,…,δh,…,δt}表示t个区间对应的可信度;
设定风力发电机的出力和光伏电站的出力在第h个小区间yh内的抽样次数为Lh,且Lh=N×δh;t个小区间内的总的抽样次数为N;并在每个小区间内进行相应次随机抽样;
将第q次抽样得到的风力发电机的出力和光伏电站的出力进行盲数计算,获得的结果代入配电网进行潮流计算,从而获得第q次抽样的风力发电机和光伏电站接入时配电网的电压、电流和功率。
与已有技术相比,本发明的有益效果体现在:
1本发明针对电网不确定信息的多样性,利用盲数理论建立了不确定信息的数学模型,由于盲数能够描述信息的多种不确定性,从而使所建立的不确定数学模型能够更加准确的描述不确定信息,能更加准确具体的反映电网的状态,更加符合实际情况。
2本发明利用盲数建立风力发电机出力和光伏电站出力的的不确定性模型,根据实际经验和测量的数据,得到风力发电机出力和光伏电站出力的小区间和可信度,避免了由于分布估计时引入的不确定性,使得结果更加准确的反映实际的情况。
3本发明对盲数表示的风力发电机出力和光伏电站出力的各个小区间进行多次抽样,结合重复潮流计算法,进行盲数潮流计算,简化了计算过程,可以快速地得到潮流结果。
4本发明针对风力发电机和光伏电站同时接入系统的情况,利用多个盲数的计算法则对风力发电机的出力和光伏电站的出力进行计算得到一个表示风力发电机和光伏电站出力的新盲数,并结合重复潮流法进行潮流计算,简化了计算过程。
5本发明的盲数潮流是由小区间和可信度组成,其包含的信息量比较大,能够给出较为丰富和准确的线路潮流信息。
附图说明
图1为本发明方法流程示意图。
具体实施方式
本实施例中,如图1所示,配电网接入有风力发电机和光伏电站;一种考虑风光出力不确定性的配电网潮流计算方法是按如下步骤进行:
步骤一、利用矩阵判断法得到盲数各区间的可信度:
步骤1.1、将风力发电机的出力分为m个小区间,记为{p1,p2,…,pi,…,pm},pi表示第i个小区间;1≤i≤m;
步骤1.2、对风力发电机的出力进行统计,获得第i个小区间pi对第j个小区间pj的可能性程度以及第j个小区间pj对第i个小区间pi的可能性程度1≤j≤m;
步骤1.3、将可能性程度进行比较,若则令判断值则令判断值a∈(1,9];
步骤1.4、利用式(1)构造判断矩阵S,其中第i行第j列的元素Sij为:
步骤1.5、求取判断矩阵S的特征值,并求出最大特征值所对应的特征向量,对最大特征值所对应的特征向量进行归一化处理,获得的归一化处理结果作为m个小区间的可信度,记为{α12,…,αi,…,αm};αi表示第i个小区间pi所对应的可信度;
根据矩阵判断法确定的含分布式电源的配电网不确定因素的盲数模型能够更加准确客观的反映信息的不确定性。
步骤二、将风力发电机的出力用盲数表示;
利用式(2)将风力发电机的出力PWTG用盲数表示为:
式(2)中,m为盲数阶数,当i≠j时,有αi≠αj,且
在实际工程应用中,信息的不确定性往往不是单一,而是具有多种不确定性,通常包含灰性、模糊性、随机性等两种或两种以上的不确定性。盲数是一种新的不确定性信息的数学表达形式,其优势在于能够描述信息的多种不确定性,可以最大程度地表示电网规划中的不确定因素。m越大,盲数的阶数越大,对不确定信息的描述也越准确,信息的不确定性程度高。
步骤三、按照步骤一和步骤二获得光伏电站的出力PPVG用盲数表示为:
式(2)中,n为盲数阶数,{x1,x2,…,xk,…,xn}表示光伏电站的出力所划分出的n个小区间,xk表示第k个小区间,βk表示第k个小区间所对应的可信度;1≤k≤n;
步骤四、求解配电网的不确定性潮流:
步骤4.1、设定在风力发电机的第i个小区间pi内的抽样次数为li,且li=N×αi;m个小区间内的总的抽样次数为N;并在每个小区间内进行li次随机抽样;定义抽样次数为q,并初始化q=1;
步骤4.2、设定在光伏电站的第k个小区间xk内的抽样次数为Lk,且Lk=N×βk;n个小区间内的总的抽样次数为N;并在每个小区间内进行Lk次随机抽样;
步骤4.3、将第q次抽样得到的风力发电机的出力代入配电网进行潮流计算,获得第q次抽样的风力发电机接入配电网时配电网的电压、电流和功率,;
进行潮流计算时需要给出一定的约束条件保证电网的正常运行。式(4)、(5)、(6)、(7)、(8)给出了潮流约束条件:
PGfmin<PGf<PGfmax (6)
QGfmin<QGf<QGfmax (7)
Ufmin<Uf<Ufmax (8)
式(4)、(5)、(6)、(7)、(8)、(9)中,分别表示接在节点f的分布式电源的有功和无功出力的盲数表达形式;分别表示f节点有功负荷和无功负荷的盲数表达形式;M表示配电网节点的个数,Uf和Ug表示f和g节点的电压值,Gfg和Bfg分别表示电导和电纳,δfg表示节点f和g之间的相角差;PGfmin和PGfmax分别表示分布式电源允许发出的最小有功功率和最大有功功率;QGfmin和QGfmax分别表示分布式电源允许发出的最小无功功率和最大无功功率;Ufmin和Ufmax分别表示f节点的电压下限和电压上限;
步骤4.4、将第q次抽样得到的光伏电站的出力代入配电网进行潮流计算,,获得第q次抽样的光伏电站接入配电网时配电网的电压、电流和功率;
步骤4.5、将第q次抽样得到的风力发电机的出力和光伏电站的出力代入配电网进行潮流计算,获得第q次抽样的风力发电机和光伏电站同时接入配电网时配电网的电压、电流和功率;
步骤4.6、将q+1赋值给q,并返回步骤4.3执行,直到q=N为止,从而获得风力发电机接入时配电网的电压、电流和功率、光伏电站接入时配电网的电压、电流和功率以及风力发电机和光伏电站同时接入时配电网的电压、电流和功率,即为潮流计算结果;
步骤五、利用盲数的求均值计算将盲数潮流的小区间用相应的数值代替,从而对潮流计算结果进行简化处理。
由于盲数潮流由大量的小区间和可信度组成,会存在大量的区间重叠,如果对各个区间之间的重叠部分进行合成,其工作量很大,需要对潮流结果进行简化处理。首先求得盲数潮流区间的最大和最小值,在区间上通过等距取点,得到很多小区间。通过统计得到各个小区间的可信度,根据盲数的求均值计算法则,利用点代替小区间,从而可以较容易地得到该线路的潮流可能分布情况。
具体实施中,步骤4.5中的风力发电机和光伏电站同时接入配电网时,风力发电机的出力和光伏电站的出力由两个不同的盲数表示,同时对风力发电机出力和光伏电站出力进行区间进行抽样时,利用多个盲数的运算法则对风力发电机的出力和光伏电站的出力进行运算得到新的盲数{y1,y2,…,yh,…,yt},表示风力发电机和光伏电站总的出力,yh表示第h个小区间,{δ12,…,δh,…,δt}表示t个区间对应的可信度,设定风力发电机的出力和光伏电站的出力在第h个小区间yh内的抽样次数为Lh,且Lh=N×δh;t个小区间内的总的抽样次数为N;并在每个小区间内进行Lh次随机抽样;将第q次抽样得到的风力发电机的出力和光伏电站的出力进行盲数计算,从而获得的结果代入配电网进行潮流计算。
当多个盲数进行潮流计算时,假设A和B为两个盲数,盲数的计算法则如下:
当A和B进行加、减、乘、除四则运算中的某一种运算时,其运算结果用A关于B的可能值带“*”矩阵和A关于B的可信度带边积矩阵进行表示,分别表一和表二所示。A关于B的可能值带“*”矩阵中的“*”表示加、减、乘、除四则运算中的一种。A与B的可信度带边积矩阵中的“*”为乘积运算,在可能值带“*”矩阵中,相同的元素需要进行合并,对应的可信度带边积矩阵中的可信度相累加。
表一 盲数A关于盲数B的可能值带“*”矩阵表
表二 盲数A关于盲数B的可信度带边积矩阵表
当多个盲数进行计算时会使盲数的阶数快速的增加,给计算带来了极大的不便。为了简化计算过程,需要在计算过程中对盲数进行降阶处理,盲数的降阶一般有两种方法:合并相交区间和舍弃小可信度区间。

Claims (2)

1.一种考虑风光出力不确定性的配电网潮流计算方法,所述配电网接入有风力发电机和光伏电站;其特征是,所述配电网潮流计算方法是按如下步骤进行:
步骤一、利用矩阵判断法得到盲数各区间的可信度:
步骤1.1、将所述风力发电机的出力分为m个小区间,记为{p1,p2,…,pi,…,pm},pi表示第i个小区间;1≤i≤m;
步骤1.2、对所述风力发电机的出力进行统计,获得第i个小区间pi对第j个小区间pj的可能性程度以及第j个小区间pj对第i个小区间pi的可能性程度1≤j≤m;
步骤1.3、将可能性程度进行比较,若则令判断值则令判断值a∈(1,9];
步骤1.4、利用式(1)构造判断矩阵S,其中第i行第j列的元素Sij为:
<mrow> <msub> <mi>S</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>=</mo> <mfrac> <msubsup> <mi>K</mi> <mi>j</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> </msubsup> <msubsup> <mi>K</mi> <mi>i</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>)</mo> </mrow> </msubsup> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
步骤1.5、求取所述判断矩阵S的特征值,并求出最大特征值所对应的特征向量,对最大特征值所对应的特征向量进行归一化处理,获得的归一化处理结果作为m个小区间的可信度,记为{α12,…,αi,…,αm};αi表示第i个小区间pi所对应的可信度;
步骤二、将风力发电机的出力用盲数表示;
利用式(2)将所述风力发电机的出力PWTG用盲数表示为:
<mrow> <msub> <mi>P</mi> <mrow> <mi>W</mi> <mi>T</mi> <mi>G</mi> </mrow> </msub> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <msub> <mi>p</mi> <mn>1</mn> </msub> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>p</mi> <mn>2</mn> </msub> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mn>...</mn> </mtd> <mtd> <mn>...</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>p</mi> <mi>i</mi> </msub> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mi>i</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>...</mn> </mtd> <mtd> <mn>...</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>p</mi> <mi>m</mi> </msub> </mtd> <mtd> <mrow> <mi>&amp;alpha;</mi> <mi>m</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
式(2)中,m为盲数阶数,当i≠j时,有αi≠αj,且
步骤三、按照步骤一和步骤二获得所述光伏电站的出力PPVG用盲数表示为:
<mrow> <msub> <mi>P</mi> <mrow> <mi>P</mi> <mi>V</mi> <mi>G</mi> </mrow> </msub> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <msub> <mi>x</mi> <mn>1</mn> </msub> </mtd> <mtd> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>x</mi> <mn>2</mn> </msub> </mtd> <mtd> <msub> <mi>&amp;beta;</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mn>...</mn> </mtd> <mtd> <mn>...</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>x</mi> <mi>k</mi> </msub> </mtd> <mtd> <msub> <mi>&amp;beta;</mi> <mi>k</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>...</mn> </mtd> <mtd> <mn>...</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>x</mi> <mi>n</mi> </msub> </mtd> <mtd> <mrow> <mi>&amp;beta;</mi> <mi>n</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
式(2)中,n为盲数阶数,{x1,x2,…,xk,…,xn}表示光伏电站的出力所划分出的n个小区间,xk表示第k个小区间,βk表示第k个小区间所对应的可信度;1≤k≤n;
步骤四、求解所述配电网的不确定性潮流:
步骤4.1、设定在所述风力发电机的第i个小区间pi内的抽样次数为li,且li=N×αi;m个小区间内的总的抽样次数为N;并在每个小区间内进行li次随机抽样;定义抽样次数为q,并初始化q=1;
步骤4.2、设定在所述光伏电站的第k个小区间xk内的抽样次数为Lk,且Lk=N×βk;n个小区间内的总的抽样次数为N;并在每个小区间内进行Lk次随机抽样;
步骤4.3、将第q次抽样得到的风力发电机的出力代入配电网进行潮流计算,获得第q次抽样的风力发电机接入时配电网的电压、电流和功率;
步骤4.4、将第q次抽样得到的光伏电站的出力代入配电网进行潮流计算,获得第q次抽样的光伏电站接入时配电网的电压、电流和功率;
步骤4.5、将第q次抽样得到的风力发电机的出力和光伏电站的出力代入配电网则获得第q次抽样的配电网的电压、电流和功率;
步骤4.6、将q+1赋值给q,并返回步骤4.3执行,直到q=N为止,从而获得风力发电机接入时配电网的电压、电流和功率、光伏电站接入时配电网的电压、电流和功率以及风力发电机和光伏电站同时接入时配电网的电压、电流和功率,即为潮流计算结果;
步骤五、利用盲数的求均值计算将盲数潮流的小区间用相应的数值代替,从而对潮流计算结果进行简化处理。
2.根据权利要求1所述的配电网潮流计算方法,其特征是,所述步骤4.5中,风力发电机和光伏电站同时接入配电网时,风力发电机的出力和光伏电站的出力由两个不同的盲数表示,同时对风力发电机出力和光伏电站出力进行区间抽样时,利用多个盲数的运算法则对风力发电机的出力和光伏电站的出力进行运算,得到新盲数{y1,y2,…,yh,…,yt}用于表示风力发电机和光伏电站总的出力,其中,yh表示第h个小区间;令{δ12,…,δh,…,δt}表示t个区间对应的可信度;
设定风力发电机的出力和光伏电站的出力在第h个小区间yh内的抽样次数为Lh,且Lh=N×δh;t个小区间内的总的抽样次数为N;并在每个小区间内进行相应次随机抽样;
将第q次抽样得到的风力发电机的出力和光伏电站的出力进行盲数计算,获得的结果代入配电网进行潮流计算,从而获得第q次抽样的风力发电机和光伏电站接入时配电网的电压、电流和功率。
CN201610230133.7A 2016-04-13 2016-04-13 一种考虑风光出力不确定性的配电网潮流计算方法 Active CN105743093B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610230133.7A CN105743093B (zh) 2016-04-13 2016-04-13 一种考虑风光出力不确定性的配电网潮流计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610230133.7A CN105743093B (zh) 2016-04-13 2016-04-13 一种考虑风光出力不确定性的配电网潮流计算方法

Publications (2)

Publication Number Publication Date
CN105743093A CN105743093A (zh) 2016-07-06
CN105743093B true CN105743093B (zh) 2018-05-18

Family

ID=56254284

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610230133.7A Active CN105743093B (zh) 2016-04-13 2016-04-13 一种考虑风光出力不确定性的配电网潮流计算方法

Country Status (1)

Country Link
CN (1) CN105743093B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107145707B (zh) * 2017-04-01 2019-12-27 三峡大学 一种计及光伏出力不确定性和全寿命周期成本的配电网变压器规划方法
CN109146249B (zh) * 2018-07-23 2021-02-09 国网江苏省电力有限公司 一种配电网可靠性预估方法、装置和设备
CN113725896B (zh) * 2021-08-11 2023-03-17 武汉大学 基于不确定性风险计算的电网新能源接入方法
CN114336637B (zh) * 2022-01-10 2023-04-07 合肥工业大学 一种考虑风光出力相关性的配电网证据理论潮流计算方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009065778A (ja) * 2007-09-06 2009-03-26 Tokyo Electric Power Co Inc:The 電力系統監視制御システム
CN201286019Y (zh) * 2008-09-19 2009-08-05 上海万德风力发电股份有限公司 一种与电网互联式110kw风光互补发电系统

Also Published As

Publication number Publication date
CN105743093A (zh) 2016-07-06

Similar Documents

Publication Publication Date Title
CN107425520B (zh) 一种含节点注入功率不确定性的主动配电网三相区间状态估计方法
CN108365822B (zh) 光伏电池在不同光照温度条件下工作特性的预测方法及系统
CN106874630A (zh) 基于电量消纳的区域电网新能源开发潜力评估方法
CN105743093B (zh) 一种考虑风光出力不确定性的配电网潮流计算方法
Fan et al. Review of uncertainty modeling for optimal operation of integrated energy system
CN108306303A (zh) 一种考虑负荷增长和新能源出力随机的电压稳定评估方法
CN109802394A (zh) 一种计及分布式电源与电动汽车接入的概率潮流计算方法
CN107145707A (zh) 一种计及光伏出力不确定性和全寿命周期成本的配电网变压器规划方法
CN106548410B (zh) 一种含分布式电源的配电网电压不平衡概率评估方法
CN105656031A (zh) 基于高斯混合分布特征的含风电电力系统安全风险评估方法
CN113361946B (zh) 一种基于分布式光伏并网系统的电能质量评估方法及装置
CN107453396A (zh) 一种分布式光伏电源出力的多目标优化调度方法
CN108599239A (zh) 一种下垂控制型孤岛微电网电压质量概率评估方法
Wu et al. Reliability evaluation of the solar power system based on the Markov chain method
CN105512775A (zh) 一种光伏发电系统的功率预测方法
Zhang et al. Data-driven sizing planning of renewable distributed generation in distribution networks with optimality guarantee
CN107171365A (zh) 基于场景解耦和异步迭代的多目标随机动态经济调度方法
CN106058863A (zh) 一种基于随机响应面法的随机最优潮流计算方法
CN107039981A (zh) 一种拟直流线性化概率最优潮流计算方法
Ni et al. A review of line loss analysis of the low-voltage distribution system
Shendryk et al. Decision Support System for Efficient Energy Management of MicroGrid with Renewable Energy Sources
CN107276093A (zh) 基于场景削减的电力系统概率潮流计算方法
CN110380419A (zh) 一种混合随机和区间变量的不确定潮流计算方法
CN107918920B (zh) 多光伏电站的出力关联性分析方法
CN116137445A (zh) 基于分布式电源的孤岛配网方法、装置、设备及介质

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