Gas-liquid two-phase metering method, system and equipment based on fusion algorithm

Abstract

The invention belongs to the field of oil-gas exploration, particularly relates to a gas-liquid two-phase metering method, system and equipment based on a fusion algorithm, and aims to solve the problem that gas-liquid two phases are difficult to meter. The method comprises the following steps: collecting original data in a target pipeline as input data through a multi-sensor module; converting the input data through an A / D converter to obtain digital input data; filtering and denoising the digital input data to obtain preprocessed input data; performing feature extraction on the preprocessed input data to obtain a sensor measurement value; performing fusion calculation on the sensor measurement values through a multi-source data fusion algorithm to obtain an initial flow value; based on a sensor measurement value, an initial flow value is corrected through a gas-liquid two-phase proportion measured by a precession vortex sensor and a high-precision matrix data iterative transmission algorithm between differential pressure flowmeters, and a corrected two-phase flow value is obtained.

Description

Background Art

[0002] Gas liquid two phase flow in pipes refers to the situation where free gas and liquid flow simultaneously in a pipeline. During this process, the free gas and liquid coexist in a gas-liquid state and transform into each other.

[0003] Among the four types of two-phase pipe flows (gas-liquid, gas-solid, liquid-liquid, and liquid-solid), the gas-liquid two-phase flow is the most complex. This is because there is a deformable two-phase interface in the gas-liquid two-phase flow, where the gas phase has high compressibility, and there is a dynamic process of mutual transformation between gas and liquid. The phenomenon that the two-phase interface is distributed in different geometric forms or different flow structure forms is called the flow pattern of two-phase flow, simply referred to as the flow state or flow pattern. Gas-liquid mixtures with different flow patterns follow their respective different flow laws.

[0004] In the prior art, the electrical signals output from temperature transmitters, pressure transmitters, differential pressure flowmeters, and swirl meters / microwave sensors need to be sent to a color screen recorder for processing and display. These electrical signals represent raw data such as temperature, pressure, local working condition flow rate, and working condition volume ratio, which cannot be directly used and have certain measurement errors affected by the working conditions. The required flow rate data output needs to be obtained through data processing. Summary of the Invention

[0005] To solve the above problems in the prior art, that is, to solve the problem of difficult measurement of gas-liquid two phases, the present invention provides a gas-liquid two-phase measurement method based on a fusion algorithm, including:

[0006] Step S10, collecting raw data in the target pipeline through a multi-sensor module as input data;

[0007] The input data includes the gas single-phase flow rate collected by a swirl meter, the total gas-liquid flow rate collected by a differential pressure flowmeter, the temperature collected by a temperature sensor, and the pressure collected by a pressure sensor;

[0008] Step S20, converting the input data through an A / D converter to obtain digital input data;

[0009] Step S30, preprocessing the digital input data to obtain preprocessed input data;

[0010] Step S40, extracting features from the preprocessed input data to obtain sensor measurement values;

[0011] Step S50, performing fusion calculation on the sensor measurement values through a multi-source data fusion algorithm to obtain an initial flow rate value;

[0012] Step S60: Based on the sensor measurement values, correct the flow value to be corrected through a high-precision matrix data iterative transmission algorithm between the swirl flow sensor and the differential pressure flowmeter to obtain the corrected two-phase flow value.

[0013] Further, the multi-module sensor module specifically includes:

[0014] The multi-sensor module includes a swirl flowmeter, a differential pressure flowmeter, a temperature sensor, and a pressure sensor. Further, the filtering and noise reduction specifically includes:

[0015] Reduce the noise in the digital input data by filtering and time-domain smoothing the digital input data to obtain the preprocessed input data.

[0016] Further, the feature extraction specifically includes:

[0017] Extract the sensor measurement values through the low-frequency component extraction method and the moving average method respectively.

[0018] Further, the multi-source data fusion algorithm specifically includes: Combine the sensor measurement values with the fluid physical property parameters, substitute them into the single-phase working condition flow calculation formula, and fuse and calculate the volume flow under the working condition. For the obtained volume flow under the working condition, convert it into the volume flow under the standard state according to the temperature collected by the temperature sensor and the pressure collected by the pressure sensor and output;

[0019] The standard state is specifically: temperature 0°C, pressure 1 atm (1 atmosphere).

