Density-dependent dispersion in heterogeneous porous media Part I: A numerical study

In this paper, we describe carefully conducted numerical experiments, in which a dense salt solution vertically displaces fresh water in a stable manner. The two-dimensional porous media are weakly heterogeneous at a small scale. The purpose of these simulations, conducted for a range of density differences, is to obtain accurate concentration profiles that can be used to validate nonlinear models for high-concentration-gradient dispersion. In this part we focus on convergence of the computations, in numerical and statistical sense, to ensure that the uncertainty in the results is small enough.Concentration variances are computed, which give estimates of the uncertainty in local concentration values. These local variations decrease with increasing density contrast. For tracer transport, obtained longitudinal dispersivities are in accordance with analytical findings. In the case of high-density contrasts, stabilizing gravity forces counteract the growth of dispersive fingers, decreasing the effective width of the transition zone. For small log-permeability variances, the decrease of the apparent dispersivity that is found is in agreement with laboratory results for homogeneous columns.     

Inversion of interference hydraulic pumping tests in both homogeneous and fractal dual media

This work proposes a complete method for automatic inversion of data from hydraulic interference pumping tests based on both homogeneous and fractal dual-medium approaches. The aim is to seek a new alternative concept able to interpret field data, identify macroscopic hydraulic parameters and therefore enhance the understanding of flow in porous fractured reservoirs. Because of its much contrasted sensitivities to parameters, the dual-medium approach yields an ill-posed inverse problem that requires a specific optimization procedure including the calculation of analytical sensitivities and their possible re-scaling. Once these constraints are fulfilled, the inversion proves accurate, provides unambiguous and reliable results. In the fractal context inverting several drawdown curves from different locations at the same time reveals more accurate. Finally, hydraulic parameters drawn from inversion should be taken into account to improve in various situations the conditioning of up-scaled flow in fractured rocks.     

西安地区黄土含水率变化特点

通过对西安地区黄土剖面的野外观察、颗粒分析和含水率的测定 ,系统地研究了黄土层含水率的变化特点 ,初步认识到在同一黄土剖面上 ,黄土含水率从上到下呈波动变化 ,其峰值位于古土壤层内 ,而波谷处于黄土层内 ,并呈波动增加的趋势。在具有结核层的黄土层中 ,其含水率从上部到下部呈由大变小再变大的趋势 ,其峰值位于结核层之上 ,CaCO3 结核层具有隔水作用 ;在古土壤层中 ,从上向下 ,其含水率呈由小变大再变小的规律 ,其峰值出现在古土壤的中、下部位 ,即粘土化最强的部位。土壤质地是决定含水率大小的重要因素 ,在包气带中 ,黄土含水率的变化与粘土颗粒含量具有明显的正相关关系。     

水与流体饱和孔隙介质界面上非镜面反射声场的实验研究

用光学法和发射-接收法研究了水与流体饱和孔隙介质界面上声波的非镜面反射声场．测量了入射角θ_i为16°－80°范围内的反射系数R，记录了θ_i在Rayleigh波临界角θ_c附近±4°范围内两束反射波的反射角、波束位移△、波形和频谱．实验发现：(1)在θ_c附近±4°范围内，观察到θ_i两束反射波A和B，其中波束B沿界面有波束位移△，实验测出△随θ_i增大而减小，并非在θ_i=θ_c时最大．(2)反射角等于θ_i的波束的频谱很接近于入射波的频谱，反射角等于θ_c的波束的频谱相对于入射波有频谱的负偏移．(3)波束A的相位与发射换能器直达波的相位一致，波束B的相位与波束A相反     

盆地超压层段非幕式突破期的地热场模型数值解法

盆地超压层段的地热场模拟对超压成因研究及油气生成，排放机理分析有重要意义，但求解超压层段的热传递方程的数值解法至今仍没有被很好地解决，问题的关键在于确定超压流体速度场，由于超压流体发生幕式突破之前的排出速率极低，其速度场v=vxi vzk近似满足条件Эux/Эx Эuz/Эz＝0，可将其近似地视为稳定的不可压缩的无源流体，利用这一条件及相应的边界压力条件，可使整个计算过程得到简化，在此基础上所建立的地热场模型有限元数值解法和模拟软件，能够实现对含油气盆地在幕式突破之前的超压层段地热场进行动态模拟。     

质疑粤西北的冰川遗迹

从地貌、沉积物、第四纪地层关系、古气候、古环境等方面论证了广东封开-怀集-带所谓的冰川谷、冰蚀三角面、悬谷、角峰、刃脊、冰半、鼻山尾、冰碛砾石、蛇形丘以及近百万年来发生过3次冰期，是既不符合地学基础理论也不符合野外实际，所以是不确实的，那里没有第四纪冰川遗迹。     

浅地层剖面仪在水域工程勘察中的应用

扼要介绍了水上浅地层剖面仪的方法技术，列举了该法在浅海，内河和湖泊减种典型环境及相应的地质条件下勘测所得到的剖面记录及地质解释结果，指出了运用该方法所能解决的工程地质问题。     

膨胀土作为桩基持力层的实验研究

根据对邯郸膨胀土中桩基的实验研究,认为膨胀土具有一定埋深时,可作为良好的桩基持力层;通过3种桩型应用情况的分析比较,论证了扩底桩是膨胀土持力层上理想的桩基形式。     

Spatial averaging of hydraulic conductivity in three-dimensional heterogeneous porous media

Numerical models of groundwater flow require the assignment of hydraulic conductivities to large grid blocks discretizing the flow domain; however, conductivity data is usually available only at the much smaller scale of core samples. This paper describes a geostatistical model for hydraulic conductivity at both the core or point scale and that of grid blocks. Conductivity at the block scale is obtained empirically as a spatial power-average of point scale values. Assuming a multivariate Gaussian model for point log-conductivity, expressions are derived for the ensemble mean and variance of block conductivity. The expression for the ensemble mean of block scale conductivity is found to be similar to an expression for the ensemble effective conductivity of an infinite field derived analytically by earlier authors. Here, block conductivities obtained by power averaging are compared with effective conductivities obtained from a numerical flow model and are found to be in excellent agreement for a suitably chosen averaging exponent. This agreement deteriorates gradually as the log variance of conductivity increases beyond 2. For arbitrary flow field geometry and anisotropic conductivity covariances, the averaging exponent can be calibrated by recourse to numerical flow experiments. For cubic fields and an isotropic spatial covariance, the averaging exponent is found to be 1/3. In this particular case, it was found that flow field discretization at the block scale through local averaging of point conductivities gave similar results to those obtained directly using a point scale discretization of the flow field.     

对穿时普遍性的解释

穿时是沉积地层的基本属性,其根本原因是地层的形成作用,即沉积作用本身。有沉积作用就有穿时(侧向加积);或者小范围内不穿时但大范围内必穿时(纵向加积)。海水进退只能影响穿时的具体形式而不是穿时的本质原因或必要条件。板块的运移和增生不仅造成了沉积地层的穿时,还造成了火山地层的穿时。这种穿时与沉积作用直接引起的穿时有本质不同,本文称之为构造穿时。