CN113552622B - 断层封堵性分析方法及装置 - Google Patents

断层封堵性分析方法及装置 Download PDF

Info

Publication number
CN113552622B
CN113552622B CN202010331935.3A CN202010331935A CN113552622B CN 113552622 B CN113552622 B CN 113552622B CN 202010331935 A CN202010331935 A CN 202010331935A CN 113552622 B CN113552622 B CN 113552622B
Authority
CN
China
Prior art keywords
fault
work area
layer
mudstone
thickness
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
CN202010331935.3A
Other languages
English (en)
Other versions
CN113552622A (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.)
Petrochina Co Ltd
Original Assignee
Petrochina Co 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 Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN202010331935.3A priority Critical patent/CN113552622B/zh
Publication of CN113552622A publication Critical patent/CN113552622A/zh
Application granted granted Critical
Publication of CN113552622B publication Critical patent/CN113552622B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/282—Application of seismic models, synthetic seismograms
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/301—Analysis for determining seismic cross-sections or geostructures
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/61—Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616—Data from specific type of measurement
    • G01V2210/6169—Data from specific type of measurement using well-logging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G01V2210/622—Velocity, density or impedance
    • G01V2210/6226—Impedance
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/64—Geostructures, e.g. in 3D data cubes
    • G01V2210/642—Faults
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/66—Subsurface modeling

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