[0020] The gas-phase mass flow calculation formula is: G g = A1f + A2

[0021] In the formula, f is the main frequency of the differential pressure signal of the swirl flowmeter, and A1 and A2 are empirical coefficients determined by experiments.

[0022] The single-phase working condition flow calculation formula is:

[0023] Among them, Q V represents the volume flow under the working condition; C represents the discharge coefficient; ε represents the flow expansion coefficient; d n represents the inner diameter of the throttling element; ΔP represents the pressure difference before and after the throttling element; ρ mix is the average density of the gas-liquid mixed flow object; D represents the inner diameter of the measuring tube; α represents the ratio of the inner diameter of the throttling element to the inner diameter of the measuring tube

[0024] The two-phase mass flow calculated according to the single-phase mass flow calculation formula is: G = Q V ρ mix

[0025] The two-phase flow rate calculated according to the single-phase flow rate calculation formula needs to be corrected by the over-reading factor Φ.

[0026] Average density ρ mix Calculation formula: ρ mix = β·ρ g +(1 - β)·ρ l

[0027] In the formula: β represents the volume gas holdup, which is the proportion of the gas phase volume in the total volume flow rate of the gas-liquid two-phase fluid; ρ g is the density of the gas phase in the gas-liquid two-phase fluid; ρ l is the density of the liquid phase in the gas-liquid two-phase fluid.

[0028] Furthermore, the conversion to standard conditions for output specifically includes:

[0029] The volume flow rate calculation formula under standard conditions is:

[0030] Among them, Q V represents the volume flow rate under working conditions; P a represents the local atmospheric pressure; P represents the gauge pressure measured by the pressure tapping hole of the flowmeter; P N represents the atmospheric pressure under standard conditions; T N represents the absolute temperature under standard conditions; T represents the absolute temperature of the fluid to be measured; Z N represents the compression factor of the gas under standard conditions; Z represents the compression factor of the gas under working conditions.

[0031] Furthermore, the high-precision matrix data iterative transmission algorithm between the swirling vortex sensor and the differential pressure flowmeter specifically includes:

[0032] The two-phase flow correction factor (over-reading factor) Φ based on differential pressure g Calculation formula:

[0033] The two-phase flow correction factor (under-reading factor) Φ' based on the swirling frequency g Calculation formula:

[0034] Where ΔP tp Indicates the pressure drop generated when two phases flow, ΔP g is the pressure drop when only single-phase gas flows in the pipeline, f tp represents the precession frequency of two-phase flow, f g It represents the precession frequency when assuming there is only single-phase gas in the tube; Φ g The smaller the under-reading factor is, the greater the deviation from the true frequency is; g is the acceleration of gravity; A is the cross-sectional area of ​​the pipe; D is the inner diameter of the pipe; G g is the gas phase mass flow rate; G l is the liquid mass flow rate.

[0035] Combined single-phase flow calculation formula, over-reading factor definition and calculation formula, under-reading factor definition and calculation formula, gas Froude number Fr g Calculation formula, Lockhart-Martinelli parameter X LM The liquid phase flow rate is calculated using the formula.

[0036] In another aspect of the present invention, a gas-liquid two-phase metering system based on a fusion algorithm is proposed, comprising a multi-sensor module, an A / D conversion module, a preprocessing module, a feature extraction module, a fusion calculation module, a data iteration module and a result output module;

[0037] A multi-sensor module, used to collect raw data in the target pipeline as input data;

[0038] An A / D conversion module is used to convert the input data through an A / D converter to obtain digital input data;

[0039] A preprocessing module, used for preprocessing the digital input data to obtain preprocessed input data;

[0040] A feature extraction module performs feature extraction on the preprocessed input data to obtain sensor measurement values;

[0041] A fusion calculation module, used for performing fusion calculation on the sensor measurement values ​​through a multi-source data fusion algorithm to obtain the two-phase flow value to be corrected;

[0042] A high-precision matrix data iterative module between a swirling vortex sensor and a differential pressure flowmeter is used to correct the flow value to be corrected based on the sensor measurement value through a high-precision matrix data iterative transmission algorithm between the swirling vortex sensor and the differential pressure flowmeter to obtain a corrected two-phase flow value;

[0043] A result output module is used to output the corrected two-phase flow value.

[0044] In a third aspect of the present invention, an electronic device is proposed, including:

[0045] At least one processor; and

[0046] A memory communicatively connected to at least one of the processors; wherein,

[0047] The memory stores instructions executable by the processor, and the instructions are used to be executed by the processor to implement the above-mentioned method for gas-liquid two-phase metering based on a fusion algorithm.

[0048] In a fourth aspect of the present invention, a computer-readable storage medium is proposed, and the computer-readable storage medium stores computer instructions, and the computer instructions are used to be executed by the computer to implement the above-mentioned method for gas-liquid two-phase metering based on a fusion algorithm.

[0049] Advantages of the present invention:

[0050] (1) The present invention integrates feature extraction and a fusion algorithm into the instrument to realize on-site gas-liquid two-phase metering;

[0051] (2) The present invention extracts more accurate and richer gas-liquid two-phase flows by integrating different data information. Description of the Drawings

[0052] By reading the detailed description of the non-restrictive embodiments with reference to the following drawings, other features, objects, and advantages of the present application will become more apparent:

[0053] Figure 1 is a flowchart of a method for gas-liquid two-phase metering based on a fusion algorithm of the present invention;

[0054] Figure 2 is a schematic diagram of a gas-liquid two-phase metering system based on a fusion algorithm of the present invention;

[0055] Figure 3 is a schematic diagram of the structure of a computer system of a server for implementing the method, system, and device embodiments of the present application. Detailed Embodiments

[0056] The present application will be further described in detail below in conjunction with the accompanying drawings and embodiments. It can be understood that the specific embodiments described herein are only used to explain the relevant invention, rather than limiting the invention. Additionally, it should be noted that for ease of description, only the parts related to the relevant invention are shown in the drawings.

[0057] It should be noted that, without conflict, the embodiments in the present application and the features in the embodiments can be combined with each other. The present application will be described in detail below with reference to the drawings and embodiments.

[0058] The present invention provides a gas-liquid two-phase metering method based on a fusion algorithm, including:

[0059] Step S10, collecting raw data in the target pipeline through a multi-sensor module as input data;

[0060] The input data includes the gas single-phase flow collected by a swirling vortex flowmeter, the total gas-liquid flow collected by a differential pressure flowmeter, the temperature collected by a temperature sensor, and the pressure collected by a pressure sensor;

[0061] Step S20, converting the input data through an A / D converter to obtain digital input data;

[0062] Step S30, preprocessing the digital input data to obtain preprocessed input data;

[0063] Step S40, extracting features from the preprocessed input data to obtain sensor measurement values;

[0064] Step S50, performing fusion calculation on the sensor measurement values through a multi-source data fusion algorithm to obtain an initial flow value;

[0065] Step S60, based on the sensor measurement values, correcting the flow value to be corrected through a high-precision matrix data iterative transmission algorithm between the swirling vortex sensor and the differential pressure flowmeter to obtain a corrected two-phase flow value.

[0066] To more clearly illustrate a gas-liquid two-phase metering method based on a fusion algorithm of the present invention, the following combines Figure 1 Details of each step in the embodiments of the present invention are elaborated.

[0067] A gas-liquid two-phase metering method based on a fusion algorithm according to the first embodiment of the present invention is described in detail as follows for each step:

[0068] Step S10, collecting raw data in the target pipeline through a multi-sensor module as input data;

[0069] The input data includes the gas single-phase flow rate collected by a swirl flowmeter, the total gas-liquid flow rate collected by a differential pressure flowmeter, the temperature collected by a temperature sensor, and the pressure collected by a pressure sensor;

[0070] In this embodiment, the multi-module sensor module specifically includes:

[0071] The multi-sensor module includes a swirl flowmeter, a differential pressure flowmeter, a temperature sensor, and a pressure sensor.

[0072] Step S20: Convert the input data through an A / D converter to obtain digital input data;

[0073] Step S30: Preprocess the digital input data to obtain preprocessed input data;

[0074] In this embodiment, the filtering and noise reduction specifically includes:

[0075] Reduce the noise in the digital input data by filtering and time-domain smoothing of the digital input data to obtain preprocessed input data.

[0076] Step S40: Extract features from the preprocessed input data to obtain sensor measurement values;

[0077] In this embodiment, the feature extraction specifically includes:

[0078] Extract the sensor measurement values through a low-frequency component extraction method and a moving average method respectively.

[0079] Step S50: Perform fusion calculation on the sensor measurement values through a multi-source data fusion algorithm to obtain an initial flow rate value;

[0080] In this embodiment, the multi-source data fusion algorithm specifically includes: Combine the sensor measurement values with the fluid physical property parameters, substitute them into the single-phase working condition flow rate calculation formula, perform fusion calculation to obtain the volume flow rate under the working condition, and for the obtained volume flow rate under the working condition, convert it into the volume flow rate under the standard state as compensation according to the temperature collected by the temperature sensor and the pressure collected by the pressure sensor and output;

[0081] The standard state is specifically: temperature 0°C, pressure 1 atm (1 atmospheric pressure).

[0082] The gas-phase mass flow rate calculation formula is: G g =A1f + A2

[0083] In the formula, f is the main frequency of the differential pressure signal of the swirl flowmeter, and A1 and A2 are empirical coefficients determined by experiments.

[0084] The volume flow rate calculation formula under single-phase conditions is as follows:

[0085] Among them, Q V represents the volume flow rate under the working conditions; C represents the discharge coefficient; ε represents the flow expansion coefficient; d n represents the inner diameter of the throttling element; ΔP represents the pressure difference before and after the throttling element; ρ mix is the average density of the gas-liquid mixed flowing object; D represents the inner diameter of the measuring pipe; α represents the ratio of the inner diameter of the throttling element to the inner diameter of the measuring pipe

[0086] The two-phase mass flow rate calculated according to the single-phase mass flow rate calculation formula is: G = Q V ρ mix

[0087] The two-phase flow rate calculated according to the single-phase flow rate calculation formula needs to be corrected by the over-reading factor Φ.

[0088] The calculation formula for the average density ρ mix is: ρ mix = β·ρ g +(1 - β)·ρ l

[0089] In the formula: β represents the volume gas content, which is the volume ratio of the gas phase in the total volume flow rate of the gas-liquid two-phase fluid; ρ g is the density of the gas phase in the gas-liquid two-phase fluid; ρ l is the density of the liquid phase in the gas-liquid two-phase fluid.

[0090] In this embodiment, the temperature collected by the temperature sensor and the pressure collected by the pressure sensor are used as compensation and converted into standard conditions for output, specifically including:

[0091] The volume flow rate calculation formula under standard conditions is:

[0092] Among them, Q V represents the volume flow rate under the working conditions; P a represents the local atmospheric pressure; P represents the gauge pressure measured by the pressure tapping hole of the flowmeter; P N represents the atmospheric pressure under standard conditions; T N represents the absolute temperature under standard conditions; T represents the absolute temperature of the fluid to be measured; Z N represents the compression coefficient of the gas under standard conditions; Z represents the compression coefficient of the gas under the working conditions;

[0093] When calibrating with a bell jar or negative pressure, take is the super compression factor, calculated according to the formula in the standard SY / T6143-1996 of China National Petroleum Corporation.

[0094] In this embodiment, the high-precision matrix data iterative transmission algorithm between the swirl sensor and the differential pressure flowmeter specifically includes:

[0095] Two-phase flow correction factor (overreading factor) based on differential pressure Φ g Calculation formula:

[0096] Two-phase flow correction coefficient (under-reading factor) Φ′ based on precession frequency g Calculation formula:

[0097] In the formula, ΔP tp Indicates the pressure drop generated when two phases flow, ΔP g is the pressure drop when only single-phase gas flows in the pipeline, f tp represents the precession frequency of two-phase flow, f g It represents the precession frequency when assuming there is only single-phase gas in the tube; Φ g The smaller the under-reading factor is, the greater the deviation from the true frequency is; g is the acceleration of gravity; A is the cross-sectional area of ​​the pipe; D is the inner diameter of the pipe; G g is the gas phase mass flow rate; G l is the liquid mass flow rate.

[0098] Combined single-phase flow calculation formula, over-reading factor definition and calculation formula, under-reading factor definition and calculation formula, gas Froude number Fr g Calculation formula, Lockhart-Martinelli parameter X LM The liquid phase flow rate is calculated using the formula.

