基于地震初至的剩余静校正量处理方法及系统.pdf

《基于地震初至的剩余静校正量处理方法及系统.pdf》由会员分享，可在线阅读，更多相关《基于地震初至的剩余静校正量处理方法及系统.pdf（28页完成版）》请在专利查询网上搜索。

1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 202010703938.5 (22)申请日 2020.07.21 (71)申请人 中国石油天然气集团有限公司 地址 100007 北京市东城区东直门北大街9 号 申请人 中国石油集团东方地球物理勘探有 限责任公司 (72)发明人 冯发全 (74)专利代理机构 北京三友知识产权代理有限 公司 11127 代理人 王天尧吴学锋 (51)Int.Cl. G01V 1/36(2006.01) (54)发明名称 一种基于地震初至的剩余静校正量处理方 法及系统 (57)摘要 本发明提出了一。

2、种基于地震初至的剩余静 校正量处理方法及系统， 该方法包括： 获取初始 地震初至数据； 将基准面静校正施加到初始地震 初至数据， 并转换为三维共中心点地震初至数 据； 设定偏移距范围， 收集偏移距和地震初至； 对 偏移距和地震初至进行拟合， 得到每一个检波点 和炮点的剩余静校正初始值； 根据设定的迭代顺 序、 迭代次数及每次迭代的误差阈值进行剩余静 校正的迭代处理， 在迭代处理完成后， 将检波点 和炮点的剩余静校正量， 累加至基准面静校正， 得到每个检波点和炮点的基准面静校正量。 该方 法及系统可以改善基准面静校正的耦合能力， 准 确得到较大或超大剩余静校正量， 提高静校正的 精度， 改善地震。

3、资料的成像质量， 确保地震资料 的处理效果。 权利要求书5页 说明书12页 附图10页 CN 111965702 A 2020.11.20 CN 111965702 A 1.一种基于地震初至的剩余静校正量处理方法， 其特征在于， 该方法包括： 获取初始地震初至数据； 将基准面静校正施加到所述初始地震初至数据， 并转换为三维共中心点地震初至数 据； 设定每一个共中心点初至数据的偏移距范围， 在所述三维共中心点地震初至数据中， 分别收集每一个共中心点地震初至道集在所述偏移距范围内的偏移距和地震初至； 对所述偏移距和地震初至进行拟合， 得到每一个检波点和炮点的剩余静校正初始值； 根据设定的检波点与炮。

4、点的迭代顺序、 迭代次数及每次迭代的误差阈值进行剩余静校 正的迭代处理， 其中， 在每一次迭代处理时， 按照检波点与炮点的迭代顺序， 根据上一次迭代得到的所述检 波点和炮点的剩余静校正量， 分别对地震初至加上对应的检波点和炮点的剩余时差， 将偏 移距及加上剩余时差后的地震初至进行曲线拟合得到地震道拟合初至； 其中， 第一次迭代 处理时采用检波点和炮点的剩余静校正初始值； 根据所述地震道拟合初至获取地震道的初至拟合时差； 选取所述初至拟合时差小于本次迭代的误差阈值对应的地震道， 分别得到地震道对应 的本次迭代中检波点和炮点的累积剩余静校正量及累加次数； 根据检波点和炮点的累积静校正量及累加次数，。

5、 分别得到检波点和炮点的静校正平均 值作为本次迭代处理的检波点和炮点的剩余静校正量； 在迭代处理完成后， 将最后一次迭代处理后得到的检波点和炮点的剩余静校正量， 累 加至每个检波点和炮点的基准面静校正， 得到每个检波点和炮点的基准面静校正量。 2.根据权利要求1所述的基于地震初至的剩余静校正量处理方法， 其特征在于， 获取初 始地震初至数据， 包括： 整理野外生产资料， 获取三维SPS文件， 其中包括： 检波点文件、 炮点文件和关系文件； 根据所述SPS文件， 拾取每炮的初至时间， 获取初始地震初至数据。 3.根据权利要求1所述的基于地震初至的剩余静校正量处理方法， 其特征在于， 将基准 面静。

6、校正施加到所述初始地震初至数据， 并转换为三维共中心点地震初至数据， 包括： 将基准面静校正施加到所述初始地震初至数据， 得到新的地震初至数据； 采用人机交互及多项式拟合剔除所述新的地震初至数据中的异常地震初至； 根据预设的地震面元， 将剔除完成后的新的地震初至数据， 转换为三维共中心点地震 初至数据。 4.根据权利要求1所述的基于地震初至的剩余静校正量处理方法， 其特征在于， 设定每 一个共中心点的偏移距范围， 包括： 按一定空间间隔指定或分段指定共中心点初至数据的偏移距范围数据， 采用空间内插 方法得到每一个共中心点初至数据的偏移距范围。 5.根据权利要求1所述的基于地震初至的剩余静校正量。

