A New Methodology for Singularity Modelling in Flow Simulations in Reservoir Engineering

In this paper, numerical methods for flow simulation in the well vicinity are discussed. Flux truncation errors are analyzed on gridblocks near the well. Due to the singularity of the well, standard numerical schemes are not efficient for near well flow simulations. To improve simulation accuracy, a new methodology, using a change of coordinates to make near well pressure linear, is presented for the singularity modelling. Based on this new approach, two numerical schemes, a two-point flux approximation scheme and a multi-point one, are proposed. These schemes are accurate for the near well modelling and are suitable for any kinds of gridblocks.     

土钉支护结构优化设计的有限元分析

考虑土钉-土体的相互作用,建立三维有限元模型,分析计算土钉支护结构的内力和位移,提出了变形控制和应力控制相结合的优化设计方法。采用遗传算法,以应力和位移为约束条件,以工程造价为目标函数,建立了土钉支护结构的优化设计数学模型,并编制了相应的计算程序,通过算例分析,并基于规范法的优化设计结果相比较,得出本文的设计方法更为合理的参考结论。     

Nuclear Waste Disposal Simulations: Couplex Test Cases

We give some results obtained for the Couplex test cases proposed by the ANDRA. In this paper our aim is twofold. Firstly, to compute the release of nuclides out of the repository by concentrating on the 3D near field (Couplex 2). The simulation of the transport phenomena takes into account the dissolution of the glass containers and congruent emissions of the radio-nuclides including filiation chains and some simplified chemistry. Secondly, it is to use the near field computations in order to simulate the nuclide migrations in a 2D far field (Couplex 3). Coupling in between the two simulations takes into consideration the periodicity of the disposal modules and the geometry of the repository described in Couplex 1. The mixed finite element and discontinuous Galerkin methods are used to solve the convection–diffusion equations. In order to handle the nonlinear precipitation/dissolution term, we developed a new iterative technique that combines Picard and Newton–Raphson methods.     

鞍山市西郊区水位预测

根据鞍山市西郊区水文地质条件，建立了双层渗流二维平面系统的数学模型：应用Galerkin有限元法求解模型，从模型的计算结果可见该区观测点水位绝对误差较小，达到了规范要求的计算精度，同时，在此区内有一个大漏斗，在东南的铁西一带也有一个小漏斗。由于铁西一带是主要开采水源区，开采水量较大，所以漏斗逐渐向外扩大且水位有明显下降的趋势。这一变化规律与实际水位一致，说明计算水位基本上反映了实际水位空间变化的规律。因此根据模拟参数利用此模型对该地区进行了水位预测，共预测了四年的水位状况，并依据水位的预测结果对该区的水资源开发利用提出了合理建议。     

Experimental and numerical studies of rock salt strain hardening

The work presented in this paper comes as part of a research program dealing with the thermomechanical behaviour of rock salt. It aims to study laboratory and in-situ long-term behaviour by means of creep tests with deviator and temperature changes. The laboratory results, using a triaxial multi-stages creep tests, highlighted the strain hardening character of rock salt. Furthermore, the in-situ results, using a borehole dilatometer multi-step creep test, have shown that the drilling is carried out in a weakly stressed pillar. The interpretation of the laboratory results, using the J.LEMAITRE law, did not indicate full agreement with all the test results. As a result a ‘double’ J.LEMAITRE model, which takes into account a double strain hardening variable, has been put forward. The validation of this model on the laboratory creep tests is very satisfactory. Furthermore, the activation energy seems satisfactory to represent the influence of the temperature. The in-situ behaviour modelling is clearly more complex than the modelization based on laboratory tests. In fact, it seems that if the rock salt behaviour is maintained by J.LEMAITRE law, it is necessary to vary with the stress, at least, one of the parameters assumed constant in the basic law.     

深海管道钢的腐蚀疲劳裂纹扩展特性试验研究

