CN109709605B - 火成岩复杂区碳酸盐岩断溶体地震数据成像方法及设备 - Google Patents
火成岩复杂区碳酸盐岩断溶体地震数据成像方法及设备 Download PDFInfo
- Publication number
- CN109709605B CN109709605B CN201910149938.2A CN201910149938A CN109709605B CN 109709605 B CN109709605 B CN 109709605B CN 201910149938 A CN201910149938 A CN 201910149938A CN 109709605 B CN109709605 B CN 109709605B
- Authority
- CN
- China
- Prior art keywords
- seismic
- velocity
- establishing
- velocity model
- imaging
- 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
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明属于地震信号处理领域,特别涉及一种火成岩复杂区碳酸盐岩断溶体地震数据成像方法及设备。所述地震数据成像方法包括:对地震数据进行预处理,获得保幅的地震道集数据;利用地质约束下火成岩的精细雕刻技术建立二叠系火成岩地层速度模型;利用分频处理方法得到优势频带的地震数据,建立奥陶系地层和寒武系地层的背景速度模型;在地震属性约束下进行断控高斯束层析,更新深层断裂处的速度模型获得高精度速度模型;通过宽频RTM成像技术获得高精度的成像结果。本发明提高了沙漠地表火成岩下覆深层断溶体地震数据处理的成像精度,为顺北油气田和国内外相似背景地区断溶体的地震成像提供技术支撑。
Description
技术领域
本发明属于地震信号处理领域,特别涉及一种火成岩复杂区碳酸盐岩断溶体地震数据成像方法及设备。
背景技术
在塔里木盆地沙漠腹地的火成岩复杂区,受二叠系火山活动影响,区域内火成岩非常发育,二叠系火成岩地层存在岩性变化大(英安岩、玄武岩、火山碎屑岩等岩性)、横向展布不连续、纵向厚度不均,由于在不同岩性的火成岩中地震波传播速度变化大,造成时间域的地震数据成像中火成岩下覆出现假背斜、假向斜或者假的断裂,如果速度建模不准,即使叠前深度偏移也不能解决火成岩对下覆地层成像的影响,导致油气勘探开发井位部署失败、钻井目标不落实等问题。
断溶体是深层碳酸盐岩中以断裂、洞穴和溶蚀缝洞为油气储集空间的地震成像目标,内部结构复杂,储集空间纵横向非均质性强,埋深通常超过7000m。发育于奥陶系和寒武系的深层断裂是规模断溶体的主要控制因素。地表沙漠对地震信号产生吸收衰减作用,奥陶系碳酸盐岩地层地震数据有效信号弱,寒武系分辨率低不利于隐蔽断裂、次级断裂的识别。
因此,二叠系火成岩地层、奥陶系地层、寒武系地层和深层断裂带的高精度速度建模和高质量的地震成像是沙漠地表火成岩复杂区深层碳酸盐岩断溶体有效识别的基础。
在常规的地震成像中,所采用的技术方法和工作流程并未对影响断溶体成像质量的“三层一带”采取针对性的处理措施和流程优化,主要技术为:
1)初始速度模型建立技术
利用叠前时间偏移最终的时间均方根速度用DIX公式转换成深度层速度体求取初始深度层速度体。均方根速度转换到深度层速度后会产生许多畸变,使得深度层速度场变化很大,必须进行编辑及平滑,得到初始深度层速度体。
2)模型层析修正背景速度场
基于模型层析成像主要考虑构造模型大套层位的平均层速度,对层间的层速度具有平均效应。沿着构造层位的规则网格进行全局寻优走时算法迭代修正深度域速度模型,垂向上是不规则的,依赖于构造模型。
3)网格层析反演速度建模技术
网格层析成像是基于模型的层析成像技术的拓展,利用CIP道集内同相轴是否拉平作为判别标准。网格层析反演基于数据驱动算法,根据目标线叠前深度偏移得到的共成像点道集,拾取深度剩余速度,形成深度剩余速度体,结合需要更新的深度速度体、拾取的内部层位等,利用最小二乘法,求解网格层析成像矩阵,得到优化后的深度域层速度体。该方法适避免了传统沿层位拾取方法的人为干扰因素,提高了深度域速度建模的精度。
4)Kirchhoff积分算法叠前深度偏移技术
Kirchhoff积分算法叠前深度偏移是当前应用和解决复杂构造和速度横向变化剧烈区域的地震成像的常规技术。深度层速度模型建立和叠前深度偏移是其关键的两大步骤。在获得较为准确的深度层速度模型以及经过合适的偏移参数测试后,对整个数据体进行叠前深度偏移。
发明内容
针对现有技术的缺陷,提供了一种火成岩复杂区深层碳酸盐岩断溶体地震数据成像方法及设备,提高了沙漠地表火成岩下覆深层断溶体地震数据处理的成像精度,为顺北油气田和国内外相似背景地区断溶体的地震成像提供技术支撑。
具体的,本发明提供了一种火成岩复杂区碳酸盐岩断溶体地震数据成像方法,包括:
对地震数据进行预处理,获得保幅的地震道集数据;
利用地质约束下火成岩的精细雕刻技术建立二叠系火成岩地层速度模型;
利用分频处理方法得到优势频带的地震数据,建立奥陶系地层和寒武系地层的背景速度模型;
在地震属性约束下进行断控高斯束层析,更新深层断裂处的速度模型获得高精度速度模型;
通过宽频RTM成像技术获得高精度的成像结果。
优选的,所述预处理包括:低频保护处理和弱信号恢复处理。
优选的,所述建立二叠系火成岩地层速度模型包括:
通过层位拾取建立层位模型;
根据时间偏移速度和测井速度建立初始速度模型;
以1000米间隔进行目标线偏移,根据模型层析建立较为准确的速度模型;
针对二叠系火成岩地层,以细网格的间隔偏移数据体进行网格层析速度模型修正。
进一步优选的,所述细网格的间隔为100m×100m及其以下。
优选的,所述建立所述背景速度模型之后,还包括:利用分频处理方法得到中低频数据,建立奥陶系地层和寒武系地层的背景速度模型,在此基础上利用分频处理方法得到高频数据,基于高频信息对背景速度模型迭代更新获得精细的寒武系地层速度模型。
优选的,所述宽频RTM成像技术选择预设的波场进行成像。
另一方面,本发明提供了一种火成岩复杂区碳酸盐岩断溶体地震信号成像服务器,包括:
地震道集数据获取模块,用于对所述地震数据进行预处理,获得保幅的地震道集数据;
火成岩地层速度模型建立模块,用于利用地质约束下所述火成岩的精细雕刻技术建立二叠系火成岩地层速度模型;
背景速度模型建立模块,用于利用分频处理方法得到优势频带的地震数据,建立奥陶系地层和寒武系地层的背景速度模型;
高精度速度模型生成模块,用于在地震属性约束下进行断控高斯束层析更新深层断裂处的速度模型获得高精度速度模型;
图像生成模块,用于通过宽频RTM成像技术获得高精度的成像结果。
优选的,所述服务器还包括:
地震数据获取模块,用于利用地震数据处理方法获得所述地震数据。
优选的,所述建立二叠系火成岩地层速度模型包括:
通过层位拾取建立层位模型;
根据时间偏移速度和测井速度建立初始速度模型;
以1000米间隔进行目标线偏移,根据模型层析建立较为准确的速度模型;
针对二叠系火成岩地层,以细网格的间隔偏移数据体进行网格层析速度模型修正。
进一步优选的,所述细网格的间隔为100m×100m及以下。
又一方面,本发明提供了一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现上述火成岩复杂区碳酸盐岩断溶体地震数据成像方法。
又一方面,本发明提供了一种计算机设备,包括存储器,处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现上述火成岩复杂区碳酸盐岩断溶体地震数据成像方法。
本发明技术方案的有益效果为:
(1)本发明面向火成岩下覆深层碳酸盐岩断溶体这一成像目标,聚焦关键地层“三层一带”(二叠系火成岩地层、奥陶系碳酸盐岩地层、寒武系地层和深层断裂带)的速度建模,在深入研究其地震信号特征的基础上,通过精细雕刻火成岩速度建模技术、密网格层析成像速度更新方法以及地震属性约束下的断控高斯束层析速度建模方法,结合信号分频技术分别利用低频、高频优势频带逐步获得稳定背景速度模型和相对高频速度模型,可以稳健而细致地建立高精度的速度模型,在此基础上通过宽频RTM成像技术获得高精度的地震数据成像,为断溶体油气藏的勘探开发提供可靠的地震数据处理成果;
(2)本发明的火成岩复杂区深层碳酸盐岩断溶体的地震数据精细地震成像的方法,适用于塔里木盆地火成岩发育区以及国内外其他类似地质背景的地震数据精细地震成像;方法可操作性强,丰富了塔里木盆地地震数据处理的理论与技术方法,为后续顺北油气田断溶体的有效识别、综合地质分析和目标靶区的优选提供了技术支撑,对提高地震勘探技求水平、促进油田增储上产具有重要的意义和良好的应用前景。
附图说明
通过阅读下文优选实施方式的详细描述,各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的,而并不认为是对本发明的限制。而且在整个附图中,用相同的参考符号表示相同的部件。在附图中:
图1为本发明一个实施例提供的沙漠地表火成岩复杂区碳酸盐岩断溶体地震数据成像方法的流程图;
图2为本发明另一个实施例提供的沙漠地表火成岩复杂区深层碳酸盐岩断溶体地震数据成像方法的流程图;
图3为本发明的火成岩复杂区深层碳酸盐岩断溶体地震数据成像服务器示意图;
图4A为常规速度建模技术建立的速度模型的剖面图;
图4B为本发明的地质约束下火成岩精细雕刻技术建立的火成岩高精度速度模型的剖面图;
图4C为常规速度建模技术建立的速度模型的切片图;
图4D为本发明的地质约束下火成岩精细雕刻技术建立的火成岩高精度速度模型的切片图;
图5为基于常规速度建模技术的速度模型的火成岩层地震成像剖面;
图6为本发明的基于地质约束下火成岩精细雕刻技术建立的速度模型的地震成像剖面;
图7为常规技术的断溶体成像剖面;
图8为本发明的利用“三层一带”速度建模技术和宽频RTM成像的断溶体成像剖面。
具体实施方式
下面将参照附图更详细地描述本发明的示例性实施例。虽然附图中显示了本发明的示例性实施例,然而应当理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更透彻地理解本发明,并且能够将本发明的范围完整的传达给本领域的技术人员。
聚焦“三层一带”速度建模是提高火成岩下覆深层碳酸盐岩断溶体成像精度的重要思路。速度建模的精度和成像算法的可靠性是影响地震成像质量的关键。针对重点层位和目的层,在深入研究其地震信号特征的基础上,通过精细雕刻速度建模技术、密网格层析成像速度修正更新方法以及地震属性约束下的断控高斯束层析更新速度模型方法,结合信号分频技术分别利用低频、高频优势频带逐步获得稳定背景速度模型和相对高频速度模型,可以稳健而细致地建立高精度的速度模型。基于高精度的速度模型,RTM技术可以对火成岩地层、奥陶系碳酸盐岩地层、寒武系基底地层和深层断裂带进行高保真的地震成像。
具体的,图1示出了本发明一个实施例提供的火成岩复杂区碳酸盐岩断溶体地震数据成像方法的流程图,如图1所示,本发明的火成岩复杂区碳酸盐岩断溶体地震数据成像方法包括:
步骤S110:对地震数据进行预处理,获得保幅的地震道集数据。
其中,所述地震数据通过利用地震勘探方法获得,所述预处理包括低频保护处理和弱信号恢复处理。
对于沙漠区深层碳酸盐岩断溶体,由于地表沙漠和浅中层地层对地震信号的吸收衰减作用,地震信号传播到深层后有效信号能量弱,尤其奥陶系碳酸盐岩地层非均质性强,波组特征和异常振幅响应更不明显,需要针对奥陶系碳酸盐岩目的层进行弱信号提取和恢复。相对于高频的吸收衰减作用更强,低频有效信号对深层地质目标体的识别更有效,所以需要在预处理过程的噪声衰减、反褶积、滤波等关键步骤中注重低频信号的保护,以获得相对保幅的地震道集数据。
步骤S120:利用地质约束下火成岩的精细雕刻技术建立二叠系火成岩地层速度模型。
基于步骤S110产生的高质量的地震道集数据,利用速度建模及其迭代更新方法,建立二叠系火成岩地层速度模型。包括:通过层位拾取建立层位模型;根据时间偏移速度和测井速度建立初始速度模型;以1000米间隔进行目标线偏移,根据模型层析建立较为准确的速度模型;针对二叠系火成岩地层,以细网格的间隔偏移数据体进行网格层析速度模型修正。
由于火成岩的横向边界变化属于一种高频变化,完全数据驱动的网格层析不能完全刻画出岩相变化的边界,需要加入地质的约束。在成像数据体上拾取二叠系火成岩的顶底地层界面,结合地震属性和地质认识,基于不同火成岩岩性的地震相雕刻火成岩岩相变化的边界,在测井数据的约束下,对火成岩地层充填不同的速度信息直至消除火成岩底界明显的扭曲,并且使火成岩底界的CIP道集拉平。
步骤S130:利用分频处理方法得到优势频带的地震数据,建立奥陶系地层和寒武系地层的背景速度模型。
由于低频数据有很好的穿透性,低频信息在深层属于优势频带,因此,优势频带的地震数据为低频数据。由S110的地震道集数据和S120的速度模型数据经过叠前深度偏移,产生CIP道集数据,针对奥陶系的网格层析成像,需要正确拾取CIP道集的剩余延迟信息,而通常奥陶系地层的地震数据信噪比低,不利于在CIP道集上正确识别奥陶系内幕地震响应的同相轴信息,利用分频处理方法优选获得高信噪比的低频数据,便于CIP道集同相轴的自动识别和正确追踪剩余延迟,从而更新得到更稳定合理的速度模型(或者背景速度模型)。
步骤S140:在地震属性约束下进行断控高斯束层析,更新深层断裂处的速度模型获得高精度速度模型。
在层位/断层几何特征及其相关地震属性信息的分析研究基础上,利用复杂地质特征约束算子和不等式约束正则化数学建模方法,针对深层断裂进行断控高斯束层析更新深层断裂处的速度模型,生成高精度速度模型。
步骤S150:通过宽频RTM成像技术获得高精度的成像结果。
常规的Kirchhoff积分算法叠前深度偏移技术不能满足横向速度剧烈变化的火成岩复杂区的地震成像,仍然消除不了火成岩对下覆地层造成的假构造,基于双程波波动方程算法的RTM技术更适合火成岩复杂区的地震成像。宽频RTM技术对预设的波场进行成像,可以去除干扰,对碳酸盐岩断溶体获得高精度的成像结果。
其中,所述预设的波场是对成像结果有较大贡献的波场,包括地震波场中的低频部分。
通过图7和图8框线示出的部分可以清楚看出,与常规技术的断溶体成像(图7)相比,本发明的利用“三层一带”速度建模技术和宽频RTM技术的断溶体成像方法提高了断溶体的反射特征,能量更突出,边界更清晰,深层断裂的成像(图8)质量也得到了提高。
图2示出了本发明另一个实施例提供的沙漠地表火成岩复杂区深层碳酸盐岩断溶体地震数据成像方法的流程图,如图2所示,本发明的沙漠地表火成岩复杂区碳酸盐岩断溶体地震数据成像方法,包括:
步骤S210:对地震数据进行预处理,获得保幅的地震道集数据。
其中,所述地震数据通过利用地震勘探方法获得,所述预处理包括低频保护处理和弱信号恢复处理。具体的,所述预处理为在去噪和反褶积等处理环节中尽可能保护更多的低频成分,并针对奥陶系碳酸盐岩地层的弱信号进行有效性号提取和恢复。
其中,所述保幅的地震道集数据包括但不限于共中心点道集。
对于沙漠区深层碳酸盐岩断溶体,由于地表沙漠和浅中层地层对地震信号的吸收衰减作用,地震信号传播到深层后有效信号能量弱,尤其奥陶系碳酸盐岩地层非均质性强,波组特征和异常振幅响应更不明显,需要针对奥陶系碳酸盐目的层进行弱信号提取和恢复。相对于高频的吸收衰减作用更强,低频有效信号对深层地质目标体的识别更有效,所以需要在预处理过程的噪声衰减、反褶积、滤波等关键步骤中注重低频信号的保护,以获得相对保幅的地震道集数据。
步骤S220:利用地质约束下火成岩的精细雕刻技术建立二叠系火成岩地层速度模型。
基于步骤S210产生的高质量的地震道集数据,利用速度建模及其迭代更新方法,建立二叠系火成岩地层速度模型。其中,所述建立二叠系火成岩地层速度模型包括:步骤S221:通过层位拾取建立层位模型;步骤S222:由测井速度约束建立火成岩地层初始速度模型;步骤S223:以1000米间隔进行目标线偏移,根据模型层析建立较为准确的速度模型;步骤S224:针对二叠系火成岩地层,以细网格的间隔偏移数据体进行网格层析速度模型修正。其中,所述细网格的间隔为100m×100m及以下。
由于火成岩的横向边界变化属于一种高频变化,完全数据驱动的网格层析不能完全刻画出岩相变化的边界,需要加入地质的约束。在成像数据体上拾取二叠系火成岩的顶底地层界面,结合地震属性和地质认识,基于不同火成岩岩性的地震相雕刻火成岩岩相变化的边界,在测井数据的约束下,对火成岩地层充填不同的速度信息直至消除火成岩底界明显的扭曲,并且使火成岩底界的CIP道集拉平。
图4A为常规速度建模技术建立的速度模型的剖面图;图4B为本发明地质约束下火成岩精细雕刻技术建立的火成岩高精度速度模型的剖面图;图4C为常规速度建模技术建立的速度模型的切片图;图4D为本发明地质约束下火成岩精细雕刻技术建立的火成岩高精度速度模型的切片图。通过图4A-图4D对比常规速度建模技术建立的速度模型和本发明的地质约束下火成岩精细雕刻技术建立的火成岩高精度速度模型,无论从纵向剖面还是横向切片来看,后者火成岩边界的形态清晰,火成岩层速度模型的精度得到提高。
图5为基于常规速度建模技术的速度模型的火成岩层地震成像剖面,其中,框线示出的部分为火成岩引起的假断裂;图6为本发明的基于地质约束下火成岩精细雕刻技术建立的速度模型的地震成像剖面,通过框线示出的部分可以清楚的看出,本发明消除了火成岩对下覆地层造成的假断裂的影响。
步骤S230:利用分频处理方法得到优势频带的地震数据,建立奥陶系地层和寒武系地层的背景速度模型。
其中,优势频带主要为中低频,地震数据具有高信噪比、高保真度,便于特征识别的特点。
由于低频数据有很好的穿透性,低频信息在深层属于优势频带。由S210的地震道集数据和S220的速度模型数据经过叠前深度偏移,产生CIP道集数据,针对奥陶系的网格层析成像,需要正确拾取CIP道集的剩余延迟信息,而通常奥陶系地层的地震数据信噪比低,不利于在CIP道集上正确识别奥陶系内幕地震响应的同相轴信息,利用分频处理方法获得的由食品袋的地震数据,信噪比高,便于CIP道集同相轴的自动识别和正确追踪剩余延迟,从而更新得到更稳定合理的速度模型(或者背景速度模型)。
步骤S240:利用分频处理方法得到高频数据,基于高频信息对背景速度模型迭代更新获得精细的寒武系地层速度模型。
相对于奥陶系地层,寒武系地层的同相轴波组特征更明显,信噪比更高,提取地震数据的相对高频信息,能显示地震信号的细节信息,方便对上步骤得到的速度模型进行精细的补充和修正,得到更高精度的速度模型。
步骤S250:在地震属性约束下进行断控高斯束层析更新深层断裂处的速度模型获得高精度速度模型。
在层位/断层几何特征及其相关地震属性信息的分析研究基础上,利用复杂地质特征约束算子和不等式约束正则化数学建模方法,针对深层断裂进行断控高斯束层析更新深层断裂处的速度模型。
步骤S260:通过宽频RTM成像技术获得高精度的成像结果。
常规的Kirchhoff积分算法叠前深度偏移技术不能满足横向速度剧烈变化的火成岩复杂区的地震成像,仍然消除不了火成岩对下覆地层造成的假构造,基于双程波波动方程算法的RTM技术更适合火成岩复杂区的地震成像。宽频RTM技术对预设的波场进行成像,可以去除干扰,对碳酸盐岩断溶体获得高精度的成像结果。
其中,所述预设的波场是对成像结果有贡献的波场,包括:地震波场中的低频部分。
通过图7和图8框线示出的部分可以清楚看出,与常规技术的断溶体成像(图7)相比,本发明的利用“三层一带”速度建模技术和宽频RTM技术的断溶体成像方法提高了断溶体的反射特征,能量更突出,边界更清晰,深层断裂的成像(图8)质量也得到了提高。
图3示出了本发明的一种火成岩复杂区碳酸盐岩断溶体地震信号成像服务器,包括:地震道集数据获取模块101,用于对所述地震数据进行预处理,获得保幅的地震道集数据;火成岩地层速度模型建立模块102,用于利用地质约束下所述火成岩的精细雕刻技术建立二叠系火成岩地层速度模型;背景速度模型建立模块103,用于利用分频处理方法得到优势频带的地震数据,建立奥陶系地层和寒武系地层的背景速度模型;高精度速度模型生成模块104,用于在地震属性约束下进行断控高斯束层析更新深层断裂处的速度模型获得高精度速度模型;图像生成模块105,用于通过宽频RTM成像技术获得高精度的成像结果。
优选的,在背景速度模型建立模块103和高精度速度模型生成模块104之间还包括寒武系地层速度速度模型建立模块(未示出),其利用分频处理方法得到高频数据,基于高频信息对背景速度模型迭代更新获得精细的寒武系地层速度模型。
其中,所述服务器还包括地震数据获取模块(未示出),用于利用地震数据处理方法获得所述地震数据。
其中,所述地震道集数据获取模块中所述预处理包括:低频保护处理和弱信号恢复处理。
其中,所述建立二叠系火成岩地层速度模型包括:通过层位拾取建立层位模型;根据时间偏移速度和测井速度建立初始速度模型;以1000米间隔进行目标线偏移,根据模型层析建立较为准确的速度模型;针对二叠系火成岩地层,以细网格的间隔偏移数据体进行网格层析速度模型修正。优选的,所述细网格的间隔为100m×100m及其以下。
此外,本发明提供了一种计算机可读存储介质,其上存储有计算机程序,所述计算机存储介质存储有至少一可执行指令,该计算机可执行指令可执行上述任意方法实施例中的火成岩复杂区碳酸盐岩断溶体地震数据成像方法。
最后,本发明还提供了一种计算机设备,包括存储器,处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现上述任意方法实施例中的火成岩复杂区碳酸盐岩断溶体地震数据成像方法。
根据本发明提供的沙漠地表火成岩复杂区碳酸盐岩断溶体地震数据成像的方法及其设备,在地震数据预处理中,通过低频保护处理和弱信号恢复等预处理技术获得相对保幅的地震道集;在前期叠前时间偏移产生的时间域层速度模型和时间域构造解释成果生成的初始深度域层速度模型基础上,利用地质约束下火成岩精细雕刻技术建立二叠系火成岩地层速度模型;基于保幅地震道集利用分频手段在低频数据上修正更新奥陶系地层速度模型和寒武系地层速度模型,生成相对低频速度模型;以此相对低频速度模型为基础利用优势信号分选技术在高频数据上对其雕刻获得精细的寒武系地层速度模型;在地震属性约束下进行断控高斯束层析获得断裂速度模型;最后,通过宽频RTM成像技术获得高精度的成像结果。本发明消除了火成岩对下覆地层成像的构造假象,提高了深层碳酸盐岩断溶体的成像精度。
在此提供的算法和显示不与任何特定计算机、虚拟系统或者其它设备固有相关。各种通用系统也可以与基于在此的示教一起使用。根据上面的描述,构造这类系统所要求的结构是显而易见的。此外,本发明也不针对任何特定编程语言。应当明白,可以利用各种编程语言实现在此描述的本发明的内容,并且上面对特定语言所做的描述是为了披露本发明的最佳实施方式。
在此处所提供的说明书中,说明了大量具体细节。然而,能够理解,本发明的实施例可以在没有这些具体细节的情况下实践。在一些实例中,并未详细示出公知的方法、结构和技术,以便不模糊对本说明书的理解。
类似地,应当理解,为了精简本发明并帮助理解各个发明方面中的一个或多个特征,在上面对本发明的示例性实施例的描述中,本发明的各个特征有时被一起分组到单个实施例、图、或者对其的描述中。然而,并不应将该公开的方法解释成反映如下意图:即所要求保护的本发明要求比在每个权利要求中所明确记载的特征更多的特征。更确切地说,如下面的权利要求书所反映的那样,发明方面在于少于前面公开的单个实施例的所有特征。因此,遵循具体实施方式的权利要求书由此明确地并入该具体实施方式,其中每个权利要求本身都作为本发明的单独实施例。
本领域那些技术人员可以理解,可以对实施例中的设备中的模块进行自适应性地改变并且把它们设置在与该实施例不同的一个或多个设备中。可以把实施例中的模块或单元或组件组合成一个模块或单元或组件,以及此外可以把它们分成多个子模块或子单元或子组件。除了这样的特征和/或过程或者单元中的至少一些是相互排斥之外,可以采用任何组合对本说明书(包括伴随的权利要求、摘要和附图)中公开的所有特征以及如此公开的任何方法或者设备的所有过程或单元进行组合。除非另外明确陈述,本说明书(包括伴随的权利要求、摘要和附图)中公开的每个特征可以由提供相同、等同或相似目的的替代特征来代替。
此外,本领域的技术人员能够理解,尽管在此的一些实施例包括其它实施例中所包括的某些特征而不是其它特征,但是不同实施例的特征的组合意味着处于本发明的范围之内并且形成不同的实施例。例如,在下面的权利要求书中,所要求保护的实施例的任意之一都可以以任意的组合方式来使用。
应该注意的是上述实施例对本发明进行说明而不是对本发明进行限制,并且本领域技术人员在不脱离所附权利要求的范围的情况下可设计出替换实施例。在权利要求中,不应将位于括号之间的任何参考符号构造成对权利要求的限制。单词“包含”不排除存在未列在权利要求中的元件或步骤。位于元件之前的单词“一”或“一个”不排除存在多个这样的元件。本发明可以借助于包括有若干不同元件的硬件以及借助于适当编程的计算机来实现。在列举了若干装置的单元权利要求中,这些装置中的若干个可以是通过同一个硬件项来具体体现。单词第一、第二、以及第三等的使用不表示任何顺序。可将这些单词解释为名称。
Claims (12)
1.一种火成岩复杂区碳酸盐岩断溶体地震数据成像方法,其特征在于,包括:
对地震数据进行预处理,获得保幅的地震道集数据;
利用地质约束下火成岩的精细雕刻技术建立二叠系火成岩地层速度模型;
利用分频处理方法得到优势频带的地震数据,建立奥陶系地层和寒武系地层的背景速度模型;
在地震属性约束下进行断控高斯束层析,更新深层断裂处的速度模型获得高精度速度模型;
通过宽频RTM成像技术获得高精度的成像结果。
2.根据权利要求1所述的方法,其特征在于,所述预处理包括:低频保护处理和弱信号恢复处理。
3.根据权利要求2所述的方法,其特征在于,所述建立二叠系火成岩地层速度模型包括:
通过层位拾取建立层位模型;
根据时间偏移速度和测井速度建立初始速度模型;
以1000米间隔进行目标线偏移,根据模型层析建立较为准确的速度模型;
针对二叠系火成岩地层,以细网格的间隔偏移数据体进行网格层析速度模型修正。
4.根据权利要求3所述的方法,其特征在于,所述细网格的间隔为100m×100m及以下。
5.根据权利要求1所述的方法,其特征在于,所述建立所述背景速度模型之后,还包括:用分频处理方法得到高频数据,基于高频信息对背景速度模型迭代更新获得精细的寒武系地层速度模型。
6.根据权利要求1所述的方法,其特征在于,所述宽频RTM成像技术选择预设的波场进行成像。
7.一种火成岩复杂区碳酸盐岩断溶体地震信号成像服务器,其特征在于,包括:
地震道集数据获取模块,用于对地震数据进行预处理,获得保幅的地震道集数据;
火成岩地层速度模型建立模块,用于利用地质约束下所述火成岩的精细雕刻技术建立二叠系火成岩地层速度模型;
背景速度模型建立模块,用于利用分频处理方法得到优势频带的地震数据,建立奥陶系地层和寒武系地层的背景速度模型;
高精度速度模型生成模块,用于在地震属性约束下进行断控高斯束层析更新深层断裂处的速度模型获得高精度速度模型;
图像生成模块,用于通过宽频RTM成像技术获得高精度的成像结果。
8.根据权利要求7所述的服务器,其特征在于,所述服务器还包括:
地震数据获取模块,用于利用地震数据处理方法获得所述地震数据。
9.根据权利要求7所述的服务器,其特征在于,所述建立二叠系火成岩地层速度模型包括:
通过层位拾取建立层位模型;
根据时间偏移速度和测井速度建立初始速度模型;
以1000米间隔进行目标线偏移,根据模型层析建立较为准确的速度模型;
针对二叠系火成岩地层,以细网格的间隔偏移数据体进行网格层析速度模型修正。
10.根据权利要求9所述的服务器,其特征在于,所述细网格的间隔为100m×100m及以下。
11.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现如权利要求1-6中任一项所述的方法。
12.一种计算机设备,包括存储器,处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现如权利要求1-6中任一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910149938.2A CN109709605B (zh) | 2019-02-28 | 2019-02-28 | 火成岩复杂区碳酸盐岩断溶体地震数据成像方法及设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910149938.2A CN109709605B (zh) | 2019-02-28 | 2019-02-28 | 火成岩复杂区碳酸盐岩断溶体地震数据成像方法及设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109709605A CN109709605A (zh) | 2019-05-03 |
CN109709605B true CN109709605B (zh) | 2020-09-01 |
Family
ID=66265539
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910149938.2A Active CN109709605B (zh) | 2019-02-28 | 2019-02-28 | 火成岩复杂区碳酸盐岩断溶体地震数据成像方法及设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109709605B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112649862B (zh) * | 2019-10-12 | 2024-06-18 | 中国石油化工股份有限公司 | 基于地层结构信息分离的断溶体识别方法和装置 |
CN111596365B (zh) * | 2020-06-18 | 2023-11-21 | 中国海洋石油集团有限公司 | 针对盐下湖相碳酸盐岩储层段的火山喷发岩地震解释方法 |
CN113466936B (zh) * | 2021-07-01 | 2023-07-18 | 中海石油(中国)有限公司上海分公司 | 一种断层阴影区crp道集获取方法、装置、设备及存储介质 |
CN117908128B (zh) * | 2024-03-08 | 2024-07-30 | 成都星辉科瑞科技有限责任公司 | 一种波数域优势频带识别岩性体断溶体边界的检测方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6278948B1 (en) * | 1999-04-02 | 2001-08-21 | Conoco Inc. | Method for gravity and magnetic data inversion using vector and tensor data |
WO2009152341A2 (en) * | 2008-06-11 | 2009-12-17 | Baker Hughes Incorporated | Vector migration of virtual source vsp data |
CN102977870A (zh) * | 2012-12-27 | 2013-03-20 | 中国石油化工股份有限公司 | 堵漏水泥浆及火成岩裂缝发育井堵漏方法 |
CN104375172A (zh) * | 2013-08-15 | 2015-02-25 | 中国石油集团东方地球物理勘探有限责任公司 | 火山岩下伏地层的构造形态正确成像的方法 |
CN104570125A (zh) * | 2014-09-26 | 2015-04-29 | 郭平 | 一种利用井数据提高成像速度模型精度的方法 |
CN104635269A (zh) * | 2013-11-13 | 2015-05-20 | 中国石油化工股份有限公司 | 一种基于叠前方位角道集预测火成岩裂缝型储层的方法 |
CN108333627A (zh) * | 2018-01-16 | 2018-07-27 | 成都理工大学 | 火成岩地区断裂真伪的识别方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105089657B (zh) * | 2015-06-15 | 2018-05-04 | 中国石油天然气股份有限公司 | 缝洞型碳酸盐岩储层油气充注的物理模拟方法及实验装置 |
-
2019
- 2019-02-28 CN CN201910149938.2A patent/CN109709605B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6278948B1 (en) * | 1999-04-02 | 2001-08-21 | Conoco Inc. | Method for gravity and magnetic data inversion using vector and tensor data |
WO2009152341A2 (en) * | 2008-06-11 | 2009-12-17 | Baker Hughes Incorporated | Vector migration of virtual source vsp data |
CN102977870A (zh) * | 2012-12-27 | 2013-03-20 | 中国石油化工股份有限公司 | 堵漏水泥浆及火成岩裂缝发育井堵漏方法 |
CN104375172A (zh) * | 2013-08-15 | 2015-02-25 | 中国石油集团东方地球物理勘探有限责任公司 | 火山岩下伏地层的构造形态正确成像的方法 |
CN104635269A (zh) * | 2013-11-13 | 2015-05-20 | 中国石油化工股份有限公司 | 一种基于叠前方位角道集预测火成岩裂缝型储层的方法 |
CN104570125A (zh) * | 2014-09-26 | 2015-04-29 | 郭平 | 一种利用井数据提高成像速度模型精度的方法 |
CN108333627A (zh) * | 2018-01-16 | 2018-07-27 | 成都理工大学 | 火成岩地区断裂真伪的识别方法及装置 |
Non-Patent Citations (1)
Title |
---|
"基于多种地震反演方法的哈拉哈塘地区火成岩识别及速度建模";崔永福 等;《东北石油大学学报》;20160831;第40卷(第4期);第54-62页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109709605A (zh) | 2019-05-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109709605B (zh) | 火成岩复杂区碳酸盐岩断溶体地震数据成像方法及设备 | |
CN111123359B (zh) | 随钻测井与地层格架约束的井周地震成像探测方法及装置 | |
Zhao et al. | Key seismic survey technologies for deep complex geological structures: A case study of the northern section of the Longmenshan Fault Fold Belt in the Sichuan Basin | |
CN116520419B (zh) | 一种热流体裂缝通道识别方法 | |
CN110095811B (zh) | 膏岩层速度模型构建处理方法及装置 | |
CN111983683A (zh) | 一种少井条件下盐下湖相灰岩储层预测方法和系统 | |
Vigh et al. | Sparse-node regional study from acquisition to imaging | |
CN103076628A (zh) | 一种孔径优化的叠前时间偏移的处理方法 | |
CN104614762B (zh) | 疏松砂岩气藏边界确定方法及装置 | |
CN105425300B (zh) | 一种剩余静校正方法 | |
CN112147686B (zh) | 多期发育火成岩叠前深度偏移成像速度建模方法及系统 | |
CN109085644A (zh) | 基于双射线走时的真地表成像方法 | |
CN112630841B (zh) | 一种微地震事件检测分析方法 | |
CN105527648A (zh) | 用于各向异性参数反演的敏感度矩阵的计算方法及系统 | |
CN105487120A (zh) | 用于确定山前带的激发井位的方法 | |
CN112031755B (zh) | 一种压裂监测装置、方法、系统、电子设备及存储介质 | |
Melosh et al. | Seismic reflection data and conceptual models for geothermal development in Nevada | |
Huang et al. | Delineating Faults beneath Basalt at the Soda Lake Geothermal Field | |
CN111310361A (zh) | 基于随钻地震的钻井指导方法、系统、设备及存储介质 | |
Walmsley et al. | Electromagnetic borehole mapping: North Sea case study | |
Zarubov* et al. | Optimal near-surface characterization—Case studies from onshore Abu Dhabi, UAE | |
CN113514882B (zh) | 一种地震勘探数据的采集方法、装置、设备及介质 | |
Zhou et al. | Anisotropic model building with well control | |
Maili et al. | Maximizing the value of existing seismic data in the Bahrain Field with wide-azimuth processing technology | |
CN110941029B (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 |