%%?------------------------------------------------------------?%%
%?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）;???