CN105403924A - 黄土层静校正方法及装置 - Google Patents

Abstract

本发明提供了一种黄土层静校正方法及装置，其中，该方法包括：获取电磁数据、微测井数据和深微测井数据；获取反演电阻率；利用这些数据，得到黄土层的厚度、速度、电阻率以及反演厚度、反演速度、反演电阻率；利用上述数据，对地震资料进行静校正。本发明解决了现有技术中如果想要保证静校正的精度，就必然要增加大量成本的技术问题，达到了在保证静校正结果准确性的同时降低成本的目的。

Description

黄土层静校正方法及装置

技术领域

本发明涉及石油勘探技术领域，特别涉及一种黄土层静校正方法及装置。

背景技术

静校正问题作为地质勘探中经常出现的问题，一直是地质勘探的难点，目前常有的静校正方法主要有：小折射法、微测井法和折射静校正方法。

然而，上述几种静校正方法都存在不同层面的问题，例如，对于小折射法，在低速带调查中，由于排列的长度有限，导致观测到的初至折射波的偏移距较小，难以全面反映低速带底界的形态，因此，对地形起伏剧烈或者是表层速度变化较大的地区并不适用；对于微测井方法，能够较准确地确定测点处表层速度随深度变化的规律，但对于低速带厚度较大或者是速度变化较大的工区，必须增加测点数量，测点数量的增多必然导致成本升高；对于折射静校正方法，仅适用于地表起伏不大、表层速度横向均质性较好、存在明显折射界面且界面比较平滑的地区。

此外，还有一些其它的静校正方法，例如：基于浅层折射和统计分析的静校正方法，具体地，包括：交互迭代静校正、多域迭代静校正、模拟退火静校正等。然而，这些静校正方法对野外数据的精度有很高的要求，因此也无法广泛应用。

针对如何进行低成本准确性高的静校正，目前尚未提出有效的解决方案。

发明内容

本发明实施例提供了一种黄土层静校正方法及装置，以达到在保证静校正结果准确性的同时降低成本的目的。

本发明实施例提供了一种黄土层静校正方法，包括：采集工区的电磁数据，其中，所述工区中有M个电磁测点，所述M个电磁测点中有N个电磁测点为微测井点，有Z个电磁测点为深微测井点，M、N、Z为正整数；对所述电磁数据进行反演，得到各个电磁测点的黄土层的反演厚度和反演电阻率；根据N个微测井点的微测井数据，得到各个微测井点的黄土层厚度和黄土层速度；对Z个深微测井点进行黄土层电阻率标定；根据所述黄土层电阻率的标定结果，对各个电磁测点的反演电阻率进行追踪，得到所述各个电磁测点的黄土层厚度；检测各个微测井点的电阻率，并根据各个微测井点的黄土层速度和电阻率，建立速度与电阻率拟合关系式，并求取各个微测井点不同阶的多项式拟合系数；将各个微测井点的黄土层的反演厚度、反演电阻率和不同阶的多项式拟合系数代入所述速度与电阻率拟合关系式，得到各个微测井点不同阶的多项式拟合系数所对应的拟合速度，并根据所述拟合速度确定统一多项式拟合阶数，以求取各个微测井点在该统一多项式拟合阶数下的拟合系数；根据所述各个微测井点在该统一多项式拟合阶数下的拟合系数，进行线性插值得到所述M个电磁测点中各个电磁测点的拟合系数；将各个电磁测点的拟合系数、黄土层的反演厚度和反演电阻率，代入所述拟合关系式，得到各个电磁测点的速度；根据所述各个电磁测点的黄土层厚度和各个电磁测点的速度，进行黄土层静校正。

在一个实施例中，所述速度与电阻率拟合关系式为：

${\mathrm{Vlg}}_i=a_i\left(0\right)+a_i\left(j\right)\underset{j=1}{\overset{n_i}{\Σ}}{\left(\log 10\left({\mathrm{Rlg}}_i*D\lg \right)\right)}^j,1\≤n_i\≤5,i=1,N_s$

其中，Vlg_i表示第i个微测井点在深度Dlg处的速度，a_i表示拟合系数，Rlg_i表示第i个微测井点在深度Dlg处的电阻率，n_i表示第i个微测井点的多项式拟合总阶数，N_s表示微测井点的个数。

在一个实施例中，将各个微测井点的黄土层的反演厚度、反演电阻率和不同阶的多项式拟合系数代入所述速度与电阻率拟合关系式，得到各个微测井点不同阶的多项式拟合系数所对应的拟合速度，并根据所述拟合速度确定统一多项式拟合阶数，以求取各个微测井点在该统一多项式拟合阶数下的拟合系数，包括：将各个微测井点的黄土层的反演厚度作为微测井点的深度、反演电阻率作为在微测井点的深度处的电阻率、不同阶的多项式拟合系数作为拟合系数，代入所述速度与电阻率拟合关系式，得到各个微测井点不同阶的多项式拟合系数所对应的拟合速度；根据各个微测井点不同阶的多项式拟合系数所对应的拟合速度，计算各个微测井点不同多项式拟合阶数下的相关系数；对每个微测井，选取最大相关系数所对应的多项式拟合阶数作为该微测井的最大多项式拟合阶数，以得到N个最大多项式拟合阶数；根据所述N个最大多项式拟合阶数计算得到统一多项式拟合阶数：将所述统一多项式拟合系数作为所述N个微测井点的多项式拟合阶数，代入所述速度与电阻率拟合关系式，得到所述各个微测井点在该多项式拟合阶数下的拟合系数。

在一个实施例中，按照以下公式计算各个微测井点不同多项式拟合阶数下的相关系数：

$C_i=\frac{\underset{k=1}{\overset{Nk}{\Σ}}(V_{Ri}^k-{\stackrel{\‾}{V}}_{Ri})({\mathrm{Vlg}}_i^k-{\stackrel{\‾}{V\lg }}_i)}{\sqrt{\underset{k=1}{\overset{Nk}{\Σ}}{(V_{Ri}^k-{\stackrel{\‾}{V}}_{Ri})}^2}\sqrt{\underset{i=1}{\overset{Nk}{\Σ}}{({\mathrm{Vlg}}_i^k-{\stackrel{\‾}{V\lg }}_i)}^2}}$

其中，C_i为第i个微测井点的相关系数，Nk为每个微测井点的深度点数，为第i个微测井点第k个深度点拟合后的速度，为第i个微测井点所有深度点拟合后的速度平均值，为第i个微测井点第k个深度点的速度，为第i个微测井点所有深度点的速度平均值。

在一个实施例中，按照以下公式计算得到统一多项式拟合阶数：

$n=\underset{i=1}{\overset{N_s}{\Σ}}\frac{n_i}{N_s}$

其中，n表示所述统一多项式拟合阶数，N_s表示微测井点的个数，n_i表示第i个微测井点的最大多项式拟合阶数。

在一个实施例中，所述电磁数据包括：二维电磁数据，或者三维电磁数据；相应地，当所述电磁数据为二维电磁数据时，进行一维线性插值得到所述M个电磁测点中各个电磁测点的拟合系数；当所述电磁数据为三维电磁数据时，进行二维线性插值得到所述M个电磁测点中各个电磁测点的拟合系数。

本发明实施例还提供了一种黄土层静校正装置，包括：采集单元，采集工区的电磁数据，其中，所述工区中有M个电磁测点，所述M个电磁测点中有N个电磁测点为微测井点，有Z个电磁测点为深微测井点，M、N、Z为正整数；反演单元，对所述电磁数据进行反演，得到各个电磁测点的黄土层的反演厚度和反演电阻率；黄土层数据获取单元，根据N个微测井点的微测井数据，得到各个微测井点的黄土层厚度和黄土层速度；标定单元，对Z个深微测井点进行黄土层速度标定和电阻率标定；追踪单元，根据所述黄土层电阻率的标定结果，对各个电磁测点的反演电阻率进行追踪，得到所述各个电磁测点的黄土层厚度；拟合关系式建立单元，检测各个微测井点的电阻率，并根据各个微测井点的黄土层速度和电阻率，建立速度与电阻率拟合关系式，并求取各个微测井点不同阶的多项式拟合系数；第一拟合系数获取单元，将各个微测井点的黄土层的反演厚度、反演电阻率和不同阶的多项式拟合系数代入所述速度与电阻率拟合关系式，得到各个微测井点不同阶的多项式拟合系数所对应的拟合速度，并根据所述拟合速度确定统一多项式拟合阶数，以求取各个微测井点在该统一多项式拟合阶数下的拟合系数；第二拟合系数获取单元，根据所述各个微测井点在该统一多项式拟合阶数下的拟合系数，进行线性插值得到所述M个电磁测点中各个电磁测点的拟合系数；速度获取单元，将各个电磁测点的拟合系数、黄土层的反演厚度和反演电阻率，代入所述拟合关系式，得到各个电磁测点的速度；静校正单元，根据所述各个电磁测点的黄土层厚度和各个电磁测点的速度，进行黄土层静校正。

在一个实施例中，所述速度与电阻率拟合关系式为：

${\mathrm{Vlg}}_i=a_i\left(0\right)+a_i\left(j\right)\underset{j=1}{\overset{n_i}{\Σ}}{\left(\log 10\left({\mathrm{Rlg}}_i*D\lg \right)\right)}^j,1\≤n_i\≤5,i=1,N_s$

其中，Vlg_i表示第i个微测井点在深度Dlg处的速度，a_i表示拟合系数，Rlg_i表示第i个微测井点在深度Dlg处的电阻率，n_i表示第i个微测井点的多项式拟合总阶数，N_s表示微测井点的个数。

在一个实施例中，第一拟合系数获取单元包括：拟合速度获取模块，将各个微测井点的黄土层的反演厚度作为微测井点的深度、反演电阻率作为在微测井点的深度处的电阻率、不同阶的多项式拟合系数作为拟合系数，代入所述速度与电阻率拟合关系式，得到各个微测井点不同阶的多项式拟合系数所对应的拟合速度；相关系数计算模块，根据各个微测井点不同阶的多项式拟合系数所对应的拟合速度，计算各个微测井点不同多项式拟合阶数下的相关系数；最大多项式拟合阶数获取模块，对每个微测井，选取最大相关系数所对应的多项式拟合阶数作为该微测井的最大多项式拟合阶数，以得到N个最大多项式拟合阶数；统一多项式拟合阶数获取模块，根据所述N个最大多项式拟合阶数计算得到统一多项式拟合阶数：拟合系数获取模块，将所述统一多项式拟合系数作为所述N个微测井点的多项式拟合阶数，代入所述速度与电阻率拟合关系式，得到所述各个微测井点在该多项式拟合阶数下的拟合系数。

