[发明专利]一种地震地层厚度变化分析方法

权利要求书

1.一种地震地层厚度变化分析方法，其特征在于，该方法包括以下步骤：

对每一道地震信号重复步骤01-步骤04：

步骤01：对每一道地震信号进行经验模式分解过程，分解为n级内蕴模式函数分量，n＝1，2，3，…N，N为总分解级数；

对每一级内蕴模式函数重复步骤02-步骤04：

步骤02：对内蕴模式函数，计算它的连续小波变换；

步骤03：对内蕴模式函数构建对应的解析信号；

步骤04：计算内蕴模式函数对应的瞬时频率；

步骤05：根据得到的各级内蕴模式函数的瞬时频率变化，分析地震的地层厚度变化值。

2.如权利要求1所述一种地震地层厚度变化分析方法，其特征在于，所述步骤01中的经验模式分解过程包括如下步骤：

首先将原始地震信号作为输入信号，进行如下筛选过程；

(A).基于输入信号的局部极大值和极小值点，采用三次B样条函数插值拟合得到平均值；

(B).从输入信号中减去平均值，得到一个新的信号；

(C).步骤(B)得到的新信号如果满足给定的条件，则可以将其视作一级内蕴模式函数，并执行步骤(D)，否则，将步骤(B)得到的新信号作为输入信号，回到步骤(A)；

D.提取内蕴模式函数分量，将原始地震信号减去已得各级内蕴模式函数的信号记为残余信号，若其呈现单调趋势或者值很小时，经验模式分解过程结束，否则将其视作为下一次筛选过程的输入信号，回到步骤(A)；

所述步骤(A)中，对输入信号的局部极大值、极小值点采用三次B样条函数插值拟合平均值时，采用四点镜像取反的方法处理端点处的平均值拟合，以输入信号左右两个时间端点为镜像点，对输入信号及其局部极大值、极小值点在左右两端分别进行镜像延拓，并对延拓的信号和局部极大值、极小值点的值取反，然后再将原输入信号的局部极大值、极小值点序列左右各增加四个局部极值点，做为新的局部极大值、极小值点序列，利用该序列拟合平均值；

所述步骤(C)中，新信号满足的给定条件是指新信号和对应输入信号的标准方差小于0.05。

3.如权利要求1所述一种地震地层厚度变化分析方法，其特征在于，所述步骤02中，内蕴模式函数的连续小波变换由以下关系式决定：

Sn(b,a)=1a∫T1T2sn(t)g‾(t-ba)dt,]]>

式中，s_n(t)表示第n级内蕴模式函数，n＝1，2，3，…N，t为时间，t属于实数集合，最小值为T₁，最大值为T₂，尺度因子a＞0，b为时移，b属于实数集合，b的最小值为T₁，最大值为T₂；S_n(b，a)为s_n(t)的连续小波变换，为小波函数的复共轭，g(t-ba)=eiω0t-bae-(t-b)22a2]]>是对基本小波函数g(t)=eiω0te-t22]]>经过伸缩a和时移b得到的，基本频率ω₀＝6.0，i=-1]]>为虚数单位。