完全匹配层吸收边界在孔隙介质弹性波模拟中的应用

模拟弹性波在孔隙介质中传播,需要稳定有效的吸收边界来消除或尽可能的减小由人工边界引起的虚假反射. 本文在前人工作基础上,首次建立了弹性孔隙介质情况下完全匹配层吸收边界的高阶速度-应力交错网格有限差分算法,并详细讨论了完全匹配层的构建及其有限差分算法实现. 首先,本文通过均匀孔隙模型的数值解与解析解的对比,验证所提出的数值方法的正确性;然后,本文考察了完全匹配层对不同入射角度入射波和自由表面上的瑞利波的吸收性能,将完全匹配层与廖氏和阻尼吸收边界进行了对比,研究了这三种吸收边界在不同吸收厚度情况下对弹性波吸收能力. 数值结果表明,在孔隙介质中,完全匹配层作为吸收边界能十分有效地吸收衰减外行波,无论对体波还是面波,是一种高效边界吸收算法.     

北太平洋对人为二氧化碳吸收的数值模拟

采用一个开边界海盆尺度环流模式研究人为CO2在北太平洋的吸收和分布,并与闭边界模式的模拟结果进行了比较.模拟结果表明,西北太平洋和赤道东太平洋是两个重要的人为CO2汇.使用较大的等密度面扩散系数使得西北太平洋人为CO2通量增大,而赤道海区通量减小（RUN2）.与闭边界模式相比,开边界模式中该两个区域的人为CO2通量都增加了.1800～1997年间,北太平洋共吸收人为CO2 2375GtC(1Gt=1×1015 g, RUN1).副极地海区是人为CO2的一个重要输出区,能输出人为CO2吸收量的38％～54％,而20°N～30°N海区是人为CO2的一个重要贮存区,占整个北太平洋的24％.开边界对于10°N以南太平洋吸收和贮存人为CO2有很大影响.与基于观测资料的估计相比,虽然模式低估了西北太平洋的人为CO2的穿透,高估了东北太平洋的人为CO2的穿透,但总的说来,模式比较好地估计了人为CO2在北太平洋的贮存.     

国内外岩溶塌陷监测方法综述

岩溶塌陷是全球广泛分布的地质灾害问题,严重影响和制约了岩溶地区社会、经济发展,因此,开展岩溶塌陷监测工作,对塌陷预测预报具有及其重要的作用。岩溶塌陷的监测技术与塌陷形成的机理关系密切。该文分析了岩溶塌陷的分布规律和形成的条件,系统总结了国内外岩溶塌陷的监测内容、技术方法、各自优缺点,如BOTDR监测、TDR监测、地质雷达监测、水动力条件监测及地-沉降等方法,为系统性地开展岩溶塌陷监测工作提供技术和方法参考,展望了未来的岩溶塌陷监测工作,认为岩溶塌陷监测需结合塌陷的形成机理、塌陷地区覆盖层特征及水动力条件,综合利用多种监测手段联合进行,并形成完善的监测网络。     

亚太区域不同遮挡情形下多系统组合PPP定位效能分析

利用多测站实测GNSS观测数据,模拟四周遮挡（城市环境）、单侧遮挡（峡谷环境）和顶空遮挡（大型桥梁路基环境）3种遮挡环境,从可用卫星数、PDOP值、可用历元率、定位精度和收敛时间5个方面,综合分析多系统PPP技术的定位服务效能。结果显示,相比GPS单系统PPP定位结果,在无遮挡环境下,四系统融合PPP技术在可用卫星数、PDOP值、可用历元率、定位精度和收敛时间方面分别改进300%、40%、2%、20%和50%;在遮挡环境下,分别改进300%、60%、25%、39%和52%。     

中条裂谷的递进演化与王屋山地区铜矿成矿地质条件分析

对中条—王屋山地区前寒武纪沉积建造、岩浆建造、构造形态的分析后,反演出中条裂谷的递进演化过程,认为王屋山地区在中条裂谷和区域地壳的演化过程中形成了极为有利的成矿条件:①形成多层含矿建造;②通达地幔的深大断裂提供了丰富的矿源和成矿空间;③频繁的岩浆活动补充了矿质并提供了成矿热能;④剥离断层系统进一步提供了扩容空间;⑤以活化转移为主要机制的强烈变形变质作用对已形成的含矿建造(矿源层)进行改造。因此,王屋山地区具有巨大的找矿潜力。     

