华北地台中元古代碳酸盐岩中的微生物成因构造及其生烃潜力

华北地台中元古界主要由环潮坪石英砂岩、浅海碳酸盐岩和浅海—泻湖相暗色页岩3种沉积相组合构成,以陆表海浅水碳酸盐岩占主导。碳酸盐岩中除含有丰富的微古植物、宏观藻类和微生物建隆外,还发育大量的微生物成因构造(MISS)和微生物诱发的碳酸盐沉淀(MMCP)。微生物席和MISS构造在高于庄组上部(约1.6 Ga)和雾迷山组下部(约1.45 Ga)碳酸盐岩中尤为发育,表明活跃的微生物活动和高有机质产量。在石化微生物席中,发现有丝状、球状细菌化石和草莓状黄铁矿;围岩中发现有针状文石、花瓣状重晶石、放射状菱铁矿、铁白云石和葡萄状碳酸盐胶结物等多种自生碳酸盐矿物,指示甲烷厌氧氧化(AOM)导致的自生碳酸盐沉淀。中元古代的温暖气候和海洋分层、缺氧、硫化条件有利于微生物的高生产量和高有机质埋藏率。气隆构造和核形石状碳酸盐结核反映浅埋藏条件下活跃的成烷作用和甲烷排放,围岩和MMCP中富沥青质。华北地台中元古界富微生物席碳酸盐岩有良好的生烃潜力,有可能形成重要的烃源岩。据微生物席、MISS构造及MMCP的研究,初步估算华北地台中元古代碳酸盐岩的概略生烃潜力约为10×108t石油当量。     

Integration of local–global upscaling and grid adaptivity for simulation of subsurface flow in heterogeneous formations

We propose a methodology, called multilevel local–global (MLLG) upscaling, for generating accurate upscaled models of permeabilities or transmissibilities for flow simulation on adapted grids in heterogeneous subsurface formations. The method generates an initial adapted grid based on the given fine-scale reservoir heterogeneity and potential flow paths. It then applies local–global (LG) upscaling for permeability or transmissibility [7], along with adaptivity, in an iterative manner. In each iteration of MLLG, the grid can be adapted where needed to reduce flow solver and upscaling errors. The adaptivity is controlled with a flow-based indicator. The iterative process is continued until consistency between the global solve on the adapted grid and the local solves is obtained. While each application of LG upscaling is also an iterative process, this inner iteration generally takes only one or two iterations to converge. Furthermore, the number of outer iterations is bounded above, and hence, the computational costs of this approach are low. We design a new flow-based weighting of transmissibility values in LG upscaling that significantly improves the accuracy of LG and MLLG over traditional local transmissibility calculations. For highly heterogeneous (e.g., channelized) systems, the integration of grid adaptivity and LG upscaling is shown to consistently provide more accurate coarse-scale models for global flow, relative to reference fine-scale results, than do existing upscaling techniques applied to uniform grids of similar densities. Another attractive property of the integration of upscaling and adaptivity is that process dependency is strongly reduced, that is, the approach computes accurate global flow results also for flows driven by boundary conditions different from the generic boundary conditions used to compute the upscaled parameters. The method is demonstrated on Cartesian cell-based anisotropic refinement (CCAR) grids, but it can be applied to other adaptation strategies for structured grids and extended to unstructured grids.     

Accurate local upscaling with variable compact multipoint transmissibility calculations

We propose a new single-phase local upscaling method that uses spatially varying multipoint transmissibility calculations. The method is demonstrated on two-dimensional Cartesian and adaptive Cartesian grids. For each cell face in the coarse upscaled grid, we create a local fine grid region surrounding the face on which we solve two generic local flow problems. The multipoint stencils used to calculate the fluxes across coarse grid cell faces involve the six neighboring pressure values. They are required to honor the two generic flow problems. The remaining degrees of freedom are used to maximize compactness and to ensure that the flux approximation is as close as possible to being two-point. The resulting multipoint flux approximations are spatially varying (a subset of the six neighbors is adaptively chosen) and reduce to two-point expressions in cases without full-tensor anisotropy. Numerical tests show that the method significantly improves upscaling accuracy as compared to commonly used local methods and also compares favorably with a local–global upscaling method.     

Estimation of hydraulic permeability considering the micro morphology of rocks of the borehole YAXCOPOIL-1 (Impact crater Chicxulub,Mexico)

Internal surface, formation factor, Nuclear Magnetic Resonance (NMR)-T2 relaxation times and pore radius distributions were measured on representative core samples for the estimation of hydraulic permeability. Permeability is estimated using various versions of the classic Kozeny–Carman-equation (K–C) and a further development of K–C, the fractal PaRiS-model, taking into account the internal surface. In addition to grain and pore size distribution, directly connected to permeability, internal surface reflects the internal structure (“micro morphology”). Lithologies could be grouped with respect to differences in internal surface. Most melt rich impact breccia lithologies exhibit large internal surfaces, while Tertiary post-impact sediments and Cretaceous lithologies in displaced megablocks display smaller internal surfaces. Investigations with scanning electron microscopy confirm the correlation between internal surface and micro morphology. In addition to different versions of K–C, estimations by means of NMR, pore radius distributions and some gas permeability measurements serve for cross-checking and calibration. In general, the different estimations from the independent methods and the measurements are in satisfactory accordance. For Tertiary limestones and Suevites bulk with very high porosities (up to 35%) permeabilites between 10⁻¹⁴ and 10⁻¹⁶ m² are found, whereas in lower Suevite, Cretaceous anhydrites and dolomites, bulk permeabilites are between 10⁻¹⁵ and 10⁻²³ m².     

Permeability development in vesiculating magmas: implications for fragmentation

 Fragmentation, or the "coming apart" of magma during a plinian eruption, remains one of the least understood processes in volcanology, although assumptions about the timing and mechanisms of fragmentation are key parameters in all existing eruption models. Despite evidence to the contrary, most models assume that fragmentation occurs at a critical vesicularity (volume percent vesicles) of 75–83%. We propose instead that the degree to which magma is fragmented is determined by factors controlling bubble coalescence: magma viscosity, temperature, bubble size distribution, bubble shapes, and time. Bubble coalescence in vesiculating magmas creates permeability which serves to connect the dispersed gas phase. When sufficiently developed, permeability allows subsequent exsolved and expanded gas to escape, thus preserving a sufficiently interconnected region of vesicular magma as a pumice clast, rather than fully fragmenting it to ash. For this reason pumice is likely to preserve information about (a) how permeability develops and (b) the critical permeability needed to insure clast preservation. We present measurements and calculations that constrain the conditions (vesicularity, bubble size distribution, time, pressure difference, viscosity) necessary for adequate permeability to develop. We suggest that magma fragments explosively to ash when and where, in a heterogeneously vesiculating magma, these conditions are not met. Both the development of permeability by bubble wall thinning and rupture and the loss of gas through a permeable network of bubbles require time, consistent with the observation that degree of fragmentation (i.e., amount of ash) increases with increasing eruption rate. Received: 5 July 1995 / Accepted: 27 December 1995     

Surface renewal analysis for sensible and latent heat flux density

High frequency temperature measurements were recorded at five heights and surface renewal (SR) analysis was used to estimate sensible heat flux density (H) over 0.1 m tall grass. Traces of the temperature data showed ramp-like structures, and the mean amplitude and duration of these ramps were used to calculate H using structure functions. Data were compared with H values measured with a sonic anemometer. Latent heat flux density (E) was calculated using an energy balance and the results were compared with E computed from the sonic anemometer data. SR analysis provided good estimates of H for data recorded at all heights but the canopy top and at the highest measurement level, which was above the fully adjusted boundary layer.     

饱和软土压缩曲线的“反常”现象

饱和软土压缩试验时,经常出现压缩曲线“反常”、压缩系数“倒大”的现象。它是土的原始结构发生破坏前后的不同压缩性的客观反映。饱和软土的结构力很微弱,唯质量好的原状土才能见到这种“反常”。应从成孔、取样、测试及资料整理等多个环节保证其工程意义。     

Fault zone geometry of a mature active normal fault: A potential high permeability channel (Pirgaki fault, Corinth rift, Greece)

We present results from petrophysical analysis of a normal fault zone with the aim of defining the flow pathways and their behavior during seismic and interseismic periods. Data are obtained on porosity geometry, strain structure and mineralogy of different domains of a normal fault zone in the Corinth rift. Data point out a close relationship between mineralogy of the clayey minerals, porosity network and strain structures and allow definition of a macroscopic anisotropy of the flow parameters with a strong control by microscopic ultracataclasite structures. The Pirgaki fault zone, developed within pelagic limestone, has a sharp asymmetric porosity profile, with a high porosity volume in the fault core and in the damage zone of the hanging wall. From porosity volumes and threshold measurements, a matrix permeability variation of 6 orders of magnitude could be expected between the protolith and the fault core. Modifications of this pathway during seismic and interseismic phases are depicted. Healing of cracks formed during seismic slip events occurred in the fault core zone and the porous network in the damage zone is sealed in a second step. The lens geometry of the fault core zone is associated with dissolution surfaces and open conduits where dissolved matter could move out of the fault core zone. These elementary processes are developed in particular along Riedel's structures and depend on the orientation of the strain surfaces relative to the local stress and depend also on the roughness of each surface type. P-surfaces are smooth and control shearing process. R-surfaces are rough and present two wavelengths of roughness. The long one controls localization of dissolution surfaces and conduits; the short one is characteristic of dissolution surfaces. The dissolved matter can precipitate in the open structures of the hanging wall damage zone, decreasing the connectivity of the macroscopic conduit developed within this part of the fault zone.     

Early Jurassic rift structures associated with the Soapaga and Boyacá faults of the Eastern Cordillera, Colombia: Sedimentological inferences and regional implications

The NW-trending Bucaramanga fault links, at its southern termination, with the Soapaga and Boyacá faults, which by their NW trend define an ample horsetail structure. As a result of their Neogene reactivation as reverse faults, they bound fault-related anticlines that expose the sedimentary fill of two Early Jurassic rift basins. These sediments exhibit the wedge-like geometry of rift fills related to west-facing normal faults. Their structural setting was controlled further by segmentation of the bounding faults at approximately 10 km intervals, in which each segment is separated by a transverse basement high. Isopach contours and different facies associations suggest these transverse anticlines may have separated depocenters of their adjacent subbasins, which were shaped by a slightly different subsidence history and thereby decoupled. The basin fill of the relatively narrow basin associated with the Soapaga fault is dominated by fanglomeratic successions organized in two coarsening-upward cycles. In the larger basin linked to the Boyacá fault, the sedimentary fill consists of two coarsening-upward sequences that, when fully developed, vary from floodplain to alluvial fan deposits. These Early Jurassic rift fills temporally constrain the evolution of the Bucaramanga fault, which accommodated right-lateral displacement during the early Mesozoic rift event.     

Radiomagnetotelluric mapping, groundwater numerical modelling and 18-Oxygen isotopic data as combined tools to determine the hydrogeological system of a landslide prone area

Three methods were combined to determine the groundwater recharge and transfer processes of a landslide prone area. First, the radiomagnetotelluric method was used to investigate the distribution of electrical resistivity (ρ) of the subsurface and build a three-dimensional model of permeability (k), through an experimental relation between ρ and k. Second, this structural model of permeability and additional climatologic data were used to fix boundary and recharge conditions to perform a three-dimensional and transient numerical simulation of the groundwater flow. Finally 18-Oxygen time series observed at the main springs were used to validate the model. This association of methods led to an improved characterization of the groundwater flow system at local scale and a better understanding of the role of this system on the landslide phenomenon. This structured approach is thought to be useful to design specific remediation strategies to drain the unstable mass.