CN111353204B - 一种岩心孔隙度参数场的处理方法 - Google Patents

一种岩心孔隙度参数场的处理方法 Download PDF

Info

Publication number
CN111353204B
CN111353204B CN201811562408.2A CN201811562408A CN111353204B CN 111353204 B CN111353204 B CN 111353204B CN 201811562408 A CN201811562408 A CN 201811562408A CN 111353204 B CN111353204 B CN 111353204B
Authority
CN
China
Prior art keywords
matrix
dimensional
core porosity
parameter field
processing
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
Application number
CN201811562408.2A
Other languages
English (en)
Other versions
CN111353204A (zh
Inventor
李玥洋
卢晓敏
魏星
张春
唐刚
王蓓
王娟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201811562408.2A priority Critical patent/CN111353204B/zh
Publication of CN111353204A publication Critical patent/CN111353204A/zh
Application granted granted Critical
Publication of CN111353204B publication Critical patent/CN111353204B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry

Landscapes

  • Chemical & Material Sciences (AREA)
  • Dispersion Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种岩心孔隙度参数场的处理方法，属于油层物理领域。该处理方法包括：对目标岩心孔隙度参数场进行格式转换，获取具有CMG数值模拟软件可识别格式的二维参数矩阵。对二维参数矩阵进行处理，获得具有CMG数值模拟软件可识别文件格式的可执行文件。利用CMG数值模拟软件，对可执行文件进行处理。该处理方法实现了实验理论参数与商用数值模拟软件的对接，将三维参数转化成可被直接编辑或计算的二维参数，如此，可利用CMG数值模拟软件对转化后的岩心孔隙度参数场进行实验模拟及预测，并且在实验过程中，可针对同一批岩心孔隙度参数场，设置不同的操作参数进行多次实验，提高对岩心孔隙度参数场的可利用性。

Description

一种岩心孔隙度参数场的处理方法
技术领域
本发明涉及油层物理领域，特别涉及一种岩心孔隙度参数场的处理方法。
背景技术
数字岩心是基于二维扫描电镜图像或三维CT扫描图像，运用计算机图像处理技术，通过一定的算法完成数字岩心重构，其能够最大程度地反映出地层信息，对于岩石物理的研究具有重要的意义。目前已公开了岩心孔隙度参数场的获取方法，其包括：对岩性样品进行CT扫描及处理，获得岩心二值化图像，对岩心二值化图像进行处理获得二维矩阵，将二维矩阵按三维顺序进行叠加，形成三维矩阵。将三维矩阵封装至单元数组的各个元胞数组中进行局部反演，获得每个元胞数组对应的孔隙度，进而获得单元数组的孔隙度参数场。根据单元数组的孔隙度参数场计算三维切割样品的平均孔隙度模拟值，并与岩心样品的平均孔隙度实测值进行比较，当两者的误差小于5％时，确定上述获得的单元数组的孔隙度参数场为岩心样品的岩心孔隙度参数场。
然而，目前所获得的岩心孔隙度参数场为三维数据，其无法被直接编辑或计算，须对应转化成可被直接编辑或计算的二维数据。但是，目前还未有将岩心孔隙度参数场转化成可编辑的二维数据的方法。
发明内容
鉴于此，本发明提供一种岩心孔隙度参数场的处理方法，可解决上述技术问题。具体而言，包括以下的技术方案：
一种岩心孔隙度参数场的处理方法，所述处理方法包括：对目标岩心孔隙度参数场进行格式转换，获取具有CMG数值模拟软件可识别格式的二维参数矩阵；
对所述二维参数矩阵进行处理，获得具有CMG数值模拟软件可识别文件格式的可执行文件；
利用CMG数值模拟软件，对所述可执行文件进行处理。
在一种可能的实现方式中，所述处理方法还包括：利用所述CMG数值模拟软件对所述可执行文件进行处理，获取所述二维参数矩阵对应的矩阵图像；
获取所述目标岩心孔隙度参数场对应的二值化矩阵图像；
将所述二维参数矩阵对应的矩阵图像与所述二值化矩阵图像进行对比，确定所述二维参数矩阵的可靠性。
在一种可能的实现方式中，所述对目标岩心孔隙度参数场进行格式转换，获取具有CMG数值模拟软件可识别格式的二维参数矩阵，包括：
对所述目标岩心孔隙度参数场的数据导入顺序进行定义；
按照定义的数据导入顺序，将所述目标岩心孔隙度参数场的各参数依次导入至平面二维矩阵中，得到所述二维参数矩阵。
在一种可能的实现方式中，选取所述目标岩心孔隙度参数场的任意点为起始原点，按照从左至右以及从上至下的顺序作为所述数据导入顺序。
在一种可能的实现方式中，所述目标岩心孔隙度参数场中，将相邻网格的N个近似数值定义成a*N，其中所述a为所述N个近似数值的取值。
在一种可能的实现方式中，所述近似数值包括：小数点后四位数值相等的数值。
在一种可能的实现方式中，通过MATLAB软件，将所述目标岩心孔隙度参数场的各参数依次导入至平面二维矩阵中，得到所述二维参数矩阵。
在一种可能的实现方式中，所述对所述二维参数矩阵进行处理，获得具有CMG数值模拟软件可识别文件格式的可执行文件，包括：
利用所述CMG数值模拟软件新建一空白文件夹；
在首次运行所述空白文件夹之前，将所述二维参数矩阵复制到所述空白文件夹内，获得所述可执行文件。
在一种可能的实现方式中，所述二维参数矩阵复制到所述空白文件夹内，获得所述可执行文件，包括：
在所述空白文件夹内找到关键字*POR；
将所述关键字*POR后面的数据删除；
将所述二维参数矩阵复制到所述关键字*POR后面，获得所述可执行文件。
本发明实施例提供的技术方案的有益效果至少包括：
本发明实施例提供的岩心孔隙度参数场的处理方法，通过将目标岩心孔隙度参数场对应的三维参数矩阵转化成能够被CMG数值模拟软件识别的二维参数矩阵，并将该二维参数矩阵转换成能够被CMG数值模拟软件识别的可执行文件，随后利用CMG数值模拟软件即可对该可执行文件进行处理。可见，本发明实施例提供的方法，实现了实验理论参数与商用数值模拟软件的对接，将三维参数转化成可被直接编辑或计算的二维参数，如此，可利用CMG数值模拟软件对转化后的岩心孔隙度参数场进行实验模拟及预测，并且在实验过程中，可针对同一批岩心孔隙度参数场，设置不同的操作参数进行多次实验，提高对岩心孔隙度参数场的可利用性。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为以左上角的起点作为原点定义数据导入方向的示意图；
图2为以中部某一点作为原点定义数据导入方向的示意图；
图3为一种目标岩心孔隙度参数场的示意图；
图4为一种二维参数矩阵的示意图；
图5为二维参数矩阵对应的矩阵图像与岩心孔隙度参数场的二值化矩阵图像的对比示意图。
其中，图5中位于左侧的图像为岩心孔隙度参数场的二值化矩阵图像，位于右侧的图像为二维参数矩阵对应的矩阵图像。
具体实施方式
为使本发明的技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。
需要说明的是，本发明实施例所涉及的岩心孔隙度参数场可以通过现有技术(例如，中国专利申请：岩心孔隙度参数场的获取方法，申请号：201610884265.1)提供的方法获取得到。该岩心孔隙度参数场为三维参数矩阵(即，由多个在三维空间内按顺序叠加的二维矩阵组成)。本发明实施例所提供的处理方法的目的是将该三维参数矩阵处理并转换成可直接编辑或计算的二维参数，具体涉及的步骤如下所示：
本发明实施例提供的岩心孔隙度参数场的处理方法包括以下步骤：对目标岩心孔隙度参数场进行格式转换，获取具有CMG数值模拟软件可识别格式的二维参数矩阵。
对二维参数矩阵进行处理，获得具有CMG数值模拟软件可识别文件格式的可执行文件。
利用CMG数值模拟软件，对可执行文件进行处理。
本发明实施例提供的岩心孔隙度参数场的处理方法，通过将目标岩心孔隙度参数场对应的三维参数矩阵转化成能够被CMG数值模拟软件识别的二维参数矩阵，并将该二维参数矩阵转换成能够被CMG数值模拟软件识别的可执行文件，随后利用CMG数值模拟软件即可对该可执行文件进行处理。可见，本发明实施例提供的方法，实现了实验理论参数与商用数值模拟软件的对接，将三维参数转化成可被直接编辑或计算的二维参数，如此，可利用CMG数值模拟软件对转化后的岩心孔隙度参数场进行实验模拟及预测，并且在实验过程中，可针对同一批岩心孔隙度参数场，设置不同的操作参数进行多次实验，提高对岩心孔隙度参数场的可利用性。
可以理解的是，本发明实施例所涉及的CMG数值模拟软件为本领域公知的商用数值模拟软件，其又称为CMG油藏数值模拟软件，本发明实施例在此对其不作更具体的描述。
其中，对目标岩心孔隙度参数场进行格式转换，获取具有CMG数值模拟软件可识别格式的二维参数矩阵，包括：对目标岩心孔隙度参数场的数据导入顺序进行定义。按照定义的数据导入顺序，将目标岩心孔隙度参数场的各参数依次导入至平面二维矩阵中，得到二维参数矩阵。
示例地，选取目标岩心孔隙度参数场的任意点为起始原点，按照从左至右以及从上至下的顺序作为数据导入顺序。
以下结合相关附图来进一步描述上述过程：
如附图1所示，对于该岩心孔隙度参数场来说，可以左上角的起点作为原点，按照从左至右(沿x方向)，以及从上至下(沿y方向)的顺序，将这些参数依次导入至平面二维矩阵中。
如附图2所示，对于该岩心孔隙度参数场来说，可以中部任一点作为原点，按照从左至右，以及从上至下的顺序，将这些参数依次导入至平面二维矩阵中。
进一步地，为了简化表达格式以节省运算空间，对于目标岩心孔隙度参数场中，将相邻网格的N个近似数值，定义成a*N，其中a为N个近似数值的取值。
举例来说，如附图3所示，第2部分的矩阵中，有三个相邻网格的数值均显示为0.065，以及，第4部分的矩阵中，有三个相邻网格的数值均显示为0.059，所以，当将目标岩心孔隙度参数场封装至平面二维矩阵的过程中，可以将这两部分相邻网格的三个数值分别定义成0.065*3和0.059*3(参见图4)。
根据所设定的数值精度来确定相邻网格的数值是否近似，本发明实施例中，上述近似数值包括：小数点后四位数值相等的数值。举例来说，对于以下三个数值来说，0.01752、0.01753、0.01751即可作为近似数值，且取值可选择为0.0175。
在将目标岩心孔隙度参数场导入至平面二维矩阵的过程中，为了提高导入效率，可以通过MATLAB软件，将目标岩心孔隙度参数场的各参数依次导入至平面二维矩阵中，得到二维参数矩阵。具体地，可以采用MATLAB软件进行程序编写(包括：上述就近似数值的定义，以及封装时的导入顺序等)，来将三维参数矩阵转化为能够被CMG数值模拟软件识别的二维参数矩阵。
在获得二维参数矩阵后，由于后续要采用CMG数值模拟软件进行数据处理，本发明实施例提供的方法需要对二维参数矩阵进行进一步处理，获得具有CMG数值模拟软件可识别文件格式的可执行文件，示例地，该过程可包括：利用CMG数值模拟软件新建一空白文件夹；在首次运行空白文件夹之前，将二维参数矩阵复制到空白文件夹内，获得可执行文件。
举例来说，CMG数值模拟软件的可执行文件为以.dat后缀结尾的文件，当建立一个新的空白文件夹时会自动生成，即，该空白文件夹以.dat后缀结尾，其包含了CMG数值模拟软件运行时所需要的全部数据体，在将二维参数矩阵复制到该空白文件夹之前，在该空白文件夹中找到关键字*POR，然后将该关键字后面的数据全部删除，随后将二维参数矩阵复制到该关键字后面即可。
进一步地，在利用本发明实施例提供的处理方法获得可被CMG数值模拟软件处理的可执行文件后，为了确定所获得二维参数矩阵是否准确(即可靠)，还可以利用CMG数值模拟软件对可执行文件进行处理，获取二维参数矩阵对应的矩阵图像；并且，获取目标岩心孔隙度参数场对应的二值化矩阵图像，将二维参数矩阵对应的矩阵图像与二值化矩阵图像进行对比，确定二维参数矩阵的可靠性。
附图5示例了二维参数矩阵对应的矩阵图像与岩心孔隙度参数场的二值化矩阵图像的对比示意图，由图5可知，两者达到了高度吻合，表明了本发明实施例提供的处理方法的可行性。
以上所述仅是为了便于本领域的技术人员理解本发明的技术方案，并不用以限制本发明。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (7)

1.一种岩心孔隙度参数场的处理方法，其特征在于，所述处理方法包括：对目标岩心孔隙度参数场进行格式转换，获取具有CMG数值模拟软件可识别格式的二维参数矩阵；
对所述二维参数矩阵进行处理，获得具有CMG数值模拟软件可识别文件格式的可执行文件；
利用CMG数值模拟软件，对所述可执行文件进行处理；
所述对目标岩心孔隙度参数场进行格式转换，获取具有CMG数值模拟软件可识别格式的二维参数矩阵，包括：
对所述目标岩心孔隙度参数场的数据导入顺序进行定义；
按照定义的数据导入顺序，将所述目标岩心孔隙度参数场的各参数依次导入至平面二维矩阵中，得到所述二维参数矩阵；
所述对所述二维参数矩阵进行处理，获得具有CMG数值模拟软件可识别文件格式的可执行文件，包括：
利用所述CMG数值模拟软件新建一空白文件夹；
在首次运行所述空白文件夹之前，将所述二维参数矩阵复制到所述空白文件夹内，获得所述可执行文件。
2.根据权利要求1所述的岩心孔隙度参数场的处理方法，其特征在于，所述处理方法还包括：利用所述CMG数值模拟软件对所述可执行文件进行处理，获取所述二维参数矩阵对应的矩阵图像；
获取所述目标岩心孔隙度参数场对应的二值化矩阵图像；
将所述二维参数矩阵对应的矩阵图像与所述二值化矩阵图像进行对比，确定所述二维参数矩阵的可靠性。
3.根据权利要求1所述的岩心孔隙度参数场的处理方法，其特征在于，选取所述目标岩心孔隙度参数场的任意点为起始原点，按照从左至右以及从上至下的顺序作为所述数据导入顺序。
4.根据权利要求1所述的岩心孔隙度参数场的处理方法，其特征在于，所述目标岩心孔隙度参数场中，将相邻网格的N个近似数值定义成a*N，其中所述a为所述N个近似数值的取值。
5.根据权利要求4所述的岩心孔隙度参数场的处理方法，其特征在于，所述近似数值包括：小数点后四位数值相等的数值。
6.根据权利要求3所述的岩心孔隙度参数场的处理方法，其特征在于，通过MATLAB软件，将所述目标岩心孔隙度参数场的各参数依次导入至平面二维矩阵中，得到所述二维参数矩阵。
7.根据权利要求1所述的岩心孔隙度参数场的处理方法，其特征在于，所述二维参数矩阵复制到所述空白文件夹内，获得所述可执行文件，包括：
在所述空白文件夹内找到关键字*POR；
将所述关键字*POR后面的数据删除；
将所述二维参数矩阵复制到所述关键字*POR后面，获得所述可执行文件。
CN201811562408.2A 2018-12-20 2018-12-20 一种岩心孔隙度参数场的处理方法 Active CN111353204B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811562408.2A CN111353204B (zh) 2018-12-20 2018-12-20 一种岩心孔隙度参数场的处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811562408.2A CN111353204B (zh) 2018-12-20 2018-12-20 一种岩心孔隙度参数场的处理方法

Publications (2)

Publication Number Publication Date
CN111353204A CN111353204A (zh) 2020-06-30
CN111353204B true CN111353204B (zh) 2023-10-31

Family

ID=71191958

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811562408.2A Active CN111353204B (zh) 2018-12-20 2018-12-20 一种岩心孔隙度参数场的处理方法

Country Status (1)

Country Link
CN (1) CN111353204B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106226831A (zh) * 2015-08-25 2016-12-14 中国石油天然气股份有限公司 一种岩心的孔隙度矩阵反演方法
CN106525680A (zh) * 2016-10-10 2017-03-22 中国石油天然气股份有限公司 岩心孔隙度参数场的获取方法
CN107341848A (zh) * 2017-07-08 2017-11-10 青岛科技大学 一种将岩心ct图像处理为商用cfd软件可读文件的新方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10465483B2 (en) * 2013-08-13 2019-11-05 Schlumberger Technology Corporation Digital core sensitivity analysis

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106226831A (zh) * 2015-08-25 2016-12-14 中国石油天然气股份有限公司 一种岩心的孔隙度矩阵反演方法
CN106525680A (zh) * 2016-10-10 2017-03-22 中国石油天然气股份有限公司 岩心孔隙度参数场的获取方法
CN107341848A (zh) * 2017-07-08 2017-11-10 青岛科技大学 一种将岩心ct图像处理为商用cfd软件可读文件的新方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
一种基于数字岩心技术的岩石等效电参数计算方法;靳钊;《物探与化探》;第42卷(第6期);全文 *
数字岩心建模方法及应用综述;郑启斌;《能源与环保》;第39卷(第12期);全文 *
数字岩心建模方法研究现状及展望;林承焰;《地球物理学进展》;第33卷(第2期);全文 *

Also Published As

Publication number Publication date
CN111353204A (zh) 2020-06-30

Similar Documents

Publication Publication Date Title
CN110609693B (zh) 基于数据标准化的代码更新方法、装置及终端设备
KR100928382B1 (ko) 메인프레임 시스템의 데이터베이스를 개방형 시스템에적합한 데이터베이스로 변환하는 마이그레이션 장치 및 그방법
EP1989685B1 (en) A method for comparing a first computer-aided 3d model with a second computer-aided 3d model
US20180129633A1 (en) Copying text format pattern
US11003635B2 (en) Database scheme for storing generic data
Rizzi et al. Simulation of time evolution with multiscale entanglement renormalization ansatz
Bahubalendruni et al. An efficient stable subassembly identification method towards assembly sequence generation
CN104572153A (zh) 一种用于车辆刷新的刷新数据转换方法
US20190220441A1 (en) Method, device and computer program product for data migration
CN111158680A (zh) 一种页面构建的方法、装置、设备及存储介质
CN111694811B (zh) 一种批量数据入库方法及装置
CN114564446A (zh) 一种文件存储的方法、装置、系统及存储介质
CN115795629A (zh) 数据转换方法、数据转换系统以及电子设备
CN111353204B (zh) 一种岩心孔隙度参数场的处理方法
CN109472364B (zh) 量子程序的处理方法及装置、存储介质和电子装置
CN113609736B (zh) 一种基于孔裂隙结构数字图像的数值计算模型构建方法
CN116894272B (zh) 一种基于高速加密技术的云计算系统数据处理方法
CN116595918A (zh) 快速逻辑等价性验证方法、装置、设备及存储介质
CN111176901A (zh) 一种hdfs删除文件恢复方法、终端设备及存储介质
CN112329267B (zh) 一种基于特征线法的组合几何中子输运处理方法及装置
Dinklage et al. Bidirectional text compression in external memory
CN114265958A (zh) Kv ssd的映射管理方法、装置及存储介质
CN115221133A (zh) 一种面向3d打印的stl文件处理方法
Maabreh et al. A multithreading and hashing technique for indexing Target‐Decoy peptides databases
Murtiyoso et al. Automatic training data generation in Deep Learning-aided semantic segmentation of Heritage buildings

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