function plot_segy(tr,SI,Tn,k)
%input:
%   tr:地震数据
%   SI:采样点数
%   Tn:地震道集道数
%   k:地震道集序号
%email:2906310084@qq.com
%2020-6-4
    x=1:SI;
    figure
    for i=(k-1)*Tn+1:Tn*k
        for n=1:size(tr,2)
            traces(:,n)=tr(:,n)/max(abs(tr(:,n)));
        end
    plot(traces(:,i)*0.5+i,x,'b');
    fill(traces(:,i)*0.5+i,x,'b');hold on
    end
    xlabel(['chen ',num2str(k)])
    ylabel('t/µs')
    set(gca,'Ydir','reverse');%翻转
    set(gca,'XAxislocation','top');
    axis([Tn*(k-1),Tn*k+1,0,2000])
end