In Situ Methane Hydrate Morphology Investigation: Natural Gas

Subscriber access provided by University of Sussex Library

Article

In situ Methane Hydrate Morphology Investigation: Natural Gas Hydratebearing Sediment Recovered from the Eastern Nankai Trough Area Yusuke Jin, Yoshihiro Konno, Jun Yoneda, Masato Kida, and Jiro Nagao Energy Fuels, Just Accepted Manuscript • DOI: 10.1021/acs.energyfuels.6b00762 • Publication Date (Web): 26 Jun 2016 Downloaded from http://pubs.acs.org on July 3, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Energy & Fuels is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

In situ Methane Hydrate Morphology Investigation: Natural Gas Hydrate-bearing Sediment Recovered from the Eastern Nankai Trough Area Yusuke Jin,1* Yoshihiro Konno, 1 Jun Yoneda, 2 Masato Kida,1 and Jiro Nagao1

1

Methane Hydrate Production Technology Research Group, Research Institute of Energy

Frontier, Department of Energy and Environment, National Institute of Advanced Industrial Science and Technology (AIST), Tsukisamu-Higashi, Toyohira-Ku, Sapporo 062-8517, Japan 2

Methane Hydrate Geo-mechanics Research Group, Research Institute of Energy Frontier,

Department of Energy and Environment, National Institute of Advanced Industrial Science and Technology (AIST), 16-1 Onogawa, Tsukuba, Ibaraki, 305-8569, Japan

*Corresponding author Email: [email protected] Address: 2-17-2-1, Tsukisamu-Higashi, Toyohira-Ku, Sapporo 062-8517, JAPAN Tel: +81-11-857-8526, Fax: +81-11-857-8417

ACS Paragon Plus Environment

1

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 33

ABSTRACT

The hydrate morphology of natural gas hydrate-bearing (GH) sediments recovered from the eastern Nankai trough area was investigated under hydrostatic pressurized conditions that prevent dissociation of gas hydrates in a sediment. We developed a novel X-ray computed tomography (CT) system and an attenuated total reflection infrared (ATR-IR) probe for use in the Instrumented Pressure Testing Chamber (IPTC) for our set of Pressure-core Nondestructive Analysis Tools (PNATs), which can measure the sediment structure, primary wave velocity (PWV), density, and shear strength under pressurized conditions. The hydrate saturation values estimated using the ATR-IR absorption bands of H2O molecules strongly correlate with PWV. Assuming homogeneity of hydrate distribution in the planes perpendicular to the sample depth direction, the hydrate morphology of natural GH sediments in the eastern Nankai trough area demonstrated a load-bearing morphology type. The predicted hydrate morphology results are in good agreement with data reported in the literature. The combination of PNATs including ATRIR spectroscopy can be used to estimate the properties of GH sediments without the release of pressure to atmospheric conditions in order to model gas hydrate reservoirs for natural gas production.

KEYWORDS natural gas; clathrate; methane; spectroscopy


2

Page 3 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

INTRODUCTION Gas clathrate hydrates (gas hydrates) are clathrate compounds that crystallize by enclosing guest molecules in H2O structures.1 Because gas hydrates are formed under high-pressure, lowtemperature conditions, natural gas hydrates are found in oceanic and permafrost sediments.1-2 Guest molecules in natural gas hydrates are almost exclusively methane (CH4); CH4 obtained from natural gas hydrates is considered to be a potential new unconventional energy resource. Gas production behavior from gas hydrates has already been studied under laboratory conditions. Some methods of gas production from natural gas hydrate reservoirs have recently been reported using large-scale apparatuses.3-6 Several gas production tests on natural gas hydrate reservoirs have been performed. The first test was conducted at the Messoyakha sites (Russia) in 1970,7-8 following which other tests were conducted at the permafrost fields of the Mallik sites (Canada) and the North Slope of Alaska (USA) in 2002 and 2007 and 2011/2012, respectively.9-11 The latest test was performed in 2013 at an offshore field in the eastern Nankai area (Japan).12 The knowledge obtained from natural gas hydrate-bearing (GH) sediments is crucial for understanding natural GH deposit layers and for modeling gas hydrate reservoirs. The coring of natural GH sediments has generally been conducted using a core sampler that maintains in situ pressure conditions in order to avoid dissociation of gas hydrates in the recovered cores. 13-16 Until 2004, after the release of the core sampler pressure to atmospheric conditions, the natural GH sediments were treated with liquid nitrogen (LN2) and were analyzed in onshore laboratories. However, the sediment structure of natural GH sediments is damaged by this pressure release and the LN2 treatment.17-20 Natural GH sediments damaged in such a manner can lead to over-/under-estimation and misunderstanding of sediment properties.21-23


3

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 33

More recently, pressurized-core analysis systems that can measure pressurized core samples without pressure release and LN2 treatment have been developed.24-26 The Pressure Core Characterization Tools (PCCTs) system, which was established by the Georgia Institute of Technology (Georgia Tech), is an innovative system for measuring the properties of natural GH sediments under pressurized conditions.26 The Instrumented Pressure Testing Chamber (IPTC) included in the PCCTs is a pressure vessel equipped with multiple sensor probe ports that can measure the primary wave (P-wave) velocity (PWV), secondary wave (S-wave) velocity (SWV), shear strength, and resistivity of pressurized GH sediments using various sensors through the probe ports.25, 27-28 The properties obtained from the IPTC are crucial for a better understanding of natural GH sediments and the characterization of gas hydrate reservoirs. The Research Consortium for Methane Hydrate Resources in Japan (MH21),29 which initiated Japan’s Methane Hydrate R&D Program (which is now managed by the Ministry of Economy, Trade, and Industry (METI) ) is currently planning a second offshore gas production test in the eastern Nankai area. To aid in this effort, MH21 is developing its own IPTC-based pressurized-core analysis system in strong collaboration with Georgia Tech and the U.S. Geological Survey (USGS).

The permeability, stiffness, and thermal conductivity of GH sediments are important characteristics determining the gas production from a gas hydrate reservoir. These properties strongly relate to the gas hydrate morphology in sandy sediments, which is mainly classified into three models: cementing, grain-coating, and load-bearing.30 The P-wave traveling through the GH sediment is one indicator of the stiffness of sediment samples. PWV is determined in part by the hydrate saturation, Sh (%), which is the volume ratio of hydrates in the pores.31-34 As gas hydrate morphology can affect the PWV behavior through alterations in the Sh,35 the


4

Page 5 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

morphology can be predicted through the correlation between PWV and Sh. At present, PWV can be measured under pressurized conditions using the IPTC system, whereas Sh in natural GH sediments is estimated directly from the gas volume released by hydrate dissociation and by the hydrate number of the gas hydrate.18, 36-37 Since there are few methods for Sh measurement under pressurized conditions, PWV–Sh correlations in pressurized GH sediments have not yet been established.

Although the absorption of the O–H stretching IR bands of H2O are too strong to measure via a normal IR systems, the O–H stretching infrared (IR) band of H2O molecules measured by attenuated total reflection–IR (ATR-IR) is suitable for monitoring solid–liquid water mixture. Jin et al.38-39 observed hydrate formation behavior and ice formation behavior in pore spaces during hydrate dissociation from the variation of the O–H stretching bands using an ATR-IR system. ATR-IR measurements could be productive for in situ observation of gas hydrates coexisting with liquid H2O, i.e., the pore spaces in GH sediments. In this study, we developed an ATR-IR probe designed for use in the IPTC system to measure the Sh of GH sediments under pressurized conditions. We also developed a large-scale X-ray computed tomography (CT) system for longlength core imaging and a bulk density measurement system for use under pressurized conditions. Furthermore, we demonstrated the nondestructive evaluation of gas hydrate morphology in pressurized natural GH sediments recovered from the eastern Nankai area using this pressurizedcore analysis system.

EXPERIMENTAL SECTION


5

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 33

Natural Gas Hydrate-bearing Sediment. The natural GH sediment sample used in this study was recovered from a coring well (hereafter AT1-C) in the eastern Nankai area during the “Daini-Atusmi Knoll” drilling campaign of June–July 2012.40 Core recovery was carried out by the deep sea drilling vessel CHIKYU (owned by the Center for Deep Earth Exploration, CDEX, of the Japan Agency for Marine Earth Science and Technology, JAMSTEC). In order to avoid gas hydrate dissociation and maintain the sediment structure, the core samples were recovered using a Hybrid PCS core sampler.16 The natural GH sediment samples were transferred to our laboratory (Sapporo, Japan), and stored at 278 K and 15 MPa (hydrostatic pressure). From the logging wired drilling (LWD) and seismic reflection results in the AT1 region,41-42 the methane hydrate concentrated zone (MHCZ) is approximately 277–318 meters below seafloor (mbsf). The MHCZ comprises two zones: alternations of sand and mud layers (upper-depth zone), and sand-dominant sequences (lower-depth zone). Further details on the “Daini-Atusmi Knoll” drilling campaign have been described in the literature.40, 42-43 The sample used in this study was recovered from 286.65–286.77 mbsf in AT1-C. We called the sample “AT1-C-10P, ” following the “Daini-Atusmi Knoll” drilling campaign.42 In the DainiAtusmi Knoll area, instead of nodule- and massive-type hydrates, fine hydrate particles exist in the pore spaces of sandy sediments, the so-called pore-filling type.44 AT1-C-10P (286.65–286.77 mbsf) was recovered from a sand layer in an alternation of sand and mud layers. The grain size of AT1-C-10P was 84.5 µm (median). The index properties of the AT1-C-10P core are listed in Table 1. The grain size distribution of AT1-C-10P is shown in Figure S1 of the Supporting Information. Pressure-core Nondestructive Analysis Tools (PNATs). Our developed Pressure-core Nondestructive Analysis Tools (PNATs) were designed for analyzing hydrostatic pressurized


6

Page 7 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

natural GH sediments. The PNATs essentially comprised a manipulator for core-cutting and transferring (PNATs-MAN) and the IPTC (PNATs-AIST-IPTC) systems. The PNATs-MAN and PNATs-AIST-IPTC were installed in a cold room kept at a controlled temperature of 278–283 K. The basic PNATs components were developed in collaboration with Georgia Tech and USGS; these components and the original IPTC/PCCTs are described in the literature.17-18, 25-26 We supplemented the already-developed PNATs components by developing a set of nondestructive analysis tools for pressurized natural GH sediments; these are described below and their functions are summarized in Table 2. Nondestructive Imaging of Sediment Structure (PNATs-X). The sediment structure of the pressurized natural GH sediment sample was nondestructively measured using an X-ray CT system (PNATs-X, TOSCANER-32300µFD, Toshiba IT control systems corporation). The PNATs-X was equipped with a flat panel detector (1024 × 1024 pixels) and a microfocus X-ray tube (maximum source voltage and current: 230 kV and 1000 µA, respectively; minimum focus size: 4 µm; white-beam X-rays). To obtain CT images of the pressurized natural GH sediment, the samples were stored in a high-pressure vessel made of duralumin (Al alloy containing Cu, Mn, and Mg), using PNATs-MAN. The vessel was contained in a cold box in which the temperature was maintained at approximately 276 K by circulating cold air during CT collection. The inner pressure of the vessel was maintained at 10 MPa using a syringe pump. Collection of the raw CT data of the sample was performed every 0.15° (2400 views) at 180 kV and 600 µA. A schematic of PNATs-X is shown in Figure S2 of the Supporting Information. P-wave Velocity Profiles (PNATs-AIST-IPTC-PWV). P-wave profiles of the sediment sample were measured using a contact pulse transmission system with piezoelectric PZT ceramic sensors (OYO Co., Japan) through probe ports in the PNATs-AIST-IPTC. A schematic of the P-wave


7

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 33

measurement is shown in Figure S3 of the Supporting Information. Details of the P-wave measurements via the IPTC system are described in the literature.25-26 The PWV can be estimated from the arrival time of a P-wave transmitted through the object from an oscillator to a receiver. We estimated the PWV, VP (m/s), of the core sample using the following equation:

V = L⁄(t − t )

(1)

where L, tArrival, and tDelay are the path length of the transmitting P-wave, the arrival time of the transmitted P-wave, and the signal delay in the measurement system, respectively. L was estimated using the installed distance between the P-wave sensors. In our measurement system, tDelay was 1.8 µs. The P-wave profiles were collected from the core at intervals of 1 cm. In general, pressurized natural GH sediments were stored in a plastic liner that protects the recovered sediment (Figure S3); therefore, the core sample was cut to a designated length and transferred from the original liner to a dedicated IPTC liner with holes every 1 cm using the PNATs-MAN after the PNATs–X measurement. The P-wave measurement probes were brought into direct contact with the core surface. Infrared Spectral Profiles (PNATs-AIST-IPTC-IR). IR spectra of sediment samples were measured using a Fourier transform IR spectrometer (VIR-200, Jasco Co., Japan) with an indium antimonide (InSb) detector and an ATR-IR system with a probe custom-designed for the IPTC system. A schematic of the ATR-IR measurement process used by the PNATs-AIST-IPTC and an image of the edge of an ATR-IR probe are shown in Figure 1. The ATR-IR probe is a coneshaped diamond crystal (incidence angle θ = 45°) with a surface area of approximately 10 mm2. As done for the P-wave measurement, the ATR-IR probe was brought into direct contact with the core surfaces. The ATR−IR spectra of the samples were measured using 32 scans at the same


8

Page 9 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

core positions as the P-wave measurements and at a resolution of 4 cm−1 over an observed range of 2220–5000 cm−1. A background spectrum was collected at atmospheric pressure.

RESULTS AND DISCUSSION Sediment Structure. Figure 2 shows a cross-sectional X-ray CT image of AT1-C-10P (286.65– 286.77 mbsf). The spatial resolution in the obtained CT image is approximately 65 µm/pixel. The core length and diameter are approximately 12 and 5.36 cm, respectively. The gray scale in Figure 2 reflects the relative linear X-ray mass attenuation coefficients of each component, with the bright gray regions showing material with high X-ray attenuation; such regions represent the natural GH sediment. The low gray region outside the core shows the original plastic liner (0.3 cm thickness). The cross-sectional CT image shows that the core sample used in this study comprises three regions: the first (bottom) is 0–3 cm in core length (286.77–286.74 mbsf), the second (middle) is 3–6 cm (286.74–286.71 mbsf) in length, and the third (top) is 6–12 cm (286.71–286.65 mbsf) in length. This separation into three regions was created during the core recovery operation.42 Using the PNATs-MAN, the AT1-C-10P sediment was divided into 2 pieces (286.77–286.695 mbsf and 286.695–286.65 mbsf). In this study, one piece from 286.77– 286.695 mbsf (0–7.5 cm as indicated by the square in Figure 2) was used for measuring the Pwave profiles and ATR-IR spectra. Hereafter, the measured points are described using the core length (0–7.5 cm), as indicated in Figure 2. P-wave Velocity. Figure 3 shows the measured P-wave profiles of AT1-C-10P (286.77–286.695 mbsf) at 283 K and 10 MPa. The x- and y-axes of Figure 3 show the elapsed time from the Pwave oscillation starting point and the amplitude of the P-wave signal, respectively. Closed circles indicate the arrivals of the first P-wave signals through each measurement point. The


9

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 33

arrival times of the P-waves are delayed by the stiffness of the sample; in general, soft materials show longer delays than hard materials. The first P-waves arriving at 2 cm were much faster than those at other points, as shown in Figure 3. The actual arrival time at 2 cm can be estimated to be 16.5 µs by subtracting tDelay = 1.8 µs from tArrival = 18.2 µs. The actual arrival times at each point are listed in Table 3. The 0 cm in the cross-sectional CT image is considered to be the measurement point for the P-wave of water, as shown in Figure 2. In the P-wave profile at 1 cm, P-waves with high amplitude are observed at 35 µs after the first P-wave arrivals at 27.3 µs. These high-amplitude waves are similar to the P-wave profile of water at 0 cm. There are no such waves at 2–7 cm. The high-amplitude wave at 1 cm is considered to be a surface wave that travels near the sample surface; at this position, there may have been a slight poor contact between the PWV probe and the core sample. In the case of poor probe contact, a surface wave traveling along the water layer near the core surface would be observed in addition to the P-wave throughout the sample. The PWVs estimated using eq. (1) are also summarized in Table 3. The estimated PWV at 0 cm is very close to that of water (approximately 1480 m/s).33 The PWVs at 2–7 cm were distributed from approximately 2000 to 3150 m/s. The PWV value of the pure methane hydrate is approximately 3800 m/s,45 and the PWVs of sandy GH sediments depend on the volume of gas hydrates in the pores.46-47 Variation of the PWV would reveal differences in the Sh at each point. Because the PWV at 2 cm has the highest value at approximately 3150 m/s among those in the 1–7 cm region, the sample region at around 2 cm is considered to contain the heaviest hydrate volume in the pores, i.e., the highest Sh. Based on the previous literature,36 the values of Sh at 2 cm can be predicted to be approximately 30%, 50%, and 95% in the cementing, grain-coating, and load-bearing cases, respectively.


10

Page 11 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

ATR-IR Spectra. Figure 4 shows the ATR-IR spectra in AT1-C-10P (286.77–286.695 mbsf) at 283 K and 10 MPa. Considering the ATR window area and the grain size in AT1-C-10P, the ATR-IR spectra in Figure 4 would show the average state of H2O molecules over a large amount of pore spaces. For spectral comparison, the ATR-IR spectra of liquid water (283 K and 10 MPa) and ice (atmospheric conditions) are also displayed as blue and light blue lines, respectively, in Figure 3. To measure the ATR-IR spectra of fresh sample surfaces, prior to spectral measurements, we bored a hole (of a depth of approximately 5 mm) through the silt layers covering the core surfaces using a drill included in the PNATs-AIST-IPTC. The ATR-IR bands in the range 2800–3600 cm−1 are designated as the O–H stretching vibration bands of H2O molecules; a change in the shape of these can be seen in Figure 4. In particular, the band shape at 2 cm differs from those at other points. Furthermore, the C–H stretching peak of the CH3-group (at approximately 2910 cm−1) can be observed at 2 cm. A weak C–H stretching peak was also detected at 5 cm. The peak observed at around 2910 cm−1 can be assigned to CH4 molecules encaged in the gas hydrate structure. The O–H stretching vibration bands of H2O molecules comprise hydrogen bonded (HBd)/nonHBd symmetric and HBd/non-HBd asymmetric O–H stretching bands. It is generally difficult to decompose the complex O–H stretching band to liquid and solid H2O bands in order to directly estimate the ratio of liquid/solid H2O in the mixture. However, an O–H stretching vibration band observed in a liquid–solid H2O mixture shows an absorption band that is a composite of liquid and solid H2O bands. The shape of the composite absorption band depends on the ratio of liquid– solid H2O areas contacted by the ATR-IR window. As in the ice state, H2O molecules in a hydrate state are bonded to four neighboring H2O molecules by four hydrogen bonds. Therefore, the shapes of the O–H stretching bands of GH sediments depend on the water–hydrate ratio in


11

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 33

the pore spaces, i.e., they depend on the hydrate saturation.38-39 We estimated the hydrate saturation by evaluating the band shapes of the ATR-IR spectra, which were evaluated as a first moment (gravity center) of the observed O–H stretching band. The first moment is obtained using

̅ = ( /)⁄ ( /)

(2)

where v, I0, and I are the wavenumber, the intensities of the blank IR spectrum, and the sample IR spectrum, respectively. The first moments were estimated in the wavenumber range 2500– 4000 cm−1. The first moments of water and ice can be estimated as approximately 3257.9 and 3275.6 cm−1, respectively, and a high saturation point is expected to show a high first moment value. Using the first moments of liquid (water) and solid (ice) H2O molecules, we could then establish a correlation between the first moment and the hydrate saturation via ATR-IR, ShIR:

ℎ = 0.056941 × ̅ − 185.5

(3)

Equation (3) is only valid for our ATR-IR system; systems with different combinations of ATR windows and detectors will show different correlations between the first moment and the Sh.39 The first moment and the ShIR values at each point were estimated using eqs. (2) and (3) and are listed in Table 3. The value of ShIR was widely distributed from 30 to over 90% in the sample, except at 0 cm. It is clear that the measurement points at 2 and 5 cm have high ShIR (or first moment) values. The intensity of the C–H stretching peaks seems to correlate with the ShIR value, as shown in Figure 4. Because the measurement point at 0 cm, a nearly completely water region as described above, the value of ShIR at 0 cm is reasonably close to 0% (estimated value: 2%). Considering the ShIR value at 0 cm, the uncertainty of the ShIR estimation is expected to be ±5 %.


12

Page 13 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Hydrate Morphology in Nankai GH Sediment. Figure 5 shows the profiles of the PWV, ShIR (derived from the first moment), the density, and the porosity in AT1-C-10P (286.77–286.695 mbsf). Here, the density data were measured using a gamma-ray density meter (PNATs-PGGamma, PM-1000, Nanogray Inc., Japan). Details of the PNATs-PG-Gamma system are described in the Supporting Information. The PWV and the ShIR were obtained by line (plane) and spot (surface) measurements, respectively. As shown in Figures 5(a) and 5(b), there is a strong correlation between the PWV and the ShIR values. The mechanical behavior of the pressurized AT1-C core sample can be well understood by assuming a homogeneous hydrate distribution in the planes perpendicular to the depth direction (hereafter the horizontal planes).18 The AT1-C-10P (286.77–286.695 mbsf) sample is considered to be relatively homogeneous in the horizontal planes (PWV measurement planes), while the PWV and ShIR profiles vary along the depth direction. The bulk density appears to be approximately 1.93–2.01 g/cm3, except at 0 cm, as shown in Figure 5(c). In the AT1-C core samples, the bulk densities of cores taken from the sand and mud layer alternation zones are in the range 1.8–2.1 g/cm3.42 The porosity distribution in the core sample was estimated as shown in Figure 5(d). The porosity, Φ, can be estimated using the relation

* = *+ × (1 − ,) + (*. × ℎ/0 + *1 × (1 – ℎ/0 )) × ,

(4)

where ρ, ρs, ρh, and ρw are the bulk, the grain (2.65 g/cm3), the hydrate (0.911 g/cm3),2 and liquid water densities (1.005 g/cm3),48 respectively. The porosity varies in the range 38–43% and the average porosity is ~40%. The estimated porosity is in close agreement with the porosities of other pressurized AT1-C cores taken from sand and mud layer alternating zones, which are in the


13

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 33

range 35–46%.18, 36, 49 Unlike the PWV and ShIR variations, the porosity is distributed through a narrow range in the core sample. It is unclear what causes the relationship between the PWV– density (porosity) and the ShIR– density (porosity) in the sample. Figure 6 shows the relationship between the PWV and the ShIR in AT1-C-10P (286.77–286.695 mbsf). Other pressurized AT1-C sample results are also plotted in Figure 6.18, 36, 49 It is clear that the PWV increases with the Sh value. The PWV–ShIR relationship in this study is consistent with that found in the literature for the ShIR range up to 80%. Konno et al.36 calculated the PWV–Sh correlation for the three hydrate morphology models (cementing, grain coating, and loadbearing) in AT1-C sediments. The predicted PWV–Sh curves are shown as dashed lines in Figure 6. The PWV–ShIR correlation, including the highest ShIR point, is in close agreement with the predicted load-bearing curve. If the Sh in a load-bearing GH sediment approaches 100%, the PWV approaches values corresponding to the cementing and grain coating models. The PWV at the highest ShIR point is approximately 3150 m/s, and tentative P-wave data from CHIKYU have shown a value of approximately 3000 m/s near the same position.42 Because there have been few reports of sandy natural core samples with high Sh (>90%), this maximum ShIR value might be slightly over-estimated, and it will be necessary to obtain a more accurate ShIR estimate in order to improve eq. (3). Sediments with a low Sh value in a load-bearing type may present insignificant inhomogeneity of hydrate distribution in the pores (patchy type). In this case, a slight gap between the prediction and the experimental data for a ShIR value of 30–50% would be observed owing to the inhomogeneity of the Sh value. The properties of the pressurized natural GH sediment can be deduced using a combination of our PNATs system. The Sh values are widely distributed in the range 30% to over 90% in the


14

Page 15 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

AT1-C-10P (286.77–286.695 mbsf), as shown in Figure 5(b). Furthermore, our analysis reveals that the pressurized AT1-C-10P sample corresponds to a load-bearing GH sediment. CONCLUSIONS We report an in situ hydrate morphology estimation of a natural gas hydrate-bearing (GH) sediment recovered from the eastern Nankai trough area, with the first gas production test from the hydrate layer carried out on March 12, 2013. As an addition to our Pressured Non-destructive Analysis Tools (PNATs) that can measure the properties of GH sediments without releasing the sample pressure, we developed an X-ray computed tomography (CT) system (PNATs-X). PNATs-X can deduce sediment structural profiles in natural gas hydrate-bearing (GH) sediment cores (maximum length: ~1.2 m). We also developed an attenuated total reflection IR (ATR-IR) probe for the Instrumented Pressure Testing Chamber (PNATs-AIST-IPTC). The hydrate saturation, i.e., the Sh value at a measured point, can be evaluated using ATR-IR measurements through PNATs-AIST-IPTC. In addition to ATR-IR measurements, we measured the P-wave velocity of the natural GH sediment using the PNATs-AIST-IPTC system. We found a strong correlation between the Sh estimated by the ATR-IR band and the PWV, in which points with higher PWVs showed higher Sh values. This correlation is in close agreement with the literature data. Assuming that hydrates are homogeneously distributed in planes perpendicular to the sample depth direction, the hydrate morphology of natural GH sediments in the eastern Nankai trough area can be predicted to be of a load-bearing type. The combination of PNATs including ATR-IR spectroscopy could be used to nondestructively evaluate important properties of gas hydrate-bearing sediments in order to model gas hydrate reservoirs and predict their gas production behavior.


15

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 33

ACKNOWLEDGMENTS This study was supported by funding from the Research Consortium for Methane Hydrate Resources in Japan (MH21 Research Consortium), planned by the Ministry of Economy, Trade, and Industry (METI), Japan. The authors thank Dr. K Egawa of INPEX Corporation, Dr. T. Ito of Research institute of Innovative Technology for the Earth (RITE), Dr. H. Haneda, Mr. K. Shinjo, Mr. T. Uchiumi, Me. H. Kaneko of AIST for their experimental support. The authors also thank all the members of the shipboard team of the drilling campaign “Daini-Atusmi Knoll” in June–July 2012.

Supporting Information.

The grain-size distribution of the AT1-10P sample used in this

paper, schematics of PNATs-X (X-ray CT system), the PWV measurement system through PNATs-AIST-IPTC, and the PNATs-PG-Gamma (a gamma density measurement system).

Funding Sources.

The Research Consortium for Methane Hydrate Resources in Japan

(MH21 Research Consortium)


16

Page 17 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

REFERENCES 1.

Franks, F., Water: A Comprehensive Treatise. Plenum Press: London, 1973; Vol. 2.

2.

Sloan, E. D.; Koh, C. A., Clathrate Hydrates of Natural Gasses. 3rd ed.; CRC Press:

New York, 2007. 3.

Konno, Y.; Jin, Y.; Shinjou, K.; Nagao, J., Experimental evaluation of the gas recovery

factor of methane hydrate in sandy sediment. Rsc Adv 2014, 4 (93), 51666-51675. 4.

Wang, Y.; Li, X. S.; Li, G.; Huang, N. S.; Feng, J. C., Experimental study on the hydrate

dissociation in porous media by five-spot thermal huff and puff method. Fuel 2014, 117, 688696. 5.

Wang, Y.; Feng, J. C.; Li, X. S.; Zhang, Y.; Li, G., Analytic modeling and large-scale

experimental study of mass and heat transfer during hydrate dissociation in sediment with different dissociation methods. Energy 2015, 90, 1931-1948. 6.

