prototipe untuk program downhole seismik //prototipe for downhole seismic program

% u/ downhole seismic

% look for velocity in vertical position

load tz.dat

load vp2a.dat

load vs2.dat

tp=tz(:,2);

ts=tz(:,3);

z=tz(:,1);

figure(1)

plot(z,tp,’*b’),title(‘waktu tempuh gel P terhadap kedalaman’),xlabel(‘depth (m)’),ylabel(‘tp (ms)’)

figure(2)

plot(z,ts,’*g’),title(‘waktu tempuh gel s terhadap kedalaman’),xlabel(‘depth (m)’),ylabel(‘ts (ms)’)

% mencari kecepatan gel P dan gel S

n=length(z)-1;

for i=1:1:n

vp(i)=(z(i+1)-z(i))/(tp(i+1)-tp(i));

vs(i)=(z(i+1)-z(i))/(ts(i+1)-ts(i));

r(i,:)=z(i):z(i+1)

end

% V as depth function with range

k=[r(:,1)

r(:,2)]

g=sort(k)

% to be continue……………………

vp2=vp2a;

figure(3)

plot(g-23,vp2),xlabel(‘kedalaman (m)’),ylabel(‘vs’),title(‘Vs terhadap kedalaman’)

vpt=vp’;

vst=vs’;

vs2=vs2;

figure(4)

plot(g,vs2),xlabel(‘kedalaman (m)’),ylabel(‘vs’),title(‘Vs terhadap kedalaman’)


(program ini sebagian hanya dummy dan tidak terpakai)

Advertisement