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

说明书

技术领域

本发明涉及地球物理成像领域，具体涉及一种地球物理成像方法及装置。

背景技术

地球结构较为复杂，为了方便人类执行任务，通常会将执行任务所在的物理结构绘制成像，以供参考。

大陆表层的物理结构是可见的，但是，地球还具有海洋、冰山、地心熔浆层等多种结构，并不是所有的物理结构均是可见的。但是，人类需要到各种环境下执行任务，例如，海底光缆的铺设，就需要明确海底的结构。因此，需要对地球进行物理成像。

目前采用的较多的地球物理成像方法是地震波成像方式，采用地震勘探收集地震波的参数，其原理为：波遇到速度差异界面产生反射和透射，震源激发后，波场的任何一个点都产生向上和向下两个方向传播的波。在地下所有的点，如果能得到两个方向波，就能求出所有点的反射系数，将各点的反射系数绘成剖面，就是一张成像剖图。

图1示出了根据背景技术的反射系数的获取方式的示意图。参见图1，每个点都能够得到向上的波Pb以及向下的波Pf，该点的反射系数R=Pb/Pf。

由于检波器只能放置在地面，因此在一定的高度下可能会检测不到向下的波Pf的情况，此时，Pf为0，则该点的R就无法获取，进而无法生成相应的图像。

发明内容

鉴于上述问题，提出了本发明以便提供一种克服上述问题或者至少部分地解决上述问题的地球物理成像方法和相应的装置。

依据本发明的一个方面，提供了一种地球物理成像方法，应用于利用声波进行地球物理成像的情况下，包括：

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

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

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

可选地，构建所述震源函数的震源倒数函数，包括：若所述震源函数为f(t)，则所述震源倒数函数F(t)=1/f(t)。

可选地，所述震源函数与所述震源倒数函数两者呈线性关系。

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

可选地，利用所述震源倒数函数以及收集的各点的波场参数，确定各点的反射系数，包括：

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

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,]]>