Method, device and medium for estimating sound wave velocity and density of multi-component mixed gas

Abstract

The application provides a method, device and medium for estimating the sound wave speed and density of a multi-component mixed gas, which comprises the following steps: introducing a temperature-dependent function according to the Wilke gas mixing rule to obtain the absorption coefficient and repulsion coefficient of each component gas; defining the temperature-dependent function to obtain the absorption coefficient and repulsion coefficient of the mixed gas; improving the Peng-Robinson state equation based on the temperature-dependent function to obtain the density and bulk modulus of the multi-component mixed gas, and then calculating the sound wave speed of the multi-component mixed gas; and calculating the viscosity of the multi-component mixed gas according to the Lohrenz-Bray-Clark model. The Peng-Robinson state equation is improved by introducing the temperature-dependent function and the Wilke gas mixing rule to calculate the density and sound wave speed of the multi-component mixed gas under different temperature-pressure conditions, and the viscosity of the multi-component mixed gas is obtained by combining the Lohrenz-Bray-Clark model, so that the accuracy of the estimation of the density, sound wave speed and viscosity of the multi-component mixed gas is improved.

Description

Technical Field

[0001] This invention belongs to the field of physical parameter estimation technology for multi-component mixed gases, and specifically relates to a method, device and medium for estimating the acoustic velocity and density of multi-component mixed gases. Background Technology

[0002] Multi-component gas mixtures are widely found in nature and in various industrial sectors such as petroleum, natural gas, chemical production, and environmental monitoring. The velocity of sound and density of these mixtures are key physical properties that directly affect process control, equipment safety, resource exploration accuracy, and environmental gas monitoring accuracy. For example, in natural gas transportation, the velocity of sound and density can be used to calculate gas volume, determine gas quality changes, and detect leaks; in chemical reactors, gas density and velocity of sound have a significant impact on reaction efficiency and process safety.

[0003] Currently, the estimation of acoustic velocity and density of mixed gases mainly relies on empirical formulas or derivations based on the ideal gas law, such as using generalized sound velocity models or component-based weighted average methods. While these methods are applicable under normal pressure, temperature, and composition conditions, they often lead to significant deviations in estimation results under actual operating conditions with high pressure, high temperature, or complex compositions (such as those containing many non-hydrocarbon impurities). This is because they do not fully consider the effects of real gases, intermolecular interactions, and the nonlinear coupling effects of the physical properties of each component. Especially under conditions where hydrocarbon gases coexist with non-hydrocarbon components such as carbon dioxide and nitrogen, existing models have poor adaptability and are difficult to meet the needs of high-precision measurement and real-time monitoring.

[0004] Furthermore, while existing gas component detection methods such as gas chromatography and mass spectrometry offer high accuracy, they suffer from drawbacks including significant detection delays, high equipment costs, and complex maintenance, making it impossible to achieve real-time inversion and online prediction of sound wave velocity and density. This is particularly true in scenarios involving dynamically changing components, such as associated gas from oil fields and chemical tail gas, where a universal method is lacking to rapidly and accurately calculate sound wave velocity and density parameters based on limited component information.

[0005] Therefore, how to provide a method for estimating the acoustic velocity and density of multi-component mixed gases, so that it can be estimated accurately and efficiently over a wide temperature and pressure range, is a problem that urgently needs to be solved by those skilled in the art. Summary of the Invention

[0006] The purpose of this invention is to provide a method, apparatus, and medium for estimating the acoustic velocity and density of a multi-component mixed gas, so as to at least solve one of the above-mentioned technical problems.

[0007] To achieve the above objectives, the first aspect of the present invention provides a method for estimating the acoustic velocity and density of a multi-component mixed gas. The method includes: introducing a temperature-dependent function based on Wilke's gas mixing rules to obtain the absorption and repulsion coefficients of each component gas; defining the temperature-dependent function to obtain the absorption and repulsion coefficients of the mixed gas; improving the Peng-Robinson equation of state based on the temperature-dependent function to obtain the density and bulk modulus of the multi-component mixed gas; calculating the acoustic velocity of the multi-component mixed gas based on the density and bulk modulus; and calculating the viscosity of the multi-component mixed gas according to the Lohrenz-Bray-Clark model.

