一种地震散射Ｐ－Ｓ转换波成像速度分析方法技术

一种地震散射Ｐ－Ｓ转换波成像速度分析方法，包括如下步骤：第一步：将地震散射Ｐ－Ｐ波地震数据读取到二维数组Ｆ中，计算出散射点位置和坐标；第二步：任意选择一个速度ｖＰｋ，对各炮检距上的散射振幅进行叠加或相关，对同一散射点的能量加权求和，便得到相应的平均振幅；第三步，制作地震散射Ｐ－Ｐ波速度谱；第四步：将地震散射Ｐ－Ｓ转换波地震数据读取到二维数组Ｅ中，计算出散射点位置和坐标；第五步：带入上边以获取的地震散射Ｐ波速度ｖＰ，对各炮检距上的散射振幅进行叠加或相关，对同一散射点的能量加权求和，便得到相应的平均振幅：第六步，制作地震散射Ｐ－Ｓ转换波速度谱，具有速度分析技术的叠加次数提高，使有效散射波能量聚焦能力更强的特点。

【技术实现步骤摘要】

本专利技术属于地震波速度参数提取与成像处理
，特别涉及一种地震散射P-S转换波成像速度分析方法。
技术介绍
速度是地震波传播的基础，是成像的关键因素之一，速度获取的准确与否，直接关系到地震波成像精确程度。反射转换波时距曲线对共转换点来说既不是对称的，也不是双曲线，传统双曲线近似对转换波来说已根本不适用，主要问题来自于CCP道集的分选和转换波时距关系的固有特性第一，即便是均匀介质的水平反射层，转换波的转换反射点的水平位置随着炮检距、反射层的深度以及纵横波速比的变化而变化；第二，反射转换波的时距曲线不是双曲线，若用双曲线的时距关系进行速度分析，不会得到正确的横波速度值，特别是对中-浅层及大炮检距会有较大误差；另外，共转换点CCP道集处理技术是针对水平层状均匀介质模型发展起来的一套处理技术，但还没有实现真正的共转换点道集选排和速度分析以及叠加成像，只是目的层共转换点选排处理，严格来说，该技术只适用于水平层状均匀介质，而实际中，这种理想的情况几乎找不到，因此，最终获取的速度参数和成像效果都存在很大的误差；另外，反射波在地震波成份中是非常有限的，而反射波与其它成份地震波传播的规律是不同的，把其它成份的地震波愣是用反射理论来处理，所获得的速度参数必定不精确。
技术实现思路
为了克服上述现有技术的不足，本专利技术的目的在于提出了一种地震散射P-S转换波成像速度分析方法，该方法无须进行道集选排，基于点散射地质模型，依据散射转换波时距双曲线规律，对具有更为广泛意义的地震散射波进行成像速度分析，与传统地震波速度分析技术相比较，在传统覆盖次数一定的情况下，由于充分利用上了所有采集到的地震信息进行速度分析，本专利技术比传统速度分析技术的叠加次数大大提高，这样不但使有效散射波能量聚焦能力更强，而且能有效提高信噪比，同时利用上了传统意义上的绕射波、断面波等异常波信息进行速度分析，因此，该方法能更精确地获取成像速度，尤其是对构造复杂地区和低信噪比地区的地震资料，能有效避免传统技术中异常波的伪能量团影响的弊端，同时可减小多次波的影响。 为了实现上述目的，本专利技术采用的技术方案是一种地震散射P-S转换波成像速度分析方法，包括如下步骤 第一步将地震散射P-P波地震数据读取到二维数组F中，同时将观测系统参数加载到原始地震数据道头中，并依据观测系统和采集参数计算出散射点位置和坐标； 第二步依据地震散射P-P波时距双曲线方程，在炮集上，固定t0PPi的情况下，任意选择一个速度vPk，该速度可以确定一条双曲线轨迹，并沿该双曲线轨迹对各炮检距上的散射振幅进行叠加或相关，然后再对同一散射点的能量加权求和，便得到相应的平均振幅 式中rji＝tji/Δτ，Δτ为采样率，tji为延迟时，m为地震记录总道数，j为道序号(j＝1，2……m)，i为采样点序号(i＝1，2……n)。 地震散射P-P波时距双曲线方程 地震散射P-P波正常时差为 地震散射P-P波成像速度 式中，j＝1，2，…，m为地震道号，tji为延迟时间，z0i为散射点距地表视深度，vP为地震散射P波传播速度，Lj为炮散距，xj为炮检距。 当vPk为正确的散射P波成像速度时，不同炮检距地震道上的振幅同相叠加或相关，叠加振幅最大，沿不同速度定义的双曲线轨迹计算叠加振幅就是对双曲线轨迹上的地震道进行相关性度量，通过不同速度对应的叠加振幅分析，便可提取地震散射P波传播速度； 第三步，制作地震散射P-P波速度谱。首先选定一系列双程垂直旅行时间 t0PP1，t0PP2，t0PP3，Λ，t0PPi，Λt0PPn 对于每一个双程散射时间t0PPi，再选定一系列的散射成像速度，其对应关系为一对多 vP1，vP2，vP3，ΛvPk，ΛvPM t0PPi扫描时间和正常时差校正扫描速度应该包括所有旅行时间和可能的散射波成像速度，计算每个网格点(t0PPi，vPk)上的平均振幅 以等值线的形式显示出来，就得到了用于地震散射P-P波速度分析的速度谱。 为了使速度谱中的平均振幅更加稳定、突出，实际用时一般选择以(t0PPi，vPk)所定义的散射双曲线为中心、宽度为一个子波长度为时窗，计算这个时窗内的平均能量或平均振幅，据此计算结果绘制速度谱。 常用的判别准则有如下几种 平均振幅能量准则 平均振幅准则 相似系数准则 式中，rji＝tji/Δτ，Δτ为采样率，tji为延迟时，N为信号延续时窗内的采样点数，m为地震记录总道数，j为道序号(j＝1，2……m)，i为采样点序号(i＝1，2……n)； 第四步将地震散射P-S转换波地震数据读取到二维数组E中，同时将观测系统参数加载到原始地震数据道头中，并依据观测系统和采集参数计算出散射点位置和坐标； 第五步依据地震散射P-S转换波时距双曲线方程，在炮集上，固定t0PSi的情况下，并带入上边以获取的地震散射P波速度vP，同时任意选择一个速度vSk，该速度可以确定一条双曲线轨迹，并沿该双曲线轨迹对各炮检距上的散射振幅进行叠加或相关，然后再对同一散射点的能量加权求和，便得到相应的平均振幅 式中rji＝tji/Δτ，Δτ为采样率，tji为延迟时，m为地震记录总道数，j为道序号(j＝1，2……m)，i为采样点序号(i＝1，2……n)。 地震散射P-S转换波时距双曲线方程 地震散射P波的垂直旅行时间 地震散射S波的垂直旅行时间 地震散射P-S转换波自激自收时间 地震散射P-S转换波的正常时差 式中，j＝1，2，…，m为地震道号，tPSj为P-S转换散射波旅行时，t0PSi为散射点垂直地表地震散射P-S转换波双程旅行时间，xj为炮检距，z0i为散射点距地表视深度，vP为地震散射P波传播速度，vS为地震散射S波传播速度，Lj为炮散距。 当vSk为正确的散射S波成像速度时，不同炮检距地震道上的振幅同相叠加或相关，叠加振幅最大，沿不同速度定义的双曲线轨迹计算叠加振幅就是对双曲线轨迹上的地震道进行相关性度量，通过不同速度对应的叠加振幅分析，便可提取地震波传播速度； 第六步，制作地震散射P-S转换波速度谱。首先选定一系列双程垂直旅行时间 t0PS1，t0PS2，t0PS3，Λ，t0PSi，Λt0PSn 对于每一个双程散射时间t0PSi，在选定一系列的散射成像速度，其对应关系为一对多 vS1，vS2，vS3，ΛvSk，ΛvSM t0PSi扫描时间和正常时差校正扫描速度应该包括所有旅行时间和可能的散射成像速度，计算每个网格点(t0PSi，vSk)上的平均振幅 以等值线的形式显示出来，就得到了用于速度分析的速度谱。 为了使速度谱中的平均振幅更加稳定、突出，实际用时一般选择以(t0PSi，vSk)所定义的散射双曲线为中心、宽度为一个子波长度为时窗，计算这个时窗内的平均能量或平均振幅，据此计算结果绘制速度谱。 常用的判别准则有如下几种 平均振幅能量准则 平均振幅准则 相似系数准则 式中，rji＝tPSji/Δτ，Δτ为采样本文档来自技高网...

