首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到18条相似文献,搜索用时 62 毫秒
1.
褶积微分算子法是一种全新的数值模拟方法,已被广泛应用于复杂介质的地震波场数值模拟,但是其边界反射问题一直没有解决。这里将最佳匹配层(PML)吸收边界条件引入到褶积微分算子法中,此方法是在研究区域的边界上加入吸收层,使边界上传入吸收层的波,随传播距离按指数规律衰减,不产生任何反射,以达到消除边界反射的目的。构造不同的模型,通过对比分析证明,PML吸收边界条件能比较好地解决褶积微分算子法的边界问题,从而验证了完全匹配层吸收效果的优越性。  相似文献   

2.
地震波场模拟中的褶积微分算子法   总被引:7,自引:1,他引:7  
褶积微分算子是将偏导数算子的带限傅氏谱进行反傅氏变换得到的,具有对称或反对称和衰减速度快的特点。根据频谱分析理论,选用序列在结束时能光滑趋于零的截断窗口将褶积算子截断,结合考虑微分算子守恒定律,将褶积算子中心值进行修正,得到适合于波场模拟的最优权系数。给出用褶积微分算子模拟声波波场的实现途径与实例。理论分析与模型试算证明了褶积微分算子法的正确性与有效性。  相似文献   

3.
地震波混合阶褶积算法模拟   总被引:6,自引:0,他引:6  
基于正反傅立叶变换,提出了地震波模拟的混合阶褶积算法。该方法原理简单、易于实现,结合了有限差分与伪谱法的优点,具有较高的精度和计算效率,适用于地震波场的正演计算,同时还给出了二维地震波场的理论计算实例。计算结果表明,此算法模拟结果正确、精度高、速度快、能适应较为复杂地质模型,并且易于推广到各向异性介质中去。  相似文献   

4.
地震波场数值模拟方法多种多样,各种方法都有各自的特点.这里推出一种全新地震波场模拟方法--基于Forsyte广义正交多项式的迭积微分算子法,该方法将计算数学中的Forsyte多项式,应用到地震波传播的数值模拟中,它同时具有广义正交多项式迭积微分算子的高精度和有限差分短算子算法的高速度.通过对算子长度的调节及算子系数的优化,可同时兼顾波场解的全局信息与局部信息.复杂非均匀介质模型数值模拟结果说明了该方法的可行性.将该方法的计算结果与傅氏变换伪谱法、错格高阶有限差分法相比较,结果说明,该方法在波场模拟方面具有较好的发展潜力,并具有自身独到的优越性.  相似文献   

5.
我们提供了一种在松散2D网格上模拟地震波传播的数值方法。该算法以弹性波方程的速度——应力公式为基础,因此利用一种交错排列网格法。不像有限单元或频谱单元法,这种方法也能处理柔性松散网格,我们利用显函数微分算子计算每一时间步的空间导数。  相似文献   

6.
笔者基于离散Shannon奇异核褶积微分算子(GDSCD)计算地震波速度应力方程的空间导数,推导了一阶GDSCD的具体形式,并提出了优化方法,即在频率域逼近平面波的真实导数,得到了不同半径和采样下限的最优窗函数系数,通过滤波响应分析算子精度,与多种数值方法对比以及模型测试表明,笔者构造的最优化GDSCD模拟地震波具有较高的计算效率和精度.  相似文献   

7.
单程波法地震波衰减特征数值模拟分析   总被引:1,自引:1,他引:1  
在非均匀粘弹性介质中,地震波衰减研究能够有效揭示地层吸收特性对地震波的影响。这里从单程波方程出发,结合粘弹性介质中的复波数理论,考虑介质的横向变化,采用分步傅立叶波场延拓技术,推导出适用于横向变速介质的带Q地震记录计算公式,最后将此方法应用于四个不同含Q地层模型。试验表明:单程波法计算速度快,波场信息简单清晰,易于深入分析和反演算法的建立。并且,Q值对地震波波形、能量、频率等都产生了较为明显的影响,对实际资料能量补偿、地震资料分辨率的提高,以及利用衰减特征进行储层识别有一定指导价值。  相似文献   

