CN111999331B - 基于核磁测试的毛管压力确定方法及装置 - Google Patents
基于核磁测试的毛管压力确定方法及装置 Download PDFInfo
- Publication number
- CN111999331B CN111999331B CN201910444319.6A CN201910444319A CN111999331B CN 111999331 B CN111999331 B CN 111999331B CN 201910444319 A CN201910444319 A CN 201910444319A CN 111999331 B CN111999331 B CN 111999331B
- Authority
- CN
- China
- Prior art keywords
- capillary pressure
- core
- nuclear magnetic
- test data
- magnetic test
- 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
- 238000012360 testing method Methods 0.000 title claims abstract description 215
- 230000005291 magnetic effect Effects 0.000 title claims abstract description 164
- 238000000034 method Methods 0.000 title claims abstract description 94
- 239000011148 porous material Substances 0.000 claims abstract description 126
- 238000004364 calculation method Methods 0.000 claims abstract description 75
- 238000006243 chemical reaction Methods 0.000 claims description 57
- 239000011435 rock Substances 0.000 claims description 33
- 238000005481 NMR spectroscopy Methods 0.000 claims description 32
- 238000004590 computer program Methods 0.000 claims description 17
- 230000035699 permeability Effects 0.000 claims description 14
- 238000003860 storage Methods 0.000 claims description 13
- 239000011159 matrix material Substances 0.000 claims description 11
- 239000002131 composite material Substances 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 21
- 230000006870 function Effects 0.000 description 19
- 238000009826 distribution Methods 0.000 description 17
- 238000004891 communication Methods 0.000 description 15
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 13
- 229910052753 mercury Inorganic materials 0.000 description 13
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 13
- 238000011161 development Methods 0.000 description 11
- 238000012545 processing Methods 0.000 description 10
- 230000000694 effects Effects 0.000 description 9
- 238000012795 verification Methods 0.000 description 9
- 238000001228 spectrum Methods 0.000 description 7
- 239000012530 fluid Substances 0.000 description 6
- 230000005311 nuclear magnetism Effects 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000009792 diffusion process Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 238000003825 pressing Methods 0.000 description 5
- 238000004457 water analysis Methods 0.000 description 5
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 4
- 239000000243 solution Substances 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 229920006395 saturated elastomer Polymers 0.000 description 3
- 241000272525 Anas platyrhynchos Species 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 239000004927 clay Substances 0.000 description 2
- 239000013078 crystal Substances 0.000 description 2
- 229910000514 dolomite Inorganic materials 0.000 description 2
- 239000010459 dolomite Substances 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000000047 product Substances 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 238000009736 wetting Methods 0.000 description 2
- 208000010392 Bone Fractures Diseases 0.000 description 1
- 206010017076 Fracture Diseases 0.000 description 1
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 241000282320 Panthera leo Species 0.000 description 1
- 108010001267 Protein Subunits Proteins 0.000 description 1
- 239000012223 aqueous fraction Substances 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000002902 bimodal effect Effects 0.000 description 1
- 238000005119 centrifugation Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 239000003921 oil Substances 0.000 description 1
- 230000005298 paramagnetic effect Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 102000030938 small GTPase Human genes 0.000 description 1
- 108060007624 small GTPase Proteins 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
- G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
- G01N24/081—Making measurements of geologic samples, e.g. measurements of moisture, pH, porosity, permeability, tortuosity or viscosity
-
- 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
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Pathology (AREA)
- Immunology (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- High Energy & Nuclear Physics (AREA)
- Dispersion Chemistry (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Bioinformatics & Computational Biology (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Geology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Geochemistry & Mineralogy (AREA)
- Measuring Magnetic Variables (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种基于核磁测试的毛管压力确定方法及装置,该方法包括:测量待计算毛管压力地区的核磁测试数据;利用所述待计算毛管压力地区的核磁测试数据及预先建立的多个不同孔隙结构类别岩心的毛管压力计算模型,生成所述待计算毛管压力地区的毛管压力。本发明可以提供一种精度高、且无需进行岩心毛管压力测试即可获得毛管压力的方法。
Description
技术领域
本发明涉及石油勘探领域,尤其是核磁技术评价孔隙结构领域,具体涉及一种基于核磁测试的毛管压力确定方法及装置。
背景技术
储层孔隙结构研究多利用毛管压力曲线测试,通过孔隙结构特征参数来评价孔隙结构,但对于相对致密低渗储层,常规压汞测试进汞饱和度普遍较低,压汞曲线不能反映进汞压力高的微小孔喉特征。
核磁被广泛用来研究孔隙结构特征,目前的做法是将核磁弛豫时间T2谱转换为毛管压力曲线来进行孔隙结构的研究。由核磁弛豫机制可知,在均匀磁场中测量的横向弛豫时间T2为:
式中,T2B为流体的体积(自由)弛豫时间,ms;D为扩散系数,μm2/ms;G为磁场梯度,gauss/cm;TE为回波间隔,ms;S为孔隙的表面积;V为孔隙的体积;ρ2为岩石的横向表面弛豫强度,m/ms。
T2B的数值通常在2~3s,要比T2大得多,即T2B>T2,因此式(1)中右边的第一项可忽略,当磁场很均匀时(对应G很小),且TE足够小时,式(1)中右边的第三项也可忽略,于是有:
得到T2与孔径rc的关系式为:
式中,Fs称为几何形状因子,对球状孔隙,Fs=3;对柱状管道,Fs=2,由式(3)可见,孔隙内流体的弛豫时间和孔隙空间大小及形状有关由物理学可知,毛管压力与毛管孔径之间的关系为:
式中,pc为毛管压力,MPa;σ为流体界面张力;θ为润湿接触角;rc为毛管半径,μm由式(3)和式(4),得:
于是
目前各种利用T2分布评价毛管压力曲线的方法大都是基于式(6)进行的一种线性转换。鉴于线性方式转换效果不理想,何雨丹等人实际地层中孔隙结构很复杂,比表面与孔径成非线性关系提出了用幂函数构造毛管压力曲线的新方法。
此外Volokitin(2003)提出了适用于中低孔渗储层的毛管压力曲线和核磁T2分布之间的经验公式:
Pc毛管压力,MPa;T2核磁弛豫时间,ms;K渗透率;m,n,A,B,C,D为常数。
核磁弛豫时间转换为毛管压力的主要的方法包括:自由水分析法,分段幂函数刻度法,二维等面积刻度转换系数法等。何雨丹博士通过对压汞法、隔板法获取毛管压力曲线的实验机理以及储集层孔隙空间束缚水模型进行深入分析认为,构造核磁毛管压力曲线与实际压汞毛管压力曲线大孔隙部分吻合较好而在小孔隙部分出现分叉的现象;对于采用双组分分析的NMR测井而言,大孔隙表面吸附的束缚水膜和孔隙中间的自由流体(自由水)的弛豫速率完全不同,通过对实际岩心的实验分析发现,当薄膜束缚水含量较多时,大孔隙空间的薄膜束缚水的弛豫信号会累加到T2分布中较小的T2弛豫组分中,基于以上分析,何雨丹等人提出了一种利用隔板法获取的毛管压力曲线与100%饱含水的核磁测井T2分布中的自由水部分进行匹配以构造核磁毛管压力曲线的方法(自由水分析法),先测量两种状态下(分别为离心束缚水状态和100%饱含水状态)的核磁测井T2分布,并将离心T2分布从100%饱含水状态下的T2分布中消除掉,以得到自由水部分的T2分布;接着将自由水部分的T2分布从T2弛豫时间最大值向最小值进行反向累加,得到一条自由水部分的T2谱积分曲线,构造核磁毛管压力曲线。
分段幂函数刻度法是何雨丹等人基于自由水分析法的基础上提出的一种核磁毛管压力曲线构造方法的改进方法,该方法主要考虑到自由水分析法只构造了大孔隙自由水部分的毛管压力曲线,而忽略掉了小孔隙部分的孔隙结构特征的问题,分段幂函数刻度法主要考虑到当储集层孔径分布范围较宽时,地层中可能存在着属于快速及受限扩散特性的小孔隙,其体积弛豫可以忽略,也可能存在接近于慢扩散特性的较大孔隙,其体积弛豫不可忽略,因此在利用核磁测井T2分布构造毛管压力曲线的过程中应对T2时间采用分段转换刻度,即将T2分布分成2个区域,设区域分界点为体积弛豫可否忽略的分界点,对快速及受限扩散特性的小孔隙部分与接近于慢扩散特性的考虑体积弛豫影响的大孔隙部分,分别采用不同的转换刻度函数来构造核磁毛管压力曲线,通过同时进行了核磁测井和压汞实验的岩心岩心的分析,何雨丹等人提出了分段幂函数构造核磁毛管压力曲线的方法,该方法指出,核磁测井T2谱积分曲线与压汞毛管压力曲线之间的转换刻度关系并不是线性刻度关系,而是幂函数刻度关系,同时,在对实际核磁测井资料进行处理的过程中,对于T2分布呈单峰分布时,采用单一幂函数来构造核磁毛管压力曲线,而对于T2分布呈双峰分布时,大孔隙部分和小孔隙部分应分别采用不同幂函数来分段构造核磁毛管压力曲线。
邵维志等人(2009)针对自由水分析法和分段幂函数刻度方法存在的问题,提出了利用二维等面积刻度转换系数法来构造核磁毛管压力曲线,利用二维等面积刻度转换系数法构造核磁毛管压力曲线。
丛云海等(2013)选取苏里格气田致密砂岩岩样进行核磁T2谱测量与压汞毛细管压力曲线测试,寻找T2弛豫时间与平均孔喉半径的关系,并以平均孔喉半径为桥梁,通过建立核磁孔隙结构模型,对核磁T2谱进行了三组分分解,得到了岩石中不同形态孔隙的T2谱贡献值;其认为致密砂岩的孔隙由球型、管型与槽型孔隙共同构成,不同的组合关系反映着岩石的不同流动形态,其关系直接制约着致密砂岩孔隙结构的评价,并建立了反应球型、管型与槽型特征的孔喉半径与岩石的T2弛豫时间幂指数函数关系。
上述毛管压力测试方法均存在着精度差,且不能应用于岩心较少的勘探地区的问题。
发明内容
针对现有技术中的问题,本发明能够建立一种在复杂储层中同时考虑孔隙度及孔隙结构的准确计算含油饱和度的方法及装置。大大提高了饱和度计算精度和解释符合率,以达到常规测井系列高精度计算复杂储层饱和度的目的,有效解决储层饱和度定量评价问题。
为解决上述技术问题,本发明提供以下技术方案:
第一方面,本发明提供一种基于核磁测试的毛管压力确定方法,包括:
测量待计算毛管压力地区的核磁测试数据;
利用待计算毛管压力地区的核磁测试数据及预先建立的多个不同孔隙结构类别岩心的毛管压力计算模型,生成待计算毛管压力地区的毛管压力。
优选地,基于核磁测试的毛管压力确定方法,还包括:
根据孔隙结构特征参数,对岩心进行分类;
对分类后的岩心进行核磁测试及毛管压力测试,得到岩心分类后的核磁测试数据及毛管压力测试数据;
根据岩心分类后的核磁测试数据及毛管压力测试数据,利用核磁测试数据与毛管压力转换公式,生成不同类别岩心的毛管压力计算模型。
优选地,根据岩心分类后的核磁测试数据及毛管压力测试数据,利用核磁测试数据与毛管压力转换公式,生成不同类别岩心的毛管压力计算模型,包括:
将岩心分类后的核磁测试数据及毛管压力测试数据进行拟合,获得核磁测试数据与毛管压力转换公式的变量系数;
根据变量系数及核磁测试数据与毛管压力转换公式,生成岩心的毛管压力计算模型。
优选地,孔隙结构特征参数包括:孔隙类型、孔隙度及渗透率。
优选地,核磁测试数据包括:核磁共振弛豫时间。
第二方面,本发明提供一种基于核磁测试的毛管压力确定装置,该装置包括:
核磁测试数据测量单元,用于测量待计算毛管压力地区的核磁测试数据;
毛管压力生成单元,用于利用待计算毛管压力地区的核磁测试数据及预先建立的多个不同孔隙结构类别岩心的毛管压力计算模型,生成待计算毛管压力地区的毛管压力。
优选地,毛管压力确定装置,还包括:
岩心分类单元,用于根据孔隙结构特征参数,对岩心进行分类;
测试单元,用于对分类后的岩心进行核磁测试及毛管压力测试,得到岩心分类后的核磁测试数据及毛管压力测试数据;
毛管压力计算模型生成单元,用于根据岩心分类后的核磁测试数据及毛管压力测试数据,利用核磁测试数据与毛管压力转换公式,生成不同类别岩心的毛管压力计算模型。
优选地,毛管压力计算模型生成单元包括:
变量系数获得模块,用于将岩心分类后的核磁测试数据及毛管压力测试数据进行拟合,获得核磁测试数据与毛管压力转换公式的变量系数;
毛管压力计算模型模块,用于根据变量系数及核磁测试数据与毛管压力转换公式,生成岩心的毛管压力计算模型。
第三方面,本发明提供一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,处理器执行程序时实现基于核磁测试的毛管压力确定方法的步骤。
第四方面,本发明提供一种计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现基于核磁测试的毛管压力确定方法的步骤。
从上述描述可知,本发明提供基于核磁测试的毛管压力确定方法及装置,首先在岩心分类的基础上,对不同孔隙类型的典型岩心分别测试其核磁共振和毛管压力,将核磁测试数据与毛管压力进行拟合,进而确定不同孔隙类型的核磁测试数据与毛管压力转换公式的变量系数;最终确定不同孔隙类型的毛管压力计算模型。该方法主要是针对岩心取心有限,或者没有取心的地区,无法通过岩心的毛管压力曲线测试结果来分析孔隙结构特征的情况下,直接将岩心或测井核磁共振测试结果通过核磁测试数据与毛管压力转换公式转换为毛管曲线直接来分析孔隙结构特征。综上,本发明可以提供一种精度高、且无需进行岩心毛管压力测试即可获得毛管压力的方法。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本申请实施例的一种基于核磁测试的毛管压力确定系统的第一种结构示意图;
图2为本申请实施例的一种基于核磁测试的毛管压力确定系统的第二种结构示意图;
图3为本发明的实施例中方法10的流程示意图;
图4为本发明的实施例中方法20的流程示意图;
图5为本发明的实施例中步骤S300的流程示意图;
图6为本发明的具体应用实例中基于核磁测试的毛管压力确定方法的流程示意图;
图7为本发明的具体应用实例中第一类岩心铸体薄片;
图8为本发明的具体应用实例中第一类岩心压汞孔喉半径直方图;
图9为本发明的具体应用实例中第一类岩心核磁共振T2分布图;
图10为本发明的具体应用实例中第二类岩心示意图;
图11为本发明的具体应用实例中第二类岩心铸体薄片;
图12为本发明的具体应用实例中第二类岩心压汞孔喉半径直方图;
图13为本发明的具体应用实例中第二类岩心核磁共振T2分布图;
图14为本发明的具体应用实例中第三类岩心铸体薄片;
图15为本发明的具体应用实例中第三类岩心T2弛豫谱示意图;
图16为本发明的具体应用实例中第三类岩心弛豫时间直方图;
图17-图24为本发明的具体应用实例三类岩心的毛管压力与核磁共振弛豫时间T2拟合结果示意图;
图25为本发明的具体应用实例中第一类岩心毛管压力计算模型校验结果示意图;
图26为本发明的具体应用实例中第二类岩心毛管压力计算模型校验结果示意图;
图27为本发明的具体应用实例中第三类岩心毛管压力计算模型校验结果示意图;
图28为本发明的具体应用实例中第一类岩心最终毛管压力计算模型校验结果示意图;
图29为本发明的具体应用实例中第二类岩心最终毛管压力计算模型校验结果示意图;
图30为本发明的具体应用实例中第三类岩心最终毛管压力计算模型校验结果示意图;
图31为本发明的具体应用实例中三类岩心最终毛管压力计算模型针对狮40井13号岩心校验结果示意图;
图32为本发明的具体应用实例中三类岩心最终毛管压力计算模型针对鸭K1-47井1号岩心校验结果示意图;
图33为本发明的实施例中的基于核磁测试的毛管压力确定装置的结构示意图;
图34为本发明的实施例中的电子设备的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
鉴于现有技术中针对毛管压力的测量方法存在着精度差,且不能应用于岩心较少的勘探地区的问题,本发明的实施例提供一种基于核磁测试的毛管压力确定系统,参见图1,该装置可以为一种服务器A1,该服务器A1可以与多个测试设备端B1通信连接,服务器A1还可以与多个数据库分别通信连接,或者如图2所示,这些数据库也可以之间设置在服务器A1中。其中,用户可以通过登录测试设备端B1的网站或者APP等,与服务器A1建立连接。服务器A1根据数据库中存储的内容获取多个岩心对应的测试指令,其中,测试指令中包含有:核磁测试指令、毛管压力测试指令及铸体薄片测试指令,根据各个测试设备标示,将对应的各个测试指令分别发送至各个测试设备端B1。
在实际应用中,进行基于核磁测试的毛管压力确定的部分可以在如上述内容的服务器A1侧执行,即,如图1或图2所示的架构。
上述的客户端设备可以具有通信模块(即通信单元),可以与远程的服务器进行通信连接,实现与服务器的数据传输。服务器可以包括任务调度中心一侧的服务器,其他的实施场景中也可以包括中间平台的服务器,例如与任务调度中心服务器有通信链接的第三方服务器平台的服务器。的服务器可以包括单台计算机设备,也可以包括多个服务器组成的服务器集群,或者分布式装置的服务器结构。
服务器与客户端设备之间可以使用任何合适的网络协议进行通信,包括在本申请提交日尚未开发出的网络协议。网络协议例如可以包括TCP/IP协议、UDP/IP协议、HTTP协议、HTTPS协议等。当然,网络协议例如还可以包括在上述协议之上使用的RPC协议(RemoteProcedure Call Protocol,远程过程调用协议)、REST协议(Representational StateTransfer,表述性状态转移协议)等。
基于上述内容,本发明的实施例还提供一种基于核磁测试的毛管压力确定方法的具体实施方式,参见图3该方法10具体包括如下内容:
步骤100:测量待计算毛管压力地区的核磁测试数据。
核磁共振测井是一种适用于裸眼井的测井新技术,是目前唯一可以直接测量任意岩性储集层自由流体(油、气、水)渗流体积特性的测井方法。核磁共振技术是利用原子核的顺磁性以及与它们相互作用的外加磁场来测量地层中的氢核在磁场中的自由旋进特性,以获取岩石内自由流体含量的测井方法。可以理解的是,步骤100中待计算毛管压力地区是指于岩心所在地区为同一区块,但没有岩心数据或者岩心数据很少的地区。
步骤200:利用待计算毛管压力地区的核磁测试数据及预先建立的多个不同类别岩心的毛管压力计算模型,生成待计算毛管压力地区的毛管压力。
将待计算毛管压力地区的核磁测试数据输入到预设的毛管压力计算模型中,以获取待计算毛管压力地区的毛管压力。需要注意的是,步骤200中的毛管压力计算模型含有多个,其数量具体依据岩心的孔隙结构来划分。
从上述描述可知,本发明提供基于核磁测试的毛管压力确定方法,首先在岩心分类的基础上,对不同孔隙类型的典型岩心分别测试其核磁共振和毛管压力,将核磁测试数据与毛管压力进行拟合,进而确定不同孔隙类型的核磁测试数据与毛管压力转换公式的变量系数;最终确定不同孔隙类型的毛管压力计算模型。该方法主要是针对岩心取心有限,或者没有取心的地区,无法通过岩心的毛管压力曲线测试结果来分析孔隙结构特征的情况下,直接将岩心或测井核磁共振测试结果通过核磁测试数据与毛管压力转换公式转换为毛管曲线直接来分析孔隙结构特征。综上,本发明可以提供一种精度高、且无需进行岩心毛管压力测试即可获得毛管压力的方法。
一实施例中,参见图4,基于测试的毛管压力确定方法还包括方法20:
步骤S100:根据孔隙结构特征参数,对岩心进行分类。
优选地,步骤S100中孔隙结构特征参数可以为,孔隙类型(例如晶间孔、粒间孔、基质微孔、溶孔及裂缝等)、孔隙度及渗透率。
步骤S200:对分类后的岩心进行核磁测试及毛管压力测试,得到岩心分类后的核磁测试数据及毛管压力测试数据。
需要注意的是,步骤S200中对分类后的岩心进行核磁测试与步骤100中的测量待计算毛管压力地区的核磁测试数据,其测量方法、测量参数以及所需获取的测量结果是一致的,只不过对象及测量环境不同。步骤100是指在井筒环境中对地层进行核磁测试;而步骤S200是指在实验室下对岩心进行核磁测试。
地下岩石的孔隙结构较为复杂,可以看做一系列相互连通的毛细管网络、不亲汞的岩石孔隙(针对汞,岩石孔隙为不润湿),在外加压力作用下,汞克服毛管力可进入岩石孔隙。随压力增加,汞一次由大到小进入岩石孔隙,岩心中的汞饱和度不断增加,注入压力与岩心中汞饱和度的关系即为毛管压力曲线。
步骤S300:根据岩心分类后的核磁测试数据及毛管压力测试数据,利用核磁测试数据与毛管压力转换公式,生成不同类别岩心的毛管压力计算模型。
现有技术在利用核磁测试数据与毛管压力转换公式时,对不同类型的储层采用一种方法进行转换,或对大孔隙和小孔隙分段采用不同的公式进行转换,但是实际岩心的孔隙类型相对实验岩心类型丰富,因此本方法在考虑转换公式的建立时,在对岩心进行详细的孔隙分类的基础上,对不同类型建立不同的转换方法。
一实施例中,参见图5步骤S300包括:
步骤S301:将岩心分类后的核磁测试数据及毛管压力测试数据进行拟合,获得核磁测试数据与毛管压力转换公式的变量系数。
岩心分类后的核磁测试数据及毛管压力测试数据进行拟合,优选的,可以使用多元非线性拟合,进而确定不同类型孔隙发育的岩心的核磁测试数据与毛管压力转换公式的变量系数;以确定不同孔隙结构类型的岩心的核磁测试数据与毛管压力转换公式。
步骤S302:根据变量系数及核磁测试数据与毛管压力转换公式,生成岩心的毛管压力计算模型。
将每类岩心的由步骤S302所得到的变量系数带入到通用的核磁测试数据与毛管压力转换公式中,以获得每类岩心各自的岩心的毛管压力计算模型。
一实施例中,孔隙结构特征参数包括:孔隙类型、孔隙度及渗透率。
一实施例中,核磁测试数据包括:核磁共振弛豫时间。
为进一步地说明本方案,本发明以某油田区块H为例,提供基于核磁测试的毛管压力确定方法的具体应用实例,该具体应用实例具体包括如下内容60,参见图6
S0:对区块H的岩心进行铸体薄片、毛管压力及核磁共振测试。
分别选取不同类型孔隙发育岩心,进行铸体薄片、毛管压力及核磁共振测试来研究其孔隙特征及孔隙结构。
S1:根据孔隙类型、孔隙度及渗透率参数,对区块H的岩心进行分类。
依据S0中测试结果,将岩心进行如下分类,第一类岩心为湖相碳酸盐岩(YH106,3294.2m,E3 2,泥晶白云岩,图7-图9),发育白云石晶间孔,孔隙尺度以纳米-微米级为主,孔隙度15.5%,渗透率0.31mD;第二类岩心为湖相碳酸盐岩(Shi40,3146.87m,E3 2,泥晶灰云岩,图10-图13),该岩心同时发育基质微孔、溶孔及裂缝,从压汞孔喉分布及核磁共振T2弛豫谱均显示其发育三类孔隙,孔隙度7.59%,渗透率0.6mD;第三类岩心为碎屑岩岩心(鸭K1-47,2805.58m,N1,细砂岩,图14及图16),主要发育粒间孔,孔隙度16.13%,渗透率6.73mD。
可以理解的是,可以在步骤S0及S1中,可以先对岩心进行铸体薄片、毛管压力及核磁共振测试,之后再将岩心进行分类,依据岩心分类结果即可获取岩心测试结果分类;但同样可以先进行岩心分类,然后对每类岩心进行铸体薄片、毛管压力及核磁共振测试,同样可以得到每类岩心的测试结果,本发明不以此为限制。
S2:将S1中三类岩心的测试结果分别带入多个经验公式中,以获取三类岩心各自的毛管压力计算模型。
本具体应用实例中核磁共振测试结果选取核磁共振弛豫时间(T2谱),即将核磁共振弛豫时间(T2)分别带入下述经验公式通式中,并对两变量进行拟合,通过拟合公式后的结果确定经验公式中的常数(如经验公式①常数m,n;和经验公式②常数a,b,c,d,k)。优选地,可以采用数据处理软件Grapher进行拟合。参见图17-图24。
经验公式①:
经验公式②:
S3:毛管压力计算模型优选。
步骤S3可以理解为岩心测试数据转换效果校验,确定不同类型孔隙适用的经验公式,并根据检验效果对经验公式做适当修改,确定不同类型孔隙适用的经验公式;
具体实施时,将不同类型孔隙的核磁弛豫时间带入图11-图18中所示的具体经验公式中,计算毛管压力,并对计算出毛管压力曲线与实测毛管曲线的符合程度来确定具体岩心应采用的经验公式,通过转换毛管曲线与实测毛管压力曲线对比发现基质微孔发育岩心,经验公式①转换毛管曲线符合程度较好(图25);裂缝发育岩心经验公式①和经验公式②的转换效果均不理想(图26);粒间孔发育岩心经验公式②转换毛管曲线与实测毛管压力曲线符合程度较高(图27)。
S4:毛管压力计算模型优化。
对比不同类型孔隙岩心的孔隙度和渗透率分布后发现,孔隙度和渗透率能较好的区分三类岩心,故在经验公式添加反映孔隙度和渗透率特征的计算项,来改善裂缝类储层的转换毛管曲线与实测毛管曲线的符合程度,通过拟合计算确定添加渗透率和孔隙度乘积的平方根做为计算项添加到经验公式①中有效的改善了转换毛管曲线与实测毛管曲线的符合程度(图28-图30)。
通过以上拟合计算确定了核磁共振T2弛豫时间转换为毛管压力曲线的复合公式:
符合公式中计算项0.06×T2 1.56代表基质微孔转换公式,P1为该类孔隙对应权重;计算项(KQ)1/2×0.01×T2 1.61为裂缝转换公式,P2为对应权重;
为粒间孔转换公式,P3为对应权重。
S5:最终毛管压力计算模型校验。
分别选取了基质微孔发育的湖相碳酸盐岩岩心和粒间孔发育的细砂岩岩心,代入公式(11)进行效果验证,发现对于基质微孔发育岩心小孔喉与基质微孔公式及裂缝公式转换效果对应较好,大孔喉出转换效果较差;对粒间孔发育岩心小孔喉与基质微孔公式转换效果较好,大孔喉处与粒间孔转换公式对应较好(图31-图32)。
S6:最终毛管压力计算模型应用。
可以理解的是,步骤S5确定了最终的毛管压力计算模型之后,可将待计算毛管压力地区(无岩心或岩心很少的地区)的核磁测试数据带入该模型,以获取待计算毛管压力地区的毛管压力。
S7:生成毛管压力曲线。
由步骤S6中的多个毛管压力值,生成待计算毛管压力地区的毛管压力曲线,以分析待计算毛管压力地区,即无岩心或岩心很少的地区的储层物性特征。
从上述描述可知,本发明提供基于核磁测试的毛管压力确定方法,首先在岩心分类的基础上,对不同孔隙类型的典型岩心分别测试其核磁共振和毛管压力,将核磁测试数据与毛管压力进行拟合,进而确定不同孔隙类型的核磁测试数据与毛管压力转换公式的变量系数;最终确定不同孔隙类型的毛管压力计算模型。该方法主要是针对岩心取心有限,或者没有取心的地区,无法通过岩心的毛管压力曲线测试结果来分析孔隙结构特征的情况下,直接将岩心或测井核磁共振测试结果通过核磁测试数据与毛管压力转换公式转换为毛管曲线直接来分析孔隙结构特征。综上,本发明可以提供一种精度高、且无需进行岩心毛管压力测试即可获得毛管压力的方法。
基于同一发明构思,本申请实施例还提供了基于核磁测试的毛管压力确定装置,可以用于实现上述实施例所描述的方法,如下面的实施例。由于基于核磁测试的毛管压力确定装置解决问题的原理与基于核磁测试的毛管压力确定方法相似,因此基于核磁测试的毛管压力确定装置的实施可以参见基于核磁测试的毛管压力确定方法实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的系统较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
本发明的实施例提供一种能够实现基于核磁测试的毛管压力确定方法的基于核磁测试的毛管压力确定装置的具体实施方式,参见图33,基于核磁测试的毛管压力确定装置具体包括如下内容:
核磁测试数据测量单元10,用于测量待计算毛管压力地区的核磁测试数据;
毛管压力生成单元20,用于利用待计算毛管压力地区的核磁测试数据及预先建立的多个不同孔隙结构类别岩心的毛管压力计算模型,生成待计算毛管压力地区的毛管压力。
优选地,毛管压力确定装置,还包括:
岩心分类单元,用于根据孔隙结构特征参数,对岩心进行分类;
测试单元,用于对分类后的岩心进行核磁测试及毛管压力测试,得到岩心分类后的核磁测试数据及毛管压力测试数据;
毛管压力计算模型生成单元,用于根据岩心分类后的核磁测试数据及毛管压力测试数据,利用核磁测试数据与毛管压力转换公式,生成不同类别岩心的毛管压力计算模型。
优选地,毛管压力计算模型生成单元包括:
变量系数获得模块,用于将岩心分类后的核磁测试数据及毛管压力测试数据进行拟合,获得核磁测试数据与毛管压力转换公式的变量系数;
毛管压力计算模型模块,用于根据变量系数及核磁测试数据与毛管压力转换公式,生成岩心的毛管压力计算模型。
从上述描述可知,本发明提供基于核磁测试的毛管压力确定装置,首先在岩心分类的基础上,对不同孔隙类型的典型岩心分别测试其核磁共振和毛管压力,将核磁测试数据与毛管压力进行拟合,进而确定不同孔隙类型的核磁测试数据与毛管压力转换公式的变量系数;最终确定不同孔隙类型的毛管压力计算模型。该方法主要是针对岩心取心有限,或者没有取心的地区,无法通过岩心的毛管压力曲线测试结果来分析孔隙结构特征的情况下,直接将岩心或测井核磁共振测试结果通过核磁测试数据与毛管压力转换公式转换为毛管曲线直接来分析孔隙结构特征。综上,本发明可以提供一种精度高、且无需进行岩心毛管压力测试即可获得毛管压力的方法。
本申请的实施例还提供能够实现上述实施例中的基于核磁测试的毛管压力确定方法中全部步骤的一种电子设备的具体实施方式,参见图34,电子设备具体包括如下内容:
处理器(processor)1201、存储器(memory)1202、通信接口(CommunicationsInterface)1203和总线1204;
其中,处理器1201、存储器1202、通信接口1203通过总线1204完成相互间的通信;通信接口1203用于实现服务器端设备、测试设备以及用户端设备等相关设备之间的信息传输。
处理器1201用于调用存储器1202中的计算机程序,处理器执行计算机程序时实现上述实施例中的基于核磁测试的毛管压力确定方法中的全部步骤,例如,处理器执行计算机程序时实现下述步骤:
步骤100:测量待计算毛管压力地区的核磁测试数据。
步骤200:利用待计算毛管压力地区的核磁测试数据及预先建立的多个不同类别岩心的毛管压力计算模型,生成待计算毛管压力地区的毛管压力。
从上述描述可知,本申请实施例中的电子设备,首先在岩心分类的基础上,对不同孔隙类型的典型岩心分别测试其核磁共振和毛管压力,将核磁测试数据与毛管压力进行拟合,进而确定不同孔隙类型的核磁测试数据与毛管压力转换公式的变量系数;最终确定不同孔隙类型的毛管压力计算模型。该方法主要是针对岩心取心有限,或者没有取心的地区,无法通过岩心的毛管压力曲线测试结果来分析孔隙结构特征的情况下,直接将岩心或测井核磁共振测试结果通过核磁测试数据与毛管压力转换公式转换为毛管曲线直接来分析孔隙结构特征。综上,本发明可以提供一种精度高、且无需进行岩心毛管压力测试即可获得毛管压力的方法。
本申请的实施例还提供能够实现上述实施例中的基于核磁测试的毛管压力确定方法中全部步骤的一种计算机可读存储介质,计算机可读存储介质上存储有计算机程序,该计算机程序被处理器执行时实现上述实施例中的基于核磁测试的毛管压力确定方法的全部步骤,例如,处理器执行计算机程序时实现下述步骤:
步骤100:测量待计算毛管压力地区的核磁测试数据。
步骤200:利用待计算毛管压力地区的核磁测试数据及预先建立的多个不同类别岩心的毛管压力计算模型,生成待计算毛管压力地区的毛管压力。
从上述描述可知,本申请实施例中的计算机可读存储介质,首先在岩心分类的基础上,对不同孔隙类型的典型岩心分别测试其核磁共振和毛管压力,将核磁测试数据与毛管压力进行拟合,进而确定不同孔隙类型的核磁测试数据与毛管压力转换公式的变量系数;最终确定不同孔隙类型的毛管压力计算模型。该方法主要是针对岩心取心有限,或者没有取心的地区,无法通过岩心的毛管压力曲线测试结果来分析孔隙结构特征的情况下,直接将岩心或测井核磁共振测试结果通过核磁测试数据与毛管压力转换公式转换为毛管曲线直接来分析孔隙结构特征。综上,本发明可以提供一种精度高、且无需进行岩心毛管压力测试即可获得毛管压力的方法。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于硬件+程序类实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下,在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外,在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中,多任务处理和并行处理也是可以的或者可能是有利的。
虽然本申请提供了如实施例或流程图的方法操作步骤,但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式,不代表唯一的执行顺序。在实际中的装置或客户端产品执行时,可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。
虽然本说明书实施例提供了如实施例或流程图所述的方法操作步骤,但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式,不代表唯一的执行顺序。在实际中的装置或终端产品执行时,可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境,甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下,并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。
为了描述的方便,描述以上装置时以功能分为各种模块分别描述。当然,在实施本说明书实施例时可以把各模块的功能在同一个或多个软件和/或硬件中实现,也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
本领域技术人员也知道,除了以纯计算机可读程序代码方式实现控制器以外,完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件,而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。
本领域技术人员应明白,本说明书的实施例可提供为方法、系统或计算机程序产品。因此,本说明书实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且,本说明书实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本说明书实施例可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本说明书实施例,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书实施例的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
以上所述仅为本说明书实施例的实施例而已,并不用于限制本说明书实施例。对于本领域技术人员来说,本说明书实施例可以有各种更改和变化。凡在本说明书实施例的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本说明书实施例的权利要求范围之内。
Claims (8)
1.一种基于核磁测试的毛管压力确定方法,其特征在于,包括:
测量待计算毛管压力地区的核磁测试数据;
利用所述待计算毛管压力地区的核磁测试数据及预先建立的多个不同孔隙结构类别岩心的毛管压力计算模型,生成所述待计算毛管压力地区的毛管压力,包括:
根据孔隙结构特征参数,对所述岩心进行分类;
对分类后的岩心进行核磁测试及毛管压力测试,得到岩心分类后的核磁测试数据及毛管压力测试数据;
根据所述岩心分类后的核磁测试数据及毛管压力测试数据,利用核磁测试数据与毛管压力转换公式,生成不同类别岩心的毛管压力计算模型;
基于每一个所述类别岩心的毛管压力计算模型及所述岩心分类在待计算毛管压力地区的权重,建立复合毛管压力计算模型;
利用所述待计算毛管压力地区的核磁测试数据及所述复合毛管压力计算模型生成所述待计算毛管压力地区的毛管压力;
所述毛管压力计算模型为:
计算项0.06×T2 1.56代表基质微孔转换公式,P1为该类孔隙对应权重;计算项(KQ)1/2×0.01×T2 1.61为裂缝转换公式,P2为对应权重;
为粒间孔转换公式,P3为对应权重。
2.如权利要求1所述的毛管压力确定方法,其特征在于,所述根据所述岩心分类后的核磁测试数据及毛管压力测试数据,利用核磁测试数据与毛管压力转换公式,生成不同类别岩心的毛管压力计算模型,包括:
将所述岩心分类后的核磁测试数据及毛管压力测试数据进行拟合,获得所述核磁测试数据与毛管压力转换公式的变量系数;
根据所述变量系数及所述核磁测试数据与毛管压力转换公式,生成所述岩心的毛管压力计算模型。
3.如权利要求1所述的毛管压力确定方法,其特征在于,所述孔隙结构特征参数包括:孔隙类型、孔隙度及渗透率。
4.如权利要求1所述的毛管压力确定方法,其特征在于,所述核磁测试数据包括:核磁共振弛豫时间。
5.一种基于核磁测试的毛管压力确定装置,其特征在于,包括:
核磁测试数据测量单元,用于测量待计算毛管压力地区的核磁测试数据;
毛管压力生成单元,用于利用所述待计算毛管压力地区的核磁测试数据及预先建立的多个不同孔隙结构类别岩心的毛管压力计算模型,生成所述待计算毛管压力地区的毛管压力,所述毛管压力生成单元包括:
岩心分类单元,用于根据孔隙结构特征参数,对所述岩心进行分类;
测试单元,用于对分类后的岩心进行核磁测试及毛管压力测试,得到岩心分类后的核磁测试数据及毛管压力测试数据;
毛管压力计算模型生成单元,用于根据所述岩心分类后的核磁测试数据及毛管压力测试数据,利用核磁测试数据与毛管压力转换公式,生成不同类别岩心的毛管压力计算模型;
基于每一个所述类别岩心的毛管压力计算模型及所述岩心分类在待计算毛管压力地区的权重,建立复合毛管压力计算模型;
利用所述待计算毛管压力地区的核磁测试数据及所述复合毛管压力计算模型生成所述待计算毛管压力地区的毛管压力;
所述毛管压力计算模型为:
计算项0.06×T2 1.56代表基质微孔转换公式,P1为该类孔隙对应权重;计算项(KQ)1/2×0.01×T2 1.61为裂缝转换公式,P2为对应权重;
为粒间孔转换公式,P3为对应权重。
6.如权利要求5所述的毛管压力确定装置,其特征在于,所述毛管压力计算模型生成单元包括:
变量系数获得模块,用于将所述岩心分类后的核磁测试数据及毛管压力测试数据进行拟合,获得所述核磁测试数据与毛管压力转换公式的变量系数;
毛管压力计算模型模块,用于根据所述变量系数及所述核磁测试数据与毛管压力转换公式,生成所述岩心的毛管压力计算模型。
7.一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现权利要求1至4任一项所述基于核磁测试的毛管压力确定方法的步骤。
8.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该计算机程序被处理器执行时实现权利要求1至4任一项所述基于核磁测试的毛管压力确定方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910444319.6A CN111999331B (zh) | 2019-05-27 | 2019-05-27 | 基于核磁测试的毛管压力确定方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910444319.6A CN111999331B (zh) | 2019-05-27 | 2019-05-27 | 基于核磁测试的毛管压力确定方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111999331A CN111999331A (zh) | 2020-11-27 |
CN111999331B true CN111999331B (zh) | 2023-10-31 |
Family
ID=73461321
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910444319.6A Active CN111999331B (zh) | 2019-05-27 | 2019-05-27 | 基于核磁测试的毛管压力确定方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111999331B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113075102B (zh) * | 2021-03-11 | 2022-02-01 | 中国地质大学(北京) | 一种建立多孔介质自发渗吸量与时间关系数学模型的方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105353419A (zh) * | 2015-11-25 | 2016-02-24 | 中国地质大学(北京) | 一种基于储集层分类的核磁毛管压力曲线构造方法 |
CN105386753A (zh) * | 2015-10-28 | 2016-03-09 | 中国地质大学(北京) | 一种利用核磁共振测井构造伪毛管压力曲线方法 |
-
2019
- 2019-05-27 CN CN201910444319.6A patent/CN111999331B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105386753A (zh) * | 2015-10-28 | 2016-03-09 | 中国地质大学(北京) | 一种利用核磁共振测井构造伪毛管压力曲线方法 |
CN105353419A (zh) * | 2015-11-25 | 2016-02-24 | 中国地质大学(北京) | 一种基于储集层分类的核磁毛管压力曲线构造方法 |
Non-Patent Citations (3)
Title |
---|
云冈石窟不同类型砂岩的核磁共振T_2谱――压汞毛管压力换算C值研究;周华;高峰;周萧;刘化冰;郭葆鑫;王昌燧;;地球物理学进展(第05期);2759-2766 * |
致密砂岩储层核磁测井资料计算毛管压力曲线方法研究;朱建华;刘跃辉;;国外测井技术(第04期);17-19 * |
酒西盆地间泉子段储层流体赋存及渗流特征;张世铭;王建功;张小军;张婷静;曹志强;杨麟科;;天然气地球科学(第08期);11-15 * |
Also Published As
Publication number | Publication date |
---|---|
CN111999331A (zh) | 2020-11-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105866009B (zh) | 一种计算致密油储层有效孔隙度的方法及装置 | |
US9140117B2 (en) | Method for evaluating relative permeability for fractional multi-phase, multi-component fluid flow through porous media | |
CN110927035A (zh) | 一种低渗致密砂岩束缚水饱和度计算方法 | |
CN105891089B (zh) | 储层渗透率确定方法及装置 | |
US9835762B2 (en) | Petrophysical rock characterization | |
CN103884633A (zh) | 一种确定岩石渗透率的方法及装置 | |
US10976274B2 (en) | Measuring the wettability of porous media based on the temperature sensitivity of nuclear magnetic resonance relaxation time | |
Zapata et al. | Modeling adsorption–desorption hysteresis in shales: Acyclic pore model | |
CN113825996B (zh) | 用于在脉冲衰减实验中确定岩心渗透率的方法和系统 | |
CN109932297A (zh) | 一种致密砂岩储层渗透率的计算方法 | |
CN106355571A (zh) | 一种白云岩储层质量的确定方法及装置 | |
CN110702580A (zh) | 一种基于信息熵的致密砂岩储层孔喉非均质性表征方法 | |
CN108268712B (zh) | 核磁共振确定孔隙介质毛管压力的方法和装置 | |
CN111999331B (zh) | 基于核磁测试的毛管压力确定方法及装置 | |
CN105844011B (zh) | 一种基于毛管模型的渗透率计算方法 | |
Cui et al. | Apparent permeability and representative size of shale: a numerical study on the effects of organic matter | |
Zan et al. | Experimental investigation of spontaneous water imbibition into methane-saturated shales under different methane pressures | |
CN108681793B (zh) | 深层油藏采油指数预测方法及装置 | |
Ji et al. | Rapid evaluation of capillary pressure and relative permeability for oil–water flow in tight sandstone based on a physics-informed neural network | |
CN116738794A (zh) | 孔裂隙介质的岩石物理数值模拟方法、装置、设备及介质 | |
CN112049623B (zh) | 一种中基性火山岩束缚水出水的判定方法和装置 | |
CN113945497A (zh) | 油气藏储层流体的可动性的评价方法 | |
CN103675945B (zh) | 一种测定孔洞型储层的饱和度的方法及设备 | |
CN112464587B (zh) | 一种预测储层毛管压力曲线的方法、装置及设备 | |
Babamahmoudi et al. | Absolute permeability assessment of porous structures under different boundary conditions using lattice Boltzmann method |
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 |