CN104122590B - 一种基于电磁勘探的油气检测方法及系统 - Google Patents

Abstract

本发明公开了一种基于电磁勘探的油气检测方法及系统，该方法包括：获取实验测线的界面分层信息及深测井电阻率数据，计算获得加权平均值电阻率；建立分层电阻率二维模型，设定实验测试频率，并计算电磁响应，获取分层电阻率二维模型的原始振幅数据以及原始相位数据；进行去噪和滤波后，获得去噪后的振幅数据和去噪后的相位数据；进行差分处理，获得差分振幅数据或差分相位数据；选取并计算差分振幅数据或差分相位数据中高频段部分对应的平均值，进行归零处理获得岩性异常振幅数据或岩性异常相位数据；计算获得振幅剩余电磁数据或相位剩余电磁数据，并以此绘制实验测线的拟深度域异常等值线断面图，根据拟深度域异常等值线断面图生成油气检测结果。

Description

一种基于电磁勘探的油气检测方法及系统

技术领域

本发明涉及地球物理油气勘探领域，尤指一种基于电磁勘探的油气检测方法及系统。

背景技术

在现今的油气勘探领域中，特别是隐蔽油气藏的勘探，油气直接检测方法一直被油气勘探家所推崇和重视。地震勘探中就有很多识别含油气异常的地震属性，众所周知的“亮点”技术一直应用于地震勘探用于识别油气目标，但是地震属性研究也不能十拿九稳，钻探成功率依然没有因此而得到明显提高，主要还在于油气与水的波阻抗以及充填油气与没有充填波阻抗差别不明显。地球化学方法被喻为唯一的直接找油方法，其通过直接测量烃类成分含量，但由于只是蜻蜓点水，在地表采集，与地下深处的油气藏距离太远，效果往往不如人意；其它方法诸如遥感、微磁等也都类似。

由于岩石电性与油气关系密切，油气与水的电阻率差异可达上1000倍，充填油气与没有充填电阻率差别非常明显，因此，电法在油气检测中占用重要地位，电测井中的电阻率、自然电位和极化率是识别评价含油气储层的主要参数，所以地面电法的油气检测技术也一直备受关注，电磁油气检测方法一直被广泛重视，方法不断推陈出新，主要方法有双频振幅和双频相位，但由于受背景异常影响，不能有效识别目标异常，而极化率反演方法由于不能同时反演电阻率和极化率，同时受反演方法和拟合精度影响，效果不能保障，难以有效提取反映油气的信息。

发明内容

针对现有技术中油气检测时会受背景异常影响，不能有效识别目标异常，使得效果不能保障，难以有效提取反映油气信息的问题，本发明对电磁油气检测方法进行了改进，通过剔除构造背景，计算并绘制振幅剩余电磁和相位剩余电磁的拟深度域异常等值线断面图，生成油气检测结果，推断油气有利区，实现圈闭的油气检测。

为达到上述目的，本发明提出了一种基于电磁勘探的油气检测方法，包括以下步骤：步骤1，获取实验测线相应的地震深度偏移剖面的界面分层信息，以及所述实验测线的深测井电阻率数据，并计算获得加权平均值电阻率；步骤2，根据所述界面分层信息以及所述深测井电阻率数据建立分层电阻率二维模型，设定实验测试频率，并计算所述分层电阻率二维模型的电磁响应，获取所述分层电阻率二维模型的原始振幅数据以及原始相位数据；步骤3，对所述原始振幅数据及原始相位数据进行去噪和滤波，获得去噪后的振幅数据和去噪后的相位数据；步骤4，对所述原始振幅数据、去噪后的振幅数据在频率对数域进行差分处理，获得差分振幅数据，并对所述原始相位数据、去噪后的相位数据在线性域进行差分处理后，获得差分相位数据；步骤5，选取并计算所述差分振幅数据或所述差分相位数据中高频段部分对应的平均值，并根据所述差分振幅数据或所述差分相位数据，进行归零处理获得岩性异常振幅数据或岩性异常相位数据；步骤6，根据所述加权平均值电阻率、所述岩性异常振幅数据或岩性异常相位数据计算获得振幅剩余电磁数据或相位剩余电磁数据，并根据所述振幅剩余电磁数据及相位剩余电磁数据绘制所述实验测线的拟深度域异常等值线断面图，根据拟深度域异常等值线断面图生成油气检测结果。

为达到上述目的，本发明还提出了一种基于电磁勘探的油气检测系统，包括：数据采集模块，用于获取实验测线相应的地震深度偏移剖面的界面分层信息，以及所述实验测线的深测井电阻率数据，并计算获得加权平均值电阻率；原始数据获取模块，用于根据所述界面分层信息以及所述深测井电阻率数据建立分层电阻率二维模型，设定实验测试频率，并计算所述分层电阻率二维模型的电磁响应，获取所述分层电阻率二维模型的原始振幅数据以及原始相位数据；去噪滤波模块，用于对所述原始振幅数据及原始相位数据进行去噪和滤波，获得去噪后的振幅数据和去噪后的相位数据；差分处理模块，用于对所述原始振幅数据、去噪后的振幅数据在频率对数域进行差分处理，获得差分振幅数据，并对所述原始相位数据、去噪后的相位数据在线性域进行差分处理后，获得差分相位数据；岩性异常数据获取模块，用于选取并计算所述差分振幅数据或所述差分相位数据中高频段部分对应的平均值，并根据所述差分振幅数据或所述差分相位数据，进行归零处理获得岩性异常振幅数据或岩性异常相位数据；油气检测模块，用于根据所述加权平均值电阻率、所述岩性异常振幅数据或岩性异常相位数据计算获得振幅剩余电磁数据或相位剩余电磁数据，并根据所述振幅剩余电磁数据及相位剩余电磁数据绘制所述实验测线的拟深度域异常等值线断面图，根据拟深度域异常等值线断面图生成油气检测结果。

本发明的基于电磁勘探的油气检测方法及系统可以有效提取反映油气的信息，通过在探区地震构造信息和探区电测井综合建模基础上，利用二维正演模拟计算出由精确构造产生的电磁异常，通过去除构造信息和近地表岩性影响获得剩余异常，主要反映地层岩性变化的信息，再采取差分处理进一步提取油气相关异常信息，最后绘制测线二维剩余差分的振幅异常和相位异常，根据异常分布圈定油气目标区，准确的实现圈闭的油气检测。

附图说明

此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中：

图1为本发明一实施例的基于电磁勘探的油气检测方法流程图。

图2为本发明一实施例的基于电磁勘探的油气检测系统的结构示意图。

图3为本发明一具体实施例的深度域振幅油气异常等值线断面图。

图4为本发明一具体实施例的深度域相位油气异常等值线断面图。

具体实施方式

以下配合图式及本发明的较佳实施例，进一步阐述本发明为达成预定发明目的所采取的技术手段。

图1为本发明一实施例的基于电磁勘探的油气检测方法流程图。如图1所示，该方法包括：

步骤S101，首先选择要处理的测线，获取该实验测线的地震资料和相应电测井资料。其中，地震资料主要是地震深度偏移剖面的界面分层信息，即界面深度和每层厚度H_w等数据。实验测线的电测井资料主要是深测井电阻率数据，并按地震分层统计出每一层的电阻率平均值ρ_w。

在本步骤中，还根据上述界面分层信息中各个分层的厚度H_w以及深测井电阻率中各个分层的电阻率平均值ρ_w计算获得加权平均值电阻率ρ₀，利用的公式为：

${\mathrm{\ρ}}_0=\left(\mathrm{\Σ}\frac{{\mathrm{\ρ}}_w}{H_w}\right)\×\left(\mathrm{\Σ}{H}_w\right);---\left(1\right)$

其中，ρ₀为加权平均值电阻率；

H_w为界面分层信息中各个分层的厚度；

ρ_w为深测井电阻率中各个分层的电阻率平均值；

角标w为各个分层的层序号。

步骤S102，根据步骤S101获得的界面分层信息以及深测井电阻率数据建立分层电阻率二维模型，设定实验测试频率f_i，并利用电磁法二维正演软件计算分层电阻率二维模型的电磁响应，获取分层电阻率二维模型的原始振幅数据E^A _i以及原始相位数据F^A _i。

在本步骤中，设定的实验测试频率f_i可以包括M个频率，M值大小和步长与实际测试频率和步长一致；角标i为实验测试频率f_i的序号，取值为1至M，与其相应的每一原始振幅数据E^A _i以及原始相位数据F^A _i等数据也分别有M个。

步骤S103，对原始振幅数据E^A _i及原始相位数据F^A _i进行去噪和滤波，获得去噪后的振幅数据E^R _i和去噪后的相位数据F^R _i。

步骤S104，对原始振幅数据E^A _i、去噪后的振幅数据E^R _i在频率对数域进行差分处理，获得差分振幅数据E^d _i，并对原始相位数据F^A _i、去噪后的相位数据F^R _i在线性域进行差分处理后，获得差分相位数据F^d _i。

在本步骤中，利用的公式有：

E^d _i＝lg(E^A _i-E^R _i)； (2)

F^d _i＝lg(F^A _i-F^R _i)； (3)

其中，E^d _i为差分振幅数据；

E^A _i为步骤S102获取的原始振幅数据；

E^R _i为步骤S103获取的去噪后的振幅数据；

F^d _i为差分相位数据；

F^A _i为步骤S102获取的原始相位数据；

F^R _i为步骤S103获取的去噪后的相位数据。

步骤S105，选取并计算差分振幅数据E^d _i或差分相位数据F^d _i中高频段部分对应的平均值，并根据差分振幅数据E^d _i或差分相位数据F^d _i，进行归零处理获得岩性异常振幅数据E^c _i或岩性异常相位数据F^c _i。

在本步骤中，选取并计算差分振幅数据E^d _i或差分相位数据F^d _i中高频段部分对应的平均值利用的公式有：

$E^h=\frac{\mathrm{\Σ}{E^d}_g}{n};---\left(4\right)$

$F^h=\frac{\mathrm{\Σ}{F^d}_g}{n};---\left(5\right)$

其中，E^h为差分振幅数据E^d _i中高频段部分对应的平均值；

E^d _g为差分振幅数据中高频段部分对应的各个差分振幅数据；

F^h为差分相位数据中高频段部分对应的平均值；

F^d _g为差分相位数据F^d _i中高频段部分对应的各个差分相位数据；

n为高频段部分对应的差分振幅数据或差分相位数据的个数；

角标g为高频段部分对应的各个差分振幅数据或差分相位数据的序号，取值为1至n，也就是从高频段部分对应的各个数据是从序号1开始的前n个。n与M还有以下关系： $n<\frac{M}{3}.$

在本步骤中，进行归零处理获得岩性异常振幅数据或岩性异常相位数据利用的公式为：

E^c _i＝E^d _i-E^h； (6)

F^c _i＝F^d _i-F^h； (7)

其中，E^c _i为岩性异常振幅数据；

E^d _i为公式(2)得到的差分振幅数据；

E^h为公式(4)差分振幅数据中高频段部分对应的平均值；

F^c _i为岩性异常相位数据；

F^d _i为公式(3)得到的差分相位数据；

F^h为公式(5)差分相位数据中高频段部分对应的平均值。

步骤S106，根据加权平均值电阻率ρ₀、岩性异常振幅数据E^c _i或岩性异常相位数据F^c _i计算获得振幅剩余电磁数据DualE_i或相位剩余电磁数据DualF_i，并根据振幅剩余电磁数据DualE_i及相位剩余电磁数据DualF_i绘制实验测线的拟深度域异常等值线断面图，根据拟深度域异常等值线断面图生成油气检测结果。

在本步骤中，计算获得振幅剩余电磁数据DualE_i或相位剩余电磁数据DualF_i利用的公式为：

${\mathrm{DualE}}_i=\frac{({E^c}_i-{E^{c3}}_i)}{({E^c}_i+{E^{c3}}_i)}\&times;356\&times;\sqrt{\left(\frac{{\mathrm{\&rho;}}_0}{f_i}\right)};---\left(8\right)$

${\mathrm{DualF}}_i=\frac{({F^c}_i-{F^{c3}}_i)}{({F^c}_i+{F^{c3}}_i)}\&times;356\&times;\sqrt{\left(\frac{{\mathrm{\&rho;}}_0}{f_i}\right)};---\left(9\right)$

其中，DualE_i为振幅剩余电磁数据；

E^c _i为公式(6)得到的岩性异常振幅数据；

E^c3 _i为每一岩性异常振幅数据E^c _i对应的3次谐波的岩性异常振幅数据；

ρ₀为公式(1)得到的加权平均值电阻率；

f_i为步骤2设定的实验测试频率；

DualF_i为相位剩余电磁数据；

F^c _i为公式(7)得到的岩性异常相位数据；

F^c3 _i为每一岩性异常相位数据F^c _i对应的3次谐波的岩性异常相位数据。

在本步骤中，3次谐波(频率是基波3倍)的数据是直接采用分频技术得到的，属于信号分析中是常用的技术手段。

在本步骤中，最后生成油气检测结果主要包括：利用得到的振幅剩余电磁数据DualE_i和相位剩余电磁数据DualF_i，绘制实验测线的拟深度域异常等值线断面图。如果有已知探井数据，根据探井数据的油气情况确定背景值；如果没有已知探井数据，选择剩余电磁的平均值作为背景值；最后，根据实验测线的拟深度域异常等值线断面图，当振幅剩余电磁数据或相位剩余电磁数据高于背景值则为异常，两者都为异常的区段为油气有利区。

基于同一发明构思，本发明实施例中还提供了一种基于电磁勘探的油气检测系统，如下面的实施例所述。

由于基于电磁勘探的油气检测系统解决问题的原理与基于电磁勘探的油气检测方法相似，因此基于电磁勘探的油气检测系统的实施可以参见基于电磁勘探的油气检测方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的系统较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。

图2为本发明一实施例的基于电磁勘探的油气检测系统的结构示意图。如图2所示，该系统包括：

数据采集模块11，用于获取实验测线相应的地震深度偏移剖面的界面分层信息，以及实验测线的深测井电阻率数据，并根据界面分层信息及深测井电阻率数据计算获得加权平均值电阻率ρ_w；

原始数据获取模块12，用于根据界面分层信息以及深测井电阻率数据建立分层电阻率二维模型，设定实验测试频率f_i，并利用电磁法二维正演软件计算分层电阻率二维模型的电磁响应，获取分层电阻率二维模型的原始振幅数据E^A _i以及原始相位数据F^A _i；

去噪滤波模块13，用于对原始振幅数据E^A _i及原始相位数据F^A _i进行去噪和滤波，获得去噪后的振幅数据E^R _i和去噪后的相位数据F^R _i；

差分处理模块14，用于对原始振幅数据E^A _i、去噪后的振幅数据E^R _i在频率对数域进行差分处理，获得差分振幅数据E^d _i，并对原始相位数据F^A _i、去噪后的相位数据F^R _i在线性域进行差分处理后，获得差分相位数据F^d _i；

岩性异常数据获取模块15，用于选取并计算差分振幅数据E^d _i或差分相位数据F^d _i中高频段部分对应的平均值，并根据差分振幅数据E^d _i或差分相位数据F^d _i，进行归零处理获得岩性异常振幅数据E^c _i或岩性异常相位数据F^c _i；

油气检测模块16，用于根据加权平均值电阻率ρ₀、岩性异常振幅数据E^c _i或岩性异常相位数据F^c _i计算获得振幅剩余电磁数据DualE_i或相位剩余电磁数据DualF_i，并根据振幅剩余电磁数据DualE_i及相位剩余电磁数据DualF_i绘制实验测线的拟深度域异常等值线断面图，根据拟深度域异常等值线断面图生成油气检测结果。

在本实施例中，数据采集模块11，计算获得加权平均值电阻率ρ₀利用的公式为：

${\mathrm{\&rho;}}_0=\left(\mathrm{\&Sigma;}\frac{{\mathrm{\&rho;}}_w}{H_w}\right)\&times;\left(\mathrm{\&Sigma;}{H}_w\right);---\left(1\right)$

其中，ρ₀为加权平均值电阻率；

H_w为界面分层信息中各个分层的厚度；

ρ_w为深测井电阻率中各个分层的电阻率平均值。

差分处理模块14，获得差分振幅数据E^d _i以及获得差分相位数据F^d _i利用的公式为：

E^d _i＝lg(E^A _i-E^R _i)； (2)

F^d _i＝lg(F^A _i-F^R _i)； (3)

其中，E^d _i为差分振幅数据；

E^A _i为原始数据获取模块12获取的原始振幅数据；

E^R _i为去噪滤波模块13获取的去噪后的振幅数据；

F^d _i为差分相位数据；

F^A _i为原始数据获取模块12获取的原始相位数据；

F^R _i为去噪滤波模块13获取的去噪后的相位数据。

在本实施例中，岩性异常数据获取模块15，选取并计算差分振幅数据E^d _i或差分相位数据F^d _i中高频段部分对应的平均值利用的公式为：

$E^h=\frac{\mathrm{\&Sigma;}{E^d}_g}{n};---\left(4\right)$

$F^h=\frac{\mathrm{\&Sigma;}{F^d}_g}{n};---\left(5\right)$

其中，E^h为差分振幅数据中高频段部分对应的平均值；

E^d _g为差分振幅数据E^d _i中高频段部分对应的各个差分振幅数据；

F^h为差分相位数据中高频段部分对应的平均值；

F^d _g为差分相位数据F^d _i中高频段部分对应的各个差分相位数据；

n为高频段部分对应的差分振幅数据或差分相位数据的个数。

在本实施例中，岩性异常数据获取模块15，根据差分振幅数据E^d _i或差分相位数据F^d _i，进行归零处理获得岩性异常振幅数据E^c _i或岩性异常相位数据F^c _i利用的公式为；

E^c _i＝E^d _i-E^h； (6)

F^c _i＝F^d _i-F^h； (7)

其中，E^c _i为岩性异常振幅数据；

E^d _i为公式(2)得到的差分振幅数据；

E^h为公式(4)差分振幅数据中高频段部分对应的平均值；

F^c _i为岩性异常相位数据；

F^d _i为公式(3)得到的差分相位数据；

F^h为公式(5)差分相位数据中高频段部分对应的平均值。

在本实施例中，油气检测模块16，根据加权平均值电阻率、岩性异常振幅数据或岩性异常相位数据计算获得振幅剩余电磁数据或相位剩余电磁数据利用的公式为；

${\mathrm{DualE}}_i=\frac{({E^c}_i-{E^{c3}}_i)}{({E^c}_i+{E^{c3}}_i)}\&times;356\&times;\sqrt{\left(\frac{{\mathrm{\&rho;}}_0}{f_i}\right)};---\left(8\right)$

${\mathrm{DualF}}_i=\frac{({F^c}_i-{F^{c3}}_i)}{({F^c}_i+{F^{c3}}_i)}\&times;356\&times;\sqrt{\left(\frac{{\mathrm{\&rho;}}_0}{f_i}\right)};---\left(9\right)$

其中，DualE_i为振幅剩余电磁数据；

E^c _i为公式(6)得到的岩性异常振幅数据；

E^c3 _i为每一岩性异常振幅数据E^c _i对应的3次谐波的岩性异常振幅数据；

ρ₀为公式(1)得到的加权平均值电阻率；

f_i为原始数据获取模块12中设定的实验测试频率；

DualF_i为相位剩余电磁数据；

F^c _i为公式(7)得到的岩性异常相位数据；

F^c3 _i为每一岩性异常相位数据F^c _i对应的3次谐波的岩性异常相位数据。

在本实施例中，油气检测模块16，根据拟深度域异常等值线断面图生成油气检测结果包括：

如果有已知探井数据，根据探井数据的油气情况确定背景值；如果没有已知探井数据，选择剩余电磁的平均值作为背景值；根据实验测线的拟深度域异常等值线断面图，当振幅剩余电磁数据或相位剩余电磁数据高于背景值则为异常，两者都为异常的区段为油气有利区。

为了对上述弹性波正演模拟方法进行更为清楚的解释，下面结合一个具体的实施例来进行说明，然而值得注意的是该实施例仅是为了更好地说明本发明，并不构成对本发明不当的限定。

结合步骤S101，选择要处理的测线，收集该测线的地震深度偏移剖面的界面分层信息，即界面深度和每层厚度H_w，一共20层，收集测井电阻率数据，并按地震分层统计出每一层的电阻率平均值ρ_w，同时，计算出加权平均值电阻率ρ₀。

结合步骤S102，根据前一步获得的界面分层信息以及深测井电阻率数据建立分层电阻率二维模型，并根据实验需求设定40个实验测试频率(步长与实测频率和步长一致)，利用电磁法二维正演软件计算分层电阻率二维模型的电磁响应，获取分层电阻率二维模型的原始振幅数据E^A _i以及原始相位数据F^A _i，这些数据对应有40个。

结合步骤S103、步骤S104，对前一部数据进行去噪滤波处理后，分别计算出差分振幅数据E^d _i以及差分相位数据F^d _i，这样就去除了剖面的构造背景，获得剖面的岩性异常。

结合步骤S105，由于浅地表岩性变化常常影响整条曲线，因此需要消除。为了去除近地表影响，将高频段数据做归零处理：

首先，求取高频段5个频率数据的平均值，利用公式(4)、(5)，其中，n取5，g为1至5，计算获得差分振幅数据E^d _i或差分相位数据F^d _i中高频段部分对应的平均值。

然后利用公式(6)、(7)，进行归零处理，获得纯反应深部岩性异常的岩性异常振幅数据E^c _i或岩性异常相位数据F^c _i。

计算出上述数据后，再结合步骤S106，利用公式(8)、(9)计算获得振幅剩余电磁数据DualE_i及相位剩余电磁数据DualF_i。

最后，根据上述数据，绘制该二维测线的振幅剩余电磁(DualE_i)和相位剩余电磁(DualF_i)的拟深度域异常等值线断面图，图3为本实施例的深度域振幅油气异常等值线断面图，图4为本实施例的深度域相位油气异常等值线断面图。结合图3及图4所示，如果有已知探井，则根据探井的油气情况确定背景值，没有已知探井则选择平均值作为背景值，高于背景值得则为异常，两者都为高异常的区段即可为油气有利区。

