Groundwater recharge of carbonate aquifers of the Silesian-Cracow Triassic (southern Poland) under human impact

A Triassic carbonate unit has been intensively drained by zinc and lead ore mines and numerous borehole fields since the nineteenth century. Its groundwater recharge has increased due to: pumping of water from boreholes, mining activity, and urbanization. An approach to determine the amounts of the recharge at a variety of spatial scales is presented in the paper. Different methods were used to identify and quantify recharge components on a regional and local scale: mathematical modelling was performed for four aquifers included in an aquifer system, an analytical estimation based on the assumption that an average recharge is equal to the average discharge of the hydrogeological system—for six man-made drainage centres, and the method of water level fluctuation (WLF) was applied in one observation borehole. Results of modelling have been supplemented by observation of environmental tracers (δ¹⁸O, δ²H, ³H), noble gases temperatures, and ⁴He_exc in groundwater. The regional aquifer’s current recharge according to estimations performed by means of modelling varies from 39 to 101 mm/year on average. Depending on the aquifer site the average precipitation ranges from 779 to 864 mm/year. In the confined part of the aquifer average recharge ranges from 26 to 61 mm/year. Within outcrops average recharge varies from 96 to 370 mm/year. Current recharge estimated by the analytical method for man-made drainage centres varies from 158 up to 440 mm/year. High values are caused by different recharge sources like precipitation, induced leakage from shallow aquifers, and water losses from streams, water mains and sewer systems. Pumping of water, mining and municipal activities constitute additional factors accounting for the intensified recharge.     

雄性三倍体长牡蛎繁殖潜力的初步研究

国内外很多人做过三倍体长牡砺的诱导工作，但在三倍体性腺发育及配子活力方面的研究，目前仅见Ａｌｌｅｎ等１９８６年和１９９０年对三倍体的海螂和长蛎异常配子的发生报道。Ｃｏｘ１９９６年研究了悉尼岩牡蛎三倍体性腺发育情况，姜卫国等对三倍体合浦珠母贝的生殖腺进行了观察，李霞［１］对三倍体长牡蛎的性腺发育做了较系统的分析和总结。作者用流式细胞仪对长牡蛎Ｃｒａｓｓｏｓｔｒｅａ　ｇｉｇａｓ二、三倍体的性腺细胞进行了初步分析，回答了三倍体长牡蛎生殖潜力及性比方面的一些问题，为进一步研究三倍体长牡蛎的生长、繁殖、杂…     

Relative molar mass distributions of chromophoric colloidal organic matter in coastal seawater determined by Flow Field-Flow Fractionation with UV absorbance and fluorescence detection

A new method for the characterization of chromophoric colloidal organic matter in seawater has been applied to samples from the Baltic Sea, Kattegatt and Skagerrak seas. Size fractionation of the sample by Flow Field-Flow Fractionation and measurement of the fluorescent and UV absorbing properties of the individual size fractions result in a relative molar mass distribution (RMM) of the optical properties. The RMM distributions have been used to estimate number and weight average relative molar masses, and polydispersity indices. At least two sources of coloured organic matter were identified from the ratio of fluorescence to UV: the Baltic surface water and the Skagerrak deep water. The dominating processes were mixing and dilution, but processes such as photo bleaching of fluorescence are also believed to be important. The RMM distribution derived from UV detection (1150–1300 Dalton) increased with increasing salinity while that derived for fluorescence (1500–1250 Dalton) decreased with increasing salinity. The specific UV absorbance taken as a proxy of the aromaticity of the chromophoric organic material showed decreasing trend with both increasing salinity and increasing UV derived weight average relative molar mass. Increasing polydispersity of the colloidal material was also observed as a function of salinity.     

Standard facies models to incorporate all heterogeneity levels in a reservoir model

The construction of reservoir models is frustrated by the fact that core and well cover only a fraction of the reservoir volume and it is therefore difficult to determine features like facies shape, -size, and -distribution, inter- and intra-facies boundaries and lateral trends from them. These features are, however, critical to fluid flow and they should necessarily be incorporated in the reservoir model and we therefore propose to systematically describe geometry and distribution of facies. To this end we make use of “standard facies models” that a priori contain all elements and boundaries of facies for a number of typical depositional environments.     

Inverse reconstructions of ecosystem flows in investigating regime shifts: impact of the choice of objective function

Inverse analysis is increasingly used in ecosystem modelling to objectively reconstruct a large number of unknown flows or interactions from a small number of observations. This type of analysis may be useful in relating observed regime shifts in ecosystem structure to underlying processes. Inversions of ecosystem flow networks currently use a constrained least-squares solution which at the same time minimizes the squared norm (the sum of squares) of the reconstructed flows. This minimum norm (MN) inversion is thought to be a parsimonious solution to the ecosystem flow inverse problem, but it may well not reflect how ecosystems are organised. It has been proposed instead that ecosystems evolve to maximize energy/mass flows or that they maximize the information content of the network weighted by ecosystem flows (ascendancy). We used simulated inverse experiments, where inverse analyses are applied to simulations of flow networks, to explore objective functions different than the MN generally used. We could not compute inverse solutions that maximize ascendancy because the objective function is unbounded. We could calculate inversions that maximize flows; however, these generally overestimated the simulated flows, even though the simulations were designed to maximize flows. It appears that the ecosystem flow inverse problem is too under-determined (too few data relative to the number of unknowns) to allow the use of these maximizing goal functions. We introduce a new minimization that simultaneously minimizes the squared flows and the squared differences between flows. This smoothing minimization makes the inverse flows as even as possible and it helps with some technical issues with MN inversions. The simulated inverse experiments indicated that this smoothed norm (SM) is the most robust in comparative analyses of contrasting ecosystem states, such as those that can be associated with regime shifts. Like the MN inversion, the SM inversion has no ecological basis. However, it is a conservative norm that is less likely to produce false differences between the dynamics of regimes.     

A moving PIV system for ship model test in a towing tank

A moving particle image velocimetry (PIV) system was successfully developed and used in a large towing tank for ship model tests to observe velocity fields near ship models. The experimental method involved adjustable optical devices for various test conditions and a special particle-seeding device. The streamwise and cross-streamwise flow fields of a yacht model and a tanker model were measured. Ship type, bottom shape, and towing speed were found to be the causes of problems affecting optical access and image quality. Possible solutions, deeper optical ducts, dark painting color, and pre-processed analysis method, were proposed and discussed.     

The dimensions of sand ripples in full-scale oscillatory flows

New large-scale experiments have been carried out in two oscillatory flow tunnels to study ripple regime sand suspension and net sand transport processes in full-scale oscillatory flows. The paper focuses on ripple dimensions and the new data are combined with existing data to make a large dataset of ripple heights and lengths for flows with field-scale amplitudes and periods. A feature of the new experiments is a focus on the effect of flow irregularity. The combined dataset is analysed to examine the range of hydraulic conditions under which oscillatory flow ripples occur, to examine the effects of flow irregularity and ripple three-dimensionality on ripple dimensions and to test and improve existing methods for predicting ripple dimensions.The following are the main conclusions. (1) The highest velocities in a flow time-series play an important role in determining the type of bedform occurring in oscillatory flow. Bedform regime is well characterised by mobility number based on maximum velocity in the case of regular flow and based on the mean of the highest one tenth peak velocities in the case of irregular flow. (2) For field-scale flows, sand size is the primary factor determining whether equilibrium ripples will be 2D or 3D. 2D ripples occur when the sand D₅₀ ≥ 0.30 mm and 3D ripples occur when D₅₀ ≤ 0.22 mm (except when the flow orbital diameter is low). (3) Ripple type (2D or 3D) is the same for regular and irregular flows and ripple dimensions produced by equivalent regular and irregular flows follow a similar functional dependence on mobility number, with mobility number based on maximum velocity in the case of regular flow and based on the mean of the highest one tenth velocities in the case of irregular flow. For much of the ripple regime, ripple dimensions have weak dependency on mobility number and ripple dimensions are similar for regular and irregular flows with the same flow orbital amplitude. However, differences in ripples produced by equivalent regular and irregular flows become significant at the high mobility end of the ripple regime. (4) Ripple dimensions predicted using the Wiberg and Harris formulae are in poor agreement with measured ripple dimensions from the large-scale experiments. Predictions based on the Mogridge et al. and the Nielsen formulae show better overall agreement with the data but also show systematic differences in cases of 3D ripples and ripples generated by irregular flows. (5) Based on the combined large-scale data, modifications to the Nielsen ripple dimension equations are proposed for the heights and lengths of 2D ripples. The same equations apply to regular and irregular flows, but with mobility number appropriately defined. 3D ripples are generally smaller than 2D ripples and estimates of 3D ripple height and length may be obtained by applying multipliers of 0.55 and 0.73 respectively to the 2D formulae. The proposed modified Nielsen formulae provide an improved fit to the large-scale data, accounting for flow irregularity and ripple three-dimensionality.     

泥浆性能对定向钻进铺管孔壁稳定性的影响

本文通过对泥浆性能的分析,特别是对泥浆流变参数的控制与调节在定向钻进非开挖铺管中的作用及其护壁原理的分析,简要阐述了泥浆性能对定向钻进铺管孔壁稳定性的影响,对非开挖定向钻进铺管施工有一定的指导意义。     

关于城市交通地理信息系统模型设计的思考与改进

GIS—T模型是建立城市交通地理信息系统的主体,本文对现有GIS—T模型进行了简单的分析与评价,并且为了增加政府与公众做出决策的正确性,提出了路径畅通可靠度模型、城市交通综合指数模型、交通出行指数模型三个数学模型,以达到量化城市交通与完善GIS—T模型的目的。     

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.