[发明专利]一种块状介质的宽频反演方法

权利要求书

1.一种块状介质的宽频反演方法，其特征在于，包括以下步骤：

步骤1：分析地震频带范围，根据地震数据的频带提取实际子波用于压缩感知的拓频；

步骤2：利用压缩感知理论建立地震的褶积模型，利用压缩感知理论进行拓频获得地震信号中全频带的反射系数，地震数据表示为反地震子波和地下反射系数的褶积，利用快速阈值迭代方法计算地震全频带反射系数；

步骤3：地震道积分是通过地震反射系数的积分求和，实现地震数据转换为地震相对阻抗剖面，采用多重积分技术对拓低频的反射系数进行奇数次多重积分，获得不同主频的相对阻抗体；

步骤4：扫描奇数次多重积分的主频，同时扫描测井阻抗频谱在积分剖面主频下的振幅相对关系；

步骤5：以测井阻抗的不同主频下的相对振幅关系为准则，计算奇数次多重积分的加权系数，利用加权相加方式重构获得地震相对阻抗反演体。

2.根据权利要求1所述的块状介质的宽频反演方法，其特征在于：所述步骤2中，地震数据用s表示，反地震子波用w表示，地下反射系数用r表示，地震数据s表示为反地震子波w和地下反射系数r的褶积s＝w*r+n，频率域表示为S＝WR＝WFr，F是Fourier变换矩阵，利用压缩感知理论，通过非线性重构来重建信号，令A＝WF,y＝S,x＝r，获得成本函数如下式：

并使成本函数的值最小，利用快速阈值迭代方法计算地震反射系数，其中||.||₂和||.||₁分别为L₂范数和L₁范数，是L₂范数和L₁范数权重的调节因子，即越大，L₁范数占的权重就越大。

3.根据权利要求2所述的块状介质的宽频反演方法，其特征在于：所述步骤3中，地震道积分的公式如下：

其中s(t)为地震反射系数，m为积分重数，s_m(t)为积分结果，在上式中，N为数据总样点数，n为当前位置计算时间样点数，Δt为采样间隔，k为正整数。

4.根据权利要求1所述的块状介质的宽频反演方法，其特征在于：所述奇数次多重积分为一重积分、三重积分、五重积分、七重积分。