陈爱华， 许鹤华

（1. 国土资源部海底矿产资源重点实验室 广州海洋地质调查局， 广州 510760； 2. 中国科学院边缘海地质重点实验室 中国科学院南海海洋研究所， 广州 510301）

0 引言

通过分析俯冲板块年龄和热流值， 马尼拉俯冲带属于热俯冲[5-11]。 陈爱华等[12]尝试利用热俯冲探讨马尼拉俯冲带的地震分布规律。 但马尼拉俯冲带的俯冲速度和俯冲角度比其它热俯冲大得多，会对俯冲洋壳的热结构分布产生很大影响[13]。 马尼拉俯冲带经历了复杂的地质演化过程， 俯冲板块深度较深， 同时存在陆壳、 洋陆过渡带和洋壳，且不同区域的地质条件相差较大[14]。 高翔等[15]选取马尼拉俯冲带南、 北两条测线， 根据两处俯冲板块的角度、年龄和运动特征的差异， 分析了摩擦热和剪切热对俯冲带热结构的影响， 讨论俯冲板块发生脱水和部分熔融的位置。 由于研究的问题不同， 建立的模型和考虑的因素也不同， 必然会导致模拟结果产生差异。 近年来， 通过对马尼拉俯冲带以及吕宋岛弧地球物理测量， 获得了马尼拉俯冲带各构造单元莫霍面的形态、 走向、 梯度和地壳性质， 为本文研究马尼拉俯冲带的热结构提供了最新的地球物理资料[16]。 本文综合马尼拉俯冲带的地震、 古地磁、 地热与层析成像等多种最新地球物理资料， 选取3 条典型剖面（位置见图1）， 构建地质模型和二维数值模型， 根据剖面板块俯冲角度、 速度等， 应用有限元方法求解非线性瞬态热方程， 得到马尼拉俯冲带的热结构特征，为进一步探讨马尼拉俯冲带地震成因机制提供热力学依据。……