CN102096106A - 基于空间网格化算法的含逆断层等值线绘图方法 - Google Patents
基于空间网格化算法的含逆断层等值线绘图方法 Download PDFInfo
- Publication number
- CN102096106A CN102096106A CN 201010582907 CN201010582907A CN102096106A CN 102096106 A CN102096106 A CN 102096106A CN 201010582907 CN201010582907 CN 201010582907 CN 201010582907 A CN201010582907 A CN 201010582907A CN 102096106 A CN102096106 A CN 102096106A
- Authority
- CN
- China
- Prior art keywords
- trap
- tomography
- point
- interpolation
- interpolation point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Abstract
本发明公开了一种基于空间网格化算法的含逆断层等值线绘图方法,包括如下步骤:含断层和逆断层的空间网格化:待插值点的初始化、网格化寻点、插值;含断层和逆断层的等直线追踪:对网格边缘进行标记、连线、绘图处理终止于断层。本发明采用空间网格化算法,以断层多边形和逆断层层次关系作为约束条件,对层位数据进行网格化;对空间网格化数据进行追踪,实现含断层和逆断层的等值线绘图。
Description
技术领域
本发明涉及一种基于空间网格化算法的含逆断层等值线绘图方法,属于地震勘探资料处理与解释领域。
背景技术
等值线图是在石油勘探开发、采矿、地质、地球物理.地球化学和气象等领域等工程和技术领域内应用极广的一种图形,是众多领域成果表示的重要图件之一。所谓绘制等值线图就是在二维平面上,把大量离散的、又具有一定规律的几何量值或物理量值,用数学的方法插值并将具有相同量值的点转换成图形的过程。
等值线图是最主要也是最常用的一种图件, 国内外许多学者对此进行了大量的研究,但目前各种方法获得的等值线图还是常常令人难以满意。 这一方面是由于已知信息量有限, 且分布不均。另一方面由于研究区常常存在许多的断层特别是逆断层,使得构造非常复杂。这都使得计算机自动绘制等值线图变得很困难,例如,林建民在《西安矿业学院学报》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个逆断层的下盘相关。
所述b步骤②中,当等值线从前一网格进入,需要从当前网格寻找出口,在当前除入口边之外其余网格边寻找是否有标记边,若有标记边则以该标记边作为该网格的出口。
所述b步骤③中,当等值线到某一网格没有找到出口,并且该网格是区域的内部网格,则等值线在断层结束。
所述断层多边形是指对于断层数据按照一系列的点来存储,把这些点在平面上按照顺序连接起来获得的多边形。
所述逆断层上层的数据称为逆断层的上盘,逆断层下层的数据称为断层的下盘。
采用本发明的优点在于:
一、本发明采用空间网格化算法,以断层多边形和逆断层层次关系作为约束条件,对层位数据进行网格化;对空间网格化数据进行追踪,实现含断层和逆断层的等值线绘图。
二、本发明与文献《三维地质数据的网格值生成和等值线绘图系统》相比,本发明可适用于数据含有断层和逆断层时的等值线绘图,且精度极高,有效支持了有断层和逆断层时的等值线绘图。
附图说明
图1为本发明含有断层和逆断层的网格化图
图2为本发明点与断层关联图
图3为网格化寻点图
图4为数据中没有断层和逆断层,待插值点与原始数据点在同一层的结构示意图
图5为本发明一层数据中没有断层和逆断层与多层数据中有断层或者逆断层,待插值点在一层,已知点在多层的示意图
图6为本发明数据中有断层或者逆断层,待插值和已知点都在多层的示意图
图7为本发明一个网格的示意图
图8为本发明等值线追踪的两种结束方式图 。
具体实施方式
一种基于空间网格化算法的含逆断层等值线绘图方法,包括如下步骤:
a、含断层和逆断层的空间网格化:
①待插值点的初始化:根据待插值点和断层多边形的包含关系,获得该待插值点有多少层数据,并把每一层的数据和逆断层的上盘、下盘关联起来;
②网格化寻点:寻找和待插值点在同一层的已知点;
③插值:通过步骤①和②找出可以计算任一待插值点高度值的已知点,把这些待插值点作为克里金算法的输入;
b、含断层和逆断层的等直线追踪:
①对网格边缘进行标记:根据输入的待插值点的高度值,对所有网格边缘作等值线可以经过该网格的标记;
②连线:在等直线追踪的时候通过标记对网格边缘进行连线;
③终止:绘图处理终止于断层。
以下对本发明作展开说明:
含断层和逆断层的空间网格化算法
含有断层的网格化算法可以实现含有正断层和逆断层的数据的网格化。为了说明含有断层的网格化算法,先对如下术语进行解释:
待插值点:需要插值的点,一般平面上的坐标(x,y)已知,高度z需要计算;
已知点:已知平面坐标(x,y)和高度z;
断层多边形:对于断层数据的存储,一般是按照一系列的点来进行存储,把这些点在平面上按照顺序连起来,就可以获得一个多边形,这就是断层多边形,图1(a)中看到的多边形F1和F2就是断层多边形;
逆断层的上盘和下盘:逆断层的存在使得逆断层内的数据具有层次性,把逆断层上层的数据称为逆断层的上盘,下层的数据称为断层的下盘。图1(a)中F2是一个逆断层,断层多边形的实现部分表示逆断层的上盘,虚线部分表示逆断层的下盘。
在一般的网格化算法中,对平面上任意的一个进行插值时,在其周围8个方位上寻找一直的点来计算待插值点的值。但是当平面含有断层的时候,根据以上方法找到的点并不能全部都作为计算待插值点值的已知点。一个已知点可以作为计算待插值点高度值的条件是两个点都在一块连续的区域。如图1所示:(a)是含断层的平面,(b)是虚线上的剖面图。图1中,b是需要插值的待插值点,a、c、d都是已知点,F1是正断层,F2逆断层。对于正断层,已知点在正断层断层多边形内或待插值点和已知点的连线和断层多边形有两个交点,这样的已知点不能作为用来计算待插值点值的已知点;对于逆断层,c的上层点可以作为计算b高度值的已知点,但是c的下层点则不满足条件,点d也不满足条件。
为了能够判断两个点是否在同一块连续的区域(中间没有断层断开),则把点和断层关联起来。点和断层关联:如图1(a)中c,c的上层点和F2的上盘在同一层,c的上层点就和逆断层F2上盘关联,同样对于c的下层点和F2的下盘相关联。并且可以知道,一个点的一层数据最多和两个断层相关联。
根据关联关系,定义算法已知点输入格式:
X Y Layer1 Layer2 Layer3…
X是点在平面上的X值;Y是点在平面上的Y值;Layer指一层数据,包括高度值和关联关系,每一个Layer格式如下:
Z faultage1 part1 faultage2 part2
Z是指高度值,faultage1指关联的断层1的断层名,part1指的是关联的断层1的上盘或下盘,faultage2指关联的断层2的断层名,part2指的是关联的断层2的上盘或下盘。当和逆断层的上盘相关联part值为“2”,当和断层的下盘相关联part值为“1”。一个点的一层数据最多可以和两个断层关联,也可以关联一个断层,也可以不和断层关联。当不和断层相关联,faultage设为“null”,part设为“-1”。图2是点和断层关联的例子,其中,a点:x y z1 F1 2 null -1 z2 F1 1 F2 2 z3 F2 1 F3 1;b点:x y z1 F3 2 null -1 z2 F3 1 F2 1;c点:x y z1 F3 2 null -1;没有和逆断层关联的点:x y z null -1 null -1。
逆断层网格化包括三个步骤:①待插值点的初始化;②网格化寻点;③克里金插值,其中①、②是含有逆断层网格化的重点。①对待插值点进行初始化,根据点和断层多边形的包含关系,获得该点有多少层数据,并把每一层的数据和逆断层的上下盘关联起来。②寻找和待插值点在同一层的已知点。
待插值点初始化:
这一步主要实现对待插值点进行初始化:根据点和断层多边形的包含关心,获得该点有多少层数据,并把每一层的数据和逆断层的上下盘关联起来。
首先,说明逆断层群:逆断层群反映的是平面上逆断层从平面上看到的逆断层多边形的相交关系。设表示逆断层群包含的所有断层,是所有断层的集合。断层群通过以下方式获得:
①任意选一逆断层F1,把F1加入;
②在中的所有逆断层中,若中的某断层和中的任一逆断层相交,将该逆断层加入;(的意思是F至S)
③循环上述步骤,当不再变化时停止循环;
上述步骤获得了包含F1的逆断层群,在剩下的断层中重复上述步骤找出所有逆断层群。逆断层群反映的是逆断层之间的相交关系,逆断层群内的所有逆断层都直接或间接的有相交关系。获得了所有的逆断层群,还需要获得逆断层群内逆断层之间的层次关系。断层的层次关系是指断层的上下关系,如图2中断层F1在断层F2之上。断层的层次关系在加载测线数据时可以获得。
对于平面上任意的待插值点具有如下特点:
①某一待插值点被n个断层包含,那么这个点对应有n+1层;
②对于第一层一定只和一断层的上盘相关,对于最下层一定只和一断层的上盘相关,对于中间的点,必同时和一断层的上盘及另一断层的下盘相关;
③相关断层的位置:可以是包含的断层,或被包含断层层次最近的断层,或上/下一包含断层到本包含断层最远的不相关断层。其关系如表1所示。
表1 断层包含关系
关联断层位置 | 距离包含断层最远 | 距离包含断层最近 |
包含断层之上 | 上盘相关 | 下盘相关 |
④对于奇数层,必与某一断层的上盘相关,对于偶数层必与某一断层的下盘相关。
利用上述特点,对待插值点进行初始化:
①判断待插值点是否被某逆断层群中的逆断层的断层多边形包含,若不被任何的逆断层多边形包含,那么该点只有一层,关联的断层置为空。
②若被逆断层多边形包含,那么该待插值点包含多层数据,需要和相关的断层关联起来。按层次由上到下记录包含该待插值点的所有断层,并先把所有层次上的值和关联的断层都置为空。对于第i个包含待插值点的逆断层,寻找所在逆断层群中包含待插值点第i-1和i之间的逆断层(若i-1和i之间没有其他逆断层,则第i层数据和包含待插值点的第i个逆断层的上盘相关),取出这些逆断层中层次上离i最近和最远的逆断层(与第i层数据关联的逆断层只能是最近和/或最远),计算平面上待插值点和这两个断层之间的最近距离是上盘还是下盘。若待插值点距离最远断层的部分是上盘,那么第i层和第i个包含待插值点的断层的上盘相关;若待插值点距离最近断层的部分是下盘,那么第i层和第i个包含待插值点的断层的下盘相关。寻找所在逆断层群中包含待插值点第i和i+1之间的逆断层(若i和i+1之间没有其他逆断层,则第i+1层数据和包含待插值点的第i个逆断层的下盘相关),取出这些逆断层中层次上离i最近和最远的逆断层,计算平面上待插值点和这两个断层之间的最近距离是上盘还是下盘。若待插值点距离最远断层的部分是上盘,那么第i层和第i个包含待插值点的断层的下盘相关;若待插值点距离最近断层的部分是下盘,那么第i层和第i个包含待插值点的断层的上盘相关。
③需要对步骤②关联各层断层进行修正:若除最底层数据外的其他层的数据(设为第i层)没有和某逆断层的上盘相关联,将第i层数据和包含待插值点的i个逆断层的上盘关联;若除最顶层数据外的其他层的数据(设为第i层)没有和某逆断层的下盘相关联,将第i层数据和包含待插值点的i-1个逆断层的下盘关联;若第i层数据同时和同一逆断层的同一上盘或下盘关联,若关联的都是下盘,取其中之一关联包含待插值点第i个逆断层的上盘,若关联的都是下盘,取其中之一关联包含待插值点第i-1个逆断层的下盘。
网格化寻点:
网格化寻点是对于任意的待插值点,在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
{
W1;
W2;
S1;
S2;
};
其中W1,W2表示对W边标记(ix,iy)和(ix,iy+1)的层次;S1,S2表示对S边标记(ix,iy)和(ix+1,iy)的层次
取相邻的两个端点,对其相连网格边缘进行标记,可分为三种情况进行讨论,假设需要追踪的等值线的高程值是Z,并进行标记的边都是W边:
①初始化所有标记,将所有网格标记的W1,W2,S1,S2都设为-1;
②两个端点都是单层点:其高度值分别是Z1、Z2,若两点都不被正断层的断层多边形包含且B(Z;Z1,Z2),那么标记该网格边缘。设置相应网格的W1=0,W2=0;
③一点是单层,一点是多层:单层点只能和多层点的最上层和最下层在同一层。取多层点的最上层,判断两点在平面上的连线是否和多层点最上层关联的断层相交,若不相交则在同一层,并且B(Z;Z1,Z2),标记该条边,设置相应网格的W1=0,W2=0;若不能和多层点的最上层点标记,多层点的最下层点,判断平面上两点的连线是否和最下层点关联的断层相交,若不相交并且B(Z;Z1,Z2),标记该条边,设置相应网格的W1=0,W2=多层点的层次;
④两点都只多层:若两点中有同一层关联同一个断层的上盘或下盘,并且B(Z;Z1,Z2),那么把相应网格的W1,W2设为相应的层次值。
连线:
当等值线从前一网格进入,需要从当前网格寻找出口,在当前除入口边之外其余网格边寻找是否有标记边,若有标记边则以该标记边作为该网格的出口。在等值线追踪过程中,需要记录每个点的层次,图7中网格等值线的出口边是E边,则该点的层次就就是对E边进行标记时点(ix+1,iy)的层次。
处理终止于断层:
当等值线到某一网格没有找到出口,并且该网格是区域的内部网格,就认为等值线将在断层结束。等值线结束与断层有两种方式:一种是延长线方式,另一种是结束于网格。
延长线方式是取当前追踪的最后两点,作延长线,寻找延长线和断层多边形的交点,把该交点作为结束点。若交点处于逆断层的下盘,设置该点为下层点(等值线中点的层次用于等值线的绘制,实际上有两种层次:上层和下层,下层点的值可以设置随意非第一层的值),否则设置为上层点。
一般情况下当某一内部网格找不到出口时,等值线要在断层结束,即当前等值线最后点靠近断层并且该断层就是等值线结束的断层。但是有些情况下延长线方式得到的等值线结束点和等值线最后点两点之间的价格太远,这样的结束点就不满足条件。在实现中,把两个网格对角线的长度作为等值线最后点和结束点之间的距离限制,即当等值线最后点和结束点之间的距离超过两个网格对角线的长度时,这种结束方式得到的结束点不满足条件,需要利用另一种结束方式就是结束于网格。
结束于网格方式首先找等值线最后点所在网格边对面网格边的两个端点,如图7中若最后点所在网格边是E边,那么找W边的两个端点就是需要寻找的端点,分别计算等值线最后点与寻找到两点的连线段与断层多边形的交点(有其中的一个连线段与断层没有交点),取其中的一个交点作为等值线结束点。结束点层次的获得和第一种方式一样。
图8展示了等值线结束于断层的两种方式。其中(a)是延长线结束方式,A1和A2是等值线最后点的两种情况。若A1是等值线的最后点B是延长线和断层的交点,这种情况下B可以作为等值线的结束点;若A2是等值线的最后点C是延长线和断层的交点,这种情况下A2C的长度超出了限制,所以C不能作为等值线结束点。图8(b)是网格结束方式,这里等值线最后点和网格的两个端点有两个交点B、C,B、C都可以作为等值线的结束点。
Claims (10)
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-1逆断层和i逆断层之间没有其他逆断层,则第i层数据和包含待插值点的i逆断层的上盘相关。
7.根据权利要求5或6所述的基于空间网格化算法的含逆断层等值线绘图方法,其特征在于:所述关联时,对于第i个包含待插值点的逆断层,寻找所在逆断层群中包含待插值点第i和i+1之间的逆断层,取出这些逆断层中层次上离i最近和最远的逆断层,计算平面上待插值点和这两个断层之间的最近距离是上盘还是下盘,若待插值点距离最远断层的部分是上盘,那么第i层和第i个包含待插值点的断层的下盘相关;若待插值点距离最近断层的部分是下盘,那么第i层和第i个包含待插值点的断层的上盘相关。
8.根据权利要求7所述的基于空间网格化算法的含逆断层等值线绘图方法,其特征在于:所述关联时,对于第i个包含待插值点的逆断层,寻找所在逆断层群中包含待插值点第i和i+1之间的逆断层,若i和i+1之间没有其他逆断层,则第i+1层数据和包含待插值点的第i个逆断层的下盘相关。
9.根据权利要求1—6或8中任一所述的基于空间网格化算法的含逆断层等值线绘图方法,其特征在于:所述b步骤②中,当等值线从前一网格进入,需要从当前网格寻找出口,在当前除入口边之外其余网格边寻找是否有标记边,若有标记边则以该标记边作为该网格的出口。
10.根据权利要求9中任一所述的基于空间网格化算法的含逆断层等值线绘图方法,其特征在于:所述b步骤③中,当等值线到某一网格没有找到出口,并且该网格是区域的内部网格,则等值线在断层结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010582907 CN102096106B (zh) | 2010-12-10 | 2010-12-10 | 基于空间网格化算法的含逆断层等值线绘图方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010582907 CN102096106B (zh) | 2010-12-10 | 2010-12-10 | 基于空间网格化算法的含逆断层等值线绘图方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102096106A true CN102096106A (zh) | 2011-06-15 |
CN102096106B CN102096106B (zh) | 2013-05-01 |
Family
ID=44129274
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010582907 Expired - Fee Related CN102096106B (zh) | 2010-12-10 | 2010-12-10 | 基于空间网格化算法的含逆断层等值线绘图方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102096106B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102243680A (zh) * | 2011-07-21 | 2011-11-16 | 中国科学技术大学 | 一种网格划分方法及系统 |
CN102385067A (zh) * | 2011-08-03 | 2012-03-21 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 含逆断层等值线绘图方法 |
CN102622514A (zh) * | 2012-02-17 | 2012-08-01 | 中国科学院地理科学与资源研究所 | 一种海洋温度场数据处理方法 |
CN103424773A (zh) * | 2012-05-25 | 2013-12-04 | 中国石油化工股份有限公司 | 一种基于层位拉平法的古地貌恢复方法 |
CN104570757A (zh) * | 2013-10-10 | 2015-04-29 | 罗伯特·博世有限公司 | 用于运行集成控制组件的方法和装置 |
CN104635262A (zh) * | 2013-11-13 | 2015-05-20 | 中国石油天然气集团公司 | 一种基于增强型矩形网格的正逆断层等值线自动生成方法 |
CN105139452A (zh) * | 2015-09-01 | 2015-12-09 | 电子科技大学 | 一种基于图像分割的地质曲线重建方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105528522A (zh) * | 2015-12-15 | 2016-04-27 | 中国石油天然气股份有限公司 | 一种基于空间网格计算连续性油气藏资源量的方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6191787B1 (en) * | 1998-02-10 | 2001-02-20 | Schlumberger Technology Corporation | Interactively constructing, editing, rendering and manipulating geoscience models |
CN101582173A (zh) * | 2009-06-24 | 2009-11-18 | 中国石油集团川庆钻探工程有限公司 | 复杂地质构造块状模型构建方法 |
-
2010
- 2010-12-10 CN CN 201010582907 patent/CN102096106B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6191787B1 (en) * | 1998-02-10 | 2001-02-20 | Schlumberger Technology Corporation | Interactively constructing, editing, rendering and manipulating geoscience models |
CN101582173A (zh) * | 2009-06-24 | 2009-11-18 | 中国石油集团川庆钻探工程有限公司 | 复杂地质构造块状模型构建方法 |
Non-Patent Citations (1)
Title |
---|
《中国优秀硕士学位论文全文数据库信息科技辑(月刊)》 20090915 黄洲 地质构造等值线的插值和绘制方法研究 31-61 1-10 , 第9期 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102243680A (zh) * | 2011-07-21 | 2011-11-16 | 中国科学技术大学 | 一种网格划分方法及系统 |
CN102243680B (zh) * | 2011-07-21 | 2013-01-16 | 中国科学技术大学 | 一种网格划分方法及系统 |
CN102385067A (zh) * | 2011-08-03 | 2012-03-21 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 含逆断层等值线绘图方法 |
CN102385067B (zh) * | 2011-08-03 | 2013-09-11 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 含逆断层等值线绘图方法 |
CN102622514A (zh) * | 2012-02-17 | 2012-08-01 | 中国科学院地理科学与资源研究所 | 一种海洋温度场数据处理方法 |
CN103424773A (zh) * | 2012-05-25 | 2013-12-04 | 中国石油化工股份有限公司 | 一种基于层位拉平法的古地貌恢复方法 |
CN103424773B (zh) * | 2012-05-25 | 2016-02-03 | 中国石油化工股份有限公司 | 一种基于层位拉平法的古地貌恢复方法 |
CN104570757A (zh) * | 2013-10-10 | 2015-04-29 | 罗伯特·博世有限公司 | 用于运行集成控制组件的方法和装置 |
CN104570757B (zh) * | 2013-10-10 | 2020-03-17 | 罗伯特·博世有限公司 | 用于运行集成控制组件的方法和装置 |
CN104635262A (zh) * | 2013-11-13 | 2015-05-20 | 中国石油天然气集团公司 | 一种基于增强型矩形网格的正逆断层等值线自动生成方法 |
CN104635262B (zh) * | 2013-11-13 | 2017-06-06 | 中国石油天然气集团公司 | 一种基于增强型矩形网格的正逆断层等值线自动生成方法 |
CN105139452A (zh) * | 2015-09-01 | 2015-12-09 | 电子科技大学 | 一种基于图像分割的地质曲线重建方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102096106B (zh) | 2013-05-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102096106B (zh) | 基于空间网格化算法的含逆断层等值线绘图方法 | |
CN109447317B (zh) | 一种基于多人同时拣货的仓库最短路径优化方法 | |
US7801924B2 (en) | Decision tree construction via frequent predictive itemsets and best attribute splits | |
US9279905B2 (en) | Method for extracting a thumbnail image from a training image so as to constrain the multipoint geostatistical modeling of the subsoil | |
US8566069B2 (en) | Method for geologically modeling seismic data by trace correlation | |
US20100185692A1 (en) | System and method for determining intervals of a space filling curve in a query box | |
CN102385067B (zh) | 含逆断层等值线绘图方法 | |
CN101980054A (zh) | 一种在高密度地震静校正处理中建立近地表速度模型的方法 | |
CN103226209A (zh) | 具有波状航线的多船地震采集 | |
CN111047700B (zh) | 一种基于断层产生时序的地质构造三维建模方法及系统 | |
CN107783185B (zh) | 一种层析静校正的处理方法及装置 | |
CN109001813A (zh) | 一种压制多次波的方法、装置及系统 | |
CN112797999B (zh) | 一种多无人艇协同遍历路径规划方法及系统 | |
CN112859852A (zh) | 一种路径规划方法和装置 | |
CN107316341A (zh) | 一种多点地质统计学沉积相建模方法 | |
AU2009223731B2 (en) | Method for determing distinct alternative paths between two object sets in 2-D and 3-D heterogeneous data | |
CN113449910B (zh) | 一种基于顺序存储二叉树的航线自动生成方法 | |
Hwang et al. | Sequencing picking operations and travel time models for man-on-board storage and retrieval warehousing system | |
CN110097636B (zh) | 一种基于可视域分析的选址规划方法 | |
Elfick | Contouring by use of a triangular mesh | |
CN102253411A (zh) | 含逆断层的复杂地质构造三维地震解释层位插值的方法 | |
CN113985899B (zh) | 基于区间多目标优化的水下机器人全局路径规划方法 | |
CN102353987B (zh) | 含逆断层的复杂地质构造三维地震解释数据管理组织方法 | |
CN114610027A (zh) | 一种船舶航海路径规划方法 | |
US11454734B2 (en) | Automated extraction of horizon patches from seismic data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130501 Termination date: 20181210 |
|
CF01 | Termination of patent right due to non-payment of annual fee |