MatLab -et krasjkurs. gjerding/data.

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]

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

Vp rho

1 1

2 2

. .

. .

n n

Bunn Ty

Top Ty

25.31.0 pV

25.paV

Gardner et al. (1974):

Brønn: a?

Vp

rho

min225. paV

i

ii

ip

ip

V

Va

41

,

41

,

close all; clearload_well;

frm = 4;id = find( well.depth >= well.top(frm) & well.depth <= well.bottom(frm) ) ;

% ...the regression coefficientsa = 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 onplot( 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 tightgrid

str = num2str(a); text(2380,2.20,strcat('Gardner: rho = ',str,'Vp^0^.^2^5'))

# depth Vp Rho

1 . . .

2 . . .

. .

.

. .

n n n n

id=

Bunn ty

Top ty

i

ii

ip

ip

V

Va

41

,

41

,

Plotte seismikkdata = load('someData');Seiswiggle(data)

mesh(data)view(60,50)View(0,90)

Fourier Transform

Andre ting

SEGY => MatLab www.stud.ntnu.no/~gjerding

MatLab – Fortran / C++ MatLab – unix

Perl, awk, sed etc. unix