[发明专利]基于自动识别折射波起点的单支折射数据解译方法和系统

权利要求书

1.一种基于自动识别折射波起点的单支折射数据解译方法，其特征在于，首先①在时距曲线上选择离炮点最近的一个接收点(T₁，X₁)，将该点(T₁，X₁)作为假想的折射波起点，将该点以后的各接收点(T₂，X₂)、(T₃，X₃)...(T_i，X_i)均作为下卧目标层的折射波点；

②根据点(T₁，X₁)的时距关系求得表层的地震波有效速度V₁；

③用最小二乘法对点(T₂，X₂)以后的各点(T_i，X_i)进行线性拟合，求得以点(T₂，X₂)为起点的下卧目标层界面速度V₂；

④用最小二乘法对点(T₃，X₃)以后的各点(T_i，X_i)进行线性拟合，求得以点(T₃，X₃)为起点的目标层界面速度V₂＇；

⑤以V₂≥2V₁和|V₂－V₂＇|/V₂≤0.1为判断条件，若不能同时满足这两个条件，则自动将下一个点(T₂，X₂)作为假想的折射波起点重复上述过程，直至同时满足上述两个判断条件，此时查找到的点(T_k，X_k)即为折射波起点；

然后自动计算表层地震波有效速度和下卧目标层的界面速度，最后计算出各测点的表层厚度。

2.根据权利要求1所述的基于自动识别折射波起点的单支折射数据解译方法，其特征在于，具体包括以下步骤：

S1，拾取地震波初至走时，在地震波形数据中识别、拾取地震波初至走时；

S2，绘制单支或相遇时距曲线，建立分析计算所需的各测点地震波时距关系曲线图；

S3，自动识别时距曲线中的折射波起点，根据所建立的地震波时距关系曲线图自动识别和查找折射波起点；

S4，分析计算并输出解译成果，自动分析计算各测点表层厚度并输出最终成果。

3.根据权利要求2所述的基于自动识别折射波起点的单支折射数据解译方法，其特征在于，S1中所述拾取地震波初至走时，具体为：先将地震仪器记录的地震数据转换为RPS格式数据，然后编辑测线号、炮点桩号、接收桩号文件头信息，根据相邻记录道地震波初至的相关性进行识别，拾取地震波初至走时。

4.根据权利要求2所述的基于自动识别折射波起点的单支折射数据解译方法，其特征在于，S2中所述绘制单支或相遇时距曲线，具体为：选择纵横比例均为1:1000的纵横绘图比例，根据所获得的地震波数据文件中的地震波初至走时，以小桩号边炮点为坐标原点、以接收桩号为横坐标、以地震波初至走时为纵坐标绘制时距曲线。

5.根据权利要求2所述的基于自动识别折射波起点的单支折射数据解译方法，其特征在于，S3中所述自动识别时距曲线中的折射波起点，具体为：采用循环过程在时距曲线上自动查找折射波起点；

每个所述循环过程均包括三个子过程：

S31，自动分析计算表层地震波有效速度；

S32，用最小二乘法分析计算下卧目标层的界面速度；

S33，用条件判断所查找到的点是否为折射波起点；

其中S31子过程是根据所选点的时距关系(T，X)来求取表层地震波有效速度V₁，S32子过程是用最小二乘法线性拟合各折射波点(T，X)来求取目标层的界面速度V₂，S33子过程是用限制条件来判断所查找到的点是否为折射波起点，也是用来终止循环的判断条件。