[发明专利]一种深度域地震子波提取与地震记录合成的一体化方法

摘要

本发明提供一种利用深度域地震资料提取深度域地震子波，并将它用于直接合成深度域地震记录的技术。该方法既考虑了地震子波在空间中的变化特点，又可以得到包含整个地层范围内有效信息的深度域地震记录。本发明包括以下主要步骤：(1)由某口井的测井数据获得深度、速度和密度信息，由此计算出对应的反射系数r；(2)为满足线性时不变条件提取子波，以常速度v_c为参考速度，将该井的井旁地震道S以及反射系数r进行常速度深度变换，得到变换后的道S₁以及反射系数r₁，并基于Gibbs抽样方法提取深度域地震子波w；(3)直接利用测井速度v、反射系数r以及提取的深度域子波w作为输入制作深度域合成地震记录。

主权项

1.一种深度域地震子波提取与地震记录合成的一体化方法，其包括以下主要步骤：1)由某口井的测井数据获得深度、速度和密度信息，由此计算出对应的反射系数r；2)为满足线性时不变条件提取子波，以常速度v_c为参考速度，将该井的井旁地震道S以及反射系数r进行常速度深度变换，得到变换后的道S₁以及反射系数r₁，并基于Gibbs抽样方法提取深度域地震子波w；3)直接利用测井速度v、反射系数r以及提取的深度域子波w作为输入制作深度域合成地震记录；其特征在于步骤3)采用了以下深度域合成地震记录制作方法：①计算的反射系数和速度场作为深度域地震记录制作的输入，将反射系数r作为初始地震记录p(z)，并开拓p(z)为p(z；j)；对深度域子波w进行傅立叶变换得到波数域子波W；②取频率j，当深度为n×Δd时，计算式中，式中，v(n×Δd)为深度为n×Δd时的速度，T为时间采样点数，Δd为深度采样间隔；将与p(n×Δd；j)相乘，得到p((n+1)×Δd；j)：将此时的地震记录p((n+1)×Δd；j)作为向深度(n+2)×Δd递推时的输入；③返回步骤②，继续深度循环，直到达到最大深度为止；对应此时的频率j得到：p(j)＝p(D×Δd；j)×W(j)，式中D为最大深度采样点数；④另取频率j+1，重复步骤②和③，直到所有频率循环完毕为止，得到地震记录p[p(1),p(2),…,p(j‑1),p(j)]；⑤取频率j，当深度为n×Δd时，计算并与p(n×Δd；j)相乘，得到深度为(n+1)×Δd的地震记录p((n+1)×Δd；j)：⑥返回步骤⑤，继续深度循环，直到达到最大深度为止，保留此时的地震记录p[p(Δd；j),p(2×Δd；j),…,p((D‑1)×Δd；j),p(D×Δd；j)]；⑦另取频率j+1，重复步骤⑤和⑥，直到所有频率循环完毕为止，最终得到深度域合成地震记录P(z)。