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

摘要

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

主权项

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