[发明专利]基于三维稀疏聚焦域反演的噪声型数据一次波反演方法

说明书

本发明涉及一种基于三维稀疏聚焦域反演的噪声型数据一次波反演方法，先用线性算子相乘的形式表示地下全体波场、一次波波场、表层相关算子和多次波波场地下的物理传播关系；以线性算子构建基于L1范数一次波反演过程；将噪声型被动源数据一次波和表层相关算子初始值带入双凸优化反演过程，获得一次波估计初值；由一次波初值和波场间的关系来求取表层相关算子的初值；构建三维聚焦稀疏算子引入双凸优化一次波反演过程求取噪声型被动源数据一次波估计值；通过表层相关算子和噪声型被动源数据一次波交替迭代反演，获得噪声型被动源数据一次波估计值。本方法能在一次波估计过程中压制噪声型被动源数据的相干噪声、纠正远偏移距同相轴误差的一体化进程。

技术领域

本发明属于地球物理勘探技术领域，具体涉及一种噪声型被动源数据的一次波反演方法，特别涉及一种基于三维稀疏聚焦域反演的噪声型数据一次波反演方法。

背景技术

在地球物理勘探中，被动源地震勘探是一个新兴领域。它的提出与应用极大地的丰富了人们对地震波传播规律及性质的认识和研究，被认为是地球物理勘探的一大进步。不同于人工激发震源的主动源勘探，被动源地震勘探不需要人工震源，而是将地下的微振动或背景噪声作为震源。因此，被动源地震勘探能够极大的节约生产成本。

噪声型被动源地震探勘是被动源地震勘探中极为重要的一种勘探方式，通常背景噪声被视为噪声型被动源数据的震源。通过将检波器放置在地表，接收来自地下的波场从而形成被动源数据地震记录。

针对被动源数据地震记录，互相关算法是最常用的一种处理被动源数据的算法。该算法能够将被动源数据进行转化，形成地表激发，地表接收的虚拟炮集记录，从而获得地下的构造信息。同时，该虚拟炮集同传统主动源地震记录一样，不仅包含一次波信息也包含表层相关多次波信息。但是，这种方法依赖于地下均匀照明的假设，在激发角度和强度方面需要具有一致性。在实际采集中，这种假设很难实现，因此互相关算法很难提供具有准确振幅的近地表响应，所以获得结果也无法为后续的地震数据处理解释以及地下成像提供准确的地震资料。

为了克服上述问题，Van Groenestijn等人(2010.Estimation of primaries bysparse inversion from passive seismic data.GEOPHYSICS 75(4).SA61–SA69)改进了针对传统主动源地震记录的一次波估计方法EPSI(Estimation of primaries by sparseinversion)，使其能够使用于被动源地震数据多次波的压制。改进后的EPSI不需要满足地下均匀照明的假设就能获得具有准确振幅的一次波响应。然而，该方法需要人为设置时窗，事先设置准确的一次波响应初始值。程浩等人(2015.Estimating primaries frompassive seismic data.EXPLORATION GEOPHYSICS 46(2).184-191)提出将L1范数稀疏反演与EPSI相结合，避免了传统EPSI 需要人为设置时窗的缺点，同时也提高了一次波响应估计结果的精度。然而，现有的这些方法基本只适用于脉冲型被动源数据，同时也无法解决被动源数据远偏移距振幅不准的问题。 Gabriel等人(2015.Closed-loop surface-relatedmultiple elimination and its application to simultaneous datareconstruction.GEOPHYSICS 80(6).V189-V199)提出Closed-loop SRME，该方法用一种更为稳定的目标函数进行主动源地震记录一次波的估计，但是该方法并不能直接应用于噪声型被动源地震数据的处理。

噪声型被动源数据中存在大量的相干噪声，会对一次波反演产生严重的干扰，因而这些方法如果直接应用在噪声型被动源数据上难以获得理想的效果。因此，针对噪声型被动源数据一次波估计的研究也鲜有报道。

发明内容