CN-122016017-A - Method for conveniently improving natural gas metering accuracy

CN122016017ACN 122016017 ACN122016017 ACN 122016017ACN-122016017-A

Abstract

The invention relates to the technical field of natural gas metering, in particular to a method convenient for improving the accuracy of natural gas metering, which comprises the steps of obtaining original metering data in the same metering period and forming an original metering sequence, extracting working condition disturbance characteristics of characteristic pressure regulating pulsation and flow fluctuation to obtain a working condition disturbance vector, extracting component related information and determining component agent quantity, mapping to obtain a component characterization result, determining a current working condition based on the component characterization result and the working condition disturbance vector, correcting the original metering sequence according to the current working condition to obtain a first corrected metering result, further screening anchoring segments, carrying out bounded updating on the corrected rule, and correcting subsequent original metering data to obtain a second corrected metering result, thereby improving the continuity, stability and accuracy of the corrected metering result.

Inventors

ZHANG CHUNYAN
WANG YU
Zang Luyang
ZHANG WEISHEN
LAI RONGBIN
JI GANG

Assignees

江西省天然气管道有限公司运营分公司

Dates

Publication Date: 20260512
Application Date: 20260416

Claims (10)

1. A method for facilitating improved accuracy in natural gas metering, comprising: acquiring original metering data of an acquisition unit in the same metering period to form an original metering sequence; Extracting working condition disturbance characteristics of the characteristic pressure regulating pulsation and the flow fluctuation based on the original metering sequence to obtain a working condition disturbance vector; extracting component related information based on the original metering sequence, determining component proxy quantity based on the component related information, and carrying out mapping treatment on the component proxy quantity to obtain a corresponding component characterization result; Determining the current working condition based on the component characterization result and the working condition disturbance vector together, and correcting the original metering sequence according to a correction rule corresponding to the current working condition to obtain a first corrected metering result; screening the anchoring segments according to screening conditions based on the original metering sequence, the current working condition, the component representation result, the working condition disturbance vector and the first correction metering result; and carrying out bounded updating on the correction rule based on the anchoring segment, and correcting the follow-up original metering data according to the updated correction rule to obtain a second correction metering result.
2. The method for facilitating the improvement of the accuracy of natural gas metering according to claim 1, wherein the specific step of obtaining the raw metering data of the acquisition unit in the same metering period and forming the raw metering sequence comprises the following steps: acquiring original metering data, and respectively carrying out buffering and time sequence arrangement according to data types, wherein the original metering data comprises acquisition records and pressure fluctuation data; determining sequence time according to the acquisition time corresponding to the acquisition record; Extracting pressure fluctuation segments in time intervals between adjacent sequence moments based on the pressure fluctuation data; And merging the acquired record and the corresponding pressure fluctuation segment to generate an original metering sequence.
3. The method for facilitating improvement of natural gas metering accuracy according to claim 1, wherein the specific step of extracting the operating condition disturbance features of the characteristic pressure regulating pulsation and the flow fluctuation based on the original metering sequence to obtain the operating condition disturbance vector comprises the following steps: performing time sequence segmentation processing on the original metering sequence to obtain a plurality of time windows; extracting working condition disturbance characteristics in each time window based on the acquired record and the pressure fluctuation data; Carrying out unified dimension processing on the working condition disturbance characteristics, and merging the results after the unified dimension processing into working condition disturbance candidate vectors corresponding to the time window; Performing difference comparison on the working condition disturbance candidate vector of the current time window and the working condition disturbance vector output by the previous time window, and determining the working condition disturbance vector of the current time window based on a comparison result; The calculation formula of the working condition disturbance candidate vector is as follows: , , , ; ; Wherein, the Representing the working condition disturbance candidate vector corresponding to the nth sliding window, Respectively representing the pressure difference fluctuation amplitude, the main fluctuation frequency, the variation coefficient of the instantaneous volume flow and the vector component of the sound velocity after the drift rate of the nth sliding window is subjected to unified dimension processing, For the magnitude of the pressure difference fluctuation, Representing the dominant pulsing frequency within the nth sliding window, A coefficient of variation representing the instantaneous volumetric flow rate within the nth sliding window, The drift rate of sound velocity in the nth sliding window relative to the previous window is represented, and the drift rate of sound velocity in the nth sliding window relative to the previous window is represented by 2.0 kilopascals, 1.5 hertz, 0.20 and 0.05, which respectively represent the fluctuation amplitude of the pressure difference, the main pulsation frequency, the variation coefficient of the instantaneous volume flow and the normalized upper limit of the sound velocity relative to the drift rate of the previous sliding window.
4. The method of claim 3, wherein the acquisition record includes an original cis-trans propagation time difference parameter, an original sound velocity parameter, an instantaneous volume flow parameter, an original pressure parameter, and an original temperature parameter, the original sound velocity parameter, the original pressure parameter, and the original temperature parameter being used to characterize the composition-related information.
5. The method for facilitating the improvement of the metering accuracy of the natural gas according to claim 4, wherein the specific steps of extracting the component related information based on the original metering sequence, determining the component proxy amount based on the component related information, and performing mapping processing on the component proxy amount to obtain the corresponding component characterization result include: a standard gas sample library is constructed in advance, and the standard gas sample library comprises reference samples corresponding to the gas components and the metering data under different working conditions; Based on the original sound velocity parameter, the original pressure parameter and the original temperature parameter in the current time window, a corresponding reference sample is called in the standard gas sample library, interpolation processing is carried out, and an intermediate component value corresponding to the current time window is obtained; merging the intermediate component values to obtain component proxy quantities corresponding to the current time window; dividing the component proxy quantity into intervals to obtain corresponding component interval numbers, and taking the component interval numbers as component characterization results; the calculation formula of the intermediate component value is as follows: ; In the formula, Representing the value of the intermediate component at the pressure-temperature combination point corresponding to the a-th pressure level and the b-th temperature level, Represents the average value of sound velocity, a represents the pressure level number, b represents the temperature level number, Represents the average value of the speed of sound within the current sliding window, Represents the standard sound velocity value for the lower hydrogen volume fraction level samples at the a-th pressure level and the b-th temperature level, Representing the standard sound velocity values for the higher hydrogen volume fraction level samples at the a-th pressure level and the b-th temperature level, Representing the hydrogen gas volume fraction values for the lower hydrogen gas volume fraction level samples at the a-th pressure level and the b-th temperature level, Representing hydrogen volume fraction values for the higher hydrogen volume fraction level samples at the a-th pressure level and the b-th temperature level; The calculation formula of the component agent amount is as follows: ; wherein G is the component proxy corresponding to the current sliding window, 、、、 Respectively represent the intermediate component values corresponding to the 4 pressure temperature combination points, Representing the position coefficient of the current sliding window pressure average between two adjacent pressure levels, Representing the coefficient of position of the current sliding window temperature average between two adjacent temperature levels.
6. The method for facilitating the improvement of the accuracy of natural gas metering according to claim 4, wherein the specific step of determining the current working condition based on the component characterization result and the working condition disturbance vector together comprises: Pre-constructing a working condition mapping table, wherein the working condition mapping table is used for storing the corresponding relation among the group interval number, the disturbance grade number, the flow subsection number and the current working condition; determining a current component interval number based on the component characterization result, determining a current disturbance level based on the working condition disturbance vector, and determining a current flow section based on instantaneous volume flow parameters in a current time window; And carrying out joint indexing based on the current component interval number, the current disturbance grade number and the current flow segmentation number, and searching and determining the current working condition from the working condition mapping table.
7. The method for facilitating the improvement of the accuracy of natural gas metering according to claim 6, wherein the specific step of correcting the original metering sequence according to the correction rule corresponding to the current working condition to obtain the first corrected metering result comprises the following steps: pre-constructing a platform calibration matrix, wherein the platform calibration matrix is used for storing correction rules corresponding to all current working conditions; based on the current working condition, a corresponding correction rule is called in a platform calibration matrix, and the correction rule is used as a correction basis of the instantaneous volume flow in the current time window; correcting the instantaneous volume flow in the current time window based on the correction rule to obtain an instantaneous volume flow correction value; determining an accumulated gas quantity corresponding to a current time window based on the instantaneous volume flow correction value, and generating a first correction metering result according to the instantaneous volume flow correction value and the accumulated gas quantity; Wherein, the calculation formula of the instantaneous volume flow correction value and the accumulated gas volume is as follows: , , In the formula, Representing the instantaneous volumetric flow correction value corresponding to the ith sample point in the nth window, For the original instantaneous volume flow corresponding to the ith sample point in the nth window, Representing the number of the current group interval corresponding to the nth window Current disturbance level numbering Current flow segment numbering The only corresponding metering correction parameter is used for metering the metering correction parameter, Represents the window accumulated gas quantity corresponding to the nth window, Representing the number of samples in the nth window, Representing the time interval between adjacent sample points.
8. The method for facilitating improved accuracy of natural gas metering of claim 7, further comprising: Executing the working condition boundary judgment on the current time window; And when the judgment result is that the current time window is at the working condition switching boundary, executing consistency processing on the correction rule of the current time window.
9. A method for facilitating improved accuracy of natural gas metering according to claim 3, wherein the specific step of screening the anchoring fragments according to screening conditions comprises: screening the continuous time window group to determine candidate anchor windows meeting the continuous consistency condition and the stability condition of the current working condition; And merging candidate anchor windows which are continuously adjacent in time and have no screening condition in the middle to obtain the anchor fragment.
10. The method for facilitating the improvement of the accuracy of natural gas metering according to claim 9, wherein the specific steps of performing a bounded update on the correction rule based on the anchor segment and correcting the subsequent original metering data according to the updated correction rule to obtain the second corrected metering result include: Determining adjacent component intervals based on the anchoring fragments, and calling correction rules corresponding to the adjacent component intervals under the same disturbance level and the same flow segmentation based on a platform calibration matrix to determine candidate correction rules corresponding to the anchoring fragments; performing clipping processing on the candidate correction rule and the correction rule before updating to obtain a correction rule; Determining an update uncertainty based on the correction rule, and determining whether the correction rule is effective according to the update uncertainty; When the correction and correction rule is effective, correcting the instantaneous volume flow in the follow-up original metering data based on the correction and correction rule to obtain a second correction and measurement result; wherein, the calculation formula of uncertainty is: ; In the formula, Representing the update uncertainty corresponding to the a-th anchor segment, Indicating the corresponding reference uncertainty of the current working condition in the platform calibration matrix, Representing the uncertainty corresponding to the dispersion of the instantaneous volumetric flow correction values within the anchor segment, Indicating the uncertainty corresponding to the single update step.

