基于横向各向同性BISQ方程的弹性波传播数值模拟-地震学报.PDFVIP

第24卷 第6期 地 震 学 报 Vol.24，No.6 2002 年11 月 (599～606) ACTA SEISMOLOGICA SINICA Nov.，2002 文章编号：0253-3782(2002)06-0599-0008 基于横向各向同性BISQ方程的 弹性波传播数值模拟* 1,2) 1) 1) 杨宽德 杨顶辉 王书强 1) 中国北京100084 清华大学数学科学系 2) 中国云南蒙自市 661100 蒙自师范高等专科学校 摘要 Biot流动和喷射流动是含流体多孔隙介质中流体流动的两种重要力学机制. 近年来，利 用同时处理这两种力学机制的 BISQ(Biot-Squirt)模型，弹性波衰减和频散的问题已被广泛研 究；然而基于BISQ方程的波场数值模拟尚未见到公开的报道．本文从BISQ方程出发，利用交 错网格方法对横向各向同性孔隙介质中不同频率和相界情况，以及双层介质中的弹性波传播进 行数值模拟，研究了在同时考虑两种流动机制作用情况下地震波和声波的传播特性及传播过程 中出现的各种波动现象． 关键词 BISQ方程 波场模拟 横向各向同性 交错网格法 中图分类号：P315.3+ 1 文献标识码：A 引言 近代地震学正在对更逼真的模型进行研究．众所周知，地球介质并不是一简单的弹性 体，油气储层、地表、海岸带、海底层等都是含流体的多孔隙固体介质．已有研究表明， 在这些介质中，流体的流动主要有两种方式(Dvorkin, Nur, 1993)，即基于宏观性质的Biot 流动(Biot, 1956, 1962)和基于微观形式的喷射流动(Mavko，Nur, 1979)．基于Biot流动 的力学机制，Biot（1956，1962）建立了多孔隙介质中的弹性波传播理论，发现了慢速P 波的存在，并得到了实验的验证（Plona,1980）．但在很多情况下这种理论预测的能量耗散 和速度频散要比实际的低，因而对许多实验结果和观测数据不能作出合理的解释．而基于 喷射流动机制的喷射流动理论，由于其理论是通过单个孔隙中流体流动的力学机理而建立， 因而其应用受到了很大的限制．Dvorkin和Nur（1993）针对一维双相系统提出了同时处理 Biot流动和喷射流动两种力学机制的统一的BISQ（Biot-Squirt）模型，并将其理论预测结 果与实验数据进行比较，获得了比Biot理论更准确的结果．之后，Parra（1997, 2000）将 这一理论推广到横向各向同性双相介质情况；杨顶辉① 以及Yang和Zhang(2002)提出了固- 流耦合各向异性的思想，并建立了一般含流体多孔隙各向异性介质中同时包含两种力学机 制的弹性波传播方程（简称BISQ方程），进一步还将其推广到了粘弹性情况(Cheng et al., * 国家自然科学基金和中国石油天然气集团公司物探重点实验室开放基金（GPKL01104）共同资助项目． 2002-01-06 收到初稿，2002-04-08 收到修改稿并决定采用． ① 杨顶辉．1998．孔隙各向异性介质中基于BISQ 模型的弹性波传播理论及有限元方法．石油大学博士后研究报告． 600 地 震 学 报 24卷 2002)．另外，基于微观流场固-流相对运动速度的各向异性，杨顶辉和张中杰(2000)研究 了Biot流动和