储层建模研究进展及发展趋势

油气田开发的后期进入高含水阶段,为了更加经济准确地进行油气开发,有必要采用储层地质建模的方法对老油气田进行储层研究。详细阐述了国内外储层地质建模的发展史,对储层地质建模的方法进行了细致的分类及论述。方法分类包括确定性建模方法和随机性建模方法。每一种建模方法又有多种子方法。提出了目前储层地质建模研究尚未很好解决的一些问题,如建模的对象局限于常规的碎屑岩储层。虽然有学者对火成岩、裂缝碳酸盐岩进行了相关的探究工作,但目前对非常规储层涉及较少,储层建模的精度也有待提高。     

Anisotropic dual-continuum representations for multiscale poroelastic materials: Development and numerical modelling

Dual-continuum (DC) models can be tractable alternatives to explicit approaches for the numerical modelling of multiscale materials with multiphysics behaviours. This work concerns the conceptual and numerical modelling of poroelastically coupled dual-scale materials such as naturally fractured rock. Apart from a few exceptions, previous poroelastic DC models have assumed isotropy of the constituents and the dual-material. Additionally, it is common to assume that only one continuum has intrinsic stiffness properties. Finally, little has been done into validating whether the DC paradigm can capture the global poroelastic behaviours of explicit numerical representations at the DC modelling scale. We address the aforementioned knowledge gaps in two steps. First, we utilise a homogenisation approach based on Levin's theorem to develop a previously derived anisotropic poroelastic constitutive model. Our development incorporates anisotropic intrinsic stiffness properties of both continua. This addition is in analogy to anisotropic fractured rock masses with stiff fractures. Second, we perform numerical modelling to test the DC model against fine-scale explicit equivalents. In doing, we present our hybrid numerical framework, as well as the conditions required for interpretation of the numerical results. The tests themselves progress from materials with isotropic to anisotropic mechanical and flow properties. The fine-scale simulations show that anisotropy can have noticeable effects on deformation and flow behaviour. However, our numerical experiments show that the DC approach can capture the global poroelastic behaviours of both isotropic and anisotropic fine-scale representations.     

The fate of Chernobyl radiocesium in the River Öre catchment,northern Sweden

A¹³⁷Cs-balance for the catchment of the River Öre in central northern Sweden which received about 30 kBq m^–2 of radiocesium from the Chernobyl accident was calculated for the period 1986–1991. Altogether, slightly less than 10% of the total deposition in the catchment was estimated to be exported from the terrestrial parts during this period of time. More than 90% of this loss is transported with the River Öre to the outer sea of the Gulf of Bothnia. The retention in Lake Örträsket which is the only lake along the river course and the Öre Estuary outside the river mouth was thus slightly less than 10%. Nearly all of the radiocesium deposited in the lake is permanently retained in the sediments and successively covered with less radioactive sediment. A considerable export of radiocesium from the estuary to the outer sea takes place due to resuspension and subsequent transport by wind and wave generated currents.     

Numerical modeling for analyzing thermal surface anomalies induced by underground coal fires

Coal seams burning underneath the surface are recognized all over the world and have drawn increasing public attention in the past years. Frequently, such fires are analyzed by detecting anomalies like increased exhaust gas concentrations and soil temperatures at the surface. A proper analysis presumes the understanding of involved processes, which determine the spatial distribution and dynamic behavior of the anomalies.In this paper, we explain the relevance of mechanical and energy transport processes with respect to the occurrence of temperature anomalies at the surface. Two approaches are presented, aiming to obtain insight into the underground coal fire situation: In-situ temperature mapping and numerical simulation. In 2000 to 2005, annual temperature mapping in the Wuda (Inner Mongolia, PR China) coal fire area showed that most thermal anomalies on the surface are closely related to fractures, where hot exhaust gases from the coal fire are released. Those fractures develop due to rock mechanical failure after volume reduction in the seams. The measured signals at the surface are therefore strongly affected by mechanical processes.More insight into causes and effects of involved energy transport processes is obtained by numerical simulation of the dynamic behavior of coal fires. Simulations show the inter-relation between release and transport of thermal energy in and around underground coal fires. Our simulation results show a time delay between the coal fire propagation and the observed appearance of the surface temperature signal. Additionally, the overall energy flux away from the burning coal seam into the surrounding bedrock is about 30-times higher than the flux through the surface. This is of particular importance for an estimation of the energy released based on surface temperature measurements. Finally, the simulation results also prove that a fire propagation rate estimated from the interpretation of surface anomalies can differ from the actual rate in the seam.     

A Holistic Approach to Upwelling and Downwelling along the South-West Coast of India

An attempt has been made to develop a holistic understanding of upwelling and downwelling along the south-west coast of India. The main objective was to elucidate the roles of different forcings involved in the vertical motion along this coast. The south-west coast of India was characterized by upwelling during the south-west monsoon (May to September) and by downwelling during the north-east monsoon and winter (November to February). The average vertical velocity calculated along the south-west coast from the vertical shift of the 26?°C isotherm is 0.57?m/day during upwelling and 0.698?m/day during downwelling. It was concluded that upwelling along the south-west coast of India is driven by offshore Ekman transport due to the alongshore wind, Ekman pumping, horizontal divergence of currents and by the propagation of coastally trapped waves. Whereas downwelling along the coast is driven only by convergence of currents and the propagation of coastally trapped Kelvin waves. Along the west coast of India, the downwelling-favorable Kelvin waves come from the equator and upwelling-favorable waves come from the Gulf of Mannar region.     

Phase equilibrium constraints on a deep crustal metamorphic field gradient: metapelitic rocks from the Ivrea Zone (NW Italy)

