CN110674606A - 一种计算矿体产状的方法 - Google Patents

一种计算矿体产状的方法 Download PDF

Info

Publication number
CN110674606A
CN110674606A CN201910958770.XA CN201910958770A CN110674606A CN 110674606 A CN110674606 A CN 110674606A CN 201910958770 A CN201910958770 A CN 201910958770A CN 110674606 A CN110674606 A CN 110674606A
Authority
CN
China
Prior art keywords
ore body
thickness
ore
line
inclination angle
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
Application number
CN201910958770.XA
Other languages
English (en)
Other versions
CN110674606B (zh
Inventor
唐忠伟
甘怀营
施耘
崔丽娜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hunan Changye Construction Engineering Construction Drawing Review Co ltd
Original Assignee
CINF Engineering Corp Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by CINF Engineering Corp Ltd filed Critical CINF Engineering Corp Ltd
Priority to CN201910958770.XA priority Critical patent/CN110674606B/zh
Publication of CN110674606A publication Critical patent/CN110674606A/zh
Application granted granted Critical
Publication of CN110674606B publication Critical patent/CN110674606B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/02—Agriculture; Fishing; Mining

Abstract

本发明公开了一种计算矿体产状的方法，可自动计算矿体产状数据的方法，解决矿体产状数据计算过程繁琐、工作量大、设计效率低的问题。本发明的方法计算矿体的产状数据非常简单，可大量节省计算时间；同时在计算面积时，其结果更为精确。本发明的方法中结合CAD图元反应器，可方便数据的修改。

Description

一种计算矿体产状的方法
技术领域
本发明属于矿山资源评价、矿山储量计算工程领域，具体涉及一种计算矿体产状的方法。
背景技术
采矿设计过程中，地质将地质资源数据和矿体各勘探线剖面线提交给采矿专业。采矿专业根据矿山具体条件划分中段后，根据矿体类型和条件进行采矿方法设计。在这个过程中需要计算矿体产状，需要计算的矿体产状的信息包括：矿体平均倾角、矿体平均厚度和各矿体类型所占的比例。目前，采矿工程师计算矿体产状的步骤为：在CAD软件中，在勘探线剖面图上切分出单元矿体，根据单元矿体轮廓线，人工计算出单元矿体倾角、厚度和面积，然后将这些结果依次按矿体类型誊抄到电子表格文件中，再在电子表格文件中编辑统计计算公式，最终得到各类型矿体的平均倾角、平均厚度及其所占的比例。采用人工计算与统计的方式存在如下缺点：1)计算数据量大、过程繁琐、效率低，在一个项目中要计算得到矿体产状结果，一般需要3到5个工作日时间；2)在计算单元矿体面积时，仅得出单元矿体面积的近似值，结果不够准确；3)计算结束后，数据结果的检验、修改极不方便。同时，目前市场上还没有一款能够有效解决矿体产状计算与统计的矿业软件或程序。
发明内容
本发明的目的是提供一种可自动计算矿体产状的方法，解决矿体产状计算过程繁琐、工作量大、设计效率低的问题。
本发明这种计算矿体产状的方法，包括以下步骤：
1)根据获取各勘探线剖面图上的各矿体的矿体剖面轮廓线，建立各勘探线剖面图上的各矿体的矿体剖面轮廓线集
Figure BDA0002228242180000011
根据设计给出的各中段标高值，建立标高直线集H(h1,…,hk,…,hl)；
其中：Si第i个矿体在各勘探线上的剖面集合，i∈[1,m]，m为矿体数量，
Figure BDA0002228242180000012
为第i个矿体在第j条勘探线上的剖面轮廓线，j∈[1,n]，n为勘探线数量，i、j、m、n都为正整数；hk为第k条标高直线，k∈[1,l]，l为标高线的条数，k、l都为正整数；
2)依次用标高直线集中两条相邻的标高线hk、hk+1分别切割所述Si集合中各
Figure BDA0002228242180000021
轮廓线，得到单元矿体轮廓线集
Figure BDA0002228242180000022
其中：U为单元矿体轮廓线构成的集合，
Figure BDA0002228242180000023
为第i个矿体在第j条勘探线剖面上被相邻的两条标高线hk、hk+1切割形成的轮廓线；
3)依次根据所述每个
Figure BDA0002228242180000024
轮廓线的拐点坐标和第i个矿体的下盘方向，分别绘制每个单元矿体的倾角图元和厚度图元；
4)依次计算所述每个
Figure BDA0002228242180000025
轮廓线的倾角厚度
Figure BDA0002228242180000027
面积
Figure BDA0002228242180000028
并将所述的倾角、厚度、面积值保存在该所述的
Figure BDA0002228242180000029
轮廓线图元的图元词典中；
5)编辑CAD图元反应器，所述反应器内容为：当厚度图元或倾角图元变化时，程序将自动重新计算厚度图元和倾角图元，并同时更新相应
Figure BDA00022282421800000210
轮廓线图元的图元词典数据，以便用户检查修改各单元矿体的产状数据值(在此步骤中，CAD图元反应器用于CAD软件的二次开发中，当某个图元被修改时，自动响应某个函数，如例为：当厚度图元或倾角图元被修改后，软件自动重新计算和编辑厚度图元和倾角图元)。
6)在CAD图形数据库中依次检索所述各矿体轮廓线Si上的
Figure BDA00022282421800000211
轮廓线图元，并从该图元的图元词典中获取
Figure BDA00022282421800000212
轮廓线的倾角、厚度、面积值，然后按矿体厚度类型分别统计出所述各矿体的平均倾角、平均厚度、所占比例；
7)在CAD图形数据库中依次检索所述各矿体轮廓线Si上的
Figure BDA00022282421800000213
轮廓线图元，并从该图元的图元词典中获取
Figure BDA00022282421800000214
轮廓线的倾角、厚度、面积值，按矿体倾角类型分别统计出所述各矿体的平均倾角、平均厚度、所占比例；
8)根据所述各矿体的计算统计结果，输出所述各矿体的产状数据表。
所述步骤1)中，标高直线为水平直线，标高直线的Y坐标值等于各中段标高值的Y坐标；标高直线集H(h1,…,hk,…,hl)中，hk按标高直线的Y坐标值从小到大的排列形成的集合。
所述步骤2)中，获得单元矿体轮廓线
Figure BDA0002228242180000031
的具体步骤如下：
2.1分别计算标高线hk、hk+1与
Figure BDA0002228242180000032
轮廓线的交点，并把所述交点插入到对应的轮廓线的拐点坐标序列中；
2.2遍历所述步骤2.1得到的新
Figure BDA0002228242180000039
轮廓线的拐点坐标，依次提取出Y坐标值不大于所述标高线hk+1的Y坐标值且不小于所述标高线hk的Y坐标值的拐点坐标集合，按集合坐标依次绘制多段线，形成的图元即为所述单元矿体轮廓线
Figure BDA0002228242180000034
所述步骤3)中，绘制每个矿体倾角图元和矿体厚度图元的具体步骤为：
3.1根据所述
Figure BDA0002228242180000035
轮廓线，得到所述轮廓线上Y坐标值等于所述hk的Y坐标值的的拐点坐标集合；
3.2如果步骤3.1所述的拐点坐标集合不为空，则设置P1为所述拐点坐标集合中X坐标值最小的拐点，设置P2为所述拐点坐标集合中X坐标值最大的拐点；如果步骤3.1所述的拐点坐标集合为空，则设置P1、P2都为轮廓线拐点坐标中Y坐标值最小的拐点；
3.3根据所述
Figure BDA0002228242180000037
轮廓线，得到所述轮廓线上Y坐标值等于所述hk+1的Y坐标值的的拐点坐标集合；
3.4如果步骤3.3所述的拐点坐标集合不为空，则设置P3为所述拐点坐标集合中X坐标值最小的拐点，设置P4为所述拐点坐标集合中X坐标值最大的拐点；如果步骤3.3所述的拐点坐标集合为空，则设置P3、P4都为
Figure BDA0002228242180000038
轮廓线拐点坐标中Y坐标值最大的拐点；
3.5如果所述下盘方向为“右边”，则以P2为起点、以P4为终点连接成线段图元，即为倾角图元；如果所述下盘方向为“左边”，则以P1为起点、以P3为终点连接成线段图元，即为倾角图元；
3.6根据步骤3.5得到的所述倾角图元，得到倾角图元的中心点，得到倾角图元的中垂线与所述
Figure BDA0002228242180000041
轮廓的交点集合，以所述中心点为起点，以所述交点集合中距离所述中心点最远的点为终点，连接成线段图元，即为厚度图元。
所述步骤4)中，如果下盘方向为“右边”，则所述倾角为所述倾角图元的角度值
Figure BDA0002228242180000042
如果下盘方向为“左边”，则所述倾角为180°减去所述倾角图元的角度值所得的结果
Figure BDA0002228242180000043
所述厚度为所述厚度图元的长度值
Figure BDA0002228242180000044
所述面积为所述
Figure BDA0002228242180000045
轮廓线所围成区域的面积值
Figure BDA0002228242180000046
所述步骤6)中，矿体的厚度类型分为薄矿体、厚矿体、极厚矿体；其划分标准为：薄矿体的矿体厚度不大于3米；厚矿体的矿体厚度在3米到15米之间；极厚矿体的矿体厚度在15米以上。
所述步骤6)中，所述按矿体厚度类型分别统计出所述第i个矿体的平均倾角、平均厚度、所占比例的计算，包括以下步骤：
6.1将矿体进行分类：按照α＝1为薄矿体类型，α＝2为厚矿体类型，α＝3为极厚矿体类型，接着统计α＝1时，各倾角
Figure BDA0002228242180000047
厚度
Figure BDA0002228242180000048
和面积
Figure BDA0002228242180000049
的值，如果计算过程中
Figure BDA00022282421800000410
属于α＝2或3矿体类型，则各倾角
Figure BDA00022282421800000411
厚度
Figure BDA00022282421800000412
和面积
Figure BDA00022282421800000413
记为0；然后按照α＝1分类计算方法，统计α＝2和3的各倾角
Figure BDA00022282421800000414
厚度
Figure BDA00022282421800000415
和面积
Figure BDA00022282421800000416
的值；具体计算公式如公式(1)所示；
6.2将6.1步骤中的
Figure BDA0002228242180000052
和
Figure BDA0002228242180000053
按照公式(2)中的公式进行计算，得到α类型矿体的平均角度α类型矿体的平均厚度
Figure BDA0002228242180000055
α类型矿体所占的比例ηα，矿体的平均角度
Figure BDA0002228242180000056
矿体的平均厚度
其中：Sak、Aak、Bak、Sa均为计算过程中的中间参数。
所述步骤(7)中，矿体倾斜类型分为缓倾斜矿体、倾斜矿体、极倾斜矿体；其中：缓倾斜矿体的矿体倾角在3度到30度之间；倾斜矿体的矿体倾角在30度到50度之间；极倾斜矿体的矿体倾角大于50度。
所述步骤7)中，按矿体倾角类型分别统计出所述各矿体的平均倾角、平均厚度、所占比例和面积，具体包括以下步骤：
7.1将矿体进行分类：按照β＝1为缓倾斜矿体，β＝2为倾斜矿体，β＝3为极倾斜矿体，接着统计β＝1时，各倾角
Figure BDA0002228242180000062
厚度
Figure BDA0002228242180000063
和面积
Figure BDA0002228242180000064
的值，如果计算过程中
Figure BDA0002228242180000065
属于β＝2或3矿体类型，则各倾角
Figure BDA0002228242180000066
厚度
Figure BDA0002228242180000067
和面积
Figure BDA0002228242180000068
记为0；然后按照β＝1分类计算方法，统计β＝2和3的各倾角
Figure BDA0002228242180000071
厚度
Figure BDA0002228242180000072
和面积
Figure BDA0002228242180000073
的值；具体计算公式如公式(3)所示；
Figure BDA0002228242180000074
7.2将7.1步骤中的和
Figure BDA0002228242180000076
按照公式(4)中的公式进行计算，得到β类型矿体的平均角度
Figure BDA0002228242180000077
β类型矿体的平均厚度β类型矿体所占的比例ηβ；
Figure BDA0002228242180000081
其中：Sβk、Aβk、Bβk、Sβ均为计算过程中的中间参数。
本发明的有益效果：1)本发明的方法计算矿体的产状数据非常简单，可大量节省计算时间；同时在计算面积时，其结果更为精确。2)本发明的方法中结合CAD图元反应器，可方便数据的修改。
附图说明
图1为本发明实施例提供的一种矿体产状自动计算与统计方法流程图；
图2为本发明实施例提供的一种相邻标高直线hk、hk+1切分第i个矿体在第j个勘探线剖面的轮廓线得到单元矿体
Figure BDA0002228242180000082
的示意图；
图3(a)为本发明实施例提供的一种下盘方向为“右边”时的单元矿体的倾角图元和矿体厚度图元绘制示意图；(b)为本发明实施例提供的一种下盘方向为“左边”时的单元矿体的倾角图元和矿体厚度图元绘制示意图；
图4为本发明实施例提供的一种单元矿体的面积计算区域示意图；
图5为本发明实施例提供的一种某个矿体产状统计计算结果图；
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合发明实施例中的附图，对本发明实施例中的技术方案进行清楚的描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例1
本实施例基于CAD软件及其二次开发语言，结合办公软件，完成矿体产状自动计算与统计系统的开发。使用该系统可在几个小时内快速完成矿体产状的计算与统计，大大提高了数据准确性和设计效率，其具体的流程如图1所示，包括以下步骤：
1、根据获取各勘探线剖面图上的各矿体的矿体剖面轮廓线，建立各勘探线剖面图上的各矿体的矿体剖面轮廓线集
Figure BDA0002228242180000091
根据设计给出的各中段标高值，建立标高直线集H(h1,…,hk,…,hl)；标高直线为水平直线，标高直线的Y坐标值等于各中段标高值的Y坐标，标高直线集中的hk按标高直线的Y坐标值按从小到大在集合中排列。
其中：Si第i个矿体在各勘探线上的剖面集合，i∈[1,m]，m为矿体数量，
Figure BDA0002228242180000092
为第i个矿体在第j条勘探线上的剖面轮廓线，j∈[1,n]，n为勘探线数量，i、j、m、n都为正整数；hk为第k条标高直线，k∈[1,l]，l为标高线的条数，k、l都为正整数。
2、依次用标高值线集中两条相邻的标高线hk、hk+1分别切割所述Si集合中各
Figure BDA0002228242180000093
轮廓线，如图2所示，分别计算标高线hk、hk+1与
Figure BDA0002228242180000094
轮廓线的交点，并把所述交点插入到对应的轮廓线的拐点坐标序列中；遍历所得到的新
Figure BDA0002228242180000102
轮廓线的拐点坐标，依次提取出Y坐标值且不大于标高线hk+1的Y坐标值且不小于所述标高线hk的Y坐标值的拐点坐标集合，按集合坐标依次绘制多段线，形成的图元即为所述单元矿体轮廓线
Figure BDA0002228242180000103
最后将单元矿体轮廓线整理成集合，即得到单元矿体轮廓线集
Figure BDA0002228242180000105
其中：U为单元矿体轮廓线构成的集合，
Figure BDA0002228242180000106
为第i个矿体在第j条勘探线剖面上被相邻的两条标高线hk、hk+1切割形成的轮廓线。
3、依次根据所述每个
Figure BDA0002228242180000107
轮廓线的拐点坐标和第i个矿体的下盘方向，分别绘制每个单元矿体的倾角图元和厚度图元，具体步骤为：
3.1根据所述
Figure BDA0002228242180000108
轮廓线，得到所述轮廓线上Y坐标值等于所述hk的Y坐标值的拐点坐标集合。
3.2如果步骤3.1所述的拐点坐标集合不为空，则设置P1为所述拐点坐标集合中X坐标值最小的拐点，设置P2为所述拐点坐标集合中X坐标值最大的拐点；如果步骤3.1所述的拐点坐标集合为空，则设置P1、P2都为
Figure BDA0002228242180000109
轮廓线拐点坐标中Y坐标值最小的拐点。
3.3根据所述
Figure BDA00022282421800001010
轮廓线，得到所述轮廓线上Y坐标值等于所述hk+1的Y坐标值的的拐点坐标集合。
3.4如果步骤3.3所述的拐点坐标集合不为空，则设置P3为所述拐点坐标集合中X坐标值最小的拐点，设置P4为所述拐点坐标集合中X坐标值最大的拐点；如果步骤3.3所述的拐点坐标集合为空，则设置P3、P4都为
Figure BDA00022282421800001011
轮廓线拐点坐标中Y坐标值最大的拐点。
3.5如果所述下盘方向为“右边”，则以P2为起点、以P4为终点连接成线段图元，即为倾角图元；如果所述下盘方向为“左边”，则以P1为起点、以P3为终点连接成线段图元，即为倾角图元；如图3a所示。
3.6根据步骤3.5得到的所述倾角图元，得到倾角图元的中心点，得到倾角图元的中垂线与所述
Figure BDA0002228242180000111
轮廓的交点集合，以所述中心点为起点，以所述交点集合中距离所述中心点最远的点为终点，连接成线段图元，即为厚度图元如图3b所示。
4、依次计算所述每个
Figure BDA0002228242180000112
轮廓线的倾角
Figure BDA0002228242180000113
厚度
Figure BDA0002228242180000114
面积
Figure BDA0002228242180000115
并将所述的倾角、厚度、面积值保存在该所述的
Figure BDA0002228242180000116
轮廓线图元的图元词典中；
其中：如果下盘方向为“右边”，则所述倾角为所述倾角图元的角度值
Figure BDA0002228242180000117
(如图3a所示)；如果下盘方向为“左边”，则所述倾角为180°减去所述倾角图元的角度值所得的结果
Figure BDA0002228242180000118
(如图3b所示)；所述厚度为所述厚度图元的长度值
Figure BDA0002228242180000119
所述面积为所述
Figure BDA00022282421800001110
轮廓线所围成区域的面积值
Figure BDA00022282421800001111
(如图4所示)。
5、编辑CAD图元反应器，所述反应器内容为：当厚度图元或倾角图元变化时，程序将自动重新计算厚度图元和倾角图元，并同时更新相应
Figure BDA00022282421800001112
轮廓线图元的图元词典数据，以便用户检查修改各单元矿体的产状数据值(在此步骤中，CAD图元反应器用于CAD软件的二次开发中，当某个图元被修改时，自动响应某个函数；如例为：当厚度图元或倾角图元被修改后，软件自动重新计算和编辑厚度图元和倾角图元)。
6、在CAD图形数据库中依次检索所述各矿体轮廓线Si上的轮廓线图元，并从该图元的图元词典中获取
Figure BDA00022282421800001114
轮廓线的倾角、厚度、面积值，然后按矿体厚度类型分别统计出所述各矿体的平均倾角、平均厚度、所占比例；其中：矿体的厚度类型分为薄矿体、厚矿体、极厚矿体；其划分标准为：薄矿体的矿体厚度不大于3米；厚矿体的矿体厚度在3米到15米之间；极厚矿体的矿体厚度在15米以上；
具体方法包括以下步骤：
6.1、将矿体进行分类：按照α＝1为薄矿体类型，α＝2为厚矿体类型，α＝3为极厚矿体类型，接着统计α＝1时，各倾角
Figure BDA0002228242180000121
厚度
Figure BDA0002228242180000122
和面积
Figure BDA0002228242180000123
的值，如果计算过程中属于α＝2或3矿体类型，则各倾角
Figure BDA0002228242180000125
厚度
Figure BDA0002228242180000126
和面积
Figure BDA0002228242180000127
记为0；然后按照α＝1分类计算方法，统计α＝2和3的各倾角
Figure BDA0002228242180000128
厚度
Figure BDA0002228242180000129
和面积
Figure BDA00022282421800001210
的值；具体计算公式如公式(1)所示；
6.2、将a3步骤中的
Figure BDA00022282421800001212
和
Figure BDA00022282421800001213
按照公式(2)中的公式进行计算，得到α类型矿体的平均角度
Figure BDA00022282421800001214
α类型矿体的平均厚度
Figure BDA00022282421800001215
α类型矿体所占的比例ηα，矿体的平均角度
Figure BDA00022282421800001216
矿体的平均厚度
Figure BDA0002228242180000131
其中：Sak、Aak、Bak、Sa均为计算过程中的中间参数。
7、在CAD图形数据库中依次检索所述各矿体轮廓线Si上的轮廓线图元，并从该图元的图元词典中获取
Figure BDA0002228242180000133
轮廓线的倾角、厚度、面积值，按矿体倾角类型分别统计出所述各矿体的平均倾角、平均厚度、所占比例；其中：矿体倾斜类型分为缓倾斜矿体、倾斜矿体、极倾斜矿体；其中：缓倾斜矿体的矿体倾角在3度到30度之间；倾斜矿体的矿体倾角在30度到50度之间；极倾斜矿体的矿体倾角大于50度。
具体方法包括以下步骤：
7.1:按照β＝1为缓倾斜矿体，β＝2为倾斜矿体，β＝3为极倾斜矿体，接着统计β＝1时，各倾角
Figure BDA0002228242180000141
厚度和面积的值，如果计算过程中
Figure BDA0002228242180000144
属于β＝2或3矿体类型，则各倾角
Figure BDA0002228242180000145
厚度
Figure BDA0002228242180000146
和面积
Figure BDA0002228242180000147
记为0；然后按照β＝1分类计算方法，统计β＝2和3的各倾角
Figure BDA0002228242180000148
厚度
Figure BDA0002228242180000149
和面积
Figure BDA00022282421800001410
的值；具体计算公式如公式(3)所示；
Figure BDA00022282421800001411
7.2：将7.1步骤中的和
Figure BDA00022282421800001413
按照公式(4)中的公式进行计算，得到β类型矿体的平均角度
Figure BDA00022282421800001414
β类型矿体的平均厚度
Figure BDA00022282421800001415
β类型矿体所占的比例ηβ。
Figure BDA0002228242180000151
其中，Sβk、Aβk、Bβk、Sβ均为计算过程中的中间参数。
8、根据所述各矿体的计算统计结果，并将所述各矿体的产状数据表输出到电子表格文件中，如图5所示。

Claims (9)

1.一种计算矿体产状的方法，包括以下步骤：
1)根据获取各勘探线剖面图上的各矿体的矿体剖面轮廓线，建立各勘探线剖面图上的各矿体的矿体剖面轮廓线集
Figure FDA0002228242170000011
根据设计给出的各中段标高值，建立标高直线集H(h1,…,hk,…,hl)；
其中：Si第i个矿体在各勘探线上的剖面集合，i∈[1,m]，m为矿体数量，
Figure FDA0002228242170000012
为第i个矿体在第j条勘探线上的剖面轮廓线，j∈[1,n]，n为勘探线数量，i、j、m、n都为正整数；hk为第k条标高直线，k∈[1,l]，l为标高线的条数，k、l都为正整数；
2)依次用标高直线集中两条相邻的标高线hk、hk+1分别切割所述Si集合中各
Figure FDA0002228242170000013
轮廓线，得到单元矿体轮廓线集
Figure FDA0002228242170000014
其中：U为单元矿体轮廓线构成的集合，
Figure FDA0002228242170000015
为第i个矿体在第j条勘探线剖面上被相邻的两条标高线hk、hk+1切割形成的轮廓线；
3)依次根据所述每个
Figure FDA0002228242170000016
轮廓线的拐点坐标和第i个矿体的下盘方向，分别绘制每个单元矿体的倾角图元和厚度图元；
4)依次计算所述每个轮廓线的倾角厚度
Figure FDA0002228242170000019
面积并将所述的倾角、厚度、面积值保存在该所述的
Figure FDA00022282421700000111
轮廓线图元的图元词典中；
5)编辑CAD图元反应器，所述反应器内容为：当厚度图元或倾角图元变化时，程序将自动重新计算厚度图元和倾角图元，并同时更新相应
Figure FDA00022282421700000112
轮廓线图元的图元词典数据，以便用户检查修改各单元矿体的产状数据；
6)在CAD图形数据库中依次检索所述各矿体轮廓线Si上的
Figure FDA00022282421700000113
轮廓线图元，并从该图元的图元词典中获取
Figure FDA00022282421700000114
轮廓线的倾角、厚度、面积值，然后按矿体厚度类型分别统计出所述各矿体的平均倾角、平均厚度、所占比例；
7)在CAD图形数据库中依次检索所述各矿体轮廓线Si上的
Figure FDA00022282421700000115
轮廓线图元，并从该图元的图元词典中获取
Figure FDA0002228242170000021
轮廓线的倾角、厚度、面积值，按矿体倾角类型分别统计出所述各矿体的平均倾角、平均厚度、所占比例；
8)根据所述各矿体的计算统计结果，输出所述各矿体的产状表。
2.根据权利要求1所述的计算矿体产状数据的方法，其特征在于，所述步骤1)中，标高直线为水平直线，标高直线的Y坐标值等于各中段标高值的Y坐标；标高直线集H(h1,…,hk,…,hl)中，hk按标高直线的Y坐标值从小到大的排列形成的集合。
3.根据权利要求1或2所述的计算矿体产状的方法，其特征在于，所述步骤2)中，获得单元矿体轮廓线
Figure FDA0002228242170000022
的具体步骤如下：
2.1分别计算标高线hk、hk+1与
Figure FDA0002228242170000023
轮廓线的交点，并把所述交点插入到对应的
Figure FDA0002228242170000024
轮廓线的拐点坐标序列中；
2.2遍历所述步骤2.1得到的新
Figure FDA0002228242170000025
轮廓线的拐点坐标，依次提取出Y坐标值不大于所述标高线hk+1的Y坐标值且不小于所述标高线hk的Y坐标值的拐点坐标集合，按集合坐标依次绘制多段线，形成的图元即为所述单元矿体轮廓线
Figure FDA0002228242170000026
4.根据权利要求1或2所述的计算矿体产状数据的方法，其特征在于，所述步骤3)中，绘制每个矿体倾角图元和矿体厚度图元的具体步骤为：
3.1根据所述
Figure FDA0002228242170000027
轮廓线，得到所述轮廓线上Y坐标值等于所述hk的Y坐标值的的拐点坐标集合；
3.2如果步骤3.1所述的拐点坐标集合不为空，则设置P1为所述拐点坐标集合中X坐标值最小的拐点，设置P2为所述拐点坐标集合中X坐标值最大的拐点；如果步骤3.1所述的拐点坐标集合为空，则设置P1、P2都为
Figure FDA0002228242170000028
轮廓线拐点坐标中Y坐标值最小的拐点；
3.3根据所述
Figure FDA0002228242170000029
轮廓线，得到所述轮廓线上Y坐标值等于所述hk+1的Y坐标值的的拐点坐标集合；
3.4如果步骤3.3所述的拐点坐标集合不为空，则设置P3为所述拐点坐标集合中X坐标值最小的拐点，设置P4为所述拐点坐标集合中X坐标值最大的拐点；如果步骤3.3所述的拐点坐标集合为空，则设置P3、P4都为
Figure FDA0002228242170000031
轮廓线拐点坐标中Y坐标值最大的拐点；
3.5如果所述下盘方向为“右边”，则以P2为起点、以P4为终点连接成线段图元，即为倾角图元；如果所述下盘方向为“左边”，则以P1为起点、以P3为终点连接成线段图元，即为倾角图元；
3.6根据步骤3.5得到的所述倾角图元，得到倾角图元的中心点，得到倾角图元的中垂线与所述
Figure FDA0002228242170000032
轮廓的交点集合，以所述中心点为起点，以所述交点集合中距离所述中心点最远的点为终点，连接成线段图元，即为厚度图元。
5.根据权利要求1所述的计算矿体产状数据的方法，其特征在于，所述步骤4)中，如果下盘方向为“右边”，则所述倾角为所述倾角图元的角度值
Figure FDA0002228242170000033
如果下盘方向为“左边”，则所述倾角为180°减去所述倾角图元的角度值所得的结果
Figure FDA0002228242170000034
所述厚度为所述厚度图元的长度值
Figure FDA0002228242170000035
所述面积为所述
Figure FDA0002228242170000036
轮廓线所围成区域的面积值
Figure FDA0002228242170000037
6.根据权利要求1所述的计算矿体产状数据的方法，其特征在于，所述步骤6)中，矿体的厚度类型分为薄矿体、厚矿体、极厚矿体；其划分标准为：薄矿体的矿体厚度不大于3米；厚矿体的矿体厚度在3米到15米之间；极厚矿体的矿体厚度在15米以上。
7.根据权利要求6所述的计算矿体产状数据的方法，其特征在于，所述步骤6)中，所述按矿体厚度类型分别统计出所述第i个矿体的平均倾角、平均厚度、所占比例的计算，包括以下步骤：
6.1将矿体进行分类：按照α＝1为薄矿体类型，α＝2为厚矿体类型，α＝3为极厚矿体类型，接着统计α＝1时，各倾角
Figure FDA0002228242170000041
厚度
Figure FDA0002228242170000042
和面积的值，如果计算过程中
Figure FDA0002228242170000044
属于α＝2或3矿体类型，则各倾角厚度
Figure FDA0002228242170000046
和面积
Figure FDA0002228242170000047
记为0；然后按照α＝1分类计算方法，统计α＝2和3的各倾角
Figure FDA0002228242170000048
厚度
Figure FDA0002228242170000049
和面积
Figure FDA00022282421700000410
的值；具体计算公式如公式(1)所示；
Figure FDA00022282421700000411
6.2将6.1步骤中的
Figure FDA00022282421700000412
和按照公式(2)中的公式进行计算，得到α类型矿体的平均角度
Figure FDA00022282421700000414
α类型矿体的平均厚度
Figure FDA00022282421700000415
α类型矿体所占的比例ηα，矿体的平均角度
Figure FDA00022282421700000416
矿体的平均厚度
Figure FDA00022282421700000417
Figure FDA0002228242170000051
其中：Sak、Aak、Bak、Sa均为计算过程中的中间参数。
8.根据权利要求1所述的计算矿体产状数据的方法，其特征在于，所述步骤(7)中，矿体倾斜类型分为缓倾斜矿体、倾斜矿体、极倾斜矿体；其中：缓倾斜矿体的矿体倾角在3度到30度之间；倾斜矿体的矿体倾角在30度到50度之间；极倾斜矿体的矿体倾角大于50度。
9.根据权利要求8所述的计算矿体产状数据的方法，其特征在于，所述步骤7)中，按矿体倾角类型分别统计出所述各矿体的平均倾角、平均厚度、所占比例，具体包括以下步骤：
7.1按照β＝1为缓倾斜矿体，β＝2为倾斜矿体，β＝3为极倾斜矿体，接着统计β＝1时，各倾角
Figure FDA0002228242170000061
厚度
Figure FDA0002228242170000062
和面积
Figure FDA0002228242170000063
的值，如果计算过程中
Figure FDA0002228242170000064
属于β＝2或3矿体类型，则各倾角
Figure FDA0002228242170000065
厚度
Figure FDA0002228242170000066
和面积记为0；然后按照β＝1分类计算方法，统计β＝2和3的各倾角
Figure FDA0002228242170000068
厚度
Figure FDA0002228242170000069
和面积
Figure FDA00022282421700000610
的值；具体计算公式如公式(3)所示；
Figure FDA00022282421700000611
7.2将7.1步骤中的
Figure FDA00022282421700000612
和
Figure FDA00022282421700000613
按照公式(4)中的公式进行计算，得到β类型矿体的平均角度
Figure FDA00022282421700000614
β类型矿体的平均厚度
Figure FDA00022282421700000615
β类型矿体所占的比例ηβ；
Figure FDA0002228242170000071
其中：Sβk、Aβk、Bβk、Sβ均为计算过程中的中间参数。
CN201910958770.XA 2019-10-10 2019-10-10 一种计算矿体产状的方法 Active CN110674606B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910958770.XA CN110674606B (zh) 2019-10-10 2019-10-10 一种计算矿体产状的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910958770.XA CN110674606B (zh) 2019-10-10 2019-10-10 一种计算矿体产状的方法

