西电数模A题风功率

1、承诺书我们完全明白，在竞赛开始后参赛队员不能以任何方式包括 、电子邮件、网上咨询等与队外的任何人研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规那么的, 如果引用别人的成果或其他公开的资料包括网上查到的资料，必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规那么，以保证竞赛的公正、公平性。如有违反竞赛规那么的行为，我们将受到严肃处理。我们授权西安电子科技大学2021年数学建模校内赛竞赛组委会，可将我们的论文以任何形式进行公开展示包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等。我们参赛选择的题号是从A/B/C中选择一项填

2、写： A 参赛报名号为 物光051 报名时所属学院请填写完整的全名： 物理与光电工程学院 参赛队员与学号 (打印，用二号字，并签名) ：1. 姜垚2. 魏春泉3. 曾增峰日期： 年 月 日西安电子科技大学2021年大学生数学建模校内赛评阅专用页评阅人1评阅人2评阅人3总评成绩实用文档风电场功率预测与规划摘 要我们用数学建模方法解决风电场功率预测问题。然后选取一个可建立风电场的地区,对其进行规划,使所建风电场能产出最大风电总功率。 首先对于问题一，由于风具有节歇性特点，且短时间内风速具有一定的稳定性，我们决定采取灰色预

3、测模型，选用当日最后24个数据作为历史数据进行预测。先通过构造关联矩阵，从而对其进行累计加权，最终得到了合理的灰色预测模型。通过对实际功率和预测功率进行相关性分析，计算出准确率和合格率。将实际数据和预测数据进行比拟得出灰色预测根本能够预测未来四个小时内的风电功率。 其次对于问题二，是对未来两天的预测，我们取预测点前两天的数据作为一组时间序列，通过对初始数据的零均值转换，使得数据趋于平稳，然后用统计软件对数据分析计算得到相关系数，进而对模型确定和定阶，经过拟合最终得到合理的时间序列ARMA模型。通过对预测数据进行相对误差分析发现时间序列ARMA模型能够准确的预测出未来48个小时的风电功率。最后对

4、于问题三，是对风电场的选址和规划，使所建风电场能产出最大的风电总功率。我们首先运用宏观选址的方法，通过对我国各地区风能资源分布情况和风电场运行数据的搜集，决定选取福建东山建立风电场。最后通过风机尾流模型，风机本钱-效益模型和风电场增量装机效益评估模型来实现风电场风机的优化布置。最后使所建风电场能产出最大风电总功率。以上求得的结果具有较高的精确性和可行性，对于风电功率的预测，风电场的选址和规划具有一定的指导意义。 关键词：风电功率 灰色预测 时间序列 宏观选址 优化布置一、问题的重述1.1背景资料与条件对于没有风电资源的电力系统来说,电网会根据日负荷曲线制定运行方案，以满足次日的电力需求。对于风

5、电富有地区的电力系统来说,由于风电具有很大的随机性、间歇性和不可控性,同时在实际运行中还常常具有反调峰特性，这给电网调度,维持发电、输电、用电之间的功率平衡以及电力市场管理带来了极大的压力。要解决大规模风电对电网功率平衡所带来的这种挑战,目前可以做的工作之一就是对风电场的输出功率进行预测,通过建立物理、数学、组合多种方法进行预测，提高预测结果的精度，以便将风电功率合理的纳入电网的调度方案中。 1.2需要解决的问题1利用附件1.1中的数据建立风电功率的超短期预测模型,能够预测未来4小时内的风电功率,给出具体预测结果,并将预测结果与附件1.2中的实际功率比拟。2利用附件2中的数据建立短期风电功率预

6、测模型,模型能够预测未来48小时内的风电功率,时间间隔为15分钟,给出预测结果和误差分析。3在进行风电入网规划时,需要将一个风电场所产出的总风电功率看作整体进行调度,在每组风电功率预测的根底上,预测风电场全区的风电总功率。根据我国风力资源、地理环境等信息,确定出一个可以建立风电场的地区,并规划该地区的风电机组位置、规模等,使所建的风电场能产出最大的风电总功率,给出具体结果。二、问题分析2.1问题的重要性分析可再生能源的开发和利用可以有效保护生态环境、保证能源平安。风能是目前最成熟、最具开发规模和商业化开展前景的可再生能源之一。我国自1986年在山东建立了第一座并网型风力发电场后，风力发电开始逐

