用户名: 密码: 验证码:
复杂矿体结构三维建模与储量计算方法研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
在各种矿产资源中,铀矿资源既是一种战略资源,又是重要的能源矿产,在国防和国民经济建设中具有非常重要的地位。和世界上其他国家比较,我国是铀矿资源不甚丰富的一个国家,勘探手段和程度都远低于发达国家水平,急需加强资源勘探与开发综合利用技术研究。随着计算机信息技术特别是地理信息技术与三维图形技术的发展,将地理信息技术、三维地质建模技术及三维可视化技术应用到固体矿产勘探、开发与管理中,是目前矿产资源勘探开发中的重要研究方向。
     论文以某地区铀矿床为研究对象,通过分析该区域的矿床结构特征和勘探数据特点,建立了一套针对复杂铀矿床的三维地质建模及储量计算方法技术体系,并开发了铀矿三维建模及储量计算可视化软件系统。通过对该矿区的矿床建模和储量计算实验,验证了方法技术的正确性,对于指导该区域的进一步勘探与开发具有重要的实用意义。
     论文针对铀矿结构特征和勘探数据特点,提出了采用三角网构模和广义三棱柱体元构模方法,分别用于构建矿体的表面模型和建立矿体内部属性模型;在三维建模中,采取了基于剖面重构的三维建模方法。考虑铀矿特殊的成矿条件,在勘探线上矿体形态呈现复杂,边缘尖锐,以及在相邻勘探线上矿体形态差异较大等特点,论文提出了一种基于矿体轮廓多边形变形技术的三维剖面重建方法。首先从勘探线中提取出矿体的轮廓边缘,经过将复杂多边形简化处理成简单多边形,然后采用同构三角网剖分及多边形插值技术生成中间多边形,最后采用等分线法将中间多边形对应连接生成三维矿体的表面。在矿体内部属性建模中,为了更好地适应矿体的边界,论文采用了基于广义三棱柱体元建模方法,并实现了三维网格剖分、网格节点插值,以及矿体边界处理等关键技术问题,成功构建了矿体内部属性模型。在对矿区储量估算中,论文提出了结合三维地质模型的Kriging储量计算方法。通过计算研究区样本数据的实验变异函数,利用Kriging插值算法生成了矿体三维空间的品位分布;然后利用矿体表面模型对矿体品位边界进行约束处理,最后通过统计矿体内体元数目和资源量得到矿体的储量。论文还研究并实现了在建模过程中涉及到的各种三维可视化与交互技术,如MC算法、LOD算法,三维图形的选择、旋转、缩放,切片制作、图形切割及虚拟漫游等,提高了建模过程的自动化,增强了三维操作的互动性与实境感。
In various minerals, uranium resources is not only a strategic resource, but alsois an important energy resources, it has a very important position in the nationaldefense and national economy construction. Compared with other countries in theworld, China is not rich of uranium resources, the degree of exploration is far belowthe level of the developed countries, and needs to strengthen the resources explorationand the technology research of comprehensive utilization. With the development ofcomputer information technology,especially the development of geographicinformation technology and3D graphics technology, geographic informationtechnology,3D geological modeling technology and3D visualization technology areapplied to the development and management of solid mineral exploration, which isone of the important research direction of reservation and calculation of mineralresources, it can improve the utilization efficiency of resources, enhance theexploration level and reduce the cost of exploration and management.
     This paper takes uranium ore field of Ruergai area in Sichuan Province as theresearch object. Through the analysis of survey data and deposit structure features inthis area, established a set of3D geological modeling of complex uranium deposit andreserves estimation methods technology system, and also developed the3D modelingand reserves of uranium calculation visualization software system. Through thecomparision of the ore deposit modeling and reserves estimation experiments, verifiedthe correctness of the method and technology, which has an important practicalsignificance for further exploration and development in this area.
     This paper analysis the structure characteristics of mining uranium andexploration data features, proposes the triangular mesh model and the generalizedthree prism element modeling method, which were used to establish surface model orebody and internal attribute distribution of ore body; In the3D modeling, takes the3Dmodeling method based on the profile reconstruction. Considering the specialmetallogenic conditions of the uranium, complex, in the exploration line the edgesharp,as well as the great difference in characteristics of the adjacent line shape of oreexploration, this paper presents a method for3D reconstruction based on contourdeformation technology. First extracted from contour edge orebody exploration line inpolygon, after simplified into a simple polygon, and then the compatible triangulationand polygonal interpolation technique to generate the intermediate polygon, thecorresponding connection to generate3D orebody surface by the bisector method. Inthe internal property modeling, in order to better adapt to the boundary of ore body, the volume element modeling method based on generalized three prism, and therealization of the technology of3D mesh, grid interpolation and boundary treatmentmethods, to construct the internal property model of ore. In the estimation of miningarea reserves,this paper proposed the Kriging reserve estimation method of thecombination of3D geological model. Through the experimental variation functioncalculation of the sample data, using Kriging interpolation algorithm to generate3Dorebody distribution; Then the constraint on the boundary of the orebody tasteorebody surface model, finally orebody reserve produced by the statistical elementnumber and the amount of resources. This paper also study and realize the3Dvisualization and interaction techniques in the process of modeling, such as MCalgorithm, LOD algorithm,3D graphics options, rotation, scaling, paraffin sections,graph cut and virtual roaming.
