【技术实现步骤摘要】
一种叠前地震直接反演的方法
本专利技术涉及叠后地震反演
，特别涉及一种叠前地震直接反演的方法。
技术介绍
与叠后地震反演不同，叠前地震反演利用未叠加或部分叠加的地震数据，除纵波阻抗信息外，可以反演出横波阻抗、密度、岩石模量、流体敏感参数、物性参数、各向异性参数、吸收参数等信息。目前，叠前地震反演主要分为基于波动方程的叠前反演和基于地震波反射系数及其近似公式的叠前反演。其中，基于波动方程的叠前反演受计算效率和稳定性的影响，离实际应用还有一定的距离。目前，叠前地震反演主要以Zoeppritz公式及其近似式为理论基础，可分为AVA/AVO直接反演和弹性阻抗反演两类。其中，AVA/AVO直接反演是指利用CMP道集、CRP道集或者角度部分叠加道集直接进行叠前参数反演，该反演方法已经得到了广泛的应用。但AVA/AVO直接反演方法对噪声敏感，为了提高反演结果的稳定性，许多学者开始对使用角度部分叠加的弹性阻抗反演方法进行研究。弹性阻抗的概念由Connolly(1999)提出，之后由Whitcombe(2002)对其扩展。利用这...
【技术保护点】
1.一种叠前地震直接反演的方法，其特征在于，包括：/n在考虑多个界面和N个角度的情况下，通过对纵波阻抗、横波阻抗及密度反射系数分别进行处理，得到纵波阻抗反射系数序列、横波阻抗反射系数序列及密度反射系数序列；/n利用所述纵波阻抗反射系数序列、横波阻抗反射系数序列及密度反射系数序列，分别建立纵波阻抗低频模型、横波阻抗低频模型以及密度低频模型；/n根据拉格朗日乘子法和所述纵波阻抗低频模型、横波阻抗低频模型以及密度低频模型，构建纵波阻抗、横波阻抗和密度同步约束的反演目标函数：/n利用凸函数差异算法DCA和交替方法乘子法ADMM对所构建的纵波阻抗、横波阻抗和密度同步约束的反演目标函数...
【技术特征摘要】
1.一种叠前地震直接反演的方法，其特征在于，包括：
在考虑多个界面和N个角度的情况下，通过对纵波阻抗、横波阻抗及密度反射系数分别进行处理，得到纵波阻抗反射系数序列、横波阻抗反射系数序列及密度反射系数序列；
利用所述纵波阻抗反射系数序列、横波阻抗反射系数序列及密度反射系数序列，分别建立纵波阻抗低频模型、横波阻抗低频模型以及密度低频模型；
根据拉格朗日乘子法和所述纵波阻抗低频模型、横波阻抗低频模型以及密度低频模型，构建纵波阻抗、横波阻抗和密度同步约束的反演目标函数：
利用凸函数差异算法DCA和交替方法乘子法ADMM对所构建的纵波阻抗、横波阻抗和密度同步约束的反演目标函数进行求解，得到纵波阻抗、横波阻抗和密度；
其中，所述N为正整数。
2.根据权利要求1所述的方法，其特征在于，所述纵波阻抗反射系数序列、横波阻抗反射系数序列及密度反射系数序列包括：
rp＝D′mp
rs＝D′ms
rρ＝D′mρ
式中，所述mp、ms、mρ是纵波阻抗、横波阻抗、密度的自然对数序列，所述D′为一阶差分矩阵，
3.根据权利要求2所述的方法，其特征在于，所述利用所述纵波阻抗反射系数序列、横波阻抗反射系数序列及密度反射系数序列，分别建立纵波阻抗低频模型、横波阻抗低频模型以及密度低频模型包括：
根据所述纵波阻抗反射系数序列和纵波阻抗与反射系数的关系式，构建纵波阻抗的自然对数序列与纵波阻抗序列的关系式；
根据低通滤波矩阵和所述纵波阻抗的自然对数序列与纵波阻抗序列的关系式，建立纵波阻抗低频模型。
4.根据权利要求2所述的方法，其特征在于，所述利用所述纵波阻抗反射系数序列、横波阻抗反射系数序列及密度反射系数序列，分别建立纵波阻抗低频模型、横波阻抗低频模型以及密度低频模型包括：
根据所述横波阻抗反射系数序列和横波阻抗与反射系数的关系式，构建横波阻抗的自然对数序列与横波阻抗序列的关系式；
根据低通滤波矩阵和所述横波阻抗的自然对数序列与横波阻抗序列的关系式，建立横波阻抗低频模型。
5.根据权利要求2所述的方法，其特征在于，所述利用所述纵波阻抗反射系数序列、横波阻抗反射系数序列及密度反射系数序列，分别建立纵波阻抗低频模型、横波阻抗低频模型以及密度低频模型包括：
根据所述密度反射系数序列和密度与反射系数的关系式，构建密度的自然对数序列与密度序列的关系式；
根据低通滤波矩阵和所述...
【专利技术属性】
技术研发人员：聂文亮，相菲，曹湘，冯丽源，文晓涛，聂祥飞，
申请(专利权)人：重庆三峡学院，
类型：发明
国别省市：重庆;50
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。