7、步开展。对于没有风电资源的电力系统来说，电网会根据日负荷曲线制定运行方案，以满足次日的电力需求。但是对于风电富有地区的电力系统来说，由于风电具有很大的随机性、间歇性和不可控性，同时在实际运行中还常常具有反调峰特性即负荷大时风电出力小，负荷小时风电出力反而大，这给电网调度，维持发电、输电、用电之间的功率平衡以及电力市场管理带来了极大的压力。要解决大规模风电对电网功率平衡所带来的这种挑战，目前可以做的工作之一就是对风电场的输出功率进行预测，以便将风电功率合理的纳入电网的调度方案中。风电场输出功率可以建立物理、数学等方法进行预测，也可以组合多种方法进行预测，提高预测结果的精度是风电功率预测永远追求的

8、目标之一2.2问题的思路分析此题研究的是风电功率预测和建模问题对短期个超短期风电功率进行准确的预测，然后根据预测结果进行选址和优化布置，实现风电场功率最大输出。2.2.1问题一该题是以通过建模进行超短期功率预测为研究目的。该题属于预测类数学问题可以直接利用历史数据进行短期预测，而常见方法有人工神经网络法、时间序列法、灰色分析预测法等。基于数据量小等应选择灰色分析预测法建立数学模型。2.2.2问题二该题是通过短期预测获得功率数据并对其进行误差分析，由于数据量相对于第一题稍多故采用时间序列分析预测法来建立数学模型。2.2.3问题三该题本质是对风电机组进行最正确选址，基于影响因素众多，应选用宏观选址

9、并进行配置优化以到达目的。 三、模型假设3.1假设一:题目中所给的数据合理正确3.2假设二:所有发电机性能相同3.3假设三:测数据时发电机都正常工作3.4假设四:软件求解过程中系统不会产生误差3.5假设五:天气因素不影响风电功率四、符号说明符号符号说明预测序列样本序列差分处理后的序列准确率合格率自相关系数偏自相关系数预测值平均绝对误差五、模型的建立与求解5.1、问题一(对未来4小时超短期风功率预测)5.1.1模型的分析1原始数据，原始数据选取风电机组后6个小时的的24个数据点；2计算生成序列，用GM(1,1)建模时，首先对原始数据进行依次累加，得到递增型序列；3得到模型的白化方程，首先，对计算

10、紧邻均值生成；然后根据GM(1,1)建模，写出灰色函数；4根据最小二乘参数估计法估计灰参数，再利用离散数据系列建立金丝的微分方程模型，得到GM(1,1)的白化方程；5白化方程的求解，得到预测值表达式，其解为时间响应函数；我们选用预测点前二十四个数据作为历史数据，运用灰色预测法进行预测。6得到灰色优化模型，通过matlab编程得出结果。5.1.2模型的建立首先，进行一次累加其次，累加数列的灰色预测模型由此，求原始数据的灰色预测模型5.1.3模型的求解预测结果:通过matlab编程对模型预测结果。得到真实观测与拟合曲线的比照用excel工具对真实值与预测值进行相关分析，得到相关系数。按照准确率计算

11、公式 时间相关系数准确率合格率四个小时内0.7878.31%87.5%将预测值与实际值进行比拟发现灰度预测模型根本能够预测未来四个小时内的风电功率趋势。5.2、问题二(对未来48小时的短期预测)5.2.1模型的分析平稳时间序列，其特性是序列的统计特性不随时间平移而变化，即均值和协方差随时间的平移而变化。运用平稳时间序列分析方法预测风电功率。而时间序列法中最常用的是累积式自回归-滑动平均模型。通常历史数据序列都是非平稳序列。因此需要进行差分变换来调整成平稳时间序列数据。建立ARMA模型的一般步骤包括模型识别、参数计算、模型检验、模型矫正和预测5.2.2模型的建立(1)数据预处理首先取预测点前两天

12、的192个样本数据序列，得到如下列图1图1序列曲线由图一说明：该样本构成的时间序列为非平稳时间序列。由此对改样本值进行有序差分变换差分算子对192个样本值进行一阶差分可得到如下序列图2图2数据差分处理序列曲线有图二说明：该序列已平稳，那么原时间序列可表示为即自回归-滑动平均模型ARIMAp,d,q(2)平稳随机时间序列模型的识别。(3)计算偏自相关系数及自相关系数。(4)自相关系数自相关函数(2.1.2)偏自相关函数当l=k时，第k系数为平稳序列的偏自相关系数运用统计学软件我们得到该平稳序列PACF/ACF (5)模型识别与模型参数估计(6)模型识别我们运用经典的Box-Jenkins模型识别