Wang, Y.; Feng, J. C.; Li, X. S.; Zhang, Y.; Li, G., Large scale experimental evaluation

to methane hydrate dissociation below quadruple point in sandy sediment. Applied Energy 2016, 162, 372-381. 7.

Makogon, Y. F.; Omelchenko, R. Y., Commercial gas production from Messoyakha

deposit in hydrate conditions. Journal of Natural Gas Science and Engineering 2013, 11, 1-6. 8.

Makogon, Y. F., Natural gas hydrates – A promising source of energy. Journal of

Natural Gas Science and Engineering 2010, 2 (1), 49-59. 9.

Dallimore, S. R.; Yamamoto, K.; Wright, J. F.; G., B., Scientific Results from the

JOGMEC/NCRCan/Aurora Mallik 2007-2008 Gas Hydrate Production Research Well Program, Mackenzie Delta, Northwest Territories, Canada. Natural Resources Canada, Earth Sciences Sector, Geological Survey of Canada: 2012; Vol. 601.


17

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

10.

Page 18 of 33

Dallimore, S. R.; Collett, T. S.; Taylor, A. E., Scientific results from the Mallik 2002 gas

hydrate production research well program, Mackenzie Delta, Northwest Territories, Canada. Geological Survey of Canada Ottawa, Ontario, Canada: 2005. 11.

The National Energy Technology Laboratory (NETL), Gas Hydrate Production Trial

Using CO2 / CH4 Exchange. http://www.netl.doe.gov/research/oil-and-gas/projectsummaries/methane-hydrate/de-nt0006553. 12.

Japan Oil, Gas, and Metals National Corporation (JOGMEC) Gas Produced from

Methane Hydrate. http://www.jogmec.go.jp/english/news/release/content/300101080.pdf. 13.

Dickens, G. R.; Paull, C. K.; Wallace, P., Direct measurement of in situ methane

quantities in a large gas-hydrate reservoir. Nature 1997, 385, 426-428. 14.

Kawasaki, M.; Umezu, S.; Yasuda, M., Pressure Temperature Core Sampler (PTCS). J.

Jpn. Assoc. Petrol. Technol. 2006, 71 (1), 139-147. 15.

Ruppel, C.; Boswell, R.; Jones, E., Scientific results from Gulf of Mexico gas hydrates

Joint Industry Project Leg 1 drilling: introduction and overview. Marine and Petroleum Geology 2008, 25 (9), 819-829. 16.

Kubo, Y.; Mizuguchi, Y.; Inagaki, F.; Yamamoto, K., A new hybrid pressure-coring

system for the drilling vessel Chikyu. Scientific Drilling 2014, 17, 37-43. 17.

Jin, Y.; Konno, Y.; Nagao, J., Pressurized subsampling system for pressured gas-hydrate-

bearing sediment: Microscale imaging using X-ray computed tomography. Review of Scientific Instruments 2014, 85 (9). 18.

Yoneda, J.; Masui, A.; Konno, Y.; Jin, Y.; Egawa, K.; Kida, M.; Ito, T.; Nagao, J.;

Tenma, N., Mechanical behavior of hydrate-bearing pressure-core sediments visualized under triaxial compression. Marine and Petroleum Geology 2015, 66 (2), 451-459.


18

Page 19 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

19.

Konno, Y.; Yoneda, J.; Egawa, K.; Ito, T.; Jin, Y.; Kida, M.; Suzuki, K.; Fujii, T.; Nagao,

J., Permeability of sediment cores from methane hydrate deposit in the Eastern Nankai Trough. Marine and Petroleum Geology 2015, 66, 487-495. 20.

Kneafsey, T. J.; Lu, H. L.; Winters, W.; Boswell, R.; Hunter, R.; Collett, T. S.,

Examination of core samples from the Mount Elbert Gas Hydrate Stratigraphic Test Well, Alaska North Slope: Effects of retrieval and preservation. Marine and Petroleum Geology 2011, 28 (2), 381-393. 21.

Jin, Y.; Hayashi, J.; Nagao, J.; Suzuki, K.; Minagawa, H.; Ebinuma, T.; Narita, H., New

method of assessing absolute permeability of natural methane hydrate sediments by microfocus X-ray computed tomography. Jpn J Appl Phys 1 2007, 46 (5a), 3159-3162. 22.

Oyama, H.; Konno, Y.; Suzuki, K.; Nagao, J., Depressurized dissociation of methane-

hydrate-bearing natural cores with low permeability. Chemical Engineering Science 2012, 68 (1), 595-605. 23.

Konno, Y.; Oyama, H.; Nagao, J.; Masuda, Y.; Kurihara, M., Numerical Analysis of the

Dissociation Experiment of Naturally Occurring Gas Hydrate in Sediment Cores Obtained at the Eastern Nankai Trough, Japan. Energy & Fuels 2010, 24 (12), 6353-6358. 24.

Schultheiss, P.; Holland, M.; Humphrey, G., Wireline Coring and Analysis under

Pressure: Recent Use and Future Developments of the HYACINTH System. Scientific Drilling 2009, (7, March 2009). 25.

Yun, T. S.; Narsilio, G. A.; Santamarina, J. C.; Ruppel, C., Instrumented pressure testing

chamber for characterizing sediment cores recovered at in situ hydrostatic pressure. Marine Geology 2006, 229 (3-4), 285-293.


19

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

26.

Page 20 of 33

Santamarina, J. C.; Dai, S.; Jang, J.; Terzariol, M., Pressure Core Characterization Tools

for Hydrate-Bearing Sediments. Scientific Drilling 2012, 14, 44-48. 27.

Yun, T. S.; Lee, C.; Lee, J.-S.; Bahk, J. J.; Santamarina, J. C., A pressure core based

characterization of hydrate-bearing sediments in the Ulleung Basin, Sea of Japan (East Sea). Journal of Geophysical Research 2011, 116 (B2). 28.

Yun, T. S.; Fratta, D.; Santamarina, J. C., Hydrate-Bearing Sediments from the Krishna−

Godavari Basin: Physical Characterization, Pressure Core Testing, and Scaled Production Monitoring. Energy & Fuels 2010, 24 (11), 5972-5983. 29.

The Research Consortium for Methane Hydrate Resources in Japan (MH21).

http://www.mh21japan.gr.jp/english/. 30.

Waite, W. F.; Santamarina, J. C.; Cortes, D. D.; Dugan, B.; Espinoza, D. N.; Germaine,

J.; Jang, J.; Jung, J. W.; Kneafsey, T. J.; Shin, H.; Soga, K.; Winters, W. J.; Yun, T. S., Physical properties of hydrate-bearing sediments. Reviews of Geophysics 2009, 47 (4). 31.

Yun, T. S.; Francisca, F. M.; Santamarina, J. C.; Ruppel, C., Compressional and shear

wave velocities in uncemented sediment containing gas hydrate. Geophysical Research Letters 2005, 32 (10). 32.

Berge, L. I.; Jacobsen, K. A.; Solstad, A., Measured acoustic wave velocities of R11

(CCl3F) hydrate samples with and without sand as a function of hydrate concentration. J Geophys Res-Sol Ea 1999, 104 (B7), 15415-15424. 33.

Lee, M. W.; Hutchinson, D. R.; Collett, T. S.; Dillon, W. P., Seismic velocities for

hydrate-bearing sediments using weighted equation. J Geophys Res-Sol Ea 1996, 101 (B9), 20347-20358.


20

Page 21 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

34.

Priest, J. A.; Best, A. I.; Clayton, C. R. I., A laboratory investigation into the seismic

velocities of methane gas hydrate-bearing sand. J Geophys Res-Sol Ea 2005, 110 (B4). 35.

Dovrkin, J.; Helgerud, M. B.; Waite, W. F.; Kirby, S. H.; Nur, A., Introduction to

physical properties and elasticity models. Natural Gas Hydrate In Oceanic and Permafrost Environments, Kluwer Academic Publishers: Dordrecht, Netherlands, 2000; p 245-260. 36.

Konno, Y.; Jin, Y.; Yoneda, J.; Kida, M.; Egawa, K.; Ito, T.; Suzuki, K.; Nagao, J., Effect

of methane hydrate morphology on compressional wave velocity of sandy sediments: Analysis of pressure cores obtained in the Eastern Nankai Trough. Marine and Petroleum Geology 2015, 66, 425-433. 37.

Kida, M.; Jin, Y.; Watanabe, M.; Konno, Y.; Yoneda, J.; Egawa, K.; Ito, T.; Nakatsuka,

Y.; Suzuki, K.; Fujii, T.; Nagao, J., Chemical and crystallographic characterizations of natural gas hydrates recovered from a production test site in the eastern Nankai Trough. Marine and Petroleum Geology 2015, 66, 396-403. 38.

Jin, Y.; Oyama, H.; Nagao, J., Infrared Spectroscopy of Gas Hydrate Dissociation

Behavior During Depressurization. Japanese Journal of Applied Physics 2009, 48 (10), 108001. 39.

Jin, Y.; Konno, Y.; Nagao, J., Growth of Methane Clathrate Hydrates in Porous Media.

Energy & Fuels 2012, 26 (4), 2242-2247. 40.

Yamamoto, K., Overview and introduction: Pressure core-sampling and analyses in the

2012e2013 MH21 offshore test of gas production from methane hydrates in the eastern Nankai Trough. Marine and Petroleum Geology 2015, 66 (2), 296-309. 41.

Fujii, T.; Suzuki, K.; Takayama, T.; Tamaki, M.; Komatsu, Y.; Konno, Y.; Yoneda, J.;

Yamamoto, K.; Nagao, J., Geological setting and characterization of a methane hydrate reservoir


21

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 33

distributed at the first offshore production test site on the Daini-Atsumi Knoll in the eastern Nankai Trough, Japan. Marine and Petroleum Geology 2015, 66, 310-322. 42.

Suzuki, K.; Schultheiss, P.; Nakatsuka, Y.; Ito, T.; Egawa, K.; Holland, M.; Yamamoto,

K., Physical properties and sedimentological features of hydrate-bearing samples recovered from the first gas hydrate production test site on Daini-Atsumi Knoll around eastern Nankai Trough. Marine and Petroleum Geology 2015, 66 (2), 346-357. 43.

Inada, N.; Yamamoto, K., Data report: Hybrid Pressure Coring System tool review and

summary of recovery result from gas-hydrate related coring in the Nankai Project. Marine and Petroleum Geology 2015, 66 (2), 323-345. 44.

Kida, M.; Suzuki, K.; Kawamura, T.; Oyama, H.; Nagao, J.; Ebinuma, T.; Narita, H.;

Suzuki, H.; Sakagami, H.; Takahashi, N., Characteristics of Natural Gas Hydrates Occurring in Pore-Spaces of Marine Sediments Collected from the Eastern Nankai Trough, off Japan. Energy & Fuels 2009, 23, 5580-5586. 45.

Helgerud, M. B.; Waite, W. F.; Kirby, S. H.; Nur, A., Elastic wave speeds and moduli in

polycrystalline ice Ih, sI methane hydrate, and sII methane-ethane hydrate. J Geophys Res-Sol Ea 2009, 114. 46.

Waite, W. F.; Winters, W. J.; Mason, D. H., Methane hydrate formation in partially

water-saturated Ottawa sand. American Mineralogist 2004, 89, 1202-1207. 47.

