共查询到20条相似文献,搜索用时 62 毫秒
1.
2.
余涛朱自强鲁光银曹书锦 《物探化探计算技术》2013,(4):446-451
重力张量是重力位的二阶空间导数,对密度体的变化和细节部份反映更为灵敏。对于复杂的重力密度模型来说,张量的解析公式是很难推导的。为了模拟这种情况,将复杂模型进行有限元剖分,计算每个单元的重力异常对测点张量的影响,最后叠加得到整个复杂模型在测点处的重力张量。采用Delaunay四面体非结构化网格对密度体进行剖分,经分析表明,这种剖分方式具有较高的精度。有限元剖分的方式拟合复杂形体缺点之一是计算量比较大,而并行化能有效解决这个问题。这里基于OpenMP并行模型编写了重力张量并行正演程序,分析了不同情况下并行的执行性能,为大规模复杂模型的重力张量正演提供了一种并行策略和思路。 相似文献
3.
重力张量是重力位的二阶空间导数,对密度体的变化和细节部份反映更为灵敏.对于复杂的重力密度模型来说,张量的解析公式是很难推导的.为了模拟这种情况,将复杂模型进行有限元剖分,计算每个单元的重力异常对测点张量的影响,最后叠加得到整个复杂模型在测点处的重力张量.采用Delaunay四面体非结构化网格对密度体进行剖分,经分析表明,这种剖分方式具有较高的精度.有限元剖分的方式拟合复杂形体缺点之一是计算量比较大,而并行化能有效解决这个问题.这里基于OpenMP并行模型编写了重力张量并行正演程序,分析了不同情况下并行的执行性能,为大规模复杂模型的重力张量正演提供了一种并行策略和思路. 相似文献
4.
在有限单元法矩形~三角形单元剖分的基础上,提出了适合复杂地形条件下的网格剖分方法:平行四边形~三角形剖分方法。计算结果表明,该方法具有快速、高效的特点。只要对原有程序稍加修改,即可实现对复杂地形条件下的复杂地电断面进行自动剖分和网格节点自动编号。因此,对于复杂地形条件下的地电模型,该方法能快速地实现较为精确的数值模拟。 相似文献
5.
《物探化探计算技术》2019,(5)
有限元求解大地电磁正演问题时,对研究区域的剖分常规做法是采用规则化的网格。但规则化网格在剖分地形、断层、褶皱等复杂模型会产生较大的几何离散误差。针对上述情况,这里采用非结构化的四边形网格对二维地电模型进行网格剖分,并与自适应有限元相结合,由剖分的粗网格出发,利用每个频点下网格单元的后验误差估计值指导网格的局部加密,优化网格质量和数量,从而提高正演模拟的精度。通过一维K型地电模型利用本文方法算出数值解与解析解进行对比分析,验证了自适应有限元法在求解大地电磁二维正演模拟中的有效性。并通过对断层、褶皱模型的正演模拟,分析了其大地电磁正演响应特征。 相似文献
6.
为了更好地模拟地下介质连续变化及开展连续介质的反演,对二维电导率分块线性变化的线源频率域可控源电磁法进行了有限元正演模拟,在剖分单元内同时对电场及电导率参数线性插值,使电导率参数在剖分单元之间保持连续变化。首先,提出有限元正演模拟的边值问题及变分问题,并详细论述了有限元的剖分、插值、单元分析及总体合成的各个步骤;其次,采用稀疏存储及基于不完全LU分解的BICGSTAB算法求解复系数方程组,节省了内存并提高了计算速度;然后,对一个均匀半空间模型进行模拟,计算结果表明,低频及高频的有限元数值解都与解析解吻合,证明了算法的正确性;最后,对水平层状模型及垂直断层模型进行正演计算,视电阻率及相位的等值线图均较好地反映出了异常体,说明文中算法能够对电导率连续变化的线源可控源电磁法进行有效地模拟。 相似文献
7.
8.
传统的MODFLOW采用结构化网格进行剖分,但是对刻画形状复杂的边界等实际问题存在局限性。最新的非结构化网格MODFLOW(MODFLOW-USG)能够提供更加灵活的网格剖分方法,并能够刻画形状不规则的边界、非均质地质体和不连续层等实际问题,提高模型仿真度与计算结果精度。利用MODFLOW-USG方法对黑河中游盆地地下水流动模型进行了研究,针对形状复杂、尺度不一的河流边界进行了精细化剖分,提高了河流及附近重点关注区模拟结果的准确性。 相似文献
9.
网格剖分是数值模拟中有限元方法的基础和关键之一 ,剖分的好坏直接影响模拟结果的精度。数值模拟目前已应用于许多领域 ,如地球化学、水文地质、工程地质、矿床地质、油气地质、油藏工程、土力学和岩石力学等。在众多的网格剖分中 ,三角剖分最具优势 ,可以很好的拟合任意形状的边界 ;三角平面的组合可以较好地逼近任意曲面 ;三角形元素的局部逼近多项式可以采用完全多项式。以往三角剖分多采用手工剖分或者设计复杂的自动剖分 ,具有速度慢、扩大误差、增加计算周期和不能充分拟合边界等缺点。笔者在工作中发现了两种三角剖分自动生成方法—… 相似文献
10.
《物探化探计算技术》2017,(6)
为满足电磁法二维数值模拟解释对交互建模与可视化网格剖分的需求,利用计算机图形学、人机交互、拓扑关系学等技术,设计了基于测线、测点等信息的二维矢量建模、四边形和三角形网格剖分方法与流程,并开发形成了软件模块。该方法与软件可提供起伏地形下大地电磁、可控源音频大地电磁法二维数值模拟的三角形、四边形网格(下边界随地形变化或平地形)剖分实用化工具。将该方法应用到多种模型试验中,取得了良好的建模与网格剖分效果。 相似文献
11.
《物探与化探》2015,(3)
介绍了有限差分算法在大地电磁测深法中的应用,推导了一维层状介质和连续介质模型下的有限差分算法,提出了一种新的网格剖分方法,并通过两层和三层介质模型的正演计算,验证了算法的正确性以及网格剖分方法的合理性。设计了一种连续介质模型,与两种不同程度近似的多层层状介质模型进行了正演结果的对比,指出了研究连续介质模型正演的必要性。对于连续介质模型,计算中剖分得越详细,其结果越接近真实值,而随着误差的逐渐减小,精确度的提高不再明显,但资源与时间的消耗将大幅增加。因此,无限制地减小网格步长是不必要的。同时,为了兼顾误差大小与网格步长的关系,采用了一种将一次网格剖分进行二次差分计算的方法得到新的结果,使计算结果明显得到了改善。 相似文献
12.
研究了不同测量方式(井-地,地-井,井-井)下点源场井中电法的三维有限元数值模拟。考虑到深度方向上大范围的网格剖分和井眼的影响及井-井测量等因素,采用放射状三棱柱单元的网格剖分方式,以提高网格质量,减少剖分单元数;给出了三棱柱单元的坐标变换公式,进行精确的单元积分,减少了单元积分时间;结合非结构化网格技术,实现了复杂模型的模拟;开发出相应的程序实现了复杂条件下(如考虑井眼影响、井井测量、倾斜井情形、地形起伏等)电法测井的三维有限元模拟,数值算例验证了方法的可靠性及计算效率,并对不同情形下的异常响应进行了分析,为进一步的反演工作奠定了基础。 相似文献
13.
14.
利用两阶段建模思路建立条件随机模拟模型。第一阶段在地质体表面利用Denaunay三角网剖分算法实现三角网网格化,第二阶段用Denaunay四面体剖分算实现地质体的四面体网格化,并对四面体单元进行三角面插值,进而得到三角网,最后用OpenGL实现地质体三维模型的可视化。 相似文献
15.
以双异性点电流源电流场二维地电断面为例,着重讨论了电阻率展析成像反演中有限单元法数值模拟边界条件的选择、网格剖分以及模拟精度等问题,同时也讨论了成像反演中最优化及演方法及拟合参数的评价等。文中提出的新的网格剖分方法和选择的边界条件可有效地提高计算速度和数值模拟的精度。 相似文献
16.
针对地下水有限元数值模拟中区域三角网格剖分复杂难以处理的情况,提出适合其特点的Delaunay三角网格自动剖分方法,并对含有多个参数分区的含水层进行网格剖分,最后利用遗传算法和有限元程序相耦合来反演含水层水文地质参数。结果表明此方法可大大简化地下水数值模拟的前处理工作,并能提高有限元网格剖分的有效性和准确性,从而得到令人满意的数值模拟结果。 相似文献
17.
18.
多极距中梯测深法是一种效率高、信息丰富的新型体积测深技术,而当前常规多极距中梯测深应用与研究多数局限于二维空间,其三维正反演一直是人们研究的热点。本文在四面体交错网格剖分基础上,采用有限单元法对球体等三维模型进行多极距中梯测深数值模拟计算。研究表明,三维数值模拟可有效反映模型的空间分布;在不增加太多工作量的情况下,实现高效体积测深;电阻率和极化率参数模拟相结合分析,模拟精度高,稳定性好;开展倾斜铜板、铜块水槽物理模拟,模拟结果较好反映了模型的空间及电性分布,验证了多极距中梯测深的可行性和正确性,为该方法的实际应用及后续的反演工作提供了重要依据。 相似文献
19.
20.
有限元法是求解地下水流和溶质运移对流-弥散方程的常用数值方法,它可以精确高效地处理以弥散为主的问题,但求解以对流为主的问题易引起显著的数值振荡。通过Galerkin有限元法对变异Henry问题进行模拟求解,得到了用不同的剖分网格及水动力弥散系数时,在特选节点处的浓度穿透曲线,分析并找到了浓度振荡的原因及合适的消除方法,即若出现浓度数值解在某值附近振荡,可以通过加密网格或增加水动力弥散系数将其消除。模拟结果及其分析表明:即使是研究区域相同,不同的边界条件、不同的水动力弥散系数对网格精度的要求不同;换言之,同一网格对不同模型参数的有效性也不同。网格Peclet数能够有效地判定给定的网格剖分是否会引起浓度振荡,对有限元法数值计算的网格剖分具有指导意义。 相似文献