CN108873085B - 一种沉积岩薄层电磁数据处理方法及装置 - Google Patents

Abstract

本申请提供一种沉积岩薄层电磁数据处理方法及装置。所述方法包括：获取采集到的电磁数据对应的三维网格数据；从三维网格数据中获取待解释剖面的剖面网格数据，获取剖面网格数据对应的异常剖面网格数据；利用电测井数据对异常剖面网格数据进行标定，获得沉积岩薄层在异常剖面网格数据中对应的电阻率异常区域；获取电阻率异常区域对应的深度值和电阻率异常平均值；根据深度值和电阻率平均值对电阻率异常区域进行搜索，将小于电阻率平均值的网格点的厚度进行叠加，获得电阻率异常区域对应的厚度数据。利用本申请中各实施例实现了沉积岩薄层识别和解释，提高了电磁数据对于高阻电性地层内夹的低阻沉积岩薄层目标的反映能力。

Description

一种沉积岩薄层电磁数据处理方法及装置

技术领域

本说明书实施例属于地球物理勘探技术领域，尤其涉及一种沉积岩薄层电磁数据处理方法及装置。

背景技术

电磁勘探是地震勘探的一种方法，通过对电磁勘探获得的电磁数据进行处理可以获得勘探地区的地震解释数据，为后续油气的开发提供理论依据。

对大地电磁数据电阻率反演结果进行地层结构解释和低阻层追踪是目前的比较常规技术手段，常规技术是基于对反演数据的剖面直接进行解释，这种方法在低阻地层厚度较大(如：厚度大于300m以上)、电阻率差异大的地区一般可以取得一定效果。但是，对于火山岩盆地、碳酸盐岩发育地区，针对低阻沉积岩薄层(如：小于200m)的识别和定量描述则存在较大困难。约束反演方法可能也存在初始模型难以建立的突出问题，目前，对于电磁数据的常规处理和解释方法往往无法识别这种低阻沉积岩薄层，这是目前油气勘探给大地电磁勘探数据处理解释提出的新挑战。因此，业内亟需一种能够满足识别火山岩盆地或碳酸盐岩发育区等低阻薄层沉积岩地层解释和目标评价的实施方案。

发明内容

本说明书实施例目的在于提供一种沉积岩薄层电磁数据处理方法及装置，提高了电磁数据对于高阻电性地层内夹的低阻沉积岩薄层目标的反映能力，为火山岩盆地及碳酸盐岩发育区的区带评价及目标优选提供了新方法。

一方面本说明书实施例提供了一种沉积岩薄层电磁数据处理方法，包括：

获取采集到的电磁数据对应的三维网格数据；

从所述三维网格数据中获取待解释剖面的剖面网格数据，获取所述剖面网格数据对应的异常剖面网格数据；

利用电测井数据对所述异常剖面网格数据进行标定，获得沉积岩薄层在所述异常剖面网格数据中对应的电阻率异常区域；

确定所述电阻率异常区域对应的深度值和电阻率异常平均值；

根据所述深度值和所述电阻率平均值对所述电阻率异常区域进行搜索，将小于所述电阻率平均值的网格点的厚度进行叠加，获得所述电阻率异常区域对应的厚度数据。

进一步地，所述方法的另一个实施例中，

所述根据所述深度值和所述电阻率异常平均值对所述电阻率异常区域进行搜索，将小于所述电阻率平均值的网格点的厚度进行叠加，获得所述电阻率异常区域对应的厚度数据，包括：

从所述深度值对应的测试网格点沿所述电阻率异常区域向下搜索，若搜索到的网格点的电阻率异常值小于所述电阻率异常平均值，则将所述网格点对应的厚度进行累加，直至搜索到的网格点的电阻率异常值等于所述电阻率异常平均值，确定出所述测试网格点对应的厚度数据；

对所述电阻率异常区域中各网格点依次进行纵向搜索，获取所述电阻率异常区域中各网格点对应的厚度数据。

进一步地，所述方法的另一个实施例中，所述方法还包括：

获取所述三维网格数据中各待解释剖面对应的厚度数据；

将所述各待解释剖面对应的厚度数据进行剖面网格化，获得所述沉积岩薄层对应的地层厚度图。

进一步地，所述方法的另一个实施例中，所述获取采集到的电磁数据对应的三维网格数据，包括：

将获取到的电磁数据进行反演，获得三维电阻率数据体；

将所述三维电阻率数据体中的电阻率数值进行对数运算；

将对数运算后的三维电阻率数据体进行网格化，获得所述三维网格数据。

进一步地，所述方法的另一个实施例中，所述获取所述剖面网格数据对应的异常剖面网格数据，包括：

