火山喷发必然伴随着高温岩浆侵入体把热量从岩石圈深部携带到地壳浅层的过程,因而火山区往往也是地热资源富集区.火山岩浆侵入体热扩散过程的数值模拟不仅是火山区构造热演化研究的重要手段,而且也是火山区地热资源评估的重要方法.本文利用在工程热物理领域应用甚广、但在地热地球物理领域鲜有报道的有限元软件ANSYS WORKBENCH对火山活动岩浆冷却过程及其对围岩的热扰动过程进行了数值模拟.首先本文建立一系列简单的柱状模型,模拟了不同规模、不同顶板埋深的侵入体的热扩散过程,引入冷却指数α=(T侵入-T)/(T侵入-T稳态)×100％及干扰温差△T=T-T稳态分别来定量分析岩浆侵入体冷却过程及围岩受侵入体热扩散扰动程度.结果表明,半径为0.5km的柱状侵入体的冷却时间约为50ka左右,对围岩热扰动持续时间不到0.5Ma,最大横向热扰动范围约为3km.而半径1km的柱状侵入体的冷却时间约为200ka左右,对围岩温度场影响可持续1.4Ma,热扰动范围可至7km.顶板埋深越大,侵入体对围岩热扰动的持续时间越短.在上述模拟实验的基础上,本文结合最新的相关研究成果,模拟计算了上地壳岩浆房和上中地壳岩浆房两种情景下长白山天池火山千年大喷发对于温度场的可能影响.结果表明,无论其岩浆房位于上地壳或中地壳上部,千年大喷发对于围岩的热扰动目前还局限于500m左右的范围内,这可能是区内尚无大规模地热显示的原因.
岩浆侵入体 数值模拟 冷却时间 热扰动