一种具有较小数值耗散的对称TVD格式

为了减小TVD格式的数值粘性,提高TVD格式的分辨度,提出了一种具有较小数值耗散的Harten型TVD格式.数值模拟结果表明,这种Harten型TVD格式比原来的Harten格式^[1]及其改进形式Harten-Yee^[2]的TVD格式能更好地模拟出边界层及涡流.     

Biot／Squirt Model in Viscoelastic Porous Media

Viscoelastic wave equations are derived for both the Biot flow and squirt-flow involved.The relationship between phase velocity or attenuation of viscoelastic wave and macro parameters of the reservoir is clarified in a porous extensive dilatancy anisotropy medium.Numerical models clearly demonstrate that the viscoelastic property of rocks,not the squirt-flow,causes the dispersion and attenuation in the low-frequency range.The attenuation of the quasi SV-wave,SH-wave and the quasi P-wave depend strongly on the directions of permeabilities.Furthermore,the trends of attenuation of the quasi P-wave fast and the quasi SV-wave are inversed in both the high-frequency and low-frequency ranges.Our numerical result is in agreement with the experimental result that was carried out for the Biot/Squirt model by Jorge [Geophys.65(2000)202].     

弹性波方程的紧致差分方法

在对弹性波方程进行数值模拟时 ,低阶差分格式往往产生严重的数值频散 ,高阶显示差分格式需要用较多的网格点 ,不利于边界的处理。而紧致差分格式吸收了它们的优点 ,弥补了它们的不足。为此该文应用紧致差分格式的思想 ,发展了二维情况下弹性波方程初值问题的紧致差分方法 ,研究了它的稳定性 ,并用 Fourier方法分析了显示差分格式和紧致差分格式的相速度误差 ,最后利用紧致差分方法在粗网格条件下对地震波传播进行了数值模拟 ,并同五点四阶中心差分方法的计算结果进行了对比。结果表明 ,求解弹性波方程的紧致差分方法有效 ,且具有比同网格点差分格式更高的计算精度和较小的数值频散。     

Viscoelastic BISQ Model for Low-Permeability Sandstone with Clay

A modified BISQ （Blot/Squirt） model for wave propagation in low-permeability sandstone is developed by introducing the viscoelastic mechanism of a porous skeleton into Dvorkin＇s model. The linear viscoelasticity of the Kelvin Voigt constitutive law is employed to describe the stress-strain relation of a solid frame with clay while the ultrasonic waves propagate through the fluid-saturated sandstone. The phase velocity and attenuation of two p-waves are given based on the present BISQ model. The comparisons between numerical results and experimental data indicate that our viscoelastic model is more realistic and feasible for wave propagation in the low-permeability sandstone, especially with clay, than traditional BISQ models.