在一个实施例中，所述相关系数计算模块具体用于按照以下公式计算各个微测井点不同多项式拟合阶数下的相关系数：

$C_i=\frac{\underset{k=1}{\overset{Nk}{\Σ}}(V_{Ri}^k-{\stackrel{\‾}{V}}_{Ri})({\mathrm{Vlg}}_i^k-{\stackrel{\‾}{V\lg }}_i)}{\sqrt{\underset{k=1}{\overset{Nk}{\Σ}}{(V_{Ri}^k-{\stackrel{\‾}{V}}_{Ri})}^2}\sqrt{\underset{i=1}{\overset{Nk}{\Σ}}{({\mathrm{Vlg}}_i^k-{\stackrel{\‾}{V\lg }}_i)}^2}}$

其中，C_i为第i个微测井点的相关系数，Nk为每个微测井点的深度点数，为第i个微测井点第k个深度点拟合后的速度，为第i个微测井点所有深度点拟合后的速度平均值，为第i个微测井点第k个深度点的速度，为第i个微测井点所有深度点的速度平均值。

在上述实施例中，在计算不同深度的黄土塬对应地层速度的过程中，在原有的电磁测点的基础上布设了微测井点和深微测井点，然后基于电磁测点的电磁数据和微测井点的微测井数据，确定各个微测井点在统一多项式拟合阶数下的拟合系数，然后基于这些微测井点在统一多项式拟合阶数下的拟合系数进行插值以得到各个电磁测点的拟合系数，从而计算出各个电测测点的厚度和速度信息，以实现对黄土层的静校正。通过上述方式解决了如有想要保证静校正的精度，就必然要增加大量成本的技术问题，达到了在保证静校正结果准确性的同时降低成本的目的。

附图说明

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

图1是本发明实施例的黄土层静校正的方法流程图；

图2是本发明实施例的二维音频大地电磁勘探数据和微测井数据结合进行二维地震剖面的黄土层静校正的方法流程图；

图3是本发明实施例的黄土层静校正的装置结构框图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施方式和附图，对本发明做进一步详细说明。在此，本发明的示意性实施方式及其说明用于解释本发明，但并不作为对本发明的限定。

考虑到现有的基于微测井方法进行静校正时，如果需要较准确地确定测点处表层速度随深度变化的规律，就需要增加大量测点，如果测点太少，则无法为地震资料提供可靠的表层模型。在本例中，主要是考虑到可以按照一定的规则增加测点，从而使得在保证准确度的情况下，尽可能少的增加测点，从而降低成本。具体地，通过在已有的微测井点的基础上，引入了电磁测点，并进行一系列处理后，可以对工区内的地震资料进行更加准确的黄土层静校正。具体地，在本例中，提供了一种黄土层静校正方法，如图1所示，可以包括以下步骤：

步骤101：采集工区的电磁数据，其中，所述工区中有M个电磁测点，所述M个电磁测点中有N个电磁测点为微测井点，有Z个电磁测点为深微测井点，M、N、Z为正整数；

在实际应用中，可以根据地震工区的勘探范围和地震资料对黄土层静校正数据密度的要求，设计并采集相应密度的电磁数据。其中，相应密度的电磁数据可以包括：音频大地电磁、瞬变电磁、大地电磁、连续电磁剖面，这些数据可以是二维电磁数据，也可以是三维电磁数据。

上述工区中的电磁测点可以按照以下规则布设：根据电磁测点的十分之一到二十分之一的密度，在地震工区内布设常规微测井，在布设的过程中，各个微测井的深度可以根据实际要求确定。布设的这些常规微测井在工区内相对均匀分布。具体地，当电磁数据为二维测线时，可以在测线的起点和终点处各布设一口微测井，当电磁数据为三维测网时，可以在测网的拐角和边界处布设一定数量的微测井。进一步的，还可以在工区内布设2口钻穿黄土层的深微测井。

步骤102：对所述电磁数据进行反演，得到各个电磁测点的黄土层的反演厚度和反演电阻率；

具体地，当电磁数据为二维电磁数据时，对该数据进行二维反演处理，当电磁数据为三维电磁数据时，对该数据进行三维反演处理，从而获得地震工区内各个电磁测点的黄土层的反演厚度和反演电阻率。

步骤103：根据N个微测井点的微测井数据，得到各个微测井点的黄土层厚度和黄土层速度；

具体地，可以采用时距图法求取各个微测井点的黄土层厚度和黄土层速度。

时距图表示地震波到达各观测点的时间与震源至观测点间距离关系的图形。它是运动地震学研究的基本内容。对时距图的分析和处理，可以对地质得出定性定量的解释。

步骤104：对Z个深微测井点进行黄土层电阻率标定；其中，这Z个深微测井点是钻穿黄土层的深微测井。

步骤105：根据所述黄土层电阻率的标定结果，对各个电磁测点的反演电阻率进行追踪，得到所述各个电磁测点的黄土层厚度；

步骤106：检测各个微测井点的电阻率，并根据各个微测井点的黄土层速度和电阻率，建立速度与电阻率拟合关系式，并求取各个微测井点不同阶的多项式拟合系数；

具体的，根据各个微测井点的黄土层速度和电阻率，所建立的速度与电阻率拟合关系式如下：

${\mathrm{Vlg}}_i=a_i\left(0\right)+a_i\left(j\right)\underset{j=1}{\overset{n_i}{\Σ}}{\left(\log 10\left({\mathrm{Rlg}}_i*D\lg \right)\right)}^j,1\≤n_i\≤5,i=1,N_s$ (公式1)

其中，Vlg_i表示第i个微测井点在深度Dlg处的速度，a_i表示拟合系数，Rlg_i表示第i个微测井点在深度Dlg处的电阻率，n_i表示第i个微测井点的多项式拟合总阶数，N_s表示微测井点的个数。

将步骤103所得到的各个微测井点的黄土层速度和黄土层厚度，以及步骤104中深微测井点所得到的各个微测井点的黄土层的电阻率代入所建立的速度与电阻率拟合关系式，可以得到地震工区内每个微测井点所对应的1至5阶多项式拟合系数a_i。

步骤107：将各个微测井点的黄土层的反演厚度、反演电阻率和不同阶的多项式拟合系数代入所述速度与电阻率拟合关系式，得到各个微测井点不同阶的多项式拟合系数所对应的拟合速度，并根据所述拟合速度确定统一多项式拟合阶数，以求取各个微测井点在该统一多项式拟合阶数下的拟合系数；

具体地，可以按照以下方式求取各个微测井点在统一多项式拟合阶数下的拟合系数，包括：

步骤1：将各个微测井点的黄土层的反演厚度作为微测井点的深度、反演电阻率作为在微测井点的深度处的电阻率、不同阶的多项式拟合系数作为拟合系数，代入速度与电阻率拟合关系式(即上述公式1)，得到各个微测井点不同阶的多项式拟合系数所对应的拟合速度；

通过将N个微测井点在步骤102中所得到的黄土层反演厚度和反演电阻率以及步骤106中所得到的地震工区内各个微测井点所对应的1至5阶多项式拟合系数a_i代入公式1，可以得到各个微测井点1至5阶多项式拟合系数所对应的拟合速度。

步骤2：根据各个微测井点不同阶的多项式拟合系数所对应的拟合速度，计算各个微测井点不同多项式拟合阶数下的相关系数；

具体地，可以按照以下公式计算各个微测井点不同多项式拟合阶数下的相关系数：

$C_i=\frac{\underset{k=1}{\overset{Nk}{\Σ}}(V_{Ri}^k-{\stackrel{\‾}{V}}_{Ri})({\mathrm{Vlg}}_i^k-{\stackrel{\‾}{V\lg }}_i)}{\sqrt{\underset{k=1}{\overset{Nk}{\Σ}}{(V_{Ri}^k-{\stackrel{\‾}{V}}_{Ri})}^2}\sqrt{\underset{i=1}{\overset{Nk}{\Σ}}{({\mathrm{Vlg}}_i^k-{\stackrel{\‾}{V\lg }}_i)}^2}}$ (公式2)

其中，C_i为第i个微测井点的相关系数，Nk为每个微测井点的深度点数，为第i个微测井点第k个深度点拟合后的速度，为第i个微测井点所有深度点拟合后的速度平均值，为第i个微测井点第k个深度点的速度，为第i个微测井点所有深度点的速度平均值。

通过将步骤1所得到的拟合速度代入公式2中，可以得到各个微测井点1至5阶多项式拟合阶数下的相关系数。

步骤3：对每个微测井，选取最大相关系数所对应的多项式拟合阶数作为该微测井的最大多项式拟合阶数，以得到N个最大多项式拟合阶数；

在本实施例中，对于N个微测井中的每个微测井均有1至5阶共5个多项式拟合阶数，对于每个微测井都选取其最大相关系数所对应的多项式的拟合阶数，作为该微测井的最大多项式拟合阶数，也就相当于在该步骤中会得到N个最大多项式拟合阶数，每个微测井有一个最大多项式拟合阶数。

步骤4：根据所述N个最大多项式拟合阶数计算得到统一多项式拟合阶数；

具体地，可以按照以下公式计算统一多项式拟合阶数：

$n=\underset{i=1}{\overset{N_s}{\Σ}}\frac{n_i}{N_s}$ (公式3)

其中，n表示所述统一多项式拟合阶数，N_s表示微测井点的个数，n_i表示第i个微测井点的最大多项式拟合阶数。

步骤5：将所述统一多项式拟合阶数作为所述N个微测井点的多项式拟合阶数，代入所述速度与电阻率拟合关系式，得到所述各个微测井点在该多项式拟合阶数下的拟合系数。

具体的，将步骤4中所得到的统一多项式拟合阶数作为一个固定值代入公式1，并分别将步骤103所得到的各个微测井点的黄土层速度和黄土层厚度，以及步骤104中深微测井点所得到的各个微测井点的黄土层的电阻率代入公式1即速度与电阻率拟合关系式，可以得到地震工区内每个微测井点所对应的在该统一多项式拟合阶数下的拟合系数a_i。

步骤108：根据所述各个微测井点在该统一多项式拟合阶数下的拟合系数，进行线性插值得到所述M个电磁测点中各个电磁测点的拟合系数；

例如，N取值为10，M取值为100，Z的取值为2，就相当于通过2个深微测井点，以及微测井数据，计算得到10个微测井点的拟合系数，然后基于这10个微测井点在统一的多项式阶数下的拟合系数进行插值，得到100个电磁测点的拟合系数。之所采用这种方式，主要是考虑到仅采用微测井方法进行静校正时，由于测点较少，因而无法为地震资料提供可靠的表层模型。本实施例在已有的N个微测井点的基础上，引入了Z个深微测井点，根据已有的N个微测井点在统一多项式拟合阶数下的拟合系数，对其进行线性插值，从而可以得到本实施例中的M个电磁测点中各个电磁测点的拟合系数。

具体地，在插值的过程中，如果电磁数据为二维电磁数据，那么就可以进行一维线性插值得到M个电磁测点中各个电磁测点的拟合系数，如果电磁数据为三维电磁数据，那么就可以进行二维线性插值得到M个电磁测点中各个电磁测点的拟合系数。

步骤109：将各个电磁测点的拟合系数、黄土层的反演厚度和反演电阻率，代入所述拟合关系式，得到各个电磁测点的速度；

即，在得到M个电磁测点中各个电磁测点的拟合系数的基础上，将由步骤102所得到的黄土层的反演厚度和反演电阻率，代入公式1即拟合关系式，从而可以得到M个电磁测点的速度；

步骤110：根据所述各个电磁测点的黄土层厚度和各个电磁测点的速度，进行黄土层静校正。

通过对M个电磁测点的黄土层厚度和黄土层速度进行分析，从而可以对工区内的地震资料进行更加准确的黄土层静校正。

下面以一个上述黄土静校正的方法的具体应用实例为例来进行说明，在本例中，以二维音频大地电磁(Audio-frequencyMagnetotellurics，简称为AMT)勘探数据和微测井数据结合来进行二维地震剖面的黄土静校正为例进行说明。然而值得注意的是，该具体实施例仅是为了更好地说明本发明，并不构成对本发明的不当限定。

如图2所示，该具体实例包括以下步骤：

步骤201：在二维地震测线上设计并采集相应密度的二维AMT电磁测点；

步骤202：在该二维测线上布设常规微测井点，常规微测井点的数量可以是二维AMT电磁测点的十分之一，常规微测井点的深度可以根据具体要求设计，且布设的常规微测井点在工区内相对均匀分布并与AMT电测测点位置重合。进一步的，还需在二维测线的起点和终点各布设一口常规微测井，同时在该二维测线的三分之一和三分之二处分别布设2口钻穿黄土层的深微测井；

步骤203：反演处理二维AMT电磁数据，得到工区内二维测线的黄土层的反演厚度和反演电阻率；

步骤204：利用微测井数据，得到各个微测井点的黄土层厚度和黄土层速度；

步骤205：对深微测井点进行黄土层电阻率标定；

步骤206：根据所述黄土层电阻率的标定结果，对二维测线上各个电磁测点的反演电阻率进行追踪，得到所述各个电磁测点的黄土层厚度；

步骤207：检测各个微测井点的电阻率，利用各个微测井点的黄土层速度、电阻率，反演厚度、反演电阻率，计算得到各个微测井点在统一多项式拟合阶数下的拟合系数；

步骤208：沿着二维AMT电磁测点位置，对步骤207中所得到的拟合系数进行一维线性插值，可以得到各个电磁测点的拟合系数；

步骤209：将各个电磁测点的拟合系数、黄土层的反演厚度和反演电阻率，代入拟合关系式，可以得到各个电磁测点的速度；

步骤210：根据各个电磁测点的黄土层厚度和各个电磁测点的速度，对该电磁测线相应的地震资料进行黄土层静校正。

基于同一发明构思，本发明实施例中还提供了一种黄土层静校正的装置，如下面的实施例所述。由于黄土层静校正的装置解决问题的原理与黄土层静校正的方法相似，因此黄土层静校正的装置的实施可以参见黄土层静校正的方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。图3是本发明实施例的黄土层静校正的装置的一种结构框图，如图3所示，包括：采集单元301、反演单元302、黄土层数据获取单元303、标定单元304、追踪单元305、拟合关系式建立单元306、第一拟合系数获取单元307、第二拟合系数获取单元308、速度获取单元309、静校正单元310，下面对该结构进行说明。

采集单元，用于采集工区的电磁数据，其中，所述工区中有M个电磁测点，所述M个电磁测点中有N个电磁测点为微测井点，有Z个电磁测点为深微测井点，M、N、Z为正整数；

反演单元，用于对所述电磁数据进行反演，得到各个电磁测点的黄土层的反演厚度和反演电阻率；

黄土层数据获取单元，用于根据N个微测井点的微测井数据，得到各个微测井点的黄土层厚度和黄土层速度；

标定单元，用于对Z个深微测井点进行黄土层速度标定和电阻率标定；

追踪单元，用于根据所述黄土层电阻率的标定结果，对各个电磁测点的反演电阻率进行追踪，得到所述各个电磁测点的黄土层厚度；

拟合关系式建立单元，用于检测各个微测井点的电阻率，并根据各个微测井点的黄土层速度和电阻率，建立速度与电阻率拟合关系式，并求取各个微测井点不同阶的多项式拟合系数；

第一拟合系数获取单元，用于将各个微测井点的黄土层的反演厚度、反演电阻率和不同阶的多项式拟合系数代入所述速度与电阻率拟合关系式，得到各个微测井点不同阶的多项式拟合系数所对应的拟合速度，并根据所述拟合速度确定统一多项式拟合阶数，以求取各个微测井点在该统一多项式拟合阶数下的拟合系数；

第二拟合系数获取单元，用于根据所述各个微测井点在该统一多项式拟合阶数下的拟合系数，进行线性插值得到所述M个电磁测点中各个电磁测点的拟合系数；

速度获取单元，用于将各个电磁测点的拟合系数、黄土层的反演厚度和反演电阻率，代入所述拟合关系式，得到各个电磁测点的速度；

静校正单元，用于根据所述各个电磁测点的黄土层厚度和各个电磁测点的速度，进行黄土层静校正。

在一个实施方式中，所述速度与电阻率拟合关系式为：

${\mathrm{Vlg}}_i=a_i\left(0\right)+a_i\left(j\right)\underset{j=1}{\overset{n_i}{\Σ}}{\left(\log 10\left({\mathrm{Rlg}}_i*D\lg \right)\right)}^j,1\≤n_i\≤5,i=1,N_s$

其中，Vlg_i表示第i个微测井点在深度Dlg处的速度，a_i表示拟合系数，Rlg_i表示第i个微测井点在深度Dlg处的电阻率，n_i表示第i个微测井点的多项式拟合总阶数，N_s表示微测井点的个数。

在一个实施方式中，第一拟合系数获取单元包括：拟合速度获取模块，用于将各个微测井点的黄土层的反演厚度作为微测井点的深度、反演电阻率作为在微测井点的深度处的电阻率、不同阶的多项式拟合系数作为拟合系数，代入所述速度与电阻率拟合关系式，得到各个微测井点不同阶的多项式拟合系数所对应的拟合速度；相关系数计算模块，用于根据各个微测井点不同阶的多项式拟合系数所对应的拟合速度，计算各个微测井点不同多项式拟合阶数下的的相关系数；最大多项式拟合阶数获取模块，用于对每个微测井，选取最大相关系数所对应的多项式拟合阶数作为该微测井的最大多项式拟合阶数，以得到N个最大多项式拟合阶数；统一多项式拟合阶数获取模块，用于根据所述N个最大多项式拟合阶数计算得到统一多项式拟合阶数：拟合系数获取模块，用于将所述统一多项式拟合系数作为所述N个微测井点的多项式拟合阶数，代入所述速度与电阻率拟合关系式，得到所述各个微测井点在该多项式拟合阶数下的拟合系数。

在一个实施方式中，所述相关系数计算模块具体用于按照以下公式计算各个微测井点不同多项式拟合阶数下的相关系数：

$C_i=\frac{\underset{k=1}{\overset{Nk}{\Σ}}(V_{Ri}^k-{\stackrel{\‾}{V}}_{Ri})({\mathrm{Vlg}}_i^k-{\stackrel{\‾}{V\lg }}_i)}{\sqrt{\underset{k=1}{\overset{Nk}{\Σ}}{(V_{Ri}^k-{\stackrel{\‾}{V}}_{Ri})}^2}\sqrt{\underset{i=1}{\overset{Nk}{\Σ}}{({\mathrm{Vlg}}_i^k-{\stackrel{\‾}{V\lg }}_i)}^2}}$

其中，C_i为第i个微测井点的相关系数，Nk为每个微测井点的深度点数，为第i个微测井点第k个深度点拟合后的速度，为第i个微测井点所有深度点拟合后的速度平均值，为第i个微测井点第k个深度点的速度，为第i个微测井点所有深度点的速度平均值。

在一个实施方式中，所述统一多项式拟合阶数获取模块具体用于按照以下公式计算得到统一多项式拟合阶数：

$n=\underset{i=1}{\overset{N_s}{\Σ}}\frac{n_i}{N_s}$

其中，n表示所述统一多项式拟合阶数，N_s表示微测井点的个数，n_i表示第i个微测井点的最大多项式拟合阶数。

在一个实施方式中，所述电磁数据包括：二维电磁数据，或者三维电磁数据；相应地，当所述电磁数据为二维电磁数据时，进行一维线性插值得到所述M个电磁测点中各个电磁测点的拟合系数；当所述电磁数据为三维电磁数据时，进行二维线性插值得到所述M个电磁测点中各个电磁测点的拟合系数。

在另外一个实施例中，还提供了一种软件，该软件用于执行上述实施例及优选实施方式中描述的技术方案。

在另外一个实施例中，还提供了一种存储介质，该存储介质中存储有上述软件，该存储介质包括但不限于：光盘、软盘、硬盘、可擦写存储器等。

从以上的描述中，可以看出，本发明实施例实现了如下技术效果：在原有的电磁测点的基础上布设了微测井点和深微测井点，然后基于电磁测点的电磁数据和微测井点的微测井数据，确定各个微测井点在统一多项式拟合阶数下的拟合系数，然后基于这些微测井点在统一多项式拟合阶数下的拟合系数进行插值以得到各个电磁测点的拟合系数，从而计算出各个电测测点的厚度和速度信息，以实现对黄土层的静校正，从而有效解决了现有技术中基于微测井方法进行静校正时，由于测点较少所产生的无法为地震资料提供可靠的表层模型的问题，进一步的，通过将电磁数据和微测井数据结合，有效地解决了黄土塬发育地区的黄土层厚度分布及速度的探测问题，为地震资料提供了更加可靠的表层模型，最终得到的黄土层静校正结果也更为符合地质构造。

