Modelling Two-Phase Incompressible Flow in Porous Media Using Mixed Hybrid and Discontinuous Finite Elements

In this paper, we present a numerical model for simulating two-phase (oil–water and air–water) incompressible and immiscible flow in porous media. The mathematical model which is based on a fractional flow formulation is formed of two nonlinear partial differential equations: a mean pressure equation and a water saturation equation. These two equations can be solved in a sequential manner. Two numerical methods are used to discretize the equations of the two-phase flow model: mixed hybrid finite elements are used to treat the pressure equation, h-based Richards' equation and the diffusion term in the saturation equation, the advection term in the saturation equation is treated with the discontinuous finite elements. We propose a better way to calculate the nonlinear coefficients contained in our equations on each element of the discretized domain. In heterogeneous porous media, the saturation becomes discontinuous at the interface between two porous media. We show in this paper how to use the capillary pressure–saturation relationship in order to handle the saturation jump in the mixed hybrid finite element method. The two-phase flow simulator is verified against analytical solutions for some flow problems treated by other authors.     

323.15K下LiCl-Li2SO4-H2O体系热力学性质的等压研究

在50℃下用等压法研究了LiCl和Li2SO4的纯盐水溶液(浓度范围分别为0 8～8 1mol·kg-1,0 7～3 1mol·kg-1)以及LiCl-Li2SO4-H2O混合体系(离子强度范围为1 0～8 6mol·kg-1)的渗透系数和水活度等热力学性质,并与0℃,25℃的数据做了对比,揭示了渗透系数随离子强度、温度及Li2SO4离子浓度分数之间的变化规律。Zdanovskii规则非理想溶液表达式和扩展式可以用来描述该体系等压平衡浓度之间的关系。根据Pitzer离子相互作用模型对实验数据进行了理论分析,拟合求取了50℃下该体系的纯盐参数和混合参数,计算值和实验值相吻合,表明Pitzer模型能够较好的描述50℃下该体系的热力学性质。     

不连续岩体-拱坝系统的动力过程数值模拟方法

提出了一种分析三维不连续变形块体系统动力学过程和建坝的数值方法,该方法不但能够模拟位移不连续面在地震载荷作用下的张开和大尺度的摩擦滑移,而且还能给出不连续面上安全系数的分布和随时间的变化.通过有2个块体组成的系统验证了所提出的方法的正确性,同时用它分析了一个具有断层和伸缩缝的拱坝系统模型在地震载荷下的动力反应.     

A discontinuous Galerkin method for seismic soil–structure interaction analysis in the time domain

A technique for modeling transient wave propagation in unbounded media is extended and applied to seismic soil–structure interaction analysis in the time domain. The technique, based on the discontinuous Galerkin method, requires lower computational cost and less storage than the boundary element method, and the time‐stepping scheme resulting from Newmark's method in conjunction with the technique is unconditionally stable, allowing for efficient and robust time‐domain computations. To extend the technique to cases characterized by seismic excitation, the free‐field motion is used to compute effective forces, which are introduced on the boundary of the computational domain containing the structure and the soil in the vicinity of the structure. A numerical example on a dam–foundation system subjected to seismic excitation demonstrates the performance of the method. Copyright © 2003 John Wiley & Sons, Ltd.     

川滇地区现今块体水平运动变形特征分析

依据川滇地区1991，1999和2001年GPS测量数据，利用DDA正分析基本原理，结合现有的地质资料，对该地区近年来的水平构造运动特征进行了模拟分析，认为川滇地区水平运动整体速率在1999～2001年期间明显减弱；川滇决体主压应变以顺时针旋转为主，方向为SE—SEE；其南部是运动变形相对剧烈、易于积累能量的区域；未来一两年内，该地区地震活动水平不会太高。     

Evaluation of seismic site factors in the Mississippi Embayment. II. Probabilistic seismic hazard analysis with nonlinear site effects

An integrated probabilistic seismic hazard analysis procedure that incorporates nonlinear site effects, PSHA-NL, is developed and used to characterize the influence of thick deposits of the upper Mississippi Embayment (ME) on seismic site coefficients. PSHA-NL follows the methodology of the 2002 USGS hazard maps and generates a compatible set of ground motion records. The motions are propagated using nonlinear and equivalent linear site response analyses and ME properties developed in a companion paper and used to derive surface uniform hazard response spectra. A set of generic site coefficients are derived and summarized in a format similar to NEHRP site coefficients, with an added dimension of ME deposits thickness to the Paleozoic rock, a physically meaningful impedance boundary. These coefficients compare well with NEHRP site coefficients for 30 m profiles. For thicker soil profiles, developed site coefficients are lower at short periods and higher at long periods than NEHRP site coefficients.     

岩土力学数值分析方法研究

随着岩土工程领域的不断扩展与延伸,岩土力学数值分析方法得到了迅速发展,出现了各种各样的数值分析方法。将岩土力学常用的数值分析方法分为连续变形分析方法与非连续变形分析方法两大类,对各种方法的研究现状和最新进展进行了评述,着重介绍了各种方法的优缺点,并提出了解决问题的思路、方法和建议。     

基于奇异值分解和小波阈值方法的地震资料处理方法

奇异值分解是一种基于代数特征值的提取方法,小波变换是一种时间频率域的去噪方法,两者在去噪方面各有特点。将奇异值分解和小波阈值去噪的方法有机地结合起来,用于消除地震勘探资料中的随机噪声。仿真实验显示对于较低信噪比资料仍有很好的处理效果。     

我国太阳总辐射月总量资料的均一性检验及订正

采用建站在30a以上的月总辐射资料，以相关性最好的月份为代表月份，以经过均一性检验的日照资料为参考序列，用Potter法进行总辐射资料的均一性检验。结果表明，存在1个或2个间断点的测站占总数的58．6％；无间断点的测站占总数的34．3％。出现间断点的年份以20世纪70年代最多，其次是90年代。90年代初全国辐射仪器的换型对辐射资料没有普遍明显的影响。最后用回归法对存在间断点的绝大多数测站进行了均一性订正，订正后辐射和目照的年平均相关系拊兽谝增大平均增大近5％．     

Identification of the diffusion coefficient for a semi-linear hyperbolic-parabolic problem with discontinuous nonlinear terms

Numerical computations concerning geochemical processes in porous media require a knowledge of the values of the macroscopic parameters. In the context of the interpretation of leaching tests one such parameter is the apparent macroscopic diffusivity. In this paper, a mathematical model for dissolution in the presence of diffusion is considered. It is shown that the constant diffusion coefficient can be identifiable when the time integral of the flux is known. Moreover, a practical algorithm for computing the diffusion coefficient is proposed. This revised version was published online in August 2006 with corrections to the Cover Date.