The metamorphic rocks of the Ivrea Zone in NW Italy preserve a deep crustal metamorphic field gradient. Application of quantitative phase equilibria methods to metapelitic rocks provides new constraints on the P–T conditions recorded in Val Strona di Omegna, Val Sesia and Val Strona di Postua. In Val Strona di Omegna, the metapelitic rocks show a structural and mineralogical change from mica‐schists with the common assemblage bi–mu–sill–pl–q–ilm ± liq at the lowest grades, through metatexitic migmatites (g–sill–bi–ksp–pl–q–ilm–liq) at intermediate grades, to complex diatexitic migmatites (g–sill–ru–bi–ksp–pl–q–ilm–liq) at the highest grades. Partial melting in the metapelitic rocks is consistent with melting via the breakdown of first muscovite then biotite. The metamorphic field gradient in Val Strona di Omegna is constrained to range from conditions of ～3.5–6.5 kbar at ≈650 °C to ～10–12 kbar at >900 °C. The peak P–T estimates, particularly for granulite facies conditions, are significantly higher than those of most earlier studies. In Val Sesia and Val Strona di Postua, cordierite‐bearing rocks record the effects of contact metamorphism associated with the intrusion of a large mafic body (the Mafic Complex). The contact metamorphism occurred at lower pressures than the regional metamorphic peak and overprints the regional metamorphic assemblages. These relationships are consistent with the intrusion of the Mafic Complex having post dated the regional metamorphism and are inconsistent with a model of magmatic underplating as the cause of granulite facies metamorphism in the region.     

Inverse modelling of aquifer parameters in basaltic rock with the help of pumping test method using MODFLOW software

The origin and movement of groundwater are the fundamental questions that address both the temporal and spatial aspects of ground water run and water supply related issues in hydrological systems. As groundwater flows through an aquifer, its composition and temperature may variation dependent on the aquifer condition through which it flows. Thus, hydrologic investigations can also provide useful information about the subsurface geology of a region. But because such studies investigate processes that follow under the Earth's shallow, obtaining the information necessary to answer these questions is not continuously easy. Springs, which discharge groundwater table directly, afford to study subsurface hydrogeological processes.The present study of estimation of aquifer factors such as transmissivity (T) and storativity (S) are vital for the evaluation of groundwater resources. There are several methods to estimate the accurate aquifer parameters (i.e. hydrograph analysis, pumping test, etc.). In initial days, these parameters are projected either by means of in-situ test or execution test on aquifer well samples carried in the laboratory. The simultaneous information on the hydraulic behavior of the well (borehole) that provides on this method, the reservoir and the reservoir boundaries, are important for efficient aquifer and well data management and analysis. The most common in-situ test is pumping test performed on wells, which involves the measurement of the fall and increase of groundwater level with respect to time. The alteration in groundwater level (drawdown/recovery) is caused due to pumping of water from the well. Theis (1935) was first to propose method to evaluate aquifer parameters from the pumping test on a bore well in a confined aquifer. It is essential to know the transmissivity (T = Kb, where b is the aquifer thickness; pumping flow rate, Q = TW (dh/dl) flow through an aquifer) and storativity (confined aquifer: S = bS_s, unconfined: S = S_y), for the characterization of the aquifer parameters in an unknown area so as to predict the rate of drawdown of the groundwater table/potentiometric surface throughout the pumping test of an aquifer. The determination of aquifer's parameters is an important basis for groundwater resources evaluation, numerical simulation, development and protection as well as scientific management. For determining aquifer's parameters, pumping test is a main method. A case study shows that these techniques have been fast speed and high correctness. The results of parameter's determination are optimized so that it has important applied value for scientific research and geology engineering preparation.     

井水温度微动态及其形成机制

介绍了我国地震地下水温度观测概况及已获得的微动态现象, 即水温固体潮、 水温同震效应、 水温前兆异常等, 并对其形成机理做了理论分析。 提出水温微动态形成机理研究应从二个地温梯度、 二个对流、 一个传导、 一个扩散等热系统与热过程去考虑。     

Cross-shelf transport into nearshore waters due to shoaling internal tides in San Pedro Bay, CA

In the summer of 2001, a coastal ocean measurement program in the southeastern portion of San Pedro Bay, CA, was designed and carried out. One aim of the program was to determine the strength and effectiveness of local cross-shelf transport processes. A particular objective was to assess the ability of semidiurnal internal tidal currents to move suspended material a net distance across the shelf. Hence, a dense array of moorings was deployed across the shelf to monitor the transport patterns associated with fluctuations in currents, temperature and salinity. An associated hydrographic program periodically monitored synoptic changes in the spatial patterns of temperature, salinity, nutrients and bacteria. This set of measurements show that a series of energetic internal tides can, but do not always, transport subthermocline water, dissolved and suspended material from the middle of the shelf into the surfzone. Effective cross-shelf transport occurs only when (1) internal tides at the shelf break are strong and (2) subtidal currents flow strongly downcoast. The subtidal downcoast flow causes isotherms to tilt upward toward the coast, which allows energetic, nonlinear internal tidal currents to carry subthermocline waters into the surfzone. During these events, which may last for several days, the transported water remains in the surfzone until the internal tidal current pulses and/or the downcoast subtidal currents disappear. This nonlinear internal tide cross-shelf transport process was capable of carrying water and the associated suspended or dissolved material from the mid-shelf into the surfzone, but there were no observation of transport from the shelf break into the surfzone. Dissolved nutrients and suspended particulates (such as phytoplankton) transported from the mid-shelf into the nearshore region by nonlinear internal tides may contribute to nearshore algal blooms, including harmful algal blooms that occur off local beaches.