Piping stress-aware mass flowmeter measurement compensation method and system

Abstract

The application provides a pipeline stress-aware mass flowmeter measurement compensation method and system, relates to the technical field of instrument correction, and the method comprises the following steps: obtaining an original stress data set; constructing a stress characteristic parameter set; constructing a sensor parameter set; calculating stress change amounts of the pipeline in axial, radial and torsional directions based on the stress characteristic parameter set and the sensor parameter set, and obtaining a stress change data set; establishing a coupling model of pipeline stress change on sensor vibration frequency and phase signal by using the stress change data set and the sensor parameter set, and calculating a corrected mass flow value; and according to the comparison result of the corrected mass flow value and a calibration point in a long-term operation process, applying a parameter optimization algorithm to adaptively adjust the coupling coefficient, so that the coupling model is optimized. The technical problem of measurement deviation caused by external stress, structural deformation and dynamic load in the prior art is solved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of instrument calibration, in particular to a mass flowmeter measurement compensation method and system combined with pipeline stress sensing. BACKGROUND

[0002] The mass flowmeter has been widely applied to the flow monitoring and metering links of process industries such as petroleum chemical industry, energy, metallurgy, food and medicine, because it has the advantages of high measurement accuracy, wide applicable medium range, and being not affected by fluid temperature and pressure changes, etc. In a complex pipeline system, the mass flowmeter is usually directly installed on a pipe section bearing mechanical stress, fluid impact and environmental vibration, and its output signal including parameters such as vibration frequency and phase difference is used to calculate instantaneous mass flow and cumulative flow, and is a key measurement unit in industrial flow metering system.

[0003] However, under long-term running conditions, the measurement accuracy of the mass flowmeter is affected by various non-ideal factors such as pipeline stress, installation strain and external vibration. The change of pipeline stress will cause the inherent frequency of the sensing pipe to deviate, resulting in distortion of the phase difference of Coriolis vibration response, thereby generating systematic measurement error. In addition, external stress often has randomness and time variability, and the traditional flowmeter compensation method mainly relies on indirect correction of temperature, pressure or vibration signals, which is difficult to effectively reflect the real stress state of the sensing pipe. At the same time, the existing device usually does not establish a physical coupling model between the pipeline stress and the change of the sensing signal, and also lacks a compensation mechanism for the structural characteristics and attitude differences of the flowmeter, resulting in insufficient stability and repeatability of the measurement results under complex load conditions.

[0004] In summary, although the existing mass flowmeter has high measurement accuracy under normal working conditions, it still has measurement deviation problems caused by external stress, structural deformation and dynamic load in a complex pipeline environment. The existing technology lacks a measurement correction method that can sense the change of pipeline stress in real time and realize coupling compensation based on stress characteristics and structural parameters, which is difficult to ensure the measurement accuracy and long-term reliability of the flowmeter under the coupling action of multiple stress fields. Therefore, there is an urgent need for a method to solve the existing technical limitations. SUMMARY

[0005] The present disclosure provides a mass flowmeter measurement compensation method and system combined with pipeline stress sensing to solve the technical problem of measurement deviation caused by external stress, structural deformation and dynamic load in the prior art.

[0006] According to a first aspect of the present disclosure, a mass flowmeter measurement compensation method combined with pipeline stress sensing is provided, comprising:

[0007] A stress sensor array is arranged at a key position of a pipe segment of the mass flow meter, and stress distribution data in axial, radial and torsional directions of the pipeline are collected in real time to obtain an original stress data set;

[0008] The original stress data set is preprocessed by filtering and time series smoothing, average stress, stress gradient, strain rate and vibration main frequency components are extracted, and a stress characteristic parameter set is constructed;

[0009] The geometric size and material performance of the mass flow meter, and the spatial attitude information of the stress sensor are obtained, and a sensor parameter set is constructed;

[0010] Based on the stress characteristic parameter set and the sensor parameter set, the stress change amount of the pipeline in the axial, radial and torsional directions is calculated, and a stress change data set is obtained;

[0011] The stress change data set and the sensor parameter set are used to establish a coupling model of the pipeline stress change to the vibration frequency and phase signal of the sensing pipe, and a corrected mass flow value is calculated;

[0012] According to the comparison result of the corrected mass flow value and the calibration point in the long-term operation process, the coupling coefficient is adaptively adjusted by applying a parameter optimization algorithm, and the coupling model is optimized.

[0013] According to the second aspect of the present disclosure, a mass flow meter measurement compensation system combined with pipeline stress sensing is provided, comprising:

[0014] The data acquisition module is used for arranging a stress sensor array at a key position of a pipe segment of the mass flow meter, collecting stress distribution data in axial, radial and torsional directions of the pipeline in real time, and obtaining an original stress data set;

[0015] The stress feature extraction module is used for filtering and time series smoothing preprocessing of the original stress data set, extracting average stress, stress gradient, strain rate and vibration main frequency components, and constructing a stress characteristic parameter set;

[0016] The sensor parameter construction module is used for obtaining the geometric size and material performance of the mass flow meter, and the spatial attitude information of the stress sensor, and constructing a sensor parameter set;

[0017] The stress change calculation module is used for calculating the stress change amount of the pipeline in the axial, radial and torsional directions based on the stress characteristic parameter set and the sensor parameter set, and obtaining a stress change data set;

[0018] The coupling model correction module is used for establishing a coupling model of the pipeline stress change to the vibration frequency and phase signal of the sensing pipe by using the stress change data set and the sensor parameter set, and calculating a corrected mass flow value.

[0019] The model self-adaptive optimization module is configured to apply a parameter optimization algorithm to adaptively adjust the coupling coefficient according to a comparison result of the corrected mass flow value and the calibration point during long-term operation, so as to realize optimization of the coupling model.

[0020] The one or more technical solutions provided in the disclosure have at least the following technical effects or advantages: stress sensor arrays are arranged at key positions of a pipe section of a mass flowmeter, stress distribution data in axial, radial and torsional directions of the pipeline are collected in real time to obtain an original stress data set; the original stress data set is preprocessed by filtering and time series smoothing, average stress, stress gradient, strain rate and vibration main frequency components are extracted to construct a stress characteristic parameter set; geometric dimensions and material properties of the mass flowmeter and spatial posture information of the stress sensor are obtained to construct a sensor parameter set; based on the stress characteristic parameter set and the sensor parameter set, stress variation amounts of the pipeline in the axial, radial and torsional directions are calculated to obtain a stress variation data set; the stress variation data set and the sensor parameter set are used to establish a coupling model of stress variation of the pipeline to vibration frequency and phase signal of the sensing pipe, and a corrected mass flow value is calculated; based on a comparison result of the corrected mass flow value and the calibration point during long-term operation, a parameter optimization algorithm is applied to adaptively adjust the coupling coefficient, so as to realize optimization of the coupling model. The technical problem of measurement deviation caused by external stress, structural deformation and dynamic load in the prior art is solved. The technical effect of realizing dynamic perception of the influence of external stress, calculation modeling and flow correction is achieved, and the measurement consistency and precision of the mass flowmeter under complex working conditions are improved.

[0021] The above description is only a summary of the technical solutions of the present application. In order to more clearly understand the technical means of the present application, the specific embodiments of the present application can be implemented according to the content of the description, and in order to make the above and other purposes, features and advantages of the present application more obvious and easy to understand, the following specific embodiments of the present application are described. BRIEF DESCRIPTION OF DRAWINGS

[0022] In order to more clearly illustrate the technical solutions in the disclosure or the prior art, the following will briefly introduce the drawings needed in the embodiments or the prior art description. Obviously, the drawings in the following description are only exemplary, and those skilled in the art can obtain other drawings according to the provided drawings without creating any creative labor.

[0023] Figure 1 The flowchart of the mass flowmeter measurement compensation method combined with pipeline stress perception provided in the embodiments of the present application;

[0024] Figure 2A structural schematic diagram of a quality flow meter measurement compensation system combined with pipeline stress perception provided by the embodiment of the present application is provided.

[0025] Label explanation: data acquisition module 11, stress feature extraction module 12, sensor parameter construction module 13, stress change calculation module 14, coupling model correction module 15, model adaptive optimization module 16. DETAILED DESCRIPTION

[0026] Exemplary embodiments of the present disclosure are described below with reference to the accompanying drawings, which include various details of the embodiments of the present disclosure to assist in understanding, and should be considered as merely exemplary. Therefore, those skilled in the art should recognize that various changes and modifications can be made to the embodiments described herein without departing from the scope and spirit of the present disclosure. Also, descriptions of well-known functions and structures are omitted in the following description for the sake of clarity and conciseness.

[0027] Embodiment one, the quality flow meter measurement compensation method combined with pipeline stress perception provided by the embodiment of the present disclosure is described with reference to Figure 1 for illustration, the method comprises:

[0028] S1: A stress sensor array is arranged at a key position of the quality flow meter installation pipe section, real-time acquisition of stress distribution data in axial, radial and torsional directions of the pipeline is performed, and an original stress data set is obtained;

[0029] Specifically, the key position of the quality flow meter installation pipe section is determined, including the upstream connecting section of the pipeline, the flow meter sensing pipe section, and the downstream connecting section. The upstream connecting section is located within one to two pipe diameters in front of the flow meter, mainly bearing the fluid inlet pulsating pressure and installation stress; the flow meter sensing pipe section is the measurement core part of the installed quality flow meter body, directly reflecting the dynamic coupling state between the fluid and the sensing pipe; the downstream connecting section is located behind the flow meter outlet, used to capture the influence of fluid inertia and installation fixed constraints on overall stress transmission.

[0030] Stress sensors are installed at equal intervals along the axial direction of the pipeline at each key position. The interval of the measuring points is determined according to the length of the pipe section and the structural characteristics, and is preferably in the range of 1 / 20 to 1 / 50 of the total length of the pipeline to ensure the spatial resolution of the stress distribution. Each measuring point includes a three-dimensional direction measuring unit that can simultaneously collect stress signals in the axial, radial and torsional directions. The axial direction is along the center line of the pipeline, and the axial stress is used to represent the elongation or compression of the pipe section caused by fluid pressure, temperature expansion or installation expansion. The radial direction is perpendicular to the pipe wall, and the radial stress reflects the fluid internal pressure and the structural bulging effect. The torsional direction is the tangential direction around the pipe axis, and the torsional stress is used to describe the torsional response caused by asymmetric flow or support eccentricity. The stress sensor uses a strain gauge type stress sensor or a fiber Bragg grating sensor, wherein the former is suitable for low-temperature and static measurement environment, and the latter has the characteristics of strong anti-electromagnetic interference ability and high long-term stability. The sensor sampling frequency is set to 1000 Hz, which is triggered by a high-precision synchronous clock or a GPS time tag to ensure that the stress signals of each measuring point are time-aligned, and the time synchronization error is not more than 1 ms. The data acquisition module stores data with time stamp, measuring point number, direction component and stress value as index to form the original stress data set, which is defined as follows:

[0031] ;

[0032] wherein represents the original stress data set, i is the measuring point number, N represents the total number of measuring points, that is, the total number of stress sensors, j represents the stress direction component, x represents the axial direction, y represents the radial direction, z represents the torsional direction, and t is the sampling time, represents the stress value of the j direction of the i measuring point, and the data set provides a basis for subsequent stress characteristic parameter calculation.

[0033] S2: filtering and time series smoothing preprocessing of the original stress data set, extracting average stress, stress gradient, strain rate and vibration main frequency component, and constructing a stress characteristic parameter set;

[0034] Further, the step S2 further comprises:

[0035] The original stress data set is preprocessed, and the preprocessing includes filtering and time series smoothing;

[0036] Based on the preprocessed original stress data set, key characteristic parameters are extracted, including average stress, stress gradient, strain rate and vibration main frequency component, wherein the specific calculation formula of the average stress is:

[0037] ;

[0038] wherein This represents the average stress at measuring point i in the θ direction, where i represents the measuring point number and θ represents the stress direction. dt represents the instantaneous stress at measuring point i in the θ direction at time t, T is the length of the time window, and dt represents the time derivative.

[0039] The specific formula for calculating the stress gradient is as follows:

[0040] ;

[0041] in, This represents the stress gradient along the pipe axis at measuring point i, where i represents the measuring point and x represents the axial direction. and These represent the axial average stress at two adjacent measuring points. This represents the distance between two adjacent measuring points;

[0042] The specific formula for calculating strain rate is as follows:

[0043] ;

[0044] in, The value represents the strain rate at measuring point i in the θ direction, where i represents the measuring point number, θ represents the stress direction, and E represents the elastic modulus of the flowmeter sensor tube material. For the rate of change of stress, This represents the direction of measurement point i at time t. Instantaneous stress, Represents the degree of stress change. For time differentiation;

[0045] The specific formula for calculating the dominant frequency component of vibration is as follows:

[0046] ;

[0047] in, Let represent the dominant frequency component of the vibration at measurement point i, and argmax() denote the independent variable that makes the function reach its maximum value. The spectral amplitude represents the stress signal, where This represents the synthesized result of the three-directional stress signals at measurement point i;

[0048] A set of stress characteristic parameters is constructed by integrating all key stress characteristic parameters with the preprocessed original stress data.

[0049] Specifically, in order to eliminate noise and non-structural interference signals, the original stress data set is filtered and time series smoothed. The filtering process adopts a band-pass filtering method, in which the low-frequency cutoff frequency is set to 0.1 Hz, which is used to eliminate low-frequency trends caused by temperature drift or slow deformation of the base, and the high-frequency cutoff frequency is set to 50 Hz, which is used to suppress high-frequency noise and high-frequency vibration of the motor or the pipe wall. The filtered original stress data set is then time series smoothed by a sliding average algorithm, and the local average value is calculated with a window length of 5 to 10 sampling points, so as to effectively eliminate transient fluctuations and isolated pulses. The processed data form a smooth stress sequence, which has a more continuous waveform and can accurately reflect the actual trend of the stress change of the pipeline.

[0050] The original stress data set after preprocessing is obtained, and on this basis, key feature parameters reflecting the mechanical state of the pipeline are extracted, and a stress feature parameter set is constructed. The key feature parameters include average stress, stress gradient, strain rate and vibration main frequency component.

[0051] The average stress is used to describe the average force of each measuring point in the stable running stage, which is obtained by integrating the instantaneous stress value in the selected time window T and averaging, and the specific calculation formula is wherein, represents the average stress of measuring point i in the direction θ, i represents the measuring point number, and θ represents the stress direction, represents the instantaneous stress of measuring point i in the direction θ at time t, T is the length of the time window, which is determined by the actual sampling frequency and the analysis requirement, and 1s is selected here, and dt represents the time differential, which reflects the integration of the stress signal in continuous time.

[0052] The stress gradient is used to describe the stress change trend of the measuring point along the axial direction of the pipeline, which is calculated based on the average stress of adjacent measuring points in the axial direction, and the specific calculation formula is wherein, represents the stress gradient of measuring point i along the axial direction of the pipeline, i represents the measuring point, and x represents the axial direction, and respectively represent the axial average stress of two adjacent measuring points, represents the position interval of two adjacent measuring points, which is directly obtained by the previously set measuring point interval.

[0053] The strain rate is used to represent the dynamic response speed of the stress of each measuring point in three directions with time, which is calculated by the time derivative of the instantaneous stress, and the discrete form is wherein, represents the strain rate of measuring point i in the direction θ, i represents the measuring point number, θ represents the stress direction, and E represents the elastic modulus of the material of the flowmeter sensing pipe, which is provided by the manufacturer, is the stress change rate, representing the instantaneous stress of the measuring point i in the direction at the time t, representing the degree of stress change, is the time differential.

[0054] The vibration main frequency component is used to characterize the frequency with the most concentrated energy in the stress signal. The frequency spectrum is obtained by performing FFT fast Fourier transform on the stress signal of each measuring point in the selected time window, and the frequency corresponding to the maximum amplitude is extracted. The formula is represented as wherein, represents the vibration main frequency component of the measuring point i, and argmax() represents the independent variable that makes the function take the maximum value, represents the amplitude of the frequency spectrum of the stress signal, wherein represents the synthesis result of the three-direction stress signal at the measuring point i, and is defined as wherein x represents the axial direction, y represents the radial direction, and z represents the torsion direction.

[0055] Finally, the above four types of key feature parameters are uniformly organized according to the measuring point number, direction component and time stamp, and are integrated with the original stress data after preprocessing to obtain a stress feature parameter set.

[0056] S3: Obtain the geometric size and material performance of the mass flowmeter, and the spatial attitude information of the stress sensor, and construct a sensor parameter set;

[0057] Further, the step S3 further comprises:

[0058] According to the structural design drawing of the mass flowmeter, the geometric size of the sensing pipe is obtained in combination with the actual measurement on site, and the geometric size includes the sensing pipe wall thickness, the sensing pipe curvature radius, the fixed end length, and the effective measurement length;

[0059] The material performance of the sensing pipe is obtained through the product specification of the mass flowmeter, and the material performance includes the elastic modulus and the density;

[0060] The pipe segment coordinate system is established with reference to the pipeline structure, the installation coordinate system is established with reference to the stress sensor installation structure, the spatial mapping relationship between the pipe segment coordinate system and the installation coordinate system is established, the spatial coordinates, the direction angle and the coordinate conversion matrix of each stress sensor are recorded and associated, and the final spatial attitude information is obtained;

[0061] The geometric size, the material performance of the sensing pipe of the mass flowmeter and the spatial attitude information of the stress sensor are integrated to construct the sensor parameter set.

[0062] Specifically, to ensure that the stress characteristic parameters can accurately map the structural response characteristics of the mass flow meter, the geometric dimensions of the mass flow meter sensing pipe, material properties, and spatial attitude information of the stress sensor need to be systematically acquired. This process aims to establish a sensor parameter set, providing structural constraints and spatial geometry benchmarks for subsequent flow correction calculations.

[0063] First, for the acquisition of geometric dimension information of the sensing pipe, it mainly includes sensing pipe wall thickness, sensing pipe curvature radius, fixed end length, and effective measurement length. Accurate determination of geometric dimensions is the premise of stress distribution calculation and pipeline vibration characteristic modeling. The above geometric parameters can be obtained in two ways: one is to directly read the nominal size parameters of the sensing pipe by consulting the structural design drawings or product technical documents provided by the mass flow meter manufacturer; the second is to use a three-dimensional laser range finder, ultrasonic thickness gauge, or vernier caliper to measure the actual installation state of the pipe shape and obtain the geometric dimensions in the actual installation state.

[0064] Second, for the acquisition of material performance parameters of the sensing pipe, the focus includes elastic modulus and density. The elastic modulus is usually calibrated by the manufacturer in the product manual, and its accuracy has a direct impact on the flow correction calculation. The density parameter can be obtained by consulting the standard material performance manual.

[0065] Third, the acquisition of spatial attitude information of the stress sensor is crucial for establishing the coordinate correspondence between the external stress direction and the internal response of the flow meter. For this purpose, a unified pipe segment coordinate system and installation coordinate system need to be established at the stress sensor installation location. The pipe segment coordinate system takes the center of the flange face as the origin, defines the axial direction as the direction along the pipe centerline, the radial direction as the direction perpendicular to the pipe wall, and the torsion direction according to the right-hand rule, which is used to describe the external stress measurement direction. The installation coordinate system is established with the stress sensor factory reference surface or installation bracket plane as the reference, which is used to describe the spatial orientation of the sensor and internal measurement elements. The coordinate positions of each stress sensor in space can be obtained through three-dimensional laser scanning, total station measurement, or laser ranging equipment. At the same time, the pitch angle, yaw angle, and roll angle of the sensor in the installation state are measured using an electronic level, a three-axis attitude sensor, or a laser attitude measurement instrument, thereby establishing the spatial mapping relationship between the installation coordinate system and the pipe segment coordinate system. This mapping relationship is realized through a rotation transformation matrix calculated from the attitude angle parameters. After completing the above measurements, the spatial coordinates, direction angles, and coordinate conversion matrices of each measurement point are recorded and correlated to obtain the final spatial attitude information.

[0066] Finally, the geometric dimension information, material performance parameters and spatial attitude information are unified, sorted and standardized to form a complete sensor parameter set. The parameter set not only reflects the inherent structural characteristics of the sensor pipe of the mass flow meter, but also provides accurate physical basis for subsequent calculation of the influence of pipe stress changes on the vibration frequency and phase signal of the sensor pipe, thereby ensuring the scientificity and verifiability of the compensation model.

[0067] S4: Based on the stress characteristic parameter set and the sensor parameter set, the stress change amount of the pipe in the axial, radial and torsional directions is calculated to obtain a stress change data set;

[0068] Further, the step S4 further includes:

[0069] The stress characteristic parameter set is obtained, and the stress change amount of each measuring point in three directions is calculated by using the difference method from the static reference state. The specific calculation formula is:

[0070] ;

[0071] Wherein, represents the stress change amount of the measuring point i in the direction d at time t, i is the measuring point index, t is the time stamp, d is the direction index, including the axial direction x, the radial direction y and the torsional direction z, represents the instantaneous stress of the measuring point i in the direction d at time t, represents the average stress of the measuring point i in the direction d at the last sampling time, and Δt is the sampling time interval;

[0072] The weight coefficients of each measuring point are set based on the sensor parameter set;

[0073] The stress change amount of each measuring point in three directions is obtained, and the average stress change amount is obtained by combining the weight coefficients of each measuring point. The specific calculation formula is:

[0074] ;

[0075] Wherein, represents the average stress change amount of the pipe section in the direction d, and N represents the total number of measuring points, represents the stress change amount of the i-th measuring point in the direction d, represents the weight coefficient of the i-th measuring point, d is the direction index, and t represents the time stamp;

[0076] The stress change amount and the average stress change amount of all measuring points in three directions are integrated to form a stress change data set.

[0077] Further, the step S4 further includes:

[0078] The weight coefficients of each measuring point are set based on the sensor parameter set, and the specific setting formula is:

[0079] ;

[0080] wherein, represents the elastic modulus of the sensing tube material where the i-th measuring point is located, represents the effective measuring length of the sensing tube corresponding to the measuring point, represents the radial distance of the measuring point, represents the radial attenuation index, represents a small constant.

[0081] Specifically, after the construction of the sensor parameter set is completed, the stress change data set is formed by calculating the stress change amount of the pipeline in the axial, radial and torsional directions in combination with the stress characteristic parameter set obtained in the foregoing steps. The core purpose of this step is to convert the instantaneous stress characteristics at different measuring points into a global stress distribution expression that affects the structure of the flowmeter, providing dynamic input for the subsequent flow correction model.

[0082] First, to realize the spatial mapping and physical consistency of multi-measuring point stress data, the time series data collected by each stress sensor needs to be projected in a unified pipe segment coordinate system. Through the installation direction mapping relationship established in the previous step, the local stress components of each stress sensor are converted into stress components in the pipe segment coordinate system, so that the axial stress corresponds to the load distribution in the direction of the pipeline center line, the radial stress corresponds to the stress distribution in the normal direction of the pipe wall, and the torsional stress reflects the shear component in the rotation direction along the center axis. This process ensures that the stresses measured by each stress sensor are spatially superimposable and comparable, so that unified stress change analysis can be performed.

[0083] Secondly, for the stress change amount of each measuring point in three directions, the difference method with the static reference state is used to calculate, and the formula is:

[0084] ;

[0085] wherein, represents the stress change amount of measuring point i in direction d at time t, i is the measuring point index, t is the time stamp, d is the direction index, including the axial direction x, the radial direction y and the torsional direction z, represents the instantaneous stress of measuring point i in direction d at time t, represents the average stress of measuring point i in direction d at the last sampling time, representing the average stress in the initial stage of steady-state operation or the calibration stage, Δt is the sampling time interval, determined by the sensor sampling frequency, and the sampling frequency of the system is set to 1000 Hz, so Δt is 1 ms.

[0086] Again, after obtaining the stress variation of each measuring point in three directions, the stress distribution is weighted and averaged and spatially interpolated using the geometric size information and material performance parameters in the sensor parameter set, so as to derive the stress variation field of the whole pipe section. The weighting process comprehensively considers the distribution density of the measuring points in the spatial position, the radial distance of the measuring points to the center of the pipeline and the elastic modulus of the material and other factors, so that the calculation result can accurately reflect the overall deformation trend of the pipeline. The specific calculation formula is:

[0087]

[0088] represents the average stress variation of the pipe section in direction d, N represents the total number of measuring points, represents the stress variation of the i-th measuring point in direction d, which is derived from the calculation of the foregoing formula, represents the weight coefficient of the i-th measuring point, which is used to reflect the contribution of the measuring point to the stress field, d is the direction index, including the axial direction x, the radial direction y and the torsion direction z, and t represents the time stamp.

[0089] Specifically, the weight coefficient of the measuring point can be determined according to the radial position of the measuring point to the center of the pipeline, the spatial distribution density of the stress sensor and the elastic modulus of the material, and is used to reflect the difference in the contribution of each measuring point to the stress. The specific determination method is:

[0090]

