Matlab l.jpg
This presentation is the property of its rightful owner.
Sponsored Links
1 / 11

MatLab PowerPoint PPT Presentation


  • 121 Views
  • Uploaded on
  • Presentation posted in: General

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].

Download Presentation

MatLab

An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Presentation Transcript


Matlab l.jpg

MatLab

-et krasjkurs


Www stud ntnu no gjerding data l.jpg

www.stud.ntnu.no/~gjerding/data


Slide5 l.jpg

>> 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]


Slide6 l.jpg

  • Plott Vp mot rho for ’Ty’ formasjonen og trekk en regresjonslinje

Top Ty

Bunn Ty


Slide7 l.jpg

rho

Vp

Gardner et al. (1974):

a?

Brønn:


Slide8 l.jpg

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 l.jpg

Plotte seismikk

data = load('someData');

Seiswiggle(data)

mesh(data)

view(60,50)

View(0,90)


Fourier transform l.jpg

Fourier Transform


Andre ting l.jpg

Andre ting

  • SEGY => MatLab

    • www.stud.ntnu.no/~gjerding

  • MatLab – Fortran / C++

  • MatLab – unix

    • Perl, awk, sed etc.

    • unix


  • Login