《测井曲线自动分层问题的解决.pptVIP

测井曲线自动分层问题 贺 衎 太原理工大学 基本内容 问题的理解和分析 模型的假设和建立 计算流程 数据处理 综合曲线的合成（包括计算权重因子） 模型求解的过程 利用趋势分析法划分井层 斜率的计算 层界面归并 模型的求解和检验 模型评价 问题的理解和分析 人工分层与自动分层的优缺点 问题主要是对数据进行处理 数据存在异常的原因： 1. 仪器设备干扰； 2. 地质结构复杂； 3：刻度和量纲所产生的影响.（见附图） 问题主要任务 根据1号井的数据建立数学模型，再根据由数学模型计算得出的2号—7号井自动分层与人工分层的结果比较，再对数学模型进行适当的调整，将数学模型进行优化处理，然后根据8号—13号井的各类曲线数据，运用所建立的数学模型将其自动分层后再分析结果. 模型的假设 假设测井数据真实可靠，不考虑人为因素引起的误差； 假设以1号井位标准，它包含所有的地质分层情况 假设不考虑地层变化引起的分层波动； 假设2号—7号井所给的分层和命名是准确的； 假设题目测量附件中所给的数据的仪器的精度满足要求； 假设根据某一个条或某几个测井曲线也可较为精确的进行自动分层 模型的建立计算流程 本文通过数据处理、分层处理两个阶段实现自动分层。 数据处理阶段的任务 剔除一些指标项，采用中值滤波的方法消除因仪器设备产生的测量数据异常；采用物理属性分析和主成分分析找出最能反映该井的测井曲线；首先统一量纲，采用加权因子分析合成综合曲线. 分层处理 自动分层的核心部分，给出了聚类分析、极值方差分析方法、趋势分析等等常用分层方法. 层界面归并和计算测井值. 数据处理 中值滤波 原始的数据矩阵（矩阵的行是条测井曲线; 矩阵的列是某一条测井曲线上的个观测点，即统一量纲） 归一化处理（由于不同的测井数据有不同的量纲, 为了统一量纲, 要进行归一化的数据变换） 综合曲线的合成（包括计算权重因子） 模型求解的过程利用趋势分析法划分井层 趋势分析是一种既直观又有效的分层方法，特别是对于层中心对称的曲线有很好的分层效果，通过计算测井曲线上各点的切线斜率，可以判断层界面的大体位置。 斜率的计算 层界面归并 为了消除因仪器设备造成的测量数据扰动，采用中值滤波对数据进行处理，但中值滤波也无法保证滤掉所有由仪器测量误差所产生的尖峰，这将在分层计算中分出一些不必要的小薄层，这时需要利用归并薄层界面来消除。 模型的求解和检验数据的处理阶段 对样本进行分层处理 对一号井样本进行聚类分析：根据权重分析原理得出四个参考指标SP1,SP,GR,AC加权因子0.1,0.3,0.5,0.1，再计算出一系列综合指标。 自动分层后与1号井手工分层进行对比校正 根据对1号井建立的模型，对2-7号井进行自动分层并与手工分层结果进行比较验证 自动分层结果 模型评价 附录 10.1去除无效点 1)[m,n]=size(data) for i=1:1:m for j=1:1:n if(data(i,j)==-9999||data(i,j)==+9999) data(i,j)=0; end end end 10.2中值滤波和极值归一化的matlab代码 1)function d=midfilt(x,n) p=size(x); x1=x; x2=x1; for j=1:p(2)-n+1 for i=1:p(1)-n+1 c=x1(i:i+(n-1),j); e=c(1,:); for u=2:n e=[e,c(u,:)]; end mm=median(e); x2(i+(n-1)/2,j)=mm; end end d=double(x2); 2)b=a; FlattenedData = b(:); MappedFlattened = mapminmax(FlattenedData, 0, 1); MappedData = reshape(MappedFlattened, size(b)); c=MappedData; * Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. Evaluation only. C