二维物性界面深度重磁反演的一种非线性迭代解法──样条函数法  

NONLINEAR ITERATIVE INVERSION OF GRAVITY AND MAGNETIC INTERFACES TO TWO DIMENSIONAL ANOMALIES BASED ON B SPLINE FUNCTION

作　　者：王硕儒[1] 于增慧[2] 

机构地区：[1]青岛海洋大学地球科学学院,青岛266003 [2]中国科学院海洋研究所,青岛266071

出　　处：《海洋与湖沼》2000年第3期302-308,共7页Oceanologia Et Limnologia Sinica

基　　金：国家自然科学基金资助项目! 49274212

摘　　要：提出一种重磁异常非线性反演的直接解法，它是基于迭代法求解一个非线性积分方 程。当物性界面起伏分布△ｈ满足｜△ｈ｜＜（２－１）ｈ时，积分的被积函数可展成西△ｈ的幂级数， 考虑△ｈ高次项的影响，应用简单迭代法可将问题化成相应的线性积分方程，并可用Ｂ样条函 数求解。在不同磁化方向的△Ｔ和Ｚａ异常上和变密度界面引起的重力异常上实现了模拟反 演。精度显著高于线性直接解。在南海东沙群岛地区进行了莫霍面的反演计算，效果好。Inversion of gravity and magnetic interfaces by gravity and magnetic anomalies, which demonstrates the distribution of subterranean materials, is an important realm in geophysics. In this realm, the common methods are based on Parker formula in frequency domain. For these methods there is the problem of choosing the power number of the relief of interfaces △h, and to magnetic interfaces they are only applicable under the condition of perpendicular magnetization. In order to avoid the errors caused by frequent Fourier transferring in frequency domain and by omitting the higher power items for linear solutions in space domain, in mis article a new nonlinear iterative solution in space domain is presented. It solves the integral equation based on simple and flexible B spline. When the relief of interfaces ah is in ｜△h＋ < (2-1)h, the effect of the higher power items is considered commonly, whereas it is omitted by other methods. Models and examples show that this method is better than Parker's formula method and the linear solution in space domain.

关 键 词：生磁异常反演 B样条函数 物性界深度 迭代解法 

分 类 号：P631.24[天文地球—地质矿产勘探]