CN105627965B - 基于真实岩心流动试验的储层液膜厚度的计算方法 - Google Patents

基于真实岩心流动试验的储层液膜厚度的计算方法 Download PDF

Info

Publication number
CN105627965B
CN105627965B CN201410638651.3A CN201410638651A CN105627965B CN 105627965 B CN105627965 B CN 105627965B CN 201410638651 A CN201410638651 A CN 201410638651A CN 105627965 B CN105627965 B CN 105627965B
Authority
CN
China
Prior art keywords
thickness
liquid film
core
permeability
capillary
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
Application number
CN201410638651.3A
Other languages
English (en)
Other versions
CN105627965A (zh
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 Petroleum and Chemical Corp
Sinopec Shengli Geological Scientific Reserch Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Shengli Geological Scientific Reserch Institute
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 Petroleum and Chemical Corp, Sinopec Shengli Geological Scientific Reserch Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201410638651.3A priority Critical patent/CN105627965B/zh
Publication of CN105627965A publication Critical patent/CN105627965A/zh
Application granted granted Critical
Publication of CN105627965B publication Critical patent/CN105627965B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本发明提供一种基于真实岩心流动试验的储层液膜厚度的计算方法,该方法包括:步骤1,获得真实岩心的有效渗流能力;步骤2,计算毛细管模型的理想渗透率;步骤3,建立当膜厚度为h时,毛细管模型有效渗透率的计算公式;步骤4,通过研究岩心的液膜厚度h与有效渗透率的关系,建立有效渗透率与液膜厚度的关系式;以及步骤5,计算出毛细管模型的有效渗透率,并根据有效渗透率与液膜厚度的关系式,计算出有效渗透率与液膜厚度的关系式。该方法能够定量表征液体进入岩心后,在岩石矿物表面吸附产生的液膜的厚薄,为深入剖析低渗砂岩油藏有效渗流能力的作用机理明确了切入点,为致密油藏的合理、有效开发奠定基础。

Description

基于真实岩心流动试验的储层液膜厚度的计算方法
技术领域
本发明涉及石油、天然气等地下流体在多孔介质内的渗流理论研究,特别是涉及到一种储层孔喉液膜厚度的间接计算方法。
背景技术
储层中的液体在与储层固相颗粒接触时,存在一种表面现象,有一层液体直接贴附在颗粒表面,这层液体即为液膜。液膜的存在降低了储层孔喉的有效半径。尤其是对于低渗油藏,液膜的厚度对储层评价、储量计算以及储层的渗流均有明显影响。为此我们发明了一种新的基于真实岩心流动试验的储层液膜厚度的计算方法,解决了以上技术问题。
发明内容
本发明的目的是提供一种能够定量计算液膜,为深入剖析低渗砂岩油藏有效渗流能力的作用机理明确了切入点的基于真实岩心流动试验的储层液膜厚度的计算方法。
本发明的目的可通过如下技术措施来实现:基于真实岩心流动试验的储层液膜厚度的计算方法,该基于真实岩心流动试验的储层液膜厚度的计算方法包括:步骤1,通过测试具体研究区块真实岩心的气体渗透率和液体渗透率,获得真实岩心的有效渗流能力;步骤2,通过真实岩心的毛管压力曲线,建立不等径毛管束模型,计算毛细管模型的理想渗透率;步骤3,建立当膜厚度为h时,毛细管模型有效渗透率的计算公式;步骤4,通过研究岩心的液膜厚度h与有效渗透率的关系,建立有效渗透率与液膜厚度的关系式;以及步骤5,利用真实岩心与毛管束模型流动等效,计算出毛细管模型的有效渗透率,并根据有效渗透率与液膜厚度的关系式,计算出有效渗透率与液膜厚度的关系式。
本发明的目的还可通过如下技术措施来实现:
在步骤1中,计算真实岩心的有效渗流能力的公式为:
PA岩心=KL/Ka (1)
其中:PA岩心为岩心的有效渗流能力;KL为岩心的液体渗透率,单位表示为10-3μm2;Ka为气体渗透率,单位表示为10-3μm2
在步骤2中,通过真实岩心的毛管压力曲线,把真实岩心假设为一组等长、不同直径的毛管束所组成,建立岩石-不等径毛细管模型,认为真实岩心的孔隙体积与毛管束体积相同,流体在单根毛管内的流动均遵循“泊稷叶”公式,流体在岩石内的流动遵循达西公式,假设真实岩心与毛管束模型在外观尺寸、流体性质、作用压差均相同,那么单根毛管内的流量与真实岩心的流量应该相等,毛细管模型与真实岩心流动等效。
在步骤2中,计算毛细管模型的理想渗透率的公式为:
其中:K理想为毛细管模型的理想渗透率,单位为10-3μm2;λ为常数;φ为孔隙度,用百分数表示;r为孔喉半径,单位为μm;s为饱和度,用百分数表示;sD为对应最大半径的汞饱和度,用百分数表示。
在步骤3中,建立的当膜厚度为h时,毛细管模型有效渗透率的计算公式:
其中:K有效为毛细管模型的有效渗透率,单位为10-3μm2;h为储层孔喉液膜厚度,单位为μm。
在步骤4中,根据毛细管模型有效渗透率的计算公式,给定一个液膜厚度h,就可以得到一个K有效值,拟合得到有效渗透率与液膜厚度的关系式:
h=aK有效 2+b K有效+c (5)
式中,a、b、c均为常数。
在步骤5中,计算毛细管模型的有效渗透率的公式为:
将计算出的毛细管模型的K有效值代入公式(5)中,求出储层孔喉液膜厚度。
本发明中的储层孔喉液膜厚度计算方法,涉及微观渗流实验分析,是研究石油天然气等地下流体在多孔介质内渗流规律的重要方法。在真实岩心有效渗流能力关系建立及基于真实岩心的毛管束模型的有效渗流能力建立的基础上,计算储层孔喉液膜的厚度。随着油田开发水平的不断提高,致密油已经成为最现实的待开发油气资源,搞清储层孔喉液膜厚度对于致密油的有效开发奠定了理论基础。该基于真实岩心流动试验的储层液膜厚度的计算方法是以液体在导管中作粘滞性流动的“泊稷叶”公式为基础,首先测试天然岩心的气体渗透率和岩心的液体渗透率;然后测试天然岩心的毛管压力曲线获取孔隙结构数据建立岩石-毛管束模型的理想渗透率和有效渗透率。该方法能够定量表征液体进入岩心后,在岩石矿物表面吸附产生的液膜的厚薄,为深入剖析低渗砂岩油藏有效渗流能力的作用机理明确了切入点,为致密油藏的合理、有效开发奠定基础。
附图说明
图1为本发明的基于真实岩心流动试验的储层液膜厚度的计算方法的一具体实施例的流程图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合所附图式,作详细说明如下。
如图1所示,图1为本发明的基于真实岩心流动试验的储层液膜厚度的计算方法的一具体实施例的流程图。
在步骤101,通过测试具体研究区块真实岩心的气体渗透率和液体渗透率,获得真实岩心的有效渗流能力:
PA岩心=KL/Ka (1)
其中:PA岩心为岩心的有效渗流能力;KL为岩心的液体渗透率,单位表示为10-3μm2;Ka为气体渗透率,单位表示为10-3μm2。流程进入到步骤102。
在步骤102,通过真实岩心的毛管压力曲线,把真实岩心假设为一组等长、不同直径的毛管束所组成,建立岩石-不等径毛细管模型,认为真实岩心的孔隙体积与毛管束体积相同,流体在单根毛管内的流动均遵循“泊稷叶”公式,流体在岩石内的流动遵循达西公式,假设真实岩心与毛管束模型在外观尺寸、流体性质、作用压差均相同,那么单根毛管内的流量与真实岩心的流量应该相等,毛细管模型与真实岩心流动等效。
通过测试具体研究区块真实岩心的毛管压力曲线,得到岩心的孔喉大小r及分布s,计算毛细管模型的理想渗透率:
其中:K理想为毛细管模型的理想渗透率,单位为10-3μm2;λ为常数;φ为孔隙度,用百分数表示;r为孔喉半径,单位为μm;s为饱和度,用百分数表示;sD为对应最大半径的汞饱和度,用百分数表示。流程进入到步骤103。
在步骤103,由于液体进入储层,在固-液界面张力的作用下,产生液膜,减小了有效渗流半径,因此根据理想渗透率的计算方法,可以建立当膜厚度为h时,毛细管模型有效渗透率的计算公式:
其中:K有效为毛细管模型的有效渗透率,单位为10-3μm2;h为储层孔喉液膜厚度,单位为μm。流程进入到步骤104。
在步骤104,通过研究岩心的液膜厚度h与有效渗透率K有效的拟合函数得:
h=f(K有效) (4)
在一实施例中,给定一个液膜厚度h,就可以得到一个K有效值,这样就可以建立有效渗透率与液膜厚度的关系式:h=aK有效 2+b K有效+c(5),式中,a、b、c均为常数。流程进入到步骤105。
在步骤105,利用真实岩心与毛管束模型流动等效,计算出毛细管模型的有效渗透率:
将上步计算出的毛细管模型的K有效值代入公式(5)中,即可求出对应实际液体渗透率的岩心孔喉中液膜厚度h。流程结束。
在应用本发明的一具体实施例中,测试真实岩心的气体渗透率Ka=1.47×10-3μm2,孔隙度φ=20.4%、液体渗透率KL=0.881×10-3μm2,得到真实岩心的有效渗流能力PA岩心=0.599。
测试真实岩心的毛管压力曲线,得到岩心的孔喉分布情况(表1),利用进行计算,得到岩心的理想渗透率K理想=9.9323×10-3μm2
表1 毛管压力试验数据
利用表1中的数据,根据有效渗透率的计算公式,可以算出固定膜厚对应的有效渗透率值(表2)。利用表2中的数据进行数据拟合,得到岩心的液膜厚度h与有效渗透率K有效的拟合函数h=0.0138K有效 2-0.2828K有效+1.502(3)。
表2 膜厚与有效渗透率的对应表
膜厚度μm 有效渗透率10-3μm2 膜厚度μm 有效渗透率10-3μm2
0.05 9.1416 0.1 8.4286
0.2 7.1847 0.3 6.1209
0.4 5.2271 0.5 4.4865
0.6 3.8605 0.7 3.3176
0.8 2.8517 0.9 2.4315
1.0 2.0530
利用真实岩心与毛管束模型流动等效,即计算出岩心-毛细管模型的有效渗透率K有效=5.9494×10-3μm2。将毛细管模型的有效渗透率K有效值代入公式(3)中,即可计算储层孔喉液膜厚度h=0.31μm。

Claims (7)

1.基于真实岩心流动试验的储层液膜厚度的计算方法,其特征在于,该基于真实岩心流动试验的储层液膜厚度的计算方法包括:
步骤1,通过测试具体研究区块真实岩心的气体渗透率和液体渗透率,获得真实岩心的有效渗流能力;
步骤2,通过真实岩心的毛管压力曲线,建立不等径毛细管模型,计算毛细管模型的理想渗透率;
步骤3,建立当液膜厚度为h时,毛细管模型有效渗透率的计算公式;
步骤4,通过研究岩心的液膜厚度h与有效渗透率的关系,建立有效渗透率与液膜厚度的关系式;以及
步骤5,利用真实岩心与毛细管模型流动等效,计算出毛细管模型的有效渗透率,并根据有效渗透率与液膜厚度的关系式,计算出储层孔喉液膜厚度。
2.根据权利要求1所述的基于真实岩心流动试验的储层液膜厚度的计算方法,其特征在于,在步骤1中,计算真实岩心的有效渗流能力的公式为:
PA岩心=KL/Ka (1)
其中:PA岩心为岩心的有效渗流能力;KL为岩心的液体渗透率,单位表示为10-3μm2;Ka为气体渗透率,单位表示为10-3μm2
3.根据权利要求2所述的基于真实岩心流动试验的储层液膜厚度的计算方法,其特征在于,在步骤2中,通过真实岩心的毛管压力曲线,把真实岩心假设为一组等长、不同直径的毛管束所组成,建立岩石-不等径毛细管模型,认为真实岩心的孔隙体积与毛管束体积相同,流体在单根毛管内的流动均遵循“泊稷叶”公式,流体在岩石内的流动遵循达西公式,假设真实岩心与毛细管模型在外观尺寸、流体性质、作用压差均相同,那么单根毛管内的流量与真实岩心的流量应该相等,毛细管模型与真实岩心流动等效。
4.根据权利要求3所述的基于真实岩心流动试验的储层液膜厚度的计算方法,其特征在于,在步骤2中,计算毛细管模型的理想渗透率的公式为:
其中:K理想为毛细管模型的理想渗透率,单位为10-3μm2;λ为常数;φ为孔隙度,用百分数表示;r为孔喉半径,单位为μm;s为饱和度,用百分数表示;sD为对应最大半径的汞饱和度,用百分数表示。
5.根据权利要求4所述的基于真实岩心流动试验的储层液膜厚度的计算方法,其特征在于,在步骤3中,建立的当膜厚度为h时,毛细管模型有效渗透率的计算公式:
其中:K有效为毛细管模型的有效渗透率,单位为10-3μm2;h为储层孔喉液膜厚度,单位为μm。
6.根据权利要求5所述的基于真实岩心流动试验的储层液膜厚度的计算方法,其特征在于,在步骤4中,根据毛细管模型有效渗透率的计算公式,给定一个液膜厚度h,就可以得到一个K有效值,拟合得到有效渗透率与液膜厚度的关系式:
h=a K有效 2+b K有效+c (5)
式中,a、b、c均为常数。
7.根据权利要求6所述的基于真实岩心流动试验的储层液膜厚度的计算方法,其特征在于,在步骤5中,计算毛细管模型的有效渗透率的公式为:
将计算出的毛细管模型的K有效值代入公式(5)中,求出储层孔喉液膜厚度。
CN201410638651.3A 2014-11-07 2014-11-07 基于真实岩心流动试验的储层液膜厚度的计算方法 Expired - Fee Related CN105627965B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410638651.3A CN105627965B (zh) 2014-11-07 2014-11-07 基于真实岩心流动试验的储层液膜厚度的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410638651.3A CN105627965B (zh) 2014-11-07 2014-11-07 基于真实岩心流动试验的储层液膜厚度的计算方法

Publications (2)

Publication Number Publication Date
CN105627965A CN105627965A (zh) 2016-06-01
CN105627965B true CN105627965B (zh) 2018-06-26

Family

ID=56043118

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410638651.3A Expired - Fee Related CN105627965B (zh) 2014-11-07 2014-11-07 基于真实岩心流动试验的储层液膜厚度的计算方法

Country Status (1)

Country Link
CN (1) CN105627965B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106383133B (zh) * 2016-10-21 2019-03-05 大连理工大学 一种基于x射线ct成像技术的co2-盐水-岩芯体系毛细管压力测量方法
CN110286067B (zh) * 2019-07-04 2021-11-30 中国石油天然气股份有限公司 一种定量表征多孔介质内等效边界层厚度的方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101592584A (zh) * 2008-05-16 2009-12-02 道达尔公司 一种估量地质构造的物理参数的方法
CN102455277A (zh) * 2010-10-22 2012-05-16 中国石油化工股份有限公司 一种高压下岩石气测渗透率的装置及方法
CN103759680A (zh) * 2013-12-31 2014-04-30 中国石油天然气股份有限公司 一种致密储层微纳米孔喉中油膜赋存厚度的测量方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005022924A (ja) * 2003-07-02 2005-01-27 Japan Fine Ceramics Center 細孔基材及びその製造方法並びにガス分離材用細孔基材

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101592584A (zh) * 2008-05-16 2009-12-02 道达尔公司 一种估量地质构造的物理参数的方法
CN102455277A (zh) * 2010-10-22 2012-05-16 中国石油化工股份有限公司 一种高压下岩石气测渗透率的装置及方法
CN103759680A (zh) * 2013-12-31 2014-04-30 中国石油天然气股份有限公司 一种致密储层微纳米孔喉中油膜赋存厚度的测量方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
岩石全应力-应变过程渗透变化规律分析;江东辉等;《金属矿山》;20120229(第2期);第22-24、44页 *

Also Published As

Publication number Publication date
CN105627965A (zh) 2016-06-01

Similar Documents

Publication Publication Date Title
CN105574320B (zh) 低渗砂岩储层有效渗流能力的评价方法
Wang Performance of multiple fractured horizontal wells in shale gas reservoirs with consideration of multiple mechanisms
Zhao et al. “Triple porosity” modeling of transient well test and rate decline analysis for multi-fractured horizontal well in shale gas reservoirs
Gao et al. Pore structure characterization, permeability evaluation and enhanced gas recovery techniques of tight gas sandstones
Clarkson et al. Production analysis of tight-gas and shale-gas reservoirs using the dynamic-slippage concept
CN110598167B (zh) 低渗透油藏油水相对渗透率实验数据的处理方法
CN107290259B (zh) 低渗砂岩储层有效渗流孔喉半径的计算方法
Alharthy et al. Physics and modeling of gas flow in shale reservoirs
Gao et al. Apparent permeability and gas flow behavior in carboniferous shale from the Qaidam Basin, China: an experimental study
Zhang et al. Wettability alteration to intermediate gas-wetting in low-permeability gas-condensate reservoirs
Cheng et al. A fractal irreducible water saturation model for capillary tubes and its application in tight gas reservoir
Ye et al. A unified method to evaluate shale gas flow behaviours in different flow regions
Chen et al. Gas slippage in anisotropically-stressed shale: An experimental study
CN105627965B (zh) 基于真实岩心流动试验的储层液膜厚度的计算方法
Sakhaee-Pour et al. Effective flow properties for cells containing fractures of arbitrary geometry
Guo et al. Pressure transient and rate decline analysis for hydraulic fractured vertical wells with finite conductivity in shale gas reservoirs
Su et al. A new method for continental shale oil enrichment evaluation
Meng et al. Study of water Huff-n-Puff in low-permeability oil reservoirs with horizontal fractures: a case study of Chang 6 reservoir in Yanchang, China
CN106769745B (zh) 单相流体通过低渗透岩心的非线性渗流特征的表征方法
Shrey et al. Modifying proppant surface with nano-roughness coating to enhance fracture conductivity
Zhang et al. A novel quantitative petrophysical model for the stress sensitivity of tight sandstones
CN108920824B (zh) 一种窄河道储层的产量递减分析方法
Huang et al. A nonlinear seepage model of gas and water transport in multi-scale shale gas reservoirs based on dynamic permeability
Zhichun et al. Numerical well test analysis for polymer flooding considering the non-Newtonian behavior
Hu et al. Acidizing flowback optimization for tight sandstone gas reservoirs

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into 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: 20180626

Termination date: 20191107

CF01 Termination of patent right due to non-payment of annual fee