【技术实现步骤摘要】
本专利技术属于地震勘探，更具体地，涉及一种地震数据处理中多层位约束计算剩余静校正的方法及介质。
技术介绍
1、在地震数据处理中，准确的静校正计算结果对地下构造成像和求取地震波速度的精度有着非常重要的作用。
2、目前存在很多计算剩余静校正的方法，包括相关法、广义线性反演法、非线性反演法等，但是这些方法都存在各种缺陷。最大能量法是地震数据处理中常用的方法，他的基本原理是，在应用正确的静校正量后，地震剖面的叠加能量达到极大值。根据这一原理，把地表一致性剩余静校正量的计算作为一个最优化问题考虑。吴波等考虑到来自同一层位相邻cmp道集的反射波存在相关性，所以提出了加权混波模型道。潘树林等在使用最大能量法计算静校正量时采用了一个计算时窗，该计算时窗沿着连续的反射同相轴追踪。即先选出地震数据的一个层位，根据层位扩展出计算时窗。该方法能够充分的利用同向轴提供的信息，计算静校正量。但是仅根据一个层位中的信息计算所有地震数据道的静校正量，会使得该层位的局部信息对整个数据的影响过大。
3、在地震数据处理中，准确的静校正计算结果对地下构造成像和求取地震波速度的精度有着非常重要的作用，因此，本领域亟需一种地震数据处理中多层位约束计算剩余静校正的方法，通过多个层位约束静校正量很有必要。
技术实现思路
1、针对现有技术存在的问题，本专利技术提供一种地震数据处理中多层位约束计算剩余静校正的方法及介质，通过拾取多个同相轴清晰的层位，由层位向上向下扩展时间窗口，计算得到受各层位信息约束的静校正
2、为实现上述目的，本专利技术提供一种地震数据处理中多层位约束计算剩余静校正量的方法，包括：
3、步骤s1：拾取多个同相轴的层位，在每个所述层位上选取时间窗口，根据预设的收敛条件，迭代计算每个层位的每个炮点和每个检波点的静校正量；
4、步骤s2：对所有层位的同一炮点和同一检波点的静校正量进行加权平均处理，得到多层位约束的每个共炮点和每个共检波点的剩余静校正量。
5、优选的，截取出所述时间窗口内的地震数据作为时窗数据用于迭代计算该层位的每个炮点和每个检波点的静校正量，包括：
6、步骤s11：基于所述时窗数据计算得出所述互相关系数，在所述时间窗口的长度范围内将所述互相关系数分别累加至该道所在多个炮点组成的炮点数组和多个检波点组成的检波点数组；
7、步骤s12：从所述炮点数组和所述检波点数组中选取最大互相关系数累加和对应的时移量分别作为炮点的静校正量和检波点的静校正量，并作为输入数据按照预设的收敛条件进行迭代处理，在迭代处理完成后，将最后一次迭代处理后得到的每个炮点的静校正量和每个检波点的静校正量作为该层最终的炮点的静校正量和检波点的静校正量。
8、优选的，所述预设的收敛条件为前后两次迭代的叠加道累加能量的差值在预设的收敛阈值内或到达预设的迭代次数。
9、优选的，对于所述时窗数据截取的所有cmp道集，所述叠加道累加能量通过下述步骤得出：
10、步骤a：对每个所述cmp道集，将该cmp道集上的所有道数据进行累加得到叠加道；
11、步骤b：将所有所述cmp道集的叠加道进行能量累加，得到所述叠加道累加能量。
12、优选的，从所述时窗数据中获取某个cmp道集的叠加道，所述步骤s11包括：
13、步骤s111：从所述叠加道中减去该cmp道集的每一道数据得到道数据集合；
14、步骤s112：从所述道数据集合中获取某一道数据的最大样点值并计算该道数据对应的标准道；
15、步骤s113：将所述cmp道集中的每一道数据与所述标准道在所述时间窗口的长度范围内计算不同时移量下的互相关系数；
16、步骤s114：将所述互相关系数分别累加至该道所在炮点数组中对应的炮号和所在检波点数组中对应的检波点号。
17、优选的，所述步骤s12中，从多个所述炮点数组和多个所述检波点数组中选取最大互相关系数累加和对应的时移量分别作为炮点的静校正量和检波点的静校正量，包括：
18、遍历完所有的所述cmp道集，得到炮点号和检波点号相同的道在相同时移量下的互相关系数的累加和，从对应的炮点数组中找出最大互相关系数累加和对应的时移量作为此次迭代计算该炮点在该层位得到的静校正量，以及从对应的检波点数组中找出最大互相关系数累加和对应的时移量作为此次迭代计算该检波点在该层位得到的静校正量。
19、优选的，所述步骤s12还包括：保存每个层位最后一次迭代计算得到的每个炮点的静校正量和该炮点的静校正量对应的互相关系数累加和，以及最后一次迭代计算得到的每个检波点的静校正量和该检波点的静校正量对应的互相关系数累加。
20、优选的，所述步骤s2包括：
21、步骤s21：将每个炮点的互相关系数累加和作为该层该炮点的能量以及每个检波点的互相关系数累加和作为该层该检波点的能量，计算所有层位的同一炮点和同一检波点的总能量；
22、步骤s22：将每层每个炮点的能量占同一炮点总能量的比重作为第一加权系数，计算出多层位约束的每个共炮点的剩余静校正量，以及，
23、将每层每个检波点的能量占同一检波点总能量的比重作为第二加权系数，计算出多层位约束的每个共检波点的剩余静校正量。
24、本专利技术还提供一种电子设备，所述电子设备包括：
25、存储器，存储有可执行指令；
26、处理器，所述处理器运行所述存储器中的所述可执行指令，以实现上述任一项所述的一种地震数据处理中多层位约束计算剩余静校正量的方法。
27、本专利技术还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述任一项所述的一种地震数据处理中多层位约束计算剩余静校正量的方法。
28、本专利技术的技术方案的有益效果在于：
29、相较于现有技术的最大能量法计算得到剩余静校正量无法实现多个地震层位的约束，本专利技术通过人工拾取多个同相轴清晰的层位，对于每个层位，首先根据层位向上向下扩展一个时间窗口，然后在在该时间窗口中迭代计算静校正量，得到各层位上每个炮点和每个检波点对应的静校正量后，根据炮点或检波点上的互相关系数累加和的能量占总能量的比重，计算各个层位静校正量的加权系数，得到最终的多层位约束的每个共炮点和每个共检波点的剩余静校正量，多个同相轴的约束使得计算的静校正量能够充分考虑多个层位的信息，能够更好地反映地下构造。
本文档来自技高网...【技术保护点】
1.一种地震数据处理中多层位约束计算剩余静校正量的方法，其特征在于，包括：
2.根据权利要求1所述的地震数据处理中多层位约束计算剩余静校正量的方法，其特征在于，所述步骤S1中，截取出所述时间窗口内的地震数据作为时窗数据用于迭代计算该层位的每个炮点和每个检波点的静校正量，包括：
3.根据权利要求2所述的地震数据处理中多层位约束计算剩余静校正量的方法，其特征在于，所述预设的收敛条件为前后两次迭代的叠加道累加能量的差值在预设的收敛阈值内或到达预设的迭代次数。
4.根据权利要求3所述的地震数据处理中多层位约束计算剩余静校正量的方法，其特征在于，对于所述时窗数据截取的所有cmp道集，所述叠加道累加能量通过下述步骤得出：
5.根据权利要求4所述的地震数据处理中多层位约束计算剩余静校正量的方法，其特征在于，从所述时窗数据中获取某个cmp道集的叠加道，所述步骤S11包括：
6.根据权利要求5所述的地震数据处理中多层位约束计算剩余静校正量的方法，其特征在于，所述步骤S12中，从多个所述炮点数组和多个所述检波点数组中选取最大互相关系数累加和
7.根据权利要求2所述的地震数据处理中多层位约束计算剩余静校正量的方法，其特征在于，所述步骤S12还包括：保存每个层位最后一次迭代计算得到的每个炮点的静校正量和该炮点的静校正量对应的互相关系数累加和，以及最后一次迭代计算得到的每个检波点的静校正量和该检波点的静校正量对应的互相关系数累加。
8.根据权利要求7所述的地震数据处理中多层位约束计算剩余静校正量的方法，其特征在于，所述步骤S2包括：
9.一种电子设备，其特征在于，所述电子设备包括：
10.一种非暂态计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1-8任一项所述的一种地震数据处理中多层位约束计算剩余静校正量的方法。
...【技术特征摘要】
1.一种地震数据处理中多层位约束计算剩余静校正量的方法，其特征在于，包括：
2.根据权利要求1所述的地震数据处理中多层位约束计算剩余静校正量的方法，其特征在于，所述步骤s1中，截取出所述时间窗口内的地震数据作为时窗数据用于迭代计算该层位的每个炮点和每个检波点的静校正量，包括：
3.根据权利要求2所述的地震数据处理中多层位约束计算剩余静校正量的方法，其特征在于，所述预设的收敛条件为前后两次迭代的叠加道累加能量的差值在预设的收敛阈值内或到达预设的迭代次数。
4.根据权利要求3所述的地震数据处理中多层位约束计算剩余静校正量的方法，其特征在于，对于所述时窗数据截取的所有cmp道集，所述叠加道累加能量通过下述步骤得出：
5.根据权利要求4所述的地震数据处理中多层位约束计算剩余静校正量的方法，其特征在于，从所述时窗数据中获取某个cmp道集的叠加道，所述步骤s11包括：
6.根据权利要求5所述的地震数据处理中...
【专利技术属性】
技术研发人员：郭清华，杨祥森，杨子兴，
申请(专利权)人：中国石油化工股份有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。