CN107220449A - 一种基于分段非线性函数的风机功率曲线建模方法 - Google Patents

一种基于分段非线性函数的风机功率曲线建模方法 Download PDF

Info

Publication number
CN107220449A
CN107220449A CN201710425152.XA CN201710425152A CN107220449A CN 107220449 A CN107220449 A CN 107220449A CN 201710425152 A CN201710425152 A CN 201710425152A CN 107220449 A CN107220449 A CN 107220449A
Authority
CN
China
Prior art keywords
power
data
piecewise
fitting
function
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
Application number
CN201710425152.XA
Other languages
English (en)
Inventor
魏善碧
柴毅
李枫
廖瑞勇
冯川
何馨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing University
Original Assignee
Chongqing University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chongqing University filed Critical Chongqing University
Priority to CN201710425152.XA priority Critical patent/CN107220449A/zh
Publication of CN107220449A publication Critical patent/CN107220449A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Structures Of Non-Positive Displacement Pumps (AREA)

Abstract

本发明公开了一种基于分段非线性函数的风机功率曲线建模方法,属于风机功率曲线分析和建模领域。包括以下步骤:首先确定获取实际数据的时间尺度。通过循环控制将数据进行预处理,接着进行数据清洗,确立所分段数,使用各拟合函数形成混合整数规划模型,通过计算机自动优化获得分界点形成分段函数,在这同时对清洗后的数据进行分段拟合。经过分析验证后,若误差较大,可调整段数和模型再次拟合,若满足要求,则显示结果,同时保存方便下次使用。本发明采用混合整数规划自动得到分界点,并对实际数据进行分段拟合,能有效减少低风速段和高风速段的直线数据段拟合带来的误差,简化了风机功率曲线拟合的步骤,提高风机功率曲线的精度和可靠性。

Description

