1.一种反演裂缝密度的方法,所述方法包括以下步骤:
基于地震资料获取相交测线的纵波反射振幅差值,并基于拟合模型获得纵波反射系数差值;
基于所述纵波反射系数差值和所述纵波反射振幅差值,构造目标函数fitness;以及
基于所述目标函数fitness进行反演计算以获得裂缝密度,
其中,所述目标函数fitness表示为:
(i=1,2,...,I;j=1,2,...,J;k=1,2,...,K),
其中,表示方位角;ηj表示相交测线的方位角差值;θk表示入射角;表示方位角分别为和入射角为θk时的纵波反射振幅差值;表示的最大值;表示反射系数差值或振幅差值最大时对应的入射角;Δγ、Δδ(V)与Δε(V)表示Thomsen参数;表示拟合模型中两条相交测线的纵波反射系数差值;表示的最大值;ωij表示权系数;I表示地震资料中方位角的测线总数量;J表示地震资料中方位角为时与其相交的测线的总数量;以及K表示地震资料中入射角θ的总数量。
2.根据权利要求1所述的反演裂缝密度的方法,其中,使用Rüger反射系数近似公式来计算所述纵波反射系数差值。
3.根据权利要求1所述的反演裂缝密度的方法,其中,使用观测地震数据中目的层位上的相交测线的反射振幅来计算所述纵波反射振幅差值。
4.根据权利要求1所述的反演裂缝密度的方法,其中,所述目标函数fitness中的权系数用于控制一组不同方位角的纵波反射振幅差值的可信程度。
5.根据权利要求1所述的反演裂缝密度的方法,其中,基于所述目标函数fitness进行反演计算以获得裂缝密度包括:
通过非线性反演算法获得目标函数fitness的最小值;
基于目标函数fitness的最小值获得Thomsen参数;以及
基于Thomsen参数与裂缝密度之间的关系和所获得的Thomsen参数进行计算以获得裂缝密度。
6.一种反演裂缝密度的系统,所述系统包括以下单元:
用于基于地震资料获取相交测线的纵波反射振幅差值,并基于拟合模型获得纵波反射系数差值的单元;
用于基于所述纵波反射系数差值和所述纵波反射振幅差值,构造目标函数fitness的单元;以及
用于基于所述目标函数fitness进行反演计算以获得裂缝密度的单元,
其中,所述目标函数fitness表示为:
(i=1,2,...,I;j=1,2,...,J;k=1,2,...,K),
其中,表示方位角;ηj表示相交测线的方位角差值;θk表示入射角;表示方位角分别为和反射角入射角为θk时的纵波反射振幅差 值;表示值的最大值;表示反射系数差值或振幅差值最大时对应的入射角;Δγ、Δδ(V)与Δε(V)表示Thomsen参数;表示拟合模型中两条相交测线的纵波反射系数差值;表示的最大值;ωij表示权系数;I表示地震资料中方位角的测线总数量;J表示地震资料中方位角为时与其相交的测线的总数量;以及K表示地震资料中入射角θ的总数量。
7.根据权利要求6所述的反演裂缝密度的系统,其中,使用Rüger反射系数近似公式来计算所述纵波反射系数差值。
8.根据权利要求6所述的反演裂缝密度的系统,其中,使用观测地震数据中目的层位上的相交测线的反射振幅来计算所述纵波反射振幅差值。
9.根据权利要求6所述的反演裂缝密度的系统,其中,所述目标函数fitness中的权系数用于控制一组不同方位角的纵波反射振幅差值的可信程度。
10.根据权利要求6所述的反演裂缝密度的系统,其中,基于所述目标函数fitness进行反演计算以获得裂缝密度包括:
通过非线性反演算法获得目标函数fitness的最小值;
基于目标函数fitness的最小值获得Thomsen参数;以及
基于Thomsen参数与裂缝密度之间的关系和所获得的Thomsen参数进行计算以获得裂缝密度。