CN103364834A - 一种利用叠前地震频散分析预测储层渗透率的方法 - Google Patents
一种利用叠前地震频散分析预测储层渗透率的方法 Download PDFInfo
- Publication number
- CN103364834A CN103364834A CN2013103216959A CN201310321695A CN103364834A CN 103364834 A CN103364834 A CN 103364834A CN 2013103216959 A CN2013103216959 A CN 2013103216959A CN 201310321695 A CN201310321695 A CN 201310321695A CN 103364834 A CN103364834 A CN 103364834A
- Authority
- CN
- China
- Prior art keywords
- frequency
- different
- reflection coefficient
- permeability
- utilize
- 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.)
- Granted
Links
- RXAHGFBMEOJVKL-UHFFFAOYSA-N C(C12)C3N1C2NC3 Chemical compound C(C12)C3N1C2NC3 RXAHGFBMEOJVKL-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
一种利用叠前地震频散分析预测储层渗透率的方法,包括根据反射系数方程:求取纵波V p随频率f的变化的情况,求取的方法包括利用反射系数方程;对不同频率f,测量不同的反射系数,从随频率不同而不同的反射系数值拟合推导出
Description
技术领域
本发明属于石油勘探领域,涉及一种利用叠前地震频散分析预测储层渗透率的方法。
背景技术
在孔隙弹性理论中,渗透率在地震勘探直接找油气技术中有着重要的地位,它一方面是油气藏描述必不可少的实质性参数,另一方面又是甜点预测的标志性参数。由地震数据中提取渗透率信息是最重要的任务之一。
用地震资料预测渗透率目前尚无成熟方法。生产上大多通过渗透率和孔隙度的回归关系预测。也有文献报道用纵横波速度比和孔隙度预测,驰豫频率预测渗透率等,在生产上尚未见到太多的成功案例。
在孔隙弹性地震参数中流动性是一个重要参数,流动性和渗透率关系密切,所以我们可以探索通过预测流动性参数来预测渗透率。而地震波在孔隙介质传播时,会造成地震波的频散和吸收效应,地震波的频散和流动性有关。
在固流双相介质中,流体的存在导致地震波发生不同程度的频散和衰减[3]。孔隙填充介质的喷射流动(或局部粘滞流动)是引起地震波发生频散和衰减的主要原因。Batzle etal.研究表明,与频率有关的速度变化受岩石的渗透率影响,并定义了流体的流动性,即岩石渗透率与流体粘滞性的比值,来描述P波频散和流体的流动性的关系。Chapman等人在《The influence of fluid –sensitive dispersion and attenuation on AVO analysis》( Chapman, M, Liu, E. and Li, X-Y Geophys. J. Int. 2006, 167, 89-105)一文中从理论上研究了两个弹性层之间夹一个填充流体导致衰减和频散的频散层组成模型的AVO特征,认为弹性层和频散层之间的界面会导致地震反射系数随频率的变化而变化,这种变化与AVO的类型有关。对于第一类AVO,速度频散导致反射波能量集中在高频段;而对于第三类AVO,速度频散导致反射波能量集中在低频段。这一研究表明可以根据反射波能量随频率的变化来研究流体的性质。这里AVO(Amplitude Versus Offset,振幅随偏移距的变化)是一项应用在地质勘探,特别是地下油气资源勘探领域的岩体特性分析技术。
反射系数与入射角的关系符合Zoeppritz(佐普里兹)方程,Zoeppritz角度域各近似方法,都采用了与S波有关的 、泊松比及等参量来表示反射系数。尽管纵波的速度信息可以利用测井数据准确地获得,但对于S波信息未知的地区,利用AVO技术进行岩性参数反演会带来很大误差和不确定性。由于流体对流体的反射系数在垂直入射时S波速度为零,而非垂直入射时反射系数的其他项可以用射线参数来替代,因此Wang Yanghua等人在《Approximations to the Zeoppritz equations and their use in AVO analysis》(Wang Yanghua. Geophysics, 64, 1920-1927)提出用流体因子及射线参数表示近似反射系数。
发明内容
为对叠前地震反演进行渗透率分析,本发明利用频变AVO理论分析纵横波参数随频率变化情况,进而分析流动性和渗透率,提出了一种利用叠前地震频散分析预测储层渗透率的方法。
本发明所述一种利用叠前地震频散分析预测储层渗透率的方法,包括根据反射系数方程:
上述方法利用从射线域对Zoeppritz方程进行简化后的公式进行频率分析,可以得到反射系数随频率变化的规律,从而对岩体渗透率进行分析。
--------② ;
通过上述处理,简化了测量和计算过程。
采用本发明所述的利用叠前地震频散分析预测储层渗透率的方法,可以反映流体速度的频散,通过分析流体速度的频变程度来衡量流体的流动性,进而达到预测渗透率的目的。
本发明与常规AVO分析不同,在射线参数域进行分析,不但考虑了炮检距的变化情况,同时还考虑了频率域的变化特征。
具体实施方式
下面对本发明的具体实施方式作进一步的详细说明。
精确求取反射系数需要借助Zeoppritz方程,而完整的Zeoppritz方程形式非常复杂,难以进行有效的数据分析,前人对该方程进行了各种简化,现有Zeoppritz方程简化多从角度域着手,而Wang Yanghua提出了从射线域简化Zeoppritz方程的思想,本发明中,优选使用Wang Yanghua等人在《Approximations to the Zeoppritz equations and their use in AVO analysis》(Wang Yanghua. Geophysics, 64, 1920-1927)提出的方程:
对反射系数进行简化处理。 相对于角度域的简化方式,射线域在远道由于实现了宽角度和宽范围的检测,得到的信息量更丰富,射线域简化的方程在进行频散分析时可以利用的信息量更多。
尽管纵波V p的速度信息可以利用测井数据准确地获得,但对于横波信息未知的地区,利用AVO技术进行岩性参数反演会带来很大误差和不确定性。由于流体对流体的反射系数在垂直入射时横波速度为零,而非垂直入射时反射系数的其他项可以用射线参数P来替代。①式能够较好的反应上述情况,对横波参数以射线参数代替,简化了反射系数的求取过程,以①式对不同频率测量不同的反射系数,利用测量结果可以观察出流体反射系数随频率变化的规律,、、P、均可以测量得到。
而根据 ----④
可见是的函数,从而可以从随频率不同而不同的反射系数值拟合推导出随频率f变化的规律。
包含纵、横波慢度参数和的项在工程应用中对反射系数影响很小,可以忽略不计,因此在简化过程中去掉,利用Gardener近似方程将参数合并成其他现有参数。使②式中不再出现在①式中出现的纵、横波慢度参数和及,简化了测量和计算过程。
将②式对以一选定的频率作频域展开,由于地震能量主要集中在地震主频带附近,在主频两侧随着频率的变化逐渐减小。频率可以为测量层段地震能量峰值频率为中心,左右各扩展ΔF1的频带内的任意频率值,其中ΔF1为预先设定的扩展频宽:例如对测量层段地震能量峰值频率为75兆,设置ΔF1为5兆,则可以在70-80兆的频带内选择任一频率。
为简化计算,对②式展开后仅取一阶项,可得
本发明利用从射线域简化的Zeoppritz方程形式,从频率域对反射系数进行测量分析,利用反射系数方程得到反射系数方程中能够反应储层渗透率的流体速度参数随频率的变化规律,从而分析储层渗透率。由于地震波的横波部分随频率的变化幅度相对纵波极小,因此流体速度参数只关注纵波的频散程度,只求取纵波的频率方程得到随频率f变化的规律客观上也简化了计算过程。
本发明从各种Zeoppritz方程的简化形式中选取了从射线域简化的WANG YANG HUA 方程①式,并对其进一步进行了简化处理,使之可以被用来进行频散分析流体速度的相关参数。利用简化后的方程进行展开并只对一阶项讨论,得到流体速度随频率的变化规律。
采用本发明所述的利用叠前地震频散分析预测储层渗透率的方法,可以反映流体速度的频散,通过分析流体速度的频变程度来衡量流体的流动性,进而达到预测渗透率的目的。
本发明与常规AVO分析不同,在射线参数域进行分析,不但考虑了炮检距的变化情况,同时还考虑了频率域的变化特征。
本发明中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块,或者二者的结合来实施。软件模块可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。
前文所述的为本发明的各个优选实施例,各个优选实施例中的优选实施方式如果不是明显自相矛盾或以某一优选实施方式为前提,各个优选实施方式都可以任意叠加组合使用,所述实施例以及实施例中的具体参数仅是为了清楚表述发明人的发明验证过程,并非用以限制本发明的专利保护范围,本发明的专利保护范围仍然以其权利要求书为准,凡是运用本发明的说明书内容所作的等同结构变化,同理均应包含在本发明的保护范围内。
Claims (5)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310321695.9A CN103364834B (zh) | 2013-07-29 | 2013-07-29 | 一种利用叠前地震频散分析预测储层渗透率的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310321695.9A CN103364834B (zh) | 2013-07-29 | 2013-07-29 | 一种利用叠前地震频散分析预测储层渗透率的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103364834A true CN103364834A (zh) | 2013-10-23 |
CN103364834B CN103364834B (zh) | 2016-08-10 |
Family
ID=49366601
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310321695.9A Active CN103364834B (zh) | 2013-07-29 | 2013-07-29 | 一种利用叠前地震频散分析预测储层渗透率的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103364834B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105510965A (zh) * | 2015-12-29 | 2016-04-20 | 中国石油天然气股份有限公司 | 流体识别方法和装置 |
CN111399044A (zh) * | 2020-04-13 | 2020-07-10 | 中国石油大学(北京) | 一种储层渗透率预测方法、装置及存储介质 |
CN112558153A (zh) * | 2019-09-25 | 2021-03-26 | 中国石油天然气股份有限公司 | 一种双相介质的油气储层预测方法及装置 |
CN113031068A (zh) * | 2021-02-24 | 2021-06-25 | 浙江大学 | 一种基于反射系数精确式的基追踪叠前地震反演方法 |
CN114185097A (zh) * | 2021-08-16 | 2022-03-15 | 中国石油天然气集团有限公司 | 预测钻前潜在井漏风险层的方法及装置 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5555218A (en) * | 1995-02-27 | 1996-09-10 | Western Atlas International, Inc. | Computation of Q-derived static corrections from seismic data |
CN1237256A (zh) * | 1997-01-31 | 1999-12-01 | 菲利浦石油公司 | 用于测定油层渗透率、孔隙率和伪相对渗透率的分布的方法 |
CN101109821A (zh) * | 2007-08-16 | 2008-01-23 | 中国石化集团胜利石油管理局 | 基于系统辨识提高地震资料分辨率的方法 |
US7672784B2 (en) * | 2006-06-21 | 2010-03-02 | Baker Hughes Incorporated | Using statistics of a fitting function for data-driven dispersion slowness processing |
CN101738637A (zh) * | 2008-11-06 | 2010-06-16 | 北京北方林泰石油科技有限公司 | 一种基于速度随频率变化信息的油气检测方法 |
CN102478667A (zh) * | 2010-11-30 | 2012-05-30 | 中国石油天然气集团公司 | 频散avo反演方法 |
CN102590863A (zh) * | 2012-01-29 | 2012-07-18 | 中国石油天然气股份有限公司 | 确定反射系数频散、反射能量及识别气层的方法及装置 |
-
2013
- 2013-07-29 CN CN201310321695.9A patent/CN103364834B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5555218A (en) * | 1995-02-27 | 1996-09-10 | Western Atlas International, Inc. | Computation of Q-derived static corrections from seismic data |
CN1237256A (zh) * | 1997-01-31 | 1999-12-01 | 菲利浦石油公司 | 用于测定油层渗透率、孔隙率和伪相对渗透率的分布的方法 |
US7672784B2 (en) * | 2006-06-21 | 2010-03-02 | Baker Hughes Incorporated | Using statistics of a fitting function for data-driven dispersion slowness processing |
CN101109821A (zh) * | 2007-08-16 | 2008-01-23 | 中国石化集团胜利石油管理局 | 基于系统辨识提高地震资料分辨率的方法 |
CN101738637A (zh) * | 2008-11-06 | 2010-06-16 | 北京北方林泰石油科技有限公司 | 一种基于速度随频率变化信息的油气检测方法 |
CN102478667A (zh) * | 2010-11-30 | 2012-05-30 | 中国石油天然气集团公司 | 频散avo反演方法 |
CN102590863A (zh) * | 2012-01-29 | 2012-07-18 | 中国石油天然气股份有限公司 | 确定反射系数频散、反射能量及识别气层的方法及装置 |
Non-Patent Citations (3)
Title |
---|
YANGHUA WANG: "Approximations to the Zoeppritz equations and their use in AVO analysis", 《GEOPHYSICS》, vol. 64, no. 6, 30 November 1999 (1999-11-30), pages 1920 - 1927 * |
刘力辉等: "一种改进的射线弹性阻抗公式及弹性参数反演", 《石油物探》, vol. 50, no. 4, 25 July 2011 (2011-07-25), pages 331 - 335 * |
刘力辉等: "叠前AVO属性的地震岩性学探索与实践研究", 《石油物探》, vol. 52, no. 3, 25 May 2013 (2013-05-25), pages 247 - 252 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105510965A (zh) * | 2015-12-29 | 2016-04-20 | 中国石油天然气股份有限公司 | 流体识别方法和装置 |
CN105510965B (zh) * | 2015-12-29 | 2018-01-05 | 中国石油天然气股份有限公司 | 流体识别方法和装置 |
CN112558153A (zh) * | 2019-09-25 | 2021-03-26 | 中国石油天然气股份有限公司 | 一种双相介质的油气储层预测方法及装置 |
CN111399044A (zh) * | 2020-04-13 | 2020-07-10 | 中国石油大学(北京) | 一种储层渗透率预测方法、装置及存储介质 |
CN111399044B (zh) * | 2020-04-13 | 2021-05-25 | 中国石油大学(北京) | 一种储层渗透率预测方法、装置及存储介质 |
CN113031068A (zh) * | 2021-02-24 | 2021-06-25 | 浙江大学 | 一种基于反射系数精确式的基追踪叠前地震反演方法 |
CN113031068B (zh) * | 2021-02-24 | 2022-05-27 | 浙江大学 | 一种基于反射系数精确式的基追踪叠前地震反演方法 |
CN114185097A (zh) * | 2021-08-16 | 2022-03-15 | 中国石油天然气集团有限公司 | 预测钻前潜在井漏风险层的方法及装置 |
CN114185097B (zh) * | 2021-08-16 | 2024-03-12 | 中国石油天然气集团有限公司 | 预测钻前潜在井漏风险层的方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN103364834B (zh) | 2016-08-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Foti et al. | Application of surface-wave methods for seismic site characterization | |
Pitarka et al. | Simulation of near-fault strong-ground motion using hybrid Green's functions | |
Roux et al. | P-waves from cross-correlation of seismic noise | |
Zerva et al. | Spatial variation of seismic ground motions: An overview | |
Zheng et al. | Seismic characterization of fractured reservoirs by focusing Gaussian beams | |
Hailemikael et al. | Ground-motion amplification at the Colle di Roio ridge, central Italy: a combined effect of stratigraphy and topography | |
Picotti et al. | Anisotropy and crystalline fabric of Whillans Ice Stream (West Antarctica) inferred from multicomponent seismic data | |
Pitarka et al. | Analysis of ground motion from an underground chemical explosion | |
Fehler et al. | Envelope broadening of outgoing waves in 2D random media: a comparison between the Markov approximation and numerical simulations | |
Titov et al. | Modeling and interpretation of scattered waves in interstage distributed acoustic sensing vertical seismic profiling survey | |
Baird et al. | Frequency-dependent seismic anisotropy due to fractures: Fluid flow versus scattering | |
CN104155693A (zh) | 储层流体流度的角道集地震响应数值计算方法 | |
CN103364834A (zh) | 一种利用叠前地震频散分析预测储层渗透率的方法 | |
Coutel et al. | Simulation-based comparison of four site-response estimation techniques | |
CN103487831A (zh) | Avo地震正演计算方法 | |
CN110618450B (zh) | 基于岩石物理建模的致密储层智能化含气性预测方法 | |
Saccorotti et al. | A review of seismic methods for monitoring and understanding active volcanoes | |
Li et al. | Identification of fractured carbonate vuggy reservoirs in the S48 well area using 3D 3C seismic technique: A case history from the Tarim Basin | |
Ziane et al. | The contribution of multiple scattering to Love wave generation in the secondary microseism | |
Zhang et al. | Retrieval of shallow S-wave profiles from seismic reflection surveying and traffic-induced noise | |
Picotti et al. | P-wave seismic attenuation by slow-wave diffusion: Numerical experiments in partially saturated rocks | |
Wang et al. | Anisotropic characteristics of mesoscale fractures and applications to wide azimuth 3D P-wave seismic data | |
Tang et al. | Single-point dispersion measurement of surface waves combining translation, rotation and strain in weakly anisotropic media: theory | |
Li et al. | Relative elastic interferometric imaging for microseismic source location | |
Kosloff et al. | Two-dimensional simulation of Rayleigh waves with staggered sine/cosine transforms and variable grid spacing |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |