CN113237917B - Method and device for calculating methane adsorption quantity in coal at different temperatures - Google Patents

Method and device for calculating methane adsorption quantity in coal at different temperatures Download PDF

Info

Publication number
CN113237917B
CN113237917B CN202110426819.4A CN202110426819A CN113237917B CN 113237917 B CN113237917 B CN 113237917B CN 202110426819 A CN202110426819 A CN 202110426819A CN 113237917 B CN113237917 B CN 113237917B
Authority
CN
China
Prior art keywords
adsorption
temperature
under
isothermal
isothermal adsorption
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110426819.4A
Other languages
Chinese (zh)
Other versions
CN113237917A (en
Inventor
邹卓
宁树正
张莉
刘亢
黄少青
孙杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
General Survey and Research Institute of China Coal Geology Bureau
Original Assignee
General Survey and Research Institute of China Coal Geology Bureau
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by General Survey and Research Institute of China Coal Geology Bureau filed Critical General Survey and Research Institute of China Coal Geology Bureau
Priority to CN202110426819.4A priority Critical patent/CN113237917B/en
Publication of CN113237917A publication Critical patent/CN113237917A/en
Application granted granted Critical
Publication of CN113237917B publication Critical patent/CN113237917B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N25/00Investigating or analyzing materials by the use of thermal means
    • G01N25/20Investigating or analyzing materials by the use of thermal means by investigating the development of heat, i.e. calorimetry, e.g. by measuring specific heat, by measuring thermal conductivity
    • G01N25/48Investigating or analyzing materials by the use of thermal means by investigating the development of heat, i.e. calorimetry, e.g. by measuring specific heat, by measuring thermal conductivity on solution, sorption, or a chemical reaction not involving combustion or catalytic oxidation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • G01N2015/0866Sorption

Landscapes

  • Chemical & Material Sciences (AREA)
  • Biochemistry (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Dispersion Chemistry (AREA)
  • Engineering & Computer Science (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Combustion & Propulsion (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

The invention provides a method and a device for calculating methane adsorption capacity in coal at different temperatures based on an adsorption heat theory, wherein the method comprises the following steps: acquiring isothermal adsorption data of methane in coal under a group of first temperature conditions; respectively inputting isothermal adsorption quantity value and adsorption pressure in isothermal adsorption data into an equal adsorption heat and isothermal adsorption quantity relation model and an adsorption pressure and fugacity relation model to obtain a group of equal adsorption heat value and first fugacity, and inputting into two fugacity and equal adsorption heat relation models at different temperatures to obtain second fugacity under a second temperature condition; taking the isothermal adsorption quantity corresponding to the first loss as the isothermal adsorption quantity of the second loss at the second temperature corresponding to the first loss; and performing curve fitting by taking the second loss and the isothermal adsorption quantity as an abscissa and an ordinate respectively to obtain an isothermal adsorption curve, and further determining isothermal adsorption data. The invention can effectively understand the adsorption characteristics of coal to methane at different temperatures and assist in judging and evaluating the coal seam.

Description

Method and device for calculating methane adsorption quantity in coal at different temperatures
Technical Field
The invention relates to the technical field of mineral exploitation, in particular to a method and a device for calculating methane adsorption quantity in coal at different temperatures.
Background
The coal is a porous medium with a pore-crack dual structure, gas in the coal mainly exists in an adsorption mode, and the adsorbed gas content in the coal can be more than 80-90% of the total adsorption quantity. The adsorption characteristic of the coal on methane gas can be known by calculating the methane content in the coal, and the gas storage capacity of the palm holding coal plays an important role in judging phenomena such as coal bed gas enrichment, gas outburst and the like, evaluating the gas content of the coal bed and predicting the capacity of the coal bed gas well.
At present, the method for calculating the methane content in the coal is to perform isothermal adsorption experiments, and because the adsorption performance of the coal on methane at different temperatures is different, a plurality of groups of isothermal adsorption experiments for methane at different temperatures are required to be performed for calculating the methane content adsorbed by the coal at different temperatures.
The isothermal adsorption experiment at different temperatures is very long in time consumption, large in workload and easy to generate experimental errors, and various cost charges are high.
Disclosure of Invention
The invention provides a method and a device for calculating methane adsorption capacity in coal at different temperatures, which are used for solving the defects of long time consumption, large workload, easy experimental error and high cost of isothermal adsorption experiments performed at different temperatures in the prior art, saving time, reducing workload, improving accuracy of determining methane adsorption capacity in coal, reducing experimental error and reducing cost.
In a first aspect, the invention provides a method for calculating methane adsorption capacity in coal at different temperatures based on a heat of adsorption theory, comprising the following steps:
acquiring isothermal adsorption data of methane in coal under a group of first temperature conditions;
inputting a set of adsorption pressure in isothermal adsorption data of methane in coal under the first temperature condition into an adsorption pressure and fugacity relation model to obtain a first fugacity under the first temperature condition;
inputting a set of isothermal adsorption values in isothermal adsorption data of methane in the coal under the first temperature condition into an equivalent adsorption heat and isothermal adsorption quantity relation model to obtain a set of equivalent adsorption heat values corresponding to the set of isothermal adsorption values;
inputting the first loss degree and the set of equivalent adsorption heat values under the first temperature condition into a relation model between the loss degrees and equivalent adsorption heat under two different temperatures to obtain a set of second loss degrees under the second temperature condition corresponding to the first loss degree under the first temperature condition;
taking the isothermal adsorption quantity corresponding to any first loss degree under the first temperature condition as the isothermal adsorption quantity under the second temperature condition corresponding to the second loss degree under the second temperature condition corresponding to the first loss degree under the first temperature condition;
and performing curve fitting by taking the second loss degree at the second temperature as an abscissa and the isothermal adsorption quantity at the second temperature as an ordinate to obtain an isothermal adsorption curve at the second temperature, and further determining isothermal adsorption data at the second temperature according to the isothermal adsorption curve.
Further, according to the method for calculating the methane adsorption quantity in the coal at different temperatures based on the adsorption heat theory, the relation model of the adsorption pressure and the fugacity is as follows:
Ln(f/P)=Z-1-Ln(Z-B)-(A/B)Ln(1+B/Z)
wherein Z is a compression factor, dimensionless; p is the pressure, MPa; A. b is a coefficient related to the pressure P, and is specifically expanded as follows: a=c 1 P/(R 2 T 2 ),B=c 2 P/(R 2 T 2 ),
Figure BDA0003029895770000021
c 2 =0.0866RT c /P c , m=0.48+1.534ω-0.176ω 2 Wherein Pc is the critical pressure of gas, and the unit is MPa; tc is the critical temperature of the gas, and the unit is K; TR is the comparative temperature, tr=t/Tc, T representing temperature; c1 and c2 are constants; omega is an eccentric factor and is dimensionless.
Further, according to the method for calculating the methane adsorption capacity in the coal at different temperatures based on the adsorption heat theory provided by the invention, the relation model of the equivalent adsorption heat and isothermal adsorption capacity is as follows:
q st =C 5 n+C 6
wherein q st Is equal in adsorption heat, n is isothermal adsorption quantity, C5 and C6 are q st Fitting parameters to n.
Further, according to the method for calculating the methane adsorption capacity in the coal at different temperatures based on the adsorption heat theory provided by the invention, a relation model between the fugacity at the two different temperatures and the equivalent adsorption heat is as follows:
Figure RE-GDA0003118729320000031
wherein q st Represents equivalent adsorption heat, and the unit is kJ/mol; t (T) 1 And T 2 Temperature is represented in K; r is a gas constant, f 1 And f 2 The loss was measured.
In a second aspect, the invention provides a device for calculating methane adsorption capacity in coal at different temperatures based on a heat of adsorption theory, comprising:
the first processing module is used for acquiring isothermal adsorption data of methane in the coal under a group of first temperature conditions;
the second processing module is used for inputting a set of adsorption pressure in isothermal adsorption data of methane in the coal under the first temperature condition into an adsorption pressure and fugacity relation model to obtain a first fugacity under the first temperature condition;
the third processing module is used for inputting a group of isothermal adsorption quantity values in isothermal adsorption data of methane in the coal under the group of first temperature conditions into an equivalent adsorption heat and isothermal adsorption quantity relation model to obtain a group of equivalent adsorption heat values corresponding to the group of isothermal adsorption quantity values;
the fourth processing module is used for inputting the first loss degree and the set of equivalent adsorption heat values under the first temperature condition into a relation model between the loss degrees and the equivalent adsorption heat under two different temperatures to obtain a set of second loss degrees under the second temperature condition corresponding to the first loss degrees under the first temperature condition;
a fifth processing module, configured to use the isothermal adsorption amount corresponding to any one of the first fugacity under the first temperature condition as the isothermal adsorption amount under the second temperature condition corresponding to the second fugacity under the second temperature condition corresponding to the first fugacity under the first temperature condition;
and the sixth processing module is used for performing curve fitting by taking the second loss degree at the second temperature as an abscissa and the isothermal adsorption quantity at the second temperature as an ordinate to obtain an isothermal adsorption curve at the second temperature, and further determining isothermal adsorption data at the second temperature according to the isothermal adsorption curve.
Further, according to the device for calculating the adsorption amount of methane in coal at different temperatures based on the adsorption heat theory provided by the invention, the relation model of adsorption pressure and fugacity in the second processing module is as follows:
Ln(f/P)=Z-1-Ln(Z-B)-(A/B)Ln(1+B/Z)
wherein Z is a compression factor, dimensionless; p is the pressure, MPa; A. b is a coefficient related to the pressure P, and is specifically expanded as follows: a=c 1 P/(R 2 T 2 ),B=c 2 P/(R 2 T 2 ),
Figure BDA0003029895770000041
c 2 =0.0866RT c /P c , m=0.48+1.534ω-0.176ω 2 Wherein Pc is the critical pressure of gas, and the unit is MPa; tc is the critical temperature of the gas, and the unit is K; TR is the comparative temperature, tr=t/Tc, T representing temperature; c1 and c2 are constants; omega is an eccentric factor and is dimensionless.
Further, according to the device for calculating the adsorption capacity of methane in coal at different temperatures based on the adsorption heat theory provided by the invention, the relation model of the equivalent adsorption heat and isothermal adsorption capacity in the third processing module is as follows:
q st =C 5 n+C 6
wherein q st Is equal in adsorption heat, n is isothermal adsorption quantity, C5 and C6 are q st Fitting parameters to n.
Further, according to the device for calculating the adsorption amount of methane in coal at different temperatures based on the adsorption heat theory provided by the invention, a relation model between the loss degree and the equivalent adsorption heat at the two different temperatures in the fifth module is as follows:
Figure RE-GDA0003118729320000051
wherein q st Represents equivalent adsorption heat, and the unit is kJ/mol; t (T) 1 And T 2 The temperature is indicated as a function of the temperature,the unit is K; r is a gas constant, f 1 And f 2 The loss was measured.
In a third aspect, the present invention provides an electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor, when executing the program, implements the steps of any of the methods for calculating methane adsorption in coal at different temperatures based on heat of adsorption theory described above.
In a fourth aspect, the present invention provides a non-transitory computer readable storage medium having stored thereon a computer program, characterized in that the computer program, when executed by a processor, implements the steps of the method for calculating methane adsorption in coal at different temperatures based on heat of adsorption theory as described in any of the above.
According to the method and the device for calculating the methane adsorption capacity in the coal under different temperatures, isothermal adsorption data of methane in the coal under a group of first temperature conditions are obtained; the first escape degree of the group of first temperature conditions is obtained by inputting a group of adsorption pressures in isothermal adsorption data of methane in coal under the group of first temperature conditions into an adsorption pressure and escape degree relation model; in addition, a group of isothermal adsorption values in isothermal adsorption data of methane in the coal under the group of first temperature conditions are input into an equivalent adsorption heat and isothermal adsorption quantity relation model, and a group of equivalent adsorption heat values corresponding to the group of isothermal adsorption values are obtained; then inputting the obtained first loss degree and the obtained set of equivalent adsorption heat values under the first temperature condition into a relation model between the loss degrees and the equivalent adsorption heat at two different temperatures to obtain a set of second loss degrees under the second temperature condition corresponding to the first loss degree under the first temperature condition; in addition, the isothermal adsorption amount corresponding to any one of the first fugacity under the first temperature condition is taken as the isothermal adsorption amount corresponding to the second fugacity under the second temperature condition corresponding to the first fugacity under the first temperature condition; and then, performing curve fitting by taking the second loss at the second temperature as an abscissa and the isothermal adsorption at the second temperature as an ordinate to obtain an isothermal adsorption curve at the second temperature, and further determining isothermal adsorption data at the second temperature according to the isothermal adsorption curve. According to the method, the adsorption characteristics of the coal to methane at different temperatures are effectively known, the phenomena of coal bed gas enrichment, gas outburst and the like are facilitated to be judged, the gas content of the coal bed is evaluated, the capacity of the coal bed gas well is predicted, meanwhile, the experimental time can be greatly saved, the experimental workload is greatly reduced, and the experimental cost is reduced.
Drawings
In order to more clearly illustrate the invention or the technical solutions of the prior art, the following description will briefly explain the embodiments or the drawings needed in the description of the prior art, and it is obvious that the drawings in the following description are some embodiments of the invention and that other drawings can be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a schematic flow chart of a method for calculating methane adsorption capacity in coal at different temperatures based on the adsorption heat theory;
FIG. 2 is a graph of Lnf-n fit of coal provided in an example of the invention at different temperatures;
FIG. 3 is a graph of Lnf-1/T fit for coal provided in the examples of the invention at different fixed adsorption capacities;
FIG. 4 is a graph showing isothermal adsorption capacity versus corresponding equivalent adsorption heat for coal provided in the examples of the present invention;
FIG. 5 is a schematic structural diagram of the device for calculating the methane adsorption capacity in coal at different temperatures based on the adsorption heat theory;
fig. 6 is a schematic structural diagram of an electronic device provided by the present invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the present invention more apparent, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
In the following, referring to fig. 1, a method for calculating methane adsorption amount in coal at different temperatures based on the theory of heat of adsorption is described in the embodiment of the present invention, which includes:
step 100: acquiring isothermal adsorption data of methane in coal under a group of first temperature conditions;
step 200: inputting a set of adsorption pressure in isothermal adsorption data of methane in coal under the first temperature condition into an adsorption pressure and fugacity relation model to obtain a first fugacity under the first temperature condition;
step 300: inputting a set of isothermal adsorption values in isothermal adsorption data of methane in the coal under the first temperature condition into an equivalent adsorption heat and isothermal adsorption quantity relation model to obtain a set of equivalent adsorption heat values corresponding to the set of isothermal adsorption values;
step 400: inputting the first loss degree and the set of equivalent adsorption heat values under the first temperature condition into a relation model between the loss degrees and the equivalent adsorption heat under two different temperatures to obtain a set of second loss degrees under the second temperature condition corresponding to the first loss degrees under the first temperature condition;
step 500: taking the isothermal adsorption quantity corresponding to any first loss degree under the first temperature condition as the isothermal adsorption quantity under the second temperature condition corresponding to the second loss degree under the second temperature condition corresponding to the first loss degree under the first temperature condition;
step 600: and performing curve fitting by taking the second loss degree at the second temperature as an abscissa and the isothermal adsorption quantity at the second temperature as an ordinate to obtain an isothermal adsorption curve at the second temperature, and determining isothermal adsorption data at the second temperature according to the isothermal adsorption curve.
Specifically, in the theoretical model of heat of adsorption, the equivalent heat of adsorption, also called differential heat of adsorption, refers to the amount of heat released when an infinitely small amount of gas molecules is adsorbed at a constant adsorption amount, and is the instantaneous enthalpy change in the adsorption process. The equivalent heat of adsorption in the adsorption process can be calculated using the Claus ius-Clapeyron equation, i.e.
Figure RE-GDA0003118729320000081
Wherein q is st kJ/mol is equal amount of adsorption heat; t represents temperature, and the unit is K; p is pressure in MPa; r is the gas constant, 8.314J/(mol.K).
In the embodiment of the present invention, in step 100, obtaining a set of isothermal adsorption data of methane in coal under a first temperature condition refers to obtaining, by an experimental method, a set of adsorption amounts of methane in coal under the same temperature condition and different pressures, or knowing a set of isothermal adsorption data. The isothermal adsorption data of methane in coal under the first temperature condition at least comprise pressure data and adsorption amount data. Meanwhile, the adsorption amount data of the pressure data spear are corresponding, namely one pressure corresponds to one adsorption amount data. And a set of data does not refer to one pair of pressure and adsorption amount data, but to a plurality of pairs of pressure and coefficient amount data.
In step 200, the pressure in the isothermal adsorption data at the first temperature acquired in step 100 is converted to a fugacity, because methane is a non-ideal gas in the supercritical temperature region, and therefore, in order to reduce the error, the pressure P is replaced with a fugacity f, wherein the fugacity in chemico-thermo-mechanical represents the effective pressure of the actual gas, which is equal to the pressure of the ideal gas having the same chemical potential under the same conditions. The activity is replaced by the fugacity in the thermodynamic equation describing the ideal gas property related to chemical potential, and the corresponding relation describing the actual gas property can be obtained. f dimensionless:
Figure RE-GDA0003118729320000082
in step 300, the relationship model between the equivalent adsorption heat and the isothermal adsorption amount is constructed in advance, and then the isothermal adsorption amount in isothermal adsorption data of methane in coal at the first temperature obtained in step 100 is combined, so that the equivalent adsorption heat corresponding to each isothermal adsorption amount can be obtained.
In step 400, through a pre-constructed relation model between the two loss values at different temperatures and the equivalent adsorption heat, the equivalent adsorption heat value corresponding to the isothermal adsorption amount at the first temperature obtained in step 300 and the first loss value at the first temperature obtained in step 200 is brought into the pre-constructed model, so as to obtain a set of second loss values at the second temperature corresponding to the first loss values at the first temperature.
In step 500, the isothermal adsorption amount corresponding to any one of the first fugacity under the first temperature condition is set as the isothermal adsorption amount under the second temperature condition corresponding to the second fugacity under the second temperature condition corresponding to the first fugacity under the first temperature condition; the isothermal adsorption amount at the first loss of the first temperature is equal to the isothermal adsorption amount corresponding to the second loss of the second temperature corresponding to the first loss of the first temperature, and the isothermal adsorption amounts are the same or close to each other when the isothermal adsorption amounts are the same or have small differences, so that the gas adsorption amounts corresponding to the first loss and the second loss of the second temperature under different temperature conditions can be considered to be equal. By this conversion, isothermal adsorption values corresponding to different fugacity values at the second temperature condition can be obtained.
Further, in step 600, the second fugacity obtained in step 400 is plotted on the abscissa, and the adsorption amount corresponding to the second fugacity under the second temperature condition obtained in step 500 is plotted on the ordinate, i.e. the plotted fit curve is the isothermal adsorption curve under the second temperature condition to be predicted. The adsorption curve is a relation curve representing the fugacity and the adsorption amount, and the relationship curve can be converted between the fugacity and the pressure. This curve can be seen as a curve between pressure and the amount of adsorption of the coal under the second temperature condition. The relationship between the pressure and the adsorption amount under the second temperature condition can be predicted by the curve, and the adsorption data under the second temperature condition can be obtained by the obtained fitting curve because the adsorption data comprise the pressure index and the adsorption amount index.
In the embodiment, the invention obtains a group of isothermal adsorption data of methane in coal under a first temperature condition; inputting a set of adsorption pressure in isothermal adsorption data of methane in coal under the first temperature condition into an adsorption pressure and fugacity relation model to obtain a first fugacity under the first temperature condition; in addition, a group of isothermal adsorption values in isothermal adsorption data of methane in the coal under the group of first temperature conditions are input into an equal adsorption heat and equal adsorption quantity relation model, and a group of equal adsorption heat values corresponding to the group of isothermal adsorption values are obtained; then inputting the obtained first loss degree and the obtained set of equivalent adsorption heat values under the first temperature condition into a relation model between the loss degrees and the equivalent adsorption heat at two different temperatures to obtain a set of second loss degrees under the second temperature condition corresponding to the first loss degrees under the first temperature condition; in addition, the isothermal adsorption amount corresponding to any one of the first fugacity under the first temperature condition is taken as the isothermal adsorption amount corresponding to the second fugacity under the second temperature condition corresponding to the first fugacity under the first temperature condition; and then, performing curve fitting by taking the second loss at the second temperature as an abscissa and the isothermal adsorption at the second temperature as an ordinate to obtain an isothermal adsorption curve at the second temperature, and further determining isothermal adsorption data at the second temperature according to the isothermal adsorption curve. According to the method, the adsorption characteristics of the coal to methane at different temperatures are effectively known, the phenomena of coal bed gas enrichment, gas outburst and the like are facilitated to be judged, the gas content of the coal bed is evaluated, the capacity of the coal bed gas well is predicted, meanwhile, the experimental time can be greatly saved, the experimental workload is greatly reduced, and the experimental cost is reduced.
Further, the present invention provides a method for calculating the amount of methane adsorbed in coal at different temperatures based on the heat of adsorption theory in the examples, wherein the substitution of the fugacity for the pressure is indicated in the foregoing, and specifically the conversion equation between the fugacity and the pressure is:
Ln(f/P)=Z-1-Ln(Z-B)-(A/B)Ln(1+B/Z) (3)
wherein Z is a compression factor, dimensionless; p is the pressure, MPa; A. b is a coefficient related to the pressure P, and is specifically expanded as follows:
A=c 1 P/(R 2 T 2 ),B=c 2 P/(R 2 T 2 ),
c 1 =0.4275R 2 T c 2 [1+m(1-T R 0.5 )] 2 /P c ,c 2 =0.0866RT c /P c , m=0.48+1.534ω-0.176ω 2
wherein P is c Is the critical pressure of gas, MPa; t (T) c Is the critical temperature of the gas, K; t (T) R For comparison of temperature, T R =T/T c ;c 1 、c 2 Is a constant; omega is an eccentric factor, dimensionless, omega takes 0.008.
By the above formula, the pressure can be converted into the corresponding fugacity. The conversion model of escape degree and pressure thus constructed can be expressed as:
Ln(f/P)=Z-1-Ln(Z-B)-(A/B)Ln(1+B/Z)
further, according to the method for calculating the methane adsorption capacity in the coal at different temperatures based on the adsorption heat theory provided by the invention, the relation model of the equivalent adsorption heat and isothermal adsorption capacity is as follows:
q st =C 5 n+C 6
wherein q st Is equal in adsorption heat, n is isothermal adsorption quantity, C5 and C6 are q st Fitting parameters to n.
Specifically, the method for constructing the relation model between the equivalent adsorption heat and the isothermal adsorption quantity comprises the following steps:
given isothermal adsorption data (i.e., adsorption pressure p—adsorption amount n) at any three temperatures for a set of samples, pressure P can be converted to fugacity f according to equation (3) above, and a linear fit can be made to Lnf and adsorption amount n as shown in fig. 2 ( temperatures 298K, 303K, 308K are examples). Lnf and isothermal adsorption n at different temperatures form a better linear relationship, and the fitting relation is as follows:
Lnf=C 1 n+C 2 (4)
wherein C is 1 、C 2 Is a fitting parameter for Lnf and n.
At 298K, 303K and 308K, 8 adsorption values are respectively set, 0.1-0.8mmol/g is taken and brought into formula (4), lnf at different temperatures is obtained, and Lnf and 1/T are linearly fitted, as shown in FIG. 3. Lnf has a better linear relationship with 1/T at different isothermal adsorption amounts. The fitting formula is as follows:
Lnf=C 3 /T+C 4 (5)
wherein C is 3 、C 4 Is a fitting parameter.
Integrating the two ends of equation (2) above yields a relationship between Lnf and 1/T:
Lnf=-q st /(RT)+C (6)
from the above formula (5) and formula (6), it can be obtained:
q st =-RC 3 (7)
therefore, the equivalent heat of adsorption q st Can be obtained by the slope of the Lnf-1/T relationship.
And drawing a relation curve (figure 4) of the equivalent adsorption heat and the isothermal adsorption amount of different coal rank coals according to the obtained equivalent adsorption heat and isothermal adsorption amount, wherein the fitting relation is as follows:
q st =C 5 n+C 6 (8)
wherein C is 5 、C 6 Is q st Fitting parameters to n.
Through the steps, the relation between the equivalent adsorption heat and the isothermal adsorption quantity can be obtained, namely, the relation model of the equivalent adsorption heat and the isothermal adsorption quantity is determined as follows:
q st =C 5 n+C 6
further, the method for calculating the adsorption capacity of methane in coal at different temperatures based on the adsorption heat theory provided by the embodiment of the invention, wherein a relation model between the fugacity at the two different temperatures and the equivalent adsorption heat is as follows:
Figure RE-GDA0003118729320000121
wherein q st Represents equivalent adsorption heat, and the unit is kJ/mol; t (T) 1 And T 2 Temperature is represented in K; r is a gas constant, f 1 And f 2 The loss was measured.
Specifically, when the known coal sample is T at any temperature 1 Isothermal adsorption data at the time of the process and the equivalent heat of adsorption obtained by the method can be calculated to obtain any temperature T 2 Isothermal adsorption curve under. Transforming the formula (6) to obtain a temperature T 1 、T 2 The relationship between the following fugacity and the equivalent heat of adsorption is as follows:
Figure BDA0003029895770000122
the method comprises the following steps of:
Figure RE-GDA0003118729320000123
wherein q is st Represents equivalent heat of adsorption, kJ/mol; t represents temperature, K; r is the gas constant, 8.314J/(mol.K), and f is the fugacity. Thus when knowing the temperature T 1 The fugacity is f 1 Adsorption quantity V 1 In the isothermal adsorption data of (2), the isothermal adsorption curve can be determined, and isothermal adsorption data at the second temperature can be obtained.
That is, through the above steps, the relationship model between the fugacity at two different temperatures and the equivalent heat of adsorption can be determined as:
Figure RE-GDA0003118729320000131
referring to fig. 5, an apparatus for calculating methane adsorption capacity in coal at different temperatures based on a heat of adsorption theory according to an embodiment of the present invention includes:
a first processing module 51, configured to obtain isothermal adsorption data of methane in coal under a set of first temperature conditions;
the second processing module 52 is configured to input a set of adsorption pressures in isothermal adsorption data of methane in the coal under the set of first temperature conditions into an adsorption pressure and fugacity relation model, and obtain a first fugacity under the set of first temperature conditions;
a third processing module 53, configured to input a set of isothermal adsorption quantity values in isothermal adsorption data of methane in the coal under the set of first temperature conditions into an equivalent adsorption heat and isothermal adsorption quantity relation model, to obtain a set of equivalent adsorption heat values corresponding to the set of isothermal adsorption quantity values;
a fourth processing module 54, configured to input the first fugacity under the set of first temperature conditions and the set of equivalent adsorption heat values into a relationship model between the fugacity under two different temperatures and the equivalent adsorption heat, to obtain a set of second fugacity under a second temperature condition corresponding to the first fugacity under the first temperature condition;
a fifth processing module 55, configured to use the isothermal adsorption amount corresponding to any one of the first fugacity under the first temperature condition as the isothermal adsorption amount under the second temperature condition corresponding to the second fugacity under the second temperature condition corresponding to the first fugacity under the first temperature condition;
and a sixth processing module 56, configured to perform curve fitting with the second loss at the second temperature as an abscissa and the isothermal adsorption amount at the second temperature as an ordinate to obtain an isothermal adsorption curve at the second temperature, and further determine isothermal adsorption data at the second temperature according to the isothermal adsorption curve.
Since the apparatus provided in the embodiments of the present invention may be used to perform the method described in the above embodiments, the working principle and the beneficial effects thereof are similar, and therefore, the details will not be described herein, and the specific content may refer to the description of the above embodiments.
Further, the invention provides an apparatus for calculating methane adsorption capacity in coal at different temperatures based on a heat of adsorption theory in the embodiment, wherein the adsorption pressure and the fugacity relation model in the second processing module is as follows:
Ln(f/P)=Z-1-Ln(Z-B)-(A/B)Ln(1+B/Z)
wherein Z is a compression factor, dimensionless; p is the pressure, MPa; A. b is a coefficient related to the pressure P, and is specifically expanded as follows: a=c 1 P/(R 2 T 2 ),B=c 2 P/(R 2 T 2 ),
Figure BDA0003029895770000141
c 2 =0.0866RT c /P c , m=0.48+1.534ω-0.176ω 2 Wherein Pc is the critical pressure of gas, and the unit is MPa; tc is the critical temperature of the gas, and the unit is K; TR is the comparative temperature, tr=t/Tc, T representing temperature; c1 and c2 are constants; omega is an eccentric factor and is dimensionless.
Further, the device for calculating the adsorption capacity of methane in coal at different temperatures based on the adsorption heat theory provided in the embodiment of the invention, wherein the relation model of the equivalent adsorption heat and isothermal adsorption capacity in the third processing module is as follows:
q st =C 5 n+C 6
wherein q st Is equal in adsorption heat, n is isothermal adsorption quantity, C5 and C6 are q st Fitting parameters to n.
The mode of construction is the same as above, so that details will not be described here, and reference will be made to the description of the above embodiments.
Further, the device for calculating the adsorption capacity of methane in coal at different temperatures based on the adsorption heat theory provided in the embodiment of the invention, wherein a relation model between the two different temperatures in the fifth module and the equivalent adsorption heat is as follows:
Figure RE-GDA0003118729320000142
wherein q st Represents equivalent adsorption heat, and the unit is kJ/mol; t (T) 1 And T 2 Temperature is represented in K; r is a gas constant, f 1 And f 2 The loss was measured.
The mode of construction is the same as above, so that details will not be described here, and reference will be made to the description of the above embodiments.
Fig. 6 illustrates a physical schematic diagram of an electronic device, as shown in fig. 6, which may include: processor 610, communication interface (Communications Interface) 620, memory 630, and communication bus 640, wherein processor 610, communication interface 620, and memory 630 communicate with each other via communication bus 640. The processor 610 may invoke logic instructions in the memory 630 to perform a method of calculating methane adsorption in coal at different temperatures based on heat of adsorption theory, the method comprising: acquiring isothermal adsorption data of methane in coal under a group of first temperature conditions; inputting a set of adsorption pressure in isothermal adsorption data of methane in coal under the first temperature condition into an adsorption pressure and fugacity relation model to obtain a first fugacity under the first temperature condition; inputting a group of isothermal adsorption values in isothermal adsorption data of methane in the coal under the group of first temperature conditions into an isothermal adsorption heat and isothermal adsorption quantity relation model to obtain a group of equivalent adsorption heat values corresponding to the group of isothermal adsorption values; inputting the first loss degree and the set of equivalent adsorption heat values under the first temperature condition into a relation model between the loss degrees and the equivalent adsorption heat values under two different temperatures to obtain a set of second loss degrees under the second temperature condition corresponding to the first loss degrees under the first temperature condition; taking the isothermal adsorption quantity corresponding to any first loss degree under the first temperature condition as the isothermal adsorption quantity under the second temperature condition corresponding to the second loss degree under the second temperature condition corresponding to the first loss degree under the first temperature condition; and performing curve fitting by taking the second loss degree at the second temperature as an abscissa and the isothermal adsorption quantity at the second temperature as an ordinate to obtain an isothermal adsorption curve at the second temperature, and further determining isothermal adsorption data at the second temperature according to the isothermal adsorption curve.
Further, the logic instructions in the memory 630 may be implemented in the form of software functional units and stored in a computer-readable storage medium when sold or used as a stand-alone product. Based on this understanding, the technical solution of the present invention may be embodied essentially or in a part contributing to the prior art or in a part of the technical solution, in the form of a software product stored in a storage medium, including instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a random access Memory (RAM, random Access Memory), a magnetic disk, or an optical disk, or other various media capable of storing program codes.
In another aspect, the present invention also provides a computer program product comprising a computer program stored on a non-transitory computer readable storage medium, the computer program comprising program instructions which, when executed by a computer, are capable of performing the methods provided by the methods described above to perform a method of calculating methane adsorption in coal at different temperatures based on heat of adsorption theory, the method comprising: acquiring isothermal adsorption data of methane in coal under a group of first temperature conditions; inputting a set of adsorption pressure in isothermal adsorption data of methane in coal under the first temperature condition into an adsorption pressure and fugacity relation model to obtain a first fugacity under the first temperature condition; inputting a group of isothermal adsorption values in isothermal adsorption data of methane in the coal under the group of first temperature conditions into an equivalent adsorption heat and isothermal adsorption quantity relation model to obtain a group of equivalent adsorption heat values corresponding to the group of isothermal adsorption values; inputting the first loss degree and the set of equivalent adsorption heat values under the first temperature condition into a relation model between the loss degrees and equivalent adsorption heat under two different temperatures to obtain a set of second loss degrees under the second temperature condition corresponding to the first loss degree under the first temperature condition; taking the isothermal adsorption quantity corresponding to any first loss degree under the first temperature condition as the isothermal adsorption quantity under the second temperature condition corresponding to the second loss degree under the second temperature condition corresponding to the first loss degree under the first temperature condition; and performing curve fitting by taking the second loss degree at the second temperature as an abscissa and the isothermal adsorption quantity at the second temperature as an ordinate to obtain an isothermal adsorption curve at the second temperature, and determining isothermal adsorption data at the second temperature according to the isothermal adsorption curve.
In yet another aspect, the present invention also provides a non-transitory computer readable storage medium having stored thereon a computer program which when executed by a processor is implemented to perform the methods provided above to perform the calculation of methane adsorption amounts in coal at different temperatures based on heat of adsorption theory, the method comprising: acquiring isothermal adsorption data of methane in coal under a group of first temperature conditions; inputting a set of adsorption pressure in isothermal adsorption data of methane in coal under the first temperature condition into an adsorption pressure and fugacity relation model to obtain a first fugacity under the first temperature condition; inputting a set of isothermal adsorption values in isothermal adsorption data of methane in the coal under the first temperature condition into an equivalent adsorption heat and isothermal adsorption quantity relation model to obtain a set of equivalent adsorption heat values corresponding to the set of isothermal adsorption values; inputting the first loss degree and the set of equivalent adsorption heat values under the first temperature condition into a relation model between the loss degrees and the equivalent adsorption heat under two different temperatures to obtain a set of second loss degrees under the second temperature condition corresponding to the first loss degrees under the first temperature condition; taking the isothermal adsorption quantity corresponding to any first loss degree under the first temperature condition as the isothermal adsorption quantity under the second temperature condition corresponding to the second loss degree under the second temperature condition corresponding to the first loss degree under the first temperature condition; and performing curve fitting by taking the second loss degree at the second temperature as an abscissa and the isothermal adsorption quantity at the second temperature as an ordinate to obtain an isothermal adsorption curve at the second temperature, and determining isothermal adsorption data at the second temperature according to the isothermal adsorption curve.
The apparatus embodiments described above are merely illustrative, wherein the elements illustrated as separate components may or may not be physically separate, and the components shown as elements may or may not be physical elements, may be located in one place, or may be distributed over a plurality of network elements. Some or all of the modules can be selected according to actual needs to achieve the purpose of the scheme of the embodiment. Those of ordinary skill in the art will understand and practice the present invention without undue burden.
From the above description of the embodiments, it will be apparent to those skilled in the art that the embodiments may be implemented by means of software plus necessary general hardware platforms, or of course may be implemented by means of hardware. Based on this understanding, the foregoing technical solution may be embodied essentially or in a part contributing to the prior art in the form of a software product, which may be stored in a computer readable storage medium, such as a ROM/RAM, a magnetic disk, an optical disk, etc., including several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to perform the method described in the various embodiments or some parts of the embodiments.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solution of the present invention, and are not limiting; although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced equivalently; such modifications and substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims (6)

1. A method for calculating methane adsorption capacity in coal at different temperatures based on a heat of adsorption theory, comprising the steps of:
acquiring isothermal adsorption data of methane in coal under a group of first temperature conditions;
inputting a set of adsorption pressure in isothermal adsorption data of methane in coal under the first temperature condition into an adsorption pressure and fugacity relation model to obtain a first fugacity under the first temperature condition;
inputting a set of isothermal adsorption values in isothermal adsorption data of methane in the coal under the first temperature condition into an equivalent adsorption heat and isothermal adsorption quantity relation model to obtain a set of equivalent adsorption heat values corresponding to the set of isothermal adsorption values; wherein, the relation model of the equivalent adsorption heat and isothermal adsorption amount is as follows:
q st =C 5 n+C 6
wherein q st Is equal in adsorption heat, n is isothermal adsorption quantity, C5 and C6 are q st Fitting parameters with n;
inputting the first loss degree and the set of equivalent adsorption heat values under the first temperature condition into a relation model between the loss degrees and the equivalent adsorption heat under two different temperatures to obtain a set of second loss degrees under the second temperature condition corresponding to the first loss degrees under the first temperature condition; the relation model between the fugacity at the two different temperatures and the equivalent adsorption heat is as follows:
Figure FDA0004044663940000011
wherein q st Represents equivalent adsorption heat, and the unit is kJ/mol; t (T) 1 And T 2 Temperature is represented in K; r is a gas constant, f 1 And f 2 Is the fugacity;
taking the isothermal adsorption quantity corresponding to any first loss degree under the first temperature condition as the isothermal adsorption quantity under the second temperature condition corresponding to the second loss degree under the second temperature condition corresponding to the first loss degree under the first temperature condition;
and performing curve fitting by taking the second loss degree at the second temperature as an abscissa and the isothermal adsorption quantity at the second temperature as an ordinate to obtain an isothermal adsorption curve at the second temperature, and determining isothermal adsorption data at the second temperature according to the isothermal adsorption curve.
2. The method for calculating the adsorption amount of methane in coal at different temperatures based on the adsorption heat theory according to claim 1, wherein the relation model of adsorption pressure and fugacity is as follows:
Ln(f/P)=Z-1-Ln(Z-B)-(A/B)Ln(1+B/Z)
wherein Z is a compression factor, dimensionless; p is the pressure, MPa; A. b is a coefficient related to the pressure P, and is specifically expanded as follows: a=c 1 P/(R 2 T 2 ),B=c 2 P/(R 2 T 2 ),
Figure FDA0004044663940000021
c 2 =0.0866RT c /P c ,m=0.48+1.534ω-0.176ω 2 Wherein Pc is the critical pressure of gas, and the unit is MPa; tc is the critical temperature of the gas, and the unit is K; TR is the comparative temperature, tr=t/Tc, T representing temperature; c1 and c2 are constants; omega is an eccentric factor and is dimensionless.
3. An apparatus for calculating methane adsorption capacity in coal at different temperatures based on a heat of adsorption theory, comprising:
the first processing module is used for acquiring isothermal adsorption data of methane in the coal under a group of first temperature conditions;
the second processing module is used for inputting a set of adsorption pressure in isothermal adsorption data of methane in the coal under the first temperature condition into an adsorption pressure and fugacity relation model to obtain a first fugacity under the first temperature condition;
the third processing module is used for inputting a group of isothermal adsorption quantity values in isothermal adsorption data of methane in the coal under the group of first temperature conditions into an equivalent adsorption heat and isothermal adsorption quantity relation model to obtain a group of equivalent adsorption heat values corresponding to the group of isothermal adsorption quantity values; the relationship model of the equivalent adsorption heat and isothermal adsorption amount in the third processing module is as follows:
q st =C 5 n+C 6
wherein q st Is equal in adsorption heat, n is isothermal adsorption quantity, C5 and C6 are q st Fitting parameters with n;
the fourth processing module is used for inputting the first loss degree and the set of equivalent adsorption heat values under the first temperature condition into a relation model between the loss degrees and the equivalent adsorption heat under two different temperatures to obtain a set of second loss degrees under the second temperature condition corresponding to the first loss degrees under the first temperature condition; the relation model between the fugacity at the two different temperatures and the equivalent adsorption heat is as follows:
Figure FDA0004044663940000031
wherein q st Represents equivalent adsorption heat, and the unit is kJ/mol; t (T) 1 And T 2 Temperature is represented in K; r is a gas constant, f 1 And f 2 Is the fugacity;
a fifth processing module, configured to use the isothermal adsorption amount corresponding to any first loss under the first temperature condition as the isothermal adsorption amount under the second temperature condition corresponding to the second loss under the second temperature condition corresponding to the first loss under the first temperature condition;
and the sixth processing module is used for performing curve fitting by taking the second loss degree at the second temperature as an abscissa and the isothermal adsorption quantity at the second temperature as an ordinate to obtain an isothermal adsorption curve at the second temperature, and further determining isothermal adsorption data at the second temperature according to the isothermal adsorption curve.
4. The apparatus for calculating methane adsorption capacity in coal at different temperatures based on heat of adsorption theory of claim 3, wherein the adsorption pressure versus fugacity model in the second treatment module is:
Ln(f/P)=Z-1-Ln(Z-B)-(A/B)Ln(1+B/Z)
wherein Z is a compression factor, dimensionless; p is the pressure, MPa; A. b is a coefficient related to the pressure P, and is specifically expanded as follows: a=c 1 P/(R 2 T 2 ),B=c 2 P/(R 2 T 2 ),
Figure FDA0004044663940000032
c 2 =0.0866RT c /P c ,m=0.48+1.534ω-0.176ω 2 Wherein Pc is the critical pressure of gas, and the unit is MPa; tc is the critical temperature of the gas, and the unit is K; TR is the comparative temperature, tr=t/Tc, T representing temperature; c1 and c2 are constants; omega is an eccentric factor and is dimensionless.
5. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the processor implements the steps of the method according to any one of claims 1 or 2 when the program is executed.
6. A non-transitory computer readable storage medium, on which a computer program is stored, characterized in that the computer program, when being executed by a processor, implements the steps of the method according to any one of claims 1 or 2.
CN202110426819.4A 2021-04-20 2021-04-20 Method and device for calculating methane adsorption quantity in coal at different temperatures Active CN113237917B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110426819.4A CN113237917B (en) 2021-04-20 2021-04-20 Method and device for calculating methane adsorption quantity in coal at different temperatures

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110426819.4A CN113237917B (en) 2021-04-20 2021-04-20 Method and device for calculating methane adsorption quantity in coal at different temperatures

Publications (2)

Publication Number Publication Date
CN113237917A CN113237917A (en) 2021-08-10
CN113237917B true CN113237917B (en) 2023-05-16

Family

ID=77128636

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110426819.4A Active CN113237917B (en) 2021-04-20 2021-04-20 Method and device for calculating methane adsorption quantity in coal at different temperatures

Country Status (1)

Country Link
CN (1) CN113237917B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113392567B (en) * 2021-08-16 2021-11-09 中国矿业大学(北京) Method and system for predicting coal bed gas of double-pore coal body

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6860147B2 (en) * 2002-09-30 2005-03-01 Alberta Research Council Inc. Process for predicting porosity and permeability of a coal bed
KR101782294B1 (en) * 2017-04-20 2017-10-23 한국가스공사 Modeling method for gas production of cbm reservoir rocks
CN108120657B (en) * 2017-12-08 2019-12-17 中国矿业大学 device for automatically and synchronously measuring high-pressure gas adsorption capacity and adsorption heat and measuring method
CN108120656B (en) * 2017-12-20 2020-05-12 辽宁工程技术大学 System and method for measuring temperature and heat in coal bed gas adsorption and desorption processes
US20220341861A1 (en) * 2019-06-12 2022-10-27 Texas Tech University System Thermodynamic Formulation for Langmuir Adsorption Isotherms
CN111896421B (en) * 2020-07-29 2021-05-28 西南石油大学 Method for calculating true adsorption capacity of methane in shale based on adsorption potential theory
CN112067783A (en) * 2020-09-03 2020-12-11 中国石油大学(北京) Method and device for determining gas-water adsorption quantity of coal reservoir

Also Published As

Publication number Publication date
CN113237917A (en) 2021-08-10

Similar Documents

Publication Publication Date Title
Otero-de-la-Roza et al. Gibbs2: A new version of the quasi-harmonic model code. I. Robust treatment of the static data
Herslund et al. Phase equilibrium modeling of gas hydrate systems for CO2 capture
Seo et al. Multiple-phase hydrate equilibria of the ternary carbon dioxide, methane, and water mixtures
Shen et al. Modeling thermodynamic derivative properties of ionic liquids with ePC-SAFT
CN113237917B (en) Method and device for calculating methane adsorption quantity in coal at different temperatures
CN110534161B (en) Adsorbate gas adsorption phase density model construction and absolute adsorption amount calculation method
Bahadori et al. Rapid estimation of equilibrium water dew point of natural gas in TEG dehydration systems
CN107153755B (en) Solving method for shale gas well numerical simulation
CN108593493B (en) Method and device for determining gas content of tight reservoir
CN111507537A (en) Method and device for predicting reserves and water invasion of water-drive gas reservoir
CN112084726B (en) Shale gas adsorption desorption diffusion kinetic model and interpretation method
Zeng et al. A time filter method for solving the double-diffusive natural convection model
CN107239648B (en) Shale gas well yield composition determination method and device
Zhang et al. One‐level and multilevel space‐time finite element method for the viscoelastic Kelvin‐Voigt model
JP6964365B2 (en) Calculation and analysis method of hydrate saturation in experimental and numerical simulation processes
Schlijper Simulation of compositional processes: The use of pseudocomponents in equation-of-state calculations
Fernandez-Prini et al. Aqueous solubility of volatile nonelectrolytes
CN112014287A (en) Method and device for determining adsorption mode of adsorbate gas
Pandey et al. Impact of bioconversion on matrix strain response of coal reservoirs: Part 2-Reservoir insights
CN109254113B (en) Method and system for measuring adsorption capacity of multi-component mixed gas in gas-solid adsorption process
CN114251075B (en) Reservoir exploitation scheme determination method, device and equipment based on multi-target parameters
CN115169253A (en) Thermodynamic two-temperature model non-equilibrium energy system improvement method
Thakre et al. Proposing ab initio assisted lattice distortion theory for phase equilibrium: Pure and mixed refrigerant gas hydrates
He et al. Error estimates of fully discrete finite element solutions for the 2D Cahn–Hilliard equation with infinite time horizon
CN112131513B (en) Fluid unbalanced phase calculation method and system

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant