Computing derivative information of sequentially coupled subsurface models

A generic framework for the computation of derivative information required for gradient-based optimization using sequentially coupled subsurface simulation models is presented. The proposed approach allows for the computation of any derivative information with no modification of the mathematical framework. It only requires the forward model Jacobians and the objective function to be appropriately defined. The flexibility of the framework is demonstrated by its application in different reservoir management studies. The performance of the gradient computation strategy is demonstrated in a synthetic water-flooding model, where the forward model is constructed based on a sequentially coupled flow-transport system. The methodology is illustrated for a synthetic model, with different types of applications of data assimilation and life-cycle optimization. Results are compared with the classical fully coupled (FIM) forward simulation. Based on the presented numerical examples, it is demonstrated how, without any modifications of the basic framework, the solution of gradient-based optimization models can be obtained for any given set of coupled equations. The sequential derivative computation methods deliver similar results compared to FIM methods, while being computationally more efficient.     

Participatory Vulnerability Reduction (PVR): an urban community-based approach for earthquake management

Participation of local people during any disaster is enormous. They possess better knowledge and information about their own community than anyone else from the outside. This study proposes Participatory Vulnerability Reduction (PVR), a community-based approach for disaster management. The concept of PVR was applied to an urban community of Dhaka city (Ward no. 06 of Dhaka North City Corporation) which has been identified by the Comprehensive Disaster Management Programme as one of the most vulnerable areas of the city for earthquake. PVR consists three steps, and in each step, different participatory urban appraisal tools were used. In the first step, the community people assessed the earthquake vulnerability. It was found that some certain parts of the study area are highly vulnerable due to lack of accessibility to the critical facilities, inadequacy of open space, poor construction practice and unsuitable soil condition for building construction. This was followed by analyzing the root causes and effects of these problems. Structural fragility of the buildings, construction of settlements by filling the low-lying areas and development of slums beside taller buildings are the three major causes behind the above vulnerable issues. In the second step, capacity of the community was assessed in terms of resources and their organizational structure. In the final step, local people developed the strategies to overcome the vulnerability and a community-based organizational set up was proposed to coordinate the collective actions. Although developed in local context, application of PVR is not limited for earthquake and it can be replicated for other communities as well.     

An orthotropic interface damage model for simulating drying processes in soils

The study of drying process in soils has received an increased attention in the last few years. This is very complex phenomenon that generally leads to the formation and propagation of desiccation cracks in the soil mass. In recent engineering applications, high aspect ratio elements have proved to be well suited to tackle this type of problem using finite elements. However, the modeling of interfaces between materials with orthotropic properties that generally exist in this type of problem using standard (isotropic) constitutive model is very complex and challenging in terms of the mesh generation, leading to very fine meshes that are intensive CPU demanding. A novel orthotropic interface mechanical model based on damage mechanics and capable of dealing with interfaces between materials in which the strength depends on the direction of analysis is proposed in this paper. The complete mathematical formulation is presented together with the algorithm suggested for its numerical implementation. Some simple yet challenging synthetic benchmarks are analyzed to explore the model capabilities. Laboratory tests using different textures at the contact surface between materials were conducted to evaluate the strengths of the interface in different directions. These experiments were then used to validate the proposed model. Finally, the approach is applied to simulate an actual desiccation test involving an orthotropic contact surface. In all the application cases the performance of the model was very satisfactory.     

Early Permian transgressive–regressive cycles: Sequence stratigraphic reappraisal of the coal-bearing Barakar Formation,Raniganj Basin,India

The present research is an attempt to assess the Barakar Formation of the Raniganj Gondwana Basin, India, in the frame of fluvio-marine (estuarine) depositional systems using sequence stratigraphic elements. Analysis of predominant facies associations signify deposition in three sub-environments: (i) a river-dominated bay-head delta zone in the inner estuary, with transition from braided fluvial channels (FA-B1) to tide-affected meandering fluvial channels and flood plains (FA-B2) in the basal part of the succession; (ii) a mixed energy central basin zone, which consists of transitional fluvio-tidal channels (FA-B2), tidal flats, associated with tidal channels and bars (FA-B3) in the middle-upper part of the succession; and (iii) a wave-dominated outer estuary (coastal) zone (FA-B4 with FA-B3) in the upper part of the succession. Stacked progradational (P1, P2)–retrogradational (R1, R2) successions attest to one major base level fluctuation, leading to distinct transgressive–regressive (T–R) cycles with development of initial falling stage systems tract (FSST), followed by lowstand systems tract (LST) and successive transgressive systems tracts (TST-1 and TST-2). Shift in the depositional regime from regressive to transgressive estuarine system in the early Permian Barakar Formation is attributed to change in accommodation space caused by mutual interactions of (i) base level fluctuations in response to climatic amelioration and (ii) basinal tectonisms (exhumation/sagging) related to post-glacial isostatic adjustments in the riftogenic Gondwana basins.     

Development of extended WRF variational data assimilation system (WRFDA) for WRF non-hydrostatic mesoscale model

The paper intends to present the development of the extended weather research forecasting data assimilation (WRFDA) system in the framework of the non-hydrostatic mesoscale model core of weather research forecasting system (WRF-NMM), as an imperative aspect of numerical modeling studies. Though originally the WRFDA provides improved initial conditions for advanced research WRF, we have successfully developed a unified WRFDA utility that can be used by the WRF-NMM core, as well. After critical evaluation, it has been strategized to develop a code to merge WRFDA framework and WRF-NMM output. In this paper, we have provided a few selected implementations and initial results through single observation test, and background error statistics like eigenvalues, eigenvector and length scale among others, which showcase the successful development of extended WRFDA code for WRF-NMM model. Furthermore, the extended WRFDA system is applied for the forecast of three severe cyclonic storms: Nargis (27 April–3 May 2008), Aila (23–26 May 2009) and Jal (4–8 November 2010) formed over the Bay of Bengal. Model results are compared and contrasted within the analysis fields and later on with high-resolution model forecasts. The mean initial position error is reduced by 33% with WRFDA as compared to GFS analysis. The vector displacement errors in track forecast are reduced by 33, 31, 30 and 20% to 24, 48, 72 and 96 hr forecasts respectively, in data assimilation experiments as compared to control run. The model diagnostics indicates successful implementation of WRFDA within the WRF-NMM system.     

An ichnofossil assemblage from the fluvial deposits of the Upper Pliocene–Pleistocene Pinjor Formation (Siwalik Group), northwestern Himalayas,India: Palaeoenvironmental implications

Continental (fluvial) strata of the Pinjor Formation (Siwalik Group), northwestern Himalayas, India, contain an invertebrate trace fossil assemblage containing Planolites beverleyensis, Palaeophycus isp., Scoyenia gracilis, Taenidium barretti and other undifferentiated traces. The traces are found in an \(\sim \)26 m thick interval of alternating pinkish red siltstone, which is intercalated with mudstone, and thickly-bedded buff and greenish coloured sandstone. These sediments are interpreted as the deposits of floodplains and channel-bars of fluvial environments and low-energy overbank floodplain deposits. The trace fossils studied here are the first well documented ichnofossil assemblage from the vast, late Cenozoic Siwalik depositional system. They are not only of palaeoenvironmental significance, but they add to the growing ichnofossil database in facies of fluvial origin and should be an impetus to further ichnological studies of the Siwalik Group.     

Modulus of subgrade reaction that varies with magnitude of displacement of cohesionless soil

Modulus of subgrade reaction is one of the required design parameters in any structural analyses of shallow foundations. However, the constant values of modulus of subgrade reaction that are determined from either from literature studies or the results of plate load tests, regardless of magnitude of soil’s displacement under design loads, have been used in structural designs of foundations. In this study, the results of 43 published full-scale field plate load tests in cohesionles soils were gathered to expose any variation in the values of modulus of subgrade reaction as soil’s displacement increases. Extensive finite element (FE) analyses were carried out while the results of FE analyses were compared with actual measured results of field load tests. The results of analyses indicated that both the modulus of subgrade reaction decrease with increase in magnitude of displacement of soils and internal forces of a design of structural frame is higher with the values of modulus of subgrade reaction that are sensitive to soil’s displacement are implemented into analyses. Therefore, structural dimensions of any structure with the constant values of modulus of subgrade reaction would not be a precise engineering solution.     

Maturity modeling and burial history reconstruction for Garau and Sargelu formations in Lurestan basin,south Iran

The Upper Jurassic Sargelu and the Cretaceous Garau formations are important source rocks in the Lurestan basin, southwest Iran. In this study, the maturity evolution and burial history of these source rocks are investigated using 1D modeling in 15 wells within the Lurestan basin. These models are calibrated using temperature and vitrinite reflectance readings at well locations. Results indicate that thermal maturity of studied source rocks increases towards the northwest of the Lurestan basin and is the lowest over the Anaran High. The present-day level of maturity suggests that these formations are suitable candidates for shale gas resources in the Lurestan basin.     

2D seismic reflection contribution to structural and geometric study of Cenozoic aquifer systems in the Sisseb El-Alem basin,central-eastern Tunisia

The Sisseb-ElAlem basin in central-eastern Tunisia contains several aquifers intensely used in agriculture. Highly exploited shallow aquifers are vulnerable to pollution and drought and reveal high salinity. Consequently, the new hydrogeological surveys were oriented to deep aquifers. The present work corresponds to study of 2D seismic sections, calibrated by water and petroleum wells, and geological outcrops. It highlights geometry and faulting of buried and preserved structures in depth and leads to identify the spatial distribution of Cenozoic potential reservoir layers in the different paleogeographic domains. The Sisseb-ElAlem basin is affected mainly by near E- and N-trending faults. Iso-depth and iso-thickness mapping were used to follow lateral and vertical variations of sedimentary series related to tectonics. The study area consists of Sisseb, Sbikha and El Ktifa uplifts, and El Alem, and Bled Saadia subsiding synclines. The Oligocene-Miocene reservoir series are exploited locally at the high borders. They reveal thickening and deepening southward in the El Alem-Bled Saadia domain, where they can be exploited at about 400 to 2000 m depths.This study should contribute to the orientation of water research to deeper aquifers, and therefore, ensure sufficient drinking and irrigation water supply mainly in areas called for sustainable development.