Elemental sulfur solubility prediction method, system, electronic device, and storage medium

Abstract

This invention discloses a method for predicting the solubility of elemental sulfur, comprising the following steps: obtaining parameter information of the target gas and calculating the target gas density ρ; selecting a prediction model based on the target gas temperature T, and determining the model parameters based on the H2S mole fraction and the target gas density ρ: when 373.15K≤T≤433.15K, the prediction model is c=ρ k exp(a / T+b); When 303.15K≤T<373.15K, the prediction model is: where c is the solubility of elemental sulfur in the target gas, k, a, and b are model parameters, model parameter k is a function of T, and the low-temperature correction coefficient d; Based on the selected prediction model and the determined model parameters, the solubility of elemental sulfur in the target gas is calculated. This invention also discloses an elemental sulfur solubility prediction system, electronic equipment, and storage medium. The elemental sulfur solubility prediction method of this invention is applicable to the temperature span in the development process of acidic gas reservoirs, introducing a low-temperature coefficient in the low-temperature range to ensure prediction accuracy throughout the entire temperature range.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of acid gas reservoir development and sulfur deposition, and particularly relates to an elemental sulfur solubility prediction method and system, an electronic device and a storage medium. BACKGROUND

[0002] Elemental sulfur deposition often occurs in acid gas reservoir development, which causes wellbore or gathering pipeline plugging and aggravates corrosion of the pipe material. Severe sulfur plugging will cause pipeline pressure to be blocked, deformed, ruptured and corroded, which may lead to leakage of H2S-containing acid gas and cause serious production safety problems. The essence of sulfur deposition is that sulfur content exceeds its solubility in acid natural gas and is precipitated, so the prerequisite for predicting the sulfur deposition site and degree is to obtain the solubility of elemental sulfur in acid gas, which requires a method for accurately predicting the solubility of sulfur.

[0003] In 1982, Chrastil derived a general model for the solubility of solid or liquid solutes in a gas. Subsequently, many scholars have conducted research on the prediction of elemental sulfur solubility in acid gas based on the Chrastil model. Roberts first applied the Chrastil model to fit sulfur solubility experimental data in 1997, obtained an empirical formula for calculating sulfur solubility, and achieved certain results when applied to the Waterton gas reservoir conditions in Canada. Since the data used to fit the model parameters are too few, the application conditions of Roberts' empirical formula are limited, and when the actual conditions deviate too much from the Waterton gas reservoir, the prediction accuracy decreases significantly.

[0004] In 2006, Qiao Haibo et al. pointed out in the paper "Prediction Model of Elemental Sulfur Solubility in Sulfur-Containing Gas" that when the temperature and pressure range is too large, the linear relationship between the logarithm of sulfur solubility and the logarithm of acid gas density deviates. Therefore, the idea of setting a density inflection point for segmented fitting is proposed, that is, when applying the Chrastil model to fit experimental data, the acid gas density is divided into high-density and low-density regions, and the experimental data in the two regions are fitted respectively to obtain two sets of model parameter values, which improves the prediction accuracy of the model to some extent, but does not consider the influence of temperature and gas composition on the model parameter values.

[0005] Patent document CN102998422A discloses a method for predicting the solubility of elemental sulfur in sulfur-containing natural gas, sets a density inflection point, and considers the influence of H2S content on the Chrastil model parameter values. Under different H2S contents, two sets of model parameter values are obtained for predicting the elemental sulfur solubility in the high / low density region, but the influence of temperature on the model parameter values is not considered.

[0006] As acidic gases flow from underground along well shafts and gathering pipelines, the temperature range spans from tens to over one hundred degrees Celsius. However, most studies on sulfur solubility prediction only cover a relatively narrow temperature range. Therefore, there is a need for a method to predict elemental sulfur solubility that considers the influence of multiple factors on the parameters of the Chrastil model, is applicable to a wide temperature range, and maintains high prediction accuracy.

[0007] The information disclosed in this background section is intended only to enhance the understanding of the overall background of the invention and should not be construed as an admission or in any way implying that the information constitutes prior art known to those skilled in the art. Summary of the Invention

[0008] One of the objectives of this invention is to provide a method, system, electronic device, and storage medium for predicting the solubility of elemental sulfur, thereby improving the problem that existing models for predicting the solubility of elemental sulfur are only applicable within a narrow temperature range.

[0009] Another objective of this invention is to provide a method, system, electronic device, and storage medium for predicting the solubility of elemental sulfur, thereby ensuring the accuracy of the prediction of the solubility of elemental sulfur over a wide temperature range.

[0010] To achieve the above objectives, according to a first aspect of the present invention, the present invention provides a method for predicting the solubility of elemental sulfur, comprising the following steps:

[0011] Obtain the parameter information of the target gas and calculate the density ρ of the target gas;

[0012] Based on the target gas temperature T, a prediction model is selected, and based on the H2S mole fraction y H2S Given the target gas density ρ, determine the model parameters:

[0013] When 373.15K≤T≤433.15K, the prediction model is c=ρ k exp(a / T+b);

[0014] When 303.15K ≤ T < 373.15K, the prediction model is:

[0015] In the formula, c represents the solubility of sulfur in the target gas, in g / m³. 3 k, a, and b are model parameters, where model parameter k is a function of T, and the low-temperature correction coefficient d = 0.0079264·T-2.2623;

[0016] Based on the selected prediction model and the determined model parameters, the solubility of elemental sulfur in the target gas is calculated.

[0017] Furthermore, in the above technical solution, when the prediction model is When determining model parameters, the following are included:

[0018] when When the model parameters are

[0019] k=-0.00558·T+3.8555, a=-9425.1, b=14.654;

[0020] when When the model parameters are

[0021] k=-0.0016125·T+1.9862, a=-5973.8, b=7.0301;

[0022] when When the model parameters are

[0023] k=-0.004768·T+3.2202, a=-9150.6, b=15.231;

[0024] when When the model parameters are

[0025] k=-0.003158·T+2.636, a=-8453.2, b=12.89.

[0026] Furthermore, in the above technical solution, when the prediction model is c = ρ k When exp(a / T+b), based on the mole fraction of H2S Given the target gas density ρ, the model parameters are determined, including those based on the H2S mole fraction. Calculate the density inflection point ρ', in kg / m³. 3 ; and compare the target gas density ρ and the density inflection point ρ', and determine the model parameters based on the comparison results.

[0027] Furthermore, in the above technical solution,

[0028] when At that time, the density inflection point ρ' = 0.0029063·T was calculated. 2 -3.2056·T+1064.7,

[0029] If ρ < ρ', then the model parameters are:

[0030] k=-0.00558·T+3.8555, a=-9425.1, b=14.654,

[0031] If ρ≥ρ', then the model parameters are:

[0032] k=-0.02158·T+12.535, a=-24839, b=40.732;

[0033] when At that time, the density inflection point ρ' = 0.0029312·T 2 -3.1846·T+1052.9,

[0034] If ρ < ρ', then the model parameters are:

[0035] k=-0.0016125·T+1.9862, a=-5973.8, b=7.0301,

[0036] If ρ≥ρ', then the model parameters are:

[0037] k=-0.017925·T+10.571, a=-21886, b=35.535;

[0038] when At that time, the density inflection point ρ' = 0.00235·T 2 -2.5418·T+845.11,

[0039] If ρ < ρ', then the model parameters are:

[0040] k=-0.029008·T+13.651, a=-30849, b=65.672,

[0041] If ρ≥ρ', then the model parameters are:

[0042] k=-0.02987·T+15.652, a=-32563, b=61.158;

[0043] when At that time, the density inflection point ρ' = 0.0023812·T 2 -2.5722·T+856.21,

[0044] If ρ < ρ', then the model parameters are:

[0045] k=-0.035523·T+16.397, a=-36200, b=78.168,

[0046] If ρ≥ρ', then the model parameters are:

[0047] k=-0.012888·T+8.4103, a=-17570, b=25.81.

[0048] Furthermore, in the above technical solution, the parameter information of the target gas includes the target gas temperature T, the target gas pressure p, the acidic gas components in the target gas, and the mole fraction of each component.

[0049] Furthermore, in the above technical solution, calculating the target gas density ρ includes: calculating the apparent critical pressure p based on the parameter information of the target gas. pc Apparent critical temperature T pc and apparent relative molecular mass M g According to the mole fraction of H2S and CO2 mole fraction For apparent critical pressure p pc and apparent critical temperature T pc Perform correction; based on the target gas parameters, the corrected apparent critical pressure, and the apparent critical temperature, calculate the apparent corresponding pressure p. pr and the corresponding apparent temperature T pr According to the apparent pressure p pr and the corresponding apparent temperature T pr Consult the natural gas compressibility factor chart to obtain the corresponding compressibility factor value Z for the target gas; and use the target gas pressure p, target gas temperature T, and apparent relative molecular mass M. g Calculate the target gas density based on the corresponding compressibility factor value Z.

[0050]

[0051] In the formula, R = 8.314 J / (mol·K) is the universal gas constant.

[0052] Furthermore, in the above technical solution, the apparent critical pressure p pc Apparent critical temperature T pc and apparent relative molecular mass M g The calculation formulas are as follows:

[0053]

[0054] In the formula, p ci T represents the critical pressure of component i, in MPa. ci M is the critical temperature of component i, in K; i y is the relative molecular mass of component i; i Let be the mole fraction of component i.

[0055] Furthermore, in the above technical solution, the corrected apparent critical pressure p pc 'and the corrected apparent critical temperature T pc 'They are respectively:

[0056]

[0057]

[0058] In the formula, ε is the correction coefficient; A is the sum of the mole fractions of H2S and CO2; and B is the mole fraction of H2S.

[0059] Furthermore, in the above technical solution, depending on the corresponding pressure p pr and the corresponding apparent temperature T pr The calculation formulas are as follows:

[0060]

[0061] Furthermore, in the above technical solution, the pressure range of the target gas is 10-60 MPa, the temperature range is 303.15 K-433.15 K, the mole fraction of H2S is 1%-25%, and the mole fraction of CO2 is 0-25%.

[0062] According to a second aspect of the present invention, the present invention provides an elemental sulfur solubility prediction system, comprising:

[0063] The data acquisition unit is used to acquire parameter information of the target gas;

[0064] A calculation unit, used to calculate the target gas density ρ;

[0065] The model acquisition unit is used to select a prediction model based on the target gas temperature T and the H2S mole fraction. Given the target gas density ρ, determine the model parameters:

[0066] When 373.15K≤T≤433.15K, the prediction model is c=ρ k exp(a / T+b);

[0067] When 303.15K ≤ T < 373.15K, the prediction model is: In the formula, c represents the solubility of sulfur in the target gas, in g / m³. 3 k, a, and b are model parameters, where model parameter k is a function of T, and the low-temperature correction coefficient d = 0.0079264·T⁻².2623; and

