重力归一化总梯度稳定解法_肖鹏飞.pdf

　2006 年 10 月 　　　　　　　　　　　　　　　 　　　　　　　　　　　　　第 4 1 卷 　第 5 期 　 ·非地震 · 重力归一化总梯度的稳定解法 ①② ③ ① ② ① 肖鹏飞 　李 　明 　徐世浙 　孟令顺 　陈生昌 ( ①浙江大学理学院地球科学系 ; ②吉林大学地球探测科学与技术学院; ③中石油勘探开发研究院地质研究所) 摘 　　　要 肖鹏飞 ,李明 ,徐世浙 ,孟令顺 ,陈生昌. 重力归一化总梯度的稳定解法. 石油地球物理勘探 ,2006 ,41 (5) :596 ～ 600 　　重力归一化总梯度法是 20 世纪 80 年代兴起的一种求解重力反演的方法 ,其关键运算是向下延拓及求导 。 原有的向下延拓方法大多数是在波数域实现的 , 由于稳定性差 ,影响了归一化总梯度法的效果 ,从而难以推广 应用 。本文以原先的希尔伯特变换法求解重力归一化总梯度值为基础 ,采用了一种新的、稳定的向下延拓算 法 ———迭代法 ,替代了原有波数域的向下延拓算子 ,提高了计算的精度和稳定性 。模型实验和实测资料应用显 示效果明显改善 。 关键词 　归一化总梯度 　稳定解法 　向下延拓 延拓因子 ,可使 G 法的精度和稳定性均能得到明 H 1 　引言 显提高 。 重力归一化总梯度法是 20 世纪 80 年代兴起的 2 　方法原理 一种求解重力异常的反演方法 。近年来关于归一化 ( ) [ 1 ] [2～5 ] 总梯度法 G 法 的研究越来越深入 ,计算方 2 . 1 　希尔伯特变换 H 法多种多样 。这主要得益于归一化总梯度法的计算 别列兹金归一化总梯度表达式为[ 1 ] 不需要附加条件 ,应用比较简便 ,并能半定量地确定 2 2 V xz ( x , z) + V zz ( x , z) ( ) ( ) G x , z = 1 地质体的位置 ,所以广泛用于研究油气构造或进行 H 1 M - 1 2 2 油气检测 [2 ] 。 M ∑ V xz ( x , z) + V zz ( x , z) j = 0 G 法的出发点在于剩余质量的引力位及导数 式中: x = j d x , d x 为采样间隔; M 为 x 方向采样点 H