共查询到19条相似文献，搜索用时 109 毫秒
1.
从三维非均匀介质中的波动方程出发，利用拟微分算子理论，Ｐａｄｅ逼近方法及因式分解技巧，获得了非均匀介质的三维高阶深度偏移方程，相应地提出了逐次低阶方程方法、低阶方程组方法及分裂方法等３种求解方法．与二维情形不同，以上每一种方法在数值求解时均存在由测线坐标ｙ的出现而带来的困难．为了克服这一困难，我们提出了差分算子分解方法，避免了近年来人们竞相研究的ｘ－ｙ方向微分算子分裂带来的分裂误差，保持了应有的相容性，解决了这一令人烦恼的问题． 相似文献
2.
利用Claerbout方程进行地震资料偏移,只适用于小倾角的情况。为了克服这一限制,R.Stolt和A.Berkhout等人导出了高阶近似的单程波方程,它们是比较复杂的高阶偏微分方程,在数值求解上存在一定的困难。本文讨论了低阶方程组形式的高阶近似,对它们构造了一些合适的差分格式。提出了求解这些差分方程的具体算法,并与15°差分偏移算法相比较,分析了此算法的计算工作量。本文提出的大倾角差分偏移方法十分有效且容易实现。 相似文献
3.
利用传统有限差分方法对基于Biot理论的双相介质波动方程进行数值求解时,由于慢纵波的存在,数值频散效应较为明显,影响模拟精度.相对于声学近似方程及普通弹性波方程,Biot双相介质波动方程在同等数值求解算法和精度要求条件下,其地震波场正演模拟需要更多的计算时间.本文针对Biot一阶速度-应力方程组发展了一种变阶数优化有限差分数值模拟方法,旨在同时提高其正演模拟的精度和效率.首先结合交错网格差分格式推导Biot方程的数值频散关系式.然后基于Remez迭代算法求取一阶空间偏导数的优化差分系数,并用于Biot方程的交错网格有限差分数值模拟.在此基础上把三类波的平均频散误差参数限制在给定的频散误差阈值和频率范围内,此时优化有限差分算子的长度就能自适应非均匀双相介质模型中的不同速度区间.数值频散曲线分析表明:基于Remez迭代算法的优化有限差分方法相较传统泰勒级数展开方法在大波数范围对频散误差的压制效果更明显;可变阶数的优化有限差分方法能取得与固定阶数优化有限差分方法相近的模拟精度.在均匀介质和河道模型的数值模拟实验中将本文变阶数优化有限差分算法与传统泰勒展开算法、最小二乘优化算法进行比较,进一步证明其在复杂地下介质中的有效性和适用性. 相似文献
4.
5.
从麦克斯韦旋度方程出发可以直接导出瞬变电磁场扩散方程,然而扩散方程不含电场对时间的一阶导数,不能构成显式的时域有限差分方程,借鉴du Fort-Frankel有限差分离散方法引入虚拟位移电流项构建显式时域有限差分方程.对Wang和Hohmann的经典时域算法进行了两点改进:第一,通过将矩形回线源电流密度加入麦克斯韦方程组的安培环路定理方程,实现回线源瞬变电磁激发源加入;第二,在计算中考虑关断时间.第一点改进使时域有限差分方程考虑了一次场的计算,并且源的计算不再依赖均匀半空间模型响应作为初始条件,使算法能够适应表层电阻率不均匀时的三维复杂模型.由于实际观测中不可能出现阶跃电流的关断形式,第二点改进可以方便设置发射电流下降沿.采用改进的三维时域有限差分正演算法对均匀半空间模型、四类三层模型、均匀半空间中含有低阻块体模型进行了计算并分别与解析解、线性数字滤波解、积分方程解和Wang的三维时域有限差分解进行了对比验证.以H模型为例,采用建立的三维时域有限差分正演算法计算了不同关断时间的斜阶跃脉冲回线源瞬变电磁中心点感应电动势衰减曲线.以实际地质资料为基础,构建包含两层采空区的三维复杂模型,以1 μs的极短关断时间进行了复杂模型定回线源瞬变电磁响应计算,并计算了该复杂模型的视电阻率曲线. 相似文献
6.
波动方程有限差分法是一种使用广泛的地震波数值模拟方法.但是有限差分法本身固有存在着数值频散问题，会降低地震波场模拟的精度与分辨率.为了克服常规有限差分算子的数值频散，本文针对VTI介质地震波数值模拟问题，构造了频率-空间域qP波波动方程高精度有限差分优化算子，根据最优化理论中高斯-牛顿法确定了高精度有限差分算子的优化系数.利用常规差分算子和高精度优化差分算子对归一化相速度的频散关系精度进行了对比分析，并对均匀各向同性介质和均匀VTI介质中的qP波地震波场进行了有限差分数值模拟，通过频散关系精度分析和波场数值模拟结果表明：有限差分优化算子具有较高的波场数值模拟精度，有效压制了传统有限差分算子数值模拟中的数值频散现象，提高了有限差分算子精度，为VTI介质频率-空间域qP波正演模拟奠定了基础. 相似文献
7.
有限差分方法是波场数值模拟的一个重要方法,交错网格差分格式比规则网格差分格式稳定性更好,但方法本身都存在因网格化而形成的数值频散效应,这会降低波场模拟的精度与分辨率.为了缓解有限差分算子的数值频散效应,精确求解空间偏导数,本文把求解波动方程的线性化方法推广到用于求解弹性波方程交错网格有限差分系数;同时应用最大最小准则作为模拟退火(SA)优化算法求解差分系数的数值频散误差判定标准来求解有限差分系数.通过上述两种方法,分别利用均匀各向同性介质和复杂构造模型进行了数值正演模拟和数值频散分析,并与传统泰勒展开算法、最小二乘算法进行比较,验证了线性化方法和模拟退火方法都能有效压制数值频散,并比较了各个算法的特点. 相似文献
8.
地震波场数值模拟是理解地震波在地下介质中的传播特点，帮助解释观测数据的有效手段，而提高计算精度和运算效率是所有波场数值模拟方法研究所追求的目标.有限差分技术是求解波动方程计算效率最高、应用最为广泛的方法之一.但传统的有限差分技术计算过程中的数值频散问题影响了该技术的计算精度与计算效率.本文通过交错网格高阶有限差分技术与通量校正传输方法（Flux|corrected transport method，FCT）相结合, 对横向各向同性介质（Transverse isotropic medium，TI）一阶速度|应力弹性波动方程组进行了数值求解研究.波场快照数值模拟结果表明，本文研究的数值模拟方法与波动方程二阶有限差分方法、交错网格四阶有限差分方法相比，在压制网格数值频散方面有明显的优势，计算精度提高，而且可以利用较大的空间步长，提高计算效率. 相似文献
9.
10.
11.
12.
波动方程有限差分法正演模拟,对认识地震波传播规律、进行地震属性研究、地震资料地质解释、储层评价等,均具有重要的理论和实际意义.但有限差分法本身固有存在着数值频散问题,数值频散在正演模拟中是一种严重的干扰,会降低波场模拟的精度与分辨率.针对TI介质波场模拟的交错网格有限差分方法,本文从空间网格离散、时间网格离散和算子近似等三个方面对其产生的数值频散进行了分析,并结合其他学者的研究成果给出了TI介质波场模拟中压制数值频散的方法与策略：在已知介质频散关系时,对差分算子可实施算子校正;通过提高差分方程的阶数来提高波场模拟精度;采用流体力学中守恒式方程的通量校正传输方法来压制波场模拟中的数值频散;在实际正演模拟时,采用交错网格高阶有限差分方程,不仅在空间上采用高阶差分,而且在时间上也要采用高阶差分,否则只在单一方向上（空间或时间）提高方程的阶数对压制数值频散也不会取得理想的效果. 相似文献
13.
横向各向同性介质紧致交错网格有限差分波场模拟(英文) 总被引：4，自引：2，他引：2
针对有限差分数值模拟的频散问题,本文将交错网格技术和紧致差分格式相结合,推导了横向各向同性介质一阶速度一应力波动方程的紧致交错网格差分格式;对比分析了紧致交错网格差分格式、交错网格差分格式以及紧致差分格式的截断误差主项,并利用Fourier误差分析方法分析了上述三种差分格式的近似精度;在此基础上,分别采用上述三种差分格式进行了波场数值模拟。结果表明,当差分方程阶数相同时,紧致交错网格差分格式截断误差最小,数值频散最弱,差分精度最高,证实了该方法的有效性。 相似文献
14.
The accuracy of finite difference computations for the dynamic motion of cylindrical shells, including transverse shear and rotatory inertia, has been assessed by comparison with Fourier series solutions. The finite difference models were based upon either the differential equations or upon control volume concepts, with the results of the latter being less sensitive to mesh size and in better agreement with the Fourier series results. A frequency analysis of all of the finite difference algorithms and the Fourier series demonstrated that for shells with clamped ends the finite difference spectra had a very considerable gap which closed slowly as the mesh size decreased, with the spectra being most complete for the best algorithm. A new algorithm was based upon these spectra and is shown to yield good results while permitting increases in the computational speed by the factor of mesh size/shell thickness. 相似文献
15.
A new upstream weighting finite element technique is developed for improved solution of the two-phase immiscible flow equations. Unlike the upstream weighting technique used by previous investigators, the new technique does not employ finite difference concepts to achieve the required upstream weighting of relative permeabilities or mobilities. Instead, upstream weighting is achieved by (1) representing the relative permeabilities or mobilities as continuous functions expressed in terms of the shape functions and nodal values (2) using asymmetric weighting functions to weight the spatial terms in the flow equations. These weighting functions are constructed such that they are dependent on the flow direction along each side of an element.In conjunction with the proposed technique, two solution schemes for treating the resulting set of non-linear algebraic equations are presented. These are the fully-implicit chord slope incremental solution scheme and the Newton-Raphson solution scheme. Both schemes allow the use of large time steps without being unstable.The proposed numerical technique is applied to two problems (1) the one-dimensional Buckley-Leverett problem (2) the two-dimensional five-spot well flow problem. Results indicate that this technique is superior to not only earlier finite element schemes but also five-point upstream finite difference formulae. 相似文献
16.
有限差分方法是波场数值模拟的一个重要方法,但常规的有限差分法本身存在着数值频散问题，会降低波场模拟的精度与分辨率，为了克服常规差分算子的数值频散，本文采用25点优化差分算子，再根据最优化理论求取的优化系数，建立了频率空间域中弹性波波动方程的差分格式；为了消除边界反射，引入最佳匹配层，构造了各向同性介质中弹性波方程在不同边界和角点处的边界条件. 最后由弹性波波动方程和边界条件，通过频率域有限差分法，分别利用不同震源对弹性波在均匀各向同性介质、层状介质及凹陷模型中的传播过程进行了数值正演模拟，得到了单频波波场、时间切片和共炮点道集，为下一步的研究工作（如成像、反演）提供了研究基础. 相似文献
17.
质点的轨迹计算是半拉格朗日模式的重要基础,传统的数值计算方法由于采用时间差分代替微分,只能得到质点运动轨迹终点的速度,因此质点的移动轨迹（位移）只能靠风速外推的方法计算,导致了模式计算不稳定等问题.借鉴精细积分法中使用半解析解的思路,利用正压原始方程研究了用运动方程的半解析解构建数值模式的可能性.求解了运动方程的一阶和二阶微分方程组的半解析解,通过时间积分半解析解计算质点运动轨迹.数值试验表明,一阶微分方程组的半解析解比差分解略有优势.二阶微分方程组的半解析解在时间步长增大时优势非常明显,而且在保证计算精度的前提下,节省计算时间,这对提高模式性能有重要作用. 相似文献
18.
本文以基于改进BISQ模型的二维双相各向同性介质一阶速度-应力方程为基础,推导出了曲线坐标系下对应的方程,然后采用低频散、低耗散的同位网格MacCormack有限差分法来离散方程,并采用紧致的单边MacCormack差分格式结合牵引力镜像法来施加自由地表边界条件,实现了地震波场数值模拟.曲线网格有限差分法采用贴体网格来描述自由表面,地表的网格线紧贴地形,避免了台阶近似造成的数值散射.数值模拟结果表明,在双相介质起伏自由地表和分界面处,各类波型复杂的反射透射规律可以清晰展现,曲线网格有限差分法可以精确地解决地震波在含起伏地表的双相各向同性介质中的传播问题. 相似文献
19.
The stability problem is a very important aspect in seismic wave numerical modeling. Based on the theory of seismic waves and constitutive equations of viscoelastic models, the stability problems of finite difference scheme for KelvinVoigt and Maxwell models with rectangular grids are analyzed. Expressions of stability conditions with arbitrary spatial accuracies for two viscoelastic models are derived. With approximation of quality factor Q≥?5, simplified expressions are developed and some numerical models are given to verify the validity of the corresponding theoretical results. Then this paper summarizes the influences of seismic wave velocity, frequency, size of grid and difference coefficients, as well as quality factor on stability condition. Finally the prerequisite conditions of the simplified stability equations are given with error analysis. 相似文献