用户名: 密码: 验证码:
松科1井复杂地层取心钻具及泥浆优化设计和应用研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
随着钻探技术的发展,无论是岩心钻探、科学钻探、还是工程勘查皆对岩心采取率及岩心质量提出了更高的要求,因此,现在的取心质量包括岩心采取率和岩心的原状结构程度。本文合松科1井对极松散地层与强塑性松软地层进行了隔液保真取样钻具及防泥包钻头的优化设计,并研究了复杂地层防塌水基泥浆体系。
     本文分析了在极为复杂地层中钻进岩心采取率极低的原因,得出在极松散地层钻进,钻具内腔的液压力需恒大于或等于内、外腔连通处的液压力;在强塑性松软地层中钻进,依靠布齿形式、排屑空间与水路结构的合理设计与钻进参数的控制防止钻头泥包;为此,开辟了复杂地层取心钻具优化设计的新思路。
     从取心钻具与水力学特性相互作用的逻辑关系入手,进行取心钻具的结构强度分析,确定了取心钻具结构空间的设计依据;然后结合水力学进行了CFD模拟;泄流槽和钻头内的挡水环是该套钻具的核心技术,能够满足保心效果;仿真分析与模拟试验表明,参数吻合性好;通过对防泥包钻头的CFD模型的分析,得出孔底流场分布均匀、效果良好,具有一定的科学性。
     优化了螺纹锥度、牙高、螺距、基面距等参数,螺纹强度受钻杆局部内径变化影响显著,内径减小,螺纹整体强度提高;螺纹锥度变化对钻杆强度的影响较大,当锥度达到1:25附近时,螺纹强度最大;螺纹强度随牙高变化极其复杂。随着螺距的增大,钻杆管体螺纹处结构内部应力分布愈不均衡,从而使螺纹的应力集中现象更加显著;螺纹强度随基面距的变化影响不大。
     结合计算流体力学的基本知识及经验,完成了隔液保真钻具和防泥包钻头的设计和模拟,掌握了其孔底流场分布规律,通过大量的实验测试:结果显示与模拟结果相符。
     隔液保真钻具的优化结果:仿真分析与CFD模拟试验表明,参数吻合性好,但泄流槽的深度对井底流域的压力场和速度场影响较为明显。50mm槽深井底流速场变化要大,对井底扰动明显,压力梯度较为明显,而70mm深泄流槽的孔底流动相对平缓,可以预测保心效果好,但估计会影响钻速。防泥包钻头优化结果:钻头切削刃前斜面的倾角对孔底流场的分布影响很大,但对孔底压力损失影响不大。钻头在孔底总的压降,30°斜面最高,更容易取心与排屑,25°斜面次之,35°斜面最小。
     以优化设计的结构参数为依据,分别针对松科一井、内蒙查干温都尔矿区、汶川科学钻井使用的取心钻具进行了优化设计,制造出隔液保真钻具、防泥包钻具;结合工程实际优选泥浆配方,设计了全新测试实验系统,开发了新型简易仪器,测试与观测泥浆的性能,研制松科一井防塌水基钻探泥浆体系。实践证明:两种钻具明显提高了复杂地层中钻进岩心采取率的提高,满足了地质的要求;防塌水基钻探泥浆在松科一井中应用成功,起到了很好的护壁效果;验证了理论设计与优化方案的可行性。
