CN105488349B - 页岩油气储层的微观地质特征确定方法和装置 - Google Patents

页岩油气储层的微观地质特征确定方法和装置 Download PDF

Info

Publication number
CN105488349B
CN105488349B CN201510862762.7A CN201510862762A CN105488349B CN 105488349 B CN105488349 B CN 105488349B CN 201510862762 A CN201510862762 A CN 201510862762A CN 105488349 B CN105488349 B CN 105488349B
Authority
CN
China
Prior art keywords
axis
area
dimensional geologic
image
shale oil
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
CN201510862762.7A
Other languages
English (en)
Other versions
CN105488349A (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 University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
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 University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN201510862762.7A priority Critical patent/CN105488349B/zh
Publication of CN105488349A publication Critical patent/CN105488349A/zh
Application granted granted Critical
Publication of CN105488349B publication Critical patent/CN105488349B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/05—Geographic models
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00—Subject matter not provided for in other main groups of this subclass

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Remote Sensing (AREA)
  • Computer Graphics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

本发明提供一种页岩油气储层的微观地质特征确定方法和装置，通过对页岩油气储层数字化三维地质体图像进行分轴处理后，确定X轴和Y轴所在水平方向上三维地质体图像的单元面积，并确定Z轴方向上三维地质体图像的单元高度；将单元面积和单元高度的乘积作为页岩油气储层的表征单元体积，采用根据页岩油气储层三维空间上水平方向单元面积和垂直方向单元高度相结合的方式，使得所确定的表征单元体积能够准确反映页岩油气储层的微观地质特征，提高了页岩油气储层的微观地质特征计算的准确度，解决了现有技术中所确定的页岩油气储层的微观地质特征参数不准确的技术问题。

Description

页岩油气储层的微观地质特征确定方法和装置
技术领域
本发明涉及油气勘探技术，尤其涉及一种页岩油气储层的微观地质特征确定方法和装置。
背景技术
页岩油气是非常规油气资源的重要组成部分，也是目前油气勘探的重点和热点。页岩油气主要富集在页岩油气储层中，不同的页岩油气储层由于微观地质特征的不同，导致其资源潜力和开发价值差别很大，页岩受沉积环境的影响，微观地质特征在垂向上和水平上存在显著的差异，这种差异在微观上对页岩油气的富集及开采具有十分重要的控制作用，但是这种差异在研究页岩油气储层的开发特征时往往容易被忽略，导致目前很多勘探井的失利。因此需要在页岩油气勘探与开发前准确确定页岩油气储层的微观地质特征。
现有技术中，对页岩油气储层的微观地质特征偏向于基于感性认识的定性描述，忽略了其在水平方向和垂直方向上的差异，往往导致所确定的页岩油气储层的微观地质特征不准确。也就是说，现有技术中，缺乏对页岩油气储层的微观地质特征的定量确定方法，存在所确定的页岩油气储层的微观地质特征准确度不高的技术问题。
发明内容
针对现有技术的上述缺陷，本发明提供一种页岩油气储层的微观地质特征确定方法和装置，用于提高所确定的页岩油气储层的微观地质特征的准确度。
本发明提供一种页岩油气储层的微观地质特征确定方法，包括：
对页岩油气储层数字化三维地质体图像进行分轴处理，其中，Z轴与地表垂直，X轴和Y轴与地表平行；
确定X轴和Y轴所在水平方向上，满足在不同第一矩形区域内，三维地质体图像中各矿物质含量比值之间的差异小于第一预设阈值条件下的最小第一矩形面积，并将最小第一矩形面积作为三维地质体图像的单元面积；
确定Z轴方向上，满足在不同高度下，三维地质体图像中各矿物质含量比值之间的差异小于第二预设阈值条件下的最小高度，并将最小高度作为三维地质体图像的单元高度；
将单元面积和单元高度的乘积作为页岩油气储层的表征单元体积。
在本发明的一实施例中，确定X轴和Y轴所在水平方向上，满足在不同第一矩形区域内，三维地质体图像中各矿物质含量比值之间的差异小于第一预设阈值条件下的最小第一矩形面积，并将最小第一矩形面积作为三维地质体图像的单元面积，包括：
确定X轴方向上，满足在不同第二矩形区域内，三维地质体图像中各矿物质含量比值之间的差异小于第三预设阈值条件下的最小第二矩形面积，并将最小第二矩形面积作为三维地质体图像的X轴单元面积；其中，各第二矩形区域与X轴平行的边长均等于三维地质体图像在X轴方向上的长度；
确定Y轴方向上，满足在不同第三矩形区域内，三维地质体图像中各矿物质含量比值之间的差异小于第四预设阈值条件下的最小第三矩形面积，并将最小第三矩形面积作为三维地质体图像的Y轴单元面积；其中，各第三矩形区域与Y轴平行的边长均等于三维地质体图像在Y轴方向上的长度；
将X轴单元面积和Y轴单元面积的均值作为三维地质体图像的单元面积。
在本发明的一实施例中，三维地质体图像为多个，将单元面积和单元高度的乘积作为页岩油气储层的表征单元体积，包括：
将各三维地质体图像的单元面积和单元高度的乘积的均值作为页岩油气储层的表征单元体积。
本发明还提供一种页岩油气储层的微观地质特征确定装置，包括：
分轴处理模块，用于对页岩油气储层数字化三维地质体图像进行分轴处理，其中，Z轴与地表垂直，X轴和Y轴与地表平行；
面积确定模块，用于确定X轴和Y轴所在水平方向上，满足在不同第一矩形区域内，三维地质体图像中各矿物质含量比值之间的差异小于第一预设阈值条件下的最小第一矩形面积，并将最小第一矩形面积作为三维地质体图像的单元面积；
高度确定模块，用于确定Z轴方向上，满足在不同高度下，三维地质体图像中各矿物质含量比值之间的差异小于第二预设阈值条件下的最小高度，并将最小高度作为三维地质体图像的单元高度；
体积确定模块，用于将单元面积和单元高度的乘积作为页岩油气储层的表征单元体积。
在本发明的一实施例中，面积确定模块，具体用于：
确定X轴方向上，满足在不同第二矩形区域内，三维地质体图像中各矿物质含量比值之间的差异小于第三预设阈值条件下的最小第二矩形面积，并将最小第二矩形面积作为三维地质体图像的X轴单元面积；其中，各第二矩形区域与X轴平行的边长均等于三维地质体图像在X轴方向上的长度；
确定Y轴方向上，满足在不同第三矩形区域内，三维地质体图像中各矿物质含量比值之间的差异小于第四预设阈值条件下的最小第三矩形面积，并将最小第三矩形面积作为三维地质体图像的Y轴单元面积；其中，各第三矩形区域与Y轴平行的边长均等于三维地质体图像在Y轴方向上的长度；
将X轴单元面积和Y轴单元面积的均值作为三维地质体图像的单元面积。
在本发明的一实施例中，三维地质体图像为多个，体积确定模块，具体用于：
将各三维地质体图像的单元面积和单元高度的乘积的均值作为页岩油气储层的表征单元体积。
本发明实施例提供的页岩油气储层的微观地质特征确定方法和装置，对页岩油气储层数字化三维地质体图像进行分轴处理后，确定X轴和Y轴所在水平方向上，满足在不同第一矩形区域内，三维地质体图像中各矿物质含量比值之间的差异小于第一预设阈值条件下的最小第一矩形面积，并将最小第一矩形面积作为三维地质体图像的单元面积；并确定Z轴方向上，满足在不同高度下，三维地质体图像中各矿物质含量比值之间的差异小于第二预设阈值条件下的最小高度，并将最小高度作为三维地质体图像的单元高度；将单元面积和单元高度的乘积作为页岩油气储层的表征单元体积。由于采用了根据页岩油气储层三维空间上水平方向单元面积和垂直方向单元高度相结合的方式，使得所确定的表征单元体积能够准确反映页岩油气储层的微观地质特征，提高了页岩油气储层的微观地质特征计算的准确度，解决了现有技术中所确定的页岩油气储层的微观地质特征参数不准确的技术问题。
附图说明
图1为本发明实施例一提供的页岩油气储层的微观地质特征确定方法的流程示意图；
图2为本发明中进行分轴处理后的三维地质体图像的结构示意图；
图3为本发明实施例二提供的页岩油气储层的微观地质特征确定方法的流程示意图；
图4为本发明提供的页岩油气储层的微观地质特征确定装置的结构示意图。
附图标记说明：
10-分轴处理模块；
20-面积确定模块；
30-高度确定模块；
40-体积确定模块。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
图1为本发明实施例一提供的页岩油气储层的微观地质特征确定方法的流程示意图，本实施例的执行主体可以是页岩油气储层的微观地质特征确定装置，该装置可以集成在计算机等处理设备上。如图1所示，本实施例提供的方法包括：
步骤S101、对页岩油气储层数字化三维地质体图像进行分轴处理。
具体的，对于采集的页岩油气储层样品，可以用扫描电镜获取其三维地质体图像，然后进行数字化图像预处理获取数字化三维地质体图像，该数字化三维地质体图像中包含不同大小的各种矿物质。为了便于说明，以下将数字化三维地质体图像简称为图像。
对于获取到的图像，可以按照样品的三维采集方向进行分轴处理，将垂直地表方向作为Z轴方向，与地表平行方向作为X、Y轴方向。图2为本发明中进行分轴处理后的三维地质体图像的结构示意图，如图2所示，划分出的Z轴与地表垂直，X轴和Y轴与地表平行。需要说明的是，图2中三维地质体图像与X、Y、Z三轴的位置关系只是作为一种示例，X、Y、Z三轴也可以穿过三维地质体图像，具体可根据需要确定。
步骤S102、确定X轴和Y轴所在水平方向上，满足在不同第一矩形区域内，三维地质体图像中各矿物质含量比值之间的差异小于第一预设阈值条件下的最小第一矩形面积，并将最小第一矩形面积作为三维地质体图像的单元面积。
具体的，可以在X轴和Y轴所在水平方向上(即与地表平行的方向上)，任意选取图像的一个截面，在截面内选取一个矩形区域，逐渐扩大该矩形区域的面积，并计算各矩形区域内页岩油气储层所含的各矿物质含量比值；当各矿物质含量比值随着矩形区域的扩大而趋于稳定，不再发生剧烈变化时，将稳定状态下的最小矩形区域面积，即在不同矩形区域内三维地质体图像中各矿物质含量比值之间的差异首次小于第一预设阈值的矩形区域面积，作为三维地质体图像的单元面积。
以图2所示图像为例，假设图像大小为X*Y*Z＝2cm*1cm*1cm，在水平方向上选取图像的一个截面，在截面内选取一个1mm*1mm的矩形区域，计算该矩形区域内页岩油气储层的各矿物质含量比值；然后逐渐扩大该矩形区域的面积，例如每次任意选取矩形区域一个边长扩大1mm，或者按照其他方式扩大矩形区域，计算每次扩大后的矩形区域内页岩油气储层的各矿物质含量比值；当各矿物质含量比值随着矩形区域的扩大而趋于稳定，不再发生剧烈变化时，例如矩形区域在大小为16mm*7mm时与在大小为16mm*8mm以及更大时，各矿物质含量比值相似，此时可将稳定状态下的最小矩形区域面积，即16mm*7mm的矩形区域面积作为三维地质体图像的单元面积。
对于第一预设阈值，可以根据经验值选取，在该阈值下，各矿物质含量比值差异较小，趋于稳定。另外，在计算矩形区域内所述三维地质体图像中各矿物质含量比值时，可以根据各矿物质在矩形区域内所占的面积进行计算。
步骤S103、确定Z轴方向上，满足在不同高度下，三维地质体图像中各矿物质含量比值之间的差异小于第二预设阈值条件下的最小高度，并将最小高度作为三维地质体图像的单元高度。
Z轴方向上的图像与地层垂直，该方向上沉积时间较长，矿物组成特征与水平方向不同，因此有必要对Z轴方向上的矿物质含量进行计算。具体的，可以在Z轴方向上(即与地表垂直的方向上)，任意选取图像的一条直线，逐渐扩大该直线的高度，并计算各直线上页岩油气储层的各矿物质含量比值；当各矿物质含量比值随着矩形区域的扩大而趋于稳定，不再发生剧烈变化时，将稳定状态下的最小直线高度，即在不同直线上三维地质体图像中各矿物质含量比值之间的差异首次小于第二预设阈值的直线高度，作为三维地质体图像的单元面积。
继续以图2所示图像为例，图像大小为X*Y*Z＝2cm*1cm*1cm，在Z轴方向上选取图像的一条直线，例如距离图像底端高度为1mm的直线，计算该直线上页岩油气储层的各矿物质含量比值；然后逐渐增加该直线的高度，例如每次将直线向上增高1mm，计算每次增高后的直线上页岩油气储层的各矿物质含量比值；当各矿物质含量比值随着直线高度的增加而趋于稳定，不再发生剧烈变化时，例如直线高度在7mm时与在8mm以及更大时，各矿物质含量比值相似，此时可将稳定状态下的最小直线高度，即7mm的高度作为三维地质体图像的单元高度。
对于第二预设阈值，可以根据经验值选取，在该阈值下，各矿物质含量比值差异较小，趋于稳定。另外，在计算直线上三维地质体图像中各矿物质含量比值时，可以根据各矿物质在直线上所占的长度进行计算。
步骤S104、将单元面积和单元高度的乘积作为页岩油气储层的表征单元体积。
根据在水平方向得到的单元面积S和垂直方向上得到的单元高度H，两者的乘积即为该页岩油气储层的表征单元体积V，即：
V＝S*H (1)
该表征单元体积V即可反映页岩油气储层的微观地质特征。
需要说明的是，步骤S102和步骤S103之间没有严格的时序关系，步骤S102也可以在步骤S103之后执行，还可以与步骤S103同时执行，具体执行顺序本实施例不做限制。
本实施例提供的页岩油气储层的微观地质特征确定方法，对页岩油气储层数字化三维地质体图像进行分轴处理后，确定X轴和Y轴所在水平方向上，满足在不同第一矩形区域内，三维地质体图像中各矿物质含量比值之间的差异小于第一预设阈值条件下的最小第一矩形面积，并将最小第一矩形面积作为三维地质体图像的单元面积；并确定Z轴方向上，满足在不同高度下，三维地质体图像中各矿物质含量比值之间的差异小于第二预设阈值条件下的最小高度，并将最小高度作为三维地质体图像的单元高度；将单元面积和单元高度的乘积作为页岩油气储层的表征单元体积。由于采用了根据页岩油气储层三维空间上水平方向单元面积和垂直方向单元高度相结合的方式，使得所确定的表征单元体积能够准确反映页岩油气储层的微观地质特征，提高了页岩油气储层的微观地质特征计算的准确度，解决了现有技术中所确定的页岩油气储层的微观地质特征参数不准确的技术问题。
图3为本发明实施例二提供的页岩油气储层的微观地质特征确定方法的流程示意图，本实施例是对上述图1所示实施例中步骤S102的一种具体实现方式的详细说明。在上述图1所示实施例的基础上，如图3所示，本实施例中，步骤S102确定X轴和Y轴所在水平方向上，满足在不同矩形区域内，三维地质体图像中各矿物质含量比值之间的差异小于第一预设阈值条件下的最小矩形面积，并将最小矩形面积作为水平方向上的单元面积，具体包括：
步骤S201、确定X轴方向上，满足在不同第二矩形区域内，三维地质体图像中各矿物质含量比值之间的差异小于第三预设阈值条件下的最小第二矩形面积，并将最小第二矩形面积作为三维地质体图像的X轴单元面积。
具体的，可以在水平方向上任意选取图像的一个截面，在截面内选取一个矩形区域，该矩形区域与X轴平行的边长固定，等于三维地质体图像在X轴方向上的长度；逐渐扩大该矩形区域的面积，即扩大矩形区域与Y轴平行的边长，并计算各矩形区域内页岩油气储层的各矿物质含量比值；当各矿物质含量比值随着矩形区域的扩大而趋于稳定，不再发生剧烈变化时，将稳定状态下的最小矩形区域面积，即在不同矩形区域内三维地质体图像中各矿物质含量比值之间的差异首次小于第三预设阈值的矩形区域面积，作为三维地质体图像的X轴单元面积。
继续以图2所示图像为例，图像大小为X*Y*Z＝2cm*1cm*1cm，在水平方向上选取图像的一个截面，在截面内靠近原点处选取一个2cm*1mm的矩形区域，计算该矩形区域内页岩油气储层的各矿物质含量比值；然后逐渐扩大该矩形区域的另一个边长，例如每次扩大1mm，计算每次扩大后的矩形区域内页岩油气储层的各矿物质含量比值；当各矿物质含量比值随着矩形区域的扩大而趋于稳定，不再发生剧烈变化时，例如矩形区域在大小为2cm*7mm时与在大小为2cm*8mm以及更大时，各矿物质含量比值相似，此时可将稳定状态下的最小矩形区域面积，即2cm*7mm的矩形区域面积作为三维地质体图像的X轴单元面积。
步骤S202、确定Y轴方向上，满足在不同第三矩形区域内，三维地质体图像中各矿物质含量比值之间的差异小于第四预设阈值条件下的最小第三矩形面积，并将最小第三矩形面积作为三维地质体图像的Y轴单元面积。
在水平方向上任意选取图像的一个截面，确定X轴单元面积后，可以继续确定该截面中的Y轴单元面积。Y轴单元面积的计算原理与X轴单元面积的计算原理类似，与计算X轴单元面积不同之处在于，选取的矩形区域与Y轴平行的边长固定，等于三维地质体图像在Y轴方向上的长度；在逐渐扩大矩形区域的面积时，扩大的是矩形区域与X轴平行的边长。具体获取Y轴单元面积的其他过程与获取X轴单元面积时类似，此处不再赘述。
步骤S203、将X轴单元面积和Y轴单元面积的均值作为三维地质体图像的单元面积。
X、Y轴方向上的图像与地层平行，两个方向上的矿物均于同期发生沉积，具有随机性，因此，该类图像的单元面积为两个方向上求取的单元面积的平均值。即，三维地质体图像的单元面积S为X轴单元面积Sx和Y轴单元面积Sy的均值，用公式表示为：
S＝(Sx+Sy)/2 (2)
需要说明的是，步骤S201和步骤S202之间没有严格的时序关系，步骤S201也可以在步骤S202之后执行，还可以与步骤S202同时执行，具体执行顺序本实施例不做限制。
对于上述第三预设阈值和第四预设阈值，都可以根据经验值选取，在该阈值下，各矿物质含量比值差异较小，趋于稳定。第三预设阈值和第四预设阈值可以相等，也可以不等。
作为本发明一种可选的实施方式，三维地质体图像可以为多个，则步骤S104中将单元面积和单元高度的乘积作为页岩油气储层的表征单元体积，包括：将各三维地质体图像的单元面积和单元高度的乘积的均值作为页岩油气储层的表征单元体积。
具体的，为了提高所确定的页岩油气储层的微观地质特征的准确度，对于页岩油气储层，可以采集多份样品，对每份样品采用上述方法获取其对应的三维地质体图像的单元面积和单元高度；将各三维地质体图像的单元面积和单元高度的乘积作为各三维地质体图像的表征单元体积，然后取各三维地质体图像的表征单元体积的均值作为页岩油气储层的表征单元体积。
当然，也可以对各三维地质体图像的表征单元体积进行其他处理获取页岩油气储层的表征单元体积，例如：观察各三维地质体图像的表征单元体积的分布特征，对取值在某个区间段内的三维地质体图像的表征单元体积进行加权平均或者取最大或最小值，作为页岩油气储层的表征单元体积。
本实施例提供的页岩油气储层的微观地质特征确定方法，通过获取三维地质体图像的X轴单元面积和Y轴单元面积，来确定三维地质体图像的单元体积，能够减小计算误差，进一步提高所确定的页岩油气储层的微观地质特征的准确度。
图4为本发明提供的页岩油气储层的微观地质特征确定装置的结构示意图，该装置可以集成在计算机等处理设备中，本实施例提供的装置包括：分轴处理模块10、面积确定模块20、高度确定模块30和体积确定模块40，其中，
分轴处理模块10，用于对页岩油气储层数字化三维地质体图像进行分轴处理，其中，Z轴与地表垂直，X轴和Y轴与地表平行；
面积确定模块20，用于确定X轴和Y轴所在水平方向上，满足在不同第一矩形区域内，三维地质体图像中各矿物质含量比值之间的差异小于第一预设阈值条件下的最小第一矩形面积，并将最小第一矩形面积作为三维地质体图像的单元面积；
高度确定模块30，用于确定Z轴方向上，满足在不同高度下，三维地质体图像中各矿物质含量比值之间的差异小于第二预设阈值条件下的最小高度，并将最小高度作为三维地质体图像的单元高度；
体积确定模块40，用于将单元面积和单元高度的乘积作为页岩油气储层的表征单元体积。
本实施例提供的装置可以执行上述方法实施例，具体实现原理和技术效果类似，此处不再赘述。
在本发明的一实施例中，面积确定模块20，具体用于：
确定X轴方向上，满足在不同第二矩形区域内，三维地质体图像中各矿物质含量比值之间的差异小于第三预设阈值条件下的最小第二矩形面积，并将最小第二矩形面积作为三维地质体图像的X轴单元面积；其中，各第二矩形区域与X轴平行的边长均等于三维地质体图像在X轴方向上的长度；
确定Y轴方向上，满足在不同第三矩形区域内，三维地质体图像中各矿物质含量比值之间的差异小于第四预设阈值条件下的最小第三矩形面积，并将最小第三矩形面积作为三维地质体图像的Y轴单元面积；其中，各第三矩形区域与Y轴平行的边长均等于三维地质体图像在Y轴方向上的长度；
将X轴单元面积和Y轴单元面积的均值作为三维地质体图像的单元面积。
在本发明的另一实施例中，三维地质体图像为多个，体积确定模块40，具体用于：将各三维地质体图像的单元面积和单元高度的乘积的均值作为页岩油气储层的表征单元体积。
本实施例提供的装置可以执行上述方法实施例，具体实现原理和技术效果类似，此处不再赘述。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。

Claims (4)

1.一种页岩油气储层的微观地质特征确定方法，其特征在于，包括：
对页岩油气储层数字化三维地质体图像进行分轴处理，其中，Z轴与地表垂直，X轴和Y轴与地表平行；
确定X轴和Y轴所在水平方向上，满足在不同第一矩形区域内，所述三维地质体图像中各矿物质含量比值之间的差异小于第一预设阈值条件下的最小第一矩形面积，并将所述最小第一矩形面积作为所述三维地质体图像的单元面积；
确定Z轴方向上，满足在不同高度下，所述三维地质体图像中各矿物质含量比值之间的差异小于第二预设阈值条件下的最小高度，并将所述最小高度作为所述三维地质体图像的单元高度；
将所述单元面积和所述单元高度的乘积作为所述页岩油气储层的表征单元体积；
所述确定X轴和Y轴所在水平方向上，满足在不同第一矩形区域内，所述三维地质体图像中各矿物质含量比值之间的差异小于第一预设阈值条件下的最小第一矩形面积，并将所述最小第一矩形面积作为所述三维地质体图像的单元面积，包括：
确定X轴方向上，满足在不同第二矩形区域内，所述三维地质体图像中各矿物质含量比值之间的差异小于第三预设阈值条件下的最小第二矩形面积，并将所述最小第二矩形面积作为所述三维地质体图像的X轴单元面积；其中，各所述第二矩形区域与X轴平行的边长均等于所述三维地质体图像在X轴方向上的长度；
确定Y轴方向上，满足在不同第三矩形区域内，所述三维地质体图像中各矿物质含量比值之间的差异小于第四预设阈值条件下的最小第三矩形面积，并将所述最小第三矩形面积作为所述三维地质体图像的Y轴单元面积；其中，各所述第三矩形区域与Y轴平行的边长均等于所述三维地质体图像在Y轴方向上的长度；
将所述X轴单元面积和所述Y轴单元面积的均值作为所述三维地质体图像的单元面积。
2.根据权利要求1所述的方法，其特征在于，所述三维地质体图像为多个，所述将所述单元面积和所述单元高度的乘积作为所述页岩油气储层的表征单元体积，包括：
将各所述三维地质体图像的单元面积和单元高度的乘积的均值作为所述页岩油气储层的表征单元体积。
3.一种页岩油气储层的微观地质特征确定装置，其特征在于，包括：
分轴处理模块，用于对页岩油气储层数字化三维地质体图像进行分轴处理，其中，Z轴与地表垂直，X轴和Y轴与地表平行；
面积确定模块，用于确定X轴和Y轴所在水平方向上，满足在不同第一矩形区域内，所述三维地质体图像中各矿物质含量比值之间的差异小于第一预设阈值条件下的最小第一矩形面积，并将所述最小第一矩形面积作为所述三维地质体图像的单元面积；
高度确定模块，用于确定Z轴方向上，满足在不同高度下，所述三维地质体图像中各矿物质含量比值之间的差异小于第二预设阈值条件下的最小高度，并将所述最小高度作为所述三维地质体图像的单元高度；
体积确定模块，用于将所述单元面积和所述单元高度的乘积作为所述页岩油气储层的表征单元体积；
所述面积确定模块，具体用于：
确定X轴方向上，满足在不同第二矩形区域内，所述三维地质体图像中各矿物质含量比值之间的差异小于第三预设阈值条件下的最小第二矩形面积，并将所述最小第二矩形面积作为所述三维地质体图像的X轴单元面积；其中，各所述第二矩形区域与X轴平行的边长均等于所述三维地质体图像在X轴方向上的长度；
确定Y轴方向上，满足在不同第三矩形区域内，所述三维地质体图像中各矿物质含量比值之间的差异小于第四预设阈值条件下的最小第三矩形面积，并将所述最小第三矩形面积作为所述三维地质体图像的Y轴单元面积；其中，各所述第三矩形区域与Y轴平行的边长均等于所述三维地质体图像在Y轴方向上的长度；
将所述X轴单元面积和所述Y轴单元面积的均值作为所述三维地质体图像的单元面积。
4.根据权利要求3所述的装置，其特征在于，所述三维地质体图像为多个，所述体积确定模块，具体用于：
将各所述三维地质体图像的单元面积和单元高度的乘积的均值作为所述页岩油气储层的表征单元体积。
CN201510862762.7A 2015-12-01 2015-12-01 页岩油气储层的微观地质特征确定方法和装置 Active CN105488349B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510862762.7A CN105488349B (zh) 2015-12-01 2015-12-01 页岩油气储层的微观地质特征确定方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510862762.7A CN105488349B (zh) 2015-12-01 2015-12-01 页岩油气储层的微观地质特征确定方法和装置

Publications (2)

Publication Number Publication Date
CN105488349A CN105488349A (zh) 2016-04-13
CN105488349B true CN105488349B (zh) 2018-06-08

Family

ID=55675323

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510862762.7A Active CN105488349B (zh) 2015-12-01 2015-12-01 页岩油气储层的微观地质特征确定方法和装置

Country Status (1)

Country Link
CN (1) CN105488349B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008041885A1 (ru) * 2006-10-06 2008-04-10 Zakrytoe Aktsionernoe Obschestvo 'modelirovanie I Monitoring Geologicheskikh Obiektov Im. V.A. Dvurechenskogo' Способ размещения поисковых, разведочных и эксплуатационных скважин на месторождениях нефти и газа на основе трехмерной геологической модели
CN101661514A (zh) * 2008-05-21 2010-03-03 中国石化股份胜利油田分公司地质科学研究院 一种油藏黑油模型数值模拟系统
CN102612682A (zh) * 2009-11-12 2012-07-25 埃克森美孚上游研究公司 用于储层建模和模拟的方法和设备
CN104573339A (zh) * 2014-12-24 2015-04-29 中国石油大学(北京) 页岩气储层的地质参数确定方法和装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008041885A1 (ru) * 2006-10-06 2008-04-10 Zakrytoe Aktsionernoe Obschestvo 'modelirovanie I Monitoring Geologicheskikh Obiektov Im. V.A. Dvurechenskogo' Способ размещения поисковых, разведочных и эксплуатационных скважин на месторождениях нефти и газа на основе трехмерной геологической модели
CN101661514A (zh) * 2008-05-21 2010-03-03 中国石化股份胜利油田分公司地质科学研究院 一种油藏黑油模型数值模拟系统
CN102612682A (zh) * 2009-11-12 2012-07-25 埃克森美孚上游研究公司 用于储层建模和模拟的方法和设备
CN104573339A (zh) * 2014-12-24 2015-04-29 中国石油大学(北京) 页岩气储层的地质参数确定方法和装置

Also Published As

Publication number Publication date
CN105488349A (zh) 2016-04-13

Similar Documents

Publication Publication Date Title
Bazaikin et al. Effect of CT image size and resolution on the accuracy of rock property estimates
Weinberg et al. What controls gold distribution in Archean terranes?
Andrew et al. Pore-scale contact angle measurements at reservoir conditions using X-ray microtomography
Hajek et al. Significance of channel-belt clustering in alluvial basins
Chamberlin et al. Interpreting paleo-avulsion dynamics from multistory sand bodies
Adams et al. Digital characterization of thrombolite-stromatolite reef distribution in a carbonate ramp system (terminal Proterozoic, Nama Group, Namibia)
Lia et al. Uncertainties in reservoir production forecasts
East et al. Fluvial–eolian interactions in sediment routing and sedimentary signal buffering: an example from the Indus Basin and Thar Desert
Taylor et al. Estimation of in situ permeability of deformation bands in porous sandstone, Valley of Fire, Nevada
CN107367762B (zh) 一种确定储层参数的方法及装置
CN107655908A (zh) 一种构建数字岩心的方法及装置
Zuchuat et al. Composite and diachronous stratigraphic surfaces in low-gradient, transitional settings: The J-3 “unconformity” and the Curtis Formation, east-central Utah, USA
CN111983671B (zh) 基于微古地貌恢复的浅水湖盆储层预测方法和装置
Osman et al. Is cell-to-cell scale variability necessary in reservoir models?
CN104360412A (zh) 致密深盆气成藏预测方法和装置
Fang et al. Fine description of structure and sedimentary microfacies of Li32 block of Lijin oilfield, Dongying depression, China
CN105488349B (zh) 页岩油气储层的微观地质特征确定方法和装置
Orellana et al. Influence of variograms in 3D reservoir-modeling outcomes: An example
CN107368688B (zh) 海相单砂体的确定方法和装置
Le Roux et al. Determining the Neogene behavior of the Nazca plate by geohistory analysis
Svancara Approximate method for direct interpretation of gravity anomalies caused by surface three-dimensional geologic structures
CN111208558A (zh) 超深低幅度三维地质构造的建立方法及装置
Chang et al. Facies parameterization and estimation for complex reservoirs-the Brugge field
Székely et al. Automated recognition of quasi‐planar ignimbrite sheets as paleosurfaces via robust segmentation of digital elevation models: an example from the Central Andes
Shao et al. 3D Geological Modeling under Extremely Complex Geological Conditions.

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Tang Xianglu

Inventor after: Jiang Zhenxue

Inventor after: Feng Jie

Inventor after: Li Zhuo

Inventor after: Yuan Yuan

Inventor after: Huang Hexin

Inventor after: Wang Pengfei

Inventor after: Wen Nuan

Inventor before: Feng Jie

Inventor before: Jiang Zhenxue

Inventor before: Tang Xianglu

Inventor before: Li Zhuo

Inventor before: Yuan Yuan

Inventor before: Huang Hexin

Inventor before: Wang Pengfei

Inventor before: Wen Nuan

GR01 Patent grant
GR01 Patent grant