Method and device for determining mass flow rate of liquid or supercritical carbon dioxide

Abstract

The invention relates to the technical field of carbon dioxide flow calculation, and provides a method and device for determining the mass flow of liquid or supercritical carbon dioxide. The method comprises the following steps: acquiring real-time parameters of liquid and / or supercritical target carbon dioxide under real-time temperature and real-time pressure; calculating the real-time density of the target carbon dioxide at the real-time temperature and the real-time pressure based on the real-time parameters; and calculating the mass flow of the target carbon dioxide by combining the real-time density according to a flow parameter measured by a preset flow meter. According to the embodiment of the invention, the calculation precision of the mass flow rate of the carbon dioxide in the liquid state and / or the supercritical state can be greatly improved on the premise of not increasing a large amount of cost.

Description

Technical Field

[0001] The present disclosure relates to the technical field of carbon dioxide flow calculation, and in particular to a method and device for determining the mass flow of liquid or supercritical carbon dioxide. Background Art

[0002] With the development of CCUS (Carbon Capture, Utilization and Storage), my country has also begun to vigorously support the development and promotion of CCUS. my country's coal-based energy structure will put my country under great pressure to reduce carbon emissions. In the next few decades, when international climate policies become clearer, CCUS technology will begin to be applied on a large scale, forming a large application market. The important issue that follows is the accurate measurement of carbon dioxide flow.

[0003] Liquid and / or supercritical carbon dioxide has special properties. Its density is a function of temperature and pressure. Its density varies greatly at different temperatures and pressures, and the density variation range at the same temperature and different pressures is also very large. The measurement error of general volume flow meters is large. Only Coriolis mass flow meters can directly detect mass flow and density, and their measurement results are not affected by density changes. Other flow meters detect volume flow and convert mass flow through density. However, Coriolis mass flow meters are expensive and cannot be widely used and promoted.

[0004] For conventional gas flow measurement, the ideal gas state equation can be used to convert the working volume flow into the standard volume flow (standard condition refers to 20°C, 1 standard atmospheric pressure), and the standard volume flow is used as the standard to avoid the impact of changes in temperature, pressure and density on the volume flow. For crude oil and refined oil flow measurement, the temperature correction coefficient and pressure correction coefficient can be queried or obtained in accordance with the national standard GB / T 1885 "Oil Volume Correction Coefficient Table", thereby calculating the net standard volume. The net standard weight of the oil in the air can also be obtained as the measurement result through the oil density and air buoyancy correction coefficient.

[0005] As for dense phase and supercritical carbon dioxide, it is not in gaseous state and cannot be calculated using the ideal gas state equation. At the same time, there is a lack of corresponding standard query temperature and pressure correction coefficients. In addition, its density is greatly affected by temperature and pressure, resulting in volume measurement cannot accurately represent its true flow rate. Current flow meters do not support dense phase and supercritical carbon dioxide mass flow density compensation calculations, but use density as a constant in the calculation, resulting in large measurement errors. Summary of the invention

[0006] In view of this, the embodiments of the present disclosure provide a method and device for determining the mass flow rate of liquid or supercritical carbon dioxide to solve the problem that flow meters in the prior art do not support density compensation calculation of dense phase and supercritical carbon dioxide mass flow, but instead use density as a constant to participate in the calculation, resulting in large measurement errors.

[0007] According to a first aspect of an embodiment of the present disclosure, a method for determining a mass flow rate of liquid or supercritical carbon dioxide is provided, comprising:

[0008] Acquire real-time parameters of target carbon dioxide in liquid state and / or supercritical state at real-time temperature and real-time pressure;

[0009] Calculate the real-time density of the target carbon dioxide at the real-time temperature and the real-time pressure based on the real-time parameters;

[0010] The mass flow rate of the target carbon dioxide is calculated according to the flow parameter measured by the preset flow meter in combination with the real-time density.

[0011] In some embodiments, the real-time parameters include a gas constant, a critical temperature, a critical pressure, an eccentricity factor, and a pseudo-comparison temperature, wherein the real-time density of the target carbon dioxide at the real-time temperature and the real-time pressure is calculated based on the real-time parameters, including:

[0012] Calculate the first coefficient and the second coefficient respectively according to the gas constant, the critical temperature and the critical pressure in the real-time parameters;

[0013] Calculating a simulated comparison temperature according to the real-time temperature and a critical temperature in the real-time parameters;

[0014] Calculating a third coefficient according to the proposed comparison temperature and an eccentricity factor in the real-time parameter;

[0015] Calculate a fourth coefficient based on the first coefficient and the third coefficient;

[0016] The real-time density of the target carbon dioxide is calculated according to the second coefficient, the fourth coefficient, the real-time temperature, and the gas constant in the real-time parameters.

[0017] In some embodiments, the calculating the first coefficient and the second coefficient respectively according to the gas constant, the critical temperature and the critical pressure in the real-time parameters includes:

[0018] The first coefficient is calculated according to a preset first calculation formula, and the gas constant, critical temperature and critical pressure in the real-time parameters, wherein the first calculation formula includes:

[0019]

[0020] The second coefficient is calculated according to a preset second calculation formula, and the gas constant, critical temperature and critical pressure in the real-time parameters, wherein the second calculation formula includes:

