notebook.community

Edit and run



In [1]:

    
using Seismic, SeismicPlotting









    



INFO: Loading help data...



In [2]:

    
download("http://seismic.physics.ualberta.ca/data/data_with_noise.su","data_with_noise.su");
SegyToSeis("data_with_noise.su","data_with_noise",["format"=>"su"]);
d,h = SeisRead("data_with_noise");
plotpar = ["style"=>"wiggles",
           "vmin"=>-5,"vmax"=>5,
           "aspect"=>"auto",
           "xlabel"=>"X","xunits"=>"","ox"=>0,"dx"=>1,
           "ylabel"=>"Time","yunits"=>"(seconds)","oy"=>0,"dy"=>h[1].d1,
           "cmap"=>"gray",
           "title"=>"data_with_noise.su"];
SeisPlot(d[:,:],plotpar);









    



  % Total    % Received % Xferd  Average Speed   Time    Time     Time  Current
                                 Dload  Upload   Total   Spent    Left  Speed
100  109k  100  109k    0     0  14.9M      0 --:--:-- --:--:-- --:--:-- 26.7M






    



nx=50
nt=501






    





Figure(PyObject <matplotlib.figure.Figure object at 0x11fc8c410>)



In [3]:

    
download("http://seismic.physics.ualberta.ca/data/syn_cmp.su","syn_cmp.su");
SegyToSeis("/Users/astanton/Seismic/data/syn_cmp.su","syn_cmp",["format"=>"su"]);
d,h = SeisRead("syn_cmp");
plotpar = ["style"=>"wiggles",
           "vmin"=>-2,"vmax"=>2,
           "aspect"=>"auto",
           "xlabel"=>"X","xunits"=>"","ox"=>0,"dx"=>1,
           "ylabel"=>"Time","yunits"=>"(seconds)","oy"=>0,"dy"=>h[1].d1,
           "cmap"=>"gray",
           "title"=>"syn_cmp.su"];
SeisPlot(d[:,:],plotpar);









    



  % Total    % Received % Xferd  Average Speed   Time    Time     Time  Current
                                 Dload  Upload   Total   Spent    Left  Speed
  0     0    0     0    0     0      0      0 --:--:-- --:--:-- --:--:--     0





    



nx=38
nt=522






    



100 88464  100 88464    0     0  13.5M      0 --:--:-- --:--:-- --:--:-- 28.1M






    





Figure(PyObject <matplotlib.figure.Figure object at 0x104344b10>)



In [4]:

    
download("http://seismic.physics.ualberta.ca/data/syn_cmp_mult.su","syn_cmp_mult.su");
SegyToSeis("syn_cmp_mult.su","syn_cmp_mult",["format"=>"su"]);
d,h = SeisRead("syn_cmp_mult");
plotpar = ["style"=>"overlay",
           "vmin"=>-1,"vmax"=>1,
           "aspect"=>"auto",
           "xlabel"=>"X","xunits"=>"","ox"=>0,"dx"=>1,
           "ylabel"=>"Time","yunits"=>"(seconds)","oy"=>0,"dy"=>h[1].d1,
           "cmap"=>"seismic",
           "title"=>"syn_cmp_mult.su"];
SeisPlot(d[:,:],plotpar);









    



  % Total    % Received % Xferd  Average Speed   Time    Time     Time  Current
                                 Dload  Upload   Total   Spent    Left  Speed
  0     0    0     0    0     0      0      0 --:--:-- --:--:-- --:--:--     0





    



nx=49
nt=1001






    



100  203k  100  203k    0     0  20.0M      0 --:--:-- --:--:-- --:--:-- 28.3M






    





Figure(PyObject <matplotlib.figure.Figure object at 0x121427b10>)



