Study on Mineral Chemistry of Liziping-Fulongchang Pb-Zn-Ag Deposit in Western Yunnan and Its Metallogenic Significance
-
摘要:
通过光学显微镜、粉晶X射线衍射分析仪(XRD)、场发射扫描电子显微镜(FESEM)、电子探针(EPMA)等方法查明了李子坪和富隆厂铅锌银矿床的矿物组成和矿物化学成分。研究表明：李子坪矿床的硫化物/硫盐矿物主要有闪锌矿、细硫砷铅锌、方铅矿、含砷黄铁矿；富隆厂矿床的硫化物/硫盐矿物主要有闪锌矿、方铅矿、灰硫砷铅锌、黄铁矿。矿物化学特征指示两矿床经历了含Pb-Zn成矿流体对富As黄铁矿角砾的交代，As元素的活化；闪锌矿、方铅矿沉淀；富As、(Sb、Cu)残余流体对方铅矿的交代形成硫盐的过程。各矿段硫盐类型的差异由成矿流体组分和物理化学条件的共同控制。针对富含As、Sb元素的特殊型沉积岩容矿铅锌银成矿体系开展的矿物学和矿物化学研究，有助于促进类似矿床成矿作用的理论认识。
-
关键词:
- 兰坪盆地 /
- 沉积岩容矿铅锌矿床 /
- 矿物化学 /
- As-Sb-Pb硫盐
Abstract:The mineral composition and mineral chemical composition of Liziping and Fulongchang lead-zinc-silver deposits have been ascertained in detail by means of optical microscopy, XRD, FESEM and EMPA. The mineralo-gical study shows that the sulfide/sulfosalt minerals in Liziping deposit mainly include sphalerite, gratonite, galena, arsenic-bearing pyrite. The sulfide/sulfosalt minerals in the Fulongchang deposit mainly include sphalerite, galena, jordanite and pyrite. The mineral chemical characteristics indicate that the two deposits have undergone metasomatism and activation of As elements by Pb-Zn ore-forming fluids to As-rich pyrite breccia; precipitation of sphale-rite and galena; the process of formation of sulfosalts by replacement of galena by As-(Sb, Cu) rich residual fluids. The difference of sulfur salt type in each ore section is controlled by ore-forming fluid composition and physicochemical conditions. In this paper, the mineralogy and mineral chemistry of the ore-bearing Pb-Zn-Ag metalloge-nic system of sediment-hosted rich in As and Sb elements will help to promote the theoretical understanding of mineralization of similar deposits.
-
传统的投资组合绩效评价均采用方差度量风险[1]，如夏普指数、特雷诺指数和詹森指数.但随着方差的广泛运用，学者们发现方差并不是一个良好的测量风险的方式，方差衡量的是随机变量与期望值之间的上下偏差，而投资者更关心下偏差[2].而风险价值(VaR)作为一种下偏差风险度量方法，在风险测量和绩效评价中得到了广泛的应用，如：VaR风险度量方法[3]；在均值-方差模型基础上的均值-VaR模型[4-5]；不允许卖空情况下的均值-方差模型与均值-VaR模型[6]；具有基数约束、上下界限制和交易成本的多阶段均值-VaR模型[7].
在投资组合模型的建构方面，一些学者考虑实际投资市场的诸多摩擦因素，构建了不同的模型并求解.如：研究了不同借贷利率下以VaR为风险度量方法的投资组合效用最大化模型[8]；在其均值-方差投资组合选择中考虑借款约束[9].交易成本同样对投资组合影响颇大，BEST和HLOUSKOVA[10]提出了分段线性交易成本的投资组合模型.另外，学者们研究了具有基数约束(限制最优资产组合中所含资产的数量)的投资组合模型，使该模型更加符合实际市场上的投资情况.如：在其均值-方差动态投资组合和市场时机选择中考虑了基数约束这一要素[11]；在研究一系列具有基数约束的投资组合选择问题时提出了一个新的混合整数二次约束的二次规划算法[12]；研究了含区间系数的基数约束投资组合模型，并运用一种新的动态差分进化算法与自适应控制参数求解转换后的确定性模型[13].
实际上，具有基数约束的投资组合模型是一个复杂的非光滑的混合整数规划问题，计算十分复杂，模型真实前沿面的解析解难以得到.因此，许多学者运用数据包络分析(Data Envelopment Analysis，DEA)法[14].如：运用传统的DEA效率评价模型优化投资组合选择[15]；在基金市场上，对30只对冲基金在规模报酬不变的情况下进行了DEA效率估计排名[16]；提出了多元一致性数据包络分析效率分析模型，在传统DEA模型的基础上考虑了投资组合的多样化[17]；将改进的DEA交叉效率模型应用于韩国股票市场投资组合选择，得到了优于传统组合选择方法的结论[18]；构建了多阶段的DEA效率估计模型，在股票市场上对10只风险证券进行了效率估计[19]；比较了前沿面方法和传统评价体系比率在投资决策设计应用中的效果，得到了前沿面评价方法更注重持续性的结论[20]；提出了一种多目标投资组合选择模型，并使用DEA模型进行了交叉效率评价[21].
本文使用BCC-DEA绩效评价方法，对具有交易成本、上下界约束、借款约束和基数约束的均值-VaR模型进行绩效评价，通过上海证券市场的真实数据进行效率估计，构建不同样本数据下的DEA效率估计模型，并构建前沿面与真实前沿面进行对比，分析DEA模型对投资组合绩效评价的有效性.
1. 均值-VaR投资组合优化模型
1.1 符号描述和定义说明
假设可供选择的风险资产有n种，第i种风险资产的收益率(随机变量)表示为Ri，其均值ri=E(Ri), 协方差矩阵G=(σij)n×n，σij=Cov(Ri, Rj) (i, j=1, 2, …, n)；xi (i=1, 2, …, n)表示第i种风险资产投资比例；xn+1=1-(x1+x2+…+xn)表示无风险资产投资比例；xn+10表示无风险资产投资比例下界限制，xn+1≥xn+10；记投资组合x=(x1, x2, …, xn, xn+1)T；rb、rl分别表示无风险资产的借款利率、贷款利率；rp、rN分别表示投资组合x的期望收益率、净收益率；Ci表示第i种风险资产的单位交易成本；r0表示投资者要求的最低期望收益率；ui、li分别表示xi的上界、下界；K为投资组合中风险资产投资比例为非负的个数.投资组合的期望收益率为rp=rTx，方差为σp2=xTGx，标准差为σp=√xTGx, 其中，r=(r1, r2, …, rn)′, ri表示第i种资产的期望收益率.
定义1 风险(VaR)是指在一定的置信度下，某一投资组合在未来特定时间内最大的可能损失.设投资组合的期望收益率为rp, 则
P(rp<−VaR)≤1−c (1) 为投资组合的VaR(VaR取正数)，其中，c(1/2≤c≤1)为常数.
式(1)表示投资组合的收益率超过-VaR的概率不低于c.假设其收益率服从正态分布N(rp, σp2)，根据中心极限定理可以得到如下定理:
定理1[6] 当投资组合中n种资产的收益率服从正态分布时，式(1)可转换为
rp≥Φ−1(c)σp−VaR, (2) 其中，Φ(·)是标准正态分布函数，Φ-1(c)是置信度为c的正态分布函数的下分位点.
随机市场情况下, 由式(2)可得VaR风险表达式为Φ-1(c)σp-r Tx.
1.2 均值-VaR投资组合模型
假设第i种资产的交易成本函数C(xi)是投资比例的分段线性函数，且为凸函数(图 1)，即
C(xi)={k1⋅xi(0≤xi≤a),k2⋅xi−(k2−k1)⋅a(a<xi≤1,k2>k1>0). (3) 借款约束是投资组合考虑的另一个重要因素.本文考虑了无风险资产借款，则投资组合x=(x1, x2, …, xn, xn+1)T的期望收益率可表达为
rp=n∑i=1rixi+rn+1(1−n∑i=1xi), (4) 其中，
rn+1={r1(1−n∑i=1xi≥0),rb(1−n∑i=1xi<0), rb≥rl.若1−n∑i=1xi⩾0, 表示存在无风险资产贷款；若1−n∑i=1xi⩽0, 表示存在无风险资产借款.
此时投资组合x=(x1, x2, …, xn, xn+1)T的净收益率为
rN=n∑i=1rixi+rn+1(1−n∑i=1xi)−C(xi). (5) 假设K≥0，且为整数，zi={0, 1}，则投资组合x=(x1, x2, …, xn, xn+1)T的基数约束为
lizi≤xi≤uizi(i=1,2,⋯,n),z1+z2+⋯+zn≤K. (6) 投资组合模型的上下限约束为：li≤xi≤ui (i=1, 2, …, n).
在文献[22]的基础上，考虑交易成本、上下界约束、借款约束和基数约束等条件，提出以下3个模型.
(1) 风险最小化的均值-VaR投资组合模型为：
minΦ−1(c)σp−rTxs.t{n∑i=1rixi+rn+1(1−n∑i=1xi)−n∑i=1C(xi)≥r0,1−n∑i=1xi≥x0n+1,lizi≤xi≤uizi(zi=0,1,i=1,2,⋯,n),n∑i=1zi≤k, (7) 其中，r0表示投资者预期给定收益率.模型共有4个约束条件：投资组合的净收益率不低于给定的期望值r0、无风险资产有借贷约束、第i种风险资产的投资比例有上下界约束和风险资产投资比例非负个数不能超过K.模型的经济意义是指在满足上述4个约束条件下，如何使投资组合风险(VaR)最小.
(2) 收益最大化的均值-VaR投资组合模型为：
maxn∑i=1rixi+rn+1(1−n∑i=1xi)−n∑i=1C(xi)s.t.{Φ−1(c)σp−rTx≤VAR0,1−n∑i=1xi≥x0n+1,lizi≤xi≤uizi(zi=0,1,i=1,2,⋯,n),n∑i=1zi≤k, (8) 其中，VaR0表示投资者预期给定的风险.模型共有4个约束条件：投资组合的风险不高于预期风险VaR0、无风险资产有借贷约束、第i种风险资产的投资比例有上下界约束和风险资产投资比例非负个数不能超过K.模型的经济意义是指在满足上述4个约束条件下，如何使投资组合收益(r)最大.
(3) 效用最大化的均值-VaR投资组合模型为：
max(1−θ)[n∑i=1rixi+rn+1(1−n∑i=1xi)−n∑i=1C(xi)]−θ[Φ−1(c)σp−rTx]s.t.{1−n∑i=1xi≥x0n+1,lizi≤xi≤uizi(zi=0,1,i=1,2,⋯,n),n∑i=1zi≤k, (9) 其中，θ (0≤θ≤1)表示投资者风险偏好系数.当θ=1时，投资者完全厌恶风险.当θ=0时，投资者只关心收益，不管风险，此时，模型(9)可以转变为：
maxn∑i=1rixi+rn+1(1−n∑i=1xi)−n∑i=1C(xi)s.t.{1−n∑i=1xi≥x0n+1,lizi≤xi≤uizi(zi=0,1,i=1,2,⋯,n),n∑i=1zi≤k, (10) 模型(10)属于非光滑的混合整数规划问题，可以运用遗传算法[23]得到其最优解，从而得到该模型目标函数最大值(即r0max)，也可得到VaR0max.
2. 均值-VaR投资组合的DEA评价方法
2.1 投资组合效率的定义
为更好地评估投资组合的绩效表现，文献[24]提出了真实前沿面的概念，并基于真实前沿面定义了投资组合效率.本文根据模型(10)对应的真实前沿面来定义n只风险资产构成的投资组合效率(图 2). 图中曲线MNW是n只风险资产构成的投资组合真实前沿面，Q是某一个投资组合.
收益导向的投资组合效率(PEr)是指相同风险下投资组合实际期望收益率与理想期望收益率的比值；风险导向的投资组合效率(PEv)是指相同期望收益率下理想风险与投资组合实际风险的比值.投资组合Q的收益导向、风险导向的效率可以分别表示为：
PEr=rqrw,PEv=VmVq. 2.2 均值-VaR的DEA效率评价模型
由于金融市场上的诸多摩擦因素，大规模投资组合难以计算，所以，当投资组合真实前沿面为凹函数时，本文利用DEA模型对决策单元进行投资组合绩效评价.
LIU等[25]运用BCC-DEA模型进行投资组合绩效评价，并指出：若r=h(VaR)是投资组合的前沿面，则r′=h′m(VaR)是m个样本的投资组合的BCC-DEA模型的前沿面，当m→+∞时，h′m(VaR)依概率收敛于h(VaR).另外，若投资组合的风险函数VaR(x)(x为投资组合比例)是凸函数，收益函数r(x)是凹函数且可行域Ω是凸集，则由模型(7)得到的投资组合的真实前沿面是凹函数，可以采用基于数据的BCC-DEA模型对决策单元进行绩效评价.随着决策单元数量的增加，其前沿面能较好地逼近真实的前沿面.
假设有m个投资组合可供评价，其中第j个投资组合xj=(x1j, x2j, …, xnj, xn+1j)′，j=1, 2, …, m，初始投资组合x0=(x10, x20, …, xn0, xf0)′，第j个投资组合的收益率为
rj=n∑i=1rixji+rn+1(1−n∑i=1xji)−n∑i=1C(xji), (11) 第j个投资组合的VaR值为
VaR(xj)=Φ−1(c)σp−[rixji+rn+1(1−xji)], (12) r0为预先给定的可能性收益率、VaR0为可能性VaR值.可使用以下3个BCC-DEA评价投资组合的有效性.
(1) 风险导向下均值-VaR投资组合的BCC-DEA评价模型：
minθs.t.{m∑j=1λjrj≥r0,m∑j=1λjVaRj≤θVaR0,m∑j=1λj=1,λj≥0(j=1,2,⋯,m), (13) 其中，λj表示第j种投资组合的比例，θ表示投资组合的效率值.
(2) 收益导向下均值-VaR投资组合的BCC-DEA评价模型：
maxθs.t.{m∑j=1λjrj≥θr0,m∑j=1λjVaRj≤VaR0,m∑j=1λj=1,λj≥0(j=1,2,⋯,m). (14) (3) 效用最大化的均值-VaR投资组合的BCC-DEA评价模型：
min(1−s+/VaR0)(1+s−/r0)s.t.{m∑j=1λjrj−s−≥r0,m∑j=1λjVaRj+s+≤VaR0,m∑j=1λj=1,λj≥0(j=1,2,⋯,m), (15) 其中，s+代表正偏差，s-代表负偏差.
以上3个模型都是线性规划模型，可以运用线性规划的旋转算法[22]求解.
3. 实证研究
本文选择上海A股市场表现较好的30只股票进行分析，股票代码分别为：S1(600000)、S2(600010)、S3(600016)、S4(600029)、S5(600030)、S6(600036)、S7(600048)、S8(600104)、S9(600109)、S10(600111)、S11(600518)、S12(600519)、S13(600637)、S14(600887)、S15(600893)、S16(601166)、S17(601186)、S18(601318)、S19(601336)、S20(601377)、S21(601390)、S22(601601)、S23(601668)、S24(601669)、S25(601688)、S26(601766)、S27(601800)、S28(601818)、S29(601985)、S30(601989).选择从2010年1月至2019年5月的周收益数据作为样本数据.
鉴于股票市场的高波动性，金融市场上往往通过构造股票投资组合来分散投资风险，而其构成投资组合的股票之间的相关系数越低，分散效果越好.由30只股票中每一只的股票收益及其与其他29只股票的相关系数(表 1)可知：样本中仅3只股票(S5、S21、S25)的相关性系数大于0.5，呈现中度相关性；20只股票的相关性系数介于0.3~0.5之间，具有较弱的相关性；其余7只股票的相关性系数介于0~0.3之间，具有极弱的相关性，可视为不相关.由此可知，本文所选取的股票样本的总体相关性较弱，可以较好地达到分散风险的效果.
表 1 股票收益率均值及其相关性Table 1. The expected return of stock and its correlation股票代码 S1 S2 S3 S4 S5 S6 S7 S8 S9 S10 收益率均值 0.002 53 0.002 70 0.002 93 0.002 64 0.001 71 0.001 64 0.002 71 0.002 32 0.002 54 0.004 37 相关系数 0.259 42 0.418 26 0.259 31 0.424 06 0.507 93 0.250 49 0.431 61 0.324 97 0.489 24 0.403 57 股票代码 S11 S12 S13 S14 S15 S16 S17 S18 S19 S20 收益率均值 0.005 49 0.003 91 0.006 47 0.005 35 0.005 41 0.002 04 0.002 42 0.001 88 0.002 98 0.001 83 相关系数 0.249 43 0.161 08 0.384 82 0.223 37 0.382 54 0.308 39 0.496 61 0.338 87 0.345 69 0.469 35 股票代码 S21 S22 S23 S24 S25 S26 S27 S28 S29 S30 收益率均值 0.003 11 0.002 04 0.002 71 0.003 11 0.002 24 0.003 92 0.004 15 0.001 83 0.001 53 0.002 97 相关系数 0.524 21 0.334 95 0.396 09 0.438 01 0.510 86 0.390 27 0.431 78 0.331 44 0.033 17 0.406 50 3.1 计算结果
假设分段线性凸交易成本为
C(xi)={0.001⋅xi(0≤xi≤0.3),0.0015⋅xi−0.00015(0.3<xi≤1), 其中，i=1, …, 30，下界限制li=-0.5，上界限制ui=0.5 (i=1, …, 30)，无风险资产借款利率rb=0.000 9，贷款利率rl=0.000 4，无风险资产投资比例下限xn+10=-0.5，风险资产投资比例非负的个数不超过6，置信度c=95%.根据实际约束生成上述30只股票的非等比例投资组合.
本文考虑无风险资产借贷利率不同的情况下，根据模型(10)得到目标函数最大值，即r0max=0.006 6和VaR0max=0.09;当资产全部投资于无风险资产时，VaR0min=0，此时r0min=0.000 4.
当VaR0在区间[0, 0.09]之间变动时，投资组合的真实前沿面如图 3所示.
当VaR0在区间[0, 0.089]之间等距离变动时，净收益率变动情况见表 2.
表 2 投资组合的净收益率(VaR0∈[0, 0.089])Table 2. The portfolio return when VaR0 ∈[0, 0.089]VaR0 0 0.003 0.006 0.009 0.012 0.015 0.018 0.021 0.024 0.027 0.030 净收益率 0.000 40 0.000 66 0.000 93 0.001 19 0.001 44 0.001 69 0.001 94 0.002 19 0.002 43 0.002 67 0.002 91 VaR0 0.033 0.036 0.039 0.042 0.045 0.048 0.051 0.054 0.057 0.060 0.063 净收益率 0.003 15 0.003 39 0.003 63 0.003 86 0.004 09 0.004 32 0.004 54 0.004 76 0.004 98 0.005 18 0.005 38 VaR0 0.066 0.069 0.072 0.075 0.078 0.081 0.084 0.087 0.089 净收益率 0.005 57 0.005 75 0.005 92 0.006 07 0.006 21 0.006 33 0.006 43 0.006 51 0.006 55 由于交易成本、借款约束、上下界约束和基数约束等摩擦因素符合模型(9)的约束条件，所以，本文构造样本数据量m分别为100、200、500、1 000、2 000、4 000的随机投资组合，并将此投资组合的样本数据量代入收益导向的BCC-DEA模型(14)中，从而得到该模型的θ，并根据θ得到有效VaR.
根据表 2和表 3，可以得到前沿面的对比(图 4)，由图可以看出：(1)随着样本数据量m的增大，DEA前沿面逐渐接近于真实前沿面.根据DEA逼近性原理[25]，当样本数据量足够大时，所得到的前沿面与真实前沿面几近重合，因此，可以近似地把样本数据量足够大的DEA前沿面视为真实前沿面. (2)当风险值小于0.03时，m=200, 500, 1 000, 2 000, 4 000的DEA前沿面差距不大，甚至几乎重合，而m=100时前沿面较低；当风险值大于0.07时，m=2 000, 4 000的前沿面出现拐点并呈现水平状，而风险值大于0.075时，m=100, 200, 500, 1 000的前沿面也出现拐点并呈现水平状，并不能很有效地接近理论上的真实前沿面.
表 3 模型(14)的不同风险和样本数据量下的净收益率Table 3. The net return according to different risk and sample sizes in Model(14)VaR0 m/只 100 200 500 1 000 2 000 4 000 0 0.000 34 0.000 34 0.000 34 0.000 34 0.000 34 0.000 34 0.003 0.000 52 0.000 54 0.000 57 0.000 57 0.000 57 0.000 58 0.006 0.000 67 0.000 78 0.000 80 0.000 81 0.000 82 0.000 82 0.009 0.000 86 0.000 97 0.000 98 0.000 99 0.001 03 0.001 05 0.012 0.001 04 0.001 17 0.001 21 0.001 22 0.001 25 0.001 28 0.015 0.001 22 0.001 38 0.001 43 0.001 45 0.001 48 0.001 50 0.018 0.001 41 0.001 57 0.001 62 0.001 65 0.001 72 0.001 73 0.021 0.001 59 0.001 76 0.001 81 0.001 86 0.001 94 0.001 96 0.024 0.001 77 0.001 95 0.001 99 0.002 06 0.002 18 0.002 19 0.027 0.001 95 0.002 13 0.002 17 0.002 26 0.002 38 0.002 41 0.030 0.002 14 0.002 30 0.002 36 0.002 47 0.002 61 0.002 64 0.033 0.002 32 0.002 47 0.002 54 0.002 67 0.002 86 0.002 87 0.036 0.002 50 0.002 65 0.002 83 0.002 87 0.003 10 0.003 10 0.039 0.002 66 0.002 82 0.002 91 0.003 08 0.003 31 0.003 32 0.042 0.002 81 0.002 99 0.003 09 0.003 28 0.003 52 0.003 55 0.045 0.002 97 0.003 16 0.003 28 0.003 48 0.003 75 0.003 78 0.048 0.003 12 0.003 33 0.003 46 0.003 69 0.003 97 0.004 01 0.051 0.003 28 0.003 51 0.003 65 0.003 89 0.004 21 0.004 24 0.054 0.003 43 0.003 68 0.003 83 0.004 09 0.004 38 0.004 46 0.057 0.003 59 0.003 85 0.004 02 0.004 30 0.004 65 0.004 69 0.060 0.003 74 0.004 01 0.004 20 0.004 50 0.004 90 0.004 92 0.063 0.003 90 0.004 15 0.004 38 0.004 67 0.005 13 0.005 15 0.066 0.004 05 0.004 30 0.004 57 0.004 84 0.005 36 0.005 37 0.069 0.004 21 0.004 44 0.004 75 0.005 02 0.005 54 0.005 57 0.072 0.004 36 0.004 58 0.004 94 0.005 19 0.005 54 0.005 57 0.075 0.004 52 0.004 73 0.005 12 0.005 36 0.005 54 0.005 57 0.078 0.004 67 0.004 80 0.005 22 0.005 42 0.005 54 0.005 57 0.081 0.004 70 0.004 80 0.005 22 0.005 42 0.005 54 0.005 57 0.084 0.004 70 0.004 80 0.005 22 0.005 42 0.005 54 0.005 57 0.087 0.004 70 0.004 80 0.005 22 0.005 42 0.005 54 0.005 57 0.089 0.004 70 0.004 80 0.005 22 0.005 42 0.005 54 0.005 57 3.2 结果分析
针对上述计算结果所得到的结论，从样本数据量出发，对样本收益率均值、不同样本数据量下的净收益率的分布和收敛效果、相关性以及增加样本数据量的代价方面进行讨论.
(1) 样本收益率均值的分布分析.由随机选取的30只股票数据的实际收益率均值分布与正态收益率分布之间的差别(图 5)及统计所得到的偏度值(1.101 6)和峰度值(3.352 1)可以看出：至少有18只股票的收益率均值集中分布在[0.001 80, 0.003 80]区间内，占比60%以上，这使得在考虑交易成本、借款约束、上下界约束和基数约束等实际约束条件下，后期计算不同样本数据量的净收益时，净收益率会集中在这个区间.
(2) DEA前沿面的收敛效果分析.本文从3个角度验证不同样本数据量下的DEA前沿面的收敛效果.
其一，PEr、PEv和DEA效率值θ间的相关系数与收敛效果.由真实前沿面得到的收益导向投资组合效率PEr、风险导向投资组合效率PEv与DEA效率值θ间的相关系数(表 4)可知：随着样本数据量的增加，PEr、PEv与θ间的相关系数越来越大, 两者之间差值越来越小，而由图 4可知构造的前沿面越来越逼近真实前沿面，表明m越大，其收敛效果越好.
表 4 收益、风险导向效率PEr、PEv与DEA效率值θ间的相关系数Table 4. The correlation coefficents between PEr or PEv and θ比较因子 m/只 100 200 500 1 000 2 000 4 000 PEr与θ 0.944 3 0.964 2 0.981 0 0.993 5 0.994 0 0.994 4 PEv与θ 0.969 1 0.982 4 0.987 2 0.991 1 0.993 0 0.993 6 其二，收益率的差值与收敛效果.由不同风险和不同样本数据量下的净收益率与理想期望收益率的差值(表 5)可以看出：在样本数据量相同的情况下，随着风险水平的增加，净收益率与理想期望收益率的差值逐渐增大；在风险水平相同的情况下，随着样本数据量的增加，净收益率与理想期望收益率的差值逐渐缩小，其中，在m=2 000, 4 000时，此差值很小，且在2种风险水平下差值相同，可见，DEA前沿面有逐渐接近真实前沿面的趋势.
表 5 不同风险和不同样本数据量下的净收益率与理想期望收益率的差值Table 5. The distance value between the rate of net return and the rate of ideal expected return according to different risk and sample sizesm/只 VaR0 0.036 0.051 0.075 0.089 100 0.000 77 0.001 26 0.001 55 0.001 85 200 0.000 61 0.001 03 0.001 34 0.001 75 500 0.000 56 0.000 90 0.000 95 0.001 33 1 000 0.000 52 0.000 65 0.000 71 0.001 12 2 000 0.000 32 0.000 33 0.000 50 0.000 98 4 000 0.000 29 0.000 31 0.000 50 0.000 98 其三，净收益率区间与收敛效果.由表 3和图 4可知：当净收益率rN在区间[0, 0.003]区间内取值时，不同样本数据量下的分布区别不大，几乎重合；在[0.004, 0.006]区间内取值时，m=200, 500, 1 000, 2 000, 4 000的净收益率随着样本数据量的增加而逐渐增加，DEA前沿面逐渐接近真实前沿面；对比m=1 000与m=2 000、4 000的前沿面，随着样本数据量的增加，同一风险水平所对应的净收益值会增多，其对应的风险更为连续，即可以在更低的风险水平下达到最大值，这使得m=2 000、4 000时前沿面出现拐点并呈现水平状.
净收益率rN在m=2 000, 4 000时取得最大值，但仅为0.005 57，这与真实前沿面理论下的最大收益率(rN=0.006 55)差距较大，即图 4中DEA前沿面的后段与真实前沿面差距较大，并不能很有效地接近理论真实的前沿面.
(3) 不同样本数据量运行时间分析.本文采用Windows7旗舰版系统，CPU型号为Inter(R) Core(TM) i7-870，2.93 Hz.由程序的运行时间(表 6)可知：随着样本数据量的翻倍增加，运行时间呈几何式增加.在m=2 000, 4 000时，相同的风险水平下，其净收益率与理想期望收益率的差值非常小，而运行时间由8 502.8 s增加到35 470.8 s，增长了约4.2倍，这种时间的代价和模拟效果的不匹配，给工作效率带来巨大的困难.因此，在研究诸如多种限制条件下的前沿面拟合效率问题时，可以基于合理假设适当地修正投资组合对应取值的股票投资比例分布，以减少工作量，提高工作效率.
表 6 不同样本数据量下的运行时间Table 6. The operation time for different sizes of sample datas m/只 100 200 500 1 000 2 000 4 000 运行时间 24.8 62.7 367.7 1 070.8 8 502.8 35 470.8 4. 结论
考虑交易成本、借款约束、上下界约束和基数约束等实际约束条件，本文提出了均值-VaR投资组合优化模型.在投资组合真实前沿面为凹函数的前提下，将BCC-DEA模型作为投资组合的绩效评价方法，根据逼近原理，将DEA前沿面与真实前沿面做对比.结果显示：(1)尽管投资组合样本数据量不同时得到的前沿面的估计值也不同，但随着样本数据量的增大，DEA前沿面逐渐接近于真实前沿面. (2)当样本数据量足够大时，所得到的前沿面与真实前沿面几近重合，可以近似地把样本数据量足够大的DEA前沿面作为真实前沿面. (3)在逼近过程中，原有的股票样本组合的收益率均值分布呈现出不均匀性，使得不同样本数据量下的DEA模型整体净收益较大值的出现频度进一步减少，整体净收益较小的值差距不大且频度较大，从而导致所得前沿面与真实前沿面差距较大，不能有效地接近真实前沿面.
另外，由收益、风险导向效率PEr、PEv与DEA效率值θ间的相关系数可知：(1)随着样本数据量的增大，两者的相关系数越来越高，证明了本文构造的DEA模型的有效性，样本数据量和逼近效果的关系证明了其收敛性. (2)前沿面拐点较早地出现，DEA前沿面靠近真实前沿面的趋势变慢也进一步证明了收敛效果.实证研究表明：当样本数据量足够大的时候，构造的样本组合前沿面可以拟合达到理论最优的真实前沿面；随着数据的增加，程序耗时呈几何式的增长，对工作效率影响较大.
-
![]()
图 3 李子坪和富隆厂铅锌矿床矿石手标本照片
注：Sp-闪锌矿；Grn-细硫砷铅矿；Anh-硬石膏；Orp-雌黄；Cal-方解石；Bit-沥青；Jor-灰硫砷铅矿；Gn-方铅矿。
Figure 3. Hand-specimen photos of the Liziping and Fulongchang deposits
![]()
图 4 李子坪和富隆厂矿床矿石光学显微镜下照片
注：Sp代表闪锌矿，Grn代表细硫砷铅矿，Anh代表硬石膏，Py代表黄铁矿，Bnn代表车轮矿，Ccp代表黄铜矿, Gn代表方铅矿, Brt代表重晶石, Cal代表方解石, Jor代表灰硫砷铅矿。
Figure 4. Reflected-light microphotographs for the Liziping and Fulongchang deposits
![]()
图 5 李子坪矿床典型样品中脉石矿物的XRD图谱
Figure 5. XRD patterns of gangue minerals in typical samples of the Liziping deposit
![]()
图 6 李子坪和富隆厂矿床部分脉石和含铅矿物的BSE图像和EDS谱图
Figure 6. BSE images and EDS spectra of some gangues and lead-bearing minerals in the Liziping and Fulongchang deposits
![]()
![]()
图 9 李子坪矿床成矿前黄铁矿和成矿期胶状含砷黄铁矿中As、Pb、Cu、Ag质量分数的变化
Figure 9. Changes of As, Pb, Cu and Ag contents in brecciated pyrite and colloidal arsenic bearing pyrite in the Liziping deposit
![]()
图 10 李子坪矿床胶状含砷黄铁矿核部到边缘的元素变化
Figure 10. Element changes from core to margin of colloidal arsenic bearing pyrite in Liziping deposit
表 1 李子坪矿床中部分典型样品的XRD矿物学组成
Table 1 XRD mineral compositions of typical samples from the Liziping deposit
样品号 矿物类型 矿物组成及其含量占比(%) LZP-1-2 脉石矿物 白云石(95%)、石英(2%)、方解石(2%)、高岭石(1%) LZP-1-3 脉石矿物 重晶石(35%)、硬石膏(63%)、石膏(2%) LZP-1-4 脉石矿物 硬石膏(93%)、石膏(7%) LZP-1-5 脉石矿物 方解石(32%)、白云石(19%)、硬石膏(39%)、石膏(10%) LZP-2-2 脉石矿物 重晶石(26%)、硬石膏(73%)、石膏(2%) LZP-2-3 脉石矿物 硬石膏(98%)、石膏(2%) LZP-3-2 脉石矿物 硬石膏(84%)、重晶石(9%)、石英(5%)、石膏(2%) LZP-3-3 脉石矿物 重晶石(46%)、硬石膏(50%)、石膏(4%) LZP-3-4 脉石矿物 重晶石(61%)、硬石膏(36%)、石膏(4%) LZP-4-5 脉石矿物 硬石膏(97%)、石膏(3%) LZP-4-7 脉石矿物 硬石膏(87%)、方解石(12%)、石膏(1%) LZP-4-8 脉石矿物 硬石膏(96%)、石膏(4%) LZP-4-9 脉石矿物 硬石膏(100%) LZP-4-11 脉石矿物 硬石膏(100%) LZP-5-3 脉石矿物 硬石膏(97%)、石膏(3%) LZP-5-6 脉石矿物 硬石膏(95%)、石膏(5%) LZP-7-2 脉石矿物 硬石膏(100%) LZP-7-3 脉石矿物 硬石膏(92%)、石膏(8%) LZP-8 脉石矿物 硬石膏(94%)、石膏(6%) 
下载: 导出CSV
表 2 不同矿物的元素含量
Table 2 Element content of different minerals
矿物 元素 质量分数/% 摩尔分数/% 矿物 元素 质量分数/% 摩尔分数/% 细硫砷铅矿 As 11.97 14.70 灰硫砷铅矿 As 7.59 9.70 S 19.04 54.60 S 17.97 53.57 Pb 68.98 30.64 Pb 67.37 31.12 Sb 7.13 5.61
下载: 导出CSV
表 3 李子坪矿床矿石矿物电子探针分析结果
Table 3 EPMA data of ore minerals from the Liziping deposit
% 点位编号 矿物 S Pb Sb As Fe Hg Ag Cu Zn Au 总计 LZP-10.1 细硫砷铅矿 19.2 69.7 0.875 9.39 0 0 0 0 0 0 99.2 LZP-10.4 细硫砷铅矿 20.0 70.5 0.907 9.64 0.034 0.020 0 0 0.228 0 101.3 LZP-10.2 细硫砷铅矿 20.3 71.8 1.552 9.71 0 0 0 0 0 0.018 103.4 LZP-10.3 细硫砷铅矿 19.6 70.2 1.814 9.30 0.021 0.045 0.062 0 0.519 0 101.5 LZP-10.5 细硫砷铅矿 19.8 69.8 0.595 9.99 0.032 0 0.060 0 0.076 0 100.4 LZP-8-2.2 细硫砷铅矿 19.3 70.6 0.076 10.10 0.041 0.116 0 0.023 0 0.017 100.3 LZP-8-2.4 细硫砷铅矿 18.8 67.2 0.249 9.45 0.061 0 0 0 0.029 0 95.8 LZP-8-2.5 细硫砷铅矿 20.1 70.8 0.288 10.30 0.026 0 0 0 0.064 0 101.6 LZP-8-2.6 细硫砷铅矿 20.4 69.4 0.376 11.13 0.033 0.092 0.065 0.016 0 0 101.4 LZP-3-1.6 细硫砷铅矿 20.0 71.1 0.405 9.56 0 0.061 0.057 0.022 0 0 101.2 LZP-3-1.5 细硫砷铅矿 19.6 69.4 0.314 9.62 0.024 0.045 0 0 0.299 0.018 99.3 LZP-3-1.3 细硫砷铅矿 20.3 70.6 0.512 9.93 0 0 0 0 0 0 101.3 LZP-3-1.2 细硫砷铅矿 20.3 71.4 0.255 10.50 0.105 0.068 0 0.011 0 0 102.7 LZP-3-1.1 细硫砷铅矿 20.2 71.6 0.252 10.10 0.036 0.048 0.025 0 0 0 102.2 LZP-1-4.1 细硫砷铅矿 19.7 68.9 0.961 9.44 0.047 0 0 0.007 0 0 99.1 LZP-1-4.2 细硫砷铅矿 19.8 69.9 0.592 10.10 0 0 0 0.003 0 0 100.4 LZP-1-4.3 细硫砷铅矿 19.8 67.8 0.514 9.58 0.017 0.009 0 0 0.325 0 98.1 LZP-1-4.4 细硫砷铅矿 20.0 69.5 0.447 9.96 0 0.080 0 0.030 0 0 100.0 LZP-4-1.1 细硫砷铅矿 19.1 69.5 1.895 9.26 0.034 0.061 0 0 0 0.029 99.9 LZP-4-1.2 细硫砷铅矿 19.8 69.1 0.488 9.72 0.031 0 0 0 0 0 99.1 LZP-4-1.3 细硫砷铅矿 19.8 71.3 0.131 10.20 0.018 0 0 0.021 0 0 101.4 LZP-4-3.1 细硫砷铅矿 19.5 70.6 0.476 10.20 0.037 0 0 0 0.090 0 100.9 LZP-4-3.3 细硫砷铅矿 19.8 69.8 1.558 9.50 0.099 0.142 0 0 0 0 103.0 LZP-4-3.5 细硫砷铅矿 19.8 69.1 0.208 11.90 0.550 0.152 0.001 0.004 0 0 101.7 LZP-4-3.6 细硫砷铅矿 19.6 69.6 1.525 9.38 0 0 0 0 0 0.042 100.1 LZP-8-2.3 方铅矿 13.3 84.0 0.410 0 0.033 0.054 0 0 0 0 97.8 LZP-8-2.6a 方铅矿 13.3 81.8 0.608 0 0.057 0.090 0 0.021 0 0 95.8 LZP-3-1.4 方铅矿 14.3 83.8 0.678 0 0.025 0 0 0.032 0 0 98.9
下载: 导出CSV
表 4 富隆厂矿床矿石矿物电子探针分析结果
Table 4 EPMA data of ore minerals from Fulongchang deposit
% 点编号 矿物 S Pb Sb As Fe Hg Ag Cu Zn Au 总计 DPZ-29.6 灰硫砷铅矿 18.8 69.9 4.01 7.555 0.002 0.101 0.118 0 0 0.002 100.6 DPZ-18.1 灰硫砷铅矿 18.9 68.3 3.43 8.005 0.016 0 0.003 0 0 0.043 98.7 DPZ-18.2 灰硫砷铅矿 19.0 69.1 3.34 7.909 0 0 0 0 0 0.004 99.4 DPZ-18.3 灰硫砷铅矿 18.6 68.0 6.36 5.343 0.005 0.002 0 0.005 0.263 0.002 98.6 DPZ-18.4 灰硫砷铅矿 19.0 68.2 3.81 7.138 0.048 0 0 0 0.058 0 98.3 DPZ-18.5 灰硫砷铅矿 19.2 69.8 3.49 8.264 0.031 0.095 0 0 0.395 0 101.3 DPZ-29.1 灰硫砷铅矿 18.8 69.1 3.69 7.905 0.010 0.023 0.135 0.020 0 0 99.7 DPZ-29.7 灰硫砷铅矿 19.0 70.1 3.25 8.641 0.017 0.081 0 0.005 0 0 101.1 DPZ-16.1 灰硫砷铅矿 19.0 70.9 3.78 6.808 0.013 0.191 0 0 1.560 0.033 102.3 DPZ-16.3 灰硫砷铅矿 19.8 70.2 3.72 8.334 0.031 0.206 0 0 0.024 0 102.3 DPZ-16.4 灰硫砷铅矿 19.5 69.0 2.95 8.217 0.031 0.099 0 0.008 0 0 99.8 DPZ-23.3 灰硫砷铅矿 19.7 68.8 5.48 6.846 0.024 0.175 0.161 0 0.241 0.022 101.4 DPZ-14-1.6 灰硫砷铅矿 19.5 66.9 5.65 6.443 0 0.018 0.127 0 0.122 0.040 98.9 DPZ-29.2 方铅矿 13.0 85.3 0.755 0 0.002 0.029 0.011 0.008 0 0.011 99.1 DPZ-29.3 方铅矿 13.2 85.2 0.831 0 0 0.110 0.005 0 0 0 99.3 DPZ-29.4 方铅矿 13.3 84.8 0.499 0 0 0.052 0 0 0 0 98.6 DPZ-16.2 方铅矿 13.7 85.6 0.201 0 0.002 0.127 0 0 0 0 99.6 DPZ-16.5 方铅矿 13.7 86.3 0.234 0 0.018 0.007 0 0 0.803 0 101.0 DPZ-8.1 方铅矿 13.9 84.7 0.913 0 0.020 0.098 0.067 0 0.289 0 100.0 DPZ-8.2 方铅矿 12.3 82.2 0.165 0 0.332 0.046 0 0 0.882 0 95.9 DPZ-8.3 方铅矿 13.2 83.4 0.678 0 0.046 0 0.091 0 0.046 0 97.9 DPZ-23.5 方铅矿 13.5 83.6 1.010 0 0.013 0.117 0 0 0.004 0 98.3 DPZ-23.4 方铅矿 14.1 83.6 1.810 0 0.024 0 0 0.003 0 0 99.6 DPZ-23.2 方铅矿 13.7 85.9 0.810 0 0.007 0.056 0 0 0.354 0 100.8 DPZ-23.1 方铅矿 13.8 87.4 0.321 0 0.053 0.057 0 0 1.770 0 103.3 DPZ-14-1.3 方铅矿 13.6 85.1 0.652 0 0.012 0 0.150 0.030 0 0.045 99.6 DPZ-14-1.5 方铅矿 13.6 85.1 0.770 0 0.041 0.046 0.337 0.002 0 0 99.9 DPZ-14-1.7 车轮矿 23.8 42.8 16.0 5.624 0.017 0 0 13.0 0 0.006 101.3 DPZ-14-1.4 车轮矿 24.3 43.9 11.9 8.229 0.018 0.140 0 13.3 0 0.010 101.8 DPZ-14-1.2 车轮矿 23.9 44.6 14.1 7.142 0.019 0.027 0 13.5 0 0 103.3
下载: 导出CSV
-
[1] 田洪亮. 兰坪白秧坪铜银多金属矿床地质特征[J]. 云南地质, 1997, 16: 105-108. [2] 胡茂德, 张洪瑞, 贾敬伍, 等. 云南兰坪李子坪铅锌矿床的控矿构造、碳和氧同位素及稀土元素地球化学特征[J]. 矿床地质, 2015, 34: 1057-1071. [3] XUE C J, ZENG R, LIU S W, et al. Geologic, fluid inclusion and isotopic characteristics of the Jinding Zn-Pb deposit, western Yunnan, South China: a review[J]. Ore Geology Reviews, 2007, 31: 337-359. doi: 10.1016/j.oregeorev.2005.04.007
[4] 薛伟, 薛春纪, 池国祥, 等. 滇西北兰坪盆地白秧坪多金属矿床流体包裹体研究[J]. 岩石学报, 2010, 26: 1773-1784. [5] WANG X H, HOU Z Q, SONG Y C, et al. Geological, fluid inclusion and isotopic studies of the Baiyangping Pb-Zn-Cu-Ag polymetallic deposit, Lanping basin, Yunnan province, China[J]. Journal of Asian Earth Sciences, 2015, 111: 853-871. doi: 10.1016/j.jseaes.2015.08.009
[6] 刘家军, 翟德高, 李志朋, 等. 兰坪盆地白秧坪银铜多金属矿集区中银、钴、铋、镍的赋存状态与成因意义[J]. 岩石学报, 2010, 26: 1646-1660. [7] 王晓虎, 侯增谦, 宋玉财, 等. 兰坪盆地白秧坪铅锌铜银多金属矿床成矿流体及成矿物质来源[J]. 地球科学(中国地质大学学报), 2012, 37: 1015-1028. [8] 侯增谦, 宋玉财, 李政, 等. 青藏高原碰撞造山带Pb-Zn-Ag-Cu矿床新类型: 成矿基本特征与构造控矿模型[J]. 矿床地质, 2008: 27(2): 123-144. doi: 10.3969/j.issn.0258-7106.2008.02.001 [9] ZOU Z C, HU R Z, BI X W, et al. Noble gas and stable isotopic constraints on the origin of the Ag-Cu polymetallic Ore deposits in the Baiyangping Area, Yunnan Province, SW China[J]. Resource Geology, 2016, 66: 183-198. doi: 10.1111/rge.12096
[10] 王长明, 陈晶源, 杨立飞, 等. 三江特提斯兰坪盆地构造-流体-成矿系统[J]. 岩石学报, 2017, 33: 1957-1977. [11] 宋玉财, 侯增谦, 杨天南, 等. "三江"喜马拉雅期沉积岩容矿贱金属矿床基本特征与成因类型[J]. 岩石矿物学杂志, 2011, 30(3): 355-380. doi: 10.3969/j.issn.1000-6524.2011.03.002 [12] DENG J, WANG Q F, LI G J, et al. Cenozoic tectono-magmatic and metallogenic processes in the Sanjiang region, southwestern China[J]. Earth-Science Reviews, 2014, 138: 268-299. doi: 10.1016/j.earscirev.2014.05.015
[13] 毕献武, 唐永永, 陶琰, 等. 西南三江碰撞造山带沉积岩容矿Pb-Zn-Ag-Cu贱金属复合成矿与深部过程[J]. 岩石学报, 2019, 35: 1341-1371. doi: 10.18654/1000-0569/2019.05.03 [14] 曾荣. 兰坪盆地流体大规模成矿过程[D]. 西安: 长安大学, 2007. [15] 孔云丽, 戚林坤. 云南兰坪中排李子坪铅锌矿矿床成因[J]. 云南地质, 2009, 28(3): 275-279. doi: 10.3969/j.issn.1004-1885.2009.03.009 [16] 邓霜岭. 云南李子坪铅锌矿段地质特征及成因[J]. 四川地质学报, 2011, 31(3): 323-325;328. doi: 10.3969/j.issn.1006-0995.2011.03.018 [17] LEACH D, SANGSTER D, KELLEY K, et al. Sediment-hosted lead-zinc deposits: a global perspective[J]. Economic Geology, 2005, 100: 561-607.
[18] 王晓虎. 兰坪盆地白秧坪铅锌铜银多金属矿床地质与成因[D]. 中国地质科学院, 2011. [19] 梁维, 侯增谦, 杨竹森, 等. 藏南扎西康大型铅锌银锑多金属矿床叠加改造成矿作用初探[J]. 岩石学报, 2013, 29(11): 3828-3842. [20] JOHAN Z, MANTIENNE J. Thallium-rich mineralization at Jas Roux, Hautes-Alpes, France: a complex epithermal, sediment-hosted, Ore-forming system[J]. Journal of the Czech Geological Society, 2000, 45: 63-77.
[21] WALIA D S, CHANG L L Y. Investigations in the systems PbS-Sb2S3-As2S3 and PbS-Bi2S3-As2S3[J]. The Canadian Mineralogist, 1973, 12: 113-119.
[22] ROLAND G W. The system Pb-As-S[J]. Mineralium Deposita, 1968, 3: 249-260.
计量
- 文章访问数: 64
- HTML全文浏览量: 12
- PDF下载量: 55
