Synthetic_Signal.rar_synthetic_初至检测_地震信号_地震模型_微地震信号

%% 微地震信号检测-初至模型生成 % 地震道最后一个数为地震记录信号的初值时间（初至时间存到地震道最后一个值中） close all clear clc %% 1.0 输入控制参数 m=2000; % 道长 n=15 ; % 道数 snr=20 ; % 信噪比 maxNum=500; % 生成模型个数 R=zeros(m+1,n); for j=1:maxNum for i=1:n T1=200+fix(abs(5*rand(1))); T2=220+fix(abs(5*rand(1))); TN=fix(50*abs(rand(1))); T1=T1+TN; T2=T2+TN; % 模型 r=zeros(m,1); r(T1)=0.6+0.4*abs(rand(1)); r(T2)=0.6+0.4*abs(rand(1)); seismic=s_data;seismic.step=1; minwav=s_seismic2wavelet(seismic,{'color','blue'},{'type','min-phase'},{'wlength',200}); % s_wplot(minwav) w=minwav.traces; % figure;plot(w); d=conv(r,w,'same'); d(abs(d(1:T1))<10^(-3)*max(abs(d)))=0; out = awgn(d',snr); % figure;plot(d); d_NorZero=find(d,1,'first'); R(:,i)=[out(:);d_NorZero]; end %% 2.0 写数据 write_segy_file(s_convert(R,0,1,[],'ms'),['seismic_',num2str(j)','.segy']); % figure;imagesc(R(1:end-1,:)); % figure;plot(R(1:end-1,:)); end