Numerical estimation of REV and permeability tensor for fractured rock masses by composite element method

The Monte Carlo method is used to generate parent stochastic discrete fracture network, from which a series of fractured rock samples of different sizes and orientations are extracted. The fracture network combined with a regular grid forms composite element mesh of the fractured rock sample, in which each composite element is composed of sub‐elements incised by fracture segments. The composite element method (CEM) for the seepage is implemented to obtain the nodal hydraulic potential as well as the seepage flow rates through the fractured rock samples. The application of CEM enables a large quantity of stochastic tests for the fractured rock samples because the pre‐process is facilitated greatly. By changing the sizes and orientations of the samples, the analysis of the seepage characteristics is realized to evaluate the variation of the permeability components, the existence of the permeability tensor and the representative element volume. The feasibility and effectiveness are illustrated in a numerical example. Copyright © 2008 John Wiley & Sons, Ltd.     

Biaxial test simulations using a packing of polygonal particles

The mechanical response of cohesionless granular materials under monotonic loading is studied by performing molecular dynamic simulations. The diversity of shapes of soil grains is modelled by using randomly generated convex polygons as granular particles. Results of the biaxial test obtained for dense and loose media show that samples achieve the same void ratio at large strains independent of their initial density state. This limit state resembles the so‐called critical state of soil mechanics, except for some stress fluctuations, which remain for large deformations. These fluctuations are studied at the micro‐mechanical level, by following the evolution of the co‐ordination number, force chains and the fraction of the sliding contacts of the sample. Copyright © 2007 John Wiley & Sons, Ltd.     

渤海湾盆地新生界生油岩系底界面温度分布

依据渤海湾盆地2000余口测温井的地温梯度数据、地层岩性描述、分层数据以及钻井资料,计算了该盆地各生油岩系底界面的温度。统计结果表明:渤海湾盆地沙河街组大部分凹陷区地层底界面温度介于90℃至150℃,目前仍具有大量生油的温度条件,而在隆起或一些凸起地区,该地层组段温度普遍小于90℃,未能达到生油的温度指标。东营组和孔店组地层也仍具有一定的生油温度条件。研究还表明:地层温度与地层界面埋深密切相关,温度随界面埋深的增大而升高,沉积厚度大的凹陷区地层界面温度大于沉积厚度小的凸起区或斜坡地带,说明地层界面埋深是决定地温高低的主要控制因素,而地温梯度对地层界面温度的影响相对较小。     

Flexural isostatic response of the Alps to increased Quaternary erosion recorded by foreland basin remnants,SE France

We test the hypothesis that flexural isostatic compensation of the mass removed by enhanced Quaternary erosion is responsible for uplift of the Western European Alps and their forelands. We use two well‐preserved and well‐dated (1.8 Ma) abandonment surfaces of foreland basin remnants in SE France (the Chambaran and Valensole plateaux) as passive benchmarks for tilting of the foreland. Estimating their initial slope from morphometric scaling relationships, we determine bulk post‐depositional tilting of 0.5–0.8% for these surfaces. The calculated isostatic response of the Alpine lithosphere to erosional unloading, using the method recently proposed by Champagnac et al. [Geology 35 (2007) 195–198] , yields a predicted tilting of 0.3–0.4% in the considered areas, explaining approximately half of the determined post‐depositional tilting. Such long‐term deformation being insensitive to cyclic loading/unloading because of glaciations, we suspect the other half to be related to as yet undetermined long‐wavelength and long‐lived tectonic process(es).     

冻融循环对塑料土工格栅拉伸性能的影响

在-25～ 25 ℃环境温度条件下, 采用水中冻融和空气中冻融两种方式, 经过50次冻融循环后, 在15 ℃、 0 ℃、 -10 ℃、 -20 ℃和-35 ℃工况下对土工格栅试样进行了拉伸性能的试验研究, 获得了聚乙烯塑料土工格栅的极限承载力、最大延伸率、 不同延伸率对应的应力值等数据. 试验结果表明, 聚乙烯塑料土工格栅在低温时抗拉性能明显提高, 从 15 ℃降至-35 ℃时, 5%延伸率对应的抗拉强度提高107%, 极限承载力提高19%. 同时, 聚乙烯塑料土工格栅极限延伸率随温度降低呈线形关系降低, 从 15 ℃降至-35 ℃时, 极限延伸率降低80%.     

墨西哥湾GC238区冷泉碳酸盐岩的微结构与石化微生物特征

采于墨西哥湾GC238海底天然气渗漏区浅表层的冷泉碳酸盐岩呈结核状产出,由方解石微晶和胶结物及少量的黄铁矿构成。胶结物由直径为0.1～0.5 m的方解石化的球体、卵形体、棒状体组成,充填于方解石晶体之间。冷泉碳酸盐岩结核下表面发育有由方解石化的球体、卵形体、棒状体组成的薄层,其中的一些球状集合体（约5 m）断面显示发育有核和外壳的层圈结构。黄铁矿呈草莓状,也具有相似的层圈结构。这种层状结构与活体古细菌被硫酸盐还原细菌包裹的层圈结构相似。样品中所保存的球体、卵形体、棒状体及其所组成的层圈结构可能是石化的甲烷氧化古细菌和硫酸盐还原细菌。     

柴达木盆地北缘侏罗系两类不同有机质丰度泥岩的 生物标志物特征对比研究

在柴达木盆地北缘地区,分别选取有机碳含量很低和较高的侏罗系泥岩样品,对比分析了它们在生物标志物组成上的差异。结果发现,高有机质丰度泥岩的生标组成与我国西北地区侏罗纪煤系有机质的特征差异不大,相比而言,低有机质丰度泥岩的正烷烃以前主峰为特征,Pr/Ph比值在1.0左右,三环萜烷和伽马蜡烷丰度较高,并在部分样品中检出了25 降藿烷系列。结合泥岩的有机岩石学特征,认为这些差异可能反映了泥岩沉积环境和生烃母质的不同:高有机质丰度泥岩的有机显微组分以相对弱还原条件下的形态有机质为主,包括藻类体、孢子体和角质体等,而低有机质丰度泥岩的有机显微组分以相对强还原条件下的矿物沥青基质为主,其母质可能来源于低等显微菌藻类。进一步通过对比不同有机质丰度泥岩,以及区内原油生标组成之间的相互关系,讨论了研究区的油源问题。     

基于MapGIS的土地利用数据库建设实践——以赣州市为例

随着"数字国土工程"的启动,国土资源信息化管理得到了各级政府的高度重视,而土地利用数据库的建设是"数字国土工程"中基础数据库的重要组成部分.文章通过分析桌面地理信息系统MapGIS的功能,结合赣州市土地利用现状图的数字化实践,总结出用MapGIS进行土地利用数据库建库的流程,并对数据库建设中应特别注意的质量问题进行了探讨.     

松散砂粒孔隙结构、孔隙分形特征及渗透率研究

多孔介质孔隙结构特征是决定流体微观运移机制的重要方面。文章用环氧树脂固结松散石英砂磨制出二维多孔介质薄片,并用孔隙结构图像处理方法提取出孔隙结构参数信息;根据孔隙结构参数,建立描述孔隙结构的分形几何模型,分析松散地层孔隙结构的分形特征;进行多孔介质渗透率的实验测量和分形计算。结果表明:多孔介质的孔隙率和孔径与粒径之间存在线性相关性,拟合度分别达到-0.976 5和0.996 6;在方格边长ε和含有孔隙的格子总数N(ε)的双对数坐标系中,lnε-lnN(ε)数据点近似成直线,且不同样品的孔隙分布分维数值比较接近;而在孔径r与大于该孔径的孔隙总数N(r)的双对数坐标系中,lnr-lnN(r)的数据点必须进行分段回归分析,以便能更好地反映孔隙结构的实际情况;孔隙分布分维数和孔径分维数与粒径也存在较高的对应关系;中砂渗透率的实验测得值是5.19×10-5mm2,孔径分维数分段回归法计算的大孔隙多孔介质渗透率为5.75×10-5mm2,测量值与计算值较为接近,孔径分维数可以计算多孔介质的渗透率。     

基于Mindlin理论Winkler地基参数的Kalman滤波识别

提出了Winkler地基参数的Kalman滤波识别方法。引入Mindlin理论后,推导了Winkler地基上板的控制微分方程。利用Fourier变换,推求了Winkler地基上简支板的Fourier闭式解。推导了Winkler地基参数的Kalman滤波方程,研究了Winkler地基参数的Kalman滤波识别的具体计算步骤。研究表明,运用Kalman滤波理论进行Winkler地基参数的识别,能有效地估计Winkler地基参数;Winkler地基参数的滤波收敛速度和精度受地基参数初始值选取和位移实测资料的影响,且Kalman滤波理论也可用于其他地基模型地基参数的识别。