共查询到20条相似文献，搜索用时 231 毫秒
1.
《物探化探计算技术》2015,(6)
在瞬变电磁法正演计算中,常常需要频时转换来获得垂直磁场随时间的变化率。目前已经提出了多种频时转换方法,其中G-S变换法只涉及实数计算,具有较快的速度,因此已经成为航空瞬变电磁正演模拟的候选方法。由于航空瞬变电磁数据巨大,时间数较多,因此在保证精度的条件下怎样提高计算速度成为当前亟待解决的问题。鉴于大量频点的正弦变换和余弦变换具有较高的精度,这里使用三层和四层典型地电模型,通过比较余弦变换、正弦变换和不同点数G-S变换的正演计算结果来确定最佳的G-S变换系数。实验结果比较和分析认为,采用8个系数的G-S变换法不仅有较快的计算速度,而且也有较高的计算精度。 相似文献
2.
3.
《物探化探计算技术》2016,(4)
在前人的工作基础上,给出一种新的快速计算框线源激发的瞬变电磁场的数值滤波算法,其中内层积分的Hankel积分式采用47点J1型线性滤波系数,外层积分采用"n点式"Gauss-Legendre数值积分法。通过线性迭加技术,将多个单线源组合成任意多边形的瞬变电磁的发射框源。算例分析表明,数值计算结果和解析解的最大相对误差小于0.007%,单测点、单时间道的计算耗时为0.1s。该方法可应用于瞬变电磁法生产实际中确定出最小延迟时间,以避免"框线效应"对瞬变电磁场的畸变影响。 相似文献
4.
5.
国外瞬变电磁法中正反演问题研究现状 总被引：2，自引：0，他引：2
国外在瞬变电磁法正反演问题研究中,运用有限差分法、有限单元法及积分方程法直接在时间域中求解以及运用富氏变换法、衰减谱法和数字滤波法从频率域变换等数值计算方法,并对其在瞬变电磁资料解释中的效果作了评述. 相似文献
6.
汉克尔变换的数值计算方法是计算电磁响应的重要工具.这里重点对比了五种已知汉克尔变换数值滤波方法在不同特性积分函数下的计算精度,以及在电磁法探测计算中的应用.结果表明,五种滤波方法均能较好地应用于电磁探测正演计算,计算精度属Anderson 801滤波较高. 相似文献
7.
《物探化探计算技术》2021,(3)
地空瞬变电磁法(Ground-airborne transient electromagnetic method)是一种新兴的物探方法,其探测范围大,工作效率高,在地形复杂的区域有着传统物探方法无法比拟的优势。这里通过对地空瞬变电磁的一维正演问题的研究,采用Hankel变换、Gauss积分和频-时转换实现了地空瞬变电磁时间域响应正演计算,通过解析解和数值解的对比验证了正演方法的准确性和稳定性,并研究了其响应特征。利用最小构造模型约束下的自适应正则化反演方法,进行了地空瞬变电磁法反演模拟研究,证实了该方法对低阻层的反映比较真实,对高阻层的反映略有误差的特点。最后用该方法对某滑坡实测数据进行了反演解释,反演结果与工区地质情况符合较好。 相似文献
8.
体积分方程法模拟复电阻率三维体电磁响应 总被引：1，自引：0，他引：1
利用体积分方程法计算了均匀半空间中复电阻率(激电)效应和电磁效应同时存在时的三维体响应。在计算中,对于需要计算三重积分的张量格林函数电荷项一次部分,应用一种差分近似的方法求解,这种方法在保证计算精度的同时更加便于计算机实现;采用二次剖分的算法解决了计算张量格林函数时的奇异值问题;计算含有贝塞尔函数的积分项时,利用一种结合连分式展开的高斯求积代替常规的快速汉克尔变换方法。验证了计算结果并分析了三维体复电阻率模型(Cole-Cole)参数对正演结果的影响,为三维体复电阻率及其参数反演提供了正演依据。 相似文献
9.
《物探化探计算技术》2017,(2)
针对航空瞬变电磁法正演算法结构与GPU的并行性运算相结合,推导出了航空瞬变电磁法一维正演公式,在2007年NVIDIA公司推出的CUDA平台上实现了一维正演算法在GPU上的调用。为验证CUDA C语言程序的正确性,设计了几个模型,并与C代码模拟的响应做了对比,结果表明,基于GPU的CUDA C代码所编的程序在航空瞬变电磁法的模拟中拥有明显优势,速度提高非常明显,在复杂的地球物理学模型计算和反演中,拥有巨大的发展潜力。 相似文献
10.
关于瞬变电磁法2.5维正演中的几个问题 总被引：3，自引：0，他引：3
熊彬 《物探化探计算技术》2006,28(2):124-128
这里讨论了瞬变电磁法2.5维正演模拟中的两个积分变换:傅里叶逆变换和拉普拉斯逆变换。针对发收距为零的中心回线方式瞬变电磁法,提出了在傅氏域中考察傅氏变换函数随波数的变化规律,进而根据曲线首尾支渐近线来划定波数覆盖的范围,然后同解析解对比,确定出最少个数的傅氏域波数。另外,介绍了只需对较少的拉氏变换变量作纯实数运算的拉普拉斯数值反演计算方法。通过对均匀半空间表面上垂直磁偶极子源形成的瞬变电磁场进行正演模拟,结果表明,关于傅氏域波数的范围划定原则及个数选取方案是合理的、拉氏逆变换算法也是切实可行的。 相似文献
11.
强震作用下土石坝极易出现失稳破坏,从而造成人员伤亡和较大的社会经济损失.由于地震的不确定性,强震作用下土石坝失稳分析通常采用失稳概率表示,目前常用方法是地震易损性分析方法,主要有云图法和增量动力分析（incremental dynamic analysis,IDA）两种方法.IDA方法计算结果准确,但计算效率低,云图法计算效率虽高,但计算精度无法得到有效保证.基于上述问题,提出了一种基于云图法和IDA方法的地震易损性快速精准分析方法（CIHA,cloud-IDA hybrid approach）.CIHA方法可兼顾计算效率和计算精度,该方法基于云图法的对数线性回归假设,通过非线性时程分析,并对地震波进行一次放缩来计算相应损伤指标下的地震动强度值,利用地震动强度值得到的均值和方差生成土石坝在各个损伤等级下的易损性曲线.通过对Lower San Fernando土石坝的地震易损性分析,将所提CIHA方法与IDA方法的计算结果进行了对比.结果表明,在计算精度方面,CIHA方法可以获得与IDA方法相近的结果,在计算效率方面,CIHA方法相比IDA方法计算效率有显著提高. 相似文献
12.
13.
双重贝寒尔函数积分的数值计算 总被引：2，自引：1，他引：2
双重贝塞尔函数积分则由于贝塞尔函数互乘项的强振荡与慢衰减特性而以以应用通常的数值积分算法。本文将被积区间[0,∞]划分为[0，λ0]、[λ0，∞]两部分，应用贝塞尔函数的克尔函数表述式及后者的大宗量渐近特性，区间[λ0，∞]的双重贝寒尔函数积分可被转化为Fourier正（余）弦变换，并可利用各种快速算法对其进行数值计算；区间[0，λ0]上的双重贝塞尔函数积分的计算可直接应用一般的数值积分算法并能获得较高的计算精度；当需大量计算有共同参量的双重贝塞 尔函数积分时，其计算效率仍显不足。此时，可应用贝塞尔函数的导数关系式对[0,λ0]内的双重贝塞尔函数积分进行恒等变换，再用差商近似导数，将其转化为对贝塞尔函数本身的积分，而该积分又仅需计算一次。故本算法对双重贝塞尔函数积分的计算效率有明显提高。 相似文献
14.
Dealing with Zeros and Missing Values in Compositional Data Sets Using Nonparametric Imputation 总被引：5，自引：0，他引：5
The statistical analysis of compositional data based on logratios of parts is not suitable when zeros are present in a data set. Nevertheless, if there is interest in using this modeling approach, several strategies have been published in the specialized literature which can be used. In particular, substitution or imputation strategies are available for rounded zeros. In this paper, existing nonparametric imputation methods—both for the additive and the multiplicative approach—are revised and essential properties of the last method are given. For missing values a generalization of the multiplicative approach is proposed. 相似文献
15.
随着城市地铁工程的高速发展,地面建筑物基坑的施工将必然会对邻近的地铁隧道造成一定的影响,运用Midas GTS三维数值模拟计算软件分析基坑开挖对邻近地铁隧道的影响,对地层自重固结、基坑开挖施工的整个过程进行模拟分析,计算的结果与实际监测数据进行对比表明,该方法对实际工程有一定的指导意义。 相似文献
16.
17.
数值积分是基于节点计算的数值分析方法的关键技术之一。针对完全无网格法的计算要求，提出了一种能自动适应任意计算域上各种节点分布方式的数值积分算法。该算法能随计算点的位置不同，自动确定积分域及积分域内的求积点，计算出数值积分结果，给出了临空边坡滑移体和圆孔方板2个实例。结果表明，按文中算法编制的程序具有良好的边界再生能力和较高的计算精度，可作为基于节点数值分析方法的一个程序模块，用于实际计算。底层算法思想，使得程序可方便地转换为其他计算语言。 相似文献
19.
结合边界元法和离散元法,提出一种可以进行计算颗粒内部应力和破碎路径的方法。该方法利用离散元法求解颗粒的相互作用和每个颗粒上的荷载。然后利用边界元法计算颗粒的应力分布,为了实现动态平衡,将颗粒的加速度视为恒定大小的体力。但体力导致边界积分方程中出现域积分,故采用直线积分法将域积分转化为边界积分,以保证边界元法降维的优势。为了提高边界元的计算效率,对于几何形状相似的颗粒,以其中一个颗粒作为模板颗粒,只需要计算模板颗粒在局部坐标系中的系数矩阵,其他相似颗粒可以通过局部和全局坐标系之间的映射获得。在得到应力后,基于Hoek-Brown准则来判断颗粒是否破碎。此外,将破坏路径简化为直线,并采用最小二乘法拟合得到破坏路径。 相似文献
20.
地空瞬变电磁法融合了地面与航空瞬变电磁的方法原理与观测模式,将发射源置于地表,并利用航空器在空中接收响应信号。因此,地空瞬变电磁法同时具有地面装置大测深、高信噪的优势以及航空装置中采集高效的优势。先有研究偏重于单一辐射源,忽略了多辐射场源在勘探深度与信噪比上的优势。因此需要研究适用于多辐射场源地空装置的数据解释方法。在解决多辐射场源地空瞬变电磁正演问题之后,以Occam反演理论为基础,对多辐射场源地空瞬变电磁一维反演进行研究。首先要构建目标函数,本文通过引入拉格朗日乘子,将数据拟合差及模型粗糙度组合起来进而构建目标函数;其次解决灵敏度矩阵的构建问题。模型计算结果表明本文多辐射场源地空瞬变电磁法一维反演算法的有效性,可为三维反演与精细解释提供良好的基础。通过对典型模型的反演计算,从其理论模型及反演结果的对比可以看出,Occam反演方法可以应用到多辐射场源地空瞬变电磁系统中,为多辐射场源地空瞬变电磁解释增加了新途径。 相似文献