基于改进Bingham模型的软岩参数非定常三维非线性黏弹塑性蠕变本构研究

传统线性元件组合模型难以描述岩石非线性加速蠕变阶段特性,给工程应用带来了巨大困难。提出一种新的应变触发非线性黏壶元件,并将其与Bingham模型串联,建立了可以描述岩石等速和加速蠕变特性的一维本构模型,进而推导出相应的三维蠕变本构模型,再将三轴压缩过程中岩石弹性模量的衰减方程引入该三维蠕变本构方程,得到一个能反映岩石蠕变全过程特性的三维非线性黏弹塑蠕变本构模型,根据泥岩蠕变试验数据对该模型参数进行了辨识。模型理论计算值与泥岩蠕变试验数据的对比表明改进的Bingham模型不仅可以充分反映软岩的等速蠕变特性,还可以很好地描述软岩的衰减和加速蠕变规律,且模型元件少,组合形式简单,为软岩非线性蠕变本构模型研究提供了新的借鉴。     

陡倾层状斜向岩层视向溃屈机制力学分析——以陕西山阳滑坡为例

以陕西山阳滑坡为例,分析了陡倾层状斜向岩层岩质滑坡的视向滑动特征、滑移-溃屈破坏模式与机制。基于梁板理论、层状板裂结构岩体弯曲-溃屈破坏的力学模型在考虑自重、地下水静水压力作用及斜倾层状山体视向滑动侧向摩阻力作用的影响下,采用岩体结构力学分析的方法建立了相应的力学模型;经过力学分析,推导出基于斜坡自重、地下水静水压力与侧向摩阻力作用下的陡倾层状斜向岩层斜坡溃屈段长度条件方程。为验证条件方程的正确性,以山阳滑坡为例进行了验算,最后得出与实际调查较一致的结果,为防御陡倾层状斜向岩层斜坡产生视向溃屈破坏提供依据。      

非均匀介质地震AVO响应模拟及分析

地震AVO （amplitude versus offset）技术是一项利用振幅信息研究岩性、检测油气的地震勘探技术。常规方法基于Zoeppritz方程计算模型界面处的反射系数,而实际地下非均匀介质中地震反射特征不仅与入射角度、物性差异有关,还与入射波频率、地层厚度、薄互层结构等因素有关。为此,应用传播矩阵理论充分考虑与这些因素有关的调谐干涉等传播效应,针对实际数据计算非均匀地下介质的高精度合成地震记录,对比Zoeppritz方程、Shuey二项近似方程、Shuey三项近似方程以及传播矩阵算法的模拟效果。研究发现：在小角度入射时Zoeppritz方程、Shuey二项近似方程和Shuey三项近似方程的反射振幅和波形基本一致,大角度入射时Zoeppritz方程与Shuey三项近似方程接近;Zoeppritz算法的模拟结果在小角度入射和浅层情况下与传播矩阵算法差别较小,而在大角度入射和深层情况下与传播矩阵差别较大,说明层间多次波的调谐干涉等传播效应不可忽略。     

A modified method to estimate chain inverse catenary profile in clay based on chain equation and chain yield envelope

Offshore floating facilities are fixed by anchoring systems embedded in seabed soils through chains or ropes. The chain inverse catenary profile embedded in soils influences both the anchor failure mechanism and the anchor holding capacity. The chain mobilizes varying soil normal and tangential resistances during motion, hence it is with difficulty to depict the chain profile. The present work proposed a modified method to estimate the chain inverse catenary profile with high accuracy based on the chain equations and the chain yield envelope. A testing arrangement with three load cells and two MEMS (Micro-electromechanical systems) accelerometers included was designed in model tests. By model tests, the loading combinations of the soil tangential and normal resistances on the chain were obtained and the yield envelopes for both chain and rope were determined. In addition, supplemental model tests were performed to validate the modified method proposed in this study, and the testing results indicated that the estimated chain inverse catenary profile was in good agreement with the actual one. Moreover, the testing arrangement is beneficial in investigating the chain-soil-anchor interaction.     

An energy-controlling boundary condition for partial wave reflections in the mild slope equation

An energy-controlling technique to actively manage the reflective property of waves from solid boundary is presented. As linear waves propagate through an energy-controlling area, a reduction in wave heights occurs due to energy dissipation, which can be placed under direct control through the imaginary part of the wavenumber and phase velocity. Based on this relationship, the present study investigates a new method to control reflected waves with desired heights in the mild slope equation model. The method is validated through numerical tests for various reflection coefficients and the results confirm the promising use of energy-controlling boundary condition for partial wave reflections.     

Experimental and theoretical methods on determination of stress-dependent natural frequency of marine sedimentary clay

Excited by the vibration sources in dynamic engineering, the natural frequency and damping factor of the saturated marine sedimentary clay are key dynamic parameters that influence the responses under cyclic loads. Experimental and theoretical methods are proposed in this paper to analyze the natural frequency and the stress-dependent nonlinearity. The experimental method shows that the natural frequency of soils with specific stress state subjected to large cyclic shear strain can be estimated from the data of dynamic triaxial tests based on the amplitude–frequency response curve. Trial and error by the criterion from the half-power bandwidth method is used to determine the optimal fitting. The results of a theoretical study on the free vibration of soil layers are then presented to derive the analytic solution of natural frequency. In addition to the two frequency-independent elements (a lumped mass matrix and a stiffness matrix), the system’s equivalent damping coefficient matrix is iteratively determined based upon the forced vibration experimentally. Finally, the impacts of the resonance phenomenon on the dynamic shear modulus and hysteretic loop are discussed.     

基于广州实证的后改革时代中国城市住房权问题

社会公正是当前地理学的一个热点也是传统研究课题。住房问题往往是社会公正研究和实践的重点案例。中国住房制度改革之后,中国住房拥有率及人均住房面积均大幅度提高。然而在住房分配的过程中,不平等和不公正现象时有发生。如何区别平均和公正,如何看待中国住房改革中的社会公平与不公共存的现象,是了解中国新城市现象的一个重要途径。本文从近期西方政治学和地理学关于住房权（housing rights）、社会公正以及城市权利等文献的讨论出发,通过在广州的实证研究,运用结构方程模型和Pratt指数对涉及住房权主要方面的住房分异的发生途径及其公正性进行综合分析。实证研究说明在住房分配领域,存在着制度因素延续下来所导致的分异。一些制度因素已经转化为市场因素,直接或间接影响个人居住条件和住房所有权的获得,并通过后者进一步影响住房分异。与此同时,纯粹市场因素的影响则相对薄弱。上述结论有助于了解当前中国城市住房权的核心问题,认识社会公正的发展前景,以及探寻减少社会不公正的制度利器。     

基于Gassmann方程的流体替换方法在珊溪水库地震研究中的应用

引入基于Gassmann方程的流体替换方法,在分析地震波P波速度、波速比与岩石孔隙度和饱和度关系的基础上,应用于珊溪水库地震波速比和P波速度变化特征研究,得到:(1)珊溪水库震中区岩石始终处于接近水饱和的饱水状态,波速比和P波速度"下降-回升"的变化实质上反映了震中区岩石"孔隙度增大(饱和度减小)-饱和度增大"的变化,每一丛地震的波速比由极小值逐渐增大为极大值是由于岩石从不饱和状态变化到饱和状态;(2)根据每一丛地震波速比的变化,计算得到珊溪水库流体扩散率αs=1.06×104 cm2 s-1,该数值与美国南卡罗莱纳水库、巴西Acu水库、广东新丰江水库的流体扩散率基本一致;(3)震源区岩石孔隙度上限值为8.7%~2.0%,该数值与华东勘测设计研究院通过室内岩石物理力学性质试验测定的珊溪水库坝址区新鲜流纹斑岩的孔隙度平均值一致。     

First-principles derivation of reactive transport modeling parameters for particle tracking and PDE approaches

Both Eulerian and Lagrangian reactive transport simulations in natural media require selection of a parameter that controls the “promiscuity” of the reacting particles. In Eulerian models, measurement of this parameter may be difficult because its value will generally differ between natural (diffusion-limited) systems and batch experiments, even though both are modeled by reaction terms of the same form. And in Lagrangian models, there previously has been no a priori way to compute this parameter. In both cases, then, selection is typically done by calibration, or ad hoc. This paper addresses the parameter selection problem for Fickian transport by deriving, from first principles and D (the diffusion constant) the reaction-rate-controlling parameters for particle tracking (PT) codes and for the diffusion–reaction equation (DRE). Using continuous time random walk analysis, exact reaction probabilities are derived for pairs of potentially reactive particles based on D and their probability of reaction provided that they collocate. Simultaneously, a second PT scheme directly employing collocation probabilities is derived. One-to-one correspondence between each of D, the reaction radius specified for a PT scheme, and the DRE decay constant are then developed. These results serve to ground reactive transport simulations in their underlying thermodynamics, and are confirmed by simulations.