In [5]:

    
download("http://seismic.physics.ualberta.ca/data/small_stack.su","small_stack.su");
SegyToSeis("small_stack.su","small_stack",["format"=>"su"]);
d,h = SeisRead("small_stack");
nt = size(d,1);
dt = h[1].d1;
download("http://seismic.physics.ualberta.ca/data/wavelet_for_small_stack.su","wavelet_for_small_stack.su");
SegyToSeis("wavelet_for_small_stack.su","wavelet_for_small_stack",["format"=>"su"]);
w,h_wav = SeisRead("wavelet_for_small_stack");
nt_wav = size(w,1);
plotpar = ["style"=>"color",
           "vmin"=>-1,"vmax"=>1,
           "aspect"=>"auto",
           "xlabel"=>"X","xunits"=>"","ox"=>0,"dx"=>1,
           "ylabel"=>"Time","yunits"=>"(seconds)","oy"=>0,"dy"=>dt,
           "cmap"=>"PuOr","fignum"=>1,
           "title"=>"small_stack.su"];

SeismicPlotting.subplot(121);
SeisPlot(d[:,:],plotpar);

SeismicPlotting.subplot(122);
w = [w[:]' zeros(nt-nt_wav)'];
SeismicPlotting.plot(w[:],[0:nt-1]*dt);
SeismicPlotting.title("wavelet_for_small_stack.su");
SeismicPlotting.suptitle("Stack and Wavelet");
ax = SeismicPlotting.gca();
ax[:invert_yaxis](); 
SeismicPlotting.show();









    



  % Total    % Received % Xferd  Average Speed   Time    Time     Time  Current
                                 Dload  Upload   Total   Spent    Left  Speed
100 20800  100 20800    0     0  5366k      0 --:--:-- --:--:-- --:--:-- 19.8M






    



nx=20





    












    



nt=200
nx=1
nt=45






    



  % Total    % Received % Xferd  Average Speed   Time    Time     Time  Current
                                 Dload  Upload   Total   Spent    Left  Speed
105   420  105   420    0     0   142k      0 --:--:-- --:--:-- --:--:--  410k



In [6]:

    
download("http://seismic.physics.ualberta.ca/data/gom_cdp_nmo.su","gom_cdp_nmo.su");
SegyToSeis("gom_cdp_nmo.su","gom_cdp_nmo",["format"=>"su"]);
d,h = SeisRead("gom_cdp_nmo");
nx = length(h);
dx = abs(h[end].h - h[1].h)/nx;
plotpar = ["style"=>"color",
           "vmin"=>-5,"vmax"=>5,
           "aspect"=>5000,
           "xlabel"=>"Offset","xunits"=>"(ft)","ox"=>0,"dx"=>dx,
           "ylabel"=>"Time","yunits"=>"(seconds)","oy"=>0,"dy"=>h[1].d1,
           "cmap"=>"gray",
           "title"=>"gom_cdp_nmo.su"];
SeisPlot(d[:,:],plotpar);
SeismicPlotting.xticks(vec([0 4000 8000 12000]));









    



  % Total    % Received % Xferd  Average Speed   Time    Time     Time  Current
                                 Dload  Upload   Total   Spent    Left  Speed
  0     0    0     0    0     0      0      0 --:--:-- --:--:-- --:--:--     0





    



nx=92
nt=1751






    



100  650k  100  650k    0     0  29.8M      0 --:--:-- --:--:-- --:--:-- 33.4M






    





Figure(PyObject <matplotlib.figure.Figure object at 0x121babd10>)



In [7]:

    
download("http://seismic.physics.ualberta.ca/data/min_phase_wavelet.su","min_phase_wavelet.su");
SegyToSeis("min_phase_wavelet.su","min_phase_wavelet",["format"=>"su"]);
d,h = SeisRead("min_phase_wavelet");
nt = size(d,1);
dt = h[1].d1;

SeismicPlotting.plot(d[:],[0:nt-1]*dt,color="k");
SeismicPlotting.title("min_phase_wavelet.su");
ax = SeismicPlotting.gca();
ax[:invert_yaxis]();
SeismicPlotting.axis([-20,20,(nt-1)*dt,0.])
SeismicPlotting.ylabel("Time (s)");









    



  % Total    % Received % Xferd  Average Speed   Time    Time     Time  Current
                                 Dload  Upload   Total   Spent    Left  Speed
126   380  126   380    0     0   126k      0 --:--:-- --:--:-- --:--:--  371k






    



nx=1
nt=35



In [ ]: