example_nmo.m

著名的seismiclab的代码 是地震学研究人员必备的工具

%
%  example_nmo
%
%  Doing and Undoing NMO (constant velocity)
%
%  SeismicLab
%  Version 1
%
%  written by M.D.Sacchi, last modified Sept., 2001  
%  [e-mail]
%
%  Copyright (C) 1998 Seismic Processing and Imaging Group
%


% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Change to cdp.sun.su for big endians
% Read data
 
[D1,H]=readsegy('cdp.linux.su');
 
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%Get the offset and dt

offset = [H.offset];
dt = H(1).dt/1000/1000;

% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Apply NMO and NMO^-1 with constant v 
v = 3600;

D2 = cvnmo(D1,offset,dt,v, 1);
D3 = cvnmo(D2,offset,dt,v,-1);

% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


% Display

nt = max(size(D1));
t = 0:dt:(nt-1)*dt;

figure(1); clf;
subplot(131); wigb(D1,0.9,offset,t); title('Original')
subplot(132); wigb(D2,0.9,offset,t); title('NMO')
subplot(133); wigb(D3,0.9,offset,t); title('NMO \times NMO^{-1}')