[发明专利]一种地震资料反演的方法及装置、设备

摘要

本发明涉及一种地震资料反演的方法及装置、设备。该方法建立基于l2范数约束的弹性阻抗反演最优化模型。提出先借助阈值滤波子和L‑曲线准则确定最优正则参数，然后利用Tikhonov滤波子求解最优化模型的极小解的方法，能有效地提高计算精度以及解的稳定性。通过低信噪比理论模型和实际地震资料对算法进行测试。低信噪比理论模型的反演结果与模型正演情况基本吻合；而基于实际地震资料的弹性阻抗反演结果与该地区岩石物理分析结果相符合，进一步获得的脆性指数分布的预测结果也与钻井吻合。该算法具有良好的稳定性和抗噪性，有实际应用推广的前景。

主权项

1.一种地震资料反演的双滤子正则化方法，具体步骤如下：步骤1：对投影模型空间到观测空间的算子矩阵W进行奇异值分解，[U，∑，V]＝svd(W)，其中Σ＝diag(σ₁，σ₂，…σ_n)是矩阵奇异值，呈降序排列；U、V分别是m×m和n×n阶正交矩阵，U＝(u₁，u₂，…，u_m)，V＝(v₁，v₂，…，v_n)，即秩为r的m*n阶矩阵，r是W的数值秩，W的奇异系统可以记为{σ_i，u_i，v_j}；步骤2：计算解的范数和残差的范数其中，σ是矩阵W的奇异值，α是正则参数，d＝WR+N，式中，d＝(d₁，d₂，…，d_m)^T，R＝(R₁，R₂，…，R_n)^T，N＝(N₁，N₂，…N_m)^T；步骤3：根据步骤2中解的范数和残差的范数(log||Rα||2，log||WRα‑d||2)绘制关系曲线，计算最优正则解。