Theory for the semi-analytical calculation of oil recovery and effective relative permeabilities using streamtubes

A semi-analytical method has been developed for calculating oil recovery in two and three dimensions, and for calculating effective relative permeabilities for coarse grids. The calculations are based on the assumption that the effects of a changing mobility field can be accounted for by using fixed streamtube geometries with flowrates updated to account for the changing mobility distribution. The single-phase pressure distribution from a numerical solution of Laplace's equation is used to calculate the pressure distribution for a two-phase flow based on a mapping of the solution of the Buckley-Leverett equation onto the streamtubes derived from the single-phase solution. The displacement calculations for oil recovery are based on theory previously developed by Dykstra and Parsons, extended to include the effects of spatially varying permeability and continuously changing mobilities, as occurs in solutions of the Buckley-Leverett equation for typical values of the mobility ratio. This idea has also been extended to the calculation of effective relative permeabilities for coarse-grid simulation and finally establishes the proper rules for averaging the results of fine-grid numerical simulations of two-phase flow for the definition of effective two-phase flow properties on coarse grids. These calculations have been generalized to three-dimensional flows by the simple device of conceptually inserting a gridded plane across the flow and defining each streamtube at that location as those streamlines which pass through any one of the grid cells. When combined with time-of-flight calculations from the gridded plane to both the producer and injector, the distribution of pore volume along each streamtube can be calculated. This information, combined with a tabulation of the single-phase, steady-state pressure distribution along each streamtube, provides all of the information needed for the semi-analytical calculation of oil recovery and effective flow properties in three-dimensional flows. © 1997 Elsevier Science Ltd. All rights reserved     

Channel bed adjustments following major aggradation in a steep headwater setting: findings from Oyabu Creek, Kyushu, Japan

A typhoon in 1993 induced major aggradation along Oyabu Creek, a steep, gravel bed mountain stream in Kyushu, Japan. Processes of sediment reworking are inferred from a 7-year monitoring program that measured adjustments to channel cross-sections, the longitudinal profile, and the extent/distribution of bedrock outcrops along a 3-km study reach. Over time, the reach adopted a riffle and pool structure, with notable increase in the area of exposed bedrock on the bed. This adjustment process was characterised by progressive reduction in sediment storage change per unit flow. The relaxation pathway following disturbance induced by the typhoon was shaped by the magnitude and frequency of subsequent rainfall events, the capacity of these events to transport available sediments, and physical linkages between reaches. Adjacent subreaches demonstrated differing relaxation pathways in response to these influences, induced by spatial and temporal variability in threshold conditions along the channel. Longer-term evidence indicates that responses to major disturbance, such as the 1993 typhoon, occur as ‘cycles’ of around 20-year duration. A relaxation period of 7 years is required to attain a quasi-equilibrium bed configuration and rate of sediment flux. The timeframe of cycles is considered to reflect changes to hillslope–channel bed coupling, marking the period required to generate sufficient sediment stores to reactivate phases of aggradation and subsequent degradation.     

Estimation of temporally averaged sediment delivery ratio using aggradational terraces in headwater catchments of the Waipaoa River,North Island,New Zealand

