技術(shù)特征：

1.一種薄儲層的識別方法，其特征在于，所述方法包括：

獲取待識別儲層的原始地震道信號S；

分別計算所述原始地震道信號S的第2階至第n階導(dǎo)數(shù)，得到n-1個導(dǎo)數(shù)：S²，...Sⁱ，...，Sⁿ，其中，Sⁱ表示所述原始地震道信號S的第i階導(dǎo)數(shù)，所述n為根據(jù)所述原始地震道信號S的導(dǎo)數(shù)確定的疊加閾值，所述i為大于等于2，且小于等于所述n的整數(shù)；

根據(jù)所述原始地震道信號S和所述n-1個導(dǎo)數(shù)：S²，...Sⁱ，...，Sⁿ，通過總體經(jīng)驗?zāi)B(tài)分解CEEMD算法確定目標(biāo)地震道信號S’，所述CEEMD算法為：

$S^{\′}=CEEMD\left\{Nor\left(Nor\right(S)+Gau(Nor\left(S\right))-\underset{i=1}{\overset{M}{\Σ}}Nor\left(S^{4i-2}\right)+\underset{k=1}{\overset{N}{\Σ}}Nor(S^{4k}\left)\right)\right\};$

其中，Gau( )表示對括號內(nèi)的內(nèi)容進行高斯變換，Nor( )表示對括號內(nèi)的內(nèi)容進行歸一化處理，M為預(yù)設(shè)的第一閾值，N為預(yù)設(shè)的第二閾值，所述第一閾值M為滿足條件：4M-2≤n中的最大值，所述第二閾值其中表示向下取整；

根據(jù)所述目標(biāo)地震道信號S’，識別所述待識別儲層中的薄儲層。

2.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述分別計算所述原始地震道信號S的第2階至第n階導(dǎo)數(shù)，得到n-1個導(dǎo)數(shù)：S²，...Sⁱ，...，Sⁿ，包括：

當(dāng)j≥2時，判斷待分析導(dǎo)數(shù)對應(yīng)的主頻D^j是否滿足循環(huán)終止條件，所述待分析導(dǎo)數(shù)包括：所述原始地震道信號S的第j階導(dǎo)數(shù)S^j，第j+1階導(dǎo)數(shù)S^j+1和第j-1階導(dǎo)數(shù)S^j-1，所述循環(huán)終止條件為：

$\frac{D^{j+1}-D^j}{D^j-D^{j-1}}\≤\mathrm{\ϵ};$

其中，ε為預(yù)設(shè)的終止參數(shù)，D^j，D^j-1和D^j+1分別為根據(jù)主頻計算公式確定的所述原始地震道信號S的第j階導(dǎo)數(shù)S^j的主頻，所述原始地震道信號S的第j+1階導(dǎo)數(shù)S^j+1的主頻，以及所述原始地震道信號S的第j-1階導(dǎo)數(shù)S^j-1的主頻，所述主頻計算公式為：

$D^j=\frac{{\∫}_0^{\mathrm{\∞}}{\mathrm{fS}}^j\left(f\right)df}{{\∫}_0^{\mathrm{\∞}}{S}^j\left(f\right)df};$

其中，f為頻率變量，S^j(f)為所述原始地震道信號S的第j階導(dǎo)數(shù)S^j的頻譜函數(shù)；

當(dāng)所述述待分析導(dǎo)數(shù)對應(yīng)的主頻D^j滿足循環(huán)終止條件時，將j+2得到新的j，重復(fù)執(zhí)行上述判斷過程，直至所述待分析導(dǎo)數(shù)對應(yīng)的主頻D^j不滿足所述循環(huán)終止條件，將不滿足所述循環(huán)終止條件的待分析導(dǎo)數(shù)對應(yīng)的j確定為所述疊加閾值。

3.根據(jù)權(quán)利要求2所述的方法，其特征在于，

所述預(yù)設(shè)的終止參數(shù)ε的取值范圍為0.2至0.3。

4.根據(jù)權(quán)利要求1至3任一所述的方法，其特征在于，

所述Nor( )表示對括號內(nèi)的內(nèi)容進行歸一化處理后得到的數(shù)值在區(qū)域[-1，1]內(nèi)。

5.一種薄儲層的識別裝置，其特征在于，所述裝置包括：

獲取模塊，用于獲取待識別儲層的原始地震道信號S；

計算模塊，用于分別計算所述原始地震道信號S的第2階至第n階導(dǎo)數(shù)，得到n-1個導(dǎo)數(shù)：S²，...Sⁱ，...，Sⁿ，其中，Sⁱ表示所述原始地震道信號S的第i階導(dǎo)數(shù)，所述n為根據(jù)所述原始地震道信號S的導(dǎo)數(shù)確定的疊加閾值，所述i為大于等于2，且小于等于所述n的整數(shù)；

確定模塊，用于根據(jù)所述原始地震道信號S和所述n-1個導(dǎo)數(shù)：S²，...Sⁱ，...，Sⁿ，通過總體經(jīng)驗?zāi)B(tài)分解CEEMD算法確定目標(biāo)地震道信號S’，所述CEEMD算法為：

$S^{\′}=CEEMD\left\{Nor\left(Nor\right(S)+Gau(Nor\left(S\right))-\underset{i=1}{\overset{M}{\Σ}}Nor\left(S^{4i-2}\right)+\underset{k=1}{\overset{N}{\Σ}}Nor(S^{4k}\left)\right)\right\};$

其中，Gau( )表示對括號內(nèi)的內(nèi)容進行高斯變換，Nor( )表示對括號內(nèi)的內(nèi)容進行歸一化處理，M為預(yù)設(shè)的第一閾值，N為預(yù)設(shè)的第二閾值，所述第一閾值M為滿足條件：4M-2≤n中的最大值，所述第二閾值其中表示向下取整；

識別模塊，用于根據(jù)所述目標(biāo)地震道信號S’，識別所述待識別儲層中的薄儲層。

6.根據(jù)權(quán)利要求5所述的裝置，其特征在于，所述計算模塊，還用于：

當(dāng)j≥2時，判斷待分析導(dǎo)數(shù)對應(yīng)的主頻D^j是否滿足循環(huán)終止條件，所述待分析導(dǎo)數(shù)包括：所述原始地震道信號S的第j階導(dǎo)數(shù)S^j，第j+1階導(dǎo)數(shù)S^j+1和第j-1階導(dǎo)數(shù)S^j-1，所述循環(huán)終止條件為：

$\frac{D^{j+1}-D^j}{D^j-D^{j-1}}\≤\mathrm{\ϵ};$

其中，ε為預(yù)設(shè)的終止參數(shù)，D^j，D^j-1和D^j+1分別為根據(jù)主頻計算公式確定的所述原始地震道信號S的第j階導(dǎo)數(shù)S^j的主頻，所述原始地震道信號S的第j+1階導(dǎo)數(shù)S^j+1的主頻，以及所述原始地震道信號S的第j-1階導(dǎo)數(shù)S^j-1的主頻，所述主頻計算公式為：

$D^j=\frac{{\∫}_0^{\mathrm{\∞}}{\mathrm{fS}}^j\left(f\right)df}{{\∫}_0^{\mathrm{\∞}}{S}^j\left(f\right)df};$

其中，f為頻率變量，S^j(f)為所述原始地震道信號S的第j階導(dǎo)數(shù)S^j的頻譜函數(shù)；

當(dāng)所述述待分析導(dǎo)數(shù)對應(yīng)的主頻D^j滿足循環(huán)終止條件時，將j+2得到新的j，重復(fù)執(zhí)行上述判斷過程，直至所述待分析導(dǎo)數(shù)對應(yīng)的主頻D^j不滿足所述循環(huán)終止條件，將不滿足所述循環(huán)終止條件的待分析導(dǎo)數(shù)對應(yīng)的j確定為所述疊加閾值。

7.根據(jù)權(quán)利要求6所述的裝置，其特征在于，

所述預(yù)設(shè)的終止參數(shù)ε的取值范圍為0.2至0.3。

8.根據(jù)權(quán)利要求5至7任一所述的裝置，其特征在于，

所述Nor( )表示對括號內(nèi)的內(nèi)容進行歸一化處理后得到的數(shù)值在區(qū)域[-1，1]內(nèi)。