0　引言

区域矿产资源可供性评价，是指在某一确定价格状态下对区域矿产资源开发利用的经济状况进行分析测评，以确定地区矿产资源储量、生产能力及矿产品产量可供应程度的分析评价方法。由此可知，区域矿产资源可供性评价是以地区已探明的矿产资源储量为基础，以当前或未来一定时期内矿产品市场某一价格水平为依据，通过合理的技术经济分析评价，对本地区矿产资源储量转化为矿产品的现实可供应能力和水平所作出的分析与判断。在当今世界形势风云变幻、国际矿产品市场剧烈动荡的背景下，开展矿产资源可供性研究和风险分析，满足国家和地区矿产资源战略、规划及相关政策体系的制定，对保障国家资源安全和实现经济持续稳定发展具有重要的现实意义。

关于矿产资源可供性研究，Mackay et al.(2001)通过建立供给和需求模型，预测了2010年英国油气资源的可供性；Sohn(2006)研究了全球对镁、铜和铁等金属矿产的供需情况，分析了收入、技术、公共政策等因素对需求的影响；Mennenga et al.(2012)建立一个定量预测模型，对澳大利亚不可再生资源的可供性进行了预测；曾澜等(1995)、章少华等(2003)、李海洋(2012)、李裕伟(2015)分别从市场供需条件出发，对不同类型矿产资源的可供性进行了分析。在蒙特卡洛风险分析方面，Bastante et al.(2008)将价格和成本引入分析模型，显示如何集成地质、技术和经济等不确定性因素开展风险分析；何巍等(2014)结合具体矿山实例总结了蒙特卡洛模拟方法的优缺点和常见随机数产生方式。

总结国内外矿产资源可供性及风险分析方面的相关文献，主要表现出两方面特征：一是现有研究成果主要在确定价格的基础上开展资源可供性研究；二是主要针对特定矿山企业生产经营活动开展风险分析。本文以河北省区域内的铁矿资源为研究对象，在分别对13个典型矿区确定价格条件下资源可供性分析评价及不同价格情景状态下区域铁矿资源总体可供性分析测评的基础上，从市场需求动态变化的角度出发，对矿产品价格不确定背景下区域铁矿资源可供性风险进行研究，以便对矿产品市场急剧变动状况下区域矿产资源可供性及其风险程度做出客观的分析和判断。

要知道，减重并不等于减脂，目前市场上的很多智能体脂称是可以针对人体的核心数据做出深度解析，科学指导，让你能够更快、更准确地掌握自己的身体健康状况，不再盲目运动，该怎么减，一目了然，让减脂塑形事半功倍。还可以定制你的私人运动计划，根据你的运动习惯及身体近况推荐可达成的目标。此外，通过手机连接上与之匹配的APP，可在APP中看到你的身体状况并综合分析提供定制级膳食运动计划，健身瘦身其实可以很简单。

1　确定价格状态下单矿区资源可供性评价

1.1　评价方法与模型

对于单个矿区而言，探明资源储量能否成功转化为现实可供的矿产品，取决于矿山生产经营企业能否获得满意的收益水平。如果实际或未来预期收益水平能够达到目标收益要求，矿山企业就会持续开展生产和经营，矿区资源的潜在可供性就会成功转化为现实可供性；否则，当收益水平达不到目标收益要求，矿山企业就会放弃矿山投资项目或终止生产经营活动，矿区资源储量就不能转化为现实矿产品可供应量。对于具体矿区而言，在已探明资源储量规模和品质、资源开采和加工技术等因素确定的条件下(即矿山企业开采加工成本确定的条件下)，矿山企业能否获得预期收益，决定于矿产品市场的价格水平，如果矿山企业的可供价格低于市场价格，则矿区的储量、产能和产量都是可供的；如果矿山企业的可供价格高于市场价格，矿区的储量、产能和产量都是不可供的。因而，单个矿区资源的可供性评价就转化为特定价格状态下对矿山企业动态经济效益的分析与测评。具体评价模型为：

/rern)

式中：NPV为矿山整个寿命期内预期获得的财务净现值(寿命期总收益水平)，cn为第n年的净现金流量(第n年现金流入与现金流出的差额)，n为折现年序号，l为矿山服务年限，r为折现率，e为常数2.72。

1.2　基础数据来源及可供性评价

蒙特卡洛模拟技术是利用服从某种分布的随机数对现实系统中可能出现的随机现象进行模拟，该方法假设独立的随机变量为xii=1,2,…,k，其对应的概率密度函数为fxq，fx2，…，fxk，功能函数为y=gx1,x2,…,xk。根据随机参数xi的分布，利用抽样技术产生随机数x1，x2，…，xk，将其代入功能函数，经过N次模拟计算，可以产生 N个功能函数yj=gx1,x2,…,xk的随机量，其中j=1,2,…,N。如果这 N个随机量中有 M个的功能函数值大于等于项目预期目标或规定的值 x0，当N充分大时，根据概率论大数定理或中心极限定理，可以测得项目的风险率 p=limpgx1,x2,…,xk≥x0=M/N，依据此模拟计算结果，便可得到地区矿产资源可供性概率分布，从而对地区资源可供性风险程度作出具体的分析和判断。

表1　矿区基础数据汇总表

Table 1　Basic data of mining areas

矿区编号开发状态矿山投资(万元)建设期(年)生产期(年)生产规模(万吨/年)单位成本(元/吨)销售价格(元/吨)一号矿区正在开采2089．47471．34863．92二号矿区正在开采2055．92464．52863．92三号矿区正在开采17493．96512．21863．92四号矿区正在开采1914．33520863．92五号矿区规划36000．0023340．00500863．92六号矿区规划86654．223．51835．91400863．92七号矿区规划2754．9131315．00520863．92八号矿区规划5268．1731010．00520863．92九号矿区规划79147．95515100．00520863．92十号矿区规划24564．744916．35400863．92十一号矿区规划25790．5531540．00520863．92十二号矿区规划12172．723733．33570863．92十三号矿区(露天开采部分)规划十三号矿区(地下开采部分)规划2307．2921416．67520863．92261．67520863．92

利用上述分析评价模型，通过试错迭代计算，可以计算出各矿区矿产品可供价格(表3)，依据各矿区资源可供价格系列分布，测算在不同价格背景下区域矿产资源可供应量分级状态，并绘制出区域铁矿资源可供性状况分布图(图1)。

表2　单矿区可供性评价结果汇总表

Table 2　Availability evaluation of mining areas

矿区编号财务净现值内部收益率所得税前(万元)所得税后(万元)所得税前(%)所得税后(%)最敏感因素盈亏平衡点产能可供程度产能利用率(%)平衡点产量(万吨)可供性状态可供量(万吨/年)可供性排序一号矿区241803．71193361．70--价格40．2936．05可供89．472二号矿区117434．0193625．93--价格36．2220．25可供55．953三号矿区1276631．301123371．09--价格59．15292．20可供493．961四号矿区20764．4416178．92--价格42．246．05可供14．339五号矿区42185．3430088．7329．7225．19价格39．1415．66可供40．004六号矿区-20．24-9539．4113．0010．91价格36．6513．16不可供0．00-七号矿区2166．161479．6827．3023．22价格44．412．22可供15．0012八号矿区5538．923806．1133．9527．97价格48．384．84可供10．0010九号矿区35547．6520495．0120．9117．87价格26．3526．35可供100．005十号矿区4534．84157．6116．6213．13价格24．403．99可供16．3511十一号矿区21539．7113730．2626．3521．88价格33．2713．31可供40．008十二号矿区24125．0418737．84158．0294．12价格40．6013．53可供33．336十三号矿区(露天开采部分)23848．4919716．60166．82145．89价格41．356．89可供16．677十三号矿区(地下开采部分)1614．741237．3241．0434．84价格33．270．56可供1．6713

从以上评价结果可以看出，除第六号矿区以外，其它矿区的财务净现值均大于零、财务内部收益率均大于13%(达到了国家规定的铁矿资源开发项目最低收益水平和标准要求)；各个矿区的盈亏平衡点较低，较容易达到盈亏平衡产量，说明各矿区资源开发利用投资具有较强的抗风险能力。评价结果表明：①在确定价格状态下(863.92元/吨)，河北省铁矿资源的整体可供性程度较高，最大可供量达到926.73万吨/年；②地区铁矿资源可供性程度重点集中在大中型矿区，排名前五位矿区的可供量达到779.35万吨/年，占区域年度总体可供量的84.10%，说明大中型矿区的成本优势，对有效提升矿区铁矿资源可供性具有明显作用；③矿产品销售价格成为影响各矿区可供性程度高低的重要风险因素，当铁精粉市场价格下降时，将会导致地区总体资源可供性程度大幅度降低。

