[发明专利]一种地震散射P-P波成像速度分析方法

说明书

技术领域

本发明属于地震波速度参数提取与成像处理技术领域，特别是一种地震散射P-P波成像速度分析方法。 

技术背景

速度是地震波传播的基础，是成像的关键因素之一，速度获取的准确与否，直接关系到地震波成像精确程度。传统地震波成像技术的速度获取出发点均源于反射波时距方程速度分析技术，这是传统地震波成像技术弊端的根源。反射波在地震波成份中是非常有限的，而反射波与其它成份地震波传播的规律是不同的，把其它成份的地震波愣是用反射理论来处理，所获得的速度参数必定不精确；另外，共中心点(CMP)道集处理技术是针对水平层状均匀介质模型发展起来的一套处理技术，严格来说，该技术只适用于水平层状均匀介质，而实际中，这种理想的情况几乎找不到。随后，在反射理论的基础上，发展起来了叠后偏移成像技术和叠前部分偏移+叠后偏移技术，成像效果比常规水平叠加技术大大改观，但还是难以跳出基于反射理论处理非反射波成份存在弊端的怪圈；紧接着发展了叠前偏移成像技术，在某种程度上是以增大运算成本为代价的一种成像技术，速度模型还是该成像技术的关键，虽然用修正速度模型的手段弱化了基于CMP道集处理获取初始速度模型的弊端，但如何修正速度模型才能正确成像的问题又摆在眼前，因为地下未知的地质情况正是想要通过处理地震数 据所想获取的，也就是速度模型预先不知道，速度模型修正准确与否直接关系到地震波成像的效果，所以叠前偏移成像技术还是跳进了一个如何修正速度模型与地震波精确成像矛盾的怪圈。 

发明内容

为了克服上述现有技术的不足，本发明的目的在于提出了一种地震散射P-P波成像速度分析方法，与传统地震波速度分析技术相比较，在传统覆盖次数一定的情况下，由于充分利用上了所有采集到的地震信息进行速度分析，本发明比传统速度分析技术的叠加次数大大提高，这样不但使有效散射波能量聚焦能力更强，而且能有效提高信噪比，同时利用上了传统意义上的绕射波、断面波等异常波信息进行速度分析，因此，该方法能更精确地获取成像速度，尤其是对构造复杂地区和低信噪比地区的地震资料，能有效避免传统技术中异常波的伪能量团影响的弊端，同时可减小多次波的影响。 

为了实现上述目的，本发明采用的技术方案是：一种地震散射P-P波成像速度分析方法，包括如下步骤： 

第一步：将地震数据读取到二维数组F中，同时将观测系统参数加载到原始地震数据道头中，并依据观测系统和采集参数计算出散射点位置和坐标； 

第二步：依据散射波时距双曲线方程，在炮集上，固定t_0i的情况下，任意选择一个速度v_k，该速度可以确定一条双曲线轨迹，并沿该双曲线轨迹对各炮检距上的散射振幅进行叠加或相关，然后再对同一散射点的能量加权求和，便得到相应的平均振幅 

F‾(i)=1mΣj=1mFj,i+rj]]>

式中：r_j＝t_j/Δτ，Δτ为采样率，t_j为延迟时，m为地震记录总道数，j为道序号(j＝1，2……m)，i为采样点序号(j＝1，2……n)。 

依据下列公式，提取地震散射P-P波传播速度。 

散射P-P波时距双曲线方程： 

散射P-P波正常时差： 

散射P-P波成像速度v_P： 

式中，j＝1，2，…，m为地震道号，t_ji为延迟时间，z_0j为散射点距地表视深度，v_P为地震波传播速度，L_j为炮散距，X_j为炮检距。