7、处理方法， 其特征在于， 对所述 偏移距和地震初至进行拟合， 得到每一个检波点和炮点的剩余静校正初始值， 包括： 对偏移距和地震初至进行多项式拟合， 得到地震道拟合初至； 根据地震道拟合初至获取地震道的初至拟合时差； 选取初至拟合时差小于预设的误差阈值的地震道， 采用二分法分别计算地震道对应的 权利要求书 1/5 页 2 CN 111965702 A 2 检波点和炮点的累积时差， 同时记录检波点和炮点的累积时差的累加次数； 根据检波点的累积时差和累加次数， 计算时差平均值作为检波点剩余静校正的初始 值； 根据炮点的累积时差和累加次数， 计算时差平均值作为炮点剩余静校正的初始值。 6.根据权利要。

8、求5所述的基于地震初至的剩余静校正量处理方法， 其特征在于， 拟合的 次数为1至6次。 7.根据权利要求5所述的基于地震初至的剩余静校正量处理方法， 其特征在于， 根据设 定的检波点与炮点的迭代顺序、 迭代次数及每次迭代的误差阈值进行剩余静校正的迭代处 理， 包括： 在设定的检波点与炮点的迭代顺序为先炮点后检波点时， 利用下式计算得到第k次迭 代时检波点的初始剩余静校正量： 其中， Tri(k)是第k次迭代处理时检波点i的初始剩余静校正量； k1、 2、 、 N， N是设定的迭代次数； Tr(k-1)是上一次迭代得到的检波点的剩余静校正量， 在k1时， Tr(0)采用检波点的剩余 静校正初始值。

9、； 是第k次迭代时地震道的检波点i及炮点j的初至拟合时差； 根据本次迭代得到的检波点i的初始剩余静校正量Tri(k)， 进一步得到每个检波点的累 积初始剩余静校正量和累加次数， 并根据每个检波点的累积初始剩余静校正量和累加次 数， 得到静校正平均值作为第k次迭代时检波点的初始剩余静校正量Tr(k)； 利用下式计算得到第k次迭代时炮点的剩余静校正量： 其中，是炮点j第k次迭代的剩余静校正量； Ts(k-1)是上一次迭代得到的炮点的剩余静校正量， 在k1时， Ts(0)采用炮点的剩余静校 正初始值； 是第k次迭代时地震道的检波点i及炮点j的初至拟合时差； Tr(k)是第k次迭代时检波点的初始剩余静。

10、校正量； Tr(k-1)是上一次迭代得到的检波点的剩余静校正量， 在k1时， Tr0采用检波点的剩余 静校正初始值； 根据本次迭代得到的炮点j的剩余静校正量进一步得到每个炮点的累积剩余静校 正量和累加次数， 并根据每个炮点的累积剩余静校正量和累加次数， 得到静校正平均值作 为第k次迭代时炮点的剩余静校正量Ts(k)； 利用下式计算得到第k次迭代时检波点的剩余静校正量： 其中， Tri(k)是检波点i第k次迭代的剩余静校正量； 权利要求书 2/5 页 3 CN 111965702 A 3 Tr(k-1)是上一次迭代得到的检波点的剩余静校正量， 在k1时， Tr(0)采用检波点的剩余 静校正初始值。

11、； 是第k次迭代时地震道的检波点i及炮点j的初至拟合时差； Ts(k)是第k次迭代时炮点的初始剩余静校正量； Ts(k-1)是上一次迭代得到的炮点的剩余静校正量， 在k1时， Ts0采用炮点的剩余静校 正初始值； 根据本次迭代得到的检波点i的剩余静校正量Tri(k)， 进一步得到每个检波点的累积剩 余静校正量和累加次数， 并根据每个检波点的累积剩余静校正量和累加次数， 得到静校正 平均值作为第k次迭代时检波点的剩余静校正量Tr(k)。 8.根据权利要求5所述的基于地震初至的剩余静校正量处理方法， 其特征在于， 根据设 定的检波点与炮点的迭代顺序、 迭代次数及每次迭代的误差阈值进行剩余静校正的迭。

12、代处 理， 包括： 在设定的检波点与炮点的迭代顺序为先检波点后炮点时， 利用下式计算得到第k次迭 代时炮点的初始剩余静校正量： 其中， Tsj(k)是第k次迭代处理时炮点j的初始剩余静校正量； k1、 2、 、 N， N是设定的迭代次数； Ts(k-1)是上一次迭代得到的炮点的剩余静校正量， 在k1时， Ts(0)采用炮点的剩余静校 正初始值； 是第k次迭代时地震道的检波点i及炮点j的初至拟合时差； 根据本次迭代得到的炮点j的初始剩余静校正量Tsj(k)， 进一步得到每个炮点的累积初 始剩余静校正量和累加次数， 并根据每个炮点的累积初始剩余静校正量和累加次数， 得到 静校正平均值作为第k次迭代。