断层封堵性分析方法及装置
技术领域
本发明涉及地质研究技术领域，特别涉及一种断层封堵性分析方法及装置。
背景技术
断层的演化过程分为活动期与静止期，活动期一般都是短暂的，断层主要为静止期。断层形成期或活动期为开启状态，可以作为油气运移的通道。断层封闭性和启闭性能评价主要针对处于静止期的断层，确定断层对油气是启闭还是封闭。断层封闭性评价历经了定性与定量两个阶段，20世纪50年代后主要通过分析断层封闭的形成机理，提出了砂泥对接、泥岩涂抹、成岩特征等均可产生封闭条件，使用“封闭”与“不封闭”定性描述断层的封闭能力；断层启闭性定量地分析方法的提出和应用始于上世纪80年代末，目前主要评价方法包括多层泥岩涂抹系数(SGR，Shale Gouge Ratio)、页岩涂抹系数(SSF，Shale SmearFactor)、Allen图解法上下盘关系、断面正应力等。
虽然目前针对断层平面封堵性定量地分析的方法有很多，但是存在一定局限性，如页岩涂抹系数(SSF)方法，仅考虑了断距和泥岩层厚度，没有考虑不同位置泥岩涂抹量不同，且页岩涂抹系数(SSF)的计算公式过于简单化，无法真正反映泥岩沿断层的纵向涂抹情况；再如Allen图解法，分析上下盘关系时，需要对断层上下盘的构造形态进行恢复，在构造复杂的情况下比较困难。断面正应力需要分析计算断面处的地层压力、断层倾角和垂直应力等，是比较困难的。而多层泥岩涂抹评价方法(SGR)是目前接受程度最高、使用最为广泛的断层封堵性评价方法，但目前的方法多集中与单剖面或者多个单剖面的评价中，通常直接利用自然伽马曲线进行泥岩含量反演，不能准确地对断层封堵性进行定量地分析。
综上所述，虽然目前针对断层封堵性定量地分析的方法有很多，但并未存在一种可准确地对断层封堵性进行分析的方法。
发明内容
本发明实施例提供了一种断层封堵性分析方法，用以准确地对断层封堵性进行定量分析，该方法包括：
根据预反演的工区相对阻抗数据体，建立工区三维构造模型；
根据工区三维构造模型，对单井岩相数据进行属性建模，得到工区相控数据体；
在工区相控数据体约束下，对泥质含量数据进行地质统计学建模，得到泥质含量模型；
根据泥质含量模型，计算出工区断层的泥岩涂抹系数；
结合区域地质数据，根据工区断层的泥岩涂抹系数，对断层封堵性进行分析。
本发明实施例还提供了一种断层封堵性分析装置，用以准确地对断层封堵性进行定量分析，该装置包括：
工区三维构造模型建立模块，用于根据预反演的工区相对阻抗数据体，建立工区三维构造模型；
相控数据体建模模块，用于根据工区三维构造模型，对单井岩相数据进行属性建模，得到工区相控数据体；
泥质含量模型建模模块，用于在工区相控数据体约束下，对泥质含量数据进行地质统计学建模，得到泥质含量模型；
泥岩涂抹系数计算模块，用于根据泥质含量模型，计算出工区断层的泥岩涂抹系数；
断层封堵性分析模块，用于结合区域地质数据，根据工区断层的泥岩涂抹系数，对断层封堵性进行分析。
本发明实施例还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述所述方法。
本发明实施例还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述所述方法的计算机程序。
本发明实施例中，根据预反演的工区相对阻抗数据体，建立工区三维构造模型；根据工区三维构造模型，对单井岩相数据进行属性建模，得到工区相控数据体；在工区相控数据体约束下，对泥质含量数据进行地质统计学建模，得到泥质含量模型；根据泥质含量模型，计算出工区断层的泥岩涂抹系数；结合区域地质数据，根据工区断层的泥岩涂抹系数，对断层封堵性进行分析，从而在工区相控数据体约束下，建模得到泥质含量模型，与现有技术对比，对泥岩分布有了以相控数据体实现的地质约束，可更加真实地以泥质含量模型反映泥岩的实际展布，可准确地对断层封堵性进行定量地分析。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种断层封堵性分析方法的流程示意图；
图2是本发明实施例提供的一种断层封堵性分析方法的实例的流程示意图；
图3是本发明实施例提供的一种断层封堵性分析方法的流程示意图；
图4是本发明实施例提供的一种断层封堵性分析方法的流程示意图；
图5是本发明实施例提供的一种断层封堵性分析装置的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
目前针对断层封堵性定量地分析的方法有很多，但并未存在一种可准确地对断层封堵性进行分析的方法。而在本发明实施例中，提供了一种断层封堵性分析方法，用以准确地对断层封堵性进行定量分析，如图1所示，该方法包括：
步骤101：根据预反演的工区相对阻抗数据体，建立工区三维构造模型；
步骤102：根据工区三维构造模型，对单井岩相数据进行属性建模，得到工区相控数据体；
步骤103：在工区相控数据体约束下，对泥质含量数据进行地质统计学建模，得到泥质含量模型；
步骤104：根据泥质含量模型，计算出工区断层的泥岩涂抹系数；
步骤105：结合区域地质数据，根据工区断层的泥岩涂抹系数，对断层封堵性进行分析。
本发明实施例中，根据预反演的工区相对阻抗数据体，建立工区三维构造模型；根据工区三维构造模型，对单井岩相数据进行属性建模，得到工区相控数据体；在工区相控数据体约束下，对泥质含量数据进行地质统计学建模，得到泥质含量模型；根据泥质含量模型，计算出工区断层的泥岩涂抹系数；结合区域地质数据，根据工区断层的泥岩涂抹系数，对断层封堵性进行分析，从而在工区相控数据体约束下，建模得到泥质含量模型，与现有技术对比，对泥岩分布有了以相控数据体实现的地质约束，可更加真实地以泥质含量模型反映泥岩的实际展布，可准确地对断层封堵性进行定量地分析。
具体实施时，先根据预反演的工区相对阻抗数据体，建立工区三维构造模型。实施例中，预反演的工区相对阻抗数据体，可以按如下方式反演：对加载到反演及建模系统的多种地震数据，利用波阻抗反演，得到工区相对阻抗数据体；所述多种地震数据包括岩相数据，地震解释数据，测井数据及分层数据等其中之一或任意组合。
实施例中，可以先统计整理单井岩相数据，多种地震数据可以包括统计整理后的岩相数据、地震解释数据、测井数据及分层数据的其中之一或任意组合。将多种地震数据加载到反演及建模系统，利用波阻抗反演，获得工区相对阻抗数据体。在一实施例中，工区相对阻抗数据体可为反演波阻抗数据体。
具体实施时，在建立工区三维构造模型后，根据工区三维构造模型，对单井岩相数据进行属性建模，得到工区相控数据体，实施例中可以包括：根据工区三维构造模型，对单井岩相数据进行属性建模中的岩相建模，得到工区相控数据体。
举一例，如图2所示，图2是本发明实施例提供的一种断层封堵性分析方法的实例的流程示意图，包括：
步骤201：基于工区若干单井声波数据、密度数据，以及工区三维地震资料，通过常规宽方位、稀疏脉冲反演等方式可得到工区波阻抗数据；
步骤202：基于工区若干单井测井解释数据、录井数据、岩性描述数据等可得到单井岩相数据，分别通过代码1、2、3分别表示泥岩、砂岩、灰岩等；
步骤203：基于建模系统或软件，以在波阻抗数据为约束下，将得到的岩相数据进行建模，最终得到相控模型，继而得到相控数据体。
具体实施时，在得到工区相控数据体后，在工区相控数据体约束下，对泥质含量数据进行地质统计学建模，得到泥质含量模型，实施例中可以包括：在工区相控数据体约束下，对泥质含量曲线进行地质统计学建模中的自然伽马地质统计学建模，得到泥质含量模型。
泥质含量数据可包括泥质含量Vsh(Shale Volume)数据曲线，和/或自然伽马GR(Natural Gamma Ray)数据曲线。在工区相控数据体约束下，利用泥质含量数据，基于地质统计学概念实现自然伽马地质统计学建模，并最终获得准确的泥质含量模型。与现有技术对比，目前断层封堵性分析方法中，通常直接利用自然伽马曲线进行泥岩含量反演得到泥质含量模型，未考虑沉积地质背景的控制，而本发明实施例中以工区相控数据体作为约束体，建模得到泥质含量模型，很好地实现了相控，对泥岩分布有了地质约束，可更加真实的反映泥岩的实际展布。
具体实施时，根据泥质含量模型，计算出工区断层的泥岩涂抹系数的方法有多种，下面给出其中一个实施例。
在一实施例中，根据泥质含量模型，计算出工区断层的泥岩涂抹系数，如图3所示，包括：
步骤301：根据泥质含量模型，计算断层每一层的泥岩厚度和断移地层厚度；
步骤302：根据断层每一层的泥岩厚度和断移地层厚度，计算断层每一层的泥岩涂抹系数。
实施例中，根据泥质含量模型计算出工区断层的泥岩涂抹系数中，对于任意断层，可截取过目的层的上盘泥质含量模型和下盘泥质含量模型，从而针对目的层的断层封堵性进行评价，举一例，如截取过目的层的以上下盘顶底为边界，上下额外各截取所建工区模型总厚度1/4大小的泥质含量模型。
实施例中，可以按如下公式，计算断层每一层的泥岩厚度：
Pn＝Vsh×ΔZn
其中，Pn表示断层地震剖面第n层的泥岩厚度，Vsh为泥质含量模型，ΔZn表示地震剖面第n层的地层厚度。
具体实施时，所述根据泥质含量模型，计算断层每一层的断移地层厚度的方法有多种，下面给出其中一个实施例。
在一实施例中，根据泥质含量模型，计算断层每一层的断移地层厚度，如图4所示，包括：
步骤401：根据泥质含量模型，获取断层每一层地震剖面的时间厚度，地层速度，以及断层倾角；
步骤402：根据断层每一层地震剖面的时间厚度和地层速度，计算断层每一层地震剖面的垂直断距；
步骤403：根据断层每一层地震剖面的垂直断距和断层倾角，计算断层每一层地震剖面的断移地层厚度。
实施例中，可以按如下公式，计算断层每一层地震剖面的断移地层厚度：
其中，Tn表示断层第n层地震剖面的断移地层厚度，tn表示断层第n层地震剖面的时间厚度，V表示断层地震剖面的地层速度，θ表示断层地震剖面的断层倾角。
实施例中，可以按如下公式，计算断层每一层的泥岩涂抹系数：
其中，SGRn为断层第n层的泥岩涂抹系数，Tn表示断层第n层地震剖面的断移地层厚度，Pn表示断层地震剖面第n层的泥岩厚度，θ表示断层地震剖面的断层倾角，ΔZn表示地震剖面第n层的地层厚度。
实施例中，根据泥质含量模型，可以针对每个小层进行断层泥岩涂抹系数评价，也可以同时评价小层内多条断层泥岩涂抹系数。
具体实施时，在计算出工区断层的泥岩涂抹系数后，结合区域地质数据，根据工区断层的泥岩涂抹系数，对断层封堵性进行分析。
实施例中，结合区域地质数据可以分别评价研究区任意目的层多条断层封堵性，从而对断层封堵性进行分析，实现了对断层封堵性的定量分析。
本发明实施例中，根据预反演的工区相对阻抗数据体，建立工区三维构造模型；根据工区三维构造模型，对单井岩相数据进行属性建模，得到工区相控数据体；在工区相控数据体约束下，对泥质含量数据进行地质统计学建模，得到泥质含量模型；根据泥质含量模型，计算出工区断层的泥岩涂抹系数；结合区域地质数据，根据工区断层的泥岩涂抹系数，对断层封堵性进行分析，从而在工区相控数据体约束下，建模得到泥质含量模型，与现有技术对比，对泥岩分布有了以相控数据体实现的地质约束，可更加真实地以泥质含量模型反映泥岩的实际展布，可准确地对断层封堵性进行定量地分析。
本发明实施例操作简单、实用性高，结果准确，与现有技术相比，利用岩相数据进行相控泥质含量建模，可获得全工区泥质含量模型，使得泥质含量模型更加精确且符合区域地质认识，可为后续断层泥岩涂抹势计算提供可靠的数据，又可简化计算流程。
基于同一发明构思，本发明实施例中还提供了一种断层封堵性分析装置，如下面的实施例所述。由于断层封堵性分析装置解决问题的原理与断层封堵性分析方法相似，因此断层封堵性分析装置的实施可以参见断层封堵性分析方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
图5是本发明实施例的断层封堵性分析装置的一种结构框图，如图5所示，包括：
工区三维构造模型建立模块01，用于根据预反演的工区相对阻抗数据体，建立工区三维构造模型；
相控数据体建模模块02，用于根据工区三维构造模型，对单井岩相数据进行属性建模，得到工区相控数据体；
泥质含量模型建模模块03，用于在工区相控数据体约束下，对泥质含量数据进行地质统计学建模，得到泥质含量模型；
泥岩涂抹系数计算模块04，用于根据泥质含量模型，计算出工区断层的泥岩涂抹系数；
断层封堵性分析模块05，用于结合区域地质数据，根据工区断层的泥岩涂抹系数，对断层封堵性进行分析。
在一实施例中，工区三维构造模型建立模块具体用于：按如下方式反演工区相对阻抗数据体：
对加载到反演及建模系统的多种地震数据，利用波阻抗反演，得到工区相对阻抗数据体；所述多种地震数据包括岩相数据，地震解释数据，测井数据及分层数据的其中之一或任意组合。
在一实施例中，相控数据体建模模块，具体用于：根据工区三维构造模型，对单井岩相数据进行属性建模中的岩相建模，得到工区相控数据体。
在一实施例中，泥质含量模型建模模块，具体用于：在工区相控数据体约束下，对泥质含量曲线进行地质统计学建模中的自然伽马地质统计学建模，得到泥质含量模型。
在一实施例中，泥岩涂抹系数计算模块，具体用于：根据泥质含量模型，计算断层每一层的泥岩厚度和断移地层厚度；
根据断层每一层的泥岩厚度和断移地层厚度，计算断层每一层的泥岩涂抹系数。
在一实施例中，泥岩涂抹系数计算模块，具体用于：按如下公式，计算断层每一层的泥岩厚度：
Pn＝Vsh×ΔZn
其中，Pn表示断层地震剖面第n层的泥岩厚度，Vsh为泥质含量模型，ΔZn表示地震剖面第n层的地层厚度。
在一实施例中，泥岩涂抹系数计算模块，具体用于：根据泥质含量模型，获取断层每一层地震剖面的时间厚度，地层速度，以及断层倾角；根据断层每一层地震剖面的时间厚度和地层速度，计算断层每一层地震剖面的垂直断距；根据断层每一层地震剖面的垂直断距和断层倾角，计算断层每一层地震剖面的断移地层厚度。
在一实施例中，泥岩涂抹系数计算模块，具体用于：按如下公式，计算断层每一层地震剖面的断移地层厚度：
其中，Tn表示断层第n层地震剖面的断移地层厚度，tn表示断层第n层地震剖面的时间厚度，V表示断层地震剖面的地层速度，θ表示断层地震剖面的断层倾角。
在一实施例中，泥岩涂抹系数计算模块，具体用于：按如下公式，计算断层每一层的泥岩涂抹系数：
其中，SGRn为断层第n层的泥岩涂抹系数，Tn表示断层第n层地震剖面的断移地层厚度，Pn表示断层地震剖面第n层的泥岩厚度，θ表示断层地震剖面的断层倾角，ΔZn表示地震剖面第n层的地层厚度。
本发明实施例还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述所述方法。
本发明实施例还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述所述方法的计算机程序。
本发明实施例中，根据预反演的工区相对阻抗数据体，建立工区三维构造模型；根据工区三维构造模型，对单井岩相数据进行属性建模，得到工区相控数据体；在工区相控数据体约束下，对泥质含量数据进行地质统计学建模，得到泥质含量模型；根据泥质含量模型，计算出工区断层的泥岩涂抹系数；结合区域地质数据，根据工区断层的泥岩涂抹系数，对断层封堵性进行分析，从而在工区相控数据体约束下，建模得到泥质含量模型，与现有技术对比，对泥岩分布有了以相控数据体实现的地质约束，可更加真实地以泥质含量模型反映泥岩的实际展布，可准确地对断层封堵性进行定量地分析。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (5)

1.一种断层封堵性分析方法，其特征在于，包括：
根据预反演的工区相对阻抗数据体，建立工区三维构造模型；
根据工区三维构造模型，对单井岩相数据进行属性建模，得到工区相控数据体；
在工区相控数据体约束下，对泥质含量数据进行地质统计学建模，得到泥质含量模型；
根据泥质含量模型，计算出工区断层的泥岩涂抹系数；
结合区域地质数据，根据工区断层的泥岩涂抹系数，对断层封堵性进行分析；
根据工区三维构造模型，对单井岩相数据进行属性建模，得到工区相控数据体，包括：根据工区三维构造模型，对单井岩相数据进行属性建模中的岩相建模，得到工区相控数据体；在工区相控数据体约束下，对泥质含量数据进行地质统计学建模，得到泥质含量模型，包括：在工区相控数据体约束下，对泥质含量曲线进行地质统计学建模中的自然伽马地质统计学建模，得到泥质含量模型；
根据泥质含量模型，计算出工区断层的泥岩涂抹系数，包括：根据泥质含量模型，计算断层每一层的泥岩厚度和断移地层厚度；根据断层每一层的泥岩厚度和断移地层厚度，计算断层每一层的泥岩涂抹系数；
按如下公式，计算断层每一层的泥岩厚度：
Pn＝Vsh×ΔZn
其中，Pn表示断层地震剖面第n层的泥岩厚度，Vsh为泥质含量模型，ΔZn表示地震剖面第n层的地层厚度；
所述根据泥质含量模型，计算断层每一层的断移地层厚度，包括：根据泥质含量模型，获取断层每一层地震剖面的时间厚度，地层速度，以及断层倾角；根据断层每一层地震剖面的时间厚度和地层速度，计算断层每一层地震剖面的垂直断距；根据断层每一层地震剖面的垂直断距和断层倾角，计算断层每一层地震剖面的断移地层厚度；
按如下公式，计算断层每一层地震剖面的断移地层厚度：
其中，Tn表示断层第n层地震剖面的断移地层厚度，tn表示断层第n层地震剖面的时间厚度，V表示断层地震剖面的地层速度，θ表示断层地震剖面的断层倾角；
按如下公式，计算断层每一层的泥岩涂抹系数：
其中，SGRn为断层第n层的泥岩涂抹系数，Tn表示断层第n层地震剖面的断移地层厚度，Pn表示断层地震剖面第n层的泥岩厚度，θ表示断层地震剖面的断层倾角，ΔZn表示地震剖面第n层的地层厚度。
2.如权利要求1所述的方法，其特征在于，所述预反演的工区相对阻抗数据体，按如下方式反演：
对加载到反演及建模系统的多种地震数据，利用波阻抗反演，得到工区相对阻抗数据体；所述多种地震数据包括岩相数据，地震解释数据，测井数据及分层数据的其中之一或任意组合。
3.一种基于权利要求1至2任一所述的断层封堵性分析方法的断层封堵性分析装置，其特征在于，包括：
工区三维构造模型建立模块，用于根据预反演的工区相对阻抗数据体，建立工区三维构造模型；
相控数据体建模模块，用于根据工区三维构造模型，对单井岩相数据进行属性建模，得到工区相控数据体；
泥质含量模型建模模块，用于在工区相控数据体约束下，对泥质含量数据进行地质统计学建模，得到泥质含量模型；
泥岩涂抹系数计算模块，用于根据泥质含量模型，计算出工区断层的泥岩涂抹系数；
断层封堵性分析模块，用于结合区域地质数据，根据工区断层的泥岩涂抹系数，对断层封堵性进行分析；
相控数据体建模模块，具体用于：根据工区三维构造模型，对单井岩相数据进行属性建模中的岩相建模，得到工区相控数据体；泥质含量模型建模模块，具体用于：在工区相控数据体约束下，对泥质含量曲线进行地质统计学建模中的自然伽马地质统计学建模，得到泥质含量模型；
泥岩涂抹系数计算模块，具体用于：根据泥质含量模型，计算断层每一层的泥岩厚度和断移地层厚度；根据断层每一层的泥岩厚度和断移地层厚度，计算断层每一层的泥岩涂抹系数。
4.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至2任一所述方法。
5.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至2任一所述方法的计算机程序。
CN202010331935.3A 2020-04-24 2020-04-24 断层封堵性分析方法及装置 Active CN113552622B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010331935.3A CN113552622B (zh) 2020-04-24 2020-04-24 断层封堵性分析方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010331935.3A CN113552622B (zh) 2020-04-24 2020-04-24 断层封堵性分析方法及装置

Publications (2)

Publication Number Publication Date
CN113552622A CN113552622A (zh) 2021-10-26
CN113552622B true CN113552622B (zh) 2024-05-28

Family

ID=78101247

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010331935.3A Active CN113552622B (zh) 2020-04-24 2020-04-24 断层封堵性分析方法及装置

Country Status (1)

Country Link
CN (1) CN113552622B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104375183A (zh) * 2014-11-20 2015-02-25 中国石油天然气股份有限公司 一种获得断层平面封堵性的方法和装置
CN110632657A (zh) * 2019-08-16 2019-12-31 中国石油天然气股份有限公司 泥岩涂抹型断层封闭性分析方法及装置
CN110632659A (zh) * 2019-08-16 2019-12-31 中国石油天然气股份有限公司 断层空间封闭性分析方法、装置及系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104375183A (zh) * 2014-11-20 2015-02-25 中国石油天然气股份有限公司 一种获得断层平面封堵性的方法和装置
CN110632657A (zh) * 2019-08-16 2019-12-31 中国石油天然气股份有限公司 泥岩涂抹型断层封闭性分析方法及装置
CN110632659A (zh) * 2019-08-16 2019-12-31 中国石油天然气股份有限公司 断层空间封闭性分析方法、装置及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
断层封堵性定量研究在复杂断块油藏中的应用;邢卫新 等;《石油物探》;第45卷(第3期);250-255 *
测井-地震-油藏模拟一体化技术及其在老油田挖潜中的应用;甘利灯 等;《石油物探》;第55卷(第5期);617-639、691 *