[0091] represents the elastic modulus of the material of the i-th measuring point, which is derived from the material performance parameter acquisition in step S3, represents the effective measuring length of the measuring point corresponding to the sensor pipe, which is derived from the geometric size information in step S3, represents the radial distance of the measuring point, i.e. the radial distance of the measuring point to the center axis of the pipeline, which is derived from the three-dimensional spatial position measurement in step S3, represents the radial attenuation index, which is used to describe the attenuation of the stress contribution to the center axis, and can be set according to the structural characteristics of the pipeline and the experience of finite element analysis, generally taking 1~2, and here it is set to 1.5, represents a small constant to prevent division by zero error, which can be taken as 10 -6 , The physical meaning of is that the closer the measuring point is to the center of the pipeline, the more rigid the material is or the longer the measuring point pipe section is, the greater the contribution to the overall stress field is, so that the proportion in the weighted average is higher.

[0092] ​​​​Finally, all the stress change values of the measuring points in three directions and the corresponding time series are sorted out, and the average stress change values at each time are integrated to form a stress change data set. The data set not only contains the stress change values of multiple measuring points in three directions, but also retains the corresponding time stamp and spatial coordinate index information, forming a four-dimensional data structure that can reflect the dynamic load distribution evolution of the pipeline. The stress change data set is the basic data input for establishing the influence model of stress change on the vibration frequency and phase signal of the sensor pipe, and its physical meaning is to quantitatively describe the real structural stress dynamic response of the pipeline system caused by pressure pulsation, fluid inertia and mechanical constraint changes during operation, providing physical support for realizing accurate compensation of flow output.

[0093] S5: using the stress change data set and the sensor parameter set, establishing a coupling model of pipeline stress change on the vibration frequency and phase signal of the sensor pipe, and calculating the corrected mass flow value;

[0094] Further, the step S5 further comprises:

[0095] The stress change data set is obtained, and the average stress change values in three directions are converted into equivalent stress, and the specific calculation formula is:

[0096] ;

[0097] wherein, , and are the average stress change values in the axial direction x, the radial direction y and the torsional direction z respectively, represents the equivalent stress;

[0098] The stress characteristic parameter set is introduced, the coupling model is constructed, and the frequency offset and the phase offset are calculated, and the specific calculation formulas are respectively:

[0099] ;

[0100] ;

[0101] wherein, represents the frequency offset, represents the phase offset, is the frequency coupling coefficient, is the phase coupling coefficient, represents the equivalent stress, represents the frequency strain rate coupling coefficient, represents the phase strain rate coupling coefficient, represents the average strain rate, which is obtained by calculating the average value of the strain rates of each measuring point in three directions in the stress characteristic parameter set at time t, and T is a single vibration period, is the frequency stress gradient coupling coefficient, is the phase stress gradient coupling coefficient, represents the average stress gradient, which is obtained by calculating the average value of the stress gradient of each measuring point along the axial direction of the pipeline at time t, is the frequency vibration main frequency component coupling coefficient, is the phase vibration main frequency component coupling coefficient, is the average vibration main frequency component, which is obtained by calculating the average value of the vibration main frequency component of each measuring point at time t,

[0102] According to the sensor parameter set, the frequency structure offset and the phase structure offset are obtained, and the specific calculation formulas are as follows:

[0103] ;

[0104] ;

[0105] wherein, represents the frequency structure offset, the phase structure offset, and are the elastic modulus and the density under the calibration reference, and are the elastic modulus and the density under the current stress state, is the frequency structure sensitivity coefficient, is the phase structure sensitivity coefficient, is the structure coupling factor;

[0106] The frequency offset, the phase offset, the frequency structure offset and the phase structure offset are applied to the Coriolis flow formula to calculate the corrected mass flow value, and the specific calculation formula is as follows:

[0107] ;

[0108] wherein, is the corrected mass flow value, and are the flow meter original output frequency and phase signals, which are derived from the measurement data of the mass flow meter, represents the frequency offset, represents the phase offset, represents the frequency structure offset, represents the phase structure offset, is the instrument constant.

[0109] Specifically, after completing the multi-point stress measurement in the pipe segment and calculating the three-direction average stress variation along the axial, radial and torsional directions, this step aims to establish the physical coupling relationship between the stress disturbance and the mass flowmeter output signal, realize the dynamic compensation of frequency and phase, and thus obtain the corrected mass flow value. First, the average stress variation along the three directions obtained in step S4 is integrated into an equivalent stress, and the specific calculation formula is:

[0110]

[0111] are the average stress variations in the axial x, radial y and torsional z directions, respectively, with the unit of Pa, represents the equivalent stress, representing the overall stress disturbance intensity of the sensing pipe, and is the main driving factor of the frequency and phase shift of the flowmeter.

[0112] On this basis, the stress characteristic parameter set obtained before is introduced to construct the coupling model, and the frequency shift and phase shift are calculated, and the specific calculation formulas are:

[0113]

[0114]

[0115] represents the frequency shift, represents the phase shift, is the frequency coupling coefficient, is the phase coupling coefficient, both of which are obtained through calibration experiments and finite element simulation, represents the equivalent stress, represents the frequency strain rate coupling coefficient, represents the phase strain rate coupling coefficient, both of which are dimensionless, represents the average strain rate, which is obtained by averaging the strain rates of each measuring point in the stress characteristic parameter set at time t in three directions, with the unit of 1 / s, reflecting the response speed of the sensing pipe under dynamic load, and T is a single vibration period, used to normalize the strain rate to form a dimensionless term, is the frequency stress gradient coupling coefficient, is the phase stress gradient coupling coefficient, both of which have the unit of m / Pa, represents the average stress gradient, which is obtained by averaging the stress gradients of each measuring point along the pipeline axial direction at time t, is the frequency vibration main frequency component coupling coefficient, is the phase vibration main frequency component coupling coefficient, both of which have the unit of S, ​​​​​​​The average vibration main frequency component is obtained by calculating the average value of the vibration main frequency component of each measuring point at time t.

[0116] Specifically, in the above frequency and phase offset calculation model, each coupling coefficient introduced is used to represent the quantitative coupling relationship between stress disturbance, structural characteristics and flow signal response. To ensure the traceability and engineering realizability of the model parameters, the value range and calculation basis of the coefficients 、 、 are determined by combining theoretical analysis, finite element simulation and experimental calibration.

[0117] Firstly, the frequency coupling coefficient and the phase coupling coefficient are obtained by finite element simulation and static and dynamic combined calibration experiment. Specifically, in the simulation stage, according to the sensor parameter set obtained in step S3, a three-dimensional finite element model is established, and the fundamental frequency and phase response curve are calculated under different stress levels. By linear regression fitting the relationship between the frequency change rate and the equivalent stress, the preliminary estimated value of is obtained. Similarly, by fitting the phase drift and stress level, the is obtained. Then in the experimental calibration, the known axial load and torsional load are applied to the mass flowmeter, and the output signal frequency and phase are monitored to correct the simulation results. In practical engineering applications, usually .

[0118] Secondly, the frequency strain rate coupling coefficient and the phase strain rate coupling coefficient respectively reflect the dynamic influence of strain rate change on frequency response and phase response. Although they are similar in form, there are differences in mechanism. It mainly reflects the transient modulation effect of load change rate on stress stiffness, which will cause a slight drift of frequency. It mainly affects the signal propagation time delay and phase integral error, and is more sensitive to phase response. Therefore, in theory, they should not be completely combined, although a unified strain rate coupling coefficient can be used under the condition of approximate linearity, but in order to maintain the model interpretability and adapt to the nonlinear compensation stage, they are still separated here. The calibration method uses dynamic loading experiment: record the instantaneous response change of frequency and phase under different strain rates, and fit the values of the two respectively. The typical value range is between 0.02 and 0.06, between 0.03 and 0.08.

[0119] Thirdly, the frequency stress gradient coupling coefficient Phase stress gradient coupling coefficient Reflects the coupling effect of the non-uniformity of the spatial stress distribution on the vibration mode of the sensing tube. The coefficient is obtained by combining finite element modal analysis and measured modal calibration. Specifically, the modal frequency variation of the sensing tube under a unit stress gradient (1 Pa / m) is calculated by simulation. Then, using modal experiments, a controllable distributed load is applied to the pipeline, and the frequency and phase changes are measured to correct the simulation results. In practical applications, Typical values .

[0120] Finally, the frequency vibration main frequency component coupling coefficient , the phase vibration main frequency component coupling coefficient is used to represent the direct coupling effect of the vibration main frequency component on the frequency and phase shift, and its unit is seconds, which ensures and are dimensionless correction quantities. The coefficient is obtained by combining finite element modal analysis and vibration experiments: record the main frequency component change and its corresponding frequency and phase response under different excitation frequencies and amplitudes, and determine and numerical values by regression fitting combined with simulation results. Typical values range from . Both reflect the modulation effect of the main modal frequency change on the signal response, and are used in the model to modify the stress and strain rate coupling term to achieve comprehensive dynamic compensation of the frequency and phase shift.

[0121] Through the above calibration process, a set of coupling parameter library matched with a specific type of flowmeter can be formed, and in field applications, the corresponding parameter group is automatically called according to the flowmeter structure type, sensing tube material and installation form, thereby ensuring the accuracy and reproducibility of the compensation model.

[0122] In addition, to consider the correction effect of the sensing tube's geometric size and material properties on the frequency and phase, a structure correction term is introduced, and the specific calculation formulas are as follows:

[0123] ;

[0124] ;

[0125] wherein represents the frequency structure offset, the phase structure offset, and are the elastic modulus and density under the calibration reference, and are the elastic modulus and density under the current stress state, is the frequency structure sensitivity coefficient, is the phase structure sensitivity coefficient, is the structural coupling factor, determined by the wall thickness, curvature radius and fixed end length of the sensing tube.

[0126] Specifically, and respectively represent the elastic modulus and density of the sensing tube material under the factory calibration reference condition, usually provided by the manufacturer or obtained from the standard material database, and are stored in the sensor parameter set. and are the material parameters under the current operating state, obtained by temperature and stress dependent correction models, and the specific calculation formulas are respectively:

[0127] ;

[0128] ;

[0129] wherein, is the stress softening coefficient, reflecting the tendency of the elastic modulus of the material to decrease under high stress, is the thermal expansion density coefficient, describing the influence of temperature change on the density of the material, and the typical experimental value range of the two is , , is the current equivalent stress, Tem is the current room temperature, which can be directly monitored by a temperature sensor, is the stress under the calibration reference, indicating that the sensing tube is in a steady mechanical state under no flow or calibration flow when the flowmeter is calibrated by standard medium after leaving the factory, which corresponds to the static stress distribution of the instrument under zero flow or standard flow conditions. In general, it is in the range of tens to hundreds of kPa, is the temperature under the calibration reference, indicating the reference temperature when the instrument is calibrated or the system is initialized, which is generally 293.15K.

[0130] Parameter and are the frequency structural sensitivity coefficient and the phase structural sensitivity coefficient, respectively, used to characterize the sensitivity of the natural frequency and phase response of the sensing tube to changes in material properties. The determination method is as follows: in the finite element modal analysis, a small perturbation is applied to the elastic modulus and density of the sensing tube model, and the relative change rate of the natural frequency and phase response is calculated. Then in the experimental stage, the frequency and phase response are measured and corrected by adjusting the temperature and loading conditions. Through the verification of multiple flowmeters, the typical range is 0.8s -1 to 1.5s -1 , the typical range is 0.5 rad to 1.2 rad.

[0131] Parameter is the structural coupling factor, which reflects the comprehensive influence of the sensing tube's geometric characteristics on the sensitivity, and its expression is where h is the sensing tube wall thickness, is the sensing tube curvature radius, is the fixed end length, and all of the above data are obtained from the sensor parameter set. The formula shows that the greater the curvature and the shorter the fixed end, the greater the overall stiffness and the weaker the structural coupling. According to the range of conventional flowmeter structural parameters, a typical value of is 0.05 to 0.25.

[0132] Through the joint determination of the above parameters, a complete structural coupling chain from geometric size and material performance to frequency and phase correction is established, so that the structural correction terms and can accurately reflect the influence of the structural characteristics of flowmeters of different types and in different installation states, thereby ensuring the engineering consistency and physical interpretability of the calculation results of the coupling model.

[0133] Finally, the frequency offset, phase offset, frequency structural offset, and phase structural offset are applied to the Coriolis flow formula to calculate the corrected mass flow value:

[0134]

[0135] where is the corrected mass flow value, in kg / s, and are the original output frequency and phase signals of the flowmeter, which are derived from the measurement data of the mass flowmeter, represents the frequency offset, represents the phase offset, represents the frequency structural offset, represents the phase structural offset, is the instrument constant, which is used to map the frequency and phase signals into the mass flow value, and its typical value range is .

[0136] Through the above model, step S5 realizes the comprehensive mapping of the three-direction average stress change obtained in step S4, the stress characteristic parameters extracted in step S2, and the sensor parameter set constructed in step S3 to the flowmeter output signal, realizes the dynamic compensation of the frequency and phase, and thus obtains a high-precision and reliable corrected mass flow output.

[0137] S6: According to the comparison results of the corrected mass flow value and the calibration point during long-term operation, a parameter optimization algorithm is applied to adaptively adjust the coupling coefficients, and the coupling model is optimized.

[0138] ​Specifically, after the real-time flow correction calculation is completed and the corrected mass flow rate value is generated in step S5, to ensure the accuracy and robustness of the coupled model in long-term operation, the coupled model parameters are adaptively optimized and self-corrected using the calibration point data obtained periodically. The calibration point data are reference flow rate values measured at several calibration times using a high-precision comparison flowmeter or a standard volume sampling device. To facilitate parameter identification and engineering implementation, the coupled coefficient set is defined as .

[0139] The optimization process can be implemented by an iterative algorithm, such as gradient descent or least squares fitting. The optimization objective function is defined as:

[0140] ;

[0141] where N is the total number of calibration point data, is the standard flow rate value corresponding to the kth calibration point data, is the corrected mass flow rate value calculated at the corresponding time point . By minimizing the objective function J, the coupled coefficients are gradually adjusted to make the flow correction model better reflect the physical relationship between the pipeline stress and the flow signal.

[0142] In the optimization process, boundary constraints and historical parameter records can be combined to ensure model stability and physical interpretability. After iterative optimization, the new self-corrected coupled coefficient set is output, which can be applied in real time in subsequent flow correction, so that the mass flowmeter maintains high-precision measurement in long-term operation and stress disturbance, and forms a traceable historical compensation parameter record, providing a basis for subsequent operation analysis and further model optimization.

[0143] Embodiment Two, based on the same inventive concept as the mass flowmeter measurement compensation method combined with pipeline stress sensing in the aforementioned embodiments, the present application also provides a mass flowmeter measurement compensation system combined with pipeline stress sensing. Please refer to FIG. 2, Figure 2 , which comprises:

[0144] A data acquisition module 11 is used to arrange an array of stress sensors at key positions of the mass flowmeter installation pipe section, to acquire real-time stress distribution data in the axial, radial and torsional directions of the pipeline, and to obtain an original stress data set;

[0145] A stress feature extraction module 12 is used to perform filtering and time series smoothing preprocessing on the original stress data set, to extract the average stress, stress gradient, strain rate and vibration main frequency components, and to construct a stress feature parameter set;

[0146] A sensor parameter construction module 13 is configured to acquire geometric dimensions and material properties of the mass flow meter and spatial posture information of the stress sensor, and construct a sensor parameter set;

[0147] A stress change calculation module 14 is configured to calculate stress change amounts of the pipeline in axial, radial and torsional directions based on the stress characteristic parameter set and the sensor parameter set, and obtain a stress change data set;

[0148] A coupling model correction module 15 is configured to establish a coupling model of stress change of the pipeline to vibration frequency and phase signals of the sensor pipe by using the stress change data set and the sensor parameter set, and calculate a corrected mass flow value;

[0149] A model adaptive optimization module 16 is configured to apply a parameter optimization algorithm to adaptively adjust the coupling coefficient according to comparison results of the corrected mass flow value and the calibration point in a long-term operation process, and realize optimization of the coupling model.

[0150] The technical features of the above embodiments can be combined in any manner. To make the description concise, all possible combinations of the technical features in the above embodiments are not described, but as long as the combinations of the technical features do not exist contradictions, they should be considered as the scope of the present disclosure.

[0151] The above description of the disclosed embodiments enables a person skilled in the art to implement or use the present application. Various modifications to the embodiments will be apparent to those skilled in the art, and the general principles defined herein can be implemented in other embodiments without departing from the spirit or scope of the application. Therefore, the present application will not be limited to the embodiments shown herein, but will conform to the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A method of measurement compensation for a pipe stress-aware mass flowmeter, characterized by, The method comprises: A stress sensor array is arranged at a key position of a mass flow meter installation pipe section, stress distribution data in axial, radial and torsional directions of the pipeline are collected in real time, and an original stress data set is obtained; The original stress data set is preprocessed by filtering and time series smoothing, average stress, stress gradient, strain rate and vibration main frequency components are extracted, and a stress characteristic parameter set is constructed; The geometric size and material performance of the mass flow meter are obtained, and the spatial attitude information of the stress sensor is obtained, and a sensor parameter set is constructed; Based on the stress characteristic parameter set and the sensor parameter set, the stress change amount of the pipeline in the axial, radial and torsional directions is calculated, and a stress change data set is obtained; The stress change data set and the sensor parameter set are used to establish a coupling model of the pipeline stress change to the vibration frequency and phase signal of the sensing pipe, and a corrected mass flow value is calculated, wherein the corrected mass flow value further comprises: The stress change data set is obtained, and the average stress change amount in three directions is converted into equivalent stress, and the specific calculation formula is: ； wherein , and are the average stress variations in axial x, radial y and torsional z directions, respectively, represents the equivalent stress; The stress characteristic parameter set is introduced to construct the coupling model, and the frequency offset and the phase offset are calculated, and the specific calculation formulas are respectively: ； ； wherein, represents a frequency offset, represents a phase offset, is a frequency coupling coefficient, is a phase coupling coefficient, represents an equivalent stress, represents a frequency strain rate coupling coefficient, represents a phase strain rate coupling coefficient, represents an average strain rate, obtained by taking the average of strain rates in three directions of each measuring point at time t, and T is a single vibration period, is a frequency stress gradient coupling coefficient, is a phase stress gradient coupling coefficient, represents an average stress gradient, obtained by taking the average of stress gradients in the axial direction of each measuring point at time t, is a frequency vibration main frequency component coupling coefficient, is a phase vibration main frequency component coupling coefficient, is an average vibration main frequency component, obtained by taking the average of vibration main frequency components of each measuring point at time t; The frequency structure offset and the phase structure offset are obtained according to the sensor parameter set, and the specific calculation formulas are respectively: ； ； wherein, represents a frequency structure offset, represents a phase structure offset, and E0and p0are the elastic modulus and density, respectively, at a calibration reference, and Eand p are the elastic modulus and density, respectively, at a current stress state, is a frequency structure sensitivity coefficient, is a phase structure sensitivity coefficient, is a structure coupling factor; The frequency offset, the phase offset, the frequency structure offset and the phase structure offset are applied to the Coriolis flow formula to calculate the corrected mass flow value, and the specific calculation formula is: ； wherein, is the corrected mass flow value, and are the flow meter raw output frequency and phase signals, derived from the measurement data of the mass flow meter, represents the frequency offset, represents the phase offset, represents the frequency structure offset, represents the phase structure offset, is the meter constant; According to the comparison result of the corrected mass flow value and the calibration point in the long-term operation process, a parameter optimization algorithm is applied to adaptively adjust the coupling coefficient, and the coupling model is optimized.

2. The method of claim 1, wherein, The stress characteristic parameter set is constructed, comprising: The original stress data set is preprocessed, and the preprocessing comprises filtering and time series smoothing; Based on the preprocessed original stress data set, key characteristic parameters are extracted, the key characteristic parameters comprise average stress, stress gradient, strain rate and vibration main frequency component, wherein the specific calculation formula of the average stress is: ； wherein, represents the average stress of measuring point i in the direction of θ, i represents the serial number of the measuring point, and θ represents the stress direction, represents the instantaneous stress of measuring point i in the direction of θ at time t, T is the length of the time window, and dt represents the time differential; The specific calculation formula of the stress gradient is: ； wherein, represents the axial stress gradient of the measuring point i along the pipe, i represents the measuring point, and x represents the axial direction, respectively represent the axial average stresses of two adjacent measuring points, represents the position distance between the two adjacent measuring points; The specific calculation formula of the strain rate is: ； wherein, represents the strain rate of the measuring point i in the direction of θ, i represents the serial number of the measuring point, θ represents the stress direction, and E represents the elastic modulus of the material of the flowmeter sensor pipe, is the stress change rate, represents the instantaneous stress of the measuring point i in the direction of θ at the moment t, represents the stress change degree, is the time differential;​ The specific calculation formula of the vibration main frequency component is: ； wherein, represents the dominant frequency component of the vibration of the measuring point i, argmax() represents taking the argument of the function that takes the maximum value, represents the spectral amplitude of the stress signal, wherein represents the synthesis result of the three-direction stress signals at the measuring point i; All key stress characteristic parameters and preprocessed original stress data are integrated to construct the stress characteristic parameter set.

3. The method of claim 1, wherein, The sensor parameter set is constructed, comprising: According to the mass flow meter structure design drawing, the geometric size of the sensing pipe is obtained in combination with the actual measurement, and the geometric size comprises the sensing pipe wall thickness, the sensing pipe curvature radius, the fixed end length and the effective measurement length; The material performance of the sensing pipe is obtained through the mass flow meter product specification, and the material performance comprises the elastic modulus and the density; The spatial mapping relationship between the pipe section coordinate system and the installation coordinate system is established by taking the pipeline structure as a reference to establish the pipe section coordinate system and taking the stress sensor installation structure as a reference to establish the installation coordinate system, the spatial coordinates, the direction angle and the coordinate conversion matrix of each stress sensor are recorded and associated, and the final spatial attitude information is obtained; The geometric size, the material performance of the mass flow meter sensing pipe and the spatial attitude information of the stress sensor are integrated to construct the sensor parameter set.

4. The method of claim 1, wherein, The stress change data set is obtained, comprising: The stress characteristic parameter set is acquired, and the stress change amount of each measuring point in three directions is calculated by using the difference value method of the static reference state, and the specific calculation formula is: ； wherein, represents the stress change amount of the measurement point i in the direction d at the time t, i is the measurement point index, t is the time stamp, d is the direction index, including the axial direction x, the radial direction y and the torsional direction z, represents the instantaneous stress of the measurement point i in the direction d at the time t, represents the average stress of the measurement point i in the direction d at the last sampling time, Δt is the sampling time interval; The weight coefficients of each measuring point are set based on the sensor parameter set; The stress change amount of each measuring point in three directions is acquired, and the average stress change amount is acquired by combining the weight coefficients of each measuring point, and the specific calculation formula is: ； wherein, represents the average stress change in direction d of the pipe section, N represents the total number of measurement points, represents the stress change in direction d of the i-th measurement point, represents the weight coefficient of the i-th measurement point, d is the direction index, and t represents the time stamp; The stress change data set is formed by integrating the stress change amount of all measuring points in three directions and the average stress change amount.

5. The method of claim 4, wherein, The weight coefficients of each measuring point are set based on the sensor parameter set, including: The weight coefficients of each measuring point are set based on the sensor parameter set, and the specific setting formula is: ； wherein, represents the elastic modulus of the sensor tube material at the i-th measuring point, represents the effective measuring length of the sensor tube corresponding to the measuring point, represents the radial distance of the measuring point, represents the radial attenuation exponent, represents a small constant.

6. A mass flow meter measurement compensation system incorporating pipe stress awareness, characterized by, The system is used to implement the quality flow meter measurement compensation method combined with pipeline stress sensing in any one of claims 1 to 5, and the system comprises: A data acquisition module, which is used to arrange a stress sensor array at a key position of a quality flow meter installation pipe section, to acquire stress distribution data in the axial, radial and torsional directions of the pipeline in real time, and to obtain an original stress data set; A stress feature extraction module, which is used to perform filtering and time sequence smoothing preprocessing on the original stress data set, to extract average stress, stress gradient, strain rate and vibration main frequency components, and to construct a stress characteristic parameter set; A sensor parameter construction module, which is used to acquire the geometric size and material performance of the quality flow meter and the spatial posture information of the stress sensor, and to construct a sensor parameter set; A stress change calculation module, which is used to calculate the stress change amount of the pipeline in the axial, radial and torsional directions based on the stress characteristic parameter set and the sensor parameter set, and to obtain a stress change data set; A coupling model correction module, which is used to establish a coupling model of the pipeline stress change on the vibration frequency and phase signal of the sensing pipe by using the stress change data set and the sensor parameter set, and to calculate a corrected quality flow value; A model adaptive optimization module, which is used to adaptively adjust the coupling coefficient by applying a parameter optimization algorithm according to the comparison result of the corrected quality flow value and the calibration point in the long-term operation process, to realize the optimization of the coupling model.

Images