基于多条测井曲线的沉积旋回划分方法技术

本发明专利技术公开了一种基于多条测井曲线的沉积旋回划分方法，具体划分步骤为：获取测井曲线，对测井曲线进行预处理后得到曲线集合A，从曲线集合A中提取最具有特征性曲线，对最具有特征性曲线进行奇异谱分解得到重建曲线集合B，利用重建曲线的信号对沉积旋回进行划分。本发明专利技术的沉积旋回划分方法，克服了现有技术中单一测井曲线划分沉积旋回效率低且干扰因素多的问题。的问题。的问题。

【技术实现步骤摘要】
基于多条测井曲线的沉积旋回划分方法


[0001]本专利技术属于物理勘探
，涉及一种基于多条测井曲线的沉积旋回划分方法。

技术介绍

[0002]在油气勘探开发中，沉积地层是油气资源的富集区域和研究重点,受到沉积构造运动周期性的影响，沉积地层的岩石物理、化学性质呈现规律性的交替变化，表现为地层沉积具有旋回性。研究地层沉积旋回，有助于实现精确的地层对比，构建地层格架，分析地层沉积特征，进而发现潜在的有利储层,测井曲线作为记录地层岩性、物性变化的基础资料，因此，地层沉积旋回的划分尤为重要。
[0003]目前，沉积旋回曲线，常用的方法是对沉积旋回性敏感的测井曲线(如自然伽马、自然电位、电阻率等)进行划分，或直接划分沉积旋回，或利用傅里叶变化、小波变化、经验模态分解等变换方法，转换至时—频域提取频率域特征划分沉积旋回。这些沉积旋回划分方法往往都是首先针对单一曲线，然后进行综合分析，最终确定各级沉积旋回，具体实施时，往往会存在工作量巨大，且受个人经验等因素影响，干扰因素多的问题。

技术实现思路

