砂岩孔渗参数的确定方法和装置制造方法及图纸

本发明专利技术提供了一种砂岩孔渗参数的确定方法和装置，其中，该方法可以包括：根据待测砂岩弛豫时间下的核磁共振数据，建立满足预定维数的候选矩阵；建立满足预定维数的二维矩阵，根据二维矩阵中各个待测元素得到关联维数的表达式；将表达式的斜率作为关联维数；基于关联维数，计算确定待测砂岩的孔渗参数。在本发明专利技术实施例中，基于待测砂岩多孔介质中的流体、孔隙表面相互作用均与流体的核磁共振特性有关的特性，通过测定流体的核磁共振特性数据来辨识不同大小的孔隙，利用由表征核磁共振数据之间关系的参数所组成的矩阵确定核磁共振数据关联维数，并确定待测砂岩的孔渗参数，提供了一种确定孔渗参数的新方法，并提高了砂岩孔渗参数的识别精度。

Method and device for determining porosity and permeability parameters of sandstone

The present invention provides a method and apparatus for determining a parameter, permeability sandstone hole, the method can include: according to the NMR data of sandstone under relaxation time tested, built to meet the predetermined candidate matrix dimension; set up to meet the predetermined two-dimensional matrix dimension, according to various two-dimensional matrix elements to be detected to obtain expressions for the correlation dimension the slope as the expression; correlation dimension; based on correlation dimension, calculate the measured porosity and permeability of sandstone. In the embodiment of the invention, the characteristics of sandstone porous medium fluid and pore surface interaction are related to NMR properties of the fluid to be measured based on the characteristics of the data by NMR determination of fluid to identify the pores of different sizes to determine the NMR data correlation dimension using matrix composed of NMR spectra according to the number of parameters the relationship between the measured and determined the sandstone porosity and permeability parameters, provides a new method to determine the porosity and permeability, and improve the recognition accuracy of parameters of sandstone pore permeability.

