发明名称 一种对电力系统微电网的全分布式二次调频方法
摘要 本发明涉及一种对电力系统微电网的全分布式二次调频方法,属于电力系统运行和控制技术领域。本发明方法首先设计了传统能源分布式电源和可再生能源分布式电源的发电成本函数,并通过求导得到其微增率。提出了分布式的数据采集和信息交互结构,采用基于次梯度的一致性算法,根据分布式电源采集的数据和交换的信息迭代计算得到输出功率的调整量,从而对分布式电源的输出功率进行实时控制。本发明在使得频率实现快速恢复的同时,保证二次调频过程中的发电成本最低,并在分布式电源之间合理分配功率。本发明的控制方法能够避免微电网的频率震荡,利用本发明方法,能够实现微电网全分布式的二次调频,并提高微电网运行的经济性。
申请公布号 CN104638643B 申请公布日期 2017.01.11
申请号 CN201510075150.3 申请日期 2015.02.12
申请人 清华大学 发明人 吴文传;张伯明;孙宏斌;王中冠;郭庆来;王彬
分类号 H02J3/00(2006.01)I;H02J3/46(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 北京清亦华知识产权代理事务所(普通合伙) 11201 代理人 罗文群
主权项 一种对电力系统微电网的全分布式二次调频方法,其特征在于该方法包括以下步骤:(1)设定微电网中包含有可再生能源分布式电源和传统能源分布式电源,设计一个微电网中可再生能源分布式电源的发电成本函数如下:<maths num="0001"><math><![CDATA[<mrow><msub><mi>C</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>P</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><msub><mi>P</mi><mi>i</mi></msub><mo>-</mo><msubsup><mi>P</mi><mi>i</mi><mi>max</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>/</mo><msubsup><mi>P</mi><mi>i</mi><mi>max</mi></msubsup><mo>=</mo><mfrac><mn>1</mn><msubsup><mi>P</mi><mi>i</mi><mi>max</mi></msubsup></mfrac><msubsup><mi>P</mi><mi>i</mi><mn>2</mn></msubsup><mo>-</mo><mn>2</mn><msub><mi>P</mi><mi>i</mi></msub><mo>+</mo><msubsup><mi>P</mi><mi>i</mi><mi>max</mi></msubsup></mrow>]]></math><img file="FDA0001014961610000011.GIF" wi="987" he="134" /></maths>其中P<sub>i</sub>是第i个可再生能源分布式电源输出的有功功率,P<sub>i</sub><sup>max</sup>是第i个可再生能源分布式电源预测最大发电容量,C<sub>i</sub>(P<sub>i</sub>)是第i个可再生能源分布式电源的发电成本;微电网中传统能源分布式电源的发电成本函数为:C<sub>j</sub>(P<sub>j</sub>)=a<sub>j</sub>P<sub>j</sub><sup>2</sup>+b<sub>j</sub>P<sub>j</sub>+c<sub>j</sub>其中,P<sub>j</sub>是第j个传统能源分布式电源输出的有功功率,a<sub>j</sub>、b<sub>j</sub>、c<sub>j</sub>均为发电成本函数系数,取值范围均大于零;分别对上述传统能源分布式电源的发电成本函数和可再生能源分布式电源的发电成本函数进行求导,得到第i个可再生能源分布式电源的发电成本微增率ICR<sub>i</sub>和第j个传统能源分布式电源的发电成本微增率ICR<sub>j</sub>:<maths num="0002"><math><![CDATA[<mrow><msub><mi>ICR</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>dC</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>P</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>dP</mi><mi>i</mi></msub></mrow></mfrac><mo>=</mo><mn>2</mn><mfrac><msub><mi>P</mi><mi>i</mi></msub><msubsup><mi>P</mi><mi>i</mi><mi>max</mi></msubsup></mfrac><mo>-</mo><mn>2</mn></mrow>]]></math><img file="FDA0001014961610000012.GIF" wi="539" he="133" /></maths><maths num="0003"><math><![CDATA[<mrow><msub><mi>ICR</mi><mi>j</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>dC</mi><mi>j</mi></msub><mrow><mo>(</mo><msub><mi>P</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>dP</mi><mi>j</mi></msub></mrow></mfrac><mo>=</mo><mn>2</mn><msub><mi>a</mi><mi>j</mi></msub><msub><mi>P</mi><mi>j</mi></msub><mo>+</mo><msub><mi>b</mi><mi>j</mi></msub></mrow>]]></math><img file="FDA0001014961610000013.GIF" wi="549" he="141" /></maths>(2)建立微电网中所有分布式电源之间的信息交换及迭代控制方式,实现对分布式电源的输出功率调整量的实时控制,具体步骤如下:(2‑1)微电网中各分布式电源分别采集各自的实时输出有功功率以及分布式电源与微电网接口处的频率值;(2‑2)所有分布式电源中任意相邻两个分布式电源之间交换各自的由上述步骤(1)得到的发电成本微增率;(2‑3)采用基于次梯度的一致性算法,计算微电网中各分布式电源的输出有功功率的调整量,具体过程如下:(2‑3‑1)将上述步骤(1)中得到的各分布式电源的发电成本微增率作为一致性变量,所有分布式电源输出的有功功率的控制目标函数为:<maths num="0004"><math><![CDATA[<mrow><mi>min</mi><mi> </mi><mi>f</mi><mrow><mo>(</mo><msub><mi>P</mi><mi>l</mi></msub><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mrow><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>P</mi><mi>l</mi></msub><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo><mo>-</mo><msub><mi>P</mi><mi>D</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow>]]></math><img file="FDA0001014961610000021.GIF" wi="655" he="127" /></maths>其中,P<sub>l</sub>[k]为微电网中第l个分布式电源在第k步迭代时输出的有功功率,P<sub>D</sub>是微电网的负荷总功率,将上述目标函数改写为如下的发电成本微增率的目标函数:<maths num="0005"><math><![CDATA[<mrow><mi>min</mi><mi> </mi><mi>f</mi><mrow><mo>(</mo><msub><mi>ICR</mi><mi>l</mi></msub><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mrow><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>r</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mfrac><mrow><msub><mi>ICR</mi><mi>r</mi></msub><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo><mo>-</mo><msub><mi>b</mi><mi>r</mi></msub></mrow><mrow><mn>2</mn><msub><mi>a</mi><mi>r</mi></msub></mrow></mfrac><mo>-</mo><msub><mi>P</mi><mi>D</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow>]]></math><img file="FDA0001014961610000022.GIF" wi="863" he="133" /></maths>其中,ICR<sub>l</sub>[k]为微电网中第l个分布式电源在第k步迭代时的发电成本微增率,a<sub>r</sub>、b<sub>r</sub>为微电网中第r个分布式电源发电成本函数系数;(2‑3‑2)根据上述步骤(2‑3‑1)发电成本微增率的目标函数,采用次梯度算法得到第k+1次发电成本微增率迭代计算式如下:<maths num="0006"><math><![CDATA[<mrow><msub><mi>ICR</mi><mi>l</mi></msub><mo>&lsqb;</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>&rsqb;</mo><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>r</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>&mu;</mi><mrow><mi>l</mi><mi>r</mi></mrow></msub><msub><mi>ICR</mi><mi>r</mi></msub><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo><mo>-</mo><msub><mi>d</mi><mi>l</mi></msub><mfrac><mrow><mo>&part;</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>ICR</mi><mi>l</mi></msub><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo><mo>)</mo></mrow></mrow><mrow><mo>&part;</mo><msub><mi>ICR</mi><mi>l</mi></msub><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></mfrac><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>r</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>&mu;</mi><mrow><mi>l</mi><mi>r</mi></mrow></msub><msub><mi>ICR</mi><mi>r</mi></msub><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo><mo>-</mo><msub><mi>d</mi><mi>l</mi></msub><mo>&CenterDot;</mo><mfrac><mn>1</mn><mrow><mn>2</mn><msub><mi>a</mi><mi>l</mi></msub></mrow></mfrac><mrow><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>r</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>P</mi><mi>r</mi></msub><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo><mo>-</mo><msub><mi>P</mi><mi>D</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001014961610000023.GIF" wi="1668" he="141" /></maths>其中,d<sub>l</sub>是中间变量,μ<sub>lr</sub>为第l个分布式电源与第r个分布式电源之间的通信系数,其数值由下式决定:<img file="FDA0001014961610000024.GIF" wi="618" he="431" />其中N<sub>l</sub>为所有与第l个分布式电源直接相连的分布式电源的集合,n<sub>l</sub>为与第l个分布式电源相连的分布式电源的个数,n<sub>r</sub>为与第r个分布式电源相连的分布式电源的个数;(2‑3‑3)采用频率偏差,替换第k+1次发电成本微增率迭代计算式中的功率不平衡量<img file="FDA0001014961610000025.GIF" wi="310" he="128" />得到的第k+1次发电成本微增率迭代计算式的表达式为:<maths num="0007"><math><![CDATA[<mrow><msub><mi>ICR</mi><mi>l</mi></msub><mo>&lsqb;</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>&rsqb;</mo><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>r</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>&mu;</mi><mrow><mi>l</mi><mi>r</mi></mrow></msub><msub><mi>ICR</mi><mi>r</mi></msub><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo><mo>-</mo><msub><mi>&lambda;</mi><mi>l</mi></msub><mrow><mo>(</mo><msubsup><mi>f</mi><mi>l</mi><mi>m</mi></msubsup><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo><mo>-</mo><msup><mi>f</mi><mo>*</mo></msup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001014961610000026.GIF" wi="939" he="143" /></maths>其中,λ<sub>l</sub>是迭代步长,f<sub>l</sub><sup>m</sup>[k]是采集的第l个分布式电源与微电网接口处的频率值,f<sup>*</sup>是微电网的频率额定值;(2‑3‑4)利用上述步骤(1)的发电成本微增率计算公式,根据(2‑3‑3)的第k+1次发电成本微增率迭代计算式,得到第j个传统能源分布式电源在第k+1步迭代计算时的输出功率调整量ΔP<sub>j</sub>[k+1]:<maths num="0008"><math><![CDATA[<mrow><msub><mi>&Delta;P</mi><mi>j</mi></msub><mo>&lsqb;</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>&rsqb;</mo><mo>=</mo><mfrac><mrow><msub><mi>ICR</mi><mi>j</mi></msub><mo>&lsqb;</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>&rsqb;</mo><mo>-</mo><msub><mi>b</mi><mi>j</mi></msub></mrow><mrow><mn>2</mn><msub><mi>a</mi><mi>j</mi></msub></mrow></mfrac><mo>-</mo><msubsup><mi>P</mi><mi>j</mi><mi>m</mi></msubsup><mo>&lsqb;</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0001014961610000031.GIF" wi="779" he="138" /></maths>其中,<img file="FDA0001014961610000032.GIF" wi="190" he="63" />是第j个传统能源分布式电源在第k+1次迭代时采集的输出有功功率值;(2‑3‑5)利用上述步骤(1)的发电成本微增率计算公式,根据(2‑3‑3)的第k+1次发电成本微增率迭代计算式,得到第i个可再生能源分布式电源在第k+1次迭代时输出功率的调整量ΔP<sub>i</sub>[k+1]:<maths num="0009"><math><![CDATA[<mrow><msub><mi>&Delta;P</mi><mi>i</mi></msub><mo>&lsqb;</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>&rsqb;</mo><mo>=</mo><msubsup><mi>P</mi><mi>i</mi><mi>max</mi></msubsup><mfrac><mrow><msub><mi>ICR</mi><mi>i</mi></msub><mo>&lsqb;</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>&rsqb;</mo><mo>+</mo><mn>2</mn></mrow><mn>2</mn></mfrac><mo>-</mo><msubsup><mi>P</mi><mi>i</mi><mi>m</mi></msubsup><mo>&lsqb;</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0001014961610000033.GIF" wi="821" he="119" /></maths>其中,<img file="FDA0001014961610000034.GIF" wi="189" he="55" />是第i个可再生能源分布式电源在第k+1次迭代时采集的输出有功功率值;(2‑4)根据上述计算得到的微电网中各分布式电源输出有功功率的调整量,实时控制分布式电源的输出功率调整量;(3)分别采集各分布式电源与微电网接口处的频率值,根据采集的频率值与微电网额定频率的偏差,实现对电力系统微电网的调频控制,包括以下步骤:(3‑1)分别采集微电网中各分布式电源与微电网接口处的频率值;(3‑2)分别计算上述频率值与微电网额定频率值之间的频率偏差,得到各分布式电源与微电网接口处的频率偏差值;(3‑3)设定一个频率偏差阈值ε,ε为正数,设定一个频率偏差维持时间阈值t,t为正数,将上述各分布式电源与微电网接口处的频率偏差值分别与频率偏差阈值进行比较,若频率偏差值小于ε,且频率偏差维持时间大于t,则进行步骤(3‑4),若频率偏差值小于ε,且频率偏差维持时间小于或等于t,则进行步骤(3‑5),若频率偏差值大于或等于ε,则进行步骤(3‑5);(3‑4)使上述步骤(2‑3‑3)中第k+1次发电成本微增率迭代计算式的表达式成为:<img file="FDA0001014961610000035.GIF" wi="565" he="126" />返回步骤(2);(3‑5)返回步骤(2)。
地址 100084 北京市海淀区清华园1号
您可能感兴趣的专利