地震数据处理方法和装置制造方法及图纸

本发明专利技术公开了一种地震数据处理方法和装置，属于地质勘测领域。所述方法包括：对动校正处理后的S个地震道集中预设范围内的地震道集进行叠加得到模型道，S为整数；计算各个地震道集与模型道的相关系数，并选出相关系数最大的第K地震道集；在预设时间范围内滑动时窗，计算第K‑1地震道集与第K地震道集的相关系数，并选取相关系数最大时的时窗作为最优点，将最优点对应的时移量作为第K‑1地震道集的剩余动校正量；在预设时间范围内滑动时窗，计算第M地震道集与第M+1地震道集的相关系数，并选取相关系数最大时的时窗作为最优点，将最优点对应的时移量作为第M地震道集的剩余动校正量，K‑2≥M≥1；采用同样方式计算第K+1至第S地震道集的剩余动校正量。

Seismic data processing method and device

The invention discloses a method and a device for processing seismic data, which belongs to the field of geological survey. The method includes: centralized preset seismic range set of S superposition model dynamic correction after the earthquake, S is an integer; calculated the correlation coefficient of each seismic trace set and model, and select the K seismic correlation coefficient of the largest set; at a preset time range when sliding the window, calculate the K 1 sets of correlation coefficients and the K seismic gathers the seismic time window, and select the maximum correlation coefficient as the most advantages will most advantages corresponding to the shift as the remaining K 1 seismic gathers the amount of correction; at a preset time within the sliding window. The correlation coefficient calculation section M seismic gathers and M+1 seismic data, the time window and select the maximum correlation coefficient as the most advantages, will has the corresponding residual time shift as the M seismic gathers correction, K 2 = M = 1; by the same way. Residual dynamic correction of seismic channel set from K+1 to S.

【技术实现步骤摘要】

