CN109710887B - 一种基于回归分析的地化录井烃类损失校正方法 - Google Patents
一种基于回归分析的地化录井烃类损失校正方法 Download PDFInfo
- Publication number
- CN109710887B CN109710887B CN201811320271.XA CN201811320271A CN109710887B CN 109710887 B CN109710887 B CN 109710887B CN 201811320271 A CN201811320271 A CN 201811320271A CN 109710887 B CN109710887 B CN 109710887B
- Authority
- CN
- China
- Prior art keywords
- core
- wall
- rock
- rock debris
- hydrocarbon
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analysing Biological Materials (AREA)
Abstract
本发明公开一种基于回归分析的地化录井烃类损失校正方法,步骤1,对岩屑与壁心的热解参数通过散点图进行相关性分析,优选拟合函数模型;步骤2,对壁心与岩心的热解参数通过散点图进行相关性分析,优选拟合函数模型等步骤;该方法通过使用基于最小二乘法和最小均方算法的多元线性回归分析及非线性回归分析方法,建立岩屑与壁心、壁心与岩心的关系式,从而建立岩屑与岩心对应热解参数的烃损恢复函数关系式,最终得到的岩心的烃类含量即为地层原始烃类含量,使烃类损失恢复校正更全面准确和科学合理。
Description
技术领域
本发明涉及地化录井烃类损失校正方法,特别是一种基于回归分析的地化录井烃类损失校正方法,适用于岩屑录井。
背景技术
地化录井主要包括岩屑录井、井壁取心及钻井取心三种取样方式,其中岩屑录井最为实用以及成本低,同时岩屑在返回井口的过程中受钻井液冲刷等因素影响,烃类损失最大;井壁取心烃类损失程度次之;钻井取心烃类损失最少,最为接近原始地层烃类含量,但钻井取心成本最高。为了节约成本,需要建立一种模型,岩屑热解数据通过此模型计算可以大致恢复地层原始烃类含量。
目前,对地化录井烃类损失校正的研究,大多以经验图版、临界点分析及一元线性回归的方法为主,这些方法多受人为因素影响较大,校正效果不好。
发明内容
本发明的目的是为了提供一种更为合理有效的地化录井烃类损失校正方法,本发明使用基于最小二乘法和最小均方算法的多元线性回归和非线性回归方法对烃类损失恢复校正,与现有一元线性回归的方法相比,对地层原始烃类含量的数据恢复更科学准确。
本发明的技术方案是:
一种基于回归分析的地化录井烃类损失校正方法,其特征在于,包括如下步骤:
步骤1,分别以岩屑和壁心热解参数为自变量和因变量,对岩屑与壁心的热解参数通过散点图进行相关性分析,优选拟合函数模型;
步骤2,分别以壁心和岩心热解参数为自变量和因变量,对壁心与岩心的热解参数通过散点图进行相关性分析,优选拟合函数模型;
步骤3,通过多元线性回归分析或非线性回归分析方法,建立岩屑到壁心的对应热解参数的烃损恢复函数关系式;
步骤4,通过多元线性回归分析或非线性回归分析方法,建立壁心到岩心的对应热解参数的烃损恢复函数关系式;
步骤5,根据步骤3、步骤4得到的关系式,建立岩屑与岩心的对应热解参数的烃损恢复函数关系式;
步骤6,对未校正岩屑的热解参数值通过步骤5得到的函数式进行计算,得到岩心的烃类含量即为地层原始烃类含量。
上述的基于回归分析的地化录井烃类损失校正方法,所述岩屑与壁心的热解参数散点关系或壁心与岩心的热解参数散点关系中,若因变量与自变量散点图呈线性关系时,采用多元线性回归表达式作为烃损恢复函数关系式;若因变量与自变量散点图呈非线性关系时,选取多个非线性函数模型并计算每个非线性函数模型的相关系数及残差平方和,选择相关系数最大及残差平方和最小的函数模型作为烃损恢复函数关系式。
上述的基于回归分析的地化录井烃类损失校正方法,以热解参数-热解烃含量S2为例,进行岩屑到壁心校正时:若因变量与自变量散点图呈线性关系,使用多元线性回归分析,多元线性回归表达式为:
S2(壁心)=a*S1(岩屑)+b*S2(岩屑)+c*pg(岩屑)+d*opi(岩屑)+e*tpi(岩屑)+f
式中:S2(壁心)为壁心热解烃含量(单位mg/g);S1(岩屑)为岩屑可溶烃含量(单位mg/g);S2(岩屑)为岩屑热解烃含量(单位mg/g);pg(岩屑)为岩屑产油潜量(单位mg/g);opi(岩屑)为岩屑油产率指数;tpi(岩屑)为岩屑油气产率指数;a、b、c、d、e、f为利用样本集求解的未知系数。
上述的基于回归分析的地化录井烃类损失校正方法,以热解参数-热解烃含量S2为例,进行岩屑到壁心校正时:若因变量与自变量散点图呈非线性关系,选取y=x/(a*x+b)、y=a*(1-exp(-b*x))、y=a*exp(-b/x)以及y=a+b*log(x)等非线性函数模型,然后计算每个非线性函数模型的相关系数及残差平方和,选择相关系数最大及残差平方和最小的函数模型,最后得到烃损恢复非线性表达式,当使用双曲线y=x/(a*x+b)模型进行非线性回归时,非线性回归表达式为:
S2(壁心)=S2(岩屑)/(a*S2(岩屑)+b)
式中,S2(壁心)为壁心热解烃含量(单位mg/g);S2(岩屑)为岩屑热解烃含量(单位mg/g);a、b为利用样本集求解的未知系数。
上述的基于回归分析的地化录井烃类损失校正方法,以热解参数-热解烃含量S2为例,进行壁心到岩心校正时:若因变量与自变量散点图呈线性关系,使用多元线性回归分析,多元线性回归表达式为:
S2(岩心)=a′*S1(壁心)+b′*S2(壁心)+c′*pg(壁心)+d′*opi(壁心)+e′*tpi(壁心)+f′
式中:S2(岩心)为岩心热解烃含量(单位mg/g);S1(壁心)为壁心可溶烃含量(单位mg/g);S2(壁心)为壁心热解烃含量(单位mg/g);pg(壁心)为壁心产油潜量(单位mg/g);opi(壁心)为壁心油产率指数;tpi(壁心)为壁心油气产率指数;a′、b′、c′、d′、e′、f′为利用样本集求解的未知系数。
上述的基于回归分析的地化录井烃类损失校正方法,以热解参数-热解烃含量S2为例,进行壁心到岩心校正时:若因变量与自变量散点图呈非线性关系,选取y=x/(a′*x+b′)、y=a′*(1-exp(-b′*x))y=a′*exp(-b′/x)以及y=a′+b′*log(x)非线性函数模型,然后计算每个非线性函数模型的相关系数及残差平方和,选择相关系数最大及残差平方和最小的函数模型,最后得到烃损恢复非线性表达式,当使用双曲线y=x/(a′*x+b′)模型进行非线性回归时,非线性回归表达式为:
S2(岩心)=S2(壁心)/(a′*S2(壁心)+b′)
式中,S2(岩心)为岩心热解烃含量(单位mg/g);S2(壁心)为壁心热解烃含量(单位mg/g);a′、b′为利用样本集求解的未知系数。
本发明的有益效果是:
1、本发明按照恢复模型建立岩屑与壁心关系式,岩心与壁心关系式,联结二者关系式,形成岩屑与岩心关系式,建立关系式时采用数据挖掘中的回归分析方法,通过使用基于最小二乘法和最小均方算法的多元线性回归分析及非线性回归分析方法,建立岩屑与壁心、壁心与岩心的关系式,从而建立岩屑与岩心对应热解参数的烃损恢复函数关系式,最终得到的岩心的烃类含量即为地层原始烃类含量,实现烃类损失恢复校正,更加科学合理。
2、以往学者常用经验图版法、一元线性回归等来对烃类损失进行恢复,这样受人为影响因素较大,同时一元线性回归分析方法较为简单不全面,使用本方法建立的烃类损失恢复模型,恢复效果更加全面准确。
附图说明
图1是某地区岩屑与壁心S2散点图(实施例1);
图2是某地区岩屑与壁心S2的拟合值与实测值线性关系图,图中圆点为实测值,直线为拟合曲线(实施例1);
图3是某地区岩屑与壁心S2散点图(实施例2);
图4是某地区岩屑与壁心S2的曲线拟合图,图中圆圈为实测值,曲线为拟合曲线(实施例2)。
具体实施方式
实施例1
本实施例中,分别以岩屑和壁心热解参数为自变量和因变量,对岩屑与壁心的热解参数通过散点图进行相关性分析,由图1可知,因变量与自变量呈线性关系,所以采用多元线性回归表达式作为烃损恢复函数关系式。以热解参数-热解烃含量S2为例,进行岩屑到壁心校正时,多元线性回归表达式为:
S2(壁心)=a*S1(岩屑)+b*S2(岩屑)+c*pg(岩屑)+d*opi(岩屑)+e*tpi(岩屑)+f
式中:S2(壁心)为壁心热解烃含量(单位mg/g);S1(岩屑)为岩屑可溶烃含量(单位mg/g);S2(岩屑)为岩屑热解烃含量(单位mg/g);pg(岩屑)为岩屑产油潜量(单位mg/g);opi(岩屑)为岩屑油产率指数;tpi(岩屑)为岩屑油气产率指数;a、b、c、d、e、f为利用样本集求解的未知系数,根据实测数据计算得到图2中拟合直线,其中,a=3.8044、b=4.0407、c=3.4821、d=21.1362、e=20.0501、f=1.1071。
本实施例中,壁心到岩心的校正方法与岩屑到壁心的校正方法一样。以热解参数-热解烃含量S2为例,进行壁心到岩心校正时:若因变量与自变量散点图呈线性关系,使用多元线性回归分析,多元线性回归表达式为:
S2(岩心)=a′*S1(壁心)+b′*S2(壁心)+c′*pg(壁心)+d′*opi(壁心)+e′*tpi(壁心)+f′
式中:S2(岩心)为岩心热解烃含量(单位mg/g);S1(壁心)为壁心可溶烃含量(单位mg/g);S2(壁心)为壁心热解烃含量(单位mg/g);pg(壁心)为壁心产油潜量(单位mg/g);opi(壁心)为壁心油产率指数;tpi(壁心)为壁心油气产率指数;a′、b′、c′、d′、e′、f′为利用样本集求解的未知系数。
实施例2
本实施例中,分别以岩屑和壁心热解参数为自变量和因变量,对岩屑与壁心的热解参数通过散点图进行相关性分析,由图3可知,因变量与自变量呈非线性关系,选取多个非线性函数模型并计算每个非线性函数模型的相关系数及残差平方和,选择相关系数最大及残差平方和最小的函数模型作为烃损恢复函数关系式。以热解参数-热解烃含量S2为例,进行岩屑到壁心校正时,选取y=x/(a*x+b)、y=a*(1-exp(-b*x))、y=a*exp(-b/x)以及y=a+b*log(x)非线性函数模型,然后计算每个非线性函数模型的相关系数及残差平方和,选择相关系数最大及残差平方和最小的函数模型,最后得到烃损恢复非线性表达式,本实施例中使用双曲线y=x/(a*x+b)模型进行非线性回归,非线性回归表达式为:
S2(壁心)=S2(岩屑)/(α*S2(岩屑)+b)
式中,S2(壁心)为壁心热解烃含量(单位mg/g);S2(岩屑)为岩屑热解烃含量(单位mg/g);a、b为利用样本集求解的未知系数,根据图3中散值分布,计算得到图4中拟合曲线,其中,a=0.0318、b=0.0737。
本实施例中,壁心到岩心的校正方法与岩屑到壁心的校正方法一样。以热解参数-热解烃含量S2为例,进行壁心到岩心校正时:若因变量与自变量散点图呈非线性关系,选取y=x/(a′*x+b′)、y=a′*(1-exp(-b′*x))、y=a′*exp(-b′/x)以及y=a′+b′*log(x)非线性函数模型,然后计算每个非线性函数模型的相关系数及残差平方和,选择相关系数最大及残差平方和最小的函数模型,最后得到烃损恢复非线性表达式,当使用双曲线y=x/(a′*x+b′)模型进行非线性回归时,非线性回归表达式为:
S2(岩心)=S2(壁心)/(a′*S2(壁心)+b′)
式中,S2(岩心)为岩心热解烃含量(单位mg/g);S2(壁心)为壁心热解烃含量(单位mg/g);a′、b′为利用样本集求解的未知系数。
Claims (5)
1.一种基于回归分析的地化录井烃类损失校正方法,其特征在于,包括如下步骤:
步骤1,分别以岩屑和壁心热解参数为自变量和因变量,对岩屑与壁心的热解参数通过散点图进行相关性分析,优选拟合函数模型;
步骤2,分别以壁心和岩心热解参数为自变量和因变量,对壁心与岩心的热解参数通过散点图进行相关性分析,优选拟合函数模型;其中:
所述优选拟合函数模型进一步包括,若散点图呈线性关系时,采用多元线性回归表达式作为烃损恢复函数关系式;若呈非线性关系时,选取多个非线性函数模型并计算每个非线性函数模型的相关系数及残差平方和,选择相关系数最大及残差平方和最小的函数模型作为烃损恢复函数关系式;
步骤3,通过多元线性回归分析或非线性回归分析方法,建立岩屑到壁心的对应热解参数的烃损恢复函数关系式;
步骤4,通过多元线性回归分析或非线性回归分析方法,建立壁心到岩心的对应热解参数的烃损恢复函数关系式;
步骤5,根据步骤3、步骤4得到的关系式,建立岩屑与岩心的对应热解参数的烃损恢复函数关系式;
步骤6,对未校正岩屑的热解参数值通过步骤5得到的函数式进行计算,得到岩心的烃类含量即为地层原始烃类含量。
2.根据权利要求1所述的一种基于回归分析的地化录井烃类损失校正方法,其特征在于:所述步骤3中岩屑到壁心校正时:若因变量与自变量散点图呈线性关系,使用多元线性回归分析,多元线性回归表达式为:
S2(壁心)=a*S1(岩屑)+b*S2(岩屑)+c*pg(岩屑)+d*opi(岩屑)+e*tpi(岩屑)+f
式中:S2(壁心)为壁心热解烃含量(单位mg/g);S1(岩屑)为岩屑可溶烃含量(单位mg/g);S2(岩屑)为岩屑热解烃含量(单位mg/g);pg(岩屑)为岩屑产油潜量(单位mg/g);opi(岩屑)为岩屑油产率指数;tpi(岩屑)为岩屑油气产率指数;a、b、c、d、e、f为利用样本集求解的未知系数。
3.根据权利要求1所述的一种基于回归分析的地化录井烃类损失校正方法,其特征在于:所述步骤3中岩屑到壁心校正时:若因变量与自变量散点图呈非线性关系,选取y=x/(a*x+b)、y=a*(1-exp(-b*x))、y=a*exp(-b/x)以及y=a+b*log(x)非线性函数模型,然后计算每个非线性函数模型的相关系数及残差平方和,选择相关系数最大及残差平方和最小的函数模型,最后得到烃损恢复非线性表达式,当使用双曲线y=x/(a*x+b)模型进行非线性回归时,非线性回归表达式为:
S2(壁心)=S2(岩屑)/(a*S2(岩屑)+b)
式中,S2(壁心)为壁心热解烃含量(单位mg/g);S2(岩屑)为岩屑热解烃含量(单位mg/g);a、b为利用样本集求解的未知系数。
4.根据权利要求1所述的一种基于回归分析的地化录井烃类损失校正方法,其特征在于:所述步骤4中壁心到岩心校正时:若因变量与自变量散点图呈线性关系,使用多元线性回归分析,多元线性回归表达式为:
S2(岩心)=a′*S1(壁心)+b′*S2(壁心)+c′*pg(壁心)+d′*opi(壁心)+e′*tpi(壁心)+f′
式中:S2(岩心)为岩心热解烃含量(单位mg/g);S1(壁心)为壁心可溶烃含量(单位mg/g);S2(壁心)为壁心热解烃含量(单位mg/g);pg(壁心)为壁心产油潜量(单位mg/g);opi(壁心)为壁心油产率指数;tpi(壁心)为壁心油气产率指数;a′、b′、c′、d′、e′、f′为利用样本集求解的未知系数。
5.根据权利要求1所述的一种基于回归分析的地化录井烃类损失校正方法,其特征在于:所述步骤4中壁心到岩心校正时:若因变量与自变量散点图呈非线性关系,选取y=x/(a′*x+b′)、y=a′*(1-exp(-b′*x))、y=a′*exp(-b′/x)以及y=a′+b′*log(x)非线性函数模型,然后计算每个非线性函数模型的相关系数及残差平方和,选择相关系数最大及残差平方和最小的函数模型,最后得到烃损恢复非线性表达式,当使用双曲线y=x/(a′*x+b′)模型进行非线性回归时,非线性回归表达式为:
S2(岩心)=S2(壁心)/(a′*S2(壁心)+b′)
式中,S2(岩心)为岩心热解烃含量(单位mg/g);S2(壁心)为壁心热解烃含量(单位mg/g);a′、b′为利用样本集求解的未知系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811320271.XA CN109710887B (zh) | 2018-11-07 | 2018-11-07 | 一种基于回归分析的地化录井烃类损失校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811320271.XA CN109710887B (zh) | 2018-11-07 | 2018-11-07 | 一种基于回归分析的地化录井烃类损失校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109710887A CN109710887A (zh) | 2019-05-03 |
CN109710887B true CN109710887B (zh) | 2023-03-14 |
Family
ID=66254194
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811320271.XA Active CN109710887B (zh) | 2018-11-07 | 2018-11-07 | 一种基于回归分析的地化录井烃类损失校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109710887B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111122381B (zh) * | 2019-12-31 | 2022-07-15 | 中法渤海地质服务有限公司 | 一种基于回归分析的地化录井原油密度预测方法 |
CN114112651A (zh) * | 2020-08-27 | 2022-03-01 | 中国石油化工股份有限公司 | 一种用于人造岩心的岩石动静态力学参数转换方法及系统 |
CN112765765B (zh) * | 2020-12-16 | 2022-03-01 | 成都理工大学 | 一种井场油基钻井液岩屑热解结果的校正方法 |
CN115394052B (zh) * | 2022-08-30 | 2023-06-27 | 重庆地质矿产研究院 | 基于机器学习获得地质灾害预警关键参数预测值的方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103543470A (zh) * | 2013-10-18 | 2014-01-29 | 中国石油大学(华东) | 一种热解岩石中的游离烃/残留烃s1的轻、重烃校正方法 |
CN105403585A (zh) * | 2015-10-28 | 2016-03-16 | 中国石油天然气股份有限公司 | 烃源岩内分散液态烃丰度下限的确定方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ZA200509840B (en) * | 2003-05-07 | 2007-03-28 | Saudi Arabian Oil Co | Compositional modeling and pyrolysis data analysis methods |
-
2018
- 2018-11-07 CN CN201811320271.XA patent/CN109710887B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103543470A (zh) * | 2013-10-18 | 2014-01-29 | 中国石油大学(华东) | 一种热解岩石中的游离烃/残留烃s1的轻、重烃校正方法 |
CN105403585A (zh) * | 2015-10-28 | 2016-03-16 | 中国石油天然气股份有限公司 | 烃源岩内分散液态烃丰度下限的确定方法 |
Non-Patent Citations (1)
Title |
---|
地化录井技术在长庆区块的应用;李会;《大庆石油地质与开发》;20050831;第24卷;第103-106页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109710887A (zh) | 2019-05-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109710887B (zh) | 一种基于回归分析的地化录井烃类损失校正方法 | |
CN103543470B (zh) | 一种热解岩石中的游离烃/残留烃s1的轻、重烃校正方法 | |
CN107966545B (zh) | 一种任意成熟度烃源岩的有机碳恢复方法及其装置 | |
CN110735635A (zh) | 一种页岩含气量测试中损失气含量的确定方法 | |
CN113553780B (zh) | 一种基于机器学习的地层孔隙压力预测方法 | |
CN103530621A (zh) | 一种基于bp神经网络的煤岩图像识别方法 | |
CN113820750A (zh) | 基于弹塑性力学定量预测泥岩构造裂缝的方法 | |
CN103628816B (zh) | 地震属性指导水平井导向的方法及其装置 | |
CN106202748B (zh) | 一种基于渗透率和井距的煤层气采收率预测方法 | |
CA2914272C (en) | Hybrid approach to assisted history matching in large reservoirs | |
CN109063296B (zh) | 页岩含气量随钻计算方法 | |
CN107909202B (zh) | 一种基于时间序列的油井产液量集成预测方法 | |
CN113671564B (zh) | 基于narx动态神经网络的微地震有效事件自动拾取方法 | |
CN105277981A (zh) | 基于波场延拓补偿的非一致性时移地震面元匹配方法 | |
CN113393051A (zh) | 基于深度迁移学习的配电网投资决策方法 | |
CN113944459B (zh) | 一种页岩原位转化可采油气资源量的预测方法及装置 | |
CN110673211B (zh) | 一种基于测井与地震数据的品质因子建模方法 | |
CN114575827B (zh) | 一种随钻测量数据智能处理系统及方法 | |
CN111650647B (zh) | 一种基于地震数据约束的声波测井曲线重构方法 | |
Feng et al. | A new method for evaluating the unstable deliverability of gas wells in gas formation testing phase | |
CN112162317B (zh) | 一种基于地震波形横向差异预测薄储层的方法 | |
CN113917238B (zh) | 基于电磁辐射信号时域频域特征的煤岩动力灾害预警方法 | |
CN112859164B (zh) | 一种致密砂岩气页岩压裂效果预判断方法及系统 | |
CN113946928B (zh) | 一种有效烃源岩物性参数的预测方法及装置 | |
CN109143365B (zh) | 山前带地震静校正方法及计算机可读存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |