• 中文核心期刊
  • 中国科技核心期刊
  • CSCD来源期刊
  • DOAJ 收录
  • Scopus 收录

热脉冲探针-时域反射技术测量含水合物沉积物的热导率及水合物饱和度

陈强, 刁少波, 孙建业, 刘昌岭, 业渝光

陈强, 刁少波, 孙建业, 刘昌岭, 业渝光. 热脉冲探针-时域反射技术测量含水合物沉积物的热导率及水合物饱和度[J]. 岩矿测试, 2013, 32(1): 108-113.
引用本文: 陈强, 刁少波, 孙建业, 刘昌岭, 业渝光. 热脉冲探针-时域反射技术测量含水合物沉积物的热导率及水合物饱和度[J]. 岩矿测试, 2013, 32(1): 108-113.
Qiang CHEN, Shao-bo DIAO, Jian-ye SUN, Chang-ling LIU, Yu-guang YE. Measurement of Thermal Conductivity and Saturation of Gas Hydrates in Sediment by Thermal Pulse Probe-Time Domain Reflection Technique[J]. Rock and Mineral Analysis, 2013, 32(1): 108-113.
Citation: Qiang CHEN, Shao-bo DIAO, Jian-ye SUN, Chang-ling LIU, Yu-guang YE. Measurement of Thermal Conductivity and Saturation of Gas Hydrates in Sediment by Thermal Pulse Probe-Time Domain Reflection Technique[J]. Rock and Mineral Analysis, 2013, 32(1): 108-113.

热脉冲探针-时域反射技术测量含水合物沉积物的热导率及水合物饱和度

基金项目: 

国家重点基础研究发展计划(973计划)项目(2009CB219503)

国家重点基础研究发展计划(973计划)项目 2009CB219503

详细信息
    作者简介:

    陈强,博士,主要从事天然气水合物模拟实验研究。E-mail:chenqiang_hds@126.com

    通讯作者:

    业渝光,研究员,主要从事天然气水合物模拟实验研究。E-mail: yeyug@public.qd.sd.cn

  • 中图分类号: TE122;P575.6

Measurement of Thermal Conductivity and Saturation of Gas Hydrates in Sediment by Thermal Pulse Probe-Time Domain Reflection Technique

  • 摘要: 我国在海洋和冻土区都已发现天然气水合物资源区并成功获取实物样品。含水合物沉积物的热导率是估算水合物资源量、设计合理开采方案的关键性数据之一。受水合物稳定条件和测量技术的限制,水合物热导率测定尚不完善。本文通过自主研制的天然气水合物热物理参数测量系统,开展了海洋沉积物中天然气水合物热导率与饱和度测量研究。实验使用取自南海神狐海域的沉积物作为反应介质,在压力7.8 MPa、温度2℃的条件下合成甲烷水合物,并利用热脉冲探针与时域反射技术联合测量的方式获得沉积物中水合物形成过程的热导率和饱和度等实验数据。结果表明,当水合物饱和度从0增加至49%时,体系热导率出现了先升高后降低的变化趋势。分析发现体系热导率随水合物饱和度的变化特征与水合物在沉积物中的填充方式有关,在实验选用的南海沉积物中,水合物优先选择在颗粒孔隙间成核生长,并最终与沉积物颗粒胶结共存。
  • 胶东是中国最重要的金矿集区,是一个主要由前寒武纪基底岩石和超高压变质岩块组成、中生代构造-岩浆作用发育的内生热液金矿成矿集中区,其金矿床形成时代和产出背景在全球造山型金矿中独一无二。山东莱州-招远金矿区是胶东金矿集区的主要组成部分,拥有玲珑、焦家等世界级金矿。其中,焦家金矿带是莱州-招远金矿区最重要的金矿带之一。该带内已发现了3个特大型金矿床及一批大中型金矿床,累计探明金储量超过500吨。焦家断裂带在寺庄金矿区内处于与成矿关系十分密切的玲珑黑云母花岗岩内及胶东群黑云斜长角闪片麻岩的接触带上,产状复杂且变化较大[1]。由于断裂带产状的变化和地质背景的复杂化,焦家金矿控矿构造系统也随之出现了明显的多样性[2]。由于该区在采矿方面存在的科学疑难问题多、找矿难度大、隐伏矿为主、物化探异常干扰显著等原因,前期对该区矿石的工艺矿物学研究以常规为主,李德亭等[3-4]对焦家金矿深部矿石矿物种类、载金矿物种类及金赋存状态等进行了初步研究,研究表明矿物主要有黄铁矿、黄铜矿、方铅矿等硫化矿和石英、绢云母、钾长石等脉石矿物;黄铁矿、黄铜矿为主要载金矿物,金常以裂隙方式充填于黄铁矿、黄铜矿等硫化矿中。现有成果对金赋存状态特征、金矿物种类、不同大小金粒的成色等研究涉及很少,且缺乏先进测试手段进行综合表征。

    焦家金矿的主矿区是原矿主要供应区,该矿区供应给焦家金矿的矿量占到总供应量的70%。本文在主矿区采集深采矿石,采用偏光反光两用显微镜、扫描电镜[5-10]、X射线衍射仪[11-13]以及能谱分析[14-16]研究工艺矿物学特征,查清载金矿物种类及含量,进而分析金在主要载金矿物中的赋存状态及特征,同时研究金的形状、金粒度大小及分布特点,对不同粒度的金成色、金矿物类型等进行系统分析,研究成果丰富了焦家金矿矿物学研究的内容,为后续选冶工艺提供了调控依据和重要信息。

    LEICA-DMLP高级研究型偏光反光两用显微镜(德国LEICA公司):主要附件LEICA MPS30照相系统,荧光附件。

    LINKAM热台THMSE600(-196~600℃):配备工艺矿物学研究的高级软件LINKSYS。

    LEICA-MZ6高级研究型体视显微镜(德国LEICA公司):该设备自带全套摄像、自动照相、颗粒分析软件。工作条件为自动对焦;固定倍率切换9段;放大倍率15~100倍;视野范围3.3~41.2 mm;工作距离100 mm;灯源LED环形灯;物镜0.63~4 X。

    SSX-550型扫描电镜及其附带的DX-4能谱仪(日产岛津公司):工作条件为二次电子图像分辨率3.5 nm,放大倍率20~30000倍。

    Empyrean型X射线衍射仪(XRD,荷兰PANALYTICALB.V公司):配合显微镜分析鉴定矿物种类,工作条件为功率3 kW,测角仪重现性0.0001°,测角仪类型T-2T。

    对矿区井下300 m、400 m、500 m处合计取样150 kg,选取有代表性的矿样按要求制作成120块抛光片和30块矿石薄片。采用不同仪器进行观察研究。

    通过光片研究并结合化学分析及化学物相分析等方法,查明焦家金矿主矿区矿物种类主要为硫化矿,以黄铁矿、黄铜矿为主,还有少量方铅矿、闪锌矿、磁黄铁矿、黝铜矿、辉铋矿、辉钼矿。金属氧化物主要为褐铁矿及少量磁铁矿等。脉石有石英、长石、绢云母、方解石和绿泥石类矿物等;还有极小量的菱铁矿、石榴石、金红石、锆石、榍石、磷灰石等矿物。

    采用高级研究型偏光反光两用显微镜对120块抛光片研究,发现载金矿物主要为硫化物,如黄铁矿、黄铜矿、闪锌矿和方铅矿;载金脉石矿物主要是石英和长石。对30块薄片进行研究,统计发现硫化物中金粒有152粒,脉石矿物中金粒22颗。硫化物中金粒占87.37%,脉石中金粒占12.63%。金在各矿物中赋存含量如图 1所示,黄铁矿、黄铜矿是重要载金矿物。66.25%的金赋存在黄铁矿中,19.29%的金赋存在黄铜矿中,11.52%的金赋存在石英中,其他矿物及脉石中含金量很少。

    图  1  载金矿物及脉石中金的含量
    Figure  1.  Gold contents in different minerals

    黄铁矿是最重要的载金矿物。在矿石中呈自形粒状、半自形粒状和它形粒状,集合体呈致密块状、粒状或结核状、脉状等。浅黄(铜黄)色,条痕绿黑色,强金属光泽,不透明,无解理,参差状断口,在地表条件下易风化为褐铁矿。黄铁矿与黄铜矿密切共生,成矿期多见黄铁矿包裹黄铜矿。黄铁矿的裂隙往往被黄铜矿充填,表明时间上的同期性。黄铁矿的嵌布粒度最小为0.001 mm,最大为3.9 mm,一般分布于0.050~0.10 mm范围内。黄铜矿也是重要的载金矿物。黄铜矿呈现黄铜色,表面常由于氧化产生斑驳的蓝、紫、褐色的锖色,强金属光泽。黄铜矿的嵌布粒度粗细不均匀。黄铜矿呈它形粒状,以细粒、微细粒状嵌布在脉石中,黄铜矿的嵌布粒度最小为0.001 mm,最大为1.21 mm,一般分布于0.030~0.080 mm范围内。

    对主要载金矿物黄铁矿、黄铜矿、石英和长石等中的金赋存状态进行分析。

    通过对120块磨制光片进行高级研究型偏光反光两用显微镜及扫描电镜观察,查明焦家金矿中金的赋存状态主要是裂隙金、包裹金和晶隙金。其中裂隙金占金矿物总量64.82%。裂隙金(图 2)成群网状分布于硫化物裂隙或硫化物和石英裂隙中。对图 2(a)中的裂隙金(1)、(2)点进行能谱分析(图 3)。从图 3(a)可知,Au的M线在2.120 keV处出现,Ag的L线在2.992 keV处出现,Fe的K线在6.405 keV处出现,可知裂隙金含有少量的Fe,另外还发现部分裂隙金中含有Cr、Cu等,裂隙金银矿物成分比较复杂;由图 3(b)可知,石英裂隙中的金不含杂质元素,表明在热液脉Au-SiO2封闭体系中不含其他元素,二氧化硅溶胶聚集形成凝胶,金溶胶可在凝胶内呈浸染状固定下来,后来随着凝胶结晶脱水作用,使金以微粒分散或溶解状态作进一步的迁移,在二氧化硅冷却形成石英期间,被圈闭的金溶胶可以发生扩散作用和聚集作用,结果形成颗粒裂隙金。有些裂隙金含有Fe、Cr、Cu等,有些又不含任何微量元素,反映成矿热液是多期次的,金银矿物结晶也是多期次的。

    图  2  裂隙金扫描电子图片
    Figure  2.  The SEM image of fissure gold
    图  3  裂隙金能谱图
    Figure  3.  The energy spectrum images of fissure gold

    通过120块磨制光片进行电子扫描显微镜观察,金赋存状态中包裹金占金矿物总量19.29%。用电子扫描显微镜研究金颗粒的生长和金成色情况。光片显示小颗粒金分布广且多。图 4是三颗包裹金的扫面电镜图片,较大颗粒(图 4中1处)被黄铜矿包裹,金颗粒中含有孔洞,这些孔洞结构有利于溶液渗透,为氰化溶金提供了良好的条件。图 4中2处是被黄铁矿包裹的小粒金,颗粒粒度为14 μm,小部分边界与石英接触;图 4中3处是被黄铁矿包裹的小粒金,颗粒粒度为12 μm,对该晶体进行了Au、Ag和Cu的线扫描(图 5)。测定Au、Ag和Cu相对的含量。从图 5可知,金边部成色低,银略高,也有微量铜。中心金高,且有微量铜。

    图  4  包裹金的扫描电镜图片
    Figure  4.  The SEM image of inclusion gold
    图  5  包裹金的Au、Ag和Cu线扫描图
    Figure  5.  The Line scanning images of Au, Ag and Cu in inclusion gold

    晶隙金类型比较多,占金矿物总量的15.89%。对磨制光片进行电子扫描显微镜观察,晶隙金经常出现在黄铁矿/黄铜矿晶隙(占52.12%)、黄铁矿/石英晶隙(占28.45%)、黄铜矿/石英晶隙(占13.20%)及黄铜矿/斑铜矿晶隙(占6.23%)。图 6是黄铁矿/石英晶隙间的金矿物,晶隙金成群分布。能谱分析可知:金矿物晶体中心(图 6中1处):含Ag 17.544%,Au 82.456%,属于自然金;金矿物晶体边部(图 6中2处):含Ag 13.453%,Au 86.547%,属于自然金。小颗粒金矿物(8 μm,图 6中3处)赋存于黄铁矿和石英晶隙间。能谱分析表明:小颗粒金矿物含Ag 15.949%,Au 82.547%,Fe 1.504%。能谱测试粒度为3~8 μm的7颗小金颗粒,表明小粒金矿物成色高,为830.62‰~898.95‰,属于自然金。说明焦家金矿成矿期间的溶液十分丰富,在高温时是一个相,随着温度降低,金先沉淀,而银后沉淀,固溶体分离,而随着结晶到后期,形成尾部和尾梢,金的成分发生变化,形成的金颗粒中银含量幅度变化较大。

    图  6  晶隙金扫描电子图片
    Figure  6.  The SEM image of intercrystalline gold

    研究结果表明焦家金矿金赋存状态主要是裂隙金、包裹金和晶隙金。裂隙金最多占金矿物总量64.82%,其次是包裹金和晶隙金。裂隙金矿物成分复杂,有时含Fe、Cr、Cu等,包裹金体内局部有孔洞,金边部成色低,晶隙金类型丰富,形状多变,多属自然金。

    用高级研究型偏光反光两用显微镜观察,硫化物黄铁矿、黄铜矿中金颗粒较大,连群分布占多数,孤立分布较少[17-18]。脉石中金颗粒细小,孤立分布占绝大多数,极少见连群分布。金矿物粒度范围较大,大颗粒可达到90~110 μm,小颗粒只有2~3 μm。其中0.104~0.147 mm的金颗粒占1.85%,0.074~0.104 mm金颗粒占4.25%,0.043~0.074 mm金颗粒占16.24%,0.037~0.043 mm金颗粒占30.58%,小于0.037 mm金颗粒占大多数,占47.08%。不同载金矿物中金粒的粒度含量见表 1。从表 1可知细粒级载金矿物中的金含量较多,如小于0.037 mm的石英中金相对含量达56.01%,粗粒中金含量较少,如黄铁矿中0.104~0.147mm的金相对含量只有3.67%。

    表  1  载金矿物中金不同粒度的含量
    Table  1.  The gold percentage characteristics in different minerals
    载金
    矿物
    不同金粒度的相对含量w(Au)/%
    0.104~0.147
    mm
    0.074~0.104
    mm
    0.043~0.074
    mm
    0.037~0.043
    mm
    <0.037
    mm
    合计
    黄铁矿3.6710.3827.1840.7718.00-
    黄铜矿03.3415.0130.2351.42-
    石英0.912.0010.8930.1956.01-
    其他007.4928.5064.01-
    含量1.854.2516.2430.5847.08100
    下载: 导出CSV 
    | 显示表格

    用高级研究型偏光反光两用显微镜和高级研究型体视显微镜观察,焦家金矿主矿区金矿物形状主要有球形(占40.25%)、三角形(占9.28%)、矩形(占6.91%)、条形(占12.14%)、块形(占8.24%)、钩形(占11.85%)、不规则形(占11.33%)金等。金的典型形状见图 7

    图  7  金赋存典型形状图
    Figure  7.  Typical photos of gold occurrence

    用高级研究型偏光反光两用显微镜观察,结合XRD衍射仪测试结果,发现金矿物种类较丰富,有自然金、银金矿、自然银、金铜矿、含铁自然银、螺硫银矿、金银碲化物等。金矿物类型与含量见表 2。其中金银系列矿物占87.55%,金铜矿占6.58%,螺硫银矿占3.31%,碲化金矿占2.56%。金银矿物成色统计见表 3。从表 3可知,金银矿物平均成色较高,其中自然金达950.21‰,银金矿达738.92‰,金银系列矿物加权成色达到728.88‰。银金矿在金银矿系列含量中的比例(分配率)占到79.70%,是最重要的金矿物,其次为自然金,分配率占14.68%。

    表  2  金矿物类型与含量
    Table  2.  The gold mineral types and their contents
    金银矿物系列矿物名称w/%
    金银系列自然金12.85
    银金矿69.78
    自然银3.67
    含铁自然银1.25
    小计87.55
    金铜系列金铜矿6.58
    硫化物系列螺硫银矿3.31
    碲化物系列碲化金银2.56
    合计-100.00
    下载: 导出CSV 
    | 显示表格
    表  3  金银矿物成色和含量
    Table  3.  The gold percentages and contents of gold-silver minerals
    矿物名称平均成色/‰w/%分配率/%加权成色/‰
    自然金950.2112.8514.68728.88
    银金矿738.9269.7879.70
    自然银11.373.674.19
    含铁自然银01.251.430
    合计-87.55100.00-
    下载: 导出CSV 
    | 显示表格

    用高级研究型偏光反光两用显微镜观察,结合能谱分析测定不同粒度金矿物成色。由表 4的分析结果可知,60 μm和32 μm的大颗粒金成色分别为543.75‰和549.38‰,14 μm和12 μm的小颗粒金成色分别是859.72‰和856.15‰,说明大粒金成色低,为银金矿;小颗粒金与大颗粒金相比成色明显较高,小颗粒金大都为自然金。

    表  4  金矿物不同晶体颗粒的成色
    Table  4.  The gold percentages of different size crystal particles
    金颗粒粒度/μmw(Ag)/‰w(Au)/‰金含量水平定名
    60456.25543.75金低银金矿
    32450.62549.38金低银金矿
    14140.28859.72金高自然金
    12143.85856.15金高自然金
    下载: 导出CSV 
    | 显示表格

    以往研究成果对焦家金矿深部矿石的矿物种类、载金矿物种类及金赋存状态等进行工艺矿物学常规研究,查明主要载金矿石是黄铁矿、黄铜矿等硫化矿及石英、长石等脉石矿物;金多数以裂隙金方式存在,少数以包裹金、晶隙金存在。随着井下纵深采矿的进一步开拓,现有工艺矿物学结论明显难以满足后续选冶技术的调控要求。本文采用偏光反光两用显微镜、扫描电镜、X射线衍射仪及能谱分析等先进手段,综合分析了焦家金矿主矿区的金矿石样品中金的不同赋存状态、嵌布粒度特征,金的形状、种类、分布状态,不同粒度金的成色特征等。研究结果显示主要载金矿物是黄铁矿、黄铜矿、闪锌矿、方铅矿、石英和长石;金赋存状态有裂隙金、包裹金和晶隙金;金矿物形状有球形、三角形等,其中球形颗粒最多;金矿物粒度范围较大,大金颗粒可达到90~110 μm,小金颗粒只有2~3 μm;金矿物种类丰富,有自然金、银金矿、自然银、金铜矿等;金银矿物的成色普遍较高,赋存于硫化物中的大颗粒金,成色较低,小粒金成色较高。

    本文利用不同先进手段表征互相印证,提高了研究结论的可靠性与准确性,丰富了焦家金矿内生热液金矿成矿区含金黄铁矿石英脉型金矿工艺矿物学的研究内容,为后续选冶作业提供了调控依据和重要矿物信息。由于缺乏对主要含金脉石矿物(石英、长石等)构造特征,金在其中赋存状态、嵌布关系等方面的研究,应进一步对其进行系统分析,以为提高后续选矿作业脉石中金的回收效果提供重要的分析依据。

  • 图  1   天然气水合物热物理实验装置图

    Figure  1.   Schematic diagram of thermal physical measurement apparatus of gas hydrate

    图  2   琼脂-水胶体温度变化(△θ)与时间(lnt)关系图

    Figure  2.   The relationship between △θ and lnt of agar solution

    图  3   饱和水沉积物结冰过程中TDR波形的变化

    Figure  3.   TDR waveform of ice freezing process in water saturated sediment

    图  4   神狐海域沉积物样品的粒度分布

    Figure  4.   Particle size distribution of Shenhu marine sediment

    图  5   沉积物中水合物饱和度的变化趋势

    Figure  5.   Saturation trend of hydrate in sediment

    图  6   热脉冲过程中温度-时间关系图

    Figure  6.   The temperature-time curve during thermal pulse

    图  7   不同反应时间热脉冲△θ-lnt关系图

    Figure  7.   The relationship between △θ and lnt

    图  8   沉积物热导率-水合物饱和度关系图

    Figure  8.   The relationship between sediment thermal conductivity and hydrate saturation

  • 许红,黄君权,夏斌,蔡乾忠.最新国际天然气水合物研究现状与资源潜力评估(上)[J].天然气工业,2005, 25(5): 21-25. http://www.cnki.com.cn/Article/CJFDTOTAL-TRQG200505006.htm
    许红,黄君权,夏斌,蔡乾忠.最新国际天然气水合物研究现状及资源潜力评估(下)[J].天然气工业,2005, 25(6): 18-23. http://www.cnki.com.cn/Article/CJFDTOTAL-TRQG200506005.htm

    Sloan E D, Bloys J B, Hydrate Engineering[M]. USA: Society of Petroleum Engineers, 2000: 148-152.

    Winters W J, Pecher I A, Waite W F, Mason D H. Physical properties and rock physics models of sediment containing natural and laboratory-formed methane gas hydrate[J]. American Mineralogist, 2004, 89(8-9): 1221-1231. doi: 10.2138/am-2004-8-909

    Winters W J, Waite W F, Mason D H, Gilbert L Y, Pecher I A. Methane gas hydrate effect on sediment acoustic and strength properties[J]. Journal of Petroleum Science and Engineering,2007, 56(1-3): 127-135. doi: 10.1016/j.petrol.2006.02.003

    金春爽,汪集旸.天然气水合物的地热研究进展[J].地球科学进展,2001, 16(4): 540-543. http://www.cnki.com.cn/Article/CJFDTOTAL-DXJZ200104015.htm

    Stoll R D, Bryan G M. Physical properties of sediments containing gas hydrates[J]. Journal of Geophysical Research, 1979, 84(B4): 1629-1634. doi: 10.1029/JB084iB04p01629

    Waite W F. Thermal conductivity measurements in porous mixtures of methane hydrate and quartz sand[J]. Geophysical Research Letters, 2002, 29(24): 2229-2235.

    Rosenbaum E J, Niall J, Johnson J K, Shaw D W, Warzinski R P. Thermal conductivity of methane hydrate from experiment and molecular simulation[J]. The Journal of Physical Chemistry B,2007, 111(46): 13194-13205. doi: 10.1021/jp074419o

    刁少波,业渝光,岳英杰,张剑,陈强,胡高伟.多孔介质中水合物的热物理参数测量[J].岩矿测试,2008, 27(3): 165-168. http://www.cnki.com.cn/Article/CJFDTOTAL-YKCS200803006.htm
    任图生,邵明安,巨兆强, Horton R.利用热脉冲-时域反射技术测定土壤水热动态和物理参数Ⅰ.原理[J].土壤学报,2004, 41(2): 225-229. doi: 10.11766/trxb200305160210
    任图生,邵明安,巨兆强, Horton R.利用热脉冲-时域反射技术测定土壤水热动态和物理参数Ⅱ.应用[J]. 土壤学报,2004, 41(4): 523-529. doi: 10.11766/trxb200305160405

    Campbell G, Calissendorff C, Williams J. Probe for measuring soil specific heat using a heat-pulse method[J]. Soil Science Society of America Journal,1991, 55(1): 291-293. doi: 10.2136/sssaj1991.03615995005500010052x

    Tarara J M, Ham J M. Measuring soil water content in the laboratory and field with dual-probe heat-capacity sensors[J]. Agronomy Journal,1997, 89(4): 535-542. doi: 10.2134/agronj1997.00021962008900040001x

    Wright J, Nixon F, Dallimore S, Matsubayashi O.A method for direct measurement of gas hydrate amounts based on the bulk dielectric properties of laboratory test media [C]//The Fourth International Conference on Gas Hydrate, ICGH-Ⅳ. Japan: Yokohama, 2002: 745-749.

    胡高伟,业渝光,刁少波,张剑,刘昌岭.时域反射技术测量海洋沉积物含水量的研究[J].现代地质,2010, 24(3): 622-626. http://www.cnki.com.cn/Article/CJFDTOTAL-XDDZ201003030.htm
    陈强,业渝光,刘昌岭.多孔介质体系中甲烷水合物生成动力学的模拟实验[J].海洋地质与第四纪地质, 2007,27(1): 111-116. http://www.cnki.com.cn/Article/CJFDTOTAL-HYDZ200701019.htm

    Dalmazzone D, Hamed N, Dalmazzone C. DSC measurements and modelling of the kinetics of methane hydrate formation in water-in-oil emulsion[J]. Chemical Engineering Science, 2009, 64(9): 2020-2026. doi: 10.1016/j.ces.2009.01.028

    胡高伟.南海沉积物的水合物声学特性模拟实验研究[D].武汉:中国地质大学, 2010.
图(8)
计量
  • 文章访问数:  1100
  • HTML全文浏览量:  285
  • PDF下载量:  6
  • 被引次数: 0
出版历程
  • 收稿日期:  2012-03-26
  • 录用日期:  2012-08-05
  • 发布日期:  2013-01-31

目录

/

返回文章
返回