速度场构建方法技术

提供一种速度场构建方法，所述方法包括：(A)确定目标工区，并获取所述目标工区的测井数据；(B)根据所述目标工区的测井数据获取所述目标工区内的多个散点的速度值；(C)将所述多个散点作为三角网的顶点来基于所述多个散点的速度值构建所述目标工区的三角网；(D)选取多个待插值点，基于构建的三角网来计算所述多个待插值点的速度值；(E)根据所述多个待插值点的速度值来构建所述目标工区的速度场。根据上述速度场构建方法能够有效提高速度场构建的效率和准确性。

【技术实现步骤摘要】
速度场构建方法
本专利技术总体来说涉及地球物理勘探
，更具体地讲，涉及一种盆地级工区的速度场构建方法。
技术介绍
目前，现有的地震资料解释中大多都采用常规变速成图，以适当反映出地下地层构造形态，尽管这种方法能够反映出地下地层构造的部分形态，但是，难以完全落实目的层构造形态、大小及构造高点的真实位置，特别是运用常规变速成图方法进行时深转换时，误差较大。在现有的地震资料解释中除常规变速成图方法之外，还存在平均速度场构建方法，但随着地震勘探领域和程度的不断扩大，地震勘探的工区也随之扩大，现有的构建速度场方法无法构建大工区(例如，四川盆地18×104km2)的速度场。
技术实现思路
本专利技术的目的是提出一种适用于盆地级工区的速度场构建方法，以解决现有的构建速度场方法无法构建大工区速度场的技术问题。根据本专利技术示例性实施例的一方面，提供一种速度场构建方法，所述方法包括：(A)确定目标工区，并获取所述目标工区的测井数据；(B)根据所述目标工区的测井数据获取所述目标工区内的多个散点的速度值；(C)将所述多个散点作为三角网的顶点来基于所述多个散点的速度值构建所述目标工区的三角网；(D)选取多个待插值点，基于构建的三角网来计算所述多个待插值点的速度值；(E)根据所述多个待插值点的速度值来构建所述目标工区的速度场。可选地，在步骤(D)中，基于构建的三角网来计算所述多个待插值点中的任一待插值点的速度值的步骤可包括：(D1)确定所述任一待插值点所在的剖面层；(D2)将构建的三角网中的所有顶点投影到所述剖面层上，得到各顶点在所述剖面层上对应的投影点；(D3)确定在所述剖面层上各投影点形成的三角网格中包围所述任一待插值点的最小三角形；(D4)计算所述最小三角形的各顶点的速度值；(D5)根据所述最小三角形的各顶点的速度值计算所述任一待插值点的速度值。可选地，在步骤(D3)之前可还包括：确定所述任一待插值点是否在各投影点形成的三角网格的范围内；如果所述任一待插值点不在各投影点形成的三角网格的范围内，则重新构建三角网，并返回执行步骤(D2)；如果所述任一待插值点在各投影点形成的三角网格的范围内，则执行步骤(D3)。可选地，重新构建三角网的步骤可包括：在所述任一待插值点附近的预定范围内以预定步长为间隔构建矩形网络；确定在所述剖面层上，与所述任一待插值点距离最近的投影点；从所述矩形网络中任意选取一个网格节点，确定与选取的网格节点距离最近的预定个数投影点；计算所述预定个数的投影点的速度值；根据所述预定个数的投影点的速度值计算所述网格节点的速度值；根据所述网格节点的速度值计算构建的矩形网络中的所有网格节点的速度值；根据矩形网络中的所有网格节点的速度值和三角网的各顶点的速度值来重新构建三角网。可选地，计算所述预定个数中的任一投影点的速度值的步骤可包括：确定所述任一投影点所在一道深度道集；确定所述一道CDP上竖直向上方向上距离所述任一投影点最近的第一散点确定所述该一道CDP上竖直向下方向上距离所述任一投影点最近的第二散点根据第一散点的深度值和速度值，以及第二散点的深度值和速度值计算所述任一投影点的速度值。采用本专利技术示例性实施例的速度场构建方法，基于三角网进行空间速度构建，以准确求取目标工区中任一点的速度值，有效提高了速度场构建的效率和准确性。附图说明通过下面结合附图进行的详细描述，本专利技术的上述和其它目的、特点和优点将会变得更加清楚，其中：图1示出根据本专利技术的示例性实施例的速度场构建方法的流程图；图2示出根据本专利技术的示例性实施例的构建的三角网与选取的多个待插值点的空间分布示意图；图3示出根据本专利技术的示例性实施例的图1中的基于构建的三角网来计算多个待插值点的速度值的步骤的流程图；图4示出根据本专利技术的示例性实施例的图2中的重新构建三角网的步骤的流程图；图5示出根据本专利技术的示例性实施例的一条测线上的多个待插值点与构建的三角网的示意图；图6示出根据本专利技术的示例性实施例的重新构建三角网的示意图；图7示出根据本专利技术的示例性实施例的计算剖面层上的各投影点的速度值的步骤的流程图；图8示出根据本专利技术的示例性实施例的计算投影点的速度值的示意图；图9示出根据本专利技术的示例性实施例的在同一剖面层上投影点与待插值点的示意图。具体实施方式下面，将参照附图详细描述本专利技术的示例性实施例。提供参照附图的以下描述以帮助对由权利要求及其等同物限定的本专利技术的示例性实施例的全面理解。包括各种特定细节以帮助理解，但这些细节仅被视为是示例性的。因此，本领域的普通技术人员将认识到在不脱离本专利技术的范围和精神的情况下，可对描述于此的示例性实施例进行各种改变和修改。此外，为了清楚和简洁，省略对公知的功能和结构的描述。图1是示出根据本专利技术的示例性实施例的速度场构建方法的流程图。参照图1，在步骤S10中，确定目标工区，并获取所述目标工区的测井数据。作为示例，所述目标工区可为盆地级工区，这里，可利用现有的各种方法来获取该目标工区的测井数据。在步骤S20中，根据所述目标工区的测井数据获取所述目标工区内的多个散点的速度值。例如，可从所述目标工区内选取多个散点，然后在从该目标工区的测井数据中提取所述多个散点的速度值。在步骤S30中，将所述多个散点作为三角网的顶点来基于所述多个散点的速度值构建所述目标工区的三角网。这里，可利用现有的各种方法来基于多个散点的速度值构建目标工区的三角网，优选地，可利用Delaunay三角网剖分分治算法来基于多个散点的速度值构建目标工区的三角网。在步骤S40中，选取多个待插值点，基于构建的三角网来计算所述多个待插值点的速度值。这里，选取的多个待插值点可为目标工区内的任意点，基于构建的三角网可计算目标工区内的任意点的速度值。图2示出了构建的三角网与选取的多个待插值点的空间分布示意图。具体说来，针对所述多个待插值点中的每个待插值点，可通过以下步骤来计算每个待插值点的速度值：确定待插值点所在的剖面层，将构建的三角网的所有顶点投影到该剖面层上，根据三角网的所有顶点在该剖面层上对应的投影点的速度值来计算待插值点的速度值。应理解，待插值点所在的剖面层与目标工区的目的层不同，可根据经验来设定不同剖面层之间的距离。在步骤S50中，根据所述多个待插值点的速度值来构建所述目标工区的速度场。这里，可利用现有的各种方法根据多个待插值点的速度值构建速度场，本专利技术对此部分的内容不再赘述。采用本专利技术示例性实施例的速度场构建方法，基于Delaunay三角网进行空间速度构建，以求取目标工区内任意点的速度值，可有效提高速度场构建的效率和准确性。下面参照图3来详细描述步骤S40中基于构建的三角网来计算多个待插值点的速度值的步骤。图3示出根据本专利技术的示例性实施例的图1中的基于构建的三角网来计算多个待插值点的速度值的步骤的流程图。参照图3，在步骤S301中，从目标工区内选取多个待插值点，例如，可选取m个待插值点，m为大于零的自然数。在步骤S302中，确定m个待插值点中的第i个待插值点所在的剖面层。在步骤S303中，将构建的三角网中的所有顶点投影到所述剖面层上，得到各顶点在所述剖面层上对应的投影点。在步骤S304中，确定第i个待插值点是否在各投影点形成的三角网格的范围内。如果第i个待插值点不在各投影点形成的三角网格的范围内，本文档来自技高网...

【技术保护点】
一种速度场构建方法，所述方法包括：(A)确定目标工区，并获取所述目标工区的测井数据；(B)根据所述目标工区的测井数据获取所述目标工区内的多个散点的速度值；(C)将所述多个散点作为三角网的顶点来基于所述多个散点的速度值构建所述目标工区的三角网；(D)选取多个待插值点，基于构建的三角网来计算所述多个待插值点的速度值；(E)根据所述多个待插值点的速度值来构建所述目标工区的速度场。

【技术特征摘要】
1.一种速度场构建方法，所述方法包括：(A)确定目标工区，并获取所述目标工区的测井数据；(B)根据所述目标工区的测井数据获取所述目标工区内的多个散点的速度值；(C)将所述多个散点作为三角网的顶点来基于所述多个散点的速度值构建所述目标工区的三角网；(D)选取多个待插值点，基于构建的三角网来计算所述多个待插值点的速度值；(E)根据所述多个待插值点的速度值来构建所述目标工区的速度场，其中，在步骤(D)中，基于构建的三角网来计算所述多个待插值点中的任一待插值点的速度值的步骤包括：(D1)确定所述任一待插值点所在的剖面层；(D2)将构建的三角网中的所有顶点投影到所述剖面层上，得到各顶点在所述剖面层上对应的投影点；(D3)确定在所述剖面层上各投影点形成的三角网格中包围所述任一待插值点的最小三角形；(D4)计算所述最小三角形的各顶点的速度值；(D5)根据所述最小三角形的各顶点的速度值计算所述任一待插值点的速度值。2.根据权利要求1所述的方法，其中，在步骤(D3)之前还包括：确定所述任一待插值点是否在各投影点形成的三角网格的范围内；如果所述任一待插值点不在各投影点形成的三角网格...

【专利技术属性】
技术研发人员：陈小二，陶正喜，彭才，范昆，徐敏，范晓，刘开元，谭荣彪，巫盛洪，
申请(专利权)人：中国石油集团川庆钻探工程有限公司地球物理勘探公司，
类型：发明
国别省市：四川;51