the effect of hydrate content on seismic attenuation a case study for mallik 2l 38 well data mackenzie delta canada 1 observations of velocities in sediments containing gas hydrates show that the strength of sediments increases with hydrate saturation hence it is expected that the attenuation of these sediments will decrease with increasing hydrate saturation however sonic log measurements in the mallik 2l 38 well and cross hole tomography measurements in the mallik field have shown that attenuation increases with hydrate saturation we studied a range of mechanisms by which increasing hydrate saturation could cause increased attenuation we found that a difference in permeability between the host sediment and the newly formed hydrate can produce the observed effect we modelled attenuation in terms of biot and squirt flow mechanisms in composite media we have used our model to predict observed attenuations in the mallik 2l 38 well mackenzie delta canada

gas hydrate reservoir and active methane venting province in sediments on 20 ma young oceanic crust in the fram strait offshore nw svalbard seafloor pockmarks are common indicators for vertical fluid flow and frequently associated with methane discharge through the gas hydrate stability zone ghsz the present day flux through these degassing systems is presumably at a low level on most rifted continental margins a pockmark field on the nw svalbard passive margin is located on young ocean crust 20 ma and shows evidence of ongoing episodic degassing new geophysical data from the vestnesa ridge 79 n a mounded and elongated sediment drift in the eastern fram strait reveal a gas hydrate free gas and venting system that is exceptionally more dynamic than documented elsewhere along the northeastern north atlantic margin the prominent bottom simulating reflection bsr about 200 mbsf separates anomalously high p wave velocities in the ghsz from a remarkable underlying low velocity zone indicating the presence of gas hydrate and gas in the pore space inversion of p wave velocity data using the differential effective medium theory yields a two dimensional concentration model of methane hydrate and free gas the model predicts saturations of up to 11 in the hydrate reservoir which due to the seafloor topography forms a large anticlinal permeability barrier below in the low velocity zone i e 1350 1500 m s up to 3 of free gas is predicted across the apex of the vestnesa ridge and in the immediate vicinity of extensional faults a conservative estimate indicates that 225 kg m2 of pure methane is stored in hydrate and gas in the upper 230 m of the sedimentary column an elongated pockmark field consisting of 100 individual pockmarks up to 600 m wide systematically aligns the apex of the vestnesa ridge active vigorous degassing from the topography controlled pressure valve system was evident from a 750 m high and 150 m wide gas flare observed in the water column during a cruise with r v jan mayen in october 2008 the gas flare documents dynamic degassing through the corresponding chimney which penetrates the entire ghsz and into the underlying free gas zone cruises in 2006 and 2007 did not detect active gas venting above the pockmark field accordingly vigorous degassing may operate in an episodic mode where hydrothermal circulation systems through young ocean crust may play a significant role

attenuation of seismic waves in methane gas hydrate bearing sand compressional wave p wave and shear wave s wave velocities vp and vs respectively from remote seismic methods have been used to infer the distribution and volume of gas hydrate within marine sediments recent advances in seismic methods now allow compressional and shear wave attenuations q1p and q1s respectively to be measured however the interpretation of these data is problematic due to our limited understanding of the effects of gas hydrate on physical properties therefore a laboratory gas hydrate resonant column was developed to simulate pressure and temperature conditions suitable for methane gas hydrate formation in sand specimens and the subsequent measurement of both q1p and q1s at frequencies and strains relevant to marine seismic surveys 13 dry gas saturated sand specimens were investigated with different amounts of methane gas hydrate evenly dispersed throughout each specimen the results show that for these dry specimens both q1p and q1s are highly sensitive to hydrate saturation with unexpected peaks observed between 3 and 5 per cent hydrate saturation it is thought that viscous squirt flow of absorbed water or free gas within the pore space is enhanced by hydrate cement at grain contacts and by the nanoporosity of the hydrate itself these results show for the first time the dramatic effect methane gas hydrate can have on seismic wave attenuation in sand and provide insight into wave propagation mechanisms these results will aid the interpretation of elastic wave attenuation data obtained using marine seismic prospecting methods

attenuation measurements from sonic waveform logs in methane hydrate bearing sediments at the nankai trough exploratory well off tokai central japan 1 i have used full waveform logs from the nankai trough exploratory well off central japan to estimate both compressional and shear attenuation in sediments containing methane hydrate mh the attenuation estimates are based on a median frequency shift to the amplitude spectrum of the recorded waveforms this paper is concerned with attenuation at sonic frequencies of 10 20 khz for compressional waves and 500 1000 hz for shear waves i observed that the presence of mh increases the host sediments seismic attenuation and shear attenuation may be more helpful than compressional attenuation in detecting or characterizing mh bearing sediments as compressional attenuation is affected by the presence of gas moreover the ratio of compressional to shear attenuation is found to be a more sensitive indicator of the presence of low saturation gas than the corresponding velocity ratio

