首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 421 毫秒
1.
三维地震波动方程有限差分数值模拟,受到计算量大且内存消耗大地制约而计算效率较低,严重影响其发展和应用。目前,一种全新的异构众核协处理器,为三维并行计算提供了一条解决计算效率和数据存储问题的重要途径。基于三维有限差分计算的特点,介绍了面向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.
朱成  李桐林  范翠松 《世界地质》2014,33(4):916-924
在可控源电磁场三维有限差分数值模拟过程中,低频时数值模拟的收敛速度变慢,导致获得精确结果的时间变长。引入电流密度的散度校正方法可以加快低频数值模拟时的收敛速度,提高解的精确性。笔者在研究可控源电磁场三维有限差分数值模拟算法的基础上,利用电流密度的散度为零,推导了散度校正公式并利用迭代解法实现了散度校正算法。模型试验表明,在低频时,散度校正能大幅度提高三维电磁场数值模拟的速度和精度,异常体规模越大,数值模拟的频率越低,散度校正的效果越明显。  相似文献   

6.
基于球坐标系下有限差分的地磁测深三维正演   总被引:2,自引:0,他引:2  
为了计算全球尺度电磁感应的响应,本文介绍地磁测深频率域三维正演。正演算法采用球坐标系下的交错网格有限差分方法,从Maxwell方程的积分形式出发,采用PARDISO对离散后的方程组求解,避免了迭代求解的散度校正。为了验证本文结果的正确性和精度,与前人的有限元和有限差分方法进行了对比,一维层状模型的三维交错网格有限差分数值结果和解析解相对误差小于5%,双半球模型的计算结果与前人的计算结果完全吻合。三维"棋盘模型"计算表明磁场分量对异常体的大小和位置具有很好的分辨能力。  相似文献   

7.
三维定位原理与三维反射波场模拟   总被引:1,自引:0,他引:1  
在三维正演模拟定位原理基础上 ,以单程声波方程作为波场外推手段 ,实现三维叠前反射波场模拟 ,完全克服了全波方程模拟的 3个问题。在波场外推中 ,使用波动方程的混合解法 ,即以波动方程解析解 (相位移法 )为主 ,在介质速度横向变化的地方 ,用数值解 (45°有限差分 )作补充 ,使得整个算法精确、快速 ,既能处理复杂构造 ,又能适应速度的任意变化。  相似文献   

8.
介质的各向异性会引起大地电磁观测曲线的畸变,对大地电磁数据解释影响很大。这里从麦克斯韦方程组出发,推导了基于二阶插值基函数的二维大地电磁非结构有限元公式,实现了任意地形起伏条件下任意各向异性二维大地电磁响应的数值模拟。各向异性断层模型的数值模拟结果与解析解吻合一致,复杂各向异性模型的有限元模拟与已有的有限差分计算结果也相符合,验证了算法的正确性。在此基础上探讨了各向异性参数对二维大地电磁响应的影响。  相似文献   

9.
波动方程有限差分正演模拟误差来源分析   总被引:2,自引:0,他引:2  
有限差分波动方程正演数值模拟主要是解决差分格式、震源选择、边界条件处理等问题.这些问题的合理解决对得到正确和准确的波动方程数值解有着十分重要的作用.笔者通过对波动方程计算机有限差分正演模拟过程的叙述,分析了在做正演模拟时的误差来源以及如何有效地减小误差.  相似文献   

10.
利用交错网格有限差分技术实现了频率域海洋可控源电磁三维正演算法,并探讨了海底地形对海洋可控源电磁响应的影响和校正方法及效果。通过与海底二维山峰地形的二维自适应有限元解对比,验证了文中算法的正确性和模拟海底起伏地形的可行性,并分析了海底二维山峰地形对三维油气储层电场响应的影响。关于海底起伏地形对海洋可控源电场响应的影响,采用比较法对其进行了地形校正,结果表明,在深水区比较法可以很好地消除地形影响,恢复海底高阻油气储层的电场响应。  相似文献   

11.
介质含水率与探地雷达信号关系数值模拟   总被引:1,自引:0,他引:1  
采用时间域有限差分技术实现了探地雷达三维正演计算,并以此对介质含水率与探地雷达信号关系开展了数值模拟研究.物理模型实验与数值模拟计算结果的对比分析证明了数值模拟的有效性.  相似文献   

12.
分布荷载推力桩计算的p-y曲线法研究   总被引:2,自引:2,他引:0  
为提高现有水平梯形分布荷载推力桩设计计算水平,提出了地基系数按非线性土抗力.水平位移(p-y)曲线抗力模式表达的水平梯形分布荷载推力桩位移和内力计算的有限差分数值分析方法,并详细推导了桩身位移的差分格式.基于这些公式,编制了全桩位移、内力及侧土抗力的计算和图形处理程序,可适用于滑坡抗滑桩和深基坑悬臂支护桩的设计计算.结合新型桩型现浇混凝土薄壁管桩(PCC桩)作为滑坡抗滑桩进行对比分析,算例表明:该方法方便可靠.当有限差分段划分得足够小时,可使数值解接近于真实解.  相似文献   

