正异常（Rw>Rmf）
泥岩基线
自然伽马
低
高
井径
缩径
扩径
深中浅电阻率
高阻
低阻
声波
<300us/m
>300us/m
钍
低
高
铀
低
高
钾
低
高
（2）碳酸盐岩剖面的测井曲线特征
碳酸盐岩剖面的测井解释任务，就是从致密的围岩中找出孔隙性、裂缝性的储集层，并判断其含油性。
碳酸盐岩剖面电阻率一般较高，自然电位效果不好。
为了区分岩性和划分储层，一般使用自然伽马测井曲线。
储集层相对于致密的围岩具有低阻、低自然伽马以及孔隙度测井反映孔隙度较大的特点。
3．划分储集层的基本方法与原则
基本要求：
凡一切可能含油气的地层都要划出来，要适当地划分明显的水层。
具体要求为：
（1）估计为油层、气层、油水同层和含油水层的储集层都必须分层解释。
（2）厚度半米以上的电性（测井曲线）可疑层（即指从测井曲线上看有油气的地层）或录井显示为微含油级别以上的储集层必须分出。
（3）选择出作为确定地层水电阻率Rw的标准水层（厚度大、岩性纯、不含油）要划分出来。
（4）录井、气测有大段油气显示而测井曲线显示不好的储集层，应选取一定层位，尤其是该组储层的顶部层位，进行分层。
（5）当有多套油水系统，油层组包括若干水层时，只解释最靠近油层的水层。
（6）对于新区探井，应做细致工作，对各个储层均应酌情选层解释，以使不漏掉可能有油气的地层。
4．正确划分出储集层的方法
（1）砂泥岩剖面
通常是自然电位（SP）曲线的异常确定渗透层的位置，用微电极曲线确定分层界面，分层前，应将井场收集的井壁取芯、气测显示等有关油气显示的资料标注在综合测井曲线图上，并根据邻井的测井和试油等资料对本井的油水关系作出初步估计。
分层时应注意：
确定分层的界面深度时，应左右环顾，照顾到分层线对每条测井曲线的合理性。
分层的深度误差不应大于0.1m。
渗透层中，凡是0.5m以上的非渗透性夹层（泥岩或致密层），应将夹层上下的渗透层分两层解释。
岩性渐变层顶界（顶部渐变层）或底界（底部渐变层）分层深度应在岩性渐变结束处。
一个厚度较大的渗透层，如有两个以上解释结论，应按解释结论分层。
在同一解释井段，如果油气层与水层岩性、地层结构和孔隙度基本相同，则油气层是纯水层的电阻率的3-5倍。
纯水层的自然电位异常最大，油气层异常明显偏小，油水同层介于油、水层之间。
并且厚度较大的油水同层，自上而下电阻率有明显减小的趋势。
（2）碳酸盐岩剖面
碳酸岩盐剖面划分渗透性地层的
5．测井曲线读数
分层以后，要从有关的主要测井曲线将代表该储层的测井曲线读数，以便计算孔隙度、饱和度等地质参数，在厚度较大的储集层中按测井曲线变化确定几个取值区，对每个取值区对应读数计算，几种主要测井曲线取值区的最小厚度如下：
各种孔隙度测井≥0.6m。
侧向测井≥0.6m
感应测井，低阻≥0.6m，高阻层≥1.5m。
每种测井曲线分层和取值要符合其方法特点，例如声波测井扣除致密夹层，选用与渗透层相对应部分的平均值。
电阻率测井曲线则扣除致密夹层，选用与渗透层相对应部分的极大值的平均值。
另外注意孔隙度与电阻率测井曲线对应取值的原则。
因为要用两者结合计算地层的含水饱和度，两者当然应该是对应深度上同一地层或同一取值区的读数。
岩层含油性的定性判断，主要依据井曲线的测井曲线特征，而电性特征是岩石物性、岩性和含油性的综合反映。
因此在判断地层的含油性时，一般应将测量井段首先按照地层水矿化度的不同分为不同的解释井段，然后才有可能对每一个解释井段在充分考虑其岩性特点的前提下进行含油性解释。
由于地下地层复杂性，仪器的局限性，上述原则是一般性的。
要做到正确地解释，一方面应多收集资料，认真分析曲线，另一方面还要了解区域性特点和规律，要积累经验。
6．计算出孔隙度、饱和度等参数。
读数以后，还要做一些定量计算，常用的公式：
孔隙度：
含水饱和度：
上式中
为当前层的声波时差，
为地层水的声波时差,189us/ft（623us/m），
为固体骨架的声波时差，对于砂岩骨架，主要矿物为石英，其声波时差为ft。
a是常数，对于砂岩地层通常取，
为当前层的电阻率，m为胶结指数。
五、实验结果
处理结果如图2所示
具体分层如下表所示：
储层序号
顶部深度
（m）
底部深度（m）
厚度（m）
测井曲线读数
孔隙度
含水
饱和度
SP
GR
DEN
1
1142
2
1162
3
1226
4
5
由自然电位的划分原则可知，在图2的自然电位的曲线上，有五个明显的偏离泥基线的负异常，故可以划分出五个储集层。
然后就可以求出泥质含量，进而求出泥质校正后的孔隙度，然后根据阿尔奇公式，就可求出含水饱和度、渗透率，进而划分油气水层。
图1某井的综合测井曲线图
1
2
3
4
5
图2SN183井测井解释综合图
实验二 利用综合方法估计地层泥质含量
一、实验目的：
通过实际计算，巩固掌握利用多种测井资料确定泥质含量的方法。
二、实验要求
自编程序，在计算机上运算出地层泥质含量。
三、实验场地、用具与设备
计算中心，尺子、像皮和计算机；
四、实验内容：
1．什么是泥质含量：
泥质是指颗粒直径小于0.01mm的碎屑物质，泥质含量，也叫做泥质体积，是指泥质的体积占岩石总体积的比：
２．确定Vsh 的重要性
泥质含量的确定，在泥质砂岩储集层的定量解释中具有重要意义。
多年来人们提出许多计算泥质含量的理论和方法。
目前求取泥质含量的方法大致可分为两类，一类是用每种测点各求出一个泥质含量，然后求出最佳值。
当岩石含有泥质时，各种测井曲线均或多或少地受到泥质的影响，其影响的程度受Vsh 的决定，评价岩石的特性时，只有已知Vsh ，才知道由于泥质带来的影响，从而将泥质的影响校正掉。
一般而言，用自然伽马或自然伽马能谱或自然电位来求取泥质含量效果最好，但自然伽马要求储层中除了泥质外，其他物质不含放射性矿物。
自然电位要求地层水电阻率保持不变，且储层中的泥质与相邻泥岩的的成分相同。
用其他方法计算泥质含量则要求更为苛刻的条件：
如电阻率方法要求储层的孔隙度和含水饱和度均要很小。
中子和声波方法则要求孔隙度很小。
3.确定Vsh 的方法：
（1）自然伽玛法
式中，
分别是砂岩和泥岩层的自然伽马值，GCUR是与地层有关的经验系数，新地层（第三系地层）GCUR=，老地层GCUR=.
（２） 自然电位法
式中，
是当前层的自然电位读数，
和
分别是纯地层和泥质地层的自然电位读数
（３）电阻率
（b=）
（４）中子法
式中，
是当前层的视中子孔隙度读数，
是泥岩层的视中子孔隙度读数。
（5）交会图法
以中子—密度测井交会图为例，通过对图2所示的石英点（Q）、水点（W）和泥岩点（SH）构成的三角形进行分解，依据资料点所落入三角形中的位置，可以推测出来泥质含量。
或者利用下式进行计算（依据点到直线的距离计算方法）：
式中，
SH
W
Q
＝0是石英点（Q）和水点（W）连线的直线方程。
依据任意两点的直线，用石英点（
）和水点（
）两个点的参数可以推出：
，
，
∴
当然，也可以用中子—声波、声波—密度交会图的类似方法求Vsh。
如果解释层段上没有纯泥岩层时，上述交会图法所定出的泥岩点位置并不代表实际的泥岩参数，导致用交会图方法估计的泥质含量比实际值偏大。
上述几种方法计算的泥质含量往往都有一定的条件，当条件满足时，泥质含量的近似结果，当条件不满足时，计算的泥质含量均可能偏高，所以在实际处理时，最后选取其中的最小值作为接近实际的泥质含量。
五、实验结果
1.已知条件：
GRmin=50GRmax=110SPmin=0SPmax=100b=Rsh=5
2.测井数据
depth
sp
msfl
lls
lld
rhob
nphi
gr
dt
cal
3.计算结果
sp
msfl
lls
lld
rhob-nphi
nphi
gr
泥质含量
附录：
程序源代码
!
功能:
利用测井数据求泥质含量
!
时间:
2016/12/22
!
----------------------------------------
!
参数说明:
!
cmdfile:
参数文件
!
file_data:
输入测井数据文件
!
file_contant:
常数参量文件
!
file_sp:
输出自然电位估计的泥质含量文件
!
file_msfl:
输出微球型聚焦估计的泥质含量文件
!
file_lls:
输出浅侧向估计的泥质含量文件
!
file_lld:
输出深侧向估计的泥质含量文件
!
file_nphi:
输出补偿中子估计的泥质含量文件
!
file_gr:
输出自然伽马估计的泥质含量文件
!
file_n_d:
输出中子-密度交会图估计的泥质含量文件
!
n:
测井数据的个数
!
grmin,grmax:
自然伽马的最小值与最大值
!
spmin,spmax:
自然电位的最小值与最大值
!
b,r_sh:
已知常数
!
fai_nsh:
已知常数
!
p_sh,p_ma,p_f:
已知常数
!
fai_nf,fai_nma:
已知常数
!
sp（:
）:
自然电位
!
msfl（:
）:
微球型聚焦
!
lls（:
）:
浅侧向
!
lld（:
）:
深侧向
!
nphi（:
）:
补偿中子
!
gr（:
）:
自然伽马
!
rhob（:
）:
岩性密度
!
----------------------------------------
programcjsy
implicitnone
character*80cmdfile
character*80file_data,file_contant
character*80file_sp,file_msfl,file_lls,file_lld,file_nphi,file_gr,file_n_d
integern
realgrmin,grmax,spmin,spmax,b,r_sh,fai_nsh,p_sh,p_ma,p_f,fai_nf,fai_nma
real,allocatable:
:
sp（:
）,msfl（:
）,lls（:
）,lld（:
）,nphi（:
）,gr（:
）,rhob（:
）
cmdfile=''
callread_cmd（cmdfile,file_data,file_contant,file_sp,file_msfl,file_lls,file_lld,file_nphi,file_gr,file_n_d）
callread_number（file_data,n）
allocate（sp（1:
n）,msfl（1:
n）,lls（1:
n）,lld（1:
n）,nphi（1:
n）,gr（1:
n）,rhob（1:
n））
callread_data（file_data,n,sp,msfl,lls,lld,nphi,gr,rhob）
callread_contant（file_contant,grmin,grmax,spmin,spmax,b,r_sh,fai_nsh,p_sh,p_ma,p_f,fai_nf,fai_nma）
callget_sp（file_sp,n,sp,spmin,spmax）
callget_r（file_msfl,n,msfl,b,r_sh）
callget_r（file_lls,n,lls,b,r_sh）
callget_r（file_lld,n,lld,b,r_sh）
callget_nphi（file_nphi,n,nphi,fai_nsh）
callget_gr（file_gr,n,gr,grmin,grmax）
callget_n_d（file_n_d,n,nphi,rhob,fai_nsh,p_sh,p_ma,p_f,fai_nf,fai_nma）
deallocate（sp,msfl,lls,lld,nphi,gr,rhob）
end
!
------------------------------------------------------------------
!
!
功能：
读取参数文件
!
!
输入参数说明：
!
cmdfile:
参数文件
!
!
输出参数说明：
!
file_data:
输入测井数据文件
!
file_contant:
常数参量文件
!
file_sp:
输出自然电位估计的泥质含量文件
!
file_msfl:
输出微球型聚焦估计的泥质含量文件
!
file_lls:
输出浅侧向估计的泥质含量文件
!
file_lld:
输出深侧向估计的泥质含量文件
!
file_nphi:
输出补偿中子估计的泥质含量文件
!
file_gr:
输出自然伽马估计的泥质含量文件
!
file_n_d:
输出中子-密度交会图估计的泥质含量文件
!
------------------------------------------------------------------
subroutineread_cmd（cmdfile,file_data,file_contant,file_sp,file_msfl,file_lls,file_lld,file_nphi,file_gr,file_n_d）
implicitnone
character*（*）cmdfile
character*（*）file_data,file_contant
character*（*）file_sp,file_msfl,file_lls,file_lld,file_nphi,file_gr,file_n_d
open（10,file=cmdfile,status='old'）
read（10,*）file_data
read（10,*）file_contant
read（10,*）file_sp
read（10,*）file_msfl
read（10,*）file_lls
read（10,*）file_lld
read（10,*）file_nphi
read（10,*）file_gr
read（10,*）file_n_d
close（10）
endsubroutineread_cmd
!
------------------------------------------------------------------
!
!
功能：
读取测井数据的个数
!
!
输入参数说明：
!
filename:
输入测井数据文件
!
!
输出参数说明：
!
number:
测井数据的个数
!
------------------------------------------------------------------
subroutineread_number（filename,number）
implicitnone
character*（*）filename
integernumber
open（11,file=filename,status='old'）
number=0
dowhile（.（11））
read（11,*,end=100,ERR=100）
number=number+1
100enddo
close（11）
number=number-1
endsubroutineread_number
!
------------------------------------------------------------------
!
!
功能：
读取测井数据
!
!
输入参数说明：
!
file_data:
输入测井数据文件
!
n:
测井数据的个数
!
!
输出参数说明：
!
sp（:
）:
自然电位
!
msfl（:
）:
微球型聚焦
!
lls（:
）:
浅侧向
!
lld（:
）:
深侧向
!
nphi（:
）:
补偿中子
!
gr（:
）:
自然伽马
!
rhob（:
）:
岩性密度
!
------------------------------------------------------------------
subroutineread_data（file_data,n,sp,msfl,lls,lld,nphi,gr,rhob）
implicitnone
character*（*）file_data
integern
realsp（1:
n）,msfl（1:
n）,lls（1:
n）,lld（1:
n）,nphi（1:
n）,gr（1:
n）,rhob（1:
n）
integeri
realh
open（10,file=file_data,status='old'）
read（10,*）
doi=1,n,1
read（10,*）h,sp（i）,msfl（i）,lls（i）,lld（i）,rhob（i）,nphi（i）,gr（i）
enddo
close（10）
endsubroutineread_data
!
------------------------------------------------------------------
!
!
功能：
读取常数参量
!
!
输入参数说明：
!
file_contant:
常数参量文件
!
!
输出参数说明：
!
grmin,grmax:
自然伽马的最小值与最大值
!
spmin,spmax:
自然电位的最小值与最大值
!
b,r_sh:
已知常数
!
fai_nsh:
已知常数
!
p_sh,p_ma,p_f:
已知常数
!
fai_nf,fai_nma:
已知常数
!
------------------------------------------------------------------
subroutineread_contant（file_contant,grmin,grmax,spmin,spmax,b,r_sh,fai_nsh,p_sh,p_ma,p_f,fai_nf,fai_nma）
implicitnone
character*（*）file_contant
realgrmin,grmax,spmin,spmax,b,r_sh
升级会员 