一种地震散射P-S转换波成像方法

摘要

一种地震散射P-S转换波成像方法，包括如下步骤：第一步：将地震数据读取到二维数组F中，将观测系统参数加载到原始地震数据道头中，依据观测系统和采集参数计算出散射点位置和坐标；第二步：依据散射波时距双曲线方程，固定t0i的情况下，以成像速度来确定散射双曲线轨迹，然后计算地震散射P-S转换波正常时差；第三步：从地震波旅行时中减去地震散射P-S转换波正常时差；第四步：对各炮检距上减去地震散射P-S转换波正常时差的双曲线轨迹的散射振幅进行加权求和，便实现地震散射P-S转换波成像；第五步：对成像结果按地震数据格式输出，具有水平叠加成像技术的叠加次数大大提高，能有效提高信噪比，而且能有效改善成像精度，获得的地质信息更为丰富。

主权项

1.一种地震散射P-S转换波成像方法，其特征在于，包括如下步骤：第一步：将地震数据读取到二维数组F中，同时将观测系统参数加载到原始地震数据道头中，并依据观测系统和采集参数计算出散射点位置和坐标；第二步：依据散射波时距双曲线方程，在炮集上，固定t_0i的情况下，以成像速度来确定散射双曲线轨迹，包括纵波速度v_P和横波速度v_S，然后计算地震散射P-S转换波正常时差。地震散射P-S转换波时距双曲线方程：$t_{\mathrm{PSji}}=\frac{1}{v_P}\sqrt{L_j^2+z_{0i}^2}+\frac{1}{v_S}\sqrt{{(x_j-L_j)}^2+z_{0i}^2}$地震P波垂直旅行时间：地震S波垂直旅行时间：地震散射P-S转换波自激自收时间：地震散射P-S转换波正常时差：${\mathrm{\Delta t}}_{\mathrm{PSji}}=t_{\mathrm{PSji}}-t_{0\mathrm{PSi}}=\frac{1}{v_P}\sqrt{L_j^2+z_{0i}^2}+\frac{1}{v_S}\sqrt{{(L_j-x_j)}^2+z_{0i}^2}-z_{0i}(\frac{1}{v_P}+\frac{1}{v_S})$式中，j＝1，2，…，m为地震道号，t_PSji为地震散射P-S转换波旅行时，t_0PSi为散射点垂直地表地震散射P-S转换波双程旅行时间，x_j为炮检距，z_0i为散射点距地表视深度，v_P为地震散射P波传播速度，v_S为地震散射S波传播速度，L_j为炮散距；第三步：从地震波旅行时中减去地震散射P-S转换波正常时差；第四步：对各炮检距上减去地震散射P-S转换波正常时差的双曲线轨迹的散射振幅进行加权求和，便实现了地震散射P-S转换波成像。地震散射波振幅加权求和原理为：式中：E为正常时差校正后的地震记录，m为地震记录总道数，j为道序号(j＝1，2……m)，i为采样点序号(i＝1，2……n)；第五步：对成像结果按地震数据格式输出。