采用侧壁开槽型CT试样,分别在空气和3.5%NaCl溶液中进行腐蚀疲劳裂纹扩展试验,研究海底管道用钢X65在海水腐蚀环境下的腐蚀疲劳裂纹扩展特性。为保证试样既能满足腐蚀疲劳试验机的夹具要求,又可控制裂纹沿直线方向扩展,采用侧壁开槽技术对CT试样进行改进;利用有限元软件ABAQUS,建立三维模型,验证了对于侧壁开槽型CT试样,可以采用ASTM标准推荐的应力强度因子表达式对改进后CT试样开展数值计算。利用YYF-50腐蚀疲劳试验装置开展试验,试验结果表明:和惰性环境相比,腐蚀环境不仅会加速X65钢的裂纹扩展速率,还会降低腐蚀疲劳体系下裂纹扩展的门槛值;在腐蚀体系下存在裂纹扩展的稳定阶段,可采用Paris公式对其腐蚀疲劳裂纹扩展速率进行预测。     

地质条件复杂地区瑞雷波勘探正演模拟

为研究复杂地质条件对瑞雷面波频散曲线的影响,采用有限差分地震数值模拟方法,对3种模型进行正演模拟,得出了某一时刻的波场快照与单炮记录。对模拟记录进行频散分析的结果表明：单层均匀介质没有频散现象;起伏地表两层均匀介质的频散曲线相速度只和介质本身的参数有关,不受地形的影响;水平地表含软夹层介质中,瑞雷波频散现象非常明显。     

A Kinematic Thermal Model for Descending Slabs with Velocity Boundary Layers: A Case Study for the Tonga Subducting Slab

For the purpose of investigating the influence of metastable olivine (MO) phase transformations on both deep seismicity and stagnation of slabs, we constructed a 2-dimensional finite element thermal model for a 120 Ma-old 50o dipping oceanic lithosphere descending at 10 cm/yr with velocity boundary layers, which would mitigate the interference of constant velocity field for the slab. The resulting temperatures show that most of intermediate and deep earthquakes occurring within the Tonga slab are occurring inside the 800oC and 1200oC isotherm, respectively. The elevation of olivine transformation near ~410 km and respective persistence of metastable olivine and spinel within the transition zone and beneath 660 km would thus result in bimodal positive, zonal, negative density anomalies, respectively. These results together with the resulting pressure anomalies may reflect the stress pattern of the Tonga slab: (i) slab pull force exerts above a depth of ~230 km; (ii) MO existence changes the buoyancy force within the transition zone and facilitates slab stagnation at a depth of 660 km; (iii) as the subducting materials accumulated over 660 km, deepest earthquakes occur due to MO transformation; (iv) a flattened ‘slab’ may penetrate into the lower mantle due to the density increment of Sp transformation.     

基于拟静力法的弹簧刚度计算方法讨论

本文讨论了拟静力法在动力状态下的地基基础协同分析中的应用,通过分析假定动力状态下岩土体的剪切模量和剪切应变的关系,将地震力转化为等效的惯性力系。对于不同类型的岩土体,剪切模量G和剪切应变Y的对应关系是有区别的,且地震过程是不可预知的,过程无法模拟。研究表明,采用线性方法简化地震力衰减过程,仍可以得到较准确的计算结果。本文还采用FLAC-3D软件对上述关系进行了模拟,通过工程实例计算验证其可行性,提供了一套在动力状态下计算弹簧刚度参数的方法。     

PARAMETER AND VARIANCE ESTIMATION IN GEOTECHNICAL BACKANALYSIS USING PRIOR INFORMATION

A probabilistic framework to perform inverse analysis of geotechnical problems is presented. The formulation allows the incorporation of existing prior information on the parameters in a consistent way. The method is based on the maximum likelihood approach that allows a straightforward introduction of the error structure of field measurements and prior information. The difficulty of ascribing definite values to the uncertainties associated with the various types of observations is overcome by including the corresponding variances in the set of parameters to be identified. The inverse analysis results in a minimization problem that is solved by coupling the optimization technique to the finite element method. Two examples are presented to illustrate the performance of the method. The first one corresponds to a synthetic case simulating the excavation of a tunnel. Young's modulus, K₀ value and measurements variances are identified. The second case concerns the excavation of a large underground cavern in which again Young's modulus and K₀ are identified. It is shown that introduction of prior information permits the estimation of parameters more consistent with all available informations that include not only monitored displacements but also results from in situ tests carried out during the site investigation stage.