M, s |= CGφ，当且仅当，对所有和s具有可及关系的t通过有穷序列的任意的～i(i∈G)步: M, t |= φ。

2　不同价格情景下的区域矿产资源可供性测评

2.1　MINSIN系统模型构建

MINSIN系统模型是美国矿业局研发的一种用于测算矿区资源可供性分析方法，本系统模型是在输入矿山建设投资、经营成本、收益指标(预期内部收益率)等参数基础上，通过系统模型分析计算，确定矿区矿产品可供价格；利用区域内不同矿区矿产品可供价格，划分地区矿产资源不同价格背景下的可供性程度和水平。具体评价模型如下：

×Mn-COn)(er-1/rern)

式中：NPV为矿山整个寿命期内预期获得的财务净现值(寿命期总收益水平)；P为矿产品价格；Mn为第n年矿产品产量；COn为第n年总成本；n为折现年序号；l为矿山服务年限；r为设定内部收益率值；e为常数2.72。

令NPV=0，求取矿产品可供价格P，再按照某个时间节点的矿产品市场价格Pt，来分析判断区域铁矿资源可供性及可供量。

由于CATS跳变到2018年4月24日，同时网关检测到其与时钟源CATS主机相差1 000 s，立即停止了网关的NTP服务；网关没有利用主备机动态控制NTP时钟服务机制，但由于前期NTP协议工具Meinberg不太稳定，出现过自动停止的情况，因此实施过程对NTP服务创建了计划任务自启动功能；由于网关在06:00设置了自启动NTP任务，因此在06:00时网关同样跳变到2018年4月24日。

2.2　区域资源可供量测算

利用上述分析评价模型，编制矿山企业生产经营财务现金流量表，价格水平选取2014年铁精粉市场平均价格863.92元/吨。分别计算14个评价对象的财务净现值和财务内部收益率，并对每个矿区预期经营状况进行盈亏平衡分析和敏感性分析，得出13个矿区矿产资源可供性分析评价结果(表2)。

表3　各矿区矿产品可供价格及可供量汇总表

Table 3　Available prices and quantity of mining areas

矿区编号可供价格(元/吨)可供量(万吨/年)矿区编号可供价格(元/吨)可供量(万吨/年)一号矿区矿410．6689．47八号矿区690．2110．00二号矿区511．7255．92九号矿区730．60100．00三号矿区557．94493．96十号矿区766．8816．35四号矿区617．8614．33十一号矿区714．0540．00五号矿区685．6840．00十二号矿区686．5333．33六号矿区864．1235．91十三号矿区(露天开采部分)537．0116．67七号矿区693．5415．00十三号矿区(地下开采部分)632．821．67

图1　铁矿石资源可供性分布图 Fig.1　Availability distribution of iron ore resources

上图中曲线为财务净现值(所得税前)等于零时区域铁矿资源的可供性曲线。由图中可供性曲线分布状态可以看出：①当铁精粉市场价格为500元/吨时，从500元处引出一条与X轴平行的直线，该平行线与可供性曲线交于一点，若矿区位于该交点左侧的曲线上，说明该矿区的产量是可供的，若矿区位于该交点右侧的曲线上，则该矿区的产量不可供，若矿区产量位于该交点，说明矿区的产量是该价格条件下的最大可供量。由此可得，当铁精矿产品处于500元/吨价格背景下时，本地区铁矿资源最大可供量为89.47万吨/年；②当铁精矿市场价格为557.94元/吨时，本地区该价格背景下铁矿资源最大可供量为655.75万吨/年；③当铁精矿市场价格为730.60元/吨时，本地区该价格背景下铁矿资源最大可供量为910.08万吨/年；④当铁精矿市场价格为863.92元/吨时，本地区该价格背景下铁矿资源最大可供量为926.73万吨/年。

3　区域矿产资源可供性风险分析

3.1　市场价格波动状态下的资源可供性风险分析方法

其中：/rern)

