一种基于范数方程提高分辨率的反褶积方法技术

本发明专利技术提出了一种基于范数方程提高分辨率的反褶积方法，本方法包括以下步骤：设定普通地震记录普通模型中的噪音为零，得到地震记录标准模型，通过标准模型，得到频率模型，并通过频率模型得到实际振幅模型，利用L1/L2范数方程，并求取实际振幅谱与标准振幅谱的误差能量，求取误差能量最小时的实际振幅谱，对实际振幅谱进行傅里叶反变换得到对应的时间函数，将时间函数带入标准模型，得到地表一致性反褶积因子，将反褶积因子与地震道褶积完成反褶积处理。使用本发明专利技术的优点在于，可以根据不同频带的信噪比求得相对稳定的反褶积算子，在计算过程中可以针对噪音振幅应用限制算子进行压制，能够有效压制高低频噪音，拓展有效频带，得到良好的反褶积结果。

A deconvolution method based on norm equation to improve resolution

【技术实现步骤摘要】
一种基于范数方程提高分辨率的反褶积方法
本专利技术涉及地震资料处理领域，具体涉及一种基于范数方程提高分辨率的反褶积方法。
技术介绍
高分辨率地震资料在油气勘探开发中发挥着越来越重要的作用，而提高资料分辨率的重要手段是反褶积，它是处理过程中的关键环节，其作用主要是压缩地震反射脉冲的长度，提高反射地震记录的分辨能力，并进一步估计地下反射界面的反射系数。反射系数序列中有波阻抗随时间变化的信息，这就提供了速度和密度随时间变化的信息，随之就可得到地层、岩性及构造在地下中间分布的信息。在有利条件下还可得到岩石孔隙率、渗透率、孔隙流体性质(油、气、水)乃至地层压力的信息。在地震资料处理过程中，常用的反褶积方法为地表一致性反褶积，该方法采用的是高斯-赛德尔迭代法，在算子分解过程中，考虑了炮点、检波点、共中心点和偏移距多个分量，但在运算过程中没有考虑非地表一致性因素，在实际资料中存在的非地表一致性噪音会影响反褶积算子的求解，然而在处理中并没有合适的手段去压制此类噪音，因此常规反褶积方法无法得到良好的反褶积效果，从而影响后续处理成像，因此需要对该方法进行改进，提高地震资料处理成像，得到高分辨率地震剖面。在地表一致性假设中，地震道被看成是炮点算子、检波点算子、反射系数算子(按CDP排列)和与炮检距有关的算子的褶积。在对数-傅里叶域中，褶积变成求和。对于给定的频率，其信号振幅的对数是炮点、检波点、共中心点和与炮检距有关项的和。我们认为，对地表同一位置，滤波作用与地震波的入射角无关，无论是浅、中、深层反射，其滤波作用均相同。因此，我们把实现这种反滤波功能的方法，称为“地表一致性反褶积”。在常规地表一致性反褶积计算过程中，首先建立符合地震记录的模型：xijk(t)＝ξijk(t)*si(t)*gj(t)*m(i+j)/2(t)*ok(t)+nijk(t)，其中si(t)表示第i炮位置的子波，gj(t)表示第j个检波点位置的子波，m(i+j)/2(t)表示共中心点的子波，ok(t)表示第k个偏移距对应的子波，xijk(t)、ξijk(t)、nijk(t)分别为对应的地震记录、反射系数和噪音。对于一个炮集记录，因为它们来自同一个激发源，所以炮集记录自相关函数的多道设计能较好地反映炮点子波的自相关函数的特性，定义其中Cijk表示的权系数，i0表示第i0炮。同理，第j0个检波点对应的子波自相关函数为共中心点对应的子波自相关函数为第k0个偏移距对应的子波函数为上述四个分量计算之后，对任意一道地震记录就可以得到其地表一致性自相关函数依据上式就可以求得最小相位子波，再选择适当的期望输出计算反褶积因子，与地震道褶积即可完成反褶积处理。常规地表一致性反褶积方法的假设条件中没有考虑实际资料里存在的非地表一致性问题，而实际资料中存在的该类干扰无法通过有效手段进行压制和去除。且常规地表一致性反褶积是采用高斯-赛德尔迭代法，没有考虑到非地表一致性的噪音影响，如果这些干扰超过了高斯误差的分布范围，在地表一致性分解时就会产生偏差，从而不能较好地提高地震资料的分辨率。而在实际资料中，存在的此类噪音会影响反褶积算子的求解，没法得到良好的效果。
技术实现思路
针对现有技术中所存在的上述技术问题，本专利技术提出了一种基于范数方程提高分辨率的反褶积方法，包括以下步骤：步骤1、建立符合地震记录的普通模型，步骤2、设定地震记录普通模型中的噪音为零，得到地震记录标准模型，步骤3、通过标准模型，得到频率模型，并通过频率模型得到实际振幅模型和实际振幅谱，步骤4、利用L1/L2范数方程，设定标准振幅谱，并求取实际振幅谱与标准振幅谱的误差能量，步骤5、求取误差能量最小时的实际振幅谱，对实际振幅谱进行傅里叶反变换得到对应的时间函数，步骤6、将时间函数带入标准模型，得到地表一致性反褶积因子，将反褶积因子与地震道褶积完成反褶积处理。在一个实施例中，步骤1中的地震记录普通模型为：xijk(t)＝ξijk(t)*si(t)*gj(t)*m(i+j)/2(t)*ok(t)+nijk(t)其中：si(t)为第i炮位置的子波，gj(t)为第j个检波点位置的子波，m(i+j)/2(t)为共中心点的子波，ok(t)为第k个偏移距对应的子波，xijk(t)为对应的地震记录、ξijk(t)为反射系数、nijk(t)为噪音。在一个实施例中，步骤2中的地震记录标准模型为：xijk(t)＝ξijk(t)*si(t)*gj(t)*m(i+j)/2(t)*ok(t)其中：si(t)为第i炮位置的子波，gj(t)为第j个检波点位置的子波，m(i+j)/2(t)为共中心点的子波，ok(t)为第k个偏移距对应的子波，xijk(t)为对应的地震记录、ξijk(t)为反射系数。在一个实施例中，步骤3具体包括以下步骤：步骤3.1、对标准模型进行傅里叶变换，转换到频率域，得到频率模型，步骤3.2、设定该模型对应的子波为最小相位，将频率模型分解到振幅谱，得到实际振幅谱，步骤3.3、对实际振幅谱取对数得到实际振幅模型。在一个实施例中，步骤3.3中的实际振幅模型为：lnAijk＝lnAs+lnAg+lnAm+lnAo其中：Aijk为实际振幅谱，As为该炮点位置的子波的振幅谱，Ag为对应的检波点位置子波的振幅谱，Am为共中心点子波的振幅谱，Ao为偏移距对应的子波的振幅谱。在一个实施例中，步骤4中的误差能量其中，Aijk为实际振幅谱，为模型振幅谱。在一个实施例中，步骤5具体包括以下步骤：步骤5.1、令求得各分量对数谱，其中AS为中心点子波的实际振幅谱，Ag为检波点子波实际振幅，Am为共中心点子波的实际振幅，Ao为偏移距对应子波的实际振幅，步骤5.2、当误差能量最小时，采用雅克比迭代法求取中心点子波、检波点子波、共中心点的子波、偏移距对应的子波的振幅谱，步骤5.3、对中心点子波、检波点子波、共中心点的子波、偏移距对应的子波的振幅谱进行傅里叶反变换得到对应的时间函数si(t)、gj(t)、m(i+j)/2(t)、ok(t)，其中si(t)、gj(t)、m(i+j)/2(t)、ok(t)，其中：si(t)为第i炮位置的子波，gj(t)为第j个检波点位置的子波，m(i+j)/2(t)为共中心点的子波，ok(t)为第k个偏移距对应的子波。与现有技术相比，本专利技术的优点在于，可以根据不同频带的信噪比求得相对稳定的反褶积算子，在计算过程中可以针对噪音振幅应用限制算子进行压制，能够有效压制高低频噪音，拓展有效频带，得到良好的反褶积结果。附图说明下面将结合附图来对本专利技术的优选实施例进行详细地描述，在图中：图1为本专利技术实施例的基于范数方程提高分辨率的反褶积方法的流程性框图；图2显示了现有技术中反褶积方法的地震偏移剖面；图3显示了应用本专利技术后的反褶积方法的地震偏移剖面；图4显示了采用现有技术与本专利技术两种方法偏移剖面的频谱分本文档来自技高网...

