【技术实现步骤摘要】
基于γ_Zoeppritz方程的叠前反演方法及系统
本专利技术涉及油气勘探领域，更具体地，涉及一种基于γ_Zoeppritz方程的叠前反演方法及系统。
技术介绍
叠前地震数据中包含了丰富的偏移距信息，振幅随偏移距/入射角的变化揭示了地下介质的岩性变化和孔隙内流体成分变化。因此利用叠前AVA同步反演可以从角度部分叠加地震数据中提取多种岩石弹性参数，其中纵横波速度比对储层岩性及孔隙内流体的变化更为敏感，是地震解释人员偏爱的“烃类指示因子”。然而常规叠前弹性参数反演需要首先反演纵波速度\纵波阻抗、横波速度\横波阻抗，然后再转化为纵横波速度比，弹性参数的间接转化会造成误差累积，从而影响纵横波速度比的反演精度。另外，常规叠前弹性参数反演大都基于Zoeppritz方程的近似公式，需要满足中小角度假设以及垂向上弹性参数变化缓慢等条件，因此进一步限制了其应用范围与反演精度。前人对叠前地震反演问题进行了深入的研究。叠前AVA反演最早可追溯至Smith提出的加权叠加方法，该方法属于带限反演，未考虑到子波的带限效应，并且反演结果依...
【技术保护点】
1.一种基于γ_Zoeppritz方程的叠前反演方法，其特征在于，包括：/n步骤1：构建低频软约束项；/n步骤2：构建地震数据拟合差项；/n步骤3：构建反射系数稀疏约束项；/n步骤4：根据所述地震数据拟合差项、所述低频软约束项以及所述反射系数稀疏约束项，构建反演目标函数；/n步骤5：根据所述反演目标函数获取弹性参数反演结果。/n
【技术特征摘要】
1.一种基于γ_Zoeppritz方程的叠前反演方法，其特征在于，包括：
步骤1：构建低频软约束项；
步骤2：构建地震数据拟合差项；
步骤3：构建反射系数稀疏约束项；
步骤4：根据所述地震数据拟合差项、所述低频软约束项以及所述反射系数稀疏约束项，构建反演目标函数；
步骤5：根据所述反演目标函数获取弹性参数反演结果。
2.根据权利要求1所述的基于γ_Zoeppritz方程的叠前反演方法，其中，所述步骤1具体包括：
步骤11：获取纵波速度模型、纵横波速度比模型以及密度模型；
步骤12：分析地震数据频谱，确定低通频与低截频，进而计算三弹性参数自然对数的低频趋势；
步骤13：构建低通滤波算子，并与所述三弹性参数自然对数的低频趋势，构建低频软约束项。
3.根据权利要求2所述的基于γ_Zoeppritz方程的叠前反演方法，其中，所述步骤2具体包括：
步骤21：获取角度叠加数据；
步骤22：读取角度子波，构建角度子波褶积矩阵；
步骤23：根据所述三弹性参数自然对数的低频趋势，构建反演初始模型；
步骤24：根据精确Zoeppritz方程，构建γ_Zoeppritz方程；
步骤25：根据所述γ_Zoeppritz方程、所述反演初始模型以及所述角度子波褶积矩阵，并结合泰勒一阶展开公式，构建基于γ_Zoeppritz方程的AVO正演方程组；
步骤26：根据所述基于γ_Zoeppritz方程的AVO正演方程组与所述角度叠加数据，构建地震数据拟合差项。
4.根据权利要求3所述的基于γ_Zoeppritz方程的叠前反演方法，其中，所述步骤24中的所述γ_Zoeppritz方程为：
其中，分别表示上层介质的纵波速度、纵横波速度比以及密度的自然对数，分别表示下层介质的纵波速度、纵横波速度比以及密度的自然对数，α1、α2、β1、β2分别表示纵波入射角、纵波透射角、转换波反射角以及转换波透射角，Rpp、Rps、Tpp、Tps分别表示纵波反射系数、转换波反射系数、纵波透射系数以及转换波透射系数。
5.根据权利要求3所述的基于γ_Zoeppritz方程的叠前反演方法，其中，步骤25中，所述基于γ_Zoeppritz方程的AVO正演方程组为：
H＝Gm(2)
其中，H＝Δd+Gm0，Δd表示实测地震数据与合成地震数据的差，m0表示反演初始模型；m表示反演参数；
W(θk)表示入射角为θk的角度子波褶积矩阵，Lvp(θk),Lγ(θk),Lρ(θk)分别表示由入射角为θk的纵波反射系数对纵波速度自然对数、纵横波速度比自然对数以及密度自然对数的一阶偏导数组成的双对角矩阵，其中k＝1,2,…M，
其中，表示第...
【专利技术属性】
技术研发人员：张丰麒，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油勘探开发研究院，
类型：发明
国别省市：北京;11
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。