自由表面多次波的压制方法和装置制造方法及图纸

本发明专利技术提供了一种自由表面多次波的压制方法和装置，其中，该方法包括：将所获取的第一单炮地震数据划分为反射部分以及非反射部分；提取反射部分中的第一自由表面多次波；提取非反射部分中的第二自由表面多次波；根据第一单炮地震数据、第一自由表面多次波以及第二自由表面多次波，计算得到压制自由表面多次波之后的第二单炮地震数据。在本发明专利技术实施例中，所得到的压制自由表面多次波之后的第二单炮地震数据精度更高。

The suppression method and device of the free surface multiple wave

The invention provides a method and device for the suppression of a free surface multiple wave, in which the method includes: dividing the obtained first single gun seismic data into reflection and non reflective parts; extracting the first free surface multiple waves in the reflected part; extracting the second free surface multiple waves in the non reflective part; According to the seismic data of the first single gun, the first free surface multiple wave and the second free surface multiple wave, the second single gun seismic data after the suppression of the free surface multiple waves are obtained. In the embodiment of the invention, the accuracy of the second single shot seismic data obtained after suppressing the free surface multiple wave is higher.

【技术实现步骤摘要】
自由表面多次波的压制方法和装置
本专利技术涉及地质勘探
，特别涉及一种自由表面多次波的压制方法和装置。
技术介绍
在地震资料处理过程中，多次波通常被视为一种相干噪声，特别是针对海上采集的地震数据，多次波发育尤为严重，处理过程中，如果不能采取有效的手段对多次波予以消除，最终的成像结果往往会造成成果数据出现地质假象，甚至会降低地震资料信噪比以及分辨率。多次波根据其产生机理，其种类非常复杂，针对海洋资料二样，能量最强的往往是自由表面多次波。通常，通过采用多次波压制方案，可以有效的消除自由表面反射率对地震数据的影响。目前，常规多次波压制方法通常只需要提供震源处地震资料的波形特征及自由表面反射率，并不需要提取详细的地下地质信息，而随着地震勘探技术的不断进步，自由表面多次波压制方法也逐渐演变成了适应性的反演处问题，即：首先，估算震源与自由表面的反射率；其次，将不含多次的数据视为反演剩余量。但是此方法受计算成本及精度的限制，无法应用于三维地震数据。
技术实现思路
本专利技术提供了一种自由表面多次波的压制方法和装置，以解决现有技术无法满足三维地震数据压制多次波的问题。本专利技术实施例提供了一种自由表面多次波的压制方法，可以包括：将所获取的第一单炮地震数据划分为反射部分以及非反射部分；提取所述反射部分中的第一自由表面多次波；提取所述非反射部分中的第二自由表面多次波；根据所述第一单炮地震数据、所述第一自由表面多次波以及所述第二自由表面多次波，计算得到压制自由表面多次波之后的第二单炮地震数据。在一个实施例中，所述提取所述反射部分中的第一自由表面多次波包括：获取所述反射部分中入射波为反射波、出射波也为反射波的第三自由表面多次波；获取所述反射部分中入射波为反射波、出射波为非反射波的第四自由表面多次波；将所述第三自由表面多次波以及所述第四自由表面多次波的和作为所述反射波中的第一自由表面多次波。在一个实施例中，按照以下公式提取所述反射波中的第一自由表面多次波：其中，表示第n次迭代后得到的第一自由表面多次波，表示第n次迭代后得到的第三自由表面多次波，表示第n次迭代后得到的第四自由表面多次波。在一个实施例中，所述提取所述反射部分中的第一自由表面多次波包括：获取所述反射部分中第n次迭代去除自由表面多次波后的波形数据，其中，n为正整数；计算所述反射部分中第n次迭代去除自由表面多次波后的波形数据与所述反射部分的乘积，将所得到的乘积作为所述第三自由表面多次波；计算所述反射部分中第n次迭代去除自由表面多次波后的波形数据与所述非反射部分的乘积，将所得到的乘积作为所述第四自由表面多次波；计算所述第三自由表面多次波以及所述第四自由表面多次波的和，将所得到的和作为所述反射部分中的第一自由表面多次波。在一个实施例中，按照以下公式提取所述反射部分中的第一自由表面多次波：其中，表示第n次迭代后得到的第一自由表面多次波，表示所述反射部分中第n次迭代去除自由表面多次波后的波形数据，Dr表示所述反射部分，Ds表示所述非反射部分。在一个实施例中，所述提取所述非反射部分中的第二自由表面多次波包括：获取所述非反射部分中入射波为非反射波、出射波也为非反射波的第五自由表面多次波；获取所述非反射部分中入射波为非反射波、出射波为反射波的第六自由表面多次波；将所述第五自由表面多次波以及所述第六自由表面多次波的和作为所述非反射波中的第二自由表面多次波。在一个实施例中，按照以下公式提取所述非反射部分中的第二自由表面多次波：其中，表示第n次迭代后得到的第二自由表面多次波，表示第n次迭代后得到的第五自由表面多次波，表示第n次迭代后得到的第六自由表面多次波。在一个实施例中，所述提取所述非反射部分中的第二自由表面多次波包括：获取所述非反射部分中第n次迭代去除自由表面多次波后的波形数据，其中，n为正整数；计算所述非反射部分中第n次迭代去除自由表面多次波后的波形数据与所述非反射部分的乘积，将所得到的乘积作为所述第五自由表面多次波；计算所述非反射部分中第n次迭代去除自由表面多次波后的波形数据与所述非反射部分的乘积，将所得到的乘积作为所述第六自由表面多次波；计算所述第五自由表面多次波以及所述第六自由表面多次波的和，将所得到的和作为所述非反射部分中的第二自由表面多次波。在一个实施例中，按照以下公式提取所述非反射部分中的第二自由表面多次波：其中，表示第n次迭代后得到的第二自由表面多次波，表示所述非反射部分中第n次迭代去除自由表面多次波后的波形数据，Dr表示所述反射部分，Ds表示所述非反射部分。在一个实施例中，按照以下公式计算得到压制自由表面多次波之后的第二单炮地震数据：上式中，上式中，表示第n次迭代压制自由表面多次波之后的第二单炮地震数据，表示第n次迭代后得到的第一自由表面多次波，表示第n次迭代后得到的第二自由表面多次波，A表示表面算子，D表示所述第一单炮地震数据。在一个实施例中，根据所述第一单炮地震数据的同相轴形态差异，将所获取的第一单炮地震数据划分为反射部分以及非反射部分。在一个实施例中，在将所获取的第一单炮地震数据划分为反射部分以及非反射部分之前，所述方法还包括：对所述第一单炮地震数据进行预处理，所述预处理包括以下至少之一：环境校正、静校正、去噪、动校正拉伸。在一个实施例中，在计算得到压制自由表面多次波之后的第二单炮地震数据之后，所述方法还包括：根据所述第一单炮地震数据的叠加速度，绘制所述第二单炮地震数据的叠加剖面。本专利技术实施例还提供了一种自由表面多次波的压制装置，可以包括：第一数据划分模块，用于将所获取的第一单炮地震数据划分为反射部分以及非反射部分；第一波形提取模块，用于提取所述反射部分中的第一自由表面多次波；第二波形提取模块，用于提取所述非反射部分中的第二自由表面多次波；第二数据划分模块，用于根据所述第一单炮地震数据、所述第一自由表面多次波以及所述第二自由表面多次波，计算得到压制自由表面多次波之后的第二单炮地震数据。在本专利技术实施例中，先将所获取的第一单炮地震数据划分为反射部分及非反射部分等两部分；并提取反射部分中所含第一自由表面多次波以及非反射部分中的第二自由表面多次波；最后，可以根据单炮地震数据、第一自由表面多次波以及第二自由表面多次波，计算得到压制自由表面多次波之后的第二单炮地震数据。由于本申请所提出的方法先将第一单炮地震数据中的反射部分和非反射部分进行分离，再分别从反射部分及非反射部分中提取自由表面多次波，解决了现有技术中自由表面多次波压制法仅可有效预测二维地震数据反射部分多次波，而无法识别分辨三维地震数据中的非反射部分所含的多次波的问题。附图说明为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。图1是本申请提供的一种自由表面多次波的压制方法流程图；图2是本申请提供的某地区原始单炮数据示意图；图3是采用现有技术进行自由表面多次波压制后所得到的地震数据示意图；图4是采用本申请所提出的自由表面多次波的压制方法所得到的地震数据示意图；图5是采用现有技术进行自由表面多次波压制后所得本文档来自技高网...

【技术保护点】
一种自由表面多次波的压制方法，其特征在于，包括：将所获取的第一单炮地震数据划分为反射部分以及非反射部分；提取所述反射部分中的第一自由表面多次波；提取所述非反射部分中的第二自由表面多次波；根据所述第一单炮地震数据、所述第一自由表面多次波以及所述第二自由表面多次波，计算得到压制自由表面多次波之后的第二单炮地震数据。

【技术特征摘要】
1.一种自由表面多次波的压制方法，其特征在于，包括：将所获取的第一单炮地震数据划分为反射部分以及非反射部分；提取所述反射部分中的第一自由表面多次波；提取所述非反射部分中的第二自由表面多次波；根据所述第一单炮地震数据、所述第一自由表面多次波以及所述第二自由表面多次波，计算得到压制自由表面多次波之后的第二单炮地震数据。2.如权利要求1所述的方法，其特征在于，所述提取所述反射部分中的第一自由表面多次波包括：获取所述反射部分中入射波为反射波、出射波也为反射波的第三自由表面多次波；获取所述反射部分中入射波为反射波、出射波为非反射波的第四自由表面多次波；将所述第三自由表面多次波以及所述第四自由表面多次波的和作为所述反射波中的第一自由表面多次波。3.如权利要求2所述的方法，其特征在于，按照以下公式提取所述反射波中的第一自由表面多次波：其中，表示第n次迭代后得到的第一自由表面多次波，表示第n次迭代后得到的第三自由表面多次波，表示第n次迭代后得到的第四自由表面多次波。4.如权利要求2所述的方法，其特征在于，所述提取所述反射部分中的第一自由表面多次波包括：获取所述反射部分中第n次迭代去除自由表面多次波后的波形数据，其中，n为正整数；计算所述反射部分中第n次迭代去除自由表面多次波后的波形数据与所述反射部分的乘积，将所得到的乘积作为所述第三自由表面多次波；计算所述反射部分中第n次迭代去除自由表面多次波后的波形数据与所述非反射部分的乘积，将所得到的乘积作为所述第四自由表面多次波；计算所述第三自由表面多次波以及所述第四自由表面多次波的和，将所得到的和作为所述反射部分中的第一自由表面多次波。5.如权利要求4所述的方法，其特征在于，按照以下公式提取所述反射部分中的第一自由表面多次波：其中，表示第n次迭代后得到的第一自由表面多次波，表示所述反射部分中第n次迭代去除自由表面多次波后的波形数据，Dr表示所述反射部分，Ds表示所述非反射部分。6.如权利要求1所述的方法，其特征在于，所述提取所述非反射部分中的第二自由表面多次波包括：获取所述非反射部分中入射波为非反射波、出射波也为非反射波的第五自由表面多次波；获取所述非反射部分中入射波为非反射波、出射波为反射波的第六自由表面多次波；将所述第五自由表面多次波以及所述第六自由表面多次波的和作为所述非反射波中的第二自由表面多次波。7.如权利要求6所述的方法，其特征在于，按照以下公式提取所述非反射部分中的第二自由表面多次波：其中，表示...

【专利技术属性】
技术研发人员：苏勤，徐兴荣，曾华会，张小美，王建华，肖明图，吴杰，张涛，赵伟，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：北京,11