IODP Proceedings    Volume contents     Search

doi:10.2204/iodp.proc.320321.217.2014

Methods

Permeability tests

Permeability tests were conducted using the Trautwein Soil Testing Equipment Company’s DigiFlow K (Fig. F2). The equipment consists of a cell (to contain the sample and provide isostatic effective stress) and three pumps (sample top pump, sample bottom pump, and cell pump). American Society for Testing and Materials (ASTM) designation D5084-90 (ASTM International, 1990) was used as a guideline for general procedures. Deionized water was used as the fluid in the pumps, whereas an idealized solution of seawater (25 g NaCl and 8 g MgSO4 per liter of water) permeates the sample. Pressure is transmitted from the deionized water to the permeant across a rubber membrane in an interface chamber (Fig. F2).

The retrieved core samples from Expeditions 320 and 321 were stored in plastic core liners and sealed in polyester film bags to prevent moisture loss. The sealed samples were stored in the refrigerator at 4°C until immediately prior to sample preparation. All testing was conducted with flow in the vertical direction (along the axis of the core). The permeability testing apparatus accommodates the whole-round core. As a result, disturbance of the sample is minimal relative to testing in which plugs or subsamples are removed from the core. The samples were carefully inspected for cracking or disturbance. Fractured or damaged samples were not tested. To provide freshly exposed surfaces, cores were trimmed on both ends using an Exacto knife, wire saw, or utility knife, depending on core properties. After trimming the ends of the sample, the diameter and height of the sample were measured. The sample was then placed in a flexible wall membrane and fitted with saturated porous disks and end caps. Next, the sample was placed in the cell, which was filled with deionized water. The membrane-encased sample was subjected to the applied pressure of the water in the cell, but fluid exchange occurs only through the flow lines connecting the end caps to the top and bottom pumps. A small confining pressure of ~0.03 MPa (5 psi) was applied, and flow lines were flushed to remove any trapped air bubbles. After flushing the flow lines, the sample was backpressured to ~0.28 MPa (40 psi) in order to ensure full saturation of the sample. Backpressure was achieved by concurrently ramping the cell pressure and the sample pressure to maintain a steady effective stress of 0.03 MPa. Once the sample reached saturation, the cell fluid pressure was increased while the sample backpressure was maintained, thus increasing the effective stress on the sample. This effective stress both consolidates the sample and pushes the flexible membrane against the sample to prevent flow bypassing the sample. For comparison with applied effective stresses, the in situ effective stress was computed from the shipboard density measurements (Table T1).

Once the target effective stress was achieved, cell pressure and backpressure were maintained. The sample was allowed to equilibrate for at least 12 h and generally 24 h. Throughout testing, inflows and outflows to the cell fluid were monitored to assess changes in sample volume, and sample data were recorded every minute. Because fluid pressure in the closed hydraulic system was affected by temperature changes, testing was conducted within a closed cabinet with a fan to keep the internal temperature uniform. Testing temperatures were 28°C ± 1°C. Three or more flow tests were performed at each effective stress level, with flow direction varied between tests. Flow tests were run by specifying pressures of the top and bottom pump and allowing flow rates into and out of the sample to equilibrate with time. Equilibration was indicated by inflow equaling outflow.

We used the measured flow rate, cross-sectional area of the sample, and the head difference between the top and bottom of the sample to calculate the hydraulic conductivity using Darcy’s law,

Q = –KAhl),

where

  • Q = measured flow rate (in cubic meters per second),
  • K = hydraulic conductivity (in meters per second),
  • A = the cross-sectional area of the sample (in square meters),
  • Δh = the difference in head across the sample (in meters), and
  • Δl = the length of the sample (in meters).

The hydraulic conductivity values were then converted to permeability (k, in square meters) using the following equation:

k = (Kµ)/(ρg),

where

  • ρ = fluid density (1023 kg/m3),
  • g = the gravitational constant (9.81 m/s2), and
  • µ = viscosity (0.000857 Pa·s).

The density value was estimated based on literature values for a temperature of 30°C and a salinity of 33 kg/m3 (Chemical Rubber Company, 2012). Assuming a reasonable water compressibility, density change due to the applied pressure is minor (<0.1%). The viscosity value was obtained from a synthesis of previous relationships (Sharqawy et al., 2010) for water at 30°C and 33 kg/m3 salinity. A 30 min or 1 h interval of stable flow rates was averaged for the permeability calculations, and the standard deviation of the permeability during that interval was calculated to assess uncertainty. The fluctuations in the calculated permeability are likely caused by slight temperature variations. The resulting volume changes would cause temporary changes in measured flow rates. The time interval was selected based on where inflow best matched outflow, indicating steady-state conditions and where the standard deviation was minimized.

For any given sample, permeability was determined at one, two, or three effective stress steps. The corresponding porosity for each effective stress was calculated using the change in volume of fluid (mL) contained in the cell during each consolidation step. Total sample volume (VT(0)) was calculated using πr2h, where r is the radius of the core sample and h is the height of the sample. Initial porosities (n0) for volume calculations were obtained from the shipboard moisture and density results of samples taken adjacent to each permeability sample. We assumed that the porosity of the sample at the end of backpressure is similar to the initial porosity (n0) of the sample because of the small change in effective stress (0.03 MPa).

Using the initial porosity (n0), the volume of voids before the testing (Vv(0)) was calculated using

Vv(0) = n0 × VT(0).

Volume of solids (Vs) was calculated using

Vs = VT(0)Vv(0).

Using the difference of cell volumes between two consecutive steps (e.g., cell volume at backpressure and cell volume at first consolidation), the change in volume of water in the cell (ΔVT(1)) was calculated. The new total volume of the sample (VT(1)) after pore spaces were reduced during the consolidation process was determined by subtracting the change in cell volume at the end of the consolidation step (ΔVT(1)) from the total sample volume (VT(0)) using

VT(1) = VT(0) – ΔVT(1).

Using the calculated new total volume of the sample (VT(1)), the new porosity at the end of the consolidation is calculated. The new porosity (n1) at the end of the consolidation is

n1 = (1 – Vs)/VT(1).

Grain size analyses

Subsamples for quantitative grain size analyses were extracted in 1.5 cm thick intervals from the sample after completion of the permeability tests. The subsamples were homogenized and disaggregated. A small aliquot of the homogenized sample was dried to determine water content, which was then used to establish the equivalent dry mass used in the particle size analysis. A sample was disaggregated and wet sieved at 63 µm to determine sand-sized fraction. A separate sample was wet sieved at 53 µm, and material <53 µm was analyzed on a 5100 Micrometrics SediGraph (Coakley and Syvitski, 1991). The SediGraph data were combined with the wet sieve results to normalize the mud and sand fraction to their relative masses to determine the proportion of sand-, silt-, and clay-sized particles.