seislab 3.02

最新matlab工具 解决 地震数据读取

% Examples4LogOnSeismicPlot
% Plot a log curve at a specified CDP on top of a seismic section

keep WF
presets

% global S4M

%     Get test data
seismic=s_data;
wlog=l_data;

%%     Convert log from depth to time
%     Add curve with two-way time （TWT）
wlog=l_depth2time（wlog{‘depth_time‘8000100}）; 
%     Remove all rows for which TWT has null values
wlog=l_rm_nulls（wlog‘anywhere‘‘twt‘）;
%      Make two-way time the first column
wlog=l_switch_depth（wlog‘twt‘）;
%     Select the desired curve （density）; remove all others
wlog=l_select（wlog{‘curves‘‘rho‘}）;
%     Create equi-distantly samples log curves in two-way travel time
step=2;      % Sample interval in ms of the depth-to-time converted well log
wlog=l_resample（wlogstep）;

%     Convert log structure into a seismic structure for the requested curve
tlog=s_log2seismic（wlog）;
%     Select the same time interval as the seismic
tlog=s_select（tlog{‘times‘seismic.firstseismic.last}{‘null‘NaN}）;
%     Remove the mean so that the curve will be plotted centered about the
%     trace number
tlog=tlog-nan_mean（tlog.traces）;

%%     Create a copy of the seismic dataset with all trace values set to
%      NaNs and sampled with the sample interval selected for the log curve;
slog=s_resample（seismicstep）;
slog.traces（1:end1:end）=NaN;

%     Find the trace number with the requested CDP （107） and copy the log curve 
%     to that trace; the only trace with valid values will be the one with the
%     log curve
loc=s_trace_numbers（seismic‘CDP‘107）; % Find the trace number with CDP == 107
slog.traces（:loc）=tlog.traces;         % Copy the log curve to that trace


%%     Plot seismic
s_wplot（seismic）                        % Plot the seismic in wiggle trace format

%     Plot log curve on top of the seismic data
s_wplot（slog{‘figure‘‘old‘}{‘wiggle‘‘red‘}{‘peak_fill‘‘‘} ...
   {‘interp‘‘linear‘}{‘wiggle_width‘1.5}{‘annotation‘‘CDP‘} ...
   {‘deflection‘1}）


%%     Plot seismic in color
s_cplot（seismic{‘shading‘‘interp‘}）    % Plot the seismic as a color image

%     Plot log curve on top of the seismic data
s_wplot（slog{‘figure‘‘old‘}{‘wiggle‘‘white‘}{‘peak_fill‘‘‘} ...
   {‘interp‘‘linear‘}{‘wiggle_width‘1.5}{‘annotation‘‘CDP‘} ...
   {‘deflection‘1}）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1005  2015-09-20 21:12  SeisLab 3.02\ACKNOWLEDGEMENT.txt
     文件        1548  2015-09-20 21:12  SeisLab 3.02\Description.txt
     文件        1861  2015-09-20 21:12  SeisLab 3.02\license.txt
     文件        1806  2015-09-20 21:12  SeisLab 3.02\README.txt
     文件        2334  2015-09-20 21:12  SeisLab 3.02\S4M\Examples\Examples4LogOnSeismicPlot.m
     文件        2048  2015-09-20 21:12  SeisLab 3.02\S4M\Examples\Examples4SeismicColorPlots.m
     文件        4765  2015-09-20 21:12  SeisLab 3.02\S4M\Examples\Examples4SeismicWigglePlots.m
     文件        3335  2015-09-20 21:12  SeisLab 3.02\S4M\Examples\Examples4WienerFilter.m
     文件        2030  2015-09-20 21:12  SeisLab 3.02\S4M\Examples\Log_examples1.m
     文件         759  2015-09-20 21:12  SeisLab 3.02\S4M\Examples\Log2seismic_examples1.m
     文件         165  2015-09-20 21:12  SeisLab 3.02\S4M\Examples\RunSeisLabHelpExamples.m
     文件        1820  2015-09-20 21:12  SeisLab 3.02\S4M\Examples\Seismic_examples1.m
     文件        1589  2015-09-20 21:12  SeisLab 3.02\S4M\Examples\Seismic_examples2.m
     文件        1525  2015-09-20 21:12  SeisLab 3.02\S4M\Examples\Seismic_examples3.m
     文件        4129  2015-09-20 21:12  SeisLab 3.02\S4M\Examples\Seismic_log_examples1.m
     文件        2314  2015-09-20 21:12  SeisLab 3.02\S4M\Examples\Seismic_phase_rotation.m
     文件        2328  2015-09-20 21:12  SeisLab 3.02\S4M\Examples\Seismic_principal_components.m
     文件        1392  2015-09-20 21:12  SeisLab 3.02\S4M\Examples\WF_Seislab_examples.m
     文件         398  2015-09-20 21:12  SeisLab 3.02\S4M\Geophysics_3.0\@struct\abs.m
     文件         394  2015-09-20 21:12  SeisLab 3.02\S4M\Geophysics_3.0\@struct\cumsum.m
     文件        1002  2015-09-20 21:12  SeisLab 3.02\S4M\Geophysics_3.0\@struct\diff.m
     文件        1371  2015-09-20 21:12  SeisLab 3.02\S4M\Geophysics_3.0\@struct\displayTools.m
     文件        1288  2015-09-20 21:12  SeisLab 3.02\S4M\Geophysics_3.0\@struct\double.m
     文件         680  2015-09-20 21:12  SeisLab 3.02\S4M\Geophysics_3.0\@struct\exp.m
     文件         740  2015-09-20 21:12  SeisLab 3.02\S4M\Geophysics_3.0\@struct\flipud.asv
     文件         854  2015-09-20 21:12  SeisLab 3.02\S4M\Geophysics_3.0\@struct\flipud.m
     文件         400  2015-09-20 21:12  SeisLab 3.02\S4M\Geophysics_3.0\@struct\imag.m
     文件         674  2015-09-20 21:12  SeisLab 3.02\S4M\Geophysics_3.0\@struct\log.m
     文件         389  2015-09-20 21:12  SeisLab 3.02\S4M\Geophysics_3.0\@struct\minus.m
     文件       11298  2015-09-20 21:12  SeisLab 3.02\S4M\Geophysics_3.0\@struct\moperation.m
     文件         212  2015-09-20 21:12  SeisLab 3.02\S4M\Geophysics_3.0\@struct\mpower.m
............此处省略407个文件信息