CN105510965B - 流体识别方法和装置 - Google Patents
流体识别方法和装置 Download PDFInfo
- Publication number
- CN105510965B CN105510965B CN201511006094.4A CN201511006094A CN105510965B CN 105510965 B CN105510965 B CN 105510965B CN 201511006094 A CN201511006094 A CN 201511006094A CN 105510965 B CN105510965 B CN 105510965B
- Authority
- CN
- China
- Prior art keywords
- frequency
- mrow
- carried out
- data
- obtains
- 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
Links
- 239000012530 fluid Substances 0.000 title claims abstract description 74
- 238000000034 method Methods 0.000 title claims abstract description 21
- 239000006185 dispersion Substances 0.000 claims abstract description 62
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 27
- 238000001615 p wave Methods 0.000 claims abstract description 27
- 238000011160 research Methods 0.000 claims abstract description 26
- 238000012545 processing Methods 0.000 claims abstract description 18
- 230000000295 complement effect Effects 0.000 claims abstract description 14
- 238000010030 laminating Methods 0.000 claims description 7
- 230000014509 gene expression Effects 0.000 claims description 6
- 239000003921 oil Substances 0.000 description 10
- 238000005516 engineering process Methods 0.000 description 9
- 230000002159 abnormal effect Effects 0.000 description 6
- 238000001514 detection method Methods 0.000 description 6
- 230000000644 propagated effect Effects 0.000 description 5
- 239000011435 rock Substances 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 229920006395 saturated elastomer Polymers 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 230000002547 anomalous effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 230000002062 proliferating effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
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. for interpretation or for event detection
- 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
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- 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
-
- 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/624—Reservoir parameters
-
- 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/64—Geostructures, e.g. in 3D data cubes
- G01V2210/647—Gas hydrates
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
本发明提供了一种流体识别方法和装置,其中,该方法包括:获取研究区的叠前角道集地震数据体;对所述叠前角道集地震数据进行分角度叠加,得到分角度叠加数据体;采用所述分角度叠加数据体进行叠前AVO反演,得到纵波阻抗体和密度体;对所述纵波阻抗体进行时频分解处理,得到分频波阻抗体;结合所述密度体,对所述分频波阻抗体进行反演运算,得到所述研究区的储层地震频散属性;根据所述储层地震频散属性进行流体识别。本发明解决了现有技术中无法准确有效地对流体进行识别的技术问题,达到了准确识别流体的目的。
Description
技术领域
本发明涉及地质勘探技术领域,特别涉及一种流体识别方法和装置。
背景技术
随着勘探开发程度的不断提高,勘探目标逐渐由最初的构造油气藏转变为岩性油气藏和隐蔽性油气藏,勘探策略也由最开始的寻找构造圈闭转变为寻找岩性异常体或者直接进行流体检测。近年来,基于地震资料反演储层弹性参数的方法已被广泛应用于储层预测当中,并且储层预测的目的也不再局限于构造、岩性圈闭识别,而是要对储层含流体性(包括流体类型以及流体饱和度)进行检测,这就对流体识别技术提出了更高的要求。
上世纪70年代,基于地震资料振幅信息的流体识别技术作为首个直接进行地下含油气性检测技术被广泛应用在实际的勘探开发中。随后,依据地震波在地下介质中传播时会在储层油气水分界面处产生沿水平方向振幅强度稳定的反射波同相轴的平点勘探方式也开始被使用。80年代提出的AVO技术,主要是利用地震波反射振幅随偏移距或者入射角的变化规律判断含气砂岩储层,使得利用AVO分析技术检测流体成为可能。流体因子是利用纵横波速度的加权组合来预测储层流体,随后将流体因子引入到反射系数域,计算纵横波反射系数的加权组合。1999年之后,对流体因子的研究逐渐从反射系数域转变到波阻抗域,其中λρ、μρ属性是阻抗域最主要且应用效果最显著的流体因子,λρ_μρ属性交会图要比常规的纵横波阻抗交会图对储层流体有更好的识别能力。
随后,Russell等对流体因子进行了进一步的改进得到了适用性更强的流体因子。近年来,对于流体因子的研究不断深入,应用领域逐渐从叠后发展到叠前、从各向同性介质发展到各向异性介质、从时间域发展到频率域。
然而,如何准确有效地对流体进行识别,目前尚未提出有效的解决方案。
发明内容
本发明实施例提供了一种流体识别方法,以达到准确识别流体的目的,该方法包括:
获取研究区的叠前角道集地震数据体;
对所述叠前角道集地震数据进行分角度叠加,得到分角度叠加数据体;
采用所述分角度叠加数据体进行叠前AVO反演,得到纵波阻抗体和密度体;
对所述纵波阻抗体进行时频分解处理,得到分频波阻抗体;
结合所述密度体,对所述分频波阻抗体进行反演运算,得到所述研究区的储层地震频散属性;
根据所述储层地震频散属性进行流体识别。
在一个实施方式中,按照以下公式结合所述密度体,对所述分频波阻抗体进行反演运算,得到所述研究区的储层地震频散属性:
其中,z(t,f)表示特定频率时的波阻抗体,z(t,f0)表示参考频率时的波阻抗体,ρ表示所述密度体,D表示反演得到的速度频散值。
在一个实施方式中,所述参考频率为所述叠前角道集地震数据体的主频率。
在一个实施方式中,对所述叠前角道集地震数据进行分角度叠加,得到分角度叠加数据体,包括:
确定分角度叠加的多个角度区间;
根据所述多个角度区间对所述叠前角道集地震数据进行分角度叠加,得到所述多个角度区间中各个角度区间对应的部分叠加数据体。
在一个实施方式中,对所述纵波阻抗体进行时频分解处理,得到分频波阻抗体,包括:
通过S变换对所述纵波阻抗体进行时频分解处理,得到分频波阻抗体。
本发明实施例还提供了一种流体识别装置,以达到准确识别流体的目的,该装置包括:
获取模块,用于获取研究区的叠前角道集地震数据体;
叠加模块,用于对所述叠前角道集地震数据进行分角度叠加,得到分角度叠加数据体;
反演模块,用于采用所述分角度叠加数据体进行叠前AVO反演,得到纵波阻抗体和密度体;
时频分解模块,用于对所述纵波阻抗体进行时频分解处理,得到分频波阻抗体;
运算模块,用于结合所述密度体,对所述分频波阻抗体进行反演运算,得到所述研究区的储层地震频散属性;
识别模块,用于根据所述储层地震频散属性进行流体识别。
在一个实施方式中,所述运算模块具体用于按照以下公式结合所述密度体,对所述分频波阻抗体进行反演运算,得到所述研究区的储层地震频散属性:
其中,z(t,f)表示特定频率时的波阻抗体,z(t,f0)表示参考频率时的波阻抗体,ρ表示所述密度体,D表示反演得到的速度频散值。
在一个实施方式中,所述参考频率为所述叠前角道集地震数据体的主频率。
在一个实施方式中,所述叠加模块包括:
角度区间确定单元,用于确定分角度叠加的多个角度区间;
叠加单元,用于根据所述多个角度区间对所述叠前角道集地震数据进行分角度叠加,得到所述多个角度区间中各个角度区间对应的部分叠加数据体。
在一个实施方式中,所述时频分解模块具体用于通过S变换对所述纵波阻抗体进行时频分解处理,得到分频波阻抗体。
在本发明实施例中,基于研究区叠前AVO反演得到的纵波阻抗数据体与密度体,利用时频分析技术对波阻抗数据体进行时频分解得到不同瞬时频率的阻抗体,然后结合地震波阻抗依赖分频反演公式提取地震频散数据体,最后根据频散属性实现流体识别,从而实现对储层流体分布规律的有效预测,提高了流体识别的准确性。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,并不构成对本发明的限定。在附图中:
图1是根据本发明实施例的流体识别方法流程图;
图2是根据本发明实施例的储层饱和不同流体时纵波速度值随频率变化关系图;
图3是根据本发明实施例的储层饱和不同流体时纵波频散值随频率变化关系图;
图4是根据本发明实施例的利用地震波阻抗依赖分频反演得到的地震频散属性剖面示意图;
图5是根据本发明实施例的利用地震波阻抗依赖分频反演得到的地震频散属性沿目的层的切片示意图;
图6是根据本发明实施例的流体识别装置的结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施方式和附图,对本发明做进一步详细说明。在此,本发明的示意性实施方式及其说明用于解释本发明,但并不作为对本发明的限定。
鉴于地震波在地下含流体介质中传播时会在不同频带范围和不同尺度上发生反射波能量衰减与速度频散,而这种现象主要是由储层孔隙流体的存在导致的。频散属性作为常规储层弹性信息的一个延伸和补充,从另一个角度展示了地震数据中被主频分量所掩盖的高频和低频信息,对于岩性及隐蔽油气藏的识别尤其是对流体检测具有独特的优势。与频率有关的速度变化主要受到储层孔隙流体的影响,并且流体的流动性(即岩石渗透率与流体粘滞性的比值)是影响地震波速度频变特征的主要因素。
发明人考虑到可以基于波阻抗依赖频散属性反演进行流体识别,具体的,可以综合利用波阻抗反演与时频分解技术来提取地震波频散属性,进而利用得到的地震波频散属性进行储层流体识别,如图1所示,该流体识别方法包括以下步骤:
步骤101:获取研究区的叠前角道集地震数据体;
步骤102:对所述叠前角道集地震数据进行分角度叠加,得到分角度叠加数据体;
在实现的时候,可以形成多个叠加道集,即形成多个分角度叠加数据体。例如,可以确定分角度叠加的多个角度区间,根据所述多个角度区间对叠前角道集地震数据进行分角度叠加,以得到多个角度区间中各个角度区间对应的部分叠加数据体。假设一共得到三个角度区间的部分叠加数据体,可以是第一角度区间、第二角度区间和第三角度区间,第一角度区间可以是5度至13度(近分角度叠加数据体),第二角度区间可以是14度至24度(中分角度叠加数据体),第三角度区间可以是25度至33度(远分角度叠加数据体)。一般情况下,不同工区和不同数据体的角度区间范围是不一样的,具体的角度大小可以根据实际情况确定。
步骤103:采用所述分角度叠加数据体进行叠前AVO反演,得到纵波阻抗体和密度体;
其中,AVO(Amplitude Versus Offset,振幅随偏移距的变化)主要是用于研究地震反射振幅随炮点与接收器之间的距离即炮检距(或入射角)的变化特征来探讨反射系数响应随炮检距(或入射角)的变化,进而确定出反射界面上覆、下伏介质的岩性特征及物性参数。AVO反演属于叠前反演,一般是利用Zoeppritz方程的近似公式,从叠前地震数据中估算岩石的弹性参数,用于岩性分析和油气预测。
在本例中是对分角度叠加数据体进行叠前AVO反演,以得到后续需要用到的纵波阻抗体和密度体。
步骤104:对所述纵波阻抗体进行时频分解处理,得到分频波阻抗体;
在具体实现的时候,考虑到S变换具有高时频分辨率,因此可以采用S变换对纵波阻抗体进行时频分解处理,以得到分频波阻抗体。
步骤105:结合所述密度体,对所述分频波阻抗体进行反演运算,得到所述研究区的储层地震频散属性;
下面对反演运算的公式进行推导,以便更好地说明反演运算的原理:
波阻抗可以表示:Z(t)=ρV(t) (公式1)
考虑到速度是频变的,根据公式1可知波阻抗也是频变的。假定密度ρ是不随频率变化的,则时间-频率域波阻抗可以表达为:
Z(t,f)=ρV(t,f) (公式2)
上式中只有波阻抗与速度是频率的函数,对其在参考频率处进行一阶泰勒展开,可以得到:
(公式3)
在参考频率处有:
Z(t,f0)=ρV(t,f0) (公式4)
将公式4代入上述公式3,可以得到:
(公式5)
从而可以推导出表征地震波频散属性的公式:
(公式6)
其中,z(t,f)表示特定频率时的波阻抗体,z(t,f0)表示参考频率时的波阻抗体,ρ表示所述密度体,D表示反演得到的速度频散值,可以表示为:
其中,上述的参考频率可以近似等于叠前角道集地震数据体的主频率。
步骤106:根据所述储层地震频散属性进行流体识别。
在上述步骤106中可以按照以下规则进行流体识别:当储层含不同流体时会表现出不同程度的频散异常,一般情况下,储层含气时的频散幅度要强于含油时,同时含油时的频散异常幅度要强于含水储层的频散异常幅度,根据反演得到的频散属性体,可以确定出高频散异常区带,从而达到流体识别的目的,一般原则是:频散异常幅度越强,则含油气的概率越大。
在上例中,基于研究区叠前AVO反演得到的纵波阻抗数据体与密度体,利用时频分析技术对波阻抗数据体进行时频分解得到不同瞬时频率的阻抗体,然后结合地震波阻抗依赖分频反演公式提取地震频散数据体,最后根据频散属性实现流体识别,从而实现对储层流体分布规律的有效预测,提高了流体识别的准确性。
在本发明实施例中还提供了一个具体的实施例对上述流体识别方法进行说明,然而值得注意的是,该具体实施例仅是为了更好地说明本发明,并不构成对本发明的不当限定。
具体的,该流体识别方法可以包括以下步骤:
S1:获取研究区叠前角道集地震数据体;
S2:对得到的角道集数据体进行分角度叠加,形成近、中、远分角度叠加数据;
S3:利用分角度叠加数据进行常规叠前AVO反演,得到纵波阻抗体以及密度体;
S4:利用广义S变换对得到的纵波阻抗数据体进行时频分解处理,得到分频波阻抗体;
S5:结合AVO反演得到的密度体对分频处理后的纵波阻抗体进行反演运算,得到研究区储层地震频散属性;
根据岩石物理理论的研究可知:当地震波在地下介质中传播时,会发生速度频散与能量衰减,使得含油气储层的地震响应具有明显的频变特征。然而,导致地震响应出现频散及衰减特征的主要因素是地下介质中孔隙流体的流动,其理论基础是岩石物理的速度频变特征。如图2所示,表明了地震波速度频变特征与地层含流体性存在一定的关系,具体的,当地层饱和气时所引起的速度频散程度要强于地层饱和油时所引起的速度频散程度,而含水时速度频变特征最弱。在地震频带范围内由于不同类型流体的流动性存在较大差异导致地震波在饱和不同流体的多孔介质中传播时表现出不同的速度频变特征。
如图3所示,为地震波在饱和不同流体地层中传播时的频散特征曲线。在地震频带范围内,饱和不同流体的地震波频散值存在明显差异:含气时最强,油次之,水最弱。
S6:利用地震频散属性进行流体识别。
选取一针对陆相沉积盆地含油气储层的三维地震数据进行流体识别研究,该地区储层特征以缝洞型为主,为典型的碳酸盐岩溶蚀孔洞储层。
根据上述流体识别方法,首先,对角度域数据体进行分角度叠加,得到三个角度(5-13度、14-24度、25-33度)的部分叠加数据体;然后,利用常规的叠前AVO反演得到纵波阻抗和密度数据体,接着对纵波阻抗数据进行时频分解处理,选取的频率为10Hz、20Hz、30Hz、40Hz、50Hz和60Hz,其中,地震数据主频接近30Hz,因此,参考频率选为30Hz。
然后,根据反演得到的密度体以及分频率阻抗体结合反演公式反演得到纵波频散属性体。如图4所示,为反演得到的频散属性剖面,在图4中,白色部分为频散属性高值区域,黑色部分为频散属性低值区域。测井与试井数据显示M1、M2和M3井均为高产油气井,而M4井为失利井仅见弱油气显示,因此,纵波频散剖面很好的吻合了井的实际生产情况,验证了频散属性对于流体检测的有效性。如图5所示为反演得到的频散属性沿目的层的切片显示图,其中,黑色表示频散属性异常高值,白色表示频散属性异常低值。由图5可以看出,在井点位置上反演的频散属性与测井、试井结果高度吻合,从而说明本例所提供的流体识别方法可以有效检测储层的含油气性,通过该流体识别方法进行流体识别降低了油气勘探的风险,提高了储层预测的精度。
基于同一发明构思,本发明实施例中还提供了一种流体识别装置,如下面的实施例所述。由于流体识别装置解决问题的原理与流体识别方法相似,因此流体识别装置的实施可以参见流体识别方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。图6是本发明实施例的流体识别装置的一种结构框图,如图6所示,可以包括:获取模块601、叠加模块602、反演模块603、时频分解模块604、运算模块605和识别模块606,下面对该结构进行说明。
获取模块601,用于获取研究区的叠前角道集地震数据体;
叠加模块602,用于对所述叠前角道集地震数据进行分角度叠加,得到分角度叠加数据体;
反演模块603,用于采用所述分角度叠加数据体进行叠前AVO反演,得到纵波阻抗体和密度体;
时频分解模块604,用于对所述纵波阻抗体进行时频分解处理,得到分频波阻抗体;
运算模块605,用于结合所述密度体,对所述分频波阻抗体进行反演运算,得到所述研究区的储层地震频散属性;
识别模块606,用于根据所述储层地震频散属性进行流体识别。
在一个实施方式中,运算模块605具体可以用于按照以下公式结合所述密度体,对所述分频波阻抗体进行反演运算,得到所述研究区的储层地震频散属性:
其中,z(t,f)表示特定频率时的波阻抗体,z(t,f0)表示参考频率时的波阻抗体,ρ表示所述密度体,D表示反演得到的速度频散值。
在一个实施方式中,参考频率可以是叠前角道集地震数据体的主频率。
在一个实施方式中,叠加模块602可以包括:角度区间确定单元,用于确定分角度叠加的多个角度区间;叠加单元,用于根据所述多个角度区间对所述叠前角道集地震数据进行分角度叠加,得到所述多个角度区间中各个角度区间对应的部分叠加数据体。
在一个实施方式中,时频分解模块604具体可以用于通过S变换对所述纵波阻抗体进行时频分解处理,得到分频波阻抗体。
从以上的描述中,可以看出,本发明实施例实现了如下技术效果:基于研究区叠前AVO反演得到的纵波阻抗数据体与密度体,利用时频分析技术对波阻抗数据体进行时频分解得到不同瞬时频率的阻抗体,然后结合地震波阻抗依赖分频反演公式提取地震频散数据体,最后根据频散属性实现流体识别,从而实现对储层流体分布规律的有效预测,提高了流体识别的准确性。
显然,本领域的技术人员应该明白,上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,并且在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明实施例不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种流体识别方法,其特征在于,包括:
获取研究区的叠前角道集地震数据体;
对所述叠前角道集地震数据进行分角度叠加,得到分角度叠加数据体;
采用所述分角度叠加数据体进行叠前AVO反演,得到纵波阻抗体和密度体;
对所述纵波阻抗体进行时频分解处理,得到分频波阻抗体;
结合所述密度体,对所述分频波阻抗体进行反演运算,得到所述研究区的储层地震频散属性;
根据所述储层地震频散属性进行流体识别;
其中,按照以下公式结合所述密度体,对所述分频波阻抗体进行反演运算,得到所述研究区的储层地震频散属性:
<mrow>
<mi>D</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>z</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>z</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&rho;</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>-</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
其中,z(t,f)表示特定频率时的波阻抗体,z(t,f0)表示参考频率时的波阻抗体,ρ表示所述密度体,D表示反演得到的速度频散值,f表示特定频率,f0表示参考频率,ρ(f-f0)表示ρ与(f-f0)相乘;
其中,所述参考频率为所述叠前角道集地震数据体的主频率。
2.如权利要求1所述的方法,其特征在于,对所述叠前角道集地震数据进行分角度叠加,得到分角度叠加数据体,包括:
确定分角度叠加的多个角度区间;
根据所述多个角度区间对所述叠前角道集地震数据进行分角度叠加,得到所述多个角度区间中各个角度区间对应的部分叠加数据体。
3.如权利要求1或2所述的方法,其特征在于,对所述纵波阻抗体进行时频分解处理,得到分频波阻抗体,包括:
通过S变换对所述纵波阻抗体进行时频分解处理,得到分频波阻抗体。
4.一种流体识别装置,其特征在于,包括:
获取模块,用于获取研究区的叠前角道集地震数据体;
叠加模块,用于对所述叠前角道集地震数据进行分角度叠加,得到分角度叠加数据体;
反演模块,用于采用所述分角度叠加数据体进行叠前AVO反演,得到纵波阻抗体和密度体;
时频分解模块,用于对所述纵波阻抗体进行时频分解处理,得到分频波阻抗体;
运算模块,用于结合所述密度体,对所述分频波阻抗体进行反演运算,得到所述研究区的储层地震频散属性;
识别模块,用于根据所述储层地震频散属性进行流体识别;
其中,所述运算模块具体用于按照以下公式结合所述密度体,对所述分频波阻抗体进行反演运算,得到所述研究区的储层地震频散属性:
<mrow>
<mi>D</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>z</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>z</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&rho;</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>-</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
其中,z(t,f)表示特定频率时的波阻抗体,z(t,f0)表示参考频率时的波阻抗体,ρ表示所述密度体,D表示反演得到的速度频散值,f表示特定频率,f0表示参考频率,ρ(f-f0)表示ρ与(f-f0)相乘;
其中,所述参考频率为所述叠前角道集地震数据体的主频率。
5.如权利要求4所述的装置,其特征在于,所述叠加模块包括:
角度区间确定单元,用于确定分角度叠加的多个角度区间;
叠加单元,用于根据所述多个角度区间对所述叠前角道集地震数据进行分角度叠加,得到所述多个角度区间中各个角度区间对应的部分叠加数据体。
6.如权利要求4或5所述的装置,其特征在于,所述时频分解模块具体用于通过S变换对所述纵波阻抗体进行时频分解处理,得到分频波阻抗体。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511006094.4A CN105510965B (zh) | 2015-12-29 | 2015-12-29 | 流体识别方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511006094.4A CN105510965B (zh) | 2015-12-29 | 2015-12-29 | 流体识别方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105510965A CN105510965A (zh) | 2016-04-20 |
CN105510965B true CN105510965B (zh) | 2018-01-05 |
Family
ID=55719073
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201511006094.4A Active CN105510965B (zh) | 2015-12-29 | 2015-12-29 | 流体识别方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105510965B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105938203B (zh) * | 2016-06-24 | 2018-07-10 | 中国石油天然气股份有限公司 | 一种储层特性的检测方法及装置 |
CN112684503B (zh) * | 2020-11-16 | 2022-08-05 | 中国石油天然气股份有限公司 | 一种基于叠后地震属性的储层流体识别方法及装置 |
CN112748467B (zh) * | 2020-12-24 | 2022-10-14 | 中国海洋石油集团有限公司 | 利用地震资料确定Russell流体因子的方法及系统 |
CN116931080B (zh) * | 2022-04-12 | 2024-06-11 | 中国石油化工股份有限公司 | 一种基于叠前频变的流体性质检测方法 |
CN116699684A (zh) * | 2023-06-09 | 2023-09-05 | 中海石油(中国)有限公司深圳分公司 | 面向储层预测的地震数据拼接方法、装置、设备及介质 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102478667A (zh) * | 2010-11-30 | 2012-05-30 | 中国石油天然气集团公司 | 频散avo反演方法 |
CN103235339A (zh) * | 2013-04-09 | 2013-08-07 | 中国石油大学(北京) | 一种时频分解地震流体识别方法 |
CN103364834A (zh) * | 2013-07-29 | 2013-10-23 | 成都晶石石油科技有限公司 | 一种利用叠前地震频散分析预测储层渗透率的方法 |
CN103713321A (zh) * | 2014-01-08 | 2014-04-09 | 孙赞东 | 基于纵波频率依赖方位avo的裂缝流体类型识别方法 |
CN103792574A (zh) * | 2012-10-30 | 2014-05-14 | 中国石油化工股份有限公司 | 一种储层频变含气性检测方法 |
CN103984010A (zh) * | 2014-04-16 | 2014-08-13 | 孙赞东 | 一种基于三项频率依赖avo反演的流体识别方法 |
CN104237936A (zh) * | 2013-06-06 | 2014-12-24 | 中国石油化工股份有限公司 | 一种油气检测的频变反演方法 |
CN104297791A (zh) * | 2014-09-25 | 2015-01-21 | 中国石油天然气股份有限公司 | 一种基于地震优势频率的反演方法及系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9482771B2 (en) * | 2012-11-05 | 2016-11-01 | Fugro Marine Geoservices, Inc. | Method of indicating the presence of gas hydrate and shallow gas in deepwater environment |
-
2015
- 2015-12-29 CN CN201511006094.4A patent/CN105510965B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102478667A (zh) * | 2010-11-30 | 2012-05-30 | 中国石油天然气集团公司 | 频散avo反演方法 |
CN103792574A (zh) * | 2012-10-30 | 2014-05-14 | 中国石油化工股份有限公司 | 一种储层频变含气性检测方法 |
CN103235339A (zh) * | 2013-04-09 | 2013-08-07 | 中国石油大学(北京) | 一种时频分解地震流体识别方法 |
CN104237936A (zh) * | 2013-06-06 | 2014-12-24 | 中国石油化工股份有限公司 | 一种油气检测的频变反演方法 |
CN103364834A (zh) * | 2013-07-29 | 2013-10-23 | 成都晶石石油科技有限公司 | 一种利用叠前地震频散分析预测储层渗透率的方法 |
CN103713321A (zh) * | 2014-01-08 | 2014-04-09 | 孙赞东 | 基于纵波频率依赖方位avo的裂缝流体类型识别方法 |
CN103984010A (zh) * | 2014-04-16 | 2014-08-13 | 孙赞东 | 一种基于三项频率依赖avo反演的流体识别方法 |
CN104297791A (zh) * | 2014-09-25 | 2015-01-21 | 中国石油天然气股份有限公司 | 一种基于地震优势频率的反演方法及系统 |
Non-Patent Citations (1)
Title |
---|
"基于广义S变换的分频AVO研究";杨鑫;《中国优秀硕士学位论文全文数据库 基础科学辑》;20140515(第2014(05)期);第1-69页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105510965A (zh) | 2016-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105510965B (zh) | 流体识别方法和装置 | |
CN104502969B (zh) | 河道砂岩性油藏的识别方法 | |
US9279897B2 (en) | Locating oil or gas passively by observing a porous oil and gas saturated system giving off its characteristic resonance response to ambient background noise, including optional differentiation of oil, locatinggas and water | |
CN101551466B (zh) | 一种利用与偏移距有关的地震属性提高油气储层预测精度的方法 | |
US9176243B2 (en) | Locating oil or gas actively by exciting a porous oil and gas saturated system to give off its characteristic resonance response, with optional differentiation of oil, gas and water | |
US6941227B2 (en) | Frequency-dependent processing and interpretation (FDPI) of seismic data for identifying, imaging and monitoring fluid-saturated underground reservoirs | |
CN104407378B (zh) | 一种各向异性参数反演方法及装置 | |
CN105425292B (zh) | 一种油气预测方法及装置 | |
CN103487831B (zh) | Avo地震正演计算方法 | |
CN104155693A (zh) | 储层流体流度的角道集地震响应数值计算方法 | |
Chen et al. | Seismic inversion for underground fractures detection based on effective anisotropy and fluid substitution | |
CN104656136A (zh) | 基于实际模型地震模拟指导下的油气储层低频阴影识别技术 | |
CA2678289C (en) | Passively locating a porous oil, gas or water saturated system giving off its characteristic resonance response to ambient background noise | |
CN106405639B (zh) | 一种叠前地震储层岩性参数的反演方法 | |
CN105938203A (zh) | 一种储层特性的检测方法及装置 | |
CN106353807A (zh) | 裂缝识别方法和装置 | |
Çemen et al. | Correlating porosity with acoustic impedance in sandstone gas reservoirs: examples from the Atokan sandstones of the Arkoma Basin, South Eastern Oklahoma | |
GB2473607A (en) | Locating hydrocarbons passively by observing a porous oil and gas saturated system giving off its characteristic resonance response to ambient background nois | |
CN110456412A (zh) | 一种基于叠后地震数据识别碳酸盐岩储层流体饱和度的方法 | |
CN104820237A (zh) | 计算地层的方差体的方法 | |
CN106324663A (zh) | 一种品质因子的获取方法 | |
CN107797144A (zh) | 基于横波分裂振幅比属性的流体检测方法 | |
Paris et al. | Predicting reservoir quality in the Bakken Formation, North Dakota, using petrophysics and 3C seismic data | |
GB2476788A (en) | Locating fluid saturated zones by applying low frequency excitation and analysing a characteristic resonant response | |
CA2687389C (en) | Locating and differentiating oil, gas or water by actively exciting a porous oil, gas or water saturated system to give off its characteristic resonance response |
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 |