岛状多年冻土管道周围热-水-应力耦合数值模拟

考虑大气温度以及管内流动介质温度的影响,建立由多个土层构成的岛状冻土物理模型,研究塔河地区埋地输油管道周围热-水-应力（T-H-M）耦合作用,分析水-热-应力耦合作用下管道的受力情况并进行管道强度验算。采用多物理场耦合软件COMSOL Multiphysics进行数值模拟,计算在T-H-M三场耦合及非耦合两种工况下岛状多年冻土管道周围应力和沉降曲线等。分析比较表明,岛状多年冻土区的温度,冰水相变以及水分迁移对应力场存在较大的影响;水分场、温度场以及应力场三者之间耦合作用明显;岛状冻土随季节呈周期性冻融,管道周围出现较大应力集中和沉降现象,会对管道安全产生一定影响。     

岩体结构与岩体水力耦合计算模型

水力耦合作用是岩体中力学过程与渗流过程相互作用的物理过程。水力耦合机理的理解是水力耦合分析的关键问题,其耦合机理是由岩体结构特性决定的。在分析岩体水力耦合过程基础上,根据岩体的基本结构及代表性单元体（REV）是否存在提出了建立水力耦合模型的方法。当裂隙岩体中不存在代表性单元体（REV）时,提出了裂隙岩体多重介质流固耦合分析的全耦合数学模型,给出流固耦合模型数值方法求解的数学模型及有限元计算表达式。程序编制和验证工作正在进行中。     

Interactive effects of the type of milling media and CuSO4 addition on the flotation performance of sulphide minerals from Merensky ore Part II: Froth stability

The limitations of pulp chemistry measurements in the flotation of a platinum group mineral (PGM) bearing Merensky ore were demonstrated in Part 1 of this article. In this paper the importance of the contribution of the froth structure due to changing froth stability is analysed using the batch flotation data. The effects of mild steel (MS) and stainless steel (SS) milling media and the addition of copper sulphate on the flotation performance of the sulphide minerals in Merensky ore have been evaluated in relation to the changes in stability of the froth phase. The effects of pulp chemistry and froth stability on the flotation of sulphide minerals were distinguished by using two different rate constants (k_t and k_w). The rate constant (k_w) calculated as a function of cumulative water recovery was used to describe characteristics of froth phase and k_t was calculated as a function of flotation time. The results revealed that the type of grinding media and copper sulphate addition had an interactive effect on the froth stability. While mild steel (MS) milling increased the froth stability due to the presence of hydrophilic iron hydroxides and colloidal metallic iron, the addition of copper sulphate reduced the stability, especially with stainless steel (SS) milling. Copper sulphate addition had a dual role in the flotation of Merensky ore in that it caused destabilisation of the froth zone as well as activation of selected sulphide minerals. The dominant effect was found to depend on the type of milling media and floatability of the mineral in question and this work has demonstrated the importance of using a combination of measurements to evaluate flotation performance holistically.     

多孔黏弹介质平面波理论及波场数值模拟

本文将黏弹固体骨架与各向异性多孔弹性模型相结合,并将地层黏弹性的特性考虑其中,更好地刻画了快P波和横波在地震频带的衰减现象。将多孔黏弹模型模拟结果同传统弹性模型模拟结果对比显示,多孔黏弹模型在模拟真实储层方面更具优势。     

线性粘弹介质中地震波场数值模拟

本文将以往两种粘弹介质中地震波模拟方法的优点结合起来，以模型理论和积分本构方程为基础，从理论上分析了模型对地震波场的影响；采用交错网格有限差分法对粘弹介质中的地震波进行数值模拟.数值计算结果表明该方法不仅便于计算，同时也便于从力学的角度来分析地震波的传播.数值计算结果与理论分析一致，说明这种方法可以更为有效地模拟粘弹介质中地震波的传播.     

垂向非均匀介质中转换波转换点计算

转换点位置的计算是转换波资料处理中的一个关键问题. 本文提出了分别基于速度随深度线性变化、速度随垂直走时线性变化、慢度随深度线性变化和慢度随垂直走时线性变化四种等效垂向非均匀介质情况下转换点位置的计算方法. 研究了通过速度拟合、走时近似和相似系数谱三种方式选择合适的等效速度方法. 结合理论模型对非均匀介质转换点计算方法、渐进转换点计算方法、Thomsen近似公式和均匀介质解析计算方法的误差进行了分析,结果表明非均匀介质转换点计算方法能更准确地计算转换点位置.     

周期性层状含孔隙、裂隙介质模型纵波衰减特征

地震波在含孔隙、裂隙斑块饱和介质传播过程中会诱发多个尺度孔隙流体流动而产生衰减和速度频散.在含有宏观尺度“Biot流”和介观尺度“局域流”衰减诱导机制的周期性层状孔隙介质模型基础上,引入了微观尺度硬币型和尖灭型裂隙“喷射流”的影响,构建了周期性层状含孔隙、裂隙介质模型.利用双解耦弹性波动方程的方法数值计算了该模型地震频带的纵波衰减和速度频散并与周期性层状孔隙介质模型做了对比研究.分析了该模型在不同裂隙参数（裂隙密度、裂隙纵横比）及裂隙体积含量下的纵波衰减和频散特征,裂隙密度越高对于纵波衰减和频散的影响越大,裂隙纵横比越小,由裂隙引起的纵波衰减部分向高频段移动,裂隙体积含量越少,纵波衰减先降低后小幅增加再降低,频散速度增加,并逐渐接近于周期性层状孔隙介质模型的纵波衰减和频散速度曲线.最后研究了周期性层状含孔隙、裂隙介质模型有效平面波模量的高低频极限以及流固相对位移在该模型中的分布特征.     

Another look at the conceptual fundamentals of porous media upscaling

We suggest a critical look at the epistemic foundations of the porous media upscaling problem that focuses on conceptual processes at work and not merely on form manipulations. We explore the way in which critical aspects of scientific methodology make their appearance in the upscaling context, thus generating useful effective parameters in practice. The fons et origo of our approach is a conceptual blending of knowledge states that requires the revision of the traditional method of scientific argument underlying most upscaling techniques. By contrast to previous techniques, the scientific reasoning of the proposed upscaling approach is based on a stochastic model that involves teleologic solutions and stochastic logic integration principles. The syllogistic form of the approach has important advantages over the traditional reasoning scheme of porous media upscaling, such as: it allows the rigorous derivation of the joint probability distributions of hydraulic gradients and conductivities across space; it imposes no restriction on the functional form of the effective parameters or the shape of the probability laws governing the random media (non-Gaussian distributions, multiple-point statistics and non-linear models are automatically incorporated); it relies on sound methodological principles rather than being ad hoc; and it offers the rational means for integrating the multifarious core knowledge bases and uncertain site-specific information sources about the subsurface system. Previous upscaling results are derived as special cases of the proposed upscaling approach under limited conditions of porous media flow, a fact that further demonstrates the generalization power of the approach. Our hope is that looking at the upscaling problem in this novel way will direct further attention to the methodological exploration of the problem at the length and the detail that it deserves.I would like to thank Drs. A. Kolovos and D.T. Hristopulos for their valuable comments. The work was supported by grants from the Army Research Office (Grant no. DAAG55–98–1-0289) and the National Institute of Environmental Health Sciences (P42-ES05948 & P30-ES10126).     

含泥质低孔渗各向异性黏弹性介质中的波频散和衰减研究

本文定义了各向异性黏弹性参数修正因子,并将其引入到黏弹性模型中以体现泥质含量对黏弹性机制的影响,同时将波传播过程中孔隙介质骨架黏弹性力学机制与两种孔隙流体流动力学机制（Biot流动和喷射流动机制）有机地统一起来处理,从而给出了描述含泥质低孔渗孔隙各向异性介质中波传播规律的黏弹性Biot/squirt (BISQ)模型.数值计算结果表明,入射波的方位角、各向异性渗透率以及泥质含量等对含流体复杂孔隙介质中波频散和衰减的影响具有显著的方位各向异性特征,在低频范围内（地震波勘探频率）黏弹性力学机制对波传播能量的衰减起主导作用.     

含气介质临界点、流体和骨架弹性参数的数值计算方法

目前临界孔隙度及孔隙介质的研究主要为测试分析的实验方法,如何运用数值计算方法求取孔隙流体介质临界点和组分弹性参数一直是备受关注的课题.本文提出了求取孔隙介质临界点、孔隙内流体和骨架弹性参数的数值计算公式及方法,并结合含气样品测试数据实现了这种数值计算.文章首先给出了3个包含孔隙度的线性方程,通过有机组合每个线性方程中的两个系数得到求取有关弹性参数的数值计算公式.然后,详细阐述了计算的方法和步骤以及需要注意的问题,把室内含气砂岩样品实测的综合介质密度、纵波速度和横波速度作为数值计算的输入数据,求得了临界点和流体及骨架等相关弹性参数的具体数值.通过比较分析本文数学方法求得的计算数据与实验方法测得的实测数据,表明了该数值计算公式的正确性和数值实现方法的有效性.