共查询到20条相似文献，搜索用时 875 毫秒
1.
本文利用第二代小波多尺度分解和快速变换的特点,构造自适应计算网格.对初始计算网格上的数值解进行第二代小波变换,得到数值解对应的小波系数空间.小波系数的大小表示相邻网格上数值变化率,小波系数大的区域网格点上的数值解变化梯度大.当小波系数大于等于预设的阈值时,在小波系数对应的网格点周围插入新的计算网格点,通过阈值可以实现网格的细化,得到多尺度下层层嵌套的细化自适应网格;由有限差分法得到相应网格点的空间导数.比较数值算例得到的波场快照和计算时间,验证了该方法的有效性. 相似文献
2.
本文将面向目标的自适应算法应用于三维大地电磁数值模拟.使用基于非结构网格的矢量有限单元法对起伏地表大地电磁正演模拟问题进行求解.使用利用垂向电流密度在物性界面上的连续性对后验误差进行估算的算法指导网格优化.由于全局自适应算法针对观测点优化网格的能力较差,本文通过求解正演问题的对偶问题计算后验误差的加权系数,并对相关加权系数进行改进,从而实现了面向目标的自适应算法.与传统基于结构化网格的电磁正演算法相比,采用非结构网格能够更好地拟合起伏地表和地下不规则异常体.由于使用了面向目标的自适应算法,本文能够使用更少的网格达到较高的计算精度.通过对比本文模拟结果与半空间响应和全局自适应算法计算结果,并通过对比使用改进前和改进后加权系数得到的网格剖分结果验证了本文算法的有效性. 相似文献
3.
二维频率空间域的数值模拟方法具有以下的优势:多炮模拟时,计算成本比时间域方法低;无累计误差;在地震反演中处理多震源模拟时,只需要有限的几个频率就可以得到好的反演结果.差分离散化形成的稀疏系数矩阵,需要求解一个巨大规模的线性方程组,最大瓶颈是需要海量的计算机内存,导致计算量庞大.本文在前人研究的基础上,采用嵌套剖分网格排序法,极大限度减少对计算机内存的需求,从而减少了计算量.针对弹性波数值模拟的特征,提出二维频率空间域弹性波多炮模拟的快速计算流程.数值模拟试验证明使用嵌套剖分排序法的弹性波多炮数值模拟比压缩存储法具有节省存储量、计算效率高等优势,为后续的二维频率空间域弹性波全波形反演奠定了很好的基础. 相似文献
4.
本文研究了二维多孔弹性波方程的多尺度波场数值模拟方法.该多尺度方法可采用较粗的网格计算,同时又能反映细尺度上物性参数的变化信息.文中详细阐述了多尺度模拟方法与算法,并推导了相应的计算格式.基本思想是建立粗细两套网格,在粗网格上,基于有限体积方法计算更新波场;在细网格上,计算多尺度基函数,这基于有限元方法通过求解一个局部化问题得到.对含有随机分布散射体的多孔介质模型进行了数值计算,计算中应用了完全匹配层(PML)吸收边界条件,数值结果验证了本文方法和算法的正确性和有效性. 相似文献
5.
将区域气候模式RegCM2与中国科学院大气物理研究所的9层全球格点大气环流模式IAP AGCM单向嵌套,对东亚现代气候进行数值模拟研究,同时检验和分析该嵌套模式的性能.已完成的10年积分结果表明,单向嵌套RegCM2由于具有较高分辨率和较完善的物理过程,因此对地面气温和降水的空间分布形势和季节变化趋势都有较好的模拟能力,且较与之嵌套的IAP AGCM的模拟效果有较大改善,如在中国区域,它模拟的年均地面气温与实况的空间相关系数由全球环流模式的092提高到094,模拟的年均降水由05提高到07. 这与嵌套RegCM2能模拟出IAP AGCM所不能分辨的中尺度信号有很大关系. 相似文献
6.
在采用有限差分方法开展探地雷达复杂目标体精细结构模拟时,为了提高计算精度,常采用非均匀网格对目标区域划分小尺寸的网格,以压制离散网格频散现象和保证有限差分方法的稳定性.常规非均匀网格和自适应亚网格技术在网格剖分数量和粗细网格边界处理上难以达到计算效率和计算精度的均衡.本文根据隐形斗篷(invisible cloak)理论,将基于变换光学(Transformation optics)理论应用于有限差分探地雷达数值计算中.该理论的主要思想是基于目标参数变化而保持电磁场的传播不变性,在坐标变换后,Maxwell方程的形式可以维持不变,而使得相对介电常数与磁导率的表达式变得复杂.通过这种方式可以虚拟地扩大目标体所占的网格节点数,减少背景介质区域的网格数,不增加模型空间的网格总数.另外,这种网格划分方式不但提高了计算效率,同时也可以克服亚网格技术边界反射误差的影响.本文推导实现了基于变换光学的二维有限差分方法,通过典型探地雷达模型测试,对比分析了该方法与常规有限差分、变网格有限差分和自适应亚网格有限差分的优缺点.计算结果验证了基于变换光学的有限差分可用于探地雷达目标精细结构模拟,具有较高的计算精度和计算效率. 相似文献
7.
本文提出了磁总场异常垂直梯度三维相关成像方法,用于成像地下等效磁源分布.它首先将地下待成像空间划分为三维规则网格,然后直接计算每个网格节点磁偶板子在观测面理论磁总场异常垂直梯度与实测磁总场异常垂直梯度的互相关,其相关系数值表征等效磁偶极子分布(即磁偶极子发生的概率).理论长方体组合模型数据和实际某矿区磁测资料试验结果表明本文方法计算得到的相关系数值能基本反映地下的磁源分布,且分辨率明显高于磁总场异常三维相关成像的分辨率,也高于基于熵滤波分离异常的磁总场异常三维相关成像的分辨率. 相似文献
8.
加权近似解析离散化(WNAD) 方法是近年发展的一种在粗网格步长条件下能有效压制数值频散的数值模拟技术. 在地震勘探的实际应用中, 不是所有情况都适合使用空间大网格步长. 为适应波场模拟的实际需要, 本文给出了求解波动方程的非一致网格上的WNAD算法. 这种方法在低速区、介质复杂区域使用细网格, 在其他区域采用粗网格计算. 在网格过渡区域, 根据近似解析离散化方法的特点, 采用了新的插值公式, 使用较少的网格点得到较高的插值精度. 数值算例表明, 非一致网格上的WNAD方法能够有效压制数值频散, 显著减少计算内存需求量和计算时间, 进一步提高了地震波场的数值模拟效率. 相似文献
9.
《地球物理学进展》2017,(3)
复杂地表条件下的有限差分地震波场的数值模拟,由于受到低速层和地表起伏的限制,模型速度分布范围变大,一般使用精细的差分网格来抑制频散,提高模拟分辨率,但精细网格会显著增加计算成本.为了能有效地解决这一问题,本文提出一种步长自适应有限差分波动方程数值模拟方法.(1)新方法根据模型中的介质速度分布,对不同的速度区域采用与该速度匹配的空间步长,实现对模型空间网格的步长自适应精细划分.对于速度分布范围大的复杂地表模型,新方法不仅能够极大地减少模型的网格节点数,同时又能提高波场的时间采样步长,减少时间采样数,提高计算效率.(2)推导了不同步长边界网格节点Laplace算子的二阶有限差分表达式,避免了在这些结点进行插值计算产生的假扰动和数值不稳定问题.(3)为了降低有限差分产生的数值频散,本文在常规的差分方程中增加了一频散校正项,能有效地衰减了高波数成分,抑制了数值频散.对复杂近地表的波场数值模拟结果表明,本文提出的步长自适应新方法能够有效减少网格节点数和时间采样数,极大地提高计算效率,计算量比常规粗网格增加一些,但效果能够达到了常规精细网格的模拟结果. 相似文献
10.
11.
12.
全波形反演是地震资料处理中速度建模的有力工具,相比层析成像等速度建模方法它能够得到速度场的更高频成分.本文给出了基于声波方程格子法正演的时间域全波形反演方法,该方法用非规则、非结构化的三角网格来离散计算区域及模型参数,能实现网格粒度与反演分辨率在空间上的自动匹配,内存需求少,计算效率高;采用L-BFGS优化方法,以分频段变网格的方式实施多尺度反演.以二维Overthrust模型进行了速度反演数值测试,显示了该方法的高效性和潜力. 相似文献
13.
有限差分法和有限单元法在大地电磁场数值模拟中已经得到了广泛的应用,但其数值结果的精度在很大程度上依赖于网格的离散程度.当模拟起伏地形、弯曲界面等复杂地电模型大地电磁场响应时,常常需要花费大量的时间以便得到较合理的离散网格.无网格局部Petrov-Galerkin法(MLPG)不同于有限差分法和有限元法,其形函数和权函数脱离了网格的束缚.本文详细推导了二维大地电磁场边值问题的弱式形式,并将其离散为局部积分域内的表达形式.通过模拟二维海洋地电模型大地电磁场响应,并与结构网格有限元结果进行对比,验证了本文算法和程序的正确性及精度.设计了一个含有弯曲界面的二维地电模型,讨论了不同离散网格对MLPG无网格法模拟结果的影响,并与结构有限元法结果进行了比较,结果表明MLPG无网格法模拟结果受离散网格影响较小.最后利用MLPG无网格法计算了两个海洋起伏地形模型的大地电磁响应,讨论了海底起伏地形对大地电磁响应的影响. 相似文献
14.
在高频面波方法中,水平分辨率是指水平方向上分辨异常体的能力.异常体在水平方向上的长度可用水平方向上横波速度的异常尺度来确定.面波多道分析（MASW）方法被广泛应用于浅地表横波速度结构的探测,然而该方法确定的横波速度是整个检波器排列的平均计算结果,因此水平分辨率较差.另外,采用共中心点（CMP）多次覆盖的方式采集数据亦增加了野外的工作量.我们在MASW方法的基础上,应用面波层析成像方法,提出一套提高面波勘探水平分辨率的完整方法的技术流程.首先,利用波场分离技术获得准确的基阶或高阶模式面波,采用相位扫描的互相关方法测量多道面波记录中任意两道之间的面波走时;然后根据面波层析成像方法,获得高分辨率的各目标网格内的纯路径相速度频散曲线;最后反演所有目标网格内的纯路径相速度频散曲线,得到研究区域的拟二维横波速度结构.这套方法具有一定的抗噪能力,理论上它可以准确地提取相邻两道之间面波的相速度频散曲线;同时由于该方法最少只需要1个排列就可以获得拟二维横波速度结构,因此它显著减小了野外工作量.理论模型和实际资料都证实了这套方法可有效提高面波勘探的水平分辨率. 相似文献
15.
坐标变换法通过将物理空间的曲网格映射为计算空间的矩形网格,将起伏地表转化为水平地表,同时将物理空间的波动方程转化为计算空间的波动方程,在计算空间完成数值模拟,坐标变换的方法对处理起伏自由边界具有较好的适应性和应用效果。本文在传统坐标变换方法的基础上,根据计算区域速度差异采用不同的网格大小和采样时间步长,提出了一种基于时空双变网格的起伏地表坐标变换正演模拟方法。在编程实现算法的基础上,通过典型模型波场模拟试算结果分析可知:(1)变网格方法与常规方法波场模拟误差在0.5%左右;(2)变网格方法计算效率视不同的变网格区域面积及变网格大小可提高几倍量级,在本文模型和计算参数下提高约5倍。(3)在满足模拟精度及频散条件要求下,变网格方法较全局细网格算法能显著节约计算内存。为此,针对起伏地表数值模拟,本文方法具有较高的模拟计算精度和一定的适应性。 相似文献
16.
本文基于自适应非结构有限元算法实现海洋电磁起伏海底地形三维正演模拟.通过采用隐式后推欧拉时间离散技术,保证在较大的时间步长条件下获得正确结果.为获得多时间道海洋电磁正演模拟的有效网格,我们采用基于法向电流连续的后验误差估计的自适应方法和网格融合技术;同时为了控制网格数量和保证正演模拟稳定性,我们还在网格融合过程中应用了随机网格挑选技术.对于方程组求解我们使用MUMPUS直接求解器.当时间步长不变时,只需对系数矩阵进行一次分解,大大提高计算速度.将本文计算结果与半空间模型解析解进行对比,验证了本文算法精度.针对海洋电磁半拖曳式和双船拖曳式工作方式,我们通过典型模型的模拟计算,研究海底地形影响及海底高阻层识别特征. 相似文献
17.
基于第二代小波变换的提升方案构造了插值小波,将雷达波场函数进行了二维小波变换,得到所有尺度上与计算网格相联系的小波系数和尺度系数.对所有尺度上的小波系数进行分析,根据解的局部性与小波系数阈值的控制,实现网格压缩和配点的自适应调节.保留大于给定阈值的小波系数及对应网格点,令小于给定阈值的小波系数为零,并舍弃其对应网格点.达到光滑区域采用较少的计算网格点,在奇异性较大的区域采用较多的计算网格点的目的.通过对自适应网格进行邻域校正、重构检查等附加修正,推导了场值更新的显式时间迭代方案.最后,以均匀、阶梯与复杂三个典型GPR模型为例,与常规数值计算结果对比表明：自适应小波配点法（AWCM）利用第二代小波的多尺度分解和快速变换的特点,可以使计算网格随着时间步适应解的移动和变化,允许计算资源更有效地使用,具有高压缩率,达到跟踪奇异性的目的,特别适合于探地雷达正演中波传问题的模拟. 相似文献
18.
本文开发了基于非结构化四面体网格的三维矢量有限元法,实现了固体矿床井眼中多分量感应测井响应的数值模拟,并分析了多分量感应测井仪器在复杂矿床模型中的响应特征.本文通过采用几何因子背景场,有效地避免了源的奇异性问题;同时,在井眼边界采用非均质网格并用Gauss-Legendre积分计算四面体单元的等效电导率.利用LU分解求解线性方程组,实现了一次网格划分多点的数值计算,提高了计算效率,从而实现快速连续的多分量感应测井模拟.非结构化的四面体网格确保了该方法可以模拟实际问题中所能遇到的复杂的矿体模型.基于水平三层分层和径向分层模型,验证了算法在各向同性和各向异性介质中的可靠性.我们还以三个不同的矿床模型为例,研究了多分量感应测井仪的不同分量的探测特性,结果表明,结合九个分量的信息,可以探测矿体的深度,也可以识别矿体的方位和走向,为精确地描述矿体的三维分布特征打下了基础. 相似文献
19.
三角网格有限元法具有网格剖分的灵活性,能有效模拟地震波在复杂介质中的传播.但传统有限元法用于地震波场模拟时计算效率较低,消耗较大计算资源.本文采用改进的核矩阵存储(IKMS)策略以提高有限元法的计算效率,该方法不用组合总体刚度矩阵,且相比于常规有限元法节省成倍的内存.对于时间离散,将有限元离散后的地震波运动方程变换至Hamilton体系,在显式二阶辛Runge-Kutta-Nystr9m(RKN)格式的基础之上加入额外空间离散算子构造修正辛差分格式,通过Taylor展开式得到具有四阶时间精度时间格式,且辛系数全为正数.本文从理论上分析了时空改进方法相比传统辛-有限元方法在频散压制、稳定性提升等方面的优势.数值算例进一步证实本方法具有内存消耗少、稳定性强和数值频散弱等优点. 相似文献
20.
为了实现大地电磁测深三维高效快速正演算法,本文研究了大地电磁测深积分方程三维正演模拟的理论与方法.首先详细地推导了Fredholm型积分方程;然后给出积分方程的离散方法、并矢格林函数积分的求解方法以及大型矩阵方程的存储与求解方法;最后,通过理论地电模型的计算检验了方法和程序编制的正确性,并展示了几种理论地电模型的数值计算结果,分析其异常特征.算例结果表明,该正演算法在CPU主频为1.81GHz、内存为1GB的计算机上能正确地模拟多个异常体的异常特征;地表电性不均匀体会使地表视电阻率曲线产生十分严重的畸变;TE模式相对TM模式纵向分辨率比较好,而TM模式的横向分辨率较好;随着异常体埋深的增加,反映异常的分辨能力进一步降低. 相似文献