CN113009115B - 一种多孔缝介质储层连通孔缝度的定量计算方法 - Google Patents

一种多孔缝介质储层连通孔缝度的定量计算方法 Download PDF

Info

Publication number
CN113009115B
CN113009115B CN202110264810.8A CN202110264810A CN113009115B CN 113009115 B CN113009115 B CN 113009115B CN 202110264810 A CN202110264810 A CN 202110264810A CN 113009115 B CN113009115 B CN 113009115B
Authority
CN
China
Prior art keywords
degree
medium reservoir
total
seam degree
seam
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.)
Active
Application number
CN202110264810.8A
Other languages
English (en)
Other versions
CN113009115A (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.)
Beijing Research Center of CNOOC China Ltd
CNOOC China Ltd
Original Assignee
Beijing Research Center of CNOOC China Ltd
CNOOC China 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 Beijing Research Center of CNOOC China Ltd, CNOOC China Ltd filed Critical Beijing Research Center of CNOOC China Ltd
Priority to CN202110264810.8A priority Critical patent/CN113009115B/zh
Publication of CN113009115A publication Critical patent/CN113009115A/zh
Application granted granted Critical
Publication of CN113009115B publication Critical patent/CN113009115B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/24Earth materials
    • 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/02Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
    • G01N23/04Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
    • G01N23/046Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N24/00Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
    • G01N24/08Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
    • G01N24/081Making measurements of geologic samples, e.g. measurements of moisture, pH, porosity, permeability, tortuosity or viscosity
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/02Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
    • G01N27/04Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance
    • G01N27/041Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance of a solid body
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/15Correlation function computation including computation of convolution operations
    • 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)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Pathology (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Food Science & Technology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Pulmonology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Computing Systems (AREA)
  • Medicinal Chemistry (AREA)
  • Remote Sensing (AREA)
  • Dispersion Chemistry (AREA)
  • Algebra (AREA)
  • Electrochemistry (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)

Abstract

本发明涉及一种多孔缝介质储层连通孔缝度的定量计算方法,包括:S1、基于岩心样品,测量得到岩心样品的总孔缝度;S2、开展电阻率实验,测量得到地层因素;S3、计算得到岩心样品的连通孔缝度;S4、基于岩心样品的影像实验,测量得到岩心样品的总孔缝度和连通孔缝度;S5、基于实验测量得到的总孔缝度和连通孔缝度,综合建立总孔缝度和连通孔缝度之间的函数关系;S6、通过测井数据的计算得到多孔缝介质储层的总孔缝度,得到使得计算的总孔缝度与实验测量得到的总孔缝度之间的误差尽量小的计算方式;S7、测量得到一系列测井数据,计算得到一口井中全井段连续的连通孔缝度。本发明可以准确计算多孔缝介质储层的连通孔缝度,具有极大的推广应用价值。

Description

一种多孔缝介质储层连通孔缝度的定量计算方法
技术领域
本发明涉及油气勘探技术领域,具体涉及一种多孔缝介质储层连通孔缝度的定量计算方法。
背景技术
连通孔缝度是岩石中连通孔隙及裂缝的孔隙体积占岩石总体积的百分比。在不发育裂缝的均质储层,孔隙类型较单一,孔隙之间的连通性较好,孤立的孔隙相对较少,连通孔缝度即为连通孔隙度,且与总孔隙度相近;因此,基于总孔隙度,可以准确计算均质储层的渗透率、饱和度等一系列的储层参数。
但在非均质储层,尤其是在多孔缝介质储层,除了发育不同类型的孔隙之外,还发育不同类型的裂缝;从孔隙及裂缝是否连通的角度,可以将孔隙及裂缝分为连通孔隙、孤立孔隙、连通裂缝、孤立裂缝等。由于孤立的孔隙及裂缝所形成的孤立孔缝度在储层中占到一定的比例,从而导致连通的孔隙及裂缝所形成的连通孔缝度与总孔缝度存在较大差异。因此,基于总孔缝度计算多孔缝介质储层的渗透率、饱和度等一系列的储层参数时会存在较大误差;只有基于连通孔缝度,才能准确计算多孔缝介质储层的渗透率、饱和度等一系列的储层参数,从而才能精确表征多孔缝介质储层。
随着油气勘探开发不断向深层、超深层发展,多孔缝介质油气藏逐渐成为增储上产的主力,而连通孔缝度的准确计算是多孔缝介质储层精确表征的关键参数,直接影响勘探、开发效果和经济效益。目前业界尚未形成连通孔缝度的定量计算方法,因此,针对多孔缝介质储层形成一种连通孔缝度的定量计算方法,能够有效地提高多孔缝介质油气藏勘探、开发的效益。
发明内容
针对上述现有技术的不足,本发明旨在提供一种连通孔缝度的定量计算方法,以便能够准确计算多孔缝介质储层的连通孔缝度。
本发明提出一种多孔缝介质储层连通孔缝度的定量计算方法,所述方法包括以下步骤:
S1、基于多孔缝介质储层的岩心样品,开展物性分析实验,测量得到岩心样品的总孔缝度;
S2、开展岩心样品的电阻率实验,测量得到所述岩心样品的地层因素;
S3、基于测量得到的所述总孔缝度和地层因素,计算得到所述岩心样品的连通孔缝度;
S4、基于所述岩心样品的影像实验,测量得到所述岩心样品的总孔缝度和连通孔缝度;
S5、基于多孔缝介质储层不同尺度岩心样品的物性分析及电阻率实验、影像实验测量得到的总孔缝度和连通孔缝度,综合建立连通孔缝度和总孔缝度之间的函数关系;
S6、通过测井数据的计算得到多孔缝介质储层的总孔缝度,得到使得计算的总孔缝度与实验测量得到的总孔缝度之间的误差尽量小的计算方式;
S7、基于油气田一口井中多孔缝介质储层测井仪器测量得到的一系列测井数据,根据建立的所述连通孔缝度和总孔缝度之间的函数关系及所述计算方式,计算得到一口井中多孔缝介质储层全井段连续的连通孔缝度。
根据本发明的一种实施方式,所述岩心样品的影像实验,包括CT扫描实验和/或核磁共振实验。
根据本发明的一种实施方式,开展所述岩心样品的电阻率实验时,使所述岩心样品100%饱和该油气田地层条件下的水,得到100%含水纯岩石电阻率R0;100%含水纯岩石电阻率R0与地层水电阻率Rw的比值,得到地层因素F,
Figure BDA0002971965040000021
式中:F为地层因素,无量纲;R0为100%含水纯岩石电阻率,Ω·m;Rw为地层水电阻率,Ω·m。
根据本发明的一种实施方式,在步骤S3中,测量得到油气田多孔缝介质储层一定数量岩心样品的总孔缝度
Figure BDA0002971965040000022
和地层因素F,以/>
Figure BDA0002971965040000023
为横坐标,总孔缝度/>
Figure BDA0002971965040000024
为纵坐标,以表征岩石的孔隙及裂缝的几何形状,线性回归可以得到直线的斜率x。
根据本发明的一种实施方式,针对油气田多孔缝介质储层的每一块岩心样品,计算得到x/(F+x-1),即得到该油气田多孔缝介质储层每一块岩心样品的连通孔缝度。
根据本发明的一种实施方式,步骤S6中,以不同尺度岩心样品实验测量得到的总孔缝度为标准,使计算的总孔缝度与实验测量得到的总孔缝度之间的误差尽量小。
根据本发明的一种实施方式,所述测井数据包括自然伽马、自然电位、深电阻率、中电阻率、浅电阻率、体积密度、中子孔隙度、纵波时差中的全部或多个的组合。
根据本发明的一种实施方式,在步骤S6中,使得计算的总孔缝度与实验测量得到的总孔缝度之间的误差尽量小时,根据步骤S5中建立的总孔缝度和连通孔缝度之间的函数关系,相应得到计算的总孔缝度与连通孔缝度之间的函数关系。
根据本发明的一种实施方式,在步骤S7中,基于油气田一口井中多孔缝介质储层测井仪器测量得到的一系列测井数据,根据所述计算方式,计算得到这口井中多孔缝介质储层的总孔缝度,并基于计算的所述总孔缝度与连通孔缝度之间的函数关系,计算得到这口井中多孔缝介质储层的连通孔缝度。
根据本发明的一种实施方式,所述方法还包括:通过Geolog、GeoFrame、Forward、Lead中的一种或多种软件来定量评价油气田井中多孔缝介质储层的连通孔缝度。
本发明由于采取以上技术方案,其具有以下优点:本发明的方法避免开展大量的岩心实验,能够有效地节约成本,具有较强的经济性。本发明的方法可以准确计算多孔缝介质储层的连通孔缝度,在保证多孔缝介质储层连通孔缝度定量计算结果准确性的同时,提供了一种简单、具有极大推广应用价值的新方法。
附图说明
图1为本发明一实施例M油田多孔缝介质储层岩心柱子显示所发育的孔隙示意图;
图2为本发明一实施例M油田多孔缝介质储层岩心铸体薄片显示所发育的裂缝示意图;
图3为本发明一实施例多孔缝介质储层发育连通孔隙、连通裂缝、孤立孔隙、孤立裂缝的理论模型示意图;
图4为本发明一实施例M油田多孔缝介质储层(1-总孔缝度)/(地层因素-1)与总孔缝度的交会示意图;
图5为本发明一实施例M油田多孔缝介质储层岩心CT扫描成果图(总孔缝度为4.20%、连通孔缝度为1.25%)示意图;
图6为本发明一实施例M油田多孔缝介质储层岩心CT扫描成果图(总孔缝度为3.75%、连通孔缝度为0.89%)示意图;
图7为本发明一实施例M油田多孔缝介质储层岩心核磁共振实验成果图(总孔缝度为4.38%、连通孔缝度为2.63%)示意图;
图8为本发明一实施例M油田多孔缝介质储层岩心核磁共振实验成果图(总孔缝度为4.36%、连通孔缝度为0.69%)示意图;
图9为本发明一实施例M油田多孔缝介质储层岩心电阻率实验、CT扫描实验和核磁共振实验测量得到的总孔缝度与连通孔缝度的交会示意图;
图10为本发明一实施例M油田A井多孔缝介质储层连通孔缝度测井定量计算成果示意图;
图11为本发明一实施例多孔缝介质储层连通孔缝度的定量计算方法流程图。
具体实施方式
以下将结合附图对本发明的较佳实施例进行详细说明,以便更清楚理解本发明的目的、特点和优点。应理解的是,附图所示的实施例并不是对本发明范围的限制,而只是为了说明本发明技术方案的实质精神。
本发明实施方式一种多孔缝介质储层连通孔缝度的定量计算方法,如图11所示,主要包括以下步骤:
S1、基于多孔缝介质储层的岩心样品,开展物性分析实验,测量得到岩心样品的总孔缝度;
S2、开展岩心样品的电阻率实验,测量得到所述岩心样品的地层因素;
S3、基于测量得到的所述总孔缝度和地层因素,计算得到所述岩心样品的连通孔缝度;
S4、基于所述岩心样品的影像实验,测量得到所述岩心样品的总孔缝度和连通孔缝度;
S5、基于多孔缝介质储层不同尺度岩心样品的物性分析及电阻率实验、影像实验测量得到的总孔缝度和连通孔缝度,综合建立连通孔缝度和总孔缝度之间的函数关系;
S6、通过测井数据的计算得到多孔缝介质储层的总孔缝度,得到使得计算的总孔缝度与实验测量得到的总孔缝度之间的误差尽量小的计算方式;
S7、基于油气田一口井中多孔缝介质储层测井仪器测量得到的一系列测井数据,根据建立的所述连通孔缝度和总孔缝度之间的函数关系及所述计算方式,计算得到一口井中多孔缝介质储层全井段连续的连通孔缝度。
本发明的上述方法通过实验和计算相互综合,可以准确计算多孔缝介质储层的连通孔缝度,在保证多孔缝介质储层连通孔缝度定量计算结果准确性的同时,提供了一种简单、具有极大推广应用价值的新方法。
本发明一实施方式提出一种多孔缝介质储层连通孔缝度的定量计算方法,主要包括以下步骤:1)基于多孔缝介质储层的岩心样品,开展物性分析实验,可以测量得到岩心样品的总孔缝度;2)开展岩心样品的电阻率实验,可以测量得到岩心样品的地层因素;3)基于步骤1)得到的总孔缝度和步骤2)得到的地层因素,可以计算得到岩心样品的连通孔缝度;4)基于岩心样品的CT扫描实验,可以测量得到总孔缝度和连通孔缝度;5)基于岩心样品的核磁共振实验,可以测量得到总孔缝度和连通孔缝度;6)基于多孔缝介质储层不同尺度岩心样品的物性分析及电阻率实验、CT扫描实验和核磁共振实验测量得到的总孔缝度和连通孔缝度,综合建立总孔缝度和连通孔缝度之间的函数关系;7)测井仪器可以测量得到一系列的测井数据,通过测井数据的计算,可以计算得到多孔缝介质储层的总孔缝度,使计算的总孔缝度与实验测量得到的总孔缝度之间的误差尽量小;8)采用步骤6)中总孔缝度和连通孔缝度的函数关系,从而可以计算得到一口井中多孔缝介质储层全井段连续的连通孔缝度。
需要说明的是,上述步骤并非都是必要设置,可根据具体情况选择合适步骤及顺序应用。
根据本发明的一种实施方式,多孔缝介质储层连通孔缝度的定量计算方法,包括以下具体步骤:
1)基于某油气田多孔缝介质储层一定深度下的岩心样品,开展岩心样品的物性分析实验,从而可以测量得到每一块岩心样品中孔隙及裂缝的孔隙体积占岩心样品总体积的百分比,即岩心样品的总孔缝度
Figure BDA0002971965040000051
2)基于该油气田多孔缝介质储层一定深度下的岩心样品,开展岩心样品的电阻率实验,使岩心样品100%饱和该油气田地层条件下的水,并测量得到100%含水纯岩石电阻率R0;100%含水纯岩石电阻率R0与地层水电阻率Rw的比值,可以得到地层因素F,见式(1)。
Figure BDA0002971965040000052
式中:F为地层因素,无量纲;R0为100%含水纯岩石电阻率,Ω·m;Rw为地层水电阻率,Ω·m。
3)基于步骤1)和步骤2)所测量得到该油气田多孔缝介质储层一定数量岩心样品的总孔缝度
Figure BDA0002971965040000053
和地层因素F,以/>
Figure BDA00029719650400000511
为横坐标,总孔缝度/>
Figure BDA0002971965040000054
为纵坐标,线性回归可以得到直线的斜率x。
4)基于步骤3)所得到直线的斜率x,针对该油气田多孔缝介质储层的每一块岩心样品,可以计算得到x/(F+x-1);因此,基于岩心样品的物性分析及电阻率实验,可以计算得到该油气田多孔缝介质储层每一块岩心样品的连通孔缝度
Figure BDA0002971965040000055
5)基于该油气田多孔缝介质储层一定深度下的岩心样品,开展岩心样品的CT扫描实验,从而可以测量得到岩心样品的总孔缝度
Figure BDA0002971965040000056
和连通孔缝度/>
Figure BDA0002971965040000057
6)基于多孔缝介质储层一定深度下的岩心样品,开展岩心样品的核磁共振实验,从而可以测量得到岩心样品的总孔缝度
Figure BDA0002971965040000058
和连通孔缝度/>
Figure BDA0002971965040000059
7)基于多孔缝介质储层不同尺度岩心样品的物性分析及电阻率实验、CT扫描实验和核磁共振实验测量得到的总孔缝度和连通孔缝度,综合不同类型的实验数据,可以建立连通孔缝度和总孔缝度之间的函数关系,见式(2)。
Figure BDA00029719650400000510
8)针对该油气田多孔缝介质储层,测井仪器可以测量得到一系列的测井数据(自然伽马、自然电位、深电阻率、中电阻率、浅电阻率、体积密度、中子孔隙度、纵波时差);通过测井数据的计算,可以计算得到多孔缝介质储层的总孔缝度
Figure BDA0002971965040000061
并以不同尺度岩心样品实验测量得到的总孔缝度/>
Figure BDA0002971965040000062
为标准,使计算的总孔缝度/>
Figure BDA0002971965040000063
与实验测量得到的总孔缝度/>
Figure BDA0002971965040000064
之间的误差尽量小,即计算的总孔缝度/>
Figure BDA0002971965040000065
约等于实验测量得到的总孔缝度/>
Figure BDA0002971965040000066
从而将式(2)变换为式(3)。
Figure BDA0002971965040000067
9)基于该油气田一口井中多孔缝介质储层测井仪器测量得到的一系列测井数据(自然伽马、自然电位、深电阻率、中电阻率、浅电阻率、体积密度、中子孔隙度、纵波时差),可以计算得到这口井中多孔缝介质储层的总孔缝度
Figure BDA0002971965040000068
并基于式(3),可以计算得到这口井中多孔缝介质储层的连通孔缝度/>
Figure BDA0002971965040000069
实施例
下面以M油田多孔缝介质储层为例,结合附图详细阐述多孔缝介质储层连通孔缝度的定量计算过程。
1)M油田多孔缝介质储层岩心柱子显示所发育的孔隙和铸体薄片显示所发育的裂缝分别见图1、图2;基于M油田多孔缝介质储层一定深度下的岩心样品,开展岩心样品的物性分析实验,从而可以测量得到每一块岩心样品中孔隙及裂缝的孔隙体积占岩心样品总体积的百分比,即7块岩心样品的总孔缝度
Figure BDA00029719650400000610
见表1。
表1M油田多孔缝介质储层岩心样品物性分析实验测量得到的总孔缝度
Figure BDA00029719650400000611
2)基于M油田多孔缝介质储层一定深度下的岩心样品,开展岩心样品的电阻率实验,使岩心样品100%饱和M油田地层条件下的水,并测量得到100%含水纯岩石电阻率R0;100%含水纯岩石电阻率R0与地层水电阻率Rw的比值,基于式(1)可以计算得到地层因素F,见表2。
Figure BDA00029719650400000612
式中:F为地层因素,无量纲;R0为100%含水纯岩石电阻率,Ω·m;Rw为地层水电阻率,Ω·m。
表2M油田多孔缝介质储层岩心样品物性分析及电阻率实验测量得到的总孔缝度与地层因素
Figure BDA0002971965040000071
3)多孔缝介质储层发育连通孔隙、连通裂缝、孤立孔隙、孤立裂缝的理论模型见图3;基于步骤1)和步骤2)所测量得到的M油田多孔缝介质储层7块岩心样品的总孔缝度
Figure BDA0002971965040000072
和地层因素F,以/>
Figure BDA0002971965040000073
为横坐标,总孔缝度/>
Figure BDA0002971965040000074
为纵坐标,线性回归可以得到直线的斜率x为18.25,见表3和图4。
表3M油田多孔缝介质储层岩心样品的
Figure BDA0002971965040000075
Figure BDA0002971965040000076
4)基于步骤3)所得到直线的斜率x为18.25,针对M油田多孔缝介质储层的每一块岩心样品,可以计算得到x/(F+x-1),即每一块岩心样品的连通孔缝度
Figure BDA0002971965040000077
见表4。
表4M油田多孔缝介质储层岩心样品连通孔缝度的计算结果
Figure BDA0002971965040000078
5)基于M油田多孔缝介质储层一定深度下的岩心样品,开展岩心样品的CT扫描实验,从而可以测量得到2块岩心样品的总孔缝度
Figure BDA0002971965040000079
和连通孔缝度/>
Figure BDA00029719650400000710
图5中CT扫描实验测定的总孔缝度/>
Figure BDA00029719650400000711
为4.20%、连通孔缝度/>
Figure BDA00029719650400000712
为1.25%;图6中CT扫描实验测定的总孔缝度
Figure BDA0002971965040000081
为3.75%、连通孔缝度/>
Figure BDA0002971965040000082
为0.89%。
6)基于M油田多孔缝介质储层一定深度下的岩心样品,开展岩心样品的核磁共振实验,从而可以测量得到2块岩心样品的总孔缝度
Figure BDA0002971965040000083
和连通孔缝度/>
Figure BDA0002971965040000084
图7中核磁共振实验测定的总孔缝度/>
Figure BDA0002971965040000085
为4.38%、连通孔缝度/>
Figure BDA0002971965040000086
为2.63%;图8中核磁共振实验测定的总孔缝度/>
Figure BDA0002971965040000087
为4.36%、连通孔缝度/>
Figure BDA0002971965040000088
为0.69%。
7)基于M油田多孔缝介质储层不同尺度岩心样品的物性分析及电阻率实验、CT扫描实验和核磁共振实验测量得到的总孔缝度和连通孔缝度,见表5;
表5M油田多孔缝介质储层岩心样品不同类型实验测量得到的总孔缝度和连通孔缝度
Figure BDA0002971965040000089
综合不同类型的实验数据,可以建立连通孔缝度和总孔缝度之间的函数关系,见图9和式(2)。
Figure BDA00029719650400000810
其中,
Figure BDA00029719650400000811
为连通孔缝度。
8)针对M油田多孔缝介质储层,测井仪器可以测量得到一系列的测井数据(自然伽马、自然电位、深电阻率、中电阻率、浅电阻率、体积密度、中子孔隙度、纵波时差);通过测井数据的计算,可以计算得到多孔缝介质储层的总孔缝度
Figure BDA00029719650400000812
并以不同尺度岩心样品实验测量得到的总孔缝度/>
Figure BDA00029719650400000813
为标准,使计算的总孔缝度/>
Figure BDA00029719650400000814
与实验测量得到的总孔缝度/>
Figure BDA00029719650400000815
之间的误差尽量小,即计算的总孔缝度/>
Figure BDA00029719650400000816
约等于实验测量得到的总孔缝度/>
Figure BDA00029719650400000817
从而将式(2)变换为式(3)。
Figure BDA00029719650400000818
9)基于M油田一口井中多孔缝介质储层测井仪器测量得到的一系列测井数据(自然伽马、自然电位、深电阻率、中电阻率、浅电阻率、体积密度、中子孔隙度、纵波时差),可以计算得到这口井中多孔缝介质储层的总孔缝度
Figure BDA00029719650400000819
并基于式(3),可以计算得到这口井中多孔缝介质储层的连通孔缝度/>
Figure BDA00029719650400000820
本发明可以通过Geolog、GeoFrame、Forward、Lead等软件来定量评价M油田A井多孔缝介质储层的连通孔缝度。
如图10所示,在A井多孔缝介质储层连通孔缝度测井定量计算成果图中,第1道是地层测量深度;第2道是自然伽马和自然电位,表示该地层的岩性特征;第3道是深、浅电阻率测井曲线,刻画该地层的电性特征;第4道是体积密度、中子孔隙度和声波时差,反映该地层的物性特征;第5道是基于上述测井曲线所计算的总孔缝度与岩心分析总孔缝度,二者之间的平均绝对误差为1.21%;第6道是测井计算连通孔缝度与岩心分析连通孔缝度,二者之间的平均绝对误差为0.83%。
上述各实施例仅用于说明本发明,其中方法的实施步骤等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。

Claims (8)

1.一种多孔缝介质储层连通孔缝度的定量计算方法,其特征在于,所述方法包括以下步骤:
S1、基于多孔缝介质储层的岩心样品,开展物性分析实验,测量得到岩心样品的总孔缝度;
S2、开展岩心样品的电阻率实验,测量得到所述岩心样品的地层因素;
S3、基于测量得到的所述总孔缝度和地层因素,计算得到所述岩心样品的连通孔缝度;
S4、基于所述岩心样品的影像实验,测量得到所述岩心样品的总孔缝度和连通孔缝度;
S5、基于多孔缝介质储层不同尺度岩心样品的物性分析及电阻率实验、影像实验测量得到的总孔缝度和连通孔缝度,综合建立连通孔缝度和总孔缝度之间的函数关系;
S6、通过测井数据的计算得到多孔缝介质储层的总孔缝度,得到使得计算的总孔缝度与实验测量得到的总孔缝度之间的误差尽量小的计算方式;
S7、基于油气田一口井中多孔缝介质储层测井仪器测量得到的一系列测井数据,根据建立的所述连通孔缝度和总孔缝度之间的函数关系及所述计算方式,计算得到一口井中多孔缝介质储层全井段连续的连通孔缝度;
在步骤S6中,使得计算的总孔缝度与实验测量得到的总孔缝度之间的误差尽量小时,根据步骤S5中建立的总孔缝度和连通孔缝度之间的函数关系,相应得到计算的总孔缝度与连通孔缝度之间的函数关系;
在步骤S7中,基于油气田一口井中多孔缝介质储层测井仪器测量得到的一系列测井数据,根据所述计算方式,计算得到这口井中多孔缝介质储层的总孔缝度,并基于计算的所述总孔缝度与连通孔缝度之间的函数关系,计算得到这口井中多孔缝介质储层的连通孔缝度。
2.根据权利要求1所述的多孔缝介质储层连通孔缝度的定量计算方法,其特征在于,所述岩心样品的影像实验,包括CT扫描实验和/或核磁共振实验。
3.根据权利要求1或2所述的多孔缝介质储层连通孔缝度的定量计算方法,其特征在于,开展所述岩心样品的电阻率实验时,使所述岩心样品100%饱和该油气田地层条件下的水,得到100%含水纯岩石电阻率R0;100%含水纯岩石电阻率R0与地层水电阻率Rw的比值,得到地层因素F,
Figure FDA0004064597140000011
式中:F为地层因素,无量纲;R0为100%含水纯岩石电阻率,Ω·m;Rw为地层水电阻率,Ω·m。
4.根据权利要求3所述的多孔缝介质储层连通孔缝度的定量计算方法,其特征在于,在步骤S3中,测量得到油气田多孔缝介质储层一定数量岩心样品的总孔缝度和地层因素F,以
Figure FDA0004064597140000021
为横坐标,总孔缝度为纵坐标,以表征岩石的孔隙裂缝的几何形状,线性回归得到直线的斜率x。
5.根据权利要求4所述的多孔缝介质储层连通孔缝度的定量计算方法,其特征在于,针对油气田多孔缝介质储层的每一块岩心样品,计算得到x/(F+x-1),即得到该油气田多孔缝介质储层每一块岩心样品的连通孔缝度。
6.根据权利要求1或2或4或5所述的多孔缝介质储层连通孔缝度的定量计算方法,其特征在于,步骤S6中,以不同尺度岩心样品实验测量得到的总孔缝度为标准,使计算的总孔缝度与实验测量得到的总孔缝度之间的误差尽量小。
7.根据权利要求1或2或4或5所述的多孔缝介质储层连通孔缝度的定量计算方法,其特征在于,所述测井数据包括自然伽马、自然电位、深电阻率、中电阻率、浅电阻率、体积密度、中子孔隙度、纵波时差中的全部或多个的组合。
8.根据权利要求1或2或4或5所述的多孔缝介质储层连通孔缝度的定量计算方法,其特征在于,所述方法还包括:通过Geolog、GeoFrame、Forward、Lead中的一种或多种软件来定量评价油气田井中多孔缝介质储层的连通孔缝度。
CN202110264810.8A 2021-03-11 2021-03-11 一种多孔缝介质储层连通孔缝度的定量计算方法 Active CN113009115B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110264810.8A CN113009115B (zh) 2021-03-11 2021-03-11 一种多孔缝介质储层连通孔缝度的定量计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110264810.8A CN113009115B (zh) 2021-03-11 2021-03-11 一种多孔缝介质储层连通孔缝度的定量计算方法

Publications (2)

Publication Number Publication Date
CN113009115A CN113009115A (zh) 2021-06-22
CN113009115B true CN113009115B (zh) 2023-07-14

Family

ID=76405066

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110264810.8A Active CN113009115B (zh) 2021-03-11 2021-03-11 一种多孔缝介质储层连通孔缝度的定量计算方法

Country Status (1)

Country Link
CN (1) CN113009115B (zh)

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102979517B (zh) * 2012-12-04 2015-06-17 中国海洋石油总公司 一种复杂油气储层饱和度定量评价方法
CN104727813B (zh) * 2015-02-12 2017-08-29 中国石油大学(华东) 一种源储共生型地层的孔隙度测量方法
DE102016201422A1 (de) * 2016-01-29 2017-08-17 Rheinisch-Westfälische Technische Hochschule (Rwth) Aachen System und Verfahren zur Messung der Porosität und Porenraumkonnektivität einer Gesteinsprobe
CN108802073B (zh) * 2017-05-05 2021-03-12 中国石油化工股份有限公司 一种基于数字岩心的岩石电性参数获取方法及装置
CN111751259A (zh) * 2020-05-12 2020-10-09 中国石油天然气股份有限公司 确定不规则岩石样品有效孔隙度的方法及装置

Also Published As

Publication number Publication date
CN113009115A (zh) 2021-06-22

Similar Documents

Publication Publication Date Title
EP1729151B1 (en) Method and apparatus for measuring the wettability of geological formations
Glover et al. Streaming potential coupling coefficient of quartz glass bead packs: dependence on grain diameter, pore size, and pore throat radius
CN105659074B (zh) 三轴核磁共振测试仪器
CN104213899B (zh) 一种地层岩石骨架的测井识别方法
CN111007230B (zh) 定量评价陆相湖盆低孔隙度致密油储层含油量的方法
CN111337408B (zh) 一种利用低场核磁共振设备测试岩石裂缝孔隙度的方法
CN102175832A (zh) 一种确定典型储层最佳饱和度计算模型的方法
US8005619B2 (en) Method of determining reservoir parameters
CN110320139A (zh) 缝洞型储层裂缝孔隙度定量评价方法及系统
CN103792575B (zh) 一种去骨架影响的致密砂岩气层识别评价方法
CN111220522A (zh) 一种高泥质细粒沉积物水合物饱和度的岩心刻度测井计算方法
CN112145165B (zh) 一种微裂缝-孔隙型储层动静态渗透率转换方法
CN110700820A (zh) 松辽盆地北部致密油储层甜点分类方法
CN106897531A (zh) 一种低渗透石灰岩储层渗透率的定量评价方法
CN110593857B (zh) 油润湿致密砂岩核磁共振测井流体识别方法及其图版
CN106777649B (zh) 一种裂缝型储集层孔隙结构定量评价方法
Liu et al. Determination of the pore-throat limits for water imbibition in tight sandstone reservoirs through NMR analysis
CN111381292A (zh) 一种预测砂岩含烃储层的测井解释方法与装置
CN113009115B (zh) 一种多孔缝介质储层连通孔缝度的定量计算方法
CN112115592A (zh) 一种基于岩石物理实验的大斜度/水平井电阻率校正方法
He et al. Radial resistivity measurement method for cylindrical core samples
Wang et al. Measurements of the seismoelectric responses in a synthetic porous rock
Ouzzane et al. Application of NMR T2 relaxation to drainage capillary pressure in vuggy carbonate reservoirs
CN112878999B (zh) 各向异性地层含水饱和度的计算方法及装置
Ou Experimental study on the anisotropy of unconventional tight reservoirs: joint ultrasonic and electrical measurements under pressure

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