[go: up one dir, main page]

CN116822999A - Method and system for predicting monitoring density of oil product of oil mixing interface of finished oil pipeline - Google Patents

Method and system for predicting monitoring density of oil product of oil mixing interface of finished oil pipeline Download PDF

Info

Publication number
CN116822999A
CN116822999A CN202311107413.5A CN202311107413A CN116822999A CN 116822999 A CN116822999 A CN 116822999A CN 202311107413 A CN202311107413 A CN 202311107413A CN 116822999 A CN116822999 A CN 116822999A
Authority
CN
China
Prior art keywords
oil
monitoring
density
downstream
upstream
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.)
Granted
Application number
CN202311107413.5A
Other languages
Chinese (zh)
Other versions
CN116822999B (en
Inventor
袁子云
陈雷
刘刚
姬浩洋
邵伟明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum East China
Original Assignee
China University of Petroleum East China
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 China University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN202311107413.5A priority Critical patent/CN116822999B/en
Publication of CN116822999A publication Critical patent/CN116822999A/en
Application granted granted Critical
Publication of CN116822999B publication Critical patent/CN116822999B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N9/00Investigating density or specific gravity of materials; Analysing materials by determining density or specific gravity
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/25Fusion techniques
    • G06F18/254Fusion techniques of classification results, e.g. of results related to same input data
    • G06F18/256Fusion techniques of classification results, e.g. of results related to same input data of results relating to different input data, e.g. multimodal recognition
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/27Regression, e.g. linear or logistic regression
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0637Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Forestry; Mining
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Data Mining & Analysis (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Human Resources & Organizations (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Strategic Management (AREA)
  • General Engineering & Computer Science (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Economics (AREA)
  • Evolutionary Biology (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Health & Medical Sciences (AREA)
  • Marketing (AREA)
  • General Health & Medical Sciences (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Educational Administration (AREA)
  • Immunology (AREA)
  • Quality & Reliability (AREA)
  • Animal Husbandry (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Mining & Mineral Resources (AREA)
  • Primary Health Care (AREA)
  • Biochemistry (AREA)
  • Agronomy & Crop Science (AREA)
  • Pathology (AREA)
  • Analytical Chemistry (AREA)
  • Operations Research (AREA)
  • Game Theory and Decision Science (AREA)
  • Development Economics (AREA)
  • Chemical & Material Sciences (AREA)
  • Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)
  • Pipeline Systems (AREA)

Abstract

本发明属于数据处理技术领域,为解决无法对后行油品监测密度的可靠预测的问题,提出了成品油管道混油界面后行油品监测密度预测方法及系统,选取历史数据中成品油管道上下游站场水力热力信息,前后行油品上游站场监测密度以及前行油品上下游站场监测密度差值作为关键输入特征变量,选取后行油品上下游站场监测密度差值作为输出变量,利用高斯混合回归算法对关键输入特征变量和输出变量进行多模态识别,根据模态识别结果,利用最大期望算法进行训练得到模态对应的后行油品上下游站场监测密度差值修正模型;并结合已获取的成品油管道后行油品上游站场监测密度值,得到预测结果。可以准确预测混油界面后行油品监测密度。

The invention belongs to the field of data processing technology. In order to solve the problem of being unable to reliably predict the density of subsequent oil product monitoring, a method and system for predicting the density of subsequent oil product monitoring at the mixed interface of a refined oil pipeline is proposed. The refined oil pipeline selected from historical data The hydraulic and thermal information of the upstream and downstream stations, the monitoring density of the upstream and downstream stations for oil products in the preceding row, and the monitoring density difference between the upstream and downstream stations of the preceding row of oil products are used as key input feature variables. The difference in monitoring density between the upstream and downstream stations of the trailing oil products is selected as For output variables, the Gaussian mixture regression algorithm is used to perform multi-modal identification of key input feature variables and output variables. Based on the modal identification results, the maximum expectation algorithm is used for training to obtain the monitoring density difference between upstream and downstream oil products corresponding to the modal. value correction model; and combined with the obtained monitoring density value of the oil product upstream station behind the product oil pipeline, the prediction results are obtained. The density of oil monitoring can be accurately predicted after the mixed oil interface.

Description

成品油管道混油界面后行油品监测密度预测方法及系统Method and system for density prediction of oil product monitoring at the mixed oil interface of refined oil pipelines

技术领域Technical field

本发明属于数据处理技术领域,尤其涉及成品油管道混油界面后行油品监测密度预测方法及系统。The invention belongs to the field of data processing technology, and in particular relates to a density prediction method and system for oil product monitoring after the mixed oil interface of a product oil pipeline.

背景技术Background technique

本部分的陈述仅仅是提供了与本发明相关的背景技术信息,不必然构成在先技术。The statements in this section merely provide background technical information related to the present invention and do not necessarily constitute prior art.

出于经济性考虑,多种成品油通常按照一定批次顺序,在相同成品油管道内连续输送。管输过程中相邻批次间将产生混油界面,站场安装的密度计虽难以获取混油界面前后行油品的真实密度值,但其实时反馈的混油界面监测密度随时间变化情况是操作人员处理混油界面的核心数据。具体而言,现场一般将混油界面密度监测值随时间变化曲线转为混油浓度分布随时间变化曲线,当混油浓度值降至一定阈值时采用批次切割方法处理混油界面。然而密度计仅能感知当前到站混油界面密度监测值,无法预知尚未到站的纯净后行油品密度监测值。若能精准把控后行油品密度监测值,可为准确指导现场油品批次切割工作提供关键数据支撑。然而管输过程水力热力条件复杂多变,基于现有油品密度计算公式无法准确获取后行油品监测密度;硬件设备普遍面临零点漂移现象,导致相同油品在上下游站场的监测密度值存在明显不同,困扰油品批次切割工作。此外,管输工况频繁变化导致数据存在多模态特性,若直接采用数据驱动建模方法预测后行油品密度,模型容易陷入过拟合误区,无法提供后行油品监测密度的可靠预测结果。For economic reasons, multiple refined oil products are usually transported continuously in the same refined oil pipeline in a certain batch sequence. During the pipeline transportation process, there will be an oil-mixed interface between adjacent batches. Although it is difficult for the density meter installed at the station to obtain the true density value of the oil products before and after the oil-mixed interface, its real-time feedback of the oil-mixed interface monitors the density changes over time. It is the core data for operators to handle the mixed oil interface. Specifically, on-site, the monitoring value of the mixed oil interface density is generally converted into a time change curve of the mixed oil concentration distribution. When the mixed oil concentration value drops to a certain threshold, the batch cutting method is used to process the mixed oil interface. However, the density meter can only sense the current density monitoring value of the mixed oil interface at the station, and cannot predict the density monitoring value of the pure oil that has not yet arrived at the station. If the oil density monitoring value can be accurately controlled, key data support can be provided to accurately guide on-site oil batch cutting work. However, the hydraulic and thermal conditions in the pipeline transportation process are complex and changeable. Based on the existing oil density calculation formula, it is impossible to accurately obtain the monitoring density of downstream oil products. Hardware equipment generally faces zero-point drift, resulting in the monitoring density values of the same oil products at upstream and downstream stations. There are obvious differences that trouble the oil batch cutting work. In addition, frequent changes in pipeline transportation conditions lead to multi-modal characteristics in the data. If the data-driven modeling method is directly used to predict the density of downstream oil products, the model will easily fall into the error of overfitting and cannot provide reliable predictions of the monitoring density of downstream oil products. result.

发明内容Contents of the invention

为克服上述现有技术的不足,本发明提供了成品油管道混油界面后行油品监测密度预测方法及系统,通过高斯混合回归算法对成品油管道数据进行多模态的识别,利用最大期望算法结合多模态识别结果进行训练得到的用于预测的后行油品上下游站场监测密度差值修正模型,从而可以准确预测混油界面后行油品监测密度。In order to overcome the shortcomings of the above-mentioned prior art, the present invention provides a method and system for predicting the oil product density after the mixed interface of a product oil pipeline. The Gaussian mixture regression algorithm is used to perform multi-modal identification of the product oil pipeline data and utilize the maximum expectation. The algorithm is trained based on multi-modal recognition results to obtain a correction model for predicting the monitoring density difference of upstream and downstream oil products at upstream and downstream stations, so that the monitoring density of downstream oil products at the mixed oil interface can be accurately predicted.

为实现上述目的,本发明的第一个方面提供成品油管道混油界面后行油品监测密度预测方法,包括:In order to achieve the above object, the first aspect of the present invention provides a density prediction method for oil product monitoring after the oil mixed interface of a product oil pipeline, including:

获取成品油管道混油界面的历史数据;Obtain historical data of the mixed oil interface of product oil pipelines;

选取历史数据中成品油管道上下游站场水力热力信息,前后行油品上游站场监测密度以及前行油品上下游站场监测密度差值作为关键输入特征变量,选取后行油品上下游站场监测密度差值作为输出变量,利用高斯混合回归算法对关键输入特征变量和输出变量进行多模态识别,得到各模态下的模态识别结果;Select the hydraulic and thermal information of the upstream and downstream stations of the refined oil pipeline in the historical data, the monitoring density of the upstream and downstream stations of the preceding oil products, and the monitoring density difference between the upstream and downstream stations of the preceding oil products as the key input feature variables, and select the upstream and downstream stations of the subsequent oil products. The station monitoring density difference is used as the output variable, and the Gaussian mixture regression algorithm is used to perform multi-modal identification of key input feature variables and output variables, and the modal identification results in each mode are obtained;

根据模态识别结果,利用最大期望算法进行训练得到模态对应的后行油品上下游站场监测密度差值修正模型;According to the mode identification results, the maximum expectation algorithm is used for training to obtain the density difference correction model for upstream and downstream oil product monitoring at the upstream and downstream stations corresponding to the mode;

利用后行油品上下游站场监测密度差值修正模型对待预测成品油管道混油界面进行预测,结合已获取的成品油管道后行油品上游站场监测密度值,得到待预测成品油管道下游站场的后行油品监测密度值。The correction model of the density difference between upstream and downstream oil product monitoring stations is used to predict the mixed oil interface of the product oil pipeline to be predicted. Combined with the obtained monitoring density value of the oil product upstream station for the product oil pipeline, the product oil pipeline to be predicted is obtained. Monitoring density value of downstream oil products at downstream stations.

本发明的第二个方面提供成品油管道混油界面后行油品监测密度预测系统,包括:The second aspect of the invention provides an oil product monitoring and density prediction system behind the mixed interface of a refined oil pipeline, including:

获取模块:获取成品油管道混油界面的历史数据;Acquisition module: obtain historical data of the mixed oil interface of the refined oil pipeline;

高斯模块:选取历史数据中成品油管道上下游站场水力热力信息,前后行油品上游站场监测密度以及前行油品上下游站场监测密度差值作为关键输入特征变量,选取后行油品上下游站场监测密度差值作为输出变量,利用高斯混合回归算法对关键输入特征变量和输出变量进行多模态识别,得到各模态下的模态识别结果;Gaussian module: Select the hydraulic and thermal information of the upstream and downstream stations of the refined oil pipeline in the historical data, the monitoring density of the upstream stations of the preceding oil products, and the monitoring density difference between the upstream and downstream stations of the preceding oil products as the key input feature variables, and select the following oil products The difference in monitoring density between the upstream and downstream stations of the product is used as the output variable, and the Gaussian mixture regression algorithm is used to perform multi-modal identification of key input feature variables and output variables, and the modal identification results in each mode are obtained;

训练模块:根据模态识别结果,利用最大期望算法进行训练得到模态对应的后行油品上下游站场监测密度差值修正模型;Training module: Based on the modal identification results, use the maximum expectation algorithm for training to obtain the density difference correction model for upstream and downstream oil product monitoring at the upstream and downstream stations corresponding to the modal;

预测模块:利用后行油品上下游站场监测密度差值修正模型对待预测成品油管道混油界面进行预测,结合已获取的成品油管道后行油品上游站场监测密度值,得到待预测成品油管道下游站场的后行油品监测密度值。Prediction module: Use the correction model of the monitoring density difference between the upstream and downstream stations of the downstream oil products to predict the mixed oil interface of the refined oil pipeline to be predicted, and combine it with the obtained monitoring density value of the upstream oil products in the refined oil pipeline to obtain the prediction Monitoring density value of downstream oil products at downstream stations of refined oil pipelines.

以上一个或多个技术方案存在以下有益效果:One or more of the above technical solutions have the following beneficial effects:

在本发明中,获取成品油管道的历史数据,从历史数据中选取上下游站场水力热力信息,前后行油品上游站场监测密度以及前行油品上下游站场监测密度差值作为关键输入特征变量,以后行油品上下游站场监测密度差值作为输出变量,利用高斯回归算法进行模态的识别,对模态识别结果基于最大期望算法进行训练,得到后行油品上下游站场监测密度差值修正模型,利用所得到的后行油品上下游站场监测密度差值修正模型对待预测的后行油品上下游站场监测密度差值进行预测;通过高斯回归算法进行多模态特征的识别,提高混油界面后行油品监测密度预测的准确性,对指导现场开展油品批次管理具有重要意义。In the present invention, the historical data of the refined oil pipeline is obtained, the hydraulic and thermal information of the upstream and downstream stations is selected from the historical data, and the monitoring density of the upstream and downstream stations of the preceding oil products and the difference in monitoring density of the upstream and downstream stations of the preceding oil products are used as the key Input the characteristic variables, and use the monitoring density difference between the upstream and downstream stations of the downstream oil products as the output variable. The Gaussian regression algorithm is used to identify the mode. The modal recognition results are trained based on the maximum expectation algorithm, and the upstream and downstream stations of the downstream oil products are obtained. The on-site monitoring density difference correction model uses the obtained upstream and downstream oil product monitoring density difference correction model to predict the to-be-predicted on-offline oil product upstream and downstream station monitoring density difference; a Gaussian regression algorithm is used to perform multiple The identification of modal characteristics improves the accuracy of oil monitoring density prediction after the mixed oil interface, and is of great significance in guiding on-site oil batch management.

本发明附加方面的优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。Additional advantages of the invention will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention.

附图说明Description of the drawings

构成本发明的一部分的说明书附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。The description and drawings that constitute a part of the present invention are used to provide a further understanding of the present invention. The illustrative embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute an improper limitation of the present invention.

图1为本发明实施例一中混油界面测量运移示意图;Figure 1 is a schematic diagram of the migration measured at the oil-mixed interface in Embodiment 1 of the present invention;

图2为本发明实施例一中常见的混油密度分布曲线;Figure 2 is a common mixed oil density distribution curve in Embodiment 1 of the present invention;

图3为本发明实施例一中由后行油品监测密度预测偏差带来的混油浓度分布曲线畸变展示图;Figure 3 is a diagram showing the distortion of the mixed oil concentration distribution curve caused by the density prediction deviation of subsequent oil monitoring in Embodiment 1 of the present invention;

图4为本发明实施例一中后行油品密度监测值预测模型流程图;Figure 4 is a flow chart of the subsequent oil product density monitoring value prediction model in Embodiment 1 of the present invention;

图5为本发明实施例一中后行油品上下游站场密度监测差值修正模型示意图。Figure 5 is a schematic diagram of the correction model for the density monitoring difference between upstream and downstream oil product stations in Embodiment 1 of the present invention.

具体实施方式Detailed ways

应该指出,以下详细说明都是示例性的,旨在对本发明提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本发明所属技术领域的普通技术人员通常理解的相同含义。It should be noted that the following detailed description is exemplary and is intended to provide further explanation of the present invention. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs.

需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本发明的示例性实施方式。It should be noted that the terms used herein are for the purpose of describing specific embodiments only, and are not intended to limit the exemplary embodiments according to the present invention.

在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。The embodiments and features of the embodiments of the present invention may be combined with each other without conflict.

实施例一Embodiment 1

本实施例公开了成品油管道混油界面后行油品监测密度预测方法,包括:This embodiment discloses a density prediction method for oil product monitoring behind the mixed oil interface of a product oil pipeline, including:

获取成品油管道混油界面的历史数据;Obtain historical data of the mixed oil interface of product oil pipelines;

选取历史数据中成品油管道上下游站场水力热力信息,前后行油品上游站场监测密度以及前行油品上下游站场监测密度差值作为关键输入特征变量,选取后行油品上下游站场监测密度差值作为输出变量,利用高斯混合回归算法对关键输入特征变量和输出变量进行多模态识别,得到各模态下的模态识别结果;Select the hydraulic and thermal information of the upstream and downstream stations of the refined oil pipeline in the historical data, the monitoring density of the upstream and downstream stations of the preceding oil products, and the monitoring density difference between the upstream and downstream stations of the preceding oil products as the key input feature variables, and select the upstream and downstream stations of the subsequent oil products. The station monitoring density difference is used as the output variable, and the Gaussian mixture regression algorithm is used to perform multi-modal identification of key input feature variables and output variables, and the modal identification results in each mode are obtained;

根据模态识别结果,利用最大期望算法进行训练得到模态对应的后行油品上下游站场监测密度差值修正模型;According to the mode identification results, the maximum expectation algorithm is used for training to obtain the density difference correction model for upstream and downstream oil product monitoring at the upstream and downstream stations corresponding to the mode;

利用后行油品上下游站场监测密度差值修正模型对待预测成品油管道混油界面进行预测,结合已获取的成品油管道后行油品上游站场监测密度值,得到待预测成品油管道下游站场的后行油品监测密度值。The correction model of the density difference between upstream and downstream oil product monitoring stations is used to predict the mixed oil interface of the product oil pipeline to be predicted. Combined with the obtained monitoring density value of the oil product upstream station for the product oil pipeline, the product oil pipeline to be predicted is obtained. Monitoring density value of downstream oil products at downstream stations.

本实施例提供的成品油管道混油界面后行油品监测密度预测方法主要包括三部分:建立成品油管道上下游站场混油数据库、建立后行油品上下游站场监测密度差值修正模型,混油界面后行油品监测密度预测。The method for predicting the density of oil products monitored at the mixed oil interface of the refined oil pipeline provided by this embodiment mainly includes three parts: establishing a mixed oil database at the upstream and downstream stations of the refined oil pipeline, and establishing the density difference correction for monitoring the density of the upstream and downstream oil products at the upstream and downstream stations. Model, oil product monitoring density prediction after mixed oil interface.

(1)建立成品油管道上下游站场混油数据库(1) Establish a mixed oil database at upstream and downstream stations of refined oil pipelines

建立成品油管道上下游站场混油数据库是基于上下游站场测量设备,获取管输水力热力信息及前后行油品监测密度并建立数据库,为后续建立后行油品监测密度差值修正模型提供数据支撑。其中,包括确定一定时段内上游站场出口温度压力、下游站场进口温度压力与上游站场前后行油品监测密度、下游站场前后行油品监测密度。The establishment of a mixed oil database at the upstream and downstream stations of the refined oil pipeline is based on the measurement equipment at the upstream and downstream stations. It obtains the pipeline water and thermal information and the front and rear oil product monitoring densities and establishes a database for subsequent establishment of the rear oil product monitoring density difference correction. The model provides data support. This includes determining the outlet temperature and pressure of the upstream station, the inlet temperature and pressure of the downstream station, the density of oil product monitoring before and after the upstream station, and the density of oil product monitoring before and after the downstream station within a certain period of time.

确定上游站场出口温度压力具体为:收集一定时段内成品油管道上游站场出口温度压力,其中温度监测样本量为,压力监测样本量为/>,取平均处理后可获得上游站场出口时均温度/>,℃;时均压力/>,单位:MPa:Determining the temperature and pressure at the outlet of the upstream station is specifically: collecting the temperature and pressure at the outlet of the upstream station of the refined oil pipeline within a certain period of time, where the temperature monitoring sample size is , the pressure monitoring sample size is/> , after average processing, the hourly average temperature of the upstream station outlet can be obtained/> , ℃; time average pressure/> , unit: MPa:

(1) (1)

(2) (2)

其中,为上游站场收集到的第/>个温度数据点,℃;/>为上游站场收集到的第/>个压力数据点,MPa。in, Collected for the upstream station/> temperature data points, ℃;/> Collected for the upstream station/> pressure data points, MPa.

确定下游站场进口温度压力具体为:收集一定时段内成品油管道下游站场进口温度压力,其中,温度监测样本量为,压力监测样本量为/>,取平均处理后可获得下游站场进口时均温度/>,单位:℃;时均压力/>,单位MPa:Determining the temperature and pressure at the downstream station inlet is specifically: collecting the temperature and pressure at the downstream station inlet of the refined oil pipeline within a certain period of time, where the temperature monitoring sample size is , the pressure monitoring sample size is/> , after average processing, the hourly average temperature of the downstream station inlet can be obtained/> , unit: ℃; hourly average pressure/> , unit MPa:

(3) (3)

(4) (4)

其中,为下游站场收集到的第/>个温度数据点,℃,/>为下游站场收集到的第/>个压力数据点,MPa。in, Collected for downstream stations/> temperature data points, ℃,/> Collected for downstream stations/> pressure data points, MPa.

确定上下游站场获取的前行油品监测密度差值与后行油品监测密度差值具体为:基于上下游站场密度计,获取上游站场前行油品密度监测值/>,单位:kg/m3,上游站场后行油品密度监测值/>,单位:kg/m3;下游站场前行油品密度监测值/>,单位:kg/m3;下游站场后行油品密度监测值/>,单位:kg/m3;据此计算上下游前行油品监测密度差值/>与上下游后行油品监测密度差值/>Determine the density difference of forward oil monitoring obtained from upstream and downstream stations Density difference from subsequent oil monitoring Specifically: based on the upstream and downstream station density meters, obtain the forward oil density monitoring value of the upstream station/> , unit: kg/m3, oil density monitoring value behind the upstream station/> , unit: kg/m3; oil density monitoring value in front of downstream station/> , unit: kg/m3; oil density monitoring value behind the downstream station/> , unit: kg/m3; based on this, calculate the density difference of upstream and downstream oil product monitoring/> Difference in density from upstream and downstream oil product monitoring/> :

(5) (5)

(6) (6)

(2)建立后行油品上下游站场监测密度差值修正模型(2) Establish a density difference correction model for upstream and downstream oil product monitoring stations

建立后行油品上下游站场监测密度差值修正模型是基于已有管输混油界面历史数据,结合高斯混合回归算法(Gaussian Mixture Regression Model, GMR),形成后行油品监测密度差值修正模型,具体包括:The establishment of the density difference correction model for upstream and downstream oil product monitoring at upstream and downstream stations is based on the historical data of the existing pipeline mixed oil interface, combined with the Gaussian Mixture Regression Model (GMR), to form the density difference for downstream oil product monitoring. Modify the model, specifically including:

选取模态识别关键特征变量:令和/>分别表征混油数据的输入变量矩阵与输出变量矩阵,/>为训练集样本量,/>为矩阵转置操作,表征第/>个训练集样本输入变量向量。Select key feature variables for modal recognition: let and/> Respectively represent the input variable matrix and output variable matrix of mixed oil data,/> is the sample size of the training set,/> For the matrix transpose operation, Characterization/> training set sample input variable vector.

为准确探索管输混油数据潜藏的多模态信息,选取上下游站场获取的温度、压力/>与/>以及上游站场获取的前行油品监测密度/>、后行油品监测密度/>及上下游站场获取的前行油品监测密度差值/>作为GMR算法识别数据多模态的关键输入特征变量。In order to accurately explore the multi-modal information hidden in the pipeline mixed oil data, the temperatures obtained from the upstream and downstream stations were selected. and , pressure/> with/> And the forward oil monitoring density obtained from the upstream station/> , rear oil monitoring density/> And the forward oil monitoring density difference obtained from upstream and downstream stations/> As the key input feature variable for the GMR algorithm to identify multi-modal data.

GMR中假定输入变量由个高斯分布组成。输入变量在第/>个分布,即第/>个模态下的边缘概率密度函数表达式/>如下:It is assumed in GMR that the input variable is given by consists of a Gaussian distribution. Input variables are at/> distribution, that is, the /> The expression of the marginal probability density function under each mode/> as follows:

(7) (7)

其中,为高斯分布,/>,/>分别指代第/>个模态下高斯分布中的均值向量与协方差矩阵。in, is a Gaussian distribution,/> ,/> Respectively refer to No./> The mean vector and covariance matrix in the Gaussian distribution under each mode.

确定回归关键输入变量:将上下游站场获取的前行油品上下游站场监测密度差值作为GMR算法回归过程的关键输入变量/>,因后行油品上下游站场监测密度差值为待预测变量,二者关系可表示为:Determine the key input variables of the regression: the difference in monitoring density of the upstream and downstream oil products obtained from the upstream and downstream stations As a key input variable in the regression process of GMR algorithm/> , due to the difference in density monitored by upstream and downstream oil product stations is the variable to be predicted, and the relationship between the two can be expressed as:

+/>(8) +/> (8)

其中,表示载第/>个模态下的高斯白噪声,其服从均值为0,方差为/>的高斯分布,T表示转置操作。在GMR中,/>个模态由对应/>个局部模型表征,/>为第/>个局部模型中的回归系数。in, Indicates the location/> Gaussian white noise in each mode has a mean value of 0 and a variance of/> Gaussian distribution, T represents the transpose operation. In GMR,/> Each modality is represented by the corresponding/> A local model representation,/> For the first/> regression coefficients in a local model.

训练后行油品监测密度差值修正模型:基于已有管输混油界面历史数据,结合最大期望算法(Expectation Maximization, EM)可估计模型参数。采用隐变量表征每个样本的模态分属情况,当/>时表明第/>个训练集样本分属于第/>个模态。计算/>属于第/>个模态的后验概率/>,并定义统计量/>Density difference correction model for oil product monitoring after training: Model parameters can be estimated based on the existing historical data of the mixed oil interface in pipelines and combined with the Expectation Maximization (EM) algorithm. Use latent variables Characterizes the modal belonging of each sample, when/> Time indicates the number/> The training set samples belong to the /> a modal. Calculate/> Belongs to/> Posterior probability of a mode/> , and define statistics/> :

(9) (9)

(10) (10)

其中,为给定第/>个样本的输入输出特征信息时,第/>个样本分属第/>个模态的条件概率;/>表征第/>个样本分属第/>个模态时,输出变量关于输入变量的条件概率;/>为第/>个样本分属第/>个模态时,输入变量的条件概率分布。in, For the given number/> When inputting and outputting feature information of a sample, the Samples belong to/> Conditional probability of a mode;/> Characterization/> Samples belong to/> mode, the conditional probability of the output variable with respect to the input variable;/> For the first/> Samples belong to/> mode, the conditional probability distribution of the input variable.

基于已获取的后验概率,可计算各模态高斯分布权重/>,均值/>及协方差/>,表达式如下:Based on the obtained posterior probability , can calculate the Gaussian distribution weight of each mode/> ,mean/> and covariance/> , the expression is as follows:

(11) (11)

(12) (12)

(13) (13)

(14) (14)

(15) (15)

其中,;/>为对角矩阵/>in, ;/> is a diagonal matrix , /> .

成品油管道混油界面后行油品监测密度预测是对后行油品监测密度待预测的混油界面,基于上下游站场测量设备获取的水力热力信息及油品监测密度信息,预测后行油品监测密度差值,结合上游站场获取的后行油品监测密度值/>,预测后行油品下游站场监测密度值/>,具体包括:The density prediction of the downstream oil product monitoring at the mixed oil interface of the refined oil pipeline is to predict the mixed oil interface where the downstream oil product monitoring density is to be predicted. Based on the hydrothermal information and oil product monitoring density information obtained by the upstream and downstream station measurement equipment, the prediction of the downstream oil product monitoring density is Oil monitoring density difference , combined with the downstream oil monitoring density value obtained from the upstream station/> , predict the monitoring density value of downstream oil products/> , specifically including:

确定待预测样本的关键变量信息:以第个待预测样本为例,结合上下游测量设备,获取上下游站场温度/>与/>、压力/>与/>以及上游站场获取的前行油品监测密度/>、后行油品监测密度/>及获取的前行油品上下游站场监测密度差值/>,组成GMR算法中的模态识别变量/>,前行油品上下游站场监测密度差值/>组成GMR算法中的回归变量Determine the key variable information of the sample to be predicted: Take the following Take a sample to be predicted as an example, and combine the upstream and downstream measurement equipment to obtain the temperature of the upstream and downstream stations/> with/> , pressure/> with/> And the forward oil monitoring density obtained from the upstream station/> , rear oil monitoring density/> And the obtained density difference between the upstream and downstream stations of Qianxing Oil Products/> , which constitutes the modal identification variable in the GMR algorithm/> , monitor the density difference between the upstream and downstream stations of Qianxing Oil/> Regressors that make up the GMR algorithm .

预测后行油品监测密度差值:对于第/>个待预测样本/>,/>为已知量,为待预测量,可基于式(16)计算其分属/>个模态的后验概率/>Oil product monitoring density difference after prediction : For the first/> samples to be predicted/> ,/> is a known quantity, is the quantity to be predicted, its distribution can be calculated based on equation (16)/> Posterior probability of a mode/> :

(16) (16)

其中,为给定第/>个样本的输入输出特征信息时,第/>个样本分属第/>个模态的条件概率,/>为第/>个样本分属第/>个模态的概率,/>为第/>个样本分属第/>个模态时,输入变量的条件概率分布。in, For the given number/> When inputting and outputting feature information of a sample, the Samples belong to/> The conditional probability of a mode,/> For the first/> Samples belong to/> The probability of a mode,/> For the first/> Samples belong to/> mode, the conditional probability distribution of the input variable.

相应地,可得到待预测量的条件概率表达式/>Correspondingly, the quantity to be predicted can be obtained Conditional probability expression/> :

(17) (17)

其中,给定输入变量时,待预测量/>的条件概率表达式;为输出变量/>服从均值为/>,方差为/>的高斯分布。in, Given input variables, the quantity to be predicted/> conditional probability expression; For output variables/> Obey the mean as/> , the variance is/> Gaussian distribution.

考虑将的期望作为预测结果/>,因此后行油品监测密度差值/>预测值/>表达式为:Consider expectations as predicted results/> , so the subsequent oil monitoring density difference/> Predicted value/> The expression is:

(18) (18)

其中,为/>在给定/>下的条件期望。in, for/> In given/> Conditional expectations below.

预测后行油品监测密度值:结合上游站场获取的后行油品监测密度值/>,预测下游站场后行油品监测密度值/>Predict the oil product monitoring density value : Combined with the downstream oil monitoring density value obtained from the upstream station/> , predict the oil product monitoring density value behind the downstream station/> :

(19) (19)

图1是混油界面测量运移示意图,多种成品油产品按照一定批次,基于成品油管道从上游厂商泵送至下游用户。混油界面过站过程中,准确估计后行油品监测密度,确定混油浓度分布信息,是实现油品批次精准切割,提高成品油管网运行经济效益的关键前提。Figure 1 is a schematic diagram of measured migration at the mixed oil interface. Various refined oil products are pumped from upstream manufacturers to downstream users based on refined oil pipelines in certain batches. During the passing process of the mixed oil interface, accurately estimating the subsequent oil monitoring density and determining the mixed oil concentration distribution information is a key prerequisite for achieving accurate oil batch cutting and improving the economic benefits of the refined oil pipeline network operation.

图2展示了常见的混油密度分布曲线,当混油界面经过上游站场时,站场安装的密度传感器最先监测到批次顺序在前的前行油品监测密度时序信息,之后逐渐感应到后行油品密度监测信息。混油浓度分布信息需基于前后行油品监测密度信息进行计算,然而混油界面在相邻两站间管道运移时,下游站场仅能感应到前行油品监测密度。目前尚缺乏直接监测后行油品密度方法,因此有必要结合现场数据,建立后行油品监测密度的高精度预测方法,为油品批次切割工作提供数据支撑。Figure 2 shows a common mixed oil density distribution curve. When the mixed oil interface passes through the upstream station, the density sensor installed at the station first detects the density timing information of the preceding oil product in the batch sequence, and then gradually senses it. After receiving the oil density monitoring information. The mixed oil concentration distribution information needs to be calculated based on the preceding and following oil product monitoring density information. However, when the mixed oil interface migrates in the pipeline between two adjacent stations, the downstream station can only sense the preceding oil product monitoring density. At present, there is still a lack of methods for directly monitoring the density of oil products in the pipeline. Therefore, it is necessary to combine on-site data to establish a high-precision prediction method for monitoring the density of oil products in the pipeline to provide data support for oil batch cutting work.

图3展示了由后行油品监测密度预测偏差带来的混油浓度分布信息误差,其中δ表征了后行油品监测密度预测误差,δ=0代表实际混油浓度分布曲线。随着后行油品监测密度预测误差逐渐增大,相应混油浓度分布曲线与真实曲线偏离程度越来越明显:混油浓度分布理应介于[0,1]区间,然而当后行油品监测密度预测结果出现误差时,混油浓度值可能出现负数,计算结果违背物理规律。此外,若将混油浓度值为1%定义为油品批次切割阈值点,当δ从0增长至1,3,5时,表征混油界面的油品批次切割阈值点对应时间偏差分别为216s,450s与624s。上述结果再次说明后行油品监测密度的精准预测对准确开展油品批次切割工作的重要意义。Figure 3 shows the mixed oil concentration distribution information error caused by the subsequent oil product monitoring density prediction error, where δ represents the subsequent oil product monitoring density prediction error, and δ=0 represents the actual mixed oil concentration distribution curve. As the monitoring density prediction error of subsequent oil products gradually increases, the corresponding mixed oil concentration distribution curve deviates more and more from the true curve: the mixed oil concentration distribution should be in the [0,1] interval, but when the subsequent oil products When there is an error in the monitoring density prediction result, the mixed oil concentration value may appear negative, and the calculation result violates the laws of physics. In addition, if the mixed oil concentration value of 1% is defined as the oil batch cutting threshold point, when δ increases from 0 to 1, 3, and 5, the corresponding time deviations of the oil batch cutting threshold points that represent the mixed oil interface are respectively for 216s, 450s and 624s. The above results once again illustrate the importance of accurate prediction of subsequent oil monitoring density for accurate oil batch cutting work.

图4展示了后行油品密度监测值预测模型流程,首先基于上下游站场安装的测量设备,分别获取温度压力及油品密度信息。需要注意的是,为精准建模,上游站场密度计监测到前后行油品密度后即可知油品批次顺序,如汽油推送柴油,简称汽柴界面;柴油推送汽油,简称柴汽界面,以此类推,据此建立相对应的后行油品密度监测差值修正模型。Figure 4 shows the downstream oil density monitoring value prediction model process. First, based on the measurement equipment installed at upstream and downstream stations, temperature, pressure and oil density information are obtained respectively. It should be noted that for accurate modeling, the upstream station density meter can know the order of oil batches after monitoring the density of oil products in the front and rear rows. For example, gasoline pushes diesel, which is referred to as the gasoline-diesel interface; diesel pushes gasoline, which is referred to as the diesel-gasoline interface. By analogy, a corresponding subsequent oil product density monitoring difference correction model is established accordingly.

图5展示了后行油品密度监测差值修正模型示意图;困扰后行油品监测密度准确预测的原因包括:柴油和汽油的密度随温度和压力的变化而变化,复杂的操作条件导致了密度与水力热力,即压力温度数据之间的高度非线性关系;硬件设备普遍面临零点漂移现象,这意味着即使在相似的温度压力条件下,上下游密度计针对相同油品可能会得到不同的测量结果。因此若建模时直接预测后行油品监测密度,水力热力条件与不同硬件设备固有的测量误差可能会导致不可靠的预测结果。Figure 5 shows a schematic diagram of the difference correction model for downstream oil product density monitoring; the reasons that trouble accurate prediction of downstream oil product density monitoring include: the density of diesel and gasoline changes with changes in temperature and pressure, and complex operating conditions lead to density changes. It has a highly nonlinear relationship with hydrothermal power, that is, pressure and temperature data; hardware equipment generally faces zero point drift, which means that even under similar temperature and pressure conditions, upstream and downstream density meters may obtain different measurements for the same oil product. result. Therefore, if the subsequent oil monitoring density is directly predicted during modeling, the measurement errors inherent in hydrothermal conditions and different hardware equipment may lead to unreliable prediction results.

然而,前后行油品在管输过程经历了相似的水力热力过程,而且两者在上下游均采用了相同设备测量其密度信息。故此前行油品上下游站场监测密度差值与后行油品上下游站场监测密度差值/>间必然高度关联。采用数据驱动模型捕捉间的函数关系,即对前行油品上下游站场监测密度差值/>与后行油品上下游站场监测密度差值/>间的依赖关系进行建模,修正后行油品上下游站场监测密度差值/>,是一种更合理的建模方法。However, the oil products from the front and back have experienced similar hydrothermal processes during pipeline transportation, and the same equipment is used to measure their density information both upstream and downstream. Therefore, it is currently necessary to monitor the density difference at upstream and downstream stations of oil products. Difference in density monitoring from upstream and downstream oil product stations/> There must be a high degree of correlation between them. Use a data-driven model to capture the functional relationship between them, that is, monitor the density difference between the upstream and downstream stations of Qianxing Oil Products/> Difference in density monitoring from upstream and downstream oil product stations/> Model the dependence between the oil products and monitor the density difference at upstream and downstream oil product stations after correction/> , is a more reasonable modeling method.

在此过程中,通过结合GMR算法引入模态识别功能,可使模型能在不同模态条件下对后行油品上下游站场监测密度差值与前行油品上下游站场监测密度差值/>间的函数关系进行精准建模,保证模型预测精度。In this process, the modal identification function is introduced by combining the GMR algorithm, so that the model can monitor the density difference of upstream and downstream oil stations under different modal conditions. Difference in density monitoring from upstream and downstream stations of Qianxing Oil Products/> Accurately model the functional relationship between them to ensure the accuracy of model prediction.

公式(20)展示了油品密度随温度变化的理论公式,其中表征了20摄氏度条件下油品的标准密度,单位为kg/m3。为展示本实施例提出的建模方法的有效性与优越性,采用理论公式与两种常见的机器学习算法,即梯度提升决策树(Gradient Boosted DecisionTree)与人工神经网络算法(Artificial Neural network,ANN)以及直接采用GMR算法直接预测后行油品监测密度得到的预测结果。此外,为尽可能全面地评估模型性能,将121个样本分为训练集与预测集,其中测试集比例从0.3逐步升至0.6。Equation (20) shows the theoretical formula for the change of oil density with temperature, where Characterizes the standard density of oil at 20 degrees Celsius, in kg/m3. In order to demonstrate the effectiveness and superiority of the modeling method proposed in this embodiment, theoretical formulas and two common machine learning algorithms are used, namely gradient boosted decision tree (Gradient Boosted DecisionTree) and artificial neural network algorithm (Artificial Neural network, ANN). ) and the prediction results obtained by directly using the GMR algorithm to directly predict the subsequent oil monitoring density. In addition, in order to evaluate the model performance as comprehensively as possible, 121 samples were divided into training sets and prediction sets, with the test set proportion gradually increasing from 0.3 to 0.6.

表1、2分别展示了汽柴界面与柴汽界面条件下基于均方根误差RMSE表征的各模型后行油品监测密度预测精度,误差指标计算式如下,其中与/>分别表征第/>个样本的真实值与拟合值,/>为测试样本量。RMSE越小,整体误差越低。Tables 1 and 2 respectively show the subsequent oil monitoring density prediction accuracy of each model based on the root mean square error RMSE under the conditions of gasoline-diesel interface and diesel-gasoline interface. The error index calculation formula is as follows, where with/> Respectively represent the first/> The true values and fitted values of the samples,/> is the test sample size. The smaller the RMSE, the lower the overall error.

(20) (20)

(21) (twenty one)

表2 汽柴界面下各模型预测结果:Table 2 Prediction results of each model under the gasoline and diesel interface:

表3柴汽界面下各模型预测结果:Table 3 Prediction results of each model under the diesel-gasoline interface:

由表2、3可得,本实施例提出的后行油品监测密度预测方法具备明显优势,证明本发明提出的方法对提高成品油管道混油界面定位具有重要意义。It can be seen from Tables 2 and 3 that the downstream oil monitoring density prediction method proposed in this embodiment has obvious advantages, proving that the method proposed by the present invention is of great significance in improving the positioning of the mixed oil interface in product oil pipelines.

实施例二Embodiment 2

本实施例的目的是提供成品油管道混油界面后行油品监测密度预测系统,包括:The purpose of this embodiment is to provide an oil product monitoring and density prediction system behind the mixed oil interface of a product oil pipeline, including:

成品油管道混油界面后行油品监测密度预测系统,其特征在于,包括:The oil product monitoring and density prediction system behind the mixed oil interface of the refined oil pipeline is characterized by:

获取模块:获取成品油管道混油界面的历史数据;Acquisition module: obtain historical data of the mixed oil interface of the refined oil pipeline;

高斯模块:选取历史数据中成品油管道上下游站场水力热力信息,前后行油品上游站场监测密度以及前行油品上下游站场监测密度差值作为关键输入特征变量,选取后行油品上下游站场监测密度差值作为输出变量,利用高斯混合回归算法对关键输入特征变量和输出变量进行多模态识别,得到各模态下的模态识别结果;Gaussian module: Select the hydraulic and thermal information of the upstream and downstream stations of the refined oil pipeline in the historical data, the monitoring density of the upstream stations of the preceding oil products, and the monitoring density difference between the upstream and downstream stations of the preceding oil products as the key input feature variables, and select the following oil products The difference in monitoring density between the upstream and downstream stations of the product is used as the output variable, and the Gaussian mixture regression algorithm is used to perform multi-modal identification of key input feature variables and output variables, and the modal identification results in each mode are obtained;

训练模块:根据模态识别结果,利用最大期望算法进行训练得到模态对应的后行油品上下游站场监测密度差值修正模型;Training module: Based on the modal identification results, use the maximum expectation algorithm for training to obtain the density difference correction model for upstream and downstream oil product monitoring at the upstream and downstream stations corresponding to the modal;

预测模块:利用后行油品上下游站场监测密度差值修正模型对待预测成品油管道混油界面进行预测,结合已获取的成品油管道后行油品上游站场监测密度值,得到待预测成品油管道下游站场的后行油品监测密度值。Prediction module: Use the correction model of the monitoring density difference between the upstream and downstream stations of the downstream oil products to predict the mixed oil interface of the refined oil pipeline to be predicted, and combine it with the obtained monitoring density value of the upstream oil products in the refined oil pipeline to obtain the prediction Monitoring density value of downstream oil products at downstream stations of refined oil pipelines.

在本实施例中,训练模块包括:In this embodiment, the training module includes:

第一计算模块:计算每个训练样本中关键输入特征变量的模态分属的后验概率;The first calculation module: calculates the posterior probability of the modal affiliation of the key input feature variables in each training sample;

第二计算模块:基于所获取的每个训练样本中关键特征变量的模态分属的后验概率,计算各模态的高斯分布的方差、均值和协方差;The second calculation module: Based on the obtained posterior probability of the modal affiliation of the key feature variables in each training sample, calculate the variance, mean and covariance of the Gaussian distribution of each mode;

第三计算模块:基于所计算得到的高斯分布的方差,建立关键输入特征变量与输出变量的关;The third calculation module: Based on the calculated variance of the Gaussian distribution, establish the relationship between the key input feature variables and the output variables;

第四计算模块:基于所计算得到的高斯分布的均值和协方差,建立关键输入特征变量的边缘概率密度函数。The fourth calculation module: Based on the calculated mean and covariance of the Gaussian distribution, establish the marginal probability density function of the key input feature variables.

在本实施例中,预测模块包括:In this embodiment, the prediction module includes:

确定模块:将待预测成品油管道上下游站场的温度、压力,上游站场获取的前后行油品监测密度作为高斯混合回归算法中的模态识别变量,将上下游站场获取的前行油品监测密度差值作为高斯混合回归算法中的回归变量;Determination module: Use the temperature and pressure of the upstream and downstream stations of the refined oil pipeline to be predicted, and the front and rear oil product monitoring densities obtained by the upstream stations as the modal identification variables in the Gaussian mixture regression algorithm, and use the forward and backward oil product monitoring densities obtained by the upstream and downstream stations as the modal identification variables in the Gaussian mixture regression algorithm. The oil product monitoring density difference is used as the regressor variable in the Gaussian mixture regression algorithm;

概率计算模块:计算模态识别变量模态分属的后验概率;Probability calculation module: calculates the posterior probability of the modal classification of the modal identification variable;

差值预测模块:基于后行油品监测密度差值修正模型、模态识别变量模态分属的后验概率以及回归变量,得到后行油品监测密度差值预测值。Difference prediction module: Based on the subsequent oil monitoring density difference correction model, the posterior probability of the modal classification of the modal identification variables, and the regression variables, the predicted value of the subsequent oil monitoring density difference is obtained.

在本实施例中,预测模块还包括:In this embodiment, the prediction module also includes:

相加模块:将所述后行油品监测密度差值预测值与所获取的成品油管道后行油品上游站场监测密度值相加,得到待预测成品油管道后行油品下游站场监测密度值。Addition module: Add the predicted density difference prediction value of the downstream oil product to the obtained monitoring density value of the downstream oil product upstream station of the refined oil pipeline to obtain the downstream oil product downstream station of the refined oil pipeline to be predicted Monitor the density value.

上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。Although the specific embodiments of the present invention have been described above in conjunction with the accompanying drawings, they do not limit the scope of the present invention. Those skilled in the art should understand that based on the technical solutions of the present invention, those skilled in the art do not need to perform creative work. Various modifications or variations that can be made are still within the protection scope of the present invention.

Claims (10)

1.成品油管道混油界面后行油品监测密度预测方法,其特征在于,包括:1. The density prediction method of oil product monitoring behind the mixed interface of refined oil pipelines is characterized by: 获取成品油管道混油界面的历史数据;Obtain historical data of the mixed oil interface of product oil pipelines; 选取历史数据中成品油管道上下游站场水力热力信息,前后行油品上游站场监测密度以及前行油品上下游站场监测密度差值作为关键输入特征变量,选取后行油品上下游站场监测密度差值作为输出变量,利用高斯混合回归算法对关键输入特征变量和输出变量进行多模态识别,得到各模态下的模态识别结果;Select the hydraulic and thermal information of the upstream and downstream stations of the refined oil pipeline in the historical data, the monitoring density of the upstream and downstream stations of the preceding oil products, and the monitoring density difference between the upstream and downstream stations of the preceding oil products as the key input feature variables, and select the upstream and downstream stations of the subsequent oil products. The station monitoring density difference is used as the output variable, and the Gaussian mixture regression algorithm is used to perform multi-modal identification of key input feature variables and output variables, and the modal identification results in each mode are obtained; 根据模态识别结果,利用最大期望算法进行训练得到模态对应的后行油品上下游站场监测密度差值修正模型;According to the mode identification results, the maximum expectation algorithm is used for training to obtain the density difference correction model for upstream and downstream oil product monitoring at the upstream and downstream stations corresponding to the mode; 利用后行油品上下游站场监测密度差值修正模型对待预测成品油管道混油界面进行预测,结合已获取的成品油管道后行油品上游站场监测密度值,得到待预测成品油管道下游站场的后行油品监测密度值。The correction model of the density difference between upstream and downstream oil product monitoring stations is used to predict the mixed oil interface of the product oil pipeline to be predicted. Combined with the obtained monitoring density value of the oil product upstream station for the product oil pipeline, the product oil pipeline to be predicted is obtained. Monitoring density value of downstream oil products at downstream stations. 2.如权利要求1所述的成品油管道混油界面后行油品监测密度预测方法,其特征在于,所述上下游站场水力热力信息包括上游站场出口温度、压力,下游站场进口温度、压力;基于成品油管道前行油品上游站场密度监测值与前行油品下游站场密度监测值确定前行油品上下游站场监测密度差值;基于成品油管道后行油品上游站场密度监测值与后行油品下游站场密度监测值确定后行油品上下游站场监测密度差值。2. The method for predicting the density of oil product monitoring at the mixed interface of a refined oil pipeline as claimed in claim 1, wherein the hydraulic and thermal information of the upstream and downstream stations includes the temperature and pressure of the outlet of the upstream station and the inlet of the downstream station. Temperature and pressure; determine the monitoring density difference between the upstream and downstream stations of the forward oil product based on the density monitoring value of the upstream oil product pipeline and the density monitoring value of the forward oil product downstream station; based on the rear oil product pipeline density monitoring value The density monitoring value of the upstream station of the product and the density monitoring value of the downstream station of the subsequent oil product determine the monitoring density difference between the upstream and downstream stations of the subsequent oil product. 3.如权利要求1所述的成品油管道混油界面后行油品监测密度预测方法,其特征在于,根据模态识别结果,利用最大期望算法进行训练得到模态对应的后行油品上下游站场监测密度差值修正模型,具体为:3. The method for monitoring and predicting the density of downstream oil products at the mixed interface of a refined oil pipeline as claimed in claim 1, characterized in that, based on the modal recognition results, the maximum expectation algorithm is used for training to obtain the downstream oil product density corresponding to the modality. Downstream station monitoring density difference correction model, specifically: 计算每个训练样本中关键输入特征变量的模态分属的后验概率;Calculate the posterior probability of the modal affiliation of the key input feature variables in each training sample; 基于所获取的每个训练样本中关键特征变量的模态分属的后验概率,计算各模态的高斯分布的方差、均值和协方差;Based on the obtained posterior probability of the modal affiliation of the key feature variables in each training sample, calculate the variance, mean and covariance of the Gaussian distribution of each modality; 基于所计算得到的高斯分布的方差,建立关键输入特征变量与输出变量的关系;Based on the calculated variance of the Gaussian distribution, establish the relationship between the key input feature variables and the output variables; 基于所计算得到的高斯分布的均值和协方差,建立关键输入特征变量的边缘概率密度函数。Based on the calculated mean and covariance of the Gaussian distribution, the marginal probability density function of the key input feature variables is established. 4.如权利要求1所述的成品油管道混油界面后行油品监测密度预测方法,其特征在于,利用后行油品监测密度差值修正模型对待预测成品油管道混油界面进行预测,具体为:4. The method for predicting the density of oil product monitoring at the mixed interface of a refined oil pipeline as claimed in claim 1, characterized in that the mixed oil interface of the refined oil pipeline to be predicted is predicted using the density difference correction model of the oil product monitored at the rear. Specifically: 将待预测成品油管道上下游站场的温度、压力,获取的前后行油品上游站场监测密度作为高斯混合回归算法中的模态识别变量,将获取的前行油品上下游站场监测密度差值作为高斯混合回归算法中的回归变量;The temperature and pressure of the upstream and downstream stations of the refined oil pipeline to be predicted, and the obtained monitoring density of the upstream and downstream oil products are used as the modal identification variables in the Gaussian mixture regression algorithm, and the obtained upstream and downstream station monitoring densities of the preceding oil products are The density difference is used as a regressor in the Gaussian mixture regression algorithm; 计算模态识别变量模态分属的后验概率;Calculate the posterior probability of the modal affiliation of the modal identification variable; 基于后行油品监测密度差值修正模型、模态识别变量模态分属的后验概率以及回归变量,得到后行油品下游站场监测密度差值的预测值。Based on the correction model of the density difference of downstream oil products monitoring, the posterior probability of the modal classification of the modal identification variables, and the regression variables, the predicted value of the monitoring density difference of downstream oil products downstream stations is obtained. 5.如权利要求4所述的成品油管道混油界面后行油品监测密度预测方法,其特征在于,将所述后行油品监测密度差值预测值与所获取的成品油管道后行油品上游站场监测密度值相加,得到待预测成品油管道后行油品下游站场监测密度值。5. The method for predicting the density of oil product monitoring at the mixed interface of a refined oil pipeline as claimed in claim 4, characterized in that the density difference prediction value of the oil product monitoring at the rear line is compared with the obtained density difference of the oil product at the mixed interface of the refined oil pipeline. The monitoring density values of the upstream oil product stations are added up to obtain the monitoring density values of the downstream oil products downstream of the refined oil pipeline to be predicted. 6.如权利要求4所述的成品油管道混油界面后行油品监测密度预测方法,其特征在于,根据模态识别变量模态分属的后验概率和回归变量,得到后行油品下游站场监测密度差值的期望,将所得到的期望作为后行油品下游站场监测密度差值的预测值。6. The density prediction method for monitoring and predicting the density of downstream oil products at the oil-mixing interface of a refined oil pipeline as claimed in claim 4, characterized in that, according to the posterior probability and regression variables of the modal assignment of the modal identification variables, the downstream oil products are obtained The expectation of the density difference monitored by the downstream station is used as the predicted value of the density difference monitored by the downstream oil product downstream station. 7.成品油管道混油界面后行油品监测密度预测系统,其特征在于,包括:7. The oil product monitoring and density prediction system behind the mixed oil interface of the refined oil pipeline is characterized by: 获取模块:获取成品油管道混油界面的历史数据;Acquisition module: obtain historical data of the mixed oil interface of the refined oil pipeline; 高斯模块:选取历史数据中成品油管道上下游站场水力热力信息,前后行油品上游站场监测密度以及前行油品上下游站场监测密度差值作为关键输入特征变量,选取后行油品上下游站场监测密度差值作为输出变量,利用高斯混合回归算法对关键输入特征变量和输出变量进行多模态识别,得到各模态下的模态识别结果;Gaussian module: Select the hydraulic and thermal information of the upstream and downstream stations of the refined oil pipeline in the historical data, the monitoring density of the upstream stations of the preceding oil products, and the monitoring density difference between the upstream and downstream stations of the preceding oil products as the key input feature variables, and select the following oil products The difference in monitoring density between the upstream and downstream stations of the product is used as the output variable, and the Gaussian mixture regression algorithm is used to perform multi-modal identification of key input feature variables and output variables, and the modal identification results in each mode are obtained; 训练模块:根据模态识别结果,利用最大期望算法进行训练得到模态对应的后行油品上下游站场监测密度差值修正模型;Training module: Based on the modal identification results, use the maximum expectation algorithm for training to obtain the density difference correction model for upstream and downstream oil product monitoring at the upstream and downstream stations corresponding to the modal; 预测模块:利用后行油品上下游站场监测密度差值修正模型对待预测成品油管道混油界面进行预测,结合已获取的成品油管道后行油品上游站场监测密度值,得到待预测成品油管道下游站场的后行油品监测密度值。Prediction module: Use the correction model of the monitoring density difference between the upstream and downstream stations of the downstream oil products to predict the mixed oil interface of the refined oil pipeline to be predicted, and combine it with the obtained monitoring density value of the upstream oil products in the refined oil pipeline to obtain the prediction Monitoring density value of downstream oil products at downstream stations of refined oil pipelines. 8.如权利要求7所述的成品油管道混油界面后行油品监测密度预测系统,其特征在于,所述训练模块,包括:8. The oil product monitoring and density prediction system behind the mixed interface of the refined oil pipeline according to claim 7, characterized in that the training module includes: 第一计算模块:计算每个训练样本中关键输入特征变量的模态分属的后验概率;The first calculation module: calculates the posterior probability of the modal affiliation of the key input feature variables in each training sample; 第二计算模块:基于所获取的每个训练样本中关键特征变量的模态分属的后验概率,计算各模态的高斯分布的方差、均值和协方差;The second calculation module: Based on the obtained posterior probability of the modal affiliation of the key feature variables in each training sample, calculate the variance, mean and covariance of the Gaussian distribution of each mode; 第三计算模块:基于所计算得到的高斯分布的方差,建立关键输入特征变量与输出变量的关;The third calculation module: Based on the calculated variance of the Gaussian distribution, establish the relationship between the key input feature variables and the output variables; 第四计算模块:基于所计算得到的高斯分布的均值和协方差,建立关键输入特征变量的边缘概率密度函数。The fourth calculation module: Based on the calculated mean and covariance of the Gaussian distribution, establish the marginal probability density function of the key input feature variables. 9.如权利要求7所述的成品油管道混油界面后行油品监测密度预测系统,其特征在于,所述预测模块,包括:9. The oil product monitoring and density prediction system behind the mixed interface of the refined oil pipeline according to claim 7, characterized in that the prediction module includes: 确定模块:将待预测成品油管道上下游站场的温度、压力,上游站场获取的前后行油品监测密度作为高斯混合回归算法中的模态识别变量,将上下游站场获取的前行油品监测密度差值作为高斯混合回归算法中的回归变量;Determination module: The temperature and pressure of the upstream and downstream stations of the refined oil pipeline to be predicted, and the front and rear oil product monitoring densities obtained by the upstream stations are used as modal identification variables in the Gaussian mixture regression algorithm, and the forward and backward oil product monitoring densities obtained by the upstream and downstream stations are used. The oil product monitoring density difference is used as the regressor variable in the Gaussian mixture regression algorithm; 概率计算模块:计算模态识别变量模态分属的后验概率;Probability calculation module: calculates the posterior probability of the modal classification of the modal identification variable; 差值预测模块:基于后行油品监测密度差值修正模型、模态识别变量模态分属的后验概率以及回归变量,得到后行油品监测密度差值的预测值。Difference prediction module: Based on the subsequent oil product monitoring density difference correction model, the posterior probability of modal identification variables and regression variables, the predicted value of the subsequent oil product monitoring density difference is obtained. 10.如权利要求7所述的成品油管道混油界面后行油品监测密度预测系统,其特征在于,所述预测模块,还包括:相加模块:将所述后行油品监测密度差值预测值与所获取的成品油管道后行油品上游站场监测密度值相加,得到待预测成品油管道后行油品下游场站监测密度值。10. The density prediction system for monitoring the downstream oil product at the mixed interface of the refined oil pipeline as claimed in claim 7, characterized in that the prediction module further includes: an addition module: adding the density difference of the downstream oil product monitoring The predicted value is added to the obtained upstream station monitoring density value of oil products flowing behind the refined oil pipeline to obtain the monitoring density value of the downstream station monitoring density of oil products flowing behind the refined oil pipeline to be predicted.
CN202311107413.5A 2023-08-31 2023-08-31 Method and system for predicting monitoring density of oil product of oil mixing interface of finished oil pipeline Active CN116822999B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311107413.5A CN116822999B (en) 2023-08-31 2023-08-31 Method and system for predicting monitoring density of oil product of oil mixing interface of finished oil pipeline

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311107413.5A CN116822999B (en) 2023-08-31 2023-08-31 Method and system for predicting monitoring density of oil product of oil mixing interface of finished oil pipeline

Publications (2)

Publication Number Publication Date
CN116822999A true CN116822999A (en) 2023-09-29
CN116822999B CN116822999B (en) 2023-12-05

Family

ID=88127906

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311107413.5A Active CN116822999B (en) 2023-08-31 2023-08-31 Method and system for predicting monitoring density of oil product of oil mixing interface of finished oil pipeline

Country Status (1)

Country Link
CN (1) CN116822999B (en)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2416153A1 (en) * 2010-08-02 2012-02-08 ENI S.p.A. Method for predicting the quality and yields of a crude oil
US20160140448A1 (en) * 2014-11-17 2016-05-19 Massachusetts Institute Of Technology Systems and methods for improving petroleum fuels production
RU2613385C1 (en) * 2016-03-10 2017-03-16 Федеральное государственное бюджетное образовательное учреждение высшего образования "Тихоокеанский государственный университет" Automated control system of oil quality
US20170145823A1 (en) * 2014-05-23 2017-05-25 Landmark Graphics Corporation Robust viscosity estimation methods and systems
WO2020261940A1 (en) * 2019-06-27 2020-12-30 積水ポリマテック株式会社 Battery module
CN112182885A (en) * 2020-09-29 2021-01-05 中国民用航空飞行学院 Fuel consumption deviation prediction method and system based on Gaussian mixture model
CN113196550A (en) * 2018-07-30 2021-07-30 凯尊创新有限公司 Housing for rechargeable battery
CN114861759A (en) * 2022-04-06 2022-08-05 中国石油大学(华东) A Distributed Training Method for Linear Dynamic System Models
CN116307303A (en) * 2023-05-24 2023-06-23 中国石油大学(华东) Mechanism-data dual-drive oil mixing length prediction method and system for finished oil pipeline

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2416153A1 (en) * 2010-08-02 2012-02-08 ENI S.p.A. Method for predicting the quality and yields of a crude oil
US20170145823A1 (en) * 2014-05-23 2017-05-25 Landmark Graphics Corporation Robust viscosity estimation methods and systems
US20160140448A1 (en) * 2014-11-17 2016-05-19 Massachusetts Institute Of Technology Systems and methods for improving petroleum fuels production
RU2613385C1 (en) * 2016-03-10 2017-03-16 Федеральное государственное бюджетное образовательное учреждение высшего образования "Тихоокеанский государственный университет" Automated control system of oil quality
CN113196550A (en) * 2018-07-30 2021-07-30 凯尊创新有限公司 Housing for rechargeable battery
WO2020261940A1 (en) * 2019-06-27 2020-12-30 積水ポリマテック株式会社 Battery module
CN113906611A (en) * 2019-06-27 2022-01-07 积水保力马科技株式会社 battery module
CN112182885A (en) * 2020-09-29 2021-01-05 中国民用航空飞行学院 Fuel consumption deviation prediction method and system based on Gaussian mixture model
CN114861759A (en) * 2022-04-06 2022-08-05 中国石油大学(华东) A Distributed Training Method for Linear Dynamic System Models
CN116307303A (en) * 2023-05-24 2023-06-23 中国石油大学(华东) Mechanism-data dual-drive oil mixing length prediction method and system for finished oil pipeline

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
ALESSIO SCANZIANI等人: "Automatic method for estimation of in situ effective contact angle from X-ray micro tomography images of two-phase flow in porous media", 《JOURNAL OF COLLOID AND INTERFACE SCIENCE》 *
袁子云等: "融合机制与高斯混合回归算法的成品油管道顺序输送混油长度预测模型", 《中国石油大学学报(自然科学版)》 *
黄懿雪;李国栋;郭长滨;: "成品油混输监测系统的关键技术问题及解决办法", 自动化博览, no. 07 *

Also Published As

Publication number Publication date
CN116822999B (en) 2023-12-05

Similar Documents

Publication Publication Date Title
Kou et al. Pairwise comparison matrix in multiple criteria decision making
CN106572493B (en) Rejecting outliers method and system in LTE network
CN110334726A (en) A method for identifying and repairing abnormal data of electric load based on density clustering and LSTM
CN111900731A (en) PMU-based power system state estimation performance evaluation method
CN115983721B (en) Cable production quality data management system based on big data of Internet of things
CN107942873B (en) An intelligent calculation and monitoring method for operating costs of a furniture manufacturing production line
CN102880809A (en) Polypropylene melt index on-line measurement method based on incident vector regression model
CN112348290A (en) River water quality prediction method, device, storage medium and equipment
CN102357934B (en) Quality monitor soft sensing method based on rubber mixing process
CN109816139B (en) Load subdivision based differential user baseline load prediction method, device and system
CN108664706B (en) Semi-supervised Bayesian Gaussian mixture model-based online estimation method for oxygen content of one-stage furnace in ammonia synthesis process
CN114782099A (en) Budget estimation and analysis method for construction cost
CN114189047A (en) False data detection and correction method for active power distribution network state estimation
CN108171002B (en) Polypropylene melt index prediction method based on semi-supervised hybrid model
CN105118293A (en) Road section traffic speed short-term prediction method in view of long-term equilibrium relationship
CN116307145B (en) Steel plate defect number prediction method based on Gauss-Poisson mixed variation self-encoder
CN110309873B (en) Integrated transient stability evaluation method and system based on Bayesian classification of power flow vectors
CN115542236B (en) Electric energy meter operation error estimation method and device
CN116822999B (en) Method and system for predicting monitoring density of oil product of oil mixing interface of finished oil pipeline
CN105939026B (en) Wind power undulate quantity probability Distribution Model method for building up based on mixing Laplace distributions
CN116307303A (en) Mechanism-data dual-drive oil mixing length prediction method and system for finished oil pipeline
CN105203327B (en) A kind of gas circuit measurement parameter selection method applied to engine air passage analysis
Mahmoud et al. Inferences for new Weibull-Pareto distribution based on progressively Type-II censored data
CN101281182B (en) Soft measurement method of component concentration of sodium aluminate solution
CN116151158B (en) A Method and System for Tracking Mixed Oil Interface of Refined Oil Pipeline Integrating Prior Knowledge

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