首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 62 毫秒
1.
Predicting the shear‐wave (S‐wave) velocity is important in seismic modelling, amplitude analysis with offset, and other exploration and engineering applications. Under the low‐frequency approximation, the classical Biot–Gassmann theory relates the Biot coefficient to the bulk modulus of water‐saturated sediments. If the Biot coefficient under in situ conditions can be estimated, the shear modulus or the S‐wave velocity can be calculated. The Biot coefficient derived from the compressional‐wave (P‐wave) velocity of water‐saturated sediments often differs from and is less than that estimated from the S‐wave velocity, owing to the interactions between the pore fluid and the grain contacts. By correcting the Biot coefficients derived from P‐wave velocities of water‐saturated sediments measured at various differential pressures, an accurate method of predicting S‐wave velocities is proposed. Numerical results indicate that the predicted S‐wave velocities for consolidated and unconsolidated sediments agree well with measured velocities.  相似文献   

2.
在实验室对5块储层砂岩进行了模拟地层压力条件下的超声波速度测试。砂岩样品采自WXS凹陷的W地层,覆盖了从低到高的孔隙度和渗透率范围。实验选用了卤水和4种不同密度油作为孔隙流体,结合温度变化,实现了对流体粘度引致的速度频散研究。对实验结果的分析表明:(1)对于高孔隙度和渗透率的样品,无论是哪种流体饱和,观察到的超声波速度测试值和零频率Gassmann预测值的差异较小(约2-3%),基本上可以用Biot模型解释;对于中等孔隙度和渗透率的样品,低粘度流体(<约3mP?S)的频散效应也可以用Biot模型得到合理解释;(2)对于低、中孔隙度和渗透率样品,当流体粘度增加时,喷射流机制起主导作用,导致严重的速度频散(可达8%)。对储层砂岩的微裂隙纵横比进行了估计并用于喷射流特征频率的计算,当高于该特征频率时,Gassmann理论的假设条件受到破坏,实验室测得的高频速度不能直接用于地震低频条件下的W地层砂岩的Gassmann流体替换研究。  相似文献   

3.
This study presents the results of experimental compaction while measuring ultrasonic velocities of sands with different grain size, shape, sorting and mineralogy. Uniaxial mechanical compaction tests up to a maximum of 50 MPa effective stress were performed on 29 dry sand aggregates derived from eight different sands to measure the rock properties. A good agreement was found between the Gassmann saturated bulk moduli of dry and brine saturated tests of selected sands. Sand samples with poor sorting showed low initial porosity while sands with high grain angularity had high initial porosity. The sand compaction tests showed that at a given stress well‐sorted, coarse‐grained sands were more compressible and had higher velocities (Vp and Vs) than fine‐grained sands when the mineralogy was similar. This can be attributed to grain crushing, where coarser grains lead to high compressibility and large grain‐to‐grain contact areas result in high velocities. At medium to high stresses the angular coarse to medium grained sands (both sorted sands and un‐sorted whole sands) showed high compaction and velocities (Vp and Vs). The small grain‐to‐grain contact areas promote higher deformation at grain contacts, more crushing and increased porosity loss resulting in high velocities. Compaction and velocities (Vp and Vs) increased with decreasing sorting in sands. However, at the same porosity, the velocities in whole sands were slightly lower than in the well‐sorted sands indicating the presence of loose smaller grains in‐between the framework grains. Quartz‐poor sands (containing less than 55% quartz) showed higher velocities (Vp and Vs) compared to that of quartz‐rich sands. This could be the result of sintering and enlargement of grain contacts of ductile mineral grains in the quartz‐poor sands increasing the effective bulk and shear stiffness. Tests both from wet measurements and Gassmann brine substitution showed a decreasing Vp/Vs ratio with increasing effective stress. The quartz‐rich sands separated out towards the higher side of the Vp/Vs range. The Gassmann brine substituted Vp and Vs plotted against effective stress provide a measure of the expected velocity range to be found in these and similar sands during mechanical compaction. Deviations of actual well log data from experimental data may indicate uplift, the presence of hydrocarbon, overpressure and/or cementation. Data from this study may help to model velocity‐depth trends and to improve the characterization of reservoir sands from well log data in a low temperature (<80–100o C) zone where compaction of sands is mostly mechanical.  相似文献   

4.
Quantifying the effects of pore-filling materials on elastic properties of porous rocks is of considerable interest in geophysical practice. For rocks saturated with fluids, the Gassmann equation is proved effective in estimating the exact change in seismic velocity or rock moduli upon the changes in properties of pore infill. For solid substance or viscoelastic materials, however, the Gassmann theory is not applicable as the rigidity of the pore fill (either elastic or viscoelastic) prevents pressure communication in the pore space, which is a key assumption of the Gassmann equation. In this paper, we explored the elastic properties of a sandstone sample saturated with fluid and solid substance under different confining pressures. This sandstone sample is saturated with octadecane, which is a hydrocarbon with a melting point of 28°C, making it convenient to use in the lab in both solid and fluid forms. Ultrasonically measured velocities of the dry rock exhibit strong pressure dependency, which is largely reduced for the filling of solid octadecane. Predictions by the Gassmann theory for the elastic moduli of the sandstone saturated with liquid octadecane are consistent with ultrasonic measurements, but underestimate the elastic moduli of the sandstone saturated with solid octadecane. Our analysis shows that the difference between the elastic moduli of the dry and solid-octadecane-saturated sandstone is controlled by the squirt flow between stiff, compliant, and the so-called intermediate pores (with an aspect ratio larger than that of compliant pore but much less than that of stiff pores). Therefore, we developed a triple porosity model to quantify the combined squirt flow effects of compliant and intermediate pores saturated with solid or viscoelastic infill. Full saturation of remaining stiff pores with solid or viscoelastic materials is then considered by the lower embedded bound theory. The proposed model gave a reasonable fit to the ultrasonic measurements of the elastic moduli of the sandstone saturated with liquid or solid octadecane. Comparison of the predictions by the new model to other solid substitution schemes implied that accounting for the combined effects of compliant and intermediate pores is necessary to explain the solid squirt effects.  相似文献   

5.
Investigating seismic dispersion and attenuation characteristics of loosely compacted marine sandstone is essential in reconciling different geophysical measurements (surface seismic, well logging and ultrasonic) for better characterization of a shallow marine sandstone reservoir. We have experimented with a typical high-porosity and high-permeability sandstone sample, extracted from the Paleogene marine depositional setting in the Gulf of Mexico, at the low-frequency band (2–500 Hz) as well as ultrasonic point (106 Hz), to investigate the effects of varying saturation levels on a rock's elasticity. The results suggest that the Young's modulus of the measured sample with adsorbed moisture at laboratory conditions (room temperature, 60%–90% humidity) exhibits dispersive behaviours. The extensional attenuation can be as high as 0.038, and the peak frequency occurs around 60 Hz. The extensional attenuation due to moisture adsorption can be dramatically mitigated with the increase of confining pressure. For partial saturation status, extensional attenuation increases as increasing water saturation by 79% with respect to the measured frequencies. Additionally, the results show that extensional attenuation at the fully water-saturated situation is even smaller than that at adsorbed moisture conditions. The Gassmann–Wood model can overall capture the P-wave velocity-saturation trend of measured data at seismic frequencies, demonstrating that the partially saturated unconsolidated sandstone at the measured seismic frequency range is prone to be in the relaxed status. Nevertheless, the ultrasonic velocities are significantly higher than the Gassmann–Wood predictions, suggesting that the rocks are in the unrelaxed status at the ultrasonic frequency range. The poroelastic modelling results based on the patchy saturation model also indicate that the characteristic frequency of the partially saturated sample is likely beyond the measured seismic frequency range.  相似文献   

6.
CO2 saturations are estimated at Sleipner using a two-step imaging workflow. The workflow combines seismic tomography (full-waveform inversion) and rock physics inversion and is applied to a two-dimensional seismic line located near the injection point at Sleipner. We use baseline data (1994 vintage, before CO2 injection) and monitor data that was acquired after 12 years of CO2 injection (2008 vintage). P-wave velocity models are generated using the Full waveform inversion technology and then, we invert selected rock physics parameters using an rock physics inversion methodology. Full waveform inversion provides high-resolution P-wave velocity models both for baseline and monitor data. The physical relations between rock physics properties and acoustic wave velocities in the Utsira unconsolidated sandstone (reservoir formation) are defined using a dynamic rock physics model based on well-known Biot–Gassmann theories. For data prior to injection, rock frame properties (porosity, bulk and shear dry moduli) are estimated using rock physics inversion that allows deriving physically consistent properties with related uncertainty. We show that the uncertainty related to limited input data (only P-wave velocity) is not an issue because the mean values of parameters are correct. These rock frame properties are then used as a priori constraint in the monitor case. For monitor data, the Full waveform inversion results show nicely resolved thin layers of CO2–brine saturated sandstones under intra-reservoir shale layers. The CO2 saturation estimation is carried out by plugging an effective fluid phase in the rock physics model. Calculating the effective fluid bulk modulus of the brine–CO2 mixture (using Brie equation in our study) is shown to be the key factor to link P-wave velocity to CO2 saturation. The inversion tests are done with several values of Brie/patchiness exponent and show that the CO2 saturation estimates are varying between 0.30 and 0.90 depending on the rock physics model and the location in the reservoir. The uncertainty in CO2 saturation estimation is usually lower than 0.20. When the patchiness exponent is considered as unknown, the inversion is less constrained and we end up with values of exponent varying between 5 and 20 and up to 33 in specific reservoir areas. These estimations tend to show that the CO2–brine mixing is between uniform and patchy mixing and variable throughout the reservoir.  相似文献   

7.
Velocities of compressional and shear waves in limestones   总被引：2，自引：1，他引：2  
Carbonate rocks are important hydrocarbon reservoir rocks with complex textures and petrophysical properties (porosity and permeability) mainly resulting from various diagenetic processes (compaction, dissolution, precipitation, cementation, etc.). These complexities make prediction of reservoir characteristics (e.g. porosity and permeability) from their seismic properties very difficult. To explore the relationship between the seismic, petrophysical and geological properties, ultrasonic compressional‐ and shear‐wave velocity measurements were made under a simulated in situ condition of pressure (50 MPa hydrostatic effective pressure) at frequencies of approximately 0.85 MHz and 0.7 MHz, respectively, using a pulse‐echo method. The measurements were made both in vacuum‐dry and fully saturated conditions in oolitic limestones of the Great Oolite Formation of southern England. Some of the rocks were fully saturated with oil. The acoustic measurements were supplemented by porosity and permeability measurements, petrological and pore geometry studies of resin‐impregnated polished thin sections, X‐ray diffraction analyses and scanning electron microscope studies to investigate submicroscopic textures and micropores. It is shown that the compressional‐ and shear‐wave velocities (Vp and Vs, respectively) decrease with increasing porosity and that Vp decreases approximately twice as fast as Vs. The systematic differences in pore structures (e.g. the aspect ratio) of the limestones produce large residuals in the velocity versus porosity relationship. It is demonstrated that the velocity versus porosity relationship can be improved by removing the pore‐structure‐dependent variations from the residuals. The introduction of water into the pore space decreases the shear moduli of the rocks by about 2 GPa, suggesting that there exists a fluid/matrix interaction at grain contacts, which reduces the rigidity. The predicted Biot–Gassmann velocity values are greater than the measured velocity values due to the rock–fluid interaction. This is not accounted for in the Biot–Gassmann velocity models and velocity dispersion due to a local flow mechanism. The velocities predicted by the Raymer and time‐average relationships overestimated the measured velocities even more than the Biot model.  相似文献   

8.
Elastic behaviour of North Sea chalk: A well-log study   总被引：1，自引：1，他引：0  
We present two different elastic models for, respectively, cemented and uncemented North Sea chalk well‐log data. We find that low Biot coefficients correlate with anomalously low cementation factors from resistivity measurements at low porosity and we interpret this as an indication of cementation. In contrast, higher Biot coefficients and correspondingly higher cementation factors characterize uncemented chalk for the same (low) porosity. Accordingly, the Poisson's ratio–porosity relationship for cemented chalk is different from that of uncemented chalk. We have tested the application of the self‐consistent approximation, which here represents the unrelaxed scenario where the pore spaces of the rock are assumed to be isolated, and the Gassmann theory, which assumes that pore spaces are connected, as tools for predicting the effect of hydrocarbons from the elastic properties of brine‐saturated North Sea reservoir chalk. In the acoustic impedance–Poisson's ratio plane, we forecast variations in porosity and hydrocarbon saturation from their influence on the elastic behaviour of the chalk. The Gassmann model and the self‐consistent approximation give roughly similar predictions of the effect of fluid on acoustic impedance and Poisson's ratio, but we find that the high‐frequency self‐consistent approach gives a somewhat smaller predicted fluid‐saturation effect on Poisson's ratio than the low‐frequency Gassmann model. The Gassmann prediction for the near and potentially invaded zone corresponds more closely to logging data than the Gassmann prediction for the far, virgin zone. We thus conclude that the Gassmann approach predicts hydrocarbons accurately in chalk in the sonic‐frequency domain, but the fluid effects as recorded by the acoustic tool are significantly affected by invasion of mud filtrate. The amplitude‐versus‐angle (AVA) response for the general North Sea sequence of shale overlying chalk is predicted as a function of porosity and pore‐fill. The AVA response of both cemented and uncemented chalk generally shows a declining reflectivity coefficient versus offset and a decreasing normal‐incidence reflectivity with increasing porosity. However, for the uncemented model, a phase reversal will appear at a relatively lower porosity compared to the cemented model.  相似文献   

9.
We conducted a laboratory study of the joint elastic‐electrical properties of sixty‐three brine‐saturated sandstone samples to assess the likely impact of differential pressure (confining minus pore fluid pressures) in the range 8–60 MPa on the joint interpretation of marine seismic and controlled‐source electromagnetic survey data. The samples showed a wide range of petrophysical properties representative of most sandstone reservoirs. We found that a regression equation comprising both a constant and an exponential part gave a good fit to the pressure dependence of all five measured geophysical parameters (ultrasonic P‐ and S‐wave velocity, attenuation and electrical resistivity). Electrical resistivity was more pressure‐sensitive in clay‐rich sandstones with higher concentrations of low aspect ratio pores and micropores than in clean sandstones. Attenuation was more pressure‐sensitive in clean sandstones with large open pores (macropores) than in clay‐rich sandstones. Pore shape did not show any influence on the pressure sensitivity of elastic velocity. As differential pressure increases, the effect of the low aspect ratio pores and micropores on electrical resistivity becomes stronger than the effect of the macropores on attenuation. Further analysis of correlations among the five parameters as a function of pressure revealed potentially diagnostic relationships for geopressure prediction in reservoir sandstones.  相似文献   

10.
Synthetic rock samples can offer advantages over natural rock samples when used for laboratory rock physical properties studies, provided their success as natural analogues is well understood. The ability of synthetic rocks to mimic the natural stress dependency of elastic wave, electrical and fluid transport properties is of primary interest. Hence, we compare a consistent set of laboratory multi-physics measurements obtained on four quartz sandstone samples (porosity range 20–25%) comprising two synthetic and two natural (Berea and Corvio) samples, the latter used extensively as standards in rock physics research. We measured simultaneously ultrasonic (P- and S-wave) velocity and attenuation, electrical resistivity, permeability and axial and radial strains over a wide range of differential pressure (confining stress 15–50 MPa; pore pressure 5–10 MPa) on the four brine saturated samples. Despite some obvious physical discrepancies caused by the synthetic manufacturing process, such as silica cementation and anisotropy, the results show only small differences in stress dependency between the synthetic and natural sandstones for all measured parameters. Stress dependency analysis of the dry samples using an isotropic effective medium model of spheroidal pores and penny-shaped cracks, together with a granular cohesion model, provide evidence of crack closure mechanisms in the natural sandstones, seen to a much lesser extent in the synthetic sandstones. The smaller grain size, greater cement content, and cementation under oedometric conditions particularly affect the fluid transport properties of the synthetic sandstones, resulting in lower permeability and higher electrical resistivity for a similar porosity. The effective stress coefficients, determined for each parameter, are in agreement with data reported in the literature. Our results for the particular synthetic materials that were tested suggest that synthetic sandstones can serve as good proxies for natural sandstones for studies of elastic and mechanical properties, but should be used with care for transport properties studies.  相似文献   

11.
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.  相似文献   

12.
含流体砂岩地震波频散实验研究   总被引：1，自引：0，他引：1       下载免费PDF全文
为了研究孔隙流体对不同渗透率岩石地震波速度的影响,在实验室利用跨频带岩石弹性参数测试系统得到了应变幅值10-6的2~2000Hz频段下的地震波速度和1 MHz频率下的超声波速度,利用差分共振声谱法得到了频率600Hz岩石干燥和完全饱水情况下岩石声学参数.实验表明,在低饱和度下,致密砂岩在地震和超声频段下没有明显的频散;在高饱和度下纵波速度的频散变得明显.从干燥到完全水饱和条件,不同频率测量的致密砂岩的体积模量随岩石孔隙度增高而降低,且体积模量的变化量受岩石微观孔隙结构的影响较大.高孔、高渗砂岩无论在低含水度下还是在高含水饱和度下频散微弱,并且在地震频段下围压对于岩石纵横波速度的影响要大于频率的影响.高孔、高渗砂岩和致密砂岩不同含水饱和度下的频散差异可应用于储层预测,油气检测等方面,同时该研究可以更好地帮助理解岩石的黏弹性行为,促进岩石物理频散理论的发展,提高地震解释的精度.  相似文献   

13.
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.  相似文献   

14.
We generalize the classical theory of acoustoelasticity to the porous case (one fluid and a solid frame) and finite deformations. A unified treatment of non‐linear acoustoelasticity of finite strains in fluid‐saturated porous rocks is developed on the basis of Biot’s theory. A strain‐energy function, formed with eleven terms, combined with Biot’s kinetic and dissipation energies, yields Lagrange’s equations and consequently the wave equation of the medium. The velocities and dissipation factors of the P‐ and S‐waves are obtained as a function of the 2nd‐ and 3rd‐order elastic constants for hydrostatic and uniaxial loading. The theory yields the limit to the classical theory if the fluid is replaced with a solid with the same properties of the frame. We consider sandstone and obtain results for open‐pore jacketed and closed‐pore jacketed ‘gedanken’ experiments. Finally, we compare the theoretical results with experimental data.  相似文献   

15.
Most sedimentary rocks are anisotropic, yet it is often difficult to accurately incorporate anisotropy into seismic workflows because analysis of anisotropy requires knowledge of a number of parameters that are difficult to estimate from standard seismic data. In this study, we provide a methodology to infer azimuthal P‐wave anisotropy from S‐wave anisotropy calculated from log or vertical seismic profile data. This methodology involves a number of steps. First, we compute the azimuthal P‐wave anisotropy in the dry medium as a function of the azimuthal S‐wave anisotropy using a rock physics model, which accounts for the stress dependency of seismic wave velocities in dry isotropic elastic media subjected to triaxial compression. Once the P‐wave anisotropy in the dry medium is known, we use the anisotropic Gassmann equations to estimate the anisotropy of the saturated medium. We test this workflow on the log data acquired in the North West Shelf of Australia, where azimuthal anisotropy is likely caused by large differences between minimum and maximum horizontal stresses. The obtained results are compared to azimuthal P‐wave anisotropy obtained via orthorhombic tomography in the same area. In the clean sandstone layers, anisotropy parameters obtained by both methods are fairly consistent. In the shale and shaly sandstone layers, however, there is a significant discrepancy between results since the stress‐induced anisotropy model we use is not applicable to rocks exhibiting intrinsic anisotropy. This methodology could be useful for building the initial anisotropic velocity model for imaging, which is to be refined through migration velocity analysis.  相似文献   

16.
Improvements in the joint inversion of seismic and marine controlled source electromagnetic data sets will require better constrained models of the joint elastic‐electrical properties of reservoir rocks. Various effective medium models were compared to a novel laboratory data set of elastic velocity and electrical resistivity (obtained on 67 reservoir sandstone samples saturated with 35 g/l brine at a differential pressure of 8 MPa) with mixed results. Hence, we developed a new three‐phase effective medium model for sandstones with pore‐filling clay minerals based on the combined self‐consistent approximation and differential effective medium model. We found that using a critical porosity of 0.5 and an aspect ratio of 1 for all three components, the proposed model gave accurate model predictions of the observed magnitudes of P‐wave velocity and electrical resistivity and of the divergent trends of clean and clay‐rich sandstones at higher porosities. Using only a few well‐constrained input parameters, the new model offers a practical way to predict in situ porosity and clay content in brine saturated sandstones from co‐located P‐wave velocity and electrical resistivity data sets.  相似文献   

17.
Wave‐induced oscillatory fluid flow in the vicinity of inclusions embedded in porous rocks is one of the main causes for P‐wave dispersion and attenuation at seismic frequencies. Hence, the P‐wave velocity depends on wave frequency, porosity, saturation, and other rock parameters. Several analytical models quantify this wave‐induced flow attenuation and result in characteristic velocity–saturation relations. Here, we compare some of these models by analyzing their low‐ and high‐frequency asymptotic behaviours and by applying them to measured velocity–saturation relations. Specifically, the Biot–Rayleigh model considering spherical inclusions embedded in an isotropic rock matrix is compared with White's and Johnson's models of patchy saturation. The modeling of laboratory data for tight sandstone and limestone indicates that, by selecting appropriate inclusion size, the Biot‐Rayleigh predictions are close to the measured values, particularly for intermediate and high water saturations.  相似文献   

