海洋小比例尺矿产储量计算方法研究

邵长高　梁建　陈宏文　曾文娟

奉新网站建设公司成都创新互联,奉新网站设计制作，有大型网站制作公司丰富经验。已为奉新上千提供企业网站建设服务。企业网站搭建\成都外贸网站制作要多少钱，请找那个售后服务好的奉新做网站的公司定做！

（广州海洋地质调查局 广州 510760）

第一作者简介：邵长高（1983—），男，硕士，主要研究方向：3S技术在资源调查和生态环境动态监测、数字海洋中的应用和开发。E-mail:zkyscg@ yahoo.com.cn。

摘要　传统矿产储量计算模型基于欧式测量，应用于小比例尺海洋矿产储量计算时存在精度差的问题，论文通过对WGS1984投影、墨卡托投影、兰勃托投影以及阿尔伯斯投影等特性的研究，提出将矿产实体进行切片处理，计算切片间矿物实体的体积累加和的方法，实现了海洋小比例尺地图投影下储量的精确量测及体积计算，系统地论述了在不规则地球椭球体下如何实现海洋矿产储量计算，为我国海洋资源探测和军事战略方面提供基础服务。

关键词　海洋量测　地理信息系统　地图投影　储量计算

1　前言

近年来资源勘探已经覆盖大部分陆地区域。越来越多的国家把目光投向海洋。海洋作为一个巨大的能源和资源宝库在国民经济、军事战略等的重要性也日益显现。各个国家竞相制定海洋科技开发规划、战略计划，优先发展海洋新技术[1]。海洋研究成为一个热点，技术的革新也日新月异。

由于海洋是一个大面积的区域，其与陆地的资源勘探技术存在较大区别，尤其在大范围海洋区域的矿产储量计算方面区别甚大。地球是一个不规则的椭球体，采用传统基于平面的欧式测量方法进行小比例尺海洋地图测量时，由于地图投影等方面的原因将会导致变形，严重影响储量计算的精确度[2]。包括欧洲石油勘探组织在内的国内外机构为了消除这种影响建立了一系列的投影转换公式。这些投影转换应用到二维投影当中一定程度上提高了地图量算的精确度。但是对于地球变形引起三维储量计算方面的误差目前并未提供行之有效的方法。本文在前人研究的基础上通过引入基于投影转换的方法，通过对WGS1984投影、墨卡托投影、兰勃托投影以及阿尔伯斯投影等特性的研究，提出将矿产实体进行切片处理，计算切片间矿物实体的体积累加和的方法，实现了海洋小比例尺地图投影下储量的精确量测及体积计算，系统地论述了在不规则地球椭球体下如何实现海洋矿产储量计算，为我国海洋资源探测和军事战略方面提供基础服务。

2　海洋投影概述

我国的海洋基本比例尺地形图中，海区小于1:50万的地形图多用等角正轴圆柱投影，又叫墨卡托投影（Mercator）[1]。现在我国企事业单位科研人员用的海图大部分为墨卡托投影。但是在海洋小比例尺下计算矿物储量时必须消除墨卡托投影引起的地图变形误差。论文引入了阿尔伯斯投影，利用其在投影变换中面积不变的特性计算储量来消除误差。在矿物深度方向上，切片间距离值取深度值的差值。

3　技术路线

海洋大面积矿产实体，跨度大，地图投影变形明显，形状不规则，因此大大增加了计算储量的难度。论文引入切片技术把矿产实体切成实体面，利用切片间实体的累加和计算实体面之间体积的总和即得矿产实体储量。示意图（图1）如下：

图1 矿物实体切片Fig.1 The slice of the mineral reserve

图1中海洋矿物实体被分割为n个切面，切面间体积和相加即为整个实体的体积。当n趋向于无穷大时则与实际体积越接近。n的值取决于实测数据的精度，也就是经纬度和深度的值的精度。

3.1　数据预处理

3.1.1　数据来源

1）多波速水深数据：多波束数据经常应用于湖泊盆地等的体积运算。多波束水深其工作原理是通过声波发射与接收换能器阵进行声波广角度定向发射、接收，通过各种传感器（卫星定位系统、运动传感器、电罗经、声速剖面仪等）对各个波束测点的空间位置归算，从而获取在与航向垂直的条带式高密度水深数据[6]