8.
谱元法是有限元法的一个重要分支,具有对模型适应能力强、计算精度较高和易于实现等优点.有别于有限元的近似计算,谱元法的质量矩阵采用数值积分得到的对角阵,既简化了计算,又可避免计算精度的下降.以用于通用计算领域的GPU为例,介绍了在CUDA编程平台下实现谱元法地震波数值模拟并行化的方法.在实例应用中为避免浮点数据的原子操作,提高计算效率,模型被采样为串行执行的四个并行单元集合体.  相似文献   

9.
针对陷落柱与周围煤体之间存在明显波阻抗差异的特点,建立巷道内地震波超前探测陷落柱数值模型。采用波场正演模拟方法,分析巷道前方存在陷落柱时的地震波场特征,总结出基于绕射偏移的陷落柱边界探测方法。在新元矿南二正巷实际超前探测中,该方法准确圈定了巷道前方陷落柱的范围。探测结果和实际揭露情况对比表明:地震波超前探测技术可以对巷道前方120m范围内的陷落柱构造进行有效预测。  相似文献   

10.
复杂地表地质条件下地震波数值模拟综述   总被引:2,自引:1,他引:2  
通过对复杂地质地表条件的分析研究,总结了近年来在这方面地震波场数值模拟的新方法。复杂地表的类型很多,对地震资料的影响程度不同,要想完全描述近地表的影响因素是不可能的。文中所介绍的模拟方法在处理这些问题上都是各有侧重、各有优势,目的是寻找合适的方法解决适合的问题。  相似文献   

11.
The perfectly matched layer(PML) was first introduced by Berenger as an absorbing boundary condition for electromagnetic wave propagation.In this article,a method is developed to ex-tend the PML to simulating seismic wave propagation in fluid-saturated porous medium.This non-physical boundary is used at the computational edge of a Forsyte polynomial convolutional differenti-ator(FPCD) algorithm as an absorbing boundary condition to truncate unbounded media.The incor-poration of PML in Biot's equations is gi...  相似文献   

12.
Doklady Earth Sciences - The problem of tsunami warning and the use of the grid-characteristic method developed for studying wave processes in heterogeneous media is addressed. A multilayer...  相似文献   

13.
This paper presents a simple and fast method to calculate flow through a dike breach. The approach was based on two-dimensional numerical simulations of idealized dike breakages at straight river-sections. As a result, computation of discharge through a breach can be achieved by use of the new developed formula (denoted as dike break formula). Furthermore, a methodology that combines one-dimensional hydrodynamic modelling, the dike break formula and a simple GIS-based method to estimate inundation areas is described. This fast and easy-to-handle tool can be utilized for near real-time forecasting or evacuation decisions. Detailed predictions were made for a number of flood and dike break scenarios at the River Rhine to prove the accuracy of the new method compared with two-dimensional numerical models.  相似文献   

14.
多波地震勘探技术能够同时获得纵波和转换波资料,与常规的纵波勘探相比,能够提供更多的地下介质信息。且转换波对于埋深较浅的小断层、小幅度构造有更高的分辨率,充分利用多分量地震资料可以有效的提高地震勘探的精度。本文将多波地震技术引入到陆域冻土区天然气水合物勘探中,对陆域冻土区水合物进行多波地震数值模拟研究。采用弹性波有限差分方法进行多分量正演模拟,利用仿射坐标转换法将原始的多分量地震数据进行纵横波叠前分离,并通过多波高斯束偏移方法分别对波场分离后的PP波和转换PS波进行叠前深度偏移成像。根据我国祁连山冻土区水合物实际地震地质资料,建立逼近真实地质情况的天然气水合物数值模型。对上述水合物模型进行多波地震数值试算表明,多波地震勘探技术是一种有效的陆域冻土区水合物探测方法,充分利用多波地震资料有利于查明陆域冻土区水合物精细的地质构造,取得更好的勘探效果。  相似文献   

