The Development of an Optimal Grid Coarsening Scheme: Interplay of Fluid Forces and Higher Moments of Fine-scale Flow Data

2005 ◽  
Author(s):  
N.H. Darman ◽  
K.S. Sorbie ◽  
G.E. Pickup
SPE Journal ◽  
2006 ◽  
Vol 11 (03) ◽  
pp. 317-327 ◽  
Author(s):  
Martin Mlacnik ◽  
Louis J. Durlofsky ◽  
Zoltan E. Heinemann

Summary A technique for the sequential generation of perpendicular-bisectional (PEBI) grids adapted to flow information is presented and applied. The procedure includes a fine-scale flow solution, the generation of an initial streamline-isopotential grid, grid optimization, and upscaling. The grid optimization is accomplished through application of a hybrid procedure with gradient and Laplacian smoothing steps, while the upscaling is based on a global-local procedure that makes use of the global solution used in the grid-determination step. The overall procedure is successfully applied to a complex channelized reservoir model involving changing well conditions. The gridding and upscaling procedures presented here may also be suitable for use with other types of structured or unstructured grid systems. Introduction Modern geological and geostatistical tools provide highly detailed descriptions of the spatial variation of reservoir properties, resulting in fine-grid models consisting of 107 to 108 gridblocks. As a consequence of this high level of detail, these models cannot be used directly in numerical reservoir simulators, but need to be coarsened significantly. Coarsening requires the averaging of rock parameters from the fine scale to the coarse scale. This process is referred to as upscaling. For simulation of flow in porous media, the upscaling of permeability is of particular interest. A large body of literature exists on this topic; for a comprehensive review of existing techniques, see Durlofsky (2005). To preserve as much of the geological information of the fine grid as possible, the grid coarsening should not be performed uniformly, but with more refinement in areas that are expected to have large impact on the flow, including structural features, such as faults. Although grid-generation techniques based on purely static, nonflow-based considerations have been shown to produce reasonable results(Garcia et al. 1992), the application of flow-based grids is often preferable. Flow-based grids require the solution of some type of fine-scale problem. They are then constructed by exploiting the information obtained from streamlines (and possibly isopotentials) either directly or indirectly. Depending on the type of grid used, points will be defined as cell vertices or nodes, resulting in either a corner-point geometry or point-distributed grid. Several gridding techniques for reservoir simulation have been introduced along these lines, as we now discuss.


2021 ◽  
Vol 17 (2) ◽  
pp. e1008588
Author(s):  
Constanze Ciavarella ◽  
Neil M. Ferguson

The spatial dynamics of epidemics are fundamentally affected by patterns of human mobility. Mobile phone call detail records (CDRs) are a rich source of mobility data, and allow semi-mechanistic models of movement to be parameterised even for resource-poor settings. While the gravity model typically reproduces human movement reasonably well at the administrative level spatial scale, past studies suggest that parameter estimates vary with the level of spatial discretisation at which models are fitted. Given that privacy concerns usually preclude public release of very fine-scale movement data, such variation would be problematic for individual-based simulations of epidemic spread parametrised at a fine spatial scale. We therefore present new methods to fit fine-scale mathematical mobility models (here we implement variants of the gravity and radiation models) to spatially aggregated movement data and investigate how model parameter estimates vary with spatial resolution. We use gridded population data at 1km resolution to derive population counts at different spatial scales (down to ∼ 5km grids) and implement mobility models at each scale. Parameters are estimated from administrative-level flow data between overnight locations in Kenya and Namibia derived from CDRs: where the model spatial resolution exceeds that of the mobility data, we compare the flow data between a particular origin and destination with the sum of all model flows between cells that lie within those particular origin and destination administrative units. Clear evidence of over-dispersion supports the use of negative binomial instead of Poisson likelihood for count data with high values. Radiation models use fewer parameters than the gravity model and better predict trips between overnight locations for both considered countries. Results show that estimates for some parameters change between countries and with spatial resolution and highlight how imperfect flow data and spatial population distribution can influence model fit.


SPE Journal ◽  
2001 ◽  
Vol 6 (01) ◽  
pp. 47-56 ◽  
Author(s):  
N.H. Darman ◽  
L.J. Durlofsky ◽  
K.S. Sorbie ◽  
G.E. Pickup
Keyword(s):  

SPE Journal ◽  
2018 ◽  
Vol 23 (02) ◽  
pp. 614-624 ◽  
Author(s):  
Shouhong Du ◽  
Larry S. Fung ◽  
Ali H. Dogru

