CN117950021A - 一种超深层碳酸盐岩断溶体储层地震相控反演方法及装置 - Google Patents

一种超深层碳酸盐岩断溶体储层地震相控反演方法及装置 Download PDF

Info

Publication number
CN117950021A
CN117950021A CN202211296792.2A CN202211296792A CN117950021A CN 117950021 A CN117950021 A CN 117950021A CN 202211296792 A CN202211296792 A CN 202211296792A CN 117950021 A CN117950021 A CN 117950021A
Authority
CN
China
Prior art keywords
seismic
inversion
model
phase control
initial model
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
Application number
CN202211296792.2A
Other languages
English (en)
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.)
Sinopec Petroleum Geophysical Exploration Technology Research Institute Co ltd
China Petroleum and Chemical Corp
Original Assignee
Sinopec Petroleum Geophysical Exploration Technology Research Institute Co ltd
China Petroleum and Chemical Corp
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 Sinopec Petroleum Geophysical Exploration Technology Research Institute Co ltd, China Petroleum and Chemical Corp filed Critical Sinopec Petroleum Geophysical Exploration Technology Research Institute Co ltd
Priority to CN202211296792.2A priority Critical patent/CN117950021A/zh
Publication of CN117950021A publication Critical patent/CN117950021A/zh
Pending legal-status Critical Current

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. for interpretation or for event detection
    • G01V1/282Application of seismic models, synthetic seismograms
    • 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. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/306Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface
    • G01V2210/624Reservoir parameters

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开一种超深层碳酸盐岩断溶体储层地震相控反演方法及装置,该方法包括:采用测井数据、地震层位和地震偏移速度体构建混合速度模型;基于混合速度模型构建反演初始模型;采用地震张量和相干能量作为表征断裂带和溶洞的先验地震相信息,并提取三维断裂目标作为相控建模的空间边界约束;利用线性弱相控约束建模的方式,向反演初始模型依次加入先验地震相信息,对反演初始模型更新迭代,获取相控约束反演初始模型;基于相控约束反演初始模型进行约束稀疏脉冲反演,获得地震相控反演结果。本发明实现了对断溶体储层的高精度地震相控约束反演,可有效提升断溶体储层识别刻画的精度和可靠性。

Description

一种超深层碳酸盐岩断溶体储层地震相控反演方法及装置
技术领域
本发明属于地震反演技术领域,更具体地,涉及一种超深层碳酸盐岩断溶体储层地震相控反演方法及装置。
背景技术
波阻抗是油气勘探领域储层识别刻画、物性预测、储量估算的重要参数。在实际生产过程中,广泛采用地震波阻抗反演技术获取地下目标地质体的阻抗参数。地震波阻抗反演可以综合利用地震、测井、地质等资料,计算出地质体的阻抗参数。然而,对于超深层碳酸盐岩油气藏断溶体储层的识别,常规地震波阻抗反演技术在实际生产应用中存在较大的局限性。这是因为,该类储层埋深较大、非均质性强,断控岩溶作用较为复杂,断裂带内部多为杂乱反射,断溶体有效储层的振幅异常不够突出,常规波阻抗反演方法难以有效预测出该类储层,影响储层识别刻画和储量计算的精度,给井位部署和开发带来较大的风险。
针对超深碳酸盐岩油气藏断溶体储层的反演识别问题,很多学者开展了相关的研究以提升储层识别的精度,包括采用高分辨率地震成像资料、高精度地震反演、地震相控约束反演等技术。其中,地震相控约束反演是在实际生产中提升此类储层识别精度较为有效且实用的技术手段。常规波阻抗反演方法一般是通过融合测井插值初始模型的低频信息和从地震数据中提取出的高频信息,获得地层的绝对波阻抗。但是,这种由测井插值的初始模型往往难以表达超深碳酸盐岩的断控岩溶特征,从而降低了这类储层反演的精度。地震相控反演技术主要是通过在初始模型中增加可以表征特定勘探目标的地质、地震等先验信息,可有效提升对目标地质体的反演识别效果。需要指出的是,当前实际生产中所采用的地震相控约束反演技术仍然面临很大的挑战,主要有三个方面的原因:
(1)常规测井插值的初始模型建模方法难以适应于超深碳酸盐岩储层。与碎屑岩油气藏不同的是,超深碳酸盐岩油气藏断溶体储层段常缺乏足够深的声波时差、密度等测井曲线。此外,断控岩溶储层的纵横向非均质性较强,导致常规测井插值建模方法难以适应于碳酸盐储层,很容易造成低频模型中出现条带状假异常。尽管有的学者尝试利用地震偏移速度构建初始模型,以解决井插值建模存在的问题,但这种方法建立的初始模型存在纵向分辨率低的问题。
(2)在初始模型的基础上采用全区整体相控建模会破坏反演结果的可靠性。常规相控约束方法是往整个初始模型中添加断裂和溶洞的先验地震相信息,以迭代更新初始模型,这种整体添加地震相的方式会破坏非断裂带初始模型的可靠性,造成反演结果失真。
(3)对地震相采用强相控约束的初始建模方式,容易造成反演结果存在明显的边界效应。常规相控约束的初始建模方法一般给不同的地震相赋予相应的固定阻抗值。这种赋值方式的相控约束性过强,导致反演结果出现严重的横向边界效应,造成反演结果失真。
发明内容
有鉴于此,本发明的目的提供一种超深层碳酸盐岩断溶体储层地震相控反演方法及装置,解决反演结果易失真的技术问题,提升断溶体储层识别刻画的精度。
为实现上述目的,一方面,本发明提供一种超深层碳酸盐岩断溶体储层地震相控反演方法,包括:
S1:采用测井数据、地震层位和地震偏移速度体构建混合速度模型;
S2:基于所述混合速度模型构建反演初始模型;
S3:采用地震张量和相干能量作为表征断裂带和溶洞的先验地震相信息,并提取三维断裂目标作为相控建模的空间边界约束;
S4:利用线性弱相控约束建模的方式,向所述反演初始模型依次加入所述断裂带和内部溶洞的所述先验地震相信息,对所述反演初始模型更新迭代,获取相控约束反演初始模型;
S5:基于所述相控约束反演初始模型进行约束稀疏脉冲反演,获得地震相控反演结果。
优选地,步骤S1中,构建的所述混合速度模型的公式为:
Vmix=f(Dtop,Dtd,Dhor,Vwell,Vseis)
式中,Vmix为混合速度模型,Dtop为测井分层数据,Dtd为标定的测井时深数据,Dhor为地震层位数据,Vwell为测井速度数据,vseis为地震偏移速度体。
优选地,步骤S2中,构建的所述反演初始模型的公式为:
PImodel=g(Vmix)
式中,PImodel为反演初始模型,Vmix为混合速度模型。
优选地,步骤S3中,所述三维断裂目标提取的方法包括:
S31:获取整个地震数据体的结构张量属性;
S32:沿目的层提取大时窗的结构张量平面属性;
S33:根据所述结构张量平面属性的断裂带分布特征,完成三维断裂目标的边界标注和空间提取。
优选地,步骤S4中,获取所述相控约束初始模型的公式为:
式中,PImodel为相控约束初始模型,E为地震属性值,Emax为地震属性的最大值,Emin为地震属性的最小值,PImax为通过测井统计的阻抗最大值,PImin通过测井统计的阻抗最小值。
另一方面,本发明还提供一种超深层碳酸盐岩断溶体储层地震相控反演装置,包括:
混合速度模型建立模块,用于采用测井数据、地震层位和地震偏移速度体构建混合速度模型;
反演初始模型建立模块,用于基于所述混合速度模型构建反演初始模型;
提取模块,用于采用地震张量和相干能量作为表征断裂带和溶洞的先验地震相信息,并提取三维断裂目标作为相控建模的空间边界约束;
相控约束反演初始模型获取模块,用于利用线性弱相控约束建模的方式,向所述反演初始模型依次加入所述断裂带和内部溶洞的所述先验地震相信息,对所述反演初始模型更新迭代,获取相控约束反演初始模型;
地震相控反演结果获取模块,用于基于所述相控约束反演初始模型进行约束稀疏脉冲反演,获得地震相控反演结果。
优选地,构建的所述混合速度模型的公式为:
Vmix=f(Dtop,Dtd,Dhor,Vwell,Vseis)
式中,Vmix为混合速度模型,Dtop为测井分层数据,Dtd为标定的测井时深数据,Dhor为地震层位数据,Vwell为测井速度数据,vseis为地震偏移速度体。
优选地,构建的所述反演初始模型的公式为:
PImodel=g(Vmix)
式中,PImodel为反演初始模型,Vmix为混合速度模型。
优选地,所述三维断裂目标提取的方法包括:
获取整个地震数据体的结构张量属性;
沿目的层提取大时窗的结构张量平面属性;
根据所述结构张量平面属性的断裂带分布特征,完成三维断裂目标的边界标注和空间提取。
优选地,获取所述相控约束反演初始模型的公式为:
式中,PImodel为相控约束初始模型,E为地震属性值,Emax为地震属性的最大值,Emin为地震属性的最小值,PImax为通过测井统计的阻抗最大值,PImin通过测井统计的阻抗最小值。
本发明的有益效果在于:
本发明通过采用基于混合速度模型构建的反演初始模型、三维断裂目标提取以及弱相控约束建模技术,在反演过程中加入与断裂带、溶洞有关的先验地震相信息,实现对断溶体储层的高精度地震波阻抗反演,有效提升断溶体储层识别刻画的精度和可靠性。
本发明的其它特征和优点将在随后具体实施方式部分予以详细说明。
附图说明
通过结合附图对本发明示例性实施方式进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显,其中,在本发明示例性实施方式中,相同的参考标号通常代表相同部件。
图1示出了本发明一种超深碳酸盐岩断溶体储层地震相控反演方法的流程图。
图2示出了本发明的一种典型的超深碳酸盐岩断溶体储层地震剖面图。
图3示出了本发明采用混合速度建模构建反演初始模型的示意图。。
图4示出了本发明的三维断裂目标提取示意图。
图5示出了本发明的线性弱相控约束下的反演初始模型示意图。
图6示出了顺北某工区断溶体储层反演结果剖面对比示意图。
图7示出了托甫台某井区单元断溶体储层反演结果剖面对比示意图。
图8示出了托甫台某井区单元断溶体储层反演结果平面对比示意图。
具体实施方式
为具体说明本发明的技术方案,本发明通过以下示例性实施方式进行更详细的描述。
实施例1
参考图1,本发明提供一种超深层碳酸盐岩断溶体储层地震相控反演方法,包括:
S1:采用测井数据、地震层位和地震偏移速度体构建混合速度模型;
S2:基于所述混合速度模型构建反演初始模型;
S3:采用地震张量和相干能量作为表征断裂带和溶洞的先验地震相信息,并提取三维断裂目标作为相控建模的空间边界约束;
S4:利用线性弱相控约束建模的方式,向所述反演初始模型依次加入所述断裂带和内部溶洞的所述先验地震相信息,对所述反演初始模型更新迭代,获取相控约束反演初始模型;
S5:基于所述相控约束反演初始模型进行约束稀疏脉冲反演,获得地震相控反演结果。
以下对上述各步骤做进一步的具体说明。
执行步骤S1和S2,采用测井数据、地震层位、地震偏移速度体构建混合速度模型;基于所述混合速度模型构建反演初始模型。
本实施例中,断溶体储层的非均质性较强,目的层段缺少深部的测井曲线,采用常规井插方法构建的初始模型,虽然纵向分辨率高,但是横向约束能力较差,很容易出现条带状假异常,且需要人为补充深部的初始模型数据。实际上,地震偏移速度体具有良好的横向约束能力,但是纵向上的分辨率不够高。本技术方案采用测井数据、地震层位和地震偏移速度体构建混合速度模型,再将混合速度模型转换成反演初始模型,这种方式既能充分利用井资料纵向分辨率高的优势,也可以利用地震偏移速度体的横向约束能力,最终获取较为准确的反演初始模型。在构建混合速度模型的过程中涉及到测井数据、时深数据、地震层位和地震偏移速度体等数据。
本实施例中,构建混合速度模型的公式为:
Vmix=f(Dtop,Dtd,Dhor,Vwell,Vseis)
式中,Vmix为混合速度模型,Dtop为测井分层数据,Dtd为标定的测井时深数据,Dhor,为地震层位数据,Vwell为测井速度数据,Vseis为地震偏移速度体。
本实施例中,构建反演初始模型的公式为:
PImodel=g(Vmix)
式中,PImodel为反演初始模型,Vmix为混合速度模型。
执行步骤S3,采用地震张量和相干能量作为表征断裂带和溶洞的先验地震相信息,并提取三维断裂目标作为相控建模的空间边界约束。
本实施例中,采用地震属性优选及三维断裂目标提取的方法,优选出可以表征断裂带和溶洞的地震属性作为先验地震相信息,并通过提取三维断裂目标作为相控建模的空间边界约束。通过这种方式,既能针对性的突出断裂带内部的储层信息,又不破坏断裂带外部的初始模型的特征。这一步包含地震属性优选、3D断裂目标提取;地震属性优选主要是优选能够有效表征断裂、溶洞的地震属性作为地震相,经过反复对比测试,本实施例采用结构张量和相干能量梯度两种地震属性分别表征断裂和溶洞信息。3D断裂目标提取是非常关键的一步,整个相控反演的核心就是围绕提取的断裂目标进行的。提取过程如下:
S31:计算整个地震数据体的结构张量属性;
S32:沿目的层提取大时窗的结构张量平面属性;
S33:根据所述结构张量平面属性的断裂带分布特征,完成三维断裂目标的边界标注和空间提取。
执行步骤S4和S5,利用线性弱相控约束建模的方式,向所述反演初始模型依次加入所述断裂带和内部溶洞的所述先验地震相信息,对所述反演初始模型更新迭代,获取相控约束反演初始模型;基于所述相控约束反演初始模型进行约束稀疏脉冲反演,获得地震相控反演结果。
本实施例中,采用线性弱相控约束建模的方式构建反演所需的初始模型。把提取的三维断裂目标作为空间边界约束,向反演初始模型依次加入断裂带和内部溶洞的先验地震相信息,实现对反演初始模型的更新迭代,建立最终的相控约束反演初始模型。常规相控约束的方法是给每种地震相赋予单一的阻抗值以更新初始模型,但这种强相控的方式会导致反演结果出现较为严重的边界效应。本实施例采用线性赋值的弱相控方式对每个地震相进行赋值,首先进行测井统计分析,在井震标定的基础上,研究不同地质目标的阻抗值及表征地震相属性值的分布规律,为每种地质目标建立阻抗与地震属性之间的线性关系。根据线性关系赋予该地震相不同的阻抗值,以更新迭代反演初始模型,获取相控约束初始模型,相控约束初始模型可用公式表示如下:
式中,PImodel为相控约束初始模型,E为地震属性值,Emax为地震属性的最大值,Emin为地震属性的最小值,PImax为通过测井统计的阻抗最大值,PImin通过测井统计的阻抗最小值。
本实施例通过统计分析建立阻抗与地震属性的线性关系,采用线性赋值的方式添加地震相信息,实现初始模型进行更新迭代,最终建立弱相控约束下的反演初始模型,可以有效避免常规相控反演存在的边界效应和反演结果失真的问题。在实际应用中可根据统计分析结果对地震属性的最大值和最小值进行截断,以避免异常极端值的影响。在此基础上,对添加地震相信息后的相控约束初始模型进行空间平滑,再次降低地震相边界约束的强度,避免反演结果出现横向边界效应,以提高断溶体储层地震反演的精度。
在相控约束反演初始模型构建好的基础上执行约束稀疏脉冲反演,即可获得最终的地震相控反演结果。
实施例2
在一个具体应用场景中,将实施例1的方法应用于塔河油田顺北某工区超深层碳酸盐岩断溶体储层的识别预测,该地区典型的断溶体储层地震响应特征如图2中的阴影范围所示,可以看出该地区断溶体埋深较大,断裂带内储层的地震响应特征与围岩的振幅差异不太明显。分别按照图3、图4和图5所示的技术路线,依次完成基于混合速度模型的反演初始模型建模、三维断裂目标提取、线性弱相控约束建模,最终构建相控约束下的反演初始模型,并完成地震波阻抗反演;其中,图3中的图3(a)为测井资料:波阻抗曲线示意图;图3(b)为井震标定示意图;图3(c)为地震层位示意图;图3(d)为地震偏移速度体示意图;图3(e)为混合速度模型构建的反演初始模型;图4中的图4(a)为原始地震数据示意图;图4(b)为提取的三维断裂目标体示意图;图5中的图5(a)为混合速度模型构建的反演初始模型示意图;图5(b)为张量属性示意图;图5(c)为相干能量梯度属性示意图;图5(d)为线性弱相控约束下的反演初始模型示意图;为了测试本发明提出的反演方法有效性,分别采用了约束稀疏脉冲反演、常规相控约束反演进行对比,结果如图6所示,其中图6中的图6(a)为地震数据图示意图,图6(b)为约束稀疏脉冲反演结果示意图,图6(c)为常规相控约束反演结果示意图,图6(d)为本发明反演方法的结果示意图。可以看出,由于断裂带内部储层的振幅异常不够明显,导致常规约束稀疏脉冲反演难以识别出该断裂带内部低阻抗的有效储层。采用常规相控约束反演方法可一定程度上识别出断裂带内储层,但由于这种方法的相控约束作用过强,且是对整个工区不区分进行相控约束建模,导致在反演剖面上断裂带内部和外部均出现严重的竖条状边界效应,造成反演结果失真,影响储层的识别刻画。本实施例提出的反演方法不仅有效识别出了断裂带内部的储层,消除了常规相控反演方法的边界效应,反演结果更合理,而且对于断裂带外部的反演,则充分利用了地震数据的保幅性,反演结果较可靠。
实施例3
在另一个具体应用场景中,将实施例1的方法用于塔河油田托甫台某井区单元深层碳酸盐岩断溶体储层的识别预测,该井区单元断裂带内部的溶蚀作用更强。图7为常规反演与本发明反演方法的剖面对比示意图。其中,图7中的图7(a)为地震数据示意图,图7(b)为常规反演结果示意图,从图7(b)中可以看出,常规波阻抗反演只能识别断裂带浅部地震响应较为明显的部分溶洞,但难以识别出断裂带中深部反射特征相对较弱的溶洞储层,这将严重影响储层刻画与储量计算的精度,图7(c)为本方案提出的线性弱相控反演结果示意图,从图7(c)中可以看出,线性弱相控约束的反演结果不仅可以识别出断裂带内的溶洞储层,对断裂带中深部的储层也能有效的识别刻画,反演结果也更符合断溶体储层的地质特征。图8为实际井区单元断溶体储层反演结果的平面图对比示意图。其中,图8中的8(a)为目的层向下10-100ms均方根振幅示意图;8(b)为常规反演结果目的层向下10-100ms最小振幅属性示意图;8(c)为本方案提出的线性弱相控反演结果目的层向下10-100ms最小振幅属性示意图;从图8中的8(a)中可以看出,目的层下部10-100ms的地震数据均方根振幅属性能基本描述整个井区单元浅部的缝洞型储集体分布特征,但对单元内F1、F2两条呈X型的走滑断裂带内的断控岩溶储集体识别效果较差。实际上,位于断裂带上的井大多具有较好的产能,生产数据表明,仅左侧北东向断裂带F1上11口井的累产油已达到42万吨。从图8(b)中可以看出,常规反演方法能基本识别出整个井区单元浅部的缝洞型储集体,识别效果比振幅属性更加精细,但同样对断裂带内部的储层难以表征,导致由此刻画出的断溶体储层明显偏低,与实际生产情况不吻合。需要指出的是,该图也能在一定程度上说明对于断裂带外部的缝洞型储集体,常规方法一般能够有效识别,因此对于断裂带外部应谨慎使用相控约束。从图8(c)中可以看出,利用本发明所提出的相控约束反演方法,不仅可以识别出整个井区单元的缝洞型储集体分布特征,而且对于断裂带内的储集体同样可以有效识别,特别是在两条断裂带的交汇处,尽管此处断裂带内地层更加破碎,溶蚀作用更强,地震反射特征更为复杂,该方法仍然能够识别刻画出断裂带内部的大量有效储层,这与井的实际生产情况更加吻合。由于该方法只对断裂带内部进行相控约束,因此不会对断裂带外部的反演结果造成影响,反演结果更可靠。
实施例4
本发明提供一种超深层碳酸盐岩断溶体储层地震相控反演装置,包括:
混合速度模型建立模块,用于采用测井数据、地震层位和地震偏移速度体构建混合速度模型;
反演初始模型建立模块,用于基于所述混合速度模型构建反演初始模型;
提取模块,用于采用地震张量和相干能量作为表征断裂带和溶洞的先验地震相信息,并提取三维断裂目标作为相控建模的空间边界约束;
相控约束反演初始模型获取模块,用于利用线性弱相控约束建模的方式,向所述反演初始模型依次加入所述断裂带和内部溶洞的所述先验地震相信息,对所述反演初始模型更新迭代,获取相控约束反演初始模型;
地震相控反演结果获取模块,用于基于所述相控约束反演初始模型进行约束稀疏脉冲反演,获得地震相控反演结果。
本实施例中,构建的所述混合速度模型的公式为:
Vmix=f(Dtop,Dtd,Dhor,Vwell,Vseis)
式中,Vmix为混合速度模型,Dtop为测井分层数据,Dtd为标定的测井时深数据,Dhor为地震层位数据,vwell为测井速度数据,Vseis为地震偏移速度体。
本实施例中,构建的所述反演初始模型的公式为:
PImodel=g(Vmix)
式中,PImodel为反演初始模型,Vmix为混合速度模型。
本实施例中,所述三维断裂目标提取的方法包括:
获取整个地震数据体的结构张量属性;
沿目的层提取大时窗的结构张量平面属性;
根据所述结构张量平面属性的断裂带分布特征,完成三维断裂目标的边界标注和空间提取。
本实施例中,获取所述相控约束反演初始模型的公式为:
式中,PImodel为相控约束初始模型,E为地震属性值,Emax为地震属性的最大值,Emin为地震属性的最小值,PImax为通过测井统计的阻抗最大值,PImin通过测井统计的阻抗最小值。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。

Claims (10)

1.一种超深层碳酸盐岩断溶体储层地震相控反演方法,其特征在于,包括:
S1:采用测井数据、地震层位和地震偏移速度体构建混合速度模型;
S2:基于所述混合速度模型构建反演初始模型;
S3:采用地震张量和相干能量作为表征断裂带和溶洞的先验地震相信息,并提取三维断裂目标作为相控建模的空间边界约束;
S4:利用线性弱相控约束建模的方式,向所述反演初始模型依次加入所述断裂带和内部溶洞的所述先验地震相信息,对所述反演初始模型更新迭代,获取相控约束反演初始模型;
S5:基于所述相控约束反演初始模型进行约束稀疏脉冲反演,获得地震相控反演结果。
2.根据权利要求1所述的地震相控反演方法,其特征在于,步骤S1中,构建的所述混合速度模型的公式为:
Vmix=f(Dtop,Dtd,Dhor,Vwell,Vseis)
式中,Vmix为混合速度模型,Dtop为测井分层数据,Dtd为标定的测井时深数据,Dhor为地震层位数据,Vwell为测井速度数据,Vseis为地震偏移速度体。
3.根据权利要求1所述的地震相控反演方法,其特征在于,步骤S2中,构建的所述反演初始模型的公式为:
PImodel=g(Vmix)
式中,PImodel为反演初始模型,Vmix为混合速度模型。
4.根据权利要求1所述的地震相控反演方法,其特征在于,步骤S3中,所述三维断裂目标提取的方法包括:
S31:获取整个地震数据体的结构张量属性;
S32:沿目的层提取大时窗的结构张量平面属性;
S33:根据所述结构张量平面属性的断裂带分布特征,完成三维断裂目标的边界标注和空间提取。
5.根据权利要求1所述的地震相控反演方法,其特征在于,步骤S4中,获取所述相控约束反演初始模型的公式为:
式中,PImodel为相控约束初始模型,E为地震属性值,Emax为地震属性的最大值,Emin为地震属性的最小值,PImax为通过测井统计的阻抗最大值,PImin通过测井统计的阻抗最小值。
6.一种超深层碳酸盐岩断溶体储层地震相控反演装置,其特征在于,包括:
混合速度模型建立模块,用于采用测井数据、地震层位和地震偏移速度体构建混合速度模型;
反演初始模型建立模块,用于基于所述混合速度模型构建反演初始模型;
提取模块,用于采用地震张量和相干能量作为表征断裂带和溶洞的先验地震相信息,并提取三维断裂目标作为相控建模的空间边界约束;
相控约束反演初始模型获取模块,用于利用线性弱相控约束建模的方式,向所述反演初始模型依次加入所述断裂带和内部溶洞的所述先验地震相信息,对所述反演初始模型更新迭代,获取相控约束反演初始模型;
地震相控反演结果获取模块,用于基于所述相控约束反演初始模型进行约束稀疏脉冲反演,获得地震相控反演结果。
7.根据权利要求6所述的地震相控反演装置,其特征在于,构建的所述混合速度模型的公式为:
Vmix=f(Dtop,Dtd,Dhor,Vwell,Vseis)
式中,Vmix为混合速度模型,Dtop为测井分层数据,Dtd为标定的测井时深数据,Dhor,为地震层位数据,Vwell为测井速度数据,Vseis为地震偏移速度体。
8.根据权利要求6所述的地震相控反演装置,其特征在于,构建的所述反演初始模型的公式为:
PImodel=g(Vmix)
式中,PImodel为反演初始模型,Vmix为混合速度模型。
9.根据权利要求6所述的地震相控反演装置,其特征在于,所述三维断裂目标提取的方法包括:
获取整个地震数据体的结构张量属性;
沿目的层提取大时窗的结构张量平面属性;
根据所述结构张量平面属性的断裂带分布特征,完成三维断裂目标的边界标注和空间提取。
10.根据权利要求6所述的地震相控反演装置,其特征在于,获取所述相控约束反演初始模型的公式为:
式中,PImodel为相控约束初始模型,E为地震属性值,Emax为地震属性的最大值,Emin为地震属性的最小值,PImax为通过测井统计的阻抗最大值,PImin通过测井统计的阻抗最小值。
CN202211296792.2A 2022-10-21 2022-10-21 一种超深层碳酸盐岩断溶体储层地震相控反演方法及装置 Pending CN117950021A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211296792.2A CN117950021A (zh) 2022-10-21 2022-10-21 一种超深层碳酸盐岩断溶体储层地震相控反演方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211296792.2A CN117950021A (zh) 2022-10-21 2022-10-21 一种超深层碳酸盐岩断溶体储层地震相控反演方法及装置

Publications (1)

Publication Number Publication Date
CN117950021A true CN117950021A (zh) 2024-04-30

Family

ID=90794884

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211296792.2A Pending CN117950021A (zh) 2022-10-21 2022-10-21 一种超深层碳酸盐岩断溶体储层地震相控反演方法及装置

Country Status (1)

Country Link
CN (1) CN117950021A (zh)

Similar Documents

Publication Publication Date Title
US10386531B2 (en) Geological model analysis incorporating cross-well electromagnetic measurements
CN104516018B (zh) 一种地球物理勘探中岩性约束下的孔隙度反演方法
US20180231681A1 (en) Reservoir resistivity characterization incorporating flow dynamics
CN105005097A (zh) 利用重力、磁力、电磁、地震资料综合识别火成岩方法
CN107219553B (zh) 基于gr分频反演的暗河充填预测方法
EP3253946B1 (en) Seismic attributes derived from the relative geological age property of a volume-based model
CN113296166B (zh) 裂缝模型的构建方法
CN109799540B (zh) 基于地质信息约束下的火山岩型铀矿床磁化率反演方法
CN112698392A (zh) 一种多层系储集体立体刻画与油气空间分布及定容方法
CN115327663B (zh) 深部矿产资源勘查用空-地-井立体地球物理探测方法
CN107861162B (zh) 基于微电极测井资料的天然裂缝识别方法及系统
Ray et al. Characterizing and modeling natural fracture networks in a tight carbonate reservoir in the Middle East: A methodology
CN107515430B (zh) 一种地震法探测盐湖卤水的方法
CN113779811B (zh) 一种构造变形对现今地应力方位扰动的定量分析方法
RU2757386C1 (ru) Способ проведения электромагнитного мониторинга ГРП
US11703612B2 (en) Methods and systems for characterizing a hydrocarbon-bearing rock formation using electromagnetic measurements
Zhang et al. Architecture characteristics and characterization methods of fault-controlled karst reservoirs: A case study of the Shunbei 5 fault zone in the Tarim Basin, China
CN117950021A (zh) 一种超深层碳酸盐岩断溶体储层地震相控反演方法及装置
CN116068663A (zh) 基于磁震联合低频建模的火成岩波阻抗反演方法
Li et al. Heterogeneous reservoir prediction of ultra-deep strike-slip fault-damaged zone constrained with local seismic anomaly data
CN115577616A (zh) 基于深度学习的碳酸岩缝洞地震刻画方法及装置
CN115598732A (zh) 基于重磁电震勘探技术融合的潜山构造识别方法
Saliés Study on the feasibility of using electromagnetic methods for fracture diagnostics
CN110927819B (zh) 一种裂缝发育程度表征方法
CN113514884A (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