CN111521532A - 一种致密储层含气饱和度的测定方法 - Google Patents
一种致密储层含气饱和度的测定方法 Download PDFInfo
- Publication number
- CN111521532A CN111521532A CN201910104816.1A CN201910104816A CN111521532A CN 111521532 A CN111521532 A CN 111521532A CN 201910104816 A CN201910104816 A CN 201910104816A CN 111521532 A CN111521532 A CN 111521532A
- Authority
- CN
- China
- Prior art keywords
- pore
- water
- saturation
- radius
- gas
- 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
- 238000000034 method Methods 0.000 title claims abstract description 61
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 306
- 239000011148 porous material Substances 0.000 claims abstract description 294
- 239000011435 rock Substances 0.000 claims abstract description 82
- 239000007789 gas Substances 0.000 claims description 71
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 claims description 44
- 229910052753 mercury Inorganic materials 0.000 claims description 44
- 238000002347 injection Methods 0.000 claims description 27
- 239000007924 injection Substances 0.000 claims description 27
- 239000012530 fluid Substances 0.000 claims description 20
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 18
- 239000002245 particle Substances 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 10
- 239000003345 natural gas Substances 0.000 claims description 9
- 238000009736 wetting Methods 0.000 claims description 7
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 239000008398 formation water Substances 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims description 3
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 2
- 239000011707 mineral Substances 0.000 claims description 2
- 239000011800 void material Substances 0.000 claims 2
- 238000005259 measurement Methods 0.000 abstract description 3
- 238000001179 sorption measurement Methods 0.000 abstract description 3
- 238000012821 model calculation Methods 0.000 abstract description 2
- 238000002474 experimental method Methods 0.000 description 6
- 238000005481 NMR spectroscopy Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 229920006395 saturated elastomer Polymers 0.000 description 2
- 238000009738 saturating Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000002209 hydrophobic effect Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000011068 loading method Methods 0.000 description 1
- 238000001225 nuclear magnetic resonance method Methods 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
- G01N33/241—Earth materials for hydrocarbon content
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
- E21B49/02—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells by mechanically taking samples of the soil
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
- G01N15/0886—Mercury porosimetry
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Pathology (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Analytical Chemistry (AREA)
- Geology (AREA)
- Dispersion Chemistry (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Remote Sensing (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Soil Sciences (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明提供了一种致密储层含气饱和度的测定方法。该方法包括以下步骤：测定致密储层岩石样品的孔径分布，根据孔径分布计算得到所述致密储层岩石样品的自由水饱和度；计算得到所述致密储层岩石样品的水膜水饱和度；计算得到所述致密储层岩石样品的角隅水饱和度；根据以下公式计算得到所述致密储层岩石样品的含气饱和度：Sg＝100‑Sw，其中，Sw为含水饱和度，单位为％；该含水饱和度为自由水饱和度、水膜水饱和度、角隅水饱和度之和；Sg为含气饱和度，单位为％。本发明的致密储层含气饱和度的测定方法采用模型计算，避免了实验过程中水挥发、表面吸附和水流观察等对含气饱和度测定结果造成的误差。
Description
技术领域
本发明涉及一种致密储层含气饱和度的测定方法，属于油气测定技术领域。
背景技术
地下岩石天然气含量通常采用体积法进行计算，需要含气饱和度、孔隙度、岩石密度、地层因子四个关键参数，其中，含气饱和度变化范围较大，目前实验室测定致密储层含气饱和度采用充注天然气的方法，在饱和水的岩石中注入天然气，保持一定的充注压力，观察出口不再有水流出，取出岩心，称量岩心实验前后减少的质量计算注入的气体体积，其与孔隙体积之比即为含气饱和度。改变不同注入压力，可以得到不同充注压力条件下的含气饱和度。
致密储层孔隙度低，容易受到实验条件影响，误差增大。具体来讲，该方法存在以下弊端：1、实验过程中液体计量不准确，致密储层本身孔隙度低，岩石含水量小，容易受到影响。如在装样、取样和称量过程时间较长，岩石中水分易挥发，造成含水量偏低；岩石表面吸附水容易使含水量偏大；2、以观察出口不再有水流出作为实验结束条件，不能观察管线中存在的水，充注出的水在管线中不能被看到，导致提前结束实验，含气饱和度偏小；3、测定多个压力点的含气饱和度，需要在每个压力点取出样品测定含气饱和度，然后再次饱和水，测定下一个压力点对应的饱和度，样品反复时间长、过程繁琐；4、岩石反复充注、饱和水容易改变岩石孔隙结构，如产生微裂缝等；5、润湿性岩石反复充注、饱和水容易改变岩石润湿性，如亲水变为疏水，容易形成优势通道，导致气体充注不完全；6、实验难以完全达到地下温压条件，流体压力小于70MPa，温度小于200℃；7、只能测定总含气饱和度单一数值。
发明内容
为解决上述技术问题，本发明的目的在于提供一种致密储层含气饱和度的测定方法，该方法能够克服现有技术的测定方法繁琐、结果不准确、难以满足地下高温高压地质条件等弊端。
为达到上述目的，本发明提供了一种致密储层含气饱和度的测定方法，其包括以下步骤：
步骤1：测定致密储层岩石样品的孔径分布，根据孔径分布计算得到所述致密储层岩石样品的自由水饱和度；
步骤2：计算得到所述致密储层岩石样品的水膜水饱和度；
步骤3：计算得到所述致密储层岩石样品的角隅水饱和度；
步骤4：根据以下公式计算得到所述致密储层岩石样品的含气饱和度：
Sg＝100-Sw
其中，Sw为含水饱和度，单位为％；该含水饱和度为自由水饱和度、水膜水饱和度、角隅水饱和度之和；
Sg为含气饱和度，单位为％。
本发明所提供的测定方法的前提是假设储层岩石中仅含有地层水和天然气两种流体，含气饱和度受到岩石含水饱和度的直接影响，得到储层含水饱和度，即可得到含气饱和度。形成岩石的矿物颗粒的大小及堆叠关系决定了流体储集空间的几何形态及空间分布。地层水在岩石中主要以三种形式赋存：(1)自由水：不能被天然气排驱的孔隙水，占据小孔喉不能被充注，其饱和度受天然气充注压力和孔喉毛管力的影响；(2)水膜水：以薄膜形式吸附在孔壁，其饱和度与孔隙壁表面积和水膜厚度有关；(3)角隅水：存在颗粒间不易驱替的角隅中，其饱和度受到岩石颗粒的大小和接触关系影响。
在上述方法中，测定致密储层岩石样品的孔径分布可以采用压汞法、核磁共振法等进行，但不限于以上两种方法。采用压汞法进行岩石孔径测试的基本原理是：汞不润湿岩石表面，是非润湿相，岩石孔隙中的空气是润湿相。汞注入岩石孔隙的过程是用非润湿相驱替润湿相。由于在孔隙喉道处毛管压力的作用，压力较低时，汞不能进入孔隙中。只有在注入压力超过毛管压力时，汞才可以进入岩样的孔隙喉道中。这时的注入压力即为毛管压力，毛管压力所对应的毛细管半径就相当于孔喉半径，而进入孔隙中汞的体积就是孔隙体积。把注入压力提高，使得汞可以进入更小的喉道控制的孔隙中。这样不断提高注入压力，就可以得到注入压力也就是毛管压力与孔隙体积的关系。由注入汞体积和岩样的孔隙体积可以计算出润湿相的饱和度，进而就可以得到不同孔喉半径区间的孔隙大小分布及毛管压力曲线。在步骤1中，所述孔径分布的数据包括不同孔喉半径对应的孔隙空间的数据。
在上述方法中，优选地，在步骤1中，所述不同孔喉半径对应的孔隙空间的数据包括不同半径孔隙空间对应的孔喉半径及所占的孔隙体积百分比。
在上述方法中，优选地，所述步骤1采用压汞法测定致密储层岩石样品的孔径分布按照以下方式进行：
将汞注入岩石孔隙之中，以驱替岩石孔隙中的空气；
当注入压力超过毛管压力时，汞进入岩样的孔隙喉道中，这时的注入压力即为毛管压力，压力所对应的毛细管半径就相当于孔喉半径，而进入孔隙中汞的体积就是孔隙体积；
提高注入压力，使得汞进入更小的喉道控制的孔隙，获得该更小的喉道对应的孔喉半径以及孔隙体积；其中，每一压力值下进行到汞不再进入岩石为止；
重复提高注入压力，获得不同的喉道对应的孔喉半径和孔隙体积。
在上述方法中，由于实验室的流体***与油藏***的表面张力不同，有必要将实验室的毛管压力数据转换为油藏毛管压力，假定“J”函数是岩石特性参数，从实验室到气藏都不发生变化，根据实际气藏条件，可以把压汞曲线换算成含水饱和度与充注压力关系(例如图3)。步骤1中根据孔径分布计算得到致密储层岩石样品的自由水饱和度可以按照以下方式进行：
(1)根据所述致密储层岩石样品中的不同的喉道对应的孔喉半径和孔隙体积得到注入压力与孔隙体积的关系：
其中：P汞-空气为毛细管压力，单位为106Pa；σ汞-空气为汞与空气的流体界面张力，单位为mN/m；θ汞-空气为汞与空气的润湿接触角，单位为°；r为毛细管半径，单位为10-6m；压汞实验中σ汞-空气＝480mN/m，θ＝140°；
(2)由注入汞体积和所述致密储层岩石样品的孔隙体积计算出空气的饱和度，进而得到不同孔喉半径区间的孔隙大小分布及毛管压力曲线：
σ气-水为界面张力，单位为mN/m；不同油藏温压条件不同，界面张力会有变化；
θ为气水接触角，优选地，θ一般为0，cosθ＝1；
P充注为充注压力，单位为MPa；
r为毛细管半径，单位为10-6m；
(3)将不同孔喉半径区间的孔隙大小分布及毛管压力曲线转换为自由水饱和度与毛管压力曲线，根据充注压力得到自由水饱和度，即根据充注压力计算得到相应的孔喉半径，将大于此半径的孔隙体积加和，该加和的值占所有孔隙体积的百分比即为自由水饱和度。
在上述方法中，优选地，所述步骤2按照以下方式进行：
(1)束缚水饱和度法将储集层中水膜的厚度视为一定值是不妥当的，同时油气藏中水膜的厚度也不可能成为储集层中有用孔喉半径的下限值(参见图1)。利用Dlvo理论，给出水膜厚度的计算公式，即按照以下公式计算得到不同孔喉半径对应的孔隙空间中的水膜厚度：
其中，a为常数，储集层越亲水，其值越大，优选取值为1.18×10-7；h为水膜厚度，单位为μm；σ为流体界面张力，单位为mN/m；r为孔喉半径，单位为10-6m；
(2)由上可以看出，均匀单管的油气驱水后会留下一层水膜，水膜是束缚水的一个组成部分。水膜的厚度会随着表面张力的增加而变薄，也会随孔径的变大而逐渐增大。将水膜所占体积与毛细管体积作比较，计算水膜束缚水饱和度的大小，即按照以下公式计算得到水膜饱和度：
其中，Swmi为第i种孔喉半径对应的孔隙空间的水膜水饱和度；
r为第i种孔喉半径对应的孔隙空间对应的孔喉半径；
h为第i种孔喉半径对应的孔隙空间对应的水膜厚度。
在上述方法中，储层角隅水饱和度受到孔隙形状的影响，对于理想情况，例如球形颗粒正排列、球形颗粒三角排列、正三角形孔隙、正方形孔隙、正六边形孔隙均可以计算含水饱和度数值。对于特定储层，角隅水饱和度计算还可采用较简便方法，通过观察储层孔隙结构，直接确定角隅水饱和度，或者统计不同孔隙结构类型所占比例，加权平均确定角隅水饱和度。优选地，所述步骤3按照以下方式之一进行：
(1)观察储层孔隙结构，直接确定角隅水饱和度；该方式适用于有条件进行图像观察的储层样品；
(2)统计不同孔隙结构类型对应的角隅水所占的比例，加权平均确定角隅水饱和度；
(3)颗粒表面水膜压力不稳定，将向角隅处汇集，在颗粒模型中大部分束缚水分布在颗粒间的角隅部位(参见图2)，形成角隅水。颗粒模型含水饱和度公式如下，按照以下公式可以计算得到角隅水饱和度：
其中：A2为气水界面的圆形截面积；
A3为角隅水的截面积；
r为气水界面圆弧半径，即孔喉半径，
单位为10-6m；
R为所述致密储层岩石样品的颗粒半径，单位为10-6m；
Swyi为第i种孔喉半径对应的孔隙空间对应的含水饱和度；
P充注为充注压力，单位为MPa；
σ为流体界面张力，单位为mN/m；
在上述方法中，优选地，所述步骤4按照以下方式进行：
(1)根据以下公式计算得到含水饱和度：
其中，Sw为含水饱和度；
λi为第i种孔喉半径对应的孔隙空间所占孔隙体积百分数；
Swyi为第i种孔喉半径对应的孔隙空间中角隅含水饱和度；
Swmi为第i种孔喉半径对应的孔隙空间中水膜含水饱和度；
Swz为自由水饱和度，即自由水孔隙空间所占孔隙体积比例；
(2)根据以下公式计算得到所述致密储层岩石样品的含气饱和度：
Sg＝100-Sw
其中，Sw为含水饱和度，单位为％；该含水饱和度为自由水饱和度、水膜水饱和度、角隅水饱和度之和；
Sg为含气饱和度，单位为％。
根据本发明的具体实施方案，优选地，本发明所提供的包括以下步骤：
步骤1、自由水饱和度：
采用压汞法测定致密储层岩石样品的孔径分布，当注入压力超过毛管压力时，汞进入岩样的孔隙喉道中，这时的注入压力即为毛管压力，毛管压力所对应的毛细管半径就相当于孔喉半径，而进入孔隙中汞的体积就是孔隙体积；
把注入压力提高，使得汞进入更小的喉道控制的孔隙中，得到相应的孔喉半径和孔隙体积；
不断提高注入压力，得到注入压力也就是毛管压力与孔隙体积的关系；
由注入汞体积和致密储层岩石样品的孔隙体积计算得到空气的饱和度，进而得到不同孔喉半径区间的孔隙大小分布及毛管压力曲线；
将致密储层岩石样品的孔径分布数据转换为油藏的毛管压力，结合实际气藏温压力条件的界面张力、接触角，将压汞曲线换算成含水饱和度随驱替压差变化的关系：
将致密储层岩石样品的孔径分布数据与毛管压力曲线转换为自由水饱和度与毛管压力曲线，得到充注压力，即可得到自由水饱和度，即根据充注压力计算得到相应的孔喉半径，将大于此半径的孔隙体积加和，该加和的值占所有孔隙体积的百分比即为自由水饱和度；
σ气-水为界面张力，单位为mN/m；优选地，气水界面张力取σ气-水＝20mN/m；
θ为气水接触角，优选地θ为0，cosθ＝1；
P充注为充注压力，单位为MPa；
r为毛细管半径，单位为10-6m；
步骤2、水膜水饱和度：
按照以下公式计算得到不同孔喉半径对应的孔隙空间中的水膜厚度：
其中，a为常数，储集层越亲水，其值越大，优选取值为1.18×10-7；h为水膜厚度，单位为μm；σ为流体界面张力，单位为mN/m；r为孔喉半径，单位为10-6m；
按照以下公式计算得到水膜饱和度：
其中，Swmi为第i种孔喉半径对应的孔隙空间的水膜水饱和度；
r为第i种孔喉半径对应的孔隙空间对应的孔喉半径；
h为第i种孔喉半径对应的孔隙空间对应的水膜厚度；
步骤3、角隅水饱和度：
按照以下公式计算得到角隅水饱和度：
其中：r为气水界面圆弧半径，即孔喉半径，
单位为10-6m；
R为所述致密储层岩石样品的颗粒半径，单位为10-6m；
Swyi为第i种孔喉半径对应的孔隙空间对应的含水饱和度；
P充注为充注压力，单位为MPa；
σ为流体界面张力，单位为mN/m；
步骤4、含气饱和度：
根据以下公式计算得到含水饱和度：
其中，Sw为含水饱和度；
λi为第i种孔喉半径对应的孔隙空间所占孔隙体积百分数；
Swyi为第i种孔喉半径对应的孔隙空间中角隅含水饱和度；
Swmi为第i种孔喉半径对应的孔隙空间中水膜含水饱和度；
Swz为自由水饱和度，即自由水孔隙空间所占孔隙体积比例。
再根据以下公式计算得到所述致密储层岩石样品的含气饱和度：
Sg＝100-Sw
其中，Sw为含水饱和度，单位为％；该含水饱和度为自由水饱和度、水膜水饱和度、角隅水饱和度之和；
Sg为含气饱和度，单位为％。
本发明所提供的致密储层含气饱和度的测定方法具有以下有益效果：
1、采用模型计算，避免了实验过程中水挥发、表面吸附和水流观察等对含气饱和度测定结果造成的误差；
2、避免了岩石反复充注对孔隙结构与润湿性的影响；
3、直接得到多个压力点的含气饱和度，简便、快速；
4、符合地下温压条件，流体压力可以大于70MPa，温度大于200℃；
5、可以分别得到自由水、水膜、角隅水这三种水的饱和度与比例；
6、不仅能测定总含气饱和度单一数值，还能得到不同孔径内的含气饱和度分布。
附图说明
图1为水膜水饱和度计算的示意图。
图2为颗粒模型角隅水计算的示意图。
图3为实施例1的充注压力与含气饱和度的关系图。
图4为实施例1的充注孔喉半径与含气饱和度的关系图。
图5为实施例2的充注压力与含气饱和度的关系图。
图6为实施例2的充注孔喉半径与含气饱和度的关系图。
具体实施方式
为了对本发明的技术特征、目的和有益效果有更加清楚的理解，现对本发明的技术方案进行以下详细说明，但不能理解为对本发明的可实施范围的限定。
实施例1
本实施例提供了一种致密储层含气饱和度的测定方法，该方法理论上可以用于不同类型岩石样品，本实施例采用某砂岩样品进行说明，不同样品输入参数不同，相应结果会不同。该方法包括以下步骤：
采用高压压汞测量样品的孔径分布，结果如表1所示。
表1某砂岩1样品孔径分布表
(1)自由水饱和度：
将孔径分布数据转换为油藏毛管压力，结合实际气藏温压力条件的界面张力、接触角等，把压汞曲线换算成含水饱和度随驱替压差关系：
将实验室测得的孔径分布数据与毛管压力曲线转换为自由水饱和度与毛管压力曲线，得到充注压力，即可得到自由水饱和度。
σ气-水为界面张力，单位为mN/m；本实施例中，气水界面张力取σ气-水＝20mN/m；
θ为气水接触角，θ为0，因此cosθ＝1；
P充注为充注压力，单位为MPa；
r为毛细管半径，单位为10-6m；
(2)水膜水饱和度：
按照以下公式计算得到不同孔喉半径对应的孔隙空间中的水膜厚度：
其中，hi为第i种孔喉半径对应的孔隙空间对应的水膜厚度，单位为μm；σ为流体界面张力，单位为mN/m；ri为第i种孔喉半径对应的孔隙空间对应的孔喉半径，在本实施例中，i为1-20，具体数据如表1所示；
按照以下公式计算得到水膜饱和度：
其中，Swmi为第i种孔喉半径对应的孔隙空间的水膜水饱和度；
ri为第i种孔喉半径对应的孔隙空间对应的孔喉半径，在本实施例中，i为1-20，具体数据如表1所示；
hi为第i种孔喉半径对应的孔隙空间对应的水膜厚度；
(3)角隅水饱和度：
按照以下公式计算得到角隅水饱和度(参见图2)：
其中：ri为气水界面圆弧半径，
单位为10-6m；
Ri为致密储层岩石样品的颗粒半径，单位为10-6m；
Swyi为第i种孔喉半径对应的孔隙空间对应的含水饱和度；
P充注i为第i种孔喉半径对应的孔隙空间对应的充注压力，单位为MPa；
σ为流体界面张力，单位为mN/m；
(4)含气饱和度：
根据以下公式计算得到含水饱和度：
其中，Sw为含水饱和度；
λi为第i种孔喉半径对应的孔隙空间所占孔隙体积百分数；
λz为自由水孔隙空间所占孔隙体积比例；
Swyi为第i种孔喉半径对应的孔隙空间中角隅含水饱和度；
Swmi为第i种孔喉半径对应的孔隙空间中水膜含水饱和度；
Swz为自由水饱和度，没有被充注，取100％；
再根据以下公式计算得到所述致密储层岩石样品的含气饱和度：
Sg＝100-Sw
其中，Sw为含水饱和度，单位为％；该含水饱和度为自由水饱和度、水膜水饱和度、角隅水饱和度之和；
Sg为含气饱和度，单位为％。
本实施例的各项参数如表2所示。图3为本实施例的充注压力与含气饱和度的关系图。图4为本实施例的充注孔喉半径与含气饱和度的关系图。
表2本实施例的参数表
通过以上内容可以看出，本实施例1所提供的致密储层含气饱和度的测定方法采用模型进行计算，避免了实验过程中水挥发、表面吸附和水流观察等对含气饱和度测定结果造成的误差，也避免了岩石反复充注对孔隙结构与润湿性的影响。该方法直接对21个压力点进行测定，直接得到多个压力点的含气饱和度，简便、快速。该方法不仅能测定总含气饱和度单一数值，还能得到不同孔径内的含气饱和度分布。
实施例2
本实施例提采用核磁共振测定该样品的孔径分布，通过核磁共振T2谱测定得到T2i弛豫时间及所占比例λi，通过公式计算孔喉半径ri＝3·ρ·T2i，其中，ρ为岩石表面弛豫率，单位为μm/ms，结果如表3所示。
表3某砂岩2样品孔径分布表
(1)自由水饱和度：
将孔径分布数据转换为油藏毛管压力，结合实际气藏温压力条件的界面张力、接触角等，把压汞曲线换算成含水饱和度随驱替压差关系：
将实验室测得的孔径分布数据与毛管压力曲线转换为自由水饱和度与毛管压力曲线，得到充注压力，即可得到自由水饱和度。
σ气-水为界面张力，单位为mN/m；本实施例中，气水界面张力取σ气-水＝50mN/m；
θ为气水接触角，θ为0，因此cosθ＝1；
P充注为充注压力，单位为MPa；
r为毛细管半径，单位为10-6m；
(2)水膜水饱和度：
按照以下公式计算得到不同孔喉半径对应的孔隙空间中的水膜厚度：
其中，hi为第i种孔喉半径对应的孔隙空间对应的水膜厚度，单位为μm；σ为流体界面张力，单位为mN/m；ri为第i种孔喉半径对应的孔隙空间对应的孔喉半径，在本实施例中，i为1-57，具体数据如表1所示；
按照以下公式计算得到水膜饱和度：
其中，Swmi为第i种孔喉半径对应的孔隙空间的水膜水饱和度；
ri为第i种孔喉半径对应的孔隙空间对应的孔喉半径，在本实施例中，i为1-57，具体数据如表1所示；
hi为第i种孔喉半径对应的孔隙空间对应的水膜厚度；
(3)角隅水饱和度：
按照以下公式计算得到角隅水饱和度：
(1)观察储层孔隙结构，直接确定角隅水饱和度；该方式适用于有条件进行图像观察的储层样品；
(2)根据表4，统计不同孔隙结构类型对应的角隅水所占的比例，加权平均确定角隅水饱和度为28.45：
Swyi＝31.2*75％+39.57*5％+21.5*10％+9.35*10％＝28.45％
表4某砂岩样品角隅水统计表
(4)含气饱和度：
根据以下公式计算得到含水饱和度：
其中，Sw为含水饱和度；
λi为第i种孔喉半径对应的孔隙空间所占孔隙体积百分数；
λz为自由水孔隙空间所占孔隙体积比例；
Swyi为第i种孔喉半径对应的孔隙空间中角隅含水饱和度；
Swmi为第i种孔喉半径对应的孔隙空间中水膜含水饱和度；
Swz为自由水饱和度，没有被充注，取100％；
再根据以下公式计算得到所述致密储层岩石样品的含气饱和度：
Sg＝100-Sw
其中，Sw为含水饱和度，单位为％；该含水饱和度为自由水饱和度、水膜水饱和度、角隅水饱和度之和；
Sg为含气饱和度，单位为％。
本实施例的各项参数如表5所示。图5为本实施例的充注压力与含气饱和度的关系图。图6为本实施例的充注孔喉半径与含气饱和度的关系图。
表5本实施例的参数表
通过以上内容可以看出，本实施例2所提供的致密储层含气饱和度的测定方法采用模型进行计算，与实施例1相比，砂岩样品的孔喉半径分布有核磁共振实验测得，角隅水的比例通过对孔隙类型观察结合加权计算获得。为了验证方法的有效性，与实际的物理模拟实验对比，充注结果与本方法计算得到结果一致物理模拟实验数据偏低的原因是由于实验时间有限，不能达到实际地下气藏充注的长时间，充注不完全造成结果偏低。但是由于物理模拟实验只能达到较低充注压力，故只能在较低充注压力范围验证，但是上述内容足以证明本发明的技术方案所具有的技术效果。
Claims (10)
1.一种致密储层含气饱和度的测定方法，其包括以下步骤：
步骤1：测定致密储层岩石样品的孔径分布，根据孔径分布计算得到所述致密储层岩石样品的自由水饱和度；
步骤2：计算得到所述致密储层岩石样品的水膜水饱和度；
步骤3：计算得到所述致密储层岩石样品的角隅水饱和度；
步骤4：根据以下公式计算得到所述致密储层岩石样品的含气饱和度：
Sg＝100-Sw
其中，Sw为含水饱和度，单位为％；该含水饱和度为自由水饱和度、水膜水饱和度、角隅水饱和度之和；
Sg为含气饱和度，单位为％。
2.根据权利要求1所述的方法，其中，在步骤1中，所述孔径分布的数据包括不同孔喉半径对应的孔隙空间的数据。
3.根据权利要求2所述的方法，其中，在步骤1中，所述不同孔喉半径对应的孔隙空间的数据包括不同半径孔隙空间对应的孔喉半径及所占的孔隙体积百分比。
4.根据权利要求1-3任一项所述的方法，其中，所述步骤1采用压汞法测定致密储层岩石样品的孔径分布按照以下方式进行：
将汞注入岩石孔隙之中，以驱替岩石孔隙中的空气；
当注入压力超过毛管压力时，汞进入岩样的孔隙喉道中，这时的注入压力即为毛管压力，压力所对应的毛细管半径就相当于孔喉半径，而进入孔隙中汞的体积就是孔隙体积；
提高注入压力，使得汞进入更小的喉道控制的孔隙，获得该更小的喉道对应的孔喉半径以及孔隙体积；
重复提高注入压力，获得不同的喉道对应的孔喉半径和孔隙体积。
5.根据权利要求4所述的方法，其中，所述步骤1中根据孔径分布计算得到所述致密储层岩石样品的自由水饱和度按照以下方式进行：
(1)根据所述致密储层岩石样品中的不同的喉道对应的孔喉半径和孔隙体积得到注入压力与孔隙体积的关系：
其中：P汞-空气为毛细管压力，单位为106Pa；σ汞-空气为汞与空气的流体界面张力，单位为mN/m；θ汞-空气为汞与空气的润湿接触角，单位为o；r为毛细管半径，单位为10-6m；压汞实验中σ汞-空气＝480mN/m，θ汞-空气＝140o；
(2)结合实际气藏温压力条件的界面张力、接触角等，将实验室的孔喉分布数据转换为油藏条件的气水充注压力P气-水：
σ气-水为界面张力，单位为mN/m；优选地，气水界面张力取σ气-水＝20mN/m；
θ气-水为气水接触角，优选地θ为0，cosθ＝1；
P充注为气水充注压力，单位为MPa；
r为毛细管半径，单位为10-6m；
(3)将不同孔喉半径区间的孔隙大小分布及毛管压力曲线转换为自由水饱和度与毛管压力关系曲线，根据充注压力得到自由水饱和度，即根据充注压力计算得到相应的孔喉半径，将大于此半径的孔隙体积加和，该加和的值占所有孔隙体积的百分比即为自由水饱和度。
6.根据权利要求1所述的方法，其中，所述步骤2按照以下方式进行：
(1)按照以下公式计算得到不同孔喉半径对应的孔隙空间中的水膜厚度：
其中，a为常数，储集层越亲水，其值越大，优选取值为1.18×10-7；h为水膜厚度，单位为μm；σ为流体界面张力，单位为mN/m；r为孔喉半径，单位为10-6m；
(2)按照以下公式计算得到水膜饱和度：
其中，Swmi为第i种孔喉半径对应的孔隙空间的水膜水饱和度；
r为第i种孔喉半径对应的孔隙空间对应的孔喉半径；
h为第i种孔喉半径对应的孔隙空间对应的水膜厚度。
7.根据权利要求1所述的方法，其中，所述步骤3按照以下方式之一进行：
(1)观察储层孔隙结构，直接确定角隅水饱和度；
(2)统计不同孔隙结构类型对应的角隅水所占的比例，加权平均确定角隅水饱和度；
(3)按照以下公式计算得到角隅水饱和度：
其中：r为气水界面圆弧半径，即孔喉半径，
单位为10-6m；
R为所述致密储层岩石样品的颗粒半径，单位为10-6m；
Swyi为第i种孔喉半径对应的孔隙空间对应的含水饱和度；
P充注为充注压力，单位为MPa；
σ为流体界面张力，单位为mN/m。
8.根据权利要求1所述的方法，其中，所述步骤4按照以下方式进行：
(1)根据以下公式计算得到含水饱和度：
其中，Sw为含水饱和度；
λi为第i种孔喉半径对应的孔隙空间所占孔隙体积百分数；
Swyi为第i种孔喉半径对应的孔隙空间中角隅含水饱和度；
Swmi为第i种孔喉半径对应的孔隙空间中水膜含水饱和度；
Swz为自由水饱和度，即自由水孔隙空间所占孔隙体积比例；
(2)根据以下公式计算得到所述致密储层岩石样品的含气饱和度：
Sg＝100-Sw
其中，Sw为含水饱和度，单位为％；该含水饱和度为自由水饱和度、水膜水饱和度、角隅水饱和度之和；
Sg为含气饱和度，单位为％。
9.根据权利要求1-8任一项所述的方法，其中，该方法包括以下步骤：
步骤1、自由水饱和度：
采用压汞法测定致密储层岩石样品的孔径分布，当注入压力超过毛管压力时，汞进入岩样的孔隙喉道中，这时的注入压力即为毛管压力，毛管压力所对应的毛细管半径就相当于孔喉半径，而进入孔隙中汞的体积就是孔隙体积；
把注入压力提高，使得汞进入更小的喉道控制的孔隙中，得到相应的孔喉半径和孔隙体积；
不断提高注入压力，得到注入压力也就是毛管压力与孔隙体积的关系；
由注入汞体积和致密储层岩石样品的孔隙体积计算得到空气的饱和度，进而得到不同孔喉半径区间的孔隙大小分布及毛管压力曲线；
将致密储层岩石样品的孔径分布数据转换为油藏的毛管压力，结合实际气藏温压力条件的界面张力、接触角，将压汞曲线换算成含水饱和度随驱替压差变化的关系：
将致密储层岩石样品的孔径分布数据与毛管压力曲线转换为自由水饱和度与毛管压力曲线，得到充注压力，即可得到自由水饱和度，即根据充注压力计算得到相应的孔喉半径，将大于此半径的孔隙体积加和，该加和的值占所有孔隙体积的百分比即为自由水饱和度；
σ气-水为界面张力，单位为mN/m；优选地，气水界面张力取σ气-水＝20mN/m；
θ为气水接触角，优选地θ为0，cosθ＝1；
P充注为充注压力，单位为MPa；
r为毛细管半径，单位为10-6m；
步骤2、水膜水饱和度：
按照以下公式计算得到不同孔喉半径对应的孔隙空间中的水膜厚度：
其中，a为常数，储集层越亲水，其值越大，优选取值为1.18×10-7h为水膜厚度，单位为μm；σ为流体界面张力，单位为mN/m；r为孔喉半径，单位为10-6m
按照以下公式计算得到水膜饱和度：
其中，Swmi为第i种孔喉半径对应的孔隙空间的水膜水饱和度；
r为第i种孔喉半径对应的孔隙空间对应的孔喉半径；
h为第i种孔喉半径对应的孔隙空间对应的水膜厚度；
步骤3、角隅水饱和度：
按照以下公式计算得到角隅水饱和度：
其中：r为气水界面圆弧半径，即孔喉半径，
单位为10-6m；
R为所述致密储层岩石样品的颗粒半径，单位为10-6m；
Swyi为第i种孔喉半径对应的孔隙空间对应的含水饱和度；
P充注为充注压力，单位为MPa；
σ为流体界面张力，单位为mN/m；
步骤4、含气饱和度：
根据以下公式计算得到含水饱和度：
其中，Sw为含水饱和度；
λi为第i种孔喉半径对应的孔隙空间所占孔隙体积百分数；
Swyi为第i种孔喉半径对应的孔隙空间中角隅含水饱和度；
Swmi为第i种孔喉半径对应的孔隙空间中水膜含水饱和度；
Swz为自由水饱和度，即自由水孔隙空间所占孔隙体积比例；
再根据以下公式计算得到所述致密储层岩石样品的含气饱和度：
Sg＝100-Sw
其中，Sw为含水饱和度，单位为％；该含水饱和度为自由水饱和度、水膜水饱和度、角隅水饱和度之和；
Sg为含气饱和度，单位为％。
10.根据权利要求1-9任一项所述的方法，其中，该方法基于以下前提：假设储层岩石中仅含有地层水和天然气两种流体，含气饱和度受到岩石含水饱和度的直接影响，得到储层含水饱和度，即可得到含气饱和度；
形成岩石的矿物颗粒的大小及堆叠关系决定了流体储集空间的几何形态及空间分布；
地层水在岩石中主要以三种形式赋存：(1)自由水：不能被天然气排驱的孔隙水，占据小孔喉不能被充注，其饱和度受天然气充注压力和孔喉毛管力的影响；(2)水膜水：以薄膜形式吸附在孔壁，其饱和度与孔隙壁表面积和水膜厚度有关；(3)角隅水：存在颗粒间不易驱替的角隅中，其饱和度受到岩石颗粒的大小和接触关系影响。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910104816.1A CN111521532B (zh) | 2019-02-01 | 2019-02-01 | 一种致密储层含气饱和度的测定方法 |
| US16/776,773 US11268944B2 (en) | 2019-02-01 | 2020-01-30 | Method for determining gas saturation of tight reservoir |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910104816.1A CN111521532B (zh) | 2019-02-01 | 2019-02-01 | 一种致密储层含气饱和度的测定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111521532A true CN111521532A (zh) | 2020-08-11 |
| CN111521532B CN111521532B (zh) | 2023-02-10 |
Family
ID=71835874
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910104816.1A Active CN111521532B (zh) | 2019-02-01 | 2019-02-01 | 一种致密储层含气饱和度的测定方法 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US11268944B2 (zh) |
| CN (1) | CN111521532B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113029908A (zh) * | 2021-03-16 | 2021-06-25 | 中国石油大学(华东) | 一种致密储层饱和度指数的实验室测量方法 |
| CN113514382A (zh) * | 2021-04-25 | 2021-10-19 | 西南石油大学 | 评价含膨胀粘土气藏岩石水相渗吸返排后水膜厚度的方法 |
| CN114037805A (zh) * | 2021-10-26 | 2022-02-11 | 重庆科技学院 | 一种非常规储层饱和度模型构建方法 |
| CN114235879A (zh) * | 2021-12-22 | 2022-03-25 | 重庆大学 | 应力加载条件下测试盖层润湿性与毛细管压力的方法 |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112129802B (zh) * | 2020-10-22 | 2023-08-15 | 重庆科技学院 | 一种水化页岩不同尺度孔隙体积增量的定量分析方法 |
| CN112377179A (zh) * | 2020-11-26 | 2021-02-19 | 中国石油天然气股份有限公司 | 一种基于分形理论的致密砂岩储层束缚水饱和度预测方法 |
| CN113155694A (zh) * | 2020-12-11 | 2021-07-23 | 中国石油天然气股份有限公司 | 一种基于恒压压汞实验的储层流动孔隙度确定方法 |
| CN112964604A (zh) * | 2021-03-25 | 2021-06-15 | 西南石油大学 | 一种考虑含水饱和度和吸附气的岩石接触角测试装置与方法 |
| CN114384227A (zh) * | 2021-06-16 | 2022-04-22 | 西安石油大学 | 一种致密砂岩的准连续溶蚀-胶结作用致密机理测定方法 |
| CN113484222B (zh) * | 2021-08-06 | 2022-09-13 | 中国地质大学（北京） | 致密气储层渗透率瓶颈区确定方法及相关设备 |
| CN113790997A (zh) * | 2021-09-02 | 2021-12-14 | 东北石油大学 | 基于梯度热处理技术测量岩心束缚流体饱和度的方法 |
| CN113740515B (zh) * | 2021-11-05 | 2022-03-01 | 中国科学院地质与地球物理研究所 | 综合表征深层海相碳酸盐岩油藏赋存下限临界条件的方法 |
| CN115163058A (zh) * | 2022-07-28 | 2022-10-11 | 中国地质大学（北京） | 基于核磁共振实验评价致密气储层液相滞留程度的方法 |
| CN115479881B (zh) * | 2022-10-25 | 2024-02-02 | 中国电建集团西北勘测设计研究院有限公司 | 一种测量砂卵砾石地层孔隙比的方法 |
| CN117347417B (zh) * | 2023-11-06 | 2024-05-03 | 西南石油大学 | 一种页岩表面弛豫率的计算方法 |
| CN117725445B (zh) * | 2024-02-07 | 2024-05-17 | 中国石油大学(华东) | 基于电导率的砂岩储层超临界二氧化碳饱和度计算方法 |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070089489A1 (en) * | 2005-10-11 | 2007-04-26 | Millipore Corporation | Methods and systems for integrity testing of porous materials |
| CN101929973A (zh) * | 2009-06-22 | 2010-12-29 | 中国石油天然气股份有限公司 | 裂缝储层含油气饱和度定量计算方法 |
| CN102174888A (zh) * | 2011-03-03 | 2011-09-07 | 康志勇 | 一种基于储集岩岩性参数的地层数据处理方法 |
| CN102434152A (zh) * | 2011-12-05 | 2012-05-02 | 中国石油天然气股份有限公司 | 一种储层含油饱和度的计算方法 |
| CN104729972A (zh) * | 2015-04-10 | 2015-06-24 | 中国石油大学(北京) | 确定碳酸盐岩气水相对渗透率的方法及装置 |
| US20150198036A1 (en) * | 2014-01-13 | 2015-07-16 | Schlumberger Technology Corporation | Method for estimating irreducible water saturation from mercury injection capillary pressure |
| US9139364B1 (en) * | 2012-07-25 | 2015-09-22 | Sandia Corporation | Method for carbon dioxide sequestration |
| CN105445441A (zh) * | 2015-12-03 | 2016-03-30 | 中国石油天然气股份有限公司 | 储层含气饱和度确定方法及装置 |
| CN105466830A (zh) * | 2015-11-19 | 2016-04-06 | 中国石油天然气股份有限公司 | 储层砂岩孔隙喉道尺寸分布识别方法 |
| CN105781539A (zh) * | 2016-03-15 | 2016-07-20 | 中国石油大学(华东) | 一种致密油气储层饱和度测井计算方法 |
| CN105986815A (zh) * | 2015-02-27 | 2016-10-05 | 中国石油化工股份有限公司 | 一种用于识别页岩地层地质甜点的方法 |
| CN108593493A (zh) * | 2018-01-10 | 2018-09-28 | 中国石油天然气股份有限公司 | 一种确定致密储层含气量的方法和装置 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101806215B (zh) | 2010-03-05 | 2013-05-01 | 中国石油集团川庆钻探工程有限公司 | 用束缚水饱和度数据判别储层流体类型的方法 |
| CA2862951C (en) * | 2011-10-21 | 2018-10-30 | Saudi Arabian Oil Company | Methods, computer readable medium, and apparatus for determining well characteristics and pore architecture utilizing conventional well logs |
| KR102145236B1 (ko) * | 2012-06-01 | 2020-08-18 | 인튜어티브 서지컬 오퍼레이션즈 인코포레이티드 | 영-공간을 이용하여 매니퓰레이터 암 대 환자의 충돌을 회피하는 방법 및 시스템 |
| CN103513270B (zh) | 2012-06-21 | 2017-06-09 | 中国石油天然气集团公司 | 一种基于岩石声学特性的气层识别评价方法及装置 |
| GB2521317B (en) | 2012-11-16 | 2018-10-31 | Halliburton Energy Services Inc | Triple phase evaluation of formation fluids |
| CN104358565B (zh) | 2014-10-31 | 2019-03-22 | 中国石油化工股份有限公司 | 一种气藏含气饱和度的确定方法 |
-
2019
- 2019-02-01 CN CN201910104816.1A patent/CN111521532B/zh active Active
-
2020
- 2020-01-30 US US16/776,773 patent/US11268944B2/en active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070089489A1 (en) * | 2005-10-11 | 2007-04-26 | Millipore Corporation | Methods and systems for integrity testing of porous materials |
| CN101929973A (zh) * | 2009-06-22 | 2010-12-29 | 中国石油天然气股份有限公司 | 裂缝储层含油气饱和度定量计算方法 |
| US20120109603A1 (en) * | 2009-06-22 | 2012-05-03 | Ning Li | Quantitative calculation method for oil (gas) saturation of fractured reservoir |
| CN102174888A (zh) * | 2011-03-03 | 2011-09-07 | 康志勇 | 一种基于储集岩岩性参数的地层数据处理方法 |
| CN102434152A (zh) * | 2011-12-05 | 2012-05-02 | 中国石油天然气股份有限公司 | 一种储层含油饱和度的计算方法 |
| US9139364B1 (en) * | 2012-07-25 | 2015-09-22 | Sandia Corporation | Method for carbon dioxide sequestration |
| US20150198036A1 (en) * | 2014-01-13 | 2015-07-16 | Schlumberger Technology Corporation | Method for estimating irreducible water saturation from mercury injection capillary pressure |
| CN105986815A (zh) * | 2015-02-27 | 2016-10-05 | 中国石油化工股份有限公司 | 一种用于识别页岩地层地质甜点的方法 |
| CN104729972A (zh) * | 2015-04-10 | 2015-06-24 | 中国石油大学(北京) | 确定碳酸盐岩气水相对渗透率的方法及装置 |
| CN105466830A (zh) * | 2015-11-19 | 2016-04-06 | 中国石油天然气股份有限公司 | 储层砂岩孔隙喉道尺寸分布识别方法 |
| CN105445441A (zh) * | 2015-12-03 | 2016-03-30 | 中国石油天然气股份有限公司 | 储层含气饱和度确定方法及装置 |
| CN105781539A (zh) * | 2016-03-15 | 2016-07-20 | 中国石油大学(华东) | 一种致密油气储层饱和度测井计算方法 |
| CN108593493A (zh) * | 2018-01-10 | 2018-09-28 | 中国石油天然气股份有限公司 | 一种确定致密储层含气量的方法和装置 |
Non-Patent Citations (4)
| Title |
|---|
| DE-HUA HAN 等: "Fizz water and low gas-saturated reservoirs", 《THE LEADING EDGE 》 * |
| 李昂 等: "致密砂岩气成藏机制类型及特征研究", 《吉林地质》 * |
| 陈思宇 等: "致密储层样品体积测量对孔隙度误差的影响", 《石油实验地质》 * |
| 马行陟 等: "页岩吸附气含量测定的影响因素定量分析", 《天然气地球科学》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113029908A (zh) * | 2021-03-16 | 2021-06-25 | 中国石油大学(华东) | 一种致密储层饱和度指数的实验室测量方法 |
| CN113029908B (zh) * | 2021-03-16 | 2021-11-26 | 中国石油大学(华东) | 一种致密储层饱和度指数的实验室测量方法 |
| CN113514382A (zh) * | 2021-04-25 | 2021-10-19 | 西南石油大学 | 评价含膨胀粘土气藏岩石水相渗吸返排后水膜厚度的方法 |
| CN113514382B (zh) * | 2021-04-25 | 2024-04-05 | 西南石油大学 | 评价含膨胀粘土气藏岩石水相渗吸返排后水膜厚度的方法 |
| CN114037805A (zh) * | 2021-10-26 | 2022-02-11 | 重庆科技学院 | 一种非常规储层饱和度模型构建方法 |
| CN114235879A (zh) * | 2021-12-22 | 2022-03-25 | 重庆大学 | 应力加载条件下测试盖层润湿性与毛细管压力的方法 |
| CN114235879B (zh) * | 2021-12-22 | 2024-04-30 | 重庆大学 | 应力加载条件下测试盖层润湿性与毛细管压力的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111521532B (zh) | 2023-02-10 |
| US11268944B2 (en) | 2022-03-08 |
| US20200249216A1 (en) | 2020-08-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111521532B (zh) | 一种致密储层含气饱和度的测定方法 | |
| Wyllie et al. | Fluid flow through unconsolidated porous aggregates | |
| Mauran et al. | Application of the Carman–Kozeny correlation to a high‐porosity and anisotropic consolidated medium: the compressed expanded natural graphite | |
| CN105866009B (zh) | 一种计算致密油储层有效孔隙度的方法及装置 | |
| Geffen et al. | Experimental investigation of factors affecting laboratory relative permeability measurements | |
| CN104330409B (zh) | 化学驱油用起泡剂体系的泡沫性能测量装置及其定量评价方法 | |
| CN114136861B (zh) | 一种储气库近井地带干化盐析效应实验***及评价方法 | |
| CN113431537B (zh) | 一种非稳态变流速大尺度岩心水驱气相对渗透率测试方法 | |
| CN109632594A (zh) | 一种基于高压压汞多尺度表征致密储层孔喉特征的方法 | |
| US20220291157A1 (en) | Method for establishing mathematical model of relationship between spontaneous imbibition volume and time of porous medium | |
| Liu et al. | Coupling lattice Boltzmann and finite volume methods for simulation of reactive transport | |
| CN109973063A (zh) | 用于确定碳化水驱油过程中碳化水对储层伤害程度的方法 | |
| CN113358542A (zh) | 用于测试不同孔喉尺寸范围内流体驱替效率的装置及方法 | |
| CN115559715A (zh) | 一种超高压低渗透气藏储层产水评价方法 | |
| Zhou et al. | The determination of surface relaxivity and application to coal spontaneous imbibition | |
| CN112147049B (zh) | 一种岩心水膜厚度的确定方法 | |
| CN209821099U (zh) | 基于核磁共振的多功能致密气储层动态参数联测装置 | |
| Li et al. | Quantitative experimental investigation of multiple PT effects on primary drainage process during scCO2 storage in deep saline aquifers | |
| CN112485282A (zh) | 含气体水合物沉积物土水特征曲线的测量***及其方法 | |
| CN113109234B (zh) | 重水核磁束缚水探测的低含油饱和度渗流规律修正方法 | |
| CN113984589A (zh) | 一种计算岩石迂曲度和气体扩散系数的方法 | |
| Al Sayari | The influence of wettability and carbon dioxide injection on hydrocarbon recovery | |
| CN110501276B (zh) | 一种针对相对渗透率曲线乳化失准的测定方法 | |
| CN112014288B (zh) | 气驱水相渗评估方法、装置及*** | |
| CN113916748A (zh) | 轻质油测页岩基质渗透率和采收率的装置和方法 |
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 |