收藏本站
《中国地质大学》 2012年
收藏 | 手机打开
二维码
手机客户端打开本文

不整合面型铀矿流体系统数值模拟

崔涛  
【摘要】:随着全球气候变化的加剧,核能受到越来越多的重视。当前,加拿大和澳大利亚的古元古代的沉积盆地中发育的不整合面型铀矿占据全球铀市场份额的30%多,但是他们的成矿机制却依然不是十分清楚。盆地中成矿流体的运移模式是什么?铀源是什么?流体运移的驱动力是什么?这些问题都需要进一步的研究。数值模拟可以整合已有的大量的地质、地球物理和地球化学数据,从而提供新的认识。在本论文中,作者利用有限元和有限差分技术,设计了大量的数值模拟实验来研究与不整合面型铀矿有关的成矿流体系统。首先,基于来自Athabasca盆地,Thelon盆地和Kombolgie盆地的地质、地球物理及地球化学数据,建立了一系列的概念模型。以这些概念模型为基础,进行了大量的数值实验来揭示流体流动、热传输、地形起伏和构造变形等过程的相互作用及其对成矿的影响。FEFLOW和FLAC被用来求解与这些过程相关的控制方程,比如质量守恒方程,动量守恒方程,和能量守恒方程等。根据本研究的需要,作者也对FEFLOW和FLAC的功能做了相应的扩展和改进。论文取得的主要研究结果如下: 1.热液对流在砂岩层中广泛发育。在典型的水文地质条件下,模拟了地热梯度介于25到35度/千米之间的流体流动,发现热液对流在有效孔隙度相对较高的砂岩层中广泛出现。热液对流对温度场扰动明显,最高可达40度,在上升流附近,温度升高,而在下降流附近,温度降低。具体的流体形态和流体速度受水文地质属性,地热梯度,地层几何特征及砂岩层厚度等的影响。较高的地热梯度会加速流体运移,但是盆地规模的基本对流模式不变。在盆地最边缘砂岩层低于1千米的地方,没有热液对流发育。空隙介质厚度是雷诺数重要的控制参数,当厚度变小时,雷诺数会降低,当其降低超过临界值时,对流便不会产生。在同一个对流环内部,流体的流动速度由外向内逐渐降低。高透水砂岩和盖层及基底之间的转换带为流体快速运移区,从而有利于矿体发育,这种特点与大部分矿床发育于基底不整合面是一致的。此外,有些研究者认为砂岩层被相对低透水单元的分隔会限制自由对流的产生,而根据本论文研究,当高透水的砂岩层被较薄的低透水单元隔离时,虽然流体运移速度会降低,但是自由对流依然发育。 2.地形起伏对盆地流体系统的运移模式影响明显,但它不是不整合面型铀矿流体系统的主要驱动力。当潜水面坡度较小时,流体流动被自由热对流所控制,但是随着坡度的加大,流体流动逐渐会被地形起伏控制。在坡度低于0.001时,整个流体系统会显示出混合对流的特点,在盆地中部,透水性砂岩厚的地方依然被自由热对流所控制,而盆地边缘却以地形驱动流为主。当地形起伏变为主要的流体驱动机制时,流体的属性(温度,压强,组分)通常会在流体流动方向上展现渐变和不对称的特征。比如,补给带的温度会降低,而泄水区的温度会升高。如果地形起伏是第一驱动力,那么矿床通常只在盆地的一侧发育。美国科罗拉多高原的板状铀矿(Tabular uranium deposits)便是在这种成矿流体作用下形成的。但是,对于不整合面型铀矿而言,矿体在盆地边缘成辐射状分布,没有明显的非对称特征。因此,地形起伏应该不是不整合面型铀矿系统中主要的流体驱动机制。 3.在不考虑机械变形对流体压力影响的前提下,断层对流体运移的影响主要取决于不同断层之间的连通性和断层是否出露地表。当三条相互孤立的宽400米的盲断层被加入到盆地尺度的概念模型中时,盆地规模的流体流动没有明显的变化。但是,断层的加入会打破已有的压力平衡,从而提高流体运移速度,但是经过10万年左右的时间,这种扰动会逐渐消失。当这些断层贯穿盖层并且被高透水性的砂岩连通时,对流模式将发生明显变化,断层一般会成为不同对流环的边界,从而成为优先的上升和下降流通道。地表流体能够沿着断层进入盆地流体系统,盆地流体也能够沿着断层溢出至盆地浅部。 4.基底不整合面(带)为氧化还原界面和典型的高速流体运移区。基底不整合面(带)的厚度从几厘米到几百米不等。由于盆地规模模型网格精度的限制,比较薄的基底不整合面(带)无法模拟。当一个厚250米的基底不整合面(带)被加入概念模型时,对200米/年,400米/年和1000米/年的水平导水率分别进行了敏感度测试。当独立的高透水率的基底不整合面(带)被加入概念模型时,盆地规模的自由热对流仍然发育,但是新加入的不整合面(带)为流体流动速度最高和流体通量最大区域。 5.构造运动中发生的机械变形会导致流体压力场变化,进而引发混合对流;在应力变化强烈的区域,变形驱动流为主,而在应力变化比较弱的区域,自由热对流依旧是主要的流体运移模式,具体的流体形态取决于应变速率及其各水文地质单元的属性。在相同的应力场下,不同地层表现出不同的压力场的演化,进而引发存在于不同层组之间的压力梯度。当断层或者其它的构造连接这些这些不同的水文地质单元时,跨越不同单元的流体流动将产生。在伸展变形过程中,整个系统的流体压力会降低,但是基底中压力降低的速度要比其它水文地质单元要快,所以盆地中的氧化性流体会被“吮吸”到基底中,从而有助于形成在基底中发育的矿体。挤压变形中,基底中的压力上升速度比上覆地层要快,基底中的还原性流体将被挤压到上面的砂岩层中,这种流体运移模式有助于在砂岩层中发育矿体。这种构造背景和矿床发育位置之间的对应关系被大部分已经发现的矿体所支持。澳大利亚的Kombolgie盆地形成后,经历了一些列的伸展变形构造,而在此盆地中的大部分矿体是发育在基底中的。加拿大Althabasca盆地的情况相对复杂,66个相对研究程度比较高的矿床中,32个是在基底中发现的,其它的34个跨越基底或者在离基底不整合面附近的砂岩里面,这也许是这些盆地经历了相对更复杂的构造过程的结果。 6.当盐度变化和溶质迁移被耦合进系统后,盆地规模流体运移模式和热液自由对流相似,对流环广泛发育于高透水砂岩层中。但是,借助FEFLOW6里面新的数据可视化技术,发现了跨越基底不整合面的流体运移。砂岩和基底中典型的流体运移速度分别为3.5米/年和2.5×10-6米/年,虽然流体流动速度相差若干数量级,但是这种现象为从基底中析出铀提供了一种可能的机制。通过溶质运移模拟发现,伴随着热盐对流,铀元素能够从基底通过平流(advection),离子扩散(ion diffusion)和机械散布(mechanical dispersion)慢慢迁移到盆地。主要的迁移方式为伴随流体流动的平流。 7.当没有强烈的构造运动时,在典型的水文地质条件下,热盐对流可以穿透到基底不整合面1-2千米以下。假设在基底中部有一个铀源,铀元素能够在热盐对流的作用下逐步迁移到沉积盆地中去。反之,如果初始铀源区被放在高渗透性的砂岩中,经过500万年铀元素也可以慢慢迁移到基底不整合面以下2千米以内。铀源区的相对位置对铀传输的效率也有很大影响。在盆地边缘,基底不整合面坡度较大的地方,铀迁移的效率和速度要比铀源区在盆地中间的情况下高,因为不整合面坡度的增加会提高流体流动速度。基于这些结果,在含有不整合面型铀矿的这些古元古代的盆地里面,热盐对流可以引起跨越基底不整合面的广泛的流-流交互和流-岩交互,并且热盐对流可以从基底中提取足够的铀,从而为形成大规模的不整合面型铀矿提供铀源。 8.从水动力学角度来看,基底可以作为铀源区。对于不整合面型铀矿铀源的争议,已经持续了几十年。根据当前广泛认可的成岩热液模型,矿床形成需要氧化性的盆地流体与还原性的基底流体或矿物在不整合面附近进行交互。铀元素必须由氧化性的盆地流体从盆地充填物或者基底中的副矿物(accessory minerals)或者以前形成的铀矿物里萃取。如上所述,高透水的砂岩中自由对流发育,如果有足够的副矿物,例如锆石和独居石,那么在合适的物理化学条件下,从盆地中萃取大量的铀元素是有可能的。同时,这些古元古代砂岩中缺少有机质和Fe2+,所以流体的氧化性可以很好的维持。但是,好多研究者认为从锆石中提取铀矿物的条件并不满足,而独居石含量又不够,所以铀源应该在其它区域。另一方面,基底中的太古代的花岗岩,淡色花岗岩,伟晶岩等铀矿物含量很高,并且易于萃取。但是,以前的认识认为,由于基底极低透水性的影响,从基底中萃取大量铀元素是不可能的。而我们的模拟结果说明,从水动力学角度来看,盆地流体可以渗入到盆地基底中并且返回盆地,所以基底是有可能作为铀源区的。
【学位授予单位】:中国地质大学
【学位级别】:博士
【学位授予年份】:2012
【分类号】:P619.14

手机知网App
【共引文献】
中国期刊全文数据库 前1条
1 凌洪飞;;论花岗岩型铀矿床热液来源——来自氧逸度条件的制约[J];地质论评;2011年02期
中国博士学位论文全文数据库 前1条
1 吴鹏;楚雄盆地砂岩型铜矿床地质地球化学及找矿预测研究[D];昆明理工大学;2009年
【相似文献】
中国期刊全文数据库 前10条
1 孙泽轩;吴英;解波;巫声扬;康友纯;陈进文;;江川-建水地区中上元古界铀矿化特征及其控矿因素[J];铀矿地质;2009年05期
2 闽山;;乌克兰铀矿勘探工作若干见闻[J];铀矿地质;1993年01期
3 丁万烈;沃罗涅日结晶地块和下莫斯科褐煤盆地的含铀远景[J];世界核地质科学;2003年04期
4 张书成,白云生,吴慧山,伍岳,常桂兰;康滇地轴南段不整合面型铀成矿基本条件的物化探研究[J];中国核科技报告;2000年00期
5 ;《铀矿地质》1990年第6卷总目录[J];铀矿地质;1990年06期
6 李耀菘;中国铀矿床同位素地球化学特征(英文)[J];中国核科技报告;1995年00期
7 王文广,陶铨,张守本;辽东中部地区区域铀矿地质特征及其成矿远景评价(英文)[J];中国核科技报告;1998年00期
8 程小久,匡耀求;盆地成矿流体系统研究中的几个重要问题[J];华南理工大学学报(自然科学版);1996年S1期
9 蒙义峰,杨竹森,曾普胜,徐文艺,王训成;铜陵矿集区成矿流体系统时限的初步厘定[J];矿床地质;2004年03期
10 陈祖伊,李子颖,周维勋,管太阳;产铀盆地的形成演化模式及其鉴别标志(英文)[J];中国核科技报告;2004年01期
中国重要会议论文全文数据库 前10条
1 ;Biogenic and petroleum-related ore-forming processes in Dongsheng uranium deposit,NW China[A];中国科学院地质与地球物理研究所2007学术论文汇编(第三卷)[C];2008年
2 ;Mineralogical and geochemical evidence for coupled bacterial uranium mineralization and hydrocarbon oxidation in the Shashagetai deposit,NW China[A];中国科学院地质与地球物理研究所2007学术论文汇编(第三卷)[C];2008年
3 ;Biogenic spherical magnetite and coffinite from a sandstone-type uranium deposit in the Dongsheng area[A];中国科学院地质与地球物理研究所2007学术论文汇编(第三卷)[C];2008年
4 ;Mineralogical and geochemical evidence for coupled bacterial uranium mineralization and hydrocarbon oxidation in the Shashagetai deposit,NW China[A];中国科学院地质与地球物理研究所2006年论文摘要集[C];2007年
5 张坚;荆高山;许俊泉;周玉祥;程京;;利用微泵和微阀构建的用于细胞分离的流体系统[A];中国生物工程学会第三次全国会员代表大会暨学术讨论会论文摘要集[C];2001年
6 张志成;;实现流体系统数字化的新思想——脉冲数字流法[A];中国力学学会学术大会'2005论文摘要集(上)[C];2005年
7 侯增谦;李荫清;张琦玲;曲晓明;;火山成因块状硫化物矿床的流体系统与混合过程:来自白银厂矿床和呷村矿床的流体包裹体证据[A];全国包裹体及地质流体学术研讨会论文摘要[C];2002年
8 王强;薛鑫;王长福;秦浩;;液压系统状态监测的磨粒分析方法[A];中国公路学会筑路机械学会第十四次学术年会论文集[C];2003年
9 ;“我们的目标是提供高质量流体系统配件和解决方法”世伟洛克公司[A];中国仪器仪表学会2005年学术年会测控技术与节能环保学术会议论文集[C];2005年
10 侯增谦;张绮玲;曲晓明;李振清;唐绍华;余金杰;;现代与古代海底热流系统流体地质-地球化学与动力学过程[A];“十五”重要地质科技成果暨重大找矿成果交流会材料二——“十五”地质行业获奖成果资料汇编[C];2006年
中国重要报纸全文数据库 前5条
1 刘洋;小部件有大学问一“阀”防溢保卫生[N];中国信息报;2009年
2 王秀兰;生物质合成二甲醚首获成功[N];中国化工报;2003年
3 ;世伟洛克变截面流量计提供定制校准[N];中国高新技术产业导报;2009年
4 班玮;德国开发出自动感应水龙头节电技术[N];中国高新技术产业导报;2010年
5 本报记者 赵宁;美晨科技:国内商用车非轮胎橡胶制品行业龙头[N];证券日报;2011年
中国博士学位论文全文数据库 前8条
1 崔涛;不整合面型铀矿流体系统数值模拟[D];中国地质大学;2012年
2 赵双良;对多孔介质中流体的研究[D];湖南大学;2007年
3 田建祥;类氩流体的物态方程和气液一级相变[D];大连理工大学;2006年
4 王勇;云南巍山—永平铜金多金属矿化集中区成矿流体特征及流体地质填图研究[D];中国地质大学(北京);2003年
5 赵克中;磁耦合双向驱动控制器的研究[D];东北大学;2006年
6 曹金涛;磁化等离子体中Richtmyer-Meshkov不稳定性的理论研究[D];中国科学技术大学;2009年
7 杜建国;大别造山带中生代岩浆作用与成矿地球化学研究[D];合肥工业大学;2000年
8 严森;厌氧条件下纳米铁还原水中六价铀的反应动力学和机理研究[D];中国地质大学;2010年
中国硕士学位论文全文数据库 前10条
1 赵士洋;流体系统元件损失建模及仿真软件开发[D];南京航空航天大学;2011年
2 杨艳;河南围山城金银多金属成矿带成矿流体研究[D];中国地质大学(北京);2008年
3 邵亚雄;三偏心金属密封蝶阀的力学特性分析及结构优化研究[D];西华大学;2007年
4 孙艳;纳沟道中扩散系数测量及流动特性研究[D];浙江大学;2007年
5 刘遥;开放系统化学势与压强对应关系的研究[D];安徽大学;2007年
6 杨章锋;内流对海底悬跨管线管壁应力影响的数值研究[D];大连理工大学;2008年
7 王霆;核地球物理方法在铀矿勘查中的应用[D];中国地质大学(北京);2011年
8 张容芳;铀胁迫对两种蓝藻生长及抗氧化酶活性的影响[D];南华大学;2012年
9 蒋行国;流体设备运行参数的网络化虚拟测试系统研究[D];重庆大学;2003年
10 刘昌军;非饱和水气两相渗流数值模拟研究及应用[D];河海大学;2005年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62791813
  • 010-62985026