首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Carbon capture and storage is a viable greenhouse gas mitigation technology and the Sleipner CO2 sequestration site in the North Sea is an excellent example. Storage of CO2 at the Sleipner site requires monitoring over large areas, which can successfully be accomplished with time lapse seismic imaging. One of the main goals of CO2 storage monitoring is to be able to estimate the volume of the stored CO2 in the reservoir. This requires a parametrization of the subsurface as exact as possible. Here we use elastic 2D time‐domain full waveform inversion in a time lapse manner to obtain a P‐wave velocity constrain directly in the depth domain for a base line survey in 1994 and two post‐injection surveys in 1999 and 2006. By relating velocity change to free CO2 saturation, using a rock physics model, we find that at the considered location the aquifer may have been fully saturated in some places in 1999 and 2006.  相似文献   

2.
In the Norwegian North Sea, the Sleipner field produces gas with a high CO2 content. For environmental reasons, since 1996, more than 11 Mt of this carbon dioxide (CO2) have been injected in the Utsira Sand saline aquifer located above the hydrocarbon reservoir. A series of seven 3D seismic surveys were recorded to monitor the CO2 plume evolution. With this case study, time‐lapse seismics have been shown to be successful in mapping the spread of CO2 over the past decade and to ensure the integrity of the overburden. Stratigraphic inversion of seismic data is currently used in the petroleum industry for quantitative reservoir characterization and enhanced oil recovery. Now it may also be used to evaluate the expansion of a CO2 plume in an underground reservoir. The aim of this study is to estimate the P‐wave impedances via a Bayesian model‐based stratigraphic inversion. We have focused our study on the 1994 vintage before CO2 injection and the 2006 vintage carried out after a CO2 injection of 8.4 Mt. In spite of some difficulties due to the lack of time‐lapse well log data on the interest area, the full application of our inversion workflow allowed us to obtain, for the first time to our knowledge, 3D impedance cubes including the Utsira Sand. These results can be used to better characterize the spreading of CO2 in a reservoir. With the post‐stack inversion workflow applied to CO2 storage, we point out the importance of the a priori model and the issue to obtain coherent results between sequential inversions of different seismic vintages. The stacking velocity workflow that yields the migration model and the a priori model, specific to each vintage, can induce a slight inconsistency in the results.  相似文献   

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

4.
Seismic monitoring of sequestered carbon dioxide (CO2) in underground deposits is a matter of growing importance. The subsurface monitoring of this greenhouse gas is possible due to the marked contrast between the physical properties of natural reservoir fluids and those of carbon dioxide after the injection. This technique makes necessary the investigation of appropriate seismic indicators to link seismic attributes to petrophysical properties, composition and state of the rock as well as pore-fluid type and in-situ physical conditions. With this motivation in mind, we use a Biot–Gassmann formulation to model the theoretical P-wave amplitude reflection coefficients vs. angle of incidence in the seismic range when a planar P-wave strikes the interface between a caprock and a porous sandstone which has its pore space saturated by a mixture of CO2 with brine or oil at different states (supercritical, liquid and gas). The effects of dissolution of CO2 in oil and the existence of a saturation threshold, above which a free CO2 phase develops, are included in the computations. Attention is particularly focused on the sensitivity of the classic best-fit amplitude variations with angle coefficients, to different degrees of CO2 saturation. We conclude from this analysis that the changes in seismic AVA attributes between 30 and 40 degrees can be useful to infer bounds on the CO2 saturation degree, to detect the presence of immiscible CO2 phase and, in some cases, to infer the physical state of the accumulations.  相似文献   

