Example Script for an Impulsive Source with Spiking Deconvolution

Copyright 1997, 1998, 1999, 2000, Infopipe Information Services Inc. Voice 303-973-7459 Fax 720-981-5668 Send email to info@infopipe.com
Example Script for an Impulsive Source with Spiking Deconvolution

@ D DERV$F.PRM ; Source Parameters File --- Derivatives - Signal
@ 1 DERVNF.PRM ; Source Parameters File --- Derivatives - Noise
@ Q Q$F.PRM ; Q-Filter Parameters File - Signal
@ I DFS5$F.PRM ; DFS5 Parameters File - Signal
@ G GEOP$F.PRM ; Geophone Parameters File - Signal
@ $ AS$.PRM ; Adjust Amplitude Parameters File - Signal
@ N ASN.PRM ; Adjust Amplitude Parameters File - Noise
@ X AUTO.PRM ; Autocorrelation Parameters File - Signal and Noise
@ Y DCN.PRM ; Deconvolution Parameters File
@ R ROSS.PRM ; Ross Weight Parameters File
@ B BLDSIG.BAT ; Batch File to construct Residual Wavelets and Filters
@ < IFFT.PRM ; Inverse FFT Parameters File
@ > FFT.PRM ; Forward FFT Parameters File
; ------------------------------------------------------------
; Copyright, 1998, 1999, 2000, Infopipe Information Services, Inc.
; darrell@infopipe.com
; ------------------------------------------------------------
; ---IMPULSIVE SOURCE MODEL-BASED WAVELET PROCESSING SCRIPT---
; ------------------------------------------------------------
; Comments: All comments begin with ";" and continue to end
; of line. All comments are stripped from output by list.exe.
; If a line begins with ";", then the line is ommitted from
; output.
; ------------------------------------------------------------
; ------------------------------------------------------------
; Tag-to-File specification. When |<tag> is encountered, list.exe
; ouputs following lines to <file name> until next |<tag> is
; encountered. NOTE: tags are single-character tags.
; Syntax:
;
;@ <tag1> <filename1>
;@ <tag2> <filename2>
;|<tag1>
;Contents to file 1
;|<tag2>
;Contents to file 2
; ------------------------------------------------------------
; ------------------------------------------------------------
; -NOTES FOR THIS PARTICULAR FILE-----------------------------
; In order to create a 'run', replace the variable "holders"
; below with the appropriate numbers.
;
; Run LIST.EXE, example: 'list myfile.fil'
;
; This will create all the input parameter files and
; the batch file to create the responses.
;
; List.exe strips all the comments from this file and
; produces the necessary files to create the responses.
;
; Next execute the batch file, example: 'bldsig'
; ------------------------------------------------------------
; -Batch File-------------------------------------------------
; ------------------------------------------------------------
; CONSTRUCT THE INDIVIDUAL FREQUENCY DOMAIN COMPONENTS
; ------------------------------------------------------------
|B
derv_f DERV$F.PRM DERV$F.F ;Derivatives(3)-Signal
derv_f DERVNF.PRM DERVNF.F ;Derivatives(1)-Noise
q_f Q$F.PRM Q$F.F ;Q-Filter-Signal
dfs5_f DFS5$F.PRM DFS5$F.F ;DFS 5
geop_f GEOP$F.PRM GEOP$F.F ;Geophone
; ------------------------------------------------------------
; COMBINE INTO RESPONSE
; ------------------------------------------------------------
; SIGNAL COMPONENT
f_mult_f DERV$F.F DFS5$F.F AUX1.F ;aux1=Derivatives*DFS 5
;f_mult_f DERV$F.F DFS5$.F AUX1.F ;aux1=Derivatives*DFS 5;From pmx_t
f_mult_f AUX1.F GEOP$F.F AUX2.F ;aux2=aux1*geophone
f_mult_f AUX2.F Q$F.F AUX$.F ;aux$=aux2*Q-Filter
; NOISE COMPONENT
f_mult_f DERVNF.F DFS5$F.F AUX3.F ;aux3=Derivatives*DFS 5
;f_mult_f DERVNF.F DFS5$.F AUX3.F ;aux3=Derivatives*DFS 5
f_mult_f AUX3.F GEOP$F.F AUXN.F ;auxN=aux3*geophone
;-------------------------------------------------------------
; CONSTRUCT THE RESIDUAL WAVELET
;-------------------------------------------------------------
; Adjust Signal to Noise Levels
f_as_f AUXN.F ASN.PRM ASN.F ;Adjust noise level
f_as_f AUX$.F AS$.PRM AS$.F ;Adjust signal level
; Construct time series
f_fft_f ASN.F IFFT.PRM ASN.FT ;Inverse FFT
f_fft_f AS$.F IFFT.PRM AS$.FT ;Inverse FFT
; Convert to real traces
f_real_t ASN.FT ASN.T ;Noise time series
f_real_t AS$.FT AS$.T ;Signal time series
; Construct autocorrelations
t_auto_t ASN.T AUTO.PRM AUTON.T ;Noise Autocorrelation
t_auto_t AS$.T AUTO.PRM AUTO$.T ;Signal Autocorrelation
; Add autocorrelations
t_sum_t AUTON.T AUTO$.T AUTOT.T ;Total Autocorrelation
; Perform deconvolution
t_dcn_t AUTOT.T DCN.PRM OPER.T ;Spiking Deconvolution
; Build residual wavelet time series
t_filt_t AS$.T OPER.T RES.T ;Residual Wavelet time series
; Shift Hilbert envelope peak to time zero
t_hshf_t RES.T SRES.T ;Place Peak energy at t=0
;-------------------------------------------------------------
; CONSTRUCT THE PHASE-ONLY CORRECTION FILTER
;-------------------------------------------------------------
; Phase-only inverse filter time series
t_poi_t SRES.T FINAL.T ;Phase only inverse
; Ross weight the filter
;t_ross_t POI.T ROSS.PRM FINAL.T ;Final Weighted Filter
; ProMax User-defined Series
t_pmx FINAL.T -.3 .3 FINAL.PMX ;ProMax User-defined Series
; Test operator
t_filt_t FINAL.T RES.T RESULT.T ;Test operator
;-------------------------------------------------------------
;-------------------------------------------------------------
;-------------------------------------------------------------
; AUXILLARY WAVELETS TO DISPLAY
; ------------------------------------------------------------
; CONVERT AUX3 TO TIME
; ------------------------------------------------------------
f_fft_f aux3.F ifft.prm aux3.tf
f_real_t aux3.tf aux3.t
; CONVERT AUX$ TO TIME
; ------------------------------------------------------------
f_fft_f aux$.F ifft.prm aux$.tf
f_real_t aux$.tf aux$.t
; ------------------------------------------------------------
; CONVERT AUX2 TO TIME
; ------------------------------------------------------------
f_fft_f aux2.F ifft.prm aux2.tf
f_real_t aux2.tf aux2.t
; ------------------------------------------------------------
; CONVERT AUX1 TO TIME
; ------------------------------------------------------------
f_fft_f aux1.F ifft.prm aux1.tf
f_real_t aux1.tf aux1.t
; ------------------------------------------------------------
; DERIVATIVES TO TIME
; ------------------------------------------------------------
f_fft_f DERV$F.F ifft.prm DERV$F.tf
f_real_t DERV$F.tf DERV$F.t
; ------------------------------------------------------------
; Q-FILTER TO TIME
; ------------------------------------------------------------
f_fft_f Q$F.F ifft.prm Q$F.tf
f_real_t Q$F.tf Q$F.t
; ------------------------------------------------------------
; DFS V INSTRUMENT TO TIME
; ------------------------------------------------------------
f_fft_f DFS5$F.F ifft.prm DFS5$F.tf
f_real_t DFS5$F.tf DFS5$F.t
f_db_t DFS5$F.F DFS5$f.A ;Amplitude
;f_fft_f DFS5$.F ifft.prm DFS5$.ttf
;f_real_t DFS5$.ttf DFS5$.t
;f_db_t DFS5$.F DFS5$f.A ;Amplitude
; ------------------------------------------------------------
; GEOPHONE TO TIME
; ------------------------------------------------------------
f_fft_f GEOP$F.F ifft.prm GEOP$F.tf
f_real_t GEOP$F.tf GEOP$F.t
; ------------------------------------------------------------
; -Parameters File--------------------------------------------
; ------------------------------------------------------------
; FORWARD FFT PARAMETERS FILE
; ------------------------------------------------------------
|>
1 ; forward parameter
; ------------------------------------------------------------
; -Parameters File--------------------------------------------
; ------------------------------------------------------------
; INVERSE FFT PARAMETERS FILE
; ------------------------------------------------------------
|<
-1 ; inverse parameter
; ------------------------------------------------------------
; -Parameters File--------------------------------------------
; ------------------------------------------------------------
; AUTOCORRELATION PARAMETERS FILE
; ------------------------------------------------------------
|X
0.128 ;Autocorrelation Lag (s)
0.128 ;Output Length (s)
0 ;Normalize = 1, Leave alone 0
; ------------------------------------------------------------
; -Parameters File--------------------------------------------
; ------------------------------------------------------------
; DECONVOLUTION PARAMETERS FILE
; ------------------------------------------------------------
|Y
0.256 ;Autocorrelation length (two-side) (s)
0.128 ;Operator length (s)
0.1 ;White Light (%)
; ------------------------------------------------------------
; -Parameters File--------------------------------------------
; ------------------------------------------------------------
; NOISE AMPLITUDE SPECIFICATION PARAMETERS FILE
; ------------------------------------------------------------
|N
40 ;Frequency (Hz) at which to normalize series
-30 ;Noise level (dB), i.e., -30.0
; ------------------------------------------------------------
; -Parameters File--------------------------------------------
; ------------------------------------------------------------
; SIGNAL AMPLITUDE SPECIFICATION PARAMETERS FILE
; ------------------------------------------------------------
|$
40 ;Frequency (Hz) at which to normalize series
0 ;Signal level (dB), i.e., 0.0
; ------------------------------------------------------------
; -Parameters File--------------------------------------------
; ------------------------------------------------------------
; ROSS WEIGHT PARAMETERS FILE
; ------------------------------------------------------------
|R
;_START1_,_END1_ ;Start and End times for on-ramp
-0.2,0.0
;_START2_,_END2_ ;Start and End times for off-ramp
0.0,0.2
;_M1_,_N1_ ;Ross weights for on-ramp
2,3
;_M2_,_N2_ ;Ross weights for off-ramp
2,3
; ------------------------------------------------------------
; -Parameters File--------------------------------------------
; ------------------------------------------------------------
; IMPULSIVE SOURCE SECTION - SOURCE PARAMETERS FILE
; ------------------------------------------------------------
|D
0.15 ;Time Shift (s)
0.002 ;Sample Interval (s)
500 ;Number of samples will be rounded to pow of 2 up
3 ;Number of Derivatives
; ------------------------------------------------------------
; -Parameters File--------------------------------------------
; ------------------------------------------------------------
; IMPULSIVE SOURCE SECTION - NOISE PARAMETERS FILE
; ------------------------------------------------------------
|1
0.15 ;Time Shift (s)
0.002 ;Sample Interval (s)
500 ;Number of samples will be rounded to pow of 2 up
1 ;Number of Derivatives
; ------------------------------------------------------------
; -Parameters File--------------------------------------------
; ------------------------------------------------------------
; EARTH FILTER SECTION - Q FILTER PARAMETERS FILE
; ------------------------------------------------------------
|Q
1 ;Time (s)
0.002 ;Sample Interval (s)
500 ;Number of samples will be rounded to pow of 2 up
30 ;Q (1/s)
10000 ;Pedestal Frequency (Hz)
1 ;1-Forward 0-Inverse
0.001 ;White Light
; ------------------------------------------------------------
; -Parameters File--------------------------------------------
; ------------------------------------------------------------
; DFS V INSTRUMENT FILTER PARAMETERS FILE
; ------------------------------------------------------------
|I
500 ;Number of samples will be rounded to pow of 2 up
0.002 ;Sample Interval (s)
40 ;Lo Cut Cutoff (Hz)
2 ;1-18dB/Oct 2-12dB/Oct 3-36dB/Oct 4-24dB/Oct
180 ;Hi Cut Cutoff (Hz)
0 ;Notch 0-None 1-60 Hz 2-50 Hz
40 ;0-dB Normalization Frequency (Hz)
; ------------------------------------------------------------
; -Parameters File--------------------------------------------
; ------------------------------------------------------------
; GEOPHONE RESPONSE PARAMETERS FILE
; ------------------------------------------------------------
|G
500 ;Number of samples will be rounded to pow of 2 up
0.002 ;Sample Interval (s)
28 ;Natural Frequency (Hz)
0.67 ;Damping Factor
; ------------------------------------------------------------
; --END OF SCRIPT FILE----------------------------------------
; -----------------------------------------------------