DZ／T 0338.3-2020 标准规范下载简介
DZ／T 0338.3-2020 固体矿产资源量估算规程 第3部分：地质统计学法.pdfDZ/T0338.3—2020目次前言范围2规范性引用文件，3术语和定义·4基本原理4.1克里格法原理·4.2距离幂次反比法原理5适用条件及方法选择估算参数选择.7数据准备8估算流程8.1地质解译建立矿化域或矿体模型.·8.2估算域的划分。8.3样品数据统计分析及特异值处理8.4样品等长度组合8.5建立块模型8.6建立变异函数模型·8.7搜索椭球体的设置LC8.8交叉验证118.9变量估值8.10变量估值结果验证119资源量分类1210资源量估算结果12附录A（资料性附录）建立数据库基础表13附录B（资料性附录)三维图件的编制15附录C（资料性附录)矿块指标体系及圈矿原则17附录D（资料性附录）特异值处理合理性检验18附录E（资料性附录）对数正态克里格法和指示克里格法的实施步骤19附录F（资料性附录）地质统计学法资源量估算成果20
本规程根据GB/T1.1一2009《标准化工作导则第1部分：标准的结构和编写》给出的规则起草。 DZ/T0338《固体矿产资源量估算规程》分为四个部分： 一第1部分：通则； 一第2部分：几何法； 一第3部分：地质统计学法； 一第4部分：SD法。 本部分为DZ/T0338的第3部分。 本部分由中华人民共和国自然资源部提出。 本部分由全国自然资源与国土空间规划标准化技术委员会（SAC/TC93)归口。 本部分主要起草单位：自然资源部矿产资源储量评审中心、北京科技大学、北京东澳达科技有限 公司。 本部分主要起草人：张树泉、胡建明、冯涛、高利民、赵婷钰、唐长钟
DZ/T0338.3—2020固体矿产资源量估算规程第3部分：地质统计学法1范围DZ/T0338的本部分规定了地质统计学法估算固体矿产资源量的基本原理、适用条件及方法选择、估算参数选择、数据准备、估算流程、资源量分类及估算成果。本部分适用于勘查程度达到详查阶段及以上要求的固体矿产勘查、矿山设计、开发阶段的资源量估算工作。2规范性引用文件下列文件对于本文件的应用是必不可少的。凡是注日期的引用文件，仅注日期的版本适用于本文件。凡是不注日期的引用文件，其最新版本（包括所有的修改单）适用于本文件。GB/T13908固体矿产地质勘查规范总则GB/T17766固体矿产资源储量分类GB/T33444固体矿产勘查工作规范3术语和定义下列术语和定义适用于本文件。3. 1矿化域domain矿化域是基于地质控矿因素（如岩性、地层、构造、围岩蚀变等）或矿化边界或将二者综合考虑圈连的，一个比较连续的矿化地质体。3. 2区域化变量regionalizedvariable区域化变量是一种限制在一个特定的空间内（如一个矿体内或矿化域内）的变量，既有结构性（相关性）又有随机性的特征。如在特定的空间内的样品品位之间是相关的，可以参与矿体品位的估算。超出这个空间范围的样品，由于样品品位之间相关性减弱，随机性的成分增大，一般不参与矿体品位的估算。3. 3变异函数variogram变异函数是地质统计学的应用工具，它能够反映区域化变量的空间变化特征一一相关性和随机性，反映局部范围和特定方向上地质特征的变化。3. 4实验变异函数experimental variogram变异函数对任一给定矿床2是未知的，需要通过样品值对之进行估计。设从矿床Q中获得一组样1
DZ/T0338.3—20203. 11带状异向性zonal anisotropy带状异向性与几何异向性相对应。根据变异函数曲线，区域化变量在不同方向上表现出总基台值不同，而不论变程相同或不相同，这种异向性均称为带状异向性，如图4所示。y(h)C+C.C0R,R,h图4带状异向性示意图3.12克里格法kriging克里格法是地质统计学中一种局部估计方法，以变异函数为主要工具，对区域化变量进行插值，求插值过程中的最优线性无偏估计量，再通过块体约束计算资源量的方法。根据研究目的和条件不同，分别有简单克里格法（simplekriging）、普通克里格法（ordinarykriging）、对数正态克里格法（lognormalordinarykriging）、指示克里格法（indicatorkriging）和泛克里格法（krigingwithtrend）等。其中应用最为广泛的是普通克里格法。3. 13距离幕次反比法inversedistanceweight距离幂次反比法是一种与空间距离有关的插值方法，即在估计待估点的值时，按照距离越近权重值越大的原则，利用已知点和待估点之间的距离取幕次后的倒数为权系数进行加权平均YD／T 1095-2018 通信用交流不间断电源（UPS），再通过块体药束计算资源量的方法。3. 14克里格方差krigingvariance采用克里格法估计时，估计误差的方差称为克里格方差，用来衡量估计的精度。4基本原理4.1克里格法原理应用克里格法对区域化变量进行局部估计时，将矿体划分成许多相同或相似的长方体，用在一定范围内系列样品点的品位值"。（α=1，2，，n）为每个长方体估值，以估计方差最小、权系数之和等于1为4
DZ/T0338.3—20205.5当区域化变量变化系数满足矿化在比较均匀的范围内（对应矿化均匀或比较均匀的，具体变化系数值可查阅各分矿种规范，一般变化系数小于或等于150%），可用距离幂次反比法或普通克里格法。6估算参数选择6.1克里格法估算参数包括变异函数参数、搜索椭球体参数和块体参数；距离幂次反比法估算参数包括幕次参数、搜索椭球体参数和块体参数。6.2变异函数参数包括：a）基本滞后距：是计算变异函数时分隔样品对的矢量的最小长度，一般由计算全向变异函数获得。当选定一个基本滞后距得到的全向实验变异函数最容易拟合出理论变异函数时，即为最佳滞后距。b)块金值（C。）：反映出区域化变量在小尺度上的变异程度，表明变量的随机成分。沿钻孔方向以最小滞后距计算变异函数时，所拟合的理论变异函数曲线与纵坐标的交点即为块金值，见图5。一般情况下，一个估算域内只有一个块金值。基台值（C）：表明了估算域内变量的相关程度，其数值由变量的方差减去块金值取得。由基台值与块金值构成了变异函数的总基台（C十C），见图5。d)变程（R）：表明了变量在估算域内自相关性存在的最大范围，一般是理论变异函数曲线中最初达到总基台值时所对应的距离，见图5DZy(h)C+C,R图5变异函数球状模型曲线示意图6.3幂次参数的选取：a）幂次一般在13之间取值。变量的空间变化性越大、变化速度越快，则取值越大。b)根据工程网度和样品密度取值。密度大，取大值，反之取小值。c）应通过交叉验证幂次的合理性。6.4搜索椭球体参数包括椭球体的方位角、倾伏角、倾角、主轴的搜索半径、次轴的搜索半径（或用主轴与次轴的比值表示）、短轴的搜索半径（或用主轴与短轴的比值表示）。若样品点空间分布不均匀，具有丛聚的情况，可将搜索椭球体分为多个扇区，每个扇区内可根据需要确定最少工程数、最少样品数和最多样品数。6.5块体模型（也称品位模型）的参数主要是设置块模型方向和块体尺寸（长、宽、高）以及进行次级分块，见图6。块体模型属性可记录块体所在的位置（即块质心点坐标）、尺寸大小，还可以附加不同属性如矿岩类型、矿种、品位、体积质量、数学运算赋值的属性、资源类别等。这些属性均保存6
DZ/T0338.3—2020在块模型文件中，既可以通过估值的方法赋值，也可以通过属性名称直接赋值，还可以进行相关数学运算赋值。a)地质解释b)块模型图6块体模型示意图6.6常用的克里格法及距离幂次反比法估算参数均包括椭球体参数、最多样品数、最少样品数、最少工程数，变异函数参数只在各克里格法中涉及，详见表1。表1不同估算方法参数表估值参数方法球体参数最多样品数最少样品数最少工程数变异函数参数附加特征参数简单克里格法V均值参数普通克里格法V指示克里格法指示值对数克里格法VV对数变换参数泛克里格法VJ漂移参数距离幂次反比法幕幂次V注：/表示应实际填人的具体参数值。数据准备7.1建立数据库，从地质数据库或电子表格中提取参与资源量估算的基础数据，至少包括勘查工程或生产工程定位表、工程测斜表、样品分析表和地质岩性表，参见附录A。表中所列的字段是必不可少的。7.2对数据进行纠错及完整性和逻辑性检查。7.3通过三维软件对工程数据进行位置对比查看和空间关系分析。8估算流程8.1地质解译建立矿化域或矿体模型8.1.1将地质数据库中的钻孔、探槽和坑道等工程数据绘制成图，作剖面图或中段平面图，参见附录B。矿岩界线清晰时，在剖面图或中段平面图上根据品位结合地层、岩性、构造，以及矿体产状等地质特征进7
DZ/T0338.3—2020行地质解译，圈定矿体或矿化域及夹石边界线。8.1.2根据矿化特征，矿岩界线呈过渡关系时，选取矿化品位值或低于一般工业指标的值圈连矿化域，采用矿块指标体系圈定矿体范围，参见附录C。8.1.3剖面图上矿体或矿化域边界线的圈定和外推原则，参见本规程通则部分和有关固体矿产勘查工作规范。8.1.4进行资源量估算通常建立的模型包括地形三维模型、矿体或矿化域模型、夹石模型、岩体模型和地质构造模型等。8.1.5采用矿体模型或矿化域模型作为估算范围时，需要对模型进行质量检查，确保模型的合理性。8.2估算域的划分8.2.1当矿体被断层、岩脉等线性构造切割或错断（见图7)，使得构造两侧区域化变量统计特征存在明显差异，应以此线性构造为界划分不同的区域，设置不同的估算参数进行估算。8.2.2当矿体出现褶皱，使得矿体的走向和倾角发生变化（见图8），无法使用一个搜索椭球体时，应将矿体划分成不同的估算域，分别设置不同的估算参数进行估算，或用动态椭球参数进行估值。图7矿体被断层错断示意图图8矿体的产状受褶皱影响发生改变示意图8.2.3当矿体内出现多个矿石类型或者有用元素局部富集（见图9)，或者受岩性控制（见图10），而使区域化变量统计特征存在明显差异时，应将矿体按不同矿石类型或富集特征划分成不同的区域，分别设置不同的估算参数进行估算。贫矿体富矿体图9根据矿体的富集特征不同划分估算域示意图图10根据岩性不同划分估算域示意图8.3样品数据统计分析及特异值处理8.3.1对样品统计分析及特异值处理时，应以矿化域或矿体为单位进行。8
DZ/T0338.3—20208.3.2对区域化变量进行估值前，分别绘制区域化变量的分布直方图、累计频率分布曲线或概率图等，判断其分布特征，并计算区域化变量的统计特征，如最小值、最大值、样品数、平均值、中值、方差、标准差、偏度、峰度和变化系数等。8.3.3识别和处理特异值可采用分位数法、估计邻域法、影响系数法、概率曲线法和累计频率分布曲线法等数理统计方法。特异值处理后的检验参见附录D。8.3.4如采用对数正态克里格法或指示克里格法进行变量估值时，不进行特异值的处理。8.4样品等长度组合8.4.1样品等长度加权组合应限制在相应的矿体（层)内。其长度取值应为统计中主要的样长值。8.4.2组合后的样段与原样段的长度变量统计特征应保持最大程度上的一致。8.4.3等长组合须在对原样品进行特异值处理后进行。8.5建立块模型8.5.1根据矿化域或矿体空间范围、形态和产状，采用一定长、宽、高尺寸的块为单位，形成整个估算区域的块模型（可根据矿体产状和形态创建次级分块和旋转块模型）。8.5.2影响块模型尺寸的主要因素有矿体规模、勘查工程间距、开采方式（露天开采或地下开采）和采矿工艺（矿房尺寸或露天开采台阶高度）、克里格效率等。其中克里格效率三（块方差一克里格方差）/块方差，其值介于0～1之间。通过比对不同尺寸的块模型的克里格效率值，克里格效率值最大时所对应的块模型尺寸是最优的。8.6建立变异函数模型8.6.1通过调整滞后距的大小，确定最佳基本滞后距。8.6.2通过计算沿钻孔小尺度上样品值的变异性，确定块金值。8.6.3通过计算各方向的变异函数，确定估算域的最大连续性方向（主轴）、次连续性方向（次轴）和最小连续性方向（短轴），三个方向在空间上应相互垂直。求得三个方向的实验变异函数（见图11），采用适当的数学模型进行拟合，求出每一个方向的基台值C、变程R。20.0短轴17.5 主轴15.0 12.510.0A7.5 5.0 2.5 255075100125150175hh一—距离，单位为米（m）。下同。图11三个方向的实验变异函数示意图9
DZ/T 0338.3—2020图15搜索椭球体与矿体或估算域叠加示意图间距，也可以通过变异函数模型确定。8.7.5采用可变椭球体进行估值时，应使用实体模型对搜索椭球体参数进行约束。8.8交叉验证8.8.1交叉验证的目的是检验变异函数模型参数和搜索椭球体参数的合理性。8.8.2交叉验证的理想结果是真值与估计值的误差均值趋近于0”，误差方差与克里格估计方差比值趋近于“1”。8.8.3若交叉验证结果不理想，应修改或重新拟合理论变异函数并再次交叉验证，直到取得一套合理的变异函数参数和搜索椭球体参数为止。8.9变量估值8.9.1对每一个矿化域或矿体，可根据区域化变量的分布和统计特征，选择普通克里格法或距离幂次反比法。其他克里格法如对数正态克里格法和指示克里格法的实施步骤参见附录E。8.9.2采用普通克里格法估算时，对个别矿体规模较小且计算变异函数有困难的情况下，可采用距离幂次反比法或其他方法进行估算。8.9.3使用距离幂次反比法进行变量估值时，应根据工程间距、样品数和估算域规模确定估值搜索半径；使用克里格法时，估值搜索半径的选择还须主要考虑变程。8.9.4对于采用工程指标体系圈定的矿体，最大搜索半径可超过变程，确保矿体内所有块体都能获得相应的估计值。采用矿块指标体系时，估值的最大搜索半径一般不超过变程的4倍。8.10变量估值结果验证8.10.1采用地质统计学法（含距离幂次反比法）对区域化变量进行估计后，应采用局部验证或全局验证的方法进行检验。8.10.2局部验证的方法是在剖面图或平面图上比较参与估算的组合样品值与相邻块体估计值之间的相近程度、外推的距离和方向，判断估值的可靠性。8.10.3全局验证的方法是沿某些方向按照一定的间距划分一系列的区域，分别计算落在每一个区域中参与估值的组合样品平均值和块体模型估值的平均值，绘制成两条曲线图，通过两条曲线的吻合度，判断估值的合理性和可靠性。8.10.4可采用另一种适用的估算方法对主要矿体进行验算，对比绝对误差和相对误差，判断估值参数的可靠性。11
9.1应根据搜索 空间上相互穿插或包含，应根据工程控制程度和地质可靠程度进行人为干预。 9.2应在剖面图上或平面图上勾画出不同类型轮廊界线并圈连成分类轮廓体，对块体模型的分类属性 进行赋值。 9.3采用我国的资源量分类规范进行分类
0.1资源量汇芯参 天规范行 10.2地质统计学法提交资 相关报告编写规范
附A （资料性附录） 建立数据库基础表
根据地质工程数据NB／T 42121-2017 火电机组辅变压器机变频器低电压穿越技术规范，建立地质工程数据库时的基础表有工程定位表（见表A.1）、工程测斜表 A.2)、样品分析表（见表A.3)和地质岩性表（见表A.4）。可结合实际需求增加相应的表格和字段
表A.1为工程定位表。
表A.2为工程测斜表
表A.3为样品分析表。
表A.4为地质岩性表
DZ/T0338.3—2020附录B（资料性附录）三维图件的编制B.1资源量计算的基本图件B.1.1按照DZ/T0078等规范要求，一般先制作地质勘查所需要的矿区（床）地形地质图、勘查线剖面图、水平断面图、中段地质图等。B.1.2根据GB/T33444、DZ/T0079等规范要求，制作与资源量估算有关的矿体（层）顶底板等高线图（或含资源量估算图）、资源量估算图（水平或垂直纵投影图）、砂矿或缓倾斜矿体（层）需要的顶底板等高线和矿层等厚线图（或含资源量估算图）等。B.2三维断面图件的制作B.2.1在实际工作中，通常是按照二维平面图件的方式完成相应的断面图，包括平面图和剖面图。这些图件的坐标均为标注值而不是实际的位置，比如，水平断面图中所有的标高均为0。在资源量估算中，如果采用三维软件按照实际剖面位置和工程坐标及测斜数据进行基本图元（工程、矿段、构造等）的制作，更有利于矿体的对应分析、块段划分和资源量估算。B.2.2制作时，可以按照勘查线或水平中段分别进行切制，形成相应的断面图位置。再利用三维软件中工程的位置、轨迹和矿层信息，在一定的断面范围内，按照实际的位置保留工程的位置并形成点或者线文件进行保存，从而形成了一个断面、一定范围内的三维断面图件的所有制图元素（见图B.1，图B.2)。单位为米削面3、削面4剖面2剖面1高程200图B.1实际位置剖面图图B.2实际高程的中段图B.3三维纵投影图的制作B.3.1矿体纵投影图分为水平投影图和垂直纵投影图。水平投影图可以投影到任意平面，但垂直纵投影图是可以实现与矿体走向相一致的三维纵投影图（见图B.3）。15
DZ/T0338.3—2020B.3.2制作三维纵投影图的步骤：首先按照基线（垂直勘查线)方向作纵剖面，其走向与矿体走向一致，但与矿体有一定的间距；其次将地形线、矿岩界线、探矿工程点、构造界线、块段界线等投影到该剖面上，并进行块段划分和资源量估算，同时需要保存这些文件。图B.3三维纵投影示意图16
SY／T 0089-2019 油气厂、站、库给水排水设计规范附录C （资料性附录） 矿块指标体系及圈矿原则
矿块指标体系一般根据地质矿化规律采用某一个品位（一般介于地质上的矿化品位与常用指标体系 中的边界品位之间)圈出的一个比较完整的矿化域，在矿化域内按照一定的大小划分估计品位的单元块， 继而对单元块进行品位估值，再采用边际品位界定单元块是矿石还是废石，然后统计资源量。 当品位分布均匀，矿体与围岩又是渐变关系时，采用矿块指标体系的优点更突出；相反条件下，采用 工程指标体系，有着显著的优点。