[0021]

[0022] Among them, a c represents the first coefficient, b represents the second coefficient, R represents the gas constant in the real-time parameter, T c represents the critical temperature in real-time parameters, p c Indicates the critical pressure in real-time parameters.

[0023] In some embodiments, the calculating the third coefficient according to the simulated comparison temperature and the eccentricity factor in the real-time parameter includes:

[0024] The third coefficient is calculated according to a preset third calculation formula, the simulated comparison temperature, and the eccentricity factor in the real-time parameter, wherein the third calculation formula includes:

[0025]

[0026] Wherein, α represents the third coefficient, ω represents the eccentricity factor in the real-time parameter, T r represents the proposed comparison temperature.

[0027] In some embodiments, calculating the real-time density of carbon dioxide according to the second coefficient, the fourth coefficient, the real-time temperature, and the gas constant in the real-time parameter includes:

[0028] The real-time density of carbon dioxide is calculated according to a preset fourth calculation formula, the second coefficient, the fourth coefficient, the real-time temperature, and the gas constant in the real-time parameter, wherein the fourth calculation formula includes:

[0029]

[0030] Wherein, p represents the real-time pressure, R represents the gas constant in the real-time parameters, T represents the real-time temperature, a represents the fourth coefficient, b represents the second coefficient, and ρ represents the real-time density of the carbon dioxide.

[0031] In some embodiments, the preset flow meter is a differential pressure flow meter, and the flow parameter measured by the differential pressure flow meter is the differential pressure, wherein the mass flow rate of the target carbon dioxide is calculated according to the differential pressure in combination with the real-time density, including:

[0032] According to the preset fifth calculation formula, the mass flow rate of carbon dioxide is calculated according to the real-time density, and the fifth calculation formula includes:

[0033]

[0034] Among them, q m represents the mass flow rate, C D represents the outflow coefficient, D represents the throttling aperture, ε represents the expansion coefficient, ρ represents the real-time density, Δp represents the differential pressure, β represents the ratio of the throttling aperture to the inner diameter of the measuring tube, and d represents the preset calculation coefficient.

[0035] According to a second aspect of the embodiments of the present disclosure, there is provided a device for determining the mass flow rate of liquid or supercritical carbon dioxide, comprising:

[0036] An acquisition module, used for acquiring real-time parameters of target carbon dioxide in liquid state and / or supercritical state at real-time temperature and real-time pressure;

[0037] A first calculation module, used for calculating the real-time density of the target carbon dioxide at the real-time temperature and the real-time pressure based on the real-time parameters;

[0038] The second calculation module is used to calculate the mass flow rate of the target carbon dioxide according to the flow parameter measured by the preset flow meter in combination with the real-time density.

[0039] In some embodiments, the real-time parameters include a gas constant, a critical temperature, a critical pressure, an eccentricity factor, and a pseudo-comparison temperature, wherein the first calculation module is further configured to:

[0040] Calculate the first coefficient and the second coefficient respectively according to the gas constant, the critical temperature and the critical pressure in the real-time parameters;

[0041] Calculating a simulated comparison temperature according to the real-time temperature and a critical temperature in the real-time parameters;

[0042] Calculating a third coefficient according to the proposed comparison temperature and an eccentricity factor in the real-time parameter;

[0043] Calculate a fourth coefficient based on the first coefficient and the third coefficient;

[0044] The real-time density of the target carbon dioxide is calculated according to the second coefficient, the fourth coefficient, the real-time temperature, and the gas constant in the real-time parameters.

[0045] In some embodiments, the calculating the first coefficient and the second coefficient respectively according to the gas constant, the critical temperature and the critical pressure in the real-time parameters includes:

[0046] The first coefficient is calculated according to a preset first calculation formula, and the gas constant, critical temperature and critical pressure in the real-time parameters, wherein the first calculation formula includes:

[0047]

[0048] The second coefficient is calculated according to a preset second calculation formula, and the gas constant, critical temperature and critical pressure in the real-time parameters, wherein the second calculation formula includes:

[0049]

[0050] Among them, a c represents the first coefficient, b represents the second coefficient, R represents the gas constant in the real-time parameter, T c represents the critical temperature in real-time parameters, p c Indicates the critical pressure in real-time parameters.

[0051] In some embodiments, the calculating the third coefficient according to the simulated comparison temperature and the eccentricity factor in the real-time parameter includes:

[0052] The third coefficient is calculated according to a preset third calculation formula, the simulated comparison temperature, and the eccentricity factor in the real-time parameter, wherein the third calculation formula includes:

[0053]

[0054] Wherein, α represents the third coefficient, ω represents the eccentricity factor in the real-time parameter, T r represents the proposed comparison temperature.

[0055] In some embodiments, calculating the real-time density of carbon dioxide according to the second coefficient, the fourth coefficient, the real-time temperature, and the gas constant in the real-time parameter includes:

[0056] The real-time density of carbon dioxide is calculated according to a preset fourth calculation formula, the second coefficient, the fourth coefficient, the real-time temperature, and the gas constant in the real-time parameter, wherein the fourth calculation formula includes:

[0057]

[0058] Wherein, p represents the real-time pressure, R represents the gas constant in the real-time parameters, T represents the real-time temperature, a represents the fourth coefficient, b represents the second coefficient, and ρ represents the real-time density of the carbon dioxide.

[0059] In some embodiments, the preset flow meter is a differential pressure flow meter, and the flow parameter measured by the differential pressure flow meter is the differential pressure, wherein the mass flow rate of the target carbon dioxide is calculated according to the differential pressure in combination with the real-time density, including:

[0060] According to the preset fifth calculation formula, the mass flow rate of carbon dioxide is calculated according to the real-time density, and the fifth calculation formula includes:

[0061]

[0062] Among them, q m represents the mass flow rate, C D represents the outflow coefficient, D represents the throttling aperture, ε represents the expansion coefficient, ρ represents the real-time density, Δp represents the differential pressure, β represents the ratio of the throttling aperture to the inner diameter of the measuring tube, and d represents the preset calculation coefficient.

[0063] According to a third aspect of an embodiment of the present disclosure, an electronic device is provided, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the steps of the above method when executing the computer program.

[0064] According to a fourth aspect of the embodiments of the present disclosure, a computer-readable storage medium is provided, which stores a computer program, and when the computer program is executed by a processor, the steps of the above method are implemented.

[0065] According to a fifth aspect of the embodiments of the present disclosure, a computer program product is provided. The computer program product includes a computer program or instructions. When the computer program or instructions are executed by a processor, the steps of the above method are implemented.

[0066] Beneficial Effects

[0067] Compared with the prior art, the embodiments of the present disclosure have at least the following beneficial effects: by acquiring real-time parameters of target carbon dioxide in liquid and / or supercritical state at real-time temperature and real-time pressure; calculating the real-time density of the target carbon dioxide at real-time temperature and real-time pressure based on the real-time parameters; calculating the mass flow rate of the target carbon dioxide in combination with the real-time density according to the flow parameters measured by a preset flow meter, so that the mass flow rate can be accurately calculated in real time in combination with the change of density without using a Coriolis flowmeter, and the calculation accuracy of the mass flow rate of carbon dioxide in liquid and / or supercritical state is greatly improved without increasing large costs. BRIEF DESCRIPTION OF THE DRAWINGS

[0068] In order to more clearly illustrate the technical solutions in the embodiments of the present disclosure, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present disclosure. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0069] Figure 1 It is a schematic diagram of an application scenario of the method for determining the mass flow rate of liquid or supercritical carbon dioxide provided in an embodiment of the present disclosure;

[0070] Figure 2 is a flow chart of some embodiments of a method for determining the mass flow rate of liquid or supercritical carbon dioxide provided according to an embodiment of the present disclosure;

[0071] Figure 3 is a flow chart of other embodiments of another method for determining the mass flow rate of liquid or supercritical carbon dioxide provided according to an embodiment of the present disclosure;

[0072] Figure 4 It is a simplified structural schematic diagram of a device for determining the mass flow rate of liquid or supercritical carbon dioxide provided according to an embodiment of the present disclosure;

[0073] Figure 5 is a schematic diagram of an electronic device provided according to an embodiment of the present disclosure. DETAILED DESCRIPTION

[0074] Embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. Although certain embodiments of the present disclosure are shown in the accompanying drawings, it should be understood that the present disclosure can be implemented in various forms and should not be construed as being limited to the embodiments set forth herein. On the contrary, these embodiments are provided to provide a more thorough and complete understanding of the present disclosure. It should be understood that the drawings and embodiments of the present disclosure are only for exemplary purposes and are not intended to limit the scope of protection of the present disclosure.

[0075] It should also be noted that, for ease of description, only the parts related to the present disclosure are shown in the drawings. In the absence of conflict, the embodiments and features in the embodiments of the present disclosure can be combined with each other.

[0076] It should be noted that the concepts such as "first" and "second" mentioned in the present disclosure are only used to distinguish different systems, devices, modules or units, and are not used to limit the order or interdependence of the functions performed by these systems, devices, modules or units.

[0077] It should be noted that the modifications of "one" and "plurality" mentioned in the present disclosure are illustrative rather than restrictive, and those skilled in the art should understand that unless otherwise clearly indicated in the context, it should be understood as "one or more".

[0078] The names of the messages or information exchanged between multiple devices in the embodiments of the present disclosure are only used for illustrative purposes and are not used to limit the scope of these messages or information.

[0079] The present disclosure will be described in detail below with reference to the accompanying drawings and in conjunction with embodiments.

[0080] Figure 1 It is a schematic diagram of an application scenario of a method for determining the mass flow rate of liquid or supercritical carbon dioxide according to some embodiments of the present disclosure.

[0081] exist Figure 1 In the application scenario, first, the computing device 101 can obtain real-time parameters 102 of the target carbon dioxide in liquid state and / or supercritical state at real-time temperature and real-time pressure.

[0082] Secondly, the computing device 101 may calculate a real-time density 103 of the target carbon dioxide at a real-time temperature and a real-time pressure based on the real-time parameter 102 .

[0083] Finally, the computing device 101 may calculate the mass flow rate 105 of the target carbon dioxide according to the flow parameter 104 measured by the preset flow meter and in combination with the real-time density 103 .

