IODP Proceedings    Volume contents     Search

doi:10.2204/iodp.proc.344.202.2015

Methods

Permeability tests

Permeability tests were conducted using the Trautwein Soil Testing Equipment Company’s DigiFlow K (Fig. F3). 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 while an idealized solution of 33 g NaCl 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. F3).

Retrieved core samples from Expedition 344 were stored in plastic core liners and sealed in aluminum bags to prevent moisture loss. The sealed samples were stored in the refrigerator at 4°C until immediately prior to sample preparation. Most of the tests (17 samples) were conducted with flow in the vertical direction (along the axis of the core) using the whole-round core. A single minicore sample (344-U1414A-44R-1W, 53–55 cm) was oriented perpendicular to the core axis. The samples were carefully inspected for cracking or disturbance; two 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. Diameters of the trimmed whole-round cores ranged from 5.3 to 6.6 cm. Sample heights ranged from 4.45 to 13.7 cm. The sample was then placed in a rubber 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 either ~0.28 or ~0.41 MPa (40 or 60 psi, respectively).

Backpressure was achieved by concurrently ramping the cell and sample pressure to maintain a steady effective stress of 0.03 MPa (5 psi). Backpressure was maintained at least 24 h. Subsequently, cell fluid pressure was increased while 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. If the membrane is not pressed against the sample it would be possible for water from the flow lines to bypass the sample by running through the space between the membrane and the core. This would result in an erroneously large permeability estimate.

Because the whole-round samples were sealed immediately after cutting the core liner, the samples were expected to be near saturation prior to testing. Backpressuring at 0.28 MPa (40 psi) for ~24 h is sufficient to ensure full saturation under these conditions (ASTM International, 1990). A B-test on each sample was implemented for the last 12 samples tested. In a B-test, the cell confining pressure was instantaneously increased by 10 psi and the sample response was measured. For saturated samples, this yields the Skempton B-coefficient. For soft to medium clays, the Skempton B-coefficient is nearly 1 (Wang, 2000), and a B-test result >95% is typically used to indicate saturation. This criteria is not applicable for more consolidated materials, as compiled B-coefficients for mudstone, sandstone, and limestone are 95%, 50% to 88%, and 25%, respectively (Wang, 2000). Samples with B-coefficients below 95% were either given additional time for saturation or increased backpressure. Saturation was assumed if the B-value did not change with increased time or backpressure.

For each sample location, in situ effective stress was estimated for hydrostatic conditions using the shipboard bulk density measurements of the overlying sediments. Effective stress increments between depths of shipboard measurements were calculated assuming hydrostatic fluid pressures and summed. The estimated in situ effective stress is generally much greater than what was reached in the laboratory testing. Thus these permeability values should not be assumed to reflect in situ conditions. As noted above, the purpose of the measurements was to construct permeability–porosity relationships for use in fluid-flow modeling.

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. Since fluid pressure in the closed hydraulic system was affected by temperature changes, testing was conducted within a closed cabinet to keep the internal temperature uniform. Testing temperatures were 28° ± 1°C. Two 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 recording flow rates into and out of the sample.

The pressure difference, ΔP, from the top and bottom pumps was converted to hydraulic head difference (Δh):

Δh = ΔPfg,

where ρf is fluid density (1021 kg/m3) and g is acceleration due to gravity (9.81 m/s2).

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

Q = –KA(Δh/Δl),

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 = hydraulic head difference across the sample in meters, and
  • Δl = the length of the sample in meters.

The hydraulic conductivity values were then converted to permeability (in square meters) using

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

where

  • ρ= fluid density (1021 kg/m3),
  • g = acceleration due to gravity (9.81 m/s2), and
  • µ = viscosity (0.000893 Pa·s).

For laboratory temperature (averaging 28°C) and fluid salinity (33 g/L), a fluid density of 1021 kg/m3 and viscosity of 0.000893 Pa·s were calculated based on relationships compiled by Sharqawy et al. (2010). Assuming a reasonable water compressibility, density change caused by the applied pressure is minor (<0.1%). A 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 every sample, up to three effective stress steps were performed. 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. The volume change during consolidation is assumed to be solely due to changes in sample porosity. Influences from material and apparatus stiffness are assumed to be negligible. 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 n0 of the sample due to the small change in effective stress (0.03 MPa).

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

Vv(0) = n0VT(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)):

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 = (VT(1)VS)/VT(1).

Grain size analyses

Eleven of the 17 permeability samples were selected to represent various lithologies and were analyzed for grain size distribution. 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 in a solution of sodium hexametaphosphate (Calgon) to inhibit clay flocculation. Subsamples were also immersed in an ultrasonic bath for a minimum of 2 h to assist disaggregation. 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. Once disaggregated, a subsample was wet-sieved at 63 µm to determine the 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 emits X-rays that record the settling rates of particles suspended in a Calgon solution. The principle of Stoke’s law is then used to calculate grain sizes. 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 size particles.