[发明专利]粘声介质叠前时间偏移计算方法及装置

权利要求书

1.一种粘声介质叠前时间偏移计算方法，其特征在于，包括：

根据地震波传播的最大旅行时、最小等效Q因子以及预设的变量R的计算公式计算出变量R的最大值；

对所述变量R的最大值进行离散化处理生成常数序列；

通过快速傅里叶变换把地震道数据转化到频率域得到地震道数据频率域表达式；

分别采用所述常数序列中的各常数对所述地震道数据频率域表达式进行补偿，得到每个所述常数各自对应的补偿函数表达式；

采用快速反傅里叶变换把每个所述补偿函数表达式转化到时间域，得到时间域地震道序列，其中，所述常数序列中的各常数与所述时间域地震道序列中的各时间域地震道数据一一对应；

从用户提供的Q场数据体中获取成像样点位置的等效Q因子和成像样点位置的旅行时，并计算出成像样点位置对应的变量R；

确定成像样点位置对应的变量R在所述常数序列中的位置，确定所述常数序列中与该位置相邻的两个常数，并确定该相邻的两个常数各自对应的时间域地震道数据；

根据该相邻的两个常数以及该相邻的两个常数各自对应的时间域地震道数据进行横向加权Kirchhoff积分法求和，以完成所有地震道数据的偏移计算。

2.根据权利要求1所述的粘声介质叠前时间偏移计算方法，其特征在于，所述变量R为地震波传播的旅行时与等效Q因子的商。

3.根据权利要求1所述的粘声介质叠前时间偏移计算方法，其特征在于，所述分别采用所述常数序列中的各常数对所述地震道数据频率域表达式进行补偿，得到每个所述常数各自对应的补偿函数表达式，具体为：

将各常数分别代入以下公式对所述地震道数据频率域表达式进行补偿，得到每个所述常数各自对应的补偿函数表达式：

其中，G′₁(ω)为补偿函数表达式，F(ω)为地震道数据频率域表达式，R表示所述常数序列中的各常数。

4.根据权利要求1所述的粘声介质叠前时间偏移计算方法，其特征在于，所述对所述变量R的最大值进行离散化处理生成常数序列，包括：

对所述变量R的最大值按照步长8进行离散化处理，生成常数序列。

5.一种粘声介质叠前时间偏移计算装置，其特征在于，包括：

变量R最大值计算单元，用于根据地震波传播的最大旅行时、最小等效Q因子以及预设的变量R的计算公式计算出变量R的最大值；

常数序列生成单元，用于对所述变量R的最大值进行离散化处理生成常数序列；

快速傅里叶变换单元，用于通过快速傅里叶变换把地震道数据转化到频率域得到地震道数据频率域表达式；

补偿单元，用于分别采用所述常数序列中的各常数对所述地震道数据频率域表达式进行补偿，得到每个所述常数各自对应的补偿函数表达式；

快速反傅里叶变换单元，用于采用快速反傅里叶变换把每个所述补偿函数表达式转化到时间域，得到时间域地震道序列，其中，所述常数序列中的各常数与所述时间域地震道序列中的各时间域地震道数据一一对应；

成像样点位置数据获取单元，用于从用户提供的Q场数据体中获取成像样点位置的等效Q因子和成像样点位置的旅行时，并计算出成像样点位置对应的变量R；

确定单元，用于确定成像样点位置对应的变量R在所述常数序列中的位置，确定所述常数序列中与该位置相邻的两个常数，并确定该相邻的两个常数各自对应的时间域地震道数据；

偏移计算单元，用于根据该相邻的两个常数以及该相邻的两个常数各自对应的时间域地震道数据进行横向加权Kirchhoff积分法求和，以完成所有地震道数据的偏移计算。

6.根据权利要求5所述的粘声介质叠前时间偏移计算装置，其特征在于，所述变量R为地震波传播的旅行时与等效Q因子的商。

7.根据权利要求5所述的粘声介质叠前时间偏移计算装置，其特征在于，所述补偿单元，具体用于将各常数分别代入以下公式对所述地震道数据频率域表达式进行补偿，得到每个所述常数各自对应的补偿函数表达式：

其中，G′₁(ω)为补偿函数表达式，F(ω)为地震道数据频率域表达式，R表示所述常数序列中的各常数。