The Acoustic Properties of Sandy and Clayey Hydrate-Bearing Sediments

Natural gas hydrates samples are rare and difficult to store and transport at in situ pressure and temperature conditions, resulting in difficulty to characterize natural hydrate-bearing sediments and to identify hydrate accumulation position and saturation at the field scale. A new apparatus was designed to study the acoustic properties of seafloor recovered cores with and without hydrate. To protect the natural frames of recovered cores and control hydrate distribution, the addition of water into cores was performed by injecting water vapor. The results show that hydrate saturation and types of host sediments are the two most important factors that govern the elastic properties of hydrate-bearing sediments. When gas hydrate saturation adds approximately to 5–25%, the corresponding P-wave velocity (Vp) increases from 1.94 to 3.93 km/s and S-wave velocity (Vs) increases from 1.14 to 2.23 km/s for sandy specimens; Vp and Vs for clayey samples are 1.72–2.13 km/s and 1.10–1.32 km/s, respectively. The acoustic properties of sandy sediments can be significantly changed by the formation/dissociation of gas hydrate, while these only minorly change for clayey specimens.


Introduction
The combustion of traditional fossil fuel has caused serious environmental issues, such as global warming, leading to the growing demand for low-carbon energy.Natural gas hydrates are solid crystals, which is a special form of natural gas that is widely accumulated in subsea and permafrost zones [1].It has been estimated that the gross organic carbon in the form of natural gas hydrates exceeds the reserves of all known fossil fuels on the earth, including coal, oil, and conventional gases [2,3].More importantly, the mole fraction of methane in natural gas hydrates is up to 99%, meaning that it can supply clean energy [4,5].To exploit this unconventional energy, it is necessary to accurately predict the location of hydrate reservoirs and to estimate the hydrate saturation and the amount of gas trapped in hydrates, which requires the development of detection methods.
Compared with other pore-filling fluids, hydrate-containing sediments exhibit relatively high acoustic velocities, which can be used to identify natural gas hydrate deposits.Song et al. [6] performed a series of tri-axial compression tests to study the mechanical properties of sediments with and without hydrate, finding that the presence of methane hydrate can increase the shear strength of host sediments.Priest et al. [7] measured the seismic velocities of methane hydrate-bearing sediments that were formed by tetrahydrofuran (THF) solution and porous media.The addition of surfactants is used to promote hydrate nucleation to reduce the induction time.Their results showed that when hydrate saturation Energies 2019, 12, 1825 2 of 11 increases from 0 to 35.27%, the corresponding P-wave velocity increased from 1639 to 2833 m/s.Li et al. [8] reported V p values from 1706 to 1782 m/s for THF aqueous solution-saturated sandy sediment, which then reached 3295-3984 m/s after completion of hydrate formation.Although V p increases with hydrate saturation, Zhang et al. [9] could not establish a universal function to correlate V p with hydrate saturation.A Gas Hydrate And Sediment Test Laboratory Instrument (GHASTLI), designed by the U.S. Geological Survey, was used to study the bulk properties of gas hydrate-bearing sediments.Through this facility, Winters et al. [10] measured the V p of coarse-grained sediments with different pore space occupancy, such as free gas, ice, water, and gas hydrates.The results showed a V p range from less than 1 km/s for gas-charged sediment to 1.77-1.94km/s for water-saturated sediment and 2.91-4.00km/s for sediment with varying degrees of hydrate saturation.Waite et al. [11] reported that methane hydrate can cement unconsolidated sediment when forming in an abundant gas environment, which would increase the acoustic velocity of hydrate-bearing sediments.Hu et al. [12] suggested that gas hydrate can be treated as a component within a matrix of consolidated sediments when hydrate saturation is above 30% in pore space, otherwise it's just part of pore fluid or adheres to the surrounding rocks.For instance, they found that numerous small particles start to coalesce together when the hydrate saturation is above 14% for the sediments recovered from the South China Sea [13].Konno et al. [14] also measured the P-wave velocity of natural sediment cores from the Eastern Nankai Trough.The result demonstrates that V p varies from less than 1700 m/s in the hydrate-free section to greater than 2300 m/s in the section with hydrate.At present, the relationship between P-wave velocity and hydrate saturation has been widely investigated by using artificial samples.
Except for hydrate saturation and pore space occupants, the lithologic characteristics of host sediments, test conditions (burial depth) and hydrate distribution can also affect the acoustic velocity of hydrate-bearing sediments [15,16].The reported gas hydrate drilling expeditions found natural gas hydrates are widely abundant in sandy sediments, such as reservoirs in the Nankai trough of offshore Japan [17,18], the Mallik permafrost site in the Northern Canada [19], and clayey sediments at several sites in the Gulf of Mexico [20], the Shenhu area in the South China Sea [21], the Northern Cascadia Margin [22], and the Krishna-Godavari basin in India [23].However, research on the influence of lithologic characteristics of host sediments and hydrate distribution on the acoustic velocities of hydrate-bearing sediments is limited, which are important to understand the interaction between gas hydrates and host sediments.
Owing to the high cost of recovering natural gas hydrates specimens, laboratory-formed hydrate-bearing sediments are widely used to provide calibration data for seismic detection.For preparing man-made gas hydrate-bearing sediments, one of the challenges is accurately adding predetermined water into host sediments.At present, the most commonly used method is to slowly push gas into initially fully water-saturated sandy sediment until a predetermined amount of water is pushed out of the sample [10].The second method is the partial water saturation method that porous media are mixed with a limited amount of water to form partially water-saturated sediment [24,25].These methods don't take hydrate distribution into consideration, yielding different velocities, even at the same hydrate saturation.The third method is the ice-seeding method that mixes sediments with small ice grains [26,27].However, this method is not suitable for natural recovered cores, which may change the physical properties of cores, such as their porosity, permeability, and skeleton structure.To accommodate the need for better methods to obtain elastic properties, we propose a new water-adding method for the preparation of gas hydrate in natural recovered cores, and aim at studying the influence of lithologic characteristics of host sediments, test conditions, and pore space occupancy on the acoustic properties of hydrate-bearing sediments, which are necessary for geologic exploration and hydrate development.

