log_examples1.m

实现地震勘探中

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

clear all
presets
global S4M
S4M.interactive=0;	% Run without user intervention (turn-off some pop-up 
                        % windows requesting user intervention)

%	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