一种二度磁异常分量换算的数字计算方法.pdfVIP

一种二度磁异常分量换算的 数字计 算方法 许 明 七 吉林省地质局第六地质调查所 本文从二度磁异常水平分量与垂直分量之间的基本关系式——希尔伯特积分出发， 应用数字积分的梯形公式，提出了由二度磁异常垂直分量值换算水平分量的新方法。其 换算公式为： HaCKJ=21- (ZaCK+lJ-ZaCK-1J)+-1- ．圣1 mCiJi-K “ ： l—K 1≠K 偿 (i=l或n) 式中mCiJ=) ll (i≠1或n) Ha[K]和ZaCiJ分别是第K个观测点上的磁异常水平分量值和第i个观测点上的磁异 常垂直分量观测值；n为磁异常剖面上观测点的个数。 本文所推荐的方法具有边缘损失小、计算精度高和便于在电子计算器上实现数字计 算等优点。文中还给出了水平圆柱体和有限延深薄板状体理论模型上的算例及换算结果 的反演数据，以验证方法的可靠性。 众所周知，许多磁异常反演计算，尤其是向量交绘法、奇点法等，不仅要用磁异常 垂直分量，还要用水平分量。但在实测资料中，常常只有垂直分量资料而缺乏水平分量 资料，因此，人们已提高许多由实测磁异常垂直分量换算水平分量的方法[1]。我们将空 间磁异常分量换算方法分为以下两种：一种为等间距取样数字计算法 如（文献CIJ中盼 第一种方法），另一种为不等长度区间内取平均值的量板法。无论是那一种换算方法， 均存在共同的缺点，那就是边部丢点较多。另外，第一种方法，因引进余项影响而求得 的末项系数太大，故当末项Za值取自高异常区时就出现很大误差}第二种换算方法则 在电子计算机 器（）上不宜实现数字计算。 为了克服磁异常分量换算中的上述缺点，我们从二度磁异常水平分量和垂直分量间 的基本关系式出发，提出新的换算方法。 基本原理 当观测面水平的时候，在垂直于二度磁性地质体走向的观测剖面 (X轴沿该剖面 上，磁异常垂直分量Za与水平分量Ha间的关系，用希尔伯特积分表示： Ha (Xk)：j— ocZaC 冗 一。cX一萋三【dx…………… (1) 式 (1)中，Xk为观测剖面上的某一点座标值 ． 在所能达到的观测精度前题下，我们假设观测剖面两端已进入正常磁场区域，则可 用实测区间CXi，XnJ来代替 (1)式中的无穷积分区域。于是有： 一 ／Xn ZaC Ha (Xk)=二L ~jdx… ．(2) 丌 XlX-: 设在观测剖面上的测点数为n，则这些观测点将观测区段分成n-l个小段。 因此以 这些小段分段进行积分 (2)式得出： Ha (Xk) :L n∑l Xi+iZa(鼍 dx…，． (3) “j=1 Xi X-Xl 显然，当i=K时，即当X=Xk时， (3)式中的两项积分是奇性的。故我们将 ~(3)式改写成： Ha似 … ÷ 蔗 I 鬻 ．dx+一} 譬 蔗 +1 Xkdx +jL n-l Xi+i2 蛩襞 dx．… … … “(4) 冗 i=K+l Xi 一 在 (4)式中的后两个积分之和，应用数值积分的梯形公式并令 Xi+i-Xi=l （即 令相邻观测点间距相等且以点距作为量距单位）得： 专 害 ： +1裂 蹇 dx一 去 [惫 曼 。 +专掣 嚣 ]+ _L K《 Za(Xi) … “(5， 兀 1：：