Description

Method for conveniently improving natural gas metering accuracy Technical Field The invention relates to the field of natural gas metering, in particular to a method for conveniently improving the accuracy of natural gas metering. Background In a dynamic settlement metering scene that the hydrogen-doped natural gas is conveyed to industrial and commercial users after pressure regulation, particularly in an ultrasonic metering branch at the front end of the industrial and commercial users, gas components, sound velocity, density, viscosity and flow stability can synchronously change along with the operation process, and meanwhile, the user side also has working conditions such as small flow fluctuation, frequent start and stop, short-time pulsation superposition and the like, so that deviation is easy to occur between an original metering value and an actual channel gas quantity of the metering instrument along with time change, and short-time deviation accumulation can be further formed. Therefore, how to perform online identification, constraint correction and bounded updating on the deviation without depending on continuous online chromatography, and stably improve the metering accuracy of natural gas becomes a technical problem to be solved under the scene. At present, the prior art usually uses the original metering model to settle accounts on the basis that the metering instrument has completed verification and recheck according to the period, and when the gas condition changes obviously, the metering reliability is restored by the way of supplementary detection, recheck or recalibration, that is, the existing scheme is mainly based on the existing calibration and post-correction processing path, and more depends on the existing model and the post-correction to maintain the metering result. However, in the dynamic settlement metering scene after the pressure regulation of the hydrogen-doped natural gas, as the existing correction logic is used for separately processing the component change and the flow pulsation, and a mechanism for uniformly modeling the component change and the flow pulsation by taking the component change and the flow pulsation as the chains with the same error source is lacking, the problem that the correction directions of the same metering instrument are inconsistent in different time periods is easy to occur in the short-time fluctuation scene after the pressure regulation, and the continuous, stable and consistent correction of the metering deviation under the dynamic working condition is difficult. Disclosure of Invention The invention aims to solve the problems that in a dynamic settlement metering scene after the pressure regulation of hydrogen-doped natural gas in the prior art, on one hand, the prior art separately processes component change and flow pulsation to cause inconsistent correction directions of the same metering instrument in different time periods, and on the other hand, correction rules lack of stable reference and update constraint, are easy to be frequently regulated due to short-time disturbance influence, and cause insufficient continuity and reliability of metering results, and provides a method convenient for improving the metering accuracy of the natural gas. In order to achieve the purpose, the technical scheme adopted by the invention is that the method for conveniently improving the metering accuracy of the natural gas comprises the following steps: acquiring original metering data of an acquisition unit in the same metering period to form an original metering sequence; Extracting working condition disturbance characteristics of the characteristic pressure regulating pulsation and the flow fluctuation based on the original metering sequence to obtain a working condition disturbance vector; extracting component related information based on the original metering sequence, determining component proxy quantity based on the component related information, and carrying out mapping treatment on the component proxy quantity to obtain a corresponding component characterization result; Determining the current working condition based on the component characterization result and the working condition disturbance vector together, and correcting the original metering sequence according to a correction rule corresponding to the current working condition to obtain a first corrected metering result; screening the anchoring segments according to screening conditions based on the original metering sequence, the current working condition, the component representation result, the working condition disturbance vector and the first correction metering result; and carrying out bounded updating on the correction rule based on the anchoring segment, and correcting the follow-up original metering data according to the updated correction rule to obtain a second correction metering result. Preferably, the specific step of obtaining the original metering data of the acquisition u