CN105425289B - 确定低频波阻抗的方法和装置 - Google Patents

确定低频波阻抗的方法和装置 Download PDF

Info

Publication number
CN105425289B
CN105425289B CN201510717763.2A CN201510717763A CN105425289B CN 105425289 B CN105425289 B CN 105425289B CN 201510717763 A CN201510717763 A CN 201510717763A CN 105425289 B CN105425289 B CN 105425289B
Authority
CN
China
Prior art keywords
frequency
wave impedance
impedance
section wave
frequency section
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510717763.2A
Other languages
English (en)
Other versions
CN105425289A (zh
Inventor
岳跃龙
陈洪涛
李玉海
李冰玲
杨峰
曾凡平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Petroleum Corp
BGP Inc
Original Assignee
China National Petroleum Corp
BGP Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China National Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CN201510717763.2A priority Critical patent/CN105425289B/zh
Publication of CN105425289A publication Critical patent/CN105425289A/zh
Application granted granted Critical
Publication of CN105425289B publication Critical patent/CN105425289B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30—Analysis
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G01V2210/622—Velocity, density or impedance
    • G01V2210/6226—Impedance

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供了一种确定低频波阻抗的方法和装置，其中，该方法包括：获得测井低频波阻抗模型和地震数据低频波阻抗模型；分频段比较分析两种模型，得到补偿系数并进行补偿；将不同频率段的模型进行融合，得到融合后的低频波阻抗模型。在本发明实施例中，最终得到的具有测井和地震信息的低频波阻抗模型充分结合了钻井、地质及地震信息，因此具有更高的横向分辨率。该方法解决了现有技术中建立的低频波阻抗模型不准确的技术问题，达到了有效提高低频波阻抗模型准确度的技术效果，从而可以进一步提高波阻抗反演技术在预测岩性时的合理性。

Description

确定低频波阻抗的方法和装置
技术领域
本发明涉及石油勘探技术领域，特别涉及一种确定低频波阻抗的方法和装置。
背景技术
近年来，随着对岩性油气藏研究的不断深入，如何准确预测复杂岩性成为关注的焦点。
在预测复杂岩性的过程中，地震反演技术发挥了重要的作用。从地震反演所用的地震资料来分，地震反演可分为：叠前反演和叠后反演；从反演所利用地震的信息来分，地震反演可分为：地震波旅行时反演和地震波振幅反演；从反演的地质结果来分，可分为：构造反演、波阻抗反演、储层参数反演、地质统计反演等。
目前常用的地震反演软件主要以解释层位来控制测井曲线，并沿着固定的插值算法来计算波阻抗，从而建立低频波阻抗模型，其插值结果与测井距离成反比。这种建模方法较为简单，在地质条件简单的区域可以使用。但是在复杂的地质条件下，采用这种方式建立的低频波阻抗模型的准确性不高。
针对上述难以建立准确的低频波阻抗模型的技术问题，目前尚未提出有效的解决方案。
发明内容
本发明实施例提供了一种确定低频波阻抗的方法，以达到建立准确的低频波阻抗模型的目的。该方法包括：
本发明实施例提供了一种确定低频波阻抗的方法，包括：通过对待分析区域的测井曲线进行插值，得到测井低频波阻抗模型；通过对所述待分析区域的地震数据进行反演，得到地震数据低频波阻抗模型；将所述测井低频波阻抗模型，分为频率范围不同的N个第一频率段波阻抗，其中，N为正整数；将所述地震数据低频波阻抗模型，分为频率范围不同的N个第二频率段波阻抗，其中，所述第一频率段波阻抗与所述第二频率段波阻抗的频率范围一一对应；对各个一一对应的第一频率段波阻抗和第二频率段波阻抗的数据进行比较，根据比较结果计算出所述N个第二频率段波阻抗中各个第二频率段波阻抗的补偿系数；根据所述各个第二频率段波阻抗的补偿系数，对所述N个第二频率段波阻抗中各个第二频率段波阻抗进行补偿，得到N个补偿后的第二频率段波阻抗；对所述N个第一频率段波阻抗和所述N个补偿后的第二频率段波阻抗进行融合，得到融合后的低频波阻抗模型。
在一个实施例中，对所述N个第一频率段波阻抗和所述N个补偿后的第二频率段波阻抗进行融合，得到融合后的低频波阻抗模型，包括：对所述N个第一频率段波阻抗和所述N个补偿后的第二频率段波阻抗一对一进行比较；将频率范围相同的第一频率段波阻抗和第二频率段波阻抗中精确度高的频率段波阻抗作为该频率范围的融合频率段波阻抗，得到N个融合频率段波阻抗；对所述N个融合频率段波阻抗进行融合，得到融合后的低频波阻抗模型。
在一个实施例中，对所述N个融合频率段波阻抗进行融合，得到融合后的低频波阻抗模型，包括：对所述N个融合频率段波阻抗分别进行傅里叶变换；根据傅里叶变换后的结果，在频率域进行拼接，得到融合后的低频波阻抗模型。
在一个实施例中，N的取值为3或4。
在一个实施例中，当N取值为3时，所述N个第一频率段波阻抗和所述N个第二频率段波阻抗对应的3个频率范围为：0Hz到3Hz、3Hz到6Hz、6Hz到12Hz。
在一个实施例中，按照以下公式对所述N个第二频率段波阻抗中各个第二频率段波阻抗进行补偿：
AIi＝ai*AIsi
其中，AIi表示第i个补偿后的第二频率段波阻抗，ai表示第i个第二频率段波阻抗的补偿系数，AIsi表示第i个补偿前的第二频率段波阻抗，i的取值为1到N。
本发明实施例还提供了一种确定低频波阻抗的装置，包括：插值单元，用于通过对待分析区域的测井曲线进行插值，得到测井低频波阻抗模型；反演单元，用于通过对所述待分析区域的地震数据进行反演，得到地震数据低频波阻抗模型；第一分频单元，用于将所述测井低频波阻抗模型，分为频率范围不同的N个第一频率段波阻抗，其中，N为正整数；第二分频单元，用于将所述地震数据低频波阻抗模型，分为频率范围不同的N个第二频率段波阻抗，其中，所述第一频率段波阻抗与所述第二频率段波阻抗的频率范围一一对应；补偿系数计算单元，用于对各个一一对应的第一频率段波阻抗和第二频率段波阻抗的数据进行比较，根据比较结果计算出所述N个第二频率段波阻抗中各个第二频率段波阻抗的补偿系数；补偿单元，用于根据所述各个第二频率段波阻抗的补偿系数，对所述N个第二频率段波阻抗中各个第二频率段波阻抗进行补偿，得到N个补偿后的第二频率段波阻抗；融合单元，用于对所述N个第一频率段波阻抗和所述N个补偿后的第二频率段波阻抗进行融合，得到融合后的低频波阻抗模型。
在一个实施例中，所述融合单元包括：比较模块，用于对所述N个第一频率段波阻抗和所述N个补偿后的第二频率段波阻抗一对一进行比较；第一融合模块，用于将频率范围相同的第一频率段波阻抗和第二频率段波阻抗中精确度高的频率段波阻抗作为该频率范围的融合频率段波阻抗，得到N个融合频率段波阻抗；第二融合模块，用于对所述N个融合频率段波阻抗进行融合，得到融合后的低频波阻抗模型。
在一个实施例中，第二融合模块包括：傅里叶变换单元，用于对所述N个融合频率段波阻抗分别进行傅里叶变换；拼接单元，用于根据傅里叶变换后的结果，在频率域进行拼接，得到融合后的低频波阻抗模型。
在一个实施例中，N的取值为3或4。
在上述实施例中，在利用测井曲线进行插值得到测井低频波阻抗模型之后，引入了地震数据低频波阻抗模型以及补偿系数，通过将对应频率段的测井低频波阻抗模型和补偿后的地震数据低频波阻抗模型进行融合，得到最终的具有测井和地震信息的低频波阻抗模型。该模型充分结合了钻井、地质及地震信息，因此具有更高的横向分辨率。通过上述方式解决了现有技术中建立的低频波阻抗模型不准确的技术问题，达到了有效提高低频波阻抗模型准确度的技术效果，从而可以进一步提高波阻抗反演技术在预测岩性时的合理性。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中：
图1是根据本发明实施例的确定低频波阻抗方法的流程图；
图2是根据本发明实施例的采用测井曲线插值对碳酸盐岩储层分析得到的低频波阻抗模型示意图；
图3是根据本发明实施例的采用确定低频波阻抗的方法对碳酸盐岩储层分析得到的低频波阻抗模型示意图；
图4是根据本发明实施例的确定低频波阻抗装置的一种结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施方式和附图，对本发明做进一步详细说明。在此，本发明的示意性实施方式及其说明用于解释本发明，但并不作为对本发明的限定。
考虑到现有技术中难以建立准确的低频波阻抗模型的技术问题，引入了测井波阻抗曲线分析、低频地震数据频谱分析、低频波阻抗补偿处理、测井曲线与低频地震波阻抗融合计算等技术，充分结合钻井、地质及地震信息，以得到更为准确的低频波阻抗模型。具体的，在本例中，提供了一种确定低频波阻抗的方法，如图1所示，可以包括以下步骤：
步骤101：通过对待分析区域的测井曲线进行插值，得到测井低频波阻抗模型；
波阻抗指的是地震波在介质中传播时，作用于某个面积上的压力与单位时间内垂直通过此面积的质点流量(即面积乘以质点振动速度)之比，具有阻力的含义，称为波阻抗，其数值等于介质密度p与波速V的乘积。
波阻抗反演是储层预测的有效手段之一，特别地对于含油气储层，例如岩性油气藏，有很好的对应性。波阻抗模型的质量对地震反演的效果有着至关重要的作用。将该波阻抗模型与岩性数据进行对比，研究该波阻抗数据的幅值变化规律以及对应岩性的敏感程度，对待分析区域的数据采用已有的测井曲线进行传统的插值计算，即可以利用约束拟合算法将深度域的测井资料与时间域的地震层位资料联系起来，生成一对时间域和深度域相一致的、平滑的、闭合的实体模型，从而可以得到测井低频波阻抗模型AIwell。
步骤102：通过对所述待分析区域的地震数据进行反演，得到地震数据低频波阻抗模型；
在对待分析区域的地震数据进行处理之前，可以先对该地震数据进行宽频采集处理，使得地震数据的低频信息能够得到最大保留，然后再通过低通滤波得到低频地震数据，在得到低频地震数据后就可以直接通过现有的反演算法得到地震数据低频波阻抗模型。具体地，可以利用井旁道提取的地震子波与低频地震数据进行反褶积运算，然后再根据反射系数与波阻抗的关系计算得到地震数据低频波阻抗模型AIseismic。
步骤103：将所述测井低频波阻抗模型，分为频率范围不同的N个第一频率段波阻抗，其中，N为正整数；
具体的，N的取值可以为3或4。其中，当N的取值为3时，表示将测井低频波阻抗模型分为3个不同的频率范围，例如，如果总的频率范围为0Hz到12Hz，则可以按照如下分法将整个频率范围分为如下三个小的频率范围：0Hz到3Hz、3Hz到6Hz、6Hz到12Hz。
然而，值得注意的是，上述分成的频率段的个数，以及具体每个频率段的频率范围仅是一种示意性描述，并不构成对本发明的不当限定，例如，N可以按照需要选取除了3和4之外的值，具体的各频率段的范围也可以按照实际需要选取。同时，每个频率段之间可以采用等分的方式，也可以采用非等分的方式。
步骤104：将所述地震数据低频波阻抗模型，分为频率范围不同的N个第二频率段波阻抗，其中，所述第一频率段波阻抗与所述第二频率段波阻抗的频率范围一一对应；
即，在该步骤中，需要将地震数据低频波阻抗模型按照与测井低频波阻抗模型同样频率段进行划分，例如，如果测井低频波阻抗模型被分为3段，每个频率段波阻抗的频率范围为：0Hz到3Hz、3Hz到6Hz、6Hz到12Hz，那么地震数据低频波阻抗模型也应该分为3段，每个频率段波阻抗的频率范围为：0Hz到3Hz、3Hz到6Hz、6Hz到12Hz。
步骤105：对各个一一对应的第一频率段波阻抗和第二频率段波阻抗的数据进行比较，根据比较结果计算出所述N个第二频率段波阻抗中各个第二频率段波阻抗的补偿系数；
通常情况下，地震数据在采集和处理的过程中会在低频段有一定的损失，地震反演的相对波阻抗数据不能直接准确解释岩性，因此，需要对地震反演的结果进行低频补偿，进而得出补偿后的波阻抗模型。
由于地震数据的低频段波阻抗衰减较快，并且不同频率段所对应的波阻抗差不同。为了实现补偿，在本例中，分别选取测井低频波阻抗模型以及地震数据低频波阻抗模型的相同频率段的数据进行比较，可以计算得到每个频率段的地震波阻抗的补偿系数a。具体地，可以以测井低频波阻抗模型为基础值，地震数据低频波阻抗模型减去基础值可以得到差异值，这个差异值就是补偿系数。在计算补偿系数的过程中，可以预先对参与计算的数据进行检查，以确定数据中是否有奇异值，如果有奇异值，则先对这些奇异值进行平滑处理，然后再基于平滑处理后的数据计算补偿系数。
步骤106：根据各个第二频率段波阻抗的补偿系数，对所述N个第二频率段波阻抗中各个第二频率段波阻抗进行补偿，得到N个补偿后的第二频率段波阻抗；
在实际实施中，可以按照以下公式对N个不同频率段的地震数据低频波阻抗模型进行补偿：
AIi＝ai*AIsi (公式1)
其中，AIi表示地震数据低频波阻抗模型的第i个补偿后的频率段波阻抗，ai表示地震数据低频波阻抗模型的第i个频率段波阻抗的补偿系数，AIsi表示地震数据低频波阻抗模型的第i个补偿前的频率段波阻抗，i的取值为1到N。
即，可以按照上述公式1对每个不同频率段的地震数据低频波阻抗模型都进行补偿，以最终得到补偿后的地震数据低频波阻抗模型AI。
根据实际分析可知：在幅值上，补偿后的地震数据低频波阻抗模型与测井低频波阻抗模型相近；在空间上，补偿后的地震数据低频波阻抗模型具有更强的地震横向分辨率。
步骤107：对所述N个第一频率段波阻抗和所述N个补偿后的第二频率段波阻抗进行融合，得到融合后的低频波阻抗模型。
在实际执行的过程中，在得到补偿后的地震数据低频波阻抗模型后，要得到最终的用于地震反演岩性预测的低频波阻抗模型，还可以执行以下操作：
S1：对所述N个第一频率段波阻抗和所述N个补偿后的第二频率段波阻抗一对一进行比较；
仍旧以上述的N取3，各频率段波阻抗对应的频率范围为：0Hz到3Hz、3Hz到6Hz、6Hz到12Hz，为例进行说明。在该步骤中，相应地，就是将0Hz到3Hz所对应的补偿后的地震数据低频波阻抗模型和0Hz到3Hz所对应的测井低频波阻抗模型进行对比，将3Hz到6Hz所对应的补偿后的地震数据低频波阻抗模型和3Hz到6Hz所对应的测井低频波阻抗模型进行对比，将6Hz到12Hz所对应的补偿后的地震数据低频波阻抗模型和6Hz到12Hz所对应的测井低频波阻抗模型进行对比。
S2：将频率范围相同的第一频率段波阻抗和第二频率段波阻抗中精确度高的频率段波阻抗作为该频率范围的融合频率段波阻抗，得到N个融合频率段波阻抗；
例如，通过对比分析这3个频率段所得到的测井低频波阻抗模型以及补偿后的地震数据低频波阻抗模型发现：在0Hz到3Hz的极低频段时，采用测井低频波阻抗模型可靠性更好一些，在3Hz到6Hz、6Hz到12Hz的中高频段时，采用补偿后的地震数据低频波阻抗模型岩性预测的效果更好。
那么，就可以将0Hz到3Hz所对应的测井低频波阻抗模型，和3Hz到6Hz、6Hz到12Hz所对应的补偿后的地震数据低频波阻抗模型作为融合频率段波阻抗。
步骤3：将所述N个融合频率段波阻抗拼接成的低频波阻抗模型作为融合后的低频波阻抗模型。
具体的，在得到融合后的低频波阻抗模型之后，可以对上述3个融合频率段波阻抗分别进行傅里叶变换，再根据傅里叶变换后的结果，在频率域进行拼接，得到最终的用于地震反演岩性预测的低频波阻抗模型AIseismic+well。该模型具有准确的测井曲线垂向识别能力以及地震横向预测能力。
为了对上例中的确定低频波阻抗的方法的效果进行说明，在本例中给出了一组实验验证数据，该组数据仍旧是以上述的N取3，各频率段范围为：0Hz到3Hz、3Hz到6Hz、6Hz到12Hz，为例进行的。如图2所示为采用测井曲线插值对碳酸盐岩储层分析得到的低频波阻抗模型示意图，图3所示为采用上述的确定低频波阻抗的方法得到的低频波阻抗模型示意图。对图2和图3进行对比，可以看出，两幅图的幅值相近，然而，图3相对于图2的中高频段在空间上的地震横向分辨率更强。
由此可见，采用本发明实施例得到的低频波阻抗模型中包含有更加丰富的波阻抗信息，因而具有更加精确的地震反演岩性预测能力。
基于同一发明构思，本发明实施例中还提供了一种确定低频波阻抗的装置，如下面的实施例所述。由于确定低频波阻抗的装置解决问题的原理与确定低频波阻抗的方法相似，因此确定低频波阻抗的装置的实施可以参见确定低频波阻抗的方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。图4是本发明实施例的确定低频波阻抗的装置的一种结构框图，如图4所示，包括：插值单元401、反演单元402、第一分频单元403、第二分频单元404、补偿系数计算单元405、补偿单元406和融合单元407，下面对该结构进行说明。
插值单元401，用于通过对待分析区域的测井曲线进行插值，得到测井低频波阻抗模型；
反演单元402，用于通过对所述待分析区域的地震数据进行反演，得到地震数据低频波阻抗模型；
第一分频单元403，用于将所述测井低频波阻抗模型，分为频率范围不同的N个第一频率段波阻抗，其中，N为正整数；
第二分频单元404，用于将所述地震数据低频波阻抗模型，分为频率范围不同的N个第二频率段波阻抗，其中，所述第一频率段波阻抗与所述第二频率段波阻抗的频率范围一一对应；
补偿系数计算单元405，用于对各个一一对应的第一频率段波阻抗和第二频率段波阻抗的数据进行比较，根据比较结果计算出所述N个第二频率段波阻抗中各个第二频率段波阻抗的补偿系数；
补偿单元406，用于根据所述各个第二频率段波阻抗的补偿系数，对所述N个第二频率段波阻抗中各个第二频率段波阻抗进行补偿，得到N个补偿后的第二频率段波阻抗；
融合单元407，用于对所述N个第一频率段波阻抗和所述N个补偿后的第二频率段波阻抗进行融合，得到融合后的低频波阻抗模型。
在一个实施方式中，所述融合单元包括：比较模块，用于对所述N个第一频率段波阻抗和所述N个补偿后的第二频率段波阻抗一对一进行比较；第一融合模块，用于将频率范围相同的第一频率段波阻抗和第二频率段波阻抗中精确度高的频率段波阻抗作为该频率范围的融合频率段波阻抗，得到N个融合频率段波阻抗；第二融合模块，用于对所述N个融合频率段波阻抗进行融合，得到融合后的低频波阻抗模型。
在一个实施方式中，第二融合模块包括：傅里叶变换单元，用于对所述N个融合频率段波阻抗分别进行傅里叶变换；拼接单元，用于根据傅里叶变换后的结果，在频率域进行拼接，得到融合后的低频波阻抗模型。
在一个实施方式中，N的取值为3或4。
在一个实施方式中，当N取值为3时，所述N个第一频率段波阻抗和所述N个第二频率段波阻抗对应的3个频率范围为：0Hz到3Hz、3Hz到6Hz、6Hz到12Hz。
在一个实施方式中，按照以下公式对所述N个第二频率段波阻抗中各个第二频率段波阻抗进行补偿：
AIi＝ai*AIsi (公式2)
其中，AIi表示第i个补偿后的第二频率段波阻抗，ai表示第i个第二频率段波阻抗的补偿系数，AIsi表示第i个补偿前的第二频率段波阻抗，i的取值为1到N。
在另外一个实施例中，还提供了一种软件，该软件用于执行上述实施例及优选实施方式中描述的技术方案。
在另外一个实施例中，还提供了一种存储介质，该存储介质中存储有上述软件，该存储介质包括但不限于：光盘、软盘、硬盘、可擦写存储器等。
从以上的描述中，可以看出，本发明实施例实现了如下技术效果：在利用测井曲线进行插值得到测井低频波阻抗模型之后，引入了地震数据低频波阻抗曲线以及补偿系数，通过将对应频率段的测井低频波阻抗模型和补偿后的地震数据低频波阻抗模型进行融合，得到最终的具有测井和地震信息的低频波阻抗模型。该模型充分结合了钻井、地质及地震信息，因此具有更高的横向分辨率。通过上述方式解决了现有技术中建立的低频波阻抗模型不准确的技术问题，达到了有效提高低频波阻抗模型准确度的技术效果，从而可以进一步提高波阻抗反演技术在预测岩性时的合理性。
显然，本领域的技术人员应该明白，上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，并且在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明实施例不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (7)

1.一种确定低频波阻抗的方法，其特征在于，包括：
通过对待分析区域的测井曲线进行插值，得到测井低频波阻抗模型；
通过对所述待分析区域的地震数据进行反演，得到地震数据低频波阻抗模型；
将所述测井低频波阻抗模型，分为频率范围不同的N个第一频率段波阻抗，其中，N为正整数；
将所述地震数据低频波阻抗模型，分为频率范围不同的N个第二频率段波阻抗，其中，所述第一频率段波阻抗与所述第二频率段波阻抗的频率范围一一对应；
对各个一一对应的第一频率段波阻抗和第二频率段波阻抗的数据进行比较，根据比较结果计算出所述N个第二频率段波阻抗中各个第二频率段波阻抗的补偿系数；
根据所述各个第二频率段波阻抗的补偿系数，对所述N个第二频率段波阻抗中各个第二频率段波阻抗进行补偿，得到N个补偿后的第二频率段波阻抗；
对所述N个第一频率段波阻抗和所述N个补偿后的第二频率段波阻抗进行融合，得到融合后的低频波阻抗模型；
其中，对所述N个第一频率段波阻抗和所述N个补偿后的第二频率段波阻抗进行融合，得到融合后的低频波阻抗模型，包括：
对所述N个第一频率段波阻抗和所述N个补偿后的第二频率段波阻抗一对一进行比较；
将频率范围相同的第一频率段波阻抗和第二频率段波阻抗中精确度高的频率段波阻抗作为该频率范围的融合频率段波阻抗，得到N个融合频率段波阻抗；
对所述N个融合频率段波阻抗进行融合，得到融合后的低频波阻抗模型；
其中，按照以下公式对所述N个第二频率段波阻抗中各个第二频率段波阻抗进行补偿：
AIi＝ai*AIsi
其中，AIi表示第i个补偿后的第二频率段波阻抗，ai表示第i个第二频率段波阻抗的补偿系数，AIsi表示第i个补偿前的第二频率段波阻抗，i的取值为1到N。
2.如权利要求1所述的方法，其特征在于，对所述N个融合频率段波阻抗进行融合，得到融合后的低频波阻抗模型，包括：
对所述N个融合频率段波阻抗分别进行傅里叶变换；
根据傅里叶变换后的结果，在频率域进行拼接，得到融合后的低频波阻抗模型。
3.如权利要求1至2中任一项所述的方法，其特征在于，N的取值为3或4。
4.如权利要求3所述的方法，其特征在于，当N取值为3时，所述N个第一频率段波阻抗和所述N个第二频率段波阻抗对应的3个频率范围为：0Hz到3Hz、3Hz到6Hz、6Hz到12Hz。
5.一种确定低频波阻抗的装置，其特征在于，包括：
插值单元，用于通过对待分析区域的测井曲线进行插值，得到测井低频波阻抗模型；
反演单元，用于通过对所述待分析区域的地震数据进行反演，得到地震数据低频波阻抗模型；
第一分频单元，用于将所述测井低频波阻抗模型，分为频率范围不同的N个第一频率段波阻抗，其中，N为正整数；
第二分频单元，用于将所述地震数据低频波阻抗模型，分为频率范围不同的N个第二频率段波阻抗，其中，所述第一频率段波阻抗与所述第二频率段波阻抗的频率范围一一对应；
补偿系数计算单元，用于对各个一一对应的第一频率段波阻抗和第二频率段波阻抗的数据进行比较，根据比较结果计算出所述N个第二频率段波阻抗中各个第二频率段波阻抗的补偿系数；
补偿单元，用于根据所述各个第二频率段波阻抗的补偿系数，对所述N个第二频率段波阻抗中各个第二频率段波阻抗进行补偿，得到N个补偿后的第二频率段波阻抗；
融合单元，用于对所述N个第一频率段波阻抗和所述N个补偿后的第二频率段波阻抗进行融合，得到融合后的低频波阻抗模型；
其中，所述融合单元包括：
比较模块，用于对所述N个第一频率段波阻抗和所述N个补偿后的第二频率段波阻抗一对一进行比较；
第一融合模块，用于将频率范围相同的第一频率段波阻抗和第二频率段波阻抗中精确度高的频率段波阻抗作为该频率范围的融合频率段波阻抗，得到N个融合频率段波阻抗；
第二融合模块，用于对所述N个融合频率段波阻抗进行融合，得到融合后的低频波阻抗模型；
其中，所述补偿单元按照以下公式对所述N个第二频率段波阻抗中各个第二频率段波阻抗进行补偿：
AIi＝ai*AIsi
其中，AIi表示第i个补偿后的第二频率段波阻抗，ai表示第i个第二频率段波阻抗的补偿系数，AIsi表示第i个补偿前的第二频率段波阻抗，i的取值为1到N。
6.如权利要求5所述的装置，其特征在于，第二融合模块包括：
傅里叶变换单元，用于对所述N个融合频率段波阻抗分别进行傅里叶变换；
拼接单元，用于根据傅里叶变换后的结果，在频率域进行拼接，得到融合后的低频波阻抗模型。
7.如权利要求5至6中任一项所述的装置，其特征在于，N的取值为3或4。
CN201510717763.2A 2015-10-29 2015-10-29 确定低频波阻抗的方法和装置 Active CN105425289B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510717763.2A CN105425289B (zh) 2015-10-29 2015-10-29 确定低频波阻抗的方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510717763.2A CN105425289B (zh) 2015-10-29 2015-10-29 确定低频波阻抗的方法和装置

Publications (2)

Publication Number Publication Date
CN105425289A CN105425289A (zh) 2016-03-23
CN105425289B true CN105425289B (zh) 2018-08-14

Family

ID=55503613

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510717763.2A Active CN105425289B (zh) 2015-10-29 2015-10-29 确定低频波阻抗的方法和装置

Country Status (1)

Country Link
CN (1) CN105425289B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106324675B (zh) * 2016-10-09 2018-09-07 中国石油大学(华东) 一种宽频地震波阻抗低频信息预测方法及系统
CN106772579B (zh) * 2016-12-07 2018-12-14 中国矿业大学(北京) 一种薄煤层中地震叠前反演方法和装置
CN107179547A (zh) * 2017-06-06 2017-09-19 中海石油(中国)有限公司 一种地震波阻抗反演低频模型建立方法
CN109283574B (zh) * 2017-07-19 2021-03-19 中国石油化工股份有限公司 低频模型构建方法及计算机可读存储介质
CN108680955B (zh) * 2017-12-12 2019-09-17 中国地质大学（武汉） 一种岩性识别方法及识别系统
CN108919344B (zh) * 2018-03-30 2020-04-21 北京诺克斯达石油科技有限公司 适用于层状介质的分频构形反演方法
CN110954947B (zh) * 2018-09-26 2021-09-17 中国石油化工股份有限公司 时域线性低频融合方法及系统
CN109444958B (zh) * 2018-11-01 2020-05-26 科吉思石油技术咨询(北京)有限公司 致密储层孔隙度定量预测的方法
CN112147700A (zh) * 2019-06-28 2020-12-29 中国石油化工股份有限公司 速度异常区的低频模型构建方法及系统
CN113917529A (zh) * 2020-07-07 2022-01-11 中国石油化工股份有限公司 波阻抗反演方法、装置、存储介质及电子设备

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4680741A (en) * 1981-10-13 1987-07-14 Geosource, Inc. Method and apparatus for seismic exploration using non-linear sweeps
FR2779533B1 (fr) * 1998-06-09 2000-08-25 Total Sa Dispositif d'acquisition sismique a haute resolution
CN101206264A (zh) * 2007-11-08 2008-06-25 符力耘 高分辨率非线性地震波阻抗反演方法
CN103293551B (zh) * 2013-05-24 2015-11-04 中国石油天然气集团公司 一种基于模型约束的阻抗反演方法及系统
CN103527184B (zh) * 2013-10-28 2016-04-20 北京大学 一种白云岩储层的预测方法和系统
CN104570125B (zh) * 2014-09-26 2016-09-28 中国石油天然气股份有限公司 一种利用井数据提高成像速度模型精度的方法

Also Published As

Publication number Publication date
CN105425289A (zh) 2016-03-23

Similar Documents

Publication Publication Date Title
CN105425289B (zh) 确定低频波阻抗的方法和装置
US10705238B2 (en) Method and apparatus for processing seismic data
CN109425896B (zh) 白云岩油气储层分布预测方法及装置
CN106033124B (zh) 一种基于随机最优化的多震源粘声最小二乘逆时偏移方法
CN107272058A (zh) 成像方法、成像装置以及计算机存储介质
CN111123354A (zh) 基于频变反射振幅衰减预测致密气层的方法及设备
CN113341455B (zh) 一种粘滞各向异性介质地震波数值模拟方法、装置及设备
CN103630932A (zh) 一种地震数据分形保幅方法
CN110542928A (zh) 基于vti各向异性传播矩阵的地震响应模拟方法
CN103744114B (zh) 基于零偏垂直地震剖面数据估计品质因子的方法和装置
CN105954803B (zh) 叠后地震反演方法及装置
CN105005075B (zh) 基于地震频率信息的多波匹配方法
CN111060961B (zh) 基于多信息约束反演的品质因子确定方法、装置及系统
CN104391324A (zh) 依赖频率的avo反演前的地震道集动校拉伸校正预处理技术
CN106574980A (zh) 用于地下地质体的岩石性质估计的系统和方法
CN103675904A (zh) 一种井震匹配解释性目标处理方法及装置
CN110297273A (zh) 储层的叠前反演方法及装置
CN113552624B (zh) 孔隙度预测方法及装置
CN106772599B (zh) 一种计算地层横波速度的方法及装置
CN106547019A (zh) 一种确定地层品质因子的方法
CN109471162A (zh) 层间多次波处理方法、系统、电子设备及可读介质
CN107678065B (zh) 提高地震分辨率的保构造井控空间反褶积方法和装置
Chai et al. Q-compensated acoustic impedance inversion of attenuated seismic data: Numerical and field-data experiments
CN106054245B (zh) 地震资料高分辨率处理定量评价方法
CN110988991B (zh) 一种弹性参数反演方法、装置及系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant