首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
We develop a semi‐empirical model which combines the theoretical model of Xu and White and the empirical formula of Han, Nur and Morgan in sand–clay environments. This new model may be used for petrophysical interpretation of P‐ and S‐wave velocities. In particular, we are able to obtain an independent estimation of aspect ratios based on log data and seismic velocity, and also the relationship between velocities and other reservoir parameters (e.g. porosity and clay content), thus providing a prediction of shear‐wave velocity. To achieve this, we first use Kuster and Toksöz's theory to derive bulk and shear moduli in a sand–clay mixture. Secondly, Xu and White's model is combined with an artificial neural network to invert the depth‐dependent variation of pore aspect ratios. Finally these aspect ratio results are linked to the empirical formula of Han, Nur and Morgan, using a multiple regression algorithm for petrophysical interpretation. Tests on field data from a North Sea reservoir show that this semi‐empirical model provides simple but satisfactory results for the prediction of shear‐wave velocities and the estimation of reservoir parameters.  相似文献   

2.
We design a velocity–porosity model for sand-shale environments with the emphasis on its application to petrophysical interpretation of compressional and shear velocities. In order to achieve this objective, we extend the velocity–porosity model proposed by Krief et al., to account for the effect of clay content in sandstones, using the published laboratory experiments on rocks and well log data in a wide range of porosities and clay contents. The model of Krief et al. works well for clean compacted rocks. It assumes that compressional and shear velocities in a porous fluid-saturated rock obey Gassmann formulae with the Biot compliance coefficient. In order to use this model for clay-rich rocks, we assume that the bulk and shear moduli of the grain material, and the dependence of the compliance on porosity, are functions of the clay content. Statistical analysis of published laboratory data shows that the moduli of the matrix grain material are best defined by low Hashin–Shtrikman bounds. The parameters of the model include the bulk and shear moduli of the sand and clay mineral components as well as coefficients which define the dependence of the bulk and shear compliance on porosity and clay content. The constants of the model are determined by a multivariate non-linear regression fit for P- and S-velocities as functions of porosity and clay content using the data acquired in the area of interest. In order to demonstrate the potential application of the proposed model to petrophysical interpretation, we design an inversion procedure, which allows us to estimate porosity, saturation and/or clay content from compressional and shear velocities. Testing of the model on laboratory data and a set of well logs from Carnarvon Basin, Australia, shows good agreement between predictions and measurements. This simple velocity-porosity-clay semi-empirical model could be used for more reliable petrophysical interpretation of compressional and shear velocities obtained from well logs or surface seismic data.  相似文献   

3.
The clay-sand mixture model of Xu and White is shown to simulate observed relationships between S-wave velocity (or transit time), porosity and clay content. In general, neither S-wave velocity nor S-wave transit time is a linear function of porosity and clay content. For practical purposes, clay content is approximated by shale volume in well-log applications. In principle, the model can predict S-wave velocity from lithology and any pair of P-wave velocity, porosity and shale volume. Although the predictions should be the same if all measurements are error free, comparison of predictions with laboratory and logging measurements show that predictions using P-wave velocity are the most reliable. The robust relationship between S- and P-wave velocities is due to the fact that both are similarly affected by porosity, clay content and lithology. Moreover, errors in the measured P-wave velocity are normally smaller than those in porosity and shale volume, both of which are subject to errors introduced by imperfect models and imperfect parameters when estimated from logs. Because the model evaluates the bulk and shear moduli of the dry rock frame by a combination of Kuster and Toksöz’ theory and differential effective medium theory, using pore aspect ratios to characterize the compliances of the sand and clay components, the relationship between P- and S-wave velocities is explicit and consistent. Consequently the model sidesteps problems and assumptions that arise from the lack of knowledge of these moduli when applying Gassmann's theory to this relationship, making it a very flexible tool for investigating how the vP-vs relationship is affected by lithology, porosity, clay content and water saturation. Numerical results from the model are confirmed by laboratory and logging data and demonstrate, for example, how the presence of gas has a more pronounced effect on P-wave velocity in shaly sands than in less compliant cleaner sandstones.  相似文献   

