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

上传人:小飞机 文档编号:6308110 上传时间:2023-10-15 格式:PPT 页数:26 大小:1.07MB
返回 下载 相关 举报
测井曲线自动分层问题的解决.ppt_第1页
第1页 / 共26页
测井曲线自动分层问题的解决.ppt_第2页
第2页 / 共26页
测井曲线自动分层问题的解决.ppt_第3页
第3页 / 共26页
测井曲线自动分层问题的解决.ppt_第4页
第4页 / 共26页
测井曲线自动分层问题的解决.ppt_第5页
第5页 / 共26页
点击查看更多>>
资源描述

《测井曲线自动分层问题的解决.ppt》由会员分享,可在线阅读,更多相关《测井曲线自动分层问题的解决.ppt(26页珍藏版)》请在三一办公上搜索。

1、测井曲线自动分层问题,贺 衎太原理工大学,基本内容,问题的理解和分析模型的假设和建立 计算流程 数据处理 综合曲线的合成(包括计算权重因子)模型求解的过程 利用趋势分析法划分井层 斜率的计算 层界面归并 模型的求解和检验模型评价,问题的理解和分析,人工分层与自动分层的优缺点问题主要是对数据进行处理数据存在异常的原因:1.仪器设备干扰;2.地质结构复杂;3:刻度和量纲所产生的影响.(见附图),问题主要任务 根据1号井的数据建立数学模型,再根据由数学模型计算得出的2号7号井自动分层与人工分层的结果比较,再对数学模型进行适当的调整,将数学模型进行优化处理,然后根据8号13号井的各类曲线数据,运用所建

2、立的数学模型将其自动分层后再分析结果.,模型的假设,假设测井数据真实可靠,不考虑人为因素引起的误差;假设以1号井位标准,它包含所有的地质分层情况假设不考虑地层变化引起的分层波动;假设2号7号井所给的分层和命名是准确的;假设题目测量附件中所给的数据的仪器的精度满足要求;假设根据某一个条或某几个测井曲线也可较为精确的进行自动分层,模型的建立计算流程,本文通过数据处理、分层处理两个阶段实现自动分层。数据处理阶段的任务 剔除一些指标项,采用中值滤波的方法消除因仪器设备产生的测量数据异常;采用物理属性分析和主成分分析找出最能反映该井的测井曲线;首先统一量纲,采用加权因子分析合成综合曲线.分层处理 自动分

3、层的核心部分,给出了聚类分析、极值方差分析方法、趋势分析等等常用分层方法.层界面归并和计算测井值.,开 始,结束,数据与处理,数据滤波器(中值滤波,归一化处理),主成分分析合成综合曲线,趋势分析,分层计算阶段,数据处理阶段,一次处理离散曲线点,确定上下界面,处理完毕,否,层界面归并,计算测井值,是,数据处理 中值滤波,原始的数据矩阵(矩阵的行是条测井曲线;矩阵的列是某一条测井曲线上的个观测点,即统一量纲),归一化处理(由于不同的测井数据有不同的量纲,为了统一量纲,要进行归一化的数据变换),式中:,代表第 条测井曲线;代表某一条测井曲线的第 个采样点;为第 条测井曲线的第 个采样点的数据;为第

4、条测井曲线数据的最大值;为 第 条测井曲线数据的最小值。归一化后数据的量纲一致,并且都在01 之间,综合曲线的合成(包括计算权重因子),模型求解的过程利用趋势分析法划分井层,趋势分析是一种既直观又有效的分层方法,特别是对于层中心对称的曲线有很好的分层效果,通过计算测井曲线上各点的切线斜率,可以判断层界面的大体位置。,斜率的计算,层界面归并,为了消除因仪器设备造成的测量数据扰动,采用中值滤波对数据进行处理,但中值滤波也无法保证滤掉所有由仪器测量误差所产生的尖峰,这将在分层计算中分出一些不必要的小薄层,这时需要利用归并薄层界面来消除。,模型的求解和检验数据的处理阶段,相对稳定又有适度波动的因素作为

5、主要考虑因素。:SP1,GR,AC,SP,对样本进行分层处理 对一号井样本进行聚类分析:根据权重分析原理得出四个参考指标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 endend,10.2中值滤波和极值归一化的m

6、atlab代码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;,

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 生活休闲 > 在线阅读


备案号:宁ICP备20000045号-2

经营许可证:宁B2-20210002

宁公网安备 64010402000987号