地震初至波走时层析方法及装置制造方法及图纸

本发明专利技术提供了一种地震初至波走时层析方法及装置，其中，该方法包括：获取地震初至波的初至时间；建立初始速度模型；执行第一迭代步骤：按炮检对进行任务分组，并行进行射线追踪、走时残差和慢度更新量的计算；将多个射线密度汇总，得到射线密度汇总值；将多个走时残差汇总，得到走时残差汇总值；将多个慢度更新量汇总，得到慢度更新量汇总值；利用射线密度汇总值、走时残差汇总值和慢度更新量汇总值，对慢度模型进行更新；在满足第一预设条件时，执行更新速度模型的步骤；执行第二迭代步骤，直至得到满足第二预设条件的近地表速度模型，第二迭代步骤包括第一迭代步骤以及更新速度模型的步骤。上述技术方案提高了初至波走时层析的效率。

【技术实现步骤摘要】
地震初至波走时层析方法及装置
本专利技术涉及地球物理地震勘探
，特别涉及一种地震初至波走时层析方法及装置。
技术介绍
在地球物理地震勘探领域中，初至波走时层析是建立近地表速度模型的主要方法之一，其建立的近地表速度模型可以用于层析静校正量的计算，也可以作为叠前深度偏移和全波形反演的近地表速度模型。在数据量较小的情况下，传统的保存射线路径方法(将射线追踪出的射线路径保存为文件)或全内存操作方法可以满足计算效率的要求。但随着地震采集密度的日益增大，有效的地震道数成倍增加，这样一方面数据量成几何倍数增加，导致了射线路径文件非常庞大，有时达到TB级，全内存操作已成为不可能的事情。而频繁的读写这一庞大的射线路径文件，给计算效率带来了很大的挑战，导致计算效率显著降低，无法满足常规生产对计算效率的要求。针对上述问题，目前尚未提出有效的解决方案。
技术实现思路
本专利技术实施例提供了一种地震初至波走时层析方法，用以提高初至波走时层析的效率，该方法包括：获取地震初至波的初至时间；建立地震初至波走时层析反演的初始速度本文档来自技高网...

【技术保护点】
1.一种地震初至波走时层析方法，其特征在于，包括：/n获取地震初至波的初至时间；/n建立地震初至波走时层析反演的初始速度模型；/n根据所述初至时间、初始速度模型，执行如下第一迭代步骤：/n按炮检对，将射线追踪、走时残差和慢度更新量的计算任务进行分组，并行进行射线追踪、走时残差和慢度更新量的计算，得到多个射线密度、多个走时残差和多个慢度更新量；/n将多个射线密度进行汇总，得到射线密度汇总值；将多个走时残差进行汇总，得到走时残差汇总值；将多个慢度更新量进行汇总，得到慢度更新量汇总值；/n利用所述射线密度汇总值、走时残差汇总值和慢度更新量汇总值，对慢度模型进行更新，得到更新后的慢度模型；/n在满足第...

【技术特征摘要】
1.一种地震初至波走时层析方法，其特征在于，包括：
获取地震初至波的初至时间；
建立地震初至波走时层析反演的初始速度模型；
根据所述初至时间、初始速度模型，执行如下第一迭代步骤：
按炮检对，将射线追踪、走时残差和慢度更新量的计算任务进行分组，并行进行射线追踪、走时残差和慢度更新量的计算，得到多个射线密度、多个走时残差和多个慢度更新量；
将多个射线密度进行汇总，得到射线密度汇总值；将多个走时残差进行汇总，得到走时残差汇总值；将多个慢度更新量进行汇总，得到慢度更新量汇总值；
利用所述射线密度汇总值、走时残差汇总值和慢度更新量汇总值，对慢度模型进行更新，得到更新后的慢度模型；
在满足第一预设条件时，根据更新后的慢度模型，执行更新速度模型的步骤；
执行第二迭代步骤，直至得到满足第二预设条件的近地表速度模型，所述第二迭代步骤包括所述第一迭代步骤以及更新速度模型的步骤；所述近地表速度模型用于地震初至波走时层析静校正量的计算。


2.如权利要求1所述的地震初至波走时层析方法，其特征在于，按炮检对，将射线追踪、走时残差和慢度更新量的计算任务进行分组，并行进行射线追踪、走时残差和慢度更新量的计算，得到多个射线密度、多个走时残差和多个慢度更新量，包括：
根据CPU核数，按炮检对，将射线追踪、走时残差和慢度更新量的计算任务进行分组，每一CPU负责一组的计算任务；每一CPU根据组内每一炮检对对应的初至时间，进行射线追踪、走时残差和慢度更新量的计算，该计算的步骤包括：
在当前速度模型下进行射线追踪，得到当前炮检对之间的射线穿过每个网格的射线路径长度；
根据所述射线路径长度，以及当前炮检对对应的地震道初至时间，计算当前慢度模型下的所述地震道的走时残差；
根据所述走时残差和射线路径长度，计算当前射线对每个网格的慢度更新量；
对每一任务组内的网格慢度模型的慢度更新量进行汇总，得到任务组内所有数据对每个网格慢度模型慢度更新量的总贡献量；对每个网格穿过的射线条数进行汇总，得到任务组内对应的射线条数；对每个网格的走时残差的平方和信息进行汇总，得到任务组内对应的走时残差的平方和信息。


3.如权利要求2所述的地震初至波走时层析方法，其特征在于，将多个慢度更新量进行汇总，得到慢度更新量汇总值，包括：将每个CPU得到的任务组内所有数据对每个网格慢度模型慢度更新量的总贡献量进行累加，得到慢度更新量汇总值；
将多个射线密度进行汇总，得到射线密度汇总值，包括：将每个CPU得到的任务组内对应的射线条数进行累加，得到射线密度汇总值；
利用所述射线密度汇总值、走时残差汇总值和慢度更新量汇总值，对慢度模型进行更新，得到更新后的慢度模型，包括：
将慢度更新量汇总值除以射线密度汇总值，得到每个网格的慢度更新量；
根据所述每个网格的慢度更新量对慢度模型进行更新，得到更新后的慢度模型；根据走时残差的平方和信息或迭代次数，确定是否继续更新慢度模型。


4.如权利要求1所述的地震初至波走时层析方法，其特征在于，在满足第一预设条件时，根据更新后的慢度模型，执行更新速度模型的步骤，包括：
在第一预设的迭代次数大于第一预设值，或前后两次迭代的走时残差的均方根误差小于第二预设值时，根据更新后的慢度模型，执行更新速度模型的步骤。


5.如权利要求1所述的地震初至...

【专利技术属性】
技术研发人员：马青坡，黄兴贵，穆志平，王彦铎，潘英杰，
申请(专利权)人：中国石油天然气集团有限公司，中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：北京;11