复杂岩性交互处理程序

复杂岩性交互处理程序
参数选择。

在程序处理过程中，由参数控制处理层段的岩性，比如是砂泥岩地层或者是碳酸盐岩地层，当YAFG=1，2时，为进行砂泥岩处理，当YAFG＝3时，进行碳酸盐岩处理。

如图分别为砂泥岩剖面和碳酸盐岩剖面。

3．输入输出曲线及参数说明
输入曲线：
AC=声波时差
CAL=井径
CNL=补偿中子
COND=感应电导
DEN=体积密度
GR=自然伽玛
K=钾
MSFL=微球聚焦
PE=光电吸收截面指数
PORS=井壁中子
RA1=电阻率？
RILD=深感应电阻率
RLL8=八侧向电阻率
RLLD=深侧向电阻率
RLLS=浅侧向电阻率
RXO=冲洗带电阻率
SP=自然电位
TC=能谱测井总计数率
TH=钍
TPI=钍-钾指数TPI=(K+3.1)*(TH+12.4)*0.1
输出曲线：
LIME=石灰岩体积
ANHY=硬石膏体积
BDEN=(PORD-POR2)/PORD
BPOR=(PORA-PORC)/PORA
C1=附加矿物1的体积（盐岩）
CARB=煤层体积
CENG=层号曲线
DOLO=白云岩体积
DREC=
EF=杨氏模量
G=
PERM=渗透率
PG=含气概率
PGW=气水概率
POR=有效孔隙度
POR2=次生孔隙度(POR-PORA)
PORB=
PORE=
PORR=
PORT=总孔隙度
PORW=地层含水孔隙度
RT=地层真电阻率
SAND=砂岩体积
SH=泥质含量
SW=地层含水饱和度
JIAO=交会法判别参数
其中：DREC=(DENW*AC2)/(ACW2*DEN)-1
DENW=CNL*DF+(1-CNL)*2.67 ACW=CNL*TF+(1-CNL)*TM
常用参数
SMN1=GR在纯地层的极小值
SMX1=GR在泥岩层的极大值
SMN2=TC在纯地层的极小值
SMX2=TC在泥岩层的极大值
SMN3=K在纯地层的极小值
SMX3=K在泥岩层的极大值
SMN4=TH在纯地层的极小值
SMX4=TH在泥岩层的极大值
SMN5=TPI在纯地层的极小值
SMX5=TPI在泥岩层的极大值
SMN6=SP在纯地层的极小值
SMX6=SP在泥岩层的极大值
SMN7=AC在纯地层的极小值
SMX7=AC在泥岩层的极大值
SMN8=CNL在纯地层的极小值
SMX8=CNL在泥岩层的极大值
GC1=计算泥质含量用的经验指数-GR
GC2=计算泥质含量用的经验指数-TC
GC3=计算泥质含量用的经验指数-K40
GC4=计算泥质含量用的经验指数-TH
GC5=计算泥质含量用的经验指数-TPI
RW=地层水电阻率
RMF=泥浆滤液电阻率
SHFG=泥质体积计算方法标志
TSH=泥岩的声波时差
TF=流体的声波时差值
TM=骨架的声波时差
AAC=AC的附加校正值
TFM=井眼未垮塌的纯地层的最大声波时差值 TMSH=干粘土骨架声波时差值
BITS=钻头直径
ACAL=CAL的附加校正值
ANSH=泥岩的中子值
DSH=泥岩的密度
DF=密度的流体值
DG=密度的骨架
DFM=井眼未垮塌的纯地层的最小密度值
SWOP=选择含水饱和度计算公式的标志
A=岩性系数
AM=胶结指数
AN=饱和度指数
ADEN=密度附加校正值
PRFG=YAFG为3时求孔隙度和矿物成分旗标
ANEU=中子附加校正值
ART=电阻率附加校正值
ARXO=RXO附加校正值
AMNE=NEU的乘法校正值
ACON=COND的附加校正值
RSH=泥岩的电阻率
RLIM=纯砂岩最大电阻率
SHCT=泥岩截止值
AM1X=附加矿物1在X的骨架值
AM1Y=附加矿物1在Y的骨架值
AM2X=附加矿物2在X的骨架值
AM2Y=附加矿物2在Y的骨架值
AM3X=附加矿物3在X的骨架值
AM3Y=附加矿物3在Y的骨架值
AM4X=附加矿物4在X的骨架值
AM4Y=附加矿物4在Y的骨架值
SAN=砂岩存在标志1-有砂岩0-无砂岩
CLIM=石灰岩存在标志1-有石灰岩0-无石灰岩 DOL=白云岩存在标志1-有白云岩0-无白云岩 ANH=石膏存在的标志1-有石膏0-无石膏
RDEP=选择RW处的深度
GRAD=地区的地温梯度
BHT=井底温度
SIRR=束缚水饱和度
SMF=泥浆滤液矿化度
DHC=油气密度
SXOL=计算SXO指数上限值
PUSH=对中子或密度或声波校正加权系数
DGMN=最小的骨架密度值
TMMX=最大的骨架声波时差值
CP=压实系数
ANFM=井眼未垮塌的纯地层的最大中子值
YAFG=地层旗标1和2-砂泥岩地层3-碳酸盐岩地层 PFG1=计算POR方法旗标1-经验系数2-常规公式 GA=当PFG和PFG1为1时求POR经验公式系数G YA=当PFG和PFG1为1时求POR经验公式系数Y POCT=分层时POR最小值
CACT=井径最大值
COAC=煤层的AC最小截至值
PFG=求POR旗标1-AC2-DEN3-CNL
GC=砂泥岩地层中子岩性校正值
GD=当PFG为2PFG1为1时求POR经验公式系数G YD=当PFG为2PFG1为1时求POR经验公式系数Y ADCF=当YAFG为2时求POR旗标
RLCT=用RLLD分层的参数
GRCT=用GR分层的参数
ACCT=用AC分层的参数
CNCT=用CNL分层的参数
C=求RT旗标1为RLLD2为RILD3为COND
FEN1=分层参数
FEN2=
ANG=中子骨架值
RXOF=RXO旗标1输入2为RLL83为RLLS4为MSFL SFG=砂泥岩剖面求PERM旗标
UF=
DOFG=
USH=
APE=
PFM=PE曲线下限截至值
CNFG=中子井眼校正旗标1校正0不校正
AMF=
A10=马五1（深度<A10）
A20=马五2（A10<深度<A20）
A30=马五3（A20<深度<A30）马五4（深度>A30） PEFG=
Z=
AMFG=
TX=分层处理旗标0-不分层1-分层
FOUT=矿物体积含量累加输出或单独输出旗标
4．旗标及公式
SHFG 1－GR,2－TC,3－AK,4－TH,5－TPI,6－SP,7－AC,8－CNL,9－RT,10－CNL-DEN,11－AC-DEN,12－AC-CNL
SWOP 1－N T
M W R POR R A SW /1)*(⋅=, 2－POR SH
SH R R SW T W *175.0*036.0/2-+=
3－利用计算M,1,/019.08.1=+=A POR M ,
若44,1.20825.0=>=>M M M POR 则若则,然后用公式1计算
4－22211221*5**4.0,,T T R R T SW RSH SH T POR
R T T W -+===
YAFG 1－用AC,CNL,DEN 单条曲线求砂泥岩地层孔隙度
2－用两条曲线求砂泥岩地层孔隙度
3－处理碳酸盐岩地层
C 1－RLL
D R T =
2－RILD R T =
3－COND R T /1000=
当YAFG=1时，
PFG = 1－用AC 求孔隙度；若PFG1=1，则YA GA CP AC POR /)/(-=;若PFG1=2，则)*)/(()(CP TM TF TM AC POR --=
2－用DEN 求孔隙度；若PFG1=1，则YD DEN GD POR /)(-=;若PFG1=2，则)/()(DF DG DEN DG POR --=
3－用CNL 求孔隙度；
为百分数，其中NSH GC CNL SH NSH GC CNL POR ,,01.0**01.0*)(-+=
当YAFG=2时，
ADCF = 1－2
2
2PORC PORA POR += 2－22
2PORC PORD POR +=
3－2/))2
((222
2PORC PORC PORA POR ++= 输出SH POR POR POR *-=
当YAFG=1,2时
若SFG ＝ 1－为小数POR SIRR
POR PERM ,*136.024
.4= 2－
221*875.31**238.0*033.01*005.4*471.0432.2SH SH POS POS SH POS PERM +++--=
其中1
11136.29*163.0SMN SMX SMN GR SH AC POS --=-=， 3－)1*6579.1)(*9908.32102.9*(*10SH AC LOG PERM -+-=
当YAFG=3时(注：PRFG 只有在YAFG=3时使用)
若PRFG = 1－用DEN,CNL 交会求POR 和矿物体积
2－用AC,CNL 交会求POR 和矿物体积
3－用DEN 求POR 和矿物体积
4－用AC 求POR 和矿物体积
5－用CNL 求POR 和矿物体积
6－用PE 求POR 和矿物体积
当yafg ＝4时，或swop=7时，
若amfg ＝1，
amfg ＝2时 m ＝1.94＋0.0062/por
amfg=3时，m ＝1.87+0.019/por
)])1((/[)(Rw Z Rmf Z Rt Rw Rmf A SW m -+⋅⋅⋅⋅⋅=φ
马五1，4
947
.2547.146-=AC φ 马五2 888
.3146-=AC φ 马五3 08.528.150*0078.0-+=
RLLD AC φ )])1((/[)(Rw Z Rmf Z Rt Rw Rmf SW m -+⋅⋅⋅⋅=φ n w b
S A K φ⋅=
马五1：a ＝0.042，b=4.2，n=2
马五2：a=0.015，b=3.46，n=2
马五3：a=0.024，b=4.3，n=1.64。