一种地震散射P-S转换波成像速度分析方法
专利汇可以提供一种地震散射P-S转换波成像速度分析方法专利检索，专利查询，专利分析的服务。并且一种 地震 散射P-S转换波成像速度分析方法，包括如下步骤：第一步：将地震散射P-P波地震数据读取到二维数组F中，计算出散射点 位置 和坐标；第二步：任意选择一个速度vPk，对各炮检距上的散射振幅进行 叠加 或相关，对同一散射点的 能量 加权求和，便得到相应的平均振幅；第三步，制作地震散射P-P 波速 度谱；第四步：将地震散射P-S转换波地震数据读取到二维数组E中，计算出散射点位置和坐标；第五步：带入上边以获取的地震散射P波速度vP，对各炮检距上的散射振幅进行叠加或相关，对同一散射点的能量加权求和，便得到相应的平均振幅：第六步，制作地震散射P-S转换波速度谱，具有速度分析技术的叠加次数提高，使有效散射波能量聚焦能 力 更强的特点。，下面是一种地震散射P-S转换波成像速度分析方法专利的具体信息内容。
1.一种地震散射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/Δτ，Δτ为采样率，tji为延迟时，N为信号延续时窗内的采样点数，m为地震记录总道数，j为道序号(j＝1，2……m)，i为采样点序号(i＝1，2……n)。
2.根据权利要求1所述的一种地震散射P-S转换波成像速度分析方法，其特征在于，包括如下步骤：
第一步：将含有146炮，每炮86道，每道2000个采样点的P-P波地震数据读取到二维数组F中，同时将观测系统参数加载到原始地震数据道头中，并依据观测系统和采集参数计算出散射点位置和坐标；
第二步：依据地震散射P-P波时距双曲线方程，在炮集上，固定t0PPi的情况下，任意选择一个速度vPk，该速度可以确定一条双曲线轨迹，并沿该双曲线轨迹对各炮检距上的散射振幅进行叠加或相关，然后再对同一散射点的能量加权求和，便得到相应的平均振幅；
第三步，制作速度谱。首先选定一系列双程垂直旅行时间：0～1000ms，步长为2.0ms；
对于每一个双程散射时间t0PPi，再选定一系列的地震散射P-P波成像速度：800～2000m/s，步长为30m/s。
t0PPi扫描时间和正常时差校正扫描速度应该包括所有旅行时间和可能的散射成像速度，计算每个网格点(t0PPi，vPk)上的平均振幅 以等值线的形式显示出来，就得到了用于速度分析的地震散射P-P波速度谱。
为了使速度谱中的平均振幅更加稳定、突出，实际用时一般选择以(t0PPi，vPk)所定义的散射双曲线为中心、宽度为一个子波长度为时窗，该时窗长度为50ms，计算这个时窗内的平均能量或平均振幅，据此计算结果绘制速度谱。
第四步：将含有146炮，每炮86道，每道2000个采样点的P-SV转换波地震数据读取到二维数组E1中，同时将观测系统参数加载到原始地震数据道头中，并依据观测系统和采集参数计算出散射点位置和坐标；
第五步：依据地震散射P-SV转换波时距双曲线方程，在炮集上，固定t0PSi的情况下，并带入上边以获取的地震散射P波速度vP，任意选择一个速度vSVk，该速度可以确定一条双曲线轨迹，并沿该双曲线轨迹对各炮检距上的散射振幅进行叠加或相关，然后再对同一散射点的能量加权求和，便得到相应的平均振幅；
第六步，制作地震散射P-SV转换波速度谱。首先选定一系列双程垂直旅行时间：0～
1000ms，步长为2.0ms；对于每一个双程散射时间t0PSi，再选定一系列的地震散射P-SV转换波成像速度：600～1200m/s，步长为30m/s。
t0PSi扫描时间和正常时差校正扫描速度应该包括所有旅行时间和可能的散射成像速度，计算每个网格点(t0PSi，vSVk)上的平均振幅 以等值线的形式显示出来，就得到了用于速度分析的地震散射P-SV转换波速度谱。
为了使速度谱中的平均振幅更加稳定、突出，实际用时一般选择以(t0PSi，vSVk)所定义的散射双曲线为中心、宽度为一个子波长度为时窗，该时窗长度为50ms，计算这个时窗内的平均能量或平均振幅，据此计算结果绘制速度谱。
第七步：将含有146炮，每炮86道，每道2000个采样点的P-SH转换波地震数据读取到二维数组E2中，同时将观测系统参数加载到原始地震数据道头中，并依据观测系统和采集参数计算出散射点位置和坐标；
第八步：依据地震散射P-SH转换波时距双曲线方程，在炮集上，固定t0PSi的情况下，并带入上边以获取的地震散射P波速度vP，任意选择一个速度vSHk，该速度可以确定一条双曲线轨迹，并沿该双曲线轨迹对各炮检距上的散射振幅进行叠加或相关，然后再对同一散射点的能量加权求和，便得到相应的平均振幅；
第九步，制作地震散射P-SH转换波速度谱，首先选定一系列双程垂直旅行时间：0～
1000ms，步长为2.0ms；对于每一个双程散射时间t0PSi，再选定一系列的地震散射P-SH转换波成像速度：600～1200m/s，步长为30m/s。
t0PSi扫描时间和正常时差校正扫描速度应该包括所有旅行时间和可能的散射成像速度，计算每个网格点(t0PSi，vSHk)上的平均振幅 以等值线的形式显示出来，就得到了用于速度分析的地震散射P-SH转换波速度谱。
为了使速度谱中的平均振幅更加稳定、突出，实际用时一般选择以(t0PSi，vSHk)所定义的散射双曲线为中心、宽度为一个子波长度为时窗，该时窗长度为50ms，计算这个时窗内的平均能量或平均振幅，据此计算结果绘制速度谱。
一种地震散射P-S转换波成像速度分析方法
技术领域
背景技术
发明内容
2000m/s，步长为30m/s。
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 一种高温介质地震波传播模拟方法 | 2020-05-12 | 242 |
| 一种地震波初至检测方法及系统 | 2020-05-12 | 52 |
| 一种地震波初至拾取方法及系统 | 2020-05-12 | 249 |
| 一种高分辨率地震波成像方法和装置 | 2020-05-13 | 644 |
| 一种地震波S波自动拾取的算法 | 2020-05-13 | 851 |
| 地震波传播展示方法 | 2020-05-11 | 592 |
| 地震波正演方法及系统 | 2020-05-11 | 820 |
| 一种利用斜井直达地震波走时计算地层地震波速度的方法 | 2020-05-12 | 626 |
| 一种增强地震波激发照明的设计方法 | 2020-05-13 | 738 |
| 用于地震波传感器的耦合装置 | 2020-05-13 | 230 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