2）地震剖面数据：海洋矿产储量数据主要来自海洋地震剖面断层数据。地震勘探方法是在地面上布置一条条的测线，沿各条测线进行地震施工采集地震信息，然后经过电子计算机处理就得出一张张地震剖面图。经过地质解释的地震剖面图就像从地面向下切了一刀，在二维空间（长度和深度方向）上显示了地下的地质构造情况（图2）[7]。海洋地震剖面中可以根据断层的层位读取炮点号，并结合导航数据读取矿产储层的坐标数据。

图2 二维地震剖面示意图Fig.1 Two dimensional seismic data

3.1.2　数据入库

从多波束或者地震剖面中提取出的位置数据，数据整理按照如下数据库格式入库：

表1 矿物储量数据结构　Table.1 The data sheet of the mineral reserve

表中数据的经度、纬度需存储经投影转换处理后变成的阿尔伯斯投影数据。

3.2　切面面积计算

3.2.1 切面绘制

运用sql语言搜索深度相同的多边形的边界值，绘制切面。方法为：

1）用sql语言搜索出数据库数据中深度值相同的数据。

2）取所有数据中一个特定数据（a1,b1），此数据需要位于所有坐标值（ax,bx）之间。

3）从（a1,b1）的0度角开始逆时针计算两者之间距离值L=sqrt[（b2-b1）2+（a2-a1）2]。同时计算角度差。如果过角度差相等则取L值较大的点。

4）把所有3）中取出数据连接成多边形即为此切面。

3.2.2　切面计算

为了保持面积计算结果不受地球椭球体影响需要将墨卡托投影转换为阿尔伯斯投影。墨卡托投影转阿尔伯斯投影在ArcEngine下方法如下[4]:

Dim pPoint As esriGeometry.IPoint

Set pPoint=New Point

pPoint.PutCoords mx,my

Set pPoint.SpatialReference=pSpRef2

pPoint.Project pSpRef1‘此处先实现由墨卡托投影到WGS1984投影中

lon=pPoint.X

lat=pPoint.Y

Set pPCS=pSpRFc.CreateProjectedCoordinateSystem（esriSRProjCS_NAD1983USA_Albers）

Set pSpRef2=pPCS

pPoint.Project pSpRef2‘实现由WGS1984投影到阿尔伯斯投影的转换

lon=pPoint.X‘lon即为在阿尔伯斯投影中的经度值

