致密砂岩储层变岩电参数法饱和度计算模型:以苏里格气田西区盒8段例

致密砂岩储层是油气勘探的重点方向之一,由于致密砂岩储层具有低孔低渗、孔隙结构复杂、非均质性强和岩电参数变化等特点,在测井解释过程中不同于常规储层。Archie公式在致密砂岩地层中使用时,强非均质性导致岩电参数a、b、m和n会随着储层的变化而变化。针对苏里格西区盒₈段致密砂岩储层提出利用变岩电参数建立饱和度模型。依据研究区盒₈段储层25块岩心岩电实验分析发现胶结系数m与lgϕ呈二次函数关系、m与岩性系数a呈幂函数关系、束缚水饱和度S_wb对饱和度指数n有较大的影响。因此,在对孔隙度、束缚水饱和度和孔隙结构指数等因素分析的基础上,综合各因素影响特征建立了变岩电参数的饱和度模型。结果表明,该方法相对传统Archie公式饱和度解释,模型精度更高,为解释复杂孔隙结构砂岩储层提供了一种有效预测方法。     

非饱和土的本构模型研究

引用平均土骨架应力的概念,研究推导出非饱和土的刚度参数随吸力变化而变化的关系式,进而推导得到用平均土骨架应力表述的非饱和土LC屈服面函数以及硬化规律。从土力学原理推导,得到土样由于在净应力和吸力作用下产生体积变形引起土样饱和度变化的关系式。由平均土骨架应力推广,得到三轴应力状态的椭圆屈服函数,这一非饱和土本构模型的优点在于考虑了应力作用后土样饱和度的变化,通过对已有试验数据的初步验证,表明提出的非饱和土本构模型的合理性和适用性。     

随钻热解地球化学录井含油饱和度的计算和应用

热解地化录井可以检测储集层含油气总量,推算原油密度值,结合储集层孔隙度值和储层岩石密度,经过换算可得出含油饱和度。热解地化录井含油饱和度是基于热解地化参数而获得,获取的机制是直接检测与合理的理论恢复相结合,能够客观的反映储层中所含油气的数量,评价储集层的性质,经试油结论验证,随钻评价结论与实际情况相符良好。     

Three-phase measurements of oil and gas trapping in sand packs

We measure the trapped saturations of oil and gas as a function of initial saturation in water-wet sand packs. We start with a water-saturated column and inject octane (oil), while water and oil are produced from the bottom. Once water production has ceased, air (gas) then enters from the top, allowing oil and gas to drain under gravity for different times. Finally water is then injected from the bottom to trap both oil and gas. The columns are sliced and the fluids analyzed using gas chromatography. We find that for high initial gas saturations more gas can be trapped in the presence of oil than in a two-phase (gas/water) system. The residual gas saturation can be over 20% compared to 14% in two-phase flow [Al Mansoori SK, Iglauer S, Pentland CH, Bijeljic B, Blunt MJ. Measurements of non-wetting phase trapping applied to carbon dioxide storage. Energy Procedia 2009;1(1):3173–80]. This is unlike previous measurements on consolidated media, where the trapped gas saturation is either similar or lower to that reached in an equivalent two-phase experiment. For lower initial gas saturation, the amount of trapping follows the initial-residual trend seen in two-phase experiments. The amount of oil trapped is insensitive to initial gas saturation or the amount of gas that is trapped, again in contrast to measurements on consolidated media. More oil is trapped than would be predicted from an equivalent two-phase (oil/water) system, although the trapped saturation is never larger than the maximum reached in two-phase flow (around 11%) [Pentland CH, Al Mansoori SK, Iglauer S, Bijeljic B, Blunt MJ. Measurement of non-wetting phase trapping in sand packs. In: SPE 115697, proceedings of the SPE annual technical conference and exhibition, Denver, Colorado, USA; 21–24 September 2008]. These initially surprising results are explained in the context of oil layer stability and the competition between snap-off and piston-like advance. In two-phase systems, displacement is principally by cooperative piston-like advance with relatively little trapping, whereas in consolidated media snap-off is generally more significant. However, oil layer collapse events during three-phase waterflooding rapidly trap the oil which acts as a barrier to direct water/gas displacement, except by snap-off, leading to enhanced gas trapping.     

Effects of random variations of soil properties on site amplification of seismic ground motions

This paper investigates the effects of random variations of soil properties on site amplification of seismic waves. First, based on attenuation laws and the filtered Tajimi–Kanai spectrum, seismic motion at the base rock of a soil site is stochastically generated according to an assumed earthquake with a given magnitude and epicentral distance. Motions on the surface of this layered random soil site are calculated by nonlinear wave propagation methods, and by assuming the incoming seismic wave consisting of SH wave or combined P and SV waves. Soil properties, including shear modulus, damping ratio and mass density, as well as ground water level are considered as random in the numerical calculation. The Rosenblueth method is used to solve the random dynamic responses of the soil site. Parametric calculations are performed to investigate the effects of various parameters on site amplification of seismic waves. The mean and maximum ground motions on surface of the site are estimated. Numerical results indicate that the estimated surface motions differ substantially if the random variations of soil properties and soil saturation level are taken into consideration in the analysis.     

Role of water in the origin of podiform chromitite deposits

We report experiments in basalt oversaturated with water to duplicate the nodular ore textures of podiform chromitite ores. In immiscible basalt-water systems saturated with olivine and chromite, olivine will reside in the melt while chromite will collect in the fluid phase. Fractionation is physical and is driven by differential wetting properties of melt and fluid against silicate and oxide surfaces. There is no need to suppress olivine from the liquidus of a primitive basaltic melt as suggested by Irvine [Irvine, T.N., Geology 5 (1977) 273-277], to achieve chromite accumulations as observed in natural podiform ore deposits. The results imply that podiform chromitite ores will form where a primitive olivine-chromite-saturated mantle melt is sufficiently water-rich to exsolve a fluid phase during passage through the uppermost mantle. The most likely geodynamic environment for podiform chromite mineralization to take place is a supra-subduction zone setting.     

Infiltration,runoff and sediment production in blanket peat catchments: implications of field rainfall simulation experiments

Blanket peat covers the headwaters of many major European rivers. Runoff production in upland blanket peat catchments is flashy with large flood peaks and short lag times; there is minimal baseflow. Little is known about the exact processes of infiltration and runoff generation within these upland headwaters. This paper presents results from a set of rainfall simulation experiments performed on the blanket peat moorland of the North Pennines, UK. Rainfall was simulated at low intensities (3–12 mm h^?1), typical of natural rainfall, on bare and vegetated peat surfaces. Runoff response shows that infiltration rate increases with rainfall intensity; the use of low‐intensity rainfall therefore allows a more realistic evaluation of infiltration rates and flow processes than previous studies. Overland flow is shown to be common on both vegetated and bare peat surfaces although surface cover does exert some control. Most runoff is produced within the top few centimetres of the peat and runoff response decreases rapidly with depth. Little vertical percolation takes place to depths greater than 10 cm owing to the saturation of the peat mass. This study provides evidence that the quickflow response of upland blanket peat catchments is a result of saturation‐excess overland flow generation. Rainfall–runoff response from small plots varies with season. Following warm, dry weather, rainfall tends to infiltrate more readily into blanket peat, not just initially but to the extent that steady‐state surface runoff rates are reduced and more flow takes place within the peat, albeit at shallow depth. Sediment erosion from bare peat plots tends to be supply limited. Seasonal weather conditions may affect this in that after a warm, dry spell, surface desiccation allows sediment erosion to become transport limited. Copyright © 2002 John Wiley & Sons, Ltd.     

流体饱和度对裂缝-孔隙岩石频变纵波速度及其各向异性的影响

流体饱和度会改变裂缝性储层的纵波速度,从而影响地层速度的频散特性及各向异性程度,导致储层地震响应特征复杂,储层预测多解性强,流体识别难度大.本文根据多相流体饱和裂缝性储层的特点,借助于Norris和KG模型,建立部分饱和裂缝-孔隙等效介质模型,给出频变地震波速度随流体饱和度变化的精确关系式.数值模拟结果表明,当气、水两相共存时,随着含水饱和度的增加,高频段纵波相速度逐渐增大,各向异性程度逐渐减小;低频段纵波相速度逐渐减小,各向异性程度不变;相速度频散及其各向异性程度逐渐增强.组合已有的孔隙弹性理论模型,对实验室人工裂缝-孔隙砂岩岩样的纵波速度进行拟合,计算得到的曲线与实验室测量散点值吻合度较高,表明组合模型在给定参数下的有效性.该研究能够为多相流体饱和裂缝性储层的地震响应特征分析奠定扎实的理论基础,为提高储层预测的确定性和流体识别的准确性提供可靠的理论依据.

     

Quantifying Damage, Saturation and Anisotropy in Cracked Rocks by Inverting Elastic Wave Velocities

Crack damage results in a decrease of elastic wave velocities and in the development of anisotropy. Using non-interactive crack effective medium theory as a fundamental tool, we calculate dry and wet elastic properties of cracked rocks in terms of a crack density tensor, average crack aspect ratio and mean crack fabric orientation from the solid grains and fluid elastic properties. Using this same tool, we show that both the anisotropy and shear-wave splitting of elastic waves can be derived. Two simple crack distributions are considered for which the predicted anisotropy depends strongly on the saturation, reaching up to 60% in the dry case. Comparison with experimental data on two granites, a basalt and a marble, shows that the range of validity of the non-interactive effective medium theory model extends to a total crack density of approximately 0.5, considering symmetries up to orthorhombic. In the isotropic case, Kachanov's (1994) non-interactive effective medium model was used in order to invert elastic wave velocities and infer both crack density and aspect ratio evolutions. Inversions are stable and give coherent results in terms of crack density and aperture evolution. Crack density variations can be interpreted in terms of crack growth and/or changes of the crack surface contact areas as cracks are being closed or opened respectively. More importantly, the recovered evolution of aspect ratio shows an exponentially decreasing aspect ratio (and therefore aperture) with pressure, which has broader geophysical implications, in particular on fluid flow. The recovered evolution of aspect ratio is also consistent with current mechanical theories of crack closure. In the anisotropic cases—both transverse isotropic and orthorhombic symmetries were considered—anisotropy and saturation patterns were well reproduced by the modelling, and mean crack fabric orientations we recovered are consistent with in situ geophysical imaging. Our results point out that: (1) It is possible to predict damage, anisotropy and saturation in terms of a crack density tensor and mean crack aspect ratio and orientation; (2) using well constrained wave velocity data, it is possible to extrapolate the contemporaneous evolution of crack density, anisotropy and saturation using wave velocity inversion as a tool; 3) using such an inversion tool opens the door in linking elastic properties, variations to permeability.     

Review Paper: Post-stack 4D seismic time-shifts: Interpretation and evaluation

A review and analysis of post-stack time-lapse time-shifts has been carried out that covers published literature supplemented by in-house datasets available to the authors. Time-shift data are classified into those originating from geomechanical effects and those due to fluid saturation changes. From these data, conclusions are drawn regarding the effectiveness of post-stack time-shifts for overburden and reservoir monitoring purposes. A variety of field examples are shown that display the range and magnitude of variation for each class of application. The underlying physical mechanisms creating these time-shifts are then described, and linked to a series of generic and field-specific rock physics calculations that predict their magnitudes. These calculations serve as a guide for practitioners wishing to utilize this information on their own datasets. Conclusions are drawn regarding the reliability of this attribute for monitoring purposes, and the extent to which further development is required and how it should be reported by authors.