[go: up one dir, main page]

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 PDF

Info

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
Application number
CN202210329151.6A
Other languages
Chinese (zh)
Other versions
CN114705591B (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.)
China University of Geosciences Beijing
Original Assignee
China University of Geosciences Beijing
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 China University of Geosciences Beijing filed Critical China University of Geosciences Beijing
Priority to CN202210329151.6A priority Critical patent/CN114705591B/en
Publication of CN114705591A publication Critical patent/CN114705591A/en
Application granted granted Critical
Publication of CN114705591B publication Critical patent/CN114705591B/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
    • G01N13/00Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N13/00Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
    • G01N13/02Investigating surface tension of liquids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/25Design optimisation, verification or simulation using particle-based methods
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N13/00Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
    • G01N13/02Investigating surface tension of liquids
    • G01N2013/0283Investigating surface tension of liquids methods of calculating surface tension
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force 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环境中吸附有气体的煤/页岩表面润湿角的表达式。本发明解决了现有技术煤/页岩润湿性测量过程中,缺少考虑气体吸附与水‑气‑煤/页岩之间相互作用对润湿性改变因素的问题。

Figure 202210329151

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.

Figure 202210329151

Description

一种高压环境中煤/页岩润湿性的计算方法A method for calculating the wettability of coal/shale in high pressure environment

所属领域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吸附气体的外表面接触有水,煤/页岩表面具有混合表面的润湿性特性;Step 1. In a high-pressure CO2 environment, CO2 gas is adsorbed on part of the coal/shale surface to replace part of the adsorbed water adsorbed on the coal/shale surface, and the CO2 gas adsorbed on the coal/shale surface is CO2 adsorbed gas, Adsorbed water and the CO 2 adsorbed gas are adsorbed on the surface of the coal/shale, the outer surface of the CO 2 adsorbed gas is in contact with water, and the surface of the coal/shale has the wettability characteristics of a mixed surface;

步骤2、由液体在混合表面的润湿性的关系式,结合步骤1的条件,得到高压CO2环境中,吸附有CO2气体的煤/页岩表面的润湿性表达式;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 surface of coal/shale with CO 2 gas adsorbed in the high-pressure CO 2 environment;

其中,液体在混合表面的润湿性的关系式为:Among them, the relationship between the wettability of the liquid on the mixed surface is:

Figure BDA0003572566430000021
Figure BDA0003572566430000021

式中,θ表示混合表面润湿角;γ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”近似定理,得到未被气体吸附的煤/页岩表面与水之间润湿角的表达式:Step 3. From the Young's equation of the wetting angle combined with the "sharp-kink" approximation theorem, the expression of the wetting angle between the surface of coal/shale that is not adsorbed by gas and water is obtained:

Figure BDA0003572566430000031
Figure BDA0003572566430000031

式中,θ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气体的煤/页岩表面润湿角的表达式为:Step 4. Combining steps 2 and 3, update the expression of the wetting angle of the coal/shale surface with CO gas adsorbed in the high pressure CO 2 environment as:

Figure BDA0003572566430000032
Figure BDA0003572566430000032

式中θ是高压CO2环境中煤/页岩表面润湿角;对于煤/页岩表面和水溶液不变的恒温情况下,已知

Figure BDA0003572566430000033
是常数,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
Figure BDA0003572566430000033
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 step 4, the method for measuring the parameters in the expression of the wetting angle of coal/shale with CO 2 gas adsorbed is:

步骤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)内,得到

Figure BDA0003572566430000034
值;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
Figure BDA0003572566430000034
value;

进一步地,不同CO2压力下煤/页岩表面润湿角的计算方法为:将步骤402中所述不同CO2压力下煤/页岩表面和水之间的接触面积占煤/页岩表面面积的比例和步骤403中得到的所述

Figure BDA0003572566430000041
值导入所述更新的高压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
Figure BDA0003572566430000041
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:

Figure BDA0003572566430000042
Figure BDA0003572566430000042

式中,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为CO21 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吸附气体的外表面接触有游离水,煤/页岩表面具有混合表面的润湿性特性。Step 1. In a high-pressure CO 2 environment, CO 2 gas is adsorbed on the surface of coal/shale to replace part of the adsorbed water, and part of the CO 2 gas adsorbed on the surface of coal/shale is CO 2 adsorbed gas. A heterogeneous mixed surface is formed on the coal/shale surface, and the mixed surface consists of the coal/shale surface that is not adsorbed by gas and the coal/shale surface that is adsorbed by gas, as shown in Figure 3. Therefore, when coal/shale is placed in a high-pressure CO2 environment, the surface moisture distribution of coal/shale changes under the action of CO2 adsorption, and the surface of coal/shale is adsorbed with adsorbed water and the CO2 adsorbed gas, the The outer surface of the CO2 adsorbed gas is in contact with free water, and the coal/shale surface has the wettability characteristics of a mixed surface.

步骤2、由液体在混合表面的润湿性的关系式,结合步骤1的条件,得到高压CO2环境中,吸附有CO2气体的煤/页岩表面的润湿性表达式。Step 2. From the relational expression of the wettability of the liquid on the mixed surface, combined with the conditions of step 1, the wettability expression of the surface of the coal/shale with CO 2 gas adsorbed in the high pressure CO 2 environment is obtained.

其中,液体在混合表面的润湿性的关系式为:Among them, the relationship between the wettability of the liquid on the mixed surface is:

Figure BDA0003572566430000061
Figure BDA0003572566430000061

式中,θ表示混合表面润湿角;L表示液体、G表示气体、S表示固体、γ表示表面张力;γLGSLSG分别表示气液、固液、气固表面张力;γ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 step 1, the surface of coal/shale has adsorbed water and the CO 2 adsorbed gas, and the wettability expression of the surface of coal/shale adsorbed with CO 2 gas is:

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也随着压力变化发生变化,由润湿角杨氏方程(

Figure BDA0003572566430000071
式中γ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的表达式:Step 3. Under different gas pressures, under the influence of gas concentration, the gas-liquid interfacial tension changes, and the wetting angle θ1 between the surface of coal/shale that is not adsorbed by gas and water also changes with the change of pressure, as shown by Wetting angle Young's equation (
Figure BDA0003572566430000071
In the formula, γ GS represents the gas-solid surface tension; γ LS represents the liquid-solid surface tension; γ LG represents the gas-liquid interfacial tension) combined with the "sharp-kink" approximation theorem (by the literature "Roshan, H., AZAl-Yaseri, M. Sarmadivalehand"). S.Iglauer (2016), On wettability of shale rocks. Journal of Colloid and Interface Science, 475, 104-111.”), the expression of wetting angle θ 1 between the surface of coal/shale and water is obtained without considering the effect of adsorbed gas formula, that is, the expression of the wetting angle θ 1 between the coal/shale surface that is not adsorbed by gas and water in the present invention:

Figure BDA0003572566430000081
Figure BDA0003572566430000081

式中,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气体的煤/页岩表面润湿角的表达式为:Step 4. For a constant temperature, the determined coal/shale surface and aqueous solution, combined with steps 2 and 3, update the expression of the wetting angle of the coal/shale surface with CO 2 gas adsorbed in the high-pressure CO 2 environment as:

Figure BDA0003572566430000082
Figure BDA0003572566430000082

式中θ是高压CO2环境中煤/页岩表面润湿角;对于煤/页岩表面和水溶液不变的恒温情况下,已知

