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.     

Gamma-ray burst neutrino background and star formation history in the universe

We estimate the flux of the gamma-ray burst (GRB) neutrino background and compute the event rate at SK and TITAND in the collapsar model, assuming that GRB formation rate is proportional to the star formation rate. We find that the predicted background neutrino flux is highly sensitive to unknown model parameters, mainly to the mass–accretion rate, to the fraction of disk energy emitted in thermal neutrinos (as opposed to emission through electromagnetic processes), and to the fraction of collapsar events leading to GRBs. The predicted neutrino flux varies over many orders of magnitude as the values of unknown model parameters are varied. We investigate the detection possibility of thermal neutrinos from collapsars which lead to GRBs by TITAND. We find that the GRB neutrino background might be detected by TITAND within 10 yrs only for the optimistic cases in which the average mass–accretion rate is high ( a few M s⁻¹), and the probability that one collapsar generates a GRB is high (f=0.5–1.0).     

Clay Minerals in an Active Hydrothermal Field at Iheya‐North‐Knoll,Okinawa Trough

Iheya‐North‐Knoll is one of the small knolls covered with thick sediments in the Okinawa Trough back‐arc basin. At the east slope of Iheya‐North‐Knoll, nine hydrothermal vents with sulfide mounds are present. The Integrated Ocean Drilling Program (IODP) Expedition 331 studied Iheya‐North‐Knoll in September 2010. The expedition provided us with the opportunity to study clay minerals in deep sediments in Iheya‐North‐Knoll. To reveal characteristics of clay minerals in the deep sediments, samples from the drilling cores at three sites close to the most active hydrothermal vent were analyzed by X‐ray diffraction, scanning electron microscope and transmission electron microscope. The sediments are classified into Layer 0 (shallow), Layer 1 (deep), Layer 2 (deeper) and Layer 3 (deepest) on the basis of the assemblage of clay minerals. Layer 0 contains no clay minerals. Layer 1 contains smectite, kaolinite and illite/smectite mixed‐layer mineral. Layer 2 contains chlorite, corrensite and chlorite/smectite mixed‐layer mineral. Layer 3 is grouped into three sub‐layers, 3A, 3B and 3C; Sub‐layer 3A contains chlorite and illite/smectite mixed‐layer mineral, sub‐layer 3B contains chlorite/smectite and illite/smectite mixed‐layer minerals, and sub‐layer 3C contains chlorite and illite. Large amounts of di‐octahedral clay minerals such as smectite, kaolinite, illite and illite/smectite mixed‐layer mineral are found in Iheya‐North‐Knoll, which is rarely observed in hydrothermal fields in mid‐ocean ridges. Tri‐octahedral clay minerals such as chlorite, corrensite and chlorite/smectite mixed‐layer mineral in Iheya‐North‐Knoll have low Fe/(Fe + Mg) ratios compared with those in mid‐ocean ridges. In conclusion, the characteristics of clay minerals in Iheya‐North‐Knoll differ from those in mid‐ocean ridges; di‐octahedral clay minerals and Fe‐poor tri‐octahedral clay minerals occur in Iheya‐North‐Knoll but not in mid‐ocean ridges.     

New buoy observation system for tsunami and crustal deformation

We have developed a new system for real-time observation of tsunamis and crustal deformation using a seafloor pressure sensor, an array of seafloor transponders and a Precise Point Positioning (PPP ) system on a buoy. The seafloor pressure sensor and the PPP system detect tsunamis, and the pressure sensor and the transponder array measure crustal deformation. The system is designed to be capable of detecting tsunami and vertical crustal deformation of ±8 m with a resolution of less than 5 mm. A noteworthy innovation in our system is its resistance to disturbance by strong ocean currents. Seismogenic zones near Japan lie in areas of strong currents like the Kuroshio, which reaches speeds of approximately 5.5 kt (2.8 m/s) around the Nankai Trough. Our techniques include slack mooring and new acoustic transmission methods using double pulses for sending tsunami data. The slack ratio can be specified for the environment of the deployment location. We can adjust slack ratios, rope lengths, anchor weights and buoy sizes to control the ability of the buoy system to maintain freeboard. The measured pressure data is converted to time difference of a double pulse and this simple method is effective to save battery to transmit data. The time difference of the double pulse has error due to move of the buoy and fluctuation of the seawater environment. We set a wire-end station 1,000 m beneath the buoy to minimize the error. The crustal deformation data is measured by acoustic ranging between the buoy and six transponders on the seafloor. All pressure and crustal deformation data are sent to land station in real-time using iridium communication.     

