CN106772680A - 基于二维叠后地震资料微观裂缝表征确定方法和装置 - Google Patents

Abstract

本申请实施例提供了一种基于二维叠后地震资料微观裂缝表征确定方法和装置。其中，该方法通过铸体薄片和扫描电镜先确定了微观裂缝发育区，然后结合测井曲线及核磁测井孔隙度曲线，研究了微观裂缝的测井响应特征，再获得岩石物理量版；其次，根据岩石物理量板得到了微观裂缝测井特征曲线，再对二维叠后地震资料进行相干属性分析，进而可以通过波形指示反演确定微观裂缝分布特征。因此解决了现有的微观裂缝确定方法中存在的使用范围有限，确定精度低的技术问题。

Description

基于二维叠后地震资料微观裂缝表征确定方法和装置

技术领域

本申请涉及油气勘探领域，特别涉及一种基于二维叠后地震资料微观裂缝表征确定方法和装置。

背景技术

随着油气勘探开发程度的不断提高，油气勘探已经逐渐开始转入对低渗透-特低渗透致密砂岩储层的研究。例如，有些已经开始对页岩油气进行勘探开发。其中，上述低渗透-特低渗透致密砂岩储层通常是指满足孔隙度小于10％，且渗透率小于1md的砂岩储层。

目前，对于大尺度的宏观断裂已经有了一些研究，但是基于二维叠后地震资料对微观尺度的微观裂缝，即微米尺度的微观裂缝的研究较少。大量的勘探实践表明，微米尺度的微观裂缝可以改善低渗透-特低渗透致密砂岩储层，甚至对于更致密的泥页岩储层的物性条件也有一定的改善作用，对提高油气产能具有重要的意义。因此，准确地预测微观裂缝分布特征在油气勘探领域中非常重要。

但是，由于受限于地震资料分辨率的影响，目前基于二维叠后地震资料对发育在低渗透-特低渗透致密砂岩储层中的微米尺度微观裂缝的地球物理表征技术方法还没有被提出。通常情况下，在地震资料上较难发现微米尺度的微观裂缝的地震响应特征，这就给表征这种微米尺度的微观裂缝造成了较大的困难。现有的微观裂缝确定方法往往只能对取芯井段的样品进行扫描分析，只能确定取井芯段附近的微观裂缝的分布。即不能应用到连井剖面，甚至平面上，只能在点上实现，且具体使用范围受到取芯井段的限制。因此，具体实施时，现有的微观裂缝确定方法常常存在使用范围有限和确定微观裂缝分布的精度不高的技术问题。

针对上述问题，目前尚未提出有效的解决方案。

发明内容

本申请实施例提供了一种基于二维叠后地震资料微观裂缝表征确定方法和装置，以解决现有的微观裂缝确定方法存在的使用范围有限和确定微观裂缝的精度不高技术问题。

本申请实施例提供了一种基于二维叠后地震资料微观裂缝表征确定方法，包括：

获取目标区域的取心段砂岩的样品、所述目标区域的测井曲线资料、所述目标区域的核磁孔隙度测井资料和所述目标区域的二维叠后地震资料；

对所述样品进行铸体薄片及扫描电镜实验分析，确定所述样品的微观裂缝发育区；

根据所述目标区域的测井曲线资料、所述目标区域的核磁孔隙度测井资料和所述样品的微观裂缝发育区，确定所述目标区域的微观裂缝的测井响应特征；

将所述目标区域按照深度分成多个连续的小层，并根据所述微观裂缝的测井响应特征，为所述多个连续的小层中的各个小层分别建立岩石物理量版；

根据多个连续的小层中各个小层的岩石物理量版，确定微观裂缝测井特征曲线；

对所述目标区域的二维叠后地震资料进行相干属性分析，得到二维叠后相干地震资料；

以所述二维叠后相干地震资料为地震波形指示反演的数据体，根据所述微观裂缝测井特征曲线，通过地震波形指示反演，得到所述目标区域的微观裂缝的分布特征。

在一个实施方式中，所述微观裂缝的测井响应特征包括：纵波速度小于预设速度，和/或，核磁测井孔隙度低于预设孔隙度。

在一个实施方式中，根据所述目标区域的测井曲线资料、所述目标区域的核磁孔隙度测井资料和所述样品的微观裂缝发育区，确定所述目标区域的微观裂缝的测井响应特征，包括：

根据所述目标区域的测井曲线资料中的声波时差曲线，确定所述目标区域的纵波速度；

根据所述目标区域的核磁孔隙度测井资料，确定所述目标区域的核磁测井孔隙度；

根据所述样品的微观裂缝发育区，确定所述预设速度和所述预设孔隙度。

在一个实施方式中，根据所述微观裂缝的测井响应特征，为所述多个连续的小层中的各个小层分别建立岩石物理量版，包括：