Apparatus
Two types of apparatus are used in this study.One is used to prepare hydrate samples, which includes a water-feeding device, a vacuum pump, a drying oven, a fridge, and a high-pressure reactor.The other is used to measure ultrasonic wave velocity, which includes a pulse sender/receiver machine (5077PR, made by Olympus NDT, Waltham, MA, USA.), a digital oscilloscope (TBS2012B, made by Tektronix Inc., Shanghai, China), two transducers, a core holder setup, a hand pump, and a cryogenic house.The water-feeding device contains a bracket to support core, under which is deionized water, as shown in Figure 1.The temperature range of the dry oven is 0-200 • C, which is used to dry the cores and supply heat for water vaporization.The high-pressure reactor where hydrates were formed is constructed from stainless steel, with an effective volume of 0.1 L (Φ40 × 80) and a maximum working pressure of 60 MPa.The experimental apparatus for measuring ultrasonic wave velocity is shown in Figure 2. The size of cryogenic house is 5 m × 2.5 m × 3 m, with working temperature of 248-323 K.The core holder setup was placed into the cryogenic house to maintain the system temperature, thus protecting methane hydrate from decomposing during acoustic test.The confining pressure was controlled by a hand pump manufactured by Jiangsu Haian petroleum research Ltd. (Haian, China).In addition, a Monitor and Control Generated System (MCGS) manufactured by Beijing Kunlun Coast Sensing Technology Co. Ltd. (Beijing, China) was used to monitor and record the temperature and pressure per minute.

Apparatus
Two types of apparatus are used in this study.One is used to prepare hydrate samples, which includes a water-feeding device, a vacuum pump, a drying oven, a fridge, and a high-pressure reactor.The other is used to measure ultrasonic wave velocity, which includes a pulse sender/receiver machine (5077PR, made by Olympus NDT, Waltham, MA, USA.), a digital oscilloscope (TBS2012B, made by Tektronix Inc., Shanghai, China), two transducers, a core holder setup, a hand pump, and a cryogenic house.The water-feeding device contains a bracket to support core, under which is deionized water, as shown in Figure 1.The temperature range of the dry oven is 0-200 °C , which is used to dry the cores and supply heat for water vaporization.The high-pressure reactor where hydrates were formed is constructed from stainless steel, with an effective volume of 0.1 L (Ф40 × 80) and a maximum working pressure of 60 MPa.The experimental apparatus for measuring ultrasonic wave velocity is shown in Figure 2. The size of cryogenic house is 5 m × 2.5 m × 3 m, with working temperature of 248-323 K.The core holder setup was placed into the cryogenic house to maintain the system temperature, thus protecting methane hydrate from decomposing during acoustic test.The confining pressure was controlled by a hand pump manufactured by Jiangsu Haian petroleum research Ltd. (Haian, China).In addition, a Monitor and Control Generated System (MCGS) manufactured by Beijing Kunlun Coast Sensing Technology Co. Ltd. (Beijing, China) was used to monitor and record the temperature and pressure per minute.Energies 2019, 12, x 3 of 11

Apparatus
Two types of apparatus are used in this study.One is used to prepare hydrate samples, which includes a water-feeding device, a vacuum pump, a drying oven, a fridge, and a high-pressure reactor.The other is used to measure ultrasonic wave velocity, which includes a pulse sender/receiver machine (5077PR, made by Olympus NDT, Waltham, MA, USA.), a digital oscilloscope (TBS2012B, made by Tektronix Inc., Shanghai, China), two transducers, a core holder setup, a hand pump, and a cryogenic house.The water-feeding device contains a bracket to support core, under which is deionized water, as shown in Figure 1.The temperature range of the dry oven is 0-200 °C, which is used to dry the cores and supply heat for water vaporization.The high-pressure reactor where hydrates were formed is constructed from stainless steel, with an effective volume of 0.1 L (Ф40 × 80) and a maximum working pressure of 60 MPa.The experimental apparatus for measuring ultrasonic wave velocity is shown in Figure 2. The size of cryogenic house is 5 m × 2.5 m × 3 m, with working temperature of 248-323 K.The core holder setup was placed into the cryogenic house to maintain the system temperature, thus protecting methane hydrate from decomposing during acoustic test.The confining pressure was controlled by a hand pump manufactured by Jiangsu Haian petroleum research Ltd. (Haian, China).In addition, a Monitor and Control Generated System (MCGS) manufactured by Beijing Kunlun Coast Sensing Technology Co. Ltd. (Beijing, China) was used to monitor and record the temperature and pressure per minute.

Materials
Four natural cores recovered from seafloor were studied in this work, two clayey cores named C1 and C2, two sandy cores named S1 and S2.All of them were manufactured in cylindrical shape and protected by heat shrink tubing.The physical parameters of these cores are shown in Table 1.Methane gas with minimum purity of 99.9% was used to prepare hydrate, supplied by Beifen Gas Industry Corporation (Beijing, China).Distilled water was made in the lab.

Experimental Procedure
The whole experiment process consists of three parts, namely water addition, preparation of gas hydrate, and acoustic velocity measurement.In the first part, all the cores were dried in an oven for 24 h, and then transferred to the water-feeding vessel.After being vacuumed, the water-feeding vessel was placed into a dry oven and heated at 363.15 K, under which deionized water would vaporize and transfer into the pore space of cores.The water addition process took 3-6 h, depending on the porosity of core and predetermined water saturation.After the water vapor injection process, the cores were sent to a fridge where free water was converted into ice.In the second part, the core was placed into a high-pressure reactor in which water and methane gas combined to form hydrate. Owing to the induction time, the hydrate formation process would take at least 12 h.Subsequently, the hydrate-containing core was transferred from the high-pressure vessel to a core holder setup under the temperature of −6 • C. Owing to the self-preservation effect of gas hydrate, the dissociation rate of methane hydrate is lowest when the system temperature is around 268 K [28,29].After that, the pore pressure was raised up to 8 MPa by injecting methane gas and the confining pressure was set at 11 MPa in the core holder setup.As gas hydrate has a memory effect [30,31], meaning that hydrate forms faster when combined with the water from decomposed hydrate, another three hours was wasted to compensate for the possible dissociated hydrate during the core transfer process.A rubber sleeve was used to separate the core from the confining liquid.In the third part, we changed the pore pressure, confining pressure, and simultaneously collected P-wave and S-wave signals to investigate the effect of effective confining pressure on acoustic velocity.In addition, the exact volume of dissociated CH 4 gas from methane hydrate was measured by a drainage method.

