发明名称 复杂空间的曲面重构方法
摘要 本发明公开了一种复杂空间的曲面重构方法,通过采用复杂空间曲面插值方法,解决了三维地质建模中复杂地质构造下的层位曲面重构问题,为等值线绘制、地质成块等提供了新的思路。本发明具有如下优点:解决了复杂地质构造下的层位插值问题,不但支持正、逆断层以及垂直断层,并且适合于多种地质建模应用;通过断层多边形边界约束下的层面拟合,实现了层面和断层无缝拟合,解决了传统的网格化插值方案层位与断层无法严格相交的难题;通过采用混合网格化的层位插值方法,既保证了数据的快速处理,又保证了拟合的精度。
申请公布号 CN102222365B 申请公布日期 2012.08.29
申请号 CN201110215270.0 申请日期 2011.07.29
申请人 电子科技大学 发明人 鲁才;陈凯云;胡光岷
分类号 G06T17/05(2011.01)I 主分类号 G06T17/05(2011.01)I
代理机构 电子科技大学专利中心 51203 代理人 周永宏
主权项 1.一种复杂地质构造的曲面重构方法,其特征在于:包括如下步骤:第一步、原始断层数据预处理:包括断层插值、拟合断层面和根据层位原始三维数据生成关联层位的断层多边形;其中断层插值和拟合断层面是先确定一个断层包络范围,利用克里金算法直接插值,然后根据插值数据拟合成断面;所述断层多边形是断层与层位在空间上相交而成的一个多边形结构,其边界由断层线组成,断层线分为上盘线和下盘线,上下盘线组成一个三维空间中的闭合多边形;第二步、原始层位数据预处理:去除违法的点数据;第三步、层位插值:(1)初始化关联的断层多边形:获取与将要插值的层位关联的所有断层多边形,初始化层位插值的约束断层多边形边界集合;(2)初始化网格属性:确定网格参数,将关联断层多边形边界按照网格离散为三维控制点,并将关联断层多边形投影到网格,确定各个网格的属性;(3)初始化种子点数据:获取将要插值的层位原始数据和关联断层多边形边界的三维控制点,转换为离散的三维种子点数据结构,获取每个种子点的断层关联属性,然后将种子点按照(x,y)坐标投影到各个网格中,定义种子点结构,获得初始化种子点的属性;(4)初始化插值点数据:在插值之前对网格的每一道进行插值点初始化,确定网格的每一道需要插值的层位点数目,以及每个插值点的断层关联属性,定义插值点结构,获得初始化插值点的属性;(5)层位插值:根据第(4)步得到的插值点数据,通过搜索插值点周围的合法种子点来对该插值点进行插值计算;在种子点稀疏的层位中,保留无法计算Z值的插值点,在所有能依靠种子点进行插值计算的插值点处理完毕后,将插值完成的点当作新的种子点来进行层位补偿插值,直到所有插值点均处理完毕;所述合法种子点的判断规则为:设待插值点为I<sub>i</sub>,关联的断层集合uF<sub>i</sub>和lF<sub>i</sub>;搜索到的种子点为S<sub>j</sub>,关联的断层集合uF<sub>j</sub>和lF<sub>j</sub>;定义与插值点和种子点的水平连线I<sub>i</sub>S<sub>j</sub>相交的断层多边形边界集合Boundary,其中uP为断层上盘边界集合,lP为断层下盘边界集合,up<sub>k</sub>为断层F<sub>k</sub>关联插值层位的多边形上盘边界,lp<sub>k</sub>为断层F<sub>k</sub>关联插值层位的多边形下盘边界,m为插值层位关联的断层多边形总数:Boundary={uP,lP}<maths num="0001"><![CDATA[<math><mfenced open='' close=''><mtable><mtr><mtd><mi>uP</mi><mo>=</mo><mo>{</mo><msub><mi>up</mi><msub><mi>i</mi><mn>1</mn></msub></msub><mo>,</mo><msub><mi>up</mi><msub><mi>i</mi><mn>2</mn></msub></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>up</mi><msub><mi>i</mi><mi>M</mi></msub></msub><mo>}</mo></mtd><mtd><mn>0</mn><mo>&le;</mo><msub><mi>i</mi><mi>M</mi></msub><mo>&le;</mo><mi>m</mi></mtd></mtr></mtable></mfenced></math>]]></maths><maths num="0002"><![CDATA[<math><mfenced open='' close=''><mtable><mtr><mtd><mi>lP</mi><mo>=</mo><mo>{</mo><msub><mi>lp</mi><msub><mi>k</mi><mn>1</mn></msub></msub><mo>,</mo><msub><mi>lp</mi><msub><mi>k</mi><mn>2</mn></msub></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>lp</mi><msub><mi>k</mi><mi>N</mi></msub></msub><mo>}</mo></mtd><mtd><mn>0</mn><mo>&le;</mo><msub><mi>k</mi><mi>N</mi></msub><mo>&le;</mo><mi>m</mi></mtd></mtr></mtable></mfenced></math>]]></maths>A.若<img file="FDA0000158725510000023.GIF" wi="146" he="38" />且<img file="FDA0000158725510000024.GIF" wi="159" he="39" />则种子点S<sub>j</sub>合法;B.若存在断层F<sub>k</sub>∈uF<sub>i</sub>且F<sub>k</sub>∈lF<sub>j</sub>,或者F<sub>k</sub>∈lF<sub>i</sub>且F<sub>k</sub>∈uF<sub>j</sub>,则种子点S<sub>j</sub>不合法;C.若Boundary中存在垂直断层边界,则种子点S<sub>j</sub>不合法;D.若Boundary中同时存在某个断层的上盘边界和下盘边界,则种子点S<sub>j</sub>不合法;E.若uF<sub>j</sub>中最后一个断层或者lF<sub>j</sub>中第一个断层为正断层,且Boundary中包含该正断层的边界,则种子点S<sub>j</sub>不合法;F.若Boundary集合不满足D和E中的条件,同时有uF<sub>i</sub>=uF<sub>j</sub>且lF<sub>i</sub>=lF<sub>j</sub>,则种子点S<sub>j</sub>合法;G.若存在断层F<sub>k</sub>∈uF<sub>i</sub>且lp<sub>k</sub>∈F<sub>k</sub>且lp<sub>k</sub>∈lP,或者F<sub>k</sub>∈lF<sub>i</sub>且up<sub>k</sub>∈F<sub>k</sub>且up<sub>k</sub>∈uP,则种子点S<sub>j</sub>不合法;H.若存在断层F<sub>k</sub>∈uF<sub>j</sub>且lp<sub>k</sub>∈F<sub>k</sub>且lp<sub>k</sub>∈lP,或者F<sub>k</sub>∈lF<sub>j</sub>且up<sub>k</sub>∈F<sub>k</sub>且up<sub>k</sub>∈uP,则种子点S<sub>j</sub>不合法;I.其它情况,则判定种子点S<sub>j</sub>合法;第四步、层位拟合成面:根据事先确定的边界约束条件将层位与断层相邻部分的边界插值点与断层多边形的上下盘线段连接起来,得到一个层面和断层完全无缝连接的面结构。
地址 611731 四川省成都市高新区(西区)西源大道2006号