根据所述微观裂缝的测井响应特征，对所述多个连续的小层中的各个小层，分别以所在小层的核磁测井孔隙度为横坐标，以纵波速度为纵坐标，建立各个小层的岩石物理量版。

在一个实施方式中，根据多个连续的小层中各个小层的岩石物理量版，确定微观裂缝测井特征曲线，包括：

从所述多个连续的小层中各个小层的岩石物理量版中逐小层获取微观裂缝的信息；

将获取的各个小层的微观裂缝的信息，作为所述微观裂缝测井特征曲线。

在一个实施方式中，以所述二维叠后相干地震资料为地震波形指示反演的数据体，根据所述微观裂缝测井特征曲线，通过地震波形指示反演，得到所述目标区域的微观裂缝的分布特征，包括：

以所述二维叠后相干地震资料为地震波形指示反演的数据体，并以所述微观裂缝测井特征曲线作为约束条件，优选高频成分，以建立反演模型；

根据所述反演模型，进行地震波形指示反演，得到所述目标区域的微观裂缝的分布特征。

在一个实施方式中，以所述二维叠后相干地震资料为地震波形指示反演的数据体，并以所述微观裂缝测井特征曲线作为约束条件，优选高频成分，以建立反演模型，包括：

以所述二维叠后相干地震资料为地震波形指示反演的数据体，并以所述微观裂缝测井特征曲线作为约束条件，优选高频成分，以建立初始模型；

根据所述初始模型，反演得到所述目标区域过井的二维叠后地震剖面上微观裂缝空间分布特征；

将所述过井二维叠后地震剖面上微观裂缝空间分布特征与所述目标区域的井上的微观裂缝测井特征曲线进行比较，得到比较误差；

如果所述比较误差小于预设阈值，则将所述初始模型作为所述反演模型。

在一个实施方式中，在将所述过井二维叠后地震剖面上微观裂缝空间分布特征与所述目标区域的井上的微观裂缝测井特征曲线进行比较，得到比较误差之后，所述方法还包括：

如果所述比较误差大于等于所述预设阈值，则重新选择高频成分，并根据重新选择的高频成分建立初始模型，直到根据建立的初始模型，反演得到的所述目标区域的过井二维叠后地震剖面上的微观裂缝空间分布特征与所述目标区域的井上的微观裂缝测井特征曲线之间的比较误差小于所述预设阈值，则将当前的初始模型作为所述反演模型。

基于相同的发明构思，本申请实施例还提供了一种基于二维叠后地震资料微观裂缝表征确定装置，包括：

获取模块，用于获取目标区域的取心段砂岩的样品、所述目标区域的测井曲线资料、所述目标区域的核磁孔隙度测井资料和所述目标区域的二维叠后地震资料；

第一确定模块，用于对所述样品进行铸体薄片及扫描电镜实验分析，确定所述样品的微观裂缝发育区；

第二确定模块，用于根据所述目标区域的测井曲线资料、所述目标区域的核磁孔隙度测井资料和所述样品的微观裂缝发育区，确定所述目标区域的微观裂缝的测井响应特征；

第三确定模块，用于将所述目标区域按照深度分成多个连续的小层，并根据所述微观裂缝的测井响应特征，为所述多个连续的小层中的各个小层分别建立岩石物理量版；

第四确定模块，用于根据多个连续的小层中各个小层的岩石物理量版，确定微观裂缝测井特征曲线；

相干属性分析模块，用于对所述目标区域的二维叠后地震资料进行相干属性分析，得到二维叠后相干地震资料；

第五确定模块，用于以所述二维叠后相干地震资料为地震波形指示反演的数据体，根据所述微观裂缝测井特征曲线，通过地震波形指示反演，得到所述目标区域的微观裂缝的分布特征。

在一个实施方式中，所述第二确定模块包括：

第一确定单元，用于根据所述目标区域的测井曲线资料中的声波时差曲线，确定所述目标区域的纵波速度；

第二确定单元，用于根据所述目标区域的核磁孔隙度测井资料，确定所述目标区域的核磁测井孔隙度；

第三确定单元，用于根据所述样品的微观裂缝发育区，确定所述预设速度和所述预设孔隙度。

在一个实施方式中，所述第四确定模块包括：

信息获取单元，用于从所述多个连续的小层中各个小层的岩石物理量版中逐小层获取微观裂缝的信息；

第四确定单元，用于将获取的各个小层的微观裂缝的信息，作为所述微观裂缝测井特征曲线。

在一个实施方式中，所述第五确定模块包括：

模型确定单元，用于以所述二维叠后相干地震资料为地震波形指示反演的数据体，并以所述微观裂缝测井特征曲线作为约束条件，优选高频成分，以建立反演模型；

第五确定单元，用于根据所述反演模型，进行地震波形指示反演，得到所述目标区域的微观裂缝的分布特征。

在本申请实施例中，通过利用铸体薄片和扫描电镜确定微观裂缝发育区，结合测井曲线和核磁测井孔隙度曲线获得微观裂缝的测井响应特征，参考建立的岩石物理量板进确定微观裂缝连续的测井特征曲线；再根据二维叠后相干地震资料和微观裂缝测井特征曲线通过地震波形指示反演确定微观裂缝的分布特征，扩大了预测的范围，解决了现有微观裂缝确定方法存在的使用范围有限和微观裂缝确定精度低的技术问题，达到了在二维叠后地震剖面上准确确定微观裂缝分布的技术效果。

附图说明

为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1是根据本申请实施例的基于二维叠后地震资料微观裂缝表征确定方法的处理流程图；

图2是根据本申请实施例的基于二维叠后地震资料微观裂缝表征确定装置的组成结构图；

图3是应用本申请实施例提供基于二维叠后地震资料微观裂缝表征确定方法/装置得到的目标区域中L1井测井资料及铸体和电镜照片综合柱状图；

图4是应用本申请实施例提供基于二维叠后地震资料微观裂缝表征确定方法/装置得到的致密砂岩储层微观裂缝岩石物理判别量版的示意图；

图5是应用本申请实施例提供基于二维叠后地震资料微观裂缝表征确定方法/装置得到的过L1井Crossline(道)方向原始二维叠后地震资料的示意图；

图6是应用本申请实施例提供基于二维叠后地震资料微观裂缝表征确定方法/装置得到的过L1井Crossline方向微观裂缝地震波形指示反演剖面的示意图。

具体实施方式

为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。

考虑到现有技术大多是利用地震资料对大尺度的裂缝进行预测，而不能直接根据地震资料确定微米尺度的微观裂缝分布特征。现有的微观裂缝的确定方法通常是对取心井段的样品进行实验分析，根据样品的分析结果确定微观裂缝发育区。但上述方法是基于取心井段的样品分析确定微观裂缝的发育区，因此受取心井段的限制，不能应用到连井剖面，甚至平面上，只能在点上实现。综上，可知现有的微观裂缝确定方法具体实施时，往往存在使用范围有限和微观裂缝分布特征确定精度低的技术问题。针对产生上述技术问题的根本原因，本申请考虑可以先通过铸体薄片和扫描电镜确定微观裂缝发育区，其次结合测井曲线及核磁测井孔隙度曲线，分析微观裂缝的测井响应特征，参考岩石物理量板，进而得到微观裂缝测井特征曲线；再将二维叠后地震资料进行相干属性分析，以使得可以通过使用二维叠后相干地震资料反应裂缝的特征；最后根据微观裂缝测井特征曲线和二维叠后相干地震资料，通过地震波形指示反演，确定微观裂缝的分布特征。从而解决了现有的微观裂缝确定方法存在的使用范围有限和微观裂缝确定精度低的技术问题，达到了在二维剖面上准确确定微观裂缝分布的技术效果。

基于上述思考思路，本申请提供了一种基于二维叠后地震资料微观裂缝表征确定方法。请参阅图1。本申请提供的基于二维叠后地震资料微观裂缝表征确定方法，具体可以包括以下步骤(步骤101至步骤107)。

步骤101：获取目标区域的取心段砂岩的样品、所述目标区域的测井曲线资料、所述目标区域的核磁孔隙度测井资料和所述目标区域的二维叠后地震资料。

在一个实施方式中，一般的探井都具有测井曲线资料，核磁孔隙度测井资料可以有针对的测得，上述测井曲线资料和核磁孔隙度测井资料常可以直接获得。目标区域的取心段砂岩的样品，通常是对取心段砂岩样品要进行取样获得。

步骤102：对所述样品进行铸体薄片及扫描电镜实验分析，确定所述样品的微观裂缝发育区。

在一个实施方式有，为了准确的确定所述样品的微观裂缝发育区，又考虑到微观裂缝的测井响应特征通常相对不明显。因此，通过同时对样品进行铸体薄片及扫描电镜分析确定样品的微观裂缝发育区。其中，通过铸体薄片可以观察到较大的微裂缝，通过扫描电镜可以观察到尺度更小级别的微裂缝。因此，将上述两种实验方式相结合可以更加准确地确定样品中的微观裂缝发育区，以便后续确定微观裂缝的测井响应特征。

步骤103：根据所述目标区域的测井曲线资料、所述目标区域的核磁孔隙度测井资料和所述样品的微观裂缝发育区，确定所述目标区域的微观裂缝的测井响应特征。

在一个实施方式中，通过分析微观裂缝发育区中的微观裂缝，可以确定所述微观裂缝的测井响应特征具体可以包括：纵波速度相对较小，核磁测井孔隙度相对较低。即纵波速度小于预设速度，和/或，核磁测井孔隙度低于预设孔隙度。

在一个实施方式中，为了确定目标区域的微观裂缝的测井响应特征，具体可以按照以下步骤执行(S1至S3)。

S1：根据所述目标区域的测井曲线资料中的声波时差曲线，确定所述目标区域的纵波速度。

S2：根据所述目标区域的核磁孔隙度测井资料，确定所述目标区域的核磁测井孔隙度。

S3：根据所述样品的微观裂缝发育区，确定所述预设速度和所述预设孔隙度。

进而可以根据上述获得的纵波速度、核磁测井孔隙度、预设速度和所述预设孔隙度，确定各个纵波速度和核磁测井孔隙度是否满足微观裂缝的测井响应特征，即是否满足纵波速度小于预设速度，和/或，核磁测井孔隙度低于预设孔隙度。

步骤104：将所述目标区域按照深度分成多个连续的小层，并根据所述微观裂缝的测井响应特征，为所述多个连续的小层中的各个小层分别建立岩石物理量版。

在一个实施方式中，为了保证结果的精度和后续确定微观裂缝测井特征曲线的准确性和连续性，可以按照预设的深度将目标区域按深度依次分为多个小层，分别进行处理。例如，目标区域的地层厚度为100m，可以按照每小层10m的预设厚度，将该地层平均分成10个小层。当然，也可以根据具体情况，将100m后的地层分成不等厚度的5个小层，其中，第一小层10m，第二小层20m，第三小层20m，第四小层40，第五小层10m。只要保证后续获得的各个小层的岩石物理量版连续准确即可。对此，本申请不作限定。

需要说明的是，为了使得后续获得岩石物理量版尽可能的准确、连续，具体实施时，还可以根据具体情况或者实施条件，将目标区域多分几个小层，以使得每个小层的厚度尽可能的小，即可以将后续获得岩石物理量版的深度范围减小，这样可以进一步减小后续确定的微裂缝测井特征曲线的误差。

在一个实施方式中，为了确定各个小层的岩石物理量版，具体实施时，可以：根据所述微观裂缝的测井响应特征，对所述多个连续的小层中的各个小层，分别以所在小层的核磁测井孔隙度为横坐标，以纵波速度为纵坐标，建立各个小层的岩石物理量版。

步骤105：根据多个连续的小层中各个小层的岩石物理量版，确定微观裂缝测井特征曲线。

在一个实施方式中，为了确定微观裂缝测井特征曲线，具体可以按照以下步骤执行(S1至S2)。

S1：从所述多个连续的小层中各个小层的岩石物理量版中逐小层获取微观裂缝的信息。

S2：将获取的各个小层的微观裂缝的信息，作为所述微观裂缝测井特征曲线。

在一个实施方式中，为了获取各个岩石物理量版中的微观裂缝的信息，具体实施时，可以通过分析落在各个岩石物理量版中左下角区域的为微观裂缝发育区域，获取该岩石物理量版中的微观裂缝的信息。因为，一般情况下，由于微观裂缝的敏感曲线特征，所以在理想状态下，岩石物理量版的微观裂缝信息主要会分布在该岩石物理量版的左下角。当然，不同研究区块由于沉积环境、成岩作用等的差异，微观裂缝数据不一定都落在左下角，具体实施时，根据实际情况进行相应调整。

步骤106：对所述目标区域的二维叠后地震资料进行相干属性分析，得到二维叠后相干地震资料。

在一个实施方式中，所述的二维叠后地震资料一般只能用于研究断裂的分布，而由于微观尺度微裂缝受普通二维叠后地震资料分辨率的限制，很难被预测出来。因此，为了后续便于研究微观裂缝在目标区域内的分布特征，具体实施时，可以先对目标区域的二维叠后地震资料进行相干属性分析，得到二维叠后相干地震资料，用以后续确定微观裂缝的分布特征。

步骤107：以所述二维叠后相干地震资料为地震波形指示反演的数据体，根据所述微观裂缝测井特征曲线，通过地震波形指示反演，得到所述目标区域的微观裂缝的分布特征。

在一个实施方式中，为了准确地确定目标区域的微观裂缝的分布特征，具体实施时，可以按照以下步骤(S1至S2)执行。

S1：以所述二维叠后相干地震资料为地震波形指示反演的数据体，并以所述微观裂缝测井特征曲线作为约束条件，优选高频成分，以建立反演模型。

S2：根据所述反演模型，进行地震波形指示反演，得到所述目标区域的微观裂缝的分布特征。

在一个实施方式中，为了建立准确的反演模型，具体实施时可以按照以下步骤(S1至S4)实施。

S1：以所述二维叠后相干地震资料为地震波形指示反演的数据体，并以所述微观裂缝测井特征曲线作为约束条件，优选高频成分，以建立初始模型，即寻找微观裂缝与相干属性波形关系。

S2：根据所述初始模型，反演得到所述目标区域过井的二维叠后地震剖面上微观裂缝空间分布特征。

S3：将所述过井二维叠后地震剖面上微观裂缝空间分布特征与所述目标区域的井上的微观裂缝测井特征曲线进行比较，得到比较误差。

S4：如果所述比较误差小于预设阈值，则将所述初始模型作为所述反演模型。

在一个实施方式中，为了对比较误差大于等于预设阈值的初始模型进行校正。在将所述过井二维叠后地震剖面上微观裂缝空间分布特征与所述目标区域的井上的微观裂缝测井特征曲线进行比较，得到比较误差之后，所述方法还可以包括：如果所述比较误差大于等于所述预设阈值，则重新选择高频成分，并根据重新选择的高频成分建立初始模型，直到根据建立的初始模型，反演得到的所述目标区域的过井二维叠后地震剖面上的微观裂缝空间分布特征与所述目标区域的井上的微观裂缝测井特征曲线之间的比较误差小于所述预设阈值，则将当前的初始模型作为所述反演模型。

在一个实施方式中，为了校正初始模型，所述重新选择高频成分具体可以包括：调整初始模型中高频成分的参数数值，例如最大频率、主频等参数。

在本申请实施例中，通过利用铸体薄片和扫描电镜确定微观裂缝发育区，参考测井曲线和核磁测井孔隙度曲线确定微观裂缝的测井响应特征，结合岩石物理量板进而确定微观裂缝测井特征曲线；再对二维叠后地震资料进行相干属性分析，使得可以利用二维叠后相干地震资料反应裂缝的特征；最后根据微观裂缝测井特征曲线和二维叠后相干地震资料，通过地震波形指示反演，确定微观裂缝分布特征，从而解决了现有的微观裂缝确定方法中存在的使用范围有限和微观裂缝确定精度低的技术问题，达到了在二维剖面上准确确定微观裂缝分布的技术效果。

基于同一发明构思，本发明实施例中还提供了一种基于二维叠后地震资料微观裂缝表征确定装置，如下面的实施例所述。由于基于二维叠后地震资料微观裂缝表征确定装置解决问题的原理与基于二维叠后地震资料微观裂缝表征确定方法相似，因此基于二维叠后地震资料微观裂缝表征确定装置的实施可以参见基于二维叠后地震资料微观裂缝表征确定方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。请参阅图2，是本发明实施例的基于二维叠后地震资料微观裂缝表征确定装置的一种组成结构图，该装置可以包括：获取模块201、第一确定模块202、第二确定模块203、第三确定模块204、第四确定模块205、相干属性分析模块206和第五确定模块207。下面对该结构进行具体说明。

获取模块201，具体可以用于获取目标区域的取心段砂岩的样品、所述目标区域的测井曲线资料、所述目标区域的核磁孔隙度测井资料和所述目标区域的二维叠后地震资料。

第一确定模块202，具体可以用于对所述样品进行铸体薄片及扫描电镜实验分析，确定所述样品的微观裂缝发育区。

第二确定模块203，具体可以用于根据所述目标区域的测井曲线资料、所述目标区域的核磁孔隙度测井资料和所述样品的微观裂缝发育区，确定所述目标区域的微观裂缝的测井响应特征。

第三确定模块204，具体可以用于将所述目标区域按照深度分成多个连续的小层，并根据所述微观裂缝的测井响应特征，为所述多个连续的小层中的各个小层分别建立岩石物理量版。

第四确定模块205，具体可以用于根据多个连续的小层中各个小层的岩石物理量版，确定微观裂缝测井特征曲线。

相干属性分析模块206，具体可以用于对所述目标区域的二维叠后地震资料进行相干属性分析，得到二维叠后相干地震资料。

第五确定模块207，具体可以用于以所述二维叠后相干地震资料为地震波形指示反演的数据体，根据所述微观裂缝测井特征曲线，通过地震波形指示反演，得到所述目标区域的微观裂缝的分布特征。

在一个实施方式中，所述微观裂缝的测井响应特征具体可以包括：纵波速度小于预设速度，和/或，核磁测井孔隙度低于预设孔隙度。

在一个实施方式中，为了确定微观裂缝的测井响应特征，所述第二确定模块203具体可以包括：

第一确定单元，用于根据所述目标区域的测井曲线资料中的声波时差曲线，确定所述目标区域的纵波速度；

第二确定单元，用于根据所述目标区域的核磁孔隙度测井资料，确定所述目标区域的核磁测井孔隙度；

第三确定单元，用于根据所述样品的微观裂缝发育区，确定所述预设速度和所述预设孔隙度。

在一个实施方式中，为了确定岩石物理量版，所述第三确定模块204具体可以用于根据所述微观裂缝的测井响应特征，对所述多个连续的小层中的各个小层，分别以所在小层的核磁测井孔隙度为横坐标，以纵波速度为纵坐标，建立各个小层的岩石物理量版。

在一个实施方式中，为了确定微观裂缝测井特征曲线，所述第四确定模块205具体可以包括：

信息获取单元，用于从所述多个连续的小层中各个小层的岩石物理量版中逐小层获取微观裂缝的信息；

第四确定单元，用于将获取的各个小层的微观裂缝的信息，作为所述微观裂缝测井特征曲线。

在一个实施方式中，为了使得获取的二维叠后地震资料可以反映微观裂缝分布特征，因此具体实施时，要使用相干分析模块206对所述目标区域的二维叠后地震资料进行相干属性分析，得到更好的反应断裂属性的二维叠后相干地震资料。

在一个实施方式中，为了确定目标区域的微观裂缝的分布特征，所述第五确定模块207具体可以包括：

模型确定单元，用于以所述二维叠后相干地震资料为地震波形指示反演的数据体，并以所述微观裂缝测井特征曲线作为约束条件，优选高频成分，以建立反演模型；

第五确定单元，用于根据所述反演模型，进行地震波形指示反演，得到所述目标区域的微观裂缝的分布特征。

需要说明的是，上述模型确定单元为了建立精度好的反演模型，具体实施时，可以先以所述二维叠后相干地震资料为地震波形指示反演的数据体，并以所述微观裂缝测井特征曲线作为约束条件，优选高频成分，以建立初始模型；再根据所述初始模型，反演得到所述目标区域过井的二维叠后地震剖面上微观裂缝空间分布特征；进而将所述过井二维叠后地震剖面上微观裂缝空间分布特征与所述目标区域的井上的微观裂缝测井特征曲线进行比较，得到比较误差；这时，如果所述比较误差小于预设阈值，则将所述初始模型作为所述反演模型；如果所述比较误差大于等于所述预设阈值，则重新选择高频成分，并根据重新选择的高频成分建立初始模型，直到根据建立的初始模型，反演得到的所述目标区域的过井二维叠后地震剖面上的微观裂缝空间分布特征与所述目标区域的井上的微观裂缝测井特征曲线之间的比较误差小于所述预设阈值，则将当前的初始模型作为所述反演模型。

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

需要说明的是，上述实施方式阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。为了描述的方便，在本说明书中，描述以上装置时以功能分为各种单元分别描述。当然，在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。

此外，在本说明书中，诸如第一和第二这样的形容词仅可以用于将一个元素或动作与另一元素或动作进行区分，而不必要求或暗示任何实际的这种关系或顺序。在环境允许的情况下，参照元素或部件或步骤(等)不应解释为局限于仅元素、部件、或步骤中的一个，而可以是元素、部件、或步骤中的一个或多个等。

从以上的描述中，可以看出，本申请实施例提供的基于二维叠后地震资料微观裂缝表征确定方法和装置。通过利用铸体薄片和扫描电镜确定微观裂缝发育区，参考测井曲线和核磁测井孔隙度曲线确定微观裂缝的测井响应特征，结合岩石物理量板进而确定微观裂缝测井特征曲线；再对二维叠后地震资料进行相干属性分析，使得可以利用二维叠后相干地震资料反应微观裂缝的特征；最后根据微观裂缝测井特征曲线和二维叠后相干地震资料，通过地震波形指示反演，确定微观裂缝分布特征，从而解决了现有的微观裂缝确定方法中存在的使用范围有限和微观裂缝确定精度低的技术问题，达到了在二维剖面上准确确定微观裂缝分布的技术效果；又通过利用铸体薄片和扫描电镜对样品进行分析，再根据样品的分析结果、测井曲线资料和核磁孔隙度测井资料确定微观裂缝的测井响应特征，进一步提高了微观裂缝确定的精度；还通过对二维叠后地震资料进行相干属性分析，使得获得的二维叠后相干地震资料可以反应微观裂缝的特征，扩大了微观裂缝的确定范围；又通过地震波形指示反演确定微观裂缝的分布特征，进一步改善了确定微观裂缝的精度。

在一个具体实施场景中，应用本申请实施例提供基于二维叠后地震资料微观裂缝表征确定方法/装置对某地的微观裂缝进行预测。具体实施时，可以按照以下几个步骤(步骤1至步骤3)执行。

步骤1、确定微观裂缝的测井响应特征。

首先要求探井必须要有常规测井曲线资料及核磁孔隙度测井资料。其次，对取心段砂岩样品要进行取样，进行铸体薄片及扫描电镜分析。建立测井曲线及铸体薄片和扫描电镜的综合分析柱状图，具体可以参阅附图3。从图中，可以发现，微米尺度的微观裂缝对纵波速度及核磁孔隙度反应较为敏感，具体地，在微观裂缝发育的地方，纵波速度较小，孔隙度较低。图3即为L1井测井资料及铸体和电镜照片综合柱状图，通过对该图的分析，发现微米尺度微观裂缝发育地方，纵波速度较小，测井孔隙度较低。

步骤2、放大微观裂缝测井信号，以确定微观裂缝测井特征曲线。

