[发明专利]一种自相控叠前反演方法

权利要求书

1.一种自相控叠前反演方法，其包括以下步骤：

1)从叠前地震道集中提取全井子波，从测井信息中提取纵横波速度、密度信息，利用全井子波和纵横波速度、密度信息正演共炮点道集得到井上正演道集；

2)对叠前地震道集进行分角度抽取并叠加得到分角度井旁叠加道集，对井上正演道集进行分角度抽取并叠加得到分角度正演叠加道集；

3)采用时频变换的方法从分角度井旁叠加道集和分角度正演叠加道集中提取保幅拓频因子，利用得到的保幅拓频因子对分角度井旁叠加道集进行保幅拓频处理，得到拓频保幅记录；

4)根据地震信息、地质先验信息、测井信息、叠加速度信息构建初始低频模型，并采用拓频保幅记录对初始低频模型进行修正得到自相控低频模型；

5)基于自相控低频模型再次对拓频保幅记录进行自相控叠前AVO/AVA反演，得到最终的自相控叠前AVO/AVA反演结果，也即地层弹性参数；

6)根据最终的自相控叠前AVO/AVA反演结果，对储层进行精细标定和预测。

2.如权利要求1所述的一种自相控叠前反演方法，其特征在于：所述步骤3)中，对分角度井旁叠加道集进行保幅拓频处理的具体步骤为：

①对分角度正演叠加道集z_i(t)和分角度井旁叠加道集s_i(t)分别作时频变换，得到时频谱Z_i(t,f)和S_i(t,f)；

②将时频谱S_i(t,f)的频带拓宽至Z_i(t,f)的频带，得到保幅拓频因子Φ_i(t,f)：

Φi(t,f)=B[Zi(t,f)]B[Si(t,f)]+δ,]]>

其中，B[·]为包络求取算子，B[S_i(t,f)]、B[Z_i(t,f)]分别为时频谱S_i(t,f)、Z_i(t,f)的Hilbert瞬时包络，δ为白噪因子；

③将保幅拓频因子Φ_i(t,f)与时频谱S_i(t,f)相乘，得到利用保幅拓频因子拓宽后的时频谱S′_i(t,f)，对拓宽后的时频谱S′_i(t,f)进行时频反变换，得到拓频保幅处理后的时间域信号s'_i(t)；

④将不同偏移距的保幅拓频因子分别作用于分角度井旁叠加道集s_i(t)中对应的道集，完成对分角度井旁叠加道集s_i(t)的保幅拓频处理，得到拓频保幅记录。

3.如权利要求1或2所述的一种自相控叠前反演方法，其特征在于：所述步骤4)中，所述自相控低频模型的构建包括以下步骤：

①根据地质先验信息、测井信息和地震信息，采用分级控制、旋回对比的方式构建精细地层格架；

②根据测井信息，借助步骤①中的精细地层格架，构建压实背景下的低频趋势；

③根据叠加速度信息，修正叠加速度异常，并进行光滑处理，得到超低频趋势；

④融合步骤②和步骤③中得到的低频趋势和超低频趋势构建初始低频模型；

⑤基于初始低频模型对拓频保幅记录进行初次叠前AVO/AVA反演，得到初始叠前AVO/AVA反演结果；

⑥根据初始叠前AVO/AVA反演结果，采用三维砂体雕刻方法得到砂体展布趋势，将该趋势与初始低频模型进行融合，得到自相控低频模型。