Standard Wireline Sonic Waveform Data
IODP -USIO logging contractor: LDEO-BRG
Hole: U1301B
Expedition: 301
Location: Juan de Fuca Ridge (NE Pacific)
Latitude: 47º 45.229' N
Longitude: 127º 45.826' W
Logging date: August 2004
Sea floor depth 2667.84 mbrf
Total penetration: 582.80 mbsf
Total core recovered: 69.05 m (29.81% of cored section)
Oldest sediment cored: no core taken in sediments.
Lithologies: basalt-hyaloclastite breccia, massive and pillow basalts
TOOL USED: DSI (Dipole Sonic Imager)
Recording mode: Monopole P&S, Cross-Dipole, and Stoneley mode (pass 1 and 2), Monopole P&S, Stoneley , and Upper and Lower Dipole mode (casing) .
Remarks about the recording: none.
MONOPOLE P&S MODE: measures compressional and hard-rock shear slowness. The monopole transmitter is excited by a high-frequency pulse, which reproduces conditions similar to previous sonic tools.
UPPER DIPOLE MODE: measures shear wave slowness using firings of the upper dipole transmitter.
LOWER DIPOLE MODE: measures shear wave slowness using firings of the lower dipole transmitter.
CROSS-DIPOLE MODE: uses alternate firings of upper and lower dipole transmitter, thus allowing acquisition of orthogonally polarized data for anisotropy studies.
STONELEY MODE: measures low-frequency Stoneley wave slowness. The monopole transmitter, driven by a low-frequency pulse, generates the Stoneley wave.
Acoustic data are recorded in DLIS format. Each of the eight waveforms generally consists of 512 samples (256 in cross-dipole mode), each recorded every 10 (monopole P&S) and 40 microsec (all dipole modes), at depth intervals of 15.24 cm (6 inches). The original data in DLIS format is first loaded on a Sun system using GeoFrame software. The packed waveform data files are then converted into ASCII and finally binary format.
Each line is composed of the entire waveform set recorded at each depth, preceded by the depth (multiplied by 10 to be stored as an integer). In the general case of 8 waveforms with 512 samples per waveform, this corresponds to 1 + 4x512 = 4097 columns. In this hole, the specifications of the files are:
Number of columns: 4097
Number of rows: 828 (Pass 1)
Number of columns: 4097
Number of rows: 760 (Pass 2)
Number of columns: 4097
Number of rows: 1123 (Through casing)
NOTE: the number of columns is 4098 in the case of the cross-dipole.
All values are stored as 'signed short integers' (= 16 bits signed).
Each file can be viewed directly as an image with Spyglass transform or NIH image on MAcs and PCs, or Fortner transform or Khoros on UNIX) - where it has to be opened as either 'raw' or 'binary matrix'. Any image or signal-processing package should also allow viewing it.
The following files have been loaded:
DSI from FMS/DSI/GPIT/SGT (Pass 1, casing at ~3014 mbrf)
U1301B-cdip_ldip_cl_p1.bin: 2968.3-3094.3 mbrf
U1301B-cdip_ldip_il_p1.bin: 2968.3-3094.3 mbrf
U1301B-cdip_udip_cl._p1.bin: 2968.3-3094.3 mbrf
U1301B-cdip_udip_il_p1.bin: 2968.3-3094.3 mbrf
U1301B-mono_p1.bin: 2968.3-3094.3 mbrf
U1301B-st_p1.bin: 2968.3-3094.3 mbrf
DSI from FMS/DSI/GPIT/SGT (Pass 2, casing at ~3014 mbrf)
U1301B-cdip_ldip_cl_p2.bin: 2978.6-3094.3 mbrf
U1301B-cdip_ldip_il_p2.bin: 2978.6-3094.3 mbrf
U1301B-cdip_udip_cl._p2.bin: 2978.6-3094.3 mbrf
U1301B-cdip_udip_il_p2.bin: 2978.6-3094.3 mbrf
U1301B-mono_p2.bin: 2978.6-3094.3 mbrf
U1301B-st_p2.bin: 2978.6-3094.3 mbrf
DSI from FMS/DSI/GPIT/SGT (through casing)
U1301B-udip_csg.bin: 2864.2-3035.2 mbrf
U1301B-ldip_csg.bin: 2864.2-3035.2 mbrf
U1301B-mono_csg.bin: 2864.2-3035.2 mbrf
U1301B-st_csg.bin: 2864.2-3035.2 mbrf
The sonic waveform files are not depth shifted to a reference run or to the seafloor. For depth shift to the sea floor, please refer to the DEPTH SHIFT section in the standard log documentation file.
NOTE: For users interested in reading and converting the data to a format more suitable for their own purpose, the fortran declaration used to open the file *.bin would be:
open (1, file = *.bin,access = 'direct', recl = 2*(1+nwav*nsamples))
where nwav is the number of waveforms (4 or 8) and nsamples the number of samples per waveforms (see above).
If the total number of depths where waveforms have been recorded is ndepth (for a 150 m interval with data every 15cm, ndepth would be 1000), a generic loop to read the data would be
do k = 1, ndepth
...
read(1, rec=k) iz, (data(i), i = 1, nwav*nsamples)
z(k)=float(iz)/10
...
enddo
In such a case, the file dimensions will be ndepth rows by nwav*nsamples+1 columns
For further information about the logs please contact:
Cristina Broglia
Phone: 845-365-8343
Fax: 845-365-3182
E-mail: Cristina Broglia