[0008] In the first aspect, obtaining the absorption coefficient and repulsion coefficient of each component gas includes: The formula for calculating the absorption coefficient 'a' of each component gas is: ; The formula for calculating the repulsion coefficient b of each component gas is: ; in, R is the temperature-dependent function; R is the gas constant. These are the critical temperatures of each component gas; denoted as , where is the critical pressure of each component gas.

[0009] In the first aspect, the temperature dependence function is defined as follows: ; ; in, This represents the decrease in temperature, where T is the test temperature; It is an exponential function; For the fitting function, It is the eccentricity factor.

[0010] In the first aspect, the absorption coefficient of the mixed gas The formula for calculation is: ; Repulsion coefficient of mixed gases The formula for calculation is: ; in, Let i be the mole fraction of the i-th gas. Let be the mole fraction of the j-th gas; Let be the absorption coefficient of the i-th gas; Let be the absorption coefficient of the j-th gas; The binary interaction coefficient between different gases; Let be the repulsion coefficient of the i-th gas.

[0011] In the first aspect, the improvement of the Peng-Robinson state equation includes: introducing a compression factor Z to algebraically simplify the Peng-Robinson state equation to obtain a cubic equation. ; Solving the cubic equation yields the molar volume of the gas mixture. ; The molar mass of the mixed gas is , Let be the molar mass of the i-th gas; The formula for calculating the density is: ; The formula for calculating the bulk modulus is: ; in, is the adiabatic index of the gas mixture.

[0012] In the first aspect, the formula for calculating the velocity of sound waves is: .

[0013] In the first aspect, the viscosity is calculated as follows: ; ; in, ; The critical density of a multi-component gas mixture; Let be the critical temperature of the i-th component gas. Let be the molar mass of the i-th component gas. Let be the critical pressure of the i-th component gas; This is a preliminary estimate of the relevant viscosity.

[0014] In the first aspect, the method further includes: verifying the accuracy of the density, the bulk modulus, the acoustic velocity, and the viscosity.

[0015] A second aspect of the present invention provides an electronic device comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor is configured to execute the computer program stored in the memory to implement a method for estimating the acoustic velocity and density of a multi-component mixed gas as described in any one aspect of the first aspect.

[0016] A third aspect of the present invention provides a computer-readable storage medium having a computer program stored thereon: when executed by a processor, the computer program implements a method for estimating the acoustic velocity and density of a multi-component mixed gas as described in any one of the first aspects.

[0017] Beneficial effects: This invention provides a method for estimating the acoustic velocity and density of a multi-component gas mixture. First, based on Wilke's gas mixing rules, a temperature-dependent function is introduced to obtain the absorption and repulsion coefficients of each component gas in relation to temperature changes. Second, the temperature-dependent function is defined and its expression is clarified, allowing the absorption and repulsion coefficients of the multi-component gas mixture to be calculated based on the absorption and repulsion coefficients of each component gas. Then, the absorption and repulsion coefficients of the multi-component gas mixture, along with the temperature-dependent function, are substituted into the Peng-Robinson equation of state to improve it and perform algebraic simplification, thereby obtaining the density and bulk modulus of the multi-component gas mixture. Finally, the acoustic velocity of the multi-component gas mixture is calculated based on the density and bulk modulus, and the viscosity is calculated using the Lohrenz-Bray-Clark model. This invention improves the Peng-Robinson equation of state by introducing a temperature-dependent function and Wilke's gas mixing rule to calculate the density and acoustic velocity of multi-component mixed gases under different temperature-pressure conditions. Furthermore, it utilizes the Lohrenz-Bray-Clark model to obtain the viscosity of multi-component mixed gases, thereby improving the accuracy of density, acoustic velocity, and viscosity estimations. This allows for wide application in fields such as oil and gas extraction, pipeline transportation, chemical production, environmental monitoring, safety inspection, and scientific research, making the estimation results more closely resemble real-world conditions. Attached Figure Description