Also Published As

Publication number Publication date
CN113552622A (zh) 2021-10-26

Similar Documents

Publication Publication Date Title
US10712472B2 (en) Method and system for forming and using a subsurface model in hydrocarbon operations
WO2019238451A1 (en) A method and a system for modelling and simulating a fractured geological structure
CN110632651B (zh) 断层封闭性定量分析方法及装置
CN105652323B (zh) 一种储层预测方法
CN110632657B (zh) 泥岩涂抹型断层封闭性分析方法及装置
CN105701319B (zh) 一种水平井规则开发井网下的沉积微相建模方法
CN110632652B (zh) 断块油气藏的油气运移方式确定方法及装置
CN110632655B (zh) 非均质型断层的封闭性分析方法及装置
GB2595809A (en) Method, apparatus, and system for quantitative analysis of sealing strength of trap fault
Ganguli Integrated reservoir studies for CO2-enhanced oil recovery and sequestration: Application to an Indian mature oil field
CN105242307A (zh) 复杂碳酸盐储层地震孔隙度获取方法及装置
CN110632656B (zh) 断层的横向封闭性分析方法及装置
CN106125133B (zh) 一种基于气云区约束下的精细速度建模方法
CN106842316B (zh) 裂缝确定方法和装置
CN113552622B (zh) 断层封堵性分析方法及装置
CN111812713A (zh) 断块圈闭确定方法、勘探辅助方法及装置
CN112505754A (zh) 基于高精度层序格架模型的井震协同划分沉积微相的方法
CN114063162B (zh) 一种基于小样本机器学习SVM Vp/Vs预测方法
CN107290803A (zh) 建立基于改进Pride-Lee模型的岩石物理新模型的方法
CN106291748A (zh) 一种获取烃源岩有机碳含量的方法
Gunnarsson 3D modeling in Petrel of geological CO2 storage site
CN113138412B (zh) 深层页岩气孔隙度地震预测方法及装置
CN110568149B (zh) 沉积盆地烃源岩生排烃史精细快速定量模拟方法
US10460051B2 (en) Computationally-efficient modeling of viscous fingering effect for enhanced oil recovery (EOR) agent injected at multiple injection concentrations
CN113534248A (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