110 likes | 338 Views
MatLab. -et krasjkurs. www.stud.ntnu.no/~gjerding/data. >> load_well; >> well well = depth: [18163x1 double] Vp: [18163x1 double] Vs: [18163x1 double] rho: [18163x1 double] name: {1x9 cell} top: [1x9 double] bottom: [1x9 double].
E N D
MatLab -et krasjkurs
>> load_well; >> well well = depth: [18163x1 double] Vp: [18163x1 double] Vs: [18163x1 double] rho: [18163x1 double] name: {1x9 cell} top: [1x9 double] bottom: [1x9 double]
Plott Vp mot rho for ’Ty’ formasjonen og trekk en regresjonslinje Top Ty Bunn Ty
rho Vp Gardner et al. (1974): a? Brønn:
Top ty id= Bunn ty close all; clear load_well; frm = 4; id = find( well.depth >= well.top(frm) & well.depth <= well.bottom(frm) ) ; % ...the regression coefficients a = sum( well.Vp(id).^.25.* well.rho(id) ) ... /sum(( well.Vp(id) .^.25).^2); x = [ min(well.Vp(id)) : max(well.Vp(id)) ]; subplot(2,2,1) hold on plot( well.Vp(id) , well.rho(id) ,'r.' ) legend('well',4) plot(x,a*x.^.25,'b') xlabel('Vp [m/s]'); ylabel('rho [g/cc]'); title(well.name(frm)) axis tight grid str = num2str(a); text(2380,2.20,strcat('Gardner: rho = ',str,'Vp^0^.^2^5'))
Plotte seismikk data = load('someData'); Seiswiggle(data) mesh(data) view(60,50) View(0,90)
Andre ting • SEGY => MatLab • www.stud.ntnu.no/~gjerding • MatLab – Fortran / C++ • MatLab – unix • Perl, awk, sed etc. • unix