1.一种地震波旅行时间计算方法,其特征在于,包括以下步骤:S1,确定基本旅行时间计算公式基本形式;所述基本旅行时间计算公式是根据泰勒级数原理展开式并利用各向异性介质模型得出基本旅行时间平方的表达式;所述基本旅行时间计算公式为:
2 2 4 6 8
T(x)=c0+c1x+c2x+c3x+c4x+...
2
c0=t
其中,x为偏移距;c0、c1、c2、c3、c4为泰勒级数的系数;t=T(0)为自激自收时间,v=vNMO表示NMO速度,Sk为非均质参数,k=2,3,4;
S2,根据所述泰勒级数原理和连续分式近似旅行时间平方公式对所述基本旅行时间计算公式进行更高阶的近似,得到下式:其中,参数B1是基于6阶泰勒展开和4阶泰勒展开的一种系数近似表示,参数B2是基于8阶泰勒展开和6阶泰勒展开的一种系数近似表示;
2
S3,以x为基数对更高阶的所述基本旅行时间计算公式中偏移距的4次方和6次方进行同时降幂处理,输出旅行时间;降幂处理后的所述基本旅行时间计算公式:其中,可取常用近似 来简化式中的参数c3、B1、B2,化简后的参数为:其中S2=1+8η,各向异性参数 ε表示各向异性强度、δ表示NMO速度和反射振幅的变化,参数ε和δ均为Thomsen各向异性参数。
2.根据权利要求1所述的一种地震波旅行时间计算方法,其特征在于,所述步骤S3之后还包括对降幂处理后的所述基本旅行时间计算公式结合近似补偿方法形成降幂补偿公式:其中,fb为初至时间,TDE(x)为降幂后的所述基本旅行时间计算公式计算出的旅行时间,n为模型层数。
3.一种地震波旅行时间计算系统,其特征在于,包括:计算优化模块,确定基本旅行时间计算公式基本形式;所述基本旅行时间计算公式是根据泰勒级数原理展开式并利用各向异性介质模型得出基本旅行时间平方的表达式;所述基本旅行时间计算公式为:
2 2 4 6 8
T(x)=c0+c1x+c2x+c3x+c4x+...
2
c0=t
其中,x为偏移距;c0、c1、c2、c3、c4为泰勒级数的系数;t=T(0)为自激自收时间,v=vNMO表示NMO速度,Sk为非均质参数,k=2,3,4;
据所述泰勒级数原理和连续分式近似旅行时间平方公式对所述基本旅行时间计算公式进行更高阶的近似,得到下式:其中,参数B1是基于6阶泰勒展开和4阶泰勒展开的一种系数近似表示,参数B2是基于8阶泰勒展开和6阶泰勒展开的一种系数近似表示;
2
计算输出模块,以x为基数对更高阶的所述基本旅行时间计算公式中偏移距的4次方和
6次方进行同时降幂处理,输出旅行时间;降幂处理后的所述基本旅行时间计算公式:其中,可取常用近似 来简化式中的参数c3、B1、B2,化简后的参数为:其中S2=1+8η,各向异性参数 ε表示各向异性强度、δ表示NMO速度和反射振幅的变化,参数ε和δ均为Thomsen各向异性参数。
4.根据权利要求3所述的一种地震波旅行时间计算系统,其特征在于,所述计算输出模块还包括,对降幂处理后的所述基本旅行时间计算公式结合近似补偿方法形成降幂补偿公式:其中,fb为初至时间,TDE(x)为降幂后的所述基本旅行时间计算公式计算出的旅行时间,n为模型层数。