微观溶蚀孔隙地球物理二维表征确定方法和装置与流程

本申请涉及油气勘探领域，特别涉及一种微观溶蚀孔隙地球物理二维表征确定方法和装置。

背景技术：

随着油气勘探技术的不断提高，发现深层及低渗透-超低渗透致密砂岩储层中探明的油气率越来越高，因此越来越受到重视。但是，这类储层的勘探难度也很大。其中，上述低渗透-特低渗透致密砂岩储层通常是指满足孔隙度小于10％，且渗透率小于1md的砂岩储层。

具体地，深层及低渗透-特低渗透致密砂岩储层一般具有较大的资源潜力，同时也具有较差的储集性能。研究发现，次生微观溶蚀孔隙的发育，能改善深层及低渗透-特低渗透致密砂岩储层的储集性能。但是由于受地震资料品质的影响，目前基于叠后二维地震资料去表征微观溶蚀孔隙的发育特征的研究很少。利用现有的方法，在测井资料上发现次生微观溶蚀孔隙的测井响应特征就很困难，在叠后二维地震资料上发现微观溶蚀孔隙的地震响应特征就更加困难。因此，现有方法具体实施时，往往存在不能准确确定目标区域中的微观溶蚀孔隙的分布特征的技术问题。

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

技术实现要素：

本申请实施例提供了一种微观溶蚀孔隙地球物理二维表征确定方法和装置，以解决现有的微观溶蚀孔隙确定方法存在的不能准确确定微观溶蚀孔隙在二维剖面上分布特征的技术问题。

本申请实施例提供了一种微观溶蚀孔隙地球物理二维表征确定方法，包括：

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

对所述样品进行铸体薄片及扫描电镜实验分析，得到分析结果；

根据所述分析结果和所述测井曲线资料，通过岩心归位，确定所述样品的微观溶蚀孔隙发育区；

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

基于相同的发明构思，本申请实施例还提供了一种微观溶蚀孔隙地球物理二维表征确定装置，包括：

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

分析模块，用于对所述样品进行铸体薄片及扫描电镜实验分析，得到分析结果；

第一确定模块，用于根据所述分析结果和所述测井曲线资料，通过岩心归位，确定所述样品的微观溶蚀孔隙发育区；

第二确定模块，用于根据所述分析结果和所述测井曲线资料，通过岩心归位，确定所述样品的微观溶蚀孔隙发育区；

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

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

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

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

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

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

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

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

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

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

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

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

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

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

附图说明

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

图1是根据本申请实施例的一种微观溶蚀孔隙地球物理二维表征确定方法的处理流程图；

图2是根据本申请实施例的一种微观溶蚀孔隙地球物理二维表征确定装置的组成结构图；

图3是应用本申请实施例提供一种微观溶蚀孔隙地球物理二维表征确定方法/装置得到的目标区域中W1井测井微观溶蚀孔隙发育特征综合柱状图；

图4是应用本申请实施例提供一种微观溶蚀孔隙地球物理二维表征确定方法/装置得到的致密砂岩储层微观溶蚀孔隙岩石物理判别量版的示意图；

图5是应用本申请实施例提供一种微观溶蚀孔隙地球物理二维表征确定方法/装置得到的过W1井Crossline(道)方向叠后二维地震资料的示意图；

图6是应用本申请实施例提供一种微观溶蚀孔隙地球物理二维表征确定方法/装置得到的过W1井Crossline方向微观溶蚀孔隙地震波形指示反演剖面的示意图。

具体实施方式

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

考虑到目前由于受地震资料品质的影响，基于叠后二维地震资料去表征微观溶蚀孔隙的发育特征的研究很少。由于微观溶蚀孔隙自身的特点，即：首先微观溶蚀孔隙尺度非常小，一般在常规测的井曲线上很难找到微观溶蚀孔隙的测井响应特征；其次，由于微观溶蚀孔隙成因较多，在常规二维叠后地震剖面上很难发现它的地震响应。因此，利用现有的方法，在测井资料上发现次生微观溶蚀孔隙的测井响应特征很困难，而在叠后二维地震资料上发现微观溶蚀孔隙的地震响应特征就更加困难。因此，现有方法在具体实施时，往往会存在不能准确确定目标区域中的微观溶蚀孔隙的分布特征的技术问题。针对产生上述技术问题的根本原因，本申请考虑可以先通过铸体薄片及扫描电镜，在岩心归位的基础上确定微观溶蚀孔隙发育区；其次结合测井曲线及核磁测井孔隙度曲线，分析微观溶蚀孔隙的测井响应特征，参考岩石物理量板，进而得到微观溶蚀孔隙测井特征曲线；最后根据微观溶蚀孔隙测井特征曲线和二维叠后地震资料，通过地震波形指示反演，确定微观溶蚀孔隙的分布特征。从而解决了现有的微观溶蚀孔隙确定方法存在的不能准确确定目标区域内的微观溶蚀孔隙分布的技术问题，达到了在二维剖面上准确确定微观溶蚀孔隙分布的技术效果。

基于上述思考思路，本申请提供了一种微观溶蚀孔隙地球物理二维表征确定方法。请参阅图1。本申请提供的一种微观溶蚀孔隙地球物理二维表征确定方法，具体可以包括以下步骤(步骤101至步骤107)。

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

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

步骤102：对所述样品进行铸体薄片及扫描电镜实验分析，得到分析结果。

在一个实施方式有，为了后续准确的确定所述样品的微观溶蚀孔隙发育区，又考虑到微观溶蚀孔隙的测井响应特征通常相对不明显，可以先确定样品的分析结果。即，可以先通过同时对样品进行铸体薄片及扫描电镜分析得到样品的铸体薄片及扫描电镜的分析结果，以便后续使用。

步骤103：根据所述分析结果和所述测井曲线资料，通过岩心归位，确定所述样品的微观溶蚀孔隙发育区。

在一个实施方式中，为了准确确定样品的微观溶蚀孔隙发育区，具体可以根据上述分析结果和测井资料，进行岩心归位；在岩心归位的基础上，确定样品的微观溶蚀孔隙发育区，以便后续确定微观溶蚀孔隙的测井响应特征。具体实施时，可以根据测井曲线资料和录井及取心段的岩心深度共同进行岩心归位。

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

步骤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，具体可以用于以所述二维叠后地震资料为地震波形指示反演的数据体，根据所述微观溶蚀孔隙测井特征曲线，通过地震波形指示反演，得到所述目标区域的微观溶蚀孔隙的分布特征。

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

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

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

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

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

需要说明的是，所述岩心归位，具体可以是指根据测井曲线资料和录井及取心段的岩心深度共同进行岩心归位。

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

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

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

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

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

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

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

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

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

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

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

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

在一个具体实施场景中，应用本申请实施例提供一种微观溶蚀孔隙地球物理二维表征确定方法/装置对某地的微观溶蚀孔隙进行预测。具体实施时，可以按照以下几个步骤(步骤1至步骤3)执行。

步骤1、确定微观溶蚀孔隙测井响应特征。

首先要求探井必须要有声波时差曲线,(即本申请中的测井曲线资料所包括的一种曲线资料)及核磁测井孔隙度资料(即本申请中的核磁孔隙度测井资料)。其次，对取心段砂岩样品要进行铸体薄片及扫描电镜分析。进而进行岩心归位，在岩心归位基础上，找出微观溶蚀孔隙发育深度范围(即本申请中的微观溶蚀孔隙发育区)。最后，结合纵波速度及核磁测井孔隙度曲线可以发现，微观溶蚀孔隙发育的地方，纵波速度较大，孔隙度较小(即得到本申请中微观溶蚀孔隙的测井响应特征)。具体可以参阅图3的W1井测井微观溶蚀孔隙发育特征综合柱状图。

步骤2、通过岩石物理量版，得到微观溶蚀孔隙测井特征曲线。

针对上述微观溶蚀孔隙发育范围内纵波速度及核磁孔隙度的测井响应特征(即本申请中的微观溶蚀孔隙的测井响应特征)，建立以核磁测井孔隙度为横坐标，以纵波速度为纵坐标的岩石物理量版，从量版中可以发现，微观溶蚀孔隙主要分布在岩石物理量版的左上角，即可以确定微观溶蚀孔隙具有较大的纵波速度和较小的核磁测井孔隙度。再分小层逐层建立这种岩石物理量版，从每个量版上读取微观溶蚀孔隙信息，进而可得到连续的微观溶蚀孔隙测井特征曲线。具体可以参阅图4的致密砂岩储层微观溶蚀孔隙岩石物理量版。

步骤3、通过测井联合地震波形指示反演，表征得到微观溶蚀孔隙在二维剖面上的分布特征。

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

通过上述具体应用本申请实施例提供的一种微观溶蚀孔隙地球物理二维表征确定方法和装置的示例，验证了通过本申请实施例提供的一种微观溶蚀孔隙地球物理二维表征确定方法和装置确实可以解决了现有的方法存在的不能准确确定目标区域的微观溶蚀孔隙在二维空间的分布特征的技术问题。

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

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

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

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

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

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

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

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