[0035] 一般情况下，对于井径严重扩大段，若要用本方法对曲线进行重构W矫正严重失 真的测井曲线，可选取扩径段上下井况较好的围岩地层作为标准层。而对于全井段的某系 列测井曲线失真或缺失的情况，则可W W邻井相同层位作为标准层。
[0036] 本发明实施例W距M14井最近的M7井的馆陶组地层为标准层，对M7井馆陶组进 行声波时差测井，M7井馆陶组测井曲线见图2。从图2中可W看出，M7井馆陶组各测井曲 线质量较好，曲线间的相互对应关系也较好，固可作为标准层。但需要注意的是，在M7井 1168m?1178m层段，CAL(井径）指示有明显扩径，固需要在标准层中剔除该段扩径率大于 10 %的层位。
[0037] 103 ;建立多元回归模型，并获取多元回归方程；
[003引本方法W式1所示的多元线性回归模型为基础，将式2?4所示的指数、对数、幕 函数多元非线性回归模型转化为如式5?7所示的线性回归模型，W便于在EXCEL软件中 进行回归计算。
[0039] 线性；y = a+b〇x〇+biXi+…+bnX。 式 1
[0040] 指数；J； = 06心川叶..+6"。 式 2
[0041] 对数；y = a+b〇Ln (x〇) +biLn (Xi) + …+bnLn (X。） 式 3
[00创幕函数；y =皆…疗 式4
[00创式中，y为待重构的测井曲线（即因变量），Xi(i = 0, 1. . .，n)为参与曲线重构的 母曲线（即自变量），a为回归常数，bi(i = 0、1，…，n)称为因变量y对自变量Xi的偏回归 系数。
[0044] 将指数、对数、幕函数非线性回归模型转化为线性模型分别为：
[0045] 指数；Y = A+b〇x〇+biXi+…+bnX。 式 5
[0046] 式中，Y = Ln (y)，A = Ln (a)
[0047] 对数；y = a+b〇X〇+biXi+... +bnXn 式 6
[0048] 式中，Xi= Ln(x i) (i = 0, 1,…，n)
[0049] 幕函数；Y = A+b〇X〇+biXi+... +bnXn 式 7
[0050] 式中，Y = Ln(y)，Xi= Ln(x i) (i = 0, 1,…，n)，A = Ln(a)。
[0化1] 根据上述多元回归模型，即可得到用于最终曲线重构的多元回归方程，分别为：
[005引 线性；y = a+b0x0+biXi+…+bnX。 式 8
[005引 指数；y=e片'(。柏脚+柏+'..+V。] 式9
[0054]对数；y = a+b〇Ln(x〇)+biLn(Xi)+...+bnLn(Xn) 式 10
[005引 幕函数：y=e护'("柄虹(、。)叫虹柄片哨虹w] 式11
[0化6] 对于线性函数来说，多元回归模型和多元回归方程一致，即公式1和公式8中的表 达式相同。
[0057] 104 ;最优母曲线及最优模型的选取；
[005引对用于重构测井曲线的母曲线的优选方法，首先是选取标准层目标曲线（即待重 构的曲线在标准层中所对应的相同曲线）与其余曲线（除目标曲线W外的所有常规曲线） 按照线性、指数、对数W及幕函数关系进行一一拟合，分析其相关系数。本发明实施例在计 算目标曲线与其余曲线之间的相关性时，主要通过EXC化软件一一绘制目标曲线与其余曲 线的散点图，再绘制出各自的拟合趋势线，并根据EXCEL软件自动计算出相关系数（R2)。本 发明实施例W重构AC曲线（声波时差）为例，分析AC与其余曲线的相关系数，如表1所示。
[0059] 表1中可W看出，AC与GR(伽马）、Rd(深电阻率）、C化（中子）、DEN(密度）曲线 之间的相关系数较高，固可优选该4类曲线作为参与多元回归计算的母曲线，而Depth(深 度）与SP(自然电位）曲线的相关系数较低，固不选取。通常情况下，相关系数大于0.25 的母曲线即满足选取要求，在不同工区，该标准可根据人工经验稍作调整，本发明实施例对 此不做限制。
[0060] 表1地层各曲线相关系数分析
[0061]
【主权项】
1. 一种重构测井曲线的方法，其特征在于，所述方法包括w下步骤： 对与目标井相邻的井进行常规测井资料的收集； 选取目标曲线的标准层； 建立线性、指数、对数、幕函数的多元回归模型，并获取相应的多元回归方程； 计算标准层目标曲线与其余各曲线之间的相关性，选取出用于重构目标曲线的母曲 线，在选出的母曲线的基础上，根据最大相关性原则选出最优多元回归模型，并推导出相应 的多元回归方程； 将最优多元回归模型结合优选出来的母曲线，代入EXC化软件线性多元回归模块，计 算回归模型中的各回归常数与系数； 在目标曲线缺失或失真的井段，将回归常数、系数W及母曲线值代入多元回归方程即 完成对目标曲线的重构。
2. 根据权利要求1所述的一种重构测井曲线的方法，其特征在于，所述目标曲线的标 准层具体为： 当目标曲线仅为部分缺失或失真时，选取目标曲线所在层段的上下围岩段作为标准 层；否则，选取邻井相同层位的地层作为标准层。
3. 根据权利要求1所述的一种重构测井曲线的方法，其特征在于，所述与目标井相邻 的井具体为；与目标井处于同一构造、同一工区的邻井。
4. 根据权利要求3所述的一种重构测井曲线的方法，其特征在于，所述邻井为与目标 井距离最近的井。
5. 根据权利要求1所述的一种重构测井曲线的方法，其特征在于，所述计算标准层目 标曲线与其余各曲线之间的相关性具体为： 通过EXC化软件一一绘制目标曲线与其余曲线的散点图，再绘制出各自的拟合趋势 线，并根据EXC化软件自动计算出相关系数。
6. 根据权利要求1所述的一种重构测井曲线的方法，其特征在于，所述方法还包括；通 过复相关系数和系数显著水平对所述重构目标曲线多元回归方程获得的重构目标曲线进 行评价。
【专利摘要】本发明公开了一种重构测井曲线的方法，涉及石油天然气勘探开发技术领域，所述方法包括以下步骤：对与目标井相邻的井进行常规测井资料的收集；选取目标曲线的标准层；建立线性、指数、对数、幂函数的多元回归模型，并获取相应的多元回归方程；计算标准层目标曲线与其余各曲线之间的相关性，选取出用于重构目标曲线的母曲线，在选出的母曲线的基础上，根据最大相关性原则选出最优多元回归模型，并推导出相应的多元回归方程；在目标曲线缺失或失真的井段，将回归常数、系数以及母曲线值代入多元回归方程即完成对目标曲线的重构。本发明实现了测井曲线的快速、精准重构，操作简便，易于推广，适用于现场工程、地质人员对地层信息进行快速评价。
【IPC分类】E21B47-00, E21B49-00, G06F17-50
【公开号】CN104533400
【申请号】CN201410645862
【发明人】邓强, 谭伟雄, 王俊瑞, 张志虎, 秦磊
【申请人】中国海洋石油总公司, 中海油能源发展股份有限公司
【公开日】2015年4月22日
【申请日】2014年11月12日