大陆架科学钻探CSDP-2井的垂直地震剖面测量

南黄海海相地层的地震波场特征和层位标定一直是困扰地震勘探的重要问题.为了近距离、高精度和高分辨率地观测井周围构造特征和岩石性质引起的波场变化,为地震资料的采集、处理与解释提供地震波衰减规律、速度与层位标定等信息,对大陆架科学钻探CSDP-2井实施了近零偏移距垂直地震剖面（VSP）观测.针对海相地层顶部强反射界面地震波穿透难的问题,采用了大容量气枪震源并设计了气枪阵列组合方式,提高了激发地震波的能量,获得了强反射界面之下清晰的PP、PS下行波和上行波信号.采用了三分量偏振合成、组合滤波和波场分离等处理方法,对VSP观测数据进行处理,获得了海相三叠系—志留系的精细的纵波、横波速度结构和地层吸收因子等物性数据,建立了钻井地层、测井、VSP上行波和多道地震剖面对应关系,实现了不同尺度的地质和地球物理属性资料的有效衔接,标定了钻井地质剖面上各深度地质体的地震反射特性,厘定了过井地震剖面上反射同相轴的地质属性.此次观测取得的纵波、横波速度信息,成为建立南黄海海相地层速度模型主要的资料来源,也是地震资料的岩性反演处理不可缺少的信息.     

油水两相渗流多孔介质弹性波传播动力学模型研究

弹性波在储层渗流场中的传播与衰减规律是研究波场强化采油动力学机理的重要基础.基于等效流体理论和饱和静态流体弹性波传播Biot理论,建立油水两非混相流体渗流条件下储层多孔介质中弹性波传播的动力学模型,通过算例求解与分析,发现含油水两相渗流储层多孔介质中同时存在着3种纵波P1、P2、P3和1种横波S;受频率和含水饱和度的影响,各波波速和品质因子呈现出不同变化规律,4种体波波速与频率、饱和度正相关,P1、P2波品质因子与饱和度正相关,P3和S波品质因子与饱和度负相关;最后,通过与传统静态弹性波模型结果对比,进一步分析了宏观渗流场对弹性波传播特征的影响规律,为揭示低频人工地震波辅助强化采油技术的动力学机理和工艺参数优化提供了重要理论依据.     

浅谈边坡稳定性分析方法

综述边坡稳定性分析方法及其研究的一些新进展,从定性、定量、非确定性、流变以及动力分析等方面进行了详细论述,指出各种方法的优缺点。     

非均匀弹性介质地震波动模拟的Fourier变换方法

一、前言 地震模型技术是应用地震法研究地下储油地层分布的基础工作,它可以用来研究地,震波的传播以及各种波的特性,可以通过修改地质模型,反复进行正演,实现对实际地震剖面的最优逼近,以实现地震反问题。本文在文献[1]的基础上,应用Fourier变换方法研究了用非均匀弹性介质波动方程正演问题的解来模拟地震波的传播,探索了该方法对非均匀弹性介质地震模型的有效性及实用性。     

Poynting and polarization vectors based wavefield decomposition and their application on elastic reverse time migration in 2D transversely isotropic media

With the progress in computational power and seismic acquisition, elastic reverse time migration is becoming increasingly feasible and helpful in characterizing the physical properties of subsurface structures. To achieve high-resolution seismic imaging using elastic reverse time migration, it is necessary to separate the compressional (P-wave) and shear (S-wave) waves for both isotropic and anisotropic media. In elastic isotropic media, the conventional method for wave-mode separation is to use the divergence and curl operators. However, in anisotropic media, the polarization direction of P waves is not exactly parallel to the direction of wave propagation. Also, the polarization direction of S-waves is not totally perpendicular to the direction of wave propagation. For this reason, the conventional divergence and curl operators show poor performance in anisotropic media. Moreover, conventional methods only perform well in the space domain of regular grids, and they are not suitable for elastic numerical simulation algorithms based on non-regular grids. Besides, these methods distort the original wavefield by taking spatial derivatives. In this case, a new anisotropic wave-mode separation scheme is developed using Poynting vectors. This scheme can be performed in the angle domain by constructing the relationship between group and polarization angles of different wave modes. Also, it is performed pointwise, independent of adjacent space points, suitable for parallel computing. Moreover, there is no need to correct the changes in phase and amplitude caused by the derivative operators. By using this scheme, the anisotropic elastic reverse time migration is more efficiently performed on the unstructured mesh. The effectiveness of our scheme is verified by several numerical examples.     

A triple porosity scheme for fluid/solid substitution: theory and experiment

Quantifying the effects of pore-filling materials on elastic properties of porous rocks is of considerable interest in geophysical practice. For rocks saturated with fluids, the Gassmann equation is proved effective in estimating the exact change in seismic velocity or rock moduli upon the changes in properties of pore infill. For solid substance or viscoelastic materials, however, the Gassmann theory is not applicable as the rigidity of the pore fill (either elastic or viscoelastic) prevents pressure communication in the pore space, which is a key assumption of the Gassmann equation. In this paper, we explored the elastic properties of a sandstone sample saturated with fluid and solid substance under different confining pressures. This sandstone sample is saturated with octadecane, which is a hydrocarbon with a melting point of 28°C, making it convenient to use in the lab in both solid and fluid forms. Ultrasonically measured velocities of the dry rock exhibit strong pressure dependency, which is largely reduced for the filling of solid octadecane. Predictions by the Gassmann theory for the elastic moduli of the sandstone saturated with liquid octadecane are consistent with ultrasonic measurements, but underestimate the elastic moduli of the sandstone saturated with solid octadecane. Our analysis shows that the difference between the elastic moduli of the dry and solid-octadecane-saturated sandstone is controlled by the squirt flow between stiff, compliant, and the so-called intermediate pores (with an aspect ratio larger than that of compliant pore but much less than that of stiff pores). Therefore, we developed a triple porosity model to quantify the combined squirt flow effects of compliant and intermediate pores saturated with solid or viscoelastic infill. Full saturation of remaining stiff pores with solid or viscoelastic materials is then considered by the lower embedded bound theory. The proposed model gave a reasonable fit to the ultrasonic measurements of the elastic moduli of the sandstone saturated with liquid or solid octadecane. Comparison of the predictions by the new model to other solid substitution schemes implied that accounting for the combined effects of compliant and intermediate pores is necessary to explain the solid squirt effects.     

弹性波逆散射微扰论

本文将普遍声逆散射微扰论应用于弹性波层析成像问题,在Born变换下推出了以旋转角为补偿参数的各阶微扰重建公式,实现了对非均匀各向同性散射体内3个参数（质量密度ρ和两个Lamé系数λ,μ）的同时重建. 对于层析成像问题,在弹性波的传播过程中P波与SV波有耦合,但它们不会和SH波发生耦合,于是可以得到3个形式相对简单的标量方程. 在Born变换下,在散射波中引入微扰参数,将散射体的3个参数分别按该微扰参数展开,然后利用二维自由空间的Green函数分别得到散射的P波、SV波和SH波的积分表示. 最后,经一维傅氏变换后,得到Born变换下散射体3个参数的各阶微扰重建公式.     

常熟台倾斜仪降雨干扰定量分析

通过常熟地震台形变山洞内三套倾斜仪观测资料的对比分析,找到降雨、蓄水干扰对观测的影响,并利用弹性形变理论验证其可靠性,分析表明,每次强降雨均伴随有来自西南方向的倾斜干扰,其原因为雨水向湖中汇聚,水位上升产生弹性形变,从而引起地表倾斜.     

移动状态边界面法在黏性土不排水蠕变孔压分析中的应用

在黏性土的各种流变现象的数学模拟研究中,不排水蠕变孔压方程几乎是空白。基于临界状态土力学理论,采用与时间相关的移动状态边界面法,给出了适用于正常固结黏性土的不依赖状态边界面具体形态的不排水蠕变孔压公式。正常固结黏性土的不排水蠕变孔压可以看作是体积蠕变势和剪应力水平增大造成的剪缩这两种作用的结果,前者与参数 有关（ 为次压缩指数, 为压缩指数）,后者则与状态边界面的具体形态有关。在此基础上进一步给出了可方便用于试验数据分析的双对数和单对数形式的两个公式。以修正剑桥模型椭圆形状态边界面为例,分析了双对数和单对数坐标系下蠕变孔压曲线的斜率m和m′与土的流变参数 的关系。采用提出的方法,很好地预测了Walkers给出的Leda clay黏土的蠕变孔压系数,并基于以上两种作用解释了Leda clay黏土的蠕变孔压形成的机制。