[发明专利]地球物理成像方法及装置

权利要求书

1.一种地球物理成像方法，应用于利用声波进行地球物理成像的情况下，包括：

获取震源函数，并构建所述震源函数的震源倒数函数；

利用所述震源倒数函数和空间波场参数，确定各点的反射系数；

根据确定出的各点的反射系数进行地球物理成像。

2.根据权利要求1所述的方法，其特征在于，构建所述震源函数的震源倒数函数，包括：若所述震源函数为f(t)，则所述震源倒数函数F(t)=1/f(t)。

3.根据权利要求2所述的方法，其特征在于，所述震源函数与所述震源倒数函数两者呈线性关系。

4.根据权利要求1至3任一项所述的方法，其特征在于，在任意时刻波场中，所述震源函数的下行波波场D(x,y,z;t)的取值等于所述震源倒数函数的下行波波场D'(x，y，z;t)取值的倒数。

5.根据权利要求4所述的方法，其特征在于，利用所述震源倒数函数以及收集的各点的波场参数，确定各点的反射系数，包括：

采用如下公式确定各点的反射系数：

I(x,y,z)=∫0tmaxU(x,y,z;t)D'(x,y,z;t),]]>

D'(x,y,z=0;t)=δ(x→-x→s)F(t)dt,]]>

其中，I(x,y,z)为坐标为（x，y，z）的点的反射系数，U(x,y,z;t)为所述震源函数的上行波波场，为空间向量，F(t)为所述震源函数的倒函数，dt为延拓时间步长，tmax为声波的最大旅行时间。

6.一种函数变换方法，包括：

获取震源函数；

在任意时刻波场中，查找等于所述震源函数的下行波波场D(x,y,z;t)的取值的倒数的下行波波场D'(x,y,z;t)；

根据所述D'(x,y,z;t)，构建所述震源函数的震源倒数函数。

7.一种地球物理成像装置，包括：

第一函数构建模块，配置为获取震源函数，并构建所述震源函数的震源倒数函数；

波场延拓模块，配置为利用所述震源倒数函数和空间波场参数，确定各点的反射系数；

成像模块，配置为根据确定出的各点的反射系数进行地球物理成像。