本专利技术涉及地质勘测领域，特别涉及一种地震数据处理方法和装置。
技术介绍
地震勘探技术是利用人工激发的地震波在地层内传播来勘探地下的地质情况。在地面某处激发的地震波向地下传播时，遇到不同弹性的地层分界面就会产生反射波或折射波返回地面，用专门的仪器可记录这些波，分析所记录的波的特点(如波的传播时间、振动形状等)，能较准确地测定这些界面的深度和形态，判断地层的岩性。由于地震记录是有噪音的，因此通常会通过不同位置进行激发、接收地下同一点的地震反射波，将这些来自同一点的地震反射波的集合称为地震道集。通过多次地震信号叠加，可以达到消除噪音，提高信噪比的目的。由于野外采集的地震资料统一反射点的地震反射时间是不一样的，其反射时间与偏移距、地震传播速度有关，因此，需要将不同偏移距地震记录时间校正到自激自收地震反射时间位置，这个过程叫动校正，其目的是为了便于对地震信号叠加。但是，因为地震波的传播速度是未知的，因此动校正的实现过程通常都设计一系列地震传播速度对地震道集进行动校正处理。由于野外地震采集数据量非常庞大，人们不可能对每一个地震反射点都进行速度分析。往往都是间隔几十道才对一个地震反射点进行速度分析，而中间未进行速度分析的道集采用已作速度分析道集的速度进行插值而得，由于地震速度的非均质性，这些插值的速度或多或少都会存在一定的误差，动校正结果往往很难达到理想效果，很难将所有地震道集同相轴都校正到水平状态。
技术实现思路
为了解决现有技术的问题，本专利技术实施例提供了一种地震数据处理方法和装置。所述技术方案如下：一方面，本专利技术实施例提供了一种地震数据处理方法，所述方法包括：对动校正处理后的S个地震道集中预设范围内的地震道集进行叠加得到模型道，S为整数；计算各个地震道集与所述模型道的相关系数，并选出相关系数最大的第K地震道集；在预设时间范围内滑动时窗，计算第K-1地震道集与所述第K地震道集的相关系数，并选取相关系数最大时的时窗作为最优点，将最优点对应的时移量作为第K-1地震道集的剩余动校正量，K＞1且K为整数；在所述预设时间范围内滑动时窗，计算第M地震道集与第M+1地震道集的相关系数，并选取相关系数最大时的时窗作为最优点，将最优点对应的时移量作为第M地震道集的剩余动校正量，K-2≥M≥1且M为整数；在所述预设时间范围内滑动时窗，计算第K+1地震道集与所述第K地震道集的相关系数，并选取相关系数最大时的时窗作为最优点，将最优点对应的时移量作为第K+1地震道集的剩余动校正量；在所述预设时间范围内滑动时窗，计算第N+1地震道集与第N地震道集的相关系数，并选取相关系数最大时的时窗作为最优点，将最优点对应的时移量作为第N+1地震道集的剩余动校正量，S-1≥N≥K+1且N为整数，所述最优点对应的时移量为时窗内的一个地震道集的地震反射时间与偏移距为0的地震道集的地震反射时间的差值；采用计算出的所述S个地震道集的剩余动校正量对所述S个地震道集进行剩余动校正。在本专利技术实施例的一种实现方式中，所述对动校正处理后的S个地震道集中预设范围内的地震道集进行叠加得到模型道，包括：采用如下公式对多个地震道集进行叠加得到模型道： A j = Σ i = m 1 m 2 a i j ; ]]>其中，m1、m2为地震道集序号，aij为第i道的第j采样点的振幅值，m1、m2、i和j均为正整数，1≤m1＜m2≤S，i∈[m1，m2]。在本专利技术实施例的另一种实现方式中，所述计算各个地震道集与所述模型道的相关系数，包括：采用以下公式计算模型道与地震道集的相关系数： r x ‾ = Σ j = 1 n ( a j - A a ) ( b j - B a ) Σ j = 1 n ( a j - A a ) 2 Σ j = 1 n ( b j - B a ) 2 ; ]]>其中，为第x道与模型道的相关系数，aj为模型道第j个采样点的振幅值，Aa为模型道振幅平均值，bj为地震道集第j个采样点的振幅值，Ba为地震道集振幅平均值，x和j均为正整数，且1≤x≤S。在本专利技术实施例的另一种实现方式中，所述方法还包括：在对动校正处理后的S个地震道集中预设范围内的地震道集进行叠加得到模型道之前，对所述各个地震道集进行振幅能量均衡化处理。在本专利技术实施例的另一种实现方式中，所述对所述各个地震道集进行振幅能量均衡化处理，包括：在设定时窗范围内统计振幅能量：Q为统计道数，P为采样点数，aj为振幅值；计算每一地震道集的振幅能量E1： E 1 = 本文档来自技高网...

【技术保护点】
一种地震数据处理方法，其特征在于，所述方法包括：对动校正处理后的S个地震道集中预设范围内的地震道集进行叠加得到模型道，S为整数；计算各个地震道集与所述模型道的相关系数，并选出相关系数最大的第K地震道集；在预设时间范围内滑动时窗，计算第K‑1地震道集与所述第K地震道集的相关系数，并选取相关系数最大时的时窗作为最优点，将最优点对应的时移量作为第K‑1地震道集的剩余动校正量，K＞1且K为整数；在所述预设时间范围内滑动时窗，计算第M地震道集与第M+1地震道集的相关系数，并选取相关系数最大时的时窗作为最优点，将最优点对应的时移量作为第M地震道集的剩余动校正量，K‑2≥M≥1且M为整数；在所述预设时间范围内滑动时窗，计算第K+1地震道集与所述第K地震道集的相关系数，并选取相关系数最大时的时窗作为最优点，将最优点对应的时移量作为第K+1地震道集的剩余动校正量；在所述预设时间范围内滑动时窗，计算第N+1地震道集与第N地震道集的相关系数，并选取相关系数最大时的时窗作为最优点，将最优点对应的时移量作为第N+1地震道集的剩余动校正量，S‑1≥N≥K+1且N为整数，所述最优点对应的时移量为时窗内的一个地震道集的地震反射时间与偏移距为0的地震道集的地震反射时间的差值；采用计算出的所述S个地震道集的剩余动校正量对所述S个地震道集进行剩余动校正。...

