买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

龙图腾网获悉应急管理部国家自然灾害防治研究院;清华大学申请的专利地震成像走时数据处理方法及装置获国家发明授权专利权，本发明授权专利权由国家知识产权局授予，授权公告号为：CN114398780B 。

龙图腾网通过国家知识产权局官网在2025-09-23发布的发明授权授权公告中获悉：该发明授权的专利申请号/专利号为：202210024491.8，技术领域涉及：G06F30/20；该发明授权地震成像走时数据处理方法及装置是由刘少林;杨顶辉;汪文帅;申文豪设计研发完成，并于2022-01-10向国家知识产权局提交的专利申请。

本地震成像走时数据处理方法及装置在说明书摘要公布了：本发明提供了一种地震成像走时数据处理方法及装置，地震成像走时数据处理方法包括：根据目标工区的地震成像数据建立所述目标工区的地震速度结构模型；建立所述地震成像数据的程函方程；根据所述地震速度结构模型、所述程函方程以及预先建立的迎风差分格式数据确定所述地震成像走时数据的地震初至到时。本发明提供的地震成像走时数据处理方法及装置，能够代替手动处理地震成像走时数据，实现了地震成像走时数据处理的全自动化，且具有效率高、精度优的优点。

本发明授权地震成像走时数据处理方法及装置在权利要求书中公布了：1.一种地震成像走时数据处理方法，其特征在于，包括： 根据目标工区的地震成像数据建立所述目标工区的地震速度结构模型； 建立所述地震成像数据的程函方程； 根据所述地震速度结构模型、所述程函方程以及预先建立的迎风差分格式数据确定所述地震成像走时数据的地震初至到时； 所述建立所述地震成像数据的程函方程包括： 根据所述地震成像数据的空间走时场、空间梯度算子、介质慢度以及所述目标工区的基础数据建立所述程函方程，具体地： 根据公式1建立程函方程： 其中T为空间走时场，为空间梯度算子，S为介质慢度；r，θ，为球坐标的三个方向坐标，r表示空间任意点与地心的距离，θ表示余纬，表示经度； 建立所述迎风差分格式数据的方法包括： 根据所述地震成像数据的空间走时场以及其网格表征参数确定迎风差分算子； 根据所述迎风差分算子建立所述迎风差分格式数据，具体地根据公式2实现： 其中D表示迎风差分算子，算子的下标r、θ和分别表示作用于r、θ和方向，走时场T的下标表示网格点编号；δr、δθ和为空间网格在r、θ和方向的网格离散间隔；ri表示r方向第i个网格点与地心的距离； 所述根据所述地震速度结构模型、所述程函方程以及预先建立的迎风差分格式数据确定所述地震成像走时数据的地震初至到时，包括： 根据所述目标工区的基础数据、所述迎风差分格式数据以及所述地震速度结构模型求解所述程函方程，以确定所述地震初至到时； 所述根据所述目标工区的基础数据、所述迎风差分格式数据以及所述地震速度结构模型求解所述程函方程，以确定所述地震初至到时，包括： 根据所述迎风差分格式数据确定所述地震成像数据的在台站处的地震射线出射方向，具体地： 根据公式2计算得到的走时场可进一步得到台站处的地震射线出射方向，地震射线的出射方向可表示为： 其中e±1,f±1,g±1是围绕台站的8个网格点；和分别为球坐标的三个基向量；和分别为平行于r、θ和方向的向量； 按照所述地震射线出射方向，对所述地震成像数据进行投影变换操作，具体地： 由公式3将地震三分量其中和分别为台站处沿着纬度、经度和垂向的三个单位向量；UN、UE和UZ为三个单位向量上的位移；投影至射线出射方向经过投影变换后的地震波形表示为： 其中，Ul为地震出射方向上的位移； 利用长短视窗比值法，根据投影变换之后的地震数据确定比值法到时，具体地： 若地震数据的滤波范围为[f1,f2]，地震波信号的平均周期为： 传统长短时窗比法确定初至到时公式如下： 其中0ab为常数，一般a取值为0.1，b取值为3.0，τi为地震记录的任一时刻，max为取最大值运算，t3为传统长短时窗比法标定的到时； 利用回溯走时搜索算法，根据所述基础数据、地震速度结构模型以及所述比值法到时求解所述程函方程，以确定所述地震初至到时，具体地： 根据长短时窗比得到的走时t3向前回溯，当长短时窗比的值为平均长短时窗比的值时该处对应的走时被认为是精确到时，回溯走时搜索算法用公式表示如下： 其中abs为取绝对值运算，t4为回溯搜索得到的初至到时。

如需购买、转让、实施、许可或投资类似专利技术，可联系本专利的申请人或专利权人应急管理部国家自然灾害防治研究院;清华大学，其通讯地址为：100085 北京市海淀区清河安宁庄路1号；或者联系龙图腾网官方客服，联系龙图腾网可拨打电话0551-65771310或微信搜索“龙图腾网”。