面向油藏开发地质问题的精细储层反演研究

垦71开发区块的主要产油层系为新近系的馆陶组河流相和古近系的东营组三角洲相地层,具有被断层复杂化的多油水系统。该区井网密集,储层厚度薄,横向关系复杂,储层与非储层速度差异小,后续开发方案调整所要求的精细储层描述难度较大。本文展示了以精细储层描述为目标的地震约束反演实例,利用多井地震、地质标定,结合三维构造解释成果,建立三维多断层精细约束地质模型。在此基础上,采用随机地震反演方法,进行了波阻抗、岩性等多种地层参数反演,较好地解决了该区薄互层储层精细描述问题,证实了精细储层反演技术在油田开发中的重要作用。     

Seismic driven probabilistic classification of reservoir facies for static reservoir modelling: a case history in the Barents Sea

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.     

鄂北上古生界隐蔽气藏地震储层预测技术

鄂尔多斯盆地北部上古生界陆相致密碎屑岩隐蔽气藏具有低孔、低渗、致密、薄互层、多层叠置和非均质性强的特点.其勘探开发一直是世界性难题.针对该区储层特点,展开了以三维地震为核心的储层综合地球物理预测技术研究,开发应用了储层岩石物理分析、地震属性优化、随机模拟反演、地震吸收衰减、地震频谱分析、地震波正演模拟、地震相划分和AVO等多项技术,形成了一整套适合该地区油气地质特征的方法技术系列,有效提高了储层预测的精度,提高了勘探开发井位部署的成功率.     

Integrating facies-based Bayesian inversion and supervised machine learning for petro-facies characterization in the Snadd Formation of the Goliat Field,south-western Barents Sea

Seismic petro-facies characterization in low net-to-gross reservoirs with poor reservoir properties such as the Snadd Formation in the Goliat field requires a multidisciplinary approach. This is especially important when the elastic properties of the desired petro-facies significantly overlap. Pore fluid corrected endmember sand and shale depth trends have been used to generate stochastic forward models for different lithology and fluid combinations in order to assess the degree of separation of different petro-facies. Subsequently, a spectral decomposition and blending of selected frequency volumes reveal some seismic fluvial geomorphological features. We then jointly inverted for impedance and facies within a Bayesian framework using facies-dependent rock physics depth trends as input. The results from the inversion are then integrated into a supervised machine learning neural network for effective porosity discrimination. Probability density functions derived from stochastic forward modelling of endmember depth trends show a decreasing seismic fluid discrimination with depth. Spectral decomposition and blending of selected frequencies reveal a dominant NNE trend compared to the regional SE–NW pro-gradational trend, and a local E–W trend potentially related to fault activity at branches of the Troms-Finnmark Fault Complex. The facies-based inversion captures the main reservoir facies within the limits of the seismic bandwidth. Meanwhile the effective porosity predictions from the multilayer feed forward neural network are consistent with the inverted facies model, and can be used to qualitatively highlight the cleanest regions within the inverted facies model. A combination of facies-based inversion and neural network improves the seismic reservoir delineation of the Snadd Formation in the Goliat Field.     

叠前反演和地震吸收技术在复杂天然气藏地震预测中的应用

针对某复杂断块天然气目标储层,在岩石物理分析的指导下,综合利用地质、地震、测井等资料,提出了一套面向复杂天然气藏的叠前地震预测技术.首先基于地震岩石物理分析得到的初始横波信息,采用叠前贝叶斯非线性三参数反演得到了井旁控制点处精确纵横波速度和密度信息,然后通过叠前/叠后联合反演技术实现了面向目标的弹性阻抗体反演及含气储层敏感参数直接提取,最后结合小波变换时频谱分析的方法从叠前地震资料中估算地层吸收参数值,提高天然气藏识别精度.实际应用表明,综合各种叠前地震预测技术,可以大大提高对复杂天然气藏的识别精度,降低勘探风险.     

A study and application of high-resolution methods for reef reservoir identification

Reefs represent a special type of carbonate trap that plays a key role in the migration, accumulation, and formation of a reservoir. They have commonly been the targets of exploration and development. However, reefs have complex interior structures and easily grow as thin, interbedded geological frames with reef microfacies that include the cap, core, and base of the reef. Because of the inherent drawbacks of seismic signals, including their low frequencies and narrow bandwidths, it is difficult to accurately identify reef reservoirs. Fortunately, the seismic frequency, phase, energy, waveform and other dynamic and geometrical properties can be used to compensate for the energy, expand the frequency bandwidth, and decompose and reconstruct the wavelet to obtain high-resolution seismic data. These data can highlight certain seismic responses of reefs, including boundary reflections, dome-shaped reflections from the reef outline, strong reflections from the reef cap, reflections from the reef bottom, and onlap reflections from the reef flanks. Some impedance response regularities, such as the lower impedance of the reef cap relative to the reef core and biodetritus beach and the fluctuating impedance of the reef-flat complex, are observed by combining log data with geological and high-resolution seismic data for a reef reservoir inversion. These methods were applied to the Changxing Formation in the Yuanba Gas Field. Good prediction results were obtained with a high consistency between the log and seismic data in a comparative analysis with the original seismic data and well logs.     

多属性融合技术在苏14井区的应用

In this study area the geological conditions are complicated and the effective sandstone is very heterogeneous. The sandstones are thin and lateral and vertical variations are large. We introduce multi-attribute fusion technology based on pre-stack seismic data, pre-stack P- and S-wave inversion results, and post-stack attributes. This method not only can keep the fluid information contained in pre-stack seismic data but also make use of the high SNR characteristics of post-stack data. First, we use a one-step recursive method to get the optimal attribute combination from a number of attributes. Second, we use a probabilistic neural network method to train the nonlinear relationship between log curves and seismic attributes and then use the trained samples to find the natural gamma ray distribution in the Su-14 well block and improve the resolution of seismic data. Finally, we predict the effective reservoir distribution in the Su-14 well block.     

基于统计岩石物理模型的各向异性页岩储层参数反演

提出了各向异性页岩储层统计岩石物理反演方法.通过统计岩石物理模型建立储层物性参数与弹性参数的定量关系,使用测井数据及井中岩石物理反演结果作为先验信息,将地震阻抗数据定量解释为储层物性参数、各向异性参数的空间分布.反演过程在贝叶斯框架下求得储层参数的后验概率密度函数,并从中得到参数的最优估计值及其不确定性的定量描述.在此过程中综合考虑了岩石物理模型对复杂地下介质的描述偏差和地震数据中噪声对反演不确定性的影响.在求取最大后验概率过程中使用模拟退火优化粒子群算法以提高收敛速度和计算准确性.将统计岩石物理技术应用于龙马溪组页岩气储层,得到储层泥质含量、压实指数、孔隙度、裂缝密度等物性,以及各向异性参数的空间分布及相应的不确定性估计,为页岩气储层的定量描述提供依据.     

基于线性预测倒谱系数的地震相分析

本文借鉴语音识别技术中的线性预测倒谱系数(LPCC系数)特征参数提取方法对地震数据进行分解,这种方法的优点是:可以获得将子波和反射系数信息分离的地震语音特征参数,对地质现象边界具有较好的描述能力,使我们可以从不同维度更细致地观察隐藏在地震数据中的地质特征.理论模型分析表明,基于LPCC系数的地震分析具有较高的地震相划分能力.实际地震资料应用表明,LPCC系数对储层特征的描述比常规三瞬属性更为细致,不同阶次LPCC系数在描述储层不同特征时也保持了内在的联系.采用K均值聚类方法对提取的12阶和24阶LPCC系数进行聚类分析,聚类结果与目的层段古地形较为吻合,较好地反映了研究区的断裂、礁滩相带、深水扇和储层的分布特征,说明在地震相分析中采用LPCC系数作为特征参数是可行和有效的.     

地震岩相识别概率表征方法

储层岩相分布信息是油藏表征的重要参数,基于地震资料开展储层岩相识别通常具有较强的不确定性.传统方法仅获取唯一确定的岩相分布信息,无法解析反演结果的不确定性,增加了油藏评价的风险.本文引入基于概率统计的多步骤反演方法开展地震岩相识别,通过在其各个环节建立输入与输出参量的统计关系,然后融合各环节概率统计信息构建地震数据与储层岩相的条件概率关系以反演岩相分布概率信息.与传统方法相比,文中方法通过概率统计关系表征了地震岩相识别各个环节中地球物理响应关系的不确定性,并通过融合各环节概率信息实现了不确定性传递的数值模拟,最终反演的岩相概率信息能够客观准确地反映地震岩相识别结果的不确定性,为油藏评价及储层建模提供了重要参考信息.模型数据和实际资料应用验证了方法的有效性.     

随机反演在储层预测中的应用

