CN113534241B - 一种方位频变流体因子识别储层的方法 - Google Patents
一种方位频变流体因子识别储层的方法 Download PDFInfo
- Publication number
- CN113534241B CN113534241B CN202110779032.6A CN202110779032A CN113534241B CN 113534241 B CN113534241 B CN 113534241B CN 202110779032 A CN202110779032 A CN 202110779032A CN 113534241 B CN113534241 B CN 113534241B
- Authority
- CN
- China
- Prior art keywords
- reservoir
- frequency
- formula
- fluid
- dispersion
- 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 49
- 238000000034 method Methods 0.000 title claims abstract description 17
- 230000001419 dependent effect Effects 0.000 title claims description 11
- 239000006185 dispersion Substances 0.000 claims abstract description 42
- 230000035945 sensitivity Effects 0.000 claims abstract description 6
- 230000001360 synchronised effect Effects 0.000 claims abstract description 5
- 239000011435 rock Substances 0.000 claims description 13
- 230000035699 permeability Effects 0.000 claims description 10
- 239000011159 matrix material Substances 0.000 claims description 8
- 150000001875 compounds Chemical class 0.000 claims description 6
- 239000013598 vector Substances 0.000 claims description 6
- 238000011160 research Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 7
- 230000002159 abnormal effect Effects 0.000 abstract description 3
- 238000001228 spectrum Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000005856 abnormality Effects 0.000 description 2
- 230000002547 anomalous effect Effects 0.000 description 2
- IAZDPXIOMUYVGZ-UHFFFAOYSA-N Dimethylsulphoxide Chemical compound CS(C)=O IAZDPXIOMUYVGZ-UHFFFAOYSA-N 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000010206 sensitivity analysis Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
Images
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/282—Application of seismic models, synthetic seismograms
-
- 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/301—Analysis for determining seismic cross-sections or geostructures
-
- 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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (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
技术领域
本发明属于油气地震勘探数据处理与解释领域,是一种通过新的与频率相关的方位AVO公式实现新频散流体因子的反演,利用与储层敏感的新频散参数和宽方位地震数据实现高精度储层识别的解释技术。
背景技术
频散是含流体储层的固有特性,可充分挖掘该属性进行储层识别,依赖频率的AVO反演方法为从实际地震数据中提取与含流体储层有关的频散异常提供了一种有效技术。而如今依赖频率的AVO反演技术是强行将频率引入常规AVO公式进行频散反演,具有不合理性,需要利用直接与频率相关的AVO公式进行反演,并且需要构建与储层敏感的流体因子来提高储层预测的精度。
不同频散参数对流体的敏感性存在差异,考虑到频散属性对流体的敏感性问题,基于敏感性分析的频散参数反演可获得与储层流体敏感的频散因子,进而提高储层预测精度。随着宽方位地震采集技术的发展,可充分利用宽方位叠前地震数据进行频散参数的反演进而进行储层的高精度预测。为此,本发明技术提供了一种通过构建对储层敏感的频散因子并结合宽方位地震数据进行储层识别的方法。
发明内容
本发明是要提供一种方位频变流体因子识别储层的方法,通过推导的一种与频率直接相关的方位AVO公式,实现五类频散因子的同步反演,并利用构建的与储层渗透率、流体粘度等参数有关的敏感频散流体因子直接识别储层。
本发明的目的在于提供一种方位频变流体因子识别储层的方法,所述方法主要包括以下步骤:
⑴对叠前角道集数据进行时频变换,并构建与新频散流体因子相关的目标函数;
r=eD
式中,r为m个频率下反射系数的列向量,共m×n行,e为m个频率下反射系数公式中呈m×n行、5列排列的系数矩阵,DM、Dμ、和分别为纵波模量、横波模量、法向弱度、新流体因子和切向弱度的频散程度。M和μ分别表示纵波模量和横波模量,δN和δT分别为岩石的法向和切向弱度,ψn为与裂缝弱度与松弛时间(受流体粘度、岩石渗透率、裂隙纵横比等参数的影响)的平方相关的参数,该参数即为构建的一种新的流体识别因子,表达式为:ψn=δNΓ2,
式中:
式中,ηf为流体的粘度,k为岩石的渗透率,Kf为流体的体积模量,a为裂隙纵横比,Γ为松弛时间。
⑵选取参考频率ω0,利用最小二乘求解该参考频率下五类模量对应的频散结果:P、S、N、K、L,并将该结果带入向量r中,利用下式进行五类频散因子的同步反演:
D=(eTe)-1eTr
⑷重复上述步骤⑴至⑶,直到处理完整个工区的所有叠前地震角道集,即可获得整个工区的三维储层预测结果,该预测结果具有很高的精度,能够压制与储层无关的背景异常,对储层具有很强的敏感性,可对后续的储层精细评价、钻井部署等高精度勘探提供支撑。
本发明的具体实现原理如下:
基于Hudson et al(1996)的等效介质模型,通过简化刚度矩阵参数,构建与裂缝弱度和新流体指示因子有关的新刚度矩阵参数,再利用散射函数与反射系数间的关系推导与频率相关的AVOA反射系数公式。对于HTI介质,其刚度矩阵C可表示为:
式中,M=λ+2μ,e为裂隙密度,λ和μ为岩石的拉梅系数,和U33为与流体参数(流体粘度,体积模量)和裂隙性质(渗透率,裂隙纵横比)有关的重要参数。此处,我们考虑裂隙中充填有流体,则和U33的表达式为:
经过分析表明,U33的虚部响度实部很小,可以忽略。进一步求取U33实部的近似表达式:
则可进一步得到刚度矩阵元素的表达式:
C11=M-MδN-ω2Mψn,
C12=λ-λδN-ω2λψn,
C23=λ-λ(1-2g)δN-ω2λ(1-2g)ψn,
C33=M-λ(1-2g)δN-ω2λ(1-2g)ψn,
C44=μ,C55=μ-μδT,
C13=C21=C31=C12,C32=C23,C66=C55,C22=C33
式中,ψn=δNΓ2,ψn为与裂缝弱度与松弛时间(受流体粘度、岩石渗透率、裂隙纵横比等参数的影响)的平方相关的参数,该参数可以作为一种新的流体识别因子,δN=4e/[3g(1-g)]和δT=16e/[3(3-2g)]分别为岩石的法向和切向弱度。
可进一步得到刚度矩阵元素中扰动:
ΔC11=(M+ΔM)-(M+ΔM)(δN+ΔδN)-ω2(M+ΔM)(ψn+Δψn)-(M-MδN-ω2Mψn)
≈ΔM-δNΔM-MΔδN-ω2Mψn-ω2ψnΔM
ΔC12≈Δλ-ΔλδN-λΔδN-ω2λΔψn-ω2ψnΔλ,
ΔC23≈Δλ-(1-2g)δNΔλ-λ(1-2g)ΔδN-ω2λ(1-2g)Δψn-ω2(1-2g)ψnΔλ,
ΔC33≈ΔM-(1-2g)δNΔλ-λ(1-2g)ΔδN-ω2λ(1-2g)Δψn-ω2(1-2g)ψnΔλ,
ΔC44=Δμ,
ΔC55≈Δμ-δTΔμ-μΔδT,
ΔC13=ΔC21=ΔC31=ΔC12,
ΔC32=ΔC23,ΔC66=ΔC55,ΔC22=ΔC33
式中,ΔM和Δμ分别为反射界面两侧纵波模量和剪切模量的变化量,ΔδN和Δψn为界面两侧δN和新流体因子ψn的变化量。此处假设界面两侧的扰动量较小,则可忽略ΔMΔδN、ΔMΔψn项。
散射函数的一般表达式为:
式中,Pi和Ps分别为入射波和散射波的极化矢量,可表示为:
则有上述公式可得:
上式即为推导的新的纵波反射系数公式,它同时与入射角、方位角和频率有关,包含了与流体特性直接相关的流体因子ψn。该公式可用于分析依赖频率的方位AVO反演。
再将上述推导的公式进行拓展,忽略密度项并将在参考频率ω0处进行泰勒展开,则有:
式中,
A′(θ)=A(θ),B′(θ)=B(θ),
令:
定义e为m×n行、5列的矩阵:
则可得到目标函数:
r=eD
其中,P、S、N、K、L的求解过程为:
a)将A(θi),B(θi),C(θi),D(θi)和E(θi)看成采样时间t和接收道n(每个采样点对应一个系数A,B、C、D和E)的函数,则:
A(θi)=A(t,n),B(θi)=B(t,n),C(θi)=C(t,n),D(θi)=D(t,n),E(θi)=E(t,n)
并对叠前AVO道集进行谱均衡处理:
b)对于某一参考频率ω0,有如下表达形式:
为了寻找使上式的误差达到最小值的变量P,S,N,K和L,利用上式分别对变量P,S,N,K和L,求偏导数如下:
令上式等于0,则有:
求解上式即可得P、S、N、K、L。
同样,给定不同的方位角,即可获得不同方位角时的P、S、N、K、L,进而可求得不同方位的频散因子。
本发明的一种方位频变流体因子识别储层的方法,具有如下优越性:
(1)推导的新的方位AVO反射系数公式直接与频率相关,可直接用于依赖频率的AVO反演,避免了常规依赖频率AVO反演中强行将频率引入AVO公式的不合理性;
(2)构建了直接与储层岩石渗透率、流体粘度有关的新流体因子,该因子对含流体储层具有更强的敏感性,能够用于高精度储层预测;
(3)可同时实现五类频散因子的同步反演,得到的新频散因子对储层具有很好的识别能力,能够压制与储层无关的异常干扰,直接突出储层异常,并且分辨率高,横向连续性好,能够解决薄储层难以识别的问题;
(4)可实现宽方位的频散因子反演,利用宽方位地震数据能够更好的分析裂缝型储层的特征,裂缝弱度等参数对识别致密裂缝型储层具有很好的效果。
附图说明
图1是本发明技术利用的与频率和方位相关的反射系数对比图。
图2是本发明技术反演得到的五类频散因子的单道结果对比图。
图3是本发明实施例工区反演得到的新频散因子剖面。
具体实施方式
本发明的具体实施方式如下:
⑴输入叠前角道集,并对角道集进行时频谱分析,生成不同频率下的振幅谱;
⑵对不同频率下的振幅谱进行振幅谱均衡处理,并获取地震数据的层速度;
⑶选定参考频率利用最小二乘反演该频率下模量的变化率,即求取P、S、N、K、L;
⑸重复步骤⑴至⑷即可获得所有方位时的频散因子结果;
⑹利用步骤⑸得到的新频散流体因子结果进行储层的识别;
⑺重复步骤⑴至⑹,处理研究工区所有的三维数据并对其进行储层预测并抽取沿层切片进行分析,即可得到整个工区不同目的层段的储层预测平面展布。
为了更清楚的表达本发明的技术优势,结合附图度对本发明实施例做进一步的详细说明。本发明的实施例说明:
图1为不同频率时的反射系数随入射角和方位角的变化对比。分别对频率为10Hz、20Hz、30Hz、40Hz、50Hz情况下的反射系数进行对比,由图可知,反射系数随方位角的变化在180°位置处对称,方位角在0~180°的反射系数也在方位角为90°的位置对称,此外,可以明显看出不同频率成分的反射系数随方位角和入射角变化的差异较大,说明反射系数与频率有较强的敏感性,可利用此特性进行依赖频率的方位AVO反演。
图2为利用本发明技术反演得到的五类频散因子单道结果对比图。由图可知,图a-d中的四个频散因子均受明显的无关异常干扰(黑色箭头所示),难以区分有利储层,而图e中新构建的频散因子最大限度地压制了其余层的异常干扰,仅在两个储层位置处有很强的频散异常。
图3为实施例工区1反演得到的新频散因子剖面。由测井曲线可知,在此处具有两套储层,而由预测的剖面结果可知,新的频散因子能够很好地识别出这两套储层,具有较高的分辨率,并且横向连续性很好,压制了无关的背景异常干扰。
上述实施例仅用于说明本发明,其中方法的各实施步骤等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。
Claims (4)
1.一种方位频变流体因子识别储层的方法,其特征在于采用以下具体步骤:
⑴对叠前角道集数据进行时频变换,并构建与新频散流体因子相关的目标函数;
r=eD
式中,r为m个频率下反射系数的列向量,共m×n行,e为m个频率下反射系数公式中呈m×n行、5列排列的系数矩阵,DM、Dμ、和分别为纵波模量、横波模量、法向弱度、新流体因子和切向弱度的频散程度,M和μ分别表示纵波模量和横波模量,δN和δT分别为岩石的法向和切向弱度,ψn为与裂缝弱度和松弛时间的平方相关的参数,其中,松弛时间受流体粘度、岩石渗透率、裂隙纵横比等参数的影响,该参数即为构建的一种新的流体识别因子,表达式为:ψn=δNΓ2,
式中:
式中,ηf为流体的粘度,k为岩石的渗透率,Kf为流体的体积模量,a为裂隙纵横比,Γ为松弛时间;
⑵选取参考频率ω0,利用最小二乘求解该参考频率下五类模量对应的频散结果:P、S、N、K、L,并将该结果带入向量r中,利用下式进行五类频散因子的同步反演:
D=(eTe)-1eTr
⑷重复上述步骤⑴至⑶,直到处理完整个工区的所有叠前地震角道集,即可获得整个工区的三维储层预测结果,该预测结果可为储层的高精度预测提供支撑。
3.根据权利要求1所述的一种方位频变流体因子识别储层的方法,其特征在于:构建了与流体粘度、岩石渗透率、裂缝纵横比等参数有关的新流体因子ψn,该因子可直接反映储层的特征。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110779032.6A CN113534241B (zh) | 2021-07-09 | 2021-07-09 | 一种方位频变流体因子识别储层的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110779032.6A CN113534241B (zh) | 2021-07-09 | 2021-07-09 | 一种方位频变流体因子识别储层的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113534241A CN113534241A (zh) | 2021-10-22 |
CN113534241B true CN113534241B (zh) | 2022-04-29 |
Family
ID=78098266
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110779032.6A Active CN113534241B (zh) | 2021-07-09 | 2021-07-09 | 一种方位频变流体因子识别储层的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113534241B (zh) |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6614717B1 (en) * | 2001-08-30 | 2003-09-02 | Nonlinear Seismic Imaging, Inc. | Time-lapse seismic using nonlinear seismic imaging |
CN103792578B (zh) * | 2012-10-26 | 2016-11-09 | 中国石油化工股份有限公司 | 一种频散avo相关属性反演的流体识别方法 |
CN105044780B (zh) * | 2015-08-26 | 2017-06-23 | 成都理工大学 | 依赖频率的avo频散参数反演中敏感参考频率的交会确定方法 |
CN108572389B (zh) * | 2017-03-14 | 2019-09-24 | 中国石油化工股份有限公司 | 频变粘弹性流体因子叠前地震反演方法 |
CN108897043B (zh) * | 2018-09-03 | 2020-09-04 | 中国石油天然气股份有限公司 | 基于弱度参数检测裂缝型储层流体的方法及装置 |
CN109507729A (zh) * | 2018-11-14 | 2019-03-22 | 成都理工大学 | 依赖频率的avo最优频散流体因子的反演方法 |
CN113126148A (zh) * | 2019-12-30 | 2021-07-16 | 中国石油大学(北京) | 一种裂缝型储层的流体识别方法、装置及*** |
-
2021
- 2021-07-09 CN CN202110779032.6A patent/CN113534241B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113534241A (zh) | 2021-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Daley et al. | Field testing of modular borehole monitoring with simultaneous distributed acoustic sensing and geophone vertical seismic profiles at Citronelle, Alabama | |
Dvorkin et al. | Dynamic poroelasticity: A unified model with the squirt and the Biot mechanisms | |
Gaiser | Multicomponent VP/VS correlation analysis | |
Yu et al. | On the characteristics of nonlinear soil response | |
US20190330981A1 (en) | Estimation of horizontal stresses and nonlinear constants in anisotropic formations such as interbedded carbonate layers in organic-shale reservoirs | |
CN110703354B (zh) | 一种基于井震结合的岩性边界识别方法 | |
CN111399056B (zh) | 一种基于划分方位角滤波预测裂缝强度的方法 | |
CN104237936B (zh) | 一种油气检测的频变反演方法 | |
Iyare et al. | Brittleness evaluation of Naparima Hill mudstones | |
US20120092958A1 (en) | Estimation of anisotropy from compressional waves from array sonic waveforms in well logging | |
Sleefe et al. | The use of broadband microseisms for hydraulic-fracture mapping | |
CN114137616A (zh) | 一种利用量子力学原理进行储层含气性检测的方法 | |
CN113534241B (zh) | 一种方位频变流体因子识别储层的方法 | |
CN115061202B (zh) | 一种页岩含气性地震储层直接检测方法 | |
CN111077574A (zh) | 一种确定地层弹性参数的方法、装置及*** | |
Merzlikin et al. | Unconventional reservoir characterization using azimuthal seismic diffraction imaging | |
Zhang et al. | Fluid identification based on P-wave anisotropy dispersion gradient inversion for fractured reservoirs | |
CN113640876A (zh) | 一种利用混沌体属性对陷落柱进行精细识别的方法 | |
Jamali et al. | Role of stress regime in azimuthal anisotropy of poroelastic media | |
CN113126148A (zh) | 一种裂缝型储层的流体识别方法、装置及*** | |
Li et al. | Characteristics of elastic wave dispersion and attenuation induced by microcracks in complex anisotropic media | |
Li et al. | Stress-sensitivity analysis of rock mechanical and petrophysical properties of fractured tight sandstone under true-triaxial stresses in the laboratory | |
CN108549103A (zh) | 饱和多孔介质叠后地震烃类检测方法 | |
Altowairqi et al. | A quantitative application of seismic inversion and multi-attribute analysis based on rock physics linear relationships to identify high total organic carbon shale-a case study from the perth basin, western australia | |
Yang et al. | In situ stress field detection of stress-induced strong anisotropy media based on Mohr circle theory |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |