储层参数定量计算
AIT阵列感应测井在鄂尔多斯盆地致密油储层的侵入校正
AIT阵列感应测井在鄂尔多斯盆地致密油储层的侵入校正夏宏泉;吕斯端;文晓峰【摘要】钻井过程中致密油储层受泥浆侵入影响,使得阵列感应测井(AIT)视电阻率偏离地层真电阻率.以陇东地区延长组致密油储层为研究对象,利用微分几何因子法对阵列感应测井视电阻率曲线进行泥浆侵入校正以求取地层真电阻率.经过泥浆侵入校正后的阵列感应深电阻率曲线更符合储层电性特征,且用校正后的电阻率计算的含水饱和度与岩石密闭取心实测的含水饱和度较吻合,进一步印证了该方法的可靠性.%There exits an invading zone in Yanchang group in Longdong region which makes array induction resistivity deviate from the true resistivity of reservoirs .In order to get over the effect of invading ,we use geometric factor of array induction logging to correct the array induction resistivity of reservoirs and calculate the true resistivity . We used the corrected resistivity to calculate the water saturation and found it is very close to the water saturation of actual measurement ,the true resistivity of oil reservoirs is more accurate after the correction ,the new method is feasible .【期刊名称】《测井技术》【年(卷),期】2017(041)006【总页数】5页(P696-700)【关键词】测井解释;致密油储层;泥浆侵入;阵列感应测井;几何因子;地层真电阻率;鄂尔多斯盆地【作者】夏宏泉;吕斯端;文晓峰【作者单位】油气藏地质及开发工程国家重点实验室西南石油大学,四川成都610500;油气藏地质及开发工程国家重点实验室西南石油大学,四川成都 610500;中国石油集团测井有限公司长庆事业部,陕西西安 710200【正文语种】中文【中图分类】P631.840 引言鄂尔多斯盆地陇东地区延长组致密油储层泥浆侵入严重,使得测井视电阻率偏离地层真电阻率,必须对泥浆侵入的影响进行校正[1]。
应用录井资料定量评价储层流体性质
下限值。统计葡西地区葡萄ห้องสมุดไป่ตู้油层 2 8口井 5 6层的已
试 油 、压裂 资料 。做 出 了 Q与 . 系 图 ( 1 从 . 关 同 )
图1 可以看出,该区葡萄花油层压裂后有流体产出的 储层必须满 足 Q≤0 3 、 . % 也就是 说,Q> .5 ≥6
03 、 . 5 < % 为非 产层 。 6
维普资讯
第2 5卷
第 3期
大庆石油地质与开发
PG 0 D D . ....
20 06年 6月
・ 7・ 2
文 章 编 号 :10 —74 (0 6 30 2 -2 003 5 2 0 )0 - 70 0
应 用 录 井 资 料定 量 评 价储 层 流体 性 质
统计 归纳 及对 比分 析 ,建 立 了不 同区 块 的校 正 关 系 。
如葡西地区葡萄花油层的校正关系为 : S Y =10 X .5 S
+2 2 l= 1 0 XS . 5; . 6 l+1 61 T = 1 0 T+ 1 4l . ; .6 . ;
孔 隙度 / t
JT=1 l S 0 9 。 S . X + .2 其中, T S 、 T O 、X 分别为 r
张 莉 ,夏峥 寒 ,赵淑英
1 4 3 2 大陕油 田仃 井仃公 卅地质求 片分公 ・ . 尼汀 大庆 6 5;. 3 d I3 l j 64 I ( 大庆油 }有限责仟公 勘探部 .黑龙江 大厌 1 { 1
摘要:应 用录 井资料 并依 据 录 井地 质评 价储 层方 法和 准则 ,按 照储 层 流体 性 质 的 不 同建 立 了各 类评 价
层 解释 标准 如表 1 所示 。
表 1 干 层 解 释 标 准
核磁共振测井技术
B0
B1
Fiberglass Sleeve
Antenna
井眼 MREX
地层
精品课件!
精品课件!
1.利用核磁测井进行储层划分(应用实例)
泥粉晶白云岩,见岩石裂缝
泥晶白云岩 ,见纵横交错岩石裂缝和 孔洞
均为
高电
均为
阻难
有
高GR
以识
效
难以
别流
储
划分
体性
层
储层
质
移谱 明显 长拖 曳现
37号层试油, 6mm油嘴,日产 油90.46吨,水 23.34方。
象,
反映
该井岩性复杂,包括
油气
含砾砂岩、砂砾岩、
信号
玄武岩、砂岩、粉砂
岩,地层电阻率变化
电性差异 不明显
大,从常规资料上划 分储层、识别储层流
体性质有很大困难。
长T2少,短T2多。
这样的储层其吸水能力弱,在油井中开采效果差, 剩余油饱和度高,多为中低孔渗层,其孔隙的大小 孔道中都有剩余油存在,在注水过程中既要考虑驱 替作用也要考虑渗流作用,提高注水速度,加大注 水量。
核磁共振测井技术的应用
水淹级别评价
T2长且幅度大, 强水淹含水率高
长T2少,水淹相 对弱,含水率低
4
粘土束 缚体积
2
可动流 体体积
40 50 40 30 20 10
0 0.5
0
1 2 4 8 16 32 64 128 256 512 1024 2048
T2(ms)
渗透率
k c m NMR ( FFI )n BVI
核磁测井计算物性参数与实验室测量结果对比
核磁共振测井技术的应用
第4章4 储层参数测井解释模型讲解
5.4 储层参数测井解释模型
储集层物性相互之间的关系:
储集层的孔隙度与渗透率是密切相关的,但又不是简单的关系,它受颗粒 大小、分选程度、胶结程度等因素的制约。一般中粗颗粒的砂岩孔隙度大,渗 透率也大,而微细颗粒砂岩孔隙度低,渗透率也小。在孔隙度与渗透率的关系 图上,资料点的分布与粒度大小有关,粒度中值Md≤0.2mm,资料点分布在左 下方,也就是孔隙度低,渗透率也小;MD≥0.4mm的资料点分布在右上方,也 就是孔隙度大渗透率也高;0.2<Md<0.4mm的资料点基本上分布在上述两者之间。
5.4 储层参数测井解释模型
自然伽马确定泥质含量
在沉积岩石中,除钾盐层外,其放射性的强弱与岩石中含泥 质的多少有密切的关系。岩石含泥质越多,自然放射性就越强。 这是因为构成泥质的粘土颗粒较细,有较大的比表面积,在沉 积过程中能够吸附较多的溶液中放射性元素的离子。另外,泥 质颗粒沉积时间较长(特别是深海沉积),有充分的时间同放 射性元素接触和离子交换,所以,泥质岩石就具有较强的自然 放射性。这就是我们利用自然伽马测井曲线定量计算地层泥质 含量的地质依据。
三种不同的角度上提供了地层的孔隙度信息。 经验表明,如果形成三孔隙度的测井系列,无论对于高-中
-低孔隙度的地层剖面,以及不同的储层类型,一般都具有较强 的求解能力,并能较好地提供满足于地质分析要求的地层孔隙 度数据。
5.4 储层参数测井解释模型
从前面的分析可知,残余油气特别是气层对声波、 密度以及中子测井计算的孔隙度影响是不同的。
1
Shr
Nhr Nmf
高阶神经网络在储层分布参数定量预测中的应用
理 研究 成 果 的基 础 上 提 出来 的 , 目前 在 油 气 资 源开 发 中已广 泛应 用 . 人 工 神经 网络 最 主 要 的 优 点
是不需 依 赖于 参数 的结 构 模 型 , 只要 知 道 系 统 的输
入及输 出 , 就 可 以建 立该 系统 的神经 网络模 型 , 系 统 具有 很强 的非 线性 结 构 , 神经 元 之 间 以权 值 的 方 式 加 以编码 并连 接起 来 , 神 经 网络 具 有 很 强 的适 应 性 和 自学 习能 力 , 通 过一 系列 的样 本 , 可 以对 网络加 以 训练 , 识 别 出 网络 节点 间 的连接 权值 , 建 立起 可 以预 测 的神经 网络体 系.
Ab s t r ac t T hi s p ap e r p r es e nt s t he me t hod s o f hi ghe r — or de r ne u r a l ne t wor k o n t he a pp l i c a t i on o f r e s e r vo i r di s t r i but i on
p oi nt ou t s om e e x i s t e nt pr o bl em s a nd op i ni ons t O be i m pr o ve d o n p r e d i c t i ng r e s e r v o i r di s t r i b ut i on p a r ame t e r s, a nd di s — c u s s e s s om e a dv a nt a ge s a nd d i s a dv a nt a g es w he n t o a pp l y i n pr e di c t i n g r e s e r vo i r di s t r i bu t i on p a r a me t e r s . Ke y wo r d s n e ur a l n e t wor k, r e s e r v oi r, pa r a me t e r s p r e di c t i n g
储层参数
=
Δt − Δtma Δtφ − Δtma
− Vsh Δtsh Δtφ
− Δtma − Δtma
值得注意是,利用声波时差确定孔隙度时,对非压实或疏松地层需进行压实校正。 对中子测井来说,有:
CNL = CNLφφ + CNLsh Vsh + CNLma Vma 式中,CNL为中子测井值;CNLΦ、CNL sh、CNL ma分别为孔隙流体、泥质和石
Vsh
Vsh
如果 Swb<15, 令 Swb=15
最后 Swb=Swb/100
1
2)
Swb
=
1 φt
⎜⎛ ⎝
Rwb Rt
⎟⎞ ⎠
2
3) Swb = Sw 1+ B
SP
其中:B = 7.5(10 81 − 1)
SSP
SP
(10 81 − 10 81 )
4) lg(Swb) = 0.18 - (1.5lg(Md + 3.6)lg( φ ) 0.18
1 = Vsh − 1 − Vsh
Rt Rsh
Rsd
考虑到纯砂岩部分应该满足尔奇公式,即:
1 = Rsd = Rsdφsd m
S w n Fsd Rw
aRw
将该式代入上式,并整理得:
( ) Sw n
=
aR w
1 − Vsh φm
⎛ ⎜
1
−
Vsh
⎞ ⎟
⎝ R t Rsh ⎠
式中Sw为含水饱和度;φ为有效孔隙度;m、n、a为地区经验系数,一般取 n=2,m=2,a=1。
1、一种孔隙度测井方法确定孔隙度
对泥质砂岩来说,密度测井响应方程为:
DEN = ρφφ + ρsh Vsh + ρ ma Vma
无孔隙度测井条件下储层孔隙度求取方法探讨
第2 0卷 第 3期
20 0 8年 9月
岩
性
油 气
藏
V0 .0 No3 1 . 2
LT I H0L 0GI S C RE ERV0I RS
S o. 0 8 et 2 0
文 章 编 号 :6 3 8 2 ( 0 8 0 - 0 9 0 17 — 9 6 20 )3 09 — 5
无孔隙度测 井条件 下储 层孔 隙度 求取方法探讨
王晓光 , 旷红 伟 , 泽 云 , 伍 苏 静
( 江大 学油气 资源与勘探 技 术教 育部 重点 实验 室) 长 摘 要 : 对缺 乏孔 隙度 测井 系列 的井 区 , 用标 准测 井资料 , 过 “ ・ 针 利 通 岩 刻度 法” 建立起 测 井孔 隙 可
i tr r t t n mo e hr u h c r c ln t o Th r r woc mmo r st d l i h a ee tbl he P n e p ea i d lt o g o e s ai g me h d. e e a et o o n po o i mo eswh c r sa i d byS y s r du to a t ra p a e tr ssiiy Afe o e cin f co nda p r n e it t. t rc mpa i gt e et d l t h d l sa ls e y s l o e t v rn h s wo mo e swi t e mo e tb ih d b haec ntn h e pa a t r,we f u d t a r st ac l td by s a e c ntn o o iy mo li e aie y mo e c i cde t he r mee s o n h tpo o iy c lu a e h l o e tp r st de s r ltv l r o n i d wi t h p r st h t ac l td byc r n lssa d isme n a s l ede ito s0 1 2 wh l h a ea ie d va in i o o i t a l u ae o ea a y i n t a b out v ai n i . , ie t eme n r ltv e i t s y c 6 o 5 9 4 .So i i e iv d t a p y n h l o t n o o i n e p ea in mo lt ac lt e e v i r st s .3 % t s b l e h ta pli g s ae c n e tp r st i tr r t t de o c l u ae r s r o rpo o iy i e y o r ltv l r c u aei h t d r a Th r b , a i emeho ss p id t v l a er s r o ra d q nt ai ey e aie ymo ea c r t nt esu ya e . e e y af sbl t d i u ple oe au t e e v i n ua i t l e t v c lu a er s r o r r mee sb sn pelgd t. a c l t e e v i a t r yu i gt a a pa y o Ke r s: tpelg;po o i y wo d y o r st y;s a ec n e t P r d c in f co h l o t n ;S e u t a tr;a a e tr ssiiy o pp r n e it t v
常用测井曲线含义及测井解释方法
主要测井曲线及其含义一、自然电位测井:测量在地层电化学作用下产生的电位。
自然电位极性的“正”、“负”以及幅度的大小与泥浆滤液电阻率Rmf和地层水电阻率Rw的关系一致。
Rmf≈Rw时,SP几乎是平直的;Rmf>Rw时S P为负异常;Rmf<Rw时,SP在渗透层表现为正异常。
自然电位测井SP曲线的应用:①划分渗透性地层。
②判断岩性,进行地层对比。
③估计泥质含量。
④确定地层水电阻率。
⑤判断水淹层。
⑥沉积相研究。
自然电位正异常Rmf<Rw时,SP出现正异常。
淡水层Rw很大(浅部地层)咸水泥浆(相对与地层水电阻率而言)自然电位测井自然电位曲线与自然伽马、微电极曲线具有较好的对应性。
自然电位曲线在水淹层出现基线偏移二、普通视电阻率测井(R4、R2.5)普通视电阻率测井是研究各种介质中的电场分布的一种测井方法。
测量时先给介质通入电流造成人工电场,这个场的分布特点决定于周围介质的电阻率,因此,只要测出各种介质中的电场分布特点就可确定介质的电阻率。
视电阻率曲线的应用:①划分岩性剖面。
②求岩层的真电阻率。
③求岩层孔隙度。
④深度校正。
⑤地层对比。
电极系测井2.5米底部梯度电阻率进套管时有一屏蔽尖,它对应套管鞋深度;若套管下的较深,在测井图上可能无屏蔽尖,这时可用曲线回零时的半幅点向上推一个电极距的长度即可。
底部梯度电极系分层:顶:低点;底:高值。
三、微电极测井(ML)微电极测井是一种微电阻率测井方法。
其纵向分辨能力强,可直观地判断渗透层。
主要应用:①划分岩性剖面。
②确定岩层界面。
③确定含油砂岩的有效厚度。
④确定大井径井段。
⑤确定冲洗带电阻率Rxo及泥饼厚度hmc。
微电极确定油层有效厚度微电极测井微电极曲线应能反映出岩性变化,在淡水泥浆、井径规则的条件下,对于砂岩、泥质砂岩、砂质泥岩、泥岩,微电极曲线的幅度及幅度差,应逐渐减小。
四、双感应测井感应测井是利用电磁感应原理测量介质电导率的一种测井方法,感应测井得到一条介质电导率随井深变化的曲线就是感应测井曲线。
《测井地质学》第七章 测井裂缝识别与评价
DSI图像
ARI图像
FMI图像
岩芯照片 岩芯照片
二、裂缝的测井响应---- 8. 井壁成像测井
压裂诱导无效缝储层测 井响应图版
二、裂缝的测井响应---- 8. 井壁成像测井
水平缝合线图版
二、裂缝的测井响应---- 8. 井壁成像测井
层理面和泥质条带 图版
二、裂缝的测井响应---- 测井综合响应
二、裂缝的测井响应---- 5. 地层倾角测井-
2)双井径曲线反映椭圆井眼
裂缝发育往往引起井壁岩块的崩落,形成椭圆井眼,因 此可利用地层倾角仪两对相互垂直的极板所测的双井径反映 出来。一般它不会长井段出现。
3)仪器转动差异
无裂缝段一般井壁光滑,在测量过程中地层倾角仪因受 电缆钢丝的扭力均匀转动。但在裂缝发育段,井壁沿裂缝方 向的崩落,或者较大的裂缝,使仪器转动减慢、不转、甚至反 转,出现“键槽效应”。
三、裂缝有效性的测井评价及参数计算
裂缝有效性的评价
井下裂缝有效与否,决定于它的张开程度、径向延伸和 连通情况,因此裂缝有效性的评价就是对这三个因素的描述 与评价: 1.从裂缝的张开度来评价裂缝的有效性 2.从裂缝的径向延伸特征来判断裂缝的有效性 3.从裂缝的连通性和渗滤性来判断裂缝的有效性
二、裂缝的测井响应---- 5. 地层倾角测井
二、裂缝的测井响应---- 5. 地层倾角测井-
二、裂缝的测井响应---- 5. 地层倾角测井
二、裂缝的测井响应---- 5. 地层倾角测井
二、裂缝的测井响应---- 5. 地层倾角测井
二、裂缝的测井响应---- 5. 地层倾角测井
二、裂缝的测井响应---- 5. 地层倾角测井
二、裂缝的测井响应---- 8. 井壁成像测井
LPM储层定量分析工作流程
LPM调谐频率到物性分析工作流程一、钻井物性数据加载统计井点处的厚度,形成散点文件1.txt井名 X Y 厚度孔隙度sm1 346396.70 4956587.30 11.00 9.50sm3 347812.00 4957055.00 16.00 9.30sm10 345730.60 4957655.50 6.00 9.10操作流程:1、运行Data ------ Loaders and UnloadersASCII LoadControl File : Scatter Set DataTarget Surface : 选择计算调谐频率的地震解释层位 gy_001Create Control File见下图,定义Property Code: Gross_Thickness把钻井统计结果加入到地震解释层位gy_001中二、LPM操作流程2、 Geology-----LPM ---Main Selection在Main Selection选择 Surface地震解释层位和Master Grid ,选择加入的钻井数据3、地质统计分析Data AnalysisLPM Data Analysis主要完成以下两项重要任务:•计算储层物性与地震属性之间的相关性。
•如果一种或几种地震属性与储层物性之间具有较好的相关性,则定量表示二者之间的关系。
3.1选择物性数据Gross_thickness和地震属性TUNING_FREQUENCY(或其他多种)进行地质统计分析3.2 在Quality Matrix中计算Quality Matrix回归关系和储层物性与地震属性交汇图点击compute -------Calibration Editor在RegressionLine中选择Line3.3在Calibration Function校正功能中选择挑选的属性TUNING_FREQUENCY,然后Compute得到地质统计计算公式4、统计应用成图Data Population在Data Population中一步一步向下计算就行其成图过程如下:▪Apply Calibration Function-根据标定公式直接将地震属性网格转化为物性网格▪Calculate Residual Scatter-计算井点处残差▪Calculate Residual Grid-根据残差散点数据产生残差网格数据▪Apply Residual-使用残差网格校正转化的物性网格数据,以保证最终网格与井点数据一致▪Estimate Confidence-显示井点与物性网格之间可能的误差,作为评价结果可行度的指标最终的结果推荐使用Apply Residual输出的网格5、物性分析显示IESX----BasemapPost ---InterpretationGrids: 把Model Filter去掉可选所有的网格。
基于铸体薄片的致密砂岩储层孔隙微观参数定量提取技术
基于铸体薄片的致密砂岩储层孔隙微观参数定量提取技术刘颜;谢锐杰;柴小颖;李功强;葛新民;陈雨霖【摘要】Aiming at the difficult problems on the quantitative characterization of the pore structure in tight sand with low porosity and low permeability,a modified Sigma filtering image processing technology is introduced to analyze the casting thin sections,in order to gain the nicer image with high resolution and signal-to-noise ratio.The pore is extracted and binaryzation by using the methods of multi-dimensional threshold segmentation and morphology filtering.Some micro pore parameters are quantitatively calculated by using gauge point statistics techniques and topology methods such as the average pore radius,roundness,the vertical and horizontal ratio of pore.The connotative information of the pore structure in the casting thin sections is excavated effectively,and the limitations of the qualitative and semi-quantitative evaluation are broken.It is consistent that the micro pore parameters are extracted from the casting thin sections and the mercury injection in the results of the pore structure characterization.A comparison is made among the micro pore structure parameters extraction technology,the capillary pressure curves of mercury injection,and the NMR experiments.The result shows that the extraction technology has an incomparable superiority because of the simple request of the samples,the fast testing speed,and the high degree of injection of pore space with colored resin by using the high temperature and high pressure injection techniques.%针对致密砂岩储层的孔隙度、渗透率极低,物性及孔隙结构定量表征困难等问题,将图像处理技术引入铸体薄片分析中,通过改进Sigma滤波有效去除图像噪声,得到高分辨率和信噪比的铸体薄片图像.应用多维阈值分割和形态学滤波方法实现孔隙提取和二值化.在孔隙提取的基础上,应用计点统计技术和拓扑学方法实现了孔隙微观参数如平均孔隙半径、圆度、孔隙纵横比的定量计算,突破了过去对铸体薄片以定性和半定量为主的局限性,有效地挖掘了铸体薄片的孔隙结构内涵信息.将结果与压汞毛管压力参数对比可知,基于铸体薄片的孔隙微观参数与压汞孔隙微观参数在表征孔隙结构上具有一致性.与压汞毛管压力、岩心核磁共振等实验相比,铸体薄片对样品要求简单,测试速度快,且由于高温高压注入技术能将染色树脂极大程度地注入孔隙空间,在致密砂岩储层孔隙结构定量表征中具有无可比拟的优越性.【期刊名称】《河南科学》【年(卷),期】2017(035)001【总页数】5页(P134-138)【关键词】铸体薄片;孔隙结构;图像处理;计点统计;定量计算【作者】刘颜;谢锐杰;柴小颖;李功强;葛新民;陈雨霖【作者单位】长江大学地球物理与石油资源学院,武汉430100;长江大学地球物理与石油资源学院,武汉430100;中国石油青海油田勘探开发研究院,甘肃敦煌736200;长江大学地球物理与石油资源学院,武汉430100;中石化华北石油工程有限公司测井分公司,河南新乡453700;中国石油大学地球科学与技术学院,山东青岛266580;长江大学地球物理与石油资源学院,武汉430100【正文语种】中文【中图分类】P313.1致密砂岩的孔隙结构分析对于储层分类、产能预测等具有重要意义[1].目前国内外学者已成功地应用铸体薄片、扫描电镜、核磁共振及成像、CT扫描技术进行微观孔隙结构的表征[2-4].以往的科学家针对铸体薄片主要以定性描述和半定量分析为主,漏掉了其中蕴含的重要信息[5-7].针对此,本文从铸体薄片制作及分析的原理出发,应用数字图像处理技术和计点统计法、拓扑学法等进行基于铸体薄片的孔隙微观参数定量提取,深入挖掘铸体薄片中的孔隙结构信息,以实现基于铸体薄片图像处理的致密砂岩微观孔隙结构定量评价和储层分类.铸体薄片是将染色树脂或液态胶在真空下注入孔隙并磨制成薄片,在显微镜下观察孔隙、喉道及其连通性的方法,主要包括制样、显微镜观察、图像采集等步骤[8].国内对铸体薄片的描述以目估法为主,主观性强,无法获得定量、精确的矿物和孔隙信息.碎屑岩矿物主要由石英、长石等骨架、黏土及孔隙组成,不同组分在光学显微镜下呈现不同特征,在RGB像素空间中孔隙一般为蓝色,黏土一般为红色、骨架一般为黄色,颜色的变化与骨架矿物组成及含量有关.通过多元分割技术(Ehlers,1987)可从图片中提取孔隙、黏土和骨架的含量和结构.本文主要以储集空间为研究对象,仅进行孔隙参数的定量提取[9-10].图像采集过程中仪器或人为因素的影响通常会给成像信号带来一定的噪声,影响目标分割和识别效果,需要通过适当的方法进行去噪.滤波是常用的去噪手段,但滤波方法的选择至关重要,若方法选择不合理,不仅难消除噪声,反而将带来更多的噪声.滤波按类型分为空域滤波和频域滤波,较为经典的方法有低通滤波、均值滤波、中值滤波、高斯滤波、带通滤波、维纳滤波、小波变换等,这些方法在滤去噪声的同时模糊了边缘信息,使孔隙分割更为困难,因而需要寻找一种选择性滤波器,在滤去噪声的同时保留边缘信息.Sigma滤波器是一种图像增强方法,算法简单,运行速度快,更重要的是在滤波过程中能保留照片的细节信息[12].Sigma滤波器采用块操作法,通过块内均值和方差进行逐块检测实现图像去噪和信号增强.Sigma滤波的原理为:对均值为ub、方差为σ2b且符合高斯分布的信号,有95.5%的信号落在[ub-2σb,ub+2σb]间,某点的真值是其相邻几点符合[ub-2σb,ub+2σb]条件的加权平均.设N=2m+1,对于(i,j),以该点为中心,判断其周围块域内的点与该点差值与方差的关系,若绝对差值小于方差则标记为1,否则为0,即:对块域的标记值与函数进行加权平均,得到(i,j)点经Sigma滤波后的值:式中:f(i,j)out为滤波后的真值.Radu(2005)、Gu(2010)等认为信噪比低时信号的方差与噪声相近,Sigma滤波效果差,并提出改进Sigma滤波算法[13-14].首先将信号分解并用高通和低通滤波器滤波分别得到水平和垂直方向的差异和均值,然后分别对四个信号进行Sigmad滤波,最后对滤波后分量信号进行重建(图1).铸体薄片以RGB格式保存,需分别对R、G、B域的数据进行处理再合成;扫描电镜以灰度图格式保存,可直接进行处理.图2(a)是某砂岩的铸体薄片(标尺为500 mm),蓝色部分为孔隙,黄色部分为骨架.图2(b)是经Sigma滤波后的照片,图2(c)是经改进Sigma滤波后的照片.从图中可知,滤波前样品的亮白色部分和蓝黑色部分的边缘十分模糊,图片中有斑点状噪声,滤波后边缘信息更加清晰,噪声明显降低,改进Sigma滤波比Sigma滤波效果好.对铸体薄片和扫描电镜进行滤波后,孔隙、黏土、骨架等信息明显增强.此时,通过目标分割可将孔隙等信息进行定量提取.阈值分割是一种简单有效的目标提取算法,用一个或多个阈值将图像按灰(色)度分为几个部分,将隶属于同一范围的像素视为相同目标.阈值分割法有双峰法、迭代法、最大类间方差法等.本章应用多维阈值分割法进行孔隙信息提取.对于铸体薄片,孔隙、骨架和黏土在R、G、B空间中具有不同的色度特征,各色度域的灰度值可表示为:式中:C表示为R、G、B.定义阈值分割的操作算子T为:式中:tlow和thigh分别表示某组分的色度下、上阈值,阈值采用试算法得到. 将铸体薄片进行阈值分割后得再通过形态学滤波压制阈值分割带来的斑点状噪声,即可得到表征孔隙结构的二值图.计点统计(Point Counting Technology)是一种图像目标测量的方法.通过对图像中目标边缘识别与检测,得到目标的分布状况.在点统计基础上,应用拓扑学方法对每个孔隙进行分析可得孔隙面积(pore aea)、周长(pore perimeter)、费雷特直径(Feret’s diameter)等参数.孔隙面积定义为孔隙像素点的个数,孔隙周长定义为由孔隙像素所占空间所连接的长度.费雷特直径定义为以孔隙边缘某像素点为原点,在顺时针(或逆时针)方向以某一角度为基本单元进行旋转,并测量旋转轴与孔隙边缘另一像素点之间的空间距离(图3)[15-16].以角度θ旋转可得最大费雷特直径Dmax和最小费雷特直径Dmin:定义形状因子(Shape Factor)为:式中:A为孔隙面积;P为孔隙周长.如图4所示,规则圆孔的形状因子为1,规则椭圆及不规则圆孔的形状因子小于1,不规则椭圆孔的形状因子远远小于1,孔隙的形状因子随着孔隙形状的复杂而变小.孔隙纵横比(aspect ratio)定义为最大费雷特直径与最小费雷特直径比:根据定义可知,规则圆孔的纵横比(AR)等于1.孔隙纵横比随着孔隙的不规则程度的增大而增大.孔隙的等效直径(equivalent diameter)可表示为:圆度(roundness)是描述孔隙与圆的近似程度,可表示为:孔隙比表面(surface area)可写为:美国国立卫生研究院提供了基于JAVA的图像处理软件ImageJ,计点统计及孔隙参数算法以该软件为基础[17-19].图5是铸体薄片经二值化滤波和反二值化后的孔隙统计.图中黑色部分为孔隙,白色部分为骨架,共有81个孔隙单元.图6(a)是压汞法排驱压力与铸体薄片法孔隙等效直径的关系.排驱压力随着孔隙等效直径的增大而减小,两者呈良好的幂指数关系,复相关系数达0.812:式中:Pc为排驱压力,MPa;De为孔隙等效直径,mm.图6(b)是压汞法变异系数与铸体薄片法圆度的关系.孔隙变异系数与圆度成反比,规则孔隙的变异系数低,复相关系数为0.814:以上分析可知,基于计点统计和拓扑学方法所得到的微观孔隙结构参数与压汞法微观孔隙结构参数相关性强,可通过铸体薄片处理得到微观孔隙结构参数用于孔隙结构定量表征.1)采用改进Sigmad滤波和形态学滤波方法对铸体薄片进行处理,在滤去噪声的同时提高了图像的分辨率,保留了图像的边缘信息.2)从铸体薄片和背散射扫描电镜中提取的微观孔隙参数与压汞法孔隙结构参数间存在良好的相关性,可用铸体薄片实现孔隙结构进行定量评价.【相关文献】[1]郝乐伟,王琪,唐俊.储层岩石微观孔隙结构研究方法与理论综述[J].岩性油气藏,2013(5):123-128.[2]何涛,王芳,汪伶俐.致密砂岩储层微观孔隙结构特征——以鄂尔多斯盆地延长组长7储层为例[J].岩性油气藏,2013(4):23-26.[3]盛军,孙卫,赵婷,等.致密砂岩气藏微观孔隙结构参数定量评价——以苏里格气田东南区为例[J].西北大学学报:自然科学版,2015,06:913-924.[4]刘晓鹏,刘燕,陈娟萍,等.鄂尔多斯盆地盒8段致密砂岩气藏微观孔隙结构及渗流特征[J].天然气地球科学,2016,(7):1225-1234.[5]赵明,郭志强,卿华,等.岩石铸体薄片鉴定与显微图像分析技术的应用[J].西部探矿工程,2009(3):66-68.[6]方少仙,何江,侯方浩,等.鄂尔多斯盆地中部气田区中奥陶统马家沟组马五_5-马五_1亚段储层孔隙类型和演化[J].岩石学报,2009(10):2425-2441.[7]张婷,徐守余,王子敏.储层微观孔喉网络图形识别方法[J].吉林大学学报:地球科学版,2011(5):1646-1650.[8]国家发展和改革委员会.SY/T 6103—2004岩石孔隙结构特征的测定-图像分析法[S].北京:石油工业出版社,2004.[9] Ehlers E G.Optical mineralogy,theory and technique[M].UK:Oxford Press,1987.[10]范宜仁,徐拥军,范卓颖,等.铸体薄片孔隙参数提取及其对饱和度指数的影响[J].测井技术,2014(3):257-261.[11]葛新民.非均质碎屑岩储层孔隙结构表征及测井精细评价研究[D].青岛:中国石油大学(华东),2013.[12] Fens T W.Petrophysical properties from small rock samples using image analysis techniques[M].Delft,the Netherlands:Delft University Press,2000.[13] Bilcu R C,Vehvilainen M.Modified sigma filter for noise reduction in Images[C]//Wseas International Conference on Commu⁃nications,2005.[14] Gu M R,Kang D S.Modified sigma filter using image decomposition[C]//Musp’10 Proceedings of the 10th Wseas International Conference on Multimedia Systems&Signal Processing,2010:193-198.[15]涂新斌,王思敬.图像分析的颗粒形状参数描述[J].岩土工程学报,2004(5):659-662.[16]张家发,叶加兵,陈劲松,等.碎石颗粒形状测量与评定的初步研究[J].岩土力学,2016(2):343-349.[17] Igathinathane C,Pordesimo L O,Columbus E P,et al.Shape identification and particles size distribution from basic shape parameters using image[J].Computers and Electronics in Agriculture,2008,63(2):168-182.[18] Tajima R,Katob arison of threshold algorithms for automatic image processing of rice roots using freeware image[J].Fuel and Energy Abstracts,2011,121(3):460-463.[19]宋玉丹,秦志钰,容幸福.用Imagej提取图像边缘的方法及展望[J].机械管理开发,2008(4):180-181.。
煤储层参数测井定量计算方法探讨
效厚 度等 ,是 研究 煤层 组 分和 评价 煤层 气 的地质 勘探 、T 业分 析及 有效 开发 的依 据 。上述 参数 一般 由钻
井取 芯后 对煤 层 岩心进 行 实验 测定 得 出 ,但 随着 煤层 气 区块 的规模 开发 ,大 量单 井煤 储层 参 数 的评 价 只
能 由测井 资料 进 行计算 求 取l ] 由于煤储 层 的非 均质 性 及 煤 岩成 份 的复 杂性 与可 变 性 ,煤 层低 孑 、低 _ 。 1 L
归 关 系的方 法来 解决 这 一 问题 ( 图 2 。 见 ) 表 1 示为研 究 区块某 井 的计算 实例 ,计 算结果 与 所 试验 分析数据 的误 差在 3 %以内 ,能够 满足使用 要求 。
d i 0 3 6 /.sn 1 7 —4 9 ( o :1 . 9 9 jis. 6 3 1 0 N). 0 2 0 . 1 21 .6 06
煤 储 层 参 数 测 井定 量 计 算 方 法探 讨
杨克兵 , 钱 铮 , 守 刚 , 双 志 ( 到\ 焦 中石油华北油田分公司勘探开发研究院, 北 任丘 025) 河 652
进行 统计 ,结合 测井 响应 特征 ,在 考虑 煤岩 骨架 参数 变 化 的基 础 上 ,建立 了一套煤 储层 组分 及孔 隙度 定
量计算 方法 。此外 ,对含 气量 的计算 考虑 了电阻率的影 响 ,提 出了依 据电阻率对 含气量进行 校正 的方法 。
1 煤 组 分 计 算 方 法
煤 岩组 分划 分为 固定 碳 、灰分 、挥 发分 、水 分 ,根据 试验 分析 结果 ,煤 岩组 分 固定碳 、灰 分相 互 间
第 9卷 第 6期
杨 克 兵 等 :煤储 层参 数 测 井 定 量 计 算 方 法 探 讨
储层有效性综合评价方法
储层综合评价方法储层评价是预测和评价研究区含油气有利区带的重要技术手段,是对储层研究的综合认识和评判。
针对单因素评价储层结果不惟一的缺点,本文研究了储层综合定量评价的方法(图1)。
该方法分为4个步骤:首先利用特征选择算法对评价参数进行筛选,然后根据灰关联分析来确定各影响因素的权重,进而运用最大值标准化法确定各项参数的评价分数,最后计算各项参数综合得分,在此基础上,运用聚类分析进行储层分类评价。
对储层评价结果进行统计分析,所划分的各类储层特征明显,与研究区储层实际特征具有很好的一致性。
最大限度地应用计算机手段对油气储层进行精细评价和综合解释具有定量化、地质意义明确等优点,有一定应用价值。
图1储层综合评价方法体系框架1评价参数的选择一项参数只从一个方面表征储层的特性,全面评价一个储层,需要采用多项参数,从多个方面进行综合评价。
对储层进行合理的分类是评价储层的基础。
迄今为止,国内外学者提出了许多储层分类的参数与方法,但应该用哪些参数、选用何种方法是一个很难解决的问题。
而储层综合定量评价是在选取储层评价参数的基础上,对储层多个影响因素做综合评价,最终得到一个综合评价指标,并依此对储层分类。
国内研究储层的学者在评价参数选择方面作了不少研究,所选择的评价指标也各不相同。
比如:刘吉余等[7]认为储层综合评价的参数主要为储层的有效厚度、砂体钻遇率、渗透率、有效孔隙度、泥质含量、黏土矿物类型、孔隙结构参数、层内非均质性参数及隔(夹)层的分布参数等;吕红华等[8]选择孔隙度、小层厚度、含油饱和度及小层钻遇率4个参数作为储层评价指标;张晓东等[9]选择有效厚度、沉积相、夹层频数、孔隙度和裂缝渗透率5个参数作为储层评价指标;张琴等[10]选取孔隙度、渗透率、颗粒分选、杂基含量、粒径、储层成岩相带、溶蚀作用及胶结作用共8个参数进行储层评价;马立文等[11]选取孔隙度、渗透率、泥质含量及渗透率突进系数4个参数进行储层评价。
一种浊积岩体储集层参数数学地质研究方法
一种浊积岩体储集层参数数学地质研究方法随着现代地质科技的发展,地质学家们日益关注浊积岩体储集层参数数学地质研究方法。
浊积岩体储集层参数数学地质研究方法是对储集层参数的研究,包括储层厚度、岩性、孔隙度、物性等。
这是一种比较新的研究方法,它可以帮助更加精确、准确,更有效的确定和预测浊积岩体储集层的性质和参数。
浊积岩体储集层参数数学地质研究方法包括一系列的定量模型,分析和处理技术,以及实验室测试方法和计算机仿真。
在定量模型方面,以简单统计逻辑、混合模型、灰色理论、回归分析等为主要定量模型,其主要任务是建立准确的数据模型,以预测储层参数的变化规律。
在分析与处理技术方面,以地质描述软件、地震处理软件、GIS 系统等为主要分析与处理技术,其主要任务是通过进行地震资料分析和GIS处理,帮助地质学家快速确定和预测储层参数,并进行准确识别。
在实验室测试方法方面,以英语仪器、地震仪器、油藏测试仪等为主要实验室测试方法,它们可以对储层参数进行更深入的测试,从而进一步提高学者对储层参数的准确预测。
在计算机仿真方面,以计算机图像分析、半定量地质物理模型、地震波模拟等为主要计算机仿真技术,这些技术能够更好地复制实地储层物理现象,从而提高储层参数的准确性。
浊积岩体储集层参数数学地质研究方法可以有效地提高石油勘探开发的成功率,作为一种新的研究方法,它具有很多优势,如可以更加准确准时的预测储层参数变化,更好的挖掘储层的矿产有效性,还可以帮助企业更好地实施矿产资源的开发和利用。
总之,浊积岩体储集层参数数学地质研究方法是一项新的、有效的研究方法,可以帮助地质学家们准确地预测、识别储层参数,提高石油勘探开发的成功率,并为地质科技的发展做出重要贡献。
综上所述,浊积岩体储集层参数数学地质研究方法无疑是一种重要的研究方法,它能够帮助地质学家们更好地预测和识别储集层参数,以提高石油勘探开发的成功率。
这一研究方法不仅有助于改善石油勘探开发的效率,也为地质科技的发展做出重要的贡献。
三维定量储层地质模型与参数动态预测方法
三维定量储层地质模型与参数动态预测方法
胡侃蔚;陈佩珍;王冠贵
【期刊名称】《石油天然气学报》
【年(卷),期】1998(000)004
【摘要】提出了三维定量储层地质模型及其参数动态预测方法,用该方法可预测储层开采状态,水淹区及剩余油区分布和各剩余油区的可采储量,三维定量储层地质模型的预测图件可为合理开采提供依据。
【总页数】1页(P45)
【作者】胡侃蔚;陈佩珍;王冠贵
【作者单位】上海-钢集团有限公司职工中专;上海-钢集团有限公司职工中专【正文语种】中文
【中图分类】P618.130.2
【相关文献】
1.煤层气储层参数多信息综合定量预测方法 [J], 胡朝元;彭苏萍;赵士华;程增庆
2.储层参数定量预测方法探讨 [J], 刘吉余;云金表;吕靖
3.储层三维定量地质模型的建立 [J], 汤春云;许发年
4.复杂断块油藏三维地质模型多参数定量评价 [J], 谭学群;刘云燕;周晓舟;刘建党;郑荣臣;贾超
5.数字岩心逆建模理论下的储层参数定量预测方法 [J], 印兴耀;郑颖;宗兆云;林利明
因版权原因,仅展示原文概要,查看原文内容请购买。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
进一步认识后的形象化描述。 即:
理论分析 实验研究
测井信息 地质信息
资料统计
测井信息 = F(地质信息) + 误差
(测井解释模型)
用数学式表达:y=F(x)+ε
仪器模型: 测井数据集、测井 值与记录值转换关 系、测量误差、环 境校正等
地质模型: 矿物成分、 流体类型等
数学模型: 数学关系式
评价地层
• 一、阿尔奇(Archie)公式
R01 R02 ... R0n F F R0
Rw1 Rw2
Rwn
Rw
• —称为地层电阻率因素
• —全含水地层电阻率
• —地层水电阻率
• Archie还通过实验发现:
• F与孔隙度及岩石孔隙结构有关(反映岩石本身
物理性质),与孔隙中是否含油气及Rw无关。
而Swt Swb Swf
因此Ct t2Swt SwbCwb Swt Swb Cwf
即Rt
t2
Swb Rwf
Rwf Rwb Swt Swb
Rwb
S wt
5.4 渗透率计算
• 一、影响因素
• 孔隙度:孔隙度↑,渗透空间大,渗透率↑ • 泥质含量:Vsh ↑,阻塞渗透通道,K↓ • 砂岩颗粒大小(用粒度中值度量):越细,k越小。 • 裂缝:裂缝愈发育,K ↑ ↑ • 压力对渗透率的影响
§5 储层参数定量计算
5.1 测井解释模型
一. 模型的基本含义
所谓模型:就是客观事物被认识后,经过抽象 ,再组合为易于理解的形象,即形象化的抽象。 模型的建立过程是:
实际 经 过 抽象
变成
简 化 形象
二. 测井解释模型
指测井信息与地质信息间的宏观关系。是经过理
论分析、实验研究、资料统计,对两者之间的关系
Qv
CEC(1t )G t
t 泥质砂岩的总孔隙度,小数;
G 岩石的平均颗粒密度, g / cm3
CEC 岩石的阳离子交换能力, mmol / g
Qv 岩石的阳离子交换容量, mmol / cm3
• 2. Waxman- Smits饱和度方程
Cwe Cw BQv / Swt
Cw
B 3.83(1 0.83e 2 )
X
hh
ww
ma ma
sh sh
Vma Vsh 1
Sh Sw 1
XX X V XX XX
ma
ma
X
w
sh
ma
ma
sh w
S XX XX 1 h
ma
ma
h w
1
Sh 0
不含油气时
X
X X X
ma
ma
sh
V
sh
X X X X ma
w
ma
w
T Vshsh
(泥质校正公式)
S
1 Cp
t t ma t f t ma
此处, Cp 为压实校正系数,可由平均时间公式与
真孔隙度值对比得到,即 Cp s 。 由岩
心实验测定或其它孔隙度测井得到。另外,C 深度有关,即
p
与埋藏
Cp a bH
实际资料表明,特别是在地层孔隙度较大时,声波 时差与孔隙度关系不再是线性关系,而是非线性关系 。
• 这种导电性与可发生阳离子交换的数量相关,该量称阳离子交换能力 CEC,它表示每100克干岩石可交换的阳离子数。单位为mmol/cm3 。
• 阳离子交换能力与地层的粘土矿物成分有 关。如:绿泥石为10-40,伊利石为10-40 ,高岭石为3-15,蒙脱石为80-150。
• 还有一种表示阳离子交换能力的方法,称 阳离子交换容量 Qv ,它表示岩样每单位 总孔隙体积中含有的可交换阳离子的摩尔 数,单位为mmol/cm3。
• 设自由水与束缚水混合电导率为Cwm
由阿尔奇公式(a b 1, m n 2)有
Ct S wt t 2 Cwm
而Cwm
bCwb b
f Cwf f
所以Ct
S wt t
2 bCwb f Cwf b f
上式同除以t ,则Ct
t
2
S
2 wt
S wb Cwb S wb
Swf Cwf Swf
含油气 泥质砂 岩体积 模型
砂岩骨架 Vma
泥 质 Vsh
油 气 Sh
水 Sw
V=1
设我们研究的物理参数为X,用
X ,X X X , ,
h
w
ma
sh
分别表示该
物理参数对油气、水、骨架、泥质的测井响应值,由体积物理
模型的思想,测井响应值X应为各部分贡献之和。即有:
X S X S X V X V
• B为交换阳离子的当量电导率
• 由纯砂岩阿尔奇公式
•有
S n* wt
F*
Ct Cwe
S
n w
F
Ct Cw
• 所以
Ct
S n* wt
F*
(Cw
BQv
/
Swt )
• n*相当于岩石不含粘土的饱和度指 数,常取2
• F*相当于岩石不含粘土的地层因素
F*
a
tm
• 三、双水模型
• 1. 双水模型的实质
• 泥质砂岩中有两种水:束缚水,即束缚在粘土表 面的水,又称近水;自由水,即位于连通孔隙空 间的水,或称远水。也就是说,泥质的存在使地 层水的导电性比纯地层下的导电性好,并认为这 种增加的导电性是粘土颗粒表面的束缚水所产生 的阳离子交换所引起的,泥质地层的全部导电能 力是由自由水和束缚水并联导电的结果。
2.全含水纯地层孔隙度方程
所谓纯地层,即不含泥质的地层,即 Vsh 0
,若不含油气,则以上通式可写成:
X X
ma
X X
ma
w
X X
ma
X X
ma
f
由于孔隙度系列测井探测 深度浅,在探测范围内孔 隙中流体为泥浆滤液
1)密度测井:将X写成ρ,则
D
ma ma
b f
• 2)中子孔隙度:
N
H Hma Hf Hma
100
电阻率增大系数
10
y = 1.0175x -1.8896 R2 = 0.9709
1
0.1
1
含水饱和度,小数
• 原始地层含水饱和度
Sw
n
abRw
mRt
• 在冲洗带中,Rt → Rxo,Rw → Rmf,Sw → Sxo,得:
Sxo
n
abRmf
m Rxo
• 阿尔奇公式适用于纯地层,对泥质地层,计算饱和度的 公式很多,如:
2
S
2 w
• 二、Waxman-Smits模型
• 1. Waxman- Smits模型
• Waxman- Smits认为:在粘土表面均具有负电荷,这样在粘土表 面吸附着部分阳离子,这些阳离子又能很容易地与水溶液中的其它 阳离子发生交换,粘土表面及其附近的的高浓度阳离子是造成所谓 粘土表面导电性的主要原因。
• Simandoux公式
•
1 V 常取m=n=2c;dl d=1n2—2,常取d=1 m
S w
d
R R aR (1V t
cl
w
cl
)
Swn
• “尼日利亚”公式
• α=1-2
1
Vcl
Rt Rcl
• “印度尼西亚公式”
2
aRw
S
2 w
• C=1-Vcl/2
1
Vccl
Rt Rcl
aRw
即:
F
a
m
•
取决于孔道的弯曲程度,颗粒的形状和排列方向,以及胶结
情况。常称为“胶结指数”,“结构指数”,“孔隙度指数”,
“Archie指数”。取值范围为1.3-3,随胶结程度的变好,该值增
大,常取2。
•
为与岩性有关的比例常数,变化范围为0.6~1.5,常取1。简
称为“Archie常数a”
• Archie还通过实验发现:地层电阻率与R0有关,同时随So 的增大而增大。定义:
菱面体排列,φ=25.96%,最小,排列最紧密, 最“稳定”
6.胶结作用:
总的来说,胶结作用愈强,φ越小。另一方面, 钙质胶结比泥质胶结更致密,即φ更小;硅质胶结 最致密,对孔隙度影响最大。即泥质好于钙质好于 硅质。
7.白云化作用:
白云化作用越强, φ越大。主要是因为产生白 云化作用后,骨架体积变小;灰岩被白云岩置换后 体积缩小12—13%
ma 1 2 f
Raymer公式
3.如何选择计算孔隙度的方法
按优先顺序:
a.根据以往经验
b.井眼条件好时,选密度、声波、中子孔隙度测 井中的两者或三者交会,即用最优化方法求
c.中子质量不好或有气存在,选声波、密度或其 交会图
d.用电阻率方法计算孔隙度时,只能计算水层
5.3 含水饱和度计算
Sw2b
• Swb,%;φ,%;K,10-3μm2
• 3.由粒度中值与孔隙度求K
• K=f(Md,φ) • Md一般可由GR或SP求
• 4.用裂缝宽度估计K
• K=8.3×106w2 • K:10-3μm2 • W为裂缝宽度,cm
• 5.由核磁共振测井求K
• SDR模型
K
C1
(nmr
100
)
4
T22g
• 由岩电实验作F~Φ关系图求a、m;作I~Sw关系图求b、n 。
50Mpa-90℃条件下地层因素与孔隙度关系 10000
1000
地层因素F
100
10
y = 0.8051x-1.9029
R2 = 0.9793