针对隐蔽油气藏储层预测的需要,开展了地震反演研究,根据目前的实际应用将储层预测中的基于模型的地震反演分为三个实施阶段:即构造反演、声波波阻抗或弹性波阻抗反演以及岩性反演,并对每个阶段的目的、关键技术及其原理进行了详细描述,尤其是详细描述了基于马尔科夫链的蒙特卡罗随机模拟技术.最后给出了一个综合应用测井、地质、地震资料进行反演,从而进行储层预测的实例.     

基于FFT-MA谱模拟的快速随机反演方法研究

虽然基于地质统计学的随机反演方法能够有效融合测井资料中的高频信息,但计算效率低,占用内存大,限制了它在实际资料中的应用.本文在保留传统随机反演方法优点的基础上,创造性地引入傅里叶滑动平均(Fast Fourier Transform-Moving Average,FFT-MA)谱模拟进行频率域的地质统计模拟,并利用逐步变形算法(Gradual Deformation Method,GDM)确保模拟结果与实际地震数据的匹配,构建了基于FFT-MA谱模拟的新的快速随机反演方法.与常规随机反演相比,新方法不仅分辨率高,而且能够使反演解得到快速收敛,有效提高计算效率,减少内存占用.模型试算获得了与理论模型吻合度较好的高分辨率反演结果.实际资料分析也表明新方法所得到的高分辨率反演结果能够对薄互储层进行良好的展示,为薄储层的识别提供高效可靠的技术支持.     

宽频数据在致密砂岩储层预测中的应用

西湖凹陷深层致密砂岩储层具有良好的勘探开发前景,受埋深影响,目的层地震资料品质较差.A构造通过斜缆宽频采集和处理获取宽频地震数据,提升了资料品质,然而应用常规子波提取方法对宽频数据进行子波提取并反演计算纵波阻抗,结果与井上实测数值差异较大,影响储层的定量解释.针对这一问题,提出统计性子波和确定性子波相结合的长短子波合并宽频子波提取方法,提取的宽频子波比常规子波低频丰富、旁瓣小,能更真实地反映地震信息,约束稀疏脉冲反演的纵波阻抗结果与测井曲线吻合度更高.基于宽频数据和常规数据分别进行约束稀疏脉冲弹性波阻抗反演,预测A构造优质储层分布,经已钻井证实,宽频数据比常规数据储层预测精度高,预测的储层展布特征与研究区地质沉积认识一致.结果表明:这种基于宽频子波提取的宽频资料应用方法有效降低了致密砂岩储层预测的多解性.     

储层弹性与物性参数地震叠前同步反演的确定性优化方法

储层弹性与物性参数可直接应用于储层岩性预测和流体识别,是储层综合评价和油气藏精细描述的基本要素之一.现有的储层弹性与物性参数地震同步反演方法大都基于Gassmann方程,使用地震叠前数据,通过随机优化方法反演储层弹性与物性参数;或基于Wyllie方程,使用地震叠后数据,通过确定性优化方法反演储层弹性与物性参数.本文提出一种基于Gassmann方程、通过确定性优化方法开展储层弹性和物性参数地震叠前反演的方法,该方法利用Gassmann方程建立储层物性参数与叠前地震观测数据之间的联系,在贝叶斯反演框架下以储层弹性与物性参数的联合后验概率为目标函数,通过将目标函数的梯度用泰勒公式展开得到储层弹性与物性参数联合的方程组,其中储层弹性参数对物性参数的梯度用差分形式表示,最后通过共轭梯度算法迭代求解得到储层弹性与物性参数的最优解.理论试算与实际资料反演结果证明了方法的可行性.     

Three‐dimensional modelling of magnetotelluric data to image Sehqanat hydrocarbon reservoir in southwestern Iran

A detailed magnetotelluric survey was conducted in 2013 in the Sehqanat oil field, southwestern Iran to map the geoelectrical structures of the sedimentary Zagros zone, particularly the boundary between the Gachsaran Formation acting as cap rock and the Asmari Formation as the reservoir. According to the electrical well logs, a large resistivity contrast exists between the two formations. The Gachsaran Formation is formed by tens to hundreds of metres of evaporites and it is highly conductive (ca. 1 Ωm–10 Ωm), and the Asmari Formation consists of dense carbonates, which are considerably more resistive (more than 100 Ωm). Broadband magnetotelluric data were collected along five southwest–northeast directed parallel lines with more than 600 stations crossing the main geological trend. Although dimensionality and strike analysis of the magnetotelluric transfer functions showed that overall they satisfied local 2D conditions, there were also strong 3D conditions found in some of the sites. Therefore, in order to obtain a more reliable image of the resistivity distribution in the Sehqanat oil field, in addition to standard 2D inversion, we investigated to what extent 3D inversion of the data was feasible and what improvements in the resistivity image could be obtained. The 2D inversion models using the determinant average of the impedance tensor depict the main resistivity structures well, whereas the estimated 3D model shows significantly more details although problems were encountered in fitting the data with the latter. Both approaches resolved the Gachsaran–Asmari transition from high conductivity to moderate conductivity. The well‐known Sehqanat anticline could also be delineated throughout the 2D and 3D resistivity models as a resistive dome‐shaped body in the middle parts of the magnetotelluric profiles.     

地震相控混沌反演在岩性油气藏储层预测中的应用——以RY坳陷M井区为例

RY坳陷M井区位于浅水三角洲前缘斜坡带,多发育水下分流河道沉积,砂岩储层厚度薄、横向变化快,且同一时期的砂体在平面上表现为多套油水系统,容易形成岩性油气藏,利用常规手段进行储层预测的难度较大。本文通过引入地震相控制,利用纵向及横向上的构造起伏、岩性变化、沉积环境变迁等先验信息降低了反演的多解性;同时将混沌理论与反演算法结合,根据反演迭代的实时状态自适应地修改正则化算子,调节反演稳定性和分辨率达到最佳平衡,既避免了混沌现象的出现,又最大限度地揭示薄层信息。M井区的实际应用显示,相控混沌反演可以突出地震波形的细节变化,准确刻画储层尖灭和分叉点位置,达到精准预测岩性油气藏的目的,对后续岩性油气藏的勘探开发具有较大指导意义。      

无井条件下建立碎屑岩储层地震地质模型研究

常规的储层建模以井数据为基础,建立孔隙度、渗透率和含油饱和度等储层属性参数模型,并通过油田开发生产数据进行拟合,得到最佳地质模型.在南海西北部深水区无井控制的情况下,针对中央峡谷浊积水道储层,采用确定性建模与连续型随机建模相结合的方法,建立碎屑岩储层地震地质模型,应用地震正演模拟和相似性分析方法,确定浊积水道砂体最佳模拟参数.文章首先分析了碎屑岩储层特征,应用连续型随机建模的方法模拟储层空间分布.然后分析实际地震、地质资料,应用层序地层学思想,对研究区划分沉积体系域,建立层序地层格架;基于波阻抗和多种地震属性,采用人机交互的方法对浊积水道储层框架结构进行精细解释,建立储层结构模型;将储层结构模型与砂泥互层随机介质进行"交"、"并"运算,利用连续型随机模拟方法对储层内部介质进行精细刻画,建立中央峡谷浊积水道随机介质储层模型;通过对模型正演模拟、常规处理和相似性分析等建立最佳储层模型.文章研究了在无井条件下建立储层地震地质模型的方法,揭示南海西北部深水区中央峡谷浊积水道储层的空间展布和内部结构特征.     

秦皇岛32-6油田明化镇组河流相储层地球物理参数影响因素研究

秦皇岛32-6油田明化镇组储层为河流相沉积,单砂体规模小,砂体间接触关系复杂,受到地震分辨率限制,对叠置砂体储层地震响应规律认识不明确,导致基于地震信息的储层砂体解释错误.通过数值正演模拟,可以正确认识砂体储层的地震响应规律,提高砂体储层预测的准确率.所以,在进行地震响应正演模拟前,非常有必要对该区储层地球物理参数影响因素展开研究,通过分析研究区14口测井资料,建立了岩石物理参数(孔隙度、泥质含量、合流体性)与地球物理参数(速度、密度、波阻抗等)之间的相关关系,分析了不同沉积微相对地球物理参数的控制作用,并给出了相应的数学表达式,为正确认识河流相储层砂体地震响应规律奠定了准确的地球物理参数基础.     

纵横波弹性阻抗联立反演在GD地区的应用

在GD油田复杂油藏描述中应用叠前纵横波弹性阻抗反演,精确地进行了油藏岩性的划分。利用三个或三个以上部分叠加数据,进行纵横波弹性阻抗联立反演,既克服了因叠后地震反演结果单一而不能满足复杂储层描述的需求,又避免了由于叠前道集信噪比低造成反演结果不稳定的缺陷。本文论述了叠前弹性波阻抗反演的基本原理,结合GD地区实际资料,对反演过程中涉及的角道集子波提取、层位标定、横波速度预测、弹性参数提取与解释等关键步骤进行了详细研究,指出基于测井资料分析的多种弹性参数综合解释是提高叠前地震反演应用效果的关键。