Code covered by the BSD License  

Highlights from
SeisLab 3.01

image thumbnail

SeisLab 3.01

by

 

24 Jul 2007 (Updated )

Release 10.0301 of a collection of functions for analysis/display of seismic and well-log data.

Log_examples1.m
% Log_examples1
%       Examples of the usage of log-related functions 

keep WF
presets
global S4M
S4M.interactive=false;	% Run without user intervention (turn-off some pop-up 
                        % windows requesting user input)

%     Create a synthetic log with a sonic made up from Gaussian noise;
depth=(5000:0.5:10000)';
nsamp=length(depth);
dt=150+25*randn(nsamp,1);
wlog=l_convert([depth,dt],{'depth','ft','Depth';'DTp','us/ft','P-sonic'});
wlog.name='Synthetic log';

%     Add Kelley bushing elevation as a parameter
wlog=ds_add_parameter(wlog,84,{'KBE','ft','Kelly bushing elevation'});

%     Add a trend to the sonic log
wlog=l_curve_math(wlog,'replace','dtp=dtp*(5000/depth)^0.333');

%     Extract sonic curve from log
dtp=l_gc(wlog,'dtp');

%     Use "dtp" to add a curve with P-velocity (curve mnemonic 'Vp')
wlog=l_curve(wlog,'add','Vp',1.0e6./dtp,'ft/s','P-velocity');

%     Use an alternative method to add a density curve computed 
%     via Gardner's formula
wlog=l_curve_math(wlog,'add_ne','rho=0.23*vp^0.25','g/cm3','Density');

%     Plot curves
l_plot(wlog)

%     Change the units of measurement from feet to meter
wlog=ds_unit_conversion(wlog,{'ft','m'});

%     Change the units of measurement from feet/second to meters/second
%     and from microseconds/foot to microseconds/meter.
wlog=ds_unit_conversion(wlog,{'ft/s','m/s'},{'us/ft','us/m'});

%     Display summary information about the log curves
disp(' ')
l_curve(wlog)

%     Compute density trend (requires the Optimization Toolbox)
if exist('fmincon','file') > 0
   disp(' ')
   disp('Regression:')
   [wlog,aux]=l_regression(wlog,'rho=x1*depth^x2',{'lbounds',0,0}, ...
              {'ubounds',10,1},{'mnem','rho_trend'});

%     Compare density with density trend
   l_plot1(wlog,{'curves','rho','rho_trend'})
   a=num2str(round(aux.x1*1000)/1000);
   b=num2str(round(aux.x2*10000)/10000);
   mytitle(['Compare density with computed density trend ',a,'*depth^{',b,'}'])
end

Contact us