CN1291241C - 等速砂泥岩地层的储层预测油气的方法 - Google Patents

等速砂泥岩地层的储层预测油气的方法 Download PDF

Info

Publication number
CN1291241C
CN1291241C CN 200510041911 CN200510041911A CN1291241C CN 1291241 C CN1291241 C CN 1291241C CN 200510041911 CN200510041911 CN 200510041911 CN 200510041911 A CN200510041911 A CN 200510041911A CN 1291241 C CN1291241 C CN 1291241C
Authority
CN
China
Prior art keywords
amplitude
sandstone
acceleration
reservoir
well
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.)
Expired - Fee Related
Application number
CN 200510041911
Other languages
English (en)
Other versions
CN1664616A (zh
Inventor
蒲仁海
钟红利
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwest University
Original Assignee
Northwest University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northwest University filed Critical Northwest University
Priority to CN 200510041911 priority Critical patent/CN1291241C/zh
Publication of CN1664616A publication Critical patent/CN1664616A/zh
Application granted granted Critical
Publication of CN1291241C publication Critical patent/CN1291241C/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了用振幅加速度在等速砂泥岩地层区预测油气储层的方法。步骤为向计算机中加载三维地震保幅体和偏移体;用三维保幅资料计算三维地震数据体的振幅加速度体。用偏移三维地震体解释并逐线或逐道内插勘探目的层的各个反射层位;根据钻井标定,按照地震的最大纵向分辨率提取平均每半个周期对应地质层位的地震振幅加速度;统计计算各井、各砂组的砂岩发育指数、对应地震层位的产能,提取各井点对应地震层位的振幅加速度值。将振幅加速度与产能或砂岩发育指数进行二元拟合。具有快捷、有效、经济、实用、运算速度快,成本低特点,成功率达60-82%,广泛用于相应地质条件下的油气勘探,有很好的社会和经济效益。

Description

等速砂泥岩地层的储层预测油气的方法
技术领域
本发明涉及地层中石油勘探预测油气的方法,特别适用于用等速砂泥岩地层地区的储层预测,而对于这类地区其他地震储层预测方法效果很差或没有效果。
背景技术
现有技术中,当含油气层系的砂岩速度大于或小于泥岩速度时,沿层振幅平面图或反演波阻抗剖面图及平面图即可反映出沉积体系格架和砂岩储层的展布,但当砂岩与泥岩速度相等或接近时,这种方法并不奏效,因而它成为目前地震储层预测的难题之一。N.S.Neidel最早指出砂岩与泥岩的速度差异对地震储层预测效果的影响,并把砂岩速度小于泥岩速度归于第一类地层,此时适合于亮点和AVO等储层与油气层预测;把砂岩速度大于泥岩速度归于第三类地层,此时适合于简单的砂岩含量或岩性预测;把砂岩速度等于泥岩速度归于第二类地层,此时用振幅、速度、波阻抗、频率等属性均不能有效的识别和预测储层。砂岩与泥岩的速度差异取决于压实作用,刘震把砂泥岩速度相等的情况称为等速型压实模型。
储层预测的另一个难题是薄互层中砂岩的预测,提高地震采集主频和测井约束反演均可提高或改善对薄层的识别和预测能力,但前者在生产中并未广泛使用,后者则需均匀和密集的钻井分布(1)。
松辽盆地南部的白垩系青山口组储层和鄂尔多斯盆地北部的二叠系石盒子组储层属于具有以上特点的储层类型。储层的存在与否或砂岩含量的变化往往难以引起振幅或波阻抗差的变化,因而用一般地震属性提取和波阻抗反演方法等方法均不能有效地识别和预测油气储层。
发明内容
本发明的发明目的在于对于砂泥岩等速地层提供一种经济、实用、简便有效的薄互层中砂岩储层的识别和预测,提高探井和评价井的钻井成功率的等速砂泥岩地层的储层预测油气的方法。
实现发明目的的技术方案是这样解决的:等速砂泥岩地层的储层预测油气的方法,包括①计算研究区的地震振幅加速度体;②沿目的层提取显示振幅加速度平面图;③在该平面图上解释和圈定振幅加速度异常值范围作为有利油气储层平面展布范围,一个计算软件程序,其方法步骤为:
a、向计算机中加载三维地震保幅体和偏移体;
b、用三维保幅资料计算三维地震数据体的振幅加速度体,设地震反射强度为S,振幅加速度为AA,则
AA=d2(LnS)/dS(t)2,即对反射强度A取对数,再求其二阶导数;式中d为导数,Ln-为自然对数,t-为地震反射时间,S-为地震反射强度;
c、用偏移三维地震体解释并逐线或逐道内插勘探目的层的各个反射层位;
d、根据钻井标定,按照地震的最大纵向分辨率提取平均每半个周期对应地质层位的地震振幅、振幅加速度;采用均方根振幅或绝对值值振幅法即半个周期内各采样点振幅绝对值的和;
e、对以上属性的数值归一化,彩色显示其平面图,对于等速型砂泥岩地区,振幅加速度异常高或低值的平面展布范围即为储层发育区,含油层可能为异常低值,含气砂层可能为异常高值;
f、在测井基线校正和标准化的基础上,用砂岩的自然电位负异常幅度SP乘以厚度H作为砂岩的发育指数SI,统计计算各井、各砂组的砂岩发育指数SI;即SI=SP×H,SI取值一般介于0~1000之间,其大小取决于自然电位SP基线校正和标准化后的取值,即从泥岩基线到砂岩负异常幅度,单位为毫伏,厚度H指每层砂岩的厚度,单位为米;
g、统计各井在对应地震层位的产能单井产油或气量,提取各井点对应地震层位的振幅加速度值AA;
h、将振幅加速度与产能或砂岩发育指数进行二元拟合,对于等速型砂泥岩地区,砂岩发育指数或砂岩累加厚度或单井产量与振幅加速度呈幂函数拟合关系,相关系数R在0.6以上;
i、振幅加速度平面图与用测井相平面图、油层产能平面图进行比较,高产井区或厚砂层发育区往往具有高或低的振幅加速度异常,即在平面上振幅加速度总是能把储层区和非储层区有效地识别和区分;
j、对以上属性的数值归一化,彩色显示其平面图,对于等速型砂泥岩地区,振幅加速度异常高或低值的平面展布范围即为储层发育区,含油层可能为异常低值,含气砂层可能为异常高值;
k、用三维地震体的沿层振幅加速度平面图可识别和预测油气储层是通过试错法多次反复实验发现的,并在随后多次实践中被证明是有效的,试错法实验过程为:在叠后地震数据体中,通过计算处理可以获取50多个与时间、振幅、频率、相位、波形有关的属性,大多数地质条件下,简单的振幅、时间属性即可识别和预测储层,但对于等速型砂泥岩地层,常用的这些属性均不能有效预测储层,通过在多个油气区针对多个等速型砂泥岩层计算提取了各层的各种地震属性,再根据已有钻井揭示的储层发育与分布状况,分别与每个属性平面图对比、拟合,寻找与钻井结果一致的能反映储层分布的有效属性,振幅加速度能识别和预测等速型砂泥岩储层就是这样试验出来的;在腰英台油田为了获得能够有效反映砂岩反映程度地震属性,采用试错法分别提取9个同相轴对应层位的地震振幅、频率、波形分类、振幅加速度、振幅变化率、峰度、歪度、相干体属性,并与测井相平面图、测井砂岩发育指数平面图比较,发现除低的振幅加速度与砂岩发育指数高值分布的平面位置较吻合外,其他属性变化与砂岩的发育程度之间几乎没有关系,振幅加速度对砂岩反映效果最好,大多数情况下,低的加速度对应于厚、纯、粗的高砂岩指数分布区,在各层位的振幅加速度图上勾绘了砂岩发育指数自然电位负异常幅度乘以厚度等值线图成为砂岩储层预测重要手段之一;在东北腰英台油田实例中,振幅加速度一般介于10-20之间,砂岩发育指数200-400;根据振幅加速度和砂岩发育指数的关系解释出一个振幅加速度小于5的低值异常反映出一个朵状砂体,其砂岩发育指数大于600,获得了累厚超过20m的含油砂层和高产油气流;
l、在等速型砂泥岩地区,由于泥岩普遍含碳酸钙,含钙泥岩及钙质泥岩与砂岩的自然伽玛、电阻率、声波测井响应类似,自然电位测井对砂岩的响应较明显,它成为识别砂岩的最有效曲线,用自然电位幅度乘以厚度作为砂岩的发育指数,并在基线校正和标准化的基础上,统计计算各井、各砂组的砂岩发育指数,用于与地震属性的拟合,以实现最终用地震属性识别和预测砂岩发育区;
m、腰英台油田各井点的重点产油砂组的振幅加速度与砂岩发育指数之间略呈幂函数相关性,相关系数在0.6左右,尽管相关系数较小,但与其它地震属性比较,仍算是在储层识别、预测和沉积相平面展布分析方面一个较为有效的方法;
n、等速型砂泥岩地区各层位的振幅、波形、频率的确存在着平面变化,但他们与砂岩含量关系不大,经正演模型和地震属性与钻井的对比分析这些地震属性的变化反映了低速纯泥岩含量的变化,在钻井剖面图上,这些纯泥岩的占地层总厚度不到20%,而其他80%的岩性为声波时差变化很小的砂岩、粉砂岩及含钙泥岩,所以地震强振幅平面展布范围多指示了纯泥岩的范围,但弱振幅却没有指示出砂岩的富集区,波形参数对砂岩反映的灵敏度介于振幅加速度和振幅之间,具相同波形的部分区域指示了砂岩较发育的低加速度区,部分区域却指示了砂岩不发育的高加速度区;
o、在等速型砂泥岩地区的各层位的瞬时频率平面图上,频率的平面变化很小,各层位的频率总体上均有南高北低的趋势,似乎与南砂北泥的趋势吻合,正演模型也反映含砂较厚时,频率变高的特点,对腰英台油田的DB13~DB27井产油区来说,瞬时频率几乎无变化,因而对储层预测来讲没有意义;
p、属性试验结果均表明,对目的层砂岩储层的预测效果可以用好与不好进行评价;振幅加速度对储层预测效果好,其成功之处就在于识别出了砂体沉积相类型及展布的有利位置;除振幅加速度外,包括振幅、振幅变化率和波形在内的绝大部分属性预测效果均不好;
q、振幅加速度在应用于预测的鄂尔多斯盆地北部塔巴庙气田下石盒子组气层砂岩时,振幅加速度平面图指示了沉积相、砂岩储层、含气量三种信息,它指示该区发育一个由北流向南的一个曲流河道,在河道拐弯处发育了五个点坝砂体,气藏主要位于点坝之上,钻于点坝上的qc 3、qc 4、qc 10、qc 11、qc 16、qc 17、qc 20井每日产气量超过30000m3,钻于漫滩上qc 2、qc 13、qc 21、qc 22、qc 23井的日产气量小于10000m3,产气量与振幅加速度成幂函数递增关系,相关系数达0.82;
r、鄂尔多斯盆地北部塔巴庙气田下石盒子组砂岩层速度一般为4300~4700m/s,含气后砂岩速度4000m/s左右,泥岩不纯,多由粉砂质泥岩和含钙泥岩组成,速度一般4300~4700m/s,因而与砂岩之间无明显的波阻抗差,地层中夹少量薄层的低速纯泥岩或碳质泥岩,他们的含量变化是造成波阻抗差和引起振幅变化的主要原因,因此根据波阻抗反演和与振幅有关的属性变化很难预测砂岩储层,在探井qc3井盒2-3段获得日产53000m3的天然气后,根据反演和其它地震属性在该井南部2km的范围内预测的qc 21、qc 22、qc 23井均只得到了日产数千方的气,未达到工业产能,在盒2-3段的振幅加速度平面图上,具较厚砂岩储层和较高产能的钻井均位于高振幅加速度区域内,这种黄红色的高振幅加速度平面形态指示了沿曲流河道分布的5个点坝砂体,高产井位于点坝砂体内,低产井及干井则位于决口扇或漫滩上;
塔巴庙气田砂岩储层的自然电位负异常响应反映效果不明显,而含钙粉砂质泥岩和钙质泥岩与砂岩储层均具有较低的自然伽玛测井值和曲线形态,所以,采用自然电位和自然伽玛反演识别和预测储层也不十分理想;
s、在塔巴庙气田的振幅加速度与钻井产能关系图上,振幅加速度与产能呈正相关,相关系数达0.82,也就是说振幅加速度越大,则该地区产能就越大,等速型砂泥岩地层的泥岩往往含钙质,因而其自然伽玛、电阻率、声波时差测井响应与砂岩类似,二者不宜区分,测井约束反演的波阻抗、自然伽玛的资料也无法识别和预测储层,当自然电位曲线反映砂岩效果较好时,测井约束反演的自然电位也可较好的指示储层,这时自然电位反演和振幅加速度可以互相印证,提高预测的可信度。
本发明与现有技术相比,具有快捷、有效、经济、实用、检测方法简便、运算速度快,成本低,预测效果好的特点。运算提取地震振幅加速度参数,并与测井、岩芯和沉积相分析相结合,很大程度上改善了识别和预测此类特殊储层的有效性。实验表明振幅加速度与砂岩发育指数和产能均呈正相关,相关系数达0.82。采用本发明,可以解决其他方法所不能解决的等速型砂泥岩地区的薄互层中的砂岩油气储层的预测问题。将该方法用于适合于此中地质条件的油气勘探和开发中,定有很好的社会和经济效益。据发明者所知,属于等速型砂泥岩地质条件的油气区在国内外占有一定的比例,至少在国内发明者就遇到了四个地区:分别是松辽盆地南部青山口组油层、鄂尔多斯盆地北部下石盒子组储层(其中包括中石油的苏里格庙气田和中石化的塔巴庙气田)、江苏油用的下第三系油层、秦潼洼陷的下第三系油层等。在未使用本发明之前,这些地区地震储层预测效果较差,钻井失利多,勘探开发成本很高,使用本技术会取得显著效果。国内外其他油气区的这种地质条件也一定存在。在一个油气区,只要用目的层的综合测井曲线就很容易判断这个地区是否属于等速型砂泥岩地层这个地质条件。
附图说明
附图1为本发明实施例中腰英台地区砂岩储层预测等值线图;
附图2为塔巴庙气田某三维地震区块盒2-3振幅加速度平面图;
附图3振幅加速度与钻井产能关系图。
具体实施方式
附图为本发明的实施例子。
下面结合附图对本发明的内容作进一步说明:一种等速砂泥岩地层区储层预测油气的方法,包括①计算研究区的地震振幅加速度体;②沿目的层提取显示振幅加速度平面图;③在该平面图上解释和圈定振幅加速度异常值范围作为有利油气储层平面展布范围,一个计算软件程序,其方法步骤为:
a、向计算机中加载三维地震保幅体和偏移体;
b、用三维保幅资料计算三维地震数据体的振幅加速度体,设地震反射强度为S,振幅加速度为AA,则
AA=d2(LnS)/dS(t)2,即对反射强度S(t)取对数,再求其二阶导数;式中d为导数,Ln-代表自然对数,t-为地震反射时间,S-为地震反射强度;
c、用偏移三维地震体解释并逐线或逐道内插勘探目的层的各个反射层位;
d、根据钻井标定,按照地震的最大纵向分辨率提取平均每半个周期对应地质层位的地震振幅、振幅加速度;采用均方根或总振幅法效果较好;
e、对以上属性的数值归一化,彩色显示其平面图,对于等速型砂泥岩地区,振幅加速度异常高或低值的平面展布范围即为储层发育区,含油层为异常低值,含气砂层为异常高值;
f、在测井基线校正和标准化的基础上,用自然电位幅度SP乘以厚度H作为砂岩的发育指数SI,统计计算各井、各砂组的砂岩发育指数SI;SI取值一般介于0~1000之间,其大小取决于自然电位SP基线校正和标准化后的取值,即从泥岩基线到砂岩负异常幅度,单位为毫伏,厚度H指每层砂岩的厚度,单位为米;;
g、统计各井在对应地震层位的产能单井产油或气量,提取各井点对应地震层位的振幅加速度值;
h、将振幅加速度与产能或砂岩发育指数进行二元拟合,对于等速型砂泥岩地区,一般砂岩发育指数或砂岩累加厚度或单井产量与振幅加速度拟合呈幂函数关系,相关系数R在0.6以上;
i、振幅加速度平面图与用测井相平面图、油层产能平面图进行比较,高产井区或厚砂层发育区往往具有高或低的振幅加速度异常,即在平面上振幅加速度总是能把储层区和非储层区有效地识别和区分;
j、对以上属性的数值归一化,彩色显示其平面图,对于等速型砂泥岩地区,振幅加速度异常高或低值的平面展布范围即为储层发育区,含油层可能为异常低值,含气砂层可能为异常高值;
k、用三维地震体的沿层振幅加速度平面图可识别和预测油气储层是通过试错法多次反复实验发现的,并在随后多次实践中被证明是有效的,试错法实验过程为:在叠后地震数据体中,通过计算处理可以获取50多个与时间、振幅、频率、相位、波形有关的属性,大多数地质条件下,简单的振幅、时间属性即可识别和预测储层,但对于等速型砂泥岩地层,常用的这些属性均不能有效预测储层,通过在多个油气区针对多个等速型砂泥岩层计算提取了各层的各种地震属性,再根据已有钻井揭示的储层发育与分布状况,分别与每个属性平面图对比、拟合,寻找与钻井结果一致的能反映储层分布的有效属性,振幅加速度能识别和预测等速型砂泥岩储层就是这样试验出来的;在腰英台油田为了获得能够有效反映砂岩反映程度地震属性,采用试错法分别提取9个同相轴对应层位的地震振幅、频率、波形分类、振幅加速度、振幅变化率、峰度、歪度、相干体属性,并与测井相平面图、测井砂岩发育指数平面图比较,发现除低的振幅加速度与砂岩发育指数高值分布的平面位置较吻合外,其他属性变化与砂岩的发育程度之间几乎没有关系,振幅加速度对砂岩反映效果最好,大多数情况下,低的加速度对应于厚、纯、粗的高砂岩指数分布区,在各层位的振幅加速度图上勾绘了砂岩发育指数自然电位负异常幅度乘以厚度等值线图成为砂岩储层预测重要手段之一;在东北腰英台油田实例中,振幅加速度一般介于10-20之间,砂岩发育指数200-400;根据振幅加速度和砂岩发育指数的关系解释出一个振幅加速度小于5的低值异常反映出一个朵状砂体,其砂岩发育指数大于600,获得了累厚超过20m的含油砂层和高产油气流;
l、在等速型砂泥岩地区,由于泥岩普遍含碳酸钙,含钙泥岩及钙质泥岩与砂岩的自然伽玛、电阻率、声波测井响应类似,自然电位测井对砂岩的响应较明显,它成为识别砂岩的最有效曲线,用自然电位幅度乘以厚度作为砂岩的发育指数,并在基线校正和标准化的基础上,统计计算各井、各砂组的砂岩发育指数,用于与地震属性的拟合,以实现最终用地震属性识别和预测砂岩发育区;
m、腰英台油田各井点的重点产油砂组的振幅加速度与砂岩发育指数之间略呈幂函数相关性,相关系数在0.6左右,尽管相关系数较小,但与其它地震属性比较,仍算是在储层识别、预测和沉积相平面展布分析方面一个较为有效的方法;
n、等速型砂泥岩地区各层位的振幅、波形、频率的确存在着平面变化,但他们与砂岩含量关系不大,经正演模型和地震属性与钻井的对比分析这些地震属性的变化反映了低速纯泥岩含量的变化,在钻井剖面图上,这些纯泥岩的占地层总厚度不到20%,而其他80%的岩性为声波时差变化很小的砂岩、粉砂岩及含钙泥岩,所以地震强振幅平面展布范围多指示了纯泥岩的范围,但弱振幅却没有指示出砂岩的富集区,波形参数对砂岩反映的灵敏度介于振幅加速度和振幅之间,具相同波形的部分区域指示了砂岩较发育的低加速度区,部分区域却指示了砂岩不发育的高加速度区;
o、在等速型砂泥岩地区的各层位的瞬时频率平面图上,频率的平面变化很小,各层位的频率总体上均有南高北低的趋势,似乎与南砂北泥的趋势吻合,正演模型也反映含砂较厚时,频率变高的特点,对腰英台油田的DB13~DB27井产油区来说,瞬时频率几乎无变化,因而对储层预测来讲没有意义;
p、属性试验结果均表明,对目的层砂岩储层的预测效果可以用好与不好进行评价;振幅加速度对储层预测效果好,其成功之处就在于识别出了砂体沉积相类型及展布的有利位置;除振幅加速度外,包括振幅、振幅变化率和波形在内的绝大部分属性预测效果均不好;
q、振幅加速度在应用于预测的鄂尔多斯盆地北部塔巴庙气田下石盒子组气层砂岩时,振幅加速度平面图指示了沉积相、砂岩储层、含气量三种信息,它指示该区发育一个由北流向南的一个曲流河道,在河道拐弯处发育了五个点坝砂体,气藏主要位于点坝之上,钻于点坝上的qc 3、qc 4、qc 10、qc 11、qc 16、qc 17、qc 20井每日产气量超过30000m3,钻于漫滩上qc 2、qc 13、qc 21、qc 22、qc 23井的日产气量小于10000m3,产气量与振幅加速度成幂函数递增关系,相关系数达0.82;
r、鄂尔多斯盆地北部塔巴庙气田下石盒子组砂岩层速度一般为4300~4700m/s,含气后砂岩速度4000m/s左右,泥岩不纯,多由粉砂质泥岩和含钙泥岩组成,速度一般4300~4700m/s,因而与砂岩之间无明显的波阻抗差,地层中夹少量薄层的低速纯泥岩或碳质泥岩,他们的含量变化是造成波阻抗差和引起振幅变化的主要原因,因此根据波阻抗反演和与振幅有关的属性变化很难预测砂岩储层,在探井qc3井盒2-3段获得日产53000m3的天然气后,根据反演和其它地震属性在该井南部2km的范围内预测的qc 21、qc 22、qc 23井均只得到了日产数千方的气,未达到工业产能,在盒2-3段的振幅加速度平面图上,具较厚砂岩储层和较高产能的钻井均位于高振幅加速度区域内,这种黄红色的高振幅加速度平面形态指示了沿曲流河道分布的5个点坝砂体,高产井位于点坝砂体内,低产井及干井则位于决口扇或漫滩上;
塔巴庙气田砂岩储层的自然电位负异常响应反映效果不明显,而含钙粉砂质泥岩和钙质泥岩与砂岩储层均具有较低的自然伽玛测井值和曲线形态,所以,采用自然电位和自然伽玛反演识别和预测储层也不十分理想;
s、在塔巴庙气田的振幅加速度与钻井产能关系图上,振幅加速度与产能呈正相关,相关系数达0.82,也就是说振幅加速度越大,则该地区产能就越大,等速型砂泥岩地层的泥岩往往含钙质,因而其自然伽玛、电阻率、声波时差测井响应与砂岩类似,二者不宜区分,测井约束反演的波阻抗、自然伽玛的资料也无法识别和预测储层,当自然电位曲线反映砂岩效果较好时,测井约束反演的自然电位也可较好的指示储层,这时自然电位反演和振幅加速度可以互相印证,提高预测的可信度。
实施例1
参照图1所示,腰英台油田青二IV-V砂组振幅加速度与砂岩发育指数的关系图,指示工区西北角有一由南向北的朵状三角洲砂体。
对腰英台油田青二IV-V砂的计算例子
工区正北—东北许多钻井如DB15、DB16、DB17、DB22、DB25等只钻到了累加3-10m厚的油层,该区的振幅加速度介于10-20之间,砂岩发育指数200-400;根据振幅加速度和砂岩发育指数的关系,在工区西北角解释出一个振幅加速度低值异常(加速度小于5)所反映的三角洲朵状砂体,其砂岩发育指数大于600,后来据此所钻的DB33和DB34井均获得了累厚超过20m的含油砂层和高产油气流。
腰英台油田是在松辽盆地南部下白垩统青山口组发现的许多油田之一,由于泥岩普遍含碳酸钙,含钙泥岩及钙质泥岩与砂岩的自然伽玛、电阻率、声波等测井响应类似,所以只有自然电位测井对砂岩的响应较明显,它成为识别砂岩的最有效曲线。然而,自然电位的负异常幅度在平面上变化较大,总体上为工区西南部的砂岩自然电位负异常幅度较大,东北部井的砂岩负异常幅度较小,这一变化趋势与地震振幅、波形、加速度等属性变化在北东—南西方向上的分带特征是一致的,所以,它们应反映了岩性及沉积相的变化,而不是测井随机误差。但是,自然电位的负异常幅度可能只反映了砂岩的粒度信息,却并不能完整地反映砂岩的发育程度。因此用自然电位幅度乘以厚度作为砂岩的发育指数,并在基线校正和标准化的基础上,统计计算各井、各砂组的砂岩发育指数,用于与地震属性的拟合,以实现最终用地震属性识别和预测砂岩发育区。
为了获得能够有效反映砂岩反映程度地震属性,采用试错法分别提取9个同相轴对应层位的地震振幅、频率、波形分类、振幅加速度、振幅变化率、峰度、歪度、相干体等属性,并与测井相平面图、测井砂岩发育指数平面图等比较,发现除低的振幅加速度与砂岩发育指数高值分布的平面位置较吻合外,其他属性变化与砂岩的发育程度之间几乎没有关系。振幅加速度对砂岩反映效果最好,大多数情况下,低的加速度对应于厚、纯、粗的高砂岩指数分布区。在各层位的振幅加速度图上勾绘了砂岩发育指数(自然电位负异常幅度乘以厚度)等值线图成为腰英台地区砂岩储层预测重要手段之一。
腰英台地区各井点的重点产层青二IV-V和青-IV砂组的振幅加速度与砂岩发育指数之间略呈一定的幂函数相关性,相关系数在0.6左右。尽管相关系数较小,但与其它地震属性比较,仍算是在储层识别、预测和沉积相平面展布分析方面一个较为有效的方法。
各层位的振幅、波形、频率等的确存在着平面变化,但他们与砂岩含量关系不大。经正演模型和地震属性与钻井的对比分析发现这些地震属性的变化反映了低速纯泥岩含量的变化,在钻井剖面图上,这些纯泥岩的占地层总厚度不到20%,而其他80%的岩性为声波时差变化很小的砂岩、粉砂岩及含钙泥岩。所以地震强振幅平面展布范围多指示了纯泥岩的范围,但弱振幅却没有指示出砂岩的富集区。波形参数对砂岩反映的灵敏度介于振幅加速度和振幅之间,研究发现,具相同波形的部分区域指示了砂岩较发育的低加速度区,部分区域却指示了砂岩不发育的高加速度区。
在各层位的瞬时频率平面图上,频率的平面变化很小。各层位的频率总体上均有南高北低的趋势,似乎与南砂北泥的趋势吻合。正演模型也反映含砂较厚时,频率变高的特点。但就DB13~DB27井产油区来说,瞬时频率几乎无变化,因而对储层预测来讲没有意义。
属性试验结果均表明,除振幅加速度外,包括振幅变化率和波形在内的绝大部分属性对该区目的层砂岩储层的反映效果不好。
振幅加速度对储层预测较成功之处,就是识别出了砂体沉积相类型及展布的有利位置。特别是在工区西北角发现青二段下部的一个较大型三角洲朵体,而目前的油井如DB13、DB18、DB26等井仅仅钻在该三角洲的东侧缘,以浅滩相为主。朝西1~3km有更厚的三角洲前缘席状砂、水下分流河道和河口坝等朵体,规模3km×5km,可望形成一个中型规模的油田。后来在该位置所钻的DB33和DB34井验证了以上预测。
另一个重要发现就是在青一段上部振幅加速度平面图上,在目前发现油的位置圈定了沿湖岸线展布的三角洲前缘席状砂和浅滩砂分布有利异常区,并在DB16井和DB17井之间发现一个由南向北伸长的反映砂岩加厚的朵状低加速度区,由于位于深湖区,又具有蘑菇朵的平面形态,所以很可能是浊积扇体,或者为水下分流河道,宽1.5~2km,长3km。重要的是该砂体分布异常区与DB19井北部工区边界吉林油田在该层位发现高产油流的黑148井位连成了一片,似乎预示了在DB16井和DB17井之间南向北伸长的朵状范围内,也有可能找到类似的高产油气层。
实施例2
图2所示为振幅加速度预测的鄂尔多斯盆地北部塔巴庙气田下石盒子组气层砂岩展布图。振幅加速度平面图指示了沉积相、砂岩储层、含气量三种信息,它指示该区发育一个由北北流向南的一个曲流河道,在河道拐弯处发育了五个点坝砂体。气藏主要位于点坝之上,钻于点坝上的qc3、qc4、qc 10、qc 11、qc 16、qc 17、qc 20等井每日产气量超过30000m3,钻于漫滩上qc2、qc 13、qc 21、qc 22、qc 23等井的日产气量小于10000m3,产气量与振幅加速度成幂函数递增关系,相关系数达0.82见图3所示。
鄂尔多斯盆地北部塔巴庙气田下石盒子组砂岩层速度一般为4300~4700m/s,含气后砂岩速度4000m/s左右,泥岩不纯,多由粉砂质泥岩和含钙泥岩组成,速度一般4300~4700m/s,因而与砂岩之间无明显的波阻抗差。地层中夹少量薄层的低速纯泥岩或碳质泥岩,他们的含量变化是造成波阻抗差和引起振幅变化的主要原因。因此根据波阻抗反演和与振幅有关的属性变化很难预测砂岩储层。在探井qc3井盒2-3段获得日产53000m3的天然气后,根据反演和其它地震属性在该井南部2km的范围内预测的qc21、qc22、qc23等井均只得到了日产数千方的气。在盒2-3段的振幅加速度平面图上,具较厚砂岩储层和较高产能的钻井均位于高振幅加速度区域内,这种黄红色的高振幅加速度平面形态指示了沿曲流河道分布的5个点坝砂体。高产井位于点坝砂体内,低产井及干井则位于决口扇或漫滩上。
塔巴庙气田砂岩储层的自然电位负异常响应反映效果不明显,而含钙粉砂质泥岩和钙质泥岩与砂岩储层均具有较低的自然伽玛测井值和曲线形态,所以,采用自然电位和自然伽玛反演识别和预测储层也不十分理想。
图3所示为振幅加速度与钻井产能关系图,由关系图中的曲线可知,振幅加速度与产能呈正相关,相关系数达0.82。也就是说振幅加速度越大,则该地区产能就越大。
等速型砂泥岩地层的泥岩往往含钙质,因而其自然伽玛、电阻率、声波时差等测井响应与砂岩类似,二者不宜区分,测井约束反演的波阻抗、自然伽玛的资料也无法识别和预测储层。当自然电位曲线反映砂岩效果较好时,测井约束反演的自然电位也可较好的指示储层。这时自然电位反演和振幅加速度可以互相印证,提高预测的可信度。例如图1的右下方有两个椭球形的低振幅加速度异常,其中左边的一个与自然电位反演指示的砂岩发育区吻合,可能为点砂坝砂体;右边异常则与自然电位反映的砂岩不吻合,可能为浅湖或漫湖泥沉积区。
振幅加速度和自然电位反演相比,前者不需要井的资料约束,运算速度快,成本低,预测效果好。自然电位反演首先要求砂岩有较明显的自然电位响应,其次,平面上要有较均匀的井的分布,再次,当反演前的测井基线校正或标准化存在问题时,反演结果也会出现多解性或错误。

Claims (1)

1、等速砂泥岩地层的储层预测油气的方法,其特征在于:包括①计算研究区的地震振幅加速度体;②沿目的层提取显示振幅加速度平面图;③在该平面图上解释和圈定振幅加速度异常值范围作为有利油气储层平面展布范围,一个计算软件程序,其方法步骤为:
a、向计算机中加载三维地震保幅体和偏移体;
b、用三维保幅资料计算三维地震数据体的振幅加速度体,设地震反射强度为S,振幅加速度为AA,则
AA=d2(LnS)/dS(t)2,即对反射强度A取对数,再求其二阶导数;
式中d为导数,Ln-为自然对数,t-为地震反射时间,S-为地震反射强度;
c、用偏移三维地震体解释并逐线或逐道内插勘探目的层的各个反射层位;
d、根据钻井标定,按照地震的最大纵向分辨率提取平均每半个周期对应地质层位的地震振幅、振幅加速度;采用均方根振幅或绝对值值振幅法即半个周期内各采样点振幅绝对值的和;
e、对以上属性的数值归一化,彩色显示其平面图,对于等速型砂泥岩地区,振幅加速度异常高或低值的平面展布范围即为储层发育区,含油层可能为异常低值,含气砂层可能为异常高值;
f、在测井基线校正和标准化的基础上,用砂岩的自然电位负异常幅度SP乘以厚度H作为砂岩的发育指数SI,统计计算各井、各砂组的砂岩发育指数SI;即SI=SP×H,SI取值一般介于0~1000之间,其大小取决于自然电位SP基线校正和标准化后的取值,即从泥岩基线到砂岩负异常幅度,单位为毫伏,厚度H指每层砂岩的厚度,单位为米;
g、统计各井在对应地震层位的产能单井产油或气量,提取各井点对应地震层位的振幅加速度值AA;
h、将振幅加速度与产能或砂岩发育指数进行二元拟合,对于等速型砂泥岩地区,砂岩发育指数或砂岩累加厚度或单井产量与振幅加速度呈幂函数拟合关系,相关系数R在0.6以上;
i、振幅加速度平面图与用测井相平面图、油层产能平面图进行比较,高产井区或厚砂层发育区往往具有高或低的振幅加速度异常,即在平面上振幅加速度总是能把储层区和非储层区有效地识别和区分;
j、对以上属性的数值归一化,彩色显示其平面图,对于等速型砂泥岩地区,振幅加速度异常高或低值的平面展布范围即为储层发育区,含油层可能为异常低值,含气砂层可能为异常高值;
k、用三维地震体的沿层振幅加速度平面图可识别和预测油气储层是通过试错法多次反复实验发现的,并在随后多次实践中被证明是有效的,试错法实验过程为:在叠后地震数据体中,通过计算处理可以获取50多个与时间、振幅、频率、相位、波形有关的属性,大多数地质条件下,简单的振幅、时间属性即可识别和预测储层,但对于等速型砂泥岩地层,常用的这些属性均不能有效预测储层,通过在多个油气区针对多个等速型砂泥岩层计算提取了各层的各种地震属性,再根据已有钻井揭示的储层发育与分布状况,分别与每个属性平面图对比、拟合,寻找与钻井结果一致的能反映储层分布的有效属性,振幅加速度能识别和预测等速型砂泥岩储层就是这样试验出来的;在腰英台油田为了获得能够有效反映砂岩反映程度地震属性,采用试错法分别提取9个同相轴对应层位的地震振幅、频率、波形分类、振幅加速度、振幅变化率、峰度、歪度、相干体属性,并与测井相平面图、测井砂岩发育指数平面图比较,发现除低的振幅加速度与砂岩发育指数高值分布的平面位置较吻合外,其他属性变化与砂岩的发育程度之间几乎没有关系,振幅加速度对砂岩反映效果最好,大多数情况下,低的加速度对应于厚、纯、粗的高砂岩指数分布区,在各层位的振幅加速度图上勾绘了砂岩发育指数自然电位负异常幅度乘以厚度等值线图成为砂岩储层预测重要手段之一;在东北腰英台油田实例中,振幅加速度一般介于10-20之间,砂岩发育指数200-400;根据振幅加速度和砂岩发育指数的关系解释出一个振幅加速度小于5的低值异常反映出一个朵状砂体,其砂岩发育指数大于600,获得了累厚超过20m的含油砂层和高产油气流;
l、在等速型砂泥岩地区,由于泥岩普遍含碳酸钙,含钙泥岩及钙质泥岩与砂岩的自然伽玛、电阻率、声波测井响应类似,自然电位测井对砂岩的响应较明显,它成为识别砂岩的最有效曲线,用自然电位幅度乘以厚度作为砂岩的发育指数,并在基线校正和标准化的基础上,统计计算各井、各砂组的砂岩发育指数,用于与地震属性的拟合,以实现最终用地震属性识别和预测砂岩发育区;
m、腰英台油田各井点的重点产油砂组的振幅加速度与砂岩发育指数之间略呈幂函数相关性,相关系数在0.6左右,尽管相关系数较小,但与其它地震属性比较,仍算是在储层识别、预测和沉积相平面展布分析方面一个较为有效的方法;
n、等速型砂泥岩地区各层位的振幅、波形、频率的确存在着平面变化,但他们与砂岩含量关系不大,经正演模型和地震属性与钻井的对比分析这些地震属性的变化反映了低速纯泥岩含量的变化,在钻井剖面图上,这些纯泥岩的占地层总厚度不到20%,而其他80%的岩性为声波时差变化很小的砂岩、粉砂岩及含钙泥岩,所以地震强振幅平面展布范围多指示了纯泥岩的范围,但弱振幅却没有指示出砂岩的富集区,波形参数对砂岩反映的灵敏度介于振幅加速度和振幅之间,具相同波形的部分区域指示了砂岩较发育的低加速度区,部分区域却指示了砂岩不发育的高加速度区;
o、在等速型砂泥岩地区的各层位的瞬时频率平面图上,频率的平面变化很小,各层位的频率总体上均有南高北低的趋势,似乎与南砂北泥的趋势吻合,正演模型也反映含砂较厚时,频率变高的特点,对腰英台油田的DB13~DB27井产油区来说,瞬时频率几乎无变化,因而对储层预测来讲没有意义;
p、属性试验结果均表明,对目的层砂岩储层的预测效果可以用好与不好进行评价;振幅加速度对储层预测效果好,其成功之处就在于识别出了砂体沉积相类型及展布的有利位置;除振幅加速度外,包括振幅、振幅变化率和波形在内的绝大部分属性预测效果均不好;
q、振幅加速度在应用于预测的鄂尔多斯盆地北部塔巴庙气田下石盒子组气层砂岩时,振幅加速度平面图指示了沉积相、砂岩储层、含气量三种信息,它指示该区发育一个由北流向南的一个曲流河道,在河道拐弯处发育了五个点坝砂体,气藏主要位于点坝之上,钻于点坝上的qc3、qc4、qc10、qc11、qc16、qc17、qc20井每日产气量超过30000m3,钻于漫滩上qc2、qc13、qc21、qc22、qc23井的日产气量小于10000m3,产气量与振幅加速度成幂函数递增关系,相关系数达0.82;
r、鄂尔多斯盆地北部塔巴庙气田下石盒子组砂岩层速度一般为4300~4700m/s,含气后砂岩速度4000m/s左右,泥岩不纯,多由粉砂质泥岩和含钙泥岩组成,速度一般4300~4700m/s,因而与砂岩之间无明显的波阻抗差,地层中夹少量薄层的低速纯泥岩或碳质泥岩,他们的含量变化是造成波阻抗差和引起振幅变化的主要原因,因此根据波阻抗反演和与振幅有关的属性变化很难预测砂岩储层,在探井qc3井盒2-3段获得日产53000m3的天然气后,根据反演和其它地震属性在该井南部2km的范围内预测的qc21、qc22、qc23井均只得到了日产数千方的气,未达到工业产能,在盒2-3段的振幅加速度平面图上,具较厚砂岩储层和较高产能的钻井均位于高振幅加速度区域内,这种黄红色的高振幅加速度平面形态指示了沿曲流河道分布的5个点坝砂体,高产井位于点坝砂体内,低产井及干井则位于决口扇或漫滩上;
塔巴庙气田砂岩储层的自然电位负异常响应反映效果不明显,而含钙粉砂质泥岩和钙质泥岩与砂岩储层均具有较低的自然伽玛测井值和曲线形态,所以,采用自然电位和自然伽玛反演识别和预测储层也不十分理想;
s、在塔巴庙气田的振幅加速度与钻井产能关系图上,振幅加速度与产能呈正相关,相关系数达0.82,也就是说振幅加速度越大,则该地区产能就越大,等速型砂泥岩地层的泥岩往往含钙质,因而其自然伽玛、电阻率、声波时差测井响应与砂岩类似,二者不宜区分,测井约束反演的波阻抗、自然伽玛的资料也无法识别和预测储层,当自然电位曲线反映砂岩效果较好时,测井约束反演的自然电位也可较好的指示储层,这时自然电位反演和振幅加速度可以互相印证,提高预测的可信度。
CN 200510041911 2005-04-07 2005-04-07 等速砂泥岩地层的储层预测油气的方法 Expired - Fee Related CN1291241C (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 200510041911 CN1291241C (zh) 2005-04-07 2005-04-07 等速砂泥岩地层的储层预测油气的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 200510041911 CN1291241C (zh) 2005-04-07 2005-04-07 等速砂泥岩地层的储层预测油气的方法

Publications (2)

Publication Number Publication Date
CN1664616A CN1664616A (zh) 2005-09-07
CN1291241C true CN1291241C (zh) 2006-12-20

Family

ID=35035800

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 200510041911 Expired - Fee Related CN1291241C (zh) 2005-04-07 2005-04-07 等速砂泥岩地层的储层预测油气的方法

Country Status (1)

Country Link
CN (1) CN1291241C (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101551466B (zh) * 2008-04-03 2011-04-20 中国石油天然气集团公司 一种利用与偏移距有关的地震属性提高油气储层预测精度的方法
CN105590018A (zh) * 2014-11-07 2016-05-18 中国石油化工股份有限公司 砂泥岩薄互层油藏的油水层判别方法

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101183154B (zh) * 2007-11-30 2010-04-07 辽河石油勘探局 一种地质录井解释方法
CN101726761B (zh) * 2008-10-15 2012-03-07 中国石油天然气股份有限公司 一种风险约束的油气资源空间分布预测方法
CN101906966A (zh) * 2010-07-16 2010-12-08 中国石油天然气股份有限公司 储层产能的预测方法及装置
CN102445708B (zh) * 2010-10-14 2013-10-23 中国石油大学(北京) 三维等效富泥质砂岩速度预测模型
CN101980053B (zh) * 2010-10-18 2013-12-18 中国石油化工股份有限公司 一种复杂礁滩储层预测方法
CN102033242B (zh) * 2010-10-22 2012-09-26 中国石油化工股份有限公司 一种深层倾斜裂缝储层地震振幅预测方法
CN103163555B (zh) * 2011-12-12 2015-09-09 中国石油化工股份有限公司 中浅层天然气气藏识别方法
US20140222347A1 (en) * 2012-11-03 2014-08-07 Drilling Info, Inc. Seismic Waveform Classification System And Method
CN104142519B (zh) * 2013-10-29 2017-02-08 中国石油化工股份有限公司 一种泥岩裂缝油藏预测方法
CN103711482B (zh) * 2014-01-06 2016-09-21 山东科技大学 一种基于两种测井曲线重叠图快速识别碎屑岩储层的方法
CN105467450B (zh) * 2015-11-09 2017-07-21 中国石油化工股份有限公司 一种基于亮点空间振幅变化率的油藏预测方法
CN106842355A (zh) * 2016-12-27 2017-06-13 中国石油天然气集团公司 一种有利储层的表征方法及装置
CN107807386B (zh) * 2017-09-27 2019-04-09 中国石油天然气股份有限公司 一种预测点坝砂体位置的方法及装置
CN108957582B (zh) * 2018-06-22 2020-01-07 中国石油天然气股份有限公司 砂体平面形态确定方法及系统
CN109781039B (zh) * 2018-12-21 2021-01-29 中国石油天然气股份有限公司 溶洞的纵向尺度定量预测的方法
CN111044702B (zh) * 2019-10-28 2020-09-04 中国石油大学(北京) 针对陆相频繁薄互层的细粒岩储层预测方法及设备

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101551466B (zh) * 2008-04-03 2011-04-20 中国石油天然气集团公司 一种利用与偏移距有关的地震属性提高油气储层预测精度的方法
CN105590018A (zh) * 2014-11-07 2016-05-18 中国石油化工股份有限公司 砂泥岩薄互层油藏的油水层判别方法

Also Published As

Publication number Publication date
CN1664616A (zh) 2005-09-07

Similar Documents

Publication Publication Date Title
CN1291241C (zh) 等速砂泥岩地层的储层预测油气的方法
Zhou et al. Exploration breakthroughs and geological characteristics of continental shale oil: A case study of the Kongdian Formation in the Cangdong Sag, China
Odling et al. Variations in fracture system geometry and their implications for fluid flow in fractures hydrocarbon reservoirs
CN108828666B (zh) 一种正断层附近地层埋藏史的恢复方法
CN103867194B (zh) 一种砂体结构的测井表征方法与钻井层段选择方法及装置
CN108222925A (zh) 页岩气储层分级综合评价方法
CN1304856C (zh) 油气藏三元分类方法
CN104932031A (zh) 针对湖相沉积的古水深定量计算方法
BAYET‐GOLL et al. Depositional Controls on the Ichnology of Ordovician Wave‐dominated Marine Facies: New Evidence from the Shirgesht Formation, Central Iran
Wu et al. Hierarchy modeling of subsurface palaeochannel reservoir architecture
Navidtalab et al. Challenges in deep basin sequence stratigraphy: a case study from the Early–Middle Cretaceous of SW Zagros
CN110400230B (zh) 一种基于大比例尺剖面沉积微相分析的稀土找矿方法
CN104047596A (zh) 三角洲前缘沉积小层对比中小断点识别的方法
CN109324171A (zh) 一种基于岩性统计的沉积相定量识别方法
CN103278852B (zh) 利用地震数据体波形结构特征模型预测油气的方法
Sjostrom et al. Sedimentology and provenance of Mesozoic nonmarine strata in western Mongolia: a record of intracontinental deformation
Riddell Potential for freshwater bedrock aquifers in northeast British Columbia: regional distribution and lithology of surface and shallow subsurface bedrock units (NTS 093I, O, P; 094A, B, G, H, I, J, N, O, P)
CN113281820A (zh) 周口坳陷及周缘地区山西组潮坪沉积体系非常规油气有利勘探区预测方法
Wolpert Upper Jurassic Geothermal Reservoirs of South Germany: Characterization and novel exploration strategy using an integrated sequence stratigraphic approach
Flood et al. Analysis of floodplain sedimentation, avulsion style and channelized fluvial sandbody distribution in an upper coastal plain reservoir: Middle Jurassic Ness Formation, Brent Field, UK North Sea
Svinth A Sedimentological and Petrographical Investigation of the Todalen Member and the Boundary Beds of the Endalen Member.: Within the Firkanten Formation (Paleocene) in the Central Basin of Spitsbergen, Svalbard.
CN107765329A (zh) 页岩储层地震识别方法
Gubin et al. Sayan–Yenisei Potential Petroleum Area and Its Resource Potential
Hardman et al. The Warnie Volcanic Province: A Jurassic Volcanic Province in Central Australia
Woodall et al. The Stratigraphy of the Pyrenees Member of the Ravensworth, Crosby and Stickle Fields, Exmouth Sub-basin, Northwest Australia: A wave-and storm-dominated shelf-margin delta

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
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20061220