在半球形区域上三维波动方程正演问题的研究

发布时间：2020-08-03 12:20

【摘要】：在地震勘探领域经常用到正演数值模拟技术,其在地震资料、解释以及观测系统设计等方面发挥着重要作用.波动方程正演数值模拟技术步骤是,已知数学模型、地震源和在地下的几何界面、物性参数(例如岩层密度,速度等)等条件下,模拟地震波在各种介质中的传播,形成人工地震记录.广泛研究的波动方程有两种类型:一是声波方程,二是弹性波方程.随着科技的发展和地震波在地下地质中研究的深入,地震波技术的模拟维数从二维向三维发展,地震波三维数值模拟可以更详细地研究地下介质对地震波场的影响.通常求解这类问题时最常用的求解波动方程正演问题的方法是有限差分法、有限元法和有限体积法等.本文首先介绍了波动方程正演的研究背景,研究意义及其进展情况,然后将在半圆形区域上研究二维波动方程的正演问题,首先主要针对二维波动方程在半圆形区域上的特点,并选取二维波动方程为具体的数学模型.通过使用极坐标方法将方程进行转化,得到在极坐标下的二维波动方程.然后利用有限差分方法进行离散得到离散方程,给出其边界条件,从而形成二维波动方程的正演的差分格式.为了验证算法的算法可行性,进行了MATLAB数值模拟.然后在二维波动方程正演问题的研究基础上,进行了对地震波的三维数值模拟技术的研究,利用球坐标转换方式建立了半球形区域上三维波动方程的模型。并利用有限差分方法进行离散,得到离散方程,对其边界条件进行处理,从而形成三维波动方程的正演的差分格式.通过MATLAB进行数值模拟,然后对比其位移函数的绝对误差,当极角划分越细时,位移函数的绝对误差就越小,处理的边界条件就越好。数值模拟结果充分表明了所构造的数值正演方法能比较有效地处理地震勘探正演问题,在一定程度上克服波动方程计算困难,不仅具有理论上的创新意义,而且具有较强的可行性、较高的效率.
【学位授予单位】：哈尔滨工程大学
【学位级别】：硕士
【学位授予年份】：2018
【分类号】：O241.8
【图文】：

波形,子波,函数,奇异性

况中波动方程在震源点的解具有奇异性，所以研究r 子波函数.其主要形式有以下三种[45]:( )2n 2 expftftππ ;( ))( )2022 202f tf t t eππ ;( ( ))( )222012 22 2x 01 2f t tf t t eππ . f 为主频，波形如图：

震源函数

)exp(2)sin(2)~(22~220,120,0,21,0,20,10,uRvuRvuuτ vπfnτπfnτmnmnmmnmmnm=+ + + (m = 1,2, , M 1;n=1,2, )值模拟检验构造的极坐标模型有限差分法的可行性和效率，对已构造的数值B 数值仿真实验，并对仿真效果进行分析.仿真实验步骤如下：使用多个发射器，为了得到更好的效果，考虑选雷克子波ftstfteππ2()sin(2) = 作为震源函数，其函数形状如图 3-1 所Z为震源函数的频率；取半圆形区域的半径 R =1000，极半径方向的，时间部长 τ =0.001，半圆形极角的等分个数分别为 M =10、 M =20表的每个等分点上.为了是正演模拟过程简便，假设时间步长保持不变测值均为常值 v = 2200m/s，故波长正演模拟中介质速度，二维波动值模拟如下图：

波场,节点,介质速度

生成地下各个节点处的波场的介质速度

【参考文献】