引文
[1] Simon W Houlding.3D geoscientific modeling-computer technique forgeological characterization[M].Springer-Verlag,1994.
    [2] Rongxing Li.Data structures and application issues in3D geographicinformation systems[J].GEOMATICA.1994,48(3):209~224.
    [3]王明华.工程岩体三维地质建模与可视化研究[D].中国科学院武汉岩土力学研究所(武汉)博士论文.2003.
    [4]王明华,白云.三维地质建模研究现状与发展趋势[J].土工基础.2006,8.vol.20.No.4.
    [5]吴立新,史文中,Christopher Gold.3D GIS与3D GMS中的空间构模技术[J].地理与地理信息科学,2003,19(1):5-11.
    [6] Christopher C M,Maydell U M. Triangulation and spatial ordering incomputer cartography[C].Proceedings of3rd Canadian CartographicAssociation Annual Meeting.Canada,1978.69-8l.
    [7] Nicolas S,Renato P. Delaunay Triangulation of Arbitrarily ShapedPlanar Domains [J].Computer Aided Geometric Design,1991,8:421-437.
    [8] Victor J D. Delaunay Triangulation in TIN Creation: A Overview andA Linear-Time Algorithm[J].GIS,1993,7(6):501-524.
    [9]赫尔丁S W.三维矿床的计算机构模方法[J].国外金属矿山,1989(2):95-96.
    [10]许斌,张森,历万庆.从序列切片重构三维对象的新方法[J].计算机学报,1994,17(1):64-71
    [11]惠勒A J,斯托克斯P C.块段模型和线框模型在地下采矿中的应用[J].国外金属矿山,1989(2):98-101.
    [12]韩国建,郭达志,金学林.矿体信息的八叉树存储和检索技术[J].测绘学报.1992.21(1):13-17.
    [13]赵树贤.煤矿床可视化构模技术[D].中国矿业大学(北京校区),1999.
    [14]张珊珊,刘志辉.基于多层DEM表面模型的地层结构的三维可视化[J].测绘信息与工程,2003,28(3):14-15.
    [15]吴立新,沙从术.真三维地学模拟系统与水利工程应用[J].南水北调与水利科技,2003,1(2):20-25.
    [16]刘陵,方军,陈利生,张洋.三维GIS的研究现状及其发展趋势[J].矿山测量,2011,(2):75.
    [17]贺怀建,白世伟.三维地层模型中地层划分的探讨[J].岩土力学,2002,23(5):637-639.
    [18]张煜,谭德宝,孙家柄.基于竖直钻孔数据的层状地层建模[J].长江科学院院报,2005,22(6):71-74.
    [19]陈学习,吴立新,车德福等.基于钻孔数据的含断层地质体三维建模方法[J].煤田地质与勘探,2005,33(5):5-8.
    [20]朱合华;张芳;叶勇庚;基于钻孔数据重构地层周围表面模型算法[J].计算机工程与应用,2006,(25),213-216.
    [21]张渭军,王文科.基于钻孔数据的地层三维建模与可视化研究[J].大地构造与成矿学,2006,30(1):108-113.
    [22]罗智勇,杨武年.基于钻孔数据的三维地质建模与可视化研究[J].测绘科学,2008,33(2):130-132.
    [23]王媛妮,葛非,李向.基于稀疏钻孔点的地层建模及可视化[J].测绘信息与工程,2008,33(1):18-20.
    [24]孟永东,徐卫亚,田斌,李明卫.基于带约束三角剖分的三维地质建模方法及应用[J].系统仿真学报,2009,21(19):5985-5989.
    [25]向中林,王妍,王润怀等.基于钻孔数据的矿山三维地质建模及可视化过程研究[J].地质与勘探,2009,45(1):75-80.
    [26]罗周全,罗贞焱,刘晓明,李畅.基于钻孔编录信息可视化解译的矿床三维建模[J].中南大学学报(自然科学版),2010,41(6):2367-2372.
    [27]高琴,洪振刚,李朝晖,龙桂华.基于VTK的三维地质建模方法研究[J].工程勘察,2011,(4):59-63.
    [28]Christian Le Carlier de Veslud, Michel Cuney, Guillaume Lorilleux,Jean-Jacques Royer, Michel Jébrak.3D modeling of uranium-bearingsolution-collapse breccias in Proterozoic sandstones (AthabascaBasin, Canada)—Metallogenic interpretations[J].2009,35(1).
    [29]李俊.固体矿产资源/储量估算方法研究现状及主要方法评介[J].中国非金属矿工业导刊,2005,46(2):53-55.
    [30]唐义,蓝运蓉.SD储量计算法[D],北京地质出版社,1990.
    [31]马洪滨,熊俊楠.基于地质统计学的储量估算系统[J].煤炭学报,2007,32(3):267-271.
    [32]王恩德,孙立双,蔡洪春等.矿体三维数据模型及品位插值方法研究[J],地质与资源,2007,16(3):222-224.
    [33]刘海英,刘修国,李超岭.基于地质统计学法的三维储量估算系统研究与应用[J].吉林大学学报(地球科学版),2009,39(3):541-546.
    [34]陈晓鹏,魏丽英.基于GTP—TEN模型矿体储量的计算研究[J].矿业工程,2010,8(2):60-62.
    [35]马翀,潘懋,金毅,姚凌青.三维地学建模技术与Kriging方法在固体矿产储量计算中的应用[J].北京大学学报(自然科学版),2011,47(3):483-490.
    [36]张明明,李晓晖,周涛发等.基于三维矿化域模型的泥河铁矿床动态储量估算[J].地质评论,2013,59(1):122-128.
    [37]侯恩科,吴立新.面向地质建模的三维体元拓扑数据模型研究[J].武汉大学学报信息科学版,2002,27(5):467-472.
    [38]杜培军,郭达志,田艳凤.顾及矿山特性的三维GIS数据结构与可视化[J].中国矿业大学学报(自然科学版),2001,30(3):238-243.
    [39]简兴祥.大地电磁二三维建模与数值模拟分析[D].成都理工大学博士论文,2003.
    [40] Keppel E. Approximating complex surfaces by triangulationofcontour lines[J]. IBM Journal of Research andDevelopment,1975,19(1):2-11.
    [41] Fuchs H, Kedem Z M, Uselton S P. Optimal surfaces reconstructionfrom planar contours[J]. Communication of theACM,1977,20(10):693-702.
    [42] Christiansen H N, Sederberg T W. Conversion of complex contourline definitions into polygonal element mosaics [J].ComputerGraphics,1978,12(3):187-192.
    [43]刘刚,胡远来,邓林.层迭三维体表面重构的算法探讨[J].成都理工大学学报:自然科学版,2003,30(5):537–540.
    [44]唐泽圣.三维数据场可视化[M].北京:清华大学出版社.2000.
    [45]张珊珊.多层DEM表面模型的建立与管理[D].武汉大学硕士论文,2002.
    [46]汪媛媛,杨忠芳,余涛.土壤碳储量计算中不同插值方法对比研究——以吉林省大安市为例[J].中国岩溶,2011,30(4):479-486.
    [47] Barequet G,Sharir M. Piecewise linear Interpolation betweenPolygonal Slices. Proc10th Ann ACM Symp on ComputationalGeometry(SoCG),1994,93~102.
    [48] Barequet G,Sharir M. Piecewise linear Interpolation betweenPolygonal Slices[J]. Computer Vision ImageUnderstanding,1996,63:251~272.
    [49] Barequet G,Shapiro D,Tal A.Multilevel Sensitive Reconstruction ofPolyhedral Surfaces from Parallel Slices[J].The VisualComputer,2000,16:116~133.
    [50] Meyers D,Skinner S.Surfaces from Contours[J].ACM Transactions onGraphics,1992,11(3):228~258.
    [51] Levin D.Multidimensional Reconstruction by set valuedApproximation[J].IMA J Numer Anal,1986,6:173~184.
    [52]陈学工,黄伟,季兴,李楠.一种由轮廓线重建物体表面的方法[J].计算机工程与应用,2011,47(11):157-159.
    [53]于杰,王怀明.基于DXF文件的CAD/CAM接口程序开发与应用[J].制造业自动化,2008,30(7):79-81.
    [54]李芳珍,许伦辉.DXF文件格式及其外部接口的研究[J].软件开发与应用,2008,27(7):83-85.
    [55]张公,周荣福,李祥.AutoCAD数据向MapGIS数据转换方法的研究[J].城市勘测,2011,(3):126-129.
    [56]邹福辉,李忠科.图像边缘检测算法的对比分析[J].计算机应用,2008,28(1):215-219.
    [57]孙惠,周红霞,李朝辉.图像处理中边缘脑开发与应用[J].电脑开发与应用,2002,(10):1-3.
    [58]马艳,张治辉.几种边缘检测算子的比较[J].工矿自动化,2004,1(2):54-56.
    [59]王淑琴,张金海,贺锡雷,邓飞.最简边缘检测算子及其在地质建模中的应用[J].计算机应用研究,2005,176-180.
    [60]王小华,钱月晶.一种改进的Canny边缘检测算法[J].机电工程,2008,25(12):60-63.
    [61]薛丽霞,李涛,王佐成.一种自适应的Canny边缘检测算法[J].计算机应用研究,2010,27(9):3588-3591.
    [62]王结臣,沈定涛,陈焱明,李利番.一种有效的复杂多边形裁剪算法.武汉大学学报信息科学版,2010,35(3):369-372.
    [63]叶琳,邱龙辉.多边形合并的算法研究[J].计算机应用与软件,2002,(8):57-59.
    [64]刘亚姝,严寒冰,范友贵.多边形简化算法及比较[J].计算机工程,2009,35(23):227-229.
    [65]赖志豪,康宝生.二维物体变形技术现状与发展.计算机技术与发展,2007,17(10):120-126.
    [66] Beier T,Neely S. Feature based image metamorphosis[J].ComputerGraphics,1992,26(2):35-42.
    [67] Sederberg T W,Greenwood E.A physically based approach to2-D shapeblending[J].Computer Graphics,1992,26(2):25-34
    [68] Goldstein E,Gotsman C,Polygon morphing using a multiresolutionrepresentation[C].Proceedings of Graphics Interface95.Paris:[s.n.],1995:247-254.
    [69] Shapira M,Rappoport A. Shape blending using the star skeletonrepresentation[J]. IEEE Transactions on Computer Graphics andApplication,1995,15(3):44-51.
    [70] Carmel E,Cohen-Or D.Warp-guided object-space morphing[J].TheVisual Computer,1997,13:465-578.
    [71]刘文予,万菲,朱光喜.一种全局优化的多边形变形方法及应用[J].计算机辅助设计与图形学学报,1999,11(2):129-133.
    [72] Alexa M,CohenOr D,Levin D.As-rigid-as possible shapeinterpolation[C].Proceedings of SIGGRAPH'2000,New York:ACMPress,2000:157-164.
    [73] Surazhsky V,Gotsman C.Guaranteed intersection-free polygonmorphing[J].Computer and Graphics,2001,25(1):67-75.
    [74] Floater M S,Gotsman C.How to Morph Tilings Injectively[J].Journalof Computational and Applied Mathematics,1999,101:117-129.
    [75]谷留新,刘克轩.改进的基于mean value重心坐标的多边形变形[J].计算机工程与应用,2005,(29):74-76.
    [76]宋伟杰,蒋大为,叶正麟.平面多边形的离散曲率插值变形[J].工程图学学报,2005,(1):50-56.
    [77]章银娥,车建仁.改进的向量多边形变形方法研究[J].大地测量与地球动力学,2009,29(5):94-98.
    [78]许婷,刘占一,胡玫瑰.可控制形状的多边形变形算法[J].计算机工程与应用,2010,46(36):216-218.
    [79]刘婕,基于小波变换的平面多边形变形算法[J].西安邮电学院学报,2011,16(2):65-67.
    [80]杨梦云,刘雁,徐丹.交互式的多边形变形控制技术[J].云南大学学报(自然科学版),2011,33(S2):225-229.
    [81]王刚,雷博琪,叶正寅.一种基于径向基函数的非结构混合网格变形技术[J].西北工业大学学报,2011,29(5):783-788.
    [82]任绍忠,刘利刚,王国瑾.保特征形状过渡的伸缩内在量算法[J].中国图象图形学报,2006,(16):147-154.
    [83] Floater MS.Mean vlaue coordinates[J].Computer Aided GeometryDesign,2003,(20):19-27.
    [84]张冬梅,刘利刚.多边形高质量同构三角剖分的有效算法.浙江大学学报(工学版),2008,(5):1-8.
    [85] Liang Wu.Report on Optimal Compatible Triangulations of SimplePolygons.Mar08,2005.
    [86] Barequet G,Sharir M. Piecewise-linear Interpolation betweenPolygonal Slices. Computer Vision Image Understanding,1996,63:251~272.
    [87]陈敏,鲍旭东.由任意形状轮廓线重建三维表面的方法研究[J].计算机工程与应用,2006.12,74-77.
    [88]马洪滨,郭甲腾.一种新的多轮廓线重构三维形体算法:切开-缝合法.东北大学学报.2007,28(1):111-114.
    [89]李翠平;李仲学;郝秀强;郝晋会.地矿工程可视化仿真中品位与储量的计算实现[J].北京科技大学学报,2007,29(9):859-862.
    [90]谢经平,朱怀球,吴泰然,潘懋.基于Voronoi Cells的自然邻点插值法在地形可视化中的应用[J].遥感信息,2002,(1):15-18.
    [91]侯景儒,黄竞先.实用地质统计学[M].北京:地质出版社,1998.
    [92]李章林,张夏林.距离平方反比法矿产资源储量计算模块设计与实现[J].地质与勘探,2007,43(6):92-97.
    [93]王靖波,潘懋,张绪.基于Kriging方法的空间散乱点插值[J].计算机辅助设计与图形学学报,1999,11(6):525-529.
    [94]曲鹏举,邓寅生,翟东兴,庞玉娟.基于Kriging算法的三棱柱地质体三维可视化实现[J].华北国土资源,2010,(1):52-56.
    [95]荆永滨,王李管,毕林,陈建宏.复杂矿体的块段模型建模算法[J].华中科技大学学报(自然科学版),2010,38(2):97-100.
    [96] M Lin,Gottschalk S.Collision Detection between GeometricModels:ASurvey[C].The Proceedings of MA Conference on Mathematics ofSurfaces,Birmingham,UK,1998,Birmingham:EurographicsAssociation,1998.
    [97]孙晓光,王明强.碰撞检测中的层次包围盒算法研究[J].现代制造工程,2009,(4):87-91.
    [98] Gott schalk S,Lin M C,Manocha D. OBBTree:a hierarchical st ructurefor rapid interference detection[C].Proceedings of the23rd AnnualConference on Computer Graphics and Interactive Techniques. NewYork:ACM,1996:171-180.
    [99] Gargantini H H A. Ray t racing an oct ree:evaluation of the firstintersection [J].Graphics Forum,1993,12(4):199-210.
    [100]江方记.基于空间分区二叉树实现多边形布尔运算的算法研究[J].北京电子科技学院学报,2008,16(2):80-83.
    [101] J ames A. Binary space partitioning for accelerated hidden surfaceremoval and rendering of static environments[D].School of InformationSystems,University of East Anglia,1999.
    [102] Ogayar C J,Segura R J,Feito F R.Point in solidstrategies[J].Computers&Graphics,2005,29(4):616-624.
    [103] Feito F R,Torres J C.Inclusion test for generalpolyhedra[J].Computers&Graphics,1997,21(1):23-30.
    [104]阎少宏,孙丽萍,纪楠,杨爱民.离散曲面上的单位法向量与平均曲率的计算.河北理工大学学报(自然科学版),2008,30(1):104-109.
    [105]僧德文,李仲学,李翠平,春民.Marching Cubes算法改进研究及应用[J].计算机应用研究,2006,(7):50-51,61.
    [106]邹耀斌.Gouraud明暗处理的探索[J].图形图形处理与游戏编程,2006,(9):61-67.
    [107] Clark J.Hierarchical Geometric Models for Visible SurfaceAlgorithms[J].Communications of the ACM,1976:547-554.
    [108]彭群生.计算机真实感图形的算法基础[D].科学出版社,2002.
    [109]张宛方,苏鸿根.3D模型LOD算法的研究及其OpenGL实现[J].计算机工程与设计,2004,25(10):1790-1796.
    [110]芮小平,张彦敏.一种实时连续LOD技术的改进算法[J].系统仿真学报,2004,16(11):2628-2630.
    [111] Shostko A A, Lohner R, Sandberg W C. Surface Triangulation overIntersecting Geometries[J]. International Journal for NumericalMethods in Engineering,1999,44(9):1359-1376.
    [112] Coelho L C G, Gattass M, Fiqueiredo L H. Intersecting and TrimmingParametric Meshes on Finite-elements Shells[J]. InternationalJournal for Numerical Methods in Engineering,1999,47(4):777-800.
    [113] Gottschalk S, Lin M C, Manocha D. OBBTree: A Hierarchical Structurefor Rapid Interference Detection[C]//Proceedings of ACM SIGGRAPH’99,New York.1999.
    [114]陈国良,刘修国,尚建嘎,花卫华.三维地质结构模型的切割分析技术及方法[J].计算机工程,2007,33(20):184-186.
    [115]苗小平,杨永国,奚砚涛.一种基于三棱柱的三维地质体可视化方法研究[J].中国矿业大学学报,2004,33(5):584-588.
    [116]张煜,白世伟.一种基于三棱柱体体元的三维地层建模方法及应用[J].中国图象图形学报,2001,6(3):285-290.
    [117]刘惠义.虚拟视景交互漫游过程中视点的运动控制方法[J].计算机辅助设计与图形学学报,2004,32(6):693-695.
    [118]侯守明,朱雅莉.基于OpenGL的虚拟漫游技术研究[J].图形图像,2004,257:43-45.
    [119]熊伟,毛善君,马蔼乃.基于观察者的碰撞检测技术在虚拟环境漫游中的应用[J].计算机应用,2002:7-10.
    [120]袁崎.储量计算中平行断面法公式的使用条件[J].地质与勘探,1984,
    [121] Matheron G. Principles of geostatistics. EconomicGeology,1963,58(8):1246-1266.
    [122]金毅,潘懋,姚凌青,等.一种稳健变差函数计算方法.北京大学学报:自然科学版,2009,45(6):1033-1038.
    [123] Journel A G, Huijbregts C H J. Mining geostatistics.London:Academic Press,1978.
    [124]孙英君,王劲峰,柏延臣.地质统计学方法进展研究[J].地球科学进展,2004,19(2):265-270.
    [125]谢锡林.地质统计学及其在储量计算中的应用[M].北京:科学技术出版社,1988.
    [126]朱裕生.矿产资源评价方法学导论[M].地质出版社,1984.
    [127]柳炳利.基于SD法的固体矿产资源研究[D].成都理工大学硕士论文,2009.

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700