【技术实现步骤摘要】
基于重磁梯度数据张量不变量的目标深度估计方法及系统
本专利技术属于重力勘探领域，尤其涉及基于重磁梯度数据张量不变量的目标深度估计方法及系统。
技术介绍
本部分的陈述仅仅是提供了与本公开相关的
技术介绍
信息，不必然构成在先技术。重力勘探起初是用来研究地球形状的重力学测量方法，是人类最早于用来研究地球的手段，如今重力勘探仍然是地球物理勘探中一个重要的分支，其目的是测量地下由密度分布不均和差异引起的重力或重力梯度异常，进而用来研究和推断地下地质体空间和密度分布属性等，更好地认识地下空间并服务于人类社会。重力勘探本身具有快速有效、轻便和成本低的优点，被广泛地应用于研究大地构造、勘探矿产和油气源分布，以及工程与环境调查领域。专利技术人在研究中发现，地质目标体的深度位置计算是重磁数据解释的重点问题，而位场数据具有较低的深度分辨率特点且易受磁倾角和磁偏角的影响，深度估计误差较大，是一个难点问题；而若采用二维剖面数据进行处理解释时，其具有局限性且易受噪声干扰的影响；而通常采用重力异常或垂向重力异常数据处理解释，其信息量少，数据...
【技术保护点】
1.基于重磁梯度数据张量不变量的目标深度估计方法，其特征是，包括：/n构建引力位在球体模型三个方向上的二阶梯度数据平面图，根据重力梯度张量数据，计算其张量矩阵的特征值及特征值组合的张量不变量；/n针对球体模型，根据特征值组合后的张量不变量与总水平导数，计算张量不变量与总水平导数比值的反正切值，建立水平距离与目标深度之间的关系；/n根据反正切值平面图上的两条45°等值线之间的距离，对目标进行深度估计。/n
【技术特征摘要】
1.基于重磁梯度数据张量不变量的目标深度估计方法，其特征是，包括：
构建引力位在球体模型三个方向上的二阶梯度数据平面图，根据重力梯度张量数据，计算其张量矩阵的特征值及特征值组合的张量不变量；
针对球体模型，根据特征值组合后的张量不变量与总水平导数，计算张量不变量与总水平导数比值的反正切值，建立水平距离与目标深度之间的关系；
根据反正切值平面图上的两条45°等值线之间的距离，对目标进行深度估计。
2.如权利要求1所述的基于重磁梯度数据张量不变量的目标深度估计方法，其特征是，重力梯度张量数据具体为：重力场是引力位在笛卡尔坐标系x,y,z方向上的矢量分量，重力梯度张量为地质体的引力位的二阶导数，即重力三分量Gx,Gy,Gz在x,y,z三个方向上的导数，其中Gx在x，y，z方向的一阶梯度分别为gxx,gxy,gxz，Gy在x，y，z方向的一阶梯度分别为gyx,gyy,gyz，Gz在x，y，z方向的一阶梯度分别为gzx,gzy,gzz，可构成梯度张量矩阵：
其中矩阵T为对称矩阵，且重力位满足拉普拉斯方程，即有gxy＝gyx，gxz＝gzx，gyz＝gzy，gxx+gyy+gzz＝0，即矩阵T有5个独立的分量。
3.如权利要求2所述的基于重磁梯度数据张量不变量的目标深度估计方法，其特征是，由于重力梯度张量矩阵T为一个实对称阵，通过计算特征值将其对称化，即张量矩阵T可以通过坐标转换成对称阵Λ：
ETTE＝Λ
其特征值Λ和特征向量E分别为：
E＝[e1e2e3]
其中，特征向量E表示新笛卡尔坐标系的正交基底，在任何坐标转换下，矩阵T将包含三个不变量要素：
I0＝trace(T)＝gxx+gyy+gzz＝0
可以通过三个不变量要素，得到求解特征值的参数方程：
λ3-I0λ2+I1λ-I2＝0
求其方程的根，得到三个特征值λ1,λ2,λ3为
λ1＝A+B(1)
其中，
即三个特征值可利用重力梯度张量矩阵中的各个分量表示。
4.如权利要求3所述的基于重磁梯度数据张量不变量的目标深度估计方法，其特征是，计算三个特征值组合后的不变量μ，即
其中不变量μ可利用重力梯度张量的各个分量表示，且不随坐标转换而变化的标量，标量值与坐标系的旋转无关，因此也与参考基无关。
5.如权利要求3所述的基于重磁梯度数据张量不变量的目标深度估计方法，其特征是，针对球体模型，假设均匀密度的球体在地面的投影坐标为(x0,y0,0)，球体中心埋深为h，半径为R，万有引力常数为G，剩余密度(与周围岩体密度差)为ρ，则剩余质量m＝4πR3ρ/3，在地面观测点产生的各个重力梯度张量分量异常为：
...
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。