CN113051697A - 古气藏地层孔隙压力演化史曲线生成方法及装置 - Google Patents
古气藏地层孔隙压力演化史曲线生成方法及装置 Download PDFInfo
- Publication number
- CN113051697A CN113051697A CN201911374034.6A CN201911374034A CN113051697A CN 113051697 A CN113051697 A CN 113051697A CN 201911374034 A CN201911374034 A CN 201911374034A CN 113051697 A CN113051697 A CN 113051697A
- Authority
- CN
- China
- Prior art keywords
- inclusion
- pore pressure
- gas reservoir
- stratum
- history curve
- 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.)
- Granted
Links
- 239000011148 porous material Substances 0.000 title claims abstract description 127
- 238000000034 method Methods 0.000 title claims abstract description 71
- 239000012530 fluid Substances 0.000 claims abstract description 41
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 37
- 239000007788 liquid Substances 0.000 claims abstract description 33
- 239000000203 mixture Substances 0.000 claims abstract description 33
- 238000007710 freezing Methods 0.000 claims abstract description 30
- 230000008014 freezing Effects 0.000 claims abstract description 29
- 238000002310 reflectometry Methods 0.000 claims description 20
- 239000004079 vitrinite Substances 0.000 claims description 20
- 230000015654 memory Effects 0.000 claims description 19
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 16
- 238000004088 simulation Methods 0.000 claims description 13
- 238000004590 computer program Methods 0.000 claims description 12
- 238000012937 correction Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 18
- 238000013461 design Methods 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 7
- 239000004215 Carbon black (E152) Substances 0.000 description 6
- 238000009825 accumulation Methods 0.000 description 6
- 229930195733 hydrocarbon Natural products 0.000 description 6
- 150000002430 hydrocarbons Chemical class 0.000 description 6
- 239000012071 phase Substances 0.000 description 6
- 239000002734 clay mineral Substances 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 239000012267 brine Substances 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 239000007791 liquid phase Substances 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- HPALAKNZSZLMCH-UHFFFAOYSA-M sodium;chloride;hydrate Chemical compound O.[Na+].[Cl-] HPALAKNZSZLMCH-UHFFFAOYSA-M 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000002427 irreversible effect Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 239000011259 mixed solution Substances 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Agronomy & Crop Science (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Animal Husbandry (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Measuring Fluid Pressure (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种古气藏地层孔隙压力演化史曲线生成方法及装置，该方法包括：获取选定的包裹体的气液比、均一温度和冰点温度；根据所述冰点温度确定所述包裹体的流体盐度；确定所述包裹体的真实成分；针对所述包裹体的真实成分确定所述包裹体在所述均一温度时的最小捕获压力，并根据该最小捕获压力确定所述包裹体对应的古气藏成藏时的古孔隙压力；根据所述古孔隙压力以及当前测量的所述包裹体所在地层的孔隙压力确定所述包裹体对应的古气藏的孔隙压力演化史曲线。本发明提供了一种设计合理、定量化计算古气藏地层孔隙压力演化过程的方法。
Description
技术领域
本发明涉及油气成藏动力学领域，具体而言，涉及一种古气藏地层孔隙压力演化史曲线生成方法及装置。
背景技术
成藏动力学是综合利用地质、地球物理、地球化学手段和计算机模拟技术，在盆地演化历史中和输导格架下，通过能量场演化及其控制的化学动力学、流体动力学和运动学过程分析，研究沉积盆地油气形成、演化和运移过程和聚集规律的综合性学科。成藏动力学研究的基础是盆地演化历史和流体输导格架，研究的核心是能量场(包括温度场、压力场、应力场)演化及其控制的化学动力学和流体动力学过程。近年来，由于油气勘探的深入和多学科联合研究的开展，成藏动力学在流体输导系统、盆地能量场演化与流体流动样式、油气成藏机理与充注历史分析等各个方面都取得了重要进展。
地层孔隙压力是沉积演化、构造运动、水动力和水－岩相互作用等综合作用的结果，是研究成藏过程、烃源岩演化、储层演化等成藏动力学方面必不可少的参数之一。现今地层孔隙压力可以通过钻杆测试获取，而古地层孔隙压力的求取却十分困难，这是因地质过程的不可逆特性，因此恢复地质过程相当困难。目前，可用于恢复古地层压力的方法有很多，如，根据流体包裹体均一温度和流体成分之间的平衡关系确定古压力(简称流体包裹体法)；根据黏土矿物形成温度及实际曲线估算黏土矿物的形成压力(简称黏土矿物法)；从泥岩声波时差资料出发，利用泥岩压实的不可逆原理，推导出古地层压力(简称声波时差法)以及盆地模拟法，其中盆地模拟法基于压力生成的物理-化学机理，既可开展一维单井井压力演化重建，也可开展二维剖面甚至三维立体的压力史模拟分析。
流体包裹体法、黏土矿物法和声波时差法，都可以计(估)算在某个地质时刻的古地层压力，无法恢复整个地质历史时期的演化过程。而盆地模拟法虽然可以从区域上反映古压力大小及其演化过程，但盆地模拟软件所采用的数学模型是对复杂地质情况的简化，如果不考虑流体包裹体法、黏土矿物法和声波时差法等方法确定的某个地质时刻的古压力和现今的实测压力，其结果必然与实际地质情况存在误差。
发明内容
本发明为了解决上述背景技术中的技术问题中的至少一个，提出了一种古气藏地层孔隙压力演化史曲线生成方法及装置。
为了实现上述目的，根据本发明的一个方面，提供了一种古气藏地层孔隙压力演化史曲线生成方法，该方法包括：
获取选定的包裹体的气液比、均一温度和冰点温度；
根据所述冰点温度确定所述包裹体的流体盐度；
根据所述流体盐度不断设定所述包裹体的成分，针对每次设定的成分，确定在所述均一温度条件下所述包裹体在最小捕获压力时的体积，以及确定在室温条件下通过调整压力使所述包裹体的体积与在最小捕获压力下的体积相等时所述包裹体的气液比，直至设定的所述包裹体的成分对应的体积相等时的所述气液比与所述包裹体的气液比相等；
针对最后设定的所述包裹体的成分确定所述包裹体在所述均一温度时的最小捕获压力，并根据该最小捕获压力确定所述包裹体对应的古气藏成藏时的古孔隙压力；
根据所述古孔隙压力以及当前测量的所述包裹体所在地层的孔隙压力确定所述包裹体对应的古气藏的孔隙压力演化史曲线。
可选的，所述根据所述古孔隙压力以及当前测量的所述包裹体所在地层的孔隙压力确定所述包裹体对应的古气藏的孔隙压力演化史曲线，具体包括：
根据所述包裹体所在单井的测井参数建立所述单井的地质模型，其中，所述测井参数包括：各地层的岩性、各地层的剥蚀量、镜质体反射率以及大地热流；
通过调节所述包裹体所在地层的岩性和/或所述包裹体上覆地层的岩性调整所述地质模型，使所述古孔隙压力和当前测量的所述包裹体所在地层的孔隙压力均落在所述地质模型模拟出的所述包裹体对应的古气藏的孔隙压力演化史曲线上。
可选的，该方法还包括：
通过调整所述各地层的剥蚀量和/或所述大地热流对所述地质模型进行修正，使所述地质模型模拟出的镜质体反射率模拟值与所述镜质体反射率相等。
可选的，所述根据所述冰点温度确定所述包裹体的流体盐度，具体包括：
根据所述冰点温度以及预设的包裹体的盐度与冰点温度的对应关系确定所述包裹体的流体盐度。
可选的，所述根据所述流体盐度不断设定所述包裹体的成分，具体包括：
通过设置不同的包裹体甲烷含量不断设定所述包裹体的成分。
为了实现上述目的，根据本发明的另一方面，提供了一种古气藏地层孔隙压力演化史曲线生成装置，该装置包括：
包裹体参数获取单元，用于获取选定的包裹体的气液比、均一温度和冰点温度；
流体盐度确定单元，用于根据所述冰点温度确定所述包裹体的流体盐度；
真实成分确定单元，用于根据所述流体盐度不断设定所述包裹体的成分，针对每次设定的成分，确定在所述均一温度条件下所述包裹体在最小捕获压力时的体积，以及确定在室温条件下通过调整压力使所述包裹体的体积与在最小捕获压力下的体积相等时所述包裹体的气液比，直至设定的所述包裹体的成分对应的体积相等时的所述气液比与所述包裹体的气液比相等；
古孔隙压力确定单元，用于针对最后设定的所述包裹体的成分确定所述包裹体在所述均一温度时的最小捕获压力，并根据该最小捕获压力确定所述包裹体对应的古气藏成藏时的古孔隙压力；
孔隙压力演化史曲线确定单元，用于根据所述古孔隙压力以及当前测量的所述包裹体所在地层的孔隙压力确定所述包裹体对应的古气藏的孔隙压力演化史曲线。
可选的，所述孔隙压力演化史曲线确定单元包括：
模型建立模块，用于根据所述包裹体所在单井的测井参数建立所述单井的地质模型，其中，所述测井参数包括：各地层的岩性、各地层的剥蚀量、镜质体反射率以及大地热流；
演化史曲线模拟模块，用于通过调节所述包裹体所在地层的岩性和/或所述包裹体上覆地层的岩性调整所述地质模型，使所述古孔隙压力和当前测量的所述包裹体所在地层的孔隙压力均落在所述地质模型模拟出的所述包裹体对应的古气藏的孔隙压力演化史曲线上。
可选的，所述孔隙压力演化史曲线确定单元还包括：
模型修正模块，用于通过调整所述各地层的剥蚀量和/或所述大地热流对所述地质模型进行修正，使所述地质模型模拟出的镜质体反射率模拟值与所述镜质体反射率相等。
可选的，所述流体盐度确定单元，具体用于根据所述冰点温度以及预设的包裹体的盐度与冰点温度的对应关系确定所述包裹体的流体盐度。
可选的，所述真实成分确定单元，具体用于通过设置不同的包裹体甲烷含量不断设定所述包裹体的成分。
为了实现上述目的，根据本发明的另一方面，还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述古气藏地层孔隙压力演化史曲线生成方法中的步骤。
为了实现上述目的，根据本发明的另一方面，还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序在计算机处理器中执行时实现上述古气藏地层孔隙压力演化史曲线生成方法中的步骤。
本发明的有益效果为：本发明针对现含油气盆地古压力恢复技术的缺陷和不足，提供一种设计合理、定量化的技术来计算古气藏地层孔隙压力演化过程的方法。在成藏动力学研究过程中，可以定量计算出一条连接成藏期压力值和现今压力值的单井孔隙压力演化曲线，提高了计算结果的可靠性和准确性，实用性更强。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1是本发明实施例古气藏地层孔隙压力演化史曲线生成方法的流程图；
图2是本发明实施例模拟孔隙压力演化史曲线的流程图；
图3是本发明实施例古气藏地层孔隙压力演化史曲线生成装置的结构框图；
图4是本发明实施例孔隙压力演化史曲线确定单元的组成结构框图；
图5是本发明实施例计算机设备示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
图1是本发明实施例古气藏地层孔隙压力演化史曲线生成方法的流程图，如图1所示，本实施例的古气藏地层孔隙压力演化史曲线生成方法包括步骤S101至步骤S105。
步骤S101，获取选定的包裹体的气液比、均一温度和冰点温度。
在本发明实施例中，对于古气藏来说，古油藏时期捕获的包裹体在后期埋藏过程中经历了较高温度后发生了裂解，因此油包裹体很难检测到，但是存在大量含烃盐水包裹体。常温下，这些气液两相的含烃盐水包裹体就是本发明所研究的对象，即本步骤选取的包裹体。进而本步骤可以在在冷热台上获得包裹体的均一温度和冰点温度。
在本发明实施例中，包裹体的气液比为包裹体气相和液相的面积比。在本发明的可选实施例中，确定气液比的方法可以为：利用共聚焦显微镜，在透光下找到所需要计算的包裹体，把它的图像输入计算机，利用共聚焦激光扫描显微镜附带的相关软件，计算出气相与液相的面积，两个面积的比值代表包裹体的实测气液比。
步骤S102，根据所述冰点温度确定所述包裹体的流体盐度。
在本发明可选实施例中，本步骤具体可以根据所述冰点温度以及预设的包裹体的盐度与冰点温度的对应关系确定所述包裹体的流体盐度。具体的本步骤利用Bodnar(1993)总结出的盐度-冰点温度关系表，查得包裹体流体盐度。
步骤S103，根据所述流体盐度不断设定所述包裹体的成分，针对每次设定的成分，确定在所述均一温度条件下所述包裹体在最小捕获压力时的体积，以及确定在室温条件下通过调整压力使所述包裹体的体积与在最小捕获压力下的体积相等时所述包裹体的气液比，直至设定的所述包裹体的成分对应的体积相等时的所述气液比与所述包裹体的气液比相等。
在本发明实施例中，本步骤用于较为准确的确定所选定的包裹体的成分。由于古气藏气体成分主要是甲烷，本发明可以将包裹体初步看作一定甲烷含量和一定盐度的混合液，由于包裹体流体盐度已由上述步骤S102确定，本步骤可以通过设置不同的甲烷含量来设定不同的包裹体的成分。
本步骤在第一次设定包裹体的成分(初始成分)后，可以确定该初始成分下的包裹体在均一温度条件下的最小捕获压力，以及在该最小捕获压力时的包裹体总体积(记为V0)。在本发明的可选实施例中，具体的可以利用PVTsim软件的“flash”选项中选择“multi-phase”，输入包裹体的均一温度。通过不断输入和调整压力P，使得输出结果显示的“Volume％”栏下的“Vapor”值刚好为零(表示包裹体中的气体完全转化为液态)，记录此时的压力P和包裹体总体积V0，此时的压力P为设定初始成分包裹体的最小捕获压力，可精确到0.1MPa。
进一步的，本步骤在室温条件(23至27度，优选25度)下，通过不断调整压力使该初始成分的包裹体的体积(V)与上述在最小捕获压力下的体积(V0)相等，并确定此时使该初始成分的包裹体的气液比。由于包裹体是一个封闭体系，在温度和压力变化的过程中，它的总体积不会变。当V＝V0时，这两点就会在同一条等容线上。在本发明的可选实施例中，具体的可以利用PVTsim软件的“flash”选项中选择“multi-phase”，输入室温(25℃)，通过不断输入和调整压力P，使得包裹体的体积V和上述V0相等，确定此时使该初始成分的包裹体的气液比。
进一步的，本步骤将该初始成分的包裹体对应的V＝V0时的气液比与包裹体的初始气液比(即步骤S101获取的包裹体的气液比)进行对比，若不相等，则重新设定包裹体的成分再进行本步骤中的上述步骤，直至设定的包裹体的成分对应的V＝V0时的气液比与包裹体的初始气液比相等，该相等时的包裹体的成分可以看作包裹体的“真实成分”。“真实成分”为该选定的包裹体相对准确的成分。
步骤S104，针对最后设定的所述包裹体的成分确定所述包裹体在所述均一温度时的最小捕获压力，并根据该最小捕获压力确定所述包裹体对应的古气藏成藏时的古孔隙压力。
在本发明实施例中，本步骤针对上述步骤S103确定的选定的包裹体的“真实成分”确定包裹体均一为液相时的最小压力，即包裹体的最小捕获压力。在本发明可选实施例中，该包裹体的最小捕获压力即为包裹体对应的古气藏成藏时的古孔隙压力。
步骤S105，根据所述古孔隙压力以及当前测量的所述包裹体所在地层的孔隙压力确定所述包裹体对应的古气藏的孔隙压力演化史曲线。
在本发明实施例中，本步骤在得知包裹体对应的古气藏成藏时的古孔隙压力以及当前测量的包裹体所在地层的孔隙压力后可以采用现有技术的演化史分析软件(例如PetroMod)模拟出该包裹体对应的古气藏的孔隙压力演化史曲线。
由以上描述可以看出，本发明针对现含油气盆地古压力恢复技术的缺陷和不足，提供一种设计合理、定量化的技术来计算古气藏地层孔隙压力演化过程的方法。在成藏动力学研究过程中，可以定量计算出一条连接成藏期压力值和现今压力值的单井孔隙压力演化曲线，提高了计算结果的可靠性和准确性，实用性更强。
图2是本发明实施例模拟孔隙压力演化史曲线的流程图，如图2所示，在本发明可选实施例中，上述步骤S105的根据所述古孔隙压力以及当前测量的所述包裹体所在地层的孔隙压力确定所述包裹体对应的古气藏的孔隙压力演化史曲线具体包括步骤S201至步骤S203。
步骤S201，根据所述包裹体所在单井的测井参数建立所述单井的地质模型。
在本发明可选实施例中，本步骤利用PetroMod 1D建立包裹体所在单井的地质模型。在本发明可选实施例中，测井参数具体包括：各地层界面对应的深度、各地层的岩性、烃源岩厚度、各地层界面对应的地质年龄、各地层的剥蚀厚度和剥蚀时间、各地层的剥蚀量、镜质体反射率(Ro)、大地热流以及地表温度等参数。
在本发明实施例中，通过PetroMod 1D建立包裹体所在单井的地质模型可以用于对该单井的一些测井参数进行模拟，以及对古气藏的孔隙压力演化史进行模拟，输出模拟的古气藏的孔隙压力演化史曲线。
步骤S202，通过调整所述各地层的剥蚀量和/或所述大地热流对所述地质模型进行修正，使所述地质模型模拟出的镜质体反射率模拟值与所述镜质体反射率相等。
在本发明实施例中，为了提高建立的地质模型的准确性，本发明还对地质模型进行修正。具体修正为：通过联合调整各地层的剥蚀量和大地热流对地质模型进行调整，使地质模型模拟出的镜质体反射率模拟值与镜质体反射率(Ro)相吻合。
步骤S203，通过调节所述包裹体所在地层的岩性和/或所述包裹体上覆地层的岩性调整所述地质模型，使所述古孔隙压力和当前测量的所述包裹体所在地层的孔隙压力均落在所述地质模型模拟出的所述包裹体对应的古气藏的孔隙压力演化史曲线上。
在本发明实施例中，通过PetroMod 1D建立的包裹体所在单井的地质模型用于对古气藏的孔隙压力演化史进行模拟，输出模拟的古气藏的孔隙压力演化史曲线。本步骤将包裹体对应的古气藏成藏时的古孔隙压力作为约束条件，当前测量的包裹体所在地层的孔隙压力作为边界条件，获得古气藏所在地层的孔隙压力演化史。具体的，本步骤通过不断的调节包裹体所在地层的岩性和包裹体上覆地层的岩性参数来调整地质模型，使地质模型输出的古气藏的孔隙压力演化史曲线同时穿过包裹体对应的古气藏成藏时的古孔隙压力和当前测量的包裹体所在地层的孔隙压力，此时地质模型输出的古气藏的孔隙压力演化史曲线为符合实际的压力演化史曲线。
由以上描述可以看出，本发明的有益效果为：本发明为一种基于储层含烃盐水包裹体热动力学参数结合相态模拟软件和盆地模拟技术定量恢复古气藏地层孔隙压力演化史的方法，克服了先前单一方法的局限性和不足，为更好的开展成藏动力学提供有利参数，实用性更强。本发明具有设计合理，实用性强等优点。
需要说明的是，在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。
基于同一发明构思，本发明实施例还提供了一种古气藏地层孔隙压力演化史曲线生成装置，可以用于实现上述实施例所描述的古气藏地层孔隙压力演化史曲线生成方法，如下面的实施例所述。由于古气藏地层孔隙压力演化史曲线生成装置解决问题的原理与古气藏地层孔隙压力演化史曲线生成方法相似，因此古气藏地层孔隙压力演化史曲线生成装置的实施例可以参见古气藏地层孔隙压力演化史曲线生成方法的实施例，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
图3是本发明实施例古气藏地层孔隙压力演化史曲线生成装置的结构框图，如图3所示，本发明实施例古气藏地层孔隙压力演化史曲线生成装置包括：包裹体参数获取单元1、流体盐度确定单元2、真实成分确定单元3、古孔隙压力确定单元4和孔隙压力演化史曲线确定单元5。
包裹体参数获取单元1，用于获取选定的包裹体的气液比、均一温度和冰点温度。
流体盐度确定单元2，用于根据所述冰点温度确定所述包裹体的流体盐度。
在本发明可选实施例中，流体盐度确定单元2具体用于根据所述冰点温度以及预设的包裹体的盐度与冰点温度的对应关系确定所述包裹体的流体盐度。
真实成分确定单元3，用于根据所述流体盐度不断设定所述包裹体的成分，针对每次设定的成分，确定在所述均一温度条件下所述包裹体在最小捕获压力时的体积，以及确定在室温条件下通过调整压力使所述包裹体的体积与在最小捕获压力下的体积相等时所述包裹体的气液比，直至设定的所述包裹体的成分对应的体积相等时的所述气液比与所述包裹体的气液比相等。
在本发明可选实施例中，真实成分确定单元3具体用于通过设置不同的包裹体甲烷含量不断设定所述包裹体的成分。
古孔隙压力确定单元4，用于针对最后设定的所述包裹体的成分确定所述包裹体在所述均一温度时的最小捕获压力，并根据该最小捕获压力确定所述包裹体对应的古气藏成藏时的古孔隙压力。
孔隙压力演化史曲线确定单元5，用于根据所述古孔隙压力以及当前测量的所述包裹体所在地层的孔隙压力确定所述包裹体对应的古气藏的孔隙压力演化史曲线。
图4是本发明实施例孔隙压力演化史曲线确定单元的组成结构框图，如图4所示，在本发明可选实施例中孔隙压力演化史曲线确定单元5具体包括：模型建立模块501、模型修正模块502和演化史曲线模拟模块503。
模型建立模块501，用于根据所述包裹体所在单井的测井参数建立所述单井的地质模型，其中，所述测井参数包括：各地层的岩性、各地层的剥蚀量、镜质体反射率以及大地热流。
模型修正模块502，用于通过调整所述各地层的剥蚀量和/或所述大地热流对所述地质模型进行修正，使所述地质模型模拟出的镜质体反射率模拟值与所述镜质体反射率相等。
演化史曲线模拟模块503，用于通过调节所述包裹体所在地层的岩性和/或所述包裹体上覆地层的岩性调整所述地质模型，使所述古孔隙压力和当前测量的所述包裹体所在地层的孔隙压力均落在所述地质模型模拟出的所述包裹体对应的古气藏的孔隙压力演化史曲线上。
为了实现上述目的，根据本申请的另一方面，还提供了一种计算机设备。如图5所示，该计算机设备包括存储器、处理器、通信接口以及通信总线，在存储器上存储有可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述实施例方法中的步骤。
处理器可以为中央处理器(Central Processing Unit，CPU)。处理器还可以为其他通用处理器、数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现场可编程门阵列(Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等芯片，或者上述各类芯片的组合。
存储器作为一种非暂态计算机可读存储介质，可用于存储非暂态软件程序、非暂态计算机可执行程序以及单元，如本发明上述方法实施例中对应的程序单元。处理器通过运行存储在存储器中的非暂态软件程序、指令以及模块，从而执行处理器的各种功能应用以及作品数据处理，即实现上述方法实施例中的方法。
存储器可以包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需要的应用程序；存储数据区可存储处理器所创建的数据等。此外，存储器可以包括高速随机存取存储器，还可以包括非暂态存储器，例如至少一个磁盘存储器件、闪存器件、或其他非暂态固态存储器件。在一些实施例中，存储器可选包括相对于处理器远程设置的存储器，这些远程存储器可以通过网络连接至处理器。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
所述一个或者多个单元存储在所述存储器中，当被所述处理器执行时，执行上述实施例中的方法。
上述计算机设备具体细节可以对应参阅上述实施例中对应的相关描述和效果进行理解，此处不再赘述。
为了实现上述目的，根据本申请的另一方面，还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序在计算机处理器中执行时实现上述古气藏地层孔隙压力演化史曲线生成方法中的步骤。本领域技术人员可以理解，实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述存储介质可为磁碟、光盘、只读存储记忆体(Read-OnlyMemory，ROM)、随机存储记忆体(RandomAccessMemory，RAM)、快闪存储器(Flash Memory)、硬盘(Hard Disk Drive，缩写：HDD)或固态硬盘(Solid-State Drive，SSD)等；所述存储介质还可以包括上述种类的存储器的组合。
显然，本领域的技术人员应该明白，上述的本发明的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (12)
1.一种古气藏地层孔隙压力演化史曲线生成方法，其特征在于，包括：
获取选定的包裹体的气液比、均一温度和冰点温度；
根据所述冰点温度确定所述包裹体的流体盐度；
根据所述流体盐度不断设定所述包裹体的成分，针对每次设定的成分，确定在所述均一温度条件下所述包裹体在最小捕获压力时的体积，以及确定在室温条件下通过调整压力使所述包裹体的体积与在最小捕获压力下的体积相等时所述包裹体的气液比，直至设定的所述包裹体的成分对应的体积相等时的所述气液比与所述包裹体的气液比相等；
针对最后设定的所述包裹体的成分确定所述包裹体在所述均一温度时的最小捕获压力，并根据该最小捕获压力确定所述包裹体对应的古气藏成藏时的古孔隙压力；
根据所述古孔隙压力以及当前测量的所述包裹体所在地层的孔隙压力确定所述包裹体对应的古气藏的孔隙压力演化史曲线。
2.根据权利要求1所述的古气藏地层孔隙压力演化史曲线生成方法，其特征在于，所述根据所述古孔隙压力以及当前测量的所述包裹体所在地层的孔隙压力确定所述包裹体对应的古气藏的孔隙压力演化史曲线，具体包括：
根据所述包裹体所在单井的测井参数建立所述单井的地质模型，其中，所述测井参数包括：各地层的岩性、各地层的剥蚀量、镜质体反射率以及大地热流；
通过调节所述包裹体所在地层的岩性和/或所述包裹体上覆地层的岩性调整所述地质模型，使所述古孔隙压力和当前测量的所述包裹体所在地层的孔隙压力均落在所述地质模型模拟出的所述包裹体对应的古气藏的孔隙压力演化史曲线上。
3.根据权利要求2所述的古气藏地层孔隙压力演化史曲线生成方法，其特征在于，还包括：
通过调整所述各地层的剥蚀量和/或所述大地热流对所述地质模型进行修正，使所述地质模型模拟出的镜质体反射率模拟值与所述镜质体反射率相等。
4.根据权利要求1所述的古气藏地层孔隙压力演化史曲线生成方法，其特征在于，所述根据所述冰点温度确定所述包裹体的流体盐度，具体包括：
根据所述冰点温度以及预设的包裹体的盐度与冰点温度的对应关系确定所述包裹体的流体盐度。
5.根据权利要求1所述的古气藏地层孔隙压力演化史曲线生成方法，其特征在于，所述根据所述流体盐度不断设定所述包裹体的成分，具体包括：
通过设置不同的包裹体甲烷含量不断设定所述包裹体的成分。
6.一种古气藏地层孔隙压力演化史曲线生成装置，其特征在于，包括：
包裹体参数获取单元，用于获取选定的包裹体的气液比、均一温度和冰点温度；
流体盐度确定单元，用于根据所述冰点温度确定所述包裹体的流体盐度；
真实成分确定单元，用于根据所述流体盐度不断设定所述包裹体的成分，针对每次设定的成分，确定在所述均一温度条件下所述包裹体在最小捕获压力时的体积，以及确定在室温条件下通过调整压力使所述包裹体的体积与在最小捕获压力下的体积相等时所述包裹体的气液比，直至设定的所述包裹体的成分对应的体积相等时的所述气液比与所述包裹体的气液比相等；
古孔隙压力确定单元，用于针对最后设定的所述包裹体的成分确定所述包裹体在所述均一温度时的最小捕获压力，并根据该最小捕获压力确定所述包裹体对应的古气藏成藏时的古孔隙压力；
孔隙压力演化史曲线确定单元，用于根据所述古孔隙压力以及当前测量的所述包裹体所在地层的孔隙压力确定所述包裹体对应的古气藏的孔隙压力演化史曲线。
7.根据权利要求6所述的古气藏地层孔隙压力演化史曲线生成装置，其特征在于，所述孔隙压力演化史曲线确定单元包括：
模型建立模块，用于根据所述包裹体所在单井的测井参数建立所述单井的地质模型，其中，所述测井参数包括：各地层的岩性、各地层的剥蚀量、镜质体反射率以及大地热流；
演化史曲线模拟模块，用于通过调节所述包裹体所在地层的岩性和/或所述包裹体上覆地层的岩性调整所述地质模型，使所述古孔隙压力和当前测量的所述包裹体所在地层的孔隙压力均落在所述地质模型模拟出的所述包裹体对应的古气藏的孔隙压力演化史曲线上。
8.根据权利要求7所述的古气藏地层孔隙压力演化史曲线生成装置，其特征在于，所述孔隙压力演化史曲线确定单元还包括：
模型修正模块，用于通过调整所述各地层的剥蚀量和/或所述大地热流对所述地质模型进行修正，使所述地质模型模拟出的镜质体反射率模拟值与所述镜质体反射率相等。
9.根据权利要求6所述的古气藏地层孔隙压力演化史曲线生成装置，其特征在于，所述流体盐度确定单元，具体用于根据所述冰点温度以及预设的包裹体的盐度与冰点温度的对应关系确定所述包裹体的流体盐度。
10.根据权利要求6所述的古气藏地层孔隙压力演化史曲线生成装置，其特征在于，所述真实成分确定单元，具体用于通过设置不同的包裹体甲烷含量不断设定所述包裹体的成分。
11.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至5任一项方法中的步骤。
12.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于，所述计算机程序在计算机处理器中执行时实现如权利要求1至5任意一项方法中的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911374034.6A CN113051697B (zh) | 2019-12-27 | 2019-12-27 | 古气藏地层孔隙压力演化史曲线生成方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911374034.6A CN113051697B (zh) | 2019-12-27 | 2019-12-27 | 古气藏地层孔隙压力演化史曲线生成方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113051697A true CN113051697A (zh) | 2021-06-29 |
| CN113051697B CN113051697B (zh) | 2022-10-04 |
Family
ID=76506253
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911374034.6A Active CN113051697B (zh) | 2019-12-27 | 2019-12-27 | 古气藏地层孔隙压力演化史曲线生成方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113051697B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115081358A (zh) * | 2022-07-22 | 2022-09-20 | 核工业北京地质研究院 | 确定砂岩铀矿渗出流体的迁移动力和迁移时间的方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103982179A (zh) * | 2014-05-26 | 2014-08-13 | 中国地质大学（北京） | 一种油藏储层的古压力定量反演探测方法 |
| CN107966443A (zh) * | 2017-12-25 | 2018-04-27 | 中国石油大学(华东) | 一种地层古流体盐度演化史的研究方法 |
| CN110259439A (zh) * | 2019-05-30 | 2019-09-20 | 中国石油大学(北京) | 用于评价碳酸盐岩断层封闭性演化史的方法及装置 |
| CN110517794A (zh) * | 2019-08-23 | 2019-11-29 | 长安大学 | 一种建立页岩气藏的埋藏-热演化史图的方法 |
-
2019
- 2019-12-27 CN CN201911374034.6A patent/CN113051697B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103982179A (zh) * | 2014-05-26 | 2014-08-13 | 中国地质大学（北京） | 一种油藏储层的古压力定量反演探测方法 |
| CN107966443A (zh) * | 2017-12-25 | 2018-04-27 | 中国石油大学(华东) | 一种地层古流体盐度演化史的研究方法 |
| CN110259439A (zh) * | 2019-05-30 | 2019-09-20 | 中国石油大学(北京) | 用于评价碳酸盐岩断层封闭性演化史的方法及装置 |
| CN110517794A (zh) * | 2019-08-23 | 2019-11-29 | 长安大学 | 一种建立页岩气藏的埋藏-热演化史图的方法 |
Non-Patent Citations (2)
| Title |
|---|
| 廖玲玲等: "利用包裹体PVTx和PetroMod盆地模拟系统反演鄂尔多斯盆地上古生界古流体压力", 《矿物岩石地球化学通报》 * |
| 张俊武等: "含烃盐水包裹体PVT模拟新方法及其在气藏古压力恢复中的应用", 《石油实验地质》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115081358A (zh) * | 2022-07-22 | 2022-09-20 | 核工业北京地质研究院 | 确定砂岩铀矿渗出流体的迁移动力和迁移时间的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113051697B (zh) | 2022-10-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Chai et al. | An efficient method for fractured shale reservoir history matching: The embedded discrete fracture multi-continuum approach | |
| CA3007950C (en) | Predicting and modeling changes in capillary pressure and relative permeabilities in a porous medium due to mineral precipitation and dissolution | |
| AU2013397958B2 (en) | A simulation-to-seismic workflow construed from core based rock typing and enhanced by rock replacement modeling | |
| NO335854B1 (no) | Automatisert system for modellering av flerverdihorisonter med forkastninger | |
| NO20111055A1 (no) | Interaktiv strukturell rekonstruksjon ved tolkning av seismikkvolumer for struktur og stratigrafi | |
| EP3096252A2 (en) | Adaptive multiscale multi-fidelity reservoir simulation | |
| Zhang et al. | A new practical method for upscaling in highly heterogeneous reservoir models | |
| Rupke et al. | Automated thermotectonostratigraphic basin reconstruction: Viking Graben case study | |
| WO2015038163A1 (en) | Pseudo-phase production simulation: a signal processing approach to assess quasi-multiphase flow production via successive analogous step-function relative permeability controlled models in reservoir flow simulation | |
| CN107808032A (zh) | 页岩气储层含水饱和度计算方法及系统 | |
| AU2013399193B2 (en) | Static earth model grid cell scaling and property re-sampling methods and systems | |
| CA2921390C (en) | Pseudo phase production simulation: a signal processing approach to assess quasi-multiphase flow production via successive analogous step-function relative permeability controlled models in reservoir flow simulation in order to rank multiple petro-physical realizations | |
| NO20190509A1 (en) | Global surface paleo-temperature modeling tool | |
| Onishi et al. | Fast Marching Method Based Rapid Simulation Accounting for Gravity | |
| AU2013400127B2 (en) | Relative permeability inversion from historical production data using viscosity ratio invariant step-function relative permeability approximations | |
| CN113051697B (zh) | 古气藏地层孔隙压力演化史曲线生成方法及装置 | |
| Correia et al. | Flow simulation using local grid refinements to model laminated reservoirs | |
| Couples | Phenomenological understanding of poroelasticity via the micromechanics of a simple digital-rock model | |
| Becerra et al. | Petroleum reservoir uncertainty mitigation through the integration with production history matching | |
| Dakhelpour-Ghoveifel et al. | Prediction of gas-oil capillary pressure of carbonate rock using pore network modeling | |
| Hui et al. | An Integrated EDFM+ DPDK Hybrid Workflow for Hierarchical Treatment of Fractures in Practical Field Studies | |
| CN117094152B (zh) | 一种耦合沉积作用的盆地沉降史模拟方法及装置 | |
| Lajevardi et al. | Stochastic regridding of geological models for flow simulation | |
| Katterbauer et al. | Data assimilation of InSAR measurements for large scale reservoirs | |
| Onishi | Effective Reservoir Management for Unconventional Reservoirs Using the Fast Marching Method and Machine Learning |
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 |