CN1956011A - 不规则三维地质几何体的自动建模方法 - Google Patents
不规则三维地质几何体的自动建模方法 Download PDFInfo
- Publication number
- CN1956011A CN1956011A CNA2005101145649A CN200510114564A CN1956011A CN 1956011 A CN1956011 A CN 1956011A CN A2005101145649 A CNA2005101145649 A CN A2005101145649A CN 200510114564 A CN200510114564 A CN 200510114564A CN 1956011 A CN1956011 A CN 1956011A
- Authority
- CN
- China
- Prior art keywords
- block window
- mrow
- moving
- key point
- section
- 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
- 238000000034 method Methods 0.000 title claims abstract description 55
- 230000001788 irregular Effects 0.000 title claims abstract description 27
- 230000008569 process Effects 0.000 claims description 9
- 230000003247 decreasing effect Effects 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 3
- 230000002452 interceptive effect Effects 0.000 claims description 2
- 230000006870 function Effects 0.000 description 13
- 238000005516 engineering process Methods 0.000 description 12
- 238000011160 research Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 238000000605 extraction Methods 0.000 description 3
- 239000002243 precursor Substances 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000007667 floating Methods 0.000 description 2
- 238000012800 visualization Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000013523 data management Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开了一种不规则三维地质几何体的自动建模方法,它包括以下步骤:在当前观测的二维地震波形剖面上,新建一个初始块体窗口;利用浏览工具进行二维地震波形剖面的单向、连续、单间隔浏览,在紧接着的下一个剖面上,拷贝上一个剖面的块体窗口到此剖面上作为此剖面的初始块体窗口;采用自动追踪算法,自动调整拷贝来的块体窗口中的所有关键点,并保存调整结果作为本剖面的块体窗口和下一个剖面的初始块体窗口;以此类推,浏览每一个二维地震波形剖面,为每一个二维地震波形剖面建立块体窗口;依次保存浏览过程中生成的块体窗口序列,构成一系列封闭的块体窗口集合;通过这些块体窗口集合在体数据空间内进行构建三维地质几何体。
Description
技术领域
本发明涉及一种不规则三维地质几何体的自动建模方法。
背景技术
在石油地质勘探领域,计算机已经应用到地震、测井、岩芯图像分析、数据管理等各个方面。随着勘探技术的发展和勘探范围的不断扩大,油气勘探面临着地表条件多样、地下构造复杂的困难。面对复杂、隐蔽的勘探目标,地球物理科学需要尽可能地综合所有相关学科的技术,进行多次反馈、交叉研究,而许多油田专家队伍之间的工作方式是流水线式的。这种工作方式虽然使各学科的职能非常明确,但不利于学科之间信息交流与反馈,一旦有新的资料加入时,很难返回重新工作。因此,地球科学界的专家队伍一直呼吁重视软件通讯、公共数据通道和三维可视化,这样三维地质建模技术就应运而生了。随着计算机技术的飞速发展,三维地质建模技术越来越受到人们的重视,并成为研究的一个热点。
三维地质建模技术是指运用计算机图形学和图像处理技术,在三维环境下,将地质空间信息管理、分析和预测、地质解译、地学统计、实体内容分析及图形可视化等工具结合起来,应用于地质分析的技术。这方面的研究在国外开展得较早,到目前,已经形成了相当的规模,各类软件层出不穷,比如,GOCAD、Landmark、EarthVision、GeoQuest、GRISYS等,而国内在这方面的研究进行得较晚,很多工作才刚刚起步,典型的软件主要有Gristation和GIVE等。
传统的基于物体表面表示的面元图形学技术已经广泛的应用于三维地质建模软件中,虽然面元图形学技术能够表示物体的外部形状和相互之间的拓扑关系,但是,它无法表示物体的内部结构,而三维地质建模中很重要的一项就是三维空间中地质属性的计算与分析。随着体元(构成三维形体的单位,一般使用单位立方体、四面体、三棱柱、球体等描述)图形学技术(1993年由科学家Kaufman首次正式提出)的出现,它为解决三维空间属性的表示与计算提供了新的可能,目前,基于体元的三维地质建模方法已经成为了三维地质建模技术研究的重点。
体数据是体元的集合表示,可以认为体数据是体元模型的实际表现形式。在体元图形技术中,使用0或1代表当前坐标点上是否存在物体。在实际应用中,体元上的数据是该点所具有的属性值,如地质勘探中的地震数据。所以,一般意义上的体数据指带有该体元属性值的体元模型。由于体数据是实体数据,数据空间中的每一个点都有其对应的属性值。将属性体按某个切面获取数据,并按属性值进行显示,就可以得到对应的图像数据。
在三维地质几何体(如图1)建模中,规则几何体的建模方法已经比较成熟,而不规则几何体的建模是一个技术研究的难点。现有的不规则地质几何块体主要是通过在三维地震数据场(其数据都是地震数据,记录了勘探结果,可以用来分析和解释地质曲面位置与形态,可以通过水平和垂直切面的方法,得到对应的二维图像数据,其常用的显示方法是基于振幅的波形显示,如图2所示)中,把地质曲面(共有两种:层面,由一系列离散点表示的地质曲面,曲面的走向接近于水平方向;断层,由一系列离散点表示的地质曲面,曲面的走向接近于垂直方向,如图3、4所示)映射到体数据空间中,再通过层面和断层相交将体数据空间的区域划分,从而获得体数据几何模型的。这种不规则几何体的建模方法存在一定的问题,其中最为突出的是:当曲面相交时有时会出现“洞”的现象,而无法构成封闭的三维地质几何体,为解决这一问题,采用了一些诸如曲线、曲面逼近和近似的方法,但效果不佳,且易出错。
发明内容
鉴于上述原因,针对现有的不规则三维地质几何体建模方法所存在的缺陷,本发明的主要目的是提供一种借助三维地震数据场的信息,通过连续的二维封闭块体窗口构建三维几何体的不规则三维地质几何体自动建模方法。
为实现上述目的,本发明采用以下技术方案:一种不规则三维地质几何体的自动建模方法,它包括以下步骤:
1、利用鼠标点击的方式,在当前观测的二维地震波形剖面(序号为i(i>=1))上,新建一个初始块体窗口;
2、利用浏览工具进行二维地震波形剖面的单向、连续、单间隔浏览,在紧接着的下一个剖面(i+1)上,拷贝上一个剖面(i)的块体窗口到此剖面上作为此剖面的初始块体窗口;采用自动追踪算法,自动调整拷贝来的块体窗口中的所有关键点,并保存调整结果作为本剖面(i+1)的块体窗口和下一个剖面(i+2)的初始块体窗口;
3、以此类推,浏览每一个二维地震波形剖面,为每一个二维地震波形剖面建立块体窗口;依次保存浏览过程中生成的块体窗口序列,构成一系列封闭的块体窗口集合;
4、通过这些块体窗口集合在体数据空间内进行构建三维地质几何体;
5、最终得到的结果是不规则三维地质几何体。
本发明提出的一种不规则三维地质几何体的自动建模方法,采用了以“连续的二维封闭块体窗口构建三维几何体”的方法,从而解决了现存方法中的缺陷,提出了一种解决问题的新途径。该方法可“快速、简洁、自动”地生成复杂的不规则三维地质几何体,为其它地质工作的开展奠定了重要基础。
附图说明
图1为不规则三维地质几何体示意图
图2为含层面和断层标记的二维波形显示的地震剖面示意图
图3为地质层面示意图
图4为地质断层示意图
图5为本发明波形显示方式下的块体窗口示意图
图6为本发明不规则三维地质几何体建模方法流程图
图7为本发明关键点及其前驱和后继的8种情况以及对应的调整策略示意图
图8为本发明图7细的8种情况中的最后3种情况细化示意图
具体实施方式
为了准确描述本发明不规则三维地质几何体自动建模方法,首先给出三个基本概念:
●块体窗口:按一定顺序(依次顺时针或逆时针)的二维层面(层面在二维图像上的投影)关键点序列(折线的离散点序列,能反映折线的形态和走势)和二维断层(断层在二维图像上的投影)关键点序列(折线的离散点序列,能反映折线的形态和走势)的封闭体,如图5所示。特别地,若当前关键点是层面和断层的交点或者是层面的两个端点,本文默认把它们算作断层关键点,方便表示。
●块体窗口集:由一系列单独剖面合成的封闭块体窗口形成的集合体。
●地质几何体:由一系列离散点,在三维地质空间表示的封闭几何体。图1所示为一个不规则的地质几何体。
此外,所谓“连续”的含义是指在三维地震数据场中,沿着某一个方向(如水平或垂直方向),按照一定的距离间隔,依次做切面获取二维地震图像数据。
基于上面的基本概念,本发明公开的不规则三维地质几何体的自动建模方法包括以下步骤,如图6所示:
1、利用鼠标点击的方式,在当前观测的二维地震波形剖面(序号为i(i>=1))上,新建一个初始块体窗口;
2、利用浏览工具进行二维地震波形剖面的单向、连续、单间隔浏览,在紧接着的下一个剖面(i+1)上,拷贝上一个剖面(i)的块体窗口到此剖面上作为此剖面的初始块体窗口;采用自动追踪算法,自动调整拷贝来的块体窗口中的所有关键点,并保存调整结果作为本剖面(i+1)的块体窗口和下一个剖面(i+2)的初始块体窗口;
3、以此类推,浏览每一个二维地震波形剖面,为每一个二维地震波形剖面建立块体窗口;依次保存浏览过程中生成的块体窗口序列,构成一系列封闭的块体窗口集合;
4、通过这些块体窗口集合在体数据空间内进行构建三维地质几何体;
在形成了封闭的块体窗口集后,需要在其所处的地质三维数据工区内,构建这些块体窗口集,并最终得到地质几何体。封闭的块体窗口集只表示了块体的边界数据,并没有包含块体内部的数据。而地质几何体,既包括边界数据,也包括块体窗口内部的数据。所以,需要根据对应剖面的逐个块体窗口来提取需要的数据,并把每个剖面提取的数据按照提取的顺序组织起来,形成地质几何体的最终数据。
5、最终得到的结果是不规则三维地质几何体。
从上不难看出,自动追踪算法是本建模方法的关键,下面具体介绍一下该方法的技术实现方案:
1、首先,定义块体窗口Block-Window(以下简写为BW)的描述方式:
BW=(V,E)
其中:
V={δi,j|δi,j∈ξ,i=i0,j=1,2,...,n,n≥0},
E={<δi,j,δi,j+1>|δi,j,δi,j+1∈ξ,i=i0,j=1,...,n-1,n≥0,}U{<δi,n,δi,1>}。
在上面的描述中,定义了一个含有n个关键点、剖面序号为i0的块体窗口BW。在定义中,ξ为关键点类型元素的集合,V是块体窗口中层面/断层关键点的有穷非空集合,E是两个相邻关键点之间边的集合;δi,j表示BW中序号为j的关键点元素,可以由三元组来表示δi,j=(Xj,Yj,λj),其中点的坐标为(Xj,Yj);若λj=0表示此关键点为层面关键点,若λj=1表示此关键点为断层关键点。特别地,由于块体窗口是封闭体,所以<δi,n,δi,1>∈E;且每个关键点有且仅有一个前驱节点元素和一个后继节点元素,如序号为n的关键点的后继节点为序号为1的关键点,如序号为1的关键点的前驱节点为序号为n的关键点。并且关键点和其前驱和后继的横坐标均不同,即Xj-1≠Xj≠Xj+1。
在此定义的基础上,定义对块体窗口的十个基本操作:
①INITIATE(BW):初始化操作,设定一个空的块体窗口BW;
②COUNT(BW):统计块体窗口中的关键点数目函数,函数值为关键点数目;
③GET(BW,k):取关键点元素函数,若1≤k≤COUNT(BW),则函数值为BW中序号为k的关键点元素δi,k,否则为空元素NULL;
④PRIOR(BW,k):求前驱函数,若1<k≤COUNT(BW),则函数值为BW中序号为k-1的元素δi,k-1,否则为空元素NULL;由于BW的封闭性,第一个元素的前驱指定为最后一个元素;即δi,COUNT(BW)=PRIOR(BW,1);
⑤NEXT(BW,k):求后继函数,若1≤k<COUNT(BW),则函数值为BW中序号为k+1的元素δi,k+1,否则为空元素NULL;由于BW的封闭性,最后一个元素的后继指定为第一个元素,即δi,l=NEXT(BW,COUNT(BW));
⑥LOCATE(BW,δ):定位函数,若BW中存在和δ完全一致的元素,则函数值为该元素在BW中的序号,否则为零;
⑦INSERT(BW,k,δ):前插操作,在BW中序号为k的元素之前插入一个新的元素δ,此操作仅在1≤k≤COUNT(BW)+1时可行,BW的元素总数目加1;
⑧DELETE(BW,k):删除操作,删除BW中序号为k的元素,此操作仅在1≤k≤COUNT(BW)时可行,BW的元素总数目减1;
⑨EMPTY(BW):判空函数,若BW为空,则返回布尔值TRUE,否则返回布尔值FALSE;
⑩CLEAR(BW):BW置空操作,无返回值。
利用以上的形式化定义和这些基本操作,可以组合形成多种复杂的操作和算法,具有良好的可扩展性。
2、根据调整策略,自动调整拷贝来的块体窗口(此剖面的初始块体窗口)中的所有关键点
调整策略也是自动追踪算法实现的关键依据,因为它决定着相邻剖面块体窗口之间的对应关系。
由于块体窗口是由一系列关键点组成的,所以块体窗口之间的对应关系就转化为相邻剖面关键点的对应关系集合。逐个解决块体窗口之间每个关键点的对应关系,也就实现了块体窗口之间的对应关系。由于关键点可以分为两类:层面关键点和断层关键点,所以针对这两类关键点分别制定调整策略。由于需要逐个对关键点进行调整,所以要找出其与前驱和后继的关系;这种关系的情况总共有
种,如图7所示:
①地质层面的形态在相邻剖面间具有局部相似性,变化浮动很小,且层面关键点反映了当前坐标局部范围内振幅的最大值;基于以上两点原因,相邻剖面间层面关键点的对应关系主要依据前一个剖面的坐标信息和当前剖面的振幅信息来制定,层面关键点的调整与其前驱和后继点无关,所以图6中前3种情况可以统一用一种调整策略,即“针对层面关键点的调整策略”。
②图7的第4、5种情况是块体窗口中所不可能出现的情况,地质形态不可能出现层面和断层的跳变,只有渐变过程,所以在块体窗口中某一种类型的关键点一旦出现,就至少连续出现2次。
③地质断层的形态在相邻剖面之间不具有局部相似性,且变化浮动很大,所以其调整策略与层面关键点调整策略无法统一。断层关键点与其前驱和后继的关键点有关,相邻剖面间断层关键点的对应关系主要依据其前驱和后继关键点来制定。根据前驱和后继关键点的类型对应为图7的最后3种情况,其调整策略为“针对断层关键点的调整策略”。
●针对层面关键点的调整策略(Auto-Tracing-Bedding-Surface)
设第i个剖面的层面关键点为δi,j,其坐标为(Xj,Yj),则其在第i+1个剖面所对应的关键点为δi+1,j,其坐标(Xj′,Yj′)满足以下2个约束条件:
●针对断层关键点的调整策略(Auto-Tracing-Fault)
设当前处理的断层关键点为δi,j,其前驱和后继分别为:δi,j-1和δi,j+1,它们的关系有如图7的最后3种情况,结合δi,j-1和δi,j+1的横坐标比较,进一步细化为图8的3种情况,相应的调整策略如下(其中参数
取值方法同层面关键点的调整策略):
■情况1:即图7的第1种情况,前驱为层面关键点,后继为断层关键点,即(λj-1=0)∧(λj+1=1),则调整策略由δi,j-1决定:
①若Xi-1<Xi
判断A(Xj+1,Y′)是否大于“0”?其中,
1.若成立:则向X增大方向以X递增1移动,在移动过程中调整Y′值,使每个中途移动点在
内振幅取得MAX值,移动直到振幅值为0或X=Xj+1为止,设此时坐标为(X′,Y′),则将(Xj,Yj)移动到(X′,Y′)即可;
2.若不成立:若A(Xj,Y′)>0,
则(Xi,Yj)移动到(Xj,Y′)即可;
②若Xi-1>Xi
判断A(Xj-1,Y′)是否大于“0”?其中,
1.若成立:则向X减小方向以X递减1移动,在移动过程中调整Y′值,使每个中途移动点在
内振幅取得MAX值,移动直到振幅值为0或X=Xj+1为止,设此时坐标为(X′,Y′),则将(Xi,Yj)移动到(X′,Y′)即可;
■情况2:即图8的第2种情况,(λj-1=1)∧(λj+1=0),调整策略由δi,j+1决定:
也分Xj+1<Xj和Xj+1>Xj的2个分支考虑,内部处理与情况1相同,此处不再重复。
■情况3:即图8的第3种情况,(λj-1=1)∧(λj+1=1),调整策略由δi,j-1和δi,j+1共同决定:
令
,将(Xj,Yj)移动到(X′,Y′)即可。
结合块体窗口(Block-Window)的定义、操作和调整策略,自动追踪算法可以表示为:
Function Auto-Tracing
{
//初始化块体窗口i,分配内存空间,人工调整
INITIATE(BWi);Manual-Adjust(BWi);
integer iIDBW=i+1;
while(bBrowse==1)//只要未停止就一直执行
{
INITIATE(BWiIDBW);
//若分配空间失败则跳出循环
if(TRUE==EMPTY(BWiIDBW))then break;
//拷贝前一个剖面的块体窗口到当前剖面上
Copy(BWiIDBW,BWiIDBW-I);
integer iNum=COUNT(BWiIDBW),iLoop;ξδB,δC,δN;
for(iLoop=1;iLoop<=iNum;iLoop++)
{
//由于断层关键点的调整依赖于其前驱和后继的层面关键点,所以要首先调整所
//有层面关键点
if(δCλC==0)Auto-Tracing-Bedding-Surface(δC);//层面关键点调整策略
}
for(iLoop=1;iLoop<=iNum;iLoop++)//再调整所有断层关键点
{
if(δCλC==1)
{
δB=PRIOR(BWiIDBW,iLoop);
δC=GET(BWiIDBW,iLoop);
δN=NEXT(BWiIDBW,iLoop);
Auto-Tracing-Fault(δB,δC,δN);//断层关键点调整策略
}
SET(BW,iLoop,δC);//将调整结果保存到BW中
}
iIDBW=iIDBW+1;//进入下一轮循环
}
Save(BWi,BWi+1,…,BWiIDBW);//保存所有块体窗口到文件中
CLEAR(BWi);…;CLEAR(BWiIDBW);//释放所有块体窗口所占内存空间}
算法中一个while循环中含有两个并列的for循环,由于单个地质几何体所对应的块体窗口的数量和单个块体窗口关键点数目差不多,且层面移动策略和断移动策略中涉及的关键点移动范围近似于常数C,所以,算法的时间和空间复杂性均为:O(n*n)=O(n2)。
算法目的是确保建模过程的自动性;在准确性得以保障的前提下,自动性带来了实际应用中的快速性。此外,算法基于制定的调整策略来实现,而策略本身具有良好的可管理、可扩充性,因而算法具有很大的拓展空间和实际价值。
在形成了封闭的块体窗口集后,需要在其所处的地质三维数据工区内,构建这些块体窗口集,并最终得到地质几何体。封闭的块体窗口集只表示了块体的边界数据,并没有包含块体内部的数据。而地质几何体,既包括边界数据,也包括块体窗口内部的数据。所以,需要根据对应剖面的逐个块体窗口来提取需要的数据,并把每个剖面提取的数据按照提取的顺序组织起来,形成地质几何体的最终数据。
本发明提出的一种不规则三维地质几何体的自动建模方法,采用了以“连续的二维封闭块体窗口构建三维几何体”的方法,从而解决了现存方法中的缺陷,提出了一种解决问题的新途径。该方法可“快速、简洁、自动”地生成复杂的不规则三维地质几何体,为其它地质工作的开展奠定了重要基础。
本发明提出的一种不规则三维地质几何体的自动建模方法与其它建模方法相比,具有以下三个显著特点:1、方法本身带来的自动性,缩短了层面、断层和不规则地质几何体的提取时间;2、减少人机交互的次数和频率,提高了执行效率;3、避免了现有的一些方法中无法实现的曲面相交后封闭的问题。
Claims (3)
1、一种不规则三维地质几何体的自动建模方法,其特征在于:它包括以下步骤:
第一步:利用鼠标点击的方式,在当前观测的二维地震波形剖面(序号为i(i>=1))上,新建一个初始块体窗口;
第二步:利用浏览工具进行二维地震波形剖面的单向、连续、单间隔浏览,在紧接着的下一个剖面(i+1)上,拷贝上一个剖面(i)的块体窗口到此剖面上作为此剖面的初始块体窗口;采用自动追踪算法,自动调整拷贝来的块体窗口中的所有关键点,并保存调整结果作为本剖面(i+1)的块体窗口和下一个剖面(i+2)的初始块体窗口;
第三步:以此类推,浏览每一个二维地震波形剖面,为每一个二维地震波形剖面建立块体窗口;依次保存浏览过程中生成的块体窗口序列,构成一系列封闭的块体窗口集合;
第四步:通过这些块体窗口集合在体数据空间内进行构建三维地质几何体;
第五步:最终得到的结果是不规则三维地质几何体。
2、根据权利要求1所述的一种不规则三维地质几何体的自动建模方法,其特征在于:所述块体窗口简写为BW按照以下方式描述:
BW=(V,E)
其中:
V={δi,j|δi,j∈ξ,i=i0,j=1,2,...,n,n≥0},
E={<δi,j,δi,j+1>|δi,j,δi,j+1∈ξ,i=i0,j=1,...,n-1,n≥0,}U{<δi,n,δi,l>}
该块体窗口含有n个关键点、剖面序号为i0,其中,ξ为关键点类型元素的集合,V是块体窗口中层面/断层关键点的有穷非空集合,E是两个相邻关键点之间边的集合;δi,j表示BW中序号为j的关键点元素,可以由三元组来表示δi,j=(Xj,Yj,λj),其中点的坐标为(Xj,Yj);若λj=0表示此关键点为层面关键点,若λj=1表示此关键点为断层关键点;由于块体窗口是封闭体,所以<δi,n,δi,1>∈E;且每个关键点有且仅有一个前驱节点元素和一个后继节点元素,如序号为n的关键点的后继节点为序号为1的关键点,如序号为1的关键点的前驱节点为序号为n的关键点,并且关键点和其前驱和后继的横坐标均不同,即Xj-1≠Xj≠Xj+1。
3、根据权利要求1所述的一种不规则三维地质几何体的自动建模方法,其特征在于:所述自动跟踪算法是指根据调整策略,自动调整拷贝来的块体窗口中的所有关键点;
所述调整策略分为层面关键点调整策略和断层关键点调整策略:
所述层面关键点的调整策略:
设第i个剖面的层面关键点为δi,j,其坐标为(Xj,Yj),则其在第i+1个剖面所对应的关键点为δi+1,j,其坐标(Xj′,Yj′)满足以下2个约束条件:
所述断层关键点的调整策略:
设当前处理的断层关键点为δi,j,其前驱和后继分别为:δi,j-1和δi,j+1:
■情况1:前驱为层面关键点,后继为断层关键点,即(λj-1=0)^(λj+1=1),则调整策略由δi,j-1决定:
①若Xj-1<Xj
若成立:则向X增大方向以X递增1移动,在移动过程中调整Y′值,使每个中途移动点在
内振幅取得MAX值,移动直到振幅值为O或X=Xj+1为止,设此时坐标为(X′,Y′),则将(Xj,Yj)移动到(X′,Y′)即可;
若不成立:若A(Xj,Y′)>0,
则(Xj,Yj)移动到(Xj,Y′)即可;若A(Xj,Y′)≤0,
则向X减小方向以X递减1移动,在移动过程中在
内调整Y′值,移动直到振幅值大于0或X=Xj-1为止,设此时坐标为(X′,Y′),则将(Xj,Yj)移动到(X′,Y′)即可。
②若Xj-1>Xj
判断A(Xj-1,Y′)是否大于“0”?其中,
若成立:则向X减小方向以X递减1移动,在移动过程中调整Y′值,使每个中途移动点在
内振幅取得MAX值,移动直到振幅值为0或X=Xj+1为止,设此时坐标为(X′,Y′),则将(Xl,Yj)移动到(X′,Y′)即可;
若不成立:若A(Xj,Y′)>0,
则(Xj,Yj)移动到(Xj,Y′)即可;若A(Xj,Y′)≤0,
则向X增大方向以X递增1移动,在移动过程中在
内调整Y′值,移动直到振幅值大于0或X=Xj-1为止,设此时坐标为(X′,Y′),则将(Xj,Yj)移动到(X′,Y′)即可。
■情况2:(λj-1=1)^(λj+1=0),调整策略由δi,j+1决定:
也分Xj+1<Xj和Xj+1>Xj的2个分支考虑,内部处理与情况1相同,此处不再重复。
■情况3:(λj-1=1)^(λj+1=1),调整策略由δi,j-1和δi,j+1共同决定:
令
将(Xj,Yj)移动到(X′,Y′)即可。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2005101145649A CN1956011B (zh) | 2005-10-26 | 2005-10-26 | 不规则三维地质几何体的自动建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2005101145649A CN1956011B (zh) | 2005-10-26 | 2005-10-26 | 不规则三维地质几何体的自动建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1956011A true CN1956011A (zh) | 2007-05-02 |
CN1956011B CN1956011B (zh) | 2010-08-11 |
Family
ID=38063316
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2005101145649A Expired - Fee Related CN1956011B (zh) | 2005-10-26 | 2005-10-26 | 不规则三维地质几何体的自动建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1956011B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101105864B (zh) * | 2007-08-07 | 2010-04-07 | 中国地质大学(武汉) | 一种基于序列剖面拓扑推理的三维地质体动态重构方法 |
CN101582173B (zh) * | 2009-06-24 | 2012-07-11 | 中国石油天然气集团公司 | 复杂地质构造块状模型构建方法 |
CN102592313A (zh) * | 2011-12-29 | 2012-07-18 | 中国水电顾问集团华东勘测设计研究院 | 一种地质三维模型自动建模与动态更新的方法 |
CN104252547B (zh) * | 2013-06-26 | 2017-06-06 | 中国石油化工股份有限公司 | 一种大规模油藏地质模型数据快速载入方法 |
CN106845040A (zh) * | 2017-03-21 | 2017-06-13 | 中国电建集团成都勘测设计研究院有限公司 | 不规则地质体三维精细化表达方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2833384B1 (fr) * | 2001-12-10 | 2004-04-02 | Tsurf | Procede, dispositif et produit programme de modelisation tridimensionnelle d'un volume geologique |
CN1156797C (zh) * | 2002-10-23 | 2004-07-07 | 北京航空航天大学 | 在虚拟环境多因素影响下生成逼真光照效果的方法 |
CN100395779C (zh) * | 2003-10-14 | 2008-06-18 | 武汉大学 | 基于立体影像镶嵌的可量测无缝空间立体模型的生成方法 |
CN100394449C (zh) * | 2004-06-30 | 2008-06-11 | 南京大学 | 基于gis与虚拟现实的实时三维地质建模的方法 |
-
2005
- 2005-10-26 CN CN2005101145649A patent/CN1956011B/zh not_active Expired - Fee Related
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101105864B (zh) * | 2007-08-07 | 2010-04-07 | 中国地质大学(武汉) | 一种基于序列剖面拓扑推理的三维地质体动态重构方法 |
CN101582173B (zh) * | 2009-06-24 | 2012-07-11 | 中国石油天然气集团公司 | 复杂地质构造块状模型构建方法 |
CN102592313A (zh) * | 2011-12-29 | 2012-07-18 | 中国水电顾问集团华东勘测设计研究院 | 一种地质三维模型自动建模与动态更新的方法 |
CN102592313B (zh) * | 2011-12-29 | 2015-04-01 | 中国水电顾问集团华东勘测设计研究院 | 一种地质三维模型自动建模与动态更新的方法 |
CN104252547B (zh) * | 2013-06-26 | 2017-06-06 | 中国石油化工股份有限公司 | 一种大规模油藏地质模型数据快速载入方法 |
CN106845040A (zh) * | 2017-03-21 | 2017-06-13 | 中国电建集团成都勘测设计研究院有限公司 | 不规则地质体三维精细化表达方法 |
CN106845040B (zh) * | 2017-03-21 | 2019-12-03 | 中国电建集团成都勘测设计研究院有限公司 | 不规则地质体三维精细化表达方法 |
Also Published As
Publication number | Publication date |
---|---|
CN1956011B (zh) | 2010-08-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR102040988B1 (ko) | 복수의 3d 곡선들로부터의 표면 생성 | |
CN104620136B (zh) | 组合多点统计和基于对象的方法的计算机系统和混合方法 | |
CN1734503A (zh) | 使用光谱分析的伸展驱动的网格参数化 | |
US11163915B2 (en) | Three-dimensional modeled object | |
CN1967596A (zh) | 三维扫描系统中三维散乱点集的三角剖分构造方法 | |
CN101051395A (zh) | 一种基于地球物理场数据的地质体三维可视化方法 | |
CN101051394A (zh) | 一种基于地球物理场数据的地质体三维可视化系统 | |
CN1604142A (zh) | 基于图像的着色和编辑方法及装置 | |
Ying et al. | An intrinsic algorithm for parallel poisson disk sampling on arbitrary surfaces | |
CN1956011A (zh) | 不规则三维地质几何体的自动建模方法 | |
EP3119986A2 (en) | Modeling intersecting faults and complex wellbores in reservoir simulation | |
CA3065179A1 (en) | Constructing structural models of the subsurface | |
AU2016346919B2 (en) | Interactive salt model modification | |
Marcias et al. | Data-driven interactive quadrangulation | |
CN108898670B (zh) | 一种基于实体和剖面的三维地质建模方法 | |
CN1873647A (zh) | Cad方法、cad系统以及存储其cad程序的程序存储介质 | |
Konyha et al. | Interactive visual analysis of families of curves using data aggregation and derivation | |
CN1430184A (zh) | 利用半边数据结构实现三维网格模型的简化方法 | |
CN1187697C (zh) | 硅基微器件计算机辅助工艺规划方法 | |
CN1311390C (zh) | 三维形状数据向单元内部数据的变换方法 | |
CN1234091C (zh) | 基于水平集和分水岭方法的医学图像分割方法 | |
CN113420359B (zh) | 一种基于revit的电力工程总图自动传参布井的方法 | |
CN110349262A (zh) | 一种基于bim的隧道超前地质预报三维建模方法及系统 | |
Becher et al. | Feature-based volumetric terrain generation and decoration | |
US11686870B2 (en) | Interpretive-guided velocity modeling seismic imaging method and system, medium and device |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100811 Termination date: 20121026 |