垂向非均匀介质中转换波转换点计算

转换点位置的计算是转换波资料处理中的一个关键问题. 本文提出了分别基于速度随深度线性变化、速度随垂直走时线性变化、慢度随深度线性变化和慢度随垂直走时线性变化四种等效垂向非均匀介质情况下转换点位置的计算方法. 研究了通过速度拟合、走时近似和相似系数谱三种方式选择合适的等效速度方法. 结合理论模型对非均匀介质转换点计算方法、渐进转换点计算方法、Thomsen近似公式和均匀介质解析计算方法的误差进行了分析,结果表明非均匀介质转换点计算方法能更准确地计算转换点位置.     

波形差异反演方法及其在气藏描述中的应用

在油气藏开发中后期,如何利用的已有的勘探成果和丰富的油气藏知识,加深对勘探区域油气开发潜力的认识,为油气田的稳产高产服务,是摆在油气田中后期开发人员面前的难题.油气田开发的研究人员,利用现有的资料,开发出了很多方法,结合地震测井等资料为油气田的稳产高产做出了巨大的贡献.本文在分析了地震波在地质体中传播过程中能量变化,相应的产生了振幅变化的理论,利用三维地震数据,通过对三维地震数据解释层位对应的波形差异的分析,结合工区内的生产井和开发井的对应地层中的油气信息,进行反演得到地层的波形差异数据体,利用相应的绘图软件绘制二维平面图形,结合工区内的沉积相等地质信息对其数据体进行油藏解释,最终得到波形差异反演数据体,利用波形差异反演数据体进行油气藏的识别、描述和预测,为油气田的中后期开发服务.     

完全匹配层吸收边界在孔隙介质弹性波模拟中的应用

模拟弹性波在孔隙介质中传播,需要稳定有效的吸收边界来消除或尽可能的减小由人工边界引起的虚假反射. 本文在前人工作基础上,首次建立了弹性孔隙介质情况下完全匹配层吸收边界的高阶速度-应力交错网格有限差分算法,并详细讨论了完全匹配层的构建及其有限差分算法实现. 首先,本文通过均匀孔隙模型的数值解与解析解的对比,验证所提出的数值方法的正确性;然后,本文考察了完全匹配层对不同入射角度入射波和自由表面上的瑞利波的吸收性能,将完全匹配层与廖氏和阻尼吸收边界进行了对比,研究了这三种吸收边界在不同吸收厚度情况下对弹性波吸收能力. 数值结果表明,在孔隙介质中,完全匹配层作为吸收边界能十分有效地吸收衰减外行波,无论对体波还是面波,是一种高效边界吸收算法.     

植被电磁散射的半空间模型研究

本文给出了地面植被电磁散射的半空间模型研究方法.在以往的相关文献中,均采用自由空间格林函数求解单个叶片散射体的散射场,本文利用半空间并矢格林函数求出了单个散射体的散射场,然后运用Monte Carlo方法模拟生成地面植被层,得到了在半空间下植被层的电磁散射特性,计算出了单、双站雷达散射截面,并与自由空间下的散射场做了相互对比.结果表明,在半空间格林函数下得到的散射场更为适用于描述地面植被的电磁散射特性.     

局部角度域波传播步进算法研究

本文从非均匀介质中波动方程出发,提出了基于一般标架的相空间（局部角度域）波传播的步进算法. 该方法在构造单程波的步进算法时,在选择标架或正交基等方面有更大的自由度. 我们以不随频率及深度变化且具有变尺度特性的Gabor_Daubechies紧标架为例,给出了单程波传输算子的具体形式及相应的波场步进算法;详细讨论了基于Gabor_Daubechies标架的传输算子的高频渐近展开问题,得出了在高频、小传输步长条件下传输算子的近似解析表达式,并给出使用条件. 通过模型算例,比较了精确传输算子与高频近似传输算子用于非均匀介质中波传播的结果,说明在一定条件下由两者得出的波场几乎是相同的.     

天体运行轨道的背景介质理论导引与自相似分形测度计算的分维微积分基础

通过讨论天体运行背景介质理论的连续轨道及离散轨道这二个研究方向的基础假设,介绍了天体运行轨道的具体方程形式及理论框架概要;进一步地通过讨论天体运行轨道Binet方程的一般形式及其行星近日点进动角的解,给出了连续轨道理论与Newton理论及Einstein广义相对论的联系与区别;通过讨论天体运行轨道的分维扩展方程,给出了包括太阳系行星、天王星卫星、地球卫星、绕月航天器等在内的离散轨道（稳定性轨道）方程及其预言数据.特别地,作为对天体在较为广泛区域作用曲线的初步探讨推论,指出仅由天体引力难以形成质量密度趋于无穷大的理想黑洞.通过讨论一般函数的分维导数的位置假设及幂函数的分维导数的形式假设,进一步明晰了幂函数的分维导数、分维微分及分维积分的具体方程形式,给出分维导数与分数阶导数的区别,随后讨论了基于一般分形测度的分维微积分形式定义导出的自相似分形的测度计算方程具体形式,给出了其与目前Hausdorff测度方法（覆盖方法）的区别,并对包括三分Cantor集合、Koch曲线、Sierpinski垫片及正交十字星形等自相似分形在内的测度进行了计算分析.     

VTI介质qP波方程高精度有限差分算子

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

非线性波消波及波场分布研究

本文基于雷诺平均N-S方程,并结合RNG k-e方程建立了粘性数值波浪水槽,对不同波陡、不同相对水深、不同相对波高的非线性规则波的阻尼消波问题和波场分布进行研究。文中提出了两种描述消波区内部阻尼变化的阻尼函数,分别适用于小波陡情形和高波陡情形。研究结果表明,小波陡组消波区可设为一个波长,阻尼系数取10~4～10~5即可满足消波要求,计算结果与实验结果及造波理论吻合良好;高波陡组消波区可设为两个波长,阻尼系数取10~4～10~5亦可满足消波要求,计算结果与实验结果吻合良好。此外,当波陡较小时,波场内反射情况的小幅改变即可对整个波场造成影响,特别是当水深较浅时这种影响极为明显,需谨慎考虑。当波陡较大时,水波能量较高,整个波场沿水波传播方向可观测到明显的衰减现象,在具体试验中需进行考虑。     

波浪作用下海床瞬态剪切破坏深度的计算分析

鉴于弹性本构模型的简便性,目前波浪作用下海床瞬态剪切破坏深度的计算一般都是基于海床的多孔弹性介质模型,这种方法不能考虑塑性单元不平衡应力的重分配,而如果将海床视为多孔理想弹塑性介质进行波浪作用下的瞬态响应分析,则可以克服上述缺陷。为此,以波浪作用下砂质、粉砂质海床瞬态响应模型为例,基于Mohr-Coulomb屈服准则,对不同海床计算参数条件下上述2种计算方法的结果进行了对比分析。结果表明,将海床视为多孔理想弹塑性介质时海床瞬态剪切破坏深度大于等于海床为多孔弹性介质时的计算结果;随着海床计算参数的变化,2种方法计算结果的变化趋势较为一致,且在砂质、粉砂质海床常见计算参数范围内,2种方法计算结果的差值在一个相对固定的范围内变化。