一种砂泥岩薄互层介质多波AVO联合反演方法技术

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

A Multi-wave AVO Joint Inversion Method for Thin Interbedded Sandstone and Mudstone Medium

The invention discloses a joint inversion method of multi-wave AVO for thin interbedded sandstone and mudstone media. For complex thin interbedded media, five parameters of equivalent VTI media are first inverted by using the VTI medium reflection coefficient equation proposed by Ruger, and then the parameters of sandstone and mudstone in thin interbedded media are further inverted by using the relationship between VTI medium parameters and its internal single layer parameters and N/G. Number and interlayer structure N/G show that the inversion results are satisfactory for both VTI medium and thin interlayer internal parameters.

【技术实现步骤摘要】
一种砂泥岩薄互层介质多波AVO联合反演方法
本专利技术涉及地震勘探
，具体涉及一种砂泥岩薄互层介质多波AVO联合反演方法。
技术介绍
随着油气勘探的不断深入，简单的构造型油气藏已经开发殆尽，未来的勘探重点是隐蔽的岩性油气藏。在我国东部和西部的一些盆地，绝大多数的中、新生代含油气盆地均以薄层砂岩泥岩沉积为主[1]。薄互层储层内部砂体通常呈条带状，数量多，厚度远远小于地震的勘探分辨率，并且岩性和厚度横向上变化快，针对薄互层储层的勘探存在很大的问题。由于地震频带内难以区分薄互层，地震波的波长远远大于薄互层，因此国际上常把薄互层整体等效为各向异性VTI(VerticalTransverselyIsotropic)介质进行研究[2-4]。薄互层的互层结构可以用N/G(net-to-gross)来表示，其中N表示砂岩的总厚度，G表示薄互层介质的总厚度[5-6]。Thomsen用P波和SV波的速度以及两个各向异性参数来描述VTI介质的性质[7]，所以薄互层可以通过P波和SV波的速度、两个各向异性参数和N/G来描述。AVO(Amplitudeversusoffset)技术是目前进行油气预测的主本文档来自技高网...

【技术保护点】
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)...

【技术特征摘要】
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；记初始模型参数为M0；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在M...

【专利技术属性】
技术研发人员：芦俊，安莹，王赟，
申请(专利权)人：中国地质大学北京，
类型：发明
国别省市：北京,11