漳州市山体滑坡主要影响因素分析

对漳州市区域内近十几年来所发生滑坡的地质灾害点进行调查研究,提出了漳州地区山体滑坡主要受地形、地质、气象和人类工程活动因素的综合影响。一般地形坡度在25°～50°,坡型为阶梯型、人工边坡无支挡,土质边坡坡度大于70°与坡高大于4m者,连续10d内降雨量在170mm以上,且日降雨量在100mm的强降雨时,发生滑坡的概率,达71.7%。故此,降雨量170mm以上、日降雨量大于或等于100mm时,可作为漳州市地质灾害预警预报的指标。     

A note on the consolidation settlement of a rigid circular foundation on a poroelastic halfspace

This paper uses Biot's poroelasticity approach to examine the consolidation behaviour of a rigid foundation with a frictionless base in contact with a poroelastic halfspace. The mathematical development of the mixed boundary value problem involves a set of dual integral equations in the Laplace transform domain which cannot be conveniently solved by employing conventional procedures. In this paper, a numerical solution is developed using a scheme where the contact normal stress is approximated by a discretized equivalent. The influence of limiting drainage boundary conditions at the entire surface of the halfspace on the degree of consolidation of the rigid circular foundation is investigated. The results obtained in this study are compared with the corresponding results given in the literature. Copyright © 2016 John Wiley & Sons, Ltd.     

One‐dimensional dynamics of saturated incompressible porous media: analytical solutions and influence of inertia terms

The complexity of formulations for the hydromechanical coupled mechanics of porous media is typically minimised by simplifying assumptions such as neglecting the effect of inertia terms. For example, three formulations commonly employed to model practical problems are classified as fully dynamic, simplified dynamic and quasi‐static. Thus, depending on the porous media conditions, each formulation will have advantages and limitations. This paper presents a comprehensive analysis of these limitations when solving one‐dimensional fully saturated porous media problems in addition to a new solution that considers a more general loading situation. A phase diagram is developed to assist on the selection of which formulation is more appropriate and convenient regarding particular cases of porosity and hydraulic conductivity values. Non‐dimensional formulations are proposed to achieve this goal. Results using the analytical solutions are compared against numerical values obtained with the finite element method, and the effect of porosity is investigated. Copyright © 2016 John Wiley & Sons, Ltd.     

Theoretical background for the inversion of seismic waveforms including elasticity and attenuation

To account for elastic and attenuating effects in the elastic wave equation, the stress-strain relationship can be defined through a general, anisotropic, causal relaxation function ^ijkl (x, ). Then, the wave equation operator is not necessarily symmetric (self-adjoint), but the reciprocity property is still satisfied. The representation theorem contains a term proportional to the history of strain. The dual problem consists of solving the wave equation withfinal time conditions and an anti-causal relaxation function. The problem of interpretation of seismic waveforms can be set as the nonlinear inverse problem of estimating the matter density (x) and all the functions ^ijkl (x, ). This inverse problem can be solved using iterative gradient methods, each iteration consisting of the propagation of the actual source in the current medium, with causal attenuation, the propagation of the residuals—acting as if they were sources—backwards in time, with anti-causal attenuation, and the correlation of the two wavefields thus obtained.     

A simplified computational approach for prediction of transient climatic conditions in underground mine airways

Summary A simplified approach for prediction of transient climatic conditions in an underground airway os presented. The dry surface temperature of the airway is calculated using a finite difference scheme, whereas the wet surface temperature is determined following the concept of pseudo-base-temperature of air and using the value of thermal gradient at the dry surface as obtained from the finite difference solution of the temperature profile in the dry rock. The air-current temperatures are regularly evaluated after a suitable time span so as to simulate a real-case situation, and the re-evaluated air-current temperatures are used for further calculations. For increasing age of the airway the wet surface temperature obtained by the finite difference scheme, with an appropriate value of thermal diffusivity, is compared with its value determined by the pseudo-base-temperature approach.