摘要
三维地质模拟(3DGM)是当前三维地理信息系统(3DGIS)和数字矿山(DM)理论与技术发展急需研究和解决的关键问题之一,其研究成果对基础地质研究、矿产资源勘查与开发等具有重要意义。
本文以地理信息系统理论为指导,以点集拓扑理论为基础,采用面向对象的方法研究了基于剖面的三维拓扑(OO-Solid)数据模型和数据结构,采用计算几何、地理信息系统、计算机图形图像学等方法研究了基于修正后的OO-Solid数据模型的三维地质模型构建与动态更新方法。在数据模型和数据结构研究方面,针对OO-Solid数据模型的不足,修正了OO-Solid数据模型及其构模元素,对修正后的OO-Solid数据模型进行了形式化描述,重新设计了相应的拓扑关系和数据结构。在基于修正后的OO-Solid数据模型的三维地质模型的构建方法研究方面,主要研究了三维地质建模的流程、空间剖面的定义与生成方法、体元的生成方法、三维地质模型的生成方法和关键算法,提出了“剖面构模元素拓扑关系自动构建”的新算法、“已知前后剖面推断中剖面”的新思路和构建复杂地质体的“分而构之”的新方法。在基于修正后的OO-Solid数据模型的三维地质模型的动态更新方法研究方面,提出了“由上至下”进行动态更新的思路和点、线、面、体四类对象动态更新的方法。在上述研究的基础上,以SQL Server2000、VC++6.0、OPENGL等为工具,进行了实验系统的设计与开发,并利用河北金牛能源股份有限公司章村矿地质勘探和生产数据对所开发的实验系统进行了应用研究,验证了上述几方面研究成果的实用性和有效性。
Three-Dimensional Geological Modeling (3DGM) is one of the key issues of Three-Dimensional Geological Information System (3DGIS) and Data Mine. It is significant that its’research results apply in the region of basic geological research and resource exploration and development.
The Geographic Information System (GIS) theory as a guide, point-set topology theory as basis, this paper uses object-oriented method to study on profile-based three-dimensional topological data model and data structure. By use of computational geometry, GIS, computer graphics images, 3DGM, the construction and dynamic updating methods of 3DGM are researched on the basis of amended OO-Solid data model. In research on data model and data structure, according to the limitations of modeling, data model and its elements of OO-Solid was revised. After revision, the revised OO-Solid data model was formal described, and the topological relationship and data structure was redesigned. Based on the revised OO-Solid model, the paper researched the construction methods of three dimensional geological models, includes the process of 3DGM the definition of spatial profile and generated methods, the generation methods of body element, generation methods and key algorithms of 3DGM. The new algorithm which the auto construction about topological relationship of profile modeling elements are automatic constructed, was proposed. The new idea, which was inferred the mid-profile on the basis of known fore-profile and back-profile, was given. The new method, which was divided the complex geological model and constructed later, was put forward. Based on revised OO-Solid data model in the dynamic updating methods of three-dimensional geological model, the new idea, which was updating from up to bottom, and the dynamic updating method about point, line, surface and body are proposed. On the basis of above study, SQL Server2000、VC++6.0、OPENGL as tools, the experimental system was designed and developed. It was shown that the paper verify the practicality and effectiveness of three-dimensional data model, three-dimensional modeling methods and algorithm for dynamic update by using examples of the actual coal mine exploration and production data from ZhangCun mine in HeBei Jinniu energy company limited.
引文
[1]陈述彭,鲁学军,周成虎.地理信息系统导论.北京:科学出版社,1999.
[2]李德仁.当前国际GIS的研究和应用现状, RS,GIS,GPS的集成与应用[J].北京:测绘出版社.
[3]龚健雅.当代地理信息系统进展综述[J].测绘与空间地理信息,2004,27(1):20-22.
[4]程朋根.地矿三维空间数据模型及相关算法研究[D].武汉:武汉大学博士学位论文, 2005.
[5]车德福.基于GTP的复杂地质体多尺度空间建模研究[D].北京:中国矿业大学博士学位论文, 2006.
[6]吴立新,张瑞新,戚宜欣,等.三维地学模拟与虚拟矿山系统[J] .测绘学报. 2002, 31(1): 28-33.
[7]吴立新,殷作如,钟亚平,等.再论数字矿山:特征、框架与关键技术[J].煤炭学报, 2003, 28(1)1-7.
[8]陈建宏.可视化集成采矿CAD系统研究[D].中南大学博士学位论文, 2002.
[9] Houlding S. W. 3D Geoscience Modeling: Computer Techniques for Geological Characterization [M]. Springer-Verlag, 1994, 27-36.
[10] Robert laurini and Derek. Thompson. Fundamentals of Spatial Information Systems[M]. Academic Press INC. 1992.
[11]侯恩科.三维地质模拟的若干理论、方法与应用问题研究[D].北京:中国矿业大学博士学位论文, 2002.
[12] Goodchild M F.Spatial Information Science[C]. In Proceedings of the Spatial Data Handling Symposium, International Geographical Union, Ohio, 1991, 3-14.
[13]郭薇.顾及空间剖分的三维拓扑空间数据模型[D].武汉:武汉测绘科技大学博士学位论文,1998.
[14] Kelk B. 3-D Modelling with geoscientific information systems[C]: The problem //Turner A K. Three-dimensional Modeling with Geoscientific Information System. Dordrecht: Kluwer Academic Publishers, 1992: 29-37.
[15] Christopher B. J. Data Structure for 3D Spatial Information System in Geology[J]. Int. J. GIS, 1989, 3(1): 15-31.
[16] Abbots F. V., Kuijk A. D. Using 3D geological modeling and connectivity analysis to locate remaining oil targets in the Brent reservoir of the mature Brent Field[C]. Offshore Europe Conference Proceedings. Aberdeen, UK: SPE, 1997, 107-114.
[17]王笑海.基于三维拓扑格网结构的GIS地层模型研究[D].中科院武汉岩土所博士学位论文, 1999.
[18]程朋根,龚健雅.地勘工程3维空间数据模型及其数据结构设计[J].测绘学报, 2001, 30(1): 74-81.
[19]张煜,白世伟.一种基于三棱柱体体元的三维地层建模方法及应用[J] .中国图象图形学报, 2001, 6(3): 83-88.
[20]齐安文,吴立新,李冰,等.一种新的三维地学空间构模方法——类三棱柱法[J] .煤炭学报, 2002, 27(2): 158-163.
[21]侯恩科,吴立新.面向地质建模的三维体元拓扑数据模型研究[J] .武汉大学学报(信息科学版): 2002, 27(5): 467-472.
[22]王勇.三维空间矢量数据生成算法研究及实验系统实现[D].北京:北京大学博士论文, 2003.
[23]武强,徐华.三维地质建模与可视化研究[J].中国科学(D辑)地球科学,2004,34(1): 54-56.
[24] Zhong D. H., Li M. C., Song L. G. , et al. Enhanced NURBS modeling and visualization for large 3D geoengineering applications: An example from the Jinping first-level hydropower engineering project, China[J]. Computers & Geosciences, 2006, 32(9): 1270-1282.
[25]吴立新,史文中, Christoher Gold. 3DGIS与2DGMS中的空间构模技术[J].地理与地理信息科学, 2003, 19(1):5-11.
[26]邓念东,侯恩科,张志华,等.三维拓扑关系形式化描述及拓扑关系模型研究[J].西安建筑科技大学学报(自然科学版) ,2007,39(6): 873-877.
[27]邓念东,侯恩科.三维体元拓扑数据模型的修正及其形式化描述[J].武汉大学学报(信息科学版). 2009, 34(1): 52-56 .
[28]郭薇.顾及空间剖分的三维拓扑空间数据模型[D].武汉:武汉测绘科技大学博士学位论文,1998.
[29]邓念东,侯恩科.基于GTP单纯剖分的地下实体拓扑关系形式化描述方法[J].煤炭学报. 2008, 33(5): 527-531.
[30] Deng nian-dong, hou en-ke, zhang zhi-hua. An improvement of OO-Solid model[C]. proceedings of the international archives of the photogrammetry, remotote sensing and spatial information sciences. Beijing:ISPRS2008. volume37, part B4-1, 343-378.
[31] Hou en-ke , Deng nian-dong. Research on dynamic updating of three dimensional geological modeling based on OO-Solid model. Coal Science & Engineering(China). 2008,14(3):420-424.
[32] Christopher C M, Maydell U M. Triangulation and spatial ordering in computer cartography[C]. Proceedings of 3rd Canadian Cartographic Association Annual Meeting. Canada, 1978. 69-8l.
[33] Nicolas S , Renato P. Delaunay Triangulation of Arbitrarily Shaped Planar Domains[J]. Computer Aided Geometric Design. 1991, 8: 421-437.
[34] Victor J D. Delaunay Triangulation in TIN Creation: A Overview and A Linear - Time Algorithm[J]. GIS. 1993, 7(6): 501-524.
[35] S.W.赫尔丁,龙子芳.三维矿床的计算机构模方法[J].矿业工程, 1989,(03): 95-96
[36]阿列尼切夫B M,弗拉基米罗夫A N.菱镁矿股份公司露天采矿工程计划编制自动化[J].国外金属矿山. 1995, (12): 70-74.
[37]赵树贤.煤矿床可视化构模技术[D] .中国矿业大学(北京校区)博士学位论文, 1999.
[38] Lemon A M, Jones N L.Building solid models from boreholes and user - defined cross - sections [J].Computers&Geosciences, 2003, 29:547- 555.
[39]朱良峰,吴信才,刘修国,等.基于钻孔数据的三维地层模型的构建[J].地理与地理信息科学. 2004, 20(3): 26-30.
[40]吴冲龙,毛小平,田宜平,等.三维数字盆地构造-地层格架模拟技术[J].地质科技情报. 2006, 25(4): 1-8.
[41]程朋根,王承瑞,甘卫军,等.基于多层DEM与QTPV的混合数据模型及其在地质建模中的应用[J].吉林大学学报(地球科学版). 2005, (6): 128-133.
[42]李青元.三维矢量结构GIS拓扑关系及其动态建立[J].测绘学报. 1997, 26(3): 235-240.
[43]陈军,郭薇.基于剖分的三维拓扑ER模型研究[J].测绘学报,1998,27(4):308-317.
[44]姜在炳.煤层动态建模技术及应用[J].煤炭学报, 2006, 31(1): 40-44.
[45] Homer H, Thomas S. A Survey of Construction and Manipulation of Octree[J]. CVGIP, 1998, 43: 409-431.
[46]边馥苓,傅仲良,胡自锋.面向目标的栅格矢量一体化三维数据模型[J].武汉测绘科技大学学报, 2000, 25(4): 294-298.
[47]韩国建,郭达志,金学林.矿体信息的八叉树存储和检索技术[J].测绘学报, 1992, 21(1): 13-17.
[48]肖乐斌,龚建华,谢传节.线性四叉树和线性八叉树邻域寻找的一种新算法[J].测绘学报, 1998, 27(3): 195-203.
[49]李清泉,李德仁.八叉树的三维行程编码[J].武汉测绘科技大学学报, 1997, 22(2): 12-16.
[50] Joe B. Construction of 3D Delaunay Triangulations Using Local Transformations[J].Computer Aided Geometric Design, 1991, 8: 123-142.
[51] Morackot P, Tempfli K, Molenaar M. A Tetrahedron-Based 3D Vector Data Model for Geoinformation[C]. In: Advanced Geographic Data Modelling(Ed. , Molenaar M). Sylvia De Hoop: Geodesy Press, 1994. 129-140.
[52] Victor J D , Alan P. Delaunay tetrahedral Data Modeling for 3D GIS Application[C]. Proc. GIS/LIS'93, 1993, 671-678.
[53] Shi W Z. Development of A Hybrid Model for Three - Dimensional GIS[J] . Geo-Spatial Information Science, 2000, 3(2): 6-12.
[54]张煜,白世伟.一种基于三棱柱体体元的三维地层建模方法及应用[J] .中国图象图形学报, 2001, 6(3): 83-88.
[55]戴吾蛟,邹峥嵘.基于体素的三维GIS数据模型的研究[J] .矿山测量, 2001, (1): 20-22.
[56] Gong J Y, Cheng P G. Study on 3D modeling and Visualization in Geological Exploration Engineering[C]. Proc. of ISPRS CommissionⅡSymposium: Integrated System for Spatial Data Production, Custodian and Decision Support. Xi'an, China, 2002, pp133-138.
[57]毛善君.灰色地理信息系统——动态修正地质空间数据理论与技术[J].北京大学学报(自然科学版) , 2002, 38(4): 556-562.
[58]侯恩科,赵洲.三维体元拓扑数据模型的改进与实验[J].煤田地质与勘探, 2006, 34(4): 13-16.
[59]赵洲,张恒,侯恩科,等.三维地质界面模型交互式选择方法研究与实现[J].西安科技大学学报, 2008, 28(4): 98-101.
[60] Shi W Z. A Hybrid Model for 3D GIS[J]. Geoinformatics, 1996(1): 400-409.
[61]唐泽圣.三维数据场可视化.北京:清华大学出版社,1999.
[62]李梅,毛善君,马蔼乃.平行轮廓线三维重建算法[J].计算机辅助设计与图形学学报. 2006, 18(7): 101-102.
[63]陈学工.三维矢量数据生成方法研究和实现[D].北京大学博士论文, 2002.
[64]屈红刚,潘懋,王勇,等.基于含拓扑剖面的三维地质建模[J] .北京大学学报(自然科学版), 2006, 42(06): 717-723.
[65]朱良峰,潘信,吴信才.三维地质建模及可视化系统的设计与开发[J].岩土力学, 2006, 27(5): 828-832.
[66] Raper J. F. Key 3D Modeling Concepts for Geoscientific Analysis. In: Turner A. K., ed., Three-Dimensional Modeling with Geoscientific Information Systems. Kluwer Academic Publishers. Dordrechts, Netherlands. 1992, 215～232.
[67]李培军.面向层状地质实体的三维地学信息系统研究[博士后研究工作报告].北京:北京大学,1997.
[68] Kell B. 3-D Modeling Geology in Three Dimension. Geobyte. 1988, 3 (1): 14～20.
[69] Rhind D. W. Spatial Data Handling In the Geosciences. In: Turner A. K., ed., Three-Dimensional Modeling with Geoscientific Information Systems. Dordrechts, Netherlands: Kluwer Academic Publishers, 1992, 13～28.
[70] Bonham-Carter G. F. Geographic Information Systems for Geoscientists: Modeling with Gis. Pergamon. 1994, 25～50.
[71] Bruce R. M. Report of Discussion by the Minerals Exploration Interest Group. In: Turner A. K., ed., Three-Dimensional Modeling with Geoscientific Information Systems. Dordrechts, Netherlands: Kluwer Academic Publishers, 1992, 391～396.
[72]朱志澄.构造地质学[M].中国地质大学出版社.1999.
[73]王润怀.矿山地质对象三维数据模型研究[D].西南交通大学博士学位论文, 2007.
[74] Frank A U.The National Center for Geographic Information and Analysis in the U.S.A, FIGXIX, 1990.50-64
[75] Golledge, Reginald G. Do People Understand Spatial Concepts:the Case of First-order Primitives.In: Lecture Notes in Computer Science 639.Pisa:Springer-Verlag,1992.1-21
[76] Egenhofer M, Franzosa R. On the Equivalence of Topological Relations [J]. International Journal of Geographic Information Systems , 1995 , 9(2):133～152.
[77]郭薇,陈军.基于点集拓扑学的三维拓扑空间关系形式化描述.测绘学报, 1997, 26(2): 122-127.
[78] Guo Wei, Zhan Ping, Chen Jun. Topological data model for 3D GIS[J]. APRS, 1998, 32(4): 657-661.
[79] Goodchild M F. Towards an Enumeration and Classification of GIS Functions[C]. In: Int. GIS Symposium: The research agenda, II67-II77, Arlington, VA: NASA,1987.
[80]毋河海,龚健雅. GIS空间数据结构与处理技术.测绘出版社, 1997.
[81] Molenaar M. An Introduction to the Theory of Spatial Object Modeling[C]. London: Taylor & Francis, 1998.
[82] Florence III, John and Egenhofer, Max J. Distribution of Topological Relations in Geographic Datasets[C]. In: ASPRS/ACSM. Annual Convention and Exposition Technical Papers. 1996. 315-325.
[83] Gold Christopher M. The Meaning of“Neighbor”[J]. Lecture Notes in Computer Science. Pisa: Springer-Verlag, 1992. 220-235.
[84]胡勇,陈军.基于Voronoi图的空间邻近关系表达和查询操作[C].中国GIS协会第二届年会论文集. 1997, 346-356.
[85] Zhao Ren liang, Chen Jun , Li Zhilin. Voronoi-based Generalized Spatial Adjacency[C]. In: Li De ren, et al., ed. The proceeding of RS, GPS, GIS, The integration and applications. Wuhan Technical University of Surveying and Mapping Press. 1998. 605~614.
[86]郭庆胜,杜晓初.不确定性线状目标之间拓扑关系的描述与判别[J].武汉大学学报(信息科学版), 2004, 29(09), 827-831.
[87]郭庆胜,杜晓初.模糊面元素空间拓扑关系抽象化方法研究[J].测绘学报, 2004, 33(04), 307-310.
[88]虞强源,刘大有,欧阳继红.基于区间值模糊集的模糊区域拓扑关系模型[J].电子学报, 2005, 33(01), 186-189.
[89]虞强源,刘大有,王生生.模糊栅格区域的层次拓扑关系模型[J].计算机研究与发展, 2005, 42(07), 1166-1172.
[90]高振记,高勇,张毅,等.基于粗糙模型的模糊地理对象空间拓扑关系研究[J].地理与地理信息科学, 2006, 22(06), 12-16.
[91]李大军,刘波,刘学军,等. GIS中模糊区域拓扑关系形式化描述[J].桂林工学院学报, 2006, 26(04), 518-520.
[92]李大军,刘波,程朋根,等.模糊空间对象拓扑关系的Rough描述[J].测绘学报, 2007, 36(01), 72-77.
[93] Egenhofer M. Franzosa R. Point-Set Topological Spatial Relations[C]. International Journal of Geographical Information Systems , 1991 , 5(2): 161～174.
[94] Egenhofer M. Reasoning about binary topological relations [C] //Proceedings of the 2nd International Symposium on Advances in Spatial Databases, Zurich, 1991: 143-160.
[95] Clementini E, Di F P, Oosterom P. A Small Set for Formal Topological Relationships Suitable for End-user Interaction [C]. In: David A, Beng C O, eds. Advances in Spatial Databases. New York: Springer-Verlag, 1993.
[96] Clementini E, Di Felice P. Topological Invariants for lines [J]. IEEE Transactions on Knowledge and Data Engineering , 1998 , 10 :38～54.
[97] Nian-dong Deng, En-ke Hou, Zhi-hua Zhang, et al. Formal Description and Modeling Topological Relations in 3D Based on Dimension-Extended and Euler-Poincare Characteristics[C]. 5th WSEAS Int. Conference on Environment, Ecosystems and Development.Spain. 2007: 29-33.
[98]邓敏,张雪松,林宗坚.拓扑关系形式化描述的Euler示性数模型[J].武汉大学学报(信息科学版), 2004, 29(10), 872-876.
[99]齐华.自动建立多边形拓扑关系算法步骤的优化与改进[J].测绘学报, 1997, 26(3), 254-260.
[100]杨海宏,李天宏,薛安.利用双邻点判断法优化拓扑关系自动生成算法[J].中国图象图形学报, 1998, 13(7), 583-587.
[101]闫浩文,杨维芳,陈全功,等.基于方位角计算的拓扑多边形自动构建快速算法[J].中国图象图形学报, 2000, 5(A)(7), 563-568.
[102]梁晓文,刘宗岐,陈宜金.基于夹角变化趋势的多边形自动搜索和生成算法[J].中国图象图形学报, 2005, 10(6), 785-789.
[103]李大军,刘波,赵宝贵,等.拓扑多边形自动构建的一种改进算法[J].计算机工程与应用, 2005, 42(16), 80-82.
[104]丁永祥,夏巨湛,王英,等.任意多边形的delaunay三角剖分[J].计算机学报, 1994, 17(4), 270-275.
[105]周晓云,刘慎权.实现约束delaunay三角剖分的健壮算法[J].计算机学报, 1996, 19(8): 615-624.
[106]VigoM, PlaN. Computing directional on constrained delaunay traingulations[J]. Computers & Graphics, 2000, 24: 181-190.
[107]关振群,宋超,顾元宪,等.有限元网格生成方法研究的新进展[J].计算机辅助设计与图形学学报, 2003, 15(1): 1-14.
[108]蔡强,杨钦,陈其明,等.地质结构重叠域的限定Delaunay三角剖分研究[J].计算机辅助设计与图形学学报, 2004, 16( 6) : 766- 771.
[109]王占刚,潘懋,屈红刚,等.三维折剖面的Delaunay三角剖分算法[J].地理与地理信息科学, 2008, 23(3):37-40.
[110]Piegl L A. Algorithm and Data Structure for Triangulation Multiply Connected Polygonal Domains. Computer and Graphics, 1993, 17(5): 563～574.
[111]Lee D T. Generalized Delaunay Triangulation for Planar Graphs. Discrete and Computational Geometry, 1986, (1): 201～217.
[112]Chew L P. Constrained Delaunay Triangulation, in Proceedings of Third ACM Symposium on Computational Geometry. Waterloo, 1987, 216～222.
[113]Boissinnat J D. Shape Reconstruction from Planar Sections. Comput Vision Graphics Image Process, 1988, 44: 1～29.
[114]刘学军,龚健雅.约束数据域的Delaunay三角剖分与修改算法.测绘学报,2001,30(1) :82～88.
[115]L De Floriani, B Falcidieno, C Pienovi. Delaunay based representation of surface defined over arbitrarily shaped domains[J]. Computer Vision, Graphics and ImageProcessing, 1985, 32:127-140.
[116]Schweitzer D, Cobb E S. Scanline rendering of parametric surfaces[J].Computer Graphics 1982,16(3):265-274.
[117]徐松,王剑英.曲面的自适应三角网格剖分[J].计算机辅助设计与图形学学报, 2000, 12(4): 4-8.
[118]史国友,贾传荧.基于NURBS的船体曲面自适应三角网格剖分[J].交通运输工程学报, 2006, 6(1): 84-88.