%在智能钻井研究工作中,使用数据的噪声处理,具有非常重要的现实意义,该方法使用了多种滤波方法。
%处理的数据主要是测井数据
% nepu liw 2017年9月3日
%-------------------------------------------------------------------
clc;
clear;
Logdata=importdata('z_da01.txt','\t');
[col,row]=size(Logdata);
% disp(Logdata);
%-------------------------------------------------------------------
Deep=Logdata(:,1);
Den=Logdata(:,2);
Den1=smooth(Den,30,'loess');
Den2=smooth(Den,30,'lowess');
[x1,y1]=size(Deep);
[x2,y2]=size(Den1);
disp([x1,y1,x1,y1]);
subplot(2,1,1);
% plot(Deep,Den,'r:','linewidth',1.5);xlabel('H (m)');ylabel('GR API');
plot(Deep,Den,'r','linewidth',1.5);xlabel('H (m)');ylabel('Den (g/cm3)');
hold on;
% plot(Deep,Den2,'k','linewidth',1.5);xlabel('H (m)');ylabel('GR API');
plot(Deep,Den1,'k','linewidth',2);
% legend('原始的密度测井数据')
% legend('原始的伽玛测井数据','噪声处理后的伽玛测井数据')
legend('原始的密度测井数据','噪声处理后的密度测井数据')
subplot(2,1,2);
plot(Deep,Den1,'k:','linewidth',2);
hold on;
plot(Deep,Den2,'k','linewidth',2);xlabel('H (m)');ylabel('Den (g/cm3)');
legend('loess','lowess')
set(0,'defaultfigurecolor','w') %设置MATLAB中figure的背景为白色
save z_data02.txt Den2 -ASCII -TABS
ddd=polyfit(Deep,Den2,2);
disp(ddd);
计算图例:

作者:LiW

本文探讨了测井等数据中的噪声处理技术,包括计算方法和应用,同时详细解析了数据图例的创建与解读,为地质分析提供关键支持。

9241

被折叠的 条评论
为什么被折叠?



