基于MWN3D算法的岩石层3D热结构模型的反演方法

在三维热导方程波数域方法(MWN3D)的基础上发展了一种反演岩石层3D热结构的方法.该方法的核心是在使用表面热流数据的同时,假设热生成率在深度方向上为幂指数衰减规律,再通过调节地表等效热生成率的方式,满足地表温度基本上为一常数的强限定条件,最后得到岩石层3D热结构模型以及相应的地表等效热生成率分布.数值模型实验表明这一方法计算速度快,处理问题灵活.在观测数据和测量参数(特别是热生成率)不充分的条件下,该方法不失为一种可行的估计区域岩石层的3D热结构的途径.

A method for inversing lithospheric 3D thermal structures based on MWN3D algorithm

Based on the method of wave-number in three dimensional domain(MWN3D)for the thermal conductive equation,a new method was established to inverse the lithospheric 3D thermal structure.With an assumption of an exponential decay with depth of the heat productions,this method takes the heat flow as upper boundary condition and adjusts the heat production distributions on the earth's surface to satisfy a constant temperature on it.This method can provide a most reasonable lithospheric 3D thermal structure and suggest a distribution of the equivalent heat production at the earth's surface as well.The test models show that this method is characterized by high computation speed and flexibility.In the absence of enough observed data and thermal parameters of the rocks,it is a useful way to investigate lithospheric 3D thermal structure.