With the development of exploration technology, core drilling, scientific exploration engineering and engineering investigation all improve the standards of recovery factor and core quality. This paper, in the light of the ultimate unconsolidated formation and the severely plastic soft formation, conducts optimized design of spacing fluid shape-retained drilling tool sampler and balling-control bit. It also makes researches on the anti-caving water-based drilling mud system in the complex formation based on the Well-1 of Songliao Scientific Drilling. spacing fluid
     This paper analyzes the cause of the extremely low recovery factor in the very complex formation, and concludes that when drilling in the ultimate unconsolidated formation, the fluid pressure within the drilling tools should be identically higher than or equal to that of the joint;when drilling in the severely plastic soft formation, the cutter’s distribution,cutting discharge space and the waterway structure should be rationally arranged, and at the same time, drilling parameters should be under control, with the view of preventing bit being balled. Therefore, a new method on the optimized design of core drillings in the complex stratum is put forward.
     Starting with the logic relationship of the interaction between coring drilling tools and hydraulics, the author analyzes the structure intensity of coring drilling tools, decides the design principles of their structural space, and then conducts CFD simulation based on hydraulics. The sillway and the water ring in the rig are the core of this set of drilling tools, which guarantee the core-protection effect. The simulation analysis and the CFD simulation experiment demonstrate that the parameters have good uniformity. The analysis of balling-conrol bit’CFD mould shows that the down-the–hole velocity field distributes evenly with good effect, which prove its scientificness.
     Through the optimization of structure parameters, such as the thread taper of drill tools, the tooth height, the pitch of screws and the datum distance, etc. , the author has the findings as follows. The thread strength is greatly influenced by the change of rod’s partial I.D., i.e., when the I.D. reduces, the whole thread strength rises; the rod strength is greatly influenced by the thread taper, i.e., when the taper is near 1:25, the thread strength is at peak; the relation between the thread strength and the tooth height is very complex. As the pitch of screws increases, the inside pressure on screw thread within the rod is unbalanced, thus, the pressure becomes more concentrated; the thread strength is less influenced by the datum distance.
     With the fundamental principle and experience of Computational Fluid Dynamics, the author finishes the design and the simulation of spacing fluid shape-retained drilling tool and balling-control bit, gets hold of the distribution pattern of down-the-hole drilling fluid field, and proves that the result is the same as simulation after numerous experiments and tests.
     Simulation analysis and the CFD simulation experiment both indicate that parameters have good uniformity, but the spillway depth affects the shaft bottom’s pressure field and velocity field a lot. The 50mm shaft bottom’s velocity field varies a lot, then the shaft bottom would be easy to be affected and its pressure gradient would be also relatively evident; the 70mm spillway has relatively slow velocity, so its core-protecting effect is predicted to be better, but it may affect the drilling rate. Above findings are got from the optimized result of the spacing fluid shape-retained drilling tool. And the optimized result of the balling-control bit demonstrates that the front bevel angel of bit-cutting-edge affects the velocity field structure a lot, but it has little effect on the down-the-hole pressure. The pressure drop of down-the-hole bit becomes the highest when the angle is 30°,which is the best moment to discharge cuttings. 25°bevel takes the second place, and then the 35°. According to the structure parameters of the optimized design, reinforcement theory and the optimized methods, the author designs the core drilling tools used in the Well-1 of Songliao Scientific Drilling., Inner Mongolia and We chuan satiric wells, manufactures the spacing fluid shape-retained drilling tool and the balling-control bit, and , on the basis of engineering practical optimized mud formula, develops the Well-1 of Songliao Scientific Drilling. anti-caving water-based drilling mud system. It has been found that the 2 kinds of drilling rods could highly increase the recovery factor in the relatively complex stratum, meeting the geological needs. The successful application of anti-caving water-based drilling mud, which acts as a good dado in the Well-1 of Songliao Scientific Drilling., proves the feasibility of the theoretical design and the optimized project.
