专利摘要
本发明公开提供了一种基于重磁梯度数据张量不变量的目标深度估计方法及系统。根据球体模型对象的重力梯度张量平面数据图，计算其梯度矩阵的特征值及特征值组合的张量不变量；利用张量不变量与总水平导数，通过反正切运算，建立水平距离与埋深之间的关系；根据两条45°等值线之间的距离，对目标体的深度进行估计。该方法对于单独地质体或距离较远的多个地质体情况下的深度估计问题，可以获得较高精度的计算结果。
权利要求
1.基于重磁梯度数据张量不变量的目标深度估计方法，其特征是，包括：
构建引力位在球体模型三个方向上的二阶梯度数据平面图，根据重力梯度张量数据，计算其张量矩阵的特征值及特征值组合的张量不变量；
针对球体模型，根据特征值组合后的张量不变量与总水平导数，计算张量不变量与总水平导数比值的反正切值，建立水平距离与目标深度之间的关系；
根据反正切值平面图上的两条45°等值线之间的距离，对目标进行深度估计；
所述张量矩阵T可以通过坐标转换成对称阵Λ：
E
其特征值Λ和特征向量E分别为：
E＝[e
其中，特征向量E表示新笛卡尔坐标系的正交基底，在任何坐标转换下，矩阵T将包含三个不变量要素：
I
可以通过三个不变量要素，得到求解特征值的参数方程：
λ
求其方程的根，得到三个特征值λ
λ
其中，
即三个特征值可利用重力梯度张量矩阵中的各个分量表示；
针对球体模型，假设均匀密度的球体在地面的投影坐标为(x0,y0,0)，球体中心埋深为h，半径为R，万有引力常数为G，剩余密度为ρ，则剩余质量m＝4πR
将球体重力梯度张量分量公式(5)～(10)带入三个特征值求解公式(1)～(3)中，并计算其组合后的张量不变量μ，有
则
令总水平导数为TDX，并带入公式(5)、(6)，即
构造公式ET并根据公式(11)、(12)简化，
简化后：
令
通过公式ET建立水平距离与埋深之间的关系，且与目标剩余质量无关；
当ET为45°时，即r
其方程解r
即
利用ET的两条45°等值线之间距离的0.4472倍进行深度估计。
2.如权利要求1所述的基于重磁梯度数据张量不变量的目标深度估计方法，其特征是，重力梯度张量数据具体为：重力场是引力位在笛卡尔坐标系x,y,z方向上的矢量分量，重力梯度张量为地质体的引力位的二阶导数，即重力三分量G
其中矩阵T为对称矩阵，且重力位满足拉普拉斯方程，即有g
3.如权利要求1所述的基于重磁梯度数据张量不变量的目标深度估计方法，其特征是，计算三个特征值组合后的不变量μ，即
其中不变量μ可利用重力梯度张量的各个分量表示，且不随坐标转换而变化的标量，标量值与坐标系的旋转无关，因此也与参考基无关。
4.一种计算装置,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征是，所述处理器执行所述程序时实现权利要求1-3任一项方法的步骤。
5.一种计算机可读存储介质，其上存储有计算机程序，其特征是，该程序被处理器执行时执行权利要求1-3任一项方法的步骤。
6.基于重磁梯度数据张量不变量的目标深度估计系统，其特征是，包括服务器，所述服务器被配置为执行权利要求1-3任一项方法的步骤。
基于重磁梯度数据张量不变量的目标深度估计方法及系统专利购买费用说明
Q：办理专利转让的流程及所需资料
A：专利权人变更需要办理著录项目变更手续，有代理机构的，变更手续应当由代理机构办理。
1：专利变更应当使用专利局统一制作的“著录项目变更申报书”提出。
2：按规定缴纳著录项目变更手续费。
3：同时提交相关证明文件原件。
4：专利权转移的，变更后的专利权人委托新专利代理机构的，应当提交变更后的全体专利申请人签字或者盖章的委托书。
Q:专利著录项目变更费用如何缴交
A：（1）直接到国家知识产权局受理大厅收费窗口缴纳,（2）通过代办处缴纳，（3）通过邮局或者银行汇款,更多缴纳方式
Q：专利转让变更，多久能出结果
A：著录项目变更请求书递交后，一般1-2个月左右就会收到通知，国家知识产权局会下达《转让手续合格通知书》。
动态评分
0.0