[发明专利]基于空间网格化算法的含逆断层等值线绘图方法

说明书

技术领域

本发明涉及一种基于空间网格化算法的含逆断层等值线绘图方法，属于地震勘探资料处理与解释领域。

背景技术

等值线图是在石油勘探开发、采矿、地质、地球物理．地球化学和气象等领域等工程和技术领域内应用极广的一种图形，是众多领域成果表示的重要图件之一。所谓绘制等值线图就是在二维平面上，把大量离散的、又具有一定规律的几何量值或物理量值，用数学的方法插值并将具有相同量值的点转换成图形的过程。

等值线图是最主要也是最常用的一种图件, 国内外许多学者对此进行了大量的研究，但目前各种方法获得的等值线图还是常常令人难以满意。 这一方面是由于已知信息量有限, 且分布不均。另一方面由于研究区常常存在许多的断层特别是逆断层，使得构造非常复杂。这都使得计算机自动绘制等值线图变得很困难，例如，林建民在《西安矿业学院学报》1991年第1期中公开了“三维地质数据的网格值生成和等值线绘图系统”。采用现有的技术，当研究区断层较复杂时特别是存在逆断层时，常常绘制不出等值线图，或者绘制出的等值线图效果并不理想，在含有逆断层中并不适用。文献《三维地质数据的网格值生成和等值线绘图系统》在网格数据是平面数据（每一个网格数据点上只有一个值）时很有效，但是针对数据含有断层和逆断层时（每一个网格数据点上可能有多个值），该技术不支持。

发明内容

本发明的目的在于克服现有技术中对断层较复杂时特别是存在逆断层时绘制等值线图存在的上述问题，提供一种基于空间网格化算法的含逆断层等值线绘图方法，本发明采用空间网格化算法，以断层多边形和逆断层层次关系作为约束条件，对层位数据进行网格化；对空间网格化数据进行追踪，实现含断层和逆断层的等值线绘图。

    为实现上述目的，本发明采用的技术方案如下：

    一种基于空间网格化算法的含逆断层等值线绘图方法，其特征在于，包括如下步骤：

    a、含断层和逆断层的空间网格化：

①待插值点的初始化：根据待插值点和断层多边形的包含关系，获得该待插值点有多少层数据，并把每一层的数据和逆断层的上盘、下盘关联起来；

②网格化寻点：寻找和待插值点在同一层的已知点；

③插值：通过步骤①和②找出可以计算任一待插值点高度值的已知点，把这些待插值点作为克里金算法的输入；

b、含断层和逆断层的等直线追踪：

①对网格边缘进行标记：根据输入的待插值点的高度值，对所有网格边缘作等值线可以经过该网格的标记；

②连线：在等直线追踪的时候通过标记对网格边缘进行连线；

③终止：绘图处理终止于断层。

所述a步骤①中，判断待插值点是否被某逆断层群中的逆断层的断层多边形包含，若不被任何的逆断层多边形包含，则关联的断层置为空。

所述a步骤①中，判断待插值点是否被某逆断层群中的逆断层的断层多边形包含，若被逆断层多边形包含，则需要和相关的断层关联起来。

所述关联时，按层次由上到下记录包含该待插值点的所有断层，并先把所有层次上的值和关联的断层均置为空。

所述关联时，对于第i个包含待插值点的逆断层，寻找所在逆断层群中包含待插值点第i-1和i之间的逆断层，取出这些逆断层中层次上离i逆断层最近和最远的逆断层，计算平面上待插值点和这两个断层之间的最近距离是上盘还是下盘，若待插值点距离最远断层的部分是上盘，那么第i层和第i个包含待插值点的断层的上盘相关；若待插值点距离最近断层的部分是下盘，那么第i层和第i个包含待插值点的断层的下盘相关。

所述关联时，对于第i个包含待插值点的逆断层，寻找所在逆断层群中包含待插值点第i-1逆断层和i逆断层之间的逆断层，若i-1逆断层和i逆断层之间没有其他逆断层，则第i层数据和包含待插值点的i逆断层的上盘相关。

所述关联时，对于第i个包含待插值点的逆断层，寻找所在逆断层群中包含待插值点第i和i+1之间的逆断层，取出这些逆断层中层次上离i最近和最远的逆断层，计算平面上待插值点和这两个断层之间的最近距离是上盘还是下盘，若待插值点距离最远断层的部分是上盘，那么第i层和第i个包含待插值点的断层的下盘相关；若待插值点距离最近断层的部分是下盘，那么第i层和第i个包含待插值点的断层的上盘相关。

所述关联时，对于第i个包含待插值点的逆断层，寻找所在逆断层群中包含待插值点第i和i+1之间的逆断层，若i和i+1之间没有其他逆断层，则第i+1层数据和包含待插值点的第i个逆断层的下盘相关。