储层模型正演模拟及AVO特征研究

为了在多波地震资料解释过程中解决多波信息识别难的问题,基于弹性波正演技术,针对某油田勘探区块的实际储层地质模型,采用高阶交错网格差分技术正演模拟非零偏炮集多波地震记录,并经多波资料处理,得到纵波和转换波CMP和CCP叠前含AVA的动校正道集及叠加剖面。在已知模型的AVA计算、模型层位参数与正演纵、横波波场关系的基础上,对比研究了正演记录中各储层的纵波和转换波对应关系。对比分析结果表明,数值模拟结果与AVA分析结果吻合,同一储层的纵波和转换波波场响应有明显差异。     

基于惠更斯原理的波动方程共炮点道集地震正演

地震非零炮检距正演模拟技术在叠前地震资料处理和解释中有重要意义.作者利用下行波将震源场延拓到地质介质任意深度,实现单程波动方程模拟叠前正演记录.用单程波动方程法对一个盐丘地质模型生成129个共炮点地震记录,并对所有共炮记录进行叠前深度偏移.在共炮记录中,来自盐丘地质模型的反射波和绕射波十分清晰.叠前深度偏移剖面和盐丘地质模型完全吻合.该方法生成记录质量好且计算效率高,并且叠前正演记录没有多次波、直达波和其他干扰波.     

任意裂隙方位双相HTI介质的地震波方程及正演模拟

基于BISQ方程的高频极限,导出了适用于任意方位角的双相HTI介质波动方程,推导了该方程在交错网格空间中求解的高阶有限差分格式和对应的完全匹配层(PML)吸收边界条件,实现了该类介质的地震波场正演模拟。模拟结果表明,该方法能准确模拟地震波在双相HTI介质中的传播过程,得到高精度的波场快照和合成记录。分析得出:在双相各向异性介质中,弹性波的传播仍表现出各向异性特征;孔隙度、粘滞系数和固液耦合密度主要影响慢纵波的振幅和形态。     

方位各向异性介质的多尺度有限差分法波场模拟

采用以裂缝走向为方位角的方位各向异性介质模型可以较好地描述裂缝性油气藏的实际情况,对弹性波在此介质中的传播过程进行准确的数值模拟有助于提高油气开采的准确程度。该文采用紧支集正交小波基对空间域进行多尺度离散,采用二阶精度有限差分算子对时域离散,推导得到了多尺度有限差分方法正演模拟的递推公式,并实现了相应的波传过程数值模拟。数值结果准确地反映了方位各向异性介质中波场的变化过程,可以清晰地观察到横波分裂和方位特征差异等现象。     

基于多GPU加速的各向异性弹性波正演模拟

比较分析了在不同网格大小介质模型情况下,分别采用串行计算、CPU 16个线程并行计算和4块GPU并行计算进行各向异性弹性波动方程正演模拟的执行时间差异。发现在网格点为2563的大模型上,用4块GPU的并行模拟计算相对16线程并行计算与串行计算的加速分别为30倍与156倍。表明多GPU并行算法可以显著缩短数值模拟时间,而且模型网格越大,加速效果越显著。因此,在单机环境下进行大尺度模型的各向异性弹性波正演模拟,采用多GPU并行计算方式是一个合适的加速选择。     

基于时频误差分析法的随机介质波场特征分析

设计了含裂缝的随机介质模型和含随机介质的复杂盐丘构造模型,基于声波方程有限差分法,进行了波场快照及共炮点地震记录模拟;利用连续小波变换的时频误差分析方法,分析了与所给模型相同构造的各向同性介质模型的地震记录的时频振幅误差;统计分析了各道地震记录的时频特性,所得结果与利用连续小波变换的时频误差分析法吻合较好.表明时频误差分析法能为研究随机介质波场特征提供一种可行的思路.     

三维各向异性介质中弹性波方程交错网格高阶有限差分法数值模拟

三维数值模拟是研究各向异性介质中复杂弹性波波场特征和传播规律的重要手段。根据Taylor展开式,推导出了交错网格一阶空间导数的任意偶数阶精度有限差分近似式和相应的差分系数计算式,给出了三维各向异性介质中弹性波一阶双曲型应力速度方程交错网格任意偶数阶精度差分格式和稳定性条件,并推导出了三维各向异性介质PML法吸收层边界条件公式和相应的交错网格差分格式。对方位各向异性介质模型和正交各向异性介质模型中弹性波的传播进行的三维数值模拟结果表明,弹性波在三维各向异性介质中传播时存在拟P波、拟SV波和拟SH波,并出现横波分裂、横波分裂盲区、波面三分叉等特殊现象,另外,弹性波场在空间是变化的,其拟P波、拟SV波和拟SH波的耦合关系比较复杂。     

隧道反射波超前探测有限差分正演模拟与偏移处理

针对地下工程领域隧道不良地质灾害超前探测问题,构建断层、倾斜软弱夹层等地质灾害体的隧道介质模型,研究隧道反射波超前探测的波场传播机理与偏移成像问题。利用一阶速度-应力弹性波方程和高阶交错网格差分计算方法,导出隧道工程反射波探测数值模拟的差分计算格式、数值稳定性条件以及边界条件,对上述模型进行正演模拟,得到相应模型不同时刻的波场快照与地震合成记录;通过对这些模拟结果进行分析,验证高阶交错网格差分法在隧道波场正演中的有效性。在导出弹性波方程逆时延拓交错网格差分计算格式基础上,通过求解程函数方程,获取激发时间成像条件,进而得到隧道全波场逆时偏移算法及偏移处理程序,并应用于隧道波场偏移成像中。研究结果表明:采用高阶交错网格有限差分逆时偏移算法能使反射波归位,绕射波收敛,从而大大提高了隧道工程反射波探测的分辨率,为隧道超前探测资料的解释提供理论依据。     

弹性波叠前逆时深度偏移技术及其应用

研究了二维各向同性介质中两分量记录的基于弹性波标量波动方程的叠前逆时深度偏移方法.从弹性波波动方程出发,推导各向同性介质中弹性波逆时深度偏移技术中相关算法,并探讨其中有关问题的解决方案,实现二维多分量资料的叠前逆时深度偏移,最后进行模型试算和实际资料试处理.结果表明,叠前逆时深度偏移方法不受地层倾角限制,较适用于高陡构造地区或介质横向速度变化较大地区的多分量地震资料处理.     

三维各向异性介质中弹性波方程交错网格高阶有限差分法数值模拟

三维数值模拟是研究各向异性介质中复杂弹性波波场特征和传播规律的重要手段。根据Taylor展开式，推导出了交错网格一阶空间导数的任意偶数阶精度有限差分近似式和相应的差分系数计算式，给出了三维各向异性介质中弹性波一阶双曲型应力-速度方程交错网格任意偶数阶精度差分格式和稳定性条件，并推导出了三维各向异性介质PML法吸收层边界条件公式和相应的交错网格差分格式。对方位各向异性介质模型和正交各向异性介质模型中弹性波的传播进行的三维数值模拟结果表明，弹性波在三维各向异性介质中传播时存在拟P波、拟SV波和拟SH波，并出现横波分裂、横波分裂盲区、波面三分叉等特殊现象，另外，弹性波场在空间是变化的，其拟P波、拟SV波和拟SH波的耦合关系比较复杂。     

K(n,-n,2n)方程的行波解

利用动力系统分支理论和定性理论研究了$K(n,-n,2n)$方程的行波解及其动力学性质. 结合可积系统的特点, 得到系统的孤立行波解,不可数无穷多光滑周期行波解和不光滑行波解;并根据行波解与相轨线间关系,揭示了不同类型行波解间转变与参数变化的关系.     

随机波浪对浪溅区三维结构物冲击作用试验研究

通过物理模型试验研究了随机波浪对位于浪溅区三维结构物的冲击作用.试验波浪为单向随机波.试验中波浪分别从正向(波向角β＝0°)和斜向(β＝15°、30°、45°)入射,入射波波高H1/3=7.5～20 cm,谱峰周期Tp=0.75～2.0 s,模型相对净空s与入射波高H1/3之比为0～0.5.对不同试验组次试验数据的统计分析比较得出:正向波和斜向波作用于结构物底面的冲击压力特征值,中部大、周围小;波浪冲击压力等值线基本上垂直对称分布于波向线两侧;结构物所受波浪冲击力随入射波波高增大而增大,对应相对板宽的变化出现2个峰值,在相对净空s/H1/3为0.1和0.2时到达最大值.     

摇板式造波机性能实验研究

通过改进后的水池造波系统对工程中常用的规则波和二维非规则波浪的模拟试验，研究了该造波系统的最大造波能力、波浪总谐波变形以及对随机波浪谱的模拟精度等波浪模拟性能．测试数据表明，该造波机控制系统经过技术改进，恢复了二维规则波的模拟功能，提高了造波品质，完全满足船模耐波性试验的要求，可应用于实际工程项目之中。     

单一过程波浪变形和爬高试验研究

为计算浅水区和破碎带内波浪变形的数学模型提供验证数据,介绍了关于孤立波和单波群（这里统称为单一过程波浪）的变形和爬高的试验研究成果。这些数据是水槽试验无法得到的,试验结果有可能引起本领域内科研工作者们的兴趣。     

一个非线性Drinfeld-Sokolov系统的行波解分支

用动力系统分支理论研究一个非线性Drin feld-Sokolov系统,证明该系统存在扭子波、反扭子波、孤立波和无穷多光滑周期波解,获得在不同参数条件下上述解存在的充分条件及其精确表达式.     

波在非均匀媒质中传播时的解

该文从声场方程出发推导出波在非均匀媒质中传播时的波动方程,并利用WKB近似法得出声波在空气中传播时波函数的表达形式,为研究声波在非均匀媒质中的传播提供了一种近似解。     

广义水波方程组行波解的分支

应用动力系统分支理论,研究广义水波方程组行波解的分支.在固定的参数条件下给出广义水波方程组的孤立波、扭结(反扭结)波解的精确表达式,并证明该方程组存在不可数无穷多个周期波解.     

上海海域风浪谱分析与研究

根据上海海域四个测点的11300多且实测波面记录，经初谱估计和分析，筛选出拟合样本280多个，并求得各测站的标准化平均样本谱及4个测站的综合平均样本谱。结果表明上海海域的风浪谱都呈现双峰谱型，且其低频部分的能量比重约为15％－20％。因此，经过多种拟合计算分析后认为：该海域拟采用JONSWAP谱与三参数谱联合或者五参数谱与三参数谱联合进行拟合较为合适，同时提出了适用于上海海域的双峰谱型及其相应参数，可作为上海海域沿海建筑物的不规则波设计的备选谱型，对我国其他类似海域也有一定的参考价值。     

柱坐标系下修正弹性应力波理论及数值分析

指出了传统的波动理论存在的一些局限性和不足。给出了修正弹性应力波理论的建立推导过程。完成了波动方程从张量状态到柱坐标系下的转换过程,得到了柱坐标系下的波动方程。建立了波动变量在加载面上与载荷及速度的关系,应力波在自由面上发生反射时边界条件对波动变量的影响。通过使用Matlab编程计算了应力波在几何模型为z方向无限长的空心圆柱结构中的传播过程以及应力波在遇到边界时的反射情况。计算结果显示了应力波随时间在结构中的传播方式,应力波在加载面上会产生2个波,即体积波和形变波,遇到边界后会反射出2个波。同时体积波和形变波的第一部分复合在一起形成一个复合脉冲以相同的波速运动。     

半地转近似下正压大气中的非线性波动

本文在半地转近似下讨论了正压模式大气中的非线性 Rossby 波.结果表明:这些非线性 Rossby 波可以是椭圆函数形式的波、孤立波或冲击波;产生这些波动的条件是由一个称之为“总能量”的物理量给出的,同时还敏感地依赖于相应量的初相值;文中还求出了一支没有线性近似对应的新频率的非线性波解,它纯粹是非线性效应的产物.