5.

HASIL DAN PEMBAHASAN

5.1. Flowchart/Algoritma Program
5.1.1. Flowchart Percobaan 1

Start

Input x_off1, TT_RW1

h=50;
x_off=x_off1';
TT_RW=TT_RW1';
n =length(x_off);

for i=1:n

G(k,1) = 1;

G(k,2) = x_off(k);

d =TT_RW;
m =inv(G'*G)*G'*d;

V2=1/m(2);
V1=x_off(1)/(TT_RW(1)-(x_off(1)/V2));
teta=inv(sin(V1/V2));
Xc=2*h*tan(teta);
Xcr=2*h*(sqrt(V2+V1/V2-V1));
1
 i

m,V2,V1,teta,Xc,Xcr

G1=[ones(1,(length(TT_RW1)));x_off1]';
GT1=G1';
d1=TT_RW1;
Id=eye(2);
beta=0:0.1:1;

for i=1:length(beta)

ms(:,i)=inv(GT1*G+beta(i)^2.*Id)*GT1*d1';

V2=1/ms(2,11);
V1=x_off(1)/(TT_RW(1)-(x_off(1)/V2));
teta=inv(sin(V1/V2));
Xc=2*h*tan(teta);
Xcr=2*h*(sqrt(V2+V1/V2-V1));

ms,V2,V1,teta,Xc,Xcr

End
 5.2.1. KODE-KODE PROGRAM

Kode-kode program inversi linier dan inversi linier terkonstrain

clear all
clc

x_off1=[100 139 178 217 257 296 335 374 413 452 491 530 570 609 648 687....
726 765 804 844 883 922 961 1000];
TT_RW1=[0.0696 0.0808 0.0920 0.1031 0.1143 0.1255 0.1367 0.1479 0.1590...
0.1702 0.1814 0.1926 0.2038 0.2149 0.2261 0.2373 0.2485 0.2597...
0.2708 0.2820 0.2932 0.3044 0.3156 0.3267];
plot(TT_RW1,x_off1)
grid on
title('\fontsize{14} GRAFIK T Vs X')
xlabel('T')
ylabel('X')
h=50;
x_off=x_off1';
TT_RW=TT_RW1';
n =length(x_off);
for k = 1:n
G(k,1) = 1;
G(k,2) = x_off(k);
end
d=TT_RW;
m =inv(G'*G)*G'*d;
disp('hasil inversi linier adalah :')
disp(m)

% menentukan nilai kecepatan V2

V2=1/m(2);
fprintf('nilai kecepatan V2 inversi linier adalah = %1.3f\n',V2)

% menentukan kecepatan V1
V1=x_off(1)/(TT_RW(1)-(x_off(1)/V2));
fprintf('nilai kecepatan V1 inversi linier adalah = %1.3f\n',V1)

% menentukan sudut kritis

teta=inv(sin(V1/V2));
fprintf('nilai sudut kritis inversi linier adalah = %1.3f\n',teta)
% menentukan jarak kritis
Xc=2*h*tan(teta);
fprintf('nilai jarak kritis inversi linier adalah = %1.3f\n',Xc)

% menentukan jarak crossover

Xcr=2*h*(sqrt(V2+V1/V2-V1));
fprintf('nilai crossover inversi linier adalah = %1.3f\n\n',Xcr)

% INVERSI LINIER TERKONSTRAIN

G1=[ones(1,(length(TT_RW1)));x_off1]';
GT1=G1';
d1=TT_RW1;
Id=eye(2);
beta=0:0.1:1;
for i=1:length(beta)
ms(:,i)=inv(GT1*G+beta(i)^2.*Id)*GT1*d1';
end
disp('hasil damped least square terkonstrain dengan beta 0:0.1:1 adalah :')
disp(ms)

% menentukan kecepatan lapisan V2 melalui inversi linier terkonstrain

V2=1/ms(2,11);
fprintf('nilai kecepatan V2 inversi linier terkonstrain pada beta sama dengan 1
adalah = %1.3f\n',V2)

% menentukan kecepatan lapisan V1 melalui inversi linier terkonstrain

V1=x_off(1)/(TT_RW(1)-(x_off(1)/V2));
fprintf('nilai kecepatan V1 inversi linier terkonstrain pada beta sama dengan 1
adalah = %1.3f\n',V1)

% menentukan sudut kritis

teta=inv(sin(V1/V2));
fprintf('nilai sudut kritis inversi linier terkonstrain adalah = %1.3f\n',teta)

% menentukan jarak kritis

Xc=2*h*tan(teta);
fprintf('nilai jarak kritis inversi linier terkonstrain adalah = %1.3f\n',Xc)

% menentukan jarak crossover

Xcr=2*h*(sqrt(V2+V1/V2-V1));
fprintf('nilai crossover inversi linier terkonstrain adalah = %1.3f\n\n',Xcr)
Kode-kode program function penjalaran seismic refraksi

function Tr=seis_refrak(V1,V2)
% Ray tracing dan Travel time seismik refrkasi kritis yang
% Parameter model
% Satuan yang digunakan :: MKS
nx=1;ny=2;blx=1000;bly=50;npar=nx*ny;
% Skala sumbu horisontal
x=[0:(nx*blx/5):nx*blx]; % skala panjang bentangan survey
% Skala sumbu vertikal
y=0:bly:ny*bly; % skala kedalaman target bawah permukaan
% Sources and receivers posistions
ns=1; % number of source(s)
nr=24; % number of receiver(s)
xs=[0]; % x-Sources position
ys=zeros(1,length(xs)); % y-sources position
near_offset=100;
xr=linspace(near_offset,nx*blx,nr); % x-receivers position
yr=zeros(1,length(xr)); % y-receivers position
% Mendefinisikan kecepatan pada setiap lapisan
% vel=2500*ones(ny,nx);
vel=[V1 V2]';
% for k=1:ny
% vel(k,:)=vel(k,:)+(k-1)*2000;
% end
a=figure(1);
set(a,'Position',[1 1 1160 380],'NumberTitle','off','Name','Seismik Refraksi')
axes1 = axes('Parent',a,'Position',[0.07 0.12 0.85 0.82]);
hold(axes1,'on');
imagesc(x,-y,vel);
colorbar
title('Model Penjalaran Gelombang Seismik Refraksi');
xlabel('[m]');ylabel('[m]');
colormap(cool);
colorbar('vert');
axis xy;axis normal;hold on
set(gca,'YLim',[-ny*bly 0])
set(gca,'YTick',[-ny*bly:bly:0])
set(gca,'XLim',[0 nx*blx])
set(gca,'XTick',[0:(nx*blx/5):nx*blx])
vel=reshape(vel',1,ny*nx);
tt=[]; % waktu tempuh gelombang
h=bly;
tc=asind(vel(1)/vel(2)); % sudut kritis
T0=2*h*cosd(tc)/vel(1);
for m=1:ns
plot(xs(m),-ys(m),'*k');
for k=1:nr
q=(m-1)*nr+k;
plot(xr(k),-yr(k),'vr');
X=abs(xr(k)-xs(m));X1=h*tand(tc);X2=X-2*X1;
if xs(m)-xr(k)<~0
X1=X1;X2=X2;
else
X1=-1*X1;X2=-1*X2;
end
%....ploting titik P1
xp1=xs(m)+X1;yp1=ys(m)+h;
ap=[xs(m) xp1];bp=[ys(m) yp1];
plot(ap,-bp,'k');hold on
%....ploting titik P2
xp2=xs(m)+X1+X2;yp2=ys(m)+h;
ag=[ap xp2];bg=[bp yp2];
plot(ag,-bg,'k');hold on
%....ploting titik geophone
as=[ag xr(k)];bs=[bg yr(k)];
plot(as,-bs,'k'); hold on
drawnow
% Perhitungan travel time
x_off(q)=sqrt((yr(k)-ys(m))^2+(xr(k)-xs(m))^2); % offset
Tr(q)=T0+(x_off(q)/vel(2)); % Travel Time
end
end
 Kode-kode program model penjalaran gelombang inversi linier

clear all
clc

% model penjalaran gelombang seismik refraksi inversi linier

V1=2436.914;
V2=3500.849;
Tr=seis_refrak(V1,V2);

Kode-kode program model penjalaran gelombang inversi linier

terkonstrain
clear all
clc

% model penjalaran gelombang seismik refraksi inversi linier terkonstrain

% pada beta sama dengan 1
V1=2501.307;
V2=3375.995;
Tr=seis_refrak(V1,V2);

5.3. ANALISIS HASIL

5.3.1. Analisis hasil inversi linier dan inversi linier terkonstrain

hasil inversi linier adalah :

0.0411
0.0003

nilai kecepatan V2 inversi linier adalah = 3500.849

nilai kecepatan V1 inversi linier adalah = 2436.914
nilai sudut kritis inversi linier adalah = 1.560
nilai jarak kritis inversi linier adalah = 8865.831
nilai crossover inversi linier adalah = 3262.869

hasil damped least square terkonstrain dengan beta 0:0.1:1 adalah :

Columns 1 through 7
 0.0411 0.0410 0.0407 0.0403 0.0397 0.0390 0.0381
0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003

Columns 8 through 11

0.0372 0.0361 0.0350 0.0338

0.0003 0.0003 0.0003 0.0003

nilai kecepatan V2 inversi linier terkonstrain pada beta sama dengan 1 adalah =
3375.995
nilai kecepatan V1 inversi linier terkonstrain pada beta sama dengan 1 adalah =
2501.307
nilai sudut kritis inversi linier terkonstrain adalah = 1.482
nilai jarak kritis inversi linier terkonstrain adalah = 1117.781
nilai crossover inversi linier terkonstrain adalah = 2958.765

5.3.2. Analisis hasil travel time inversi linier dan inversi linier terkonstrain
Hasil travel time inversi linier
Tr =

Columns 1 through 7

0.0580 0.0692 0.0804 0.0916 0.1027 0.1139 0.1251

Columns 8 through 14

0.1363 0.1474 0.1586 0.1698 0.1810 0.1922 0.2033

Columns 15 through 21

0.2145 0.2257 0.2369 0.2480 0.2592 0.2704 0.2816

Columns 22 through 24

0.2928 0.3039 0.3151

 Hasil travel time inversi linier terkonstrain
Tr =

Columns 1 through 7

0.0565 0.0681 0.0797 0.0912 0.1028 0.1144 0.1260

Columns 8 through 14

0.1376 0.1492 0.1608 0.1724 0.1840 0.1956 0.2072

Columns 15 through 21

0.2187 0.2303 0.2419 0.2535 0.2651 0.2767 0.2883

Columns 22 through 24

0.2999 0.3115 0.3231

5.4. Gambar
5.4.1. Gambar grafik T Vs X
 5.4.2. Gambar penjalaran gelombang inversi linier

5.4.3. Gambar penjalaran gelomang inversi linier terkonstrain

5.4. PEMBAHASAN