Robust numerical implementation of a 3D rate-dependent model for reservoir geomechanical simulations

A 3D elasto-plastic rate-dependent model for rock mechanics is formulated and implemented into a Finite Element (FE) numerical code. The model is based on the approach proposed by Vermeer and Neher (A soft soil model that accounts for creep. In: Proceedings of the International Symposium “Beyond 2000 in Computational Geotechnics,” pages 249-261, 1999). An original strain-driven algorithm with an Inexact Newton iterative scheme is used to compute the state variables for a given strain increment.The model is validated against laboratory measurements, checked on a simplified test case, and used to simulate land subsidence due to groundwater and hydrocarbon production. The numerical results prove computationally effective and robust, thus allowing for the use of the model on real complex geological settings.     

软粘土的各向异性和小应变条件下的本构模型(ASM)

以能量方程为基础推导出各向异性和小应变条件下弹塑性的本构方程,并将屈服面与修正桥模型和试验测得的屈服面进行了对比,表明新的屈服面比修正剑桥模型更接近试验结果,采用从初始应力到状态边界面的距离为参数的系数来修正硬化准则以模拟土体小应谱条件下的应力－应变特征。将新的本构模型编入有限元程序,对小应变试验进行了计算分析,计算结果表明,新的模型比修正剑桥模型能更好地反映土体的主要力学特性。     

Ethnic coexistence in a pluralistic campus environment

In Malaysia, ethnic segregation seems to have grown more and more pronounced at all levels of education, which may have in the main contributed to increasing occupational segregation by ethnicity when the graduates left to join the job market. Such trends may be disturbing given the effort the country has put in to promote interethnic understanding and reduce interethnic economic disparity since 1970. By critically investigating the dynamics of ethnic coexistence in the microcosm of the university campus environment, this paper provides statistical evidence to show how far the country has progressed in terms of ethnic relations since the watershed events of May 13, 1969; to what extent Malaysian multiethnic society is different now compared to the unmistakable racial “corporateness” and interethnic “separateness” that Furnivall observed in his classic study of 1948; and in what ways ethnic relations have been reshaped by three decades of affirmative action policies and the form of ethnic democracy adapted for this unique society.

Use of electrical resistivity methods for detecting subsurface fresh and saline water and delineating their interfacial configuration: a case study of the eastern Dead Sea coastal aquifers, Jordan

The alluvial aquifer is the primary source of groundwater along the eastern Dead Sea shoreline, Jordan. Over the last 20 years, salinity has risen in some existing wells and several new wells have encountered brackish water in areas thought to contain fresh water. A good linear correlation exists between the water resistivity and the chloride concentration of groundwater and shows that the salinity is the most important factor controlling resistivity. Two-dimensional electrical tomography (ET) integrated with geoelectrical soundings were employed to delineate different water-bearing formations and the configuration of the interface between them. The present hydrological system and the related brines and interfaces are controlled by the Dead Sea base level, presently at 410 m b.s.l. Resistivity measurements show a dominant trend of decreasing resistivity (thus increasing salinity) with depth and westward towards the Dead Sea. Accordingly, three zones with different resistivity values were detected, corresponding to three different water-bearing formations: (1) strata saturated with fresh to slightly brackish groundwater; (2) a transition zone of brine mixed with fresh to brackish groundwater; (3) a water-bearing formation containing Dead Sea brine. In addition, a low resistivity unit containing brine was detected above the 1955 Dead Sea base level, which was interpreted as having remained unflushed by infiltrating rain.     

考虑循环软化的非线性动力本构模型在FLAC3D中的实现

已有震害研究表明,震后边坡会因持续变形而破坏,且伴随着土体强度逐渐降低的现象,即土的循环软化行为。因此,有必要研究考虑循环软化的非线性动力本构模型以用于复杂条件下地震边坡稳定性分析。在已有的非线性动力本构模型基础上,提出了考虑循环软化的处理方法。同时,在FLAC3D平台上实现了本构模型二次开发,并通过了理论公式与已有文献中试验数据的验证。结果表明：计算出的骨干曲线与理论公式一致,且计算出的动剪切模量比及阻尼比与试验数据吻合较好,能够克服Hardin-Drnevich模型和Davidenkov模型在较大应变处（>0.01%）过高地估算阻尼比的缺陷;考虑了循环软化后,计算出的剪切强度有明显降低,且当遇到骨干曲线剪应力可以连续地过渡到软化后的主干曲线上,模型的收敛性较好。所开发的本构模型可为大应变条件下软土场地及边坡地震灾害评估提供支持。     

Occurrence of salt water above fresh water in dynamic equilibrium in a coastal groundwater flow system near De Panne, Belgium

A salt water lens is found above fresh water under the shore between Dunkerque (France) and Nieuwpoort (Belgium). This inverse density distribution is in a dynamic equilibrium. It develops due to the infiltration of salt water on the back shore during high tide. Under this salt water lens, water infiltrated in the adjacent dune area flows towards the sea and discharges at the seabed. This water quality distribution differs from the classic salt water wedge under fresh water described in the literature. Here, the evolution to this water quality distribution is simulated with a density dependent numerical model. A large tidal range, shore morphology and a permeable groundwater reservoir are the main conditions for the observed water quality distribution. By altering these conditions, intermediate water quality distributions between the classic salt water wedge and the one discussed here develop. Based on these simulations, it is expected that similar kinds of inverse density distribution could be present in a number of coastal areas, which have tides, a gently sloping shore and a permeable substratum.     

粗粒筑坝材料的动力变形特性

结合室内大型动三轴试验资料,对粗粒料的动力变形特性进行研究后认为,在复杂的高应力条件下振动过程中试验粗粒料由于体积缩小而变得更为密实,其动应力-应变关系与Hardin-Drnevich模型的双曲骨架曲线假定有明显差异,其阻尼特性也不符合Hardin假定;提出了反映试验过程中材料振动硬化特性的幂函数型动应力-应变关系模型和基于试验资料的粗粒料阻尼比估算公式;粗粒料的动残余变形特性与其初始密度、静应力状态、动荷载、振动持续时间等因素密切相关,建立了考虑初始固结围压影响的粗粒料残余应变与动剪应力之间的关系模型。可为强震区超高土石坝的深入研究提供可靠的理论依据。     

基于物态本构模型的路基动力反应分析

以动力固结理论和瞬态动力学理论为基础,将有效应力物态动本构关系引入到以动力反应与动力固结相耦合、静应力与动应力变化相耦合、孔压的产生、扩散和消散相耦合的动力三维瞬态动力学基本方程组中,形成了饱和土体有效应力物态地震反应分析的完整理论体系。利用有限元法对饱和路基进行了瞬态地震反应分析,较好地反映了土体在震动过程中剪 胀、剪缩的实际性态,比以往引进孔压模型和静-动交替分析计算更为合理。     

非饱和黄土的弹塑性软化本构模型

通过对陕西杨陵地区的非饱和黄土进行等应力比压缩试验和常规三轴试验,在对试验结果的应力-应变曲线特性进行分析的基础上,由试验结果直接确定屈服函数,根据Drucker公设,即相适应的流动法则,选择合适的硬化参数,建立了非饱和黄土的弹塑性本构模型,并将模型计算结果与试验结果进行比较。结果表明,该模型能够较好地模拟非饱和黄土与固结应力有关的应变软化现象,计算的变形规律与试验结果基本吻合。     

泥岩渗流-应力耦合蠕变损伤模型研究(Ⅱ):数值仿真和参数反演

进一步分析了第Ⅰ部分[1]提出的泥岩渗流-应力耦合蠕变损伤模型。在连续损伤力学理论和比奥(Biot)理论的基础上,导出了考虑渗流-应力-损伤耦合的蠕变损伤有限元格式,建立了弹性预测、塑性修正、损伤修正-渗透系数修正的数值分析框架,编制了非线性有限元分析程序。根据监测的衬砌长期变形数据,采用优化反分析法获得了蠕变损伤模型中的待定参数,并应用于比利时核废料库施工过程中泥岩巷道围岩渗流-应力耦合过程、损伤演化以及长期稳定性分析,研究结果表明,泥岩开挖后渗透性明显增大,约为原岩的120倍,蠕变效应导致泥岩裂隙和渗透性自愈合,约3.5年后渗透性基本恢复到原岩的数量级,围岩中部的蠕变明显大于顶部和底部。研究成果对软岩隧洞长期稳定性的预测与预报具有一定的参考意义。