CN115467653A - 一种测井渗透率谱的获取方法 - Google Patents

一种测井渗透率谱的获取方法 Download PDF

Info

Publication number
CN115467653A
CN115467653A CN202211047482.7A CN202211047482A CN115467653A CN 115467653 A CN115467653 A CN 115467653A CN 202211047482 A CN202211047482 A CN 202211047482A CN 115467653 A CN115467653 A CN 115467653A
Authority
CN
China
Prior art keywords
permeability
logging
resistivity
equivalent
obtaining
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
CN202211047482.7A
Other languages
English (en)
Other versions
CN115467653B (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.)
Chengdu Univeristy of Technology
Original Assignee
Chengdu Univeristy of Technology
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 Chengdu Univeristy of Technology filed Critical Chengdu Univeristy of Technology
Priority to CN202211047482.7A priority Critical patent/CN115467653B/zh
Publication of CN115467653A publication Critical patent/CN115467653A/zh
Application granted granted Critical
Publication of CN115467653B publication Critical patent/CN115467653B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • E21B47/002Survey of boreholes or wells by visual inspection
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Geophysics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
  • Measurement Of Resistance Or Impedance (AREA)

Abstract

本申请提供一种测井渗透率谱的获取方法,包括:基于微电阻率成像测井获取的第i个极板测量的电导率Ci;根据井场水分析资料获取泥浆滤液电阻率Rmf;获取孔隙度曲折度T;根据所述电导率Ci、泥浆滤液电阻率Rmf、孔隙度曲折度T确定第i个极板对应的等效渗透率Ki;根据所述第i个极板对应的等效渗透率Ki确定微电阻率成像测井统计窗长内对应的的渗透率谱。本发明给出了参数更为简洁的渗透率谱计算模型,不需常规测井等参数,适用性更广。相较于常规测井解释的渗透率谱,本发明提供的方法可利用分辨率更高的微电阻率成像测井获得更为精细的渗透率,且能进行谱显示。

Description

一种测井渗透率谱的获取方法
技术领域
本发明涉及石油勘探储层评价技术领域,尤其涉及一种测井渗透率谱的获取方法。
背景技术
石油行业中,渗透率是储层评价、勘探开发的关键参数。常规测井计算的渗透率分辨率低,在复杂储层难以适用;岩心实验毛管压力实验测量的渗透率虽然精细,但成本高昂,难以获得连续的地层渗透率。目前,急需一套连续且精细的渗透率计算方法。
发明内容
本发明的目的在于解决上述现有技术存在的缺陷,提供一种于微电阻率成像测井的渗透率谱的计算方法。
一种测井渗透率谱的获取方法,包括:
基于微电阻率成像测井获取的第i个极板测量的电导率Ci
根据井场水分析资料获取泥浆滤液电阻率Rmf
获取孔隙度曲折度T;
根据所述电导率Ci、泥浆滤液电阻率Rmf、孔隙度曲折度T确定第i个极板对应的等效渗透率Ki
根据所述第i个极板对应的等效渗透率Ki确定基于微电阻率成像测井统计窗长内对应的渗透率谱。
进一步地,如上所述的测井渗透率谱的获取方法,所述根据所述电导率Ci、泥浆滤液电阻率Rmf、孔隙度曲折度T确定第i个极板对应的等效渗透率Ki包括采用以下公式获取等效渗透率Ki
Figure BDA0003820711510000021
进一步地,如上所述的测井渗透率谱的获取方法,所述渗透率谱基于毛管物理模型构建得到。
进一步地,如上所述的测井渗透率谱的获取方法,所述基于微电阻率成像测井获取的第i个极板测量的电导率Ci包括:
根据毛管物理模型,确定孔隙流体的电阻R以及单位体积内的岩石边长Lh
根据所述孔隙流体的电阻R以及岩石边长Lh,确定微电阻率成像测井单一极板测量的电导率Ci
进一步地,如上所述的测井渗透率谱的获取方法,所述根据所述孔隙流体的电阻R以及岩石边长Lh,确定基于微电阻率成像测井单一极板测量的电导率Ci包括根据以下公式获取电导率Ci
Figure BDA0003820711510000022
进一步地,如上所述的测井渗透率谱的获取方法,所述等效渗透率Ki基于微电阻率成像测井测量范围内等效孔喉毛管的平均流量q、以及孔喉半径r获取得到:
其中,微电阻率成像测井测量范围内等效孔喉毛管的平均流量为:
Figure BDA0003820711510000023
上式中,r为孔喉半径,ΔP为压力差,μ为流体粘度;
根据Darcy定律,流过等效毛管模型的流量可表示为:
Figure BDA0003820711510000024
等效毛管模型中流量相等,联立以上两式,并带入孔喉半径表达式
Figure BDA0003820711510000025
Figure BDA0003820711510000026
消去Lh,即可获取所述等效渗透率Ki
有益效果:
本发明提供的测井渗透率谱的获取方法,通过获取电导率Ci、泥浆滤液电阻率Rmf、孔隙度曲折度T即可确定第i个极板对应的等效渗透率Ki,并根据第i个极板对应的等效渗透率Ki最终确定对应的渗透率谱,该方法给出了参数更为简洁的渗透率谱计算模型,不需常规测井等参数,适用性更广。相较于常规测井解释的渗透率谱,本发明提供的方法可利用分辨率更高的微电阻率成像测井获得更为精细的渗透率,且能进行谱显示。
附图说明
图1为本申请基于微电阻率成像测井的岩石电阻率毛管物理模型示意图;
图2为毛管物理模型的等效电阻率图;
图3为本申请基于微电阻率成像测井的成像渗透率谱图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面对本发明中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1为本申请基于微电阻率成像测井的岩石电阻率毛管物理模型示意图,图2为毛管物理模型的等效电阻率图;如图1、图2所示,基于单位体积内的毛管物理模型,从骨架与孔喉部分构建单岩石电阻率毛管物理模型,设单位体积内,岩石边长为Lh(孔喉半径单位为微米,在标准单位运算下,可消去Lh),岩石的平均孔喉半径为r,曲折度为T,喉道长度为L,则L=LhT。微电阻率成像测井电极测量的电阻率可视为井壁冲洗带电阻率。若假定:①井壁冲洗带岩石孔隙内流体完全被泥浆滤液替代,不含残余油气,泥浆滤液电阻率为Rmf;②岩石骨架不导电,且不存在泥质附加导电性。假定单位体积内模型边长为Lh(标准单位运算下,可消去Lh),岩石的平均孔喉半径为r,曲折度为T,喉道长度为L,则L=LhT。基于单位体积内的毛管物理模型,构建渗透率谱。
根据毛管物理模型,由欧姆定律可知,孔隙流体的电阻R可表示为:
Figure BDA0003820711510000041
井壁岩石电阻可视为孔喉流体电阻与骨架电阻的并联,转换为微电阻率成像测井单一极板测量的电导率Ci
Figure BDA0003820711510000042
由式(2)可知,平均孔喉半径r可表示为:
Figure BDA0003820711510000043
据泊塞叶定律(Purcell,W.R.,1949),微电阻率成像测井测量范围内等效孔喉毛管的平均流量为:
Figure BDA0003820711510000044
式中,r为孔喉半径,ΔP为压力差,μ为流体粘度。
根据Darcy定律,流过等效毛管模型的流量可表示为:
Figure BDA0003820711510000045
等效毛管模型中流量相等,联立式(4)、(5),并带入孔喉半径表达式(3),消去Lh,基于微电阻率成像测井建立的等效渗透率为:
Figure BDA0003820711510000046
式中,Ki为微电阻率成像测井极板i计算的等效渗透率;Ci为微电阻率成像测井极板i测量电导率;Rmf为泥浆滤液电阻率,可通过井场水分析资料获取;T为孔喉曲折度,可通过岩心实验获取,也可取理论值
Figure BDA0003820711510000047
其中,Φ为孔隙度,可通过微电阻率成像测井计算获取,m为孔隙度指数,可通过岩心实验获取。
由渗透率谱模型可知,确定好微电阻率成像测井统计窗长,可获得窗长内微电阻率成像测井计算得到的渗透率K,并输出得到对应的渗透率谱。如图3所示,第一道为深度道;第二道为基于微电阻率成像测井的成像图道;第三道为计算出的平均渗透率道;第四道为计算出的渗透率谱道,x轴为渗透率值分量,y轴代表对应渗透率出现的频率,渗透率谱的形态反应着储层不同渗透率值的分布。
相较于常规测井解释的渗透率,本方法可利用分辨率更高的微电阻率成像测井获得更为精细的渗透率,且能进行谱显示。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (6)

