CN108074282A - 基于三维地质建模软件的物性参数的厚度校正方法 - Google Patents
基于三维地质建模软件的物性参数的厚度校正方法 Download PDFInfo
- Publication number
- CN108074282A CN108074282A CN201610981900.8A CN201610981900A CN108074282A CN 108074282 A CN108074282 A CN 108074282A CN 201610981900 A CN201610981900 A CN 201610981900A CN 108074282 A CN108074282 A CN 108074282A
- Authority
- CN
- China
- Prior art keywords
- thickness
- data
- modeling software
- reservoir
- dimensional geological
- 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
Classifications
-
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Software Systems (AREA)
- Remote Sensing (AREA)
- Computer Graphics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Geophysics And Detection Of Objects (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种基于三维地质建模软件的物性参数的厚度校正方法,包括:步骤1,将储层厚度数据统一格式,并进行整数编码;步骤2,将厚度属性编码正确加载到三维地质建模软件中;步骤3,利用厚度属性编码计算得到准确的储层厚度属性曲线数据;步骤4,利用建立的厚度校正的解释模型,对孔隙度、渗透率等物性参数进行厚度校正计算。该基于三维地质建模软件的物性参数的厚度校正方法,通过计算和加载,准确地构建了厚度属性曲线,解决了三维地质建模中缺少储层厚度属性曲线的问题,在三维地质建模软件中实现了储层厚度的校正,计算得到了更为准确的孔隙度、渗透率等储层属性数据,并为利用厚度属性曲线进行的其他计算提供了基础数据。
Description
技术领域
本发明属于石油勘探开发的三维地质建模领域,特别是涉及到一种基于三维地质建模软件的物性参数的厚度校正方法。
背景技术
三维地质建模是对地质体的数字化表述过程,主要应用于开发阶段的油藏静态描述,借助于地质模型可以更好地进行油田储层评价、开发方案编制与调整及开发效果预测等工作。属性参数建模是三维地质建模的重要组成部分,主要包括反映地层岩性、油气资源分布等特性的多种参数,如孔隙度、渗透率、含油饱和度等。这些属性参数均由测井曲线计算得到,而测井曲线受纵向分辨率的影响,在薄层的测井响应中包含了较多围岩的信号。因此,薄层的属性参数结果应该进行储层厚度的影响校正,储层的厚度属性是进行该项校正的基础。
Petrel软件由斯伦贝谢公司开发,是目前比较成熟的三维地质建模软件之一,利用Petrel软件可以建立起定量的三维可视化油藏模型。但在Petrel软件中并不能直接获得储层厚度属性,给薄层的属性参数校正带来不便。为此,我们发明了一种基于三维地质建模软件的物性参数的厚度校正方法,解决了以上技术问题。
发明内容
本发明的目的是提供一种能够方便快捷的物性参数的厚度校正方法。
本发明的目的可通过如下技术措施来实现:基于三维地质建模软件的物性参数的厚度校正方法,该基于三维地质建模软件的物性参数的厚度校正方法包括:步骤1,将储层厚度数据统一格式,并进行整数编码;步骤2,将厚度属性编码正确加载到三维地质建模软件中;步骤3,利用厚度属性编码计算得到准确的储层厚度属性曲线数据。步骤4,利用建立的厚度校正的解释模型,对孔隙度、渗透率等物性参数进行厚度校正计算。
本发明的目的还可通过如下技术措施来实现:
步骤1包括:
1)将储层厚度数据保留小数点后1位数字;
2)构造出1套离散厚度属性数据,其数值等于储层厚度值×10;
3)对离散数据的整数编码特征需要进行初始设置,根据储层的最大厚度确定储层编码的最大值Numer_Max:
Numer_Max=max(Ti)
其中,Numer_Max表示所有小层厚度属性的最大值;Ti表示第i个小层的厚度值,i=1,2,3,…,n;n为总小层数。
在步骤2中,在地质建模软件加载属性数据时,需要满足的格式为:
起始深度 属性代码
终止深度 非代码字符
准备的小层厚度属性表格式为:
井号 起始深度 终止深度 厚度属性代码
将小层数据表快速导出为建模软件可用的数据格式,具体实现过程为:
1)读取井号列数据,并作为新建文本的文件名;
2)读取同一井号的所有行,并将井号、起始深度、终止深度、厚度属性代码分别存储;
3)将起始深度和厚度属性代码写入第一行,终止深度和非代码字符写入第一行,并用空格隔开;
4)重复以上3个步骤,将所有井的厚度属性数据加载到三维地质建模软件中。
在步骤3中,从厚度属性编码计算得到准确的井点储层厚度属性曲线数据的计算公式为:
Ti=Ti/10
其中,Ti表示第i个小层的厚度值,i=1,2,3,…,n。
在步骤4中,厚度校正的解释模型的计算公式为:
PorC=K×Por
其中,Por为厚度校正前的孔隙度值,PorC为厚度校正的孔隙度值,K为校正系数,含义是岩心声波孔隙度与测井解释声波孔隙度比值,其解释模型为:
K=a*exp(-T/b)+c,
a,b,c为模型系数,由地区资料回归得到。
本发明由于采取以上技术方案,其具有以下优点:1、该方法以地质研究小层数据表成果为基础,储层厚度参数来源准确可靠。2、在具体实现时,利用Office软件的VBA工具,可以快速地将小层数据表转换为建模软件可用的数据格式,并正确加载到地质建模软件中。3、加载后,利用软件自带的计算器进行简单计算就可以得到准确的储层厚度参数,为后续的研究提供重要的储层参数。该方法方便快捷,技巧性强。利用本发明,准确地构建了厚度属性曲线,从而解决了三维地质建模中缺少储层厚度属性曲线的问题,并为利用厚度属性曲线进行的计算提供了数据。弥补了Petrel地质建模软件无法直接提供储层厚度属性的不足,为储层厚度的影响校正,进而精确求取孔隙度等储层属性提供了数据。
附图说明
图1为本发明的基于三维地质建模软件的物性参数的厚度校正方法的一具体实施例的流程图;
图2为本发明的一具体实施例中厚度数据成果图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合附图所示,作详细说明如下。
如图1所示,图1为本发明的一种基于三维地质建模软件的物性参数的厚度校正方法的流程图。
在步骤101,具体实现过程为:
1)将储层厚度数据保留小数点后1位数字;
2)构造出1套离散厚度属性数据,其数值等于储层厚度值×10;
3)对离散数据的整数编码特征需要进行初始设置,根据储层的最大厚度确定储层编码的最大值Numer_Max:
Numer_Max=max(Ti)
其中,Numer_Max表示所有小层厚度属性的最大值;Ti表示第i个小层的厚度值,i=1,2,3,…,n;n为总小层数。
比如储层最大厚度为2.8米,则离散数据的属性编码增加至28个。流程进入到步骤103。
在步骤102,在地质建模软件加载属性数据时,一般需要满足的格式为:
起始深度 属性代码
终止深度 -9999(或其他非代码字符)
准备的小层厚度属性表格式为:
井号 起始深度 终止深度 厚度属性代码
因此,利用Office软件的VBA功能将小层数据表快速导出为建模软件可用的数据格式。其具体实现过程为:
1)读取井号列数据,并作为新建文本的文件名;
2)读取同一井号的所有行,并将井号、起始深度、终止深度、厚度属性代码分别存储;
3)将起始深度和厚度属性代码写入第一行,终止深度和-9999写入第一行,并用空格隔开;
4)重复以上3个步骤,将所有井的厚度属性数据加载到三维地质建模软件中。
实例井X1的小层数据表如表1。
表1实例井X1的小层数据表
转换后生成名为X1.txt的文档,储层厚度参数属性,其格式为:
上例中,实例X1井2914.3-2917.1米油层处加载的厚度属性为28(代表厚度2.8米×10),将X1.txt加载到Petrel模型中。流程进入到步骤103。
在步骤103,从厚度属性编码计算得到准确的井点储层厚度属性曲线数据的计算公式为:
Ti=Ti/10
其中,Ti表示第i个小层的厚度值,i=1,2,3,…,n。
利用建模软件自带的计算器功能,将厚度属性离散数据/10,得到准确的储层厚度属性曲线。
实例井X1某干层起止深度为2907.3-2908.1米,厚度为0.8米,米油层处加载的厚度属性计算方法如下:8/10,计算得到结果为0.8,代表该储层厚度属性为7.3米。
步骤104,利用建立的厚度校正的解释模型,对孔隙度、渗透率等物性参数进行厚度校正计算。厚度校正的解释模型一般计算公式为:
PorC=K×Por
其中,Por为厚度校正前的孔隙度值,PorC为厚度校正的孔隙度值,K为校正系数,含义是岩心声波孔隙度与测井解释声波孔隙度比值,其解释模型一般为:
K=a*exp(-T/b)+c,
a,b,c为模型系数,由地区资料回归得到。
图2为本发明的一具体实施例中一种基于三维地质建模软件的物性参数的厚度校正方法成果示意图。研究区a=-0.10775,b=0.48648,c=0.99914。图中左起第1道为深度道,第2道为自然伽马(GR)、自然电位(SP)曲线道,第3道为声波时差(AC)曲线道,第4道为4米梯度电阻率(R4)曲线道,第5道为感应电导率(COND)曲线道,第6道为解释结论道,第7道为储层厚度(T)曲线道,第8、9道分别为厚度校正前(POR0)、厚度校正后(PORC)孔隙度曲线道。
第1层干层起止深度为2907.3-2908.1米,厚度为0.8米,第2层干层起止深度为2909.4-2910.4米,厚度为1.0米,第1层油层起止深度为2914.3-2917.1米,厚度为2.8米。在第7道储层厚度曲线与储层厚度相符,厚度属性构建及加载正确。通过对厚度校正前、后孔隙度曲线的对比,薄层的孔隙度得到了校正。
上述各实施例仅用于说明本发明,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。
Claims (5)
1.基于三维地质建模软件的物性参数的厚度校正方法,其特征在于,该基于三维地质建模软件的物性参数的厚度校正方法包括:
步骤1,将储层厚度数据统一格式,并进行整数编码;
步骤2,将厚度属性编码正确加载到三维地质建模软件中;
步骤3,利用厚度属性编码计算得到准确的储层厚度属性曲线数据;
步骤4,利用建立的厚度校正的解释模型,对孔隙度、渗透率等物性参数进行厚度校正计算。
2.根据权利要求1所述的基于三维地质建模软件的物性参数的厚度校正方法,其特征在于,步骤1包括:
1)将储层厚度数据保留小数点后1位数字;
2)构造出1套离散厚度属性数据,其数值等于储层厚度值×10;
3)对离散数据的整数编码特征需要进行初始设置,根据储层的最大厚度确定储层编码的最大值Numer_Max:
Numer_Max=max(Ti)
其中,Numer_Max表示所有小层厚度属性的最大值;Ti表示第i个小层的厚度值,
i=1,2,3,…,n;n为总小层数。
3.根据权利要求1所述的基于三维地质建模软件的物性参数的厚度校正方法,其特征在于,在步骤2中,在地质建模软件加载属性数据时,需要满足的格式为:
起始深度属性代码
终止深度非代码字符
准备的小层厚度属性表格式为:
井号起始深度终止深度厚度属性代码
将小层数据表快速导出为建模软件可用的数据格式,具体实现过程为:
1)读取井号列数据,并作为新建文本的文件名;
2)读取同一井号的所有行,并将井号、起始深度、终止深度、厚度属性代码分别存储;
3)将起始深度和厚度属性代码写入第一行,终止深度和非代码字符写入第一行,并用空格隔开;
4)重复以上3个步骤,将所有井的厚度属性数据加载到三维地质建模软件中。
4.根据权利要求1所述的基于三维地质建模软件的物性参数的厚度校正方法,其特征在于,在步骤3中,从厚度属性编码计算得到准确的井点储层厚度属性曲线数据的计算公式为:
Ti=Ti/10
其中,Ti表示第i个小层的厚度值,i=1,2,3,…,n。
5.根据权利要求1所述的基于三维地质建模软件的物性参数的厚度校正方法,其特征在于,在步骤4中,厚度校正的解释模型的计算公式为:
PorC=K×Por
其中,Por为厚度校正前的孔隙度值,PorC为厚度校正的孔隙度值,K为校正系数,含义是岩心声波孔隙度与测井解释声波孔隙度比值,其解释模型为:
K=a*exp(-T/b)+c,
a,b,c为模型系数,由地区资料回归得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610981900.8A CN108074282B (zh) | 2016-11-08 | 2016-11-08 | 基于三维地质建模软件的物性参数的厚度校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610981900.8A CN108074282B (zh) | 2016-11-08 | 2016-11-08 | 基于三维地质建模软件的物性参数的厚度校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108074282A true CN108074282A (zh) | 2018-05-25 |
CN108074282B CN108074282B (zh) | 2021-02-23 |
Family
ID=62154145
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610981900.8A Active CN108074282B (zh) | 2016-11-08 | 2016-11-08 | 基于三维地质建模软件的物性参数的厚度校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108074282B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110599594A (zh) * | 2019-07-29 | 2019-12-20 | 成都理工大学 | 一种岩石物性结构三维建模的方法 |
CN111079217A (zh) * | 2019-12-17 | 2020-04-28 | 山东省交通规划设计院 | 一种基于bim的岩土工程综合勘察信息解译方法及*** |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102313790A (zh) * | 2011-07-19 | 2012-01-11 | 北京师范大学 | 一种海底地质体二氧化碳封存潜力评估方法 |
CN102707317A (zh) * | 2010-10-27 | 2012-10-03 | 中国石油化工股份有限公司 | 一种利用地震波吸收衰减特征进行储层分析的方法 |
CN104142516A (zh) * | 2013-10-28 | 2014-11-12 | 中国石油化工股份有限公司 | 一种薄单砂层厚度预测方法 |
CN104239673A (zh) * | 2014-06-05 | 2014-12-24 | 中国石油大学 | 一种用于含油气盆地内三维热流时空场的重建方法 |
US20150066457A1 (en) * | 2013-08-27 | 2015-03-05 | Halliburton Energy Services, Inc. | Modeling Fluid Displacement in a Well System Environment |
CN104809277A (zh) * | 2015-04-14 | 2015-07-29 | 中国石油天然气股份有限公司 | 一种超低渗致密储层水平井地质建模方法 |
CN104992468A (zh) * | 2015-07-30 | 2015-10-21 | 中国石油集团川庆钻探工程有限公司 | 缝洞型碳酸盐岩油气藏三维地质建模方法 |
CN105184867A (zh) * | 2015-09-29 | 2015-12-23 | 中国矿业大学 | 基于规则格网和角点网格技术的三维地层建模方法 |
CN105443122A (zh) * | 2015-12-28 | 2016-03-30 | 中国石油天然气股份有限公司 | 一种测井解释模型的处理方法及装置 |
CN105549082A (zh) * | 2014-10-29 | 2016-05-04 | 中国石油天然气股份有限公司 | 超深碳酸盐储层三维地质力学场建立方法和*** |
CN105626054A (zh) * | 2014-11-07 | 2016-06-01 | 中国石油化工股份有限公司 | 厚层断块油藏断棱精细刻画的方法 |
CN105911592A (zh) * | 2016-06-30 | 2016-08-31 | 中国石油化工股份有限公司 | 一种实钻约束的三维地震资料长波长静校正方法 |
-
2016
- 2016-11-08 CN CN201610981900.8A patent/CN108074282B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102707317A (zh) * | 2010-10-27 | 2012-10-03 | 中国石油化工股份有限公司 | 一种利用地震波吸收衰减特征进行储层分析的方法 |
CN102313790A (zh) * | 2011-07-19 | 2012-01-11 | 北京师范大学 | 一种海底地质体二氧化碳封存潜力评估方法 |
US20150066457A1 (en) * | 2013-08-27 | 2015-03-05 | Halliburton Energy Services, Inc. | Modeling Fluid Displacement in a Well System Environment |
CN104142516A (zh) * | 2013-10-28 | 2014-11-12 | 中国石油化工股份有限公司 | 一种薄单砂层厚度预测方法 |
CN104239673A (zh) * | 2014-06-05 | 2014-12-24 | 中国石油大学 | 一种用于含油气盆地内三维热流时空场的重建方法 |
CN105549082A (zh) * | 2014-10-29 | 2016-05-04 | 中国石油天然气股份有限公司 | 超深碳酸盐储层三维地质力学场建立方法和*** |
CN105626054A (zh) * | 2014-11-07 | 2016-06-01 | 中国石油化工股份有限公司 | 厚层断块油藏断棱精细刻画的方法 |
CN104809277A (zh) * | 2015-04-14 | 2015-07-29 | 中国石油天然气股份有限公司 | 一种超低渗致密储层水平井地质建模方法 |
CN104992468A (zh) * | 2015-07-30 | 2015-10-21 | 中国石油集团川庆钻探工程有限公司 | 缝洞型碳酸盐岩油气藏三维地质建模方法 |
CN105184867A (zh) * | 2015-09-29 | 2015-12-23 | 中国矿业大学 | 基于规则格网和角点网格技术的三维地层建模方法 |
CN105443122A (zh) * | 2015-12-28 | 2016-03-30 | 中国石油天然气股份有限公司 | 一种测井解释模型的处理方法及装置 |
CN105911592A (zh) * | 2016-06-30 | 2016-08-31 | 中国石油化工股份有限公司 | 一种实钻约束的三维地震资料长波长静校正方法 |
Non-Patent Citations (1)
Title |
---|
徐浩: "喇嘛甸油田气层厚度解释标准研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110599594A (zh) * | 2019-07-29 | 2019-12-20 | 成都理工大学 | 一种岩石物性结构三维建模的方法 |
CN110599594B (zh) * | 2019-07-29 | 2021-07-20 | 成都理工大学 | 一种岩石物性结构三维建模的方法 |
CN111079217A (zh) * | 2019-12-17 | 2020-04-28 | 山东省交通规划设计院 | 一种基于bim的岩土工程综合勘察信息解译方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN108074282B (zh) | 2021-02-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bisdom et al. | Calibrating discrete fracture-network models with a carbonate three-dimensional outcrop fracture network: Implications for naturally fractured reservoir modeling | |
CA2783977C (en) | Method and system for partitioning parallel simulation models | |
CN104136942B (zh) | 用于大规模并行储层仿真的千兆单元的线性求解方法和装置 | |
CN105913494B (zh) | 多尺度裂缝精细地质建模及数值模拟方法和装置 | |
CN107367762B (zh) | 一种确定储层参数的方法及装置 | |
CN107942404B (zh) | 一种确定裂缝的方法及装置 | |
CN105223617A (zh) | 利用地震属性控制初始波阻抗的建模方法及*** | |
King et al. | Reservoir modeling: From rescue to resqml | |
US20180038974A1 (en) | System and method for petro-elastic modeling | |
Hannes et al. | Revisiting hydraulic hysteresis based on long‐term monitoring of hydraulic states in lysimeters | |
Correia et al. | Multiscale integration for karst-reservoir flow-simulation models | |
Spooner | Lifting the fog of confusion surrounding clay and shale in petrophysics | |
CN108074282A (zh) | 基于三维地质建模软件的物性参数的厚度校正方法 | |
Martin et al. | Quantifying downstream, vertical and lateral variation in fluvial deposits: implications from the Huesca Distributive Fluvial System | |
CN110863818B (zh) | 一种剩余油/气分布的描述方法及装置 | |
Ziegler et al. | Manual of the Matlab Script FAST Calibration v2. 4 | |
CN106970423A (zh) | 一种确定页岩气地层构造模型的方法及装置 | |
Mendicelli et al. | Italian seismic amplification factors for peak ground acceleration and peak ground velocity | |
CN111472763B (zh) | 地层厚度预测方法及装置 | |
Uli et al. | Integrated carbonate field management using forward stratigraphic modelling approach. Case study of a field in Central Luconia | |
CN113550735B (zh) | 一种利用测井资料计算页岩油储层含油量的方法 | |
Korus et al. | Integrating AEM and borehole data for regional hydrogeologic synthesis: Tools and examples from Nebraska, USA | |
Heidbach et al. | Manual of the Tecplot 360 Add-on GeoStress v2. 0 | |
Wharton Jr | Isopachous maps of sand reservoirs | |
BR112014032302B1 (pt) | Método para determinação de um diagrama de truncamento para uma estimativa plurigaussiana de parâmetros geológicos associados a uma zona de um subsolo real, método para determinar os parâmetros geológicos estimados em um modelo de subsolo real, meio de armazenamento não transitório legível em computador, e dispositivo para determinação de um diagrama de truncamento para uma estimativa plurigaussiana do parâmetro geológico associado a uma zona de um subsolo real |
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 |