本发明为基于张量分解的卫星磁场数据融合地震异常提取方法，包括：多通道卫星磁场数据减去CHAOS‑7磁场模型数据去除地球主磁场和静态岩石圈磁场；利用小波分解去除剩余数据的低频成分；根据同步挤压小波变换获取的各通道磁场数据的时频幅值矩阵构造对应的三阶时频谱张量；利用非负张量分解对构造的三阶时频谱张量进行分解，分解出与地震相关的融合局部特征分量；根据能量‑熵比选取融合局部特征分量，通过阈值分割法方法提取该分量中的融合地震异常，并输出融合地震异常结果。本发明结合了多通道卫星磁场数据的信息，从中获取了与地震相关的融合局部特征分量，提取的融合地震异常更加可靠。