Taking into account three important porous media mechanisms during wave propagation （the Biot-flow, squirt-flow, and solid-skeleton viscoelastic mechanisms）, we introduce water saturation into the dynamic governing equations of wave propagation by analyzing the effective medium theory and then providing a viscoelastic Biot/squirt （BISQ） model which can analyze the wave propagation problems in a partially viscous pore fluid saturated porous media. In this model, the effects of pore fluid distribution patterns on the effective bulk modulus at different frequencies are considered. Then we derive the wave dynamic equations in the time-space domain. The phase velocity and the attenuation coefficient equations of the viscoelatic BISQ model in the frequency-wavenumber domain are deduced through a set of plane harmonic solution assumptions. Finally, by means of numerical simulations, we investigate the effects of water saturation, permeability, and frequency on compressional wave velocity and attenuation. Based on tight sandstone and carbonate experimental observed data, the compressional wave velocities of partially saturated reservoir rocks are calculated. The compressional wave velocity in carbonate reservoirs is more sensitive to gas saturation than in sandstone reservoirs.
In heterogeneous natural gas reservoirs, gas is generally present as small patchlike pockets embedded in the water-saturated host matrix. This type of heterogeneity, also called "patchy saturation", causes significant seismic velocity dispersion and attenuation. To establish the relation between seismic response and type of fluids, we designed a rock physics model for carbonates. First, we performed CT scanning and analysis of the fluid distribution in the partially saturated rocks. Then, we predicted the quantitative relation between the wave response at different frequency ranges and the basic lithological properties and pore fluids. A rock physics template was constructed based on thin section analysis of pore structures and seismic inversion. This approach was applied to the limestone gas reservoirs of the right bank block of the Amu Darya River. Based on poststack wave impedance and prestack elastic parameter inversions, the seismic data were used to estimate rock porosity and gas saturation. The model results were in good agreement with the production regime of the wells.
黏声波方程常被用于描述地下介质的黏弹性及波的传播现象,频域有限差分(finite difference frequency domain,FDFD)方法是黏声波和黏弹性波波场模拟的常用工具.目前FDFD黏声波模拟常用的二阶五点方法和优化九点方法在一个波长内的网格点数小于4时误差较大.通过令FDFD系数随一个波长内的网格点数自适应从而提高FDFD方法的精度,本文针对黏声波波场模拟发展了一种适用于不同空间采样间隔之比的通用格式自适应系数FDFD方法.同时,为了验证自适应系数FDFD方法对一般黏声波模型的有效性,本文针对三个典型的黏声波模型,分别采用解析解和基于高阶FDFD的参考解验证了所提出方法的有效性.本方法的FDFD格式通过在传统的二阶FDFD格式的基础上引入相关校正项得到,其中校正项按网格点与中心点的距离进行分类选取,同时校正项对应的自适应FDFD系数不仅和空间采样间隔之比相关,还和一个波长内的采样点数相关.所需的自适应FDFD系数可通过声波方程的数值频散关系和查找表高效给出.数值频散分析表明,在空间采样间隔相等或不等的情况下,以相速度误差不超过1%为标准,通用格式自适应系数FDFD方法所需的一个波长内的采样点数均小于2.5.数值模拟实验表明,对于不同的空间采样间隔之比,相对于常用的二阶五点FDFD方法和优化九点FDFD方法,通用格式自适应系数FDFD方法均可在相似的计算量和内存需求下,有效提高黏声波模拟的精度.
