共查询到20条相似文献，搜索用时 453 毫秒
1.
曲面位场导数的换算方法大多基于调和函数的积分方程数值解。本文提出了基于B样条函数的曲面上位场导数直接计算法。引入按弧长S作变量的三次样条插值,得出了位场水平导数的计算式。该方法还有一个特点,那就是可以用一元样条来实现曲面上的求导。这些,都决定了该方法的简便性和精确性。利用曲线上和曲面上的两种理论模型进行试算,分别计算了Zα/x和Δg/Z~2计算和理论的结果对比表明,符合很好。误差主要取决于曲面起伏大小,异常变化陡度和取样的间隔。 相似文献
2.
《物探化探计算技术》1991,(2)
CMD软件包采用模块结构,用FORTRAN77语言编写,软件包的主要功能:①高精度频率域求导;②不同精度的空间域求异;③三次样条求导;④位场横向数据的垂向导数;⑤弱异常信息的增强。为完成上述功能,模块中有快速富里埃变换、希尔伯特变换、补偿圆偿、绘曲线等模块。软件包设有HELP文件,以帮助用户正确地使用。计算结果可用打印机打出,也可用屏幕显示。更值得一提的是,该软件提供绘曲线程序,计算所得的导数值可用导函数的曲线形式在屏幕上显示,存绘图仪上绘制导数曲线图,也可以在打字机上输出导数曲线图。不仅如此,他是按照国家标准局(现国家技术监督局标准化司)已经批准的两项计算机软件国家标准: 相似文献
3.
本文用样条函数法实现了二维位场的向上延拓和向下延拓,向上延拓用样条函数法计算泊松积分,向下延拓用样条函数法求解第一类Ｆｒｅｄｈｏｌｍ积分方程。模型和实例计算表明,本文方法正确可行且精度高。 相似文献
4.
本文用样条函数法实现了二维位场的向上延拓和向下延拓，向上延拓用样条函数法计算泊松积分，向下延拓用样条函数法求解第一类Ｆｒｅｄｈｏｌｍｔｍ积分方程。模型和实例计算表明，本文方法正确可行且精度高。 相似文献
5.
导数换算作为常用的位场数据处理手段具有重要的物理意义。由于实际测量误差和噪声的影响,波数域导数换算存在计算过程不稳定、求导结果精度差的问题。为了降低噪声干扰,提高导数计算精度,提出了在波数域常规导数算子基础上附加Chebyshev低通滤波器的任意阶导数换算方法。通过分析该滤波器的滤波特性,并结合位场异常径向平均功率谱曲线特征,确定了Chebyshev低通滤波法的滤波参数,减少了人为因素对求导结果的干扰。二度体和三度体模型试验表明,与其他3种方法(常规FFT(fast Fourier transformation)法、向上延拓法和ISVD(integrated second vertical derivative)法)相比,该方法计算的导数与理论值的均方根误差最小,求导结果精度较高。在虎林盆地布格重力异常数据处理中,利用该方法计算的垂向二阶导数受噪声干扰小,结果可靠性较高,依此划分出的8条大断裂和11条小断裂在以往研究中均得到证实。 相似文献
6.
垂向导数广泛应用于位场数据处理和解释中,常规的空间域求导公式多以台劳级数为基础,这不可避免会引入截断误差,且不便于求取高阶导数。另一种应用较多的波数域算法。其求导算子在高波数成份产生振荡失真,很不稳定,在实际应用中很难取得理想结果。这里提出了一种计算位场数据垂向导数的稳定算法,在波数域求垂向积分,在空间域计算二个二阶水平方向导数,借助拉普拉斯方程,求取位场异常的各阶垂向导数。算法只在求垂向积分时引入了一次付里叶变换,其它计算全部在空间域中进行,从而有效地避免了吉布斯效应的重复累加传递,有效降低了常规波数域求导算子在高波数成份产生的振荡失真。通过理论模型对比,证明了这里所提出的算法较常规的波数域方法具有良好的稳定性,特别是在计算高阶垂向导数时有明显优势,能有效地提高与垂向导数有关的位场定量解释方法的准确性。 相似文献
7.
8.
《水文地质工程地质》2015,(5):1-5
求解含水层中溶质和热量的运移时,连续的达西速度场十分重要。传统有限单元法能够精确求解水头,但其水头一阶导数在节点上不连续,求解通过截面流量时,流入量和流出量不相等。将有限单元法与三次样条技术结合,利用三次样条函数逼近水头,保证水头一阶导数的连续性,从而利用达西公式求解出连续的达西速度场,保证流量的连续性。通过对均质和非均质介质下的二维稳定流的数值模拟,发现这种结合三次样条技术的有限单元法能够精确求解水头、达西速度与流量。 相似文献
9.
在位场数据处理中,垂向导数具有重要的物理意义。其在一定程度上可以划分不同深度和大小异常源产生的叠加异常,且导数的阶次越高,这种分辨能力就越强,但通常认为高阶导数的换算是不稳定的。本文在Tikhonov正则化求位场垂向高阶导数的基础上,结合迭代法进行逐次逼近,提出了位场高阶导数的Tikhonov正则化迭代法,并且得到Tikhonov正则化迭代法的递推公式。通过对该方法的滤波特性分析可以看出,该方法计算的位场垂向高阶导数具有一定的稳定性及保幅性。模型试验和实际数据的处理表明,该方法计算结果较常规FFT求导法有更高的稳定性和实用价值。 相似文献
10.
在地球物理参数反演中，观测的物理场量，例如波场中的位移或电流场中的电位，对于介质参数的编数的计算是非常重要的。本文采用自伴随算子和格林函数的方法给出了弹性波场和稳恒电流场问题的Ｆｒｅｃｈｅｔ导数解析式。该式表明，Ｆｒｅｃｈｅｔ导数可由格林函数Ｇ及其梯度ΔＧ求得。把任意场源正演问题的解和Ｆｒｅｃｈｅｔ导数的计算统一起来，并提出用有限元方法计算任意背景场下的格林函数，从而得到Ｆｒｅｃｈｅｔ导数的数值解 相似文献
11.
在地球物理参数反演中，观测的物理场量，例如波场中的位移或电流场中的电位，对于介质参数的偏导数的计算是非常重要的。本文采用自伴随算子和格林函数的方法给出了弹性波场和稳恒电流场问题的Ｆｒｅｃｈｅｔ导数解析式。该式表明，Ｆｒｅｃｈｅｔ导数可由格林函数Ｇ及其梯度△Ｇ求得。把任意场源正演问题的解和Ｆｅｅｃｈｅｔ导数的计算统一起来，并提出用有限元方法计算任意背景场下的格林函数，从而得到Ｆｒｅｃｈｅｔ导数的数值解。文中以稳恒电流场问题为例，给出了直流电法视电阻率偏导数的数值结果。 相似文献
12.
13.
杨辉 《物探化探计算技术》1998,20(4):347-350
根据台劳级数展开式,将窗口内网格结点位场数据在计算点展开,且对应点两两相加,由最小二乘法求得窗口内网格结点的圆滑滤波算子和二次导数算子,将此算子与观测的位场数据褶积,就可得到位场圆滑后的数据及水平、垂直二次导数。 相似文献
14.
极限平衡分析方法是斜坡稳定性评价中的常用方法,在长期的工程实践中积累了丰富的经验,但其不能考虑斜坡岩土体中实际存在的不确定性,在应用中具有一定的局限性。可靠度分析方法可有效地考虑斜坡系统内的不确定性和相关性,但因状态函数偏导数的求解比较困难,使可靠度分析方法在实际中应用不便。为解决上述问题,根据二元函数插值逼近原理,在矩形区域上构造拉格朗日不完全双二次多项式逼近状态函数,从而近似地计算状态函数的偏导数,求得状态函数的均值和方差,并利用精度较高的一次二阶矩方法来计算斜坡的可靠指标和破坏概率。据鄂西恩施地区马堡营滑坡实例分析表明,引入二元函数插值逼近的一次二阶矩方法计算结果与剩余推力法及Monte-carlo模拟方法结果一致,其精度可满足工程需求。 相似文献
15.
扶永铭 《成都理工学院学报》1995,22(3):109-114
该文研究了低频人工时变场（似稳场）中任意导电导磁矿体的边界元素法。通过对矢量位的边界积分方程，再以样条边界元素法进行离散，得到实分量的常系数方程组和虚分量的二阶微分方程组。此方程可利用时间积分的样条函数方法求解，解出适量位后，利用磁场强度与矢量位Ａ的相应关系而求得二次磁场的水平分量和垂直分量。 相似文献
16.
综合利用总场异常及其垂直梯度的快速曲化平方法 总被引：1，自引：0，他引：1
摘 要 结合观测异常及其垂直梯度的特点给出了一种新的快速曲化平方法。这种方法利
用了位场的调和性质借助了三次样条优良拟合特点。模型计算表明该方法计算速度快
精度较高充分发挥了垂直梯度反映地质体垂向变化的特点。 相似文献
17.
引言众所周知,重力垂直二次导数是布伽异常沿垂直方向的二次导数,也就是重力位沿垂直方向的三次导数,是划分区域场和局部异常的有效方法.重力垂直二次导数有许多计算公式,所得结果不同,原因何在,究竟它们的意义是什么?它们之间有什么异同?能否提出一个准则来选择某一公式?计算精度如何评定?都是值得我们探讨的.本文结合一些实际成果粗略谈几点认识,抛砖引玉,请同志们批评指正.一二次导数法的一般意义大家知道,布伽异常中包含着比局部异常大得多的区域性异常,而和矿产分布有关的局部异常,往往被其掩盖或发生畸变.过去,我们习惯于从实测布伽异常图上来分辨局部异常,就是通常所讲的封闭和凸扭.这种 相似文献
18.
根据位场理论,将观测曲面上测得的△T(或Za)视为位函数,代入偶层位公式的极限形式,求解出对应于△T(或Za)的等效偶层的磁化强度J(α,β,γ)。将J(α,β,γ,)代入偶层场强分量表达式 (α)/(α)n(1/r)dS 只要给定了方向u,就可得到△T沿该方向的一阶导数。将J(α,β,γ)代入(α)T/(α)x,(α)T/(α)y,(α)T/(α)z的积分表达式就能得到△T的二阶导数。因为使用的是二重积分形式的导数的精确表达式,它克服了空间域差商求导的致命弱点。特别是该方法能将起伏地形曲面上的场值向上延拓(包括曲化平)和求一阶(或二阶)导数两个转换过程集中在一个积分表达式中,一次计算完毕。 相似文献
19.
在电子计算机技术应用于物探资料处理以前,垂向二次导数的计算主要用于消除重力资料中的区域场以突出局部场和分解多体场源的迭加异常。随着频谱分析方法和电子计算机技术在物探资料解释中的逐步推广使用,垂向二次导数异常的计算成了磁异常数据处理和解释主要内容之一。垂向二次导数异常为磁异常的推断解释提供某些与地质体形态更为直接相关的信息,而大量的磁异常资料又为垂向二次导数方法的应用提供广阔的前景。从而,垂向二次导数方法又重新引起人们的兴趣和重视,近年来,国内就发表了不少关于垂向二次导数物 相似文献