[发明专利]一种地震反演的方法和装置

权利要求书

1.一种地震反演的方法，其特征在于，包括：

获取工区的转换波数据；

基于所述转换波数据得到所述工区中分角度叠加的地震数据；

根据所述工区中参考井的参考测井数据、所述地震数据中的入射角度信息，计算得到所述参考井的转换波弹性阻抗曲线；

基于所述工区中的层位数据以及所述转换波弹性阻抗曲线，得到所述工区的初始弹性阻抗模型；

利用所述地震数据、所述初始弹性阻抗模型、所述地震数据对应的地震子波反演得到所述工区的转换波弹性阻抗；

基于所述转换波弹性阻抗中的井旁道转换波弹性阻抗、所述参考测井数据，反演得到表征所述井旁道转换波弹性阻抗与所述参考测井数据之间关系的关系系数；

根据所述关系系数，结合所述转换波弹性阻抗，计算得到所述工区的地震属性数据；

其中，基于所述转换波数据得到所述工区中分角度叠加的地震数据，包括：

基于所述转换波数据得到多组不同入射角范围的地震数据；

相应的，根据所述工区中参考井的参考测井数据、所述地震数据中的入射角度信息，计算得到所述参考井的转换波弹性阻抗曲线，包括：

根据所述工区中参考井的参考测井数据、预设的入射角度，计算得到所述参考井中所述预设的入射角度的转换波弹性阻抗曲线，其中，根据所述地震数据的入射角范围确定所述预设的入射角度，其中，所述参考测井数据包括：所述参考井的横波速度和密度；

其中，按照以下公式根据所述工区中参考井的参考测井数据、所述地震数据中的入射角度信息，计算得到所述参考井的转换波弹性阻抗曲线：

其中，SEIO(θ)＝ρ₀^a(θ)-1β₀^b(θ)-1；

a(θ)＝sinθ[-(1+2γ)+γ(3γ/2+1)sin²θ+γ⁴sin⁴θ]；

b(θ)＝2γsinθ[-2+(1+2γ)sin²θ+γ³sin⁴θ]；

其中，SEI(i,θ)表示所述参考井的转换波弹性阻抗曲线，θ表示所述预设的入射角度，ρ(i)表示所述参考测井数据中第i个采样点的密度，β(i)表示所述参考测井数据中第i个采样点的横波速度，ρ₀表示所述参考测井数据中密度的平均值，β₀表示所述参考测井数据中横波速度的平均值，γ表示地下介质的纵横波速度比，i＝1,2,…,m。

2.如权利要求1所述的方法，其特征在于，按照以下公式基于所述转换波弹性阻抗中的井旁道转换波弹性阻抗、所述参考测井数据，反演得到表征所述井旁道转换波弹性阻抗与所述参考测井数据之间关系的关系系数：

其中，P(i,θ_j)＝ln(SEI(i,θ_j)·SEIO(θ_j))＝a(θ_j)lnρ(i)+b(θ_j)lnβ(i)；

a(θ_j)＝sinθ_j[-(1+2γ)+γ(3γ/2+1)sin²θ_j+γ⁴sin⁴θ_j]；

b(θ_j)＝2γsinθ_j[-2+(1+2γ)sin²θ_j+γ³sin⁴θ_j]；

其中，θ_j表示第j个所述预设的入射角度，ρ(i)表示所述参考测井数据中第i个采样点的密度，β(i)表示所述参考测井数据中第i个采样点的横波速度，γ表示地下介质的纵横波速度比，c_j表示第j个关系系数，i＝1,2,…,m，j＝1,2,…,n。