CN116882079A - Water pump characteristic curve self-adaptive calibration and prediction method - Google Patents

Water pump characteristic curve self-adaptive calibration and prediction method Download PDF

Info

Publication number
CN116882079A
CN116882079A CN202310697387.XA CN202310697387A CN116882079A CN 116882079 A CN116882079 A CN 116882079A CN 202310697387 A CN202310697387 A CN 202310697387A CN 116882079 A CN116882079 A CN 116882079A
Authority
CN
China
Prior art keywords
water pump
data
collected
characteristic
flow
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
CN202310697387.XA
Other languages
Chinese (zh)
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.)
Zhejiang Yuanchuang Intelligent Control Technology Co ltd
Zhejiang University ZJU
Original Assignee
Zhejiang Yuanchuang Intelligent Control Technology Co ltd
Zhejiang University ZJU
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 Zhejiang Yuanchuang Intelligent Control Technology Co ltd, Zhejiang University ZJU filed Critical Zhejiang Yuanchuang Intelligent Control Technology Co ltd
Priority to CN202310697387.XA priority Critical patent/CN116882079A/en
Publication of CN116882079A publication Critical patent/CN116882079A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/27Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/12Timing analysis or timing optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Control Of Positive-Displacement Pumps (AREA)

Abstract

本发明公开了一种水泵特性曲线自适应校准及预测方法,包括以下步骤:经外部设备实时采集水泵各个特征的数据,将采集到的数据上传到指定平台,并对当前和历史数据储存;建立基于分类改进的ARMA模型,并对水泵各个特征进行预测;将步骤S2中水泵各个特征的预测值以及当前时刻水泵各个特征的采集值,代入水泵在额定频率下的运行规律公式,对当前时刻水泵的特性曲线进行修正和预测;基于采集的管路参数通过公式以最小二乘的方式进行拟合;绘制管路特性曲线和水泵特性曲线,该两条曲线的交点就是预测水泵工作的实际工况点。本发明具有能够有效提高准确性以及扩大适用范围的特点。

The invention discloses an adaptive calibration and prediction method for water pump characteristic curves, which includes the following steps: collecting data on various characteristics of the water pump in real time through external equipment, uploading the collected data to a designated platform, and storing current and historical data; establishing Based on the classification-improved ARMA model, each feature of the water pump is predicted; the predicted values of each feature of the water pump in step S2 and the collected values of each feature of the water pump at the current moment are substituted into the operating rule formula of the water pump at the rated frequency to predict the water pump at the current moment. Correct and predict the characteristic curve; fit the collected pipeline parameters in a least squares manner through formulas; draw the pipeline characteristic curve and the water pump characteristic curve. The intersection of the two curves is to predict the actual working conditions of the water pump. point. The invention has the characteristics of effectively improving accuracy and expanding the scope of application.

Description

一种水泵特性曲线自适应校准及预测方法An adaptive calibration and prediction method for water pump characteristic curves

技术领域Technical field

本发明涉及一种水泵特性处理方法,尤其涉及一种水泵特性曲线自适应校准及预测方法。The invention relates to a water pump characteristic processing method, and in particular to a water pump characteristic curve adaptive calibration and prediction method.

背景技术Background technique

当前传统的水泵建模方式通常以机理建模为主,在模型中没有考虑同一类型的水泵在不同工况环境以及随着使用时间的增长而带来的特性曲线的变化,因此传统的模型会有较大的误差。此外,由于不同工业场景的特殊性,采集到的数据会受到一定的限制(特殊的工业场景下设备无法任意的启停),只能采集到定频或者特定频率下的部分数据点。例如,在专利名称为“变频调速水泵数学模型的建模方法和系统”的现有技术中,其基于海量的水泵历史运行数据,对其按照相同的频率段进行分组,拟合对应频率段的流量-扬程和流量-效率曲线。The current traditional water pump modeling methods are usually based on mechanism modeling. The changes in the characteristic curves of the same type of water pump in different working conditions and with the increase of use time are not considered in the model. Therefore, the traditional model will There is a large error. In addition, due to the particularity of different industrial scenarios, the collected data will be subject to certain restrictions (the equipment cannot be started and stopped arbitrarily in special industrial scenarios), and only some data points at fixed frequencies or specific frequencies can be collected. For example, in the existing technology patent titled "Modeling Method and System for Mathematical Model of Variable Frequency Adjustable Speed Water Pumps", it is based on a large amount of historical water pump operating data, grouping them according to the same frequency band, and fitting the corresponding frequency band The flow-head and flow-efficiency curves.

但仍存在以下几点不足:However, there are still several shortcomings:

1.该模型族仅适用于对同一频率带内的数据拟合水泵的流量-扬程和流量-效率曲线,当采集数据不能满足在同一频率带或在同一频率带内的数据点较少时会受到限制。1. This model family is only suitable for fitting the flow-head and flow-efficiency curves of water pumps to data in the same frequency band. It will be used when the collected data cannot meet the requirements of the same frequency band or there are few data points in the same frequency band. restricted.