[0084] It should be noted that the computing device 101 can be hardware or software. When the computing device is hardware, it can be a single server or a single terminal device, or a programmable logic controller (PLC) or a distributed control system (DCS). When the computing device is embodied as software, it can be installed in the hardware devices listed above. It can be implemented as, for example, a single software or a software module. No specific limitation is made here.

[0085] It should be understood that Figure 1 The number of computing devices in the embodiment is only illustrative. Any number of computing devices may be provided according to implementation requirements.

[0086] Continue to refer Figure 2 , shows a process 200 of some embodiments of the method for determining the mass flow rate of liquid or supercritical carbon dioxide according to the present disclosure. The method can be Figure 1 The method for determining the mass flow rate of liquid or supercritical carbon dioxide comprises the following steps:

[0087] Step 201 : obtaining real-time parameters of target carbon dioxide in liquid state and / or supercritical state at real-time temperature and real-time pressure.

[0088] In some embodiments, the execution subject of the method for determining the mass flow rate of liquid or supercritical carbon dioxide (e.g. Figure 1 The computing device 101) shown can be connected to the target device via a wired connection or a wireless connection, and then obtain real-time parameters of the target carbon dioxide in liquid and / or supercritical state at real-time temperature and real-time pressure.

[0089] In the process of oil development, injecting carbon dioxide into oil reservoirs can not only improve oil and gas recovery rates, but also reduce greenhouse gas emissions. It is an important measure for carbon emission reduction. During the distribution and transportation of carbon dioxide, it often needs to be compressed into a liquid or supercritical state (that is, the physical properties are close to those of a gas, and the density is close to that of a liquid). Since carbon dioxide has different densities at different temperatures and pressures. The target carbon dioxide is liquid and / or supercritical carbon dioxide whose mass flow rate needs to be detected. The real-time temperature may refer to the real-time temperature detected in the carbon dioxide transmission pipeline. The real-time pressure may refer to the real-time pressure detected in the carbon dioxide transmission pipeline. The real-time parameters may refer to various parameters of the detectable liquid and / or supercritical carbon dioxide.

[0090] It should be noted that the above-mentioned wireless connection methods may include but are not limited to 3G / 4G / 5G connection, WiFi connection, Bluetooth connection, WiMAX connection, Zigbee connection, UWB (ultra wideband) connection, and other wireless connection methods currently known or developed in the future.

[0091] Step 202: Calculate the real-time density of the target carbon dioxide at the real-time temperature and real-time pressure based on the real-time parameters.

[0092] In some embodiments, the execution entity may calculate the real-time density of the target carbon dioxide at the real-time temperature and the real-time pressure based on the real-time parameters.

[0093] Due to the special physical properties of liquid and / or supercritical carbon dioxide, the density is a function of temperature and pressure. Its density varies greatly at different temperatures and pressures, and the density variation range at the same temperature and different pressures is also very large. Therefore, when accurately calculating carbon dioxide, it is necessary to first determine the density of carbon dioxide that changes in real time, that is, the above-mentioned real-time density.

[0094] In some optional implementations, the real-time parameters may include at least a gas constant, a critical temperature, a critical pressure, an eccentricity factor, and a pseudo-comparison temperature. The above parameters are common parameters in the art and will not be described in detail.

[0095] In some optional implementations, the execution entity may calculate the real-time density of the target carbon dioxide at the real-time temperature and real-time pressure based on the real-time parameters according to the following steps:

[0096] In the first step, the first coefficient and the second coefficient are calculated respectively according to the gas constant, the critical temperature and the critical pressure in the real-time parameters. The first coefficient and the second coefficient may refer to two different intermediate coefficients obtained in the calculation process.

[0097] As some preferred implementations, the first coefficient and the second coefficient are calculated respectively according to the gas constant, the critical temperature and the critical pressure in the real-time parameters, including:

[0098] The first coefficient is calculated according to a preset first calculation formula, and the gas constant, critical temperature and critical pressure in the real-time parameters, wherein the first calculation formula includes:

[0099]

[0100] The second coefficient is calculated according to a preset second calculation formula, and the gas constant, critical temperature and critical pressure in the real-time parameters, wherein the second calculation formula includes:

[0101]

[0102] Among them, a c represents the first coefficient, b represents the second coefficient, R represents the gas constant in the real-time parameter, T c represents the critical temperature in real-time parameters, p c Indicates the critical pressure in real-time parameters.

[0103] The first calculation formula and the second calculation formula may refer to preset calculation formulas for calculating the first coefficient and the second coefficient, respectively.

[0104] The second step is to calculate the simulated comparison temperature according to the real-time temperature and the critical temperature in the real-time parameters. The simulated comparison temperature may refer to the temperature simulated by calculation.

[0105] For example, the pseudo-comparison temperature can be calculated by the following calculation formula:

[0106]

[0107] Wherein, Tr represents the simulated comparison temperature, T represents the real-time temperature, and Tc represents the critical temperature.

[0108] The third step is to calculate a third coefficient according to the proposed comparison temperature and the eccentricity factor in the real-time parameter. The third coefficient may refer to an intermediate coefficient obtained during the calculation process.

[0109] As some preferred implementations, the calculating the third coefficient according to the simulated comparison temperature and the eccentricity factor in the real-time parameter includes:

[0110] The third coefficient is calculated according to a preset third calculation formula, the simulated comparison temperature, and the eccentricity factor in the real-time parameter, wherein the third calculation formula includes:

[0111]

[0112] Wherein, α represents the third coefficient, ω represents the eccentricity factor in the real-time parameter, T r represents the proposed comparison temperature.

[0113] The fourth step is to calculate a fourth coefficient according to the first coefficient and the third coefficient. The fourth coefficient may refer to an intermediate coefficient obtained during the calculation process.

[0114] For example, the fourth coefficient can be calculated by the following formula:

[0115] a=a c α

[0116] Wherein, a represents the fourth coefficient, a c represents the first coefficient, α represents and belongs to the third coefficient.

[0117] The fifth step is to calculate the real-time density of the target carbon dioxide according to the second coefficient, the fourth coefficient, the real-time temperature, and the gas constant in the real-time parameters.

[0118] As some preferred implementations, the calculating the real-time density of carbon dioxide according to the second coefficient, the fourth coefficient, the real-time temperature, and the gas constant in the real-time parameter includes:

[0119] The real-time density of carbon dioxide is calculated according to a preset fourth calculation formula, the second coefficient, the fourth coefficient, the real-time temperature, and the gas constant in the real-time parameter, wherein the fourth calculation formula includes:

[0120]

[0121] Wherein, p represents the real-time pressure, R represents the gas constant in the real-time parameters, T represents the real-time temperature, a represents the fourth coefficient, b represents the second coefficient, and ρ represents the real-time density of the carbon dioxide.

[0122] Step 203: Calculate the mass flow rate of the target carbon dioxide according to the flow parameter measured by the preset flow meter in combination with the real-time density.

[0123] Existing non-Coriolis mass flowmeter flowmeters include at least: differential pressure flowmeter, velocity flowmeter, volume flowmeter, etc. Since different flowmeters measure different parameters, their corresponding calculation formulas are also different.

[0124] In some embodiments, the preset flow meter may be a differential pressure flow meter, and the flow parameter measured by the differential pressure flow meter is the differential pressure.

[0125] In some optional implementations, calculating the target carbon dioxide mass flow rate according to the differential pressure in combination with the real-time density includes:

[0126] According to the preset fifth calculation formula, the mass flow rate of carbon dioxide is calculated according to the real-time density, and the fifth calculation formula includes:

[0127]

[0128] Among them, q m represents the mass flow rate, C D represents the outflow coefficient, D represents the throttling aperture, ε represents the expansion coefficient, ρ represents the real-time density, Δp represents the differential pressure, β represents the ratio of the throttling aperture to the inner diameter of the measuring tube, and d represents the preset calculation coefficient.

[0129] In other embodiments, the preset flow meter may be a velocity flow meter or a volume flow meter, which may directly calculate (or calculate in combination with the real-time density) the volume flow rate of the target carbon dioxide in liquid and / or supercritical state, and then calculate the mass flow rate in combination with the real-time density. The corresponding calculation formula for the volume flow rate is prior art and will not be described in detail here.

[0130] The beneficial effects of one of the above-mentioned embodiments of the present disclosure include at least: obtaining real-time parameters of the target carbon dioxide in liquid and / or supercritical state at real-time temperature and real-time pressure; calculating the real-time density of the target carbon dioxide at real-time temperature and real-time pressure based on the real-time parameters; calculating the mass flow rate of the target carbon dioxide in combination with the real-time density according to the flow parameters measured by a preset flow meter, so that the mass flow rate can be accurately calculated in real time in combination with the change in density without using a Coriolis flowmeter, and the calculation accuracy of the mass flow rate of carbon dioxide in liquid and / or supercritical state is greatly improved without increasing large costs.

[0131] Continue to refer Figure 3 , shows a process 300 of another embodiment of a method for determining the mass flow rate of liquid or supercritical carbon dioxide according to the present disclosure, which can be performed by Figure 1The method for determining the mass flow rate of liquid or supercritical carbon dioxide includes:

[0132] Step 301, obtaining real-time parameters of target carbon dioxide in liquid state and / or supercritical state at real-time temperature and real-time pressure, wherein the real-time parameters include gas constant, critical temperature, critical pressure, eccentricity factor, and pseudo-comparison temperature.

[0133] Step 302, calculating the first coefficient and the second coefficient respectively according to the gas constant, critical temperature and critical pressure in the real-time parameters.

[0134] Step 303: Calculate the simulated comparison temperature according to the real-time temperature and the critical temperature in the real-time parameters.

[0135] Step 304: Calculate a third coefficient according to the simulated comparison temperature and the eccentricity factor in the real-time parameter.

[0136] Step 305: Calculate a fourth coefficient according to the first coefficient and the third coefficient.

[0137] Step 306: Calculate the real-time density of the target carbon dioxide according to the second coefficient, the fourth coefficient, the real-time temperature, and the gas constant in the real-time parameters.

[0138] Step 307: Calculate the mass flow rate of the target carbon dioxide according to the flow parameter measured by the preset flow meter in combination with the real-time density.

[0139] In some embodiments, the specific implementation of steps 301-307 and the technical effects brought about can be referred to Figure 2 The steps in the corresponding embodiments will not be repeated here.

[0140] All the above optional technical solutions can be arbitrarily combined to form optional embodiments of the present application, which will not be described one by one here.

[0141] The following are embodiments of the device disclosed herein, which can be used to execute the method embodiments disclosed herein. For details not disclosed in the device embodiments disclosed herein, please refer to the method embodiments disclosed herein.

[0142] Further references Figure 4 As an implementation of the above methods in the above figures, the present disclosure provides some embodiments of a device for determining the mass flow rate of liquid or supercritical carbon dioxide. These device embodiments are similar to Figure 2 The above method embodiments correspond.

[0143] like Figure 4 As shown, in some embodiments, the mass flow rate determination device 400 of liquid or supercritical carbon dioxide includes:

[0144] An acquisition module, used for acquiring real-time parameters of target carbon dioxide in liquid state and / or supercritical state at real-time temperature and real-time pressure;

[0145] A first calculation module, used for calculating the real-time density of the target carbon dioxide at the real-time temperature and the real-time pressure based on the real-time parameters;

[0146] The second calculation module is used to calculate the mass flow rate of the target carbon dioxide according to the flow parameter measured by the preset flow meter in combination with the real-time density.

[0147] In some optional implementations of some embodiments, the real-time parameters include a gas constant, a critical temperature, a critical pressure, an eccentricity factor, and a pseudo-comparison temperature, wherein the first calculation module is further configured to:

[0148] Calculate the first coefficient and the second coefficient respectively according to the gas constant, the critical temperature and the critical pressure in the real-time parameters;

[0149] Calculating a simulated comparison temperature according to the real-time temperature and a critical temperature in the real-time parameters;

[0150] Calculating a third coefficient according to the proposed comparison temperature and an eccentricity factor in the real-time parameter;

[0151] Calculate a fourth coefficient based on the first coefficient and the third coefficient;

[0152] The real-time density of the target carbon dioxide is calculated according to the second coefficient, the fourth coefficient, the real-time temperature, and the gas constant in the real-time parameters.

[0153] In some optional implementations of some embodiments, respectively calculating the first coefficient and the second coefficient according to the gas constant, the critical temperature and the critical pressure in the real-time parameters includes:

[0154] The first coefficient is calculated according to a preset first calculation formula, and the gas constant, critical temperature and critical pressure in the real-time parameters, wherein the first calculation formula includes:

[0155]

[0156] The second coefficient is calculated according to a preset second calculation formula, and the gas constant, critical temperature and critical pressure in the real-time parameters, wherein the second calculation formula includes:

[0157]

[0158] Among them, a crepresents the first coefficient, b represents the second coefficient, R represents the gas constant in the real-time parameter, T c represents the critical temperature in real-time parameters, p c Indicates the critical pressure in real-time parameters.

[0159] In some optional implementations of some embodiments, calculating the third coefficient according to the simulated comparison temperature and the eccentricity factor in the real-time parameter includes:

[0160] The third coefficient is calculated according to a preset third calculation formula, the simulated comparison temperature, and the eccentricity factor in the real-time parameter, wherein the third calculation formula includes:

[0161]

[0162] Wherein, α represents the third coefficient, ω represents the eccentricity factor in the real-time parameter, T r represents the proposed comparison temperature.

[0163] In some optional implementations of some embodiments, calculating the real-time density of carbon dioxide according to the second coefficient, the fourth coefficient, the real-time temperature, and the gas constant in the real-time parameter includes:

[0164] The real-time density of carbon dioxide is calculated according to a preset fourth calculation formula, the second coefficient, the fourth coefficient, the real-time temperature, and the gas constant in the real-time parameter, wherein the fourth calculation formula includes:

[0165]

[0166] Wherein, p represents the real-time pressure, R represents the gas constant in the real-time parameters, T represents the real-time temperature, a represents the fourth coefficient, b represents the second coefficient, and ρ represents the real-time density of the carbon dioxide.

[0167] In some optional implementations of some embodiments, the preset flow meter is a differential pressure flow meter, and the flow parameter measured by the differential pressure flow meter is the differential pressure, wherein the mass flow rate of the target carbon dioxide is calculated according to the differential pressure in combination with the real-time density, including:

[0168] According to the preset fifth calculation formula, the mass flow rate of carbon dioxide is calculated according to the real-time density, and the fifth calculation formula includes:

[0169]

[0170] Among them, q m represents the mass flow rate, C Drepresents the outflow coefficient, D represents the throttling aperture, ε represents the expansion coefficient, ρ represents the real-time density, Δp represents the differential pressure, β represents the ratio of the throttling aperture to the inner diameter of the measuring tube, and d represents the preset calculation coefficient.

[0171] It is understandable that the modules described in the device 400 are similar to those described in the reference Figure 2 Therefore, the operations, features and beneficial effects described above for the method are also applicable to the device 400 and the modules contained therein, and will not be described in detail here.

[0172] like Figure 5 As shown, the electronic device 500 may include a processing device (e.g., a central processing unit, a graphics processing unit, etc.) 501, which can perform various appropriate actions and processes according to a program stored in a read-only memory (ROM) 502 or a program loaded from a storage device 508 into a random access memory (RAM) 503. In the RAM 503, various programs and data required for the operation of the electronic device 500 are also stored. The processing device 501, the ROM 502, and the RAM 503 are connected to each other via a bus 504. An input / output (I / O) interface 505 is also connected to the bus 504.

[0173] Typically, the following devices may be connected to the I / O interface 505: an input device 506 including, for example, a touch screen, a touch pad, a keyboard, a mouse, a camera, etc.; an output device 507 including, for example, a liquid crystal display (LCD), a speaker, etc.; a storage device 508 including, for example, a magnetic tape, a hard disk, etc.; and a communication device 509. The communication device 509 may allow the electronic device 500 to communicate with other devices wirelessly or by wire to exchange data. Figure 5 The electronic device 500 is shown with various devices, but it should be understood that it is not required to implement or possess all the devices shown. More or fewer devices may be implemented or possessed instead. Figure 5 Each block shown in the figure may represent one device, or may represent multiple devices as required.

[0174] In particular, according to some embodiments of the present disclosure, the process described above with reference to the flowchart can be implemented as a computer software program. For example, some embodiments of the present disclosure include a computer program product, which includes a computer program or instruction carried on a computer-readable medium, and the computer program or instruction contains a program code for executing the method shown in the flowchart. In some such embodiments, the computer program can be downloaded and installed from the network through the communication device 509, or installed from the storage device 508, or installed from the ROM 502. When the computer program is executed by the processing device 501, the above-mentioned functions defined in the method of some embodiments of the present disclosure are executed.

[0175] It should be noted that the computer-readable medium in some embodiments of the present disclosure may be a computer-readable signal medium or a computer-readable storage medium or any combination of the two. The computer-readable storage medium may be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, device or device, or any combination of the above. More specific examples of computer-readable storage media may 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 some embodiments of the present disclosure, the computer-readable storage medium may be any tangible medium containing or storing a program that can be used by or in combination with an instruction execution system, device or device. In some embodiments of the present disclosure, the computer-readable signal medium may include a data signal propagated in a baseband or as part of a carrier wave, in which a computer-readable program code is carried. This propagated data signal may take a variety of forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination of the above. The computer readable signal medium may also be any computer readable medium other than a computer readable storage medium, which may 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 may be transmitted using any suitable medium, including but not limited to: wires, optical cables, RF (radio frequency), etc., or any suitable combination of the above.

[0176] In some embodiments, the client and the server may communicate using any currently known or future developed network protocol such as HTTP (HyperText Transfer Protocol), and may be interconnected with any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include a local area network ("LAN"), a wide area network ("WAN"), an internet (e.g., the Internet), and a peer-to-peer network (e.g., an ad hoc peer-to-peer network), as well as any currently known or future developed network.

[0177] The computer-readable medium may be included in the device; or it may exist independently without being assembled into the electronic device. The computer-readable medium carries one or more programs. When the one or more programs are executed by the electronic device, the electronic device: obtains real-time parameters of the target carbon dioxide in liquid and / or supercritical state at real-time temperature and real-time pressure; calculates the real-time density of the target carbon dioxide at real-time temperature and real-time pressure based on the real-time parameters; calculates the mass flow rate of the target carbon dioxide based on the flow parameters measured by the preset flow meter and the real-time density.

[0178] Computer program code for performing the operations of some embodiments of the present disclosure may be written in one or more programming languages ​​or a combination thereof, including object-oriented programming languages ​​such as Java, Smalltalk, C++, and conventional procedural programming languages ​​such as "C" or similar programming languages. The program code may be executed entirely on the user's computer, partially on the user's computer, as a separate software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In cases involving a remote computer, the remote computer may be connected to the user's computer via any type of network, including a local area network (LAN) or a wide area network (WAN), or may be connected to an external computer (e.g., via the Internet using an Internet service provider).

[0179] The flow chart and block diagram in the accompanying drawings illustrate the possible architecture, function and operation of the system, method and computer program product according to various embodiments of the present disclosure. In this regard, each square box in the flow chart or block diagram can represent a module, a program segment or a part of a code, and the module, the program segment or a part of the code contains one or more executable instructions for realizing the specified logical function. It should also be noted that in some implementations as replacements, the functions marked in the square box can also occur in a sequence different from that marked in the accompanying drawings. For example, two square boxes represented in succession can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each square box in the block diagram and / or flow chart, and the combination of the square boxes in the block diagram and / or flow chart can be implemented with a dedicated hardware-based system that performs a specified function or operation, or can be implemented with a combination of dedicated hardware and computer instructions.

[0180] The modules described in some embodiments of the present disclosure may be implemented in software or hardware. The modules described may also be set in a processor, for example, may be described as:

[0181] Acquisition module, first calculation module and second calculation module. For example, the acquisition module can also be described as a "module for acquiring real-time parameters of target carbon dioxide in liquid and / or supercritical state at real-time temperature and real-time pressure".

[0182] The functions described above herein may be performed at least in part by one or more hardware logic components. For example, without limitation, exemplary types of hardware logic components that may be used include: field programmable gate arrays (FPGAs), application specific integrated circuits (ASICs), application specific standard products (ASSPs), systems on chips (SOCs), complex programmable logic devices (CPLDs), and the like.

