【技术实现步骤摘要】
本专利技术属于油气田开发，特别涉及一种过套管电阻率测井数据预处理方法及系统。
技术介绍
1、油田开发中后期，剩余油饱和度监测对油田增储上产和稳油控水具有重要意义。过套管电阻率测井开辟了在适当条件下有效和准确地计算套后饱和度的新前景，其解释技术简单、方法成熟。由于几乎所有裸眼井都使用电阻率进行原始地层含水饱和度的计算，因此只需将原始电阻率测井曲线与过套管电阻率测井曲线进行比较，就可以非常直接地估计当前的含油饱和度水平，在国内外油田得到了广泛应用。过套管电阻率测井仪器通过检测从套管体泄漏到地层中微弱的电流信号实现对地层电阻率的测量，为测量到可靠的资料，仪器在设计及测量方式上与一般测井仪器不同。因此，原始的过套管电阻率测井资料也与常见的测井资料有明显不同。与常规测井资料相比，过套管电阻率测井资料具有以下特点：
2、(1)测量值离散。过套管电阻率测井仪器均采用点测方式进行，仪器在套管井中的每一个测量点静止测量，因为仪器在行进中引入的噪声信号将远大于要检测的地层信号。点测的特点，决定了测量值在深度上是离散的。
3、(2)一个测量点有多个测量值，还可能含有异常的测量值。多次测量方式在能重复的观测系统中被广泛使用，目的是减小误差。过套管电阻率测井仪器在目的层段应要求有3次以上的测量，但重复测量的次数受到许多方面的制约，特别是仪器保护方面。过套管电阻率测井仪器在井下高温、高压环境下作业，仪器一次下井测量时间不宜过长，以免损坏仪器或缩短仪器的使用寿命。因此，重复测量的次数一般不超过10次。这样一来，同一深度测量的地层电阻率
4、(3)测量点之间的深度间隔不等。为了节省测量时间，过套管电阻率测井的深度间隔一般为0.5m，在薄目的层处可减小到0.3m，在非目的层段测量点之间的深度间隔则相对较大，如1m。由此，造成了一口井的测量数据整体上深度间隔不等。而常规的测井资料数字化后深度间隔一般为0.125m。
5、鉴于以上特点，需要对过套管电阻率测量数据进行预处理以形成单值、等深度间距的数据，便于后续处理，从而与其他测井资料配套应用。
6、目前，过套管电阻率测量值筛选常用的方法主要是格拉布斯准则筛选法、插值方法主要有拉格朗日插值、三次样条插值和线性插值方法。格拉布斯准则一次检测一个异常值。从数据集中删除该异常值，并且迭代测试直到没有检测到异常值。但是，多次迭代会改变检测概率，并且测试不应该用于六个或更少的样本大小(n>6)，因为它经常将大多数点标记为异常值。
7、目前，常用过套管电阻率的插值方法有牛顿、拉格朗日、埃尔米特等提出的全局多项式插值法及三次样条插值法，但这些方法容易出现形态摇摆现象，计算量较大。而且上述插值方法形成的线性插值曲线为折线，不具备光滑性。
8、为了解决上述问题，需要设计一种过套管电阻率测井数据预处理方法及系统，剔除异常测量值并插值形成等深度间隔、加密的数据。
技术实现思路
1、针对上述问题，本专利技术提供一种过套管电阻率测井数据预处理方法，所述方法包括：
2、读取原始数据，提取油井的深度与过套管电阻率；
3、利用箱型图对每一深度点上过套管电阻率进行处理，得到对应深度点上过套管电阻率的合理数据的上限值与下限值；
4、依据所述上限值和下限值剔除过套管电阻率的异常值，并对同一深度点上过套管电阻率的正常值进行处理生成过套管电阻率的单一数据；
5、利用akima插值法将过套管电阻率的单一数据插值成等深度间隔的数据。
6、优选地，利用箱型图对每一深度点上过套管电阻率进行处理包括以下步骤：
7、将过套管电阻率从大到小进行排序，得到过套管电阻率序列；
8、从所述过套管电阻率序列中提取出中位数、下四分位数q1和上四分位数q3；
9、计算出过套管电阻率的合理数据的上限值和下限值。
10、优选地，所述中位数为过套管电阻率序列中最中间的数或者最中间两个数的平均值；
11、所述下四分位数q1为过套管电阻率序列平均分成四份后取(1+n)/4的数；
12、所述上四分位数q3为过套管电阻率序列平均分成四份后取(1+n)*3/4的数；
13、其中，n表示某一深度下测量的过套管电阻率数据点个数。
14、优选地，所述上限值的计算方法包括：
15、计算q3+1.5iqr，得到上限对比值，其中四分位距iqr＝q3-q1；
16、比较所述上限对比值与过套管电阻率序列中极大值的大小，将较小的值作为上限值。
17、优选地，所述下限值的计算方法包括：
18、计算q1-1.5iqr，得到下限对比值，其中四分位距iqr＝q3-q1；
19、比较所述下限对比值与过套管电阻率序列中极小值的大小，将较大的值作为下限值。
20、优选地，所述异常值为低于下限值和高于上限值的过套管电阻率的值；
21、所述正常值为大于等于下限值且小于等于上限值的过滤套管电阻率的值。
22、优选地，所述过套管电阻率的单一数据的生成方法包括：
23、求出同一深度点上的过套管电阻率的所有正常值的算数平均值，将算数平均值作为该深度点上过套管电阻率的单一数据。
24、优选地，利用akima插值法将过套管电阻率的单一数据插值成等深度间隔的数据包括以下步骤：
25、在每两个单一数据之间建立一条由三次多项式拟合而成的曲线，所述曲线保证一阶导数连续；
26、将所有n个单一数据构建为不等间距样本点(xi,yi)(i＝0,1,…,n-1)，其中，xi表示第i个深度，yi表示与深度对应的合成过套管电阻率单一数据，x0<x1<…<xn-1；若在子区间[xk,xk+1](k＝0,1,2,…,n-2)上的2个端点处满足：
27、
28、则在任意2个相邻子区间可以确定唯一的三次多项式：
29、f(x)＝a0+a1(x-xk)+a2(x-xk)2+a3(x-xk)3，
30、由所述三次多项式计算子区间中的插值点x(x∈[xk,xk+1])处的函数近似值。
31、本专利技术还提供一种过套管电阻率测井数据预处理系统，所述系统包括：
32、信息提取模块，用于提取油井的深度与过套管电阻率；
33、信息处理模块，用于对每一深度点上过套管电阻率进行处理，得到对应深度点上过套管电阻率的合理数据的上限值与下限值；
34、数据生成模块，用于剔除过套管电阻率的异常值，并对同一深度点上过套管电阻率的正常值进行处理生成过套管电阻率的单一数据；
35、数据处理模块，用于将过套管电阻率插值成等深度间隔的数据。
36、优选地，所述信息处理模块用于对每一深度点上过套管电阻率进行处本文档来自技高网...
【技术保护点】
1.一种过套管电阻率测井数据预处理方法，其特征在于，所述方法包括：
2.根据权利要求1所述的一种过套管电阻率测井数据预处理方法，其特征在于，
3.根据权利要求2所述的一种过套管电阻率测井数据预处理方法，其特征在于，
4.根据权利要求3所述的一种过套管电阻率测井数据预处理方法，其特征在于，
5.根据权利要求3所述的一种过套管电阻率测井数据预处理方法，其特征在于，
6.根据权利要求1所述的一种过套管电阻率测井数据预处理方法，其特征在于，
7.根据权利要求1所述的一种过套管电阻率测井数据预处理方法，其特征在于，
8.根据权利要求1所述的一种过套管电阻率测井数据预处理方法，其特征在于，
9.一种过套管电阻率测井数据预处理系统，其特征在于，所述系统包括：
10.根据权利要求9所述的一种过套管电阻率测井数据预处理系统，其特征在于，
11.根据权利要求10所述的一种过套管电阻率测井数据预处理系统，其特征在于，
12.根据权利要求11所述的一种过套管电阻率测井数据预
13.根据权利要求11所述的一种过套管电阻率测井数据预处理系统，其特征在于，
14.根据权利要求9所述的一种过套管电阻率测井数据预处理系统，其特征在于，
15.根据权利要求9所述的一种过套管电阻率测井数据预处理系统，其特征在于，
...【技术特征摘要】
1.一种过套管电阻率测井数据预处理方法，其特征在于，所述方法包括：
2.根据权利要求1所述的一种过套管电阻率测井数据预处理方法，其特征在于，
3.根据权利要求2所述的一种过套管电阻率测井数据预处理方法，其特征在于，
4.根据权利要求3所述的一种过套管电阻率测井数据预处理方法，其特征在于，
5.根据权利要求3所述的一种过套管电阻率测井数据预处理方法，其特征在于，
6.根据权利要求1所述的一种过套管电阻率测井数据预处理方法，其特征在于，
7.根据权利要求1所述的一种过套管电阻率测井数据预处理方法，其特征在于，
8.根据权利要求1所述的一种过套管电阻率测井数...
【专利技术属性】
技术研发人员：黄奇志，杨原军，冯敏，王瑞峰，王敏，杨轩宇，廖长霖，徐庆岩，雷诚，
申请(专利权)人：中国石油天然气集团有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。