引文
[美]陈景仁,湍流模型及有限分析法,上海交通大学出版社,1989年2月,第一版,12-45;
    [英]詹姆斯(James/C.)对比分析[M]青岛出版社2005-6-13 2-8.
    Aadnoy BS, Larsen K, Berg PC . Analysis of stuck pipe in deviated boreholes[J]. JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING, 2003,37(3-4):195-212.
    DESEVAUX P, PRENEL J P, HOSTACHE G. An optical analysis of an induced flow ejector using light polarization properties[J]. Experiments in Fluids,1994,16(3):165-170.
    Dogay S, Ozbayoglu E, Kok MV. Trajectory Estimation in Directional Drilling Using Bottom Hole Assembly Analysis[J]. ENERGY SOURCES PART A-RECOVERY UTILIZATION AND ENVIRONMENTAL EFFECTS, 2009,31(7):553-559.
    E.Tsaru,M.Orgasuwera. Fatigue of drillstring connections[J]. JOURNAL OF OFFSHORE MECHANICS AND ARCTIC ENGINEERING-TRANSACTIONS OF THE ASME,1995,117(2):126-132.
    EAMES I M, APHORNRATANA S, H IDER H. A theoretical and experimental study of a small-scale steam jet refrigerator[J]. Int. J . Refrig, 1995,18(6):18-23.
    Joern Lauterjung 2000. The International Continental Scientific Drilling Program(ICDP).ICDP Newsletter,(2):3-6.
    M.Guagliano. Fatigue failure of a bar of a twin-screw extruder for plastics[J]. Engineering Failure Analysis,2005,12( 6):128-131;
    Richard Selescu,A new type of shock-free axisymmetic supersonic flow,Nonlinear Analysis 47 2001.
    Ritto TG, Soize C, Sampaio R. Non-linear dynamics of a drill-string with uncertain model of the bit-rock interaction[J]. INTERNATIONAL JOURNAL OF NON-LINEAR MECHANICS, 2009,44(8): 865-876.
    Singh SP . The influence of geology on blasthole deviation[C]. PROCEEDINGS OF THE TWENTY-SECOND ANNUAL CONFERENCE ON EXPLOSIVES AND BLASTING TECHNIQUE, VOL II, 1996: 320-332.
    Taghavi R.,Raman G.,Flow characteristics of a rectangular multielement supersonic mixer-ejector,Journal of Propulsion and Power,v12 n5 Sep-0ct 1996.
    Vaz MA, Patel MH. Anslysis of drillstring in vertical and deviated holes using the galerkin technique[J]. ENGINEERING STRUCTURES, 1995,17(6): 437-442.
    Wang JL, Wang WZ. Research on the Impact Dynamic Simulation between Drillstring and Wellbore[C]. PROCEEDINGS OF THE SECOND INTERNATIONAL CONFERENCE ON MODELLING AND SIMULATION (ICMS2009), VOL 1,2009:262-266.Wolf S F.Field measurement of downhole drillstring vibrations[C]. SPE 14330,1985:237-242.
    艾春娟.复杂地层中深部绳索取心钻杆及钻具的改进研究: [学位论文].中国地质大学(北京), 2005.
    蔡记华,乌效鸣,杨倩云. LG植物胶的室内实验与初步应用研究.石油钻采工艺, 2006(12),28(6).
    蔡记华,乌效鸣等.LG植物胶处理剂的研究及应用《石油钻采工艺》2006-12-30
    蔡记华.暂堵型钻井液的试验研究《中国地质大学硕士论文》2003.05.01
    常学芳,薄玉成,徐健.基于CFDesign的导气孔流量系统的数值仿真[J],机械工程与自动化,2006,12(6):50-51.
    陈毓川.序(为全国深部找矿工程研讨会专辑)[J].中国地质, 2007, (增刊): 381-382.
    崔迎春,王贵和.钻井液技术发展趋势浅析.钻井液与完井液.2005(1),22(1),60~70.
    灯虎.裂缝性泥页岩的水化稳定性研究及其应用[博士学位论文].成都:西南石油学院.2004年4月.
    邓传光,PDC钻头泥包原因分析及对策[J,.钻采工艺, 2006, 29(6): 127-128.
    董智光,常喜顺.复杂地层金刚石钻进钻杆折断原因初探[J].煤田地质与勘探, 1994, (3): 33-36.
    段志东,刘克温.弹塑性变形对螺栓应力和轴力分布的影响[J].机械设计与制造, 2009, (4): 231-235.
    冯清文,唐胜利.国内外绳索取心钻进现状[J].探矿工程, 1990, (9): 57-60.
    冯清文.绳索取心钻杆螺纹锥度与强度关系及应力研究[J].探矿工程, 1995, (6): 24-27.
    甘行平.深部找矿中的关键技术问题及对策[J].中国地质, 2007, (增刊): 383-387.
    高长福,范斌,张庆昌等.钻杆失效分析与可靠性探讨[J].石油工业技术监督, 2007, (10): 47-49.
    高德利. 21世纪油气钻井技术展望[J].石油化工动态, 2000, 8(2): 31-34.
    高德利.谈谈石油钻柱失效问题[J].石油钻采工艺, 2000, (增刊): 9-12.
    高德利.油气钻井技术展望[J].石油大学学报(自然科学版), 2003, 27(2): 29-32.
    高磊.耐温耐盐滤失剂的合成与蒙脱土的相互作用研究《山东大学硕士论文》2010-05-01
    高有峰.松辽盆地陆相白垩系科学钻探—“松科1”井南孔岩石学特征研究,《吉林大学硕士论文》2007-05-20.
    耿瑞伦,陈星庆等.多工艺空气钻探[M].北京:地质出版社, 1998. 61-62.
    耿瑞伦,方啸虎.在我国危机矿山接替资源勘查工作中探矿工程新技术的应用[J].探矿工程(岩土钻掘工程), 2005, (增刊): 5-10.
    耿瑞伦,冯国强.液动锤技术及其应用前景[J].探矿工程(岩土钻掘工程), 2001, (1):32-33.
    耿瑞伦,张林霞,冉恒谦等. 21世纪地质钻探技术应用与发展趋势[J].探矿工程(岩土钻掘工程), 2001, (增刊): 9-13.
    谷芳.规整填料局部流动与传质的计算流体力学研究[D].天津:天津大学,2004:36-45.
    郝树青.贯通式潜孔锤反循环取心(取样)钻进孔底流场模拟与实验研究[D].长春:吉林大学,2007:18-22.
    何爱国等.进口大尺寸444.5mm钻头在两口科探井中的应用[J].石油钻采工艺, 1999, 21(5): 27-31.
    胡福昌.钻杆性能及其结构的发展[J].探矿工程, 1982, (5): 78-81.
    胡郁乐,张晓西,张恒春,等.CFD在保真取心钻具结构优化设计中的应用[J].地质与勘探,2009(9):626-628.
    黄纪达.金刚石钻探用绳索取心钻杆及接头的强化方法探讨[J].探矿工程, 1991, (1): 22-25;
    黄麟森.基于ANSYS的钻杆螺纹结构形式分析[J].煤矿机械, 2009, 30 (7): 3-5.
    姜光忍,李忠,王献斌.绳索取心钻探施工中钻杆折断原因分析及应对措施[J],探矿工程(岩土钻掘工程), 2009, (3): 15-18.
    蒋荣庆等.冲击回转多工艺钻进方法的新发展.长春地质学院建院40周年科学研究论文集(勘察技术与方法)[C].长春:吉林科学技术出版社, 1992. 15-18.
    兰洪波,高德利,张国辉. 149.3mm钻杆高强度接头数值分析[J].天然气工业, 2008, (9): 67-69.
    李传武,李发东,任海军.液动锤在科钻一井先导孔钻井中的应用[J].石油钻探技术, 2002, 30(5): 12-14.
    李海军.喷射器性能、结构及特殊流动现象研究[D],大连:大连理工学院,2004:36-42.
    李皓月,周田朋,刘相新. ANSYS工程计算应用教程[M].北京:中国铁道出版社, 2003. 1-5.
    李润方,林腾胶,施太和.石油钻杆连接螺纹弹塑性接触有限元分析[J].石油矿场机械, 1998, (5): 44-47.
    李晓芬.水基防塌钻井液在大庆松科1井的应用研究,《中国地质大学硕士论文》2008.
    李玉柱,范明顺.流体力学[M],北京:高等教育出版社,1998年6月,第一版,65-85.
    李玉柱.液体力学[M].北京:高等教育出版社。1998:79-184.
    李兆敏,蔡国琰.非牛顿流体力学,石油大学出版社,1998年4月,第一版,32-43.
    李子丰,马兴瑞,黄文虎.钻柱力学基本方程及其应用[J].力学学报; 1995, (4): 8-11.
    李子丰.油气井杆管柱力学[M].北京:石油工业出版社, 2008. 13-18.
    刘飞.离心压缩机的CFD分析与建模参数化[D].武汉:华中科技大学,2005:12-15.
    刘辉. S135钻杆失效分析[J].钻采工艺, 2009, (7): 63-66.
    刘希圣等.钻井工艺原理,石油工业出版社,1984年9月,第一版,46-89.
    刘远亮,乌效鸣,朱永宜,李晓芬.松科一井长裸眼防塌钻井液技术,《石油钻采工艺》.2009.08
    刘远亮.水平定向钻进铺管防冒将技术研究《中国地质大学硕士论文》2010.
    刘振铎等.刘广志文集[M].北京:地质出版社, 2003. 18-21.
    吕建国,桂暖银,胡远彪.钻杆失效理论分析方法初探[J].探矿工程(岩土钻掘工程), 2001, (5): 55-57.
    孟庆鸿,张恒春,胡郁乐,王稳石.防泥包钻头的优化设计与应用[J].煤田地质与勘探,2011.39(1),71-73
    孟庆鸿.干旱缺水地区水资源快速评价钻探技术的研究与开发,中国科技发展论坛,2005.06
    孟庆鸿.齐长缨等水资源快速钻探技术在山西探采结合孔的应用[J],探矿工程(岩土钻掘工程),2002.03
    孟庆鸿.西部特殊景观区深层地下水综合钻探技术论述[J],探矿工程(岩土钻掘工程),2007,34(9)24-27.
    冉诚实.钻挺螺纹失效分析及改进方法研究: [学位论文].西南石油大学, 2003.
    申昭熙.材料形变强化和摩擦系数对圆螺纹接头滑脱性能的影响[J].应用力学学报, 2008, (2): 293-296.
    沈玺,宋小青.工程钻进钻杆的实验模拟与仿真研究[J].地质装备, 2007, (4): 17-19.
    苏长寿.液动潜孔锤的现状及用于石油钻井应注意的几个问题[J].探矿工程(岩土钻掘工程), 2002, (6): 30-31;
    苏长寿.液动潜孔锤技术现状及发展设想[J].探矿工程(岩土钻掘工程), 2003, (1): 28-30.
    苏义脑,张国红.纵横弯曲法对变截面钻具组合的受力变形分析[J].石油钻采工艺, 1996, (2): 12-13.
    隋跃华,成效华,高建礼.强抑制性钻井完井液研究与应用.钻井泥浆与完井液,2001,18(6):11~14.
    孙丙伦.深部找矿组合钻探技术研究: [学位论文].吉林大学, 2009.
    唐继平等.塔里木山前复杂地层钻头使用技术[J].石油钻采工艺, 2003, 25(2): 39-43.
    陶勤南.农业试验设计与统计方法一百例[M]陕西科学技术出版社1987年9月473-569.
    陶兴华.提高深井钻井速度的有效技术方法[J].石油钻采工艺, 2001, 23(5): 4-8.
    王达,张伟,张晓西等中国大陆科学钻探工程科钻一井钻探工程技术,北京,科学出版社,2007,173-176.
    王关清,李克向,陈元顿等.中国陆上石油钻井科技发展回顾与展望[J].石油钻采工艺, 2000, 22(1): 1-6.
    王关清等.深探井和超深探井钻井的难点分析和对策探讨[J].石油钻采工艺, 1998, 20(1): 1-7.
    王国芳. S75绳索取心钻杆脱扣原因浅析[J].探矿工程, 1983, (2): 61.
    王国强.实用工程数值模拟技术及其在ANSYS上的实践[M].西安:西北工业大学出版社, 1999. 1-3.
    王良平,周世良.钻井泥浆处理剂及其作用原理.北京:石油工业出版社,2003.181~232.
    王人杰,苏长寿.我国液动冲击回转钻探的回顾与展望[J].探矿工程(岩土钻掘工程), 1999, (增刊): 321-323.
    王同良等.世界石油钻井科技发展水平与展望[J].石油钻采工艺, 2000, 22(2): 1-7.
    王同良等.中国钻井技术经济指标的发展规律研究[J].石油钻采工艺, 2001, 23(1): 1-8.
    王中华.钻井液技术现状及发展方向.断块油气田.2004(9),11(50):59~62;
    尉春玲,顾成林.钻杆的损坏原因及维护使用[J].西部探矿工程, 2008, (5): 59-61.
    魏存祥,刘海军,魏志.新型高强度钻杆性能及应用[J].重庆科技学院学报, 2008, (8): 33-34.
    乌效鸣,蔡记华,李云波.LG植物胶处理剂的试验研究.钻井泥浆与完井液,2005,22(1).
    乌效鸣,胡郁乐,贺冰新等.钻井泥浆与岩土工程浆液[M].武汉:中国地质大学出版社,2002.
    乌效鸣,朱永宜等.松科一井长裸眼防塌钻井液技术刘远亮[J],《石油钻采工艺》
    乌效鸣,朱永宜等.松科1井(主井)防塌钻井液技术蔡记华[J]《石油钻探技术》2008-09-29
    吴福云.S75绳索取心钻杆折断事故的原因分析[J].西部探矿工程, 2006, (10): 218-221.
    徐同台,赵忠举,徐云英.国内外钻井液技术新进展及对21世纪的展望(I).钻井液与完井液.2000,17(6),30~37.
    徐同台,赵忠举,徐云英.国内外钻井液技术新进展及对21世纪的展望(II).钻井液与完井液.2001,18(1),30~40.
    徐同台.井壁稳定技术研究现状及发展方向.钻井液与完井液.1997,14(4),36~43.
    许明财. ANSYS Workbench 10.0中文培训资料[M].北京: ANSYS-CHINA北京办事处, 2006. 3-12.
    闫洪泉.山东统计[J] 2001年第4期32-32.
    杨进能.钻杆偏磨的原因及其对疲劳寿命的影响[J].金属材料与冶金工程, 2008, (5): 16-18.
    杨自林.钻具失效事故的原因分析及对策[J].石油钻探技术, 2001, 24(5): 10-12.
    殷琨,蒋荣庆.潜孔锤反循环钻进技术及其应用[J].探矿工程, 1996, (5): 13-15.
    余伦秀.浅谈提高工勘中岩心采取率的几点措施[J]《四川地质学报》2009-06-15
    张春波.绳索取心钻杆损坏原因的分析[J].探矿工程, 1981, (5): 21-25.
    张春波等.绳索取心金刚石钻进技术[M].北京:地质出版社, 1985. 1-2.
    张洪武.有限元分析与CAE技术基础[M].北京:清华大学出版社,2004.39-63.
    张金昌,孟庆鸿.水文水井钻探技术现状及钻探新技术的应用,严重缺水地区地下水勘查论文集(2) 2004.10.
    张金昌,冉恒谦等.钻探技术面临的新形势、新机遇和新任务[J].探矿工程(岩土钻掘工程), 2007, (9): 13-15.
    张伟.2001年中国大陆科学钻探工程取得重要进展[J].探矿工程(岩土钻掘工程), 2002, (1),63.
    张伟.地质管材螺纹间隙及其计算[J].探矿工程(岩土钻掘工程), 1996, (2): 37-38.
    张晓静,乌效鸣等.LG植物胶泥浆在煤田复杂地层钻探中的应用研究探矿工程(岩土钻掘工程), 2007-12-25.
    张晓静.水敏/松散地层钻井液的护壁机理分析与应用研究,《中国地质大学硕士论文》2007.
    张炎,陈平,施太和.深井中钻挺螺纹连接的强度研究[J].探矿工程(岩土钻掘工程), 1998, (5): 21-23.
    张焱,曹里民,刘坤方等.一种计算钻柱螺纹弹塑性接触应力的新方法[J].石油钻采工艺1999,(5): 21-23.
    张毅,赵鹏,王治国.钻杆内加厚过渡区应力分布的有限元分析[J].试验与研究, 2001, (12): 11-15.
    张永勤.高效钻探技术是加速“危机矿山接替资源勘探规划”实施的最有效手段[J].探矿工程(岩土钻掘工程), 2006, (1): 6-8.
    张永勤.我国地质找矿取心(样)钻探设备现状及提高效能的分析研究[J].探矿工程(岩土钻掘工程), 2006, (8): 45-50.
    张云连.冲击旋转钻井工具的发展现状及应用研究[J].石油矿场机械, 2001, 30(5): 12-14.
    张作龙,高学仕,赵怀文等.两级柱螺纹套管接箍扣牙力计算与分析[J].石油大学学报(自然科学版); 1994, (5): 21-22.
    郑新权等.中国石油钻井技术现状及需求[J].石油钻采工艺, 2003, 25(2): 1-4.
    朱永宜,孟庆鸿等.矿产资源调查评价空气反循环钻探技术配套示范研究与开发,项目验收报告,2002.
    朱永宜,王稳石.松科一井(主井)取心钻进工艺[J].探矿工程(岩土钻掘工程),2008(9);1-5.
    朱永宜,王稳石.中国白垩纪科学钻探松科一井(主井)钻探工程概要[J].探矿工程(岩土钻掘工程), 2008.3.25.
    朱永宜,王稳石.中国陆相白垩纪科学钻探松科一井(北井)钻探工程技术配套[J],探矿工程(岩土钻掘工程), 2009增刊, (36): 388-390.

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700