[发明专利]一种砂泥岩薄互层介质多波AVO联合反演方法

摘要

本发明公开了一种砂泥岩薄互层介质多波AVO联合反演方法，针对复杂薄互层介质，利用Ruger提出的VTI介质反射系数方程，首先利用叠前AVA角道集反演出等效VTI介质的5个参数，然后进一步利用VTI介质参数与其内部单层参数和N/G的关系反演出薄互层内部砂岩和泥岩的参数以及互层结构N/G，结果表明，无论是对VTI介质还是薄互层的内部参数，反演均能取得令人满意的效果。

主权项

1.一种砂泥岩薄互层介质多波AVO联合反演方法，其特征在于，包括如下步骤：S1、抽取PP波AVA道集和PS波AVA道集，并将PS波资料压缩到PP波的旅行时上，使PS波与PP波时间剖面相匹配；S2、建立VTI介质参数反演的初始模型：2.1)利用测井资料建立纵横速度和密度的线性初始模型；2.2)根据岩石物理测井信息或者速度分析，得到初始各向异性参数模型；S3、建立薄互层内部参数反演的初始模型：3.1)根据测井资料得到的砂泥岩纵横波速度、密度以及厚度资料；3.2)利用步骤3.1)得到的资料建立砂泥岩纵横波速度、密度和N/G的线性初始模型；S4、多波AVA反演S4.1、VTI介质参数反演：4.1.1)将步骤S2中建立得到的纵横速度和密度的线性初始模型以及初始各向异性参数模型转化为以Δε,Δδ,为模型参数如下：M＝[A1,B1,C1,D1,E1,…Am,Bm,Cm,Dm,Em]T,其中，D＝Δε,E＝Δδ,；m表示反演剖面中每一道的采样点数，则模型参数的维数为5*m；记初始模型参数为M₀；4.1.2)利用式(1)计算模型参数的反射系数序列，并与提取出的子波进行褶积，得到PP波和PS波的合成记录；其中，VP0和VS0分别表示P波和SV波的垂直速度，ρ是密度，ε和δ是VTI介质的各向异性参数，i和j分别为P波的入射角和透射角；下标1和2分别表示反射界面上下层的弹性参数；4.1.3)根据式(2)、(3)计算出模型增量ΔM和目标函数Q1，如果Q1不满足精度要求，则ΔM叠加到初始参数模型上得到新的模型参数，并返回步骤4.1.2)；直到Q1满足精度要求，迭代停止；此时得到一个达到目标分辨率的纵波速度、横波速度、密度和各向异性参数的数据体；ΔM＝(JTJ)‑1JT[Y(M,θ,t)‑Y(M0,θ,t)], (2)Q1＝||(JTJ+λI)·ΔM‑JT[Y(M,θ,t)‑Y(Minv,θ,t)]||2. (3)设实测的地震记录Y(M,θ,t)＝[y(M,θ1,t),y(M,θ2,t),…,y(M,θn,t)]T，其中n表示CMP的道数，θ表示入射角，t表示时间，地震记录等于反射系数序列和子波的褶积，即：Y＝ω(t)*R(θ),式中，ω(t)是地震子波的序列，R(θ)是反射系数序列，具体表示为：对地震记录Y在M0处进行Taylor展开并略去高阶项，得到：ΔY＝JΔM,其中，ΔM＝M‑M0代表模型参数的增量，ΔY＝Y(M,θ,t)‑Y(M0,θ,t)为实测数据地震记录与初始模型参数地震记录计算得到的地震记录之差；J是Jacobian矩阵，是地震记录Y(M,θ,t)对于模型参数M的一阶导数；J表示为：S4.2、薄互层内部参数反演：4.2.1)薄互层内部参数的模型参数定义为：P＝[VPa1,VSa1,ρa1,VPb1,VSb1,ρb1,r1,…VPax,VSax,ρax,VPbx,VSbx,ρbx,rx]；r表示N/G；记薄互层内部参数的初始模型参数为P0；4.2.2)用初始模型参数为P0计算得到VTI介质模型参数S：S＝[VP1,VS1,ρ1,ε1,δ1,…VPx,VSx,ρx,εx,δx]T；式中x表示VTI介质的个数；4.2.2)求解下式，得到模型增量ΔP和目标函数Q2；如果Q2不满足精度要求，则ΔP叠加到初始模型参数上得到新的模型参数，并用新的模型参数计算VTI介质模型参数S，并再次计算下式，直到Q2满足精度要求，迭代停止；此时获得一个达到目标分辨率的砂泥岩纵波速度、横波速度和密度以及砂岩与储层厚度比N/G的数据体；ΔS＝GΔP；Q2＝||(GTG+λI)·ΔP‑GT(S‑Sinv)||2；其中，ΔP＝P‑P0代表模型参数的增量，ΔS＝S‑S0是实际VTI介质参数与通过初始模型参数计算得到的VTI参数之差；G是Jacobian矩阵，是VTI介质模型参数S对于模型参数P的一阶导数,G表示为：