Upscaling permeability for three-dimensional fractured porous rocks with the multiple boundary method

Upscaling permeability of grid blocks is crucial for groundwater models. A novel upscaling method for three-dimensional fractured porous rocks is presented. The objective of the study was to compare this method with the commonly used Oda upscaling method and the volume averaging method. First, the multiple boundary method and its computational framework were defined for three-dimensional stochastic fracture networks. Then, the different upscaling methods were compared for a set of rotated fractures, for tortuous fractures, and for two discrete fracture networks. The results computed by the multiple boundary method are comparable with those of the other two methods and fit best the analytical solution for a set of rotated fractures. The errors in flow rate of the equivalent fracture model decrease when using the multiple boundary method. Furthermore, the errors of the equivalent fracture models increase from well-connected fracture networks to poorly connected ones. Finally, the diagonal components of the equivalent permeability tensors tend to follow a normal or log-normal distribution for the well-connected fracture network model with infinite fracture size. By contrast, they exhibit a power-law distribution for the poorly connected fracture network with multiple scale fractures. The study demonstrates the accuracy and the flexibility of the multiple boundary upscaling concept. This makes it attractive for being incorporated into any existing flow-based upscaling procedures, which helps in reducing the uncertainty of groundwater models.     

The bubbles or the boiling pot? An ecosystemic approach to culture,environment and quality of life

For the diagnosis and prognosis of the problems of quality of life, a multidisciplinary ecosystemic approach encompasses four dimensions of being-in-the-world, as donors and recipients: intimate, interactive, social and biophysical. Social, cultural and environmental vulnerabilities are understood and dealt with, in different circumstances of space and time, as the conjugated effect of all dimensions of being-in-the-world, as they induce the events (deficits and assets), cope with consequences (desired or undesired) and contribute for change. Instead of fragmented and reduced representations of reality, diagnosis and prognosis of cultural, educational, environmental and health problems considers the connections (assets) and ruptures (deficits) between the different dimensions, providing a planning model to develop and evaluate research, teaching programmes, public policies and field projects. The methodology is participatory, experiential and reflexive; heuristic-hermeneutic processes unveil cultural and epistemic paradigms that orient subject-object relationships; giving people the opportunity to reflect on their own realities, engage in new experiences and find new ways to live better in a better world. The proposal is a creative model for thought and practice, providing many opportunities for discussion, debate and development of holistic projects integrating different scientific domains (social sciences, psychology, education, philosophy, etc.). Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.

HP–UHP metamorphism as an indicator of slab dip variations in the Alpine arc

HP/UHP and LT metamorphic units that commonly occur in the inner parts of mountain belts result from the subduction of continental and oceanic material, most often exhumed prior to continental collision. The prograde pressure–temperature history of HP–UHP rocks strongly depends on the convergence rate and on the subduction zone geometry. The maximum pressure recorded provides a proxy for the depth of shearing off and stacking of HP metamorphic nappes. A 2-D thermal model of continental subduction at lithospheric scale is used to compute the length and pressure peak of detached HP metamorphic units as a function of the slab dip angle and the convergence rate. Model results are applied to the metamorphic nappe pile of the inner Alps. A mean convergence rate of 1 cm/year during the subduction of the Briançonnais terrane is indicated by the paleogeographic reconstructions between 46 and 38 Ma. On this basis, the available petrological data and lengths of metamorphic units are used to compute the variations of the slab dip angle. The slab dip angle is shown to increase, from the northeast to the southwest, along the Alpine arc with estimated values of 20° for Suretta, 30–45° for Monte Rosa and Gran Paradiso, and 60° for Dora Maira. From Eocene to Oligocene times, the increase in slab dip angle is controlled by changes of buoyancy, due to the spatial configuration of the Valaisan trough and the incoming of crustal material within the subduction zone.     

Ensemble filter methods with perturbed observations applied to nonlinear problems