河北省铁矿资源较为丰富，主要集中分布在唐山、张家口以及邯邢地区；根据当前本地区铁矿资源开发利用状况及开发利用规划方案，本次单矿区资源可供性评价包括13个矿区、14个评价对象(第十三号矿区包括露天开采和地下开采两个规划部分)，基础数据来源：正在开采矿区，取矿山企业实际生产经营数据；规划开采矿区，取矿山建设项目可行性研究报告中的规划预测数据(表1)。同时，由于各个矿区生产的铁矿石原矿品位各不相同，有些甚至相差较大，因此本次资源可供性评价统一选取品位为63%的铁精矿产品为评价对象，各矿区所计算的铁矿石产量均为铁精矿可供量。

据统计，“十二五”期间，我国由公路建设引起的地表扰动破坏总面积为119.6万公顷，占开发建设项目干扰总面积的21.6%。其中，山地、丘陵地区75.3万公顷，占公路干扰破坏总面积的62.97%。且山地丘陵区自然条件较差，一旦被破坏，极难被恢复。

3.2　河北省铁矿资源矿产品价格分布特征

通过河北省主要铁矿资源价格信息报告以及国内外铁矿石价格行情分析，可得到河北省2006年7月至2015年6月铁精矿平均价格发展变化分布特征(图2)。由图中可以看出，2006年7月以来，河北省铁精矿价格经历了一个跌宕起伏的变化；在2008年6至7月达到最高值1270.21元/吨，最低价格跌至489.84元/吨；2011年至今又出现缓慢的下降趋势，在2015年6月已下降至499.28元/吨。

图2　2006～2015年铁精矿平均价格曲线图 Fig.2　Average prices of iron ore in 2006～2015

3.3　各矿区铁矿资源可供性风险分析

本次模拟分析以矿区生产经营财务净现值(FNPV)为功能函数，以地区2006年7月～2015年6月每月铁精矿平均价格变化序列为随机变量，利用蒙特卡洛模拟方法crystal ball软件包进行动态模拟计算。先后对13个矿区、14个评价对象分别运行1000、5000、10000次模拟，将输出结果进行对比分析发现，当模拟次数达到5000 时，13个矿区财务净现值的均值和标准差均已达到收敛状态，由此可以得到各矿区资源可供性模拟分析结果及可供性状态概率分布(图3～16、表4)。

图3　一号矿区分析结果 Fig.3　Result of No.1 mining area

图4　二号矿区分析结果 Fig.4　Result of No.2 mining area

图5　三号矿区分析结果 Fig.5　Result of No.3 mining area

图6　四号矿区分析结果 Fig.6　Result of No.4 mining area

图7　五号矿区分析结果 Fig.7　Result of No.5 mining area

图8　六号矿区分析结果 Fig.8　SResult of No.6 mining area

图9　七号矿区分析结果 Fig.9　Result of No.7 mining area

图10　八号矿区分析结果 Fig.10　Result of No.8 mining area

图11　九号矿区分析结果 Fig.11　Result of No.9 mining area

图12　十三号矿区(露天)分析结果 Fig.12　Result of No.13 mining area(open)

图13　十三号矿区(地下)分析结果 Fig.13　Result of No.13 mining area(underground)

图14　十号矿区分析结果 Fig.14　Result of No.10 mining area

图15　十一号矿区分析结果 Fig.15　Result of No.11 mining area

图16　十二号矿区分析结果 Fig.16　Result of No.12 mining area

3.4　区域总体资源可供性风险分析

(ⅱ) 假设(u,v)是系统(3)的一个共存解,则a>λ1,因此θa存在且唯一,从而u<θa,又因d>λ1,则v>θd,故

×ki

由以上资源可供性分析测评结果可以看出，区域矿产资源可供性程度与矿产品市场价格具有紧密的相关关系，而矿产品市场由于受到政治经济形势、供需状况、资源保有储量、开发技术水平、矿产品贸易政策等众多因素的影响，矿产品市场价格经常会出现较大幅度的变化，这就要求必须开展价格波动状态下区域矿产资源可供性分析，以便对地区一定时期内矿产资源可供性风险程度做出科学全面的分析和判断。

3) 遥控启动。当主机处于备车完毕状态、控制部位在集控台时，将集控台车钟手柄从停车位置扳向正车(或倒车)任意位置，系统向电喷控制系统发出正车(或倒车)启动指令，由电喷控制系统控制主机正转(反转)启动。