1.一种测井渗透率谱的获取方法,其特征在于,包括:
基于微电阻率成像测井获取的第i个极板测量的电导率Ci
根据井场水分析资料获取泥浆滤液电阻率Rmf
获取孔隙度曲折度T;
根据所述电导率Ci、泥浆滤液电阻率Rmf、孔隙度曲折度T确定第i个极板对应的等效渗透率Ki
根据所述第i个极板对应的等效渗透率Ki确定基于微电阻率成像测井统计窗长内对应的渗透率谱。
2.根据权利要求1所述的测井渗透率谱的获取方法,其特征在于,所述根据所述电导率Ci、泥浆滤液电阻率Rmf、孔隙度曲折度T确定第i个极板对应的等效渗透率Ki包括采用以下公式获取等效渗透率Ki
Figure FDA0003820711500000011
3.根据权利要求1所述的测井渗透率谱的获取方法,其特征在于,所述渗透率谱基于毛管物理模型构建得到。
4.根据权利要求3所述的测井渗透率谱的获取方法,其特征在于,所述基于微电阻率成像测井获取的第i个极板测量的电导率Ci包括:
根据毛管物理模型,确定孔隙流体的电阻R以及单位体积内的岩石边长Lh
根据所述孔隙流体的电阻R以及岩石边长Lh,确定基于微电阻率成像测井单一极板测量的电导率Ci
5.据权利要求4所述的测井渗透率谱的获取方法,其特征在于,所述根据所述孔隙流体的电阻R以及岩石边长Lh,确定基于微电阻率成像测井单一极板测量的电导率Ci包括根据以下公式获取电导率Ci
Figure FDA0003820711500000012
6.根据权利要求4所述的测井渗透率谱的获取方法,其特征在于,所述等效渗透率Ki基于微电阻率成像测井测量范围内等效孔喉毛管的平均流量q、以及孔喉半径r获取得到:
其中,微电阻率成像测井测量范围内等效孔喉毛管的平均流量为:
Figure FDA0003820711500000021
上式中,r为孔喉半径,ΔP为压力差,μ为流体粘度;
根据Darcy定律,流过等效毛管模型的流量可表示为:
Figure FDA0003820711500000022
等效毛管模型中流量相等,联立以上两式,并带入孔喉半径表达式
Figure FDA0003820711500000023
Figure FDA0003820711500000024
消去Lh,即可获取所述等效渗透率Ki
CN202211047482.7A 2022-08-29 2022-08-29 一种测井渗透率谱的获取方法 Active CN115467653B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211047482.7A CN115467653B (zh) 2022-08-29 2022-08-29 一种测井渗透率谱的获取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211047482.7A CN115467653B (zh) 2022-08-29 2022-08-29 一种测井渗透率谱的获取方法

Publications (2)

Publication Number Publication Date
CN115467653A true CN115467653A (zh) 2022-12-13
CN115467653B CN115467653B (zh) 2024-06-18

Family

ID=84368564

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211047482.7A Active CN115467653B (zh) 2022-08-29 2022-08-29 一种测井渗透率谱的获取方法

