Identification of diffusion parameters in a nonlinear convection–diffusion equation using the augmented Lagrangian method

Numerical identification of diffusion parameters in a nonlinear convection–diffusion equation is studied. This partial differential equation arises as the saturation equation in the fractional flow formulation of the two-phase porous media flow equations. The forward problem is discretized with the finite difference method, and the identification problem is formulated as a constrained minimization problem. We utilize the augmented Lagrangian method and transform the minimization problem into a coupled system of nonlinear algebraic equations, which is solved efficiently with the nonlinear conjugate gradient method. Numerical experiments are presented and discussed. This work was partially supported by the Research Council of Norway (NFR), under grant 128224/431.     

基于宽巷载波的GPS周跳探测与修复的算法改进

鉴于双频P码伪距结合相位观测值用于探测周跳存在的不足,提出了结合DCPC图像和相应的算法来分析周跳发生大小以及发生位置的方法。尤其是在载波L1与L2出现大小相同的周跳且宽巷载波失效的情况下,需结合DCPC图像来分析周跳发生的历元,通过M-W宽巷载波组合及其电离层残差组合,用两历元间的二元一次方程来修复周跳。该方法能够快速地探测及修复周跳,具有实用性。     

GPS软件接收机精确载波频率捕获算法分析与实现

采用基于快速傅里叶变换（FFT）的并行码相位搜索捕获算法实现GPS信号的捕获,根据相位关系对载波频率进行精化,使其满足跟踪模块的需求。在分析算法的基础上,通过对模拟信号和实测数据进行捕获实验,实现了对载波频率的精化,获得精确频率可直接用于后续跟踪环路。     

DuMu: DUNE for multi-{phase, component, scale, physics, …} flow and transport in porous media

DuMu^x is a free and open-source simulator for flow and transport processes in porous media, based on the Distributed and Unified Numerics Environment (DUNE). Its main intention is to provide a sustainable and consistent framework for the implementation and application of model concepts, constitutive relations, discretizations, and solvers. The paper provides an overview of DuMu^x with the focus on software-related aspects. Selected examples highlight the multi-scale and the parallel capabilities.     

Divergence of solutions to perturbation-based advection-dispersion moment equations

Moment equation methods are popular and powerful tools for modeling transport processes in randomly heterogeneous porous media, but the application of these methods to advection-dispersion equations often leads to erroneous oscillations. Perturbative methods, required to close systems of moment equations, become inaccurate for large perturbations; however, little quantitative theory exists for determining when this occurs for advection-dispersion equations. We consider three different methods (asymptotic approximation, Eulerian truncation, and iterative solution) for closing and solving advection-dispersion moment equations describing transport in stratified porous media with random permeability. We obtain approximate analytical expressions for time above which the asymptotic approximation to the mean diverges, in particular quantifying the impact that dispersion has on delaying—but not eliminating—divergence. We demonstrate that Eulerian truncation and iterative solution methods do not eliminate divergent behavior either. Our divergence criteria provide a priori estimates that signal a warning to the practitioner of stochastic advection-dispersion equations to carefully consider whether to apply perturbative approaches.     

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.     

湖泊沉积物的矿物组成、成因、环境指示及研究进展

湖泊沉积物是不同地质、气候、水文条件下各类碎屑、黏土、自生/生物成因矿物以及有机物质等的综合体。沉积矿物蕴含着丰富区域和全球环境演变信息,如湖水的化学组成、流域构造、气候、水文以及人类活动的相互作用等。相关信息可以赋存在矿物外部微形貌、内部微结构、化学组成、物理和化学性质、同位素组成、谱学特征、成因以及共生组合等方面。因此,湖泊科学的许多关键课题都离不开矿物学,特别是在利用湖泊沉积物进行区域及过去全球变化研究中,深入的机理研究归根到底都要涉及矿物学,如流域化学风化作用、粒度组成、生物壳体化学组成、测年材料的选择等。然而,由于湖泊沉积物中矿物的多源性、复杂性,如何有效提取和解译其中的环境信息,是一项长期困扰研究者的课题,湖泊沉积矿物学的研究往往被许多研究者所忽视,中国的相关研究也较为薄弱。笔者综述了湖泊沉积物中碎屑、黏土、自生/生物矿物的矿物组合、特征、成因在(古)环境反演中的作用及最新研究进展,提出除了继续加强对湖泊沉积物中矿物来源、成因和古环境示踪的深入研究以外,矿物相间的转变及其对湖水和孔隙水组成的响应、一些非晶质或隐晶质及低丰度矿物相在湖泊化学和动力学中的作用也是很有潜力的研究领域,最后提出了研究中存在的一些问题、面对的挑战以及对研究前景的展望。     

油藏条件下孔隙岩样毛管和电学性质研究

本文系统总结了国内外关于油藏条件下岩石物理基础实验研究的成果、现状及发展方向。油藏条件主要是指油藏实际的温度、压力、毛管及电性平衡、饱和顺序和润湿性质等。实验研究及现场资料均表明，油藏条件和地面实际室条件下饱和流体孔隙岩样毛管及电学性质有很大差惜同藏条件下毛管及电学性质实际研究对准确评价油气储量及促进岩石物理学基础理论发展具有重要意义。     

Nonequilibrium effects in models of three-phase flow in porous media

In this paper we extend to three-phase flow the nonequilibrium formalism proposed by Barenblatt and co-workers for two-phase porous media flow. The underlying idea is to include nonequilibrium effects by introducing a pair of effective water and gas saturations, which are linked to the actual saturations by a local evolution equation. We illustrate and analyze how nonequilibrium effects lead to qualitative and quantitative differences in the solution of the three-phase flow equations.     

Unsaturated water flow across soil aggregate contacts

Unsaturated water flow through soil aggregates is controlled by the contacts between aggregates. The contacts are highly conductive when wet and become bottle-necks for flow when drained. We postulate that the hydraulic conductivity of the contacts is in first place determined by the water-filled contact area. The objective of this study was to measure and model the water-filled contact area and to relate it to the conductivity of a series of aggregates. We performed microscopic tomography of an aggregate pair equilibrated at different water potentials. By means of image analysis and a morphological pore network model, the water-filled contact area was calculated. We found that the aggregate surface is rough and the contact region contains macropores which are rapidly drained. As a consequence the water-filled contact area dramatically decreases as the water potential is diminished. We modeled this process by describing the aggregates as spheres covered by much smaller spheres representing the roughness. The water-filled contact was analytically calculated from this model. Knowing the water-filled contact area we up-scale the hydraulic conductivity of a series of aggregates. This is calculated as the harmonic mean of the contact and aggregate conductivities. The contact conductivity is calculated from the water-filled contact area. Near saturation the conductivity of a series of aggregates is close to the conductivity of a single aggregate, and, when further drained, it rapidly decreases as the water-filled contact area. The model matches the experimental data well.