Three dimensional simulation of fluid flow in X‐ray CT images of porous media

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical manifold method for dynamic nonlinear analysis of saturated porous media

It is well known that the Babuska–Brezzi stability criterion or the Zienkiewicz–Taylor patch test precludes the use of the finite elements with the same low order of interpolation for displacement and pore pressure in the nearly incompressible and undrained cases, unless some stabilization techniques are introduced for dynamic analysis of saturated porous medium where coupling occurs between the displacement of solid skeleton and pore pressure. The numerical manifold method (NMM), where the interpolation of displacement and pressure can be determined independently in an element for the solution of u–p formulation, is derived based on triangular mesh for the requirement of high accurate calculations from practical applications in the dynamic analysis of saturated porous materials. The matrices of equilibrium equations for the second‐order displacement and the first‐order pressure manifold method are given in detail for program coding. By close comparison with widely used finite element method, the NMM presents good stability for the coupling problems, particularly in the nearly incompressible and undrained cases. Numerical examples are given to illustrate the validity and stability of the manifold element developed. Copyright © 2006 John Wiley & Sons, Ltd.     

滇西动态重力网的分形特征及空间分辨力

介绍了用分形学描述地球物理观测网非均匀分布和联合应用分形特征与Ｓｈａｎｎｏｎ定理确定最佳网格距的方法。研究了滇西动态重力网的分形特征，它的分形维Ｄｆ＝１．２７，最佳网格距ｒ０＝２８ｋｍ。讨论了１９８６年至１９９４年间网区发生的８次Ｍｓ＞５．０级地震前重力场变化的可信性，指出这些变化的波长均大于１００ｋｍ，至少２倍于网的最小空间分辨力，因而是可信的     

About Darcy's law in non-Galilean frame

This paper is aimed towards investigating the filtration law of an incompressible viscous Newtonian fluid through a rigid non-inertial porous medium (e.g. a porous medium placed in a centrifuge basket). The filtration law is obtained by upscaling the flow equations at the pore scale. The upscaling technique is the homogenization method of multiple scale expansions which rigorously gives the macroscopic behaviour and the effective properties without any prerequisite on the form of the macroscopic equations. The derived filtration law is similar to Darcy's law, but the tensor of permeability presents the following remarkable properties: it depends upon the angular velocity of the porous matrix, it verifies Hall–Onsager's relationship and it is a non-symmetric tensor. We thus deduce that, under rotation, an isotropic porous medium leads to a non-isotropic effective permeability. In this paper, we present the results of numerical simulations of the flow through rotating porous media. This allows us to highlight the deviations of the flow due to Coriolis effects at both the microscopic scale (i.e. the pore scale), and the macroscopic scale (i.e. the sample scale). The above results confirm that for an isotropic medium, phenomenological laws already proposed in the literature fails at reproducing three-dimensional Coriolis effects in all types of pores geometry. We show that Coriolis effects may lead to significant variations of the permeability measured during centrifuge tests when the inverse Ekman number Ek⁻¹ is 𝒪(1). These variations are estimated to be less than 5% if Ek⁻¹<0.2, which is the case of classical geotechnical centrifuge tests. We finally conclude by showing that available experimental data from tests carried out in centrifuges are not sufficient to determining the effective tensor of permeability of rotating porous media. Copyright © 2004 John Wiley & Sons, Ltd.     

地球化学景观吸引子的发现及其意义

勘查地球化学对异常评价的最大不确定性源于对地球化学景观特性的认识不够。我们应用构造相空间的方法分析了浙江儒岙幅（1:5万）水系沉积物测量的Ag含量数据，揭示了地球化学景观是具有低维吸引子（DZ≈2．9）的混饨系统。这一发现将对勘查地球化学中采样密度的选择、地球化学场模拟、异常评价等根本性问题产生重大影响。     

On the practice of estimating fractal dimension

Coastlines epitomize deterministic fractals and fractal (Hausdorff-Besicovitch) dimensions; a divider [compass] method can be used to calculate fractal dimensions for these features. Noise models are used to develop another notion of fractals, a stochastic one. Spectral and variogram methods are used to estimate fractal dimensions for stochastic fractals. When estimating fractal dimension, the objective of the analysis must be consistent with the method chosen for fractal dimension calculation. Spectal and variogram methods yield fractal dimensions which indicate the similarity of the feature under study to noise (e.g., Brownian noise). A divider measurement method yields a fractal dimension which is a measure of complexity of shape.     

Influence of Structural Non-Stationarity of Surface Roughness on Morphological Characterization and Mechanical Deformation of Rock Joints

Summary Structural non-stationarity of surface roughness affects accurate morphological characterization as well as mechanical behaviour of rock joints at the laboratory scale using samples with a size below the stationarity threshold. In this paper, the effect of structural non-stationarity of surface roughness is investigated by studying the scale dependence of surface roughness and mechanical behaviour of rock joints. The results show that the structural non-stationarity mainly affects the accurate characterization of the surface roughness of the fracture samples. It also controls the amount and location of the contact areas during shear tests, which in turn affects the mechanical properties and asperity degradation of the samples. It is concluded that for accurate determination of the morphological and mechanical properties of rock joints at laboratory and field scales, samples with size equal to or larger than the stationarity threshold are required. Author’s address: Nader Fardin, Rock Mechanics Group, Department of Mining Engineering, Faculty of Engineering, University of Tehran, P.O. Box: 11365/4563, Tehran, Iran     

多波束数据的海底数字地形模型构建

提出大批量多格式原始多波束数据的DTM构建方法，以满足大区域、大比例尺海底地形制图的需要。在多种格式原始多波束数据接口和系统内部标准数据结构的基础上，通过对数据文件、数据种类和数据运算量的有效组织管理，实施边读入、边权重配赋的网格插值，分析了高斯、指数和距离平方反比权重函数的适用性，及最小值、最大值和平均值的实用性。在权重配赋网格插值基础上，提出分形fBM和张力样条配合使用的方程式插值方法，保证DTM数据的有效外延和地形分辨。整套算法效率高，并能有效保证DTM的精度，对存储在外部介质的数据遍历一次即可完成网格插值。     

声波在水-多孔介质海底界面上的反射与透射

声波在海底界面的反射和透射是海底散射、海底混响、海底目标探测的重要问题。利用Biot多孔介质声传播理论对声波在水-多孔海底界面上的反射和透射进行了分析,具体给出反射声波的反射系数,3种透射声波的透射系数以及声能透射系数随入射波入射角和频率(10~40 kHz)的变化关系,分析了各种透射波对透射声能的贡献。多孔海底介质参数分别采用Stoll和Chotiros给出的2组参数进行理论计算。     

富钴结壳分布的分形特征

利用我国在西太平洋国际海底两个海山获得的富钴结壳资料进行分形研究后发现，富钴结壳的分布具有变维分形分布的特点，结壳丰度呈二阶累计和分形分布，结壳厚度呈多阶不定累计和分形分布，通过计算得到了它们各自的分维数，进而对两个海山的结壳资源进行了初步评价。