本申请实施例提供的一种压制虚反射信号的方法及装置。所述方法包括：获取频率-波数域的海上地震数据；根据所述海上地震数据的分析处理结果，计算虚反射压制因子；根据所述虚反射压制因子计算虚反射压制算子；根据所述海上地震数据以及所述虚反射压制算子计算压制虚反射信号后的海上地震数据。利用所述本发明专利技术实施例提供的技术方案，计算得到比较准确的虚反射压制因子和稳定性较强的压制虚反射信号后的海上地震数据，可以压制虚反射信号，消除虚反射信号的干扰。
【技术实现步骤摘要】
本专利技术涉及海上地震资料处理
，尤其涉及一种压制虚反射信号的方法及 装置。
技术介绍
地球物理勘探
中，虚反射往往是影响海上地震资料品质的一个重要因 素。在海上进行地震勘探时，一般将震源和检波器放置于海面以下一定深度处，由于海水与 空气的波阻抗差异较大，往往在海水与空气之间形成较强的反射面，所述反射面构成自由 界面。震源激发的地震波经过所述自由界面反射后可以形成虚反射。所述虚反射的波形、 频率、视速度等参数一般和一次反射波相似，因此，所述虚反射可能会干扰一次反射波，造 成地震记录的低频响应，降低地震剖面的辨率。所述虚反射有时甚至会产生虚假同相轴，给 地震资料反演与解释造成较大的困难。因此，压制虚反射信号是海上地震资料处理最重要 的步骤之一。 现有技术中常用反褶积方法压制虚反射信号，具体可以设计虚反射滤波器，在时 间域利用预测反褶积来压制虚反射信号。所述反褶积方法假设地震记录是反射波与虚反射 滤波器褶积以及所述虚反射滤波器为最小相位，实际上所述虚反射滤波器的振幅谱存在一 系列的"零值点"，所以所述虚反射滤波器不是最小相位的，不满足反褶积的条件，因此时间 域反褶积方法不能完全消除虚反射。
技术实现思路
本申请的目的在于提供一种压制虚反射信号的方法及装置，可以压制虚反射信 号，消除虚反射信号的干扰。 本申请提供的一种压制虚反射信号的方法及装置是这样实现的： -种压制虚反射信号的方法，所述方法包括： 获取频率-波数域的海上地震数据； 基于检波器在海水中的深度、所述海上地震数据的频率、波数，计算频率-波数域 的虚反射压制因子； 根据所述虚反射压制因子计算虚反射压制算子； 根据所述虚反射压制算子对所述频率-波数域的海上地震数据进行虚反射信号 压制处理，得到虚反射信号压制后的海上地震数据。 可选的，在本专利技术的一个实施例中，所述基于检波器在海水中的深度、所述海上地 震数据的频率、波数，计算频率-波数域的虚反射压制因子包括： 采用下式计算所述虚反射压制因子： 式中，G(F，K)为虚反射压制因子在频率-波数域的值，F和K分别为所述海上地 震数据的频率和波数，R为空气和海面间自由界面的反射系数，i为虚数单位，z为检波器在 海水中的深度，C为海水速度。 可选的，在本专利技术的一个实施例中，所述根据所述虚反射压制因子对所述频 率-波数域的海上地震数据进行虚反射信号压制处理包括： 采用下式计算虚反射压制后的海上地震数据： 式中，A(F，K)为虚反射压制算子在频率-波数域的值，G(F，K)为虚反射压制因子 在频率-波数域的值，G丨厂，幻为 G(F，κ)的共辄。 可选的，在本专利技术的一个实施例中，所述根据所述虚反射压制算子对所述频 率-波数域的海上地震数据进行虚反射信号压制处理，包括： 将所述频率-波数域的海上地震数据和所述虚反射压制算子做乘积运算。 可选的，在本专利技术的一个实施例中，所述根据所述虚反射压制因子对所述频 率-波数域的海上地震数据进行虚反射信号压制处理包括： 采用下式计算虚反射压制后的海上地震数据： P (F, Κ) = D (F, k) 式中，P(F，K)为压制虚反射信号后的频率-波数域的海上地震数据，D(F，K)为压 制虚反射信号前的频率-波数域的海上地震数据，A (F，K)为虚反射压制算子在频率-波数 域的值，white为白噪声系数。 可选的，在本专利技术的一个实施例中，所述白噪声系数的计算公式如下式所示： 式中，white为白噪声系数，F为所述海上地震数据的频率。 可选的，在本专利技术的一个实施例中，所述获取频率-波数域的海上地震数据，包 括： 获取时间-空间域的海上地震数据； 对所述时间-空间域的海上地震数据进行二维傅里叶变换。 -种压制虚反射信号的装置，所述装置包括： 地震数据获取单元，用于获取频率-波数域的海上地震数据； 压制因子计算单元，用于根据所述海上地震数据的分析处理结果，计算虚反射压 制因子； 虚反射压制算子计算单元，用于根据所述虚反射压制因子计算虚反射压制算子； 压制后地震数据计算单元，用于根据所述海上地震数据以及所述虚反射压制算子 计算压制虚反射信号后的海上地震数据。 可选的，在本专利技术的一个实施例中，所述压制因子计算单元包括： 所述虚反射压制因子的计算公式包括： 式中，G(F，K)为虚反射压制因子在频率-波数域的值，F和K分别为海上地震数 据的频率和波数，R为空气和海面间自由界面的反射系数，i为虚数单位，ζ为检波器在海水 中的深度，c为海水速度。 可选的，在本专利技术的一个实施例中，所述虚反射压制算子计算单元包括： 所述虚反射压制算子的计算公式包括： 式中，A(F，K)为虚反射压制算子在频率-波数域的值，G(F，K)为虚反射压制因子 在频率-波数域的值，Z 7. /〇为G (F，K)的共辄。 可选的，在本专利技术的一个实施例中，所述压制后地震数据计算单元包括： 计算子单元，用于将所述频率-波数域的海上地震数据和所述虚反射压制算子做 乘积运算。 可选的，在本专利技术的一个实施例中，所述压制后地震数据计算单元包括： 所述虚反射压制因子的计算公式包括： P (F, K) = D (F, k) 式中，P(F，K)为压制虚反射信号后的频率-波数域的海上地震数据，D(F，K)为压 制虚反射信号前的频率-波数域的海上地震数据，A (F，K)为虚反射压制算子在频率-波数 域的值，white为白噪声系数。 可选的，在本专利技术的一个实施例中，所述白噪声系数的计算公式如下式所示： 式中，white为白噪声系数，F为所述海上地震数据的频率。 可选的，在本专利技术的一个实施例中，所述地震数据获取单元，包括： 时深域地震数据获取单元，用于获取时间-空间域的海上地震数据； 傅里叶变换单元，用于对所述时间-空间域的海上地震数据进行二维傅里叶变 换。 由此可见，本专利技术的一种压制虚反射信号的方法和装置的技术方案通过对频 率-波数域内的海当前第1页1 2 3 4 本文档来自技高网...
【技术保护点】
一种压制虚反射信号的方法，其特征在于，所述方法包括：获取频率‑波数域的海上地震数据；基于检波器在海水中的深度、所述海上地震数据的频率、波数，计算频率‑波数域的虚反射压制因子；根据所述虚反射压制因子计算虚反射压制算子；根据所述虚反射压制算子对所述频率‑波数域的海上地震数据进行虚反射信号压制处理，得到虚反射信号压制后的海上地震数据。
【技术特征摘要】
【专利技术属性】
技术研发人员：马光凯，周铮铮，赵波，钱忠平，任晓乔，李合群，
申请(专利权)人：中国石油天然气集团公司，中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：北京;11
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。