Country Status (1)

Country Link
CN (1) CN115467653B (zh)

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5193059A (en) * 1990-06-06 1993-03-09 Western Atlas International Inc. Method for identifying and characterizing hydraulic units of saturated porous media: tri-kappa zoning process
US5335542A (en) * 1991-09-17 1994-08-09 Schlumberger Technology Corporation Integrated permeability measurement and resistivity imaging tool
US5463549A (en) * 1993-10-15 1995-10-31 Schlumberger Technology Corporation Method and apparatus for determining permeability of subsurface formations
US20040055745A1 (en) * 2001-07-20 2004-03-25 Baker Hughes Incorporated Method and apparatus for combined NMR and formation testing for assessing relative permeability with formation testing and nuclear magnetic resonance testing
US7054749B1 (en) * 2000-11-13 2006-05-30 O'meara Jr Daniel J Method for determining reservoir fluid volumes, fluid contacts, compartmentalization, and permeability in geological subsurface models
US20130308831A1 (en) * 2012-05-18 2013-11-21 Ingrain, Inc. Method And System For Estimating Rock Properties From Rock Samples Using Digital Rock Physics Imaging
CN104712330A (zh) * 2015-01-30 2015-06-17 中国地质大学(武汉) 一种测井渗透率解释方法
CN104912547A (zh) * 2014-03-11 2015-09-16 中国石油化工集团公司 应用电阻率成像测井资料连续定量评价储层非均质特征的方法
CN110439547A (zh) * 2019-08-15 2019-11-12 中国海洋石油集团有限公司 储集层微电阻成像生成孔隙度谱的方法
CN110955982A (zh) * 2019-12-18 2020-04-03 长江大学 一种变质岩储层渗透率计算方法、装置及计算机存储介质
US20200271820A1 (en) * 2017-08-22 2020-08-27 Halliburton Energy Services, Inc. A new porosity independent methodology for permeability prediction based on micro-resistivity images and laterolog resistivities
US20210223192A1 (en) * 2018-06-01 2021-07-22 Board Of Regents, The University Of Texas System Method for assessment of pore-throat size distribution and permeability in porous media
CN113640890A (zh) * 2021-08-31 2021-11-12 成都理工大学 基于核磁共振测井的储层气-水相对渗透率曲线构建方法、系统
CN113818867A (zh) * 2021-09-17 2021-12-21 成都理工大学 一种构建伪毛管压力曲线方法、系统、介质、设备及应用

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5193059A (en) * 1990-06-06 1993-03-09 Western Atlas International Inc. Method for identifying and characterizing hydraulic units of saturated porous media: tri-kappa zoning process
US5335542A (en) * 1991-09-17 1994-08-09 Schlumberger Technology Corporation Integrated permeability measurement and resistivity imaging tool
US5463549A (en) * 1993-10-15 1995-10-31 Schlumberger Technology Corporation Method and apparatus for determining permeability of subsurface formations
US7054749B1 (en) * 2000-11-13 2006-05-30 O'meara Jr Daniel J Method for determining reservoir fluid volumes, fluid contacts, compartmentalization, and permeability in geological subsurface models
US20040055745A1 (en) * 2001-07-20 2004-03-25 Baker Hughes Incorporated Method and apparatus for combined NMR and formation testing for assessing relative permeability with formation testing and nuclear magnetic resonance testing
US20130308831A1 (en) * 2012-05-18 2013-11-21 Ingrain, Inc. Method And System For Estimating Rock Properties From Rock Samples Using Digital Rock Physics Imaging
CN104912547A (zh) * 2014-03-11 2015-09-16 中国石油化工集团公司 应用电阻率成像测井资料连续定量评价储层非均质特征的方法
CN104712330A (zh) * 2015-01-30 2015-06-17 中国地质大学(武汉) 一种测井渗透率解释方法
US20200271820A1 (en) * 2017-08-22 2020-08-27 Halliburton Energy Services, Inc. A new porosity independent methodology for permeability prediction based on micro-resistivity images and laterolog resistivities
US20210223192A1 (en) * 2018-06-01 2021-07-22 Board Of Regents, The University Of Texas System Method for assessment of pore-throat size distribution and permeability in porous media
CN110439547A (zh) * 2019-08-15 2019-11-12 中国海洋石油集团有限公司 储集层微电阻成像生成孔隙度谱的方法
CN110955982A (zh) * 2019-12-18 2020-04-03 长江大学 一种变质岩储层渗透率计算方法、装置及计算机存储介质
CN113640890A (zh) * 2021-08-31 2021-11-12 成都理工大学 基于核磁共振测井的储层气-水相对渗透率曲线构建方法、系统
CN113818867A (zh) * 2021-09-17 2021-12-21 成都理工大学 一种构建伪毛管压力曲线方法、系统、介质、设备及应用

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
ARTUR POSENATO GARCIA, ETC: "A New Multiphysics Method for Simultaneous Assessment of Permeability and Saturation-Dependent Capillary Pressure in Hydrocarbon-Bearing Rocks", SOCIETY OF PETROLEUM ENGINEERS, 6 June 2019 (2019-06-06), pages 1 - 17 *
王亮等: "泥质砂岩液相渗透率计算新方法", 地球物理学报, vol. 58, no. 10, 31 October 2015 (2015-10-31), pages 3837 - 3844 *
郭志华: "骨架导电泥质岩石油水相对渗透率与电阻率关系实验研究", 中国优秀博士学位论文全文数据库(工程科技Ⅰ辑), 15 March 2021 (2021-03-15), pages 1 - 126 *

