共查询到19条相似文献,搜索用时 156 毫秒
1.
三维层状地质结构下整体迭代射线追踪法 总被引:2,自引:1,他引:2
作者在文中对三维任意界面情况下的两点间射线追踪问题进行了研究,从Fermat原理出发,对射线路径上的旅行时,采用一阶Taylor不完全展开等办法,推出了三维情况下计算射线路径的整体迭代算法。并通过模型计算证明了本方法的稳健性以及计算的高精确度和高效率。为三维结构下地震偏移和层析成像提供了一种可行的射线追踪方法。 相似文献
2.
作者在文中对三维任意界面情况下的两点间射线追踪问题进行了研究,从Fermat原理出发,对射线路径上的旅行时,采用一阶Taylor不完全展开等办法,推出了三维情况下计算射线路径的整体迭代算法。并通过模型计算证明了本方法的稳健性以及计算的高精度和高效率,为三维结构下地震偏移和层析成像提供了一种可行的射线追踪方法。 相似文献
3.
基于F erm at原理和Sne ll定理对三维任意界面情况下的两点间GPR波传播射线追踪问题进行了研究,采用了一个适用于任意界面情况下计算反(折)射点的一阶近似公式,从任意给定的初始路径出发进行逐段迭代计算,给出了三维逐段迭代射线追踪算法的计算格式,并进行了三维介质下射线追踪计算。结果表明计算速度相当快,且其计算精度可以根据需要达到任意要求。在此基础上,根据雷达波传播规律,考虑每一段的能量损失。根据地质雷达的发射子波,最终获得三维介质中雷达波传播和接收记录。并将结果应用于实际剖面拟合中,提高了解释的精度。 相似文献
4.
二维声波方程稳健迭代速度反演 总被引:1,自引:0,他引:1
从二维声波方程初、边值问题出发,通过引进一个广义幂指数误差分布函数,以及对波动方程的Lippman-Schwinger方程的解进行Bon近似,导出一个残差加权迭代最小二乘(稳健迭代)算法,进而实现对地下介质剖面的速度结构以演计算。模拟计算表明,本方法具有较高的精度和较强的抗干扰能力,是求解声波方程反演问题的一种有效方法,同时也为地下介质的速度分布成像提供一种新的技术。 相似文献
5.
6.
7.
首次提出了用于解决塔里木盆地地震资料数据处理中浅部地震速度分析的折射层析成象方法。该方法可以通过求解任意变速介质中的二维程函方程、射线追踪和SIRT迭代解法实现。通过对塔里木盆地地震资料的层析反演取得了较好的结果,揭示了浅部地震速度具有较大的横向变化。这一变化对建立地下速度模型及反射构造成象有较大影响,钻井资料证实了这一结果。这项研究的成功对研究塔里木盆地地震波速度的变化规律及进一步提高该区整体地震资料处理与解释的质量具有重要的实际意义。 相似文献
8.
二维介质速度成像的新方法是基于地震波的旅行时间,由射线路径沿线的速度分布特征所决定的,将二维介质中的速度分布用两个分别以二维空间坐标为自变量的勒让德多项式的线性组合来描述。根据走时积分方程,即可建立地震波走时与其射线路径沿线的速度分布之间的对应关系,进而重建出地震波射线覆盖区域的速度结构。该成像技术是一个迭代过程,且适用于折射波及反射波。 相似文献
9.
针对混凝土质量检测的特点,采用弯曲射线进行超声层析成像。首先建立波速分布的网格模型,并引入最短路径法进行射线追踪,然后采用联合迭代重建技术(SIRT)反演成像。为使算法适应弯曲射线,对SIRT的计算步骤进行调整,增加外层迭代,以便对修正后的路径进行下一轮反演;为提高反演精度,根据射线传播距离、幅度衰减、频率突变等信息建立和调整信号可信度系数,并依据该系数在求取网格慢度修正值时进行加权平均;为缩短计算时间,动态调整网格节点个数,并进行双向最短路径追踪。经过数值模拟及对混凝土基桩的实测,证明该方法在计算准确性和精度上有明显改进,计算速度也可以满足工程应用的需要。 相似文献
10.
论述了用数字方法设计出地层弹性波速度模型的原理和方法,阐述了地震射线追踪法的算法并用其对数字模型的反射、折射波时距曲线做了计算,通过实例分析了复杂速度结构与地震剖面的对应关系,可指导解释人员根据实际地质情况任意设计地层模型进行计算,很快得出结果。 相似文献
11.
12.
为了实现适应崎岖海底、大陡坡、海底火山等复杂海底地质条件且灵活、稳定、精度高的射线追踪方法,并基于该方法详细分析复杂海底条件下各种地震波型的运动学特征,综合多种算法的优势,实现了一种快速推进迎风插值射线追踪方法。首先,采用混合网格法剖分复杂海底地质模型;其次,通过融入迎风差分思想的线性插值策略来构建精度高、无条件稳定且灵活的局部走时和射线路径计算公式;然后,综合应用这些公式和多级次快速推进法,灵活计算整个模型各种类型地震波的走时,并基于逆向追踪方法计算射线路径;最后,对该方法进行了精度分析,发现其能够获得相对高的走时和射线路径计算精度,且反射波的计算精度远高于入射波。此外,计算实例分析显示,初至波中富含折射波和陡倾构造的反射波在很大偏移距还能被接收,崎岖海底各种波型传播路径复杂;基于此提出加大采集排列长度和采用直达波走时可提高复杂构造成像质量等复杂海底地震数据采集与处理方面的思考与建议。 相似文献
13.
Investigation of highly efficient algorithms for solving linear equations in the discontinuous deformation analysis method 下载免费PDF全文
Large‐scale engineering computing using the discontinuous deformation analysis (DDA) method is time‐consuming, which hinders the application of the DDA method. The simulation result of a typical numerical example indicates that the linear equation solver is a key factor that affects the efficiency of the DDA method. In this paper, highly efficient algorithms for solving linear equations are investigated, and two modifications of the DDA programme are presented. The first modification is a linear equation solver with high efficiency. The block Jacobi (BJ) iterative method and the block conjugate gradient with Jacobi pre‐processing (Jacobi‐PCG) iterative method are introduced, and the key operations are detailed, including the matrix‐vector product and the diagonal matrix inversion. Another modification consists of a parallel linear equation solver, which is separately constructed based on the multi‐thread and CPU‐GPU heterogeneous platforms with OpenMP and CUDA, respectively. The simulation results from several numerical examples using the modified DDA programme demonstrate that the Jacobi‐PCG is a better iterative method for large‐scale engineering computing and that adoptive parallel strategies can greatly enhance computational efficiency. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
14.
15.
16.
Astronomy Reports - An algorithm for the numerical solution of Kepler’s equation with machine precision is presented. The convergence of the iterative sequence of Newton’s method is... 相似文献
17.
针对岩土工程材料应变软化问题及有限元对其数值计算时切线刚度矩阵负定造成求解困难的问题进行研究。建立了基于Drucker-Prager(D-P)强度准则的岩石弹塑性应变软化本构模型,本构积分算法采用一种完全隐式返回映射算法,它具有无条件稳定和精确的特点,详细论述了如何进行本构模型的程序化求解;考虑弧长法在判断切线刚度矩阵正定性导致效率低的缺点,在弹塑性增量有限元方程的迭代计算中尝试采用Newton-Raphson法和arc-length法(NR-AL法)联合迭代求解的思路,即在结构未达到极限荷载前采用NR迭代法,而当结构接近极限荷载时转换为AL法控制迭代,从而使结构越过峰值点进入软化区直至破坏,NR-AL法汲取了2者迭代求解中具有的优势;利用C++语言对所建应变软化模型的本构求解和弹塑性增量有限元方程迭代求解过程给予程序实现,应用所编程序进行数值计算,分析了D-P理想弹塑性模型、应变软化模型、应变硬化模型计算的应力-应变曲线的区别,同时将应变软化模型计算结果与试验数据进行了对比。研究结果表明:所建应变软化本构模型可以较好地模拟岩石材料的峰后软化特性,能够揭示峰后应变软化特性和破坏机制,同时NR-AL法能够求解由于应变软化造成的负刚度问题,也克服了单独使用弧长法时判断切线刚度矩阵正定性效率低的缺点。 相似文献
18.
A parallelizable, semi‐implicit numerical method is proposed for the study of naturally‐fractured reservoir systems. It has proved to be computationally efficient in producing accurate numerical solutions for the dual‐porosity model for immiscible, two‐phase flow in such reservoirs. The method combines hybridized mixed finite elements, a new version of the modified method of characteristics, a sophisticated operator‐splitting procedure for separating the pressure calculation in the fractures from that of the saturation, another operator splitting to handle the interaction of the matrix blocks and the fractures, and domain decomposition iterative procedures for both the pressure and the saturation. It permits moderately long time steps for the pressure and the saturation in the fractures and matrix blocks by using short, inexpensive microsteps to treat the transport portion of the saturation equation in the fractures. This paper is devoted to the formulation of the method and a discussion of numerical results for five‐spot and vertical cross‐section examples. 相似文献
19.
《Computers and Geotechnics》1988,5(2):151-170
Mathematical and computational aspects of a three-dimensional boundary element formulation designed specifically for application to mining rock mechanics problems are described. Stress analysis of mining geometries differs in many respects from equivalent analyses of civil or mechanical engineering. These differences are highlighted and the manner in which these differences are accomodated by the formulation is described. The formulation has been implemented into a program BEAP (Boundary Element Analysis Package). A test problem is presented and compared with other formulations.Design goals for the program include the following: ease of use, ability to run on personal computers, use of an iterative equation solver to reduce total CPU for large problems, use of single precision to store matrix coefficients thereby reducing disk storage requirements, limited application to multi-subregion problems and easy inter-facing to third party software for pre- and post-processing.Finally, a practical mining example is presented. It is concluded that the formulation is very well suited to analysis of mining rock mechanics problems, providing a sound basis for future developments in this field. 相似文献