Some recent advances in the modelling of soft rock joints in direct shear

This paper presents a review of recent developments made by the authors into the modelling of rock joints in direct shear. Careful observation of laboratory direct shear testing on concrete/rock joints containing two-dimensional roughness has allowed theoretical models of behaviour to be developed. The processes modelled include asperity sliding, asperity shearing, post-peak behaviour, asperity deformation and distribution of stresses on the joint interface. Model predictions compare extremely well with laboratory test results. These models were then applied to direct shear tests on rock/rock joints, and although behaviour in general was well predicted, the strength of rock/rock joints was over-predicted. Direct shear tests have also been carried out on samples containing both two- and three-dimensional roughness to test the accuracy of the two-dimensional approximation to roughness adopted in the theoretical models.     

Theoretical Prediction of Gibbs Free Energies of Formation for Crystalline α-MOOH and α-M2O3 Based on a Linear Free-Energy Relationship

In the present study, the modified Sverjensky–Molling equation, derived from a linear-free energy relationship, is used to predict the Gibbs free energies of formation of crystalline phases of α-MOOH (with a goethite structure) and α-M2O3 (with a hematite structure) from the known thermodynamic properties of the corresponding aqueous trivalent cations (M3+). The modified equation is expressed as ΔG0f,MVX=aMVXΔG0n,M3++bMVX+βMVXγM3+, where the coefficients aMVX, bMVX, and βMVX characterize a particular structural family of MvX (M is a trivalent cation [M3+] and X represents the remainder of the composition of solid); γ3+ is the ionic radius of trivalent cations (M3+); ΔG0f,MVX is the standard Gibbs free energy of formation of MvX; and ΔG0n,M3+ is the non-solvation energy of trivalent cations (M3+). By fitting the equation to the known experimental thermodynamic data, the coefficients for the goethite family (α-MOOH) are aMVX=0.8838, bMVX=?424.4431 (kcal/mol), and βMVX=115 (kcal/mol.?), while the coefficients for the hematite family (α-M2O3) are aMVX=1.7468, bMVX=?814.9573 (kcal/mol), and βMVX=278 (kcal/mol.?). The constrained relationship can be used to predict the standard Gibbs free energies of formation of crystalline phases and fictive phases (i.e. phases that are thermodynamically unstable and do not occur at standard conditions) within the isostructural families of goethite (α-MOOH) and hematite (α-M2O3) if the standard Gibbs free energies of formation of the trivalent cations are known.     

现代地球化学的理论思维结构

地球科学和地球化学的进展,不但加深和扩大了对地球的认识和理论成果,同时也总结提炼出了全新的科学思维方法。复杂性理论和系统思维是科学认识论的一次飞跃。现代地球化学建立了较完整的理论和方法体系,给地质科学认知地球提供了依据和支撑。作者论证了系统思维和逻辑演绎法相结合的思维方式,是研究和认识地球复杂巨系统以及推动地球科学理论创新的指导思想,并提出了运用现代地球化学理论构架和方法体系,综合地球科学系统思维与逻辑思维相结合方向的地球化学认知模型,为深入研究地球物质能量系统的行为和解决人类面临的资源、环境和灾害问题提供新的思路。     

氢氧化物族矿物的氧同位素分馏

应用增量方法计算了氢氧化物族矿物的氧同位素分馏,得到常见氢氧化物的18O富集顺序为：褐铁矿＞三水铝石＞针铁矿＞水镁石＞硬水铝石。氢氧化物与其对应的氧化物相比显著地富集18O。三价阳离子的氢氧化物和氧化物的18O富集顺序为：M（OH）3＞MO（OH）＞M2O3。Al（OH）3同质多象变体之间也存在一定的分馏。对于石英-氢氧化物、方解石-氢氧化物和氢氧化物-水体系,本文计算提供了在0-1200℃温度范围内三组内部一致的分馏系数方程。这些理论校准与合成实验结果和／或地表温度下的天然样品相吻合,特别针铁矿、勃姆石和硬水铝石与水之间的氧同位素分馏关系能够满足地质测温的要求。因此,对氢氧化物-水体系的氧同位素分析可望提供表生环境下可靠的地质温度计。     

The Mw = 6.0, 7 September 1999 Athens Earthquake

On 7 September 1999 at 11:56 GMT a destructive earthquake (M_w = 6.0) occurred close to Athens (Greece). The rupture process is examined using data from the Cornet local permanent network, as well as teleseismic recordings. Data recorded by a temporary seismological network were analyzed to study the aftershock sequence. The mainshock was relocated at 38.105°N, 23.565°E, about 20 km northwest of Athens. Four foreshocks were also relocated close to the mainshock. The modeling of teleseismic P and SH waves provides a well-constrained focal mechanism of the mainshock (strike = 105°, dip = 55° and rake = -80°) at a depth of 8 km and a seismic moment M₀ = 1.010²⁵ dyn·cm. The obtained fault plane solution represents normal faulting indicating an almost north-south extension. More than 3500 aftershocks were located, 1813 of which present RMS < 0.1 s and ERH, ERZ < 1.0 km. Two main clusters were distinguished, while the depth distribution is concentrated between 2 and 11 km. Over 1000 fault plane solutions of aftershocks were constrained, the majority of which also correspond to N–S extension. No surface breaks were observed but the fault plane solution of the mainshock is in agreement with the tectonics of the area and with the focal mechanisms obtained by aftershocks. The hypocenter of the mainshock is located on the deep western edge of the fault plane. The relocated epicenter coincides with the fringe that represents the highest deformation observed on the differential interferometric image. The calculated source duration is 5 sec, while the estimated dimensions of the fault are 15 km length and 10 km width. The source process is characterized by unilateral eastward rupture propagation, towards the city of Athens. An evident stop phase observed in the recordings of the Cornet local stations is interpreted as a barrier caused by the Aegaleo Mountain.     