Also Published As

Publication number Publication date
CN115467653B (zh) 2024-06-18

Similar Documents

Publication Publication Date Title
US7363164B2 (en) Method of evaluating fluid saturation characteristics in a geological formation
US8301381B2 (en) Method for the characterization of geological formations
Wu et al. Numerical simulation of mud-filtrate invasion in deviated wells
CN106468172A (zh) 一种超低渗砂岩油藏低阻储层测井解释方法
CN104213899B (zh) 一种地层岩石骨架的测井识别方法
EP3108099B1 (en) Measuring behind casing hydraulic conductivity between reservoir layers
CN102540265B (zh) 一种基于网络模拟的多孔介质含水饱和度计算方法
CN110967771B (zh) 一种利用常规测井计算裂缝、溶蚀孔洞孔隙度的方法
WO2016115197A1 (en) Measuring inter-reservoir cross flow rate between adjacent reservoir layers form transient pressure tests
CN110043254B (zh) 一种基于电缆地层测试资料地层有效渗透率的获取方法
Yamada et al. Revisiting porosity analysis from electrical borehole images: integration of advanced texture and porosity analysis
CN112363242B (zh) 基于测录井融合的储层流体识别方法与装置
CN110554441A (zh) 一种基于钙质含量校正的致密砂岩含气饱和度计算方法
US8005619B2 (en) Method of determining reservoir parameters
CN112145165B (zh) 一种微裂缝-孔隙型储层动静态渗透率转换方法
CN106777649B (zh) 一种裂缝型储集层孔隙结构定量评价方法
Delfiner Three statistical pitfalls of Phi-K transforms
RU2580177C1 (ru) Способ определения изменений параметров пористой среды под действием загрязнителя
Hördt et al. Case histories of hydraulic conductivity estimation with induced polarization at the field scale
CN112115592A (zh) 一种基于岩石物理实验的大斜度/水平井电阻率校正方法
CN110017136B (zh) 一种基于视水层电阻率的水淹层识别与产水率预测方法
CN115467653A (zh) 一种测井渗透率谱的获取方法
Beik et al. Effects of highly laminated reservoirs on the performance of wireline and while-drilling formation-tester sampling with oval, focused, and conventional probe types
CN108828687A (zh) 一种基于电成像面孔率的渗透率计算方法
CN115421203A (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