[0068] The prediction unit calculates the solubility of elemental sulfur in the target gas based on the selected prediction model and the determined model parameters.

[0069] According to a third aspect of the present invention, an electronic device is provided, comprising: at least one processor; and a memory communicatively connected to the at least one processor; wherein the memory stores instructions executable by the at least one processor, the instructions being executed by the at least one processor to cause the at least one processor to perform the elemental sulfur solubility prediction method as described in any of the above technical solutions.

[0070] According to a fourth aspect of the present invention, the present invention provides a non-transitory computer-readable storage medium storing computer-executable instructions for causing a computer to execute the elemental sulfur solubility prediction method of any of the above-described technical solutions.

[0071] Compared with the prior art, the present invention has one or more of the following beneficial effects:

[0072] 1. The elemental sulfur solubility prediction method of the present invention has a wide applicable temperature range (303.15K~433.15K), which can cover most of the conventional temperature range in the development of acidic gas reservoirs. By introducing a low temperature coefficient in the low temperature range, high prediction accuracy is guaranteed throughout the entire temperature range.

[0073] 2. The elemental sulfur solubility prediction method of the present invention comprehensively considers the influence of temperature, acid gas density and gas composition on model parameters, and divides multiple condition intervals to further improve prediction accuracy.

[0074] 3. This invention establishes a functional relationship between model parameter k and T, which makes the prediction model continuous across various temperature ranges and can predict the solubility of elemental sulfur at any temperature within the range.

[0075] 4. This invention establishes a functional relationship between density inflection point and temperature, which is more consistent with experimental laws and provides more accurate prediction results.

[0076] 5. The elemental sulfur solubility prediction method of the present invention has a wide range of applications and can predict the elemental sulfur solubility of target gases with a pressure range of 10 to 60 MPa, a temperature range of 303.15 K to 433.15 K, an H2S mole fraction range of 1% to 25%, and a CO2 mole fraction range of 0% to 25%. The process is simple and the results are highly accurate.

[0077] The above description is merely an overview of the technical solution of the present invention. In order to better understand the technical means of the present invention and to implement it according to the contents of the specification, and to make the above and other objects, technical features and advantages of the present invention easier to understand, one or more preferred embodiments are listed below and described in detail with reference to the accompanying drawings. Attached Figure Description

[0078] Figure 1 This is a flowchart illustrating a method for predicting the solubility of elemental sulfur according to an embodiment of the present invention.

[0079] Figure 2 This is a schematic diagram of the structure of an elemental sulfur solubility prediction system according to an embodiment of the present invention.

[0080] Figure 3This is a schematic diagram of the hardware structure of an electronic device for performing an elemental sulfur solubility prediction method according to an embodiment of the present invention. Detailed Implementation

[0081] The specific embodiments of the present invention will now be described in detail with reference to the accompanying drawings, but it should be understood that the scope of protection of the present invention is not limited to the specific embodiments.

[0082] Unless otherwise expressly stated, throughout the specification and claims, the term "comprising" or its variations such as "including" or "comprises" shall be understood to include the stated elements or components without excluding other elements or other components.

[0083] In this document, for ease of description, spatial relative terms such as “below,” “under,” “down,” “above,” “above,” “upper,” etc., are used to describe the relationship of one element or feature to another element or feature in the accompanying drawings. It should be understood that spatial relative terms are intended to encompass different orientations of an object in use or operation, in addition to those depicted in the figures. For example, if an object in the figure is flipped, an element described as “below” or “under” another element or feature would be oriented “above” that element or feature. Thus, the exemplary term “below” can encompass both the downward and upward orientations. An object may also have other orientations (rotated 90 degrees or other orientations), and the spatial relative terms used herein should be interpreted accordingly.

[0084] In this document, the terms "first," "second," etc., are used to distinguish two different elements or parts, and are not used to define specific positions or relative relationships. In other words, in some embodiments, the terms "first," "second," etc., can also be used interchangeably.

[0085] like Figure 1 As shown, the elemental sulfur solubility prediction method according to a specific embodiment of the present invention includes the following steps:

[0086] S100 acquires the parameter information of the target gas and calculates the target gas density ρ.

[0087] S200 selects a prediction model based on the target gas temperature T and the H2S mole fraction. Given the target gas density ρ, determine the model parameters:

[0088] When 373.15K≤T≤433.15K, the prediction model is:

[0089] c = ρ k exp(a / T+b),

[0090] In the formula, c represents the solubility of sulfur in the target gas, in g / m³.3 k, a, and b are model parameters, and the model parameter k is a function of T.

[0091] S211 when When calculating the density inflection point (unit: kg / m³), 3 ):

[0092] ρ' = 0.0029063·T 2 -3.2056·T+1064.7,

[0093] If ρ < ρ', then the model parameters are:

[0094] k=-0.00558·T+3.8555, a=-9425.1, b=14.654,

[0095] If ρ≥ρ', then the model parameters are:

[0096] k=-0.02158·T+12.535, a=-24839, b=40.732;

[0097] S212 when At that time, the density inflection point (unit: kg / m³) 3 ):

[0098] ρ' = 0.0029312·T 2 -3.1846·T+1052.9,

[0099] If ρ < ρ', then the model parameters are:

[0100] k=-0.0016125·T+1.9862, a=-5973.8, b=7.0301,

[0101] If ρ≥ρ', then the model parameters are:

[0102] k=-0.017925·T+10.571, a=-21886, b=35.535;

[0103] S213 when At that time, the density inflection point (unit: kg / m³) 3 ):

[0104] ρ' = 0.00235·T 2 -2.5418·T+845.11,

[0105] If ρ < ρ', then the model parameters are:

[0106] k=-0.029008·T+13.651, a=-30849, b=65.672,

[0107] If ρ≥ρ', then the model parameters are:

[0108] k=-0.02987·T+15.652, a=-32563, b=61.158;

[0109] S214 when At that time, the density inflection point (unit: kg / m³) 3 ):

[0110] ρ' = 0.0023812·T 2 -2.5722·T+856.21,

[0111] If ρ < ρ', then the model parameters are:

[0112] k=-0.035523·T+16.397, a=-36200, b=78.168,

[0113] If ρ≥ρ', then the model parameters are:

[0114] k=-0.012888·T+8.4103, a=-17570, b=25.81.

[0115] When 303.15K ≤ T < 373.15K, the prediction model is:

[0116]

[0117] d = 0.0079264·T - 2.2623,

[0118] In the formula, c represents the solubility of sulfur in the target gas, in g / m³. 3 k, a, and b are model parameters, where k is a function of T and d is the low-temperature correction coefficient.

[0119] S221 when At that time, the model parameters are:

[0120] k=-0.00558·T+3.8555, a=-9425.1, b=14.654.

[0121] S222 when When the model parameters are

[0122] k=-0.0016125·T+1.9862, a=-5973.8, b=7.0301.

[0123] S223 when When the model parameters are

[0124] k=-0.004768·T+3.2202, a=-9150.6, b=15.231.

[0125] S224 when When the model parameters are

[0126] k=-0.003158·T+2.636, a=-8453.2, b=12.89.

[0127] S300 calculates the solubility of elemental sulfur in the target gas based on the selected prediction model and the determined model parameters.

[0128] Furthermore, in one or more exemplary embodiments of the present invention, the parameter information of the target gas includes the target gas temperature T, the target gas pressure p, the acidic gas components in the target gas, and the mole fraction of each acidic gas component.

[0129] Furthermore, in one or more exemplary embodiments of the present invention, calculating the target gas density ρ may include the following steps:

[0130] S110 calculates the apparent critical pressure p based on the parameter information of the target gas. pc Apparent critical temperature T pc and apparent relative molecular mass M g The calculation formulas are as follows:

[0131]

[0132] In the formula, p ci T represents the critical pressure of component i, in MPa. ci M is the critical temperature of component i, in K; i y is the relative molecular mass of component i; i Let represent the mole fraction of component i. The critical pressure and critical temperature of each component can be obtained from existing data.

[0133] S120 is based on the mole fraction of H2S and CO2 mole fraction For apparent critical pressure p pc and apparent critical temperature T pc Perform correction. Obtain the corrected apparent critical pressure p. pc 'and the corrected apparent critical temperature T pc 'They are respectively:

[0134]

[0135]

[0136] In the formula, ε is the correction coefficient; A is the sum of the mole fractions of H2S and CO2; and B is the mole fraction of H2S.

[0137] S130 calculates the apparent pressure p based on the target gas parameters, the corrected apparent critical pressure, and the apparent critical temperature. pr and the corresponding apparent temperature T pr :

[0138]

[0139] S140 according to the apparent corresponding pressure p pr and the corresponding apparent temperature T pr Query the natural gas compressibility factor chart to obtain the corresponding compressibility factor value Z of the target gas.

[0140] S150 uses target gas pressure p, target gas temperature T, and apparent relative molecular mass M. g Calculate the target gas density (in kg / m³) based on the corresponding compressibility factor value Z. 3 ):

[0141]

[0142] In the formula, R = 8.314 J / (mol·K) is the universal gas constant.

[0143] Furthermore, in one or more exemplary embodiments of the present invention, the target gas has a pressure range of 10 to 60 MPa, a temperature range of 303.15 K to 433.15 K, an H2S mole fraction range of 1% to 25%, and a CO2 mole fraction range of 0% to 25%.

[0144] The method, system, electronic device and storage medium for predicting the solubility of elemental sulfur of the present invention are described in more detail below by way of specific embodiments. It should be understood that the embodiments are merely exemplary and the present invention is not limited thereto.

[0145] Example 1

[0146] refer to Figure 2 As shown, the elemental sulfur solubility prediction system includes: a data acquisition unit 10, which acquires parameter information of the target gas; a calculation unit 20, which calculates the target gas density ρ; and a model acquisition unit 30, which selects a prediction model based on the target gas temperature T and the H2S mole fraction. Given the target gas density ρ, determine the model parameters:

[0147] When 373.15K≤T≤433.15K, the prediction model is c=ρ k exp(a / T+b);

[0148] When 303.15K ≤ T < 373.15K, the prediction model is:

[0149] In the formula, c represents the solubility of sulfur in the target gas, in g / m³. 3 k, a, and b are model parameters, where model parameter k is a function of T, and the low-temperature correction coefficient d = 0.0079264·T-2.2623; and prediction unit 40, which calculates the elemental sulfur solubility in the target gas based on the selected prediction model and the determined model parameters.

[0150] Example 2

[0151] The elemental sulfur solubility prediction method according to the present invention was used to predict the elemental sulfur solubility of the target gas. In this embodiment, the composition (mole fraction) of the target gas is 20% H2S, 10% CO2, 66% CH4 and 4% N2, and the temperature and pressure of the target gas are shown in Table 1.

[0152] According to the elemental sulfur solubility prediction method of the present invention, the target gas density is calculated according to steps S110 to S150: based on the parameter information of the target gas, the apparent critical pressure p is calculated. pc Apparent critical temperature T pc and apparent relative molecular mass M g According to the mole fraction of H2S and CO2 mole fraction For apparent critical pressure p pc and apparent critical temperature T pc Perform correction; based on the target gas parameters, the corrected apparent critical pressure, and the apparent critical temperature, calculate the apparent corresponding pressure p. pr and the corresponding apparent temperature T pr According to the apparent pressure p pr and the corresponding apparent temperature T pr Consult the natural gas compressibility factor chart to obtain the corresponding compressibility factor value Z for the target gas; and use the target gas pressure p, target gas temperature T, and apparent relative molecular mass M. g Calculate the target gas density based on the corresponding compressibility factor value Z.

[0153]

[0154] In the formula, R = 8.314 J / (mol·K) is the universal gas constant. The calculated target gas density is shown in Table 1.

[0155] In this embodiment, the target gas temperature is 373.15K ≤ T ≤ 433.15K, therefore the prediction model is chosen as c = ρ k exp(a / T+b), mole fraction of H2S The density is 20%. According to step S211, the density inflection point is calculated:

[0156] ρ' = 0.0029063·T 2 -3.2056·T+1064.7, compare the target gas density ρ and the density inflection point ρ'. If ρ < ρ', the model parameters are: k = -0.00558·T+3.8555, a = -9425.1, b = 14.654; if ρ ≥ ρ', the model parameters are: k = -0.02158·T+12.535, a = -24839, b = 40.732.

[0157] Based on the selected prediction model and the determined model parameters, the solubility of elemental sulfur in the target gas is calculated, the predicted value of elemental sulfur solubility is obtained, and the relative error is calculated by comparing it with the measured value, as shown in Table 1.

[0158] Table 1. Experimental conditions and comparison of measured and predicted values ​​of elemental sulfur solubility.

[0159]

[0160]

[0161] As shown in Table 1, the predicted values ​​of elemental sulfur solubility according to the present invention are very close to the measured values, with an average relative error of only 5.51%.

[0162] Example 3

[0163] The elemental sulfur solubility prediction method according to the present invention was used to predict the elemental sulfur solubility of the target gas. In this embodiment, the composition (molar fraction) of the target gas is 7% H2S, 20% CO2, 65% CH4 and 8% N2, and the temperature and pressure of the target gas are shown in Table 2.

[0164] According to the elemental sulfur solubility prediction method of the present invention, the target gas density is calculated according to steps S110 to S150, and the calculated target gas density is shown in Table 2.

[0165] In this embodiment, the target gas temperature is 373.15K ≤ T ≤ 433.15K, therefore the prediction model is chosen as c = ρ k exp(a / T+b), mole fraction of H2S The value is 7%. According to step S212, the density inflection point is calculated as: ρ' = 0.0029312·T 2-3.1846·T+1052.9, compare the target gas density ρ and the density inflection point ρ'. If ρ < ρ', the model parameters are: k = -0.0016125·T+1.9862, a = -5973.8, b = 7.0301; if ρ ≥ ρ', the model parameters are: k = -0.017925·T+10.571, a = -21886, b = 35.535.

[0166] Based on the selected prediction model and the determined model parameters, the solubility of elemental sulfur in the target gas is calculated, the predicted value of elemental sulfur solubility is obtained, and the relative error is calculated by comparing it with the measured value, as shown in Table 2.

[0167] Table 2 Comparison of experimental conditions and measured and predicted values ​​of elemental sulfur solubility.

[0168]

[0169]

[0170] As shown in Table 2, the predicted values ​​of elemental sulfur solubility according to the present invention are very close to the measured values, with an average relative error of only 5.56%.

[0171] Example 4

[0172] The elemental sulfur solubility prediction method according to the present invention was used to predict the elemental sulfur solubility of the target gas. In this embodiment, the composition (mole fraction) of the target gas is 6% H2S, 9% CO2, 81% CH4 and 4% N2, and the temperature and pressure of the target gas are shown in Table 3.

[0173] According to the elemental sulfur solubility prediction method of the present invention, the target gas density is calculated according to steps S110 to S150, and the calculated target gas density is shown in Table 3.

[0174] In this embodiment, the target gas temperature is 373.15K ≤ T ≤ 433.15K, therefore the prediction model is chosen as c = ρ k exp(a / T+b), mole fraction of H2S The value is 6%. According to step S213, the density inflection point is calculated as: ρ' = 0.00235·T 2 -2.5418·T+845.11, compare the target gas density ρ and the density inflection point ρ'. If ρ < ρ', the model parameters are: k = -0.029008·T+13.651, a = -30849, b = 65.672; if ρ ≥ ρ', the model parameters are: k = -0.02987·T+15.652, a = -32563, b = 61.158.

[0175] Based on the selected prediction model and the determined model parameters, the solubility of elemental sulfur in the target gas is calculated, the predicted value of elemental sulfur solubility is obtained, and the relative error is calculated by comparing it with the measured value, as shown in Table 3.

[0176] Table 3 Comparison of experimental conditions and measured and predicted values ​​of elemental sulfur solubility.

[0177]

[0178]