5.
Time-lapse seismic data is useful for identifying fluid movement and pressure and saturation changes in a petroleum reservoir and for monitoring of CO2 injection. The focus of this paper is estimation of time-lapse changes with uncertainty quantification using full-waveform inversion. The purpose of also estimating the uncertainty in the inverted parameters is to be able to use the inverted seismic data quantitatively for updating reservoir models with ensemble-based methods. We perform Bayesian inversion of seismic waveform data in the frequency domain by combining an iterated extended Kalman filter with an explicit representation of the sensitivity matrix in terms of Green functions (acoustic approximation). Using this method, we test different strategies for inversion of the time-lapse seismic data with uncertainty. We compare the results from a sequential strategy (making a prior from the monitor survey using the inverted baseline survey) with a double difference strategy (inverting the difference between the monitor and baseline data). We apply the methods to a subset of the Marmousi2 P-velocity model. Both strategies performed well and relatively good estimates of the monitor velocities and the time-lapse differences were obtained. For the estimated time-lapse differences, the double difference strategy gave the lowest errors.  相似文献   

6.
CO2 has been injected into the saline aquifer Utsira Fm at the Sleipner field since 1996. In order to monitor the movement of the CO2 in the sub‐surface, the seventh seismic monitor survey was acquired in 2010, with dual sensor streamers which enabled optimal towing depths compared to previous surveys. We here report both on the time‐lapse observations and on the improved resolution compared to the conventional streamer surveys. This study shows that the CO2 is still contained in the subsurface, with no indications of leakage. The time‐lapse repeatability of the dual sensor streamer data versus conventional data is sufficient for interpreting the time‐lapse effects of the CO2 at Sleipner, and the higher resolution of the 2010 survey has enabled a refinement of the interpretation of nine CO2 saturated layers with improved thickness estimates of the layers. In particular we have estimated the thickness of the uppermost CO2 layer based on an analysis of amplitude strength together with time‐separation of top and base of this layer and found the maximum thickness to be 11 m. This refined interpretation gives a good base line for future time‐lapse surveys at the Sleipner CO2 injection site.  相似文献   

7.
The effect of sub‐core scale heterogeneity on fluid distribution pattern, and the electrical and acoustic properties of a typical reservoir rock was studied by performing drainage and imbibition flooding tests with CO2 and brine in a laboratory. Moderately layered Rothbach sandstone was used as a test specimen. Two core samples were drilled; one perpendicular and the other parallel to the layering to allow injection of fluids along and normal to the bedding plane. During the test 3D images of fluid distribution and saturation levels were mapped by an industrial X‐ray CT‐scanner together with simultaneous measurement of electrical resistivity, ultrasonic velocities as well as amplitudes. The results showed how the layering and the flooding direction influenced the fluid distribution pattern and the saturation level of the fluids. For a given fluid saturation level, the measured changes in the acoustic and electrical parameters were affected by both the fluid distribution pattern and the layering orientation relative to the measurement direction. The P‐wave amplitude and the electrical resistivity were more sensitive to small changes in the fluid distribution patterns than the P‐wave velocity. The change in amplitude was the most affected by the orientation of the layering and the resulting fluid distribution patterns. In some instances the change due to the fluid distribution pattern was higher than the variation caused by the change in CO2 saturation. As a result the Gassmann relation based on ‘uniform' or ‘patchy' saturation pattern was not suitable to predict the P‐wave velocity variation. Overall, the results demonstrate the importance of core‐imaging to improve our understanding of fluid distribution patterns and the associated effects on measured rock‐physics properties.  相似文献   

8.
The vertical stratification of carbon dioxide (CO2) injected into a deep layered aquifer made up of high-permeability and low-permeability layers, such as Utsira aquifer at Sleipner site in Norway, is investigated with a Buckley–Leverett equation including gravity effects. In a first step, we study both by theory and simulation the application of this equation to the vertical migration of a light phase (CO2), in a denser phase (water), in 1D vertical columns filled with different types of porous media: homogeneous, piecewise homogeneous, layered periodic and finally heterogeneous. For each case, we solve the associated Riemann problems and propose semi-analytical solutions describing the spatial and temporal evolution of the light phase saturation. These solutions agree well with simulation results. We show that the flux continuity condition at interfaces between high-permeability and low-permeability layers leads to CO2 saturation discontinuities at these interfaces and, in particular, to a saturation increase beneath low-permeability layers. In a second step, we analyze the vertical migration of a CO2 plume injected into a 2D layered aquifer. We show that the CO2 vertical stratification under each low-permeability layer is induced, as in 1D columns, by the flux continuity condition at interfaces. As the injection takes place at the bottom of the aquifer the velocity and the flux function decrease with elevation and this phenomenon is proposed to explain the stratification under each mudstone layer as observed at Sleipner site.  相似文献   

9.
Shear-wave velocity is a key parameter for calibrating monitoring time-lapse 4D seismic data during CO2-EOR (Enhanced Oil Recovery) and CO2 sequestration. However, actual S-wave velocity data are lacking, especially in 4D data for CO2 sequestration because wells are closed after the CO2 injection and seismic monitoring is continued but no well log data are acquired. When CO2 is injected into a reservoir, the pressure and saturation of the reservoirs change as well as the elastic parameters of the reservoir rocks. We propose a method to predict the S-wave velocity in reservoirs at different pressures and porosities based on the Hertz–Mindlin and Gassmann equations. Because the coordination number is unknown in the Hertz–Mindlin equation, we propose a new method to predict it. Thus, we use data at different CO2 injection stages in the Gao89 well block, Shengli Oilfield. First, the sand and mud beds are separated based on the structural characteristics of the thin sand beds and then the S-wave velocity as a function of reservoir pressure and porosity is calculated. Finally, synthetic seismic seismograms are generated based on the predicted P- and S-wave velocities at different stages of CO2 injection.  相似文献   

10.
引入基于Gassmann方程的流体替换方法,在分析地震波P波速度、波速比与岩石孔隙度和饱和度关系的基础上,应用于珊溪水库地震波速比和P波速度变化特征研究,得到:(1)珊溪水库震中区岩石始终处于接近水饱和的饱水状态,波速比和P波速度"下降-回升"的变化实质上反映了震中区岩石"孔隙度增大(饱和度减小)-饱和度增大"的变化,每一丛地震的波速比由极小值逐渐增大为极大值是由于岩石从不饱和状态变化到饱和状态;(2)根据每一丛地震波速比的变化,计算得到珊溪水库流体扩散率αs=1.06×104 cm2 s-1,该数值与美国南卡罗莱纳水库、巴西Acu水库、广东新丰江水库的流体扩散率基本一致;(3)震源区岩石孔隙度上限值为8.7%~2.0%,该数值与华东勘测设计研究院通过室内岩石物理力学性质试验测定的珊溪水库坝址区新鲜流纹斑岩的孔隙度平均值一致。  相似文献   

11.
二氧化碳地质封存是减少温室气体排放和减缓温室效应的重要手段.二氧化碳封存的一个重要组成部分是地震监测,即用地震的方法监测封存后的二氧化碳的分布变化.为了实现这个目标,需要建立储层参数与地震性质之间的关系(岩石物理模型)和从地震监测数据中反演获得储层流体的饱和度等参数.首先,本文以Biot理论为基础,结合多相流模型研究了多个物理参数(孔隙度、二氧化碳饱和度、温度和压力等)对同时含有二氧化碳和水的孔隙介质的波速和衰减等属性的影响.结果表明:孔隙度和二氧化碳饱和度对岩石的频散和衰减属性影响强烈,而温度和压力通过孔隙流体性质对岩石的波速产生影响.然后,本文基于含多相流的Biot理论,应用抗干扰能力强、且具有更好的局部搜索能力和抗早熟能力的自适应杂交遗传算法对实际数据进行了反演研究.对岩心实验数据的反演研究表明了算法的有效性,而且表明含多相流的Biot理论能够很好地解释水和二氧化碳饱和岩石的波速特征.最后,我们将自适应杂交遗传算法应用于实际封存项目的地震监测数据,获得了封存后不同时期的二氧化碳饱和度,达到了用地震方法监测二氧化碳分布的目的.  相似文献   

