声波全波计算方法的研究及应用.pdfVIP

徐佳等：声波全渡计算方法的研究及应嗣 声波全波计算方法的研究及应用 TheResearchandApplicationofAcousticFull——waveCalculating 徐 佳 庞 巨丰 范 斐 董兰屏 王迎辉 (西安石油大学电子工程学院，陕西 西安 710065) 摘 要：随着测井技术的发展，产生了声波全波列测井技术，声波全波测井比声速测井能得到更丰富的地层信息，这些信息对油田的勘探与开发有着 广泛的应用。本文主要介绍了声波全波的计算方法，实轴积分法适用于井外为均匀各向同性弹性介质的情况，而对于复杂介质等井眼物理模型则主 要应用有限差分法 。声波全波信息在确定岩性、岩层孔隙度和含油饱和度等参数的方面有着重要意义。 关键词：声波全波；实轴积分；有限差分 1 声波全波数值计算方法 1．1 实轴积分(离散波数法) = 案 第一类数值方法是计算接收器所接收到的全波波形 和某一种波形(如纵波和斯通利波)。这种方法主要采用 …(走、 ■ 实轴积分法和割线积分法(或极点留数计算)，用时域波形和 频域特性曲线来表示不同地层~．JC-T的全波列曲线。 c(志) 考虑裸眼井的理想模型为无限大弹性介质内的一个 半径为 口的无限长圆孔，即井半径为 a，井筒在上、下两 Q=2 101 个方向上无限延伸；井 内流体和井外地层的参数分别以 下标 1、2区分；令介质密度为p，纵波速度为 Vp，横波速 {11[2k2-1)K0(ffva)+ K丽0(,usa)} 度为 s，流体速度 ；并假设井 内流体和固体都是理想 将A(k)和Q代人 1，得井内任意点处的声压 的弹性介质，建立柱坐标系(rOz)，z轴与井轴重合，声源 r∞ roo 位于坐标系远点。 l=I s(∞)I √ 一 。。 √ 一 OO 考虑到泥浆 中有入射波和反射波，地层中只有折射 纵波和折射横波。且入射波在井轴处幅度大，远离井轴 [ 十 训 幅度逐渐减小；反射波在井壁处幅度最大，向井内传播幅 1．2 有限差分法 度逐渐减小；折射波在井壁处幅度最大，向地层传播幅度 传统的解析方法如离散波数法适用于介质轴对称柱 逐渐减小 。 状分层且声源居 中情况。对含水平、倾斜分层的介质或 用分离变量法求解 ，得 测井声系偏离井轴情况下的测井问题则需半解析或纯数 roo rco 值方法进行研究。有限差分法作为一种经典的数值方 1=l S(CO)I c[Ko(1r)+A(k)jro(1r)] 法，能够对场域的几何形状不规则的问题给出全波列的 (“一如) ∞ (1) 解 ，在声波测井中得到广泛的应用。 roo rco = I s((U)l B(志)厂0( )(“一垃)(z (2) 文献[3]提出了一种声压一速度有限差分的方法，用 r∞ r 速度矢量描述弹性固体介质。这样选择变量给整个差分 ，=l s()I c(k)K1( )(一砬 (