[0179] As shown in Table 3, the predicted values ​​of elemental sulfur solubility according to the present invention are very close to the measured values, with an average relative error of only 6.13%.

[0180] Example 5

[0181] The elemental sulfur solubility prediction method according to the present invention was used to predict the elemental sulfur solubility of the target gas. In this embodiment, the composition (mole fraction) of the target gas is 1% H2S, 14% CO2, 81% CH4 and 4% N2, and the temperature and pressure of the target gas are shown in Table 4.

[0182] According to the elemental sulfur solubility prediction method of the present invention, the target gas density is calculated according to steps S110 to S150, and the calculated target gas density is shown in Table 4.

[0183] In this embodiment, the target gas temperature is 373.15K ≤ T ≤ 433.15K, therefore the prediction model is chosen as c = ρ k exp(a / T+b), mole fraction of H2S Given a density of 1%, according to step S214, the density inflection point is calculated as: ρ' = 0.0023812·T 2 -2.5722·T+856.21, compare the target gas density ρ and the density inflection point ρ'. If ρ < ρ', the model parameters are: k = -0.035523·T+16.397, a = -36200, b = 78.168; if ρ ≥ ρ', the model parameters are: k = -0.012888·T+8.4103, a = -17570, b = 25.81.

[0184] Based on the selected prediction model and the determined model parameters, the solubility of elemental sulfur in the target gas is calculated, the predicted value of elemental sulfur solubility is obtained, and the relative error is calculated by comparing it with the measured value, as shown in Table 4.

[0185] Table 4. Comparison of experimental conditions and measured and predicted values ​​of elemental sulfur solubility.

[0186]

[0187] As shown in Table 4, the predicted values ​​of elemental sulfur solubility according to the present invention are very close to the measured values, with an average relative error of only 9.64%.

[0188] Example 6

[0189] The elemental sulfur solubility of the target gas was predicted using the elemental sulfur solubility prediction method according to the present invention. In this embodiment, the composition (mole fraction) of the target gas was 9.93% H2S, 7.16% CO2, and 82.91% CH4, and the temperature and pressure of the target gas are shown in Table 5.

[0190] According to the elemental sulfur solubility prediction method of the present invention, the target gas density is calculated according to steps S110 to S150, and the calculated target gas density is shown in Table 5.

[0191] In this embodiment, the target gas temperature is 303.15K ≤ T < 373.15K, therefore the prediction model is selected as follows. d = 0.0079264·T - 2.2623,

[0192] H2S mole fraction y H2S The percentage is 9.93%. According to step S222, the model parameters are:

[0193] k=-0.0016125·T+1.9862, a=-5973.8, b=7.0301.

[0194] Based on the selected prediction model and the determined model parameters, the solubility of elemental sulfur in the target gas is calculated, the predicted value of elemental sulfur solubility is obtained, and the relative error is calculated by comparing it with the measured value, as shown in Table 5.

[0195] Table 5. Comparison of experimental conditions and measured and predicted values ​​of elemental sulfur solubility.

[0196]

[0197] As shown in Table 5, the predicted values ​​of elemental sulfur solubility according to the present invention are very close to the measured values, with an average relative error of only 8.12%.

[0198] Example 7

[0199] The elemental sulfur solubility prediction method according to the present invention was used to predict the elemental sulfur solubility of the target gas. In this embodiment, the composition (molar fraction) of the target gas is 17.71% H2S, 6.81% CO2 and 75.48% CH4, and the temperature and pressure of the target gas are shown in Table 6.

[0200] According to the elemental sulfur solubility prediction method of the present invention, the target gas density is calculated according to steps S110 to S150, and the calculated target gas density is shown in Table 6.

[0201] In this embodiment, the target gas temperature is 303.15K ≤ T < 373.15K, therefore the prediction model is selected as follows. d = 0.0079264·T - 2.2623,

[0202] H2S mole fraction The percentage is 17.71%. According to step S221, the model parameters are:

[0203] k=-0.00558·T+3.8555, a=-9425.1, b=14.654.

[0204] Based on the selected prediction model and the determined model parameters, the solubility of elemental sulfur in the target gas is calculated, the predicted value of elemental sulfur solubility is obtained, and the relative error is calculated by comparing it with the measured value, as shown in Table 6.

[0205] Table 6. Comparison of experimental conditions and measured and predicted values ​​of elemental sulfur solubility.

[0206]

[0207] As shown in Table 6, the predicted values ​​of elemental sulfur solubility according to the present invention are very close to the measured values, with an average relative error of only 10.21%.

[0208] Example 8

[0209] The elemental sulfur solubility prediction method according to the present invention was used to predict the elemental sulfur solubility of the target gas. In this embodiment, the composition (mole fraction) of the target gas is 10% H2S, 0.86% CO2 and 89.14% CH4, and the temperature and pressure of the target gas are shown in Table 7.

[0210] According to the elemental sulfur solubility prediction method of the present invention, the target gas density is calculated according to steps S110 to S150, and the calculated target gas density is shown in Table 7.

[0211] In this embodiment, the target gas temperature is 303.15K ≤ T < 373.15K, therefore the prediction model is selected as follows. d = 0.0079264·T - 2.2623,

[0212] H2S mole fraction The value is 10%. According to step S222, the model parameters are:

[0213] k=-0.0016125·T+1.9862, a=-5973.8, b=7.0301.