12.
Between the years 2008 and 2013, approximately 67 kilotons of CO2 have been injected at the Ketzin site, Germany. As part of the geophysical monitoring programme, time‐lapse electrical resistivity tomography has been applied using crosshole and surface‐downhole measurements of electrical resistivity tomography. The data collection of electrical resistivity tomography is partly based on electrodes that are permanently installed in three wells at the site (one injection well and two observation wells). Both types of ERT measurements consistently show the build‐up of a CO2‐related resistivity signature near the injection point. Based on the imaged resistivity changes and a petrophysical model, CO2 saturation levels are estimated. These CO2 saturations are interpreted in conjunction with CO2 saturations inferred from neutron‐gamma loggings. Apart from the CO2–brine substitution response in the observed resistivity changes, significant imprints from the dynamic behaviour of the CO2 in the reservoir are observed.  相似文献   

13.
随着石油勘探工业的持续发展和技术水平的日益提高,低孔低渗油气藏已成为我国油气勘探开发的重要领域之一,但是该类油气藏的储层岩石物理关系复杂,对其评价也相对较难。本文针对中深层气藏低孔低渗储层评价存在的困难,采用修正的White气体包裹体模型开展岩石物理研究。首先对气体包裹体模型进行分析,得到纵横波速度的计算公式,进而理论计算并分析纵横波速度与孔隙度、饱和度、压力及温度参数的变化关系,最后结合岩心实验对理论计算结果进行了验证。研究结果表明基于气体包裹体模型的理论计算结果与实测数据吻合较好,可以较好地为低孔低渗复杂储层声波速度测量及解释提供技术支持。   相似文献   

14.
致密砂岩气储层的岩石物理模型研究   总被引：3，自引：1，他引：2       下载免费PDF全文
王大兴 《地球物理学报》2016,59(12):4603-4622
根据鄂尔多斯盆地苏里格气田以往实测和新测的共17口井51块岩样超声波实验数据,得到304组不同孔隙度和不同含水饱和度下对应的纵横波速度、泊松比等弹性参数.重新优选计算体积模量和泊松比与含气饱和度的关系,表明苏里格气田上古生界二叠系石盒子组盒8致密砂岩储层的模型与Brie模型(e=2)相似度最高.由此建立的苏里格气田储层岩石物理模型,更好的表征了致密岩石储层物理参数随含气饱和度变化规律,为该区储层预测提供了理论依据.致密储层岩石物理模型研究成果应用于苏里格气田多波地震资料气水预测中,实际例子表明该模型适用于该区的储层和含气性预测,并取得了较好的效果.  相似文献   

15.
The phase velocity of seismic waves varies with the propagation frequency, and thus frequency-dependent phenomena appear when CO2 gas is injected into a reservoir. By dynamically considering these phenomena with reservoir conditions it is thus feasible to extract the frequency-dependent velocity factor with the aim of monitoring changes in the reservoir both before and after CO2 injection. In the paper, we derive a quantitative expression for the frequency-dependent factor based on the Robinson seismic convolution model. In addition, an inversion equation with a frequency-dependent velocity factor is constructed, and a procedure is implemented using the following four processing steps: decomposition of the spectrum by generalized S transform, wavelet extraction of cross-well seismic traces, spectrum equalization processing, and an extraction method for frequency-dependent velocity factor based on the damped least-square algorithm. An attenuation layered model is then established based on changes in the Q value of the viscoelastic medium, and spectra of migration profiles from forward modeling are obtained and analyzed. Frequency-dependent factors are extracted and compared, and the effectiveness of the method is then verified using a synthetic data. The frequency-dependent velocity factor is finally applied to target processing and oil displacement monitoring based on real seismic data obtained before and after CO2 injection in the G89 well block within Shengli oilfield. Profiles and slices of the frequency-dependent factor determine its ability to indicate differences in CO2 flooding, and the predicting results are highly consistent with those of practical investigations within the well block.  相似文献   

