CN106249294A - 一种储层烃类检测方法 - Google Patents
一种储层烃类检测方法 Download PDFInfo
- Publication number
- CN106249294A CN106249294A CN201510320573.7A CN201510320573A CN106249294A CN 106249294 A CN106249294 A CN 106249294A CN 201510320573 A CN201510320573 A CN 201510320573A CN 106249294 A CN106249294 A CN 106249294A
- Authority
- CN
- China
- Prior art keywords
- reservoir
- theta
- frequency
- wave
- incidence
- 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.)
- Pending
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种储层烃类检测方法,其特征在于,包括;获得地下介质的粘弹性参数,包括:纵波速度、横波速度、密度、纵波衰减以及横波衰减参数;根据粘弹性介质AVAF关系近似式,计算出反射系数随入射角、频率的变化关系;根据所述反射系数随入射角、频率的变化关系,获得储层发育的详细信息,由此推测出储层的烃类含量。本发明在常规的振幅随入射角变化关系(AVA)的基础上引入了频率对地震振幅的影响,能更加方便、精确地描述发生在界面处的反射,易于实现油气勘探领域中储层的精细描述,可以准确推测出储层的烃类含量。
Description
技术领域
本发明属于石油地球物理勘探领域,具体涉及一种使用粘弹性介质振幅随入射角及频率变化关系的储层烃类检测方法。
背景技术
可靠的烃类检测可降低勘探风险与投入成本,对储层评价和油气勘探具有重要意义。自20世纪70年代以来,基于地震振幅差异和储层弹性参数的烃类检测技术一直是主导应用技术,如最初的“亮点”和AVO分析技术,以及后来的LMR(拉梅参数)技术、Russell流体因子和泊松阻抗方法等。这些技术虽然在提高油气勘探成功率上做出了巨大的贡献,但随着勘探主体变为隐蔽岩性油气藏或深海油气藏,单纯依靠振幅异常和储层弹性参数的烃类检测技术在实际储层预测中逐渐显露出局限性。
首先,隐蔽岩性油气藏的气水关系复杂、储层非均质性强、地震振幅响应与弹性参数异常特征不明显,常规方法的多解性严重;另外,海上钻井较少,叠前地震反演的实施难度较大,影响了弹性参数的精确求取。
发明内容
本发明的目的在于解决上述现有技术中存在的难题,提供了一种储层烃类检测方法,可以利用粘弹性介质振幅随入射角及频率变化关系(AVAF),更好地描述反射波随频率的变化关系,实现对储层烃类的有效检测,提高储层含油气性预测的准确性。
本发明是通过以下技术方案实现的:
获得地下介质的粘弹性参数,包括:纵波速度、横波速度、密度、纵波衰减以及横波衰减参数;
根据下列公式,计算出反射系数随入射角、频率的变化关系:
其中,Vp表示介质的纵波速度,Vs表示介质的横波速度,ρ表示介质的密度,θ分别表示地震波由上覆地层进入下方地层的入射角,ω表示地震波的频率,为频率相关项,αp表示纵波衰减因子,αs表示横波衰减因子;
根据所述反射角随入射角、频率的变化关系,获得储层发育的详细信息,由此推测出储层的烃类含量。
与现有技术相比,本发明的有益效果是:本发明在常规的振幅随入射角变化关系(AVA)的基础上引入了频率对地震振幅的影响,即AVAF。首先,粘弹性介质的假设更符合地下介质的实际情况,这对理论进行了完善;其次,频率的引入将为基于地震数据的反演和分析带来更多有用的信息。粘弹性介质AVAF近似式将更加方便、精确地描述发生在界面处的反射,为在粘弹性介质反演与分析奠定良好的理论基础,易于实现油气勘探领域中储层的精细描述,可以准确推测出储层的烃类含量。
附图说明
图1为本发明的储层烃类检测方法的流程图。
图2为含气砂岩储层模型示意图;
图3为储层顶界面精确反射系数随入射角和频率的变化;
图4为储层顶界面近似反射系数随入射角和频率的变化;
图5为储层顶界面反射系数误差分析;
具体实施方式
下面结合附图对本发明作进一步详细描述。
发明人通过深入的研究发现,地震波在含油气储层中传播时会发生能量衰减异常,可利用衰减特征参数进行储层烃类检测,这样不仅能较好地弥补无法实施传统方法的不足,而且可与其他相关方法相互补充验证,降低预测结果的多解性。
基于此,本发明提出了一种储层烃类检测方法,图1示出了本发明的储层烃类检测方法的流程图。
参考图1,在步骤101中,获得地下介质的粘弹性参数,包括:纵波速度、横波速度、密度、纵波衰减以及横波衰减参数。
在步骤102中,根据粘弹性介质AVAF关系近似式,计算出反射角随入射角、频率的变化关系。
首先,由运动方程和虎克定律可推导出弹性波方程,在此基础上结合连续性方程、平面波的解、斯奈尔定律可得到Zoeppritz方程:
其中,VP,VS,ρ依次分别表示介质的纵波速度,横波速度和密度,θ表示入射角,RPP,RPS,TPP,TPS依次分别表示PP波反射系数、PS波反射系数,PP波透射系数和PS波透射系数。
对于粘弹性介质,相速度V(ω)是包含频率信息的复数,写为:
其中,H(ω)和F(ω)分别为新定义参数,V(ω)是相速度,ωr是参考频率,Vr是参考频率的相速度,Q为地下介质的衰减,请发明人补充方程(2)中参数的含义。需要注意的是,式中的速度和衰减即包含了纵波也包含了横波。当反射界面上下的速度差不是很大时,由传统的Aki-Richards近似可以得到:
其中,VP,VS,ρ依次分别表示介质的纵波速度、横波速度和密度,θ表示入射角,Δ表示上下介质的差,
将(2)式代入(3),
其中,R(θ,ω)表示随入射角和频率变换的反射系数,Vp表示介质的纵波速度,Vs表示介质的横波速度,ρ表示介质的密度,θ分别表示地震波由上层介质进入下层介质的入射角,ω表示地震波的频率,下标1和2分别表示上下介质,下标r表示参数H(ω)的实部(关于H(ω)的定义见式2)
对(4)式中的速度项进行泰勒(Taylor)展开,展开条件如下:
其中符号“→”表示“趋近于”。
在(5)条件下将纵波速度反射和横波速度反射展开成:
将(6)式代入Aki-Richards近似式,可以得到粘弹性介质下线性化的AVAF关系式:
其中,HP(ω)和HS(ω)分别表示与P波和S波相关的参数,H(ω)的定义仍然参照式(2)。
进一步线性化,得到:
其中α=Q-1,由此(7)式写为:
方程式(9)即为粘弹性介质AVAF关系近似式,其中,Vp表示介质的纵波速度,Vs表示介质的横波速度,ρ表示介质的密度,θ分别表示入射角,ω表示频率,为频率相关项,αp表示纵波衰减因子,αs表示横波衰减因子。
利用在步骤101中获得的纵波速度、横波速度、密度、纵波衰减以及横波衰减参数,根据公式(9),可以计算出反射系数随入射角、频率的变化关系,即得到一个三维数据体(即由入射角度域、频率域和时间域组成的三维空间,在这个空间内的每一个点都对应一个数据值,即三维数据体),由此可获得储层发育的详细信息。根据反射系数随入射角和频率的变化关系可以推知是否有储层发育以及是何种类型储层,例如:如看到反射系数随入射角增大而逐渐增大,那么该层段极有可能是含气性储层;同理,如看到反射系数随入射角增大而出现先增大后减小态势,则为含油储层,由此即实现了对储层烃类的检测作用。
下面通过一个含气砂岩储层模型来说明本发明的效果。模型为泥岩层中间夹有含气砂岩储层(如图1所示),三层界面的参数分别设为:
Vp1=3000m/s,Vs1=1500m/s,ρ1=2700kg/m3,Qp1=100,Qs1=150;
Vp2=2800m/s,Vs2=1200m/s,ρ2=2650kg/m3,Qp2=10,Qs2=20;
Vp3=3500m/s,Vs3=1800m/s,ρ3=2750kg/m3,Qp3=100,Qs3=150;
其中,Vp,Vs,ρ,Qp和Qs依次分别代表介质中纵波速度、横波速度、密度、纵波品质因子和横波品质因子,下标代表不同的层。
储层顶界面位于100ms处,实验模拟的入射角度范围从0°—50°变化,频率范围为:0hz—100hz。
图2和图3分别为储层顶界面精确AVAF关系与近似AVAF关系,仅从AVO分类上二者都体现出III类AVO效应,如果同时考虑频率引起的变换,可以看到当入射角增大时,反射系数随频率增加而减小,绝对值增大。
图4为本发明近似式与精确反射系数之间的误差,在常规勘探的入射角和频带范围内,二者的误差很小,几乎可以忽略。说明本发明在地球物理勘探领域具有良好的适用性。
本发明针对传统的AVO或AVA分析存在的缺陷,专利提出了一种粘弹性介质振幅随入射角和和频率变化的近似技术。该技术基于粘弹性介质理论,可有效分析来自储层顶、底界面反射系数随入射角和频率的变化规律,实现储层的精细刻画与描述。同时,近似式的表达具有较高的准确性,为通过利用角度和频率信息进行反演,实现储层定量化描述奠定了良好的基础。
上述技术方案只是本发明的一种实施方式,对于本领域内的技术人员而言,在本发明公开了应用方法和原理的基础上,很容易做出各种类型的改进或变形,而不仅限于本发明上述具体实施方式所描述的方法,因此前面描述的方式只是优选的,而并不具有限制性的意义。
Claims (4)
1.一种储层烃类检测方法,其特征在于,包括以下步骤;
获得地下介质的粘弹性参数,包括:纵波速度、横波速度、密度、纵波衰减以及横波衰减参数;
根据粘弹性介质AVAF关系近似式,计算出反射系数随入射角、频率的变化关系:
其中,Vp表示介质的纵波速度,Vs表示介质的横波速度,ρ表示介质的密度,θ分别表示地震波由上覆地层进入下方地层的入射角,ω表示地震波的频率,为频率相关项,αp表示纵波衰减因子,αs表示横波衰减因子;
根据所述反射系数随入射角、频率的变化关系,获得储层发育的详细信息,由此推测出储层的烃类含量。
2.根据权利要求1所述的储层烃类检测方法,其特征在于,所述根据所述反射系数随入射角、频率的变化关系,获得储层发育的详细信息可以包括:所述根据所述反射系数随入射角、频率的变化关系,得到三维数据体,利用该三维数据体,获得储层发育的详细信息。
3.根据权利要求1所述的储层烃类检测方法,其特征在于,所述储层发育的详细信息包括:储层发育的孔隙度,以及所含的流体及储层是否属于优质储层。
4.根据权利要求1所述的储层烃类检测方法,其特征在于,所述推测出储层的烃类含量包括:
根据所述反射系数随入射角、频率的变化关系,反演出地震波在储层中传播时的能量衰减异常,由此推测出储层的烃类含量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510320573.7A CN106249294A (zh) | 2015-06-12 | 2015-06-12 | 一种储层烃类检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510320573.7A CN106249294A (zh) | 2015-06-12 | 2015-06-12 | 一种储层烃类检测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106249294A true CN106249294A (zh) | 2016-12-21 |
Family
ID=57626299
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510320573.7A Pending CN106249294A (zh) | 2015-06-12 | 2015-06-12 | 一种储层烃类检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106249294A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107479104A (zh) * | 2017-08-04 | 2017-12-15 | 成都晶石石油科技有限公司 | 泊松反射率及泊松流体反射率的求取方法及地层分析方法 |
CN108957556A (zh) * | 2018-08-13 | 2018-12-07 | 西北大学 | 非气藏地震波衰减异常的地质成因排除方法 |
CN111427086A (zh) * | 2020-04-02 | 2020-07-17 | 同济大学 | 利用纵横波振幅频率分量预测地层最大衰减参数的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2138135A (en) * | 1983-04-06 | 1984-10-17 | Chevron Res | Interpretation of seismic records |
CN102749645A (zh) * | 2012-03-14 | 2012-10-24 | 中国石油天然气股份有限公司 | 一种利用角度阻抗梯度进行储层烃类检测的方法及装置 |
CN103257361A (zh) * | 2013-05-24 | 2013-08-21 | 中国石油天然气集团公司 | 基于Zoeppritz方程近似式的油气预测方法及系统 |
CN103643949A (zh) * | 2013-12-20 | 2014-03-19 | 中国石油天然气集团公司 | 一种储层含油气的定量预测方法及装置 |
-
2015
- 2015-06-12 CN CN201510320573.7A patent/CN106249294A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2138135A (en) * | 1983-04-06 | 1984-10-17 | Chevron Res | Interpretation of seismic records |
CN102749645A (zh) * | 2012-03-14 | 2012-10-24 | 中国石油天然气股份有限公司 | 一种利用角度阻抗梯度进行储层烃类检测的方法及装置 |
CN103257361A (zh) * | 2013-05-24 | 2013-08-21 | 中国石油天然气集团公司 | 基于Zoeppritz方程近似式的油气预测方法及系统 |
CN103643949A (zh) * | 2013-12-20 | 2014-03-19 | 中国石油天然气集团公司 | 一种储层含油气的定量预测方法及装置 |
Non-Patent Citations (5)
Title |
---|
HENG LUO 等: "A linearized AVAF equation and Q impedance in anelastic media", 《CPS/SEG BEIJING 2014 INTERNATIONAL GEOPHYSICAL CONFERENCE》 * |
乔玉雷: "《三维垂直地震剖面资料处理与解释技术-以胜利油田垦71井区为例》", 31 May 2011, 石油工业出版社 * |
朱光亚 等: "《中国科学技术文库》", 31 January 1998, 科学技术文献出版社 * |
李雪英 等: "薄层时频特征的正演模拟", 《地球物理学报》 * |
陈怀震 等: "基于各向异性岩石物理的缝隙流体因子AVAZ反演", 《地球物理学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107479104A (zh) * | 2017-08-04 | 2017-12-15 | 成都晶石石油科技有限公司 | 泊松反射率及泊松流体反射率的求取方法及地层分析方法 |
CN108957556A (zh) * | 2018-08-13 | 2018-12-07 | 西北大学 | 非气藏地震波衰减异常的地质成因排除方法 |
CN111427086A (zh) * | 2020-04-02 | 2020-07-17 | 同济大学 | 利用纵横波振幅频率分量预测地层最大衰减参数的方法 |
CN111427086B (zh) * | 2020-04-02 | 2021-03-26 | 同济大学 | 利用纵横波振幅频率分量预测地层最大衰减参数的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104570072B (zh) | 一种粘弹性介质中的球面pp波反射系数建模方法 | |
Hornby et al. | Interferometric imaging of a salt flank using walkaway VSP data | |
CN106842326B (zh) | 无横波速度测井时砂泥互层co2地质封存时移地震正演模拟方法 | |
CN110618450B (zh) | 基于岩石物理建模的致密储层智能化含气性预测方法 | |
Zeng et al. | Recent progress in analysis of seismically thin beds | |
CN110187389B (zh) | 一种基于薄层反射理论的ava反演方法 | |
Leiceaga et al. | Enhanced density estimation from prestack inversion of multicomponent seismic data | |
CN106249294A (zh) | 一种储层烃类检测方法 | |
CN103472482B (zh) | 基于基因排序体系的多波地震资料时间域匹配方法及系统 | |
Kadkhodaie et al. | Pore throat size characterization of carbonate reservoirs by integrating core data, well logs and seismic attributes | |
CN105093291A (zh) | 一种恢复油气储层地震反射特征的方法 | |
Landrø et al. | Using diving waves for detecting shallow overburden gas layers | |
Wallick et al. | Interpretation of full-azimuth broadband land data from Saudi Arabia and implications for improved inversion, reservoir characterization, and exploration | |
Gray et al. | Fractured reservoir characterization using avaz on the Pinedale anticline, Wyoming | |
Naseer et al. | Delineation of stratigraphic traps within the basin floor fans of Miocene sedimentary sequences, offshore Indus, Pakistan using inverted acoustic impedance simulations | |
CN109507725B (zh) | 富砂层段有效储层地震属性预测方法及系统 | |
Holman | Seismic characterization of fractured rock fabric in Mississippian limestone, Payne county Oklahoma | |
Strecker et al. | Seismic attribute analysis in hydrothermal dolomite, Devonian slave point formation, Northeast British Columbia, Canada | |
Sumantri et al. | Seismic Qualitative and Quantitative Analysis for Reservoir Characterization of Globigerina BH Gas Field, North East Java Basin. | |
Hussein et al. | Prospects delineation and evaluation using the CDP gathers characterization, onshore Nile Delta, Egypt | |
Xu | Oil reservoir assessment using multicomponent seismic data | |
Wang et al. | An introduction to seismic exploration of the Micang-Dabashan foothill belt in the Sichuan Basin | |
Ren et al. | AVO influence factors exploration and case analysis | |
Wei et al. | High quality reservoir prediction of sublacustrine fan in Liaozhong Depression | |
Zeng et al. | Research and application of integrated technology of tight sandstone gas acquisition processing and interpretation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20161221 |
|
RJ01 | Rejection of invention patent application after publication |