[发明专利]一种LogR-ΔT测井评价的泥页岩油含量预测方法

权利要求书

1.一种LogR-ΔT测井评价的泥页岩油含量预测方法，其特征在于采用以下步骤：

步骤①、测井数据的预处理；

步骤②、极小值节点的识别；

步骤③、测井曲线拉平处理；

步骤④、标定预测泥页岩地层含油量S₁和氯仿沥青“A”的公式参数；

步骤⑤、定量评价非实测数据泥页岩地层页岩油S₁和氯仿沥青“A”含量，实现无取样泥页岩段地层页岩油含量S₁和氯仿沥青“A”连续预测评价。

2.根据权利要求1所述的一种LogR-ΔT测井评价的泥页岩油含量预测方法，其特征在于泥页岩地层含油量S₁的评价预测的步骤是：

步骤1、测井参数预处理：

①提取声波时差测井数据ΔT与电阻率测井数据R，ΔT的单位为毫秒每英尺(μs/ft),R的单位为欧姆米(Ω·m)，

②对电阻率R取对数处理得到logR；

步骤2、节点识别：

从对电阻率R测井曲线上和声波时差ΔT测井曲线上识别出局部范围的极小值，去除非含残留烃层段的影响；数学上求极小值通常用求导的方法，因此需对测井曲线求导；

第i个深度点的logR的一阶导数计算公式为：

logR′(d_i)=[logR(d_i+1)-logR(d_i)]/H              (8)

d_i为第i点的深度值，d_i的单位为米(m)；logR(d_i)为深度d_i处的logR，无量纲；logR(d_i+1)为深度d_i+1处的logR，无量纲；测井数据记录间隔为H(取0.125m)；

第i个深度点的声波时差值ΔT的一阶导数计算公式为：

ΔT′(d_i)=[ΔT(d_i+1)-ΔT(d_i)]/H            (9)

d_i为第i点的深度值，d_i的单位为米(m)；ΔT(d_i)为深度d_i处的声波时差值，ΔT(d_i)的单位为毫秒每英尺(μs/ft)；ΔT(d_i+1)为深度d_i+1处的声波时差值，ΔT(d_i+1)的单位为毫秒每英尺(μs/ft)；测井数据记录间隔为H(取0.125m)；

但由于测井曲线具有一定的采样间隔，在一定尺度上讲属不连续曲线，计算得到的极值点和实际的极值点不一定吻合，采用导数符号变化的点作为极值点；

logR′(d_i)<=0,logR′(d_i+1)>0;则logR(d_i)为极小值，记录d_m=d_i；

同样，如果ΔT′(d_i)<=0，ΔT′(d_i+1)>0;则ΔT(d_i)为极小值，记录d_n=d_i；

步骤3、测井曲线拉平：

在d_i深度段泥岩层中，d_m<d_i<d_m+1,则：

LlogR(d_i)=logR(d_i)-[logR(d_m)+logR(d_m-1)]/2      (10)

d_i为第i点的深度值，d_i的单位为米(m)；logR(d_i)为深度d_i处的logR，无量纲；logR(d_m)为深度d_m处的logR，无量纲；logR(d_m-1)为深度d_m-1处的logR，无量纲；LlogR(d_i)为拉平处理后的logR(d_i)值，无量纲；测井数据记录间隔为H(取0.125m)；

同样，在d_i深度段泥岩层中，d_n<d_i<d_n+1,则：

LΔT(d_i)=ΔT(d_i)-[ΔT(d_n)+ΔT(d_n-1)]/2        (11)

d_i为第i点的深度值，单位：米(m)；ΔT(d_i)为深度d_i处的声波时差值，单位为毫秒每英尺(μs/ft)；ΔT(d_n)为深度d_n处的声波时差值，单位为毫秒每英尺(μs/ft)；ΔT(d_n-1)为深度d_n-1处的声波时差值，单位为毫秒每英尺(μs/ft)；LΔT(d_i)为拉平处理后的声波时差值，单位为毫秒每英尺(μs/ft)；测井数据记录间隔为H(取0.125m)；

步骤4、含油量S₁的预测模型的标定：

S₁=a×LlogR+b×LΔT+c               (12)

对于一系列(n个点，n>=3)残留烃S₁的实测点：

都有对应的LlogR和LΔT值：

(LlogR_i,LΔT_i,S1_i)i=0,1,..,n-1

LlogR(_i)为第i个点拉平处理后的logR(d_i)值，无量纲；LΔT(_i)为第i个点拉平处理后的声波时差值，单位为毫秒每英尺(μs/ft)；S1_i为第i个点的热解含油量，单位为毫克每克(mg/g);

要使实测的含油量S₁与计算的含油量S₁最接近，即需要满足：

S=Σi=0n-1(a×LlogRi+b×LΔTi+c-Sli)2---(13)]]>

最小；

要使得S最小，应满足：

∂S∂a=0,∂S∂b=0,∂S∂c=0---(14)]]>

即：Σ2(a×LlogRi+b×LΔTi+c-Sli)LlogRi=0Σ2(a×LlogRi+b×LΔTi+c-Sli)LΔTi=0Σ2(a×LlogRi+b×LΔTi+c-Sli)=0---(15)]]>

有：

aΣ(LlogRi)2+bΣ(LlogRi)×(LΔTi)+cΣ(LlogRi)=Σ(LlogRi)×(Sli)bΣ(LlogRi)×(LΔTi)+bΣ(LΔTi)2+cΣ(LΔTi)=Σ(LΔTi)×(Sli)cΣ(LlogRi)+bΣ(LΔTi)+cn=Σ(Sli)---(16)]]>

解上述线形方程组，得：a，b，c;

即：S₁=a×LlogR+b×LΔT+c                 (12)

步骤5、模拟计算：

依据S₁=a×LlogR+b×LΔT+c，结合处理的测井数据推广计算无实测数据段的S₁值。