[go: up one dir, main page]

CN110966000A - 一种用于表征砂砾岩储层层状指数的方法及系统 - Google Patents

一种用于表征砂砾岩储层层状指数的方法及系统 Download PDF

Info

Publication number
CN110966000A
CN110966000A CN202010000508.7A CN202010000508A CN110966000A CN 110966000 A CN110966000 A CN 110966000A CN 202010000508 A CN202010000508 A CN 202010000508A CN 110966000 A CN110966000 A CN 110966000A
Authority
CN
China
Prior art keywords
index
reservoir
conductivity
column
row
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
CN202010000508.7A
Other languages
English (en)
Other versions
CN110966000B (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.)
North China Branch Of Cnpc Logging Co ltd
China National Petroleum Corp
China Petroleum Logging Co Ltd
Original Assignee
North China Branch Of Cnpc 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 North China Branch Of Cnpc Logging Co Ltd filed Critical North China Branch Of Cnpc Logging Co Ltd
Priority to CN202010000508.7A priority Critical patent/CN110966000B/zh
Publication of CN110966000A publication Critical patent/CN110966000A/zh
Application granted granted Critical
Publication of CN110966000B publication Critical patent/CN110966000B/zh
Expired - Fee Related 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
    • E21B49/00Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling 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

Landscapes

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

Abstract

本发明公开了一种用于表征砂砾岩储层层状指数的方法及系统,包括:利用电成像测井技术,采集砂砾岩储层图像;按照预设的行间隔和列间隔,将当前图像进行网格化划分,由每行的行电导率值的标准方差,得到当前储层横向连续性指数,以及由每列的列电导率值的标准方差,得到相应的纵向连续性指数;根据横向连续性指数和纵向连续性指数,计算当前储层的层状指数,用以定量表征当前储层成层状性的强弱。本发明实现了储层层状指数的量化表征,并确定了有效储层的层状指数界限,为基于电成像的有效储层识别提供了有效手段。

Description

一种用于表征砂砾岩储层层状指数的方法及系统
技术领域
本发明涉及石油测井工程领域,尤其是涉及一种用于表征砂砾岩储层层状指数的方法及系统。
背景技术
近些年,随着我国油气勘探开发领域的不断拓展,砂砾岩油气藏成为了深入勘探的对象之一。电成像测井能够提供井筒周围地层直观的、清晰的、分辨率高的图像,其电成像图像包含了地层中丰富的岩性、岩石结构等地质信息。在砂砾岩储层段,电成像图像能清晰地反映储层的层状结构特征。
目前,通过对研究区大量实例的定性分析,发现砂砾岩储层的成层状性与储层产能密切相关,储层成层状性的强弱是判断储层有效性的关键指标。前人在利用电成像图像评价砂砾岩储层结构方面做了大量研究,但是对砂砾岩成层状性的储层结构研究较少,且都停留在对储层成层状性的定性描述上,缺乏对储层成层状性进行量化表征的有效手段,也就无法基于储层成层状性来建立有效储层的定量判识标准。
发明内容
为了解决上述技术问题,本发明实施例提供了一种用于表征砂砾岩储层层状指数的方法,包括:步骤一、利用电成像测井技术,采集砂砾岩储层图像;步骤二、按照预设的行间隔和列间隔,将所述图像进行网格化划分,由每行的行电导率值的标准方差,得到当前储层横向连续性指数,以及由每列的列电导率值的标准方差,得到相应的纵向连续性指数;步骤三、根据所述横向连续性指数和所述纵向连续性指数,计算当前储层的层状指数,用以定量表征当前储层成层状性的强弱。
优选地,所述方法,还包括:根据待研究区块内不同储层对应的所述层状指数,并结合所述不同储层的产能信息,确定用于识别所述待研究区块内储层是否有效的层状指数标准,以利用所述层状指数标准判断所述待研究区块内的储层有效性。
优选地,在所述步骤二中,构建关于所述砂砾岩储层图像的电导率矩阵;计算所述图像中每行的行电导率标准方差,并计算各行的行电导率标准方差的平均值,得到所述横向连续性指数;计算所述图像中每列的列电导率标准方差,并计算各列的列电导率标准方差的平均值,得到所述纵向连续性指数。
优选地,将所述纵向连续性指数与所述横向连续性指数进行除法运算,得到所述层状指数。
优选地,所述方法,还包括:按照预设窗口,对经过网格化划分处理的所述图像进行扫描,构建关于每个扫描窗口的电导率矩阵;通过执行所述步骤二和所述步骤三,计算每个扫描窗口对应的所述横向连续性指数、所述纵向连续性指数和相应的所述层状指数;将所有所述层状指数按照储层深度进行连接,生成针对当前储层的层状指数曲线。
另一方面,提供了一种用于表征砂砾岩储层层状指数的系统,包括:图像获取模块,其配置为利用电成像测井技术,采集砂砾岩储层图像;连续性指数生成模块,其配置为按照预设的行间隔和列间隔,将所述图像进行网格化划分,由每行的行电导率值的标准方差,得到当前储层横向连续性指数,以及由每列的列电导率值的标准方差,得到相应的纵向连续性指数;层状指数生成模块,其配置为根据所述横向连续性指数和所述纵向连续性指数,计算当前储层的层状指数,用以定量表征当前储层成层状性的强弱。
优选地,所述系统,还包括:标准生成模块,其配置为根据待研究区块内不同储层对应的所述层状指数,并结合所述不同储层的产能信息,确定用于识别所述待研究区块内储层是否有效的层状指数标准,以利用所述层状指数标准判断所述待研究区块内的储层有效性。
优选地,所述连续性指数生成模块,包括:电导率矩阵构建单元,其配置为构建关于所述砂砾岩储层图像的电导率矩阵;横向连续性指数计算单元,其配置为计算所述图像中每行的行电导率标准方差,并计算各行的行电导率标准方差的平均值,得到所述横向连续性指数;纵向连续性指数计算单元,其配置为计算所述图像中每列的列电导率标准方差,并计算各列的列电导率标准方差的平均值,得到所述纵向连续性指数。
优选地,所述层状指数生成模块,其进一步配置为将所述纵向连续性指数与所述横向连续性指数进行除法运算,得到所述层状指数。
优选地,所述系统还包括曲线绘制模块,其中,所述曲线绘制模块具备:窗口扫描单元,其配置为按照预设窗口,对经过网格化划分处理的所述图像进行扫描,构建关于每个扫描窗口的电导率矩阵;窗口层状指数计算单元,其配置为执行所述连续性指数生成模块和所述层状指数生成模块,计算每个扫描窗口对应的所述横向连续性指数、所述纵向连续性指数和相应的所述层状指数;层状指数曲线生成单元,其配置为将所有所述层状指数按照储层深度进行连接,生成针对当前储层的层状指数曲线。
与现有技术相比,上述方案中的一个或多个实施例可以具有如下优点或有益效果:
本发明提出了一种用于表征砂砾岩储层层状指数的方法及系统。该方法及系统基于电成像测井技术,利用储层电导率图像,计算横向连续性、以及纵向连续性,进一步得到表征储层成层状性强弱的量化指标(层状指数)。另外,利用构建储层层状指数这一量化指标的计算过程,结合待研究区块内不同储层样本的产能情况,确定用于评价当前待研究区块内储层有效性的关键层状指数标准,用以自动判断并识别出待研究区块内储层的有效性。本发明实现了储层成层状性的量化表征,并确定了有效储层的层状指数界限,为基于电成像的有效储层识别提供了有效手段,有效储层识别准确率达84.2%。
本发明的其它特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例共同用于解释本发明,并不构成对本发明的限制。在附图中:
图1为本申请实施例的用于表征砂砾岩储层层状指数的方法的步骤图。
图2为本申请实施例的用于表征砂砾岩储层层状指数的方法的具体流程图。
图3为本申请实施例的用于表征砂砾岩储层层状指数的方法中生成层状指数曲线的具体流程图。
图4为本申请实施例的用于表征砂砾岩储层层状指数的方法中砂砾岩储层图像的两个示例。
图5为本申请实施例的用于表征砂砾岩储层层状指数的方法中横向连续性指数和纵向连续性指数的计算示意图。
图6为本申请实施例的用于表征砂砾岩储层层状指数的方法中层状指数计算实例示意图。
图7为本申请实施例的用于表征砂砾岩储层层状指数的方法中确定层状指数标准的实例示意图。
图8为本申请实施例的用于表征砂砾岩储层层状指数的系统的模块框图。
具体实施方式
以下将结合附图及实施例来详细说明本发明的实施方式,借此对本发明如何应用技术手段来解决技术问题,并达成技术效果的实现过程能充分理解并据以实施。需要说明的是,只要不构成冲突,本发明中的各个实施例以及各实施例中的各个特征可以相互结合,所形成的技术方案均在本发明的保护范围之内。
近些年,随着我国油气勘探开发领域的不断拓展,砂砾岩油气藏成为了深入勘探的对象之一。电成像测井能够提供井筒周围地层直观的、清晰的、分辨率高的图像,其电成像图像包含了地层中丰富的岩性、岩石结构等地质信息。在砂砾岩储层段,电成像图像能清晰地反映储层的层状结构特征。
目前,通过对研究区大量实例的定性分析,发现砂砾岩储层的成层状性与储层产能密切相关,储层成层状性的强弱是判断储层有效性的关键指标。前人在利用电成像图像评价砂砾岩储层结构方面做了大量研究,但是对砂砾岩成层状性的储层结构研究较少,且都停留在对储层成层状性的定性描述上,缺乏对储层成层状性进行量化表征的有效手段,也就无法基于储层成层状性来建立有效储层的定量判识标准。
为了解决上述现有技术中的问题,本发明提出了一种用于表征砂砾岩储层层状指数的方法及系统。该方法及系统基于针对砂砾岩储层的电成像测井资料,通过量化表征储层图像在横向和纵向上的连续性,在此基础上构建出储层的层状指数及相应的层状指数连续曲线,实现了储层成层状性强弱的量化表征方法,以此来预测或评价砂砾岩储层的有效性。另外,本发明还能够利用上述层状指数量化表征方法,在待研究区域内,建立起基于层状指数的有效储层判识标准(有效储层识别阈值参数)。这样,本发明不仅满足了砂砾岩储层图像的成层状性量化表征的迫切需求,也形成了一种储层有效性判识的有效手段。
图1为本申请实施例的用于表征砂砾岩储层层状指数的方法的步骤图。图2为本申请实施例的用于表征砂砾岩储层层状指数的方法的具体流程图。结合图1和图2,步骤S110利用电成像测井技术,采集砂砾岩储层图像。在步骤S110中,对待研究区域内的砂砾岩储层进行电成像测井施工,获取能够反映井筒周围地层岩性、岩石结构等地质信息的高分辨率静态图像或动态图像。图4为本申请实施例的用于表征砂砾岩储层层状指数的方法中砂砾岩储层图像的两个示例。图4展示了两例砂砾岩储层图像,图4a展示了成层状性较弱的砂砾岩储层的图像,左图为静态图像,右图为动态图像;图4b展示了成层状性较强的砂砾岩储层的图像,左图为静态图像,右图为动态图像。需要说明的是,在利用电成像测井技术采集到的砂砾岩储层图像中,能够计算出图像中每个像素点的地层电导率值,基于此,利用后续的步骤S120和步骤S130来对当前储层的成层状性进行量化评价。
而后,步骤S120按照预设的行间隔和列间隔,将当前获得的砂砾岩储层图像进行网格化划分,由每行(网格)的行电导率值的标准方差,得到当前储层的横向连续性指数,以及由每列(网格)的列电导率值的标准方差,得到当前储层的纵向连续性指数。具体地,首先,步骤S1201按照预设的行间隔和列间隔,将当前获得的砂砾岩储层静态图像或动态图像进行网格化划分,从而将当前砂砾岩储层图像划分为若干个网格,每个网格的横向宽度与列间隔的长度相符、并且每个网格的纵向宽度与行间隔的长度相符,从而进入到步骤S1202中。在本发明实施例中,对行间隔和列间隔的长度不作具体限定,本领域技术人员可根据实际的层状指数计算精度进行设置,每个网格可以包括由多个像素点形成的像素块构成,也可以为单独一个像素点。
步骤S1202基于经过网格化划分处理的图像,构建关于当前砂砾岩储层静态图像或动态图像的电导率矩阵。其中,利用如下矩阵表达式表示电导率矩阵:
Figure BDA0002353089880000051
在上述矩阵表达式中,C表示电导率矩阵,m表示电导率矩阵C的行数,n表示电导率矩阵C的列数。其中,电导率矩阵中的每个矩阵元素代表一个网格,每个矩阵元素的值为该网格的电导率值。若一个网格中仅含有一个像素点,那么该网格的电导率值即为该像素点的电导率值;若网格为像素块,那么该网格的电导率值优选为当前块内各个像素点的电导率的平均值。
在完成电导率矩阵构建后,进入到步骤S1203中。步骤S1203沿着地层水平方向,计算当前砂砾岩储层静态图像或动态图像中每行网格的行电导率标准方差,并计算各行的行电导率标准方差的平均值,利用该平均值表示当前储层在横向上的连续性指数,即当前储层横向连续性指数。图5为本申请实施例的用于表征砂砾岩储层层状指数的方法中横向连续性指数和纵向连续性指数的计算示意图。如图5所示,右向箭头所示的方向即为横向连续性表征的方向。
具体地,首先,对每行中所有网格的电导率值进行平均值运算,利用行电导率计算式,计算图像中每行的行电导率;而后,根据每行的行电导率,利用行电导率方差计算式,计算相应行的行电导率的标准方差;最后,根据每行的行电导率标准方差值,利用横向连续性指数计算式,得到关于当前图像反映出的当前储层在横向上的电导率连续性特征,即横向连续性指数。其中,上述行电导率计算式利用下述表达式表示:
Figure BDA0002353089880000061
式(1)中,i表示行序号,j表示列序号,MEAN_H(i)表示第i行的行电导率,C(i,j)表示第i行第j列网格的电导率值。进一步,行电导率方差计算式利用如下表达式表示:
Figure BDA0002353089880000062
式(2)中,STDD_H(i)表示第i行的行电导率标准方差。进一步,横向连续性指数计算式利用如下表达式表示:
Figure BDA0002353089880000063
式(3)中,CONT_H表示当前图像的横向连续性指数。图像的横向连续性指数表示储层电导率在水平横向上的连续性,横向连续性指数值越大,则表明储层在横向上连续性越差,也就是说,在沿着地层水平方向上,储层连续性越弱;反之,横向连续性指数值越小,则表明储层在横向上连续性越好,也就是说,在沿着地层水平方向上,储层连续性越强。
在完成横向连续性指数计算后,进入到步骤S1204中。步骤S1204沿着地层水平的垂直方向,计算当前砂砾岩储层静态图像或动态图像中每列网格的列电导率标准方差,并计算各列的列电导率标准方差的平均值,利用该平均值表示当前储层纵向上的连续性,即当前储层纵向连续性指数。如图5所示,指向向上的箭头所示的方向即为纵向连续性表征的方向。
具体地,首先,对每列中所有网格的电导率值进行平均值运算,利用列电导率计算式,计算图像中每列的列电导率;而后,根据每列的列电导率,利用列电导率方差计算式,计算相应列的列电导率的标准方差;最后,根据每列的列电导率标准方差值,利用纵向连续性指数计算式,得到关于当前图像反映出的当前储层在纵向上的电导率连续性特征,即纵向连续性指数。其中,上述列电导率计算式利用下述表达式表示:
Figure BDA0002353089880000071
式(4)中,MEAN_V(i)表示第j列的列电导率值。进一步,列电导率方差计算式利用如下表达式表示:
Figure BDA0002353089880000072
式(5)中,STDD_V(j)表示第j列的列电导率标准方差。进一步,纵向连续性指数计算式利用如下表达式表示:
Figure BDA0002353089880000073
式(6)中,CONT_V表示当前图像的纵向连续性指数。图像的纵向连续性指数表示储层电导率在垂直、纵向上的连续性,纵向连续性指数值越大,则表明储层在纵向上连续性越差,也就是说,在地层纵向方向上,储层成层状性可能越强;反之,纵向连续性指数值越小,则表明储层在纵向上连续性越连续,也就是说,在地层纵向方向上,储层成层状性可能越弱。
而后,在完成纵向连续性指数的计算后,上述步骤S120结束,进入到步骤S130中。步骤S130根据上述步骤S120得到的横向连续性指数和纵向连续性指数,计算当前砂砾岩储层的层状指数,以利用该层状指数来定量表征当前储层成层状性的强弱。进一步,步骤S1301利用层状指数计算式,将上述纵向连续性指数与上述横向连续性指数进行除法运算,用二者的比值计算结果来表示量化的层状指数。其中,层状指数计算式利用如下表达式表示:
Figure BDA0002353089880000081
式(7)中,LI表示层状指数。在实际应用过程中,层状指数数值越小,则表明储层的成层状性越弱,当前储层越接近于成层状性弱的储层;反之,层状指数数值越大,则表明储层的成层状性越强,当前储层越接近于成层状性强的储层。
这样,本发明实施例利用上述步骤S110~步骤S130计算出了可用于量化表征当前储层成层状性强弱的层状指数参数。
在另一实施例中,也能够基于上述所述的方法,生成针对当前储层的层状指数曲线。图3为本申请实施例的用于表征砂砾岩储层层状指数的方法中生成层状指数曲线的具体流程图。如图3所述,首先,步骤S301利用电成像测井技术,采集砂砾岩储层的静态或动态图像。步骤S302按照预设的行间隔和列间隔,将当前砂砾岩储层的静态或动态图像进行网格化划分处理。需要说明的是,由于步骤S301与上述步骤S110类似,步骤S302与上述步骤S201类似,故在此不作赘述。
而后,进入到步骤S303中。步骤S303按照预设窗口、以及预设深度步长,将经过网格化划分处理的当前砂砾岩储层的静态或动态图像进行窗口扫描处理,构建关于每个扫描窗口的电导率矩阵。扫描窗口的横向宽度与当前静态或动态图像的宽度相符,预设的深度步长的长度一般小于扫描窗口的纵向宽度。其中,在当前砂砾岩储层静态或动态图像中,被扫描窗口所框定的部分即为待处理图像。在步骤S303中,需要按照预先设定的扫描窗口,并按照预设的深度步长,沿着深度方向(向上或向下)移动扫描窗口进行扫描,以构建每个扫描窗口的电导率矩阵,从而进入到步骤S304中。步骤S304计算每个扫描窗口(将每个扫描窗口作为待处理的图像)的横向连续性指数,而后,步骤S305计算每个扫描窗口的纵向连续性指数,接着,步骤S306计算每个扫描窗口的层状指数。需要说明的是,由于步骤S304与上述步骤S1203类似,步骤S305与上述步骤S1204类似,步骤S306与上述步骤S130类似,故在此不作赘述。最后,步骤S307将所有层状指数按照每个扫描窗口对应的储层深度范围进行连接,生成针对当前储层的层状指数曲线。
在本发明实施例中,需要利用电成像测井技术采集待研究区域内不同储层的图像。在砂砾岩储层中,利用电成像图像资料量化表征储层的成层状性,通过计算图像横向连续性和纵向连续性,在此基础上构建砂砾岩储层的层状指数。下面以EL油田腾一段砂砾岩储层为例结合附图对本发明生成层状指数参数的实施例做进一步说明,包括以下步骤:
(1)基于电成像测井资料,选择出研究区块中多个具有代表性的不同成层状性的砂砾岩储层图像,作为该区块储层的研究样本。如图4所示为两例不同成层状性的砂砾岩储层图像。
(2)获取储层的动态或静态图像电导率矩阵。
(3)依次利用式(1)、式(2)和式(3)计算当前储层的横向连续性指数。
(4)依次利用式(4)、式(5)和式(6)计算当前储层的纵向连续性指数。
(6)利用纵向连续性指数CONT_V与横向连续性CONT_H指数的比值,构建储层层状指数LI。层状指数越小,表明储层的成层状性越弱,越接近于成层状性弱的储层;反之,越接近成层状性越强的储层。
进一步,下面以EL油田腾一段砂砾岩储层为例结合附图对本发明生成层状指数曲线的实施例做进一步说明,包括以下步骤:
(1)采集某深度段内储层的电成像图像;
(2)基于预设窗口,按照不同深度位置自上而下,对储层图像进行窗口扫描;
(3)构建每个扫描窗口的电导率矩阵;
(4)根据每个电导率矩阵依次计算横向连续性指数、纵向连续性指数和层状指数;
(5)按照扫描窗口对应的深度位置(范围),将层状指数进行连接,绘制该深度段储层的层状指数曲线。
图6为本申请实施例的用于表征砂砾岩储层层状指数的方法中层状指数计算实例示意图。如图6所示,Tt_X_Continuity_d、Tt_Y_Continuity_d和Tt_XY_LayerIndex_d曲线为基于动态图像的储层横向连续性、储层纵向连续性和储层的层状指数;Tt_X_Continuity、Tt_Y_Continuity和Tt_XY_LayerIndex曲线为基于静态图像的储层横向连续性、储层纵向连续性和储层的层状指数。在实际处理过程中,使用静态图像和动态图像均可,但一般使用动态图像。
此外,在完成层状性指数计算后,由于储层的层状指数能够反映出该储层的产能情况,因此,为了建立层状指数与产能情况之间的关系,进一步更加方便、快捷的识别出待研究区域内不同位置的储层是否为有效储层,在本发明所述的量化表征方法中,还能够利用不同储层对应的层状指数,确定出用于区分待研究区域内的储层是否有效的层状指数标准。
进一步,本发明所述的量化表征方法还包括:根据待研究区块内不同储层对应的层状指数,并结合不同储层对应的产能信息,确定用于识别待研究区块内储层是否有效的层状指数标准,以对待研究区块的储层是否有效来进行自主识别。
具体的,首先,收集待研究区块内多个具有代表性的储层的产能情况信息,作为标准确定样本集合,并对样本集合中每个样本进行储层属性归类。其中,储层属性包括:有效储层或无效储层。有效储层表示试油或投产情况达到工业产能,无效储层表示试油或投产情况属于干层。而后,利用上述步骤S110~步骤S130,分别计算上述样本集合中每个样本的层状指数。最后,将每个样本储层的储层属性与相应的层状指数进行关联,确定待研究储层的层状指数标准,从而利用这一层状指标准判断待研究区块内的储层有效性。其中,层状指数标准为能够自动识别待研究区块的储层是否有效的指标参数。
进一步,在判断当前待研究区块内的(待评价)某一储层是否有效时,将上述层状指数标准作为有效性评价阈值,并将待评价储层对应的层状指数与层状指数标准进行对比,利用二者的比较结果,判断当前待评价储层的有效性。其中,在当前待评价储层对应的层状指数大于当前层状指数标准时,判断当前储层为无效储层。另外,在当前待评价储层对应的层状指数小于当前层状指数标准时,判断当前储层为有效储层。
在一个具体示例中,基于电成像测井资料,选择出待研究区块中多个具有代表性的砂砾岩储层图像,作为该区块储层的研究样本集合。第一,依据当前待研究区试油和/或试采资料,确定各个储层样本的储层产能情况,将各样本储层分为有效储层和无效储层;第二,利用上述步骤S110~步骤S130建立不同砂砾岩储层样本的层状指数;第三,通过比较有效储层和无效储层的层状指数,并结合每个储层样本相对应的层状指数,快速建立有效储层的判识标准。图7为本申请实施例的用于表征砂砾岩储层层状指数的方法中确定层状指数标准的实例示意图。如图7所示,样本集合中选取19个样本层,其中,第1-8为无效储层,第9-19为有效储层,通过统计图可知,该区域有效储层和无效储层的层状指数界限约为1.5(如图7中虚线所示)。在待研究区域内储层的层状指数大于1.5时,判别当前储层为无效储层;在待研究区域内储层的层状指数小于1.5时,判别当前储层为有效储层,整体储层有效性的判识准确率为84.2%。
另一方面,基于上述用于表征砂砾岩储层层状指数的方法,本发明还提出了一种用于表征砂砾岩储层层状指数的系统。图8为本申请实施例的用于表征砂砾岩储层层状指数的系统的模块框图。如图8所示,该系统包括:图像获取模块81、连续性指数生成模块82和层状指数生成模块83。其中,图像获取模块81按照上述步骤S110所述的方法实施,配置为利用电成像测井技术,采集砂砾岩储层图像。连续性指数生成模块82按照上述步骤S120所述的方法实施,配置为按照预设的行间隔和列间隔,将当前图像进行网格化划分,由每行的行电导率值的标准方差,得到当前储层横向连续性指数,以及由每列的列电导率值的标准方差,得到相应的纵向连续性指数。层状指数生成模块83按照上述步骤S130所述的方法实施,配置为根据横向连续性指数和纵向连续性指数,计算当前储层的层状指数,用以定量表征当前储层成层状性的强弱。
进一步,上述连续性指数生成模块82包括:电导率矩阵构建单元821、横向连续性指数计算单元822和纵向连续性指数计算单元823。其中,电导率矩阵构建单元821,其配置为构建关于当前砂砾岩储层图像的电导率矩阵。横向连续性指数计算单元822,其配置为计算当前图像中每行的行电导率标准方差,并计算各行的行电导率标准方差的平均值,得到横向连续性指数。纵向连续性指数计算单元823,其配置为计算当前图像中每列的列电导率标准方差,并计算各列的列电导率标准方差的平均值,得到纵向连续性指数。
进一步,上述层状指数生成模块83,其配置为将上述纵向连续性指数与上述横向连续性指数进行除法运算,得到层状指数。
本发明所述的量化表征系统,还包括:曲线绘制模块84。曲线绘制模块84包括:窗口扫描单元841、窗口层状指数计算单元842和层状指数曲线生成单元843。其中,窗口扫描单元841,其配置为按照预设窗口,对经过网格化划分处理的当前图像进行扫描,构建关于每个扫描窗口的电导率矩阵。窗口层状指数计算单元842,其配置为通过执行上述连续性指数生成模块82和层状指数生成模块83,计算每个扫描窗口对应的横向连续性指数、纵向连续性指数和相应的层状指数。层状指数曲线生成单元843,其配置为将所有扫描窗口对应的层状指数按照储层深度进行连接,生成针对当前储层的层状指数曲线。
此外,本发明所述的量化表征系统,还包括:标准生成模块85。其中,标准生成模块85,其配置为根据待研究区块内不同储层对应的层状指数,并结合不同储层的产能信息,确定用于识别所述待研究区块内储层是否有效的层状指数标准,以利用这一层状指标准判断待研究区块内的储层有效性。
本发明公开了一种用于表征砂砾岩储层层状指数的方法及系统。该方法及系统基于电成像测井技术,利用储层电导率图像,计算横向连续性、以及纵向连续性,进一步得到表征储层成层状性强弱的量化指标(层状指数)。另外,利用构建储层层状指数这一量化指标的计算过程,结合待研究区块内不同储层样本的产能情况,确定用于评价当前待研究区块内储层有效性的关键层状指数标准,用以自动判断并识别出待研究区块内储层的有效性。本发明实现了储层成层状性的量化表征,并确定了有效储层的层状指数界限,为基于电成像的有效储层识别提供了有效手段,有效储层识别准确率达84.2%。
虽然本发明所揭露的实施方式如上,但所述的内容只是为了便于理解本发明而采用的实施方式,并非用以限定本发明。任何本发明所属技术领域内的技术人员,在不脱离本发明所揭露的精神和范围的前提下,可以在实施的形式上及细节上作任何的修改与变化,但本发明的专利保护范围,仍须以所附的权利要求书所界定的范围为准。

Claims (10)

1.一种用于表征砂砾岩储层层状指数的方法,其特征在于,包括:
步骤一、利用电成像测井技术,采集砂砾岩储层图像;
步骤二、按照预设的行间隔和列间隔,将所述图像进行网格化划分,由每行的行电导率值的标准方差,得到当前储层横向连续性指数,以及由每列的列电导率值的标准方差,得到相应的纵向连续性指数;
步骤三、根据所述横向连续性指数和所述纵向连续性指数,计算当前储层的层状指数,用以定量表征当前储层成层状性的强弱。
2.根据权利要求1所述的方法,其特征在于,所述方法,还包括:
根据待研究区块内不同储层对应的所述层状指数,并结合所述不同储层的产能信息,确定用于识别所述待研究区块内储层是否有效的层状指数标准,以利用所述层状指数标准判断所述待研究区块内的储层有效性。
3.根据权利要求1或2所述的方法,其特征在于,在所述步骤二中,
构建关于所述砂砾岩储层图像的电导率矩阵;
计算所述图像中每行的行电导率标准方差,并计算各行的行电导率标准方差的平均值,得到所述横向连续性指数;
计算所述图像中每列的列电导率标准方差,并计算各列的列电导率标准方差的平均值,得到所述纵向连续性指数。
4.根据权利要求1~3中任一项所述的方法,其特征在于,
将所述纵向连续性指数与所述横向连续性指数进行除法运算,得到所述层状指数。
5.根据权利要求1所述的方法,其特征在于,所述方法,还包括:
按照预设窗口,对经过网格化划分处理的所述图像进行扫描,构建关于每个扫描窗口的电导率矩阵;
通过执行所述步骤二和所述步骤三,计算每个扫描窗口对应的所述横向连续性指数、所述纵向连续性指数和相应的所述层状指数;
将所有所述层状指数按照储层深度进行连接,生成针对当前储层的层状指数曲线。
6.一种用于表征砂砾岩储层层状指数的系统,其特征在于,包括:
图像获取模块,其配置为利用电成像测井技术,采集砂砾岩储层图像;
连续性指数生成模块,其配置为按照预设的行间隔和列间隔,将所述图像进行网格化划分,由每行的行电导率值的标准方差,得到当前储层横向连续性指数,以及由每列的列电导率值的标准方差,得到相应的纵向连续性指数;
层状指数生成模块,其配置为根据所述横向连续性指数和所述纵向连续性指数,计算当前储层的层状指数,用以定量表征当前储层成层状性的强弱。
7.根据权利要求6所述的系统,其特征在于,所述系统,还包括:
标准生成模块,其配置为根据待研究区块内不同储层对应的所述层状指数,并结合所述不同储层的产能信息,确定用于识别所述待研究区块内储层是否有效的层状指数标准,以利用所述层状指数标准判断所述待研究区块内的储层有效性。
8.根据权利要求6或7所述的系统,其特征在于,所述连续性指数生成模块,包括:
电导率矩阵构建单元,其配置为构建关于所述砂砾岩储层图像的电导率矩阵;
横向连续性指数计算单元,其配置为计算所述图像中每行的行电导率标准方差,并计算各行的行电导率标准方差的平均值,得到所述横向连续性指数;
纵向连续性指数计算单元,其配置为计算所述图像中每列的列电导率标准方差,并计算各列的列电导率标准方差的平均值,得到所述纵向连续性指数。
9.根据权利要求6~8中任一项所述的系统,其特征在于,
所述层状指数生成模块,其进一步配置为将所述纵向连续性指数与所述横向连续性指数进行除法运算,得到所述层状指数。
10.根据权利要求6所述的系统,其特征在于,所述系统还包括曲线绘制模块,其中,所述曲线绘制模块具备:
窗口扫描单元,其配置为按照预设窗口,对经过网格化划分处理的所述图像进行扫描,构建关于每个扫描窗口的电导率矩阵;
窗口层状指数计算单元,其配置为执行所述连续性指数生成模块和所述层状指数生成模块,计算每个扫描窗口对应的所述横向连续性指数、所述纵向连续性指数和相应的所述层状指数;
层状指数曲线生成单元,其配置为将所有所述层状指数按照储层深度进行连接,生成针对当前储层的层状指数曲线。
CN202010000508.7A 2020-01-02 2020-01-02 一种砂砾岩储层成层状性的评价方法及系统 Expired - Fee Related CN110966000B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010000508.7A CN110966000B (zh) 2020-01-02 2020-01-02 一种砂砾岩储层成层状性的评价方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010000508.7A CN110966000B (zh) 2020-01-02 2020-01-02 一种砂砾岩储层成层状性的评价方法及系统

