1.一种重力数据与大地电磁数据联合反演方法，其特征在于，包括如下步骤：

步骤1)，数据采集步骤：采集工作区所有预先选定的测量点的重力数据和大地电磁数据，其中，所述重力数据是指剩余重力数据，所述大地电磁数据是指电场数据和磁场数据；

步骤2)，数据筛选步骤：分别计算出所述重力数据和大地电磁数据中跳变数据的数据量占各自数据总量的百分比r1和r2，其中，所述跳变数据是指与相邻采集点的数据误差均超过误差阈值的采集点的数据；

步骤3)，初始模型确定步骤：采集工作区的岩石样本至少100个，分别测量样本的密度和电阻率，使用三次样条函数拟合二者之间的关系，由此得到岩石样本的密度模型与岩石样本的电阻率模型之间的对应关系；比较r1与r2，若r1≥r2，则使用均匀半空间密度模型作为重力数据反演初始模型，使用第一次重力数据反演结果的密度模型所对应的电阻率模型作为大地电磁数据反演初始模型，并进行大地电磁数据的第一次反演；若r1＜r2，则使用均匀半空间电阻率模型作为大地电磁反演初始模型，使用第一次大地电磁反演结果的电阻率模型所对应的密度模型作为重力数据反演的初始模型，并进行重力数据的第一次反演；

步骤4)，反演步骤：进行重力数据和大地电磁的联合反演，每进行一次重力数据、大地电磁数据的反演计算并更新密度和电阻率模型后，计算重力数据拟合差Grms和大地电磁数据拟合差MTrms，判断是否同时小于对应的数据拟合差阈值Gth和MTth，若成立，则退出循环，若不成立，则判断迭代次数n是否达到最大的迭代次数Nmax，若成立，则退出循环，若不成立，则继续进行重力数据、大地电磁数据的迭代计算；

所述步骤3)中的岩石样本的密度模型与岩石样本的电阻率模型之间的对应关系根据对所述工作区的岩石样本统计分析确定。

2.如权利要求1所述的重力数据与大地电磁数据联合反演方法，其特征在于，所述步骤

1)中的测量点个数为100-2500个，相邻测量点之间的距离为50-500米。

3.如权利要求1所述的重力数据与大地电磁数据联合反演方法，其特征在于，所述步骤

2)中的误差阈值为20％-40％。

4.如权利要求1所述的重力数据与大地电磁数据联合反演方法，其特征在于，所述步骤

4)的反演步骤具体为：按顺序分别进行重力数据和大地电磁数据的反演迭代计算，其中，从第二次反演迭代开始，利用最新的电阻率模型约束重力数据的反演，利用最新的密度模型约束大地电磁数据的反演。

5.如权利要求4所述的重力数据与大地电磁数据联合反演方法，其特征在于，所述约束为交叉梯度约束。