浅层气对下伏地层深度预测影响定量研究及在渤海B油田中的应用

浅层气的存在造成下伏地层地震反射出现假象,使得深度预测不准确,增大了油田开发钻井风险。为提高浅层气下伏地层深度预测精度,利用理论模型分析量化描述浅层气对下伏地层的影响范围和影响程度,同时通过正演模拟手段构建时差定量校正量版,指导储层高精度深度预测。根据研究成果和认识,成功解释了渤海B油田两口已钻井深度误差出现的原因,并进一步利用时差校正量版提高了该油田区目标储层深度预测的精度,成功指导了后续多口开发井的准确钻探。     

不同含水条件下板岩力学实验研究与理论分析

以南水北调西线的板岩为研究对象,进行了不同泡水条件下的三轴压缩实验,通过引入相对吸水率的定义,分析了吸水率变化对岩石强度、弹性模量及泊松比等力学参数的影响,并从微观的角度分析了试样发生变化的机理,得出如下结论：(1)板岩在泡水后三轴抗压强度发生降低,最大降低辐度达46 %左右,抗压强度随相对吸水率的变化可以用函数式来有效的描述;(2)在泡水后,板岩的弹性模量与泊松比随着吸水率的增大而增大,但由于板岩自身各向异性明显,规律性不是特别明显,只能近似用关系式来模拟;(3)理论分析证实：实际泡水条件下的岩样,在受力挤压后,会经历排水与不排水两种过程,随着吸水率的提高,弹性模量与泊松比增大,验证了实验结果;(4)利用微观力学模型初步证实：泡水后板岩颗粒的吸水性与颗粒间毛细管力的减小是造成板岩泡水后发生软化的原因。     

煤层水压致裂后煤岩应力解析

水压致裂后煤岩应力分布规律对水压致裂防冲效果起关键性作用。采用理论研究方法得出高压注水压致裂后及卸水后水区和气区的孔隙、瓦斯压力和煤体应力解析解。结果表明,致裂后水区孔隙压力沿径向变化不大,与注水压力接近;气区瓦斯压力沿径向呈递减趋势;在水区外围一定范围内形成瓦斯压力升高区;水区煤体环向应力将会减小,直到变为拉应力;气区煤体径向应力沿径向递减。卸水后水区孔隙压力、煤体径向应力沿径向呈递增趋势;气区煤体径向应力沿径向呈递增趋势,趋近于原始煤体应力;气区煤体环向应力沿径向呈递减趋势;气区孔隙压力沿径向呈递减趋势。这为煤层水压致裂预防冲击地压提供理论基础。     

海洋沉积物动力学的示踪物方法

本文的论题是示踪物方法在海洋物源追踪和沉积物输运率计算中的应用 ,以及示踪物方法的普适性理论框架。定量的物源追踪需要适当的示踪标记和物质混合模型 ,有必要根据沉积动力过程的研究来确定示踪标记的变换函数 ,从而将改进的示踪标记用于混合模型分析。人工示踪物实验的现有方法以空间积分法最为常用 ,其关键是示踪物质心运动和沉积物活动层的界定。天然示踪物的质心位置是无法定义的 ,因此不能借用传统人工示踪物方法来估算物质输运率 ;在某些特殊情形下 ,可以利用天然示踪物的质量守衡原理来获得物质输运信息。今后 ,以示踪物质的连续方程为基础 ,有可能建立一种同时适用于人工和天然示踪物的普适理论框架。其中需解决的问题包括示踪物与现场物质的差异、沉积速率和活动层厚度的时间尺度、沉积物扩散过程等。     

Directions and possibilities of mathematical geology

This paper considers the present state of mathematical geology. Three directions are recognized: applied, theoretical, and mathematical. Applied mathematical geology includes formal use of mathematics to solve problems and computer processing of data. Success is achieved by a correspondence of mathematical methods used to the nature of geological data. This correspondence can be demonstrated by purely mathematical means. Theoretical mathematical geology uses mathematics as a language of geology; however, a number of methodological problems must be solved: formalization of initial geological concepts and creation of a strict conceptual basis, substantiation of initial principles of mathematical simulation, creation of theoretical geological models, problems of elementary and coincidence in geology, and methodological substantiations of possibilities of any mathematical model to approximate geological models. The essense and significance of these problems are considered. The main task of mathematical geology is to prove its correspondence to the nature of the geological objects studied, geological data obtained, and geological problems solvable. Finally, the main problems of mathematical geology are not so much mathematical as geological and methodological.