【技术保护点】
一种地震散射Ｐ－Ｓ转换波成像速度分析方法，其特征在于，包括如下步骤：第一步：将地震散射Ｐ－Ｐ波地震数据读取到二维数组Ｆ中，同时将观测系统参数加载到原始地震数据道头中，并依据观测系统和采集参数计算出散射点位置和坐标；第二步：依据地震散射Ｐ－Ｐ波时距双曲线方程，在炮集上，固定ｔ↓［０ＰＰｉ］的情况下，任意选择一个速度ｖ↓［Ｐｋ］，该速度可以确定一条双曲线轨迹，并沿该双曲线轨迹对各炮检距上的散射振幅进行叠加或相关，然后再对同一散射点的能量加权求和，便得到相应的平均振幅：＊（ｉ）＝１为采样率，ｔｊｉ为延迟时，Ｎ为信号延续时窗内的采样点数，ｍ为地震记录总道数，ｊ为道序号（ｊ＝１，２……ｍ），ｉ为采样点序号（ｉ＝１，２……ｎ）；第四步：将地震散射Ｐ－Ｓ转换波地震数据读取到二维数组Ｅ中，同时将观测系统参数加载到原始地震数据道头中，并依据观测系统和采集参数计算出散射点位置和坐标；第五步：依据地震散射Ｐ－Ｓ转换波时距双曲线方程，在炮集上，固定ｔ↓［０ＰＳｉ］的情况下，并带入上边以获取的地震散射Ｐ波速度ｖ↓［Ｐ］，同时任意选择一个速度ｖ↓［Ｓｋ］，该速度可以确定一条双曲线轨迹，并沿该双曲线轨迹对各炮检距上的散射振幅进行叠加或相关，然后再对同一散射点的能量加权求和，便得到相应的平均振幅：＊（ｉ）＝１／ｍ＊Ｅ↓［ｊ，ｉ＋ｒ↓［ｊｉ］］式中：ｒ↓［ｊｉ］＝ｔ↓［ｊｉ］／Δτ，Δτ为采样率，ｔ↓［ｊｉ］为延迟时，ｍ为地震记录总道数，ｊ为道序号（ｊ＝１，２……ｍ），ｉ为采样点序号（ｉ＝１，２……ｎ）。地震散射Ｐ－Ｓ转换波时距双曲线方程：ｔ↓［ＰＳｊｉ］＝１／ｖ↓［Ｐ］＊＊＊＋１／ｖ↓［Ｓ］＊＊＊地震散射Ｐ波的垂直旅行时间：ｔ↓［０Ｐｉ］＝ｚ↓［０ｉ］／ｖ↓［Ｐ］地震散射Ｓ波的垂直旅行时间：ｔ↓［０Ｓｉ］＝ｚ↓［０ｉ］／ｖ↓［Ｓ］地震散射Ｐ－Ｓ转换波自激自收时间：ｔ↓［０ＰＳｉ］＝ｔ↓［０Ｐｉ］＋ｔ↓［０Ｓｉ］＝ｚ↓［０ｉ］／ｖ↓［Ｐ］＋ｚ↓［０ｉ］／ｖ↓［Ｓ］地震散射Ｐ－Ｓ转换波的正常时差：Δｔ↓［ＰＳｊｉ］＝ｔ↓［ＰＳｊ］－ｔ↓［０ＰＳｉ］＝１／ｖ↓［Ｐ］＊＊＊＋１／ｖ↓［Ｓ］＊＊＊－ｚ↓［０ｉ］（１／ｖ↓［Ｐ］＋１／ｖ↓［Ｓ］）式中，ｊ＝１，２，…，ｍ为地震道号，ｔ↓［ＰＳｊ］为Ｐ－Ｓ转换散射波旅行时，ｔ↓［０ＰＳｉ］为散射点垂直地表地震散射Ｐ－Ｓ转换波双程旅行时间，ｘ↓［ｊ］为炮检距，ｚ↓［０ｉ］为散射点距地表视深度，ｖ↓［Ｐ］为地震...

【技术特征摘要】

【专利技术属性】
技术研发人员：沈鸿雁，
申请(专利权)人：西安石油大学，
类型：发明
国别省市：87[中国|西安]