Integration of local–global upscaling and grid adaptivity for simulation of subsurface flow in heterogeneous formations

We propose a methodology, called multilevel local–global (MLLG) upscaling, for generating accurate upscaled models of permeabilities or transmissibilities for flow simulation on adapted grids in heterogeneous subsurface formations. The method generates an initial adapted grid based on the given fine-scale reservoir heterogeneity and potential flow paths. It then applies local–global (LG) upscaling for permeability or transmissibility [7], along with adaptivity, in an iterative manner. In each iteration of MLLG, the grid can be adapted where needed to reduce flow solver and upscaling errors. The adaptivity is controlled with a flow-based indicator. The iterative process is continued until consistency between the global solve on the adapted grid and the local solves is obtained. While each application of LG upscaling is also an iterative process, this inner iteration generally takes only one or two iterations to converge. Furthermore, the number of outer iterations is bounded above, and hence, the computational costs of this approach are low. We design a new flow-based weighting of transmissibility values in LG upscaling that significantly improves the accuracy of LG and MLLG over traditional local transmissibility calculations. For highly heterogeneous (e.g., channelized) systems, the integration of grid adaptivity and LG upscaling is shown to consistently provide more accurate coarse-scale models for global flow, relative to reference fine-scale results, than do existing upscaling techniques applied to uniform grids of similar densities. Another attractive property of the integration of upscaling and adaptivity is that process dependency is strongly reduced, that is, the approach computes accurate global flow results also for flows driven by boundary conditions different from the generic boundary conditions used to compute the upscaled parameters. The method is demonstrated on Cartesian cell-based anisotropic refinement (CCAR) grids, but it can be applied to other adaptation strategies for structured grids and extended to unstructured grids.     

In situ estimation of the effective chemical diffusion coefficient of a rock matrix in a fractured aquifer

An in situ method of estimating the effective diffusion coefficient for a chemical constituent that diffuses into the primary porosity of a rock is developed by abruptly changing the concentration of the dissolved constituent in a borehole in contact with the rock matrix and monitoring the time-varying concentration. The experiment was conducted in a borehole completed in mudstone on the campus of the University of the Free State in Bloemfontein, South Africa. Numerous tracer tests were conducted at this site, which left a residual concentration of sodium chloride in boreholes that diffused into the rock matrix over a period of years. Fresh water was introduced into a borehole in contact with the mudstone, and the time-varying increase of chloride was observed by monitoring the electrical conductivity (EC) at various depths in the borehole. Estimates of the effective diffusion coefficient were obtained by interpreting measurements of EC over 34 d. The effective diffusion coefficient at a depth of 36 m was approximately 7.8×10^?6 m²/d, but was sensitive to the assumed matrix porosity. The formation factor and mass flux for the mudstone were also estimated from the experiment.     

Hydrogeochemistry and isotope analyses used to determine groundwater recharge and flow in low-gradient catchments of the province of Buenos Aires,Argentina

The region known as Pampa Plain, in Argentina, is a vast area characterized by slopes of less than 0.05%. The surface sediment is silty sand, mainly Aeolian, referred to as Pampean Loess, which is a phreatic aquifer unit of utmost importance for the water supply of the region. On account of the slight gradient, hydrogeological analyses using only hydraulic measurements are difficult to perform, often leading to confusing results. Thus, the study presented relies on hydrochemical modeling and isotopic determinations as well. The study area comprises three catchments in the inter-mountainous area corresponding to El Moro, Tamangueyú and Seco creeks, and covering 2,570 km² in the province of Buenos Aires. Measurements of piezometric levels and samples of groundwater, surface water and rainwater were carried out between January and March 2005. Major ion results were analyzed by means of hydrochemical graphs and hydrogeochemical modelling using NETPATH. The resulting data show that it is possible to identify local changes in recharge, flow direction and stream/groundwater relationship by using hydrochemical and isotopic information, which may become a useful and more precise tool for the study of particularly flat landscapes.     

Estimates of heat flow and heat production and a thermal model of the São Francisco craton

An updated analysis of geothermal data from the highland area of eastern Brazil has been carried out and the characteristics of regional variations in geothermal gradients and heat flow examined. The database employed includes results of geothermal measurements at 45 localities. The results indicate that the Salvador craton and the adjacent metamorphic fold belts northeastern parts of the study area are characterized by geothermal gradients in the range of 6–17°C/km. The estimated heat flow values fall in the range of 28–53 mW/m², with low values in the cratonic area relative to the fold belts. On the other hand, the São Francisco craton and the intracratonic São Francisco sedimentary basin in the southwestern parts are characterized by relatively higher gradient values, in the range of 14–42°C/km, with the corresponding heat flow values falling in the range of 36–89 mW/m². Maps of regional variations indicate that high heat flow anomaly in the São Francisco craton is limited to areas of sedimentary cover, to the west of the Espinhaço mountain belt. Crustal thermal models have been developed to examine the implications of the observed intracratonic variations in heat flow. The thermal models take into consideration variation of thermal conductivity with temperature as well as change of radiogenic heat generation with depth. Vertical distributions of seismic velocities were used in obtaining estimates of radiogenic heat production in crustal layers. Crustal temperatures are calculated based on a procedure that makes simultaneous use of the Kirchoff and Generalized Integral Transforms, providing thereby analytical solutions in 2D and 3D geometry. The results point to temperature variations of up to 300°C at the Moho depth, between the northern Salvador and southern São Francisco cratons. There are indications that differences in rheological properties, related to thermal field, are responsible for the contrasting styles of deformation patterns in the adjacent metamorphic fold belts.     

The role of overwash in the evolution of mixing zone morphology within barrier islands