[0018] To more clearly illustrate the technical solutions in the embodiments of this specification or the prior art, the drawings used in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0019] Figure 1 A schematic flowchart illustrating a method for estimating the acoustic velocity and density of a multi-component mixed gas provided by the present invention; Figure 2 A comparison chart of the bulk moduli of hydrogen, helium, methane, and carbon dioxide predicted by different methods; Figure 3 A comparison chart of the densities of hydrogen, helium, methane, and carbon dioxide predicted by different methods; Figure 4 A comparison chart of the acoustic velocities of hydrogen, helium, methane, and carbon dioxide predicted by different methods; Figure 5 A comparison chart of the viscosities of hydrogen, helium, methane, and carbon dioxide predicted by different methods; Figure 6 A schematic diagram of the structure of an electronic device provided by the present invention; Figure 7 This is a schematic diagram of the structure of a computer-readable storage medium provided by the present invention. Detailed Implementation

[0020] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative effort are within the scope of protection of the present invention.

[0021] The terms "first," "second," etc., used in the specification, claims, and accompanying drawings of this application are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that the embodiments described herein can be implemented in a sequence other than that illustrated or described herein. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover non-exclusive inclusion; for example, a process, method, system, product, or device that includes a series of steps or modules is not necessarily limited to those explicitly listed, but may include other steps or modules not explicitly listed or inherent to such processes, methods, products, or devices. The naming or numbering of steps appearing in this application does not imply that the steps in the method flow must be performed in the chronological / logical order indicated by the naming or numbering. The execution order of named or numbered process steps can be changed according to the desired technical purpose, as long as the same or similar technical effect is achieved.

[0022] The module division described in this application is a logical division. In practical applications, there may be other division methods. For example, multiple modules may be combined or integrated into another system, or some features may be ignored or not executed. In addition, the shown or discussed mutual coupling, direct coupling, or communication connection may be through some interface, and the indirect coupling or communication connection between modules may be electrical or other similar forms, none of which are limited in this application. Furthermore, the modules or sub-modules described as separate components may or may not be physically separated, may or may not be physical modules, or may be distributed in multiple circuit modules. Some or all of the modules can be selected to achieve the purpose of the solution in this application according to actual needs.

[0023] Example 1 Please see Figure 1 This embodiment provides a method for estimating the acoustic velocity and density of a multi-component gas mixture. The method includes: introducing a temperature-dependent function based on Wilke's gas mixing rule to obtain the absorption and repulsion coefficients of each component gas; defining the temperature-dependent function to obtain the absorption and repulsion coefficients of the gas mixture; improving the Peng-Robinson equation of state based on the temperature-dependent function to obtain the density and bulk modulus of the multi-component gas mixture; calculating the acoustic velocity of the multi-component gas mixture based on the density and bulk modulus; and calculating the viscosity of the multi-component gas mixture based on the Lohrenz-Bray-Clark model.

[0024] Specifically, this invention provides a method for estimating the acoustic velocity and density of a multi-component gas mixture. First, based on Wilke's gas mixing rules, a temperature-dependent function is introduced to obtain the absorption and repulsion coefficients of each component gas that are related to temperature changes. Second, the temperature-dependent function is defined and its expression is clarified, allowing the absorption and repulsion coefficients of the multi-component gas mixture to be calculated based on the absorption and repulsion coefficients of each component gas. Then, the absorption and repulsion coefficients of the multi-component gas mixture, along with the temperature-dependent function, are substituted into the Peng-Robinson equation of state to improve it and perform algebraic simplification, thereby obtaining the density and bulk modulus of the multi-component gas mixture. Finally, the acoustic velocity of the multi-component gas mixture is calculated based on the density and bulk modulus, and the viscosity is calculated using the Lohrenz-Bray-Clark model. This invention improves the Peng-Robinson equation of state by introducing a temperature-dependent function and Wilke's gas mixing rule to calculate the density and acoustic velocity of multi-component mixed gases under different temperature-pressure conditions. Furthermore, it utilizes the Lohrenz-Bray-Clark model to obtain the viscosity of multi-component mixed gases, thereby improving the accuracy of density, acoustic velocity, and viscosity estimations. This allows for wide application in fields such as oil and gas extraction, pipeline transportation, chemical production, environmental monitoring, safety inspection, and scientific research, making the estimation results more closely resemble real-world conditions.

