塔河奥陶系薄层状孔洞储集体定量预测方法研究

本文以塔河油田非"串珠"薄层状储集体作为研究对象,结合钻井、测井资料,指出薄层状储集体储集空间主要为溶蚀孔洞,主要位于T47不整合面下0~60 m.正演模拟结果证实T47不整合面低频强反射屏蔽了下部薄层状孔洞储集体弱信号,利用多子波地震道分解技术剥离了T47旁瓣强能量,重构数据恢复了薄层状溶蚀孔洞储集体弱信号.采用多子波分解和多属性线性回归融合技术建立了薄层状储集体厚度定量预测方法,降低了单一属性预测的多解性,提高了薄层状储集体厚度预测的准确性.塔河十区东地区应用结果表明,薄层状孔洞储集体预测厚度与实钻厚度的相对误差小于10%,说明该方法取得了较好的应用效果.     

辽河油田薄层砂体的预测与追踪

辽河盆地及其外围盆地在漫长的地质历史时期形成了浊积扇、河流相、扇三角洲等多种类型的砂体,单层砂体较薄.为了准确地部署井位,应用井约束反演,地震属性分析和分频分析等技术,对不同类型的薄层砂体及含气砂层进行了预测,较准确地刻画了薄砂层的分布及形态,为部署井位提供了依据.经钻探,预测比较准确,效果良好.     

塔北隆起玉东地区盐下薄砂岩高分辨地震预测

古近系薄砂岩是塔里木盆地玉东地区盐下主力产气储层段,薄砂岩储层的准确预测是玉东地区勘探成功的关键.本文通过岩心、测井、地震等资料综合研究及储层岩石物理分析建立了薄砂岩正演地质模型,正演模拟与地震响应特征分析发现玉东三维8种敏感地震属性与薄砂岩厚度呈4类(线性正相关、线性负相关、指数负相关、对数正相关)强相关,其中交流分量、道积分与薄层指示因子等敏感属性经RGB-HIS信息融合具有最佳的薄砂岩分辨能力.经地震属性信息融合、地质统计学高分辨反演和地震沉积学预测,古近系薄砂岩储层砂体主要分布于玉东北部的YT5井至YD1井、YD4井一带,薄砂岩发育三期透镜状水下分流河道和砂坝叠加砂体,砂体相互叠置连通.研究形成了一套薄砂岩储层预测新方法,为塔里木盆地古近系盐下薄砂岩储层预测和岩性圈闭勘探提供技术指导.     

场地砂土地震液化的模糊聚类分析

鉴于场地砂土地震液化的概念以及某些主要影响因子都是不确切的,具有中间过渡性质,无明确的内涵和外延。故文中提出了一个新的分析途径,即模糊聚类分圻的方法,来预测现场砂层的液化。根据地震的基本数据和土层的性状,应用模糊聚类的分析方法将50个场地砂层地震液化的历史资料划分为5类,作为预测场地砂层液化的模糊模式。这种分析方法除具有研究模糊问题的功用外,还县有在分析时既不需要任何假定,也不引入任何试验误差等优点。 文中最后基于所建立的模糊模式,预测了90个场地液化的历史资料,表明所建议的方法是有效而可行的。     

淤砂层对重力坝地震响应的影响分析

在大坝地震响应分析中如何考虑淤砂层,以及淤砂层对大坝的地震响应影响如何,是目前混凝土坝抗震分析和设计中需要进一步研究的问题。将库水-淤砂层-坝体-坝基动力相互作用当做波动散射问题,采用广义饱和多孔介质统一计算框架分析库水(流体)、淤砂层(饱和多孔介质)、坝体与坝基(固体)间的耦合。研究了一种考虑库水可压缩性与淤砂层影响的库水-淤砂层-坝体-坝基动力相互作用的高效分析方法,并通过自编程序实现了该方法。以Koyna重力坝为对象,以脉冲P波和SV波为输入,给出了库水-坝体-坝基体系不同时刻的波场快照,分析了其波场特征。设计了不同淤砂层厚度的4种工况,对比分析了4种工况下坝体的位移、加速度和最大应力,研究结果表明：库底淤砂层对坝体位移影响很小,但对坝体各点加速度和最大应力有减小作用,并且随着淤砂层厚度的增加,作用会增强。     

对称型隔夹层与单砂层时频特征对比分析

隔夹层作为单砂层与薄互层、薄互层与薄互层之间的过渡型地层,对不同薄互层单元具有重要的分隔作用.厘清单砂层与隔夹层之间时频特征的区别与联系,对后续的薄互层单元划分与预测具有重要意义.本文基于波动理论,利用深度域相移法对不同厚度的单砂层、对称型隔夹层进行正演模拟.采用广义S变换分析零偏移距地震道反射复合波的瞬时频谱,提取单砂层与对称型隔夹层时域、频域属性参数进行对比分析.研究发现:当对称型隔夹层顶部砂岩与中部泥岩厚度之和与顶底反射系数同号的单砂层厚度相同时,二者频谱特征极为相似,时域波形存在镜像关系.当对称型隔夹层顶部砂岩与中部泥岩厚度之和与顶底反射系数异号的单砂层厚度相同时,二者频谱特征存在明显差异,这为砂泥岩剖面中单砂层与对称型隔夹层的定性区分奠定了理论基础.此外,对称型隔夹层的第一陷频频率与顶部砂岩和中部泥岩厚度之和具有近似的反比关系,这为对称型隔夹层厚度预测提供了重要技术手段.     

碳酸盐岩溶洞厚度定量预测方法研究与应用

针对地震反演薄储层的厚度预测多解性强,且难以适应碳酸盐岩溶洞储层地层横向突变的状况,难以满足该类储层定量描述和储量计算的要求等问题.在综合考虑到溶洞型储层顶底地震反射对偶响应的特点、以及深入分析了储层厚度变化与合成地震记录波形和能量的变化特征,而提出了一种储层地震响应特征建立与厚度定量预测的方法即为建立储层厚度谱.该方法主要关键技术是:地震子波精确提取、一组不同储层厚度变化与合成地震响应模型库的建立、在镶边余弦平方衰减梯度函数约束下构建地震波形和能量特征双重约束的储层厚度谱、预测误差分析以及储层厚度三维定量预测.该方法的特点是充分利用储层厚度变化引起的波形和能量的差异,不需要井约束和层位约束,克服了部分人为干扰,减少了预测多解性.预测结果的垂向分辨率突破了传统的1/4波长的限制,一般可达1/20～1/40波长,对于最难预测的薄互层也能达到1/8波长.实际资料表明该方法对于非均质储层等特殊储层预测有较好的推广和应用价值.     

倾斜地层地震液化和滑移的有限元分析

应用饱和多孔介质动力学分析倾斜地层的土壤动力线性反应、液化和液化滑移问题，地下水位上的地层简化为单相介质层，饱和夹砂层看作是两相介质，水是可压缩的，采用双曲线非线性本构关系，考虑了砂土的剪胀性、刚度退化、滞回特性和土水相对运动等因素。基于土力学模型，建立了适用于分析非自由场地液化的动力方程组，基于是否考虑发生渗流问题，同时建立了两种离散形式：一种是以土骨架位移和水位移为未知量的矩阵方程，另一种是以土骨架位移、水位移和孔隙水压力为未知量的矩阵方程，初步分析了适用于多孔介质波动模拟的离散模型的人工边界问题，形成的方法将有助于问题的解决。     

渤海西南部新近系超薄储层定量预测技术研究与应用

渤海西南部新近系明化镇组主要为极浅水三角洲和曲流河沉积,受特定沉积环境的控制,发育大量厚度小于λ/8(λ为视波长)的超薄储层,利用常规地球物理手段难以进行定量预测.为定量预测此类超薄储层,本文通过对多项地球物理技术改造升级与优化组合,建立了一套完整的超薄储层定量预测技术,具体包括两项关键子技术.(1)基于薄储层的地震高分辨率处理技术.首先优选能较好继承地震有效信息的部分地震道进行叠加,然后开展基于ARMA的谱蓝化处理,提高地震资料分辨率.(2)基于薄储层的解释性优化处理技术.在地震资料高分辨率处理基础上进行基于匹配追踪算法的子波重构处理,将薄储层的地震反射信息进行小幅度放大,然后进行倾角中值滤波和倾角扩散滤波处理,有效压制数据中的随机噪声和突出储层边界,使得地震资料可解释性增强.技术成果应用于渤海A油田,成效显著,同时取得两点主要认识:(1)方法可识别厚度小于λ/8的超薄储层,在渤海西南部海域取得了成功应用,对陆相断陷盆地浅层超薄储层预测具有较好借鉴意义;(2)方法针对薄储层进行地震优化处理,为凸显薄储层地震响应,一定程度破坏了适中厚度和厚储层的地震响应,不能做到薄、厚储层同时预测.     

利用测井资料评价有利沉积相带及其储量分布

针对雁木西油田白垩系低阻油层受沉积相带及其断裂、构造和岩性等多种因素控制影响,测井信息应用难度增大的特点,提出利用多种测井响应提取单渗砂层能量厚度信息,使其曲线形态、特征和数据大小能够集中反映相对高渗的单渗砂层最大沉积能量及厚度变化,控制划分扇三角洲前缘水下辫状河道、浅水河道骨架砂体有利微相带的发育、规模及分布范围,克...     

基于频率域峰值属性的河道砂体定量预测及应用(英文)

河道砂体是陆相含油气盆地最重要的储集类型之一,其边界识别和厚度定量预测是储层预测的热点难题。本文在总结现有方法技术的基础上,提出一种利用频率域峰值属性进行河道砂体边界识别和厚度定量预测的新方法。对典型河道薄砂体地震反射进行了正演模拟,构造了一种新的地震属性——峰值频率-振幅比,研究表明:峰值频率属性对地层厚度变化敏感,振幅属性对地层岩性变化敏感,两者比值突出河道砂体的边界,同时,借助峰值频率与薄层厚度间存在的定量关系进行薄砂体厚度计算。实际数据应用表明,地震峰值频率属性可以较好的刻画河道的平面展布特征;峰值频率-振幅比属性可以提高对河道砂体边界的识别能力;利用频率域地震属性进行砂体边界识别及厚度定量预测是可行的。     

等厚薄互层时频特征的正演模拟

本文基于波动理论,利用深度域相移法对不同互层数、不同单层厚度的等厚薄互层进行正演模拟.采用广义S变换分析零偏移距地震道反射复合波的瞬时频谱.研究发现当单层厚度大于3/16波长时,等厚薄互层时域特征为中、高频等幅振荡,瞬时频谱为与单层厚度成反比的单一峰值梳状函数;当单层厚度小于3/16波长时,震源子波的有限带宽对瞬时频谱的高频频率峰有压制作用,导致等厚薄互层中部的时域振幅和频域最大幅度均为低值,薄互层的结构特征不易分辨,说明宽高频的震源子波是提高薄互层分辨能力的关键.此外,理论推导和实验分析均证明:当地层厚度大于1/8波长时,其峰值频率与薄互层单层厚度存在定量解析关系,这为等厚薄互层单层厚度定量预测提供了重要的技术手段.     

地震资料预测薄层厚度方法研究与应用

薄层厚度预测一直是公认的难题之一,其难度就在于如何准确地识别和提取薄层的地震属性.常规方法是利用时间域或频率域地震属性与薄层厚度的线性关系计算.但是理论与实际资料表明,不同的薄层和地层组合对地震波的动力学的信息影响很大,各种参数与薄层厚度成非线性关系,使用单一的信息不可能准确预测薄层厚度.本文利用三种线性预测原理（模型）,经数学变换为属性参数,采用非线性BP网络预测薄层厚度,取得了令人满意的效果.     

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

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

Crust-mantle transitional zone of Tianshan orogenic belt and Junggar Basin and its geodynamic implication

The traveling time of the reflection waves of each shot point from the crust-mantle transitional zone has been obtained by data processing using wavelet transform to the waves reflected from the crust-mantle transitional zone. The crust-mantle transitional zone of the Xayar-Burjing geoscience transect can be divided into three sections: the northern margin of the Tarim Basin, the Tianshan orogenic belt and Junggar Basin. The crust-mantle transitional zone is composed mainly of first-order discontinuity in the Tarim Basin and the Junggar Basin, but in the Tianshan orogenic belt, it is composed of 7–8 thin layers which are 2-3 km in thickness and high and low alternatively in velocity, with a total thickness of about 20km. The discovery of the crust-mantle transitional zone of the Tianshan orogenic belt and Junggar Basin and their differences in tectonic features provide evidence for the creation of the geodynamic model “lithospheric subduction with intrusion layers in crust” for the Tianshan orogenic belt.     

薄层时频特征的正演模拟

基于正演模拟的薄层时频特征响应分析是薄层定性识别与厚度定量预测的基础与关键.本文基于波场延拓理论,利用深度域相移法对不同厚度薄层进行了正演模拟.采用广义S变换对零偏移距地震道的反射复合波进行瞬时频谱分析,发现当薄层顶底反射系数极性相反时,复合波的干涉作用具有升频降幅作用,而极性相同时具有升幅降频作用;理论推导和实验分析均证明主峰值频率与薄层双程旅行时间厚度存在定量解析关系,这为薄层厚度定量预测提供了重要的技术手段;由于峰值频率对反射系数大小与极性变化不敏感,因此,峰值频率方法与时域振幅方法相比更具稳定性.     

饱和砂土透镜体液化对建筑物地震反应的影响

采用一种能分析有结构物存在的场地地震液化问题的二维有效应力有限元分析方法,研究饱和砂土透镜体液化对建筑物地震反应的影响。计算中采用了更为合理的迭代方式处理土的非线性,考虑了Ｋｃ对孔压的作用,引入了透射边界。取建筑物为短周期结构。考虑了透镜化宽度、厚度、埋深以及输入地震动类型幅值对结构加速度反应的影响。计算结果表明：（１）所采用的方法与已有的模型实验结果有很好的对应关系,可用于招考莪存在下的场地砂土     

基于反射系数谱理论的薄层多波AVO

AVO分析是目前地震勘探潜在油气储层的一个重要方法。体散射信息包含了层结构、岩性和孔隙流体信息,对地震勘探非常有用。但是基于 Zoeppritz方程的传统AVO分析之只包含了单层信息。薄层厚度定量解释对构造解释、储层描述和储层横向预测都非常重要。本文阐述的基于频率域弹性传播矩阵反射系数谱方法既考虑了层界面引起的振幅变化(Zoeppritz方程),也考虑了层内传播引起的振幅变化。因此该反射系数谱既包括单一层界面信息,也包括层内体散射信息。该反射系数谱是层厚和频率的连续函数,便于分析频率和层厚对反射系数谱的影响。可分析的薄层厚度可以无限小,直至消失。可分析的频率是任意的和连续的。这是对时间域反射系数做傅里叶变换无法实现的。地震波的传播是复杂的,各种波型是同时存在而且相互转换的,该反射系数谱考虑了各种波型在传播过程中的相互转换以及多次波。与比射线方法比更便于正演薄层多波多分量AVO响应。     

Effect of anelastic patchy saturated sand layers on the reflection and transmission responses of a periodically layered medium

Based on analytic relations, we compute the reflection and transmission responses of a periodically layered medium with a stack of elastic shales and partially saturated sands. The sand layers are considered anelastic (using patchy saturation theory) or elastic (with effective velocity). Using the patchy saturation theory, we introduce a velocity dispersion due to mesoscale attenuation in the sand layer. This intrinsic anelasticity is creating frequency dependence, which is added to the one coming from the layering (macroscale). We choose several configurations of the periodically layered medium to enhance more or less the effect of anelasticity. The worst case to see the effect of intrinsic anelasticity is obtained with low dispersion in the sand layer, strong contrast between shales and sands, and a low value of the net‐to‐gross ratio (sand proportion divided by the sand + shale proportion), whereas the best case is constituted by high dispersion, weak contrast, and high net‐to‐gross ratio. We then compare the results to show which dispersion effect is dominating in reflection and transmission responses. In frequency domain, the influence of the intrinsic anelasticity is not negligible compared with the layering effect. Even if the main resonance patterns are the same, the resonance peaks for anelastic cases are shifted towards high frequencies and have a slightly lower amplitude than for elastic cases. These observations are more emphasized when we combine all effects and when the net‐to‐gross ratio increases, whereas the differences between anelastic and elastic results are less affected by the level of intrinsic dispersion and by the contrast between the layers. In the time domain, the amplitude of the responses is significantly lower when we consider intrinsic anelastic layers. Even if the phase response has the same features for elastic and anelastic cases, the anelastic model responses are clearly more attenuated than the elastic ones. We conclude that the frequency dependence due to the layering is not always dominating the responses. The frequency dependence coming from intrinsic visco‐elastic phenomena affects the amplitude of the responses in the frequency and time domains. Considering intrinsic attenuation and velocity dispersion of some layers should be analyzed while looking at seismic and log data in thin layered reservoirs.     

Side lobes of wavelets impact identification of thin sand bodies

Recently, we found that side lobes of wavelets have a large impact on the identification of thin sand reservoirs when studying some gas fields in a basin in Northwest China. Reflections from the top of the H Formation, in which there are gas-bearing thin sand bodies, have the main wavelet lobe between two weak peak side lobes. The lower one always mixes with another peak reflected from the top of a thin sand reservoir. That makes it difficult to identify the sand reservoir. In order to solve this, many forward models were set up using typical well logs. 2D synthetic profiles were produced using Ricker wavelets to study the relationships between the effects of wavelet side lobes and thin sand position and frequency and between amplitude and the thin sand body. We developed the following conclusions： First, it is easier to identify thin sands in a shallower position. Second, a good way to tell sand body reflections from side lobes is by comparing profiles with different frequency windows. Third, it is helpful and effective to describe sand extent using amplitude attributes.