Calculation Equations
The hydrate saturation was calculated based on following two equations.
where S h is the volume saturation of gas hydrate in pore space, m h is the mass of gas hydrate, g; ρ h is the density of methane hydrate, 0.910 g/cm 3 ; V pore is the total pore volume of specimen; m w is the mass of distilled water, g; n is the hydrate number, which is equal to 6 in this work [32].
Energies 2019, 12, 1825 5 of 11 The P and S-wave velocities can be calculated by the length of core and travel time of acoustic wave.
where L is the length of specimen, m; t is the travel time of P and S-wave, as shown in Figure 3; t 0 is the travel time when two detectors tightly touched together, the value is 5.704 × 10 −6 s.
the travel time when two detectors tightly touched together, the value is 5.704×10 −6 s.
The Vp/Vs ratio is calculated by Equation ( 4): Combined with Vp and Vs, the elastic modulus of specimen can be obtained by Equation (5).The effective confining pressure is the difference between confining pressure and pore pressure.

Acoustic Signal
Elastic modulus, shear modulus, and Poisson's ratio are the key parameters for developing the seismic inversion model, and can be figured out by Vp and Vs.Hence, it is important to acquire the Vp and Vs of hydrate-bearing sediments at the same time, rather than Vp only.On the basis of other seismic instruments [33], we built a new apparatus that not only collects P and S-wave signals, but also supplies confining pressure and pore pressure during the experiment process, allowing simulation of hydrate in situ accumulation conditions.One of the wave diagrams collected in the experiment is shown in Figure 3, from which the travel time can be achieved.Then, acoustic velocity and other parameters of hydrate-bearing sediments were calculated based on Equations (3)-(5), and are listed in Table 2.The V p /V s ratio is calculated by Equation ( 4): Combined with V p and V s , the elastic modulus of specimen can be obtained by Equation ( 5).
The effective confining pressure is the difference between confining pressure and pore pressure.
δ e f f = P C −P P (6)

Acoustic Signal
Elastic modulus, shear modulus, and Poisson's ratio are the key parameters for developing the seismic inversion model, and can be figured out by V p and V s .Hence, it is important to acquire the V p and V s of hydrate-bearing sediments at the same time, rather than V p only.On the basis of other seismic instruments [33], we built a new apparatus that not only collects P and S-wave signals, but also supplies confining pressure and pore pressure during the experiment process, allowing simulation of hydrate in situ accumulation conditions.One of the wave diagrams collected in the experiment is shown in Figure 3, from which the travel time can be achieved.Then, acoustic velocity and other parameters of hydrate-bearing sediments were calculated based on Equations (3)-(5), and are listed in Table 2.

Effect of Pore Occupants
The most relevant physical properties for predicting the behavior of gas hydrate-bearing sediments are pore space occupants, effective stress, and lithologic characteristics of host sediments.The influence of pore occupants on the elastic properties was investigated by measuring the acoustic velocities of recovered cores with water saturation of 20%, hydrate saturation of 20%, and a completely dry core.As shown in Figure 4, both V p and V s of hydrate-containing cores are much higher than that of water-containing and dry cores, indicating that gas hydrate is the most important pore occupant that could significantly change the elastic properties of host sediments.Gas hydrate not only can occupy the pore space, but also works as an adhesive to cement surrounding rock particles, which would further strengthen the skeleton structure of sediments [10,12].Hence, the increase of V p and V s are comprehensive results caused by the increase of hydrate concentration in the pore space.

