空间域位场延拓新方法研究

doi:10.6038/j.issn.1004-2903.2012.04.026

摘要/Abstract

摘要： 在空间域进行位场延拓,需要数值求解第一类Fredholm积分方程,由于所得方程组系数矩阵不是稀疏矩阵,求解该方程组需要的计算机内存大,计算量大,导致延拓算法在一般计算机上难以实现,阻碍了对空间域位场延拓方法的研究.在分析系数矩阵结构特征的基础上,本文证明了方程组系数矩阵是对称的分块Toeplitz型矩阵.利用系数矩阵的对称性和分块Toeplitz型矩阵与向量相乘的快速算法,解决了系数矩阵的存储和计算问题,使得空间域位场延拓成为可能,为研究新的位场延拓方法和分析延拓误差提供了一条新的途径.利用模型数据和实测资料,对空间域位场向上延拓、空间域积分迭代法向下延拓进行了检验,结果证实了空间域位场延拓的可行性和正确性.

Abstract: To apply continuation of potential field in spatial domain, Fredholm integral equation of the first kind has to be solved numerically. Because the coefficient matrix that comes from discretizing integral equation is not sparse and has very large scale, it requires large computer memory to save the matrix, and operation on it consumes very long time. The continuation method in spatial domain is not developed in practice. In this paper, based on the analysis of structure of coefficient matrix, we find an important fact that the coefficient matrix is a symmetric, block Toeplitz type matrix. A mathematical proof is given in the paper. By using fast algorithm that already exists to solve block Toeplitz system, we realize spatial domain upward continuation and downward continuation. Model data and aeromagnetic data test shows the feasibility and correctness of spatial continuation of potential field.

中图分类号:

P631

陈龙伟, 胡小平, 吴美平, 马涛. 2012. 空间域位场延拓新方法研究. 地球物理学进展, (4): 1509-1518. doi: 10.6038/j.issn.1004-2903.2012.04.026 .

CHEN Long-wei, HU Xiao-ping, WU Mei-ping, MA Tao. 2012. Research on spatial domain continuation method for potential field. Progress in Geophysics. (4): 1509-1518. doi: 10.6038/j.issn.1004-2903.2012.04.026 .

0
/ / 推荐

导出引用管理器 EndNote|Reference Manager|ProCite|BibTeX|RefWorks

链接本文: http://www.progeophys.cn/CN/10.6038/j.issn.1004-2903.2012.04.026

http://www.progeophys.cn/CN/Y2012/V/I4/1509

参考文献

[1] 侯重初. 位场的频率域向下延拓方法. 物探与化探, 1982, 1: 33-40. Hou Z C. Frequency domain downward continuation method for potential field. Geophysical & Geochemical Exploration, 1982, 1: 33-40.
[2] 栾文贵. 场位解析延拓的稳定化算法. 地球物理学报, 1983, 26 (3): 263-274. Luan W G. The stabilized algorithm of the analytic continuation for the potential field. Chinese J. Geophys. (in Chinese), 1983, 26(3): 263-274.
[3] 王延忠, 熊光楚. 位场向下延拓组合滤波器的设计和应用. 地球物理学报, 1985, 28 (5): 537-543. Wang Y Z, Xiong G C. A combined filter used for the downward continuation of the potential field. Chinese J. Geophys. (in Chinese), 1985, 28(5): 537-543.
[4] 梁锦文. 位场向下延拓的正则化方法. 地球物理学报, 1989, 32(3): 600-608. Liang J W. Regularization methods for potential field downward continuation.Chinese J. Geophys. (in Chinese), 1989, 32(3): 600-608.
[5] Xia J, Sprowl D R, Adkins-Heljeson D. Correction of Topographic Distortions in Potential-Field Data: A Fast and Accurate Approach. Geophysics, 1993, 58(4): 515-523.
[6] Naidu P S, Mathew M P. Fast Reduction of Potential Fields Measured over an Uneven Surface to a Plane Surface. IEEE Transactions on Geoscience and Remote Sensing, 1994, 32(3): 508-512.
[7] 毛小平, 吴蓉元, 曲赞. 频率域位场下延的振荡机制及消除方法. 石油地球物理勘探, 1998, 33(2): 230-237. Mao X P, Wu R Y, Qu Z. Oscillation mechanism of potential field downward continuation in frequency domain and its elimination. OGP, 1998, 33(2): 230-237.
[8] Fedi M, Florio G. A Stable Downward Continuation by Using the ISVD Method. Geophysical Journal International, 2002, 151: 146-156.
[9] 徐世浙. 迭代法与FFT法位场向下延拓效果的比较. 地球物理学报,2007, 50(1): 285-289. Xu S Z. A comparison of effects between the iteration method and FFT for downward continuation of potential fields. Chinese J. Geophys. (in Chinese), 2007, 50(1): 285-289.
[10] 陈生昌, 肖鹏飞. 位场向下延拓的波数域广义逆算法. 地球物理学报, 2007, 50(6): 1816-1822. Chen S C, Xiao P F. Wavenumber domain generalized inverse algorithm for potential field downward continuation. Chinese J. Geophys. (in Chinese), 2007, 50(6): 1816-1822.
[11] 张辉, 陈龙伟, 任治新等. 位场向下延拓迭代法收敛性分析及稳健向下延拓方法研究. 地球物理学报, 2009, 52(4): 1107-1113. Zhang H, Chen L W, Ren Z X, et al. Analysis on convergence of iteration method for potential fields downward continuation and research on robust downward continuation method. Chinese J. Geophys. (in Chinese), 2009, 52(4): 1107-1113.
[12] 刘东甲, 洪天求, 贾志海等. 位场向下延拓的波数域迭代法及其收敛性. 地球物理学报, 2009, 52(6): 1599-1605. Liu D J, Hong T Q, Jia Z H, et al. Wave number domain iteration method for downward continuation of potential fields and its convergence. Chinese J. Geophys. (in Chinese), 2009, 52 (6): 1599-1605.
[13] Bath M 著, 郑治真等译. 地球物理学中的谱分析. 北京: 地震出版社, 1978. Bath M, translated by Zheng Z Z,et al. Spectral analysis in Geophysics. Beijing: Seism Press, 1978.
[14] 熊光楚. 磁(重力)异常的变换及滤波技术. 北京: 冶金工业出版社, 1990. Xiong G C. Transform and filtering technique for Magnetic anomaly, Beijing: Metallurgical Press, 1990.
[15] 熊光楚. 有限长及离散化数据对重、磁异常在频率域作变换处理的影响. 桂林冶金地质学院学报, 1982, (1): 23-34. Xiong G C. The effect of frequency domain transformation of discrete gravity and magnetic anomalies with finite length. Journal of Guilin College of Geology, 1982, (1): 23-34.
[16] 熊光楚. 在地球物理测量的数据处理中值得进一步研究的一些问题. 物探与化探, 1985, 9(2): 100-107.. Xiong G C. Some problems worthy of further study in geophysical data processing. Geophysical & Geochemical Exploration, 1985, 9(2): 100-107.
[17] 吴宣志, 薛光琦. 谱分析方法正演位场的若干问题. 物探与化探, 1987, 11 (4): 298-306. Wu X Z, Xue G Q. Some aspects on calculation of potential fields using Fourier transformation. Geophysical & Geochemical Exploration, 1987, 11 (4): 298-306.
[18] Ricard Y, Blakely R J. A Method to Minimize Edge Effects in Two-Dimensional Discrete Fourier Transform. Geophysics, 1988, 53(8): 1113-1117.
[19] 柴玉璞. 位场dft算法研究. 地球物理学报, 1988, 31(2): 211-224. Chai Y P. Algorithm investigation of DFT of potential field. Chinese J. Geophys. (in Chinese), 1988, 31(2): 211-224.
[20] 柴玉璞. Fourier变换数值计算的偏移抽样理论. 中国科学(E辑), 1996, 26 (5): 450-456. Chai Y P. Shift sampling theory about Fourier transform computation. Science in China (Series E), 1996, 26(5): 450-456.
[21] Huestis S P, Parker R L. Upward and Downward Continuation as Inverse Problems. Geophysical Journal Research, 1979, 57: 171-188.
[22] 杨文采. 用于位场数据处理的广义反演技术. 地球物理学报, 1986, 29 (3): 283-291. Yang W C. A generalized inversion technique for potential field data processing. Chinese J. Geophys. (in Chinese), 1986, 29 (3): 283-291.
[23] Huestis S P. The Continuation Inverse Problems Revisited. Geophysical Journal International, 1998, 133: 705-712.
[24] Curtis R. Vogel. Computational Methods for Inverse Problems. SIAM, Philadelphia, 2002.
[25] Zhou J, Shi G G, Ge Z L. Study of Iterative Regularization Methods for Potential Field Downward Continuation in Geophysical Navigation. journal of Astronautics, 2011, 32(4): 787-794.
[26] 曾小牛, 李夕海, 韩绍卿等. 位场向下延拓三种迭代方法之比较. 地球物理学进展, 2011, 26(3): 908-915. Zeng X N, Li X H, Han S Q, et al. A comparison of three iteration methods for downward continuation of potential fields. Progress in Geophys. (in Chinese), 2011, 26(3): 908-915.
[27] 陈龙伟, 张辉, 郑志强等. 水下地磁辅助导航中地磁场延拓方法. 中国惯性技术学报, 2007, 15(6): 693-697. Chen L W, Zhang H, Zheng Z Q, et al. Technique of geomagnetic field continuation in underwater geomagnetic aided navigation. Journal of Chinese Inertial Technology, 2007, 15(6): 693-697.
[28] 陈龙伟, 徐世浙, 胡小平等. 位场向下延拓的迭代最小二乘法. 地球物理学进展, 2011, 26(3): 894-901. Chen L W, Xu S Z, Hu X P, et al. The iterative least square method for downward continuation of potential fields, Progress in Geophys. (in Chinese), 2011, 26(3): 894-901.
[29] M. Hanke. Conjugate Gradient Type Methods for Ill-Posed Problems. Longman Scientific & Technical, Essex, UK, 1995.
[30] 程方道. 波数域重磁数据处理中的随机误差分析. 石油物探, 1985, 24(1): 70-76. Cheng F D. Stochastical error analysis for the gravitational and the magnetic data processing in wave-number domain. Geophysical Prospecting for Petroleum, 1985, 24(1): 70-76.
[31] 程方道. 重磁异常波数域变换中误差传递系数的计算及其应用. 矿产与地质, 1987, 1(2): 37-41. Cheng F D. Calculation of error transmission coefficient for gravitational and magnetic anomaly transform in wavenumber domain and its application. Mine and Geology, 1987, 1(2): 37-41.
[32] 熊光楚. 随机误差在频率域的传播. 物探与化探, 1988, 12 (1): 23-27. Xiong G C. The propagation of random errors in frequency domain. Geophysical & Geochemical Exploration, 1988, 12(1): 23-27.
[33] 管志宁, 张昌达等. 磁法勘探重要问题理论分析与应用. 北京:地质出版社,1993. Guan Z N, Zhang C D, et al. Theoretical analysis of important problems in magnetic exploration and its application. Beijing: Geology Press, 1993.
[34] Chai Y. A-E Equation of Potential Field Transformations in the Wavenumber Domain and Its Application. Applied Geophysics, 2009, 6(3): 205-216. 