4.
This study devises a new analytical relationship to determine the porosity of water-saturated soils at shallow depth using seismic compressional and shear wave velocities. Seismic refraction surveys together with soil sample collection were performed in selected areas containing water-saturated clay–silt, sand and gravely soils. Classification of clay–silt, sand and gravel dense soils provided the coefficient of experimental equation between the data sets, namely, Poisson's ratio, shear modulus and porosity values. This study presents a new analytical relationship between Poisson's ratio and shear modulus values, which are obtained from seismic velocities and porosity values of water-saturated material computed from water content and grain densities, which are determined by laboratory analysis of disturbed samples. The analytical relationship between data sets indicates that when the shear modulus of water-saturated loose soil increases, porosity decreases logarithmically. If shear modulus increases in dense or solid saturated soils, porosity decreases linearly.  相似文献   

5.
—The feasibility of modeling elastic properties of a fluid-saturated sand-clay mixture rock is analyzed by assuming that the rock is composed of macroscopic regions of sand and clay. The elastic properties of such a composite rock are computed using two alternative schemes.¶The first scheme, which we call the composite Gassmann (CG) scheme, uses Gassmann equations to compute elastic moduli of the saturated sand and clay from their respective dry moduli. The effective elastic moduli of the fluid-saturated composite rock are then computed by applying one of the mixing laws commonly used to estimate elastic properties of composite materials.¶In the second scheme which we call the Berryman-Milton scheme, the elastic moduli of the dry composite rock matrix are computed from the moduli of dry sand and clay matrices using the same composite mixing law used in the first scheme. Next, the saturated composite rock moduli are computed using the equations of Brown and Korringa, which, together with the expressions for the coefficients derived by Berryman and Milton, provide an extension of Gassmann equations to rocks with a heterogeneous solid matrix.¶For both schemes, the moduli of the dry homogeneous sand and clay matrices are assumed to obey the Krief’s velocity-porosity relationship. As a mixing law we use the self-consistent coherent potential approximation proposed by Berryman.¶The calculated dependence of compressional and shear velocities on porosity and clay content for a given set of parameters using the two schemes depends on the distribution of total porosity between the sand and clay regions. If the distribution of total porosity between sand and clay is relatively uniform, the predictions of the two schemes in the porosity range up to 0.3 are very similar to each other. For higher porosities and medium-to-large clay content the elastic moduli predicted by CG scheme are significantly higher than those predicted by the BM scheme.¶This difference is explained by the fact that the BM model predicts the fully relaxed moduli, wherein the fluid can move freely between sand and clay regions. In contrast, the CG scheme predicts the no-flow or unrelaxed moduli. Our analysis reveals that due to the extremely low permeability of clays, at seismic and higher frequencies the fluid has no time to move between sand and clay regions. Thus, the CG scheme is more appropriate for clay-rich rocks.  相似文献   

6.
The effect of clay distribution on the elastic properties of sandstones   总被引：1，自引：0，他引：1  
The shape and location of clay within sandstones have a large impact on the P‐wave and S‐wave velocities of the rock. They also have a large effect on reservoir properties and the interpretation of those properties from seismic data and well logs. Numerical models of different distributions of clay – structural, laminar and dispersed clay – can lead to an understanding of these effects. Clay which is located between quartz grains, structural clay, will reduce the P‐wave and S‐wave velocities of the rock. If the clay particles become aligned or form layers, the velocities perpendicular to the alignment will be reduced further. S‐wave velocities decrease more rapidly than P‐wave velocities with increasing clay content, and therefore Poisson's ratios will increase as the velocities decrease. These effects are more pronounced for compacted sandstones. Small amounts of clay that are located in the pore space will have little effect on the P‐wave velocity due to the competing influence of the density effect and pore‐fluid stiffening. The S‐wave velocity will decrease due to the density effect and thus the Poisson's ratio will increase. When there is sufficient clay to bridge the gaps between the quartz grains, P‐wave and S‐wave velocities rise rapidly and the Poisson's ratios decrease. These effects are more pronounced for under‐compacted sandstones. These general results are only slightly modified when the intrinsic anisotropy of the clay material is taken into account. Numerical models indicate that there is a strong, nearly linear relationship between P‐wave and S‐wave velocity which is almost independent of clay distribution. S‐wave velocities can be predicted reasonably accurately from P‐wave velocities based on empirical relationships. However, this does not provide any connection between the elastic and petrophysical properties of the rocks. Numerical modelling offers this connection but requires the inclusion of clay distribution and anisotropy to provide a model that is consistent with both the elastic and petrophysical properties. If clay distribution is ignored, predicting porosities from P‐wave or S‐wave data, for example, can result in large errors. Estimation of the clay distribution from P‐wave and S‐wave velocities requires good estimates of the porosity and clay volume and verification from petrographic analyses of core or cuttings. For a real data example, numerical models of the elastic properties suggest the predominance of dispersed clay in a fluvial sand from matching P‐wave and S‐wave velocity well log data using log‐based estimates of the clay volume and porosity. This is consistent with an interpretation of other log data.  相似文献   

7.
We obtain the wave velocities and quality factors of clay‐bearing sandstones as a function of pore pressure, frequency and partial saturation. The model is based on a Biot‐type three‐phase theory that considers the coexistence of two solids (sand grains and clay particles) and a fluid mixture. Additional attenuation is described with the constant‐Q model and viscodynamic functions to model the high‐frequency behaviour. We apply a uniform gas/fluid mixing law that satisfies the Wood and Voigt averages at low and high frequencies, respectively. Pressure effects are accounted for by using an effective stress law. By fitting a permeability model of the Kozeny– Carman type to core data, the model is able to predict wave velocity and attenuation from seismic to ultrasonic frequencies, including the effects of partial saturation. Testing of the model with laboratory data shows good agreement between predictions and measurements.  相似文献   

8.
地震波本征衰减反映了地层及其所含流体的一些特性,对油气勘探开发有重要意义.已有的理论研究与实验发现,地震频带内的衰减主要与中观尺度(波长与颗粒尺度之间)的斑状部分饱和、完全饱和岩石弹性非均匀性情况下波诱导的局部流体流有关.这种衰减与岩石骨架、孔隙度及充填流体的性质密切相关.本文着重讨论均匀流体分布、斑状或非均匀流体分布两种情况下部分饱和岩石的纵波模量差异.以经典岩石物理理论和衰减机制认识为基础,通过分析低频松弛状态、高频非松弛状态岩石的弹性模量,讨论储层参数(如孔隙度、泥质含量以及含水饱和度等)与纵波衰减之间的确定性关系.上述方法与模型在陆相砂泥岩地层与海相碳酸盐岩地层中的适用性通过常规测井资料得到了初步验证.  相似文献   

9.
We measured in the laboratory ultrasonic compressional and shear‐wave velocity and attenuation (0.7–1.0 MHz) and low‐frequency (2 Hz) electrical resistivity on 63 sandstone samples with a wide range of petrophysical properties to study the influence of reservoir porosity, permeability and clay content on the joint elastic‐electrical properties of reservoir sandstones. P‐ and S‐wave velocities were found to be linearly correlated with apparent electrical formation factor on a semi‐logarithmic scale for both clean and clay‐rich sandstones; P‐ and S‐wave attenuations showed a bell‐shaped correlation (partial for S‐waves) with apparent electrical formation factor. The joint elastic‐electrical properties provide a way to discriminate between sandstones with similar porosities but with different clay contents. The laboratory results can be used to estimate sandstone reservoir permeability from seismic velocity and apparent formation factor obtained from co‐located seismic and controlled source electromagnetic surveys.  相似文献   

10.
Wave velocities in haloanhydrites are difficult to determine and significantly depend on the mineralogy. We used petrophysical parameters to study the wave velocity in haloanhydrites in the Amur Darya Basin and constructed a template of the relation between haloanhydrite mineralogy (anhydrite, salt, mudstone, and pore water) and wave velocities. We used the relation between the P-wave moduli ratio and porosity as constraint and constructed a graphical model (petrophysical template) for the relation between wave velocity, mineral content and porosity. We tested the graphical model using rock core and well logging data.  相似文献   

11.
We have studied 56 unfractured chalk samples of the Upper Cretaceous Tor Formation of the Dan, South Arne and Gorm Fields, Danish North Sea. The samples have porosities of between 14% and 45% and calcite content of over 95%. The ultrasonic compressional‐ and shear‐wave velocities (VP and VS) for dry and water‐saturated samples were measured at up to 75 bar confining hydrostatic pressure corresponding to effective stress in the reservoir. The porosity is the main control of the ultrasonic velocities and therefore of the elastic moduli. The elastic moduli are slightly higher for samples from the South Arne Field than from the Dan Field for identical porosities. This difference may be due to textural differences between the chalk at the two locations because we observe that large grains (i.e. filled microfossils and fossil fragments) that occur more frequently in samples from the Dan Field have a porosity‐reducing effect and that samples rich in large grains have a relatively low porosity for a given P‐wave modulus. The clay content in the samples is low and is mainly represented by either kaolinite or smectite; samples with smectite have a lower P‐wave modulus than samples with kaolinite at equal porosity. We find that ultrasonic VP and VS of dry chalk samples can be satisfactorily estimated with Gassmann's relationships from data for water‐saturated samples. A pronounced difference between the VP/VS ratios for dry and water‐saturated chalk samples indicates promising results for seismic amplitude‐versus‐offset analyses.  相似文献   

12.
We study the propagation of elastic waves that are generated in a fluid‐filled borehole surrounded by a cracked transversely isotropic medium. In the model studied the anisotropy and borehole axes coincide. To obtain the effective elastic moduli of a cracked medium we have applied Hudson's theory that enables the determination of the overall properties as a function of the crack orientation in relation to the symmetry axis of the anisotropic medium. This theory takes into account the hydrodynamic mechanism of the elastic‐wave attenuation caused by fluid filtration from the cracks into a porous matrix. We have simulated the full waveforms generated by an impulse source of finite length placed on the borehole axis. The kinematic and dynamic parameters of the compressional, shear and Stoneley waves as functions of the matrix permeability, crack orientation and porosity were studied. The modelling results demonstrated the influence of the crack‐system parameters (orientation and porosity) on the velocities and amplitudes of all wave types. The horizontally orientated cracks result in maximal decrease of the elastic‐wave parameters (velocities and amplitudes). Based on the fact that the shear‐ and Stoneley‐wave velocities in a transversely isotropic medium are determined by different shear moduli, we demonstrate the feasibility of the acoustic log to identify formations with close to horizontal crack orientations.  相似文献   

13.
In impure chalk, the elastic moduli are not only controlled by porosity but also by contact‐cementation, resulting in relatively large moduli for a given porosity, and by admixtures of clay and fine silica, which results in relatively small moduli for a given porosity. Based on a concept of solids suspended in pore fluids as well as composing the rock frame, we model P‐wave and S‐wave moduli of dry and wet plug samples by an effective‐medium Hashin–Shtrikman model, using chemical, mineralogical and textural input. For a given porosity, the elastic moduli correspond to a part of the solid (the iso‐frame value) forming the frame of an Upper Hashin–Shtrikman bound, whereas the remaining solid is modelled as suspended in the pore fluid. The iso‐frame model is thus a measure of the pore‐stiffness or degree of cementation of the chalk. The textural and mineralogical data may be assessed from logging data on spectral gamma radiation, density, sonic velocity and water saturation in a hydrocarbon zone, whereas the iso‐frame value of a chalk may be assessed from the density and acoustic P‐wave logs alone. The iso‐frame concept may thus be directly used in conventional log‐analysis and is a way of incorporating sonic‐logging data. The Rigs‐1 and Rigs‐2 wells in the South Arne field penetrate the chalk at the same depth but differ in porosity and in water saturation although almost the entire chalk interval has irreducible water saturation. Our model, combined with petrographic data, indicates that the difference in porosity is caused by a higher degree of pore‐filling cementation in Rigs‐1. Petrographic data indicate that the difference in water saturation is caused by a higher content of smectite in the pores of Rigs‐1. In both wells, we find submicron‐size diagenetic quartz.  相似文献   