Choi, J. H.; Dai, S.; Cha, J. H.; Seol, Y., Laboratory formation of noncementing hydrates

in sandy sediments. Geochemistry, Geophysics, Geosystems 2014, 15 (4), 1648-1656. 48.

Lemmon, E. W.; McLinden, M. O.; Friend, D. G., "Thermophysical Properties of Fluid

Systems" in NIST Chemistry WebBook, NIST Standard Reference Database Number 69, Eds.


22

Page 23 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

P.J. Linstrom and W.G. Mallard, National Institute of Standards and Technology, Gaithersburg MD, 20899, http://webbook.nist.gov, (retrieved February 29, 2016). ed. 49.

Santamarina, J. C.; Dai, S.; Terzariol, M.; Jang, J.; Waite, W. F.; Winters, W. J.; Nagao,

J.; Yoneda, J.; Konno, Y.; Fujii, T.; Suzuki, K., Hydro-bio-geomechanical properties of hydratebearing sediments from Nankai Trough. Marine and Petroleum Geology 2015, 66, 434-450.


23

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 33

FIGURE CAPTIONS Figure 1. ATR-IR probe developed for IPTC. (a) Schematic of the measurement using the ATRIR probe. (b) Image of the edge of the ATR-IR probe. Figure 2. Cross-sectional X-ray CT image of AT1-C-10P (286.65–286.77 mbsf). The core length and diameter are 12 cm and 5.36 cm, respectively. The core top and bottom are 286.64 and 286.77 mbsf in depth, respectively. The core sample used for the PWV and ATR-IR measurements is shown within the square. Figure 3. P-wave profiles for AT1-C-10P (286.77–286.695 mbsf). The closed circles at each measurement point indicate the determined arrival times. P-wave measurements were carried out at 278 K and 10 MPa. Figure 4. O–H stretching ATR-IR bands in AT1-C-10P (286.77–286.70 mbsf). ATR-IR spectroscopic measurements were carried out at 283 K and 10 MPa. Figure 5. Property profiles for AT1-C-10P (286.77–286.695 mbsf). (a) P-wave velocity. (b) Hydrate saturation estimated from the ATR-IR spectra. (c) Density obtained using a gamma-ray density system. (d) Porosity estimated from the hydrate saturation and density. Each property at 0 cm was obtained from liquid water, not from the sediment. Figure 6. Correlation between P-wave velocity and hydrate saturation in AT1-C-10P (286.77– 286.695 mbsf). Closed circles: this study; open circles: Konno et al.;36 open squares: Santamarina et al.;49 open diamond: Yoneda et al.18 The dashed lines are PWV–Sh correlations predicted from the properties of the AT1-C samples and the hydrate morphology types (red: cementing; black: grain-coating; blue: load-bearing).36 The PWV–Sh correlation at 0 cm is not shown.


24

Page 25 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

FIGURES

Figure 1.


25

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 33

Figure 2.


26

Page 27 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Figure 3.


27

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 33

Figure 4.


28

Page 29 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Figure 5.


29

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 30 of 33

Figure 6.


30

Page 31 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

TABLES Table 1. Indexes of the natural GH sediment sample AT-C-10P used in this study

Median

Depth (mbsf) Core length

grain size top

bottom

(cm)

286.65

286.77

12

Classification

(µm) 84.5

Silty sand


31

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 32 of 33

Table 2. Function summary of the developed PNAT family

PNATs family (PNATs-) MAN

X AIST-IPTC

-PWV -IR

PG

-Gamma

-PWV TACCT SUB

Descriptions Manipulating sediment samples for sample cutting in a designated length, sample transfer to high-pressure vessel/PNATs apparatus, and sample moving to change point of measurement by using PNATs apparatus. Micro-focus X-ray CT system for measuring 3dimensional sediment structure. Core analysis system developed with the strong cooperation of Georgia Tech and the USGS. Core properties can be measured using several measurement probes. Use with PNATs-MAN. Probe designed for IPTC. P-wave velocity measurements for estimating hydrate saturation. AIST original probe designed for IPTC. ATRIR spectra measurement for estimating hydrate saturation. Gamma-ray density system for measuring bulk density of sediment samples. Use with PNATsMAN

Note Ref. 18.

New development. Basic concept is described in Ref. 25 New development. New development.

New development. Details in the Supporting Information. Non-contact P-wave velocity measurement Under developing. sytem. Use with PNATs-MAN Tri-axial testing system for measuring strength, Ref. 18 deformation behavior, and/or permeability. Re-sampling of small core from sediment Ref. 17 samples.


32

Page 33 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Table 3. P-wave and ATR-IR measurement results for AT1-C-10P (286.77–286.695 mbsf) Core

Actual

Wave path

Depth

First

Hydrate

moment

saturation

(cm-1)**

(%)

P-wave position

arrival time

length

(mbsf)

velocity (m/s) (cm)

(µs)*

(mm)

286.77

0

30.0

44.6

1486.7

3258.1

(2.0)

286.76

1

25.5

49.6

1945.1

3263.6

33.3

286.75

2

16.4

51.6

3146.3

3275.0

98.2

286.74

3

20.6

49.6

2407.8

3266.8

51.5

286.73

4

24.6

52.1

2117.9

3263.7

33.8

286.72

5

21.0

53.6

2481.0

3270.4

72.0

286.71

6

23.2

52.6

2267.2

3262.8

28.7

286.70

7

24.0

52.6

2191.7

3266.5

49.8

*Actual arrival time was estimated by subtracting the system delay (1.8 µs) **The first moment values of water and ice Ih are 3257.9 and 3275.6 cm−1, respectively.


33

In Situ Methane Hydrate Morphology Investigation: Natural Gas

Recommend Documents