CN107423475A - 砂岩百分含量约束的岩相建模方法 - Google Patents
砂岩百分含量约束的岩相建模方法 Download PDFInfo
- Publication number
- CN107423475A CN107423475A CN201710315195.2A CN201710315195A CN107423475A CN 107423475 A CN107423475 A CN 107423475A CN 201710315195 A CN201710315195 A CN 201710315195A CN 107423475 A CN107423475 A CN 107423475A
- Authority
- CN
- China
- Prior art keywords
- percent content
- petrofacies
- constraint
- sandstone percent
- sandstone
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种砂岩百分含量约束的岩相建模方法,该砂岩百分含量约束的岩相建模方法包括:步骤1,对研究区的实钻井进行统计分析,将沉积相图转化为研究区砂岩百分含量分带图;步骤2,将位图格式的砂岩百分含量分带图进行矢量化;步骤3,建立模型初步格架;步骤4,在岩相建模的时候进行约束;步骤5,在建模的同时启用平面约束机制;步骤6,完成岩相建模过程,输出各期次岩相展布特征图。该砂岩百分含量约束的岩相建模方法基于实钻井数据统计,获得工区砂地比数据,工区实钻井数量越多获得的数据越详细,对于岩性模型的约束更加准确。
Description
技术领域
本发明涉及岩性建模领域,特别是涉及到一种砂岩百分含量约束的岩相建模方法。
背景技术
油气地质模型研究是一项用于油气描述和油气分布预测的复合学科理论和方法体系。它是集沉积学、储层地质学、构造地质学和石油地质学等地质理论,数学地质、地质统计学和油层物理学等方法为一体的,最大限度应用计算机技术进行模型预测,揭示岩性展布规律,表征地下地质特征和各种油藏参数三维空间分布,建立能描述岩性展布特征的、地质的、岩石物理的、成岩的、构造的、流体及工程等意义的参数地质模型。
在实际的建模过程中,要建立一个合理的、科学的、完整的地质模型必须根据具体的地质情况,在一定的建模策略指导下,通过一定的建模数学工具,采用合理的建模方法,建立符合或者尽可能逼近符合地下实际的地质模型,而最终建立的地质模型是否可靠,则要依靠模型的前期约束方法和后期的优选。
在进行建模的时候在算法上常采用变差函数方法进行约束,同时采用相关性高的属性体进行约束,能在一定程度上增加预测的准确性,但是与实钻井对比时仍然存在不小误差,需要更加准确的约束方法。为此我们发明了一种新的砂岩百分含量约束的岩相建模方法,解决了以上技术问题。
发明内容
本发明的目的是提供一种基于实钻井数据统计,获得工区砂地比数据,工区实钻井数量越多获得的数据越详细,对于岩性模型的约束更加准确的砂岩百分含量约束的岩相建模方法。
本发明的目的可通过如下技术措施来实现:砂岩百分含量约束的岩相建模方法,该砂岩百分含量约束的岩相建模方法包括:步骤1,对研究区的实钻井进行统计分析,将沉积相图转化为研究区砂岩百分含量分带图;步骤2,将位图格式的砂岩百分含量分带图进行矢量化;步骤3,建立模型初步格架;步骤4,在岩相建模的时候进行约束;步骤5,在建模的同时启用平面约束机制;步骤6,完成岩相建模过程,输出各期次岩相展布特征图。
本发明的目的还可通过如下技术措施来实现:
在步骤1,对研究区的实钻井进行岩性统计分析,获得目的层砂地比数据,结合研究区已有的沉积相平面图,将不同的沉积相带范围赋值为不同的砂岩百分含量数值,将沉积相图转化为研究区砂岩百分含量分带图。
在步骤2,将位图格式的砂岩百分含量分带图进行矢量化,并将图件的四个角位置赋值上坐标信息,获得具有坐标信息的砂岩百分含量分带图的矢量图格式。
在步骤3,将目的层顶底面以及目的层解释层位数据、井数据、曲线数据、时深关系数据作为输入,建立模型空间格架。
在步骤4,选取合适的精度将模型进行网格,选取合适的变差函数参数以及敏感属性体作为约束条件对模型进行算法上、立体上的约束。
在步骤5,在平面上采用已处理好的砂岩百分含量分带图进行约束,使用已钻井资料,对模型进行最大程度的约束,使结果更加符合地质规律,与实钻井更加吻合。
本发明中的砂岩百分含量约束的岩相建模方法,涉及岩性建模领域,特别是模型约束方法方面,利用实钻井统计获得的砂地比数据对岩性模型进行约束,使得模型更加准确且复合实际。首先将工区内的井尽可能多的对各期次的砂地比分别进行统计分析,也可辅助使用该区较成熟的沉积相图(根据统计将不同的沉积相带赋予不同的砂岩百分含量数值)获得研究区不同区域的砂岩百分含量分布图;将精度较高的砂岩百分含量的位图转化为矢量图,并且图件的四个角位置赋值实际坐标,供后期模型约束使用;选取典型井、地震层位、岩相曲线等输入数据,建立岩相模型格架,利用变差函数、敏感属性体对模型进行约束,重点是针对不同的三角洲期次采用对应的砂岩百分含量图进行约束,经约束过的岩相模型更加符合地质规律,与钻井实际吻合的更好。该方法基于实钻井数据统计,获得工区砂地比数据,工区实钻井数量越多获得的数据越详细,对于岩性模型的约束更加准确。
附图说明
图1为本发明的砂岩百分含量约束的岩相建模方法的一具体实施例的流程图;
图2为本发明的一具体实施例中研究区沉积相图坐标赋值的示意图;
图3为本发明的一具体实施例中研究区沉积相图数字化的示意图;
图4为本发明的一具体实施例中岩相建模结果分析的示意图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合附图所示,作详细说明如下。
如图1所示,图1为本发明的砂岩百分含量约束的岩相建模方法的流程图。
步骤101,对研究区的实钻井尽可能多的进行岩性统计分析,获得目的层砂地比数据,可结合研究区已有的较为成熟的沉积相平面图,将不同的沉积相带(可细分至亚相、微相)范围赋值为不同的砂岩百分含量数值,这样沉积相图就转化为研究区砂岩百分含量分带图。
步骤102,将位图格式的砂岩百分含量分带图进行矢量化,并将图件的四个角位置赋值上坐标信息,这样就获得了具有坐标信息的砂岩百分含量分带图的矢量图格式。
步骤103,将目的层顶底面以及目的层解释层位数据、井数据、曲线数据、时深关系数据等作为输入,建立好模型空间格架。
步骤104,选取合适的精度将模型进行网格,选取合适的变差函数参数以及敏感属性体作为约束条件对模型进行算法上、立体上的约束。
步骤105,在平面上则采用已处理好的砂岩百分含量分带图进行约束,可最大限度的使用已钻井资料,对模型进行最大程度的约束,使结果更加符合地质规律,与实钻井更加吻合。
步骤106,完成岩相建模过程,输出各期次岩相展布特征图。
在应用本发明的一具体实施例中,东营凹陷东部沙三段为三角洲沉积类型,作为该区重点层系具有较大的勘探潜力。发育9个期次三角洲,本研究以第二期次为例来进行说明。
第一步:对研究区的实钻井进行统计分析,结合沉积相图,将研究区不同的沉积相带赋值对应的砂地比数值,三角洲平原相砂地比赋值50%,三角洲前缘砂地比赋值70%,湖泊相砂地比赋值0,浊积岩砂地比赋值80%,分流间湾相砂地比赋值20%,最终将沉积相图转化为砂地比图;
第二步:将位图格式的砂岩百分含量分带图进行矢量化,并将图件的四个角位置赋值上坐标信息(图2、图3);
第三步:将解释层位、典型井、测井曲线等作为输入数据,建立模型初步格架;
第四步:选取合适的精度将模型进行网格,选取合适的变差函数参数以及敏感属性体作为约束条件之一,在岩相建模的时候进行约束;
第五步:在建模的同时启用平面约束机制,针对目的层第二期次,在建模的同时采用砂岩百分含量平面分布图(砂地比)对岩相模型进行约束;
第六步:完成岩相建模过程,输出各期次岩相展布特征图。
建模方法常用于开发上,也仅仅是采用常规的建模方法,约束因素不够导致了精度有所欠缺。而基于砂地比约束的岩相建模方法,充分利用了研究区钻井较多的特点,明确了不同地区砂岩百分含量分布特征,依次对岩相模型过程进行约束,使结果更加符合地质规律(图4)。通过与沉积相图(图2)进行对比,两者吻合程度也更高。
Claims (6)
1.砂岩百分含量约束的岩相建模方法,其特征在于,该砂岩百分含量约束的岩相建模方法包括:
步骤1,对研究区的实钻井进行统计分析,将沉积相图转化为研究区砂岩百分含量分带图;
步骤2,将位图格式的砂岩百分含量分带图进行矢量化;
步骤3,建立模型初步格架;
步骤4,在岩相建模的时候进行约束;
步骤5,在建模的同时启用平面约束机制;
步骤6,完成岩相建模过程,输出各期次岩相展布特征图。
2.根据权利要求1所述的砂岩百分含量约束的岩相建模方法,其特征在于,在步骤1,对研究区的实钻井进行岩性统计分析,获得目的层砂地比数据,结合研究区已有的沉积相平面图,将不同的沉积相带范围赋值为不同的砂岩百分含量数值,将沉积相图转化为研究区砂岩百分含量分带图。
3.根据权利要求1所述的砂岩百分含量约束的岩相建模方法,其特征在于,在步骤2,将位图格式的砂岩百分含量分带图进行矢量化,并将图件的四个角位置赋值上坐标信息,获得具有坐标信息的砂岩百分含量分带图的矢量图格式。
4.根据权利要求1所述的砂岩百分含量约束的岩相建模方法,其特征在于,在步骤3,将目的层顶底面以及目的层解释层位数据、井数据、曲线数据、时深关系数据作为输入,建立模型空间格架。
5.根据权利要求1所述的砂岩百分含量约束的岩相建模方法,其特征在于,在步骤4,选取合适的精度将模型进行网格,选取合适的变差函数参数以及敏感属性体作为约束条件对模型进行算法上、立体上的约束。
6.根据权利要求1所述的砂岩百分含量约束的岩相建模方法,其特征在于,在步骤5,在平面上采用已处理好的砂岩百分含量分带图进行约束,使用已钻井资料,对模型进行最大程度的约束,使结果更加符合地质规律,与实钻井更加吻合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710315195.2A CN107423475A (zh) | 2017-05-05 | 2017-05-05 | 砂岩百分含量约束的岩相建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710315195.2A CN107423475A (zh) | 2017-05-05 | 2017-05-05 | 砂岩百分含量约束的岩相建模方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107423475A true CN107423475A (zh) | 2017-12-01 |
Family
ID=60425778
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710315195.2A Pending CN107423475A (zh) | 2017-05-05 | 2017-05-05 | 砂岩百分含量约束的岩相建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107423475A (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2761953A1 (fr) * | 2010-12-15 | 2012-06-15 | IFP Energies Nouvelles | Procede d'exploitation d'un gisement petrolier a partir d'un modele d'ecoulement cale au moyen d'une mise a l'echelle de cartes de pression et de saturation |
CN104657523A (zh) * | 2013-11-20 | 2015-05-27 | 中国石油化工股份有限公司 | 砂砾岩综合地质建模方法 |
CN105136628A (zh) * | 2015-07-31 | 2015-12-09 | 中国石油天然气股份有限公司 | 三角洲沉积遥感探测方法及装置 |
CN105701319A (zh) * | 2016-03-08 | 2016-06-22 | 中国石油大学(华东) | 一种水平井规则开发井网下的沉积微相建模方法 |
-
2017
- 2017-05-05 CN CN201710315195.2A patent/CN107423475A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2761953A1 (fr) * | 2010-12-15 | 2012-06-15 | IFP Energies Nouvelles | Procede d'exploitation d'un gisement petrolier a partir d'un modele d'ecoulement cale au moyen d'une mise a l'echelle de cartes de pression et de saturation |
CN104657523A (zh) * | 2013-11-20 | 2015-05-27 | 中国石油化工股份有限公司 | 砂砾岩综合地质建模方法 |
CN105136628A (zh) * | 2015-07-31 | 2015-12-09 | 中国石油天然气股份有限公司 | 三角洲沉积遥感探测方法及装置 |
CN105701319A (zh) * | 2016-03-08 | 2016-06-22 | 中国石油大学(华东) | 一种水平井规则开发井网下的沉积微相建模方法 |
Non-Patent Citations (1)
Title |
---|
赵鹏大 等: "《临清坳陷东部油气地质异常研究与资源综合评价》", 31 March 2010, 中国地质大学出版社 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7113869B2 (en) | Method for prospect identification in asset evaluation | |
US8694261B1 (en) | 3D-well log invention | |
CN108931814A (zh) | 一种基于多属性融合的基岩裂缝预测的方法 | |
Kupfersberger et al. | Methodology for integrating analog geologic data in 3-D variogram modeling | |
CN109425896A (zh) | 白云岩油气储层分布预测方法及装置 | |
NO326598B1 (no) | Tre-dimensjonal geologisk modellering | |
US20120296618A1 (en) | Multiscale Geologic Modeling of a Clastic Meander Belt Including Asymmetry Using Multi-Point Statistics | |
CN105116449B (zh) | 一种弱反射储层的识别方法 | |
CN104834003B (zh) | 一种非常规致密储层的相控压缩系数地震预测方法 | |
CN107807410A (zh) | 储层地球物理敏感参数优选方法及系统 | |
CN105301644A (zh) | 基于多参数梯度向量和海色矩阵的油气检测方法及装置 | |
CN106405639A (zh) | 一种叠前地震储层岩性参数的反演方法 | |
CN108629459A (zh) | 储层含烃孔隙的检测方法及装置 | |
CN107942405A (zh) | 预测砂泥岩薄互层砂体累积厚度的方法 | |
CN114609675A (zh) | 基于高频旋回对碳酸盐岩地层沉积微地貌的定量恢复方法 | |
CN105986819A (zh) | 用于测井资料自动处理与综合解释的方法和装置 | |
CN112505754B (zh) | 基于高精度层序格架模型的井震协同划分沉积微相的方法 | |
CN104820237B (zh) | 计算地层的方差体的方法 | |
Caf et al. | Dolomitization geometry and reservoir quality from supervised Bayesian classification and probabilistic neural networks: Midland Basin Leonardian Wichita and Clear Fork Formations | |
Pranter et al. | Sequence-stratigraphic, petrophysical, and multicomponent seismic analysis of a shelf-margin reservoir: San Andres Formation (Permian), Vacuum field, New Mexico, United States | |
CN109283577B (zh) | 一种地震层位标定方法 | |
Mehdipour et al. | The Best Scenario for Geostatistical Modeling of Porosity in the Sarvak Reservoir in an Iranian Oil Field, Using Electrofacies, Seismic Facies, and Seismic Attributes | |
Chowdhury et al. | Geometric shapes derived from airborne gravity gradiometry data: New tools for the explorationist | |
CN105242317B (zh) | 一种纵波速度的确定方法及装置 | |
CN107423475A (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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20171201 |
|
WD01 | Invention patent application deemed withdrawn after publication |