对所述剖面网格数据进行向上延拓计算，获得所述待解释剖面的背景场数据；

将所述剖面网格数据与所述背景场数据的差值作为所述异常剖面网格数据。

进一步地，所述方法的另一个实施例中，所述确定所述电阻率异常区域对应的深度值和电阻率异常平均值，包括：

沿所述待解释剖面横向搜索所述电阻率异常区域，获取所述电阻率异常区域对应的深度值和所述沉积岩薄层底界向上预设厚度内的电阻率异常平均值。

另一方面，本说明书实施例还提供了一种沉积岩薄层电磁数据处理装置，包括：

三维网格数据获取模块，用于获取采集到的电磁数据对应的三维网格数据；

异常剖面网格获取模块，用于从所述三维网格数据中获取待解释剖面的剖面网格数据，获取所述剖面网格数据对应的异常剖面网格数据；

异常区域确定模块，用于利用电测井数据对所述异常剖面网格数据进行标定，获得沉积岩薄层在所述异常剖面网格数据中对应的电阻率异常区域；

电阻率平均值处理模块，用于确定所述电阻率异常区域对应的深度值和电阻率异常平均值；

网格搜索模块，用于根据所述深度值和所述电阻率平均值对所述电阻率异常区域进行搜索，将小于所述电阻率平均值的网格点的厚度进行叠加，获得所述电阻率异常区域对应的厚度数据。

进一步地，所述装置的另一个实施例中，所述网格搜索模块具体用于：

从所述深度值对应的测试网格点沿所述电阻率异常区域向下搜索，若搜索到的网格点的电阻率异常值小于所述电阻率异常平均值，则将所述网格点对应的厚度进行累加，直至搜索到的网格点的电阻率异常值等于所述电阻率异常平均值，确定出所述测试网格点对应的厚度数据；

对所述电阻率异常区域中各网格点依次进行纵向搜索，获取所述电阻率异常区域中各网格点对应的厚度数据。

进一步地，所述装置的另一个实施例中，所述装置还包括平面处理模块用于：

获取所述三维网格数据中各待解释剖面对应的厚度数据；

将所述各待解释剖面对应的厚度数据进行剖面网格化，获得所述沉积岩薄层对应的地层厚度图。

进一步地，所述装置的另一个实施例中，所述三维网格数据获取模块具体用于：

将获取到的电磁数据进行反演，获得三维电阻率数据体；

将所述三维电阻率数据体中的电阻率数值进行对数运算；

将对数运算后的三维电阻率数据体进行网格化，获得所述三维网格数据。

进一步地，所述装置的另一个实施例中，所述异常剖面网格获取模块具体用于：

对所述剖面网格数据进行向上延拓计算，获得所述待解释剖面的背景场数据；

将所述剖面网格数据与所述背景场数据的差值作为所述异常剖面网格数据。

进一步地，所述装置的另一个实施例中，所述电阻率平均值处理模块具体用于：

沿所述待解释剖面横向搜索所述电阻率异常区域，获取所述电阻率异常区域对应的深度值和所述沉积岩薄层底界向上预设厚度内的电阻率异常平均值。

再一方面，本说明书实施例提供了一种计算机存储介质，其上存储有计算机程序，所述计算机程序被执行时，实现上述沉积岩薄层电磁数据处理方法。

还一方面，本说明书实施例提供了沉积岩薄层电磁数据处理系统，包括处理器以及用于存储处理器可执行指令的存储器，所述处理器执行所述指令时实现如上述沉积岩薄层电磁数据处理方法。

本说明书实施例提供的沉积岩薄层电磁数据处理方法及装置，通过对电磁数据进行三维网格化处理，根据获得的三维网格数据获得待解释剖面的异常剖面网格数据。对异常剖面网格数据进行标定，获得沉积岩薄层在剖面网格中的深度以及对应的电阻率异常区域。对获得的电阻率异常区域进行追踪搜索，获得待解释剖面的厚度数据，获得了低阻沉积岩薄层的厚度信息特征。实现了沉积岩薄层的识别和地层解释，提高了电磁数据对于高阻电性地层内夹的低阻沉积岩薄层目标的反映能力，为火山岩盆地及碳酸盐岩发育区的区带评价及目标优选提供了新方法，克服了约束反演依赖于初始模型的缺点。

附图说明

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

图1是本说明书一个实施例中的一种沉积岩薄层电磁数据处理方法流程示意图；

图2是本说明书又一个实施例中沉积岩薄层电磁数据的处理方法流程示意图；

图3是本说明书一个实施例中三维电阻率数据体的立体示意图；

图4是本说明书一个实施例中获得的沉积岩薄层地层厚度图的示意图；

图5是本说明书提供的沉积岩薄层电磁数据处理装置一个实施例的模块结构示意图；

图6是本说明书提供的一种沉积岩薄层电磁数据处理系统实施例的模块结构示意图。

具体实施方式

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

在地震勘探发生产过程中，利用不同的勘探方法可以获得相应的地震数据，对地震数据进行处理，可以得到勘探区域的地质结构解释等信息，为油气的开采提供理论依据。电磁勘探是一种比较常见的地震勘探方法，可以对电磁勘探中采集到的电磁数据进行处理，获得勘探区域的地层解释。

本申请实施例提出沉积岩薄层电磁数据处理方法，通过对采集到的电磁数据进行三维网格化处理，对三维网格数据进行解释追踪，可以获得沉积岩薄层的厚度信息特征，实现了沉积岩薄层的识别和地层解释，提高了电磁数据对低阻沉积岩薄层的反映能力，获取低阻沉积岩地层分布和厚度信息，为区带评价和目标优选提供资料依据。

具体地，图1是本申请一个实施例中的一种沉积岩薄层电磁数据处理方法流程示意图，本申请提供的沉积岩薄层电磁数据处理方法包括：

S1、获取采集到的电磁数据对应的三维网格数据。

利用电磁勘探可以采集到电磁数据，对电磁数据进行处理可以获得对应的三维网格数据。三维网格数据可以包括勘探区地下三维空间内各个节点对应的信息如：电阻率值、深度等。可以通过对采集到的电磁数据进行网格化处理，获取到三维网格数据。

本说明书一个实施例中，所述获取采集到的电磁数据对应的三维网格数据，可以包括：

将获取到的电磁数据进行反演，获得三维电阻率数据体；

将所述三维电阻率数据体中的电阻率数值进行对数运算；

将对数运算后的三维电阻率数据体进行网格化，获得所述三维网格数据。

在具体的实施过程中，可以对野外采集的三维大地电磁数据进行三维电阻率反演，在进行三维电阻率反演时，可以将地下空间剖分在低阻薄层发育深度段的深度剖分间距设置不大于薄层厚度，即可以按工区地质情况和勘探需要，对沉积岩薄层的厚度进行定义，比如不大于150m。可以采用共轭梯度法反演，获得三维电阻率数据体。三维电阻率反演是一种比较常用的大地电磁资料的反演方法，本说明书实施例中三维电阻率反演的具体过程可以包括：加载三维大地电磁资料数据、剖分地下空间、共轭梯度法反演地下空间电阻率分布、反演终止、输出反演结果，获得三维电阻率数据体。三维电阻率数据体可以包括地下空间各节点的电阻率值，具体可以包括每个节点的空间坐标(x，y，z)和节点处的电阻率值。

对获得的三维电阻率数据体的电阻率数值取对数运算，将对数运算后的三维电阻率数据体进行三维网格化，获得对数电阻率的三维网格数据。三维网格化可以表示对不规则空间分布的数据网格化为规则性网格(长、宽、高等距)的节点数据。取对数运算可以降低采集到的电磁数据的数值变化率，使得网格数据比较平滑自然，方便后续对数据进行追踪解释。

S2、从所述三维网格数据中获取待解释剖面的剖面网格数据，获取所述剖面网格数据对应的异常剖面网格数据。

获得三维网格数据后，可以从三维网格数据中指定需要解释的剖面即待解释剖面的剖面网格数据，如可以从三维网格数据中提取出待解释剖面的剖面网格数据。待解释剖面可以根据实际需要进行指定，也可以依次将各个剖面进行解释。获得待解释剖面的剖面网格数据后，可以对剖面网格数据进行处理，获得剖面网格数据中的异常剖面网格数据，异常剖面网格数据可以表示采集的电磁数据中有波动的数据，这些波动可以表征地下构造特征。可以通过去除背景场数据获得异常剖面网格数据，背景场也叫区域背景场，是场中的平缓变化部分或缓慢变化成分。

本申请一个实施例中，可以采用下述方法获得所述剖面网格数据对应的异常剖面网格数据：

对所述剖面网格数据进行向上延拓计算，获得所述待解释剖面的背景场数据；

将所述剖面网格数据与所述背景场数据的差值作为所述异常剖面网格数据。

在具体的实施过程中，向上延拓是位场计算的一种方法，可以表示计算高度大于已知场高度的位场的过程，可以通过向上延拓计算获得背景场。可以对获得的待解释剖面的剖面网格数据进行向上延拓计算，向上延拓的距离可以根据实际需要进行设置，本申请实施例不作具体限定。如：可以向上延拓500米，可以将向上延拓500米的计算结果作为待解释剖面的背景场数据。可以根据剖面网格数据和背景场数据计算出异常剖面网格数据，如：可以将获得的剖面网格数据与背景场数据的差值作为所述异常剖面网格数据。通过向上延拓计算，确定出待解释剖面的背景场数据和异常剖面网格数据，方法简单，为后续沉积岩薄层的识别以及厚度数据的获取奠定了理论基础。

S3、利用电测井数据对所述异常剖面网格数据进行标定，获得沉积岩薄层在所述异常剖面网格数据中对应的电阻率异常区域。

获得待解释剖面的异常剖面网格数据后，可以利用电测井数据对异常剖面网格数据进行标定，通过电测井数据可以确定出大地电磁剖面数据对应的地质含义、层位、构造、深度等，实现异常剖面网格数据的标定。根据标定后的异常剖面网格数据可以获得异常剖面网格数据中各个节点表示的地质含义，确定出沉积岩薄层在异常剖面网格数据中对应的电阻率异常区域。如：可以以电测井数据的沉积岩薄层层段底界深度对应的剖面网格数据的深度为标志，确定沉积岩薄层在剖面网格上的深度和对应的电阻率异常区域。

S4、确定所述电阻率异常区域对应的深度值和电阻率异常平均值。

在具体的实施过程中，可以利用电测井数据，确定出电阻率异常区域对应的深度值，以及区域内的电阻率异常平均值。

本说明书一个实施例中，可以沿待解释剖面横向追踪该电阻率异常区域，记录该电阻率异常区域对应的深度值和沉积岩薄层底界以上预设厚度以内的电阻率异常平均值。预设厚度的取值可以根据实际需要进行设置，本说明书实施例不作具体限定。

S5、根据所述深度值和所述电阻率平均值对所述电阻率异常区域进行搜索，将小于所述电阻率平均值的网格点的厚度进行叠加，获得所述电阻率异常区域对应的厚度数据。

在具体的实施过程中，可以从所述深度值对应的测试网格点对电阻率异常区域进行搜索，如：可以先纵向搜索，再横向搜索，或者可以采用其他搜索规则进行搜索。将小于电阻率平均值的网格点的厚度进行叠加，获得电阻率异常区域对应的厚度数据。

在本说明书一个实施例中，可以采用下述方法对电阻率异常区域进行搜索：

从所述深度值对应的测试网格点沿所述电阻率异常区域向下搜索，若搜索到的网格点的电阻率异常值小于所述电阻率异常平均值，则将所述网格点对应的厚度进行累加，直至搜索到的网格点的电阻率异常值等于所述电阻率异常平均值，确定出所述测试网格点对应的厚度数据；

采用相同的方法，对所述电阻率异常区域中其他网格点依次进行纵向搜索，如沿测试网格点的横向扫描，依次对各个网格点进行纵向搜索，获取所述电阻率异常区域中其他网格点对应的厚度数据，确定出待解释剖面对应的厚度数据。通过对电阻率异常区域中的网格点依次进行搜索解释，确定出待解释剖面的厚度数据，可以获得沉积岩薄层的厚度信息特征，为沉积岩薄层的地层解释提供了理论依据。

本说明书一个实施例中，还可以依次获取三维网格数据中各待解释剖面对应的厚度数据，将所述各待解释剖面对应的厚度数据进行剖面网格化，获得所述沉积岩薄层对应的地层厚度图。对待解释剖面完成上述操作后，可以得到待解释剖面这一层的厚度数据，它在平面上是一条线上的各个网格点的厚度值。不同的待解释剖面的网格点的厚度值，构成了平面的厚度数据分布，可以对这一平面分布的厚度数据(二维的平面的各网格点厚度值组成的数据体)进行网格化(也称规则化)，展示出厚度数据的平面分布成果(如，勾绘等值线图)。如：可以对所有待解释剖面的厚度数据结果进行平面网格化，再进行厚度平面等值线图输出，获得低阻沉积岩薄层的地层厚度图。通过对各个待解释剖面进行搜索解释，获得各剖面对应的厚度数据，并对获得的厚度数据进行平面网格化处理，获得的沉积岩薄层的平面分布，为沉积岩薄层的地层解释提供了理论依据。

本说明书实施例提供了一种适用于沉积岩薄层的电磁数据处理方法，通过对电磁数据进行三维网格化处理，根据获得的三维网格数据获得待解释剖面的异常剖面网格数据。对异常剖面网格数据进行标定，获得沉积岩薄层在剖面网格中的深度以及对应的电阻率异常区域。对获得的电阻率异常区域进行追踪搜索，获得待解释剖面的厚度数据，获得了低阻沉积岩薄层的厚度信息特征。实现了沉积岩薄层的识别和地层解释，提高了电磁数据对于高阻电性地层内夹的低阻沉积岩薄层目标的反映能力，为火山岩盆地及碳酸盐岩发育区的区带评价及目标优选提供了新方法，克服了约束反演依赖于初始模型的缺点。

图2是本申请又一个实施例中沉积岩薄层电磁数据的处理方法流程示意图，如图2所示，本说明书一个实施例中，在某盆地进行了三维电法勘探，应用上述实施例中的方法对三维电法反演后的三维电阻率数据实施了沉积岩薄层解释工作，实施例的具体过程和步骤如下：

1)对野外采集的某盆地三维大地电磁资料进行三维电阻率反演，地下空间剖分在低阻薄层发育深度段的深度剖分间距不大于薄层厚度(如：50m)，采用共轭梯度法反演，得到三维电阻率数据体，得到的三维电阻率数据体如图3所示，图3是本说明书一个实施例中三维电阻率数据体的立体示意图。

2)对三维电阻率数据体的电阻率数值取对数运算，并进行三维网格化，形成对数电阻率的三维网格数据。

3)从对数电阻率的三维网格数据中指定需要解释的剖面网格数据，对该剖面网格数据进行向上延拓500m计算，以向上延拓500m计算结果作为背景场而求取电阻率剩余异常的异常剖面网格数据。

4)用YC1电测井数据标定电阻率剩余异常的异常剖面网格数据，以YC1电测井数据的火山岩层段底界深度对应的异常剖面网格数据的深度为标志，确定火山岩岩层在异常剖面网格数据上的深度和对应的电阻率剩余异常带即电阻率异常区域，沿剖面横向追踪该电阻率剩余异常带，记录该电阻率剩余异常带对应的深度值和火山岩层段底界以上50m以内的电阻率异常平均值。

5)以记录的电阻率剩余异常带对应的深度值向下搜索，当搜索到的网格节点的电阻率剩余异常数值小于记录的电阻率异常平均值时就累加厚度，直至搜索电阻率剩余异常值等于记录的电阻率异常平均值的网格节点时，结束累加厚度计算。沿剖面横向推广进行搜索和记录各剖面上各测点下的最终累加厚度数值。

6)按步骤3)、4)、5)方法处理其余待解释剖面，然后对所有待解释剖面的厚度数据结果进行平面网格化，再进行厚度平面等值线图输出，等值线间距取5m，即为低阻沉积岩薄层的地层厚度图，地层厚度图的具体形式如图4所示，图4是本说明书一个实施例中获得的沉积岩薄层地层厚度图的示意图。

本说明书实施例通过对大地电磁三维电阻率数据进行数据处理和解释追踪，获得了低阻沉积岩薄层的平面分布和厚度信息特征，低阻沉积岩薄层深度约2500m-4000m，而可识别的低阻薄层的厚度在约50m-200m，该方法明显地提高电磁数据对于高阻电性地层内夹的低阻沉积岩薄层目标的反映能力，为火山岩盆地及碳酸盐岩发育区的区带评价及目标优选提供了新方法，克服了约束反演依赖于初始模型的缺点。

基于上述所述的沉积岩薄层电磁数据处理方法，本说明书一个或多个实施例还提供一种沉积岩薄层电磁数据处理装置。所述的装置可以包括使用了本说明书实施例所述方法的系统(包括分布式系统)、软件(应用)、模块、组件、服务器、客户端等并结合必要的实施硬件的装置。基于同一创新构思，本说明书实施例提供的一个或多个实施例中的装置如下面的实施例所述。由于装置解决问题的实现方案与方法相似，因此本说明书实施例具体的装置的实施可以参见前述方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。

具体地，图5是本说明书实施例提供的沉积岩薄层电磁数据处理装置一个实施例的模块结构示意图，如图5所示，本说明书实施例中提供的沉积岩薄层电磁数据处理装置包括：三维网格数据获取模块51，异常剖面网格获取模块52，异常区域确定模块53，电阻率平均值处理模块54，网格搜索模块55。

三维网格数据获取模块，可以用于获取采集到的电磁数据对应的三维网格数据；

异常剖面网格获取模块，可以用于从所述三维网格数据中获取待解释剖面的剖面网格数据，获取所述剖面网格数据对应的异常剖面网格数据；

异常区域确定模块，可以用于利用电测井数据对所述异常剖面网格数据进行标定，获得沉积岩薄层在所述异常剖面网格数据中对应的电阻率异常区域；

电阻率平均值处理模块，可以用于确定所述电阻率异常区域对应的深度值和电阻率异常平均值；