【技术特征摘要】
1.一种地震数据处理方法，其特征在于，所述方法包括：对动校正处理后的S个地震道集中预设范围内的地震道集进行叠加得到模型道，S为整数；计算各个地震道集与所述模型道的相关系数，并选出相关系数最大的第K地震道集；在预设时间范围内滑动时窗，计算第K-1地震道集与所述第K地震道集的相关系数，并选取相关系数最大时的时窗作为最优点，将最优点对应的时移量作为第K-1地震道集的剩余动校正量，K＞1且K为整数；在所述预设时间范围内滑动时窗，计算第M地震道集与第M+1地震道集的相关系数，并选取相关系数最大时的时窗作为最优点，将最优点对应的时移量作为第M地震道集的剩余动校正量，K-2≥M≥1且M为整数；在所述预设时间范围内滑动时窗，计算第K+1地震道集与所述第K地震道集的相关系数，并选取相关系数最大时的时窗作为最优点，将最优点对应的时移量作为第K+1地震道集的剩余动校正量；在所述预设时间范围内滑动时窗，计算第N+1地震道集与第N地震道集的相关系数，并选取相关系数最大时的时窗作为最优点，将最优点对应的时移量作为第N+1地震道集的剩余动校正量，S-1≥N≥K+1且N为整数，所述最优点对应的时移量为时窗内的一个地震道集的地震反射时间与偏移距为0的地震道集的地震反射时间的差值；采用计算出的所述S个地震道集的剩余动校正量对所述S个地震道集进行剩余动校正。2.根据权利要求1所述的方法，其特征在于，所述对动校正处理后的S个地震道集中预设范围内的地震道集进行叠加得到模型道，包括：采用如下公式对多个地震道集进行叠加得到模型道： A j = Σ i = m 1 m 2 a i j ; ]]>其中，m1、m2为地震道集序号，aij为第i道的第j采样点的振幅值，m1、m2、i和j均为正整数，1≤m1＜m2≤S，i∈[m1，m2]。3.根据权利要求1所述的方法，其特征在于，所述计算各个地震道集与所述模型道的相关系数，包括：采用以下公式计算模型道与地震道集的相关系数： r x ‾ = Σ j = 1 n ( a j - A a ) ( b j - B a ) Σ j = 1 n ( a j - A a ) 2 Σ j = 1 n ( b j - B a ) 2 ; ]]>其中，为第x道与模型道的相关系数，aj为模型道第j个采样点的振幅值，Aa为模型道振幅平均值，bj为地震道集第j个采样点的振幅值，Ba为地震道集振幅平均值，x和j均为正整数，且1≤x≤S。4.根据权利要求1至3任一项所述的方法，其特征在于，所述方法还包括：在对动校正处理后的S个地震道集中预设范围内的地震道集进行叠加得到模型道之前，对所述各个地震道集进行振幅能量均衡化处理。5.根据权利要求4所述的方法，其特征在于，所述对所述各个地震道集进行振幅能量均衡化处理，包括：在设定时窗范围内统计振幅能量：Q为统计道数，P为采样点数，aj为振幅值；计算每一地震道集的振幅能量E1： E 1 = 1 P Σ j = 1 P a j 2 ; ]]>计算每一地震道集的振幅均衡系数：k＝E0/E1；采用所述振幅均衡系数对每一地震道集进行振幅均衡处理：其中为每一地震道集振幅值a组成的序列，为振...

【专利技术属性】
技术研发人员：张光荣，杨跃明，肖富森，张世荣，喻颐，冉崎，廖奇，周肖，马波，谢冰，梁翰，赖强，陈康，张旋，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：北京;11