Effect of Pore Occupants
The most relevant physical properties for predicting the behavior of gas hydrate-bearing sediments are pore space occupants, effective stress, and lithologic characteristics of host sediments.The influence of pore occupants on the elastic properties was investigated by measuring the acoustic velocities of recovered cores with water saturation of 20%, hydrate saturation of 20%, and a completely dry core.As shown in Figure 4, both Vp and Vs of hydrate-containing cores are much higher than that of water-containing and dry cores, indicating that gas hydrate is the most important pore occupant that could significantly change the elastic properties of host sediments.Gas hydrate not only can occupy the pore space, but also works as an adhesive to cement surrounding rock particles, which would further strengthen the skeleton structure of sediments [10,12].Hence, the increase of Vp and Vs are comprehensive results caused by the increase of hydrate concentration in the pore space.As well as pore occupants, the influence of effective confining pressure on acoustic velocities was also studied in this work.During the acoustic velocity test, the pore pressure was set at 8 MPa to maintain the stability of methane hydrate and the confining pressure ranged from 11 to 23 MPa, resulting in the maximum effective confining pressure reaches 15 MPa.As shown Figure 4, the acoustic Energies 2019, 12, 1825 7 of 11 velocity of hydrate-containing core linearly increases with effective confining pressure.When the effective confining pressure increases from 3 MPa to 15 MPa, the corresponding P and S-wave velocities would increase about 20% and 10%, respectively.On the other hand, when hydrate saturation changes from 0 to 20%, the corresponding acoustic velocity would increase by about 90%.This means that the effect of buried depth of natural gas hydrates on the seismic velocity can be neglected if the thickness of the hydrate layer is within 100 m.

Effects of Host Sediments
The host sediments of natural gas hydrates are mainly composed of gravel, sand, silt, clay, and mud [23].Here, we studied the acoustic properties of sandy and clayey hydrate-bearing sediments.If pore water content is below 30% in sandy sediments, most of the water is bounded in pore space [34].In addition, disseminated gas hydrate saturations in the pore space of fine-grained clayey sediments are typically <10% [23].Thereby, the volume fraction of methane hydrate distributed in the cores was controlled below 30% in this work.As shown in Figure 5, when hydrate saturation increases approximately from 5% to 25%, the corresponding V p ranges from 1.94 to 3.93 km/s, V s ranges from 1.14 to 2.23 km/s for sandy specimens; V p and V s for clayey specimens are 1.72-2.13km/s and 1.10-1.32km/s, respectively.When at the same hydrate saturation, the V p and V s of sandy specimens are higher than that of clayey specimens, and the difference enlarges with hydrate saturation.Winters et al. [10] suggested that the P-wave velocity of fine-grained sediment is higher than that of coarser-grained sediment, because their bury depths under the sea are different.However, the effective confining pressure was kept the same in our experiment, indicating the difference of acoustic velocity is mainly owing to the lithologic difference of host sediments and the interaction between hydrate sandy particles.
S-wave velocities would increase about 20% and 10%, respectively.On the other hand, when hydrate saturation changes from 0 to 20%, the corresponding acoustic velocity would increase by about 90%.This means that the effect of buried depth of natural gas hydrates on the seismic velocity can be neglected if the thickness of the hydrate layer is within 100 m.