Broad-band power-law spectra of well-log data in Japan

For the purpose of revealing the statistical characteristics of P -wave velocity, S -wave velocity and density in the uppermost part of the crust, we analysed well-log data obtained from five deep wells in different tectonic regions in Japan: three wells through the mainly sedimentary rocks in the Kanto plain and two wells in the Kuju volcano group in Kyushu Island. In the Kanto plain, the power spectral density of fractional fluctuation of P -wave velocity and that of density are proportional to a power of the spatial wavelength from a few metres to 100 m. where the power index (slope of the power spectral density at double logarithmic scale) is 1.1-1.3. At the Kuju volcano group, that of P - and S -wave velocity and density also obey a power law, with a power index of 1.3-1.6 for wavelengths from a few metres to few hundred metres. Correcting the effect of the moving box-car observation window which corresponds to the separation of two receivers of the logging tool, we find that the power-law characteristics hold for wavelengths down to a few tens of centimetres. The 1-D sections of the elastic inhomogeneities follow a kind of band-limited self-affine random process. Comparing the power spectral densities, we find smaller values of the power index in stable areas and larger values in tectonically active areas. The difference in the power index arises from long-wavelength components.     

Mercury concentration in the ocean

Seventy percent of 342 seawater samples collected in the Bering Sea, North and South Pacific, Japan Sea, East and South China Seas, and Indian Ocean had concentrations of “total” mercury ranging from 3 to 6 ng Hg l^?1 with an arithmetic mean of 5.3 ng l^?1 and a geometric mean of 5.0 ng l^?1. In some cases, a higher concentration was observed at the surface, at the halocline or thermocline, or in the bottom water. But in general, there was no consistent correlation between mercury concentration and depth, except for a statistical tendency for mercury concentration to be slightly higher in the surface water. This tendency suggests that mercury in the ocean is supplied from the atmosphere by rain washout. The latitudinal variation of surface mercury concentrations showed that the maximum concentration at each latitude decreased from 40°N to 30°S. This variation provides evidence that atmospheric mercury is emitted mainly from continental areas naturally or anthropogenically.     

In situ filtration sampler for taking filtered water samples

A new type ofin situ filtration sampler is described. The all Teflon sampler is hung on a nylon rope attached to a steel wire (Fig. 3). By dropping a messenger, the sampler piston is pulled down by a released weight of 50 kg, and the sampler is filled with about 400 ml of filtered water. The sample is not contaminated by metals because there is no metal around the sampler.     

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.     

Forced vibration test of a building with semi‐active damper system

The authors developed a semi‐active hydraulic damper (SHD) and installed it in an actual building in 1998. This was the first application of a semi‐active structural control system that can control a building's response in a large earthquake by continuously changing the device's damping coefficient. A forced vibration test was carried out by an exciter with a maximum force of 100 kN to investigate the building's vibration characteristics and to determine the system's performance. As a result, the primary resonance frequency and the damping ratio of a building that the SHDs were not jointed to, decreased as the exciting force increased due to the influence of non‐linear members such as PC curtain walls. The equivalent damping ratio was estimated by approximating the resonance curves using the steady‐state response of the SDOF bilinear hysteretic system. After the eight SHDs were jointed to the building, the system's performance was identified by a response control test for steady‐state vibration. The elements that composed the semi‐active damper system demonstrated the specified performance and the whole system operated well. Copyright © 2000 John Wiley & Sons, Ltd.