共查询到20条相似文献,搜索用时 421 毫秒
1.
《物探化探计算技术》2017,(5)
三维地震波动方程有限差分数值模拟,受到计算量大且内存消耗大地制约而计算效率较低,严重影响其发展和应用。目前,一种全新的异构众核协处理器,为三维并行计算提供了一条解决计算效率和数据存储问题的重要途径。基于三维有限差分计算的特点,介绍了面向MIC(Many Integrated Core)平台的三维瑞雷波有限差分并行算法移植和性能优化,采用ofload模式将计算核函数加载到MIC上,在MIC协处理上使用openMP多线程并行方式,并通过循环合并、nocopy数据传输、SIMD向量化和CPU+MIC协同计算方式进行优化。通过模型计算可知:三维数值计算程序在MIC上具有近线性加速比,MIC+CPU协同并行计算的性能可以达到纯CPU节点的4.5倍,大大提高了计算效率。通过对模拟结果与解析解的比较,以及对频散曲线特征的分析,验证了数值模拟的正确性。 相似文献
2.
瞬变电磁法是探测煤田采空区的一种有效方法。瞬变电磁拟地震处理解释方法研究日渐成为电磁勘探领域的热点并且取得了突破性进展,所以可以把瞬变电磁拟地震处理解释运用到煤田采空区的探测中来提高勘探精度。文中在瞬变电磁三维时域有限差分数值模拟两层含水采空区模型的基础上,运用预条件正则化共轭梯度算法实现了瞬变电磁扩散场到虚拟波场的稳定转化,并对虚拟波场进行了脉冲压缩处理;最后通过克希霍夫积分解,对虚拟波场进行延拓成像,实现了两层含水采空区模型的三维解释。 相似文献
3.
根据汉口新业大厦水文地质与工程地质条件,建立开采地下水非稳定流定解问题三维数学模型,利用现有承压含水层非稳定流定解问题数值方案求解,重点讨论了源汇项的具体数学表达式。又简介根据达西定律和水均衡原理,建立了三维有限差分数值方程及其解法。根据土的有效应力原理和压密原理推导出了抽水沉降非线性压密方程,并研制了相应的三维有限元和有限差分计算程序。最后,对该环境岩土工程问题进行了分析和评价,计算结果与实测结果符合较好。 相似文献
4.
为高效求解扩展型Boussinesq水波方程,建立了基于有限差分和有限体积方法的混合数值格式。将一维控制方程写为守恒形式,方程中通量部分采用有限体积方法求解,剩余部分采用有限差分方法求解。其中,有限体积方法采用Godunov类高分辨率格式,并结合HLL(Harten-Lax and van Leer)式黎曼问题近似解求界面数值通量,黎曼问题界面左右变量通过高精度状态插值方法(MUSCL)构筑。有限差分方法则采用具有二阶精度的中心差分公式进行。采用具有TVD(Total Variation Diminishing)性质的三阶龙格-库塔多步积分法进行时间积分。对数值模式进行了验证,数值结果同解析解或实验数据吻合良好。 相似文献
5.
6.
基于球坐标系下有限差分的地磁测深三维正演 总被引:2,自引:0,他引:2
为了计算全球尺度电磁感应的响应,本文介绍地磁测深频率域三维正演。正演算法采用球坐标系下的交错网格有限差分方法,从Maxwell方程的积分形式出发,采用PARDISO对离散后的方程组求解,避免了迭代求解的散度校正。为了验证本文结果的正确性和精度,与前人的有限元和有限差分方法进行了对比,一维层状模型的三维交错网格有限差分数值结果和解析解相对误差小于5%,双半球模型的计算结果与前人的计算结果完全吻合。三维"棋盘模型"计算表明磁场分量对异常体的大小和位置具有很好的分辨能力。 相似文献
7.
8.
介质的各向异性会引起大地电磁观测曲线的畸变,对大地电磁数据解释影响很大。这里从麦克斯韦方程组出发,推导了基于二阶插值基函数的二维大地电磁非结构有限元公式,实现了任意地形起伏条件下任意各向异性二维大地电磁响应的数值模拟。各向异性断层模型的数值模拟结果与解析解吻合一致,复杂各向异性模型的有限元模拟与已有的有限差分计算结果也相符合,验证了算法的正确性。在此基础上探讨了各向异性参数对二维大地电磁响应的影响。 相似文献
9.
10.
《物探化探计算技术》2017,(1)
利用交错网格有限差分技术实现了频率域海洋可控源电磁三维正演算法,并探讨了海底地形对海洋可控源电磁响应的影响和校正方法及效果。通过与海底二维山峰地形的二维自适应有限元解对比,验证了文中算法的正确性和模拟海底起伏地形的可行性,并分析了海底二维山峰地形对三维油气储层电场响应的影响。关于海底起伏地形对海洋可控源电场响应的影响,采用比较法对其进行了地形校正,结果表明,在深水区比较法可以很好地消除地形影响,恢复海底高阻油气储层的电场响应。 相似文献
11.
12.
分布荷载推力桩计算的p-y曲线法研究 总被引:2,自引:2,他引:0
为提高现有水平梯形分布荷载推力桩设计计算水平,提出了地基系数按非线性土抗力.水平位移(p-y)曲线抗力模式表达的水平梯形分布荷载推力桩位移和内力计算的有限差分数值分析方法,并详细推导了桩身位移的差分格式.基于这些公式,编制了全桩位移、内力及侧土抗力的计算和图形处理程序,可适用于滑坡抗滑桩和深基坑悬臂支护桩的设计计算.结合新型桩型现浇混凝土薄壁管桩(PCC桩)作为滑坡抗滑桩进行对比分析,算例表明:该方法方便可靠.当有限差分段划分得足够小时,可使数值解接近于真实解. 相似文献
13.
ZTEM(Z轴倾子电磁法)是一种天然场源的频率域航空电磁法,其特点是磁场垂直分量在空中机载平台测量,磁场水平分量在地面的固定基站测量,具有勘探深度大、速度快、成本低、覆盖面积大等技术优势。本文实现了ZTEM二维有限差分正演和二维非线性共轭梯度(NLCG)反演算法。研究对象是倾子资料,反演过程中通过解"拟正演"问题来避开雅克比矩阵的直接计算。通过理论模型合成数据反演试算,验证了ZTEM倾子资料二维NLCG反演算法的稳定性与可靠性。与大地电磁(MT) TE模式阻抗资料反演结果进行对比,发现在异常体横向边界的约束方面,ZTEM倾子反演比MT阻抗反演更具优越性。 相似文献
14.
15.
为了研究复杂地电模型的航空瞬变电磁法全波形响应特征,需要开发考虑发射波形的三维数值模拟算法。本研究基于非结构四面体网格和位移逆Krylov子空间(Shift-and-Invert Krylov,简称SAI Krylov)方法,采用基于电偶极子离散的场源处理方法模拟场源,在时间域进行计算实现了全波形航空瞬变电磁法矢量有限元三维数值模拟。使用均匀半空间模型在阶跃波、半正弦波、三角波和梯形波激发下的全波形解析解、VTEM实际激发波形的后推欧拉算法计算结果,检验了本研究开发的数值模拟算法的正确性。设计地表起伏异常体模型,计算和分析了航空瞬变电磁响应特征。开发的基于位移逆Krylov子空间的全波形航空瞬变电磁法三维数值模拟算法适合模拟复杂地电模型的响应,具有较高的计算精度。 相似文献
16.
17.
18.
采用有限差分方法,对水平电偶极子源在海洋二维介质情况下电磁响应进行数值模拟。利用油气高阻异常体模型的数值模拟结果,研究MVO(Magnitude Versus Offset)曲线解释方法的有效性,进而分析高阻油气异常体的水平位置的变化,对电场和磁场幅值的MVO曲线产生的影响,得出了通过MVO曲线定性分析高阻油气异常体相对水平位置的方法。本研究可为海洋可控源电磁法数据的定性分析提供部份理论依据。 相似文献
19.
殷文 《吉林大学学报(地球科学版)》2008,38(1):144-151
在弹性波频率空间域有限差分数值模拟方面,差分网格及边界条件是影响弹性波模拟成功与否的关键,为了压制数值模拟中的网格频散,采用25点有限差分算子,建立了有限差分矩阵方程,且借鉴匹配层衰减边界条件思想,设计了弹性波频率空间域有限差分数值模拟算法。由于采用高阶有限差分法来提高差分格式的精度,将会导致计算量显著增加,为此,对频率空间域有限差分弹性波数值模拟方法,采用流水线技术与分治策略进行了并行算法研究,提高了计算效率,使得在合理的计算时间内更精确地模拟弹性波在弹性介质中的传播过程。 相似文献
20.
从回线源瞬变电磁法三维有限差分正演算法出发,与均匀半空间解析解作比较,验证了数值模拟算法的正确性。利用该正演算法分别研究了巷道、背景地层电阻率、掘进机对矿井瞬变电磁超前探测响应的影响。结果表明,背景地层电阻率越低,响应幅值越强,巷道影响越显著;掘进机距探测线框超过10m时,影响幅度明显变小;矿井瞬变电磁法的方向性具体表现为在不同的探测角度上测量的异常幅度不同。通过研究这些影响因素,对分析矿井瞬变电磁超前探测资料具有重要的借鉴作用。 相似文献