网格搜索模块，可以用于根据所述深度值和所述电阻率平均值对所述电阻率异常区域进行搜索，将小于所述电阻率平均值的网格点的厚度进行叠加，获得所述电阻率异常区域对应的厚度数据。

本说明书实施例提供的沉积岩薄层电磁数据处理装置，通过对电磁数据进行三维网格化处理，根据获得的三维网格数据获得待解释剖面的异常剖面网格数据。对异常剖面网格数据进行标定，获得沉积岩薄层在剖面网格中的深度以及对应的电阻率异常区域。对获得的电阻率异常区域进行追踪搜索，获得待解释剖面的厚度数据，获得了低阻沉积岩薄层的厚度信息特征。提高了电磁数据对于高阻电性地层内夹的低阻沉积岩薄层目标的反映能力，为火山岩盆地及碳酸盐岩发育区的区带评价及目标优选提供了新方法，克服了约束反演依赖于初始模型的缺点。

在上述实施例的基础上，本说明书一个实施例中，所述网格搜索模块具体用于：

从所述深度值对应的测试网格点沿所述电阻率异常区域向下搜索，若搜索到的网格点的电阻率异常值小于所述电阻率异常平均值，则将所述网格点对应的厚度进行累加，直至搜索到的网格点的电阻率异常值等于所述电阻率异常平均值，确定出所述测试网格点对应的厚度数据；

对所述电阻率异常区域中各网格点依次进行纵向搜索，获取所述电阻率异常区域中各网格点对应的厚度数据。

本说明书实施例，通过向上延拓计算，确定出待解释剖面的背景场数据和异常剖面网格数据，方法简单，为后续沉积岩薄层的识别以及厚度数据的获取奠定了理论基础。

在上述实施例的基础上，本说明书一个实施例中，所述装置还包括平面处理模块用于：

获取所述三维网格数据中各待解释剖面对应的厚度数据；

将所述各待解释剖面对应的厚度数据进行剖面网格化，获得所述沉积岩薄层对应的地层厚度图。

本说明书实施例，通过对各个待解释剖面进行搜索解释，获得各剖面对应的厚度数据，并对获得的厚度数据进行平面网格化处理，获得的沉积岩薄层的平面分布，为沉积岩薄层的地层解释提供了理论依据。

在上述实施例的基础上，本说明书一个实施例中，所述三维网格数据获取模块具体用于：

将获取到的电磁数据进行反演，获得三维电阻率数据体；

将所述三维电阻率数据体中的电阻率数值进行对数运算；

将对数运算后的三维电阻率数据体进行网格化，获得所述三维网格数据。

本说明书实施例，通过取对数运算可以降低采集到的电磁数据的数值变化率，再对取对数运算后的三维电阻率数据体进行网格化，使得网格数据比较平滑自然，方便后续对数据进行追踪解释。

在上述实施例的基础上，本说明书一个实施例中，所述异常剖面网格获取模块具体用于：

对所述剖面网格数据进行向上延拓计算，获得所述待解释剖面的背景场数据；

将所述剖面网格数据与所述背景场数据的差值作为所述异常剖面网格数据。

本说明书实施例，通过向上延拓计算，确定出待解释剖面的背景场数据和异常剖面网格数据，方法简单，为后续沉积岩薄层的识别以及厚度数据的获取奠定了理论基础。

在上述实施例的基础上，本说明书一个实施例中，所述电阻率平均值处理模块具体用于：

沿所述待解释剖面横向搜索所述电阻率异常区域，获取所述电阻率异常区域对应的深度值和所述沉积岩薄层底界向上预设厚度内的电阻率异常平均值。

本申请提供的沉积岩薄层电磁数据处理装置，通过对大地电磁三维电阻率数据进行数据处理和解释追踪，获得了低阻沉积岩薄层的平面分布和厚度信息特征，提高了电磁数据对于高阻电性地层内夹的低阻沉积岩薄层目标的反映能力，为火山岩盆地及碳酸盐岩发育区的区带评价及目标优选提供了新方法，克服了约束反演依赖于初始模型的缺点。

需要说明的，上述所述的装置根据方法实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。

上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。

本说明书提供的上述实施例所述的方法或装置可以通过计算机程序实现业务逻辑并记录在存储介质上，所述的存储介质可以计算机读取并执行，实现本说明书实施例所描述方案的效果。

本说明书实施例还可以提供一种计算机存储介质，其上存储有计算机程序，所述计算机程序被执行时，可以实现上述实施例中的沉积岩薄层电磁数据处理方法，例如：

获取采集到的电磁数据对应的三维网格数据；

从所述三维网格数据中获取待解释剖面的剖面网格数据，获取所述剖面网格数据对应的异常剖面网格数据；

利用电测井数据对所述异常剖面网格数据进行标定，获得沉积岩薄层在所述异常剖面网格数据中对应的电阻率异常区域；

确定所述电阻率异常区域对应的深度值和电阻率异常平均值；

根据所述深度值和所述电阻率平均值对所述电阻率异常区域进行搜索，将小于所述电阻率平均值的网格点的厚度进行叠加，获得所述电阻率异常区域对应的厚度数据。

本说明书实施例提供的上述沉积岩薄层电磁数据处理方法或装置可以在计算机中由处理器执行相应的程序指令来实现，如使用windows操作系统的c++语言在PC端实现、linux系统实现，或其他例如使用android、iOS系统程序设计语言在智能终端实现，以及基于量子计算机的处理逻辑实现等。本说明书提供的一种沉积岩薄层电磁数据处理系统的一个实施例中，图6是本说明书实施例提供的一种沉积岩薄层电磁数据处理系统实施例的模块结构示意图，如图6所示，本说明书另一实施例提供的沉积岩薄层电磁数据处理系统可以包括处理器61以及用于存储处理器可执行指令的存储器62，

处理器61和存储器62通过总线63完成相互间的通信；

所述处理器61用于调用所述存储器62中的程序指令，以执行上述各沉积岩薄层电磁数据处理方法实施例所提供的方法，例如包括：获取采集到的电磁数据对应的三维网格数据；从所述三维网格数据中获取待解释剖面的剖面网格数据，获取所述剖面网格数据对应的异常剖面网格数据；利用电测井数据对所述异常剖面网格数据进行标定，获得沉积岩薄层在所述异常剖面网格数据中对应的电阻率异常区域；确定所述电阻率异常区域对应的深度值和电阻率异常平均值；根据所述深度值和所述电阻率平均值对所述电阻率异常区域进行搜索，将小于所述电阻率平均值的网格点的厚度进行叠加，获得所述电阻率异常区域对应的厚度数据。

上述沉积岩薄层电磁数据处理系统还可以实现如下方法：利用可控震源配备的电台接收其他可控震源发布的状态询问信息，所述电台的通讯距离大于或等于可控震源的可同时激发距离；根据所述状态询问信息，向所述其他可控震源反馈状态信息，以使得所述其他可控震源根据所述状态信息进行扫描控制，所述状态信息包括：可控震源的位置信息、激发状态、将要激发的时间。

需要说明的是说明书上述所述的装置根据相关方法实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照方法实施例的描述，在此不作一一赘述。本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于硬件+程序类实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。

本说明书实施例并不局限于必须是符合行业通信标准、标准计算机数据处理和数据存储规则或本说明书一个或多个实施例所描述的情况。某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取、存储、判断、处理方式等获取的实施例，仍然可以属于本说明书实施例的可选实施方案范围之内。

在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)，如调制的数据信号和载波。

本领域技术人员应明白，本说明书一个或多个实施例可提供为方法、系统或计算机程序产品。因此，本说明书一个或多个实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本说明书一个或多个实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。

以上所述仅为本说明书一个或多个实施例的实施例而已，并不用于限制本本说明书一个或多个实施例。对于本领域技术人员来说，本说明书一个或多个实施例可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在权利要求范围之内。

Claims (12)

1.一种沉积岩薄层电磁数据处理方法，其特征在于，包括：

获取采集到的电磁数据对应的三维网格数据；

从所述三维网格数据中获取待解释剖面的剖面网格数据，获取所述剖面网格数据对应的异常剖面网格数据，所述异常剖面网格数据表示采集的电磁数据中有波动的数据；

利用电测井数据对所述异常剖面网格数据进行标定，获得沉积岩薄层在所述异常剖面网格数据中对应的电阻率异常区域；

确定所述电阻率异常区域对应的深度值和电阻率异常平均值；

根据所述深度值和所述电阻率平均值对所述电阻率异常区域进行搜索，将小于所述电阻率平均值的网格点的厚度进行叠加，获得所述电阻率异常区域对应的厚度数据；

其中，所述获取所述剖面网格数据对应的异常剖面网格数据，包括：

对所述剖面网格数据进行向上延拓计算，获得所述待解释剖面的背景场数据；

将所述剖面网格数据与所述背景场数据的差值作为所述异常剖面网格数据。

2.如权利要求1所述的一种沉积岩薄层电磁数据处理方法，其特征在于，所述根据所述深度值和所述电阻率异常平均值对所述电阻率异常区域进行搜索，将小于所述电阻率平均值的网格点的厚度进行叠加，获得所述电阻率异常区域对应的厚度数据，包括：