velocity structure of a gas hydrate reflector seismic reflection profiles across many continental margins have imaged bottom simulating reflectors bsrs parallel to the seabed these are often interpreted as the base of a zone in which methane hydrate ice is stable waveform inversion of seismic reflection data can be used to estimate from seismic data worldwide the velocity structure of a bsr and its thickness a test of this method at a drill site of the ocean drilling program predicts that sediment pores beneath the bsr contain free methane for approximately 30 meters the hydrate and underlying gas represent a large global reservoir of methane which may have economic importance and may influence global climate

influence of gas hydrate morphology on the seismic velocities of sands this paper reports the results of a series of resonant column tests on specimens where gas hydrate has been formed in sands using an excess water technique in these specimens the amount of hydrate formed is restricted by the amount of gas in the specimen and with an excess of water being present in the pore space results of resonant column tests carried out to determine compressional and shear wave velocities suggest that gas hydrate formed in this way are frame supporting in contrast the behavior observed in sands where the hydrate is formed from finite water where the remaining pore space is saturated with methane gas termed in this paper the excess gas method exhibits a cementing behavior while tetrahydrofuran hydrate sands or where the hydrate is formed from dissolved methane within the pore water exhibit a pore filling behavior for hydrate saturations less than 40 for sands where the hydrate is formed using the excess water method much larger volumes of hydrate are required before a significant increase in shear wave velocity occurs although increases in compressional wave velocity are seen at lower hydrate contents these results suggest that hydrate interaction with the sediment is strongly dependent on morphology and that natural hydrate may exhibit contrasting seismic signatures depending upon the geological environment in which it forms

mechanical properties of sand silt and clay containing tetrahydrofuran hydrate 1 the mechanical behavior of hydrate bearing sediments subjected to large strains has relevance for the stability of the seafloor and submarine slopes drilling and coring operations and the analysis of certain small strain properties of these sediments for example seismic velocities this study reports on the results of comprehensive axial compression triaxial tests conducted at up to 1 mpa confining pressure on sand crushed silt precipitated silt and clay specimens with closely controlled concentrations of synthetic hydrate the results show that the stress strain behavior of hydrate bearing sediments is a complex function of particle size confining pressure and hydrate concentration the mechanical properties of hydrate bearing sediments at low hydrate concentration probably 50 of pore space the behaviorbecomesmoreindependentofstressbecausethehydratescontrolbothstiffnessand strength and possibly the dilative tendency of sediments by effectively increasing interparticle coordination cementing particles together and filling the pore space the cementation contribution to the shear strength of hydrate bearing sediments decreases with increasing specific surface of soil minerals the lower the effective confining stress the greater the impact of hydrate formation on normalized strength

characterization of a stratigraphically constrained gas hydrate system along the western continental margin of svalbard from ocean bottom seismometer data the ongoing warming of bottom water in the arctic region is anticipated to destabilize some of the gas hydrate present in shallow seafloor sediment potentially causing the release of methane from dissociating hydrate into the ocean and the atmosphere ocean bottom seismometer obs experiments were conducted along the continental margin of western svalbard to quantify the amount of methane present as hydrate or gas beneath the seabed p and s wave velocities were modeled for five sites along the continental margin using ray trace forward modeling two southern sites were located in the vicinity of a 30 km long zone where methane gas bubbles escaping from the seafloor were observed during the cruise the three remaining sites were located along an e w orientated line in the north of the margin at the deepest northern site vp anomalies indicate the presence of hydrate in the sediment immediately overlying a zone containing free gas up to 100 m thick the acoustic impedance contrast between the two zones forms a bottom simulating reflector bsr at approximately 195 m below the seabed the two other sites within the gas hydrate stability zone ghsz do not show the clear presence of a bsr or of gas hydrate however anomalously low vp indicating the presence of free gas was modeled for both sites the hydrate content was estimated from vp and vs using effective medium theory at the deepest northern site modeling suggests a pore space hydrate concentration of 7 12 if hydrate forms as part of a connected framework and about 22 if it is pore filling at the two other northern sites located between the deepest site and the landward limit of the ghsz we suggest that hydrate is present in the sediment as inclusions hydrate may be present in small quantities at these two sites 4 5 of the pore space the variation in lithology for the three sites indicated by high resolution seismic profiles may control the distribution concentration and formation of hydrate and free gas

physical properties and rock physics models of sediment containing natural and laboratory formed methane gas hydrate this paper presents results of shear strength and acoustic velocity p wave measurements performed on 1 samples containing natural gas hydrate from the mallik 2l 38 well mackenzie delta northwest territories 2 reconstituted ottawa sand samples containing methane gas hydrate formed in the laboratory and 3 ice bearing sands these measurements show that hydrate increases shear strength and p wave velocity in natural and reconstituted samples the proportion of this increase depends on 1 the amount and distribution of hydrate present 2 differences in sediment properties and 3 differences in test conditions stress strain curves from the mallik samples suggest that natural gas hydrate does not cement sediment grains however stress strain curves from the ottawa sand containing laboratory formed gas hydrate do imply cementation is present acoustically rock physics modeling shows that gas hydrate does not cement grains of natural mackenzie delta sediment natural gas hydrates are best modeled as part of the sediment frame this finding is in contrast with direct observations and results of ottawa sand containing laboratory formed hydrate which was found to cement grains waite et al 2004 it therefore appears that the microscopic distribution of gas hydrates in sediment and hence the effect of gas hydrate on sediment physical properties differs between natural deposits and laboratory formed samples this difference may possibly be caused by the location of water molecules that are available to form hydrate models that use laboratory derived properties to predict behavior of natural gas hydrate must account for these differences

characterization of in situ elastic properties of gas hydrate bearing sediments on the blake ridge during ocean drilling program leg 164 shear sonic velocity and other geophysical logs were acquired in gas hydrate bearing sediments on the blake ridge to characterize the very distinct seismic signature of such formations anomalous low amplitudes overlying a strong bottom simulating reflector bsr a comparison of the bulk moduli derived from the logs to standard elastic consolidation models shows that the sediments are overconsolidated above the bsr at 440 meters below seafloor mbsf because of the presence of hydrates below the bottom of the thermodynamic hydrate stability zone at 520 mbsf the high compressibility of the formation and the attenuation of the monopole sonic waveforms are typical of sediments partially saturated with free gas between these two depths gas hydrate and free gas seem to coexist within the gas hydrate stability zone we estimate the amount of gas hydrates using different models based on theories for wave scattering in multiphase media and for grain cementation in close agreement with measurements made on discrete in situ samples the latter describes most accurately the interactions between the matrix the pore fluids and the hydrates this model indicates that 5 to 10 of the pore space is occupied by hydrates deposited uniformly on the surface of the grains the comparison with gassmann s model also show that the amount of free gas below the bsr never exceeds 5 of the pore space but is high enough to generate the bsr the coexistence of free gas and gas hydrates below the bsr may be explained by capillary effects in the smaller pores or by remaining crystalline structures after partial hydrate decomposition

development of three dimensional basement structure in taiwan deduced from past plate motion consistency with the present seismicity 1 using colored clay we examined geometrical evolution of the three dimensional basement structure in taiwan due to relative plate motion under the following assumptions the motion of the philippine sea plate relative to the eurasian plate has been constant during the last 15 myr the ryukyu and manila trenches had been connected by a transform fault before the collision between the luzon arc and eurasian continental margin and the collision started at 5 ma the basement structure obtained from the clay model was fully consistent with observed deep seismicity which led us to the following conclusions there are no oceanic slabs beneath central taiwan instead the eurasian continental margin is underthrusting beneath it the philippine sea slab extends northwestward from the ryukyu trench while the south china sea slab extends southeastward from the manila trench the downdip length of the south china sea slab becomes shorter from south to north

effective medium modeling of gas hydrate filled fractures using the sonic log in the krishna godavari basin offshore eastern india 1 indian national gas hydrate program expedition 01 has established that clay rich marine sediment from the krishna godavari kg basin in the eastern indian margin hosts one of the richest gas hydrate deposits in the world resistivity at bit images and pressure cores reveal that the gas hydrate morphology in clay rich sediment varies from complex vein structures grain displacing to invisible pore filling existing rock physics models which relate acoustic data to in situ gas hydrate concentrations generally assume isotropic pore filling gas hydrate which yields misleading concentration estimates for fractured fine grained sediments the anisotropic kg basin sediment presents additional complications here we apply differential effective medium theory to incorporate grain displacing morphologies by which gas hydrate is included as vertical ellipsoids with aspect ratios ranging from those of thin veins up to those of nodules in an elastic anisotropic background we have estimated gas hydrate concentrations from sonic velocities at hole 10d in the kg basin considering three basic gas hydrate morphologies i pore filling ii grain displacing and iii a combination of grain displacing and pore filling average gas hydrate saturations for these three cases are 35 42 27 30 and 33 41 of the total porosity respectively in the depth range 60 140 m below seafloor mbsf saturation is highest at 67 mbsf for any morphology but the values differ between morphologies for the pore filling morphology the maximum gas hydrate saturation of 56 is 18 22 higher than the grain displacing morphology and 2 9 higher than the combined morphology estimates differ by 6 of the sediment volume with rotations of gas hydrate veins from vertical to horizontal

