[go: up one dir, main page]

CN118857235B - Accurate detection method for installation angle of rotor flowmeter - Google Patents

Accurate detection method for installation angle of rotor flowmeter Download PDF

Info

Publication number
CN118857235B
CN118857235B CN202411320045.7A CN202411320045A CN118857235B CN 118857235 B CN118857235 B CN 118857235B CN 202411320045 A CN202411320045 A CN 202411320045A CN 118857235 B CN118857235 B CN 118857235B
Authority
CN
China
Prior art keywords
data point
installation inclination
smoothing
inclination angle
sequence
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
Application number
CN202411320045.7A
Other languages
Chinese (zh)
Other versions
CN118857235A (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.)
Ningbo Dongchi Measurement And Control Technology Co ltd
Original Assignee
Ningbo Dongchi Measurement And Control Technology Co ltd
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 Ningbo Dongchi Measurement And Control Technology Co ltd filed Critical Ningbo Dongchi Measurement And Control Technology Co ltd
Priority to CN202411320045.7A priority Critical patent/CN118857235B/en
Publication of CN118857235A publication Critical patent/CN118857235A/en
Application granted granted Critical
Publication of CN118857235B publication Critical patent/CN118857235B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C9/00Measuring inclination, e.g. by clinometers, by levels
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F25/00Testing or calibration of apparatus for measuring volume, volume flow or liquid level or for metering by volume
    • G01F25/10Testing or calibration of apparatus for measuring volume, volume flow or liquid level or for metering by volume of flowmeters
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Operations Research (AREA)
  • Fluid Mechanics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Indication And Recording Devices For Special Purposes And Tariff Metering Devices (AREA)
  • Measuring Volume Flow (AREA)

Abstract

本申请涉及倾角检测技术领域,具体涉及一种用于转子流量计安装倾角的精确检测方法。该方法包括:采集每个时刻的流量值和安装倾角,构成序列,对于安装倾角序列,通过平滑前后安装倾角差异和偏离程度差异获取平滑度;通过相邻数据点的安装倾角差异和偏离程度差异获取干扰度,获取流量序列的干扰度,基于两者的干扰度和序列值获取流量影响度,对干扰度优化;基于优化后干扰度和平滑度获取平滑因子的处理效果,并以此获取调节后的平滑因子,通过调节后的平滑因子对安装倾角进行平滑,完成精确检测。本申请提高了转子流量计安装倾角的检测精度与效率。

The present application relates to the technical field of inclination detection, and specifically to a method for accurately detecting the installation inclination of a rotor flowmeter. The method comprises: collecting the flow value and installation inclination at each moment to form a sequence, and for the installation inclination sequence, obtaining the smoothness by smoothing the installation inclination difference and the deviation difference before and after; obtaining the interference by the installation inclination difference and the deviation difference of adjacent data points, obtaining the interference of the flow sequence, obtaining the flow influence based on the interference and sequence value of the two, and optimizing the interference; obtaining the processing effect of the smoothing factor based on the optimized interference and smoothness, and thereby obtaining the adjusted smoothing factor, smoothing the installation inclination by the adjusted smoothing factor, and completing the accurate detection. The present application improves the detection accuracy and efficiency of the installation inclination of the rotor flowmeter.

Description

Accurate detection method for installation inclination angle of rotameter
Technical Field
The application relates to the technical field of inclination angle detection, in particular to an accurate detection method for a rotor flowmeter installation inclination angle.
Background
The rotameter measures the flow rate of fluid according to the throttling principle, and keeps the differential pressure on the upper and lower sides of the rotor constant by changing the flow area of the fluid, so the rotameter is also called a variable flow area constant differential pressure flowmeter. The traditional rotameter has higher requirement on installation during installation, and vertical installation is usually required, however, in the practical application process, the traditional rotameter is often limited by installation conditions, and the installation accuracy cannot meet the conditions, so that the measurement of the rotameter has deviation and even cannot work normally. Therefore, it is generally required to detect the installation inclination angle of the rotameter and compensate the flow measurement value according to the installation inclination angle of the rotameter and corresponding software, wherein the detection accuracy of the installation inclination angle of the rotameter directly affects the compensation amount of the flow measurement value, that is, affects the accuracy and efficiency of measuring the flow of the rotameter.
When the horizontal sensor is used for acquiring the installation dip angle of the rotameter, the acquired installation dip angle data is lower in accuracy due to the influence of surrounding vibration, temperature change and the like of the rotor flow, and further the compensation accuracy of the rotameter is lower. Therefore, a filtering algorithm is generally required to perform filtering processing on the inclination angle data acquired by the level sensor, so that the measurement accuracy and efficiency of the level sensor are improved. An exponentially weighted moving average filter is a common filtering algorithm, where the choice of a smoothing factor is critical to the performance of the filter, too large a smoothing factor would make the filter too sensitive to noise, and too small a smoothing factor would make the filter too slow to react to changes in the data.
Disclosure of Invention
In order to solve the technical problem of lower installation inclination angle measurement precision, the application provides an accurate detection method for the installation inclination angle of a rotameter, which adopts the following technical scheme:
The application provides an accurate detection method for the installation inclination angle of a rotameter, which comprises the following steps:
collecting the flow value and the installation inclination angle at each moment;
Respectively forming a flow sequence and an installation dip sequence by the flow value and the installation dip, and smoothing the installation dip sequence by a preset smoothing factor in the installation dip sequence to obtain an installation dip after the data point is smoothed;
obtaining a neighborhood data point of each data point, and obtaining the interference degree of the data point according to the difference of the installation inclination angle of the data point and the neighborhood data point and the difference of the deviation degree of the data point and the average value of the deviation degree of the neighborhood data point; in the flow sequence, the interference degree of the data points in the flow sequence is obtained according to the deviation degree of the data points and the neighborhood data points and the difference of the flow values;
In the installation dip sequence, comparing the smoothness of the data points with the optimized interference degree and the smoothness of the rest data points in the installation dip sequence with the optimized interference degree to obtain a processing effect value of the smoothing factor;
And smoothing the installation dip angle according to the regulated smoothing factor, and accurately detecting the installation dip angle.
In the scheme, the smoothness of the installation dip angle data point is built through the deviation degree of the data points in the installation dip angle sequence under the action of the smoothing factor and the change of the installation dip angle, meanwhile, the interference degree of the installation dip angle data point is built through the local data change of the installation dip angle data point and the deviation degree distribution characteristics, the flow influence degree of the installation dip angle data point and the flow data is built based on the influence of the fluid flow data measured by the rotameter on the installation dip angle, the interference degree is optimized based on the flow influence degree, the processing effect value of the smoothing factor is built in a self-adaptive mode based on the optimized interference degree and the smoothness, the acquisition of the optimal smoothing factor is completed based on the processing effect value, the optimized interference degree and the smoothness self-adaptive mode, the filtering effect of the exponentially weighted moving average filter is improved, and the detection precision and the efficiency of the installation dip angle of the rotameter are further improved.
In one embodiment, the method for obtaining the smoothness of each data point under different smoothness factors according to the difference of the front-back installation inclination angles and the difference of the deviation degree of the same data point is as follows:
Calculating the Z fraction difference of each data point before and after smoothing in the installation dip sequence as the difference of the deviation degree of the installation dip sequence;
Acquiring the smoothness of the data points under a smoothing factor according to the Z fraction difference before and after smoothing and the installation inclination difference before and after smoothing;
The smoothness of the data points under the smoothing factor is positively correlated with the Z-score difference before and after smoothing and the installation inclination angle difference before and after smoothing respectively, and is negatively correlated with the Z-score after smoothing.
In one embodiment, the method for obtaining the interference degree of each data point according to the difference of the installation inclination angle of the data point and the neighborhood data point and the difference of the deviation degree of the data point and the mean value of the deviation degree of the neighborhood data point is as follows:
In the installation dip sequence, selecting a plurality of data points closest to the data point as neighborhood data points of the data point;
Obtaining the interference degree of the data points according to the difference of the mean values of the Z scores of the data points and the Z scores of the neighborhood data points and the difference of the installation inclination angles of the data points and the neighborhood data points;
The disturbance degree of the data points is positively correlated with the difference of the installation inclination angles, and the difference of the mean values of the Z scores of the data points in the neighborhood is negatively correlated.
In one embodiment, the method for obtaining the flow influence degree of the data points according to the interference degree difference of the same data point in different sequences and the difference of the installation inclination angle and the flow value is as follows:
, indicating the degree of interference of the q-th data point in the installation dip sequence, Indicating the degree of interference of the q-th data point in the traffic sequence,The normalized value of the mounting tilt for the q-th data point in the mounting tilt sequence is shown,Represents the normalized value of the flow value for the q-th data point in the flow sequence,Represents an exponential function with a base of a natural constant,The flow effect of the qth data point is shown.
In one embodiment, in the installation inclination sequence, the method for optimizing the interference degree according to the influence of the flow to obtain the optimized interference degree is as follows:
, Indicating the degree of interference of the q-th data point of the installation dip sequence after optimization, Indicating the flow effect of the qth data point,The degree of interference of the q-th data point of the installation dip sequence is shown.
In one embodiment, the method for comparing the smoothness of the data point with the optimized disturbance degree and the smoothness of the rest data points in the installation inclination angle sequence with the optimized disturbance degree to obtain the processing effect value of the smoothing factor is as follows:
, A smoothing effect value representing the q-th data point of the installation dip sequence, Representing the maximum value of the processing effect of all data points in the installation dip sequence at the smoothing factor,Representing the minimum value of the processing effect of all data points in the installation dip sequence at the smoothing factor,A linear normalization function is represented and,The treatment effect value of the q-th data point of the installation dip sequence under the smoothing factor is represented.
In one embodiment, the method for adjusting the smoothing factor according to the processing effect value, the smoothness and the optimized interference degree to obtain the adjusted smoothing factor includes:
In the formula, Representing the smoothness of the q-th data point of the installation dip sequence at the current smoothing factor,Indicating the degree of interference of the q-th data point in the installation dip sequence after optimization,Representing the effect of processing the qth data point in the installation dip sequence at the current smoothing factor,As the current smoothing factor is used,Is the adjusted smoothing factor.
In one embodiment, the method for smoothing the installation inclination angle according to the adjusted smoothing factor is as follows:
Obtaining the installation inclination angle after the smoothing treatment by an exponential moving weighted average method;
The algorithm is input into the adjusted smoothing factors, the collected installation inclination angle and the installation inclination angle after the previous data point smoothing treatment.
In one embodiment, the method for obtaining the smoothness of the data points under the smoothing factor according to the Z-score difference before and after smoothing and the installation inclination angle difference before and after smoothing is as follows:
g represents a data point before smoothing, E represents a data point after smoothing, Representing the Z-score of the data point before smoothing,Representing the Z-score of the smoothed data point,Representing the installation tilt of the data point before smoothing,Representing the installation tilt angle of the smoothed data point,Representing the smoothness of the data point at the current smoothing factor.
In one embodiment, the method for obtaining the interference degree of the data points according to the difference of the mean values of the Z scores of the data points and the Z scores of the neighborhood data points and the difference of the installation inclination angles of the data points and the neighborhood data points is as follows:
, The Z-score representing the q-th data point, The mean of the Z-scores of all neighbor data points representing the qth data point,Represents the mounting tilt of the q-th data point,Representing the mounting inclination of the ith neighbor data point of the qth data point,Representing the number of data points of the field,The degree of interference of the q-th data point is shown.
The method has the advantages that the smoothness of the installation dip data point is built through the deviation degree of the data points in the installation dip sequence under the action of the smoothness factor and the change of the installation dip, meanwhile, the interference degree of the installation dip data point is built through the local data change of the installation dip data point and the deviation degree distribution characteristics, the flow influence degree of the installation dip data point and the flow data is built based on the influence of the fluid flow data measured by the rotameter on the installation dip, the interference degree is optimized based on the flow influence degree, the processing effect value of the smoothness factor is built in a self-adaptive mode based on the optimized interference degree and the smoothness, the acquisition of the optimal smoothness factor is completed in a self-adaptive mode based on the processing effect value, the optimized interference degree and the smoothness degree, the filtering effect of the exponentially weighted moving average filter is improved, and the detection accuracy and the detection efficiency of the installation dip of the rotameter are further improved.
Drawings
In order to more clearly illustrate the embodiments of the application or the technical solutions and advantages of the prior art, the following description will briefly explain the drawings used in the embodiments or the description of the prior art, and it is obvious that the drawings in the following description are only some embodiments of the application, and other drawings can be obtained according to the drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a method for accurately detecting a rotameter mounting tilt angle according to one embodiment of the present application;
FIG. 2 is a schematic illustration of the installation of a rotameter with a conical tube;
Fig. 3 is a comparison of the installation inclination of the application in different situations.
Detailed Description
In order to further describe the technical means and effects adopted by the present application to achieve the preset purpose, the following detailed description refers to the specific implementation, structure, characteristics and effects of the accurate detection method for the installation inclination angle of the rotameter according to the present application, which is provided by the present application, with reference to the accompanying drawings and the preferred embodiments. In the following description, different "one embodiment" or "another embodiment" means that the embodiments are not necessarily the same. Furthermore, the particular features, structures, or characteristics of one or more embodiments may be combined in any suitable manner.
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 application belongs.
An embodiment of a precise detection method for the installation inclination angle of a rotameter:
The following specifically describes a specific scheme of the accurate detection method for the installation inclination angle of the rotameter provided by the application with reference to the accompanying drawings.
Referring to fig. 1, a flowchart of a method for precisely detecting a mounting inclination angle of a rotameter according to an embodiment of the present application is shown, and the method includes the following steps:
and S001, collecting the flow value and the installation inclination angle at each moment.
The rotor flowmeter is placed in a conical tube universal to the rotor flowmeter, a horizontal sensor is prevented from being arranged on the bottom surface of the conical tube, the center of gravity direction of the rotor flowmeter is coincident with the tube axis direction of the conical tube, as shown in fig. 2, 1 in fig. 2 is the rotor universal to the rotor flowmeter, 2 is the horizontal sensor, and 3 is the conical tube universal to the rotor flowmeter. The installation inclination angle of the rotameter can be obtained through the horizontal sensor, and the flow value is obtained in real time through the rotameter. Wherein the flow value at this time is an uncompensated flow value. The installation tilt angle and flow rate values are collected every five seconds.
So far, the flow value and the installation inclination angle at each moment are obtained.
Step S002, the flow value and the installation dip angle respectively form a flow sequence and an installation dip angle sequence, the installation dip angle sequence is smoothed by a preset smoothing factor to obtain the smoothed installation dip angles of the data points in the installation dip angle sequence, and the smoothness of each data point in the installation dip angle sequence under different smoothing factors is obtained according to the difference of the installation dip angles before and after the same data point in the installation dip angle sequence is smoothed and the difference of the deviation degree.
And for the collected installation dip angle and flow rate values, taking one day as a period, so that all the installation dip angle and flow rate values collected from the starting time of each day to the current moment respectively form a sequence, and respectively recording the sequence as an installation dip angle sequence and a flow rate sequence.
Because the fluid flow is overlarge, the pipeline pressure is overlarge and local deformation or pipeline vibration is caused, the installation dip angle can fluctuate near a certain fixed value, the collected installation dip angle is inaccurate, the installation dip angle needs to be smoothed, the corresponding smoothness can be adaptively constructed through the change of the installation dip angle before and after the smoothing of the installation dip angle to represent the change of the installation dip angle under the action of the filter, and the smoothing factor parameters are regulated and controlled based on the smoothness.
Because the installation dip angle sequence and the flow sequence are both time sequences, an initial smoothing factor of the exponentially weighted moving average filter is set first, and in this embodiment, the value of the initial smoothing factor is 0.3, and the initial smoothing factor is the smoothing factor of the first data point of the sequence, and the smoothing factor of each data point is obtained by adjusting the smoothing factor of the last data point. The installation dip sequence is smoothed by an exponential weighted moving average method, and the smoothed installation dip of each data point in the installation dip sequence is obtained after the smoothing, wherein the exponential weighted moving average method is a known technology and the application is not described in detail.
The Z fraction of each data point before and after the smoothing in the installation dip sequence is calculated, the calculation of the Z fraction is a known technology, the Z fraction is used for obtaining the deviation degree of each data point relative to the average value, the greater the deviation degree difference before and after the smoothing is, the more the installation dip with larger fluctuation is described to be corrected, the better the smoothing effect is, and the greater the installation dip difference before and after the smoothing is, the greater the smoothing effect is described. Based on this, the smoothness of each data point under the smoothing factor is obtained.
The smoothness of the data points under the smoothing factor is positively correlated with the Z-score difference before and after smoothing and the installation inclination angle difference before and after smoothing respectively, and is negatively correlated with the Z-score after smoothing.
It should be noted that positive correlation indicates that one variable grows and the other variable also grows, and that the directions of the two variables change are the same, and that when one variable changes from large to small or from small to large, the other variable also changes from large to small or from small to large, and the specific relation is determined by practical application, and the application is not limited in particular.
It should be noted that, the negative correlation indicates that one variable increases and the other variable decreases, and the directions of the two variable changes are opposite, when one variable changes from large to small or from small to large, the other variable also changes from small to large or from large to small, and the specific relation is determined by practical application, and the application is not limited in particular.
Preferably, in the present embodiment, the expression of the smoothness under the current smoothing factor is:
Where G represents the data point before smoothing, E represents the data point after smoothing, Representing the Z-score of the data point before smoothing,Representing the Z-score of the smoothed data point,Representing the installation tilt of the data point before smoothing,Representing the installation tilt angle of the smoothed data point,Representing the smoothness of the data point at the current smoothing factor.
Preferably, in another embodiment of the present application, the expression of the smoothness under the current smoothing factor is:
The larger the Z fraction of the smoothed data point in the installation dip sequence, the farther the deviation degree of the smoothed data point is, the worse the smoothing effect is, the larger the Z fraction difference of the smoothed data point is, the better the smoothing effect of the installation dip under the current smoothing factor is, the larger the difference of the installation dip is, the larger the change of the smoothed data point is, and the better the smoothing effect is.
Thus, the smoothness of each data point under the current smoothing factor is obtained.
Step S003, obtaining a neighborhood data point of each data point, obtaining the interference degree of the data point according to the difference of the installation inclination angle of the data point and the neighborhood data point and the difference of the deviation degree of the data point and the deviation degree mean value of the neighborhood data point, obtaining the interference degree of the data point in the flow sequence according to the deviation degree of the data point and the neighborhood data point and the difference of the flow value in the flow sequence, obtaining the flow influence degree of the data point according to the interference degree difference of the same data point in different sequences and the difference of the installation inclination angle and the flow value, and obtaining the optimized interference degree according to the flow influence degree in the installation inclination angle sequence.
The smoothness of each data point in the installation dip sequence under the current smoothing factor is obtained through the steps. Since the smoothness is not as high as better for different data points, when the data points in the installation dip sequence are normal data points (i.e. data points with smaller interference), the greater the smoothness is, the greater the influence weight of the historical data in the current filter is indicated, and the current normal data points are excessively smoothed, and when the data points in the installation dip sequence are abnormal data points (i.e. data points with larger interference), the greater the smoothness is, the better the noise removal effect of the current filter is indicated.
Therefore, data points in the installation dip sequence need to be analyzed, and the interference degree of each data point in the installation dip sequence is obtained according to the difference of the installation dip of the data point and the data point in the neighborhood of the data point and the Z fraction difference.
For each data point in the installation dip sequence, selecting M data points closest to the data point in the installation dip sequence as neighborhood data points, and if M data points do not exist in the installation dip sequence, taking all the data points, wherein the value of M is 11.
The degree of disturbance of the data points is positively correlated with the difference in the mounting inclination angles and with the difference in the Z fractions, respectively.
Preferably, in this embodiment, the interference of the data points in the installation dip sequence is:
, The Z-score representing the q-th data point, The mean of the Z-scores of all neighbor data points representing the qth data point,Represents the mounting tilt of the q-th data point,Representing the mounting inclination of the ith neighbor data point of the qth data point,Representing the number of data points of the field,The degree of interference of the q-th data point is shown.
Preferably, in another embodiment of the present application, the degree of interference of the data points in the installation dip sequence is:
The larger the Z score of the data point is, the smaller the average value of the Z scores of all the neighborhood data points is, namely, the larger the difference between the Z scores of the data point and the neighborhood data point is, the larger the deviation of the data point relative to the neighborhood data point is, namely, the larger the independence of the data point is, the more likely the change of the installation inclination angle of the data point is an independent change, and the larger the disturbance degree of the data point is. The larger the difference between the installation inclination angle of the data point and the installation inclination angle of the neighborhood data point is, the abrupt change of the installation inclination angle of the data point is indicated, the larger the installation inclination angle is compared with the normal fluctuation difference, and the interference degree of the data point is larger.
The interference degree of each data point in the installation dip angle sequence is obtained through the steps, and then for each flow sequence, the interference degree of each data point in the flow sequence is obtained through the differences of the Z fraction and the flow value of the data point in the flow sequence and the flow value and the Z fraction average value of the neighborhood data point.
Because the installation inclination can receive the influence of the flow of fluid when gathering, when the flow is too big, the pressure of fluid to the pipeline changes, and then leads to the local load of pipeline too big, appears the phenomenon of local deformation, so when the flow is too big, can make the installation inclination take place great change, and then appear the too big phenomenon of interference degree of data point in the installation inclination sequence.
Based on the analysis, the flow value acquired by the rotameter is further analyzed, and because the flow value and the installation dip angle are acquired at each moment, the data points of the installation dip angle sequence and the data points in the flow sequence are in one-to-one correspondence at the same moment, so that the flow influence degree of each data point is acquired according to the interference degree difference of the corresponding data points of different sequences and the difference of the installation dip angle and the flow value of the data points.
Preferably, in this embodiment, first, a normalized value of the installation inclination angle and the flow rate value of each data point is obtained, and the flow rate influence is expressed as follows:
, indicating the degree of interference of the q-th data point in the installation dip sequence, Indicating the degree of interference of the q-th data point in the traffic sequence,The normalized value of the mounting tilt for the q-th data point in the mounting tilt sequence is shown,Represents the normalized value of the flow value for the q-th data point in the flow sequence,Represents an exponential function with a base of a natural constant,The flow effect of the qth data point is shown.
When the flow data is too large, the installation inclination angle can be changed due to local deformation caused by the excessive local load of the pipeline, so that the flow influence degree reflects the possibility that the abnormality of the installation inclination angle is influenced by the flow data. If the difference between the disturbance degrees of the data points corresponding to the installation inclination angle sequence and the flow sequence is smaller, the flow influence degree is larger, the change of the installation inclination angle and the flow value is reflected on the whole, the change of the numerical value is ignored, and therefore the flow influence degree of the data points is obtained by combining the difference of the numerical values.
Optimizing the interference degree of the data points in the installation dip angle sequence through the flow influence of the data points, and obtaining the optimized interference degree, wherein the expression is as follows:, the disturbance degree after the q-th data point of the installation dip sequence is optimized is shown.
So far, the interference degree of each data point in the installation dip sequence after optimization is obtained.
And S004, comparing the smoothness of the data points with the optimized interference degree and the smoothness of the rest data points in the installation dip angle sequence with the optimized interference degree in the installation dip angle sequence to obtain a processing effect value of the smoothing factor, and adjusting the smoothing factor according to the processing effect value, the smoothness and the optimized interference degree to obtain the adjusted smoothing factor.
Through the steps, the smoothness of the data points and the optimized interference degree are obtained, and the processing effect of the smoothing factor is constructed based on the smoothness and the optimized interference degree.
For each data point, the smoothness of the data point indicates that the data point has larger change under the action of the smoothness factor, the optimized interference degree of the data point indicates that the data point is more likely to be a data point acquired for receiving interference, the product of the smoothness of each data point and the optimized interference degree is calculated and recorded as a smoothness effect value, the smoothness effect value characterizes the processing effect of the smoothness factor, when the smoothness effect value is larger, the smoothness effect of the installation dip angle data point is larger when the smoothness effect value is larger, and when the smoothness effect value is smaller, the smoothness effect of the installation dip angle is smaller when the smoothness effect value is not received, and the content remained after the smoothness is more.
And obtaining the processing effect value of the smoothing factor by installing the smoothing effect value of each data point in the dip angle sequence.
Preferably, in the present embodiment, the expression of the processing effect value of the smoothing factor is:
, A smoothing effect value representing the q-th data point of the installation dip sequence, Representing the maximum value of the processing effect of all data points in the installation dip sequence at the smoothing factor,Representing the minimum value of the processing effect of all data points in the installation dip sequence at the smoothing factor,A linear normalization function is represented and,The treatment effect value of the q-th data point of the installation dip sequence under the smoothing factor is represented.
The processing effect value of each data point of the installation dip angle sequence under the smoothing factor is obtained through the steps, and the smoothing factor is adaptively adjusted based on the processing effect value, and the adjustment mode is as follows:
In the formula, Representing the smoothness of the q-th data point of the installation dip sequence at the current smoothing factor,Indicating the degree of interference of the q-th data point in the installation dip sequence after optimization,Representing the effect of processing the qth data point in the installation dip sequence at the current smoothing factor,As the current smoothing factor is used,Is the adjusted smoothing factor.
When the processing effect of the data point under the smoothing factor is weaker, the smoothing factor is corrected, and when the interference degree after optimization is greater than or equal to the smoothness, the interference suffered by the data point in the installation dip angle sequence is larger, and at the moment, the installation dip angle after the adjustment of the last data point is more important, and the current smoothing factor needs to be reduced. And when the interference degree after optimization is smaller than the smoothness, the data points in the current installation dip sequence are excessively smoothed, and more attention is needed to be paid to the installation dip acquired by the data points, so that the smoothing factor is needed to be increased.
Thus, the adjusted smoothing factor is obtained.
And step S005, smoothing the installation inclination angle according to the regulated smoothing factor, and accurately detecting the installation inclination angle.
For each data point in the installation dip angle sequence, the smoothing factor of the previous data point is used as the current smoothing factor, the current smoothing factor is regulated through the steps, the installation dip angle after the smoothing treatment is obtained through the regulated smoothing factor through an exponential moving weighted average method, and the accurate detection of the installation dip angle is completed. The input is the adjusted smoothing factor, the collected installation inclination angle and the installation inclination angle after the previous data point smoothing treatment. Unfiltered mounting angles, original filtered mounting angles and optimized filtered mounting angles of the present application are shown in fig. 3.
It should be noted that the foregoing embodiments are merely illustrative of the technical solutions of the present application, and not restrictive, and although the present application has been described in detail with reference to the foregoing embodiments, it should be understood by those skilled in the art that the technical solutions described in the foregoing embodiments may be modified or some of the technical features may be replaced equivalently, and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the scope of the technical solutions of the embodiments of the present application and are included in the protection scope of the present application.
In this specification, each embodiment is described in a progressive manner, and identical and similar parts of each embodiment are all referred to each other, and each embodiment mainly describes differences from other embodiments.

