三、磁异常处理转换.ppt

《三、磁异常处理转换.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、1、加深对磁性体磁异常在空间域处理转换原理与作用的认识2、用Matlab语言编程实现球体或水平圆柱体磁异常(包括Za、Ha、Δt)的向上延拓和分量转换。一、实验目的二、实验内容利用两个大小与埋深不同的球体或水平圆柱体产生的磁异常（ΔT、Za、Ha），进行上延计算；利用单个球体或水平圆柱体产生的磁异常（ΔT、Za、Ha），进行分量转换分析磁异常上延计算与分量转换的作用对计算结果进行误差分析，分析计算异常与理论异常误差的主要来源数据准备模型（球体或水平圆柱体）示意图数据准备设在地下有两个磁性体，一个浅而小（半径r=5m，埋深R=

2、25m），一个深而大（半径r=50m，埋深R=100m），有效磁化强度Ms=0.2A/m,is=90°情况下，正演其磁异常（ΔT、Za、Ha）。三、实验要求一、上延计算利用两个球体或水平圆柱体的正演叠加磁异常（Za、Ha或△T）进行上延计算，画出ΔT、Za、Ha各分量向上延拓5m、10m的异常图；比较上延前后各磁异常分量曲线的不同及不同地质体的上延后效果的不同，分析上延计算的作用。对计算结果进行误差分析，分析计算异常与理论异常误差的主要来源三、实验要求二、分量转换利用单个球体或水平圆柱体的正演磁异常（Za、Ha）进行分量转换

3、，画出Za、Ha分量转换图；比较分量转换前后磁异常曲线形态，分析分量转换的作用。对计算结果进行误差分析，分析计算异常与理论异常误差的主要来源四、实验原理1、空间域向上延拓：换算平面位于实测平面之上。主要用途是削弱局部异常干扰，反映深部异常。设坐标原点位于计算点下方实测剖面上，延拓高度为一个点距h，则原点的向上延拓公式：四、实验原理剖面向上延拓计算原理剖面向上延拓计算原理经过复杂的求解：剖面向上延拓计算原理：为延拓高度剖面向上延拓计算原理为之间某值（中值定理）～剖面向上延拓计算原理怎么选择？剖面向上延拓计算原理(1).～等间

4、距划分剖面向上延拓计算原理(1).～等间距划分剖面向上延拓计算原理(1).～等间距划分剖面向上延拓计算原理(1).～等间距划分剖面向上延拓计算原理(2)不等间距划分怎么选择？剖面向上延拓计算原理(2)不等间距划分怎么选择？剖面向上延拓计算原理(2)不等间距划分取：剖面向上延拓计算原理取：则：剖面向上延拓计算原理取：则：P剖面向上延拓计算原理取：则：P量板法(等权)剖面向上延拓计算原理量板法怎么提高精度？P剖面向上延拓计算原理四、实验原理2、原点处的磁异常的分量换算公式：ZaHa：HaZa：磁异常分量换算原理经过求解：

5、拉普拉斯方程的诺伊曼问题（p166)：互算：得：得：互算：得：互算：互算：得：在原平面上换算时,z=0,又设在坐标原点为计算点，则x=0：互算：将积分区间分成5个部分：互算：在小积分区间内线性变化：互算：在小积分区间内线性变化：互算：在小积分区间内线性变化：互算：在小积分区间内线性变化：互算：在积分区间分成一系列等间隔的小区间：互算：在积分区间分成一系列等间隔的小区间：在每个小区间内，磁异常近似线性变化。利用积分中值定理：互算：在积分区间分成一系列等间隔的小区间：在每个小区间内，磁异常近似线性变化。利用积分中值定理：互算：在

6、积分区间没有值，对于二维磁异常，设按照距离平方反比衰减：互算：在积分区间没有值，对于二维磁异常，设按照距离平方反比衰减：互算：其中：互算：此外，可得：五、程序编写与图示结果例子：1）数据准备定义球体或水平圆柱体的参数r=10;%浅水平圆柱体截面半径，单位：mr0=30;%深水平圆柱体截面半径，单位：mR=20;%浅水平圆柱体柱轴的埋深,单位：mR0=50;%深水平圆柱体柱轴的埋深,单位：m正演计算两个（单个）球体或两个（单个）水平圆柱体的磁异常（Za、Ha、△T）五、程序编写与图示结果例子：2）计算上延异常（Za、….)%Z

7、a向上延拓m=2米t=10;%每个点延拓采用的测点数nx=length(x);%X方向测点数dx=1;%X方向测点间距C0=(1/pi)*(atan((1/2)*m)-atan((-1/2)*m))五、程序编写与图示结果例子：fori=1:length(x)-t*2forj=1:tC(j)=(1/pi)*(atan((j+1/2)*m)-atan((j-1/2)*m));Z(j)=C(j)*(Za(i+t-j)+Za(i+j+t));endZmh(i)=C(j)*Za(i+t)+sum(Z);endx1=xmin+t*dx:

8、dx:(xmin+(nx-1)*dx-t*dx);3）绘图六、结果分析1、球体或水平圆柱体磁异常任一分量上延前后的变化特征2、球体或水平圆柱体磁异常任意两个分量转换前后的变化特征3、总结磁异常上延与分量转换的作用4、对计算结果进行误差分析，分析计算异常与理论异常误差的主要来源七、实验报告报