13、时炮点的初始剩余静校正量Ts(k)； 利用下式计算得到第k次迭代时检波点的剩余静校正量： 其中， Tri(k)是检波点i第k次迭代的剩余静校正量； Tr(k-1)是上一次迭代得到的检波点的剩余静校正量， 在k1时， Tr(0)采用检波点的剩余 静校正初始值； 是第k次迭代时地震道的检波点i及炮点j的初至拟合时差； Ts(k)是第k次迭代时炮点的初始剩余静校正量； Ts(k-1)是上一次迭代得到的炮点的剩余静校正量， 在k1时， Ts0采用炮点的剩余静校 正初始值； 根据本次迭代得到的检波点i的剩余静校正量Tri(k)， 进一步得到每个检波点的累积剩 余静校正量和累加次数， 并根据每个检波点的累。

14、积剩余静校正量和累加次数， 得到静校正 平均值作为第k次迭代时检波点的剩余静校正量Tr(k)； 利用下式计算得到第k次迭代时炮点的剩余静校正量： 权利要求书 3/5 页 4 CN 111965702 A 4 其中，是炮点j第k次迭代的剩余静校正量； Ts(k-1)是上一次迭代得到的炮点的剩余静校正量， 在k1时， Ts(0)采用炮点的剩余静校 正初始值； 是第k次迭代时地震道的检波点i及炮点j的初至拟合时差； Tr(k)是第k次迭代时检波点的初始剩余静校正量； Tr(k-1)是上一次迭代得到的检波点的剩余静校正量， 在k1时， Tr0采用检波点的剩余 静校正初始值； 根据本次迭代得到的炮点j的。

15、剩余静校正量进一步得到每个炮点的累积剩余静校 正量和累加次数， 并根据每个炮点的累积剩余静校正量和累加次数， 得到静校正平均值作 为第k次迭代时炮点的剩余静校正量Ts(k)。 9.一种基于地震初至的剩余静校正量处理系统， 其特征在于， 该系统包括： 数据获取模块， 用于获取初始地震初至数据； 数据处理模块， 用于将基准面静校正施加到所述初始地震初至数据， 并转换为三维共 中心点地震初至数据； 偏移距和地震初至收集模块， 设定每一个共中心点的偏移距范围， 在所述三维共中心 点地震初至数据中， 分别收集每一个共中心点地震初至道集在所述偏移距范围内的偏移距 和地震初至； 拟合模块， 用于对所述偏移距。

16、和地震初至进行拟合， 得到每一个检波点和炮点的剩余 静校正初始值； 迭代处理模块， 用于根据设定的检波点与炮点的迭代顺序、 迭代次数及每次迭代的误 差阈值进行剩余静校正的迭代处理， 其中， 包括： 拟合单元， 用于在每一次迭代处理时， 按照检波点与炮点的迭代顺序， 根据上一次迭代 得到的所述检波点和炮点的剩余静校正量， 分别对地震初至加上对应的检波点和炮点的剩 余时差， 将偏移距及加上剩余时差后的地震初至进行曲线拟合得到地震道拟合初至； 其中， 第一次迭代处理时采用检波点和炮点的剩余静校正初始值； 初至拟合时差获取单元， 用于根据所述地震道拟合初至获取地震道的初至拟合时差； 累积剩余静校正量计。

17、算单元， 用于选取所述初至拟合时差小于本次迭代的误差阈值对 应的地震道， 分别得到地震道对应的本次迭代中检波点和炮点的累积剩余静校正量及累加 次数； 剩余静校正量计算单元， 用于根据检波点和炮点的累积静校正量及累加次数， 分别得 到检波点和炮点的静校正平均值作为本次迭代处理的检波点和炮点的剩余静校正量； 基准面静校正量计算模块， 用于在迭代处理完成后， 将最后一次迭代处理后得到的检 波点和炮点的剩余静校正量， 累加至每个检波点和炮点的基准面静校正， 得到每个检波点 和炮点的基准面静校正量。 10.一种计算机设备， 包括存储器、 处理器及存储在存储器上并可在处理器上运行的计 算机程序， 其特征在。

18、于， 所述处理器执行所述计算机程序时实现权利要求1至8任一所述方 权利要求书 4/5 页 5 CN 111965702 A 5 法。 11.一种计算机可读存储介质， 其特征在于， 所述计算机可读存储介质存储有计算机程 序， 所述计算机程序被处理器执行时实现权利要求1至8任一所述方法。 权利要求书 5/5 页 6 CN 111965702 A 6 一种基于地震初至的剩余静校正量处理方法及系统 技术领域 0001 本发明涉及地球物理勘探技术领域， 尤指一种基于地震初至的剩余静校正量处理 方法及系统。 背景技术 0002 静校正是目前地震资料处理中要做的一个重要基础工作， 对确保地震成像质量至 关重。

