基于波动方程正演模拟分析薄砂层含不同流体的AVO特征

塔河油田碎屑岩储层厚度较小,薄层AVO特征受薄层内多次波反射以及薄层调谐效应和频率滤波效应影响,导致薄层AVO特征与厚层AVO特征差异较大。为此,笔者设计了不同含油气性薄层砂体模型来进行波动方程正演模拟,得到了CRP道集记录,提取了目的层反射振幅,基于入射波振幅计算了反射系数,从而获得了薄层地震正演AVO曲线,分析了塔河油田碎屑岩储层含不同流体时的AVO特征,找出了实际地质条件下的含油气性与AVO属性特征的对应关系,以指导地震剖面上的含流体预测工作。     

地震多属性分析在混积环境砂体识别中的应用——以桑塔木南地区石炭系Ⅳ砂组油为例

地震多属性分析是储层预测的重要手段之一。以塔河油田桑塔木南地区石炭系Ⅳ砂组油为例,在单井沉积相分析的基础上,优选振幅属性信息,进行AVF分频反演,建立振幅属性与砂体之间的耦合关系及砂体识别模式,刻画沉积相平面分布。将地震反射结构类型细分为4种类型;以Ⅳ砂组油层顶面为基准,上下浮动10 ms进行平均波峰振幅属性提取,对灰质含量重的西南地区进行AVF分频反演;振幅属性与砂体之间总体存在正相关关系,砂含量越高,振幅值越大。水下分流河道主要分布在S1072井区及研究区西北缘;砂坪主要分布在AT10及S111井区,呈片状展布;潮砂脊主要分布在南部,可作为研究区岩性圈闭的重点目标区。      

基于振幅谱互相关的薄层厚度计算

反射波振幅和峰值频率两种属性当前被广泛应用于预测薄层砂体厚度。然而,实际地震资料的振幅通常不保幅会降低振幅计算薄层厚度的可靠性;峰值频率法计算薄层厚度具有泰勒展开近似前提,会降低其计算薄层厚度的精度。为了提高薄层厚度估计的精度,提出了一种利用振幅谱互相关函数计算薄层厚度的新方法。首先,构建目的层的振幅谱与子波振幅谱和反射系数褶积后振幅谱之间的相关系数函数,其最大值所对应的地层厚度为目的层砂体的厚度;然后,利用高斯迭代方法求取相关系数函数最大值,从而计算出薄层厚度。通过楔状体模型测试证实新方法是可行、有效的。相较于峰值频率方法,新方法能够更精确地计算薄层厚度。实际资料的应用进一步证实新方法适用于薄层厚度计算。     

河流相砂泥岩薄互层预测方法研究与应用

由于河道变迁频繁,河流相储层单砂体厚度薄、横向变化快,利用井网预测砂体分布的方法无法满足老油田中后期剩余油挖潜的需要。针对大港港东油田河流相储层发育的特点,建立三维地质模型进行储层预测方法实验,探索预测单期河道砂体空间展布的方法,实验表明:通过离散合成地震记录的方式,分别分析每套砂体地震反射对地震记录能量的贡献,发现在其中一套砂体地震反射能量的零点提取的振幅切片能够反映另一套单砂体的空间展布形态,利用此方法可达到预测单砂体的目的。零值法提取振幅切片能够有效地刻画井间砂体的横向变化特征,为厘定港东油田的油水关系,调整井网部署提供有力的依据。     

薄储层段地震响应特征分析——以大庆西部斜坡区某研究区块为例

正研究区位于松辽盆地西部斜坡区,主要目的层为萨尔图油层,其萨二油层组为主要储层发育位置。对萨二层段提取地震振幅类属性,发现强振福区域并不能反映河道砂体发育位置。通过对区内振幅强弱与储层发育位置进行研究后,认为利用地震振幅属性不能完全将砂泥区分。因此从地震波波形的角度出发,分析地震波波波形与薄储层发育情况之间存在对应关系,最终发现通波形聚类分析的办法可以对区内储层发育区域作出了预测。1区域地质概况     

塔河油田叠后地震资料高保真处理技术

塔河油田岩溶储集体(洞穴、礁滩、孔洞、裂缝或裂隙等)的地震波场特征是弱反射背景之上的强弱不同的绕射波,因此,绕射波准确归位成像和振幅保持是研究岩溶储集体的基础。笔者从塔河油田地震资料高保真处理实践中,总结出了一套适应该区碳酸盐岩储层预测要求的以叠后高精度偏移算法和振幅保持处理为主的叠后高保真处理技术,取得显著的效果。     

塔河油田奥陶系缝洞型碳酸盐岩储层预测研究

针对塔河油田奥陶系碳酸盐岩缝洞型储层的结构特点,以地质、钻井、测井及生产测试等资料分析为出发点,以各类储层地震响应为基础,分析了有利储层的地震属性特征、地震反射结构及波形相关特点,最终将有利储层综合划分为两类,并建立了各自的地震综合预测模式:Ⅰ类储层表现为串珠状、条带状振幅变化率异常和残丘高点+杂乱反射特征,相位角变化大于20°;Ⅱ类储层为弱或强振幅、高振幅变化率、弱相干、低速度和多种反射结构。提高了塔河油田及其周缘的钻井成功率。     

地震分频解释技术在河道砂预测中的应用

地震资料分频处理技术是一种比常规地震属性识别储层更有效的方法。该方法将一定时窗内地震反射波通过离散傅里叶变从时间域换转换到频率域,用转换后的振幅谱和相位谱识别薄储层的厚度变化和边界。通过该技术可以对地震资料中所有单频信号进行有针对性的分析,已经证明是一种具有较高可信度的技术。针对胜利油田老河口油田老168井区馆陶组河道砂油藏厚度薄、横向变化大的地质特征,本文运用分频处理技术对该井区馆陶组地层砂体分布进行预测,并对储层非均质性进行了分析,经钻井证实,预测结果与实际情况十分吻合。     

薄互层干涉对叠前AVO属性的影响分析

常规AVO分析是以单一界面为前提,用来分析薄层的AVO响应特征存在误差。Brekhovski给出的层状传播矩阵方程考虑了层内波的反射和透射以及层厚等影响因素,更适合用于讨论薄互层的地震响应特征。对薄层进行Brekhovski正演,发现不同入射角下振幅与厚度关系不同,与叠后全叠加振幅相比,截距P属性和薄层厚度之间的线性关系更强,说明P属性比全叠加振幅更有利于薄层厚度预测。不同薄互层模型的Brekhovski正演结果表明:净厚度和泊松比是影响AVO曲线特征的两个主要因素,净厚度越小,泊松比对截距P和梯度G属性的影响越小;当泊松比一定时,净厚度和P、G属性呈近似线性关系,P属性与净厚度关系更敏感,可用于薄互层净厚度预测。应用实例证明在薄互层区P属性用于砂地比预测效果比较好,叠后振幅和叠前反演的弹性阻抗预测的净砂岩厚度精度较差。     

塔河油田卡拉沙依组砂组沉积相与储层研究

卡拉沙依组砂岩是塔河油田石炭系重要储层，从上到下分为4个砂组，即Ⅰ砂组、Ⅱ砂组、Ⅲ砂组及Ⅳ砂组，根据岩石组合特征、岩石粒度、分选性、以及沉积构造等特征，可得出Ⅰ砂组为较典型的扇三角洲沉积环境；Ⅱ-Ⅳ砂组为潮坪环境，其中Ⅲ砂组沉积时的古环境水动力相对最强，以潮间带下部为主，储层砂体最发育，单个砂体的厚度也最大，为纵向上最好储层段。     

数字地震台网地震数据库管理系统

为了能较好地满足台网编目室日常工作的需要 ,我们开发了数字台网地震数据库管理系统。该系统满足了编目室日常编报工作的需要 ,实现了对地震信息的数据库管理 ,为用户提供了在矢量电子地图上进行任意多边形地震信息查询的功能。     

地震反演——煤田地震勘探的新进展

地震反演根据数据来源、信息类型、反演结果及实现方法可划分不同种类,地震反演的基本方法可分为基于波动方程的反演和以地震褶积模型为基础的反演二大类。目前国内建立在地震记录褶积模型基础的主要反演软件有STRATA、JASON、ISIS。依托国家重大产业技术开发专项“西部煤炭资源高精度三维地震勘探技术”,在阳泉二矿7个钻孔中．应用STRATA和JASON软件对其三维地震资料进行波阻抗反演,其结果不仅可提高纵向分辨率,准确界定煤层顶、底板,提高下组弱反射煤层的连续性,而且通过波阻抗切片,可获得煤层及顶板的岩性信息。另外通过对比STRATA和JASON软件,发现前者的分辨率较高,而后者在波形连续性表现较好。     

合成地震记录在煤炭地震勘探中的应用

人工 合 成 地 震 记 录 是 利 用 已 知 测 井 资 料 ,通 过 数 学 方 法 模 拟 实 际 地 震 记 录 的 一 项 技 术 ,其 主 要 用 途 是 :鉴别 反 射 波,识 别多 次 波 ,建 立 地震 模 型 。该 技 术 在 铁法 矿 区 地 震勘 探 中 的 应用 ,优化 了 地 震 勘探 设 计 、提 高 了 采 集 及解 释 精 度、准 确标 定 地 质 层位 。     

地震危险性、地震危害性和地震易损性