sonic waveform attenuation in gas hydrate bearing sediments from the mallik 2l 38 research well mackenzie delta canada 1 the mallik 2l 38 research well was drilled to 1150 m under the mackenzie delta canada and penetrated a subpermafrost interval where methane hydrate occupies up to 80 of the pore space a suite of high quality downhole logs was acquired to measure in situ the physical properties of these hydrate bearing sediments similar to other hydrate deposits resistivity and compressional and shear sonic velocity data increase with higher hydrate saturation owing to electrical insulation of the pore space and stiffening of the sediment framework in addition sonic waveforms show strong amplitude losses of both compressional and shear waves in intervals where methane hydrate is observed we use monopole and dipole waveforms to estimate compressional and shear attenuation comparing with hydrate saturation values derived from the resistivity log we observe a linear increase in both attenuation measurements with increasing hydrate saturation which is not intuitive for stiffening sediments numerical modeling of the waveforms allows us to reproduce the recorded waveforms and illustrate these results we also use a model for wave propagation in frozen porous media to explain qualitatively the loss of sonic waveform amplitude in hydrate bearing sediments we suggest that this model can be improved and extended allowing hydrate saturation to be quantified from attenuation measurements in similar environments and providing new insight into how hydrate and its sediment host interact

hydrate morphology physical properties of sands with patchy hydrate saturation 1 the physical properties of gas hydrate bearing sediments depend on the volume fraction and spatial distribution of the hydrate phase the host sediment grain size and the state of effective stress determine the hydrate morphology in sediments this information can be used to significantly constrain estimates of the physical properties of hydrate bearing sediments including the coarse grained sands subjected to high effective stress that are of interest as potential energy resources reported data and physical analyses suggest hydrate bearing sands contain a heterogeneous patchy hydrate distribution whereby zones with 100 pore space hydrate saturation are embedded in hydrate free sand accounting for patchy rather than homogeneous hydrate distribution yields more tightly constrained estimates of physical properties in hydrate bearing sands and captures observed physical property dependencies on hydrate saturation for example numerical modeling results of sands with patchy saturation agree with experimental observation showing a transition in stiffness starting near the series bound at low hydrate saturations but moving toward the parallel bound at high hydrate saturations the hydrate patch size itself impacts the physical properties of hydrate bearing sediments for example at constant hydrate saturation we find that conductivity electrical hydraulic and thermal increases as the number of hydrate saturated patches increases this increase reflects the larger number of conductive flow paths that exist in specimens with many small hydrate saturated patches in comparison to specimens in which a few large hydrate saturated patches can block flow over a significant cross section of the specimen

poroelastic model to relate seismic wave attenuation and dispersion to permeability anisotropy a transversely isotropic model with a horizontal axis of symmetry based on the biot and squirt flow mechanisms predicts seismic waves in poroelastic media the model estimates velocity dispersion and attenuation of waves propagating in the frequency range of crosswell and high resolution reverse vertical seismic profiling vsp 250 1250 hz for vertical permeability values much greater than horizontal permeability parameters the model assumes the principal axes of the stiffness constant tensor are aligned with the axes of the permeability and squirt flow tensors in addition the unified biot and squirt flow mechanism bisq model is adapted to simulate cracks in permeable media under these conditions the model simulations demonstrate that the preferential direction of fluid flow in a reservoir containing fluid filled cracks can be determined by analyzing the phase velocity and attenuation of seismic waves propagating at different azimuth and incident angles as a result the fast compressional wave

an effective medium inversion algorithm for gas hydrate quantification and its application to laboratory and borehole measurements of gas hydrate bearing sediments the presence of gas hydrate in marine sediments alters their physical properties in some circumstances gas hydrate may cement sediment grains together and dramatically increase the seismic p and s wave velocities of the composite medium hydrate may also form a load bearing structure within the sediment microstructure but with different seismic wave attenuation characteristics changing the attenuation behaviour of the composite here we introduce an inversion algorithm based on effective medium modelling to infer hydrate saturations from velocity and attenuation measurements on hydrate bearing sediments the velocity increase is modelled as extra binding developed by gas hydrate that strengthens the sediment microstructure the attenuation increase is modelled through a difference in fluid flow properties caused by different permeabilities in the sediment and hydrate microstructures we relate velocity and attenuation increases in hydrate bearing sediments to their hydrate content using an effective medium inversion algorithm based on the self consistent approximation sca differential effective medium dem theory and biot and squirt flow mechanisms of fluid flow the inversion algorithm is able to convert observations in compressional and shear wave velocities and attenuations to hydrate saturation in the sediment pore space we applied our algorithm to a data set from the mallik 2l 38 well mackenzie delta canada and to data from laboratory measurements on gas rich and water saturated sand samples predictions using our algorithm match the borehole data and water saturated laboratory data if the proportion of hydrate contributing to the load bearing structure increases with hydrate saturation the predictions match the gas rich laboratory data if that proportion decreases with hydrate saturation we attribute this difference to differences in hydrate formation mechanisms between the two environments

