期刊文献+
共找到64篇文章
< 1 2 4 >
每页显示 20 50 100
双相介质地震波场数值模拟的迭积微分算子及其PML边界条件 被引量:4
1
作者 贺同江 刘红艳 李小凡 《地球物理学进展》 CSCD 北大核心 2010年第4期1180-1188,共9页
将基于计算数学中F0rsyte广义正交多项式的迭积微分算子引入到地震波动方程的一阶速度——应力方程的空间微分运算中去,并采用时间错格有限差分算子替代传统的差分算子以匹配高精度的空间迭积微分算子,从而发展一种全新的地震波场正演... 将基于计算数学中F0rsyte广义正交多项式的迭积微分算子引入到地震波动方程的一阶速度——应力方程的空间微分运算中去,并采用时间错格有限差分算子替代传统的差分算子以匹配高精度的空间迭积微分算子,从而发展一种全新的地震波场正演模拟方法,来解决复杂非均匀介质模型中的波场传播问题.为了大幅衰减人工边界引起的反射,本文将完全匹配层(Perfectly Matched Layer,PML)吸收边界条件引入到所构建的方法中,以解决迭积微分算子法的边界问题.以二维波动方程为例,用迭积微分算子法实现了双相介质的地震波场正演模拟,模拟结果表明,双相介质模型较好地解释了合流体孔隙特性.同时也表明迭积微分算子法是一种非常实用、有效的数值模拟方法. 展开更多
关键词 地震波 数值模拟 迭积微分算子 pml边界条件 双相介质
在线阅读 下载PDF
基于PML边界的变网格高阶有限差分声波方程逆时偏移 被引量:20
2
作者 郭念民 吴国忱 《石油地球物理勘探》 EI CSCD 北大核心 2012年第2期256-265,182,共10页
叠前逆时深度偏移采用全声波方程求解,不受介质横向速度变化和高陡倾角的影响,具有成像精度高、相位准确、实现回转波成像等优点。逆时偏移利用双程波动方程构造波场延拓算子,正向延拓时间域震源点波场,逆时反向外推时间域检波点波场,... 叠前逆时深度偏移采用全声波方程求解,不受介质横向速度变化和高陡倾角的影响,具有成像精度高、相位准确、实现回转波成像等优点。逆时偏移利用双程波动方程构造波场延拓算子,正向延拓时间域震源点波场,逆时反向外推时间域检波点波场,然后利用互相关成像条件实现成像,因此正演模拟技术是其成功与否的关键。当浅层为海水或者低速层时,常规的有限差分方法必须采用小网格才能有效压制频散,得到高质量的波场记录,从而保证成像精度。但是若对整个区域都用小网格和小的时间采样间隔进行波场计算,势必造成计算量的增加。本文给出了声波方程变网格算法的差分格式,推导了基于PML边界条件的变网格高阶有限差分方程,将可变网格和可变时间步长算法应用于逆时偏移的波场外推,既保证了波场外推计算的精度和最终逆时偏移的成像效果,同时又提高了计算效率,并通过数值算例试算和逆时偏移成像的应用,说明了该方法的有效性和可行性。 展开更多
关键词 逆时偏移 pml边界条件 变网格 高阶有限差分 正演模拟
在线阅读 下载PDF
耦合PML边界条件的三维大地电磁二次场有限差分法 被引量:7
3
作者 薛帅 白登海 +1 位作者 唐静 马晓冰 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2017年第1期337-348,共12页
本文将大地电磁场分解为一次场和二次场,应用交错网格有限差分法模拟计算大地电磁二次场,并引入各向异性最佳匹配层(PML)吸收边界条件作为二次场边界条件,实现了耦合PML吸收边界条件的三维大地电磁二次场有限差分正演模拟.为了确保正演... 本文将大地电磁场分解为一次场和二次场,应用交错网格有限差分法模拟计算大地电磁二次场,并引入各向异性最佳匹配层(PML)吸收边界条件作为二次场边界条件,实现了耦合PML吸收边界条件的三维大地电磁二次场有限差分正演模拟.为了确保正演的稳定性和效率,QMR求解器和磁感应矢量散度校正技术被用于PML吸收边界条件下系数矩阵的快速求解.三维模型正演响应表明,基于二次场的三维大地电磁有限差分算法具有较高的计算精度和可靠性.通过计算分析不同PML吸收因子条件的大地电磁正演结果,显示在适当的吸收因子下,PML吸收边界条件可较大幅度的减小外边界距离,从而有效的压缩模型求解空间,最终提高三维大地电磁正演模拟的效率. 展开更多
关键词 大地电磁法 二次场 有限差分交错网格法 pml吸收边界
在线阅读 下载PDF
PML边界条件下孔隙介质弹性波可变网格正演模拟方法研究 被引量:6
4
作者 孙林洁 刘春成 张世鑫 《石油物探》 EI CSCD 北大核心 2015年第6期652-664,共13页
孔隙介质弹性波正演模拟对于储层预测和油气检测具有重要的指导意义。为提高孔隙介质波场正演模拟的精度,将可变网格方法应用于基于Biot理论的双相介质弹性波高阶有限差分正演模拟方法中,有效解决了二阶混合偏导数差分近似的可变网格过... 孔隙介质弹性波正演模拟对于储层预测和油气检测具有重要的指导意义。为提高孔隙介质波场正演模拟的精度,将可变网格方法应用于基于Biot理论的双相介质弹性波高阶有限差分正演模拟方法中,有效解决了二阶混合偏导数差分近似的可变网格过渡问题,推导了基于PML边界条件的可变网格差分格式,实现了空间与时间步长任意整数倍变化的双相介质弹性波可变网格正演模拟。对于细微地质结构或横向变化剧烈的介质区域,利用可变网格对区域精细重采样可有效提高正演的模拟精度;对于介质相对稳定的区域,采用大步长采样可保证正演计算效率。数值试验结果表明,与常规正演模拟方法相比,双相介质可变网格正演模拟方法的精度、效率有明显提高,对精细研究孔隙介质地震响应有着重要意义。 展开更多
关键词 BIOT理论 孔隙介质 pml边界条件 可变网格 有限差分
在线阅读 下载PDF
稀疏存储的显式有限元三角网格地震波数值模拟及其PML吸收边界条件 被引量:30
5
作者 刘有山 滕吉文 +1 位作者 刘少林 徐涛 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2013年第9期3085-3099,共15页
有限元法是复杂介质地震模拟的有力工具,它能比较客观地反映地震波的传播,比较细致地再现地震图像.但是,为了获得较精确的结果,有限元法模拟地震波的传播需要的网格点数多,具有计算量大和消耗内存多的缺点.针对上述缺点,本文对刚度矩阵... 有限元法是复杂介质地震模拟的有力工具,它能比较客观地反映地震波的传播,比较细致地再现地震图像.但是,为了获得较精确的结果,有限元法模拟地震波的传播需要的网格点数多,具有计算量大和消耗内存多的缺点.针对上述缺点,本文对刚度矩阵采用压缩存储行(CSR)格式,以减少计算量并节省内存;采用集中质量矩阵得到对角的质量矩阵以提高有限元法(显式有限元)的计算效率;时间离散采用保能量的Newmark算法以提高有限元法的计算精度;采用变分形式(弱形式)的PML吸收边界条件对人工截断边界进行处理.通过与高精度的数值方法——谱元法的数值试验的对比表明,上述方法的引入可使有限元法在计算精度和计算效率方面均可取得比较显著的改进.为了获得相当的计算精度,相比于7阶谱元法,显式有限元法需要更精细的网格.然而,显式有限元法的计算速度比前者快近2倍,而内存需求仅为谱元法的1/4~1/6. 展开更多
关键词 显式有限元法 压缩存储行(CSR)格式 集中质量矩阵 Newmark算法 pml吸收边界条件 谱元法
在线阅读 下载PDF
PML吸收边界条件影响因素分析 被引量:8
6
作者 丁科 《物探与化探》 CAS CSCD 2012年第4期623-627,共5页
影响完全匹配层方法吸收效果的主要因素有吸收层厚度和衰减系数,笔者通过地震波场数值模拟,讨论了吸收层厚度和衰减系数对吸收效果的影响。研究表明,在衰减系数一定时,吸收层厚度越大,边界反射越弱。吸收层的厚度一般在12~20道较为适当... 影响完全匹配层方法吸收效果的主要因素有吸收层厚度和衰减系数,笔者通过地震波场数值模拟,讨论了吸收层厚度和衰减系数对吸收效果的影响。研究表明,在衰减系数一定时,吸收层厚度越大,边界反射越弱。吸收层的厚度一般在12~20道较为适当,在吸收层的厚度比较小时,随着衰减系数的增大,边界反射逐渐减弱甚至完全没有边界反射,但是若再进一步增大衰减系数则又会逐步出现边界反射。因此,在实际应用中应该注意衰减系数的选择,衰减系数一般在500~2 000之间,衰减系数较大反而影响其效果。 展开更多
关键词 pml吸收边界 数值模拟 影响因素 波动方程
在线阅读 下载PDF
PML边界条件在空腔噪声计算中的应用
7
作者 陈维 《战术导弹技术》 北大核心 2018年第4期18-24,83,共8页
通过数值模拟气动噪声对空腔进行外形优化,以达到降低空腔气动噪声的目的。为准确计算流场声压脉动,数值方法采用了大涡模拟(LES)技术,远场边界采用了完全匹配层边界条件(PML),空间格式采用色散相关保持格式(DRP),同时时间格式采用低耗... 通过数值模拟气动噪声对空腔进行外形优化,以达到降低空腔气动噪声的目的。为准确计算流场声压脉动,数值方法采用了大涡模拟(LES)技术,远场边界采用了完全匹配层边界条件(PML),空间格式采用色散相关保持格式(DRP),同时时间格式采用低耗散低色散5级Runge-Kutta(LDDRK5)格式以使得流场耗散和色散最小。通过对三种外形的空腔在来流马赫数0.41时监测点以及脉动流场的分析和判断,采用PML边界条件对内部流场的反射误差很小,可以稳定的吸收出流波动。结果表明,该计算方法可以较好的计算空腔的声压,给出气动噪声最小的最优外形,为空腔降噪优化和气动噪声的计算提供了参考。 展开更多
关键词 空腔 pml边界条件 气动噪声 DRP格式 LDDRK
原文传递
PML吸收边界条件在波导不连续性问题中的应用 被引量:1
8
作者 邢立丹 徐立 +1 位作者 朱小芳 李斌 《真空电子技术》 2010年第5期25-29,共5页
采用矢量有限元方法结合PML吸收边界条件分析波导不连续性问题。利用PML吸收边界条件截断计算区域,可以减少内存消耗;采用模式正交性计算波导结构的S参数,可以获得更加精确的解。在结果验证阶段,分别对PML参数设置、厚度选择、波导介质... 采用矢量有限元方法结合PML吸收边界条件分析波导不连续性问题。利用PML吸收边界条件截断计算区域,可以减少内存消耗;采用模式正交性计算波导结构的S参数,可以获得更加精确的解。在结果验证阶段,分别对PML参数设置、厚度选择、波导介质不连续性以及波导结构不连续性模型进行分析,将计算结果与文献对比,验证了PML吸收边界条件的有限元法分析波导不连续性问题的准确性。 展开更多
关键词 有限元 pml吸收边界条件 波导不连续性 S参数
在线阅读 下载PDF
时变媒质边界条件与局部吸收负载的设置与应用 被引量:4
9
作者 任武 刘波 +1 位作者 高本庆 杨仕明 《微波学报》 CSCD 北大核心 2003年第3期1-6,共6页
首先对媒质匹配条件进行了分析 ,得出时变媒质中电磁波无反射传播所需的参数关系 ,并将普通的Mur和PML边界拓展到时变媒质的FDTD分析中 ;根据时变媒质的特性 ,提出了由时变媒质构造的局部吸收负载 ,并对其应用进行了探讨。结果表明 ,此... 首先对媒质匹配条件进行了分析 ,得出时变媒质中电磁波无反射传播所需的参数关系 ,并将普通的Mur和PML边界拓展到时变媒质的FDTD分析中 ;根据时变媒质的特性 ,提出了由时变媒质构造的局部吸收负载 ,并对其应用进行了探讨。结果表明 ,此吸收负载具有灵活多变、简单便捷。 展开更多
关键词 时变媒质 媒质匹配条件 pml边界 FDTD 局部吸收负载 电磁场
在线阅读 下载PDF
交错网格有限差分正演模拟的联合吸收边界 被引量:7
10
作者 胡建林 宋维琪 +2 位作者 张建坤 邢文军 徐文会 《石油地球物理勘探》 EI CSCD 北大核心 2018年第5期914-920,共7页
三维声波方程交错网格有限差分正演模拟中的边界问题一直是热点问题。完全匹配层吸收边界(PML)具有较强且稳定的吸收效果,但必须具有一定的边界厚度才能吸收干净,这就增大了三维正演模拟的模型空间,即增加了运算量;Higdon边界能消除任... 三维声波方程交错网格有限差分正演模拟中的边界问题一直是热点问题。完全匹配层吸收边界(PML)具有较强且稳定的吸收效果,但必须具有一定的边界厚度才能吸收干净,这就增大了三维正演模拟的模型空间,即增加了运算量;Higdon边界能消除任意角度入射波的边界反射,也具有较强稳定性,但该高阶吸收边界离散化后过于复杂,而低阶时吸收效果不如PML边界。因此,基于对PML吸收层中的平面波传播规律的研究,重新推导PML最外层的Higdon吸收边界条件,得到含PML吸收系数的新的Higdon吸收边界条件。联合吸收边界不仅可使用较小厚度(相对于单纯PML边界)的PML层对分量进行衰减,而且在PML边界外层,能应用新推导的Higdon吸收边界条件对反射波进行匹配吸收。在相同吸收效果下,联合吸收边界大幅度降低了PML厚度,减小了运算量,得到精确的模拟结果。 展开更多
关键词 三维声波方程 交错网格有限差分 正演 pml边界 Higdon边界 联合吸收边界
在线阅读 下载PDF
VTI介质一阶准P波方程正演模拟及边界条件 被引量:10
11
作者 韩令贺 何兵寿 《石油地球物理勘探》 EI CSCD 北大核心 2010年第6期819-825,共7页
本文根据具有垂直对称轴的横向各向同性(VTI)介质中的一阶准P波方程,构造了该方程在交错网格中求解的高阶有限差分格式,给出了其稳定性条件;重点研究了适用于该方程的完全匹配层(PML)吸收边界条件,在此基础上实现了VTI介质中一阶准P波... 本文根据具有垂直对称轴的横向各向同性(VTI)介质中的一阶准P波方程,构造了该方程在交错网格中求解的高阶有限差分格式,给出了其稳定性条件;重点研究了适用于该方程的完全匹配层(PML)吸收边界条件,在此基础上实现了VTI介质中一阶准P波方程的正演模拟,得到了高精度的波场快照和合成记录,并采用叠前逆时深度偏移算法对复杂模型的合成记录进行偏移成像,得到了较好的成像效果。数值算例表明:准P波方程正演模拟能克服常规弹性波方程正演模拟的某些缺陷,便于研究各向异性介质中纵波的传播机理;叠前逆时深度偏移算法能够对复杂模型准确成像;通过对比分析各向同性与各向异性Marmousi模型的合成记录以及两种偏移算法的成像结果,显示各向异性参数对波场特征(旅行时和振幅)会产生明显影响,在地震资料处理、反演和解释中不可忽略。 展开更多
关键词 VTI介质 准P波方程 正演 有限差分 交错网格 逆时偏移 pml吸收边界条件
在线阅读 下载PDF
基于PML井间地震纵横波分离的弹性波数值模拟
12
作者 李洋森 王涛 +1 位作者 孟凡顺 张亮 《物探化探计算技术》 CAS CSCD 2012年第6期735-740,624,共6页
井间地震方法观测到的波场信息丰富且复杂,其难点是如何识别和分离这些波场信息。这里从纵波波场为无旋场,横波波场为无散场出发,利用一阶速度~应力方程的交错网格高阶有限差分法,对井间弹性波纵、横波场分离问题进行了数值模拟。经数... 井间地震方法观测到的波场信息丰富且复杂,其难点是如何识别和分离这些波场信息。这里从纵波波场为无旋场,横波波场为无散场出发,利用一阶速度~应力方程的交错网格高阶有限差分法,对井间弹性波纵、横波场分离问题进行了数值模拟。经数值模拟试验表明,同二阶的位移方程相比,使用等价的波场分离交错网格的一阶速度应力方程,无论是模拟精度还是模拟的稳定性方面都有了较大的提高,同时该横波资料的获取方法为井间地震资料的解释应用,提供了一种新的技术思路。 展开更多
关键词 井间地震 波场分离 交错网格 pml吸收边界条件
在线阅读 下载PDF
褶积微分算子法数值模拟的边界问题研究
13
作者 贺同江 刘红艳 李小凡 《物探化探计算技术》 CAS CSCD 2010年第6期571-577,共7页
褶积微分算子法是一种全新的数值模拟方法,已被广泛应用于复杂介质的地震波场数值模拟,但是其边界反射问题一直没有解决。这里将最佳匹配层(PML)吸收边界条件引入到褶积微分算子法中,此方法是在研究区域的边界上加入吸收层,使边界上传... 褶积微分算子法是一种全新的数值模拟方法,已被广泛应用于复杂介质的地震波场数值模拟,但是其边界反射问题一直没有解决。这里将最佳匹配层(PML)吸收边界条件引入到褶积微分算子法中,此方法是在研究区域的边界上加入吸收层,使边界上传入吸收层的波,随传播距离按指数规律衰减,不产生任何反射,以达到消除边界反射的目的。构造不同的模型,通过对比分析证明,PML吸收边界条件能比较好地解决褶积微分算子法的边界问题,从而验证了完全匹配层吸收效果的优越性。 展开更多
关键词 褶积微分算子 数值模拟 pml吸收边界
在线阅读 下载PDF
色散媒质中的完全匹配层吸收边界
14
作者 李莉 高攸纲 《电波科学学报》 EI CSCD 1999年第1期1-6,共6页
将A.P.Zhao提出的与媒质无关的完全匹配层吸收边界推广到色散媒质中,推导出了用于色散媒质的无反射完全匹配条件。并对导电媒质进行了数值实验,结果十分理想。与Q.H.Liu的方法相比,该方法可以应用到更复杂的媒质中。... 将A.P.Zhao提出的与媒质无关的完全匹配层吸收边界推广到色散媒质中,推导出了用于色散媒质的无反射完全匹配条件。并对导电媒质进行了数值实验,结果十分理想。与Q.H.Liu的方法相比,该方法可以应用到更复杂的媒质中。而与S.D.Gedney的各向异性媒质吸收层相比,在角区的处理要简单得多,更适合于时域有限差分法。对于三维问题更是如此。 展开更多
关键词 时域有限差分法 pml吸收边界 色散媒质
在线阅读 下载PDF
Acoustic reflection well logging modeling using the frequency-domain finite-element method with a hybrid PML 被引量:6
15
作者 Wang Bing Kuo Zhang +2 位作者 Guo Tao He Liu Zhang Xiao-Liang 《Applied Geophysics》 SCIE CSCD 2018年第1期35-45,148,共12页
In this paper, we propose a hybrid PML (H-PML) combining the normal absorption factor of convolutional PML (C-PML) with tangential absorption factor of Mutiaxial PML (M-PML). The H-PML boundary conditions can be... In this paper, we propose a hybrid PML (H-PML) combining the normal absorption factor of convolutional PML (C-PML) with tangential absorption factor of Mutiaxial PML (M-PML). The H-PML boundary conditions can better suppress the numerical instability in some extreme models, and the computational speed of finite-element method and the dynamic range are greatly increased using this HPML. We use the finite-element method with a hybrid PML to model the acoustic reflection of the interface when wireline and well logging while drilling (LWD), in a formation with a reflector outside the borehole. The simulation results suggests that the PS- and SP- reflected waves arrive at the same time when the inclination between the well and the outer interface is zero, and the difference in arrival times increases with increasing dip angle. When there are fractures outside the well, the reflection signal is clearer in the subsequent reflection waves and may be used to identify the fractured zone. The difference between the dominant wavelength and the model scale shows that LWD reflection logging data are of higher resolution and quality than wireline acoustic reflection logging. 展开更多
关键词 Frequency domain finite-element method hybrid pml boundary wireline well logging logging while drilling
在线阅读 下载PDF
VSP逆时偏移及其存储策略研究 被引量:5
16
作者 王维红 郭雪豹 +1 位作者 石颖 刘诗竹 《石油地球物理勘探》 EI CSCD 北大核心 2015年第1期77-83,16,共7页
在逆时偏移算法中,应用常规的随机边界方法虽可节约存储空间,但浅层常伴有随机噪声。应用PML边界方法虽可改善上述状况,却又面临地震波场存储的压力。为此,本文采用优化系数的高阶有限差分方法实现VSP数据逆时深度偏移,采用拉普拉斯去... 在逆时偏移算法中,应用常规的随机边界方法虽可节约存储空间,但浅层常伴有随机噪声。应用PML边界方法虽可改善上述状况,却又面临地震波场存储的压力。为此,本文采用优化系数的高阶有限差分方法实现VSP数据逆时深度偏移,采用拉普拉斯去噪方法压制低频噪声,并兼顾考虑精度和存储,在PML边界震源波场正传过程中保存部分波场,进而利用保存的信息与检波点波场同步反传,不仅可有效地节约存储空间,也确保了替代波场信息的可靠性。断层模型测试表明,本文方法能够以低存储实现高精度的VSP逆时偏移,相比于地面地震偏移,断层成像更清晰、准确。 展开更多
关键词 VSP 逆时偏移 随机边界 pml边界
在线阅读 下载PDF
基于双相介质BISQ模型的地震波正演模拟 被引量:4
17
作者 胡英 陈辉 +2 位作者 苏云 刘忠亮 张红霖 《石油地球物理勘探》 EI CSCD 北大核心 2012年第6期901-907,1024+839,共7页
Biot流动机制和喷射流动机制是含流体孔隙介质中的两种重要的力学机制,同时包含这两种机制的BISQ模型能够更好地反映孔隙介质中弹性波的传播规律。本文从BISQ模型的基本方程出发,推导了完全匹配层(PML)吸收边界的交错网格高阶有限差分... Biot流动机制和喷射流动机制是含流体孔隙介质中的两种重要的力学机制,同时包含这两种机制的BISQ模型能够更好地反映孔隙介质中弹性波的传播规律。本文从BISQ模型的基本方程出发,推导了完全匹配层(PML)吸收边界的交错网格高阶有限差分方法的表达式,对双相各向同性介质中地震波场进行了数值模拟;然后通过调整孔隙度、渗透率等储层参数,对比分析了储层参数对双相介质波场特征的影响。数值模拟结果表明,基于BISQ模型的正演模拟正确反映了双相各向同性介质的波场特征;储层参数取值变化对波场特征的影响也十分明显。 展开更多
关键词 BISQ模型 双相介质 pml边界条件 波场模拟
在线阅读 下载PDF
基于高阶谱差分的CAA模型预测高速列车过隧微气压波 被引量:3
18
作者 杨志刚 谭晓明 +1 位作者 梁习锋 任鑫 《铁道学报》 EI CAS CSCD 北大核心 2014年第7期85-89,共5页
高速列车以大于300 km/h的速度穿越隧道时,在隧道出口形成的微气压波具有较强的非线性。该微气压波幅值的数量级达到102 Pa,难以用传统的CFD方法(2阶精度)精确模拟。基于高阶谱差分的CAA模型是一种值得尝试的方法。CAA方法能在强非线... 高速列车以大于300 km/h的速度穿越隧道时,在隧道出口形成的微气压波具有较强的非线性。该微气压波幅值的数量级达到102 Pa,难以用传统的CFD方法(2阶精度)精确模拟。基于高阶谱差分的CAA模型是一种值得尝试的方法。CAA方法能在强非线性区域保持流体变量和声学变量的传递精度,高阶谱差分方法能够融入非结构网格,获得较高的计算效率和几何灵活性。本文在动模型试验平台开展200~350 km/h速度下的微气压波缩比试验研究;采用基于3阶谱差分的CAA方法,对轴对称圆管模型,数值研究了管口气压爆波的产生和近场辐射特性。数值结果和缩比试验结果的一致性,验证了该数值方法的准确性和可靠性。 展开更多
关键词 微气压波 pml边界条件 高阶谱差分 计算气动声学 非线性特性
在线阅读 下载PDF
基于高阶有限差分纵横波分解的弹性波数值模拟 被引量:4
19
作者 赵茂强 吴国忱 《物探化探计算技术》 CAS CSCD 2010年第5期487-494,共8页
常规弹性波数值模拟方法得到的VSP波场信息中,P波和S波是相互耦合在一起的,波场信息复杂,在进行多波多分量地震波场特征分析时需要进行波场分离。这里采用纵波、横波解耦的弹性波波动方程高阶有限差分算法进行数值模拟,并采用PML边界条... 常规弹性波数值模拟方法得到的VSP波场信息中,P波和S波是相互耦合在一起的,波场信息复杂,在进行多波多分量地震波场特征分析时需要进行波场分离。这里采用纵波、横波解耦的弹性波波动方程高阶有限差分算法进行数值模拟,并采用PML边界条件处理人工边界条件,得到纵波、横波解耦的独立VSP纵波、横波记录,为进一步研究弹性波VSP波场特征奠定了基础。 展开更多
关键词 弹性波 高阶有限差分 纵横波解耦 pml边界
在线阅读 下载PDF
基于LS-RSGFD方法优化的横向各向同性(TI)介质一阶qP波高精度数值模拟 被引量:12
20
作者 李振春 杨富森 王小丹 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2016年第4期1477-1490,共14页
为克服各向异性介质弹性波数值模拟中存在着计算量大和波场分离困难等局限,研究了声学近似的VTI介质和TTI介质一阶qP波数值模拟方法.首先对VTI介质弹性波方程进行声学近似,推导了VTI介质一阶qP波方程;然后基于精确的TTI介质频散关系,引... 为克服各向异性介质弹性波数值模拟中存在着计算量大和波场分离困难等局限,研究了声学近似的VTI介质和TTI介质一阶qP波数值模拟方法.首先对VTI介质弹性波方程进行声学近似,推导了VTI介质一阶qP波方程;然后基于精确的TTI介质频散关系,引入一个包含各向异性控制参数σ的新辅助波场,推导了稳定的TTI介质二阶耦合qP波波动方程,并通过引入波场的伪速度分量,推导了等价的一阶应力-速度形式.结合旋转交错网格有限差分(RSGFD)和基于最小二乘优化的有限差分(LS-FD)两种各具优势的方法,研究了最小二乘旋转交错网格有限差分(LS-RSGFD)方法,并用其数值求解VTI和TTI介质一阶qP波方程,然后通过构造其LS-RSGFD格式,实现了高精度的各向异性介质qP波波场数值模拟.数值模拟结果表明:TI介质一阶qP波方程能够准确地模拟各向异性介质中qP波的运动学特征,引入控制参数σ能够有效地减弱不稳定性问题,保证非均匀TTI介质中qP波场的稳定传播;利用优化的LS-RSGFD方法可以得到高精度的合成地震记录,同时还可以相对地提高计算效率. 展开更多
关键词 TI介质 LS-RSGFD方法 一阶qP波 数值模拟 pml边界条件 控制参数 不稳定性
在线阅读 下载PDF
上一页 1 2 4 下一页 到第
使用帮助 返回顶部