共查询到20条相似文献，搜索用时 859 毫秒
1.
Jean‐Paul van Gestel Kevin D. Best Olav I. Barkved Jan H. Kommedal 《Geophysical Prospecting》2011,59(4):673-681
The frequent time‐lapse observations from the life of field seismic system across the Valhall field provide a wealth of information. The responses from the production and injection wells can be observed through time‐shift and amplitude changes. These observations can be compared to modelled synthetic seismic responses from a reservoir simulation model of the Valhall Field. The observed differences between the observations and the modelling are used to update and improve the history match of the reservoir model. The uncertainty of the resulting model is reduced and a more confident prediction of future reservoir performance is provided. A workflow is presented to convert the reservoir model to a synthetic seismic response and compare the results to the observed time‐lapse responses for any time range and area of interest. Correlation based match quality factors are calculated to quantify the visual differences. This match quality factor allows us to quantitatively compare alternative reservoir models to help identify the parameters that best match the seismic observations. Three different case studies are shown where this workflow has helped to reduce the uncertainty range associated with specific reservoir parameters. By updating various reservoir model parameters we have been able to improve the match to the observations and thereby improve the overall reservoir model predictability. The examples show positive results in a range of different reservoir modelling issues, which indicates the flexibility of this workflow and the ability to have an impact in most reservoir modelling challenges. 相似文献
2.
3.
早期的褶积微分算子都是基于正反傅立叶变换而实现的,其精度比四阶有限差分的精度稍高,本文将计算数学中的Forsyte广义正交多项式微分算子与褶积算子相结合,构建了一个新的快速、高精度褶积微分算子,其计算结果非常接近实验函数微分的精确值,精度与16阶有限差分的精度相当,远优于错格伪谱法的精确度.另外,2.5维数值模拟比二维模拟可以更真实地模拟三维介质的臬个剖面的波场,并且2.5维地震波模拟的计算量比三维模拟的计算量及计算耗时要大大减少.本文利用基于Forsyte广义正交多项式褶积微分算子法计算2.5维非均匀介质地震波场,模拟结果表明,该算法的计算速度快,计算精度高,能够直观、高效地反映复杂介质中波场的传播规律,并且2.5维波场数值模拟具有更高的计算效率,是一种非常值得深入研究并广泛应用的方法. 相似文献
4.
基于共聚焦点技术的叠前AVP（振幅随射线参数变化）分析与常规叠后反演方法相比优势明显,但传统通过褶积和互相关运算来实现的方法依赖于聚焦算子,而在复杂构造区走时计算困难且子波难以精确提取,从而导致了聚焦算子不准确,而且褶积和互相关运算会影响信噪比和分辨率,基于此,本文提出了基于保真振幅单程波延拓算法的叠前AVP成像方法.该方法利用保真振幅傅里叶有限差分延拓算法实现两步聚焦,分别生成共聚焦点道集和网格点道集,既充分利用了保真振幅延拓算法在振幅保持方面的优势,也可以发挥傅里叶有限差分方法对复杂构造区横向变速适应性强的优势,而且两步聚焦过程都不需要聚焦算子,从而解决了传统方法中走时计算和子波提取的问题.模型试算结果表明了方法的正确性和可行性,而针对实际地震资料的试处理结果与传统方法相比具有更高的信噪比和分辨率,表明了方法的有效性.该方法为复杂构造区油气检测提供了一种新的地球物理依据. 相似文献
5.
地震波传播有限差分模拟的人工边界问题 总被引：6，自引：3，他引：6
在用有限差分法模拟地震波传播中，差分网格人工边界的反射问题是至关的。本文简要介绍了几种消除人工边界反射的方法及其主要优缺点，联立这些不同的吸收边界条件于实际的有限差分计算可以给出好的结果。 相似文献
6.
A magnetotelluric finite-element modelling algorithm is developed, which is capable of handling three-dimensional conductive and magnetic anisotropic anomalies. Different from earlier three-dimensional magnetotelluric anisotropic modelling methods, the algorithm we presented has taken the magnetic anisotropy into consideration. The variational equations are produced by the Galerkin method and the governing equations are solved using a hexahedral vector edge finite-element method. The accuracy of this algorithm is firstly validated by comparing its solutions with the results of finite-difference method for a three-dimensional conductive arbitrary anisotropic model, and then validated by comparing with analytical solutions for a one-dimensional magnetic model. The responses of four kinds of models under different conditions are studied, and some conclusions are obtained. It shows that for materials with a high magnetic permeability, its influence on magnetotelluric responses cannot be ignored in some circumstances. Especially, if the magnetic susceptibility is exceptionally high, it may really distort the apparent resistivities of lower resistive anomalies. These conclusions are also beneficial for magnetotelluric survey. 相似文献
7.
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. 相似文献
8.
Many of the numerical techniques used for seismic zonation studies treat one-dimensional structural models and/or the incidence of plane polarized body waves. These techniques are often not adequate for laterally heterogeneous structures and for sources that are not located beneath the site of interest. In such cases a more rigorous treatment of the combined effects of the source, the path and the site response is needed. This can be accomplished with a hybrid approach combining modal summation and the finite-difference technique. To demonstrate the differences between these techniques, the ground motion in the city of Benevento (Italy) is modelled. We first compare the results obtained with one-and two-dimensional structural models for vertical incidence of plane polarized body waves. These results are then compared with those obtained with the hybrid approach for two-dimensional structural models.The comparisons have allowed us to find important differences in the response obtained with the different modelling techniques. For the same site, these differences consist of strong variations in amplitude and in the shape of the spectral amplifications. For a seismic source which is not located beneath the site, vertical incidence of waves significantly overestimates the local hazard in a laterally homogeneous structure. For a laterally heterogeneous area, we can conclude that one-dimensional modelling fails to estimate the seismic hazard, whereas for a seismic source which is not located beneath the site of interest, two-dimensional modelling with vertical incidence of plane polarized body waves may not allow reliable estimates to be made of the frequency bands at which amplifications occur. The results obtained for two-dimensional structural models are used for a zonation of the city of Benevento. 相似文献
9.
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. 相似文献
10.
The effects of source and receiver motion on seismic data are considered using extensions of the standard convolutional model. In particular, receiver motion introduces a time-variant spatial shift into data, while source motion converts the effect of the source signature from a single-channel convolution in time to a multichannel convolution in time and space. These results are consistent with classical Doppler theory and suggest that Doppler shifting can introduce distortions into seismic data even at relatively slow acquisition speeds. It is shown that, while both source and receiver motion are known to be important for marine vibroseis acquisition, receiver motion alone can produce significant artifacts in marine 3D data. Fortunately, the convolutional nature of the distortions renders them amenable to correction using simple deconvolution techniques. Specifically, the effects of receiver motion can be neutralized by applying an appropriate reverse time-variant spatial shift, while those due to source motion can be addressed by introducing time-variant spatial shifts both before and after standard, deterministic, signature deconvolution or correlation. 相似文献
11.
In geophysical inverse problems, the posterior model can be analytically assessed only in case of linear forward operators, Gaussian, Gaussian mixture, or generalized Gaussian prior models, continuous model properties, and Gaussian-distributed noise contaminating the observed data. For this reason, one of the major challenges of seismic inversion is to derive reliable uncertainty appraisals in cases of complex prior models, non-linear forward operators and mixed discrete-continuous model parameters. We present two amplitude versus angle inversion strategies for the joint estimation of elastic properties and litho-fluid facies from pre-stack seismic data in case of non-parametric mixture prior distributions and non-linear forward modellings. The first strategy is a two-dimensional target-oriented inversion that inverts the amplitude versus angle responses of the target reflections by adopting the single-interface full Zoeppritz equations. The second is an interval-oriented approach that inverts the pre-stack seismic responses along a given time interval using a one-dimensional convolutional forward modelling still based on the Zoeppritz equations. In both approaches, the model vector includes the facies sequence and the elastic properties of P-wave velocity, S-wave velocity and density. The distribution of the elastic properties at each common-mid-point location (for the target-oriented approach) or at each time-sample position (for the time-interval approach) is assumed to be multimodal with as many modes as the number of litho-fluid facies considered. In this context, an analytical expression of the posterior model is no more available. For this reason, we adopt a Markov chain Monte Carlo algorithm to numerically evaluate the posterior uncertainties. With the aim of speeding up the convergence of the probabilistic sampling, we adopt a specific recipe that includes multiple chains, a parallel tempering strategy, a delayed rejection updating scheme and hybridizes the standard Metropolis–Hasting algorithm with the more advanced differential evolution Markov chain method. For the lack of available field seismic data, we validate the two implemented algorithms by inverting synthetic seismic data derived on the basis of realistic subsurface models and actual well log data. The two approaches are also benchmarked against two analytical inversion approaches that assume Gaussian-mixture-distributed elastic parameters. The final predictions and the convergence analysis of the two implemented methods proved that our approaches retrieve reliable estimations and accurate uncertainties quantifications with a reasonable computational effort. 相似文献
12.
Maryam Sadat Mirkamali Abdolrahim Javaherian Hossein Hassani Mohammad Reza Saberi Seyed Abolfazl Hosseini 《Geophysical Prospecting》2020,68(7):2195-2216
Quality, availability and consistency of the measured and interpreted well log data are essential in the seismic reservoir characterization methods, and seismic petrophysics is the recommended workflow to achieve data consistency between logs and seismic domains. This paper uses seismic petrophysics workflow to improve well logs and pore geometry interpretations for an oil carbonate reservoir in the Fahliyan Formation in the southwest of Iran. The petrophysical interpreted well logs, rock physics and well-to-seismic tie analysis are integrated into the proposed workflow. Our implementation incorporates revising petrophysical well log interpretations and updating pore geometry characteristics to obtain a better well-tie quality. We first propose an improved pore-type characterization approach based on both P- and S-wave velocities for quantifying pore geometry. Then, seismic logs are estimated accordingly, and the results are used in the well-to-seismic analysis. The quality of the well-tie is improved, furthermore, by iterating on the petrophysical interpreted well logs as well as the calculated pore geometries. For the intervals with high-quality data, our workflow improves the consistency between the results of measured and modelled seismic logs. For the intervals with problematic well logs, the application of our proposed workflow results in the successful replacement of the poor data and subsequently leads to an improved wavelet estimation and well-tie results. In both cases, a higher quantification of pore geometries is achieved, which in turn is confirmed by the core images and formation micro-imager analysis. 相似文献
13.
传统的完全匹配层技术是一种能够较为有效地消除边界反射的边界条件,但是当表层为泊松比较高的自由表面时,该技术可能会产生不稳定的现象.针对传统的完全匹配层技术固有的不稳定和掠射情况下吸收效果不佳等缺陷,发展了多轴完全匹配层、卷积完全匹配层以及将两者结合的多轴卷积完全匹配层等3种边界条件.本文介绍了水平自由表面的不同处理方法以及传统、多轴、卷积和多轴卷积等4种完全匹配层条件的原理,通过二维半无限空间模型的交错网格有限差分正演模拟对比,分析了几种自由边界实施方法在这几种完全匹配层条件下的稳定性,并通过提取单道波形与解析解进行对比,定性分析了水平自由表面几种不同处理方法的准确性以及各自的适用条件. 结果表明,泊松比和水平自由表面实施方法对波场模拟效果及其稳定性有重要影响. 相似文献
14.
L. Alvarez G. F. Panza F. Vaccari B. E. González 《Pure and Applied Geophysics》2001,158(9-10):1763-1782
— We present the results of complete P-SV and SH waves modelling, up to a maximum frequency of 1?Hz, along two profiles in Santiago de Cuba city. The seismic sources are located in the depth range from 10 to 40?km on the Oriente fault zone at distances of several tens of kilometres from the city. The calculation has been made by a hybrid method: Modal summation in the regional anelastic model (one-dimensional) where the source is buried, and finite differences in the local sedimentary anelastic models (two-dimensional). The analysis of the influence of the depth and of the distance of the source on the site effects shows that standard traditional methods, based on the deconvolution analysis of the rock outcrop motion, can lead to erroneous results. 相似文献
15.
传统的高阶有限差分波动方程数值模拟方法采用高阶差分算子近似空间偏导数,能有效抑制空间频散.然而,传统的有限差分法仅采用二阶差分算子近似时间偏导数,这使得地震波场沿时间外推的精度较低.当采用较大的时间采样间隔,传统的有限差分法模拟波场会出现明显的时间频散,甚至不稳定.本文基于新的差分结构和中心网格剖分,发展了一种空间任意偶数阶精度、时间四阶和六阶精度的时空域有限差分方法.基于对离散后的频散关系进行泰勒展开,本文推导了时空域高阶有限差分算子的差分系数.相速度分析表明时间四阶、六阶精度的差分方法能显著地减小传统时间二阶精度差分方法的时间频散.在相同的精度下与传统差分法比较,本文发展的时间四阶、六阶有限差分方法的计算效率比传统方法高.均匀和非匀均介质中的波场数值模拟实验进一步证实本文研究的时空高阶有限差分方法的优越性. 相似文献
16.
We explore the link between basin modelling and seismic inversion by applying different rock physics models. This study uses the E‐Dragon II data in the Gulf of Mexico. To investigate the impact of different rock physics models on the link between basin modelling and seismic inversion, we first model relationships between seismic velocities and both (1) porosity and (2) effective stress for well‐log data using published rock physics models. Then, we build 1D basin models to predict seismic velocities derived from basin modelling with different rock physics models, in a comparison with average sonic velocities measured in the wells. Finally, we examine how basin modelling outputs can be used to aid seismic inversion by providing constraints for the background low‐frequency model. For this, we run different scenarios of inverting near angle partial stack seismic data into elastic impedances to test the impact of the background model on the quality of the inversion results. The results of the study suggest that the link between basin modelling and seismic technology is a two‐way interaction in terms of potential applications, and the key to refine it is establishing a rock physics models that properly describes changes in seismic signatures reflecting changes in rock properties. 相似文献
17.
为了研究松南地区陆相砂泥岩薄互层储层岩石物性横向变化、微构造等对油气聚集的影响和探索利用高分辨率井间地震技术直接指示油气分布的可能性,开展了井间地震及地面小三维地震、VSP、测井等立体地震观测．使用了自行研制的井中炸药震源和常规地震记录仪器．采用层析成像方法、反射资料叠前偏移成像方法等对采集到的资料进行研究,得到了不同分辨率的地下构造和岩石物性信息,展示出不同方法技术在分辨薄储层能力上的差异．高分辨率的井间地震成果较好地实现了井间地层对比,查清了两井间储层横向变化和油气分布情况．层析速度与构造信息等的综合应用解释了两井油气产能差异的原因,建立了油气地质模型．所得成果显示出井间地震技术在油气开发中的潜力． 相似文献
18.
Selecting the correct resolution in distributed hydrological modelling at the watershed scale is essential in reducing scale-related errors. The work presented herein uses information content (entropy) to identify the resolution which captures the essential variability, at the watershed scale, of the infiltration parameters in the Green and Ampt infiltration equation. A soil map of the Little Washita watershed in south-west Oklahoma, USA was used to investigate the effects of grid cell resolution on the distributed modelling of infiltration. Soil-derived parameters and infiltration exhibit decreased entropy as resolutions become coarser. This is reflected in a decrease in the maximum entropy value for the reclassified/derived parameters vis a vis the original data. Moreover, the entropy curve, when plotted against resolution, shows two distinct segments: a constant section where no entropy was lost with decreasing resolution and another part which is characterized by a sharp decrease in entropy after a critical resolution of 1209 m is reached. This methodology offers a technique for assessing the largest cell size that captures the spatial variability of infiltration parameters for a particular basin. A geographical information system (GIS) based rainfall-runoff model is used to simulate storm hydrographs using infiltration parameter maps at different resolutions as inputs. Model results up to the critical resolution are reproducible and errors are small. However, at resolutions beyond the critical resolution the results are erratic with large errors. A major finding of this study is that a large resolution (1209 m for this basin) yields reproducible model results. When modelling a river basin using a distributed model, the resolution (grid cell size) can drastically affect the model results and calibration. The error structure attributable to grid cell resolution using entropy as a spatial variability measure is shown. 相似文献
19.
地震储层预测是油气勘探的重要组成部分,但完成该项工作往往需要经历多个环节,而多工序或长周期的研究分析降低了勘探效率.基于油气藏分布规律及其在地震响应上所具有的特点,本文引入卷积神经网络深度学习方法,用于智能提取、分类并识别地震油气特征.卷积神经网络所具有的强适用性、强泛化能力,使之可以在小样本条件下,对未解释地震数据体进行全局优化提取特征并加以分类,即利用有限的已知含油气井段信息构建卷积核,以地震数据为驱动,借助卷积神经网络提取、识别蕴藏其中的地震油气特征.将本方案应用于模型数据及实际数据的验算,取得了预期效果.通过与实际钻井信息及基于多波地震数据机器学习所预测结果对比,本方案利用实际数据所演算结果与实际情况有较高的吻合度.表明本方案具有一定的可行性,为缩短相关环节的周期提供了一种新的途径. 相似文献
20.
非稳态地震稀疏约束反褶积研究（英文） 总被引：1，自引：1，他引：0
传统Robinson褶积模型主要受缚于三种不合理的假设,即白噪反射系数、最小相位地震子波与稳态假设,而现代反射系数反演方法（如稀疏约束反褶积等）均在前两个假设上寻求突破的同时却忽视了一个重要事实：实际地震信号具有典型的非稳态特征,这直接冲击着反射系数反演中地震子波不随时间变化的这一基础性假设。本文首先通过实际反射系数测试证实,非稳态效应造成重要信息无法得到有效展现,且对深层影响尤为严重。为校正非稳态影响,本文从描述非稳态方面具有普适性的非稳态褶积模型出发,借助对数域的衰减曲线指导检测非稳态影响并以此实现对非稳态均衡与校正。与常规不同,本文利用对数域Gabor反褶积仅移除非稳态影响,而将分离震源子波和反射系数的任务交给具有更符合实际条件的稀疏约束反褶积处理,因此结合两种反褶积技术即可有效解决非稳态特征影响,又能避免反射系数和地震子波理想化假设的不利影响。海上地震资料的应用实际表明,校正非稳态影响有助于恢复更丰富的反射系数信息,使得与地质沉积和构造相关的细节特征得到更加清晰的展现。 相似文献