CN117192601A - 一种基于时间域层速度的时深转换速度场建立方法 - Google Patents
一种基于时间域层速度的时深转换速度场建立方法 Download PDFInfo
- Publication number
- CN117192601A CN117192601A CN202311116356.7A CN202311116356A CN117192601A CN 117192601 A CN117192601 A CN 117192601A CN 202311116356 A CN202311116356 A CN 202311116356A CN 117192601 A CN117192601 A CN 117192601A
- Authority
- CN
- China
- Prior art keywords
- time
- layer
- speed
- velocity
- well
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 59
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 32
- 238000012937 correction Methods 0.000 claims abstract description 40
- 238000001914 filtration Methods 0.000 claims abstract description 8
- 230000008859 change Effects 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 9
- 238000010276 construction Methods 0.000 claims description 7
- 230000007704 transition Effects 0.000 claims description 4
- 238000012800 visualization Methods 0.000 claims description 3
- 238000005553 drilling Methods 0.000 claims description 2
- 230000015572 biosynthetic process Effects 0.000 abstract description 3
- 238000011161 development Methods 0.000 abstract description 3
- 239000003208 petroleum Substances 0.000 abstract description 3
- 238000003786 synthesis reaction Methods 0.000 abstract 1
- 238000011160 research Methods 0.000 description 9
- 230000000694 effects Effects 0.000 description 7
- 230000008901 benefit Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 238000013508 migration Methods 0.000 description 5
- 230000005012 migration Effects 0.000 description 5
- 238000009933 burial Methods 0.000 description 4
- 238000004587 chromatography analysis Methods 0.000 description 4
- 238000005056 compaction Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000007480 spreading Effects 0.000 description 4
- 230000006835 compression Effects 0.000 description 3
- 238000007906 compression Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000000265 homogenisation Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000005498 polishing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种基于时间域层速度的时深转换速度场建立方法,属于石油勘探开发技术领域。本发明的方法包括以下步骤:1)通过井的时间域合成记录标定获取时深曲线和层速度曲线;2)对井的层速度曲线进行低通滤波处理;3)计算出井与地震之间的层速度误差曲线,对误差进行构造约束的克里金插值;4)原始地震层速度场减去层速度误差场,得到第1次校正后的时间域层速度场;5)采用井控的浅层地震层位时间值横向趋势约束预测、克里金误差校正的方法,建立浅层水平控制时间层以上的填充速度体;6)将第1次校正后的层速度场的浅层部分用填充速度体进行替换,得到最终速度场。本发明提出的速度场建立方法,能够改善时深转换速度场的可靠性。
Description
技术领域
本发明涉及一种基于时间域层速度的时深转换速度场建立方法,属于石油勘探开发领域。
背景技术
随着地震资料处理技术的发展进步,近年来叠前偏地震资料处理技术得到了越来越广泛的应用,地震层速度在地震资料处理中得到了比较普遍的应用。目前主流的地震资料解释技术,通常在时间域内进行,对于在时间域解释得到的构造数据,需要采用较为准确的时间域速度场进行时深转换,才能得到相对真实的深度域构造数据。因此,如何建立比较可靠的时间域速度场,成为构造研究中的一个关键环节。
以往建立时深转换速度场的技术,通常以地震资料处理过程中得到的地震叠加速度为基础,将其转换到平均速度域,然后使用井的平均速度资料进行校正,再用于时深转换研究。但是对于叠前偏处理相关的时间域地震层速度数据,如果采用以往的先转化为平均速度场再使用井速度进行校正的流程,表现出明显的不适应性。主要问题在于:对于时间域的地震层速度,虽然可以采用公式计算的方法将其转换为地震平均速度,但是由于地震层速度在浅层的误差相对较大(一方面由于地震资料通常在浅层的覆盖次数低、信噪比差,地震资料处理过程中难以获得准确的浅层层速度;另一方面由于地震浅层速度在地震资料处理过程中往往采用固定速度填充的简化处理方法,造成浅层速度与实际地层速度的差异较大),而在从地震层速度向地震平均速度进行转换的计算过程中,需要自浅而深累计平均,那么浅层不准确的地震层速度将会对中深层的地震平均速度计算结果产生持续性的不良影响,造成转换后的地震平均速度在中深层(石油勘探的目标层一般位于中深层)与井的平均速度资料缺乏准确的对应关系,使得利用井速度校正地震速度缺乏可靠的基础条件,进而影响到最终速度场的精确性与可靠性。
通过查询现有的基于层速度的速度建场技术,大多属于地震资料处理领域,研究的侧重点主要在于提高地震资料的成像精度,而在时间域地震数据的时深转换速度场建立方面,目前利用井震层速度数据开展的研究相对较少。
例如公开号为CN105277977A的中国专利申请公布文本公开了一种叠前深度偏移速度场模型建立方法和装置,在利用初至层析的方法反演近地表速度模型的过程中,引入了射线密度,通过射线密度的垂向展布来确定初至层析反演模型的垂向有效性,再通过将有效的初至层析的近地表速度模型与深层反射层析的速度模型进行拼接,得到更加精确的浅层速度模型,从而提高深度偏移的成像质量。
2015年的硕士论文《地震波逆时偏移深度域层速度建模方法研究》,提出了构建深度域层速度初始模型及后续修正的方法,首先从共炮点道集抽取CMP道集,进行叠加速度分析,将叠加速度转变为均方根速度,再根据DIX公式转为时间域层速度;然后对时间域层速度进行时深转换,得到深度域初始层速度模型;在深度域通过构建层析速度反演方程组求取速度更新量,自上而下逐层修正初始速度模型,当逆时偏移角道集的同向轴水平时,得到适合于逆时偏移的深度域层速度模型。
发明内容
本发明的目的是提供一种基于时间域层速度的时深转换速度场建立方法,用于建立可靠性更佳的时深转换速度场。
为实现上述目的,本发明的方案包括:
本发明的一种基于时间域层速度的时深转换速度场建立方法的技术方案,包括如下步骤:
1)根据已钻井的时深曲线,得到井的时间域层速度曲线;
2)提取井旁地震时间域层速度曲线,与所述井的时间域层速度曲线进行对比,并对井的时间域层速度曲线进行处理;
3)使用井旁地震时间域层速度曲线减去处理后的井的时间域层速度曲线,得到层速度误差曲线;根据所述层速度误差曲线得到时间域层速度误差场;
4)使用原始地震时间域层速度场减去时间域层速度误差场,得到第1次校正后的时间域层速度场;
5)在浅层选取一个井控程度高的水平控制时间层,采用浅层地震层位时间值横向趋势约束预测、克里金误差校正的方法,得到浅层横向变化、垂向均质的填充速度体;
6)通过所述填充速度体对第1次校正后的时间域层速度场进行调整,得到第2次校正后的最终时间域层速度场;基于DIX公式将所述最终时间域层速度场转换为平均速度场,完成建立过程。
本发明采用浅层地震层位时间值横向趋势约束预测、克里金误差校正的方法,能够使稀疏井点的浅层填充速度的平面插值结果与实际地震时间层位的平面趋势较为相似。这是有其合理性的。从理论上分析,位于水平控制时间面到零时间面之间的浅层地层的填充速度的横向分布趋势,与地层的新老、岩性的分布、埋深压实作用等因素存在着相关性;而地层的新老、岩性的分布、埋深压实作用等因素在水平方向的展布趋势,与附近非水平分布的实际地震层位的时间值的平面趋势具有一定程度的相关性;因此浅层地层的填充速度的横向分布趋势,应该与附近非水平分布的实际地震层位的时间值的横向分布趋势具有相关性。基于这种认识,可以采用浅层地震层位时间值横向趋势约束预测、克里金插值误差校正的方法,改善稀疏井点浅层填充速度在平面上插值结果的可靠性。
本发明的基于时间域层速度的时深转换速度场建立方法,能够利用地震速度的趋势约束优势和井速度的可靠性好及细节丰富优势,建立综合性的速度场;在时间域层速度场的校正方面,直接在层速度域内进行而不是在平均速度域内进行,避免了地震浅层层速度的不准确性带来的在平均速度域地震速度与井速度存在较大的趋势误差问题,有助于改善利用井速度校正地震速度的匹配性;对于浅层水平控制时间层以上的速度场填充速度,采用浅层地震层位时间值横向趋势约束预测、克里金误差校正的方法进行获取,改善浅层填充速度在稀疏井插值条件下的可靠性,从而得到较为可靠的速度场,应用于时深转换研究。
进一步地,所述步骤5)中,通过如下步骤计算得到所述填充速度体:
在浅层选取一个水平控制时间层位,所述层位满足使设定数量的已钻井在水平控制时间面位置具有时深曲线值;
使用在水平控制时间层位具有井的时深曲线值的已钻井,根据水平控制时间层位的时间值在井的时深曲线上查找,得到该水平控制时间值在井上所对应的深度值;
根据已钻井相应于水平控制时间层位的时间值、深度值,得到已钻井位置的浅层填充速度;
在水平控制时间层位附近,使用地震波数据体追踪产生一个实际地震时间层位,统计已钻井点坐标位置在该实际地震时间层位的时间值;以该时间值为横轴,以所述浅层填充速度为纵轴,绘制出层位时间值与浅层填充速度值的交会图,根据所述交会图拟合得到从层位时间值到浅层填充速度值的转换关系式,使用该关系式对所述实际地震时间层位进行转换,得到初步浅层填充速度平面图;统计初步浅层填充速度平面图与已钻井的浅层填充速度的误差,使用克里金插值方法对误差进行平面插值分配;使用初步浅层填充速度平面图减去误差平面分配结果,得到最终的浅层平面填充速度;
将所述浅层平面填充速度三维立体化,三维体的垂向范围采用从水平控制时间面至零时间面的范围,平面范围采用地震速度场的范围,得到所述填充速度体。
进一步地,所述步骤1)中,根据所述时深曲线上的数据点,采用“相邻点深度差÷相邻点时间差=相邻点之间的层速度”的计算方法,得到井的时间域层速度曲线。
进一步地,在时间域对已钻井进行合成记录标定,根据波形匹配情况进行拉伸压缩,得到已钻井的时深曲线。
进一步地,所述步骤2)中,对井的时间域层速度曲线进行低通滤波处理,使处理后的井的时间域层速度曲线的变化细节略多于地震时间域层速度曲线的变化细节。
进一步地,所述步骤3)中,使用时间域构造格架做为约束,对所述层速度误差曲线进行克里金插值,得到时间域层速度误差场。
进一步地,所述步骤6)中,所述调整包括:根据第1次校正后的时间域层速度场,将其浅层水平控制时间层以上至零时间面之间的部分用所述填充速度体进行替换。
本发明首先在时间域对已钻井进行合成记录标定,根据波形匹配情况进行适当的拉伸压缩,得到已钻井的时深曲线;根据该时深曲线上的数据点,采用“相邻点深度差÷相邻点时间差=相邻点之间的层速度”的计算方法,得到井的时间域层速度曲线;
然后提取井旁地震时间域层速度曲线,与上述井的层速度曲线进行对比;对井的层速度曲线进行低通滤波处理,使处理后的井的层速度曲线的变化细节略多于地震层速度的变化细节;
然后使用井旁地震时间域层速度曲线减去滤波处理后的井的时间域层速度曲线,得到层速度误差曲线;使用时间域构造格架做为约束,对层速度误差曲线进行克里金插值,得到时间域层速度误差场;
然后使用原始地震时间域层速度场减去时间域层速度误差场,得到第1次校正后的时间域层速度场;
然后在浅层选取一个井控程度较高的水平控制时间层,采用浅层地震层位时间值横向趋势约束预测、克里金误差校正的方法,得到浅层横向变化、垂向均质的填充速度体;
最后根据第1次校正后的时间域层速度场,将其浅层水平控制时间层以上至零时间面之间的部分用上述浅层横向变化、垂向均质的填充速度体进行替换,得到第2次校正后的时间域层速度场,也就是最终的层速度场;基于DIX公式将该层速度场转换为平均速度场,用于时深转换研究。
其中,计算得到浅层横向变化的填充速度体的方法包括以下步骤:
在浅层选取一个水平控制时间层位,该层位需要满足使大约80%以上的已钻井在水平控制时间面位置具有时深曲线值(在水平控制时间层满足高于所有目标层时间段的前提下,此百分比数值越高越好,亦即井控程度越高越好);
使用在水平控制时间面具有井的时深曲线值的已钻井,根据水平控制时间层的时间值在井的时深曲线上查找,得到该水平控制时间值在井上所对应的深度值;
根据已钻井相应于水平控制时间面的时间值、深度值,采用“深度÷时间=速度”的计算方法,得到已钻井位置的浅层填充速度;
在水平控制时间层附近,使用地震波数据体追踪产生一个实际地震时间层位(实际地震的时间层位必然具有非水平分布特征),统计已钻井点坐标位置在该实际地震时间层位的时间值;以该时间值为横轴,以前面得到的井点的浅层填充速度为纵轴,绘制出层位时间值与浅层填充速度值的交会图,在该图上拟合得到从层位时间值到浅层填充速度值的转换关系式,使用该关系式对上述浅层实际地震时间层位进行转换,得到初步的浅层填充速度平面图;统计初步浅层填充速度平面图与已钻井的浅层填充速度的误差,使用克里金插值方法对误差进行平面插值分配;使用初步的浅层填充速度平面图减去误差平面分配结果,得到最终的浅层平面填充速度;
将上述的浅层平面填充速度三维立体化,三维体的垂向范围采用从水平控制时间面至零时间面的范围,平面范围采用地震速度场的范围,三维体中每个样点的属性值为浅层填充速度值,该值在浅层三维体内具有横向变化、垂向均质的特点。
附图说明
图1是本发明的基于时间域层速度的时深转换速度场建立方法流程图;
图2(a)是本发明实施例中井的时深曲线图;
图2(b)是本发明实施例中井的层速度曲线图;
图3(a)是本发明实施例中地震层速度曲线图;
图3(b)是本发明实施例中井的滤波后层速度曲线图;
图4是本发明实施例中井震层速度误差曲线与误差体剖面叠合图;
图5(a)是本发明实施例中井旁地震平均速度曲线图;
图5(b)是本发明实施例中井的平均速度曲线图;
图6是本发明实施例中第1次校正后的层速度场图;
图7是本发明实施例中的浅层实际地震时间层位在平面的分布图;
图8是本发明实施例中的浅层实际地震层位时间值与浅层填充速度的交会图;
图9是本发明实施例中的初步浅层填充速度平面图;
图10是本发明实施例中的浅层填充速度误差的克里金插值平面图;
图11是本发明实施例中的最终浅层填充速度平面图;
图12是本发明实施例中的浅层填充速度在剖面的分布图;
图13是本发明实施例中仅采用井的浅层填充速度进行克里金插值(未用浅层地震层位时间值横向趋势约束)得到的填充速度平面图;
图14是本发明实施例中采用浅层填充速度替换后的层速度体剖面图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明了,下面结合附图及实施例,对本发明作进一步的详细说明。
研究区面积约150km2,是有利的油气勘探区,构造趋势的高低对油气成藏综合研究具有重要作用。该区已钻井11口,所存在的技术问题是,如何合理采用井震结合的研究方法,建立可靠性较好的时深转换速度场,为构造研究奠定良好的速度基础。
本实施例的一种基于时间域层速度的时深转换速度场建立方法,流程如图1所示,包括以下步骤:
1)在时间域对已钻井进行合成记录标定,根据波形匹配情况进行适当的拉伸压缩,得到已钻井的时深曲线;根据该时深曲线上的数据点,采用“相邻点深度差÷相邻点时间差=相邻点之间的层速度”的计算方法,得到井的时间域层速度曲线。
以实施例中的W9井为例,井的时深曲线、层速度曲线分别如图2(a)、图2(b)所示。其中井的层速度曲线不是直接根据声波曲线计算得到,而是根据井的合成记录时深曲线计算得到。
2)提取井旁地震时间域层速度曲线,与上述井的层速度曲线进行对比;对井的层速度曲线进行低通滤波处理,使处理后的井的层速度曲线的变化细节略多于地震层速度的变化细节。
以实施例中的W9井为例,井旁地震层速度曲线、井的低通滤波后层速度曲线分别如图3(a)、图3(b)所示。
3)使用井旁地震时间域层速度曲线减去滤波处理后的井的时间域层速度曲线,得到层速度误差曲线;使用时间域构造格架做为约束,对层速度误差曲线进行克里金插值,得到层速度误差场。
在本实施例中,以过W1至W5井的层速度误差曲线与层速度误差场的剖面叠合图为例,对层速度误差场的插值效果进行说明,如图4所示,图中A层位、B层位是根据地震波数据体追踪得到的时间域地震层位,层位的空间高低趋势控制了井震层速度误差在三维空间的插值趋势。
这里说明一下为什么要在层速度域对井震速度进行误差曲线提取,而不是采用常规的在平均速度域进行误差曲线提取。图5(a)、图5(b)展示了使用图3(a)、图3(b)的层速度曲线依据DIX公式转换得到的地震平均速度曲线、井的平均速度曲线图,从图5(a)、图5(b)可以看出,在平均速度域,从浅层至中深层的井震平均速度存在较大的误差,产生这种现象的原因主要是由于图3(a)所示的浅层地震层速度并不准确(可以看出图3(a)的最浅层地震层速度使用了固定速度2000m/s进行填充,然后向下逐渐过渡),造成转换得到的地震平均速度曲线(图5(a))从浅层至中深层都不够准确可靠,因此缺乏在平均速度域内进行井震速度对比和校正的条件;而对比观察图3(a)、图3(b)所示的地震层速度曲线、井的层速度曲线特征,两者在中深层的层速度趋势符合较好,且层速度误差明显较小,说明在层速度域内井震速度之间具有较好的匹配关系,具备良好的井震速度对比和校正条件。
4)使用原始地震时间域层速度场减去时间域层速度误差场,得到第1次校正后的层速度场。
在本实施例中,图6为过W1至W5井的第1次校正后的层速度场剖面,可以看出层速度分布具有比较和谐的渐变特征,没有发生畸变。
5)在浅层选取一个井控程度较高的水平控制时间层,采用浅层地震层位时间值横向趋势约束预测、克里金误差校正的方法,计算得到浅层横向变化的填充速度体。
在步骤5)中,计算得到浅层横向变化的填充速度体的方法包括以下步骤:
在浅层选取一个水平控制时间层位,该层位需要满足使大约80%以上的已钻井在水平控制时间面位置具有时深曲线值(在水平控制时间层满足高于所有目标层时间段的前提下,此百分比数值越高越好,亦即井控程度越高越好)。根据上述原则,本实施例中优选出的水平控制时间层位为时间值等于1500ms的水平时间层位,该层位能够满足使100%的已钻井在时间值为1500ms的水平时间面上具有时深曲线值,且该1500ms水平时间层高于该研究区有利目的层段所在的1650-3000ms的时间分布范围。
使用在水平控制时间面具有井的时深曲线值的已钻井,根据水平控制时间层的时间值在井的时深曲线上查找,得到该水平控制时间值在井上所对应的深度值。
根据已钻井相应于水平控制时间面的时间值、深度值,采用“深度÷时间=速度”的计算方法,得到已钻井位置的浅层填充速度。
在水平控制时间层附近,使用地震波数据体追踪产生一个实际地震时间层位(实际地震的时间层位必然具有非水平分布特征),采用地震层位时间值横向趋势约束预测、克里金误差校正的方法,得到浅层平面填充速度。
在本实施例中,通过在图6的剖面图中进行观察,发现距离1500ms的水平控制时间层较近的实际地震层位是A层位,因此选择A层位做为浅层填充速度的横向趋势约束层位,A层位在平面上的分布如图7所示。统计已钻井点位置在A层位的时间值,以该时间值为横轴,以前面得到的井点的浅层填充速度为纵轴,绘制出层位时间值与浅层填充速度值的交会图,如图8所示,在该图上拟合得到从层位时间值到浅层填充速度值的转换关系式为“填充速度=1.8069×层位时间-919.1865”;使用该关系式对A层位进行转换,得到初步的浅层填充速度平面图,如图9所示;统计初步浅层填充速度平面图与已钻井的浅层填充速度的误差,使用克里金插值方法对误差进行平面插值分配,如图10所示;使用初步的浅层填充速度平面图减去误差平面分配结果,得到最终的浅层平面填充速度,如图11所示。对比最终的浅层平面填充速度图(图11)与A层位的时间值平面分布图(图7)可以看出,采用上述地震层位时间值横向趋势约束预测、克里金误差校正的方法,能够使稀疏井点的浅层填充速度的平面插值结果的趋势,与浅层实际地震时间层位的平面趋势较为相似,从而使浅层填充速度的平面插值结果较为可靠。
从理论上分析,位于水平控制时间面到零时间面之间的浅层地层的填充速度的横向分布趋势,与地层的新老、岩性的分布、埋深压实作用等因素存在着相关性;而地层的新老、岩性的分布、埋深压实作用等因素在水平方向的展布趋势,与附近非水平分布的实际地震层位的时间值的平面趋势具有一定程度的相关性;因此浅层地层的填充速度的横向分布趋势,应该与附近非水平分布的实际地震层位的时间值的横向分布趋势具有相关性。基于这种认识,可以采用浅层地震层位时间值横向趋势约束预测、克里金插值误差校正的方法,改善稀疏井点浅层填充速度在平面上插值结果的可靠性。
将上述平面填充速度三维立体化,三维体的垂向范围采用从水平控制时间面至零时间面的范围,平面范围采用地震速度场的范围,三维体中每个样点的属性值为浅层填充速度值,该值在浅层三维体内具有横向变化、垂向均质的特点。选取过W1至W5井的填充速度剖面,对其空间分布特征进行说明,如图12所示。
在浅层填充速度的获取方面,如果不采用本发明提出的地震层位时间值横向趋势约束预测、克里金误差校正的方法,而仅采用井点的浅层速度直接进行克里金插值,得到的浅层填充速度平面图如图13所示,可以看出,该图的填充速度插值结果趋势受已钻井点的平面分布影响较大,容易产生局部的插值畸变,且这种填充速度插值结果没有考虑速度与宏观地层展布趋势的关联性,得到的浅层填充速度分布结果不够合理可靠。
6)根据第1次校正后的时间域层速度场,将其浅层水平控制时间层以上至零时间面之间的部分用上述浅层横向变化、垂向均质的填充速度体进行替换,得到第2次校正后的时间域层速度场,也就是最终的层速度场,层速度场的剖面效果如图14所示;基于DIX公式将该层速度场转换为平均速度场,用于时深转换研究。
从上述实施例可以看出,本发明的基于时间域层速度的时深转换速度场建立方法,能够利用地震速度的趋势约束优势和井速度的可靠性好及细节丰富优势,建立综合性的速度场;在地震速度场的校正方面,直接在层速度域内进行而不是在平均速度域内进行,避免了地震浅层层速度的不准确性带来的在平均速度域地震速度与井速度存在较大的趋势误差问题,从而改善利用井速度校正地震速度的匹配性;对于浅层水平控制时间层以上的速度场填充速度,采用附近地震层位时间值横向趋势约束预测、克里金误差校正的方法获取,有助于改善浅层填充速度在稀疏井插值条件下的可靠性,最终得到较为可靠的时深转换速度场。
针对如何更有效地利用时间域地震层速度数据、井的层速度数据进行联合分析,建立可靠性较好的时深转换速度场的问题,本发明专利提出了一种解决方案。首先利用时间域井的合成记录层速度资料对时间域地震层速度进行构造约束的克里金误差校正,得到初步的层速度场;然后利用已钻井在浅层水平控制时间面的层速度资料,使用浅层地震层位时间值横向趋势约束预测、克里金误差校正的方法,得到速度场的浅层填充速度;对时间域层速度场的水平控制时间层以上的部分使用横向变化、垂向均质的速度插值充填,得到最终的时间域层速度场,将其转换为时间域平均速度场,用于时深转换研究。
Claims (7)
1.一种基于时间域层速度的时深转换速度场建立方法,其特征在于,包括如下步骤:
1)根据已钻井的时深曲线,得到井的时间域层速度曲线;
2)提取井旁地震时间域层速度曲线,与所述井的时间域层速度曲线进行对比,并对井的时间域层速度曲线进行处理;
3)使用井旁地震时间域层速度曲线减去处理后的井的时间域层速度曲线,得到层速度误差曲线;根据所述层速度误差曲线得到时间域层速度误差场;
4)使用原始地震时间域层速度场减去时间域层速度误差场,得到第1次校正后的时间域层速度场;
5)在浅层选取一个井控程度高的水平控制时间层,采用浅层地震层位时间值横向趋势约束预测、克里金误差校正的方法,得到浅层横向变化、垂向均质的填充速度体;
6)通过所述填充速度体对第1次校正后的时间域层速度场进行调整,得到第2次校正后的最终时间域层速度场;基于DIX公式将所述最终时间域层速度场转换为平均速度场,完成建立过程。
2.根据权利要求1所述的基于时间域层速度的时深转换速度场建立方法,其特征在于,所述步骤5)中,通过如下步骤计算得到所述填充速度体:
在浅层选取一个水平控制时间层位,所述层位满足使设定数量的已钻井在水平控制时间面位置具有时深曲线值;
使用在水平控制时间层位具有井的时深曲线值的已钻井,根据水平控制时间层位的时间值在井的时深曲线上查找,得到该水平控制时间值在井上所对应的深度值;
根据已钻井相应于水平控制时间层位的时间值、深度值,得到已钻井位置的浅层填充速度;
在水平控制时间层位附近,使用地震波数据体追踪产生一个实际地震时间层位,统计已钻井点坐标位置在该实际地震时间层位的时间值;以该时间值为横轴,以所述浅层填充速度为纵轴,绘制出层位时间值与浅层填充速度值的交会图,根据所述交会图拟合得到从层位时间值到浅层填充速度值的转换关系式,使用该关系式对所述实际地震时间层位进行转换,得到初步浅层填充速度平面图;统计初步浅层填充速度平面图与已钻井的浅层填充速度的误差,使用克里金插值方法对误差进行平面插值分配;使用初步浅层填充速度平面图减去误差平面分配结果,得到最终的浅层平面填充速度;
将所述浅层平面填充速度三维立体化,三维体的垂向范围采用从水平控制时间面至零时间面的范围,平面范围采用地震速度场的范围,得到所述填充速度体。
3.根据权利要求1所述的基于时间域层速度的时深转换速度场建立方法,其特征在于,所述步骤1)中,根据所述时深曲线上的数据点,采用“相邻点深度差÷相邻点时间差=相邻点之间的层速度”的计算方法,得到井的时间域层速度曲线。
4.根据权利要求3所述的基于时间域层速度的时深转换速度场建立方法,其特征在于,在时间域对已钻井进行合成记录标定,根据波形匹配情况进行拉伸压缩,得到已钻井的时深曲线。
5.根据权利要求1所述的基于时间域层速度的时深转换速度场建立方法,其特征在于,所述步骤2)中,对井的时间域层速度曲线进行低通滤波处理,使处理后的井的时间域层速度曲线的变化细节略多于地震时间域层速度曲线的变化细节。
6.根据权利要求1所述的基于时间域层速度的时深转换速度场建立方法,其特征在于,所述步骤3)中,使用时间域构造格架做为约束,对所述层速度误差曲线进行克里金插值,得到时间域层速度误差场。
7.根据权利要求1所述的基于时间域层速度的时深转换速度场建立方法,其特征在于,所述步骤6)中,所述调整包括:根据第1次校正后的时间域层速度场,将其浅层水平控制时间层以上至零时间面之间的部分用所述填充速度体进行替换。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311116356.7A CN117192601A (zh) | 2023-08-31 | 2023-08-31 | 一种基于时间域层速度的时深转换速度场建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311116356.7A CN117192601A (zh) | 2023-08-31 | 2023-08-31 | 一种基于时间域层速度的时深转换速度场建立方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117192601A true CN117192601A (zh) | 2023-12-08 |
Family
ID=88999068
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311116356.7A Pending CN117192601A (zh) | 2023-08-31 | 2023-08-31 | 一种基于时间域层速度的时深转换速度场建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117192601A (zh) |
-
2023
- 2023-08-31 CN CN202311116356.7A patent/CN117192601A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7400978B2 (en) | Method for updating a geologic model by seismic data | |
CN109738945B (zh) | 一种利用叠前深度偏移成果直接生成构造图的方法 | |
CN108196305B (zh) | 一种山地静校正方法 | |
CN110361781B (zh) | 一种基于地震反射振幅分析的沉积期古地貌恢复方法 | |
CN111722284B (zh) | 一种基于道集数据建立速度深度模型的方法 | |
CN105089652A (zh) | 一种拟声波曲线重构与稀疏脉冲联合反演方法 | |
CN107942379A (zh) | 一种提高复杂断块速度模型精度的方法 | |
CN113031068B (zh) | 一种基于反射系数精确式的基追踪叠前地震反演方法 | |
CN109444959A (zh) | 全频高精度层速度场建立方法 | |
CN113093286B (zh) | 一种冷泉发育区储层非均质性的反演方法 | |
CN111077578B (zh) | 岩层分布预测方法和装置 | |
CN109143335B (zh) | 一种合成地震记录的制作方法、系统、介质及设备 | |
CN117192601A (zh) | 一种基于时间域层速度的时深转换速度场建立方法 | |
CN111060967B (zh) | 一种井控速度场建模方法 | |
CN116840909A (zh) | 一种针对潜山内幕的联合深度域建模方法 | |
CN113589375B (zh) | 基于倾斜层约束旅行时计算的vsp层速度反演方法 | |
CN112946742A (zh) | 一种拾取精确叠加速度谱的方法 | |
CN111310361A (zh) | 基于随钻地震的钻井指导方法、系统、设备及存储介质 | |
CN110244383A (zh) | 基于近地表数据的地质岩性综合模型创建方法 | |
CN109061737A (zh) | 一种基于合成地震记录的储层预测方法及装置 | |
CN114814949B (zh) | 一种浅层逆vsp初至层析及地层预测方法 | |
CN114966826B (zh) | 基于野外近地表调查资料的高低频融合近地表建模方法 | |
CN113589366B (zh) | 基于全波形反演的宽频融合建模方法 | |
CN114076982B (zh) | 一种基于波形特征差异的岩溶古地貌恢复方法及装置 | |
CN112526610B (zh) | 一种约束表层建模的三维地震采集激发井深设计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication |