简介:地震波场数值模拟方法对理解和分析地震波的传播规律具有着重要的意义。弹性波动方程能够模拟地下介质的实际情况,为偏移和成像提供有效的依据。在弹性波波场数值模拟中,旋转交错网格数值模拟(RSM)修改了标准交错网格数值模拟(SSM)方法,将同类的参数定义在同样的节点上,拓宽了稳定性条件的约束,但在低速区会出现较严重的频散。变阶数差分方法是自适应空间算子长度方法的一种变化和推广。它以理论频散误差研究为基础,结合实际波场传播的情况进行误差计算,对不同速度匹配不同的差分阶数。本文研究了变阶数旋转交错网格数值模拟(VRSM),即是籽变阶数方法应用到RSM中,它可以很好地解决RSM在低速区域的数值频散问题,以及减少不必要的时间损耗;同时讨论了旋转交错网格的理论频散特性,并基于波场分离的方法分析了实际波场传播的频散误差,将原方法的应用范围由声波推广到剪切波,由理论值推广到时变值。在数值模拟试验中,VRSM将被应用于水平层模型和Overthrust模型。通过阶数分配以及相应波场传播效果和计算时间的分析,验证了该方法应用于复杂介质波场模拟中的实用性和有效性。实验的结果表明VRSM能够合理分配不同速度所对应的差分阶数,能保证计算的精确性,并合理控制计算的时间。
简介:孔隙结构特征及类型划分对低渗透储层勘探开发至关重要,基于多重分形理论与核磁共振实验,对东营凹陷南坡沙河街组沙四段(ES4)复杂低渗透砂岩进行孔隙结构研究。首先,根据岩心物性、铸体薄片、压汞等资料所反映的孔隙结构参数差异,将研究区的岩石孔隙结构类型划分了3大类5小类;然后,针对不同类型岩石的核磁共振T2谱进行插值并计算其对数坐标下的一维、三维分形维数以及多重分形谱,并提取多重分形参数奇异性强度α、分布稠密度f(a),结果显示孔隙结构类型不同,盒维数尤其是多重分形参数值差异明显,孔隙结构好,其a、f(a)偏向高值,以此划分孔隙结构类型与压汞、薄片分析结果基本一致;最后,将该方法应用到核磁共振测井剖面上,应用效果较好,表明多重分形是核磁共振T2谱的一种属性,利用核磁测井T2谱多重分形特征及参数能够连续较好地评价低渗透砂岩孔隙结构类型与预测有效储层。
简介:有限差分方法广泛应用于求解许多科技领域所涉及的偏微分方程,高阶显式有限差分方法通常用来提高求解精度,已经提出的高阶隐式有限差分方法和截断高阶显式有限差分方法可用来进一步提高模拟精度而不增加计算量。本文首先计算了针对常规网格上的一阶导数和二阶导数、交错网格上的一阶导数的有限差分系数,发现高阶隐式有限差分系数中存在一些小的系数。频散分析结果表明:忽略这些小的差分系数能够近似维持有限差分的精度,但是显著减小了计算量。然后,引入镜像对称边界条件来提高隐式有限差分方法的精度和稳定性,采用混合吸收边界条件来减小来自模型边界所不需要的反射。最后,给出了针对均匀和非均匀介质模型的弹性波模拟例子,表明了本文方法的优点。