开云网址完整)电力系统潮流计算方法分析1、电力系统潮流分析基于牛拉法和保留非线性的随机潮流姓名:*学号:*1潮流算法简介1.1常规潮流计算常规的潮流计算是在确定的状态下.即:通过已知运行条件(比如节点功率或网络结构等)得到系统的运行状态(比如所有节点的电压值与相角、所有支路上的功率分布和损耗等)。常规潮流算法中的一种普遍采用的方法是牛顿-拉夫逊法.当初始值和方程的精确解足够接近时,该方法可以在很短时间内收敛.下面简要介绍该方法。1.1。1牛顿拉夫逊方法原理对于非线T),在待求量x初次的估计值x(0)附近,用泰勒级数(忽略二阶和以上的高阶项)表示它,可获得如式(1-2)的线性化变换后的方程组,该方程组被称为修正方程组。
3、设系统的网络中一共设有n个节点当采用直角坐标系解决潮流方程,此时待解电压和导纳如下式:(1-7)1-8)平衡节点的电压是已知的,平衡节点表示如下.V=e+jfnnn(完整)电力系统潮流计算方法分析除了平衡节点以外的所有2(n-1)个节点是需要求解的量。每个节点可列出两个方程式假定系统中前m个节点为P一Q节点,第m+1到n1个节点为P-V节点。对于PQ节点,P和Q的值是ii固定的,对于PV节点,P和V的值是固定的。iiaP=Pe工(GeB/)/工(G/+Be)=0,.iisijjjjjjjjj(1-9)Q=Q-f工(GeBf)+e乞(Gf+Be)=0i=1,2,mijjijjjijjijjji-
9、标形式的公式不存在截断误差,因此为了减小计算误差,本文以直角坐标形式的牛拉法为基础编写了保留非线性潮流计算方法的程序。迭代公式为:x(k+i)=J-iy(x(o)ys+y(Ax(k)(114)迭代过程和牛拉法相类似,流程图如下所示:图1。2保留非线蒙特卡罗模拟原理(完整)电力系统潮流计算方法分析蒙特卡罗模拟方法的思想是,是当求解问题是一不确定事件的平均值时,我们通过构建模型并采用某特定的“实验”,就可以实验中此事件发生的频率去估算概率.1.2。2蒙特卡罗模拟步骤1)根据不同新能源的特点建立新能源输出功率的样本,规模为N;2)将得到的N个样本值带入对应接
10、入新能源的各节点,得到接入光伏后的各节点的值。3)按照1。1所述的牛顿拉夫逊法进行确定性潮流计算,得到N组关于节点的电压,支路功率与网损的数据等。4)运用数学上的统计原理,可以求出输出变量的分布情况.1。3拉丁超立方采样法1。3。1拉丁超立方采样原理拉丁超立方采样由M.D.McKay、R。J.Beckman和W。J。Conover在1979年提出,它通过分层采样使采样点能够覆盖到整个随机变量的分布范围。该方法分成两步:1)采样:所有的输入变量可以通过分层采样,使得样本点更加准确均匀的分布;2)排列:改变初次采样得到的样本数据的顺序,令变量数据之间的关联程度最小,或者通过排序达到指定的相关系数。
11、1.3.2拉丁超立方采样优点1)可以使采样得到的数据较为全面地覆盖变量所分布的范围,同时分层使得采样时不会再采到一样或相似的数据,更准确地体现变量的总体情况,同时减小了样本规模。一些文献证明了拉丁超立方采样与简单随机采样在采样规模同是M时,两种方法抽取到的变量假设是独立的,那么它们的联合覆盖空间百分比平均值表示如下:M-12P=xlOO%116)1LM_M-112P=xlOO%m_M+1_可以看出,当M大于等于2时,一式大于二式,表明拉丁超立方采样比随机采样覆盖的范围(完整)电力系统潮流计算方法分析(完整)电力系统潮流计算方法分析a一rand、x=F-1(z)=Ft(),a=1,2,Nkak
12、aA(1-19)(1-20)大比如当M=20时,按式(1T6)计算得:P=90.25%,P=81.86%。lm2)拉丁超立方采样的稳健性好开云网址。假设一输出随机变量Y满足下式Y=cX(117)iii=1c是常数,丫是输入随机变量X的线性函数。在相同采样规模下,进行一定次数的蒙特卡罗ii模拟,每一次都能获得一个关于丫的分布情况.由每个丫的分布的期望值可以得到一个新的分布。用方差O表示这个分布的离散程度。若O越大,表明不同仿真间的差异越大,算法的稳健性越ZZ不好。文献指出通过拉丁超立方采样法得到的方差要比随机采样得到的方差小1/N2。表明一共进行总数为N3的随机采样得到的方差a与只需进行N次拉丁超立方采
13、样得到的方差a相同。ZZ1。3.3拉丁超立方采样步骤1)采样假设X,X,,X是随机潮流计算的N个输入变量。X的累积概率分布是:12NkZ=F(X),k=1,2,N(1-18)kkk取采样规模为A,采样步骤为:将Z的取值范围】0,1均匀分为A等份,即0,丄,丄,?,.,匕!,1;kAAAA从所有区间内依次抽取一个值作为一个采样值,区间内的抽取是随机的;c。由累积概率分布Zk的反函数变换后,便能得到输入变量Xk的样本数据.第a个区间z的采样值和X的第n个采样值如下:kka一randZ=,a=1,2,NkaA总共有N个输入变量,每个随机变量采样规模为A,假设将随机变量的数据以行为单位依次排列,那么最
14、终可以得到N*A阶的样本矩阵2)排序在求解随机潮流时,往往假设输入随机变量是独立的,但是按照上述方法得到的样本矩阵具有一定的相关性。我们需要分析和处理样本矩阵的关联性。使得变量数据值之间的关联性最小或者通过排序达到指定的相关系数。2系统模型建立光伏接入后的配电网系统主要由光伏发电系统、负荷和发电机三部分组成。太阳能光伏发电利用光伏电池可将光照转变为电动势的原理。在研究光伏并网后的随机潮流计算等有关问题时,首先要确定的是光伏发电的输出功率的随机特性,而此出力与太阳的光照强度密切相关,所以要想得到出力情况,必须先求出光照强度的随机分布30-3。4本次光伏发电,采用的是典型的Beta分布。此时我们可
15、以得到光照强度的概率密度函数为:r(a)r(p)Is、1max1S丿max(2-1)其中S是指光照强度统计时间内的实际值,s是指最大值.r是Gamma函数.a和p是形状max参数,将一段时间里太阳光照强度的期望值卩和方差5进行下式的变换便能得到形状参数35-36。卩(1-卩)1p=(1卩)卩(1-卩)(2-2)(2-3)假设光伏发电所用的电池方阵中有N个电池组,每个电池组的面积为A,光电转换效率为n耳n二1,2,N。n那么电池方阵总体的光电之间转化效率耳和方阵总的面积A分别是:A耳nn2-4)(2-5)门=-n=1AA亠nn=l此时这个电池方阵总的输出功率为:P二耳SA(26)N通过(24)-
16、(26),在光照强度的概率密度函数基础上,便能推导出光伏输出功率的概率密度函数为:f(P)=r(a+p)、a1了DAp1r(a)r(p)IP丿ma1P丿max(2-7)其中,P二耳AS,为光伏出力的最大值.mama当a=0.8,p=2时,光照强度的概率分布曲线为:概率密度函数(IIIIIP0,10.20J04CL50.6OJ0,8C191S/Smax图2。1形状参数为0。8和2时光照强度的概率分布图配电网中可以将接入光伏的节点视为PQ节点,主要由于通过调节电容器可以使得功率因数恒定.3IEEE30节点算例3.1IEEE30节点系统介绍IEEE-30节点系统包括6台发电机,30个节点与41条支路
22、用蒙特卡罗模拟法得到节点1电压的PDF与CDF如图3。1和3.2所示。可以看出两种算法还是存在差异的。(完整)电力系统潮流计算方法分析1(a)保留非线.0381.041.042U(b)牛顿拉夫逊图3.2两种算法下电压1的PDF图U1860000FDC-B(完整)电力系统潮流计算方法分析642000FDC-N1%25104TU(b)牛顿拉夫逊图3。3两种算法下电压1的CDF图3.3两种随机潮流算法的比较将以简单随机采样为基础的蒙特卡罗模拟法(MCSRS)和以拉丁超立方采样为基础的模拟法(MCLHS)得出的数据从准
23、确性和性能等方面做一个评估,全面比较两种随机潮流算法。3。3。1模型的准确性评估通过对输入随机变量的概率分布参数拟合,来分析所建立的模型的有效性和正确性。拟合的效果用相对误差指标来表示,表明分布情况的参数x的相对误差指标计算公式如下:x100%3-1)ccxfxbcxbc和c分别为参数x的样本拟合值和给定值.xfxb对光伏的输出功率采用Beta分布模型进行评估。Beta分布的两个形状参数的选取值为:a=0.9,卩二0.85.在一定规模下,根据光伏采样样本得到样本的平均值和方差,得到形状参数a,P的拟合值。并根据式(3-1)与实际的给定值0.9、0.85相比较得到误差.不同规模下分别采样50次后
25、的影响。由表可以看出,相同规模下,MCLHS比MCSRS的误差更小,用MCLHS生成的样本准确性更高。随着规模的增加,MCLHS和MCSRS生成的样本数据的正确性都有很大的提高。3。3.2性能评估通过算出的输出变量的平均值与标准差去评估MCLHS与MCSRS两种方法的计算精确度。计算公式如下:Euxcucuxfxbcuxbx100%Edxcdcdxfxbcdxbx100%32)33)上面两个式子式分别用来表示平均值与标准差的相对误差指标。采样规模为N时,一类输出变量便有N个数值,输出变量相对误差指标用这N个值的期望值表示.X分为mean、std、max和min四类。为减小随机性对结果产生的影响
26、,对两种方法在不同规模下分别采样50次,最后输出变量误差指标用50次误差的平均值mean表示,将这50次误差计算的标准差std、最大值max与最小值min用来评估上述方法收敛性与稳健性。cu和cd是误差计算的参考值.分别选取用xbxb20000次蒙特卡罗模拟得到的所选取的电压、功率和网损值来作为参考值。本次算例以节点18电压值、支路编号为3(3-4)的功率值与网损值作为研究对象.1)选取采样规模为500,以节点18电压值开云网址,支路3的功率值与网损值为研究对象,将得到的(完整)电力系统潮流计算方法分析平均值和标准差与参考值比较得到误差。两种方法均在此规模下进行50次仿线次计算结果的平均值、
29、CSRS的稳定性和收敛性较差一点。2)以支路3(P3-4)的有功功率为研究对象,在不同采样规模下,利用两种方法进行50次仿真,将得到的平均值、标准差和最大值绘制成图,如下:(完整)电力系统潮流计算方法分析由图可以看出,不同图模下MiCLiHS方计在不准确开云网址、稳健的性和收敛性都要优于mcsrs。且随着采较图样规模的增加,上述方法的计算准确性、稳键性与收敛性都有提高。个别出现误差增大的情况与蒙特卡罗模拟存在随机性有关,选取的20000次的参考值仍有一定的误差。3)以支路3(P34)的有功功率为研究对象。选取N=20000次的MCSRS与N=800次的MCLHS仿线的有功功率的分布情况.概率
30、分布函数图累积分布函数图3.6两种方法下累积分布函数对比图由图可以看出规模为800次的MCLHS基本可以达到规模为20000次的MCSRS的的水平。4总结(完整)电力系统潮流计算方法分析本文在讨论牛顿拉夫逊法和保留非线性法两种常规潮流算法的基础上,又研究了蒙特卡罗模拟和拉丁超立方采样两种随机潮流方法。牛拉法和保留非线性各有优缺点。比如牛拉法的迭代次数少,保留非线性的总计算时间少,这与书本上的知识也相符合.蒙特卡罗模拟通过采样化不确定为确定,而拉丁超立方采样通过改变采样方式,提高了样本质量,减小了采样规模。一个网络,可以使用的潮流计算的方法多种多样,针对不同的网络特点,针对不同的要求选择合适的计算方法。
1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。