19、要。 0003 目前， 基准面静校正是先建立近地表模型然后根据给定的相关参数进行计算， 由 于在复杂地区， 模型精度较低， 在此模型上计算的基准面静校正的解在空间上并不耦合， 导 致地震资料成像较差或不能成像。 目前对基准面静校正的解是否在空间上耦合并不检查分 析， 而只是采用剩余静校正进行解决， 当基准面静校正在空间上耦合较差时， 剩余静校正较 大而很难求准， 目前的剩余静校正技术无法满足大剩余静校正的要求。 0004 综上来看， 亟需一种可以解决上述问题的剩余静校正量处理技术方案。 发明内容 0005 为解决上述问题， 本发明提出了一种基于地震初至的剩余静校正量处理方法及系 统。 通过在成。

20、像域进行静校正耦合求解， 得到大剩余静校正量， 改善基准面静校正的耦合能 力， 确保地震资料的处理效果， 以满足地震资料处理对静校正高精度的要求。 0006 在本发明实施例的第一方面， 提出了一种基于地震初至的剩余静校正量处理方 法， 该方法包括： 0007 获取初始地震初至数据； 0008 将基准面静校正施加到所述初始地震初至数据， 并转换为三维共中心点地震初至 数据； 0009 设定每一个共中心点初至数据的偏移距范围， 在所述三维共中心点地震初至数据 中， 分别收集每一个共中心点地震初至道集在所述偏移距范围内的偏移距和地震初至； 0010 对所述偏移距和地震初至进行拟合， 得到每一个检波点。

21、和炮点的剩余静校正初始 值； 0011 根据设定的检波点与炮点的迭代顺序、 迭代次数及每次迭代的误差阈值进行剩余 静校正的迭代处理， 其中， 0012 在每一次迭代处理时， 按照检波点与炮点的迭代顺序， 根据上一次迭代得到的所 述检波点和炮点的剩余静校正量， 分别对地震初至加上对应的检波点和炮点的剩余时差， 将偏移距及加上剩余时差后的地震初至进行曲线拟合得到地震道拟合初至； 0013 根据所述地震道拟合初至获取地震道的初至拟合时差； 其中， 第一次迭代处理时 采用检波点和炮点的剩余静校正初始值； 0014 选取所述初至拟合时差小于本次迭代的误差阈值对应的地震道， 分别得到地震道 对应的本次迭代。

22、中检波点和炮点的累积剩余静校正量及累加次数； 说明书 1/12 页 7 CN 111965702 A 7 0015 根据检波点和炮点的累积静校正量及累加次数， 分别得到检波点和炮点的静校正 平均值作为本次迭代处理的检波点和炮点的剩余静校正量； 0016 在迭代处理完成后， 将最后一次迭代处理后得到的检波点和炮点的剩余静校正 量， 累加至每个检波点和炮点的基准面静校正， 得到每个检波点和炮点的基准面静校正量。 0017 在本发明实施例的第二方面， 还提出了一种基于地震初至的剩余静校正量处理系 统， 该系统包括： 0018 数据获取模块， 用于获取初始地震初至数据； 0019 数据处理模块， 用于。

23、将基准面静校正施加到所述初始地震初至数据， 并转换为三 维共中心点地震初至数据； 0020 偏移距和地震初至收集模块， 设定每一个共中心点的偏移距范围， 在所述三维共 中心点地震初至数据中， 分别收集每一个共中心点地震初至道集在所述偏移距范围内的偏 移距和地震初至； 0021 拟合模块， 用于对所述偏移距和地震初至进行拟合， 得到每一个检波点和炮点的 剩余静校正初始值； 0022 迭代处理模块， 用于根据设定的检波点与炮点的迭代顺序、 迭代次数及每次迭代 的误差阈值进行剩余静校正的迭代处理， 其中， 包括： 0023 拟合单元， 用于在每一次迭代处理时， 按照检波点与炮点的迭代顺序， 根据上一。

24、次 迭代得到的所述检波点和炮点的剩余静校正量， 分别对地震初至加上对应的检波点和炮点 的剩余时差， 将偏移距及加上剩余时差后的地震初至进行曲线拟合得到地震道拟合初至； 其中， 第一次迭代处理时采用检波点和炮点的剩余静校正初始值； 0024 初至拟合时差获取单元， 用于根据所述地震道拟合初至获取地震道的初至拟合时 差； 0025 累积剩余静校正量计算单元， 用于选取所述初至拟合时差小于本次迭代的误差阈 值对应的地震道， 分别得到地震道对应的本次迭代中检波点和炮点的累积剩余静校正量及 累加次数； 0026 剩余静校正量计算单元， 用于根据检波点和炮点的累积静校正量及累加次数， 分 别得到检波点和炮。

25、点的静校正平均值作为本次迭代处理的检波点和炮点的剩余静校正量； 0027 基准面静校正量计算模块， 用于在迭代处理完成后， 将最后一次迭代处理后得到 的检波点和炮点的剩余静校正量， 累加至每个检波点和炮点的基准面静校正， 得到每个检 波点和炮点的基准面静校正量。 0028 在本发明实施例的第三方面， 还提出了一种计算机设备， 包括存储器、 处理器及存 储在存储器上并可在处理器上运行的计算机程序， 所述处理器执行所述计算机程序时实现 基于地震初至的剩余静校正量处理方法。 0029 在本发明实施例的第四方面， 还提出了一种计算机可读存储介质， 所述计算机可 读存储介质存储有计算机程序， 所述计算机。

26、程序被处理器执行时实现基于地震初至的剩余 静校正量处理方法。 0030 本发明提出的基于地震初至的剩余静校正量处理方法及系统可以改善基准面静 校正的耦合能力， 准确得到较大或超大剩余静校正量， 提高静校正的精度， 改善地震资料的 成像质量， 确保地震资料的处理效果。 说明书 2/12 页 8 CN 111965702 A 8 附图说明 0031 为了更清楚地说明本申请实施例技术方案， 下面将对实施例描述中所需要使用的 附图作简单地介绍， 显而易见地， 下面描述中的附图是本申请的一些实施例， 对于本领域普 通技术人员来讲， 在不付出创造性劳动的前提下， 还可以根据这些附图获得其它的附图。 003。

27、2 图1是本发明一实施例的基于地震初至的剩余静校正量处理方法流程示意图。 0033 图2是本发明一实施例的剩余静校正迭代处理的详细流程示意图。 0034 图3是本发明一具体实施例的基于地震初至的剩余静校正量处理方法的流程示意 图。 0035 图4是本发明一具体实施例的LX地区三维某检波点线基准面静校正曲线示意图。 0036 图5是本发明一具体实施例的应用本发明方法后的剩余静校正静校正曲线示意 图。 0037 图6是本发明一具体实施例的LX地区三维应用剩余静校正前的单炮记录示意图。 0038 图7是本发明一具体实施例的应用本发明方法后的剩余静校正的单炮记录示意 图。 0039 图8是本发明一具体。

28、实施例的应用本发明方法后的剩余静校正的叠加剖面示意 图。 0040 图9是本发明一具体实施例的LX地区三维应用剩余静校正前的叠加剖面示意图。 0041 图10是本发明一实施例的基于地震初至的剩余静校正量处理系统的结构示意图。 0042 图11是本发明一实施例的迭代处理模块的结构示意图。 0043 图12是本发明一实施例的数据获取模块的结构示意图。 0044 图13是本发明一实施例的数据处理模块的结构示意图。 0045 图14是本发明一实施例的拟合模块的结构示意图。 0046 图15是本发明一实施例的计算机设备结构示意图。 具体实施方式 0047 下面将参考若干示例性实施方式来描述本发明的原理和。

29、精神。 应当理解， 给出这 些实施方式仅仅是为了使本领域技术人员能够更好地理解进而实现本发明， 而并非以任何 方式限制本发明的范围。 相反， 提供这些实施方式是为了使本公开更加透彻和完整， 并且能 够将本公开的范围完整地传达给本领域的技术人员。 0048 本领域技术人员知道， 本发明的实施方式可以实现为一种系统、 装置、 设备、 方法 或计算机程序产品。 因此， 本公开可以具体实现为以下形式， 即： 完全的硬件、 完全的软件 (包括固件、 驻留软件、 微代码等)， 或者硬件和软件结合的形式。 0049 根据本发明的实施方式， 提出了一种基于地震初至的剩余静校正量处理方法及系 统。 该方法及系统。

30、在进行剩余静校正前剔除较大的异常初至， 对于不同共中心点使用的偏 移距是空变的， 一个共中心点道集上使用的偏移距有时是多个偏移距范围， 并且不同偏移 距范围可以部分重合； 在处理过程中， 使用较高次数的多项式拟合， 提高地震时差的计算精 度， 应用改进的高斯赛德尔方法进行迭代计算检波点、 炮点的剩余静校正量； 整体方案处理 速度快， 能较好解决较大(超大)的剩余静校正， 改善基准面静校正的耦合能力， 同时便于处 理系统剩余静校正的发挥， 进一步提高静校正的精度， 改善地震资料的成像质量， 确保地震 说明书 3/12 页 9 CN 111965702 A 9 资料的处理效果； 在复杂地区的表层建。

31、模和静校正工作中具有广泛的适用性和通用性。 0050 下面参考本发明的若干代表性实施方式， 详细阐释本发明的原理和精神。 0051 图1是本发明一实施例的基于地震初至的剩余静校正量处理方法流程示意图。 如 图1所示， 该方法包括： 0052 步骤S101， 获取初始地震初至数据。 0053 步骤S102， 将基准面静校正施加到所述初始地震初至数据， 并转换为三维共中心 点地震初至数据。 0054 步骤S103， 设定每一个共中心点初至数据的偏移距范围， 在所述三维共中心点地 震初至数据中， 分别收集每一个共中心点地震初至道集在所述偏移距范围内的偏移距和地 震初至。 0055 步骤S104， 对。

32、所述偏移距和地震初至进行拟合， 得到每一个检波点和炮点的剩余 静校正初始值。 0056 步骤S105， 根据设定的检波点与炮点的迭代顺序、 迭代次数及每次迭代的误差阈 值进行剩余静校正的迭代处理。 0057 其中， 参考图2， 为本发明一实施例的剩余静校正迭代处理的详细流程示意图。 如 图2所示， 步骤S105包括： 0058 步骤S1051， 在每一次迭代处理时， 按照检波点与炮点的迭代顺序， 根据上一次迭 代得到的所述检波点和炮点的剩余静校正量， 分别对地震初至加上对应的检波点和炮点的 剩余时差， 将偏移距及加上剩余时差后的地震初至进行曲线拟合得到地震道拟合初至。 其 中， 第一次迭代处理。

33、时采用检波点和炮点的剩余静校正初始值。 0059 步骤S1052， 根据所述地震道拟合初至获取地震道的初至拟合时差。 0060 步骤S1053， 选取所述初至拟合时差小于本次迭代的误差阈值对应的地震道， 分别 得到地震道对应的本次迭代中检波点和炮点的累积剩余静校正量及累加次数。 0061 步骤S1054， 根据检波点和炮点的累积静校正量及累加次数， 分别得到检波点和炮 点的静校正平均值作为本次迭代处理的检波点和炮点的剩余静校正量。 0062 步骤S106， 在迭代处理完成后， 将最后一次迭代处理后得到的检波点和炮点的剩 余静校正量， 累加至每个检波点和炮点的基准面静校正， 得到每个检波点和炮点。

34、的基准面 静校正量。 0063 为了对上述基于地震初至的剩余静校正量处理方法进行更为清楚的解释， 下面结 合具体的实施例来进行说明。 0064 如图3所示， 为本发明一具体实施例的基于地震初至的剩余静校正量处理方法的 流程示意图。 如图3所示， 处理流程为： 0065 S301： 0066 整理野外生产资料， 获取三维SPS文件， 其中包括： 检波点文件、 炮点文件和关系文 件， 检波点文件按检波点桩号升序排序、 炮点文件按炮点桩号升序排序。 0067 S302： 0068 根据所述SPS文件， 拾取每炮的初至时间， 获取初始地震初至数据。 0069 S303： 0070 将基准面静校正施加到。

35、所述初始地震初至数据， 得到新的地震初至数据。 说明书 4/12 页 10 CN 111965702 A 10 0071 S304： 0072 采用人机交互及多项式拟合剔除所述新的地震初至数据中的异常地震初至， 以提 高基于地震初至数据进行剩余静校正处理的结果的准确性。 0073 S305： 0074 设定地震面元， 根据地震面元将剔除完成后的新的地震初至数据转换为三维共中 心点地震初至数据。 0075 S306： 0076 按一定空间间隔指定或分段指定共中心点初至数据的偏移距范围数据， 采用空间 内插方法得到每一个共中心点初至数据的偏移距范围。 0077 对于分段指定的偏移距范围数据可以有部。

36、分重合， 例如， 500m-2500m， 2000m- 4000m。 0078 S307： 0079 指定要使用的共中心点地震初至的起始和终止CMP线号、 CMP线号增量、 起始和终 止CMP点号、 CMP点号增量， 设定初始误差阈值。 0080 S308： 0081 在选定的一个共中心点地震初至道集上， 收集或分段收集在指定偏移距范围内的 偏移距和地震初至。 0082 S309： 0083 对收集的偏移距和地震初至进行多项式拟合， 得到各地震道的拟合初至， 并计算 各地震道的初至拟合时差， 其中， 拟合的次数为1至6次。 通过使用较高次数的多项式拟合， 可以提高地震时差的计算精度。 0084。

37、 S310： 0085 对在小于误差阈值范围内的初至拟合时差采用二分法分别计算该地震道对应的 检波点和炮点的累积时差， 同时记录检波点和炮点的累积时差的累加次数。 0086 S311： 0087 根据检波点的累积时差和累加次数， 计算时差平均值作为检波点剩余静校正的初 始值。 0088 S312： 0089 根据炮点的累积时差和累加次数， 计算时差平均值作为炮点剩余静校正的初始 值。 0090 在得到检波点和炮点剩余静校正的初始值后， 本发明将利用改进的高斯赛德尔方 法进行迭代计算检波点、 炮点的剩余静校正量， 具体过程为： 0091 S313： 0092 设置迭代次数、 每次迭代误差阈值， 。

