CN106844858A - 地层裂缝发育区带预测方法及装置 - Google Patents
地层裂缝发育区带预测方法及装置 Download PDFInfo
- Publication number
- CN106844858A CN106844858A CN201611190644.7A CN201611190644A CN106844858A CN 106844858 A CN106844858 A CN 106844858A CN 201611190644 A CN201611190644 A CN 201611190644A CN 106844858 A CN106844858 A CN 106844858A
- Authority
- CN
- China
- Prior art keywords
- structural
- recent structural
- model
- construction
- recent
- 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
- 208000010392 Bone Fractures Diseases 0.000 title claims abstract description 77
- 206010017076 Fracture Diseases 0.000 title claims abstract description 77
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 76
- 238000011161 development Methods 0.000 title claims abstract description 69
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000010276 construction Methods 0.000 claims abstract description 124
- 238000003325 tomography Methods 0.000 claims description 41
- 239000011435 rock Substances 0.000 claims description 23
- 230000008859 change Effects 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 238000010998 test method Methods 0.000 claims description 4
- 238000011156 evaluation Methods 0.000 abstract description 3
- 230000000694 effects Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 230000033001 locomotion Effects 0.000 description 5
- 241001045988 Neogene Species 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000013508 migration Methods 0.000 description 4
- 230000005012 migration Effects 0.000 description 4
- 101150091879 neo gene Proteins 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 241000208340 Araliaceae Species 0.000 description 3
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 3
- 235000003140 Panax quinquefolius Nutrition 0.000 description 3
- 230000009471 action Effects 0.000 description 3
- 230000008602 contraction Effects 0.000 description 3
- 235000008434 ginseng Nutrition 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 239000004215 Carbon black (E152) Substances 0.000 description 2
- 229930195733 hydrocarbon Natural products 0.000 description 2
- 150000002430 hydrocarbons Chemical class 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000010008 shearing Methods 0.000 description 2
- 239000004575 stone Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000003556 assay Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000010429 evolutionary process Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000001125 extrusion Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000012395 formulation development Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000005226 mechanical processes and functions Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000037361 pathway Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- 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
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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/06—Energy or water supply
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- Computer Hardware Design (AREA)
- Theoretical Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Game Theory and Decision Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Quality & Reliability (AREA)
- Geometry (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本申请提供一种地层裂缝发育区带预测方法及装置，其中，该方法包括：建立待研究区域内现今构造的现今地质模型和现今构造所经历的各个演化阶段构造的各个古构造地质模型，获取现今构造的现今地应力和各个演化阶段构造的各个古地应力，获取现今构造以及各个演化阶段构造分别对应构造特征的物理力学参数，进而将现今地质模型和各个古构造地质模型分别转换为现今构造力学模型和各个古构造力学模型，并根据现今构造力学模型和各个古构造力学模型，预测待研究区域内现今构造的地层裂缝发育区带。该技术方案综合考虑现今地应力和各个古地应力，能够完全反映地层裂缝的实际发育状况，提高了预测地层裂缝发育区带的精度。
Description
技术领域
本申请涉及石油勘探技术领域，尤其涉及一种地层裂缝发育区带预测方法及装置。
背景技术
裂缝是油气的存储空间和运移通道，其分布与发育状况对油气成藏和运移有着深刻的影响。在石油勘探开发前期，准确预测地层裂缝的分布和发育状况有助于确定勘探靶区和制定开发方案，所以，地层裂缝发育区带预测属于油气勘探中的一项重要研究内容。
目前，地层裂缝区带预测方法主要通过对待研究区域的构造建立地质模型，并根据该构造的特征和现今地应力对已建立的地质模型进行分析处理，最后根据分析处理的结果来预测待研究区域内构造的地层裂缝发育情况。
然而，由于地层裂缝的空间展布规律性差，待研究区域内构造的现今地应力并不能完全反映地层裂缝的实际发育状况，故上述仅考虑现今地应力的地层裂缝预测方法的预测准确率低。
发明内容
本申请提供一种地层裂缝发育区带预测方法及装置，以解决现有地层裂缝预测方法的预测准确率低的问题。
本申请提供的一种地层裂缝发育区带预测方法，包括：
建立待研究区域内现今构造的现今地质模型和所述现今构造所经历的各个演化阶段构造的各个古构造地质模型；
获取所述现今构造的现今地应力和所述各个演化阶段构造的各个古地应力；
获取所述现今构造以及所述各个演化阶段构造分别对应构造特征的物理力学参数；
根据所述现今构造以及所述各个演化阶段构造分别对应构造特征的物理力学参数以及所述现今地应力和所述各个古地应力，将所述现今地质模型和所述各个古构造地质模型分别转换为现今构造力学模型和各个古构造力学模型；
根据所述现今构造力学模型和所述各个古构造力学模型，预测所述待研究区域内所述现今构造的地层裂缝发育区带。
本申请还提供一种地层裂缝发育区带预测装置，包括：
模型建立模块，用于建立待研究区域内现今构造的现今地质模型和所述现今构造所经历的各个演化阶段构造的各个古构造地质模型；
地应力获取模块，用于获取所述现今构造的现今地应力和所述各个演化阶段构造的各个古地应力；
力学参数获取模块，用于获取所述现今构造以及所述各个演化阶段构造分别对应构造特征的物理力学参数；
模型转换模块，用于根据所述现今构造以及所述各个演化阶段构造分别对应构造特征的物理力学参数以及所述现今地应力和所述各个古地应力，将所述现今地质模型和所述各个古构造地质模型分别转换为现今构造力学模型和各个古构造力学模型；
裂缝预测模块，用于根据所述现今构造力学模型和所述各个古构造力学模型，预测所述待研究区域内所述现今构造的地层裂缝发育区带。
本申请提供的地层裂缝发育区带预测方法及装置，通过建立待研究区域内现今构造的现今地质模型和该现今构造所经历的各个演化阶段构造的各个古构造地质模型，获取现今构造的现今地应力和各个演化阶段构造的各个古地应力，获取现今构造以及各个演化阶段构造分别对应构造特征的物理力学参数，以及根据现今构造以及各个演化阶段构造分别对应构造特征的物理力学参数以及现今地应力和各个古地应力，将现今地质模型和各个古构造地质模型分别转换为现今构造力学模型和各个古构造力学模型，根据现今构造力学模型和各个古构造力学模型，预测待研究区域内现今构造的地层裂缝发育区带。本申请实施例的技术方案综合考虑现今地应力和各个古地应力，能够完全反映地层裂缝的实际发育状况，提高了预测地层裂缝发育区带的精度。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本申请提供的地层裂缝发育区带预测方法实施例一的流程示意图；
图2为本申请提供的地层裂缝发育区带预测方法实施例二的流程示意图；
图3A至图3E分别为四个演化阶段构造以及现今构造的地质模型示意图；
图4为本申请提供的地层裂缝发育区带预测方法实施例三的流程示意图；
图5为本申请提供的地层裂缝发育区带预测方法实施例四的流程示意图；
图6为本申请提供的地层裂缝发育区带预测方法实施例五的流程示意图；
图7A至图7E为各个古构造力学模型和现今构造力学模型中地层裂缝发育区带分布的示意图；
图8为本申请提供的地层裂缝发育区带预测装置实施例一的结构示意图；
图9为本申请提供的地层裂缝发育区带预测装置实施例二的结构示意图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
近年来，随着油气田勘探程度的深入，构造简单并且容易开采的油气藏逐渐减少，人们不得不去开采那些构造复杂的油气藏。但是这些构造复杂的油气藏经历了多次构造运动，地层变形严重，许多地层倾角变陡，并且产生了大量裂缝。对于油气藏而言，这些地区裂缝不仅是重要的油气储集空间和渗流通道，同时还是优质储层发育带区的主控因素之一，因此，准确预测地层裂缝的空间展布是油气勘探中一项重要的研究内容。
由于地应力、岩性、断层，以及沉积环境等因素均会不同程度地影响地层裂缝的发育状况。地层历经漫长的地质年代，遭受了极为复杂的变形和重组，裂缝作为这种变化的产物而存在于地层中，影响油气的运移和油气藏的分布。
此外，从力学机制上讲，裂缝是地质材料在地应力作用下发生变形破坏的产物。所以，现今构造中的裂缝发育带，不仅受控于现今地应力场的作用，而且具有历史继承性，同样受各演化阶段构造的地应力作用影响，而现有地层裂缝预测方法没有考虑地层裂缝发育的地质历史继承性，致使预测准确率低。
本申请提供一种地层裂缝发育区带预测方法及装置，用于解决现有地层裂缝预测方法的预测准确率低的问题。下面，通过具体实施例对本申请的技术方案进行详细说明。
需要说明的是，下面这几个具体的实施例可以相互结合，对于相同或相似的概念或过程可能在某些实施例中不再赘述。
图1为本申请提供的地层裂缝发育区带预测方法实施例一的流程示意图。如图1所示，本申请实施例提供的地层裂缝发育区带预测方法包括以下步骤：
步骤11：建立待研究区域内现今构造的现今地质模型和该现今构造所经历的各个演化阶段构造的各个古构造地质模型。
在本实施例中，若要预测待研究区域内地层裂缝发育区带的分布情况，首先确定出待研究区域、该待研究区域内的现今构造以及该现今构造在形成过程中所经历的各个演化阶段对应的各个演化阶段构造。其次，根据实践应用中收集的构造解释资料，针对现今构造建立现今地质模型、针对该现今构造所经历的各个演化阶段构造，分别建立对应的各个古构造地质模型。
步骤12：获取上述现今构造的现今地应力和各个演化阶段构造的各个古地应力。
地应力是指存在于地壳中的应力，具体是由于岩石形变而引起的介质内部单位面积上的作用力。地应力在不同的阶段地应力不同，其中，现今地应力是指待研究区域中存在的或正在活动的地应力，古地应力是指某一地质时期或某一重要地质事件以前的地应力，在本实施例中，现今地应力指现今构造的地应力，各个古地应力分别与各个演化阶段相对应。
可选的，为了准确分析待研究区域内地层裂缝发育区带的分布情况，本申请实施例分别获取现今构造的现今地应力和各个演化阶段构造的各个古地应力，具体的，包括获取现今地应力的方向和大小以及各个古地应力的方向和大小。
步骤13：获取现今构造以及各个演化阶段构造分别对应构造特征的物理力学参数。
构造特征指构造内部构造层的叠加数以及每个构造层的构造形态特征。具体的，本申请实施例中构造特征主要指地层层序、断层分布和断层大小。
可选的，在本申请实施例中，待研究区域的构造在不同演化阶段的构造特征不同，所以，首先确定出现今构造和各个演化阶段构造包括的构造特征，并且测出各个构造特征的物理力学参数。
具体的，本实施例主要是获取现今构造对应地层岩石的物理力学参数以及断层介质的物理力学参数，以及各个演化阶段构造对应地层岩石的物理力学参数以及断层介质的物理力学参数。
步骤14：根据现今构造以及各个演化阶段构造分别对应构造特征的物理力学参数以及现今地应力和各个古地应力，将现今地质模型和各个古构造地质模型分别转换为现今构造力学模型和各个古构造力学模型。
可选的，为了便于分析待研究区域的地层裂缝发育状态，本申请实施例将采用上述获取到的现今构造以及各个演化阶段构造分别对应构造特征的物理力学参数对现今地质模型和各个古构造地质模型进行离散化处理，将现今地质模型离散化为现今构造力学模型以及将各个古构造地质模型离散化为各个古构造力学模型，该现今构造力学模型和各个古构造力学模型均利用后续分析计算。
此外，将上述现今地应力的方向和大小作为现今构造力学模型的边界条件加载到现今构造力学模型上，将各个古地应力的方向和大小作为各个古构造力学模型的边界条件加载到各个古构造力学模型上，实现现今地质模型和各个古构造地质模型的离散化处理。
步骤15：根据现今构造力学模型和各个古构造力学模型，预测该待研究区域内现今构造的地层裂缝发育区带。
在本实施例中，由于现今构造力学模型综合了现今构造对应构造特征的物理力学参数和现今地应力，各个古构造力学模型综合了各个演化阶段构造对应构造特征的物理力学参数和各个古地应力，故综合分析现今构造力学模型和各个古构造力学模型能够准确预测待研究区域内的现今构造的演化发展历程，进而可预测出该待研究区域内现今构造的地层裂缝发育区带。
本申请实施例提供的地层裂缝发育区带预测方法，通过建立待研究区域内现今构造的现今地质模型和该现今构造所经历的各个演化阶段构造的各个古构造地质模型，获取现今构造的现今地应力和各个演化阶段构造的各个古地应力，获取现今构造以及各个演化阶段构造分别对应构造特征的物理力学参数，以及根据现今构造以及各个演化阶段构造分别对应构造特征的物理力学参数以及现今地应力和各个古地应力，将现今地质模型和各个古构造地质模型分别转换为现今构造力学模型和各个古构造力学模型，根据现今构造力学模型和各个古构造力学模型，预测待研究区域内现今构造的地层裂缝发育区带。本申请实施例的技术方案综合考虑现今地应力和各个古地应力，能够完全反映地层裂缝的实际发育状况，提高了预测地层裂缝发育区带的精度，可用于辅助判断油气的运移方向和勘探靶区的确定。
进一步的，在上述实施例的基础上，图2为本申请提供的地层裂缝发育区带预测方法实施例二的流程示意图。如图2所示，本申请实施例提供的地层裂缝发育区带预测方法，在上述步骤11(建立待研究区域内现今构造的现今地质模型和现今构造所经历的各个演化阶段构造的各个古构造地质模型)之前，还包括：
步骤21：解析待研究区域，确定待研究区域的构造特征以及现今构造所经历的多个演化阶段。
具体的，由于地层经历了漫长的地质年代，所以，待研究区域经历了多个不同的演化阶段，并且在不同的演化阶段内具有不同的构造特征。
举例来说，对于待研究区域内的构造A，经过分析可知，该构造A长期受构造挤压作用，主要以其经历五个构造演化过程的影响进行说明，故待研究区域内的构造共经历侏罗纪末期构造、白垩纪末期构造、古近纪构造、新近纪构造、第四纪构造(现今构造)等过程。
相应的，图3A至图3E分别为四个演化阶段构造以及现今构造的地质模型示意图。
在本实施例中，侏罗纪末期受印支运动影响构造A在近南北向区域挤压背景下发生区域隆升，形成侏罗纪末期构造，建立的地质模型如图3A所示。
白垩纪受燕山中期运动影响，区域收缩构造变形导致地层发生强烈挤压收缩变形，构造A进一步隆升，形成白垩纪末期构造，建立的地质模型如图3B所示。
古近纪受燕山晚期运动影响，持续受到近南北向的挤压作用，形成古近纪构造，建立的地质模型如图3C所示。
古近纪受喜马拉雅早期运动影响，区域收缩构造变形导致地层发生较为强烈的挤压收缩变形，形成新近纪构造，建立的地质模型如图3D所示。
第四纪构造受喜马拉雅晚期运动影响，在区域近南北向挤压作用下，构造A发生收缩构造变形，背斜曲率增大，建立的现今地质模型如图3E所示。
值得说明的是，在图中，字母Q表示第四纪构造(现今构造)的地层，N表示新近纪构造的地层，E表示古近纪构造的地层，K表示白垩纪末期构造的地层，J表示侏罗纪末期构造的地层，以及T表示三叠系，P表示二叠系，C表示石炭系，D表示泥盆系。
此外，图3A-图3E为构造A南北走向剖面，在图3A中，南北方向长38.5km，深度8km；在图3B中，南北方向长38km，深度9km；在图3C中，南北方向长37km，深度11km；在图3D中，南北方向长36km，深度14km；在图3E中，南北方向长36km，深度15km。
步骤22：根据待研究区域的构造特征和现今构造所经历的多个演化阶段，确定影响待研究区域内地应力分布的地层界面和断层。
本实施例中，若要建立现今地质模型和各个古构造地质模型，首先解析待研究区域的构造特征和演化过程，分析待研究区域的构造特征和现今构造所经历的多个演化阶段，进而确定出能够影响待研究区域内地应力分布的地层界面和断层。此处，地层界面和断层的分布和大小在每个演化阶段均可能不同，因此，本申请实施例确定的是待研究区域内每个演化阶段影响的地应力分布的地层界面和断层。
本申请实施例提供的地层裂缝发育区带预测方法，通过解析待研究区域，确定待研究区域的构造特征以及现今构造所经历的多个演化阶段，并根据待研究区域的构造特征和现今构造所经历的多个演化阶段，确定影响待研究区域内地应力分布的地层界面和断层，这样为后续建立地质模型和获取构造特征的物理学参数奠定了基础。
进一步的，在图2所示实施例的基础上，上述步骤12(获取上述现今构造的现今地应力和各个演化阶段构造的各个古地应力)可通过如下可能实现方式实现，具体参见图4所示实施例。
图4为本申请提供的地层裂缝发育区带预测方法实施例三的流程示意图。如图4所示，上述步骤12(获取上述现今构造的现今地应力和各个演化阶段构造的各个古地应力)，包括：
步骤41：分别获取现今构造的现今地层岩芯和各个演化阶段构造的各个古地层岩芯。
步骤42：利用地层岩芯的古地磁定向和岩石声发射特性，分别测定现今构造的现今地应力和各个演化阶段构造的各个古地应力。
可选的，若想获取现今构造的现今地应力和各个演化阶段构造的各个古地应力，首先获取待研究区域内现今构造和各个演化阶段构造的地层岩芯，利用古地磁定向和岩石声发射特性测定现今构造的现今地应力的方向和大小以及各个演化阶段构造对应古地应力的方向和大小。
表1为构造A采用声发射特性测量出的最大有效地应力值。本实施例中利用地层岩芯古地磁定向结合岩石声发射特性测定的现今地应力和古地应力的大小如表1所示。可选的，现今地应力和各个古地应力的方向均为近南北向挤压应力。
表1构造A采用声发射特性测量出的最大有效地应力值
| 地层时代 | 最大有效地应力/MPa |
| Q | 58.3 |
| N | 65.2 |
| E | 48.6 |
| K | 53.6 |
| J | 38.6 |
本申请实施例提供的地层裂缝发育区带预测方法，在获取现今地应力和各个古地应力的方向和大小时，通过获取现今构造的现今地层岩芯和各个演化阶段构造的各个古地层岩芯，进而利用地层岩芯的古地磁定向和岩石声发射特性，分别测定现今构造的现今地应力和各个演化阶段构造的各个古地应力，测定方法简单、准确率高。
进一步的，图5为本申请提供的地层裂缝发育区带预测方法实施例四的流程示意图。如图5所示，上述步骤13(获取现今构造以及各个演化阶段构造分别对应构造特征的物理力学参数)，包括：
步骤51：分别获取现今构造和各个演化阶段构造在地层界面的地层岩芯以及现今构造和各个演化阶段构造对应断层的断层介质。
步骤52：利用岩石力学测试方法，获取现今构造和各个演化阶段构造分别对应的各个地层界面的地层岩芯的物理力学参数以及现今构造和各个演化阶段构造对应断层的断层介质的物理力学参数。
可选的，通过获取不同地质年代对应构造的地层岩芯，具体的，分别利用钻孔设备在现今构造和各个演化阶段构造对应地层界面的地层岩芯或者对应断层的断层介质，通过开展岩石静态物理力学参数测试，或利用测井资料计算获取现今构造和各个演化阶段构造对应断层的断层介质的物理力学参数。可选的，岩石静态物理力学参数测试方式的具体实现步骤参考国际岩石力学学会《岩石力学试验建议方法》，此处不再赘述。利用测井资料获取时，首先获取纵横波声波时差，通过一系列对应的计算公式用来计算岩石和断层介质的物理力学参数。
举例来说，表2现今构造和各个演化阶段构造对应的地层岩芯和断层介质的物理力学参数。使用岩石静态物理力学参数测试测定现今构造和各个演化阶段构造对应断层的断层介质的物理力学参数，如表2所示。
表2现今构造和各个演化阶段构造对应的地层岩芯和断层介质的物理力学参数
| 地层时代 | 岩性 | 弹性模量/Pa | 泊松比 | 体积模量/Pa | 剪切模量/Pa | 黏聚力/Pa | 内摩擦角/° | 抗拉强度/Pa | |
| Q | 粉砂岩 | 2460 | 0.200 | 38 | |||||
| N | 泥岩 | 2461 | 0.260 | 30 | 6.05×105 | ||||
| E | 砂质泥岩 | 2510 | 0.147 | 36 | |||||
| K | 细砂岩 | 2873 | 0.235 | 42 | |||||
| J | 砂岩 | 2487 | 0.123 | 40 | |||||
| T | 泥岩 | 2483 | 0.204 | 32 | |||||
| P | 粉砂岩 | 2460 | 0.200 | 38 | |||||
| C | 砂岩 | 2580 | 0.159 | 42 | |||||
| D | 砂质泥岩 | 2530 | 0.147 | 40 | |||||
| 1号断层 | 断层介质1 | 1320 | 0.178 | 30 | |||||
| 2号断层 | 断层介质2 | 1300 | 0.182 | 30 |
如表2所示，物理力学参数包括每个地层的岩性、岩芯密度、弹性模量、泊松比、体积模量、剪切模量、黏聚力、内摩擦角和抗拉强度等。
可选的，图6为本申请提供的地层裂缝发育区带预测方法实施例五的流程示意图。如图6所示，上述步骤14(根据现今构造以及各个演化阶段构造分别对应构造特征的物理力学参数以及现今地应力和各个古地应力，将现今地质模型和各个古构造地质模型分别转换为现今构造力学模型和各个古构造力学模型)，包括：
步骤61：利用现今构造对应构造特征的物理力学参数将现今地质模型离散为现今构造力学模型，上述现今地应力作为现今构造力学模型的边界条件。
步骤62：利用各个演化阶段构造分别对应构造特征的物理力学参数将各个古构造地质模型离散为各个古构造力学模型，上述各个古地应力分别作为各个古构造力学模型的边界条件。
作为一种示例，利用数值方法将图1所示实施例得到的现今地质模型和各个古地质模型离散化为现今构造力学模型和各个古构造力学模型。根据上述步骤得到的现今构造和各个演化阶段构造分别对应的各个地层界面的地层岩芯和对应断层的断层介质的物理力学参数，通过有限元法将上述现今地质模型和古地质模型进行离散化处理，分别离散为力学计算模型，以及将上述现今地应力作为现今构造力学模型的边界条件，将各个古地应力分别作为各个古构造力学模型的边界条件。
在本实施例中现今构造力学模型和各个古构造力学模型在南北方向受构造挤压作用，在竖直方向受重力作用。
对于上述构造A，图7A至图7E为各个古构造力学模型和现今构造力学模型中地层裂缝发育区带分布的示意图。如图7A至图7E所示，由于待研究区域受到多期南北向挤压作用，造成了多期地层岩石的变形破坏，形成了丰富的裂缝发育区带，图中横向线条表示地层界面，图中填充斜线的长条矩形、填充点、线条的不规则形状等均表示裂缝发育区带的示意图。
此外，本申请实施例并不限定裂缝发育区带的标识形状，只要是能够表示裂缝发育区带的表示方式均可以被采用，此处不再赘述。
值得说明的是，在图中，字母Q表示第四纪构造(现今构造)的地层，N表示新近纪构造的地层，E表示古近纪构造的地层，K表示白垩纪末期构造的地层，J表示侏罗纪末期构造的地层，以及T表示三叠系，P表示二叠系，C表示石炭系，D表示泥盆系。
本申请实施例提供的地层裂缝发育区带预测方法，通过模拟古地应力与现今地应力对各个演化阶段构造与现今构造的力学作用，考察不同演化阶段构造与现今构造地层中裂缝的发育状况，能够准确预测地层裂缝发育区带的分布，填补了现有地层裂缝发育区带预测方法的不足，为油气勘探开发提供技术支持。
下述为本申请装置实施例，可以用于执行本申请方法实施例。对于本申请装置实施例中未披露的细节，请参照本申请方法实施例。
图8为本申请提供的地层裂缝发育区带预测装置实施例一的结构示意图。如图8所示，本申请实施例提供的地层裂缝发育区带预测装置，包括：
模型建立模块81，用于建立待研究区域内现今构造的现今地质模型和所述现今构造所经历的各个演化阶段构造的各个古构造地质模型.
地应力获取模块82，用于获取所述现今构造的现今地应力和所述各个演化阶段构造的各个古地应力。
力学参数获取模块83，用于获取所述现今构造以及所述各个演化阶段构造分别对应构造特征的物理力学参数。
模型转换模块84，用于根据所述现今构造以及所述各个演化阶段构造分别对应构造特征的物理力学参数以及所述现今地应力和所述各个古地应力，将所述现今地质模型和所述各个古构造地质模型分别转换为现今构造力学模型和各个古构造力学模型。
裂缝预测模块85，用于根据所述现今构造力学模型和所述各个古构造力学模型，预测所述待研究区域内所述现今构造的地层裂缝发育区带。
本申请实施例提供的地层裂缝发育区带预测装置，可用于执行如图1所示方法实施例的技术方案，其实现原理和技术效果类似，此处不再赘述。
进一步的，在上述实施例的基础上，图9为本申请提供的地层裂缝发育区带预测装置实施例二的结构示意图。如图9所示，本申请实施例提供的地层裂缝发育区带预测装置，还包括：解析模块91和确定模块92。
该解析模块91，用于在模型建立模块81建立待研究区域内现今构造的现今地质模型和所述现今构造所经历的各个演化阶段构造的各个古构造地质模型之前，解析待研究区域，确定所述待研究区域的构造特征以及所述现今构造所经历的多个演化阶段。
该确定模块92，用于根据所述待研究区域的构造特征和所述现今构造所经历的多个演化阶段，确定影响所述待研究区域内地应力分布的地层界面和断层。
可选的，上述地应力获取模块82，具体用于分别获取所述现今构造的现今地层岩芯和所述各个演化阶段构造的各个古地层岩芯，利用地层岩芯的古地磁定向和岩石声发射特性，分别测定所述现今构造的现今地应力和所述各个演化阶段构造的各个古地应力。
在一实施例中，力学参数获取模块83，具体用于分别获取所述现今构造和所述各个演化阶段构造在所述地层界面的地层岩芯以及所述现今构造和所述各个演化阶段构造对应断层的断层介质，利用岩石力学测试方法，获取所述现今构造和所述各个演化阶段构造分别对应的各个地层界面的地层岩芯的物理力学参数以及所述现今构造和所述各个演化阶段构造对应断层的断层介质的物理力学参数。
在上述任一实施例中，模型转换模块84，具体用于利用所述现今构造对应构造特征的物理力学参数将所述现今地质模型离散为所述现今构造力学模型，所述现今地应力作为所述现今构造力学模型的边界条件，利用各个演化阶段构造分别对应构造特征的物理力学参数将所述各个古构造地质模型离散为所述各个古构造力学模型，所述各个古地应力分别作为所述各个古构造力学模型的边界条件。
本申请实施例提供的地层裂缝发育区带预测装置，可用于上述方法实施例的技术方案，其实现原理和技术效果类似，此处不再赘述。
本领域普通技术人员可以理解：实现上述各方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成。前述的程序可以存储于一计算机可读取存储介质中。该程序在执行时，执行包括上述各方法实施例的步骤；而前述的存储介质包括：ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
最后应说明的是：以上各实施例仅用以说明本申请的技术方案，而非对其限制；尽管参照前述各实施例对本申请进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本申请各实施例技术方案的范围。
Claims (10)
1.一种地层裂缝发育区带预测方法，其特征在于，包括：
建立待研究区域内现今构造的现今地质模型和所述现今构造所经历的各个演化阶段构造的各个古构造地质模型；
获取所述现今构造的现今地应力和所述各个演化阶段构造的各个古地应力；
获取所述现今构造以及所述各个演化阶段构造分别对应构造特征的物理力学参数；
根据所述现今构造以及所述各个演化阶段构造分别对应构造特征的物理力学参数以及所述现今地应力和所述各个古地应力，将所述现今地质模型和所述各个古构造地质模型分别转换为现今构造力学模型和各个古构造力学模型；
根据所述现今构造力学模型和所述各个古构造力学模型，预测所述待研究区域内所述现今构造的地层裂缝发育区带。
2.根据权利要求1所述的方法，其特征在于，在所述建立待研究区域内现今构造的现今地质模型和所述现今构造所经历的各个演化阶段构造的各个古构造地质模型之前，还包括：
解析待研究区域，确定所述待研究区域的构造特征以及所述现今构造所经历的多个演化阶段；
根据所述待研究区域的构造特征和所述现今构造所经历的多个演化阶段，确定影响所述待研究区域内地应力分布的地层界面和断层。
3.根据权利要求2所述的方法，其特征在于，所述获取所述现今构造的现今地应力和所述各个演化阶段构造的各个古地应力，包括：
分别获取所述现今构造的现今地层岩芯和所述各个演化阶段构造的各个古地层岩芯；
利用地层岩芯的古地磁定向和岩石声发射特性，分别测定所述现今构造的现今地应力和所述各个演化阶段构造的各个古地应力。
4.根据权利要求2所述的方法，其特征在于，所述获取所述现今构造以及所述各个演化阶段构造分别对应构造特征的物理力学参数，包括：
分别获取所述现今构造和所述各个演化阶段构造在所述地层界面的地层岩芯以及所述现今构造和所述各个演化阶段构造对应断层的断层介质；
利用岩石力学测试方法，获取所述现今构造和所述各个演化阶段构造分别对应的各个地层界面的地层岩芯的物理力学参数以及所述现今构造和所述各个演化阶段构造对应断层的断层介质的物理力学参数。
5.根据权利要求1～4任一项所述的方法，其特征在于，所述根据所述现今构造以及所述各个演化阶段构造分别对应构造特征的物理力学参数以及所述现今地应力和所述各个古地应力，将所述现今地质模型和所述各个古构造地质模型分别转换为现今构造力学模型和各个古构造力学模型，包括：
利用所述现今构造对应构造特征的物理力学参数将所述现今地质模型离散为所述现今构造力学模型，所述现今地应力作为所述现今构造力学模型的边界条件；
利用各个演化阶段构造分别对应构造特征的物理力学参数将所述各个古构造地质模型离散为所述各个古构造力学模型，所述各个古地应力分别作为所述各个古构造力学模型的边界条件。
6.一种地层裂缝发育区带预测装置，其特征在于，包括：
模型建立模块，用于建立待研究区域内现今构造的现今地质模型和所述现今构造所经历的各个演化阶段构造的各个古构造地质模型；
地应力获取模块，用于获取所述现今构造的现今地应力和所述各个演化阶段构造的各个古地应力；
力学参数获取模块，用于获取所述现今构造以及所述各个演化阶段构造分别对应构造特征的物理力学参数；
模型转换模块，用于根据所述现今构造以及所述各个演化阶段构造分别对应构造特征的物理力学参数以及所述现今地应力和所述各个古地应力，将所述现今地质模型和所述各个古构造地质模型分别转换为现今构造力学模型和各个古构造力学模型；
裂缝预测模块，用于根据所述现今构造力学模型和所述各个古构造力学模型，预测所述待研究区域内所述现今构造的地层裂缝发育区带。
7.根据权利要求6所述的装置，其特征在于，还包括：解析模块和确定模块；
所述解析模块，用于在所述模型建立模块建立待研究区域内现今构造的现今地质模型和所述现今构造所经历的各个演化阶段构造的各个古构造地质模型之前，解析待研究区域，确定所述待研究区域的构造特征以及所述现今构造所经历的多个演化阶段；
所述确定模块，用于根据所述待研究区域的构造特征和所述现今构造所经历的多个演化阶段，确定影响所述待研究区域内地应力分布的地层界面和断层。
8.根据权利要求7所述的装置，其特征在于，所述地应力获取模块，具体用于分别获取所述现今构造的现今地层岩芯和所述各个演化阶段构造的各个古地层岩芯，利用地层岩芯的古地磁定向和岩石声发射特性，分别测定所述现今构造的现今地应力和所述各个演化阶段构造的各个古地应力。
9.根据权利要求7所述的装置，其特征在于，所述力学参数获取模块，具体用于分别获取所述现今构造和所述各个演化阶段构造在所述地层界面的地层岩芯以及所述现今构造和所述各个演化阶段构造对应断层的断层介质，利用岩石力学测试方法，获取所述现今构造和所述各个演化阶段构造分别对应的各个地层界面的地层岩芯的物理力学参数以及所述现今构造和所述各个演化阶段构造对应断层的断层介质的物理力学参数。
10.根据权利要求6～9任一项所述的装置，其特征在于，所述模型转换模块，具体用于利用所述现今构造对应构造特征的物理力学参数将所述现今地质模型离散为所述现今构造力学模型，所述现今地应力作为所述现今构造力学模型的边界条件，利用各个演化阶段构造分别对应构造特征的物理力学参数将所述各个古构造地质模型离散为所述各个古构造力学模型，所述各个古地应力分别作为所述各个古构造力学模型的边界条件。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611190644.7A CN106844858B (zh) | 2016-12-21 | 2016-12-21 | 地层裂缝发育区带预测方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611190644.7A CN106844858B (zh) | 2016-12-21 | 2016-12-21 | 地层裂缝发育区带预测方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106844858A true CN106844858A (zh) | 2017-06-13 |
| CN106844858B CN106844858B (zh) | 2020-02-18 |
Family
ID=59135107
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611190644.7A Active CN106844858B (zh) | 2016-12-21 | 2016-12-21 | 地层裂缝发育区带预测方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106844858B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111060416A (zh) * | 2019-12-23 | 2020-04-24 | 中国石油大学(北京) | 一种地表露头与岩心岩石力学层划分方法及系统 |
| CN114153002A (zh) * | 2020-09-04 | 2022-03-08 | 中国石油化工股份有限公司 | 储层天然裂缝三维地质建模方法、装置、电子设备及介质 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070073542A1 (en) * | 2005-09-23 | 2007-03-29 | International Business Machines Corporation | Method and system for configurable allocation of sound segments for use in concatenative text-to-speech voice synthesis |
| US20130275101A1 (en) * | 2007-08-24 | 2013-10-17 | Sheng-Yuan Hsu | Method For Modeling Deformation In Subsurface Strata |
| CN104200039A (zh) * | 2014-09-17 | 2014-12-10 | 中国石油大学(华东) | 一种构造裂缝产状定量预测方法 |
| CN105319603A (zh) * | 2015-11-06 | 2016-02-10 | 中国石油大学(华东) | 致密砂岩储层复杂网状裂缝的预测方法 |
-
2016
- 2016-12-21 CN CN201611190644.7A patent/CN106844858B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070073542A1 (en) * | 2005-09-23 | 2007-03-29 | International Business Machines Corporation | Method and system for configurable allocation of sound segments for use in concatenative text-to-speech voice synthesis |
| US20130275101A1 (en) * | 2007-08-24 | 2013-10-17 | Sheng-Yuan Hsu | Method For Modeling Deformation In Subsurface Strata |
| CN104200039A (zh) * | 2014-09-17 | 2014-12-10 | 中国石油大学(华东) | 一种构造裂缝产状定量预测方法 |
| CN105319603A (zh) * | 2015-11-06 | 2016-02-10 | 中国石油大学(华东) | 致密砂岩储层复杂网状裂缝的预测方法 |
Non-Patent Citations (2)
| Title |
|---|
| 丁文龙等: "页岩储层构造应力场模拟与裂缝分布预测方法及应用", 《地学前缘》 * |
| 唐湘蓉，李晶: "构造应力场有限元数值模拟在裂缝预测中的应用", 《特种油气藏》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111060416A (zh) * | 2019-12-23 | 2020-04-24 | 中国石油大学(北京) | 一种地表露头与岩心岩石力学层划分方法及系统 |
| CN114153002A (zh) * | 2020-09-04 | 2022-03-08 | 中国石油化工股份有限公司 | 储层天然裂缝三维地质建模方法、装置、电子设备及介质 |
| CN114153002B (zh) * | 2020-09-04 | 2024-05-07 | 中国石油化工股份有限公司 | 储层天然裂缝三维地质建模方法、装置、电子设备及介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106844858B (zh) | 2020-02-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103258091B (zh) | 非常规油气藏水平井段三维岩体力学模型建立的方法及装置 | |
| US8768672B2 (en) | Method for predicting time-lapse seismic timeshifts by computer simulation | |
| CN103256046B (zh) | 非常规油气藏水平井全缝长压裂参数模拟的方法及装置 | |
| US10712472B2 (en) | Method and system for forming and using a subsurface model in hydrocarbon operations | |
| Angus et al. | Integrated hydro-mechanical and seismic modelling of the Valhall reservoir: a case study of predicting subsidence, AVOA and microseismicity | |
| MX2014009200A (es) | Pronostico de produccion de recursos. | |
| CN102656480A (zh) | 建模储油层盆地的方法 | |
| CN102636812B (zh) | 一种获得碳酸盐岩储层储集空间体积的方法 | |
| CN105527653B (zh) | 一种基于地质信息的虚拟井构建方法 | |
| CN101231346A (zh) | 应用地震波速度预测煤、岩体物理力学参数方法 | |
| CN106368691A (zh) | 基于岩石物理地震信息三维异常孔隙压力预测方法 | |
| Elenius et al. | Assessment of CO2 storage capacity based on sparse data: Skade Formation | |
| CN105629308A (zh) | 相控非均质力学参数地应力方法 | |
| CN110501744A (zh) | 烃源岩有机碳地球物理定量预测方法、装置、设备及存储介质 | |
| CN109446735A (zh) | 一种模拟测井数据的生成方法、设备以及系统 | |
| CN105242307A (zh) | 复杂碳酸盐储层地震孔隙度获取方法及装置 | |
| CN105866855B (zh) | 一种地质构造演化和变形过程的分析方法 | |
| CN107728205A (zh) | 一种地层压力预测方法 | |
| Bigi et al. | Discrete fracture network of the Latemar carbonate platform | |
| CN106844858A (zh) | 地层裂缝发育区带预测方法及装置 | |
| CN107817524A (zh) | 三维地震层析成像的方法和装置 | |
| CN111008472A (zh) | 一种基于离散元的裂谷盆地伸展过程模拟方法 | |
| CN106568918A (zh) | 页岩有机碳含量toc预测方法 | |
| Felletti | Statistical modelling and validation of correlation in turbidites: an example from the Tertiary Piedmont Basin (Castagnola Fm., Northern Italy) | |
| CN104345337B (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 |