[發(fā)明專利]基于測量的重力異常數(shù)據(jù)提取局部重力異常的方法

在審

• 申請?zhí)枺?/td>: 202111055317.1
• 申請日: 2021-09-09
• 公開（公告）號: CN113761457A
• 公開（公告）日: 2021-12-07
• 發(fā)明（設計）人: 孟慶奎;舒晴;高維;周堅鑫;徐光晶;張文志
• 申請（專利權）人: 中國自然資源航空物探遙感中心
• 主分類號: G06F17/10
• 分類號: G06F17/10;G06F17/13;G01V7/00
• 代理公司: 北京云嘉湃富知識產權代理有限公司 11678
• 代理人: 程凌軍

1.一種基于測量的重力異常數(shù)據(jù)提取局部重力異常的方法，其特征在于包括以下步驟：

通過航空、航天、船載或地面平臺測量重力異常數(shù)據(jù)Z(x,y)，該重力異常數(shù)據(jù)Z(x，y)由區(qū)域重力異常R(x，y)和局部重力異常L(x，y)組成；

確定切割半徑r；

根據(jù)重力異常數(shù)據(jù)Z(x，y)和切割半徑r計算與點(x，y)相距r的某4點重力異常的平均值A(x，y)；

計算加權參數(shù)e；

比較區(qū)域重力異常與局部重力異常的曲線彎曲方向；

對于彎曲方向相反的情形，區(qū)域重力異常R_n(x，y)確定為：

ΔZ(x，y)為切割區(qū)域重力異常修正量，第1次切割的區(qū)域重力異常用R₁(x，y)表示，令R₁(x，y)作為下次迭代的重力異常數(shù)據(jù)Z(x，y)，通過該迭代的重力異常數(shù)據(jù)Z(x，y)確定第2次切割的區(qū)域重力異常R₂(x，y)；依次迭代下去，最終有

于是切割半徑r區(qū)域重力異常

計算并提取局部重力異常L(x，y)＝Z(x，y)-R(x，y)。

2.根據(jù)權利要求1所述的基于測量的重力異常數(shù)據(jù)提取局部重力異常的方法，其特征在于計算加權參數(shù)e具體為；

e＝c+d (0≤e≤2)

(0≤c≤1；當時，c＝1)

(0≤d≤1；當時，d＝1)

ΔB_x＝Z(x+r，y)-Z(x-r，y)

ΔB_y＝Z(x，y+r)-Z(x，y-r)。

3.根據(jù)權利要求2所述的基于測量的重力異常數(shù)據(jù)提取局部重力異常的方法，其特征在于所述比較區(qū)域重力異常與局部重力異常的曲線彎曲方向具體包括：計算比較參數(shù)

B_x1＝Z(x+r，y+r)-[Z(x+2r，y+r)+Z(x，y+r)]/2

B_y1＝Z(x+r，y+r)-[Z(x+r，y+2r)+Z(x+r，y)]/2

B_x2＝Z(x+2r，y+r)-[Z(x+4r，y+r)+Z(x，y+r)]/2

B_y2＝Z(x+r，y+2r)-[Z(x+r，y+4r)+Z(x+r，y)]/2

當比較參數(shù)同時滿足下列條件時，表明區(qū)域重力異常和局部重力異常的曲線彎曲方向相反：B_y1×B_x1＞0，B_y1×B_y2＜0，B_x1×B_x2＜0，B_y2×B_x2＞0。

4.根據(jù)權利要求3所述的基于測量的重力異常數(shù)據(jù)提取局部重力異常的方法，其特征在于當比較參數(shù)不能同時滿足下列條件時，表明區(qū)域重力異常和局部重力異常的曲線彎曲方向相同：B_y1×B_x1＞0，B_y1×B_y2＜0，B_x1×B_x2＜0，B_y2×B_x2＞0。

5.根據(jù)權利要求4所述的基于測量的重力異常數(shù)據(jù)提取局部重力異常的方法，其特征在于對于彎曲方向相同的情形，區(qū)域重力異常