X-ray CT based numerical analysis of fracture flow for core samples under various confining pressures

The X-ray CT based numerical analysis of fracture flow for core samples, recently developed by the authors, was applied to two granite core samples having either a mated artificial or a mated natural fracture at confining pressures of 5 to 50 MPa. A third-generation medical X-ray CT scanner was used to image the samples within a core holder consisting of an aluminum liner and a carbon fiber overwrap. Fracture models (i.e., aperture distributions) were obtained by the CT images, the resolution of which was coarser than the apertures, and a single-phase flow simulation was performed using a local cubic law-based fracture flow model. Numerical results were evaluated by a fracture porosity measurement and a solution displacement experiment using NaCl and NaI aqueous solutions. These numerical results coincided only qualitatively with the experimental results, primarily due to image noise from the aluminum liner of the core holder. Nevertheless, the numerical results revealed flow paths within the fractures and their changes with confining pressure, whereas the experimental results did not provide such results. Different stress-dependencies in the flow paths were observed between the two samples despite the similar stress-dependency in fracture porosity and permeability. The changes in total area of the flow paths with confining pressure coincided qualitatively with changes in breakthrough points in the solution displacement experiment. Although the data is limited, the results of the present study suggest the importance of analyzing fluid flows within naturally fractured core samples under in situ conditions in order to better understand the fracture flow characteristics in a specific field. As demonstrated herein, X-ray CT-based numerical analysis is effective for addressing this concern. Using a multi-phase flow model, as well as a core holder constructed of an engineered plastic, should provide a useful, non-destructive, and non-contaminative X-ray CT-based fracture flow analysis for core samples under in situ conditions in future studies.     

Culture-dependent and independent analyses of subsurface microbial communities in oil-bearing strata of the Sagara oil reservoir

Abstract   Culture-dependent and independent methods were used to evaluate the microbial communities in cores collected at depths up to 200 m in oil-bearing and oil-free strata near the Sagara oil reservoir near Shizuoka, Japan. Direct microbial counts revealed much higher numbers (2.2 × 10⁴−7.9 × 10⁶/g) of microbes in the oil-bearing strata than in the oil-free zones, where counts were uniform at approximately 1.0 × 10⁴/g. Molecular taxonomic analyses via 16S rRNA gene sequence comparisons showed that the oil-free strata were dominated by members of the γ-Proteobacteria including Pseudomonas , Stenotrophomonas and Sphingomonas , whereas the oil-bearing strata were dominated by a single species closely related to Pseudomonas stutzeri . All archaeal clones were phylogenetically affiliated with the uncultured soil group in Crenarchaeota with the exception of a single phylotype that belonged to the genus Thermococcus . Culture-dependent analysis was carried out by most-probable-number culturing as well as direct plating to determine viable cell counts, using both complex organic substrates or native oil and autotrophic media. Both culture-dependent and independent methods revealed the abundant cultivable member was the aerobic oil-degrading bacterium Pseudomonas stutzeri and neither autotrophs nor anaerobic heterotrophs could be detected in the oil-bearing strata.     

A structure of the Kuroshio and its related upwelling on the East China Sea shelf slope

The ADCP on an advanced towed fish with controllable main and tail wings, called DRAKE measured a detailed sectional structure of the Kuroshio flowing to the NE along the East China Sea shelf slope west of Okinawa. At the observation period, a countercurrent directed to the SW formed in near-bottom water on the shelf slope. The horizontal flow perpendicular to the stream axis of the Kuroshio constructed a convergence zone around the boundary between the Kuroshio and the countercurrent. An intensive upwelling with the maximum velocity of 2.8 cm s^–1 was found to distribute on the shelf slope around the convergence zone. A dynamic cause of this intensive upwelling is discussed carefully.     

Plutonium diffusivity in compacted bentonite

Measurements on plutonium diffusivity in water-saturated compacted bentonite were carried out. Representative specimens of sodium bentonite were taken from the Tsukinuno and Kuroishi mines situated in northeast Japan. Tsukinuno bentonite was divided into three types: raw type, purified Na-type, and H-type which was prepared by treating Na-type bentonite with hydrochloric acid. Kuroishi bentonite contained chlorite as impurity. H-type bentonite was used as reference for the convenience of profile measurement in bentonite, since plutonium diffusivity in H-type bentonite was considered to be larger than that in Na-type bentonite because of low pH and low swelling pressure of H-type bentonite.

Sampled bentonite was compacted into pellets of 20 mm in diameter and 20 mm in height. Bulk densities of these specimens were 1200–1800 kg/m³ for purified Na-type and H-type bentonite and 1600 kg/m³ for raw type bentonite.

Plutonium profiles obtained in H-type bentonite can be explained by diffusion equation with constant concentration source. Diffusivity ranges from 10^-13 to 10¹² m²/s for H-type and Kuroishi impure sodium bentonite. Diffusivity in both raw type and purified Tsukinuno bentonite was was estimated to less than 10^-14 M²/s. Diffusivity in H-type bentonite showed a tendency to decrease with increasing density. Influence of in bentonite was also studied. Quartz content up to 50% or hematite content up to 1% did not influence diffusivity significantly in H-type bentonite.

The chemical species of plutonium in pore water of Na-type and H-type were estimated Pu(OH)₃^-; and PuO₂^- , respectively.     

Spatial delineation of riparian groundwater within alluvium deposit of mountainous region using Laplace equation

The interaction between surface and groundwater plays a key role in a riparian ecosystem while the size of riparian groundwater has not been typically incorporated into hydrological modelling systems. An extensive geophysical survey composed of 25 individual DC electrical resistivity profiles was conducted at the Blair–Wallis site in Wyoming. The observed resistivity images show a near‐surface aquifer interpreted as the saturated alluvium deposit along the channel, rather than the geological bedrock. Based on the electrical resistivity images, it can be inferred that only the near‐surface portion of the groundwater actively interacts with the stream flow in the mountainous and hilly watershed. This study attempted the spatial extrapolation of the measured riparian aquifer depths by means of fitting functions based on the surface topography. The analysis indicated that the boundary of the riparian aquifer well corresponds to the topographical inflexion point of the hill slope profile. It was also demonstrated that the extent of alluvium deposit, where the area of riparian aquifer is indicated, can be delineated using the slope and curvature maps in the geographic information system. Then, the parabolic and biharmonic functions were tested for the groundwater depth estimation using the developed alluvium deposit map. The proposed methodology was effective if geological diffusion processes by wind and water dominated the topography. The spatial map of the active aquifer will be useful in hydrological drought analysis because it is considered to be a main source of baseflow during dry seasons.     

Resolution Improvement of Solar Images

A method for the improvement of resolution in an observed solar image is proposed. A blind deconvolution method is used for restoration of an atmospherically-degraded solar image, and a super-resolution method is applied to its restored image to improve the resolution. It is confirmed that a blind deconvolution process can restore fine structures that are blurred in an observed image, and that the super-resolution process can make a cutoff frequency in a blind-deconvolved image higher. A time series of super-resolved images of a sunspot observed with the 70-cm Vacuum Tower Telescope at Teide Observatory is presented.     

Reciprocal sound transmission measurement of mean current and temperature variations in the central part (Aki-nada) of the Seto Inland Sea,Japan

A 30 km-range reciprocal sound transmission experiment was carried out on the line connecting Honshu and Shikoku (the first and fourth biggest main Japanese islands, respectively) in the central part (Aki-nada) of the Seto Inland Sea, Japan, during March–May 2010 to measure the mean current and temperature variations over the sea. The range-averaged current along the sound transmission line was estimated to have a mean and standard deviation of (3.8–4.4) ± (1.7–1.8) cm/s after converting the travel time difference data into currents, including a fortnightly tidal variation in the range of ±30 cm/s. The positive mean current implies slow water movement from the west to east through Aki-nada. The range-averaged speed of sound was estimated by converting from the mean travel time or one-way travel time into the speed of sound, and further converted into temperature for fixed values of salinity and depth, according to the standard speed of sound formula. Besides the precise measurement (to an accuracy of 0.01°C) of semidiurnal and diurnal tidal variations and seasonal warming, the temperature data showed periodic variations with periods of 7.0 and 21.1 days that had never been observed in Aki-nada before. This study suggests that reciprocal sound transmission is a powerful technique for the long-term accurate measurement of mean current and temperature variations in coastal and inland seas.     

Optimum design of the ocean acoustic tomography system for the Sea of Japan

An ocean acoustic tomography system covering the region of 800×1000 km with the spatial resolution of eddy-resolving scales has been designed on the basis of computer experiments using the hydrographic data collected in the Sea of Japan. The optimum number of acoustic sources required for 20 acoustic receivers was estimated as 13 by changing the source number. The spatial resolution for the optimum system was 41 km smaller than the dominant size of meso-scale eddies in the Sea of Japan. The effect of travel-time errors on tomographic maps is also quantified.     

A coastal acoustic tomography experiment in the Tokyo Bay

Eight sets of coastal acoustic tomography (CAT) systems were deployed during November 29 to December 10, 2002 at the coasts on both sides of Tokyo Bay to measure tidal current structures at 15-rain interval.Sound transmission across the Tokyo Bay (between Yokohama and Chiba)was successfillly traced,even under severe interference from ship generated wakes and bubbles.Tidal current fields changing from northward to southward flow are well reconstructed by the inverse analysis of travel-time difference data for a period with the best sound transmission condition. It is suggested that the CAT is the most powerful tool to continuously map tidal current fields in the coastal seas with heavy shipping traffic and fisheries activity.     

Monitoring Surface Velocity from Repeated ADCP Observations and Satellite Altimetry

A method has been developed to monitor the surface velocity field by combining repeated acoustic Doppler current profiler (ADCP) observations and satellite altimetry data. The geostrophic velocity anomaly is calculated from the sea surface height anomaly field estimated from the altimetry data by an optimal interpolation. It has been confirmed that this accurately observes the smoothed velocity anomaly field when the interpolation scales are set according to the spatio-temporal sampling pattern of the altimeter used. The velocity anomaly obtained from the altimetry data is subtracted from the repeated ADCP observations to estimate temporal mean velocity along the ship tracks. Regularly sampled, nine-year time series of surface velocity can then be obtained by adding the computed mean velocity and the altimetry anomaly components. This clearly illustrates surface velocity fluctuations such as the movement of the Kuroshio axis due to its meandering and an increase of the interannual variability of the Subtropical Countercurrent toward its downstream region. This revised version was published online in July 2006 with corrections to the Cover Date.