黏滞声学介质地震波吸收补偿叠前时间偏移方法

为了在叠前偏移中沿波场传播路径补偿黏滞性介质的吸收衰减效应,提高地震分辨率,首先导出黏滞声学介质条件下的波场频散关系,引入等效Q值实现频率域波场延拓,基于稳相原理求出波场延拓的渐近解,利用反褶积成像条件得到成像结果,并给出了偏移算法稳定性和数据高频折叠的处理方法,研制出工业化应用软件,完成多个三维工区的地震数据处理。另外,根据VSP资料,正演分析了松辽盆地的地层吸收对地震波场传播的衰减效应。与现有常规叠前时间偏移结果相比,黏滞声学介质吸收补偿地震叠前时间偏移在保护地震低频成分的前提下展宽频带15Hz左右,大幅提高了地震成像分辨率。     

基于格子Boltzmann方法的非常规天然气微尺度流动基础模型

基于相似准则和气体真实物性参数提出一种格子Boltzmann模型中选取无因次松弛时间的新方法,同时考虑边界努森层的影响对无因次松弛时间进行修正,结合壁面二阶滑移边界条件推导出相应组合反弹/镜面反射边界条件中的关键参数,从而建立了基于格子Boltzmann方法的适用于高温高压条件的非常规天然气微尺度流动模拟新模型.将甲烷...     

基于格子Boltzmann方法的3D数字岩心渗流特征分析

页岩气储层物性差、各向异性突出且孔隙结构存在多尺度特性,内部流动非常复杂,由于纳米孔隙的存在,宏观连续理论不适用。基于页岩二维扫描电镜图片和计算机处理技术,采用随机生长方法分别重建大孔隙和小孔隙的三维数字岩心,将大孔隙和小孔隙初始模型组合在一起构成多尺度的页岩气储层三维数字岩心组合模型,在此基础上,应用格子Boltzmann方法研究不同尺度孔隙的渗流行为。四参数随机方法重建的图像能较好体现原图的性质,储层的孔隙结构特征对其渗透率具有重要影响,渗透率随着孔隙度的增大而增大,组合模型数字岩心的渗透率大于大孔隙数字岩心和小孔隙数字岩心的渗透率,并且大于二者的总和,说明小孔隙的存在可以极大提高储层连通性。     

基于格子Boltzmann方法颗粒团聚物曳力系数的数值模拟及应用研究

应用气固两相流动的格子Boltzmann-离散颗粒运动模型,气体流动采用格子Boltzmann方法的D2Q9模型,固体颗粒运动采用牛顿第二运动定律计算,壁面采用具有二阶计算精度的半步长反弹边界条件,数值模拟了循环流化床提升管内颗粒团聚物的运动过程。模拟结果表明,颗粒团聚物的气固相间曳力系数随颗粒团聚物的空隙率和雷诺数的增大而减小;拟合得到基于格子Boltzmann方法的气固相间曳力模型;采用基于格子Boltzmann方法的气固相间曳力模型的欧拉-欧拉双流体模型模拟循环流化床提升管内颗粒相密度和颗粒质量流量径向分布,模拟结果与实验结果很接近。     

黏滞介质吸收补偿叠前时间偏移的倾角道集孔径优选

黏滞介质吸收补偿叠前偏移是提高地震成像分辨率的合理方法,考虑了地震波传播过程的高频损失,但在补偿地震高频能量的同时也放大了噪声,如何有效压制偏移噪声是实现高分辨率成像必须解决的问题。为此,基于稳相原理,提出了黏滞介质吸收补偿叠前时间偏移过程中生成倾角道集及自动拾取菲涅耳带的方法,认为倾角道集的同相轴以稳相点为中心向两侧上翘,存在明显的道集时差和振幅变化,同相轴集中弯曲的顶点为稳相点,稳相点邻域的同相轴近平直部分为菲涅耳带,即为最佳偏移孔径。经对拾取结果检查、修改、插值,形成时、空变偏移孔径。实际数据处理结果表明,与常规叠前时间偏移相比,变偏移孔径与黏滞介质吸收补偿叠前时间偏移相结合,具有保证陡倾反射成像、压制偏移噪声、改善成像分辨率和提高偏移计算效率的优势。     

声格固体模型地震波场模拟

ＰｅｔｅｒＭｏｒａ为模拟复杂介质中的纵波，提出了声格固体模型，即在格子Ｂｏｌｔｚｍａｎ模型基础上引入变速粒子而建立的细胞自动机模型，从而为利用这种方法模拟地震波场奠定了理论基础。文中简要介绍了Ｐｅｔｅｒ Ｍｏｒａ的声格固体模型的基本理论框架及差分算法，同时还介绿了黄联捷的插值算法。     

Numerical simulation of non-Archie electrophysical property of saturated rock with lattice Boltzmann method

The electrophysical property of saturated rocks is very important for reservoir identification and evaluation. In this paper, the lattice Boltzmann method （LBM） was used to study the electrophysical property of rock saturated with fluid because of its advantages over conventional numerical approaches in handling complex pore geometry and boundary conditions. The digital core model was constructed through the accumulation of matrix grains based on their radius distribution obtained by the measurements of core samples. The flow of electrical current through the core model saturated with oil and water was simulated on the mesoscopic scale to reveal the non-Archie relationship between resistivity index and water saturation （I-Sw）. The results from LBM simulation and laboratory measurements demonstrated that the I-Sw relation in the range of low water saturation was generally not a straight line in the log-log coordinates as described by the Archie equation. We thus developed a new equation based on numerical simulation and physical experiments. This new equation was used to fit the data from laboratory core measurements and previously published data. Determination of fluid saturation and reservoir evaluation could be significantly improved by using the new equation.     

优化系数的吸收边界条件方程

在波动方程有限差分的正演及偏移成像过程中，由于引入人为的边界条件，必然会导致产生边界反射，从而影响正演结果及偏移成像的效果。为此，人们提出了各种吸收边界条件方程，但其吸收效果均不理想。本文在最小二乘意义下，用共轭梯度导出了优化求解吸收边界条件方程中的系数，使得吸收边界条件方程尽可能与内部方程相逼近，从而使得边界反射波被全部或部分吸收掉。理论分析和偏移成像实践表明，优化系数的吸收边界条件方程对以各种     

基于REV尺度格子Boltzmann方法的页岩气藏渗流规律

页岩储层一般天然微裂缝发育,基质孔隙结构复杂,使得页岩气渗流过程呈现出多尺度多场耦合的特征。为了研究页岩气藏复杂的渗流规律,重构了天然微裂缝发育的页岩储层多孔介质模型,并围绕页岩气多重运移机制对广义格子Boltzmann模型进行了修正,建立了适用于模拟页岩气渗流特征的表征单元体（REV）尺度格子Boltzmann模型（LB模型）,并基于天然微裂缝物性特征以及气体滑脱、吸附/解吸、表面扩散效应等渗流特征对该模型进行了敏感性参数分析。结果表明：当页岩储层天然微裂缝较发育时,微裂缝为气体在基质中流动的主要通道;其中裂缝密度是影响储层表观渗透率的主要参数,裂缝密度增大3～4倍,储层表观渗透率可增大10倍以上,而裂缝长度以及裂缝开度的影响程度均次之;努森数（Kn）是影响页岩气渗流的主要参数,随着Kn增大,克氏效应愈显著,特别当Kn > 0.1时,多孔介质表观渗透率增幅显著增大;页岩储层多孔介质表观渗透率会随着吸附气量的增大而减小,特别是当储层压力较低时,该现象更为显著;气体表面扩散效应对页岩气渗流过程的影响程度大,同等条件下考虑气体表面扩散效应的储层表观渗透率较忽略该效应可提升2～5倍,但提升作用受制于储层吸附气量的多少。该研究成果为页岩气微观渗流理论研究提供了新思路,为页岩气藏高效勘探开发提供了技术支撑。     

Hudson理论中裂隙参数的适用性研究

