非混溶饱和两相渗流与孔隙介质耦合作用的理论研究——数学模型

基于连续介质力学的一般理论,在黑油多相渗流理论模型的基础上,推导了非混溶饱和两相渗流与变形孔隙介质耦合作用的数学模型方程,并对流固耦合效应进行了初步的讨论     

Seismic Wave Propagation and Excitation in Multi-layered Media with Irregular Interfaces Part (Ⅰ): SH Case

Study of seismic wave excitation and propagation in laterally heterogeneous media was an active and important subject in seismology in the past two decades, numerous analytical and numerical efforts have been made in this research field. In this article, I have, first, made a brief review on those developments and then introduced and summarized a unified and efficient method, global generalized reflection-transmission (abbreviated to R/T thereafter) matrices method, for synthetic seismograms in multi-layered media with irregular interfaces developed by the author [24~26]. As demonstrated in this article, this method could be regarded as an extension of the generalized R/T coefficients method for the horizontally layered case [2,5] to the layered media with irregularly shaped interfaces by incorporating the T matrices technique [27,28]. Because of the use of a recursive scheme in computing the global generalized R/T matrices, this method is efficient, particularly for the case with a large number of irre     

混合ADI-FDTD亚网格技术在探地雷达频散媒质中的高效正演

提出混合ADI-FDTD亚网格技术开展频散介质GPR正演,即在物性参数变化剧烈局部区域采用细网格剖分ADI-FDTD计算,其他的区域采用粗网格剖分常规FDTD计算,ADI-FDTD突破了CFL条件的限制,可选取与粗网格一致的大时间步长,有效地提高了计算效率.本文首先基于Debye方程,推导了粗网格FDTD及细网格ADI-FDTD频散介质差分格式,着重对粗细两种网格结合的场值交换方式进行了深入探讨,给出了该算法的计算流程.然后以一个薄层模型为例,分别应用粗网格、细网格、混合ADI-FDTD亚网格算法对该模型进行正演,计算资源的占用及模拟精度说明了混合ADI-FDTD亚网格算法的优势.最后,建立频散介质与非频散介质的组合模型,应用3种方法对该模型进行正演,对比3种方法优劣,分析雷达剖面中非频散介质及频散介质中波形特征,有效地指导雷达资料的精确解释.     

Refraction Effect in Geothermal Heat Flow Due to a 3-D Prism in a Two-Layered Earth

A mathematical model of the disturbance of the stationary geothermal field due to a three-dimensional perturbing body embedded near the surface of a two layered earth. The theoretical analysis is based on the generalized theory of the double-layer potential, similar to the boundary integral method likely to direct current geoelectricity problems. Special attention is paid to cases in which a pan of the body's surface is tangential to the planar boundaries of the substratum or of the surface of the earth The numerical calculations are performed for a three-dimensional prism using methods similar to the boundary element method (BEM). Numerous graphs are shown for the disturbance of the temperature field and of the heat flow on the surface of the earth and along vertical profiles (boreholes).     

Numerical modeling of two-phase flow in heterogeneous permeable media with different capillarity pressures

Contrast in capillary pressure of heterogeneous permeable media can have a significant effect on the flow path in two-phase immiscible flow. Very little work has appeared on the subject of capillary heterogeneity despite the fact that in certain cases it may be as important as permeability heterogeneity. The discontinuity in saturation as a result of capillary continuity, and in some cases capillary discontinuity may arise from contrast in capillary pressure functions in heterogeneous permeable media leading to complications in numerical modeling. There are also other challenges for accurate numerical modeling due to distorted unstructured grids because of the grid orientation and numerical dispersion effects. Limited attempts have been made in the literature to assess the accuracy of fluid flow modeling in heterogeneous permeable media with capillarity heterogeneity. The basic mixed finite element (MFE) framework is a superior method for accurate flux calculation in heterogeneous media in comparison to the conventional finite difference and finite volume approaches. However, a deficiency in the MFE from the direct use of fractional flow formulation has been recognized lately in application to flow in permeable media with capillary heterogeneity. In this work, we propose a new consistent formulation in 3D in which the total velocity is expressed in terms of the wetting-phase potential gradient and the capillary potential gradient. In our formulation, the coefficient of the wetting potential gradient is in terms of the total mobility which is smoother than the wetting mobility. We combine the MFE and discontinuous Galerkin (DG) methods to solve the pressure equation and the saturation equation, respectively. Our numerical model is verified with 1D analytical solutions in homogeneous and heterogeneous media. We also present 2D examples to demonstrate the significance of capillary heterogeneity in flow, and a 3D example to demonstrate the negligible effect of distorted meshes on the numerical solution in our proposed algorithm.     

Nonequilibrium effects in models of three-phase flow in porous media

In this paper we extend to three-phase flow the nonequilibrium formalism proposed by Barenblatt and co-workers for two-phase porous media flow. The underlying idea is to include nonequilibrium effects by introducing a pair of effective water and gas saturations, which are linked to the actual saturations by a local evolution equation. We illustrate and analyze how nonequilibrium effects lead to qualitative and quantitative differences in the solution of the three-phase flow equations.     

Application of implicit sub-time stepping to simulate flow and transport in fractured porous media

In general, the accuracy of numerical simulations is determined by spatial and temporal discretization levels. In fractured porous media, the time step size is a key factor in controlling the solution accuracy for a given spatial discretization. If the time step size is restricted by the relatively rapid responses in the fracture domain to maintain an acceptable level of accuracy in the entire simulation domain, the matrix tends to be temporally over-discretized. Implicit sub-time stepping applies smaller sub-time steps only to the sub-domain where the accuracy requirements are less tolerant and is most suitable for problems where the response is high in only a small portion of the domain, such as within and near the fractures in fractured porous media. It is demonstrated with illustrative examples that implicit sub-time stepping can significantly improve the simulation efficiency with minimal loss in accuracy when simulating flow and transport in fractured porous media. The methodology is successfully applied to density-dependent flow and transport simulations in a Canadian Shield environment, where the flow and transport is dominated by discrete, highly conductive fracture zones.     

水介质对水下目标体磁场的影响

根据扬子江航道上一系列磁测及验证结果,指出在不同介质中的相似磁性物体所引起的磁场强度会产生一定的变化,实例说明：这一现象可能对磁法勘查带来的影响.文章没有对不同介质中磁性体磁场的变化机制进行讨论.     

胆脂瘤型中耳炎CT表现及临床价值

目的: 研究胆脂瘤型中耳炎CT表现,并探讨其临床价值。 方法: 15例19耳胆脂瘤型中耳炎,均经乳突根治术及病理证实,结合手术对照分析其CT表现。 结果: ①发生部位:主要位于中耳腔,占93%(14/15)。位于中耳腔及乳突窦,占27%(4/15)。位于岩锥,占7%(1/15)。②CT表现:窦腔内软组织肿块,骨质破坏,prussaks间隙增宽,窦腔破坏区周边硬化。 结论: CT检查能准确显示胆脂瘤型中耳炎影像特征,对胆脂瘤型中耳炎诊断和鉴别诊断有重要作用,对手术具有十分重要的指导意义。     

考虑土体固液二相性质的地下管线地震反应研究

本文将管线周围的土体视为固-液二相介质，采用考虑地震过程中土体内孔隙水压力的增长和消散的有效应力分析方法，并以非线性的土体本构模型模拟土体动力特性的改变以及在孔隙水压作用下的变形，以求接近土体的真实地震反应。研究表明随着土体软化程度的不同，受到土体约束的管线应力亦发生相应变化。文中详细研究了地震动峰值、管线直径、行波波速和入射角度等对管线应力的影响。