[发明专利]一种地震波场坡印廷矢量确定方法及系统

权利要求书

1.一种地震波场坡印廷矢量确定方法，其特征在于，包括：

对地震波场的时间导数的幅值进行调整，计算调整后的地震波场的时间导数的坡印廷矢量，得到第一坡印廷矢量；所述第一坡印廷矢量的幅值与第二坡印廷矢量的幅值的差值处于设定范围内；所述第二坡印廷矢量为所述地震波场的坡印廷矢量；

对所述第二坡印廷矢量和所述第一坡印廷矢量进行运算，得到所述地震波场的最终坡印廷矢量；

所述对地震波场的时间导数的幅值进行调整，计算调整后的地震波场的时间导数的坡印廷矢量，得到第一坡印廷矢量，具体包括：

计算所述地震波场的一阶时间导数；

计算所述地震波场的一阶时间导数的幅值调整因子；

采用所述幅值调整因子调整所述地震波场的一阶时间导数的幅值，得到调整时间导数；

计算所述调整时间导数的坡印廷矢量，得到所述第一坡印廷矢量。

2.根据权利要求1所述的一种地震波场坡印廷矢量确定方法，其特征在于，所述对所述第二坡印廷矢量和所述第一坡印廷矢量进行运算，得到所述地震波场的最终坡印廷矢量，具体包括：

将所述第二坡印廷矢量和所述第一坡印廷矢量相加，得到所述最终坡印廷矢量。

3.根据权利要求1所述的一种地震波场坡印廷矢量确定方法，其特征在于，所述计算所述地震波场的一阶时间导数，具体包括：

采用有限差分法计算所述地震波场的一阶时间导数。

4.根据权利要求3所述的一种地震波场坡印廷矢量确定方法，其特征在于，所述幅值调整因子的计算公式为：

其中，a为幅值调整因子；f_d为地震波场的主频；Δt为有限差分的时间间隔。

5.根据权利要求1所述的一种地震波场坡印廷矢量确定方法，其特征在于，所述调整时间导数的计算公式为：

其中，为t时刻的调整时间导数；a为幅值调整因子；t为时间；u^t为t时刻的地震波场；为对u^t求时间偏导数得到的一阶时间导数。

6.根据权利要求1所述的一种地震波场坡印廷矢量确定方法，其特征在于，所述第一坡印廷矢量的计算公式为：

其中，为第一坡印廷矢量；为t时刻的调整时间导数；t为时间；为的梯度。

7.根据权利要求1所述的一种地震波场坡印廷矢量确定方法，其特征在于，所述第一坡印廷矢量的计算公式为：

其中，为第一坡印廷矢量；为t时刻的调整时间导数；t为时间；为对求梯度；u^t为t时刻的地震波场；f_d为地震波场的主频；Δt为有限差分的时间间隔。

8.根据权利要求3所述的一种地震波场坡印廷矢量确定方法，其特征在于，所述一阶时间导数的计算公式为：

其中，u^t为t时刻的地震波场；t为时间；Δt为有限差分的时间间隔；u^t+Δt为t+Δt时刻的地震波场；u^t-Δt为t-Δt时刻的地震波场。

9.一种地震波场坡印廷矢量确定系统，其特征在于，包括：

幅值调整模块，用于对地震波场的时间导数的幅值进行调整，计算调整后的地震波场的时间导数的坡印廷矢量，得到第一坡印廷矢量；所述第一坡印廷矢量的幅值与第二坡印廷矢量的幅值的差值处于设定范围内；所述第二坡印廷矢量为所述地震波场的坡印廷矢量；

坡印廷矢量运算模块，用于对所述第二坡印廷矢量和所述第一坡印廷矢量进行运算，得到所述地震波场的最终坡印廷矢量；

所述对地震波场的时间导数的幅值进行调整，计算调整后的地震波场的时间导数的坡印廷矢量，得到第一坡印廷矢量，具体包括：

计算所述地震波场的一阶时间导数；

计算所述地震波场的一阶时间导数的幅值调整因子；

采用所述幅值调整因子调整所述地震波场的一阶时间导数的幅值，得到调整时间导数；

计算所述调整时间导数的坡印廷矢量，得到所述第一坡印廷矢量。