A coarse-scale compositional model

In subsurface flow modeling, compositional simulation is often required to model complex recovery processes, such as gas/CO ₂ injection. However, compositional simulation on fine-scale geological models is still computationally expensive and even prohibitive. Most existing upscaling techniques focus on black-oil models. In this paper, we present a general framework to upscale two-phase multicomponent flow in compositional simulation. Unlike previous studies, our approach explicitly considers the upscaling of flow and thermodynamics. In the flow part, we introduce a new set of upscaled flow functions that account for the effects of compressibility. This is often ignored in the upscaling of black-oil models. In the upscaling of thermodynamics, we show that the oil and gas phases within a coarse block are not at chemical equilibrium. This non-equilibrium behavior is modeled by upscaled thermodynamic functions, which measure the difference between component fugacities among the oil and gas phases. We apply the approach to various gas injection problems with different compositional features, permeability heterogeneity, and coarsening ratios. It is shown that the proposed method accurately reproduces the averaged fine-scale solutions, such as component overall compositions, gas saturation, and density solutions in the compositional flow.     

Intensity measures for the seismic response evaluation of buried steel pipelines under near-field pulse-like ground motions

Ground-motion Intensity Measures (IMs) are used to quantify the strength of ground motions and evaluate the response of structures. IMs act as a link between seismic demand and seismic hazard analysis and therefore, have a key role in performance-based earthquake engineering. Many studies have been carried out on the determination of suitable IMs in terms of efficiency, sufficiency and scaling robustness. The majority of these investigations focused on ordinary structures such as buildings and bridges, and only a few were about buried pipelines. In the current study, the optimal IMs for predicting the seismic demand of continuous buried steel pipelines under near-field pulse-like ground motion records is investigated. Incremental dynamic analysis is performed using twenty ground motion records. Using the results of the regression analysis, the optimality of 23 potential IMs are studied. It is concluded that specific energy density (SED) followed by \(\sqrt {VSI[{\omega _1}(PGD + RM{S_d})]} \) are the optimal IMs based on efficiency, sufficiency and scaling robustness for seismic response evaluation of buried pipelines under near-field ground motions.

     

A multiple-point statistics algorithm for 3D pore space reconstruction from 2D images

Fluid flow behavior in a porous medium is a function of the geometry and topology of its pore space. The construction of a three dimensional pore space model of a porous medium is therefore an important first step in characterizing the medium and predicting its flow properties. A stochastic technique for reconstruction of the 3D pore structure of unstructured random porous media from a 2D thin section training image is presented. The proposed technique relies on successive 2D multiple point statistics simulations coupled to a multi-scale conditioning data extraction procedure. The Single Normal Equation Simulation Algorithm (SNESIM), originally developed as a tool for reproduction of long-range, curvilinear features of geological structures, serves as the simulation engine. Various validating criteria such as marginal distributions of pore and grain, directional variograms, multiple-point connectivity curves, single phase effective permeability and two phase relative permeability calculations are used to analyze the results. The method is tested on a sample of Berea sandstone for which a 3D micro-CT scanning image is available. The results confirm that the equi-probable 3D realizations obtained preserve the typical patterns of the pore space that exist in thin sections, reproduce the long-range connectivities, capture the characteristics of anisotropy in both horizontal and vertical directions and have single and two phase flow characteristics consistent with those of the measured 3D micro-CT image.     

Comparison of fundamental properties of new types of fiber‐mesh‐reinforced seismic isolators with conventional isolators

New types of fiber‐reinforced rubber‐based seismic isolators have been a research interest for a number of engineers in the past decade. These new types of isolators can have similar seismic performances compared with the conventional ones. In most of the previous researches, the fiber‐reinforced rubber‐based isolators is usually manufactured with placing fiber sheets between precut rubber layers with the use of a bonding agent. This research differs from the previous researches in terms of manufacturing process, use of fiber mesh instead of fiber sheets, and use of lead in the core for some of the bearings. The aim of this research is to provide comparisons in fundamental seismic response properties of the new type of fiber mesh reinforced isolators and conventional isolators. In this scope, four pairs of fiber mesh reinforced elastomeric bearings and four pairs of steel‐reinforced elastomeric bearings are subjected to various levels of compression stresses and cyclic shear strains under constant vertical pressure. The tested types of isolators are fiber mesh reinforced elastomeric bearing, fiber mesh reinforced elastomeric bearing with lead core, steel‐reinforced elastomeric bearings, and steel‐reinforced elastomeric bearings with lead core. In this research, steel‐reinforced bearings are called conventional isolators. The major advantage for fiber mesh reinforced bearings observed during the tests is that these isolators can develop a considerable low horizontal stiffness compared with the conventional isolators. The damping characteristics of the new and conventional types are similar to each other. Copyright © 2013 John Wiley & Sons, Ltd.     

