首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
常规滤波方法对浅海多次波压制不能取得理想效果.通过改变模型表面边界条件,求解波动方程实现多情景地震数值模拟.利用自由边界条件模拟时获得含海平面强反射多次波地震记录,而利用吸收边界条件时得到理想条件下无多次波地震记录.对比分析多次波对有效地震信号的影响得到海平面强反射多次波改变了有效地震信号的振幅、频率和相位.利用归一化方法计算得到压制多次波滤波器算子,并作用于含海平面强反射多次波地震记录.叠加地震剖面对比分析表明,归一化滤波器算子较好的压制了地震记录中的多次波,提高地震资料分辨率.最后,根据多次波产生条件,给出了利用实际区域简化地质模型的地震数值模拟结果计算归一化滤波器算子的方法,从而利用该方法压制实际地震数据中多次波.  相似文献   

2.
在线数值模拟方法是一种基于模型更新的混合试验方法,其是不完整边界条件混合试验的通用解决方案。在该方法中:结构反应由结构整体数值模型计算得到,边界条件自然得以满足;物理子结构的恢复力不再参与运动方程的求解,而是用于估计数值模型的本构参数,提高了数值模型的可靠性。论文提出多本构模型参数同步更新的在线数值模拟方法,以提高多种本构模型参数不确定性同时存在时在线数值模拟方法的精度,并以Levenberg-Marquardt算法为参数估计方法,实现了采用恢复力对钢筋和混凝土材料本构模型参数的同步估计。研究表明:以恢复力来估计材料本构模型参数是可行的,且Levenberg-Marquardt算法具有较高的收敛速度和计算精度;在线数值模拟方法有效地提高了混合试验的模拟精度。  相似文献   

3.
各向异性介质地震波场的优化褶积微分算子法数值模拟   总被引:1,自引:1,他引:0  
在前人工作基础上,通过对窗函数参数进行优化实现了对基于Shannon奇异核理论的交错网格褶积微分算子的优化过程.应用这种优化褶积微分算子方法对各向异性介质进行了数值模拟,讨论了优化褶积微分算子法模拟的PML吸收边界条件以及稳定性条件,分析了弹性波在此类介质中的传播特征,并与高阶交错网格有限差分方法进行了对比.数值实验结果表明,该方法适用于各向异性介质中弹性波场模拟,精度高,稳定性好,是一种研究复杂介质中地震波传播的有效数值方法.  相似文献   

4.
通过对数值实验的比较和分析, 提出了一种旨在改善区域耦合模式中ENSO模拟的回归修正方法. 该方法主要用于修正耦合模式中海气间交换的通量. 具体步骤如下: 首先, 利用多年的观测资料计算得到驱动海洋模式所需的动量及热量通量, 驱动海洋模式进行长期积分; 其次, 用海洋模式模拟的SST作为大气模式的边界条件, 相应积分大气模式; 再利用大气模式模拟变量和相应观测资料建立线性关系, 通过线性拟合得到修正系数; 最后, 利用随时间和空间变化的回归修正系数修正计算动量及热量通量的变量, 并用修正后的变量计算海气交换通量, 进行耦合模式积分. 同时利用一个热带太平洋-全球大气耦合模式对该方案及常用的“距平耦合”方案进行了检验. 结果表明, 该方案优于“距平耦合”方案, 不仅可以更好的控制气候“漂移”现象, 而且, 能够改善区域耦合模式在热带太平洋区域的ENSO模拟.  相似文献   

5.
将基于计算数学中Forsyte 广义正交多项式的迭积微分算子引入到地震波动方程的一阶速度--应力方程的空间微分运算中去,并采用时间错格有限差分算子替代传统的差分算子以匹配高精度的空间迭积微分算子,从而发展一种全新的地震波场正演模拟方法,来解决复杂非均匀介质模型中的波场传播问题.为了大幅衰减人工边界引起的反射,本文将完全匹配层(Perfectly Matched Layer,PML)吸收边界条件引入到所构建的方法中,以解决迭积微分算子法的边界问题.以二维波动方程为例,用迭积微分算子法实现了双相介质的地震波场正演模拟,模拟结果表明,双相介质模型较好地解释了含流体孔隙特性.同时也表明迭积微分算子法是一种非常实用、有效的数值模拟方法.  相似文献   

6.
数值频散和边界反射是频率域模拟时需要解决的两个重要问题.然而,受计算效率和分解阻抗矩阵时的内存占用量的制约,提高有限差分算子长度或增加有限差分网格数目均不是提高频率域模拟精度的最优解决方案.本文首先分析了数值频散产生的理论机制,在此基础上,推导了一种“波数补偿”的声波方程表达式来压制数值频散,并给出其物理意义,有效地改善了数值频散问题,提高了模拟精度;在边界问题上,本文采用多轴卷积完全匹配层(MCPML)边界条件代替传统的完全匹配层(PML)边界条件,快速吸收边界内的残余能量,压制边界反射.结合改进声波方程和MCPML边界条件,给出了一种高精度的频率域声波方程有限差分格式.数值模拟结果表明,在不增加计算量和内存占用量的前提下,本文研究的方法、正演精度高、波场模拟清晰、无干扰反射,是一种可靠高效的频率域模拟方法.  相似文献   

7.
利用偏移进行视反射率估计的初步研究   总被引:1,自引:0,他引:1       下载免费PDF全文
视反射率估计是地震数据处理解释中的一项重要内容,通常采用反演的方法得到.本文以地震偏移和地震线性反演理论相结合为基础,并利用保幅单程波传播算子和保幅波动方程叠前偏移算法以及成像空间中的角度域波动方程偏移成像和照明补偿等方法技术,提出了一种利用单程波波动方程偏移进行地下反射面视反射率估计方法,并进行了理论模型的数值试验.这种估计方法得到的视反射率估计是一种近法向入射的小角度反射率.  相似文献   

8.
VTI介质qP波方程高精度有限差分算子   总被引:6,自引:4,他引:6       下载免费PDF全文
波动方程有限差分法是一种使用广泛的地震波数值模拟方法.但是有限差分法本身固有存在着数值频散问题,会降低地震波场模拟的精度与分辨率.为了克服常规有限差分算子的数值频散,本文针对VTI介质地震波数值模拟问题,构造了频率-空间域qP波波动方程高精度有限差分优化算子,根据最优化理论中高斯-牛顿法确定了高精度有限差分算子的优化系数.利用常规差分算子和高精度优化差分算子对归一化相速度的频散关系精度进行了对比分析,并对均匀各向同性介质和均匀VTI介质中的qP波地震波场进行了有限差分数值模拟,通过频散关系精度分析和波场数值模拟结果表明:有限差分优化算子具有较高的波场数值模拟精度,有效压制了传统有限差分算子数值模拟中的数值频散现象,提高了有限差分算子精度,为VTI介质频率-空间域qP波正演模拟奠定了基础.  相似文献   

9.
反褶积是提高地震资料分辨率的重要方法,子波估计与反褶积算子设计是反褶积方法的两个重要方面.本文在分析谱模拟法与自相关法各自特点的基础上,将两种方法的优势相结合,提出了自适应谱模拟方法,提高了谱模拟方法的适应性和准确性.在反褶积算子计算过程中对误差能量进行加权,减少了反褶积过程中处理噪音的产生,同时引入信噪比谱约束,以达到保证信噪比的前提下合理提升地震资料分辨率的目的.  相似文献   

10.
在地震波传播数值模拟的过程中,需要使用吸收边界条件达到衰减人为边界反射的目的,其中完全匹配层(PML)吸收边界条件被认为是目前最理想的吸收边界条件.但在三维地震波传播数值模拟中直接应用常规的PML吸收边界条件存在编程实现复杂、计算资源需求高等不足之处.本文基于传统的NPML吸收边界条件,提出了一种改进的方法.该方法通过引入伸展函数和辅助变量,把PML介质中的三维弹性波动方程分解成正常计算项和衰减项,其中衰减项的计算只在PML区域内部进行,其值只与坐标轴的方向相关,可通过离散数值迭代求解.分析和试验结果表明,该方法避免了卷积运算,减少了辅助变量的个数,具有理论方程简洁、编程实现简单、计算资源需求少等优点,能更有效地解决三维地震波传播数值模拟中的吸收边界问题.  相似文献   

11.
波场模拟中的数值频散分析与校正策略   总被引:22,自引:5,他引:17       下载免费PDF全文
波动方程有限差分法正演模拟,对认识地震波传播规律、进行地震属性研究、地震资料地质解释、储层评价等,均具有重要的理论和实际意义.但有限差分法本身固有存在着数值频散问题,数值频散在正演模拟中是一种严重的干扰,会降低波场模拟的精度与分辨率.针对TI介质波场模拟的交错网格有限差分方法,本文从空间网格离散、时间网格离散和算子近似等三个方面对其产生的数值频散进行了分析,并结合其他学者的研究成果给出了TI介质波场模拟中压制数值频散的方法与策略:在已知介质频散关系时,对差分算子可实施算子校正;通过提高差分方程的阶数来提高波场模拟精度;采用流体力学中守恒式方程的通量校正传输方法来压制波场模拟中的数值频散;在实际正演模拟时,采用交错网格高阶有限差分方程,不仅在空间上采用高阶差分,而且在时间上也要采用高阶差分,否则只在单一方向上(空间或时间)提高方程的阶数对压制数值频散也不会取得理想的效果.  相似文献   

12.
用伴随方法求解多个工业污染源优化布局问题   总被引:5,自引:1,他引:5  
刘峰  胡非  朱江 《中国科学D辑》2005,35(1):64-71
发展了一种伴随方法, 引入罚函数来处理环境约束不等式, 构造优化问题的Lagrange函数, 从理论上推导了伴随方程和目标函数梯度公式. 通过求解伴随方程得到目标函数的梯度, 利用梯度信息使得目标函数下降, 用迭代的方法逼近最优解. 用二维的简化模式进行了数值试验, 结果验证了理论的正确性. 该方法收敛速度快, 效率高, 与空气质量数值模式和计算机技术的飞速发展相适应, 为工业规划和污染防治提供了有力的研究工具.  相似文献   

13.
复杂介质地震波传播的褶积微分算子数值模拟h   总被引:1,自引:1,他引:0       下载免费PDF全文
李信富  李小凡 《地震学报》2008,30(4):377-382
将佛尔塞(Forsyte)广义正交多项式微分算子地震波场正演模拟算法运用于复杂非均匀介质模型的波场数值模拟中,并比较了该方法在计算效率和计算精度方面与有限差分方法和伪谱法的差异. 数值结果表明,这种广义正交多项式微分算子法计算速度快、精度高,对计算资源需求低,是一种颇具潜力的数值模拟方法.   相似文献   

14.
本文导出了一种由单程波方程利用反射/透射算子的可分表示方法模拟复杂介质中一次反射地震波的数值算法. 文中利用算子可分表示理论将反射/透射算子分解成适合于双域(空间域和波数域)运算的表达形式,使得本文得到的地震波数值模拟算法可适应于一定程度横向非均匀介质和界面起伏情况,在入射角小于45°时能够准确模拟振幅随入射角(AVA)的关系. 就模拟一次反射地震波而言,与前人研究的双程波动方程伪谱法地震波模拟相比,本文算法具有足够高的模拟精度,且计算效率成倍地提高.  相似文献   

15.
本文以多孔介质中大尺度传热问题为基础,结合热平衡理论分析与数值计算,探讨了上通流对大陆岩石圈地幔-地壳热结构模式的潜在影响.根据大陆岩石圈中孔隙波传热概念模型的初步理论分析结果,指出了采用理论分析和数值模拟相结合的方法在研究大陆岩石圈地幔-地壳热结构模式时的重要性.理论分析方法可用来确定岩石圈尺度范围内大陆岩石圈的厚度和大陆地壳相关的边界条件,从而为地壳范围内数值模型的建立提供一些重要信息.数值模拟方法可以用来模拟地壳尺度范围内地壳的详细结构和复杂几何形状.如果地壳内的热分布是所考虑的主要因素,采用具有地壳尺度的合理数值模型可以有效减少计算机工作量.利用理论分析方法求出的岩石圈尺度范围内大陆岩石圈厚度与地幔传导热流之间关系的理论解,不仅可以用来验证模拟大陆岩石圈内传热问题所采用的数值方法, 而且可以用来初步研究大陆岩石圈内热分布的基本规律,为研究岩石圈地幔热事件中大陆岩石圈热减薄过程提供相应的边界条件.本文从理论分析的观点初步探讨了中国大陆不同构造背景下大陆岩石圈的热结构模式,其结果与从地球物理和地质资料中获得的大陆岩石圈热结构模式十分吻合.研究结果表明由大陆岩石圈中孔隙波传播所导致的上通流是影响大陆岩石圈地幔-地壳热结构模式及大陆岩石圈地幔与地壳之间物质和能量交换的可能机制之一.  相似文献   

