CN105048499B - 基于模型预测控制的风电并网实时调度方法及系统 - Google Patents
基于模型预测控制的风电并网实时调度方法及系统 Download PDFInfo
- Publication number
- CN105048499B CN105048499B CN201510441034.9A CN201510441034A CN105048499B CN 105048499 B CN105048499 B CN 105048499B CN 201510441034 A CN201510441034 A CN 201510441034A CN 105048499 B CN105048499 B CN 105048499B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- power
- mtd
- wind
- 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
Links
Classifications
-
- 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
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A30/00—Adapting or protecting infrastructure or their operation
- Y02A30/60—Planning or developing urban green infrastructure
-
- 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
Abstract
本发明公开一种基于模型预测控制的风电并网实时调度方法及系统，能够提高风电并网实时调度的准确度，从而能够提升风电并网的稳定性。所述方法包括：计算大电网中各台风机的发电量指标；建立风机功率预测模型，以及风机出力计划与所述风机功率预测模型的功率预测结果的关系，通过对所述风机功率预测模型进行滚动优化，对风机出力计划与所对应的发电量指标之间的误差进行反馈校正，使风机出力计划与所对应的发电量指标之间的误差在预设的误差范围内。
Description
技术领域
本发明涉及电力系统运行与控制技术领域，具体涉及一种基于模型预测控制的风电并网实时调度方法及系统。
背景技术
新能源发电在不可再生资源日益枯竭的状况下发展迅速，其中风能的发展更为迅速，近几年风电场装机容量在逐渐增大，如何提高风电场有功功率预测的准确度和并网的稳定性以满足大电网调度要求是影响风电并网普及的难题。
风能由于不确定性和间歇性阻碍了对其的利用和推广。目前的风电场功率预测，大多数是在预测算法上进行改进，其中用到了机器学习和概率分布等对风电功率进行预测，选取预测精度高的模型。但是由于预测模型的固化以及在线数据的实时变化，使得单纯使用功率预测模型进行调度的准确度较低，从而导致风电并网的稳定性较差。
发明内容
本发明的目的在于，提供一种基于模型预测控制的风电并网实时调度方法及系统，能够提高风电并网实时调度的准确度，从而能够提升风电并网的稳定性。
为此目的，一方面，本发明提出一种基于模型预测控制的风电并网实时调度方法，包括：
计算大电网中各台风机的发电量指标；
建立风机功率预测模型，以及风机出力计划与所述风机功率预测模型的功率预测结果的关系，通过对所述风机功率预测模型进行滚动优化，对风机出力计划与所对应的发电量指标之间的误差进行反馈校正，使风机出力计划与所对应的发电量指标之间的误差在预设的误差范围内。
另一方面，本发明提出一种基于模型预测控制的风电并网实时调度系统，包括：
发电量指标计算单元，用于计算大电网中各台风机的发电量指标；
模型优化单元，用于建立风机功率预测模型，以及风机出力计划与所述风机功率预测模型的功率预测结果的关系，通过对所述风机功率预测模型进行滚动优化，对风机出力计划与所对应的发电量指标之间的误差进行反馈校正，使风机出力计划与所对应的发电量指标之间的误差在预设的误差范围内。
本发明实施例所述的基于模型预测控制的风电并网实时调度方法及系统，延续了对功率预测模型的算法改进，加上对其中误差、负荷调度的考虑，通过预测控制技术即功率预测模型的滚动优化和误差的反馈校正，使得功率预测较之现有的预测方法在精度上有所提高，相较于传统的风电并网仅利用预测模型进行调度，具有超前消除预测偏差、降低风电并网调度不确定性因素影响的优点，能够提高风电并网实时调度的准确度，从而能够提升风电并网的稳定性。
附图说明
图1为本发明基于模型预测控制的风电并网实时调度方法一实施例的流程示意图；
图2为本发明基于模型预测控制的风电并网实时调度方法另一实施例的流程示意图；
图3为本发明基于模型预测控制的风电并网实时调度系统一实施例的方框结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1所示，本实施例公开一种基于模型预测控制的风电并网实时调度方法，包括：
S1、计算大电网中各台风机的发电量指标；
S2、建立风机功率预测模型，以及风机出力计划与所述风机功率预测模型的功率预测结果的关系，通过对所述风机功率预测模型进行滚动优化，对风机出力计划与所对应的发电量指标之间的误差进行反馈校正，使风机出力计划与所对应的发电量指标之间的误差在预设的误差范围内。
本发明实施例所述的基于模型预测控制的风电并网实时调度方法，延续了对功率预测模型的算法改进，加上对其中误差、负荷调度的考虑，通过预测控制技术即功率预测模型的滚动优化和误差的反馈校正，使得功率预测较之现有的预测方法在精度上有所提高，相较于传统的风电并网仅利用预测模型进行调度，具有超前消除预测偏差、降低风电并网调度不确定性因素影响的优点，能够提高风电并网实时调度的准确度，从而能够提升风电并网的稳定性。
可选地，在本发明基于模型预测控制的风电并网实时调度方法的实施例中，所述计算大电网中各台风机的发电量指标，包括：
获取大电网的历史负荷数据，根据所述历史负荷数据建立负荷预测模型，对未来预设时段的负荷进行预测，得到所述未来预设时段大电网的负荷预测值，并确定出分配给风电场发电机组的发电量指标；
将所述风电场发电机组的发电量指标分配给风电场的各台风机，得到所述各台风机的发电量指标。
本发明实施例中，在得到未来预设时段大电网的负荷预测值之后，可以依照一定的比例给传统发电机组和风力发电机组分配发电量指标。在将风电场发电机组的发电量指标分配给风电场的各台风机时，可以依据当前各个风机的指标(比如风机所处环境的风速、湿度等)进行分配。
可选地，在本发明基于模型预测控制的风电并网实时调度方法的实施例中，所述建立风机功率预测模型，包括：
计算非线性风电功率预测拓扑在希尔伯特空间的线性系数α(*)，计算公式为
其中，l为待建模的风机的历史数据数量，i为2和l之间的整数，α＝(α1,α2,...,αl)，xi(i∈(1,2,...,l))为此风机第i个历史指标值，yi(i∈(1,2,...,l))为此风机对应第i个历史指标值的历史真实功率输出值，K(x,x)为映射核函数，ε为常数；
确定出风机功率预测模型为其中，或者x为当前预测时间段数值天气预报所预测的此风机的指标值，P为此风机的功率预测结果。
本发明实施例中，和α＝(α1,α2,...,αl)是非线性风电功率预测拓扑在希尔伯特空间的两组线性系数，当选择线性系数时，当选择线性系数α＝(α1,α2,...,αl)时，
可选地，在本发明基于模型预测控制的风电并网实时调度方法的实施例中，所述风机出力计划与所述风机功率预测模型的功率预测结果的关系为
其中，为第i个风机的出力计划，Pi为第i个风机的功率预测结果，kmax、kmin分别为风电场的机组调整裕度的上、下限的系数，Cwp是风电场机组的容量。
可选地，在本发明基于模型预测控制的风电并网实时调度方法的实施例中，所述对所述风机功率预测模型进行滚动优化，包括：
按照公式计算所述风机功率预测模型的功率预测结果的控制增量，计算公式为
其中，是风机功率预测模型的功率预测结果；N1、N2分别为优化时域的始值和终值，NU是控制时域，u为控制变量，Δ是差分算子，Δu是控制增量，yE是风机的期望输出功率，y是风机的实际输出功率，α和c为常数；
将所述控制增量和所述风机功率预测模型的功率预测结果进行加法运算，并将所述加法运算的结果作为所述风机功率预测模型的功率预测结果的修正值。
本发明利用控制增量对风机功率预测模型的功率预测结果进行修正，以实现对风机功率预测模型的滚动优化，保证了利用风机功率预测模型进行风机功率预测的准确度，有助于实现风电并网的稳定运行。
可选地，在本发明基于模型预测控制的风电并网实时调度方法的实施例中，所述对风机出力计划与所对应的发电量指标之间的误差进行反馈校正，包括：
计算风机出力计划与所述风机功率预测模型的功率预测结果之间的误差，根据计算出的所述误差计算出风机的指标调整量，并将风机的指标按照所述指标调整量进行相应的调整。
本发明实施例中，可以通过风机控制器对风机的指标进行调整，已知风电场出力计划，反馈校正风机控制器，从而达到满足整体发电计划的要求，通过对模型的滚动优化、风机控制器的反馈校正来进一步提高各个风电场风机输出功率预测的准确度并且满足大电网调度需求。
图2为本发明的一具体流程图，下面结合图2对本发明的具体实施例进行详细描述。
参看图2，本发明基于模型预测控制的风电并网实时调度方法，具体包括如下步骤：
(1)在当前时刻采样时刻t0，搜集某地历史负荷数据，建立超短期负荷预测模型，计算得到大电网的未来设定时间段的负荷输出；据此来分配传统发电机组和风机发电机组的发电指标；
(2)计算得到区域各个风电场各台风机的发电分配目标值；
(3)以风电场各台风机的历史指标数据作为训练集，构建风机的功率预测模型，将依据数值天气预报数据得到的各台风机的指标代入风机的功率预测模型，得到各台风机的功率预测值，在此过程中可配合风电场监测设备监测得到的风机的实际输出功率对各台风机的功率预测值进行优化；
(4)根据各台风机的功率预测值计算出各台风机的出力计划；
(5)计算各台风机的发电量指标和出力计划之间的误差，根据误差计算风机的指标调整量，通过对风机的指标进行调整(此调整过程可配合风电场监测设备进行指标监测)使各台风机的发电量指标和出力计划之间的误差尽量的小。
参看图3，本实施例公开一种基于模型预测控制的风电并网实时调度系统，包括：
发电量指标计算单元1，用于计算大电网中各台风机的发电量指标；
模型优化单元2，用于建立风机功率预测模型，以及风机出力计划与所述风机功率预测模型的功率预测结果的关系，通过对所述风机功率预测模型进行滚动优化，对风机出力计划与所对应的发电量指标之间的误差进行反馈校正，使风机出力计划与所对应的发电量指标之间的误差在预设的误差范围内。
本发明实施例所述的基于模型预测控制的风电并网实时调度系统，延续了对功率预测模型的算法改进，加上对其中误差、负荷调度的考虑，通过预测控制技术即功率预测模型的滚动优化和误差的反馈校正，使得功率预测较之现有的预测方法在精度上有所提高，相较于传统的风电并网仅利用预测模型进行调度，具有超前消除预测偏差、降低风电并网调度不确定性因素影响的优点，能够提高风电并网实时调度的准确度，从而能够提升风电并网的稳定性。
可选地，在本发明基于模型预测控制的风电并网实时调度系统的实施例中，所述发电量指标计算单元，包括：
计算子单元，用于获取大电网的历史负荷数据，根据所述历史负荷数据建立负荷预测模型，对未来预设时段的负荷进行预测，得到所述未来预设时段大电网的负荷预测值，并确定出分配给风电场发电机组的发电量指标；
分配子单元，用于将所述风电场发电机组的发电量指标分配给风电场的各台风机，得到所述各台风机的发电量指标。
可选地，在本发明基于模型预测控制的风电并网实时调度系统的实施例中，所述模型优化单元，包括：
线性系数计算子单元，用于计算非线性风电功率预测拓扑在希尔伯特空间的线性系数α(*)，计算公式为
其中，l为待建模的风机的历史数据数量，i为2和l之间的整数，α＝(α1，α2，...，αl)，xi(i∈(1,2,...,l))为此风机第i个历史指标值，yi(i∈(1，2，...，l))为此风机对应第i个历史指标值的历史真实功率输出值，K(x,x)为映射核函数，ε为常数；
确定出风机功率预测模型为其中，或者x为当前预测时间段数值天气预报所预测的此风机的指标值，p为此风机的功率预测结果。
可选地，在本发明基于模型预测控制的风电并网实时调度系统的实施例中，所述模型优化单元，用于计算风机出力计划与所述风机功率预测模型的功率预测结果之间的误差，根据计算出的所述误差计算出风机的指标调整量，并将风机的指标按照所述指标调整量进行相应的调整。
虽然结合附图描述了本发明的实施方式，但是本领域技术人员可以在不脱离本发明的精神和范围的情况下做出各种修改和变型，这样的修改和变型均落入由所附权利要求所限定的范围之内。
Claims (8)
1.一种基于模型预测控制的风电并网实时调度方法，其特征在于，包括：
计算大电网中各台风机的发电量指标；
建立风机功率预测模型，以及风机出力计划与所述风机功率预测模型的功率预测结果的关系，通过对所述风机功率预测模型进行滚动优化，对风机出力计划与所对应的发电量指标之间的误差进行反馈校正，使风机出力计划与所对应的发电量指标之间的误差在预设的误差范围内；
所述建立风机功率预测模型，包括：
计算非线性风电功率预测拓扑在希尔伯特空间的线性系数α(*)，计算公式为
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<munder>
<mi>min</mi>
<msup>
<mi>&alpha;</mi>
<mrow>
<mo>(</mo>
<mo>*</mo>
<mo>)</mo>
</mrow>
</msup>
</munder>
</mtd>
<mtd>
<mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>l</mi>
</munderover>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>&alpha;</mi>
<mi>i</mi>
<mo>*</mo>
</msubsup>
<mo>-</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>&alpha;</mi>
<mi>j</mi>
<mo>*</mo>
</msubsup>
<mo>-</mo>
<msub>
<mi>&alpha;</mi>
<mi>j</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mi>K</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&epsiv;</mi>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>l</mi>
</munderover>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>&alpha;</mi>
<mi>i</mi>
<mo>*</mo>
</msubsup>
<mo>+</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>-</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>l</mi>
</munderover>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>&alpha;</mi>
<mi>i</mi>
<mo>*</mo>
</msubsup>
<mo>-</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>l</mi>
</munderover>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msubsup>
<mi>&alpha;</mi>
<mi>i</mi>
<mo>*</mo>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msubsup>
<mi>&alpha;</mi>
<mi>i</mi>
<mo>*</mo>
</msubsup>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>l</mi>
<mo>.</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
其中，l为待建模的风机的历史数据数量，α＝(α1,α2,…,αl)，xi(i∈(1,2,…,l))为此风机第i个历史指标值，yi(i∈(1,2,…,l))为此风机对应第i个历史指标值的历史真实功率输出值，K(x,x)为映射核函数，ε为常数；
确定出风机功率预测模型为其中，或者x为当前预测时间段数值天气预报所预测的此风机的指标值，P为此风机的功率预测结果。
2.根据权利要求1所述的基于模型预测控制的风电并网实时调度方法，其特征在于，所述计算大电网中各台风机的发电量指标，包括：
获取大电网的历史负荷数据，根据所述历史负荷数据建立负荷预测模型，对未来预设时段的负荷进行预测，得到所述未来预设时段大电网的负荷预测值，并确定出分配给风电场发电机组的发电量指标；
将所述风电场发电机组的发电量指标分配给风电场的各台风机，得到所述各台风机的发电量指标。
3.根据权利要求1所述的基于模型预测控制的风电并网实时调度方法，其特征在于，所述风机出力计划与所述风机功率预测模型的功率预测结果的关系为
<mrow>
<msubsup>
<mi>P</mi>
<mi>i</mi>
<mrow>
<mi>w</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>k</mi>
<mi>min</mi>
</msub>
<msup>
<mi>C</mi>
<mrow>
<mi>w</mi>
<mi>p</mi>
</mrow>
</msup>
<mo>&le;</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<msup>
<mi>C</mi>
<mrow>
<mi>w</mi>
<mi>p</mi>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<msup>
<mi>C</mi>
<mrow>
<mi>w</mi>
<mi>p</mi>
</mrow>
</msup>
<mo>,</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo><</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<msup>
<mi>C</mi>
<mrow>
<mi>w</mi>
<mi>p</mi>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>k</mi>
<mi>max</mi>
</msub>
<msup>
<mi>C</mi>
<mrow>
<mi>w</mi>
<mi>p</mi>
</mrow>
</msup>
<mo>,</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>></mo>
<msub>
<mi>k</mi>
<mi>max</mi>
</msub>
<msup>
<mi>C</mi>
<mrow>
<mi>w</mi>
<mi>p</mi>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
其中，为第i个风机的出力计划，Pi为第i个风机的功率预测结果，kmax、kmin分别为风电场的机组调整裕度的上、下限的系数，Cwp是风电场机组的容量。
4.根据权利要求1所述的基于模型预测控制的风电并网实时调度方法，其特征在于，所述对所述风机功率预测模型进行滚动优化，包括：
按照公式计算所述风机功率预测模型的功率预测结果的控制增量，计算公式为
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>u</mi>
</mrow>
</munder>
<mi>E</mi>
<mo>{</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>m</mi>
<mo>=</mo>
<msub>
<mi>N</mi>
<mn>1</mn>
</msub>
</mrow>
<msub>
<mi>N</mi>
<mn>2</mn>
</msub>
</munderover>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mo>(</mo>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mi>m</mi>
</mrow>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>y</mi>
<mi>E</mi>
</msub>
<mo>(</mo>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mi>m</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>N</mi>
<mi>U</mi>
</mrow>
</munderover>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>&Delta;u</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mi>E</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>&alpha;</mi>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mi>c</mi>
<mo>,</mo>
<mi>t</mi>
<mo>></mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mi>E</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
其中，是风机功率预测模型的功率预测结果；N1、N2分别为优化时域的始值和终值，NU是控制时域，u为控制变量，Δ是差分算子，Δu是控制增量，yE是风机的期望输出功率，y是风机的实际输出功率，α和c为常数，λ(m)是指控制加权系数；
将所述控制增量和所述风机功率预测模型的功率预测结果进行加法运算，并将所述加法运算的结果作为所述风机功率预测模型的功率预测结果的修正值。
5.根据权利要求1所述的基于模型预测控制的风电并网实时调度方法，其特征在于，所述对风机出力计划与所对应的发电量指标之间的误差进行反馈校正，包括：
计算风机出力计划与所述风机功率预测模型的功率预测结果之间的误差，根据计算出的所述误差计算出风机的指标调整量，并将风机的指标按照所述指标调整量进行相应的调整。
6.一种基于模型预测控制的风电并网实时调度装置，其特征在于，包括：
发电量指标计算单元，用于计算大电网中各台风机的发电量指标；
模型优化单元，用于建立风机功率预测模型，以及风机出力计划与所述风机功率预测模型的功率预测结果的关系，通过对所述风机功率预测模型进行滚动优化，对风机出力计划与所对应的发电量指标之间的误差进行反馈校正，使风机出力计划与所对应的发电量指标之间的误差在预设的误差范围内；
所述模型优化单元，包括：
线性系数计算子单元，用于计算非线性风电功率预测拓扑在希尔伯特空间的线性系数α(*)，计算公式为
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<munder>
<mi>min</mi>
<msup>
<mi>&alpha;</mi>
<mrow>
<mo>(</mo>
<mo>*</mo>
<mo>)</mo>
</mrow>
</msup>
</munder>
</mtd>
<mtd>
<mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>l</mi>
</munderover>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>&alpha;</mi>
<mi>i</mi>
<mo>*</mo>
</msubsup>
<mo>-</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>&alpha;</mi>
<mi>j</mi>
<mo>*</mo>
</msubsup>
<mo>-</mo>
<msub>
<mi>&alpha;</mi>
<mi>j</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mi>K</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&epsiv;</mi>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>l</mi>
</munderover>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>&alpha;</mi>
<mi>i</mi>
<mo>*</mo>
</msubsup>
<mo>+</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>-</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>l</mi>
</munderover>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>&alpha;</mi>
<mi>i</mi>
<mo>*</mo>
</msubsup>
<mo>-</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>l</mi>
</munderover>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msubsup>
<mi>&alpha;</mi>
<mi>i</mi>
<mo>*</mo>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msubsup>
<mi>&alpha;</mi>
<mi>i</mi>
<mo>*</mo>
</msubsup>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>l</mi>
<mo>.</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
其中，l为待建模的风机的历史数据数量，α＝(α1,α2,…,αl)，xi(i∈(1,2,…,l))为此风机第i个历史指标值，yi(i∈(1,2,…,l))为此风机对应第i个历史指标值的历史真实功率输出值，K(x,x)为映射核函数，ε为常数；
确定出风机功率预测模型为其中，或者x为当前预测时间段数值天气预报所预测的此风机的指标值，P为此风机的功率预测结果。
7.根据权利要求6所述的基于模型预测控制的风电并网实时调度装置，其特征在于，所述发电量指标计算单元，包括：
计算子单元，用于获取大电网的历史负荷数据，根据所述历史负荷数据建立负荷预测模型，对未来预设时段的负荷进行预测，得到所述未来预设时段大电网的负荷预测值，并确定出分配给风电场发电机组的发电量指标；
分配子单元，用于将所述风电场发电机组的发电量指标分配给风电场的各台风机，得到所述各台风机的发电量指标。
8.根据权利要求6所述的基于模型预测控制的风电并网实时调度装置，其特征在于，所述模型优化单元，用于计算风机出力计划与所述风机功率预测模型的功率预测结果之间的误差，根据计算出的所述误差计算出风机的指标调整量，并将风机的指标按照所述指标调整量进行相应的调整。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510441034.9A CN105048499B (zh) | 2015-07-24 | 2015-07-24 | 基于模型预测控制的风电并网实时调度方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510441034.9A CN105048499B (zh) | 2015-07-24 | 2015-07-24 | 基于模型预测控制的风电并网实时调度方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105048499A CN105048499A (zh) | 2015-11-11 |
| CN105048499B true CN105048499B (zh) | 2017-11-03 |
Family
ID=54454821
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510441034.9A Active CN105048499B (zh) | 2015-07-24 | 2015-07-24 | 基于模型预测控制的风电并网实时调度方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105048499B (zh) |
Families Citing this family (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105375479B (zh) * | 2015-12-14 | 2017-11-17 | 东南大学 | 一种基于模型预测控制的分布式能源能量管理方法 |
| WO2018006920A1 (en) * | 2016-07-06 | 2018-01-11 | Vestas Wind Systems A/S | A wind power plant having a plurality of wind turbine generators and a power plant controller |
| CN106786807B (zh) * | 2016-12-15 | 2019-05-28 | 电子科技大学 | 一种基于模型预测控制的风电场有功功率控制方法 |
| CN107276221B (zh) * | 2017-05-23 | 2019-12-03 | 中国农业大学 | 一种优化风电消纳的电力系统调度方法 |
| CN107508327B (zh) * | 2017-10-20 | 2020-07-24 | 云南电网有限责任公司 | 基于模型预测控制的配电网自治区域并网功率控制方法 |
| CN107979112B (zh) * | 2017-11-30 | 2021-08-27 | 全球能源互联网研究院有限公司 | 一种风机控制方法、系统、终端及可读存储介质 |
| CN109376426B (zh) * | 2018-10-23 | 2023-04-07 | 株洲中车时代电气股份有限公司 | 一种风电并网功率调度方法及装置 |
| CN109657867B (zh) * | 2018-12-25 | 2023-04-07 | 广州汇电云联互联网科技有限公司 | 一种a类机组出力曲线分解方法 |
| CN109685271B (zh) * | 2018-12-25 | 2020-11-17 | 广州汇电云联互联网科技有限公司 | 一种统调负荷曲线分解方法 |
| CN112103939B (zh) * | 2019-06-18 | 2024-01-26 | 上海电机学院 | 一种电力系统稳定性分析方法及装置 |
| CN110474328A (zh) * | 2019-08-20 | 2019-11-19 | 杭州市电力设计院有限公司余杭分公司 | 一种新能源电站的并网控制方法及系统 |
| US10954919B1 (en) * | 2019-10-02 | 2021-03-23 | General Electric Company | Machine-learning model-based analytic for monitoring wind farm power performance |
| CN111342499B (zh) * | 2020-03-05 | 2023-09-08 | 宁夏嘉泽新能源股份有限公司 | 一种基于风功率预测数据的风电场实时调度方法 |
| CN112039132B (zh) * | 2020-08-14 | 2021-11-30 | 中国农业大学 | 一种不确定性电源的有功控制和精细指令分配方法 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102075014B (zh) * | 2011-01-06 | 2013-02-13 | 清华大学 | 消纳风电接入的大电网实时调度方法 |
| CN102170170B (zh) * | 2011-04-02 | 2013-02-13 | 清华大学 | 消纳风电接入的大电网调度滚动计划方法 |
| CN102738833B (zh) * | 2012-06-20 | 2014-07-09 | 湖北省电力公司 | 一种含风电的电力系统多时间尺度滚动协调调度方法 |
| CN103427444B (zh) * | 2013-05-16 | 2015-07-08 | 新疆大学 | 一种减小风电并网调度计划误差的控制方法 |
| CN104181895A (zh) * | 2014-08-28 | 2014-12-03 | 国家电网公司 | 适应新能源接入的短期与超短期协调滚动调度优化策略 |
-
2015
- 2015-07-24 CN CN201510441034.9A patent/CN105048499B/zh active Active
Non-Patent Citations (1)
| Title |
|---|
| 大规模风电接入电网的有功分层模型预测控制方法;张伯明等;《电力系统自动化》;20140510;第38卷(第9期);6-14 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105048499A (zh) | 2015-11-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105048499B (zh) | 基于模型预测控制的风电并网实时调度方法及系统 | |
| CN105046374B (zh) | 一种基于核极限学习机模型的功率区间预测方法 | |
| CN105375479B (zh) | 一种基于模型预测控制的分布式能源能量管理方法 | |
| CN103390116B (zh) | 采用分步方式的光伏电站发电功率预测方法 | |
| CN106451550B (zh) | 一种基于改进次梯度粒子群的微电网并网优化调度方法 | |
| CN102102626B (zh) | 风电场短期功率预测方法 | |
| CN103401236B (zh) | 一种基于风电场流动相关性的风电场机组分组方法 | |
| CN105811407B (zh) | 一种基于分布式牛顿法的微电网一次调频控制方法 | |
| CN107769254A (zh) | 一种风电集群轨迹预测与分层控制方法 | |
| CN102856917B (zh) | 一种配电网无功优化方法 | |
| CN103208803A (zh) | 风电和光电集中并网的无功电压优化控制方法 | |
| CN109523084A (zh) | 一种基于主元分析和机器学习的风电功率超短期预测方法 | |
| CN104978608A (zh) | 一种风电功率预测装置及预测方法 | |
| CN105207272B (zh) | 基于通用分布的电力系统动态随机经济调度方法及装置 | |
| CN105354620A (zh) | 一种风机发电功率的预测方法 | |
| CN105024398A (zh) | 一种基于最优风电置信度的优化调度方法 | |
| CN103489046A (zh) | 风电场短期功率预测方法 | |
| CN104500336B (zh) | 一种基于Hammerstein‑Wiener模型的风电机组恒功率广义预测控制方法 | |
| CN106953363A (zh) | 一种风电场限功率运行状态下电网旋转备用优化配置方法 | |
| CN107204615A (zh) | 一种实现功率预测的方法及系统 | |
| CN112510703B (zh) | 一种考虑煤耗曲线修正的多能源接入电网优化调度方法 | |
| CN103887813B (zh) | 基于风功率预测不确定度的风电系统运行的控制方法 | |
| CN105425591A (zh) | 基于风电场区域实时风速分析调控风电场有功的方法和装置 | |
| CN102509026A (zh) | 基于最大信息熵原理的风电场短期输出功率预测综合模型 | |
| CN107706932A (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 |