13、方法对于AR模型，其偏自相关函数满足下式对于MA模型，其相关函数满足下式对于ARMA模型，同时具备AR和MA模型的特征。由此我们给出三种模型的根本特征如下表一表一 三种根本模型特征类别模型名称 ARp MA(q) ARMA(p，q)自相关函数拖尾截尾拖尾偏自相关函数截尾拖尾拖尾有上面、统计特性，我们可以判断该模型属于ARMA(p，q)(7)模型参数估计有2，2.1我们得到该模型属于ARMA(p，q)模型，因此我们对该模型进行参数估计，由该型可推得我们令，那么这样就将原ARMAp，q变成MA(q)模型，我们根据MAq的阶数较低我们可直接求解。假设MAq的阶数较高，可运用线性迭代法求解。(8)模型

14、定阶根据AIC最小信息准那么发进行模型定阶，经过逐步的模型拟合，矩估计模型参数估计，我们最终得到当模型的阶数为AR6模型，拟合效果到达最优。由此我们确定模型为ARMA(6,1)(9)建立预测模型对于ARMA(6,1,0)模型，由于样本个数h?p，因此残差项，这里运用统计spss软件可求得预测方程为：5.2.3模型的求解运用预测方程对未来两天的数据进行预测得到如下曲线图预测数据与真实值得比拟曲线图曲线图说明根本上符合了未来两天的风电功率趋势，但预测明显存在延时性，现在对其进行误差分析：平均绝对误差：平均相对误差：时间绝对误差相对误差未来48h2.27850.0825.3问题三(风电场的选址与优化

15、布置)5.3.1风电场的选址分析风电场宏观选址的相关知识如下：宏观选址的概念：风电场宏观选址是在认真研究国家和地区风电开展规划的根底上，详细调查地区风能资源分布情侣，广泛收集区域风电场运行数据，通过对假设干场址的风能资源、电网接入和其它建设条件的分析和比拟，确定风电场的建设地点，开发价值、开发策略和开发步骤的过程，是保证公司风电产业又好又快开展的关键。宏观选址的根本原那么：1风能资源丰富、风能质量好拟选场址年平均风速一般应大于6m/s，风功率密度一般应大于200W/m2盛行风向稳定，风切变较小，端流强度较小无破坏性风速。2满足联网要求研究电网网架结构和规划开展情况，保证风电场接得进、送得出、落

16、得下。3具备交通运输和施工安装条件拟选场址周围港口、公路、铁路等交通运输条件满足风电机组、施工机械、吊装设备和其它设备、材料的进场要求。4保证工程平安拟选场址应防止洪水、潮水、地震、火灾和其它地质灾害、气象灾害等对工程造成破坏性的影响。5满足环境保护要求避开鸟类的迁徙路径、留鸟和其它动物的繁殖区。和居民区保持一定距离。尽量减少耕地等的占用。6规划装机规模满足经济型开发要求，工程满足投资回报；要求。风电场选址需要考虑的因素如下：1风能资源和气象条件相关气象条件要提出温度、沙尘、雷电、雨雪、台风对风电机组、发电量、工程施工等的影响。还应初步选择一种机型，初步比拟各选址的年发电量。2各场址的地形和交

17、通条件地形平坦单一，有利于减少湍流强度，有利于风电机组的场内运输、摆放，有利于吊装机械和其它施工机械的作业。3各场址的工程地质条件在风电场选址时，应尽量选择地震强度小、工程地质和水文地质条件较好的场址。4规划装机容量风电场装机容量取决于风电场可利用土地面积、地形、风资源分布，主导风向等因素。5接入系统初步方案场址和现有变电站的距离，是否需要新建、改建变电站，线路的电压等级，电网的结构和容量等。6环境影响初步评价风电场的建设对场址及周围环境有无影响或者影响的程度，是否涉及到自然保护区、文化遗产、风景名胜、对动植物、居民有无影响。7其他内容当地政府和居民对在该地区建风电场的态度：1用地方面：当地是

18、否已将场址规划为其他用途，或附近有无和建风电场冲突的工程。是否涉及到建筑迁拆，鱼塘、耕地、林地的占用，地下有无矿产，风电场用地是否涉及两个不同的地市、省区；2环境保护方面：是否涉及到自然保护区、文化遗产、风景名胜、对动植物、居民有无影响；3总装机容量：较大的装机容量可以摊抵道路、接入系统等固定本钱；4投资：初步做出各场址的投资估算；5电价：预计工程能够争取到的电价。5.3.2风电场优化布置1.模型的建立风电场风机优化布置的数学模型主要包括以下三个模型:1风机尾流模2风机本钱-效益模型；3风电场增量装机效益评估模型。风机尾流效应的模拟是整个风电场发电量效益的根底，因此，风机尾流模型的合理性将直接

