发明名称 一种电力系统潮流无解调整方法
摘要 本发明公开了一种电力系统潮流无解调整方法,涉及电力系统的运行与调度。该方法利用拉格朗日函数、牛顿法和对数障碍函数法,先将潮流无解问题转化为非线性规划问题,以当前状态为初始点,通过最优化理论确定迭代步长及最优搜索方向,再根据迭代步长和搜索方向,更新原对偶变量,确定下一次迭代的出发点,直到满足潮流收敛条件为止。本发明的优点是:所建立的非线性规划模型准确、灵活,目标函数可根据需要有多种选取方法,同时可以保证恢复后的潮流解的合理性;在迭代计算过程中不需要求解潮流方程式,也不依赖于潮流算法的收敛性。方程的迭代次数与系统规模关系不大,在大型电力系统中计算速度比现有手工方法有较大提高。
申请公布号 CN101976837A 申请公布日期 2011.02.16
申请号 CN201010508547.4 申请日期 2010.10.17
申请人 华中电网有限公司;华中科技大学 发明人 陈金富;奚江惠;洪峰;李勇;郭崇军;姜曼;徐友平
分类号 H02J3/00(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 武汉荆楚联合知识产权代理有限公司 42215 代理人 王健
主权项 1.一种电力系统潮流无解调整方法,其特征在于,该方法包含下列步骤:a)在给定的有m+n+1个节点的电力系统中,PQ节点为m个,PV节点为n个,平衡节点1个,其中第1到第<img file="2010105085474100001DEST_PATH_IMAGE001.GIF" wi="12" he="23" />节点为可切除的负荷节点,引入方程组:<img file="367247DEST_PATH_IMAGE002.GIF" wi="507" he="180" />(*)其中,<img file="2010105085474100001DEST_PATH_IMAGE003.GIF" wi="164" he="28" />表示第k个节点注入功率的切除比例,<img file="776231DEST_PATH_IMAGE004.GIF" wi="82" he="23" />为含参数<img file="2010105085474100001DEST_PATH_IMAGE005.GIF" wi="20" he="25" />的潮流方程,<img file="113059DEST_PATH_IMAGE006.GIF" wi="16" he="17" />为各个节点的电压幅值和相角,<img file="2010105085474100001DEST_PATH_IMAGE007.GIF" wi="72" he="28" />为系统状态变量,<img file="265692DEST_PATH_IMAGE008.GIF" wi="81" he="31" />为松弛变量,<img file="2010105085474100001DEST_PATH_IMAGE009.GIF" wi="66" he="28" />为<img file="571908DEST_PATH_IMAGE010.GIF" wi="36" he="25" />的拉格朗日乘子,<img file="2010105085474100001DEST_PATH_IMAGE011.GIF" wi="80" he="28" />为<img file="964844DEST_PATH_IMAGE012.GIF" wi="80" he="51" />的拉格朗日乘子,<img file="2010105085474100001DEST_PATH_IMAGE013.GIF" wi="134" he="50" />,<img file="173321DEST_PATH_IMAGE014.GIF" wi="128" he="46" />,<img file="2010105085474100001DEST_PATH_IMAGE015.GIF" wi="74" he="25" />,<img file="677115DEST_PATH_IMAGE016.GIF" wi="21" he="20" />为对数障碍参数,<img file="2010105085474100001DEST_PATH_IMAGE017.GIF" wi="17" he="21" />、<img file="837838DEST_PATH_IMAGE018.GIF" wi="20" he="21" />是分别以<img file="2010105085474100001DEST_PATH_IMAGE019.GIF" wi="22" he="26" />、<img file="588625DEST_PATH_IMAGE020.GIF" wi="26" he="25" />为元素的对角矩阵,<img file="2010105085474100001DEST_PATH_IMAGE021.GIF" wi="18" he="21" />、<img file="834799DEST_PATH_IMAGE022.GIF" wi="23" he="24" />是分别以<img file="2010105085474100001DEST_PATH_IMAGE023.GIF" wi="23" he="25" />、<img file="394480DEST_PATH_IMAGE024.GIF" wi="26" he="26" />为元素的对角矩阵;b)将方程组(*)中的<img file="2010105085474100001DEST_PATH_IMAGE025.GIF" wi="69" he="26" />,<img file="675289DEST_PATH_IMAGE026.GIF" wi="54" he="26" />进行赋值,<img file="2010105085474100001DEST_PATH_IMAGE027.GIF" wi="176" he="27" />;<img file="410027DEST_PATH_IMAGE028.GIF" wi="136" he="27" />,其中<img file="2010105085474100001DEST_PATH_IMAGE029.GIF" wi="73" he="28" />为原变量,<img file="143497DEST_PATH_IMAGE030.GIF" wi="60" he="28" />为对偶变量;c)对<img file="2010105085474100001DEST_PATH_IMAGE031.GIF" wi="42" he="26" />,<img file="238360DEST_PATH_IMAGE032.GIF" wi="43" he="26" />,<img file="2010105085474100001DEST_PATH_IMAGE033.GIF" wi="49" he="26" />,<img file="394183DEST_PATH_IMAGE034.GIF" wi="50" he="26" />,<img file="2010105085474100001DEST_PATH_IMAGE035.GIF" wi="48" he="28" />,<img file="299822DEST_PATH_IMAGE036.GIF" wi="50" he="28" />,<img file="2010105085474100001DEST_PATH_IMAGE037.GIF" wi="50" he="28" />,<img file="520588DEST_PATH_IMAGE038.GIF" wi="51" he="28" />,<img file="2010105085474100001DEST_PATH_IMAGE039.GIF" wi="46" he="28" />,<img file="419143DEST_PATH_IMAGE040.GIF" wi="48" he="28" />,<img file="DEST_PATH_IMAGE041.GIF" wi="49" he="27" />,<img file="674544DEST_PATH_IMAGE042.GIF" wi="51" he="27" />,<img file="2010105085474100001DEST_PATH_IMAGE043.GIF" wi="51" he="27" />,<img file="751084DEST_PATH_IMAGE044.GIF" wi="53" he="27" />,<img file="2010105085474100001DEST_PATH_IMAGE045.GIF" wi="34" he="24" />,<img file="196497DEST_PATH_IMAGE046.GIF" wi="35" he="24" />,<img file="2010105085474100001DEST_PATH_IMAGE047.GIF" wi="44" he="27" />,<img file="836425DEST_PATH_IMAGE048.GIF" wi="46" he="27" />,<img file="2010105085474100001DEST_PATH_IMAGE049.GIF" wi="45" he="25" />,<img file="8650DEST_PATH_IMAGE050.GIF" wi="45" he="25" />,<img file="2010105085474100001DEST_PATH_IMAGE051.GIF" wi="41" he="25" />进行赋值;令M = 1;d)计算方程组(*),求出各原变量和对偶变量的迭加量<img file="443042DEST_PATH_IMAGE052.GIF" wi="253" he="26" />;e)计算<img file="DEST_PATH_IMAGE053.GIF" wi="20" he="21" />,<img file="185870DEST_PATH_IMAGE054.GIF" wi="27" he="19" />和<img file="DEST_PATH_IMAGE055.GIF" wi="29" he="23" />,其中:<img file="360981DEST_PATH_IMAGE056.GIF" wi="102" he="45" /><img file="DEST_PATH_IMAGE057.GIF" wi="479" he="40" /><img file="325395DEST_PATH_IMAGE058.GIF" wi="460" he="39" />f)当下列不等式全部满足时,执行步骤i),否则,执行步骤g);<img file="DEST_PATH_IMAGE059.GIF" wi="380" he="123" /><img file="930688DEST_PATH_IMAGE060.GIF" wi="124" he="26" />,<img file="DEST_PATH_IMAGE061.GIF" wi="134" he="25" />,<img file="613342DEST_PATH_IMAGE062.GIF" wi="143" he="25" />g)对<img file="DEST_PATH_IMAGE063.GIF" wi="75" he="28" />;<img file="332424DEST_PATH_IMAGE064.GIF" wi="59" he="27" />重新进行赋值:<img file="DEST_PATH_IMAGE065.GIF" wi="340" he="28" />,<img file="416924DEST_PATH_IMAGE066.GIF" wi="273" he="26" />;h)令M= M+1,当M<<img file="255436DEST_PATH_IMAGE051.GIF" wi="41" he="25" />时,返回步骤d),当M≥<img file="972856DEST_PATH_IMAGE051.GIF" wi="41" he="25" />时,执行步骤j);i)计算结束,输出<img file="368065DEST_PATH_IMAGE025.GIF" wi="69" he="26" />和<img file="DEST_PATH_IMAGE067.GIF" wi="53" he="25" />的值,这些参数可以用于电网的运行调度;j)结束计算,用本方法无法得到原变量<img file="307071DEST_PATH_IMAGE029.GIF" wi="73" he="28" />和对偶变量<img file="313555DEST_PATH_IMAGE068.GIF" wi="56" he="27" />的合理值。
地址 430077 湖北省武汉市武昌区徐东大街47号