井间地震资料的层位标定.docxVIP

井间地震资料的层位标定 1 合成记录的生成 在时间间隔合成地震资料的生产中，主要有三个问题：（1）反射系数r（n）提取；（2）地震波w（n）3）的提取；（3）反射系数序列r（n）t）和地震波w（n）t对应于井外地震记录。 井间地震给出的一般是深度域的剖面,所以在井间地震资料的地质解释中,用于层位标定的合成地震记录必须是深度域的合成记录。然而,深度域合成记录的制作不同于时间域的制作方法,目前还没有直接适于深度域地震记录的数学模型。据此,要制作深度域的合成记录有两条可能的途径:一是推导出深度域的合成记录公式;另一是利用现有的理论制作时间域的合成记录,然后把它转化到深度域。 在上述前一途径中,由时域褶积模型s(nΔt)=r(nΔt)*w(nΔt)=Μ-1∑m=0r(m)w(n-m)s(nΔt)=r(nΔt)*w(nΔt)=∑m=0M?1r(m)w(n?m),根据时深关系得到s(nΔzv(z))=r(nΔzv(z))*w(nΔzv(z))s(nΔzv(z))=r(nΔzv(z))*w(nΔzv(z)),由于v(z)本身就是深度z的函数,所以无法直接得到类似于时间域褶积公式的深度域公式s(nΔz)=r(nΔz)*w(nΔz),要根据复合函数推导出以nΔz为自变量的深度域公式,还需要知道速度和深度的函数关系v(z)。有人试图把时间域子波转化到深度域,或者从深度偏移的深度域剖面中提取地震子波后,直接利用时间域褶积公式得到合成记录,这种方法无论在理论上还是物理意义上都是讲不通的;还有人仿照时间域褶积公式,考虑到界面上、下速度的不同,在深度域中用子波算子代替时域子波,用深度域不对称子波算子求取合成记录。该方法似乎可行,但是有些理论问题尚待进一步证实。 因此,由时间域合成记录再转化到深度域是现阶段最符合理论根据的一种深度域合成地震记录制作方法。 2 深度合成数据集的准备和生产方法 2.1 多层介质反射系数的频域特征分析 根据声波、密度等测井曲线,求取每个采样点的v(z),v(z)=Δz/Δt(z),并建立时深关系t(z)?t(z)=∑iΔt(zi)。根据速度剖面进行层位划分,具体实现时取一个门槛值Th,如果 |v(zi+1)-v(zi)|Τh 或 |ρi+1v(zi+1)-ρiv(zi)|Τh, 就属于同一个层位,这一层的速度就等于同一层所有v(zi)的平均值。通过分层得到采样点数n和每一层的密度ρi、速度vi和厚度Δzi值。 在根据分层结果计算反射系数序列时,可利用多层介质反射系数的波动方程的平面简谐波解和边界条件的递推关系,求出每一个波阻抗界面上每一频率(可以给出频率范围)的反射系数Ri(kΔf),通过叠加可以得到该分界面的频域反射系数,所有层都如法炮制就得到频率域的反射系数R(kΔf),即反射系数的频谱。采样间隔为Δt,Δf的选择可以通过傅里叶变换的关系ΔtΔf=1/N得到。 2.2 子波的数谱密度 地震子波的求取有多种途径,此研究是通过井间地震资料的井旁道记录采用同态反褶积方法获取。具体过程为:首先,要对深度域的井旁道记录进行深时转换,得到时间域的井旁道记录;然后用复赛谱方法求取子波的频谱W(ω),即时间域的褶积模型为s(t)=r(t)*w(t),其频谱关系为S(ω)=R(ω)W(ω),求其自然对数,得到lnS(ω)=lnR(ω)+lnW(ω),即为地震记录s(t),反射系数r(t)和地震子波w(t)的对数谱,简记为?S(ω)=?R(ω)+?W(ω);对此作反傅里叶变换,得到s(t),r(t),w(t)的复赛谱,记为(s(t)=?r(t)+?w(t),利用复赛谱的性质可以分解出所求的子波复赛谱?w(t),对它作傅里叶变换后再求指数,得到子波的频谱W(ω)。 2.3 井旁记录st的特征 利用S(ω)=R(ω)W(ω)得到合成记录的频谱,作反傅里叶变换就得到时间域合成地震记录s(nΔt)。根据时深关系,求出每个时域采样点对应的深度值,然后按照所要求的深度采样间隔Δz进行重采样,得到深度域合成地震记录。 为使合成记录与井旁地震道有比较好的吻合,可以适当损失测井资料的分辨率,对r(t)作一次带通滤波,通频带一般为10～400 Hz,使测井与井间地震资料的频带相近。 为了监控或比较井旁记录S′(t)与合成记录S(t)的精度或吻合程度,通常采用如下定量指标: (1) 归一化的相关系数Rss′ Rss′=∫Τ0S(t)S′(t)dt√∫Τ0S2(t)dt∫Τ0S′2(t)dt, 利用Rss′的大小判断两记录的相似程度。 (2)均方误差m2 m2=1Τ∫Τ0[S(t)-S′(t)]2dt=2σ2-2Τ∫Τ0S(t)S′(t)dt=2σ2[1-1Τ∫Τ0S(t)S′(t)dtσ2]=2σ2[1-Rss′]. 式中,方差σ2=1Τ∫Τ0S2(t)dt=1Τ