2.该模型的基于历史数据对水泵特性曲线的拟合,可以一定程度上应对水泵长期使用而带来特性曲线偏移的问题,但是在水泵发生故障,采集数据偏差较大时无法及时发现。2. The model's fitting of the water pump characteristic curve based on historical data can, to a certain extent, deal with the problem of characteristic curve deviation caused by long-term use of the water pump. However, when the water pump fails and the collected data deviates greatly, it cannot be discovered in time.

因此,现有的技术存在着准确性较差以及适用范围较窄的问题。Therefore, the existing technology has problems of poor accuracy and narrow scope of application.

发明内容Contents of the invention

本发明的目的在于,针对现有技术的不足,提供一种水泵特性曲线自适应校准及预测方法。本发明具有能够有效提高准确性以及扩大适用范围的特点。The purpose of the present invention is to provide an adaptive calibration and prediction method for water pump characteristic curves in view of the shortcomings of the existing technology. The invention has the characteristics of effectively improving accuracy and expanding the scope of application.

本发明的技术方案:一种水泵特性曲线自适应校准及预测方法,包括以下步骤:The technical solution of the present invention: an adaptive calibration and prediction method for water pump characteristic curves, including the following steps:

S1、经外部设备实时采集水泵和管路两部分各个特征的数据,将采集到的数据上传到指定平台,并对当前和历史数据储存;S1. Collect data on various characteristics of the water pump and pipeline in real time through external equipment, upload the collected data to the designated platform, and store current and historical data;

S2、建立基于分类改进的ARMA模型,并对水泵各个特征进行预测;S2. Establish an ARMA model based on classification improvement and predict various characteristics of the water pump;

S3、将步骤S2中水泵各个特征的预测值以及当前时刻水泵各个特征的采集值,代入水泵在额定频率下的运行规律公式包括流量-扬程曲线方程、流量-功率曲线方程以及效率曲线方程,对当前时刻水泵的特性曲线进行修正,并得到预测的下一个时间单位内水泵运行时的流量-扬程和流量-效率曲线;S3. Substitute the predicted values of each characteristic of the water pump in step S2 and the collected values of each characteristic of the water pump at the current moment into the operating law formula of the water pump at the rated frequency, including the flow-head curve equation, the flow-power curve equation and the efficiency curve equation. The characteristic curve of the water pump at the current moment is corrected, and the predicted flow-head and flow-efficiency curves when the water pump is running in the next time unit are obtained;

S4、基于采集的管路参数以最小二乘的方式拟合得到管路特性曲线,拟合公式为其中/>为管路的扬程,Q为管路的流量,D,S为需要拟合的参数;绘制管路特性曲线和水泵特性曲线,两条曲线的交点为预测水泵工作的实际工况点。S4. Based on the collected pipeline parameters, the pipeline characteristic curve is obtained by fitting the least squares method. The fitting formula is: Among them/> is the head of the pipeline, Q is the flow rate of the pipeline, D and S are the parameters that need to be fitted; draw the pipeline characteristic curve and the water pump characteristic curve, and the intersection of the two curves is the actual working condition point for predicting the work of the water pump.

进一步地,步骤S1中,每日采集一次水泵和管路各个特征的数据,每条数据对应水泵和管路的流量、扬程;Further, in step S1, data on various characteristics of the water pump and pipeline are collected once a day, and each piece of data corresponds to the flow rate and head of the water pump and pipeline;

对于每一个采集的数据,通过多次采集取平均值的方式,作为该条数据的测量值。For each collected data, the average value is taken through multiple collections as the measured value of this piece of data.

进一步地,步骤S2中,基于分类改进的ARMA模型通过以下步骤进行实现:Further, in step S2, the ARMA model based on classification improvement is implemented through the following steps:

1)对步骤S1中采集的水泵各个特征的数据进行分类,分为同一频率下的数据即定频数据和不同频率下的数据即变频数据;1) Classify the data of each characteristic of the water pump collected in step S1 into data at the same frequency, that is, fixed frequency data, and data at different frequencies, that is, variable frequency data;

2)对定频数据和变频数据分别按季度分类,并提取同一类中的各个特征数据信息,包括各个特征的最大值、最小值、均值、截尾均值和中位数,选取特征的其中一个数据信息来描述该类的该特征,记为表示第i个类下的第j个特征点采集数据的中位数;2) Classify fixed-frequency data and variable-frequency data by quarter, and extract each feature data information in the same category, including the maximum value, minimum value, mean, censored mean and median of each feature, and select one of the features Data information is used to describe the characteristics of this class, recorded as Represents the median of the j-th feature point collection data under the i-th class;

3)将每一类中各个特征的数据信息构成向量,根据Legendre最佳逼近准则,找到与特征数据信息向量最接近的采集数据向量,并记录下对应的本类中的日期di,t3) Construct the data information of each feature in each category into a vector, and according to the Legendre best approximation criterion, find the collection data vector closest to the feature data information vector, and record the corresponding date d i,t in this category;

4)在分类后的数据集中取前N个类下的m个特征,基于ARMA模型对当前类即第N+1类的这m个特征点分别进行预测,具体公式如下:4) Take m features from the first N categories in the classified data set, and predict the m feature points of the current class, that is, the N+1th class, based on the ARMA model. The specific formula is as follows:

其中,p为自回归项AR的阶数,αi(i=0,1,K,p)自回归系数;q为移动平均项MA的阶数,βj(j=1,2,K,q)为移动平均系数,ξt为白噪声,限定了噪声的边界;E表示数学期望,Var表示方差,/>表示任意采集的同一类型的水泵数据。Among them, p is the order of the autoregressive term AR, α i (i=0,1,K,p) autoregressive coefficient; q is the order of the moving average term MA, β j (j=1,2,K, q) is the moving average coefficient, ξ t is white noise, Defines the boundary of noise; E represents mathematical expectation, Var represents variance,/> Represents any collected water pump data of the same type.

5)用采集得到的各个特征变量的值代入上述步骤中得到预测的下一个类的各个特征的数据信息,对应的日期为 5) Substitute the collected values of each feature variable into the above steps to obtain the predicted data information of each feature of the next class. The corresponding date is

6)基于测量的第N个类中的最后kt日的特征值和预测的未来的特征数据信息向量,拟合各个特征的变化曲线,用函数fi来表述,以15日为周期进行预测,实现对下一季度每15日进行一次水泵特性曲线的更新和迭代;具体公式如下:6) Based on the measured characteristic values of the last kt day in the Nth class and the predicted future characteristic data information vector, fit the change curve of each characteristic, expressed by the function f i , and predict with a period of 15 days, Realize the update and iteration of the water pump characteristic curve every 15 days in the next quarter; the specific formula is as follows:

用最小二乘的方式求解公式,得到函数fi,其中,函数fi表示第i个特征拟合的预测函数,基于预测函数得到下一个类中任意时刻的预测值。Use the least squares method to solve the formula to obtain the function f i , where the function f i represents the prediction function of the i-th feature fitting, and the predicted value at any time in the next class is obtained based on the prediction function.

进一步地,步骤S3中,流量-扬程曲线方程、流量-功率曲线方程以及效率曲线方程,分别如下所示:Further, in step S3, the flow-head curve equation, flow-power curve equation and efficiency curve equation are as follows:

H=aQ2+bQ+cH=aQ 2 +bQ+c

P=jQ2+kQ+lP=jQ 2 +kQ+l

其中,H是水泵扬程,Q是水泵流量,η是水泵效率,P是水泵功率,a,b,c,j,k,l是拟合二次曲线的参数。Among them, H is the water pump head, Q is the water pump flow rate, eta is the water pump efficiency, P is the water pump power, a, b, c, j, k, l are the parameters of fitting the quadratic curve.

进一步地,同一频率下采集的数据,通过流量-扬程曲线方程、流量-功率曲线方程能够得到特定频率下的流量-扬程曲线方程和流量-功率曲线方程;Furthermore, for the data collected at the same frequency, the flow-head curve equation and the flow-power curve equation at a specific frequency can be obtained through the flow-head curve equation and the flow-power curve equation;

不同频率下采集的数据,通过对采集数据以最小二乘的处理方式结合机理模型进行拟合,如下所示:The data collected at different frequencies are fitted by combining the least squares processing method with the mechanism model, as shown below:

H0=aQ0 2+bQ0i+ci2 H 0 =aQ 0 2 +bQ 0 i+ci 2

P0=jQ0 2i+kQ0i2+li3 P 0 =jQ 0 2 i+kQ 0 i 2 +li 3

其中,H0,Q0,P0为在水泵运行频率为f0时的扬程、流量以及功率,f为额定频率,i为运行频率与额定频率的比例系数。Among them, H 0 , Q 0 , P 0 are the head, flow rate and power when the pump operating frequency is f 0 , f is the rated frequency, and i is the proportional coefficient between the operating frequency and the rated frequency.

本发明的有益效果:Beneficial effects of the present invention:

与现有技术相比,本发明为了提升水泵的数学模型的准确性,对于特定频率下的水泵的特性曲线采用两方面分别建模进行数据融合:一是变频调速下的水泵特性曲线拟合,采集不同频率下的水泵的流量、扬程、功率的参数,通过机理模型将这些数据转化为同一频率下的数据,进一步给出该对应频率的特性曲线;二是定频时的水泵特性曲线的拟合,通过采集的是同一频率下的水泵特性参数,并拟合出该频率下的特性曲线;最后是管路特性曲线的拟合,通过测量管路在不同流量下时的管路扬程,拟合该场景下的管路特性曲线。本发明将数据采集、数据处理、数据检验、水泵特性曲线的实时更新集成化、自动化,不再需要研发人员进行额外的工作;通过采集的同一频率和不同频率下的数据能够预测未来一个或多个时间单位内水泵的流量-扬程和流量效率曲线,该方法降低了工作复杂度的同时保持了特性曲线的精度。而且,应用场景可以推广,并不只适用于某一特定的工业场景。Compared with the existing technology, in order to improve the accuracy of the mathematical model of the water pump, the present invention adopts two aspects of modeling for data fusion for the characteristic curve of the water pump at a specific frequency: First, the fitting of the water pump characteristic curve under variable frequency speed regulation , collect the flow, head, and power parameters of the water pump at different frequencies, convert these data into data at the same frequency through the mechanism model, and further give the characteristic curve of the corresponding frequency; the second is the characteristic curve of the water pump at fixed frequency Fitting, by collecting the characteristic parameters of the water pump at the same frequency, and fitting the characteristic curve at that frequency; finally, the fitting of the pipeline characteristic curve, by measuring the pipeline head at different flow rates, Fit the pipeline characteristic curve in this scenario. This invention integrates and automates data collection, data processing, data inspection, and real-time updating of water pump characteristic curves, and no longer requires R&D personnel to perform additional work; it can predict one or more future events through collected data at the same frequency and different frequencies. The flow-head and flow efficiency curve of the water pump within a time unit. This method reduces the work complexity while maintaining the accuracy of the characteristic curve. Moreover, the application scenarios can be generalized and are not only applicable to a specific industrial scenario.

本发明通过实时采集数据,并基于实时采集的数据对水泵和管路的特性曲线进行校准,模型自适应更新,并对采集的历史数据进行分析,根据改进的时间序列的ARMA模型,判断当前时刻采集的数据是否异常,进一步预测未来的流量-扬程和流量-效率曲线,对水泵的异常检测等问题可以作为一个判断条件;适用于不同工业环境下,可以同时对定频和变频采集的数据分别拟合指定频率下的水泵特性曲线。综上所述,本发明具有能够有效提高准确性以及扩大适用范围的特点。This invention collects data in real time and calibrates the characteristic curves of water pumps and pipelines based on the real-time collected data, updates the model adaptively, analyzes the collected historical data, and determines the current moment based on the improved time series ARMA model. Whether the collected data is abnormal, further predict the future flow-head and flow-efficiency curves, and can be used as a judgment condition for abnormal detection of water pumps. It is suitable for different industrial environments and can simultaneously analyze data collected at fixed frequency and variable frequency. Fit the water pump characteristic curve at the specified frequency. To sum up, the present invention has the characteristics of effectively improving accuracy and expanding the scope of application.

附图说明Description of the drawings

图1本发明实施例提供的一种水泵特性曲线自适应校准及预测方法的流程图。Figure 1 is a flow chart of an adaptive calibration and prediction method for a water pump characteristic curve provided by an embodiment of the present invention.

具体实施方式Detailed ways

以下结合附图对本发明具体实施方式作进一步详细说明。The specific embodiments of the present invention will be described in further detail below with reference to the accompanying drawings.

本发明提供的一种水泵特性曲线自适应校准及预测方法,The invention provides a method for adaptive calibration and prediction of water pump characteristic curves.

本发明实施的是一种水泵特性曲线自适应检验、校准以及预测的方法,包括数据采集、数据处理与合理性分析、水泵模型校准、水泵特性曲线预测四个方面,通过以下步骤实现,如图1所示为具体的流程图:What the present invention implements is a method for adaptive testing, calibration and prediction of water pump characteristic curves, which includes four aspects: data collection, data processing and rationality analysis, water pump model calibration, and water pump characteristic curve prediction. It is implemented through the following steps, as shown in the figure 1 shows the specific flow chart:

1.外部设备实时采集数据,每日采集一次数据,每条数据对应水泵/管路的流量、扬程等,将外部设备实时采集到的数据上传到指定平台,并对当前和历史数据储存。1. External equipment collects data in real time, and collects data once a day. Each piece of data corresponds to the flow rate, head, etc. of the water pump/pipeline. The data collected in real time by the external equipment is uploaded to the designated platform, and the current and historical data are stored.

对于每一个采集的数据点,通过多次采集取平均值的方式,作为该条数据的测量值。For each collected data point, the average value is taken through multiple collections as the measured value of the data.

2.基于分类改进的ARMA模型的建立。传统的ARMA模型基于每日采集的数据可以对未来多日的数据进行预测,然而,对于水泵特性曲线采集的数据点这类在短期内变化幅度小,需要在更长时间尺度上进行预测的数据,预测数据趋势时会存在趋势不清晰,并浪费了较多的计算资源。因此对传统的模型进行改进,充分挖掘数据的信息,提高数据预测的整体精度。改进的模型通过以下步骤进行实现:2. Establishment of ARMA model based on classification improvement. The traditional ARMA model can predict multi-day data in the future based on daily collected data. However, for data points collected from the water pump characteristic curve, which have small changes in the short term, predictions need to be made on a longer time scale. , when predicting data trends, the trend will be unclear and a lot of computing resources will be wasted. Therefore, the traditional model should be improved to fully mine the information of the data and improve the overall accuracy of data prediction. The improved model is implemented through the following steps:

对1中采集的数据分类,分为同一频率下的数据(定频)和不同频率下的数据(变频)两部分。The data collected in 1 is classified into two parts: data at the same frequency (fixed frequency) and data at different frequencies (variable frequency).

对两部分的数据分别按季度分类,并提取同一类中的各个特征(流量、扬程等)的信息,包括各个特征的最大值、最小值、均值、截尾均值(在同一类中去掉最大值和最小值后的均值),中位数,选取特征的中位数来描述该类的该特征,记为表示第i个类下的第j个特征点采集数据的中位数。Classify the two parts of data by quarter, and extract the information of each feature (flow rate, head, etc.) in the same category, including the maximum value, minimum value, mean value, and censored mean value of each feature (removing the maximum value in the same category and the mean after the minimum value), median, select the median of the feature to describe the feature of this class, recorded as Represents the median of the j-th feature point collection data under the i-th class.

对于每一类中各个特征的中位数可以构成向量对每一日采集的数据用向量可以表示为/>其中/>表示第i类下第t日的第q个特征,按照Legendre最佳逼近准则,以式(1-1)为目标,找到与数据中位数向量最接近的采集数据向量,并记录下对应的本类中的日期di,tThe median of each feature in each category can form a vector The data collected on each day can be expressed as a vector as/> Among them/> Represents the q-th feature on the t-th day under the i-th category. According to Legendre's best approximation criterion and taking formula (1-1) as the goal, find the collected data vector closest to the data median vector, and record the corresponding Date d i,t in this class.

在分类后的数据集中取前N个类下的m个特征,基于ARMA模型对当前类(第N+1类)的这m个特征点分别进行预测。具体公式如下:Take m features from the first N categories in the classified data set, and predict the m feature points of the current class (N+1th class) based on the ARMA model. The specific formula is as follows:

其中,p为自回归项AR的阶数,可由PACF图确定,αi(i=0,1,K,p)自回归系数;q为移动平均项MA的阶数,可由ACF图确定,βj(j=1,2,K,q)为移动平均系数,ξt为白噪声。Among them, p is the order of the autoregressive term AR, which can be determined from the PACF chart, α i (i=0,1,K,p) autoregressive coefficient; q is the order of the moving average term MA, which can be determined from the ACF chart, β j (j=1,2,K,q) is the moving average coefficient, and ξ t is white noise.

用测量得到的各个特征变量的值带入上述步骤中可以得到预测的下一个类的各个特征的中位数,对应的日期为 By bringing the measured values of each feature variable into the above steps, you can get the predicted median of each feature of the next class, and the corresponding date is