The ensemble Kalman filter (EnKF) has been shown repeatedly to be an effective method for data assimilation in large-scale problems, including those in petroleum engineering. Data assimilation for multiphase flow in porous media is particularly difficult, however, because the relationships between model variables (e.g., permeability and porosity) and observations (e.g., water cut and gas–oil ratio) are highly nonlinear. Because of the linear approximation in the update step and the use of a limited number of realizations in an ensemble, the EnKF has a tendency to systematically underestimate the variance of the model variables. Various approaches have been suggested to reduce the magnitude of this problem, including the application of ensemble filter methods that do not require perturbations to the observed data. On the other hand, iterative least-squares data assimilation methods with perturbations of the observations have been shown to be fairly robust to nonlinearity in the data relationship. In this paper, we present EnKF with perturbed observations as a square root filter in an enlarged state space. By imposing second-order-exact sampling of the observation errors and independence constraints to eliminate the cross-covariance with predicted observation perturbations, we show that it is possible in linear problems to obtain results from EnKF with observation perturbations that are equivalent to ensemble square-root filter results. Results from a standard EnKF, EnKF with second-order-exact sampling of measurement errors that satisfy independence constraints (EnKF (SIC)), and an ensemble square-root filter (ETKF) are compared on various test problems with varying degrees of nonlinearity and dimensions. The first test problem is a simple one-variable quadratic model in which the nonlinearity of the observation operator is varied over a wide range by adjusting the magnitude of the coefficient of the quadratic term. The second problem has increased observation and model dimensions to test the EnKF (SIC) algorithm. The third test problem is a two-dimensional, two-phase reservoir flow problem in which permeability and porosity of every grid cell (5,000 model parameters) are unknown. The EnKF (SIC) and the mean-preserving ETKF (SRF) give similar results when applied to linear problems, and both are better than the standard EnKF. Although the ensemble methods are expected to handle the forecast step well in nonlinear problems, the estimates of the mean and the variance from the analysis step for all variants of ensemble filters are also surprisingly good, with little difference between ensemble methods when applied to nonlinear problems.     

Precise ionosphere-free single-frequency GNSS positioning

Ionospheric delays can be efficiently eliminated from single-frequency data using a combination of carrier phases and code ranges. Unfortunately, GPS and GLONASS ranges are relatively noisy which can limit the use of the positioning method. Nevertheless, position standard deviations are in the range of 6–8 cm (horizontal) and 7–9 cm (3d) obtained from diurnal data batches from selected IGS reference stations can be further reduced to 2–3 cm (3d) for weekly smoothed averages. GPS data sets collected in Ghana (Africa) reveal a typical level of 10 cm of deviation that must be anticipated under average conditions. Looking at the future of GNSS, the European Galileo system will, in contrast to GPS, provide the broadband signal E5 that is by far less affected by multipath thus providing rather precise range measurements. Simulated processing runs featuring both high ionospheric and tropospheric delay variations show a 3d position precision of 4 cm even for a data batch as short as just 1 h, whereas GPS L1/Galileo E1 performance is close to 13 cm for the same data set.     

中国东北地区林地面积变化的动态模拟

There is plenty of forests in Northeast China which contributes a lot to the conservation of water and land resources, produces timber products, and provides habitats for a huge number of wild animals and plants. With changes of socio-economic factors as well as the geophysical conditions, there are dramatic changes on the spatial patterns of forest area. In this sense, it is of great significance to shed light on the dynamics of forest area changes to find the underlining reasons for shaping the changing patterns of forest area in Northeast China. To explore the dynamics of forest area change in Northeast China, an econometric model is developed which is composed of three equations identifying forestry production, conversion from open forest to closed forest and conversion from other land uses to closed forest so as to explore the impacts on the forest area changes from demographic, social, economic, location and geophysical factors. On this basis, we employ the Dynamics of Land System (DLS) model to simulate land-use conversions between forest area and non-forest cover and the land-use conversions within the sub-classes of forest area for the period 2000–2020 under business as usual scenario, environmental protection scenario and economic growth scenario. The simulation results show that forest area will expand continuously and there exist various kinds of changing patterns for the sub-classes of forest area, for example, closed forest will expand continuously and open forest and shrub will decrease a little bit, while area of other forest will keep intact. The research results provide meaningful decision-making information for conserving and exploiting the forest resources and making out the planning for forestry production in the Northeast China region.     

Seismic isolation analysis of FPS bearings in spatial lattice shell structures

