CN106845863A - 一种分布式风力发电出力与供暖负荷同步指数预测方法 - Google Patents
一种分布式风力发电出力与供暖负荷同步指数预测方法 Download PDFInfo
- Publication number
- CN106845863A CN106845863A CN201710098969.0A CN201710098969A CN106845863A CN 106845863 A CN106845863 A CN 106845863A CN 201710098969 A CN201710098969 A CN 201710098969A CN 106845863 A CN106845863 A CN 106845863A
- Authority
- CN
- China
- Prior art keywords
- vgn
- tgn
- ign
- pgn
- wind power
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000010438 heat treatment Methods 0.000 claims abstract description 63
- 238000010248 power generation Methods 0.000 claims abstract description 55
- 238000013528 artificial neural network Methods 0.000 claims abstract description 32
- 239000002245 particle Substances 0.000 claims description 59
- 230000001360 synchronised effect Effects 0.000 claims description 25
- 230000006870 function Effects 0.000 claims description 24
- 238000012549 training Methods 0.000 claims description 22
- 238000013178 mathematical model Methods 0.000 claims description 6
- 238000011946 reduction process Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 4
- 238000012544 monitoring process Methods 0.000 abstract description 4
- 230000009897 systematic effect Effects 0.000 abstract 1
- 230000008901 benefit Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005338 heat storage Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/043—Architecture, e.g. interconnection topology based on fuzzy logic, fuzzy membership or fuzzy inference, e.g. adaptive neuro-fuzzy inference systems [ANFIS]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- H02J3/386—
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/003—Load forecast, e.g. methods or systems for forecasting future load demand
-
- 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
-
- 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
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P80/00—Climate change mitigation technologies for sector-wide applications
- Y02P80/10—Efficient use of energy, e.g. using compressed air or pressurized fluid as energy carrier
- Y02P80/15—On-site combined power, heat or cool generation or distribution, e.g. combined heat and power [CHP] supply
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Health & Medical Sciences (AREA)
- Software Systems (AREA)
- General Health & Medical Sciences (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Quality & Reliability (AREA)
- Pure & Applied Mathematics (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Water Supply & Treatment (AREA)
- Power Engineering (AREA)
- Automation & Control Theory (AREA)
- Computational Mathematics (AREA)
- Fuzzy Systems (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Primary Health Care (AREA)
- Public Health (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
Abstract
本发明提供一种分布式风力发电出力与供暖负荷同步指数预测方法,涉及风力发电技术领域。该方法对风电供暖系统的系统参数的非线性时间序列进行相空间重构后,作为样本集进行分布式风力发电出力与供暖负荷同步指数数学建模,利用模糊神经网络求解该数学模型,针对相空间重构后的数据对下一时刻的系统同步性进行预测,得到分布式风力发电出力与供暖负荷同步指数预测值。本发明针对热电联产系统进行实时监测,测量系统运行参数以及系统的地理气象环境参数,对分布式风力发电出力与供暖负荷同步指数进行预测,根据计算结果实时地对热电联产系统进行控制,能够有效利用风能,显著提高热电联产系统运行的可靠性与经济性。
Description
技术领域
本发明涉及风力发电技术领域,尤其涉及一种分布式风力发电出力与供暖负荷同步指数预测方法。
背景技术
从能源产业发展的角度看,为解决“三北”地区供暖季弃风难题,通过电蓄热锅炉替代燃煤小锅炉,增加冬季低谷负荷时段的用电负荷,热电联产系统日益推广。以往分布式风力发电出力与供暖负荷同步性计算方法的特点是忽略网架结构与地理气象条件间的相互作用关系,由单一的气象条件对热电联产系统的同步性进行计算,不能有效利用电网运行状态数据及地理气象综合数据资源,计算准确度和可靠度不高。
发明内容
针对现有技术的缺陷,本发明提供一种分布式风力发电出力与供暖负荷同步指数预测方法,根据分布式风力发电与供暖联合系统的网架结构、地理气象条件、电网和负荷运行条件,对分布式风力发电出力与供暖负荷同步指数进行有效预测,对热电联产系统的风险及效益做出评估,使热电联产系统能安全、稳定、高效运行。
一种分布式风力发电出力与供暖负荷同步指数预测方法,包括如下步骤:
步骤1:建立分布式风力发电出力与供暖负荷同步预测系统参数非线性时间序列,具体方法为:
步骤1.1:根据时间顺序按相同的时间间隔测量风电供暖系统的系统参数,所述系统参数包括:分布式风电系统所在环境的风速v、温度T、气压P、供暖负荷总负荷电流I及系统母线电压V;
步骤1.2:根据测量的系统参数构建如下非线性时间序列:
其中,vgn1,vgn2,...,vgnn表示风速时间序列,Tgn1,Tgn2,....,Tgnn表示温度时间序列,Pgn1,Pgn2,....,Pgnn表示气压时间序列,Ign1,Ign2,....,Ignn表示总负荷电流时间序列,Vgn1,Vgn2,...,Vgnn表示系统母线电压平均值时间序列;
步骤2:采用坐标延迟法对构建的非线性时间序列进行相空间事构,具体如下:
根据风速时间序列vgn1,vgn2,....,vgnn,在n维状态空间中重构风速状态矢量vgni′为:vgni′={vgni,vgni+τ,...,vgni+(m-1)τ},其中,i=1,2,...,n,τ为延迟时间,m为嵌入维数;
根据温度时间序列Tgn1,Tgn2,....,Tgnn,在n维状态空间中重构温度状态矢量Tgni′为:Tgni′={Tgni,Tgni+τ,…,Tgni+(m-1)τ};
根据气压时间序列Pgn1,Pgn2,....,Pgnn,在n维状态空间中重构气压状态矢量Pgni′为:Pgni′={Pgni,Pgni+τ,…,Pgni+(m-1)τ};
根据供暖负荷总负荷电流的时间序列为Ign1,Ign2,....,Ignn,在n维状态空间中重构总负荷电流状态矢量Igni′为:Igni′={Igni,Igni+τ,...,Igni+(m-1)τ};
根据母线电压的时间序列为Vgn1,Vgn2,...,Vgnn,在n维状态空间中重构母线电压状态矢量Vgni′为:Vgni′={Vgni,Vgni+τ,...,Vgni+(m-1)τ};
步骤3:将步骤2重构的相空间中的相点作为样本集,建立分布式风力发电出力与供暖负荷同步指数数学模型,如下式所示;
其中,ygn为分布式风力发电出力与供暖负荷同步指数,fmb(vgni′,Tgni′,Pgni′,Igni′,Vgni′)为目标函数,其中gni′为公式简化过程的中间符号,gcf(vgni′,Tgni′,Pgni′)为分布式风力发电出力与供暖负荷同步指数目标函数的惩罚因子,rys(Igni′,Vgni′)为分布式风力发电出力与供暖负荷同步指数目标函数的约束项,
步骤4:利用模糊神经网络求解分布式风力发电出力与供暖负荷同步指数数学模型,针对相空间重构后的数据对下一时刻的系统同步性进行预测,得到分布式风力发电出力与供暖负荷同步指数预测值,具体方法为:
步骤4.1:采用粒子群聚类算法对模糊神经网络隐含层节点数进行求取;
利用粒子群聚类算法对重构后的相空间相点所代表的系统运行工况进行分类,得到分类数mgn,并将mgn作为神经网络隐含层的节点数,其中粒子群聚类算法的输入数据为相空间重构后的相点,输出为分类数mgn,输入数据个数为Ngn,最大迭代次数为τmax;
步骤4.2:根据隐含层节点数mgn建立多输入单输出的模糊神经网络系统,根据模糊神经网络系统计算分布式风力发电出力与供暖负荷同步指数预测值。
进一步地,所述步骤4.1的具体方法为:
步骤4.1.1:分类数初始化,mgn=1;
步骤4.1.2:对粒子群参数进行初始化;相空间相点中随机选择mgn个值作为聚类中心值,并将这mgn个聚类中心值作为粒子群初始值;
计算粒子的适应度,粒子群适应度用如下函数计算;
其中,ffit为粒子群适应度函数,Jm为各个粒子到聚类中心的欧氏距离之和;
步骤4.1.3:进行粒子群迭代搜索,得到运行方式分类数mgl下的最优聚类中心,最大迭代次数为τmax,粒子群迭代更新方程如下式所示;
其中,为第k次迭代下,第i个粒子的速度;表示第k次迭代为止第i个粒子的最佳位置;表示第k次迭代为止粒子种群中所有粒子位置中的最佳位置,表示第i个粒子第k次迭代所在位置,c1、c2表示学习常数,rand()为0~1之间的随机数,w为惯性权重,用于权衡局部最优能力和全局最优能力;
粒子群适应度用如下函数计算;
其中,ffit为粒子群适应度函数,Jm为各个粒子到聚类中心的欧氏距离之和;
步骤4.1.4:计算聚类中心值间的欧氏距离之和JPphm,其中phm=1,2,…,Ngn;
步骤4.1.5:将mgn加1,若mgn<Ngn,则重复步骤4.1.2至4.1.5,否则,执行步骤4.1.6;
步骤4.1.6:比较分类数mgn从1到Ngn对应的欧氏距离之和将JPphm最小值所对应的mgn作为神经网络隐含层的节点数。
进一步地,所述步骤4.2的具体方法为:
步骤4.2.1:神经网络目标函数如下式所示;
其中,为神经网络第m次训练后系统期望输出,ygn(m)为神经网络第m次训练后系统实际输出;
步骤4.2.2:进行网络权值训练,计算输出层权值ω,权值更新公式如下式所示;
ω(m)=ω(m-1)+△ω(m)+a(ω(m-1)-ω(m-2))
其中,ω(m)为第m次训练后的输出层权值,η为学习速率,a为动量因子,η,a∈[0,1];设定阈值当时,网络权值训练结束;
步骤4.2.3:预测分布式风力发电出力与供暖负荷同步指数;
将重构相空间中的相点vgni′,Tgni′,Pgni′,Igni′和Vgni′输入权值训练结束的神经网络,网络实际输出的ygn即为分布式风力发电出力与供暖负荷同步指数的预测值。
由上述技术方案可知,本发明的有益效果在于:本发明提供的一种分布式风力发电出力与供暖负荷同步指数预测方法,针对热电联产系统进行实时监测,选取系统运行参数——供暖负荷总负荷电流、系统母线电压平均值,以及系统所在环境地理气象环境参数——风速、温度、气压;并根据监测参数对分布式风力发电出力与供暖负荷同步指数进行预测,根据计算结果实时地对热电联产系统进行控制,能够有效利用风能,显著提高热电联产系统运行的可靠性与经济性。
附图说明
图1为本发明实施例提供的分布式风力发电出力与供暖负荷同步指数预测方法流程图。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
一种分布式风力发电出力与供暖负荷同步指数预测方法,如图1所示,具体步骤如下。
步骤1:建立分布式风力发电出力与供暖负荷同步预测系统参数非线性时间序列;
步骤1.1:根据时间顺序按相同的时间间隔测量风电供暖系统的一系列系统参数,其中系统参数包括:分布式风电系统所在环境的风速v、温度T、气压P、供暖负荷总负荷电流I及系统母线电压V;
步骤1.2:根据测量的系统参数构建如下非线性时间序列:
其中,vgn1,vgn2,....,vgnn表示风速时间序列、Tgn1,Tgn2,....,Tgnn表示温度时间序列、Pgn1,Pgn2,....,Pgnn表示气压时间序列、Ign1,Ign2,....,Ignn表示总负荷电流时间序列、Vgn1,Vgn2,...,Vgnn表示系统母线电压平均值时间序列。
步骤2:采用坐标延迟法对构建的非线性时间序列进行相空间重构,具体如下:
在n维状态空间中重构的风速状态矢量vgni′、温度状态矢量Tgni′、气压状态矢量Pgni′、总负荷电流状态矢量Igni′、母线电压状态矢量Vgni′分别表示为:
vgni′={vgni,vgni+τ,...,vgηi+(m-1)τ};
Tgni′={Tgni,Tgni+τ,...,Tgni+(m-1)τ};
Pgni′={Pgni,Pgni+τ,…,Pgni+(m-1)τ};
Igni′={Igni,Igni+τ,...,Igni+(m-1)τ};
Vgni′={Vgni,Vgni+τ,...,Vgni+(m-1)τ};
其中,i=1,2,...,n,τ为延迟时间,m为嵌入维数。在本实施例中,选择m=10,τ=105ms。
步骤3:将步骤2重构的相空间中的相点作为样本集,建立分布式风力发电出力与供暖负荷同步指数数学模型,如下式所示;
其中,ygn为分布式风力发电出力与供暖负荷同步指数,fmb(vgni′,Tgni′,Pgni′,Igni′,Vgni′)为目标函数,其中gni′为公式简化过程的中间符号,gcf(vgni′,Tgni′,Pgni′)为分布式风力发电出力与供暖负荷同步指数目标函数的惩罚因子,rys(Igni′,Vgni′)为分布式风力发电出力与供暖负荷同步指数目标函数的约束项,max()为相空间重构后数据中的最大值,min()为相空间重构后数据中的最小值。
步骤4:利用模糊神经网络求解分布式风力发电出力与供暖负荷同步指数数学模型,针对相空间重构后的数据对下一时刻的系统同步性进行预测,得到分布式风力发电出力与供暖负荷同步指数预测值,分布式风力发电出力与供暖负荷同步指数越大则分布式风力发电出力与供暖负荷同步可能性越大,具体方法为:
步骤4.1:采用粒子群聚类算法对模糊神经网络隐含层节点数进行求取;
利用粒子群聚类算法对重构后的相空间相点所代表的系统运行工况进行分类,得到分类数mgn,并将mgn作为神经网络隐含层的节点数,其中粒子群聚类算法的输入数据为相空间重构后的相点,输出为分类数mgn,输入数据个数Ngn=500,具体方法如下:
步骤4.1.1:分类数初始化,mgn=1;
步骤4.1.2:对粒子群参数进行初始化;相空间相点中随机选择mgn个值作为聚类中心值,并将这mgn个聚类中心值作为粒子群初始值;
计算粒子的适应度,粒子群适应度用如下函数计算;
其中,ffit为粒子群适应度函数,Jm为各个粒子到聚类中心的欧氏距离之和;
步骤4.1.3:进行粒子群迭代搜索,得到运行方式分类数mgl下的最优聚类中心,最大迭代次数为τmax=50,粒子群迭代更新方程如下式所示;
其中,为第k次迭代下,第i个粒子的速度;表示第k次迭代为止第i个粒子的最佳位置;表示第k次迭代为止粒子种群中所有粒子位置中的最佳位置,表示第i个粒子第k次迭代所在位置,c1、c2表示学习常数,rand()为0~1之间的随机数,w为惯性权重,用于权衡局部最优能力和全局最优能力。本实施例中,c1=1,c2=2,w=0.6,当随机数rand()分别为0.2或0.5时,粒子群算法寻优效果最好。
粒子群适应度用如下函数计算;
其中,ffit为粒子群适应度函数,Jm为各个粒子到聚类中心的欧氏距离之和。
步骤4.1.4:计算聚类中心值间的欧氏距离之和JPphm,其中phm=1,2,…,Ngn;
步骤4.1.5:将mgn加1,若mgn<Ngn,则重复步骤4.1.2至4.1.5,否则,执行步骤4.1.6;
步骤4.1.6:比较分类数mgn从1到Ngn对应的欧氏距离之和将JPphm最小值所对应的mgn作为神经网络隐含层的节点数。
本实施例中,设最大迭代次数为500,得到的神经网络隐含层的节点数为11。
步骤4.2:根据隐含层节点数mgn建立多输入单输出的模糊神经网络系统,根据模糊神经网络系统计算分布式风力发电出力与供暖负荷同步指数预测值,具体方法为:
步骤4.2.1:神经网络目标函数如下式所示;
其中,为神经网络第m次训练后系统期望输出,ygn(m)为神经网络第m次训练后系统实际输出;
步骤4.2.2:进行网络权值训练,计算输出层权值ω,权值更新公式如下式所示;
ω(m)=ω(m-1)+Δω(m)+a(ω(m-1)-ω(m-2))
其中,ω(m)为第m次训练后的输出层权值,η为学习速率,a为动量因子,η,a∈[0,1];设定阈值当时,网络权值训练结束;本实施例中,设定阈值为训练结束后fa=0.048。
步骤4.2.3:预测分布式风力发电出力与供暖负荷同步指数;
将重构相空间的相点vgni′,Tgni′,Pgni′,Igni′和Vgni′输入权值训练结束的神经网络,网络实际输出的ygn即为分布式风力发电出力与供暖负荷同步指数的预测值。本实施例中,最终得到的分布式风力发电出力与供暖负荷同步指数预测值为ygn=0.8。
在分布式风力发电系统的控制中,根据预测得到的同步指数来进行风机出力大小的控制,若得到的预测值ygn小于0.7,则增加风机出力,若得到的预测值ygn处于0.7~1.2时,则保持风机现有出力,若得到的预测值ygn大于1.2时,则减小风机出力,使分布式风力发电出力与供暖负荷同步指数保持在0.7~1.2之间,从而提高系统的稳定性。
本实施例提供的一种分布式风力发电出力与供暖负荷同步指数预测方法,根据分布式风力发电与供暖联合系统的网架结构、地理气象条件、电网和负荷运行条件,针对热电联产系统进行实时监测,选取系统运行参数——供暖负荷总负荷电流、系统母线电压平均值,以及系统所在环境地理气象环境参数——风速、温度、气压,并根据监测参数对分布式风力发电出力与供暖负荷同步指数进行预测,根据计算结果实时地对热电联产系统进行控制,对热电联产系统的风险及效益做出评估,能够有效利用风能,使热电联产系统能安全、稳定、高效运行。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明权利要求所限定的范围。
Claims (3)
1.一种分布式风力发电出力与供暖负荷同步指数预测方法,其特征在于:该方法包括如下步骤:
步骤1:建立分布式风力发电出力与供暖负荷同步预测系统参数非线性时间序列,具体方法为:
步骤1.1:根据时间顺序按相同的时间间隔测量风电供暖系统的系统参数,所述系统参数包括:分布式风电系统所在环境的风速v、温度T、气压P、供暖负荷总负荷电流I及系统母线电压V;
步骤1.2:根据测量的系统参数构建如下非线性时间序列:
其中,vgn1,vgn2,....,vgnn表示风速时间序列,Tgn1,Tgn2,....,Tgnn表示温度时间序列,Pgn1,Pgn2,....,Pgnn表示气压时间序列,Ign1,Ign2,....,Ignn表示总负荷电流时间序列,Vgn1,Vgn2,...,Vgnn表示系统母线电压平均值时间序列;
步骤2:采用坐标延迟法对构建的非线性时间序列进行相空间重构,具体如下:
根据风速时间序列vgn1,vgn2,...,vgnn,在n维状态空间中重构风速状态矢量vgni′为:vgni′={vgni,vgni+τ,...,vgni+(m-1)τ},其中,i=1,2,...,n,τ为延迟时间,m为嵌入维数;
根据温度时间序列Tgn1,Tgn2,....,Tgnn,在n维状态空间中重构温度状态矢量Tgni′为:Tgni′={Tgni,Tgni+τ,...,Tgni+(m-1)τ};
根据气压时间序列Pgn1,Pgn2,....,Pgnn,在n维状态空间中重构气压状态矢量Pgni′为:Pgni′={Pgni,Pgni+τ,…,Pgni+(m-1)τ};
根据供暖负荷总负荷电流的时间序列为Ign1,Ign2,....,Ignn,在n维状态空间中重构总负荷电流状态矢量Igni′为:Igni′={Igni,Igni+τ,...,Igni+(m-1)τ};
根据母线电压的时间序列为Vgn1,Vgn2,...,Vgnn,在n维状态空间中重构母线电压状态矢量Vgni′为:Vgni′={Vgni,Vgni+τ,...,Vgni+(m-1)τ};
步骤3:将步骤2重构的相空间中的相点作为样本集,建立分布式风力发电出力与供暖负荷同步指数数学模型,如下式所示;
其中,ygn为分布式风力发电出力与供暖负荷同步指数,fmb(vgni′,Tgni′,Tgni′,Igni′,Vgni′)为目标函数,其中gni′为公式简化过程的中间符号,gcf(vgni′,Tgni′,Pgni′)为分布式风力发电出力与供暖负荷同步指数目标函数的惩罚因子,rys(Igni′,Vgni′)为分布式风力发电出力与供暖负荷同步指数目标函数的约束项,
步骤4:利用模糊神经网络求解分布式风力发电出力与供暖负荷同步指数数学模型,针对相空间重构后的数据对下一时刻的系统同步性进行预测,得到分布式风力发电出力与供暖负荷同步指数预测值,具体方法为:
步骤4.1:采用粒子群聚类算法对模糊神经网络隐含层节点数进行求取;
利用粒子群聚类算法对重构后的相空间相点所代表的系统运行工况进行分类,得到分类数mgn,并将mgn作为神经网络隐含层的节点数,其中粒子群聚类算法的输入数据为相空间重构后的相点,输出为分类数mgn,输入数据个数为Ngn,最大迭代次数为τmax;
步骤4.2:根据隐含层节点数mgn建立多输入单输出的模糊神经网络系统,根据模糊神经网络系统计算分布式风力发电出力与供暖负荷同步指数预测值。
2.根据权利要求1所述的一种分布式风力发电出力与供暖负荷同步指数预测方法,其特征在于:所述步骤4.1的具体方法为:
步骤4.1.1:分类数初始化,mgn=1;
步骤4.1.2:对粒子群参数进行初始化;相空间相点中随机选择mgn个值作为聚类中心值,并将这mgn个聚类中心值作为粒子群初始值;
计算粒子的适应度,粒子群适应度用如下函数计算;
其中,ffit为粒子群适应度函数,Jm为各个粒子到聚类中心的欧氏距离之和;
步骤4.1.3:进行粒子群迭代搜索,得到运行方式分类数mgl下的最优聚类中心,最大迭代次数为τmax,粒子群迭代更新方程如下式所示;
其中,为第k次迭代下,第i个粒子的速度;表示第k次迭代为止第i个粒子的最佳位置;表示第k次迭代为止粒子种群中所有粒子位置中的最佳位置,表示第i个粒子第k次迭代所在位置,c1、c2表示学习常数,rand()为0~1之间的随机数,w为惯性权重,用于权衡局部最优能力和全局最优能力;
粒子群适应度用如下函数计算;
其中,ffit为粒子群适应度函数,Jm为各个粒子到聚类中心的欧氏距离之和;
步骤4.1.4:计算聚类中心值间的欧氏距离之和JPphm,其中phm=1,2,…,Ngn;
步骤4.1.5:将mgn加1,若mgn<Ngn,则重复步骤4.1.2至4.1.5,否则,执行步骤4.1.6;
步骤4.1.6:比较分类数mgn从1到Ngn对应的欧氏距离之和将JPphm最小值所对应的mgn作为神经网络隐含层的节点数。
3.根据权利要求2所述的一种分布式风力发电出力与供暖负荷同步指数预测方法,其特征在于:所述步骤4.2的具体方法为:
步骤4.2.1:神经网络目标函数如下式所示;
其中,为神经网络第m次训练后系统期望输出,ygn(m)为神经网络第m次训练后系统实际输出;
步骤4.2.2:进行网络权值训练,计算输出层权值ω,权值更新公式如下式所示;
ω(m)=ω(m-1)+△ω(m)+a(ω(m-1)-ω(m-2))
其中,ω(m)为第m次训练后的输出层权值,η为学习速率,a为动量因子,η,a∈[0,1];设定阈值当时,网络权值训练结束;
步骤4.2.3:预测分布式风力发电出力与供暖负荷同步指数;
将重构相空间中的相点vgni′,Tgni′,Pgni′,Igni′和Vgni′输入权值训练结束的神经网络,网络实际输出的ygn即为分布式风力发电出力与供暖负荷同步指数的预测值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710098969.0A CN106845863A (zh) | 2017-02-23 | 2017-02-23 | 一种分布式风力发电出力与供暖负荷同步指数预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710098969.0A CN106845863A (zh) | 2017-02-23 | 2017-02-23 | 一种分布式风力发电出力与供暖负荷同步指数预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106845863A true CN106845863A (zh) | 2017-06-13 |
Family
ID=59133163
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710098969.0A Pending CN106845863A (zh) | 2017-02-23 | 2017-02-23 | 一种分布式风力发电出力与供暖负荷同步指数预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106845863A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107575928A (zh) * | 2017-09-27 | 2018-01-12 | 赫普科技发展(北京)有限公司 | 一种农村分布式风电供暖系统及方法 |
CN108376262A (zh) * | 2018-02-23 | 2018-08-07 | 新疆大学 | 一种风电出力典型特性的分析模型构建方法 |
CN109086805A (zh) * | 2018-07-12 | 2018-12-25 | 华南理工大学 | 一种基于深度神经网络和成对约束的聚类方法 |
CN109539359A (zh) * | 2018-10-15 | 2019-03-29 | 国网辽宁省电力有限公司 | 分工况pid+自适应前馈补偿的相变电蓄热供暖系统及方法 |
CN111798034A (zh) * | 2020-06-10 | 2020-10-20 | 河海大学 | 基于风流场时空图像学习的风电功率预测系统及预测方法 |
CN113404655A (zh) * | 2020-08-24 | 2021-09-17 | 湖南科技大学 | 一种基于ps0-anfis的风力发电机传感器状态诊断系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102082433A (zh) * | 2010-11-29 | 2011-06-01 | 沈阳工业大学 | 一种风电并网系统电压稳定预测装置及方法 |
US20120035898A1 (en) * | 2005-05-28 | 2012-02-09 | Carlos Repelli | Method and apparatus for generating an environmental element prediction for a point of interest |
CN102915511A (zh) * | 2012-09-21 | 2013-02-06 | 四川大学 | 一种电力负荷的混沌相空间神经网络模型安全监测方法 |
CN103301926A (zh) * | 2013-05-22 | 2013-09-18 | 西安交通大学 | 火电厂球磨机制粉系统运行工况划分方法 |
CN104361393A (zh) * | 2014-09-06 | 2015-02-18 | 华北电力大学 | 基于粒子群优化算法改进的神经网络模型用于数据预测方法 |
CN105354635A (zh) * | 2015-10-30 | 2016-02-24 | 国家电网公司 | 一种配电网不平衡谐波负荷指数预测方法 |
CN106355299A (zh) * | 2016-10-31 | 2017-01-25 | 国家电网公司 | 一种分布式光伏并网节点薄弱性指数预测方法 |
CN106374531A (zh) * | 2016-10-31 | 2017-02-01 | 国家电网公司 | 一种主动配电网光储能量交换指数评估方法 |
-
2017
- 2017-02-23 CN CN201710098969.0A patent/CN106845863A/zh active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120035898A1 (en) * | 2005-05-28 | 2012-02-09 | Carlos Repelli | Method and apparatus for generating an environmental element prediction for a point of interest |
CN102082433A (zh) * | 2010-11-29 | 2011-06-01 | 沈阳工业大学 | 一种风电并网系统电压稳定预测装置及方法 |
CN102915511A (zh) * | 2012-09-21 | 2013-02-06 | 四川大学 | 一种电力负荷的混沌相空间神经网络模型安全监测方法 |
CN103301926A (zh) * | 2013-05-22 | 2013-09-18 | 西安交通大学 | 火电厂球磨机制粉系统运行工况划分方法 |
CN104361393A (zh) * | 2014-09-06 | 2015-02-18 | 华北电力大学 | 基于粒子群优化算法改进的神经网络模型用于数据预测方法 |
CN105354635A (zh) * | 2015-10-30 | 2016-02-24 | 国家电网公司 | 一种配电网不平衡谐波负荷指数预测方法 |
CN106355299A (zh) * | 2016-10-31 | 2017-01-25 | 国家电网公司 | 一种分布式光伏并网节点薄弱性指数预测方法 |
CN106374531A (zh) * | 2016-10-31 | 2017-02-01 | 国家电网公司 | 一种主动配电网光储能量交换指数评估方法 |
Non-Patent Citations (1)
Title |
---|
谷彦章: "风电供热综合效益评价研究", 《中国优秀硕士学位论文全文数据库经济与管理科学辑(月刊)》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107575928A (zh) * | 2017-09-27 | 2018-01-12 | 赫普科技发展(北京)有限公司 | 一种农村分布式风电供暖系统及方法 |
CN108376262A (zh) * | 2018-02-23 | 2018-08-07 | 新疆大学 | 一种风电出力典型特性的分析模型构建方法 |
CN108376262B (zh) * | 2018-02-23 | 2021-08-10 | 新疆大学 | 一种风电出力典型特性的分析模型构建方法 |
CN109086805A (zh) * | 2018-07-12 | 2018-12-25 | 华南理工大学 | 一种基于深度神经网络和成对约束的聚类方法 |
CN109086805B (zh) * | 2018-07-12 | 2020-07-28 | 华南理工大学 | 一种基于深度神经网络和成对约束的聚类方法 |
CN109539359A (zh) * | 2018-10-15 | 2019-03-29 | 国网辽宁省电力有限公司 | 分工况pid+自适应前馈补偿的相变电蓄热供暖系统及方法 |
CN111798034A (zh) * | 2020-06-10 | 2020-10-20 | 河海大学 | 基于风流场时空图像学习的风电功率预测系统及预测方法 |
CN111798034B (zh) * | 2020-06-10 | 2021-09-07 | 河海大学 | 基于风流场时空图像学习的风电功率预测系统及预测方法 |
CN113404655A (zh) * | 2020-08-24 | 2021-09-17 | 湖南科技大学 | 一种基于ps0-anfis的风力发电机传感器状态诊断系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106845863A (zh) | 一种分布式风力发电出力与供暖负荷同步指数预测方法 | |
Azman et al. | A unified online deep learning prediction model for small signal and transient stability | |
CN112733462B (zh) | 一种结合气象因素的超短期风电场功率预测方法 | |
CN111179121B (zh) | 基于专家系统与深度逆向强化学习的电网紧急控制方法 | |
Kusiak et al. | Short-horizon prediction of wind power: A data-driven approach | |
JP5888640B2 (ja) | 太陽光発電予測装置、太陽光発電予測方法及び太陽光発電予測プログラム | |
CN108090621B (zh) | 一种基于分阶段整体优化的短期风速预测方法与系统 | |
Kassa et al. | Short term wind power prediction using ANFIS | |
CN105046374A (zh) | 一种基于核极限学习机模型的功率区间预测方法 | |
Zhang et al. | Application of simulated annealing genetic algorithm-optimized back propagation (BP) neural network in fault diagnosis | |
Gallego et al. | Maintaining flexibility in smart grid consumption through deep learning and deep reinforcement learning | |
CN116722545B (zh) | 基于多源数据的光伏发电预测方法和相关设备 | |
Annoni et al. | Sparse-sensor placement for wind farm control | |
CN109242200A (zh) | 一种贝叶斯网络预测模型的风电功率区间预测方法 | |
CN115688574A (zh) | 一种电力系统的暂态稳定在线评估方法及系统 | |
Zhang et al. | Safe building HVAC control via batch reinforcement learning | |
Sun et al. | Hybrid reinforcement learning for power transmission network self-healing considering wind power | |
Chen et al. | Distributed hierarchical deep reinforcement learning for large-scale grid emergency control | |
Naug et al. | Data driven methods for energy reduction in large buildings | |
Kumar et al. | Forecasting indoor temperature for smart buildings with ARIMA, SARIMAX, and LSTM: A fusion approach | |
CN105303319A (zh) | 一种基于运行数据的风电场动态分群等值方法 | |
Meng et al. | Hybrid neural network based on GRU with uncertain factors for forecasting ultra-short-term wind power | |
CN113784380A (zh) | 一种采用图注意力网络与融合邻域的拓扑预测方法 | |
CN105447658A (zh) | 一种含风电随机模糊注入电力系统波动的电压崩溃点求取方法 | |
CN113962357B (zh) | 基于gwo-wnn的分布式光伏功率数据虚拟采集方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170613 |