岩浆侵入体冷却过程数值模拟研究.pdfVIP

中国地球科学联合学术年会 2014 ·429 · 岩浆侵入体冷却过程的数值模拟研究 段文涛 1 黄少鹏 1,2* 时庆金 1 陈 希 1 1.西安交通大学全球环境变化研究院地热与环境研究实验室，西安，710049；2.Department of Earth and Environmental Sciences, University of Michigan, Ann Arbor, MI 48109-1005, USA 地热能作为新能源家族中的重要一员，在经济性、环保性和稳定性上具有其他新能源无法比拟的优 势。理论上在任何一个地方达到一定深度都可以获得所需要的地热能，新生代火山构造区更是寻找高温 地热资源的首选区，因为火山喷发必然伴随着高温岩浆侵入体把热量从岩石圈深部携带到地壳浅层的过 程，我们利用有限元数值模拟软件 ANSYS 建立了火山岩浆活动概念模型，模拟计算岩浆侵入体的冷却 过程，为在火山活动构造区地热能源开发提供参考数据。所选取地质数值模型直径 40km 、深度40km， 假设火山活动的岩浆通道为位于模型中心的直立圆柱体，岩浆侵入体正上方另设一底部直径 2.4km、顶 部直径 2km、高 1km 的圆台以模拟火山锥。模型下边界传导热流边界设定为40mW/m2 ，上边界地表温 度为 10℃，地面温度随高程的变化率为 5℃/km。计算模型上（0-15km ）、中（15-30km）、下（30-40km ） -1 -1 -1 -1 地壳岩石和侵入体的热导率分别为 2.6 、2.5 、2.8 和 2.8 W·m ·K , 比热 1050、970、950 和 1100 J·kg ·K ， 密度 2750 、2900 、3200 和 3200 kg·m-3 ，生热率 1.5、0.4、0.02 和 0.02 μW·m-3 。 本文着重分析不同规模和埋深的岩浆侵入体的冷却过程，计算了埋藏深度分别为 0km、1km、3km、 6km 和 11km，岩浆通道半径分别为 0.5km 和 1km 的案例，对比了 10 组模拟结果。为了描述方便，我 们引入冷却因子α= （T -T)/T -T )×100% 来衡量岩浆侵入体逐渐冷却最终趋于与围岩达到温度平衡的程 1 1 0 度，其中 T1 为岩浆侵入体初始温度，设定为 1000℃，T0 为稳态温度，T 为冷却过程温度。由于径向的 传热方式是由中心向外的，因而在确定侵入体某一深度的冷却时间时，我们以该深度岩浆侵入体中心处 衰减因子达到 90%为准。 模拟计算结果（见表 1）表明，岩浆的冷却时间与侵入体的规模（半径）平方成正比，与顶板埋深 成反比，与侵入深度成反比，但到一定深度趋于定值。比如，对于半径为 500m、顶板埋深为 1km 的侵 入体，深度 5.0km、6.0km、8.0km 和 11km 处的岩浆冷却时间分别大致为 20.5ka 、21.2ka 、21.5ka 和 21.5ka ； 半径为 1km、顶板埋深同样为 1km 的侵入体，深度 5.0km、6.0km、8.0km 和 11km 处的岩浆冷却时间 分别约为 83.9ka、86.0ka、87.2ka 和 87.3ka；而对于半径为 1km、侵入深度为 3km 的侵入体，深度 5.0km、 6.0km、8.0km 和 11km 处的岩浆冷却时间分别为 57.2ka、64.8ka、86.6ka 和 86.9ka。 表 1 火山活动侵入岩浆达到冷却因子为 90%所需要的时间（单位 ka ） 岩浆侵入柱体顶板埋深/km 深度 侵入体半径 0.5km 侵入体半径 1.0km