Publications (2)

Publication Number Publication Date
CN110674606A true CN110674606A (zh) 2020-01-10
CN110674606B CN110674606B (zh) 2022-10-14

Family

ID=69081559

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910958770.XA Active CN110674606B (zh) 2019-10-10 2019-10-10 一种计算矿体产状的方法

Country Status (1)

Country Link
CN (1) CN110674606B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103824329A (zh) * 2014-01-20 2014-05-28 中国地质科学院矿产资源研究所 一种地质勘探三维可视化储量估算方法
CN109356585A (zh) * 2018-11-16 2019-02-19 青海省第三地质矿产勘查院 一种判断金属硫化物矿体产状的方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103824329A (zh) * 2014-01-20 2014-05-28 中国地质科学院矿产资源研究所 一种地质勘探三维可视化储量估算方法
CN109356585A (zh) * 2018-11-16 2019-02-19 青海省第三地质矿产勘查院 一种判断金属硫化物矿体产状的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
VIKTOR A. KHAKULOV 等: ""Creation of Automated Lump Sorting of Molybdenum-Containing Ores in the Technology of Forming the Quality of Ore Mass Flows at the Mining Stage"", 《2018 IEEE INTERNATIONAL CONFERENCE "QUALITY MANAGEMENT, TRANSPORT AND INFORMATION SECURITY, INFORMATION TECHNOLOGIES" (IT&QM&IS)》 *
徐帅等: "基于多方法联合的露天转地下境界矿柱厚度优化", 《东北大学学报(自然科学版)》 *
罗晓峰: "微机AUTO CAD在矿山地质设计及其储量计算中的应用", 《新疆有色金属》 *

Also Published As

Publication number Publication date
CN110674606B (zh) 2022-10-14

Similar Documents

Publication Publication Date Title
CN101673410B (zh) 一种基于矢量建筑图纸重建三维模型的方法
Farouki The characterization of parametric surface sections
CN102136218B (zh) 一种计算机地图制图中的地图注记自动配置方法
CN110633518B (zh) 基于bim的简易参数化单元式幕墙材质分类统计方法
CN110765523A (zh) 一种基于bim技术的深基坑围护结构快速构建方法
CN111310278A (zh) 一种基于仿真的船舶自动化建模方法
CN103150428A (zh) 一种基于地质信息的地下洞室自动建模方法
CN103150328A (zh) 一种利用数据库自动制图的集成方法
CN101930483A (zh) 应用参数化设计模型化简数字地图居民地多边形的方法
CN109918760A (zh) 基于Revit互用性的钢筋建模方法
CN102360436B (zh) 一种基于部件的联机手写藏文字符的识别方法
CN106844966B (zh) 一种螺旋桨叶面叶背精确建模方法
CN110674606B (zh) 一种计算矿体产状的方法
CN103106313B (zh) 轧后件顺序重构方法
US7225396B2 (en) System and method for generating a report on an object
CN108595792B (zh) 一种飞机油箱燃油测量传感器布局优化方法
CN108121853A (zh) 一种基于AutoCAD计算挖填方工程量的系统与方法
CN116485102A (zh) 基于bim技术的建筑表面积材料量精确算量方法
CN113820756B (zh) 实时日照检测方法及系统
CN115690340A (zh) 一种基于二、三维实时联动的交叉地质剖面界线调整方法
CN113888486A (zh) 一种基于三维模型辅助矢量化、辅助dlg质检的方法
CN115578530A (zh) 基于Dynamo创建三维隧道地质模型的方法
CN114463494A (zh) 一种地形特征线自动提取算法
JP2002259473A (ja) ３次元ソリッドモデルのモデリングに必要な基礎デ−タの受け渡し方式
CN108197360B (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
TR01 Transfer of patent right

Effective date of registration: 20240417

Address after: 410000, No. 199 Jiefang Middle Road, Furong District, Changsha City, Hunan Province

Patentee after: Hunan Changye Construction Engineering Construction Drawing Review Co.,Ltd.

Country or region after: China

Address before: 410019 No. 299, Mu Lin East Road, Yuhua District, Changsha, Hunan

Patentee before: CHANGSHA ENGINEERING & RESEARCH INSTITUTE OF NONFERROUS METALLURGY Co.,Ltd.

Country or region before: China