CN106569271B - 储层裂缝显示方法及装置 - Google Patents

Abstract

本发明公开了一种储层裂缝显示方法及装置，其涉及石油勘探技术领域，所述储层裂缝显示方法包括以下步骤：获取目标储层的具有三维参数的地震资料；基于所述地震资料获得所述目标储层的裂缝发育特征，所述裂缝发育特征至少包括裂缝密度矩阵、裂缝发育方向矩阵和裂缝纵横比矩阵；基于所述裂缝密度矩阵、所述裂缝发育方向矩阵和所述裂缝纵横比矩阵得到裂缝延伸长度和张开度；所述裂缝延伸长度和所述张开度得到裂缝延伸权重矩阵；通过所述裂缝密度矩阵中各点的裂缝密度与所述裂缝延伸权重矩阵得到裂缝密度。应用本储层裂缝显示方法及装置形成的裂缝成图结果能够有效展现目标储层裂缝张开度。

Description

储层裂缝显示方法及装置

技术领域

本发明涉及石油勘探技术领域，特别涉及一种储层裂缝显示方法及装置。

背景技术

裂缝是岩石受力发生破裂、沿破裂面两侧的岩石没有发生明显位移的一种断裂构造现象。它既是油、气、水等地下流体的运移通道，也是流体的储集空间。通过裂缝可以极大地提升储层的渗透率，所以裂缝探测对于石油勘探和开发均有重要意义。

关于裂缝的裂缝密度和发育方向是描述目标储层裂缝的两个较为重要参数。裂缝密度表示裂缝在储层中的线、面积或体积百分比，而裂缝方向表示裂缝平面的延伸方向。为了同时表达裂缝密度和方向，人们通常用裂缝矢量图显示裂缝发育特征，其中矢量长度与裂缝密度成正比，矢量方向则表示裂缝发育方向。但是上述显示方式没有考虑裂缝的延伸特征，未能表现出裂缝的张开度。

发明内容

为了克服现有技术的上述缺陷，本发明实施例中提供了一种储层裂缝显示方法及装置，其能够解决裂缝成图结果不能表现裂缝张开度的问题。

本发明实施例的具体技术方案是：

一种储层裂缝显示方法，所述储层裂缝显示方法包括以下步骤：

获取目标储层的具有三维参数的地震资料；

基于所述地震资料获得所述目标储层的裂缝发育特征，所述裂缝发育特征至少包括裂缝密度矩阵、裂缝发育方向矩阵和裂缝纵横比矩阵；

基于所述裂缝密度矩阵、所述裂缝发育方向矩阵和所述裂缝纵横比矩阵得到裂缝延伸长度和张开度；

基于所述裂缝延伸长度和所述张开度得到裂缝延伸权重矩阵；

通过所述裂缝密度矩阵中各点的裂缝密度与所述裂缝延伸权重矩阵得到新的裂缝密度。

在一个优选地实施方式中，在所述基于所述地震资料获得所述目标储层的裂缝发育特征的步骤中，具体为对所述地震资料进行各向异性分析进而获得所述目标储层的裂缝发育特征。

在一个优选地实施方式中，所述裂缝密度矩阵中至少包括各点的裂缝密度。

在一个优选地实施方式中，在所述基于所述裂缝密度矩阵、所述裂缝发育方向矩阵和所述裂缝纵横比矩阵得到裂缝延伸长度和张开度的步骤中，其具体的计算过程如下：

ρ_max＝max(P_ni×nx)

ρ_o＝P(i,j)

r_s＝r_lα

φ_o＝Φ_ni×nx(i,j)

其中，r_l表示裂缝延伸长度，r_s表示张开度，ρ_max表示裂缝密度的最大值点，r_max表示根据地质调查和测井裂缝研究得到裂缝密度最大值点的裂缝延伸长度，ρ_max表示裂缝密度最大值点，α表示裂缝纵横比矩阵中的纵横比，ρ_O表示裂缝密度矩阵中某一点的裂缝密度，ρ_O中O点对应的椭圆盘长轴方向与裂缝发育方向相同，P_ni×nx表示多个点的裂缝密度，ni为目标储层工区内inline线的数目，nx为目标储层工区内xline线的数目，i，j为点O对应的线道号，i＝1,2,...,ni，j＝1,2,...,nx，Φ_ni×nx(i,j)表示目标储层工区内不同的点的椭圆盘长轴方向，φ_o表示裂缝发育方向矩阵中的裂缝发育方向。

在一个优选地实施方式中，所述裂缝延伸权重矩阵中包括各点的权值，权值的具体计算过程如下：

γ_M＝d_M/r_M

ω_M＝f(γ_M)

f(γ_M＝0)＝1

f(γ_M≥1)＝0

其中，ω_M表示权值，权值是该权值位置的点与椭圆盘中心O之间的距离d_M与椭圆盘该方向极径r_M比值γ_M的函数，椭圆盘内各点权值从中心等于1递减至边界处等于0，权重等值线为椭圆。

在一个优选地实施方式中，在所述通过所述裂缝密度矩阵中各点的裂缝密度与所述裂缝延伸权重矩阵得到新的裂缝密度的步骤中，将所述裂缝密度矩阵中各点的裂缝密度与所述裂缝延伸权重矩阵的乘积进行叠加以对储层内每点的裂缝进行延伸扩展进而生成新的裂缝密度，其具体计算公式如下：

或

或

其中，W_O表示裂缝延伸权重矩阵，ρ_O表示裂缝密度矩阵中各点的裂缝密度，P_ni×nx表示裂缝密度，ni为工区内inline线的数目，nx为工区内xline线的数目。

一种储层裂缝显示装置，所述储层裂缝显示装置包括：

地震资料获取模块，其用于获取目标储层的具有三维参数的地震资料；

裂缝发育特征生成模块，其用于对所述地震资料进行各向异性分析进而获得所述目标储层的裂缝发育特征；

裂缝延伸长度和张开度计算模块，其用于根据所述裂缝密度矩阵、所述裂缝发育方向矩阵和所述裂缝纵横比矩阵计算得到裂缝延伸长度和张开度；

裂缝延伸权重矩阵生成模块，其用于根据所述裂缝延伸长度和所述张开度得到裂缝延伸权重矩阵；

新的裂缝密度生成模块，其用于根据所述裂缝密度矩阵中各点的裂缝密度与所述裂缝延伸权重矩阵得到新的裂缝密度。

在一个优选地实施方式中，在所述裂缝延伸长度和张开度计算模块中，根据所述裂缝密度矩阵、所述裂缝发育方向矩阵和所述裂缝纵横比矩阵计算得到裂缝延伸长度和张开度的计算过程如下，

ρ_max＝max(P_ni×nx)

ρ_o＝P(i,j)

r_s＝r_lα

φ_o＝Φ_ni×nx(i,j)

其中，r_l表示裂缝延伸长度，r_s表示张开度，ρ_max表示裂缝密度的最大值点，r_max表示根据地质调查和测井裂缝研究得到裂缝密度最大值点的裂缝延伸长度，ρ_max表示裂缝延伸长度，α表示裂缝纵横比矩阵中的纵横比，ρ_O表示裂缝密度矩阵中某一点的裂缝密度，ρ_O中O点对应的椭圆盘长轴方向与裂缝发育方向相同，P_ni×nx表示多个点的裂缝密度，ni为工区内inline线的数目，nx为工区内xline线的数目，i，j为点O对应的线道号，i＝1,2,...,ni，j＝1,2,...,nx，Φ_ni×nx(i,j)表示目标储层工区内不同的点的椭圆盘长轴方向，φ_o表示裂缝发育方向矩阵中的裂缝发育方向。

在一个优选地实施方式中，在所述裂缝延伸权重矩阵生成模块中，所述裂缝延伸权重矩阵中包括各点的权值，权值的具体计算过程如下：

γ_M＝d_M/r_M

ω_M＝f(γ_M)

f(γ_M＝0)＝1

f(γ_M≥1)＝0

其中，ω_M表示权值，权值是该权值位置的点与椭圆盘中心O之间的距离d_M与椭圆盘该方向极径r_M比值γM的函数，椭圆盘内各点权值从中心等于1递减至边界处等于0，权重等值线为椭圆。

在一个优选地实施方式中，在所述新的裂缝密度生成模块中，将所述裂缝密度矩阵中各点的裂缝密度与所述裂缝延伸权重矩阵的乘积进行叠加以对储层内每点的裂缝进行延伸扩展进而生成新的裂缝密度，其具体计算公式如下：

或

或

其中，W_O表示裂缝延伸权重矩阵，ρ_O表示裂缝密度矩阵中各点的裂缝密度，P_ni×nx表示裂缝密度，ni为目标储层工区内inline线的数目，nx为目标储层工区内xline线的数目。

本发明实施例中的储层裂缝显示方法及装置，其通过裂缝密度矩阵、裂缝发育方向矩阵和裂缝纵横比矩阵得到裂缝延伸长度和张开度，进而根据裂缝延伸长度和张开度形成裂缝延伸权重矩阵，裂缝密度矩阵中各点的裂缝密度与裂缝延伸权重矩阵的乘积进行叠加以对储层内每点的裂缝进行延伸扩展最终生成新的裂缝密度，考虑到了裂缝的延伸特征，新的裂缝密度形成的沿层切片图像能有效展现了目标储层的裂缝张开度、裂缝密度和裂缝发育方向矩阵，该新的裂缝密度形成的切片图像可以用于分析储层裂缝发育形态和连通关系。

附图说明

在此描述的附图仅用于解释目的，而不意图以任何方式来限制本发明公开的范围。另外，图中的各部件的形状和比例尺寸等仅为示意性的，用于帮助对本发明的理解，并不是具体限定本发明各部件的形状和比例尺寸。本领域的技术人员在本发明的教导下，可以根据具体情况选择各种可能的形状和比例尺寸来实施本发明。

图1是本发明实施例提供的一种目标储层裂缝显示方法的流程图。

图2是本发明实施例提供的一种计算裂缝延伸权重矩阵参数的示意图。

图3是本发明实施例提供的一种计算裂缝延伸权重矩阵内某一点权值的示意图。

图4是本发明实施例提供的一种裂缝延伸权重矩阵的示意图。

图5是本发明实施例提供的裂缝延伸权重矩阵加权叠加后的效果图。

图6是本发明实施例提供的一种储层裂缝显示装置的结构框图。

具体实施方式

结合附图和本发明具体实施方式的描述，能够更加清楚地了解本发明的细节。但是，在此描述的本发明的具体实施方式，仅用于解释本发明的目的，而不能以任何方式理解成是对本发明的限制。在本发明的教导下，技术人员可以构想基于本发明的任意可能的变形，这些都应被视为属于本发明的范围。

图1是本发明实施例提供的一种目标储层裂缝显示方法的流程图，如图1所示，在本申请中提出了一种储层裂缝显示方法，储层裂缝显示方法包括以下步骤：

S101：获取目标储层的具有三维参数的地震资料。

S102：基于地震资料获得目标储层的裂缝发育特征，裂缝发育特征至少包括裂缝密度矩阵、裂缝发育方向矩阵和裂缝纵横比矩阵。

由于地震资料具有三维参数，所以可以对获取的目标储层的地震资料进行各向异性分析，如此可以获得目标储层的裂缝发育特征。在裂缝发育特征中至少需要包括裂缝密度矩阵、裂缝发育方向矩阵和裂缝纵横比矩阵等参数。为了下面计算的需要，将目标储层的裂缝密度以裂缝密度矩阵的方式表示，裂缝密度矩阵可以表示为P_ni×nx，其具体为ni乘以nx的矩阵，其中，ni为目标储层工区内inline线的数目，其中nx为目标储层工区内xline线的数目，所以裂缝密度矩阵中包括各点的裂缝密度，即裂缝密度矩阵中的元素为各点的裂缝密度。

S103：基于裂缝密度矩阵、裂缝发育方向矩阵和裂缝纵横比矩阵得到裂缝延伸长度和张开度。

根据裂缝密度矩阵、裂缝发育方向矩阵和裂缝纵横比矩阵计算得到裂缝延伸长度和张开度的，其具体的计算过程如下：

ρ_max＝max(P_ni×nx)

ρ_o＝P(i,j)

r_s＝r_lα

φ_o＝Φ_ni×nx(i,j)

其中，r_l表示裂缝延伸长度，r_s表示张开度，ρ_max表示裂缝密度的最大值点，r_max表示根据地质调查和测井裂缝研究得到的最大裂缝延伸长度，ρ_max表示裂缝密度最大值点，α表示裂缝纵横比矩阵中的纵横比，ρ_O表示裂缝密度矩阵中某一点的裂缝密度，ρ_O中O点对应的椭圆盘长轴方向与裂缝发育方向相同，P_ni×nx表示裂缝密度，ni为目标储层工区内inline线的数目，nx为目标储层工区内xline线的数目，i，j为点O对应的线道号，i＝1,2,...,ni，j＝1,2,...,nx，Φ_ni×nx(i,j)表示目标储层工区内不同的点的椭圆盘长轴方向，φ_o表示裂缝发育方向矩阵中裂缝发育方向。

在上述过程中，先求出目标储层工区所有点的裂缝密度，进而得到其中裂缝密度的最大值点ρ_max，根据裂缝密度的最大值点ρ_max和根据地质调查、测井裂缝研究得到的最大裂缝延伸长度r_max和裂缝密度矩阵中某一点的裂缝密度ρ_O可以计算得到该点的裂缝延伸长度r_l，再通过裂缝纵横比矩阵中的纵横比α计算得到该点的r_s。通过上述方法，计算得到目标储层工区内各点裂缝的延伸长度和张开度。

S104：基于裂缝延伸长度、张开度和裂缝发育方向矩阵得到裂缝延伸权重矩阵。

图2是本发明实施例提供的一种计算裂缝延伸权重矩阵参数的示意图，图3是本发明实施例提供的一种计算裂缝延伸权重矩阵内某一点权值的示意图，如图2、图3，根据目标储层工区内各点裂缝的延伸长度和张开度得到裂缝延伸权重矩阵W_O。该裂缝延伸权重矩阵内的非零值呈椭圆盘形状，该椭圆盘的尺寸参数由裂缝延伸长度和张开度确定，裂缝延伸权重矩阵是ni乘以nx大小的矩阵，如图3所示，当O点固定时，网格nx乘以ni上的每一个网格点(例如，iM与jM位置的黑点M)，就可以得到该点的权值ω_M。如图2所示，权值ω_M是点M与圆盘中心O之间的距离d_M与椭圆盘该方向极径r_M比值γ_M的函数，其计算过程如下：

γ_M＝d_M/r_M

ω_M＝f(γ_M)

f(γ_M＝0)＝1

f(γ_M≥1)＝0

椭圆盘内各点权值从中心等于1递减至边界处等于0，权重等值线均为椭圆，如公式所示，裂缝延伸权重矩阵W_O在椭圆内是有值的，中心处等于1，边界处等于0，如此，计算得到裂缝延伸权重矩阵W_O内每一个元素的值，图4是本发明实施例提供的一种裂缝延伸权重矩阵的示意图，如图4所示，图中椭圆盘的不同颜色表示裂缝延伸权重矩阵内元素的大小。

S105：通过裂缝密度矩阵中各点的裂缝密度与裂缝延伸权重矩阵得到新的裂缝密度。

将裂缝密度矩阵中各点的裂缝密度与裂缝延伸权重矩阵的乘积进行叠加以对储层内每点的裂缝进行延伸扩展进而生成新的裂缝密度，其具体计算公式如下：

或

或

其中，W_O表示裂缝延伸权重矩阵，ρ_O表示裂缝密度矩阵中各点的裂缝密度，P_ni×nx表示多个点下裂缝密度，ni为工区内inline线的数目，nx为工区内xline线的数目，P_new表示新的裂缝密度。

图5是本发明实施例提供的裂缝延伸权重矩阵加权叠加后的效果图，如图5所示，当形成新的裂缝密度后，可以通过沿层切片和伪彩图方式成图以显示该新的裂缝密度，由于考虑到了裂缝的延伸特征，新的裂缝密度形成的沿层切片图像能有效展现了目标储层的裂缝张开度、裂缝密度和裂缝发育方向矩阵，该新的裂缝密度形成的切片图像可以用于分析储层裂缝发育形态和连通关系。

在本发明的实施例中还提出了一种储层裂缝显示装置，图6是本发明实施例提供的一种储层裂缝显示装置的结构框图，如图6所示，储层裂缝显示装置包括：地震资料获取模块，裂缝发育特征生成模块，裂缝延伸长度和张开度计算模块，裂缝延伸权重矩阵生成模块，新的裂缝密度生成模块。

地震资料获取模块用于获取目标储层的具有三维参数的地震资料。

裂缝发育特征生成模块用于对地震资料进行各向异性分析进而获得目标储层的裂缝发育特征。

裂缝延伸长度和张开度计算模块用于根据裂缝密度矩阵、裂缝发育方向矩阵和裂缝纵横比矩阵计算得到裂缝延伸长度和张开度，计算过程如下，

ρ_max＝max(P_ni×nx)

ρ_o＝P(i,j)

r_s＝r_lα

φ_o＝Φ_ni×nx(i,j)

其中，r_l表示裂缝延伸长度，r_s表示张开度，ρ_max表示裂缝密度的最大值点，r_max表示根据地质调查和测井裂缝研究得到裂缝密度最大值点的裂缝延伸长度，ρ_max表示裂缝延伸长度，α表示裂缝纵横比矩阵中的纵横比，ρ_O表示裂缝密度矩阵中某一点的裂缝密度，ρ_O中O点对应的椭圆盘长轴方向与裂缝发育方向相同，P_ni×nx表示多个点的裂缝密度，ni为工区内inline线的数目，nx为工区内xline线的数目，i，j为点O对应的线道号，i＝1,2,...,ni，j＝1,2,...,nx，Φ_ni×nx(i,j)表示目标储层工区内不同的点的椭圆盘长轴方向，φ_o表示裂缝发育方向矩阵中的裂缝发育方向。

裂缝延伸权重矩阵生成模块用于根据裂缝延伸长度和张开度得到裂缝延伸权重矩阵，裂缝延伸权重矩阵中包括各点的权值，权值的具体计算过程如下：

γ_M＝d_M/r_M

ω_M＝f(γ_M)

f(γ_M＝0)＝1

f(γ_M≥1)＝0

其中，ω_M表示权值，权值是该权值位置的点与椭圆盘中心O之间的距离d_M与椭圆盘该方向极径r_M比值γ_M的函数，椭圆盘内各点权值从中心等于1递减至边界处等于0，权重等值线为椭圆。

新的裂缝密度生成模块用于根据裂缝密度矩阵中各点的裂缝密度与裂缝延伸权重矩阵得到新的裂缝密度，具体为将裂缝密度矩阵中各点的裂缝密度与裂缝延伸权重矩阵的乘积进行叠加以对储层内每点的裂缝进行延伸扩展进而生成新的裂缝密度，其具体计算公式如下：

或

或

其中，W_O表示裂缝延伸权重矩阵，ρ_O表示裂缝密度矩阵中各点的裂缝密度，P_ni×nx表示裂缝密度，ni为目标储层工区内inline线的数目，nx为目标储层工区内xline线的数目。

在20世纪90年代，对于一个技术的改进可以很明显地区分是硬件上的改进(例如，对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而，随着技术的发展，当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此，不能说一个方法流程的改进就不能用硬件实体模块来实现。例如，可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable GateArray，FPGA))就是这样一种集成电路，其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上，而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且，如今，取代手工地制作集成电路芯片，这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现，它与程序开发撰写时所用的软件编译器相类似，而要编译之前的原始代码也得用特定的编程语言来撰写，此称之为硬件描述语言(Hardware Description Language，HDL)，而HDL也并非仅有一种，而是有许多种，如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware DescriptionLanguage)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(RubyHardware Description Language)等，目前最普遍使用的是VHDL(Very-High-SpeedIntegrated Circuit Hardware Description Language)与Verilog。本领域技术人员也应该清楚，只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中，就可以很容易得到实现该逻辑方法流程的硬件电路。

控制器可以按任何适当的方式实现，例如，控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application Specific Integrated Circuit，ASIC)、可编程逻辑控制器和嵌入微控制器的形式，控制器的例子包括但不限于以下微控制器：ARC 625D、Atmel AT91SAM、Microchip PIC18F26K20以及Silicone Labs C8051F320，存储器控制器还可以被实现为存储器的控制逻辑的一部分。

本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。

上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。

为了描述的方便，描述以上装置时以功能分为各种单元分别描述。当然，在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。

通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。该计算机软件产品可以包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。该计算机软件产品可以存储在内存中，内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括短暂电脑可读媒体(transitory media)，如调制的数据信号和载波。

本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。

本申请可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。

本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。

虽然通过实施例描绘了本申请，本领域普通技术人员知道，本申请有许多变形和变化而不脱离本申请的精神，希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。

Claims (6)

1.一种储层裂缝显示方法，其特征在于，所述储层裂缝显示方法包括以下步骤：

获取目标储层的具有三维参数的地震资料；

基于所述地震资料获得所述目标储层的裂缝发育特征，所述裂缝发育特征至少包括裂缝密度矩阵、裂缝发育方向矩阵和裂缝纵横比矩阵；

基于所述裂缝密度矩阵、所述裂缝发育方向矩阵和所述裂缝纵横比矩阵得到裂缝延伸长度和张开度，其具体的计算过程如下：

ρ_max＝max(P_ni×nx)

ρ_o＝P(i,j)

r_s＝r_lα

φ_o＝Φ_ni×nx(i,j)

其中，r_l表示裂缝延伸长度，r_s表示张开度，ρ_max表示裂缝密度的最大值点，r_max表示根据地质调查和测井裂缝研究得到裂缝密度最大值点的裂缝延伸长度，α表示裂缝纵横比矩阵中的纵横比，ρ_O表示裂缝密度矩阵中某一点的裂缝密度，ρ_O中O点对应的椭圆盘长轴方向与裂缝发育方向相同，P_ni×nx表示多个点的裂缝密度，ni为目标储层工区内inline线的数目，nx为目标储层工区内xline线的数目，i，j为点O对应的线道号，i＝1,2,...,ni，j＝1,2,...,nx，Φ_ni×nx(i,j)表示目标储层工区内不同的点的椭圆盘长轴方向，φ_o表示裂缝发育方向矩阵中的裂缝发育方向；

基于所述裂缝延伸长度和所述张开度得到裂缝延伸权重矩阵，所述裂缝延伸权重矩阵内的非零值呈椭圆盘形状，该椭圆盘的尺寸参数由裂缝延伸长度和张开度确定，所述裂缝延伸权重矩阵中包括各点的权值，权值的具体计算过程如下：

γ_M＝d_M/r_M

ω_M＝f(γ_M)

f(γ_M＝0)＝1

f(γ_M≥1)＝0

其中，ω_M表示权值，权值是该权值位置的点与椭圆盘中心O之间的距离d_M与椭圆盘该方向极径r_M比值γ_M的函数，椭圆盘内各点权值从中心等于1递减至边界处等于0，权重等值线为椭圆；

通过所述裂缝密度矩阵中各点的裂缝密度与所述裂缝延伸权重矩阵得到新的裂缝密度。

2.根据权利要求1所述的储层裂缝显示方法，其特征在于，在所述基于所述地震资料获得所述目标储层的裂缝发育特征的步骤中，具体为对所述地震资料进行各向异性分析进而获得所述目标储层的裂缝发育特征。

3.根据权利要求1所述的储层裂缝显示方法，其特征在于，所述裂缝密度矩阵中至少包括各点的裂缝密度。

4.根据权利要求1所述的储层裂缝显示方法，其特征在于，在所述通过所述裂缝密度矩阵中各点的裂缝密度与所述裂缝延伸权重矩阵得到新的裂缝密度的步骤中，将所述裂缝密度矩阵中各点的裂缝密度与所述裂缝延伸权重矩阵的乘积进行叠加以对储层内每点的裂缝进行延伸扩展进而生成新的裂缝密度，其具体计算公式如下：

或

或

其中，W_O表示裂缝延伸权重矩阵，ρ_O表示裂缝密度矩阵中各点的裂缝密度，P_ni×nx表示裂缝密度，ni为工区内inline线的数目，nx为工区内xline线的数目。

5.一种储层裂缝显示装置，其特征在于，所述储层裂缝显示装置包括：

地震资料获取模块，其用于获取目标储层的具有三维参数的地震资料；

裂缝发育特征生成模块，其用于对所述地震资料进行各向异性分析进而获得所述目标储层的裂缝发育特征；

裂缝延伸长度和张开度计算模块，其用于根据所述裂缝密度矩阵、所述裂缝发育方向矩阵和所述裂缝纵横比矩阵计算得到裂缝延伸长度和张开度，计算过程如下，

ρ_max＝max(P_ni×nx)

ρ_o＝P(i,j)

r_s＝r_lα

φ_o＝Φ_ni×nx(i,j)

其中，r_l表示裂缝延伸长度，r_s表示张开度，ρ_max表示裂缝密度的最大值点，r_max表示根据地质调查和测井裂缝研究得到裂缝密度最大值点的裂缝延伸长度，ρ_max表示裂缝延伸长度，α表示裂缝纵横比矩阵中的纵横比，ρ_O表示裂缝密度矩阵中某一点的裂缝密度，ρ_O中O点对应的椭圆盘长轴方向与裂缝发育方向相同，P_ni×nx表示多个点的裂缝密度，ni为工区内inline线的数目，nx为工区内xline线的数目，i，j为点O对应的线道号，i＝1,2,...,ni，j＝1,2,...,nx，Φ_ni×nx(i,j)表示目标储层工区内不同的点的椭圆盘长轴方向，φ_o表示裂缝发育方向矩阵中的裂缝发育方向；

裂缝延伸权重矩阵生成模块，其用于根据所述裂缝延伸长度和所述张开度得到裂缝延伸权重矩阵，所述裂缝延伸权重矩阵内的非零值呈椭圆盘形状，该椭圆盘的尺寸参数由裂缝延伸长度和张开度确定，所述裂缝延伸权重矩阵中包括各点的权值，权值的具体计算过程如下：

γ_M＝d_M/r_M

ω_M＝f(γ_M)

f(γ_M＝0)＝1

f(γ_M≥1)＝0

其中，ω_M表示权值，权值是该权值位置的点与椭圆盘中心O之间的距离d_M与椭圆盘该方向极径r_M比值γ_M的函数，椭圆盘内各点权值从中心等于1递减至边界处等于0，权重等值线为椭圆；

新的裂缝密度生成模块，其用于根据所述裂缝密度矩阵中各点的裂缝密度与所述裂缝延伸权重矩阵得到新的裂缝密度。

6.根据权利要求5所述的储层裂缝显示装置，其特征在于，在所述新的裂缝密度生成模块中，将所述裂缝密度矩阵中各点的裂缝密度与所述裂缝延伸权重矩阵的乘积进行叠加以对储层内每点的裂缝进行延伸扩展进而生成新的裂缝密度，其具体计算公式如下：

或

或

其中，W_O表示裂缝延伸权重矩阵，ρ_O表示裂缝密度矩阵中各点的裂缝密度，P_ni×nx表示裂缝密度，ni为目标储层工区内inline线的数目，nx为目标储层工区内xline线的数目。