Analysis of rock drainage and cooling experiments for underground cryogenic LNG storage

Hydrogeological monitoring was conducted around a pilot cavern for underground cryogenic LNG (Liquefied Natural Gas) storage. The monitoring was mainly focused on the operation of a drainage and recharge system. After the operation of the drainage system commenced, the drainage rate decreased rapidly in the initial stages and then decreased gradually. Hydrogeological monitoring revealed that the rock drainage system operated effectively. During drainage, the water table was maintained below the cavern roof. The recharge for ice-ring formation was performed in two phases. The first phase involved the cessation of pumping in downward-drainage holes and the second involved the closure of upward boreholes. Since the water table was maintained below the cavern roof, artificial recharge was planned at first. However, it was not implemented due to heavy rainfall in the recharge stage. On the basis of hydrogeological monitoring and hydraulic tests, it was found that the fractures above the roof and on the right wall of the pilot cavern mainly affected seepage into the cavern and thermal variation due to the storage of liquid nitrogen. Thermal variation was examined by the thermometers installed around the pilot cavern. The cooling and thawing processes reveal the characteristics of thermal distribution in the rock and the 0 °C isotherm. The cooling phase lasted for six months, and the 0 °C isotherm progressed in time after the injection of liquid nitrogen into the cavern. The isotherm propagated up to about 4 m from the floor and the sidewall of the cavern and about 3 m from the cavern roof. The cooling rate of the rock mass above the cavern roof was lower than that of the other cavern sides due to the gaseous space in the upper part of the containment. The fractures were analyzed and considered for thermal modeling. A two-dimensional finite element analysis was performed to compare the field monitoring at the pilot cavern. The numerical modeling shows the distance between the ice ring and heat transfer pattern of the fractures around the pilot cavern. The propagation of the measured and calculated 0 °C isotherm reveals that the water-conveying joint on the right wall might affect thermal propagation through a thermal pipe.     

Engineering geological aspects of replacing a solid waste disposal site with a sanitary landfill

The current solid waste disposal site in the Mamak district of Ankara is being engulfed by the growing city. All varieties of solid wastes, including medical wastes, are stored at the present site in an irregular manner. Topographical and geological conditions at Mamak waste site are favorable for constructing a sanitary landfill. Located at the edge of a topographical depression, the site is underlain by the natural hydraulic barriers such as clay and metagreywacke. The terrestrial clay has a permeability of 10⁻⁷ to 10⁻⁸ cm/s and low to moderate values of CEC. The proposed sanitary landfill to replace the present solid waste site has a capacity of storing solid waste over 50 years. The details of base liner, final cover, toe embankment, and drainage of leachate and gas are presented in the paper.     

Gum rosin (colophony): A suitable material for thermomechanical modelling of the lithosphere

We investigate the use of a ductile material with temperature-sensitive viscosity for thermomechanical modelling of the lithosphere. First, we consider the scaling of mechanical and thermal properties. For a normal field of gravity, the balance of stresses and body forces sets the stress scale, in proportion to the linear dimensions and the densities. The equation of thermal conduction sets the time scale. The activation enthalpy for creep sets the temperature scale; but the thermal expansivity provides an additional constraint on this temperature scale.

Gum rosin appears to be a suitable material for lithospheric modelling. We have measured its flow properties, at various temperatures, in a specially designed rotary viscometer with unusually low machine friction. The rosin is almost Newtonian. Strain rate depends upon stress to the power n, where 1.0 <n < 1.14. The viscosity varies over 5 orders of magnitude, from about 10² Pa s at 80°C, to about 10⁷ Pa s at 40°C. The activation enthalphy is thus about 250 kJ/mol. Measured with a needle probe, the thermal conductivity is 0.113 ± 0.001 W m⁻¹K⁻¹; the thermal diffusivity, (6±3) ×10⁻⁷ m² s⁻¹. Calculated from X-ray profiles, the thermal expansivity is about 3 × 10⁻⁴ K⁻¹. These thermal and mechanical properties make gum rosin suitable for thermomechanical models, where linear dimensions scale down by a factor of 10⁶; time, by 10¹¹; viscosity, by 10¹⁷; and temperature change, by 10¹.     

