Methodology of hydrogeological characterization of deep carbonate aquifers as potential reservoirs of groundwater. Case of study: the Jurassic aquifer of <Emphasis Type="Italic">El Maestrazgo</Emphasis> (Castellón,Spain)

A methodology for the characterization of deep carbonate aquifers has been developed and applied to El Maestrazgo Jurassic aquifer in Castellón, Spain. Characterization of these aquifer formations, located at more than 300 m deep, consisted of a previous phase of compilation, analysis and synthesis of the existing information about the area, followed by a coordinated combination of different speciality studies: geology, stratigraphy, structural analysis, hydrogeology, hydrochemistry, geophysics and remote sensing. Geological studies included geological mapping, definition of stratigraphical units and facies and structural analysis. The aim of the hydrogeology study was to define aquifer formations, recharge area, aquifer points inventory and groundwater flow directions for the establishment of piezometric and water quality observation nets. Special techniques were applied, like thermal infrared aerial images and the evaluation of submarine groundwater discharge by means of natural radium isotopes. Hydrochemical techniques, including majority elements characterization and stable isotopes (¹⁸O, ²H and ³H) determination, allowed classifying hydrochemical facies and establishing a renewal pattern for water within the system. Geophysics was useful in determining the aquifer geometry, the features of the basement and the petrophysical characteristics of the geological formations. Preliminary results show an important tectonic complexity and the possibilities for groundwater uses in the area of study.     

影响低渗透油藏低速非线性渗流的实验研究

以膨胀性的粘土为主要研究对象,通过粘土浆体的流变性实验和填砂管渗流实验,并使用一种新研制的微流量测量工具,研究了低渗岩心在较低流速下的渗流特性;利用微尺度流动理论,Einstein粘度定律等进行了理论分析,对引起低渗透油藏低速非线性渗流的影响因素进行了初步的探索.结论认为,粘土的膨胀性以及其水化后形成浆体的流变性是引起低速非线性渗流的主要因素之一.     

Isotopic variation of fishes in freshwater and estuarine zones of a large subtropical coastal lagoon

We used stable C and N isotope ratios of tissues from 29 fish species from a large subtropical lagoon in southern Brazil to examine spatial variability in isotopic composition and vertical trophic structure across freshwater and estuarine habitats. Nitrogen isotope ratios indicated a smooth gradation in trophic positions among species, with most fishes occupying the secondary and tertiary consumer level. Fish assemblages showed a significant shift in their carbon isotopic signatures between freshwater and estuarine sites. Depleted carbon signatures (from −24.7‰ to −17.8‰) were found in freshwater, whereas more enriched signatures (from −19.1‰ to −12.3‰) were obtained within the estuarine zone downstream. Based on our survey of the C₃ and C₄ plants and isotopic values for phytoplankton and benthic microalgae reported for ecosystems elsewhere, we hypothesized that the observed δ¹³C differences in the fish assemblage between freshwater and estuarine sites is due to a shift from assimilating organic matter ultimately derived from C₃ freshwater marsh vegetation and phytoplankton at the freshwater site (δ¹³C ranging from −25‰ to −19‰), to C₄ salt-marsh (e.g. Spartina) and widgeon grass (Ruppia maritima), benthic microalgae and marine phytoplankton at the estuarine sites (from −18‰ to −12‰). Our results suggested that fish assemblages are generally supported by autochthonous primary production. Freshwater fishes that likely were displaced downstream into the estuary during periods of high freshwater discharge had depleted δ¹³C values that were characteristic of the upper lagoon. These results suggest that spatial foodweb subsidies can occur within the lagoon.     

Numerical Experiment on Two-Dimensional Line Thermal

The time evolution of a two-dimensional line thermal-a turbulent flow produced by an initial element with signifi-cant buoyancy released in a large water body, is numerically studied with the two-equation k - e model for turbulence closure. The numerical results show that the thermal is characterized by a vortex pair flow and a kidney shaped concentra-tion structure with double peak maxima; the computed flow details and scalar mixing characteristics can be described by self-similar relations beyond a dimensionless time around 10. There are two regions in the flow field of a line thermal: a mixing region where the concentration of tracer fluid is high and the flow is turbulent and rotational with a pair of vortex eyes, and an ambient region where the concentration is zero and the flow is potential and well-described by a model of doublet with strength very close to those given by early experimental and analytical studies. The added virtual mass coeffi-cient of the thermal motion is found to be approximat     

秦山核电站邻近水域的流况分析

本文分析了秦山核电站邻近水域各测站的实测流、潮流和佘流分布及季节变化。结果表明:测区内潮流为非正规丰日浅海潮流,其运动形式呈往复流。核电站二期工程取水口附近(L_1,18站)为一大流速区,而一期工程出水口周围(L_2,L_2′站)有一相对小的流速区。大流速区落潮流历时大于涨潮流历时,而相对小流速区则相反;大流速区的余流较小,而小流速区的余流大。在四个季节中,涨潮流速秋季最大,春季次之,夏季最小;落潮流速秋季最大,夏季次之,春季最小。余流冬、秋两季较大。在实测流中以潮流为主,径流和季风的作用也相当显著。     

How tides and river flows determine estuarine bathymetries

For strongly tidal, funnel-shaped estuaries, we examine how tides and river flows determine size and shape. We also consider how long it takes for bathymetric adjustment, both to determine whether present-day bathymetry reflects prevailing forcing and how rapidly changes might occur under future forcing scenarios.Starting with the assumption of a 'synchronous' estuary (i.e., where the sea surface slope resulting from the axial gradient in phase of tidal elevation significantly exceeds the gradient in tidal amplitude ), an expression is derived for the slope of the sea bed. Thence, by integration we derive expressions for the axial depth profile and estuarine length, L, as a function of and D, the prescribed depth at the mouth. Calculated values of L are broadly consistent with observations. The synchronous estuary approach enables a number of dynamical parameters to be directly calculated and conveniently illustrated as functions of and D, namely: current amplitude Û, ratio of friction to inertia terms, estuarine length, stratification, saline intrusion length, flushing time, mean suspended sediment concentration and sediment in-fill times.Four separate derivations for the length of saline intrusion, L_I, all indicate a dependency on (U_o is the residual river flow velocity and f is the bed friction coefficient). Likely bathymetries for `mixed' estuaries can be delineated by mapping, against and D, the conditions L<sub>I</sub>/L<1,E<sub>X</sub>/L<1 (E<sub>X</sub> is the tidal excursion) alongside the Simpson-Hunter criteria D/U<sup>3</sup><50 m<sup>−2</sup> s<sup>3</sup>. This zone encompasses 24 out of 25 `randomly' selected UK estuaries.However, the length of saline intrusion in a funnel-shaped estuary is also sensitive to axial location. Observations suggest that this location corresponds to a minimum in landward intrusion of salt. By combining the derived expressions for L and L_I with this latter criterion, an expression is derived relating D_i, the depth at the centre of the intrusion, to the corresponding value of U_o. This expression indicates U_o is always close to 1 cm s⁻¹, as commonly observed. Converting from U_o to river flow, Q, provides a morphological expression linking estuarine depth to Q (with a small dependence on side slope gradients).These dynamical solutions are coupled with further generalised theory related to depth and time-mean, suspended sediment concentrations (as functions of and D). Then, by assuming the transport of fine marine sediments approximates that of a dissolved tracer, the rate of estuarine supply can be determined by combining these derived mean concentrations with estimates of flushing time, F_T, based on L_I. By further assuming that all such sediments are deposited, minimum times for these deposition rates to in-fill estuaries are determined. These times range from a decade for the shortest, shallowest estuaries to upwards of millennia in longer, deeper estuaries with smaller tidal ranges.     

河口环流和盐水入侵Ⅱ--径流量和海平面上升的影响

本文利用本序列上篇论文建立的理想河口数值模式,研究径流量和海平面上升对河口环流和盐水入侵的影响.在径流量增大的情况下,口门内表层向海的流速增大,底层向陆的密度流减弱,滞流点下移.口门外侧向口门的密度流增大,上升流趋于增强.口门内盐水入侵减弱,口外盐度减小、冲淡水扩展范围增大.在口门上游北岸底层盐度下降明显,口门处南岸表层盐度下降明显.径流量变化对盐水入侵影响十分巨大.在海平面上升的情况下,拦门沙区域向陆的密度流增强,滞流点上移,表层向海的流动增大.口门内盐水入侵增强,口外盐度增大,冲淡水扩展范围减小.海平面上升对盐水入侵影响十分明显,北岸底层盐度增大尤为特出.     

一个南海热带气旋发展的环流特征分析

本文用ECMWF2.5×2.5网格点资料通过一个例子,对有利于南海热带气旋发展的环流场作了详细的分析,发现几点事实:(1)扰动在风的垂直切变较大的环境中仍可发展。(2)扰动位于中、低空西南急流左侧,高空南支东风急流右侧时,有利发展。(3)中,低空西南急流形成了扰动环流的强风潮,强风潮从扰动南侧入角,并包围扰动的东半圆。(4)扰动中、低层的入流和高层的外流呈现严重的不对称性。(5)当扰动发展到一定强度时,对南海热带季风环流圈的加强起到正反馈作用。     

Fluid flow equations governing primary oil migration as a separate phase

A mathematical model of primary oil migration as a separate phase out of compacting shales is presented. During burial and oil generation, source rock porosity decreases and oil saturation increases until residual oil saturation is reached. At this stage oil is expelled out by capillary and excess fluid pressure gradients. The model is a system of differential equations which relate changes in oil and water saturation in time to water and oil flow out of the source rock during burial. An additional set of equations for periods of erosion of overburden are also provided. The equations can be numerically solved by finite difference method. If oil and water flow is to be simulated during oil generation, then at each time step, changes by oil generation in oil and water saturations and porosity must be calculated. The solution procedure is briefly outlined.     

Simulation of green water flow on deck using non-linear dam breaking theory

Experiments carried out with models of floating production, storage and offloading platforms (FPSOs) showed that the flow of water over the deck edge, onto the deck resembled a suddenly released wall of water rather than a breaking wave. Therefore green water flow onto the deck was simulated using dam breaking theory, but the theory’s shallow-water assumptions may be limiting. In this paper a non-linear dam breaking problem is formulated. Equations of motion in the Lagrangian form are used and the solution is sought as an infinite series in time. Comparisons with the shallow water approximation are carried out.