Publications (2)

Publication Number Publication Date
CN110966000A true CN110966000A (zh) 2020-04-07
CN110966000B CN110966000B (zh) 2021-06-29

Family

ID=70037957

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010000508.7A Expired - Fee Related CN110966000B (zh) 2020-01-02 2020-01-02 一种砂砾岩储层成层状性的评价方法及系统

Country Status (1)

Country Link
CN (1) CN110966000B (zh)

Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040225441A1 (en) * 2003-05-08 2004-11-11 Schlumberger Technology Corporation, A Texas Corporation Automated borehole geology and petrophysics interpretation using image logs
CN102200008A (zh) * 2010-03-26 2011-09-28 中国石油天然气股份有限公司 一种基于电成像测井的储层有效性识别方法
CN103115820A (zh) * 2013-03-07 2013-05-22 北京科技大学 一种节理岩体各向异性的确定方法
WO2014116896A1 (en) * 2013-01-25 2014-07-31 Services Petroliers Schlumberger Pressure transient testing with sensitivity analysis
WO2014130342A1 (en) * 2013-02-20 2014-08-28 Apache Corporation Methods for determining well log attributes for formation characterization
US20140297237A1 (en) * 2011-11-30 2014-10-02 Baker Hughes Incorporated Use of Monte Carlo Simulations For Establishing Correlations and Their Uncertainties
US20150098627A1 (en) * 2013-10-08 2015-04-09 Shin-Ju Chu Ye Automatic Dip Picking From Wellbore Azimuthal Image Logs
CN104732562A (zh) * 2015-03-24 2015-06-24 中国石油集团川庆钻探工程有限公司 一种电成像测井图像相似度比对方法
CN104793254A (zh) * 2015-04-28 2015-07-22 中国科学院地质与地球物理研究所 一种地下电性薄层的识别方法和装置
CN104808248A (zh) * 2015-04-22 2015-07-29 长江大学 基于电成像测井全井眼图像的地层产状检测方法
CN104849766A (zh) * 2015-05-19 2015-08-19 中国石油集团川庆钻探工程有限公司 电成像测井图像的地层产状识别方法
CN104912547A (zh) * 2014-03-11 2015-09-16 中国石油化工集团公司 应用电阻率成像测井资料连续定量评价储层非均质特征的方法
US20160146970A1 (en) * 2014-11-20 2016-05-26 Apache Corporation Automated method for selecting positions within formations from which to extract samples thereof
CN106447670A (zh) * 2016-08-31 2017-02-22 中国石油集团川庆钻探工程有限公司 基于电成像测井图像的孔洞参数自动计算方法
WO2017058692A1 (en) * 2015-09-30 2017-04-06 Schlumberger Technology Corporation Methods and systems to analyze bed boundary detection
US20170103554A1 (en) * 2013-03-28 2017-04-13 Reeves Wireline Technologies Limited Borehole Log Data Processing Methods
CN107219159A (zh) * 2017-04-19 2017-09-29 中国石油天然气股份有限公司 一种确定储层非均质性的方法和装置
CN109085663A (zh) * 2018-07-20 2018-12-25 中国矿业大学 一种致密砂岩储层层理缝识别方法

Patent Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040225441A1 (en) * 2003-05-08 2004-11-11 Schlumberger Technology Corporation, A Texas Corporation Automated borehole geology and petrophysics interpretation using image logs
CN102200008A (zh) * 2010-03-26 2011-09-28 中国石油天然气股份有限公司 一种基于电成像测井的储层有效性识别方法
US20140297237A1 (en) * 2011-11-30 2014-10-02 Baker Hughes Incorporated Use of Monte Carlo Simulations For Establishing Correlations and Their Uncertainties
WO2014116896A1 (en) * 2013-01-25 2014-07-31 Services Petroliers Schlumberger Pressure transient testing with sensitivity analysis
WO2014130342A1 (en) * 2013-02-20 2014-08-28 Apache Corporation Methods for determining well log attributes for formation characterization
CN103115820A (zh) * 2013-03-07 2013-05-22 北京科技大学 一种节理岩体各向异性的确定方法
US20170103554A1 (en) * 2013-03-28 2017-04-13 Reeves Wireline Technologies Limited Borehole Log Data Processing Methods
US20150098627A1 (en) * 2013-10-08 2015-04-09 Shin-Ju Chu Ye Automatic Dip Picking From Wellbore Azimuthal Image Logs
CN104912547A (zh) * 2014-03-11 2015-09-16 中国石油化工集团公司 应用电阻率成像测井资料连续定量评价储层非均质特征的方法
US20160146970A1 (en) * 2014-11-20 2016-05-26 Apache Corporation Automated method for selecting positions within formations from which to extract samples thereof
CN104732562A (zh) * 2015-03-24 2015-06-24 中国石油集团川庆钻探工程有限公司 一种电成像测井图像相似度比对方法
CN104808248A (zh) * 2015-04-22 2015-07-29 长江大学 基于电成像测井全井眼图像的地层产状检测方法
CN104793254A (zh) * 2015-04-28 2015-07-22 中国科学院地质与地球物理研究所 一种地下电性薄层的识别方法和装置
CN104849766A (zh) * 2015-05-19 2015-08-19 中国石油集团川庆钻探工程有限公司 电成像测井图像的地层产状识别方法
WO2017058692A1 (en) * 2015-09-30 2017-04-06 Schlumberger Technology Corporation Methods and systems to analyze bed boundary detection
CN106447670A (zh) * 2016-08-31 2017-02-22 中国石油集团川庆钻探工程有限公司 基于电成像测井图像的孔洞参数自动计算方法
CN107219159A (zh) * 2017-04-19 2017-09-29 中国石油天然气股份有限公司 一种确定储层非均质性的方法和装置
CN109085663A (zh) * 2018-07-20 2018-12-25 中国矿业大学 一种致密砂岩储层层理缝识别方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JOSSELIN KHERROUBI等: "Lamination Analysis from Electrical Borehole Images A Quantitative Workflow", 《SPWLA 57TH ANNUAL LOGGING SYMPOSIUM》 *
RICHARD BOOTLE等: "Laminated Sand Shale Formation Evaluation Using Azimuthal LWD Resistivity", 《THE 2009 SPE ANNUAL TECHNICAL CONFERENCE AND EXHIBITION》 *
杨玉卿等: "基于电成像测井图像识别求取地层产状的新方法", 《石油物探》 *

