CN107688200B - 针对低幅度构造畸变的速度校正方法 - Google Patents

针对低幅度构造畸变的速度校正方法 Download PDF

Info

Publication number
CN107688200B
CN107688200B CN201710237118.XA CN201710237118A CN107688200B CN 107688200 B CN107688200 B CN 107688200B CN 201710237118 A CN201710237118 A CN 201710237118A CN 107688200 B CN107688200 B CN 107688200B
Authority
CN
China
Prior art keywords
velocity
speed
sandstone
well
correction
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710237118.XA
Other languages
English (en)
Other versions
CN107688200A (zh
Inventor
刘玉娟
李梅
邱津
张志让
张梦萦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
New JCM Group Co.,Ltd.
Original Assignee
Landocean Group Ltd By Share Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Landocean Group Ltd By Share Ltd filed Critical Landocean Group Ltd By Share Ltd
Priority to CN201710237118.XA priority Critical patent/CN107688200B/zh
Publication of CN107688200A publication Critical patent/CN107688200A/zh
Application granted granted Critical
Publication of CN107688200B publication Critical patent/CN107688200B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy

Abstract

本发明提出了一种针对低幅度构造畸变的速度校正方法,包括:根据井上精确的时深标定关系和处理得到的叠加速度构建初始的空变速度场;通过对井上实测速度曲线的分析,通过敏感地震属性的优选勾勒高速含钙砂岩异常体的平面分布边界;进行纯砂岩段速度和含钙后整体岩性速度的统计;拟合构建高速含钙砂岩校正量随时间的变化关系的速度校正关系式;利用拟合的速度校正关系式得到相应的速度校正量体,与初始的空变速度场相结合得到异常体范围内的校正后的速度场;对初始的空变速度场进行相应异常体范围内的速度替换,形成高速含钙砂岩异常速度校正后的空变速度场。本发明解决了因特殊速度体存在导致的低幅度构造特征不真实的现象。

Description

针对低幅度构造畸变的速度校正方法
技术领域
本发明涉及HSE管理技术领域,特别涉及一种针对低幅度构造畸变的速度校正方法。
背景技术
低幅度构造,即在地震资料上表现为反射同相轴平直而变化幅度很小,一般在几米到几十米左右。在低幅度构造主体部位上实施的两口钻井A和B,对同一目的层顶部,在时间剖面上A井的时间小于B井(即,理论上该目的层顶部在A井的实钻深度要浅于B井),如果在A井附近出现局部速度异常(砂岩储层内出现含钙砂岩,含钙砂岩速度高于砂岩速度),将造成下伏构造形态产生畸变,出现B井实钻深度要小于A井的反转现象,这种速度异常引起的构造畸变尤其对低幅度构造的影响是非常大的。
基于地震剖面解释得到的构造图属于时间域,新钻井用到的构造图实为深度域,这就需要构建精细的速度场将时间域构造转成深度域构造。准确的速度求取,一直是勘探和开发的核心问题。最常用的速度场构建信息来源于井上时深标定关系以及地震处理得到的层速度。
如果速度异常体分布范围较广,在常规的地震速度分析中就能识别,能够通过地震资料得出大致的速度异常体的平面分布范围。但往往这种特殊的速度异常体在空间上局部发育,呈现横向上延展范围小,纵向上厚度薄的特点。受限于地震速度的分辨率,这种异常体是很难在地震资料处理过程中识别出来。与地震资料相比,井上速度信息分辨率高,能够较好的指示出井点处速度异常体的纵向分布,但缺乏速度异常体平面展布的认识。因此可看出,地震资料和井资料在识别异常体上各有优缺点,单凭任一资料都很难相对准确的勾勒出速度异常体的空间展布。
井震结合构建速度场技术,即以井上标定的速度为硬数据控制,以地震处理得到的平均速度(层速度经DIX公式转换)为软数据作为横向趋势约束,该技术较传统的单纯基于井上速度和基于地震速度的构建方法在准确度上有了较大的提高。该方法在地下构造相对简单、地层横向速度变化不大的情况下能够很好的满足勘探和开发的需求,但在发育特殊地质体而导致的速度横向变化,尤其对低幅度构造地区,该方法仍具有一定的缺陷。因为,在实际地震资料处理中,拾取的速度谱在横、纵向上一般间隔较大,难以反映叠加速度的横向变化细节。即便是采用缩小速度谱间隔的方法,尽管横向分辨率有所提高,但缺乏针对性,也会产生大量的工作量。
发明内容
本发明的目的旨在至少解决所述技术缺陷之一。
为此,本发明的目的在于提出一种针对低幅度构造畸变的速度校正方法,解决因特殊速度体存在导致的低幅度构造特征不真实的现象。
为了实现上述目的,本发明的实施例提供一种针对低幅度构造畸变的速度校正方法,包括如下步骤:
步骤S101,根据井上精确的时深标定关系和处理得到的叠加速度构建初始的空变速度场,其中,所述初始的空变速度场包括:井上高纵向分辨精度和地震横向分辨率;
步骤S102,通过对井上实测速度曲线的分析,确定高速含钙砂岩异常体在纵向上发育的层段及个数,获取含钙砂岩体对应的地震反射特征,并通过敏感地震属性勾勒高速含钙砂岩异常体的平面分布边界;
步骤S103,以井上的实测速度曲线为分析对象,针对整套目的层选取多个含钙砂岩发育段,进行纯砂岩段速度和含钙后整体岩性速度的统计;
步骤S104,根据砂岩段速度和含钙后整体岩性速度的统计结果,拟合构建高速含钙砂岩校正量随时间的变化关系的速度校正关系式;
步骤S105,在高速含钙砂岩平面分布范围内,利用拟合的速度校正关系式得到相应的速度校正量体,与所述步骤S101中初始的空变速度场相结合得到异常体范围内的校正后的速度场;
步骤S106,对所述步骤S101中初始的空变速度场进行相应异常体范围内的速度替换,形成高速含钙砂岩异常速度校正后的空变速度场。
进一步,在所述步骤S101中,在构建初始的空变速度场之前,进行井震时深标定,以从中选取曲线质量和标定效果优的井,用于初步速度场构建的纵向约束;
从地震处理输出的速度场需要经过DIX公式转成层速度及平均速度,以为所述空变速度场的构建提供横向变化趋势。
进一步,在所述步骤S102中,
选取纯波地震数据体,在构造地层格架的约束下结合井震标定结果分析速度异常体对应的地震反射特征。
进一步,在所述步骤S103中,对井上实测速度曲线进行环境校正以及归一化,以确保除地质因素外,不同井对应大致深度的相似岩性段具有接近的速度范围。
进一步,在所述步骤S106中,替换速度场具有与原始速度场相同的时间范围、数据体采样率以及近似的量纲范围。
根据本发明实施例的针对低幅度构造畸变的速度校正方法,解决了因特殊速度体存在导致的低幅度构造特征不真实的现象。在井震综合构建的初步速度场的基础上,针对性的从井上特殊地质体的岩性和速度出发,结合地震反射特征进行其平面分布范围的刻画以及纵向上速度校正量的建立,经特殊地质体速度校正后的速度场在横向及纵向分辨率精度上进一步得到了提高,具有快速、高效的特点。
本发明附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1为根据本发明实施例的针对低幅度构造畸变的速度校正方法的流程图;
图2为根据本发明实施例的基于井和地震信息的初始空变速度场的示意图;
图3为根据本发明实施例的井上高速异常体的识别及相应的地震反射特征标定的示意图;
图4为根据本发明实施例的研究区速度校正样本选取及速度校正量随时间变化关系的构建示意图;
图5为根据本发明实施例的校正前后速度对比图;
图6为根据本发明实施例的目标区速度校正前后构造对比图。
具体实施方式
下面详细描述本发明的实施例,实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
如图1所示,本发明实施例的针对低幅度构造畸变的速度校正方法,包括如下步骤:
步骤S101,根据井上精确的时深标定关系和处理得到的叠加速度构建初始的空变速度场。
其中,初始的空变速度场包括:井上高纵向分辨精度和地震横向分辨率。
在本步骤中,在构建初始的空变速度场之前,要进行精细的井震时深标定,以便从中选取曲线质量好、标定效果高的井,用于初步速度场构建的纵向约束。
从地震处理输出的速度场需要经过DIX公式转成层速度及平均速度,以为空变速度场的构建提供横向变化趋势。
步骤S102,通过对井上实测速度曲线的分析,确定高速含钙砂岩异常体在纵向上发育的层段及个数,在井震资料标定的基础上,在层序地层解释格架的约束下,获取含钙砂岩体对应的地震反射特征,并通过敏感地震属性勾勒高速含钙砂岩异常体的平面分布边界。
具体地,选取的地震数据体必须是纯波地震数据体,在构造地层格架的约束下结合井震标定结果,从振幅、能量、相位以及频率等方面要充分分析速度异常体对应的地震反射特征,挖掘异常体典型的响应特征,避免与其他岩性混淆。
步骤S103,以井上的实测速度曲线为分析对象,针对整套目的层,从含钙砂岩的厚度、深度等方面出发,选取多个含钙砂岩发育段,进行纯砂岩段速度和含钙后整体岩性速度的统计。
在本步骤之前,还需要对井上实测速度曲线进行环境校正以及归一化,以确保除地质因素外,不同井对应大致深度的相似岩性段具有接近的速度范围。
步骤S104,根据砂岩段速度和含钙后整体岩性速度的统计结果,形成相应的包含目的层深度、时间值、纯砂岩速度、含钙后整体砂岩速度、样本的时间厚度、样本的深度厚度以及速度校正量,拟合构建高速含钙砂岩校正量随时间的变化关系的速度校正关系式,作为研究区速度场精细空间异常体速度校正的依据。
在本步骤中,高速含钙砂岩样本的选取要具有典型性,充分考虑到不同含钙砂岩厚度、不同深度以及不同含钙砂岩与纯砂岩组合对速度变化的影响,样本选取要尽量均匀分布在目的层段。
步骤S105,在高速含钙砂岩平面分布范围内,利用拟合的速度校正关系式得到相应的速度校正量体,与步骤S101中初始的空变速度场相结合得到异常体范围内的校正后的速度场。
通过对研究区内不同井选取的样本拟合得到的速度校正量随时间的变化关系要具有较高的相关系数,提高速度校正的准确度。
步骤S106,对步骤S101中初始的空变速度场进行相应异常体范围内的速度替换,形成高速含钙砂岩异常速度校正后的空变速度场,包括:形成最终的融合井上纵向分辨率、地震横向变化趋势以及特殊异常体的高精度速度场。
其中,替换速度场具有与原始速度场相同的时间范围、数据体采样率以及近似的量纲范围,从而保证校正后的速度横、纵向变化合理。
下面参考图2至图6对本发明实施例的针对低幅度构造畸变的速度校正方法进行说明。
(1)首先把处理提供的叠加速度做DIX转换,转换成平均速度,再利用井的时深速度对平均速度体进行校正,构建初始的空变速度场。该速度场纵向上具有高精度的井上速度信息,横向上反映了速度的变化趋势,如图2所示。
(2)参考图3,从层速度对比看,5口井整体的速度低频背景趋势比较一致,但Well2井存在若干薄的高速层,如900m-1200m,1300m-1700m深度范围。经井震标定得出,浅层高速含钙砂岩段与地震剖面局部反射轴上拉相对应,并呈杂乱、中弱反射,越靠近速度变化层,地震反射特征畸变越大(箭头所指)。基于地震反射特征的识别刻画异常速度体的平面分布范围。
(3)以发育速度异常体的井为基础,从浅层到深层,选取不同厚度的高速含钙砂岩和其相近砂岩的速度为样本,利用井上时深关系得出样本点的时间深度。根据含钙砂岩和不含钙砂岩的速度相比得到该时间深度对应的速度校正量,如图4所示。拟合关系式是速度校正量与高速含钙砂岩厚度以及时间深度的综合关系表征,可看出速度校正量随时间深度的增加而递减,可作为地震速度谱校正的有利依据。
(4)黑色边界为基于地震剖面及属性识别的高速含钙砂岩体分布范围,对于这一层段速度校正量大致分布在1.0342-1.0360之间。如图5所示,通过校正前后的速度对比可看出,黑色边界内以及下方的速度分布趋势有所不同,校正后的速度数值变大,分布趋势横向更合理。
(5)图6为根据本发明实施例的目标区速度校正前后构造对比图。针对目标区,两次成图的构造形态有很大改变,含钙砂岩校正后的构造特征表现为构造压缩,长轴为近北西向,而非近南北向,且构造高点向主控断层走向方向延伸,与区域构造特征和后续新钻井相吻合。
根据本发明实施例的针对低幅度构造畸变的速度校正方法,解决了因特殊速度体存在导致的低幅度构造特征不真实的现象。在井震综合构建的初步速度场的基础上,针对性的从井上特殊地质体的岩性和速度出发,结合地震反射特征进行其平面分布范围的刻画以及纵向上速度校正量的建立,经特殊地质体速度校正后的速度场在横向及纵向分辨率精度上进一步得到了提高,具有快速、高效的特点。
本发明实施例的针对低幅度构造畸变的速度校正方法,具有以下优点:
(1)在利用井上精度和地震横向速度趋势信息构建空变速度场基础上,又综合考虑了井上已知异常速度体与地震反射特征的关系,进而在全区利用敏感地震属性精细的确定了其平面分布范围,弥补了地震速度精度低、难以表征速度异常体的不足。
(2)本发明是一种空变的速度校正方法,通过在目的层段内选取的多个典型样本构建的速度校正量与时间的变化关系不仅仅反映了某一段包含异常体造成的速度变化,也同时体现了浅层异常速度对深层速度分析的累积效应,更精确、更合理的解决了因异常速度体带来的下伏低幅度构造畸变的现象,较真实的还原了下伏的构造形态,降低了后期钻井勘探的风险。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。本发明的范围由所附权利要求及其等同限定。

Claims (5)

1.一种针对低幅度构造畸变的速度校正方法,其特征在于,包括如下步骤:
步骤S101,根据井上精确的时深标定关系和处理得到的叠加速度构建初始的空变速度场,其中,所述初始的空变速度场包括:井上高纵向分辨精度和地震横向分辨率;
步骤S102,通过对井上实测速度曲线的分析,确定高速含钙砂岩异常体在纵向上发育的层段及个数,获取含钙砂岩体对应的地震反射特征,并通过敏感地震属性勾勒高速含钙砂岩异常体的平面分布边界;
步骤S103,以井上的实测速度曲线为分析对象,针对整套目的层选取多个含钙砂岩发育段,进行纯砂岩段速度和含钙后整体岩性速度的统计;
步骤S104,根据砂岩段速度和含钙后整体岩性速度的统计结果,拟合构建高速含钙砂岩校正量随时间的变化关系的速度校正关系式;
步骤S105,在高速含钙砂岩平面分布范围内,利用拟合的速度校正关系式得到相应的速度校正量体,与所述步骤S101中初始的空变速度场相结合得到异常体范围内的校正后的速度场;
步骤S106,对所述步骤S101中初始的空变速度场进行相应异常体范围内的速度替换,形成高速含钙砂岩异常速度校正后的空变速度场。
2.如权利要求1所述的针对低幅度构造畸变的速度校正方法,其特征在于,在所述步骤S101中,在构建初始的空变速度场之前,进行井震时深标定,以从中选取曲线质量和标定效果优的井,用于初步速度场构建的纵向约束;
从地震处理输出的速度场需要经过DIX公式转成层速度及平均速度,以为所述空变速度场的构建提供横向变化趋势。
3.如权利要求1所述的针对低幅度构造畸变的速度校正方法,其特征在于,在所述步骤S102中,
选取纯波地震数据体,在构造地层格架的约束下结合井震标定结果分析速度异常体对应的地震反射特征。
4.如权利要求1所述的针对低幅度构造畸变的速度校正方法,其特征在于,在所述步骤S103中,对井上实测速度曲线进行环境校正以及归一化,以确保除地质因素外,不同井对应大致深度的相似岩性段具有接近的速度范围。
5.如权利要求1所述的针对低幅度构造畸变的速度校正方法,其特征在于,在所述步骤S106中,替换速度场具有与原始速度场相同的时间范围、数据体采样率以及近似的量纲范围。
CN201710237118.XA 2017-04-12 2017-04-12 针对低幅度构造畸变的速度校正方法 Active CN107688200B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710237118.XA CN107688200B (zh) 2017-04-12 2017-04-12 针对低幅度构造畸变的速度校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710237118.XA CN107688200B (zh) 2017-04-12 2017-04-12 针对低幅度构造畸变的速度校正方法

Publications (2)

Publication Number Publication Date
CN107688200A CN107688200A (zh) 2018-02-13
CN107688200B true CN107688200B (zh) 2019-09-03

Family

ID=61152428

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710237118.XA Active CN107688200B (zh) 2017-04-12 2017-04-12 针对低幅度构造畸变的速度校正方法

Country Status (1)

Country Link
CN (1) CN107688200B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111208558B (zh) * 2018-11-22 2022-05-06 中国石油天然气股份有限公司 超深低幅度三维地质构造的建立方法及装置
CN112835095B (zh) * 2019-11-25 2023-08-22 中国石油天然气股份有限公司 地震资料低幅度构造成图方法及装置
CN112859167B (zh) * 2019-11-28 2023-08-22 中国石油天然气股份有限公司 畸变地质体速度场的校正方法及装置
CN113655525A (zh) * 2020-05-12 2021-11-16 中国石油天然气集团有限公司 基于高速砾岩的构造高点偏移的校正方法及装置
CN113009577B (zh) * 2021-03-23 2023-07-21 中海石油深海开发有限公司 一种低幅构造校正方法、装置、计算机设备及存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101750628A (zh) * 2008-12-11 2010-06-23 中国石油天然气股份有限公司 二维叠加速度及均方根速度场闭合差校正方法
CN105353412A (zh) * 2015-12-14 2016-02-24 中国石油大学(华东) 一种井震联合平均速度场的计算方法及系统
CN105717540A (zh) * 2016-03-14 2016-06-29 中国海洋石油总公司 一种微幅度构造精细预测方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9091786B2 (en) * 2010-10-15 2015-07-28 Oz YILMAZ Image based effective medium modeling of the near surface earth formation

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101750628A (zh) * 2008-12-11 2010-06-23 中国石油天然气股份有限公司 二维叠加速度及均方根速度场闭合差校正方法
CN105353412A (zh) * 2015-12-14 2016-02-24 中国石油大学(华东) 一种井震联合平均速度场的计算方法及系统
CN105717540A (zh) * 2016-03-14 2016-06-29 中国海洋石油总公司 一种微幅度构造精细预测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
三维地震速度场建立技术_以饶阳凹陷河间南地区为例;易远元 等;《石油学报》;20150731;第36卷(第7期);第820-826页
井震时深转换技术在低幅度构造评价中的应用;梁卫 等;《中国海上油气》;20140630;第26卷(第3期);第61-64、71页

Also Published As

Publication number Publication date
CN107688200A (zh) 2018-02-13

Similar Documents

Publication Publication Date Title
CN107688200B (zh) 针对低幅度构造畸变的速度校正方法
CN108957549B (zh) 一种辫状河沉积非均质致密砂岩气藏地质建模方法
CN104047597B (zh) 含油气泥页岩地层测井曲线标准化方法
CN101551466B (zh) 一种利用与偏移距有关的地震属性提高油气储层预测精度的方法
CN108680951A (zh) 一种基于地震信息判断煤层气富集沉积控制作用的方法
CN106054248A (zh) 一种基于大面积致密储层地震岩石物理反演方法
Zhang et al. Brittleness evaluation of resource plays by integrating petrophysical and seismic data analysis
CN109388817A (zh) 一种储层裂缝三维建模方法
CN103454685A (zh) 利用测井约束波阻抗反演预测砂体厚度的方法和装置
CN103670388A (zh) 一种泥页岩有机碳含量的评价方法
CN109633750A (zh) 基于测井相波阻抗与地震波形的非线性映像关系反演方法
CN107065013B (zh) 一种地震尺度下的层速度确定方法及装置
CN104570067A (zh) 一种地球物理勘探中相控地震反演方法
CN104200115A (zh) 一种基于地质统计学模拟的全地层速度建模方法
CN110954943A (zh) 被动源地震频率谐振勘探方法
Chen et al. The application study on the multi-scales integrated prediction method to fractured reservoir description
CN103777245A (zh) 基于地震资料的油气成藏条件定量评价方法
Zhou et al. Sweet spot prediction in tight sandstone reservoir based on well-bore rock physical simulation
Liu et al. Combined petrophysics and 3D seismic attributes to predict shale reservoirs favourable areas
CN102914797B (zh) 一种获得地层各向异性系数的方法及装置
Sun et al. Shale rocks brittleness index prediction method using extended elastic impedance inversion
Christensen et al. Seismically driven reservoir characterization using an innovative integrated approach: Syd Arne Field
CN110244383A (zh) 基于近地表数据的地质岩性综合模型创建方法
Fernandez-Concheso Characterizing an unconventional reservoir with conventional seismic data: A case study using seismic inversion for the Vaca Muerta Formation, Neuquen Basin, Argentina
CN109283577A (zh) 一种地震层位标定方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CP01 Change in the name or title of a patent holder

Address after: 100094 Room 401, building 4, courtyard 3, fengxiu Middle Road, Haidian District, Beijing

Patentee after: New JCM Group Co.,Ltd.

Address before: 100094 Room 401, building 4, courtyard 3, fengxiu Middle Road, Haidian District, Beijing

Patentee before: LANDOCEAN ENERGY SERVICES CO.,LTD.

CP01 Change in the name or title of a patent holder