双相介质中储层波场特征数值模拟研究  

RESERVOIR WAVE FIELD CHARACTERISTIC NUMERICAL SIMULATION IN TWO-PHASE MEDIA

作　　者：李勇[1,2] 肖富森[3] 马廷虎[3] 李义明[3] 谭勇[3] 

机构地区：[1]成都理工大学"地球探测与信息技术"教育部重点实验室,四川成都610059 [2]成都理工大学信息工程学院,四川成都610059 [3]中国石油西南油气田分公司川东开发公司,重庆400021

出　　处：《西南石油大学学报（自然科学版）》2009年第3期13-15,176,共3页Journal of Southwest Petroleum University(Science & Technology Edition)

基　　金：国家自然科学基金项目(40874053/D0409);成都理工大学地球探测与信息技术教育部重点实验室开放基金(2008DTKF002)

摘　　要：根据双相介质弹性波方程,并对弹性波方程进行分离得到了膨胀波方程。利用高阶差分对膨胀波方程进行数值模拟,边界处理采用完全匹配层的吸收边界条件,算法实现了任意高阶差分和自动加载边界条件的纵波方程正演模拟。计算结果表明,在双相各向同性介质中存在快纵波、慢纵波,两种纵波有明显的区别,快纵波的速度远大于慢纵波的速度。在分界面上,快纵波要产生透射快纵波、透射转换慢纵波,反射快纵波、反射转换慢纵波。慢纵波具有很强的衰减性,耗散系数越小,慢纵波越明显,反之亦然。该算法的高阶差分形式可以显著地降低数值频散,有效提高地震波正演计算的精度,适应性好,操作简便灵活。Based on the equation of elastic wave propagation in two-phase media and by separating the equation of elastic wave,the dilatational wave equation in two-phase in media is derived.The higher finite-difference (FD) and the perfectly matched layer boundary conditions are used to numerical simulation.The numerical simulation of the any higher-order difference and automatic loading boundary conditions is realized.The results of the calculation show that in the two-phase media,there are the fast dilatational wave and slow dilatational wave and they are different obviously.The velocity of the fast P wave is higher.On the interface,incident fast P-wave will generate two kinds of transmission and two kinds of reflection waves.The slow P-wave has intensive attenuation.The dissipation modulus is smaller,then the slow P-wave is more obvious,vice versa.The high-order FD pattern of this method can reduce the numerical dispersion markedly,and improve the accuracy of numerical simulation effectively and with good adaptability.The operation is convenient and flexible.

关 键 词：双相介质 波动方程 波场特征 有限差分 数值模拟 

分 类 号：TE122.2[石油与天然气工程—油气勘探]