[发明专利]基于三角网格的三维表层模型构建方法

说明书

技术领域

本发明涉及地球物理勘探技术，更具体地讲，涉及一种能够近乎准确反映近地表结构的基于三角网格的三维表层模型构建方法。

背景技术

目前，我国地震勘探的主战场主要集中在山地、沙漠、黄土塬等复杂地表区，所面临的地表条件异常复杂，不同区域有不同的表层特点，大量存在高速层顶界不平坦和表层速度横向变化大的情况，表层问题是地震勘探的首要问题。

地震勘探解释的理论都假定激发点与接收点是在一个水平面上，并且地层速度是均匀的。但实际上地面常常不平坦，各个激发点深度也可能不同，低速带中的波速与地层中的波速又相差悬殊，所以必将影响实测的时距曲线形状。为了消除这些影响，对原始地震数据要进行地形校正、激发深度校正、低速带校正等，这些校正对同一观测点的不同地震界面都是不变的，因此统称静校正。

另外，随着各油田勘探程度的提高，勘探目标向低幅度、小断块和岩性圈闭的勘探方向转移，勘探难度不断加大，精度要求不断提高。为了准确地分辨出地下复杂圈闭，静校正的问题日显突出。许多专家认为，静校正已经成为制约勘探的一个瓶颈，要想获得地震勘探的进一步突破，首先要解决静校正问题。

此外，由静校正引起的闭合问题是复杂的，它是由基准面、高速层顶界面、表层调查点、交点、表层模型、时深曲线、替换速度和计算方法等多种因素综合形成的，在静校正计算中，如果上述因素不可靠或使用不当，就会造成剖面闭合问题或构造形态失真问题。这是勘探中长期存在的问题，特别是对低幅度构造的勘探尤为突出。所以，必需努力解决静校正计算的基础控制问题，从而保证在控制点、交点上各种数据的受控和闭合，进而保证剖面的闭合。要完成上述工作，建立准确的表层模型是后续数据处理的基础。

利用建立近地表模型来计算静校正量的方法从原理上来讲是非常可靠和真实的，能够反映近地表地质的变化规律。只是如何求取到符合表层地质规律的模型是非常困难的工作，控制点数量的多少、对该地区沉积规律的了解，均影响表层点的解释成果。另外表层控制点处的模型即便可以真实的确定，但表层点间低降速界面形态，往往得不到准确的确认，这样就会使模型造成偏差，进而影响校正量精度和处理精度。

以往建模的方法多采用距离加权厚度线性内插法，然后将逐段进行调整厚度，这种工作不但繁琐，效率低，且不容易反应在整个面上的表层模型变化规律，受人为因素影响较大，所解释的结果带有不确定性。由于地震采集与处理技术的迅速发展，人们对如何准确求取表层结构的重视程度越来越高，这些均推动了建模技术的发展同时又对表层结构的建模准确性提出了更高的要求。

发明内容

本发明的目的在于提供一种能够提高建模准确性的基于三角网格的三维表层模型构建方法。

为了实现上述目的，提供一种基于三角网格的三维表层模型构建方法，所述方法包括如下步骤：对某区域的表层调查点进行测量，以直接或间接获得每个表层调查点的地理位置信息和表层模型参数；获取所述区域所有炮点检波点的地理位置信息；对由表层调查点的地理位置信息构成的二维点集进行三角剖分，以获得Delaunay三角网格；找到每个炮点检波点所在的三角形，并以该三角形的三个顶点位置的表层调查点作为控制点，并基于面积加权插值法计算每个炮点检波点的表层模型参数Vp，由此获得所述区域的三维表层模型。

在所述方法中，利用Bowyer-Watson算法对由表层调查点构成的二维点集进行三角剖分。

所述面积加权插值法的计算公式为：

Vp=Σn=13An*VnAt,]]>

其中，At表示由所三个控制点的地理位置构成的三角形的面积；

An表示被计算的炮点检波点的地理位置与所述三个控制点形成的第n个小三角形的面积；

Vn表示所述三角形的第n个控制点处的表层调查点的表层模型参数，