利用测井数据进行正演

%% ------------------------------------------------------------ %%
% DAQING Well logs plotting 1-3-30
% 
% - Written by Zhiqi Guo on 05 April 2012
% - during Marathon Shale Gas projects in Ediburgh Anisotropy Project （EAP）
% ----------------------------------------------------------------------- %


%% ------------------------------------------------------------ %%
clc; 
close all;
clear all;

%% ------------------------------------------------------------ %%                             
%--------------------- 1-Setting and data loading ------------------------%
 well          =  load（‘data\WEL1-3-30.txt‘）;    % well
%    well          =  importdata（‘input_well_jiaoye1_forBackus.txt‘）; 
  n              =    6;
  LineWidth      =    2;
  FontSize       =   12;
  AspectRatio_1  =  4/1;
  
target_well   =   1:1913;
% target_well   =   2242:3522;
 depth_range  =   length（target_well）;
  

interval = 1; 
sampling = 1:interval:depth_range;   % interval = 1 for 0.125m; 
                                     % interval = 8 for 1m; 
                               
   depth_well  =              well（target_well（sampling）1）;
      DT_well  =              well（target_well（sampling）2）;
      GR_well  =              well（target_well（sampling）3）;
     LLD_well  =              well（target_well（sampling）4）;
    NPHI_well  =              well（target_well（sampling）5）;
    RHOB_well  =              well（target_well（sampling）6）;
      SP_well  =              well（target_well（sampling）7）;
      
%    depth_well  =              well.data（target_well（sampling）1）;
%       vp_well  =              well.data（target_well（sampling）3）;
%       GR_well  =              well.data（target_well（sampling）2）;
%       vs_well  =              well.data（target_well（sampling）4）;
%      rho_well  =              well.data（target_well（sampling）5）;
      vp_well  =  1e3*0.3048./well（target_well（sampling）2）;    % from DT to Vp           
      vs_well  =  0.79*vp_well-0.79;                            % Vs from Vp  
     rho_well  = RHOB_well;
            
% -- from depth to time using Vp
time_well = zeros（1length（rho_well））;                           
time_well（1） = 2*0.125/（vp_well（1））;                            
for ii = 2:1:length（rho_well）
      time_well（ii） = time_well（ii-1） + 2*0.125/（vp_well（ii））;
end

% -- Ip and Rpp
Ip = vp_well.*rho_well;                                        
Ip（length（Ip）+1） = Ip（length（Ip））;
RefleCoeff = zeros（1length（rho_well））;
for ii = 1:1:length（rho_well）
    RefleCoeff（ii） = （Ip（ii+1） - Ip（ii）） / （Ip（ii+1） + Ip（ii））;
end

% -- Media filter  
  interval_medfilt = 1;
   GR_well_medfilt = medfilt1（GR_wellinterval_medfilt）;         
   vp_well_medfilt = medfilt1（vp_wellinterval_medfilt）;   
   vs_well_medfilt = medfilt1（vs_wellinterval_medfilt）;    
  rho_well_medfilt = medfilt1（rho_wellinterval_medfilt）;       
RefleCoeff_medfilt = medfilt1（RefleCoeffinterval_medfilt）;