吉林大学学报(地球科学版) ›› 2015, Vol. 45 ›› Issue (3): 934-940.doi: 10.13278/j.cnki.jjuese.201503302

• 地球探测与信息技术 • 上一篇    下一篇

对迭代法位场向下延拓方法的剖析

于德武, 龚胜平   

  1. 中国地质科学院地球物理地球化学勘查研究所, 河北 廊坊 065000
  • 收稿日期:2014-06-25 发布日期:2015-05-26
  • 作者简介:于德武(1955),男,教授级高级工程师,主要从事重磁勘探数据处理和解释研究,E-mail:yu_dewu@aliyun.com。
  • 基金资助:

    国家重大科学仪器设备开发专项(2011YQ050060,2011YQ05006011)

Analysis of the Potential Field Downward Continuation Iteration Method

Yu Dewu, Gong Shengping   

  1. Institute of Geophysical & Geochemical Exploration, Chinese Academy of Geological Sciences, Langfang 065000, Hebei, China
  • Received:2014-06-25 Published:2015-05-26

摘要:

位场向下延拓迭代法的实质内容是"向上延拓而不是向下延拓"和以"迭代的结果趋近于观测值"为标准的操作过程。根据数据操作流程剖析了位场向下延拓迭代法的运行机制,得到了迭代数据在空间域的变化规律,即用迭代法将观测高度的位场向下延拓一个深度h。这实际上是通过不同高度的向上延拓来实现的。也就是说,迭代次数增加一次,涉及的上延平面就增大一个h的高度。一般地,迭代次数n与上延高度h的关系为n~(n+1)h。在空间域中,初值、上延结果、差以及每一次校正后的结果都能用满足莱布尼兹定理的交错级数表示,从而得出了迭代法能够收敛的结论;或者,以"观测高度上的实测值与计算值的差值小到可以忽略"为标准,从数学上也能证明迭代法能够收敛。数学推论和模型试验结果说明了迭代的位场初值可以任意给定。在实际操作中,迭代误差标准的影响和由于迭代误差标准不恰当可能出现不能达到迭代标准的情况,需引起注意,也值得进一步研究。

关键词: 位场, 向下延拓, 收敛性, 空间域, 迭代法, 频率域, 快速傅里叶变换

Abstract:

The substance of the potential field downward continuation iteration method is an operating process of "upward continuation rather than downward continuation" based on the criterion of "iteration result tend to be the observation data". We have analyzed the operational mechanism of the method by chasing its operation routine, and obtained the dynamic rule of the data flow in spatial domain. It turned out that the iteration method for a downward continuation leads to an upward continuation of the potential field on the individual levels. The relationship between the iteration period n and the height h of the upward continuation is n~(n+1)h.Shown in an alternating serie of satisfying Leibniz theorem to each iteration result in spatial domain, the convergence of this method is obtained. The convergence can also be verified by using the criterion of minimizing the difference between measured value and calculated data. The mathematical reasoning and numerical computation results show that an initial value of the iteration method can be arbitrary. Also, we have studied the effect of the iteration tolerance criterion, and pointed out the possibility that a tolerance might never be satisfied in practice if the criterion is not fitted.

Key words: potential field, downward continuation, convergence, spatial domain, iteration method, frequency domain, fast Fourier transformation (FFT)

中图分类号: 

  • P631.2

[1] Fedi M,Florio G.A Stable Downward Continuation by Using the ISVD Method[J].Geophysical Journal International,2002,151(1):146-156.

[2] Cooper G.The Stable Downward Continuation of Potential Field Data[J].Exploration Geophysics,2004,35(2):260-265.

[3] 刘保华,焦湘恒,王重德.位场向下延拓的边界单元法[J].石油地球物理勘探,1990,25(4):495-499. Liu Baohua,Jiao Xiangheng,Wang Chongde.Boundary Element Method for Continuation of Potential Field[J]. Oil Geophysical Prospecting, 1990,25(4):495-499.

[4] 毛小平,吴蓉元,曲赞.频率域位场下延的振荡机制及消除方法[J].石油地球物理勘探,1998,33(2):230-237. Mao Xiaoping,Wu Rongyuan,Qu Zan. Oscillation Mechanism of Potential Field Downward Continuation in Frequency Domain and Its Elimination[J]. Oil Geophysical Prospecting, 1998,33(2):230-237.

[5] 徐世浙,沈晓华,邹乐君,等.将航磁异常从飞行高度向下延拓至地形线[J].地球物理学报,2004,47(6):1127-1130. Xu Shizhe,Shen Xiaohua,Zou Lejun,et al.Downward Continuation of Aeromagnetic Anomaly from Flying Altitude to Terrain[J].Chinese Journal of Geophysics,2004,47(6):1127-1130.

[6] 徐世浙.位场延拓的积分-迭代法[J].地球物理学报,2006,49(4):1176-1182. Xu Shizhe.The Integral-Iteration Method for Continuation of Potential Fields[J].Chinese Journal of Geophysics,2006,49(4):1176-1182.

[7] 徐世浙.迭代法与FFT法位场向下延拓效果的比较[J].地球物理学报,2007,50(1):285-289. Xu Shizhe.A Comparison of Effects Between the Iteration Method and FFT for Downward Continuation of Potential Fields[J].Chinese Journal of Geophysics,2007,50(1):285-289.

[8] 徐世浙,余海龙.位场曲化平的插值-迭代法[J].地球物理学报,2007,50(6):1811-1815. Xu Shizhe, Yu Hailong.The Interpolation-Iteration Method for Potential Field Continuation from Undulating Surface to Plane[J].Chinese Journal of Geophysics,2007,50(6):1811-1815.

[9] 王彦国,王祝文,张凤旭,等. 位场向下延拓的导数迭代法[J].吉林大学学报:地球科学版,2012,42(1): 240-245. Wang Yanguo,Wang Zhuwen,Zhang Fengxu,et al. Derivative-Iteration Method for Downward Continuation of Potential Fields[J]. Journal of Jilin University:Earth Science Edition, 2012,42(1): 240-245.

[10] 张志厚,吴乐园. 位场向下延拓的相关系数法法[J].吉林大学学报:地球科学版,2012,42(6): 1912-1919. Zhang Zhihou,Wu Leyuan. Corelation Coefficient Method for Downward Continuation of Potential Fields[J]. Journal of Jilin University:Earth Science Edition, 2012,42(6): 1912-1919.

[11] 刘东甲,洪天求,贾志海,等.位场向下延拓的波数域迭代法及其收敛性[J].地球物理学报,2009,52(6):1599-1605. Liu Dongjia,Hong Tianqiu,Jia Zhihai,et al.Wave Number Domain Iteration Method for Downward Continuation of Potential Fields and Its Convergence[J].Chinese Journal of Geophysics,2009,52(6):1599-1605.

[12] 骆遥.位场迭代法向下延拓的地球物理含义:以可下延异常逐次分离过程为例[J].地球物理学进展,2011,26(4):1197-1200. Luo Yao.The Integral-Iteration Method for Continuation of Potential Fields: A Case by Computing Regional Anomaly Can be Continued Downward[J].Progress in Geophysics,2011,26(4):1197-1200.

[13] 姚长利,李宏伟,郑元满,等.重磁位场转换计算中迭代法的综合分析与研究[J].地球物理学报,2012,55(6):2062-2078. Yao Changli, Li Hongwei, Zheng Yuanman, et al. Research on Iteration Method Using in Potential Field Transformations[J]. Chinese Journal of Geophysics,2012,55(6):2062-2078.

[1] 刘长胜, 高菲. 基于地空频率域电磁法的深部地层极化效应探测方法[J]. 吉林大学学报(地球科学版), 2026, 56(1): 366-376.
[2] 丁乾龙, 沈金松, 陈双全, 冉尚, 龙刚. 频率域传播矩阵法地震AVA多参数反演[J]. 吉林大学学报(地球科学版), 2025, 55(3): 957-969.
[3] 罗新刚, 王万银, .

基于正则化思想的tilt-Euler法在边缘深度反演中的应用 [J]. 吉林大学学报(地球科学版), 2024, 54(2): 633-646.

[4] 陈亮, 张田玉, 王言章, 周海根, 蒋川东, .

基于同步提取变换的地空频率域电磁信号幅度提取方法 [J]. 吉林大学学报(地球科学版), 2024, 54(2): 647-654.

[5] 张建民, 胡英飒, 曾昭发, 杜威, 张领, 槐楠, 胡志鹏. 基于数学形态学的位场数据边缘增强方法[J]. 吉林大学学报(地球科学版), 2019, 49(3): 845-856.
[6] 范雕, 李姗姗, 孟书宇, 邢志斌, 冯进凯. 利用重力异常反演马里亚纳海沟海底地形[J]. 吉林大学学报(地球科学版), 2018, 48(5): 1483-1492.
[7] 王泰涵, 黄大年, 马国庆, 李野, 林松. 基于并行预处理算法的三维重力快速反演[J]. 吉林大学学报(地球科学版), 2018, 48(2): 384-393.
[8] 曾昭发, 霍祉君, 李文奔, 李静, 赵雪宇, 何荣钦. 任意各向异性介质三维有限元航空电磁响应模拟[J]. 吉林大学学报(地球科学版), 2018, 48(2): 433-444.
[9] 杜威, 许家姝, 吴燕冈, 郝梦成. 位场垂向高阶导数的Tikhonov正则化迭代法[J]. 吉林大学学报(地球科学版), 2018, 48(2): 394-401.
[10] 李光, 渠晓东, 黄玲, 方广有. 基于磁偶极子的频率域电磁系统几何误差分析[J]. 吉林大学学报(地球科学版), 2017, 47(4): 1255-1267.
[11] 冯晅, 鲁晓满, 刘财, 周超, 金泽龙, 张明贺. 基于逐减随机震源采样法的频率域二维黏滞声波方程全波形反演[J]. 吉林大学学报(地球科学版), 2016, 46(6): 1865-1873.
[12] 高成, 孙建国. 不同域的局部平面波分解应用与对比[J]. 吉林大学学报(地球科学版), 2015, 45(5): 1523-1529.
[13] 王卫平, 曾昭发, 李静, 吴成平. 频率域航空电磁法地形影响和校正方法[J]. 吉林大学学报(地球科学版), 2015, 45(3): 941-951.
[14] 郭振波,李振春. 起伏地表条件下频率-空间域声波介质正演模拟[J]. 吉林大学学报(地球科学版), 2014, 44(2): 683-693.
[15] 马国庆,黄大年,杜晓娟,李丽丽. Hartley变换在位场(重、磁)异常导数计算中的应用[J]. 吉林大学学报(地球科学版), 2014, 44(1): 328-335.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 于 平,李瑞磊,付 雷,郝 雪,张向军,廉国芬. 松辽盆地滨北地区区域构造特征及意义--地震长剖面给出的证据[J]. J4, 2005, 35(05): 611 -615 .
[2] 张原庆, 宋炳忠, 王玉福, 张宁. 鲁西铜石岩体金成矿规律和成矿预测[J]. J4, 2010, 40(6): 1287 -1294 .
[3] 李建平,李桐林,张 辉,徐凯军. 不规则回线源层状介质瞬变电磁场正反演研究及应用[J]. J4, 2005, 35(06): 790 -0795 .
[4] 高松, 宋莹, 王琳,蒋步新. 水体中秋兰姆特效降解菌的筛选及其降解性能研究[J]. J4, 2006, 36(03): 455 -457 .
[5] 张凤君,李 卿,马玖彤,于广菊. 膜蒸馏处理糠醛废水的实验研究[J]. J4, 2006, 36(02): 270 -0273 .
[6] 卢双舫,李吉君,薛海涛,徐立恒. 油成甲烷碳同位素分馏的化学动力学及其初步应用[J]. J4, 2006, 36(05): 825 -829 .
[7] 丁志宏,冯平,毛慧慧. 考虑径流年内分布影响的丰枯划分方法及其应用[J]. J4, 2009, 39(2): 276 -0280 .
[8] 任何军, 刘娜,高松,张兰英,张玉玲,周睿. 假单胞菌DN2对多氯联苯的降解及bphA1核心序列测定[J]. J4, 2009, 39(2): 312 -0316 .
[9] 黄奇波, 覃小群, 刘朋雨, 康志强, 唐萍萍. 半干旱区岩溶碳汇原位监测方法适宜性研究[J]. 吉林大学学报(地球科学版), 2015, 45(1): 240 -246 .
[10] 熊晓亮,孙红月,张世华,蔡岳良. 高扬程虹吸保障条件分析与合理管径选择数值模拟[J]. 吉林大学学报(地球科学版), 2014, 44(5): 1595 -1601 .