式中：Q为区域矿产资源总体可供量；qi为第i个矿区资源可供量；NPVi为第i个矿区整个寿命期内某一随机价格状态下的财务净现值(寿命期总收益水平)；ki为矿区不可分割性参数，当NPVi≥0时，令ki=1，当NPVi<0时，ki=0；i为矿区分析对象编号(i=1，2，3…，14)。

表4　各矿区矿产品可供概率及可供量汇总表

Table 4　Availability probability and available quantity of mining areas

矿区编号可供概率(%)可供量(万吨/年)矿区编号可供概率(%)可供量(万吨/年)一号矿区100．0089．47八号矿区63．5010．00二号矿区89．7555．92九号矿区63．86100．00三号矿区82．01493．96十号矿区60．4316．35四号矿区71．4314．33十一号矿区57．5240．00五号矿区67．2840．00十二号矿区61．6833．33六号矿区42．7435．91十三号矿区(露天开采部分)86．0516．67七号矿区63．1115．00十三号矿区(地下开采部分)69．791．67

通过crystal ball软件对各矿区进行蒙特卡洛模拟分析运算，得出每个矿区的可供概率，与此同时，抽取的5000个随机价格样本数据，计算得出相应的财务净现值，再利用以上模型的综合模拟计算，便可得到本地区铁矿资源总体可供性风险分析结果(表5)。

在以上各单个矿区资源可供性风险分析的基础上，可以构建以下区域总体资源可供性风险分析模型，来分析判断地区资源可供性总体风险程度。

从以上分析结果可以看出，河北省铁矿资源可供量处于900.35～952.61万吨/年的概率为60.44%；供应量为800.35万吨/年的概率为1.26%；可供量处于712.02～760.35万吨/年的概率为5.60%；可供量处于656.02～672.02万吨/年的概率为14.72%；可供量处于145.39～162.06万吨/年的概率为7.74%；可供量处于89.47万吨/年的概率为10.24%。表明在当前铁矿石市场价格发展变化状态下，河北省铁矿资源总体可供性状态处于900～950万吨/年的水平，可供性保证程度只有60%，说明地区铁矿资源可供性存在着较大的风险性。

表5　区域铁矿资源可供性状态概率分布汇总表

Table 5　Availability probability distribution of regional iron ore resources

区域资源可供量(万吨/年)可供频率可供概率(%)区域资源可供量(万吨/年)可供频率可供概率(%)89．4751210．24745．35180．36145．391853．70755．35190．38162．062024．04760．35711．42656．0252910．58800．35631．26670．35821．64900．351462．92672．021252．50916．7073814．76712．021723．44952．61213842．76

4　研究结论及对策建议

本文依据河北省铁矿资源13个矿区的资源特征及开发利用条件，开展了确定价格状态和不同价格情景下地区铁矿资源可供性评价，并结合一定时期内铁矿石市场价格波动变化形成的价格序列分布，对地区铁矿资源可供性进行了风险分析。研究结果表明：河北省铁矿资源具备一定的可供性，可供规模为900～950万吨/年，但可供性概率比较低，可供性保证程度只有60%左右，标志着地区铁矿资源可供性存在着较大的风险性。结合地区铁矿资源地质特征、开发利用条件及矿产品市场的发展变化，提出以下对策建议：

①九麦2号、②中麦895、③小偃22（CK）、④秦农578、⑤西农223、⑥陕农33、⑦武农6号、⑧凳峰168。

(1)建立铁矿资源可供性预警系统，及时跟踪和预测地区资源可供性发展变化态势。矿产资源可供性评价是一个包含资源特征、开发利用条件、开采加工技术工艺、投资规模、成本水平、税费标准、矿产品价格等众多影响因素和参数的综合评价过程；一次可供性评价往往只能对特定时间节点和确定评价参数状态下地区资源可供性作出分析和判断，随着时间的变化，各项评价参数、特别是矿产品市场价格的波动性变化，将直接影响和改变资源可供性状态。因此，必须建立地区资源可供性评价预警系统，根据资源开发利用条件、技术工艺及矿产品市场的供需条件变化，适时开展资源可供性分析研究，对资源供需形势及时作出科学的分析、预警和判断，为地区矿产资源战略规划的制定、矿产资源政策的实施与调整提供决策参考。