The sediment delivery ratio was estimated for two periods (28 years and eight years) following reforestation of seven tributary catchments (0·33 to 0·49 km²) in the headwaters of the Waipaoa River basin, North Island, New Zealand. In these catchments, gully erosion, which largely resulted from clearance of the natural forest between 1880 and 1920, is the main source of sediment to streams. Reforestation commenced in the early 1960s in an attempt to stabilize hillslopes and reduce sediment supply. Efforts have been partially successful and channels are now degrading, though gully erosion continues to supply sediment at accelerated rates in parts of the catchment. Data from the area indicate that the sediment delivery ratio (SDR) can be estimated as a function of two variables, ψ (the product of catchment area and channel slope) and A _g (the temporally averaged gully area for the period). Sediment input from gullies was determined from a well defined relationship between sediment yield and gully area. Sediment scoured from channels was estimated from dated terrace remnants and the current channel bed. Terrace remnants represent aggradation during major floods. This technique provides estimates of SDR averaged over periods between large magnitude terrace‐forming events and with the present channel bed. The technique averages out short‐term variability in sediment flux. Comparison of gully area and sediment transport between two periods (1960–1988 and 1988–1996) indicates that the annual rate of sediment yield from gullies for the later period has decreased by 77 per cent, sediment scouring in channels has increased by 124 per cent, and sediment delivered from catchments has decreased by 78 per cent. However, average SDR for the tributaries was found to be not significantly different between these periods. This may reflect the small number of catchments examined. It is also due to the fact that the volume of sediment scoured from channels was very small relative to that produced by gullies. According to the equation for SDR determined for the Waipaoa headwaters, SDR increases with increasing catchment area in the case where A _g and channel slope are fixed. This is because the amount of sediment produced from a channel by scouring increases with increasing catchment area. However, this relationship does not hold for the main stem of the study catchments, because sediment delivered from its tributaries still continues to accumulate in the channel. Higher order channels are, in effect, at a different stage in the aggradation/degradation cycle and it will take some time until a main channel reflects the effects of reforestation and its bed adjusts to net degradation. Results demonstrate significant differences among even low order catchments, and such differences will need to be taken into consideration when using SDR to estimate sediment yields. Copyright © 2001 John Wiley & Sons, Ltd.     

Analysis and control of vibration of ropes in a high-rise elevator under earthquake excitation

In this research, the vibration of elevator ropes, including the main rope and compensation rope are investigated simultaneously in a high-rise elevator system under earthquake excitation. Moreover, the paper presents a new control method to restrain the sway of both ropes. This study considers varying rope lengths during elevator operation which cause other system parameters such as natural frequency, and damping ratio to be time-variant variables. The dynamics of the ropes are analyzed by solving the governing non-stationary, nonlinear equation numerically. In order to mitigate the vibration of ropes in several motion conditions, particularly upwards movement, downward movement, stopped at the lowest position, and stopped at the highest position, an active equipment is installed at the compensation sheave. The stability of the system using the controller is analyzed at four states: without disturbance and static car, without disturbance and mobile car, including disturbance and static car, and including disturbance and mobile car. The efficiency of the controller used for dampening the vibration of elevator ropes is validated by numerical simulation results.     

Influence of submarine fumaroles on the distribution of mercury in the sediment of Kagoshima Bay,Japan

To estimate the influence of mercury emitted from submarine fumaroles, the horizontal and vertical distribution of mercury in sediment of Kagoshima Bay was studied. The fumaroles are located in the northern bay head area, and the sediment samples had been taken from 52 points throughout the bay with a gravity core sampler. The core samples obtained were cut at a thickness of 1–2 cm and used for measurements. The total concentration of mercury in surface sediment in the northern and central areas of the bay was 51–679 μg kg^− 1 (average 199 μg kg^− 1, n = 22) and 23–100 μg kg^− 1 (average 55 μg kg^− 1, n = 30), respectively. The highest value was obtained in the vicinity of the fumaroles. The mercury concentration in sediment near the fumaroles varied with depth, which may reflect the variation in fumarolic activity. A successive extraction method was applied to the speciation of mercury in the sediment. The results showed that sediment taken in the vicinity of submarine fumaroles contained a higher percentage of mercury bound with organic matter.     

Gully systems under undisturbed indigenous forest, East Coast Region, New Zealand

Gully erosion is commonly associated with agricultural landscapes where vegetation clearance and farming practices increase runoff, leading to fluvial incision. However, gully erosion can also occur in forests that have undergone some form of disturbance, either natural or resulting from human impacts. This paper reports on recent gully development within areas of undisturbed indigenous forest as a result of a high magnitude rainfall event on the East Coast of New Zealand's North Island. This region, through a combination of crushed and sheared rock types, steep topography, and tectonic and climatic setting, has high natural rates of erosion, exacerbated by European deforestation in the late 19th and early 20th centuries.Sequential air photographs, spanning a 58 year period between 1939 and 1997 were used to classify and document the growth and recovery of gully systems in the 14.1 km² headwaters of the Mangaoporo catchment. Following a severe cyclone in 1988, with a rainfall of 535 mm, there were 21 active gully systems within the indigenous forest. On photography prior to 1988 only four gully systems were present. During this period there were 8 major rainfall events (150–250 mm). Despite further 5 rainfall events of 150–250 mm between 1988 and 1997 all gully systems showed signs of recovery, with a combined reduction in active area of 37%. The nature and location of these features is strongly influenced by lithology (orientation of jointing and bedding), and to a topographic threshold defined by catchment slope and catchment area.     

Characterization of slip zone associated with the 1999 Taiwan Chi-Chi earthquake: X-ray CT image analyses and microstructural observations of the Taiwan Chelungpu fault

To characterize the fault-related rocks within the Chelungpu fault, we performed X-ray computed tomography (CT) image analyses and microstructural observations of Hole B core samples from the Taiwan Chelungpu-fault Drilling Project. We identified the slip zone associated with the 1999 Chi-Chi earthquake, within the black gouge zone in the shallowest major fault zone, by comparison with previous reports. The slip zone was characterized by low CT number, cataclastic (or ultracataclastic) texture, and high possibility to have experienced a mechanically fluidized state. Taking these characteristics and previous reports of frictional heating in the slip zone into consideration, we suggested that thermal pressurization was the most likely dynamic weakening mechanism during the earthquake.     

Two modes of climatic control in the Holocene stalagmite record from the Japan Sea side of the Japanese islands

The Holocene stalagmite FG01 collected at the Fukugaguchi Cave in Itoigawa, central Japan provides a unique high‐resolution record of the East Asian winter monsoon. Because of the climate conditions on the Japan Sea side of the Japanese islands, the volume of precipitation during the winter is strongly reflected in the stalagmite δ¹⁸O signal. Examination of the carbon isotopes and the Mg/Ca ratio of FG01 provided additional information on the Holocene climate in Itoigawa, which is characterized by two different modes separated at 6.4 ka. Dripwater composition and the correlation between the δ¹³C and Mg/Ca data of FG01 indicate the importance of prior calcite precipitation (PCP), a process that selectively eliminated ¹²C and calcium ions from infiltrating water from CO₂ degassing and calcite precipitation. In an earlier period (10.0–6.4 ka), an increase in soil pCO₂ associated with warming and wetting climate trends was a critical factor that enhanced PCP, and resulted in an increasing trend in the Mg/Ca and δ¹³C data and a negative correlation between the δ¹³C and δ¹⁸O profiles. A distinct peak in the δ¹³C age profile at 6.8 ka could be a response to an increase of approximately 10% in C4 plants in the recharge area. At 6.4 ka, the climate mode changed to another, and correlation between δ¹⁸O and δ¹³C became positive. In addition, a millennial‐scale variation in δ¹⁸O and pulsed changes in δ¹³C and Mg/Ca became distinct. Assuming that δ¹⁸O and PCP were controlled by moisture in the later period, the volume of precipitation was high during 6.0–5.2, 4.4–4.0, and 3.0–2.0 ka. In contrast, the driest interval in Itoigawa was during 0.2–0.4 ka, and broadly corresponds to the Little Ice Age.     

Isotopic geochemistry of Hualalai shield-stage tholeiitic basalts from submarine North Kona region,Hawaii

The geochemical data of Hualalai tholeiitic basalts allow extension of the temporal variations established at Mauna Loa back in time, and provide important information for the long-term temporal variation of the Hawaiian lavas. We report new Hf, Pb, Nd, and Sr isotope compositions for 32 Hualalai tholeiitic basalts collected from deep submarine portions of the North Kona region. The samples were collected from the lower section of the North Kona bench (dives K218 and K219), a submarine stratigraphic section at Hualalai volcano's northwest rift zone (dive S690), and an elongate ridge outboard of the central section of the bench (dive S692), during two JAMSTEC Hawaii cruises in 2001 and 2002. The Hualalai shield-stage tholeiitic basalts have magma source isotopic signatures similar to Mauna Loa. The new data shows temporal Pb and Sr isotope trends that correspond to the long-term temporal variations in Loa-trend lavas, and the Hualalai–Mauna Loa lavas seem to show inter-shield geochemical excursions. Variation in Pb and Sr isotopes at Hualalai appears to take place over a longer time scale than at Mauna Loa. The merged Hualalai–Mauna Loa isotopic trends support models where heterogeneous material in the plume conduit is distributed chaotically, with variable cross-sectional density and length scale.