显然，本领域的技术人员应该明白，上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，并且在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明实施例不限制于任何特定的硬件和软件结合。

以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种黄土层静校正方法，其特征在于，包括：

采集工区的电磁数据，其中，所述工区中有M个电磁测点，所述M个电磁测点中有N个电磁测点为微测井点，有Z个电磁测点为深微测井点，M、N、Z为正整数；

对所述电磁数据进行反演，得到各个电磁测点的黄土层的反演厚度和反演电阻率；

根据N个微测井点的微测井数据，得到各个微测井点的黄土层厚度和黄土层速度；

对Z个深微测井点进行黄土层电阻率标定；

根据所述黄土层电阻率的标定结果，对各个电磁测点的反演电阻率进行追踪，得到所述各个电磁测点的黄土层厚度；

检测各个微测井点的电阻率，并根据各个微测井点的黄土层速度和电阻率，建立速度与电阻率拟合关系式，并求取各个微测井点不同阶的多项式拟合系数；

将各个微测井点的黄土层的反演厚度、反演电阻率和不同阶的多项式拟合系数代入所述速度与电阻率拟合关系式，得到各个微测井点不同阶的多项式拟合系数所对应的拟合速度，并根据所述拟合速度确定统一多项式拟合阶数，以求取各个微测井点在该统一多项式拟合阶数下的拟合系数；

根据所述各个微测井点在该统一多项式拟合阶数下的拟合系数，进行线性插值得到所述M个电磁测点中各个电磁测点的拟合系数；

将各个电磁测点的拟合系数、黄土层的反演厚度和反演电阻率，代入所述拟合关系式，得到各个电磁测点的速度；

根据所述各个电磁测点的黄土层厚度和各个电磁测点的速度，进行黄土层静校正。

2.如权利要求1所述的方法，其特征在于，所述速度与电阻率拟合关系式为：

${\mathrm{Vlg}}_i=a_i\left(0\right)+a_i\left(j\right)\underset{j=1}{\overset{n_i}{\Σ}}{\left(\log 10\left({\mathrm{Rlg}}_i*D\lg \right)\right)}^j,1\≤n_i\≤5,i=1,N_s$

其中，Vlg_i表示第i个微测井点在深度Dlg处的速度，a_i表示拟合系数，Rlg_i表示第i个微测井点在深度Dlg处的电阻率，n_i表示第i个微测井点的多项式拟合总阶数，N_s表示微测井点的个数。

3.如权利要求1所述的方法，其特征在于，将各个微测井点的黄土层的反演厚度、反演电阻率和不同阶的多项式拟合系数代入所述速度与电阻率拟合关系式，得到各个微测井点不同阶的多项式拟合系数所对应的拟合速度，并根据所述拟合速度确定统一多项式拟合阶数，以求取各个微测井点在该统一多项式拟合阶数下的拟合系数，包括：

将各个微测井点的黄土层的反演厚度作为微测井点的深度、反演电阻率作为在微测井点的深度处的电阻率、不同阶的多项式拟合系数作为拟合系数，代入所述速度与电阻率拟合关系式，得到各个微测井点不同阶的多项式拟合系数所对应的拟合速度；

根据各个微测井点不同阶的多项式拟合系数所对应的拟合速度，计算各个微测井点不同多项式拟合阶数下的相关系数；

对每个微测井，选取最大相关系数所对应的多项式拟合阶数作为该微测井的最大多项式拟合阶数，以得到N个最大多项式拟合阶数；

根据所述N个最大多项式拟合阶数计算得到统一多项式拟合阶数：

将所述统一多项式拟合系数作为所述N个微测井点的多项式拟合阶数，代入所述速度与电阻率拟合关系式，得到所述各个微测井点在该多项式拟合阶数下的拟合系数。

4.如权利要求3所述的方法，其特征在于，按照以下公式计算各个微测井点不同多项式拟合阶数下的相关系数：

$C_i=\frac{\underset{k=1}{\overset{Nk}{\Σ}}\left(V_{Ri}^k-{\stackrel{\‾}{V}}_{Ri}\right)\left({\mathrm{Vlg}}_i^k-{\stackrel{\‾}{V\lg }}_i\right)}{\sqrt{\underset{k=1}{\overset{Nk}{\Σ}}{\left(V_{Ri}^k-{\stackrel{\‾}{V}}_{Ri}\right)}^2}\sqrt{\underset{i=1}{\overset{Nk}{\Σ}}{\left({\mathrm{Vlg}}_i^k-{\stackrel{\‾}{V\lg }}_i\right)}^2}}$

其中，C_i为第i个微测井点的相关系数，Nk为每个微测井点的深度点数，为第i个微测井点第k个深度点拟合后的速度，为第i个微测井点所有深度点拟合后的速度平均值，为第i个微测井点第k个深度点的速度，为第i个微测井点所有深度点的速度平均值。

5.如权利要求3所述的方法，其特征在于，按照以下公式计算得到统一多项式拟合阶数：

$n=\underset{i=1}{\overset{N_s}{\Σ}}\frac{n_i}{N_s}$

其中，n表示所述统一多项式拟合阶数，N_s表示微测井点的个数，n_i表示第i个微测井点的最大多项式拟合阶数。

6.如权利要求1至5中任一项所述的方法，其特征在于，所述电磁数据包括：二维电磁数据，或者三维电磁数据；

相应地，

当所述电磁数据为二维电磁数据时，进行一维线性插值得到所述M个电磁测点中各个电磁测点的拟合系数；

当所述电磁数据为三维电磁数据时，进行二维线性插值得到所述M个电磁测点中各个电磁测点的拟合系数。

7.一种黄土层静校正装置，其特征在于，包括：

采集单元，用于采集工区的电磁数据，其中，所述工区中有M个电磁测点，所述M个电磁测点中有N个电磁测点为微测井点，有Z个电磁测点为深微测井点，M、N、Z为正整数；

反演单元，用于对所述电磁数据进行反演，得到各个电磁测点的黄土层的反演厚度和反演电阻率；

黄土层数据获取单元，用于根据N个微测井点的微测井数据，得到各个微测井点的黄土层厚度和黄土层速度；

标定单元，用于对Z个深微测井点进行黄土层速度标定和电阻率标定；

追踪单元，用于根据所述黄土层电阻率的标定结果，对各个电磁测点的反演电阻率进行追踪，得到所述各个电磁测点的黄土层厚度；

拟合关系式建立单元，用于检测各个微测井点的电阻率，并根据各个微测井点的黄土层速度和电阻率，建立速度与电阻率拟合关系式，并求取各个微测井点不同阶的多项式拟合系数；

第一拟合系数获取单元，用于将各个微测井点的黄土层的反演厚度、反演电阻率和不同阶的多项式拟合系数代入所述速度与电阻率拟合关系式，得到各个微测井点不同阶的多项式拟合系数所对应的拟合速度，并根据所述拟合速度确定统一多项式拟合阶数，以求取各个微测井点在该统一多项式拟合阶数下的拟合系数；

第二拟合系数获取单元，用于根据所述各个微测井点在该统一多项式拟合阶数下的拟合系数，进行线性插值得到所述M个电磁测点中各个电磁测点的拟合系数；

速度获取单元，用于将各个电磁测点的拟合系数、黄土层的反演厚度和反演电阻率，代入所述拟合关系式，得到各个电磁测点的速度；

静校正单元，用于根据所述各个电磁测点的黄土层厚度和各个电磁测点的速度，进行黄土层静校正。

8.如权利要求7所述的装置，其特征在于，所述速度与电阻率拟合关系式为：

${\mathrm{Vlg}}_i=a_i\left(0\right)+a_i\left(j\right)\underset{j=1}{\overset{n_i}{\Σ}}{\left(\log 10\left({\mathrm{Rlg}}_i*D\lg \right)\right)}^j,1\≤n_i\≤5,i=1,N_s$

其中，Vlg_i表示第i个微测井点在深度Dlg处的速度，a_i表示拟合系数，Rlg_i表示第i个微测井点在深度Dlg处的电阻率，n_i表示第i个微测井点的多项式拟合总阶数，N_s表示微测井点的个数。

9.如权利要求7所述的装置，其特征在于，第一拟合系数获取单元包括：

拟合速度获取模块，用于将各个微测井点的黄土层的反演厚度作为微测井点的深度、反演电阻率作为在微测井点的深度处的电阻率、不同阶的多项式拟合系数作为拟合系数，代入所述速度与电阻率拟合关系式，得到各个微测井点不同阶的多项式拟合系数所对应的拟合速度；

相关系数计算模块，用于根据各个微测井点不同阶的多项式拟合系数所对应的拟合速度，计算各个微测井点不同多项式拟合阶数下的相关系数；

最大多项式拟合阶数获取模块，用于对每个微测井，选取最大相关系数所对应的多项式拟合阶数作为该微测井的最大多项式拟合阶数，以得到N个最大多项式拟合阶数；

统一多项式拟合阶数获取模块，用于根据所述N个最大多项式拟合阶数计算得到统一多项式拟合阶数；

拟合系数获取模块，用于将所述统一多项式拟合系数作为所述N个微测井点的多项式拟合阶数，代入所述速度与电阻率拟合关系式，得到所述各个微测井点在该多项式拟合阶数下的拟合系数。

10.如权利要求9所述的装置，其特征在于，所述相关系数计算模块具体用于按照以下公式计算各个微测井点不同多项式拟合阶数下的相关系数：

$C_i=\frac{\underset{k=1}{\overset{Nk}{\Σ}}\left(V_{Ri}^k-{\stackrel{\‾}{V}}_{Ri}\right)\left({\mathrm{Vlg}}_i^k-{\stackrel{\‾}{V\lg }}_i\right)}{\sqrt{\underset{k=1}{\overset{Nk}{\Σ}}{\left(V_{Ri}^k-{\stackrel{\‾}{V}}_{Ri}\right)}^2}\sqrt{\underset{i=1}{\overset{Nk}{\Σ}}{\left({\mathrm{Vlg}}_i^k-{\stackrel{\‾}{V\lg }}_i\right)}^2}}$

其中，C_i为第i个微测井点的相关系数，Nk为每个微测井点的深度点数，为第i个微测井点第k个深度点拟合后的速度，为第i个微测井点所有深度点拟合后的速度平均值，为第i个微测井点第k个深度点的速度，为第i个微测井点所有深度点的速度平均值。