A numerical study of mechanical behavior of a cement paste under mechanical loading and chemical leaching

Based on relevant experimental data of a petroleum cement paste under mechanical loading and chemical leaching, an elastic‐plastic model is first proposed by taking into account plastic shearing and pore collapse. The degradation of mechanical properties induced by the chemical leaching is characterized by a chemical damage variable which is defined as the increase of porosity. Both elastic and plastic properties of the cement paste are affected by the chemical damage. The proposed model is calibrated from and applied to describe mechanical responses in triaxial compression tests respectively on sound and fully leached samples. In the second part, a phenomenological chemical model is defined to establish the relationship between porosity change and calcium dissolution process. The dissolution kinetics is governed by a diffusion law taking into account the variation of diffusion coefficient with calcium concentration. The chemical model is coupled with the mechanical model, and both are applied to describe mechanical response of cement paste samples subjected to progressive chemical leaching and compressive stresses. Comparisons between experimental data and numerical results are presented.     

“Transverse Action” – A model benchmark exercise for numerical analysis of the Callovo-Oxfordian claystone hydromechanical response to excavation operations

The Callovo-Oxfordian claystone (COx) is considered as a potential geological formation to host an industrial radioactive waste repository in France. A detailed understanding of the thermo-hydro-mechanical (THM) behavior of the COx is a key issue for design of different repository structures and the safety calculations of the project. More particularly, numerical modeling of induced fracture networks around drifts excavated at the main level of the Andra’s Meuse/Haut-Marne Underground Research Laboratory (M/HM URL) and short and long term behavior of the COx around these drifts are of great interest. Several constitutive models have been developed/used in the framework of the R&D and simulation programs of Andra. A model benchmark exercise has been launched since 2012 to provide an overall view of the developed models regarding the in situ observations. In this view, two series of test cases, respectively at material point scale and at drift excavation scale are defined. Different kinds of constitutive models based on the elasto-visco-plasticity concept, continuum damage mechanics, the rigid block spring method and two-scale computational homogenized model (CHM) are used within this exercise. The obtained results show that accounting for material anisotropy and strain localization treatment techniques can improve the obtained results when elasto-visco-plastic models are used. Damage mechanics based approaches and methods accounting for discontinuities through discrete elements provide also interesting insights especially when fracturing processes must be modeled. However, more efforts are necessary to improve the robustness of these kinds of approaches in the complex context of COx response to excavation works.     

基于四面体格网的三维地质体建模算法

构建真三维模型是地质学、三维GIS、科学可视化及工程应用领域共同关心的课题。提出了基于TEN的三维离散数据的构模,详细说明了从离散的三维初始数据点建Delaunay四面体格网的算法和步骤。算例表明,该算法能够较好地生成任意形状三维地质体。     

安塞油田坪桥水平井区沉积微相三维建模研究

以安塞油田坪桥水平井区为例,探讨了应用水平井和露头资料及标点过程(示性点过程)方法进行沉积微相三维建模的思路与方法,同时对研究区三角洲前缘分流河道的定量结构模式进行了深入的研究。研究表明,坪桥水平井区三角洲分流河道单砂体宽度一般为40～100m左右。砂体结构呈两种模式,其一为“迷宫状”,砂体离散地分布于湖相泥岩中,其二为“拼合板状”,多个砂体在侧向上叠置而形成连片分布的砂体,但部分单砂体间可被泥质或钙质薄层侧积层所隔挡。应用水平井资料和露头资料,可有效地提取三维建模所必须的地质统计学参数。研究表明,采用标点过程方法,可建立符合地质实际的沉积微相三维模型。     

基于共反射点叠加的观测参数选择

在地震勘探数据处理中，通常采用的共中心点(CMP)叠加技术对于复杂的地质构造难以取得好的效果。若通过从道间距、覆盖次数等观测参数与目标层地层倾角、共反射点离散度的关系出发，并结合模型正演，对复杂构造地层的观测参数选取进行优化，使之能对目标层尤其是逆断层下降盘的阴影层段进行较充分的“照明”，即可得到共反射点叠加剖面，这对共反射点叠加是一项有益的研究。     

均匀介质中F-K反偏移法正演模拟

F-K反偏移法正演模拟是在给定地质模型后通过构建人工偏移场,再对人工偏移场进行F-K反偏移来实现。人工偏移场的构建是在地质模型上直接放上经正确拉伸和幅度校正后的地震脉冲子波。F-K反偏移中插值映射采用两点sinc插值,其能避免假同相轴的出现且具有较快的速度。奇异点处理使用泰勒近似展开。最后对凹陷模型进行了F-K反偏移法正演模拟实现,得到了预期的效果。     

起伏地表条件下地震波场数值模拟有限积分变换有限差分方法

针对起伏地表条件下的地震波数值模拟问题,提出将起伏的地表映射到一个规则的长方形网格坐标系中,并在此基础上推导出变换域中的波动方程。根据导出来的波动方程应用有限余弦变换有限差分方法进行地震波场的数值模拟,以解决起伏地表条件下的数值模拟。     

北京平原区土壤地球化学特征及影响因素分析

基于北京平原区土壤地球化学调查数据,运用数理统计、因子分析、相关性分析、回归分析等方法研究了土壤元素的横向和纵向分布特征,并探讨其影响因素.结果 表明:北京平原区土壤CaO、MgO、Na20、Hg、Cd、Sr、Ba、Sn等含量偏高,Sb、As、Th、W、Br、U、I、Mo、有机质等含量偏低.R型因子分析结果显示,F1因...     

基于GIS的三维水文地质建模及应用——以咸阳市为例

为了提升城市地质调查工作的总体服务水平,通过搭建城市地质综合GIS平台,采用人机交互式单元分割-拼接建模方法,以水文地质剖面为基础数据,加入水文地质钻孔、水文地质图、地表高程等多源数据作为约束,首次在咸阳地区构建了建模深度300 m、建模面积3714 km2的水文地质三维结构模型.本次建模的工作顺利完成为研究咸阳地区地下水文地质结构、地下水资源赋存特征提供了三维可视化的数据支撑,也为其他城市开展类似建模工作提供借鉴.     

河南尉氏县西部地质地球物理综合解译及地热资源远景区预测

尉氏县地热资源开发尚处于起步阶段,为适应大规模地热资源开发需求,当前急需探明该区地热资源赋存情况。研究区位于太康隆起次级构造单元中部鼻状构造带及尉氏—薛店次凹带结合部位,本文通过尉氏县城西部所采集的220个大地电磁测深点的二维及三维反演工作,结合二维反射地震剖面,分析了研究区地层、断裂展布构造,并对地热资源远景区进行了预测。维性分析结果表明研究区主构造走向为北西向,二维及三维大地电磁测深反演结果表明研究区4 km以浅地层由上到下可分为四层：新近系—第四系、三叠系、石炭系—二叠系、寒武系—奥陶系;综合二维反射地震及大地电磁测深资料,剖面位置处可识别出三条正断层。选择断裂带周边或穿过断裂带钻井,对水量有一定保证。