[发明专利]用于各向异性参数反演的敏感度矩阵的计算方法及系统

权利要求书

1.一种用于各向异性参数反演的敏感度矩阵的计算方法，包括以下步骤：

基于实际测井数据建立纵波反射系数与纵波阻抗、横波阻抗、密度以及各向异性参数的泛函关系；

基于褶积原理，利用所建立的泛函关系与给定的地震子波得到模拟叠前方位角度道集；

根据实际叠前方位角度道集和模拟叠前方位角度道集建立误差函数；

对误差函数求偏导数得到用于各向异性参数反演的敏感度矩阵。

2.根据权利要求1所述的方法，其特征在于，根据以下表达式建立所述泛函关系：

r_pj＝c₁(lnZ_pj+1-lnZ_pj)+c₂(lnZ_sj+1-lnZ_sj)+c₃(lnρ_j+1-lnρ_j)+c₄(δ_j+1-δ_j)+c₅(ε_j+1-ε_j)+c₆(γ_j+1-γ_j)

式中，r_pj表示第j个界面的纵波反射系数，j表示界面的个数；Z_pj为第j个界面的纵波阻抗，Z_sj为第j个界面的横波阻抗，ρ_j为第j个界面的密度；δ_j、ε_j和γ_j分别为第j个界面的各向异性参数；ln表示取自然对数；系数c₁、c₂、c₃、c₄、c₅以及c₆分别为：

c₁＝(1+tan²θ)/2

c₂＝-4(ksinθ)²

c₃＝tan²θ/2-2(ksinθ)²

其中，θ和分别为入射角和方位角；k为横波速度与纵波速度的比值。

3.根据权利要求2所述的方法，其特征在于，在建立所述误差函数的步骤中，

先由模拟叠前方位角度道集中的第i个样点值和实际叠前方位角度道集中的第i个样点值建立方位角为入射角为θ时的模拟叠前方位角度道集和实际叠前方位角度道集的误差函数；

再按照先入射角后方位角的顺序将对应不同方位角和入射角的误差函数逐一累加得到实际叠前方位角度道集和模拟叠前方位角度道集的误差函数。

4.根据权利要求3所述的方法，其特征在于，根据以下表达式建立实际叠前方位角度道集和模拟叠前方位角度道集的误差函数F：

式中，w为给定的地震子波；是实际叠前方位角度道集的第i个样点值；m为界面数，n为样点数；P1和P2为方位角的起始值和终止值，Q1和Q2为入射角的起始值和终止值为方位角为入射角为θ时的模拟叠前方位角度道集和实际叠前方位角度道集的误差函数，的表达式为：