资源描述:
《应用地磁学课6new》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、第六章磁异常的处理与转换第六章磁异常的处理与转换目的:①满足解释理论要求的假设条件•根据磁异常的数学物理特征复杂→简化(曲面→平面;迭加•进行必要的数据加工处理→孤立;等)•满足某些特定需要②满足解释方法(某分量→另一分量;磁化方向换算;磁场值→频谱值)。磁异常的处理与转换是磁力勘探理论的③突出磁异常某一方面特点重要组成部分,对提高解释效果很重要(上延→压制浅部、突出深度滤波→可突出浅或深,或某一方面)第六章磁异常的处理与转换磁异常处理转换及其意义主要内容:①圆滑、划分异常1、向上延拓Î压制浅层(干扰),突出深部(趋势);②空间换算2、水平导数Î突出方向构造信息;③分量换算3、垂直导数④导数
2、换算Î突出浅层场源信息;4、化磁极⑤不同磁化方向换算Î消除斜磁化影响,简化磁场形态;10001000800800水平导数Î突出方向构造信息6006001000100040040080080020020060060000向上020040060080010000200400600800100010001000Î延拓4004008008002002006006004004000002004006008001000020040060080010002002000002004006008001000020040060080010001化磁极Î消除斜磁化影响垂直导数Î突出浅层场源信息1000100080
3、0800600600ÎÎ400400200200000200400600800100002004006008001000地磁场化磁极Î消除斜磁化影响化磁极Î消除斜磁化影响10001000800800Î6006004004002002000002004006008001000020040060080010001000I=60度I=30度10001000800800800600600600400400400Î200200200000020040060080010000200400600800100002004006008001000I=-30度I=-60度I=90度第六章磁异常的处理与转换第六章
4、磁异常的处理与转换处理与转换要遵循的原则:§1、空间域剖面磁异常处理与转换①选择合适的处理转换方法。要求使用者必须掌握各种处理与转换方法§2、频率域磁异常处理与转换的原理。②处理转换只是一种数学加工方法,能突出某些信息,而不能获得数据中不包含的信息。2§1、剖面处理、转换§1、剖面处理、转换一、磁异常空间换算一、磁异常空间换算即解析延拓(即解析延拓)根据观测面上的实测异常,换算场源以外①向上其他空间位置的磁异常①向上②向下②向下§1、剖面处理、转换§1、剖面处理、转换一、解析延拓原理一、解析延拓原理•磁异常为位函数,磁性体外部空间的磁•把观测面以上的半空间作为调和域D,在域内磁场和磁位均具有
5、调和函数的性质。场随着观测面距离的增加而衰减,到无穷远趋于•调和函数的性质:零,且域内满足拉普拉斯方程。在调和域中满足拉普拉斯方程且具有连•磁测数据通常是在地表或近地表取得,如果测区有足够的范围,在观测面上已有实测的磁异常续的一、二阶导数;及法向导数,则可求出观测面上半空间中的磁异调和域中函数值可以通过域边界上的函常,即实现了磁异常的解析延拓。数值及法向导数值来计算。§1、剖面处理、转换§1、剖面处理、转换1、向上延拓:1、向上延拓:Ux(,0)Ux(,0)xxrP(,0)xrP(,0)xvvzzQ(,)ξζQ(,)ξζ3§1、剖面处理、转换§1、剖面处理、转换⎧2⎧21、向上延拓:⎪∇=<
6、Uz0(0)1、向上延拓:⎪∇Uz=<0(0)⎪⎪⎨⎨⎪UU=(,)xo⎪UU=(,)xo⎪z=0⎪z=0⎩⎩经过复杂的求解:Ux(,0)zUd+∞(,0)ξξUxz(,)=−πξ∫−∞()−+x22zxrP(,0)xvzTd+∞Δ(,0)ξξΔ=Txz(,)−πξ∫−∞()−+x22zzQ(,)ξζ§1、剖面处理、转换§1、剖面处理、转换1、向上延拓:zUd+∞(,0)ξξzUd+∞(,0)ξξUxz(,)=−Uxz(,)=−πξ∫−∞()−x22+zπξ∫−∞()−+x22z积分中值定理设坐标原点位于计算点下方实测剖面,延拓高度为一个点距(h):∞1ξi+1h•U(ξ,0)dξUiξi+
7、1hU(0,−h)==dξ∞∑∫ξ22∑∫ξ221(ξhU•ξ,0)dξi=−∞πiξ+hiπiξ+hi+1Uh(0,−=)∑∫22i=−∞πξξi+hh:为延拓高度Ui为U(ξi,0)~U(,ξi+10)之间某值(中值定理)§1、剖面处理、转换§1、剖面处理、转换Uiξi+1h(−11等间距划分Uh(0,−=)∑∫22dξ(1).i)h~(i+)hiπξξi+h22P(0,−h)P(0,−h)hx11xP(