Summary Grid coarsening outside of the areas of interest is a common method to reduce computational cost in reservoir simulation. Aquifer regions are candidates for grid coarsening. In this situation, upscaling is applied to the fine grid to generate coarse-grid flow properties. The efficacy of the approach can be judged easily by comparing the simulation results between the coarse-grid model and the fine-grid model. For many reservoirs in the Middle East bordered by active aquifers, transient water influx is an important recovery mechanism that needs to be modeled correctly. Our experience has shown that the standard grid coarsening and upscaling method do not produce correct results in this situation. Therefore, the objective of this work is to build a method that retains the fine-scale heterogeneities to accurately represent the water movement, but to significantly reduce the computational cost of the aquifer grids in the model. The new method can be viewed as a modified two-level multigrid (MTL-MG) or a specialized adaptation of the multiscale method. It makes use of the vertical-equilibrium (VE) concept in the fine-scale pressure reconstruction in which it is applicable. The method differs from the standard grid coarsening and upscaling method in which the coarse-grid properties are computed a priori. Instead, the fine-scale information is restricted to the coarse grid during Newton's iteration to represent the fine-scale flow behavior. Within the aquifer regions, each column of fine cells is coarsened vertically based on fine-scale z-transmissibility. A coarsened column may consist of a single amalgamated aquifer cell or multiple vertically disconnected aquifer cells separated by flow barriers. The pore volume (PV), compressibility, and lateral flow terms of the coarse cell are restricted from the fine-grid cells. The lateral connectivity within the aquifer regions and the one between the aquifer and the reservoir are honored, inclusive of the fine-scale description of faults, pinchouts, and null cells. Reservoir regions are not coarsened. Two alternatives exist for the fine-scale pressure reconstruction from the coarse-grid solution. The first method uses the VE concept. When VE applies, pressure variation can be analytically computed in the solution update step. Otherwise, the second method is to apply a 1D z-line solve for the fine-scale aquifer pressure from the coarse-grid solution. Simulation results for several examples are included to demonstrate the efficacy and efficiency of the method. We have applied the method to several Saudi Arabian complex full-field simulation models in which the transient aquifer water influx has been identified as a key factor. These models include dual-porosity/dual-permeability (DPDP) models, as well as models with faults and pinchouts in corner-point-geometry grids, for both history match and prediction period. The method is flexible and allows for the optional selection of aquifer regions to be coarsened, either only peripheral aquifers or both the peripheral and bottom aquifers. The new method gives nearly identical results compared with the original runs without coarsening, but with significant reduction in computer time or hardware cost. These results will be detailed in the paper.


Author(s):  
Russell L. Steere

Complementary replicas have revealed the fact that the two common faces observed in electron micrographs of freeze-fracture and freeze-etch specimens are complementary to each other and are thus the new faces of a split membrane rather than the original inner and outer surfaces (1, 2 and personal observations). The big question raised by published electron micrographs is why do we not see depressions in the complementary face opposite membrane-associated particles? Reports have appeared indicating that some depressions do appear but complementarity on such a fine scale has yet to be shown.Dog cardiac muscle was perfused with glutaraldehyde, washed in distilled water, then transferred to 30% glycerol (material furnished by Dr. Joaquim Sommer, Duke Univ., and VA Hospital, Durham, N.C.). Small strips were freeze-fractured in a Denton Vacuum DFE-2 Freeze-Etch Unit with complementary replica tooling. Replicas were cleaned in chromic acid cleaning solution, then washed in 4 changes of distilled water and mounted on opposite sides of the center wire of a Formvar-coated grid.


2016 ◽  
Vol 550 ◽  
pp. 147-161 ◽  
Author(s):  
PM Zhadan ◽  
MA Vaschenko ◽  
VB Lobanov ◽  
AF Sergeev ◽  
SA Kotova

2020 ◽  
Vol 655 ◽  
pp. 185-198
Author(s):  
J Weil ◽  
WDP Duguid ◽  
F Juanes

Variation in the energy content of prey can drive the diet choice, growth and ultimate survival of consumers. In Pacific salmon species, obtaining sufficient energy for rapid growth during early marine residence is hypothesized to reduce the risk of size-selective mortality. In order to determine the energetic benefit of feeding choices for individuals, accurate estimates of energy density (ED) across prey groups are required. Frequently, a single species is assumed to be representative of a larger taxonomic group or related species. Further, single-point estimates are often assumed to be representative of a group across seasons, despite temporal variability. To test the validity of these practices, we sampled zooplankton prey of juvenile Chinook salmon to investigate fine-scale taxonomic and temporal differences in ED. Using a recently developed model to estimate the ED of organisms using percent ash-free dry weight, we compared energy content of several groups that are typically grouped together in growth studies. Decapod megalopae were more energy rich than zoeae and showed family-level variability in ED. Amphipods showed significant species-level variability in ED. Temporal differences were observed, but patterns were not consistent among groups. Bioenergetic model simulations showed that growth rate of juvenile Chinook salmon was almost identical when prey ED values were calculated on a fine scale or on a taxon-averaged coarse scale. However, single-species representative calculations of prey ED yielded highly variable output in growth depending on the representative species used. These results suggest that the latter approach may yield significantly biased results.


Sign in / Sign up

Export Citation Format

Share Document