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种方法计算结果的差值在一个相对固定的范围内变化。     

基于粒径分维值检验的松散介质渗透系数空间估值——以河北平原为例

为了克服松散含水介质非均质性各向异性对介质渗透系数的求取所带来的困难,以水动力条件为线索,利用互相关法寻找同一时代地层及同一时代地表水的主、次流方向.在主流方向上,建立了石家庄-束鹿-衡水段介质颗粒粒径随距离递变的数学模型,并以此为基础推求无钻孔地段的粒径.同时,运用分形理论,计算了中、上更新统地层主流方向平均粒径空间分布的分维值,并将其作为检验粒径估值正确与否的判据.通过不同岩性的渗透系数-深度模型,实现了用平均粒径预测渗透系数的全过程.     

双重介质污染质运移模型研究及其在岩溶区的应用

本文用双重介质理论,在天然裂隙系统中对其裂隙域和孔隙岩块域建立了二维水流和污染质运移模型。讨论了岩块为斜方体、椭园柱体和任意不规则几何形态条件下,耦合项(Γ)的表达式及计算方法。分别对流速、弥散度、孔隙度等因素对污染质运移的控制和影响进行了定量分析。最后,模拟了济南东郊约85km~2岩溶含水层中Cl~-离子的运移情况,取得了比较满意的结果。     

准各向异性粘弹介质地震波的数字仿真

本文论述准各向异性粘弹波动方程有限元数值解地震波响应的原理与计算方法,并给出地震衰减因子与方程中粘滞系数的计算关系式。列出了所研制的地面记录地震剖面与井中记录垂直地震剖面,以及波场时间切片图、散度与旋度图、应变能与动能的时间切片图。本方法适于任意分层形态且层间物性参数为强间断的实际地震剖面的仿真,旨在研究正问题以实现反问题。