s_log2reflectivity.m

地震、测井方面matlab代码,解释的比较详细

function reflect=s_log2reflectivity(wlog,step,varargin)
% Create acoustic reflection coefficient series in time from well log in 
% depth or time. An impedance curve to use can be specified via the
% "redefine curve name" facility: {'aimp',impedance_curve}  where 
% "impedance_curve" is the mnemmonic of the impedance curve to use.
% Written by: E. R.: December 17, 2003
% Last updated: July 18, 2005: Impedance curve to use for reflectivity 
%           calculation can be chosen via the "redefine curve name" facility.
%
%           reflect=s_log2reflectivity(wlog,step,varargin)
% INPUT
% wlog      log structure with curves "aImp",  or "DT" and "rho", 
%           or "Vp" and "rho".
%           Mnemonics can be changed: For example {'dt','dtp'} implies the the 
%           sonic curve is represented by "dtp" (by default curve mnemonics 
%           are not case sensitive).
% step      seismic sample interval
% varargin  one or more cell arrays; the first element of each cell array is a keyword,
%           the other elements are parameters. Presently, keywords are:
%     'depth_time'  depth value and the corresponding time to establish time to top of log
%           the depth value must satisfy wlog.first <= depth_time{1} <= wlog.last
%           Default:{'depth_time',first,0} : "first" is the depth of the first
%           non-NaN  sample of the sonic or velocity log; thus first >= wlog.first
%
% OUTPUT
% reflect   seismic structure representing reflection coefficients

%	Set default parameters
param.depth_time=[];
 
[param,cm]=l_assign_input(param,varargin);

[dummy,ierr]=curve_index1(wlog,cm.aimp);
if ierr
   disp(' Acoustic impedance computed.')
   wlog1=l_seismic_acoustic(wlog,{'Vp',cm.vp},{'DTp',cm.dtp},{'rho',cm.rho}, ...
                        {'aImp',cm.aimp},{'twt',cm.twt});
else
   disp([' Reflectivity computed from curve "',cm.aimp,'".'])
   wlog1=wlog;
end
 
if ~iscurve(wlog1,cm.twt)
   if isempty(param.depth_time)
      wlog1=l_sonic2time(wlog1,{'output',cm.twt});
   else
      wlog1=l_sonic2time(wlog1,{'depth_time',param.depth_time{1},param.depth_time{2}}, ...
         {'output','twt'});
   end
elseif ~isempty(param.depth_time)
   wlog1=l_sonic2time(wlog1,{'action','replace'},{'output',cm.twt});
end

wlog1=l_switch_depth(wlog1,cm.twt);
wlog1=l_resample(wlog1,step);

impedance=s_convert(l_gc(wlog1,cm.aimp),wlog1.first,step);
reflect=s_reflcoeff(impedance);
reflect.name=['Reflectivity (',wlog.name,')'];