[0025] In some possible implementations, obtaining the absorption coefficient and repulsion coefficient of each component gas includes: The formula for calculating the absorption coefficient 'a' of each component gas is: ; The formula for calculating the repulsion coefficient b of each component gas is: ; in, This is a temperature-dependent function; R is the gas constant, R = 8.31 J / (mol·K); These are the critical temperatures of each component gas; denoted as , where is the critical pressure of each component gas.

[0026] In some possible implementations, the temperature dependence function is defined as follows: ; ; in, This represents the decrease in temperature, where T is the test temperature; It is an exponential function; For the fitting function, This refers to the eccentricity factor, which varies depending on the gas. It can be obtained by consulting the Chemical Online Manual of the National Institute of Standards and Technology.

[0027] In some possible implementations, the absorption coefficient of the mixed gas The formula for calculation is: ; Repulsion coefficient of mixed gases The formula for calculation is: ; in, Let i be the mole fraction of the i-th gas. Let be the mole fraction of the j-th gas; Let be the absorption coefficient of the i-th gas; Let be the absorption coefficient of the j-th gas; The binary interaction coefficient between different gases can be determined by regression fitting of experimental data. Let be the repulsion coefficient of the i-th gas.

[0028] Those skilled in the art will understand that the temperature-dependent absorption and repulsion coefficients of the single-component gas obtained in the aforementioned steps are not sufficient for the absorption and repulsion coefficients of the mixed gas; the coupling effect between different gases also needs to be considered. Therefore, it is necessary to introduce... By combining the calculation formulas for the absorption coefficient and repulsion coefficient of a single component gas, the absorption coefficient and repulsion coefficient of a mixed gas can be obtained.

[0029] In some possible implementations, the improvement of the Peng-Robinson equation of state includes: introducing a compression factor Z to algebraically simplify the Peng-Robinson equation of state to obtain a cubic equation. ; Solving the cubic equation yields the molar volume of the gas mixture. ; The molar mass of the mixed gas is , Let be the molar mass of the i-th gas; The formula for calculating the density is: ; The formula for calculating the bulk modulus is: ; in, is the adiabatic index of the gas mixture.

[0030] Specifically, after introducing the temperature-dependent function mentioned above, the Peng-Robinson equation of state for a multi-component gas mixture is as follows: Where P is the total pressure of the gas mixture and V is the total volume of the gas mixture, and by algebraic simplification, we can obtain the compressibility factor. A cubic equation with variables is used to derive the equations for the compressibility factor Z of a real gas mixture under different conditions: By solving the above equations and selecting a suitable solution Z, the molar volume of the gas mixture can be obtained. The molar mass of the gas mixture is the sum of the molar masses of each component gas: ; in, Let be the mole fraction of the i-th gas; The molar mass of the i-th gas can be obtained by consulting the Chemical Online Manual of the National Institute of Standards and Technology.

[0031] According to the formula for calculating density Given the molar volume and molar mass of the gas mixture, its density can be obtained; and based on this, the bulk modulus of the gas mixture can be calculated. ; in, is the adiabatic index of the gas mixture.

[0032] The adiabatic index of monatomic gases (such as helium and neon) is typically 1.67, that of diatomic gases (such as nitrogen, oxygen, and hydrogen) is 1.4, and that of triatomic gases (such as carbon dioxide) is 1.3. However, polyatomic gases do not always have the same adiabatic index; it depends on the degrees of freedom of the molecules. For example, methane has an adiabatic index of 1.32. For gas mixtures, the isobaric specific volume of each component gas can be calculated first. and specific heat at constant volume : ; Then the adiabatic index of the single-component gas can be obtained. The formula for calculation is: ; Then, according to Wilke's gas mixing rules, the adiabatic index of the mixed gas... for: ; in, Let be the mole fraction of the i-th gas; Let be the isobaric specific volume of the i-th gas. Let be the specific heat at constant volume for the i-th gas.

[0033] In some possible implementations, the velocity of sound is calculated as follows: .

