专利名称：基于空间网格化算法的含逆断层等值线绘图方法
技术领域：
本发明涉及一种基于空间网格化算法的含逆断层等值线绘图方法，属于地震勘探 资料处理与解释领域。
背景技术：
等值线图是在石油勘探开发、采矿、地质、地球物理.地球化学和气象等领域等工 程和技术领域内应用极广的一种图形，是众多领域成果表示的重要图件之一。所谓绘制等 值线图就是在二维平面上，把大量离散的、又具有一定规律的几何量值或物理量值，用数学 的方法插值并将具有相同量值的点转换成图形的过程。等值线图是最主要也是最常用的一种图件，国内外许多学者对此进行了大量的 研究，但目前各种方法获得的等值线图还是常常令人难以满意。这一方面是由于已知信息 量有限，且分布不均。另一方面由于研究区常常存在许多的断层特别是逆断层，使得构造 非常复杂。这都使得计算机自动绘制等值线图变得很困难，例如，林建民在《西安矿业学院 学报》1991年第1期中公开了“三维地质数据的网格值生成和等值线绘图系统”。采用现有 的技术，当研究区断层较复杂时特别是存在逆断层时，常常绘制不出等值线图，或者绘制出 的等值线图效果并不理想，在含有逆断层中并不适用。文献《三维地质数据的网格值生成和 等值线绘图系统》在网格数据是平面数据(每一个网格数据点上只有一个值)时很有效，但 是针对数据含有断层和逆断层时(每一个网格数据点上可能有多个值)，该技术不支持。
发明内容
本发明的目的在于克服现有技术中对断层较复杂时特别是存在逆断层时绘制等 值线图存在的上述问题，提供一种基于空间网格化算法的含逆断层等值线绘图方法，本发 明采用空间网格化算法，以断层多边形和逆断层层次关系作为约束条件，对层位数据进行 网格化；对空间网格化数据进行追踪，实现含断层和逆断层的等值线绘图。为实现上述目的，本发明采用的技术方案如下
一种基于空间网格化算法的含逆断层等值线绘图方法，其特征在于，包括如下步骤
a、含断层和逆断层的空间网格化
①待插值点的初始化根据待插值点和断层多边形的包含关系，获得该待插值点有多 少层数据，并把每一层的数据和逆断层的上盘、下盘关联起来；
②网格化寻点寻找和待插值点在同一层的已知点；
③插值通过步骤①和②找出可以计算任一待插值点高度值的已知点，把这些待插值 点作为克里金算法的输入；
b、含断层和逆断层的等直线追踪
①对网格边缘进行标记根据输入的待插值点的高度值，对所有网格边缘作等值线可 以经过该网格的标记；
②连线在等直线追踪的时候通过标记对网格边缘进行连线；③终止绘图处理终止于断层。所述a步骤①中，判断待插值点是否被某逆断层群中的逆断层的断层多边形包 含，若不被任何的逆断层多边形包含，则关联的断层置为空。所述a步骤① 中，判断待插值点是否被某逆断层群中的逆断层的断层多边形包 含，若被逆断层多边形包含，则需要和相关的断层关联起来。所述关联时，按层次由上到下记录包含该待插值点的所有断层，并先把所有层次 上的值和关联的断层均置为空。所述关联时，对于第i个包含待插值点的逆断层，寻找所在逆断层群中包含待插 值点第i-Ι和i之间的逆断层，取出这些逆断层中层次上离i逆断层最近和最远的逆断层， 计算平面上待插值点和这两个断层之间的最近距离是上盘还是下盘，若待插值点距离最远 断层的部分是上盘，那么第i层和第i个包含待插值点的断层的上盘相关；若待插值点距离 最近断层的部分是下盘，那么第i层和第i个包含待插值点的断层的下盘相关。所述关联时，对于第i个包含待插值点的逆断层，寻找所在逆断层群中包含待插 值点第i-Ι逆断层和i逆断层之间的逆断层，若i-Ι逆断层和i逆断层之间没有其他逆断 层，则第i层数据和包含待插值点的i逆断层的上盘相关。所述关联时，对于第i个包含待插值点的逆断层，寻找所在逆断层群中包含待插 值点第i和i+ι之间的逆断层，取出这些逆断层中层次上离i最近和最远的逆断层，计算平 面上待插值点和这两个断层之间的最近距离是上盘还是下盘，若待插值点距离最远断层的 部分是上盘，那么第i层和第i个包含待插值点的断层的下盘相关；若待插值点距离最近断 层的部分是下盘，那么第i层和第i个包含待插值点的断层的上盘相关。所述关联时，对于第i个包含待插值点的逆断层，寻找所在逆断层群中包含待插 值点第i和i+ι之间的逆断层，若i和i+ι之间没有其他逆断层，则第i+ι层数据和包含待 插值点的第i个逆断层的下盘相关。所述b步骤②中，当等值线从前一网格进入，需要从当前网格寻找出口，在当前除 入口边之外其余网格边寻找是否有标记边，若有标记边则以该标记边作为该网格的出口。所述b步骤③中，当等值线到某一网格没有找到出口，并且该网格是区域的内部 网格，则等值线在断层结束。所述断层多边形是指对于断层数据按照一系列的点来存储，把这些点在平面上按 照顺序连接起来获得的多边形。所述逆断层上层的数据称为逆断层的上盘，逆断层下层的数据称为断层的下盘。采用本发明的优点在于
一、本发明采用空间网格化算法，以断层多边形和逆断层层次关系作为约束条件，对层 位数据进行网格化；对空间网格化数据进行追踪，实现含断层和逆断层的等值线绘图。二、本发明与文献《三维地质数据的网格值生成和等值线绘图系统》相比，本发明 可适用于数据含有断层和逆断层时的等值线绘图，且精度极高，有效支持了有断层和逆断 层时的等值线绘图。
图1为本发明含有断层和逆断层的网格化2为本发明点与断层关联图 图3为网格化寻点图
图4为数据中没有断层和逆断层，待插值点与原始数据点在同一层的结构示意图 图5为本发明一层数据中没有断层和逆断层与多层数据中有断层或者逆断层，待插值 点在一层，已知点在多层的示意图
图6为本发明数据中有断层或者逆断层，待插值和已知点都在多层的示意图 图7为本发明一个网格的示意图 图8为本发明等值线追踪的两种结束方式图。
具体实施例方式一种基于空间网格化算法的含逆断层等值线绘图方法，包括如下步骤
a、含断层和逆断层的空间网格化
①待插值点的初始化根据待插值点和断层多边形的包含关系，获得该待插值点有多 少层数据，并把每一层的数据和逆断层的上盘、下盘关联起来；
②网格化寻点寻找和待插值点在同一层的已知点；
③插值通过步骤①和②找出可以计算任一待插值点高度值的已知点，把这些待插值 点作为克里金算法的输入；
b、含断层和逆断层的等直线追踪
①对网格边缘进行标记根据输入的待插值点的高度值，对所有网格边缘作等值线可 以经过该网格的标记；
②连线在等直线追踪的时候通过标记对网格边缘进行连线；
③终止绘图处理终止于断层。以下对本发明作展开说明 含断层和逆断层的空间网格化算法
含有断层的网格化算法可以实现含有正断层和逆断层的数据的网格化。为了说明含有 断层的网格化算法，先对如下术语进行解释
待插值点需要插值的点，一般平面上的坐标(x，y)已知，高度ζ需要计算； 已知点已知平面坐标(χ，y)和高度ζ ；
断层多边形对于断层数据的存储，一般是按照一系列的点来进行存储，把这些点在平 面上按照顺序连起来，就可以获得一个多边形，这就是断层多边形，图1 (a)中看到的多边 形Fl和F2就是断层多边形；
逆断层的上盘和下盘逆断层的存在使得逆断层内的数据具有层次性，把逆断层上层 的数据称为逆断层的上盘，下层的数据称为断层的下盘。图1 (a)中F2是一个逆断层，断 层多边形的实现部分表示逆断层的上盘，虚线部分表示逆断层的下盘。在一般的网格化算法中，对平面上任意的一个进行插值时，在其周围8个方位上 寻找一直的点来计算待插值点的值。但是当平面含有断层的时候，根据以上方法找到的点 并不能全部都作为计算待插值点值的已知点。一个已知点可以作为计算待插值点高度值的 条件是两个点都在一块连续的区域。如图1所示(a)是含断层的平面，(b)是虚线上的剖 面图。图1中，b是需要插值的待插值点，a、c、d都是已知点，Fl是正断层，F2逆断层。对于正断层，已知点在正断层断层多边形内或待插值点和已知点的连线和断层多边形有两个 交点，这样的已知点不能作为用来计算待插值点值的已知点；对于逆断层，C的上层点可以 作为计算b高度值的已知点，但是C的下层点则不满足条件，点d也不满足条件。为了能够判断两个点是否在同一块连续的区域(中间没有断层断开)，则把点和断 层关联起来。点和断层关联如图1 (a)中c，c的上层点和F2的上盘在同一层，c的上层 点就和逆断层F2上盘关联，同样对于c的下层点和F2的下盘相关联。并且可以知道，一个 点的一层数据最多和两个断层相关联。根据关联关系，定义算法已知点输入格式 X Y Layerl Layer2 Layer3."
X是点在平面上的X值；Y是点在平面上的Y值；Layer指一层数据，包括高度值和关联 关系，每一个Layer格式如下
Z faultagel parti faultage2 part2
Z是指高度值，faultagel指关联的断层1的断层名，parti指的是关联的断层1的上 盘或下盘，faultage2指关联的断层2的断层名，part2指的是关联的断层2的上盘或下盘。 当和逆断层的上盘相关联part值为“2”，当和断层的下盘相关联part值为“1”。一个点的 一层数据最多可以和两个断层关联，也可以关联一个断层，也可以不和断层关联。当不和断 层相关联，faultage设为“null ”，part设为“_1”。图2是点和断层关联的例子，其中，a点 χ y zl Fl 2 null -1 z2 Fl 1 F2 2 z3 F2 1 F3 1；b点x y zl F3 2 null -1 z2 F3 1 F2 1 ；c 点X y zl F3 2 null -1 ；没有和逆断层关联的点x y ζ null -1 null -1。逆断层网格化包括三个步骤①待插值点的初始化；②网格化寻点；③克里金插 值，其中①、②是含有逆断层网格化的重点。①对待插值点进行初始化，根据点和断层多边 形的包含关系，获得该点有多少层数据，并把每一层的数据和逆断层的上下盘关联起来。② 寻找和待插值点在同一层的已知点。待插值点初始化
这一步主要实现对待插值点进行初始化根据点和断层多边形的包含关心，获得该点 有多少层数据，并把每一层的数据和逆断层的上下盘关联起来。首先，说明逆断层群逆断层群反映的是平面上逆断层从平面上看到的逆断层多 边形的相交关系。设表示逆断层群包含的所有断层，是所有断层的集合。断层群通过以下 方式获得
①任意选一逆断层Fl，把Fl加入；
②在中的所有逆断层中，若中的某断层和中的任一逆断层相交，将该逆断层加入；(的 意思是F至S)
③循环上述步骤，当不再变化时停止循环；
上述步骤获得了包含Fl的逆断层群，在剩下的断层中重复上述步骤找出所有逆断层 群。逆断层群反映的是逆断层之间的相交关系，逆断层群内的所有逆断层都直接或间接的 有相交关系。获得了所有的逆断层群，还需要获得逆断层群内逆断层之间的层次关系。断 层的层次关系是指断层的上下关系，如图2中断层Fl在断层F2之上。断层的层次关系在 加载测线数据时可以获得。对于平面上任意的待插值点具有如下特点①某一待插值点被η个断层包含，那么这个点对应有η+ 1层；
②对于第一层一定只和一断层的上盘相关，对于最下层一定只和一断层的上盘相关， 对于中间的点，必同时和一断层的上盘及另一断层的下盘相关；
③相关断层的位置可以是包含的断层，或被包含断层层次最近的断层，或上/下一包 含断层到本包含断层最远的不相关断层。其关系如表1所示。 表1断层包含关系
关联断层位置距离包含断层最远距离包含断层最近包含断层之上上盘相关下盘相关
④对于奇数层，必与某一断层的上盘相关，对于偶数层必与某一断层的下盘相关。利用上述特点，对待插值点进行初始化
①判断待插值点是否被某逆断层群中的逆断层的断层多边形包含，若不被任何的逆断 层多边形包含，那么该点只有一层，关联的断层置为空。②若被逆断层多边形包含，那么该待插值点包含多层数据，需要和相关的断层关 联起来。按层次由上到下记录包含该待插值点的所有断层，并先把所有层次上的值和关联 的断层都置为空。对于第i个包含待插值点的逆断层，寻找所在逆断层群中包含待插值点 第i-1和i之间的逆断层(若i-Ι和i之间没有其他逆断层，则第i层数据和包含待插值点 的第i个逆断层的上盘相关)，取出这些逆断层中层次上离i最近和最远的逆断层(与第i 层数据关联的逆断层只能是最近和/或最远)，计算平面上待插值点和这两个断层之间的最 近距离是上盘还是下盘。若待插值点距离最远断层的部分是上盘，那么第i层和第i个包 含待插值点的断层的上盘相关；若待插值点距离最近断层的部分是下盘，那么第i层和第i 个包含待插值点的断层的下盘相关。寻找所在逆断层群中包含待插值点第i和i+ι之间的 逆断层(若i和i+ι之间没有其他逆断层，则第i+ι层数据和包含待插值点的第i个逆断层 的下盘相关)，取出这些逆断层中层次上离i最近和最远的逆断层，计算平面上待插值点和 这两个断层之间的最近距离是上盘还是下盘。若待插值点距离最远断层的部分是上盘，那 么第i层和第i个包含待插值点的断层的下盘相关；若待插值点距离最近断层的部分是下 盘，那么第i层和第i个包含待插值点的断层的上盘相关。③需要对步骤②关联各层断层进行修正若除最底层数据外的其他层的数据(设 为第i层)没有和某逆断层的上盘相关联，将第i层数据和包含待插值点的i个逆断层的上 盘关联；若除最顶层数据外的其他层的数据(设为第i层)没有和某逆断层的下盘相关联， 将第i层数据和包含待插值点的i-Ι个逆断层的下盘关联；若第i层数据同时和同一逆断 层的同一上盘或下盘关联，若关联的都是下盘，取其中之一关联包含待插值点第i个逆断 层的上盘，若关联的都是下盘，取其中之一关联包含待插值点第i-ι个逆断层的下盘。网格化寻点
网格化寻点是对于任意的待插值点，在8个方位角上搜寻可以用来计算其值的已知点 (这类点称为候选点)。搜索方法和一般网格化方法相同，但是含有断层网格化搜索到的点 必须满足和已知点在同一块连续的区域，这也是含有断层网格化算法需要解决的主要问题 之一。对于每个输入的已知点，已经知道和断层的关联关系，而对于需要插值的未知点，通 过上一步的初始化，也可以和断层关联起来，利用这些关系，可以确定怎样的点可以用来作 为计算待插值点的高度值。可以根据以下特点来判断待插值点和已知点是否在同一块区域①判断是否在同一层对于被逆断层包含的点，可以和断层关联起来。如图3中，b5对 应点，被断层包含分为两层，上层点断层的上盘关联，下盘点和断层的下盘相关联；同理a6 对应点和b5对应点的情况是一样的。在判断是否同一层的时候，根据上面的层次关联关 系，b5对应点的上层和a6对应点的上层在同一层，b5对应点的下层和a6对应点的下层在
同一层。②经过相关断层上面点关联到断层还有一点就是具体和断层的上盘或下盘相 关。根据这一点，对于点a，由于a没有被断层包含，所以a不能和任何断层相关联。现在找 到了候选点a2，对于a2的上盘点，与F2的上盘相关，a与a2的上盘点的连线没有经过a2 经过的断层，所以a2上盘点可以作为网格化a的有效点；对于a2的下盘点与F2的下盘点 相关联，但是a于a2的连线经过了 a2下盘点所关联的断层，所以a2的下盘点不能作为网 格化a的点。同理若是待插值点在断层内，则待插值点也可以关联到断层的上盘或下盘，若 是与找到的候选点的连线经过了待插值点的相关联的断层，则这是找到的点也不能用来网 格化该待插值点。③经过断层若待插值点和候选点的连线同时经过了一个断层的上盘或下盘，那 么该候选点也就不能用来网格化该待插值点。利用以上特点来判断待插值点和已知点的数据是否在同一块连续的区域。
对于待插值点某一层的数据，找到其周围某个方向上的已知点
①如图4，待插值点和已知点都只有一层判断两点之间的连线是否经过断层(正断层 和逆断层)，若经过断层则舍弃该已知点，否则可以作为候选点。②如图5，待插值点只有一层，已知点多层从已知点的上层到下层依次进行判 断，若已知点某一层满足待插值点和已知点在平面上的连线没有经过已知点该层关联的断 层和任何一个正断层，那么可以作为候选点，否则舍弃该层点，继续下一层。③待插值点多层，已知点只有一层，对待插值点其中一层插值若待插值点和已知 点连线在平面上的连线没有经过待插值点该层关联的断层和任一正断层相关，那么该已知 点可以作为候选点，否则舍弃该点。④如图6，待插值多层，已知点多层，对待插值点其中一层插值从已知点的上层 到下层依次进行判断，若待插值点该层关联的断层和已知点关联的断层有同一个断层的同 一盘，那么可以把已知点该层数据作为候选点，否则继续判断。若两点在平面上的连线没有 经过待插值点该层关联的断层，没有经过已知点关联的断层和任何正断层，该层已知点也 可以作为候选点，否则舍弃已知点该层的点。克里金插值
通过以上步骤可以找出可以计算任一待插值点高度值的已知点，把这些点作为克里金 算法的输入，即可以获得待插值点的高度值。当待插值落在正断层的断层多边形区域内时， 不需要对该点进行插值，根据需要，可以把其置为便于处理的填充值。含断层和逆断层的等值线追踪
含有断层的等值线追踪需要同样需要处理正断层和逆断层的问题。一般等值线追踪 算法有两种情况的等值线起始于区域的边缘终止于区域的边缘；起始点和终止点都在区 域内部且是同一点。但是含有断层的等值线的追踪由于断层的等值线，等值线会终止于断 层，同时由于逆断层使得数据是具有多层的，应当保证一条连续的等值线在同一块连续的区域。实现含有断层的等值线追踪算法包含三个步骤①对网格边缘进行标记；②连线；③ 处理终止于断层；
对网格边缘进行标记
根据输入的高度值，对所有网格边缘进行标记。当某一网格的边缘被标记，说明等值线 可以经过该网格被标记的边缘，所以在追踪的时候只要识别标记进行连线即可。这里首先 定义一个符号B (A;B，C) :A、B、C指的是三个数值，B (A ;B，C)值A值介于B、C之间。对网格的边缘进行标记需要遍历所有的网格边缘，如图3 - 12所示黑点平面上显 示的网格化的点，称为网格端点；黑点之间的连线是网格的边缘；虚线是等值线，根据标记 穿越于网格之间。图7是一个网格，cell是当前网格标号，E，S，W，N是网格的四条边，需要对当前网 格的W和S边进行标记，并且需要记录标记是和端点的哪层数据定义的。定义记录标记层 次的数据格式如下
struct flag
{
Wl ； W2; Si;
S2; }；
其中Wl，W2表示对W边标记(ix, iy)和(ix, iy+1)的层次；Si，S2表示对S边标记(ix, iy)和(ix+1，iy)的层次
取相邻的两个端点，对其相连网格边缘进行标记，可分为三种情况进行讨论，假设需要 追踪的等值线的高程值是Z，并进行标记的边都是W边
①初始化所有标记，将所有网格标记的Wl，W2，Si，S2都设为一1 ；
②两个端点都是单层点其高度值分别是Zl、Z2，若两点都不被正断层的断层多边形 包含且B (Z;Z1，Z2)，那么标记该网格边缘。设置相应网格的Wl = 0，W2 = 0 ；
③一点是单层，一点是多层单层点只能和多层点的最上层和最下层在同一层。取多层 点的最上层，判断两点在平面上的连线是否和多层点最上层关联的断层相交，若不相交则 在同一层，并且B (Z ；Z1，Z2)，标记该条边，设置相应网格的Wl = 0，W2 = 0 ；若不能和多层 点的最上层点标记，多层点的最下层点，判断平面上两点的连线是否和最下层点关联的断 层相交，若不相交并且B (Z ；Zl, Z2)，标记该条边，设置相应网格的Wl = 0，W2 =多层点的 层次；
④两点都只多层若两点中有同一层关联同一个断层的上盘或下盘，并且B(Z;Z1, Z2)，那么把相应网格的Wl，W2设为相应的层次值。连线
当等值线从前一网格进入，需要从当前网格寻找出口，在当前除入口边之外其余网格 边寻找是否有标记边，若有标记边则以该标记边作为该网格的出口。在等值线追踪过程中， 需要记录每个点的层次，图7中网格等值线的出口边是E边，则该点的层次就就是对E边进 行标记时点(ix+l，iy)的层次。
处理终止于断层
当等值线到某一网格没有找到出口，并且该网格是区域的内部网格，就认为等值线将 在断层结束。等值线结束与断层有两种方式一种是延长线方式，另一种是结束于网格。
延长线方式是取当前追踪的最后两点，作延长线，寻找延长线和断层多边形的交 点，把该交点作为结束点。若交点处于逆断层的下盘，设置该点为下层点(等值线中点的层 次用于等值线的绘制，实际上有两种层次上层和下层，下层点的值可以设置随意非第一层 的值)，否则设置为上层点。一般情况下当某一内部网格找不到出口时，等值线要在断层结束，即当前等值线 最后点靠近断层并且该断层就是等值线结束的断层。但是有些情况下延长线方式得到的等 值线结束点和等值线最后点两点之间的价格太远，这样的结束点就不满足条件。在实现中， 把两个网格对角线的长度作为等值线最后点和结束点之间的距离限制，即当等值线最后点 和结束点之间的距离超过两个网格对角线的长度时，这种结束方式得到的结束点不满足条 件，需要利用另一种结束方式就是结束于网格。结束于网格方式首先找等值线最后点所在网格边对面网格边的两个端点，如图7 中若最后点所在网格边是E边，那么找W边的两个端点就是需要寻找的端点，分别计算等值 线最后点与寻找到两点的连线段与断层多边形的交点(有其中的一个连线段与断层没有交 点)，取其中的一个交点作为等值线结束点。结束点层次的获得和第一种方式一样。图8展示了等值线结束于断层的两种方式。其中(a)是延长线结束方式，Al和A2 是等值线最后点的两种情况。若Al是等值线的最后点B是延长线和断层的交点，这种情况 下B可以作为等值线的结束点；若A2是等值线的最后点C是延长线和断层的交点，这种情 况下A2C的长度超出了限制，所以C不能作为等值线结束点。图8 (b)是网格结束方式，这 里等值线最后点和网格的两个端点有两个交点B、C，B、C都可以作为等值线的结束点。
权利要求
1.一种基于空间网格化算法的含逆断层等值线绘图方法，其特征在于，包括如下步骤a、含断层和逆断层的空间网格化①待插值点的初始化根据待插值点和断层多边形的包含关系，获得该待插值点有多 少层数据，并把每一层的数据和逆断层的上盘、下盘关联起来；②网格化寻点寻找和待插值点在同一层的已知点；③插值通过步骤①和②找出可以计算任一待插值点高度值的已知点，把这些待插值 点作为克里金算法的输入；b、含断层和逆断层的等直线追踪①对网格边缘进行标记根据输入的待插值点的高度值，对所有网格边缘作等值线可 以经过该网格的标记；②连线在等直线追踪的时候通过标记对网格边缘进行连线；③终止绘图处理终止于断层。
2.根据权利要求1所述的基于空间网格化算法的含逆断层等值线绘图方法，其特征在 于所述a步骤①中，判断待插值点是否被某逆断层群中的逆断层的断层多边形包含，若不 被任何的逆断层多边形包含，则关联的断层置为空。
3.根据权利要求2所述的基于空间网格化算法的含逆断层等值线绘图方法，其特征在 于所述a步骤①中，判断待插值点是否被某逆断层群中的逆断层的断层多边形包含，若被 逆断层多边形包含，则需要和相关的断层关联起来。
4.根据权利要求3所述的基于空间网格化算法的含逆断层等值线绘图方法，其特征在 于所述关联时，按层次由上到下记录包含该待插值点的所有断层，并先把所有层次上的值 和关联的断层均置为空。
5.根据权利要求4所述的基于空间网格化算法的含逆断层等值线绘图方法，其特征在 于所述关联时，对于第i个包含待插值点的逆断层，寻找所在逆断层群中包含待插值点第 i-1和i之间的逆断层，取出这些逆断层中层次上离i逆断层最近和最远的逆断层，计算平 面上待插值点和这两个断层之间的最近距离是上盘还是下盘，若待插值点距离最远断层的 部分是上盘，那么第i层和第i个包含待插值点的断层的上盘相关；若待插值点距离最近断 层的部分是下盘，那么第i层和第i个包含待插值点的断层的下盘相关。
6.根据权利要求5所述的基于空间网格化算法的含逆断层等值线绘图方法，其特征在 于所述关联时，对于第i个包含待插值点的逆断层，寻找所在逆断层群中包含待插值点第 i-1逆断层和i逆断层之间的逆断层，若i-i逆断层和i逆断层之间没有其他逆断层，则第 i层数据和包含待插值点的i逆断层的上盘相关。
7.根据权利要求5或6所述的基于空间网格化算法的含逆断层等值线绘图方法，其特 征在于所述关联时，对于第i个包含待插值点的逆断层，寻找所在逆断层群中包含待插值 点第i和i+Ι之间的逆断层，取出这些逆断层中层次上离i最近和最远的逆断层，计算平面 上待插值点和这两个断层之间的最近距离是上盘还是下盘，若待插值点距离最远断层的部 分是上盘，那么第i层和第i个包含待插值点的断层的下盘相关；若待插值点距离最近断层 的部分是下盘，那么第i层和第i个包含待插值点的断层的上盘相关。
8.根据权利要求7所述的基于空间网格化算法的含逆断层等值线绘图方法，其特征在于所述关联时，对于第i个包含待插值点的逆断层，寻找所在逆断层群中包含待插值点第 i和i+ι之间的逆断层，若i和i+ι之间没有其他逆断层，则第i+ι层数据和包含待插值点 的第i个逆断层的下盘相关。
9.根据权利要求1一6或8中任一所述的基于空间网格化算法的含逆断层等值线绘图 方法，其特征在于所述b步骤②中，当等值线从前一网格进入，需要从当前网格寻找出口， 在当前除入口边之外其余网格边寻找是否有标记边，若有标记边则以该标记边作为该网格 的出口。
10.根据权利要求9中任一所述的基于空间网格化算法的含逆断层等值线绘图方法， 其特征在于所述b步骤③中，当等值线到某一网格没有找到出口，并且该网格是区域的内 部网格，则等值线在断层结束。
全文摘要
本发明公开了一种基于空间网格化算法的含逆断层等值线绘图方法，包括如下步骤含断层和逆断层的空间网格化待插值点的初始化、网格化寻点、插值；含断层和逆断层的等直线追踪对网格边缘进行标记、连线、绘图处理终止于断层。本发明采用空间网格化算法，以断层多边形和逆断层层次关系作为约束条件，对层位数据进行网格化；对空间网格化数据进行追踪，实现含断层和逆断层的等值线绘图。
文档编号G01V1/30GK102096106SQ20101058290
公开日2011年6月15日 申请日期2010年12月10日 优先权日2010年12月10日
发明者刘鸿, 吕文彪, 姚兴苗, 巫骏, 邹文, 陈三平 申请人:中国石油集团川庆钻探工程有限公司