Scribd
Upload
30 views2 pages

Ray Tracing.m

The document describes modeling seismic refraction tomography travel time data. It defines the geometry of sources and receivers, velocity models for different blocks, and performs ray tracing on the true velocity model. Travel times are calculated based on ray paths and slopes in different blocks. A graph of travel time vs distance is plotted.

Uploaded by

Sekmawati Ar
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
30 views2 pages

Ray Tracing.m

The document describes modeling seismic refraction tomography travel time data. It defines the geometry of sources and receivers, velocity models for different blocks, and performs ray tracing on the true velocity model. Travel times are calculated based on ray paths and slopes in different blocks. A graph of travel time vs distance is plotted.

Uploaded by

Sekmawati Ar
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 2

% Pemodelan Tomografi Seismik Refraksi

% Desain data Sintetik Travel Time


clear all; clc;
nx = 100; % Jumlah blok dalam arah horizontal
ny = 2; % Jumlah blok dalam arah vertikal
bly = 20; % Dimensi dalam arah vertikal
x = [ 0 5 1000]; % Mendifnisikan skala sumbu
horizontal
y = [0 5 40]; % Mendifinisikan skala sumbu vertikal

% Membuat konfigurasi geometri sumber dan receiver


ys = [0];
ns = length(ys);
xs = 10;
xr = [300:60:1000];
nr = length(xr);
yr = zeros(1,nr); % receiver

% Mendefinisikan kecepatan pada setiap blok


% jumlah blok = nx*ny
v1 = 400;
v2 = 820;
V1 = v1*ones(1,nx);
V2 = v2*ones(1,nx);
V = [V1 V2];
velo = reshape (V,nx,ny);

figure (1)
a = subplot(211);
imagesc(x,-y,velo');
shading('flat');
colormap(cool);
title('Ray Tracing pada true Model');
xlabel ('Jarak (m)');
ylabel ('Kedelaman(m)');
axis xy;
axis ([min(x) max(x) -max(y) min(y)]);
hold on
colorbar('vert');
ve1 = reshape(V',1,ny*nx);
hold on
h = bly;
for m = 1:ns
plot(xs(m),-ys(m),'*k');
hold on
for k = 1:nr
plot(xr(k),-yr(k),'ok');
hold on
plot(xr(k),-yr(k),'ok');
hold on
end
end

for m = 1:ns
plot(xs(m),-ys(m),'*k');
hold on
for k = 1:nr
plot(xr(k),-yr(k),'ok');
hold on
plot(xr(k),-yr(k),'ok');
hold on
mk = (m-1)*nr+k;
tetaC = asind(v1/v2);
X(mk) = xr(k)-xs(m);
X1 (mk) = h*tand(tetaC);
X3 (mk) = h*tand(tetaC);
X2 (mk) = X(mk)-(X1(mk)+X3(mk));
P = h/cosd(tetaC);
T1 (mk) = P/v1;
T2 (mk) = (X(mk)-2*h*tand(tetaC))/v2;
T3 (mk) = T1(mk);
Tr (mk) = T1(mk) + T2(mk) + T3(mk);
xray = [xs(m) xs(m)+X1(mk) xr(k)-X3(mk)
xr(k)];
yray = [ys(m) h h yr(k)];
plot(xray,-yray,'r');
hold on;
drawnow;
end
end

% PLOT GRAFIK T-X


x = X; t = Tr;
b =subplot(212);
plot (x,t,'-or');
% slope =(t(NR)-T(1))/(X(NR)-X(1))
title('Grafik T-X');
xlabel ('X');
ylabel ('Tr');

You might also like