低孔渗砂岩孔隙结构分类方法研究

低孔渗砂岩储层孔隙结构复杂,在二维与三维空间中具有较强非均质性,对储集空间的测井响应有极大影响。而压汞曲线、核磁T₂谱、铸体薄片3种资料在一定程度上各自反映了孔隙结构的不同特征,本文分别提取3种资料的微观孔隙结构参数进行相关性分析,最终分别利用主成分分析-模糊均值聚类方法实现孔隙结构分类,其分类结果压汞曲线分形法所得一致。     

Numerical analysis of frost effects in porous media. Benefits and limits of the finite element poroelasticity formulation

The aim of this paper is to analyse the performance of a finite element formulation usable for predicting the mechanical consequences of frost effects on porous media. It considers the characteristics of porous media and how the frost action can be assessed. The problem is then separated into two parts: thermal and poromechanical calculations. The constitutive equations developed in the framework of poromechanics are presented and the implementation in a usual finite element poroelasticity formulation based on Zuber's method is adopted. An analysis of the time‐step influence on the convergence rate is given and leads us to propose a simple method in order to obtain objectivity of the finite element response and avoid over‐long calculations. Frost effect simulations are carried out on real porous media (two fired clays) as a case study. Although the experimental behaviour of the porous media subjected to frost action is in accordance with some observations, the calculated strains appear to be overestimated compared with measurements. The problem could be largely attributable to the difficulty of assessing permeability evolution during frost development. Copyright © 2011 John Wiley & Sons, Ltd.     

Influence of pore size and geometry on peat unsaturated hydraulic conductivity computed from 3D computed tomography image analysis

In organic soils, hydraulic conductivity is related to the degree of decomposition and soil compression, which reduce the effective pore diameter and consequently restrict water flow. This study investigates how the size distribution and geometry of air‐filled pores control the unsaturated hydraulic conductivity of peat soils using high‐resolution (45 µm) three‐dimensional (3D) X‐ray computed tomography (CT) and digital image processing of four peat sub‐samples from varying depths under a constant soil water pressure head. Pore structure and configuration in peat were found to be irregular, with volume and cross‐sectional area showing fractal behaviour that suggests pores having smaller values of the fractal dimension in deeper, more decomposed peat, have higher tortuosity and lower connectivity, which influences hydraulic conductivity. The image analysis showed that the large reduction of unsaturated hydraulic conductivity with depth is essentially controlled by air‐filled pore hydraulic radius, tortuosity, air‐filled pore density and the fractal dimension due to degree of decomposition and compression of the organic matter. The comparisons between unsaturated hydraulic conductivity computed from the air‐filled pore size and geometric distribution showed satisfactory agreement with direct measurements using the permeameter method. This understanding is important in characterizing peat properties and its heterogeneity for monitoring the progress of complex flow processes at the field scale in peatlands. Copyright © 2010 John Wiley & Sons, Ltd.     

对饱和粘性土渗透规律的新认识及其应用

通过实验，对饱和粘性土的渗透规律进行研究，得到几点新认识。利用此认识建立了饱和粘性土越流补给量的计算公式，并在实际应用中得到验证。     

Dual models of pore spaces

We present a model for pore spaces that consists of two parts related by duality: (1) a decomposition of an open polyhedral pore space into open contractible pore bodies separated by relatively open interfaces and (2) a pore network that is homotopy equivalent to the pore space. The dual model is unique and free of parameters, but it relies on regularity conditions for the pore space. We show how to approximate any pore space by the interior of a polyhedral complex such that the regularity conditions are fulfilled. Thus, we are able to calculate the dual model from synthetic porous media and images of real porous media. The pore bodies are unions of relatively open Delaunay cells with respect to the corners of the pore boundary, and the pore network consists of certain at most two-dimensional (2D) Voronoi cells with respect to the corners of the pore boundary. The pore network describes the neighborhood relations between the pore bodies. In particular, any relatively open 2D Delaunay face f separating two pore bodies has a unique (relatively open) dual network edge. In our model, f is a pore throat only if it is hit by its dual network edge. Thus, as opposed to widespread intuition, any pore throat is convex, and adjacent pore bodies are not necessarily separated by pore throats. Due to the duality between the pore network and the decomposition of the pore space into pore bodies it is straightforward to store the geometrical properties of the pore bodies [pore throats] as attributes of the dual network vertices [edges]. Such an attributed network is used to perform 2D drainage simulations. The results agree very well with those from a pore-morphology based modeling approach performed directly on the digital image of a porous medium. Contractibility of the pore bodies and homotopy equivalence of the pore space and the pore network is proven using discrete Morse theory and the nerve theorem from combinatorial topology.     

远区重力数据对扰动位径向二阶梯度的影响

提出了利用地面重力异常数据计算地面扰动位径向二阶梯度时，远距离带区影响的误差公式，结果以级数形式给出，以其数学模型为基础，利用全球重力位模型GPM98，计算确定了不同平均异常网格大小情况下截断半径取一系列值时的远区影响，在将卫星重力梯度和航空重力梯度观测数据与地面重力数据进行联合处理时，文中提出了远区影响的数量关系，对于确定合理的截断半径大小具有一定的参考价值。     

双侧破裂模型对地震动空间相关性影响的定量研究

本在震源为双侧线状断层破裂模型的假定下，根据震源理论和波谱随空间坐标的展开式，得到了场地两点地震动Fourier谱的表达式Ab(f)=Aa(f)exp(a1δ0 a2δθ，中提出了将双侧破裂模型分解为两个单侧破裂模型计算双侧破裂模型a1和a2的数值方法，以一次地震（M＝6）为例，分析了双侧破裂模型对地震动空间相关性的影响，计算了 a1和a2的实部和虚部以及孔径随震中距，方位角和频率变化的三维图像。     

Streambank dewatering for increased stability

Streambank erosion is often the dominant source of sediment leaving modified watersheds. Mass failure of high, steep banks is one of the most serious forms of streambank erosion. The risk of a given bank experiencing mass failure is a function of bank height, angle, and soil strength, which is governed by soil moisture. Two methods for bank dewatering were tested in adjacent sections of streambank bordering a deeply incised channel in northern Mississippi: a low‐cost pump system and subsurface horizontal drains. Pore water pressures (both positive and negative pressures, or matric suction) were continuously monitored for 2 years at the pumped site, at an adjacent untreated control section, and for 1 year at the site stabilized with horizontal drains. Resulting data were used to calculate a time series of the factor of safety using a computer model. Over the course of two wet seasons, average bank retreats for the control and pumped plots were 0·43 and 0·21 m, respectively. More limited monitoring revealed that the site with passive drains retreated about 0·23 m. At the pumped site pore water pressure was 3–4 kPa lower than at the control site during the most critical periods. Accordingly, computed factors of safety were above the failure threshold at the pumped site, but fell below unity at the control site on 11 occasions over the period of observation. Similarly, the drained site displayed generally lower pore water pressure and higher safety factors except for two events when drains were evidently overwhelmed with the volume of local surface and subsurface flows. These results suggest, but do not prove, that bank dewatering promoted lower rates of bank retreat and higher levels of stability since the three sites had slight differences in soils, geometry and boundary conditions. Initial cost of the dewatering treatments were significantly less than orthodox bank stabilization measures, but operation and maintenance requirements may be greater. Published in 2009 by John Wiley & Sons, Ltd.