Effects of Host Sediments
The host sediments of natural gas hydrates are mainly composed of gravel, sand, silt, clay, and mud [23].Here, we studied the acoustic properties of sandy and clayey hydrate-bearing sediments.If pore water content is below 30% in sandy sediments, most of the water is bounded in pore space [34].In addition, disseminated gas hydrate saturations in the pore space of fine-grained clayey sediments are typically <10% [23].Thereby, the volume fraction of methane hydrate distributed in the cores was controlled below 30% in this work.As shown in Figure 5, when hydrate saturation increases approximately from 5% to 25%, the corresponding Vp ranges from 1.94 to 3.93 km/s, Vs ranges from 1.14 to 2.23 km/s for sandy specimens; Vp and Vs for clayey specimens are 1.72-2.13km/s and 1.10-1.32km/s, respectively.When at the same hydrate saturation, the Vp and Vs of sandy specimens are higher than that of clayey specimens, and the difference enlarges with hydrate saturation.Winters et al. [10] suggested that the P-wave velocity of fine-grained sediment is higher than that of coarser-grained sediment, because their bury depths under the sea are different.However, the effective confining pressure was kept the same in our experiment, indicating the difference of acoustic velocity is mainly owing to the lithologic difference of host sediments and the interaction between hydrate and sandy particles.
The acoustic velocities measured in our experiment were also compared with the field data of Mallik 2L-38 well in Northern Canada [35], where the occurrence of natural gas hydrates has been verified.The gravel and sand content of natural samples obtained from Mallik 2L-38 borehole is around 90%, the rest is silt and clay [10].According to Figure 5, the experimental measured Vp and Vs for sandy hydrate-bearing sediments are higher than field data, while for clayey hydrate-bearing sediments they are lower or close to field data when hydrate saturation is below 25%.Hence the experimental research can supply comparable data for the seismic wave detection of natural gas hydrates.The acoustic velocities measured in our experiment were also compared with the field data of Mallik 2L-38 well in Northern Canada [35], where the occurrence of natural gas hydrates has been verified.The gravel and sand content of natural samples obtained from Mallik 2L-38 borehole is around 90%, the rest is silt and clay [10].According to Figure 5, the experimental measured V p and V s for sandy hydrate-bearing sediments are higher than field data, while for clayey hydrate-bearing sediments they are lower or close to field data when hydrate saturation is below 25%.Hence the experimental research can supply comparable data for the seismic wave detection of natural gas hydrates.
The physical parameters, such as, elastic modulus, shear modulus, and Poisson's ratio, of recovered cores containing laboratory-formed hydrate were calculated and listed in Table 2.The V p /V s ratio ranges from 1.7 to 1.9 for sandy cores, and 1.55-1.7 for clayey cores.The elastic modulus of hydrate-bearing sediments is an indispensable parameter for seismic inversion models.As shown in Table 2, when hydrate saturation approximates 5-25%, the elastic modulus of sandy hydrate-bearing Energies 2019, 12, 1825 8 of 11 sediments ranges from 6 to 24 GPa, and ranges from 5 to 9 GPa for clayey hydrate-bearing sediments.The results demonstrate that hydrate saturation and types of host sediments are two dominant factors that affect the elastic properties of hydrate-containing sediments, especially for sandy sediments.

Effects of Hydrate Dissociation
The interaction between host sediments and gas hydrates is a quite important issue, not only for developing remote detection technology, but also for potential environment impacts caused by hydrate dissociation from the sediments.Being solid crystal, gas hydrates can be viewed as pore occupants, or be part of rocks that works as skeleton structure of the reservoir [36,37].Once the hydrate is dissociating, it may significantly affect layer mechanical strength and cause serious geological disasters, such as submarine landslide, gas release, and break well systems [38][39][40].In this work, we studied the influence of hydrate decomposition on the physical mechanism of hydrate sediments by comparing the acoustic velocity before and after hydrate decomposition.As shown in Figure 6, the V p and V s of clayey specimens barely change after hydrate decomposition, indicating the layer stability can be maintained.It is consistent with the actual result of field hydrates production at Shenhu area in the South China Sea, which was conducted in May, 2017.Hydrate-bearing sediments in the Shenhu area are mainly composed of muddy silt [2], and there were no collapses or landslides duringthe 30-day gas production period.For sandy specimens, the V p and V s decrease significantly after hydrate decomposition, especially when hydrate saturation is above 10%, as shown in Figure 6.Both the Mallik onshore production test and Japan Nankai Trough offshore trial found the fluid migration conduits were collapsed or deformed as the hydrates dissociated, causing the subsequent decline in production rate [41,42].The host sediments of those two sites are coarse sand.Therefore, the risk of submarine landslide and layer collapse is higher for the development of sandy hydrate reservoirs.
recovered cores containing laboratory-formed hydrate were calculated and listed in Table 2.The Vp/Vs ratio ranges from 1.7 to 1.9 for sandy cores, and 1.55-1.7 for clayey cores.The elastic modulus of hydrate-bearing sediments is an indispensable parameter for seismic inversion models.As shown in Table 2, when hydrate saturation approximates 5-25%, the elastic modulus of sandy hydrate-bearing sediments ranges from 6 to 24 GPa, and ranges from 5 to 9 GPa for clayey hydrate-bearing sediments.The results demonstrate that hydrate saturation and types of host sediments are two dominant factors that affect the elastic properties of hydrate-containing sediments, especially for sandy sediments.

Effects of Hydrate Dissociation
The interaction between host sediments and gas hydrates is a quite important issue, not only for developing remote detection technology, but also for potential environment impacts caused by hydrate dissociation from the sediments.Being solid crystal, gas hydrates can be viewed as pore occupants, or be part of rocks that works as skeleton structure of the reservoir [36,37].Once the hydrate is dissociating, it may significantly affect layer mechanical strength and cause serious geological disasters, such as submarine landslide, gas release, and break well systems [38][39][40].In this work, we studied the influence of hydrate decomposition on the physical mechanism of hydrate sediments by comparing the acoustic velocity before and after hydrate decomposition.As shown in Figure 6, the Vp and Vs of clayey specimens barely change after hydrate decomposition, indicating the layer stability can be maintained.It is consistent with the actual result of field hydrates production at Shenhu area in the South China Sea, which was conducted in May, 2017.Hydrate-bearing sediments in the Shenhu area are mainly composed of muddy silt [2], and there were no collapses or landslides duringthe 30-day gas production period.For sandy specimens, the Vp and Vs decrease significantly after hydrate decomposition, especially when hydrate saturation is above 10%, as shown in Figure 6.Both the Mallik onshore production test and Japan Nankai Trough offshore trial found the fluid migration conduits were collapsed or deformed as the hydrates dissociated, causing the subsequent decline in production rate [41,42].The host sediments of those two sites are coarse sand.Therefore, the risk of submarine landslide and layer collapse is higher for the development of sandy hydrate reservoirs.2.  2.

Conclusions
Seismic wave detection is an effective method to identify subsea hydrate accumulation zones and saturation, but corresponding data are limited because of the high cost of recovering natural gas hydrate cores at in situ conditions.This work performed a series of experiments to study the acoustic properties of submarine recovered cores that contain laboratory-formed hydrate.Compared with effective confining stress and other pore occupants, the amount of gas hydrate obviously plays a dominant role in the physical properties of host sediments.The influence of buried depth of natural gas hydrates on the seismic velocity can be neglected if the thickness of hydrate layer is within 100 m.The acoustic velocity of sandy specimens is much higher than that of clayey specimens when at the same hydrate saturation.Conversely, the dissociation of gas hydrate can significantly change the strength of the skeleton structure of sandy hydrate-bearing sediments.Therefore, the risk of submarine landslide and layer collapse may be higher for the development of sandy hydrate reservoirs.

Figure 1 .Figure 1 .
Figure 1.Details of the water-adding device.(a) The picture of a bracket and cylindrical vessel; (b) the bracket was placed into the cylindrical vessel for supporting the cores when adding water vapor.

Figure 1 .Figure 2 .
Figure 1.Details of the water-adding device.(a) The picture of a bracket and cylindrical vessel; (b) the bracket was placed into the cylindrical vessel for supporting the cores when adding water vapor.

Figure 3 .
Figure 3. Wave diagram of P and S-wave signals collected in the experiment.

Figure 3 .
Figure 3. Wave diagram of P and S-wave signals collected in the experiment.

Figure 4
Figure 4The relationship between acoustic velocities, P-wave velocity (a) and S-wave velocity (b), and effective confining pressure, pore occupants for core S1.

Figure 4 .
Figure 4.The relationship between acoustic velocities, P-wave velocity (a) and S-wave velocity (b), and effective confining pressure, pore occupants for core S1.

Figure 6 .
Figure 6.Comparison of the acoustic velocities of hydrate-bearing sediments before and after hydrate dissociation.(a):P-wave velocity, (b)S-wave velocity.Hydrate saturations are shown in Table2.

Figure 6 .
Figure 6.Comparison of the acoustic velocities of hydrate-bearing sediments before and after hydrate dissociation.(a) P-wave velocity, (b) S-wave velocity.Hydrate saturations are shown in Table2.

Table 1 .
The physcial parameters of natural recovered cores.

Table 2 .
Theacoustic parameters of hydrate-bearing sediments formed with natural recovered cores.

Table 2 .
Theacoustic parameters of hydrate-bearing sediments formed with natural recovered cores.
Note: the absolute errors of P and S-wave velocity are both 0.01 km/s, the combined relative errors of elastic modulus, bulk modulus, shear modulus, and Poisson's ratio are 9.6%, 4.7%, 1.8%, and 11.0%, respectively.