18.
Compressional and shear-wave velocities have been measured and a novel approach using digital processing employed to study wave attenuation for brine- and gas- saturated sandstones, over a range of effective stresses from 5 to 60 MPa. Also measured were the complex conductivity in the brine-saturated state and permeability in the gas-saturated state over the same range of stresses as for the velocity measurements. Broadband ultrasonic pulses of P- and orthogonally polarized S-waves in the frequency range 0.3–0.8 MHz are transmitted through the specimen to be characterized for comparison with a reference (aluminium) having low attenuation. The attenuation is calculated in terms of the quality factor Q from the Fourier spectral ratios, using the frequency spectral ratios technique. The corrections necessary for the effects of diffraction due to the finite size of the ultrasonic transducers have been carried out for the case of measurements under lower confining stress. To interpret the laboratory measured velocity and attenuation data under the physical conditions of this study and to estimate the effects of pore structure, numerical modelling of velocities and attenuation as functions of the confining stress have been performed, based on the MIT model. Theoretical models based on several hypothesized attenuation mechanisms are considered in relation to laboratory data of the effects of confining pressure, fluid saturation and pore structure on attenuation. Numerical calculations using these models with the experimental data indicate that friction on thin cracks and grain boundaries is the dominant attenuation mechanism for dry and brine-saturated sandstones at low effective stresses for the frequencies tested. However, for brine-saturated sandstones at moderately high effective stresses, fluid flow could play a more important role in ultrasonic S-wave attenuation, depending on the pore structure of the sample.  相似文献   

19.
The hydrodynamic characterization of the epikarst, the shallow part of the unsaturated zone in karstic systems, has always been challenging for geophysical methods. This work investigates the feasibility of coupling time‐lapse refraction seismic data with petrophysical and hydrologic models for the quantitative determination of water storage and residence time at shallow depth in carbonate rocks. The Biot–Gassmann fluid substitution model describing the seismic velocity variations with water saturation at low frequencies needs to be modified for this lithology. I propose to include a saturation‐dependent rock‐frame weakening to take into account water–rock interactions. A Bayesian inversion workflow is presented to estimate the water content from seismic velocities measured at variable saturations. The procedure is tested first with already published laboratory measurements on core samples, and the results show that it is possible to estimate the water content and its uncertainty. The validated procedure is then applied to a time‐lapse seismic study to locate and quantify seasonal water storage at shallow depth along a seismic profile. The residence time of the water in the shallow layers is estimated by coupling the time‐lapse seismic measurements with rainfall chronicles, simple flow equations, and the petrophysical model. The daily water input computed from the chronicles is used to constraint the inversion of seismic velocities for the daily saturation state and the hydrodynamic parameters of the flow model. The workflow is applied to a real monitoring case, and the results show that the average residence time of the water in the epikarst is generally around three months, but it is only 18 days near an infiltration pathway. During the winter season, the residence times are three times shorter in response to the increase in the effective rainfall.  相似文献   

20.
We present a new petro-elastical and numerical-simulation methodology to compute synthetic seismograms for reservoirs subject to CO2 sequestration. The petro-elastical equations model the seismic properties of reservoir rocks saturated with CO2, methane, oil and brine. The gas properties are obtained from the van der Waals equation and we take into account the absorption of gas by oil and brine, as a function of the in situ pore pressure and temperature. The dry-rock bulk and shear moduli can be obtained either by calibration from real data or by using rock-physics models based on the Hertz-Mindlin and Hashin-Shtrikman theories. Mesoscopic attenuation due to fluids effects is quantified by using White's model of patchy saturation, and the wet-rock velocities are calculated with Gassmann equations by using an effective fluid modulus to describe the velocities predicted by White's model. The simulations are performed with a poro-viscoelastic modeling code based on Biot's theory, where viscoelasticity is described by generalizing the solid/fluid coupling modulus to a relaxation function. Using the pseudo-spectral method, which allows general material variability, a complete and accurate characterization of the reservoir can be obtained. A simulation, that considers the Utsira sand of the North Sea, illustrates the methodology.  相似文献   

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

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