[发明专利]构建地层速度场的方法及装置

说明书

技术领域

本申请涉及地球物理勘探领域，尤其涉及一种构建地下地层速度场的方法。

背景技术

现有的地震资料解释中大多都采用常规变速成图反映地质构造，尽管这种方法也能反映出地下地层构造的部分形态，但是，难以完全反映出地层构造的形态、大小以及构造高点的真实位置，特别是，当研究区块遭受过不同时期构造运动的影响，底层不同程度的受到挤压、褶皱、剥蚀时，这种岩性的横向相变造成地下介质存在较强的各向异性，使得地震波的传播速度在纵、横方向上变化均较大。如果继续采用现有的常规变速成图方法进行深度-时间域(Z，t)转换，往往会导致错误的结果。

发明内容

本发明的目的在于提供一种构建地层速度场的方法和装置，以便科学地反映出地下地层的构造形态。

根据本发明的一方面，提供一种构建地层速度场的方法，所述方法包括：获取目标区域内沿每条地表测线剖面上的每道速度序列的每个种子点处的地层速度；按照预设的网格步长对所述目标区域的地下地层进行空间划分，以生成该地下地层的六面体网格体系；基于获取的每道速度序列的每个种子点处的地层速度，确定每道速度序列在穿越生成的六面体网格体系的每个网格平面时的地层速度；基于确定的每道速度序列在穿越生成的六面体网格体系的每个网格平面时的地层速度，确定生成的六面体网格体系的每个网格平面上的每个网格顶点处的地层速度，以便构建所述目标区域的地层速度场。

优选地，所述确定每道速度序列在穿越生成的六面体网格体系的每个网格平面时的地层速度的具体步骤为：采用线性插值法，确定每道速度序列在穿越生成的六面体网格体系的每个网格平面时的地层速度。

优选地，当每道速度序列垂直穿越每个网格平面时，所述线性插值法的计算公式为：

wk=wj(zk-zj)+wi(zi-zk)zi-zj]]>

其中，z_k表示生成的六面体网格体系的第k个网格平面的高度；z_i、z_j分别表示在同一地表测线剖面上的同一道速度序列中距离第k个网格平面最近的两个种子点的高度，并且z_i>z_k>z_j；w_i表示z_i处的地层速度；w_j表示z_j处的地层速度；w_k表示该道速度序列在穿越第k个网格平面时的地层速度。

优选地，所述确定生成的六面体网格体系的每个网格平面上的每个网格顶点处的地层速度的具体步骤为：采用距离加权反比插值法，确定生成的六面体网格体系的每个网格平面上的每个网格顶点处的地层速度。

优选地，所述距离加权反比插值法的计算公式为：