共查询到20条相似文献，搜索用时 604 毫秒
1.
《物探化探计算技术》2021,(4)
提出一种基于积分解的空间-波数混合域二度体磁异常数值模拟快速算法。该方法将磁异常二维空间域卷积问题,通过傅里叶变换转换为空间-波数混合域垂向一维积分问题,将一个复杂问题分解为多个小问题,不同波数的小问题之间具有高度并行性;保留深度方向为空间域,采用二次插值的形函数计算垂向一维积分,便于浅层网格适当加密,深层网格适当稀疏,兼顾计算精度、计算效率及模拟复杂地形。在此基础上,根据一维形函数积分的特点,提出了一种适用于起伏地形条件下的磁异常快速计算方法,核心思想是对于相同的单元积分进行存储,避免重复计算,进一步提高了计算效率,尤其适用于复杂地形条件下的模拟。模型算例中分别设计了突变介质模型、起伏地形模型和复杂模型,通过数值解与解析解对比结果表明:本算法正确、可靠,且具有数值模拟精度高、计算速度快,适用于任意复杂地形的特点。 相似文献
2.
垂向导数广泛应用于位场数据处理和解释中,常规的空间域求导公式多以台劳级数为基础,这不可避免会引入截断误差,且不便于求取高阶导数。另一种应用较多的波数域算法。其求导算子在高波数成份产生振荡失真,很不稳定,在实际应用中很难取得理想结果。这里提出了一种计算位场数据垂向导数的稳定算法,在波数域求垂向积分,在空间域计算二个二阶水平方向导数,借助拉普拉斯方程,求取位场异常的各阶垂向导数。算法只在求垂向积分时引入了一次付里叶变换,其它计算全部在空间域中进行,从而有效地避免了吉布斯效应的重复累加传递,有效降低了常规波数域求导算子在高波数成份产生的振荡失真。通过理论模型对比,证明了这里所提出的算法较常规的波数域方法具有良好的稳定性,特别是在计算高阶垂向导数时有明显优势,能有效地提高与垂向导数有关的位场定量解释方法的准确性。 相似文献
3.
利用有限元法实现了任意方向偶极子源在二维介质中频率域电磁响应的数值模拟,研究了波数取值对模拟结果的影响.通过对构造走向的Fourier变换,将全三维电磁问题,转化为一系列二维问题,并在波数域求解,极大地减小了计算工作量,导出了波数域耦合适用于二维电性介质中任何方向电或磁偶极子响应计算的电磁场方程.针对每个给定的波数,上述耦合电磁场方程用等参有限元方法在x-z平面内求解.采用Fourier逆变换,将波数域解积分,得到空间域电磁场.针对电磁模拟计算中,源点的奇异性,采用具有一定面积的伪δ函数表达源电流分布,使数值解精度得以提高.另外,采用等参有限元,使地下复杂地质体得到准确表达.利用不同波数值对均匀介质与层状介质的模拟结果与解析解的对比,验证了算法的正确性与精度.利用层状介质模型的解析解与数值计算结果的对比,分析了波数的优化取值范围及取值点数对数值模拟结果的影响,考察了算法对非均匀介质的适应性. 相似文献
4.
关于瞬变电磁法2.5维正演中的几个问题 总被引：3，自引：0，他引：3
熊彬 《物探化探计算技术》2006,28(2):124-128
这里讨论了瞬变电磁法2.5维正演模拟中的两个积分变换:傅里叶逆变换和拉普拉斯逆变换。针对发收距为零的中心回线方式瞬变电磁法,提出了在傅氏域中考察傅氏变换函数随波数的变化规律,进而根据曲线首尾支渐近线来划定波数覆盖的范围,然后同解析解对比,确定出最少个数的傅氏域波数。另外,介绍了只需对较少的拉氏变换变量作纯实数运算的拉普拉斯数值反演计算方法。通过对均匀半空间表面上垂直磁偶极子源形成的瞬变电磁场进行正演模拟,结果表明,关于傅氏域波数的范围划定原则及个数选取方案是合理的、拉氏逆变换算法也是切实可行的。 相似文献
5.
移动荷载作用下半无限体的动力响应解 总被引：3，自引：3，他引：3
以Duhamel积分为基础，采用Fourier变换和Floquet变换等方法，给出了移动荷载作用下半无限体中任一点的动力响应在时域、频域和频率-波数域内的表达式。该公式的推出，在交通荷载作用下，可为土体和结构的动力响应分析提供一定的理论基础。 相似文献
6.
动力响应问题的求解对于地基在外荷载作用下引起的弹性波动问题研究有重要的意义。本文提出了一种求解横观各向同性层状地基在施加时间简谐荷载作用下任意点的应力响应的算法。此算法利用傅里叶变换将广义平面应变问题频率-空间域的动力方程转化到频率-波数域内,结合对偶变量的引入,利用高精度的精细积分算法对状态方程进行求解,在得到频率-波数域内的位移响应的基础上,利用傅里叶逆变换得到任意点的动应力响应。简谐荷载不仅可以施加在地基表面,而且可以施加在地基内部。对比算例验证了本文算法的准确性,同时对地基各向异性特性、激励频率和阻尼比对动应力响应的影响进行了参数分析,为工程实际提供可靠的数值依据。 相似文献
7.
《物探化探计算技术》2017,(2)
解决任意密度分布复杂地质体重力异常三维正演快速、高精度计算问题,是实现重力三维反演、人机交互解释建模的关键。针对该问题,从积分方程出发,提出一种波数域重力异常三维正演方法,其关键环节包括三个方面:(1)将研究区域剖分成许多规则小棱柱体,每个小棱柱体密度值可以任意给定,以此刻画任意密度分布和起伏地形条件下的复杂地质体;(2)给出一种新的高精度均匀棱柱体重力异常二维波数域的计算公式,用于计算组合棱柱体模型的重力异常;(3)采用Gauss-FFT法将重力异常从波数域转换到空间域,保证计算效率的同时,有效克服了传统FFT法引起的边界效应问题。模型算例检验结果表明,该算法计算速度快、精度高,对于剖分为百万个棱柱体的模型,耗时只需几秒。 相似文献
8.
探讨了在波数域作静效应校正的可能性,提出利用各频点互相关矩阵作为选择滤波截止波数的理论依据,并对三维数值模拟结果作了处理,表明在波数域消除静效应是可行的。最后,对云南某地的CSAMT数据作了试算,得出了较好的效果。 相似文献
9.
虽然傅氏转换能够筒化波动方程式的解法，但在频率-波数领域中，求解波动方程式，计算震波模型，该方法不能适应传波介质之界面弯斜或传波速度横向变化。把界面视为一组不连续的绕射点，用各点反射系数作为讯号源强度，各绕射点分别在频率-波数领域中实施波场向上外插，再把各点外插波场相加，由此得到震波模型。利用相位转移函数，取代每一绕射点的傅氏转换，能节省大量电算时间。这种频率-波数领域震波模型法，比时间-空间领域法更容易处理复杂介质模型。用新的方法对几种典型之地质模型，包括断层、背斜、尖灭、多层簿地层岩相突变等进行了计算，所得震波特性相当合理，符合震波传播理论。实例表明，即使对于岩相突变，计算的震波模型与实际震波记录，亦相当吻合。频率-波数领域震波模型法，对于复杂界面和速度变化之地质模型，其数学处理上远比射线理论和克希霍夫积分方程来得可行且单纯。且比定差法和定元法显著节省网格点和容易处理边界条件。唯这种频率-波数领域绕射点法，随绕射点数之增加，其电算时间之增加幅度似乎较高。关键词震波模型傅氏转换波动方程式频率-波数领域波场外插相位转移函数 相似文献
10.
双重贝寒尔函数积分的数值计算 总被引：2，自引：1，他引：2
双重贝塞尔函数积分则由于贝塞尔函数互乘项的强振荡与慢衰减特性而以以应用通常的数值积分算法。本文将被积区间[0,∞]划分为[0，λ0]、[λ0，∞]两部分，应用贝塞尔函数的克尔函数表述式及后者的大宗量渐近特性，区间[λ0，∞]的双重贝寒尔函数积分可被转化为Fourier正（余）弦变换，并可利用各种快速算法对其进行数值计算；区间[0，λ0]上的双重贝塞尔函数积分的计算可直接应用一般的数值积分算法并能获得较高的计算精度；当需大量计算有共同参量的双重贝塞 尔函数积分时，其计算效率仍显不足。此时，可应用贝塞尔函数的导数关系式对[0,λ0]内的双重贝塞尔函数积分进行恒等变换，再用差商近似导数，将其转化为对贝塞尔函数本身的积分，而该积分又仅需计算一次。故本算法对双重贝塞尔函数积分的计算效率有明显提高。 相似文献
11.
最优化波数在点源二维有限元正演计算中的应用 总被引：3，自引：0，他引：3
文章结合地电模型给出了计算误差随最优化波数个数的变化曲线,揭示出该曲线上存在突降稳定点,这一规律使计算所需波数大大减少。研究还发现误差随网格加密参数曲线上也存在突降稳定点。上述规律的发现可以大大节省二维有限元正演的计算时间,使之在反演中的充分应用成为可能。 相似文献
12.
二阶弹性波动方程高精度交错网格波场分离数值模拟 总被引：4，自引：0，他引：4
给出了一种等价的二阶弹性波动方程,以解决弹性波场中完全弹性波动方程不能完全分离耦合的纵、横波波场问题.应用高阶交错网格有限差分法求解该波动方程,并使用通量校正技术(FCT)进一步压制频散,采用均匀介质模型和层状介质模型进行波场分离数值试验,精确得到了混合波场、完全分离的纯纵波及纯横波波场.数值结果分析表明,本文方法在均匀介质情况下准确可靠,在分离后的纯纵、横波波场中可观察到较为丰富的能量转换信息,这对认识复杂弹性波的传播规律及弹性波理论具有重要意义. 相似文献
13.
14.
Conventional investigations of waves–seabed interaction problems have been only concerned with the soil response due to two-dimensional linear progressive waves over a uniform seabed. However, the effects of non-linear waves which have been reported in the literature may be significantly different. In this paper, a finite element model is developed to investigate the non-linear wave-induced seabed response with variable permeability and shear modulus in a three-dimensional domain. The finite element formulations are fully presented in this paper. The numerical model is verified with the previous investigations through the reduced form of the present solution. The numerical results indicate that the influence of non-linear wave components cannot always be ignored without substantial error. Furthermore, the wave-induced seabed response is affected significantly by variable permeability in coarser seabeds and variable shear modulus in finer seabeds. © 1997 by John Wiley & Sons, Ltd. 相似文献
15.
有限元法是复杂介质地震模拟的有力工具,它能比较客观地反映地震波传播过程,比较细致地再现地震图像。但由于有限元法地震模型所需网格节点数多、数据规模大,使它的计算受到限制。采用对角化集中质量矩阵、刚度矩阵变二维为一维存贮,并通过使用有限差分吸收边界条件来处理地质模型的人为边界反射问题,从而使计算所需内存降低,计算时间减少,并能比较好地实现地震声波用有限元法计算,理论模型的计算结果验证了算法的正确性和可行性。 相似文献
16.
研究了不同测量方式(井-地,地-井,井-井)下点源场井中电法的三维有限元数值模拟。考虑到深度方向上大范围的网格剖分和井眼的影响及井-井测量等因素,采用放射状三棱柱单元的网格剖分方式,以提高网格质量,减少剖分单元数;给出了三棱柱单元的坐标变换公式,进行精确的单元积分,减少了单元积分时间;结合非结构化网格技术,实现了复杂模型的模拟;开发出相应的程序实现了复杂条件下(如考虑井眼影响、井井测量、倾斜井情形、地形起伏等)电法测井的三维有限元模拟,数值算例验证了方法的可靠性及计算效率,并对不同情形下的异常响应进行了分析,为进一步的反演工作奠定了基础。 相似文献
17.
杨奎 《物探化探计算技术》2012,34(4):436-443,368
面波勘探现在已经被广泛应用于工程选址及地震安全评价等方面,但其大多是针对Rayleigh波的应用,对Love波的关注较少.针对面波信息中的两种主要成份Rayleigh波和Love波对速度递增模型和存在低速层模型进行多模频散曲线的敏感性和特征分析,得到了多模式波对半空间介质模型和多层介质的响应特征,而且还对其基模式波和高模式波在高频和低频上的速度趋势进行了分析,细化了高频时的多模式波的理论最低速,得出了一些具有普遍性的结论,并为面波信息的综合利用提供了有力支持. 相似文献
18.
19.
用ＩＲ－１６００系列傅利叶变换红外线分光光度计对腐植煤系列的褐煤、长焰煤、气煤、工、焦煤、瘦煤、贫煤和无烟煤煤样的镜质作了测定，共红外光谱特征反映了它们在变质程度上的区别。得出的镜质体芳核Ｃ＝Ｃ键红外吸收波数随煤级变化曲线，反映腐植煤的芳核Ｃ＝Ｃ键红外吸收波 随煤化程度增加而降低或向低波娄迁移的规律。有理由相信、腐植煤的煤化程度或牌号煤也可以用红我光谱方法精确确定。 相似文献
20.
煤层工作面内陷落柱与煤层之间具有显著的电阻率、波速和密度差异,适宜于采用无线电波透视法和震波透视CT法来探测煤层工作面内陷落柱的边界范围;陷落柱体通常裂隙发育,与围岩相比较,其富水陷落柱为相对低电阻率值范围,而不富水陷落柱则为相对高电阻率值范围,因此,可利用煤层底板三维电法来探查底板陷落柱富水性。谢桥矿13218工作面1^#陷落柱的探测表明,该陷落柱范围具有较高的电磁波吸收系数特征和相对较高的纵波波速特征;工作面底板陷落柱向下范围为相对高电阻率值特征,为不富水陷落柱。该综合物探探测结果与实际验证资料基本吻合,有效地指导了本煤层工作面开采及底板煤层工作面的开采设计与施工。 相似文献