compressional and shear wave velocities in uncemented sediment containing gas hydrate 1 the competing hypotheses for gas hydrate formation at the particle scale in sediments describe processes of pore filling frame building or cementation new measurements of compressional vp and shear wave vs velocities in fine grained sands subjected to low confinement and monitored during formation of tetrahydrofuran hydrate indicate that hydrate nucleates in the pore space presumably at grain boundaries and grows with limited impact on the sediment shear stiffness vp and vs until crystals begin to interact with the granular skeleton at 40 hydrate concentration vs increases significantly more than vp at higher hydrate concentrations reflecting larger changes in the specimen s shear stiffness than its bulk stiffness the results indicate that seismic velocities and or their ratio vp vs have limited capability for locating hydrate or constraining hydrate concentrations

anisotropic effective medium modeling of the elastic properties of shales shales are complex porous materials normally consisting of percolating and interpenetrating fluid and solid phases the solid phase is generally comprised of several mineral components and forms an intricate and anisotropic microstructure the shape orientation and connection of the two phases control the anisotropic elastic properties of the composite solid we develop a theoretical framework that allows us to predict the effective elastic properties of shales its usefulness is demonstrated with numerical modeling and by comparison with established ultrasonic laboratory experiments the theory is based on a combination of anisotropic formulations of the self consistent sca and differential effective medium dem approximations this combination guarantees that both the fluid and solid phases percolate at all porosities our modeling of the elastic properties of shales proceeds in four steps first we consider the case of an aligned biconnected clay fluid composite composed of ellipsoidal inclusion

a laboratory investigation into the seismic velocities of methane gas hydrate bearing sand remote seismic methods which measure the compressional wave p wave velocity vp and shear wave s wave velocity vs can be used to assess the distribution and concentration of marine gas hydrates in situ however interpreting seismic data requires an understanding of the seismic properties of hydrate bearing sediments which has proved problematic because of difficulties in recovering intact hydrate bearing sediment samples and in performing valid laboratory tests therefore a dedicated gas hydrate resonant column ghrc was developed to allow pressure and temperature conditions suitable for hydrate formation to be applied to a specimen with subsequent measurement of both vp and vs made at frequencies and strains relevant to marine seismic investigations thirteen sand specimens containing differing amounts of evenly dispersed hydrate were tested the results show a bipartite relationship between velocities and hydrate pore saturation with a marked transition between 3 and 5 hydrate pore saturation for both vp and vs this suggests that methane hydrate initially cements sand grain contacts then infills the pore space these results show in detail for the first time using a resonant column how hydrate cementation affects elastic wave properties in quartz sand this information is valuable for validating theoretical models relating seismic wave propagation in marine sediments to hydrate pore saturation

elastic properties of hydrate bearing sediments using effective medium theory accurate and detailed models of the seismic velocity structure of gas hydrate bearing sediments may be determined by careful analysis of controlled source seismic data however interpretation of these velocities in terms of hydrate saturation of the pore space has hitherto relied on semiempirical formulas and or simple effective medium theory we develop a rigorous theoretical scheme to relate the seismic properties of a clay rich hydrate bearing sediment to its porosity mineralogy microstructural features and hydrate saturation we consider separately the two possible end members for the distribution of hydrate in the pore space 1 hydrates are unconnected and located in the pore voids without appreciable grain contact and 2 connected hydrates are forming cement binding around the grains the scheme is transversely isotropic to allow for anisotropy due to alignment of clay platelets and is based on a combination of a self consistent approximation a differential effective medium theory and a method of smoothing for crystalline aggregates we have applied the scheme to lithological and seismic velocity data from ocean drilling program site 995 on the blake ridge southeastern u s continental margin to make estimates of the hydrate saturation it was found that the hydrates are probably unconnected and their volume concentration varies between approx 0 at 100 m below the seabed and approx 9 at 400 m depth just above the bottom simulating reflector if the clay platelet orientation distribution resembles the function we have used

methane hydrate formation in partially water saturated ottawa sand bulk properties of gas hydrate bearing sediment strongly depend on whether hydrate forms primarily in the pore fluid becomes a load bearing member of the sediment matrix or cements sediment grains our compressional wave speed measurements through partially water saturated methane hydrate bearing ottawa sands suggest hydrate surrounds and cements sediment grains the three ottawa sand packs tested in the gas hydrate and sediment test laboratory instrument ghastli contain 38 1 porosity initially with distilled water saturating 58 31 and 16 of that pore space respectively from the volume of methane gas produced during hydrate dissociation we calculated the hydrate concentration in the pore space to be 70 37 and 20 respectively based on these hydrate concentrations and our measured compressional wave speeds we used a rock physics model to differentiate between potential pore space hydrate distributions model results suggest methane hydrate cements unconsolidated sediment when forming in systems containing an abundant gas phase