[0099] Although the steps are described in the above order in the above embodiments, those skilled in the art can understand that in order to achieve the effects of this embodiment, different steps do not have to be executed in such an order, and they can be executed simultaneously (in parallel) or in a reversed order, and these simple changes are all within the protection scope of the present invention.

[0100] A gas-liquid two-phase metering system based on a fusion algorithm according to the second embodiment of the present invention, as Figure 2 shown, includes a multi-sensor module, an A / D conversion module, a preprocessing module, a feature extraction module, a fusion calculation module, a data iteration module, and a result output module:

[0101] The multi-sensor module is used to collect the original data in the target pipeline as input data;

[0102] The A / D conversion module is used to convert the input data through an A / D converter to obtain digital input data;

[0103] The preprocessing module is used to preprocess the digital input data to obtain preprocessed input data;

[0104] The feature extraction module extracts features from the preprocessed input data to obtain sensor measurement values;

[0105] The fusion calculation module is used to perform fusion calculation on the sensor measurement values through a multi-source data fusion algorithm to obtain a two-phase flow value to be corrected;

[0106] A high-precision matrix data iteration module between a swirling vortex sensor and a differential pressure flowmeter is used to correct the flow value to be corrected based on the sensor measurement values through a high-precision matrix data iteration transmission algorithm between the swirling vortex sensor and the differential pressure flowmeter to obtain a corrected two-phase flow value;

[0107] The result output module is used to output the corrected two-phase flow value.

[0108] Those skilled in the art of the technical field can clearly understand that for the convenience and brevity of description, the specific working process and related descriptions of the above-described system can refer to the corresponding process in the foregoing method embodiment, and will not be repeated here.

[0109] It should be noted that the gas-liquid two-phase metering system based on the fusion algorithm provided in the above embodiments is only illustrated by dividing the above functional modules. In practical applications, the above functions can be allocated to different functional modules according to needs, that is, the modules or steps in the embodiments of the present invention can be further decomposed or combined. For example, the modules in the above embodiments can be combined into one module, or further split into multiple sub-modules to complete all or part of the functions described above. The names of the modules and steps involved in the embodiments of the present invention are only used to distinguish each module or step, and are not regarded as an improper limitation of the present invention.

[0110] An electronic device according to a third embodiment of the present invention includes:

[0111] At least one processor; and

[0112] A memory communicatively connected to at least one of the processors; wherein,

[0113] The memory stores instructions executable by the processor, and the instructions are used to be executed by the processor to implement the above-mentioned gas-liquid two-phase metering method based on the fusion algorithm.

[0114] A computer-readable storage medium according to a fourth embodiment of the present invention stores computer instructions, and the computer instructions are used to be executed by the computer to implement the above-mentioned gas-liquid two-phase metering method based on the fusion algorithm.

[0115] Those skilled in the art can clearly understand that for the convenience and conciseness of description, the specific working processes and related descriptions of the above-mentioned storage devices and processing devices can refer to the corresponding processes in the foregoing method embodiments, and will not be repeated here.

[0116] Those skilled in the art should be able to realize that the modules and method steps of each example described in combination with the embodiments disclosed herein can be implemented by electronic hardware, computer software, or a combination of the two. The programs corresponding to the software modules and method steps can be placed in a random access memory (RAM), memory, read-only memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disk, removable disk, CD-ROM, or any other form of storage medium well-known in the technical field. To clearly illustrate the interchangeability of electronic hardware and software, the components and steps of each example have been generally described according to functions in the above description. Whether these functions are executed in the form of electronic hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementation should not be considered to exceed the scope of the present invention.

[0117] Refer to the following Figure 3 , which shows a schematic structural diagram of a computer system of a server for implementing the method, system, and device embodiments of the present application. Figure 3 The server shown is merely an example and should not impose any limitations on the functions and usage scope of the embodiments of the present application.

[0118] As Figure 3 shown, the computer system includes a central processing unit (CPU, Central Processing Unit) 601, which can perform various appropriate actions and processes according to the program stored in the read-only memory (ROM, Read Only Memory) 602 or the program loaded from the storage section 608 into the random access memory (RAM, Random Access Memory) 603. In the RAM 603, various programs and data required for system operation are also stored. The CPU 601, ROM 602, and RAM 603 are connected to each other via a bus 604. The input / output (I / O, Input / Output) interface 605 is also connected to the bus 604.