19、影响到风电场效益估算以及风电优化布置的正确性。风电场风机本钱-效益模型用于对某一选定的风机布置方案进行其与风机相关的投资本钱核算，并结合尾流模型对评估周期内的发电量效益进行估算，该模型的合理性会影响到风机布置方案的经济性评价结论。风电场增量装机效益评价模型用于对选定的不同风机布置方案的 比照评价分析，并最终确定给定区域风电场的风机最优布置方案。2.风机尾流模型传统采用简化风机尾流线性扩张模型，即尾流影响边界距离线性增大模型。而风经过风机后的尾流影响区域直径会突然扩大。3.本钱-效益模型在风机优化布置时，只需要考虑与风机设备相关总投资本钱约占工程总投资的70%80%，包括：与风机相关固定投资本钱

20、和评估周期内的总运行维护经营本钱两项。其中，由风电场相关设备折算生产本钱模型进行和估算；与很多因素有关，如人工费、管理费和折算每年的机组大修费用等，将每年运行维护费用表示为与有关的比例函数：式中，为单位风机折算投资本钱；为风电场效益评估周期；为风电场风机总台数；为每年单位风机运行维护总本钱占的百分比系数。关于风电场效益的计算，一般主要计算其发电效率，但还要CDM清洁开展机制。在选定的风电场效益评估周期内，风电场风机总平均功率可根据品均风俗概率密度函数、风功率曲线计算得出：式中，为风电场内的平均天然风速；为受尾流效应影响下的风机工作风速。式中，风电场风速概率函数采用Weibull量参数概率分布函

21、数进行拟合求解；单台风机工作风速由尾流模型计算得到；对应风速下的风机功率由风机功率曲线插值得出。计入CDM收益的本钱折减，求出在CMD工程实施周期内的风能可利用时间：式中，、分别为风机的启动风俗和切除风速；、风别为风俗概率密度函数的Weibull分布的尺度参数和形状参数。由此可计算CMD收益的风电场风机最终折算本钱为：式中，为计入CMD收益的电量本钱折减系数。4.增量装机效益评价模型风电场风机布置的最终评价目标是：通过优化布置，使得风电场总投入与总产出之比最小，即单位投资效益最大化。评价目标函数为：采用增量评价法对风电场增装风机的本钱效益比进行评价。即在满足相同的地形、风能和风机条件下，在N台

22、风机的根底上再增装第N+1台风机后，将风机的投入和产出的增量进行比拟。如果增装的第N+1台风电机的本钱增量与产出增量之比？满足如下方程，那么说明增装该台风机经济性不好，从而即知拟建风电场布置N台风机为最优。5.计算方法对于上述的风机优化布置数学模型，可采用多种数学方法求解。采用遗传算法，按照以下步骤进行优化计算：(1)遗传编码，即对风电场流场概化，包括网格的划分、风机布点的选择等；(2)产生初始种群，即设定风电场风机数量及布置形式的初始条件；(3)适应值函数的设计，即设定表征风电场风机布置优劣的特征函数；(4)杂交算子和变异算子的设计，即如何根据已有风电场布置形式生成一种或多种新的风电场布置形

23、式；(5)控制参数的制定，即设定风电场风机布置的进化代数等控制参数，控制参数超过该设定值，计算中止；(6)计算终止条，即风电场效益评价模型是否得到满足，假设满足，那么计算完成。具体实现过程：将风电场画风为的正方形网格，假设只在每个方格的中心点布置风机，且网格边长大于风机产生的最大为刘半径。将风电场风机的每种布置方式视为一个个体，并设布置有风电机的网格点值为1，其他点为0.选择初始种群大小为n,每个个体中有N个网格点数值为1，即布置N台风机，并随机给出风电场n种风机布置形式。根据上述本钱-效益模型，建立并计算遗传算法的适应值函数？。选择适应值较小的。个个体进入交配池作为父代产生n-1个子代，与父

24、代中适应值最好的个体组成新的种群。需指出的是，为保证在个体交叉过程中的风机数目不变，算法从一个父代中两次随机选取k个基因，两次选取的基因各异，将两次选取的k个基因交叉得出新的个体。设置风电场风机优化配置的最大进化代数控制值，当进化代数超过此控制值时那么中断。5.3.3配置参考通过查资料获得我国风电场分布图如下：最终我们决定选取福建东山风电场作为参考结果，结果如下表：地区台数容量kw估计年利用小时数福建东山1060002800六、模型的评价与优化6.1模型一的评价优点：算法简单、可利用较少数据建模。缺点：对于复杂的非线性系统来说预测的效果不是很理想6.2模型二的评价优点：对于短、近期预测比拟显著

