一种浅层地震单支反射数据自动解译方法和系统的制作方法
【技术领域】
[0001]本发明属于工程物探领域，涉及工程物探浅层地震反射勘探数据分析领域，尤其涉及一种浅层地震单支反射数据自动解译方法。主要应用于水电工程、水利工程、石油勘探、公路、铁路工程勘察等领域。
【背景技术】
[0002]浅层地震反射勘探是工程物探常用的方法，经济高效，它包括单支反射、反射多次覆盖等，其中单支反射是各类工程简单场地覆盖层厚度和地质分层探测的有效物探手段之一，也是测试各反射界面上覆地层平均速度参数的重要手段。
[0003]工程实际应用中，由于浅地表地质性质的横向不均匀变化，在缺乏钻孔提供速度参数情况下，通常需要在不同地质单元使用浅层地震单支反射来进行覆盖层厚度和地质分层探测、为地震反射多次覆盖数所处理和解译提供地层物性(速度)分层标志和地层速度参数。目前，浅层地震单支反射数据解译多通过手工分析来提供参考性的地层参数，效率极低，数据解译的精度不高。
【发明内容】
[0004]本发明的目的在于提供一种浅层地震单支反射数据自动解译方法，有效地提高工作效率，同时提高浅层地震单支反射数据解译结果的精度，从而解决现有技术中存在的前述问题。
[0005]为了实现上述目的，本发明采用的技术方案如下:
[0006]一种浅层地震单支反射数据自动解译方法，包括以下步骤:
[0007]SI，数据预处理，具体为对单支反射地震数据进行初步处理；
[0008]S2，拾取各反射波走时，具体为拾取来自各反射界面的反射波走时；
[0009]S3，绘制 X2 — T2 曲线；
[0010]S4，数据解译，具体为解译单支地震反射数据；
[0011]S5，成果及输出，具体为生成解译结果并输出结果。
[0012]优选的，步骤SI包括以下步骤:
[0013]S11，数据格式转换，即将地震仪器记录的原始地震数据转换为所需要的RPS格式数据；
[0014]S12，编辑文件头，即对格式转换后RPS格式数据中的文件头信息进行编辑，所述文件头信息包括:测线号、激发粧号、接收粧号。
[0015]优选的，步骤S2具体为，在波形图中，根据来自各反射界面反射波的时距关系，通过波组对比分析，确定来自不同反射界面的反射波，并拾取反射波走时。
[0016]优选的，步骤S3具体为，先分别计算各接收点到激发点距离的平方值记为X2、各接收点来自不同反射界面的反射波走时的平方值记为T2，然后以X2为横坐标，以T 2为纵坐标绘制成X2—T2曲线。
[0017]优选的，步骤S4包括以下步骤:
[0018]S41，解译上覆地层平均波速，即选择不同反射界面的(X2，T2)坐标点分别求取各反射界面上覆地层的平均地震纵波速度；
[0019]S42，解译反射界面埋深，即根据反射波的几何传播规律解译各测点不同反射界面的埋深。
[0020]优选的，步骤S5包括以下步骤:
[0021]S51，绘制剖面图，即将根据解译得到的各测点的不同反射界面埋深按规定的比例绘制工程剖面图；
[0022]S52，标注地层平均速度，即在工程剖面图上标注各地层的平均地震波速度信息；
[0023]S53，成果输出，即在时间剖面和工程剖面中进行选择，然后输出。
[0024]优选的，反射波走时与X2 — T2曲线之间的动态关联使得修改反射波走时后及时动态更新X2—T2曲线。
[0025]优选的，拾取或修改反射波走时可及时动态绘制该反射界面的时距曲线，通过各反射界面时距曲线的双曲线形态特征检验反射走时拾取的合理性和正确性。
[0026]一种浅层地震单支反射数据自动解译系统，包括:
[0027]数据预处理模块，用于对单支反射地震数据进行初步处理；
[0028]各反射波走时拾取模块，用于拾取来自各反射界面的反射波走时；
[0029]X2-T2曲线绘制模块，用于绘制X 2— T2曲线；
[0030]数据解译模块，用于解译单支地震反射数据；
[0031]成果及输出模块，用于生成解译结果并输出结果。
[0032]本发明的有益效果是:
[0033]本发明主要应用于水电工程、水利工程、石油勘探、公路、铁路工程勘察等领域。通过本发明的浅层地震单支反射数据自动解译方法和系统，有效地提高工作效率，同时提高浅层地震单支反射数据解译结果的精度。
【附图说明】
[0034]图1是本发明实施例的流程示意图；
[0035]图2是本发明实施例中拾取各反射界面反射波走时和绘制时距曲线图；
[0036]图3是本发明实施例中的X2 — T2曲线、各反射界面上覆地层平均速度图；
[0037]图4是本发明实施例的解译成果图。
【具体实施方式】
[0038]为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图，对本发明进行进一步详细说明。应当理解，此处所描述的【具体实施方式】仅仅用以解释本发明，并不用于限定本发明。
[0039]如附图所示，本发明公开了一种浅层地震单支反射数据自动解译方法和系统，该方法包括以下步骤:
[0040]1、数据预处理A，对单支反射地震数据进行初步处理；
[0041]2、拾取各界面反射波走时B，拾取来自各反射界面的反射波走时；
[0042]3、绘制X2 — T2曲线C，绘制X 2—T2曲线；
[0043]4、数据解译D，解译单支地震反射数据；
[0044]5、成果及输出E，生成解译结果并输出结果。
[0045]数据预处理A包括数据格式转换Al、编辑文件头A2两个子过程。其中数据格式转换Al是将地震仪器记录的原始地震数据转换为“物探测试数据综合分析系统”特定的RPS格式数据；编辑文件头A2是对格式转换后RPS格式数据中的测线号、激发粧号、接收粧号等文件头信息进行编辑。
[0046]拾取反射波走时B是在波形图中，根据来自各反射界面反射波的时距关系，通过波组对比分析，确定来自不同反射界面的反射波，并拾取反射波走时。
[0047]绘制X2—T2曲线C，是先分别计算X2 (各接收点到激发点距离的平方值)、T2 (各接收点来自不同反射界的反射波走时的平方值)的值，然后以为X2横坐标，以T2为纵坐标绘制成X2 — T2曲线。
[0048]数据解译D包括解译上覆地层平均波速Dl和解译反射界面埋深D2两个子过程，其中解译上覆地层平均波速Dl是合理选择不同反射界面的(X2，T2)坐标点分别求取各反射界面上覆地层的平均地震纵波速度；解译反射界面埋深D2是根据反射波的几何传播规律解译各测点不同反射界面的埋深。当地震单支反射方法用于简单场地勘探时，解译上覆地层平均波速Dl是数据解译D的一个关键过程和步骤，目的是解译各反射点埋深；当地震单支反射方法用于提供地震多次覆盖数据处理和解译的参数时，解译上覆地层平均波速Dl也是数据解译D的目的之一。
[0049]成果及输出E包括绘制剖面图Ε1、标注地层平均速度Ε2和成果输出Ε3三个子过程，其中绘制剖面图El是将根据解译得到的各测点的不同反射界面埋深按规定的比例绘制工程剖面图；标注地层平均速度Ε2是在工程剖面图上标注各地层的平均地震波速度信息；成果输出Ε3是在时间剖面和工程剖面中进行选择，然后输出。
[0050]反射波走时与X2 — T2曲线之间的动态关联使得修改反射波走时后及时动态更新X2-T2曲线；在X 2—T2曲线上选择不同反射界面的(X2、