[0004]本专利技术的目的是提供一种基于多条测井曲线的沉积旋回划分方法，解决了现有技术中基于单一曲线的划分结果，进行人为的综合判定的沉积旋回划分方法，工作量巨大，且受个人经验等因素影响，干扰因素多的问题。
[0005]本专利技术所采用的技术方案是，一种基于多条测井曲线的沉积旋回划分方法，具体步骤为：获取测井曲线，对测井曲线进行预处理后得到曲线集合A，从曲线集合A中提取最具有特征性曲线，对最具有特征性曲线进行奇异谱分解得到重建曲线集合B，利用重建曲线的信号对沉积旋回进行划分。
[0006]本专利技术的特点还在于，
[0007]获取测井曲线的具体步骤为：收集工区范围内地质、测井资料，确定沉积旋回划分的目标层段，地层岩性等信息，挑选出对沉积旋回性敏感的M条测井曲线。
[0008]测井曲线预处理步骤为：分别对每条测井曲线进行异常值校正、归一化处理后，利用线性拟合的方法计算各曲线之间的相关系数，通过数学变换确保所有曲线之间保持正相关关系，得到测井曲线的基础曲线集合A＝{a1,a2,a3......,a
M
},M≥2。
[0009]重建曲线集合的具体步骤为：根据奇异谱各分量的贡献率，划分不同组合，进行奇异谱重建，获得含N条重建曲线的集合B＝{b1,b2,b3......,b
N
}。
[0010]沉积旋回划分的原则为：中—长周期旋回和短周期旋回。
[0011]本专利技术的有益效果是：本专利技术通过融合多条测井曲线特征，创建具有代表性的曲线，分析奇异谱分解中各分量贡献率，划分多个重建组合，进行奇异谱重建，从而突出有效信息，降低噪声影响，利用重建曲线的不同信号将沉积旋回划分为中—长周期沉积旋回和
短周期沉积旋回，降低了沉积旋回划分工作量，提高了划分的准确度。
附图说明
[0012]图1是本专利技术基于多条测井曲线的沉积旋回划分方法的流程图；
[0013]图2是本专利技术基于多条测井曲线的沉积旋回划分方法中根据地层岩性选择的测井曲线图和降维提取的新曲线图；
[0014]图3是本专利技术基于多条测井曲线的沉积旋回划分方法中奇异谱各分量贡献率图；
[0015]图4是本专利技术基于多条测井曲线的沉积旋回划分方法中基于奇异谱重建沉积旋回划分图。
具体实施方式
[0016]下面结合附图和具体实施方式对本专利技术进行详细说明。
[0017]本专利技术的一种基于多条测井曲线的沉积旋回划分方法，如图1所示具体步骤为：
[0018]步骤1、收集工区范围内地质、测井资料，确定沉积旋回划分的目标层段，地层岩性等信息，挑选出对沉积旋回性敏感的M条测井曲线；
[0019]步骤2、分别对每条测井曲线进行异常值校正、归一化处理后，利用线性拟合的方法计算各曲线之间的相关系数，通过数学变换确保所有曲线之间保持正相关关系，得到测井曲线的基础曲线集合A＝{a1,a2,a3......,a
M
},M≥2。
[0020]步骤3、采用降维算法从基础曲线集合中提取最具有特征性曲线T；
[0021]步骤4、对曲线T进行奇异谱分解，根据奇异谱各分量的贡献率，划分不同组合，进行奇异谱重建，获得含N条重建曲线的集合B＝{b1,b2,b3......,b
N
}；
[0022]步骤5、利用重建曲线集合，根据重建曲线的不同信号将沉积旋回划分为中—长周期旋回和短周期旋回。
[0023]下面结合具体实施例验证本专利技术方法的有效性：
[0024]实施例1
[0025]收集工区内地质和测井等资料，确定划分沉积旋回的目标层段和主要岩性，区内地层以砂泥岩为主，据此选择自然伽马(GR)、自然电位(SP)、声波时差(DT)和电阻率(RT)作为划分沉积旋回的基础曲线。
[0026]以Well
‑
1井为例，目标层段深度为710
‑
860m，提取4条基础曲线进行异常值校正、归一化处理，然后分别计算各曲线之间的相关系数R，其中R
GR
‑
SP
＝0.66，R
GR
‑
DT
＝0.67，R
GR
‑
RT
＝
‑
0.35，R
SP
‑
DT
＝0.38，R
SP
‑
RT
＝
‑
0.23，R
DT
‑
RT
＝
‑
0.24。由此可知，GR与SP、DT正相关，RT与其他曲线负相关。本次实例中设定RT
′
＝
‑
RT，使得RT
′
与其他曲线转变为正相关关系。图2为均呈现正相关的基础曲线集合A＝{GR,SP,DT,RT'}。
[0027]本实例中采用了主成分分析(PCA)，通过线性变换，将多维数据映射到低维子空间中，从而实现降维目标。采用PCA降维算法，从基础曲线集A中获得最大特征值，提取出一条最具特征性的新曲线T。
[0028]以790
‑
810m层段为例如图2中所示，新曲线T融合了多曲线特征，保留了具有代表性的主要信息。
[0029]对新曲线T进行奇异谱分解，分析窗长范围2
‑
5m，本实例中设定最小分析窗长为
2.5m，获得奇异谱分量总数为统计各分量方差{var1,var2,var3....var
60
}，定义各分量贡献率为其中var
i
为奇异谱分量方差，图3展示了奇异谱各分量贡献率图，其上标示出多个贡献率变化拐点。根据奇异谱定义和各分量贡献率图，将60个奇异谱分量划分为D＝1，D＝2
‑
3，D＝4
‑
13，D＝14
‑
23和D＝24
...

【技术保护点】

【技术特征摘要】
1.一种基于多条测井曲线的沉积旋回划分方法，其特征在于，具体划分步骤为：获取测井曲线，对测井曲线进行预处理后得到曲线集合A，从曲线集合A中提取最具有特征性曲线，对最具有特征性曲线进行奇异谱分解得到重建曲线集合B，利用重建曲线的信号完成对沉积旋回进行划分。2.根据权利要求1所述的一种基于多条测井曲线的沉积旋回划分方法，其特征在于，所述获取测井曲线的具体步骤为：收集工区范围内地质、测井资料，确定沉积旋回划分的目标层段，地层岩性等信息，挑选出对沉积旋回性敏感的M条测井曲线。3.根据权利要求1所述的一种基于多条测井曲线的沉积旋回划分方法，其特征在于，所述测井曲线预处理步骤为：分别对每条测井曲线进行异常值校正、...

【专利技术属性】
技术研发人员：杨兴利，柳朝阳，辛毅超，黄博，王振宇，寸少妮，
申请(专利权)人：延长油田股份有限公司，
类型：发明
国别省市：