CN111832975A - 一种用于原地浸出矿床的非均质性评价方法、装置及介质 - Google Patents
一种用于原地浸出矿床的非均质性评价方法、装置及介质 Download PDFInfo
- Publication number
- CN111832975A CN111832975A CN202010744233.8A CN202010744233A CN111832975A CN 111832975 A CN111832975 A CN 111832975A CN 202010744233 A CN202010744233 A CN 202010744233A CN 111832975 A CN111832975 A CN 111832975A
- Authority
- CN
- China
- Prior art keywords
- spatial information
- point cloud
- heterogeneity
- information entropy
- maximum outer
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Mining
Abstract
本发明公开了一种用于原地浸出矿床的非均质性评价方法、装置及介质,其中方法以矿层块段模型为硬数据,首先计算每个封闭空间各属性类别所在空间区域的最大外轮廓;然后用多点样板来扫描最大外轮廓得到局部空间信息熵,并按最大外轮廓、封闭空间、待评价区域这样从局部到整体,依次对空间信息熵进行融合,从而获得整个待评价区域的联合空间信息熵;最后,根据建立的信息熵与非均质性的映射关系,来评价采区的非均质性。本发明可以有效反映采区的非均质性,从而可实现对原地浸出矿床非均质性的有效自动化评价,且评价结果可用于指导地浸区生产钻孔布置及井场结构的优化,具有重大的实践意义。
Description
技术领域
本发明涉及原地浸出矿床储层评价技术领域,特别地,涉及一种用于原地浸出矿床的非均质性评价方法、装置及介质。
背景技术
原地浸出过程中矿化非均质性和渗透特性是决定能否实施地浸开采的关键技术,两者缺一不可。与此同时,影响地浸开采技术是否经济合理的因素有矿体品位、矿物成分、产状、埋深以及储量等。长期以来关于储层评价技术的研究多集中于油气藏领域,而与砂岩铀矿储层评价技术相关的研究少之甚少。因此,在砂岩铀矿山井场设计工作中,人们通常基于经验,考虑地下水流方向,按照统一的井型、井距等进行井场平面布置,并根据单个注入井或抽出井所揭露的岩性和品位信息,完成垂向的井构布置。这无疑将导致采区建设和生产成本增加、设备材料浪费、生产效率和效益降低等诸多问题。
因此,为了保证砂岩型铀矿床原地浸出开采设计方案的合理性、最优性,我们亟需在三维地质块段模型的基础上,建立一套能够自动化定量评价溶浸开采矿床储层非均质性的方法和指标,来指导采区井场结构的优化,达到降低建设投资、生产成本和提高目标矿物的回收率的目的。
发明内容
本发明提供一种用于原地浸出矿床的非均质性评价方法,通过计算采区的信息熵,然后根据信息熵与非均质性的映射关系,来有效评价原地浸出矿床储层各属性空间展布的非均质性,评价结果可用于指导地浸区生产钻孔布置及井场结构的优化。
为实现上述技术目的,本发明采用如下技术方案:
一种用于原地浸出矿床的非均质性评价方法,包括以下步骤:
步骤S1,根据钻孔的开口坐标,将非均质性待评价区域划分为若干封闭空间;
步骤S2,获取非均质性待评价区域对应的地质块段模型,从中提取所有点云数据,并判断每个点云所属的封闭空间;
步骤S3,针对每个封闭空间:均按照其中点云的属性类别,在封闭空间内生成与属性类别对应的最大外轮廓;
步骤S4,定义多点样板的尺寸;针对每个封闭空间的每个最大外轮廓:均使用所有的多点样板扫描最大外轮廓内的点云,获得所有多点样板的概率分布,进而计算最大外轮廓的空间信息熵;
步骤S5,针对每个封闭空间:根据其所有的最大外轮廓计算封闭空间的空间信息熵;
步骤S6,聚合所有封闭空间的空间信息熵,得到能够评价整个待评价区域的联合空间信息熵;
步骤S7,根据步骤S5和步骤S6得到的空间信息熵,以及空间信息熵与非均质性的映射关系,评价非均质性待评价区域内每个封闭空间以及整个待评价区域的非均质性。
进一步的,步骤S2中,判断每个点云所属的封闭空间的方法为:
设点云坐标为A(x,y),封闭空间为多边形P{p1,p2,…,pK},将点云与多边形的每条边进行叉乘,即:{Ap1×Ap2,Ap2×Ap3,…,ApK-1×ApK,ApK×Ap1},如果计算结果为0或者符号相同,则表示点云A(x,y)在多边形内或在多边形上;反之,表示点云A(x,y)在多边形外;
其中,P表示多边形,pi表示多边形的第i个顶点,i=1,2,…,K,多边形的所有顶点按照顺时针或逆时针方向排列依次为p1,p2,…,pK,Api为点云A到顶点pi的矢量。
进一步的,步骤S3中,按照其中点云的属性类别,生成与属性类别对应的最大外轮廓,具体过程为:
设当前封闭空间内的所有点云集合为Ω,不同属性类别的点云集合为Ωm,m=1,2,…,M;采用Alpha Shapes算法计算点云集合Ωm的最大外轮廓为Cm;
其中,计算点云集合Ωm中的每个点云A(x,y)是否在最大外轮廓Cm内的方法为:
将点云A(x,y)与最大外轮廓Cm{s1,s2,…,sN}的顶点进行叉乘,即:{As1×As2,As2×As3,…,AsN-1×AsN,AsN×As1},如果计算结果为0或者符号相同,则表示点云A(x,y)在最大外轮廓Cm的内部或在最大外轮廓Cm的面上;反之,表示点云A(x,y)在最大外轮廓Cm的外部;
式中,si表示最大外轮廓Cm的第i个顶点,i=1,2,…,N,最大外轮廓Cm的所有顶点按照顺时针或逆时针方向依次排列为s1,s2,…,sN,Asi为点云A到顶点si的矢量。
进一步的,步骤S4中,扫描最大外轮廓内的点云,是指点云集合Ω分布在最大外轮廓内的点云;
定义多点样板的网格尺寸为l*w*h,由于每个点云的属性类别包括M种,故多点样板共有Num=Ml*w*h种不同配置,定义不同配置的多点样板表示为Temk,k=1,2,…,Num;
使用第k个多点样板Temk扫描最大外轮廓Cm内的点云时,若多点样板在最大外轮廓Cm内,计为1次有效模拟;当多点样板Temk扫描完最大外轮廓Cm的所有点云后,统计多点样板Temk对应的有效模拟数量Numk;
当统计完所有多点样板的有效模拟数量后,计算所有多点样板的概率分布pk:
最大外轮廓Cm的空间信息熵hm的计算方法为:
式中,μ是一个正常数。
进一步的,所述步骤S5中,封闭空间的空间信息熵H的计算方法为:
式中,Pm为最大外轮廓Cm内的点云数量。
进一步的,设第j个封闭空间的空间信息熵为Hj,则聚合所有封闭空间的空间信息熵的具体聚合方法为:
式中,H即为聚合得到的非均质性待评价区域的空间信息熵,Wj表示第j个封闭空间的权重,J为待评价区域所包括的封闭空间数量,且
式中,Ωj为第j个封闭空间内的点云数量。
进一步的,空间信息熵与非均质性的映射关系为:若空间信息熵的范围为0~0.2,对应的封闭空间或者整个待评价区域为均质区;若空间信息熵的范围为0.2~0.4,对应的封闭空间或者整个待评价区域为弱非均质区;若空间信息熵的范围为0.4~0.6,对应的封闭空间或者整个待评价区域为非均质区;若空间信息熵的范围为0.6~1,对应的封闭空间或者整个待评价区域为强非均质区。
进一步的,所述点云为地质块段模型中每个基础块的中心点,点云的属性类别即为基础块的属性类别,所述属性类别包括高品位和低品位2种。
本发明还提供一种用于原地浸出矿床的非均质性评价装置,包括:
封闭空间划分模块,用于:根据钻孔的开口坐标,将非均质性待评价区域划分为若干封闭空间;
点云数据划分模块,用于:获取非均质性待评价区域对应的地质块段模型,从中提取所有点云数据,并判断每个点云所属的封闭空间;
最大外轮廓生成模块,用于:针对每个封闭空间均按照其中点云的属性类别,生成与属性类别对应的最大外轮廓;
第一局部空间信息熵计算模块,用于:定义多点样板的尺寸;针对每个封闭空间的每个最大外轮廓:均使用所有的多点样板扫描最大外轮廓内的点云,获得所有多点样板的概率分布,进而计算最大外轮廓的空间信息熵;
第二局部空间信息熵计算模块,用于:针对每个封闭空间,根据其所有的最大外轮廓计算封闭空间的空间信息熵;
空间信息熵聚合模块,用于:聚合所有封闭空间的空间信息熵,得到非均质性待评价区域的空间信息熵;
非均质性评价模块,用于:根据封闭空间的空间信息熵和非均质性待评价区域的空间信息熵,以及空间信息熵与非均质性的映射关系,评价非均质性待评价区域内每个封闭空间以及整个待评价区域的非均质性。
本发明还提供一种计算机存储介质,其特征在于,用于存储程序,所述程序被执行时,用于实现上述任一方法技术方案。
有益效果
本发明提供的用于原地浸出矿床的非均质性评价方法、装置及介质,以矿层块段模型为硬数据,首先计算每个封闭空间各属性类别所在空间区域的最大外轮廓;然后用多点样板来扫描最大外轮廓得到局部空间信息熵,并按最大外轮廓、封闭空间、待评价区域这样从局部到整体,依次对空间信息熵进行融合,从而获得整个待评价区域的联合空间信息熵;最后,根据建立的信息熵与非均质性的映射关系,来评价采区的非均质性。由于本发明中计算的空间信息熵是根据所有不同配置的多点样板的概率分布计算得到,可以有效反映采区的非均质性,从而可实现对原地浸出矿床非均质性的有效自动化评价,且评价结果可用于指导地浸区生产钻孔布置及井场结构的优化,具有重大的实践意义。
附图说明
构成本申请的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是本发明优选实施例的原地浸出开采矿床非均质性定量评价方法的流程示意图。
图2是本发明优选实施例的步骤S1中地质人员根据经验划分井型,井型为规则五点型或其他类型示意图。
图3是本发明优选实施例的步骤S1中根据工程布置,划分块段模型点云数据到对应封闭空间示意图。
图4是本发明优选实施例的步骤S2中矿区品位模型空间展布图。
图5是本发明优选实施例的步骤S3中Alpha Shapes算法圈定散点最大外轮廓,并填充空区示意图。
图6是本发明优选实施例的步骤S4中多点样板分布示意图。
图7是本发明优选实施例的步骤S4中多点样板扫描最大外轮廓示意图。
图8是本发明优选实施例的步骤S7中建立信息熵与非均质性的映射关系,评价各封闭空间及采区模型的非均质性示意图。
具体实施方式
下面对本发明的实施例作详细说明,本实施例以本发明的技术方案为依据开展,给出了详细的实施方式和具体的操作过程,对本发明的技术方案作进一步解释说明。
如图1所示,本发明优选实施例提供的用于原地浸出矿床非均质性评价方法,包括如下步骤:
步骤S1,根据钻孔的开口坐标,将非均质性待评价区域划分为若干封闭空间;
具体过程为:溶浸采矿中地质人员会根据经验划分井型,通常井型为规则的五点型或其他类型。提取钻孔开口坐标来组建封闭区域多边形,即封闭空间,并以此来划分地质块段模型,如图2所示。
步骤S2,获取非均质性待评价区域对应的地质块段模型,从中提取所有点云数据,并判断每个点云所属的封闭空间;
具体过程为:地质块段模型是由大量基础块组成,提取点云数据是指提取地质块段模型中每个基础块的中心点(x,y,z)、尺寸(xlen,ylen,zlen)及基础块包含的属性(property),其中的中心点即为基础块对应的点云,基础块的属性一般包括若干种不同的类别。根据点云与封闭区域的空间关系来划分点云数据,如图3所示,将地质块段模型划分为36个封闭区域;本实施例以研究区品位模型为例,该模型中包含高品位和低品位2类属性(如图4所示),根据属性类别依次给点云数据添加相应的标签{1,2};
设点云A(x,y,z)在xy平面的投影坐标为A′(x,y),封闭空间为多边形P{p1,p2,…,pK},将点云与多边形的每条边进行叉乘,即:{A′p1×A′p2,A′[2×A′p3,…,A′pK-1×A′pK,A′pK×A′p1},如果计算结果为0或者符号相同,则表示点云A′(x,y)在多边形内或在多边形上,点云A属于多边形封闭空间P;反之,表示点云A′(x,y)在多边形外,点云A不属于多边形封闭空间P;
其中,P表示多边形,pi表示多边形的第i个顶点,i=1,2,K,多边形的所有顶点按照顺时针或逆时针方向排列依次为p1,p2,…,pK,Api为点云A到顶点pi的矢量。
步骤S3,针对每个封闭空间:均按照其中点云的属性类别,生成与属性类别对应的最大外轮廓;
具体过程为:如图5所示,当前封闭空间内的所有点云集合为Ω,不同属性类别的点云集合为Ωm,m=1,2,…,M;采用Alpha Shapes算法计算点云集合Ωm的最大外轮廓为Cm,在计算中,‘Alpha’取值为Inf(正无穷)。此时,用点云集合Ω与最大外轮廓求交集来填充空区,此时最大外轮廓内的点云集合为P1。
其中,计算点云集合Ωm中的每个点云A(x,y,z)是否在最大外轮廓Cm内的方法为:
将点云A(x,y,z)与最大外轮廓Cm{s1,s2,…,sN}的顶点进行叉乘,即:{As1×As2,As2×As3,…,AsN-1×AsN,AsN×As1},如果计算结果为0或者符号相同,则表示点云A(x,y,z)在最大外轮廓Cm的内部或在最大外轮廓Cm的面上;反之,表示点云A(x,y,z)在最大外轮廓Cm的外部;
式中,si表示最大外轮廓Cm的第i个顶点,i=1,2,N,最大外轮廓Cm的所有顶点按照顺时针或逆时针方向依次排列为s1,s2,…,sN,Asi为点云A到顶点si的矢量。
步骤S4,定义多点样板的尺寸;针对每个封闭空间的每个最大外轮廓:均使用所有的多点样板扫描最大外轮廓内的点云,获得所有多点样板的概率分布,进而计算最大外轮廓的空间信息熵;
具体过程为:设定义多点样板的网格尺寸为l*w*h,且多点样板在x,y,z方向的尺寸均小于最大外轮廓的尺寸;由于每个点云的属性类别包括M种,故多点样板共有Num=Ml *w*h种不同配置,定义不同配置的多点样板表示为Temk,k=1,2,…,Num。在本实施例中,属性类别包括高品位和低品位两种,即M=2,例如定义品位值G>=0.5‰的基础块为高品位,0<G<0.5‰的基础块为低品位,该品位值可以通俗理解为基础块中有用元素的含量高低;另外本实施例定义的多点样板网格尺寸为2*2*1,即在x,y,z方向上分别有2、2、1个网格,因此共有16(24)种不同配置的多点样本,如图6所示。
使用第k个多点样板Temk扫描最大外轮廓Cm内的点云时,若多点样板在最大外轮廓Cm内,计为1次有效模拟,如图7所示;当多点样板Temk扫描完最大外轮廓Cm的所有点云后,统计多点样板Temk对应的有效模拟数量Numk;
当统计完所有多点样板的有效模拟数量后,计算所有多点样板的概率分布pk:
最大外轮廓Cm的空间信息熵hm的计算方法为:
式中,μ是一个正常数。对数底取多点样板中有效的多点样板个数,保证H∈[0,1]。
步骤S5,针对每个封闭空间:根据其所有的最大外轮廓计算封闭空间的空间信息熵;
具体地,由于最大外轮廓与按属性类别对应,因此每个封闭空间内最大外轮廓的数量与属性类别的数量相同,在本实施例中为M=2。本实施例中,针对每个封闭空间,均基于加法的方法聚合M=2个最大外轮廓的空间信息熵:
式中,Pm为最大外轮廓Cm内的点云数量。
步骤S6,聚合所有封闭空间的空间信息熵,得到非均质性待评价区域的空间信息熵;
具体地,本实施例采用基于乘法的方法聚合36个封闭空间的空间信息熵。设第j个封闭空间的空间信息熵为Hj,则聚合所有封闭空间的空间信息熵的具体聚合方法为:
式中,H即为聚合得到的非均质性待评价区域的空间信息熵,Wj表示第j个封闭空间的权重,J为待评价区域所包括的封闭空间数量,且
式中,Ωj为第j个封闭空间内的点云数量。
步骤S7,根据步骤S5和步骤S6得到的空间信息熵,以及空间信息熵与非均质性的映射关系,评价非均质性待评价区域内每个封闭空间以及整个待评价区域的非均质性;
其中,空间信息熵与非均质性的映射关系,是由地质专家根据经验预先制定,在本实施例中具体为:若空间信息熵的范围为0~0.2,对应的封闭空间或者整个待评价区域为均质区;若空间信息熵的范围为0.2~0.4,对应的封闭空间或者整个待评价区域为弱非均质区;若空间信息熵的范围为0.4~0.6,对应的封闭空间或者整个待评价区域为非均质区;若空间信息熵的范围为0.6~1,对应的封闭空间或者整个待评价区域为强非均质区。
与上述方法实施例相对应,本发明还提供一种用于原地浸出矿床的非均质性评价装置,包括:
封闭空间划分模块,用于:根据钻孔的开口坐标,将非均质性待评价区域划分为若干封闭空间;
点云数据划分模块,用于:获取非均质性待评价区域对应的地质块段模型,从中提取所有点云数据,并判断每个点云所属的封闭空间;
最大外轮廓生成模块,用于:针对每个封闭空间均按照其中点云的属性类别,生成与属性类别对应的最大外轮廓;
第一局部空间信息熵计算模块,用于:定义多点样板的尺寸;针对每个封闭空间的每个最大外轮廓:均使用所有的多点样板扫描最大外轮廓内的点云,获得所有多点样板的概率分布,进而计算最大外轮廓的空间信息熵;
第二局部空间信息熵计算模块,用于:针对每个封闭空间,根据其所有的最大外轮廓计算封闭空间的空间信息熵;
空间信息熵聚合模块,用于:聚合所有封闭空间的空间信息熵,得到非均质性待评价区域的空间信息熵;
非均质性评价模块,用于:根据封闭空间的空间信息熵和非均质性待评价区域的空间信息熵,以及空间信息熵与非均质性的映射关系,评价非均质性待评价区域内每个封闭空间以及整个待评价区域的非均质性。
本发明还提供一种计算机存储介质,用于存储程序,所述程序被执行时,用于实现上述方法实施例。
本发明以矿层块段模型为硬数据,首先计算每个封闭空间各属性类别所在空间区域的最大外轮廓;然后用多点样板来扫描最大外轮廓得到局部空间信息熵,并按最大外轮廓、封闭空间、待评价区域这样从局部到整体,依次对空间信息熵进行融合,从而获得整个待评价区域的联合空间信息熵;最后,根据建立的信息熵与非均质性的映射关系,来评价采区的非均质性。由于本发明中计算的空间信息熵是根据所有不同配置的多点样板的概率分布计算得到,可以有效反映采区的非均质性,从而可实现对原地浸出矿床非均质性的有效自动化评价,且评价结果可用于指导地浸区生产钻孔布置及井场结构的优化,具有重大的实践意义。
以上实施例为本申请的优选实施例,本领域的普通技术人员还可以在此基础上进行各种变换或改进,在不脱离本申请总的构思的前提下,这些变换或改进都应当属于本申请要求保护的范围之内。
Claims (10)
1.一种用于原地浸出矿床的非均质性评价方法,其特征在于,包括以下步骤:
步骤S1,根据钻孔的开口坐标,将非均质性待评价区域划分为若干封闭空间;
步骤S2,获取非均质性待评价区域对应的地质块段模型,从中提取所有点云数据,并判断每个点云所属的封闭空间;
步骤S3,针对每个封闭空间:均按照其中点云的属性类别,在封闭空间内生成与属性类别对应的最大外轮廓;
步骤S4,定义多点样板的尺寸;针对每个封闭空间的每个最大外轮廓:均使用所有的多点样板扫描最大外轮廓内的点云,获得所有多点样板的概率分布,进而计算最大外轮廓的空间信息熵;
步骤S5,针对每个封闭空间:根据其所有的最大外轮廓计算封闭空间的空间信息熵;
步骤S6,聚合所有封闭空间的空间信息熵,得到能够评价整个待评价区域的联合空间信息熵;
步骤S7,根据步骤S5和步骤S6得到的空间信息熵,以及空间信息熵与非均质性的映射关系,评价非均质性待评价区域内每个封闭空间以及整个待评价区域的非均质性。
2.根据权利要求1所述的方法,其特征在于,步骤S2中,判断每个点云所属的封闭空间的方法为:
设点云坐标为A(x,y),封闭空间为多边形P{p1,p2,…,pK},将点云与多边形的每条边进行叉乘,即:{Ap1×Ap2,Ap2×Ap3,…,ApK-1×ApK,ApK×Ap1},如果计算结果为0或者符号相同,则表示点云A(x,y)在多边形内或在多边形上;反之,表示点云A(x,y)在多边形外;
其中,P表示多边形,pi表示多边形的第i个顶点,i=1,2,…,K,多边形的所有顶点按照顺时针或逆时针方向排列依次为p1,p2,…,pK,Api为点云A到顶点pi的矢量。
3.根据权利要求1所述的方法,其特征在于,步骤S3中,按照其中点云的属性类别,生成与属性类别对应的最大外轮廓,具体过程为:
设当前封闭空间内的所有点云集合为Ω,不同属性类别的点云集合为Ωm,m=1,2,…,M;采用Alpha Shapes算法计算点云集合Ωm的最大外轮廓为Cm;
其中,计算点云集合Ωm中的每个点云A(x,y)是否在最大外轮廓Cm内的方法为:
将点云A(x,y)与最大外轮廓Cm{s1,s2,…,sN}的顶点进行叉乘,即:{As1×As2,As2×As3,…,AsN-1×AsN,AsN×As1},如果计算结果为0或者符号相同,则表示点云A(x,y)在最大外轮廓Cm的内部或在最大外轮廓Cm的面上;反之,表示点云A(x,y)在最大外轮廓Cm的外部;
式中,si表示最大外轮廓Cm的第i个顶点,i=1,2,…,N,最大外轮廓Cm的所有顶点按照顺时针或逆时针方向依次排列为s1,s2,…,sN,Asi为点云A到顶点si的矢量。
4.根据权利要求1所述的方法,其特征在于,步骤S4中,扫描最大外轮廓内的点云,是指点云集合Ω分布在最大外轮廓内的点云;
定义多点样板的网格尺寸为l*w*h,由于每个点云的属性类别包括M种,故多点样板共有Num=Ml*w*h种不同配置,定义不同配置的多点样板表示为Temk,k=1,2,…,Num;
使用第k个多点样板Temk扫描最大外轮廓Cm内的点云时,若多点样板在最大外轮廓Cm内,计为1次有效模拟;当多点样板Temk扫描完最大外轮廓Cm的所有点云后,统计多点样板Temk对应的有效模拟数量Numk;
当统计完所有多点样板的有效模拟数量后,计算所有多点样板的概率分布pk:
最大外轮廓Cm的空间信息熵hm的计算方法为:
式中,μ是一个正常数。
7.根据权利要求1所述的方法,其特征在于,空间信息熵与非均质性的映射关系为:若空间信息熵的范围为0~0.2,对应的封闭空间或者整个待评价区域为均质区;若空间信息熵的范围为0.2~0.4,对应的封闭空间或者整个待评价区域为弱非均质区;若空间信息熵的范围为0.4~0.6,对应的封闭空间或者整个待评价区域为非均质区;若空间信息熵的范围为0.6~1,对应的封闭空间或者整个待评价区域为强非均质区。
8.根据权利要求1所述的方法,其特征在于,所述点云为地质块段模型中每个基础块的中心点,点云的属性类别即为基础块的属性类别。
9.一种用于原地浸出矿床的非均质性评价装置,其特征在于,包括:
封闭空间划分模块,用于:根据钻孔的开口坐标,将非均质性待评价区域划分为若干封闭空间;
点云数据划分模块,用于:获取非均质性待评价区域对应的地质块段模型,从中提取所有点云数据,并判断每个点云所属的封闭空间;
最大外轮廓生成模块,用于:针对每个封闭空间均按照其中点云的属性类别,生成与属性类别对应的最大外轮廓;
第一局部空间信息熵计算模块,用于:定义多点样板的尺寸;针对每个封闭空间的每个最大外轮廓:均使用所有的多点样板扫描最大外轮廓内的点云,获得所有多点样板的概率分布,进而计算最大外轮廓的空间信息熵;
第二局部空间信息熵计算模块,用于:针对每个封闭空间,根据其所有的最大外轮廓计算封闭空间的空间信息熵;
空间信息熵聚合模块,用于:聚合所有封闭空间的空间信息熵,得到能够评价整个待评价区域的联合空间信息熵;
非均质性评价模块,用于:根据封闭空间的空间信息熵和非均质性待评价区域的空间信息熵,以及空间信息熵与非均质性的映射关系,评价非均质性待评价区域内每个封闭空间以及整个待评价区域的非均质性。
10.一种计算机存储介质,其特征在于,用于存储程序,所述程序被执行时,用于实现如权利要求1-8任一所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010744233.8A CN111832975A (zh) | 2020-07-29 | 2020-07-29 | 一种用于原地浸出矿床的非均质性评价方法、装置及介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010744233.8A CN111832975A (zh) | 2020-07-29 | 2020-07-29 | 一种用于原地浸出矿床的非均质性评价方法、装置及介质 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111832975A true CN111832975A (zh) | 2020-10-27 |
Family
ID=72920059
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010744233.8A Pending CN111832975A (zh) | 2020-07-29 | 2020-07-29 | 一种用于原地浸出矿床的非均质性评价方法、装置及介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111832975A (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105201493A (zh) * | 2015-09-17 | 2015-12-30 | 中国石油大学(华东) | 一种识别水驱油藏储层非均质性的双示踪剂方法 |
CN105651676A (zh) * | 2016-02-17 | 2016-06-08 | 中国石油大学(华东) | 一种水平井规则开发井网下的储层非均质性表征方法 |
US20170011149A1 (en) * | 2015-07-08 | 2017-01-12 | Conocophillips Company | Geobody continuity in geological models based on multiple point statistics |
CN110687603A (zh) * | 2019-11-07 | 2020-01-14 | 中海石油(中国)有限公司 | 一种海上油田储层内部渗流屏障的地质建模方法 |
CN110702580A (zh) * | 2019-10-21 | 2020-01-17 | 西安石油大学 | 一种基于信息熵的致密砂岩储层孔喉非均质性表征方法 |
-
2020
- 2020-07-29 CN CN202010744233.8A patent/CN111832975A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170011149A1 (en) * | 2015-07-08 | 2017-01-12 | Conocophillips Company | Geobody continuity in geological models based on multiple point statistics |
CN105201493A (zh) * | 2015-09-17 | 2015-12-30 | 中国石油大学(华东) | 一种识别水驱油藏储层非均质性的双示踪剂方法 |
CN105651676A (zh) * | 2016-02-17 | 2016-06-08 | 中国石油大学(华东) | 一种水平井规则开发井网下的储层非均质性表征方法 |
CN110702580A (zh) * | 2019-10-21 | 2020-01-17 | 西安石油大学 | 一种基于信息熵的致密砂岩储层孔喉非均质性表征方法 |
CN110687603A (zh) * | 2019-11-07 | 2020-01-14 | 中海石油(中国)有限公司 | 一种海上油田储层内部渗流屏障的地质建模方法 |
Non-Patent Citations (4)
Title |
---|
曹茜等: "砂砾岩储层分布非均质性和质量非均质性研究――以克拉玛依油田五2东区克上组为例", 《岩性油气藏》 * |
李月: "镇北地区长8油层组储层非均质性的研究", 《中国优秀硕士学位论文全文数据库》 * |
窦文超: "鄂尔多斯盆地西南部长6-长7段砂岩致密成因及非均质性研究", 《中国优秀硕士学位论文全文数据库》 * |
贺康: "Entropy Evaluation Method for Sandstone Uranium Reservoir Characteristics Based on Convex Hull Search", 《IEEE ACCESS(VOLUME:8)》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2018222331A1 (en) | Constructing structural models of the subsurface | |
CN106529755A (zh) | 一种矿山地质资源储量管理方法 | |
CN104809266B (zh) | 一种基于样条曲线的工作面煤层赋存情况精确预测方法 | |
CN102609982B (zh) | 空间地质数据非结构化模式的拓扑发现方法 | |
CN114663627B (zh) | 一种基于三维点云数据库的矿山数字模型建立方法 | |
CN103218850B (zh) | 一种真三维采矿爆破单元体模型建立方法 | |
CN112562078A (zh) | 一种三维地质分析预测模型构建方法 | |
Kessler et al. | Unlocking the potential of digital 3D geological subsurface models for geotechnical engineers | |
Li et al. | Uncertainty visualisation of a 3D geological geometry model and its application in GIS-based mineral resource assessment: a case study in Huayuan District, Northwestern Hunan Province, China | |
CN111832975A (zh) | 一种用于原地浸出矿床的非均质性评价方法、装置及介质 | |
CN111739165A (zh) | 一种矿体三维品位模型局部更新方法 | |
CN104462649B (zh) | 一种矿体块段模型储量的自动更新方法 | |
Zhang et al. | Method and application of urban 3D rapid modeling of geology based on CAD borehole logs | |
CN116699715A (zh) | 一种高光谱遥感信息参与的铀成矿有利区预测方法 | |
CN115035258A (zh) | 一种高效的基于cad钻孔柱状图的城市三维地质建模方法 | |
Houlding | The application of new 3-D computer modelling techniques to mining | |
Rezaei et al. | An integrated geo-statistical methodology for an optimum resource estimation of angouran underground mine | |
He et al. | Entropy evaluation method for sandstone uranium reservoir characteristics based on convex hull search | |
Sun et al. | Adaptive Interpolation Method for Generalized Triangular Prism (GTP) Geological Model Based on the Geometric Smoothness Rule | |
CN110675495B (zh) | 基于地质体建模的采煤工作面煤层三维展示方法及其应用 | |
Liu et al. | Fine geological modeling of complex fault block reservoir based on deep learning | |
CN116912534B (zh) | 自适应搜索匹配的热液矿床成矿系统空间结构识别方法 | |
CN113436327B (zh) | 一种三维地质模型建模方法、装置及可读存储介质 | |
CN106887041A (zh) | 模拟地质现象空间发育趋势的地质三维建模方法和系统 | |
Wang et al. | A Calculation Method For Mineral Reserves Based On Multi-point Geostatistics |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20201027 |
|
WD01 | Invention patent application deemed withdrawn after publication |