深源远震波形反演介质层状结构

在远场正演计算的基础上，运用３０°～９０°的四个不同深源远震记录波形，通过试错法选取台站下方的结构，借助波场反向重建算法，反演确定宜昌台站下方的速度结构。这一算法的实现有助于提高用试错法所得结果的可靠性和精确度，对于利用天然地震资料研究地球深部构造有一定的实用价值。     

Full-wave synthetic acoustic logs in porous media

In order to investigate the permeability of porous media,it is necessary to figure out the property of wave propagation in this medium. The authors calculated the full waveforms and the 2-dimensional spectrum of fluid saturated porous media,analyzed the wave propagation characteristics of different permeable stratums. From the theoretical acoustic pressure waveform received at the well-axis and the 2-D spectrum it is revealed that longitudinal-,transverse-,Stoneley-,pseudo-Rayleigh-,longitudinal mode-and transverse mode waves are different in characters in various permeable strata. Stoneley wave is sensitive to the change of permeability,and its frequency range and amplitude attenuation are influenced by permeability.     

煤系变形介质条件研究在煤层小断层预测中的应用——以下峪口矿2号煤层为例

煤系介质条件作为构造形成和发育的的物质基础，必然影响断层的发育。以陕西下峪口煤矿2号煤层为对象，研究了煤层围岩介质参数与煤层小断层间的关系，确定了控制2号煤层小断层发育的4个主要因素，并据此建立了预测小断层发育程度的数学模型。     

爆炸应力波作用下缺陷介质裂纹扩展的动态分析

采用透射式焦散线测试系统,进行爆炸应力波作用下缺陷介质裂纹扩展试验,研究了含与炮孔共线的预制裂隙介质裂纹扩展速度、加速度、裂纹尖端动态应力强度因子和动态能量释放率的变化规律以及它们之间的变化关系。试验结果表明,在爆炸应力波作用下裂隙两端产生了两条翼裂纹A、B,翼裂纹A的长度较翼裂纹B长,两条翼裂纹向相反的方向扩展;在翼裂纹扩展过程中,存在着加速与减速的过程,扩展速度瞬间达到峰值,其后逐渐振荡下降;动态应力强度因子也呈现瞬间达最大值到逐渐减小连续振荡变化的趋势,动态应力强度因子 > ;翼裂纹尖端的动态能量释放率对裂纹扩展具有驱动作用。     

Convective stability analysis of the long-term storage of carbon dioxide in deep saline aquifers

Deep saline aquifers are one of the most suitable geologic formations for carbon sequestration. The linear and global stability analysis of the time-dependent density-driven convection in deep saline aquifers is presented for long-term storage of carbon dioxide (CO₂). The convective mixing that can greatly accelerate the CO₂ dissolution into saline aquifers arises because the density of brine increases upon the dissolution of CO₂ and such a density difference may induce instability. The effects of anisotropic permeability on the stability criteria, such as the critical time for the appearance of convective phenomena and the critical wavelength of the most unstable perturbation, are investigated with linear and global stability analysis. The linear stability analysis provides a sufficient condition for instability while the global stability analysis yields a sufficient condition for stability. The results obtained from these two approaches are not exactly the same but show a consistent trend, both indicating that the anisotropic system becomes more unstable when either the vertical or horizontal permeability increases.     

二维饱和多孔介质因点汇诱发比奥固结的解析解

给出了有限二维饱和多孔介质因点汇诱发的Biot固结的一个解析解。其中假设多孔介质为均匀各向同性和线弹性,假设孔隙压力场符合第1类边界条件,数学模型采用可压缩多孔介质模型。利用傅里叶和拉普拉斯变换及相应反演获得了双重无穷项级数和形式的精确解。然后特别探讨了定流量点汇诱发的稳态解析解,并用文献现有解析解进行了验证。所提出的解析解适合于验证数值解,并可用于深入分析有限二维多孔介质的流-固耦合行为。     

基于SPH方法的非饱和多孔介质相变耦合研究

在考虑相变的热能平衡方程和非饱和水分迁移质量控制方程的基础上,建立温度场-水分场的耦合模型,并采用一种无网格粒子算法（SPH）进行数值求解。其中,耦合方程中考虑了水流传热以及温度势对水流的直接驱动,在不考虑相变的情况下,该耦合模型可退化为常温下的水-热耦合模型,故可用于模拟冻融循环的相关问题。从求解热能平衡方程中的含冰量出发,实现解耦并对半无限单向冻结条件下介质内非稳态温度场和体积含水率分布场进行模拟,将耦合作用下的温度场与不耦合的解析解进行对比,反映出水分迁移对温度场存在较大影响。最后,求解了路基边坡在季节性周期温度边界下,温度场、水分场分布的演变规律,并评估了边坡阴阳面受热不均对水热两场分布的影响。计算结果基本能反映土冻结相变的实际物理过程,光滑粒子算法可以用于尝试解决冻土领域的其他相关问题。     

Response of a tunnel deeply embedded in a viscoelastic medium

This paper presents a three‐dimensional energy‐based solution for the time‐dependent response of a deeply embedded and unsupported semi‐infinite tunnel of circular cross‐section. The tunnel is taken to be excavated quasi‐instantaneously from an infinite rock body that initially exhibits an isotropic stress state and that is made up of a homogeneous, isotropic and viscoelastic material. The viscoelastic behaviour is modelled by means of Burger's model, and the rock is taken to behave volumetrically linear elastic and to exhibit exclusively deviatoric creep. This viscoelastic problem is transformed into the Laplace domain, where it represents a quasi‐elastic problem. The displacement fields in the new solution are taken to be the products of independent functions that vary in the radial and longitudinal directions. The differential equations governing the displacements of the system and appropriate boundary conditions are obtained using the principle of minimum potential energy. The solutions for these governing equations in the Laplace domain are then obtained analytically and numerically using a one‐dimensional finite difference technique. The results are then transformed back into the time domain using an efficient numerical scheme. The accuracy of the new solution is comparable with that of a finite element analysis but requires much less computation effort. Copyright © 2011 John Wiley & Sons, Ltd.