[0119] The following components are connected to the I / O interface 605: an input section 606 including a keyboard, a mouse, etc.; an output section 607 including, for example, a cathode ray tube (CRT, Cathode Ray Tube), a liquid crystal display (LCD, Liquid Crystal Display), etc., and a speaker, etc.; a storage section 608 including a hard disk, etc.; and a communication section 609 including a network interface card such as a LAN (Local Area Network) card, a modem, etc. The communication section 609 performs communication processing via a network such as the Internet. A drive 610 is also connected to the I / O interface 605 as needed. A removable medium 611, such as a magnetic disk, an optical disk, a magneto-optical disk, a semiconductor memory, etc., is installed on the drive 610 as needed so that a computer program read from it can be installed into the storage section 608 as needed.

[0120] In particular, according to an embodiment of the present disclosure, the processes described above with reference to the flowcharts can be implemented as computer software programs. For example, an embodiment of the present disclosure includes a computer program product that includes a computer program carried on a computer-readable medium, and the computer program includes program code for performing the method shown in the flowchart. In such an embodiment, the computer program can be downloaded and installed from a network through the communication section 609 and / or installed from the removable medium 611. When the computer program is executed by the central processing unit (CPU) 601, the above-described functions defined in the method of the present application are performed. It should be noted that the computer-readable medium in the present application can be a computer-readable signal medium or a computer-readable storage medium or any combination of the two. The computer-readable storage medium can be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination of the above. More specific examples of the computer-readable storage medium can include, but are not limited to: an electrical connection with one or more wires, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the above. In the present application, the computer-readable storage medium can be any tangible medium that contains or stores a program that can be used by or in conjunction with an instruction execution system, apparatus, or device. In the present application, the computer-readable signal medium can include a data signal propagated in a baseband or as part of a carrier wave, which carries the computer-readable program code. Such a propagated data signal can take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination of the above. The computer-readable signal medium can also be any computer-readable medium other than the computer-readable storage medium, which can send, propagate, or transmit a program for use by or in conjunction with an instruction execution system, apparatus, or device. The program code contained on the computer-readable medium can be transmitted by any appropriate medium, including but not limited to: wireless, wire, optical cable, RF, etc., or any suitable combination of the above.

[0121] Computer program code for performing the operations of this application can be written in one or more programming languages or combinations thereof. The above-mentioned programming languages include object-oriented programming languages such as Java, Smalltalk, C++, and also include conventional procedural programming languages such as the "C" language or similar programming languages. The program code can be executed entirely on the user's computer, partially on the user's computer, executed as a stand-alone software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In the case of a remote computer, the remote computer can be connected to the user's computer through any kind of network, including a local area network (LAN) or a wide area network (WAN), or it can be connected to an external computer (for example, by using an Internet service provider to connect through the Internet).

[0122] The flowcharts and block diagrams in the accompanying drawings illustrate the possible architectures, functions, and operations of systems, methods, and computer program products according to various embodiments of this application. In this regard, each block in the flowchart or block diagram can represent a module, a program segment, or a part of the code, and this module, program segment, or part of the code contains one or more executable instructions for implementing the specified logical function. It should also be noted that in some alternative implementations, the functions marked in the blocks may occur in a different order than marked in the accompanying drawings. For example, two consecutive blocks shown may actually be executed substantially in parallel, and they may sometimes be executed in the reverse order, depending on the functions involved. It should also be noted that each block in the block diagram and / or flowchart, and the combination of blocks in the block diagram and / or flowchart, can be implemented by a dedicated hardware-based system for performing the specified functions or operations, or can be implemented by a combination of dedicated hardware and computer instructions.

[0123] The terms "first", "second", etc. are used to distinguish similar objects and are not used to describe or represent a specific order or sequence.

[0124] The term "comprising" or any other similar term is intended to cover non-exclusive inclusion, so that a process, method, article, or device / equipment that includes a series of elements not only includes those elements, but also includes other elements not explicitly listed, or also includes elements inherent in these processes, methods, articles, or devices / equipment.

[0125] So far, the technical solution of the present invention has been described in conjunction with the preferred embodiments shown in the accompanying drawings. However, it is easily understood by those skilled in the art that the protection scope of the present invention is obviously not limited to these specific embodiments. Without departing from the principle of the present invention, those skilled in the art can make equivalent changes or substitutions to the relevant technical features, and the technical solutions after these changes or substitutions will all fall within the protection scope of the present invention.

Claims

1. A gas-liquid two-phase metering method based on a fusion algorithm, characterized in that, The method includes: Step S10, collecting raw data in the target pipeline through a multi-sensor module as input data; The input data includes the gas single-phase flow rate collected by a swirl flowmeter, the total gas-liquid flow rate collected by a differential pressure flowmeter, the temperature collected by a temperature sensor, and the pressure collected by a pressure sensor; Step S20, converting the input data through an A / D converter to obtain digital input data; Step S30, preprocessing the digital input data to obtain preprocessed input data; Step S40, extracting features from the preprocessed input data, and obtaining the low-frequency component of the preprocessed input data as the sensor measurement value; Step S50, performing fusion calculation on the sensor measurement values through a multi-source data fusion algorithm to obtain an initial flow rate value; Step S60, based on the sensor measurement values, correcting the flow rate value to be corrected through a high-precision matrix data iterative transmission algorithm between the swirl sensor and the differential pressure flowmeter to obtain a corrected two-phase flow rate value.

2. The gas-liquid two-phase metering method based on a fusion algorithm according to claim 1, wherein The multi-module sensor module specifically includes: The multi-sensor module includes a swirl flowmeter, a differential pressure flowmeter, a temperature sensor, and a pressure sensor.

3. A gas-liquid two-phase metering method based on a fusion algorithm according to claim 1, characterized in that The filtering and noise reduction specifically includes: Reducing the noise in the digital input data by filtering and time-domain smoothing of the digital input data to obtain preprocessed input data.

4. A gas-liquid two-phase metering method based on a fusion algorithm according to claim 1, characterized in that The feature extraction specifically includes: Respectively extracting the sensor measurement values, including the temperature value of the temperature sensor, the pressure value of the pressure sensor, the differential pressure value of the orifice flowmeter, and the frequency value of the swirl flowmeter, through a low-frequency component extraction method and a moving average method.

5. A gas-liquid two-phase metering method based on a fusion algorithm according to claim 1, characterized in that The multi-source data fusion algorithm specifically includes: Substituting the sensor measurement values and fluid physical property parameters into the single-phase working condition flow rate calculation formula, and performing fusion calculation to obtain the volume flow rate under the working condition. For the obtained volume flow rate under the working condition, it is compensated and converted into the volume flow rate under the standard state according to the temperature collected by the temperature sensor and the pressure collected by the pressure sensor and output; the specific standard state is: temperature 0°C, pressure 1 atm; The gas-phase mass flow rate calculation formula is: G g = A1f + A2 In the formula, f is the main frequency of the differential pressure signal of the swirl flowmeter, and A1 and A2 are empirical coefficients determined by experiments; The gas phase flow rate in the gas-liquid two-phase flow calculated according to the gas phase mass flow rate calculation formula needs to be corrected with the under-reading factor Φ g ′; The two-phase volume flow rate calculated according to the single-phase working condition volume flow rate calculation formula is: Among them, Q V represents the volumetric flow rate under the operating conditions; C represents the discharge coefficient; ε represents the expansion coefficient of the flow stream; d n represents the inner diameter of the throttling element; ΔP represents the pressure difference across the throttling element; ρ mix is the average density of the gas-liquid mixed flowing object; D represents the inner diameter of the measuring tube; α represents the ratio of the inner diameter of the throttling element to the inner diameter of the measuring tube The two-phase mass flow rate calculated according to the single-phase mass flow rate calculation formula is: G = Q V ρ mix The two-phase flow rate calculated according to the single-phase flow rate calculation formula needs to be corrected by the over-reading factor Φ g Correction; Average density ρ mix Calculation formula: ρ mix = β·ρ g + (1 - β)·ρ l In the formula: β represents the volume gas holdup, which is the volume proportion of the gas phase in the total volume flow rate of the gas-liquid two-phase fluid; ρ g is the density of the gas phase in the gas-liquid two-phase fluid; ρ l is the density of the liquid phase in the gas-liquid two-phase fluid.

