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.     

天文(重力)水准的新手段及其计算方法

介绍当今确定天文大地垂偏差的新仪器－CCD自动天体测量仪和确定两点间高程异常差的新方法，该方法用单极坐标代替繁琐的双极坐标进行计算，讨论了天文重力水准的误差及垂线偏差非线性影响等问题，估计在不久的将来，用这一手段施测山区似大地水准面的精度可望达到厘米级。     

地壳演化与成矿作用——以川滇地洼系“四层楼”铜矿床序列为例

提出川滇地洼系“四层楼”铜矿床序列的形成与陆壳演化的成生联系,是与本区陆壳由前地槽—地槽—地台—地洼演化各阶段与之相匹配的成矿作用的产物.与此同时,并总结了本区“四层楼”铜矿床序列的成矿作用具有明显的继承性、新生性、旋回性及层控性四大特点和多因复成矿床的成矿模式.     

用随机介质模型方法描述孔洞型油气储层

孔洞型储集层具有强烈的非均质性,在空间分布上变化剧烈,规律性差,其地球物理参数(如速度、密度、弹性参数等)在空间上的变化很难用传统的层状介质模型进行刻画.本文借助随机模型方法,研究孔洞型储层中不同孔洞尺度、分布密度与模型参数的关系,建立反映实际孔洞油气储层介质地球物理参数空间分布统计特征的随机介质模型,以更方便地对孔洞型油气储层进行描述,并为碳酸盐岩孔洞型储层地震反射特征的研究提供有效的方法途径.     

The mitigation of the vibration effects caused by tunnel blasts in urban areas: a case study in Istanbul

Ground vibrations generated by commercial explosives in tunnel construction may cause structural damage in urban areas. Therefore, suppressing the vibration effects and mitigating the possible hazard after blasting is important. We present a new method of controlled blasting that is environmentally friendly, and easy to utilize for tunnel construction. Small charges in this method are detonated sequentially to produce minimum side effects. The efficiency of the charges may be increased based on the previously monitored shots. This method is utilized in a tunnel construction in Istanbul with five experimental shots. In these experiments, the duration and also the quantity of explosives were carefully controlled. We were able to obtain better results with short durations (480 ms) instead of long durations (9,000 ms) although the vibration levels defined as peak particle velocity (PPV) became bigger while the quantity of the explosive charge increased from 3.088 to 9.264 kg.     

库车坳陷大北地区下白垩统巴什基奇克组储层成因地质分挝

塔里木盆地库车坳陷大北地区白垩系巴什基奇克组储层以细一中粒岩屑砂岩为主夹少量长石岩屑砂岩,杂基为铁泥质和泥质(2%～10%).胶结物以方解石(3%～15%)为主.粘土矿物以伊利石、绿泥石和伊/蒙混层为主,不含高岭石,伊/蒙混层中的蒙皂石含量一般为(15～20)%;主要孔隙组合为残余原生粒问孔-溶蚀孔-微孔隙,占储集空间总量的(50～90)%,其次为构造缝-溶蚀孔,主要的孔隙类型为溶蚀孔和裂缝,总体属于低孔低渗-特低孔特低渗储层.综合分析认为:阶段性的前陆构造逆冲推覆作用对沉积旋回性、水体盐度变化及物源区距离变化的控制至关重要;构造活跃期扇三角洲前缘环境中形成的差分选中-粗砂岩、近源快速堆积造成的高岩屑含量、干旱成湖导致的高碳酸盐胶结和早成岩期相对深埋-后期持续快速深埋压实作用,是特低孔特低渗储层的主要成因;而较优质储层足构造平缓期辫状河三角洲前缘环境中形成的分选好的中-细砂岩、早期表生溶蚀及弱碳酸盐胶结作用、早期构造破裂的复合作用及其叠加效应的结果.     

多传感器观测下带乘性噪声系统的逆向滤波与反褶积融合算法

针对多传感器观测环境下带乘性噪声系统的逆向最优滤波与反褶积融合估计问题 ,本文提出了 1种基于极大似然准则的最优融合算法。该算法中各单传感器间并行计算 ,并且融合中心与单传感器处理中心间无反向通讯 ,因而执行效率较高。仿真表明 ,该融合算法产生的逆向滤波与反褶积比单传感器处理结果有较明显提高     

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

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