【技术实现步骤摘要】
砂岩孔渗参数的确定方法和装置
本专利技术涉及地质勘探
，特别涉及一种砂岩孔渗参数的确定方法和装置。
技术介绍
岩石孔渗参数包括岩石的孔隙度和渗透率，可以指的是岩石所具有的孔隙和吼道的几何形状、大小、分布及其互连通关系。岩石的孔隙度和渗透率是岩石储层结构的重要组成部分，也是影响油气储集与油田开发效果的内在因素。在油气开发微观实验中测定不同尺度岩心样品的孔渗参数时，可以通过测定结果了解待测储层的物性特征，为制定合理的油气开采方式提供直接依据。砂岩气藏具有独特的开发特征，储层均质性良好，具备可观的储量和开发条件。目前，对砂岩气藏储层孔渗参数的测定主要是采用常规的岩心压汞或者CT扫描岩心等技术手段来分析砂岩的孔隙结构，并测定砂岩的渗透率和孔隙度，然而，使用上述常规方法测得的砂岩孔渗参数准确度不高。针对上述问题，目前尚未提出有效的解决方案。
技术实现思路
本专利技术提供了一种砂岩孔渗参数的确定方法和装置，以达到提高砂岩孔渗参数准确度的目的。本专利技术实施例提供了一种砂岩孔渗参数的确定方法，可以包括：根据获取的待测砂岩弛豫时间下的核磁共振数据，建立满足预定维数的候选矩阵，所述预定维数大于等于1且小于等于所述核磁共振数据的个数；建立满足预定维数的二维矩阵，所述二维矩阵中的待测元素是根据所述候选矩阵中与所述待测元素的行数和列数对应的核磁共振数据之间的欧式距离确定的；计算所述待测元素中满足小于第一阈值的第一距离数目、所述待测元素中满足大于第一阈值的第二距离数目，将所述第一距离数目与所述第二距离数目相加得到距离总数，计算所述第一距离数目和所述距离总数的比值；基于所述比值，得到关联维数的表达式，所述表达式表征的是在第二阈值范围内，所述比值的对数等于所述第一阈值的对数与所述关联维数的乘积；计算所述比值的对数与所述第一阈值的对数成线性关系时所述表达式的斜率，将所述斜率作为所述关联维数；基于所述关联维数，计算确定所述待测砂岩的孔渗参数。在一个实施例中，所述孔渗参数可以包括：孔隙度；相应的，可以按照以下公式计算确定所述待测砂岩的孔渗参数：上式中，φ表示所述待测砂岩的孔隙度，Dr表示所述关联维数，a为第一实数，b为第二实数。在一个实施例中，所述孔渗参数可以包括：渗透率；相应的，可以按照以下公式计算确定所述待测砂岩的孔渗参数：上式中，k表示所述待测砂岩的渗透率，Dr表示所述关联维数，c为第三实数，d为第四实数。在一个实施例中，计算所述比值的对数与所述第一阈值的对数成线性关系时所述表达式的斜率，将所述斜率作为所述关联维数，可以包括：分别计算所述表达式关于所述第一阈值范围的一阶导数和二阶导数；将所述一阶导数为常数并且所述二阶导数为0时，所述表达式中第一阈值对应的区间，作为所述第二阈值范围；在所述第二阈值范围中，对所述表达式进行拟合，并确定拟合时所对应的拟合度；选取所述拟合度最大时的区间作为最佳区间，与所述最佳区间对应的一阶导数作为所述待测砂岩的关联维数。在一个实施例中，可以按照以下方式确定所述预定维数：从1逐渐增加维数，根据获取的待测砂岩弛豫时间下的核磁共振数据，分别建立列数满足所述维数的矩阵；计算所述矩阵对应的伪邻点值；计算相邻维数对应的伪邻点值之间的差值；当所述差值满足预设的第三阈值范围时，将所述差值对应的矩阵作为候选矩阵，将所述差值对应的维数作为预定维数。在一个实施例中，根据获取的待测砂岩弛豫时间下的核磁共振数据，分别建立列数满足所述维数的矩阵，可以包括：从所述核磁共振数据中依次选取和所述维数相等数目的数据作为第一行向量；步骤1：以所述第一行向量中的第二个数据作为第一个数据，从所述核磁共振数据中依次选取和所述维数相等数目的数据作为第二行向量；步骤2：将所述第二行向量作为第一行向量；重复执行步骤1至步骤2，直至所述核磁共振数据均存在于所述矩阵中。在一个实施例中，可以按照以下公式计算所述矩阵对应的伪邻点值：上式中，Rm表示所述待测砂岩的伪邻点值，N表示所述核磁共振数据的个数，m表示所述维数，m为正整数，1≤m≤N，Yi(m)表示列数为m的矩阵中第i个行向量，Yn(m)表示所述矩阵中的所有m维行向量中与Yi(m)欧式距离最小的行向量。本专利技术实施例还提供了一种砂岩孔渗参数的确定装置，可以包括：矩阵建立模块，可以用于根据获取的待测砂岩弛豫时间下的核磁共振数据，建立满足预定维数的候选矩阵，所述预定维数大于等于1且小于等于所述核磁共振数据的个数；二维矩阵建立模块，可以用于建立满足预定维数的二维矩阵，所述二维矩阵中的待测元素是根据所述候选矩阵中与所述待测元素的行数和列数对应的核磁共振数据之间的欧式距离确定的；比值计算模块，可以用于计算所述待测元素中满足小于第一阈值的第一距离数目、所述待测元素中满足大于第一阈值的第二距离数目，将所述第一距离数目与所述第二距离数目相加得到距离总数，计算所述第一距离数目和所述距离总数的比值；表达式确定模块，可以用于基于所述比值，得到关联维数的表达式，所述表达式表征的是在第二阈值范围内，所述比值的对数等于所述第一阈值的对数与所述关联维数的乘积；关联维数计算模块，可以用于计算所述比值的对数与所述第一阈值的对数成线性关系时所述表达式的斜率，将所述斜率作为所述关联维数；孔渗参数计算模块，可以用于基于所述关联维数，计算确定所述待测砂岩的孔渗参数。在一个实施例中，所述孔渗参数可以包括：孔隙度；相应的，所述孔渗参数计算模块具体可以用于按照以下公式计算确定所述待测砂岩的孔渗参数：上式中，φ表示所述待测砂岩的孔隙度，Dr表示所述关联维数，a为第一实数，b为第二实数。在一个实施例中，所述孔渗参数可以包括：渗透率；相应的，所述参数计算模块具体可以用于按照以下公式计算确定所述待测砂岩的孔渗参数：上式中，k表示所述待测砂岩的渗透率，Dr表示所述关联维数，c为第三实数，d为第四实数。本专利技术实施例还提供了一种砂岩孔渗参数的确定装置，可以包括处理器以及用于存储处理器可执行指令的存储器，所述处理器执行所述指令时可以实现：根据获取的待测砂岩弛豫时间下的核磁共振数据，建立满足预定维数的候选矩阵，所述预定维数大于等于1且小于等于所述核磁共振数据的个数；建立满足预定维数的二维矩阵，所述二维矩阵中的待测元素是根据所述候选矩阵中与所述待测元素的行数和列数对应的核磁共振数据之间的欧式距离确定的；计算所述待测元素中满足小于第一阈值的第一距离数目、所述待测元素中满足大于第一阈值的第二距离数目，将所述第一距离数目与所述第二距离数目相加得到距离总数，计算所述第一距离数目和所述距离总数的比值；基于所述比值，得到关联维数的表达式，所述表达式表征的是在第二阈值范围内，所述比值的对数等于所述第一阈值的对数与所述关联维数的乘积；计算所述比值的对数与所述第一阈值的对数成线性关系时所述表达式的斜率，将所述斜率作为所述关联维数；基于所述关联维数，计算确定所述待测砂岩的孔渗参数。本专利技术实施例还提供了一种计算机可读存储介质，其上存储有计算机指令，所述指令被执行时可以实现以下步骤：根据获取的待测砂岩弛豫时间下的核磁共振数据，建立满足预定维数的候选矩阵，所述预定维数大于等于1且小于等于所述核磁共振数据的个数；建立满足预定维数的二维矩阵，所述二维矩阵中的待测元素是根据所述候选矩阵中与所述待测元素的本文档来自技高网...

