共查询到20条相似文献，搜索用时 593 毫秒
1.
位场垂向二次导数是基于拉氏方程求得的。1969年B.K.Bhattacharyya实现了二元三次样条函数计算位场导数的方法。本文采用了非等距结点B样条函数求导的方法,它具有灵活统一的优点。取作为二次微商计算式。由于用曲线插值求导取代曲面插值求导,程序更便于编写。在计算机上分别用三种模型作了试算,效果很好,精度比其它方法(包括上述二元三次样条函数法)提高了一个级次。 相似文献
2.
根据位场理论,将观测曲面上测得的△T(或Za)视为位函数,代入偶层位公式的极限形式,求解出对应于△T(或Za)的等效偶层的磁化强度J(α,β,γ)。将J(α,β,γ,)代入偶层场强分量表达式 (α)/(α)n(1/r)dS 只要给定了方向u,就可得到△T沿该方向的一阶导数。将J(α,β,γ)代入(α)T/(α)x,(α)T/(α)y,(α)T/(α)z的积分表达式就能得到△T的二阶导数。因为使用的是二重积分形式的导数的精确表达式,它克服了空间域差商求导的致命弱点。特别是该方法能将起伏地形曲面上的场值向上延拓(包括曲化平)和求一阶(或二阶)导数两个转换过程集中在一个积分表达式中,一次计算完毕。 相似文献
3.
起伏地形位场快速延拓新方法 总被引：2，自引：0，他引：2
将起伏地形上的测点垂直投影至水平面上,并将地形上的位场u放在水平面的对应点上,这样就有一个平面上的假想位场u。用快速傅立叶变换（FFT）计算平面位场的垂向导数δu/δn,将这个垂向导数移至地形的对应点上,近似作为地形上的位的垂向导数δu/δn。根据Green公式,从地形上的位u和位的法向导数δu/δn,用积分的方法可计算地形上部任一点的位。在计算出地形上部一个平面上的位后,可用FFT迅速计算其它平面的位,包括与地形相交的平面上的位。在通常的位场延拓方法中,如用等效源法解一个线形代数方程组要需用很大内存量,要耗去大量计算时间。本方法的主要特点是省略了解线形代数方程组的步骤,从而大大节约计算时间。模型计算表明,本方法的计算精度是令人满意的。用本方法对航磁实际资料进行延拓,取得较好的效果。 相似文献
4.
在位场数据处理中,垂向导数具有重要的物理意义。其在一定程度上可以划分不同深度和大小异常源产生的叠加异常,且导数的阶次越高,这种分辨能力就越强,但通常认为高阶导数的换算是不稳定的。本文在Tikhonov正则化求位场垂向高阶导数的基础上,结合迭代法进行逐次逼近,提出了位场高阶导数的Tikhonov正则化迭代法,并且得到Tikhonov正则化迭代法的递推公式。通过对该方法的滤波特性分析可以看出,该方法计算的位场垂向高阶导数具有一定的稳定性及保幅性。模型试验和实际数据的处理表明,该方法计算结果较常规FFT求导法有更高的稳定性和实用价值。 相似文献
5.
本文从单层位出发,导出了曲面上位场转换所需的基本公式:单层场强分量及单层场强分量的一阶与二阶导数表达式。再从联系磁位和引力位的泊松公式出发,导出了“磁化方向与磁场分量方向互换定理“,从而建立了一个多功能的曲面上的位场转换解释系统,给出了使用该系统处理与解释实际资料的例子。 相似文献
6.
7.
在双正交小波分析下,实现了重磁异常数据高精度的导数计算。该算法是在小波分析得到的系数基础上,采用Monte-Carlo方法进行噪声的估计,根据获得参数分别进行去噪处理。理论计算结果表明:算法的计算效率高,计算的结果精度高,抗干扰能力强。实际数据处理结果表明:采用基于小波变换的方法完成导数处理,结果对地质填图中地质体边界的划分和磁性地层参数分析提供了重要的资料,也为采用基于导数的位场数据反演解释方法提供了优良的基础数据。 相似文献
8.
位场场源边界识别是位场数据分析解释中的重点内容之一。相邻异常源由于场的叠加而难以区分,尤其是对于叠加异常中不同深度的场源边界的识别。倾斜角法是位场数据处理中常用的边界增强方法之一,能有效突出低振幅、短波长的异常特征。本文基于倾斜角法定义了新的位场边界探测方法,该方法利用倾斜角方程中位场二阶垂直导数与垂向导数水平导数之比进行计算,能有效区分相邻异常源产生的叠加异常,突出叠加异常中不同深度场源的边界特征。通过综合模型和实测数据,并与现有的基于导数的位场边界探测方法进行比较,评价了该方法的应用效果。结果表明,该方法对场源深度的变化不太敏感,它可以更准确地确定场源边界,区分相邻异常源的叠加异常,尤其对密集型组合地质构造的识别有较好的效果,即使是在场源埋深较大情况下,也能突出叠加异常中不太明显的微弱变化的线性异常特征,以及弱异常源的形态特征,可为精细研究局部地质构造提供技术支持。 相似文献
9.
位场垂向导数广泛应用于位场数据处理与解释当中,目前求取位场各阶垂向导数方法大致可以分为两类：①在频率域进行;②在空间域进行。针对这两类方法存在着稳定性较差的问题,这里介绍了计算位场各阶垂向导数的ISVD算法,并通过对单一地质体模型和地质体组合模型产生的布格重力异常进行处理分析,证明了ISVD算法在计算位场各阶垂向导数时,特别是在计算位场高阶垂向导数时,比在频率域及空间域计算结果具有较好的稳定性,ISVD算法能够有效地识别出浅层小尺度地质体边界。为了验证该方法对实际资料的处理效果,对某区块实际布格重力异常数据应用ISVD算法试算各阶垂向导数,结合地质背景等资料划分出断裂体系。结果证明,应用ISVD算法计算的位场各阶垂向导数,对于确定浅层小尺度地质体边界有较好的效果。 相似文献
10.
导数换算作为常用的位场数据处理手段具有重要的物理意义。由于实际测量误差和噪声的影响,波数域导数换算存在计算过程不稳定、求导结果精度差的问题。为了降低噪声干扰,提高导数计算精度,提出了在波数域常规导数算子基础上附加Chebyshev低通滤波器的任意阶导数换算方法。通过分析该滤波器的滤波特性,并结合位场异常径向平均功率谱曲线特征,确定了Chebyshev低通滤波法的滤波参数,减少了人为因素对求导结果的干扰。二度体和三度体模型试验表明,与其他3种方法(常规FFT(fast Fourier transformation)法、向上延拓法和ISVD(integrated second vertical derivative)法)相比,该方法计算的导数与理论值的均方根误差最小,求导结果精度较高。在虎林盆地布格重力异常数据处理中,利用该方法计算的垂向二阶导数受噪声干扰小,结果可靠性较高,依此划分出的8条大断裂和11条小断裂在以往研究中均得到证实。 相似文献
11.
在地球物理参数反演中，观测的物理场量，例如波场中的位移或电流场中的电位，对于介质参数的偏导数的计算是非常重要的。本文采用自伴随算子和格林函数的方法给出了弹性波场和稳恒电流场问题的Ｆｒｅｃｈｅｔ导数解析式。该式表明，Ｆｒｅｃｈｅｔ导数可由格林函数Ｇ及其梯度△Ｇ求得。把任意场源正演问题的解和Ｆｅｅｃｈｅｔ导数的计算统一起来，并提出用有限元方法计算任意背景场下的格林函数，从而得到Ｆｒｅｃｈｅｔ导数的数值解。文中以稳恒电流场问题为例，给出了直流电法视电阻率偏导数的数值结果。 相似文献
12.
在地球物理参数反演中，观测的物理场量，例如波场中的位移或电流场中的电位，对于介质参数的编数的计算是非常重要的。本文采用自伴随算子和格林函数的方法给出了弹性波场和稳恒电流场问题的Ｆｒｅｃｈｅｔ导数解析式。该式表明，Ｆｒｅｃｈｅｔ导数可由格林函数Ｇ及其梯度ΔＧ求得。把任意场源正演问题的解和Ｆｒｅｃｈｅｔ导数的计算统一起来，并提出用有限元方法计算任意背景场下的格林函数，从而得到Ｆｒｅｃｈｅｔ导数的数值解 相似文献
13.
导数计算是位场（重磁）数据处理中必不可少的技术手段,现今大多采用Fourier变换来进行。Hartley变换是在Fourier变换基础上定义的一种实数域运算,比Fourier变换更加对称,所需的运算量更少。笔者推导出基于Hartley变换的重磁异常导数计算公式,计算结果与理论值之间误差小于5%,通过理论模型证明Hartley变换可替代Fourier变换进行位场异常的导数计算,且受噪音干扰较小。将Hartley变换用于位场边界识别滤波器的计算,获得了清晰的断裂分布。 相似文献
14.
15.
本文用样条函数法实现了二维位场的向上延拓和向下延拓,向上延拓用样条函数法计算泊松积分,向下延拓用样条函数法求解第一类Ｆｒｅｄｈｏｌｍ积分方程。模型和实例计算表明,本文方法正确可行且精度高。 相似文献
16.
本文用样条函数法实现了二维位场的向上延拓和向下延拓，向上延拓用样条函数法计算泊松积分，向下延拓用样条函数法求解第一类Ｆｒｅｄｈｏｌｍｔｍ积分方程。模型和实例计算表明，本文方法正确可行且精度高。 相似文献
17.
针对向下延拓的不适定问题,提出了一种新的位场向下延拓方法-导数迭代法。从位场垂向一阶导数的定义出发,将观测面和向上延拓、向下延拓同等高度平面上的位场值近似联系起来,采用迭代法进行逐次逼近,推导出空间域位场向下延拓的导数迭代法的递推公式。考虑到空间域向上延拓积分方程实现的复杂性,对递推公式进行快速傅里叶变换,整理得到波数域中的导数迭代公式,同时证明了该方法的收敛性。模型检验和实例分析均表明迭代法相对直接FFT法具有稳定性强和下延深度大的优点。 相似文献
18.
位场数据边缘增强对于地质构造研究和矿床勘查具有重要意义。现有的边缘增强方法多数基于位场数据的导数构建,往往存在对噪声敏感的缺点,且处理结果中常有虚假边界的出现。本文基于数学形态学良好的形状分析和噪声压制性能,通过数学形态学基本算子的组合叠加,提出了一种新的位场数据边缘增强方法。该方法容易实现,且不需要计算位场数据的水平和垂向导数。板状体模型试验结果表明,该方法能有效地对重力异常和化极后的磁异常数据进行边缘增强。针对三维模型试验中正负异常同时存在的情况,相比以往几种方法,该方法输出的边缘更清晰且没有额外边缘的产生,并表现出更强的抗噪能力。最后将该方法和以往几种方法应用到老挝万象地区实际重力异常数据的处理中,结果显示该方法刻画出的边缘更加连续清晰;同时结合以往地质资料,推断了成矿的有利部位。 相似文献
19.
《物探化探计算技术》1991,(2)
CMD软件包采用模块结构,用FORTRAN77语言编写,软件包的主要功能:①高精度频率域求导;②不同精度的空间域求异;③三次样条求导;④位场横向数据的垂向导数;⑤弱异常信息的增强。为完成上述功能,模块中有快速富里埃变换、希尔伯特变换、补偿圆偿、绘曲线等模块。软件包设有HELP文件,以帮助用户正确地使用。计算结果可用打印机打出,也可用屏幕显示。更值得一提的是,该软件提供绘曲线程序,计算所得的导数值可用导函数的曲线形式在屏幕上显示,存绘图仪上绘制导数曲线图,也可以在打字机上输出导数曲线图。不仅如此,他是按照国家标准局(现国家技术监督局标准化司)已经批准的两项计算机软件国家标准: 相似文献
20.
综合利用总场异常及其垂直梯度的快速曲化平方法 总被引：1，自引：0，他引：1
摘 要 结合观测异常及其垂直梯度的特点给出了一种新的快速曲化平方法。这种方法利
用了位场的调和性质借助了三次样条优良拟合特点。模型计算表明该方法计算速度快
精度较高充分发挥了垂直梯度反映地质体垂向变化的特点。 相似文献