过现有的反演算法得到地震数据低频波阻抗模型。具体地，可以利用井旁道提取的地震子波与低频地震数据进行反褶积运算，然后再根据反射系数与波阻抗的关系计算得到地震数据低频波阻抗模型AIMlsnu。。
[0032]步骤103:将所述测井低频波阻抗模型，分为频率范围不同的N个第一频率段波阻抗，其中，N为正整数；
[0033]具体的，N的取值可以为3或4。其中，当N的取值为3时，表示将测井低频波阻抗模型分为3个不同的频率范围，例如，如果总的频率范围为OHz到12Hz，则可以按照如下分法将整个频率范围分为如下三个小的频率范围:0Hz到3Hz、3Hz到6Hz、6Hz到12Hz。
[0034]然而，值得注意的是，上述分成的频率段的个数，以及具体每个频率段的频率范围仅是一种示意性描述，并不构成对本发明的不当限定，例如，N可以按照需要选取除了 3和4之外的值，具体的各频率段的范围也可以按照实际需要选取。同时，每个频率段之间可以采用等分的方式，也可以采用非等分的方式。
[0035]步骤104:将所述地震数据低频波阻抗模型，分为频率范围不同的N个第二频率段波阻抗，其中，所述第一频率段波阻抗与所述第二频率段波阻抗的频率范围一一对应；
[0036]S卩，在该步骤中，需要将地震数据低频波阻抗模型按照与测井低频波阻抗模型同样频率段进行划分，例如，如果测井低频波阻抗模型被分为3段，每个频率段波阻抗的频率范围为:0Hz到3Hz、3Hz到6Hz、6Hz到12Hz，那么地震数据低频波阻抗模型也应该分为3段，每个频率段波阻抗的频率范围为:0Hz到3Hz、3Hz到6Hz、6Hz到12Hz。
[0037]步骤105:对各个一一对应的第一频率段波阻抗和第二频率段波阻抗的数据进行比较，根据比较结果计算出所述N个第二频率段波阻抗中各个第二频率段波阻抗的补偿系数；
[0038]通常情况下，地震数据在采集和处理的过程中会在低频段有一定的损失，地震反演的相对波阻抗数据不能直接准确解释岩性，因此，需要对地震反演的结果进行低频补偿，进而得出补偿后的波阻抗模型。
[0039]由于地震数据的低频段波阻抗衰减较快，并且不同频率段所对应的波阻抗差不同。为了实现补偿，在本例中，分别选取测井低频波阻抗模型以及地震数据低频波阻抗模型的相同频率段的数据进行比较，可以计算得到每个频率段的地震波阻抗的补偿系数a。具体地，可以以测井低频波阻抗模型为基础值，地震数据低频波阻抗模型减去基础值可以得到差异值，这个差异值就是补偿系数。在计算补偿系数的过程中，可以预先对参与计算的数据进行检查，以确定数据中是否有奇异值，如果有奇异值，则先对这些奇异值进行平滑处理，然后再基于平滑处理后的数据计算补偿系数。
[0040]步骤106:根据各个第二频率段波阻抗的补偿系数，对所述N个第二频率段波阻抗中各个第二频率段波阻抗进行补偿，得到N个补偿后的第二频率段波阻抗；
[0041]在实际实施中，可以按照以下公式对N个不同频率段的地震数据低频波阻抗模型进行补偿:
[0042]AIi= a ^AIsi(公式 I)
[0043]其中，AI1表示地震数据低频波阻抗模型的第i个补偿后的频率段波阻抗，a i表示地震数据低频波阻抗模型的第i个频率段波阻抗的补偿系数，々^表示地震数据低频波阻抗模型的第i个补偿前的频率段波阻抗，i的取值为I到N。
[0044]S卩，可以按照上述公式I对每个不同频率段的地震数据低频波阻抗模型都进行补偿，以最终得到补偿后的地震数据低频波阻抗模型Al。
[0045]根据实际分析可知:在幅值上，补偿后的地震数据低频波阻抗模型与测井低频波阻抗模型相近；在空间上，补偿后的地震数据低频波阻抗模型具有更强的地震横向分辨率。
[0046]步骤107:对所述N个第一频率段波阻抗和所述N个补偿后的第二频率段波阻抗进行融合，得到融合后的低频波阻抗模型。
[0047]在实际执行的过程中，在得到补偿后的地震数据低频波阻抗模型后，要得到最终的用于地震反演岩性预测的低频波阻抗模型，还可以执行以下操作:
[0048]S1:对所述N个第一频率段波阻抗和所述N个补偿后的第二频率段波阻抗一对一进行比较；
[0049]仍旧以上述的N取3，各频率段波阻抗对应的频率范围为:0Hz到3Hz、3Hz到6Hz、6Hz到12Hz，为例进行说明。在该步骤中，相应地，就是将OHz到3Hz所对应的补偿后的地震数据低频波阻抗模型和OHz到3Hz所对应的测井低频波阻抗模型进行对比，将3Hz到6Hz所对应的补偿后的地震数据低频波阻抗模型和3Hz到6Hz所对应的测井低频波阻抗模型进行对比，将6Hz到12Hz所对应的补偿后的地震数据低频波阻抗模型和6Hz到12Hz所对应的测井低频波阻抗模型进行对比。
[0050]S2:将频率范围相同的第一频率段波阻抗和第二频率段波阻抗中精确度高的频率段波阻抗作为该频率范围的融合频率段波阻抗，得到N个融合频率段波阻抗；
[0051]例如，通过对比分析这3个频率段所得到的测井低频波阻抗模型以及补偿后的地震数据低频波阻抗模型发现:在OHz到3Hz的极低频段时，采用测井低频波阻抗模型可靠性更好一些，在3Hz到6Hz、6Hz到12Hz的中高频段时，采用补偿后的地震数据低频波阻抗模型岩性预测的效果更好。
[0052]那么，就可以将OHz到3Hz所对应的测井低频波阻抗模型，和3Hz到6Hz、6Hz到12Hz所对应的补偿后的地震数据低频波阻抗模型作为融合频率段波阻抗。
[0053]步骤3:将所述N个融合频率段波阻抗拼接成的低频波阻抗模型作为融合后的低频波阻抗模型。
[0054]具体的，在得到融合后的低频波阻抗模型之后，可以对上述3个融合频率段波阻抗分别进行傅里叶变换，再根据傅里叶变换后的结果，在频率域进行拼接，得到最终的用于地震反演岩性预测的低频波阻抗模型AIMlsnM+irell。该模型具有准确的测井曲线垂向识别能力以及地震横向预测能力。
[0055]为了对上例中的确定低频波阻抗的方法的效果进行说明，在本例中给出了一组实验验证数据，该组数据仍旧是以上述的N取3，各频率段范围为:0Hz到3Hz、3Hz到6Hz、6Hz到12Hz，为例进行的。如图2所示为采用测井曲线插值对碳酸盐岩储层分析得到的低频波阻抗模型示意图，图3所示为采用上述的确定低频波阻抗的方法得到的低频波阻抗模型示意图。对图2和图3进彳丁对比，可以看出，两幅图的幅值相近，然而，图3相对于图2的中尚频段在空间上的地震横向分辨率更强。
[0056]由此可见，采用本发明实施例得到的低频波阻抗模型中包含有更加丰富的波阻抗信息，因而具有更加精确的地震反演岩性预测能力。
[0057]基于同一发明构思，本发明实施例中还提供了一种确定低频波阻抗的装置，如下面的实施例所述。由于确定低频波阻抗的装置解决问题的原理与确定低频波阻抗的方法相似，因此确定低频波阻抗的装置的实施可以参见确定低频波阻抗的方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。图4是本发明实施例的确定低频波阻抗的装置的一种结构框图，如图4所示，包括:插值单元401、反演单元402、第一分频单元403、第二分频单元404、补偿系数计算单元405、补偿单元406和融合单元407，下面对该结构进行说明。
[0058]插值单元401，用于通过对待分析区域的测井曲线进行插值，得到测井低频波阻抗丰旲型；
[0059]反演单元402，用于通过对所述待分析区域的地震数据进行反演，得到地震数据低频波阻抗模型；
[0060]第一分频单元403，用于将所述测井低频波阻抗模型，分为频率范围不同的N个第一频率段波阻抗，其中，N为正整数；
[0061]第二分频单元404，用于将所述地震数据低频波阻抗模型，分为频率范围不同的N个第二频率段波阻抗，其中，所述第一频率段波阻抗与所述第二频率段波阻抗的频率范围
——对应；
[0062]补偿系数计算单元405，用于对各个一一对应的第一频率段波阻抗和第二频率段波阻抗的数据进行比较，根据比较结果计算出所述N个第二频率段波阻抗中各个第二频率段波阻抗的补偿系数；
[0063]补偿单元406，用于根据所述各个第二频率段波阻抗的补偿系数，对所述N个第二频率段波阻抗中各个第二频率段波阻抗进行补偿，得到N个补偿后的第二频率段波阻抗；
[0064]融合单元407，用于对所述N个第