38、并确定迭检波点和炮点的迭代顺序。 例如， 迭 代次数可以设置为3次， 相应的， 每次迭代误差阈值分别设置为30ms、 20ms、 10ms， 迭代顺序 可以为： 先迭代炮点、 后迭代检波点(或者先迭代检波点、 后迭代炮点)。 0093 S314： 0094 进行剩余静校正的迭代处理， 具体过程为： 0095 在每一次迭代处理时， 按照检波点与炮点的迭代顺序， 根据上一次迭代得到的所 说明书 5/12 页 11 CN 111965702 A 11 述检波点和炮点的剩余静校正量， 分别对地震初至加上对应的检波点和炮点的剩余时差， 将偏移距及加上剩余时差后的地震初至进行曲线拟合得到地震道拟合初至； 。

39、(拟合的次数 为1至6次)。 其中， 第一次迭代处理时采用检波点和炮点的剩余静校正初始值。 0096 根据所述地震道拟合初至获取地震道的初至拟合时差。 0097 选取所述初至拟合时差小于本次迭代的误差阈值对应的地震道， 分别得到地震道 对应的本次迭代中检波点和炮点的累积剩余静校正量及累加次数。 0098 根据检波点和炮点的累积静校正量及累加次数， 分别得到检波点和炮点的静校正 平均值作为本次迭代处理的检波点和炮点的剩余静校正量。 0099 对于两种迭代顺序， 以下分别说明计算检波点和炮点的剩余静校正量的详细过 程： 0100 在设定的检波点与炮点的迭代顺序为先炮点后检波点时， 首先， 利用下式。

40、计算得 到第k次迭代时检波点的初始剩余静校正量： 0101 0102其中，是第k次迭代处理时检波点i的初始剩余静校正量； 0103 k1、 2、 、 N， N是设定的迭代次数； 0104 Tr(k-1)是上一次迭代得到的检波点的剩余静校正量， 在k1时， Tr(0)采用检波点的 剩余静校正初始值； 0105是第k次迭代时地震道的检波点i及炮点j的初至拟合时差； 0106 根据本次迭代得到的检波点i的初始剩余静校正量Tri(k)， 进一步得到每个检波 点的累积初始剩余静校正量和累加次数， 并根据每个检波点的累积初始剩余静校正量和累 加次数， 得到静校正平均值作为第k次迭代时检波点的初始剩余静校正。

41、量Tr(k)； 0107 利用下式计算得到第k次迭代时炮点的剩余静校正量： 0108 0109其中，是炮点j第k次迭代的剩余静校正量； 0110 Ts(k-1)是上一次迭代得到的炮点的剩余静校正量， 在k1时， Ts(0)采用炮点的剩余 静校正初始值； 0111是第k次迭代时地震道的检波点i及炮点j的初至拟合时差； 0112 Tr(k)是第k次迭代时检波点的初始剩余静校正量； 0113 Tr(k-1)是上一次迭代得到的检波点的剩余静校正量， 在k1时， Tr0采用检波点的 剩余静校正初始值； 0114根据本次迭代得到的炮点j的剩余静校正量进一步得到每个炮点的累积剩 余静校正量和累加次数， 并根。

42、据每个炮点的累积剩余静校正量和累加次数， 得到静校正平 均值作为第k次迭代时炮点的剩余静校正量Ts(k)； 0115 利用下式计算得到第k次迭代时检波点的剩余静校正量： 0116 说明书 6/12 页 12 CN 111965702 A 12 0117 其中， Tri(k)是检波点i第k次迭代的剩余静校正量； 0118 Tr(k-1)是上一次迭代得到的检波点的剩余静校正量， 在k1时， Tr(0)采用检波点的 剩余静校正初始值； 0119是第k次迭代时地震道的检波点i及炮点j的初至拟合时差； 0120 Ts(k)是第k次迭代时炮点的初始剩余静校正量； 0121 Ts(k-1)是上一次迭代得到的。

43、炮点的剩余静校正量， 在k1时， Ts0采用炮点的剩余 静校正初始值； 0122 根据本次迭代得到的检波点i的剩余静校正量Tri(k)， 进一步得到每个检波点的累 积剩余静校正量和累加次数， 并根据每个检波点的累积剩余静校正量和累加次数， 得到静 校正平均值作为第k次迭代时检波点的剩余静校正量Tr(k)。 0123 在设定的检波点与炮点的迭代顺序为先检波点后炮点时， 首先， 利用下式计算得 到第k次迭代时炮点的初始剩余静校正量： 0124 0125 其中， Tsj(k)是第k次迭代处理时炮点j的初始剩余静校正量； 0126 k1、 2、 、 N， N是设定的迭代次数； 0127 Ts(k-1)。

44、是上一次迭代得到的炮点的剩余静校正量， 在k1时， Ts(0)采用炮点的剩余 静校正初始值； 0128是第k次迭代时地震道的检波点i及炮点j的初至拟合时差； 0129 根据本次迭代得到的炮点j的初始剩余静校正量Tsj(k)， 进一步得到每个炮点的 累积初始剩余静校正量和累加次数， 并根据每个炮点的累积初始剩余静校正量和累加次 数， 得到静校正平均值作为第k次迭代时炮点的初始剩余静校正量Ts(k)； 0130 利用下式计算得到第k次迭代时检波点的剩余静校正量： 0131 0132 其中， Tri(k)是检波点i第k次迭代的剩余静校正量； 0133 Tr(k-1)是上一次迭代得到的检波点的剩余静校。

45、正量， 在k1时， Tr(0)采用检波点的 剩余静校正初始值； 0134是第k次迭代时地震道的检波点i及炮点j的初至拟合时差； 0135 Ts(k)是第k次迭代时炮点的初始剩余静校正量； 0136 Ts(k-1)是上一次迭代得到的炮点的剩余静校正量， 在k1时， Ts0采用炮点的剩余 静校正初始值； 0137 根据本次迭代得到的检波点i的剩余静校正量Tri(k)， 进一步得到每个检波点的累 积剩余静校正量和累加次数， 并根据每个检波点的累积剩余静校正量和累加次数， 得到静 校正平均值作为第k次迭代时检波点的剩余静校正量Tr(k)； 0138 利用下式计算得到第k次迭代时炮点的剩余静校正量： 0。

46、139 说明书 7/12 页 13 CN 111965702 A 13 0140其中，是炮点j第k次迭代的剩余静校正量； 0141 Ts(k-1)是上一次迭代得到的炮点的剩余静校正量， 在k1时， Ts(0)采用炮点的剩余 静校正初始值； 0142是第k次迭代时地震道的检波点i及炮点j的初至拟合时差； 0143 Tr(k)是第k次迭代时检波点的初始剩余静校正量； 0144 Tr(k-1)是上一次迭代得到的检波点的剩余静校正量， 在k1时， Tr0采用检波点的 剩余静校正初始值； 0145根据本次迭代得到的炮点j的剩余静校正量进一步得到每个炮点的累积剩 余静校正量和累加次数， 并根据每个炮点的累。

47、积剩余静校正量和累加次数， 得到静校正平 均值作为第k次迭代时炮点的剩余静校正量Ts(k)。 0146 S315： 0147 在最后一次迭代完成后， 各检波点和炮点的剩余静校正量即为各检波点和炮点的 最终剩余静校正量， 将剩余静校正量累加到各检波点和炮点的基准面静校正量上， 作为各 检波点和炮点新的基准面静校正量； 该新的基准面静校正量可以应用于地震资料处理， 也 可以供其它方法进一步计算剩余静校正量。 0148 由于上述过程增加的计算量较低， 可以保证较快的处理速度， 通过静校正耦合求 解得到大剩余静校正量， 便于处理系统剩余静校正的发挥， 提高基准面静校正的耦合能力， 同时可以提高静校正的。

48、精度， 改善地震资料的成像质量， 确保地震资料的处理效果； 本方案 在复杂地区的表层建模和静校正工作中具有广泛的适用性和通用性。 0149 需要说明的是， 尽管在上述实施例及附图中以特定顺序描述了本发明方法的操 作， 但是， 这并非要求或者暗示必须按照该特定顺序来执行这些操作， 或是必须执行全部所 示的操作才能实现期望的结果。 附加地或备选地， 可以省略某些步骤， 将多个步骤合并为一 个步骤执行， 和/或将一个步骤分解为多个步骤执行。 0150 为了对上述基于地震初至的剩余静校正量处理方法进行更为清楚的解释， 下面结 合一个具体的实施例来进行说明， 然而值得注意的是该实施例仅是为了更好地说明本。

49、发 明， 并不构成对本发明不当的限定。 0151 以LX地区为例， 该地区属于复杂地区， 在该地区表层结构复杂， 由于表层建模困 难、 模型精度较差， 基准面静校正问题严重， 静校正耦合较差， 导致出现剩余静校正量较大， 常规剩余静校正无法解决， 影响地震资料的成像质量等问题， 三维基准面静校正存在较大 剩余时差不能满足地震资料的处理要求。 0152 对此， 结合图1至图3的步骤流程， 可以较好地获取剩余静校正量， 改善基准面静校 正的耦合能力， 同时便于处理系统剩余静校正的发挥， 进一步提高静校正的精度， 改善地震 资料的成像质量， 确保地震资料的处理效果。 0153 首先， 拾取地震初至并。

50、施加基准面静校正得到新的地震初至。 0154 设置地震面元大小20cm20cm， 按地震面元生成三维共中心点的地震初至文件。 0155 指定或分段指定每一个共中心点要使用的偏移距范围， 例如， 要使用的共中心点 地震初至的起始CMP线号1001、 终止CMP线号2185、 CMP线号增量2、 起始CMP点号1020、 终止 CMP点号3150、 CMP点号增量2， 初始误差阈值30ms。 说明书 8/12 页 14 CN 111965702 A 14 0156 在选定的一个共中心点地震初至道集上， 收集或分段收集在指定偏移距范围内的 偏移距和地震初至， 对收集的偏移距和地震初至进行多项式拟合得。