【技术保护点】
1.一种基于范数方程提高分辨率的反褶积方法，其特征在于，包括以下步骤：/n步骤1、建立符合地震记录的普通模型，/n步骤2、设定所述地震记录普通模型中的噪音为零，得到地震记录标准模型，/n步骤3、通过所述标准模型，得到频率模型，并通过所述频率模型得到实际振幅模型和实际振幅谱，/n步骤4、利用L1/L2范数方程，设定标准振幅谱，并求取所述实际振幅谱与所述标准振幅谱的误差能量，/n步骤5、求取误差能量最小时的实际振幅谱，对所述实际振幅谱进行傅里叶反变换得到对应的时间函数，/n步骤6、将所述时间函数带入所述标准模型，得到地表一致性反褶积因子，将所述反褶积因子与地震道褶积完成反褶积处理。/n

【技术特征摘要】
1.一种基于范数方程提高分辨率的反褶积方法，其特征在于，包括以下步骤：
步骤1、建立符合地震记录的普通模型，
步骤2、设定所述地震记录普通模型中的噪音为零，得到地震记录标准模型，
步骤3、通过所述标准模型，得到频率模型，并通过所述频率模型得到实际振幅模型和实际振幅谱，
步骤4、利用L1/L2范数方程，设定标准振幅谱，并求取所述实际振幅谱与所述标准振幅谱的误差能量，
步骤5、求取误差能量最小时的实际振幅谱，对所述实际振幅谱进行傅里叶反变换得到对应的时间函数，
步骤6、将所述时间函数带入所述标准模型，得到地表一致性反褶积因子，将所述反褶积因子与地震道褶积完成反褶积处理。


2.根据权利要求1所述的反褶积方法，其特征在于，步骤1中所述的地震记录普通模型为：
xijk(t)＝ξijk(t)*si(t)*gj(t)*m(i+j)/2(t)*ok(t)+nijk(t)
其中：si(t)为第i炮位置的子波，gj(t)为第j个检波点位置的子波，m(i+j)/2(t)为共中心点的子波，ok(t)为第k个偏移距对应的子波，xijk(t)为对应的地震记录、ξijk(t)为反射系数、nijk(t)为噪音。


3.根据权利要求2所述的反褶积方法，其特征在于，步骤2中所述的地震记录标准模型为：
xijk(t)＝ξijk(t)*si(t)*gj(t)*m(i+j)/2(t)*ok(t)
其中：si(t)为第i炮位置的子波，gj(t)为第j个检波点位置的子波，m(i+j)/2(t)为共中心点的子波，ok(t)为第k个偏移距对应的子波，xijk(t)为对应的地震记录、ξijk(t)为反射系数。

【专利技术属性】
技术研发人员：金红娣，刘晨，白帆，李芦茜，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京;11