seislab 3.02

% 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.02ACKNOWLEDGEMENT.txt
     文件        1548  2015-09-20 21:12  SeisLab 3.02Description.txt
     文件        1861  2015-09-20 21:12  SeisLab 3.02license.txt
     文件        1806  2015-09-20 21:12  SeisLab 3.02README.txt
     文件        2334  2015-09-20 21:12  SeisLab 3.02S4MExamplesExamples4LogOnSeismicPlot.m
     文件        2048  2015-09-20 21:12  SeisLab 3.02S4MExamplesExamples4SeismicColorPlots.m
     文件        4765  2015-09-20 21:12  SeisLab 3.02S4MExamplesExamples4SeismicWigglePlots.m
     文件        3335  2015-09-20 21:12  SeisLab 3.02S4MExamplesExamples4WienerFilter.m
     文件        2030  2015-09-20 21:12  SeisLab 3.02S4MExamplesLog_examples1.m
     文件         759  2015-09-20 21:12  SeisLab 3.02S4MExamplesLog2seismic_examples1.m
     文件         165  2015-09-20 21:12  SeisLab 3.02S4MExamplesRunSeisLabHelpExamples.m
     文件        1820  2015-09-20 21:12  SeisLab 3.02S4MExamplesSeismic_examples1.m
     文件        1589  2015-09-20 21:12  SeisLab 3.02S4MExamplesSeismic_examples2.m
     文件        1525  2015-09-20 21:12  SeisLab 3.02S4MExamplesSeismic_examples3.m
     文件        4129  2015-09-20 21:12  SeisLab 3.02S4MExamplesSeismic_log_examples1.m
     文件        2314  2015-09-20 21:12  SeisLab 3.02S4MExamplesSeismic_phase_rotation.m
     文件        2328  2015-09-20 21:12  SeisLab 3.02S4MExamplesSeismic_principal_components.m
     文件        1392  2015-09-20 21:12  SeisLab 3.02S4MExamplesWF_Seislab_examples.m
     文件         398  2015-09-20 21:12  SeisLab 3.02S4MGeophysics_3.0@structabs.m
     文件         394  2015-09-20 21:12  SeisLab 3.02S4MGeophysics_3.0@structcumsum.m
     文件        1002  2015-09-20 21:12  SeisLab 3.02S4MGeophysics_3.0@structdiff.m
     文件        1371  2015-09-20 21:12  SeisLab 3.02S4MGeophysics_3.0@structdisplayTools.m
     文件        1288  2015-09-20 21:12  SeisLab 3.02S4MGeophysics_3.0@structdouble.m
     文件         680  2015-09-20 21:12  SeisLab 3.02S4MGeophysics_3.0@structexp.m
     文件         740  2015-09-20 21:12  SeisLab 3.02S4MGeophysics_3.0@structflipud.asv
     文件         854  2015-09-20 21:12  SeisLab 3.02S4MGeophysics_3.0@structflipud.m
     文件         400  2015-09-20 21:12  SeisLab 3.02S4MGeophysics_3.0@structimag.m
     文件         674  2015-09-20 21:12  SeisLab 3.02S4MGeophysics_3.0@structlog.m
     文件         389  2015-09-20 21:12  SeisLab 3.02S4MGeophysics_3.0@structminus.m
     文件       11298  2015-09-20 21:12  SeisLab 3.02S4MGeophysics_3.0@structmoperation.m
     文件         212  2015-09-20 21:12  SeisLab 3.02S4MGeophysics_3.0@structmpower.m
............此处省略407个文件信息