做好城市防震减灾工作 ,关键是对未来可能遭遇地震灾害定量化预测。本文论述地震灾害定量化的三要素 :地震危险性 ,地震危害性和地震易损性 ,它们概念完全不同 ,且很容易混淆 ,但又存在因与果的关系。本文还简要介绍厦门市地震科技工作者开展的“闽南地区综合防震减灾示范工程”。     

隐伏活断层地震勘探正演模拟方法

浅层高分辨率地震勘探是城市活断层探测手段中最有效、最可靠的方法之一,但是由于隐伏活断层勘探时存在着勘探深度浅、地层的属性差异小以及环境噪声大等因素,导致地震勘探数据信噪比低、反射信号弱等问题,严重影响地震数据的采集、处理和解释。采用一阶应力-速度方程的高阶交错网格有限差分方法, 结合完全匹配层(PML)的边界处理方法和通量校正(FCT)法消除频散的技术手段,对活断层模型地震勘探进行正演模拟,清晰直观地表现了断层处的细微差别,有助于地震勘探中对信号的分析与理解。结果表明,这些处理手段结合波动方程有限差分算法,能够有效实现薄层、弱反射的隐伏活断层模型的波场正演模拟,可有效识别隐伏活断层的断层倾角、断距、规模、跨越范围、断层深度以及断裂带影响范围,为分析这些参数对地震记录的影响提供有力帮助,同时结合实测地震数据的处理、解释,为防灾减灾提供具体有力的判据。     

Seismic Hazard of Honduras

In this paper we have described the proceduresused, input data applied and results achieved in ourefforts to develop seismic hazard maps of Honduras.The probabilistic methodology of Cornell is employed.Numerical calculations were carried out by making useof the computer code SEISRISK III. To examine theimpact of uncertainties in seismic and structuralcharacteristics, the logic tree formalism has beenused. We compiled a de-clustered earthquake cataloguefor the region comprising 1919 earthquakes occurringduring the period from 1963 to 1997. Unified momentmagnitudes were introduced. Definition of aseismotectonic model of the whole region under review,based on geologic, tectonic and seismic information,led to the definition of seven seismogenetic zones forwhich seismic characteristics were determined. Fourdifferent attenuation models were considered. Resultsare expressed in a series of maps of expected PGA for60% and 90% probabilities of nonexceedence in a50-year interval which corresponds to return periodsof 100 and 475 years, respectively. The highest PGAvalues of about 0.4g (90% probability ofnon-exceedence) are expected along the borders withGuatemala and El Salvador.     

Seismic hazard of Egypt

Earthquake hazard parameters such as maximum expected magnitude,M _max, annual activity rate,, andb value of the Gutenberg-Richter relation have been evaluated for two regions of Egypt. The applied maximum likelihood method permits the combination of both historical and instrumental data. The catalogue used covers earthquakes with magnitude 3 from the time interval 320–1987. The uncertainties in magnitude estimates and threshold of completeness were taken into account. The hazard parameter determination is performed for two study areas. The first area, Gulf of Suez, has higher seismicity level than the second, all other active zones in Egypt.b-values of 1.2 ± 0.1 and 1.0 ± 0.1 are obtained for the two areas, respectively. The number of annually expected earthquakes with magnitude 3 is much larger in the Gulf of Suez, 39 ± 2 than in the other areas, 6.1 ± 0.5. The maximum expected magnitude is calculated to be 6.5 ± 0.4 for a time span of 209 years for the Gulf of Suez and 6.1 ± 0.3 for a time span of 1667 years for the remaining active areas in Egypt. Respective periods of 10 and 20 years were reported for earthquakes of magnitude 5.0 for the two subareas.     

Seismic Zoning of Croatia

The territory of Croatia and neighboring regions is divided into 17 seismic source zones, considering available seismological and geological data. On this basis, seismic hazard elements (seismicity rate, maximum magnitude, b-value, probabilities of exceedance and return periods for a predefined set of magnitudes) are computed using the maximum likelihood method appropriate for treating data-sets with variable completeness thresholds. The values of long term expected peak horizontal acceleration obtained by using a combination of the deterministic and the probabilistic procedure are the highest in the Dubrovnik zone, while the Zagreb zone has the highest earthquake hazard in the continental part of the country. This revised version was published online in June 2006 with corrections to the Cover Date.     

Seismic hazard of Bulgaria

A simplified seismotectonic model is proposed for Bulgaria by introducing generalized seismogenic areas containing systems of complex geometry faults. A tectonic scheme, which considers the main faults only, is then derived from this. The assessment of the regional seismic hazard is done using different approaches: the Gumbel, the Cornell, and the fault rupture model methods. A series of relations among seismological parameters are derived from the available data. The results obtained by the different approaches are similar: the most dangerous area is the Struma zone, located in southwestern Bulgaria.