本发明的基于电磁勘探的油气检测方法及系统可以有效提取反映油气的信息，通过在探区地震构造信息和探区电测井综合建模基础上，利用二维正演模拟计算出由精确构造产生的电磁异常，通过去除构造信息和近地表岩性影响获得剩余异常，主要反映地层岩性变化的信息，再采取差分处理进一步提取油气相关异常信息，最后绘制测线二维剩余差分的振幅异常(图3)和相位异常(图4)，根据异常分布圈定油气目标区，准确的实现圈闭的油气检测。

以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (12)

1.一种基于电磁勘探的油气检测方法，其特征在于，包括以下步骤：

步骤1，获取实验测线相应的地震深度偏移剖面的界面分层信息，以及所述实验测线的深测井电阻率数据，并计算获得加权平均值电阻率；其中，计算获得加权平均值电阻率利用的公式为：

${\mathrm{\&rho;}}_0=\left(\mathrm{\&Sigma;}\frac{{\mathrm{\&rho;}}_w}{H_w}\right)\&times;\left({\mathrm{\&Sigma;H}}_w\right);---\left(1\right)$

其中，ρ₀为所述加权平均值电阻率；

H_w为所述界面分层信息中各个分层的厚度；

ρ_w为所述深测井电阻率中各个分层的电阻率平均值；

步骤2，根据所述界面分层信息以及所述深测井电阻率数据建立分层电阻率二维模型，设定实验测试频率，并计算所述分层电阻率二维模型的电磁响应，获取所述分层电阻率二维模型的原始振幅数据以及原始相位数据；

步骤3，对所述原始振幅数据及原始相位数据进行去噪和滤波，获得去噪后的振幅数据和去噪后的相位数据；

步骤4，对所述原始振幅数据、去噪后的振幅数据在频率对数域进行差分处理，获得差分振幅数据，并对所述原始相位数据、去噪后的相位数据在线性域进行差分处理后，获得差分相位数据；

步骤5，选取并计算所述差分振幅数据或所述差分相位数据中高频段部分对应的平均值，并根据所述差分振幅数据或所述差分相位数据，进行归零处理获得岩性异常振幅数据或岩性异常相位数据；

步骤6，根据所述加权平均值电阻率、所述岩性异常振幅数据或岩性异常相位数据计算获得振幅剩余电磁数据或相位剩余电磁数据，并根据所述振幅剩余电磁数据及相位剩余电磁数据绘制所述实验测线的拟深度域异常等值线断面图，根据拟深度域异常等值线断面图生成油气检测结果。

2.根据权利要求1所述的基于电磁勘探的油气检测方法，其特征在于，所述步骤4中获得差分振幅数据以及获得差分相位数据利用的公式为：

E^d _i＝lg(E^A _i-E^R _i) ； (2)

F^d _i＝lg(F^A _i-F^R _i)； (3)

其中，E^d _i为所述差分振幅数据；

E^A _i为所述步骤2获取的原始振幅数据；

E^R _i为所述步骤3获取的去噪后的振幅数据；

F^d _i为所述差分相位数据；

F^A _i为所述步骤2获取的原始相位数据；

F^R _i为所述步骤3获取的去噪后的相位数据。

3.根据权利要求2所述的基于电磁勘探的油气检测方法，其特征在于，所述步骤5中，选取并计算所述差分振幅数据或所述差分相位数据中高频段部分对应的平均值利用的公式为：

$E^h=\frac{{{\mathrm{\&Sigma;E}}^d}_g}{n};---\left(4\right)$

$F^h=\frac{{{\mathrm{\&Sigma;F}}^d}_g}{n};---\left(5\right)$

其中，E^h为所述差分振幅数据中高频段部分对应的平均值；

E^d _g为所述差分振幅数据中高频段部分对应的各个差分振幅数据；

F^h为所述差分相位数据中高频段部分对应的平均值；

F^d _g为所述差分相位数据中高频段部分对应的各个差分相位数据；

n为所述高频段部分对应的差分振幅数据或差分相位数据的个数。

4.根据权利要求3所述的基于电磁勘探的油气检测方法，其特征在于，所述步骤5中，根据所述差分振幅数据或所述差分相位数据，进行归零处理获得岩性异常振幅数据或岩性异常相位数据利用的公式为；

E^c _i＝E^d _i-E^h； (6)

F^c _i＝F^d _i-F^h； (7)

其中，E^c _i为所述岩性异常振幅数据；

E^d _i为公式(2)得到的所述差分振幅数据；

E^h为公式(4)所述差分振幅数据中高频段部分对应的平均值；

F^c _i为所述岩性异常相位数据；

F^d _i为公式(3)得到的所述差分相位数据；

F^h为公式(5)所述差分相位数据中高频段部分对应的平均值。

5.根据权利要求4所述的基于电磁勘探的油气检测方法，其特征在于，所述步骤6中，根据所述加权平均值电阻率、所述岩性异常振幅数据或岩性异常相位数据计算获得振幅剩余电磁数据或相位剩余电磁数据利用的公式为；

${\mathrm{DualE}}_i=\frac{({E^c}_i-{E^{c3}}_i)}{({E^c}_i+{E^{c3}}_i)}\&times;356\&times;\sqrt{\left(\frac{{\mathrm{\&rho;}}_0}{f_i}\right)};---\left(8\right)$

${\mathrm{DualF}}_i=\frac{({F^c}_i-{F^{c3}}_i)}{({F^c}_i+{F^{c3}}_i)}\&times;356\&times;\sqrt{\left(\frac{{\mathrm{\&rho;}}_0}{f_i}\right)};---\left(9\right)$

其中，DualE_i为所述振幅剩余电磁数据；

E^c _i为公式(6)得到的所述岩性异常振幅数据；

E^c3 _i为每一所述岩性异常振幅数据E^c _i对应的3次谐波的岩性异常振幅数据；

ρ₀为公式(1)得到的所述加权平均值电阻率；

f_i为所述步骤2设定的所述实验测试频率；

DualF_i为所述相位剩余电磁数据；

F^c _i为公式(7)得到的所述岩性异常相位数据；

F^c3 _i为每一所述岩性异常相位数据F^c _i对应的3次谐波的岩性异常相位数据。

6.根据权利要求5所述的基于电磁勘探的油气检测方法，其特征在于，所述步骤6中，根据所述振幅剩余电磁数据及相位剩余电磁数据绘制所述实验测线的拟深度域异常等值线断面图，根据拟深度域异常等值线断面图生成油气检测结果包括：

如果有已知探井数据，根据所述探井数据的油气情况确定背景值；

如果没有已知探井数据，选择剩余电磁的平均值作为背景值；

根据所述实验测线的拟深度域异常等值线断面图，当所述振幅剩余电磁数据或相位剩余电磁数据高于背景值则为异常，两者都为异常的区段为油气有利区。

7.一种基于电磁勘探的油气检测系统，其特征在于，包括：

数据采集模块，用于获取实验测线相应的地震深度偏移剖面的界面分层信息，以及所述实验测线的深测井电阻率数据，并计算获得加权平均值电阻率；其中，计算获得加权平均值电阻率利用的公式为：

${\mathrm{\&rho;}}_0=\left(\mathrm{\&Sigma;}\frac{{\mathrm{\&rho;}}_w}{H_w}\right)\&times;\left({\mathrm{\&Sigma;H}}_w\right);---\left(1\right)$

其中，ρ₀为所述加权平均值电阻率；

H_w为所述界面分层信息中各个分层的厚度；

ρ_w为所述深测井电阻率中各个分层的电阻率平均值；

原始数据获取模块，用于根据所述界面分层信息以及所述深测井电阻率数据建立分层电阻率二维模型，设定实验测试频率，并计算所述分层电阻率二维模型的电磁响应，获取所述分层电阻率二维模型的原始振幅数据以及原始相位数据；

去噪滤波模块，用于对所述原始振幅数据及原始相位数据进行去噪和滤波，获得去噪后的振幅数据和去噪后的相位数据；

差分处理模块，用于对所述原始振幅数据、去噪后的振幅数据在频率对数域进行差分处理，获得差分振幅数据，并对所述原始相位数据、去噪后的相位数据在线性域进行差分处理后，获得差分相位数据；

岩性异常数据获取模块，用于选取并计算所述差分振幅数据或所述差分相位数据中高频段部分对应的平均值，并根据所述差分振幅数据或所述差分相位数据，进行归零处理获得岩性异常振幅数据或岩性异常相位数据；

油气检测模块，用于根据所述加权平均值电阻率、所述岩性异常振幅数据或岩性异常相位数据计算获得振幅剩余电磁数据或相位剩余电磁数据，并根据所述振幅剩余电磁数据及相位剩余电磁数据绘制所述实验测线的拟深度域异常等值线断面图，根据拟深度域异常等值线断面图生成油气检测结果。

8.根据权利要求7所述的基于电磁勘探的油气检测系统，其特征在于，所述差分处理模块，获得差分振幅数据以及获得差分相位数据利用的公式为：

E^d _i＝lg(E^A _i-E^R _i) ； (2)

F^d _i＝lg(F^A _i-F^R _i)； (3)

其中，E^d _i为所述差分振幅数据；

E^A _i为所述原始数据获取模块获取的原始振幅数据；

E^R _i为所述去噪滤波模块获取的去噪后的振幅数据；

F^d _i为所述差分相位数据；

F^A _i为所述原始数据获取模块获取的原始相位数据；

F^R _i为所述去噪滤波模块获取的去噪后的相位数据。

9.根据权利要求8所述的基于电磁勘探的油气检测系统，其特征在于，所述岩性异常数据获取模块，选取并计算所述差分振幅数据或所述差分相位数据中高频段部分对应的平均值利用的公式为：

$E^h=\frac{{{\mathrm{\&Sigma;E}}^d}_g}{n};---\left(4\right)$

$F^h=\frac{{{\mathrm{\&Sigma;F}}^d}_g}{n};---\left(5\right)$

其中，E^h为所述差分振幅数据中高频段部分对应的平均值；

E^d _g为所述差分振幅数据中高频段部分对应的各个差分振幅数据；

F^h为所述差分相位数据中高频段部分对应的平均值；

F^d _g为所述差分相位数据中高频段部分对应的各个差分相位数据；

n为所述高频段部分对应的差分振幅数据或差分相位数据的个数。

10.根据权利要求9所述的基于电磁勘探的油气检测系统，其特征在于，所述岩性异常数据获取模块，根据所述差分振幅数据或所述差分相位数据，进行归零处理获得岩性异常振幅数据或岩性异常相位数据利用的公式为；

E^c _i＝E^d _i-E^h； (6)

F^c _i＝F^d _i-F^h； (7)

其中，E^c _i为所述岩性异常振幅数据；

E^d _i为公式(2)得到的所述差分振幅数据；

E^h为公式(4)所述差分振幅数据中高频段部分对应的平均值；

F^c _i为所述岩性异常相位数据；

F^d _i为公式(3)得到的所述差分相位数据；

F^h为公式(5)所述差分相位数据中高频段部分对应的平均值。

11.根据权利要求10所述的基于电磁勘探的油气检测系统，其特征在于，所述油气检测模块，根据所述加权平均值电阻率、所述岩性异常振幅数据或岩性异常相位数据计算获得振幅剩余电磁数据或相位剩余电磁数据利用的公式为；

${\mathrm{DualE}}_i=\frac{({E^c}_i-{E^{c3}}_i)}{({E^c}_i+{E^{c3}}_i)}\&times;356\&times;\sqrt{\left(\frac{{\mathrm{\&rho;}}_0}{f_i}\right)};---\left(8\right)$

${\mathrm{DualF}}_i=\frac{({F^c}_i-{F^{c3}}_i)}{({F^c}_i+{F^{c3}}_i)}\&times;356\&times;\sqrt{\left(\frac{{\mathrm{\&rho;}}_0}{f_i}\right)};---\left(9\right)$

其中，DualE_i为所述振幅剩余电磁数据；

E^c _i为公式(6)得到的所述岩性异常振幅数据；

E^c3 _i为每一所述岩性异常振幅数据E^c _i对应的3次谐波的岩性异常振幅数据；

ρ₀为公式(1)得到的所述加权平均值电阻率；

f_i为所述原始数据获取模块中设定的所述实验测试频率；

DualF_i为所述相位剩余电磁数据；

F^c _i为公式(7)得到的所述岩性异常相位数据；

F^c3 _i为每一所述岩性异常相位数据F^c _i对应的3次谐波的岩性异常相位数据。

12.根据权利要求11所述的基于电磁勘探的油气检测系统，其特征在于，所述油气检测模块，根据所述振幅剩余电磁数据及相位剩余电磁数据绘制所述实验测线的拟深度域异常等值线断面图，根据拟深度域异常等值线断面图生成油气检测结果包括：

如果有已知探井数据，根据所述探井数据的油气情况确定背景值；

如果没有已知探井数据，选择剩余电磁的平均值作为背景值；

根据所述实验测线的拟深度域异常等值线断面图，当所述振幅剩余电磁数据或相位剩余电磁数据高于背景值则为异常，两者都为异常的区段为油气有利区。