(2)强化矿产资源开发利用及生产经营过程技术进步，不断降低矿产品生产加工成本。针对本地区铁矿资源贫矿多、富矿少、小型矿区多、大中型矿区少、开采深度大、选矿难度大等自然特征和地区铁矿资源开采加工技术落后、开发成本偏高的现实状况。努力通过资源开发利用新技术工艺的研究与开发，实现资源开采加工技术的新突破；积极引进国内外先进的生产加工技术和工艺，不断改进矿山企业的生产技术水平；有效利用高新技术产业的成熟技术成果，加快对资源开发利用产业的技术改造和升级；强化对资源开发利用产业的技术整合，利用少量增量资产投入带动存量资源技术水平大幅度提升等方法、途径和措施，大力改进地区资源开发利用的技术水平，不断降低矿产品生产加工成本，有效提升地区铁矿资源可供性的现实能力和水平。

(3)加强区域铁矿资源勘查，不断增加地区资源可供基地和可供程度。河北省位于华北地台北缘，出露地层齐全，构造和岩浆活动频繁，具备良好的铁矿资源成矿地质条件，形成一系列不同类型及具有地域特点的铁矿资源成矿带和成矿区域，铁矿资源异常区和矿化点分布广泛，潜在资源量可观，仅沉积变质型铁矿的潜在资源量保守预测有21亿吨。这些潜在资源量通过深入的地质勘查工作，转化为现实可开发的资源量，将会大大提高地区资源的可供能力。

(4)广泛构筑外部资源供应体系，提高地区铁矿资源的保障程度。河北省地处东北亚经济圈的有利地带，与俄罗斯、朝鲜、蒙古等这些铁矿资源丰富的国家相毗邻；同时，优越的港口条件，为有效利用澳大利亚、印度、巴西、南非等国家优质、丰富的铁矿资源提供了有利条件。可以通过获取境外资源勘查权和采矿权，从源头上保障资源的可供性；通过重点地区选择与宏观布局相结合，实现资源供给的全球优化配置；通过以资源开发企业为主体的勘查、开发、加工和贸易一体化模式，有效降低地区资源可供性风险等方式与途径，广泛开拓外部资源市场，构筑完善的资源供应和保障体系。

[References]

Bastante F G,Taboada J, Alejano L, Alonso E. 2008. Optimization tools and simulation methods for designing and evaluating a mining operation[J].Stochastic Environmental Research & Risk Assessment, 22(6): 727-735

He wei,Huan Bing-han,Wu Xiao-ming.2014.Application of Monte Carlo method in the risk analysis of mining investment[J].Value Engineering,(2):324-325(in Chinese with English abstract)

Li Hai-yang.2012.Demand prediction and availability analysis of mineral resources in Liaoyang[J].Science & Technology Information,(18):243(in Chinese)

Li Yu-wei.2015.The principle and methodology of mineral availability analysis[J].Nature Resource Economics of China,28(2):8-13(in Chinese with English abstract)

Maekay R M, Probert S D. 2001. Forecasting the United Kingdom’s supplies and demands for fluid fossil-fuels[J].Applied Energy, 69: 161-189

Mennenga M,Thiede S, Beier J,Dettmer T,Kara S,Herrmann C. 2012. A forecasting model for the evaluation of future resource availability [M]. Leveraging Technology for a Sustainable World: 449-454

Sohn I. 2006. Long-term projections of non-fuel minerals: We were wrong, but why?[J].Resources Policy, 30: 259-284

Zeng Lan,Song Jian-jun.1995.Analysis of the availability of regional resources in China in 1996-2010[J].Geology in China,22(9):14-16(in Chinese)

Zhang Shao-hua,Li Chang-you.2003.Analysis of the availability of some non-metallic mineral resources[J].China Non-metallic Minerals Industry,(2):7-11(in Chinese)

[附中文参考文献]

何　巍,宦秉炼,吴晓明.2014.蒙特卡洛模拟方法在矿业投资风险分析中的应用[J].价值工程,(2):324-325

李海洋.2012.辽阳市矿产资源需求预测及可供性分析[J].科技资讯,(18):243

李裕伟.2015.矿产资源可供性分析的原理与方法[J].中国国土资源经济,28(2):8-13

曾　澜,宋建军.1995.1996-2010年我国区域资源可供性分析[J].中国地质,22(9):14-16

章少华,李常有.2003.若干非金属矿资源可供性分析[J].中国非金属矿工业导刊,(2):7-11