一种基于分段非线性函数的风机功率曲线建模方法
技术领域
本发明属于风机功率曲线分析和建模领域,涉及一种基于风机功率曲线实际数据的,将混合整数规划引入以实现自动分段函数拟合功率曲线降低误差的方法。
背景技术
风力发电在提倡新能源和清洁能源的二十一世纪,正在逐渐被各国政府和企业所重视。而在风力发电的领域内,风力发电机组的功率曲线作用举足轻重。功率曲线描述风速与机组输出功率之间关系,其不仅是风机设计过程中进行控制系统设计时的重要依据,更是考核风力发电机组性能、评估风力发电机组发电能力的一项重要指标。风机功率曲线对风机寿命、风机机组发电量以及风机性能检测等都有着举足轻重的作用。现有风电功率曲线的机理模型,但机理模型往往难以运用。于是现如今大多数功率模型是通过对风机实际数据进行模型拟合。
但是,基于实际数据拟合模型也面临一些问题:
(1)现如今常见的精度较高的四参数、五参数和七参数拟合方法在功率曲线的起始段和末尾段拟合效果较差,实际数据中,起始段有一段功率为零的数据带,而高风速的末尾段则会到达功率最大值呈现一段直线带,一头一尾的数据影响着中间数据的拟合;
(2)为了更好的拟合效果,通常在取线时将头尾去除,当作分段函数处理,此时需要提前增加程序步骤,同时人为的操作不能准确地找出分段点,影响曲线精度;
分段函数是一个函数在不同自变量定义的时候拥有不同的对应法则,在风机功率曲线上,低风速段的零点段和高风速段的满发功率直线段与中间部分的“S”形曲线在分界点处不连续,所以不存在只一个普通函数式就可以完美地拟合出一条完整的功率曲线,应用分段函数,有效地降低拟合带来的误差,不仅使中间部分的拟合曲线更贴近原效果,同时减少了在分界点处不连续带来的误差。一般的拟合函数如logistic函数对水平直线的拟合效果较差,采用分段函数可以有效地解决这个问题。
使用分段函数会产生上述的第二个问题,现今暂时解决去除头尾的问题有两个方法,一是通过人为地设定分界点,即通过对图形的大体印象,设定一个不变的常数作为两个分界点,但除了人为设定不精确之外,不变的常数不能适用于不同的风机,若不同的风机设定不同的常数又增加了工作量;而另一种方法是先采用类似四参数logistic方程进行拟合,将拟合出的曲线穿过横坐标轴和功率最大值的点设为分界点,但此时不仅增加了拟合的次数,同时这种分界点容易受到拟合方程原函数的约束,准确性下降。本发明采用整数规划来解决以上问题。整数规划是指一类要求问题中的全部或一部分变量为整数的数学规划。是近三十年来发展起来的、规划论的一个分支。整数规划问题是要求决策变量取整数值的线性规划或非线性规划问题。在整数规划中,如果所有变量都限制为整数,则称为纯整数规划;如果仅一部分变量限制为整数,则称为混合整数规划。整数规划的一种特殊情形是01规划,它的变数仅限于0或1。
发明内容
有鉴于此,本发明的目的在于提供一种将混合整数规划和分段函数应用于大量风机运行实际数据相结合的模型计算技术,该方法采用了分段函数,通过设定两个分界点将风机功率曲线分为三段,有效地降低了风机功率曲线由于头尾的直线段对最后拟合效果的影响,同时直接采用风机的实际数据拟合,降低了由于实际环境、不同风机对风机功率曲线的影响。并利用混合整数规划,通过0-1规划将分段函数当做一个整体进行拟合,通过计算机自动优化寻找分界点,进一步提高了整个方法的精度和可靠性。
基于分段非线性函数的风机功率曲线模型,包括以下模块:数据采集模块,数据预处理模块;数据清洗模块;混合整数规划模块;分段拟合模块;数据分析处理模块。
基于分段非线性函数的风机功率曲线建模技术,它包含以下步骤:
步骤一:将数据按所需时间尺度进行读入,并通过简单的循环处理去除丢失参数和严重故障的数据,再将预处理后的数据进行数据清洗,清洗方法可采用OIV算法或者多倍均方差算法;
步骤二:设定所需的分段数,建立混合整数规划模型,将分界点当作变量列入模型之内,同时将零值、中间段拟合方程和最大值等加入到模型内,对同一台风机的大量数据使用整个模型进行自定义拟合,通过计算机自动优化寻找到分界点形成分段函数,同时每一段都进行非线性拟合,中间段拟合采用四参数logistic方程进行非线性拟合;
步骤三:对非线性拟合的风机功率曲线进行误差检验,通过拟合的函数方程计算如拐点等所需的数据,通过屏幕输出结果,同时保存拟合的结果。
本发明的有益技术效果为:本发明可以充分利用分段函数的特性,将不同的数据段用不同的函数进行拟合,有效地降低了不同风速段的数据集采用相同模型拟合带来的误差,同时采用混合整数规划将整个模型放入大量的实际数据进行拟合,能够一次性得到所有的分界点,降低了人为确认分段函数分界点的不确定性,相对于传统的使用实际数据进行取线拟合的方法,节省了采用bin法、最大似然法或者最大值法取线的步骤,同时也简化了数据后期处理分析的步骤,其任意分段的可扩展性也提高了该模型在不同风机的功率曲线建模下的适用性和可靠性。
附图说明
为了使本发明的目的、技术方案和有益效果更加清楚,本发明提供如下附图进行说明:
图1为本发明所述基于分段非线性函数的风机功率曲线建模技术的系统结构图。
图2为本发明所述基于分段非线性函数的风机功率曲线建模技术的系统模块图。
图3为本发明所述基于分段非线性函数的风机功率曲线建模技术的系统流程图。
具体实施方式
下面将结合附图,对本发明的优选实施例进行详细的描述。
图2为本发明所述基于分段非线性函数的风机功率曲线建模技术的系统结构图,包括以下模块:数据采集模块,数据预处理模块;数据清洗模块;混合整数规划模块;分段拟合模块;数据分析处理模块。数据采集模块根据预先设定的时间尺度将数据读取,数据预处理模块将大量选取的数据进行简单的处理例如去除负数点等,随后数据清洗模块对预处理后的数据进行清洗,同时混合整数规划模块根据所选的段数构建整数规划模型。随后分段拟合模块将混合整数模块自动优化得到的分界点进行分段拟合,数据分析处理模块对拟合结果进行误差分析,若误差较大,则重新选取段数拟合,若误差符合精度要求,则进行数据的后期分析处理、显示以及存储。
图3为本发明所述基于分段非线性函数的风机功率曲线建模技术的系统流程图,具体实现步骤如下:
(1)通过数据采集模块,设定风机实际运行数据的时间尺度(1分钟平均、5分钟平均、10分钟平均),将大量实际数据读取到所设定的变量之中。
(2)数据预处理模块将一定时间尺度的数据读取进来后,通过设定的简单的逻辑循环结构,将数据分成多个风速比恩,对每一个比恩的数据进行数据筛选,将类似高风速零功率点的故障点以及与在同一比恩内平均风速误差过大的湍流点等先行去除。随后将处理过的数据保存在预定的变量或者文件之中。
(3)数据清洗模块对预处理后的数据进行清洗,采用改良后的滑差算法对预处理模块无法清洗的边缘点进行清洗,清洗过程为将数据等比例分为多个比恩,对每一个比恩内的数据进行排序,依次计算前i个坐标点中功率的方差值,并作为该点对应的滑差值,由于这种计算方式类似于滑动,故将计算结果形象地称为滑差值,所得到的数据集称为滑差集;最后,通过与初始阈值S比较分析实现数据清洗的功能。清洗过后,对数据进行检测,若故障点已基本消除,则继续运行,若仍存在大量故障点,则需重复清洗直到满足要求。
(4)数据清洗完成后,确定分段函数所分的段数,下面的步骤以分3段为例。
将曲线分为J1、J2、J3三段,同时设两个分界点的横坐标为F1、F2(F1<F2),其中F1是风机功率曲线的切入点,而F2是风机功率曲线的切出点,x是自变量风速,同时设两个状态量δ1、δ2
则令:
则有非线性规划模型:
J=δ1δ2J1+(1-δ12J2+(1-δ1)(1-δ2)J3
(5)分段拟合模块将预先设定的拟合函数代入进行拟合。
其中起始段J1和末尾段J3的拟合方程如下:
J1=||0-y||
J3=||pmax-y||
其中y为实际数据点的功率值,pmax是实际数据中功率的最大值,可以设pmax=f(F2),即最大功率值为切出点处的功率值。
中间段J2采用四参数的logistic方程,其式子为:
其中p、m、n、r为logistic方程的四个待定参数,x是自变量风速。
所以:
J2=||f(x)-y||
将这些函数代入混合整数规划模型,使用自定义非线性拟合函数fittype对实际数据进行拟合,将优化计算得到的结果保存到设定的变量中。
(6)最后,数据分析处理模块将拟合的结果与原实际数据进行误差分析,若误差不满足实际要求,则重新设置段数或者重新设置不同的拟合函数进行拟合,直到误差达到系统的精度要求为止。
当误差达到所需精度后,即拟合效果满足要求,将处理过后的数据显示在屏幕上,并存储到预先设定好的变量或者文件之中,等待之后的处理。
最后说明的是,以上优选实施例仅用以说明本发明的技术方案而非限制,尽管通过上述优选实施例已经对本发明进行了详细的描述,但本领域技术人员应当理解,可以在形式上和细节上对其作出各种各样的改变,而不偏离本发明权利要求书所限定的范围。

Claims (4)

1.一种基于分段非线性函数的风机功率曲线建模技术,其特征在于:包括以下数据处理步骤:
步骤一:将数据按所需时间尺度进行读入,并通过循环处理去除丢失参数和严重故障的数据,再将预处理后的数据进行数据清洗,设定所需的分段数;
步骤二:建立混合整数规划模型,将分界点当作变量列入模型之内,通过计算机自动优化寻找到分界点形成分段函数,同时使用不同函数对每一段进行非线性拟合;
步骤三:对非线性拟合的风机功率曲线进行误差检验,通过拟合的函数方程计算如拐点等所需的数据,同时保存拟合的结果。
2.根据权利要求1所述的基于分段非线性函数的风机功率曲线建模技术,其特征在于:步骤一中所述的设定分段数,本次建模技术可以将功率曲线模型分为任意段,可以从不同段数的分段模型去分析系统模型,对不同的风机在不同的环境以及所要求的功率曲线精度不同,不同的段数使功率曲线模型适用性更强。
3.根据权利要求1所述的基于分段非线性函数的风机功率曲线建模技术,其特征在于:步骤二所述的混合整数规划模型,将分段函数不同定义段的不同函数形成一个整体,将整体模型进行拟合简化了分段函数拟合的步骤,同时也省去了原步骤中的比恩取线步骤。
4.根据权利要求1所述的基于分段非线性函数的风机功率曲线建模技术,其特征在于:步骤二所述的混合整数规划模型,通过计算机自动优化得到分段函数分界点,降低了人为直接设定或先拟合再定点带来的影响,简化了原风机功率曲线后期数据处理的步骤。
CN201710425152.XA 2017-06-07 2017-06-07 一种基于分段非线性函数的风机功率曲线建模方法 Pending CN107220449A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710425152.XA CN107220449A (zh) 2017-06-07 2017-06-07 一种基于分段非线性函数的风机功率曲线建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710425152.XA CN107220449A (zh) 2017-06-07 2017-06-07 一种基于分段非线性函数的风机功率曲线建模方法

Publications (1)

Publication Number Publication Date
CN107220449A true CN107220449A (zh) 2017-09-29

Family

ID=59947796

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710425152.XA Pending CN107220449A (zh) 2017-06-07 2017-06-07 一种基于分段非线性函数的风机功率曲线建模方法

Country Status (1)

Country Link
CN (1) CN107220449A (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111523707A (zh) * 2020-04-08 2020-08-11 江苏方天电力技术有限公司 基于分段数据拟合的综合能源系统设备变工况特性建模方法
CN112231182A (zh) * 2020-12-15 2021-01-15 长沙树根互联技术有限公司 物联设备工况数据分析方法、装置和计算机设备
CN112462312A (zh) * 2020-11-04 2021-03-09 胜达克半导体科技(上海)有限公司 一种用于芯片测试机的自动校准方法及其应用
CN113064077A (zh) * 2021-02-20 2021-07-02 山东科技大学 一种基于分段拟合的锂电池寿命衰减曲线拟合方法
CN114472765A (zh) * 2021-11-23 2022-05-13 重庆大学 一种环件轧制过程环件圆度误差理论计算的方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103020377A (zh) * 2012-12-26 2013-04-03 国电联合动力技术有限公司 一种直驱风力发电机的螺栓强度分析方法
CN104573172A (zh) * 2014-11-19 2015-04-29 新疆金风科技股份有限公司 风力发电机组中的结构件的疲劳分析方法和疲劳分析装置
CN106447210A (zh) * 2016-10-10 2017-02-22 国家电网公司 一种计及可信度评价的配网设备健康度动态诊断方法
CN106529075A (zh) * 2016-11-28 2017-03-22 西安交通大学 一种考虑分时段的非线性模拟风速方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103020377A (zh) * 2012-12-26 2013-04-03 国电联合动力技术有限公司 一种直驱风力发电机的螺栓强度分析方法
CN104573172A (zh) * 2014-11-19 2015-04-29 新疆金风科技股份有限公司 风力发电机组中的结构件的疲劳分析方法和疲劳分析装置
CN106447210A (zh) * 2016-10-10 2017-02-22 国家电网公司 一种计及可信度评价的配网设备健康度动态诊断方法
CN106529075A (zh) * 2016-11-28 2017-03-22 西安交通大学 一种考虑分时段的非线性模拟风速方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
郭晓杰: ""监督学习方法在短期风电功率预测中的应用研究"", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111523707A (zh) * 2020-04-08 2020-08-11 江苏方天电力技术有限公司 基于分段数据拟合的综合能源系统设备变工况特性建模方法
CN111523707B (zh) * 2020-04-08 2022-06-07 江苏方天电力技术有限公司 基于分段数据拟合的综合能源系统设备变工况特性建模方法
CN112462312A (zh) * 2020-11-04 2021-03-09 胜达克半导体科技(上海)有限公司 一种用于芯片测试机的自动校准方法及其应用
CN112231182A (zh) * 2020-12-15 2021-01-15 长沙树根互联技术有限公司 物联设备工况数据分析方法、装置和计算机设备
CN113064077A (zh) * 2021-02-20 2021-07-02 山东科技大学 一种基于分段拟合的锂电池寿命衰减曲线拟合方法
CN114472765A (zh) * 2021-11-23 2022-05-13 重庆大学 一种环件轧制过程环件圆度误差理论计算的方法

Similar Documents

Publication Publication Date Title
CN107220449A (zh) 一种基于分段非线性函数的风机功率曲线建模方法
Yang et al. Netadapt: Platform-aware neural network adaptation for mobile applications
CN109543559A (zh) 基于孪生网络和动作选择机制的目标跟踪方法及系统
CN104794492B (zh) 基于功率特征模型的机床设备加工运行状态在线识别方法
CN103838820B (zh) 基于近邻传播的进化多目标优化社区检测方法
CN112131206B (zh) 一种多模型数据库OrientDB参数配置自动调优方法
CN110674120A (zh) 一种风电场数据清洗方法及装置
CN110727670B (zh) 基于流程图的数据结构预测传递及自动化数据处理方法
CN115438714A (zh) 分簇式联邦学习驱动的机械故障诊断方法、装置及介质
CN112733455A (zh) 水轮机建模数据的修正方法、系统、装置和存储介质
CN111160170B (zh) 一种自学习的人体行为识别与异常检测方法
CN110348540B (zh) 基于聚类的电力系统暂态功角稳定故障筛选方法及装置
CN111209914A (zh) 一种异常风功率数据的剔除方法
CN108763584B (zh) 一种风功率曲线散点过滤的方法及其系统
CN108198268A (zh) 一种生产设备数据标定方法
CN113537245B (zh) 一种基于特征图的神经网络剪枝方法
CN110210657B (zh) 一种基于单机模型的风机功率预测方法、系统及计算机存储介质
CN109035311A (zh) 一种弯骨骨折自动配准及内固定钢板预弯建模方法
CN114871858B (zh) 一种机器人全口径研磨抛光系统及方法
CN116662804A (zh) 基于数据驱动的烘丝机烘丝质量预测模型建立方法及装置
CN110517226A (zh) 基于双边滤波的多特征纹理图像融合的烟梗区域提取方法
CN104134091B (zh) 一种神经网络训练方法
CN108573275A (zh) 一种在线分类微服务的构建方法
CN114070438A (zh) 一种5g基站电磁信号的检测方法及系统
CN114077924A (zh) 考虑极端大风天气的风电功率预测大风数据增强方法

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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20170929