CN111693425B - 一种基于压汞曲线的岩心薄膜束缚水饱和度的测定方法 - Google Patents
一种基于压汞曲线的岩心薄膜束缚水饱和度的测定方法 Download PDFInfo
- Publication number
- CN111693425B CN111693425B CN201910192711.6A CN201910192711A CN111693425B CN 111693425 B CN111693425 B CN 111693425B CN 201910192711 A CN201910192711 A CN 201910192711A CN 111693425 B CN111693425 B CN 111693425B
- Authority
- CN
- China
- Prior art keywords
- mercury
- saturation
- film
- pore throat
- formula
- 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.)
- Expired - Fee Related
Links
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 title claims abstract description 170
- 229910052753 mercury Inorganic materials 0.000 title claims abstract description 170
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 97
- 239000011435 rock Substances 0.000 title claims abstract description 31
- 238000000034 method Methods 0.000 title claims abstract description 25
- 239000011148 porous material Substances 0.000 claims abstract description 69
- 238000002474 experimental method Methods 0.000 claims abstract description 9
- 238000012360 testing method Methods 0.000 claims description 41
- 238000009736 wetting Methods 0.000 claims description 8
- 238000003556 assay Methods 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 21
- 238000011156 evaluation Methods 0.000 abstract description 3
- 238000011160 research Methods 0.000 abstract description 3
- 238000005259 measurement Methods 0.000 abstract description 2
- 230000007246 mechanism Effects 0.000 abstract description 2
- 239000010408 film Substances 0.000 description 44
- 239000012528 membrane Substances 0.000 description 5
- 239000010409 thin film Substances 0.000 description 4
- 239000002245 particle Substances 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000009969 flowable effect Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
Images
Classifications
-
- 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
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了基于压汞曲线的岩心薄膜束缚水饱和度的测定方法。所述方法包括如下步骤:对岩心进行压汞实验,得到不同进汞压力下的进汞饱和度,进而得到岩心的孔喉半径和孔吼分布频率;根据进汞压力得到岩心的薄膜束缚水厚度;根据孔喉半径得到岩心的不同进汞区间孔喉平均半径;根据孔喉分布频率、薄膜束缚水厚度和不同进汞区间孔喉平均半径得到岩心的不同进汞区间薄膜束缚水饱和度,进而得到岩心的薄膜束缚水饱和度。本发明可以根据实验压汞曲线数据,快速得到岩心薄膜束缚水饱和度,具有简便实用、可操作性强的特点。本发明薄膜束缚水饱和度的测定结果可应用于储层评价、产能预测、储量计算以及油气藏出水机理研究等众多方面,具有广阔的应用前景。
Description
技术领域
本发明涉及一种基于压汞曲线的岩心薄膜束缚水饱和度的测定方法,属于油气田生产技术领域。
背景技术
束缚水是滞留在储层岩石微小孔隙或吸附于岩石颗粒表面的不可流动的水。束缚水饱和度为束缚水所占孔隙的体积与岩石总孔隙体积之比。束缚水饱和度是油藏单相与油水两相渗流的临界参数,也是进行储层评价、产能预测、油水渗流规律分析和储量计算的基础参数。
束缚水包括两种形式:一种是滞留在储层岩石微小孔隙中的毛细管束缚水(如图1a所示);另一种是吸附于岩石颗粒表面(或大孔隙壁上)的薄膜束缚水(如图1b所示)。因此,束缚水饱和度也可分为两种:即毛细管束缚水饱和度和薄膜束缚水饱和度。目前,束缚水饱和度主要通过压汞、离心毛管力、油水相渗、核磁共振等流体实验方法测定。由于这些实验手段不能对毛细管束缚水和薄膜束缚水进行有效区分,致使目前薄膜束缚水饱和度尚无法准确定量的测定,这在很大程度上制约了薄膜束缚水的理论研究和实践应用。因此,建立一种薄膜束缚水饱和度的测定方法具有重要的理论意义和实践价值。
发明内容
本发明的目的是提供一种基于压汞曲线的岩心薄膜束缚水饱和度的测定方法,具有简便实用、可操作性强的特点。
本发明所提供的岩心薄膜束缚水饱和度的测定方法,包括如下步骤:
1)对岩心进行压汞实验,得到不同进汞压力下的进汞饱和度,进而得到岩心的孔喉半径和孔吼分布频率;
2)根据所述进汞压力得到岩心的薄膜束缚水厚度;
3)根据所述孔喉半径得到岩心的不同进汞区间孔喉平均半径;
4)根据所述孔喉分布频率、所述薄膜束缚水厚度和所述不同进汞区间孔喉平均半径得到岩心的不同进汞区间薄膜束缚水饱和度,进而得到岩心的薄膜束缚水饱和度。
上述的测定方法中,步骤1)中,所述压汞实验过程中,通过测量不同进汞压力下的进汞饱和度,便可绘制压汞曲线,如图2所示,进而确定岩心孔喉分布,具体地,可按照下述公式得到:
根据所述进汞压力,结合式(1)得到所述孔喉半径;
式(1)中,ri表示第i个进汞压力对应的孔喉半径,μm;σ表示汞表面张力,为0.48N/m;θ表示汞润湿角,为140°;Pi表示第i个测试点进汞压力,MPa;n表示进汞压力测试总点数。
根据所述进汞饱和度,结合式(2)得到所述孔喉分布频率;
fi=ΔSi=Si+1-Si (i=1,2,……,n-1;) (2)
式(2)中,fi表示第i个区间孔喉分布频率,%;ΔSi表示第i个区间汞饱和度增量,%;Si+1表示第i+1个测试点进汞压力对应的汞饱和度,%;Si表示第i个测试点进汞压力对应的汞饱和度,%;n表示进汞压力测试总点数。
上述的测定方法中,步骤2)中,根据式(3)得到所述薄膜束缚水厚度;
式(3)中,Pn表示第n个测试点进汞压力,MPa;rn表示束缚水薄膜厚度,μm;σ表示汞表面张力,为0.48N/m;θ表示汞润湿角,为140°。
上述的测定方法中,步骤3)中,根据式(4)得到所述不同进汞区间孔喉平均半径;
式(4)中,rci表示第i个进汞区间孔喉平均半径,μm;ri表示第i个测试点进汞压力对应的孔喉半径,μm;ri+1表示第i+1个测试点进汞压力对应的孔喉半径,μm;n表示进汞压力测试总点数。
上述的测定方法中,步骤4)中,建立了薄模束缚水同心毛细管模型,如图3所示,进而根据式(5)得到所述不同进汞区间薄膜束缚水饱和度;
式(5)中,Swni表示第i个进汞区间薄膜束缚水饱和度,%;rci表示第i个进汞区间孔喉平均半径,μm;rn表示束缚水薄膜厚度,μm;ΔSi表示第i个进汞区间汞饱和度增量,即第i个区间孔喉分布频率,%;n表示进汞压力测试总点数。
上述的测定方法中,步骤4)中,根据式(6)得到所述薄膜束缚水饱和度;
式(6)中,Swni表示第i个进汞区间薄膜束缚水饱和度,%;Swn:岩心薄膜束缚水饱和度,%,n表示进汞压力测试总点数。
本发明由于采取以上技术方案,其具有以下优点:
1、本发明可以根据实验压汞曲线数据,快速得到岩心薄膜束缚水饱和度,具有简便实用、可操作性强的特点。
2、本发明有效弥补了目前实验方法不能对薄膜束缚水饱和度进行定量测定的缺点。
3、基于本发明方法,可以进一步对薄膜束缚水饱和度和毛细管束缚水饱和度构成进行定量评价。
4、本发明薄膜束缚水饱和度的测定结果可应用于储层评价、产能预测、储量计算以及油气藏出水机理研究等众多方面,具有广阔的应用前景。
附图说明
图1为束缚水的两种形式示意图,其中,图1a表示毛细管束缚水,图1b表示薄膜束缚水。
图2为由压汞曲线确定孔喉分布的示意图。
图3为薄模束缚水同心毛细管的模型图。
图4为本发明方法的流程图。
图5为本发明实施例中的压汞曲线确定孔喉分布图。
具体实施方式
下述实施例中所使用的实验方法如无特殊说明,均为常规方法。
下述实施例中所用的材料、试剂等,如无特殊说明,均可从商业途径得到。
如图4所示,本发明提供的基于压汞曲线的岩心薄膜束缚水饱和度的测定方法包括以下步骤:
1、根据压汞曲线数据,确定岩心孔喉分布
由压汞曲线数据,确定岩心孔喉分布计算公式如下:
1.1孔喉半径计算
式中,Pi:第i个测试点进汞压力(以下简称第i个进汞压力),MPa;ri:第i个进汞压力对应的孔喉半径,μm;σ:汞表面张力,0.48mN/m;θ:汞润湿角,140°;n:进汞压力测试总点数。
1.2孔喉分布频率计算
压汞实验中两相邻测试点进汞压力(Pi,Pi+1)构成一个进汞区间i(简称区间)。其中,孔喉分布频率(fi)等于区间汞饱和度增量ΔSi,其计算公式为:
fi=ΔSi=Si+1-Si(i=1,2,……,n-1;)
式中,fi:第i个区间孔喉分布频率,%;ΔSi:第i个区间汞饱和度增量,%;Si+1:第i+1个进汞压力Pi+1对应的汞饱和度,%;Si:第i个进汞压力Pi对应的汞饱和度,%;n:进汞压力测试总点数。
2、确定薄膜束缚水厚度
式中,Pn:第n个进汞压力,MPa;rn:束缚水薄膜厚度,μm;σ:汞表面张力,0.48N/m;θ:汞润湿角,140°。
3、求取不同进汞区间薄膜束缚水饱和度
3.1计算不同进汞区间孔喉平均半径
式中,rci:第i个进汞区间孔喉平均半径,μm;ri:第i个进汞压力对应的孔喉半径,μm;ri+1:第i+1个进汞压力对应的孔喉半径,μm;n:进汞压力测试总点数。
3.2计算不同进汞区间薄膜束缚水饱和度
建立了薄模束缚水同心毛细管模型,如图3所示,设区间薄膜束缚水饱和度为Swni,其计算公式为:
式中,rci:第i个进汞区间孔喉平均半径,μm;rn:束缚水薄膜厚度,μm;Swni:第i个进汞区间薄膜束缚水饱和度,%;ΔSi:第i个进汞区间汞饱和度增量,%;n:进汞压力测试总点数。
4、确定岩心薄膜束缚水饱和度
式中,Swni:第i个进汞区间薄膜束缚水饱和度,%;Swn:岩心薄膜束缚水饱和度,%,n表示进汞压力测试总点数。
下面通过一个具体的实施例,进一步说明本发明的技术效果。
南梁油田位于鄂尔多斯盆地陕北斜坡西南部,地处甘肃省华池、庆阳县境内,面积约2600km2。主力层为三叠系上统延长组长6油组长63砂组,属于深湖-半深湖相重力流沉积,储层平均孔隙度9.09%、平均渗透率0.213×10-3μm2,属于特低孔-超低渗储层。以南梁油田长6油组山156井18#岩心(2060.1m)压汞曲线为实施例,该岩样压汞实验基础数据如表1所示,压汞实验测试数据如表2所示。
表1压汞实验基础数据表
表2压汞实验测试数据表
1)根据压汞曲线数据,确定岩心孔喉分布
由压汞曲线数据,确定岩心孔喉分布计算公式如下:
①孔喉半径计算
式中,Pi:第i个测试点进汞压力(以下简称第i个进汞压力),MPa;ri:第i个进汞压力对应的孔喉半径,μm;σ:汞表面张力,0.48mN/m;θ:汞润湿角,140°;n:进汞压力测试总点数,计算结果如表3所示。
②孔喉分布频率计算
压汞实验中两相邻测试点进汞压力(Pi,Pi+1)构成一个进汞区间i(简称区间)。其中,孔喉分布频率(fi)等于区间汞饱和度增量ΔSi,其计算公式为:fi=ΔSi=Si+1-Si(i=1,2,……,n-1;)
式中,fi:第i个区间孔喉分布频率,%;ΔSi:第i个区间汞饱和度增量,%;Si+1:第i+1个进汞压力Pi+1对应的汞饱和度,%;Si:第i个进汞压力Pi对应的汞饱和度,%;n:进汞压力测试总点数,计算结果如表3所示,由18#岩心压汞曲线确定孔喉分布如图5。
2)确定薄膜束缚水厚度
式中,Pn:第n个进汞压力,MPa;rn:束缚水薄膜厚度,μm;σ:汞表面张力,0.48N/m;θ:汞润湿角,140°,取n=29,计算结果如表3所示。
3)求取不同进汞区间薄膜束缚水饱和度
①计算不同进汞区间孔喉平均半径
式中,rci:第i个进汞区间孔喉平均半径,μm;ri:第i个进汞压力对应的孔喉半径,μm;ri+1:第i+1个进汞压力对应的孔喉半径,μm;n:进汞压力测试总点数,计算结果如表3所示。
②计算不同进汞区间薄膜束缚水饱和度
式中,rci:第i个进汞区间孔喉平均半径,μm;rn:束缚水薄膜厚度,μm;Swni:第i个进汞区间薄膜束缚水饱和度,%;ΔSi:第i个进汞区间汞饱和度增量,%;n:进汞压力测试总点数,计算结果如表3所示。
4)确定岩心薄膜束缚水饱和度
式中,Swni:第i个进汞区间薄膜束缚水饱和度,%;Swn:岩心薄膜束缚水饱和度,%,n:进汞压力测试总点数,计算结果如表3所示。
表3薄膜束缚水饱和度计算结果表
以上描述仅为本申请的实施例以及对所运用技术原理的说明。本领域技术人员应当理解,本申请中所涉及的发明范围,并不限于上述技术特征的特定组合而成的技术方案,同时也应涵盖在不脱离所述发明构思的情况下,由上述技术特征或其等同特征进行任意组合而形成的其它技术方案。例如上述特征与本申请中公开的(但不限于)具有类似功能的技术特征进行互相替换而形成的技术方案。
Claims (7)
1.一种岩心薄膜束缚水饱和度的测定方法,包括如下步骤:
1)对岩心进行压汞实验,得到不同进汞压力下的进汞饱和度,进而得到岩心的孔喉半径和孔喉分布频率;
2)根据所述进汞压力得到岩心的薄膜束缚水厚度;
3)根据所述孔喉半径得到岩心的不同进汞区间孔喉平均半径;
4)根据所述孔喉分布频率、所述薄膜束缚水厚度和所述不同进汞区间孔喉平均半径得到岩心的不同进汞区间薄膜束缚水饱和度,进而得到岩心的薄膜束缚水饱和度。
3.根据权利要求1或2所述的测定方法,其特征在于:步骤1)中,根据所述进汞饱和度,结合式(2)得到所述孔喉分布频率;
fi=ΔSi=Si+1-Si (2)
式(2)中,fi表示第i个区间孔喉分布频率,%;ΔSi表示第i个区间汞饱和度增量,%;Si+1表示第i+1个测试点进汞压力对应的汞饱和度,%;Si表示第i个测试点进汞压力对应的汞饱和度,%;i=1,2,……,n-1,n表示进汞压力测试总点数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910192711.6A CN111693425B (zh) | 2019-03-14 | 2019-03-14 | 一种基于压汞曲线的岩心薄膜束缚水饱和度的测定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910192711.6A CN111693425B (zh) | 2019-03-14 | 2019-03-14 | 一种基于压汞曲线的岩心薄膜束缚水饱和度的测定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111693425A CN111693425A (zh) | 2020-09-22 |
CN111693425B true CN111693425B (zh) | 2021-05-11 |
Family
ID=72475089
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910192711.6A Expired - Fee Related CN111693425B (zh) | 2019-03-14 | 2019-03-14 | 一种基于压汞曲线的岩心薄膜束缚水饱和度的测定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111693425B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114112841B (zh) * | 2021-11-19 | 2023-06-09 | 中国石油大学(华东) | 一种致密砂岩储层束缚水饱和度的计算方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2341372A1 (en) * | 2009-12-16 | 2011-07-06 | BP Exploration Operating Company Limited | Method for measuring rock wettability |
CN101799559B (zh) * | 2010-02-09 | 2011-12-07 | 康志勇 | 基于砂岩油藏油层电阻率的地层数据处理方法 |
CN103278436B (zh) * | 2013-02-01 | 2015-02-04 | 西安石油大学 | 特低渗透双重介质砂岩油藏微观孔隙结构的定量表征方法 |
CN104458525B (zh) * | 2013-09-13 | 2017-04-05 | 中国石油天然气股份有限公司 | 一种常规压汞曲线表征微观孔隙喉道特征的采集处理方法 |
CN103900942B (zh) * | 2014-01-10 | 2016-06-01 | 中国石油大学(华东) | 基于高压压汞分析的储层微观孔喉参数的连续表征方法 |
US10495774B2 (en) * | 2014-01-13 | 2019-12-03 | Schlumberger Technology Corporation | Method for estimating irreducible water saturation from mercury injection capillary pressure |
CN104200105A (zh) * | 2014-09-05 | 2014-12-10 | 中国石油大学(华东) | 确定致密砂岩气充注物性下限的方法 |
CN107944671B (zh) * | 2017-11-07 | 2021-11-02 | 中国石油天然气股份有限公司 | 一种确定地层束缚水饱和度的方法 |
-
2019
- 2019-03-14 CN CN201910192711.6A patent/CN111693425B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN111693425A (zh) | 2020-09-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110346258B (zh) | 一种致密岩石油相相对渗透率测定的方法 | |
CN107063967B (zh) | 一种基于微观孔隙通道的饱和/非饱和渗透系数预测方法 | |
CN107525746B (zh) | 一种表征致密油藏低速非达西渗流特征的方法及装置 | |
CN110598167B (zh) | 低渗透油藏油水相对渗透率实验数据的处理方法 | |
Muskat et al. | Flow of gas through porous materials | |
Wan et al. | Chloride content and pH value in the pore solution of concrete under carbonation | |
CN114720331B (zh) | 岩土体迂曲度幂律浆液时变性的柱半球渗透半径确定方法 | |
CN110516322B (zh) | 一种不同水力梯度下黏土饱和非线性渗透系数的预测方法 | |
CN111693425B (zh) | 一种基于压汞曲线的岩心薄膜束缚水饱和度的测定方法 | |
CN108444890B (zh) | 测试中、高渗岩心液测渗透率的非稳态滴定装置及方法 | |
Lu et al. | Constant flow method for concurrently measuring soil-water characteristic curve and hydraulic conductivity function | |
CN115618643B (zh) | 宾汉浆液时水效应岩土体迂曲度柱半球渗透半径确定方法 | |
CN106777515B (zh) | 一种基于岩心实验资料分析致密气井产能的方法 | |
CN108133086A (zh) | 一种应力敏感储层中产水压裂气井裂缝半长反演方法 | |
Yanjie et al. | Lower limit of tight oil flowing porosity: Application of high-pressure mercury intrusion in the fourth Member of Cretaceous Quantou Formation in southern Songliao Basin, NE China | |
CN109520891A (zh) | 过滤材料接触角的测定方法 | |
CN116383573B (zh) | 一种基于多区相变传质渗流耦合的凝析气产能评价方法 | |
CN109932296B (zh) | 一种定量表征贾敏效应动态变化的方法 | |
CN111208053A (zh) | 一种基于恒压压汞实验的储层微毛细管孔隙度确定方法 | |
CN114112841B (zh) | 一种致密砂岩储层束缚水饱和度的计算方法 | |
CN115083543A (zh) | 一种超大掺量矿物掺合料水泥石毛细吸水性能的预测方法 | |
CN115559715A (zh) | 一种超高压低渗透气藏储层产水评价方法 | |
CN207457017U (zh) | 一种用于测试多孔工程材料孔隙分布特征的试验装置 | |
CN109490162A (zh) | 考虑应力敏感条件下的致密油藏毛细管压力计算方法 | |
CN113945497B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210511 |
|
CF01 | Termination of patent right due to non-payment of annual fee |