[发明专利]一种稀疏正则化的岩石物理弹性参数提取方法

权利要求书

1.一种稀疏正则化的岩石物理弹性参数提取方法，包括以下步骤：

1)叠前地震道集分小、中、大三个角度范围进行叠加，获得小、中两个角度部分叠加道集，经过小、中两个角度部分叠加道集反演得到小、中两个角度对应的弹性阻抗；

2)将小、中两个角度及小、中两个角度对应的弹性阻抗代入弹性阻抗方程得到弹性阻抗方程组：

1+tan2θ1-8Ksin2θ11-4Ksin2θ11+tan2θ2-8Ksin2θ21-4Ksin2θ2In(Vp)In(Vs)In(ρ)=In(EI(θ1))In(EI(θ2));]]>

3)令A=1+tan2θ1-8Ksin2θ11-4Ksin2θ11+tan2θ2-8Ksin2θ21-4Ksin2θ2;]]>

x=In(Vp)In(Vs)In(ρ);b=In(EI(θ1))In(EI(θ2));]]>则弹性阻抗方程组简化为欠定方程组：

Ax＝b；

4)对欠定方程组进行推导获得待求解线性方程组：

F(x,z,s)=Ax-bATz+s-C0Ds-βμe=000;]]>

式中，F表示对x,z,s的运算法则，_T表示转秩，z为x的对偶解，s为松弛因子，C₀＝[1…1]，D＝diag(x₁,x₂…x_N)，diag表示对角矩阵；β为衰减变量，0＜β＜1；μ也为引入的松弛因子，计算过程中趋向于0；e为单位矩阵；

5)对待求解线性方程组求解得到(x_k,z_k,s_k)^T；

6)根据In(Vp)In(Vs)In(ρ)=xk]]>计算获得弹性参数纵波速度Vp、横波速度Vs和密度ρ。