16.
The injection of CO2 at the Ketzin pilot site commenced in June 2008 and was terminated in August 2013 after 67 kT had been injected into a saline formation at a depth of 630–650 m. As part of the site monitoring program, four 3D surface seismic surveys have been acquired to date, one baseline and three repeats, of which two were conducted during the injection period, and one during the post‐injection phase. The surveys have provided the most comprehensive images of the spreading CO2 plume within the reservoir layer. Both petrophysical experiments on core samples from the Ketzin reservoir and spectral decomposition of the 3D time‐lapse seismic data show that the reservoir pore pressure change due to CO2 injection has a rather minor impact on the seismic amplitudes. Therefore, the observed amplitude anomaly is interpreted to be mainly due to CO2 saturation. In this study, amplitude versus offset analysis has been applied to investigate the amplitude versus offset response from the top of the sandstone reservoir during the injection and post‐injection phases, and utilize it to obtain a more quantitative assessment of the CO2 gaseous saturation changes. Based on the amplitude versus offset modelling, a prominent decrease in the intercept values imaged at the top of the reservoir around the injection well is indeed associated solely with the CO2 saturation increase. Any change in the gradient values, which would, in case it was positive, be the only signature induced by the reservoir pressure variations, has not been observed. The amplitude versus offset intercept change is, therefore, entirely ascribed to CO2 saturation and used for its quantitative assessment. The estimated CO2 saturation values around the injection area in the range of 40%–60% are similar to those obtained earlier from pulsed neutron‐gamma logging. The highest values of 80% are found in the second seismic repeat in close vicinity to the injection and observation wells.  相似文献   

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

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

19.
Understanding multiphase transport within saline aquifers is necessary for safe and efficient CO2 sequestration. To that end, numerous full‐physics codes exist for rigorously modeling multiphase flow within porous and permeable rock formations. High‐fidelity simulation with such codes is data‐ and computation‐intensive, and may not be suitable for screening‐level calculations. Alternatively, under conditions of vertical equilibrium, a class of sharp‐interface models result in simplified relationships that can be solved with limited computing resources and geologic/fluidic data. In this study, the sharp‐interface model of Nordbotten and Celia (2006a,2006b) is evaluated against results from a commercial full‐physics simulator for a semi‐confined system with vertical permeability heterogeneity. In general, significant differences were observed between the simulator and the sharp‐interface model results. A variety of adjustments were made to the sharp‐interface model including modifications to the fluid saturation and effective viscosity in the two‐phase region behind the CO2‐brine interface. These adaptations significantly improved the predictive ability of the sharp interface model while maintaining overall tractability.  相似文献   

20.
傅磊  刘四新 《地球物理学报》2016,59(12):4464-4472
本文提出了一种初至纵波(P波)与瑞雷面波的交叉梯度联合反演策略.通过对初至P波进行全波形反演可以获得近地表P波速度结构;通过对仅含瑞雷面波信息的地震数据转换到频率-波数域进行加窗振幅波形反演(Windowed-Amplitude Waveform Inversion,w-AWI)可获得近地表横波(S波)速度结构.在二者反演的目标函数中均加入P波速度和S波速度的交叉梯度作为正则化约束项,使得在反演过程中P波速度和S波速度相互制约,相互约束,从而实现对地震初至P波与瑞雷面波的联合反演.数值模拟结果表明交叉梯度联合反演可以提高S波速度反演分辨率,而P波速度反演结果并没有得到提高.实际资料的反演结果表明,交叉梯度联合反演能够获得更加可信的近地表速度结构.  相似文献   

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

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