14.
To investigate the effect of permeability on the propagation of seismo-acoustic waves through marine sediments, a theoretical model based on Biot's equations is established which relates the compressional wave velocity measured at a fixed frequency to computed velocities at zero and infinite frequencies in terms of sediment porosity and permeability. The model is examined experimentally in a standard soil mechanics consolidation test (itself dependent, among other things, on sediment porosity and permeability) which has been modified to include measurements of compressional wave velocity at 1 MHz and shear-wave velocity at 5 kHz. This test allows the elastic modulus of the sediment frame to be assessed under different load conditions simultaneous with the velocity determinations. From a number of tests on different samples, five samples are chosen to typify the range of sediment sizes. The results show that the difference between the measured velocity at 1 MHz and the model-derived velocity at zero frequency increases with increasing particle size (from clays to fine sand), with decreasing porosity, and with increasing permeability. For sediments coarser than fine sand the simple model breaks down, possibly because of the dominance of scattering/diffraction effects at the high frequency of the experiment. Within this limitation the model seems satisfactory to offer a capability of predicting the permeability of a sea floor sediment to an order of magnitude by the in situ measurement of seismic velocities over a wide range of frequencies; the prediction process requires a good in situ determination of sediment porosity such as that offered by electrical formation factor measurements.  相似文献   

15.
成岩作用是影响砂岩声波速度的地质因素之一,定量表征压实和胶结作用的砂岩声波速度岩石物理模型具有重要的理论和实践应用意义.选取视压实率和视胶结率定量表征砂岩成岩作用,通过建立视压实率与颗粒配位数的关系将压实作用的影响引入修正的定量表征胶结作用的CCT模型,最终建立了一种能够定量表征压实和胶结作用对砂岩声波速度影响的岩石物理模型.理论考察发现,随胶结率的增大,岩石声波速度首先迅速增大,随后趋于稳定;随视压实率增大,岩石声波速度同样逐渐增大,当胶结率较大时声波速度变化更为明显.为了验证该声波速度模型,分别对人造砂岩和天然样品进行了声波速度实验观测,结果表明：实验结果与理论分析的趋势吻合良好.该模型易于使用,能够为应用地震和测井资料识别有利储层、定量评价孔隙度以及开展横波速度预测等应用提供理论基础.  相似文献   

