[发明专利]共中心点初至道集生成方法、计算剩余静校正方法及装置

说明书

本发明公开了一种共中心点初至道集生成方法、计算剩余静校正方法及装置。本发明将各共中心点的位置网格化，确定每个初至地震道在该共中心点中的顺序，并以此统计每个共中心点所在面元的初至地震道的数量，以及每个面元第一个初至地震道的存储位置，这样，在生成共中心点初至道集文件时，可利用面元第一个初至地震道的存储位置，以及自身的顺序，可快速确定每个面元所有初至地震道的地址的偏移量，从而快速完成各面元对应的各共中心点的所有初至地震道的数据的写入，以生成共中心点初至道集文件。本发明缩短了初至地震道数据的处理时间，提高了共中心点初至道集的生成效率，较好地满足了地震资料处理的时效要求。

技术领域

本发明涉及地球物理勘探数据处理技术领域，特别涉及一种共中心点初至道集生成方法、计算剩余静校正方法及装置。

背景技术

在地震数据中，多个炮点-检波点对，共享炮点和检波点的中心点，此类中心点称之为共中心点。

共中心点初至道集是地震资料处理常用的一个数据体，野外数据是共炮点道集，因此共中心点初至道集需要根据SPS文件(包含炮点文件、检波点文件及以及定义炮点-检波点关系的关系文件)自主生成。

最简单的生成方法就是从炮点中逐一寻找，提取相关信息，然后生成共中心点初至道集，但由于炮点数量多，这样做其计算量也是最大的，尤其是三维，其计算量更是十分惊人，计算时间较长。

目前在公开的文献中，未能检索到生成共中心点初至道集的公开方法。因此依照上述的简单分析，生成三维共中心点初至道集的计算量十分巨大，在普通微机上所需时间相当长，实际生产中很难满足地震资料处理的时效要求。

发明内容

鉴于上述问题，提出了本发明以便提供一种克服上述问题或者至少部分地解决上述问题的一种共中心点初至道集生成方法、计算剩余静校正方法及装置。

第一方面，本发明实施例提供一种共中心点初至道集的生成方法，包括：

根据测线的方位角，将各检波点和各炮点的第一坐标进行旋转，得到以所述测线为正东或正北的新坐标系中的第二坐标；

根据预设的面元大小，对所述新坐标系进行网格状划分，根据各检波点和各炮点的第二坐标，确定各检波点-炮点对的共中心点所在的面元；

确定每炮的初至地震道所在的共中心点，以及所述初至地震道在共享所述共中心点的所有初至地震道中的位置顺序；

根据每炮的初至地震道的共中心点所在的面元的标识，统计每个网格面元的共中心点包含的初至地震道的总道数；

根据每个面元中共中心点包含的初至地震道的总道数，按网格中各面元的顺序，依次为每个面元的共中心点的第一个初至地震道，分配其在共中心点初至道集文件中的记录位置，作为所述共中心点的起始记录位置；

针对每个初至地震道，按照所述初至地震道所在共中心点的所述起始记录位置和所述初至地震道在共中心点中的所述位置顺序，确定每个初至地震道在共中心点初至道集文件中的实际记录位置，并将所述初至地震道的相关信息，写入所述实际记录位置，以生成所述共中心点初至道集文件。

进一步地，根据测线的方位角，将各检波点和各炮点的第一坐标进行旋转，得到以所述测线为正东或正北的新坐标系中的第二坐标，具体包括：

将第一坐标按照顺时针方向进行旋转；

得到所述第二坐标为：

X_R＝(X-X₀)*cos(90-θ)+(Y-Y₀)*sin(90-θ)+X₀

Y_R＝(Y-Y₀)*cos(90-θ)-(X-X₀)*sin(90-θ)+Y₀