铜陵冬瓜山铜矿床的三维形态及成矿动力学计算模拟
|
摘要
成矿系统是一个复杂的时空体系,计算模拟是揭示其空间形态和结构特征以及时间进程动力学的最有效手段。
铜陵冬瓜山铜矿床是长江中下游成矿带最大的铜矿床,本文选择冬瓜山铜矿床为对象进行三维形态建模及成矿动力学计算模拟研究。以GOCAD为建模工具,通过离散光滑插值和Delaunay三角剖分等建立矿体及其它相关地质体的三维形态模型。通过编程实现GOCAD的SGrid模型向FLAC3D的六面体网格模型的转换,在地质研究的基础上建立三维动力学模型进行力-热-流体的耦合动力学计算模拟实验,根据计算模拟实验的结果分析成矿的动力学机制以及控矿因素。
1)通过三维建模及可视化,展示了矿体的三维形态及其与其他地质要素之间的空间关系,发现了矿体受背斜构造和岩体的双重控制。
2)系统总结了三维网格剖分通用算法,分析了网格剖分技术在地学领域的研究现状和发展前景。
3)编制了从GOCAD的SGrid网格转化成FLAC3D的六面体剖分网格接口程序,实现了复杂实体地质模型向FLAC3D模型的转化,并将此方法成功地应用于冬瓜山铜矿床的动力学建模。
4)动力学计算模拟结果显示,冬瓜山矿床的形成与岩体同构造冷却过程耦合动力学作用密切相关,力-热-流耦合作用在中-上石炭统的碳酸盐岩中形成汇流扩容带,为层控矿体的形成创造了有利的条件。
The metallogenic system is complex in spatial structure and temporal process. The computational modeling is the most effective method for revealing its shape and structure and replaying its geodynamic processes. In this thesis, the Dongguashan copper deposit, as the largest copper deposit in the Yangtze River Metallogenic Belt, was studied through 3D shape simulation and computational modeling of its metallogenic processes. By using of the GOCAD code,3D geological models are constructed through discrete smooth interpolation and triangulated irregular network. By programming, the SGrid model of GOCAD is transferred to the hexahedron-mesh model of FLAC3D. Based on the geological investigation, the 3D geodynamic model involved the coupled mechano-thermo-hydrological processes is constructed and computed by numerical method, of which the results are analyzed to extract the implication for dynamic mechanism of ore-formation and ore-controlling factors. The main outcomes in this thesis include the following 4 aspects.
1) The 3D shape of the main ore body and the spatial relationship of the main ore body with other geological factors are revealed by 3D geological modeling, and it is discovered that the main ore body is controlled by both the anticline and the intrusion.
2) By systematically summarizing the algorithms of 3D mesh generation, the current situation and development trend of the mesh generation and its application in the geosciences are analyzed.
3) The complex 3D geological model is transferred to the FLAC3D model through programming the transferring of the SGrid model of GOCAD into the hexahedron-mesh model of FLAC3D. It is successfully used to model the geodynamic processes of the metallogenic system in the Dongguashan copper deposit.
4) The results of the geodynamic computational modeling experiments demonstrate that, the ore formation is closely associated with coupled geodynamic processes during the syn-tectonic cooling of the intrusion, and the coupled mechano-thermo-hydrological processes result in the generation of the fluid-focus dilation zone in the carbonate of the medium to upper Carboniferous, favorable for forming the stratum-bound ore body.
铜陵冬瓜山铜矿床是长江中下游成矿带最大的铜矿床,本文选择冬瓜山铜矿床为对象进行三维形态建模及成矿动力学计算模拟研究。以GOCAD为建模工具,通过离散光滑插值和Delaunay三角剖分等建立矿体及其它相关地质体的三维形态模型。通过编程实现GOCAD的SGrid模型向FLAC3D的六面体网格模型的转换,在地质研究的基础上建立三维动力学模型进行力-热-流体的耦合动力学计算模拟实验,根据计算模拟实验的结果分析成矿的动力学机制以及控矿因素。
1)通过三维建模及可视化,展示了矿体的三维形态及其与其他地质要素之间的空间关系,发现了矿体受背斜构造和岩体的双重控制。
2)系统总结了三维网格剖分通用算法,分析了网格剖分技术在地学领域的研究现状和发展前景。
3)编制了从GOCAD的SGrid网格转化成FLAC3D的六面体剖分网格接口程序,实现了复杂实体地质模型向FLAC3D模型的转化,并将此方法成功地应用于冬瓜山铜矿床的动力学建模。
4)动力学计算模拟结果显示,冬瓜山矿床的形成与岩体同构造冷却过程耦合动力学作用密切相关,力-热-流耦合作用在中-上石炭统的碳酸盐岩中形成汇流扩容带,为层控矿体的形成创造了有利的条件。
The metallogenic system is complex in spatial structure and temporal process. The computational modeling is the most effective method for revealing its shape and structure and replaying its geodynamic processes. In this thesis, the Dongguashan copper deposit, as the largest copper deposit in the Yangtze River Metallogenic Belt, was studied through 3D shape simulation and computational modeling of its metallogenic processes. By using of the GOCAD code,3D geological models are constructed through discrete smooth interpolation and triangulated irregular network. By programming, the SGrid model of GOCAD is transferred to the hexahedron-mesh model of FLAC3D. Based on the geological investigation, the 3D geodynamic model involved the coupled mechano-thermo-hydrological processes is constructed and computed by numerical method, of which the results are analyzed to extract the implication for dynamic mechanism of ore-formation and ore-controlling factors. The main outcomes in this thesis include the following 4 aspects.
1) The 3D shape of the main ore body and the spatial relationship of the main ore body with other geological factors are revealed by 3D geological modeling, and it is discovered that the main ore body is controlled by both the anticline and the intrusion.
2) By systematically summarizing the algorithms of 3D mesh generation, the current situation and development trend of the mesh generation and its application in the geosciences are analyzed.
3) The complex 3D geological model is transferred to the FLAC3D model through programming the transferring of the SGrid model of GOCAD into the hexahedron-mesh model of FLAC3D. It is successfully used to model the geodynamic processes of the metallogenic system in the Dongguashan copper deposit.
4) The results of the geodynamic computational modeling experiments demonstrate that, the ore formation is closely associated with coupled geodynamic processes during the syn-tectonic cooling of the intrusion, and the coupled mechano-thermo-hydrological processes result in the generation of the fluid-focus dilation zone in the carbonate of the medium to upper Carboniferous, favorable for forming the stratum-bound ore body.
引文
[1]於崇文.地质系统的复杂性[M].北京：地质出版社,2003
[2]刘亮明.成矿理论的预测能力及其改善途径[J].地学前缘,2007,14(5)：82-91
[3]Houlding S W.3D geoscience modeling:computer techniques for geological characterization[M]. Berlin Heidelberg:Springer Verlag,1994.
[4]李清泉,杨必胜,史文中,等.三维空间数据的实时获取、建模与可视化[M].武汉：武汉大学出版社,2003
[5]吴立新,史文中.地理信息系统原理与算法[M].北京：科学出版社,2003
[6]Mallet J.L. Discrete smooth interpolation in geometric modeling[J]. Computer Aided Design,1992,24(4):178-190
[7]Mallet J.L. Discrete modeling for nature objects[J]. Mathematical Geology,1997, 29(2):199-218
[8]毛先成.三维数字矿床与隐伏矿体立体定量预测研究[博士学位论文][D].中南大学,2006
[9]朱良峰,潘信,吴信才,等.地质断层三维可视化模型的构建方法与实现技术[J].软件学报,2008,19(8)：2004-2017
[10]曾新平,吴健生,杨自安,等.三维GIS环境下的地质体可视化和特征分析[J].地质与勘探,2005,41(1)：72-76
[11]王纯祥,白世伟,贺怀建.三维地层可视化中地质建模研究[J].岩石力学与工程学报,2003,22(10)：1722-1726
[12]张宝一,尚建嘎,吴鸿敏,等.三维地质建模及可视化技术在固体矿产储量估算中的应用[J].地质与勘探,2007,43(2)：76-81
[13]李培军层状地质体的三维模拟与可视化[J].地学前缘,2000,7(suppl)271-277
[14]吴立新,陈学习,车德福,等.一种基于GTP的地下真3D集成表达的实体模型[J].武汉大学学报(信息科学版),2007,32(4)：331-335
[15]吴立新,陈学习,史文中.基于GTP的地下工程与围岩一体化真三维空间构模[J].地理与地理信息科学,2003,19(6)：l-6
[16]武强,关文革,贾丽萍,等.面向矿区复杂地质体的四面体生成算法[J].中国矿业大学学报,2005,34(5)：617-621
[17]Gocad research group. What is Gocad[EB/OL].2010. http://www.gocad.org /www/gocad/index.xhtml
[18]Dynamic Graphics, Inc. Dynamic Graphics Earthvision Technology[EB/OL]. 2010. http://www.dgi.com/earthvision/evmain.html
[19]Schlumberger limited. GeoFrame Reservoir Characterization Software[EB/OL]. 2010. http://www.slb.com/services/software/geo/geoframe.aspx
[20]徐浩,王庆飞,杨立强,等.成矿多过程耦合模型数值模拟研究态势[J].矿产与地质,2006,20(4-5)：494-497
[21]刘亮明,吕俊武,彭省临,等.成熟勘探矿集区新一轮找矿：勘查战略创新及铜陵矿集区找矿实例[J].地质论评,2005,51(3)：325-333
[22]Hobbs B E, Zhang Y, Ord A, et al. Application of coupled deformation, fluid flow, thermal and chemical modeling to predictive mineral exploration[J]. Journal of Geochemical Exploration,2000,69(70):505-509
[23]Liu L M, Peng S L, Zhang Y H. Numerical geodynamic modeling of coupled mechano-thermo-hydrological processes and its application in predictive exploration in the Fenghuangshan ore field, Tongling[J]. Geotectonica et Metallogenia,2005,29(2):164-173
[24]古成中,吴新跃.有限元网格剖分及发展趋势[J].计算机科学与探索,2008,2(3)：248-259
[25]Cundall P A. A computer model for simulating progressive large scale movement in block rock system[C]. Symposium ISRM,1971, Proc2:129-136
[26]Ransing R S, Gethin D T, Khoei A R, et al. Powder compaction modelling via the discrete and finite element method[J]. Materials and Design,2000,21 (4): 263-269
[27]胥建龙,唐志平.离散元与有限元结合的多尺度方法及其应用[J].计算物理,2003,20(6)：477-482
[28]President's information technology advisory committee. Computational scie nce:ensuring America's competitiveness[EB/OL].2005. http://www.nitrd.gov /pitac/reports/20050609_computational/computational.pdf
[29]赵崇斌,B.E.Hobbs, A.Ord,等.用计算地球科学研究方法探讨地质现象的动力学机制—以断层中等距成矿分布为例[J].中国科学D辑：地球科学,2008,38(5)：646-652
[30]Chongbin Zhao, B.E.Hobbs, H.B.Muhlhaus, et al. Computer simulations of coupled problems in geological and geochemical systems[J]. Computer Methods in Applied Mechanics and Engineering,2002,191(29-30): 3137-3152
[31]刘亮明,疏志明,赵崇斌,等.矽卡岩矿床的汇流扩容空间控矿机制及其对深部找矿的意义：以铜陵-安庆地区为例[J].岩石学报,2008,24(8)：1848-1856
[32]谭凯旋,谢焱石,赵志忠,等.构造-流体-成矿体系的反应-输运-力学耦合模型和动力学模拟[J].地学前缘,2001,8(4)：311-321
[33]邓军,翟裕生,杨立强,等.剪切带构造-流体-成矿系统动力学模拟[J].地学前缘,1999,6(1)：115-127
[34]Itasca Consulting Group, Inc. Fast Lagrangian Analysis of Continua in 3 Dimensions, version 3.0, user's manual[R]. USA:Itasca Consulting Group, Inc FLAC3D,2005
[35]Itasca Consulting Group, Inc. The overview, features, options of FLAC3D [EB/OL].2010. http://www.itascacg.com/flac3d/index.php
[36]ANSYS, inc. The introduction of ANSYS products[EB/OL].2010. http://w ww.ansys.com/products/default.asp
[37]Simulia, inc. The introduction of Abaqus FEA[EB/OL].2010. http://www.s imulia.com/products/abaqus_fea.html
[38]李进文,裴荣富,梅燕雄,等.安徽铜陵狮子山铜(金)矿田成矿流体地球化学研究[J].矿床地质,2006,25(4)：427-437
[39]姜章平,张国培,姚孝德,等.安徽铜陵狮子山矿区冬瓜山铜矿床南段勘探地质报告[R].安徽省地质矿产局321地质队,1995
[40]毛先成,戴塔根.地质体几何造型的计算机实现[J].中南工业大学学报,1997,28(5)：419-422
[41]戴晟.地质体三维可视化研究与系统实现[硕士学位论文][D].华东师范大学,2008
[42]董辉.地质体三维可视化原理与方法研究[硕士学位论文][D].中南大学,2004
[43]易善桢,李琦.3D-GIS数据表示和空间插值方法研究[J].中国图形图像学报,1999,4(8)：698-701
[44]邢延涛.复杂地质岩体三维可视化关键技术研究[硕士学位论文][D].华中科技大学,2006
[45]朱海龙.Gocad项目组简介[J].勘探地球物理进展,2008,31(1)：80-81
[46]董梅,慎乃齐,胡辉,等.基于GOCAD的三维地质模型构建方法[J].桂林工学院学报,2008,28(2)：188-192
[47]Olivier Kaufmann, Thierry Martin.3D geological modeling from boreholes, cross-sections and geological maps, application over former natural gasstorages in coal mines[J]. Computers&Geosciences,2008,34(3):278-290
[48]Andrea Bistacchi, Matteo Massironi, Giorgio V.DalPiaz, et al.3D fold and fault reconstruction with an uncertainty model:An example from an Alpine tunnel case study[J]. Computers&Geosciences,2008,34(4):351-372
[49]Kevin Sprague, Eric de Kemp, WinstonWong, et al. Spatial targeting using queries in a 3-D GIS environment with application to mineral exploration[J]. Computers&Geosciences,2006,32(3):396-418
[50]武晓波,王世新,肖春生.Delaunay三角网的生成算法研究[J].测绘学报,1999,28(1)：28-35
[51]彭仪普.地形三维可视化及其实时绘制技术研究[博士学位论文][D].西南交通大学,2002
[52]刘怀辉.平面区域有限元三角网格剖分算法研究[硕士学位论文][D].山东大学,2007
[53]黄晓东,杜群贵,叶邦彦.三维实体有限元自适应网格规划生成[J].计算机辅助设计与图形学学报,2005,17(7)：1446-1451
[54]金世双.离散数据三角网格剖分研究与实现[J].上海海事大学学报,2005,26(3)：85-87
[55]徐永安,赵静,沈玲玲.约束Delaunay三角网格生成地质构造模型[J].南京师范大学学报,2008,8(4)：155-158
[56]李刚,赵玉新.基于边界特征点提取的约束Delaunay三角剖分算法[J].系统仿真学报,2007,19(16)：3734-3738
[57]刘少兵.带断层地震数据的Delaunay三角剖分算法[J].煤田地质与勘探,2008,36(6)：70-72
[58]陈涛,李光耀.STL格式文件的四边形网格剖分与网格光顺[J].中国机械工程,2009,20(5)：522-528
[59]杨伟军,包忠诩,扶名福,等.映射法在三维六面体有限元网格生成中的应用[J].南昌大学学报(工科版),1999,21(4)：39-43
[60]王东风,翟建军,陈文亮.基于映射法的六面体网格生成算法[J].现代设计与先进制造技术,2009,38(5)：25-27
[61]贾虹,卢炎麟,高发兴.高品质全四边形有限元网格生成的铺砌法[J].浙江工业大学学报,2000,28(4)：254-257
[62]余学进,彭迎风.自动生成有限元四边形网格的铺筑法[J].南昌航空工业学院学报(自然科学版),17(4)：24-26
[63]杜丽惠,江春波.多种介质区域有限元网格自动生成技术[J].水利发电,2000,(2)：55-56
[64]吴宝海,王尚锦.参数曲面网格生成的改进波前法[J].计算机辅助设计与图形学学报,2005,17(8)：1686-1690
[65]Wang Desheng, Wan Shui. Robust Delaunay Tetrahedronal Meshing Coupled with Advancing Front Method[J]. Journal of Southeast University(English Edition),2002,18(2):131-135
[66]杨长伟,孙立镌.结点连元法的一种改进算法[J].哈尔滨理工大学学报,2000,5(5)：97-100
[67]李华,程耿东,顾元宪.一种新的全四边形网格快速生成方法——模板法[J].计算结构力学及其应用,1996,13(1)：25-33
[68]张玉峰,朱以文.有限元网格自动生成的典型方法与研究前瞻[J].武汉大学学报(工学版),2005,38(2)：54-59
[69]徐能雄,武雄,汪小刚,等.基于三维地质建模的复杂构造岩体六面体网格剖分方法[J].岩土工程学报,2006,28(8)：957-961
[70]关振群,宋超,顾元宪,等.有限元网格生成方法研究的新进展[J].计算机辅助设计与图形学学报,2003,15(1)：1-13
[71]吕军,王忠金,王仲仁.有限元六面体网格的典型生成方法及发展趋势[J].哈尔滨工业大学学报,2001,33(4)：485-490
[72]Pinto V, Font X, Salgot M, et al. Using 3D Structures and Their Virtual Representation as a Tool for Restoring Opencast Mines and Quarries[J]. Engineering Geology,2002,63(1-2):121-129
[73]王长海,许国,范孝锋,等,GOCAD与ABAQUS的数据接口实现[J].测绘信息与工程,2008,33(5)：41-42
[74]陈邦国,姜章平,张卫平,等.安徽冬瓜山叠生式层状铜矿热液改造型流体研究[J].江苏地质,2002,26(2)：65-69
[75]黄顺生,徐兆文,倪培,等.安徽铜陵冬瓜山热液叠加改造型铜矿床流体包裹体地球化学特征[J].地质找矿论丛,2003,18(1)：34-38
[2]刘亮明.成矿理论的预测能力及其改善途径[J].地学前缘,2007,14(5)：82-91
[3]Houlding S W.3D geoscience modeling:computer techniques for geological characterization[M]. Berlin Heidelberg:Springer Verlag,1994.
[4]李清泉,杨必胜,史文中,等.三维空间数据的实时获取、建模与可视化[M].武汉：武汉大学出版社,2003
[5]吴立新,史文中.地理信息系统原理与算法[M].北京：科学出版社,2003
[6]Mallet J.L. Discrete smooth interpolation in geometric modeling[J]. Computer Aided Design,1992,24(4):178-190
[7]Mallet J.L. Discrete modeling for nature objects[J]. Mathematical Geology,1997, 29(2):199-218
[8]毛先成.三维数字矿床与隐伏矿体立体定量预测研究[博士学位论文][D].中南大学,2006
[9]朱良峰,潘信,吴信才,等.地质断层三维可视化模型的构建方法与实现技术[J].软件学报,2008,19(8)：2004-2017
[10]曾新平,吴健生,杨自安,等.三维GIS环境下的地质体可视化和特征分析[J].地质与勘探,2005,41(1)：72-76
[11]王纯祥,白世伟,贺怀建.三维地层可视化中地质建模研究[J].岩石力学与工程学报,2003,22(10)：1722-1726
[12]张宝一,尚建嘎,吴鸿敏,等.三维地质建模及可视化技术在固体矿产储量估算中的应用[J].地质与勘探,2007,43(2)：76-81
[13]李培军层状地质体的三维模拟与可视化[J].地学前缘,2000,7(suppl)271-277
[14]吴立新,陈学习,车德福,等.一种基于GTP的地下真3D集成表达的实体模型[J].武汉大学学报(信息科学版),2007,32(4)：331-335
[15]吴立新,陈学习,史文中.基于GTP的地下工程与围岩一体化真三维空间构模[J].地理与地理信息科学,2003,19(6)：l-6
[16]武强,关文革,贾丽萍,等.面向矿区复杂地质体的四面体生成算法[J].中国矿业大学学报,2005,34(5)：617-621
[17]Gocad research group. What is Gocad[EB/OL].2010. http://www.gocad.org /www/gocad/index.xhtml
[18]Dynamic Graphics, Inc. Dynamic Graphics Earthvision Technology[EB/OL]. 2010. http://www.dgi.com/earthvision/evmain.html
[19]Schlumberger limited. GeoFrame Reservoir Characterization Software[EB/OL]. 2010. http://www.slb.com/services/software/geo/geoframe.aspx
[20]徐浩,王庆飞,杨立强,等.成矿多过程耦合模型数值模拟研究态势[J].矿产与地质,2006,20(4-5)：494-497
[21]刘亮明,吕俊武,彭省临,等.成熟勘探矿集区新一轮找矿：勘查战略创新及铜陵矿集区找矿实例[J].地质论评,2005,51(3)：325-333
[22]Hobbs B E, Zhang Y, Ord A, et al. Application of coupled deformation, fluid flow, thermal and chemical modeling to predictive mineral exploration[J]. Journal of Geochemical Exploration,2000,69(70):505-509
[23]Liu L M, Peng S L, Zhang Y H. Numerical geodynamic modeling of coupled mechano-thermo-hydrological processes and its application in predictive exploration in the Fenghuangshan ore field, Tongling[J]. Geotectonica et Metallogenia,2005,29(2):164-173
[24]古成中,吴新跃.有限元网格剖分及发展趋势[J].计算机科学与探索,2008,2(3)：248-259
[25]Cundall P A. A computer model for simulating progressive large scale movement in block rock system[C]. Symposium ISRM,1971, Proc2:129-136
[26]Ransing R S, Gethin D T, Khoei A R, et al. Powder compaction modelling via the discrete and finite element method[J]. Materials and Design,2000,21 (4): 263-269
[27]胥建龙,唐志平.离散元与有限元结合的多尺度方法及其应用[J].计算物理,2003,20(6)：477-482
[28]President's information technology advisory committee. Computational scie nce:ensuring America's competitiveness[EB/OL].2005. http://www.nitrd.gov /pitac/reports/20050609_computational/computational.pdf
[29]赵崇斌,B.E.Hobbs, A.Ord,等.用计算地球科学研究方法探讨地质现象的动力学机制—以断层中等距成矿分布为例[J].中国科学D辑：地球科学,2008,38(5)：646-652
[30]Chongbin Zhao, B.E.Hobbs, H.B.Muhlhaus, et al. Computer simulations of coupled problems in geological and geochemical systems[J]. Computer Methods in Applied Mechanics and Engineering,2002,191(29-30): 3137-3152
[31]刘亮明,疏志明,赵崇斌,等.矽卡岩矿床的汇流扩容空间控矿机制及其对深部找矿的意义：以铜陵-安庆地区为例[J].岩石学报,2008,24(8)：1848-1856
[32]谭凯旋,谢焱石,赵志忠,等.构造-流体-成矿体系的反应-输运-力学耦合模型和动力学模拟[J].地学前缘,2001,8(4)：311-321
[33]邓军,翟裕生,杨立强,等.剪切带构造-流体-成矿系统动力学模拟[J].地学前缘,1999,6(1)：115-127
[34]Itasca Consulting Group, Inc. Fast Lagrangian Analysis of Continua in 3 Dimensions, version 3.0, user's manual[R]. USA:Itasca Consulting Group, Inc FLAC3D,2005
[35]Itasca Consulting Group, Inc. The overview, features, options of FLAC3D [EB/OL].2010. http://www.itascacg.com/flac3d/index.php
[36]ANSYS, inc. The introduction of ANSYS products[EB/OL].2010. http://w ww.ansys.com/products/default.asp
[37]Simulia, inc. The introduction of Abaqus FEA[EB/OL].2010. http://www.s imulia.com/products/abaqus_fea.html
[38]李进文,裴荣富,梅燕雄,等.安徽铜陵狮子山铜(金)矿田成矿流体地球化学研究[J].矿床地质,2006,25(4)：427-437
[39]姜章平,张国培,姚孝德,等.安徽铜陵狮子山矿区冬瓜山铜矿床南段勘探地质报告[R].安徽省地质矿产局321地质队,1995
[40]毛先成,戴塔根.地质体几何造型的计算机实现[J].中南工业大学学报,1997,28(5)：419-422
[41]戴晟.地质体三维可视化研究与系统实现[硕士学位论文][D].华东师范大学,2008
[42]董辉.地质体三维可视化原理与方法研究[硕士学位论文][D].中南大学,2004
[43]易善桢,李琦.3D-GIS数据表示和空间插值方法研究[J].中国图形图像学报,1999,4(8)：698-701
[44]邢延涛.复杂地质岩体三维可视化关键技术研究[硕士学位论文][D].华中科技大学,2006
[45]朱海龙.Gocad项目组简介[J].勘探地球物理进展,2008,31(1)：80-81
[46]董梅,慎乃齐,胡辉,等.基于GOCAD的三维地质模型构建方法[J].桂林工学院学报,2008,28(2)：188-192
[47]Olivier Kaufmann, Thierry Martin.3D geological modeling from boreholes, cross-sections and geological maps, application over former natural gasstorages in coal mines[J]. Computers&Geosciences,2008,34(3):278-290
[48]Andrea Bistacchi, Matteo Massironi, Giorgio V.DalPiaz, et al.3D fold and fault reconstruction with an uncertainty model:An example from an Alpine tunnel case study[J]. Computers&Geosciences,2008,34(4):351-372
[49]Kevin Sprague, Eric de Kemp, WinstonWong, et al. Spatial targeting using queries in a 3-D GIS environment with application to mineral exploration[J]. Computers&Geosciences,2006,32(3):396-418
[50]武晓波,王世新,肖春生.Delaunay三角网的生成算法研究[J].测绘学报,1999,28(1)：28-35
[51]彭仪普.地形三维可视化及其实时绘制技术研究[博士学位论文][D].西南交通大学,2002
[52]刘怀辉.平面区域有限元三角网格剖分算法研究[硕士学位论文][D].山东大学,2007
[53]黄晓东,杜群贵,叶邦彦.三维实体有限元自适应网格规划生成[J].计算机辅助设计与图形学学报,2005,17(7)：1446-1451
[54]金世双.离散数据三角网格剖分研究与实现[J].上海海事大学学报,2005,26(3)：85-87
[55]徐永安,赵静,沈玲玲.约束Delaunay三角网格生成地质构造模型[J].南京师范大学学报,2008,8(4)：155-158
[56]李刚,赵玉新.基于边界特征点提取的约束Delaunay三角剖分算法[J].系统仿真学报,2007,19(16)：3734-3738
[57]刘少兵.带断层地震数据的Delaunay三角剖分算法[J].煤田地质与勘探,2008,36(6)：70-72
[58]陈涛,李光耀.STL格式文件的四边形网格剖分与网格光顺[J].中国机械工程,2009,20(5)：522-528
[59]杨伟军,包忠诩,扶名福,等.映射法在三维六面体有限元网格生成中的应用[J].南昌大学学报(工科版),1999,21(4)：39-43
[60]王东风,翟建军,陈文亮.基于映射法的六面体网格生成算法[J].现代设计与先进制造技术,2009,38(5)：25-27
[61]贾虹,卢炎麟,高发兴.高品质全四边形有限元网格生成的铺砌法[J].浙江工业大学学报,2000,28(4)：254-257
[62]余学进,彭迎风.自动生成有限元四边形网格的铺筑法[J].南昌航空工业学院学报(自然科学版),17(4)：24-26
[63]杜丽惠,江春波.多种介质区域有限元网格自动生成技术[J].水利发电,2000,(2)：55-56
[64]吴宝海,王尚锦.参数曲面网格生成的改进波前法[J].计算机辅助设计与图形学学报,2005,17(8)：1686-1690
[65]Wang Desheng, Wan Shui. Robust Delaunay Tetrahedronal Meshing Coupled with Advancing Front Method[J]. Journal of Southeast University(English Edition),2002,18(2):131-135
[66]杨长伟,孙立镌.结点连元法的一种改进算法[J].哈尔滨理工大学学报,2000,5(5)：97-100
[67]李华,程耿东,顾元宪.一种新的全四边形网格快速生成方法——模板法[J].计算结构力学及其应用,1996,13(1)：25-33
[68]张玉峰,朱以文.有限元网格自动生成的典型方法与研究前瞻[J].武汉大学学报(工学版),2005,38(2)：54-59
[69]徐能雄,武雄,汪小刚,等.基于三维地质建模的复杂构造岩体六面体网格剖分方法[J].岩土工程学报,2006,28(8)：957-961
[70]关振群,宋超,顾元宪,等.有限元网格生成方法研究的新进展[J].计算机辅助设计与图形学学报,2003,15(1)：1-13
[71]吕军,王忠金,王仲仁.有限元六面体网格的典型生成方法及发展趋势[J].哈尔滨工业大学学报,2001,33(4)：485-490
[72]Pinto V, Font X, Salgot M, et al. Using 3D Structures and Their Virtual Representation as a Tool for Restoring Opencast Mines and Quarries[J]. Engineering Geology,2002,63(1-2):121-129
[73]王长海,许国,范孝锋,等,GOCAD与ABAQUS的数据接口实现[J].测绘信息与工程,2008,33(5)：41-42
[74]陈邦国,姜章平,张卫平,等.安徽冬瓜山叠生式层状铜矿热液改造型流体研究[J].江苏地质,2002,26(2)：65-69
[75]黄顺生,徐兆文,倪培,等.安徽铜陵冬瓜山热液叠加改造型铜矿床流体包裹体地球化学特征[J].地质找矿论丛,2003,18(1)：34-38