A novel finite element double porosity model for multiphase flow through deformable fractured porous media

Based on the theory of double-porosity, a novel mathematical model for multiphase fluid flow in a deforming fractured reservoir is developed. The present formulation, consisting of both the equilibrium and continuity equations, accounts for the significant influence of coupling between fluid flow and solid deformation, usually ignored in the reservoir simulation literature. A Galerkin-based finite element method is applied to discretize the governing equations both in the space and time domain. Throughout the derived set of equations the solid displacements as well as the fluid pressure values are considered as the primary unknowns and may be used to determine other reservoir parameters such as stresses, saturations, etc. The final set of equations represents a highly non-linear system as the elements of the coefficient matrices are updated during each iteration in terms of the independent variables. The model is employed to solve a field scale example where the results are compared to those of ten other uncoupled models. The results illustrate a significantly different behaviour for the case of a reservoir where the impact of coupling is also considered. © 1997 by John Wiley & Sons, Ltd.     

Timing of Late Quaternary glaciations in the Himalayas of northern Pakistan

Optically stimulated luminescence dating of Late Quaternary glaciogenic sediments was undertaken in critical areas of the Himalayas of northern Pakistan in order to examine the timing of glaciation. The dates demonstrate that several glaciations occurred during the last glacial cycle. In Swat, the Grabral 2 Stade and the Kalam I Stade were dated at ca. 77 ka and ca. 38 ka, respectively. The error on the former date is large and it is conceivable that the moraines may have formed during the early part of Oxygen Isotope Stage 3 rather than during Oxygen Isotope Stage 4. The Kalam I Stade, however, clearly represents a glaciation during Oxygen Isotope Stage 3. The oldest moraines and those at the lowest altitude in the Indus valley at Shatial have an age of ca. 60 ka. These also relate to a major glacial advance during Oxygen Isotope Stage 3. A younger series of moraines, the Jalipur Tillite, and glaciofluvial sands at Liachar in the Indus valley, and moraines at Rampur–Tarshing have ages of ca. 27 ka, ca. 21–23 ka and ca. 15 ka, respectively. These dates show that glaciers also occupied parts of the Indus valley during Oxygen Isotope Stage 2. These dates and the morphostratigraphy show that glaciation in the Pakistani Himalaya was more extensive during the early part of the last glacial cycle and that the local last glacial maximum in Pakistan was asynchronous with the maximum extent of Northern Hemisphere ice sheets. Copyright © 2000 John Wiley & Sons, Ltd.     

Preliminary soil micromorphology studies of landscape and occupation history at Tabon Cave,Palawan, Philippines

Preliminary soil micromorphology study of cultural sediments at Tabon Cave, Philippines, supports interpretations of sporadic occupation in the Paleolithic. The presence in some deposits of authigenic minerals potentially related to altered cultural materials, such as ash, needs further investigation. Later in the sequence there is a marked change in local depositional processes, with the onset of significant quartz sand deposition in layers dating from the Middle to Late Holocene. This could relate to beach development in the area. Future sedimentological study and dating would confirm this interpretation, which suggests that in appropriate settings, stratified cave sediments could be useful for the study of regional sea level rise. © 2007 Wiley Periodicals, Inc.     

Plateau reduction by drainage divide migration in the Eastern Cordillera of Colombia defined by morphometry and 10Be terrestrial cosmogenic nuclides

Catchment‐wide erosion rates were defined using ¹⁰Be terrestrial cosmogenic nuclides for the Eastern Cordillera of the Colombian Andes to help determine the nature of drainage development and landscape evolution. The Eastern Cordillera, characterized by a smooth axial plateau bordered by steep flanks, has a mean erosion rate of 11 ± 1 mm/ka across the plateau and 70 ± 10 mm/ka on its flanks, with local high rates >400 mm/ka. The erosional contrast between the plateau and its flanks was produced by the increase in the orogen regional slope, derived from the progressive shortening and thickening of the Eastern Cordillera. The erosion rates together with digital topographic analysis show that the drainage network is dynamic and confirms the view that drainage divides in the Eastern Cordillera are migrating towards the interior of the mountain belt resulting in progressive drainage reorganization from longitudinal to transverse‐dominated rivers and areal reduction of the Sabana de Bogotá plateau. Copyright © 2016 John Wiley & Sons, Ltd.     

Seismic vulnerability scenarios of Unreinforced Masonry churches in New Zealand

A vulnerability analysis of c.300 unreinforced Masonry churches in New Zealand is presented. The analysis uses a recently developed vulnerability index method (Cattari et al. in Proceedings of the New Zealand Society for Earthquake Engineering NZSEE 2015 conference, Rotorua, New Zealand, 2015a; b; SECED 2015 conference: earthquake risk and engineering towards a Resilient World, Cambridge; Goded et al. in Vulnerability analysis of unreinforced masonry churches (EQC 14/660)—final report, 2016; Lagomarsino et al. in Bull Earthq Eng, 2018), specifically designed for New Zealand churches, based on a widely tested approach for European historical buildings. It consists of a macroseismic approach where the seismic hazard is defined by the intensity and correlated to post seismic damage. The many differences in typologies of New Zealand and European churches, with very simple architectural designs and a majority of one nave churches in New Zealand, justified the need to develop a method specifically created for this country. A statistical analysis of the churches damaged during the 2010–2011 Canterbury earthquake sequence was previously carried out to develop the vulnerability index modifiers for New Zealand churches. This new method has been applied to generate seismic scenarios for each church, based on the most likely seismic event for 500 years return period, using the latest version of New Zealand’s National Seismic Hazard Model. Results show that highly vulnerable churches (e.g. stone churches and/or with a weak structural design) tend to produce higher expected damage even if the intensity level is lower than for less vulnerable churches in areas with slightly higher seismicity. The results of this paper provide a preliminary tool to identify buildings requiring in depth structural analyses. This paper is considered as a first step towards a vulnerability analysis of all the historical buildings in the country, in order to preserve New Zealand’s cultural and historical heritage.     

Variability of the useful life of reservoirs in tropical locations: A case study from the Burdekin Falls Dam,Australia

Large dam construction for irrigation, hydropower, water supply, and flood control in tropical to sub-tropical areas increased markedly after 1950. Many of the dams built during this period have filled with sediment and no longer perform to their original design specifications; in many cases, forecast dam life was greatly overestimated. This study investigates the useful life of Burdekin Falls Dam (BFD), northeastern Australia and compares the findings with other tropical reservoirs. Using two independent methods it is estimated that between 61 and 65 million m³ of sediment has been deposited in the reservoir over the 24 years of operation through 2011. This sediment volume equates to an average of 0.15% of capacity lost per year since construction was completed. If current sediment loads/climate regimes persist, reservoir capacity will be reduced by 50% after 345 years. However, the useful life of the BFD reduces to just 276 years when drawdown data are considered; these data show reservoir use would be affected once 40% of storage was filled with sediment, with a 60% drawdown return period of 1 in 15 years. When compared to similar large tropical to sub-tropical reservoirs, the BFD has a slightly longer reservoir useful life than dams in India and a much longer half-life than for both similar-sized and larger dams in China, Brazil, and Iran. Properties of the BFD that promote a longer useful life include a lower trap efficiency, relatively low annual sediment load delivered to the reservoir, limited sediment deposition behind the dam wall (and uniform distribution of deposited sediment), and the export of highly turbid annual floodwaters before settling and deposition of any remaining sediment within the reservoir.