资源描述:
文章编号1001-6112201201-0019-06 川川东北元坝地区中生代构造与动态热演化史 磷灰石、锆石U-Th / He 定年分析 王摇 杰,秦建中,刘文汇,陶摇 成,腾格尔 1. 中国石油化工集团公司 油气成藏重点实验室,江苏 无锡摇 214151; 2. 中国石油化工股份有限公司 石油勘探开发研究院 无锡石油地质研究所,江苏 无锡摇 214126 摘要通过对川东北元坝地区须家河组T3x嘉定组K1j钻井岩屑样品镜质体反射率和锆石、磷灰石U-Th / He 定年分析, 建立了该区 He 年龄深度/ 温度动态演化模式,推断出元坝地区磷灰石 He 封闭温度为 95 益左右。 元坝地区 T3x-K1j 中生代地 层基本都经历了磷灰石 He 封闭温度95 益;所有样品未经历锆石 He 封闭温度,T3x2-J1z 地层部分样品可能经历了约 170 益的 最高古地温。 元坝地区中生代地层在古近纪新近纪0.2 36. 4 Ma发生重大冷却抬升剥蚀,剥蚀速率约为 109. 9 m/ Ma,K1j 及以上地层最大抬升剥蚀厚度约为 4 000 m。 系统揭示了该区动态热演化历史,中生代地层最高古地温接近于 170 190 益,随 后地层发生抬升,古地温下降;在 36 176 Ma 之间时,古地温在 95 170 益之间;在 0 36 Ma 时,现今地温小于 95 益。 关键词镜质体反射率;磷灰石和锆石 U-Th/ He 定年;剥蚀厚度;热演化史;元坝地区;川东北 中图分类号TE121.1摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 文献标识码A Mesozoic tectonics and dynamic thermal history in Yuanba area of northeastern Sichuan Basin application of U-Th / He dating of apatite and zircon Wang Jie, Qin Jianzhong, Liu Wenhui, Tao Cheng, Tenger 1. SINOPEC Key Laboratory of Petroleum Accumulation Mechanisms, Wuxi, Jiangsu 214151, China; 2. Wuxi Institute of Petroleum Geology, SINOPEC, Wuxi, Jiangsu 214151,China Abstract In Yuanba area of the northeastern Sichuan Basin, vitrinite reflectance and U-Th / He dating analy鄄 ses of apatite and zircon in chip samples from the Xujiahe T3x-Jiading K1j ations were carried out. A dynamic evolution model based on He age and depth/ temperature was set up, by means of which we concluded that the He closure temperature of apatite was about 95 益. Most of the Mesozoic strata T3x-K1j in the study area experienced the He closure temperature of apatite 95 益 while no sample experienced the He closure tem鄄 perature of zircon, and some samples from T3x2-J1z underwent the maxium palaeogeothermal of about 170 益. From Paleogene to Neogene 0. 2-36. 4 Ma, the Mesozoic strata in Yuanba area went through a severe cooling uplift and denudation with the denudation rate of 109. 9 m/ Ma. In K1j and the overlying strata, the maximum de鄄 nudation thickness was about 4 000 m. The dynamic thermal evolution history in the study area was illustrated systematically as followed. In Mesozoic, the maximum palaeogeothermal of strata was 170-190 益. As strata up鄄 raise, palaeogeothermal decreased. From 36 to 176 Ma, palaeogeothermal was 95-170 益. From 0 to 36 Ma, the present geotemperature is below 95 益. Key words vitrinite reflectance; U-Th/ He dating of apatite and zircon; denudation thickness; thermal evolu鄄 tion history; Yuanba area; northeastern Sihuan Basin 摇 摇 国内外对沉积盆地热演化历史的恢复方法大 体上有 3 种1正演技术,利用盆地演化的热动力 学模型恢复热史[1-2];2反演技术,利用各种古温 标来恢复古地温[3-6];3将正演技术与反演技术 结合起来的综合法[7]。 古温标方法被认为是研究 精度较高且切实可行的方法,前人已进行了一些有 益的探索,取得了重要进展[8-12]。 近年来,逐渐发 展起来的U-Th / He 低温热定年技术为盆地热历 收稿日期2011-01-18;修订日期2011-12-29。 作者简介王杰1975,男,博士,高级工程师,从事油气地质与有机地球化学研究工作。 E-mailwangjie. syky sinopec. com。 基金项目国家重点基础研究973 项目2005CB422102和国家重大专项2008ZX05005-01-02。 摇 摇 第 34 卷第 1 期 2012 年 1 月 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 石摇 油摇 实摇 验摇 地摇 质 PETROLEUM GEOLOGY 该区位于 上扬子地块北缘,是由古生代和中生代地层组成的 一个复杂变形系统。 自晋宁期基底形成以后,该区 在震旦纪中三叠世演化为被动大陆边缘盆地;晚 三叠世演化为前陆盆地;侏罗纪早白垩世发展成 陆内坳陷盆地;晚白垩世沉积盆地萎缩衰亡,开始进 入以抬升剥蚀为主的构造演化阶段[19]。 研究区震 旦系至中三叠统为海相地层,发育碳酸盐岩夹碎屑 岩沉积;上三叠统至下白垩统为陆相地层,发育碎屑 岩沉积。 前人对川东北地区热历史做过部分研究工 作[20-23],但对于地温随时间演化历史的研究则不 多见。 目前元坝地区是中国石化在川东北地区的 重点勘探区域,研究该区构造及动态热演化历史, 对天然气资源潜力评价及勘探具有重要意义。 1摇 样品与实验 获取合适的磷灰石和锆石样品是得到准确 He 年龄 数 据 的 关 键。 由 于 元 坝 地 区 须 家 河 组 T3x嘉定组K1j为碎屑岩沉积,砂岩或含砂 泥岩发育,可以获取合适的磷灰石和锆石颗粒进行 U-Th / He 年龄测试。 本研究样品取自元坝 3 井 T3x-K1j 钻井岩屑,岩屑深度数据根据钻井地质分 层数据获得,挑选出合适的磷灰石、锆石颗粒矿物; U-Th / He 年龄测试在澳大利亚联邦科学与工业 研究院CSIRO进行,样品 He 年龄的测试过程详 见文献[13],所有样品 He 年龄最后均进行 琢 发射 效应校正方法见 Farley 等,1996[24]。 本次测试 He 年龄数据见表 1。 2摇 构造与动态热演化史反演 由于不同矿物颗粒的U-Th / He 定年体系封 闭温度差异较大,所以根据不同矿物颗粒U-Th / He 体系的 He 年龄差、封闭温度差以及地温梯度 图 1,恢复地层的剥蚀速率和剥蚀厚度,进而可 以反演研究区时间温度的动态热演化历史[13]。 2. 1摇 U-Th / He 封闭温度和年龄 元 坝地区T3x2-K1j砂岩样品磷灰石He年龄 表 1摇 川东北元坝 3 井磷灰石和锆石U-Th / He 年龄 Table 1摇 U-Th/ He age of apatite and zircon from well Yuanba 3, northeastern Sichuan Basin 层位深度/ m海拔/ m 磷灰石 He 年龄/ Ma 单颗粒平均值 锆石 He 年龄/ Ma 锆石经历 温度 现今 地温/ 益 K1j87390. 3 J1z4 154.5-3 677. 2 T3x44 372-3 894. 7 T3x44 476.5-3 999. 2 T3x24 776.5-4 299. 2 35. 7 13. 1 60. 5 14. 9 20. 9 17. 8 5.9 0.2 1.2 2.8 1.3 2.4 0.2 0.1 0.2 摇36. 4 Ma,接 近于封闭温度 小于 100 益 摇14. 9 Ma,经 历封闭温度大 于 100 益 摇 1. 6 Ma,经历 封闭温度大于 100 益 摇0. 17 Ma,经 历封闭温度大 于 100 益 157.5 178.6 194.5 135.2 194.1 168.1 127.0 111.9 220.2 97.3 112.0 182.1 149.1 222.2 181.0 163.6 202.0 摇未 经 历 锆 石 封 闭 温度 25 摇接 近 于 锆 石 封 闭 温度 95 摇接 近 于 锆 石 封 闭 温度 98 摇接 近 于 锆 石 封 闭 温度 100 摇接 近 于 锆 石 封 闭 温度 105 摇 摇 摇 摇 摇 摇 摇 摇 摇注实测温度数据依据地温梯度得到,取 1.7 益 / hm,地表温度均定为 24 益。 02摇 摇 摇 摇 摇 摇 石摇 油摇 实摇 验摇 地摇 质摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇第 34 卷摇 摇 图 1摇 磷灰石和锆石 U-Th/ He 定年体系及定年时钟效应 Fig. 1摇 U-Th/ He dating system and clock effect of apatite and zircon 均远小于样品沉积年龄,表明该区上述层位样品都 进入了磷灰石 He 封闭温度和封闭体系。 元坝地区 磷灰石 He 年龄随现今埋藏深度或实测地层温度增 大而逐渐减小图2,表现出地层沉积年龄越小,磷 灰石 He 年龄越大的特征。 He 年龄随温度/ 深度演 化剖面反映出该区磷灰石U-Th/ He 体系封闭温 度约为95 益左右,此时埋藏深度约为4 500 m。 元坝地区 K1j 样品磷灰石 He 年龄相对分散, 分布在13. 1 60. 5 Ma 之间,处于 He 部分保留区, 磷灰石中 He 部分热重置,说明该样品比较接近于 封闭温度和封闭体系,该地层沉积时其最高古地温 曾接近磷灰石封闭温度,其它样品都经历了磷灰石 封闭温度。 由表 1 和图 2 可见,元坝地区 T3x2-K1j 钻井砂岩样品磷灰石 He 年龄介于 0. 1 60. 5 Ma, 均值为 0. 2 36. 4 Ma,说明该区中生代地层在古 近纪新近纪0 60 Ma 发生主要抬升剥蚀事 件,比中生代普光地区抬升剥蚀时间略微要晚普 光地区 1. 5 25. 9 Ma,比通南巴地区抬升剥蚀时 间要早通南巴地区 9. 2 25. 8 Ma [13],古地温逐 渐降低到磷灰石 He 封闭温度,直至现今地温。 元坝地区 K1j 样品锆石U-Th / He 年龄为 157. 5 Ma 和 178. 6 Ma,均大于样品地层沉积年龄 表 1,图 2,反映了该样品曾经经历的最高古地 温不超过锆石 He 封闭温度。 该样品的锆石颗粒 没有经历盆地最后一期热事件的重置或重置不完 全,说明元坝地区 K1j 地层锆石颗粒来源于在侏罗 纪经历了锆石U-Th / He 封闭温度的抬升剥蚀 区,表明该区在侏罗纪发生抬升剥蚀冷却。 T3x2- J1z 部分样品锆石 He 年龄大于地层沉积年龄,且分 布相对分散,表明上述样品可能比较接近于锆石封 闭温度,曾经历了接近于 170 益 左右的最高古地 温,该区埋藏史图印证了这一点。 由图2 可见,随现今地温或埋藏深度的增加,磷 灰石He 年龄及变化范围逐渐变小,说明该地区磷灰 石经历了封闭温度,He 年龄反映了地层抬升剥蚀的 主要时期,该区在古近纪新近纪0 60 Ma遭受 强烈抬升剥蚀事件;而锆石He年龄及变化范围随 图 2摇 川东北元坝 3 井磷灰石和锆石U-Th / He 年龄随埋藏深度和现今地温的关系 Fig.2摇 U-Th/ He age of apatite and zircon vs. burial depth and present geotemperature, well Yuanba 3, northeastern Sichuan Basin 12摇 第 1 期摇 摇 摇 王杰,等. 川东北元坝地区中生代构造与动态热演化史 磷灰石、锆石U-Th / He 定年分析 现今地温或埋藏深度的增加有变大的趋势,表明该 区中生代地层没有经历锆石封闭温度,只是处于 He 部分保留区域内,其经历的最高古地温曾接近 于锆石封闭温度170 190 益,所以样品锆石 He 年龄可以在一定程度上反映样品沉积时间。 2. 2摇 剥蚀历史恢复 从磷灰石的 He 年龄数据可以得到元坝地区 在古近纪新近纪0. 1 60. 5 Ma处于冷却抬升 剥蚀的主要时期。 假使该区在地质历史时期古地 温梯度恒定1. 7 益 / hm,从磷灰石 He 年龄与埋 藏深度的关系图 2来看,利用年龄高程法[25]推 算出此期间的平均剥蚀速率约为 109. 9 m/ Ma,大 致相当于 1. 9 益 / Ma 平均降温速率。 鉴于 K1j 样 品磷灰石处于 He 部分保留区,可能接近于封闭温 度,其最高古地温曾接近磷灰石封闭温度,K1j 及 以上地层最大抬升剥蚀厚度约为 4 000 m 左右。 根据锆石 He 年龄数据可以得到该区在侏罗世期 间处于抬升剥蚀冷却的状态。 元坝 3 井实测 RoRb值随埋藏深度增加呈半 对数直线变化的方程为D 4 305. 8 lnVRo或 Rb1 655. 3R2 0. 89。 式中D 为埋藏深度, VRo或 Rb为镜质体反射率或沥青反射率。 其直线 斜率较高,说明古地温梯度相对较低图 3。 当接 近地表 Ro值约为 0. 25 左右时,深度约为-4 300 m,推算该井最大剥蚀厚度约为 4 300 m 左右该 井未有火成岩活动,与根据磷灰石U-Th / He 年 龄与封闭温度得出的剥蚀深度基本上是一致的。 综合来说,该区 K1j 及以上地层抬升剥蚀厚度约为 4 000 m 左右。 2. 3摇 热史动态演化历程 由表 1、图 2 和图 4 可见,川东北元坝地区 T3x2-K1j 样品磷灰石 He 年龄随埋藏深度或实测 地层温度增大逐渐减小,说明该区样品都经历了磷 灰石 He 封闭温度。 该区埋藏深度达到 4 800 m 左 右时,样品锆石 He 年龄小于地层年龄,可能接近 于锆石封闭温度。 所以,元坝地区深度大于 4 800 m 的样品最高古地温接近于 170 190 益;随后地 层抬升,古地温下降,当样品在 35. 7 194 Ma 之间 时,古地温在 95 170 益之间;当样品在 0 35. 7 Ma 时,进入磷灰石封闭温度,现今地温小于95 益。 根据不同层位样品进入磷灰石封闭温度的时间对 埋藏史图进行限制,使之更接近于实际地质条件。 3摇 结果与讨论 1元坝地区磷灰石 He 年龄随样品沉积年龄 变小而越大,反映了典型的磷灰石经历封闭温度后 的 He 年龄变化特征,建立磷灰石 He 年龄深度/ 温度动态演化模式,该区磷灰石 He 年龄的封闭温 度为 95 益。 2元坝地区中生代地层在古近纪新近纪 0. 2 36. 4 Ma发生重大冷却抬升剥蚀事件,剥 蚀速率约为109. 9 m/ Ma,K1j 及以上地层最大剥蚀 厚度为 4 000 m 左右。 3元坝地区 T3x-K1j 中生代地层基本上都经 历了磷灰石 He 封闭温度95 益,K1j 地层接近于 封闭温度,其沉积时最高古地温曾接近磷灰石封闭 图 3摇 川东北元坝 3 井镜质体反射率与磷灰石、锆石U-Th / He 年龄以及现今埋藏深度的关系 Fig.3摇 Vitrinite reflectance, U-Th/ He age of apatite and zircon vs. present burial depth, well Yuanba 3, northeastern Sichuan Basin 22摇 摇 摇 摇 摇 摇 石摇 油摇 实摇 验摇 地摇 质摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇第 34 卷摇 摇 图 4摇 川东北元坝 3 井地层埋藏史和热演化史 Fig. 4摇 Burial and thermal evolution histories in well Yuanba 3, northeastern Sichuan Basin 温度;所有样品未经历锆石 He 封闭温度,T3x2-J1z 地层部分样品可能曾经历了 170 益左右的最高古 地温。 根据磷灰石、锆石封闭温度和 He 年龄,揭 示了该区动态热演化历史,该区中生代地层最高古 地温接近于 170 190 益,随后地层发生抬升,古地 温下降;在 36 176 Ma 之间时,古地温在 95 170 益之间;在 0 36 Ma 时,现今地温小于 95 益。 参考文献 [1]摇 何丽娟. 沉积盆地构造热演化模拟的研究进展[J]. 地球科学 进展,2000,156661-665. [2]摇 Mckenzie D. Some remarks on the development of sedimentary ba鄄 sins[J]. Earth and Planetary Science Letters,1978,4025-32. [3]摇 秦建中,饶丹,蒋宏. 高演化海相碳酸盐岩层系古温标的直接指 标包裹体均一温度[J].石油实验地质,2008,305494-498. [4]摇 袁玉松,郑和荣,涂伟. 沉积盆地剥蚀量恢复方法[J]. 石油实 验地质,2008,306636-642. [5]摇Qin Jianzhong, Wang Jie, Qiu Nansheng. Evidence of thermal evolution history of northeast Sichuan basinU-Th / He low temperature thermochronometry of apatite and zircon[J]. Journal of China University of Geosciences,2008,196591-601. [6]摇徐二社,邱楠生,秦建中,等. 沉积盆地有机质自由基热演化 特征及其作为古温标的探索[J]. 地质学报,2008,283 413-419. [7]摇 任战利,刘丽,崔军平. 盆地构造热演化史在油气成藏期次研 究中的应用[J]. 石油与天然气地质,2008,294502-506. [8]摇 曹长群,尚庆华,方一亭. 探讨笔石反射率对奥陶系、志留系烃源 岩成熟度的指示作用[J].古生物学报,2000,391151-156. [9]摇 邱楠生,李慧莉,金之钧. 碳酸盐岩层系热历史恢复的有机质 自由基古温标研究[J]. 地质学报,2006,803390-397. [10]摇 何谋春,吕新彪,刘艳荣.激光拉曼光谱在油气勘探中的应用研 究初探[J]. 光谱学与光谱分析,2004,24111363-1366. [11]摇 Orchard M J A. Conodont based temporal and spatial framework for the Triassic of North America[J]. AAPG Bulletin,2002, 8661156. [12]摇 王飞宇,刘德汉. 镜状体反射率可作为下古生界高过成熟烃 源岩成熟度标尺[J]. 天然气工业,1996,16414-19. [13]摇 秦建中,王杰,邱楠生,等. 反演南方海相层系热史动态演化 的新温标 磷灰石、锆石U-Th / He 年龄和封闭温度[J]. 石 油与天然气地质,2010,313277-287. [14]摇 Evans N J,Wilson N S F,Cline J S,et al. Fluorite U-Th / He thermochronologyConstraints on the low temperature history of Yucca Mountain,Nevada[J]. Applied Geochemistry,2005,20 1099-1105. [15]摇 Wolf R A,Farley K A,Silver L A. Helium diffusion and low temperature thermochronometry of apatite [ J]. Geochimica et Cosmochimica Acta,1996,60214231-4240. [16]摇Reiners P W,Farley K A,Hickes H J. Helium diffusion and 32摇 第 1 期摇 摇 摇 王杰,等. 川东北元坝地区中生代构造与动态热演化史 磷灰石、锆石U-Th / He 定年分析 U-Th / He thermochronometry of zircon Initial results from Fish Canyon Tuff and Gold Butte[J]. Tectonophysics,2002, 349247-308. [17]摇Reiners P W,Farley K A. Helium diffusion and U-Th / He thermochronometry titanite [ J]. Geochimica et Cosmochimica Acta,1999,63223845-3859. [18]摇 Wolf R,Farley K A,Kass D. A sensitivity analysis of the apatite U-Th / He thermochronometer[J]. Chemical Geology,1998, 1481-2105-114. [19]摇 郭正吾,邓康龄,韩永辉,等. 四川盆地形成与演化[M]. 北 京地质出版社,199648-88. [20]摇 韩永辉,吴春生. 四川盆地地温梯度及几个深井的热流值[J]. 石油与天然气地质,1993,14180-84. [21]摇 王一刚,余晓锋,杨雨,等. 流体包裹体在建立四川盆地古地温 剖面研究中的应用[J].地球科学,1998,233285-288. [22]摇 卢庆治,胡圣标,郭彤楼,等. 川东北地区异常高压形成的地 温场背景[J]. 地球物理学报,2005,4851110-1116. [23]摇 沈传波,梅廉夫,郭彤楼. 川东北地区中、新生代热历史的裂 变径迹分析[J]. 天然气工业,2007,27724-26. [24]摇 Farley K A,Wolf R A,Silver L T. The effects of long alpha stop鄄 ping distances on U-Th / He ages[J]. Geochimica et Cosmo鄄 chimica Acta,1996,60214223-4229. [25]摇 Gleadow A J W,Fitzgerald P G. Uplift history and structure of the Transantarctic mountainsNew evidence from fission track dating of basement apatite in the dry valleys area, southern Victoria land[J]. Earth and Planetary Science Letters,1987,821-14. 编辑摇徐文明 势势势势势势势势势势势势势势势势势势势势势势势势势势势势势势势势势势势势势势势势势势势势势势 上接第 18 页 [17]摇 Nunn J A,Meulbroek P. Kilometer-scale upward migration of hydrocarbons in geopressured sediments by buoyancy - driven propagation of methane - flied fractures [ J]. AAPG Bulletin, 2002,865907-918. [18]摇 潘钟祥. 不整合对于油气运移聚集的重要性及寻找不整合 面下的某些油气藏[J]. 地质论评,1983,294374-381. [19]摇李坤,沈忠民,肖宗林. 卡塔克古隆起油气运移输导体系演 化特征[J]. 断块油气田,2011,183317-320. [20]摇 Bethke C M,Reed J D,Oltz D F. Long-range petroleum migration in the Illinois basin[J].AAPG Bulletin,1991,755925-945. [21]摇 吕修样.塔里木盆地不整合面石油运移散失量的实验研究[J]. 石油大学学报自然科学版,2000,244112-114. [22]摇 李朝阳,刘玉平,管太阳,等. 不整合面中的成矿机制与找矿 研究[J]. 地学前缘,2004,112353- 360. [23]摇 Nie F J,Li S T,Wang H,et al.Lateral migration pathways of petro鄄 leum in the Zhu III subbasin, Pearl river mouth basin, South Chi鄄 na Sea[J]. Marine and Petroleum Geology,2001,18561-575. [24]摇 杨宁,吕修祥,陈梅涛,等. 英买力地区碳酸盐岩古风化壳油气 聚集模式探讨[J]. 天然气地球科学,2005,32135-39. [25]摇 Xie Xinong,Li Sitian,Hu Xiangyun,et al. Conduit systems and ation mechanism of heat fluids in diapiric belt of Ynggehai basin, China[J]. Science in ChinaSeries D,1999,426 561-571. [26]摇王振峰,何家雄,解习农. 莺歌海盆地泥一流体底辟带热流 体活动对天然气运聚成藏的控制作用[J]. 地球科学 中 国地质大学学报,2004,292203-210. [27]摇 沈朴,张善文,林会喜,等. 油气输导体系研究综述[J]. 油气 地质与采收率,2010,1744-8. [28]摇 Hindle A D. Petroleum migration pathway and charge concentra鄄 tionA three - imensional model [ J]. AAPG Bulletin,1997, 8191451-1481. [29]摇 谢泰俊,潘祖荫,杨学昌. 油气运移动力及通道体系[M] / / 龚再升,李思田. 南海北部大陆边缘盆地分析与油气聚集. 北京科学出版社,1997385-405. [30]摇顾乔元,韩剑发,张耀堂. 轮南地区奥陶系断裂及储集层裂 缝分形特征[J]. 新疆石油地质,1998,194287-290. [31]摇 陈强路,黄继文. 塔河油田油气输导体系与运聚模式[J]. 新 疆石油地质,2004,252143-146. [32]摇 吕修祥,杨宁,解启来,等. 塔中地区深部流体对碳酸盐岩 储层的改造作用[J]. 石油与天然气地质,2005,263 284-289. [33]摇 金之钧. 塔里木盆地热液活动地质地球化学特征及其对储 层影响[J]. 地质学报,2006,802245-253. [34]摇杨宁,吕修祥,陈梅涛. 塔里木盆地塔河油田奥陶系碳酸盐 岩油气成藏特征[J]. 西安石油大学学报自然科学版, 2008,2331-6. [35]摇顾忆,黄继文,邵志兵. 塔河油田奥陶系油气地球化学特征 与油气运移[J]. 石油实验地质,2003,256746-750. [36]摇 饶丹,马绪杰,贾存善,等. 塔河油田主体区奥陶系缝洞系统 与油气分布[J]. 石油实验地质,2007,296589-592. 编辑摇徐文明 42摇 摇 摇 摇 摇 摇 石摇 油摇 实摇 验摇 地摇 质摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇第 34 卷摇 摇
展开阅读全文