Claims (10)

1.用于转子流量计安装倾角的精确检测方法,其特征在于,该方法包括以下步骤:1. A method for accurately detecting the installation inclination angle of a rotor flowmeter, characterized in that the method comprises the following steps: 采集每个时刻的流量值和安装倾角;Collect the flow value and installation inclination at each moment; 将流量值和安装倾角分别构成流量序列和安装倾角序列,在安装倾角序列中,通过预设平滑因子对安装倾角序列进行平滑获取数据点平滑后的安装倾角;根据同一个数据点平滑前后安装倾角的差异和偏离程度的差异获取每个数据点在不同平滑因子下的平滑度;The flow value and the installation inclination angle constitute a flow sequence and an installation inclination angle sequence respectively. In the installation inclination angle sequence, the installation inclination angle sequence is smoothed by a preset smoothing factor to obtain the installation inclination angle after smoothing of the data point; the smoothness of each data point under different smoothing factors is obtained according to the difference in the installation inclination angle before and after smoothing and the difference in the degree of deviation of the same data point; 获取每个数据点的邻域数据点,根据数据点与其邻域数据点的安装倾角差异以及数据点的偏离程度与邻域数据点偏离程度均值的差异获取数据点的干扰度;在流量序列中,根据数据点和邻域数据点偏离程度和流量值的差异获取流量序列中数据点的干扰度;根据同一个数据点在不同序列的干扰度差异以及安装倾角与流量值的差异获取数据点的流量影响度;在安装倾角序列中,根据流量影响度对干扰度优化获取优化后的干扰度;Obtain the neighborhood data points of each data point, and obtain the interference degree of the data point according to the difference in installation inclination between the data point and its neighborhood data points and the difference between the deviation degree of the data point and the average deviation degree of the neighborhood data points; in the flow sequence, obtain the interference degree of the data point in the flow sequence according to the difference in deviation degree and flow value between the data point and the neighborhood data points; obtain the flow influence degree of the data point according to the difference in interference degree of the same data point in different sequences and the difference between the installation inclination angle and the flow value; in the installation inclination sequence, optimize the interference degree according to the flow influence degree to obtain the optimized interference degree; 在安装倾角序列中,将数据点的平滑度和优化后的干扰度和安装倾角序列中其余数据点的平滑度和优化后的干扰度比较获取平滑因子的处理效果值;根据处理效果值、平滑度以及优化后的干扰度对平滑因子调节获取调节后的平滑因子;In the installation inclination sequence, the smoothness and optimized interference of the data point are compared with the smoothness and optimized interference of the remaining data points in the installation inclination sequence to obtain a processing effect value of the smoothing factor; the smoothing factor is adjusted according to the processing effect value, the smoothness and the optimized interference to obtain an adjusted smoothing factor; 根据调节后的平滑因子对安装倾角进行平滑,对安装倾角进行精确检测。The installation inclination angle is smoothed according to the adjusted smoothing factor, and the installation inclination angle is accurately detected. 2.如权利要求1所述的用于转子流量计安装倾角的精确检测方法,其特征在于,所述根据同一个数据点平滑前后安装倾角的差异和偏离程度的差异获取每个数据点在不同平滑因子下的平滑度的方法为:2. The accurate detection method for the installation inclination angle of the rotor flowmeter according to claim 1 is characterized in that the method for obtaining the smoothness of each data point under different smoothing factors according to the difference in the installation inclination angle before and after smoothing and the difference in the degree of deviation of the same data point is: 计算平滑前后每个数据点在安装倾角序列的Z分数差异作为安装倾角序列的偏离程度的差异;The difference in the Z score of each data point in the installation inclination sequence before and after smoothing is calculated as the difference in the degree of deviation of the installation inclination sequence; 根据平滑前后的Z分数差异和平滑前后的安装倾角差异获取数据点在平滑因子下的平滑度;The smoothness of the data points under the smoothing factor is obtained according to the difference in Z scores before and after smoothing and the difference in the installation inclination angle before and after smoothing; 数据点在平滑因子下的平滑度分别与平滑前后的Z分数差异和平滑前后的安装倾角差异呈正相关,与平滑后的Z分数呈负相关。The smoothness of the data points under the smoothing factor is positively correlated with the difference in Z scores before and after smoothing and the difference in mounting inclination angles before and after smoothing, and negatively correlated with the Z score after smoothing. 3.如权利要求2所述的用于转子流量计安装倾角的精确检测方法,其特征在于,所述获取每个数据点的邻域数据点,根据数据点与其邻域数据点的安装倾角差异以及数据点的偏离程度与邻域数据点偏离程度均值的差异获取数据点的干扰度的方法为:3. The accurate detection method for the installation inclination angle of the rotor flowmeter according to claim 2 is characterized in that the method of obtaining the neighborhood data points of each data point and obtaining the interference degree of the data point according to the difference in the installation inclination angle between the data point and its neighborhood data points and the difference between the deviation degree of the data point and the mean deviation degree of the neighborhood data points is: 在安装倾角序列中,对于数据点选取距离其最近的若干个数据点作为其的邻域数据点;In the installation inclination sequence, for a data point, several data points closest to it are selected as its neighborhood data points; 根据数据点的Z分数和邻域数据点Z分数的均值的差异以及数据点和邻域数据点的安装倾角差异获取数据点的干扰度;Obtain the interference degree of the data point according to the difference between the Z score of the data point and the mean Z score of the neighboring data points and the difference in the installation inclination angle between the data point and the neighboring data points; 数据点的干扰度分别与安装倾角差异呈正相关,Z分数和邻域数据点Z分数的均值的差异呈负相关。The interference degree of the data points was positively correlated with the difference in installation inclination angles, and the difference in Z scores and the mean Z scores of neighboring data points was negatively correlated. 4.如权利要求1所述的用于转子流量计安装倾角的精确检测方法,其特征在于,所述根据同一个数据点在不同序列的干扰度差异以及安装倾角与流量值的差异获取数据点的流量影响度的方法为:4. The accurate detection method for the installation inclination angle of a rotor flowmeter according to claim 1 is characterized in that the method for obtaining the flow influence of a data point according to the difference in interference degree of the same data point in different sequences and the difference between the installation inclination angle and the flow value is: 表示安装倾角序列中第q个数据点的干扰度,表示流量序列中第q个数据点的干扰度,表示安装倾角序列中第q个数据点的安装倾角归一化后的值,表示流量序列中第q个数据点的流量值归一化后的值,表示以自然常数为底的指数函数,表示第q个数据点的流量影响度。 , represents the interference degree of the qth data point in the installation inclination sequence, represents the interference degree of the qth data point in the traffic sequence, represents the normalized value of the installation inclination angle of the qth data point in the installation inclination angle sequence, Represents the normalized value of the flow value of the qth data point in the flow sequence, represents an exponential function with a natural constant as base, Indicates the traffic influence of the qth data point. 5.如权利要求1所述的用于转子流量计安装倾角的精确检测方法,其特征在于,所述在安装倾角序列中,根据流量影响度对干扰度优化获取优化后的干扰度的方法为:5. The accurate detection method for the installation inclination angle of the rotor flowmeter according to claim 1 is characterized in that, in the installation inclination angle sequence, the method for optimizing the interference degree according to the flow influence degree to obtain the optimized interference degree is: 表示安装倾角序列第q个数据点优化后的干扰度,表示第q个数据点的流量影响度,表示安装倾角序列第q个数据点的干扰度。 , It represents the interference degree after optimization of the qth data point in the installation inclination sequence, represents the traffic impact of the qth data point, Represents the interference degree of the qth data point in the installation inclination sequence. 6.如权利要求1所述的用于转子流量计安装倾角的精确检测方法,其特征在于,所述将数据点的平滑度和优化后的干扰度和安装倾角序列中其余数据点的平滑度和优化后的干扰度比较获取平滑因子的处理效果值的方法为:6. The accurate detection method for the installation inclination angle of a rotor flowmeter according to claim 1 is characterized in that the method of comparing the smoothness and optimized interference degree of a data point with the smoothness and optimized interference degree of other data points in the installation inclination angle sequence to obtain the processing effect value of the smoothing factor is: 表示安装倾角序列的第q个数据点的平滑效果值,表示安装倾角序列中所有数据点在平滑因子下的处理效果的最大值, 表示安装倾角序列中所有数据点在平滑因子下的处理效果的最小值,表示线性归一化函数,表示安装倾角序列的第q个数据点在平滑因子下的处理效果值。 , represents the smoothing effect value of the qth data point of the installed inclination sequence, It represents the maximum value of the processing effect of all data points in the installation inclination sequence under the smoothing factor. It represents the minimum value of the processing effect of all data points in the installation inclination sequence under the smoothing factor. represents the linear normalization function, It represents the treatment effect value of the qth data point of the installation inclination angle sequence under the smoothing factor. 7.如权利要求1所述的用于转子流量计安装倾角的精确检测方法,其特征在于,所述根据处理效果值、平滑度以及优化后的干扰度对平滑因子调节获取调节后的平滑因子的方法为:7. The accurate detection method for the installation inclination angle of the rotor flowmeter according to claim 1 is characterized in that the method of adjusting the smoothing factor according to the processing effect value, the smoothness and the optimized interference degree to obtain the adjusted smoothing factor is: 式中,表示安装倾角序列第q个数据点在当前平滑因子下的平滑度,表示安装倾角序列中第q个数据点优化后的干扰度,表示安装倾角序列中第q个数据点在当前平滑因子下的处理效果,为当前平滑因子,为调节后的平滑因子。In the formula, Indicates the smoothness of the qth data point in the installation inclination sequence under the current smoothing factor, It represents the interference degree after optimization of the qth data point in the installation inclination sequence, It indicates the processing effect of the qth data point in the installation inclination sequence under the current smoothing factor. is the current smoothing factor, is the adjusted smoothing factor. 8.如权利要求1所述的用于转子流量计安装倾角的精确检测方法,其特征在于,所述根据调节后的平滑因子对安装倾角进行平滑的方法为:8. The accurate detection method for the installation inclination angle of the rotor flowmeter according to claim 1 is characterized in that the method for smoothing the installation inclination angle according to the adjusted smoothing factor is: 通过指数移动加权平均法获取平滑处理后的安装倾角;The installation inclination angle after smoothing is obtained by using the exponential moving weighted average method; 算法的输入为调节后的平滑因子,采集的安装倾角以及前一个数据点平滑处理后的安装倾角。The inputs to the algorithm are the adjusted smoothing factor, the collected installation inclination angle, and the installation inclination angle after smoothing of the previous data point. 9.如权利要求2所述的用于转子流量计安装倾角的精确检测方法,其特征在于,所述根据平滑前后的Z分数差异和平滑前后的安装倾角差异获取数据点在平滑因子下的平滑度的方法为:9. The accurate detection method for the installation inclination angle of the rotor flowmeter according to claim 2, characterized in that the method for obtaining the smoothness of the data point under the smoothing factor according to the difference in Z scores before and after smoothing and the difference in the installation inclination angle before and after smoothing is: ,G表示平滑前的数据点,E表示平滑后的数据点,表示平滑前数据点的Z分数,表示平滑后数据点的Z分数,表示平滑前数据点的安装倾角,表示平滑后数据点的安装倾角,表示数据点在当前平滑因子下的平滑度。 , G represents the data points before smoothing, and E represents the data points after smoothing. represents the Z score of the data point before smoothing, represents the Z score of the smoothed data point, represents the installation inclination angle of the data point before smoothing, represents the installation inclination angle of the smoothed data point, Indicates the smoothness of the data point under the current smoothing factor. 10.如权利要求3所述的用于转子流量计安装倾角的精确检测方法,其特征在于,所述根据数据点的Z分数和邻域数据点Z分数的均值的差异以及数据点和邻域数据点的安装倾角差异获取数据点的干扰度的方法为:10. The accurate detection method for the installation inclination angle of the rotor flowmeter according to claim 3, characterized in that the method for obtaining the interference degree of the data point according to the difference between the Z score of the data point and the mean value of the Z score of the neighboring data points and the difference between the installation inclination angle of the data point and the neighboring data points is: 表示第q个数据点的Z分数,表示第q个数据点的所有邻域数据点的Z分数的均值,表示第q个数据点的安装倾角,表示第q个数据点的第i个邻域数据点的安装倾角,表示领域数据点的数量,表示第q个数据点的干扰度。 , represents the Z score of the qth data point, represents the mean of the Z scores of all neighboring data points of the qth data point, represents the installation inclination angle of the qth data point, represents the installation inclination of the ith neighboring data point of the qth data point, represents the number of domain data points, Represents the interference degree of the qth data point.
CN202411320045.7A 2024-09-23 2024-09-23 Accurate detection method for installation angle of rotor flowmeter Active CN118857235B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202411320045.7A CN118857235B (en) 2024-09-23 2024-09-23 Accurate detection method for installation angle of rotor flowmeter

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202411320045.7A CN118857235B (en) 2024-09-23 2024-09-23 Accurate detection method for installation angle of rotor flowmeter

Publications (2)

Publication Number Publication Date
CN118857235A CN118857235A (en) 2024-10-29
CN118857235B true CN118857235B (en) 2025-02-28

Family

ID=93169450

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202411320045.7A Active CN118857235B (en) 2024-09-23 2024-09-23 Accurate detection method for installation angle of rotor flowmeter

Country Status (1)

Country Link
CN (1) CN118857235B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA261957A (en) * 1926-06-22 Brodhead Linderman Garrett Meter mechanism
CN107255438A (en) * 2017-04-25 2017-10-17 中国计量大学 A kind of micro- angle change measurement apparatus of simple ultrasonic water meter reflector plate

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4675523B2 (en) * 2001-09-05 2011-04-27 日新製鋼株式会社 Method for calculating the terrace length of the raw material deposition layer at the top of the blast furnace furnace
JP5397836B2 (en) * 2009-12-10 2014-01-22 横河電機株式会社 Ultrasonic flow meter
CN107830903B (en) * 2017-10-25 2020-02-18 傅古月 Microwave solid flowmeter detecting system and device
CN212133746U (en) * 2020-04-28 2020-12-11 宁波东驰测控技术有限公司 Corrosion-resistant metal rotor flowmeter
CN114529521B (en) * 2022-01-26 2024-06-28 湖北飞来钟粮油设备有限公司 Method and device for controlling inclination angle of grain husking machine, grain husking machine and storage medium
CN115067931B (en) * 2022-06-28 2025-04-22 北京巨驰医药技术有限公司 Valsalva test automatic detection method, device, system and storage medium

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA261957A (en) * 1926-06-22 Brodhead Linderman Garrett Meter mechanism
CN107255438A (en) * 2017-04-25 2017-10-17 中国计量大学 A kind of micro- angle change measurement apparatus of simple ultrasonic water meter reflector plate

Also Published As

Publication number Publication date
CN118857235A (en) 2024-10-29

Similar Documents

Publication Publication Date Title
CN118954657B (en) Mine water micro-sand separation system and method based on dynamic filtering algorithm
CN118211082B (en) Oil filter element residual life prediction method and system based on data analysis
CN115290316B (en) Fault diagnosis method for eccentric rotary valve
CN114817228A (en) A method for removing outliers from long-term monitoring data of lake water level
CN116242443B (en) Dynamic metering method of ultrasonic metering instrument
CN118348814B (en) High-precision proportional valve control method and system based on robust control
CN118502255B (en) Remote monitoring method for multi-axis numerical control equipment
CN117848438A (en) High-precision measuring method for gas flow and thermal gas mass flowmeter
CN114638435B (en) Diesel engine security parameter prediction method based on data driving
CN118857235B (en) Accurate detection method for installation angle of rotor flowmeter
CN119848436B (en) Hydraulic and hydroelectric engineering data management method based on BIM information processing
CN118629160A (en) A landslide early warning method based on adaptive local regression
CN113901379A (en) A dynamic online fast processing method for real-time data at the edge
CN117171685A (en) Operation monitoring method of turbidity measurement system
CN112003891B (en) Multi-sensor data fusion method for intelligent networked vehicle controller
CN118499233A (en) Vacuum pump fault alarm method and system
CN118981737A (en) IoT device operation monitoring method and system based on edge computing
CN117268743B (en) Fault diagnosis method for proportional flow valve
CN115015810B (en) An electrical connection defect diagnosis system for inductive capacitors
CN110852397A (en) An adaptive signal fusion method and system based on relative fluctuation
CN119293429B (en) A method for improving the accuracy of gyro array data for stable platforms
CN119062426B (en) Differential pressure sensor control method for measuring engine aftertreatment differential pressure
CN119916698B (en) Anti-interference six-dimensional sensor acquisition system
CN119861649A (en) Intelligent control system for switching pulse signals
KR102718498B1 (en) Signal processing device and signal processing method

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