一阶速度-应力Biot双相各向同性介质弹性波 First-order Velocity-Stress Elastic Wave Field Separation Scheme for Biot Two-phase Isotropic Medium.pdfVIP

第28卷第3期 计 算 物 理 V01．28．No．3 OFCOMPUTATl0NALPHYSICS 2011年5月 CHINESEJOURNAL May，201l 文章编号：100l一246X(2011)03m404舶 一阶速度一应力Biot双相各向同性介质弹性波 波场分离数值模拟 陈可洋 (中国石油大庆油田有限责任公司勘探开发研究院，黑龙江大庆163712) 摘要：提出～种等价的一阶双曲型速度一应力Biot双相各向同性介质弹性波波动方程，以实现双相介质混合 波场中纯快慢纵波和纯横波波场分离的问题．应用散度和旋度理论证明双相介质等价方程波场分离的可行性， 采用高阶交错网格有限差分法构建高精度正演算子，推导其PML吸收边界条件和稳定性条件，并对均匀双相介 质和层状非均匀双相介质模型进行数值模拟试验，准确得到固流相混合弹性波场、被完全分离的纯纵横波波场， 同时边界吸收效果良好，数值模拟精度较高．计算结果还表明，固流相中的快慢纵波相互伴生因而无法实现分 离，且归属于纯纵波波场，流相慢纵波能量比固相慢纵波能量强，这对认识双相介质弹性波的传播规律以及完 善双相介质理论具有重要意义． 关键词：Biot双相各向同性介质；等效波场分离数值模拟方程；纯快慢纵波和纯横波；PML吸收边界条件；固相 和流相 中图分类号：P315．3 文献标识码：A O 引言 目前，地震波场正演数值模拟技术在石油等地球物理勘探和开发、天然地震观测以及超声检测中发挥着 重要的作用，由人工合成的弹性波数值模拟数据有助于实际地震资料的处理和解释，可以深化对弹性波场在 复杂介质中传播的机理的认识，并检验地震资料处理方法和观测系统设计的合理性等…．近些年来对复杂 孔隙介质的研究已成为当今前沿热点内容，其应用领域包括岩土工程、地震灾害预防、石油勘探与开发、解决 如封存cO：的动态监测、提高石油采收率以及评价洋底似反射层下的天然气储量等．通常认为孔隙介质是 由固体颗粒和孔隙流体构成，由此构建的理论包括Biot双相介质理论，其他理论大多在此基础上进行改进 和应用分析，如考虑流体的喷射机制后，形成的BIsQ双相介质理论；考虑多相介质组合情况形成的Santos理 论模型等心。]．到目前为止，已形成了相应孔隙介质理论的一整套高精度波场正演数值模拟方法及其吸收边 界条件，但仍处于进一步完善过程中． 弹性波正演数值模拟技术是认识弹性波传播规律的一种重要手段，其中研究弹性波的波场分离数值模 拟方法逐步兴起，开始作为数值模拟理论和反演问题新的研究方向．马德堂等Mo首次提出满足纵波为无旋 场、横波为无散场的等价方程思路，采用虚谱法求解等价的二阶弹性波动方程，实现了单相介质波场分离数 值模拟，获得较好的数值模拟结果．但是由于伪谱法的人工边界处理困难，计算量较大，使得该方法的推广应 用受到限制．陈可洋等“1提出了基于二阶单相各向同性介质弹性波的高精度正演数值模拟方法和基于散 度和旋度的一阶双曲型单相各向同性介质纵横波波场分离数值方法，取得了高质量的数值模拟结果，且其计 Jianlei等一。1叫提出了基于一阶速度一应力单相各向同性介质 算速度快，波场的信噪比较高．李振春和zhang 弹性波波场分离等价的弹性波动方程，并采用高阶交错网格有限差分法实现波场分离数值模拟，其精度较 高，且不需要对弹性常数(如速度和密度)求导数．目前国内外对双相介质的波场分离数值试验还未见相关 报道．本文提出了一种基于无耗散型Bjot双相各向同性介质的高精度波场分离数值模拟方法，通过对固流 相任一分量的波场正演递推算子进行分裂运算，得到了任一分量各自的纯纵波和纯横波波场，并通过构建高 收稿日期：2010—05一04；修回日期：2010一09一13 作者简介：陈可洋(1983一)，男，浙江诸暨，助理工程师，硕士，主要从事地震渡正演数值模拟和逆时成像方法研究，大庆油田勘探开 发研究院地震处理二室163712，E·mail：key“gchen@163。com． 万方数据