CN112613174B - Evaluation method for methane adsorption capacity of shale considering multiple adsorption mechanisms - Google Patents
Evaluation method for methane adsorption capacity of shale considering multiple adsorption mechanisms Download PDFInfo
- Publication number
- CN112613174B CN112613174B CN202011506504.2A CN202011506504A CN112613174B CN 112613174 B CN112613174 B CN 112613174B CN 202011506504 A CN202011506504 A CN 202011506504A CN 112613174 B CN112613174 B CN 112613174B
- Authority
- CN
- China
- Prior art keywords
- adsorption
- methane
- shale
- isothermal
- model
- 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
Links
- 238000001179 sorption measurement Methods 0.000 title claims abstract description 219
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 title claims abstract description 214
- 230000007246 mechanism Effects 0.000 title claims abstract description 38
- 238000011156 evaluation Methods 0.000 title claims description 13
- 238000000034 method Methods 0.000 claims abstract description 19
- 239000011148 porous material Substances 0.000 claims abstract description 17
- 239000002052 molecular layer Substances 0.000 claims abstract description 11
- 238000002474 experimental method Methods 0.000 claims description 15
- 238000012360 testing method Methods 0.000 claims description 12
- 238000009835 boiling Methods 0.000 claims description 9
- 229920006395 saturated elastomer Polymers 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 6
- 238000005303 weighing Methods 0.000 claims description 6
- 238000005516 engineering process Methods 0.000 claims description 4
- 239000000126 substance Substances 0.000 claims description 4
- 238000005553 drilling Methods 0.000 claims description 3
- 238000001035 drying Methods 0.000 claims description 3
- 239000012535 impurity Substances 0.000 claims description 3
- 238000004377 microelectronic Methods 0.000 claims description 3
- 239000000843 powder Substances 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 239000010419 fine particle Substances 0.000 claims description 2
- 238000012216 screening Methods 0.000 claims description 2
- 238000011161 development Methods 0.000 abstract description 8
- 239000003345 natural gas Substances 0.000 abstract description 5
- 238000012512 characterization method Methods 0.000 abstract description 4
- 239000003245 coal Substances 0.000 abstract description 4
- 239000010410 layer Substances 0.000 description 16
- 239000007789 gas Substances 0.000 description 12
- 238000011160 research Methods 0.000 description 5
- 230000007547 defect Effects 0.000 description 4
- 238000000329 molecular dynamics simulation Methods 0.000 description 4
- 239000002356 single layer Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Solid-Sorbent Or Filter-Aiding Compositions (AREA)
Abstract
本发明公开了一种考虑多重吸附机制的页岩吸附甲烷能力评估方法，基于不同尺寸纳米孔道中吸附机制的差异，以及不同吸附机制下吸附相甲烷密度的不同，对微孔充填形式吸附和多分子层形式吸附进行联合表征，所建立的DA‑TBET等温吸附模型具有清晰的物理背景，克服了当前吸附模型的不足，评估所得页岩吸附甲烷能力结果更加合理准确。此外，本发明除了用于评估页岩吸附甲烷能力，还可用于评估煤层样品吸附甲烷能力，在非常规天然气勘探开发中具有重要应用价值。
The invention discloses a method for evaluating the methane adsorption capacity of shale considering multiple adsorption mechanisms. Based on the differences in adsorption mechanisms in nano-pore channels of different sizes and the difference in adsorption phase methane density under different adsorption mechanisms, the adsorption and multi-adsorption in the form of micropore filling are evaluated. Combined characterization of adsorption in molecular layer form, the established DA-TBET adsorption isotherm model has a clear physical background, overcomes the shortcomings of the current adsorption model, and evaluates the methane adsorption capacity of the obtained shale The results are more reasonable and accurate. In addition, the present invention can be used for evaluating the methane adsorption capacity of shale, and can also be used for evaluating the methane adsorption capacity of coal seam samples, which has important application value in the exploration and development of unconventional natural gas.
Description
Technical Field
The invention relates to the technical field of unconventional natural gas exploration and development, in particular to a shale methane adsorption capacity evaluation method considering multiple adsorption mechanisms.
Background
The reasonable and accurate evaluation of the methane adsorption capacity of the shale is an important basis for realizing economic and efficient development of the shale gas reservoir, and plays a significant role in gas reservoir development. Through reasonably and accurately evaluating the methane adsorption capacity of the shale, the gas reservoir geological reserve can be calculated, and the method has important guiding significance for making and adjusting a shale gas reservoir development scheme, and related researches are widely concerned.
The indoor isothermal adsorption experiment is an important method for researching the methane adsorption capacity of the shale, and researchers develop a large amount of experimental improvement innovation and theoretical derivation work based on the experiment, wherein the research on the isothermal adsorption model representation of the methane adsorption capacity of the shale is particularly hot. According to different assumed adsorption mechanisms, the isothermal adsorption model mainly comprises a monolayer model, a polymolecular layer model and a micropore filling model. At present, although there are many isothermal adsorption models for methane in shale, no model considers the difference of adsorption phase methane density in pore channels with different sizes, and the adsorption phase methane density in different isothermal adsorption models is considered to be the same. According to the conclusion of molecular dynamics simulation research, methane molecules in micropores are influenced by the adsorption potential of the wall surface of a pore passage in an overlapping manner, so that the adsorption phase molecules are arranged more closely and the density of the adsorption phase molecules is higher, which is not distinguished in the conventional research on an isothermal adsorption model. In addition, although researchers consider microporous adsorption mode as filling adsorption mode and monolayer adsorption mode in mesopores and macropores and provide corresponding combined models (such as DA-Langmuir model and DA-LF model), the researchers find that methane is actually adsorbed in the mesopores as a polymolecular layer based on molecular dynamics simulation. Therefore, corresponding research work is necessary to be carried out, an isothermal adsorption model for effectively representing multiple adsorption mechanisms of shale gas is established, and support is provided for reasonably and accurately evaluating methane adsorption capacity of shale.
Disclosure of Invention
The invention mainly overcomes the defects in the prior art, and aims to provide the shale methane adsorption capacity evaluation method considering multiple adsorption mechanisms.
In order to achieve the technical purpose, the invention adopts the following technical scheme:
an evaluation method for methane adsorption capacity of shale considering multiple adsorption mechanisms is characterized in that the calculation method comprises the following steps:
s1: drilling a shale core, crushing the core into powder in a laboratory, and screening out shale fine particles with the size of 60-80 meshes by using a screen;
s2: drying the shale sample prepared in the step in a thermostat at 80 ℃, and performing mass weighing record once by using a micro electronic balance every 12 hours until the mass change of the sample obtained by two adjacent times of weighing is less than or equal to 2%;
s3: loading the sample into a sample cylinder of a volumetric or gravimetric isothermal adsorption instrument, continuously vacuumizing for 24 hours, and further removing impurity adsorption gas and water vapor in the sample;
s4: carrying out isothermal adsorption experiments on the shale samples by using a volumetric method or gravimetric method isothermal adsorption instrument at least 3 different temperatures, and recording the change data of the excess adsorption quantity along with the pressure;
s5: fitting the experimental data by adopting a DA-TBET isothermal adsorption model considering a multiple adsorption mechanism to obtain isothermal adsorption model parameters at different temperatures;
s6: and calculating the adsorption capacity curve of the shale to the methane at different temperatures and pressures by using the isothermal adsorption model parameters at different temperatures obtained in the steps.
Further, the method for crushing the shale sample in the step S1 can refer to GB/T19560-2008 standard.
Further, the DA-TBET isothermal adsorption model expression considering the multiple adsorption mechanism established in step S5 is as follows:
in the formula, VexFor actually measuring the excess adsorption capacity, cm3/g；ρad1Adsorption phase methane density in cm in micropore filling mode3/g；ρad2Adsorption phase methane density in cm in the form of medium-large pore polymolecular layer3/g；ρgIs free phase methane density, cm3/g；VwThe volume of the adsorption space in cm is filled in the form of micropores in the sample3(ii)/g; r is a general gas constant, and the value here is 8.314J/(mol.K); t is adsorption experiment test temperature, K; e is the adsorption characteristic energy, J/mol; p is the experimental test pressure, MPa; p is a radical of0Is the methane virtual saturated vapor pressure, MPa; q is a parameter for representing the heterogeneity of the shale adsorption pores; vmIs the maximum adsorption capacity of the first molecular layer in the mesopores and macropores, cm3(ii)/g; c is a constant related to the heat of adsorption; n represents the number of adsorption layers; free phase methane density ρgCan be directly obtained by searching in chemical databases of national institute of standards and technology.
Further, the key parameters in the DA-TBET isothermal adsorption model considering the multiple adsorption mechanism established in step S5 are selected as follows:
adsorption phase methane density rho for micropore filling formad1Because of the superposition of adsorption potentials in narrow pore canals, methane molecules are closely arranged, and the density of the methane molecules is considered to be equal to the density of a boiling point, and the value is 0.425g/cm3；
Adsorption phase methane density rho for multi-molecular layer form in mesopore and macroporead2Using an adsorbed phase methane densitometer as proposed by Ozawa et al (1976)Calculating an empirical formula, wherein the expression is as follows:
ρad2＝ρbexp[-0.0025(T-Tb)]
in the formula, ρad2Adsorption phase methane density in cm in the form of medium-large pore polymolecular layer3/g；ρbThe boiling point density of methane is 0.425g/cm3(ii) a T is adsorption experiment test temperature, K; t isbThe boiling point temperature of methane under the standard atmospheric pressure is 111.7K;
for the number n of adsorbed layers of the polymolecular layer in the mesopores and macropores, the study conclusion of Zhongshang et al (shale gas adsorption mechanism study based on modified BET model, 2018) and the study conclusion of the applicant based on molecular dynamics simulation should be combined, and the number of adsorbed layers of the polymolecular layer in the shale should be 2.
Further, the virtual saturated vapor pressure calculation formula in the DA-TBET isothermal adsorption model considering the multiple adsorption mechanism established in step S5 is as follows:
in the formula, p0Is the methane virtual saturated vapor pressure, MPa; p is a radical ofcThe critical pressure of methane is 4.5992 MPa; t iscThe critical temperature of methane is 190.564K; t is adsorption experiment test temperature, K; k is a parameter characterizing the adsorption system.
Further, when the DA-TBET isothermal adsorption model is adopted to fit the experimental data at different temperatures in step S5, it is required to satisfy the parameter V in the DA-TBET modelwE, k and q remain the same.
Further, the adsorption capacity curve expression of the shale to methane at different temperatures and pressures in the step S6 is as follows:
in the formula, VabIs the real adsorption capacity of the shale to the methane, cm3(ii)/g; the remaining parameters are the values fitted in step S5.
According to the shale methane adsorption capacity evaluation method considering multiple adsorption mechanisms, based on the difference of adsorption mechanisms in nanopores with different sizes and the difference of adsorption phase methane densities under different adsorption mechanisms, the combined characterization is carried out on the adsorption in the micropore filling form and the adsorption in the multi-molecular-layer form, the established DA-TBET isothermal adsorption model has a clear physical background, the defects of the current adsorption model are overcome, and the evaluation result of the shale methane adsorption capacity is more reasonable and accurate. In addition, the method can be used for evaluating the methane adsorption capacity of the shale and the methane adsorption capacity of the coal bed sample, and has important application value in unconventional natural gas exploration and development.
Has the advantages that:
compared with the prior art, the invention has the following beneficial effects:
based on the difference of adsorption mechanisms in nanopores with different sizes and the difference of adsorption phase methane densities under different adsorption mechanisms, the combined characterization of micropore filling type adsorption and multi-molecular-layer type adsorption is carried out, the established DA-TBET isothermal adsorption model has a clear physical background, the defects of the current adsorption model are overcome, and the evaluation result of the methane adsorption capacity of the shale is more reasonable and accurate. In addition, the method can be used for evaluating the methane adsorption capacity of the shale and the methane adsorption capacity of the coal bed sample, and has important application value in unconventional natural gas exploration and development.
Drawings
FIG. 1 is measured isothermal adsorption data for shale samples 4-08 at different temperatures;
FIG. 2 is a comparison graph of measured isothermal adsorption data of shale samples 4-08 at different temperatures and fitting results of a DA-TBET model;
FIG. 3 is a graph of the methane adsorption capacity of shale samples 4-08 evaluated at different temperatures based on the DA-TBET model.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail with reference to the following embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
Example (b):
a shale methane adsorption capacity evaluation method considering multiple adsorption mechanisms comprises the following specific calculation methods:
s1: drilling a shale core, crushing the core into powder according to GB/T19560-;
s2: drying the shale sample prepared in the step in a thermostat at 80 ℃, and performing mass weighing record once by using a micro electronic balance every 12 hours until the mass change of the sample obtained by two adjacent times of weighing is less than or equal to 2%;
s3: loading the sample into a sample cylinder of a volumetric or gravimetric isothermal adsorption instrument, continuously vacuumizing for 24 hours, and further removing impurity adsorption gas and water vapor in the sample;
s4: carrying out isothermal adsorption experiments on the shale samples by using a volumetric method or gravimetric method isothermal adsorption instrument at least 3 different temperatures, and recording the change data of the excess adsorption quantity along with the pressure;
s5: fitting the experimental data by adopting a DA-TBET isothermal adsorption model considering a multiple adsorption mechanism to obtain isothermal adsorption model parameters at different temperatures;
the DA-TBET isothermal adsorption model expression taking the multiple adsorption mechanism into consideration established in the step S5 is as follows:
in the formula, VexFor actually measuring the excess adsorption capacity, cm3/g；ρad1Adsorption phase methane density in cm in micropore filling mode3/g；ρad2Adsorption phase methane density in cm in the form of medium-large pore polymolecular layer3/g；ρgIs a free phaseDensity of methane in cm3/g；VwThe volume of the adsorption space in cm is filled in the form of micropores in the sample3(ii)/g; r is a general gas constant, and the value here is 8.314J/(mol.K); t is adsorption experiment test temperature, K; e is the adsorption characteristic energy, J/mol; p is the experimental test pressure, MPa; p is a radical of0Is the methane virtual saturated vapor pressure, MPa; q is a parameter for representing the heterogeneity of the shale adsorption pores; vmIs the maximum adsorption capacity of the first molecular layer in the mesopores and macropores, cm3(ii)/g; c is a constant related to the heat of adsorption; n represents the number of adsorption layers; free phase methane density ρgCan be directly obtained by searching in chemical databases of national institute of standards and technology.
The key parameters in the DA-TBET isothermal adsorption model considering the multiple adsorption mechanism established in the step S5 are selected as follows:
adsorption phase methane density rho for micropore filling formad1Because of the superposition of adsorption potentials in narrow pore canals, methane molecules are closely arranged, and the density of the methane molecules is considered to be equal to the density of a boiling point, and the value is 0.425g/cm3；
Adsorption phase methane density rho for multi-molecular layer form in mesopore and macroporead2An empirical formula was calculated using the adsorption phase methane density proposed by Ozawa et al (1976) and expressed as:
ρad2＝ρbexp[-0.0025(T-Tb)]
in the formula, ρad2Adsorption phase methane density in cm in the form of medium-large pore polymolecular layer3/g；ρbThe boiling point density of methane is 0.425g/cm3(ii) a T is adsorption experiment test temperature, K; t isbThe boiling point temperature of methane under the standard atmospheric pressure is 111.7K;
for the number n of adsorbed layers of the polymolecular layer in the mesopores and macropores, the study conclusion of Zhongshang et al (shale gas adsorption mechanism study based on modified BET model, 2018) and the study conclusion of the applicant based on molecular dynamics simulation should be combined, and the number of adsorbed layers of the polymolecular layer in the shale should be 2.
The virtual saturated vapor pressure calculation formula in the DA-TBET isothermal adsorption model considering the multiple adsorption mechanism established in step S5 is as follows:
in the formula, p0Is the methane virtual saturated vapor pressure, MPa; p is a radical ofcThe critical pressure of methane is 4.5992 MPa; t iscThe critical temperature of methane is 190.564K; t is adsorption experiment test temperature, K; k is a parameter characterizing the adsorption system.
When the DA-TBET isothermal adsorption model is adopted to fit the experimental data at different temperatures in the step S5, the parameter V in the DA-TBET model needs to be satisfiedwE, k and q remain the same.
S6: calculating the adsorption capacity curve of the shale to the methane at different temperatures and pressures by using the isothermal adsorption model parameters at different temperatures obtained in the step;
the expression of the adsorption capacity curve of the shale to methane at different temperatures and pressures in the step S6 is as follows:
in the formula, VabIs the real adsorption capacity of the shale to the methane, cm3(ii)/g; the remaining parameters are the values fitted in step S5.
Example 1:
the data used in this example are from 4-08 adsorption data for shale samples taken from the Shirai reef dam in the open literature (Tian et al, 2016) and the isothermal adsorption data obtained using a volumetric instrumentation at 35.4 deg.C, 50.4 deg.C and 65.4 deg.C, respectively, are shown in FIG. 1 and Table 1.
Free phase methane densities at different temperature and pressure conditions were obtained from queries in the national institute of standards and technology chemical database, as shown in table 1.
Adopting the established DA-TBET isothermal adsorption model to perform isothermal adsorption experiments on the samples measured at different temperaturesFitting data, calculating model parameters at different temperatures and fitting R2As shown in table 2. And substituting the parameters obtained by fitting into a DA-TBET model, and comparing the calculation result with the actually measured experimental data, as shown in figure 2. As can be seen from Table 2, the fitting R at different temperatures2All exceed 0.999, which shows that the DA-TBET model has good fitting effect. It can also be seen from FIG. 2 that the model fitting results are in good agreement with the measured isothermal adsorption experimental data. Based on the fitted model parameters, the methane adsorption capacity curve of the shale sample at different temperatures can be calculated, as shown in fig. 3.
In addition, by the established DA-TBET isothermal adsorption model, the total adsorption capacity of the shale to methane can be evaluated, and the methane gas quantity adsorbed by different adsorption modes in the shale pore canal can be respectively calculated by the first half part (micropore filling adsorption) and the second half part (multi-molecular layer adsorption) of the model, so that the method has important significance for further clarifying the adsorption capacity of the shale in different pore canals.
TABLE 1 measured isothermal adsorption data and free phase methane Density at different temperatures and pressures
TABLE 2 model parameters and fitting R calculated at different temperatures2
| T | Vw | E | k | q | Vm | C | n | ρad1 | ρad2 | R2 |
| 308.55 | 0.00114 | 13721.5 | 4.71078 | 3.3193 | 2.10586 | 9.49687 | 2.0 | 0.425 | 0.25822 | 0.99926 |
| 323.55 | 0.00114 | 13721.5 | 4.71078 | 3.3193 | 1.98224 | 11.11999 | 2.0 | 0.425 | 0.24872 | 0.99968 |
| 338.55 | 0.00114 | 13721.5 | 4.71078 | 3.3193 | 1.96559 | 11.34600 | 2.0 | 0.425 | 0.23957 | 0.99903 |
According to the shale methane adsorption capacity evaluation method considering multiple adsorption mechanisms, based on the difference of adsorption mechanisms in nanopores with different sizes and the difference of adsorption phase methane densities under different adsorption mechanisms, the combined characterization is carried out on the adsorption in the micropore filling form and the adsorption in the multi-molecular-layer form, the established DA-TBET isothermal adsorption model has a clear physical background, the defects of the current adsorption model are overcome, and the evaluation result of the shale methane adsorption capacity is more reasonable and accurate. In addition, the method can be used for evaluating the methane adsorption capacity of the shale and the methane adsorption capacity of the coal bed sample, and has important application value in unconventional natural gas exploration and development.
Although the present invention has been described with reference to the above embodiments, it should be understood that the present invention is not limited to the above embodiments, and those skilled in the art can make various changes and modifications without departing from the scope of the present invention.
Claims (5)
1. A shale methane adsorption capacity evaluation method considering multiple adsorption mechanisms is characterized by comprising the following steps:
s1: drilling a shale core, crushing the core into powder in a laboratory, and screening out shale fine particles with the size of 60-80 meshes by using a screen;
s2: drying the shale sample prepared in the step in a thermostat at 80 ℃, and performing mass weighing record once by using a micro electronic balance every 12 hours until the mass change of the sample obtained by two adjacent times of weighing is less than or equal to 2%;
s3: loading the sample into a sample cylinder of a volumetric or gravimetric isothermal adsorption instrument, continuously vacuumizing for 24 hours, and further removing impurity adsorption gas and water vapor in the sample;
s4: carrying out isothermal adsorption experiments on the shale samples by using a volumetric method or gravimetric method isothermal adsorption instrument at least 3 different temperatures, and recording the change data of the excess adsorption quantity along with the pressure;
s5: fitting the experimental data by adopting a DA-TBET isothermal adsorption model considering a multiple adsorption mechanism to obtain isothermal adsorption model parameters at different temperatures; the DA-TBET isothermal adsorption model expression taking the multiple adsorption mechanism into consideration established in the step S5 is as follows:
in the formula, VexFor actually measuring the excess adsorption capacity, cm3/g；ρad1Adsorption phase methane density in cm in micropore filling mode3/g；ρad2Adsorption phase methane density in cm in the form of medium-large pore polymolecular layer3/g；ρgIs free phase methane density, cm3/g；VwThe volume of the adsorption space in cm is filled in the form of micropores in the sample3(ii)/g; r is a general gas constant, and the value here is 8.314J/(mol.K); t is adsorption experiment test temperature, K; e is adsorptionCharacteristic energy, J/mol; p is the experimental test pressure, MPa; p is a radical of0Is the methane virtual saturated vapor pressure, MPa; q is a parameter for representing the heterogeneity of the shale adsorption pores; vmIs the maximum adsorption capacity of the first molecular layer in the mesopores and macropores, cm3(ii)/g; c is a constant related to the heat of adsorption; n represents the number of adsorption layers; free phase methane density ρgCan be directly obtained by inquiring in a chemical database of the national institute of standards and technology;
s6: calculating the adsorption capacity curve of the shale to the methane at different temperatures and pressures by using the isothermal adsorption model parameters at different temperatures obtained in the step;
the expression of the adsorption capacity curve of the shale to methane at different temperatures and pressures in the step S6 is as follows:
in the formula, VabIs the real adsorption capacity of the shale to the methane, cm3(ii)/g; the remaining parameters are the values fitted in step S5.
2. The method for evaluating the methane adsorption capacity of shale considering multiple adsorption mechanisms as claimed in claim 1, wherein the method for crushing the shale sample in step S1 is according to GB/T19560-2008 standard.
3. The method for evaluating the methane adsorption capacity of shale considering multiple adsorption mechanisms according to claim 1, wherein the key parameters in the DA-TBET isothermal adsorption model considering multiple adsorption mechanisms established in step S5 are selected as follows:
adsorption phase methane density rho for micropore filling formad1Because of the superposition of adsorption potentials in narrow pore canals, methane molecules are closely arranged, and the density of the methane molecules is considered to be equal to the density of a boiling point, and the value is 0.425g/cm3；
Adsorption phase methane density rho for multi-molecular layer form in mesopore and macroporead2Which isThe expression is as follows:
ρad2＝ρbexp[-0.0025(T-Tb)]
in the formula, ρad2Adsorption phase methane density in cm in the form of medium-large pore polymolecular layer3/g；ρbThe boiling point density of methane is 0.425g/cm3(ii) a T is adsorption experiment test temperature, K; t isbThe boiling point temperature of methane under the standard atmospheric pressure is 111.7K;
and for the number n of the adsorption layers of the polymolecular layers in the mesopores and the macropores, the number of the adsorption layers of the polymolecular layers in the shale is 2.
4. The method for evaluating the methane adsorption capacity of shale in consideration of multiple adsorption mechanisms according to claim 1, wherein the virtual saturated vapor pressure calculation formula in the DA-TBET isothermal adsorption model in consideration of multiple adsorption mechanisms established in the step S5 is as follows:
in the formula, p0Is the methane virtual saturated vapor pressure, MPa; p is a radical ofcThe critical pressure of methane is 4.5992 MPa; t iscThe critical temperature of methane is 190.564K; t is adsorption experiment test temperature, K; k is a parameter characterizing the adsorption system.
5. The method for evaluating the methane adsorption capacity of shale considering multiple adsorption mechanisms as claimed in claim 1, wherein the parameter V in the DA-TBET model is satisfied when the DA-TBET isothermal adsorption model is used to fit experimental data at different temperatures in step S5wE, k and q remain the same.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011506504.2A CN112613174B (en) | 2020-12-18 | 2020-12-18 | Evaluation method for methane adsorption capacity of shale considering multiple adsorption mechanisms |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011506504.2A CN112613174B (en) | 2020-12-18 | 2020-12-18 | Evaluation method for methane adsorption capacity of shale considering multiple adsorption mechanisms |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112613174A CN112613174A (en) | 2021-04-06 |
| CN112613174B true CN112613174B (en) | 2022-03-15 |
Family
ID=75240630
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011506504.2A Active CN112613174B (en) | 2020-12-18 | 2020-12-18 | Evaluation method for methane adsorption capacity of shale considering multiple adsorption mechanisms |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112613174B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113670960B (en) * | 2021-07-23 | 2021-12-28 | 西南石油大学 | A real shale gas adsorption prediction method based on molecular simulation |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107106665A (en) * | 2014-06-06 | 2017-08-29 | 纪念斯隆-凯特琳癌症中心 | Target Chimeric antigen receptor of mesothelin and application thereof |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106251232B (en) * | 2016-08-12 | 2019-09-06 | 中国石油天然气股份有限公司 | Method and device for determining gas content of shale |
| CN111896421B (en) * | 2020-07-29 | 2021-05-28 | 西南石油大学 | A method for calculating the real adsorption amount of methane in shale based on adsorption potential theory |
-
2020
- 2020-12-18 CN CN202011506504.2A patent/CN112613174B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107106665A (en) * | 2014-06-06 | 2017-08-29 | 纪念斯隆-凯特琳癌症中心 | Target Chimeric antigen receptor of mesothelin and application thereof |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112613174A (en) | 2021-04-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Cho et al. | Isotherms of individual pores by gas adsorption crystallography | |
| CN110210460A (en) | A kind of shale gas apparent permeability calculation method for considering multiple factors and influencing | |
| Alnoaimi et al. | Experimental and numerical analysis of gas transport in shale including the role of sorption | |
| CN113670960B (en) | A real shale gas adsorption prediction method based on molecular simulation | |
| CN105866009B (en) | A kind of method and device calculating fine and close oily reservoir effecive porosity | |
| CN110672813B (en) | Shale gas content calculation method | |
| CN109490168B (en) | Quantitative evaluation method for the contribution rate of different mass transfer and diffusion mechanisms to reservoir seepage capacity in shale gas reservoirs | |
| CN105866002A (en) | Method for accurately measuring nuclear magnetic resonance porosity of oil-containing shale | |
| CN112613174B (en) | Evaluation method for methane adsorption capacity of shale considering multiple adsorption mechanisms | |
| Vandamme | Coupling between adsorption and mechanics (and vice versa) | |
| CN111896421B (en) | A method for calculating the real adsorption amount of methane in shale based on adsorption potential theory | |
| CN111537544B (en) | Improve nuclear magnetic resonance T 2 Conversion method for spectral characterization of pore size distribution precision of tight reservoir | |
| CN113075102A (en) | Method for establishing mathematical model of relation between spontaneous imbibition amount of porous medium and time | |
| CN102706778B (en) | Measuring device of effective diffusion coefficient of weak adsorption gas in bulk nano-porous material under low pressure | |
| CN107271314A (en) | A kind of method for measuring coal petrography adsorption swelling coefficient | |
| Fan et al. | A novel experimental system for accurate gas sorption and its application to various shale rocks | |
| Wang et al. | The characteristic development of micropores in deep coal and its relationship with adsorption capacity on the eastern margin of the Ordos Basin, China | |
| Zhao et al. | Experimental investigation on effect of water film thickness in unsaturated sandstone cores on CO2 transport during geologic storage | |
| Zhou et al. | Investigation of microscopic pore structure and permeability prediction in sand-conglomerate reservoirs | |
| Kolesnikov et al. | Pore size distribution of MCM-41-type silica materials from pseudomorphic transformation-A minimal input data approach based on excess surface work | |
| Zhang et al. | A review of molecular models for gas adsorption in shale nanopores and experimental characterization of shale properties | |
| CN110715879B (en) | Evaluation method of gas adsorption in micropores of highly evolved shale reservoirs based on gas-water distribution | |
| Hu et al. | Experimental and numerical study on scale effects of gas emission from coal particles | |
| CN101266237A (en) | A Disc Type Unsaturated Soil Moisture Diffusion Instrument | |
| Kreisberg et al. | Influence of the acid concentration on the morphology of micropores and mesopores in porous glasses |
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 |