一种由地震数据道头炮检点坐标重建SPS文件的方法技术

本发明专利技术涉及一种由地震数据道头炮检点坐标重建SPS文件的方法，属于地震资料数字处理技术领域。本发明专利技术包括如下步骤：由原始地震数据逐炮逐道提取道头中的炮点和检波点坐标信息；判断是否需要进行坐标旋转，若是，则进行坐标旋转；若否，则执行下一步骤；可视化炮检点坐标，分析判断观测系统的类型；根据观测系统的类型将炮检点所在区域网格化，并逐炮逐道计算炮检点所在网格点位置，赋予相应的线号、桩号，生成SPS文件。本发明专利技术解决了地震资料数字处理过程中SPS文件缺失的问题。过程中SPS文件缺失的问题。过程中SPS文件缺失的问题。

【技术实现步骤摘要】
一种由地震数据道头炮检点坐标重建SPS文件的方法


[0001]本专利技术涉及一种由地震数据道头炮检点坐标重建SPS文件的方法，属于地震资料数字处理


技术介绍

[0002]在地震资料数字处理的过程中，观测系统的定义非常重要，通常采用SPS标准格式的文件定义观测系统。SPS(Shell Processing Support Format for 3D Surveys)格式标准由英国SHELL公司制定，已被SEG推荐为地震资料处理辅助数据的标准格式。其作用在于将地震勘探野外采集的炮检点坐标以及有关的辅助数据，以一种通用的标准格式传送到处理中心。在2011年国家能源局发布的《地震勘探辅助数据SPS格式》中，已将其作为我国石油天然气行业标准在全国行业范围内执行。
[0003]目前，国内外实际地震资料采集均要求附带SPS文件，但实际地震资料处理中经常会遇到SPS文件缺失的情况。一方面，随着时移地震(4D地震)的发展，一些老工区的地震资料需要重新处理，这些数据有可能因为采集比较早而没有记录SPS文件，也有可能因为后期资料的整理与保存出现问题而导致SPS文件丢失。另一方面，在实际地震资料野外采集的过程中，当采集信号传输出现问题时，则无法从采集仪器上下载SPS文件。针对这些问题，我们需要用特定的方法恢复或者重建SPS文件。
[0004]SPS文件包含四个部分：其一，炮点(S)文件，定义炮点的线号、桩号、坐标等信息；其二，检波点(R)文件，定义检波点的线号、桩号、坐标等信息；其三，炮点和检波点的关系(X)文件，定义每个地震记录的文件号和对应的炮点信息以及与其对应的检波点信息，即观测系统。其四，格式(C)文件，定义炮点文件、检波点文件和炮检点关系文件的格式，通常用标准格式SPS Rev2.1即可。在SPS文件缺失的情况下，炮点和检波点的坐标信息可从地震数据道头中提取得到，往往是缺失炮点和检波点的线号和桩号信息，从而难以建立炮检点关系文件。传统的SPS文件恢复方法需要已有野外采集的班报或记录文件，并借助专门的文件整理软件或者观测系统设计软件，其过程繁琐，容易出错。
[0005]本专利技术提出一种由地震数据道头炮检点坐标重建SPS文件的方法，主要用于生成炮点文件、检波点文件和炮检点关系文件。该方法通过网格化炮点、检波点所在区域，通过逐炮逐道计算炮点、检波点所在网格点的位置以赋予相应的线号、桩号，并确定炮检点的关系，从而重建得到新的SPS文件，用于定义观测系统，以便后期的地震资料数字处理。

技术实现思路

[0006]针对上述
技术介绍
和思路，本专利技术提出了一种由地震数据道头炮检点坐标重建SPS文件的方法，以解决地震资料数字处理过程中SPS文件缺失的问题。
[0007]本专利技术是采用以下的技术方案实现的：本专利技术所述一种由地震数据道头炮检点坐标重建SPS文件的方法，包括如下步骤：步骤一：由原始地震数据逐炮逐道提取道头中的炮点和检波点坐标信息；
步骤二：判断是否需要进行坐标旋转，若是，则进行坐标旋转；若否，则执行步骤三；步骤三：可视化炮检点坐标，分析判断观测系统的类型；步骤四：根据观测系统的类型将炮检点所在区域网格化，并逐炮逐道计算炮检点所在网格点位置，赋予相应的线号、桩号，生成SPS文件；进一步地，步骤二中所述判断是否需要进行坐标旋转，可利用如下方法：当X方向最大值与最小值之差大于检波点线距的一半时，则进行坐标旋转；进一步地，步骤二中所述坐标旋转，可具体包括如下步骤：(1)计算测线方位角：
ꢀꢀꢀꢀꢀꢀ
(1)其中，a为测线方位角，arctan表示反正切函数，min_y0表示Y方向最小值，min_x0表示Y方向最小值min_y0对应的X坐标，max_x0表示X方向最大值，max_y0表示X方向最大值max_x0对应的Y坐标；(2)进行坐标旋转：
ꢀꢀꢀꢀꢀꢀꢀ
(2)其中，(x0,y0)表示坐标原点，即分别表示X方向最小和Y方向最小值，(x,y)表示旋转之前的坐标，(x
’
,y
’
)表示旋转之后的坐标，sin表示正弦函数，cos表示余弦函数。
[0008]进一步地，步骤三所述观测系统类型，包括如下两大类：(1)平行型观测系统，指炮线和检波点线平行，又可具体分为两种：
①
炮线和检波点线均平行于X轴；
②
炮线和检波点线均平行于Y轴；(2)正交型观测系统，指炮线和检波点线互相垂直，又可具体分为两种：
①
炮线平行于X轴，检波点线平行于Y轴；
②
炮线平行于Y轴，检波点线平行于X轴；进一步地，步骤四所述将炮检点所在区域进行网格化，逐炮逐道计算炮检点所在网格点位置，并赋予相应的线号和桩号，线号的给定根据步骤三中判定的观测系统类型具体有以下四种：(1)平行型观测系统，炮线和检波点线均平行于X轴：
ꢀꢀꢀꢀꢀ
(3)(2)平行型观测系统，炮线和检波点线均平行于Y轴：
ꢀꢀꢀꢀꢀꢀ
(4)(3)正交型观测系统，炮线平行于X轴，检波点线平行于Y轴：
ꢀꢀꢀꢀꢀ
(5)(4)正交型观测系统，炮线平行于Y轴，检波点线平行于X轴：
ꢀꢀꢀꢀꢀ
(6)以上各式中，source_line表示炮线号，sx表示炮点X坐标，sy表示炮点Y坐标，min_sx表示炮点X坐标最小值，min_sy表示炮点Y坐标最小值，dsx表示网格化炮点坐标时给定的沿X方向的网格大小，dsy表示网格化炮点坐标时给定的沿Y方向的网格大小；rec_line表示检波点线号，rx表示检波点X坐标，ry表示检波点Y坐标，min_rx表示检波点X坐标最小值，min_ry表示检波点Y坐标最小值，drx表示网格化检波点坐标时给定的沿X方向的网格大小dry表示网格化检波点坐标时给定的沿Y方向的网格大小；而桩号的给定则可以在一开始赋予起始的炮点桩号的检波点桩号，待计算到新的炮点或检波点时，将已有的炮点桩号和检波点桩号分别赋予当前炮点或检波点，再分别进行自增运算，以计算下一个新的炮点或检波点。需要注意的是，我们通过规则网格化方法赋予炮检点相应的线号和桩号，所以应保存落入每一个网格内的炮点和检波点，计算下一个炮点或检波点时，应与落入同一个网格内的炮点或检波点进行遍历比较，以确定是否新的炮点或检波点。
[0009]本专利技术的有益效果是：采用本专利技术所述的一种由地震数据道头炮检点坐标重建SPS文件的方法，仅利用地震数据道头炮检点坐标便可通过网格化的方法生成SPS文件，解决了实际地震数据处理过程中SPS文件缺失的问题。
附图说明
[0010]图1为本专利技术的流程图；图2为实际资料原始炮检点位置坐标分布图；图3为经过坐标旋转后的炮检点坐标图；图4为炮点文件说明图；图5为检波点文件说明图；图6为炮检点关系文件说明图；图7为共炮点道集示意图；图8为共检波点道集示意图。
具体实施方式
[0011]为了使本专利技术目的、技术方案更加清楚明白，下面结合实施例，对本专利技术作进一步详细说明。
...

【技术保护点】

【技术特征摘要】
1.一种由地震数据道头炮检点坐标重建SPS文件的方法，其特征在于，包括如下步骤：步骤一：由原始地震数据逐炮逐道提取道头中的炮点和检波点坐标信息；步骤二：判断是否需要进行坐标旋转，若是，则进行坐标旋转；若否，则执行步骤三；步骤三：可视化炮检点坐标，分析判断观测系统的类型；步骤四：根据观测系统的类型将炮检点所在区域网格化，并逐炮逐道计算炮检点所在网格点位置，赋予相应的线号、桩号，生成SPS文件。2.根据权利要求1所述的一种由地震数据道头炮检点坐标重建SPS文件的方法，其特征在于，所述步骤二中，判断是否需要进行坐标旋转，可利用如下方法：当X方向最大值与最小值之差大于检波点线距的一半时，则进行坐标旋转；进一步地，步骤二中所述坐标旋转，可具体包括如下步骤：(1)计算测线方位角：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)其中，a为测线方位角，arctan表示反正切函数，min_y0表示Y方向最小值，min_x0表示Y方向最小值min_y0对应的X坐标，max_x0表示X方向最大值，max_y0表示X方向最大值max_x0对应的Y坐标；(2)进行坐标旋转：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)其中，(x0,y0)表示坐标原点，即分别表示X方向最小值和Y方向最小值，(x,y)表示旋转之前的坐标，(x
’
,y
’
)表示旋转之后的坐标，sin表示正弦函数，cos表示余弦函数。3.根据权利要求1所述的一种由地震数据道头炮检点坐标重建SPS文件的方法，其特征在于，所述步骤三中，观测系统类型包括如下两大类：(1)平行型观测系统，指炮线和检波点线平行，又可具体分为两种：
①
炮线和检波点线均平行于X轴；
②
炮线和检波点线均平行于Y轴；(2)正交型观测系统，指炮线和检波点线互相垂直，又可具体分为两种：
①
炮线平行于X轴，检波点线平行于Y轴；
②
炮线平行于Y轴，检波...

【专利技术属性】
技术研发人员：何万利，张建中，杨华臣，
申请(专利权)人：中国海洋大学，
类型：发明
国别省市：