基于FLAC3D成矿过程数值模拟:以南泥湖钼矿床为例
|
摘要
南泥湖钼矿床是一个斑岩-矽卡岩二位一体的矿床,成矿的三大核心条件是控岩构造条件、岩浆条件和物源条件。本文首先通过GOCAD对南泥湖钼矿床建立理想的成矿模式,并运用C#编写的程序把模型的格式转换成FLAC3D可以识别的格式,导入FLAC3D中。然后通过FLAC3D对南泥湖矿床岩体侵入后固化冷却过程中进行力-热-流耦合数值模拟,通过数值模拟结果显示流体流向岩体接触带附近,并在南泥湖组出现最大体应变带形成扩容空间,扩容空间与流体交汇处形成汇流区是成矿的有利部位。
The Nannihu Mo deposit is a porphyry-skarn the two incorporated into one of the deposits.Metallogenic conditions of three core is accused rock structure conditions,magmatic conditions and sediment source.In this paper,first of all,established the ideal metallogenic model of Nannihu Mo deposit by GOCAD software,then converted into a format that can be imported into the FLAC3 D,by means of C+ +programming.the coupled mechanical—thermal-hydrological numerical modeling of hydrothermal metallogenic processes for the Nannihu Mo deposit was performed in FLAC3 D.The results show that,fluid flow focusing near the contact rock zone of the Nannihu units.These areas with the maximum dilation.Maximum dilation and fluid interchange form confluence area,are favorable for mineralization.
The Nannihu Mo deposit is a porphyry-skarn the two incorporated into one of the deposits.Metallogenic conditions of three core is accused rock structure conditions,magmatic conditions and sediment source.In this paper,first of all,established the ideal metallogenic model of Nannihu Mo deposit by GOCAD software,then converted into a format that can be imported into the FLAC3 D,by means of C+ +programming.the coupled mechanical—thermal-hydrological numerical modeling of hydrothermal metallogenic processes for the Nannihu Mo deposit was performed in FLAC3 D.The results show that,fluid flow focusing near the contact rock zone of the Nannihu units.These areas with the maximum dilation.Maximum dilation and fluid interchange form confluence area,are favorable for mineralization.
引文
[1]Liu L M,Zhang Y H.Numerical modeling of the coupled mechanical and hydrological processes during deformation and mineralization in the Mount Isa Block,Australia[R].2007.
[2]Lin G,Zhang Y,Guo F,et al.Numerical modelling of lithosphere evolution in the North China Block:Thermal versus tectonic thinning[R].2005.
[3]Mallet J.L.Discrete modeling for nature objects[J].Mathematical Geology,1997,29(2):199-218.
[4]Norton D L,Dutrow B L.Complex behavior of magjmahydrothermal Proeess:role of supercritical fluid[J].Geochimiea et CosmoehimieaAeta,2001,65:4009-4017.
[5]England P C,Mckenzie D A.Thin viscous sheetmodel for conti2 nental deformation[J].Journal of Geophysics Research,1982,70:2952321.
[6]龚纪文,崔建军,席先武,等.FLAC数值模拟软件及其在地学中的应用[J].大地构造与成矿学,2002,26(3):321-325.
[7]陈庆宣,王维襄,孙叶,等.岩石力学与构造应力场[M].北京:地质出版社,1998.
[8]陈繁荣,裘愉卓.江西贵溪冷水坑多金属矿床成矿过程流体地球化学模拟及其地质意义[J].地球化学,1995,24(zk):24-32.
[9]邓军,翟裕生,杨立强,等.剪切带构造流体成矿系统动力学模拟[J].地学前缘,1999,6(l):115-128.
[10]邓军,杨立强,翟裕生,等.构造流体成矿系统及其动力学的理论格架与方法体系[J].地球科学,2000,25(1):71-78.
[11]粱正召,丙勇勤,唐春安,等.露天矿潜滑体边坡变形、失稳模拟分析[J].地质灾害与环境保护,2000,11(3):235-237.
[12]孙叶,谭成轩.构造应力场研究与实践[J].地质力学学报,2001,7(3):54-58.
[13]任启江,郭国章,冯祖钧,等.陕西金堆城斑岩钼矿成矿过程中热及流体传输的计算模拟[J].矿床地质,1994,13(1):88-94.
[14]刘亮明,疏志明,赵崇斌,等.矽卡岩矿床的汇流扩容空间控矿机制及其对深部找矿的意义:以铜陵一安庆地区为例[J].岩石学报,2008(8):1545-1556.
[15]刘亮明,周瑞超,赵崇斌.构造应力环境对浅成岩体成矿系统的制约:从安庆月山岩体冷却过程动力学计算模拟结果来分析[J].岩石学报,2010,26(9):2869-2878.
[16]龚纪文,崔建军,席先武,等.FLAC数值模拟软件及其在地学中的应用[J].大地构造与成矿学,2002,26(3):321-325.
[17]赵义来,刘亮明.复杂形态岩体接触带成矿祸合动力学三维数值模拟:以安庆铜矿为例[J].大地构造与成矿学,2011,125(35):128-136.
[18]孙涛.铜陵冬瓜山铜矿床的三维形态及成矿动力学计算模拟[D].长沙:中南大学,2010.
[19]孙涛,刘亮明,赵义来,等.基于GOCAD平台的复杂地质体系的动力学建模研究[J].矿产与地质,2011,25(2):163-167.
[20]杨立强,张中杰,林炯,等.FIAC基本原理及其在地学中的应用[C]//中国科学院地质与地球物理研究所二○○三学术论文汇编·第二卷(青藏高原),2003.
[21]于淼.热液成矿过程数值模拟研究—以个旧锡矿为例[D].北京:中国地质大学(北京),2014.
[22]岑况,於崇文.铜陵地区硫化物矿床成矿过程的热传导和物质输运动力学[J]地球化学,2001,30(6):533-539.
[23]刘孝善,吴澄宇,严正富,等.河南栾川南泥湖-三道庄钼(钨)矿床围岩蚀变的地质地球化学研究[J].地质找矿论丛,1986,1(4):27-39.
[24]汪江河.河南三道庄钨钼矿区深部下层矿的找矿意义[J].地质调查与研究,2009,32(4):291-296.
[25]毛景文,李晓峰,张荣华,等.深部流体成矿系统[M].北京:中国大地出版社,2005.
[26]翁纪昌,付治国,黄超勇,等.南泥湖钼矿田成矿物质组成研究[J].地质调查与研究,2008(2):97-106.
[27]刘永春,付治国,高飞,等.河南栾川南泥湖特大型钼矿床成矿母岩地质特征研究[J].中国钼业,2006(3):13-17.
[2]Lin G,Zhang Y,Guo F,et al.Numerical modelling of lithosphere evolution in the North China Block:Thermal versus tectonic thinning[R].2005.
[3]Mallet J.L.Discrete modeling for nature objects[J].Mathematical Geology,1997,29(2):199-218.
[4]Norton D L,Dutrow B L.Complex behavior of magjmahydrothermal Proeess:role of supercritical fluid[J].Geochimiea et CosmoehimieaAeta,2001,65:4009-4017.
[5]England P C,Mckenzie D A.Thin viscous sheetmodel for conti2 nental deformation[J].Journal of Geophysics Research,1982,70:2952321.
[6]龚纪文,崔建军,席先武,等.FLAC数值模拟软件及其在地学中的应用[J].大地构造与成矿学,2002,26(3):321-325.
[7]陈庆宣,王维襄,孙叶,等.岩石力学与构造应力场[M].北京:地质出版社,1998.
[8]陈繁荣,裘愉卓.江西贵溪冷水坑多金属矿床成矿过程流体地球化学模拟及其地质意义[J].地球化学,1995,24(zk):24-32.
[9]邓军,翟裕生,杨立强,等.剪切带构造流体成矿系统动力学模拟[J].地学前缘,1999,6(l):115-128.
[10]邓军,杨立强,翟裕生,等.构造流体成矿系统及其动力学的理论格架与方法体系[J].地球科学,2000,25(1):71-78.
[11]粱正召,丙勇勤,唐春安,等.露天矿潜滑体边坡变形、失稳模拟分析[J].地质灾害与环境保护,2000,11(3):235-237.
[12]孙叶,谭成轩.构造应力场研究与实践[J].地质力学学报,2001,7(3):54-58.
[13]任启江,郭国章,冯祖钧,等.陕西金堆城斑岩钼矿成矿过程中热及流体传输的计算模拟[J].矿床地质,1994,13(1):88-94.
[14]刘亮明,疏志明,赵崇斌,等.矽卡岩矿床的汇流扩容空间控矿机制及其对深部找矿的意义:以铜陵一安庆地区为例[J].岩石学报,2008(8):1545-1556.
[15]刘亮明,周瑞超,赵崇斌.构造应力环境对浅成岩体成矿系统的制约:从安庆月山岩体冷却过程动力学计算模拟结果来分析[J].岩石学报,2010,26(9):2869-2878.
[16]龚纪文,崔建军,席先武,等.FLAC数值模拟软件及其在地学中的应用[J].大地构造与成矿学,2002,26(3):321-325.
[17]赵义来,刘亮明.复杂形态岩体接触带成矿祸合动力学三维数值模拟:以安庆铜矿为例[J].大地构造与成矿学,2011,125(35):128-136.
[18]孙涛.铜陵冬瓜山铜矿床的三维形态及成矿动力学计算模拟[D].长沙:中南大学,2010.
[19]孙涛,刘亮明,赵义来,等.基于GOCAD平台的复杂地质体系的动力学建模研究[J].矿产与地质,2011,25(2):163-167.
[20]杨立强,张中杰,林炯,等.FIAC基本原理及其在地学中的应用[C]//中国科学院地质与地球物理研究所二○○三学术论文汇编·第二卷(青藏高原),2003.
[21]于淼.热液成矿过程数值模拟研究—以个旧锡矿为例[D].北京:中国地质大学(北京),2014.
[22]岑况,於崇文.铜陵地区硫化物矿床成矿过程的热传导和物质输运动力学[J]地球化学,2001,30(6):533-539.
[23]刘孝善,吴澄宇,严正富,等.河南栾川南泥湖-三道庄钼(钨)矿床围岩蚀变的地质地球化学研究[J].地质找矿论丛,1986,1(4):27-39.
[24]汪江河.河南三道庄钨钼矿区深部下层矿的找矿意义[J].地质调查与研究,2009,32(4):291-296.
[25]毛景文,李晓峰,张荣华,等.深部流体成矿系统[M].北京:中国大地出版社,2005.
[26]翁纪昌,付治国,黄超勇,等.南泥湖钼矿田成矿物质组成研究[J].地质调查与研究,2008(2):97-106.
[27]刘永春,付治国,高飞,等.河南栾川南泥湖特大型钼矿床成矿母岩地质特征研究[J].中国钼业,2006(3):13-17.