CN109598068B - 古构造约束建模方法、装置和设备 - Google Patents
古构造约束建模方法、装置和设备 Download PDFInfo
- Publication number
- CN109598068B CN109598068B CN201811487404.2A CN201811487404A CN109598068B CN 109598068 B CN109598068 B CN 109598068B CN 201811487404 A CN201811487404 A CN 201811487404A CN 109598068 B CN109598068 B CN 109598068B
- Authority
- CN
- China
- Prior art keywords
- dimensional
- model
- target
- attribute
- 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 56
- 238000004088 simulation Methods 0.000 claims description 31
- 230000006870 function Effects 0.000 claims description 24
- 238000010276 construction Methods 0.000 claims description 18
- 230000008859 change Effects 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 3
- 238000011161 development Methods 0.000 abstract description 9
- 238000011160 research Methods 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 12
- 238000004590 computer program Methods 0.000 description 4
- 230000008021 deposition Effects 0.000 description 4
- 230000035699 permeability Effects 0.000 description 4
- 239000011148 porous material Substances 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000002948 stochastic simulation Methods 0.000 description 3
- 208000035126 Facies Diseases 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000009731 jinlong Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- 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
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Acoustics & Sound (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Geophysics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种古构造约束建模方法、装置和设备，其中，该方法包括：获取目标井区的三维地震资料、测井资料以及岩芯数据；根据所述目标井区的测井资料、所述岩芯数据以及所述三维地震资料，针对目标属性构建所述目标井区的三维属性模型；根据所述三维地震资料，确定所述目标井区的古构造特征；根据所述古构造特征构建所述目标井区的三维古构造模型；基于所述三维古构造模型对所述三维属性模型进行校正。在本发明实施例中，通过建立目标井区的三维古构造模型对三维属性模型进行校正，使得校正后的三维属性模型可以用于分析当前实际地质条件对断层发育样式的影响，从而为地质研究提供理论依据。
Description
技术领域
本发明涉及油藏建模技术领域，特别涉及一种古构造约束建模方法、装置和设备。
背景技术
火山岩具有很好的储集物性，蕴藏着丰富的油气资源，在对火山岩油藏的开发过程中需要根据相关地震数据解释、火山岩相分布数据、单井物性数据等基础地质数据对火山岩油藏进行三维属性建模，从而对火山岩油藏进行储层表征，为油藏开发提供依据。
由于火山岩油藏在成因类型、储层类型、储渗类型以及气水分布上均较为复杂，并且容易受上覆风化壳的淋滤溶蚀形成溶蚀孔隙，导致了火山岩油藏在建模表征上具有一定的难度。而现有的三维属性建模的方法在地质构造成因在理论与实际之间缺乏有效的关联手段，从而无法实现对火山岩油藏储集体上的属性特征进行定性描述。
针对上述问题，目前尚未提出有效的解决方案。
发明内容
本发明实施例提供了一种古构造约束建模方法、装置和设备，以解决现有技术中的三维属性建模方法无法实现对火山岩油藏储集体上的属性特征进行定性描述、地质结构的成因分析过程得不到验证的问题。
本发明实施例提供了一种古构造约束建模方法，包括：根据所述目标井区的测井资料、所述岩芯数据以及所述三维地震资料，针对目标属性构建所述目标井区的三维属性模型；根据所述三维地震资料，确定所述目标井区的古构造特征；根据所述古构造特征构建所述目标井区的三维古构造模型；基于所述三维古构造模型对所述三维属性模型进行校正。
在一个实施例中，根据所述古构造特征构建所述目标井区的三维古构造模型，包括：根据所述三维地震资料的精细解释，构建所述目标井区的三维构造模型；基于所述古构造特征和所述三维构造模型确定所述目标井区的三维古构造模型。
在一个实施例中，基于所述古构造特征和所述三维构造模型确定所述三维古构造模型，包括：基于所述古构造特征，对所述三维构造模型进行层拉平。
在一个实施例中，根据所述目标井区的测井资料、所述岩芯数据以及所述三维地震资料构建所述目标井区的三维属性模型，包括：根据所述目标井区的测井资料、所述岩芯数据以及所述三维地震资料构建所述目标工区的三维岩相模型；基于所述三维岩相模型通过井控的克里金模拟确定所述目标属性的所述三维属性模型。
在一个实施例中，基于所述三维古构造模型对所述三维属性模型进行校正，包括：通过所述三维古构造模型，获取古构造面的深度数据；根据所述三维属性模型的单井解释获取所述目标属性在多个方向上的变差函数；基于所述目标井区的岩相属性，确定所述古构造面的深度数据与所述单井解释中目标属性之间的相关系数；根据所述古构造面的深度数据、所述变差函数以及所述相关系数，通过同位协同克里金随机模拟法对所述三维属性模型进行校正。
在一个实施例中，所述目标井区的油藏为火山岩油藏。
在一个实施例中，在基于所述三维古构造模型对所述三维属性模型进行校正之后，所述方法还包括：根据校正后的三维属性模型进行地质勘探。
本发明实施例还提供了一种古构造约束建模装置，包括：获取模块，用于获取目标井区的三维地震资料、测井资料以及岩芯数据；
第一构建模块，用于根据所述目标井区的测井资料、所述岩芯数据以及所述三维地震资料针对目标属性构建所述目标井区的三维属性模型；
确定模块，用于根据所述三维地震资料，确定所述目标井区的古构造特征；
第二构建模块，用于根据所述古构造特征构建所述目标井区的三维古构造模型；
校正模块，用于基于所述三维古构造模型对所述三维属性模型进行校正。
在一个实施例中，所述第二构建模块，包括：第一构建单元，用于根据所述三维地震资料的精细解释，构建所述目标井区的三维构造模型；第一确定单元，用于基于所述古构造特征和所述三维构造模型确定所述目标井区的三维古构造模型。
在一个实施例中，所述第一确定单元，包括：处理单元，用于基于所述古构造特征，对所述三维构造模型进行层拉平。
在一个实施例中，所述第一构建模块，包括：第二构建单元，用于根据所述目标井区的测井资料、所述岩芯数据以及所述三维地震资料构建所述目标工区的三维岩相模型；第二确定单元，用于基于所述三维岩相模型通过井控的克里金模拟确定所述目标属性的所述三维属性模型。
在一个实施例中，所述校正模块包括：第一获取单元，用于通过所述三维古构造模型，获取古构造面的深度数据；第二获取单元，用于根据所述三维属性模型的单井解释获取所述目标属性在多个方向上的变差函数；第三确定单元，用于基于所述目标井区的岩相属性，确定所述古构造面的深度数据与所述单井解释中目标属性之间的相关系数；校正单元，用于根据所述古构造面的深度数据、所述变差函数以及所述相关系数，通过同位协同克里金随机模拟法对所述三维属性模型进行校正。
在一个实施例中，所述目标井区的油藏为火山岩油藏。
本发明实施例还提供一种古构造约束建模设备，包括处理器以及用于存储处理器可执行指令的存储器，所述处理器执行所述指令时实现所述古构造约束建模方法的步骤。
本发明实施例还提供一种计算机可读存储介质，其上存储有计算机指令，所述指令被执行时实现所述古构造约束建模方法的步骤。
本发明实施例提供了一种古构造约束建模方法，通过建立目标井区的三维古构造模型对三维属性模型进行校正，使得校正后的三维属性模型可以用于分析当前实际地质条件对断层发育样式的影响，从而为地质研究提供理论依据。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中：
图1是根据本发明实施例提供的古构造约束建模方法步骤示意图；
图2是根据本发明具体实施例提供的古构造约束建模方法的示意图；
图3是根据本发明具体实施例提供的地层模型示意图；
图4是根据本发明具体实施例提供的三维古构造模型示意图；
图5是根据本发明具体实施例提供的三维属性模型示意图；
图6是根据本发明具体实施例提供的储层古构造面示意图；
图7是根据本发明具体实施例提供的储层古构造面的三维网格模型示意图；
图8是根据本发明具体实施例提供的古构造面深度值与井筒解释孔隙度的统计关系示意图；
图9是根据本发明具体实施例提供的单井孔隙度在不同方向的变差函数示意图；
图10是根据本发明具体实施例提供的一个孔隙度模型实现示意图；
图11是根据本发明实施例提供的古构造约束建模装置示意图；
图12是根据本发明实施例提供的古构造约束建模电子设备示意图。
具体实施方式
下面将参考若干示例性实施方式来描述本发明的原理和精神。应当理解，给出这些实施方式仅仅是为了使本领域技术人员能够更好地理解进而实现本发明，而并非以任何方式限制本发明的范围。相反，提供这些实施方式是为了使本申请公开更加透彻和完整，并且能够将本公开的范围完整地传达给本领域的技术人员。
本领域的技术人员知道，本发明的实施方式可以实现为一种系统、装置设备、方法或计算机程序产品。因此，本申请公开可以具体实现为以下形式，即：完全的硬件、完全的软件(包括固件、驻留软件、微代码等)，或者硬件和软件结合的形式。
考虑到现有的三维属性建模的方法在地质构造成因在理论与实际之间缺乏有效的关联手段，无法实现对火山岩油藏储集体上的特征进行定性描述。
基于以上问题，本发明实施例提供了一种古构造约束建模方法，如图1所示，可以包括以下步骤：
步骤S101：获取目标井区的三维地震资料、测井资料以及岩芯数据。
步骤S102：根据目标井区的测井资料、岩芯数据以及三维地震资料，针对目标属性构建目标井区的三维属性模型。
考虑到储层属性的分布依赖于岩相分布，其中岩相是沉积物的沉积环境和表明沉积环境的岩性特征和生物特征的总和。因此在建立目标井区的三维属性模型前可以整合上述测井资料、岩芯数据以及三维地震资料建立三维岩相模型，在上述三维岩相模型的约束下，根据测井资料的解释中的目标属性数据，应用井控的克里金模拟的算法建立目标属性的三维属性模型。其中，上述目标属性包括但不限于：孔隙度、渗透率。
进一步的，在建立的三维岩相模型中的每一个岩相范围内，可以分别应用井控克里金模拟，产生井间属性数据，并根据每个岩相范围内的井间属性数据建立目标井区的三维属性模型。
步骤S103：根据三维地震资料，确定目标井区的古构造特征。
考虑到常规的地震解释实质上是对三维资料的二维平面解释，使用该解释方法难以发现小的构造和地层特征，而采用三维地震资料的精细解释技术可以提高地震资料解释的准确性。因此，可以对上述三维地震资料进行精细解释，以获取目标井区地层解释数据、断层解释数据、速度模型以及井点分层数据。
在区域构造的大背景下，可以将上述断层解释数据中各断层进行合理组合，建立初始三维断层模型。为了对上述初始三维断层模型进行进一步优化，可以利用上述速度模型将前期地震解释形成的时间域的三维断层模型转化为深度域，结合区域构造背景、断层的几何趋势以及单井上的断点解释结果，对深度域的三维断层模型中的断层进行几何学上的调整与简化，并将上述调整与简化后的各个断层在三维空间中进行合理组合，建立目标井区三维断层模型。
根据上述地层解释数据建立时间域的三维层面模型，将建立的时间域的三维层面模型转化到深度域，对目标井区边界范围内的三维空间进行网格化，其中，平面网格的精度可以设置为20m*20m，可以理解的是，平面网格的精度可以根据实际需要设定为其它值，本申请对此不作限定。然后，在此基础上，根据井点分层数据和深度域地震层面模型，建立目标层面的地层模型。
根据上述目标井区三维断层模型和地层模型，确定目标井区的古构造特征。
步骤S104：根据古构造特征构建目标井区的三维古构造模型。
对上述三维地震资料进行精细解释，以获取目标井区构造的断裂和地层分布特征信息。基于所述三维构造建模的基本原则，根据上述目标井区构造的断裂和地层分布特征信息，结合地质模型的相似性原则，构建目标井区的三维构造模型。其中，上述地质模型相似性原则主要是建立的三维构造模型中断层的空间形态和尺度与实际地震剖面解释得到的断层产状要素保持一致。
基于上述目标井区的三维构造模型以及古构造特征，对三维构造模型中不整合面进行层拉平，以恢复形成不整合面活动之前的三维古构造模型。其中，上述不整合面指曾经沉积区遭受区域抬升后，发生沉积间断-剥蚀，后期又沉降发生沉积的作用面。这个面与上覆地层呈不整合接触(角度不整合、假整合)。常见不整合面识别标志包括但不限于：风化面、底砾、孔洞缝发育。
步骤S105：基于三维古构造模型对三维属性模型进行校正。
为了获取古构造面的深度数据，可以将上述三维古构造模型赋值到三维网格模型中，得到的三维网格模型中垂直方向上所有网格值均相等，且等于该点平面坐标对应的古构造面的深度数据。
基于上述三维岩相模型以及测井资料，针对不同岩相分别进行单井解释，获取目标属性在不同方向上的变差函数，定量表征目标属性的三维变化规律，其中，上述变差函数是用来描述油藏属性空间变化的一种方法，可以定量的描述区域化变量的空间相关项，这种空间上的相关性是各向异性的，因此需要从不同方向上描述某个属性的变差函数。
针对上述不同岩相，分别统计上述古构造面的深度数据与单井解释中目标属性的相关关系和相关系数。并在上述三维属性模型中赋予上述古构造面的深度数据，以单井解释的目标属性为硬数据、古构造面的深度数据为协同变量，应用同位协同克里金随机模拟的方法对目标属性进行随机模拟，模拟可产生多个实现，每个实现均不同但均为等可能的。根据同位协同克里金随机模拟过程中目标参数变化与模拟结果的对应关系，对三维属性模型进行校正，根据校正后的三维属性模型可以分析当前的实际地质条件对断层发育样式的影响，从而为地质研究提供理论依据。
下面结合一个具体实施例对上述方法进行说明，然而，值得注意的是，该具体实施例仅是为了更好地说明本申请，并不构成对本申请的不当限定。
本发明实施提供了一种古构造约束建模方法，以金龙10井区火山岩油藏储层建模为例，如图2所示，可以包括：
步骤S201：基于三维地震资料精细解释确定古构造特征。
获取井区的三维地震资料、测井资料以及岩芯数据，在地震剖面上对井区主要断层进行精细解释，在区域构造的大背景下，将上述断层解释数据中各断层进行合理组合，建立初始三维断层模型。为了对上述初始三维断层模型进行进一步优化，可以利用上述速度模型将前期地震解释形成的时间域的三维断层模型转化为深度域，结合区域构造背景、断层的几何趋势以及单井上的断点解释结果，对深度域的三维断层模型中的断层进行几何学上的调整与简化，并将上述调整与简化后的各个断层在三维空间中进行合理组合，建立目标井区三维断层模型。
根据上述地层解释数据建立时间域的三维层面模型，将建立的时间域的三维层面模型转化到深度域，对目标井区边界范围内的三维空间进行网格化，其中，将平面网格的精度设置为20m*20m，然后，在此基础上，根据井点分层数据和深度域地震层面模型，建立目标层面的地层模型，上述地层模型如图3所示，其中，图中所示类似JL10的图标表示目标井区的井名。
根据上述井区三维断层模型和地层模型，确定井区的古构造特征。
步骤S202：基于构造的断裂和地层分层特征构建三维古构造模型。
对上述三维地震资料进行精细解释，以获取井区构造的断裂和地层分布特征信息。基于所述三维构造建模的基本原则，根据上述构造的断裂和地层分布特征信息，结合地质模型的相似性原则，构建井区的三维构造模型。
基于上述三维构造模型以及古构造特征，将风化壳进行层拉平，以恢复储层在遭受次生溶蚀时的顶面构造面，即三维古构造模型，如图4所示，其表示储层到风化壳的距离，古构造面越高，储层到风化壳距离则越近。其中，图中所示类似JI122的图标表示目标井区的井名。
步骤S203：基于储层属性特征建立三维属性模型。
由于储层属性的分布依赖于火成岩岩相分布，因而，需要首先整合测井资料、岩心资料、三维地震资料以及地质模式信息，确定储层的属性特征，并建立基质火成岩岩相三维模型；在基质火成岩岩相三维模型；约束下，以测井解释的孔渗数据为硬数据，应用高斯克里金随机模拟和高斯协同克里金随机模拟的算法分别建立孔隙度和渗透率的三维属性模型，参见图5，其中，图中所示类似JL1010的图标表示目标井区的井名。
步骤S204：基于三维古构造模型校正三维属性模型。
将三维古构造模型中储层古构造面(图6)赋值到三维网格模型中，参见图7，模型垂向所有网格值均相等且等于该点平面坐标对应的古构造面的深度值。其中，图6中所示古高、古低代表着古构造面中地形的高低。图7中图例所示为古构造面的深度，图中等值线代表现今构造面的深度，二者不一致。
针对不同火成岩相，分别统计古构造面深度值与井筒解释孔隙度的相关关系及相关系数，古构造面深度值与井筒解释孔隙度的统计关系如图8所示，其中，图中直线所示为古构造面深度值与井筒解释孔隙度二者的线性关系，二者的相关系数R为0.4135。该相关系数将作为后续古构造协同相控孔隙度随机模拟的一个输入数据。
在不同火成岩相内，分别统计单井孔隙度在不同方向的变差函数，参见图9，以此定量表征孔隙度的三维变化规律，为下一步古构造协同相控孔隙度随机模拟提供输入数据。其中，由于各向异性，变差函数需要从不同的方向上进行计算，通常需要从主变程、次变程以及垂直变程方向上进行计算。
将单井解释的孔隙度数据作为硬数据，以储层顶面古构造模型为协同变量，输入前述已经统计得到的古构造与孔隙度的相关系数及变差函数，应用同位协同克里金随机模拟的方法，来对井间孔隙度值进行随机模拟，模拟可产生多个实现。图10为模拟产生的一个孔隙度模型实现，根据同位协同克里金随机模拟过程中目标参数变化与模拟结果的对应关系，对三维属性模型进行校正。
基于同一发明构思，本发明实施例中还提供了一种古构造约束建模装置，如下面的实施例所述。由于古构造约束建模装置解决问题的原理与古构造约束建模方法相似，因此古构造约束建模装置的实施可以参见古构造约束建模方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。图11是本发明实施例的古构造约束建模装置的一种结构框图，如图11所示，包括：获取模块111、第一构建模块112、确定模块113、第二构建模块114以及校正模块115，下面对该结构进行说明。
获取模块111，用于获取目标井区的三维地震资料、测井资料以及岩芯数据。
第一构建模块112，用于根据所述目标井区的测井资料、所述岩芯数据以及所述三维地震资料针对目标属性构建所述目标井区的三维属性模型。
考虑到储层属性的分布依赖于岩相分布，其中岩相是沉积物的沉积环境和表明沉积环境的岩性特征和生物特征的总和。因此，在一个实施方式，上述可以包括第二构建单元，用于在建立目标井区的三维属性模型前可以整合上述测井资料、岩芯数据以及三维地震资料建立三维岩相模型，在上述三维岩相模型的约束下，在第二确定单元中根据测井资料的解释中的目标属性数据，应用井控的克里金模拟的算法建立目标属性的三维属性模型。其中，上述目标属性包括但不限于：孔隙度、渗透率。
进一步的，在建立的三维岩相模型中的每一个岩相范围内，可以分别应用井控克里金模拟，产生井间属性数据，并根据每个岩相范围内的井间属性数据建立目标井区的三维属性模型。
确定模块113，用于根据所述三维地震资料，确定所述目标井区的古构造特征。
考虑到常规的地震解释实质上是对三维资料的二维平面解释，使用该解释方法难以发现小的构造和地层特征，而采用三维地震资料的精细解释技术可以提高地震资料解释的准确性。因此，可以对上述三维地震资料进行精细解释，以获取目标井区地层解释数据、断层解释数据、速度模型以及井点分层数据。
在区域构造的大背景下，可以将上述断层解释数据中各断层进行合理组合，建立初始三维断层模型。为了对上述初始三维断层模型进行进一步优化，可以利用上述速度模型将前期地震解释形成的时间域的三维断层模型转化为深度域，结合区域构造背景、断层的几何趋势以及单井上的断点解释结果，对深度域的三维断层模型中的断层进行几何学上的调整与简化，并将上述调整与简化后的各个断层在三维空间中进行合理组合，建立目标井区三维断层模型。
根据上述地层解释数据建立时间域的三维层面模型，将建立的时间域的三维层面模型转化到深度域，对目标井区边界范围内的三维空间进行网格化，其中，平面网格的精度可以设置为20m*20m，可以理解的是，平面网格的精度可以根据实际需要设定为其它值，本申请对此不作限定。然后，在此基础上，根据井点分层数据和深度域地震层面模型，建立目标层面的地层模型。
根据上述目标井区三维断层模型和地层模型，确定目标井区的古构造特征。
第二构建模块114，用于根据所述古构造特征构建所述目标井区的三维古构造模型。
在一个实施方式中，上述第二构建模块114可以包括：第一构建单元，用于对上述三维地震资料进行精细解释，以获取目标井区构造的断裂和地层分布特征信息。基于所述三维构造建模的基本原则，根据上述目标井区构造的断裂和地层分布特征信息，结合地质模型的相似性原则，构建目标井区的三维构造模型。其中，上述地质模型相似性原则主要是建立的三维构造模型中断层的空间形态和尺度与实际地震剖面解释得到的断层产状要素保持一致。
上述第二构建模块114还可以包括：第一确定单元，用于基于上述目标井区的三维构造模型以及古构造特征，其中，第一确定单元中包括处理单元，用于对三维构造模型中不整合面进行层拉平，以恢复形成不整合面活动之前的三维古构造模型。其中，上述不整合面指曾经沉积区遭受区域抬升后，发生沉积间断-剥蚀，后期又沉降发生沉积的作用面。这个面与上覆地层呈不整合接触(角度不整合、假整合)。常见不整合面识别标志包括但不限于：风化面、底砾、孔洞缝发育。
校正模块115，用于基于所述三维古构造模型对所述三维属性模型进行校正。
具体的，该校正模块可以包括：第一获取单元，可以用于将上述三维古构造模型赋值到三维网格模型中，得到的三维网格模型中垂直方向上所有网格值均相等，且等于该点平面坐标对应的古构造面的深度数据；第二获取单元，可以用于基于上述三维岩相模型以及测井资料，针对不同岩相分别进行单井解释，获取目标属性在不同方向上的变差函数，定量表征目标属性的三维变化规律，其中，上述变差函数是用来描述油藏属性空间变化的一种方法，可以定量的描述区域化变量的空间相关项，这种空间上的相关性是各向异性的，因此需要从不同方向上描述某个属性的变差函数；第三确定单元，可以用于针对上述不同岩相，分别统计上述古构造面的深度数据与单井解释中目标属性的相关关系和相关系数。还包括校正单元，可以用于在上述三维属性模型中赋予上述古构造面的深度数据，以单井解释的目标属性为硬数据、古构造面的深度数据为协同变量，应用同位协同克里金随机模拟的方法对目标属性进行随机模拟，模拟可产生多个实现，每个实现均不同但均为等可能的。根据同位协同克里金随机模拟过程中目标参数变化与模拟结果的对应关系，对三维属性模型进行校正，根据校正后的三维属性模型可以分析当前的实际地质条件对断层发育样式的影响，从而为地质研究提供理论依据。
从以上的描述中，可以看出，本发明实施例实现了如下技术效果：通过对目标井区的三维地震资料进行精细解释获取实际的地址条件参数，确定目标井区的三维古构造模型，并通过三维古构造模型对三维属性模型进行校正，使得校正后的三维属性模型可以用于分析当前实际地质条件对断层发育样式的影响，从而为地质研究提供理论依据。
本申请实施方式还提供了一种电子设备，具体可以参阅图12所示的基于本申请实施例提供的古构造约束建模方法的电子设备组成结构示意图，电子设备具体可以包括输入设备121、处理器122、存储器123。其中，输入设备121具体可以用于输入标井区的三维地震资料、测井资料以及岩芯数据。处理器122具体可以用于根据目标井区的测井资料、岩芯数据以及三维地震资料，针对目标属性构建目标井区的三维属性模型；根据三维地震资料，确定目标井区的古构造特征；根据古构造特征构建目标井区的三维古构造模型；基于三维古构造模型对三维属性模型进行校正。存储器123具体可以用于存储目标井区的测井资料、岩芯数据、三维地震资料以及校正后的三维属性模型等参数。
在本实施方式中，所述输入设备具体可以是用户和计算机系统之间进行信息交换的主要装置之一。所述输入设备可以包括键盘、鼠标、摄像头、扫描仪、光笔、手写输入板、语音输入装置等；输入设备用于把原始数据和处理这些数的程序输入到计算机中。所述输入设备还可以获取接收其他模块、单元、设备传输过来的数据。所述处理器可以按任何适当的方式实现。例如，处理器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application Specific Integrated Circuit，ASIC)、可编程逻辑控制器和嵌入微控制器的形式等等。所述存储器具体可以是现代信息技术中用于保存信息的记忆设备。所述存储器可以包括多个层次，在数字系统中，只要能保存二进制数据的都可以是存储器；在集成电路中，一个没有实物形式的具有存储功能的电路也叫存储器，如RAM、FIFO等；在系统中，具有实物形式的存储设备也叫存储器，如内存条、TF卡等。
在本实施方式中，该电子设备具体实现的功能和效果，可以与其它实施方式对照解释，在此不再赘述。
本申请实施方式中还提供了一种基于古构造约束建模方法的计算机存储介质，所述计算机存储介质存储有计算机程序指令，在所述计算机程序指令被执行时实现：基于所述三维古构造模型对三维属性模型进行校正。
在本实施方式中，上述存储介质包括但不限于随机存取存储器(Random AccessMemory,RAM)、只读存储器(Read-Only Memory,ROM)、缓存(Cache)、硬盘(Hard DiskDrive,HDD)或者存储卡(Memory Card)。所述存储器可以用于存储计算机程序指令。网络通信单元可以是依照通信协议规定的标准设置的，用于进行网络连接通信的接口。
在本实施方式中，该计算机存储介质存储的程序指令具体实现的功能和效果，可以与其它实施方式对照解释，在此不再赘述。
显然，本领域的技术人员应该明白，上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，并且在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明实施例不限制于任何特定的硬件和软件结合。
应该理解，以上描述是为了进行图示说明而不是为了进行限制。通过阅读上述描述，在所提供的示例之外的许多实施方式和许多应用对本领域技术人员来说都将是显而易见的。因此，本申请的范围不应该参照上述描述来确定，而是应该参照前述权利要求以及这些权利要求所拥有的等价物的全部范围来确定。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (13)
1.一种古构造约束建模方法，其特征在于，包括：
获取目标井区的三维地震资料、测井资料以及岩芯数据；
根据所述目标井区的测井资料、所述岩芯数据以及所述三维地震资料，针对目标属性构建所述目标井区的三维属性模型；
根据所述三维地震资料，确定所述目标井区的古构造特征；
根据所述古构造特征构建所述目标井区的三维古构造模型；
基于所述三维古构造模型对所述三维属性模型进行校正；
其中，基于所述三维古构造模型对所述三维属性模型进行校正，包括：将所述三维古构造模型赋值到三维网格模型中，得到古构造面的深度数据；根据所述三维属性模型的单井解释获取所述目标属性在多个方向上的变差函数，所述目标属性在多个方向上的变差函数用于定量表征所述目标属性的三维变化规律；基于所述目标井区的岩相属性，确定所述古构造面的深度数据与所述单井解释中目标属性之间的相关系数；将所述单井解释的目标属性作为硬数据，以所述古构造面的深度数据为协同变量，根据所述变差函数和所述相关系数应用同位协同克里金随机模拟法对所述目标属性进行随机模拟，得到模拟结果；根据同位协同克里金随机模拟过程中目标参数变化与所述模拟结果的对应关系，对所述三维属性模型进行校正。
2.根据权利要求1所述的方法，其特征在于，根据所述古构造特征构建所述目标井区的三维古构造模型，包括：
根据所述三维地震资料的精细解释，构建所述目标井区的三维构造模型；
基于所述古构造特征和所述三维构造模型确定所述目标井区的三维古构造模型。
3.根据权利要求2所述的方法，其特征在于，基于所述古构造特征和所述三维构造模型确定所述三维古构造模型，包括：
基于所述古构造特征，对所述三维构造模型进行层拉平。
4.根据权利要求1所述的方法，其特征在于，根据所述目标井区的测井资料、所述岩芯数据以及所述三维地震资料构建所述目标井区的三维属性模型，包括：
根据所述目标井区的测井资料、所述岩芯数据以及所述三维地震资料构建所述目标工区的三维岩相模型；
基于所述三维岩相模型通过井控的克里金模拟确定所述目标属性的所述三维属性模型。
5.根据权利要求1所述的方法，其特征在于，所述目标井区的油藏为火山岩油藏。
6.根据权利要求1所述的方法，其特征在于，在基于所述三维古构造模型对所述三维属性模型进行校正之后，所述方法还包括：
根据校正后的三维属性模型进行地质勘探。
7.一种古构造约束建模装置，其特征在于，包括：
获取模块，用于获取目标井区的三维地震资料、测井资料以及岩芯数据；
第一构建模块，用于根据所述目标井区的测井资料、所述岩芯数据以及所述三维地震资料针对目标属性构建所述目标井区的三维属性模型；
确定模块，用于根据所述三维地震资料，确定所述目标井区的古构造特征；
第二构建模块，用于根据所述古构造特征构建所述目标井区的三维古构造模型；
校正模块，用于基于所述三维古构造模型对所述三维属性模型进行校正；其中，基于所述三维古构造模型对所述三维属性模型进行校正，包括：将所述三维古构造模型赋值到三维网格模型中，得到古构造面的深度数据；根据所述三维属性模型的单井解释获取所述目标属性在多个方向上的变差函数，所述目标属性在多个方向上的变差函数用于定量表征所述目标属性的三维变化规律；基于所述目标井区的岩相属性，确定所述古构造面的深度数据与所述单井解释中目标属性之间的相关系数；将所述单井解释的目标属性作为硬数据，以所述古构造面的深度数据为协同变量，根据所述变差函数和所述相关系数应用同位协同克里金随机模拟法对所述目标属性进行随机模拟，得到模拟结果；根据同位协同克里金随机模拟过程中目标参数变化与所述模拟结果的对应关系，对所述三维属性模型进行校正。
8.根据权利要求7所述的装置，其特征在于，所述第二构建模块，包括：
第一构建单元，用于根据所述三维地震资料的精细解释，构建所述目标井区的三维构造模型；
第一确定单元，用于基于所述古构造特征和所述三维构造模型确定所述目标井区的三维古构造模型。
9.根据权利要求8所述的装置，其特征在于，所述第一确定单元，包括：
处理单元，用于基于所述古构造特征，对所述三维构造模型进行层拉平。
10.根据权利要求7所述的装置，其特征在于，所述第一构建模块，包括：
第二构建单元，用于根据所述目标井区的测井资料、所述岩芯数据以及所述三维地震资料构建所述目标工区的三维岩相模型；
第二确定单元，用于基于所述三维岩相模型通过井控的克里金模拟确定所述目标属性的所述三维属性模型。
11.根据权利要求7所述的装置，其特征在于，所述目标井区的油藏为火山岩油藏。
12.一种古构造约束建模设备，包括处理器以及用于存储处理器可执行指令的存储器，所述处理器执行所述指令时实现权利要求1至6中任一项所述方法的步骤。
13.一种计算机可读存储介质，其上存储有计算机指令，所述指令被执行时实现权利要求1至6中任一项所述方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811487404.2A CN109598068B (zh) | 2018-12-06 | 2018-12-06 | 古构造约束建模方法、装置和设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811487404.2A CN109598068B (zh) | 2018-12-06 | 2018-12-06 | 古构造约束建模方法、装置和设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109598068A CN109598068A (zh) | 2019-04-09 |
| CN109598068B true CN109598068B (zh) | 2021-06-18 |
Family
ID=65962301
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811487404.2A Active CN109598068B (zh) | 2018-12-06 | 2018-12-06 | 古构造约束建模方法、装置和设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109598068B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112419493B (zh) * | 2019-08-21 | 2024-04-30 | 中国石油天然气股份有限公司 | 页岩储层三维属性模型建立方法及装置 |
| CN113933899B (zh) * | 2021-11-10 | 2023-07-28 | 中国石油大学(华东) | 基于地震属性约束指导下的砂砾岩储层建模方法和系统 |
| CN115393528B (zh) * | 2022-09-20 | 2023-04-18 | 中国地震局地球物理研究所 | 一种古地震探槽的三维建模方法及其系统 |
Family Cites Families (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101906965B (zh) * | 2010-07-14 | 2011-12-21 | 中国矿业大学(北京) | Vr地质环境下的地层识别与分析方法 |
| CN102313790A (zh) * | 2011-07-19 | 2012-01-11 | 北京师范大学 | 一种海底地质体二氧化碳封存潜力评估方法 |
| US20150153476A1 (en) * | 2012-01-12 | 2015-06-04 | Schlumberger Technology Corporation | Method for constrained history matching coupled with optimization |
| CN104657523A (zh) * | 2013-11-20 | 2015-05-27 | 中国石油化工股份有限公司 | 砂砾岩综合地质建模方法 |
| EP3146146A2 (en) * | 2014-05-07 | 2017-03-29 | King Abdullah University Of Science And Technology | Multi data reservior history matching and uncertainty quantification framework |
| CN104183018B (zh) * | 2014-08-24 | 2017-06-23 | 西南石油大学 | 用于带水碳酸盐岩气藏气水分布表征的六阶段建模方法 |
| CN104156617B (zh) * | 2014-08-24 | 2017-06-23 | 西南石油大学 | 用于多层砂岩气藏气层品质分类表征的六阶段建模方法 |
| CN104453834A (zh) * | 2014-10-31 | 2015-03-25 | 中国石油化工股份有限公司 | 一种井组注采关系优化调整方法 |
| FR3034894B1 (fr) * | 2015-04-09 | 2018-08-10 | Services Petroliers Schlumberger | Reservoir de saturation de gisement de petrole et modelisation de permeabilite |
| CN104809277B (zh) * | 2015-04-14 | 2017-09-01 | 中国石油天然气股份有限公司 | 一种超低渗致密储层水平井地质建模方法 |
| FR3038408B1 (fr) * | 2015-06-30 | 2017-08-25 | Services Petroliers Schlumrberger | Modelisation de la saturation et permeabilite de reservoir de champ petrolifere |
| CN107288633A (zh) * | 2016-03-31 | 2017-10-24 | 中国石油天然气股份有限公司 | 具有重力分异现象的凝析气藏地质储量计算方法 |
-
2018
- 2018-12-06 CN CN201811487404.2A patent/CN109598068B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN109598068A (zh) | 2019-04-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104992468B (zh) | 缝洞型碳酸盐岩油气藏三维地质建模方法 | |
| CN102147479B (zh) | 一种储层空间物性参数的建模方法 | |
| CN103077558B (zh) | 碳酸盐岩缝洞型油藏大型溶洞储集体分布模型的建模方法 | |
| CN109116428B (zh) | 一种缝洞型碳酸盐岩油藏不确定性建模方法及其装置 | |
| CN104850682B (zh) | 基于位置的多点地质统计学建模方法 | |
| EP2869096B1 (en) | Systems and methods of multi-scale meshing for geologic time modeling | |
| CN109598068B (zh) | 古构造约束建模方法、装置和设备 | |
| EP3293552B1 (en) | System and method for editing geological models by switching between volume-based models and surface-based structural models augmented with stratigraphic fiber bundles | |
| CN103116192A (zh) | 一种缝洞型碳酸盐岩油藏储集体建模方法 | |
| EP3329307A1 (en) | Assignment of systems tracts | |
| US10387583B2 (en) | Rotations from gradient directions | |
| EP2956912B1 (en) | Gridless simulation of a fluvio-deltaic environment | |
| CN104504754A (zh) | 一种多点统计建模的方法及装置 | |
| CN111815773A (zh) | 适用于机器学习算法的三维复杂地质模型标签制作方法 | |
| CN105549084A (zh) | 一种三维高精度速度建模方法和系统 | |
| CN112925017A (zh) | 一种断缝型储层地质建模方法 | |
| CN107609265B (zh) | 一种基于蚂蚁追踪的地层应力场有限元模拟方法及系统 | |
| CN106815412B (zh) | 一种构造应力场的模拟方法和装置 | |
| US12105241B2 (en) | Conditioning method and system for channel lobe deposition environment | |
| CN105301638B (zh) | 一种提取风化层底界面的方法和装置 | |
| CN111768503B (zh) | 一种基于三维地质模型的海砂资源量估算方法 | |
| CN109472867B (zh) | 一种将钻井位置信息影响范围定量化的方法 | |
| CN107368688B (zh) | 海相单砂体的确定方法和装置 | |
| CN113742962B (zh) | 一种基于水平井的页岩储层属性三维建模方法 | |
| CA3019487A1 (en) | Systems and methods for unfaulting point clouds |
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 |