CN112836393A - 一种基于多尺度熵分析储层非均质性的方法 - Google Patents
一种基于多尺度熵分析储层非均质性的方法 Download PDFInfo
- Publication number
- CN112836393A CN112836393A CN202110252785.1A CN202110252785A CN112836393A CN 112836393 A CN112836393 A CN 112836393A CN 202110252785 A CN202110252785 A CN 202110252785A CN 112836393 A CN112836393 A CN 112836393A
- Authority
- CN
- China
- Prior art keywords
- logging
- reservoir
- scale
- time
- entropy
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/12—Timing analysis or timing optimisation
-
- 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
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Abstract
本发明公开了一种基于多尺度熵分析储层非均质性的方法,包括:采集待测井中储层在不同井段的测井信号,并得到待测井不同井段的测井曲线;将测井曲线进行归一化处理和粗粒化处理,并计算得到测井曲线所表征的储层非均质性在不同时间尺度下对应的样本熵值;依据时间尺度与其对应的样本熵值之间的关系绘制函数图,并依据函数图构建不同井段的储层的多尺度熵曲线;比对分析不同井段的多尺度熵曲线,来评价不同井段的储层非均质性强弱。以归一化和粗粒化的处理方式,将测井信号根据不同尺度得到多个不同的信息,多方位、多角度地提高待处理数据的信息量,同时构建多尺度熵曲线,实现全方位的比对,提高分析信息的准确可靠性。
Description
技术领域
本发明涉及储层非均质性的分析方法领域,具体涉及一种基于多尺度熵分析储层非均质性的方法。
背景技术
在过去的几十年中,各国致力于开发传统的石油和天然气,但是,随着全球能源格局的巨大变化,当前许多国家正在将其勘探和开发目标转向非常规资源。然而,在实际开采过程中,会受到很多因素的影响,特别是非均质性问题在低孔隙度和低渗透率的非常规油藏开发中的影响则更为突出,因此,大量学者对储层非均质性研究给予了极大的关注。
然而,目前主要还是针对孔隙度以及渗透率的静态资料进行统计计算来分析非均质性。虽然有一些学者从测井信号的复杂度研究,表明测井信号度越高,地层的非均质性越强,但是,这些方法受到算法和成本的限制,导致结果不准确和不恰当,并且,由于他们大多数单纯从信号的复杂度出发,没有考虑到地层系统的固有的复杂波动以及多期性。
例如:2016年,D.Subhakar应用去趋势波动分析法以及多重分形去趋势波动方法分析测井曲线的自相似性并结合赫斯特指数分析测井曲线的复杂程度;2017年,GauravS.Gairola应用希尔伯特黄变换,分析通过分析测井信号分解除的固有模态函数分量的个数以及所含能量的大小来分析原始测井信号的复杂程度,进而评价储层的非均质性;2019年,肖子亢通过R/S分形理论表分析测井曲线的长期自相似性表征测井曲线的复杂程度,然后通过应用经验模态分解方法分析井曲线的短期自相似性表征测井曲的复杂程度,最后将两个复杂程度结合提出新的复杂度评价指标。
上述几种方式中,去趋势波动分析与R/S分形理论的核心都是通过分析测井信号的自相关性,希尔伯特黄变换以及经验模态分解方法通过分析分解后所得信号的个数以及所含的能量高低作为指标来评价原始信号的复杂程度,这些方法都忽视了信号内在的固有联系以及地层单元的复杂情况,均是在单一尺度(分辨率)下进行的复杂度分析,因此,会造成分析结果偏差较大。
尽管人们正在花费很多精力来克服这些缺点,但尚未开发出一种高效有效的方法。
熵是描述系统在物理领域中的复杂性的度量(Grassberger,1991)。后来,学者们将熵理论引入信号分析领域,以描述信号的复杂性,并提出了香农(Shannon)熵,Kolmogorov熵,KS熵,近似熵和样本熵(Benettin等人,1976年;Pincus,1911年;Richman和Moorman,2000;香农,1948年)。与其他熵分析相比,样本熵分析(Sampleentropyanalysis,SEA)具有较强的抗干扰能力和抗噪声性能。样本熵分析提供了一种在时间尺度上量化数据复杂性的工具。样本熵分析有两个有点:1)对时间序列长度的依赖性低2)复杂度的一致性,也就是如果某一段数据的样本熵值高于另一段数据的样本熵值,则无论如何改变算法参数都不会改变最终的分析结果(Richman and Moorman,2000)。由于测井信号的采样频率固定,因此可以将测井信号视为一个时间序列(Subhakar和Chandrasekhar,2016年;Xiao等人,2019年)。
尽管样本熵分析可以分析时间序列的复杂性,但它也有局限性。仅基于单尺度样本熵分析,得出的复杂度信息不足。
发明内容
为此,本发明提供一种基于多尺度熵分析储层非均质性的方法,以解决现有技术中尽管样本熵分析可以分析时间序列的复杂性,但它也有局限性,仅基于单尺度样本熵分析,得出的复杂度信息不足提出的问题。
为了实现上述目的,本发明的实施方式提供如下技术方案:
本发明提供了一种基于多尺度熵分析储层非均质性的方法,包括如下步骤:
步骤100、采集待测井中储层在不同井段的测井信号,并得到所述待测井不同井段的测井曲线;
步骤200、将所述测井曲线进行归一化处理和粗粒化处理,并计算得到所述测井曲线所表征的储层非均质性在不同时间尺度下对应的样本熵值;
步骤300、依据所述时间尺度与其对应的样本熵值之间的关系绘制函数图,并依据所述函数图构建不同井段的储层的多尺度熵曲线;
步骤400、比对分析不同井段的多尺度熵曲线,来评价不同井段的储层非均质性强弱。
作为本发明一种优选地方案,所述步骤S200中所述样本熵值的计算包括:
归一化处理,将所述测井曲线中的振幅以0~1或者-1~1之间的小数表示进行数据化,得到归一化测井曲线;
粗粒化处理,测井曲线为N个元素组成的时间序列,将所述时间尺度的数值设定为τ,对不同个数的时间尺度τ内的连续点取平均值,以创建粗粒化后不同时间尺度的时间序列,即为:
当时间尺度τ=1时,粗粒化数据的结果就是原始的时间序列;
当时间尺度τ=2时,粗粒化数据的结果是通过计算两个连续元素的平均值形成的,变为长度为N/2的时间序列;
当时间尺度τ=3时,粗粒化数据的结果是通过计算两个连续元素的平均值形成的,变为长度为N/3的时间序列;
以此类推,以获得粗粒化后的时间序列,即为粗粒化测井曲线;
由此,归纳总结得到:
每个粗粒化后时间序列中的元素为:
作为本发明一种优选地方案,步骤200中,步骤200中,样本熵值的计算具体包括:
步骤201、所述粗粒化曲线组构建成N1个数据点组成的时间序列,记为:
{x}={x(1),x(2),x(3),…,x(N1)},其中,N1=N/τ;
步骤202、自i点开始选择m个连续的x值形成一组维数为m的向量序列,记为:
Xm(i)={x(i),x(i+1),x(i+2),…,x(i+m-1)},1≤i≤N1-m+1;
步骤203、定义向量Xm(i)与Xm(j)之间的最大距离d[Xm(i),Xm(j)]为两者对应元素中最大差值的绝对值,i,j均为元素的位置,得到:
d[Xm(i),Xm(j)]=maxk=0,…,m-1(|x(i+k)-x(j+k)|),0≤k≤m-1,
其中,1≤j≤N1-m+1,1≤i≤N1-m+1,且j≠i,k为递增个数;
其中,r为约束条件,B(m)(r)为Xm(i)与Xm(j)在r下匹配m个点的概率,r为相似容限;
其中,B(m+1)(r)为Xm+1(i)与Xm+1(j)在r下匹配m+1个点的概率;
作为本发明一种优选地方案,时间尺度τ对应连续粗粒化时间序列为;依据步骤S201至步骤S206计算每个时间尺度下粗粒化的时间序列的样本熵值:、
MSE(x,τ,m,r)=SE(y(τ),m,r)。
作为本发明一种优选地方案,所述相似容限r=0.1~0.25std,其中,std为所述测井曲线对应时间序列的标准差。
作为本发明一种优选地方案,还包括步骤500,依据不同井段的储层的测井曲线的信息复杂度,评估不同测井的储存非均质性强弱,包括如下步骤:
步骤501、计算每个井段的测井信号在每个时间尺度下的样本熵值;
步骤502、比较每个井段的所有样本熵值;
步骤503、当一个井段的样本熵值超过一半大于另一个测井的样本熵值时,则说明前一个井段的测井信号的复杂度高,则储存非均质性强;反之亦然。
作为本发明一种优选地方案,步骤100还包括:
步骤101、将多个测井仪器按照标记顺序装配在测量通道的多个均匀分布的测量节点上,并将所述测量通道以固定速率下放,并在测井中标注采集测井信号的起始线,依据所述固定速率和相邻所述测量节点的单元距离设定所述测井仪器采集并记录一次测井信号的单元时间;
步骤102、依据所述测量通道上多个所述测量节点的分布,将测井内待分析储层分为多个储层井段,设定的储层井段的初始单元厚度为所述相邻所述测量节点的单元距离,在第一个所述测井仪器到达所述起始线的同时启动第一个所述测井仪器进行采集;
步骤103、以所述固定速率将所述测量通道继续下放,并按照顺序且以所述单元时间为间隔依次开启每个所述测井仪器进行采集,且已开启的所述测井仪器继续采集下一个所述储层井段的测井信号;
步骤104、对在每个所述储层井段采集所述测井信号进行整合处理,剔除存在明显测量误差的测井信号并获得每个储层井段的储层优化数据,并对所述储层优化数据进行处理以获得每个所述储层井段的测井曲线。
作为本发明一种优选地方案,还包括:将所述测井仪器的采集记录一次测井信号的时间设定为至少2倍的单元时间,以获取不同厚度的所述储层井段的测井曲线,即获取不同时间尺度的井段的测井曲线。
作为本发明一种优选地方案,所述步骤100还包括:
在测量通道上设置用于连通每个所述测井仪器并接收采集数据的数据处理模块,并按照所述测量仪器启动采集的顺序进行编号;
在所述数据处理模块中设定储层井段的初始单元厚度u,设定时间尺度τ为所述测量仪器采集记录一次测井信号的时间T的调节因子,记固定速率为v,则T=τu/v,则采集的实际井段厚度H=τu;
在所述数据处理模块中以所述实际井段厚度H为第一参数,以测量仪器的编号为第二参数,以测井信号为第三参数存储每一种时间尺度τ下在每一个储层井段中每个所述测井仪器采集的测井信号的数据集合。
作为本发明一种优选地方案,将所述测量通道按照所述测井仪器的编号的正反顺序分别在所述测井内采集一次测井信号,以使得每个所述储层井段的数据采集次数相同。
本发明的实施方式具有如下优点:
以归一化和粗粒化的处理方式,将原始的测井信号根据不同尺度得到多个不同的信息,多方位、多角度地提高待处理数据的信息量;将多组信息数据分别计算得到样本熵值并采用一一对应的方式构建多尺度熵曲线,在完成整体信息的多元化构建的基础上,实现全方位的比对,提高分析信息的准确可靠性。
附图说明
为了更清楚地说明本发明的实施方式或现有技术中的技术方案,下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍。显而易见地,下面描述中的附图仅仅是示例性的,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图引伸获得其它的实施附图。
本说明书所绘示的结构、比例、大小等,均仅用以配合说明书所揭示的内容,以供熟悉此技术的人士了解与阅读,并非用以限定本发明可实施的限定条件,故不具技术上的实质意义,任何结构的修饰、比例关系的改变或大小的调整,在不影响本发明所能产生的功效及所能达成的目的下,均应仍落在本发明所揭示的技术内容得能涵盖的范围内。
图1为本发明实施方式中基于多尺度熵分析储层非均质性的方法的流程图;
图2为本发明实施方式中τ=2时,粗粒化处理的过程;
图3为本发明实施方式中τ=3时,粗粒化处理的过程;
图4为本发明实施方式中A、B、C、D四个油井的尺度与样本熵值的对应函数曲线及产量图;
图5为本发明实施方式中C、E、F、G四个油井的尺度与样本熵值的对应函数曲线及月产量值图;
图6为本发明实施方式中同一油井不同位置的尺度与样本熵值的对应函数曲线图。
具体实施方式
以下由特定的具体实施例说明本发明的实施方式,熟悉此技术的人士可由本说明书所揭露的内容轻易地了解本发明的其他优点及功效,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明提供了一种基于多尺度熵分析储层非均质性的方法,包括如下步骤:
步骤100、采集待测井中储层在不同井段的测井信号,并得到所述待测井不同井段的测井曲线;
步骤200、将所述测井曲线进行归一化处理和粗粒化处理,并计算得到所述测井曲线所表征的储层非均质性在不同时间尺度下对应的样本熵值;
步骤300、依据所述时间尺度与其对应的样本熵值之间的关系绘制函数图,并依据所述函数图构建不同井段的储层的多尺度熵曲线;
步骤400、比对分析不同井段的多尺度熵曲线,来评价不同井段的储层非均质性强弱。
本发明以以归一化和粗粒化的处理方式,将原始的测井信号根据不同尺度得到多个不同的信息,多方位、多角度地提高待处理数据的信息量;将多组信息数据分别计算得到样本熵值并采用一一对应的方式构建多尺度熵曲线,在完成整体信息的多元化构建的基础上,实现全方位的比对,提高分析信息的准确可靠性。
其中,时间尺度是指完成某一种物理过程所花费时间的平均度量,物理过程的演变越慢,其时间尺度越长,物理过程涉及的空间范围越大,其时间尺度也越长。在储层的应用中可以理解为不同深度的储层。
其中,如图2和图3所示,所述步骤S200中所述样本熵值的计算包括:
归一化处理,将所述测井曲线中的振幅以0~1或者-1~1之间的小数表示进行数据化,得到归一化测井曲线;
粗粒化处理,测井曲线为N个元素组成的时间序列,将所述时间尺度的数值设定为τ,对不同个数的时间尺度τ内的连续点取平均值,以创建粗粒化后不同时间尺度的时间序列,即为:
当时间尺度τ=1时,粗粒化数据的结果就是原始的时间序列;
当时间尺度τ=2时,粗粒化数据的结果是通过计算两个连续元素的平均值形成的,变为长度为N/2的时间序列;
当时间尺度τ=3时,粗粒化数据的结果是通过计算两个连续元素的平均值形成的,变为长度为N/3的时间序列;
以此类推,以获得粗粒化后的时间序列,即为粗粒化测井曲线;
由此,归纳总结得到:
每个粗粒化后时间序列中的元素为:
其中,骤200中,样本熵值的计算具体包括:
步骤201、所述粗粒化曲线组构建成N1个数据点组成的时间序列,记为:
{x}={x(1),x(2),x(3),…,x(N1)},其中,N1=N/τ;
步骤202、自i点开始选择m个连续的x值形成一组维数为m的向量序列,记为:
Xm(i)={x(i),x(i+1),x(i+2),…,x(i+m-1)},1≤i≤N1-m+1;
步骤203、定义向量Xm(i)与Xm(j)之间的最大距离d[Xm(i),Xm(j)]为两者对应元素中最大差值的绝对值,i,j均为元素的位置,得到:
d[Xm(i),Xm(j)]=maxk=0,…,m-1(|x(i+k)-x(j+k)|),0≤k≤m-1,
其中,1≤j≤N1-m+1,1≤i≤N1-m+1,且j≠i,k为递增个数;
其中,r为约束条件,B(m)(r)为Xm(i)与Xm(j)在r下匹配m个点的概率,r为相似容限;
其中,B(m+1)(r)为Xm+1(i)与Xm+1(j)在r下匹配m+1个点的概率;
时间尺度τ对应连续粗粒化时间序列为;依据步骤201至步骤206计算每个时间尺度下粗粒化的时间序列的样本熵值:
MSE(x,τ,m,r)=SE(y(τ),m,r)。
相似容限r=0.1~0.25std,其中,std为所述测井曲线对应时间序列的标准差。
为了使得分析比对结果与实际差距更小,提高结果的准确度,步骤S400具体包括:
S401、将两个不同的待分析储层中在同一尺度下得到的样本熵值一一对应进行比对,设定比对组数为n;
S402、对于样本熵值较高占比超过n/2的,则判定其测井信号复杂度高,非均质性强;
S403、重复步骤S401和S402,完成多个待分析储层之间的比对。
本实施方式还包括步骤500,依据不同井段的储层的测井曲线的信息复杂度,评估不同测井的储存非均质性强弱,包括如下步骤:
步骤501、计算每个井段的测井信号在每个时间尺度下的样本熵值;
步骤502、比较每个井段的所有样本熵值;
步骤503、当一个井段的样本熵值超过一半大于另一个测井的样本熵值时,则说明前一个井段的测井信号的复杂度高,则储存非均质性强;反之亦然。
同样的,利用这种方式也可以用于判断不同测井的储层非均质性强弱。
上述样本熵值较高的理解为:设定其中一个待分析储层在尺度1下的样本熵值为a1,在尺度2下的样本熵值为a2…依此类推,直至an;另一个待分析储层在尺度1下的样本熵值为b1,在尺度2下的样本熵值为b2…依此类推,直至bn;
将a1与b1进行比对,a2与b2进行比对,依此进行一一比对,每组比对结果对数值较大的分别计数在对应的待分析储层下,当数值较大的样本熵值占比超过n/2的,则判定其测井信号复杂度高,非均质性强。相应地,能够根据非均质性强弱判定其能源的产量。
上述方式从测井曲线出发,通过多维度多角度的分析方式,减少了对岩心、录井等资料依赖,降低储层非均质性分析的成本,为储层非均质性评价以及后期的产能预测与对比提供一个新手段。
采用上述方式分别对A、B、C、D四个油井为一组构建其函数图,具体如图4所示,同时右侧为其产量,通过图的左侧可以直观地看出,C、D的样本熵值的整体趋势在高位,而A、B的样本熵值的整体趋势在低位,因此,在此基础上,则可以看出A、B产量应大于C、D,比对右侧产量实际值,完全符合本申请中的判定方法得到的结论。
同时采用上述方式分别对C、E、F、G四个油井为一组构建其函数图,具体如图5所示,同时右侧为其月产量值,可以直观看出,E油井的复杂程度最低,非均质性弱,产量最高;F油井的样本熵值较高的最多,复杂程度最高,非均质性高,产量也最低;C油井和G油井虽部分样本熵值较为接近,但G油井部分样本熵值明显低于C油井,很直观地能够看出其复杂程度整体相对低于C油井,而其月产量值也很明显地给出了这个结果。
当然,这里可以为不同的油井之间的对比,也可以为同一油井不同位置的对比,例如,如图6所示,为单井中不同深度的函数图的构建,分别对应1067.8-1077.2m深度和1077.2-1087.4m深度下的储层。
上述步骤100还包括如下步骤:
步骤101、将多个测井仪器按照标记顺序装配在测量通道的多个均匀分布的测量节点上,并将所述测量通道以固定速率下放,并在测井中标注采集测井信号的起始线,依据所述固定速率和相邻所述测量节点的单元距离设定所述测井仪器采集并记录一次测井信号的单元时间;
步骤102、依据所述测量通道上多个所述测量节点的分布,将测井内待分析储层分为多个储层井段,设定的储层井段的初始单元厚度为所述相邻所述测量节点的单元距离,在第一个所述测井仪器到达所述起始线的同时启动第一个所述测井仪器进行采集;
步骤103、以所述固定速率将所述测量通道继续下放,并按照顺序且以所述单元时间为间隔依次开启每个所述测井仪器进行采集,且已开启的所述测井仪器继续采集下一个所述储层井段的测井信号;
步骤104、对在每个所述储层井段采集所述测井信号进行整合处理,剔除存在明显测量误差的测井信号并获得每个储层井段的储层优化数据,并对所述储层优化数据进行处理以获得每个所述储层井段的测井曲线。
本步骤100的实施例主要是能够充分利用测量仪器在测井内上下移动的时间获得更全面的数据,提高对测井内不同储层井段的测井信号的准确性,进而提高储层非均质性强弱的判断准确性。
在本实施例中,随着测井通道不断的下放,测井仪器不断的间歇性的采集和记录每一个按照预设要求而被区分的多个储层井段的测井信号,从而在一个测井仪器中记录了多个不同储层井段的测井信号,按照上述方法,第一个测井仪器记录全部的储层井段的测井信号,其他测井仪器采集的测井信号数量依次递减,(例如第二个测井仪器无法采集最底部的一个储层井段,依次类推)。
通过这种方式,在一次下放测井仪器的单次过程中,能够最大化的获得最多的数据,为了进一步抬高每一个储层井段的测井信号准确性,可以将所述测量通道按照所述测井仪器的编号的正反顺序分别在所述测井内采集一次测井信号,以使得每个所述储层井段的数据采集次数相同。
进一步的,为了分析不同厚度的储层井段,将所述测井仪器的采集记录一次测井信号的时间设定为至少2倍的单元时间,以获取不同厚度的所述储层井段的测井曲线,即获取不同时间尺度的井段的测井曲线。
上述步骤100还包括:
在测量通道上设置用于连通每个所述测井仪器并接收采集数据的数据处理模块,并按照所述测量仪器启动采集的顺序进行编号;
在所述数据处理模块中设定储层井段的初始单元厚度u,设定时间尺度τ为所述测量仪器采集记录一次测井信号的时间T的调节因子,记固定速率为v,则T=τu/v,则采集的实际井段厚度H=τu;
在所述数据处理模块中以所述实际井段厚度H为第一参数,以测量仪器的编号为第二参数,以测井信号为第三参数存储每一种时间尺度τ下在每一个储层井段中每个所述测井仪器采集的测井信号的数据集合。
通过这种方式,能够以时间尺度为唯一的调节因子,通过数据处理模块控制每一个测量仪器的动作(采集起始动作,采集频率等)。
本实施例可以通过人工控制测量通道的下方速率,也可以通过机械设备更为准确的控制测量通道的下放速率。
虽然,上文中已经用一般性说明及具体实施例对本发明作了详尽的描述,但在本发明基础上,可以对之作一些修改或改进,这对本领域技术人员而言是显而易见的。因此,在不偏离本发明精神的基础上所做的这些修改或改进,均属于本发明要求保护的范围。
Claims (10)
1.一种基于多尺度熵分析储层非均质性的方法,其特征在于:包括如下步骤:
步骤100、采集待测井中储层在不同井段的测井信号,并得到所述待测井不同井段的测井曲线;
步骤200、将所述测井曲线进行归一化处理和粗粒化处理,并计算得到所述测井曲线所表征的储层非均质性在不同时间尺度下对应的样本熵值;
步骤300、依据所述时间尺度与其对应的样本熵值之间的关系绘制函数图,并依据所述函数图构建不同井段的储层的多尺度熵曲线;
步骤400、比对分析不同井段的多尺度熵曲线,来评价不同井段的储层非均质性强弱。
2.根据权利要求1所述的一种基于多尺度熵分析储层非均质性的方法,其特征在于,所述步骤S200中所述样本熵值的计算包括:
归一化处理,将所述测井曲线中的振幅以0~1或者-1~1之间的小数表示进行数据化,得到归一化测井曲线;
粗粒化处理,测井曲线为N个元素组成的时间序列,将所述时间尺度的数值设定为τ,对不同个数的时间尺度τ内的连续点取平均值,以创建粗粒化后不同时间尺度的时间序列,即为:
当时间尺度τ=1时,粗粒化数据的结果就是原始的时间序列;
当时间尺度τ=2时,粗粒化数据的结果是通过计算两个连续元素的平均值形成的,变为长度为N/2的时间序列;
当时间尺度τ=3时,粗粒化数据的结果是通过计算两个连续元素的平均值形成的,变为长度为N/3的时间序列;
以此类推,以获得粗粒化后的时间序列,即为粗粒化测井曲线;
由此,归纳总结得到:
每个粗粒化后时间序列中的元素为:
3.根据权利要求2所述的一种基于多尺度熵分析储层非均质性的方法,其特征在于,步骤200中,样本熵值的计算具体包括:
步骤201、所述粗粒化曲线组构建成N1个数据点组成的时间序列,记为:
{x}={x(1),x(2),x(3),…,x(N1)},其中,N1=N/τ;
步骤202、自i点开始选择m个连续的x值形成一组维数为m的向量序列,记为:
Xm(i)={x(i),x(i+1),x(i+2),…,x(i+m-1)},1≤i≤N1-m+1;
步骤203、定义向量Xm(i)与Xm(j)之间的最大距离d[Xm(i),Xm(j)]为两者对应元素中最大差值的绝对值,i,j均为元素的位置,得到:
d[Xm(i),Xm(j)]=maxk=0,…,m-1(|x(i+k)-x(j+k)|),0≤k≤m-1,
其中,1≤j≤N1-m+1,1≤i≤N1-m+1,且j≠i,k为递增个数;
其中,r为约束条件,B(m)(r)为Xm(i)与Xm(j)在r下匹配m个点的概率,r为相似容限;
其中,B(m+1)(r)为Xm+1(i)与Xm+1(j)在r下匹配m+1个点的概率;
4.根据权利要求3所述的一种基于多尺度熵分析储层非均质性的方法,其特征在于,
时间尺度τ对应连续粗粒化时间序列为;依据步骤201至步骤206计算每个时间尺度下粗粒化的时间序列的样本熵值:
MSE(x,τ,m,r)=SE(y(τ),m,r)。
5.根据权利要求3或4所述的一种基于多尺度熵分析储层非均质性的方法,其特征在于,所述相似容限r=0.1~0.25std,其中,std为所述测井曲线对应时间序列的标准差。
6.根据权利要求4所述的一种基于多尺度熵分析储层非均质性的方法,其特征在于,还包括步骤500,依据不同井段的储层的测井曲线的信息复杂度,评估不同测井的储存非均质性强弱,包括如下步骤:
步骤501、计算每个井段的测井信号在每个时间尺度下的样本熵值;
步骤502、比较每个井段的所有样本熵值;
步骤503、当一个井段的样本熵值超过一半大于另一个测井的样本熵值时,则说明前一个井段的测井信号的复杂度高,则储存非均质性强;反之亦然。
7.根据权利要求6所述的一种基于多尺度熵分析储层非均质性的方法,其特征在于,步骤100还包括:
步骤101、将多个测井仪器按照标记顺序装配在测量通道的多个均匀分布的测量节点上,并将所述测量通道以固定速率下放,并在测井中标注采集测井信号的起始线,依据所述固定速率和相邻所述测量节点的单元距离设定所述测井仪器采集并记录一次测井信号的单元时间;
步骤102、依据所述测量通道上多个所述测量节点的分布,将测井内待分析储层分为多个储层井段,设定的储层井段的初始单元厚度为所述相邻所述测量节点的单元距离,在第一个所述测井仪器到达所述起始线的同时启动第一个所述测井仪器进行采集;
步骤103、以所述固定速率将所述测量通道继续下放,并按照顺序且以所述单元时间为间隔依次开启每个所述测井仪器进行采集,且已开启的所述测井仪器继续采集下一个所述储层井段的测井信号;
步骤104、对在每个所述储层井段采集所述测井信号进行整合处理,剔除存在明显测量误差的测井信号并获得每个储层井段的储层优化数据,并对所述储层优化数据进行处理以获得每个所述储层井段的测井曲线。
8.根据权利要求7所述的一种基于多尺度熵分析储层非均质性的方法,其特征在于,还包括:将所述测井仪器的采集记录一次测井信号的时间设定为至少2倍的单元时间,以获取不同厚度的所述储层井段的测井曲线,即获取不同时间尺度的井段的测井曲线。
9.根据权利要求8所述的一种基于多尺度熵分析储层非均质性的方法,其特征在于,所述步骤100还包括:
在测量通道上设置用于连通每个所述测井仪器并接收采集数据的数据处理模块,并按照所述测量仪器启动采集的顺序进行编号;
在所述数据处理模块中设定储层井段的初始单元厚度u,设定时间尺度τ为所述测量仪器采集记录一次测井信号的时间T的调节因子,记固定速率为v,则T=τu/v,则采集的实际井段厚度H=τu;
在所述数据处理模块中以所述实际井段厚度H为第一参数,以测量仪器的编号为第二参数,以测井信号为第三参数存储每一种时间尺度τ下在每一个储层井段中每个所述测井仪器采集的测井信号的数据集合。
10.根据权利要求7所述的一种基于多尺度熵分析储层非均质性的方法,其特征在于,将所述测量通道按照所述测井仪器的编号的正反顺序分别在所述测井内采集一次测井信号,以使得每个所述储层井段的数据采集次数相同。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110252785.1A CN112836393B (zh) | 2021-03-09 | 2021-03-09 | 一种基于多尺度熵分析储层非均质性的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110252785.1A CN112836393B (zh) | 2021-03-09 | 2021-03-09 | 一种基于多尺度熵分析储层非均质性的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112836393A true CN112836393A (zh) | 2021-05-25 |
CN112836393B CN112836393B (zh) | 2023-07-18 |
Family
ID=75930075
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110252785.1A Active CN112836393B (zh) | 2021-03-09 | 2021-03-09 | 一种基于多尺度熵分析储层非均质性的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112836393B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113869289A (zh) * | 2021-12-02 | 2021-12-31 | 西北工业大学深圳研究院 | 基于熵的多通道舰船辐射噪声特征提取方法 |
CN116663989A (zh) * | 2023-06-01 | 2023-08-29 | 中国铁道科学研究院集团有限公司 | 线状工程地形复杂性测度方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104730596A (zh) * | 2015-01-25 | 2015-06-24 | 中国石油大学(华东) | 一种基于多尺度因素约束的离散裂缝建模方法 |
CN109782341A (zh) * | 2019-02-20 | 2019-05-21 | 西安交通大学 | 一种基于信息散度的地震数据属性提取方法 |
CN110094198A (zh) * | 2019-06-05 | 2019-08-06 | 中国地质大学(北京) | 基于spss软件的油井灰质含量测定方法 |
CN110852527A (zh) * | 2019-11-20 | 2020-02-28 | 成都理工大学 | 一种结合深度学习的储层物性参数预测方法 |
US20200251201A1 (en) * | 2019-02-06 | 2020-08-06 | Oberon Sciences Ilan Ltd. | Identifying and quantifying individualized variability-patterns and implementing them for improved efficacy of systems |
-
2021
- 2021-03-09 CN CN202110252785.1A patent/CN112836393B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104730596A (zh) * | 2015-01-25 | 2015-06-24 | 中国石油大学(华东) | 一种基于多尺度因素约束的离散裂缝建模方法 |
US20200251201A1 (en) * | 2019-02-06 | 2020-08-06 | Oberon Sciences Ilan Ltd. | Identifying and quantifying individualized variability-patterns and implementing them for improved efficacy of systems |
CN109782341A (zh) * | 2019-02-20 | 2019-05-21 | 西安交通大学 | 一种基于信息散度的地震数据属性提取方法 |
CN110094198A (zh) * | 2019-06-05 | 2019-08-06 | 中国地质大学(北京) | 基于spss软件的油井灰质含量测定方法 |
CN110852527A (zh) * | 2019-11-20 | 2020-02-28 | 成都理工大学 | 一种结合深度学习的储层物性参数预测方法 |
Non-Patent Citations (2)
Title |
---|
XINYUAN JI 等: "A novel heterogeneity analysis method for clastic reservoir through multiscale entropy analysis of logging data: A case study of the Yanchang Formation, Ordos Basin, China", 《JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING》, vol. 196, pages 1 - 16 * |
张存才 等: "如何在PRMS储量分类标准下快速评估资产探讨", 《海洋石油》, vol. 40, no. 3, pages 13 - 18 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113869289A (zh) * | 2021-12-02 | 2021-12-31 | 西北工业大学深圳研究院 | 基于熵的多通道舰船辐射噪声特征提取方法 |
CN113869289B (zh) * | 2021-12-02 | 2022-03-01 | 西北工业大学深圳研究院 | 基于熵的多通道舰船辐射噪声特征提取方法 |
CN116663989A (zh) * | 2023-06-01 | 2023-08-29 | 中国铁道科学研究院集团有限公司 | 线状工程地形复杂性测度方法及装置 |
CN116663989B (zh) * | 2023-06-01 | 2024-01-23 | 中国铁道科学研究院集团有限公司 | 线状工程地形复杂性测度方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN112836393B (zh) | 2023-07-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111046341B (zh) | 一种基于主成分分析的非常规天然气压裂效果评价及产能预测方法 | |
CN112836393A (zh) | 一种基于多尺度熵分析储层非均质性的方法 | |
CN113570000A (zh) | 一种基于多模型融合的海洋单要素观测质量控制方法 | |
CN103336305B (zh) | 一种基于灰色理论划分致密砂岩储层岩石物理相的方法 | |
CN113419284B (zh) | 一种基于聚类分析测井岩石物理相双甜点识别方法 | |
CN116644284A (zh) | 一种地层分类特征因子确定方法、系统、电子设备及介质 | |
CN113283180A (zh) | 基于K-means与SVR组合的致密油藏水平井压裂产能预测方法与应用 | |
CN112305441A (zh) | 一种集成式聚类下的动力电池健康状态评估方法 | |
CN111429012A (zh) | 页岩脆性甜点评价方法 | |
CN111626377B (zh) | 一种岩相识别方法、装置、设备及存储介质 | |
CN110322363A (zh) | 页岩气储层改造体积计算方法及系统 | |
CN116562428A (zh) | 一种基于机器学习的压裂施工参数优化方法 | |
CN112282742A (zh) | 一种页岩油优质储层的预测方法 | |
CN109655394B (zh) | 一种孔喉结构参数约束下的核磁共振t2谱渗透率计算方法 | |
CN115017646A (zh) | 压裂参数优化方法、装置、电子设备和存储介质 | |
CN114912477A (zh) | 大跨度小净距隧道爆破中不同掏槽方式的识别分类方法 | |
CN113792499A (zh) | 基于沉积特征的松散承压含水层富水性动态确定方法 | |
CN114186879A (zh) | 地质参数对资源量计算误差影响的评估方法及设备 | |
CN113627640A (zh) | 一种针对缝洞型油藏油气井的产能试井预测方法及系统 | |
CN111898664B (zh) | 一种基于Block-Bootstrap和多阶段聚类的桥梁模态参数自动识别方法 | |
CN110570102A (zh) | 一种储层评价方法 | |
CN112630839A (zh) | 一种测井曲线标准化方法及系统 | |
CN111582379A (zh) | 一种基于聚类算法的岩土层智能分层方法和系统 | |
CN113449261B (zh) | 页岩气储层开采性分析方法及装置 | |
CN112835113B (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 |