首页 / 专利库 / 地球科学 / 油页岩 / 干酪根 / 一种基于数理统计的镜质体反射率校准方法

一种基于数理统计的镜质体反射率校准方法

阅读:699发布:2020-08-02

专利汇可以提供一种基于数理统计的镜质体反射率校准方法专利检索,专利查询,专利分析的服务。并且一种基于数理统计的镜质体反射率校准方法:根据来自同一口井不同深度段的若干 烃 源岩样品的 岩石 热解 数据最高裂解 温度 推算对应样品的镜质体反射率;对每一个烃源岩样品,分别取镜质体反射率实测值与镜质体反射率值之差的绝对值,计算该组绝对值的算术平均值;根据算术平均值计算烃源岩样品的样本标准偏差或平均偏差;根据Chauvenet法则计算 阀 值;分别用不同烃源岩样品的镜质体反射率实测值减去阀值;将保留镜质体反射率实测值的烃源岩样品的最高裂解温度和镜质体反射率实测值拟合一次函数求出系数;求对应烃源岩样品校正后的镜质体反射率校正值;采用统计学方法检测镜质体反射率校正值的准确性。本 发明 对于单井和区域上的热史恢复具有良好的应用前景。,下面是一种基于数理统计的镜质体反射率校准方法专利的具体信息内容。

1.一种基于数理统计的镜质体反射率校准方法,其特征在于,包括如下步骤:
1)根据来自同一口井不同深度段的若干源岩样品的岩石热解数据最高裂解温度推算对应样品的镜质体反射率Ro;
2)对每一个烃源岩样品,分别取镜质体反射率实测值与步骤1)推算的镜质体反射率Ro值之差的绝对值,得到一组绝对值,计算该组绝对值的算术平均值;
3)根据算术平均值计算烃源岩样品的样本标准偏差或平均偏差;
4)根据Chauvenet法则计算值;
5)分别用不同烃源岩样品的镜质体反射率实测值Roi减去阀值Ф:
Roi-Ф(5),
将结果大于0的烃源岩样品的镜质体反射率实测值作异常值剔除,并同时剔除最高裂解温度,保留其余烃源岩样品的镜质体反射率实测值和最高裂解温度;
6)将保留镜质体反射率实测值的烃源岩样品的最高裂解温度作为x,将对应的镜质体反射率实测值作为y,拟合x对y的一次函数:
y=ax+b(6)
求出系数a、b;
7)将步骤5)所剔除烃源岩样品的最高裂解温度作为x′、将步骤6所得系数a、b分别代入如下公式:
y′=ax′+b(7)
所得y′值为对应烃源岩样品校正后的镜质体反射率校正值;
8)若烃源岩样品有最高裂解温度以外的热应数据,则将所述的最高裂解温度以外的热应力数据分别作为x和x′,重复步骤6)~步骤7)检验镜质体反射率校正值准确性,具体为将结果与步骤7)中利用最高裂解温度计算得到的镜质体反射率校正值比较,采用统计学方法检测镜质体反射率校正值的准确性。
2.根据权利要求1所述的一种基于数理统计的镜质体反射率校准方法,其特征在于,步骤1)所述的推算烃源岩样品的镜质体反射率实测值是采用如下公式:
Ro(Tmax)=0.0180×Tmax–7.16(1)
式中,Ro表示镜质体反射率实测值,Tmax为最高裂解温度,420℃干酪根,不适用于I型干酪根。
3.根据权利要求1所述的一种基于数理统计的镜质体反射率校准方法,其特征在于,步骤2)是采用如下公式得到平均值:
μ=∑ni=1[︱Roi–Roi(Tmax)︱]/n(2)
式中,n为样品总数。
4.根据权利要求1所述的一种基于数理统计的镜质体反射率校准方法,其特征在于,步骤3)是采用如下公式进行计算:
标准偏差σ={∑ni=1[︱Roi–Roi(Tmax)︱-μ]2/(n-1)}1/2,n>5(3)平均偏差δ=∑ni=1[︱Roi–Roi(Tmax)︱-μ]/n,n≤5(3’)
式中,Roi为第i个样品的镜质体反射率实测值,Roi(Tmax)为第i个样品的镜质体反射率推算值,n为样品总数。
5.根据权利要求1所述的一种基于数理统计的镜质体反射率校准方法,其特征在于,步骤4)是采用如下公式计算:
Ф=Zτσ  (4)
式中,Ф为阀值,Zτ为一个与样品总数n相联系的系数,Zτ由样本容量决定,通过查表获得,σ为样本标准偏差。
6.根据权利要求1所述的一种基于数理统计的镜质体反射率校准方法,其特征在于,步骤8)所述最高裂解温度以外的热应力参数是4/1-甲基二苯并噻吩转化率或甾烷异构化参数、或Ts/(Ts+Tm)或C3217α,21β(H)22S/(22S+22R)升藿烷异构化参数。