基于测量的第N个类中的最后kt日(分为k段)的特征值和预测的未来的中位数向量,拟合各个特征的变化曲线(用函数fi来表述),以15日为周期进行预测,实现对下一季度每15日进行一次水泵特性曲线的更新和迭代。具体的实现过程通过公式(1-4)描述,用最小二乘的方式求解公式(1-4),得到函数fi Based on the measured characteristic values of the last kt days (divided into k segments) in the Nth class and the predicted future median vector, the change curve of each characteristic (expressed by the function f Forecast for the cycle, and update and iterate the water pump characteristic curve every 15 days for the next quarter. The specific implementation process is described by formula (1-4), and the least square method is used to solve formula (1-4) to obtain the function f i

其中函数fi表示第i个特征拟合的预测函数,基于预测函数得到下一个类中任意时刻的预测值。The function fi represents the prediction function of the i-th feature fitting, and the predicted value at any time in the next class is obtained based on the prediction function.

3.基于步骤2中各个特征的预测值以及当前时刻各数据的测量值,带入水泵在额定频率下的运行规律公式包括流量-扬程曲线方程、流量-功率曲线方程以及效率曲线方程,如下所示:3. Based on the predicted values of each feature in step 2 and the measured values of each data at the current moment, the operating formulas of the water pump at the rated frequency include the flow-head curve equation, the flow-power curve equation and the efficiency curve equation, as follows Show:

H=aQ2+bQ+c (1-5)H=aQ 2 +bQ+c (1-5)

P=jQ2+kQ+l (1-6)P=jQ 2 +kQ+l (1-6)

其中H是水泵扬程,Q是水泵流量,η是水泵效率,P是水泵功率,a,b,c,j,k,l是拟合二次曲线的参数。针对不同的工业场景以及采集的数据,分别处理同一频率以及不同频率下采集的数据。Among them, H is the water pump head, Q is the water pump flow rate, eta is the water pump efficiency, P is the water pump power, a, b, c, j, k, l are the parameters of fitting the quadratic curve. For different industrial scenarios and collected data, data collected at the same frequency and at different frequencies are processed respectively.

同一频率下采集的数据通过公式(1-5)(1-6)可以得到特定频率下的流量-扬程曲线方程和流量-功率曲线方程。The data collected at the same frequency can be used to obtain the flow-head curve equation and flow-power curve equation at a specific frequency through formulas (1-5) (1-6).

不同频率下采集的数据,通过对采集数据以最小二乘的处理方式结合机理模型进行拟合,如下所示:The data collected at different frequencies are fitted by combining the least squares processing method with the mechanism model, as shown below:

H0=aQ0 2+bQ0i+ci2 (1-9)H 0 =aQ 0 2 +bQ 0 i+ci 2 (1-9)

P0=jQ0 2i+kQ0i2+li3 (1-10)P 0 =jQ 0 2 i+kQ 0 i 2 +li 3 (1-10)

其中H0,Q0,P0为在水泵运行频率为f0时的扬程、流量以及功率,f为额定频率,i为运行频率与额定频率的比例系数。Among them, H 0 , Q 0 , P 0 are the head, flow rate and power when the pump operating frequency is f 0 , f is the rated frequency, and i is the proportional coefficient between the operating frequency and the rated frequency.

假设采集了n个不同频率时的数据,由(1-5),(1-9)和(1-6),(1-10)通过最小二乘法计算:Assume that n data at different frequencies are collected and calculated by the least squares method from (1-5), (1-9) and (1-6), (1-10):

Ax=T (1-11)Ax=T (1-11)

By=W (1-12)By=W (1-12)

其中:in:

通过该步骤不仅可以对当前时刻水泵的特性曲线进行修正,还能够得到预测的下一个时间单位内水泵运行时的流量-扬程和流量-效率曲线。This step can not only correct the characteristic curve of the water pump at the current moment, but also obtain the predicted flow-head and flow-efficiency curves when the water pump is running in the next time unit.

4.基于采集的管路参数通过公式(1-14)以最小二乘的方式进行拟合4. Based on the collected pipeline parameters, the least squares method is used to fit the formula (1-14).

其中为管路的扬程,Q为管路的流量,D,S为需要拟合的参数。in is the head of the pipeline, Q is the flow rate of the pipeline, D and S are the parameters that need to be fitted.

绘制管路特性曲线和水泵特性曲线,它们的交点就是预测水泵工作的实际工况点。Draw the pipeline characteristic curve and the water pump characteristic curve. Their intersection is the actual working condition point where the water pump is predicted to work.

本发明实施例可以实现以下功能:平台上的每个组件都应该可以独立地和平台外的系统进行交互,实时的采集和传输数据,并上传到数据处理平台;Embodiments of the present invention can achieve the following functions: each component on the platform should be able to independently interact with systems outside the platform, collect and transmit data in real time, and upload it to the data processing platform;

基于实时采集的数据,完成水泵以及管路特性曲线自动化拟合更正的脚本,计算实时的管路阻抗系数,并对水泵原先的机理模型进行实时的更新和修正;Based on the data collected in real time, complete the script for automatic fitting and correction of the water pump and pipeline characteristic curves, calculate the real-time pipeline impedance coefficient, and update and correct the original mechanism model of the water pump in real time;

基于历史的数据,对水泵特性曲线的变化情况进行预测;Based on historical data, predict the changes in the water pump characteristic curve;

对于不同的工业环境都具有适用性,对采集的数据分块处理,得到更为精确的水泵特性曲线。It is applicable to different industrial environments. The collected data is processed in blocks to obtain a more accurate water pump characteristic curve.

为了满足上述功能,需外部系统实时的采集和输入数据,并实时更新该水泵的特性曲线,而无需平台的开发人员额外进行人工处理,从而避免了平台的二次开发,提升效率。In order to meet the above functions, an external system is required to collect and input data in real time, and update the characteristic curve of the water pump in real time, without the need for additional manual processing by the platform developers, thereby avoiding secondary development of the platform and improving efficiency.

上述实施例用来解释说明本发明,而不是对本发明进行限制,在本发明的精神和权利要求的保护范围内,对本发明作出的任何修改和改变,都落入本发明的保护范围。The above embodiments are used to illustrate the present invention, rather than to limit the present invention. Within the spirit of the present invention and the protection scope of the claims, any modifications and changes made to the present invention fall within the protection scope of the present invention.

Claims (5)

1.一种水泵特性曲线自适应校准及预测方法,其特征在于,包括以下步骤:1. An adaptive calibration and prediction method for water pump characteristic curves, which is characterized by including the following steps: S1、经外部设备实时采集水泵和管路两部分各个特征的数据,将采集到的数据上传到指定平台,并对当前和历史数据储存;S1. Collect data on various characteristics of the water pump and pipeline in real time through external equipment, upload the collected data to the designated platform, and store current and historical data; S2、建立基于分类改进的ARMA模型,并对水泵各个特征进行预测;S2. Establish an ARMA model based on classification improvement and predict various characteristics of the water pump; S3、将步骤S2中水泵各个特征的预测值以及当前时刻水泵各个特征的采集值,代入水泵在额定频率下的运行规律公式包括流量-扬程曲线方程、流量-功率曲线方程以及效率曲线方程,对当前时刻水泵的特性曲线进行修正,并得到预测的下一个时间单位内水泵运行时的流量-扬程和流量-效率曲线;S3. Substitute the predicted values of each characteristic of the water pump in step S2 and the collected values of each characteristic of the water pump at the current moment into the operating law formula of the water pump at the rated frequency, including the flow-head curve equation, the flow-power curve equation and the efficiency curve equation. The characteristic curve of the water pump at the current moment is corrected, and the predicted flow-head and flow-efficiency curves when the water pump is running in the next time unit are obtained; S4、基于采集的管路参数以最小二乘的方式拟合得到管路特性曲线,拟合公式为其中/>为管路的扬程,Q为管路的流量,D,S为需要拟合的参数;绘制管路特性曲线和水泵特性曲线,两条曲线的交点为预测水泵工作的实际工况点。S4. Based on the collected pipeline parameters, the pipeline characteristic curve is obtained by fitting the least squares method. The fitting formula is: Among them/> is the head of the pipeline, Q is the flow rate of the pipeline, D and S are the parameters that need to be fitted; draw the pipeline characteristic curve and the water pump characteristic curve, and the intersection of the two curves is the actual working condition point for predicting the work of the water pump. 2.根据权利要求1所述的一种水泵特性曲线自适应校准及预测方法,其特征在于,步骤S1中,每日采集一次水泵和管路各个特征的数据,每条数据对应水泵和管路的流量、扬程;2. A water pump characteristic curve adaptive calibration and prediction method according to claim 1, characterized in that, in step S1, data on various characteristics of the water pump and pipeline are collected once a day, and each piece of data corresponds to the water pump and pipeline. The flow rate and head; 对于每一个采集的数据,通过多次采集取平均值的方式,作为该条数据的测量值。For each collected data, the average value is taken through multiple collections as the measured value of this piece of data. 3.根据权利要求1所述的一种水泵特性曲线自适应校准及预测方法,其特征在于,步骤S2中,基于分类改进的ARMA模型通过以下步骤进行实现:3. A water pump characteristic curve adaptive calibration and prediction method according to claim 1, characterized in that, in step S2, the ARMA model based on classification improvement is implemented through the following steps: 1)对步骤S1中采集的水泵各个特征的数据进行分类,分为同一频率下的数据即定频数据和不同频率下的数据即变频数据;1) Classify the data of each characteristic of the water pump collected in step S1 into data at the same frequency, that is, fixed frequency data, and data at different frequencies, that is, variable frequency data; 2)对定频数据和变频数据分别按季度分类,并提取同一类中的各个特征数据信息,包括各个特征的最大值、最小值、均值、截尾均值和中位数,选取特征的其中一个数据信息来描述该类的该特征,记为表示第i个类下的第j个特征点采集数据的中位数;2) Classify fixed-frequency data and variable-frequency data by quarter, and extract each feature data information in the same category, including the maximum value, minimum value, mean, censored mean and median of each feature, and select one of the features Data information is used to describe the characteristics of this class, recorded as Represents the median of the j-th feature point collection data under the i-th class; 3)将每一类中各个特征的数据信息构成向量,根据Legendre最佳逼近准则,找到与特征数据信息向量最接近的采集数据向量,并记录下对应的本类中的日期di,t3) Construct the data information of each feature in each category into a vector, and according to the Legendre best approximation criterion, find the collection data vector closest to the feature data information vector, and record the corresponding date d i,t in this category; 4)在分类后的数据集中取前N个类下的m个特征,基于ARMA模型对当前类即第N+1类的这m个特征点分别进行预测,具体公式如下:4) Take m features from the first N categories in the classified data set, and predict the m feature points of the current class, that is, the N+1th class, based on the ARMA model. The specific formula is as follows: 其中,p为自回归项AR的阶数,αi(i=0,1,K,p)自回归系数;q为移动平均项MA的阶数,βj(j=1,2,K,q)为移动平均系数,ξt为白噪声,限定了噪声的边界;E表示数学期望,Var表示方差,/>表示任意采集的同一类型的水泵数据;Among them, p is the order of the autoregressive term AR, α i (i=0,1,K,p) autoregressive coefficient; q is the order of the moving average term MA, β j (j=1,2,K, q) is the moving average coefficient, ξ t is white noise, Defines the boundary of noise; E represents mathematical expectation, Var represents variance,/> Represents any collected water pump data of the same type; 5)用采集得到的各个特征变量的值代入上述步骤中得到预测的下一个类的各个特征的数据信息,对应的日期为 5) Substitute the collected values of each feature variable into the above steps to obtain the predicted data information of each feature of the next class. The corresponding date is 6)基于测量的第N个类中的最后kt日的特征值和预测的未来的特征数据信息向量,拟合各个特征的变化曲线,用函数fi来表述,选定预测周期,实现对下一季度每一次预测周期进行一次水泵特性曲线的更新和迭代;具体公式如下:6) Based on the measured characteristic values of the last kt day in the Nth class and the predicted future characteristic data information vector, fit the change curve of each characteristic, express it with the function f i , select the prediction period, and realize the next prediction The water pump characteristic curve is updated and iterated for each forecast period in the first quarter; the specific formula is as follows: 用最小二乘的方式求解公式,得到函数fi,其中,函数fi表示第i个特征拟合的预测函数,基于预测函数得到下一个类中任意时刻的预测值。Use the least squares method to solve the formula to obtain the function f i , where the function f i represents the prediction function of the i-th feature fitting, and the predicted value at any time in the next class is obtained based on the prediction function. 4.根据权利要求1所述的一种水泵特性曲线自适应校准及预测方法,其特征在于,步骤S3中,流量-扬程曲线方程、流量-功率曲线方程以及效率曲线方程,分别如下所示:4. A water pump characteristic curve adaptive calibration and prediction method according to claim 1, characterized in that, in step S3, the flow-head curve equation, the flow-power curve equation and the efficiency curve equation are as follows: H=aQ2+bQ+cH=aQ 2 +bQ+c P=jQ2+kQ+lP=jQ 2 +kQ+l 其中,H是水泵扬程,Q是水泵流量,η是水泵效率,P是水泵功率,a,b,c,j,k,l是拟合二次曲线的参数。Among them, H is the water pump head, Q is the water pump flow rate, eta is the water pump efficiency, P is the water pump power, a, b, c, j, k, l are the parameters of fitting the quadratic curve. 5.根据权利要求4所述的一种水泵特性曲线自适应校准及预测方法,其特征在于,同一频率下采集的数据,通过流量-扬程曲线方程、流量-功率曲线方程能够得到特定频率下的流量-扬程曲线方程和流量-功率曲线方程;5. A water pump characteristic curve adaptive calibration and prediction method according to claim 4, characterized in that the data collected at the same frequency can be obtained at a specific frequency through the flow-head curve equation and the flow-power curve equation. Flow-head curve equation and flow-power curve equation; 不同频率下采集的数据,通过对采集数据以最小二乘的处理方式结合机理模型进行拟合,如下所示:The data collected at different frequencies are fitted by combining the least squares processing method with the mechanism model, as shown below: H0=aQ0 2+bQ0i+ci2 H 0 =aQ 0 2 +bQ 0 i+ci 2 P0=jQ0 2i+kQ0i2+li3 P 0 =jQ 0 2 i+kQ 0 i 2 +li 3 其中,H0,Q0,P0为在水泵运行频率为f0时的扬程、流量以及功率,f为额定频率,i为运行频率与额定频率的比例系数。Among them, H 0 , Q 0 , P 0 are the head, flow rate and power when the pump operating frequency is f 0 , f is the rated frequency, and i is the proportional coefficient between the operating frequency and the rated frequency.
CN202310697387.XA 2023-06-13 2023-06-13 Water pump characteristic curve self-adaptive calibration and prediction method Pending CN116882079A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310697387.XA CN116882079A (en) 2023-06-13 2023-06-13 Water pump characteristic curve self-adaptive calibration and prediction method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310697387.XA CN116882079A (en) 2023-06-13 2023-06-13 Water pump characteristic curve self-adaptive calibration and prediction method

