一种基于卫星测高重力数据的局部重力异常提取方法技术

本发明专利技术公布了一种基于卫星测高重力数据的局部重力异常提取方法，其基本步骤为：利用sobel算子处理卫星测高重力数据，得到各重力异常区域；根据各重力异常区域的边界线，提取出各重力异常区域的最小外接矩形；将各重力异常区域的最小外接矩形进行旋转坐标系得到变换矩形，利用改进后的曲面拟合法得到各变换矩形区域内的重力异常值；利用网函数插值处理各变换矩形的重力异常值，得到各重力异常区域的区域异常重力值；最终获得各重力异常区域的局部重力异常值。本发明专利技术的方法科学合理，具有较强的可操作性和实用性，利用网函数插值法对卫星测高重力数据进行处理得到局部重力异常信息。

A local gravity anomaly extraction method based on satellite altimetry gravity data

【技术实现步骤摘要】
一种基于卫星测高重力数据的局部重力异常提取方法
本专利技术涉及网函数插值方法，具体说是一种基于卫星测高重力数据的局部重力异常提取方法。
技术介绍
一般获得的重力数据为布格重力异常数据，它是当代重力研究中获取最为简单且研究最为广泛的重力异常数据。布格重力异常是由区域重力异常和局部重力异常叠加组合而成。其中，由分布范围较大且影响较深的地质因素(例如：区域地质构造)所引起的重力异常称为区域重力异常，由地质因素(例如：矿产，油气)所引起的小范围的重力异常称为局部重力异常。由此可见，局部重力异常的提取不仅是资源勘探的重要研究内容，同时也有助于地质构造的研究。在此基础下，如何有效地提取出局部重力异常也就成为重要问题。目前用于提取局部重力异常的方法有很多，其中应用最为广泛的方法就是圆周法与趋势分析法。两种方法的基本原理基本一致，皆是通过计算点位的区域重力异常值，再用重力异常值减去区域重力异常值得到局部重力异常值。这两种方法的实验原理较为简易但是在实验过程中会有较多局限性。例如，会出现边缘缺失或是边缘拟合精度低的现象，由此可见，传统提取方法很难有效地提本文档来自技高网...

【技术保护点】
1.一种基于卫星测高重力数据的局部重力异常提取方法，其特征在于，包括以下步骤：/n(1)利用重力数据G和sobel算子得到重力梯度数据G′，进而得到P个重力异常区域；/n(2)在各重力异常区域中，提取重力异常区域边界线，得到各重力异常区域的最小外接矩形ABCD和与之对应的旋转矩阵R(θ)。将各最小外接矩形ABCD经过旋转变换得到变换矩形G

【技术特征摘要】
1.一种基于卫星测高重力数据的局部重力异常提取方法，其特征在于，包括以下步骤：
(1)利用重力数据G和sobel算子得到重力梯度数据G′，进而得到P个重力异常区域；
(2)在各重力异常区域中，提取重力异常区域边界线，得到各重力异常区域的最小外接矩形ABCD和与之对应的旋转矩阵R(θ)。将各最小外接矩形ABCD经过旋转变换得到变换矩形G1G2G3G4；
(3)利用步骤(2)得到的变换矩形G1G2G3G4，通过改进的曲面拟合法计算各变换矩形G1G2G3G4内的重力异常值Gk；
(4)将步骤(3)得到的变换矩形G1G2G3G4进行网函数插值，得到插值g。利用旋转矩阵的逆矩阵R(θ)′，将变换矩形G1G2G3G4旋转回原最小外接矩形ABCD，并从插值g中提取出重力异常区域的区域重力异常值gk；
(5)利用步骤(1)中的原始重力数据G与步骤(4)得到的区域异常值gk，计算各重力异常区域的局部重力异常值△gk。


2.根据权利要求1所述的一种基于卫星测高重力数据的局部重力异常提取方法，其特征在于：所述步骤(2)中构建变换矩形G1G2G3G4过程为：
A.提取重力异常区域的边界，计算边界每两相邻点的方位角θi′，方位角θi′对取余得到倾斜角θi(i＝1，2，3...)，构建各倾斜角的倾斜矩阵R(θi)(i＝1，2，3...)；
B.根据步骤A中两相邻点，做该重力异常区域的外接矩形AiBiCiDi(i＝1，2，3...)，记录外接矩形AiBiCiDi(i＝1，2，3...)的面积Si(i＝1，2，3...)，其中面积最小的外接矩形为最小外接矩形ABCD，其对应的旋转矩阵为R(θ)；<...

【专利技术属性】
技术研发人员：刘善伟，叶传龙，孙钦婷，
申请(专利权)人：中国石油大学华东，
类型：发明
国别省市：山东;37