从所述深度值对应的测试网格点沿所述电阻率异常区域向下搜索，若搜索到的网格点的电阻率异常值小于所述电阻率异常平均值，则将所述网格点对应的厚度进行累加，直至搜索到的网格点的电阻率异常值等于所述电阻率异常平均值，确定出所述测试网格点对应的厚度数据；

对所述电阻率异常区域中各网格点依次进行纵向搜索，获取所述电阻率异常区域中各网格点对应的厚度数据。

3.如权利要求1所述的一种沉积岩薄层电磁数据处理方法，其特征在于，所述方法还包括：

获取所述三维网格数据中各待解释剖面对应的厚度数据；

将所述各待解释剖面对应的厚度数据进行剖面网格化，获得所述沉积岩薄层对应的地层厚度图。

4.如权利要求1所述的一种沉积岩薄层电磁数据处理方法，其特征在于，所述获取采集到的电磁数据对应的三维网格数据，包括：

将获取到的电磁数据进行反演，获得三维电阻率数据体；

将所述三维电阻率数据体中的电阻率数值进行对数运算；

将对数运算后的三维电阻率数据体进行网格化，获得所述三维网格数据。

5.如权利要求1所述的一种沉积岩薄层电磁数据处理方法，其特征在于，所述确定所述电阻率异常区域对应的深度值和电阻率异常平均值，包括：

沿所述待解释剖面横向搜索所述电阻率异常区域，获取所述电阻率异常区域对应的深度值和所述沉积岩薄层底界向上预设厚度内的电阻率异常平均值。

6.一种沉积岩薄层电磁数据处理装置，其特征在于，包括：

三维网格数据获取模块，用于获取采集到的电磁数据对应的三维网格数据；

异常剖面网格获取模块，用于从所述三维网格数据中获取待解释剖面的剖面网格数据，获取所述剖面网格数据对应的异常剖面网格数据，所述异常剖面网格数据表示采集的电磁数据中有波动的数据；

异常区域确定模块，用于利用电测井数据对所述异常剖面网格数据进行标定，获得沉积岩薄层在所述异常剖面网格数据中对应的电阻率异常区域；

电阻率平均值处理模块，用于确定所述电阻率异常区域对应的深度值和电阻率异常平均值；

网格搜索模块，用于根据所述深度值和所述电阻率平均值对所述电阻率异常区域进行搜索，将小于所述电阻率平均值的网格点的厚度进行叠加，获得所述电阻率异常区域对应的厚度数据；

其中，所述异常剖面网格获取模块具体用于：

对所述剖面网格数据进行向上延拓计算，获得所述待解释剖面的背景场数据；

将所述剖面网格数据与所述背景场数据的差值作为所述异常剖面网格数据。

7.如权利要求6所述的一种沉积岩薄层电磁数据处理装置，其特征在于，所述网格搜索模块具体用于：

从所述深度值对应的测试网格点沿所述电阻率异常区域向下搜索，若搜索到的网格点的电阻率异常值小于所述电阻率异常平均值，则将所述网格点对应的厚度进行累加，直至搜索到的网格点的电阻率异常值等于所述电阻率异常平均值，确定出所述测试网格点对应的厚度数据；

对所述电阻率异常区域中各网格点依次进行纵向搜索，获取所述电阻率异常区域中各网格点对应的厚度数据。

8.如权利要求6所述的一种沉积岩薄层电磁数据处理装置，其特征在于，所述装置还包括平面处理模块用于：

获取所述三维网格数据中各待解释剖面对应的厚度数据；

将所述各待解释剖面对应的厚度数据进行剖面网格化，获得所述沉积岩薄层对应的地层厚度图。

9.如权利要求6所述的一种沉积岩薄层电磁数据处理装置，其特征在于，所述三维网格数据获取模块具体用于：

将获取到的电磁数据进行反演，获得三维电阻率数据体；

将所述三维电阻率数据体中的电阻率数值进行对数运算；

将对数运算后的三维电阻率数据体进行网格化，获得所述三维网格数据。

10.如权利要求6所述的一种沉积岩薄层电磁数据处理装置，其特征在于，所述电阻率平均值处理模块具体用于：

沿所述待解释剖面横向搜索所述电阻率异常区域，获取所述电阻率异常区域对应的深度值和所述沉积岩薄层底界向上预设厚度内的电阻率异常平均值。

11.一种计算机存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被执行时，实现权利要求1-5任一项所述的方法。

12.一种沉积岩薄层电磁数据处理系统，其特征在于，包括处理器以及用于存储处理器可执行指令的存储器，所述处理器执行所述指令时实现如权利要求1至5中任意一项所述方法的步骤。

Images