【技术保护点】
一种砂岩孔渗参数的确定方法，其特征在于，包括：根据获取的待测砂岩弛豫时间下的核磁共振数据，建立满足预定维数的候选矩阵，所述预定维数大于等于1且小于等于所述核磁共振数据的个数；建立满足预定维数的二维矩阵，所述二维矩阵中的待测元素是根据所述候选矩阵中与所述待测元素的行数和列数对应的核磁共振数据之间的欧式距离确定的；计算所述待测元素中满足小于第一阈值的第一距离数目、所述待测元素中满足大于第一阈值的第二距离数目，将所述第一距离数目与所述第二距离数目相加得到距离总数，计算所述第一距离数目和所述距离总数的比值；基于所述比值，得到关联维数的表达式，所述表达式表征的是在第二阈值范围内，所述比值的对数等于所述第一阈值的对数与所述关联维数的乘积；计算所述比值的对数与所述第一阈值的对数成线性关系时所述表达式的斜率，将所述斜率作为所述关联维数；基于所述关联维数，计算确定所述待测砂岩的孔渗参数。

【技术特征摘要】
1.一种砂岩孔渗参数的确定方法，其特征在于，包括：根据获取的待测砂岩弛豫时间下的核磁共振数据，建立满足预定维数的候选矩阵，所述预定维数大于等于1且小于等于所述核磁共振数据的个数；建立满足预定维数的二维矩阵，所述二维矩阵中的待测元素是根据所述候选矩阵中与所述待测元素的行数和列数对应的核磁共振数据之间的欧式距离确定的；计算所述待测元素中满足小于第一阈值的第一距离数目、所述待测元素中满足大于第一阈值的第二距离数目，将所述第一距离数目与所述第二距离数目相加得到距离总数，计算所述第一距离数目和所述距离总数的比值；基于所述比值，得到关联维数的表达式，所述表达式表征的是在第二阈值范围内，所述比值的对数等于所述第一阈值的对数与所述关联维数的乘积；计算所述比值的对数与所述第一阈值的对数成线性关系时所述表达式的斜率，将所述斜率作为所述关联维数；基于所述关联维数，计算确定所述待测砂岩的孔渗参数。2.如权利要求1所述的方法，其特征在于，所述孔渗参数包括：孔隙度；相应的，按照以下公式计算确定所述待测砂岩的孔渗参数：上式中，φ表示所述待测砂岩的孔隙度，Dr表示所述关联维数，a为第一实数，b为第二实数。3.如权利要求1所述的方法，其特征在于，所述孔渗参数包括：渗透率；相应的，按照以下公式计算确定所述待测砂岩的孔渗参数：上式中，k表示所述待测砂岩的渗透率，Dr表示所述关联维数，c为第三实数，d为第四实数。4.如权利要求1所述的方法，其特征在于，计算所述比值的对数与所述第一阈值的对数成线性关系时所述表达式的斜率，将所述斜率作为所述关联维数，包括：分别计算所述表达式关于所述第一阈值范围的一阶导数和二阶导数；将所述一阶导数为常数并且所述二阶导数为0时，所述表达式中第一阈值对应的区间，作为所述第二阈值范围；在所述第二阈值范围中，对所述表达式进行拟合，并确定拟合时所对应的拟合度；选取所述拟合度最大时的区间作为最佳区间，与所述最佳区间对应的一阶导数作为所述待测砂岩的关联维数。5.如权利要求1所述的方法，其特征在于，按照以下方式确定所述预定维数：从1逐渐增加维数，根据获取的待测砂岩弛豫时间下的核磁共振数据，分别建立列数满足所述维数的矩阵；计算所述矩阵对应的伪邻点值；计算相邻维数对应的伪邻点值之间的差值；当所述差值满足预设的第三阈值范围时，将所述差值对应的矩阵作为候选矩阵，将所述差值对应的维数作为预定维数。6.如权利要求5所述的方法，其特征在于，根据获取的待测砂岩弛豫时间下的核磁共振数据，分别建立列数满足所述维数的矩阵，包括：从所述核磁共振数据中依次选取和所述维数相等数目的数据作为第一行向量；步骤1：以所述第一行向量中的第二个数据作为第一个数据，从所述核磁共振数据中依次选取和所述维数相等数目的数据作为第二行向量；步骤2：将所述第二行向量作为第一行向量；重复执行步骤1至步骤2，直至所述核磁共振数据均存在于所述矩阵中。7.如权利要求5所述的方法，其特征在于，按照以下公式计算所述矩阵对应的伪邻点值：上式中，Rm表示所述待测砂岩的伪邻点值，N表示所述核磁共振数据的个数，m表示所述维数，m为正整数，1≤m≤N，Yi(m)表示列数为m的矩阵中第i个行向量，Yn(m)表示所述矩阵中的所有m维行向量中与Yi(m)欧式距离最小的行向量。8.一种砂岩孔渗参数的确定装置，其特征在于，包括：矩阵建立模块，用于根据获取的待测砂岩弛豫时间下的核磁共振...

【专利技术属性】
技术研发人员：吴国铭，叶礼友，高树生，刘华勋，朱文卿，安为国，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：北京,11