夏季青藏高原低涡的切向流场及波动特征分析

从大气动力学原理出发,将高原低涡视为受热源强迫的边界层内涡旋,建立了柱坐标下满足梯度风平衡的低涡控制方程组,分析高原低涡切向流场的基本特征.在此基础上,通过求解线性化涡旋模式,得出高原低涡中各类波动的频散关系及其特征,同时定性讨论了热力作用对混合波动的影响以及混合波动与高原低涡流场特征的联系.使用中尺度数值模式WRF分...     

Non-Fickian mass transport in fractured porous media

The paper provides an introduction to fundamental concepts of mathematical modeling of mass transport in fractured porous heterogeneous rocks. Keeping aside many important factors that can affect mass transport in subsurface, our main concern is the multi-scale character of the rock formation, which is constituted by porous domains dissected by the network of fractures. Taking into account the well-documented fact that porous rocks can be considered as a fractal medium and assuming that sizes of pores vary significantly (i.e. have different characteristic scales), the fractional-order differential equations that model the anomalous diffusive mass transport in such type of domains are derived and justified analytically. Analytical solutions of some particular problems of anomalous diffusion in the fractal media of various geometries are obtained. Extending this approach to more complex situation when diffusion is accompanied by advection, solute transport in a fractured porous medium is modeled by the advection-dispersion equation with fractional time derivative. In the case of confined fractured porous aquifer, accounting for anomalous non-Fickian diffusion in the surrounding rock mass, the adopted approach leads to introduction of an additional fractional time derivative in the equation for solute transport. The closed-form solutions for concentrations in the aquifer and surrounding rocks are obtained for the arbitrary time-dependent source of contamination located in the inlet of the aquifer. Based on these solutions, different regimes of contamination of the aquifers with different physical properties can be readily modeled and analyzed.     

Accuracy of numerical simulations of contaminant transport in heterogeneous aquifers: A comparative study

This work deals with a comparison of different numerical schemes for the simulation of contaminant transport in heterogeneous porous media. The numerical methods under consideration are Galerkin finite element (GFE), finite volume (FV), and mixed hybrid finite element (MHFE). Concerning the GFE we use linear and quadratic finite elements with and without upwind stabilization. Besides the classical MHFE a new and an upwind scheme are tested. We consider higher order finite volume schemes as well as two time discretization methods: backward Euler (BE) and the second order backward differentiation formula BDF (2). It is well known that numerical (or artificial) diffusion may cause large errors. Moreover, when the Péclet number is large, a numerical code without some stabilising techniques produces oscillating solutions. Upwind schemes increase the stability but show more numerical diffusion. In this paper we quantify the numerical diffusion for the different discretization schemes and its dependency on the Péclet number. We consider an academic example and a realistic simulation of solute transport in heterogeneous aquifer. In the latter case, the stochastic estimates used as reference were obtained with global random walk (GRW) simulations, free of numerical diffusion. The results presented can be used by researchers to test their numerical schemes and stabilization techniques for simulation of contaminant transport in groundwater.     

Microtopographic properties of sparse moss vegetation in the Antarctic polar desert

The effect of topography on moss vegetation is examined to clarify the processes that affect the colonization of polar deserts on continental Antarctica. Data on the presence of the mosses Bryum pseudotriquetrum and Pottia heimii, and relative altitude were recorded. The altitude measurements were used to infer the underlying topographical attributes of the substrate in the study plots. Specifically, the local distribution of moss plants was clarified using the topographical attributes to construct generalized linear mixed models (GLMMs). The models suggested that steep slopes and convex microhabitats within areas of concave general relief (at the plot scale 4 × 4 m) promoted the establishment of moss. This correspondence to general relief was more apparent for B. pseudotriquetrum than for P. heimii. Among the study plots, general relief was found to be an important determinant of the precise spatial distribution of B. pseudotriquetrum. The standard surface estimated using the robust methods presented in this study is shown to be more accurate for describing moss distribution than the prevailing least-squares method.     

Numerical studies of methane production from Class 1 gas hydrate accumulations enhanced with carbon dioxide injection

Class 1 gas hydrate accumulations are characterized by a permeable hydrate-bearing interval overlying a permeable interval with mobile gas, sandwiched between two impermeable intervals. Depressurization-induced dissociation is currently the favored technology for producing gas from Class 1 gas hydrate accumulations. The depressurization production technology requires heat transfer from the surrounding environment to sustain dissociation as the temperature drops toward the hydrate equilibrium point and leaves the reservoir void of gas hydrate. Production of gas hydrate accumulations by exchanging carbon dioxide with methane in the clathrate structure has been demonstrated in laboratory experiments and proposed as a field-scale technology. The carbon dioxide exchange technology has the potential for yielding higher production rates and mechanically stabilizing the reservoir by maintaining hydrate saturations. We used numerical simulation to investigate the advantages and disadvantages of using carbon dioxide injection to enhance the production of methane from Class 1 gas hydrate accumulations. Numerical simulations in this study were primarily concerned with the mechanisms and approaches of carbon dioxide injection to investigate whether methane production could be enhanced through this approach. To avoid excessive simulation execution times, a five-spot well pattern with a 500-m well spacing was approximated using a two-dimensional domain having well boundaries on the vertical sides and impermeable boundaries on the horizontal sides. Impermeable over- and under burden were included to account for heat transfer into the production interval. Simulation results indicate that low injection pressures can be used to reduce secondary hydrate formation and that direct contact of injected carbon dioxide with the methane hydrate present in the formation is limited due to bypass through the higher permeability gas zone.     

位场各阶垂向导数的ISVD算法及其应用

位场垂向导数广泛应用于位场数据处理与解释当中,目前求取位场各阶垂向导数方法大致可以分为两类：①在频率域进行;②在空间域进行。针对这两类方法存在着稳定性较差的问题,这里介绍了计算位场各阶垂向导数的ISVD算法,并通过对单一地质体模型和地质体组合模型产生的布格重力异常进行处理分析,证明了ISVD算法在计算位场各阶垂向导数时,特别是在计算位场高阶垂向导数时,比在频率域及空间域计算结果具有较好的稳定性,ISVD算法能够有效地识别出浅层小尺度地质体边界。为了验证该方法对实际资料的处理效果,对某区块实际布格重力异常数据应用ISVD算法试算各阶垂向导数,结合地质背景等资料划分出断裂体系。结果证明,应用ISVD算法计算的位场各阶垂向导数,对于确定浅层小尺度地质体边界有较好的效果。     

徐州大北望地区馒头组碳酸盐-陆源碎屑混合沉积

陆源碎屑与碳酸盐混合沉积是一种沉积机理特殊而又有重要意义的沉积现象,在徐州大北望地区寒武系馒头组广泛发育。在野外观察和地质资料分析的基础上,结合薄片分析,对徐州大北望地区的混合沉积进行了详细的研究。结果表明:发育混合沉积的环境有混积潮坪、混积局限台地、混积开阔台地、混积台地边缘。在此基础上,以沉积学原理为指导,建立了徐州大北望地区馒头组混合沉积模式,并对研究区混合沉积的控制因素进行了归纳和讨论。     

Identification and quantification of mixed sources of oil spills based on distributions and isotope profiles of long-chain n-alkanes

Combined with quantitative determination of concentration and isotopic composition of petroleum hydrocarbons, weathering simulation experiments on artificially mixed oils and their two end-member oils are performed for identification and quantification of mixed sources. The >C₁₈n-alkanes show no appreciable losses during a short-term weathering process. An approach based on distribution of long-chain n-alkanes (>C₁₈) is suggested for estimating the contribution proportion of each source in mixed oils. Stable carbon isotope profile of individual n-alkanes is a powerful tool to differentiate sources of oil spills, but unavailable to accurately allocate each contribution due to a relatively large analytical error.     

Numerical modeling of two-phase flow in heterogeneous permeable media with different capillarity pressures

Contrast in capillary pressure of heterogeneous permeable media can have a significant effect on the flow path in two-phase immiscible flow. Very little work has appeared on the subject of capillary heterogeneity despite the fact that in certain cases it may be as important as permeability heterogeneity. The discontinuity in saturation as a result of capillary continuity, and in some cases capillary discontinuity may arise from contrast in capillary pressure functions in heterogeneous permeable media leading to complications in numerical modeling. There are also other challenges for accurate numerical modeling due to distorted unstructured grids because of the grid orientation and numerical dispersion effects. Limited attempts have been made in the literature to assess the accuracy of fluid flow modeling in heterogeneous permeable media with capillarity heterogeneity. The basic mixed finite element (MFE) framework is a superior method for accurate flux calculation in heterogeneous media in comparison to the conventional finite difference and finite volume approaches. However, a deficiency in the MFE from the direct use of fractional flow formulation has been recognized lately in application to flow in permeable media with capillary heterogeneity. In this work, we propose a new consistent formulation in 3D in which the total velocity is expressed in terms of the wetting-phase potential gradient and the capillary potential gradient. In our formulation, the coefficient of the wetting potential gradient is in terms of the total mobility which is smoother than the wetting mobility. We combine the MFE and discontinuous Galerkin (DG) methods to solve the pressure equation and the saturation equation, respectively. Our numerical model is verified with 1D analytical solutions in homogeneous and heterogeneous media. We also present 2D examples to demonstrate the significance of capillary heterogeneity in flow, and a 3D example to demonstrate the negligible effect of distorted meshes on the numerical solution in our proposed algorithm.     

A finite element solution for the fractional advection–dispersion equation

The fractional advection–dispersion equation (FADE) known as its non-local dispersion, has been proven to be a promising tool to simulate anomalous solute transport in groundwater. We present an unconditionally stable finite element (FEM) approach to solve the one-dimensional FADE based on the Caputo definition of the fractional derivative with considering its singularity at the boundaries. The stability and accuracy of the FEM solution is verified against the analytical solution, and the sensitivity of the FEM solution to the fractional order α and the skewness parameter β is analyzed. We find that the proposed numerical approach converge to the numerical solution of the advection–dispersion equation (ADE) as the fractional order α equals 2. The problem caused by using the first- or third-kind boundary with an integral-order derivative at the inlet is remedied by using the third-kind boundary with a fractional-order derivative there. The problems for concentration estimation at boundaries caused by the singularity of the fractional derivative can be solved by using the concept of transition probability conservation. The FEM solution of this study has smaller numerical dispersion than that of the FD solution by Meerschaert and Tadjeran (J Comput Appl Math 2004). For a given α, the spatial distribution of concentration exhibits a symmetric non-Fickian behavior when β = 0. The spatial distribution of concentration shows a Fickian behavior on the left-hand side of the spatial domain and a notable non-Fickian behavior on the right-hand side of the spatial domain when β = 1, whereas when β = −1 the spatial distribution of concentration is the opposite of that of β = 1. Finally, the numerical approach is applied to simulate the atrazine transport in a saturated soil column and the results indicat that the FEM solution of the FADE could better simulate the atrazine transport process than that of the ADE, especially at the tail of the breakthrough curves.