共查询到20条相似文献，搜索用时 31 毫秒
1.
Seismic driven probabilistic classification of reservoir facies for static reservoir modelling: a case history in the Barents Sea 总被引：2，自引：0，他引：2
Dario Grana Enrico Paparozzi Silvia Mancini Cristiano Tarchiani 《Geophysical Prospecting》2013,61(3):613-629
In this paper we present a case history of seismic reservoir characterization where we estimate the probability of facies from seismic data and simulate a set of reservoir models honouring seismically‐derived probabilistic information. In appraisal and development phases, seismic data have a key role in reservoir characterization and static reservoir modelling, as in most of the cases seismic data are the only information available far away from the wells. However seismic data do not provide any direct measurements of reservoir properties, which have then to be estimated as a solution of a joint inverse problem. For this reason, we show the application of a complete workflow for static reservoir modelling where seismic data are integrated to derive probability volumes of facies and reservoir properties to condition reservoir geostatistical simulations. The studied case is a clastic reservoir in the Barents Sea, where a complete data set of well logs from five wells and a set of partial‐stacked seismic data are available. The multi‐property workflow is based on seismic inversion, petrophysics and rock physics modelling. In particular, log‐facies are defined on the basis of sedimentological information, petrophysical properties and also their elastic response. The link between petrophysical and elastic attributes is preserved by introducing a rock‐physics model in the inversion methodology. Finally, the uncertainty in the reservoir model is represented by multiple geostatistical realizations. The main result of this workflow is a set of facies realizations and associated rock properties that honour, within a fixed tolerance, seismic and well log data and assess the uncertainty associated with reservoir modelling. 相似文献
2.
Estimating gas saturation in a thin layer by using frequency‐dependent amplitude versus offset modelling 
下载免费PDF全文
下载免费PDF全文 Zhaoyu Jin Mark Chapman Xiaoyang Wu Giorgos Papageorgiou 《Geophysical Prospecting》2017,65(3):747-765
Various models have been proposed to link partial gas saturation to seismic attenuation and dispersion, suggesting that the reflection coefficient should be frequency‐dependent in many cases of practical importance. Previous approaches to studying this phenomenon typically have been limited to single‐interface models. Here, we propose a modelling technique that allows us to incorporate frequency‐dependent reflectivity into convolutional modelling. With this modelling framework, seismic data can be synthesised from well logs of velocity, density, porosity, and water saturation. This forward modelling could act as a basis for inversion schemes aimed at recovering gas saturation variations with depth. We present a Bayesian inversion scheme for a simple thin‐layer case and a particular rock physics model and show that, although the method is very sensitive to prior information and constraints, both gas saturation and layer thickness theoretically can be estimated in the case of interfering reflections. 相似文献
3.
Underground fractures play an important role in the storage and movement of hydrocarbon fluid. Fracture rock physics has been the useful bridge between fracture parameters and seismic response. In this paper, we aim to use seismic data to predict subsurface fractures based on rock physics. We begin with the construction of fracture rock physics model. Using the model, we may estimate P-wave velocity, S-wave velocity and fracture rock physics parameters. Then we derive a new approximate formula for the analysis of the relationship between fracture rock physics parameters and seismic response, and we also propose the method which uses seismic data to invert the elastic and rock physics parameters of fractured rock. We end with the method verification, which includes using well-logging data to confirm the reliability of fracture rock physics effective model and utilizing real seismic data to validate the applicability of the inversion method. Tests show that the fracture rock physics effective model may be used to estimate velocities and fracture rock physics parameters reliably, and the inversion method is resultful even when the seismic data is added with random noise. Real data test also indicates the inversion method can be applied into the estimation of the elastic and fracture weaknesses parameters in the target area. 相似文献
4.
Time-lapse seismic data are generally used to monitor the changes in dynamic reservoir properties such as fluid saturation and pore or effective pressure. Changes in saturation and pressure due to hydrocarbon production usually cause changes in the seismic velocities and as a consequence changes in seismic amplitudes and travel times. This work proposes a new rock physics model to describe the relation between saturation-pressure changes and seismic changes and a probabilistic workflow to quantify the changes in saturation and pressure from time-lapse seismic changes. In the first part of this work, we propose a new quadratic approximation of the rock physics model. The novelty of the proposed formulation is that the coefficients of the model parameters (i.e. the saturation-pressure changes) are functions of the porosity, initial saturation and initial pressure. The improvements in the results of the forward model are shown through some illustrative examples. In the second part of the work, we present a Bayesian inversion approach for saturation-pressure 4D inversion in which we adopt the new formulation of the rock physics approximation. The inversion results are validated using synthetic pseudo-logs and a 3D reservoir model for CO2 sequestration. 相似文献
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.
基于岩石物理和地震反演理论,提出了一种同步反演储层孔隙度和含水饱和度的方法.以岩石物理为基础,建立了砂泥岩储层物性和弹性参数之间定量的关系-Simon模型,以贝叶斯理论为手段,结合不同类型的砂泥岩储层,建立了多信息联合约束的物性参数反演目标函数,并通过蒙特卡罗和遗传算法相结合的思路求解该目标函数,最终得到孔隙度和含水饱和度的同步反演结果.将该方法应用于河道砂和砂砾岩两种不同的砂泥岩储层中,孔隙度和含水饱和度数据的联合应用,进一步减少了储层预测的多解性,为石油地质综合研究提供了更加丰富准确的基础数据. 相似文献
7.
By applying seismic inversion, we can derive rock impedance from seismic data. Since it is an interval property, impedance is valuable for reservoir characterization. Furthermore, the decomposition of the impedance into two fundamental properties, i.e. velocity and density, provides a link to the currently available rock‐physics applications to derive quantitative reservoir properties. However, the decomposition is a challenging task due to the strong influence of noise, especially for seismic data with a maximum offset angle of less than 30°. We present a method of impedance decomposition using three elastic impedance data derived from the seismic inversion of angle stacks, where the far‐stack angle is 23.5°. We discuss the effect of noise on the analysis as being the most significant cause of making the decomposition difficult. As the result, the offset‐consistent component of noise mostly affects the determination of density but not the velocities (P‐ and S‐wave), whereas the effect of the random component of noise occurs equally in the determination of the velocities and density. The effect is controlled by the noise enhancement factor 1/A, which is determined by a combination of stack angles. Based on the results of the analysis, we show an innovative method of decomposition incorporating rock‐physics bounds as constraints for the analysis. The method is applied to an actual data set from an offshore oilfield; we demonstrate the result of analysis for sandbody detection. 相似文献
8.
Tight oil siltstones are rocks with complex structure at pore scale and are characterized by low porosity and low permeability at macroscale. The production of tight oil siltstone reservoirs can be increased by hydraulic fracturing. For optimal fracking results, it is desirable to map the ability to fracture based on seismic data prior to fracturing. Brittleness is currently thought to be a key parameter for evaluating the ability to fracture. To link seismic information to the brittleness distribution, a rock physics model is required. Currently, there exists no commonly accepted rock physics model for tight oil siltstones. Based on the observed correlation between porosity and mineral composition and known microstructure of tight oil siltstone in Daqing oilfield of Songliao basin, we develop a rock physics model by combining the Voigt–Reuss–Hill average, self-consistent approximation and differential effective medium theory. This rock physics model allows us to explore the dependence of the brittleness on porosity, mineral composition, microcrack volume fraction and microcrack aspect ratio. The results show that, as quartz content increases and feldspar content decreases, Young's modulus tends to increase and Poisson ratio decreases. This is taken as a signature of higher brittleness. Using well log data and seismic inversion results, we demonstrate the versatility of the rock physics template for brittleness prediction. 相似文献
9.
相对于常规砂岩,致密砂岩在岩石物理性质、力学性质等方面具有明显差异,并呈现出很强的非均质性.岩石物理模型能将储层参数与地震特性信息联系起来,因此可以作为致密砂岩储层参数与地震特性信息转换的桥梁.常规的岩石物理模型通常只考虑单一因素(例如非均匀性,单一孔隙,单一尺度等),建立的岩石物理模板并不适用于致密砂岩.本文针对高饱和气、微裂隙发育、非均质性强的致密砂岩储层,利用Voigt-Reuss-Hill模型计算混合矿物的弹性模量,采用微分等效介质(DEM)模型描述含裂隙、孔隙岩石的骨架弹性模量,基于Biot-Rayleigh波动方程构建了岩石物理弹性模板,给出了致密砂岩储层弹性参数与物性的关系.基于测井数据和实验数据对岩石物理弹性模板进行校正,并将校正后的岩石物理弹性模板结合叠前地震资料应用于川西地区储层孔隙度与裂隙含量预测.结果显示,反演裂隙含量、孔隙度与储层试气报告、测井孔隙度基本吻合,表明该模板能够较合理地应用于致密砂岩储层孔隙度及裂隙含量解释中. 相似文献
10.
岩石物理模板采用测井解释的各类岩性矿物骨架点值,选取适合该地区的岩石物理模型,模拟在不同储层组合、不同孔隙及不同饱和度情况下,储层岩石物理参数变化引起的储层测井参数及地球物理响应特征的变化,定量地建立起储层参数同地球物理弹性参数间的解释关系模版.本文根据新场JS42气藏储层参数分析结果,尝试性地将岩石物理解释模板应用于储层定量解释中,对储层高产气区、含水区域进行定量解释,并预测了该气藏的气水界面,该预测结果与实钻井测试情况吻合,证实了该方法的科学性. 相似文献
11.
By substituting rock skeleton modulus expressions into Gassmann approximate fluid equation, we obtain a seismic porosity inversion equation. However, conventional rock skeleton models and their expressions are quite different from each other, resuling in different seismic porosity inversion equations, potentially leading to difficulties in correctly applying them and evaluating their results. In response to this, a uniform relation with two adjusting parameters suitable for all rock skeleton models is established from an analysis and comparison of various conventional rock skeleton models and their expressions including the Eshelby-Walsh, Pride, Geertsma, Nur, Keys-Xu, and Krief models. By giving the two adjusting parameters specific values, different rock skeleton models with specific physical characteristics can be generated. This allows us to select the most appropriate rock skeleton model based on geological and geophysical conditions, and to develop more wise seismic porosity inversion. As an example of using this method for hydrocarbon prediction and fluid identification, we apply this improved porosity inversion, associated with rock physical data and well log data, to the ZJ basin. Research shows that the existence of an abundant hydrocarbon reservoir is dependent on a moderate porosity range, which means we can use the results of seismic porosity inversion to identify oil reservoirs and dry or water-saturated reservoirs. The seismic inversion results are closely correspond to well log porosity curves in the ZJ area, indicating that the uniform relations and inversion methods proposed in this paper are reliable and effective. 相似文献
12.
13.
Jin-Yong Gui Jian-Hu Gao Xue-Shan Yong Sheng-Jun Li Bin-Yang Liu Wan-Jin Zhao 《应用地球物理》2015,12(4):523-532
Variation of reservoir physical properties can cause changes in its elastic parameters. However, this is not a simple linear relation. Furthermore, the lack of observations, data overlap, noise interference, and idealized models increases the uncertainties of the inversion result. Thus, we propose an inversion method that is different from traditional statistical rock physics modeling. First, we use deterministic and stochastic rock physics models considering the uncertainties of elastic parameters obtained by prestack seismic inversion and introduce weighting coefficients to establish a weighted statistical relation between reservoir and elastic parameters. Second, based on the weighted statistical relation, we use Markov chain Monte Carlo simulations to generate the random joint distribution space of reservoir and elastic parameters that serves as a sample solution space of an objective function. Finally, we propose a fast solution criterion to maximize the posterior probability density and obtain reservoir parameters. The method has high efficiency and application potential. 相似文献
14.
Benyamin Khadem Mohammad Reza Saberi Mohammad Eslahati Bita Arbab 《Geophysical Prospecting》2020,68(5):1613-1632
Rock typing and flow unit detection are more challenging in clastic reservoirs with a uniform pore system. An integrated workflow based on well logs, inverted seismic data and rock physics models is proposed and developed to address such challenges. The proposed workflow supplies a plausible reservoir model for further investigation and adds extra information. Then, this workflow has been implemented in order to define different rock types and flow units in an oilfield in the Persian Gulf, where some of these difficulties have been observed. Here, rock physics models have the leading role in our proposed workflow by providing a diagnostic framework in which we successfully differentiate three rock types with variant characteristics on the given wells. Furthermore, permeability and porosity are calculated using the available rock physics models to define several flow units. Then, we extend our investigation to the entire reservoir by means of simultaneous inversion and rock physics models. The outcomes of the study suggest that in sediments with homogeneous pore size distribution, other reservoir properties such as shale content and cementation (which have distinct effects on the elastic domain) can be used to identify rock types and flow units. These reservoir properties have more physical insights for modelling purposes and can be distinguished on seismic cube using proper rock physics models. The results illustrate that the studied reservoir mainly consists of rock type B, which is unconsolidated sands and has the characteristics of a reservoir for subsequent fluid flow unit analysis. In this regard, rock type B has been divided into six fluid units in which the first detected flow unit is considered as the cleanest unit and has the highest reservoir process speed about 4800 to 5000 mD. Here, reservoir quality decreases from flow unit 1 to flow unit 6. 相似文献
15.
This paper discusses and addresses two questions in carbonate reservoir characterization: how to characterize pore‐type distribution quantitatively from well observations and seismic data based on geologic understanding of the reservoir and what geological implications stand behind the pore‐type distribution in carbonate reservoirs. To answer these questions, three geophysical pore types (reference pores, stiff pores and cracks) are defined to represent the average elastic effective properties of complex pore structures. The variability of elastic properties in carbonates can be quantified using a rock physics scheme associated with different volume fractions of geophysical pore types. We also explore the likely geological processes in carbonates based on the proposed rock physics template. The pore‐type inversion result from well log data fits well with the pore geometry revealed by a FMI log and core information. Furthermore, the S‐wave prediction based on the pore‐type inversion result also shows better agreement than the Greensberg‐Castagna relationship, suggesting the potential of this rock physics scheme to characterize the porosity heterogeneity in carbonate reservoirs. We also apply an inversion technique to quantitatively map the geophysical pore‐type distribution from a 2D seismic data set in a carbonate reservoir offshore Brazil. The spatial distributions of the geophysical pore type contain clues about the geological history that overprinted these rocks. Therefore, we analyse how the likely geological processes redistribute pore space of the reservoir rock from the initial depositional porosity and in turn how they impact the reservoir quality. 相似文献
16.
Wide-azimuth seismic data can be used to derive anisotropic parameters on the subsurface by observing variation in subsurface seismic response along different azimuths. Layer-based high-resolution estimates of components of the subsurface anisotropic elastic tensor can be reconstructed by using wide-azimuth P-wave data by combining the kinematic information derived from anisotropic velocity analysis with dynamic information obtained from amplitude versus angle and azimuth analysis of wide-azimuth seismic data. Interval P-impedance, S-impedance and anisotropic parameters associated with anisotropic fracture media are being reconstructed using linearized analysis assuming horizontal transverse anisotropy symmetry. In this paper it is shown how additional assumptions, such as the rock model, can be used to reduce the degrees of freedom in the estimation problem and recover all five anisotropic parameters. Because the use of a rock model is needed, the derived elastic parameters are consistent with the rock model and are used to infer fractured rock properties using stochastic rock physics inversion. The inversion is based on stochastic rock physics modelling and maximum a posteriori estimate of both porosity and crack density parameters associated with the observed elastic parameters derived from both velocity and amplitude versus angle and azimuth analysis. While the focus of this study is on the use of P-wave reflection data, we also show how additional information such as shear wave splitting and/or anisotropic well log data can reduce the assumptions needed to derive elastic parameter and rock properties. 相似文献
17.
Tor Arne Johansen Åsmund Drottning Isabelle Lecomte Håvar Gjøystdal 《Geophysical Prospecting》2002,50(2):119-137
The aim of seismic reservoir monitoring is to map the spatial and temporal distributions and contact interfaces of various hydrocarbon fluids and water within a reservoir rock. During the production of hydrocarbons, the fluids produced are generally displaced by an injection fluid. We discuss possible seismic effects which may occur when the pore volume contains two or more fluids. In particular, we investigate the effect of immiscible pore fluids, i.e. when the pore fluids occupy different parts of the pore volume. The modelling of seismic velocities is performed using a differential effective‐medium theory in which the various pore fluids are allowed to occupy the pore space in different ways. The P‐wave velocity is seen to depend strongly on the bulk modulus of the pore fluids in the most compliant (low aspect ratio) pores. Various scenarios of the microscopic fluid distribution across a gas–oil contact (GOC) zone have been designed, and the corresponding seismic properties modelled. Such GOC transition zones generally give diffuse reflection regions instead of the typical distinct GOC interface. Hence, such transition zones generally should be modelled by finite‐difference or finite‐element techniques. We have combined rock physics modelling and seismic modelling to simulate the seismic responses of some gas–oil zones, applying various fluid‐distribution models. The seismic responses may vary both in the reflection time, amplitude and phase characteristics. Our results indicate that when performing a reservoir monitoring experiment, erroneous conclusions about a GOC movement may be drawn if the microscopic fluid‐distribution effects are neglected. 相似文献
18.
Exploring seismic inversion methodologies for non‐stationary geological environments: a benchmark study between deterministic and geostatistical seismic inversion 下载免费PDF全文
下载免费PDF全文 This paper presents a comparison between subsurface impedance models derived from different deterministic and geostatistical seismic inversion methodologies applied to a challenging synthetic dataset. Geostatistical seismic inversion methodologies nowadays are common place in both industry and academia, contrasting with traditional deterministic seismic inversion methodologies that are becoming less used as part of the geo‐modelling workflow. While the first set of techniques allows the simultaneous inference of the best‐fit inverse model along with the spatial uncertainty of the subsurface elastic property of interest, the second family of inverse methodology has proven results in correctly predicting the subsurface elastic properties of interest with comparatively less computational cost. We present herein the results of a benchmark study performed over a realistic three‐dimensional non‐stationary synthetic dataset in order to assess the performance and convergence of different deterministic and geostatistical seismic inverse methodologies. We also compare and discuss the impact of the inversion parameterisation over the exploration of the model parameter space. The results show that the chosen seismic inversion methodology should always be dependent on the type and quantity of the available data, both seismic and well‐log, and the complexity of the geological environment versus the assumptions behind each inversion technique. The assessment of the model parameter space shows that the initial guess of traditional deterministic seismic inversion methodologies is of high importance since it will determine the location of the best‐fit inverse solution. 相似文献
19.
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. 相似文献
20.