lat=pPoint.Y‘lat即为阿尔伯斯投影中的纬度值[4’

ArcEngine是目前地理信息系统处理方面比较流行的二次开发工具。墨卡托投影转化为阿尔伯斯投影时，每一个坐标点均要做转换，通常是采用W GS1984投影作为中间转换投影。先将墨卡托投影转化到WGS1984投影，然后将转化来的WGS1984投影转化成阿尔伯斯投影。

阿尔伯斯投影最大的特点是投影前后面积保持不变，本文采用质心量算法进行面积计算，具体步骤是先寻找多边形的质心，然后由质心到各多边形顶点引直线，最后把每个多边形的面积相加即得结果。计算步骤如图3[4]。

方法为[4]:

1）首先遍历数据库，读取数据库中高程相等数据的坐标值组成平面多边形。找出多边形质心。

2）连接多边形每个点与质心。

3）计算每个小多边形的面积然后相加。S=s1+s2+s3………。其中S表示多边形面积，s1、s2、s3等表示小三角形面积[4]。

设L为边长，L两端点坐标值为（a1,b1），（a2,b2）。如图4所示：

则：L=sqrt[（b2-b1）2＋（a2-a1）2]

每个小三角形面积计算源代码为[4]:

s=（L1+L2+L3）/2

S=sqrt[s＊（s-L1））＊（s-L2）＊（s-L3）]

图3 多边形的面积量算[4]Fig.3 Area measurement of the polygon

图4 每个小三角形面积计算Fig.4 The calculation of every triangle

此处S值即为切面面积。切面面积的计算结果考虑了地球椭球体引起的误差更接近实际值。

3.3　切面间体积计算

将矿物实体分割切片后其中每个切面间体积v的计算模拟梯形计算模式，S上为上切面面积，S下为下切面面积，h为切面间高度差。如图5所示：

图5 单个切面实体Fig.5 Single slice object

则切面间体积v＝（S上＋S下）h/2。图1和5中当切面数n趋向于无穷大时，切面1和切面2之间的面积差值越小，相应的两个多边形的形状也就最接近，h值也就最小。此时可以得到误差较为小的切面体积计算结果。

3.4　矿物储量计算

将矿产实体分割成n个切面后，每个小切面的体积的累加和即为整个矿产实体的储量。切面数n的值越大所切割的体积个数越多，则切面值越接近实际值。体积值V即是每个小切面间体积v的累加和。

南海地质研究.2010

式中：V即为整个矿物储量。它累加了所有的切面间实体的体积之和，切面间实体的个数取决于n的大小。当n趋向于无穷大时最接近实际值。

4　结语

本文介绍了基于投影转换的海洋小比例尺矿产储量的计算方法，同时提供了基于Ar-cEngine的投影转换方法。矿产储量的计算模式不同于传统的计算模式，关键在于考虑到了小比例尺下由于地球椭球体变形引起的误差。所以论文引入了投影变换的方法，从一定程度上降低了地球的不规则性引起的误差。但是此方法只适应于固体矿产的储量计算，对于石油、水合物等的储量计算只能做体积计算的一个参数。

参考文献

[1]单宝强，毛永强.GIS中的坐标系定义与转换[J].黑龙江国土资源，2005,11,38～39

[2]欧洲石油勘探组.Coordinate Conversions and Transformation including Formulas[M].国际石油技术软件开放公司，2008

[3]苏国辉，戴勤奋，魏合龙.海洋地质数据库数据的存储结构[J].海洋地质动态，2003,19（6）:5～7

[4]邵长高，谭建军等.海洋小比例尺地图精确测量及计算方法[J].地理与地理信息科学，2009,25（2）:42～45

[5]

[6]

[7]

Method of Precise Measurem ent and Calculation of Small Scale Mineral Reserve Calculation

Shao Changgao,Liang Jian,Chen Hongwen,Zeng Wenjuan

（Guangzhou Marine Geological Survey,Guangzhou,510760）

Abstract:To the small-scale map in ocean mine reserve field,the traditional measurement method computes the reserve with a relatively coarse precision.In order to improve that,a new method has been provided in this study,which uses Arc Engine technology to finish the conversion between different projections and measure the earth's area as well as other information precisely.And then cut the mine reserve object into several pieces,so we can calculate the volume of the reserve by summing every piece.The different projections,such as WGS1984,Mercator,and Albers,also have been discussed,which can provide a good service for the military strategy and exploration of ocean resources.

Key words:Ocean measurements,GIS,Map projection,Reserve Calculation

矿产资源储量空间数据库建设与维护中MAPGIS 应用实践

蒋洪明 李雪洁 肖 茜

（江苏省地质资料馆）

摘 要 利用馆藏地质资料信息对全省矿产资源储量空间数据库进行实时更新维护，是馆藏地质资料服务经济社会的有效途径之一。为了提高建库与维护工作的质量、效率以及加强技术交流，作者对江苏矿产资源储量空间数据库建设应用过程中，地质资料发挥的重要作用，以及利用 MAPGIS 技术采集矿体资源储量估算范围、水平投影图形拐点坐标的方法和步骤进行了总结阐述。

关键词 地质资料 MAPGIS 矿产资源 空间数据库

随着地质资料信息化的发展，对地质资料信息需求也在不断地提高。储量空间数据库建设与数据维护工作，是我省地质资料信息资源管理的一项重要内容。在促进地质资料深化利用过程中，地质资料为矿产资源储量空间数据库建设与维护，提供了重要的基础信息资源。通过地质资料提供的可靠数据信息来源，我馆以 MAPGIS 技术为基础，采用创新方法以及现代化信息系统技术，使矿产资源储量空间数据库数据得到了不断补充、更新与完善，实现了矿产资源储量空间数据库的动态管理，同时，也进一步加强了地质资料信息集群化产业化的发展。

1 储量空间数据库的建设与维护

随着信息技术在国土资源领域的应用发展，2003 年 10 月 14 日，国土资源部办公厅下发了《关于建立矿产资源储量空间数据库的通知》（国土资厅发 [2003]324 号），要求各省国土资源行政主管部门开展矿产资源储量空间数据库建设与日常维护工作，带动矿产资源储量科学化、规范化、标准化、图形化管理发展，促进矿产资源储量管理由传统的“以数管矿”向“以图管矿”新模式的转变，增加矿产资源储量信息对探矿权、采矿权和矿产资源规划等矿政管理工作的支撑力，进一步提高矿政管理现代化水平。此项工作的开展，拓展了馆藏地质资料服务经济社会的新途径。

1.1 矿产资源储量空间数据库建设

根据国土资源部下发的矿产资源储量空间数据库建设工作技术要求，矿产资源储量空间数据库建设，是通过提取地质资料信息，进行矿区矿产特征以及资源储量估算范围水平投影图形拐点坐标数据等相关信息的采集。

储量空间数据库建库的流程为：确定入库矿产地；收集矿区资料；矿区资料预处理分析，确定作为矿产资源储量估算水平投影的图件；圈定矿体资源储量估算范围水平投影图形，量算拐点坐标；采集矿区地质特征等相关信息；坐标数据与信息采集质量检查与订正；数据录入；数据录入质量检查与订正；提供利用。

1.2 矿产资源储量空间数据库维护

矿产资源储量空间数据库维护工作就是对新发现的矿产地，按建库工作流程，采集坐标数据与信息，填加入库；同时对已提供利用的矿产资源储量空间数据库在应用过程中发现的问题，进行检查、核实与更正，从而增强矿产资源储量空间数据库的现势性、准确性。

2 矿产资源储量空间数据库维护现状

矿产资源储量空间数据库建设与维护是一项长期、持续的工作。在工作应用过程中，发现存在着以下两个方面的问题。

一是矿产资源储量数据库中的矿产基本特征信息的数据采集，来源于汇交至我馆的地质资料报告中的数据，但在实际的工作过程中，仍然有大量的成果地质资料报告并未汇交至我馆。

二是在初期建库阶段，矿体资源储量范围水平投影图形拐点坐标，采用了人工直接在纸质矿区图件上丈量、读取数据、填写数据采集表；这一方法的弊端为采集的数据精度比较低，而且差错率也较高，大大增加了数据录入及校对工作的工作量，从而降低了工作效率。

3 MAPGIS 技术在矿产资源储量空间数据库建设中的辅助作用

针对储量空间数据库应用中存在的问题，我馆总结经验，研究方法，从而采取了相应的解决措施。矿产资源储量空间数据库建设是“以图管矿”思想的体现，其技术关键点是矿体资源储量估算范围水平投影图形科学、合理、准确地圈定，拐点坐标精准地量算。从建库工作流程和此项工作的技术关键点可以看出，矿体资源储量估算范围水平投影图形的拐点坐标精准量算工作环节，存在探索研究的空间与必要性。

为了提高矿体资源储量估算范围水平投影图形拐点坐标数据采集的精度与工作效率，江苏在矿产资源储量空间数据库建设过程中进行了探索研究，引入 GIS 技术，实现计算机自动、精准采集坐标数据，很大程度上提高了拐点坐标数据采集精度与工作效率，创建了矿体资源储量范围水平投影图形拐点数据采集新方法。

由于矿体是成矿地质作用形成的可开发利用的自然地质体，且形态多样，千差万别，其水平投影范围的平面形状亦是多样、各异。因此，人工量算其坐标点数据，费时、费力，且精度与效率都低。应用MAPGIS 软件的强大功能，可精确、快捷地采集拐点数据，并能自动生成坐标点数据。具体工作步骤如下。

3.1 MAPGIS 生成标准图框

打开 MAPGIS 投影变换模块，利用键盘生成矩形图框，根据需要填写上相应的参数。坐标系通常填写“国家坐标系”，它的起始代号根据原地质图相应的带号填写；标注选择为“公里值”；网起始值是公里网从哪点开始，以光栅图内图廓左下角 X 及 Y 值作为起始公里值，以内图廓右上角 X 及 Y 坐标作为终止公里值，单位选择“公里”。根据原图比例尺，填写网格间距，例如，原图比例尺为 1∶2000，则网格间距为 0.2，网格线类型选择“绘制实线坐标线”，通过以上设置，将生成所需要的标准图框。

3.2 生成 MSI 影像文件

在 MAPGIS 主界面上，通过图像处理模块，打开图像分析，在转换数据类型中选择要转换的光栅文件的类型（如JPG、TIF、BMP等），选择需要转换的光栅文件，然后指定转换后的MSI影像文件存放目录。在影像文件中，选定生成的 MSI 影像文件，打开“镶嵌融合”功能，选择“参照点 / 线 / 区文件”。为保证配准的精度，选取多个控制点，进行配准，每确定一个控制点后，影像图下方控制点 ID 就会依次出现控制点的 X、Y 坐标，然后通过校正预览，查看校正后的影像是否准确，影像完成校正后，对形成的校正影像文件进行保存。

3.3 拐点投影及坐标转换

在 MAPGIS 中，进入实用服务模块，选择投影变换下级目录用户数据点文件投影转换，将已知的坐标点文件（TXT 文本文件）导入进来。通过“用户投影参数”设置，输入投影相应的参数，根据要求，设置相应的文件类型。按要求完成“结果投影参数”的设置。其中，注意比例尺按需要投影的比例尺确定，投影带类型及序号应根据坐标确定，对“点图元参数”进行“子图号”的设定。设置完成后，将点保存，即形成了矿体资源储量估算范围水平投影图形拐点（含标志点）坐标表，见表 1。

新建一个MAPGIS项目工程，新建一个点文件，打开点编辑，根据属性标注释。然后再新建一个线文件，将点文件按属性标上的序号顺序，进行连接，按住 Ctr 键，点鼠标右键，线将自动闭合。通过其他文字及图的编辑功能，形成矿体资源储量估算范围水平投影图形，见图 1。至此，矿体资源储量估算范围水平投影图形拐点坐标计算机软件化采集工作完毕。

表 1 矿体资源储量估算范围水平投影图形拐点坐标表

图 1 矿体资源储量估算范围水平投影图

4 地质资料在矿产资源储量空间库中的作用

从矿产资源储量空间库建设与维护过程可以看出，矿产资源储量空间数据库建设核心任务是采集矿体资源储量估算范围水平投影范围的拐点坐标和矿区地质特征等相关信息进行入库。这两方面的重点信息均蕴藏在矿区地质勘查报告中。因此，收集到合格的矿区地质勘查报告是矿产资源储量空间数据库建设的前提条件。馆藏的矿区地质勘查报告均为通过评审、正式印刷汇交的地质报告；且通过了汇交验收，报告质量与权威性毋庸置疑，其文字内容的齐全性、数据的准确性、图件的精确性，为矿产资源储量空间库的建设提供了数据与信息来源与质量保障。

5 结语

利用馆藏地质资料，通过 GIS 技术在矿产资源储量空间数据库建设中的创新应用，江苏于 2004 年率先建设了全省矿产资源储量数据库，并提供全省矿政管理和基本建设等领域利用。同时根据有关要求，在已建数据库的基础上，不断补充、维护、完善，实现对空间数据库的动态管理。江苏地质资料馆持续应用在实践中自己创建的矿体资源储量范围水平投影图形拐点坐标采集方法，按年度维护与实时维护相结合的方式，开展全省矿产资源储量空间数据库维护工作。截至 2011 年底，全省矿产资源储量空间数据库中入库的矿区达 499 个，涉及 616 个矿产地，较好地维护了江苏矿产地信息的现势性与准确性，为全省矿业权设置、矿产资源规划、建设项目压矿审批等矿政管理工作提供了有力的支撑服务。

用于矿床资源储量计算的地理信息系统

地理信息系统与计算机建模技术被矿业企业、地质机构、咨询公司、设计与研究机构所广泛使用。

以下是开发用于地质勘查活动以及矿产储量的计算与开采的计算机软件领域的全球领导者：

1）AQUILA Mining Systems Ltd.（加拿大蒙特利尔）：该公司是地质勘查软件开发领域的先驱，是世界上第一家研制和实施能独立识别岩石及其主要特性的钻探系统的公司。该公司目前是卡特彼勒公司（Gaterpillar Gorporation）旗下的一员；

2）DataMine（Datamine国际公司，英国）：该软件产品主要被咨询公司及矿业企业用作不同种类矿产的综合采矿管理系统。它是世界上最强大的多功能系统之一；

3）Vulcan（由KJRA系统公司开发）是足够强大的一体化系统，由大量模块组成，用以解决各种地质、采矿、测量、环境目标等；

4）MineScape系统（以及Ellipse、MineStar、MineMarket程序）由澳大利亚Mincom PtyLtd.开发。MineScape被开发各类固体矿产的企业所采用；

5）Lynx系统（由南非Lynx Geosystems S.A.（Pty）Ltd.开发）由一套尖端地压仪器组成，主要设计用于解决环境问题；

6）MineSight系统（由Mintec，Inc.开发）被美国和加拿大的矿业企业广泛采用。它拥有强大的地质信息处理、矿藏建模与矿产储量估算能力；

7）Gemcom系统（由加拿大Gemcom Software International Inc.开发）可运行一个完整的工作循环，最后创建出矿藏模型，完成露天矿和地下矿的设计与规划；

8）Micromine系统（由澳大利亚Micromine Pty开发）：这款软件产品包含许多用于解决地质勘查活动及矿床作业过程中所出现的问题的工具；

9）Techbase套件（由美国Minsoft Ltd.开发）主要开发用于地质方面，适合用于矿藏、地质图及储量估算的3D建模；

10）Geostat软件包（由加拿大Geostat Systems International Inc.开发）被用于矿藏建模、储量估算与采矿作业规划；

11）GDM一体化系统（由法国地质调查局BRDM开发）：该系统的主要用户包括油气田开发企业。系统拥有尖端的地压零部件及灵活的绘图功能。它具有强大的地球物理与地球化学数据处理与分析能力；

12）PetrelE＆P软件平台一体化系统（由国际公司斯伦贝谢开发）。该系统的目的是解决油气田储量建模与计算的各种相关问题。系统拥有用于任何类型油气田建模的强大数学工具、地球物理测量数据处理能力及其他邻域功能；

13）MineFrame地理信息系统（由俄罗斯科学院柯拉科学中心矿业大学开发）设计用于各种地质目标的综合性解决方案及矿藏（主要是固体矿产）的3D建模，以满足矿业企业、科研和设计机构的需求；

14）地质统计软件工具（由Volodymyr Maltsev开发）设计用于解决相关地质统计问题及矿产储量的计算；

15）K-MINE一体化地理信息系统（由乌克兰KRYVBASAKADEMINVEST开发）设计用于解决各类矿藏建模的各种相关问题。所创建的模型可在矿产开发阶段被矿业企业的作业部门所采用。

大部分计算机程序（系统）是集地质环境数据收集、处理与传输软件和设备于一体的信息分析套件。它们可完成地质构造与矿山巷道的3D显示、矿产储量的计算、地质勘查与采矿作业的规划与优化、环保措施的实施等工作。

开发和改迸地质信息解释方法作为矿产储量计算的基础，考虑了矿产资源储量分类的需求，和为便于利用计算机技术迸行矿床建模和储量计算充分发挥各信息价值，而满足改变传统（前苏联时期广泛应用，独联体国家到现在仍保留）的地质勘查信息采集、解释分析、综合整理方法的条件。

矿床计算机建模和地质统计学方法的使用，能精确地反映出矿床内矿化（流体饱和度）参数的空间分布规律，因为它们能考虑影响储量计算的全套指标。采用假设的几何形状和勘查网的密度，在此基础上构建的块模型及动力学模型能最准确地反映出矿体和构造的天然各向异性。

计算机建模的使用能帮助评价和圈出矿产资源储量的空间形态，区分不同类型和不同商业（技术）品位的矿产的储量/资源。

创建3D矿床模型的方法主要取决于矿床构造及矿产类型。建模过程具体包括以下七个步骤[1]。

1）数据库结构的开发，以储存地质勘查获取的原始信息；

2）地质成果数据的录入和分析：编辑待录入系统的地质信息；地质试验、地球物理测量信息录入；原始地质数据的统计分析，查证（纠错），数据分组，数据库整理，找出规律性；

3）地质勘查信息解释、矿床建模：模型空间内钻孔构建，按剖面线分组；根据地层和岩性定义和圈定矿与非矿层段，根据边界品位调整层段（地质数据解释）；考虑到构造变形，根据地球物理测量信息（地震、电法、磁法和重力测量），清楚划分岩石空间分界线；

4）创建地质体的线框模型：线框矿床建模（矿体和围岩、地层、异常、构造（traps）等的建模）；

5）地质统计实地考察：对空间数据、变异性、各组分地质特性的空间变异（各向异性）规律迸行地质统计学分析；水动力系统建模，计算迁移、污染、化学成分等的；

6）建立矿床的块模型：创建空块模型；采用数学方法（最近地区（多边形）法、距离反比加权插值（IDW）法、克里格法（正在修改）等迸行各成分含量的插值；根据预设的成矿条件，调整某一矿床的岩石分布等高线；矿产资源储量级别和类别的定义；

7）评价和估算储量：确定矿物成分（原矿状态）的最低边界品位；定义储量类别和级别；

不同类型矿产的实地建模技术基本相同，已投入运营的矿产地建模技术略有不同（对于这类矿产地，通常已编制了采矿图形文件（平面图、剖面图、地质图），以便根据生产勘探、采样测试以及实际的矿床开发数据，调整空间岩石分布等高线）。在勘查信息解释和储量计算阶段，不同类型矿床的模型构建步骤区别明显。

钻孔位置选择、合理确定勘查工程网度、试验测试方法和质量[5，11]，应在实地建模与地质数据解释之前迸行。