共查询到20条相似文献，搜索用时 250 毫秒
1.
储层弹性与物性参数可直接应用于储层岩性预测和流体识别,是储层综合评价和油气藏精细描述的基本要素之一.现有的储层弹性与物性参数地震同步反演方法大都基于Gassmann方程,使用地震叠前数据,通过随机优化方法反演储层弹性与物性参数;或基于Wyllie方程,使用地震叠后数据,通过确定性优化方法反演储层弹性与物性参数.本文提出一种基于Gassmann方程、通过确定性优化方法开展储层弹性和物性参数地震叠前反演的方法,该方法利用Gassmann方程建立储层物性参数与叠前地震观测数据之间的联系,在贝叶斯反演框架下以储层弹性与物性参数的联合后验概率为目标函数,通过将目标函数的梯度用泰勒公式展开得到储层弹性与物性参数联合的方程组,其中储层弹性参数对物性参数的梯度用差分形式表示,最后通过共轭梯度算法迭代求解得到储层弹性与物性参数的最优解.理论试算与实际资料反演结果证明了方法的可行性. 相似文献
2.
The added value of the joint pre-stack inversion of PP (incident P-wave and reflected P-wave) and PS (incident P-wave and reflected S-wave) seismic data for the time-lapse application is shown. We focus on the application of this technique to the time-lapse (four-dimensional) multicomponent Jubarte field permanent reservoir monitoring seismic data. The joint inversion results are less sensitive to noise in the input data and show a better match with the rock physics models calibrated for the field. Further, joint inversion improves S-impedance estimates and provides a more robust quantitative interpretation, allowing enhanced differentiation between pore pressure and fluid saturation changes, which will be extremely useful for reservoir management. Small changes in reservoir properties are expected in the short time between the time-lapse seismic acquisitions used in the Jubarte project (only 1 year apart). The attempt to recover subtle fourth-dimensional effects via elastic inversion is recurrent in reservoir characterization projects, either due to the small sensitivity of the reservoirs to fluid and pressure changes or the short interval between the acquisitions. Therefore, looking for methodologies that minimize the uncertainty of fourth-dimensional inversion outputs is of fundamental importance. Here, we also show the differences between PP only and joint PP–PS inversion workflows and parameterizations that can be applied in other projects. We show the impact of using multicomponent data as input for elastic seismic inversions in the analysis of the time-lapse differences of the elastic properties. The larger investment in the acquisition and processing of multicomponent seismic data is shown to be justified by the improved results from the fourth-dimensional joint inversion. 相似文献
3.
Anne-Kari Furre Mona Andersen Anita Smalø Moen Randi Kristin Tønnessen 《Geophysical Prospecting》2009,57(3):427-437
We have developed a method to determine the effect of pore pressure depletion on elastic framework moduli, by utilizing sonic logs from wells drilled at different locations through a reservoir at varying depletion stages. This is done by first inverting the sonic logs for elastic framework bulk and shear moduli, thus carefully removing pressure dependent fluid effects. By crossplotting these elastic framework moduli against an increase in net stress (which is directly related to depletion), we derive the stress sensitivity of the elastic framework moduli. We found that the observed stress sensitivity was consistent with time-lapse seismic results and that the sensitivity derived from the sonic logs was much smaller than predicted by hydrostatic measurements on core samples. This method is applicable to depletion scenarios in mature fields and can be used both for modelling and inverting time-lapse seismic data for effects of pore pressure depletion on seismic data. 相似文献
4.
B. Gurevich 《Pure and Applied Geophysics》2000,157(5):811-827
5.
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. 相似文献
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.
基于流体替换技术的地震AVO属性气藏识别(英文) 总被引：2，自引：1，他引：1
传统上,油藏地球物理工程师是基于测井数据进行流体替换,计算油藏饱和不同流体时的弹性参数,并通过地震正演模拟分析油藏饱和不同流体时的地震响应,从而进行油气藏识别研究。该研究方案为油藏研究提供了重要的弹性参数和地震响应信息,但这些信息仅限于井眼位置。对于实际油藏条件,地下储层参数都是随位置变化而变化的,如孔隙度、泥质含量和油藏厚度等,因此基于传统流体替换方案得到的流体变化地震响应信息对于油气藏识别具有很大的局限性。研究通过设定联系油藏弹性参数与孔隙度、矿物组分等参数的岩石物理模型,并基于三层地质模型,进行地震正演模拟与AVO属性计算。得到油藏孔隙度、泥质含量和储层厚度变化时地震AVO属性,并建立了饱和水储层和含气储层对应AVO属性(包括梯度与截距)之间的定量关系。建立的AVO属性之间的线性关系可以实现基于地震AVO属性直接进行流体替换。最后,应用建立的流体替换前后AVO属性之间线性方程,对模拟地震数据直接进行流体替换,并通过流体替换前后AVO属性交汇图分析实现了气藏识别。 相似文献
8.
地震尺度下碳酸盐岩储层的岩石物理建模方法(英文) 总被引：3，自引：3，他引：0
碳酸盐岩油藏的强非均质性以及孔隙结构的复杂性,使得作为连接油藏参数与地震参数重要桥梁的岩石物理模型,以及作为油藏预测和定量表征最有效工具的流体替换成为岩石物理建模的难点与重点。在碳酸盐岩储层复杂孔隙结构与地震尺度下碳酸盐岩储层非均质性分析基础上,研究采用岩石网格化方法,将地震尺度下非均质碳酸盐岩储层岩石划分为具有独立岩石参数的均质岩石子体,根据岩石孔隙成因与结构特征采用不同岩石物理模型分步计算岩石子块干岩石弹性模量,并根据不同孔隙连通性进行流体替换,计算饱和不同流体岩石弹性模量。基于计算的岩石子块弹性模量,采用Hashin-Shtrikman-Walpole弹性边界计算理论方法实现地震尺度下碳酸盐岩储层弹性参数计算。通过对含有不同类型孔隙组合碳酸盐岩储层模型的弹性模量进行计算与分析,明确不同孔隙对岩石弹性参数的影响特征,模拟分析结果与实际资料认识一致。 相似文献
9.
Larry Lines 《Geophysical Prospecting》2015,63(2):422-441
Unlike light oils, heavy oils do not have a well‐established scheme for modelling elastic moduli from dynamic reservoir properties. One of the main challenges in the fluid substitution of heavy oils is their viscoelastic nature, which is controlled by temperature, pressure, and fluid composition. Here, we develop a framework for fluid substitution modelling that is reliable yet practical for a wide range of cold and thermal recovery scenarios in producing heavy oils and that takes into account the reservoir fluid composition, grounded on the effective‐medium theories for estimating elastic moduli of an oil–rock system. We investigate the effect of fluid composition variations on oil–rock elastic moduli with temperature changes. The fluid compositional behaviour is determined by flash calculations. Elastic moduli are then determined using the double‐porosity coherent potential approximation method and the calculated viscosity based on the fluid composition. An increase in temperature imposes two opposing mechanisms on the viscosity behaviour of a heavy‐oil sample: gas liberation, which tends to increase the viscosity, and melting, which decreases the viscosity. We demonstrate that melting dominates gas liberation, and as a result, the viscosity and, consequently, the shear modulus of the heavy oils always decrease with increasing temperature. Furthermore, it turns out that one can disregard the effects of gas in the solution when modelling the elastic moduli of heavy oils. Here, we compare oil–rock elastic moduli when the rock is saturated with fluids that have different viscosity levels. The objective is to characterize a unique relation between the temperature, the frequency, and the elastic moduli of an oil–rock system. We have proposed an approach that takes advantage of this relation to find the temperature and, consequently, the viscosity in different regions of the reservoir. 相似文献
10.
岩石物理测量是油藏水驱开采时移地震监测的基础.在实验室对来自胜利油田的5块岩石样品模拟储层条件进行了水驱和气驱动态岩石物理弹性测量,重点分析了流体替换、温度、孔隙压力对岩石纵、横波速度的影响.实验表明,在水驱情形下,由于流体替换和温度、孔隙压力变化所引起的岩石纵横波速度的变化均很小,实施时移地震监测具有较大的风险性.相比之下,气驱可能引起较为明显的纵波速度变化,有利于时移地震监测的实施.进一步完善实验方法、丰富实验内容、是今后时移地震岩石物理实验研究的主要任务. 相似文献
11.
在双相介质理论指导下构建的流体因子Gassmann流体项f可以实现对油水层的有效判识。本文基于包含流体因子f的反射系数近似公式推导,得到相应的流体弹性阻抗方程,将流体因子纳入地震反演过程并直接参与目标储层的流体识别,实现了从弹性阻抗数据体中直接反演流体因子f的方法技术流程。流体因子f的直接反演有效减少间接计算所带来的累计误差,提高了地震流体识别的精度。在实际资料应用中,结合研究区地质背景与地球物理特征开展的岩石物理分析表明,Gassmann流体因子f对含油和含水特征有较好的区分作用。反演结果表明,流体因子f能为有效地识别目的层的油水特征,为下一步的钻井开发提供可靠的技术支持。 相似文献
12.
The elastic moduli of four sandstone samples are measured at seismic (2?2000 Hz) and ultrasonic (1 MHz) frequencies and water- and glycerin-saturated conditions. We observe that the high-permeability samples under partially water-saturated conditions and the low-permeability samples under partially glycerin-saturated conditions show little dispersion at low frequencies (2?2000 Hz). However, the high-permeability samples under partially glycerin-saturated conditions and the low-permeability samples under partially water-saturated conditions produce strong dispersion in the same frequency range (2?2000 Hz). This suggests that fluid mobility largely controls the pore-fluid movement and pore pressure in a porous medium. High fluid mobility facilitates pore-pressure equilibration either between pores or between heterogeneous regions, resulting in a low-frequency domain where the Gassmann equations are valid. In contrast, low fluid mobility produces pressure gradients even at seismic frequencies, and thus dispersion. The latter shows a systematic shift to lower frequencies with decreasing mobility. Sandstone samples showed variations in Vp as a function of fluid saturation. We explore the applicability of the Gassmann model on sandstone rocks. Two theoretical bounds for the P-velocity are known, the Gassmann–Wood and Gassmann–Hill limits. The observations confirm the effect of wave-induced flow on the transition from the Gassmann–Wood to the Gassmann–Hill limit. With decreasing fluid mobility, the P-velocity at 2–2000 Hz moves from the Gassmann–Wood boundary to the Gassmann–Hill boundary. In addition,, we investigate the mechanisms responsible for this transition. 相似文献
13.
引入基于Gassmann方程的流体替换方法,在分析地震波P波速度、波速比与岩石孔隙度和饱和度关系的基础上,应用于珊溪水库地震波速比和P波速度变化特征研究,得到:(1)珊溪水库震中区岩石始终处于接近水饱和的饱水状态,波速比和P波速度"下降-回升"的变化实质上反映了震中区岩石"孔隙度增大(饱和度减小)-饱和度增大"的变化,每一丛地震的波速比由极小值逐渐增大为极大值是由于岩石从不饱和状态变化到饱和状态;(2)根据每一丛地震波速比的变化,计算得到珊溪水库流体扩散率αs=1.06×104 cm2 s-1,该数值与美国南卡罗莱纳水库、巴西Acu水库、广东新丰江水库的流体扩散率基本一致;(3)震源区岩石孔隙度上限值为8.7%~2.0%,该数值与华东勘测设计研究院通过室内岩石物理力学性质试验测定的珊溪水库坝址区新鲜流纹斑岩的孔隙度平均值一致。 相似文献
14.
Sandra Barbouteau Patrick Rasolofosaon Noalwenn Dubos‐Sallée Virginie Vassil Jean‐François Nauroy 《Geophysical Prospecting》2017,65(1):240-250
An alternative laboratory technique to measure the elastic constants of solid samples, based on the analysis of the cross‐correlation spectra of the vibratory response of randomly excited short solid cylinders, has been recently proposed. The aim of this paper is to check the ability of the technique called passive ultrasonic interferometry to monitor fluid substitution in different rock samples. Velocity variations due to fluid substitution are easily measured if the wave attenuation in the fluid‐saturated rock is not too large (typically in rocks with few cracks or microfractures). The experimental results are in agreement with the predictions of Biot–Gassmann poroelastic theory. The effect of substituting water with a stiffer saturating fluid, such as ethylene glycol, is to increase the overall bulk modulus of the rock, without any substantial effect on shear modulus. Furthermore, the experimental results compare well with those obtained independently with conventional pulse‐transmission technique using ultrasonic transducers. However, the measured pulse‐transmission bulk moduli are slightly larger than the corresponding measured ultrasonic interferometry moduli, with the deviation increasing with increasing fluid viscosity. This can be explained by dispersion due to wave‐induced flow of the viscous fluid since pulse‐transmission experiments involve higher frequencies than ultrasonic interferometry experiments. 相似文献
15.
非重复采集时移地震正演模拟及可行性分析 总被引：1，自引：0，他引：1
本文主要研究前期为常规三维地震勘探和后期为高精度三维地震勘探的这种典型非重复采集时移地震的可行性,主要使用正演模拟和匹配处理相结合的方法进行.本文推导了基于PML边界的三维声波波动方程时空域高阶有限差分正演方程,根据岩石物理流体替代分析方法计算得到时移前后储层的弹性参数,建立了时移地震地球物理模型,设计了常规三维和高精度三维观测系统,进行地震正演模拟,获得了非重复采集时移地震理论模拟数据.利用匹配处理法对非重复采集时移地震正演模拟数据和实际资料进行可行性分析,匹配处理过程主要包括初始数据匹配、一致性处理和互均化处理等.理论模拟数据和实际资料的匹配处理结果表明,通过针对性的分析和处理,利用非重复采集的地震资料进行时移地震研究是可行的,能够得到较好的时移地震响应,在一定程度上反映了油藏的变化情况. 相似文献
16.
Maxim Lebedev Moyra E.J. Wilson Vassily Mikhaltsevitch 《Geophysical Prospecting》2014,62(6):1253-1265
Petrophysical properties of carbonate reservoirs are less predictable than that of siliciclastic reservoirs. One of the main reasons for this is the physical and chemical interactions of carbonate rocks with pore fluids. Such interactions can significantly change the elastic properties of the rock matrix and grains, making the applicability of Gassmann's fluid substitution procedure debatable. This study is an attempt to understand the mechanisms of fluid‐rock interactions and the influence of these interactions on elastic parameters of carbonates. We performed precise indentation tests on Savonnières limestone at a microscale level under dry, distilled water, and n‐Decane saturated conditions. Our experiments display softening of the rock matrix after water saturation. We have found that mainly the ooid cortices, peloid nuclei and prismatic intergranular cement are affected by water flooding. We also observed a shear modulus reduction in Savonnières limestone in an experiment performed at ultrasonic frequencies. One of the most important results obtained in our experimental study is that the Gassmann fluid substitution theory might not always be applicable to predict the elastic moduli of fluid‐saturated limestones. 相似文献
17.
Pore structure and mineral matrix elastic moduli are indispensable in rock physics models. We propose an estimation method of pore structure and mineral moduli based on Kuster-Toksöz model and Biot’s coefficient. In this technique, pore aspect ratios of five different scales from 100 to 10?4 are considered, Biot’s coefficient is used to determine bounds of mineral moduli, and an estimation procedure combined with simulated annealing (SA) algorithm to handle real logs or laboratory measurements is developed. The proposed method is applied to parameter estimations on 28 sandstone samples, the properties of which have been measured in lab. The water saturated data are used for estimating pore structure and mineral moduli, and the oil saturated data are used for testing these estimated parameters through fluid substitution in Kuster-Toksöz model. We then compare fluid substitution results with lab measurements and find that relative errors of P-wave and S-wave velocities are all less than 5%, which indicates that the estimation results are accurate. 相似文献
18.
Nisar Ahmed Perveiz Khalid Hafiz Muhammad Bilal Shafi Patrick Connolly 《Acta Geophysica》2017,65(5):991-1007
The use of seismic direct hydrocarbon indicators is very common in exploration and reservoir development to minimise exploration risk and to optimise the location of production wells. DHIs can be enhanced using AVO methods to calculate seismic attributes that approximate relative elastic properties. In this study, we analyse the sensitivity to pore fluid changes of a range of elastic properties by combining rock physics studies and statistical techniques and determine which provide the best basis for DHIs. Gassmann fluid substitution is applied to the well log data and various elastic properties are evaluated by measuring the degree of separation that they achieve between gas sands and wet sands. The method has been applied successfully to well log data from proven reservoirs in three different siliciclastic environments of Cambrian, Jurassic, and Cretaceous ages. We have quantified the sensitivity of various elastic properties such as acoustic and extended elastic (EEI) impedances, elastic moduli (K sat and K sat–μ), lambda–mu–rho method (λρ and μρ), P-to-S-wave velocity ratio (V P/V S), and Poisson’s ratio (σ) at fully gas/water saturation scenarios. The results are strongly dependent on the local geological settings and our modeling demonstrates that for Cambrian and Cretaceous reservoirs, K sat–μ, EEI, V P/V S, and σ are more sensitive to pore fluids (gas/water). For the Jurassic reservoir, the sensitivity of all elastic and seismic properties to pore fluid reduces due to high overburden pressure and the resultant low porosity. Fluid indicators are evaluated using two metrics: a fluid indicator coefficient based on a Gaussian model and an overlap coefficient which makes no assumptions about a distribution model. This study will provide a potential way to identify gas sand zones in future exploration. 相似文献
19.
天然气在开发过程中,储层有效压力和含气饱和度均会发生变化,研究有效压力和含气饱和度的变化对地震响应特征的影响,在基于时移地震的剩余气分布预测研究中具有重要意义。天然气和石油的声学性质有着明显的差异,油藏时移地震的研究成果不能直接应用于气藏,因此需要开展气藏的时移地震研究。利用Shapiro模型表征干岩石弹性模量随有效压力的变化,借助Batzle-Wang方程描述流体速度随压力的变化关系,联合Gassmann理论进行流体替代,表征饱和流体岩石速度随含气饱和度的变化,建立了饱和流体岩石速度随有效压力和饱和度变化的岩石物理模型。基于该模型,对不同含气饱和度和不同有效压力下的气藏储层模型进行了多波时移地震叠前振幅变化（AVO）模拟。结果表明多波时移地震AVO技术可以有效地区分有效压力变化和含气饱和度变化,为进一步开展气藏多波时移地震流体监测提供了理论参考依据。 相似文献
20.
Fracture detection and fluid identification are important tasks for a fractured reservoir characterization. Our goal is to demonstrate a direct approach to utilize azimuthal seismic data to estimate fluid bulk modulus, porosity, and dry fracture weaknesses, which decreases the uncertainty of fluid identification. Combining Gassmann’s (Vier. der Natur. Gesellschaft Zürich 96:1–23, 1951) equations and linear-slip model, we first establish new simplified expressions of stiffness parameters for a gas-bearing saturated fractured rock with low porosity and small fracture density, and then we derive a novel PP-wave reflection coefficient in terms of dry background rock properties (P-wave and S-wave moduli, and density), fracture (dry fracture weaknesses), porosity, and fluid (fluid bulk modulus). A Bayesian Markov chain Monte Carlo nonlinear inversion method is proposed to estimate fluid bulk modulus, porosity, and fracture weaknesses directly from azimuthal seismic data. The inversion method yields reasonable estimates in the case of synthetic data containing a moderate noise and stable results on real data. 相似文献