地球化学反应模型用于水—岩相互作用的研究
|
摘要
水-岩相互作用是地球科学领域十分关注的问题。多年来科学家们通过实验室内反复实验和野外观察,积累了大量的水-岩体系热力学和动力学参数,这为深入了解自然界水-岩作用现象和规律提供丰富的理论基础和实际数据资料。但这些工作一般都是通过化学反应测试获得,有很大的局限性。随着计算机科学的发展,通过建立适当的地球化学模型,用以模拟计算自然界中的水-岩相互作用,近年来国际上相关领域发展迅速。
本文在地球化学反应平衡模拟理论基础上,建立了矿物溶解滴定模型,选用美国地质调查局(USCS)提供的Phreeqc模拟软件,研究了不同条件下,pH、温度、CO_2分压、离子强度等条件矿物溶解作用的影响。同时,在前人实地测试勘查的数据资料基础上,选择珠海校区的地下水-岩系统作为研究对象,通过软件计算,考察多种影响因素,进行了研究区的的水-岩地球化学反应模拟。
模拟结果表明:(1)珠海校区地下水中典型矿物,如文石、方解石、白云石等都处于非饱和状态,由于这些矿物是灰岩的重要组成,其饱和指数均小于零,所以灰岩可能会持续溶解。(2)pH值、温度、二氧化碳分压、离子强度、同离子效应等因素的变化对于体系的平衡具有重要影响,文中对这些应影响作用进行了定量研究,发现温度、二氧化碳分压的影响满足一定的函数关系式。(3)研究了酸雨对水-岩作用的影响。模拟了Ca~(2+)、NH_4~+等离子对于典型矿物如文石、方解石、白云石等溶解反应平衡的影响。(4)模拟数据显示其地下水pH小于7,根据前人的实地调查数据资料,如果不加以控制处理,将会使其pH值进一步降低,将会给环境带来的潜在风险。
Water-rock interaction is an attractive subject in earth sciences. Scientists have got a lot of thermodynamics and kinetics data of water-rock systems by means of experiments and investigations, which provide the basic theory and information to well understand water-rock interaction in the nature. However, these results are often obtained through geochemical analysis during reactions and this approach has many limitations. With the development of computer science, the study on rock-water interaction by geochemical modeling has been developed in recent years.
Based on the theory of geochemical reaction equilibrium modeling, a kind of titration model was established in the thesis. The simulation software Phreeqc donated by USGS was used to study the mineral dissolution under different conditions such as pH, temperature, CO_2 partial pressure, ionic strength, etc. On the basis of data materials surveyed by others, groundwater-rock system in Zhuhai school zone was studied by geochemical reaction simulation.
The results showed as follows. (1) Some minerals, such as aragonite, calcite, dolomite, were not saturated with groundwater. It indicated that limestone might be continued to dissolve. (2) The equilibrium of these systems was influenced by the factors including pH, temperature, CO_2 partial pressure, ionic strength and common-ion effect. It also indicated that temperature and CO_2 partial pressure have a relationship of mathematic functions with species. (3) The effect of acid rain was studied on water-rock interaction. It also showed ions such as Ca~(2+)、NH_4~+ would impact aragonite, calcite and dolomite dissolution equilibrium. (4) The simulation showed pH of groundwater would be less than 7, which might be a potential risk in the environment.
本文在地球化学反应平衡模拟理论基础上,建立了矿物溶解滴定模型,选用美国地质调查局(USCS)提供的Phreeqc模拟软件,研究了不同条件下,pH、温度、CO_2分压、离子强度等条件矿物溶解作用的影响。同时,在前人实地测试勘查的数据资料基础上,选择珠海校区的地下水-岩系统作为研究对象,通过软件计算,考察多种影响因素,进行了研究区的的水-岩地球化学反应模拟。
模拟结果表明:(1)珠海校区地下水中典型矿物,如文石、方解石、白云石等都处于非饱和状态,由于这些矿物是灰岩的重要组成,其饱和指数均小于零,所以灰岩可能会持续溶解。(2)pH值、温度、二氧化碳分压、离子强度、同离子效应等因素的变化对于体系的平衡具有重要影响,文中对这些应影响作用进行了定量研究,发现温度、二氧化碳分压的影响满足一定的函数关系式。(3)研究了酸雨对水-岩作用的影响。模拟了Ca~(2+)、NH_4~+等离子对于典型矿物如文石、方解石、白云石等溶解反应平衡的影响。(4)模拟数据显示其地下水pH小于7,根据前人的实地调查数据资料,如果不加以控制处理,将会使其pH值进一步降低,将会给环境带来的潜在风险。
Water-rock interaction is an attractive subject in earth sciences. Scientists have got a lot of thermodynamics and kinetics data of water-rock systems by means of experiments and investigations, which provide the basic theory and information to well understand water-rock interaction in the nature. However, these results are often obtained through geochemical analysis during reactions and this approach has many limitations. With the development of computer science, the study on rock-water interaction by geochemical modeling has been developed in recent years.
Based on the theory of geochemical reaction equilibrium modeling, a kind of titration model was established in the thesis. The simulation software Phreeqc donated by USGS was used to study the mineral dissolution under different conditions such as pH, temperature, CO_2 partial pressure, ionic strength, etc. On the basis of data materials surveyed by others, groundwater-rock system in Zhuhai school zone was studied by geochemical reaction simulation.
The results showed as follows. (1) Some minerals, such as aragonite, calcite, dolomite, were not saturated with groundwater. It indicated that limestone might be continued to dissolve. (2) The equilibrium of these systems was influenced by the factors including pH, temperature, CO_2 partial pressure, ionic strength and common-ion effect. It also indicated that temperature and CO_2 partial pressure have a relationship of mathematic functions with species. (3) The effect of acid rain was studied on water-rock interaction. It also showed ions such as Ca~(2+)、NH_4~+ would impact aragonite, calcite and dolomite dissolution equilibrium. (4) The simulation showed pH of groundwater would be less than 7, which might be a potential risk in the environment.
引文
Anon.Visual MODFLOW V.2.8.2 User's Manual for Professional Applications in Three-Dimensional Groundwater Flow and Contaminant Transport Modeling [Z]. Ontario:Waterloo Hydrogeologic Inc.,2000. 1~3.
Garrels. R.M.&M.E.Thompson. Geochemical model for sea water at 25℃and one atmosphere total pressure atmosphere. Jour. Sci.,1962. 58~66
Hegeson.H.C. Evaluation of irreversible reations in geoehemical processes involving minerals and aqueous solutions. Thermodynamics relationship.Geochim.Cosmochim,1968. 853~877
陈家玮,鲍征宇.地球化学反应平衡模型的建模方法.地质科技情报,2001,20(3):41~45
邓海琳,涂光炽,李朝阳.水/岩比的地球化学意义.矿物学报,1999,19(3):267~272
郭乙木.《线性与非线性有限元及其应用》.北京:机械工业出版社,2005.156-158
蒋引珊,王玉洁,徐长耀等.热分析法研究矿物分解过程动力学.长春科技大学学报,2000, 30(1):91~93
蒋倩,李心清,丁文慈等.岩溶水对不同环境的响应以及其与大气CO2源汇的关系——基于热力学的研究.矿物岩石地球化学通报,2000,25(3):227~233
李义连,王焰新,周来茹,等.地下水矿物饱和度的水文地球化学模拟分析——以娘子关泉域岩溶水为例.地质科技情报,2003,21(1) :33~35
刘再华.不同CO2分压条件下的白云岩溶解动力学机理.中国科学,2001,31(4):378~384
刘再华.方解石溶解、沉积速率控制的物理、化学机制.重要会议:中国地质科学院“九五”科技成果汇编,2001,253~256
刘再华,W. Dreybrodt,李华举.灰岩和白云岩溶解速率控制机理的比较.地球科学—中国地质大学学报,2006,1(3):412~415
刘再华,W.Dreybrodt,韩军,等. CaCO3-CO2-H2O岩溶系统的平衡化学及分其析.中国岩溶,2005,24(1):1~14
李进文,裴荣富,张荣华,等.安徽铜陵胡村矽卡岩型铜矿水-岩反应动力学实验研究.高校地质学报,2007,13(2):261~271
李义连,王焰新,张江华,等.娘子关泉域岩溶水硫酸盐污染的地球化学模拟分析.地球科学——中国地质大学学报,2000,25(5):467~471
刘兴起,蔡克勤,于升松.基于Pitzer模型的茶卡盐湖卤水和盐类矿物的形成进行了地球化学模拟.中国科学.2003,33(3):247~254
马瑞.水岩相互作用研究-以山西太原为例:[博士学位论文].武汉:中国地质大学,2007
马尧.地球化学模式PHREEQC在地侵工艺中的应用.铀矿冶.2007,26(2):67~71
欧亚波.加拿大Thompson镍矿矿山废水/尾砂作用:[硕士学位论文].四川:成都理工大学,2008
潘根兴,滕永忠.土壤一灰岩岩溶系统中水文地球化学动力学过程模拟及其意义.地球化学,2000,29(3):274~276
乔淑卿,杨作升,潘燕俊,等.沉积物中碳酸盐(方解石和白云石)测定方法的研究.中国海洋大学学报,2006,36(3):484~488
齐雯,郑绵平.西藏扎布耶盐湖卤水演化的Pitzer模型模拟.地质学报,2007,81(12):1735~1741
沈照理,王焰新.水—岩相互作用的回顾与展望.地球科学——中国地质大学学报,2002,27(2):128~133
施枫芝,迟宝明,潘杰.基于化学动力学理论的中深层地下水资源评价.水土保持研究,2007,14(5):211~238
孙连环.塔里木盆地塔中碳酸盐岩储层酸岩反应动力学实验研究.石油与天然气化工.2006,25(1):51~54
陶小晚.不同二氧化碳分压环境近平衡状态下海洋碳酸盐溶解及沉淀动力学研究:[硕士学位论文].北京:中国科学院海洋研究所,2007
童海涛,宋汉周.水-岩作用系统的水文地球化学模拟.水科学进展,2004,15(2):211~215
王丽,王金生,林学钰.水文地球化学模型研究进展.水文地质工程地质,2003(6):105~109
王永新.水-岩相互作用机理及其对库岸边坡稳定性影响的研究:[硕士学位论文].重庆:重庆大学,2006
王益峰,冯祖钧.一个计算机模拟水-岩反应过程的数学模型.矿床地质,1987,6(2):77~88
肖林萍,黄思静.方解石和白云石溶蚀实验热力学模型及地质意义.矿物岩石,2003,23(1):114~116
肖林萍,黄思静.碳酸盐岩溶蚀实验热力学模型及工程地质意义.西南交通大学学报,200237(3):250~253
谢水波,陈泽昂,张晓健,等.铀尾矿库区浅层地下水中U(Ⅵ)迁移的模拟.原子能科学技术,2007,41(1):58~64
徐乐昌.地下水模拟常用软件介绍.铀冶矿,2002,21(1):34~28
于炳松,赖兴运.成岩作用中的地下水碳酸体系与方解石溶解度.沉积学报,2006,24(5):628~635
袁道先,章程.岩溶动力学的理论探索与实践.地球学报,2008,29(3):355~365
杨涛,周文斌.江西科学,2005,23(1):28~33
张国军.夹层软泥化过程及其化学动力学模拟:硕士论文.江苏,河海大学.2007
张文涛,关平等.成岩作用的热力学分析及应用实例.自然科学进展,2007,17(9):1209~1215
张建利,孙忠杰,张泽兰.碳酸盐岩油藏酸盐反应动力学实验研究.油田化学,2003,20(3):216~218
B.Jmerkel B.Planer-Friedrich.水文地球化学模拟程序.(朱义年,王焰新译).武汉:中国地质大学出版社,2005.75~94
周海燕,周训,姚锦梅.广东从化温泉的水文地球化学模拟.现代地质, 2007,21(4):620~623
张凤娥,卢耀如,郭秀红等.符合岩溶形成机理研究.地学前沿,2003,10(2):495~500
赵新锋、陈建耀、唐常源等.珠江三角洲小流域地下水化学特征及演化规律.生态环境,2007,16(6):1620~1626
周绪伦,李振伦等.钙华洞穴的形成与保护.中国岩溶,2004,23(2):107~108
Garrels. R.M.&M.E.Thompson. Geochemical model for sea water at 25℃and one atmosphere total pressure atmosphere. Jour. Sci.,1962. 58~66
Hegeson.H.C. Evaluation of irreversible reations in geoehemical processes involving minerals and aqueous solutions. Thermodynamics relationship.Geochim.Cosmochim,1968. 853~877
陈家玮,鲍征宇.地球化学反应平衡模型的建模方法.地质科技情报,2001,20(3):41~45
邓海琳,涂光炽,李朝阳.水/岩比的地球化学意义.矿物学报,1999,19(3):267~272
郭乙木.《线性与非线性有限元及其应用》.北京:机械工业出版社,2005.156-158
蒋引珊,王玉洁,徐长耀等.热分析法研究矿物分解过程动力学.长春科技大学学报,2000, 30(1):91~93
蒋倩,李心清,丁文慈等.岩溶水对不同环境的响应以及其与大气CO2源汇的关系——基于热力学的研究.矿物岩石地球化学通报,2000,25(3):227~233
李义连,王焰新,周来茹,等.地下水矿物饱和度的水文地球化学模拟分析——以娘子关泉域岩溶水为例.地质科技情报,2003,21(1) :33~35
刘再华.不同CO2分压条件下的白云岩溶解动力学机理.中国科学,2001,31(4):378~384
刘再华.方解石溶解、沉积速率控制的物理、化学机制.重要会议:中国地质科学院“九五”科技成果汇编,2001,253~256
刘再华,W. Dreybrodt,李华举.灰岩和白云岩溶解速率控制机理的比较.地球科学—中国地质大学学报,2006,1(3):412~415
刘再华,W.Dreybrodt,韩军,等. CaCO3-CO2-H2O岩溶系统的平衡化学及分其析.中国岩溶,2005,24(1):1~14
李进文,裴荣富,张荣华,等.安徽铜陵胡村矽卡岩型铜矿水-岩反应动力学实验研究.高校地质学报,2007,13(2):261~271
李义连,王焰新,张江华,等.娘子关泉域岩溶水硫酸盐污染的地球化学模拟分析.地球科学——中国地质大学学报,2000,25(5):467~471
刘兴起,蔡克勤,于升松.基于Pitzer模型的茶卡盐湖卤水和盐类矿物的形成进行了地球化学模拟.中国科学.2003,33(3):247~254
马瑞.水岩相互作用研究-以山西太原为例:[博士学位论文].武汉:中国地质大学,2007
马尧.地球化学模式PHREEQC在地侵工艺中的应用.铀矿冶.2007,26(2):67~71
欧亚波.加拿大Thompson镍矿矿山废水/尾砂作用:[硕士学位论文].四川:成都理工大学,2008
潘根兴,滕永忠.土壤一灰岩岩溶系统中水文地球化学动力学过程模拟及其意义.地球化学,2000,29(3):274~276
乔淑卿,杨作升,潘燕俊,等.沉积物中碳酸盐(方解石和白云石)测定方法的研究.中国海洋大学学报,2006,36(3):484~488
齐雯,郑绵平.西藏扎布耶盐湖卤水演化的Pitzer模型模拟.地质学报,2007,81(12):1735~1741
沈照理,王焰新.水—岩相互作用的回顾与展望.地球科学——中国地质大学学报,2002,27(2):128~133
施枫芝,迟宝明,潘杰.基于化学动力学理论的中深层地下水资源评价.水土保持研究,2007,14(5):211~238
孙连环.塔里木盆地塔中碳酸盐岩储层酸岩反应动力学实验研究.石油与天然气化工.2006,25(1):51~54
陶小晚.不同二氧化碳分压环境近平衡状态下海洋碳酸盐溶解及沉淀动力学研究:[硕士学位论文].北京:中国科学院海洋研究所,2007
童海涛,宋汉周.水-岩作用系统的水文地球化学模拟.水科学进展,2004,15(2):211~215
王丽,王金生,林学钰.水文地球化学模型研究进展.水文地质工程地质,2003(6):105~109
王永新.水-岩相互作用机理及其对库岸边坡稳定性影响的研究:[硕士学位论文].重庆:重庆大学,2006
王益峰,冯祖钧.一个计算机模拟水-岩反应过程的数学模型.矿床地质,1987,6(2):77~88
肖林萍,黄思静.方解石和白云石溶蚀实验热力学模型及地质意义.矿物岩石,2003,23(1):114~116
肖林萍,黄思静.碳酸盐岩溶蚀实验热力学模型及工程地质意义.西南交通大学学报,200237(3):250~253
谢水波,陈泽昂,张晓健,等.铀尾矿库区浅层地下水中U(Ⅵ)迁移的模拟.原子能科学技术,2007,41(1):58~64
徐乐昌.地下水模拟常用软件介绍.铀冶矿,2002,21(1):34~28
于炳松,赖兴运.成岩作用中的地下水碳酸体系与方解石溶解度.沉积学报,2006,24(5):628~635
袁道先,章程.岩溶动力学的理论探索与实践.地球学报,2008,29(3):355~365
杨涛,周文斌.江西科学,2005,23(1):28~33
张国军.夹层软泥化过程及其化学动力学模拟:硕士论文.江苏,河海大学.2007
张文涛,关平等.成岩作用的热力学分析及应用实例.自然科学进展,2007,17(9):1209~1215
张建利,孙忠杰,张泽兰.碳酸盐岩油藏酸盐反应动力学实验研究.油田化学,2003,20(3):216~218
B.Jmerkel B.Planer-Friedrich.水文地球化学模拟程序.(朱义年,王焰新译).武汉:中国地质大学出版社,2005.75~94
周海燕,周训,姚锦梅.广东从化温泉的水文地球化学模拟.现代地质, 2007,21(4):620~623
张凤娥,卢耀如,郭秀红等.符合岩溶形成机理研究.地学前沿,2003,10(2):495~500
赵新锋、陈建耀、唐常源等.珠江三角洲小流域地下水化学特征及演化规律.生态环境,2007,16(6):1620~1626
周绪伦,李振伦等.钙华洞穴的形成与保护.中国岩溶,2004,23(2):107~108