Seismic activity along the central america volcanic arc: is it related to subduction of the Cocos plate?

We determine seismic strain rate of tectonic earthquakes along the Central America Volcanic Arc. We then compare this result to those obtained from earthquakes related to the convergence of the Cocos and Caribbean plates and to earthquakes in the back-arc region of northern Central America.

The seismic strain-rate tensor for shallow-focus earthquakes along the Central America volcanic arc since 1700, has a compressive eigenvector with a magnitude of 0.7 × 10⁻⁸ year⁻¹, and oriented in a 357° azimuth. The extensive eigenvector is oriented in a 86° azimuth, with a magnitude of 0.82 × 10⁻⁸ year⁻¹. When only Centroid Moment-tensor solutions (CMT) are considered, the respective eigenvectors are 1.2 × 10⁻⁸ year⁻¹ and 1.0 × 10⁻⁸ year⁻¹.

The compressive eigenvector from the seismic strain-rate tensor for earthquakes along the Cocos-Caribbean convergent margin is 2.0 × 10⁻⁸ year⁻¹, plunging at 25°, and oriented in a 29° azimuth. Its magnitude and direction are similar to those of the compressive eigenvector for earthquakes along the volcanic arc. The extensive eigenvector along the convergent margin, on the other hand, has a large vertical component. The compressive and extensive eigevenvectors are 4.9 × 10⁻⁸ year⁻¹ and 4.6 × 10⁻⁸ year⁻¹, using only CMTs as the database.

Earthquakes along the grabens of northern Central America yield a seismic strain-rate tensor whose extensive eigenvector has a magnitude of 2.4 × 10⁻⁸ year⁻¹, oriented in a 109° azimuth. Magnitude and direction are similar to those of the extensive eigenvector for earthquakes along the volcanic arc. The compressive eigenvector along the grabens is practically vertical.

Similarities in magnitudes and directions for compressive and extensive eigenvectors suggest to us that the strain field along the Central America volcanic arc is the result of compression along the convergent Cocos-Caribbean margin, and extension in the back-arc region, along the grabens of northern Central America. This field is resolved as strike-slip faulting along the arc.     

The usability of bottom ash as an engineering material when amended with different matrices

This study aims at investigating the utilization of bottom ash obtained from four different power stations as a construction fill and landfill bottom liner. For the matrix material, commercial powdered bentonite, construction lime and natural clay were used. Compaction tests (Standard Proctor and vibratory hammer) were carried out on the different ratios of bottom ash and matrix material. The optimum water content ranged from 40 to 45% yielding a dry density mostly ca 1 Mg m⁻³. Uniaxial compressive strength of mixtures ranges from 0.1 to 0.5 kgf cm⁻² which showed a 3–20-fold increase when tested on 28-day cured specimens. Triaxial compression tests yield varying rates of shear strength which also showed as high as an 11-fold increase for cured specimens. The hydraulic conductivity of those mixtures was mostly ca 10⁻⁴ cm s⁻¹, which is not considered to be low enough for landfill lining. Leaching tests using deionized water were also performed to investigate the possible effect of leachate produced from the mixtures on the environment. In conclusion a light density and environmentally friendly mixture is determined and proposed as construction filler.     

多孔联合水幕效率试验方法及应用实例

我国地下水封洞库主要修建在结晶岩体中,岩体裂隙发育程度和空间展布的非均质各向异性使得水封洞库的水封条件十分复杂。水幕系统是保证地下水封洞库密封性的核心要素,因此所设计的水幕孔与围岩间的水力联系是否满足密封性需求,需进行检验。前人在进行相关问题研究时往往将地质体简化为均质各向同性介质,个别文献在实验室尺度对地质体各向异性进行讨论,但工程适用性不强。为分析场地尺度岩体渗透性能各向异性对水幕孔水力传输效率的影响,并探究考虑这种影响的水幕系统优化方法,本文借鉴工程经验,引入水幕效率的概念,设计一种多孔联合水幕效率试验,以烟台地下水封洞库工程项目为依托,进行场地试验并利用阈值法得到试验结果,进而对水幕孔进行相应补充以优化水幕系统网络,最终使水幕系统网络在考虑地质体各向异性的实际影响下满足气密性要求。实践证明,该方法能显著加强水幕系统在各向异性场区的水力传输能力,比现行的等间距设置水幕系统钻孔更适于工程实际。     

Strain rate effect on the compressive strength of frozen sand

Cylindrical specimens of fine Ottawa sand (A.S.T.M. designation C-109), compacted at the optimum moisture content and saturated before unidirectional freezing, have been tested in uniaxial compression at a cold room temperature of —5.5°C and strain rates between 10⁻⁷ and 10⁻² s⁻¹. The results agree with an extrapolation of data obtained by Sayles and Epanchin [1], but are much higher than those obtained by both Goughnour and Andersland [2] and Perkins and Ruedrich [3] at strain rates below 10⁻⁵ s⁻¹. There is evidence that this may be due to variation in total moisture (ice) content, the conditions under which the specimens were frozen (closed system or an open system) and to the end effects at the platen—specimen interface.     

Studying the migration behaviour of selenate in Boom Clay by electromigration

For the disposal of high-level waste (HLW) in a deep geological formation as Boom Clay, safety assessment studies have shown that long lived ⁷⁹Se is one of the more critical fission products. Therefore, the knowledge of its migration properties (diffusion, retention) through the geological barrier (Boom Clay) is of paramount importance. The migration behaviour of selenium strongly depends on its speciation. Under reducing conditions, selenide would be the dominant species and selenium migration would mainly be controlled by the low solubility of Se(−II)-bearing minerals. However Se species are often found in redox disequilibrium and more oxidized species might also coexist. Therefore, the study of selenate migration requires attention, as it might be the most mobile selenium species in the host rock. Electromigration experiments performed with a ⁷⁵Se-labeled selenate in Boom Clay indicate a high mobility for this species. The apparent diffusion coefficient (D_app) of selenate in Boom Clay is estimated from electromigration experiments performed under different electric fields. Using two independent approaches, the value of D_app for selenate is shown to fall in the range from 1.7×10⁻¹¹ to 6.2×10⁻¹¹ m² s⁻¹. Moreover, no reduction of selenate in Boom Clay was observed.     

Equilibria in the system Au–Ag–S–fluid: computed and experimental data

We present results of computations on the interaction of solid-phase electrum–argentite–pyrite (weight ratios 210⁻⁵/ 210⁻³/1 and 210⁻⁵/410⁻²/1) association with Cl-containing aqueous moderately acid solutions (0.5m NaCl, pH = 3.08) at 300 °C and 500 bars. These data are a physicochemical basis for predicting the geochemical behavior of Au and Ag during the hydrothermal-metasomatic transformation of Au-Ag-pyrite. We also propose a technique of study of this process based on the phase equilibria of the subsystem Au–Ag–S with the aqueous solution at different liquid/solid (l/s) ratios, with the use of new graphic diagrams. The relationship of the composition of the solid-phase association with l/s ratio in real boundary conditions (Au = 17 ppm, mAu/mAg = 10^–3.57–10^–2.28) is shown. The maximum l/s values for complete leaching of gold and silver (l/smax = 200–800) are estimated. It has been established that argentite is the first to dissolve when mAu/mAg(s) > mAu/mAg(sol), and electrum, when mAu/mAg(s) < mAu/mAg(sol).

The experimental results showed that at 300 °C, the conversion of electrum (N_Au = 300‰) nonequilibrated with pyrite into an Au-richer form (N_Au = 730‰) and argentite follows an intricate kinetic scheme. Using the Pilling-Bedwords kinetic equation for processing data yielded the process rate constant K = 2.8(±0.5)10⁻⁵ g²cm⁻⁴day⁻¹. With this equation, the time of the complete conversion of 200 μm thick flat gold grains is 604 days. These data evidence a significant role of kinetic factors in hydrothermal-metasomatic processes involving native gold, which requires combination of thermodynamic and kinetic approaches on the construction of geologo-genetic models for hydrothermal sulfide formation.     

Water-enhanced dynamic recrystallization and solution transfer in experimentally deformed carnallite

Cylindrical samples of polycrystalline carnallite (KMgCl₃, 6H₂O) were deformed in a triaxial apparatus at 60°C, at confining pressures between 0.1 and 31 MPa and at strain rates between 10⁻⁴ and 10⁻⁸ s⁻¹. In a number of cases, small amounts of saturated carnallite brine were added. Samples without added brine deform by intracrystalline slip, mechanical twinning, cracking, and by frictional sliding on crack surfaces. Stress-strain curves of these samples are strongly dependent on confining pressure. Addition of brine has a dramatic effect on both microstructural development and mechanical properties. Grain-boundary migration is strongly enhanced. At lower strain rates, additional intracrystalline effects start to appear, together with the onset of solution transfer. Rapid compaction in samples deformed with added brine causes high fluid pressures to develop. At higher strain rates addition of brine results in a decrease of the flow stress by a factor of two. This weakening will increase even further at strain rates below about 10⁻⁹ s⁻¹, when solution transfer becomes rate controlling. It is argued that deformation of carnallite in nature is adequately described by the flow law found for samples deformed with added brine.     

Assessment of crushed shales for use as compacted landfill liners

The possibility of using crushed shales as landfill liners is investigated in this study. Two types of shales were studied by performing the following laboratory tests: hydraulic conductivity, compaction, swelling, consolidation, X-ray diffraction (XRD), scanning electron microscope (SEM) and chemical analysis. For both compacted shales, the hydraulic conductivity was in order of 10^− 7 cm/s or less which satisfies the specifications for landfill liners. The results of XRD and SEM support the low values of the hydraulic conductivity. Because of the dominant presence of low-activity kaolinite, there was no significant change in the hydraulic conductivity when the compacted shales are exposed to calcium chloride solution. The compressibility of the compacted clay was low and no serious post-construction settlement is expected. The shear strength of the compacted shales was within the usual expected range for earthen liners and, therefore, should pose no challenges with respect to shear strength. The crushed shales also satisfy the other criteria related to Atterberg limits and grain size.     

Ice lens structures, compression strengths and creep behavior of some synthetic frozen silty soils

This research paper concerns the frozen state properties of several silty soils encountered in the construction of a large LNG storage cavern near Carlstadt, New Jersey.

The creep and compression strengths were correlated with both field and laboratory-measured mechanical properties. These cryogenic mechanical properties were related to ice lens structures, sizes and occurrences, and to heat flow and geological structure effects. Three distinct and separate types of ice lenses were identified. These cryogenic ground data pertain to the prediction of the engineering properties of frozen earth for large structures. Sample size effect problems were investigated, this in relationship to the number and size and orientations of the ice lenses present.     

Single-well and large-scale cross-hole tracer experiments in fractured rocks at two sites in Sweden

As part of the Swedish site investigations and research associated with the disposal of spent nuclear fuel, tracer experiments in deep boreholes have been employed in order to characterize solute transport and hydraulic connectivity in the fractured bedrock at two sites??Forsmark and Laxemar-Simpevarp. Performance and analytical results are presented from a suite of tracer experiments in varying scales, from single-hole injection-withdrawal tests, intermediate-scale tests with sorbing tracers, to large-scale connectivity tests over distances up to several hundred meters in major fracture zones or networks of zones. In addition to demonstrating transport connectivity over large distances, a general result is that the single-hole tests, as well as the cross-hole tests with sorbing tracers, have clearly demonstrated the process of solute retention of water-conductive features in the rock at different scales.     

Determining mineral weathering rates based on solid and solute weathering gradients and velocities: application to biotite weathering in saprolites

Chemical weathering gradients are defined by the changes in the measured elemental concentrations in solids and pore waters with depth in soils and regoliths. An increase in the mineral weathering rate increases the change in these concentrations with depth while increases in the weathering velocity decrease the change. The solid-state weathering velocity is the rate at which the weathering front propagates through the regolith and the solute weathering velocity is equivalent to the rate of pore water infiltration. These relationships provide a unifying approach to calculating both solid and solute weathering rates from the respective ratios of the weathering velocities and gradients. Contemporary weathering rates based on solute residence times can be directly compared to long-term past weathering based on changes in regolith composition. Both rates incorporate identical parameters describing mineral abundance, stoichiometry, and surface area.

Weathering gradients were used to calculate biotite weathering rates in saprolitic regoliths in the Piedmont of Northern Georgia, USA and in Luquillo Mountains of Puerto Rico. Solid-state weathering gradients for Mg and K at Panola produced reaction rates of 3 to 6×10⁻¹⁷ mol m⁻² s⁻¹ for biotite. Faster weathering rates of 1.8 to 3.6×10⁻¹⁶ mol m⁻² s⁻¹ are calculated based on Mg and K pore water gradients in the Rio Icacos regolith. The relative rates are in agreement with a warmer and wetter tropical climate in Puerto Rico. Both natural rates are three to six orders of magnitude slower than reported experimental rates of biotite weathering.     

A review and analysis of silicate mineral dissolution experiments in natural silicate melts

We present a database and a graphical analysis of published experimental results for dissolution rates of olivine, quartz plagioclase, clinopyroxene, orthopyroxene, spinel, and garnet in basaltic and andesitic melts covering a range of experimental temperatures (1100–1500°C) and pressures (10⁵ Pa-3.0 GPa). The published datasets of Donaldson (1985, 1990) and Brearly and Scarfe (1986) are the most complete. Experimental dissolution rates from all datasets are recalculated and normalized to a constant oxygen basis to allow for direct comparison of dissolution rates between different minerals. Dissolution rates (ν) range from 5·10⁻¹⁰ oxygen equivalent moles (o.e.m.) cm⁻² s⁻¹ for olivine in a basaltic melt to 1.3·10⁻⁵ o.e.m. cm⁻² s⁻¹ for garnet in a basaltic melt. Values of ln ν are Arthenian for the experiments examined and activation energies range from 118 to 1800 kJ/o.e.m. for quartz and clinopyroxene, respectively.

The relationship between calculated A/RT for the dissolution reactions, where A is the thermodynamic potential affinity, and values of ν is linear for olivine, plagioclase, and quartz. We interpret this as strong evidence in support of using calculated A as a predictor of ν for, at least, superliquidus melt conditions.     

Formation of vivianite in holocene clay sediments

Concretions of vivianite attaining 0.5 g occur in the upper 4 m of the bottom sediments of the Åsrum lake, Vestfold, S. Norway. The vivianite-bearing zone is of lacustrine origin. Vivianite is absent in the underlying lagunal and marine sediments. Electro-chemical and chemical data are given for the interstitial water and the clay.

Solubility product of vivianite is determined. The equilibrium constant is K = 10⁻⁴⁰ for the reaction 3 Fe⁺⁺ + 2 HPO₄⁻ + 2 OH⁻ = vivianite. The concretions grew by diffusion from a slightly supersaturated interstitial water phase. The degree of supersaturation is calculated. The stability relations of vivianite and hydroxyl apatite are discussed.     

Rock Cavern Stability Analysis Under Different Hydro-Geological Conditions Using the Coupled Hydro-Mechanical Model

Rock cavern stability has a close relationship with the uncertain geological parameters, such as the in situ stress, the joint configurations, and the joint mechanical properties. Therefore, the stability of the rock cavern should be studied with variable geological conditions. In this paper, the coupled hydro-mechanical model, which is under the framework of the discontinuous deformation analysis, is developed to study the underground cavern stability when considering the hydraulic pressure after excavation. Variable geological conditions are taken into account to study their impacts on the seepage rate and the cavern stability, including the in situ stress ratio, joint spacing, and joint dip angle. In addition, the two cases with static hydraulic pressure and without hydraulic pressure are also considered for the comparison. The numerical simulations demonstrate that the coupled approach can capture the cavern behavior better than the other two approaches without the coupling effects.     

An empirical model for estimating hydraulic conductivity of highly disturbed clastic sedimentary rocks in Taiwan

This paper presents the establishment of an empirical HC model for estimating rock mass hydraulic conductivity of highly disturbed clastic sedimentary rocks in Taiwan using high-resolution borehole acoustic televiewer and double packer hydraulic tests. Four geological parameters including rock quality designation (RQD), depth index (DI), gouge content designation (GCD), and lithology permeability index (LPI) were adopted for establishing the empirical HC model. To verify rationality of the proposed HC model, 22 in-situ hydraulic tests were carried out to measure the hydraulic conductivity of the highly disturbed clastic sedimentary rocks in three boreholes at two different locations in Taiwan. Besides, the model verification using another borehole data with four additional in-situ hydraulic tests from similar clastic sedimentary rocks was also conducted to further verify the feasibility of the proposed empirical HC model. The field results indicated that the rock mass in the study area has a conductivity between the order of 10^− 10 m/s and 10^− 6 m/s at the depth between 34 m and 275 m below ground surface. Results demonstrate that the empirical HC model may provide a useful tool to predict hydraulic conductivity of the highly disturbed clastic sedimentary rocks in Taiwan based on measured HC-values.     

Engineering geological characterization of low strength anisotropic rocks in the Himalayan region for assessment of tunnel support

A detailed engineering geological assessment of low strength (6–35 MPa) anisotropic rocks at an ongoing Hydroelectric Project in the Himalayan Region has been carried out. The project (the Nathpa Jhakri Hydroelectric Project) will have one of the world's largest and longest water conducting systems, consisting of a 10.15 m diameter and 27.3 km long head race tunnel, a 942 000 m³ underground desilting complex and a 20 × 49 × 216 m powerhouse cavern in the area. Because these constructions are made in low strength metamorphosed anisotropic rocks like quartz mica schists, biotite schists and muscovite schists, it seemed necessary and worthwhile to undertake a comprehensive study of such rocks. The studies include petrographic and petrofabric analyses of the rocks, geo-mechanical properties, rock stress measurements, rock mass classification using the Q-System and data concerning joint geometry, joint roughness and joint strength. The use of mapped geological and geotechnical data along with the experimental and modelling studies have helped to better understand the behaviour of these rocks in underground openings.     

Oxidative dissolution of pyritic sludge from the Aznalcóllar mine (SW Spain)

As a result of the collapse of a mine tailing dam, a large extension of the Guadiamar valley was covered with a layer of pyritic sludge. Despite the removal of most of the sludge, a small amount remained in the soil, constituting a potential risk of water contamination. The kinetics of the sludge oxidation was studied by means of laboratory flow-through experiments at different pH and oxygen pressures. The sludge is composed mainly of pyrite (76%), together with quartz, gypsum, clays, and sulphides of zinc, copper, and lead. Trace elements, such as arsenic and cadmium, also constitute a potential source of pollution. The sludge is fine grained (median of 12 μm) and exhibits a large surface (BET area of 1.4±0.2 m² g⁻¹).

The dissolution rate law of sludge obtained is r=10^{−6.1(±0.3)} [O₂(aq)]^0.41(±0.04) a_H+^0.09(±0.06) g_sludge m⁻² s⁻¹ (22 °C, pH=2.5–4.7). The dissolution rate law of pyrite obtained is r=10^{−7.8(±0.3)} [O₂(aq)]^0.50(±0.04) a_H+^0.10(±0.08) mol m⁻² s⁻¹ (22 °C, pH=2.5–4.7). Under the same experimental conditions, sphalerite dissolved faster than pyrite but chalcopyrite dissolves at a rate similar to that of pyrite. No clear dependence on pH or oxygen pressure was observed. Only galena dissolution seemed to be promoted by proton activity. Arsenic and antimony were released consistently with sulphate, except at low pH conditions under which they were released faster, suggesting that additional sources other than pyrite such as arsenopyrite could be present in the sludge. Cobalt dissolved congruently with pyrite, but Tl and Cd seemed to be related to galena and sphalerite, respectively.

A mechanism for pyrite dissolution where the rate-limiting step is the surface oxidation of sulphide to sulphate after the adsorption of O₂ onto pyrite surface is proposed.