16.
We calculate the compressional- and shear-wave velocities of permafrost as a function of unfrozen water content and temperature. Unlike previous theories based on simple slowness and/or moduli averaging or two-phase models, we use a Biot-type three-phase theory that considers the existence of two solids (solid and ice matrices) and a liquid (unfrozen water). The compressional velocity for unconsolidated sediments obtained with this theory is close to the velocity computed with Wood's model, since Biot's theory involves a Wood averaging of the moduli of the single constituents. Moreover, the model gives lower velocities than the well-known slowness averaging theory (Wyllie's equation). For consolidated Berea sandstone, the theory underestimates the value of the compressional velocity below 0°C. Computing the average bulk moduli by slowness averaging the ice and solid phases and Wood averaging the intermediate moduli with the liquid phase yields a fairly good fit of the experimental data. The proportion of unfrozen water and temperature are closely related. Fitting the wave velocity at a given temperature allows the prediction of the velocity at the whole range of temperatures, provided that the average pore radius and its standard deviation are known.  相似文献   

17.
本文定义了各向异性黏弹性参数修正因子,并将其引入到黏弹性模型中以体现泥质含量对黏弹性机制的影响,同时将波传播过程中孔隙介质骨架黏弹性力学机制与两种孔隙流体流动力学机制（Biot流动和喷射流动机制）有机地统一起来处理,从而给出了描述含泥质低孔渗孔隙各向异性介质中波传播规律的黏弹性Biot/squirt (BISQ)模型.数值计算结果表明,入射波的方位角、各向异性渗透率以及泥质含量等对含流体复杂孔隙介质中波频散和衰减的影响具有显著的方位各向异性特征,在低频范围内（地震波勘探频率）黏弹性力学机制对波传播能量的衰减起主导作用.  相似文献   

18.
Continuous, single-channel reflection profiling has been carried out in PVC-lined boreholes, primarily with the aim of ascertaining the position of an old subsurface gas storage tunnel on a proposed dam site. Tube wave reflection patterns thus generated have been interpreted in terms of sediment rigidity and shear wave velocity, and these results could be compared with some independent data. It is interesting to note that, within the well section penetrating Tertiary clays, the velocity of the hydraulic transients apparently was not affected by the PVC casing, which might be explained by a tight coupling between casing and clay wall. In such situations, tube waves turn out a straightforward tool for the determination of shear wave velocity and the derivation of dynamic elastic moduli of unconsolidated sediments. Further applications of the study of the distribution of seismic velocities on the dam site dealt with the consolidation history of the clays. A level of abnormally low P-wave velocities has been detected and interpreted as a gas-charged horizon which, by its coincidence with the base level of clay diapirs, might be considered to have contributed to clay flowage in past geological times. Data about maximum past burial depth, derived from shear wave velocities, turned out to be in agreement with results from consolidation testing.  相似文献   

19.
This paper tests the ability of various rock physics models to predict seismic velocities in shallow unconsolidated sands by comparing the estimates to P and S sonic logs collected in a shallow sand layer and ultrasonic laboratory data of an unconsolidated sand sample. The model fits are also evaluated with respect to the conventional model for unconsolidated sand. Our main approach is to use Hertz‐Mindlin and Walton contact theories, assuming different weight fractions of smooth and rough contact behaviours, to predict the elastic properties of the high porosity point. Using either the Hertz‐Mindlin or Walton theories with rough contact behaviour to define the high porosity endpoint gives an over‐prediction of the velocities. The P‐velocity is overpredicted by a factor of ～1.5 and the S‐velocity by a factor of ～1.8 for highly porous gas‐sand. The degree of misprediction decreases with increasing water saturation and porosity.Using the Hertz‐Mindlin theory with smooth contact behaviour or weighted Walton models gives a better fit to the data, although the data are best described using the Walton smooth model. To predict the properties at the lower porosities, the choice of bounding model attached to the Walton Smooth model controls the degree of fit to the data, where the Reuss bound best captures the porosity variations of dry and wet sands in this case since they are caused by depositional differences. The empirical models based on lab experiments on unconsolidated sand also fit the velocity data measured by sonic logs in situ, which gives improved confidence in using lab‐derived results.  相似文献   

20.
To simultaneously take into account the Biot-flow mechanism, the squirt-flow mechanism, and the frame-viscoelasticity mechanism, a generalized viscoelastic BISQ （Biot/squirt） model is developed for wave propagation in clay-bearing sandstones based on Dvorkin＇s elastic BISQ model. The present model is extended to a wide range of permeabilities （k 〉 0.05 mD） by introducing a dimensionless correction factor for viscoelastic parameters, defined as a function of the permeability and the clay content. We describe the frame＇s stress-strain relationship of the clay-bearing sandstones by the differential constitutive equations of generalized viscoelasticity and then derive the viscoelastic-wave dynamic equations. With the assumption of a plane-wave solution, we finally yield the phase velocities and the attenuation coefficients by solving the dynamic wave equations in the frequency and wave number domain. The comparison of numerical results and experimental data shows that the generalized viscoelastic BISQ model is applicable for modeling the wave propagation in most of the sandstones mainly bearing kaolinite clay.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号