CN116128211A - Wind-light-water combined short-term optimization scheduling method based on wind-light uncertainty prediction scene - Google Patents
Wind-light-water combined short-term optimization scheduling method based on wind-light uncertainty prediction scene Download PDFInfo
- Publication number
- CN116128211A CN116128211A CN202211612301.0A CN202211612301A CN116128211A CN 116128211 A CN116128211 A CN 116128211A CN 202211612301 A CN202211612301 A CN 202211612301A CN 116128211 A CN116128211 A CN 116128211A
- Authority
- CN
- China
- Prior art keywords
- wind
- solar
- forecast
- scenario
- scene
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; 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
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Educational Administration (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Wind Motors (AREA)
Abstract
Description
技术领域Technical Field
本发明涉及考虑风光不确定性的风光水联合短期优化调度方法,具体是一种基于风光不确定性预报场景的风光水联合短期优化调度方法。The present invention relates to a wind-solar-water joint short-term optimization scheduling method considering the uncertainty of wind and solar power, and specifically to a wind-solar-water joint short-term optimization scheduling method based on a wind-solar-water uncertainty forecast scenario.
背景技术Background Art
可再生能源是绿色低碳能源,是我国多轮驱动能源供应体系的重要组成部分,对于改善能源结构、保护生态环境、应对气候变化、实现经济社会可持续发展具有重要意义。近十年来,可再生能源已成为我国发电新增装机的主体,我国可再生能源发电量增量在全社会用电量增量中的占比将超过50%,风电和太阳能发电量将再实现翻倍。但风光与光伏发电本身存在不确定性,直接接入将会给电网带来较大冲击,但利用水电灵活出力的特点,通过将风光水打捆外送,可在一定程度上平抑风光波动性,减少其不确定性带来的影响。Renewable energy is a green and low-carbon energy source, and is an important part of my country's multi-wheel-driven energy supply system. It is of great significance to improve the energy structure, protect the ecological environment, respond to climate change, and achieve sustainable economic and social development. In the past decade, renewable energy has become the main body of my country's new installed capacity for power generation. The increase in my country's renewable energy power generation will account for more than 50% of the increase in electricity consumption in the whole society, and wind power and solar power generation will double again. However, there is uncertainty in wind, solar and photovoltaic power generation. Direct access will bring a greater impact to the power grid. However, by taking advantage of the flexible output of hydropower, by bundling wind, solar and hydropower for external transmission, the volatility of wind and solar power can be smoothed to a certain extent, reducing the impact of its uncertainty.
目前,国内外学者对于考虑风光不确定性的风光水联合短期优化调度方法,研究工作主要集中在基于贝叶斯理论和模糊数学理论的建模方法。如刘乔波(2018年)基于历史数据对风电及光伏在不同出力水平下预测误差的分布进行拟合,引入spearman相关系数分季节描述风光预测误差的互补性,建立风光联合预测误差分布模型;张俊涛(2020年)基于随机规划理论,采用耦合分位点回归将历史统计信息由确定性预测序列转化为场景集;徐野驰(2022年)基于预测误差分布的风功率随机性分析法,通过生成描述风功率随机性的风功率出力场景,建立考虑频率响应的随机优化调度模型。然而,上述方法都是基于历史数据描述单独的场景或是多个不相关的独立场景用于制定优化调度方案,并不能考虑场景之间的相关性以及随时段更新的实测值对于调度系统的优化,而刻画风光不确定性的准确程度将会影响优化调度方案的制定,最终造成调度目标的偏差。At present, domestic and foreign scholars have focused their research on the short-term optimal dispatching method of wind, solar and water combined with the consideration of wind and solar uncertainty on modeling methods based on Bayesian theory and fuzzy mathematics theory. For example, Liu Qiaobo (2018) fitted the distribution of wind power and photovoltaic prediction errors at different output levels based on historical data, introduced the spearman correlation coefficient to describe the complementarity of wind and solar prediction errors by season, and established a wind and solar joint prediction error distribution model; Zhang Juntao (2020) used coupled quantile regression based on stochastic programming theory to transform historical statistical information from a deterministic prediction sequence into a scenario set; Xu Yechi (2022) established a random optimization dispatching model considering frequency response by generating wind power output scenarios describing the randomness of wind power based on the wind power randomness analysis method of prediction error distribution. However, the above methods are all based on historical data to describe a single scenario or multiple unrelated independent scenarios for the formulation of optimal dispatching plans, and cannot consider the correlation between scenarios and the optimization of the dispatching system by the measured values updated at different time periods. The accuracy of describing wind and solar uncertainty will affect the formulation of the optimal dispatching plan, and ultimately cause a deviation in the dispatching target.
发明内容Summary of the invention
发明目的:针对现有技术存在的问题与不足,本发明提供一种基于风光不确定性预报场景的风光水联合短期优化调度方法。基于历史风光出力的预测和实测数据,建立基于场景构建的风光场景预报模型、场景削减模型以及场景迁移概率计算模型,通过逐时更新风光预报典型场景集及其迁移概率系数,从而逐时计算更新优化调度决策,进行风光水联合短期优化调度。Purpose of the invention: In view of the problems and shortcomings of the prior art, the present invention provides a wind-solar-water joint short-term optimal scheduling method based on wind-solar uncertainty forecast scenarios. Based on the forecast and measured data of historical wind-solar output, a scene-based wind-solar scenario forecast model, a scene reduction model, and a scene migration probability calculation model are established. By updating the typical scene set of wind-solar forecasts and their migration probability coefficients hourly, the optimal scheduling decision is calculated and updated hourly, and the wind-solar-water joint short-term optimal scheduling is performed.
技术方案:一种基于风光不确定性预报场景的风光水联合短期优化调度方法,包括以下步骤:Technical solution: A wind-solar-water combined short-term optimization scheduling method based on wind-solar uncertainty forecast scenario, including the following steps:
S1、选取历史风光出力数据,对缺失值和异常值进行处理,对处理后的数据进行样本归一化,并划分样本层级;S1. Select historical wind and solar power output data, process missing values and outliers, normalize the processed data, and divide the sample levels;
S2、建立基于场景构建的风光场景预报模型,采用各层级的样本数据进行层级内的样本特征分析,运用分析结果基于风光出力的日前预报数据进行风光场景预报,得到风光预报场景集;S2. Establish a wind and solar scenario forecast model based on scenario construction, use sample data of each level to analyze sample characteristics within the level, use the analysis results to forecast wind and solar scenarios based on the day-ahead forecast data of wind and solar output, and obtain a wind and solar forecast scenario set;
S3、建立场景削减模型,对初始的风光预报场景集进行场景削减,迭代计算挑选最具有代表性的风光预报典型场景,组成风光预报典型场景集;S3, establishing a scenario reduction model, performing scenario reduction on the initial wind and solar forecast scenario set, iteratively calculating and selecting the most representative wind and solar forecast typical scenarios, and forming a wind and solar forecast typical scenario set;
S4、基于风光出力历史实测值建立场景迁移概率计算模型,对实测值进行区间划分,计算存在时间相关性的实测值之间的迁移概率;S4. Establish a scenario migration probability calculation model based on the historical measured values of wind and solar power output, divide the measured values into intervals, and calculate the migration probability between the measured values with time correlation;
S5、根据风光预报典型场景集中的结果,基于风光出力逐时实测数据,逐时更新优化调度决策,按照源荷匹配度最好原则,进行风光水联合短期优化调度。S5. According to the results of typical wind and solar forecast scenarios, based on the hourly measured data of wind and solar output, the optimized scheduling decision is updated hourly, and the short-term optimized scheduling of wind, solar and water is carried out according to the principle of best source-load matching.
进一步的,步骤S1具体为:历史风光出力数据主要包括,时空相关性一致的日前预报风出力数据以及实测风出力数据,以及时空相关性一致的日前预报光伏出力数据以及实测光伏出力数据;对于缺失值的处理方法为:按照相邻两组数据进行线性插值方法进行数据补充;对于异常值的处理方法为:首先对异常值进行删除处理,其缺失数据按照缺失值的方法补充;对处理后的数据进行样本归一化,样本归一化方法为:Further, step S1 is specifically as follows: the historical wind and solar power output data mainly include the day-ahead forecast wind output data and the measured wind output data with consistent temporal and spatial correlation, as well as the day-ahead forecast photovoltaic output data and the measured photovoltaic output data with consistent temporal and spatial correlation; the method for processing missing values is: the data is supplemented by the linear interpolation method according to two adjacent groups of data; the method for processing outliers is: firstly, the outliers are deleted, and the missing data is supplemented according to the missing value method; the processed data is sample normalized, and the sample normalization method is:
s′=(s-μ)/σs′=(s-μ)/σ
式中,s′为归一化后样本值;s为样本原始数据;μ为样本原始集的均值;σ为样本原始集的标准差。Where s′ is the normalized sample value; s is the original sample data; μ is the mean of the original sample set; σ is the standard deviation of the original sample set.
进一步的,步骤S2中建立的基于场景构建的风光场景预报模型结构为:Furthermore, the structure of the scene-based forecast model established in step S2 is:
基于样本归一化处理后得到的预报与实测数据集,将其按照预报值数值大小进行分级,计算每一层级内数据特征值,特征值计算方法如下:Based on the forecast and measured data sets obtained after sample normalization, they are classified according to the numerical value of the forecast value, and the characteristic value of the data in each level is calculated. The characteristic value calculation method is as follows:
式中,ε′i,j为第i层级内的第j个样本相对误差值;sp′i,j为预测样本的归一化值;st′i,j为实测样本的归一化值。Where ε′ i,j is the relative error value of the jth sample in the i-th level; sp′ i,j is the normalized value of the predicted sample; and st′ i,j is the normalized value of the measured sample.
式中,为层级的第一特征值;ni为第i层级的样本相对误差值数量。In the formula, is the first eigenvalue of the level; ni is the number of relative error values of samples in the i-th level.
式中,为层级的第二特征值。In the formula, is the second eigenvalue of the level.
进一步的,基于该样本层级中的两个特征值,采用超拉丁方抽样生成风光出力预报场景集。Furthermore, based on the two eigenvalues in the sample level, super Latin square sampling is used to generate a set of wind and solar output forecast scenarios.
进一步的,步骤S3具体为:采用K-means随机质心聚类法进行场景削减,聚类结果控制在4-8之间。Furthermore, step S3 specifically includes: using K-means random centroid clustering method to reduce the scenes, and controlling the clustering result between 4 and 8.
进一步的,步骤S4具体为:根据风光出力实测数据归一化结果进行区间划分,建立场景迁移概率计算模型,概率计算方法如下:Further, step S4 is specifically: according to the normalized results of the wind and solar power output measured data, the interval is divided and a scene migration probability calculation model is established. The probability calculation method is as follows:
st″=(st′-μ′)/σ′st″=(st′-μ′)/σ′
式中,st″表示标准正态分布的标准分;μ′表示实测样本归一化值均值;σ′表示实测样本归一化值标准差。Where st″ represents the standard score of the standard normal distribution; μ′ represents the mean of the normalized values of the measured samples; σ′ represents the standard deviation of the normalized values of the measured samples.
式中,Pst′(u,d)表示当前实测值归一化后位于(u,d)区间时,下一实测值归一化后为st′的概率。Where Pst′ (u, d) represents the probability that the next measured value will be st′ after normalization when the current measured value is in the interval (u, d) after normalization.
进一步的,步骤S5具体为:基于风光预报典型场景集中的结果,逐时计算各个典型场景发生的概率,不断生成用于短期优化调度计算的决策场景,计算方法如下:Further, step S5 is specifically as follows: based on the results of the typical scene set of wind and solar forecast, the probability of occurrence of each typical scene is calculated hourly, and the decision scene for short-term optimization scheduling calculation is continuously generated. The calculation method is as follows:
式中,P′q为风光预报典型场景集中第q个场景的迁移概率系数。Where P′q is the migration probability coefficient of the qth scene in the typical scene set of wind and solar forecast.
式中,x′p表示优化调度期内第p个小时的风光预报场景值;xp,q表示风光预报典型场景值中第w个场景中第p个小时的值。Where x′p represents the wind and solar forecast scenario value of the pth hour during the optimization scheduling period; xp,q represents the value of the pth hour in the wth scenario of the typical wind and solar forecast scenario value.
进一步的,以源荷匹配度最优为目标构建风光水联合短期优化调度模型,基于不断生成的短期优化调度的决策场景,实现考虑风光不确定性的风光水联合短期优化调度,计算方法如下:Furthermore, a wind-solar-water joint short-term optimal scheduling model is constructed with the goal of optimizing source-load matching. Based on the continuously generated decision scenarios of short-term optimal scheduling, the wind-solar-water joint short-term optimal scheduling considering the uncertainty of wind and solar is realized. The calculation method is as follows:
Nhp=Nw+Ns+Nh Nh p =N w +N s +N h
式中,Nv为源荷匹配率;Nyp为第p个小时的系统总出力;Nhp为第p个小时的系统负荷;Nw、Ns、Nh分别为风电、光电、水电出力值;为水电站第j时的发电流量;Hj为水电站第j时的发电水头。Where, Nv is the source-load matching rate; Nyp is the total system output at the pth hour; Nhp is the system load at the pth hour; Nw , Ns , and Nh are the wind power, photovoltaic, and hydropower output values, respectively; is the power generation flow of the hydropower station at the jth time; Hj is the power generation head of the hydropower station at the jth time.
一种基于风光不确定性预报场景的风光水联合短期优化调度系统,包括:A wind-solar-water combined short-term optimization dispatching system based on wind-solar uncertainty forecasting scenarios, comprising:
数据处理模块,用于对输入的历史风光出力数据进行层级划分,并对各层级的历史数据进行特征值的计算,用于场景生成;The data processing module is used to classify the input historical wind and solar power output data and calculate the characteristic values of the historical data at each level for scene generation;
场景生成模块,用于建立风光场景预报模型与场景削减模型,基于日前预报数据采用各层级样本的特征值生成场景集,运用迭代聚类进行场景集的削减,输出风光预报典型场景集;The scenario generation module is used to establish a wind and solar scenario forecast model and a scenario reduction model. Based on the day-ahead forecast data, the characteristic values of samples at each level are used to generate a scenario set, and iterative clustering is used to reduce the scenario set, and a typical wind and solar forecast scenario set is output;
场景迁移概率计算模块,用于计算存在时间相关性的实测值之间的迁移概率,输出风光预报典型场景集中各场景的迁移概率系数;The scene migration probability calculation module is used to calculate the migration probability between measured values with time correlation and output the migration probability coefficient of each scene in the typical scene set of wind and solar forecast;
优化调度模块,用于制定风光水联合优化调度计划,基于风光预报典型场景集与各场景迁移概率系数,按照源荷匹配度最好原则,进行风光水联合短期优化调度计算。The optimization scheduling module is used to formulate a wind, solar and water joint optimization scheduling plan. Based on the typical scene set of wind and solar forecasts and the migration probability coefficient of each scene, it performs wind, solar and water joint short-term optimization scheduling calculations in accordance with the principle of best source-load matching.
一种计算机设备,该计算机设备包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,处理器执行上述计算机程序时实现如上所述的基于风光不确定性预报场景的风光水联合短期优化调度方法。A computer device includes a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the wind-solar-water combined short-term optimization scheduling method based on the wind-solar uncertainty forecast scenario as described above is implemented.
一种计算机可读存储介质,该计算机可读存储介质存储有执行如上所述的基于风光不确定性预报场景的风光水联合短期优化调度方法的计算机程序。A computer-readable storage medium stores a computer program for executing the wind-solar-water combined short-term optimization scheduling method based on the wind-solar uncertainty forecast scenario as described above.
有益效果:与现有技术相比,本发明的技术效果为:(1)本发明基于历史风光出力的预测和实测数据,建立基于场景构建的风光场景预报模型、场景削减模型以及场景迁移概率计算模型,通过生成风光典型预报场景集,制定风光水联合调度计划,并据此提出了新的考虑不确定性的风光水联合短期优化调度方法;(2)通过逐时更新风光预报典型场景集及其迁移概率系数,可以逐时计算更新优化调度决策。该方法不确定信息提取部分所需数据均为历史数据,可编写PYTHON程序应对任一调度期不同输入数据情况下风光水联合短期优化调度计算;(3)原理简单,操作简便灵活,易于实施。该技术方法基于场景构建与削减以及场景迁移概率模型进行优化调度计算,计算速度快,响应时间短;(4)可以为后续其他可再生能源的接入提供支撑;该技术为接入新能源的多能互补系统优化调度提供了新的方法和技术,也为多能互补系统新能源消纳、减少电网冲击、优化能源结构打下了较好的基础。Beneficial effects: Compared with the prior art, the technical effects of the present invention are as follows: (1) Based on the forecast and measured data of historical wind and solar output, the present invention establishes a scene-based wind and solar forecast model, a scene reduction model and a scene migration probability calculation model, and formulates a wind, solar and water joint dispatch plan by generating a typical wind and solar forecast scene set, and proposes a new wind, solar and water joint short-term optimal dispatch method considering uncertainty; (2) By updating the typical wind and solar forecast scene set and its migration probability coefficient hour by hour, the optimized dispatch decision can be calculated and updated hour by hour. The data required for the uncertain information extraction part of this method are all historical data, and a PYTHON program can be written to cope with the wind, solar and water joint short-term optimal dispatch calculation under different input data in any dispatch period; (3) The principle is simple, the operation is simple and flexible, and it is easy to implement. This technical method performs optimized dispatch calculation based on scene construction and reduction and scene migration probability model, with fast calculation speed and short response time; (4) It can provide support for the subsequent access of other renewable energy sources; This technology provides a new method and technology for the optimal dispatch of multi-energy complementary systems accessing new energy, and also lays a good foundation for the consumption of new energy in multi-energy complementary systems, reducing grid shocks and optimizing energy structure.
附图说明BRIEF DESCRIPTION OF THE DRAWINGS
图1为本发明的方法流程图;Fig. 1 is a flow chart of the method of the present invention;
图2为本发明的系统结构框图。FIG. 2 is a system structure block diagram of the present invention.
具体实施方式DETAILED DESCRIPTION
下面结合具体实施例,进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。The present invention is further explained below in conjunction with specific embodiments. It should be understood that these embodiments are only used to illustrate the present invention and are not used to limit the scope of the present invention. After reading the present invention, various equivalent forms of modifications to the present invention by those skilled in the art all fall within the scope defined by the claims attached to this application.
本发明基于历史风光出力的预测和实测数据,建立基于场景构建的风光场景预报模型、场景削减模型以及场景迁移概率计算模型,通过逐时更新风光预报典型场景集及其迁移概率系数,从而逐时计算更新优化调度决策,进行风光水联合短期优化调度。Based on the prediction and measured data of historical wind and solar power output, the present invention establishes a wind and solar scene forecast model, a scene reduction model and a scene migration probability calculation model based on scenario construction. By updating the typical scene set of wind and solar forecast and its migration probability coefficient hourly, the optimized scheduling decision is calculated and updated hourly to carry out the short-term optimal scheduling of wind, solar and water joint.
如图1所示,本发明的一种基于风光不确定性预报场景的风光水联合短期优化调度方法,包括以下步骤:As shown in FIG1 , a wind-solar-water joint short-term optimization scheduling method based on a wind-solar uncertainty forecast scenario of the present invention comprises the following steps:
S1、选取历史风光出力数据,对缺失值和异常值进行处理,对处理后的数据进行样本归一化,并划分样本层级,具体步骤如下:S1. Select historical wind and solar output data, process missing values and outliers, normalize the processed data, and divide the sample levels. The specific steps are as follows:
S1.1、首先对历史风光出力数据进行筛选,选取相关风电场与光伏电站附近的大气环境测点数据进行计算,数据精度控制在1-15min,同时保证数据样本的总体一致性,以确保其时空相关性。S1.1. First, the historical wind and solar power output data are screened, and the atmospheric environment measurement point data near the relevant wind farms and photovoltaic power stations are selected for calculation. The data accuracy is controlled within 1-15min, and the overall consistency of the data samples is guaranteed to ensure its temporal and spatial correlation.
S1.2、筛选得到历史风光出力预测和实测数据后,需要对其进行缺失值与异常值的处理,对于风光出力数据而言,将不会出现负值以及超出测量仪器量程的数据,特别的,对于光伏出力数据,其在太阳升起前的测量值应始终为零,因此,选取的历史数据需要完整且相对合理的风光出力的预测与实测数据。S1.2. After screening the historical wind and solar power output forecast and measured data, it is necessary to process the missing values and outliers. For wind and solar power output data, there will be no negative values or data exceeding the range of the measuring instrument. In particular, for photovoltaic output data, the measured value before the sun rises should always be zero. Therefore, the selected historical data requires complete and relatively reasonable wind and solar power output forecast and measured data.
S1.3、对处理后的样本进行归一化处理,样本归一化方法为:S1.3. Normalize the processed samples. The sample normalization method is:
s′=(s-μ)/σs′=(s-μ)/σ
式中,s′为归一化后样本值;s为样本原始数据;μ为样本原始集的均值;σ为样本原始集的标准差。Where s′ is the normalized sample value; s is the original sample data; μ is the mean of the original sample set; σ is the standard deviation of the original sample set.
S2、建立基于场景构建的风光场景预报模型,采用各层级的样本数据进行层级内的样本特征分析,运用分析结果基于风光出力的日前预报数据进行风光场景预报,得到风光预报场景集,具体步骤如下:S2. Establish a wind and solar scenario forecast model based on scenario construction, use sample data of each level to analyze sample characteristics within the level, use the analysis results to forecast wind and solar scenarios based on the day-ahead forecast data of wind and solar output, and obtain a wind and solar forecast scenario set. The specific steps are as follows:
S2.1、计算数据集中各时间点预报值的预报误差,根据预报误差的统计结果进行层级划分,可采用误差值进行梯级划分或按误差样本数量均分,确定层级内样本上下边界用于后续抽样进行场景生成。S2.1. Calculate the forecast error of the forecast value at each time point in the data set, and divide the layers according to the statistical results of the forecast error. The error value can be used for ladder division or the number of error samples can be evenly divided to determine the upper and lower boundaries of the samples within the layer for subsequent sampling for scene generation.
S2.2、计算每一层级内数据特征值,特征值计算方法如下:S2.2. Calculate the characteristic value of the data in each level. The characteristic value calculation method is as follows:
式中,ε′i,j为第i层级内的第j个样本相对误差值;sp′i,j为预测样本的归一化值;st′i,j为实测样本的归一化值。Where ε′ i,j is the relative error value of the jth sample in the i-th level; sp′ i,j is the normalized value of the predicted sample; and st′ i,j is the normalized value of the measured sample.
式中,为层级的第一特征值;ni为第i层级的样本相对误差值数量。In the formula, is the first eigenvalue of the level; ni is the number of relative error values of samples in the i-th level.
式中,为层级的第二特征值。In the formula, is the second eigenvalue of the level.
S2.3、基于该样本层级中的两个特征值,采用超拉丁方抽样生成风光出力预报场景集。超拉丁方抽样(LHS)是通过较少迭代次数的抽样,准确建立输入分布,与蒙特卡洛发相比,其关键在于对输入概率进行分层,并建立累计概率尺度相等的区间,通过从输入分布的每个区间或层中进行样本抽取,抽样结果被强制代表每个区间的值,因此输入概率分布被强制重建。风光与光伏预报场景集的生成策略为:整个抽样过程采用“抽样不替换”原则,累计分布的分层数目与整个执行过程中的迭代次数相同。需要注意的是,当使用该方法进行抽样时,需要保持变量之间的独立性,独立性的保持通过对每个变量随机选择抽样区间来实现,这样可以有效避免变量间的无意相关。S2.3. Based on the two eigenvalues in the sample level, super Latin square sampling is used to generate a set of wind and solar output forecast scenarios. Super Latin square sampling (LHS) accurately establishes the input distribution through sampling with a smaller number of iterations. Compared with Monte Carlo, the key is to stratify the input probability and establish intervals with equal cumulative probability scales. By extracting samples from each interval or layer of the input distribution, the sampling results are forced to represent the values of each interval, so the input probability distribution is forced to be reconstructed. The generation strategy of the wind and photovoltaic forecast scenario set is: the entire sampling process adopts the "sampling without replacement" principle, and the number of stratifications of the cumulative distribution is the same as the number of iterations in the entire execution process. It should be noted that when using this method for sampling, the independence between variables needs to be maintained. The independence is maintained by randomly selecting sampling intervals for each variable, which can effectively avoid unintentional correlation between variables.
S3、建立场景削减模型,对初始的风光预报场景集进行场景削减,迭代计算挑选最具有代表性的风光预报典型场景,组成风光预报典型场景集。S3. Establish a scenario reduction model, perform scenario reduction on the initial wind and solar forecast scenario set, iteratively calculate and select the most representative wind and solar forecast typical scenarios to form a wind and solar forecast typical scenario set.
K-means算法作为聚类中的基础算法,属于无监督学习中的一种算法,其基本原理是随机确定k个初始点作为簇质心,然后计算样本数据中各点与簇质心之间的计算距离,并以此距离为依据对样本进行分类,通过不断迭代簇质心的位置,可以得到更好的聚类结果。The K-means algorithm is a basic algorithm in clustering and an algorithm in unsupervised learning. Its basic principle is to randomly determine k initial points as cluster centroids, then calculate the distance between each point in the sample data and the cluster centroid, and classify the samples based on this distance. By continuously iterating the position of the cluster centroid, better clustering results can be obtained.
具体步骤如下:The specific steps are as follows:
S3.1、根据初始风光预报场景集,我们选择4-8个初始质心,即为典型簇质心。S3.1. Based on the initial set of scenery forecast scenes, we select 4-8 initial centroids, which are the typical cluster centroids.
S3.2、计算整个场景集中其他元素距离典型簇质心的距离,这里距离采用欧式距离进行计算。S3.2. Calculate the distance between other elements in the entire scene set and the centroid of the typical cluster. Here, the distance is calculated using Euclidean distance.
S3.3、选择新的簇质心并计算在新的簇质心下整个场景集的分类情况,通过不断重复S3.2-S3.3,指导整个场景集的分类结果收敛,则表明聚类结束。S3.3. Select a new cluster centroid and calculate the classification of the entire scene set under the new cluster centroid. By continuously repeating S3.2-S3.3, the classification results of the entire scene set converge, indicating that clustering is completed.
S4、基于风光出力历史实测值建立场景迁移概率计算模型,对实测值进行区间划分,计算存在时间相关性的实测值之间的迁移概率,具体步骤如下:S4. A scenario migration probability calculation model is established based on the historical measured values of wind and solar power output, and the measured values are divided into intervals to calculate the migration probability between the measured values with time correlation. The specific steps are as follows:
S4.1、计算典型场景集中各个场景的对应情况的正态分布标准分,计算方法如下:S4.1. Calculate the normal distribution standard score of the corresponding situation of each scenario in the typical scenario set. The calculation method is as follows:
st″=(st′-μ′)/σ′st″=(st′-μ′)/σ′
式中,st″表示标准正态分布的标准分;μ′表示实测样本归一化值均值;σ′表示实测样本归一化值标准差。Where st″ represents the standard score of the standard normal distribution; μ′ represents the mean of the normalized values of the measured samples; σ′ represents the standard deviation of the normalized values of the measured samples.
S4.2、根据风光实测数据归一化结果进行区间划分,建立场景迁移概率计算模型,概率计算方法如下:S4.2. According to the normalized results of the wind and solar measured data, the intervals are divided and a scene migration probability calculation model is established. The probability calculation method is as follows:
式中,Pst′(u,d)表示当前实测值归一化后位于(u,d)区间时,下一实测值归一化后为st′的概率。Where Pst′ (u, d) represents the probability that the next measured value will be st′ after normalization when the current measured value is in the interval (u, d) after normalization.
S5、根据风光预报典型场景集中的结果,基于风光出力逐时实测数据,逐时更新优化调度决策,按照源荷匹配度最好原则,进行风光水联合短期优化调度,具体步骤如下:S5. According to the results of the typical scene of wind and solar forecast, based on the hourly measured data of wind and solar output, the optimized dispatch decision is updated hourly, and the short-term optimized dispatch of wind, solar and water is carried out according to the principle of the best source-load matching. The specific steps are as follows:
S5.1、基于风光预测典型场景集中各个场景迁移概率值,计算各个场景的迁移概率系数,计算方法如下:S5.1. Based on the migration probability values of each scene in the typical scene set of wind and solar forecast, the migration probability coefficient of each scene is calculated. The calculation method is as follows:
式中,P′q为风光预报典型场景集中第q个场景的迁移概率系数。Where P′q is the migration probability coefficient of the qth scene in the typical scene set of wind and solar forecast.
S5.2、基于各个场景的迁移概率系数,生成用于短期优化调度的计算场景,计算方法如下:S5.2. Based on the migration probability coefficient of each scenario, generate a calculation scenario for short-term optimization scheduling. The calculation method is as follows:
式中,x′p表示优化调度期内第p个小时的风光预报场景值;xp,q表示风光预报典型场景值中第q个场景中第p个小时的值。Where x′p represents the wind and solar forecast scenario value of the pth hour during the optimization scheduling period; xp,q represents the value of the pth hour in the qth scenario of the wind and solar forecast typical scenario value.
S5.3、以源荷匹配度最优为目标构建风光水联合短期优化调度模型,基于不断生成的短期优化调度的决策场景,实现考虑风光不确定性的风光水联合短期优化调度,计算方法如下:S5.3. A wind-solar-water joint short-term optimal dispatch model is constructed with the goal of optimizing source-load matching. Based on the continuously generated decision scenarios for short-term optimal dispatch, a wind-solar-water joint short-term optimal dispatch considering the uncertainty of wind and solar is realized. The calculation method is as follows:
Nhp=Nw+Ns+Nh Nh p =N w +N s +N h
式中,Nv为源荷匹配率;Nyp为第p个小时的系统总出力;Nhp为第p个小时的系统负荷;Nw、Ns、Nh分别为风电、光电、水电出力值;为水电站第j时的发电流量;Hj为水电站第j时的发电水头。Where, Nv is the source-load matching rate; Nyp is the total system output at the pth hour; Nhp is the system load at the pth hour; Nw , Ns , and Nh are the wind power, photovoltaic, and hydropower output values, respectively; is the power generation flow of the hydropower station at the jth time; Hj is the power generation head of the hydropower station at the jth time.
本发明基于历史风光出力的预测和实测数据,建立基于场景构建的风光场景预报模型、场景削减模型以及场景迁移概率计算模型,通过逐时更新风光预报典型场景集及其迁移概率系数,从而逐时计算更新优化调度决策,进行风光水联合短期优化调度。该技术为接入新能源的多能互补系统优化调度提供了新的方法和技术,也为多能互补系统新能源消纳、减少电网冲击、优化能源结构打下了较好的基础。Based on the forecast and measured data of historical wind and solar output, the present invention establishes a scene-based wind and solar scenario forecast model, a scene reduction model, and a scene migration probability calculation model. By updating the typical scene set of wind and solar forecast and its migration probability coefficient hourly, the optimized scheduling decision is calculated and updated hourly, and the short-term optimized scheduling of wind, solar and water joint is performed. This technology provides new methods and technologies for the optimized scheduling of multi-energy complementary systems connected to new energy, and also lays a good foundation for the consumption of new energy in multi-energy complementary systems, reducing grid shocks, and optimizing energy structure.
如图2所示,基于风光不确定性预报场景的风光水联合短期优化调度系统,包括:As shown in Figure 2, the wind-solar-water combined short-term optimization dispatching system based on the wind-solar uncertainty forecast scenario includes:
数据处理模块,用于对输入的数据进行层级划分,并对各层级的历史数据进行特征值的计算,用于场景生成;The data processing module is used to divide the input data into layers and calculate the characteristic values of the historical data at each layer for scenario generation;
场景生成模块,用于建立风光场景预报模型与场景削减模型,基于日前预报数据采用各层级样本的特征值生成场景集,运用迭代聚类进行场景集的削减,输出风光预报典型场景集;The scenario generation module is used to establish a wind and solar scenario forecast model and a scenario reduction model. Based on the day-ahead forecast data, the characteristic values of samples at each level are used to generate a scenario set, and iterative clustering is used to reduce the scenario set, and a typical wind and solar forecast scenario set is output;
场景迁移概率计算模块,用于计算存在时间相关性的实测值之间的迁移概率,输出风光预报典型场景集中各场景的迁移概率系数;The scene migration probability calculation module is used to calculate the migration probability between measured values with time correlation and output the migration probability coefficient of each scene in the typical scene set of wind and solar forecast;
优化调度模块,用于制定风光水联合优化调度计划,基于风光预报典型场景集与各场景迁移概率系数,按照源荷匹配度最好原则,进行风光水联合短期优化调度计算。The optimization scheduling module is used to formulate a wind, solar and water joint optimization scheduling plan. Based on the typical scene set of wind and solar forecasts and the migration probability coefficient of each scene, it performs wind, solar and water joint short-term optimization scheduling calculations in accordance with the principle of best source-load matching.
系统的实现过程与方法一样,不再赘述。The implementation process and method of the system are the same and will not be repeated here.
显然,本领域的技术人员应该明白,上述的本发明实施例的基于风光不确定性预报场景的风光水联合短期优化调度方法各步骤基于风光不确定性预报场景的风光水联合短期优化调度系统各模块可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,并且在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明实施例不限制于任何特定的硬件和软件结合。Obviously, those skilled in the art should understand that the various steps of the wind-solar-water joint short-term optimal scheduling method based on the wind-solar uncertainty forecast scenario and the various modules of the wind-solar-water joint short-term optimal scheduling system based on the wind-solar uncertainty forecast scenario of the above-mentioned embodiment of the present invention can be implemented by a general computing device, they can be concentrated on a single computing device, or distributed on a network composed of multiple computing devices, optionally, they can be implemented with executable program codes of the computing device, so that they can be stored in a storage device and executed by the computing device, and in some cases, the steps shown or described can be executed in a different order than here, or they can be made into individual integrated circuit modules respectively, or multiple modules or steps therein can be made into a single integrated circuit module for implementation. In this way, the embodiments of the present invention are not limited to any specific combination of hardware and software.
Claims (10)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211612301.0A CN116128211A (en) | 2022-12-15 | 2022-12-15 | Wind-light-water combined short-term optimization scheduling method based on wind-light uncertainty prediction scene |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211612301.0A CN116128211A (en) | 2022-12-15 | 2022-12-15 | Wind-light-water combined short-term optimization scheduling method based on wind-light uncertainty prediction scene |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116128211A true CN116128211A (en) | 2023-05-16 |
Family
ID=86309127
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211612301.0A Pending CN116128211A (en) | 2022-12-15 | 2022-12-15 | Wind-light-water combined short-term optimization scheduling method based on wind-light uncertainty prediction scene |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116128211A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116960933A (en) * | 2023-06-13 | 2023-10-27 | 国网山东省电力公司威海供电公司 | An optimized dispatching method and system for marine ranch microgrid |
| CN117154725A (en) * | 2023-10-31 | 2023-12-01 | 长江三峡集团实业发展(北京)有限公司 | Scheduling methods, devices, computer equipment and media for water, wind and solar multi-energy complementation |
| CN119271845A (en) * | 2024-04-03 | 2025-01-07 | 河海大学 | A method for generating multi-stage scene tree of wind-solar-load based on uncertainty feature extraction of source load |
-
2022
- 2022-12-15 CN CN202211612301.0A patent/CN116128211A/en active Pending
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116960933A (en) * | 2023-06-13 | 2023-10-27 | 国网山东省电力公司威海供电公司 | An optimized dispatching method and system for marine ranch microgrid |
| CN117154725A (en) * | 2023-10-31 | 2023-12-01 | 长江三峡集团实业发展(北京)有限公司 | Scheduling methods, devices, computer equipment and media for water, wind and solar multi-energy complementation |
| CN117154725B (en) * | 2023-10-31 | 2024-01-26 | 长江三峡集团实业发展(北京)有限公司 | Water-wind-solar multi-energy complementary scheduling method, device, computer equipment and medium |
| CN119271845A (en) * | 2024-04-03 | 2025-01-07 | 河海大学 | A method for generating multi-stage scene tree of wind-solar-load based on uncertainty feature extraction of source load |
| CN119271845B (en) * | 2024-04-03 | 2025-11-07 | 河海大学 | Wind-solar-load multi-stage scene tree generation method based on source-load uncertainty feature extraction |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Cai et al. | Wind speed forecasting based on extreme gradient boosting | |
| Chen et al. | Effective load carrying capability evaluation of renewable energy via stochastic long-term hourly based SCUC | |
| CN109615146B (en) | Ultra-short-term wind power prediction method based on deep learning | |
| CN116128211A (en) | Wind-light-water combined short-term optimization scheduling method based on wind-light uncertainty prediction scene | |
| CN113255973A (en) | Power load prediction method, power load prediction device, computer equipment and storage medium | |
| CN112100911A (en) | A Solar Radiation Prediction Method Based on Deep BISLTM | |
| CN109840633B (en) | Photovoltaic output power prediction method, system and storage medium | |
| CN105631558A (en) | BP neural network photovoltaic power generation system power prediction method based on similar day | |
| CN108429256B (en) | Power system operation optimization method and terminal equipment | |
| CN115618922A (en) | Photovoltaic power prediction method, equipment, photovoltaic power generation system and storage medium | |
| CN114819374B (en) | Regional new energy ultra-short term power prediction method and system | |
| CN117713238B (en) | Stochastic Optimal Operation Strategy of Photovoltaic Power Generation and Energy Storage Microgrid | |
| Xu et al. | Adaptive feature selection and GCN with optimal graph structure-based ultra-short-term wind farm cluster power forecasting method | |
| CN116014722A (en) | Sub-solar photovoltaic power generation prediction method and system based on seasonal decomposition and convolution network | |
| Su et al. | Short-term transmission capacity prediction of hybrid renewable energy systems considering dynamic line rating based on data-driven model | |
| CN115511170A (en) | A Modeling Method for Power Prediction Error of Multiple Photovoltaic Power Stations | |
| Wang et al. | An ultra-short-term forecasting model for high-resolution solar irradiance based on SOM and deep learning algorithm | |
| Zim et al. | Short-term weather forecasting for wind energy generation using a deep learning technique | |
| Fu et al. | A spatial forecasting method for photovoltaic power generation combined of improved similar historical days and dynamic weights allocation | |
| Yang et al. | Hybrid data and model‐driven joint identification of distribution‐network topology and parameters | |
| CN119518681A (en) | Method for determining photovoltaic power generation, computer program product and electronic device | |
| CN108460501B (en) | Wind power station output power prediction method based on combined model | |
| Liu et al. | A novel photovoltaic power output forecasting method based on weather type clustering and wavelet support vector machines regression | |
| CN117808148A (en) | Method for predicting uncertainty of micro-grid based on signal decomposition and deep learning | |
| CN116644903A (en) | Method, system, equipment and storage medium for assessing peak shaving adequacy of power system |
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 |