CN107976390A - 一种基于覆压孔渗实验确定平均孔喉比的计算方法 - Google Patents

一种基于覆压孔渗实验确定平均孔喉比的计算方法 Download PDF

Info

Publication number
CN107976390A
CN107976390A CN201710954227.3A CN201710954227A CN107976390A CN 107976390 A CN107976390 A CN 107976390A CN 201710954227 A CN201710954227 A CN 201710954227A CN 107976390 A CN107976390 A CN 107976390A
Authority
CN
China
Prior art keywords
mrow
average pore
throat ratio
pressure
msub
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
CN201710954227.3A
Other languages
English (en)
Other versions
CN107976390B (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 National Petroleum Corp
China Petroleum Logging Co Ltd
Original Assignee
China National Petroleum Corp
China Petroleum Logging Co Ltd
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 National Petroleum Corp, China Petroleum Logging Co Ltd filed Critical China National Petroleum Corp
Priority to CN201710954227.3A priority Critical patent/CN107976390B/zh
Publication of CN107976390A publication Critical patent/CN107976390A/zh
Application granted granted Critical
Publication of CN107976390B publication Critical patent/CN107976390B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/082Investigating permeability by forcing a fluid through a sample
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/088Investigating 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)
  • Geophysics And Detection Of Objects (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种基于覆压孔渗实验确定平均孔喉比的计算方法,包括:建立覆压孔渗数据与平均孔喉比之间的理论关系;选取目的层岩心开展不同压力下的覆压孔渗实验,通过实验测量得到不同压力下的孔隙度、渗透率数值;利用覆压孔渗与平均孔喉比之间的数学关系与测量覆压孔渗数据,构建平均孔喉比的数学方程,通过迭代法求解平均孔喉比;将三种不同压力下孔渗数据求取的平均孔喉比进行平均,其平均值为该岩心的平均孔喉比。本发明在实际应用中,实现了在无压汞、核磁实验的情况下,通过覆压孔渗实验确定岩心孔喉比,实验成本低、快速方便。

Description

一种基于覆压孔渗实验确定平均孔喉比的计算方法
技术领域
本发明涉及石油勘探中的岩石物理实验数据深化应用技术领域,涉及一种基于覆压实验确定平均孔吼比的计算方法。
背景技术
岩石的孔隙空间可以划分为孔隙体与喉道,喉道为连通孔隙的狭小区域,控制岩石的渗流能力,孔隙体为孔隙空间较宽区域,反应岩石存储能力。因此,评价孔隙及喉道大小、连通性等孔隙结构参数是指导致密油气藏勘探开发、分析渗流与成藏机理的关键因素。目前,孔喉结构的评价手段主要有压汞、核磁共振、CT扫描等。压汞法是通过进汞压力涨落分别测量孔隙与喉道体积及数量,进而确定孔隙与喉道的大小分布及孔喉比等参数;核磁共振是通过一定刻度转换为伪毛管压力曲线后,应用与压汞法相似的分析法确定孔喉比等参数,或者通过T2几何均值与实验孔喉参数建立经验关系;CT实验是通过对岩样进行扫描切片成像,得到岩样孔喉三维立体图,直观显示孔喉的形态特征,通过三维建模最终计算得到孔吼比等参数。
以上方法中压汞实验、核磁共振实验需要的周期长而且实验费用高,确定的孔吼比是测量间接物理量压力的变化,通过一定的数学关系转换过来,精度受到一定影响;而CT扫描得到三维立体成图后,需要进行岩心的三维模型重构,数据运算量大,计算复杂,同时该实验费用较高。
发明内容
为解决现有技术中存在的上述缺陷,本发明的目的在于提供一种基于覆压孔渗实验确定平均孔喉比的计算方法,能够实现在无压汞、CT实验的情况下,通过不同压力下的覆压孔渗实验快速准确地得到岩样的平均孔喉比。
本发明是通过下述技术方案来实现的。
所述方法按照以下步骤进行操作:
步骤101,建立覆压孔渗数据与平均孔喉比之间的理论关系;
步骤102,选取目的层岩心开展不同压力下的覆压孔渗实验,通过实验测量得到不同压力下的孔隙度、渗透率数值;
步骤103,利用覆压孔渗与平均孔喉比之间的数学关系与测量覆压孔渗数据,构建平均孔喉比的数学方程,通过迭代法求解平均孔喉比;
步骤104,将三种不同压力下孔渗数据求取的平均孔喉比进行平均,其平均值为该岩心的平均孔喉比。
进一步地,所述步骤101中,依据串联毛管模型,推导覆压孔渗数据与平均孔喉比之间的计算公式如下:
式中,相关参数由下式得到:
式中:δ为平均孔喉比;A1为覆压下孔隙截面积变化率;A2为覆压下喉道截面积变化率;C为覆压前后孔隙度比值;Por为覆压为零时测量的孔隙度;Por′为一定覆压下测量的孔隙度;K为覆压为零时测量的渗透率;K′为一定覆压下测量的渗透率。
进一步地,所述步骤102中,需要进行岩心不同压力下的覆压孔渗实验,通过实验测量得到不同压力下的孔隙度、渗透率数值。
具体地,所述岩心覆压实验按照《岩心分析方法SY/T 5336-2006》标准规定的流程进行。
进一步地,所述步骤103中,利用覆压孔渗与平均孔喉比之间的数学关系与测量覆压孔渗数据,构建平均孔喉比的数学方程:
结合A1取值范围,通过迭代法求解F(δ)=0时的根δ,确定岩心平均孔喉比。
进一步地,所述步骤104中,求取岩心最终孔喉比的公式如下:
式中:为最终确定的岩心孔喉比;δ1为P1压力下实验数据确定的平均孔喉比;δ2为P2压力下实验数据确定的平均孔喉比;δ3为P3压力下实验数据确定的平均孔喉比。
进一步地,所述三种不同压力P1、P2和P3应满足逐渐递增,最小压力递增不小于5Mpa,本方法选取的实验压力值为5Mpa、15Mpa和30Mpa。
本发明的有益效果在于,实现了在无压汞与CT扫描实验的情况下,通过不用压力下的覆压孔渗实验数据,可以计算得到岩心的平均孔吼比,比以往的实验确定方法更加快速准确,而且成本较低。
附图说明
图1为本发明提供的一种基于覆压实验确定平均孔喉比方法的流程图;
图2(a)、2(b)分别为串联毛管模型的孔喉示意图;
图3为本发明专利确定的孔吼比与压汞实验确定孔吼比对比图。
具体实施方式
下面结合附图和实施例对发明作进一步的详细说明,但并不作为对发明做任何限制的依据。
如图1所示,本发明实施例提供的一种基于覆压孔渗实验确定平均孔喉比的计算方法,包括如下步骤:
步骤101:建立覆压孔渗数据与平均孔喉比之间的理论关系;
如图2(a)和图2(b)所示,依据串联毛管模型,推导覆压孔渗数据与平均孔喉比之间的计算公式如下:
式中,相关参数由下式得到:
式中:δ为平均孔喉比;A1为覆压下孔隙截面积变化率;A2为覆压下喉道截面积变化率;C为覆压前后孔隙度比值;Por为覆压为零时测量的孔隙度;Por′为一定覆压下测量的孔隙度;K为覆压为零时测量的渗透率;K′为一定覆压下测量的渗透率。
步骤102:选取目的层岩心开展不同压力下的覆压孔渗实验,通过实验测量得到不同压力下的孔隙度、渗透率数值;在进行该步骤操作之前,还需进行岩心覆压孔渗实验,并获取不同压力下岩心的孔隙度、渗透率测量值。选择某一目标研究区块,选取不同物性的岩心按照《岩心分析方法SY/T 5336-2006》标准规定,开展岩心不同压力下的覆压孔渗实验,并获取不同压力下的孔隙度、渗透率测量值,为后续的计算准备数据。
步骤103:利用覆压孔渗与平均孔喉比之间的数学关系与测量覆压孔渗数据,构建平均孔喉比的数学方程。
下面,通过对本实施例的具体实施情况做进一步详细说明,以支持本发明所要解决的技术问题。
利用平均孔喉比计算模型与覆压孔渗测量数据,确定计算模型中的相关参数,如表1所示:
表1
构建平均孔喉比数学方程为:
结合A1取值范围,通过迭代法求解F(δ)=0时的根δ,确定岩心平均孔喉比。
步骤104,通过三种不同压力下5Mpa、15Mpa和30Mpa的覆压孔渗数据分别求取的平均孔喉比,求取均值作为该岩心的平均孔喉比:
式中:为最终确定的岩心孔喉比;δ1为P1压力下实验数据确定的平均孔喉比;δ2为P2压力下实验数据确定的平均孔喉比;δ3为P3压力下实验数据确定的平均孔喉比,计算结果如表2所示。
表2
图3为利用本发明方法确定的平均孔喉比与压汞实验确定平均孔喉比对比图,从图中可以看出,本发明方法确定的平均孔喉比与压汞实验结果一致性较好,计算结果准确可靠。
最后所应说明的是,以上具体实施方式仅用以说明本发明的技术方案而非限制,尽管参照实例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (7)

1.一种基于覆压孔渗实验确定平均孔喉比的计算方法,其特征在于,包括步骤:
步骤101,建立覆压孔渗数据与平均孔喉比之间的理论关系;
步骤102,选取目的层岩心开展不同压力下的覆压孔渗实验,通过实验测量得到不同压力下的孔隙度、渗透率数值;
步骤103,利用覆压孔渗与平均孔喉比之间的数学关系与测量覆压孔渗数据,构建平均孔喉比的数学方程,通过迭代法求解平均孔喉比;
步骤104,将三种不同压力下孔渗数据求取的平均孔喉比进行平均,其平均值为该岩心的平均孔喉比。
2.如权利要求1所述的一种基于覆压孔渗实验确定平均孔喉比的计算方法,其特征在于,所述步骤101中,依据串联毛管模型,推导覆压孔渗数据与平均孔喉比之间的计算公式如下:
<mrow> <msup> <mi>&amp;delta;</mi> <mn>2</mn> </msup> <mo>=</mo> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <msub> <mi>A</mi> <mn>2</mn> </msub> <mo>&amp;CenterDot;</mo> <mi>C</mi> </mrow> <mrow> <msub> <mi>A</mi> <mn>1</mn> </msub> <mo>&amp;CenterDot;</mo> <mi>C</mi> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
式中,相关参数由下式得到:
<mrow> <mi>C</mi> <mo>=</mo> <mfrac> <mrow> <mi>P</mi> <mi>o</mi> <mi>r</mi> </mrow> <mrow> <msup> <mi>Por</mi> <mo>&amp;prime;</mo> </msup> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>A</mi> <mn>2</mn> </msub> <mo>=</mo> <mfrac> <msup> <mi>K</mi> <mo>&amp;prime;</mo> </msup> <mi>K</mi> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mfrac> <mn>1</mn> <mi>C</mi> </mfrac> <mo>&lt;</mo> <msub> <mi>A</mi> <mn>1</mn> </msub> <mo>&lt;</mo> <mn>1</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
式中:δ为平均孔喉比;A1为覆压下孔隙截面积变化率;A2为覆压下喉道截面积变化率;C为覆压前后孔隙度比值;Por为覆压为零时测量的孔隙度;Por′为一定覆压下测量的孔隙度;K为覆压为零时测量的渗透率;K′为一定覆压下测量的渗透率。
3.如权利要求1所述的一种基于覆压孔渗实验确定平均孔喉比的计算方法,其特征在于,在进行步骤102操作之前,还需进行岩心覆压孔渗实验,并获取不同压力下岩心的孔隙度、渗透率测量值。
4.如权利要求3所述的一种基于覆压孔渗实验确定平均孔喉比的计算方法,其特征在于,所述岩心覆压实验按照《岩心分析方法SY/T 5336-2006》标准规定的流程进行。
5.根据权利要求1所述的一种基于覆压孔渗实验确定平均孔喉比的计算方法,其特征在于,所述步骤103中,利用覆压孔渗与平均孔喉比之间的数学关系与测量覆压孔渗数据,构建平均孔喉比的数学方程如下:
<mrow> <mi>F</mi> <mrow> <mo>(</mo> <mi>&amp;delta;</mi> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mi>&amp;delta;</mi> <mn>2</mn> </msup> <mo>-</mo> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <msub> <mi>A</mi> <mn>2</mn> </msub> <mo>&amp;CenterDot;</mo> <mi>C</mi> </mrow> <mrow> <msub> <mi>A</mi> <mn>1</mn> </msub> <mo>&amp;CenterDot;</mo> <mi>C</mi> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
结合A1取值范围,通过迭代法求解F(δ)=0时的根δ,确定岩心平均孔喉比。
6.根据权利要求1所述的一种基于覆压孔渗实验确定平均孔喉比的计算方法,其特征在于,所述步骤104中,通过同一块岩心不同压力下测量的覆压孔渗数据,求取三种不同压力P1、P2和P3下的平均孔喉比,进行平均确定为该岩心的平均孔喉比:
<mrow> <mover> <mi>&amp;delta;</mi> <mo>&amp;OverBar;</mo> </mover> <mo>=</mo> <mfrac> <mrow> <msub> <mi>&amp;delta;</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>&amp;delta;</mi> <mn>2</mn> </msub> <mo>+</mo> <msub> <mi>&amp;delta;</mi> <mn>3</mn> </msub> </mrow> <mn>3</mn> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
式中:为最终确定的岩心孔喉比;δ1为P1压力下实验数据确定的平均孔喉比;δ2为P2压力下实验数据确定的平均孔喉比;δ3为P3压力下实验数据确定的平均孔喉比。
7.根据权利要求6所述的一种基于覆压孔渗实验确定平均孔喉比的计算方法,其特征在于,所述三种不同压力P1、P2和P3应满足逐渐递增,最小压力递增不小于5Mpa。
CN201710954227.3A 2017-10-13 2017-10-13 一种基于覆压孔渗实验确定平均孔喉比的计算方法 Active CN107976390B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710954227.3A CN107976390B (zh) 2017-10-13 2017-10-13 一种基于覆压孔渗实验确定平均孔喉比的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710954227.3A CN107976390B (zh) 2017-10-13 2017-10-13 一种基于覆压孔渗实验确定平均孔喉比的计算方法

Publications (2)

Publication Number Publication Date
CN107976390A true CN107976390A (zh) 2018-05-01
CN107976390B CN107976390B (zh) 2020-05-08

Family

ID=62012455

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710954227.3A Active CN107976390B (zh) 2017-10-13 2017-10-13 一种基于覆压孔渗实验确定平均孔喉比的计算方法

Country Status (1)

Country Link
CN (1) CN107976390B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4783751A (en) * 1983-08-17 1988-11-08 University Of South Carolina Analysis of pore complexes
US20120197526A1 (en) * 2011-01-27 2012-08-02 Instituto Mexicano Del Petroleo Procedure for the determination of effective and total porosity of carbonated sedimentary rocks, and morphology characterization of their micro and nanopores
US20130042670A1 (en) * 2010-04-30 2013-02-21 National University Of Singapore Determination of pore size in porous materials by evaporative mass loss
CN104101562A (zh) * 2013-04-15 2014-10-15 中国石油天然气集团公司 一种确定储层渗透率的方法
CN105466830A (zh) * 2015-11-19 2016-04-06 中国石油天然气股份有限公司 储层砂岩孔隙喉道尺寸分布识别方法
CN105844011A (zh) * 2016-03-22 2016-08-10 中国石油天然气集团公司 一种基于毛管模型的渗透率计算方法
CN106198345A (zh) * 2016-07-05 2016-12-07 西北大学 一种储层评价方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4783751A (en) * 1983-08-17 1988-11-08 University Of South Carolina Analysis of pore complexes
US20130042670A1 (en) * 2010-04-30 2013-02-21 National University Of Singapore Determination of pore size in porous materials by evaporative mass loss
US20120197526A1 (en) * 2011-01-27 2012-08-02 Instituto Mexicano Del Petroleo Procedure for the determination of effective and total porosity of carbonated sedimentary rocks, and morphology characterization of their micro and nanopores
CN104101562A (zh) * 2013-04-15 2014-10-15 中国石油天然气集团公司 一种确定储层渗透率的方法
CN105466830A (zh) * 2015-11-19 2016-04-06 中国石油天然气股份有限公司 储层砂岩孔隙喉道尺寸分布识别方法
CN105844011A (zh) * 2016-03-22 2016-08-10 中国石油天然气集团公司 一种基于毛管模型的渗透率计算方法
CN106198345A (zh) * 2016-07-05 2016-12-07 西北大学 一种储层评价方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
成志刚等: "基于储层孔喉特征参数计算致密砂岩渗透率", 《测井技术》 *
魏广振: "岩石孔喉平均半径SEM测定法", 《ANALYTICAL INSTRUMENTATION》 *

Also Published As

Publication number Publication date
CN107976390B (zh) 2020-05-08

Similar Documents

Publication Publication Date Title
Chen et al. Experimental study on the effect of fracture geometric characteristics on the permeability in deformable rough-walled fractures
Schelle et al. Analysis of the Agreement of Soil Hydraulic Properties Obtained from Multistep-Outflow and Evaporation MethodsAll rights reserved. No part of this periodical may be reproduced or transmitted in any form or by any means, electronic or mechanical, including photocopying, recording, or any information storage and retrieval system, without permission in writing from the publisher.
CN107340183B (zh) 结构性软土次固结系数描述方法
CN106018740B (zh) 孔压静力触探标定罐系统
CN109145473A (zh) 一种竖向排水板固结饱和软土的沉降过程预测方法
CN110927035A (zh) 一种低渗致密砂岩束缚水饱和度计算方法
CN107169271A (zh) 一种应变软化边坡的稳定性动态评价方法
da Fonseca et al. Recommended procedures to assess critical state locus from triaxial tests in cohesionless remoulded samples
CN104101562A (zh) 一种确定储层渗透率的方法
CN105606513B (zh) 确定裂隙岩体渗透率及孔隙度对有效应力敏感程度的方法
CN104237957A (zh) 一种核磁共振测井t2截止值的确定方法
CN105891089A (zh) 储层渗透率确定方法及装置
CN102915406A (zh) 径向流条件下油水相对渗透率曲线的计算方法
CN109100812A (zh) 基于核磁共振的岩石孔隙分形维数评价方法及装置
CN104297127B (zh) 基于三轴固结试验确定土体渗透系数与孔隙比关系的方法
CN105928827B (zh) 一种细粒土击实试验最佳含水率计算方法
CN104374679B (zh) 一种确定储层渗透率贡献值的方法
Zhang Limitations of suction‐controlled triaxial tests in the characterization of unsaturated soils
CN107894367A (zh) 一种输电线路湿陷性黄土本构模拟方法
CN107976390A (zh) 一种基于覆压孔渗实验确定平均孔喉比的计算方法
CN105844011A (zh) 一种基于毛管模型的渗透率计算方法
CN109612900A (zh) 一种储层岩石渗透率预测方法、装置及存储介质
Lovisa et al. An in-depth comparison of cv values determined using common curve-fitting techniques
CN112051292A (zh) 一种孔喉结构动态表征方法及装置
Cui et al. Modeling of three-dimensional single rough rock fissures: A study on flow rate and fractal parameters using the Weierstrass-Mandelbrot function

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