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

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

Info

Publication number
CN113009115A
CN113009115A CN202110264810.8A CN202110264810A CN113009115A CN 113009115 A CN113009115 A CN 113009115A CN 202110264810 A CN202110264810 A CN 202110264810A CN 113009115 A CN113009115 A CN 113009115A
Authority
CN
China
Prior art keywords
porosity
medium reservoir
total
porous
core sample
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
CN202110264810.8A
Other languages
English (en)
Other versions
CN113009115B (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

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/24—Earth materials
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00—Investigating 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/02—Investigating 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/04—Investigating 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/046—Investigating 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]
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
    • G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
    • G01N24/081—Making measurements of geologic samples, e.g. measurements of moisture, pH, porosity, permeability, tortuosity or viscosity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/02—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
    • G01N27/04—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance
    • G01N27/041—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance of a solid body
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/15—Correlation function computation including computation of convolution operations
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • Theoretical Computer Science (AREA)
  • Geology (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Environmental & Geological Engineering (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Pulmonology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Computing Systems (AREA)
  • Electrochemistry (AREA)
  • Dispersion Chemistry (AREA)
  • Medicinal Chemistry (AREA)
  • Algebra (AREA)
  • Food Science & Technology (AREA)
  • Remote Sensing (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
  • Geophysics And Detection Of Objects (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 (10)

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

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102979517A (zh) * 2012-12-04 2013-03-20 中国海洋石油总公司 一种复杂油气储层饱和度定量评价方法
CN104727813A (zh) * 2015-02-12 2015-06-24 中国石油大学(华东) 一种源储共生型地层的孔隙度测量方法
WO2017129812A1 (de) * 2016-01-29 2017-08-03 Rwth Aachen System und verfahren zur messung der porosität und porenraumkonnektivität einer gesteinsprobe
CN108802073A (zh) * 2017-05-05 2018-11-13 中国石油化工股份有限公司 一种基于数字岩心的岩石电性参数获取方法及装置
CN111751259A (zh) * 2020-05-12 2020-10-09 中国石油天然气股份有限公司 确定不规则岩石样品有效孔隙度的方法及装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102979517A (zh) * 2012-12-04 2013-03-20 中国海洋石油总公司 一种复杂油气储层饱和度定量评价方法
CN104727813A (zh) * 2015-02-12 2015-06-24 中国石油大学(华东) 一种源储共生型地层的孔隙度测量方法
WO2017129812A1 (de) * 2016-01-29 2017-08-03 Rwth Aachen System und verfahren zur messung der porosität und porenraumkonnektivität einer gesteinsprobe
CN108802073A (zh) * 2017-05-05 2018-11-13 中国石油化工股份有限公司 一种基于数字岩心的岩石电性参数获取方法及装置
CN111751259A (zh) * 2020-05-12 2020-10-09 中国石油天然气股份有限公司 确定不规则岩石样品有效孔隙度的方法及装置

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
孙泽 等: ""基于三维CT扫描图像的致密砂岩孔隙数据提取方法"", 《华北科技学院学报》 *
李雄炎 等: ""复杂储层连通孔隙度评价与渗透率定量计算方法"", 《石油地球物理勘探》 *
李雄炎等: "高精度胶结指数模型的建立与应用", 《石油学报》 *
杜玉荣等: "三孔隙模型在评价塔河油田碳酸盐岩储层的应用", 《国外测井技术》 *
闫伟林等: "基于导电孔隙的中基性火山岩储层含气饱和度解释模型", 《吉林大学学报(地球科学版)》 *

Also Published As

Publication number Publication date
CN113009115B (zh) 2023-07-14

Similar Documents

Publication Publication Date Title
US7532983B2 (en) Method and apparatus for measuring the wettability of geological formations
CN106468172A (zh) 一种超低渗砂岩油藏低阻储层测井解释方法
RU2731842C1 (ru) Способы и системы для определения объемной плотности, пористости и распределения размера пор подповерхностной формации
EA011046B1 (ru) Способ определения водонасыщенности подземной формации
US8005619B2 (en) Method of determining reservoir parameters
CN112145165B (zh) 一种微裂缝-孔隙型储层动静态渗透率转换方法
CN107688037A (zh) 一种利用核磁测井t2分布确定井下岩石粒度曲线的方法
CN110593857B (zh) 油润湿致密砂岩核磁共振测井流体识别方法及其图版
CN111381292B (zh) 一种预测砂岩含烃储层的测井解释方法与装置
CN111220522A (zh) 一种高泥质细粒沉积物水合物饱和度的岩心刻度测井计算方法
CN110344822B (zh) 含钙质致密储层含水饱和度的确定方法
CN113216945A (zh) 一种致密砂岩储层渗透性定量评价方法
CN110410058B (zh) 一种校正岩心实验结果刻度二维核磁测井的方法
CN112115592A (zh) 一种基于岩石物理实验的大斜度/水平井电阻率校正方法
CN113009115A (zh) 一种多孔缝介质储层连通孔缝度的定量计算方法
CN113006779B (zh) 一种砾岩地层岩石可钻性级值的评价方法和校正方法
CN115391739A (zh) 一种裂缝渗透率的定量计算方法和系统
CN112878999B (zh) 各向异性地层含水饱和度的计算方法及装置
US11754518B2 (en) Methods and systems to determine tortuosity of rock and fluids in porous media
CN113720745A (zh) 含碳屑碎屑岩储层地球物理测井计算孔隙度的方法
CN114428365A (zh) 一种基于含气特征的致密砂岩气层测井识别方法及装置
CN114086938B (zh) 一种非均质性砂岩储层的含气饱和度预测方法
CN114428049B (zh) 一种计算古老碳酸盐岩储层沥青含量的方法
RU2778620C1 (ru) Метод нейтрон-нейтронной цементометрии - ннк-ц для контроля качества цементирования облегченными и обычными цементами строящихся скважин и состояния цементного камня эксплуатируемых нефтегазовых скважин, заполненных любыми типами флюидов
CN112922587B (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
GR01 Patent grant