An application of Karmarkar's interior-point linear programming algorithm for multi-reservoir operations optimization

Optimization of multi-reservoir systems operations is typically a very large scale optimization problem. The following are the three types of optimization problems solved using linear programming (LP): (i) deterministic optimization for multiple periods involving fine stage intervals, for example, from an hour to a week (ii) implicit stochastic optimization using multiple years of inflow data, and (iii) explicit stochastic optimization using probability distributions of inflow data. Until recently, the revised simplex method has been the most efficient solution method available for solving large scale LP problems. In this paper, we show that an implementation of the Karmarkar's interior-point LP algorithm with a newly developed stopping criterion solves optimization problems of large multi-reservoir operations more efficiently than the simplex method. For example, using a Micro VAX II minicomputer, a 40 year, monthly stage, two-reservoir system optimization problem is solved 7.8 times faster than the advanced simplex code in MINOS 5.0. The advantage of this method is expected to be greater as the size of the problem grows from two reservoirs to multiples of reservoirs. This paper presents the details of the implementation and testing and in addition, some other features of the Karmarkar's algorithm which makes it a valuable optimization tool are illuminated.     

塔里木河流域径流变化趋势及其对气候变化的响应

This paper has studied the change of streamflow and the impact of climatic variability conditions on regional hydrological cycle in the headwater of the Tarim River Basin. This study investigates possible causes of observed trends in streamflow in an environment which is highly variable in terms of atmospheric conditions, and where snow and ice melt play an important role in the natural hydrological regime. The discharge trends of three head streams have a significant increase trend from 1957 to 2002 with the Mann-Kendall test. Complex time-frequency distributions in the streamflow regime are demonstrated especially by Morlet wavelet analysis over 40 years. The purpose is to ascertain the nature of climatic factors spatial and temporal distribution, involved the use of EOF （Empirical Orthogonal Function） to compare the dominant temperature, precipitation and evaporation patterns from normally climatic records over the Tarim＇s headwater basin. It shows that the first principal component was dominated since the 1990s for temperature and precipitation, which identifies the significant ascending trend of spatial and temporal pattern characteristics under the condition of the global warming. An exponential correlation is highlighted between surface air temperature and mean river discharge monthly, so the regional runoff increases by 10%-16% when surface air temperature rises by 1 ℃. Results suggest that headwater basins are the most vulnerable environments from the point of view of climate change, because their watershed properties promote runoff feeding by glacier and snow melt water and their fundamental vulnerability to temperature changes affects rainfall, snowfall, and glacier and ice melt.     

Reconstructing hydrological variability from testate amoebae analysis in Carpathian peatlands

Peatlands offer the potential for high resolution records of water balance over Holocene timescales, yet this potential is under-exploited in many areas of the world. Within Europe, peatlands are mostly confined to areas north of 55° N, but several areas of southern and eastern Europe contain small peatlands which may be suitable for palaeoclimatic reconstruction. In this paper we test the potential of peatlands in the Carpathian region for deriving quantified estimates of water table changes using testate amoebae analysis. A training set for palaeohydrological reconstruction from testate amoebae assemblages was obtained by collecting surface samples from 13 peatlands, including 9 from Hungary and 4 from Transylvania (Romania). Using a simple measure of mean annual water tables estimated from staining of PVC tape, we found that some peatlands were heavily influenced by runoff and groundwater, and were therefore not suitable as modern analogues of ombrotrophic climatically sensitive sites. The relationship between the testate amoebae assemblages in the modern samples and the environmental variables was explored using CCA. The CCA biplot showed that the most important variables are depth to water table and moisture content, confirming that hydrology is a key control on taxon distribution. pH was a secondary gradient. A transfer function for % moisture and depth to water table was established and applied to fossil assemblages from a sequence from Fenyves-tető, Transylvania, Romania. The reconstructed water table shows a number of variations which have parallels with other palaeoclimatic records from Europe and the North Atlantic prominent phases of higher water tables are associated with the periods 8000–8300 cal BP, 3000–2500 cal BP and after 600 cal BP. We suggest that these were periods of particular intensification of westerly airflow which affected eastern Europe as well as western and central Europe.     

Sedimentological evolution of the Lower Series formation in the southern area of the Hassi R’Mel field,Saharan Platform,Algeria

The Hassi R’Mel field, which is the largest gas field in Algeria, produces principally gas condensate. The produced gas and oil comes from four reservoirs of Triassic age (reservoir A, reservoir B, reservoir C, and the Lower Series). All the previous studies in this field focused only on the upper reservoirs (A, B, and C), which have shown an excellent petrophysical properties. However, due to the presence of shale facies as well as volcanic intrusions, the Lower Series reservoir has been less explored. Thus, our study is mainly focused on this Lower Series in order to find the best sandy levels. The aim of this paper was to determine the facies types in the area, their distribution in the southern area of the Hassi R’Mel gas field, and their depositional environment. The Triassic rocks in this field are divided into four formations. From the youngest to oldest, these are S4, Lower Shale, Horizon A, and the Lower Series. The Lower Series formation can be further divided into three members (shaly sandstone, andesite, and lower shale). The shaly sandstone member within the Lower Series is composed of four units that are present all over the study area with an important pinch-out toward the northern part. The depositional environment within the defined units is a braided fluvial system in the southern part that evolves into a meandering fluvial system toward the north. Thus, the study area is considered to be a transitional zone between two depositional systems.     

Volcanic emissions and the early Earth atmosphere

Despite uncertainties in our understanding of early Earth volcanism and atmospheric composition, thermodynamic modelling is able to offer estimates of the global production of reactive trace species (NO, OH, SO₃, Cl, Br and I) from early Earth volcanism, and thereby to shed light on processes which may have been different in Earth’s early atmosphere. Model results show that thermal decomposition of magmatic H₂O, CO₂ and SO₂ in high-T mixtures of magmatic and atmospheric gases (at T > 1400 °C) generate high levels of reactive trace gas species. Production of these reactive trace species is insensitive to atmospheric CO₂ in mixtures where the atmospheric gas is the minor component and will hence continue during periods of low atmospheric CO₂. Fluxes of NO, OH, Cl, Br and I from early Earth volcanism are predicted to exceed those from modern Earth volcanism as the higher temperature of early Earth emissions compensates for lower levels of O₂ in the atmosphere, compared to the modern Earth. Under certain conditions, the volcanic NO flux from early Earth volcanism is found to be comparable to other sources of reactive N such as lightning NO and photochemical HCN. This is one possible source of fixed nitrogen which may alleviate any postulated Archean nitrogen crisis. Our thermodynamic model reveals that production of SO₃ (a potential precursor for near-source volcanic sulphate and hence ‘primary’ volcanic aerosol) is likely to be significantly lower from early Earth volcanism. Uncertainty in the pathway to near-source sulphate in modern volcanism (i.e., the reaction of SO₃ with water or direct emission) introduces a large uncertainty into the production rate of near-source volcanic sulphate on the early Earth.     

Sr isotopic composition of Paleoproterozoic C-rich carbonate rocks: The Tulomozero Formation,SE Fennoscandian Shield

Rb–Sr systematics has been studied in ¹³C-rich carbonate rocks of the Paleoproterozoic (2.09 ± 0.07 Ga) Tulomozero Formation in the northern Onega Lake area, the SE Fennoscandian Shield. The formation is divided into eight members (A–F) consisting of greenschist-facies-grade, variegated sandstones, siltstones, mudstones, stromatolitic dolostones and subordinate crystalline limestones. Samples of carbonate rocks were obtained from two overlapping drillholes intersecting the entire thickness of the Tulomozero Formation. Prior to isotope analysis, the rocks powders were treated with 1N ammonium acetate for partial removal of the late epigenetic carbonate phases. Major resetting of the Rb–Sr systems in the Tulomozero carbonate rocks appears to take place during the Svecofennian regional metamorphic event, and it was screened by using Mn/Sr, Fe/Sr, Mg/Ca, and ¹⁸O/¹⁶O ratios. High Sr content (up to 2080 μg/g in limestones, and 530 μg/g in dolostones) coupled with low Fe/Mn (<0.40) ratios in the Tulomozero carbonate rocks of Members A, B (the lower part), D, F, and E are consistent with accumulation of original carbonate sediments in evaporitic lacustrine, playa, and sabkha environments. A decrease in the Sr content with concurrent increase in the Fe/Mn ratio (>0.40) in dolostones of the upper part of Member B, and of Members G and H is indicative of seawater influxes (sea transgression) into the Tulomozero basin. The ⁸⁷Sr/⁸⁶Sr values in the least altered (Mn/Sr < 2.0) marine dolostones are 0.70418–0.70442 and 0.70343–0.70409 for the earlier and late phases of the marine transgression, respectively. The decrease in the ⁸⁷Sr/⁸⁶Sr ratio in ca. 2.1 Ga seawater is attributable to an increase in hydrothermal flux Sr into the Palaeoproterozoic ocean.     

Behaviour of a small sedimentary volcanic aquifer receiving irrigation return flows: La Aldea,Gran Canaria,Canary Islands (Spain)

In many arid and semi-arid areas, intensive cultivation is practiced despite water commonly being a limiting factor. Often, irrigation water is from local aquifers or imported from out-of-area aquifers and surface reservoirs. Irrigation return flows become a significant local recharge source, but they may deteriorate aquifer water quality. La Aldea valley, located in the western sector of Gran Canaria Island (Atlantic Ocean), is a coastal, half-closed depression in altered, low-permeability volcanics with alluvium in the gullies and scree deposits over a large part of the area. This area is intensively cultivated. Irrigation water comes from reservoirs upstream and is supplemented (average 30 %) by local groundwater; supplementation goes up to 70 % in dry years, in which groundwater reserves are used up to exhaustion if the dry period persists. Thus, La Aldea aquifer is key to the water-supply system, whose recharge is mostly from return irrigation flows and the scarce local rainfall recharge on the scree formations, conveyed to the gully deposits. To quantify the hydrogeological conceptual model and check data coherence, a simplified numerical model has been constructed, which can be used as a tool to help in water management.     

Submersible observations of an iceberg pit and scour on the Grand Banks of Newfoundland

An elongated sea bed depression (pit), 80 by 125 m and 10 m deep located on the Grand Banks of Newfoundland was identified on high resolution surficial geophysical survey lines and investigated using a manned submersible. The pit occurs at the termination of a long (greater than 3 km) iceberg scour, 1.5 m deep, and 80 m wide in 91 m of water. It is similar to other depressions in the region that indent the sea bed up to seven times deeper than the surrounding deepest iceberg furrows. The pit is interpreted to be the sea bed response to impact and loading by a rolling and grounding iceberg. The iceberg scour and pit appear to be of recent origin.