System dynamics approach for simulating water resources of an urban water system with emphasis on sustainability of groundwater

Water resources management is an important driver in social and economic development. Water shortage is one of the most critical issues threatening human welfare, public health, and ecosystems. This issue has turned into a major challenge in many river basins all around the world due to the imbalance in water supply and demand. Use of simulation models can be effective tools in providing water managers with scientifically supported decisions in dealing with complex and uncertain water resource systems. System dynamics approach serves as a management tool and may play an important role in understanding the cause–effect in water resources systems. In the present study, system dynamics approach was applied to simulate management strategies dealing with Tehran metropolitan water resources systems. In the developed model, the trend of water storage in the next 30-year period and the effectiveness of water supply strategies were simulated. The results showed that, despite the growing shortage of the water resources, optimal use of existing resources under appropriate strategies could reduce water deficit within the next 30 years.     

Copper isotope fractionation in the Meiduk porphyry copper deposit,Northwest of Kerman Cenozoic magmatic arc,Iran

The Meiduk deposit possesses three different Cu reservoirs each with a unique Cu isotope signature. δ⁶⁵Cu for the leached cap minerals ranges from ?2.5 to +0.49‰ to ?0.45 to +0.3‰ for hypogene minerals and from +1.3 to +4.4‰ for supergene enrichment minerals. Oxidation of hypogene sulphides and effective trapping of copper (from solutions derived from the leached cap) in the supergene enrichment zone caused this relationship. A systematic pattern of low Cu isotope values close to the surface and higher isotope values with depth reveals a palaeo‐fluid pathway in the northwest–southeast direction over the deposit. Thus, the copper isotope data from leached cap and enrichment minerals can be used to monitor copper migration during supergene weathering at the Meiduk deposit.     

Enhancing PIV image and fractal descriptor for velocity and shear stresses propagation around a circular pier

In this study,the fractal dimensions of velocity fluctuations and the Reynolds shear stresses propagation for flow around a circular bridge pier are presented.In the study reported herein,the fractal dimension of velocity fluctuations(u′,v′,w′) and the Reynolds shear stresses(u′v′ and u′w′) of flow around a bridge pier were computed using a Fractal Interpolation Function(FIF) algorithm.The velocity fluctuations of flow along a horizontal plane above the bed were measured using Acoustic Doppler Velocity meter(ADV)and Particle Image Velocimetry(P1V).The PIV is a powerful technique which enables us to attain high resolution spatial and temporal information of turbulent flow using instantaneous time snapshots.In this study,PIV was used for detection of high resolution fractal scaling around a bridge pier.The results showed that the fractal dimension of flow fluctuated significantly in the longitudinal and transverse directions in the vicinity of the pier.It was also found that the fractal dimension of velocity fluctuations and shear stresses increased rapidly at vicinity of pier at downstream whereas it remained approximately unchanged far downstream of the pier.The higher value of fractal dimension was found at a distance equal to one times of the pier diameter in the back of the pier.Furthermore,the average fractal dimension for the streamwise and transverse velocity fluctuations decreased from the centreline to the side wall of the flume.Finally,the results from ADV measurement were consistent with the result from PIV,therefore,the ADV enables to detect turbulent characteristics of flow around a circular bridge pier.     

Evaluation of Land Surface Scheme SABAE-HW in Simulating Snow Depth,Soil Temperature and Soil Moisture within the BOREAS Site,Saskatchewan

The Soil Atmosphere Boundary, Accurate Evaluation of Heat and Water (SABAE-HW) model is a multilayered, one-dimensional, physically based version of the Canadian Land Surface Scheme (CLASS) and uses the same methodologies as CLASS, version 2.6. SABAE provides an improved interface for groundwater modelling to simulate soil moisture, soil temperature, energy fluxes and snow depth for a wide range of soil and vegetation. This paper reports the results of the first field comparison of SABAE-HW using an extensive ten-year dataset from the Boreal Ecosystem Atmosphere Study (BOREAS) and the Boreal Ecosystem Research and Monitoring Sites (BERMS) project, an area in central Saskatchewan, Canada, rich in terms of hydrological and meteorological data. The model is also independently tested and verified with the Simultaneous Heat and Water (SHAW) model, which is an unsaturated-zone transport model. Two boundary conditions are considered at the bottom of the soil profile: a water table boundary condition and a unit gradient boundary condition. There was substantial agreement between the results of the simulations and observations in terms of snow depth and soil temperature. Snow depth and soil temperature were simulated reasonably well by SABAE, with correlation values of 0.96 and 0.98, respectively. However, there were some discrepancies for simulated soil temperature in winter. General agreement was obtained in terms of unfrozen soil moisture results, especially at greater depths, but there were general similarities in observed and simulated soil moisture trends in winter. An average correlation of 0.55 was found for SABAE while the correlation for SHAW was much smaller (less than 0.30), which indicates a better fit between simulated and field data by SABAE. Although a unit gradient boundary condition does not influence soil moisture, it was found that unit gradient boundary runs resulted in increased bias towards overestimation of the soil temperature. Thus, a safer and more accurate approach, we believe, is to adopt a first type boundary (i.e., water table) condition at the bottom of the domain. This has implications for climate and weather modelling in general. The result of this field testing demonstrated the potential and high accuracy of SABAE-HW as a Canadian model capable of simulating snow depth, snow temperature, soil moisture, energy fluxes, and we believe it is now appropriate to include this land surface scheme with its counterparts.

R ésumé ?[Traduit par la rédaction] Le modèle Soil Atmosphere Boundary, Accurate Evaluation of Heat and Water (SABAE-HW) est une version multicouche, à une dimension, basée sur la physique du schéma CLASS (Canadian Land Surface Scheme) qui utilise les mêmes méthodologies que le CLASS version 2.6. Le SABAE offre une interface améliorée pour la modélisation des eaux sous-terraines permettant de simuler l'humidité du sol, la température du sol, les flux d'énergie et l'épaisseur de la neige pour une grande variété de sols et de végétation. Cet article présente les résultats de la première comparaison terrain du SABAE-HW en utilisant une base de données étendue de dix ans de l'Étude de l'atmosphère et des écosystèmes boréaux (BOREAS) et du projet des Sites de recherche et de surveillance des écosystèmes boréaux (BERMS), une région du centre de la Saskatchewan, au Canada, riche en données hydrologiques et météorologiques. Le modèle est aussi indépendamment testé et vérifié à l'aide du Simultaneous Heat and Water (SHAW), un modèle de transport en zone non saturée. Deux conditions aux limites sont supposées au fond du profil du sol : une condition aux limites de nappe phréatique et une condition aux limites de gradient unitaire. On a trouvé une concordance importante entre les résultats des simulations et les observations en ce qui a trait à l'épaisseur de la neige et à la température du sol. L'épaisseur de la neige et la température du sol ont été raisonnablement bien simulées par le modèle SABAE, avec des corrélations de 0,96 et 0,98, respectivement. Cependant, il y avait certaines divergences pour la température simulée du sol en hiver. Pour ce qui est des résultats concernant l'humidité du sol non gelé, ils s'accordaient généralement, surtout pour les plus grandes profondeurs, mais il y avait des similarités générales dans les tendances observées et simulées de l'humidité du sol en hiver. Nous avons trouvé une corrélation moyenne de 0,55 pour le SABAE alors que la corrélation pour le SHAW était beaucoup plus faible (0,30), ce qui indique un meilleur ajustement des données simulées aux données de terrain pour le SABAE. Même si une condition aux limites de gradient unitaire n'influence pas l'humidité du sol, il ressort que des passes faites avec un gradient unitaire aux limites ont produit un biais accru vers la surestimation de la température du sol. Donc, nous croyons qu'une approche plus sûre et plus précise serait d'adopter une condition aux limites de Dirichlet (c.-à-d. une nappe phréatique) au fond du domaine. Ceci a des répercussions sur la modélisation du climat et du temps en général. Le résultat de cet essai sur le terrain a démontré le potentiel et la grande exactitude du SABAE-HW en tant que modèle canadien capable de simuler l'épaisseur de la neige, la température du sol, l'humidité du sol et les flux d'énergie et nous croyons qu'il est maintenant approprié d'inclure ce schéma de surface avec ses contreparties.