[0034] Those skilled in the art will understand that, given the adiabatic index of a gas mixture, its bulk modulus can be calculated. Combined with the density of the gas mixture, the longitudinal wave velocity can then be obtained. In deep underground fluids, the longitudinal wave velocity represents the sound wave velocity; that is, in deep underground flow, the sound wave velocity of the gas mixture is... .

[0035] Furthermore, according to the Lohrenz-Bray-Clark model, the viscosity of the gas mixture can be expressed as: ; ; in, ; The critical density of a multi-component gas mixture; Let be the critical temperature of the i-th component gas. Let be the molar mass of the i-th component gas. Let be the critical pressure of the i-th component gas; For the preliminary estimate of the relevant viscosity; It can be expressed as a simplified density. ; in, Let be the critical volume of the i-th gas; ; According to The range of values ​​depends on the choice of different formulas for calculation. This makes the viscosity calculation of the mixed gas more accurate.

[0036] In some possible implementations, the method further includes verifying the accuracy of the density, the bulk modulus, the acoustic velocity, and the viscosity.

[0037] In this application, taking a mixed gas comprising hydrogen, helium, methane, and carbon dioxide as an example, the density, bulk modulus, sound velocity, and viscosity of the mixed gas are calculated using both conventional methods and the calculation method provided above in this invention, and compared with data from the NIST (National Institute of Standards and Technology, Chemical Handbook) database. Figure 2-5As shown in the figure, the experimental data indicates that the data obtained by using the method provided by the present invention is closer to the NIST data. That is, the calculation method provided by the present invention is closer to the real underground state and improves the accuracy of underground state prediction.

[0038] Figure 2 A comparison chart of the bulk moduli of hydrogen, helium, methane, and carbon dioxide predicted by different methods; Figure 3 A comparison chart of the densities of hydrogen, helium, methane, and carbon dioxide predicted by different methods; Figure 4 A comparison chart of the acoustic velocities of hydrogen, helium, methane, and carbon dioxide predicted by different methods; Figure 5 A comparison chart of the viscosities of hydrogen, helium, methane, and carbon dioxide predicted by different methods; in, Figure 2-5 In the figure, the red solid lines represent predictions made using the calculation method of this application, the black solid lines represent predictions made using traditional methods, and the blue dots represent data from the NIST database. As can be seen from the figure, there are some discrepancies between the data calculated by the traditional method and the NIST data, while the data of this application has a high degree of agreement with the NIST data. That is, this application provides more accurate predictions and can more accurately simulate the compression characteristics of mixed gases under high pressure, making it applicable to the simulation of states at deeper underground locations, thereby providing a theoretical basis for the development of resources at deeper underground locations.

[0039] In summary, this invention comprehensively considers the elastic parameters and viscosity changes of different component gases under different temperatures and pressures. It does not treat the mixed gas as an ideal gas, but rather performs a simple linear weighting and improves the traditional Peng-Robinson equation and Lohrenz-Bray-Clark model to make them closer to the real state of multi-component mixed gases in underground reservoirs, thereby enabling more accurate prediction of underground conditions.

[0040] Example 2 Please see Figure 6 , Figure 6 This is a schematic diagram of one embodiment of the electronic device 100 of the present invention, including: The system includes a memory 101, a processor 102, and a computer program 103 stored in the memory and executable on the processor. When the processor executes the computer program 103 stored in the memory, it implements the above-described method for estimating the acoustic velocity and density of a multi-component mixed gas.

[0041] For ease of explanation, only the parts related to the embodiments of the present invention are shown. For specific technical details not disclosed, please refer to the part on the method for estimating the acoustic velocity and density of multi-component mixed gases in the embodiments of the present invention. The memory 101 can be used to store the computer program 103, which includes software programs, modules, and data. The processor 102 executes the computer program 103 stored in the memory 101 to perform various functional applications and data processing of the electronic device.

[0042] Example 3 This invention also provides a computer-readable storage medium; please refer to [link to relevant documentation]. Figure 7 , Figure 7 This is a schematic diagram of an embodiment of a computer-readable storage medium in the present invention. The computer-readable storage medium may store a computer program, which, when executed, includes some or all of the steps of the method for estimating the acoustic velocity and density of a multi-component mixed gas described in the above method embodiments.

[0043] Those skilled in the art will readily understand that, for the sake of convenience and brevity, the specific working processes of the devices, electronic equipment, and computer-readable storage media described above can be referred to the corresponding processes of the sound wave velocity and density estimation methods for multi-component mixed gases in the foregoing method embodiments, and will not be repeated here.

[0044] In the embodiments provided in this application, it should be understood that the disclosed systems and methods can be implemented in other ways. For example, the system embodiments described above are merely illustrative; for instance, the division of units is only a logical functional division, and in actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Furthermore, the mutual coupling or direct coupling or communication connection shown or discussed may be through some interfaces; the indirect coupling or communication connection between systems or units may be electrical, mechanical, or other forms.

[0045] The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the units can be selected to achieve the purpose of this embodiment according to actual needs.

[0046] Furthermore, the functional units in the various embodiments of the present invention can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit. The integrated unit can be implemented in hardware or as a software functional unit.

[0047] If the integrated unit is implemented as a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or all or part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the method for estimating the acoustic velocity and density of multi-component mixed gases according to various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0048] The above-described embodiments are only used to illustrate the technical solutions of the present invention, and are not intended to limit it. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features. Such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for estimating the acoustic velocity and density of a multi-component mixed gas, characterized in that, The method includes: Based on Wilke's gas mixing rules, a temperature-dependent function is introduced to obtain the absorption and repulsion coefficients of each component gas. The temperature dependence function is defined to obtain the absorption coefficient and repulsion coefficient of the mixed gas; Based on the temperature dependence function, the Peng-Robinson equation of state is improved to obtain the density and bulk modulus of the multi-component mixed gas; Calculate the acoustic velocity of the multi-component mixed gas based on the density and the bulk modulus. The viscosity of the multi-component gas mixture was calculated based on the Lohrenz-Bray-Clark model. The process of obtaining the absorption coefficient and repulsion coefficient of each component gas includes: The formula for calculating the absorption coefficient 'a' of each component gas is: The formula for calculating the repulsion coefficient b of each component gas is: in, R is the temperature-dependent function; R is the gas constant. The critical temperature of each component gas; The critical pressure of each component gas; The temperature dependence function is defined as follows: ； ； in, This represents the decrease in temperature, where T is the test temperature; It is an exponential function; For the fitting function, It is the eccentricity factor; Absorption coefficient of mixed gas The formula for calculation is: Repulsion coefficient of mixed gases The formula for calculation is: in, Let i be the mole fraction of the i-th gas. Let be the mole fraction of the j-th gas; Let be the absorption coefficient of the i-th gas; Let be the absorption coefficient of the j-th gas; The binary interaction coefficient between different gases; Let be the repulsion coefficient of the i-th gas; The improvement to the Peng-Robinson state equations includes: introducing a compression factor Z, and performing algebraic simplification on the Peng-Robinson state equations to obtain cubic equations. ； Solving the cubic equation yields the molar volume of the gas mixture. ; The molar mass of the mixed gas is , Let be the molar mass of the i-th gas; The formula for calculating the density is: ； The formula for calculating the bulk modulus is: ； in, is the adiabatic index of the gas mixture.

2. The method for estimating the acoustic velocity and density of a multi-component mixed gas according to claim 1, characterized in that, The formula for calculating the velocity of sound waves is: .

3. The method for estimating the acoustic velocity and density of a multi-component mixed gas according to claim 2, characterized in that, The formula for calculating the viscosity is: ； ； in, ; The critical density of a multi-component gas mixture; Let be the critical temperature of the i-th component gas. Let be the molar mass of the i-th component gas. Let be the critical pressure of the i-th component gas; To and Preliminary estimates of the relevant viscosity.

4. The method for estimating the acoustic velocity and density of a multi-component mixed gas according to claim 3, characterized in that, The method further includes verifying the accuracy of the density, the bulk modulus, the acoustic velocity, and the viscosity.

5. An electronic device, comprising: A memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that the processor, when executing the computer program stored in the memory, implements a method for estimating the acoustic velocity and density of a multi-component mixed gas as described in any one of claims 1-4.

6. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by the processor, it implements the method for estimating the acoustic velocity and density of a multi-component mixed gas as described in any one of claims 1-4.

Images