结合微观裂缝对纵波速度和核磁孔隙度的测井响应特征，分小层逐层建立以核磁测井孔隙度为横坐标，以纵波速度为纵坐标的岩石物理量版。从量版中可以发现，微观尺度微观裂缝主要分布在岩石物理量版的左下角。从每个小层的岩石物理量版上读取微观尺度微观裂缝信息，进而可得到连续的微观裂缝测井特征曲线。具体可以参阅图4的致密砂岩储层微观裂缝岩石物理判别量版，结合微观裂缝的测井响应特征，可以认为微观裂缝主要发育在岩石物理量版的左下角。

步骤3.测井联合地震波形指示反演微观裂缝在二维叠后地震资料上的确定微观裂缝的分布特征(即确定微观裂缝的分布特征)。

首先，对过井的二维叠后地震资料进行相干属性分析。

其次，以步骤2中得到的微观裂缝特征曲线为约束条件，优选出合适的高频成分，建立合理的初始模型(即寻找微观裂缝与相干属性波形关系)，再运用波形指示反演预测连井剖面上微观裂缝空间分布特征。并与井上微观裂缝特征曲线做对比，如果波形反演结果与井上微观裂缝特征曲线的吻合率较低，则重新建立初始模型，直到反演效果较好为止。具体可以参阅图5的过L1井的Crossline(道)方向上的原始二维叠后地震资料。进一步，具体可以参阅图6的过L1井Crossline方向上的微观裂缝地震波形指示反演剖面，图中井旁边曲线为微观裂缝测井特征曲线。从该过井反演剖面可以看到反演所得的微观裂缝与井上统计结果(井柱子旁边的特征测井曲线)基本保持一致，说明结合岩石物理量版建立的微观裂缝测井特征曲线，通过波形指示反演，可以较为准确地表征微观裂缝在空间的发育程度。

通过上述具体应用本申请实施例提供的基于二维叠后地震资料微观裂缝表征确定方法和装置的示例，验证了通过本申请实施例提供的基于二维叠后地震资料微观裂缝表征确定方法和装置确实可以解决了现有的微观裂缝确定方法中存在的使用范围有限和微观裂缝确定精度低的技术问题，并且达到了在二维叠后地震剖面上准确确定微观裂缝分布的技术效果。

尽管本申请内容中提到不同的微观裂缝确定方法或装置，但是，本申请并不局限于必须是行业标准或实施例所描述的情况等，某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取、处理、输出、判断方式等的实施例，仍然可以属于本申请的可选实施方案范围之内。

虽然本申请提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。

上述实施例阐明的装置或模块等，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述模块的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个模块或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。

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

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

通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，移动终端，服务器，或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。

本说明书中的各个实施例采用递进的方式描述，各个实施例之间相同或相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。本申请可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。

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

Claims (12)

1.一种基于二维叠后地震资料微观裂缝表征确定方法，其特征在于，包括：

获取目标区域的取心段砂岩的样品、所述目标区域的测井曲线资料、所述目标区域的核磁孔隙度测井资料和所述目标区域的二维叠后地震资料；

对所述样品进行铸体薄片及扫描电镜实验分析，确定所述样品的微观裂缝发育区；

根据所述目标区域的测井曲线资料、所述目标区域的核磁孔隙度测井资料和所述样品的微观裂缝发育区，确定所述目标区域的微观裂缝的测井响应特征；

将所述目标区域按照深度分成多个连续的小层，并根据所述微观裂缝的测井响应特征，为所述多个连续的小层中的各个小层分别建立岩石物理量版；

根据多个小层中各个小层的岩石物理量版，确定微观裂缝测井特征曲线；

对所述目标区域的二维叠后地震资料进行相干属性分析，得到二维叠后相干地震资料；

以所述二维叠后相干地震资料为地震波形指示反演的数据体，根据所述微观裂缝测井特征曲线，通过地震波形指示反演，得到所述目标区域的微观裂缝的分布特征。

2.根据权利要求1所述的方法，其特征在于，

所述微观裂缝的测井响应特征包括：纵波速度小于预设速度，和/或，核磁测井孔隙度低于预设孔隙度。

3.根据权利要求2所述的方法，其特征在于，根据所述目标区域的测井曲线资料、所述目标区域的核磁孔隙度测井资料和所述样品的微观裂缝发育区，确定所述目标区域的微观裂缝的测井响应特征，包括：

根据所述目标区域的测井曲线资料中的声波时差曲线，确定所述目标区域的纵波速度；

根据所述目标区域的核磁孔隙度测井资料，确定所述目标区域的核磁测井孔隙度；

根据所述样品的微观裂缝发育区，确定所述预设速度和所述预设孔隙度。

4.根据权利要求1所述的方法，其特征在于，根据所述微观裂缝的测井响应特征，为所述多个连续的小层中的各个小层分别建立岩石物理量版，包括：

根据所述微观裂缝的测井响应特征，对所述多个连续的小层中的各个小层，分别以所在小层的核磁测井孔隙度为横坐标，以纵波速度为纵坐标，建立各个小层的岩石物理量版。

5.根据权利要求1所述的方法，其特征在于，根据多个连续的小层中各个小层的岩石物理量版，确定微观裂缝连续的测井特征曲线，包括：

从所述多个连续的小层中各个小层的岩石物理量版中逐小层获取微观裂缝的信息；

将获取的各个小层的微观裂缝的信息，作为所述微观裂缝测井特征曲线。

6.根据权利要求1所述的方法，其特征在于，以所述二维叠后相干地震资料为地震波形指示反演的数据体，根据所述微观裂缝测井特征曲线，通过地震波形指示反演，得到所述目标区域的微观裂缝的分布特征，包括：

以所述二维叠后相干地震资料为地震波形指示反演的数据体，并以所述微观裂缝测井特征曲线作为约束条件，优选高频成分，以建立反演模型；

根据所述反演模型，进行地震波形指示反演，得到所述目标区域的微观裂缝的分布特征。

7.根据权利要求6所述的方法，其特征在于，以所述二维叠后相干地震资料为地震波形指示反演的数据体，并以所述微观裂缝测井特征曲线作为约束条件，优选高频成分，以建立反演模型，包括：

以所述二维叠后相干地震资料为地震波形指示反演的数据体，并以所述微观裂缝测井特征曲线作为约束条件，优选高频成分，以建立初始模型；

根据所述初始模型，反演得到所述目标区域过井的二维叠后地震剖面上微观裂缝空间分布特征；

将所述过井二维叠后地震剖面上微观裂缝空间分布特征与所述目标区域的井上的微观裂缝测井特征曲线进行比较，得到比较误差；

如果所述比较误差小于预设阈值，则将所述初始模型作为所述反演模型。

8.根据权利要求7所述的方法，其特征在于，在将所述过井二维叠后地震剖面上微观裂缝空间分布特征与所述目标区域的井上的微观裂缝测井特征曲线进行比较，得到比较误差之后，所述方法还包括：

如果所述比较误差大于等于所述预设阈值，则重新选择高频成分，并根据重新选择的高频成分建立初始模型，直到根据建立的初始模型，反演得到的所述目标区域的过井二维叠后地震剖面上的微观裂缝空间分布特征与所述目标区域的井上的微观裂缝测井特征曲线之间的比较误差小于所述预设阈值，则将当前的初始模型作为所述反演模型。

9.一种基于二维叠后地震资料微观裂缝表征确定装置，其特征在于，包括：

获取模块，用于获取目标区域的取心段砂岩的样品、所述目标区域的测井曲线资料、所述目标区域的核磁孔隙度测井资料和所述目标区域的二维叠后地震资料；

第一确定模块，用于对所述样品进行铸体薄片及扫描电镜实验分析，确定所述样品的微观裂缝发育区；

第二确定模块，用于根据所述目标区域的测井曲线资料、所述目标区域的核磁孔隙度测井资料和所述样品的微观裂缝发育区，确定所述目标区域的微观裂缝的测井响应特征；

第三确定模块，用于将所述目标区域按照深度分成多个连续的小层，并根据所述微观裂缝的测井响应特征，为所述多个连续的小层中的各个小层分别建立岩石物理量版；

第四确定模块，用于根据多个连续的小层中各个小层的岩石物理量版，确定微观裂缝测井特征曲线；

相干属性分析模块，用于对所述目标区域的二维叠后地震资料进行相干属性分析，得到二维叠后相干地震资料；

第五确定模块，用于以所述二维叠后相干地震资料为地震波形指示反演的数据体，根据所述微观裂缝测井特征曲线，通过地震波形指示反演，得到所述目标区域的微观裂缝的分布特征。

10.根据权利要求9所述的装置，其特征在于，所述第二确定模块包括：

第一确定单元，用于根据所述目标区域的测井曲线资料中的声波时差曲线，确定所述目标区域的纵波速度；

第二确定单元，用于根据所述目标区域的核磁孔隙度测井资料，确定所述目标区域的核磁测井孔隙度；

第三确定单元，用于根据所述样品的微观裂缝发育区，确定所述预设速度和所述预设孔隙度。

11.根据权利要求9所述的装置，其特征在于，所述第四确定模块包括：

信息获取单元，用于从所述多个连续的小层中各个小层的岩石物理量版中逐小层获取微观裂缝的信息；

第四确定单元，用于将获取的各个小层的微观裂缝的信息，作为所述微观裂缝测井特征曲线。

12.根据权利要求9所述的装置，其特征在于，所述第五确定模块包括：

模型确定单元，用于以所述二维叠后相干地震资料为地震波形指示反演的数据体，并以所述微观裂缝测井特征曲线作为约束条件，优选高频成分，以建立反演模型；

第五确定单元，用于根据所述反演模型，进行地震波形指示反演，得到所述目标区域的微观裂缝的分布特征。