A Connectivity Index for Discrete Fracture Networks

Connectivity is an important measure for assessing flow transport in rock, especially through fractures. In this paper, rock fracture systems are modelled by a discrete fracture model simulated by a marked point process. A connectivity index is then introduced to quantify the connectivity between any two points in space. Monte Carlo simulation is used to evaluate the connectivity index for stationary cases and relationships between the connectivity index and the parameters of the discrete fracture model are analysed. The average number of intersections per fracture, X_f, and the fracture intensity, P12 (P32), are calculated and the relationships between these parameters and the connectivity index are investigated, concluding that X_f is the more suitable parameter for the classification of rock mass flow properties. The relationships between the connectivity index and the percolation state of the fractured medium are also discussed. An edge correction is briefly discussed and a practical example is used to demonstrate the method of computing the connectivity index.     

Short-term beach–dune sand budgets on the north sea coast of France: Sand supply from shoreface to dunes, and the role of wind and fetch

Three experimental plots, covering the transition from the upper beach to the dune, on the North Sea coast of France were monitored at various intervals over a period of 18–24 months via high resolution terrain surveys in order to determine inter-site sand budget variability, as well as patterns and processes involved in sand exchanges between the upper beach and dune. The wind regime consists of a fairly balanced mix of moderate (80% of winds are below 8 m/s) onshore, offshore and shore-parallel winds. Sustained dune accretion over several years depends on the periodic local onshore welding of shoreface tidal banks that have developed in the storm- and tide-dominated setting of the southern North Sea. The only site where this has occurred in the recent past is Calais, where bank welding has created a wide accreting upper beach sand flat. At this site, significant sand supply from the subtidal sand bank reservoir to the upper beach flat occurred only once over the 18-month survey following a major storm. The bulk of the sand deposited over this large flat is not directly integrated into the adjacent embryo dunes by onshore winds but is progressively reworked in situ into developing dunes or transported alongshore by the balanced wind regime, thus resulting in alongshore stretching of the embryo dune system. The Leffrinckoucke site near Belgium shows moderate beach–dune mobility and accretion, while the Wissant site exhibits significant upper beach bedform mobility controlled by strong longshore currents that result in large beach budget fluctuations with little net budget change, to the detriment of the adjacent dunes. Accretion at these two sites, which are representative of the rest of the North Sea coast of France, is presently constrained by the absence of a shore-attached sand bank supply reservoir, while upper beach–dune sand exchanges are further limited by the narrow wave-affected upper beach, the intertidal morphology of bars and troughs which segments the aeolian fetch, and the moderate wind energy conditions. The balanced wind regime limits net sand mobilisation in favour of either the beach or the dune, and may explain the relatively narrow longshore morphology of the dune ridges bounding this coast.     

A greenhouse study on root dynamics of cactus pears, Opuntia ficus-indica and O. robusta

Over the last 10 years a great interest in spineless cactus pear was shown in the drier areas in terms of both fresh fruit and fodder production. However, there is a lack of knowledge on quantitative data on root dynamics of these plants needed to fully understand its potential under water limiting conditions. This study aimed at quantifying the effects of water stress on the growth of tap roots, side roots and rain roots of the species Opuntia ficus-indica (L.) Miller (cultivar Morado—green cladode) and O. robusta Wendl. (cultivar Monterey—blue cladode). One-year-old cladodes were planted in root boxes and pots (2002/2003 season) that were kept in the greenhouse at day/night temperatures of 25–30 °C/15–18 °C. Placing the cladodes flat on the soil, more areoles came in contact with the soil and therefore more roots developed in both species with an average of only 3.4% areole complexes not rooting. Each areole complex formed on average 3 roots. The highest daily tap root growth was 42 and 36 mm for O. ficus-indica and O. robusta, respectively. Tap root growth increased in the morning with water stress for both species but decreased in the afternoon. Side root growth increased with water stress, with that of O. robusta more per tap root than O. ficus-indica. O. robusta showed a finer root system than O. ficus-indica. The side roots grew as much as 8 and 5 mm per day for O. ficus-indica and O. robusta, respectively. Whitish rain roots developed on the established roots within the first hour after rewetting the soil and grew for only 3 days. Rain roots grew up to 7 and 5 mm within a day for O. ficus-indica and O. robusta, respectively. Considering all studied aspects of their roots systems, O. robusta appears to be better adapted to drought (less sensitive to water stress) than O. ficus-indica.     

Optimization of SAGD Well Elevation

The application of steam-assisted gravity drainage (SAGD) to recover heavy oil sands is becoming increasingly important in the northern Alberta McMurray Formation because of the vast resources/reserves accessible with this mechanism. Selecting the stratigraphic elevations of SAGD well pairs is a vital decision for reservoir evaluation and planning. The inherent uncertainty in the distribution of geological variables significantly influences this decision. Geostatistical simulation is used to capture geological uncertainty, which is used can be transformed into a distribution of the best possible well pair elevations. A simple exhaustive calculation scheme is used to determine the optimum stratigraphic location of a SAGD well pair where the recovery R is maximized. There are three basic steps to the methodology: (1) model the uncertainty in the top continuous bitumen (TCB) and bottom continuous bitumen (BCB) surfaces, (2) calculate the recovery at all possible elevation increments within the TCB and BCB interval, and (3) identify the elevation that maximizes R. This is repeated for multiple TCB/BCB pairs of surfaces to assess uncertainty. The methodology is described and implemented on a subset of data from the Athabasca Oilsands in Fort McMurray, Alberta.     

Effect of sampling and linkage on fault length and length–displacement relationship

The power-law exponent (n) in the equation: D=cL ⁿ , with D = maximum displacement and L = fault length, would be affected by deviations of fault trace length. (1) Assuming n=1, numerical simulations on the effect of sampling and linkage on fault length and length–displacement relationship are done in this paper. The results show that: (a) uniform relative deviations, which means all faults within a dataset have the same relative deviation, do not affect the value of n; (b) deviations of the fault length due to unresolved fault tip decrease the values of n and the deviations of n increase with the increasing length deviations; (c) fault linkage and observed dimensions either increase or decrease the value of n depending on the distribution of deviations within a dataset; (d) mixed deviations of the fault lengths are either negative or positive and cause the values of n to either decrease or increase; (e) a dataset combined from two or more datasets with different values of c and orders of magnitude also cause the values of n to deviate. (2) Data including 19 datasets and spanning more than eight orders of fault length magnitudes (10⁻²–10⁵ m) collected from the published literature indicate that the values of n range from 0.55 to 1.5, the average value being 1.0813, and the peak value of n _d (double regression) is 1.0–1.1. Based on above results from the simulations and published data, we propose that the relationship between the maximum displacement and fault length in a single tectonic environment with uniform mechanical properties is linear, and the value of n deviated from 1 is mainly caused by the sampling and linkage effects.     

鄱阳盆地构造—沉积特征及其演化史

 鄱阳盆地是发育在江南—九岭和怀玉—官帽基底拆离造山带上的白垩-古近纪张扭性断陷盆地,明显受近南北走向的赣江走滑断裂系和前白垩系逆冲断裂系晚期反转的双重因素控制,呈“两坳夹一隆”的区域构造格局。受赣江、进贤—石门街走滑断裂夹持的南昌凹陷,呈近南北向凹凸相间展布格局。进贤—石门街走滑断裂以东地区,“南断北超”型断陷自南而北由北东走向渐转为北北东向展布,其源于晚印支-燕山期北东向逆冲断裂在晚期的反转作用,西缘和北部受赣江断裂系东枝走滑断裂影响明显而发生转向。鄱阳盆地经历了早白垩世晚期(冷水坞组沉积期)拉分盆地、晚白垩世早-中期(周家店组—南雄组沉积早期)走滑张扭断(坳)陷盆地、晚白垩世晚期(南雄组沉积中-晚期)走滑伸展盆地和古近纪坳陷盆地等演化阶段,东南部断陷区(江埠—二甲村凹陷)由于喜马拉雅期的隆升剥蚀而只残存周家店组早期坳陷和周家店组中期—南雄组中期断陷两个发育阶段的沉积。     

裂缝性岩石裂缝扩展的光弹构造物理模拟——以南华北盆地户部寨沙河阶组为例

以南华北盆地户部寨第三系沙河阶组粉砂岩低渗透储层为例,总结了其中发育的5种类型构造裂缝组合,即:①南北向雁形张节理组成的共轭剪节理带;②东西向雁形张节理组成的共轭剪节理带;③缠结状多组剪切构造裂缝;④组合状交切的多组剪切构造裂缝及张裂缝;⑤散状无序排列的构造裂缝。在此基础上,建立光弹性试验构造物理模型,沿现应力场96°方向,均匀地对该模型逐步增加外载力。结果表明,沿各组合裂缝尖端分别形成对称于裂缝走向、且主轴垂直于裂缝的Ⅰ型等色线条纹;对称于裂缝走向、且主轴平行于裂缝走向的Ⅱ型条纹和斜交于裂缝走向的复合型(Ⅲ型)条纹。第①类组合裂缝尖端以Ⅰ型等色线条纹的拉张应力和张剪应力为主;第②类组合裂缝尖端以复合型(Ⅲ型)等色线条纹的挤压应力为主;第③、④和⑤类组合裂缝尖端的应力分布较为复杂。当模型外载荷接近1000kg时,裂缝尖端起裂,并开始发育显微裂缝,而当外载荷达1100-1200 kg时,沿裂缝尖端部位形成宏观裂缝,并开始扩展。5种类型的裂缝组合扩展发育程度由强至弱依次为:①,③,⑤,②,④。其中第①类共轭裂缝组合起裂最早,扩展最快,裂缝间的连通趋势最好,说明规则分布的构造裂缝组合,在垂直于外施加载衙时更易达到连通的目的。     

渗透系数空间变异性研究

水文地质参数的空间变异性是随机理论研究的基础,而渗透系数是最为重要的水文地质参数。国外有关渗透参数空间变异性的研究工作已开展很多,但渗透系数究竟服从什么分布目前尚无确切答案。利用Borden含水层试验数据,对渗透系数的空间变异性进行探讨,结果表明若处理方法得当,渗透系数应服从对数正态分布。同时,还对今后野外开展含水层渗透系数空间变异性试验研究提出几点建议。     

地质雷达及浅层地震在泥质粉砂岩中断裂带上的曲线特征

探讨了如何利用地质雷达及浅层地震有关原理判释泥质粉砂岩中隐伏断裂带的性质和特征。主要的工作方法是通过现场实测，分析地质雷达反映在张性断裂带和压性断裂带上的曲线特征，分析浅层地震相遇时距曲线在断裂带上的曲线特征。后者的主要特征有：（1）在断裂带上初至时间延长；（2）在断裂带上左右支时距曲线上下跳动的幅度均等；（3）在断裂带上跳动的长度不一；（4）在断裂带左右支时距曲线上下跳动或反相；（5）时距曲线呈喇叭口。阐明了地质雷达和浅层地震的原理及其在泥质粉砂岩中勘查断裂构造的效果