共查询到19条相似文献,搜索用时 187 毫秒
1.
地震波传播的褶积微分算子法数值模拟 总被引:2,自引:0,他引:2
为了解决地震波场数值模拟中的速度与精度的匹配及局部信息与全局信息的耦合等问题, 该文借助新推出的基于Forsyte广义正交多项式的褶积微分算子, 将计算数学中的Forsyte多项式褶积微分算子应用到地震波传播的数值模拟中.复杂非均匀介质模型数值模拟结果说明了该方法的可行性和优越性.该方法同时具有广义正交多项式褶积微分算子的高精度和有限差分短算子算法的高速度.通过对算子长度的调节及算子系数的优化, 可同时兼顾波场解的全局信息与局部信息. 相似文献
2.
褶积微分算子法是一种全新的数值模拟方法,已被广泛应用于复杂介质的地震波场数值模拟,但是其边界反射问题一直没有解决。这里将最佳匹配层(PML)吸收边界条件引入到褶积微分算子法中,此方法是在研究区域的边界上加入吸收层,使边界上传入吸收层的波,随传播距离按指数规律衰减,不产生任何反射,以达到消除边界反射的目的。构造不同的模型,通过对比分析证明,PML吸收边界条件能比较好地解决褶积微分算子法的边界问题,从而验证了完全匹配层吸收效果的优越性。 相似文献
3.
地震波场数值模拟方法多种多样,各种方法都有各自的特点.这里推出一种全新地震波场模拟方法--基于Forsyte广义正交多项式的迭积微分算子法,该方法将计算数学中的Forsyte多项式,应用到地震波传播的数值模拟中,它同时具有广义正交多项式迭积微分算子的高精度和有限差分短算子算法的高速度.通过对算子长度的调节及算子系数的优化,可同时兼顾波场解的全局信息与局部信息.复杂非均匀介质模型数值模拟结果说明了该方法的可行性.将该方法的计算结果与傅氏变换伪谱法、错格高阶有限差分法相比较,结果说明,该方法在波场模拟方面具有较好的发展潜力,并具有自身独到的优越性. 相似文献
4.
杨辉 《物探化探计算技术》1998,20(4):347-350
根据台劳级数展开式,将窗口内网格结点位场数据在计算点展开,且对应点两两相加,由最小二乘法求得窗口内网格结点的圆滑滤波算子和二次导数算子,将此算子与观测的位场数据褶积,就可得到位场圆滑后的数据及水平、垂直二次导数。 相似文献
5.
杨辉 《物探化探计算技术》1997,19(3):246-251
根据台劳级数展开式,将窗口内网格结点位场数据在计算点展开,且对应点两两相减,由最小二乘法求得窗口内网格结点的水平导数算子,将此算子与观测的位场数据褶积,就可求得位场的水平导数。 相似文献
6.
7.
我们建立了一个频率域的位场转换解释系统,由于可以利用褶积定理(空间域的褶积对应于频率域的乘积)的特点,在导出一系列位场转换的频率域滤波算子的基础上,不但可以实现各种单一的位场转换,而且可以很方便地进行综合的位场转换——将若干个频率域的滤波算子连乘(滤波器的串联),从而使这个解释系统具有多功能位场转换的特点,以适应解释的要求。此外,借助于快速富里叶变换的计算机算法,使这个系统在计算机上能省时而准确地实现计算。鉴于上述两方面的优越性,近几年来这个系统在重磁资料解释实践中,特别是在解释宽缓异常时,逐步地得到了推广使用。 相似文献
8.
9.
10.
谱分析方法正演位场的若干问题 总被引:1,自引:0,他引:1
谱分析方法正演重磁场时,谱的离散化和截断产生异常的褶迭和皱波效应给计算结果带来了误差.本文提出消除办法,即先对谱进行褶迭和褶积,再作富里叶变换,得到的空间域异常就不再产生褶迭和皱波效应了.本文还对谱的解析式中某些奇异点的处理提出解决措施. 相似文献
11.
地震波混合阶褶积算法模拟 总被引:6,自引:0,他引:6
基于正反傅立叶变换,提出了地震波模拟的混合阶褶积算法。该方法原理简单、易于实现,结合了有限差分与伪谱法的优点,具有较高的精度和计算效率,适用于地震波场的正演计算,同时还给出了二维地震波场的理论计算实例。计算结果表明,此算法模拟结果正确、精度高、速度快、能适应较为复杂地质模型,并且易于推广到各向异性介质中去。 相似文献
12.
13.
基于散度和旋度的弹性波波场分离数值模拟方法 总被引:4,自引:1,他引:3
弹性波波场分离数值模拟方法是研究混合波场中纯纵波和纯横波波场传播规律的一种重要手段,通常采用算子分裂来实现,这不但增加了计
算量,而且还浪费了大量计算时间。因此提出了基于散度和旋度的弹性波波场分离数值模拟新方法,仅依赖于交错网格来计算出散度场(纯纵波)
和旋度场(纯横波),从而提高计算效率。在人工截断边界处采用完全匹配层吸收边界条件(PML)。数值模拟结果表明,该方法能够从混合波场中
准确地分离出纯纵横波波场,且精度高,边界吸收效果好,计算速度快,计算过程稳定,可用于指导实际多波多分量地震资料处理,同时指出本文
算法对相互耦合的波型(如面波)无法实现准确分离。 相似文献
14.
15.
The offset-domain prestack depth migration with optimal separable approximation, based on the double square root equation, is used to image complex media with large and rapid velocity variations. The method downward continues the source and the receiver wavefields simultaneously. The mixed domain algorithm with forward Fourier and inverse Fourier transform is used to construct the double square root equation wavefield extrapolation operator. This operator separates variables in the wave number domain and variables in the space domain. The phase operation is implemented in the wave number domain, whereas the time delay for lateral velocity variation is corrected in the space domain. The migration algorithm is efficient since the seismic data are not computed shot by shot. The data set test of the Marmousi model indicates that the offset-domain migration provides a satisfied seismic migration section on which complex geologic structures are imaged in media with large and rapid lateral velocity variations. 相似文献
16.
17.
基于一阶速度—应力弹性波波动方程,利用傅里叶变换将空间导数变换到频率域求解,对介质波场进行正演模拟。本文结合实际勘探中可能存在的地质问题,设计了多种地质模型,包括背斜模型、向斜模型、透镜体模型和断层--侵入体模型等,进行了伪谱法地震波场正演模拟、合成地震记录的实现,并分析了其波场动力学变化特征。结果表明,应用伪谱法对复杂模型的模拟有较高的精度和可靠性,可直接用于辅助地震野外工作设计、反演解释和资料解释等工作。 相似文献
18.
基于波动方程照明分析原理和组合震源激发理论,对二维声波波场进行规则网格有限差分正演模拟,并给出组合震源照明度的计算方法,对于多个震源组合激发,通过改变相邻两个震源之间的激发时间,激发所产生的地震波场会沿着某一确定的方向传播,可以增强地下特定的区域的地质体的照明分析。模型试验结果表明,震源数目越多,激发的地震波最大能量指向越集中,改变震源之间的延迟时间,可以调整波场的传播方向,该方法对于分析简单的二维模型以及复杂的逆冲推覆构造模型和marmousi模型的照明能量分布均适用,对于指导炮点的设计有借鉴意义。 相似文献
19.