Publications (1)

Publication Number Publication Date
CN116882079A true CN116882079A (en) 2023-10-13

Family

ID=88255758

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310697387.XA Pending CN116882079A (en) 2023-06-13 2023-06-13 Water pump characteristic curve self-adaptive calibration and prediction method

Country Status (1)

Country Link
CN (1) CN116882079A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117540163A (en) * 2023-11-09 2024-02-09 南栖仙策(南京)高新技术有限公司 Pump performance curve generation method, model construction and training method and device
CN118095053A (en) * 2024-01-15 2024-05-28 上海碳索能源服务股份有限公司 Indirect prediction method, system, terminal and medium for flow of cooling water pump
CN119042683A (en) * 2024-10-30 2024-11-29 浙江中广电器集团股份有限公司 Heat pump unit, pipe network calibration method thereof and actual operation environment determination method
CN117540163B (en) * 2023-11-09 2025-02-25 南栖仙策(南京)高新技术有限公司 Pump performance curve generation method, model building and training method and device

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117540163A (en) * 2023-11-09 2024-02-09 南栖仙策(南京)高新技术有限公司 Pump performance curve generation method, model construction and training method and device
CN117540163B (en) * 2023-11-09 2025-02-25 南栖仙策(南京)高新技术有限公司 Pump performance curve generation method, model building and training method and device
CN118095053A (en) * 2024-01-15 2024-05-28 上海碳索能源服务股份有限公司 Indirect prediction method, system, terminal and medium for flow of cooling water pump
CN119042683A (en) * 2024-10-30 2024-11-29 浙江中广电器集团股份有限公司 Heat pump unit, pipe network calibration method thereof and actual operation environment determination method

