地球物理学进展 ›› 2012, Vol. ›› Issue (4): 1509-1518.doi: 10.6038/j.issn.1004-2903.2012.04.026
陈龙伟, 胡小平, 吴美平, 马涛
收稿日期:
2011-08-10
修回日期:
2011-11-12
出版日期:
2012-08-20
发布日期:
2012-08-20
作者简介:
陈龙伟,男,1982年生,国防科技大学博士生,主要从事位场数据处理.(E-mail: longweichen_nudt@163.com)
基金资助:
国家自然科学基金项目(61174206)资助.
CHEN Long-wei, HU Xiao-ping, WU Mei-ping, MA Tao
Received:
2011-08-10
Revised:
2011-11-12
Online:
2012-08-20
Published:
2012-08-20
摘要: 在空间域进行位场延拓,需要数值求解第一类Fredholm积分方程,由于所得方程组系数矩阵不是稀疏矩阵,求解该方程组需要的计算机内存大,计算量大,导致延拓算法在一般计算机上难以实现,阻碍了对空间域位场延拓方法的研究.在分析系数矩阵结构特征的基础上,本文证明了方程组系数矩阵是对称的分块Toeplitz型矩阵.利用系数矩阵的对称性和分块Toeplitz型矩阵与向量相乘的快速算法,解决了系数矩阵的存储和计算问题,使得空间域位场延拓成为可能,为研究新的位场延拓方法和分析延拓误差提供了一条新的途径.利用模型数据和实测资料,对空间域位场向上延拓、空间域积分迭代法向下延拓进行了检验,结果证实了空间域位场延拓的可行性和正确性.
中图分类号:
陈龙伟, 胡小平, 吴美平, 马涛. 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 .
[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. |
[1] | 石颖,井洪亮,李莹. 反馈迭代法压制表面多次波效果分析[J]. 地球物理学进展, 2012, (4): 1493-1500. |
[2] | 宋维琪,杨晓东. 基于射线追踪的微地震多波场正演模拟[J]. 地球物理学进展, 2012, (4): 1501-1508. |
[3] | 张小莉,冯乔,冯强汉,陈凤喜,刘晶晶,孟展,刘广伟. 鄂尔多斯盆地中东部长6复杂低渗储层测井评价[J]. 地球物理学进展, 2012, (4): 1519-1524. |
[4] | 佘德平. 多震源地震正演数值模拟技术[J]. 地球物理学进展, 2012, (4): 1533-1540. |
[5] | 井涌泉,余杰,范廷恩,陈开远,胡光义. 秦皇岛32-6油田明化镇组河流相储层地球物理参数影响因素研究[J]. 地球物理学进展, 2012, (4): 1541-1547. |
[6] | 许璐,孟小红,刘国峰. 逆时偏移去噪方法研究进展[J]. 地球物理学进展, 2012, (4): 1548-1556. |
[7] | 朱海龙,张山. 高保真VSP波场分离方法研究[J]. 地球物理学进展, 2012, (4): 1557-1564. |
[8] | 杨晓利,王童奎,柳晓风,郑晓凤,刑文军. 相控地震多属性预测技术在南堡凹陷碳酸盐岩储层的应用[J]. 地球物理学进展, 2012, (4): 1565-1571. |
[9] | 刘淑芬,文瑞霞,陈昌旭,王国雷,张冬杰,梁清. 层控垂向多参数拟合变速成图方法研究——以明水绥化地区为例[J]. 地球物理学进展, 2012, (4): 1572-1580. |
[10] | 于正军,董冬冬,宋维琪,高艳珂,吴时国. 相带控制下协克里金方法孔隙度预测[J]. 地球物理学进展, 2012, (4): 1581-1587. |
[11] | 刘军迎,雍学善,杨午阳,高建虎,周春雷. 基于叠前方位振幅的大港——埕海地区奥陶系风化壳裂缝储层的叠前预测[J]. 地球物理学进展, 2012, (4): 1588-1597. |
[12] | 陈永波,潘建国,高建虎,白晓寅,尹路,许多年,陈雪珍,王斌. 岩性油气藏关键技术攻关与应用研究——以准噶尔盆地准东阜11井区为例[J]. 地球物理学进展, 2012, (4): 1598-1608. |
[13] | 石磊,郭良辉,孟小红. 磁总场异常垂直梯度三维相关成像[J]. 地球物理学进展, 2012, (4): 1609-1614. |
[14] | 王佳文,柴冬梅,凡友华,刘忠诚. 广义S变换在瑞利波频散曲线提取中的应用[J]. 地球物理学进展, 2012, (4): 1615-1623. |
[15] | 郝建飞,周灿灿,李霞,程相志,李潮流,宋连腾. 页岩气地球物理测井评价综述[J]. 地球物理学进展, 2012, (4): 1624-1632. |
阅读次数 | ||||||
全文 |
|
|||||
摘要 |
|
|||||