Figure BDA0003572566430000083
是常数(由文献《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
Figure BDA0003572566430000083
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 formula 5, the wetting angle of coal/shale surface under different gas pressures can be obtained by obtaining f 1 , ρ g , ρ lf and the surface wetting angle of coal/shale under standard conditions (ie, normal pressure).

为了准确获取高压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 above formula 5, the present invention provides a measurement method, comprising:

步骤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内,得到

Figure BDA0003572566430000091
值;Step 403: Measure the wetting angle of the coal/shale surface under standard conditions, import it into the coal/shale wettability formula 5, and obtain
Figure BDA0003572566430000091
value;

步骤404、将所述不同CO2压力下煤/页岩表面和液体之间的接触面积占煤/页岩表面面积的比例f1

Figure BDA0003572566430000092
值导入所述煤/页岩润湿性公式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
Figure BDA0003572566430000092
The value is imported into the coal/shale wettability formula 5 to obtain the coal/shale surface wettability angle under different pressures.

本发明提供的测量方法中,将低场核磁共振技术与煤/页岩表面水接触润湿角测量技术结合起来,通过建立煤/页岩的润湿性与核磁共振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。根据公式

Figure BDA0003572566430000093
结合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
Figure BDA0003572566430000093
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、获取不同压力下f13. 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值为:

Figure BDA0003572566430000101
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:
Figure BDA0003572566430000101

4、测量标准状况下煤/页岩表面润湿角。4. Measure the wetting angle of coal/shale surface under standard conditions.

取煤样/页岩粉碎至大于200目,取煤粉用加压成型模具在50MPa压力下压制成具有光滑平面的试片,可制成圆柱体或椭圆体试片。在恒湿室内,用注射器往试片上滴落不足4mg的蒸馏水水滴,用快速相机进行拍照,通过图像分析,测量水滴与表面接触的长度a,通过单圆拟合水滴形状获取圆半径r,由公式

Figure BDA0003572566430000111
求出接触角θ.即标况下煤/页岩表面润湿角。利用公式4及已知的标况下煤/页岩表面润湿角和气液密度,此处液体为蒸馏水,固体表面液体膜密度ρlf为水密度取1g/cm3,求得
Figure BDA0003572566430000112
值。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
Figure BDA0003572566430000111
Obtain the contact angle θ. That is, the wetting angle of the coal/shale surface under standard conditions. Using formula 4 and the known wetting angle and gas-liquid density of the coal/shale surface under standard conditions, where the liquid is distilled water, the liquid film density ρ lf on the solid surface is the water density and takes 1g/cm 3 to obtain
Figure BDA0003572566430000112
value.

5、利用公式5,计算获取不同压力下的煤/页岩表面润湿角。5. Using formula 5, calculate and obtain the wetting angle of coal/shale surface under different pressures.

利用获得的不同压力下的f1的值和不同压力下CO2气体密度及步骤4中的

Figure BDA0003572566430000113
值,即可获得不同压力下煤/页岩样品的表面润湿角。Using the obtained values of f1 at different pressures and CO2 gas densities at different pressures and the values in step 4
Figure BDA0003572566430000113
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、获取不同压力下f13. 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值为:

Figure BDA0003572566430000114
继续向煤样中注入CO2气体,压力达到2MPa,静置10个小时,然后测试完全饱和水样品的T2谱图,记做T2-2,其中P1谱峰信号幅度为P1-2,2MPa压力下的f1值为:
Figure BDA0003572566430000121
依次类推直至压力增大到6MPa,6MPa压力下的f1值为:
Figure BDA0003572566430000122
依次计算得到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:
Figure BDA0003572566430000114
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:
Figure BDA0003572566430000121
And so on until the pressure increases to 6MPa, the value of f 1 under 6MPa pressure is:
Figure BDA0003572566430000122
The f 1 values are calculated sequentially as follows:

表1不同压力下P1峰幅度值及f1 Table 1 P1 peak amplitude and f 1 under different pressures

Figure BDA0003572566430000123
Figure BDA0003572566430000123

4、测量标准状况下煤/页岩表面润湿角。4. Measure the wetting angle of coal/shale surface under standard conditions.

取煤样粉碎至大于200目,取煤粉用加压成型模具在50MPa压力下压制成具有光滑平面的圆柱体试片。在恒湿室内,用注射器往试片上滴落不足4mg的蒸馏水水滴,用快速相机进行拍照,得到图5,通过图像分析,测量水滴与表面接触的长度a=1.94mm和半径r=0.99mm,由公式

Figure BDA0003572566430000124
求出表面接触角78.46°.利用公式5及已知的标况下煤/页岩表面润湿角和气液密度求得
Figure BDA0003572566430000125
值为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
Figure BDA0003572566430000124
Obtain the surface contact angle of 78.46°. Use formula 5 and the known standard conditions of coal/shale surface wetting angle and gas-liquid density to obtain
Figure BDA0003572566430000125
The value is 1.2079.

5、利用公式5,计算获取不同压力下煤/页岩样品表面润湿角。5. Using formula 5, calculate and obtain the surface wetting angle of coal/shale samples under different pressures.

利用获得的不同压力下的f1的值和不同压力下CO2气体密度及上述步骤中的

Figure BDA0003572566430000126
值,即可获得不同压力下煤/页岩样品表面润湿角,如表2所示。Using the obtained values of f1 at different pressures and CO2 gas density at different pressures and the above steps
Figure BDA0003572566430000126
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

压力(MPa)Pressure(MPa) 00 11 22 33 44 55 66 润湿角(°)Wetting angle (°) 7878 86.686.6 92.992.9 96.496.4 99.699.6 102.2102.2 106.3106.3

以上实施例仅为本申请的示例性实施例,不用于限制本申请,本申请的保护范围由权利要求书限定。本领域技术人员可以在本申请的实质和保护范围内,对本申请做出各种修改或等同替换,这种修改或等同替换也应视为落在本申请的保护范围内。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)

1. A method for calculating coal/shale wettability in a high-pressure environment is characterized by comprising the following steps:
step 1, high pressure CO2In the environment, CO2The gas is adsorbed on part of the coal/shale surface to replace part of the adsorbed water adsorbed on the coal/shale surface, and the CO adsorbed on the coal/shale surface2The gas being CO2Adsorbing gas, adsorbing water and the CO on the surface of the coal/shale2Adsorbing gas, said CO2With gas-adsorbing external surfaces in contact with water and coal/shale surfaces having mixed surfacesWettability characteristics;
step 2, obtaining high-pressure CO by combining the relation of the wettability of the liquid on the mixed surface and the conditions of the step 12In the environment, CO is adsorbed2A wettability expression of the coal/shale surface of the gas;
wherein the wettability relation of the liquid on the mixing surface is as follows:
Figure FDA0003572566420000011
wherein θ represents a mixing surface wetting angle; gamma rayLGRepresents the gas-liquid surface tension; gamma rayi,SG、γi,SLRespectively representing the interfacial tension of each material in the mixing surface with the gas and the liquid; i represents the material in the mixing surface; f. ofiRepresenting the proportion of the contact area of each material with the mixing surface in the mixing surface, f1+f2+…fnN denotes the total number of materials in the mixing surface as 1; thetaiIndicating the wetting angle between the liquid and each material;
high pressure CO2In the environment, the coal/shale surface is composed of coal/shale surface which is not adsorbed by gas and coal/shale surface which is adsorbed by gas, CO is adsorbed2The wettability expression of the coal/shale surface of the gas is:
cosθ=f1cosθ1+f2cosθ2 (2)
in the formula, theta represents high-pressure CO2The wetting angle of the coal/shale surface in the environment, f1 represents the proportion of the coal/shale surface which is not adsorbed with gas, namely the proportion of the contact area of the coal/shale surface and the adsorbed water in the coal/shale surface; theta1Represents the wetting angle between the coal/shale surface and water that is not adsorbed by the gas; f2 represents the ratio of the surface of coal adsorbed by gas, i.e., the CO2The contact area of the adsorption gas and the water accounts for the proportion of the surface of the coal/shale; theta.theta.2Indicating the wetting angle between the surface of the coal adsorbed by the gas and the water, i.e. the CO2The wetting angle between the adsorbed gas and water is pi.
Step 3, combining a Young equation of the wetting angle with the approximate theorem of "sharp-kink" to obtain an expression of the wetting angle between the surface of the coal/shale which is not adsorbed by the gas and the water:
Figure FDA0003572566420000021
in the formula, theta1Represents the wetting angle between the coal/shale surface and water that is not adsorbed by the gas, I represents the van der waals potential energy integral; Δ ρ represents a gas-liquid density difference; gamma rayLGRepresents the gas-liquid interfacial tension, ρgRepresents the gas density; rholfIndicating that the coal/shale surface liquid film density is equal to the liquid density.
Step 4, combining the steps 2 and 3, updating the high-pressure CO2In the environment, CO is adsorbed2The expression for the coal/shale surface wetting angle of the gas is:
Figure FDA0003572566420000022
in which theta is high pressure CO2Coal/shale surface wetting angle in the environment; for the constant temperature of the coal/shale surface and aqueous solution, it is known
Figure FDA0003572566420000023
Is a constant, f1Represents the proportion of the coal/shale surface which is not adsorbed with gas, namely the proportion of the contact area of the coal/shale surface and the adsorbed water in the coal/shale surface, rhogRepresenting high pressure CO2Gas density in the environment; ρ is a unit of a gradientlfIndicating that the liquid film density at the coal/shale surface is equal to the aqueous solution density.
2. The method of claim 1, wherein the updated high pressure CO of step 4 is used to calculate the coal/shale wettability2In the environment, CO is adsorbed2Method for measuring parameters in coal/shale wetting angle expression of gasThe method comprises the following steps:
step 401, measuring different CO's using low-field nuclear magnetic resonance2T2 spectrum of water in a saturated water coal sample under pressure;
step 402, analyzing the data of the T2 spectrogram to obtain different CO2The contact area between the coal/shale surface and the water under the pressure is calculated to obtain different CO2The ratio of the contact area between the coal/shale surface and the water under pressure to the coal/shale surface area;
step 403, measuring the wetting angle of the coal/shale surface under standard conditions, and introducing the updated high-pressure CO2In the environment, CO is adsorbed2Within the expression (5) of the coal/shale wetting angle of the gas, the result is
Figure FDA0003572566420000024
The value is obtained.
3. The method of claim 2, wherein the differential CO is different from the CO2The method for calculating the surface wetting angle of the coal/shale under pressure comprises the following steps: combining the different COs in step 4022The ratio of the contact area between the coal/shale surface and the water under pressure to the coal/shale surface area and the ratio obtained in step 403
Figure FDA0003572566420000031
Value introduction into the renewed high pressure CO2In the environment, CO is adsorbed2Within the expression (5) for the coal/shale wetting angle of the gas, different CO's were obtained2Coal/shale surface wetting angle under pressure.
4. The method of claim 2, wherein the different CO is calculated in step 402 by using the method of calculating the coal/shale wettability in the high pressure environment2The contact area between the coal/shale surface and the water under pressure is referred to the T2 spectrogram data of the water in a saturated water coal sample under standard conditions.
5. A high voltage as claimed in claim 2Method for calculating coal/shale wettability in environment, characterized in that, in step 402, said different CO2The formula of the ratio of the contact area between the coal/shale surface and the water under pressure to the coal/shale surface area is:
Figure FDA0003572566420000032
in the formula (f)1Representing the ratio of the contact area between the coal/shale surface and the water to the coal/shale surface area, P1-0Indicating the peak signal amplitude, P, of the first peak to the left in the spectrum of T2 under standard conditions1-nIndicating the peak signal amplitude of the first left peak in the spectrum of T2 at different pressures.
6. The method for calculating the coal/shale wettability in the high-pressure environment according to claim 2, wherein the method for measuring the wetting angle of the coal/shale surface under the standard condition in the step 403 is as follows:
taking a coal/shale sample, crushing to be larger than 200 meshes, and taking coal powder to prepare a test piece with a smooth plane; dripping less than 4mg of distilled water drops on the test piece in a constant humidity chamber; taking a picture, measuring the contact length of the water drop and the surface of the test piece and the radius of the water drop to obtain a contact angle; and further obtaining the coal/shale surface wetting angle under the standard condition.
7. The method for calculating the coal/shale wettability in the high-pressure environment according to claim 2, wherein the saturated water coal sample in the step 401 is prepared by the following steps:
and crushing the coal/shale sample to prepare a sample, drying the sample to constant weight, and putting a quantitative sample into the vacuumized saturated water to obtain a saturated water-coal sample.
CN202210329151.6A 2022-03-30 2022-03-30 Method for calculating wettability of coal/shale in high-pressure environment Active CN114705591B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (3)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
陈基瑜 等: "二氧化碳和氦气对低煤阶煤煤润湿性的影响研究", 《煤炭科学技术》, vol. 43, no. 11, 30 November 2015 (2015-11-30), pages 129 - 134 *

Cited By (6)

* Cited by examiner, † Cited by third party
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