CN114384584B - 裂缝建模方法及装置 - Google Patents
裂缝建模方法及装置 Download PDFInfo
- Publication number
- CN114384584B CN114384584B CN202011122331.4A CN202011122331A CN114384584B CN 114384584 B CN114384584 B CN 114384584B CN 202011122331 A CN202011122331 A CN 202011122331A CN 114384584 B CN114384584 B CN 114384584B
- Authority
- CN
- China
- Prior art keywords
- angle
- crack
- cracks
- low
- data
- 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
- 238000000034 method Methods 0.000 title claims abstract description 69
- 238000011161 development Methods 0.000 claims abstract description 54
- 238000003384 imaging method Methods 0.000 claims abstract description 36
- 238000010276 construction Methods 0.000 claims abstract description 24
- 206010017076 Fracture Diseases 0.000 claims description 82
- 208000010392 Bone Fractures Diseases 0.000 claims description 75
- 238000012545 processing Methods 0.000 claims description 17
- 238000004590 computer program Methods 0.000 claims description 13
- 239000012530 fluid Substances 0.000 claims description 13
- 238000003860 storage Methods 0.000 claims description 11
- 238000012937 correction Methods 0.000 claims description 8
- 230000009545 invasion Effects 0.000 claims description 8
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 230000006870 function Effects 0.000 description 14
- 238000010586 diagram Methods 0.000 description 12
- 238000005516 engineering process Methods 0.000 description 12
- 238000004891 communication Methods 0.000 description 7
- 238000012512 characterization method Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 6
- 238000009826 distribution Methods 0.000 description 5
- 241001415846 Procellariidae Species 0.000 description 4
- 239000000872 buffer Substances 0.000 description 4
- 239000006185 dispersion Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 208000035126 Facies Diseases 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 238000012800 visualization Methods 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
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- 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/34—Displaying seismic recordings or visualisation of seismic data or attributes
- G01V1/345—Visualisation of seismic data or attributes, e.g. in 3D cubes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- 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/6161—Seismic or acoustic, e.g. land or sea measurements
-
- 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/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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/70—Other details related to processing
- G01V2210/74—Visualisation of seismic data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Software Systems (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Geophysics (AREA)
- Acoustics & Sound (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Biophysics (AREA)
- Data Mining & Analysis (AREA)
- Medical Informatics (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Computer Hardware Design (AREA)
- Computational Linguistics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Computer Graphics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种裂缝建模方法及装置,方法包括:根据地震资料数据和成像测井数据生成高角度及斜交缝的裂缝密度体,根据常规测井数据和预先获得的构造曲率数据建立低角度缝裂缝密度体;根据目标区域的成像测井数据统计目标区域中的高角度裂缝、斜交缝、低角度裂缝,确定各产状裂缝在不同走向的发育比例;根据各产状裂缝在不同走向的发育比例,高角度缝及斜交缝的裂缝密度体以及低角度缝裂缝密度体生成高角度及斜交缝在各裂缝走向的裂缝密度体、低角度缝在各裂缝走向的裂缝密度体;根据在各裂缝走向的裂缝密度体进行裂缝建模。本发明可以体现低角度缝对油藏开发的影响,实现多走向多倾角复杂裂缝网络发育的油田的双重介质地质建模。
Description
技术领域
本发明涉及地质勘探技术,具体的讲是一种裂缝建模方法及装置。
背景技术
目前的裂缝建模技术主要是针对高角度缝及斜交缝(裂缝>30°),采用双重介质地质建模来完成裂缝表征的目标,为了适应不同走向裂缝精细建模的需要,部分学者提出采用裂缝分组的方法进行建模。但对于多种倾角裂缝同时发育的油田目前尚缺乏有效的裂缝表征的技术。
现有技术的裂缝建模方案难以表征不同倾角裂缝的发育特征,高角度缝及斜交缝对地下流体具有重要的控制作用,而低角度缝却与高角度缝和斜交缝有本质的区别,其对流体渗流无明显控制作用,但却对地下流体渗流速度具有重要的影响。如何对同时发育高角度缝、斜交缝和低角度缝的油田的裂缝问题进行有效表征,目前仍是缺乏有效的方法。
发明内容
为了至少解决现有技术中裂缝建模的一问题,本发明实施例提供了一种裂缝建模方法,方法包括:
根据地震资料数据和成像测井数据生成高角度及斜交缝的裂缝密度体;
根据常规测井数据和预先获得的构造曲率数据建立低角度缝裂缝密度体;
根据目标区域的成像测井数据统计目标区域中的高角度裂缝、斜交缝、低角度裂缝,确定各产状裂缝在不同走向的发育比例;
根据各产状裂缝在不同走向的发育比例,所述高角度缝及斜交缝的裂缝密度体以及低角度缝裂缝密度体生成高角度及斜交缝在各裂缝走向的裂缝密度体、低角度缝在各裂缝走向的裂缝密度体;
根据高角度及斜交缝在各裂缝走向的裂缝密度体、低角度缝在各裂缝走向的裂缝密度体进行裂缝建模。
本发明实施例中,所述的根据地震资料数据和成像测井数据生成高角度及斜交缝的裂缝密度体包括:
利用蚂蚁追踪算法对所述地震资料数据进行处理生成地震属性体数据;
获取目标区域的单井成像测井高角度裂缝和斜交缝裂缝发育强度曲线及预先建立的地质模型;
将所述的地震属性体数据作为所述地质模型的约束,根据所述目标区域的成像测井高角度裂缝和斜交缝裂缝强度曲线生成高角度缝及斜交缝的裂缝密度体。
本发明实施例中,所述的根据常规测井数据和预先获得的构造曲率数据建立低角度缝裂缝密度体包括:
根据所述常规测井数据,利用电阻率侵入校正差比法确定的高角度缝及斜交缝剔除三孔隙度比值法中识别的高角度缝及斜交缝确定三孔隙度比值法识别的单井低角度缝;
根据确定的三孔隙度比值法识别的单井低角度缝及构造曲率数据建立低角度缝裂缝密度体。
本发明实施例中,所述的各产状裂缝在各裂缝走向的裂缝发育比例包括:各产状裂缝在各裂缝走向的裂缝发育比例。
同时,本发明还提供一种裂缝建模装置,包括:
高角度及斜交缝的裂缝密度体生成模块,用于根据地震资料数据和成像测井数据生成高角度及斜交缝的裂缝密度体;
低角度缝裂缝密度体生成模块,用于根据常规测井数据和构造曲率数据建立低角度缝裂缝密度体;
裂缝发育比例统计模块,用于根据目标区域的成像测井数据统计目标区域中的高角度裂缝、斜交缝、低角度裂缝,确定各产状裂缝在不同走向的发育比例;
不同走向裂缝密度体生成模块,用于根据各产状裂缝在不同走向的发育比例,所述高角度缝及斜交缝的裂缝密度体以及低角度缝裂缝密度体生成高角度及斜交缝的裂缝密度体在各裂缝走向的裂缝密度体、低角度缝在各裂缝走向的裂缝密度体;
建模模块,用于根据高角度及斜交缝的裂缝密度体在各裂缝走向的裂缝密度体、低角度缝在各裂缝走向的裂缝密度体进行裂缝建模。
本发明实施例中,所述的高角度及斜交缝的裂缝密度体生成模块包括:
属性数据处理单元,用于利用蚂蚁追踪算法对所述地震资料数据进行处理生成地震属性体数据;
模型数据获取单元,用于获取目标区域的单井成像测井高角度裂缝和斜交缝裂缝发育强度曲线及预先建立的地质模型;
密度体生成单元,用于将所述的地震属性体数据作为所述地质模型的约束,根据所述目标区域的成像测井高角度裂缝和斜交缝裂缝强度曲线生成高角度缝及斜交缝的裂缝密度体。
本发明实施例中,所述的低角度缝裂缝密度体生成模块包括:
低角度缝识别单元,用于根据所述常规测井数据,利用电阻率侵入校正差比法确定的高角度缝及斜交缝剔除三孔隙度比值法中识别的高角度缝及斜交缝确定三孔隙度比值法识别的单井低角度缝;
低角度缝裂缝密度体生成单元,用于根据确定的三孔隙度比值法识别的单井低角度缝及构造曲率数据建立低角度缝裂缝密度体。
本发明实施例中,所述的各产状裂缝在各裂缝走向的裂缝发育比例包括:各产状裂缝在各裂缝走向的裂缝发育比例。
同时,本发明还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,处理器执行计算机程序时实现上述方法。
同时,本发明还提供一种计算机可读存储介质,计算机可读存储介质存储有执行上述方法的计算机程序。
本发明提供的裂缝建模方法及装置,对于具有不同走向和倾角裂缝发育的油藏,低角度缝和高角度缝及斜交缝具有显著差异,高角度缝及斜交缝对地下流体渗流方向具有明显的控制作用,而低角度缝对地下流体流动方向无明显影响。本方案基于对裂缝走向分组建模的基础上继续根据裂缝倾角进行分组,基于裂缝的走向和倾角两个参数进行分组建模,从而体现出低角度缝对渗流方向无明显控制作用的特征,并能继续在建模中体现高角度缝对渗流方向有明显影响的特性,从而既在建模中保留了高角度缝对油藏开发影响特征,又将低角度缝对油藏开发的影响体现出来,实现了多走向多倾角复杂裂缝网络发育的油田的双重介质地质建模。
为让本发明的上述和其他目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附图式,作详细说明如下。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明提供的裂缝建模方法的流程图;
图2为本发明实施例在Petrel软件中对北东高角度裂缝、北东斜交缝、北东低角度裂缝、北西高角度裂缝、北西斜交缝和北西低角度裂缝6组裂缝建模处理界面图;
图3为本发明在Petrel软件中微体现低角度缝对渗流方向控制较弱特征而设置Fisher模型中散度的参数图;
图4为本发明不同方向各裂缝集离散网络模型图;
图5为本发明提供的裂缝建模装置的框图;
图6为本发明提供的电子设备实施例的示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
现有技术中,对于高角度缝、斜交缝和低角度缝均发育的油田,目前的裂缝建模方案缺乏有效的表征技术,难以将不同倾角的裂缝的发育特征表征出来。本发明提供一种满足多走向多倾角发育时裂缝表征的技术建模方案,通过考虑裂缝走向和倾角两个因素的影响进行裂缝分组,从而实现对多走向多倾角裂缝的有效表征,为油田数值模拟及开发方案设计奠定了基础。
如图1所示,为本发明提供的一种裂缝建模方法的流程图,其包括:
步骤S101,根据地震资料数据和成像测井数据生成高角度及斜交缝的裂缝密度体;
步骤S102,根据常规测井数据和构造曲率数据建立低角度缝裂缝密度体;
步骤S103,根据目标区域的成像测井数据统计目标区域中的高角度裂缝、斜交缝、低角度裂缝,确定各产状裂缝在不同走向的发育比例;
步骤S104,根据各产状裂缝在不同走向的发育比例,所述高角度缝及斜交缝的裂缝密度体以及低角度缝裂缝密度体生成高角度及斜交缝的裂缝密度体在各裂缝走向的裂缝密度体、低角度缝在各裂缝走向的裂缝密度体;
步骤S105,根据高角度及斜交缝的裂缝密度体在各裂缝走向的裂缝密度体、低角度缝在各裂缝走向的裂缝密度体进行裂缝建模。
本发明实施例中,根据地震资料和成像测井数据生成高角度及斜交缝的裂缝密度体包括:
利用蚂蚁追踪算法对所述地震资料数据进行处理生成地震属性体数据;
获取目标区域的单井成像测井高角度裂缝和斜交缝裂缝发育强度曲线及预先建立的地质模型;
将所述的地震属性体数据作为所述地质模型的约束,根据所述目标区域的成像测井高角度裂缝和斜交缝裂缝强度曲线生成高角度缝及斜交缝的裂缝密度体。
具体的,本发明实施例中,基于蚂蚁追踪技术对地震资料体进行处理,获取大尺度高角度缝及斜交缝分布,并且,本实施例中,通过Petrel软件将处理后的地震体转换为地震属性体加入到地质模型中,以该属性体作为约束,对成像测井高角度缝和斜交缝裂缝密度曲线进行井间插值,从而生成高角度缝及斜交缝的裂缝密度体。本实施例中,采用的Petrel软件是以三维地质模型为中心的勘探开发一体化平台,属于地球物理专业软件。
本发明实施例中,所述的根据常规测井数据和所述构造曲率体建立低角度缝裂缝密度体包括:
根据所述测井数据,利用电阻率侵入校正差比法确定的高角度缝及斜交缝剔除三孔隙度比值法中识别的高角度缝及斜交缝确定三孔隙度比值法识别的单井低角度缝;
根据确定的三孔隙度比值法识别的单井低角度缝和构造曲率数据,建立低角度缝裂缝密度体。其中,构造曲率数据可根据构造面在Petrel中直接计算得到,本实施例中构造面可根据油田地震资料解释获得,对本领域技术人员而言可清楚获知如何从油田地震资料获得构造面及跟进构造面确定构造曲率数据,本方案中不再详细赘述。
即本发明实施例中,基于常规测井资料,使用电阻率侵入校正差比法识别的高角度缝及斜交缝(>30°)约束三孔隙度比值法,从而剔除三孔隙度比值法中识别的高角度缝及斜交缝,实现三孔隙度比值法识别单井低角度缝,并且,本实施例中具体的,使用构造曲率数据作为约束在Petrel软件中通过使用序贯高斯插值实现低角度缝裂缝密度体的建立。
本发明实施例中,所述的各产状裂缝在各裂缝走向的裂缝发育比例包括:各产状裂缝在各裂缝走向的裂缝发育比例。
本发明实施例具体通过以下技术方案实现,具体步骤包括:
(1)基于蚂蚁追踪技术对地震资料体进行处理,获取大尺度高角度缝及斜交缝分布,并通过Petrel软件将处理后的地震体转换为地震属性体加入到地质模型中,以该属性体作为约束,对成像测井进行井间插值,从而生成高角度缝及斜交缝的裂缝密度体。
(2)基于常规测井资料,使用电阻率侵入校正差比法识别的高角度缝及斜交缝(>30°)约束三孔隙度比值法,从而剔除三孔隙度比值法中识别的高角度缝及斜交缝,实现三孔隙度比值法识别单井低角度缝,并在Petrel软件中使用曲率体作为约束,通过使用序贯高斯插值实现低角度缝裂缝密度体的建立。
(3)假设研究区裂缝走向以北东-南西和北西-南东方向为主,裂缝产状有高角度缝、斜交缝和低角度缝。通过使用成像测井分别统计北东高角度裂缝、北东斜交缝、北东低角度裂缝、北西高角度裂缝、北西斜交缝和北西低角度裂缝6组裂缝发育的比例,将北东高角度裂缝、北东斜交缝、北西高角度裂缝和北西斜交缝按比例劈分高角度缝及斜交缝的裂缝密度体,从而生成不同方向裂缝发育的裂缝密度体,将北东低角度裂缝和北西低角度裂缝统计数据对低角度缝裂缝密度体进行等比例劈分,从而分别生成低角度缝不同方向的裂缝密度体。
本发明实施例中,各产状裂缝在各裂缝走向的裂缝发育比例包括:各产状裂缝在各裂缝走向的裂缝发育比例,即根据目标研究区域的裂缝的主要走向,分别同各产状裂缝在各走向的裂缝发育比例。
(4)在Petrel中根据分的六组裂缝及裂缝密度体建立裂缝模型,为了体现低角度缝对地下流体渗流方向影响较弱的特征,在裂缝分组建模的过程中将Fisher模型的散度(concentration)适当减小,从而实现低角度缝对地下流体渗流方向影响较弱的目标。
Fisher模型由一个平均极方位和一个离散度参数K来定义,离散度参数K的变化范围可以从0(无平均值)到∞(所有方位均相同)之间变化。当K→0时,Fisher分布模型趋于均匀分布,即裂缝极点均匀分布在极点图上;当K→∞时,裂缝极点的分布趋于集中分布于平均矢量极点。K的值适当减小,可以使裂缝分布更离散,符合低角度缝没有主走向的特征。
本发明对于具有多走向多产状裂缝发育油田具有适应性,尤其对建模中体现低角度裂缝对渗流方向控制较弱这一特征提出了一种新的解决方案。
本发明实施例中Petrel是一套基于Windows平台的三维可视化建模软件,它集地震解释、构造建模、岩相建模、油藏属性建模和油藏数值模拟显示及虚拟现实于一体,为地质学家、地球物理学家,岩石物理学家,油藏工程人员提供了一个共享的信息平台。同时,Petrel应用了各种先进技术:强大的构造建模技术、高精度的三维网格化技术、确定性和随机性沉积相模型建立技术、科学的岩石物理建模技术、先进的三维计算机可视化和虚拟现实技术。
图2为本发明实施例在Petrel软件中对北东高角度裂缝、北东斜交缝、北东低角度裂缝、北西高角度裂缝、北西斜交缝和北西低角度裂缝6组裂缝建模处理界面图。
图3为本发明在Petrel软件中为体现低角度缝对渗流方向控制较弱特征而设置Fisher模型中散度的参数图。
图4为本发明不同方向各裂缝集离散网络模型图。
本发明对于具有不同走向和不同倾角裂缝发育的油藏,低角度缝和高角度缝及斜交缝具有显著差异,高角度缝及斜交缝对地下流体渗流方向具有明显的控制作用,而低角度缝对地下流体流动方向无明显影响。本发明基于前人对裂缝走向分组建模的基础上继续根据裂缝倾角进行分组,基于裂缝的走向和倾角两个参数进行分组建模,并通过设置Fisher模型中散度参数,从而体现出低角度缝对渗流方向无明显控制作用的特征,并能继续在建模中体现高角度缝对渗流方向有明显影响的特性,从而既在建模中保留了高角度缝对油藏开发影响特征,又将低角度缝对油藏开发的影响体现出来,实现了多走向多倾角复杂裂缝网络发育的油田的双重介质地质建模。
本发明对于具有多走向多产状裂缝发育油田具有适应性,尤其对建模中体现低角度裂缝对渗流方向控制较弱这一特征提出了一种新的解决方案。
在本发明创造中借用了前人关于高角度缝及斜交缝裂缝建模的方案以及前人对多走向裂缝问题进行分组建模的技术路线。
本发明一是对于具有多走向多倾角裂缝发育的油田,提出考虑裂缝走向和倾角两个因素的影响,从而对裂缝进行分组,从而更加准确精细的表征了油田复杂裂缝网络。二是针对高角度缝及斜交缝对地下流体渗流方向具有控制作用而低角度缝对地下流体无明显影响的问题,提出基于分组建模的结果通过设置Fisher模型中散度参数,从而实现了对低角度缝特征的表征。
如图5所示,本发明还提供一种裂缝建模装置,包括:
高角度及斜交缝的裂缝密度体生成模块501,用于根据地震资料数据和成像测井高角度及斜交缝裂缝发育强度曲线生成高角度及斜交缝的裂缝密度体;
低角度缝裂缝密度体生成模块502,用于根据常规测井数据和构造曲率数据建立低角度缝裂缝密度体;
裂缝发育比例统计模块503,用于根据目标区域的成像测井数据统计目标区域中的高角度裂缝、斜交缝、低角度裂缝,确定各产状裂缝在各裂缝走向的裂缝发育比例;
不同走向裂缝密度体生成模块504,用于根据各产状裂缝在各裂缝走向的裂缝发育比例,所述高角度及斜交缝的裂缝密度体以及低角度缝裂缝密度体生成高角度及斜交缝的裂缝密度体在各裂缝走向的裂缝密度体、低角度缝在各裂缝走向的裂缝密度体;
建模模块505,用于根据地震资料和成像测井数据得到高角度及斜交缝的裂缝密度体在各裂缝走向的裂缝密度体、以及根据常规测井资料和构造曲率约束体建立的低角度缝密度体得到的在各裂缝走向的裂缝密度体进行裂缝建模。
本发明实施例中,所述的高角度及斜交缝的裂缝密度体生成模块包括:
属性数据处理单元,用于利用蚂蚁追踪算法对所述地震资料数据进行处理生成地震属性体数据;
模型数据获取单元,用于获取目标区域的单井成像测井高角度裂缝和斜交缝裂缝发育强度曲线及预先建立的地质模型;
密度体生成单元,用于将所述的地震属性体数据作为所述地质模型的约束,根据所述目标区域的成像测井高角度裂缝和斜交缝裂缝强度曲线生成高角度缝及斜交缝的裂缝密度体。
本发明实施例中,所述的低角度缝裂缝密度体生成模块包括:
低角度缝识别单元,用于根据所述测井数据,利用电阻率侵入校正差比法确定的高角度缝及斜交缝剔除三孔隙度比值法中识别的高角度缝及斜交缝确定三孔隙度比值法识别的单井低角度缝;
低角度缝裂缝密度体生成单元,用于根据确定的三孔隙度比值法识别的单井低角度缝和构造曲率数据建立低角度缝裂缝密度体。
本实施例还提供一种电子设备,该电子设备可以是台式计算机、平板电脑及移动终端等,本实施例不限于此。在本实施例中,该电子设备可以参照前述方法及装置的实施例,其内容被合并于此,重复之处不再赘述。
图6为本发明实施例的电子设备600的系统构成的示意框图。如图6所示,该电子设备600可以包括中央处理器100和存储器140;存储器140耦合到中央处理器100。值得注意的是,该图是示例性的;还可以使用其他类型的结构,来补充或代替该结构,以实现电信功能或其他功能。
一实施例中,裂缝建模功能可以被集成到中央处理器100中。其中,中央处理器100可以被配置为进行如下控制:
根据地震资料数据和成像测井数据生成高角度及斜交缝的裂缝密度体;
根据常规测井数据和构造曲率约束体建立低角度缝裂缝密度体;
根据目标区域的成像测井数据统计目标区域中的高角度裂缝、斜交缝、低角度裂缝,确定各产状裂缝在各裂缝走向的裂缝发育比例;
根据各产状裂缝在各裂缝走向的裂缝发育比例,所述高角度及斜交缝的裂缝密度体以及低角度缝裂缝密度体生成高角度及斜交缝的裂缝密度体在各裂缝走向的裂缝密度体、低角度缝在各裂缝走向的裂缝密度体;
根据地震资料和成像测井数据得到高角度及斜交缝的裂缝密度体在各裂缝走向的裂缝密度体、以及根据常规测井资料和构造曲率约束体建立的低角度缝密度体得到的在各裂缝走向的裂缝密度体进行裂缝建模。
如图6所示,该电子设备600还可以包括:通信模块110、输入单元120、音频处理单元130、显示器160、电源170。值得注意的是,电子设备600也并不是必须要包括图6中所示的所有部件;此外,电子设备600还可以包括图6中没有示出的部件,可以参考现有技术。
如图6所示,中央处理器100有时也称为控制器或操作控件,可以包括微处理器或其他处理器装置和/或逻辑装置,该中央处理器100接收输入并控制电子设备600的各个部件的操作。
其中,存储器140,例如可以是缓存器、闪存、硬驱、可移动介质、易失性存储器、非易失性存储器或其它合适装置中的一种或更多种。可储存上述与失败有关的信息,此外还可存储执行有关信息的程序。并且中央处理器100可执行该存储器140存储的该程序,以实现信息存储或处理等。
输入单元120向中央处理器100提供输入。该输入单元120例如为按键或触摸输入装置。电源170用于向电子设备600提供电力。显示器160用于进行图像和文字等显示对象的显示。该显示器例如可为LCD显示器,但并不限于此。
该存储器140可以是固态存储器,例如,只读存储器(ROM)、随机存取存储器(RAM)、SIM卡等。还可以是这样的存储器,其即使在断电时也保存信息,可被选择性地擦除且设有更多数据,该存储器的示例有时被称为EPROM等。存储器140还可以是某种其它类型的装置。存储器140包括缓冲存储器141(有时被称为缓冲器)。存储器140可以包括应用/功能存储部142,该应用/功能存储部142用于存储应用程序和功能程序或用于通过中央处理器100执行电子设备600的操作的流程。
存储器140还可以包括数据存储部143,该数据存储部143用于存储数据,例如联系人、数字数据、图片、声音和/或任何其他由电子设备使用的数据。存储器140的驱动程序存储部144可以包括电子设备的用于通信功能和/或用于执行电子设备的其他功能(如消息传送应用、通讯录应用等)的各种驱动程序。
通信模块110即为经由天线111发送和接收信号的发送机/接收机110。通信模块(发送机/接收机)110耦合到中央处理器100,以提供输入信号和接收输出信号,这可以和常规移动通信终端的情况相同。
基于不同的通信技术,在同一电子设备中,可以设置有多个通信模块110,如蜂窝网络模块、蓝牙模块和/或无线局域网模块等。通信模块(发送机/接收机)110还经由音频处理器130耦合到扬声器131和麦克风132,以经由扬声器131提供音频输出,并接收来自麦克风132的音频输入,从而实现通常的电信功能。音频处理器130可以包括任何合适的缓冲器、解码器、放大器等。另外,音频处理器130还耦合到中央处理器100,从而使得可以通过麦克风132能够在本机上录音,且使得可以通过扬声器131来播放本机上存储的声音。
本发明实施例还提供一种计算机可读程序,其中当在电子设备中执行所述程序时,所述程序使得计算机在所述电子设备中执行如上面实施例所述的裂缝建模方法。
本发明实施例还提供一种存储有计算机可读程序的存储介质,其中所述计算机可读程序使得计算机在电子设备中执行上面实施例所述的裂缝建模功能。
以上参照附图描述了本发明的优选实施方式。这些实施方式的许多特征和优点根据该详细的说明书是清楚的,因此所附权利要求旨在覆盖这些实施方式的落入其真实精神和范围内的所有这些特征和优点。此外,由于本领域的技术人员容易想到很多修改和改变,因此不是要将本发明的实施方式限于所例示和描述的精确结构和操作,而是可以涵盖落入其范围内的所有合适修改和等同物。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (8)
1.一种裂缝建模方法,其特征在于,所述的方法包括:
根据地震资料数据和成像测井数据生成高角度及斜交缝的裂缝密度体;
根据常规测井数据和预先获得的构造曲率数据建立低角度缝裂缝密度体;
根据目标区域的成像测井数据统计目标区域中的高角度裂缝、斜交缝、低角度裂缝,确定各产状裂缝在不同走向的发育比例;
根据各产状裂缝在不同走向的发育比例,所述高角度缝及斜交缝的裂缝密度体以及低角度缝裂缝密度体生成高角度及斜交缝在各裂缝走向的裂缝密度体、低角度缝在各裂缝走向的裂缝密度体;
根据高角度及斜交缝在各裂缝走向的裂缝密度体、低角度缝在各裂缝走向的裂缝密度体进行裂缝建模;在裂缝分组建模的过程中将Fisher模型的散度减小,从而实现低角度缝对地下流体渗流方向影响较弱的目标。
2.如权利要求1所述的裂缝建模方法,其特征在于,所述的根据地震资料数据和成像测井数据生成高角度及斜交缝的裂缝密度体包括:
利用蚂蚁追踪算法对所述地震资料数据进行处理生成地震属性体数据;
获取目标区域的单井成像测井高角度裂缝和斜交缝裂缝发育强度曲线及预先建立的地质模型;
将所述的地震属性体数据作为所述地质模型的约束,根据所述目标区域的成像测井高角度裂缝和斜交缝裂缝强度曲线生成高角度缝及斜交缝的裂缝密度体。
3.如权利要求1所述的裂缝建模方法,其特征在于,所述的根据常规测井数据和预先获得的构造曲率数据建立低角度缝裂缝密度体包括:
根据所述常规测井数据,利用电阻率侵入校正差比法确定的高角度缝及斜交缝剔除三孔隙度比值法中识别的高角度缝及斜交缝确定三孔隙度比值法识别的单井低角度缝;
根据确定的三孔隙度比值法识别的单井低角度缝及构造曲率数据建立低角度缝裂缝密度体。
4.一种裂缝建模装置,其特征在于,所述的装置包括:
高角度及斜交缝的裂缝密度体生成模块,用于根据地震资料数据和成像测井数据生成高角度及斜交缝的裂缝密度体;
低角度缝裂缝密度体生成模块,用于根据常规测井数据和预先获得的构造曲率数据建立低角度缝裂缝密度体;
裂缝发育比例统计模块,用于根据目标区域的成像测井数据统计目标区域中的高角度裂缝、斜交缝、低角度裂缝,确定各产状裂缝在不同走向的发育比例;
不同走向裂缝密度体生成模块,用于根据各产状裂缝在不同走向的发育比例,所述高角度缝及斜交缝的裂缝密度体以及低角度缝裂缝密度体生成高角度及斜交缝在各裂缝走向的裂缝密度体、低角度缝在各裂缝走向的裂缝密度体;
建模模块,用于根据高角度及斜交缝在各裂缝走向的裂缝密度体、低角度缝在各裂缝走向的裂缝密度体进行裂缝建模;在裂缝分组建模的过程中将Fisher模型的散度减小,从而实现低角度缝对地下流体渗流方向影响较弱的目标。
5.如权利要求4所述的裂缝建模装置,其特征在于,所述的高角度及斜交缝的裂缝密度体生成模块包括:
属性数据处理单元,用于利用蚂蚁追踪算法对所述地震资料数据进行处理生成地震属性体数据;
模型数据获取单元,用于获取目标区域的单井成像测井高角度裂缝和斜交缝裂缝发育强度曲线及预先建立的地质模型;
密度体生成单元,用于将所述的地震属性体数据作为所述地质模型的约束,根据所述目标区域的成像测井高角度裂缝和斜交缝裂缝强度曲线生成高角度缝及斜交缝的裂缝密度体。
6.如权利要求4所述的裂缝建模装置,其特征在于,所述的低角度缝裂缝密度体生成模块包括:
低角度缝识别单元,用于根据所述常规测井数据,利用电阻率侵入校正差比法确定的高角度缝及斜交缝剔除三孔隙度比值法中识别的高角度缝及斜交缝确定三孔隙度比值法识别的单井低角度缝;
低角度缝裂缝密度体生成单元,用于根据确定的三孔隙度比值法识别的单井低角度缝及构造曲率数据建立低角度缝裂缝密度体。
7.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至3任一项所述方法。
8.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至3任一项所述方法的计算机程序。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011122331.4A CN114384584B (zh) | 2020-10-20 | 2020-10-20 | 裂缝建模方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011122331.4A CN114384584B (zh) | 2020-10-20 | 2020-10-20 | 裂缝建模方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114384584A CN114384584A (zh) | 2022-04-22 |
CN114384584B true CN114384584B (zh) | 2024-05-28 |
Family
ID=81192560
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011122331.4A Active CN114384584B (zh) | 2020-10-20 | 2020-10-20 | 裂缝建模方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114384584B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115576007B (zh) * | 2022-11-22 | 2023-03-14 | 西南石油大学 | 一种基于乱序矩阵的半确定性裂缝建模方法及系统 |
CN115657132B (zh) * | 2022-12-29 | 2023-03-10 | 成都捷科思石油天然气技术发展有限公司 | 一种基于机器学习的储层预测方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104992067A (zh) * | 2015-07-14 | 2015-10-21 | 北京博达瑞恒科技有限公司 | 一种基于裂缝产状的裂缝分析方法 |
CN109100793A (zh) * | 2017-06-20 | 2018-12-28 | 中国石油化工股份有限公司 | 一种定量分析裂缝因素对储层影响的方法 |
CN111754623A (zh) * | 2019-03-27 | 2020-10-09 | 中国石油化工股份有限公司 | 一种基于多点地质统计学的离散裂缝建模方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3568568B1 (en) * | 2017-01-13 | 2023-10-11 | Baker Hughes Holdings LLC | Near wellbore discrete fracture networks |
-
2020
- 2020-10-20 CN CN202011122331.4A patent/CN114384584B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104992067A (zh) * | 2015-07-14 | 2015-10-21 | 北京博达瑞恒科技有限公司 | 一种基于裂缝产状的裂缝分析方法 |
CN109100793A (zh) * | 2017-06-20 | 2018-12-28 | 中国石油化工股份有限公司 | 一种定量分析裂缝因素对储层影响的方法 |
CN111754623A (zh) * | 2019-03-27 | 2020-10-09 | 中国石油化工股份有限公司 | 一种基于多点地质统计学的离散裂缝建模方法 |
Non-Patent Citations (4)
Title |
---|
Reservoir architecture and fracture characterization of low-permeability sandstone reservoir: a case study of Biandong oilfield, Jinhu depression, northern Jiangsu Basin, China;Xue Li 等;Arabian Journal of Geosciences;20180718;第1-13页 * |
中国南海流花油田生物礁储层裂缝预测建模研究;张永江 等;石油地质与工程;20170531;第31卷(第03期);第66-69页 * |
基于离散裂缝网络模型裂缝储层的随机模拟――以王窑区块为例;唐永 等;地质科技情报;20100930;第29卷(第05期);第71-76页 * |
构造裂缝识别与建模研究——以鄂尔多斯盆地中西部定边—华池地区延长组为例;潘欣;中国优秀硕士学位论文全文数据库基础科学辑;20100915(第9期);第A011-70页 * |
Also Published As
Publication number | Publication date |
---|---|
CN114384584A (zh) | 2022-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109710968B (zh) | 一种基岩潜山裂缝预测方法及装置 | |
CN114384584B (zh) | 裂缝建模方法及装置 | |
CN107765298B (zh) | 一种确定近地表速度模型的方法及装置 | |
CN104504754B (zh) | 一种油气储层多点统计建模的方法及装置 | |
CN115310381A (zh) | 基于典型单元体计算的裂隙岩体渗流场模拟方法及终端 | |
CN112415591A (zh) | 绕射波的成像方法、装置、电子设备及存储介质 | |
CN113433587B (zh) | 岩石物理反演方法及系统 | |
CN113933901B (zh) | 深度域地震数据拼接方法及装置 | |
CN113586042B (zh) | 孔隙型碳酸盐岩动态渗透率确定方法及装置 | |
CN112507551B (zh) | 一种非结构化动态网格剖分方法及装置 | |
CN114114392B (zh) | 层速度模型建立方法及装置 | |
CN117270032A (zh) | 基于三维地震的砂岩型铀矿铀特征参数反演方法及装置 | |
CN114578428A (zh) | 一种速度场高效构建方法及装置 | |
CN112462426B (zh) | 横波矢量静校正方法及装置 | |
CN112696181A (zh) | 基于地应力和天然裂缝耦合关系的压裂预测方法及装置 | |
CN111751879A (zh) | 储层孔隙度预测方法及装置 | |
CN111754611B (zh) | 基于地质模型约束的井间砂体连通图自动绘制方法及系统 | |
CN112782780A (zh) | 一种基于岩石物理相的储层评价方法、装置及设备 | |
CN104834005B (zh) | 一种确定地质界面的方法 | |
CN113075747B (zh) | 储层裂缝发育区域的预测方法及装置 | |
CN104316964B (zh) | 一种近地表表层的建模方法及系统 | |
CN114442169B (zh) | 地震资料中近源信号压制方法及装置 | |
CN112505748B (zh) | 炮域反射波拾取方法及装置 | |
CN115903026B (zh) | 复合砂体构型解析方法、设备及介质 | |
CN117214956A (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 |