CN114705591A - A method for calculating the wettability of coal/shale in high pressure environment - Google Patents
A method for calculating the wettability of coal/shale in high pressure environment Download PDFInfo
- Publication number
- CN114705591A CN114705591A CN202210329151.6A CN202210329151A CN114705591A CN 114705591 A CN114705591 A CN 114705591A CN 202210329151 A CN202210329151 A CN 202210329151A CN 114705591 A CN114705591 A CN 114705591A
- Authority
- CN
- China
- Prior art keywords
- coal
- shale
- gas
- adsorbed
- water
- 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.)
- Granted
Links
- 239000003245 coal Substances 0.000 title claims abstract description 240
- 238000000034 method Methods 0.000 title claims abstract description 33
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 119
- 238000009736 wetting Methods 0.000 claims abstract description 75
- 239000007788 liquid Substances 0.000 claims abstract description 41
- 238000001179 sorption measurement Methods 0.000 claims abstract description 19
- 229920006395 saturated elastomer Polymers 0.000 claims description 29
- 238000001228 spectrum Methods 0.000 claims description 24
- 238000012360 testing method Methods 0.000 claims description 15
- 239000000463 material Substances 0.000 claims description 14
- 238000002156 mixing Methods 0.000 claims description 10
- 238000005481 NMR spectroscopy Methods 0.000 claims description 8
- 239000007864 aqueous solution Substances 0.000 claims description 7
- 239000012153 distilled water Substances 0.000 claims description 5
- 238000001035 drying Methods 0.000 claims description 4
- 238000005076 Van der Waals potential Methods 0.000 claims description 3
- 239000000843 powder Substances 0.000 claims description 2
- 238000005259 measurement Methods 0.000 abstract description 20
- 230000003993 interaction Effects 0.000 abstract description 7
- 230000008569 process Effects 0.000 abstract description 5
- 239000007789 gas Substances 0.000 description 90
- 239000011148 porous material Substances 0.000 description 10
- 239000007787 solid Substances 0.000 description 9
- 238000004364 calculation method Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 6
- 238000000691 measurement method Methods 0.000 description 6
- 230000003595 spectral effect Effects 0.000 description 5
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 238000010997 low field NMR spectroscopy Methods 0.000 description 3
- 239000004696 Poly ether ether ketone Substances 0.000 description 2
- JUPQTSLXMOCDHR-UHFFFAOYSA-N benzene-1,4-diol;bis(4-fluorophenyl)methanone Chemical compound OC1=CC=C(O)C=C1.C1=CC(F)=CC=C1C(=O)C1=CC=C(F)C=C1 JUPQTSLXMOCDHR-UHFFFAOYSA-N 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000010191 image analysis Methods 0.000 description 2
- 230000008676 import Effects 0.000 description 2
- 238000011065 in-situ storage Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 229920002530 polyetherether ketone Polymers 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 230000009919 sequestration Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 238000000685 Carr-Purcell-Meiboom-Gill pulse sequence Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000000084 colloidal system Substances 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000005213 imbibition Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012900 molecular simulation Methods 0.000 description 1
- 238000000465 moulding Methods 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 239000002352 surface water Substances 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
- G01N13/02—Investigating surface tension of liquids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- 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
- G06F30/25—Design optimisation, verification or simulation using particle-based methods
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
- G01N13/02—Investigating surface tension of liquids
- G01N2013/0283—Investigating surface tension of liquids methods of calculating surface tension
-
- 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)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Analytical Chemistry (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Computational Mathematics (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Force Measurement Appropriate To Specific Purposes (AREA)
Abstract
本发明提供了一种高压环境中煤/页岩润湿性的计算方法,包括:步骤1、高压CO2环境中,煤/页岩表面吸附有吸附水和CO2吸附气体,煤/页岩表面具有混合表面的润湿性特性;步骤2、由液体在混合表面的润湿性的关系式,结合步骤1的条件,得到吸附有气体的煤/页岩表面的润湿性表达式;步骤3、由润湿角杨氏方程结合“sharp‑kink”近似定理,得到未被气体吸附的煤/页岩表面润湿角的表达式;步骤4、结合前述更新得到高压CO2环境中吸附有气体的煤/页岩表面润湿角的表达式。本发明解决了现有技术煤/页岩润湿性测量过程中,缺少考虑气体吸附与水‑气‑煤/页岩之间相互作用对润湿性改变因素的问题。
The present invention provides a method for calculating the wettability of coal/shale in a high-pressure environment, including: step 1. In a high-pressure CO 2 environment, adsorbed water and CO 2 adsorbed gas are adsorbed on the surface of the coal/shale, and the coal/shale is adsorbed on the surface. The surface has the wettability characteristics of the mixed surface; step 2, from the wettability relationship of the liquid on the mixed surface, combined with the conditions of step 1, obtain the wettability expression of the gas-adsorbed coal/shale surface; step 3. From the Young's equation of wetting angle combined with the "sharp-kink" approximation theorem, the expression of the wetting angle of the coal/shale surface that is not adsorbed by gas is obtained ; Expression for the wetting angle of the coal/shale surface of the gas. The invention solves the problem that the wettability changing factor of the interaction between gas adsorption and water-gas-coal/shale is not considered in the prior art coal/shale wettability measurement process.
Description
所属领域Field of study
本发明属于煤层气、页岩气勘探开发以及CO2地质封存过程中储层润湿性测量评价的技术领域,具体涉及一种高压环境中煤/页岩润湿性的计算方法。The invention belongs to the technical field of reservoir wettability measurement and evaluation in the process of coalbed methane and shale gas exploration and development and CO 2 geological storage, and in particular relates to a calculation method of coal/shale wettability in high pressure environment.
背景技术Background technique
润湿性是指一种液体在一种固体表面铺展的能力,高压气体环境中煤/页岩润湿性判断对于煤层气和页岩气开发、煤层或者页岩储层的CO2封存至关重要。高压气体环境中煤/页岩润湿性是影响天然气开发过程中地层中气、水分布,气、水两相渗流过程的重要因素。同时对于页岩层,润湿性的判断是评价其作为CO2封存盖层重要指标。在高压CO2环境中,若水湿性的页岩被改变成CO2气润湿,盖层毛管力减小,则失去了其盖层的封闭作用,导致CO2泄露。Wettability refers to the ability of a liquid to spread on a solid surface. Judgment of the wettability of coal/shale in high pressure gas environments is critical for CBM and shale gas development, and CO 2 sequestration in coal or shale reservoirs. important. The wettability of coal/shale in high-pressure gas environment is an important factor affecting the distribution of gas and water in the formation and the two-phase seepage process of gas and water in the process of natural gas development. At the same time, for shale layers, the judgment of wettability is an important index to evaluate it as a CO 2 sequestration caprock. In the high-pressure CO2 environment, if the water-wet shale is changed to be wetted by CO2 gas, the capillary force of the caprock is reduced, and the sealing effect of the caprock is lost, resulting in CO2 leakage.
常规环境中润湿性定量测定的方法主要有:(1)离心机法(USBM,United StatesBureau of Mines)、Amott方法(渗吸与排驱法),这两种方法均基于常压下气、水驱替实验,由于煤/页岩样品多为低孔低渗,驱替实验较难操作。(2)接触角测量法,方法测量结果直观准确,是煤或页岩润湿性测量的主要方法。The methods for quantitative determination of wettability in conventional environments mainly include: (1) Centrifuge method (USBM, United States Bureau of Mines), Amott method (imbibition and expulsion method), both of which are based on gas under normal pressure, In the water flooding experiment, because the coal/shale samples are mostly low porosity and low permeability, the flooding experiment is difficult to operate. (2) Contact angle measurement method, the measurement results are intuitive and accurate, and it is the main method for measuring the wettability of coal or shale.
高压CO2环境中,由于CO2气体在煤/页岩表面的吸附以及气-水界面张力的改变等原因导致煤/页岩表面润湿性改变。目前高压气体环境中煤/页岩的水润湿性的测量技术发展较为缓慢。主要的测量技术是将样品置于高压气体腔内,对水接触角进行测量。此类方法存在的最大问题是由于测量时间短,很难考虑气体吸附以及与水-气-煤/页岩之间相互作用对润湿性的改变。因此亟需建立一套高压气体环境中,考虑气体吸附,水-气-煤/页岩之间相互作用及气-水界面张力变化等多种因素的润湿性测量新方法。In the high pressure CO2 environment, the wettability of the coal/shale surface changes due to the adsorption of CO2 gas on the coal/shale surface and the change of the gas-water interfacial tension. At present, the measurement technology of water wettability of coal/shale in high pressure gas environment is relatively slow. The main measurement technique is to place the sample in a high pressure gas chamber and measure the water contact angle. The biggest problem with such methods is that due to the short measurement time, it is difficult to take into account the changes in wettability caused by gas adsorption and interactions with water-gas-coal/shale. Therefore, it is urgent to establish a new wettability measurement method in a high-pressure gas environment that takes into account gas adsorption, the interaction between water-gas-coal/shale, and the change of gas-water interfacial tension.
发明内容SUMMARY OF THE INVENTION
本发明提供了一种高压环境中煤/页岩润湿性的计算方法,解决了现有技术煤/页岩润湿性测量过程中,缺少考虑气体吸附以及与水-气-煤/页岩之间相互作用对润湿性改变因素的问题。The invention provides a method for calculating the wettability of coal/shale in a high pressure environment, which solves the problem of the lack of consideration of gas adsorption and water-gas-coal/shale wettability measurement process in the prior art. Interactions between wettability altering factors.
一种高压环境中煤/页岩润湿性的计算方法,方法包括:A method for calculating the wettability of coal/shale in a high pressure environment, comprising:
步骤1、高压CO2环境中,CO2气体吸附在部分煤/页岩表面置换一部分吸附在煤/页岩表面的吸附水,吸附在煤/页岩表面的CO2气体为CO2吸附气体,煤/页岩表面吸附有吸附水和所述CO2吸附气体,所述CO2吸附气体的外表面接触有水,煤/页岩表面具有混合表面的润湿性特性;
步骤2、由液体在混合表面的润湿性的关系式,结合步骤1的条件,得到高压CO2环境中,吸附有CO2气体的煤/页岩表面的润湿性表达式;
其中,液体在混合表面的润湿性的关系式为:Among them, the relationship between the wettability of the liquid on the mixed surface is:
式中,θ表示混合表面润湿角;γLG表示气液表面张力;γi,SG、γi,SL分别表示混合表面中每一种材质与气体和液体的界面张力;i表示混合表面中的材质;fi表示每种材质与混合表面的接触面积在混合表面中所占比例,f1+f2+…fn=1,n表示混合表面中材质的总数量;θi表示液体与每种材质之间的润湿角;In the formula, θ is the wetting angle of the mixed surface; γ LG is the gas-liquid surface tension; material; f i represents the proportion of the contact area of each material and the mixing surface in the mixing surface, f 1 +f 2 +…f n =1, n represents the total number of materials in the mixing surface; θ i represents the liquid and Wetting angle between each material;
高压CO2环境中,煤/页岩表面由未被气体吸附的煤/页岩表面和被气体吸附的煤/页岩表面组成,吸附有CO2气体的煤/页岩表面的润湿性表达式为:In the high pressure CO2 environment, the coal/shale surface is composed of the coal/shale surface not adsorbed by gas and the coal/shale surface adsorbed by gas, and the wettability expression of the coal/shale surface adsorbed by CO2 gas The formula is:
cosθ=f1cosθ1+f2cosθ2 (2)cosθ=f 1 cosθ 1 +f 2 cosθ 2 (2)
式中,θ表示高压CO2环境中煤/页岩表面润湿角,f1表示未被吸附气体的煤/页岩表面所占比例,即煤/页岩表面和所述吸附水的接触面积占煤/页岩表面的比例;θ1表示未被气体吸附的煤/页岩表面与所述水之间润湿角;f2表示被气体吸附的煤表面所占比例,即所述CO2吸附气体和水的接触面积占煤/页岩表面的比例;θ2表示被气体吸附的煤表面和水之间的润湿角,即所述CO2吸附气体与水之间润湿角为π。In the formula, θ represents the wetting angle of the coal/shale surface in the high-pressure CO2 environment, and f1 represents the proportion of the coal/shale surface without adsorbed gas, that is, the contact area between the coal/shale surface and the adsorbed water accounts for The ratio of coal/shale surface; θ1 represents the wetting angle between the coal/shale surface not adsorbed by gas and the water; f2 represents the proportion of coal surface adsorbed by gas, that is, the CO2 adsorbed gas The contact area with water accounts for the proportion of coal/shale surface; θ 2 represents the wetting angle between the coal surface adsorbed by the gas and water, that is, the wetting angle between the CO 2 adsorbed gas and water is π.
步骤3、由润湿角杨氏方程结合“sharp-kink”近似定理,得到未被气体吸附的煤/页岩表面与水之间润湿角的表达式:
式中,θ1表示未被气体吸附的煤/页岩表面与水之间润湿角,I表示范德华势能积分;Δρ表示气液密度差;γLG表示气液界面张力,ρg表示气体密度;ρlf表示煤/页岩表面液体膜密度等于液体密度。where θ 1 is the wetting angle between the coal/shale surface not adsorbed by gas and water, I is the integral of the van der Waals potential energy; Δρ is the gas-liquid density difference; γ LG is the gas-liquid interfacial tension, and ρ g is the gas density ; ρ lf indicates that the coal/shale surface liquid film density is equal to the liquid density.
步骤4、结合步骤2和3,更新高压CO2环境中,吸附有CO2气体的煤/页岩表面润湿角的表达式为:
式中θ是高压CO2环境中煤/页岩表面润湿角;对于煤/页岩表面和水溶液不变的恒温情况下,已知是常数,f1表示未被吸附气体的煤/页岩表面所占比例,即煤/页岩表面和所述吸附水的接触面积占煤/页岩表面的比例,ρg表示高压CO2环境中的气体密度;ρlf表示煤/页岩表面的液体膜密度等于水溶液密度。where θ is the wetting angle of the coal/shale surface in the high pressure CO environment; for the constant temperature of the coal/shale surface and the aqueous solution, it is known that is a constant, f 1 represents the proportion of coal/shale surface without adsorbed gas, that is, the contact area of coal/shale surface and the adsorbed water accounts for the proportion of coal/shale surface, ρ g represents high pressure CO environment gas density in ; ρ lf means that the density of the liquid film on the coal/shale surface is equal to the density of the aqueous solution.
进一步地,步骤4中所述更新的高压CO2环境中,吸附有CO2气体的煤/页岩润湿角表达式内的参数的测量方法为:Further, in the updated high-pressure CO 2 environment described in
步骤401、使用低场核磁共振技术测量不同CO2压力下饱和水煤样中的水的T2谱图;Step 401, using low-field nuclear magnetic resonance technology to measure the T2 spectrum of water in the water-saturated coal sample under different CO2 pressures;
步骤402、对所述T2谱图的数据分析,获得不同CO2压力下煤/页岩表面和水之间的接触面积,计算出不同CO2压力下煤/页岩表面和水之间的接触面积占煤/页岩表面面积的比例;Step 402 , analyze the data of the T2 spectrum, obtain the contact area between the coal/shale surface and water under different CO 2 pressures, and calculate the contact between the coal/shale surface and water under different CO 2 pressures Area as a percentage of coal/shale surface area;
步骤403、测量标况下煤/页岩表面的润湿角,导入所述更新高压CO2环境中,吸附有CO2气体的煤/页岩润湿角的表达式(5)内,得到值;Step 403: Measure the wetting angle of the coal/shale surface under standard conditions, and import it into the updated high-pressure CO 2 environment, in the expression (5) of the wetting angle of the coal/shale with CO 2 gas adsorbed, to obtain value;
进一步地,不同CO2压力下煤/页岩表面润湿角的计算方法为:将步骤402中所述不同CO2压力下煤/页岩表面和水之间的接触面积占煤/页岩表面面积的比例和步骤403中得到的所述值导入所述更新的高压CO2环境中,吸附有CO2气体的煤/页岩润湿角的表达式(5)内,得到不同CO2压力下煤/页岩表面润湿角。Further, the calculation method of the wetting angle of the coal/shale surface under different CO 2 pressures is as follows: the contact area between the coal/shale surface and the water under different CO 2 pressures described in step 402 is accounted for the coal/shale surface. The ratio of the area and the one obtained in step 403 as described The value is imported into the expression (5) of the wetting angle of coal/shale adsorbed with CO 2 gas in the updated high-pressure CO 2 environment, and the surface wetting angle of coal/shale under different CO 2 pressures is obtained.
进一步地,步骤402中所述不同CO2压力下煤/页岩表面和水之间的接触面积,以标况下饱和水煤样中的水的T2谱图数据为参照。Further, the contact area between the coal/shale surface and water under different CO 2 pressures described in step 402 is based on the T2 spectrum data of water in the water-saturated coal sample under standard conditions.
进一步地,步骤402中所述不同CO2压力下煤/页岩表面和水之间的接触面积占煤/页岩表面面积的比例的公式为:Further, the formula for the ratio of the contact area between the coal/shale surface and water to the coal/shale surface area under different CO 2 pressures described in step 402 is:
式中,f1表示煤/页岩表面和水之间接触面积占煤/页岩表面面积的比例,P1-0表示标况下T2谱图中左边第一个峰的谱峰信号幅度,P1-n表示不同压力下T2谱图中左边第一个峰的谱峰信号幅度。In the formula, f 1 represents the ratio of the contact area between the coal/shale surface and water to the coal/shale surface area, P 1-0 represents the spectral peak signal amplitude of the first peak on the left in the T2 spectrum under standard conditions, P 1-n represents the spectral peak signal amplitude of the first peak on the left in the T2 spectrum under different pressures.
进一步地,步骤403中所述标况下煤/页岩表面的润湿角的测量方法为:Further, the method for measuring the wetting angle of the coal/shale surface under the standard conditions described in step 403 is:
取煤/页岩样粉碎至大于200目,取煤粉制成具有光滑平面的试片;恒湿室内,在试片上滴落小于4mg的蒸馏水水滴;拍照,测量水滴与试片表面接触的长度和水滴半径,得到接触角;进而得到标况下煤/页岩表面润湿角。Take the coal/shale sample and crush it to a size larger than 200 mesh, and take the coal powder to make a test piece with a smooth plane; in the constant humidity chamber, drop the distilled water droplets of less than 4 mg on the test piece; take a picture, and measure the length of the contact between the water droplet and the test piece surface and the water droplet radius to obtain the contact angle; and then obtain the wetting angle of the coal/shale surface under standard conditions.
进一步地,所述步骤401中饱和水煤样的制作过程为:Further, the production process of the saturated water-coal sample in the step 401 is:
粉碎煤/页岩样制成样品,干燥至恒重后,取定量样品放入抽真空饱和水中,得到饱和水煤样。Pulverized coal/shale samples were made into samples, and after drying to constant weight, quantitative samples were taken into vacuum-saturated water to obtain saturated water-saturated coal samples.
本发明与现有技术相比,具有如下有益效果:Compared with the prior art, the present invention has the following beneficial effects:
1、本发明提供的计算方法适用于在高压CO2环境中,综合考虑了气体吸附作用及水-气-煤/页岩之间相互作用对润湿性的改变等因素,得到更优化的煤/页岩润湿性测量结果。1. The calculation method provided by the present invention is suitable for high pressure CO 2 environment, comprehensively considering factors such as gas adsorption and the change of wettability caused by the interaction between water-gas-coal/shale, etc., to obtain a more optimized coal. / Shale wettability measurements.
2、本发明实现了对不同气体压力下煤/页岩润湿性的计算和测量,是一种可即时的、原位的、动态的测量煤/页岩润湿性的方法。2. The present invention realizes the calculation and measurement of the wettability of coal/shale under different gas pressures, and is a method for measuring the wettability of coal/shale in real time, in situ and dynamically.
附图说明Description of drawings
为了更清楚地说明本发明的实施方式或现有技术中的技术方案,下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍。显而易见地,下面描述中的附图仅仅是示例性的,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图引伸获得其它的实施附图。In order to illustrate the embodiments of the present invention or the technical solutions in the prior art more clearly, the following briefly introduces the accompanying drawings that are required to be used in the description of the embodiments or the prior art. Obviously, the drawings in the following description are only exemplary, and for those of ordinary skill in the art, other implementation drawings can also be obtained according to the extension of the drawings provided without creative efforts.
图1是煤/页岩孔隙表面吸附点位及分子模拟水分子在煤/页岩表面吸附的示意图;Fig. 1 is a schematic diagram of adsorption sites on coal/shale pore surface and adsorption of molecular simulated water molecules on coal/shale surface;
图2是CO2环境中CO2吸附对水分子置换的示意图。Figure 2 is a schematic diagram of the displacement of water molecules by CO adsorption in a CO environment.
图3是煤表面在CO2环境中呈现固体和吸附气组成的混合表面。Figure 3 is a coal surface showing a mixed surface of solid and adsorbed gas composition in a CO2 environment.
图4是本发明实施例中0-6MPa CO2气体中样品T2谱图变化的示意图;4 is a schematic diagram of the change of the T spectrum of the sample in the 0-6MPa CO gas in the embodiment of the present invention ;
图5是本发明实施例中样品在标准状况下润湿角的测量图。FIG. 5 is a measurement diagram of the wetting angle of the sample under standard conditions in the embodiment of the present invention.
图中标号:Labels in the figure:
1为吸附点位;2为吸附水;3为煤/页岩表面;4为游离水;5为CO2。1 is adsorption site; 2 is adsorbed water; 3 is coal/shale surface; 4 is free water; 5 is CO 2 .
具体实施方式Detailed ways
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, but not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.
煤/页岩孔隙表面存在一系列吸附点位,这些吸附点位是煤表面水分附着的主要位点,水分吸附在吸附点位上,并借助水分子间氢键逐渐形成多层水分子簇。如图1所示,煤/页岩孔隙表面吸附点位示意图及分子模拟吸附水在煤/页岩表面吸附。当CO2等吸附性气体吸附在煤/页岩表面时,由于CO2与煤/页岩表面非吸附位点间相互作用力更强,因此CO2气体吸附在煤/页岩表面置换取代一些吸附水,如图2所示,被取代的吸附水运移到孔隙中间形成游离水,图2中空白处均为游离水,游离水覆盖了吸附在煤/页岩表面的CO2气体周围。而直接吸附在吸附点位上的吸附水由于其间作用力较强难以被置换,仍然会吸附在煤/页岩孔隙表面。因此,高压CO2环境中,CO2吸附到煤/页岩表面,置换一部分表面的吸附水,如图3所示。CO2对水分子的置换能力主要与煤表面对CO2的吸附能力有关。There are a series of adsorption sites on the pore surface of coal/shale. These adsorption sites are the main sites for moisture attachment on the coal surface. As shown in Figure 1, the schematic diagram of the adsorption sites on the coal/shale pore surface and the molecular simulation of the adsorption of water on the coal/shale surface. When CO2 is adsorbed on the surface of coal/shale, due to the stronger interaction between CO2 and the non-adsorbed sites on the surface of coal/shale, the CO2 gas adsorbed on the surface of coal/shale replaces some Adsorbed water, as shown in Fig. 2, the displaced adsorbed water migrates to the middle of the pores to form free water, the blank spaces in Fig. 2 are all free water, and the free water covers the CO gas adsorbed on the surface of coal/shale. However, the adsorbed water directly adsorbed on the adsorption site is difficult to be replaced due to the strong force between them, and will still be adsorbed on the surface of coal/shale pores. Therefore, in a high-pressure CO2 environment, CO2 is adsorbed to the coal/shale surface, displacing a part of the adsorbed water on the surface, as shown in Figure 3. The displacement capacity of CO2 for water molecules is mainly related to the adsorption capacity of CO2 on the coal surface.
本发明中综合气体吸附作用及水-气-煤/页岩之间相互作用对润湿性的改变因素,提供了一种高压CO2环境中煤/页岩润湿性的计算方法,In the present invention, a method for calculating the wettability of coal/shale in a high-pressure CO 2 environment is provided by synthesizing the change factors of gas adsorption and the interaction between water-gas-coal/shale on wettability,
步骤1、高压CO2环境中,CO2气体吸附在煤/页岩表面置换部分吸附水,吸附在煤/页岩表面的部分CO2气体为CO2吸附气体。在煤/页岩表面形成非均质混合表面,混合表面由未被气体吸附的煤/页岩表面和被气体吸附的煤/页岩表面组成,如图3所示。因此,当煤/页岩置于高压CO2环境中,煤/页岩表面水分分布在CO2吸附作用下发生改变,煤/页岩表面吸附有吸附水和所述CO2吸附气体,所述CO2吸附气体的外表面接触有游离水,煤/页岩表面具有混合表面的润湿性特性。
步骤2、由液体在混合表面的润湿性的关系式,结合步骤1的条件,得到高压CO2环境中,吸附有CO2气体的煤/页岩表面的润湿性表达式。
其中,液体在混合表面的润湿性的关系式为:Among them, the relationship between the wettability of the liquid on the mixed surface is:
式中,θ表示混合表面润湿角;L表示液体、G表示气体、S表示固体、γ表示表面张力;γLG,γSL,γSG分别表示气液、固液、气固表面张力;γi,SG、γi,SL分别表示混合表面中每一种材质与气体和液体的界面张力;i表示混合表面中的材质;fi表示每种材质与混合表面的接触面积在混合表面中所占比例,f1+f2+…fn=1,n表示混合表面中材质的总数量;θi表示液体与每种材质之间的润湿角。In the formula, θ represents the wetting angle of the mixed surface; L represents liquid, G represents gas, S represents solid, and γ represents surface tension; γ LG , γ SL , γ SG represent the surface tension of gas-liquid, solid-liquid, and gas-solid, respectively; γ i, SG , γ i, SL represent the interfacial tension between each material in the mixed surface and gas and liquid, respectively; i represents the material in the mixed surface; f i represents the contact area between each material and the mixed surface in the mixed surface. Proportion, f 1 +f 2 +...f n =1, n represents the total number of materials in the mixed surface; θ i represents the wetting angle between the liquid and each material.
结合步骤1中高压CO2环境中,煤/页岩表面吸附有吸附水和所述CO2吸附气体,吸附有CO2气体的煤/页岩表面的润湿性表达式:Combined with the high-pressure CO 2 environment in
cosθ=f1cosθ1+f2cosθ2 (2)cosθ=f 1 cosθ 1 +f 2 cosθ 2 (2)
式中,θ表示高压CO2环境中煤/页岩表面润湿角,f1表示未被吸附气体的煤/页岩表面所占比例,即煤/页岩表面和所述吸附水的接触面积占煤/页岩表面的比例;θ1表示未被气体吸附的煤/页岩表面与水之间润湿角;f2表示被气体吸附的煤表面所占比例,即所述CO2吸附气体和煤/页岩表面的接触面积所占比例,由于所述CO2吸附气体吸附在煤/页岩表面上,所述CO2吸附气体的外表面接触有游离水,所以所述CO2吸附气体和煤/页岩表面的接触面积等于所述CO2吸附气体和水的接触面积;θ2表示被气体吸附的煤表面和水之间的润湿角,即所述CO2吸附气体与水之间润湿角,为π。In the formula, θ represents the wetting angle of the coal/shale surface in the high-pressure CO2 environment, and f1 represents the proportion of the coal/shale surface without adsorbed gas, that is, the contact area between the coal/shale surface and the adsorbed water accounts for The ratio of coal/shale surface; θ1 represents the wetting angle between the coal/shale surface not adsorbed by gas and water; f2 represents the proportion of coal surface adsorbed by gas, that is, the CO2 adsorbed gas and coal The ratio of the contact area to the surface of the shale/shale. Since the CO 2 adsorbed gas is adsorbed on the coal/shale surface, and the outer surface of the CO 2 adsorbed gas is in contact with free water, the CO 2 adsorbed gas and coal / The contact area of the shale surface is equal to the contact area of the CO 2 adsorbed gas and water; θ 2 represents the wetting angle between the coal surface adsorbed by the gas and water, that is, the wetting angle between the CO 2 adsorbed gas and water Wet angle, π.
由于f1+f2=1,则高压CO2环境中煤/页岩表面润湿角的计算公式为:Since f 1 +f 2 =1, the calculation formula of the wetting angle of coal/shale surface in high pressure CO 2 environment is:
cosθ=f1cosθ1-f2=f1(1+cosθ1)-1 (3)cosθ=f 1 cosθ 1 -f 2 =f 1 (1+cosθ 1 )-1 (3)
步骤3、在不同的气体压力下,受到气体浓度影响,气液界面张力发生改变,未被气体吸附的煤/页岩表面与水之间润湿角θ1也随着压力变化发生变化,由润湿角杨氏方程(式中γGS表示气固表面张力;γLS表示液固表面张力;γLG表示气液界面张力)结合“sharp-kink”近似定理(由文献《Roshan,H.,A.Z.Al-Yaseri,M.Sarmadivalehand S.Iglauer(2016),On wettability of shale rocks.Journal of Colloid andInterface Science,475,104-111.》可知),得到未考虑吸附气体作用下煤/页岩表面与水之间润湿角θ1的表达式,也就是本发明中未被气体吸附的煤/页岩表面与水之间润湿角θ1的表达式:
式中,I表示范德华势能积分;Δρ表示气液密度差;γLG表示气液界面张力,ρg表示气体密度;ρlf表示固体表面液体膜密度,一般等于液体密度。In the formula, I is the integral of van der Waals potential energy; Δρ is the gas-liquid density difference; γ LG is the gas-liquid interfacial tension, ρ g is the gas density; ρ lf is the liquid film density on the solid surface, which is generally equal to the liquid density.
步骤4、对于温度恒定,确定的煤/页岩表面和水溶液,结合步骤2和3,更新高压CO2环境中,吸附有CO2气体的煤/页岩表面润湿角的表达式为:
式中θ是高压CO2环境中煤/页岩表面润湿角;对于煤/页岩表面和水溶液不变的恒温情况下,已知是常数(由文献《Al-Yaseri,A.Z.,H.Roshan,M.Lebedev,A.Barifcani,and S.Iglauer(2016),Dependence of quartz wettability on fluiddensity.Geophysical Research Letters,43(8),3771-3776.》可知),f1表示未被吸附气体的煤/页岩表面所占比例,即煤/页岩表面和所述吸附水的接触面积也就是煤/页岩表面和水的接触面积占煤/页岩表面的比例,ρg表示高压CO2环境中的气体密度;ρlf表示煤/页岩表面的液体膜密度,一般等于水溶液密度。where θ is the wetting angle of the coal/shale surface in the high pressure CO environment; for the constant temperature of the coal/shale surface and the aqueous solution, it is known that is a constant (from the literature "Al-Yaseri, AZ, H. Roshan, M. Lebedev, A. Barifcani, and S. Iglauer (2016), Dependence of quartz wettability on fluiddensity. Geophysical Research Letters, 43(8), 3771- 3776.”), f 1 represents the proportion of the coal/shale surface without adsorbed gas, that is, the contact area between the coal/shale surface and the adsorbed water, that is, the contact area between the coal/shale surface and water. The ratio of coal/shale surface, ρ g is the gas density in high pressure CO2 environment; ρ lf is the liquid film density on the coal/shale surface, which is generally equal to the density of the aqueous solution.
因此,根据公式5,获得f1,ρg,ρlf和标况下(即常压下)的煤/页岩表面润湿角即可获取不同气体压力下煤/页岩表面润湿角。Therefore, according to
为了准确获取高压CO2环境中煤/页岩表面润湿角,根据以上高压润湿角计算原理,本发明提供一种高压环境中煤/页岩润湿性的参数测量方法,利用低场核磁共振谱图获取不同压力下f1的值,同时获取常压下煤/页岩表面润湿角,利用以上计算公式即可获得不同气体压力下煤/页岩润湿性的定量测定,是一种可即时的、原位的、动态的测量煤/页岩的润湿性的新方法。In order to accurately obtain the wetting angle of coal/shale surface in high pressure CO 2 environment, according to the above calculation principle of high pressure wetting angle, the present invention provides a method for measuring the wettability of coal/shale in high pressure environment. The value of f 1 under different pressures is obtained from the resonance spectrum, and the wetting angle of the coal/shale surface under normal pressure is obtained at the same time. Using the above calculation formula, the quantitative determination of the wettability of coal/shale under different gas pressures can be obtained. A new method for instant, in situ, and dynamic measurement of coal/shale wettability.
针对上述公式5中的未知数,本发明提供一种测量方法,包括:For the unknowns in the
步骤401、使用低场核磁共振技术测量不同CO2压力下饱和水煤样中的水的T2谱图;Step 401, using low-field nuclear magnetic resonance technology to measure the T2 spectrum of water in the water-saturated coal sample under different CO2 pressures;
步骤402、对所述T2谱图数据分析,获得不同CO2压力下煤/页岩表面和水之间的接触面积,计算出不同CO2压力下煤/页岩表面和水之间的接触面积占煤/页岩表面面积的比例;Step 402: Analyze the T2 spectrum data to obtain the contact area between the coal/shale surface and water under different CO 2 pressures, and calculate the contact area between the coal/shale surface and water under different CO 2 pressures percentage of coal/shale surface area;
步骤403、测量标况下煤/页岩表面的润湿角,导入所述煤/页岩润湿性公式5内,得到值;Step 403: Measure the wetting angle of the coal/shale surface under standard conditions, import it into the coal/
步骤404、将所述不同CO2压力下煤/页岩表面和液体之间的接触面积占煤/页岩表面面积的比例f1和值导入所述煤/页岩润湿性公式5内,得到不同压力下煤/页岩表面润湿角。Step 404: Calculate the ratio f1 of the contact area between the coal/shale surface and the liquid to the coal/shale surface area under different CO The value is imported into the coal/
本发明提供的测量方法中,将低场核磁共振技术与煤/页岩表面水接触润湿角测量技术结合起来,通过建立煤/页岩的润湿性与核磁共振T2谱图信号的标度关系,实现了对不同气体压力下煤/页岩润湿性的定量测定。In the measurement method provided by the present invention, the low-field nuclear magnetic resonance technology is combined with the coal/shale surface water contact wetting angle measurement technology, and the scale of the wettability of the coal/shale and the nuclear magnetic resonance T2 spectrum signal is established. A quantitative determination of the wettability of coal/shale under different gas pressures is achieved.
在本申请中,不同CO2压力下对饱和水煤样中的水进行低场核磁共振测量,获得煤中水的T2谱图,根据T2谱图左边第一个峰(P1)的幅度代表煤中孔隙表面吸附水。在不同的气体压力下,P1峰面积发生变化,代表煤中CO2气体吸附到煤孔隙表面,置换吸附水,这部分水不再与孔表面接触,变成游离水运移到大孔中间。因此利用NMR技术(核磁共振)对孔表面剩余吸附水的定量测量,即可获得不同气体压力下固体和液体之间接触面积即f1。根据公式结合T2谱图得到的f1值及不同压力下CO2气体密度,和标况下的煤/页岩表面润湿角,从而求得设定气体压力下煤/页岩的水润湿角。In this application, low-field NMR measurements were performed on water in water-saturated coal samples under different CO 2 pressures, and the T 2 spectrum of water in coal was obtained. According to the amplitude of the first peak (P1) on the left of the T 2 spectrum Represents the adsorbed water on the pore surface of coal. Under different gas pressures, the area of the P1 peak changes, which means that CO2 gas in the coal is adsorbed to the coal pore surface and replaces the adsorbed water. This part of the water is no longer in contact with the pore surface, and becomes free water and migrates to the middle of the large pores. Therefore, the quantitative measurement of the remaining adsorbed water on the pore surface by NMR technology (nuclear magnetic resonance) can obtain the contact area between the solid and the liquid under different gas pressures, that is, f 1 . According to the formula Combined with the f 1 value obtained from the T2 spectrum, the CO 2 gas density under different pressures, and the surface wetting angle of coal/shale under standard conditions, the water-wetting angle of coal/shale under the set gas pressure was obtained.
在一个可能的实施例中,测量方法的步骤如下:In a possible embodiment, the steps of the measurement method are as follows:
1、制样及装样1. Sample preparation and loading
将煤/页岩样粉碎取直径约1cm左右的煤样颗粒若干,放入干燥箱中干燥至恒重,称重,制成样品;将样品泡入抽真空饱和水中至48小时,保证样品完全饱和水,即得到饱和水煤样。取完全饱和水样品,擦去完全饱和水样品表面的多余水分,将完全饱和水样品装入PEEK耐高压样品腔中。Pulverize the coal/shale sample and take a number of coal sample particles with a diameter of about 1 cm, put it in a drying box to dry to constant weight, and weigh to make a sample; soak the sample in vacuum-saturated water for 48 hours to ensure that the sample is completely Saturated water, that is, saturated water-coal samples. Take the fully saturated water sample, wipe off the excess water on the surface of the fully saturated water sample, and put the fully saturated water sample into the PEEK high pressure resistant sample chamber.
2、选取低场核磁共振测量参数并设定环境温度2. Select the low-field NMR measurement parameters and set the ambient temperature
选用CPMG(Carr-Purcell-Meiboom-Gill)测量序列;测量具体参数包括:回波时间(TE)、叠加次数(SCANS)、等待时间(TW)、回波数(NECH);设定测试环境温度至适当温度并保持恒温。Select the CPMG (Carr-Purcell-Meiboom-Gill) measurement sequence; the specific measurement parameters include: echo time (TE), stacking times (SCANS), waiting time (TW), echo number (NECH); set the test ambient temperature to appropriate temperature and maintain a constant temperature.
所述核磁共振仪采用CPMG序列,测量具体参数优选如下:回波时间(TE)=0.3ms,叠加次数(SCANS)=64,等待时间(TW)=3s,回波数(NECH)=18000;所述适当环境温度为25℃。The NMR instrument adopts the CPMG sequence, and the specific measurement parameters are preferably as follows: echo time (TE)=0.3ms, stacking times (SCANS)=64, waiting time (TW)=3s, echo number (NECH)=18000; The appropriate ambient temperature is 25°C.
3、获取不同压力下f1值3. Obtain f 1 values under different pressures
测量标况下完全饱和水样品的T2谱图,记做T2-0,其中P1谱峰信号幅度为P1-0。将nMPa压力的CO2气体注入到样品腔中,静置10个小时,待CO2充分吸附到煤中,测试此时完全饱和水样品的T2谱图,记做T2-n,其中P1谱峰信号幅度为P1-n。因此,n MPa压力下的f1值为: Measure the T2 spectrum of the fully saturated water sample under standard conditions, denoted as T 2-0 , and the signal amplitude of the P1 spectrum peak is P 1-0 . Inject the CO 2 gas of nMPa pressure into the sample chamber, let it stand for 10 hours, and wait until the CO 2 is fully adsorbed into the coal. Test the T2 spectrum of the fully saturated water sample at this time, denoted as T 2-n , in which the P1 spectrum The peak signal amplitude is P 1-n . Therefore, the value of f1 at n MPa pressure is:
4、测量标准状况下煤/页岩表面润湿角。4. Measure the wetting angle of coal/shale surface under standard conditions.
取煤样/页岩粉碎至大于200目,取煤粉用加压成型模具在50MPa压力下压制成具有光滑平面的试片,可制成圆柱体或椭圆体试片。在恒湿室内,用注射器往试片上滴落不足4mg的蒸馏水水滴,用快速相机进行拍照,通过图像分析,测量水滴与表面接触的长度a,通过单圆拟合水滴形状获取圆半径r,由公式求出接触角θ.即标况下煤/页岩表面润湿角。利用公式4及已知的标况下煤/页岩表面润湿角和气液密度,此处液体为蒸馏水,固体表面液体膜密度ρlf为水密度取1g/cm3,求得值。Take the coal sample/shale and pulverize it to more than 200 mesh, take the pulverized coal and press it into a test piece with a smooth plane under a pressure of 50MPa, which can be made into a cylinder or ellipsoid test piece. In the constant humidity chamber, drop a drop of distilled water less than 4 mg onto the test piece with a syringe, take a picture with a fast camera, measure the length a of the contact between the water drop and the surface through image analysis, and obtain the circle radius r by fitting the shape of the water drop with a single circle. formula Obtain the contact angle θ. That is, the wetting angle of the coal/shale surface under standard conditions. Using
5、利用公式5,计算获取不同压力下的煤/页岩表面润湿角。5. Using
利用获得的不同压力下的f1的值和不同压力下CO2气体密度及步骤4中的值,即可获得不同压力下煤/页岩样品的表面润湿角。Using the obtained values of f1 at different pressures and CO2 gas densities at different pressures and the values in step 4 value, the surface wetting angles of coal/shale samples under different pressures can be obtained.
实施例1Example 1
为演示,本发明提供一个具体的应用实施例,操作步骤如下:For demonstration, the present invention provides a specific application embodiment, and the operation steps are as follows:
1、制样及装样1. Sample preparation and loading
选取鄂尔多斯盆地中煤阶煤样,将煤样粉碎取直径约1cm左右的煤样颗粒,放入干燥箱中干燥至恒重,称重19.88g;将样品进行抽真空饱和水48小时,保证样品完全饱和水。取完全饱和水样品,擦去完全饱和水样品表面多余水分,将完全饱和水样品装入PEEK耐高压样品腔中。Select the middle-rank coal sample in the Ordos Basin, pulverize the coal sample to obtain coal sample particles with a diameter of about 1 cm, put it in a drying oven to dry to a constant weight, and weigh 19.88 g; vacuum the sample with saturated water for 48 hours to ensure that the sample is Fully saturated water. Take the fully saturated water sample, wipe off the excess water on the surface of the fully saturated water sample, and put the fully saturated water sample into the PEEK high pressure resistant sample chamber.
2、选取低场核磁共振测量参数并设定环境温度2. Select the low-field NMR measurement parameters and set the ambient temperature
选用CPMG(Carr-Purcell-Meiboom-Gill)测量序列,测量具体参数如下:回波时间(TE)=0.3ms,叠加次数(SCANS)=64,等待时间(TW)=3s,回波数(NECH)=18000;设定适当环境温度为25℃。The CPMG (Carr-Purcell-Meiboom-Gill) measurement sequence is selected, and the specific parameters of the measurement are as follows: echo time (TE)=0.3ms, stacking times (SCANS)=64, waiting time (TW)=3s, echo number (NECH) =18000; set the appropriate ambient temperature to 25℃.
3、获取不同压力下f1值3. Obtain f 1 values under different pressures
如图4所示,图中纵坐标是“T2谱信号幅度(a.u.),测量标况下完全饱和水样品的T2谱图,记做T2-0,其中P1谱峰信号幅度为P1-0。向样品腔中注入1MPa压力的CO2气体,加压后静置10个小时,然后测试完全饱和水样品的T2谱图,记做T2-1,其中P1谱峰信号幅度为P1-1,此时,1MPa压力下的f1值为:继续向煤样中注入CO2气体,压力达到2MPa,静置10个小时,然后测试完全饱和水样品的T2谱图,记做T2-2,其中P1谱峰信号幅度为P1-2,2MPa压力下的f1值为:依次类推直至压力增大到6MPa,6MPa压力下的f1值为:依次计算得到f1数值如下:As shown in Figure 4, the ordinate in the figure is "T2 spectral signal amplitude (au), and the T2 spectral diagram of the fully saturated water sample under standard conditions is measured, marked as T 2-0 , and the P1 spectral peak signal amplitude is P 1- 0. Inject CO 2 gas with a pressure of 1MPa into the sample cavity, let it stand for 10 hours after pressurization, and then test the T2 spectrum of the fully saturated water sample, denoted as T 2-1 , where the P1 spectrum peak signal amplitude is P 1 -1 , at this time, the value of f 1 under 1MPa pressure is: Continue to inject CO 2 gas into the coal sample, the pressure reaches 2MPa, let it stand for 10 hours, and then test the T2 spectrum of the fully saturated water sample, denoted as T 2-2 , where the P1 spectrum peak signal amplitude is P 1-2 , The value of f1 under 2MPa pressure is: And so on until the pressure increases to 6MPa, the value of f 1 under 6MPa pressure is: The f 1 values are calculated sequentially as follows:
表1不同压力下P1峰幅度值及f1 Table 1 P1 peak amplitude and f 1 under different pressures
4、测量标准状况下煤/页岩表面润湿角。4. Measure the wetting angle of coal/shale surface under standard conditions.
取煤样粉碎至大于200目,取煤粉用加压成型模具在50MPa压力下压制成具有光滑平面的圆柱体试片。在恒湿室内,用注射器往试片上滴落不足4mg的蒸馏水水滴,用快速相机进行拍照,得到图5,通过图像分析,测量水滴与表面接触的长度a=1.94mm和半径r=0.99mm,由公式求出表面接触角78.46°.利用公式5及已知的标况下煤/页岩表面润湿角和气液密度求得值为1.2079。The coal sample was pulverized to a size larger than 200 mesh, and the pulverized coal was pressed into a cylindrical test piece with a smooth plane under a pressure of 50 MPa with a pressurized molding die. In the constant humidity chamber, drop the distilled water droplet of less than 4 mg on the test piece with a syringe, and take a picture with a fast camera to obtain Figure 5. Through image analysis, the length a=1.94mm and the radius r=0.99mm of the contact between the water droplet and the surface are measured. by formula Obtain the surface contact angle of 78.46°. Use
5、利用公式5,计算获取不同压力下煤/页岩样品表面润湿角。5. Using
利用获得的不同压力下的f1的值和不同压力下CO2气体密度及上述步骤中的值,即可获得不同压力下煤/页岩样品表面润湿角,如表2所示。Using the obtained values of f1 at different pressures and CO2 gas density at different pressures and the above steps value, the surface wetting angles of coal/shale samples under different pressures can be obtained, as shown in Table 2.
表2不同压力下煤/页岩样品表面润湿角Table 2 Surface wetting angles of coal/shale samples under different pressures
以上实施例仅为本申请的示例性实施例,不用于限制本申请,本申请的保护范围由权利要求书限定。本领域技术人员可以在本申请的实质和保护范围内,对本申请做出各种修改或等同替换,这种修改或等同替换也应视为落在本申请的保护范围内。The above embodiments are only exemplary embodiments of the present application, and are not intended to limit the present application. The protection scope of the present application is defined by the claims. Those skilled in the art can make various modifications or equivalent replacements to the present application within the spirit and protection scope of the present application, and such modifications or equivalent replacements should also be regarded as falling within the protection scope of the present application.
Claims (7)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210329151.6A CN114705591B (en) | 2022-03-30 | 2022-03-30 | Method for calculating wettability of coal/shale in high-pressure environment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210329151.6A CN114705591B (en) | 2022-03-30 | 2022-03-30 | Method for calculating wettability of coal/shale in high-pressure environment |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114705591A true CN114705591A (en) | 2022-07-05 |
CN114705591B CN114705591B (en) | 2024-08-02 |
Family
ID=82169824
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210329151.6A Active CN114705591B (en) | 2022-03-30 | 2022-03-30 | Method for calculating wettability of coal/shale in high-pressure environment |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114705591B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114577673A (en) * | 2022-01-18 | 2022-06-03 | 重庆大学 | A NMR characterization method for pore wettability of coal powder with equal particle size |
CN117423404A (en) * | 2023-10-18 | 2024-01-19 | 中国地质大学(北京) | Acquisition method, terminal and storage medium for wettability parameters of coal reservoir |
CN117804999A (en) * | 2024-02-29 | 2024-04-02 | 陕西延长石油(集团)有限责任公司 | Method for quantitatively characterizing influence of chemical action of carbon dioxide on wettability of sandstone |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2490614C1 (en) * | 2012-03-05 | 2013-08-20 | Учреждение Российской академии наук Институт проблем нефти и газа РАН | Method to determine differentiated wettability of minerals within reservoir formations |
CN108593697A (en) * | 2018-05-09 | 2018-09-28 | 中国地质大学(北京) | A kind of method and apparatus measuring Coal Under High Pressure wetability using low-field nuclear magnetic resonance |
CN109142154A (en) * | 2017-06-16 | 2019-01-04 | 中国石油化工股份有限公司 | A method of calculating shale gas reservoir angle of wetting |
-
2022
- 2022-03-30 CN CN202210329151.6A patent/CN114705591B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2490614C1 (en) * | 2012-03-05 | 2013-08-20 | Учреждение Российской академии наук Институт проблем нефти и газа РАН | Method to determine differentiated wettability of minerals within reservoir formations |
CN109142154A (en) * | 2017-06-16 | 2019-01-04 | 中国石油化工股份有限公司 | A method of calculating shale gas reservoir angle of wetting |
CN108593697A (en) * | 2018-05-09 | 2018-09-28 | 中国地质大学(北京) | A kind of method and apparatus measuring Coal Under High Pressure wetability using low-field nuclear magnetic resonance |
Non-Patent Citations (1)
Title |
---|
陈基瑜 等: "二氧化碳和氦气对低煤阶煤煤润湿性的影响研究", 《煤炭科学技术》, vol. 43, no. 11, 30 November 2015 (2015-11-30), pages 129 - 134 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114577673A (en) * | 2022-01-18 | 2022-06-03 | 重庆大学 | A NMR characterization method for pore wettability of coal powder with equal particle size |
CN114577673B (en) * | 2022-01-18 | 2024-02-02 | 重庆大学 | Pore wettability NMR characterization method of equal-particle-size pulverized coal |
CN117423404A (en) * | 2023-10-18 | 2024-01-19 | 中国地质大学(北京) | Acquisition method, terminal and storage medium for wettability parameters of coal reservoir |
CN117423404B (en) * | 2023-10-18 | 2024-04-26 | 中国地质大学(北京) | Method, terminal and storage medium for obtaining wettability parameters of coal reservoir |
CN117804999A (en) * | 2024-02-29 | 2024-04-02 | 陕西延长石油(集团)有限责任公司 | Method for quantitatively characterizing influence of chemical action of carbon dioxide on wettability of sandstone |
CN117804999B (en) * | 2024-02-29 | 2024-05-03 | 陕西延长石油(集团)有限责任公司 | Method for quantitatively characterizing influence of chemical action of carbon dioxide on wettability of sandstone |
Also Published As
Publication number | Publication date |
---|---|
CN114705591B (en) | 2024-08-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114705591A (en) | A method for calculating the wettability of coal/shale in high pressure environment | |
CN112378943B (en) | A shale oil saturation evaluation model, evaluation method and application | |
CN113075102B (en) | Method for establishing mathematical model of relation between spontaneous imbibition amount of porous medium and time | |
Song et al. | Multi-fractal distribution analysis for pore structure characterization of tight sandstone—A case study of the Upper Paleozoic tight formations in the Longdong District, Ordos Basin | |
US20240027379A1 (en) | Method for quantitative evaluation on sensitivity of shale oil and gas reservoir to injected fluids | |
CN107063919B (en) | The device and method of carbon dioxide and alkane competitive Adsorption amount in a kind of measurement shale | |
EP3286555B1 (en) | Nuclear magnetic resonance gas isotherm technique to evaluate reservoir rock wettability | |
CN103994960B (en) | A Coal/Shale Isothermal Adsorption Experimental Method | |
CN106153662A (en) | The measuring method of rock core stress sensitivity | |
CN103344541A (en) | Method for measuring total porosity of shale | |
CN110108616B (en) | Method for compensating signal loss of coal sample in centrifugal process in nuclear magnetic resonance experiment | |
Zhang et al. | Microscopic distribution and dynamic variation of water under stress in middle and high rank coal samples | |
CN105675468A (en) | Evaluation method of shale organic matter porosity | |
CN113607621A (en) | A method for characterizing closed pores and connectivity of coal-measure composite reservoirs | |
Plug et al. | Capillary pressure and wettability behavior of CO2 sequestration in coal at elevated pressures | |
CN108593697A (en) | A kind of method and apparatus measuring Coal Under High Pressure wetability using low-field nuclear magnetic resonance | |
Ren et al. | Influence of micro‐pore structure in tight sandstone reservoir on the seepage and water‐drive producing mechanism—a case study from Chang 6 reservoir in Huaqing area of Ordos basin | |
WO2020097037A1 (en) | Nuclear magnetic resonance gas isotherm technique to evaluate reservoir rock wettability | |
Wang et al. | Pore structure and fractal characteristics of Wufeng–Longmaxi formation shale in northern Yunnan–Guizhou, China | |
Xin et al. | Declining high-pressure sorption isotherms on shale and coal: systematic comparison of the contributing factors | |
Chen et al. | Pore and gas desorption characteristics of primary coal with different degrees of metamorphism | |
CN105547959B (en) | Analysis method for natural gas storage capacity of matrix pores in dense bedrock | |
CN115876659A (en) | A method and system for obtaining shale pore and fracture volume | |
Wang et al. | Study of fluid displacement processes and sequestration of CO2 in coal reservoirs using nuclear magnetic resonance spectroscopy | |
Liu et al. | Characterization of Micro-nano pore structure of tight sandstone based on nuclear magnetic resonance experiments (NMR) |
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 |