共查询到20条相似文献，搜索用时 703 毫秒
1.
根据位场理论,将观测曲面上测得的△T(或Za)视为位函数,代入偶层位公式的极限形式,求解出对应于△T(或Za)的等效偶层的磁化强度J(α,β,γ)。将J(α,β,γ,)代入偶层场强分量表达式 (α)/(α)n(1/r)dS 只要给定了方向u,就可得到△T沿该方向的一阶导数。将J(α,β,γ)代入(α)T/(α)x,(α)T/(α)y,(α)T/(α)z的积分表达式就能得到△T的二阶导数。因为使用的是二重积分形式的导数的精确表达式,它克服了空间域差商求导的致命弱点。特别是该方法能将起伏地形曲面上的场值向上延拓(包括曲化平)和求一阶(或二阶)导数两个转换过程集中在一个积分表达式中,一次计算完毕。 相似文献
2.
导数换算作为常用的位场数据处理手段具有重要的物理意义。由于实际测量误差和噪声的影响,波数域导数换算存在计算过程不稳定、求导结果精度差的问题。为了降低噪声干扰,提高导数计算精度,提出了在波数域常规导数算子基础上附加Chebyshev低通滤波器的任意阶导数换算方法。通过分析该滤波器的滤波特性,并结合位场异常径向平均功率谱曲线特征,确定了Chebyshev低通滤波法的滤波参数,减少了人为因素对求导结果的干扰。二度体和三度体模型试验表明,与其他3种方法(常规FFT(fast Fourier transformation)法、向上延拓法和ISVD(integrated second vertical derivative)法)相比,该方法计算的导数与理论值的均方根误差最小,求导结果精度较高。在虎林盆地布格重力异常数据处理中,利用该方法计算的垂向二阶导数受噪声干扰小,结果可靠性较高,依此划分出的8条大断裂和11条小断裂在以往研究中均得到证实。 相似文献
3.
在磁力、重力异常的解释中,常常要进行异常场的换算,例如向上延拓、向下延拓,曲化平(将起伏地形上的观测值换算到某一水平面上的场值),磁异常不同分量之间的换算,求变换磁化方向后的异常(包括计算异常的化到磁极),以及异常的一次、二次垂向导数等等。这些换算已有一些方法可以进行,在上延、不同分量的换算和化到磁极等方面均取得了好的效果,但是在许多时候,特别对于观测面起伏高差较大的情况,常常效果不佳。而且有的方法丢掉边部测点的信息较多。在这里我们介绍重磁异常换算的等效源法,这种方法可以将上述不同方面场的换算问 相似文献
4.
垂向导数广泛应用于位场数据处理和解释中,常规的空间域求导公式多以台劳级数为基础,这不可避免会引入截断误差,且不便于求取高阶导数。另一种应用较多的波数域算法。其求导算子在高波数成份产生振荡失真,很不稳定,在实际应用中很难取得理想结果。这里提出了一种计算位场数据垂向导数的稳定算法,在波数域求垂向积分,在空间域计算二个二阶水平方向导数,借助拉普拉斯方程,求取位场异常的各阶垂向导数。算法只在求垂向积分时引入了一次付里叶变换,其它计算全部在空间域中进行,从而有效地避免了吉布斯效应的重复累加传递,有效降低了常规波数域求导算子在高波数成份产生的振荡失真。通过理论模型对比,证明了这里所提出的算法较常规的波数域方法具有良好的稳定性,特别是在计算高阶垂向导数时有明显优势,能有效地提高与垂向导数有关的位场定量解释方法的准确性。 相似文献
5.
《物探化探计算技术》1991,(2)
CMD软件包采用模块结构,用FORTRAN77语言编写,软件包的主要功能:①高精度频率域求导;②不同精度的空间域求异;③三次样条求导;④位场横向数据的垂向导数;⑤弱异常信息的增强。为完成上述功能,模块中有快速富里埃变换、希尔伯特变换、补偿圆偿、绘曲线等模块。软件包设有HELP文件,以帮助用户正确地使用。计算结果可用打印机打出,也可用屏幕显示。更值得一提的是,该软件提供绘曲线程序,计算所得的导数值可用导函数的曲线形式在屏幕上显示,存绘图仪上绘制导数曲线图,也可以在打字机上输出导数曲线图。不仅如此,他是按照国家标准局(现国家技术监督局标准化司)已经批准的两项计算机软件国家标准: 相似文献
6.
《物探化探计算技术》1982,(4)
我们已知磁力线的分布及磁位法,确定磁化倾角。利用各方向的水平分量或水平导数零线的交点,确定磁性体的中心。利用各方向的水平导数极值点的联线,确定磁性体的边界范围。现以谢尔塔拉异常为例说明其实际效果。谢尔塔拉异常(图1)从异常形态及正负比来看,颇具近垂直磁化的特征。但这与本区67。的正常地磁倾角是不相符的。为此,对异常进行了不同方法的计算,以进一步探讨其磁化倾角,这对进行正演及各种计算是十分必要的。 相似文献
7.
在位场数据处理中,垂向导数具有重要的物理意义。其在一定程度上可以划分不同深度和大小异常源产生的叠加异常,且导数的阶次越高,这种分辨能力就越强,但通常认为高阶导数的换算是不稳定的。本文在Tikhonov正则化求位场垂向高阶导数的基础上,结合迭代法进行逐次逼近,提出了位场高阶导数的Tikhonov正则化迭代法,并且得到Tikhonov正则化迭代法的递推公式。通过对该方法的滤波特性分析可以看出,该方法计算的位场垂向高阶导数具有一定的稳定性及保幅性。模型试验和实际数据的处理表明,该方法计算结果较常规FFT求导法有更高的稳定性和实用价值。 相似文献
8.
9.
10.
在物化探电子计算技术79年第2期中,我们讨论了用等效源法进行重磁异常换算的效果[1]。方法的基本思想是,采用最优化方法求出观测异常的场源分布,根据这个场源分布,再来计算任意点的各种场值,如向上延拓,向下延拓,曲化平(将起伏地形上的观测值换算到某一水平面上的场值),磁异常不同分量之间的换 相似文献
11.
在野外1∶1万地质测量、井下大量坑道观察和编录基础上,对冀东峪耳崖金矿床矿区1∶1万高精度磁法测量数据进行化极、延拓、求导处理,将峪耳崖矿区的高磁异常划分为低缓正异常场、低缓负异常场和中高正异常场3种场区,解译出不同方向线性断裂构造31条,以北东向、北西向构造为主,构成矿区构造格架;磁化极异常进行4个高度50 m、100 m、200 m、500 m垂向二阶导数计算,形态受北东向、北西向及南北向构造综合控制,沿北东方向带状分布。结合已知矿体的构造控矿规律及岩体控矿特征,确定区内成矿远景区圈定依据为成矿磁性体形态、产状变化地段;北东向构造面与北西向、南北向、东西向构造交汇复合部位;黄铁矿化、硅化、大理岩化及绢云母化等蚀变发育部位。 相似文献
12.
曲面位场导数的换算方法大多基于调和函数的积分方程数值解。本文提出了基于B样条函数的曲面上位场导数直接计算法。引入按弧长S作变量的三次样条插值,得出了位场水平导数的计算式。该方法还有一个特点,那就是可以用一元样条来实现曲面上的求导。这些,都决定了该方法的简便性和精确性。利用曲线上和曲面上的两种理论模型进行试算,分别计算了Zα/x和Δg/Z~2计算和理论的结果对比表明,符合很好。误差主要取决于曲面起伏大小,异常变化陡度和取样的间隔。 相似文献
13.
位场垂向二次导数是基于拉氏方程求得的。1969年B.K.Bhattacharyya实现了二元三次样条函数计算位场导数的方法。本文采用了非等距结点B样条函数求导的方法,它具有灵活统一的优点。取作为二次微商计算式。由于用曲线插值求导取代曲面插值求导,程序更便于编写。在计算机上分别用三种模型作了试算,效果很好,精度比其它方法(包括上述二元三次样条函数法)提高了一个级次。 相似文献
14.
引言众所周知,重力垂直二次导数是布伽异常沿垂直方向的二次导数,也就是重力位沿垂直方向的三次导数,是划分区域场和局部异常的有效方法.重力垂直二次导数有许多计算公式,所得结果不同,原因何在,究竟它们的意义是什么?它们之间有什么异同?能否提出一个准则来选择某一公式?计算精度如何评定?都是值得我们探讨的.本文结合一些实际成果粗略谈几点认识,抛砖引玉,请同志们批评指正.一二次导数法的一般意义大家知道,布伽异常中包含着比局部异常大得多的区域性异常,而和矿产分布有关的局部异常,往往被其掩盖或发生畸变.过去,我们习惯于从实测布伽异常图上来分辨局部异常,就是通常所讲的封闭和凸扭.这种 相似文献
15.
利用重磁水平和垂直二阶导数确定东北地区梯度带 总被引：1，自引：0，他引：1
分析了东北地区的重、磁场特征,同时对研究区的布格重力异常和航磁异常数据进行水平和垂直求导。水平和垂直参量图能给出断裂和走向的信息。结合布格重力异常图可以得出东北地区重力梯度带的走向和边界,并对比重力和磁力的垂直和各个角度的水平二阶导数来综合确定东北地区梯度带。利用求重力和磁力的水平和垂直二阶导数来划分东北地区梯度带。 相似文献
16.
位场数据异常存在的线性特征往往对应着地下断裂构造、不同岩性地质体的边界接触带或者其他具有一定密度或磁性差异的构造特征。对这些线性特征进行增强、提取并进行半定量地解释是重磁资料处理的主要内容。然而,位场数据中若混入噪声干扰,再利用总水平导数、总水平梯度倾斜角等方法进行高次求导运算会把噪声放大,导致提取的线性构造位置发生偏离甚至出现错误。为了获取较为准确的线性构造及边界位置,选取对噪声干扰不敏感的小波模极大值方法,将该方法用于模型试验和实际矿区数据处理中,都能较好地定位出异常体模型边界和矿体的投影边界,表明小波模极大值是一种有效的重磁异常线性特征增强与提取方法。 相似文献
17.
航空电磁数据中经常存在条带状误差,为消除这种误差,通常采用窗口滤波。但是,窗口参数及滤波参数要求极为严格,且不能区分调平误差和其类似大小的地质异常,需要操作者多次尝试以寻求合适的处理参数。这里提出了导数突变点统计法,该方法首先在沿伪切割线方向（垂直飞行方向）,对网格化的航空电磁数据求取伪切割线的导数,采用高通滤波方法滤掉幅值比较小的导数突变点;然后统计每条航线上的突变点个数,根据预设的阈值,确定出存在条带误差的航线;再针对这条航线上的数据,采用均值滤波得到调平结果。导数突变点统计法能自动找出存在误差的航线,并只针对误差线进行滤波,能有效消除条带状误差,极大地保留真实地质异常。导数突变点统计方法算法简单,输入参数少,数据处理速度快,为航空电磁数据调平提供了一种新思路。 相似文献
18.
用于频率域航电数据处理的伪切割线自动调平法 总被引：1，自引：0，他引：1
针对频率域航空电磁数据中存在的区块状测量误差,提出了伪切割线自动调平法。首先利用一维移动滤波滤掉高频地质异常,得到包含背景值和调平误差的中间值;通过对伪切割线求导,由突变点区间统计法判断出的调平边界突变点位置来划分待调平区域;沿伪切割线方向由边界向外延拓一定距离,确定伪切割线范围,在此范围内,对边界突变点采用中值滤波;然后在确定的伪切割线范围内积分得到背景值,背景值与中间值作差求得调平误差;最后与原始数据结合得到调平结果。该方法输入参数少,能自动划分调平区域,极大地减少了工作强度及人为误差,同时能有效地消除区块状误差,保留真实地质异常。 相似文献
19.
基于DCT的重力归一化总梯度及相位在识别断裂构造中的研究与应用 总被引：1，自引：0，他引：1
区别于传统的别列兹金(V.M.Berezkin)法(即傅立叶级数法)和傅立叶变换法,提出在位场转换时采用离散余弦变换(DCT)计算重力归一化总梯度和归一化相位,并把二者结合起来研究断裂构造.分别研究了比较典型的无限延伸铅垂台阶、有限延伸铅垂台阶和正断层模型的GH场和相位特征,模型实验证实了方法的可行性和有效性.得出了寻找地质体断裂的空间位置和水平尺度的规律:在断裂处GH场有极大值圈闭,相位发生突变;多个断面相位有多个突变,两个突变点之间的水平距离反映地质体的水平尺度;不同深度相位突变点之间的连线大致反映了断裂的倾向.将研究方法应用于兴凯-虎林线剖面异常数据处理中,根据该线GH场和相位特征,划分出13条断裂. 相似文献
20.
为调查庐枞矿集区中东部七家山—马口地区深部矿产资源潜力,利用小波变换方法提取深源重磁场,通过对垂向一阶导数Z与水平方向导数模S进行参数变换,突出深源重磁场的局部异常; 结合频谱激电(spectrum induced polarization,SIP)、可控源音频大地电磁测深(controlled source audio-frequency magneto tellurics,CSAMT)成果,在调查区深部查明了一条长大于3.0 km、宽大于0.5 km的高重低磁、低阻高极化异常带; 根据调查区内外已知铜矿(化)体分布规律,结合已知矿(化)体在深源重磁异常中的分布及钻孔(zk4601)钻探验证结果,认为上述异常带是铜多金属矿化蚀变带的反映,是寻找铜多金属矿的有利部位,具有较好的找矿前景。 相似文献