[0214] Based on the selected prediction model and the determined model parameters, the solubility of elemental sulfur in the target gas is calculated, the predicted value of elemental sulfur solubility is obtained, and the relative error is calculated by comparing it with the measured value, as shown in Table 7.

[0215] Table 7 Comparison of experimental conditions and measured and predicted values ​​of elemental sulfur solubility.

[0216]

[0217] As shown in Table 7, the predicted values ​​of elemental sulfur solubility according to the present invention are very close to the measured values, with an average relative error of only 5.97%.

[0218] Example 9

[0219] The elemental sulfur solubility prediction method according to the present invention was used to predict the elemental sulfur solubility of the target gas. In this embodiment, the composition (mole fraction) of the target gas is 10.03% H2S, 10.39% CO2 and 79.58% CH4, and the temperature and pressure of the target gas are shown in Table 7.

[0220] According to the elemental sulfur solubility prediction method of the present invention, the target gas density is calculated according to steps S110 to S150, and the calculated target gas density is shown in Table 8.

[0221] In this embodiment, the target gas temperature is 303.15K ≤ T < 373.15K, therefore the prediction model is selected as follows. d = 0.0079264·T - 2.2623,

[0222] H2S mole fraction The value is 10.03%. According to step S222, the model parameters are:

[0223] k=-0.0016125·T+1.9862, a=-5973.8, b=7.0301.

[0224] Based on the selected prediction model and the determined model parameters, the solubility of elemental sulfur in the target gas is calculated, the predicted value of elemental sulfur solubility is obtained, and the relative error is calculated by comparing it with the measured value, as shown in Table 8.

[0225] Table 8. Comparison of experimental conditions and measured and predicted values ​​of elemental sulfur solubility.

[0226]

[0227] As shown in Table 8, the predicted values ​​of elemental sulfur solubility according to the present invention are very close to the measured values, with an average relative error of only 10.54%.

[0228] Example 10

[0229] This embodiment provides a non-transitory (non-volatile) computer storage medium that stores computer-executable instructions that can execute the methods in any of the above method embodiments and achieve the same technical effect.

[0230] Example 11

[0231] This embodiment provides a computer program product, which includes a computer program stored on a non-transitory computer-readable storage medium. The computer program includes program instructions, which, when executed by a computer, cause the computer to perform the methods described above and achieve the same technical effects.

[0232] Example 12

[0233] Figure 3 This is a schematic diagram of the hardware structure of the electronic device for performing the elemental sulfur solubility prediction method according to this embodiment. The device includes one or more processors 610 and a memory 620. Taking one processor 610 as an example, the device may also include an input device 630 and an output device 640.

[0234] The processor 610, memory 620, input device 630, and output device 640 can be connected via a bus or other means. Figure 3 Taking the example of a connection between China and Israel via a bus.

[0235] The memory 620, as a non-transitory computer-readable storage medium, can be used to store non-transitory software programs, non-transitory computer-executable programs, and modules. The processor 610 executes various functional applications and data processing of the electronic device by running the non-transitory software programs, instructions, and modules stored in the memory 620, thereby implementing the processing method of the above-described method embodiments.

[0236] The memory 620 may include a program storage area and a data storage area, wherein the program storage area may store the operating system and applications required for at least one function; the data storage area may store data, etc. Furthermore, the memory 620 may include high-speed random access memory, and may also include non-transitory memory, such as at least one disk storage device, flash memory device, or other non-transitory solid-state storage device. In some embodiments, the memory 620 may optionally include memory remotely located relative to the processor 610, and these remote memories may be connected to the processing device via a network. Examples of such networks include, but are not limited to, the Internet, corporate intranets, local area networks, mobile communication networks, and combinations thereof.

[0237] Input device 630 can receive input digital or character information and generate signal input. Output device 640 may include display devices such as a display screen.

[0238] One or more modules are stored in memory 620 and, when executed by one or more processors 610, execute:

[0239] Obtain the parameter information of the target gas and calculate the density ρ of the target gas;

[0240] Based on the target gas temperature T, a prediction model is selected, and based on the H2S mole fraction... Given the target gas density ρ, determine the model parameters:

[0241] When 373.15K≤T≤433.15K, the prediction model is c=ρ k exp(a / T+b);

[0242] When 303.15K ≤ T < 373.15K, the prediction model is:

[0243] In the formula, c is the solubility of elemental sulfur in the target gas, k, a and b are model parameters, model parameter k is a function of T, and the low temperature correction coefficient d = 0.0079264·T-2.2623;

[0244] Based on the selected prediction model and the determined model parameters, the solubility of elemental sulfur in the target gas is calculated.

[0245] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate, and 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 modules can be selected to achieve the purpose of this embodiment according to actual needs.

[0246] Through the above description of the embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus a general-purpose hardware platform, and of course, it can also be implemented by hardware. Based on this understanding, the above technical solutions, in essence or the parts that contribute to the related technology, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., including several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute the methods of various embodiments or some parts of embodiments.

[0247] The foregoing description of specific exemplary embodiments of the present invention is for illustrative and explanatory purposes. These descriptions are not intended to limit the invention to the precise forms disclosed, and it will be apparent that many changes and variations can be made in accordance with the foregoing teachings. The exemplary embodiments were chosen and described in order to explain the specific principles of the invention and its practical application, thereby enabling those skilled in the art to implement and utilize various different exemplary embodiments of the invention, as well as various different choices and variations. Any simple modifications, equivalent changes, and alterations made to the foregoing exemplary embodiments should fall within the scope of protection of the present invention.

