基于横向各向同性BISQ方程的弹性波传播数值模拟.pdf

第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)和基于微观形式的喷射流动(MavkoNur, 1979)基于 Biot流 动的力学机制 Biot立了多孔隙介质中的弹性波传播理论 发现了慢速 P波的存在 并得到了实验的验证 Plona,1980但在很多情况下这种理论预测的能量 耗散和速度频散要比实际的低 因而对许多实验结果和观测数据不能作出合理的解释 而 基于喷射流动机制的喷射流动理论 由于其理论是通过单个孔隙中流体流动的力学机理而 建立 因而其应用受到了很大的限制 Dvorkin和 Nur1993针对一维双相系统提出了 同时处理 Biot流动和喷射流动两种力学机制的统一的 BISQ Biot-Squirt模型 并将 其理论预测结果与实验数据进行比较 获得了比Biot理论更准确的结果 之后 Parra1997, 2000将这一理论推广到横向各向同性双相介质情况 杨顶辉 以及Yang和Zhang(2002) 提出了固-流耦合各向异性的思想 并建立了一般含流体多孔隙各向异性介质中同时包含 两种力学机制的弹性波传播方程 简称 BISQ方程 进一步还将其推广到了粘弹性情况 * 国家自然科学基金和中国石油天然气集团公司物探重点实验室开放基金 GPKL01104 共同资助项目 2002-01-06 收到初稿 2002-04-08 收到修改稿并决定采用 杨顶辉 1998 孔隙各向异性介质中基于BISQ 模型的弹性波传播理论及有限元方法 石油大学博士后研究报告 万方数据 600 地 震 学 报 24 卷 (Chenget al., 2002)另外 基于微观流场固-流相对运动速度的各向