CN109736793B - 一种评价储层孔隙结构的方法及设备 - Google Patents
一种评价储层孔隙结构的方法及设备 Download PDFInfo
- Publication number
- CN109736793B CN109736793B CN201910047194.3A CN201910047194A CN109736793B CN 109736793 B CN109736793 B CN 109736793B CN 201910047194 A CN201910047194 A CN 201910047194A CN 109736793 B CN109736793 B CN 109736793B
- Authority
- CN
- China
- Prior art keywords
- magnetic resonance
- nuclear magnetic
- echo data
- resonance echo
- reservoir
- 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
Images
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明实施例提供一种评价储层孔隙结构的方法及设备,该方法通过获取目标储层的核磁共振回波数据,对核磁共振回波数据进行时域变换,并计算变换后的核磁共振回波数据的振幅的平方及频率,对振幅的平方和频率进行线性拟合,得到拟合直线,获取拟合直线对应的斜率,并根据斜率计算核磁共振回波数据的分形维数,根据分形维数评价目标储层的孔隙结构,自动根据核磁共振回波数据进行分形分析,得到相应的分形维数,无需进行核磁共振回波数据反演,计算量较小,减少计算误差。
Description
技术领域
本发明实施例涉及油气田勘探开发技术领域,尤其涉及一种评价储层孔隙结构的方法及设备。
背景技术
在油气开发技术领域,岩石的孔隙结构是指岩石内的孔隙和喉道类型、大小、分布及其相互连通关系。研究油气储层的岩石孔隙结构,深入揭示油气储层的内部结构,对油气田勘探和开发有着重要的意义。
目前,为了评价油气储层的孔隙结构,获取油气储层的孔隙结构参数,一般采用的是核磁共振测井。当需要评价储层孔隙结构时,通常对核磁共振测井采集到的核磁共振回波数据进行反演得到核磁共振T2分布,然后计算核磁共振T2分布的分形维数来评价储层孔隙结构。
然而,发明人发现现有技术中至少存在如下问题:对于比常规储层更加复杂的致密砂岩储层,进行核磁共振测井时,在低信噪比下核磁共振回波数据的反演会存在较大的误差,导致计算出来的分形维数误差较大。
发明内容
本发明实施例提供一种评价储层孔隙结构的方法及设备,以减少计算分形维数的计算误差。
第一方面,本发明实施例提供一种评价储层孔隙结构的方法,包括:
获取目标储层的核磁共振回波数据;
对所述核磁共振回波数据进行时域变换,并计算变换后的核磁共振回波数据的振幅的平方及频率;
对所述振幅的平方和频率进行线性拟合,得到拟合直线;
获取所述拟合直线对应的斜率,并根据所述斜率计算所述核磁共振回波数据的分形维数;
根据所述分形维数评价所述目标储层的孔隙结构。
在一种可能的设计中,所述根据所述斜率计算所述核磁共振回波数据的分形维数的计算公式为:
其中,Dechoe为所述分形维数,k为拟合直线斜率的绝对值。
在一种可能的设计中,所述对所述核磁共振回波数据进行时域变换之前,还包括:
判断目标储层核磁共振回波数据对应的概率分布是否符合高斯分布,其中目标储层核磁共振回波数据为所述核磁共振回波数据中的任一个;
若所述目标储层核磁共振回波数据的概率分布不符合高斯分布,则对所述目标储层核磁共振回波数据进行正态变换。
在一种可能的设计中,所述对所述核磁共振回波数据进行时域变换,包括:
按照预设的傅立叶变换公式,对所述核磁共振回波数据进行傅立叶变换;
在一种可能的设计中,所述计算变换后的核磁共振回波数据的振幅的平方及频率,包括:
获取核磁共振回波数据的回波个数以及所述傅立叶变换公式中的实部和虚部;
根据所述核磁共振回波个数和所述傅立叶变换公式中的实部与虚部,计算所述变换后的核磁共振回波数据的振幅的平方;
获取所述核磁共振回波数据的回波间隔;
根据所述回波间隔和所述核磁共振回波个数计算所述核磁共振回波数据中的采样点的频率。
在一种可能的设计中,所述根据所述核磁共振回波个数和所述傅立叶变换公式中的实部与虚部计算所述变换后的核磁共振回波数据的振幅的平方的计算公式为:
其中,所述a为所述傅立叶变换公式中的实部,b为所述傅立叶变换公式中的虚部,N1为所述核磁共振回波个数,P(ω)为所述振幅的平方。
在一种可能的设计中,所述根据所述回波间隔和所述核磁共振回波个数计算所述核磁共振回波数据中的采样点的频率的计算公式为:
其中,ω为第n个采样点的频率,N1为所述核磁共振回波个数,TE为所述回波间隔。
在一种可能的设计中,所述对所述振幅的平方和频率进行线性拟合,包括:
在预设的双对数坐标系中,对所述振幅的平方和频率进行线性拟合。
第二方面,本发明实施例提供一种评价储层孔隙结构的设备,包括:
回波数据获取模块,用于获取目标储层的核磁共振回波数据;
回波数据变换模块,用于对所述核磁共振回波数据进行时域变换,并计算变换后的核磁共振回波数据的振幅的平方及频率;
线性拟合模块,用于对所述振幅的平方和频率进行线性拟合,得到拟合直线;
分形维数计算模块,用于获取所述拟合直线对应的斜率,并根据所述斜率计算所述核磁共振回波数据的分形维数;
孔隙结构评价模块,用于根据所述分形维数评价所述目标储层的孔隙结构。
在一种可能的设计中,在所述分形维数计算模块中,根据所述斜率计算所述核磁共振回波数据的分形维数的计算公式为:
其中,Dechoe为所述分形维数,k为拟合直线斜率的绝对值。
在一种可能的设计中,所述评价储层孔隙结构的设备还包括:正态变换模块;
所述正态变换模块,用于在对所述核磁共振回波数据进行时域变换之前,判断目标储层核磁共振回波数据对应的概率分布是否符合高斯分布,其中目标储层核磁共振回波数据为所述核磁共振回波数据中的任一个;
若所述目标储层核磁共振回波数据的概率分布不符合高斯分布,则对所述目标储层核磁共振回波数据进行正态变换。
在一种可能的设计中,所述回波数据变换模块包括:傅立叶变换单元;
所述傅立叶变换单元,用于按照预设的傅立叶变换公式,对所述核磁共振回波数据进行傅立叶变换。
在一种可能的设计中,所述回波数据变换模块包括:
变换数据获取单元,用于获取核磁共振回波数据的核磁共振回波个数以及所述傅立叶变换公式中的实部和虚部;
振幅平方计算单元,用于根据所述核磁共振回波个数和所述傅立叶变换公式中的实部与虚部,计算所述变换后的核磁共振回波数据的振幅的平方;
回波间隔获取单元,用于获取所述核磁共振回波数据的回波间隔;
频率计算单元,用于根据所述回波间隔和所述核磁共振回波个数计算所述核磁共振回波数据中的采样点的频率。
在一种可能的设计中,在所述振幅平方计算单元中,根据所述核磁共振回波个数和所述傅立叶变换公式中的实部与虚部计算所述变换后的核磁共振回波数据的振幅的平方的计算公式为:
其中,所述a为所述傅立叶变换公式中的实部,b为所述傅立叶变换公式中的虚部,N1为所述核磁共振回波个数,P(ω)为所述振幅的平方。
在一种可能的设计中,在所述频率计算单元中,根据所述回波间隔和所述核磁共振回波个数计算所述核磁共振回波数据中的采样点的频率的计算公式为:
其中,ω为第n个采样点的频率,N1为所述核磁共振回波个数,TE为所述回波间隔。
在一种可能的设计中,所述线性拟合模块具体用于:
在预设的双对数坐标系中,对所述振幅的平方和频率进行线性拟合。
第三方面,本发明实施例提供一种评价储层孔隙结构的设备,包括:至少一个处理器和存储器;
所述存储器存储计算机执行指令;
所述至少一个处理器执行所述存储器存储的计算机执行指令,使得所述至少一个处理器执行如第一方面任一项所述的评价储层孔隙结构的方法。
第四方面,本发明实施例提供一种计算机可读存储介质,所述计算机可读存储介质中存储有计算机执行指令,当处理器执行所述计算机执行指令时,实现如第一方面任一项所述的评价储层孔隙结构的方法。
本发明实施例提供的评价储层孔隙结构的方法及设备,该方法通过获取目标储层的核磁共振回波数据,对核磁共振回波数据进行时域变换,并计算变换后的核磁共振回波数据的振幅的平方及频率,对振幅的平方和频率进行线性拟合,得到拟合直线,获取拟合直线对应的斜率,并根据斜率计算核磁共振回波数据的分形维数,根据分形维数评价目标储层的孔隙结构,自动根据核磁共振回波数据进行分形分析,得到相应的分形维数,无需进行核磁共振回波数据反演,计算量较小,减少计算误差。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的评价储层孔隙结构的场景示意图;
图2为本发明实施例提供的评价储层孔隙结构的方法的流程图一;
图3为本发明实施例提供的线性拟合示意图;
图4为本发明实施例提供的分形维数与储层物性参数及孔隙结构参数的交会图的示意图;
图5为本发明实施例提供的评价储层孔隙结构的方法的流程图二;
图6为本发明实施例提供的评价储层孔隙结构的设备的结构示意图一;
图7为本发明实施例提供的评价储层孔隙结构的设备的结构示意图二;
图8为本发明实施例提供的评价储层孔隙结构的设备的硬件结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”、“第三”“第四”等(如果存在)是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
图1为本发明实施例提供的评价储层孔隙结构的场景示意图,如图1所示,评价储层孔隙结构包括评价储层孔隙结构的设备101和核磁共振测井102。
将核磁共振井102放至目标储层,每隔一定时间,核磁共振测井采集目标储层的核磁共振回波信号,并将采集到的核磁共振回波信号发送至评价储层孔隙结构的设备。评价储层孔隙结构的设备101接收到多段回波信号。评价储层孔隙结构的设备101基于回波信号进行分形分析,得到对应的分形维数。储层的分形维数与孔隙结构之间存在一定的关系,根据核磁共振回波信号对应的分形维数可以得到目标储层的孔隙结构。
下面以具体地实施例对本发明的技术方案进行详细说明。下面这几个具体的实施例可以相互结合,对于相同或相似的概念或过程可能在某些实施例不再赘述。
图2为本发明实施例提供的评价储层孔隙结构的方法的流程图一。本实施例的方法应用于图1所示实施例中的评价储层孔隙结构的设备,如图2所示,该方法包括:
S201:获取目标储层的核磁共振回波数据。
在本实施例中,每隔一定时间,核磁共振测井将采集到的目标储层的核磁共振回波信号实时发送至评价储层孔隙结构的设备。
在本实施例中,核磁共振回波数据是由多段回波信号对应的数据组成,一段回波信号与另一段回波信号之间存在间隔的时间,即存在回波间隔。
在本实施例中,回波信号对应的数据包含接收时刻和该时刻对应的回波幅值。
S202:对核磁共振回波数据进行时域变换,并计算变换后的核磁共振回波数据的振幅的平方及频率。
在本实施例中,将核磁共振回波数据从时域变换到频域,实现时域与频域的转换,并计算变换后的核磁共振回波数据的振幅的平方及频率。
在一种可能的设计中,在对核磁共振回波数据进行时域变换之前,判断目标储层核磁共振回波数据对应的概率分布是否符合高斯分布,其中目标储层核磁共振回波数据为核磁共振回波数据中的任一个,若目标储层核磁共振回波数据的概率分布不符合高斯分布,则对目标储层核磁共振回波数据进行正态变换。
在本实施例中,正态变换可以为BOX-COX变换(Box-Cox transformation)。
在本实施例中,目标储层核磁共振回波数据为核磁共振回波数据中的一段回波信号对应的数据。
在本实施例中,判断目标储层核磁共振回波数据是否符合高斯分布,若该目标储层核磁共振回波数据符合高斯分布,则直接对目标储层核磁共振回波数据进行时域变换。若该目标储层核磁共振回波数据不符合高斯分布,则先对不符合高斯分布的目标储层核磁共振回波数据进行BOX-COX变换,然后对进行BOX-COX变换后的目标储层核磁共振回波数据进行时域变换。
其中,对不符合高斯分布的目标储层核磁共振回波数据进行BOX-COX变换的计算公式是:其中,λ为预设数值,其可以由相关人员根据实际情况选择的合适的值,g(t)为概率分布不满足正态分布的核磁共振回波数据,H(t)为正态变换后的近似正态分布的核磁共振回波数据。
在利用BOX-COX变换时,函数g(t)的各个分量应均大于0,否则应将g(t)做平移后,使得g(t+α)的各个分量都大于0后再进行BOX-COX变换。
在本实施例中,可以根据选取的λ的具体数值,利用对应的BOX-COX变换公式,对不符合高斯分布的目标储层核磁共振回波数据进行BOX-COX变换。
S203:对振幅的平方和频率进行线性拟合,得到拟合直线。
在一种可能的设计中,步骤S203具体包括:在预设的双对数坐标系中,对振幅的平方和频率进行线性拟合。
在本实施例中,在对振幅的平方和频率进行线性拟合之前,还需要将振幅的平方和频率中的高频噪声部分去掉,从而使分形维数的计算结果更加准确。
在本实施例中,如图3所示,振幅的平方的对数与频率的对数之间存在线性关系,即log(P(ω))=-klog(ω)+C'。因此,可以在预设的双对数坐标系中,对振幅的平方和频率进行线性拟合。
S204:获取拟合直线对应的斜率,并根据斜率计算核磁共振回波数据的分形维数。
在一种可能的设计中,根据斜率计算所述核磁共振回波数据的分形维数的计算公式为:
其中,Dechoe为所述分形维数,k为拟合直线斜率的绝对值。
S205:根据分形维数评价目标储层的孔隙结构。
在本实施例中,如图4所示,随着核磁共振回波数据的分形维数(D)增大,T2几何均值(T2lm)减小,储层相对分选系数增大,储层的非均质性增加,储层孔隙结构变差,相应的储层平均孔隙半径减小,渗透率(k)降低,排驱压力(Pd)及中值压力(Pc50)增大。
在本实施例中,通过获取目标储层的核磁共振回波数据,对核磁共振回波数据进行时域变换,并计算变换后的核磁共振回波数据的振幅的平方及频率,对振幅的平方和频率进行线性拟合,得到拟合直线,获取拟合直线对应的斜率,并根据斜率计算核磁共振回波数据的分形维数,根据分形维数评价目标储层的孔隙结构,自动根据核磁共振回波数据进行分形分析,得到相应的分形维数,无需进行核磁共振回波数据反演,计算量较小,减少计算误差。
图5为本发明实施例提供的评价储层孔隙结构的方法的流程示意图二,本实施例在图2实施例的基础上,对步骤S202的具体实现过程进行了详细说明。如图5所示,该方法包括:
S501:按照预设的傅立叶变换公式,对核磁共振回波数据进行傅立叶变换。
在本实施例中,时域变换可以为傅立叶变换,获取核磁共振数据的回波个数,然后基于傅立叶变换公式对核磁共振回波数据进行傅立叶变换,傅立叶变换公式为:
其中,Gk为核磁共振回波数据经时域变换后的频率信号,gn为时域变换前的核磁共振回波数据(时间信号),N1为核磁共振回波个数。
其中,核磁共振回波个数表示核磁共振回波数据包含多少个回波信号。
S502:获取核磁共振回波数据的核磁共振回波个数以及傅立叶变换公式中的实部和虚部。
在本实施例中,Gk=a+bi,a为傅立叶变换公式中的实部,b为傅立叶变换公式中的虚部。相应的,实部和虚部为:
S503:根据核磁共振回波个数和傅立叶变换公式中的实部与虚部,计算变换后的核磁共振回波数据的振幅的平方。
在一种可能的设计中,根据回波间隔和核磁共振回波个数计算核磁共振回波数据中的采样点的频率的计算公式为:
其中,a为傅立叶变换公式中的实部,b为傅立叶变换公式中的虚部,N1为核磁共振回波个数,P(ω)为振幅的平方。
S504:获取核磁共振回波数据的回波间隔。
S505:根据回波间隔和核磁共振回波个数计算核磁共振回波数据中的采样点的频率。
在一种可能的设计中,根据回波间隔和核磁共振回波个数计算核磁共振回波数据中的采样点的频率的计算公式为:
其中,ω为第n个采样点的频率,N1为核磁共振回波个数,TE为回波间隔。
图6为本发明实施例提供的评价储层孔隙结构的设备的结构示意图一,如图6所示,本实施例提供的评价储层孔隙结构的设备600,其可以包括:回波数据获取模块610、回波数据变换模块620、线性拟合模块630、分形维数计算模块640和孔隙结构评价模块650。
其中,回波数据获取模块610,用于获取目标储层的核磁共振回波数据。
回波数据变换模块620,用于对核磁共振回波数据进行时域变换,并计算变换后的核磁共振回波数据的振幅的平方及频率。
线性拟合模块630,用于对振幅的平方和频率进行线性拟合,得到拟合直线。
分形维数计算模块640,用于获取拟合直线对应的斜率,并根据斜率计算核磁共振回波数据的分形维数。
孔隙结构评价模块650,用于根据分形维数评价目标储层的孔隙结构。
在一种可能的设计中,在分形维数计算模块中,根据斜率计算核磁共振回波数据的分形维数的计算公式为:
其中,Dechoe为分形维数,k为拟合直线斜率的绝对值。
在一种可能的设计中,评价储层孔隙结构的设备还包括:正态变换模块。
正态变换模块,用于在对核磁共振回波数据进行时域变换之前,判断目标储层核磁共振回波数据对应的概率分布是否符合高斯分布,其中目标储层核磁共振回波数据为核磁共振回波数据中的任一个。
若目标储层核磁共振回波数据的概率分布不符合高斯分布,则对目标储层核磁共振回波数据进行正态变换。
在一种可能的设计中,线性拟合模块具体用于:
在预设的双对数坐标系中,对振幅的平方和频率进行线性拟合。
本发明实施例提供的评价储层孔隙结构的设备,可以实现上述如图2所示的实施例的评价储层孔隙结构的方法,其实现原理和技术效果类似,此处不再赘述。
图7为本发明实施例提供的评价储层孔隙结构的设备的结构示意图二,如图7所示,图6中的回波数据变换模块620还可以包括:傅立叶变换单元621、变换数据获取单元622、振幅平方计算单元623、回波间隔获取单元624和频率计算单元625。
其中,傅立叶变换单元621,用于按照预设的傅立叶变换公式,对核磁共振回波数据进行傅立叶变换。
变换数据获取单元622,用于获取核磁共振回波数据的核磁共振回波个数以及傅立叶变换公式中的实部和虚部。
振幅平方计算单元623,用于根据核磁共振回波个数和傅立叶变换公式中的实部与虚部,计算变换后的核磁共振回波数据的振幅的平方。
回波间隔获取单元624,用于获取核磁共振回波数据的回波间隔。
频率计算单元625,用于根据回波间隔和核磁共振回波个数计算核磁共振回波数据中的采样点的频率。
在一种可能的设计中,在振幅平方计算单元中,根据核磁共振回波个数和傅立叶变换公式中的实部与虚部计算变换后的核磁共振回波数据的振幅的平方的计算公式为:
其中,a为傅立叶变换公式中的实部,b为傅立叶变换公式中的虚部,N1为核磁共振回波个数,P(ω)为振幅的平方。
在一种可能的设计中,在频率计算单元中,根据回波间隔和核磁共振回波个数计算核磁共振回波数据中的采样点的频率的计算公式为:
其中,ω为第n个采样点的频率,N1为核磁共振回波个数,TE为回波间隔。
本发明实施例提供的评价储层孔隙结构的设备,可以实现上述如图5所示的实施例的评价储层孔隙结构的方法,其实现原理和技术效果类似,此处不再赘述。
图8为本发明实施例提供的评价储层孔隙结构的设备的硬件结构示意图。如图8所示,本实施例提供的客户端800包括:至少一个处理器801和存储器802。其中,处理器801、存储器802通过总线803连接。
在具体实现过程中,至少一个处理器801执行所述存储器802存储的计算机执行指令,使得至少一个处理器801执行上述方法实施例中的文件格式的转换方法。
处理器801的具体实现过程可参见上述方法实施例,其实现原理和技术效果类似,本实施例此处不再赘述。
在上述的图8所示的实施例中,应理解,处理器可以是中央处理单元(英文:Central Processing Unit,简称:CPU),还可以是其他通用处理器、数字信号处理器(英文:Digital Signal Processor,简称:DSP)、专用集成电路(英文:Application SpecificIntegrated Circuit,简称:ASIC)等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合发明所公开的方法的步骤可以直接体现为硬件处理器执行完成,或者用处理器中的硬件及软件模块组合执行完成。
存储器可能包含高速RAM存储器,也可能还包括非易失性存储NVM,例如至少一个磁盘存储器。
总线可以是工业标准体系结构(Industry Standard Architecture,ISA)总线、外部设备互连(Peripheral Component,PCI)总线或扩展工业标准体系结构(ExtendedIndustry Standard Architecture,EISA)总线等。总线可以分为地址总线、数据总线、控制总线等。为便于表示,本申请附图中的总线并不限定仅有一根总线或一种类型的总线。
本发明实施例还提供一种计算机可读存储介质,所述计算机可读存储介质中存储有计算机执行指令,当处理器执行所述计算机执行指令时,实现上述图2对应的实施例的评价储层孔隙结构的方法。
上述的计算机可读存储介质,上述可读存储介质可以是由任何类型的易失性或非易失性存储设备或者它们的组合实现,如静态随机存取存储器(SRAM),电可擦除可编程只读存储器(EEPROM),可擦除可编程只读存储器(EPROM),可编程只读存储器(PROM),只读存储器(ROM),磁存储器,快闪存储器,磁盘或光盘。可读存储介质可以是通用或专用计算机能够存取的任何可用介质。
一种示例性的可读存储介质耦合至处理器,从而使处理器能够从该可读存储介质读取信息,且可向该可读存储介质写入信息。当然,可读存储介质也可以是处理器的组成部分。处理器和可读存储介质可以位于专用集成电路(Application Specific IntegratedCircuits,简称:ASIC)中。当然,处理器和可读存储介质也可以作为分立组件存在于设备中。
本领域普通技术人员可以理解:实现上述各方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成。前述的程序可以存储于一计算机可读取存储介质中。该程序在执行时,执行包括上述各方法实施例的步骤;而前述的存储介质包括:ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (9)
1.一种评价储层孔隙结构的方法,其特征在于,包括:
获取目标储层的核磁共振回波数据;
对所述核磁共振回波数据进行时域变换,并计算变换后的核磁共振回波数据的振幅的平方及频率;
对所述振幅的平方和频率进行线性拟合,得到拟合直线;
获取所述拟合直线对应的斜率,并根据所述斜率计算所述核磁共振回波数据的分形维数;
根据所述分形维数评价所述目标储层的孔隙结构;
按照预设的傅立叶变换公式,对所述核磁共振回波数据进行傅立叶变换;
所述计算变换后的核磁共振回波数据的振幅的平方及频率,包括:
获取核磁共振回波数据的核磁共振回波个数以及所述傅立叶变换公式中的实部和虚部;
根据所述核磁共振回波个数和所述傅立叶变换公式中的实部与虚部,计算所述变换后的核磁共振回波数据的振幅的平方;
获取所述核磁共振回波数据的回波间隔;
根据所述回波间隔和所述核磁共振回波个数计算所述核磁共振回波数据中的采样点的频率。
3.根据权利要求1所述的方法,其特征在于,所述对所述核磁共振回波数据进行时域变换之前,还包括:
判断目标储层核磁共振回波数据对应的概率分布是否符合高斯分布,其中目标储层核磁共振回波数据为所述核磁共振回波数据中的任一个;
若所述目标储层核磁共振回波数据的概率分布不符合高斯分布,则对所述目标储层核磁共振回波数据进行正态变换。
6.根据权利要求1所述的方法,其特征在于,所述对所述振幅的平方和频率进行线性拟合,包括:
在预设的双对数坐标系中,对所述振幅的平方和频率进行线性拟合。
7.一种评价储层孔隙结构的设备,其特征在于,包括:
回波数据获取模块,用于获取目标储层的核磁共振回波数据;
回波数据变换模块,用于对所述核磁共振回波数据进行时域变换,
并计算变换后的核磁共振回波数据的振幅的平方及频率;
线性拟合模块,用于对所述振幅的平方和频率进行线性拟合,得到拟合直线;
分形维数计算模块,用于获取所述拟合直线对应的斜率,并根据所述斜率计算所述核磁共振回波数据的分形维数;
孔隙结构评价模块,用于根据所述分形维数评价所述目标储层的孔隙结构;
所述回波数据变换模块包括:傅立叶变换单元;
所述傅立叶变换单元,用于按照预设的傅立叶变换公式,对所述核磁共振回波数据进行傅立叶变换;
所述回波数据变换模块包括:
变换数据获取单元,用于获取核磁共振回波数据的核磁共振回波个数以及所述傅立叶变换公式中的实部和虚部;
振幅平方计算单元,用于根据所述核磁共振回波个数和所述傅立叶变换公式中的实部与虚部,计算所述变换后的核磁共振回波数据的振幅的平方;
回波间隔获取单元,用于获取所述核磁共振回波数据的回波间隔;
频率计算单元,用于根据所述回波间隔和所述核磁共振回波个数计算所述核磁共振回波数据中的采样点的频率。
8.一种评价储层孔隙结构的设备,其特征在于,包括:至少一个处理器和存储器;
所述存储器存储计算机执行指令;
所述至少一个处理器执行所述存储器存储的计算机执行指令,使得所述至少一个处理器执行如权利要求1至6任一项所述的评价储层孔隙结构的方法。
9.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质中存储有计算机执行指令,当处理器执行所述计算机执行指令时,实现如权利要求1至6任一项所述的评价储层孔隙结构的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910047194.3A CN109736793B (zh) | 2019-01-18 | 2019-01-18 | 一种评价储层孔隙结构的方法及设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910047194.3A CN109736793B (zh) | 2019-01-18 | 2019-01-18 | 一种评价储层孔隙结构的方法及设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109736793A CN109736793A (zh) | 2019-05-10 |
CN109736793B true CN109736793B (zh) | 2021-01-29 |
Family
ID=66365195
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910047194.3A Active CN109736793B (zh) | 2019-01-18 | 2019-01-18 | 一种评价储层孔隙结构的方法及设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109736793B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112734143A (zh) * | 2019-10-14 | 2021-04-30 | 中国石油大学(北京) | 一种孔隙结构分形表征的方法、装置及系统 |
CN110987751B (zh) * | 2019-11-15 | 2021-05-18 | 东北石油大学 | 一种致密储层三维空间孔喉定量分级评价方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106525681A (zh) * | 2016-10-14 | 2017-03-22 | 中国石油大学(北京) | 一种确定页岩储层孔隙直径大小的方法 |
CN108596383A (zh) * | 2018-04-18 | 2018-09-28 | 中国石油大学(北京) | 储层分类的方法及装置 |
CN109100812A (zh) * | 2018-07-25 | 2018-12-28 | 中国石油大学(北京) | 基于核磁共振的岩石孔隙分形维数评价方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10228484B2 (en) * | 2015-10-30 | 2019-03-12 | Schlumberger Technology Corporation | Robust multi-dimensional inversion from wellbore NMR measurements |
-
2019
- 2019-01-18 CN CN201910047194.3A patent/CN109736793B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106525681A (zh) * | 2016-10-14 | 2017-03-22 | 中国石油大学(北京) | 一种确定页岩储层孔隙直径大小的方法 |
CN108596383A (zh) * | 2018-04-18 | 2018-09-28 | 中国石油大学(北京) | 储层分类的方法及装置 |
CN109100812A (zh) * | 2018-07-25 | 2018-12-28 | 中国石油大学(北京) | 基于核磁共振的岩石孔隙分形维数评价方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN109736793A (zh) | 2019-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109736793B (zh) | 一种评价储层孔隙结构的方法及设备 | |
CN109340586B (zh) | 一种供水管道泄露的检测方法及系统 | |
CN110288003B (zh) | 数据变化识别方法及设备 | |
CN109060592B (zh) | 一种页岩损失气量获取方法、装置及终端设备 | |
CN117590172A (zh) | 应用于变压器的局部放电声电联合定位方法、装置和设备 | |
CN114168586A (zh) | 一种异常点检测的方法和装置 | |
CN113465734B (zh) | 一种结构振动的实时估计方法 | |
CN109945075B (zh) | 一种供水管道泄漏程度检测方法及装置 | |
CN111596379A (zh) | 一种地震观测系统异常识别方法和装置 | |
CN116775409A (zh) | 功耗评估方法、装置、电子设备及计算机可读存储介质 | |
CN104832161B (zh) | 基于双尺度相关对比的自动深度校正方法及系统 | |
CN110263496B (zh) | 岩芯特征尺寸确定方法及装置 | |
CN114062190B (zh) | 表征页岩油微观赋存特征的方法、装置、终端及存储介质 | |
CN113609971B (zh) | 微地震观测设备的检查方法、装置、设备及存储介质 | |
CN113960555B (zh) | 目标太赫兹时域回波处理方法、装置、设备及存储介质 | |
CN117554593B (zh) | 一种农田土壤墒情实时监测方法 | |
CN109639318A (zh) | 航空通信干扰脉冲的消除方法、装置、设备及存储介质 | |
CN117270063B (zh) | 一种信号频率获取方法以及装置 | |
CN116318445B (zh) | 频谱感知方法、装置、电子设备和存储介质 | |
CN113870956B (zh) | 基于页岩油赋存状态表征的页岩游离油饱和度预测方法 | |
CN113433513B (zh) | 一种雷达脉冲信号信道检测方法 | |
CN110778315B (zh) | 油井井筒脱气诊断方法 | |
CN117741765A (zh) | 一种空变子波提取方法、装置、存储介质及电子设备 | |
CN118444385A (zh) | 基于等间隔分立谱特征的工频干扰噪声压制方法及系统 | |
CN113723225A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |