CN113466931B - 一种地质时期断层封闭性评价方法及系统 - Google Patents
一种地质时期断层封闭性评价方法及系统 Download PDFInfo
- Publication number
- CN113466931B CN113466931B CN202110599091.5A CN202110599091A CN113466931B CN 113466931 B CN113466931 B CN 113466931B CN 202110599091 A CN202110599091 A CN 202110599091A CN 113466931 B CN113466931 B CN 113466931B
- Authority
- CN
- China
- Prior art keywords
- fault
- point
- section
- activity rate
- geological
- 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
- 238000011156 evaluation Methods 0.000 title claims abstract description 51
- 230000000694 effects Effects 0.000 claims abstract description 84
- 238000000034 method Methods 0.000 claims abstract description 41
- 238000012937 correction Methods 0.000 claims abstract description 34
- 239000011435 rock Substances 0.000 claims abstract description 34
- 238000004364 calculation method Methods 0.000 claims abstract description 31
- 239000003208 petroleum Substances 0.000 claims abstract description 28
- 238000007789 sealing Methods 0.000 claims description 35
- 238000004590 computer program Methods 0.000 claims description 17
- 238000005553 drilling Methods 0.000 claims description 8
- 238000003860 storage Methods 0.000 claims description 7
- 238000001514 detection method Methods 0.000 claims description 6
- 230000000903 blocking effect Effects 0.000 claims description 5
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 239000004927 clay Substances 0.000 claims description 3
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 10
- 238000011161 development Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 8
- 229930195733 hydrocarbon Natural products 0.000 description 6
- 150000002430 hydrocarbons Chemical class 0.000 description 6
- 230000000875 corresponding effect Effects 0.000 description 5
- 238000013508 migration Methods 0.000 description 5
- 230000005012 migration Effects 0.000 description 5
- 239000004215 Carbon black (E152) Substances 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000009826 distribution Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000001419 dependent effect Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- BPQQTUXANYXVAA-UHFFFAOYSA-N Orthosilicate Chemical compound [O-][Si]([O-])([O-])[O-] BPQQTUXANYXVAA-UHFFFAOYSA-N 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 239000000945 filler Substances 0.000 description 1
- 210000001503 joint Anatomy 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
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/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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
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
本发明提出了一种地质时期断层封闭性评价方法及系统，其中，该方法包括：采集石油地震勘探数据；根据石油地震勘探数据，建立地下断层模型，并计算出断层的断距、断面上每个点的断层倾角值和垂直断距；根据地下断层模型确定待评价的目的层段在断面上的空间投影窗口位置，并选取断面上要进行封闭性评价的计算点；计算断面上每个点的断层倾角校正因子；计算计算点的有效断层岩的泥质含量ESGR值；根据断层活动速率计算得到在一地质时期的断面上每个点的归一化活动速率；计算代表地质时期的断层封闭能力的参数Pf；统计勘探区域上已发现油藏断层封闭下限阈值，将Pf值与封闭下限阈值进行比较，根据比较结果定量评价断层在一地质时期的封闭性能。
Description
技术领域
本发明涉及油气田勘探开发技术领域，尤指一种地质时期断层封闭性评价方法及系统。
背景技术
断层裂缝之间充填有厚度、岩性等分布不均一的断层岩，断层的侧向封闭能力不是取决于两盘岩性对接情况，而是取决于断层岩与目的盘储层之间的排替压力差异，而对于碎屑岩地层而言，断层岩的排替压力大小主要取决于断层岩中泥质含量的大小，断层岩泥质含量越大，断层侧向封闭能力越强；反之越弱。研究者们试图通过研究断层岩泥质含量的办法间接或直接评价断层的侧向封闭油气的能力。
断层岩中的泥质主要有泥页岩层的泥岩涂抹和砂泥混杂充填物。
目前，被国内外学者比较公认的评价泥岩涂抹发育程度的方法主要有两个，一个是泥岩涂抹因子SSF，另外一个是泥岩涂抹势CSP，但上述方法的不足是仅考虑了断距和泥岩层厚度，没有考虑不同位置泥岩涂抹量不同，也没有考虑到断层不同位置的倾角是变化的这一客观现象。对于断层中砂泥混杂充填物的研究，主要采用有效断层泥比率ESGR方法，即为断层在断移过程中滑过某点泥岩层累加厚度与断距的比值，该方法考虑了实际的断裂作用过程中，砂岩中的泥(硅酸盐)对断裂带中断层泥的比率的贡献，ESGR更能准确预测断裂带中断层泥含量和泥岩涂抹的程度。
国内外研究方法集中关注了现今地质条件下断层的封闭性，是对当前状态的静态评价，没有考虑断层的历史活动性，无法确定断层的古封闭性。对于断块圈闭的评价而言，断层的古封闭性是重大地质勘探风险。当前研究现状认为断层侧向封闭能力的最根本因素是断裂带泥质成分含量，用参数SGR表征，该参数是垂直断距、断裂两侧岩层厚度、地层泥质含量三个参数的函数。上述参数的弊端在于只能评价当前状态的封闭性，没有考虑断层的历史活动性和演化过程，无法确定断层的古封闭性，而且也没有考虑高陡断层和低角度平缓断层在封闭性上的差异。
综上来看，亟需一种可以克服上述缺陷，能够准确评价多种时期断层封闭性的技术方案。
发明内容
为了解决背景技术中所提到的技术问题，本发明提出了一种地质时期断层封闭性评价方法及系统，本发明适合碎屑岩地层的断层封闭性定量判定；通过建立考虑断层演化中动态活动因素、断层倾角因素的断层封闭性定量评价地质与数学模型，计算断面上某地质时期每个点的归一化活动速率Vf、求取计算点的有效断层岩的泥质含量ESGR值、计算断面上每个点的断层倾角校正因子；并利用上述3个参数所代表的地质因素求取古封闭性评价参数Pf值，再统计区域上已发现油藏断层封闭下限阈值，将计算的待预测断层的地质时期封闭性Pf值与封闭下限阈值比较，定量评价断层某地质时期封闭与否及封闭能力的大小，提高断层封闭性评价的准确度，以期降低断层圈闭钻探风险，也为开发阶段优化断块油气藏开发方案提供地质评价参考。
在本发明实施例的第一方面，提出了一种地质时期断层封闭性评价方法，该方法包括：
采集石油地震勘探数据；
根据所述石油地震勘探数据，建立地下断层模型，并计算出断层的断距、断面上每个点的断层倾角值和垂直断距；
根据地下断层模型确定待评价的目的层段在断面上的空间投影窗口位置，并选取断面上要进行封闭性评价的计算点；
根据断面上每个点的断层倾角值，计算得到断面上每个点的断层倾角校正因子；
根据石油地震勘探数据及垂直断距，计算得到所述计算点的有效断层岩的泥质含量ESGR值；
根据断层的断距得到断层活动速率，并根据所述断层活动速率计算得到在一地质时期的断面上每个点的归一化活动速率；
根据断面上每个点的断层倾角校正因子、计算点的有效断层岩的泥质含量ESGR值及在一地质时期的断面上每个点的归一化活动速率，计算得到代表地质时期的断层封闭能力的参数Pf；
统计勘探区域上已发现油藏断层封闭下限阈值，将待评价断层的地质时期封闭性Pf值与封闭下限阈值比较，根据比较结果定量评价断层在一地质时期的封闭性能。
进一步的，根据地下断层模型确定待评价的目的层段在断面上的空间投影窗口位置，并选取断面上要进行封闭性评价的计算点，包括：
根据石油地震勘探和钻井数据建立地层模型；
将所述地层模型与地下断层模型耦合，确定待评价的目的层段在断面上的空间投影窗口位置，并选取断面上要进行封闭性评价的计算点。
进一步的，根据断面上每个点的断层倾角值，计算得到断面上每个点的断层倾角校正因子，包括：
断面上每个点的断层倾角校正因子计算式为：
其中，f(θ)为断层倾角校正因子，无量纲；
θ为断层倾角值，°。
进一步的，根据石油地震勘探数据及垂直断距，计算得到所述计算点的有效断层岩的泥质含量ESGR值，包括：
其中，ESGR为有效断层岩的泥质含量，％；
i为滑过断点的第i层岩层；
Ti为滑过断点的第i层岩层厚度，m；
Vshi为滑过断点的第i层岩层泥质含量，％；
D为断层的垂直断距，m。
进一步的，根据断层的断距得到断层活动速率，并根据所述断层活动速率计算得到在一地质时期的断面上每个点的归一化活动速率，包括：
根据一地质时期断层的断距，计算断层活动后的累加断距与对应活动时期的时间长度的比值，所述比值为一地质时期的断层活动速率；
归一化活动速率的计算式为：
其中，Vf为归一化活动速率，无量纲，％；
V为一地质时期的断层活动速率，m/Ma；
Vmax和Vmin分别为待评价断层在历史活动时期所具备的最大活动速率和最小的活动速率，m/Ma。
进一步的，根据断面上每个点的断层倾角校正因子、计算点的有效断层岩的泥质含量ESGR值及在一地质时期的断面上每个点的归一化活动速率，计算得到代表地质时期的断层封闭能力的参数Pf，包括：
代表地质时期的断层封闭能力的参数Pf的计算式为：
其中，Pf为代表地质时期的断层封闭能力的参数。
进一步的，统计勘探区域上已发现油藏断层封闭下限阈值，将待评价断层的地质时期封闭性Pf值与封闭下限阈值比较，根据比较结果定量评价断层在一地质时期的封闭性能，包括：
将Pf与封闭下限阈值比较；其中，若Pf比封闭下限阈值大，断层侧向封闭，反之断层侧向开启；
封闭能力通过Pf与封闭下限阈值的差值表征，若差值越大，断层封闭能力越强，反之越弱。
在本发明实施例的第二方面，提出了一种地质时期断层封闭性评价系统，该系统包括：
数据采集模块，用于采集石油地震勘探数据；
模型建立模块，用于根据所述石油地震勘探数据，建立地下断层模型，并计算出断层的断距、断面上每个点的断层倾角值和垂直断距；
计算点选取模块，用于根据地下断层模型确定待评价的目的层段在断面上的空间投影窗口位置，并选取断面上要进行封闭性评价的计算点；
校正因子计算模块，用于根据断面上每个点的断层倾角值，计算得到断面上每个点的断层倾角校正因子；
泥质含量计算模块，用于根据石油地震勘探数据及垂直断距，计算得到所述计算点的有效断层岩的泥质含量ESGR值；
活动速率计算模块，用于根据断层的断距得到断层活动速率，并根据所述断层活动速率计算得到在一地质时期的断面上每个点的归一化活动速率；
封闭能力计算模块，用于根据断面上每个点的断层倾角校正因子、计算点的有效断层岩的泥质含量ESGR值及在一地质时期的断面上每个点的归一化活动速率，计算得到代表地质时期的断层封闭能力的参数Pf；
封闭性评价模块，用于统计勘探区域上已发现油藏断层封闭下限阈值，将待评价断层的地质时期封闭性Pf值与封闭下限阈值比较，根据比较结果定量评价断层在一地质时期的封闭性能。
在本发明实施例的第三方面，提出了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现地质时期断层封闭性评价方法。
在本发明实施例的第四方面，提出了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现地质时期断层封闭性评价方法。
本发明提出的地质时期断层封闭性评价方法及系统引入断层活动速率这一与断层历史演化相关联的因素，反映了某地质时期的断层封闭特征；而且考虑了断层倾角变化的因素给断层岩泥质含量计算带来的影响，本发明所带来的评价效果更符合实际地质规律，有效提高了断层封闭性评价的准确度，降低断层圈闭钻探风险，为开发阶段优化断块油气藏开发方案提供地质评价参考。
附图说明
为了更清楚地说明本申请实施例技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其它的附图。
图1是本发明一实施例的地质时期断层封闭性评价方法流程示意图。
图2是本发明一具体实施例的断层倾角变化对断层倾角校正因子的影响关系示意图。
图3A至图3C是本发明一具体实施例的不同地质时期的封闭能力参数Pf计算结果图。
图4是本发明一实施例的地质时期断层封闭性评价系统架构示意图。
图5是本发明一实施例的计算机设备结构示意图。
具体实施方式
下面将参考若干示例性实施方式来描述本发明的原理和精神。应当理解，给出这些实施方式仅仅是为了使本领域技术人员能够更好地理解进而实现本发明，而并非以任何方式限制本发明的范围。相反，提供这些实施方式是为了使本公开更加透彻和完整，并且能够将本公开的范围完整地传达给本领域的技术人员。
本领域技术人员知道，本发明的实施方式可以实现为一种系统、装置、设备、方法或计算机程序产品。因此，本公开可以具体实现为以下形式，即：完全的硬件、完全的软件(包括固件、驻留软件、微代码等)，或者硬件和软件结合的形式。
根据本发明的实施方式，提出了一种地质时期断层封闭性评价方法及系统。
下面参考本发明的若干代表性实施方式，详细阐释本发明的原理和精神。
图1是本发明一实施例的地质时期断层封闭性评价方法流程示意图。如图1所示，该方法包括：
步骤S1，采集石油地震勘探数据；
步骤S2，根据所述石油地震勘探数据，建立地下断层模型，并计算出断层的断距、断面上每个点的断层倾角值和垂直断距；
步骤S3，根据地下断层模型确定待评价的目的层段在断面上的空间投影窗口位置，并选取断面上要进行封闭性评价的计算点；
步骤S4，根据断面上每个点的断层倾角值，计算得到断面上每个点的断层倾角校正因子；
步骤S5，根据石油地震勘探数据及垂直断距，计算得到所述计算点的有效断层岩的泥质含量ESGR值；
步骤S6，根据断层的断距得到断层活动速率，并根据所述断层活动速率计算得到在一地质时期的断面上每个点的归一化活动速率；
步骤S7，根据断面上每个点的断层倾角校正因子、计算点的有效断层岩的泥质含量ESGR值及在一地质时期的断面上每个点的归一化活动速率，计算得到代表地质时期的断层封闭能力的参数Pf；
步骤S8，统计勘探区域上已发现油藏断层封闭下限阈值，将待评价断层的地质时期封闭性Pf值与封闭下限阈值比较，根据比较结果定量评价断层在一地质时期的封闭性能。
为了对上述地质时期断层封闭性评价方法进行更为清楚的解释，下面结合每一步骤进行详细说明。
步骤S1：
采集石油地震勘探数据。
步骤S2：
根据所述石油地震勘探数据，建立地下断层模型，并计算出断层的断距、断面上每个点的断层倾角值和垂直断距；
步骤S3：
根据石油地震勘探和钻井数据建立地层模型；
将所述地层模型与地下断层模型耦合，确定待评价的目的层段在断面上的空间投影窗口位置，并选取断面上要进行封闭性评价的计算点。
步骤S4：
根据断面上每个点的断层倾角值，计算得到断面上每个点的断层倾角校正因子。
参考图2，为本发明一具体实施例的断层倾角变化对断层倾角校正因子的影响关系示意图。
断面上每个点的断层倾角校正因子计算式为：
其中，f(θ)为断层倾角校正因子，无量纲；
θ为断层倾角值，°。
步骤S5：
根据石油地震勘探数据及垂直断距，计算得到所述计算点的有效断层岩的泥质含量ESGR值，包括：
其中，ESGR为有效断层岩的泥质含量，％；
i为滑过断点的第i层岩层；
Ti为滑过断点的第i层岩层厚度，m；
Vshi为滑过断点的第i层岩层泥质含量，％；
D为断层的垂直断距，m。
步骤S6：
根据断层的断距得到断层活动速率，并根据所述断层活动速率计算得到在一地质时期的断面上每个点的归一化活动速率，包括：
根据一地质时期断层的断距，计算断层活动后的累加断距与对应活动时期的时间长度的比值，所述比值为一地质时期的断层活动速率；
归一化活动速率的计算式为：
其中，Vf为归一化活动速率，无量纲，％；
V为一地质时期的断层活动速率，m/Ma；
Vmax和Vmin分别为待评价断层在历史活动时期所具备的最大活动速率和最小的活动速率，m/Ma。
步骤S7：
根据断面上每个点的断层倾角校正因子、计算点的有效断层岩的泥质含量ESGR值及在一地质时期的断面上每个点的归一化活动速率，计算得到代表地质时期的断层封闭能力的参数Pf，包括：
代表地质时期的断层封闭能力的参数Pf的计算式为：
其中，Pf为代表地质时期的断层封闭能力的参数。
步骤S8：
统计勘探区域上已发现油藏断层封闭下限阈值，将待评价断层的地质时期封闭性Pf值与封闭下限阈值比较，根据比较结果定量评价断层在一地质时期的封闭性能，包括：
将Pf与封闭下限阈值比较；其中，若Pf比封闭下限阈值大，断层侧向封闭，反之断层侧向开启；
封闭能力通过Pf与封闭下限阈值的差值表征，若差值越大，断层封闭能力越强，反之越弱。
需要说明的是，尽管在上述实施例及附图中以特定顺序描述了本发明方法的操作，但是，这并非要求或者暗示必须按照该特定顺序来执行这些操作，或是必须执行全部所示的操作才能实现期望的结果。附加地或备选地，可以省略某些步骤，将多个步骤合并为一个步骤执行，和/或将一个步骤分解为多个步骤执行。
为了对上述地质时期断层封闭性评价方法进行更为清楚的解释，下面结合一个具体的实施例来进行说明，然而值得注意的是该实施例仅是为了更好地说明本发明，并不构成对本发明不当的限定。
以苏丹Muglad盆地Kaikang凹陷西斜坡某断层为例，采用本发明提出的地质时期断层封闭性评价方法进行断层封闭性评价。
参考图3A，示例性展示了断层的归一化活动速率；图3B示例性展示了白垩纪末期油气关键成藏时刻断层的Pf分布情况；图3C示例性展示了古近纪末期时断层的PfPf的分布情况。
该断层的倾角在40-60度之间，对应的断层倾角校正因子f(θ)位于1.12-1.08；
另一参数SGR在下白垩统、上白垩统和古近系目的层段分别为24％～29％、28％～31％、34％～38％。
计算目的层段Vf所需的下白垩统、上白垩统和古近系所对应的地质历史时窗分别取值为17、20.6、44.8Ma，断距是根据构造图计算而来的。
该断层是北西-南东方向延展的，自北而南长约120km，在横向上具有分段生长的特征，见图3A，在平面上可分为三段：其中在第二分段内发现上白垩统的构造油藏；在第三分段中上白垩统圈闭钻探失利，在古近系中发现油气显示；第一分段尚未有钻井。
该断层地质时期历经3个演化时期：早白垩世、晚白垩世和古近纪活动时期。图3A表明整体上该断层在早白垩世末期和晚白垩世末期的活动速率普遍高于古近纪末期，这说明在该区的关键运移成藏期(晚白垩世末期)，该断层是活动较为强烈的，相应的封闭性能弱一些，有利于下白垩统源岩生成的油气沿着断面垂向运移到上白垩统的主要目的层；但是在横向分段的不同部位，其封闭能力有差异。
图3B中表明在关键油气成藏时期(晚白垩世末期)第二分段封闭性参数Pf整体小于第三分段，表明第二分段封闭性明显弱于第三分段，这说明该段是利于油气运移的，这与实际发现的油气聚集在第二分段的事实是符合的，相反，在油气运移的关键时期，第三分段封闭性高于第二分段，不利于当时油气的运移。
在图3C上整体表现为较强的封闭性，在第二分段Pf整体表现为高值，说明在古近纪该段断层具备良好了封闭性，这有利于油气的保存；但在第三分段Pf值呈现低值(根据第二分段油气分布情况，统计此处取Pf的阈值为5.2)，说明第三分段封闭性差，造成该段封闭性差的原因在于该段在古近纪的活动速率变大(见图3A)；由此可能会造成早期形成的油气藏沿着断层封闭性差的部位向上调整，最终致使油气散逸地表或者分散在古近系的圈闭中，这与在第三分段针对上白垩统断块圈闭目标钻探失利、在古近系中发现油气显示的事实是相吻合的，说明本发明的效果符合地质规律，具有可靠性。
在介绍了本发明示例性实施方式的方法之后，接下来，参考图4对本发明示例性实施方式的地质时期断层封闭性评价系统进行介绍。
地质时期断层封闭性评价系统的实施可以参见上述方法的实施，重复之处不再赘述。以下所使用的术语“模块”或者“单元”，可以是实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
基于同一发明构思，本发明还提出了一种地质时期断层封闭性评价系统，如图4所示，该系统包括：
数据采集模块410，用于采集石油地震勘探数据；
模型建立模块420，用于根据所述石油地震勘探数据，建立地下断层模型，并计算出断层的断距、断面上每个点的断层倾角值和垂直断距；
计算点选取模块430，用于根据地下断层模型确定待评价的目的层段在断面上的空间投影窗口位置，并选取断面上要进行封闭性评价的计算点；
校正因子计算模块440，用于根据断面上每个点的断层倾角值，计算得到断面上每个点的断层倾角校正因子；
泥质含量计算模块450，用于根据石油地震勘探数据及垂直断距，计算得到所述计算点的有效断层岩的泥质含量ESGR值；
活动速率计算模块460，用于根据断层的断距得到断层活动速率，并根据所述断层活动速率计算得到在一地质时期的断面上每个点的归一化活动速率；
封闭能力计算模块470，用于根据断面上每个点的断层倾角校正因子、计算点的有效断层岩的泥质含量ESGR值及在一地质时期的断面上每个点的归一化活动速率，计算得到代表地质时期的断层封闭能力的参数Pf；
封闭性评价模块480，用于统计勘探区域上已发现油藏断层封闭下限阈值，将待评价断层的地质时期封闭性Pf值与封闭下限阈值比较，根据比较结果定量评价断层在一地质时期的封闭性能。
应当注意，尽管在上文详细描述中提及了地质时期断层封闭性评价系统的若干模块，但是这种划分仅仅是示例性的并非强制性的。实际上，根据本发明的实施方式，上文描述的两个或更多模块的特征和功能可以在一个模块中具体化。反之，上文描述的一个模块的特征和功能可以进一步划分为由多个模块来具体化。
基于前述发明构思，如图5所示，本发明还提出了一种计算机设备500，包括存储器510、处理器520及存储在存储器510上并可在处理器520上运行的计算机程序530，所述处理器520执行所述计算机程序530时实现前述地质时期断层封闭性评价方法。
基于前述发明构思，本发明提出了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现前述地质时期断层封闭性评价方法。
本发明提出的地质时期断层封闭性评价方法及系统引入断层活动速率这一与断层历史演化相关联的因素，反映了某地质时期的断层封闭特征；而且考虑了断层倾角变化的因素给断层岩泥质含量计算带来的影响，本发明所带来的评价效果更符合实际地质规律，有效提高了断层封闭性评价的准确度，降低断层圈闭钻探风险，为开发阶段优化断块油气藏开发方案提供地质评价参考。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应说明的是：以上所述实施例，仅为本发明的具体实施方式，用以说明本发明的技术方案，而非对其限制，本发明的保护范围并不局限于此，尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。
Claims (7)
1.一种地质时期断层封闭性评价方法，其特征在于，该方法包括：
采集石油地震勘探数据；
根据所述石油地震勘探数据，建立地下断层模型，并计算出断层的断距、断面上每个点的断层倾角值和垂直断距；
根据地下断层模型确定待评价的目的层段在断面上的空间投影窗口位置，并选取断面上要进行封闭性评价的计算点；
根据断面上每个点的断层倾角值，计算得到断面上每个点的断层倾角校正因子；其中，包括：
断面上每个点的断层倾角校正因子计算式为：
其中，f(θ)为断层倾角校正因子，无量纲；
θ为断层倾角值，°；
根据石油地震勘探数据及垂直断距，计算得到所述计算点的有效断层岩的泥质含量ESGR值；
根据断层的断距得到断层活动速率，并根据所述断层活动速率计算得到在某一地质时期的断面上每个点的归一化活动速率Vf；其中，包括：
根据某一地质时期断层的断距，计算断层活动后的累加断距与对应活动时期的时间长度的比值，所述比值为某一地质时期的断层活动速率V；
根据断层活动速率V得到断层归一化活动速率Vf；归一化活动速率Vf的计算式为：
其中，Vf为归一化活动速率，无量纲，％；
V为某一地质时期的断层活动速率，m/Ma；
Vmax和Vmin分别为待评价断层在历史活动时期所具备的最大活动速率和最小的活动速率，m/Ma；
根据断面上每个点的断层倾角校正因子、计算点的有效断层岩的泥质含量ESGR值及在某一地质时期的断面上每个点的归一化活动速率，计算得到代表地质时期的断层封闭能力的参数Pf；其中，包括：
代表地质时期的断层封闭能力的参数Pf的计算式为：
其中，Pf为代表地质时期的断层封闭能力的参数；
f(θ)为断层倾角校正因子，无量纲；
ESGR为有效断层岩的泥质含量，％；
Vf为归一化活动速率，无量纲，％；
统计勘探区域上已发现油藏断层封闭下限阈值，将待评价断层的地质时期封闭性Pf值与封闭下限阈值比较，根据比较结果定量评价断层在某一地质时期的封闭性能。
2.根据权利要求1所述的地质时期断层封闭性评价方法，其特征在于，根据地下断层模型确定待评价的目的层段在断面上的空间投影窗口位置，并选取断面上要进行封闭性评价的计算点，包括：
根据石油地震勘探和钻井数据建立地层模型；
将所述地层模型与地下断层模型耦合，确定待评价的目的层段在断面上的空间投影窗口位置，并选取断面上要进行封闭性评价的计算点。
3.根据权利要求1所述的地质时期断层封闭性评价方法，其特征在于，根据石油地震勘探数据及垂直断距，计算得到所述计算点的有效断层岩的泥质含量ESGR值，包括：
其中，ESGR为有效断层岩的泥质含量，％；
i为滑过断点的第i层岩层；
Ti为滑过断点的第i层岩层厚度，m；
Vshi为滑过断点的第i层岩层泥质含量，％；
D为断层的垂直断距，m。
4.根据权利要求3所述的地质时期断层封闭性评价方法，其特征在于，统计勘探区域上已发现油藏断层封闭下限阈值，将待评价断层的地质时期封闭性Pf值与封闭下限阈值比较，根据比较结果定量评价断层在某一地质时期的封闭性能，包括：
将Pf与封闭下限阈值比较；其中，若Pf比封闭下限阈值大，断层侧向封闭，反之断层侧向开启；
封闭能力通过Pf与封闭下限阈值的差值表征，若差值越大，断层封闭能力越强，反之越弱。
5.一种地质时期断层封闭性评价系统，其特征在于，该系统包括：
数据采集模块，用于采集石油地震勘探数据；
模型建立模块，用于根据所述石油地震勘探数据，建立地下断层模型，并计算出断层的断距、断面上每个点的断层倾角值和垂直断距；
计算点选取模块，用于根据地下断层模型确定待评价的目的层段在断面上的空间投影窗口位置，并选取断面上要进行封闭性评价的计算点；
校正因子计算模块，用于根据断面上每个点的断层倾角值，计算得到断面上每个点的断层倾角校正因子；其中，包括：
断面上每个点的断层倾角校正因子计算式为：
其中，f(θ)为断层倾角校正因子，无量纲；
θ为断层倾角值，°；
泥质含量计算模块，用于根据石油地震勘探数据及垂直断距，计算得到所述计算点的有效断层岩的泥质含量ESGR值；
活动速率计算模块，用于根据断层的断距得到断层活动速率，并根据所述断层活动速率计算得到在某一地质时期的断面上每个点的归一化活动速率Vf；其中，包括：
根据某一地质时期断层的断距，计算断层活动后的累加断距与对应活动时期的时间长度的比值，所述比值为某一地质时期的断层活动速率V；
根据断层活动速率V得到断层归一化活动速率Vf；归一化活动速率Vf的计算式为：
其中，Vf为归一化活动速率，无量纲，％；
V为某一地质时期的断层活动速率，m/Ma；
Vmax和Vmin分别为待评价断层在历史活动时期所具备的最大活动速率和最小的活动速率，m/Ma；
封闭能力计算模块，用于根据断面上每个点的断层倾角校正因子、计算点的有效断层岩的泥质含量ESGR值及在某一地质时期的断面上每个点的归一化活动速率，计算得到代表地质时期的断层封闭能力的参数Pf；其中，包括：
代表地质时期的断层封闭能力的参数Pf的计算式为：
其中，Pf为代表地质时期的断层封闭能力的参数；
f(θ)为断层倾角校正因子，无量纲；
ESGR为有效断层岩的泥质含量，％；
Vf为归一化活动速率，无量纲，％；
封闭性评价模块，用于统计勘探区域上已发现油藏断层封闭下限阈值，将待评价断层的地质时期封闭性Pf值与封闭下限阈值比较，根据比较结果定量评价断层在某一地质时期的封闭性能。
6.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至4任一所述方法。
7.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现权利要求1至4任一所述方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110599091.5A CN113466931B (zh) | 2021-05-31 | 2021-05-31 | 一种地质时期断层封闭性评价方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110599091.5A CN113466931B (zh) | 2021-05-31 | 2021-05-31 | 一种地质时期断层封闭性评价方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113466931A CN113466931A (zh) | 2021-10-01 |
| CN113466931B true CN113466931B (zh) | 2023-08-22 |
Family
ID=77871945
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110599091.5A Active CN113466931B (zh) | 2021-05-31 | 2021-05-31 | 一种地质时期断层封闭性评价方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113466931B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118210044A (zh) * | 2024-02-29 | 2024-06-18 | 重庆地质矿产研究院 | 一种利用本底参数的断层活动风险评价方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106324676A (zh) * | 2016-10-21 | 2017-01-11 | 中国石油天然气股份有限公司 | 一种确定断层封闭性的方法及装置 |
| CN108241181A (zh) * | 2018-01-30 | 2018-07-03 | 中国石油化工股份有限公司 | 一种断层封闭性的评价方法 |
| CN111337974A (zh) * | 2020-02-17 | 2020-06-26 | 中国石油天然气股份有限公司 | 断层封闭性确定方法及装置 |
| CN111638550A (zh) * | 2019-03-01 | 2020-09-08 | 中国石油化工股份有限公司 | 一种断裂控藏的评价方法 |
| WO2021008630A1 (zh) * | 2019-07-18 | 2021-01-21 | 中国石油天然气股份有限公司 | 圈闭断层的封闭性定量分析方法、装置及系统 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20190227185A1 (en) * | 2017-05-26 | 2019-07-25 | Chevron U.S.A. Inc. | System and method for predicting fault seal from seismic data |
-
2021
- 2021-05-31 CN CN202110599091.5A patent/CN113466931B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106324676A (zh) * | 2016-10-21 | 2017-01-11 | 中国石油天然气股份有限公司 | 一种确定断层封闭性的方法及装置 |
| CN108241181A (zh) * | 2018-01-30 | 2018-07-03 | 中国石油化工股份有限公司 | 一种断层封闭性的评价方法 |
| CN111638550A (zh) * | 2019-03-01 | 2020-09-08 | 中国石油化工股份有限公司 | 一种断裂控藏的评价方法 |
| WO2021008630A1 (zh) * | 2019-07-18 | 2021-01-21 | 中国石油天然气股份有限公司 | 圈闭断层的封闭性定量分析方法、装置及系统 |
| CN111337974A (zh) * | 2020-02-17 | 2020-06-26 | 中国石油天然气股份有限公司 | 断层封闭性确定方法及装置 |
Non-Patent Citations (1)
| Title |
|---|
| 冀中坳陷大柳泉地区油源断裂优势通道输导油气能力定量评价;王浩然;《中国博士学位论文全文数据库 基础科学辑》(第01期);第1-115页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113466931A (zh) | 2021-10-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7363164B2 (en) | Method of evaluating fluid saturation characteristics in a geological formation | |
| Jacobi et al. | Integrated petrophysical evaluation of shale gas reservoirs | |
| CN110632651B (zh) | 断层封闭性定量分析方法及装置 | |
| WO2021008630A1 (zh) | 圈闭断层的封闭性定量分析方法、装置及系统 | |
| CN104678438B (zh) | 一种co2地质封存中四维地震资料co2分布预测的方法 | |
| RU2505676C2 (ru) | Способ определения коэффициента обводненности и состава притока нефтяной скважины | |
| CN109239778B (zh) | 一种断层侧向封闭性定量评价方法 | |
| CN110632655B (zh) | 非均质型断层的封闭性分析方法及装置 | |
| CN110632652B (zh) | 断块油气藏的油气运移方式确定方法及装置 | |
| US9739902B2 (en) | System and method for remaining resource mapping | |
| CN108508484A (zh) | 一种适用于张扭性断层泥岩涂抹效果的评价方法 | |
| CN111337974B (zh) | 断层封闭性确定方法及装置 | |
| CN113466931B (zh) | 一种地质时期断层封闭性评价方法及系统 | |
| CN110632656B (zh) | 断层的横向封闭性分析方法及装置 | |
| CN106772584B (zh) | 一种顾及地质构造条件约束的精细体元属性赋值方法 | |
| CN110633557B (zh) | 一种煤层气构造有利区识别方法 | |
| CN115542419A (zh) | 海洋天然气水合物钻后综合评价方法、系统和介质 | |
| CN110244357B (zh) | 一种构造破坏型古油藏的判识方法 | |
| CN113534248A (zh) | 圈闭断层的封闭性定量分析方法、装置及系统 | |
| Faraji et al. | Estimation and mapping of vitrinite reflectance from seismic data in South Pars gas field, Persian Gulf | |
| CN113031056B (zh) | 构造约束下断层封闭性分析方法及装置 | |
| CN111650644A (zh) | 盐胶结砂岩定量预测方法 | |
| CN113687412B (zh) | 盐间地层压力预测方法、装置、电子设备及介质 | |
| Haris et al. | Reservoir Compartment Assessment: A Case Study Of Bangko And Bekasap Formation, Central Sumatra Basin Indonesia | |
| CN113495293B (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 |