CN113468751B - Recursion Lasso-based flowmeter anomaly online monitoring method and system and storage medium - Google Patents
Recursion Lasso-based flowmeter anomaly online monitoring method and system and storage medium Download PDFInfo
- Publication number
- CN113468751B CN113468751B CN202110765956.0A CN202110765956A CN113468751B CN 113468751 B CN113468751 B CN 113468751B CN 202110765956 A CN202110765956 A CN 202110765956A CN 113468751 B CN113468751 B CN 113468751B
- Authority
- CN
- China
- Prior art keywords
- flow
- flowmeter
- autoregressive model
- lasso
- model
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 74
- 238000012544 monitoring process Methods 0.000 title claims abstract description 56
- 238000003860 storage Methods 0.000 title claims abstract description 13
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 39
- 230000002159 abnormal effect Effects 0.000 claims abstract description 35
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 13
- 238000004519 manufacturing process Methods 0.000 claims abstract description 12
- 241000208125 Nicotiana Species 0.000 claims abstract description 9
- 235000002637 Nicotiana tabacum Nutrition 0.000 claims abstract description 9
- 230000008569 process Effects 0.000 claims description 18
- 238000004590 computer program Methods 0.000 claims description 9
- 238000002790 cross-validation Methods 0.000 claims description 3
- 230000002547 anomalous effect Effects 0.000 claims 1
- 230000005856 abnormality Effects 0.000 description 37
- 238000010586 diagram Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 235000019504 cigarettes Nutrition 0.000 description 3
- 238000013480 data collection Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- 241000723345 Chrysophyllum Species 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000013598 vector Substances 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/20—Administration of product repair or maintenance
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/04—Manufacturing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Economics (AREA)
- Marketing (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Entrepreneurship & Innovation (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Measuring Volume Flow (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
Description
技术领域technical field
本发明涉及烟草加工技术领域,尤其涉及一种基于递推Lasso的流量计异常在线监测方法、系统和存储介质。The invention relates to the technical field of tobacco processing, in particular to a recursive Lasso-based flowmeter abnormality online monitoring method, system and storage medium.
背景技术Background technique
计量器具产品在使用的过程中会出现异常,在异常出现之后如果没有及时进行排查和维修,将会影响计量器具产品的准确性,从而对后续的工业生产过程带来巨大的隐患,也会增加企业的管理支出成本。There will be abnormalities in the use of measuring instruments. If they are not checked and repaired in time after the abnormality occurs, the accuracy of the measuring instruments will be affected, which will bring huge hidden dangers to the subsequent industrial production process and increase The cost of administrative expenses for the business.
为了确保计量器具测量数据的准确性,就需要在计量器具工作过程中实时监测计量器具测得的数据。一旦出现异常,如数据跳变等情况,需要判断是计量器具出现故障,抑或是生产条件的调整;若计量器具出现故障,则需要对计量器具进行异常排查和维修。In order to ensure the accuracy of the measurement data of the measuring instruments, it is necessary to monitor the data measured by the measuring instruments in real time during the working process of the measuring instruments. Once an abnormality occurs, such as data jumps, etc., it is necessary to judge whether the measuring instrument is faulty or the adjustment of the production conditions; if the measuring instrument is faulty, it is necessary to check and repair the abnormality of the measuring instrument.
在线流量计是卷烟制造过程溯源链中最重要的基础数据获取源之一,其检测性能在生产过程中需要始终维持在允许范围内。目前,卷烟企业生产线对流量计的异常监测方式缺乏专业性,一方面仅是简单参考常规点检规程,而忽视了生产过程中的在线监测特点;另一方面则是仅凭借技术人员的维护经验进行主观预判,具有一定随意性,且存在人为干预导致的失误风险,因此无法确保流量计是否真正出现故障。可见,由于缺乏科学依据和数据支撑,现有的流量计异常监测方法未能较好的做到故障预警的作用,不能及时发现流量计故障,因此存在较大的生产隐患。The online flowmeter is one of the most important sources of basic data acquisition in the traceability chain of the cigarette manufacturing process, and its detection performance needs to be kept within the allowable range during the production process. At present, the production line of cigarette enterprises lacks professionalism in the abnormal monitoring methods of flowmeters. On the one hand, they simply refer to the routine point inspection procedures, while ignoring the characteristics of online monitoring in the production process; on the other hand, they only rely on the maintenance experience of technicians Subjective prediction is somewhat arbitrary, and there is a risk of error caused by human intervention, so it is impossible to ensure whether the flowmeter is actually faulty. It can be seen that due to the lack of scientific basis and data support, the existing flowmeter anomaly monitoring methods fail to achieve a good fault warning function, and cannot detect flowmeter failures in time, so there are large hidden dangers in production.
因此,亟需一种基于递推Lasso的流量计异常在线监测方法、系统和存储介质。Therefore, there is an urgent need for a recursive Lasso-based flowmeter abnormality online monitoring method, system and storage medium.
发明内容Contents of the invention
本发明的目的是提供一种基于递推Lasso的流量计异常在线监测方法、系统和存储介质,以解决上述现有技术中的问题,能够及时监测到流量计的异常情况。The object of the present invention is to provide a recursive Lasso-based flowmeter abnormality online monitoring method, system and storage medium to solve the above-mentioned problems in the prior art, and to be able to monitor the abnormality of the flowmeter in time.
本发明提供了一种基于递推Lasso的流量计异常在线监测方法,其中,包括:The present invention provides a flowmeter abnormal online monitoring method based on recursive Lasso, which includes:
根据流量计的历史流量数据,建立自回归模型,并通过离线Lasso算法确定所述自回归模型的模型系数;According to the historical flow data of flow meter, set up autoregressive model, and determine the model coefficient of described autoregressive model by off-line Lasso algorithm;
利用流量计的实时流量,通过递推Lasso算法对所述自回归模型进行更新;Using the real-time flow rate of the flowmeter, the autoregressive model is updated through the recursive Lasso algorithm;
根据基于更新后的所述自回归模型得到的预测流量和流量计在下一时刻的实际流量确定流量计是否异常。Whether the flowmeter is abnormal is determined according to the predicted flow obtained based on the updated autoregressive model and the actual flow of the flowmeter at the next moment.
如上所述的基于递推Lasso的流量计异常在线监测方法,其中,优选的是,所述根据流量计的历史流量数据,建立自回归模型,并通过离线Lasso算法确定所述自回归模型的模型系数,具体包括:The above-mentioned recursive Lasso-based flowmeter abnormality online monitoring method, wherein, preferably, the autoregressive model is established according to the historical flow data of the flowmeter, and the model of the autoregressive model is determined by the offline Lasso algorithm coefficients, including:
采集流量计流量的数据样本,得到历史流量数据集合;Collect data samples of flowmeter flow to obtain historical flow data collection;
根据所述历史流量数据集合,构造自回归模型;Construct an autoregressive model according to the historical flow data set;
基于离线Lasso算法求取所述自回归模型的回归系数;Obtain the regression coefficient of described autoregressive model based on off-line Lasso algorithm;
利用交叉验证的方法确定所述自回归模型的的正则化参数。The regularization parameters of the autoregressive model are determined by means of cross-validation.
如上所述的基于递推Lasso的流量计异常在线监测方法,其中,优选的是,通过以下公式确定所述离线Lasso算法的最优条件,并根据所述离线Lasso算法的最优条件确定所述自回归模型的回归系数和正则化参数:The above-mentioned recursive Lasso-based flowmeter abnormality online monitoring method, wherein, preferably, the optimal condition of the offline Lasso algorithm is determined by the following formula, and the optimal condition of the offline Lasso algorithm is determined according to the Regression coefficients and regularization parameters for an autoregressive model:
其中,p为自回归模型的滞后系数,yn表示当前时刻待预测的流量值,yn-i表示历史时刻的流量值,α表示自回归模型的回归系数,μn-1表示自回归模型的正则化参数。Among them, p is the lag coefficient of the autoregressive model, y n represents the flow value to be predicted at the current moment, y ni represents the flow value at the historical moment, α represents the regression coefficient of the autoregressive model, μ n-1 represents the regularity of the autoregressive model parameterization.
如上所述的基于递推Lasso的流量计异常在线监测方法,其中,优选的是,所述利用流量计的实时流量,通过递推Lasso算法对所述自回归模型进行更新,具体包括:The above-mentioned recursive Lasso-based flowmeter abnormal online monitoring method, wherein, preferably, the real-time flow of the flowmeter is used to update the autoregressive model through the recursive Lasso algorithm, specifically including:
将通过离线Lasso算法所确定的所述自回归模型的模型系数作为所述自回归模型在更新过程中的初始模型系数;Using the model coefficients of the autoregressive model determined by the offline Lasso algorithm as the initial model coefficients of the autoregressive model in the update process;
每读入一个新的流量计实时流量,则根据该新的流量计实时流量对所述自回归模型的模型参数进行更新,并且通过以下公式对所述自回归模型进行更新:Whenever a new real-time flow of the flowmeter is read, the model parameters of the autoregressive model are updated according to the new real-time flow of the flowmeter, and the autoregressive model is updated by the following formula:
其中,Yn-1表示流量计测得流量的历史流量数据集合,Zn-1表示自回归模型的自回归项,yn表示当前时刻待预测的流量值,Zn表示n时刻对应的自回归项,t为递推Lasso算法的参数,其取值为0-1,在读入新数据后,t的值都会从0变到1,当t=1时表示新数据完全读入,自回归模型内的参数α和μ也更新完毕,准备读入下一个时刻的新数据,当下一时刻的新数据读入之后,t的值又从0开始变到1。Among them, Y n-1 represents the historical flow data collection of flow measured by the flowmeter, Z n-1 represents the autoregressive item of the autoregressive model, y n represents the flow value to be predicted at the current moment, and Z n represents the autoregressive value corresponding to time n. Regression item, t is the parameter of the recursive Lasso algorithm, and its value is 0-1. After the new data is read in, the value of t will change from 0 to 1. When t=1, it means that the new data is completely read in. The parameters α and μ in the regression model are also updated, and the new data at the next moment is ready to be read in. After the new data at the next moment is read in, the value of t changes from 0 to 1 again.
如上所述的基于递推Lasso的流量计异常在线监测方法,其中,优选的是,所述根据基于更新后的所述自回归模型得到的预测流量和流量计在下一时刻的实际流量确定流量计是否异常,具体包括:In the online monitoring method for flowmeter abnormalities based on recursive Lasso as described above, preferably, the flowmeter is determined according to the predicted flow rate obtained based on the updated autoregressive model and the actual flow rate of the flowmeter at the next moment. Is it abnormal, including:
根据更新后的所述自回归模型,实时预测流量计下一时刻的流量,得到预测流量;According to the updated autoregressive model, the flow rate of the flow meter at the next moment is predicted in real time to obtain the predicted flow rate;
根据得到的不同时刻对应的预测流量,得到预测流量的估计曲线图;Obtain an estimated curve diagram of the predicted flow according to the obtained predicted flow corresponding to different moments;
根据所述预测流量和流量计在下一时刻的实际流量的偏差,得到残差图;According to the deviation between the predicted flow and the actual flow of the flowmeter at the next moment, a residual map is obtained;
根据所述残差图确定流量计是否异常。Whether the flowmeter is abnormal is determined according to the residual error graph.
如上所述的基于递推Lasso的流量计异常在线监测方法,其中,优选的是,所述根据所述残差图确定流量计是否异常,具体包括:According to the recursive Lasso-based flowmeter abnormality online monitoring method as described above, preferably, the determining whether the flowmeter is abnormal according to the residual map specifically includes:
在所述残差图的纵坐标超过预设阈值时,则确定流量计的流量数据异常;When the ordinate of the residual graph exceeds a preset threshold, it is determined that the flow data of the flowmeter is abnormal;
根据生产记录确定导致流量计的流量数据异常的原因是流量计出现异常还是梗丝产品批次发生更换。According to the production records, it is determined whether the abnormality of the flow data of the flowmeter is caused by the abnormality of the flowmeter or the replacement of the shredded stem product batch.
如上所述的基于递推Lasso的流量计异常在线监测方法,其中,优选的是,所述基于递推Lasso的流量计异常在线监测方法用于监测烟草制丝过程增湿水流量是否异常,The above-mentioned recursive Lasso-based flowmeter abnormality online monitoring method, wherein, preferably, the recursive Lasso-based flowmeter abnormality online monitoring method is used to monitor whether the humidification water flow in the shredded tobacco process is abnormal,
所述基于递推Lasso的流量计异常在线监测方法,还包括:The recursive Lasso-based flowmeter abnormal online monitoring method also includes:
基于预设的梗丝产品批次的流量阈值和所述预测流量,确定当前时刻的梗丝产品规格。Based on the preset flow threshold of the cut stem product batch and the predicted flow, the specification of the cut stem product at the current moment is determined.
本发明还提供一种采用上述方法的基于递推Lasso的流量计异常在线监测系统,包括:The present invention also provides a flowmeter abnormal online monitoring system based on recursive Lasso using the above method, including:
自回归模型建立模块,用于根据流量计的历史流量数据,建立自回归模型,并通过离线Lasso算法确定所述自回归模型的模型系数;The autoregressive model building module is used to establish an autoregressive model according to the historical flow data of the flowmeter, and determine the model coefficients of the autoregressive model by an off-line Lasso algorithm;
自回归模型更新模块,用于利用流量计的实时流量,通过递推Lasso算法对所述自回归模型进行更新;The autoregressive model update module is used to update the autoregressive model through the recursive Lasso algorithm by using the real-time flow of the flowmeter;
流量监测模块,用于根据基于更新后的所述自回归模型得到的预测流量和流量计在下一时刻的实际流量确定流量计是否异常。The flow monitoring module is used to determine whether the flow meter is abnormal according to the predicted flow obtained based on the updated autoregressive model and the actual flow of the flow meter at the next moment.
本发明还提供一种计算机可读存储介质,所述计算机可读存储介质中存储有计算机程序,当其在计算机上运行时,使得计算机执行上述的基于递推Lasso的流量计异常在线监测方法。The present invention also provides a computer-readable storage medium, wherein a computer program is stored in the computer-readable storage medium, and when it is run on a computer, the computer executes the above-mentioned recursive Lasso-based flowmeter abnormality online monitoring method.
本发明还提供一种计算机程序产品,其特征在于,所述计算机程序产品在终端设备上运行时,使所述终端设备执行上述的基于递推Lasso的流量计异常在线监测方法。The present invention also provides a computer program product, which is characterized in that, when the computer program product is run on a terminal device, the terminal device is made to execute the above-mentioned recursive Lasso-based flowmeter abnormality online monitoring method.
本发明提供一种基于递推Lasso的流量计异常在线监测方法,通过递推Lasso算法对自回归模型进行更新,根据基于更新后的自回归模型得到的预测流量和实际流量来监测流量计的异常情况,具有准确度高、操作方便和实时跟踪等特点,为烟草制丝过程增湿水流量的在线监测提供了科学、客观、可靠的技术支持,进而能够保证流量计的计量性能稳定性,并且本发明也适用于其他计量器具的检测场景中,应用前景广泛。The present invention provides a flowmeter abnormality online monitoring method based on recursive Lasso, which updates the autoregressive model through the recursive Lasso algorithm, and monitors the abnormality of the flowmeter according to the predicted flow and actual flow obtained based on the updated autoregressive model It has the characteristics of high accuracy, convenient operation and real-time tracking. It provides scientific, objective and reliable technical support for the on-line monitoring of humidification water flow in the tobacco shred process, thereby ensuring the stability of the measurement performance of the flowmeter, and The invention is also applicable to detection scenarios of other measuring instruments, and has wide application prospects.
附图说明Description of drawings
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步描述,其中:In order to make the purpose, technical solutions and advantages of the present invention more clear, the present invention will be further described below in conjunction with accompanying drawing, wherein:
图1为本发明提供的基于递推Lasso的流量计异常在线监测方法的实施例的流程图;Fig. 1 is the flowchart of the embodiment of the flow meter abnormal online monitoring method based on recursive Lasso provided by the present invention;
图2为本发明实施例提供的预测流量的估计曲线图;Fig. 2 is the estimation graph of the predicted flow provided by the embodiment of the present invention;
图3为本发明实施例提供的残差图;FIG. 3 is a residual diagram provided by an embodiment of the present invention;
图4为本发明提供的基于递推Lasso的流量计异常在线监测系统的实施例的结构框图。Fig. 4 is a structural block diagram of an embodiment of an online monitoring system for flowmeter abnormality based on recursive Lasso provided by the present invention.
具体实施方式detailed description
现在将参照附图来详细描述本公开的各种示例性实施例。对示例性实施例的描述仅仅是说明性的,决不作为对本公开及其应用或使用的任何限制。本公开可以以许多不同的形式实现,不限于这里所述的实施例。提供这些实施例是为了使本公开透彻且完整,并且向本领域技术人员充分表达本公开的范围。应注意到:除非另外具体说明,否则在这些实施例中阐述的部件和步骤的相对布置、材料的组分、数字表达式和数值应被解释为仅仅是示例性的,而不是作为限制。Various exemplary embodiments of the present disclosure will now be described in detail with reference to the accompanying drawings. The description of the exemplary embodiments is illustrative only, and in no way restricts the disclosure, its application or uses. The present disclosure can be implemented in many different forms and is not limited to the embodiments described here. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art. It should be noted that relative arrangements of parts and steps, compositions of materials, numerical expressions and numerical values set forth in these embodiments should be interpreted as illustrative only and not as limiting, unless specifically stated otherwise.
本公开中使用的“第一”、“第二”:以及类似的词语并不表示任何顺序、数量或者重要性,而只是用来区分不同的部分。“包括”或者“包含”等类似的词语意指在该词前的要素涵盖在该词后列举的要素,并不排除也涵盖其他要素的可能。“上”、“下”等仅用于表示相对位置关系,当被描述对象的绝对位置改变后,则该相对位置关系也可能相应地改变。"First", "second" and similar words used in the present disclosure do not indicate any order, quantity or importance, but are only used to distinguish different parts. Words like "comprising" or "comprising" mean that the elements preceding the word cover the elements listed after the word, and do not exclude the possibility of also covering other elements. "Up", "Down" and so on are only used to indicate the relative positional relationship. When the absolute position of the described object changes, the relative positional relationship may also change accordingly.
在本公开中,当描述到特定部件位于第一部件和第二部件之间时,在该特定部件与第一部件或第二部件之间可以存在居间部件,也可以不存在居间部件。当描述到特定部件连接其它部件时,该特定部件可以与所述其它部件直接连接而不具有居间部件,也可以不与所述其它部件直接连接而具有居间部件。In the present disclosure, when it is described that a specific component is located between a first component and a second component, there may or may not be an intervening component between the specific component and the first component or the second component. When it is described that a specific component is connected to other components, the specific component may be directly connected to the other component without an intermediate component, or may not be directly connected to the other component but has an intermediate component.
本公开使用的所有术语(包括技术术语或者科学术语)与本公开所属领域的普通技术人员理解的含义相同,除非另外特别定义。还应当理解,在诸如通用字典中定义的术语应当被解释为具有与它们在相关技术的上下文中的含义相一致的含义,而不应用理想化或极度形式化的意义来解释,除非这里明确地这样定义。All terms (including technical terms or scientific terms) used in the present disclosure have the same meaning as understood by one of ordinary skill in the art to which the present disclosure belongs, unless otherwise specifically defined. It should also be understood that terms defined in, for example, general-purpose dictionaries should be interpreted as having meanings consistent with their meanings in the context of the relevant technology, and should not be interpreted in idealized or extremely formalized meanings, unless explicitly stated herein Defined like this.
对于相关领域普通技术人员已知的技术、方法和设备可能不作详细讨论,但在适当情况下,技术、方法和设备应当被视为说明书的一部分。Techniques, methods and devices known to those of ordinary skill in the relevant art may not be discussed in detail, but where appropriate, techniques, methods and devices should be considered part of the description.
本发明针对流量计这一在卷烟生产过程中使用较为频繁的传感器件,基于实际应用场景和客观科学思想,提出一种基于递推Lasso的流量计异常在线监测方法对流量计进行在线监测,及时发现故障,减小损失。Aiming at the flowmeter, which is a sensor device that is frequently used in the cigarette production process, based on the actual application scene and objective scientific ideas, the present invention proposes an online monitoring method for flowmeter abnormality based on recursive Lasso to monitor the flowmeter online, and timely Find faults and reduce losses.
本发明的基于递推Lasso的流量计异常在线监测方法用于监测烟草制丝过程增湿水流量是否异常,如图1所示,本实施例提供的基于递推Lasso的流量计异常在线监测方法在实际执行过程中,具体包括如下步骤:The online monitoring method for flowmeter abnormality based on recursive Lasso of the present invention is used to monitor whether the flow rate of humidifying water in the tobacco shred process is abnormal. As shown in Figure 1, the online monitoring method for abnormal flowmeter based on recursive Lasso provided by this embodiment In the actual implementation process, the specific steps are as follows:
步骤S1、根据流量计的历史流量数据,建立自回归模型,并通过离线Lasso算法确定所述自回归模型的模型系数。Step S1. Establish an autoregressive model according to the historical flow data of the flowmeter, and determine the model coefficients of the autoregressive model through an offline Lasso algorithm.
在本发明的基于递推Lasso的流量计异常在线监测方法的一种实施方式中,所述步骤S1具体可以包括:In one embodiment of the online monitoring method for flowmeter abnormalities based on recursive Lasso of the present invention, the step S1 may specifically include:
步骤S11、采集流量计流量的数据样本,得到历史流量数据集合。Step S11, collecting data samples of the flow rate of the flow meter to obtain a historical flow data set.
具体地,按预设时间长度,收集在梗丝加料单元中增湿水流量计所采集到的数据作为在线监测的数据源,以得到历史流量数据集合。Specifically, the data collected by the humidification water flow meter in the cut stem feeding unit is collected according to a preset time length as a data source for online monitoring, so as to obtain a set of historical flow data.
步骤S12、根据所述历史流量数据集合,构造自回归模型。Step S12, constructing an autoregressive model according to the historical flow data set.
步骤S13、基于离线Lasso算法求取所述自回归模型的回归系数。Step S13, calculating the regression coefficient of the autoregressive model based on the offline Lasso algorithm.
步骤S14、利用交叉验证的方法确定所述自回归模型的的正则化参数。Step S14, using a cross-validation method to determine the regularization parameters of the autoregressive model.
具体而言,通过以下公式确定所述离线Lasso算法的最优条件,并根据所述离线Lasso算法的最优条件确定所述自回归模型的回归系数和正则化参数:Specifically, the optimal condition of the offline Lasso algorithm is determined by the following formula, and the regression coefficient and the regularization parameter of the autoregressive model are determined according to the optimal condition of the offline Lasso algorithm:
其中,p为自回归模型的滞后系数,yn表示当前时刻待预测的流量值,yn-i表示历史时刻的流量值,α表示自回归模型的回归系数,μn-1表示自回归模型的正则化参数。Among them, p is the lag coefficient of the autoregressive model, y n represents the flow value to be predicted at the current moment, y ni represents the flow value at the historical moment, α represents the regression coefficient of the autoregressive model, μ n-1 represents the regularity of the autoregressive model parameterization.
将l1范数引入自回归模型中,假设流量计测得的流量数据为y∈Rn-1,模型的回归系数为p,其中εi表示流量计测量的噪声数据。那么Lasso的最优解就变为:The l 1 norm is introduced into the autoregressive model, assuming that the flow data measured by the flowmeter is y∈R n-1 , and the regression coefficient of the model is p, Where εi represents the noise data measured by the flowmeter. Then the optimal solution of Lasso becomes:
其中μn-1是正则化参数,公式(1)解的形式是稀疏的,只有少量的元素是非零的,Lasso算法的目的是为了在变量数较多的情况下,通过把一些不重要的变量置零来突出一些非零的重要变量,这些非零变量则是对预测有用的变量,因此可以根据这些非零变量来确定哪些历史时刻的流量对预测当前时刻的流量有帮助。Among them, μ n-1 is the regularization parameter. The form of the solution of formula (1) is sparse, and only a small number of elements are non-zero. The purpose of the Lasso algorithm is to combine some unimportant Variables are set to zero to highlight some non-zero important variables, which are useful variables for prediction, so it can be determined which historical moments of flow are helpful to predict the current moment of flow based on these non-zero variables.
将回归系数向量α中非零元素的索引定义一个“有效集”,为了使得符号变得简便,将“有效集”放在前面,例如αT=(α1 T,0T)、vT=(v1 T,v2 T)。其中,对于所有的i来说,满足越阶函数v1i=sgn(α1i);对于所有的j来说,-1≤v2j≤1,其中,i和j是针对不同的集合分开表示的索引,这样可以避免混淆。同时根据“有效集”,将Y划分成Y=(Y1 Y2),其中Y表示公式(1)中的yn-1到yn-p的集合,每个集合中的元素都对应了一个回归系数α,Y1表示回归系数α不为零的对应Y中的元素,Y2是指表示回归系数α为零的对应Y中的元素。如果求得的解是唯一的,那么Y1 TY1是可逆的,那么,最优条件为The index of the non-zero elements in the regression coefficient vector α defines an "effective set". In order to make the notation easier, the "effective set" is placed in front, for example, α T =(α 1 T ,0 T ), v T = (v 1 T ,v 2 T ). Among them, for all i, the higher-order function v 1i =sgn(α 1i ) is satisfied; for all j, -1≤v 2j ≤1, where i and j are represented separately for different sets index, which avoids confusion. At the same time, according to the "effective set", Y is divided into Y=(Y 1 Y 2 ), where Y represents the set of y n-1 to y np in formula (1), and the elements in each set correspond to a regression The coefficient α, Y1 means the element in Y corresponding to which the regression coefficient α is not zero, and Y2 refers to the element in Y corresponding to which the regression coefficient α is zero. If the obtained solution is unique, then Y 1 T Y 1 is reversible, then the optimal condition is
需要说明的是,在知道“有效集”和特征向量α内系数的符号的情况下,就可以以封闭的形式计算α。It should be noted that α can be computed in closed form given the “active set” and the sign of the coefficients in the eigenvector α.
步骤S2、利用流量计的实时流量,通过递推Lasso(Least absolute shrinkageand selection operator,最小绝对收敛和选择算子)算法对所述自回归模型进行更新。Step S2, using the real-time flow rate of the flowmeter to update the autoregressive model through a recursive Lasso (Least absolute shrinkage and selection operator) algorithm.
根据读入的新数据来微调自回归模型的模型参数,可以提高自回归模型的预测效果。在本发明的基于递推Lasso的流量计异常在线监测方法的一种实施方式中,所述步骤S2具体可以包括:Fine-tuning the model parameters of the autoregressive model according to the new data read in can improve the predictive effect of the autoregressive model. In one embodiment of the online monitoring method for flowmeter abnormality based on recursive Lasso of the present invention, the step S2 may specifically include:
步骤S21、将通过离线Lasso算法所确定的所述自回归模型的模型系数作为所述自回归模型在更新过程中的初始模型系数。Step S21, using the model coefficients of the autoregressive model determined by the offline Lasso algorithm as the initial model coefficients of the autoregressive model in the updating process.
步骤S22、每读入一个新的流量计实时流量,则根据该新的流量计实时流量对所述自回归模型的模型参数进行更新,并且通过以下公式对所述自回归模型进行更新:Step S22, each time a new real-time flow rate of the flowmeter is read, the model parameters of the autoregressive model are updated according to the new real-time flow rate of the flowmeter, and the autoregressive model is updated by the following formula:
其中,Yn-1表示流量计测得流量的历史流量数据集合,Zn-1表示自回归模型的自回归项,yn表示当前时刻待预测的流量值,Zn表示n时刻对应的自回归项,t为递推Lasso算法的参数,其取值为0-1,在读入新数据后,t的值都会从0变到1,当t=1时表示新数据完全读入,自回归模型内的参数α和μ也更新完毕,准备读入下一个时刻的新数据,当下一时刻的新数据读入之后,t的值又从0开始变到1。Among them, Y n-1 represents the historical flow data collection of flow measured by the flowmeter, Z n-1 represents the autoregressive item of the autoregressive model, y n represents the flow value to be predicted at the current moment, and Z n represents the autoregressive value corresponding to time n. Regression item, t is the parameter of the recursive Lasso algorithm, and its value is 0-1. After the new data is read in, the value of t will change from 0 to 1. When t=1, it means that the new data is completely read in. The parameters α and μ in the regression model are also updated, and the new data at the next moment is ready to be read in. After the new data at the next moment is read in, the value of t changes from 0 to 1 again.
具体地,假设已经计算出在n-1时刻公式(1)的解α(n-1)后,得到了新的观测数据yn,那么就可以利用yn和历史数据去预测n+1时刻的数据yn+1,那么目标就改为计算α(n)的值。由此引出递推Lasso的解法。记zn=(yn,yn-1,...,yn-p),Zn=(z1,z2,...,zn)T,Yn-1=(y1,y2,...yn-1)T,其中,zn表示预测当前流量所需要用到的历史时刻流量,Zn表示之前所有历史时刻流量的集合,Yn-1表示所有历史时刻流量,那么优化的目标函数就变为:Specifically, assuming that the solution α (n-1) of formula (1) at time n-1 has been calculated, new observation data y n is obtained, then y n and historical data can be used to predict time n+1 data y n+1 , then the goal is to calculate the value of α (n) . This leads to the solution of recursive Lasso. Denote z n =(y n ,y n-1 ,...,y np ), Z n =(z 1 ,z 2 ,...,z n ) T , Y n-1 =(y 1 ,y 2 ,...y n-1 ) T , where z n represents the historical time flow required to predict the current flow, Z n represents the collection of all previous historical time flow, Y n-1 represents all historical time flow, Then the optimized objective function becomes:
整个更新路径就变成α(n-1)=α(0,μn-1)到α(n)=α(1,μn),根据这个更新路径,更新方法可以分为两步:The entire update path becomes α (n-1) = α(0,μ n-1 ) to α (n) =α(1,μ n ). According to this update path, the update method can be divided into two steps:
第一步:当t=0时,从μn-1到μn更新正则化参数,这相当于采用最小角回归的方法计算两者之间的正则化路径。Step 1: when t=0, update the regularization parameters from μ n−1 to μ n , which is equivalent to calculating the regularization path between the two using the minimum angle regression method.
第二步:当μ=μn时,计算t从0到1。The second step: when μ=μ n , calculate t from 0 to 1.
上述第二步的求解过程为:首先需要证明的是α(t,μ)是针对t的分段光滑函数。为了使得符号更加简便,令α(t)=α(t,μ),若已经知道“有效集”以及α内系数的符号,就能计算出Lasso的解。在第一步中,已经计算出“有效集”以及α内系数的符号,并且当t∈[0,t*)时(其中,t*表示在t∈[0,t*)时α(t)是一条光滑的曲线),“有效集”以及α内系数的符号保持不变,Lasso的解α(t)是光滑的。将“有效集”改变的一个点称为转折点,接下来要分析怎么去计算这个点:The solution process of the above second step is as follows: firstly, it needs to be proved that α(t, μ) is a piecewise smooth function for t. In order to make the sign more convenient, let α(t)=α(t,μ), if you already know the "effective set" and the sign of the coefficient in α, you can calculate the solution of Lasso. In the first step, the "effective set" and the signs of the coefficients in α have been calculated, and when t∈[0,t * ) (where t * means at t∈[0,t * ) when α(t ) is a smooth curve), the "effective set" and the signs of the coefficients in α remain unchanged, and Lasso's solution α(t) is smooth. A point where the "effective set" changes is called a turning point. Next, we will analyze how to calculate this point:
当t=t*时,更新“有效集”以及α内系数的符号,并在下一个转折点到达之前保持不变,不断迭代这个过程,直到t=1,这样就可以计算出所需要的解α(t)。When t=t * , update the "effective set" and the sign of the coefficient in α, and keep it unchanged until the next turning point arrives, and iterate this process until t=1, so that the required solution α(t ).
从而,得到了添加观测值的Lasso的在线更新算法:Thus, the online update algorithm of Lasso with added observations is obtained:
STEP1:计算α(n-1)=α(0,μn-1)到α(n)=α(1,μn);STEP1: Calculate α (n-1) = α(0, μ n-1 ) to α (n) = α(1, μ n );
STEP2:初始化α(0,μn)的非零系数到“有效集”,令v=sgn(α(0,μn)),令v1和zn,1是v和zn根据“有效集”划分的子向量,是的子矩阵,其中的列是“有效集”,初始化初始化转折点t'=0;STEP2: Initialize the non-zero coefficients of α(0,μ n ) to the "effective set", let v=sgn(α(0,μ n )), let v 1 and z n,1 be v and z n according to the "effective set"Set" sub-vectors, yes A submatrix of , the columns of which are the "active set", initialized Initialize turning point t'=0;
STEP3:计算下一个转折点t'。如果该转折点小于之前的转折点或者转折点大于1,跳转到STEP5,STEP3: Calculate the next turning point t'. If the turning point is smaller than the previous turning point or the turning point is greater than 1, jump to STEP5,
第一种情形是:先α1(t')中第i个元素变成0;然后从“有效集”中移除i;再将vi置0;The first case is: first the i-th element in α 1 (t') becomes 0; then remove i from the "effective set"; then set v i to 0;
第二种情形是:先将ω2(t')的第j个元素的绝对值到达1;然后将j加入“有效集”;接着,如果该元素到达1(或-1),那么将vj置1(或-1)。The second case is: first, the absolute value of the jth element of ω 2 (t') reaches 1; then j is added to the "effective set"; then, if the element reaches 1 (or -1), then v j is set to 1 (or -1).
STEP4:根据更新后的“有效集”更新v1,和zn,1,更新 STEP4: Update v 1 according to the updated "effective set", and z n,1 , update
STEP5:计算当t=1时最后的结果,其中α(n)的值由的有效集给出。STEP5: Calculate the final result when t=1, where the value of α (n) is determined by The effective set of is given.
当数据集中观测值过多就会使得STEP1的计算时间变长,且在后续计算的过程中,距离当前时刻较为久远的观测点数据对预测当前时刻的值影响不大,因此考虑在一段时间后剔除最初的历史数据,使得模型的计算时间稳定在一定范围之内。When there are too many observations in the data set, the calculation time of STEP1 will be longer, and in the subsequent calculation process, the observation point data that is farther away from the current time has little effect on predicting the value at the current time, so consider after a period of time Eliminate the initial historical data, so that the calculation time of the model is stable within a certain range.
假设已经计算出在n时刻的解α(n)后,在新的观测数据读入之前,需要剔除最初的数据,由此引出递推Lasso剔除历史数据的解法。记z1=(y0,y-1,...,y-p),Z=(z2,...,zn)T,Y=(y2,y3,...,yn)T,那么优化的目标函数就变为:Assuming that after the solution α (n) at time n has been calculated, the original data needs to be eliminated before the new observation data is read in, thus leading to the solution of recursive Lasso elimination of historical data. Write z 1 =(y 0 ,y -1 ,...,y -p ), Z=(z 2 ,...,z n ) T , Y=(y 2 ,y 3 ,...,y n ) T , then the optimized objective function becomes:
整个更新路径就变成α(n)=α(1,μn)到α(n')=α(0,μn'),根据这个更新路径,更新方法可以分为两步:STEP1:当t=1时,从μn到μn'更新正则化参数,这相当于采用最小角回归的方法计算两者之间的正则化路径。The entire update path becomes α (n) = α(1,μ n ) to α (n') = α(0,μ n' ), according to this update path, the update method can be divided into two steps: STEP1: When When t=1, the regularization parameters are updated from μ n to μ n ', which is equivalent to calculating the regularization path between the two using the minimum angle regression method.
STEP2:当μ=μn'时,计算t从1到0。STEP2: When μ=μ n' , calculate t from 1 to 0.
此时STEP2的计算步骤与添加观测数据的步骤相同。At this time, the calculation steps of STEP2 are the same as the steps of adding observation data.
步骤S3、根据基于更新后的所述自回归模型得到的预测流量和流量计在下一时刻的实际流量确定流量计是否异常。Step S3. Determine whether the flowmeter is abnormal according to the predicted flow obtained based on the updated autoregressive model and the actual flow of the flowmeter at the next moment.
在本发明的基于递推Lasso的流量计异常在线监测方法的一种实施方式中,所述步骤S3具体可以包括:In one embodiment of the online monitoring method for flowmeter abnormality based on recursive Lasso of the present invention, the step S3 may specifically include:
步骤S31、根据更新后的所述自回归模型,实时预测流量计下一时刻的流量,得到预测流量。Step S31 , according to the updated autoregressive model, predict the flow rate of the flow meter at the next moment in real time, and obtain the predicted flow rate.
步骤S32、根据得到的不同时刻对应的预测流量,得到预测流量的估计曲线图。Step S32 , according to the obtained predicted flow corresponding to different time points, an estimated graph of the predicted flow is obtained.
如图2的虚线为预测流量的估计曲线图。The dotted line in Fig. 2 is the estimated curve diagram of the predicted flow.
步骤S33、根据所述预测流量和流量计在下一时刻的实际流量的偏差,得到残差图。Step S33 , according to the deviation between the predicted flow rate and the actual flow rate of the flowmeter at the next moment, a residual map is obtained.
这里的偏差为预测值与实际流量差值的绝对值。其中,残差图用于表征流量在预定的使用环境下是否发生异常。本发明在一些实施方式中,得到的残差图如图3所示,图中直线从上到下依次为50kg/h、25kg/h、10kg/h。The deviation here is the absolute value of the difference between the predicted value and the actual flow. Among them, the residual graph is used to characterize whether the flow is abnormal in a predetermined usage environment. In some embodiments of the present invention, the obtained residual diagram is shown in Figure 3, and the straight lines in the diagram are 50kg/h, 25kg/h, and 10kg/h from top to bottom.
步骤S34、根据所述残差图确定流量计是否异常。Step S34. Determine whether the flowmeter is abnormal according to the residual map.
利用偏差值对在线流量工作过程中的变化进行估计,当偏差值过大时就认定流量出现异常操作。在本发明的基于递推Lasso的流量计异常在线监测方法的一种实施方式中,所述步骤S34具体可以包括:Use the deviation value to estimate the change in the online flow working process, and when the deviation value is too large, it is determined that the flow has abnormal operation. In one embodiment of the online monitoring method for flowmeter abnormality based on recursive Lasso of the present invention, the step S34 may specifically include:
步骤S341、在所述残差图的纵坐标超过预设阈值时,则确定流量计的流量数据异常。Step S341, when the ordinate of the residual graph exceeds a preset threshold, it is determined that the flow data of the flow meter is abnormal.
在图3中,当第7000个时刻开始流量计出现故障后,预测值与实际值偏差较大,残差(残差图的纵坐标)大于50kg/h,从而被认定为流量数据异常(图中的异常点在100kg/h附近)。In Figure 3, when the flowmeter fails at the 7000th moment, the predicted value deviates greatly from the actual value, and the residual error (the vertical axis of the residual error graph) is greater than 50 kg/h, which is considered as abnormal flow data (Fig. The abnormal point in is around 100kg/h).
步骤S342、根据生产记录确定导致流量计的流量数据异常的原因是流量计出现异常还是梗丝产品批次发生更换。Step S342 , according to the production records, it is determined whether the abnormality of the flow data of the flowmeter is caused by the abnormality of the flowmeter or the replacement of the cut stem product batch.
若梗丝产品批次发生更换,则确定导致流量计的流量数据异常的原因是生产进行了调整;若梗丝产品批次未发生更换,则确定导致流量计的流量数据异常的原因是流量计出现异常,此时需要对流量计及时进行维护。在本发明中,利用残差图实时跟踪流量计的状态性能,当流量出现异常后,能够及时监测并发出警报。If the cut stem product batch is changed, it is determined that the reason for the abnormal flow data of the flowmeter is that the production has been adjusted; if the cut stem product batch has not been changed, it is determined that the cause of the abnormal flow data of the flowmeter is the If an abnormality occurs, the flowmeter needs to be maintained in time. In the present invention, the status performance of the flow meter is tracked in real time by using the residual graph, and when the flow is abnormal, it can be monitored in time and an alarm can be issued.
进一步地,在本发明的一些实施方式中,所述基于递推Lasso的流量计异常在线监测方法,还包括:Further, in some embodiments of the present invention, the recursive Lasso-based flowmeter abnormality online monitoring method also includes:
步骤S4、基于预设的梗丝产品批次的流量阈值和所述预测流量,确定当前时刻的梗丝产品规格。Step S4, based on the preset flow threshold of the cut-stem product batch and the predicted flow, determine the specification of the cut-stem product at the current moment.
基于预设的梗丝产品批次的流量阈值以及基于多个时刻的预测流量所绘制的估计曲线图,确定当前时刻进入梗丝加料单元的产品批次。由于每个梗丝产品批次不同,增湿水流量也会不同,因此通过设定增湿水流量的范围,则可根据自回归模型估计的流量大小来判断出梗丝产品规格。Based on the preset flow threshold of the cut-stem product batch and the estimated curve drawn based on the predicted flow at multiple times, the product batch entering the cut-stem feeding unit at the current moment is determined. Since each batch of cut stem products is different, the flow rate of humidification water will also be different. Therefore, by setting the range of flow rate of humidification water, the specification of cut stem products can be judged according to the flow rate estimated by the autoregressive model.
数据样品中的梗丝分为两种,一种是普通梗,另一种是黄金叶专用梗。当物料供给稳定后,普通梗的增湿水流量在110~130kg/h的范围内,黄金叶专用梗的增湿水流量在140~170kg/h的范围内。因此可以根据预测的流量值来判断梗丝的产品规格。在图2中,虚线为预测值,实曲线为实际值;实直线为140kg/h的划分线,若物料供给稳定后,预测流量小于140kg/h则判断为普通梗,测流量大于140kg/h则判断为黄金叶专用梗。The shredded stems in the data samples are divided into two types, one is ordinary stems, and the other is special golden leaf stems. When the material supply is stable, the humidifying water flow rate of common stalks is within the range of 110-130kg/h, and the humidifying water flow rate of special golden leaf stalks is within the range of 140-170kg/h. Therefore, the product specifications of the shredded stems can be judged according to the predicted flow value. In Figure 2, the dotted line is the predicted value, and the solid curve is the actual value; the solid line is the dividing line of 140kg/h. If the material supply is stable and the predicted flow rate is less than 140kg/h, it is judged as a common stem, and the measured flow rate is greater than 140kg/h. It is judged as a special stalk for golden leaves.
本发明实施例提供的基于递推Lasso的流量计异常在线监测方法,通过递推Lasso算法对自回归模型进行更新,根据基于更新后的自回归模型得到的预测流量和实际流量来监测流量计的异常情况,具有准确度高、操作方便和实时跟踪等特点,为烟草制丝过程增湿水流量的在线监测提供了科学、客观、可靠的技术支持,进而能够保证流量计的计量性能稳定性,并且本发明也适用于其他计量器具的检测场景中,应用前景广泛。The recursive Lasso-based flowmeter abnormality online monitoring method provided by the embodiment of the present invention updates the autoregressive model through the recursive Lasso algorithm, and monitors the flowmeter according to the predicted flow rate and actual flow rate obtained based on the updated autoregressive model. Abnormal conditions, with the characteristics of high accuracy, convenient operation and real-time tracking, provide scientific, objective and reliable technical support for the on-line monitoring of humidification water flow in the tobacco shred process, thereby ensuring the stability of the measurement performance of the flowmeter. Moreover, the present invention is also applicable to detection scenarios of other measuring instruments, and has wide application prospects.
相应地,如图4所示,本发明还提供一种基于递推Lasso的流量计异常在线监测系统,包括:Correspondingly, as shown in Fig. 4, the present invention also provides a flowmeter abnormality online monitoring system based on recursive Lasso, including:
自回归模型建立模块1,用于根据流量计的历史流量数据,建立自回归模型,并通过离线Lasso算法确定所述自回归模型的模型系数;The autoregressive model building module 1 is used to set up an autoregressive model according to the historical flow data of the flowmeter, and determine the model coefficients of the autoregressive model by an off-line Lasso algorithm;
自回归模型更新模块2,用于利用流量计的实时流量,通过递推Lasso算法对所述自回归模型进行更新;The autoregressive
流量监测模块3,用于根据基于更新后的所述自回归模型得到的预测流量和流量计在下一时刻的实际流量确定流量计是否异常。The
本发明实施例提供的基于递推Lasso的流量计异常在线监测系统,利用自回归模型更新模块通过递推Lasso算法对自回归模型进行更新,流量监测模块根据基于更新后的自回归模型得到的预测流量和实际流量来监测流量计的异常情况,具有准确度高、操作方便和实时跟踪等特点,为烟草制丝过程增湿水流量的在线监测提供了科学、客观、可靠的技术支持,进而能够保证流量计的计量性能稳定性,并且本发明也适用于其他计量器具的检测场景中,应用前景广泛。The recursive Lasso-based flowmeter abnormality online monitoring system provided by the embodiment of the present invention uses the autoregressive model update module to update the autoregressive model through the recursive Lasso algorithm, and the flow monitoring module is based on the prediction obtained based on the updated autoregressive model It has the characteristics of high accuracy, convenient operation and real-time tracking. It provides scientific, objective and reliable technical support for the on-line monitoring of humidification water flow in the tobacco shredded process, and can further The metering performance stability of the flowmeter is guaranteed, and the present invention is also applicable to detection scenarios of other metering instruments, and has broad application prospects.
应理解图4所示的基于递推Lasso的流量计异常在线监测系统的各个部件的划分仅仅是一种逻辑功能的划分,实际实现时可以全部或部分集成到一个物理实体上,也可以物理上分开。且这些部件可以全部以软件通过处理元件调用的形式实现;也可以全部以硬件的形式实现;还可以部分部件以软件通过处理元件调用的形式实现,部分部件通过硬件的形式实现。例如,某个上述模块可以为单独设立的处理元件,也可以集成在电子设备的某一个芯片中实现。其它部件的实现与之类似。此外这些部件全部或部分可以集成在一起,也可以独立实现。在实现过程中,上述方法的各步骤或以上各个部件可以通过处理器元件中的硬件的集成逻辑电路或者软件形式的指令完成。It should be understood that the division of the various components of the recursive Lasso-based flowmeter abnormality online monitoring system shown in Figure 4 is only a division of logical functions. In actual implementation, it can be fully or partially integrated into a physical entity, or physically separate. And these components can all be implemented in the form of software called by the processing element; they can also be implemented in the form of hardware; some components can also be implemented in the form of software called by the processing element, and some components can be implemented in the form of hardware. For example, a certain above-mentioned module may be a separately established processing element, or may be integrated into a certain chip of the electronic device for implementation. The implementation of other components is similar. In addition, all or part of these components can be integrated together, or implemented independently. In the process of implementation, each step of the above method or each of the above components can be completed by an integrated logic circuit of hardware in the processor element or an instruction in the form of software.
本发明还提供一种计算机可读存储介质,所述计算机可读存储介质中存储有计算机程序,当其在计算机上运行时,使得计算机执行上述的基于递推Lasso的流量计异常在线监测方法。The present invention also provides a computer-readable storage medium, wherein a computer program is stored in the computer-readable storage medium, and when it is run on a computer, the computer executes the above-mentioned recursive Lasso-based flowmeter abnormality online monitoring method.
本发明还提供一种计算机程序产品,其特征在于,所述计算机程序产品在终端设备上运行时,使所述终端设备执行上述的基于递推Lasso的流量计异常在线监测方法。The present invention also provides a computer program product, which is characterized in that, when the computer program product is run on a terminal device, the terminal device is made to execute the above-mentioned recursive Lasso-based flowmeter abnormality online monitoring method.
在上述实施例中,可以全部或部分地通过软件、硬件、固件或者其任意结合来实现。当使用软件实现时,可以全部或部分地以计算机程序产品的形式实现。所述计算机程序产品包括一个或多个计算机指令。在计算机上加载和执行所述计算机指令时,全部或部分地产生按照本申请实施例所述的流程或功能。所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程设备。所述计算机指令可以存储在计算机可读存储介质中,或者从一个计算机可读存储介质向另一个计算机可读存储介质传输,例如,所述计算机指令可以从一个网站站点、计算机、服务器或数据中心通过有线(例如:同轴电缆、光纤、数据用户线(Digital Subscriber Line,DSL))或无线(例如:红外、无线、微波等)方式向另一个网站站点、计算机、服务器或数据中心进行传输。所述计算机可读存储介质可以是计算机能够存取的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。所述可用介质可以是磁性介质(例如:软盘、硬盘、磁带)、光介质(例如:数字通用光盘(Digital Versatile Disc,DVD))、或者半导体介质(例如:固态硬盘(Solid State Disk,SSD))等。In the above embodiments, all or part may be implemented by software, hardware, firmware or any combination thereof. When implemented using software, it may be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer instructions. When the computer instructions are loaded and executed on the computer, the processes or functions according to the embodiments of the present application will be generated in whole or in part. The computer can be a general purpose computer, a special purpose computer, a computer network, or other programmable devices. The computer instructions may be stored in or transmitted from one computer-readable storage medium to another computer-readable storage medium, for example, the computer instructions may be transmitted from a website, computer, server or data center Transmission to another website site, computer, server or data center via wired (eg coaxial cable, optical fiber, Digital Subscriber Line (DSL)) or wireless (eg infrared, wireless, microwave, etc.). The computer-readable storage medium may be any available medium that can be accessed by a computer, or a data storage device such as a server or a data center integrated with one or more available media. The available medium may be a magnetic medium (for example: floppy disk, hard disk, magnetic tape), an optical medium (for example: Digital Versatile Disc (Digital Versatile Disc, DVD)), or a semiconductor medium (for example: Solid State Disk (Solid State Disk, SSD) )Wait.
至此,已经详细描述了本公开的各实施例。为了避免遮蔽本公开的构思,没有描述本领域所公知的一些细节。本领域技术人员根据上面的描述,完全可以明白如何实施这里公开的技术方案。So far, the embodiments of the present disclosure have been described in detail. Certain details known in the art have not been described in order to avoid obscuring the concept of the present disclosure. Based on the above description, those skilled in the art can fully understand how to implement the technical solutions disclosed herein.
虽然已经通过示例对本公开的一些特定实施例进行了详细说明,但是本领域的技术人员应该理解,以上示例仅是为了进行说明,而不是为了限制本公开的范围。本领域的技术人员应该理解,可在不脱离本公开的范围和精神的情况下,对以上实施例进行修改或者对部分技术特征进行等同替换。本公开的范围由所附权利要求来限定。Although some specific embodiments of the present disclosure have been described in detail through examples, those skilled in the art should understand that the above examples are for illustration only, rather than limiting the scope of the present disclosure. Those skilled in the art should understand that the above embodiments can be modified or some technical features can be equivalently replaced without departing from the scope and spirit of the present disclosure. The scope of the present disclosure is defined by the appended claims.
Claims (5)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110765956.0A CN113468751B (en) | 2021-07-05 | 2021-07-05 | Recursion Lasso-based flowmeter anomaly online monitoring method and system and storage medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110765956.0A CN113468751B (en) | 2021-07-05 | 2021-07-05 | Recursion Lasso-based flowmeter anomaly online monitoring method and system and storage medium |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113468751A CN113468751A (en) | 2021-10-01 |
CN113468751B true CN113468751B (en) | 2022-12-27 |
Family
ID=77878744
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110765956.0A Active CN113468751B (en) | 2021-07-05 | 2021-07-05 | Recursion Lasso-based flowmeter anomaly online monitoring method and system and storage medium |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113468751B (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105975443A (en) * | 2016-05-04 | 2016-09-28 | 西南大学 | Lasso-based anomaly detection method and system |
CN109443419A (en) * | 2018-08-31 | 2019-03-08 | 广州市世科高新技术有限公司 | A kind of rectifier on-line monitoring method based on machine learning |
CN111628961A (en) * | 2020-03-30 | 2020-09-04 | 西安交大捷普网络科技有限公司 | DNS (Domain name Server) anomaly detection method |
CN111737249A (en) * | 2020-08-24 | 2020-10-02 | 国网浙江省电力有限公司 | Abnormal data detection method and device based on Lasso algorithm |
CN112631250A (en) * | 2020-12-15 | 2021-04-09 | 中国计量大学 | Fault isolation and identification method in nonlinear process based on denoising autoencoder |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10068204B2 (en) * | 2014-07-23 | 2018-09-04 | International Business Machines Corporation | Modeling and visualizing a dynamic interpersonal relationship from social media |
CN107610464B (en) * | 2017-08-11 | 2019-02-12 | 河海大学 | A Trajectory Prediction Method Based on Gaussian Mixture Time Series Model |
US11379284B2 (en) * | 2018-03-13 | 2022-07-05 | Nec Corporation | Topology-inspired neural network autoencoding for electronic system fault detection |
-
2021
- 2021-07-05 CN CN202110765956.0A patent/CN113468751B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105975443A (en) * | 2016-05-04 | 2016-09-28 | 西南大学 | Lasso-based anomaly detection method and system |
CN109443419A (en) * | 2018-08-31 | 2019-03-08 | 广州市世科高新技术有限公司 | A kind of rectifier on-line monitoring method based on machine learning |
CN111628961A (en) * | 2020-03-30 | 2020-09-04 | 西安交大捷普网络科技有限公司 | DNS (Domain name Server) anomaly detection method |
CN111737249A (en) * | 2020-08-24 | 2020-10-02 | 国网浙江省电力有限公司 | Abnormal data detection method and device based on Lasso algorithm |
CN112631250A (en) * | 2020-12-15 | 2021-04-09 | 中国计量大学 | Fault isolation and identification method in nonlinear process based on denoising autoencoder |
Non-Patent Citations (2)
Title |
---|
PEMOGEN:automatic adaptive performance modeling during program runtime;PEMOGEN automatic adaptive performance modeling during program;《PACT "14: Proceedings of the 23rd international conference on Parallel architectures and compilation》;20140831;393-404页 * |
基于向量自回归模型的移动通信基站流量预测;何勇 等;《工业工程与管理》;20171231;第22卷(第04期);79-84页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113468751A (en) | 2021-10-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110375787B (en) | An online evaluation method for the running state of a metering instrument | |
CN111123188A (en) | Comprehensive verification method and system of electric energy meter based on improved least squares method | |
CN117606578A (en) | An intelligent gas flow monitoring method | |
CN112284440B (en) | An Adaptive Correction Method for Sensor Data Bias | |
CN117871771B (en) | Big data-based gas energy monitoring method | |
CN118211837A (en) | A performance evaluation method and device for smart electric meter | |
CN118689866B (en) | A data optimization method and system for ultrasonic flowmeter | |
CN116292367B (en) | Power plant fan system abnormal condition detection method based on one-dimensional convolution | |
CN116755411A (en) | Industrial equipment data acquisition method and device | |
CN113468751B (en) | Recursion Lasso-based flowmeter anomaly online monitoring method and system and storage medium | |
CN110781206A (en) | Method for predicting whether electric energy meter in operation fails or not by learning meter-dismantling and returning failure characteristic rule | |
CN118466648A (en) | A constant temperature and humidity control method and system for a cell culture box | |
CN110914776B (en) | Apparatus and method for identifying effects and causes on variability of model-based controller performance or control deviations | |
CN117691736A (en) | Non-invasive intelligent factory energy consumption monitoring method and system | |
CN105117606A (en) | Method for determining element fault probability change tendency | |
WO2021140942A1 (en) | Diagnosing device, diagnosing method, and program | |
CN115307670B (en) | Method, apparatus and medium for locating central air conditioning system anomaly sensors | |
Kučera et al. | The digital twin of a measuring process within the Industry 4.0 concept | |
CN112507286A (en) | State-based equipment fault probability analysis method and device and electronic equipment | |
CN118348341B (en) | Data model for intelligent monitoring of operation state of substation lightning arrester | |
CN119147141B (en) | Substation SF6 pressure gauge based data acquisition and processing method | |
CN118482817B (en) | Digital textile equipment management method and system based on data analysis | |
CN119468584A (en) | A circulating cooling water system integrating big data and artificial intelligence | |
Tong et al. | A novel adaptive extreme learning machine-autoencoder algorithm for online monitoring and fault detection | |
CN117010213A (en) | Nuclear sensor reliability assessment method based on operation data and storage medium |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20240912 Address after: No.16 Yulin South Road, Zhengdong New District, Zhengzhou City, Henan Province Patentee after: CHINA TOBACCO HENAN INDUSTRIAL Co.,Ltd. Country or region after: China Address before: No.16 Yulin South Road, Zhengdong New District, Zhengzhou City, Henan Province Patentee before: CHINA TOBACCO HENAN INDUSTRIAL Co.,Ltd. Country or region before: China Patentee before: China Jiliang University |
|
TR01 | Transfer of patent right |