CN110764163A - 基于γ总量及高密度电法测量下的伟晶岩锂矿的圈定方法 - Google Patents
基于γ总量及高密度电法测量下的伟晶岩锂矿的圈定方法 Download PDFInfo
- Publication number
- CN110764163A CN110764163A CN201911116736.4A CN201911116736A CN110764163A CN 110764163 A CN110764163 A CN 110764163A CN 201911116736 A CN201911116736 A CN 201911116736A CN 110764163 A CN110764163 A CN 110764163A
- Authority
- CN
- China
- Prior art keywords
- apparent
- gamma
- polarizability
- software
- lithium ore
- 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
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于γ总量及高密度电法测量下的伟晶岩锂矿的圈定方法,其包采集每个断面上γ照射量率及断面不同深度处的视电阻率和视极化率;标记视电阻率、视极化率和γ照射量率中的异常测量点;将视电阻率和视极化率输入反演软件得到视电阻率断面图和视极化率断面图;采用Mapgis软件生成γ总量测量的平面剖面图;将格式转换后的视电阻率断面图和视极化率断面图及平面剖面图输入Mapgis软件生成剖面综合图;采用异常测量点在剖面综合图上分别圈出第一赋存区和第二赋存区;将重合的第一赋存区和第二赋存区标记为伟晶岩脉锂矿富存区,将余下的第一赋存区和第二赋存区标记为伟晶岩脉锂矿疑似区。
Description
技术领域
本发明涉及矿床勘测技术领域,具体涉及基于γ总量及高密度电法测量下的伟晶岩锂矿的圈定方法。
背景技术
锂矿是重要的战略矿种,包括地下卤水型钾盐伴生的盐湖型锂矿和花岗伟晶岩型的硬岩型锂矿这两种,目前卤水型钾盐伴生的盐湖型锂矿的开发利用技术还未取得较好突破,新能源锂电池产业利用的锂资源还是以花岗伟晶岩型的硬岩型锂矿作为主要的材料来源。
传统伟晶岩锂矿的找矿方法是通过区域地质调查、专项地质调查、化探等手段确定靶区,然后通过大量的槽探工作进行揭露验证,在此过程中有几个问题:在四川很多地方山高林密,地形切割厉害,地质调查工作要花费大量的人力和物力,大量的槽探工作无法利用大型机械进行施工,只能人工进行探槽开挖,效率低下,工期较长,会增加更多的成本;另外,在当前形势下,在环保更严厉的要求下,提倡绿色勘查理念,传统的地质找矿通过大量的槽探揭露矿体会越来越受到约束。
发明内容
针对现有技术中的上述不足,本发明提供的基于γ总量及高密度电法测量下的伟晶岩锂矿的圈定方法在伟晶岩型锂矿探测过程中不需要开挖大量的探槽进行探测。
为了达到上发明目的,本发明采用的技术方案为:
提供一种基于γ总量及高密度电法测量下的伟晶岩锂矿的圈定方法,其包括:在目标地质体的可能走向选取若干断面;
在断面上确定若干测量点,采用高密度电法测量系统勘测断面不同深度处的视电阻率和视极化率,采用γ测量仪采集每个测量点处的γ照射量率;
将所有视电阻率中大于第一阈值、所有视极化率中大于第二阈值及所有γ照射量率中大于第三阈值的标记为异常测量点;
将所有视电阻率和视极化率单独输入高密度电法反演软件进行反演,得到视电阻率断面图和视极化率断面图;将所有γ照射量率输入Mapgis软件生成γ总量测量的平面剖面图;
将视电阻率断面图、视极化率断面图转换为能输入Mapgis软件的格式,之后将格式转换后的视电阻率断面图和视极化率断面图及γ总量测量的平面剖面图输入Mapgis软件生成剖面综合图;
根据剖面综合图及标记的异常测量点,在剖面综合图上分别圈出基于γ照射量率的第一赋存区和基于视电阻率的第二赋存区;
将第一赋存区和第二赋存区重合的部分标记为伟晶岩脉锂矿富存区,将余下的第一赋存区和第二赋存区标记为伟晶岩脉锂矿疑似区。
本发明的有益效果为:本方案通过γ测量仪和高密度电法测量系统进行目标地质体的断面(测量剖面)处的γ照射量率、视电阻率和视极化率进行采集,之后根据采集的数据与高密度电法反演软件和Mapgis软件的相互结合,能够实现基于γ总量及高密度电法测量下的目标地质体中的伟晶岩型锂矿的圈定;
通过该种方式与传统的探槽方式相比,通过仪器在目标地质体的表面进行数据采集,不需要破坏岩体,能够大幅度提高寻找伟晶岩锂矿工作效率,不存在探槽开挖可以节约探测成本;开展高密度测量及地面γ总量测量剖面工作可以较好分析判断伟晶岩锂矿深度、产状、延伸等情况;圈定伟晶岩脉锂矿的分布情况后可以进行地表异常露头查证和少量的钻探工作查证,可以很好满足绿色勘查要求。
附图说明
图1为基于γ总量及高密度电法测量下的伟晶岩锂矿的圈定方法的流程图。
图2为实施例中γ总量测量解译剖面图。
图3为实施例中γ总量剖面图。
图4为实施例中反演得到的视电阻率断面图。
图5为实施例中剖面综合图。
图6为实施例中圈定后的平面图。
图7为实施例中伟晶岩矿锂矿脉的钻孔验证图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
参考图1,图1示出了基于γ总量及高密度电法测量下的伟晶岩锂矿的圈定方法的流程图;如图1所示,该方法100包括步骤101至步骤107。
在步骤101中,在目标地质体的可能走向(勘测时获得的推测走向)选取若干断面;本方案中目标地质体的所处的工作区相对较小,一般覆盖面积在2平方公里左右;在进行断面选择时,一般在目标地质体对应的山脉走向上等间距选取,两个断面之间的距离越小,采集的γ照射量率及视电阻率和视极化率就越准确,本方案优选相邻断面的间距为200m。
在步骤102中,在断面上确定若干测量点,采用高密度电法测量系统勘测断面不同深度处的视电阻率和视极化率,采用γ测量仪采集每个测量点处的γ照射量率。
实施时,本方案优选高密度电法测量系统为DUK-2B高密度电法测量系统;γ测量仪的型号为FD-3013。
在步骤103中,将所有视电阻率中大于第一阈值、所有视极化率中大于第二阈值及所有γ照射量率中大于第三阈值的标记为异常测量点。
在本发明的一个实施例中,所述第一阈值/第二阈值/第三阈值的获取方法包括:
S1、计算所有视电阻率/视极化率/γ照射量率的平均值μ和标准偏差σ;
S2、根据视电阻率/视极化率/γ照射量率的平均值μ和标准偏差σ,计算参考阈值T=μ+3σ;
S3、判断视电阻率/视极化率/γ照射量率是否大于其对应的参考阈值;若是,则进入步骤S4,否则输出视电阻率/视极化率/γ照射量率的参考阈值T;
S4、删除大于参考阈值的视电阻率/视极化率/γ照射量率,并返回步骤S1。
在采集数据时由于操作不当或外界气体环境干扰因素的影响,使得部分数据会存在异常偏高的情况,本方法在计算第一阈值/第二阈值/第三阈值时,通过得到的对应参考值对异常偏高数据进行剔除,以使最终得到的参考阈值更能反映真实值,以保证后续伟晶岩脉锂矿赋存区圈定的准确性。
在步骤104中,将所有视电阻率和视极化率单独输入高密度电法反演软件进行反演,得到视电阻率断面图和视极化率断面图;将所有γ照射量率输入Mapgis软件生成γ总量测量的平面剖面图。
在本发明的一个实施例中,将γ照射量率输入Mapgis软件前还包括对原始数据进行预处理及数据格式转换:
将每个测量点GPS定位转换为采样坐标X、Y;
将所有测量点的γ照射量率与对应的坐标X、Y按时间顺序导入Excel表中,Excel表的第一行为标题行;
将Excel表另存为.csv格式文件,并将.csv格式文件的扩展名更改为.txt文档,得到含γ照射量率的.txt文档;
采用Mapgis软件配备的记事本打开.txt文档,并将.txt文档第一行的标题行更改为notgrid,之后将文档的扩展名更改为.dat。
实施时,本方案优选将所有γ照射量率输入Mapgis软件生成γ总量测量的平面剖面图仅有包括步骤A1至步骤A3:
在步骤A1中,采用Mapgis软件中的用户文件投影转换功能对扩展名为.dat的文档进行数据投影转换;
步骤A1具体实现步骤为:选取投影变换中的用户数据点投影变换打开“XX地区伽玛总量测量结果.dat”(只能是纯文本格式)→设置用户文件选项(选择“按指定分隔符”)→设置分隔符→指定数据起始位置→设置用户投影参数→结果投影参数→设置用户文件选项→设置生成点图元参数→投影转换→确定。如果需要将投影结果写到文本文件,那么点“写到文件”按系统提示输入投影结果文件名,输入完毕即进行转换。需要注意的是,保存的文件没有文件类型,需要手动输入.txt文件类型,才可以使用记事本打开。数据投影转换也可以采用Mapgis软件二次开发软件Section来实现。
在步骤A2中,采用Mapgis软件中的Mapgis编辑子系统根据数据投影转换后的文件中的数据生成γ总量剖面图。
步骤A2主要如下:在Mapgis软件中的主菜单采用图形处理输入编辑新建工程,打开Mapgis编辑子系统,选择“不生成可编辑项”,在“文件名”栏目下方空白处,右击鼠标,添加保存的伽玛总量测量WT、WL、WP(γ照射量率在输入Mapgis软件后,Mapgis软件能够根据输入的所有γ照射量率生成点文件WT、线文件WL和区文件WP)文件添加到系统中,复位窗口后就可以看到添加的图形内容,按照WT、WL、WP排序,点击刷新按钮。还可以按照需要添加其他已编辑好的地质、地理资料图层,通过编辑有关WT、WL、WP文件,添加文字说明、比例尺、图例、图框,就可以形成一张符合规范的地面γ总量剖面图。
在步骤A3中,在γ总量剖面图中进行地质解译,画出剖面解译剖面图,之后采用mapgis软件生成γ总量测量的平面剖面图。
在步骤105中,将视电阻率断面图、视极化率断面图转换为能输入Mapgis软件的格式,之后将格式转换后的视电阻率断面图和视极化率断面图及γ总量测量的平面剖面图输入Mapgis软件生成剖面综合图;
视电阻率断面图和视极化率断面图在输入Mapgis软件后,也会首先生成点文件WT、线文件WL和区文件WP。
在步骤106中,根据剖面综合图及标记的异常测量点,在剖面综合图上分别圈出基于γ照射量率的第一赋存区和基于视电阻率的第二赋存区;
在步骤107中,将第一赋存区和第二赋存区重合的部分标记为伟晶岩脉锂矿富存区,将余下的第一赋存区和第二赋存区标记为伟晶岩脉锂矿疑似区。
下面结合具体的实例对本方案的圈定方法的准确性进行说明:
以四川省金川县观音桥镇地区12#伟晶岩矿锂矿脉露头上采用本方案提供的方法进行伟晶岩锂矿的圈定,具体工作仪器及工作参数见表:
表1试验仪器
仪器名称 | 仪器型号 | 仪器产地 |
高密度电法测量系统 | DUK-2B | 重庆地质仪器厂 |
γ辐射仪 | FD-3013 | 上海申核仪器厂 |
将观音桥镇地区12#的长度方向等间距划分为10个断面,并在每个断面上选取10个测量点,采用本方法得到的γ总量测量解译剖面图如图2所示,采用在步骤A2中生成的γ总量剖面图如图3所示,图3中的γ1至γ4为采用步骤106圈定出来的伟晶岩脉锂矿赋存区,通过反演得到的视电阻率断面图如图4所示,通过步骤105可以得到观音桥镇地区12#处的剖面综合图,参考图5;通过步骤106的操作可以圈出第一赋存区和第二赋存区,圈定后的平面图可以参考图6,图6中重叠的区域即为伟晶岩脉锂矿富存区,圈出来的余下区域即为伟晶岩脉锂矿疑似区。
在伟晶岩脉锂矿富存区及伟晶岩脉锂矿疑似区采用进行地表异常露头查证和钻探工作查证的方式对伟晶岩脉锂矿进行验证,在伟晶岩脉锂矿富存区开5个钻孔,发现其中有4个钻孔发现了伟晶岩脉锂矿,在伟晶岩脉锂矿疑似区开设2个钻孔,发现其中有1个钻孔发现了伟晶岩脉锂矿,图7为转孔验证效果图。
通过地表异常露头查证和钻孔方式进行查证,在伟晶岩脉锂矿富存区80%的钻孔存在伟晶岩脉锂矿,在伟晶岩脉锂矿疑似区50%的钻孔存在伟晶岩脉锂矿,由此可以看出,采用本方案基于γ照射量率及视极化率和视电阻率的方式进行伟晶岩脉锂矿圈定,其准确率较高。
Claims (5)
1.基于γ总量及高密度电法测量下的伟晶岩锂矿的圈定方法,其特征在于,包括:
在目标地质体的可能走向选取若干断面;
在断面上确定若干测量点,采用高密度电法测量系统勘测断面不同深度处的视电阻率和视极化率,采用γ测量仪采集每个测量点处的γ照射量率;
将所有视电阻率中大于第一阈值、所有视极化率中大于第二阈值及所有γ照射量率中大于第三阈值的标记为异常测量点;
将所有视电阻率和视极化率单独输入高密度电法反演软件进行反演,得到视电阻率断面图和视极化率断面图;将所有γ照射量率输入Mapgis软件生成γ总量测量的平面剖面图;
将视电阻率断面图、视极化率断面图转换为能输入Mapgis软件的格式,之后将格式转换后的视电阻率断面图和视极化率断面图及γ总量测量的平面剖面图输入Mapgis软件生成剖面综合图;
根据剖面综合图及标记的异常测量点,在剖面综合图上分别圈出基于γ照射量率的第一赋存区和基于视电阻率的第二赋存区;
将第一赋存区和第二赋存区重合的部分标记为伟晶岩脉锂矿富存区,将余下的第一赋存区和第二赋存区标记为伟晶岩脉锂矿疑似区。
2.根据权利要求1所述的基于γ总量及高密度电法测量下的伟晶岩锂矿的圈定方法,其特征在于,所述第一阈值/第二阈值/第三阈值的获取方法包括:
S1、计算所有视电阻率/视极化率/γ照射量率的平均值μ和标准偏差σ;
S2、根据视电阻率/视极化率/γ照射量率的平均值μ和标准偏差σ,计算参考阈值T=μ+3σ;
S3、判断视电阻率/视极化率/γ照射量率是否大于其对应的参考阈值;若是,则进入步骤S4,否则输出视电阻率/视极化率/γ照射量率的参考阈值T;
S4、删除大于参考阈值的视电阻率/视极化率/γ照射量率,并返回步骤S1。
3.根据权利要求1所述的基于γ总量及高密度电法测量下的伟晶岩锂矿的圈定方法,其特征在于,将γ照射量率输入Mapgis软件前还包括对原始数据进行预处理及数据格式转换:
将每个测量点GPS定位转换为采样坐标X、Y;
将所有测量点的γ照射量率与对应的坐标X、Y按时间顺序导入Excel表中,Excel表的第一行为标题行;
将Excel表另存为.csv格式文件,并将.csv格式文件的扩展名更改为.txt文档,得到含γ照射量率的.txt文档;
采用Mapgis软件配备的记事本打开.txt文档,并将.txt文档第一行的标题行更改为notgrid,之后将文档的扩展名更改为.dat。
4.根据权利要求3所述的基于γ总量及高密度电法测量下的伟晶岩锂矿的圈定方法,其特征在于,所述将所有γ照射量率输入Mapgis软件生成γ总量测量的平面剖面图进一步包括:
采用Mapgis软件中的用户文件投影转换功能对扩展名为.dat的文档进行数据投影转换;
采用Mapgis软件中的Mapgis编辑子系统根据数据投影转换后的文件中的数据生成γ总量剖面图;
在γ总量剖面图中进行地质解译,画出剖面解译剖面图,之后采用mapgis软件生成γ总量测量的平面剖面图。
5.根据权利要求1所述的基于γ总量及高密度电法测量下的伟晶岩锂矿的圈定方法,其特征在于,高密度电法测量系统为DUK-2B高密度电法测量系统;高密度电法反演软件为二维高密度电法反演软件。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911116736.4A CN110764163B (zh) | 2019-11-15 | 2019-11-15 | 基于γ总量及高密度电法测量下的伟晶岩锂矿的圈定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911116736.4A CN110764163B (zh) | 2019-11-15 | 2019-11-15 | 基于γ总量及高密度电法测量下的伟晶岩锂矿的圈定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110764163A true CN110764163A (zh) | 2020-02-07 |
CN110764163B CN110764163B (zh) | 2021-06-08 |
Family
ID=69338417
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911116736.4A Active CN110764163B (zh) | 2019-11-15 | 2019-11-15 | 基于γ总量及高密度电法测量下的伟晶岩锂矿的圈定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110764163B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111398360A (zh) * | 2020-04-21 | 2020-07-10 | 山东大学 | 基于ERT和IP圈定包气带中LNAPLs的污染源区探测方法及系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2539838C1 (ru) * | 2013-07-16 | 2015-01-27 | Федеральное Государственное Унитарное Предприятие "Сибирский Научно-Исследовательский Институт Геологии, Геофизики И Минерального Сырья" | Способ прогноза скрытого оруденения, связанного с гранитоидами |
CN105700040A (zh) * | 2016-03-05 | 2016-06-22 | 西安西北有色物化探总队有限公司 | 一种铝土矿定位方法和装置 |
CN106990445A (zh) * | 2017-03-08 | 2017-07-28 | 核工业二六大队 | 放射性地质环境现状调查中放射性异常边界的确定方法 |
CN108241180A (zh) * | 2018-01-19 | 2018-07-03 | 五矿勘查开发有限公司 | 离子型稀土矿底板勘查方法 |
CN108459358A (zh) * | 2018-03-12 | 2018-08-28 | 有色金属矿产地质调查中心 | 沉积盆地基底顶面形态恢复和深部矿体定位预测的新方法 |
US10116232B2 (en) * | 2015-12-18 | 2018-10-30 | Schlumberger Technology Corporation | Collocation of radiation generator components for limited-space devices |
CN109270567A (zh) * | 2018-09-11 | 2019-01-25 | 湖南中核环保科技有限公司 | 一种在铀矿勘查中制成探槽γ照射量率等值图的方法 |
-
2019
- 2019-11-15 CN CN201911116736.4A patent/CN110764163B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2539838C1 (ru) * | 2013-07-16 | 2015-01-27 | Федеральное Государственное Унитарное Предприятие "Сибирский Научно-Исследовательский Институт Геологии, Геофизики И Минерального Сырья" | Способ прогноза скрытого оруденения, связанного с гранитоидами |
US10116232B2 (en) * | 2015-12-18 | 2018-10-30 | Schlumberger Technology Corporation | Collocation of radiation generator components for limited-space devices |
CN105700040A (zh) * | 2016-03-05 | 2016-06-22 | 西安西北有色物化探总队有限公司 | 一种铝土矿定位方法和装置 |
CN106990445A (zh) * | 2017-03-08 | 2017-07-28 | 核工业二六大队 | 放射性地质环境现状调查中放射性异常边界的确定方法 |
CN108241180A (zh) * | 2018-01-19 | 2018-07-03 | 五矿勘查开发有限公司 | 离子型稀土矿底板勘查方法 |
CN108459358A (zh) * | 2018-03-12 | 2018-08-28 | 有色金属矿产地质调查中心 | 沉积盆地基底顶面形态恢复和深部矿体定位预测的新方法 |
CN109270567A (zh) * | 2018-09-11 | 2019-01-25 | 湖南中核环保科技有限公司 | 一种在铀矿勘查中制成探槽γ照射量率等值图的方法 |
Non-Patent Citations (2)
Title |
---|
史剑锋 等: "综合物探在柳林铀矿勘查中的应用效果", 《CT理论与应用研究》 * |
张永恒 等: "四川九龙县锂铍稀有金属矿区天然放射性水平评价", 《世界核地质科学》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111398360A (zh) * | 2020-04-21 | 2020-07-10 | 山东大学 | 基于ERT和IP圈定包气带中LNAPLs的污染源区探测方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN110764163B (zh) | 2021-06-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Siler et al. | Three-dimensional geologic mapping to assess geothermal potential: examples from Nevada and Oregon | |
CN110109895A (zh) | 适用于tbm掘进隧道的围岩分级联合预测方法及应用 | |
EP2235568B1 (en) | Geophysical data processing systems | |
AU2013277928B2 (en) | Systems and methods for processing geophysical data | |
CN111045091B (zh) | 一种玄武岩覆盖下古河道的识别定位方法 | |
CN107966739B (zh) | 钍放射性异常的圈定方法 | |
Wahlgren et al. | Geology Laxemar. Site descriptive modelling SDM-Site Laxemar | |
CN108761550A (zh) | 一种南方花岗岩型铀矿找矿方法 | |
CN110764163B (zh) | 基于γ总量及高密度电法测量下的伟晶岩锂矿的圈定方法 | |
Cao et al. | Research on comprehensive detection and visualize of hidden cavity goaf | |
Gao et al. | Site-scale bedrock fracture modeling of a spent fuel reprocessing site based on borehole group in Northwest, China | |
Gunnarsson | 3D modeling in Petrel of geological CO2 storage site | |
CN110454155B (zh) | 一种应用磁化率测井方法确定第四纪地层年代的方法 | |
Hitchon et al. | Dynamic basin analysis: an integrated approach with large data bases | |
Martins et al. | Impact of geophysics in small-scale mining | |
Lane et al. | Improving geological models through statistical integration of borehole data and geologists’ cross-sections | |
Akanbi et al. | 2-D electrical resistivity survey for cassiterite potential mapping in Jos-Bukuru Area, North Central, Nigeria | |
Breytenbach et al. | Application, advantages and limitations of high-density gravimetric surveys compared with three-dimensional geological modelling in dolomite stability investigations | |
CN117251739B (zh) | 用于砂岩型铀资源定量预测的样本集的构建方法 | |
Shaw et al. | Interpretation and modelling: geology, to support the development of an Integrated Site Descriptive Model | |
Coccia et al. | Accurate 3D location of mine induced seismicity in complex near-field underground conditions | |
Wu et al. | Application of 3D Geological Model in Subway Construction in Singapore | |
Baldwin Jr | Paleoseismic investigation of the San Andreas fault on the north coast segment, near Manchester, California | |
YU et al. | For personal use only | |
Bele | 3D Geological Modeling in Mineral Deposits (Copper Ore Body Cases) |
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 |