说明书全文

一种基于数理统计的镜质体反射率校准方法

技术领域

[0001] 本发明涉及一种镜质体反射率校准方法。特别是涉及一种基于数理统计的镜质体反射率校准方法。

背景技术

[0002] 地质学认为,镜质体是由高等植物的根、茎、叶在覆的还原条件下,经凝胶化作用而形成的。由于镜质体的芳构化是一个不可逆的过程,且转化率可以达到100%,因此与源岩的热演化阶段有很好的对应关系。镜质体反射率(Ro),是有机岩石学中根据镜质体的芳构化程度判断其成熟度的一项重要参数,是指在油浸条件下煤(烃源岩)中均质镜质体或基质镜质体抛光表面的反射光强度与垂直入射光强度之比,用百分数表示,其是煤岩学中确定煤阶的最重要方法。然而常规生油岩中的分散有机质很难像煤一样划分出各个亚组分,在基质镜质体中富氢镜质体和再循环镜质体是不能作为成熟度指标的,只有正常镜质体可以用来测定烃源岩镜质体反射率。在油田地质学中,I型和II1型干酪根基质镜质体的形成有低等生物和其他类脂化合物的参与,存在大量富氢镜质体和再循环镜质体,从而导致富氢。这类镜质体呈现明显的镜质组/类脂组过渡特征,它们的Ro实测值往往比估计的要低,抑制的主要原因在于泥炭化作用阶段腐殖体形成过程中类脂物的明显输入(朱抱荃,1987)。
[0003] 渤海湾盆地各富烃凹陷内烃源岩干酪根类型以II型为主,排烃限所对应的镜质体反射率大约为0.5%,但是盆地内各凹陷烃源岩镜质体反射率在进入古~始新统较高热演化阶段的层位后不再随埋深发生明显变化,实测干酪根镜质体反射率大多保持在0.5~1.0%之间,怀疑与上述镜质体表面富氢使其芳构化作用受到抑制有关。众所周知,镜质体反射率(Ro)是恢复盆地热史的重要参数,它的准确与否关系到烃源岩评价、资源量计算,与油气勘探息息相关,因此对实测值与预期相差甚远的可疑数据加以修正十分必要。

发明内容

[0004] 本发明所要解决的技术问题是,提供一种对于单井和区域上的热史恢复具有良好的应用前景的基于数理统计的镜质体反射率校准方法。
[0005] 本发明所采用的技术方案是:一种基于数理统计的镜质体反射率校准方法,包括如下步骤:
[0006] 1)根据来自同一口井不同深度段的若干烃源岩样品的岩石热解数据最高裂解温度推算对应样品的镜质体反射率Ro;
[0007] 2)对每一个烃源岩样品,分别取镜质体反射率实测值与步骤1)推算的镜质体反射率Ro值之差的绝对值,得到一组绝对值,计算该组绝对值的算术平均值;
[0008] 3)根据算术平均值计算烃源岩样品的样本标准偏差或平均偏差;
[0009] 4)根据Chauvenet法则计算值;
[0010] 5)分别用不同烃源岩样品的镜质体反射率实测值Roi减去阀值Ф:
[0011] Roi-Ф  (5),
[0012] 将结果大于0的烃源岩样品的镜质体反射率实测值作异常值剔除,并同时剔除最高裂解温度,保留其余烃源岩样品的镜质体反射率实测值和最高裂解温度;
[0013] 6)将保留镜质体反射率实测值的烃源岩样品的最高裂解温度作为x,将对应的镜质体反射率实测值作为y,拟合x对y的一次函数:
[0014] y=ax+b  (6)
[0015] 求出系数a、b;
[0016] 7)将步骤5)所剔除烃源岩样品的最高裂解温度作为x′、将步骤6所得系数a、b分别代入如下公式:
[0017] y′=ax′+b  (7)
[0018] 所得y′值为对应烃源岩样品校正后的镜质体反射率校正值;
[0019] 8)若烃源岩样品有最高裂解温度以外的热应数据,则将所述的最高裂解温度以外的热应力数据分别作为x和x′,重复步骤6)~步骤7)检验镜质体反射率校正值准确性,具体为将结果与步骤7)中利用最高裂解温度计算得到的镜质体反射率校正值比较,采用统计学方法检测镜质体反射率校正值的准确性。
[0020] 步骤1)所述的推算烃源岩样品的镜质体反射率实测值是采用如下公式:
[0021] Ro(Tmax)=0.0180×Tmax–7.16(1)
[0022] 式中,Ro表示镜质体反射率实测值,Tmax为最高裂解温度,420℃
[0023] 步骤2)是采用如下公式得到平均值:
[0024] μ=∑ni=1[︱Roi–Roi(Tmax)︱]/n  (2)
[0025] 式中,n为样品总数。
[0026] 步骤3)是采用如下公式进行计算:
[0027] 标准偏差σ={∑ni=1[︱Roi–Roi(Tmax)︱-μ]2/(n-1)}1/2,n>5  (3)[0028] 平均偏差δ=∑ni=1[︱Roi–Roi(Tmax)︱-μ]/n,n≤5  (3’)
[0029] 式中,Roi为第i个样品的镜质体反射率实测值,Roi(Tmax)为第i个样品的镜质体反射率推算值,n为样品总数。
[0030] 步骤4)是采用如下公式计算:
[0031] Ф=Zτσ  (4)
[0032] 式中,Ф为阀值,Zτ为一个与样品总数n相联系的系数,Zτ由样本容量决定,通过查表获得,σ为样本标准偏差。
[0033] 步骤8)所述最高裂解温度以外的热应力参数是4/1-甲基二苯并噻吩转化率或甾烷异构化参数、或Ts/(Ts+Tm)或C3217α,21β(H)22S/(22S+22R)升藿烷异构化参数。
[0034] 本发明的一种基于数理统计的镜质体反射率校准方法,对于拥有纵向上连续的烃源岩热解数据或其它地化数据的单井十分有效,通常一口单井的烃源岩镜质体发射率抑制作用仅发生在某一深度段,如烃源岩热演化较高的,或干酪根类型偏腐泥型的层位,并非所有烃源岩样品的实测Ro值都不准,通过实测Ro值与其它热应力参数转化的Ro值互检的手段可以找到这些异常点,有针对性的进行纠正,且不过于依赖经验公式,而是根据每口井自身的情况建立独立的校准公式,对于单井和区域上的热史恢复具有良好的应用前景。附图说明
[0035] 图1是本发明一种基于数理统计的镜质体反射率校准方法的流程图
[0036] 图2是剔除异常数据后所展示的Tmax对Ro(实测)的拟合结果;
[0037] 图3是将所剔除样品的Tmax值带回校正公式后的计算结果;
[0038] 图4是校正前样品实测Ro与对应氢指数(HI)随深度的变化趋势;
[0039] 图5是样品校正前的Ro与校正后的Ro数据对比以及随深度的变化趋势;
[0040] 图6是样品校正前的Ro与校正后的Ro数据放大后与氢指数随深度的变化趋势;
[0041] 图7是将Tmax改为4-甲基二苯并噻吩转化率(MDR-4)后对Ro(实测)的拟合结果;
[0042] 图8是将所剔除样品的MDR-4值带回校正公式后的计算结果。

具体实施方式

[0043] 下面结合实施例和附图对本发明的一种基于数理统计的镜质体反射率校准方法做出详细说明。
[0044] 如图1所示,本发明是一种基于数理统计的镜质体反射率校准方法,包括如下步骤:
[0045] 1)根据来自同一口井不同深度段的若干烃源岩样品的岩石热解数据(Rock-Eval)最高裂解温度(Tmax)推算对应样品的镜质体反射率;
[0046] 所述的推算烃源岩样品的镜质体反射率实测值是采用如下公式:
[0047] Ro(Tmax)=0.0180×Tmax–7.16  (1)
[0048] 式中,Ro表示镜质体反射率实测值,Tmax为最高裂解温度,420℃<Tmax<500℃该公式可较好地适用于许多II型和III型干酪根,但不适用于I型干酪根。
[0049] 2)对每一个烃源岩样品,分别取镜质体反射率实测值与步骤1)推算的镜质体反射率Ro(Tmax)值之差的绝对值,得到一组绝对值,计算该组绝对值的算术平均值;是采用如下公式得到平均值:
[0050] μ=∑ni=1[︱Roi–Roi(Tmax)︱]/n  (2)
[0051] 式中,n为样品总数。
[0052] 3)根据算术平均值计算烃源岩样品的样本标准偏差或平均偏差;是采用如下公式进行计算:
[0053] 标准偏差σ={∑ni=1[︱Roi–Roi(Tmax)︱-μ]2/(n-1)}1/2,n>5  (3)[0054] 平均偏差δ=∑ni=1[︱Roi–Roi(Tmax)︱-μ]/n,n≤5  (3’)
[0055] 式中,Roi为第i个样品的镜质体反射率实测值,Roi(Tmax)为第i个样品的镜质体反射率推算值,n为样品总数。
[0056] 4)根据Chauvenet法则计算阀值;是采用如下公式计算:
[0057] Ф=Zτσ  (4),
[0058] 式中,Ф为阀值,Zτ为一个与样品总数n相联系的系数,Zτ由样本容量决定,通过查表获得,σ为样本标准偏差。
[0059] 5)分别用不同烃源岩样品的镜质体反射率实测值Roi减去阀值Ф:
[0060] Roi-Ф  (5)
[0061] 将结果大于0的烃源岩样品的镜质体反射率实测值作异常值剔除,并同时剔除最高裂解温度,保留其余烃源岩样品的镜质体反射率实测值和最高裂解温度;
[0062] 6)将保留镜质体反射率实测值的烃源岩样品的最高裂解温度作为x,将对应的镜质体反射率实测值作为y,拟合x对y的一次函数:
[0063] y=ax+b  (6)
[0064] 求出系数a、b;
[0065] 7)将步骤5)所剔除烃源岩样品的最高裂解温度作为x′、将步骤6所得系数a、b分别代入如下公式:
[0066] y′=ax′+b  (7)
[0067] 所得y′值为对应烃源岩样品校正后的镜质体反射率校正值;
[0068] 8)若烃源岩样品有最高裂解温度以外的热应力数据,则将所述的最高裂解温度以外的热应力数据分别作为x和x′,重复步骤6)~步骤7)检验镜质体反射率校正值准确性,具体为将结果与步骤7)中利用最高裂解温度计算得到的镜质体反射率校正值比较,采用统计学方法检测镜质体反射率校正值的准确性。所述最高裂解温度以外的热应力参数是4/1-甲基二苯并噻吩转化率(MDR)或甾烷异构化参数、或Ts/(Ts+Tm)或C3217α,21β(H)22S/(22S+22R)升藿烷异构化参数。
[0069] 下面结合附图和表格以QHD29-2-1井为例进行具体说明。
[0070] 1)根据岩石热解数据(Rock-Eval)中的最高裂解温度(Tmax)推算烃源岩样品的镜质体反射率值,Tmax-Ro转换公式为:Ro=Tmax*0.018-7.16(Teichmuller and Durand,1983),但注意事先剔除不适用于推导公式的数据。如表1第二、三列所示,可以看到部分推算镜质体反射率与实测镜质体反射率(表1第五列所示)有较大差距。
[0071] 2)用推算镜质体反射率检测实测镜质体反射率,本技术的关键在于将镜质体反射率的实测值与推算值之差的绝对值(如表1第六列所示)作为一系列独立的数据,计算︱Ro(实测)–Ro(Tmax)︱的算术平均值μ如表1第九列所示,计算︱实测镜质体反射率–推算镜质体反射率︱的标准偏差σ或平均偏差δ如表1第十列所示,根据Chauvenet法则计算阀值Ф=Zτσ,因所举实例中样品数为36,所以Zτ=2.39,Ф=2.39σ,计算结果如表1第十一列所示,该阀值相当于以真值为圆心,以阀值为半径画圆,若镜质体反射率的实测值与推算值之差在圆圈范围内(即小于Ф),则认为该实测镜质体反射率值可靠,予以保留,否则认为该实测镜质体反射率不可靠,予以剔除。
[0072] 3)剔除实测镜质体反射率与推算镜质体反射率差距过大(高于阀值)的样品数据点,即︱实测镜质体反射率–推算镜质体反射率︱–Ф≥0的点如表1第八列带*的数据所示。
[0073] 4)剩余样品所展示的最高裂解温度对实测镜质体反射率的变化规律被认为是可信的,以最高裂解温度为“x”,以实测镜质体反射率为“y”,拟合该井最高裂解温度与实测镜质体反射率的相关公式,如图2所示,结果为:y=0.0088x+3.2262,即:a=0.0088,b=3.2262。
[0074] 5)将所剔除样品的最高裂解温度值(Tmax)带回上一步所求得公式,计算出校正后的镜质体反射率,结果如图3、表1第十二列带*数据所示,如图4所示,校正前的镜质体反射率受到显著的抑制作用,表现为进入氢指数较高(≥350mg/g)的烃源岩层位后实测镜质体反射率不再随深度发生显著变化,如图5所示,校正镜质体反射率,要比实测镜质体发射率高出0.1-0.2%,且校正后的镜质体反射率在氢指数较高(≥350mg/g)的烃源岩层位仍随埋深升高(如图6所示),明显克服了因镜质体芳构化进程受阻对反射率的抑制作用。
[0075] 6)如表1第七列所示,将岩石最高裂解温度(Tmax)改为4-甲基二苯并噻吩转化率(MDR-4),该参数与镜质体反射率(Ro)、岩石最高热裂解温度(Tmax)同为热应力参数,如图7、图8所示,重复步骤(4)-(5),检验校正镜质体反射率的准确性。由于这类数据不多,导致将最高热裂解温度改为4-甲基二苯并噻吩转化率后样本容量骤减,但使用4-甲基二苯并噻吩转化率对镜质体反射率的校正结果与使用最高裂解温度比较接近,证明这种基于数理统计的镜质体反射率校准方法对该井的受抑制实测镜质体反射率的改善是可行的。
[0076] 上述实施例仅是用于说明本发明的技术原理和实际应用,其中方法的各实施步骤都是可以修改的,凡是在本发明技术方案的基础上进行的变化或改进,均在本发明的保护范围之内。
[0077] 表1
[0078]
[0079]
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