[发明专利]一种快速转换波直接模拟以确定地层岩性、岩相变化的方法

摘要

本发明是一种实现转换波叠后地震记录的快速直接模拟以确定地层岩性、岩相变化的方法，步骤是利用纵、横波速度得到P－SV转换波的双程垂直旅行时，将深度域测井资料转换到P－SV波的时间域；对时间域非规则采样的测井资料进行均匀重采样，确定每个新采样点上的P－SV转换波反射系数；按常规纵波处理方法，将得到的P－SV转换波反射系数与地震子波进行褶积得到P－SV转换波的模拟地震记录；按常规纵波处理方法，将模拟结果与实际叠后地震记录进行相关处理，实现P－SV转换波层位标定，确定地层岩性、岩相变化。本发明不需要射线追踪计算入射角和透射角和对测井曲线进行滤波处理，提高了叠后地震记录模拟的计算效率，可以反映测井资料的高频响应特征。

主权项

1、一种快速转换波直接模拟以确定地层岩性、岩相变化的方法，其特征在于：通过以下具体技术步骤实现：1)用全波列测井方法记录纵、横波时差和密度，将纵、横波时差转换为纵、横波速度；2)利用纵、横波速度按下式得到P-SV转换波的双程垂直旅行时，将深度域测井资料转换到P-SV波的时间域；${\mathrm{Tps}}_i=\mathrm{\Delta h}\underset{i=1}{\overset{n}{\Sigma }}(1/{\alpha }_i+1/{\beta }_i)---\left(3\right)式中，Tpsi为深度域测井资料转换到P-SV转换波时间域的双程旅行时，\Delta h为深度域测井资料的采样间隔，\alpha i和\beta i分别为深度域测井曲线上第i个采样点上的纵、横波速度值，n为深度域测井资料的采样数；3)对时间域非规则采样的测井资料进行均匀重采样，按下式确定每个新采样点上的P-SV转换波反射系数；$ {\mathrm{Rps}}_i=2\frac{{\beta }_i+{\beta }_{i+1}}{{\alpha }_i+{\alpha }_{i+1}}\frac{{\rho }_i{\beta }_i-{\rho }_{i+1}{\beta }_{i+1}}{{\rho }_{i+1}{\beta }_{i+1}+{\rho }_i{\beta }_i}\sin 2\theta ---\left(4\right)式中，\theta 为P-SV转换波叠加剖面的等效角度，\alpha i、\beta i和\rho i分别为时间域测井曲线上第i个采样点上的纵波速度、横波速度和密度，\alpha i+1、\beta i+1和\rho i+1分别为时间域测井曲线上第i+1个采样点上的纵波速度、横波速度和密度；4)按常规纵波处理方法，将步骤3)中得到的P-SV转换波反射系数与地震子波进行褶积得到P-SV转换波的模拟地震记录；5)按常规纵波处理方法，将模拟结果与实际叠后地震记录进行相关处理，实现P-SV转换波层位标定，确定地层岩性、岩相变化。$$