Claims

1. A method for predicting the solubility of elemental sulfur, characterized in that, Includes the following steps: Obtain the parameter information of the target gas and calculate the density of the target gas. ; Based on target gas temperature Select a prediction model and base it on H2S mole fraction. and target gas density Determine the model parameters: When 373.15K≤ When K ≤ 433.15, the prediction model is: When the prediction model is At that time, based on the H2S mole fraction and target gas density The model parameters were determined based on the H2S mole fraction. Calculate the density inflection point ; and compare the density of the target gas and density inflection point The model parameters are determined based on the comparison results; When 303.15K≤ When <373.15K, the prediction model is , In the formula The solubility of sulfur in the target gas. , and For model parameters, model parameters for The function, low temperature correction coefficient ;as well as Based on the selected prediction model and the determined model parameters, the solubility of elemental sulfur in the target gas is calculated.

2. The method for predicting the solubility of elemental sulfur according to claim 1, characterized in that, When the prediction model is When determining model parameters, the following are included: When 13.5%≤ When ≤25%, the model parameters are: ； When 6.5%≤ When <13.5%, the model parameters are: ； When 3.5%≤ When <6.5%, the model parameters are: ;as well as When 1%≤ When <3.5%, the model parameters are: 。 3. The method for predicting the solubility of elemental sulfur according to claim 1, characterized in that, When 13.5%≤ Calculate the density inflection point when ≤25%. , like < Then the model parameters are: ， like ≥ Then the model parameters are: ； When 6.5%≤ When <13.5%, density inflection point , like < Then the model parameters are: ， like ≥ Then the model parameters are: ； When 3.5%≤ When <6.5%, density inflection point , like < Then the model parameters are: ， like ≥ Then the model parameters are: ； as well as When 1%≤ When <3.5%, density inflection point , like < Then the model parameters are: ， like ≥ Then the model parameters are: 。 4. The method for predicting the solubility of elemental sulfur according to claim 1, characterized in that, The parameter information of the target gas includes the target gas temperature. Target gas pressure The acidic gas components in the target gas and the mole fraction of each component.

5. The method for predicting the solubility of elemental sulfur according to claim 1, characterized in that, Calculate the target gas density include: Calculate the apparent critical pressure based on the parameter information of the target gas. Apparent critical temperature and apparent relative molecular mass ; Based on H2S mole fraction and CO2 mole fraction , for visual critical pressure and apparent critical temperature Perform correction; Based on the target gas parameters, the corrected apparent critical pressure, and the apparent critical temperature, calculate the apparent corresponding pressure. and corresponding temperature ; According to the corresponding pressure and corresponding temperature Consult the natural gas compressibility factor chart to obtain the corresponding compressibility factor value for the target gas. ;as well as Using target gas pressure Target gas temperature Apparent relative molecular mass Corresponding compression factor value Calculate the density of the target gas ， In the formula, R = 8.314 J / (mol·K) is the universal gas constant.

6. The method for predicting the solubility of elemental sulfur according to claim 5, characterized in that, Apparent critical pressure Apparent critical temperature and apparent relative molecular mass The calculation formulas are as follows: ， In the formula, p ci Components i The critical pressure, in MPa; T ci Components i The critical temperature, expressed in K; M i Components i The relative molecular mass; y i Components i The mole fraction.

7. The method for predicting the solubility of elemental sulfur according to claim 6, characterized in that, Corrected apparent critical pressure and the corrected apparent critical temperature They are respectively: ， ， ， In the formula, A is the correction factor; A is the sum of the mole fractions of H2S and CO2; B is the mole fraction of H2S.

8. The method for predicting the solubility of elemental sulfur according to claim 7, characterized in that, Depending on the corresponding pressure and corresponding temperature The calculation formulas are as follows: 。 9. The method for predicting the solubility of elemental sulfur according to claim 1, characterized in that, The target gas has a pressure range of 10~60MPa, a temperature range of 303.15K~433.15K, an H2S mole fraction range of 1%~25%, and a CO2 mole fraction range of 0~25%.

10. A system for predicting the solubility of elemental sulfur, characterized in that, include: The data acquisition unit is used to acquire parameter information of the target gas; The computing unit is used to calculate the density of the target gas. ; Model acquisition unit, which is used to acquire models based on target gas temperature Select a prediction model and base it on H2S mole fraction. Determine the model parameters: When 373.15K≤ When K ≤ 433.15, the prediction model is: When the prediction model is At that time, based on the H2S mole fraction and target gas density The model parameters were determined based on the H2S mole fraction. Calculate the density inflection point ; And compare the density of the target gas and density inflection point The model parameters are determined based on the comparison results; When 303.15K≤ When <373.15K, the prediction model is , In the formula The solubility of sulfur in the target gas. , and For model parameters, model parameters for The function, low temperature correction coefficient ;as well as The prediction unit calculates the solubility of elemental sulfur in the target gas based on the selected prediction model and the determined model parameters.

11. An electronic device, characterized in that, include: At least one processor; as well as A memory communicatively connected to the at least one processor; wherein, The memory stores instructions that can be executed by the at least one processor to cause the at least one processor to perform the elemental sulfur solubility prediction method as described in any one of claims 1 to 9.

12. A non-transitory computer-readable storage medium, characterized in that, The non-transitory computer-readable storage medium stores computer-executable instructions for causing the computer to perform the elemental sulfur solubility prediction method as described in any one of claims 1 to 9.

Images