Overwash is a major controlling factor in the morphology of the mixing zone of coastal aquifers. Conceptual models of the mixing zone describe an interface controlled by tidal oscillations, wave run-up, and other factors; however, few describe the influence of large storm events. In August 1993, Hatteras Island, North Carolina, USA, experienced a 3-m storm surge due to Hurricane Emily. Sound-side flooding infiltrated a wellfield, causing a dramatic increase in TDS levels that persisted for more than 3 years. Two-dimensional simulations with SUTRA, the USGS finite-element model, are calibrated to the TDS breakthrough data of this storm to infer model dispersivity values. Simulations using the calibrated dispersivity values, predicted flooding levels, and 54 years of hurricane records to determine the influence of the overwash events suggest that it is rare for the mixing zone to approximate the conceptual morphology. Even during quiescent periods such as between 1965 and 1975, TDS levels do not return to theoretical levels before being elevated by a subsequent storm event. Thus, while tidal oscillations and other factors are important to mixing zone development, basic wind events and more severe storm events may have more influence and lasting effect on the morphology of the mixing zone.     

Surface water–groundwater interaction and chemistry in a mineral-armored hydrothermal outflow channel, Yellowstone National Park, USA

Small quantities of groundwater interact with hydrothermal surface water to drive in-stream geochemical processes in a silica-armored hot-spring outflow channel in Yellowstone National Park, USA. The objective of this study was to characterize the hydrology and geochemistry of this unique system in order to (1) learn more about the Yellowstone Plateau’s subsurface water mixing between meteoric and hydrothermal waters and (2) learn more about the chemical and physical processes that lead to accumulation of streambed cements, i.e., streambed armor. A combination of hydrological, geochemical, mineralogical, microscopic, and petrographic techniques were used to identify groundwater and surface-water exchange. Interaction could be identified in winter because of differences in surface water and groundwater composition but interaction at other times of the year cannot be ruled out. Dissolved constituents originating from groundwater (e.g., Fe(II) and Mg) were traced downstream until oxidation and/or subsequent precipitation with silica removed them, particularly where high affinity substrates like cyanobacterial surfaces were present. Because the stream lies in a relatively flat drainage basin and is fed mainly by a seasonally relatively stable hot spring, this system allowed study of the chemical processes along a stream without the obscuring effects of sedimentation.     

New basin modelling results from the Polish part of the Central European Basin system: implications for the Late Cretaceous–Early Paleogene structural inversion

Combined subsidence and thermal 1D modelling was performed on six well-sections located in the north-western Mid-Polish Trough/Swell in the eastern part of the Central European Basin system. The modelling allowed constraining quantitatively both the Mesozoic subsidence and the magnitude of the Late Cretaceous–Paleocene inversion and erosion. The latter most probably reached 2,400 m in the Mid-Polish Swell area. The modelled Upper Cretaceous thickness did not exceed 500 m, and probably corresponded to 200–300 m in the swell area as compared with more than 2,000 m in the adjacent non-inverted part of the basin. Such Upper Cretaceous thickness pattern implies early onset of inversion processes, probably in the Late Turonian or Coniacian. Our modelling, coupled with previous results of stratigraphic and seismic studies, demonstrates that the relatively low sedimentation rates in the inverted part of the basin during the Late Cretaceous were the net result of several discrete pulses of non-deposition and/or erosion that were progressively more pronounced towards the trough axis. The last phase of inversion started in the Late Maastrichtian and was responsible for the total amount of erosion, which removed also the reduced Upper Cretaceous deposits. According to our modelling results, a Late Cretaceous heat-flow regime which is similar to the present-day conditions (about 50 mW/m²) was responsible for the observed organic maturity of the Permian-Mesozoic rocks. This conclusion does not affect the possibility of Late Carboniferous–Permian and Late Permian–Early Triassic thermal events.     

Increased Terrestrial to Ocean Sediment and Carbon Fluxes in the Northern Chesapeake Bay Associated With Twentieth Century Land Alteration

We calculated Chesapeake Bay (CB) sediment and carbon fluxes before and after major anthropogenic land clearance using robust monitoring, modeling and sedimentary data. Four distinct fluxes in the estuarine system were considered including (1) the flux of eroded material from the watershed to streams, (2) the flux of suspended sediment at river fall lines, (3) the burial flux in tributary sediments, and (4) the burial flux in main CB sediments. The sedimentary maximum in Ambrosia (ragweed) pollen marked peak land clearance (~1900 a.d.). Rivers feeding CB had a total organic carbon (TOC)/total suspended solids of 0.24?±?0.12, and we used this observation to calculate TOC fluxes from sediment fluxes. Sediment and carbon fluxes increased by 138–269% across all four regions after land clearance. Our results demonstrate that sediment delivery to CB is subject to significant lags and that excess post-land clearance sediment loads have not reached the ocean. Post-land clearance increases in erosional flux from watersheds, and burial in estuaries are important processes that must be considered to calculate accurate global sediment and carbon budgets.     

印度板块和亚洲大陆在何时何地碰撞

印度板块和亚洲大陆的初始碰撞时间是所有相关的喜马拉雅-西藏造山体系演化模式的主控条件,并严重影响到对众多与青藏高原隆升和东亚大陆挤出相关的地质过程速率的解释,以及对新生代全球气候变化的理解。尽管印度板块和亚洲大陆汇聚的速率在55Ma突然减缓被广泛地认为是初始碰撞的标志,但这次碰撞所造成的主要构造效应直到20多个百万年以后才显现出来。对印度板块和亚洲大陆相对位置的重新估算,表明它们在55Ma时并没有达到可以彼此发生碰撞的距离。基于来自西藏新的野外证据和对已有数据的重新评估,认为初始碰撞发生在始新世—渐新世之交（约34Ma）,并对55Ma时发生的地质事件提出了另一种解释