[发明专利]一种地质条件勘探方法

摘要

本发明公开一种地质条件勘探方法。所述方法将笛卡尔坐标系下的地下介质传播速度v(x)和检波器排列的地震记录P(x,t)插值到扁球坐标系中，获取扁球坐标系下的地下介质传播速度模型v(ξ)和检波器排列的地震记录P(ξ,t)；将检波器排列的地震记录P(ξ,t)通过傅里叶变换至频率域的检波器排列的地震记录P(ξ,ω)；在扁球坐标系下，对单个频率的检波器波场沿ξ₃坐标轴方向向下延拓一个Δξ₃，通过速度时移和有限差分校正获取ξ₃+Δξ₃处的检波器波场直至将波场延拓到模型的底部ξ₃′处；将主频频带范围内的不同频率的波场均延拓到模型的底部ξ₃′处，叠加所有频率的偏移量，得到扁球坐标系的叠加后偏移量P(ξ)；将叠后偏移量P(ξ)插值到笛卡尔坐标系，得到最终的偏移剖面P(x)。采用本发明的方法，能够提高地质条件勘探中的成像的精度。

主权项

1.一种地质条件勘探方法，其特征在于，所述方法包括：获取笛卡尔坐标系下地下介质传播速度v(x)、检波器排列的地震记录P(x,t)和震源子波的主频频带范围；确定一个扁球坐标系，所述扁球坐标系包括ξ₁、ξ₂、ξ₃三个坐标轴，将所述笛卡尔坐标系下的地下介质传播速度v(x)和检波器排列的地震记录P(x,t)插值到所述扁球坐标系中，获取扁球坐标系下的地下介质传播速度模型v(ξ)和检波器排列的地震记录P(ξ,t)；将所述扁球坐标系下的所述检波器排列的地震记录P(ξ,t)通过傅里叶变换至频率域的检波器排列的地震记录P(ξ,ω)；在扁球坐标系下，对单个频率的检波器波场沿ξ₃坐标轴方向将波场向下延拓一个Δξ₃，通过速度时移和有限差分校正获取ξ₃+Δξ₃处的检波器波场直至将波场延拓到模型的底部ξ₃′处；将所述震源子波的主频频带范围内的不同频率的波场均延拓到模型的底部ξ₃′处，叠加所有频率的偏移量，得到扁球坐标系的叠加后偏移量P(ξ)；将扁球坐标系的叠加后偏移量P(ξ)插值到笛卡尔坐标系，得到最终的偏移剖面P(x)；所述对单个频率的检波器波场沿ξ₃坐标轴方向通过速度时移，具体包括：所述速度时移通过计算得到，其中ξ为扁球坐标系下一点的空间位置，(ξ₁，ξ₂，ξ₃)为扁球坐标系下ξ点的坐标，ξ_i(i＝1,2,3)为扁球坐标系下ξ点的坐标变量，P(ξ,ω)为ξ点的波场，ω为圆频率，v(ξ)为ξ点的速度，A为控制扁球坐标系宽度的常量，i为虚数单位；所述有限差分校正通过计算得到，其中，ξ_i(i＝1,2,3)为扁球坐标系下ξ点的坐标变量，ξ＝(ξ₁，ξ₂，ξ₃)为扁球坐标系下一点的空间位置，P(ξ,ω)为ξ点的波场，ω为圆频率，v(ξ)为ξ点的速度，A为控制扁球坐标系宽度的常量,a_i和b_i为无量纲的常数，i为虚数单位。