[0183] The above descriptions are only some preferred embodiments of the present disclosure and an explanation of the technical principles used. Those skilled in the art should understand that the scope of the invention involved in the embodiments of the present disclosure is not limited to the technical solutions formed by a specific combination of the above-mentioned technical features, but should also cover other technical solutions formed by any combination of the above-mentioned technical features or their equivalent features without departing from the above-mentioned inventive concept. For example, the above-mentioned features are replaced with the technical features with similar functions disclosed in the embodiments of the present disclosure (but not limited to) and the technical solutions formed.

Claims

1. A method for determining the mass flow rate of liquid or supercritical carbon dioxide, characterized in that: include: Acquire real-time parameters of target carbon dioxide in liquid state and / or supercritical state at real-time temperature and real-time pressure; Calculate the real-time density of the target carbon dioxide at the real-time temperature and the real-time pressure based on the real-time parameters; The mass flow rate of the target carbon dioxide is calculated according to the flow parameter measured by the preset flow meter in combination with the real-time density.

2. The method according to claim 1, characterized in that The real-time parameters include a gas constant, a critical temperature, a critical pressure, an eccentricity factor, and a pseudo-comparison temperature, wherein the real-time density of the target carbon dioxide at the real-time temperature and the real-time pressure is calculated based on the real-time parameters, including: Calculate the first coefficient and the second coefficient respectively according to the gas constant, the critical temperature and the critical pressure in the real-time parameters; Calculating a simulated comparison temperature according to the real-time temperature and a critical temperature in the real-time parameters; Calculating a third coefficient according to the proposed comparison temperature and an eccentricity factor in the real-time parameter; Calculate a fourth coefficient based on the first coefficient and the third coefficient; The real-time density of the target carbon dioxide is calculated according to the second coefficient, the fourth coefficient, the real-time temperature, and the gas constant in the real-time parameters.

3. The method according to claim 2, characterized in that The calculating the first coefficient and the second coefficient respectively according to the gas constant, the critical temperature and the critical pressure in the real-time parameters comprises: The first coefficient is calculated according to a preset first calculation formula, and the gas constant, critical temperature and critical pressure in the real-time parameters, wherein the first calculation formula includes: The second coefficient is calculated according to a preset second calculation formula, and the gas constant, critical temperature and critical pressure in the real-time parameters, wherein the second calculation formula includes: Among them, a c represents the first coefficient, b represents the second coefficient, R represents the gas constant in the real-time parameter, T c represents the critical temperature in real-time parameters, p c Indicates the critical pressure in real-time parameters.

4. The method according to claim 2, characterized in that: The calculating the third coefficient according to the proposed comparison temperature and the eccentricity factor in the real-time parameter includes: The third coefficient is calculated according to a preset third calculation formula, the simulated comparison temperature, and the eccentricity factor in the real-time parameter, wherein the third calculation formula includes: Wherein, α represents the third coefficient, ω represents the eccentricity factor in the real-time parameter, T r represents the proposed comparison temperature.

5. The method according to claim 2, characterized in that: The step of calculating the real-time density of carbon dioxide according to the second coefficient, the fourth coefficient, the real-time temperature, and the gas constant in the real-time parameter comprises: The real-time density of carbon dioxide is calculated according to a preset fourth calculation formula, the second coefficient, the fourth coefficient, the real-time temperature, and the gas constant in the real-time parameter, wherein the fourth calculation formula includes: Wherein, p represents the real-time pressure, R represents the gas constant in the real-time parameters, T represents the real-time temperature, a represents the fourth coefficient, b represents the second coefficient, and ρ represents the real-time density of the carbon dioxide.

6. The method according to claim 2, characterized in that The preset flow meter is a differential pressure flow meter, and the flow parameter measured by the differential pressure flow meter is the differential pressure. The mass flow rate of the target carbon dioxide is calculated according to the differential pressure in combination with the real-time density, including: According to the preset fifth calculation formula, the mass flow rate of carbon dioxide is calculated according to the real-time density, and the fifth calculation formula includes: Among them, q m represents the mass flow rate, C D represents the outflow coefficient, D represents the throttling aperture, ε represents the expansion coefficient, ρ represents the real-time density, Δp represents the differential pressure, β represents the ratio of the throttling aperture to the inner diameter of the measuring tube, and d represents the preset calculation coefficient.

7. A device for determining the mass flow rate of liquid or supercritical carbon dioxide, characterized in that: include: An acquisition module, used for acquiring real-time parameters of target carbon dioxide in liquid state and / or supercritical state at real-time temperature and real-time pressure; A first calculation module, used for calculating the real-time density of the target carbon dioxide at the real-time temperature and the real-time pressure based on the real-time parameters; The second calculation module is used to calculate the mass flow rate of the target carbon dioxide according to the flow parameter measured by the preset flow meter in combination with the real-time density.

8. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the steps of the method according to any one of claims 1 to 6 are implemented.

9. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 6 are implemented.

10. A computer program product, comprising a computer program or instructions, characterized in that: When the computer program or instruction is executed by a processor, the steps of the method according to any one of claims 1 to 6 are implemented.

Images