seismic characterization of hydrates in faulted fine grained sediments of krishna godavari basin full waveform inversion in fine grained faulted sediments both stratigraphic and fault induced structural variations can simultaneously determine the gas hydrate distribution insights into hydrate distribution can be obtained from p wave velocity vp and attenuation qp 1 character of the gas hydrate stability zone ghsz in this paper we apply frequency domain full waveform inversion fwi to surface towed 2d multichannel seismic data from the krishna godavari kg basin india to image the fine scale 100 x 30 m vp and qp 1 variations within the ghsz we validate the inverted vp model by reconciling it with a sonic log from a nearby 250 m well the vp model shows a patchy distribution of hydrate away from the faults dominated parts of the profile hydrates demonstrate stratigraphic control which appears to be permeability driven the qp 1 model suggests that attenuation is relatively suppressed in hydratesbearing sediments elevated attenuation in non hydrate bearing sediments could be driven by the apparent pore fluid immiscibility at seismic wavelengths the vp and the qp 1 models also suggest that fault zones within the ghsz can be hydrate or free gas rich depending on the relative supply of free gas and water from below the ghsz

characterization of seismic wave velocity and attenuation and interpretation of tetrahydrofuran hydrate bearing sand using resonant column testing abstract understanding the geophysical properties of hydrate bearing sediments is critical for hydrate utilization such as prediction of hydrate distribution and concentration as well as safety evaluation during hydrate production here a customized resonant column is employed to test tetrahydrofuran thf hydrate bearing sandy specimens the results show that wave velocities of all tested specimens increase exponentially with increasing effective confining stress v whereas wave velocity evolution becomes almost stress independent when critical hydrate saturation sh 0 53 in this study is reached which transforms the hydrate morphology from non cementing to cementing type in addition regardless of structure differences between tested specimens an increased factor and decreased exponent follows a uniform equation 0 61 0 2 log m s measured wave velocity differences between hydrate bearing and hydrate free specimens confirm the coexistence of various hydrate morphologies based on matching results of theoretical predictions from rock physics models and interpretation of volumetric proportions of hydrates with different morphologies in tested specimens wave attenuation of tested specimens generally decreases with increasing effective confining stress and increase with increasing hydrate saturation the relationship between wave attenuation and hydrate saturation qs 1 0 015 0 17 sh is confirmed to be a hydrate saturation indicator considering data from field seismic surveys this type of monotonous relationship appears to be less sensitive for hydrate morphology nevertheless through the combined interpretation of wave velocities and attenuations it is theoretically feasible to not only predict hydrate saturation but also illustrate the evolution of hydrate morphology our study contributes to a better understanding of stress related geophysical properties of hydrate bearing sediments which can be useful for the characterization of hydrate reservoirs

sediments with gas hydrates internal structure from seismic avo we interpret amplitude variation with offset avo data from a bottom simulating reflector bsr offshore florida by using rock physics based synthetic seismic models a previously conducted velocity and avo analysis of the in situ seismic data showed that the bsr separates hydrate bearing sediments from sediments containing free methane the amplitudes at the bsr are increasingly negative with increasing offset this behavior was explained by p wave velocity above the bsr being larger than that below the bsr and s wave velocity above the bsr being smaller than that below the bsr we use these avo and velocity results to infer the internal structure of the hydrated sediment to do so we examine two micromechanical models that correspond to the two extreme cases of hydrate deposition in the pore space 1 the hydrate cements grain contacts and strongly reinforces the sediment and 2 the hydrate is located away from grain contacts and does not affect the stiffness of the sediment frame only the second model can qualitatively reproduce the observed avo response thus inferred internal structure of the hydrate bearing sediment means that 1 the sediment above the bsr is uncemented and thereby mechanically weak and 2 its permeability is very low because the hydrate clogs large pore space conduits the latter explains why free gas is trapped underneath the bsr the seismic data also indicate the absence of strong reflections at the top of the hydrate layer this fact suggests that the high concentration of hydrates in the sediment just above the bsr gradually decreases with decreasing depth this effect is consistent with the fact that the low permeability hydrated sediments above the bsr prevent free methane from migrating upwards

gas hydrates free gas distribution and fault pattern on the west svalbard continental margin summary a 3 d tomographic analysis of seismic velocity and attenuation fields is presented with the purpose to analyse the intimate relation of gas hydrates and free gas distribution with the fault pattern the 3 d four component seismic data have been acquired offshore western svalbard the analysis of the subbottom topography of the base of the stability field of gas hydrates indicated by the bottom simulating reflector and the thickness of the underlying free gas bearing zone suggest a fault induced compartmentalization of the gas reservoir the proposed fluid and gas circulation scheme assumes deep sourced warm fluids moving upwards mainly along high permeability faults and fractures below the gas hydrate stability zone faults change from being a pathway of fluid flow to being barriers to fluid flow in the hydrate stability zone consequently below the hydrate stability zone free gas is transferred from the permeable faults to non faulted but fault bounded sediments where it accumulates below the base of the hydrate stability zone producing a significantly thicker free gas layer the tomographic approach is therefore proposed as an effective procedure to provide detailed information on 3 d p and s wave velocity and attenuation distribution in marine sediments and specifically to examine the details of the distribution of the gas hydrates and of the free gas for which the dependence on the vp vs and qp is now fairly well known

gas hydrate stability zone modelling in areas of salt tectonics and pockmarks of the barents sea suggests an active hydrocarbon venting system the barents sea seabed exhibits an area of major glacial erosion exposing parts of the old hydrocarbon basins in this region we modelled the gas hydrate stability field in a 3d perspective including the effects of higher order hydrocarbon gases we used 3d seismic data to analyse the linkage between fluid flow expressions and hydrate occurrences above old sedimentary basin systems and vertical faults pockmarks showed a relation to fault systems where some of them are directly connected to hydrocarbon bearing sedimentary formations the influence of bottom water temperature pore water salinity and geothermal gradient variation on gas hydrate stability zone ghsz thickness is critically analysed in relation to both geological formations and salt tectonics our analysis suggests a highly variable ghsz in the barents sea region controlled by local variations in the parameters of stability conditions recovery of gas hydrate sample from the region and presence of gas enhanced reflections below estimated bsr depths may indicate a prevalent gas hydrate stable condition

elastic wave velocity in marine sediments with gas hydrates effective medium modeling we offer a first principle based effective medium model for elastic wave velocity in unconsolidated high porosity ocean bottom sediments containing gas hydrate the dry sediment frame elastic constants depend on porosity elastic moduli of the solid phase and effective pressure elastic moduli of saturated sediment are calculated from those of the dry frame using gassmann s equation to model the effect of gas hydrate on sediment elastic moduli we use two separate assumptions a hydrate modifies the pore fluid elastic properties without affecting the frame b hydrate becomes a component of the solid phase modifying the elasticity of the frame the goal of the modeling is to predict the amount of hydrate in sediments from sonic or seismic velocity data we apply the model to sonic and vsp data from odp hole 995 and obtain hydrate concentration estimates from assumption b consistent with estimates obtained from resistivity chlorinity and evolved gas data

relationships among compressional wave attenuation porosity clay content and permeability in sandstones anelastic attenuation is the process by which rocks convert compressional waves into heat and thereby modify the amplitude and phase of the waves understanding the causes of compressional wave attenuation is important in the acquisition processing and interpretation of high resolution seismic data and in deducing the physical properties of rocks from seismic data we have measured the attenuation coefficients of compressional waves in 42 sandstones at a confining pressure of 40 mpa equivalent to a depth of burial of about 1 5 km in a frequency range from 0 5 to 1 5 mhz the compressional wave measurements were made using a pulse echo method in which the sample 5 cm diameter 1 8 cm to 3 5 cm long was sandwiched between perspex lucite buffer rods inside the high pressure rig the attenuation of the sample was estimated from the logarithmic spectral ratio of the signals corrected for beam spreading reflected from the top and base of the sample the results show that for these samples compressio

seismic velocities for hydrate bearing sediments using weighted equation a weighted equation based on the three phase time average and wood equations is applied to derive a relationship between the compressional wave p wave velocity and the amount of hydrates filling the pore space the proposed theory predicts accurate p wave velocities of marine sediments in the porosity range of 40 80 and provides a practical means of estimating the amount of in situ hydrate using seismic velocity the shear s wave velocity is derived under the assumption that the p to s wave velocity ratio of the hydrated sediments is proportional to the weighted average of the p to s wave velocity ratios of the constituent components of the sediment in the case that all constituent components are known a weighted equation using multiphase time average and wood equations is possible however this study showed that a three phase equation with modified matrix velocity compensated for the clay content is sufficient to accurately predict the compressional wave velocities for the marine sediments this theory was applied to the laboratory measurements of the p and s wave velocities in permafrost samples to infer the amount of ice in the unconsolidated sediment the results are comparable to the results obtained by repeatedly applying the two phase wave scattering theory the theory predicts that the poisson s ratio of the hydrated sediments decreases as the hydrate concentration increases and the porosity decreases in consequence the amplitude versus offset avo data for the bottom simulating reflections may reveal positive negative or no avo anomalies depending on the concentration of hydrates in the sediments

theoretical modeling insights into elastic wave attenuation mechanisms in marine sediments with pore filling methane hydrate the majority of presently exploitable marine methane hydrate reservoirs are likely to host hydrate in disseminated form in coarse grain sediments for hydrate concentrations below 25 40 disseminated or pore filling hydrate does not increase elastic frame moduli thus making impotent traditional seismic velocity based methods here we present a theoretical model to calculate frequency dependent p and s wave velocity and attenuation of an effective porous medium composed of solid mineral grains methane hydrate methane gas and water the model considers elastic wave energy losses caused by local viscous flow both i between fluid inclusions in hydrate and pores and ii between different aspect ratio pores created when hydrate grows the inertial motion of the frame with respect to the pore fluid biot s type fluid flow and gas bubble damping the sole presence of pore filling hydrate in the sediment reduces the available porosity and intrinsic permeability of the sediment affecting biot s type attenuation at high frequencies our model shows that attenuation maxima due to fluid inclusions in hydrate are possible over the entire frequency range of interest to exploration seismology 1 106 hz depending on the aspect ratio of the inclusions whereas maxima due to different aspect ratio pores occur only at sonic to ultrasound frequencies 104 106 hz this frequency response imposes further constraints on possible hydrate saturations able to reproduce broadband elastic measurements of velocity and attenuation our results provide a physical basis for detecting the presence and amount of pore filling hydrate in seafloor sediments using conventional seismic surveys

estimating the amount of gas hydrate and free gas from marine seismic data marine seismic data and well log measurements at the blake ridge offshore south carolina show that prominent seismic bottom simulating reflectors bsrs are caused by sediment layers with gas hydrate overlying sediments with free gas we apply a theoretical rock physics model to 2 d blake ridge marine seismic data to determine gas hydrate and free gas saturation high porosity marine sediment is modeled as a granular system where the elastic wave velocities are linked to porosity effective pressure mineralogy elastic properties of the pore filling material and water gas and gas hydrate saturation of the pore space to apply this model to seismic data we first obtain interval velocity using stacking velocity analysis next all input parameters to the rock physics model except porosity and water gas and gas hydrate saturation are estimated from geologic information to estimate porosity and saturation from interval velocity we first assume that the entire sediment does not contain gas hydrate o

overall properties of heterogeneous material summary small scale heterogeneity in crustal structure can be very complex and difficult to describe in detail and yet at the same time can be very important for the description of for instance tectonic stress and porosity statistical properties of such heterogeneity can be derived from the properties of waves of relatively large wavelength as they propagate through the crust the differences between measured wavelengths and attenuation coefficients and those of compact rock and the variations of these quantities in space and time provide in principle a means of determining quantities like the density and orientations of microfractures and the nature of crack infill material the theoretical basis for inferences of this kind is the concept of the effective or equivalent material based on an averaging process taken over the microstructure a number of methods have been used to calculate the properties of the effective medium several of which are described here

estimation of gas hydrate concentration from multi component seismic data at sites on the continental margins of nw svalbard and the storegga region of norway high resolution seismic experiments employing arrays of closely spaced four component ocean bottom seismic recorders were conducted at a site off western svalbard and a site on the northern margin of the storegga slide off norway to investigate how well seismic data can be used to determine the concentration of methane hydrate beneath the seabed data from p waves and from s waves generated by p s conversion on reflection were inverted for p and s wave velocity vp and vs using 3d travel time tomography 2d ray tracing inversion and 1d waveform inversion at the nw svalbard site positive vp anomalies above a sea bottom simulating reflector bsr indicate the presence of gas hydrate a zone containing free gas up to 150 m thick lying immediately beneath the bsr is indicated by a large reduction in vp without significant reduction in vs at the storegga site the lateral and vertical variation in vp and vs and the variation in amplitude and polarity of reflectors indicate a heterogeneous distribution of hydrate that is related to a stratigraphically mediated distribution of free gas beneath the bsr derivation of hydrate content from vp and vs was evaluated using different models for how hydrate affects the seismic properties of the sediment host and different approaches for estimating the background velocity of the sediment host the error in the average vp of an interval of 20 m thickness is about 2 5 at 95 confidence and yields a resolution of hydrate concentration of about 3 if hydrate forms a connected framework or about 7 if it is both pore filling and framework forming at nw svalbard in a zone about 90 m thick above the bsr a biot theory based method predicts hydrate concentrations of up to 11 of pore space and an effective medium based method predicts concentrations of up to 6 if hydrate forms a connected framework or 12 if hydrate is both pore filling and framework forming at storegga hydrate concentrations of up to 10 or 20 were predicted depending on the hydrate model in a zone about 120 m thick above a bsr with seismic techniques alone we can only estimate with any confidence the average hydrate content of broad intervals containing more than one layer not only because of the uncertainty in the layer by layer variation in lithology but also because of the negative correlation in the errors of estimation of velocity between adjacent layers in this investigation an interval of about 20 m thickness equivalent to between 2 and 5 layers in the model used for waveform inversion was the smallest within which one could sensibly estimate the hydrate content if lithological layering much thinner than 20 m thickness controls hydrate content then hydrate concentrations within layers could significantly exceed or fall below the average values derived from seismic data