16.
高精度频率域弹性波方程有限差分方法及波场模拟   总被引:14,自引:4,他引:14       下载免费PDF全文
有限差分方法是波场数值模拟的一个重要方法,但常规的有限差分法本身存在着数值频散问题,会降低波场模拟的精度与分辨率,为了克服常规差分算子的数值频散,本文采用25点优化差分算子,再根据最优化理论求取的优化系数,建立了频率空间域中弹性波波动方程的差分格式;为了消除边界反射,引入最佳匹配层,构造了各向同性介质中弹性波方程在不同边界和角点处的边界条件. 最后由弹性波波动方程和边界条件,通过频率域有限差分法,分别利用不同震源对弹性波在均匀各向同性介质、层状介质及凹陷模型中的传播过程进行了数值正演模拟,得到了单频波波场、时间切片和共炮点道集,为下一步的研究工作(如成像、反演)提供了研究基础.  相似文献   

17.
在浅水情况下,由于观测数据中缺少近偏移距信息,水层多次波的压制面临挑战.利用多道预测算子压制水层多次波是浅水环境下压制多次波的重要方法之一,这种方法先从输入数据中估计出多道预测算子,再将预测算子和输入数据做褶积预测出水层相关多次波.然而,估计的多道预测算子很容易受噪声污染,从而影响多次波模型的精度.所以,我们提出了改进的多道预测算子压制浅水多次波方法.该方法先从数据中估计出多道预测算子,并利用估计的算子构建出精确的水层模型;然后,通过计算算子的走时信息、估计振幅信息、合成新算子三个步骤来修正原始的多道预测算子.修正的算子不仅不受噪声影响,还含有精确的走时信息、可靠的振幅信息;最后,该方法用修正的算子来预测多次波,并结合自适应相减,将预测的多次波从输入数据中去除.通过合成数据和实际资料的验证表明,相比于原始的多道预测算子压制多次波方法,改进的方法能够取得更好的压制效果.  相似文献   

18.
地震波正演数值模拟中所产生的数值频散严重地影响波场模拟的精度和分辨率.基于声波方程,以8阶精度的近似解析离散算子来校正差分算子对空间偏导数进行离散,以三阶Runge-Kutta方法对时间导数进行离散,并结合通量校正传输技术以抑制数值频散,从而建立基于NAD算子和TCF技术改进的Runge-Kutta方法.波场模拟结果表明,该方法在抑制数值频散上明显优越于传统的高阶有限差分方法.同时,该方法对模拟地震波在复杂介质中的传播有着很强的适应能力.  相似文献   

19.
江文滨  林缅 《地球物理学报》2011,54(6):1679-1689
本文采用POM(Princeton Ocean Model)模式,结合多波束测量得到的高分辨率地形数据,通过网格嵌套分级计算,得到了南海北部琼西南海域近40km2研究区内的水平网格分辨率达200 m的海底流场.通过数值试验,确定了内嵌模型的最佳侧面开边界条件.计算结果表明,本文所提出超高水平分辨率嵌套模式能够模拟复杂海...  相似文献   

20.
应用波动方程有限差分方法人工合成地震记录进行城市隐伏断层的地震勘探.采用适当的差分算法、震源子波以及有效地边界条件和频散消除方法来提高正演模拟的精度和分辨率.结果表明地震数值模拟可以用来分析地层厚度、断层倾角、断层深度以及断裂带宽度等参数对地震记录的影响;通过实际试验数据和正演模拟的合成单炮地震记录的对比,可以快速判断断层的大致位置和断裂带的大致范围以及目标地层的大致深度.该方法能有效辅助实际野外工作中的数据采集参数的估计,提高勘探效率和精度.  相似文献   

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

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