共查询到20条相似文献，搜索用时 515 毫秒
1.
We obtain the wave velocities of clay-bearing sandstones as a function of clay content, porosity and frequency. Unlike previous theories, based simply on slowness and/or moduli averaging or two-phase models, we use a Biot-type three-phase theory that considers the existence of two solids (sand grains and clay particles) and a fluid. The theory, which is consistent with the critical porosity concept, uses three free parameters that determine the dependence of the dry-rock moduli of the sand and clay matrices as a function of porosity and clay content.
Testing of the model with laboratory data shows good agreement between predictions and measurements. In addition to a rock physics model that can be useful for petrophysical interpretation of wave velocities obtained from well logs and surface seismic data, the model provides the differential equation for computing synthetic seismograms in inhomogeneous media, from the seismic to the ultrasonic frequency bands. 相似文献
Testing of the model with laboratory data shows good agreement between predictions and measurements. In addition to a rock physics model that can be useful for petrophysical interpretation of wave velocities obtained from well logs and surface seismic data, the model provides the differential equation for computing synthetic seismograms in inhomogeneous media, from the seismic to the ultrasonic frequency bands. 相似文献
2.
Yuki Asahara Motohiko Murakami Yasuo Ohishi Naohisa Hirao Kei Hirose 《Earth and Planetary Science Letters》2010,289(3-4):479-485
We extended the pressure range of sound velocity measurements for liquid water to 25 GPa and 900 K along the melting curve using a laser heated diamond anvil cell with a combined system of Brillouin scattering and synchrotron X-ray diffraction. Experimental pressure and temperature were obtained by solving simultaneous equations: the melting curve of ice and the equation of state for gold. The sound velocities obtained in liquid water at high pressures and melting temperatures were converted to density using Murnaghan's equation of state by fitting a parameter of the pressure derivative of bulk modulus at 1 GPa. The results are in good agreement with the values predicted by a previously reported equation of state for water based on sound velocity measurements. The equation of state for water obtained in this study could be applicable to water released by dehydration reactions of dense hydrous magnesium silicate phases in cold subducting slabs at lower mantle conditions, although the validity of Murnaghan's equation of state for water should be evaluated in a wider pressure and temperature ranges. The present velocity data provides the basis for future improvement of the accurate thermodynamic model for water at high pressures. 相似文献
3.
A model has been developed to relate the velocities of acoustic waves Vp and Vs in unconsolidated permafrost to the porosity and extent of freezing of the interstitial water. The permafrost is idealized as an assemblage of spherical quartz grains embedded in a matrix composed of spherical inclusions of water in ice. The wave-scattering theory of Kuster and Toksoz is used to determine the effective elastic moduli, and hence the acoustic velocities. The model predicts Vp and Vs to be decreasing functions of both the porosity and the water-to-ice ratio. The theory has been applied to laboratory measurements of Vp and Vs in 31 permafrost samples from the North American Arctic. Although no direct measurements were made of the extent of freezing in these samples, the data are consistent with the predictions of the model. Electrical resistivity measurements on the permafrost samples have demonstrated their essentially resistive behaviour. The ratio of resistivity of permafrost in its frozen state to that in its unfrozen state has been related to the extent of freezing in the samples. Electromagnetic and seismic reflection surveys can be used together in areas of permafrost: firstly an EM survey to determine the extent of freezing and then the acoustic velocity model to predict the velocities in the permafrost. The necessary transit time corrections can thus be made on seismic reflection records to compensate for the presence of permafrost. 相似文献
4.
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. 相似文献
5.
We obtain the wave velocities and quality factors of gas‐hydrate‐bearing sediments as a function of pore pressure, temperature, frequency and partial saturation. The model is based on a Biot‐type three‐phase theory that considers the existence of two solids (grains and gas hydrate) and a fluid mixture. Attenuation is described with the constant‐Q model and viscodynamic functions to model the high‐frequency behaviour. We apply a uniform gas/water mixing law that satisfies Wood's and Voigt's averages at low and high frequencies, respectively. The acoustic model is calibrated to agree with the patchy‐saturation theory at high frequencies (White's model). Pressure effects are accounted by using an effective stress law for the dry‐rock moduli and permeabilities. The dry‐rock moduli of the sediment are calibrated with data from the Cascadia margin. Moreover, we calculate the depth of the bottom simulating reflector (BSR) below the sea floor as a function of sea‐floor depth, geothermal gradient below the sea floor, and temperature at the sea floor. 相似文献
6.
Yongyang Sun Boris Gurevich Maxim Lebedev Stanislav Glubokovskikh Vassili Mikhaltsevitch Junxin Guo 《Geophysical Prospecting》2019,67(4):888-899
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. 相似文献
7.
Analysis of different parameterisations of waveform inversion of compressional body waves in an elastic transverse isotropic Earth with a vertical axis of symmetry 
下载免费PDF全文
下载免费PDF全文 In a multi‐parameter waveform inversion, the choice of the parameterisation influences the results and their interpretations because leakages and the tradeoff between parameters can cause artefacts. We review the parameterisation selection when the inversion focuses on the recovery of the intermediate‐to‐long wavenumbers of the compressional velocities from the compressional body (P) waves. Assuming a transverse isotropic medium with a vertical axis of symmetry and weak anisotropy, analytical formulas for the radiation patterns are developed to quantify the tradeoff between the shear velocity and the anisotropic parameters and the effects of setting to zero the shear velocity in the acoustic approach. Because, in an anisotropic medium, the radiation patterns depend on the angle of the incident wave with respect to the vertical axis, two particular patterns are discussed: a transmission pattern when the ingoing and outgoing slowness vectors are parallel and a reflection pattern when the ingoing and outgoing slowness vectors satisfy Snell's law. When the inversion aims at recovering the long‐to‐intermediate wavenumbers of the compressional velocities from the P‐waves, we propose to base the parameterisation choice on the transmission patterns. Since the P‐wave events in surface seismic data do not constrain the background (smooth) vertical velocity due to the velocity/depth ambiguity, the preferred parameterisation contains a parameter that has a transmission pattern concentrated along the vertical axis. This parameter can be fixed during the inversion which reduces the size of the model space. The review of several parameterisations shows that the vertical velocity, the Thomsen parameter δ, or the Thomsen parameter ε have a transmission pattern along the vertical axis depending on the parameterisation choice. The review of the reflection patterns of those selected parameterisations should be done in the elastic context. Indeed, when reflection data are also inverted, there are potential leakages of the shear parameter at intermediate angles when we carry out acoustic inversion. 相似文献
8.
基于Kelvin-Voigt黏弹性骨架的含非饱和流体孔隙介质BISQ模型(英文) 总被引：4，自引：1，他引：3
本文综合考虑了在波传播过程中孔隙介质的三种重要力学机制——"Biot流动机制一squirt流动机制-固体骨架黏弹性机制",借鉴等效介质思想,将含水饱和度引入波动力学控制方程,并考虑了不同波频率下孔隙流体分布模式对其等效体积模量的影响,给出了能处理含粘滞性非饱和流体孔隙介质中波传播问题的黏弹性Biot/squirt(BISQ)模型。推导了时间-空间域的波动力学方程组,由一组平面谐波解假设,给出频率-波数域黏弹性BISQ模型的相速度和衰减系数表达式。基于数值算例分析了含水饱和度、渗透率与频率对纵波速度和衰减的影响,并结合致密砂岩和碳酸盐岩的实测数据,对非饱和情况下的储层纵波速度进行了外推,碳酸盐岩储层中纵波速度对含气饱和度的敏感性明显低于砂岩储层。 相似文献
9.
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. 相似文献
10.
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. 相似文献
11.
CLAUDIA KERNER 《Geophysical Prospecting》1990,38(2):111-137
Alekseev and Mikhailenko have developed a wavenumber-summation method which combines a finite integral transformation with a finite-difference calculation and involves no approximations other than numerical ones. However, numerical anisotropy causes velocity errors for shear waves which are unacceptable if Poisson's ratios are larger than 0.4 and unless the number of grid points per wavelength is chosen considerably higher than the value generally regarded as sufficient in finite-difference computations. To overcome this limitation in the applicability of the otherwise very powerful modelling scheme, the method is applied to the elastodynamic equations for the velocity vector. Thus, instead of solving a second-order hyperbolic system as in the case of the wave equation, solutions to a first-order hyperbolic system are computed. The finite-difference iteration is performed in a staggered grid. In addition to mastering the numerical difficulties in cases where the Poisson's ratio is unusually high, this approach results in a code which can be used for the modelling of liquid layers. With the new scheme, water reverberations are investigated in terms of normal modes. It is found that for realistic sea-bottom velocities the critical and supercritical cases exist only for P-waves. It means that compressional waves are trapped within the water layer but energy leaks into the substratum through converted shear waves. These leaky compressional normal modes attain properties similar to those of shear normal modes or Pseudo-Love waves. Due to their origin from conversion of dispersed multi-modal compressional waves the shear waves generated at the sea-bottom form a long complex wavetrain. They were found to mask the reflections from the target horizon in an offset-VSP field section. 相似文献
12.
依据研究区的地热梯度（25℃／km）,在高温高压（最高温度为1050℃,最高压力为1.2GPa）条件下系统测量了横穿红河-哀牢山断裂带的元江-墨江地质剖面上的哀牢山岩群各类变质岩（千枚岩、片岩、浅粒岩、变粒岩、大理岩和片麻岩）的纵波速度．实验结果表明,不同岩类的纵波速度随温度压力变化的趋势不同．在相当于衷牢山岩群变质岩峰期变质温度和压力条件下（P=0．4-0.8GPa,T=35-700℃）,测得大部分岩石的纵波速度为5．50－5．80km／s,这一纵波速度值与区域地球物理测深揭示的中地壳低速层的纵波速度相当因此,结合该区变质岩、地壳内热状态及地球物理测深研究成果可初步认为：组成哀牢山岩群的浅粒岩、变粒岩、酸性片麻岩以及部分千枚岩、片岩为该地区中地壳低速层的主要岩石类型． 相似文献
13.
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. 相似文献
14.
Yuriy Ivanov 《Geophysical Prospecting》2019,67(9):2287-2297
Degeneracies of the slowness surfaces of shear (and compressional) waves in low-symmetry anisotropic media (such as orthorhombic), known as point singularities, pose difficulties during modelling and inversion, but can be potentially used in the latter as model parameter constraints. I analyse the quantity and spatial arrangement of point singularities in orthorhombic media, as well as their relation to the overall strength of velocity anisotropy. A classification scheme based on the number and spatial distribution of singularity directions is proposed. In normal orthorhombic models (where the principal shear moduli are smaller than the principal compressional moduli), point singularities can only be arranged in three distinct patterns, and media with the theoretical minimum (0) and maximum (16) number of singularities are not possible. In orthorhombic models resulting from embedding vertical fractures in transversely isotropic background, only two singularity distributions are possible, in contrast to what was previously thought. Although the total number of singularities is independent of the overall anisotropy strength, for general (non-normal) orthorhombic models, different spatial distributions of singularities become more probable with increasing magnitude of anisotropy. 相似文献
15.
Elastic properties of eclogite rocks from the Bohemian massif 总被引：1，自引：0，他引：1
Vladislav Babuška Jiří Fiala D. J. Mayson R. C. Liebermann Reviewer J. Buben 《Studia Geophysica et Geodaetica》1978,22(4):348-361
Summary Compressional velocity anistropy has been studied in detail at atmospheric pressure for 78 specimens of 23 types of eclogite rocks from the Bohemian massif. For nine of these rock types, compressional and shear velocities were measured as a function of pressure to750 MPa at room temperature. The velocity anisotropy for both compressional and shear waves is less than4% at high pressure. The velocities increase with increasing garnet content and decrease with increasing symplectitization. The Moldanubian eclogites have significantly higher velocities, on the average, than the eclogites from the Kruné hory crystalline complex, although the densities of both groups are comparable. 相似文献
16.
Wyllie's time-average equation and subsequent refinements have been used for over 20 years to estimate the porosity of reservoir rocks from compressional (P)-wave velocity (or its reciprocal, transit time) recorded on a sonic log. This model, while simple, needs to be more convincingly explained in theory and improved in practice, particularly by making use of shear (S)-wave velocity. One of the most important, although often ignored, factors affecting elastic velocities in a rock is pore structure, which is also a controlling factor for transport properties of a rock. Now that S-wave information can be obtained from the sonic log, it may be used with P-waves to provide a better understanding of pore structure. A new acoustic velocities-to-porosity transform based on an elastic velocity model developed by Kuster and Toksöz is proposed. Employing an approximation to an equivalent pore aspect ratio spectrum, pore structure for reservoir rocks is taken into account, in addition to total pore volume. Equidimensional pores are approximated by spheres and rounded spheroids, while grain boundary pores and flat pores are approximated by low aspect ratio cracks. An equivalent pore aspect ratio spectrum is characterized by a power function which is determined by compressional-and shear-wave velocities, as well as by matrix and inclusion properties. As a result of this more sophisticated elastic model of porous rocks and a stricter theory of elastic wave propagation, the new method leads to a more satisfactory interpretation and fuller use of seismic and sonic log data. Calculations using the new transform on data for sedimentary rocks, obtained from published literature and laboratory measurements, are presented and compared at atmospheric pressure with those estimated from the time-average equation. Results demonstrate that, to compensate for additional complexity, the new method provides more detailed information on pore volume and pore structure of reservoir rocks. Examples are presented using a realistic self-consistent averaging scheme to consider interactions between pores, and the possibility of extending the method to complex lithologies and shaly rocks is discussed. 相似文献
17.
J.Peter Watt Richard J. O&#x;Connell 《Physics of the Earth and Planetary Interiors》1980,21(4):359-370
Existing data supporting or disputing the validity of the Hashin-Shtrikman bounds on the elastic properties of multiphase aggregates often do not consider porosity, elastic anisotropy, or experimental errors. In this experiment, two-phase aggregates of KCl + (NH4Br, TlBr, CsCl, NaCl, Cu, and LiF) at every 20% volume fraction were vacuum hot-pressed and the compressional and shear velocities were measured with a computer-controlled ultrasonic interferometer to ±0.2%. The ratio of the shear moduli, μ, (phase 2/KCl) varied from about 1 to 5, producing a range of separations between the theoretical two-phase Hashin-Shtrikman bounds for the composites. Samples were generally 99% or better of the theoretical density, with less than 1% velocity anisotropy. Porosity corrections were applied assuming spherical pores, based on the observed velocity-pressure behaviour. Velocities agreed with the HS bounds calculated from the end-member single-crystal stiffnesses when anisotropy was taken into account.The velocity data were also used to estimate the bulk modulus, K, and shear modulus of the second phase by means of the matrix method — taking the K and μ of KCl as known and calculating the moduli of the other phase assuming that the measured velocities were the two-phase Hashin-Shtrikman bounds or the Voigt-Reuss-Hill average. A narrow range of moduli estimates results only if the μ's of both phases are fairly closely matched. For μ's mismatched by a factor of 5, the theoretical uncertainty in the estimates can be 10 times larger than the experimental uncertainty. Estimates using the VRH average can lie outside the HS-based results. 相似文献
18.
Jos M. Carcione 《Geophysical Prospecting》1998,46(3):249-270
19.
Elastic and electromagnetic waves are commonly used to investigate various soil characteristics. The goal of this study is to estimate the elastic moduli and the void ratio based on both the compressional and shear wave velocities, and the electrical resistivity measured by field velocity resistivity probe (FVRP). The compressional and shear waves are measured by piezoelectric disk elements and bender elements installed at the end of the FVRP frame tip. The electrical resistivity is determined by the electrical resistivity probe installed at the tip of the FVRP frame. The FVRP tests are carried out in a clay–sand mixture prepared in a calibration chamber and in silty sand to silty clay soils in the field. The elastic waves and electrical resistivity are measured at every 1 cm. The field tests are carried out at a depth of 6–20 m, at 10 cm intervals, at the Southern coastal area of the Korean peninsula. The measured data are converted into the constraint and shear moduli based on the elastic waves. Void ratios are evaluated based on the elastic wave velocities and the electrical resistivity, and these void ratios match the volumetric void ratio well. This study suggests that the FVRP may effectively determine the elastic moduli and void ratio. 相似文献
20.
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. 相似文献