介绍了Hudson裂缝介质模型；阐述了边界元法的基本原理；建立了纵波在裂缝介质中传播的边界元法数值模型；模拟了裂缝密度、裂缝纵横比和裂缝填充物对纵波的传播速度和衰减因子的影响；分析了Hudson理论中裂隙参数的有效取值范围以及可能带来的误差。模拟结果表明，边界元法可以很好地模拟含裂缝非均匀介质中波的传播问题，裂缝密度对纵波传播速度和衰减因子有一定的影响。在地震长波长情况下，相同特征尺度的圆形干裂缝比薄形干裂缝对弹性波的散射可以进行更有效的衰减；裂缝饱和充水后各向异性程度下降，且纵横比值越小，岩石介质性质对裂缝流体越敏感。     

组合边界条件下二维三分量TTI介质波场数值模拟

从TTI介质一阶应力—速度方程出发,利用旋转交错网格高阶有限差分方法,将非分裂完全匹配层(Non-spliting Perfect Match Layer,简称NPML)边界吸收条件和自由边界条件相结合形成组合边界条件,进行了二维三分量TTI介质弹性波场数值模拟。波场快照和炮记录表明:①采用非分裂式边界条件能较好地消除近地表大角度入射波和瞬逝波;②组合边界条件与NPML边界吸收条件相比,不仅有效地压制了边界反射,同时实现了对自由地表的模拟,获得了丰富的全波场信息,其中在地表产生的PS转换横波作为一种特殊的横波现象,可为近地表结构调查以及多波波场分析等提供有益信息;③自由地表引起的面波以及多次波对偏移结果有着重要影响,因此在实际地震资料处理中应当充分考虑自由地表条件对波场的影响效应。数值模拟结果证实了组合边界条件下二维三分量TTI介质波场数值模拟方法的可行性和正确性。     

基于流线EnKF油藏自动历史拟合

作为一种油藏参数动态估计的优化算法,将集合卡尔曼滤波EnKF与流线方法相结合,在算法实现过程中,选择具有更高计算精度和更快计算速度的流线模拟器进行预测,将油藏模型自动还原为一系列沿流线的一维模型,使数值弥散和受网格划分的影响达到最小,保持了明显的驱替前缘。采用序贯高斯协模拟方法生成一组地质模型,通过实时观测数据（产量及含水率等）,连续动态更新油藏模型的静态参数（渗透率及孔隙度等）、动态参数（压力及饱和度等）,同时实现流线分布的更新,直观地反映油藏流体在注采井之间的运动轨迹。采用Bayes理论阐述了流线EnKF数学原理,并通过拟合计算一个二维水驱油藏模型,验证了方法的有效性。     

基于完全匹配层构建新方法的2.5维声波近似方程数值模拟

传统的吸收边界条件只对有限入射角和有限频率范围内的波具有很好的吸收效果,而在有效吸收范围外的边界反射将对数值模拟结果造成污染。为此,在2.5维声波数值模拟中引入完全匹配层吸收边界条件,并提出了一种完全匹配层构建新方法。该方法将完全匹配层构建在研究区域的右边界和下边界的外部,同时改变相应的衰减因子,使之关于完全匹配层的中心呈对称分布。在采用有限差分法计算波场时,将完全匹配层外边界的网格点与相应的研究区域边界的网格点连接起来,使波在传入完全匹配层之后又传入到研究区域中。由于完全匹配层良好的吸收效果,波在完全匹配层中被完全吸收,基本不会再传入到研究区域。在构建相同数量的完全匹配层的情况下,完全匹配层构建新方法使波在传播过程中比以往的构建方法多通过一倍的完全匹配层介质,从而提高了对边界反射的吸收效果。通过均匀速度模型和Marmousi速度模型证实,完全匹配层构建新方法对边界反射的吸收效果较好。     

夏盐南三维静校正方法

夏盐南三维工区地表条件复杂，低、降速带变化大，静校正问题突出，为此，进行了针对性的系统研究。分别用折射初至静校正方法和表层模型静校正方法提供了2套静校正量。由于该区折射层变化频繁，折射波难以追踪，折射初至静校正方法没能得到好的应用效果；表层模型法将微测井和小折射有机地结合起来进行静校正量的计算，加上控制点较多，较好地解决了该区的静校正问题。