25、。缺点：延伸到更远的未来就会出现很大的局限性，导致预测值偏离实际较大而使决策出错。6.3模型三的评价优点：对于广域范围内进行选址效果明显。缺点：当地域具有一定局限性时无法精确选址。七、模型的推广参考文献：1 姜启源. 数学模型第三版M.北京:高等教育出版社，1999.2 韩中庚.数学建模方法及其运用第二版M.北京：高等教育出版社,20213 丁明.给予时间序列分析的风电场风速预测模型J电气自动化设备,2002,18(4)60-724 赵建华.风电厂选址及最正确风力发电机组选型J.电力勘测设计,2021,21(1):56-58. 5 赵伟然.风电场选址与风机优化排布实用技术探讨J.电力科学与工程

26、,2021,03(4):44-476 张新房.风力发电技术的开展及假设干问题J.现代电力,2003,05(6):33-37. 7 陶奕衫.风电场宏观选址综合决策方法的研究J.四川电力技术,2021,12(1):21-23. 附录附录一：Matlab灰度预测程序：程序1：%非线性最小二乘法拟合%灰色模型的风功率预测 Grey modelclc,clear all,close all%预测未来四小时的功率t=1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24;x0=41.07 38.47 38.67 36.93 38.93

27、40.67 42.07 38.27 37.67 38.2 38.4 . 37.73 35.27 34.07 32.67 29.4 30.07 26.53 26.13 31.53 33.2 31.6 . 28.87 25.8;n=length(x0);pre_num=n;disp(级比检验)lambda=x0(1:end-1)./x0(2:end)range=minmax(lambda);x1=cumsum(x0);z=0.5*(x1(2:end)+x1(1:end-1);Y=x0(2:end);B=-z(1:end) ones(n-1,1);u=inv(B*B)*B*Y; a=u(1);b=u

28、(2);for k=1:pre_num-1 x0_pre(k+1)=(x0(1)-b/a)*(exp(-a*k)-exp(-a*(k-1);end err=x0 - x0_pre(1:n); epsilon=abs(err)./x0(1:n).*100; disp(预测值); disp(x0_pre); disp(相对误差); disp(epsilon); t=1:24; plot(t,x0,rp,t,x0_pre,g-) %原始数据与预测数据的比拟xlabel(时间);ylabel(风功率(KW); legend(原始数据,预测数据,2)%预测数据表t1=89:112; %预测值后四小时n1

29、=length(t1); pre_num1=n1;for k=1:pre_num1-1 x0_pre1(k+1)=(x0(1)-b/a)*(exp(-a*k)-exp(-a*(k-1);end figure n=97 - 89; i=0; for k=n+1:pre_num1 i=i+1; t2(i)=t1(1,k); x0_pre2(i)=x0_pre1(1,k); end plot(t2,x0_pre2,r-) %预测数据表xlabel(时间);ylabel(功率(KW);legend(预测数据);程序2：%非线性最小二乘法拟合%灰色模型的风功率预测 Grey modelclc,clear

30、 all,close all %预测未来四小时的功率t=74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97;x0=43.8 45 49.6 55.4 57.33 58.67 53.33 61 63.07 59.33. 58.07 57.73 64.93 76.2 76.4 77.6 71.8 73.27 82.2 83.07 80.47. 82.8 79 76.13;n=length(x0);pre_num=n;disp(级比检验)lambda=x0(1:end-1)./x0(2:end)range=m

31、inmax(lambda);x1=cumsum(x0);z=0.5*(x1(2:end)+x1(1:end-1);Y=x0(2:end);B=-z(1:end) ones(n-1,1);u=inv(B*B)*B*Y; a=u(1);b=u(2);for k=1:pre_num-1 x0_pre(k+1)=(x0(1)-b/a)*(exp(-a*k)-exp(-a*(k-1);end err=x0 - x0_pre(1:n); epsilon=abs(err)./x0(1:n).*100; disp(预测值); disp(x0_pre); disp(相对误差); disp(epsilon); t=74:97; plot(t,x0,o,t,x0_pre,-) %原始数据与预测数据的比拟 xlabel(时间);ylabel(风功率(KW); legend(原始数据,预测数据,2) %预测数据表t1=89:112; %预测值1四小时n1=length(t1); pre_num1=n1;for k=1:pre_num1-1 x0_pre1(k+1)=(x0(1)-b/a)*(exp(-a*k)-exp(-a*(k-1);end figure n=97 - 89; i=