15.
In the present study seismic wave propagation in heterogeneous media is numerically simulated by using the pseudospectral method with the staggered grid RFFT differentiation in order to clarify the cause for the complicated distribution characteristics of strong ground motion in regions with basin structure. The results show that the maximum amplitudes of simulated ground acceleration waveforms are closely related to the basin structure. Interference of seismic waves in the basin strongly affects the distribution of maximum seismic waveforms, which may result in peak disasters during earthquakes. Peak disasters might be away from basin boundaries or earthquake faults. Seismic energy transmitted into the basin from the bedrock can hardly penetrate the bottom of the basin and then travel back into the bedrock region. The seismic energy is absorbed by basin media, and transferred into the kinematical energy of seismic waves with great amplitude in the basin. Seismic waves between basins may result in seriou  相似文献   

16.
Homogenization in fractured rock analyses is essentially based on the calculation of equivalent elastic parameters. In this paper, a new numerical homogenization method that was programmed by means of a MATLAB code, called HLA-Dissim, is presented. The developed approach simulates a discontinuity network of real rock masses based on the International Society of Rock Mechanics (ISRM) scanline field mapping methodology. Then, it evaluates a series of classic joint parameters to characterize density (RQD, specific length of discontinuities). A pulse wave, characterized by its amplitude, central frequency, and duration, is propagated from a source point to a receiver point of the simulated jointed rock mass using a complex recursive method for evaluating the transmission and reflection coefficient for each simulated discontinuity. The seismic parameters, such as delay, velocity, and attenuation, are then calculated. Finally, the equivalent medium model parameters of the rock mass are computed numerically while taking into account the natural discontinuity distribution. This methodology was applied to 17 bench fronts from six aggregate quarries located in Tunisia, Spain, Austria, and Sweden. It allowed characterizing the rock mass discontinuity network, the resulting seismic performance, and the equivalent medium stiffness. The relationship between the equivalent Young’s modulus and rock discontinuity parameters was also analyzed. For these different bench fronts, the proposed numerical approach was also compared to several empirical formulas, based on RQD and fracture density values, published in previous research studies, showing its usefulness and efficiency in estimating rapidly the Young’s modulus of equivalent medium for wave propagation analysis.  相似文献   

17.
田小甫 《地质与勘探》2012,48(4):840-846
[摘 要] 本文采用数值模拟的方法,利用UDEC 软件,开展了结构面对岩质边坡地震动影响的数值 模拟研究,研究内容涉及结构面产状和位置变化对岩体边坡地震动放大系数空间展布的影响。根据岩 体边坡中结构面发育的一般规律,本文模型工况主要考虑与坡面相交的贯穿性结构面,模拟的主要物理 对象是层面、贯穿性长大软弱结构面等。数值模拟结果表明:结构面越靠近坡顶,边坡上部的地震动响 应越强烈;顺倾结构面会使坡肩部形成更强烈的地震动响应;边坡越陡,坡顶和坡肩部的地震动响应越 强烈。上述规律对于岩质边坡地震稳定性分析具有指导意义。  相似文献   

18.
反转构造模型的地震波模拟及成像   总被引:1,自引:0,他引:1  
为了检验由地震剖面解释得出的反转构造的可靠性,探讨反转构造地区在地震数据采集,处理中应采用的有效方法,设计了一个典型的反转构造模型-突起构造。以该模型为基础,针对复杂构造地震波模拟中常见起源频散严重、人为边界反射吸收不理想两个问题,应用高阶差分法以及一种新的边界吸收条件,对反转模型进行了地震波场数值模拟和深度偏移成像处理。采用中心放炮两边接收的观测方式,共模拟了185炮地震记录。研究表明,波动方程高阶差分法是模拟复杂构造模型中地震波传播的有效方法,除了近于直立的断层因难以接收到其反射波而较难成像外,其余的反转断层及其组合在理论上是可以通过合适方法(如Fourier有限差分深度偏移)较好地成像的。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号