[发明专利]基于混合Markov邻域的叠前地震多参数同步反演非线性方法在审
申请号: | 202010765384.1 | 申请日: | 2020-07-31 |
公开(公告)号: | CN112255670A | 公开(公告)日: | 2021-01-22 |
发明(设计)人: | 严立志;张宏兵;尚作萍;韩飞龙 | 申请(专利权)人: | 河海大学 |
主分类号: | G01V1/28 | 分类号: | G01V1/28;G01V1/30 |
代理公司: | 南京纵横知识产权代理有限公司 32224 | 代理人: | 董建林 |
地址: | 210024 江*** | 国省代码: | 江苏;32 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 基于 混合 markov 邻域 地震 参数 同步 反演 非线性 方法 | ||
1.一种基于混合Markov邻域的叠前地震多参数同步反演非线性方法,其特征在于,
获取测井数据,包括测井速度和密度;
根据测井数据和Zoeppritz方程获取反射系数;
根据测井数据、假设条件和背景趋势关系获取拟合偏差值;
根据反射系数、拟合偏差值对目标函数进行求解得到反演数据。
2.根据权利要求1所述的一种基于混合Markov邻域的叠前地震多参数同步反演非线性方法,其特征在于,所述Zoeppritz方程为:
其中,RPP代表纵波入射纵波反射的反射系数,RPS代表纵波入射横波反射的反射系数,TPP代表纵波入射纵波反射的透射系数,TPS代表纵波入射横波反射的透射系数,v是测井速度,ρ1为上覆介质的密度,ρ2为下伏介质的密度,vs1为上覆介质的横波速度,vs2为下伏介质的横波速度,vp1为上覆介质的纵波速度,vp2为下伏介质的纵波速度,θ1和θ2分别表示纵波的反射角和透射角,和分别表示横波的反射角和透射角。
3.根据权利要求1所述的一种基于混合Markov邻域的叠前地震多参数同步反演非线性方法,其特征在于,所述假设条件包括:
令VS/VP为常数;
其中,VS为测井横波速度,VP为测井纵波速度;
所述背景趋势关系为:
lnVS=k·lnVP+kC+ΔLS,
lnρ=m·lnVP+mC+ΔLD,
其中,k表示测井纵横波速度的对数的拟合曲线的斜率,m表示测井纵波速度和测井密度的对数的拟合曲线的斜率,kC表示测井纵横波速度的对数的拟合曲线的截距,mC表示测井纵波速度与测井密度的对数的拟合曲线的截距,△LS和△LD分别是横波阻抗和密度的拟合偏差值,ρ为测井密度。
4.根据权利要求1所述的一种基于混合Markov邻域的叠前地震多参数同步反演非线性方法,其特征在于,所述目标函数为:
其中,J1为数据项,J2为先验项,J3为测井数据的约束项,Z为待反演的参数,D(Z)是Z的梯度值,σ为白噪声的标准方差,Y(θ)为地震响应值,W(θ)为角度子波,R(θ)为Zoeppritz方程中所求的纵波入射纵波反射的反射系数,φ为势函数,为模型参数的梯度值,δ为刻度参数,C1k为不同类型的阶团,C2k为垂直方向上的点集,λ1为先验项的加权参数,λ2为测井数据的约束项的加权参数,为测井纵波速度的梯度值,δP为测井纵波速度的刻度参数,为测井横波速度的梯度值,δS为测井横波速度的刻度参数,为测井密度的梯度值,δρ为测井密度的刻度参数。
5.根据权利要求1所述的一种基于混合Markov邻域的叠前地震多参数同步反演非线性方法,其特征在于,所述目标函数的求解过程包括:
在求解初期通过三阶邻域对目标函数进行求解;
在求解后期通过一阶邻域与三阶邻域的差值对目标函数进行求解。
6.根据权利要求5所述的一种基于混合Markov邻域的叠前地震多参数同步反演非线性方法,其特征在于,所述一阶邻域与三阶邻域的差值为:
其中,Δφ为一阶邻域与三阶邻域的差值,为一阶邻域,为三阶邻域。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于河海大学,未经河海大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202010765384.1/1.html,转载请声明来源钻瓜专利网。