CN111856564B - 一种储层流体检测方法 - Google Patents
一种储层流体检测方法 Download PDFInfo
- Publication number
- CN111856564B CN111856564B CN201910337150.4A CN201910337150A CN111856564B CN 111856564 B CN111856564 B CN 111856564B CN 201910337150 A CN201910337150 A CN 201910337150A CN 111856564 B CN111856564 B CN 111856564B
- Authority
- CN
- China
- Prior art keywords
- waves
- seismic
- fast
- transverse
- slow
- 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 144
- 238000001514 detection method Methods 0.000 title claims abstract description 26
- 239000011435 rock Substances 0.000 claims abstract description 28
- 238000000034 method Methods 0.000 claims description 37
- 238000004458 analytical method Methods 0.000 claims description 20
- 230000008859 change Effects 0.000 claims description 16
- 238000001228 spectrum Methods 0.000 claims description 12
- 239000006185 dispersion Substances 0.000 claims description 8
- 230000010287 polarization Effects 0.000 claims description 8
- 230000000694 effects Effects 0.000 claims description 5
- 238000000926 separation method Methods 0.000 claims description 5
- 238000012360 testing method Methods 0.000 claims description 5
- 239000013598 vector Substances 0.000 claims description 4
- 238000012800 visualization Methods 0.000 claims description 4
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 3
- 238000011161 development Methods 0.000 abstract description 9
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 8
- 230000004044 response Effects 0.000 description 7
- 238000011160 research Methods 0.000 description 6
- 230000035945 sensitivity Effects 0.000 description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 6
- 238000011156 evaluation Methods 0.000 description 4
- 239000004215 Carbon black (E152) Substances 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 3
- 229930195733 hydrocarbon Natural products 0.000 description 3
- 150000002430 hydrocarbons Chemical class 0.000 description 3
- 230000007480 spreading Effects 0.000 description 3
- 230000003044 adaptive effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 229920006395 saturated elastomer Polymers 0.000 description 2
- 238000010206 sensitivity analysis Methods 0.000 description 2
- 241001442234 Cosa Species 0.000 description 1
- 244000089409 Erythrina poeppigiana Species 0.000 description 1
- 235000019738 Limestone Nutrition 0.000 description 1
- 235000009776 Rathbunia alamosensis Nutrition 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002637 fluid replacement therapy Methods 0.000 description 1
- 239000006028 limestone Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013031 physical testing Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 230000009466 transformation 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/30—Analysis
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
技术领域
本发明涉及一种储层流体检测方法,属于石油勘探开发领域。
背景技术
现阶段储层流体检测的方法,主要包括以下内容:(1)根据不同流体性质下储层岩石物理性质的差异性分析研究,分析储层环境下流体的物理性质,结合流体的分布特征及其与孔隙结构的联系,估算含流体岩石的波响应特征,并与实际问题相联系;(2)在储层地震资料叠前和叠后反演的基础上,进行储层流体检测;(3)发现并分析对储层流体敏感的地震属性,进行储层流体检测。
申请号为CN201210315777.8的中国专利文件公开了一种储层流体检测方法和储层流体检测装置,具体涉及一种储层流体检测方法,该方法包括:角道集获取步骤,对地震数据进行处理,提取角道集;参数获取步骤,根据所获得的角道集,获得纵波速度Vp、横波速度Vs和密度p;流体因子确定步骤,根据计算出的纵波速度Vp、横波速度Vs和密度p确定指数泊松比流体因子,其对数是岩石或地层的泊松比的函数,以及油气层确定步骤,根据所述指数泊松比流体因子,确定是否存在油气层。
申请号为CN201410593709.7的中国专利文件公开了一种基于等效流体声波阻抗的流体识别方法,具体涉及一种基于等效流体声波阻抗的流体识别方法,其先通过测井方法或岩石物理测试方法获得储层的数据集合,再构建饱和流体的岩石密度,然后利用反演得到饱和流体岩石的声波阻抗,将声波阻抗分解成骨架等效声波阻抗与等效流体声波阻抗之和,利用反演的声波阻抗及骨架等效声波阻抗计算等效流体声波阻抗;等效流体声波阻抗的数据集合用石油物探专业的地震显示软件进行图形显示,并依据显示图形进行流体识别。
申请号为CN201310392124.4的中国专利文件公开了一种基于信号低频强异常自适应提取的地震资料流体识别方法,具体涉及一种基于信号低频强异常自适应提取的地震资料流体识别方法,该方法包括:(1)查看钻井岩心资料确定含流体储层的地下深度位置,并对井旁地震道的地震信号做谱分解,然后观察含流体储层相比围岩是否具有低频强振幅的地震反射特征,如果具有该特征,则转入步骤2;(2)对一道地震信号做基于低频强振幅异常的谱分解得到该道地震信号的时频谱;(3)获取步骤2得到的时频谱中每个时刻的瞬时最大能量;(4)对下一道地震信号重复步骤2和步骤3,直到所有地震道部完成为止;(5)利用所述低频强振幅的地震反射特征进行地震资料流体识别。
油气藏评价与开发(2012-2)文章《岩石物理参数的流体敏感性分析》针对不同岩石物理参数对储层流体的敏感性各异,提出了岩石物理流体敏感性参数的构建方法,并对不同的岩石物理参数及组合进行流体敏感性的分析。在岩石物理参数分类基础上,构建由剪性参数和体性参数组成的组合参数,将体性参数λ表示为干岩石参量和流体参量两部分的综合贡献,参考实验数据,选择最佳系数c值,使构建的组合参数最大程度地反映流体的贡献。在东营北带,通过实验室中深层岩石样品的测定和计算,得到基本岩石物理参数,按照流体敏感参数的定义对该区的不同岩石物理参数及组合进行流体敏感性分析,结果表明,构建的组合参数对储层流体的敏感程度更高,这对寻求敏感性较好的岩石物理参数进行储层流体识别具有指示意义。
地球物理学进展(2013-6)文章《阿姆河右岸麦捷让地区碳酸盐岩储层流体检测研究》从岩石物理模型入手,分析石灰岩纵横波速度、弹性参数与岩石基本性质之间的定量关系,寻找对流体敏感的弹性属性,认为纵横波速度比能够较好地识别气水分布。通过测井弹性参数交会和流体替换,以实际资料验证岩石物理分析结果并选择气藏敏感参数,建立适合该地区的气藏预测模版。利用叠前时间偏移得到的共反射点道集,开展三维叠前同时反演,利用反演结果进行储层流体识别。
石油物探(2013-5)文章《基于叠前反演的缝洞型碳酸盐岩流体检测及评价技术研究》针对流体检测特别是缝洞型碳酸盐岩储集体含流体性质检测的地球物理难题,基于改进的适用于碳酸盐岩岩石物理分析的Xu&White岩石物理模型,建立起碳酸盐岩纵、横波速度预测方法技术流程,对弹性参数检测流体的有效性进行分析。在碳酸盐岩岩石物理分析基础上,开展叠前反演弹性参数流体检测和叠前流体反演评价技术研究,探索实现塔河油田主体目标区缝洞型碳酸盐岩储层流体检测及预测结果的概率评价分析。
现阶段储层流体检测依然是世界性难题,在储层进行开发时,储层流体的分布预测和流体性质检测是储层成功开发的重要保证。然而,储层流体检测对储层的开发意义重大却仍面临巨大的困难。首先,由于储层流体变化对地震响应的影响相对于储层变化对地震响应的影响更小也更复杂,因此,很难确定储层不同流体性质下地震响应的特点;其次,地下各种温度、(静岩)压力、孔隙流体压力、不同流体饱和度,都会对储层的岩石物理特性产生影响,强烈的非均质性和各向异性容易造成波场特征复杂,流体对地震响应的影响难以和其它因素变化造成的对地震响应的影响区分开来;再次,油和水对储层P波和S波的影响基本相同,波阻抗的变化较小,这些都会制约着储层流体检测的效果;最后,对储层流体识别技术的研究特别是针对裂缝型储层的流体识别研究依然缺乏系统性。此外,现有的岩石物理显示,随频率变化,油和水造成的剪切S波分裂。流体饱和度对P波和快S波速度的影响小于1%,而流体性质对慢(似)S波速度的影响大于4%,而这些特征在现在的碳酸盐岩储层流体检测技术中还没有很好地加以应用。
发明内容
针对上述问题,本发明的目的是提供一种储层流体检测方法,其能够进行储层流体的分布预测和流体性质检测,获取各种流体空间上的分布情况,进而能够实现油气的高效安全开发。
为实现上述目的,本发明采取以下技术方案:一种储层流体检测方法,包括以下步骤:根据地震波中横波地震数据,获取快横波和慢横波地震数据;根据地震波中纵波和快、慢横波地震数据,获取地震波中纵波和快、慢横波的分频属性和频率衰减属性;根据地震波中纵波和快、慢横波地震数据,获取地震波中纵波和快、慢横波的速度的空间展布;根据地震波中纵波和快、慢横波在各向异性介质储层中传播,确定地震波中纵波和快、慢横波的岩石物理特性的差异;根据地震波中纵波和快、慢横波的岩石物理特性的差异,结合纵波和快、慢横波的分频属性和频率衰减属性以及速度的空间展布进行储层流体检测并确定储层流体性质。
在一个优选的实施例中,根据储层流体性质,获取储层中各流体空间上的展布状况。
在一个优选的实施例中,根据横波地震数据,结合快横波和慢横波在空间偏振的差异进行地震波中快横波和慢横波的分离,获取快横波和慢横波地震数据;
确定快、慢横波的径向分量和水平分量地震记录所依据的公式为:
R(t)=S1(t)cosα+S2(t)sinα (1)
T(t)=S1(t)sinα-S2(t)cosα (2)
式中,R为快横波和慢横波的径向分量地震记录,T为快横波和慢横波的横向分量地震记录,S1为快横波地震记录,S2(t)为慢横波地震记录,α为S1与R之间的夹角。
在一个优选的实施例中,采用互相关法和正交基旋转法分别进行地震波中快横波和慢横波的分离。
在一个优选的实施例中,采用互相关法分离地震波中的快横波和慢横波所依据的公式为:
S1(t)=R(t)cosa+T(t)sinα (3)
S2(t)=R(t)sina+T(t)cosα (4)
式中,对α进行角度扫描,分别计算对应角度的快横波地震记录和慢横波地震记录,并计算快横波地震记录和慢横波地震记录之间的相关值,α等于裂缝走向与测线方向的夹角时,快横波地震记录与慢横波地震记录的互相关值达到极大值,对快横波和慢横波的径向分量和横向分量地震记录进行旋转分离出快横波和慢横波;
基于快横波和慢横波为正交偏振时,欧氏空间中两个相互正交的向量的内积为零,得到采用正交基旋转法分离地震波中快横波和慢横波所依据的公式为:
在一个优选的实施例中,获取地震波中纵波和快、慢横波的分频属性,包括利用时频分析算法将地震波中纵波和快、慢横波的地震数据由时间域变换到频率域,计算每个地震道所有时间点的频谱,获取共频率的地震数据体,并利用可视化对各分频地震数据体进行对比、分析和解释。
在一个优选的实施例中,获取地震波中纵波和快、慢横波的频率衰减属性,包括对每个地震道的地震记录进行时频分析,在时频剖面上将每个时间样点处局部频率中的能量频率的极大值作为初始的衰减频率,再计算总能量中间隔两点的斜线斜率,并确定该斜率为地震波中纵波和快、慢横波的频率衰减属性。
在一个优选的实施例中,根据纵波和快、慢横波地震数据,进行地震波的速度分析或地震速度反演,获取地震波中纵波和快、慢横波速度的空间展布状况。
在一个优选的实施例中,地震波中纵波和快、慢横波的岩石物理特性的差异包括储层中流体有和无以及流体性质不同时,纵波和快、慢横波的分频属性和频率衰减属性以及速度的差异。
在一个优选的实施例中,储层中流体有和无以及流体性质不同对纵波和快、慢横波的分频属性和频率衰减属性以及速度的影响包括储层流体的黏度增加,纵波和慢横波的频散曲线随频率变化分别形成一个对黏度敏感的频带,且纵波和慢横波的频散曲线随频率变化呈相反趋势移动,快横波速度不随频率变化而变化,纵波和慢横波速度随频率增加而增加,且频率增加对慢横波速度的影响程度大于频率增加对纵波速度的影响程度。
本发明由于采取以上技术方案,其具有以下优点:1、本发明能够简捷高效地确定储层中流体性质,进而能够准确地获取各种流体空间上的分布情况,根据储层中流体性质和各种流体空间上的分布情况,实现油气的高效安全开发。2、本发明能够利用横波地震数据进行快横波和慢横波的分离。3、本发明能够快速且准确地确定储层中流体有和无以及流体性质不同分别对纵波和快、慢横波的分频属性和频率衰减属性以及速度的影响。
本发明的其它特征和优点将在随后的说明书中阐述,并且,部分的从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
附图说明
为了更清楚的说明本发明实施例中的技术方案,下面将对实施例描述中所需要的附图做简单的介绍:
图1是本发明的一种储层流体检测方法的一个具体实施例的流程示意图;
图2是本发明的一个具体实施例的快、慢横波的径向和横向分量以及极化方向的示意图;
图3是本发明的一个具体实施例的快、慢横波径向分量时差的示意图;
图4是本发明的一个具体实施例的快、慢横波横向分量时差的示意图;
图5是本发明的一个具体实施例的储层中三个各向异性层分别重建的快、慢横波的地震记录的示意图;
图6是本发明的一个具体实施例的地震数据进行分频属性体生成的示意图;
图7是本发明的一个具体实施例的纵波速度随频率和流体性质变化而变化的示意图;
图8是本发明的一个具体实施例的快横波和慢横波的速度随频率和流体性质变化而变化的示意图;
图9是本发明的一个具体实施例的横波分裂随频率和流体性质变化而变化的示意图。
具体实施方式
以下将结合附图及实施例来详细说明本发明的实施方式,借此对本发明如何应用技术手段来解决技术问题,并达成技术效果的实现过程能充分理解并据以实施。需要说明的是,只要不构成冲突,本发明中的各个实施例以及各实施例中的各个特征可以相互结合,所形成的技术方案均在本发明的保护范围之内。
如图1所示,本发明提出的储层流体检测方法,包括以下步骤:
(1)根据横波地震数据,进行地震波中快横波和慢横波(包括慢横波和似慢横波)的分离,获取地震波中快横波和慢波地震数据
首先,获取地震数据,地震数据包括纵波和横波地震数据。接着,根据横波地震数据,结合快横波和慢横波在空间偏振的差异进行地震波中快横波和慢横波的分离,获取地震波中快横波和慢横波地震数据。
具体的,根据储层中流体的存在和流体性质的不同,地震波中P-SV横波穿越储层裂隙介质时会分裂成平行裂隙走向的快横波和垂直裂隙走向的慢横波,在某一时刻,快横波和慢横波的偏振在空间具有不同的组合形式。即地震波中横波在各向异性介质储层中会产生横波分裂现象,分裂后的快横波和慢横波在检波器的x和Y分量地震记录均能够产生能量投影(如图2所示),x和Y分量地震记录中分别包括快横波地震记录和慢横波地震记录。进一步地,根据快横波和慢横波存在的速度上的差异,进行地震波中快横波和慢横波的分离,获取地震波中快横波和慢横波地震数据。
优选的,检波器为三分量检波器。
根据投影叠加原理,确定快横波和慢横波的径向分量和水平分量地震记录所依据的公式为:
R(t)=S1(t)cosα+S2(t)sinα (1)
T(t)=S1(t)sinα-S2(t)cosα (2)
式中,R为快横波和慢横波的径向分量地震记录,T为快横波和慢横波的横向分量地震记录,S1为快横波地震记录,S2(t)为慢横波地震记录,α为S1与R之间的夹角。
进一步地,快横波和慢横波各向异性的大小可以通过两方面得到:一是基于运动学的快横波和慢横波旅行时差得到(如图3、图4所示);二是基于动力学的快横波和慢横波振幅差异得到。
进一步地,根据横波地震数据,基于快横波和慢横波在空间偏振的差异进行地震波中快横波和慢横波的分离时,可以采用互相关法和正交基旋转法分别进行地震波中快横波和慢横波的分离。
具体的,由于分离后的快横波和慢横波都是由一个横波分离出来的,快横波和慢横波仅存在极化方向和传播速度的不同,且快横波和慢横波的波形相似,因此,可以采用互相关法分离地震波中的快、慢横波。
采用互相关法分离地震波中的快、慢横波所依据的公式为:
S1(t)=R(t)cosα+T(t)sinα (3)
S2(t)=R(t)sinα+T(t)cosα (4)
式中,对α进行角度扫描,分别计算对应角度的快横波地震记录和慢横波地震记录,并计算快横波地震记录和慢横波地震记录之间的相关值;当α等于裂缝走向与测线方向的夹角时,快横波地震记录与慢横波地震记录的互相关值最大,确定此时的α角度值,接着根据公式(1)和(2)对快、慢横波的径向分量和横向分量地震记录进行旋转分离出快、慢横波;其中,公式(3)和(4)由公式(1)和(2)反解得到。
具体的,基于快横波和慢横波为正交偏振时,欧氏空间中两个相互正交的向量的内积为零,得到采用正交基旋转法(OBR)分离地震波中的快、慢横波所依据的公式为:
此时,欧氏空间中两个相互正交的向量的内积为零:
进而确定:
优选的,当数据信噪比较高时,采用正交基旋转法进行地震波中快横波和慢横波的分离。正交基旋转法具有可行性和实用性,准确度高的特点。
在一个具体实施例中,获取储层中三个各向异性层分别重建的快横波和慢横波的地震记录(如图5所示)。
(2)根据地震波中纵波和快、慢横波(快横波和慢横波)地震数据,获取地震波中纵波和快、慢横波的分频属性和频率衰减属性
分频属性是地震数据基于频谱分解得到的一种属性,通过数学变换,将时域信号转换成时频域的信息。
获取地震波中纵波和快、慢横波的分频属性,包括以下内容:利用时频分析算法将地震波中纵波和快、慢横波的地震数据由时间域变换到频率域,计算每个地震道所有时间点的频谱;接着,按照频率重新排列生成共频率的地震数据体;再利用可视化对各分频地震数据体进行对比、分析和解释。具体的,沿滑动时窗对整个地震波中纵波和快、慢横波的地震数据体进行时频分析,计算每个样本的振幅谱和相位谱,然后生成共频率的地震数据体。共频率地震数据体纵坐标是时间,与常规地震数据相同,每个共频率地震数据体中仅包括一个频率,可以从时间、空间、频率可视化共频率地震数据体进行解释。
在一个具体实施例中,如图6所示,确定地震数据的分频属性体,包括:
①获取3-D地震数据体,并进行时频分析;
②计算3-D地震数据体中每个样本的振幅谱和相位谱,生成4-D时频体数据;
③离散频率能量体,获得多个共频率的地震数据体。
获取地震波中纵波和快、慢横波的频率衰减属性,包括以下内容:首先对每个地震道的地震记录做时频分析,包括对每个地震道的地震记录进行时频分析,在时频剖面上将每个时间样点处局部频率中的能量频率的极大值作为初始的衰减频率,再计算总能量65%和85%两点之间的斜线斜率,并确定该斜率为频率衰减属性。由于频率衰减属性的计算方法在地震数据信噪比比较高,因此其能够较好地反映频率的衰减特征。
(3)根据地震波纵波和快、慢横波地震数据,获取地震波中纵波和快、慢横波的速度的空间展布
根据纵波和快、慢横波地震数据,进行地震波的速度分析或地震速度反演,获取地震波中纵波和快、慢横波速度的空间展布状况。
具体的,对地震波中纵波和快、慢横波的地震数据进行速度分析或地震速度反演,得到三维数据体,在此基础上,能够确定纵波和快、慢横波速度的空间展布情况。其中,
速度分析方法包括通过速度扫描法、常速叠加(CVS)法和速度谱法,获得叠加速度,并通过Dix公式转化为深度域层速度。
地震速度反演方法包括通过地表或井间观测到的地震运动学信息(如走时和反射路径等)或者动力学信息(如振幅和波形等)的分析,建立模型参数与观测数据之间的泛函关系,并利用最佳匹配恢复储层地下或井间速度参数信息。
(4)根据地震波中纵波和快、慢横波在各向异性介质储层中传播,确定地震波中纵波和快、慢横波的岩石物理特性的差异
根据地震波中纵波和快、慢横波在各向异性介质储层中传播时,流体有和无以及流体性质不同,表现出不同的岩石物理特性(具体指流体有和无以及流体性质不同对分频属性和频率衰减属性以及速度的影响),进一步确定流体有和无以及流体性质不同时,地震波中纵波和快、慢横波的岩石物理特性的差异。
地震波中纵波和快、慢横波的岩石物理特性差异包括储层中流体有和无以及流体性质不同而导致的纵波和快、慢横波的分频属性和频率衰减属性以及速度的差异。
储层中流体有和无以及流体性质不同对地震波中纵波和快、慢横波的分频属性和频率衰减属性以及的速度影响包括储层流体的黏度增加(气-水-油黏度增加),即流体性质变化,纵波的频散曲线向低频(非方位)移动形成一个对黏度敏感的频带(如图7所示),慢横波的频散曲线与纵波的频散曲线的移动方向相反,慢横波的频散曲线向高频(非方位)移动形成一个对黏度敏感的频带(如图8所示);同时,快横波速度不随频率变化而发生变化,即流体性质变化,对快横波速度影响不大,纵波和慢横波速度随频率增加而增加,且频率增加对慢横波速度的影响程度大于频率增加对纵波速度的影响程度(如图7、图8所示)。其中,储层流体性质变化包括储层流体中气、水和油三种流体形式的变化,气-水-油的黏度依次增加。
在一个具体实施例中,频率增加对慢横波速度的影响大于4%(如图7、图8所示)。
在一个具体的实施例中,储层液体中油和水造成的剪切横波分裂随频率和流体性质变化而变化(如图9所示)。
(5)根据地震波中纵波和快、慢横波的岩石物理特性的差异,结合纵波和快、慢横波的分频属性和频率衰减属性以及速度的空间展布进行储层流体检测并确定储层流体性质
根据储层中流体有和无以及流体性质不同时,所致的地震波中纵波和快、慢横波的岩石物理特性的差异,结合纵波和快、慢横波的分频属性和频率衰减属性进行储层流体检测并确定储层流体性质。
(6)进一步获取储层中各流体空间上的展布状况
根据储层中流体性质,获取储层中各流体空间上的展布状况。
进一步地,根据储层中各流体空间上的展布状况,实现油气的高效安全开发。
虽然本发明所公开的实施方式如上,但所述的内容只是为了便于理解本发明而采用的实施方式,并非用以限定本发明。任何本发明所属技术领域内的技术人员,在不脱离本发明所公开的精神和范围的前提下,可以在实施的形式上及细节上作任何的修改与变化,但本发明的专利保护范围,仍须以所附的权利要求书所界定的范围为准。
Claims (10)
1.一种储层流体检测方法,其特征在于,包括以下步骤:
根据地震波中横波地震数据,获取快横波和慢横波地震数据;
根据地震波中纵波和快、慢横波地震数据,获取地震波中纵波和快、慢横波的分频属性和频率衰减属性;
根据地震波中纵波和快、慢横波地震数据,获取地震波中纵波和快、慢横波的速度;
根据地震波中纵波和快、慢横波在各向异性介质储层中传播,确定地震波中纵波和快、慢横波的岩石物理特性的差异;
根据地震波中纵波和快、慢横波的岩石物理特性的差异,结合纵波和快、慢横波的分频属性和频率衰减属性以及速度进行储层流体检测并确定储层流体性质。
2.根据权利要求1所述的一种储层流体检测方法,其特征在于,根据储层流体性质,获取储层中各流体空间上的展布状况。
3.根据权利要求1或2所述的一种储层流体检测方法,其特征在于,根据横波地震数据,结合快横波和慢横波在空间偏振的差异进行地震波中快横波和慢横波的分离,获取快横波和慢横波地震数据;
确定快、慢横波的径向分量和水平分量地震记录所依据的公式为:
R(t)=S1(t)cosα+S2(t)sinα (1)
T(t)=S1(t)sinα-S2(t)cosα (2)
式中,R(t)为快横波和慢横波的径向分量地震记录,T(t)为快横波和慢横波的横向分量地震记录,S1(t)为快横波地震记录,S2(t)为慢横波地震记录,α为S1(t)与R(t)之间的夹角。
4.根据权利要求3所述的一种储层流体检测方法,其特征在于,采用互相关法和正交基旋转法分别进行地震波中快横波和慢横波的分离。
5.根据权利要求4所述的一种储层流体检测方法,其特征在于,采用互相关法分离地震波中的快横波和慢横波所依据的公式为:
S1(t)=R(t)cosα+T(t)sinα (3)
S2(t)=R(t)sinα+T(t)cosα (4)
式中,对α进行角度扫描,分别计算对应角度的快横波地震记录和慢横波地震记录,并计算快横波地震记录和慢横波地震记录之间的相关值,α等于裂缝走向与测线方向的夹角时,快横波地震记录与慢横波地震记录的互相关值达到极大值,对快横波和慢横波的径向分量和横向分量地震记录进行旋转分离出快横波和慢横波;
基于快横波和慢横波为正交偏振时,欧氏空间中两个相互正交的向量的内积为零,得到采用正交基旋转法分离地震波中快横波和慢横波所依据的公式为:
6.根据权利要求1所述的一种储层流体检测方法,其特征在于,获取地震波中纵波和快、慢横波的分频属性,包括利用时频分析算法将地震波中纵波和快、慢横波的地震数据由时间域变换到频率域,计算每个地震道所有时间点的频谱,获取共频率的地震数据体,并利用可视化对各分频地震数据体进行对比、分析和解释。
7.根据权利要求1所述的一种储层流体检测方法,其特征在于,获取地震波中纵波和快、慢横波的频率衰减属性,包括对每个地震道的地震记录进行时频分析,在时频剖面上将每个时间样点处局部频率中的能量频率的极大值作为初始的衰减频率,再计算总能量中间隔两点的斜线斜率,并确定该斜率为地震波中纵波和快、慢横波的频率衰减属性。
8.根据权利要求1所述的一种储层流体检测方法,其特征在于,根据纵波和快、慢横波地震数据,进行地震波的速度分析或地震速度反演,获取地震波中纵波和快、慢横波速度的空间展布状况。
9.根据权利要求1所述的一种储层流体检测方法,其特征在于,地震波中纵波和快、慢横波的岩石物理特性的差异包括储层中流体有和无以及流体性质不同时,纵波和快、慢横波的分频属性和频率衰减属性以及速度的差异。
10.根据权利要求1所述的一种储层流体检测方法,其特征在于,储层中流体有和无以及流体性质不同对纵波和快、慢横波的分频属性和频率衰减属性以及速度的影响包括储层流体的黏度增加,纵波和慢横波的频散曲线随频率变化分别形成一个对黏度敏感的频带,且纵波和慢横波的频散曲线随频率变化呈相反趋势移动,快横波速度不随频率变化而变化,纵波和慢横波速度随频率增加而增加,且频率增加对慢横波速度的影响程度大于频率增加对纵波速度的影响程度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910337150.4A CN111856564B (zh) | 2019-04-24 | 2019-04-24 | 一种储层流体检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910337150.4A CN111856564B (zh) | 2019-04-24 | 2019-04-24 | 一种储层流体检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111856564A CN111856564A (zh) | 2020-10-30 |
CN111856564B true CN111856564B (zh) | 2022-02-01 |
Family
ID=72952323
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910337150.4A Active CN111856564B (zh) | 2019-04-24 | 2019-04-24 | 一种储层流体检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111856564B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5835452A (en) * | 1995-10-06 | 1998-11-10 | Amoco Corporation | Reflected shear wave seismic processes |
CN102053259A (zh) * | 2009-10-28 | 2011-05-11 | 中国石油化工股份有限公司 | 一种油气藏识别方法 |
CN102879817A (zh) * | 2012-07-23 | 2013-01-16 | 北京多分量地震技术研究院 | 基于地面地震数据获取地下裂缝信息的控制方法 |
CN106353802A (zh) * | 2016-08-19 | 2017-01-25 | 中国石油化工股份有限公司石油勘探开发研究院 | 转换波信号的快慢波分离方法 |
CN107678061A (zh) * | 2017-10-24 | 2018-02-09 | 中国地质大学(北京) | 基于横波分裂振幅差属性的流体检测方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100238765A1 (en) * | 2009-03-20 | 2010-09-23 | Grechka Vladimir | Reservoir characterization from multicomponent microseismic data |
US20140334262A1 (en) * | 2013-04-22 | 2014-11-13 | Robert H. Brune | Method and Apparatus for Active Seismic Shear Wave Monitoring of Hydro-Fracturing of Oil and Gas Reservoirs Using Arrays of Multi-Component Sensors and Controlled Seismic Sources |
CN107797144A (zh) * | 2017-10-24 | 2018-03-13 | 中国地质大学(北京) | 基于横波分裂振幅比属性的流体检测方法 |
-
2019
- 2019-04-24 CN CN201910337150.4A patent/CN111856564B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5835452A (en) * | 1995-10-06 | 1998-11-10 | Amoco Corporation | Reflected shear wave seismic processes |
CN102053259A (zh) * | 2009-10-28 | 2011-05-11 | 中国石油化工股份有限公司 | 一种油气藏识别方法 |
CN102879817A (zh) * | 2012-07-23 | 2013-01-16 | 北京多分量地震技术研究院 | 基于地面地震数据获取地下裂缝信息的控制方法 |
CN106353802A (zh) * | 2016-08-19 | 2017-01-25 | 中国石油化工股份有限公司石油勘探开发研究院 | 转换波信号的快慢波分离方法 |
CN107678061A (zh) * | 2017-10-24 | 2018-02-09 | 中国地质大学(北京) | 基于横波分裂振幅差属性的流体检测方法 |
Non-Patent Citations (3)
Title |
---|
Shear Wave Splitting Analysis to Estimate Fracture Orientation and Frequency Dependent Anisotropy;Raoof GHOLAMI et al.;《Acta Geophysica》;20160229;第64卷(第1期);第76-100页 * |
云南高黎贡山典型板岩地震波各向异性研究;郭彬彬等;《地球物理学报》;20140331;第57卷(第3期);第837-846页 * |
地震资料中的多尺度断裂与流体响应特征:数值模拟分析;郭桂红等;《地球物理学报》;20130630;第56卷(第6期);第2002-2011页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111856564A (zh) | 2020-10-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Paillet et al. | Acoustic modes of propagation in the borehole and their relationship to rock properties | |
US5831935A (en) | Method for geophysical processing and interpretation using seismic trace difference for analysis and display | |
US8706420B2 (en) | Seismic fluid prediction via expanded AVO anomalies | |
US6789018B1 (en) | Mapping reservoir rocks using frequency spectral broadening and the presence of the slow-wave | |
US6631783B2 (en) | Mapping reservoir characteristics using earth's nonlinearity as a seismic attribute | |
Russell | Prestack seismic amplitude analysis: An integrated overview | |
Qi et al. | Sonic QP/QS ratio as diagnostic tool for shale gas saturation | |
Prioul et al. | Forward modeling of fracture-induced sonic anisotropy using a combination of borehole image and sonic logs | |
Naseer et al. | Porosity prediction of lower cretaceous unconventional resource play, south Indus Basin, Pakistan, using the seismic spectral decomposition technique | |
Germán Rubino et al. | Seismic characterization of thin beds containing patchy carbon dioxide-brine distributions: A study based on numerical simulations | |
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 | |
CN111399056A (zh) | 一种基于划分方位角滤波预测裂缝强度的方法 | |
Li et al. | Fractured formation evaluation by seismic attenuation derived from array acoustic log waves based on modified spectral ratio method and an extended Biot's poroelastic model | |
Burns et al. | Fracture properties from seismic scattering | |
Durrani et al. | Characterization of seismic anisotropy using azimuthal AVO analysis (AVAz)-An application case study in the deep and tight carbonate reservoirs from Potwar Basin onshore Pakistan | |
CN111856564B (zh) | 一种储层流体检测方法 | |
Wang et al. | Gas prediction using low-frequency components of variable-depth streamer seismic data applied to the deepwater area of the South China Sea | |
Luo et al. | Abnormal transmission attenuation and its impact on seismic-fracture prediction—A physical modeling study | |
CN110658553B (zh) | 一种储层流体性质检测的方法与系统 | |
Li et al. | Recent advances in application of AVO to carbonate reservoirs | |
Wang et al. | Fluids discrimination by ray-path elastic impedance inversion: A successful case from Sulige tight gas field | |
Gao et al. | Sparse spectral attributes to detect hydrocarbons: A case study from deep dolomite reservoirs | |
Strecker et al. | Seismic attribute analysis in hydrothermal dolomite, Devonian slave point formation, Northeast British Columbia, Canada | |
CN112764100B (zh) | 储层油气范围检测方法及装置 | |
Yuan et al. | Interpretation of fractures and joint inversion using multicomponent seismic data—Marcellus Shale example |
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 |