共查询到18条相似文献，搜索用时 593 毫秒
1.
利用线性旅行时插值射线追踪 计算近地表模型初至波走时 总被引：4，自引：0，他引：4
利用线性旅行时插值射线追踪对近地表模型进行正演计算,可以快捷、准确地获得初至波走时和射线路径。由于该算法计算的初至波不局限于折射波,因此很好地解决了浅层折射勘探中的低速"隐蔽层"问题;而且,由于该算法是基于网格划分和线性插值,因此它不仅可以追踪任意复杂介质的初至波,而且可以使得追踪的初至波射线路径逼于真实,避免了同类算法直接连接网络节点形成射线路径的缺陷(路径过于弯折,计算走时偏大)。将LTI算法同其他几种算法的追踪结果进行的对比和分析表明,LTI算法在计算初至波走时和射线路径方面较其他算法更为精准、稳定,是一种有效的射线追踪方法。 相似文献
2.
射线追踪是一种重要的地震数据处理工具,而地震波旅行时的计算是射线理论的核心问题。平面波近似旅行时计算方法精度低,且对于介质速度的适用性弱,因此本文在精度较高的球面波近似旅行时计算方法的基础上,对射线追踪算法提出了改进。数值模拟结果表明了该算法的正确性,能够有效地处理速度畸变界面,在不影响计算精度的情况下提高计算效率。 相似文献
3.
旅行时射线追踪的精度和效率一直是影响地震层析成像质量和效果的关键因素。这里在三维近地表速度结构层析成像应用中,结合旅行时插值算法与最短路径算法,提出了一种三维初至波射线追踪算法。通过对双线性旅行时插值算法进行改进,并运用判定条件与简化插值计算公式进行快速计算,减少了插值次数,降低了运算量,同时也保证了较高的射线精度,有效地解决了三维射线追踪算法的计算低效率,射线精度不高的问题。 相似文献
4.
利用格林公式及点源场的边界条件,将地下半空间位函数的体积积分变为地面位函教的积分方程以后,用边界元方法解此积分方程便可得到地面各节点的电位值。其计算精度和速度取决于剖分方式和节点密度。本文采用变网格的剖分方式,用较少的点实现了计算,并兼顾了精度和速度。 相似文献
5.
《物探与化探》2017,(5)
随着浅层地裂缝等地质灾害问题精细探测需求的增大,采用初至波速度层析反演技术,可以提供较高精度的浅层速度场的纵横向异常变化,可为识别地裂缝等地质灾害现象提供依据。通常,在采集参数一定的条件下,其反演成像精度受算法、网格类型和剖分尺度的限制,在网格类型不变的条件下,要提高精度就必须加密正反演计算网格,改进算法,从而实现增量而不减速且高精度的效果。文中利用OpenMP的单机多核并行技术,探讨并实现了初至波层析成像并行算法,其中正演部分使用了改进的旅行时线性插值法,针对原有的按列计算,加入了按行计算,充分考虑到来自各方向的射线,使得计算出的旅行时最小,反演部分使用了能迅速稳定收敛的LSQR法以及正则化技术,通过理论和实际资料测试,其收敛速度快,反演结果较好。同时,在加密采样和缩小网格尺度的条件下,与常规串行算法相比,其运算速度和效率有较大提高,与商业软件比较,其运算效率差异不大,但反演结果的精度和可分辨性较好。 相似文献
6.
提出了一种新的射线追踪方法——模拟退火法.新方法不仅较好地解决了传统射线追踪方法过分依赖初始模型的问题,而且对于复杂速度场模型也能保证在固定的发射与接收点之间构建令人满意的射线路径及其相应的走时,搜索到满足旅行时全局最小的射线路径.此外,新方法还可计算局部最小旅行时,并可方便地通过指定射线经过固定点来对射线路径进行限制... 相似文献
7.
针对混凝土质量检测的特点,采用弯曲射线进行超声层析成像。首先建立波速分布的网格模型,并引入最短路径法进行射线追踪,然后采用联合迭代重建技术（SIRT）反演成像。为使算法适应弯曲射线,对SIRT的计算步骤进行调整,增加外层迭代,以便对修正后的路径进行下一轮反演;为提高反演精度,根据射线传播距离、幅度衰减、频率突变等信息建立和调整信号可信度系数,并依据该系数在求取网格慢度修正值时进行加权平均;为缩短计算时间,动态调整网格节点个数,并进行双向最短路径追踪。经过数值模拟及对混凝土基桩的实测,证明该方法在计算准确性和精度上有明显改进,计算速度也可以满足工程应用的需要。 相似文献
8.
扩展有限元法的网格剖分不需要考虑裂纹的位置,但计算时需根据计算网格和裂纹位置关系确定加强节点和加强类型。一般根据加强单元类型确定加强节点和加强类型。确定加强单元类型,不仅要考虑裂纹位置,还要考虑单元形状。指出现有文献中只考虑裂纹位置根据节点水平集值判别加强单元类型的不足,提出相应的改进方法。对于复杂形状裂纹,采用水平集法确定加强单元类型并不方便;给出了一种确定加强单元类型、加强节点和加强类型的有效方法,即根据裂纹与单元边界交点数和裂纹拐点位置确定加强单元类型,然后根据加强单元类型确定加强节点及加强类型;给出了该方法详细的执行过程。算例分析表明了本文方法的正确性。 相似文献
9.
在利用有限差分等基于网格的数值分析方法解地震波走时所满足的程函方程时,由于速度模型的网格化离散等原因,会使走时在各网格节点之间不具有计算射线路径所要求的光滑性,即走时在邻近网格节点之间不具有连续的一阶导数。因此,直接利用网格节点走时计算射线路径会使最终的射线路径不光滑。为解决这个问题,已有研究者提出了基于B样条插值的逆向梯度方案（法）。然而,在速度发生突变时,B样条逆向梯度法所计算出的射线路径会具有较大的误差。针对这个问题,首先采用适合于解最小零偏差逼近及最佳平方逼近问题的Chebyshev多项式取代B样条对来自于分区多级计算方案的网格节点走时进行最佳逼近,得到在最小平方意义下的最优走时公式;然后采用与B样条逆向梯度法类似的计算过程得到光滑的射线路径。数值实验表明,利用Chebyshev多项式逼近走时可以得到具有很高精度的多次反射射线路径,在多次波偏移成像研究中具有潜在的价值。 相似文献
10.
因对初至旅行时层数的确定和拐点拾取存存不确定性．导致时间项反演方法在近地表复杂区建立的速度模型精度不高,静校正效果不理想。对此提出了利刚旅行时差值曲线方法进行层数确定和拐点拾取。理论分析认为,运用旅行时差值曲线方法时,应尽餐选择相邻的激发点。通过对三层近地表速度模型进行射线追踪。得到初至旅行时,根据其旅行时曲线初步判定层数,再运用旅行时差值曲线进一步确定层数。观察时间项反演速度模型,其与近地表速度模型较接近,经计算反演速度模型与设计的近地表速度模型之间的均方根误差仅为0．95m,可见该方法可有效提高时间项反演的精度。 相似文献
11.
地震波走时计算的逆风差分算法 总被引：5，自引：0，他引：5
快速精确地旅行时计算在地震资料的叠前偏移与层析反演中起重要作用。利用逆足差分格式求解程函方程，首先在波阵面上寻找全局极小点，然后计算该点周围各点的地震波初至走时，同时，采用局部算法以确保计算的稳定性。该算法对于任意复杂地质模型都具有较高的精度。 相似文献
12.
波前扩展有限差分地震波走时算法具有物理意义明确、因果稳定性强的特点,但每次波前扩展都要寻找波前面上的最小走时点。当计算网格点数较多,特别是涉及到三维走时计算时,寻找波前面上的最小走时点是一项十分耗时的工作。研究发现,波前扩展有限差分地震波走时算法的波前点具有两个突出特点:①波前点更新十分频繁,通常每次取出波前最小走时点后都要插入若干新的波前点;②新计算出的波前点的走时通常比较大。数据结构中的二叉树堆排序方法可以提高寻找波前面上最小走时点的效率,根据特点①,在原始二叉树堆排序方法的基础上,优化了插入新波前点和移除波前面上最小走时点的流程,实际计算结果表明,与原始的二叉树堆排序方法相比,改进后的二叉树堆排序方法可以提高大约20%的计算效率。根据特点②,将原始的二叉树堆排序方法推广到多叉树,实际计算结果表明,完全三叉树堆排序方法优于二叉树和四叉树堆排序方法,可以再提高5%的计算效率。 相似文献
13.
结合GIS(geographic information systems,地理信息系统)的空间数据处理能力,将潘家铮法从二维扩展到三维,提出一个基于栅格柱体单元的库岸滑坡滑速计算方法。建立了基于栅格柱体单元的三维滑坡滑速计算分析模型,给出了模型参数在GIS中的空间计算表达式。其次,通过对计算模型进行受力分析,结合牛顿运动定律,选取了沿滑体滑动方向和垂直方向来建立动力平衡方程,从而求解三维滑坡滑速。在GIS中开发出一个三维滑坡滑速计算扩展模块,并通过实例与潘家铮法所得结果进行对比。结果表明,该方法比潘家铮法计算的最大速度结果大15.2%,比启动加速度结果大32.8%,比达到最大速度的滑动时间短1 s。 相似文献
14.
以重力位在场源内部满足泊松方程为依据,以重力矢量满足第三类边界条件为切入点,推导了与三度体重力矢量满足的边值问题相对应的变分问题,进而利用有限单元法实现了对变分问题的求解.立方体模型试验结果表明:文中提出的新的系数矩阵存储方式较之传统方式能够更有效地节约存储空间,且为利用预条件共轭梯度技术更加快速地求解线性方程组提供了保障;重力矢量的计算精度与边界长度及单元网格的边长息息相关,其计算效率则主要取决于所要计算的节点总数和大型稀疏线性方程组求解算法的优劣;一般情况下,当单元的边长小于场源体边长的1/10、边界长度大于场源体长度的7.5倍时,能够获得理想的结果. 相似文献
15.
线源二维时间域电磁响应的有限差分解 总被引：1，自引：0，他引：1
本文以ＫｕｏａｎｄＣｈｏ（１９８０）和Ｏｒｉｓｔａｇｌｉｏ ａｎｄ Ｈｏｈｍａｎｎ（１９８３）工作基础，对时间域中直接求解无限长电流源激发下任意二维电阻率分布的麦克期韦方程的有限差分方法进行了重新研究。提出了一种改进方法，这种方法的独特的上在于用变化方向法把二维问题化为一系列一维问题来求解，并保留了前人在地－空边界使用地表向 拓的边界条件以避免求空气中的场的处理技术，从而节省了计算机内存，简化了计 相似文献
16.
Richards方程在非饱和渗流模拟及其他相关领域应用广泛。在数值求解过程中,可以采用有限差分方法进行数值离散并迭代求解,为了获得较可靠的数值解,常规的均匀网格空间步长往往是较小的。在一些不利数值条件下,如入渗于干燥土壤,迭代计算费时甚至精度也不能得到很好改善。因此,文章提出Chebyshev空间网格改进方法,结合有限差分方法对Richards方程进行数值离散以获得线性方程组,并通过经典的Picard迭代方法进行迭代求解线性方程组以得到Richards方程的数值解。通过均质土和分层土2个不利情况下的非饱和渗流算例,又结合模型解析解和软件Hydrus-1D,对比研究了改进网格方法与均匀网格方法获得数值解的精度。结果表明,提出的Chebyshev网格方法相较于传统的均匀网格,可以在较少的节点数下获得较高的数值精度,又具有较小的计算开销,有较好的应用前景。 相似文献
17.
地震波旅行时计算是射线追踪方法的核心问题,其精度直接影响着路径的精度和正演模拟的准确程度。在简要分析目前最典型的几种地震波旅行时计算方法(如有限差分法、线性插值法等)局限性的基础上,提出了一种基于非线性插值的地震波旅行时计算方法,并推导了其计算公式。通过模型的对比试算,不仅验证了所提出的旅行时非线性插值算法的正确性,而且还表明本算法计算旅行时的精度,比Asakawa[5]提出的线性插值算法的精度更高。 相似文献
18.
通常情况下,含水层的参数可分别借助于Dupuit和Theis井流公式,在局部地段利用稳定流和非稳定流抽水试验获得.在新井设计阶段,一般是依据现有参照井数据进行线性或非线性插值计算,但插值受插值表的约束,效率不高.本次研究采用2次优化方法,用网格法确定所求参数最优估计值作为下一步的叠代初值.用共轭梯度法逐步逼近所求参数的精确解,高效地利用参照井的基本数据得到了精度较高的新井参数,为准确地设计新井提供可靠的依据,实地应用效果良好. 相似文献