[发明专利]一种基于纵波和转换横波联合反演截距和梯度的方法

说明书

技术领域

本发明属于地球物理勘探领域，具体地，涉及一种基于纵波和转换横波联合反演截距和梯度的方法。

背景技术

振幅随偏移距变化(Amplitude Versus Offet)是通过研究反射振幅随偏移距变化关系来获取速度、密度等岩性参数，进而进行油气储层预测的一项技术。利用Shuey(1985)公式提取相应的截距(P)和梯度(G)或者它们的组合参数可以有效的进行油气识别和流体预测。其中，截距(P)剖面，即直垂入射反射系数剖面，具有较高的信噪比和分辨率，在近似零炮检距处垂直入射、反射，不会产生转换横波，可视为纵波剖面。梯度(G)剖面与泊松比有关，含有反射界面上下地层岩性变化的信息，截距(P)和梯度(G)是很好的油气指示和流体预测参数。

但是，传统提取截距(P)和梯度(G)的方法仅仅依赖于纵波资料，难以克服提取过程中的不稳定性和提取结果的多解性；另外一个方面，多分量地震资料记录到了丰富的地震波场信息，可以弥补单纯利用纵波方法的不足，减少反演过程的多解性，增加反演结果的稳定性，为油气预测提供可靠的方法。本发明综合利用纵波和转换横波联合提取截距(P)和梯度(G)，充分利用接收到的多分量地震波场信息，减少了反演过程的多解性，提高了反演过程的稳定性，为油气预测提供了一种可靠的方法。

发明内容

为了克服单独使用纵波资料反演截距(P)和梯度(G)参数的不足，本发明提供了一种基于纵波和转换横波联合反演截距(P)和梯度(G)的方法；充分利用纵波、转换横波波场信息，在以截距(P)和梯度(G)为参数的纵波反射系数近似式(Shuey，1985)基础上，进一步推导出了以截距(P)和梯度(G)为参数的转换横波反射系数近似式，利用获得的纵波和转换横波角道集，采用奇异值分解方法进行联合反演，获得截距(P)和梯度(G)。

为实现上述目的，本发明的技术方案如下：

基于纵波和转换横波联合反演截距和梯度的方法，包括以下步骤：

步骤1：获得纵波、转换横波角度域道集，同时获得纵波和横波速度和测井解释结果

步骤2：对纵波和转换横波角度域道集提取子波获得角度域反射系数剖面，并利用获得的纵波和横波速度求取横纵波速度比剖面

步骤3：利用推导出的以截距(P)和梯度(G)为参数的转换横波反射系数公式以及纵波反射系数近似式，采用奇异值分解方法进行联合反演，获得截距(P)和梯度(G)

步骤4：利用测井解释结果对反演出的截距(P)和梯度(G)进行标定建立流体识别模式，对未钻井区进行流体检测

相对于现有技术，本发明具有如下的有益效果：本发明综合利用纵波和转换横波联合提取截距(P)和梯度(G)，充分利用接收到的多分量地震波场信息，减少了反演过程的多解性，提高了反演过程的稳定性，为油气预测提供了一种可靠的方法。

附图说明

图1是纵波和转换横波联合提取截距(P)和梯度(G)方法流程图；

图2是进行纵波和转换横波联合提取截距(P)和梯度(G)的层状模型数据；

图3a是由Zoeppritz方程直接合成的纵波角度域道集(0-38°)；

图3b是由Zoeppritz方程直接合成的转换横波角度域道集(0-38°)；

图4a是用最小平方脉冲反褶积原理去掉图3a中纵波角度域道集的子波得到的纵波角度域反射系数；

图4b是用最小平方脉冲反褶积原理去掉图3b中转换横波角度域道集的子波得到的转换横波角度域反射系数；

图5a是纵波和转换横波联合提取的截距(P)和仅利用纵波提取的截距(P)与截距(P)真实值的比较；

图5b是利用纵波和转换横波联合提取的梯度(G)和仅利用纵波提取的梯度(G)与梯度(G)真实值的比较；

图6a是由Zoeppritz方程直接合成的纵波角度域道集(0-38°)并加有25dB背景噪声；

图6b是由Zoeppritz方程直接合成的转换横波角度域道集(0-38°)并加有25dB背景噪声；

图7a是用最小平方脉冲反褶积原理去掉图6a中纵波角度域道集的子波得到的纵波角度域反射系数；

图7b是用最小平方脉冲反褶积原理去掉图6b中转换横波角度域道集的子波得到的转换横波角度域反射系数；

图8a在输入角道集中加入随机噪声情况下，纵波和转换横波联合提取的截距(P)和仅利用纵波提取的截距(P)与截距(P)真实值的比较；