高温下钙蒙脱石膨胀特性的分子动力学模拟  

Molecular dynamics simulation of swelling properties of Ca-montmorillonite at high temperatures

作　　者：杨亚帆 王建州[1] 商翔宇[1] 王涛[1] 孙树瑜 Yang Ya-Fan;Wang Jian-Zhou;Shang Xiang-Yu;Wang Tao;Sun Shu-Yu(State Key Laboratory for Geomechanics and Deep Underground Engineering,China University of Mining and Technology,Xuzhou 221116,China;Physical Science and Engineering Division,King Abdullah University of Science and Technology,Thuwal 23955-6900,Saudi Arabia;Institute of Geophysics and Geomatics,China University of Geosciences,Wuhan 430074,China)

机构地区：[1]中国矿业大学深部岩土力学与地下工程国家重点实验室,徐州221116 [2]阿卜杜拉国王科技大学物理科学与工程学部,图瓦23955-6900 [3]中国地质大学(武汉)地球物理与空间信息学院,武汉430074

出　　处：《物理学报》2022年第4期42-53,共12页Acta Physica Sinica

基　　金：国家自然科学基金重点项目(批准号:51936001);国家自然科学基金(批准号:51874262);阿卜杜拉国王科技大学(批准号:BAS/1/1351-01,URF/1/4074-01,URF/1/3769-01);中国矿业大学引进人才科研启动经费(批准号:102521155)资助的课题。

摘　　要：高温下蒙脱石的膨胀特性在核废料深部封存、二氧化碳封存及页岩气开发等应用中有着重要影响,但相关机理尚不明确.本工作使用分子动力学模拟为技术手段计算5 MPa和298—500 K等条件下,1.40—4.00 nm晶面间距(d)的一系列饱和钙蒙脱石的膨胀压力.以模拟所得的数值结果为依据,基于水化效应、双电层效应和离子关联效应等模型推演膨胀压力随温度与d的变化规律,并与相应的实验数据进行对比.模拟结果表明,当d较小时,因为高温会弱化水化力的强度,钙蒙脱石膨胀压力震荡的幅度降低,同时水化力作用的d的范围减小.当d较大时,因为高温强化离子关联效应,膨胀压力降低,同时双电层力的作用的d的范围增加.在较高温度和较大d时,膨胀压力为收缩力,阻碍膨胀.这些膨胀压力的变化规律与前期钠蒙脱石体系的研究类似.然而,通过对比两种蒙脱石体系的模拟结果,发现两种体系存在显著的差异—钙蒙脱石比钠蒙脱石更难膨胀到较大的d.此模拟结果与前人实验观测的结果相符.我们进一步将此差异归于钙蒙脱石的离子关联效应要远大于钠蒙脱石.有别于分子模拟中对于离子关联效应的精确描述,连续化的Poisson-Boltzmann方程因为忽略了离子关联效应,从而无法表达出与两种体系模拟结果都相吻合的膨胀压力变化规律.The swelling of Ca-montmorillonite at elevated temperatures is important for many applications including geological disposal of radioactive waste,subsurface carbon sequestration,and shale gas exploration.However,the experimentally observed swelling behaviors of Ca-montmorillonite contacting liquid water and the temperature effects on the swelling pressure are not well understood.In this work,molecular dynamics simulations are carried out to study the swelling of Wyoming Ca-montmorillonite with a d-spacing(d)range of 1.40–4.00 nm at 5 MPa and various temperatures(298–500 K).The ClayFF and SPC are adopted for modeling Ca-montmorillonite and water,respectively.The simulation box is measured to be 11.15,3.66,and 28.00 nm in the x-,y-,and z-direction.Atomistic pistons are used to control the bulk pressure of the water environment,and the implicit walls are implemented for preventing the ions from leaking from the pore into the water environment.The clay atoms are fixed during the simulation and the swelling pressure is calculated through dividing the force by the area.The equilibrium time is at least 20 ns and the production time falls in a range of 50–88 ns.The swelling pressure results show that for small d,high temperature reduces the magnitude of the oscillating curve of swelling pressure and also reduces the range of d where hydration force dominates the swelling pressure.This temperature effect is due to the weakened hydration force as evidenced from the weakened water density distributions inside the pore.For large d,high temperature reduces the swelling pressure,which is consistent with the experimental result,and increases the range of d where double layer force dominates the swelling pressure.The reduction of the swelling pressure can be explained by the enhanced ion correlation that reduces the double layer force according to the strong coupling theory,given that the calculated coupling parameters at higher temperatures are smaller.The swelling pressures are negative at elevated temperatures and large

关 键 词：分子动力学模拟 膨胀压力 钙蒙脱石 离子关联效应 

分 类 号：O643.1[理学—物理化学]