Similar Documents

Publication Publication Date Title
CN111913803B (en) Service load fine granularity prediction method based on AKX hybrid model
CN116882079A (en) Water pump characteristic curve self-adaptive calibration and prediction method
CN111275288A (en) XGboost-based multi-dimensional data anomaly detection method and device
CN107765347A (en) A kind of Gaussian process returns and the short-term wind speed forecasting method of particle filter
CN111525587B (en) A method and system for grid reactive voltage control based on reactive load situation
CN117668743B (en) Time sequence data prediction method of association time-space relation
CN118469332B (en) A multi-type load collaborative control method and system for industrial parks
CN113449919A (en) Power consumption prediction method and system based on feature and trend perception
CN116646929A (en) A short-term wind power forecasting method based on PSO-CNN-BILSTM
CN113408808B (en) Training method, data generation device, electronic equipment and storage medium
CN117494888A (en) Distributed photovoltaic power prediction method and system based on optimized cascade
CN117216572A (en) Wind-solar power generation load prediction method for multi-scene application
CN118378553B (en) A digital twin model construction method and system
CN117175595B (en) Power grid regulation and control method and system based on multi-level data
CN114372615A (en) Short-term photovoltaic power prediction method and system considering spatial correlation
CN111027612B (en) Energy metering data feature reduction method and device based on weighted entropy FCM
CN118971349A (en) Intelligent monitoring method and system for power distribution network based on data center and edge computing
CN106228277B (en) Reservoir dispatching forecast information effective precision identification method based on data mining
CN113742989B (en) Combustion optimization control method and device, storage medium and electronic equipment
CN116578858A (en) Air compressor fault prediction and health degree evaluation method and system based on graphic neural network
CN111626465B (en) New energy power short-term interval prediction method and system
CN115114983B (en) Method for acquiring and analyzing electric quantity data based on big data equipment and computer system
CN118567240B (en) Adaptive control method and system for aluminum smelting process parameters based on big data modeling
CN110322063A (en) A kind of power consumption simulated prediction method and storage medium
CN109670243A (en) A kind of life-span prediction method based on lebesgue space model

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