13.
ZTEM(Z轴倾子电磁法)是一种天然场源的频率域航空电磁法,其特点是磁场垂直分量在空中机载平台测量,磁场水平分量在地面的固定基站测量,具有勘探深度大、速度快、成本低、覆盖面积大等技术优势。本文实现了ZTEM二维有限差分正演和二维非线性共轭梯度(NLCG)反演算法。研究对象是倾子资料,反演过程中通过解"拟正演"问题来避开雅克比矩阵的直接计算。通过理论模型合成数据反演试算,验证了ZTEM倾子资料二维NLCG反演算法的稳定性与可靠性。与大地电磁(MT) TE模式阻抗资料反演结果进行对比,发现在异常体横向边界的约束方面,ZTEM倾子反演比MT阻抗反演更具优越性。  相似文献   

14.
研究了电法测井不同测量方式(井—地,地—井,井—井)下点源场井中电法的三维有限差分数值模拟。采用六面体网格剖分方式来对模型进行剖分,运用一维非零元素行压缩存储模式来存储系数矩阵,减少了内存需求和计算量;采用不完全Cholesky共轭梯度(ICCG)方法来求解线性方程组,提高了求解效率;编制了相应的程序实现了井—地、地—井、井—井、倾斜井条件下的电法测井三维有限差分数值模拟。设计的算例结果验证了该算法的正确性和效率性,并且分析了各种情形下的异常特征,为进一步的反演工作打下了基础。  相似文献   

15.
董岩  谭捍东  付兴 《现代地质》2023,37(1):74-83
为了研究复杂地电模型的航空瞬变电磁法全波形响应特征,需要开发考虑发射波形的三维数值模拟算法。本研究基于非结构四面体网格和位移逆Krylov子空间(Shift-and-Invert Krylov,简称SAI Krylov)方法,采用基于电偶极子离散的场源处理方法模拟场源,在时间域进行计算实现了全波形航空瞬变电磁法矢量有限元三维数值模拟。使用均匀半空间模型在阶跃波、半正弦波、三角波和梯形波激发下的全波形解析解、VTEM实际激发波形的后推欧拉算法计算结果,检验了本研究开发的数值模拟算法的正确性。设计地表起伏异常体模型,计算和分析了航空瞬变电磁响应特征。开发的基于位移逆Krylov子空间的全波形航空瞬变电磁法三维数值模拟算法适合模拟复杂地电模型的响应,具有较高的计算精度。  相似文献   

16.
以巷道底板供电为例,采用三维有限差分数值模拟技术对矿井下有巷道存在时的全空间点源电流场的分布进行了计算,根据计算结果初步总结了巷道的影响规律.  相似文献   

17.
以体积分方程法模拟三维电磁正演问题的原理为基础,利用数值滤波方法实现了时间域水平长导线源的瞬变电磁三维问题的模拟,对电阻率、极化率、时间常数和频率相关系数的变化分别进行了三维数值模拟。为了提高计算的精度,变换过程中采用双精度数值滤波方法。利用偶极源的解析解和数值解的对比来检验数值滤波方法的可靠性。数值模拟结果表明,此种数值滤波方法精度高,结果稳定。高精度数值滤波结果为瞬变电磁法在石油、矿产、地热以及其他方面的应用提供了理论指导。  相似文献   

18.
采用有限差分方法,对水平电偶极子源在海洋二维介质情况下电磁响应进行数值模拟。利用油气高阻异常体模型的数值模拟结果,研究MVO(Magnitude Versus Offset)曲线解释方法的有效性,进而分析高阻油气异常体的水平位置的变化,对电场和磁场幅值的MVO曲线产生的影响,得出了通过MVO曲线定性分析高阻油气异常体相对水平位置的方法。本研究可为海洋可控源电磁法数据的定性分析提供部份理论依据。  相似文献   

19.
在弹性波频率空间域有限差分数值模拟方面,差分网格及边界条件是影响弹性波模拟成功与否的关键,为了压制数值模拟中的网格频散,采用25点有限差分算子,建立了有限差分矩阵方程,且借鉴匹配层衰减边界条件思想,设计了弹性波频率空间域有限差分数值模拟算法。由于采用高阶有限差分法来提高差分格式的精度,将会导致计算量显著增加,为此,对频率空间域有限差分弹性波数值模拟方法,采用流水线技术与分治策略进行了并行算法研究,提高了计算效率,使得在合理的计算时间内更精确地模拟弹性波在弹性介质中的传播过程。  相似文献   

20.
从回线源瞬变电磁法三维有限差分正演算法出发,与均匀半空间解析解作比较,验证了数值模拟算法的正确性。利用该正演算法分别研究了巷道、背景地层电阻率、掘进机对矿井瞬变电磁超前探测响应的影响。结果表明,背景地层电阻率越低,响应幅值越强,巷道影响越显著;掘进机距探测线框超过10m时,影响幅度明显变小;矿井瞬变电磁法的方向性具体表现为在不同的探测角度上测量的异常幅度不同。通过研究这些影响因素,对分析矿井瞬变电磁超前探测资料具有重要的借鉴作用。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号