CN105627965A - 基于真实岩心流动试验的储层液膜厚度的计算方法 - Google Patents
基于真实岩心流动试验的储层液膜厚度的计算方法 Download PDFInfo
- Publication number
- CN105627965A CN105627965A CN201410638651.3A CN201410638651A CN105627965A CN 105627965 A CN105627965 A CN 105627965A CN 201410638651 A CN201410638651 A CN 201410638651A CN 105627965 A CN105627965 A CN 105627965A
- Authority
- CN
- China
- Prior art keywords
- liquid film
- thickness
- core
- capillary
- permeability
- 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
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+bK有效+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+bK有效+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=aK有效 2+bK有效+c(5)
式中,a、b、c均为常数。
7.根据权利要求6所述的基于真实岩心流动试验的储层液膜厚度的计算方法,其特征在于,在步骤5中,计算毛细管模型的有效渗透率的公式为:
将计算出的毛细管模型的K有效值代入公式(5)中,求出储层孔喉液膜厚度。
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 true CN105627965A (zh) | 2016-06-01 |
CN105627965B 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) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106383133A (zh) * | 2016-10-21 | 2017-02-08 | 大连理工大学 | 一种基于x射线ct成像技术的co2‑盐水‑岩芯体系毛细管压力测量方法 |
CN110286067A (zh) * | 2019-07-04 | 2019-09-27 | 中国石油天然气股份有限公司 | 一种定量表征多孔介质内等效边界层厚度的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005022924A (ja) * | 2003-07-02 | 2005-01-27 | Japan Fine Ceramics Center | 細孔基材及びその製造方法並びにガス分離材用細孔基材 |
CN101592584A (zh) * | 2008-05-16 | 2009-12-02 | 道达尔公司 | 一种估量地质构造的物理参数的方法 |
CN102455277A (zh) * | 2010-10-22 | 2012-05-16 | 中国石油化工股份有限公司 | 一种高压下岩石气测渗透率的装置及方法 |
CN103759680A (zh) * | 2013-12-31 | 2014-04-30 | 中国石油天然气股份有限公司 | 一种致密储层微纳米孔喉中油膜赋存厚度的测量方法 |
-
2014
- 2014-11-07 CN CN201410638651.3A patent/CN105627965B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005022924A (ja) * | 2003-07-02 | 2005-01-27 | Japan Fine Ceramics Center | 細孔基材及びその製造方法並びにガス分離材用細孔基材 |
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)
Title |
---|
江东辉等: "岩石全应力-应变过程渗透变化规律分析", 《金属矿山》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106383133A (zh) * | 2016-10-21 | 2017-02-08 | 大连理工大学 | 一种基于x射线ct成像技术的co2‑盐水‑岩芯体系毛细管压力测量方法 |
CN106383133B (zh) * | 2016-10-21 | 2019-03-05 | 大连理工大学 | 一种基于x射线ct成像技术的co2-盐水-岩芯体系毛细管压力测量方法 |
CN110286067A (zh) * | 2019-07-04 | 2019-09-27 | 中国石油天然气股份有限公司 | 一种定量表征多孔介质内等效边界层厚度的方法 |
CN110286067B (zh) * | 2019-07-04 | 2021-11-30 | 中国石油天然气股份有限公司 | 一种定量表征多孔介质内等效边界层厚度的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105627965B (zh) | 2018-06-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105574320B (zh) | 低渗砂岩储层有效渗流能力的评价方法 | |
CN106599449B (zh) | 一种溶洞体积计算的试井解释方法 | |
CN110598167A (zh) | 低渗透油藏油水相对渗透率实验数据的处理方法 | |
MX2014011709A (es) | Un metodo eficiente para seleccionar volumen elemental representativo en representaciones digitales de medios porosos. | |
Meng et al. | An improved productivity model in coal reservoir and its application during coalbed methane production | |
CN107290259B (zh) | 低渗砂岩储层有效渗流孔喉半径的计算方法 | |
CN107526891B (zh) | 一种聚合物驱大孔道油藏试井分析方法 | |
You et al. | Pore structure and limit pressure of gas slippage effect in tight sandstone | |
CN109339774A (zh) | 基于边界层的低渗透变形介质非线性渗流数值模拟方法 | |
CN106897531A (zh) | 一种低渗透石灰岩储层渗透率的定量评价方法 | |
Chen | Production from a fractured well with finite fracture conductivity in a closed reservoir: an exact analytical solution for pseudosteady-state flow | |
CN105627965A (zh) | 基于真实岩心流动试验的储层液膜厚度的计算方法 | |
CN104407397B (zh) | 一种致密油气储层物性下限的确定方法 | |
CN111577264A (zh) | 裂缝孔隙型油藏水平井产能预测方法及装置 | |
CN107843531B (zh) | 低渗透低粘油油层孔隙结构的核磁共振参数表征方法 | |
CN106769745B (zh) | 单相流体通过低渗透岩心的非线性渗流特征的表征方法 | |
CN108920824B (zh) | 一种窄河道储层的产量递减分析方法 | |
Wu et al. | Pressure analysis for volume fracturing well in tight reservoirs considering dynamic threshold pressure gradient and stress sensitivity | |
Hu et al. | Acidizing flowback optimization for tight sandstone gas reservoirs | |
Zhichun et al. | Numerical well test analysis for polymer flooding considering the non-Newtonian behavior | |
CN109490162A (zh) | 考虑应力敏感条件下的致密油藏毛细管压力计算方法 | |
CN106446489B (zh) | 连续型致密碎屑岩油藏储层物性上限的计算方法 | |
Jia et al. | Dynamic Effect of capillary pressure in tight gas reservoir | |
He et al. | Waterflood direction and front characterization with multiple methods: a case study in Changqing oilfield, China | |
CN106769688B (zh) | 一种预测地质条件下致密砂岩中天然气扩散系数的方法 |
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 |