Also Published As

Publication number Publication date
CN110966000B (zh) 2021-06-29

Similar Documents

Publication Publication Date Title
CN107703560B (zh) 一种基于三重信息的泥页岩岩相精细识别方法
CN110276827B (zh) 一种基于泥页岩储层的有效性的评价方法
US10228478B2 (en) Method of searching for oil-gas reservoir based on trap-3D software
Guo et al. Multiple-point geostatistics-based three-dimensional automatic geological modeling and uncertainty analysis for borehole data
US10641923B2 (en) Method for characterizing and exploiting a subterranean formation comprising a network of fractures
EP3217192A1 (en) Method for looking for oil and gas reservoir on the basis of trap-3d software
CN110866912B (zh) 基于成像测井图像纹理的页岩纹层非均质性数据处理方法
CN106875471A (zh) 煤系含或隔水层三维可视化建模方法
CN115166853A (zh) 页岩气藏天然裂缝模型建立方法、装置、电子设备及介质
CN104749651A (zh) 相控多级重构测井定量识别泥页岩甜点方法
CN114153002A (zh) 储层天然裂缝三维地质建模方法、装置、电子设备及介质
CN118759600A (zh) 一种适用于构造蚀变岩型金矿的找矿预测方法
Zhou et al. Fracture prediction of tight sandstone reservoirs using outcrops and log curve-based extremum method: a case study of the Chang 7 Member of the Yanchang Formation in Block X, Ordos Basin
CN116559938A (zh) 一种油气藏裂缝模型建立方法及电子设备
CN110966000B (zh) 一种砂砾岩储层成层状性的评价方法及系统
Sarkheil et al. The fracture network modeling in naturally fractured reservoirs using artificial neural network based on image loges and core measurements
CN116263901A (zh) 页岩气开发评价方法及系统
CN119294030A (zh) 深层超压高温条件下数字岩心电阻率数值模拟及饱和度模型构建方法
CN119090666A (zh) 一种瓦斯储量预测及抽采效果动态评价系统
CN115266836B (zh) 一种地层流体类型识别的方法及装置
CN119378179A (zh) 一种用于反演走滑断裂形成期构造应力的方法及系统
CN119066827A (zh) 一种碳酸盐岩潜山油藏改建储气库评价方法及装置
CN117908156A (zh) 一种基于参数耦合评价模型评价页岩油地质的方法及系统
CN115831276A (zh) 一种土料的空间分类分区方法和系统
CN117687096B (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
TA01 Transfer of patent application right

Effective date of registration: 20200917

Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No.

Applicant after: CHINA NATIONAL PETROLEUM Corp.

Applicant after: CHINA PETROLEUM LOGGING Co.,Ltd.

Applicant after: North China branch of CNPC logging Co.,Ltd.

Address before: 062550 No. 002, Bohai East Road, Renqiu City, Cangzhou City, Hebei Province

Applicant before: North China branch of CNPC logging Co.,Ltd.

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210629

CF01 Termination of patent right due to non-payment of annual fee