Montecarlo DLA-type simulations of wetting effects in fluid displacement in porous media

In this work, we report diffusion-limited aggregation (DLA)-type Montecarlo computations of a stochastic model of displacement of a viscous fluid by another that preferentially wets a porous medium, for the case when both fluids are immiscible in the absence of buoyancy forces. The model has the aim to simulate cooperative invasion processes found in experiments of immiscible wetting displacement. The model considers the nonlocal effects of the Laplacian pressure field and the capillary forces via hydrodynamic equations in the Darcy regime with a boundary condition for the pressure at the interface. The boundary condition contains two different types of disorder: the capillary term, which constitutes an additive random disorder, and a term containing an effective random surface tension, which couples to a curvature (it constitutes a multiplicative random term that carries nonlocal information of the whole pressure). We generate different displacement patterns for different setting of the parameters of the model. We analyze these patterns by studying the scaling properties of the interface that separate the two fluids and calculating the fractal dimension of the interface. The results show the existence of three distinct regimes of scaling. One regime at the smallest-length scales is due to the multiplicative random disorder together with the nonlocal coupling; it reveals itself in a roughness exponent α ≈ 0.80. Additionally, we find a DLA-type scaling regime with a roughness exponent α ≈ 0.60 at the largest scales and intermediate scaling regime with α ≈ 0.70 corresponding to invasion percolation with trapping. Each regime has definite scaling ranges that depend on the capillary number and the relative wetting tendency of the fluids. The behavior of the fractal dimensions of the interfaces of the aggregates constitutes a further confirmation of the existence of three scaling regimes and the multi-self-affinity of the perimeter of the interface boundaries.     

基于岩石物理实验的时移地震研究

时移地震技术是油田开发阶段动态检测的一个重要手段.岩石物理实验室能实现模拟油气田在开采过程中其储层参数变化对岩石弹性参数的影响.本文基于岩石物理实验及地震正演模拟,分析了X气田在开采过程中其孔隙压力下降及水不断侵入情况下,进行时移地震监测的可行性及其主要影响因素.     

地震尺度下碳酸盐岩储层的岩石物理建模方法(英文)

碳酸盐岩油藏的强非均质性以及孔隙结构的复杂性,使得作为连接油藏参数与地震参数重要桥梁的岩石物理模型,以及作为油藏预测和定量表征最有效工具的流体替换成为岩石物理建模的难点与重点。在碳酸盐岩储层复杂孔隙结构与地震尺度下碳酸盐岩储层非均质性分析基础上,研究采用岩石网格化方法,将地震尺度下非均质碳酸盐岩储层岩石划分为具有独立岩石参数的均质岩石子体,根据岩石孔隙成因与结构特征采用不同岩石物理模型分步计算岩石子块干岩石弹性模量,并根据不同孔隙连通性进行流体替换,计算饱和不同流体岩石弹性模量。基于计算的岩石子块弹性模量,采用Hashin-Shtrikman-Walpole弹性边界计算理论方法实现地震尺度下碳酸盐岩储层弹性参数计算。通过对含有不同类型孔隙组合碳酸盐岩储层模型的弹性模量进行计算与分析,明确不同孔隙对岩石弹性参数的影响特征,模拟分析结果与实际资料认识一致。     

Predicting the displacement of triple pendulum™ bearings in a full‐scale shaking experiment using a three‐dimensional element

The accuracy of a series spring model to predict the peak displacement and displacement history of Triple Pendulum? (TP) bearings in a strongly shaken, full‐scale building is evaluated in this paper. The series spring model was implemented as a self‐contained three‐dimensional TP bearing element in OpenSees and is now available for general use. The TP bearing element contains the option for constant friction or a generalized friction model that accounts for the effect of instantaneous velocity and compression load on the friction coefficient. Comparison between numerical simulation and experimental data of a five‐story steel moment frame building shows that the peak displacement of isolation system can generally be predicted with confidence using a constant friction coefficient model. The friction coefficient model accounting for the effect of axial load and velocity leads to minor improvement over the constant friction coefficient models in some cases. Copyright © 2013 John Wiley & Sons, Ltd.     

Determination of scoured bridge natural frequencies with soil–structure interaction

This study developed a finite element method with the effect of soil–fluid–structure interaction to calculate bridge natural frequencies. The finite element model includes bridge girders, piers, foundations, soil, and water. The effective mass above the soil surface was then used to find the first natural frequency in each direction. A field experiment was performed to validate that the natural frequencies calculated using the proposed finite element method had acceptable accuracy. The calculated natural frequencies with the fluid–structure interaction effect are always smaller than those without this effect. However, the frequency change due to the fluid effect is not obvious, so using the soil–structure interaction model is accurate enough in the bridge natural frequency analysis. The trend of the frequency decreases with the increase of the scour depth, but the curve is not smooth because of non-uniform foundation sections and layered soils. However, when the scour depth is such that pile cap is exposed, the changes in natural frequency with the scour depth are more obvious, and this is useful for measurement of the depth using bridge natural frequencies.     

On the impossibility of mean-field dynamos with some spherical symmetry of the motions

Abstract

A spherical mean-field dynamo model is considered in which both the mean motion and the mean electro-motive force due to fluctuating motions show some spherical symmetry. It is shown that under some reasonable assumptions the magnetic field is bound to decay to zero.     

A review of: “Problems of Solar and Stellar Oscillations. Iau Colloquium No. 66. Edited by D. O. Gough. D. Reidel Publishing Company, 493 pp. Df1195.00/U.S. $84.50 (ISBN 90-277-1554-8). 1983.”

Abstract

The model equations describing two-dimensional thermohaline convection of a Boussinesq fluid in a rotating horizontal layer are known to support multiple instabilities, depending on the values of certain control parameters (Arneodo et al., 1985). Most of these multiple instabilities have already been studied for double or triple diffusive convection, where behaviours ranging from simple steady to irregular motions have been found. Here we consider the one remaining bifurcation mentioned by Arneodo et al. (1985): the interaction between a steady and an oscillatory convection roll when the linear spectrum for a single wavenumber comprises one zero and one pair of purely imaginary eigenvalues. The method of centre manifolds and normal forms is used to derive evolution equations for the amplitudes of the convection rolls close to bifurcation and the behaviours associated with the equations is discussed.     

Finite amplitude convection and magnetic field generation in a rotating spherical shell

Abstract

Finite amplitude solutions for convection in a rotating spherical fluid shell with a radius ratio of η=0.4 are obtained numerically by the Galerkin method. The case of the azimuthal wavenumber m=2 is emphasized, but solutions with m=4 are also considered. The pronounced distinction between different modes at low Prandtl numbers found in a preceding linear analysis (Zhang and Busse, 1987) is also found with respect to nonlinear properties. Only the positive-ω-mode exhibits subcritical finite amplitude convection. The stability of the stationary drifting solutions with respect to hydrodynamic disturbances is analyzed and regions of stability are presented. A major part of the paper is concerned with the growth of magnetic disturbances. The critical magnetic Prandtl number for the onset of dynamo action has been determined as function of the Rayleigh and Taylor numbers for the Prandtl numbers P=0.1 and P=1.0. Stationary and oscillatory dynamos with both, dipolar and quadrupolar, symmetries are close competitors in the parameter space of the problem.     

On topographic core-mantle coupling

Abstract

The purpose of this note is two-fold: to draw attention to a perplexing difficulty connected with topographic core-mantle coupling, and to suggest tentatively an explanation. The difficulty is an apparent conflict between the most comprehensive theory of the coupling so far attempted (Anufriev and Braginsky, 1975a, b, 1977a, b) and recent explicit calculations based on magnetic and seismic information (Speith et al., 1986). It is argued that asymmetric deviations from Anufriev and Braginsky's basically axisymmetric model of the underlying core flow are capable of resolving the difficulty.     

近断层地震动作用下LRB隔震储罐地震反应特征研究

采用Bouc-Wen模型模拟铅芯橡胶支座(LRB)非线性力学性能,建立LRB基础隔震储罐地震反应的数值模型.利用叠加原理得到人工合成近断层脉冲型地震动,从实际典型近断层地震动和人工近断层脉冲地震动输入两个角度出发,以远场地震反应为参照,系统探讨了近断层地震动作用下LRB隔震储罐地震反应特点以及近断层地震动运动特征的影响规律.研究发现,近断层地震动作用下LRB隔震储罐地震反应明显大于远场地震反应值,显著的近断层脉冲效应是隔震储罐设计不容忽视的问题.近断层地震动的PGV/PGA值是决定隔震储罐地震反应的综合指标,PGV/PGA值较大时,隔震储罐地震反应脉冲效应突出.另外,脉冲周期、脉冲参与系数以及脉冲类型也是影响LRB隔震储罐地震反应的重要因素.在脉冲参与系数越大、含有脉冲数量越多的近断层地震动作用下,隔震储罐地震反应越强烈.当脉冲周期接近储罐晃动周期时,晃动波高会出现峰值.