井间地震层析成像分辨率研究.docVIP

　第26卷第3期　2002年6月 物　探　与　化　探 GEOPHYSICALGEOCHEMICALEXPLORATION Vol.26,No.3　 Jun.,2002　 井间地震层析成像分辨率研究 裴正林1,余钦范2,狄帮让1 (1.石油大学中国石油然气总公司物探重点试验室,北京　100083;2.中国地质大学,北京　100083) 摘要:通过大量数值模型试验,并根据射线密度、射线正交性和点扩散函数,对井间地震层析成像的分辨率问题进行了研究。结果表明,井间地震层析成像的分辨率在空间上是可能性变的,垂直方向的分辨率好于水平方向,而且高速异常体和低速异常体的图像形态和分辨率是非一致的。高速异常体波形CT的图像水平方向最小分辨率距离为1/2波长,垂直分辨率最小分辨率距离为3/10波长;而低速异常体的波形CT图像的水平方向的最小分辨率约为2/5波长,垂向最小分辨率约1/5个波长;走时CT的图像水平方向最小分辨率为3个像元,垂直方向的最小分辨率距离为2个像元;当井深与井间距之比为2时,图像分辨率最高。关键词:井间地震;层析成像;分辨率 中图分类号:P631.4;P631.8　　　文献标识码:A　　　:-)0218-07 　　后,构和岩性分布。因此,对于解的可靠性和分辨率的研究就显得非常重要。在线性反演中,Backus和Gilbert(1969)提出用分辨率矩阵评价解的分辨率。冯锐(1991)提出了射线密度矢量和椭圆的方法。Humphreys和Clayton(1988)提出了检测板方法。而RectorⅢ和Washbourne(1994)、赵永贵等(1996)则采用射线密度和射线正交性评价解的分辨率。 井间地震层析成像的分辨率大小,主要受反演算法、观测系统、震源频率、目标体的形态及性质(高速体还是低速体)等方面的影响。近年来,许多研究者对井间地震走时层析成像分辨率问题进行了研究,并取得了可喜的成果。Williamson(1991,1993)从散射成像与数值模拟的角度研究,认为该方法所能分辨的最小尺度不小于第一Fresnel带;Schuster(1996)从均匀介质中的孤立散射体的散射出发,导出了井间地震走时层析成像纵横分辨率估计式:Δzx0,Δx=(4x0/H)λx0/4,x0=D/2,D是井间距,λ是波长,H是井深。曹俊兴等(1995)根= ,。笔者利用数值模拟试验研究了井间地震弯曲射线走时成 像、波形成像的分辨能力。 1　解估计的分辨率矩阵与图像分辨率评价 的点扩散函数 1.1　分辨率矩阵 井间地震层析成像的线性化反演的关键是通过研究广义逆的性质对解估计进行预测和评价,使重建的图像符合地质意义。线性化反演可以归结为解线性方程组 (1)Ax=d, 对于射线CT,式中,d为地震走时的观测值矢量,x为真实速度模型矢量。Am×n是Jacobi矩阵,是根据速度模型计算出的地震波射线参数矩阵。通常,上述方程没有严格的真解,只有近似解。设A的广义逆为A-g,方程(1)的广义解估计表示为 -g (2 )x=Ad 把(1)式代入(2)式,则有 x=A -g Ax=Rx,(3) 据费马原理,对均匀背景下速度异常体的井间地震层析成像分辨率给出了估算式。杨文采(1992)从计算像元方面来考虑,认为该方法技术在理想情况下,最佳分辨宽度约3个像元。陈国金(1996)根据直射 收稿日期:2001-05-11 式中,R=A-gA为方阵。当R=I(I为单位方阵)时,有^x=x,说明所有模型参数的估计值都是精确解。但是在实际问题中,常常是R≠I,这时,模型参数的某个估计值^xk可以表示为 J ^xk= j=1 ∑R kjxj,(4) 　3期裴正林等: 井间地震层析成像分辨率研究 ?219? 式中,J表示速度模型矢量维数。(4)式说明了当R≠I时,模型参数中的任意一个估计值都与该点附近的若干参数值有关,它实际上是真实模型{xj},j=1,…,J的空间加权平均,其中的权为矩阵R的行向量。R越接近单位阵,则^x越接近真解x,解估计的某个参数与其它参数之间的关系不密切,解估计的分辨率就越高,因此,R被称为模型的分辨矩阵(Backus和Gilbert,1968)。对^xk的分辨率的定量估计可以用分辨核表示 J 一般来说,当源位置变化时,h(x,y)也会有所变化。如果一个点源在物场中移动时,它在像场上的像只移动位置而不改变函数形态,这种成像系统称为空不变的。对于井间地震CT而言,成像系统是空可变的。 对于上述映射系统,可利用系统的点扩散函数 h来描述和预估成像图像的质量。 rk= j=1 ∑(R kj -δj) (5) 在成像理论中,空间分辨率一般定义为:若物函 数为2个相隔一定距离的点,所得像函数一般有2个峰,当2点逐渐靠近时,2峰间的下凹渐渐变小,终于合成为一个峰