<正>A theoretical model of a friction pendulum system(FPS) is introduced to examine its application for the seismic isolation of spatial lattice shell structures.An equation of motion of the lattice shell with FPS bearings is developed.Then, seismic isolation studies are performed for both double-layer and single-layer lattice shell structures under different seismic input and design parameters of the FPS.The influence of frictional coefficients and radius of the FPS on seismic performance are discussed.Based on the study,some suggestions for seismic isolation design of lattice shells with FPS bearings are given and conclusions are made which could be helpful in the application of FPS.     

Recent evolution of Lake Arreo,northern Spain: influences of land use change and climate

We present a high-resolution, multiproxy reconstruction of the depositional history of Lake Arreo, northern Spain, for the last 60 years. We conducted sedimentological, geochemical and diatom analyses in short cores and made a detailed comparison with regional instrumental climate data (1952–2007), limnological monitoring of the lake (1992–2008) and recent land use changes that affect the lake catchment. Chronology is based on “floating” discontinuous varve counts and ¹³⁷Cs and ¹⁴C dates. Four periods were identified in the Lake Arreo recent history: (1) prior to 1963, varved facies intercalated with fine turbidite deposits, and diatom assemblages dominated by Cyclotella taxa indicate predominantly meromictic conditions, (2) from 1964 to 1978, permanent anoxia persisted in bottom waters, as shown by similar facies and diatom assemblages as before, though detrital layers were coarser, (3) from 1979 to 1994, sediment delivery to the lake increased and laminated, clastic facies were deposited, and (4) from 1995 to 2008, dominance of massive facies and an increase in Fragilaria tenera and Achnanthes minutissima reflect relatively lower lake levels, less frequent bottom anoxia with more frequent water column mixing, similar to modern conditions. The period 1952–1979 was a time of meromixis and varved facies deposition, and was characterized by higher rainfall and less intense agricultural pressure in the watershed. There were two short humid periods (1992–1993 and 1996–1998) when monitoring data show more anoxic weeks per year and relatively higher lake levels. Increased cultivation of small landholdings in 1963, and particularly after 1979, caused a large increase in sediment delivery to the lake. The inferred lake evolution is in agreement with monitoring data that suggest a transition from dominantly meromictic conditions prior to 1993–1994 to a predominantly monomictic pattern of circulation since then, particularly after 2000. The synergistic effects of intensive water extraction for irrigation and lower rainfall since 1979, and particularly since 1994, brought the long period of meromictic conditions in Lake Arreo to an end. Water balance and sediment delivery to the lake are dominant factors that control the limnological and mixing conditions in Lake Arreo and they must be considered in management and restoration plans.     

Colloidal Stability and Toxicity of Gold Nanoparticles and Gold Chloride on <Emphasis Type="Italic">Chlamydomonas reinhardtii</Emphasis>

Here we have examined interactions of gold nanoparticles differing in primary particle size and coating with the green algae Chlamydomonas reinhardtii as function of the colloidal stability of the particles in the experimental media used for toxicity studies. Interactions of dissolved Au³⁺ ions with algae were also examined. Included endpoints were photosynthetic yield and algal growth. Morphological and structural effects were examined microscopically and by flow cytometry. The results indicate no significant toxicity of gold nanoparticles to C. reinhardtii. Analysis of published data suggests toxicity of gold nanoparticles on algal growth to relate rather to particular coatings than to the gold core.     

Detecting red tides in the eastern Seto inland sea with satellite ocean color imagery

A multi-spectral classification scheme is proposed to identify water with red tide(s) using satellite ocean color imagery obtained by the Sea-viewing Wide Field-of-view Sensor (SeaWiFS). The study area was the eastern Seto Inland Sea in Japan, where serious red tides frequently occur. “Background Ocean Colors” (BOCs hereafter), or colors of water around a red tide or those of the water before/after a red tide, are calculated as the monthly climatological average of normalized water-leaving radiances (nLw) with 0.01 degree spatial resolution with SeaWiFS imagery. Criteria for detecting red-tide pixels are established from analyses of characteristics of the nLws (in the 443, 490, 510, and 555 nm bands) anomalies from BOCs and the nLw spectra together with the red-tide records in Osaka Bay. The proposed scheme can efficiently indicate the presence or absence of red tides for independent match-ups with 83% accuracy. Additional validations of specific events indicate that the algorithm performed well in the study area. These results suggest that the scheme is appropriate to detect red tides in the optically complex coastal water of the eastern Seto Inland Sea.