6. The gas-liquid two-phase metering method based on a fusion algorithm according to claim 5, wherein, The conversion to standard condition output specifically includes: The volume flow rate calculation formula under the standard state is: Among them, Q V represents the volume flow rate under the operating condition; P a represents the local atmospheric pressure; P represents the gauge pressure measured by the pressure tapping hole of the flowmeter; P N represents the atmospheric pressure under standard conditions; T N represents the absolute temperature under standard conditions; T represents the absolute temperature of the fluid to be measured; Z N represents the compression factor of the gas under standard conditions; Z represents the compression factor of the gas under the operating condition.

7. A gas-liquid two-phase metering method based on a fusion algorithm according to claim 5, characterized in that The high-precision matrix data iterative transmission algorithm between the swirl sensor and the differential pressure flowmeter specifically includes: Over-reading factor Φ g is defined as the two-phase flow correction factor based on differential pressure: Where, ΔP tp represents the pressure drop generated during two-phase flow, and ΔP g is the pressure drop assuming that only single-phase gas flows in the pipeline; Φ g Calculation formula: where X LM is the Lockhart-Martinelli parameter, reflecting the relative magnitude of the gas-phase velocity and the liquid-phase velocity; Fr g is the gas Froude number, reflecting the internal relationship among the gas-phase flow rate, pressure, and density factors; the formula is: Coefficient calculation formula: where g is the acceleration due to gravity; A is the cross-sectional area of the pipeline; D is the inner diameter of the pipeline; G g is the gas-phase mass flow rate; G l is the liquid-phase mass flow rate; β represents the volume gas holdup, which is the proportion of the gas-phase volume in the total volume flow rate of the gas-liquid two-phase fluid; ρ g is the density of the gas phase in the gas-liquid two-phase fluid; ρ l is the density of the liquid phase in the gas-liquid two-phase fluid; Under-reading factor Φ′ g Defined as the two-phase flow correction coefficient based on the precession frequency: where, f tp represents the precession frequency in two-phase flow (the main frequency of differential pressure measured by the precession vortex flowmeter), and f g represents the precession frequency assuming only single-phase gas in the pipe; Φ′ g The smaller the under-reading factor, the greater the deviation from the true frequency; Φ′ g Calculation formula: Among them, X LM , Fr g are defined in the same way as those in the calculation formula of Φ g ; Simultaneously establish the gas-phase mass flow rate calculation formula, single-phase operating condition flow rate calculation formula, over-reading factor definition formula and calculation formula, under-reading factor definition formula and calculation formula, gas Froude number Fr g calculation formula, Lockhart-Martinelli parameter X LM Calculate the liquid-phase flow rate value using the calculation formula.

8. A gas-liquid two-phase metering system based on a fusion algorithm, characterized in that, The system includes a multi-sensor module, an A / D conversion module, a preprocessing module, a feature extraction module, a fusion calculation module, a data iteration module, and a result output module; The multi-sensor module is used to collect raw data in the target pipeline as input data; The A / D conversion module is used to convert the input data through an A / D converter to obtain digital input data; The preprocessing module is used to preprocess the digital input data to obtain preprocessed input data; A feature extraction module that extracts features from the preprocessed input data to obtain sensor measurement values; A fusion calculation module for performing fusion calculation on the sensor measurement values through a multi-source data fusion algorithm to obtain the two-phase flow values to be corrected; A high-precision matrix data iterative module between a swirl flow sensor and a differential pressure flowmeter, which is used to correct the flow values to be corrected based on the sensor measurement values through a high-precision matrix data iterative transmission algorithm between the swirl flow sensor and the differential pressure flowmeter to obtain the corrected two-phase flow values; A result output module for outputting the corrected two-phase flow values.

9. An electronic device, characterized in that, Comprising: At least one processor; And A memory communicatively connected to at least one of the processors; wherein, The memory stores instructions executable by the processor, and the instructions are used to be executed by the processor to implement a method for gas-liquid two-phase metering based on a fusion algorithm according to any one of claims 1-7.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores computer instructions, and the computer instructions are used to be executed by the computer to implement a method for gas-liquid two-phase metering based on a fusion algorithm according to any one of claims 1-7.

Images