Spatial distribution of block falls using volumetric GIS–decision-tree models

Block falls are considered a significant aspect of surficial instability contributing to losses in land and socio-economic aspects through their damaging effects to natural and human environments. This paper predicts and maps the geographic distribution and volumes of block falls in central Lebanon using remote sensing, geographic information systems (GIS) and decision-tree modeling (un-pruned and pruned trees). Eleven terrain parameters (lithology, proximity to fault line, karst type, soil type, distance to drainage line, elevation, slope gradient, slope aspect, slope curvature, land cover/use, and proximity to roads) were generated to statistically explain the occurrence of block falls. The latter were discriminated using SPOT4 satellite imageries, and their dimensions were determined during field surveys. The un-pruned tree model based on all considered parameters explained 86% of the variability in field block fall measurements. Once pruned, it classifies 50% in block falls’ volumes by selecting just four parameters (lithology, slope gradient, soil type, and land cover/use). Both tree models (un-pruned and pruned) were converted to quantitative 1:50,000 block falls’ maps with different classes; starting from Nil (no block falls) to more than 4000 m³. These maps are fairly matching with coincidence value equal to 45%; however, both can be used to prioritize the choice of specific zones for further measurement and modeling, as well as for land-use management. The proposed tree models are relatively simple, and may also be applied to other areas (i.e. the choice of un-pruned or pruned model is related to the availability of terrain parameters in a given area).     

CO2 saturation estimates at Sleipner (North Sea) from seismic tomography and rock physics inversion

CO₂ saturations are estimated at Sleipner using a two-step imaging workflow. The workflow combines seismic tomography (full-waveform inversion) and rock physics inversion and is applied to a two-dimensional seismic line located near the injection point at Sleipner. We use baseline data (1994 vintage, before CO₂ injection) and monitor data that was acquired after 12 years of CO₂ injection (2008 vintage). P-wave velocity models are generated using the Full waveform inversion technology and then, we invert selected rock physics parameters using an rock physics inversion methodology. Full waveform inversion provides high-resolution P-wave velocity models both for baseline and monitor data. The physical relations between rock physics properties and acoustic wave velocities in the Utsira unconsolidated sandstone (reservoir formation) are defined using a dynamic rock physics model based on well-known Biot–Gassmann theories. For data prior to injection, rock frame properties (porosity, bulk and shear dry moduli) are estimated using rock physics inversion that allows deriving physically consistent properties with related uncertainty. We show that the uncertainty related to limited input data (only P-wave velocity) is not an issue because the mean values of parameters are correct. These rock frame properties are then used as a priori constraint in the monitor case. For monitor data, the Full waveform inversion results show nicely resolved thin layers of CO₂–brine saturated sandstones under intra-reservoir shale layers. The CO₂ saturation estimation is carried out by plugging an effective fluid phase in the rock physics model. Calculating the effective fluid bulk modulus of the brine–CO₂ mixture (using Brie equation in our study) is shown to be the key factor to link P-wave velocity to CO₂ saturation. The inversion tests are done with several values of Brie/patchiness exponent and show that the CO₂ saturation estimates are varying between 0.30 and 0.90 depending on the rock physics model and the location in the reservoir. The uncertainty in CO₂ saturation estimation is usually lower than 0.20. When the patchiness exponent is considered as unknown, the inversion is less constrained and we end up with values of exponent varying between 5 and 20 and up to 33 in specific reservoir areas. These estimations tend to show that the CO₂–brine mixing is between uniform and patchy mixing and variable throughout the reservoir.     

计算机辅助药物设计中用于建模的计算方法研究进展

通过检索国内外相关领域的大量研究文献,并结合作者的研究经验,综述了计算机辅助药物设计中计算方法的发展及其应用的最新进展,并对计算方法的选用原则进行了探讨,以期为从事CADD工作者合理选择建模计算方法提供有价值的参考。     

Vertical distribution of bacteria and archaea in a CO2-rich meromictic lake: A case study of Lake Monoun

In order to assess the diversity and spatial distribution of bacterial and archaeal communities in a CO₂-rich and meromictic lake, samples were collected along the water column of Lake Monoun and analyzed using PCR-DGGE and quantitative analyses of the 16S rRNA gene. The retrieved sequences were affiliated to 6 bacterial phyla and two archaeal phyla which plausible environmental functions map with the physico-chemical parameters of the lake. Unclassified sequences were also detected. This suggests heterogeneity in community composition and existence of potential candidate divisions. For instance, amongst the bacterial sequences, 18.2% matched with methanotrophic bacteria of the order methylococcales and amongst the archaeal sequences, 16.6% matched with methanogenic species of the order methanomicrobiales. Hence, evidencing the existence of methane-related prokaryotes in the lake, a finding that would play a key role in our understanding of the methane puzzle of Lake Monoun. Other groups capable of a wide range of metal and nutrients transformations were also detected, as well as those of unknown functions. The layering of microbial communities also appeared to directly or indirectly depend on oxygen availability. DGGE and qPCR analyses both suggested a scarcity of archaea in the surface samples. Furthermore, qPCR revealed that bacteria were numerically more important than archaea in all the samples. The general distribution along the water column indicated that archaeal abundance increases with depth.     

瞬变电磁场资料的联合时-频分析解释

本文充分利用瞬变电磁（TEM）信号包含的信息,采用联合时-频分析方法解释了TEM资料。文章选取了Wigner-Ville分布和Gabor展开的联合时-频分布（或表示）,将一维时间域信号拓展到时-频二维平面上,不仅将不同地电断面的响应成功地分开来,而且所显示的信号特征与地层电性结构之间有着明确的物理意义,从时间和频率两个方面同时描述了瞬变涡流场在地层中的激发和衰减的过程。利用时间窗数据和频率窗数据所进行的联合时-频反演,缩小了拟合差,提高了定量解释精度。     

DECOVALEX Project: from 1992 to 2007

The DECOVALEX project is a unique international research collaboration, initiated in 1992, for advancing the understanding and mathematical modelling of coupled thermo-hydro-mechanical (THM) and thermo-hydro-mechanical-chemical (THMC) processes in geological systems—subjects of importance for performance assessment of radioactive waste repositories in geological formations. From 1992 up to 2007, the project has made important progress and played a key role in the development of numerical modelling of coupled processes in fractured rocks and buffer/backfill materials. The project has been conducted by research teams supported by a large number of radioactive-waste-management organizations and regulatory authorities, including those of Canada, China, Finland, France, Japan, Germany, Spain, Sweden, UK, and the USA. Through this project, in-depth knowledge has been gained of coupled THM and THMC processes associated with nuclear waste repositories, as well as numerical simulation models for their quantitative analysis. The knowledge accumulated from this project, in the form of a large number of research reports and international journal and conference papers in the open literature, has been applied effectively in the implementation and review of national radioactive-waste-management programmes in the participating countries. This paper presents an overview of the project.     

地震散射波的高精度数值模拟与振幅分析

随着地震勘探要求和技术的不断提高,散射波地震勘探方法逐渐引起人们的重视,尤其在构造复杂地区的高精度勘探中,散射波理论具有广泛的应用前景.本文在讨论了现阶段几种散射波数值模拟方法的利弊后,设计了一种基于微扰论的FK域积分法,在散射场的二次震源和空间能量衰减处理两方面做了改进,其计算精度和效率均得到了提高.此外,采用该算法针对散射波动力学进行了正演模拟试算,从扰动介质空间分布和扰动量两个角度定量地分析了散射波振幅动力学方面的特征.在此基础上得到了五个经验常数,并给出了对应的计算公式和物理意义,通过这五个常数可以划分散射波不同形态的动力学特征.     

埋藏岩溶洞穴垮塌深度定量图版及其在碳酸盐岩缝洞型储层地质评价预测中的意义

通过对岩溶洞穴垮塌的岩石力学原理及塔里木盆地奥陶系实钻资料的分析,基本明确了洞穴埋藏垮塌的控制因素是岩石抗弯强度、洞穴尺度、洞穴距风化暴露面的距离等,并得到了洞穴垮塌深度的定量图版。该图版对于钻前预测现今埋藏于地下的碳酸盐岩洞穴是否已经垮塌、洞穴埋藏演化史精确恢复以及对裂缝—洞穴型储层的评价预测具有较强的实践意义。     

定量地震地貌学研究进展

定量地震地貌学以三维地震数据为基础,通过对地貌特征参数的定量分析来了解盆地历史、沉积过程和填充构架,在沉积体系定量分析和储层表征方面为一新方向.其主要研究内容包括系统地质研究、地震属性分析、地貌特征定量分析、沉积单元交汇分析和沉积体系影响因素研究等五方面,核心为多数据融合及沉积单元形态定量数据的利用.定量地震地貌学在冲积河道、三角洲、浊积水道、水下扇储层特征及结构研究方面已取得较大突破,在地层沉积定量研究、储集体规模预测方面具有较大优势和前景.     

青藏高原北部地区地震基本活动状态定量评价

基于区域构造背景、应力场状态及动力学环境等,将青藏高原北部地区划分为4个相对独立的统计单元,即祁连山地震带、甘东南地区、柴达木一共和地块及库玛地震带.通过地震频度和b值拟合确定了各个地区不同下限震级的地震目录完整性起始时间;确定了各构造单元中强以上地震活跃与平静交替活动的特征;定量计算了在平均状态中强以上地震活动特征参数及平静阶段与活跃阶段地震活动特征参数;分别获得了各个区域在平均、活跃与平静状态下的小震活动状态参数.