摘要
本文在阅读大量国内外文献的基础上,研究了一套准确、经实践验证的混凝土防渗墙施工过程质量检测控制,成墙后垂直反射法、瞬变电磁法地面无损普查与少损层析成像详查相结合的混凝土防渗墙质量检测模式,并将研究成果应用于河北省黄壁庄水库副坝混凝土防渗墙质量控制及检测中,取得了良好的效果。为我国今后大量的混凝土防渗墙工程施工质量检测提供了技术保障。主要研究工作及成果如下:
(1)对混凝土防渗墙施工过程中常规质量检测指标和检测方法进行分析总结,进一步规范化;针对施工规范及设计提出的混凝土防渗墙垂直度的要求,分析研究施工槽孔孔斜率的检测方法并论证叠加孔斜率方法检测孔斜的严格性和可行性;根据设计提出的混凝土防渗墙厚度的要求和相交圆的性质,分析研究Ⅰ、Ⅱ期槽孔最小套接厚度的检测指标和计算方法。
(2)从弹性波传播的基本理论研究出发,推导了弹性波在类似混凝土防渗墙模型的二维垂向不均匀介质中传播的波动方程式;建立了混凝土防渗墙二维弹性薄板模型,研究弹性波在混凝土防渗墙中的传播、衰减、反射、透射规律及能量分配关系;从分辨率的概念出发,提出了垂直反射法检测混凝土防渗墙分辨率的内涵,对影响分辨率的因素进行了分析,推导出了缺陷层复合反射系数、可检测缺陷层极限厚度的表达式;建立了垂直反射法采样系统各环节之间的关系,研究分析了震源特性、接收传感器对测试结果的影响,提出了对测试仪器参数的要求及采集参数的选择原则。从理论、模型、分辨率、测试技术等方面,研究论证了垂直反射法检测混凝土防渗墙质量是一种实用、无损、效果明显的方法,作为一种质量普查方法是比较适宜的。
(3)针对混凝土防渗墙的特征,进一步完善瞬变电磁法检测防渗墙的理论体系、数据处理程序以及现场观测系统;利用等效代换的方法建立了瞬变电磁响应的概念、方法以及瞬变电磁场的扩散规律,推导计算出全区视电阻率及视深度参数的计算公式,进行混凝土防渗墙的二维电阻率成像,提高了检测结果定量分析的水平。
(4)分析总结了层析成像技术的发展及研究状况,检测混凝土防渗墙质量的物理基础及数学基础;归纳总结了各种射线追踪技术、及算法的优缺点,在此基础之上,编制层析成像直射线及弯曲射线正反演处理系统软件;依据混凝土防渗墙的特点,分别设计不同缺陷模型,进行正反演解析,确定了层析成像技术检测混凝土防渗墙缺陷的准确性、提出提高成像分辨率的措施;建立空洞模型,研究了层析成像最低分辨率的条件以及分辨能力与异常体尺寸、速度之间的关系,并用混凝土模型板试验验证了层析成像技术检测混凝土防渗墙缺陷的有效性;利用标准模块波速测试和抗压强度试验数据,建立了波速和混凝土强度之间的非线性回归方程,提高了层析成像检测混凝土防渗墙的定量化水平。
Based on abundant reading of domestic and foreign documentations, thisarticle researches on quality inspection pattern for concrete cutoff wallcombined with quality inspection control, vertical reflection, on surfacenon-destructive reconnaissance by transient electromagnetic methods and fewdestructive tomographic imaging. The research achievement has been appliedin quality control and inspection for concrete cutoff wall in attached dam ofHuangbizhuang reservoir in Hebei province, and good results have beenachieved. All this provides technical guarantees for construction qualityinspection for future large numbers of concrete cutoff wall projects. The mainresearches and achievements are as follows:
(1) Analyze and further formalize normal quality inspection index andinspection methods used in concrete cutoff wall construction; In accordancewith the verticality required by construction specification and design ofconcrete cutoff wall, test method for construction slot hole inclination areanalyzed and researched, the strictness and feasibility of hole inclinationinspection by stacked hole inclination method are demonstrated; Inspectionindex and calculation methods of minimum joint thickness of stageⅠ, stageⅡslot hole are analyzed based on design requirement of concrete cutoff wallthickness requirement and property of crossed circle.
(2) Based on basic theory of elastic wave propagation, wave equation forpropagation of elastic wave in 2 dimension vertical inhomogeneous mediumof similar concrete cutoff wall model is developed. 2 dimension elastic sheetof concrete cutoff wall is established to research law and energy distributionrelations between propagation, attenuation, reflection and transmission ofelastic wave in concrete cutoff wall; Based on the concept of resolution, thecontent of vertical reflection method used to inspect concrete cutoff wallresolution is raised and factors will have influences on resolution are analyzed,expression of defective layer composite reflection factor and detectabledefective layer maximum thickness are developed; Relations between allsteps of sampling system of vertical reflection are build, the influence ofcentrum characteristics and receiving sensor on testing results are analyzed,therefore parameter requirements for testing instrument and parameter selection principle of collected parameters are raised. Based on the theory,model, resolution and testing technology, it is demonstrated that verticalreflection method used for testing of the quality of concrete cutoff wall is apractical, nondestructive and effective method and suitable for qualityreconnaissance.
(3) In accordance with the characteristics of concrete cutoff wall, furthercomplete the theory system, data processing program and on situ observationsystem of transient electromagnetic method used in inspecting of cutoff wall;Concept and method of transient electromagnetic response and dispersion lawof transient electromagnetic field are established based on equivalentsubstation, formula of all area apparent resistivity and apparent depthparameters are developed and calculated, 2 dimensional resistivity imaging ofconcrete cutoff wall is conducted, the level of quantitative analysis of testingresults is improved.
(4) Development and researching conditions of tomographic imaging,physical and mathematic foundations of concrete cutoff wall qualityinspection are summarized; Advantages & disadvantages of different raystracking technology and calculation are summarized, and straight and curvedray forward and inversion processing software of tomographic image arecompiled based on that; In accordance with the characteristics of concretecutoff wall, models with various defects are designed for forward andinversion analysis, accuracy of tomographic image used to inspect defects inconcrete cutoff wall is determined, measures to increase image resolution areforwarded; Void model is established to research the conditions of minimumresolution of tomographic image and relations between resolution capabilityand size and speed of anomalous body, the effectiveness of tomographicimaging for defect inspection of concrete cutoff wall is verified by concretemodel test; Nonlinear regression equation of wave speed and concretestrength is established based on testing data of standard modular wave speedtesting and compressive strength test, quantification level of tomographicimaging used in inspecting concrete cutoff wall is improved.
引文
[1] 陶景良.混凝土防渗墙施工.北京:水利电力出版社,1988.1-20
[2] 高钟璞.大坝基础防渗墙.北京:水利电力出版社,2000.1-35
[3] 丛蔼森.防渗墙的设计施工与应用.北京:水利电力出版社,2001.1-60
[4] 段世超,王素静.小浪底大坝防渗墙设计.人民黄河,1996,7:55-58
[5] 中国水利水电基础工程局.水利水电工程混凝土防渗墙施工技术规范.北京:水利电力出版社,1996
[6] 许一欣.混凝土防渗墙在黑龙潭水库除险加固工程中的应用.北京:云南电力技术,2003,31(4):37-38
[7] 丁伯明.混凝土防渗墙是解决坝基渗漏的有效措施.河北水利水电技术,2001,4:37-38
[8] 陈庚仪.我国水工混凝土防渗墙技术的应用和发展.见:水利水电地基与基础工程新技术.天津:天津科学出版社,2002.3-22
[9] 宋功业,邵界立.混凝土工程施工技术与质量控制.北京:中国建材工业出版社,2003.1-28
[10] 冷元宝.工程CT技术在工程勘察中的应用.见:98水利水电地基与基础工程学术交流论文集.天津:天津科学技术出版社,1999.241-243
[11] 段世超,孙胜利,杨永叶.小浪底大坝混凝土防渗墙质量检测.西部探矿工程,2000,66(5):36-38
[12] 梁国钱,苏全.层析成像技术在拱坝质量检测中的应用.水利发电,2002,(5):20-22.
[13] 朱介寿,严忠琼,曹俊兴.井间地球物理层析成像软件系统研究.物探化探计算技术,1994,16(4):310-320
[14] 靳洪晓,赵永贵等.地面地震CT在浅层勘探中的应用.工程地质学报,2002,8(2):239-243
[15] 王辉,黄鼎成.地震层析成像方法及其在岩体结构研究中的应用.工程地质学报,2000,8(1):109-117
[16] 宋先海,黄荣荣.弹性波CT在大体积混凝土无损检测中的应用.人民长江,2003,34(11):23-25
[17] 王家映.地球物理反演理论.北京:高等教育出版社,2002.1-181
[18] LanczosC. Aniteration methodforthesolution of the eigenvalue problem olinear differential andinterral operators. J. Res. N.B.S.,1950,45:255-282
[19] SassaKoichi. Suggested methodforseismictestingand between borfholes Part4:Suggested method forseis mictomography. International Journal Rock Mechanic sand Mining sciences&Geomechanics Abstracts, 1988,25:447-472
[20] 杨文采.地球物理反演的理论与方法.北京:地质出版社,1993.1-166
[21] 王运生.最佳路径算法在计算波路中的应用.物探化探计算技术,1992,14(1):32-36
[22] LanganRT, etal. Tracingofrays through heterogeneousmedia:Anaccu2ratean deficient procedure. Geophysics, 1985,50:1456-1465
[23] 牛建军,杜立志,谷城.岩溶探测中的弹性波CT方法.吉林大学学报,2004,34(4):630-633
[24] 王兴泰.工程与环境物探新方法新技术教程.北京:地质出版社,1996.1-38
[25] 刘合,王玉普.国外井间地震技术.北京:石油工业出版社,1998.1-56
[26] 邱庆程,李伟和.跨孔地震CT层析成像在岩溶勘察中的应用.物探与化探,2001,25(3):236-240
[27] 金亚荣,陈国平.层析技术在岩溶地区岩土工程勘察中的应用.中山大学学报,2000,39(2):114-118
[28] 曹小林,严忠琼,朱介寿.井间地震层析成像技术及其在岩溶勘察中的应用.地质灾害与环境保护,1997,8(2):53-58.
[29] 黄金莉等.工程物探中层析成像的研究.地球物理学报,1994:37(2):682-686
[30] 杨文采,杜剑渊.层析成像新算法及其在工程检测中的应用.地球物理学报,1994,37(1):239-244
[31] Lanczos C. Aniteration method for the solution of the eigenvalue problem of linear differ eniatal and integral operaters. J. Res. N.B.S., 1950,45:255-281
[32] Paige C C, Sounders M A. Anal gorithm for sparse line are quation and sparse least squares, ACMTrans. Math. Softw. 1982,8:195-209
[33] 孙跃军,俞国柱,石桂.井间地震层析成像技术及其在岩溶勘察中的应用.CT理论与应用研究,2001,10(4):10-13
[34] 刘江平.浅层地震技术在大堤防渗墙质量检测中应用.物探与化探,2000,24(4):303-305
[35] 杨成林.瑞雷波勘探法原理及其应用.物探化探,1989,13(6):465-468
[36] 何樵登,熊维纲.应用地球物理教程.北京:地质出版社,1991.1-45
[37] 夏唐代,林水珍.波速法在防渗墙质量检测中应用研究.沈阳化工学院学报,2002,14(4):273-276 [38] NazarianS. Insi2tu Determination of Elastic Moduli of Soil DePosits and Pavement System by Spectral 2Analysis2 of 2surface Wave method[D].University of TexasofAustin, 1984
[39] Gucunslin, WoodsRD. Use of Rayleigh Modesin Inter pretation of SASWTest. Second International Conference on Recent Advance sin Geotechnical Earth quake Engineering and SoilDynamics[Z].St. Louis, Missout, 1991.1399-1408.
[40] 夏唐代,陈云敏,吴世明.匀质软夹层地基瑞利波弥散特性.振动工程学报,1993,6(1):42-50
[41] Eringen AC,舒胡华SA.弹性动力学,第二卷.戈革译.北京:石油工业出版社,1984.1-124
[42] 罗莎,闫海新,李广智.瑞雷面波在防渗墙无损检测中应用.地球物理学进展,2003,18(3):416-419
[43] 余才盛.堤防防渗芯墙质量检测技术及效果.地球物理学进展,2003,18(3):410-415
[44] 何继善等.可控源音频大地电磁法.长沙:中南工业大学出版社,1989.1-45
[45] 王山山,任青文.黄河大堤防渗墙质量无损检测方法研究.河海大学学报,2004,32(4):405-409
[46] 李力,侯建国,杨成林等.探地雷达和瑞雷波方法在水工混凝土质量检测中应用.水利水电技术,2001,32(5):35-39
[47] 王世恩,刘超常,锅伯强.用瑞利波检验高喷防渗墙的施工质量及地层分布特征工程勘探,1998,(6):62-64
[48] RENWei2xin, GUIDODR. Structural damage identification using modal data. 1: simulation verification. ASCE, Journal of Structural Engineering, 2002, 128(1):87-95
[49] RENWei2xin, GUIDODR. Structural damage identification using modal data. 2: test verification. ASCE, Journal of Structural Engineering, 2002,128 (1):96-104
[50] 王山山,江泉,马良筠.用振动法检测特殊桩基的质量.力学与实践,2001,23(4):19-20
[51] 蔡靖,张献民,王建华等.面波与振动映像结合的防渗墙无损检测技术.水利水电技术,2005,25(2):40-42
[52] 王建华,蔡靖,孙建运.水泥搅拌土防渗墙中瑞利面波特性的研究.第九界全国土动力学及基础工程学术会议论文集编委会.北京:清华大学出版社,2003
[53] 高建东.使用地质雷达探测水库地下防渗墙.水利水电技术1999,35(2) :115-117
[54] 张升光,侯建国,杨成林等.探地雷达和瑞雷波法在水工混凝土土质量检测中的应用.水利水电技术,2001,32(5):35-39
[55] 肖自龙,李宏展,梁娜.透地雷达在检测塑性混凝土垂直防渗墙中的应用.华北水利水电学院学报,2005,26(1):39-41
[56] 李大心.探地雷达法与应用.北京:地质出版社,1994.1-80
[57] 中国水利水电基础工程局主编.水利水电工程混凝防渗墙施工技术规范(SL174—96),北京:中国水利水电出版社,1996
[58] 中华人民共和国水利电力部主编.水工混凝土施工规范范(SDJ—207—82).北京:水利电力出版社,1982
[59] 中华人民共和国水利电力部主编.水工混凝土施工规范(SD105—82).北京:水利电力出版社,1983
[60] 谢里夫 RE,吉尔达特 LP.勘测地震学(第二版).北京:石油工业出版社.1999.1-24
[61] 王振东编.浅层地震勘探应用技术.北京:地质出版社,1988.55-67
[62] 吴小平,赵鸿儒.我国工程多波地震勘探研究与应用现状.地震学报,1996,18(2):262~268
[63] 俞寿朋.高分辨率地震勘探技术.北京:石油工业出版社,1993.1-135
[64] 徐攸在,刘兴满.桩的动测新技术,北京:中国建材工业出版社,1989.1-15
[65] 王雪峰,吴世明.基桩动测技术,北京:科学出版社,2001.8-57
[66] 王礼立,朱兆祥.应力波基础,北京:国防工业出版社,1984.1-25
[67] Rayleigh, Lord, On Waves Propagated along the Plan Surface of an Elastic Solid, Proc. London Math. SOC.,17(1887)4
[68] Rayleigh, Lord. Theory of Sound, Macmillan, London(1894)
[69] 王俊茹编著.工程与环境地震勘探技术.北京:地质出版社,2002.1-23
[70] Pao Y H. Elastic waves in solids, ASME Trans. Journal of Applied Mechanichs, 1983,50:1152~1164
[71] 田钢,刘菁华,曾绍发.环境地球物理教程.北京:地质出版社,2005.1-72
[72] Sheriff R E. Encyclopedic dictionary of exploration geophysics(3~(RD)edition).Tulsa:USA SOC. Expl. Geophys, 1991
[73] Widess M B. How thin is a thin bed? Geophysics, 1973,38(46):1176-1180
[74] Farr J B. High-resolution seismic methods improve stratigraphic exploration. Oil&Gas Journal, 1977,75(48):182-188
[75] 唐文榜.全波震相在地震勘探中的应用前景.地球物理学进展,1989,4(1):6 —8
[76] 赵鸿儒,孙进忠,唐文榜.全波震相分析的应用.地球物理学报,1990,33(1):54-61
[77] 赵鸿儒,郭铁栓,徐子君等.工程多波地震勘探.北京:地震出版社,1996.13-53
[78] 薛云峰,袁江华,孙晓暾.垂直反射法检测砼防渗墙的研究应用.物探探,2004,28(5):467-470
[79] 何汉漪主编.海上地震资料高分辨率处理技术论文集.北京:地质出版社,2001.1-16
[80] 何汉漪.高分辨率地震的优势与误区.中国海上油气,1999,13(6):439-442
[81] 李庆忠.走向精确勘探的道路.北京:石油工业出版社,1993.17-23
[82] Widess, M,B. Quantifying resolving Pover of seismic systems. Geophysics, 1982,47:1160~1173
[83] 何汉漪.海上高分辨率地震实例研究.地球物理学报,1999,42(1):121-125
[84] 王雪峰,吴世明,张岚.基桩完整性检测中的振源特性分析.工程勘察,2000,6:32-34
[85] 王雪峰,吴世明.不同传感器在桩基完整性检测中的应用.工程勘察,2000,4:35-37
[86] 罗骐先.桩基工程检测手册.北京:人民交通出版社,2003.78-83
[87] 王义江,于永志.速度计和加速度计性能参数对基桩动测结果的影响.工程质量,2001,9:27-29
[88] 牛之琏,李爱华.瞬变电磁在铝锌金矿床上的找矿效果.地质与勘探,1991,11:41-43
[89] 李貅,郭文波等.瞬变电磁法在煤田矿井涌水通道勘察中的应用.西安工程学院学报,2000,22(3):451-454
[90] 戴定旋,于进,马长金等.瞬变电磁法在工程地质勘察中的应用.物探化探,1995,19(5):385-389
[91] 葛云生,海戴嫒.珠海磨刀门大桥桥址瞬变电磁勘探.桂林冶金地质学院学报,1992,129(增刊):25-28
[92] M.N Nabighian. Quasi-static transient response of a conducting half space-An approximate representation. Geophysics, 1979,44(10)
[93] M.N Nabighian. Electromagnetic methods in applied geophysics—Theory Volume I, SeG, 1988
[94] 牛之琏编著.时间域电磁法原理.长沙:中南工业大学出版社,1992.1-28
[95] Raab P. and Frisehknecht F. Desktop Computer Processing of Concident and Central loop time domain electromagnetic data. U.S.G.S Open-File Roport: 83-240
[96] 薛云峰,孙晓暾,孙雅芳.瞬变电磁电阻率成像法探测岩体渗漏通道的研究应用.水文地质工程地质,2004,31(197):96-98
[97] 李貅.瞬变电磁测深的理论与应用.西安:陕西科学技术出版社,2002.79-84
[98] A. M Cormack. Representation of a faaction by its line ittegrals with som radiologioal application.J. Appl. Phys. 1963, 34:2722-2727
[99] 孙坚.层析成像技术在土木工程中的应用研究[博士论文].杭州,浙江大学建筑工程学院,1999
[100] Chapman C.H. Genneralimed Redom transforms and slant stackes. Geophyys, J, R, Astr, soc, 1981,66
[101] Bois P.,etal. Well-to well seismic measurements. Geophysics, 1972, 37:471-480
[102] Dines K A, Lytle RJ. Computerined geophysical tomography. Proc, IEEE, 1979,67:1065-1073
[103] Andersion D.L. Seismic tomography of the Earth's interior. American scientist, 1979,67:1063-1073
[104] Daily W D. Underground oil shale retort monitoring using geotomography. Geophysics, 1984,172:1701-1711
[105] Somerstein, Radio-frequencey geotomography for remotely.Probing the interior of operating mini and commercial-sined oil-shale retorts. Geophysics, 1984,49:1288-1300
[106] Pratt R G,Worthington M H. The application of diffraction tomography to crossbole seismic data. Geophysics, 1988,53:1284-1294
[107] Bishop T N, et al. Tomographic determination of velocity and depth in laterally varying media. Geophysics, 1985,50:903-923
[108] 杜剑渊,杨文采.跨孔地震层析成像的级联方法.地球物理学报,1991,31(6):441-443
[109] 杨国宪,张可霓,李丕武.层析成像在黄河小浪底工程地下厂房区域的应用研究.地球物理学报,1993,36(1):106-112
[110] 杨文采.应用于地质勘探和矿产开发的地表层析成像.中国地质,1991,(7):22-26
[111] 杨文采,杜剑渊.煤矿陷落柱地震层析成像试验.物探与化探,1992,16(3):174-180
[112] 王运生,冷元宝.地震折射成像在水库工程中的应用.勘察科学技术,1992(1):12-14
[113] Julian, B.R, Gubbins, D. Three-dimensional seismic ray tracing.J. Geophysics. 1997,43:95-114
[114] Um, J. and Thurber,C.A fast algorithm for two Points seismic ray tracing.Bull. seis. SOC. am.,1987,79:927-986
[115] Lee, J.M. and Shalev, E, On the stability of P-wave tomography at Loma Prieta:A comparison of parameterinations linear and nonlionear inversions. Bull. Seis. Soc. Am.,1992,82:1821-1839
[116] Vidal, J.,Finite-difference calculation of travel times. Bull. Seis. Soc. Am.,1988,78:2062-2076
[117] Qin, F.H. etal.,Solution of the eironal equation by a finite differencemethod, 6oth SEG, 1990
[118] Vetle Vinje. Einar Iversen, Havar Gjoystaal. Traveltime and amplitude estmatio using wavefront construction. Geophysics, 1993,58(8):1157-1166
[119] 张建中,陈世军,余大祥。最短路径射线追踪方法及其改进,地球物理学进展,2003,18(1):146-150
[120] 高尔根,徐果明,蒋先艺等.三维结构下逐段选代射线追踪方法.石油地球物理勘探,2002,37(1):11-16
[121] 聂建新,杨慧珠.地震波旅行时二次1线性联合插值法.清华大学学报(自然科学版),2003,43(11):1495-1498
[122] 宋炜.球面波近似快速波前射线追踪.勘探地球物理进展,2003,26(4):255-259
[123] 刘光鼎,常旭.浅层地球物理工程.科学,1993,45(4):16-19
[124] Youli Quan, Jerry M. Harris, Seismic attenuastion tomography using the freauency shift method, Geophysics, 1997,63(3):895-905
[125] Kjartansson E., Constant Q-wave Propagation and attenuation.,J..Geophys. Res.,1979,82(7): 4737-4748
[126] D.H.Johnston等.有关衰减的基本知识.石油物探译丛,1983,38(9):845-854
[127] 左德.地震衰减层析成像理论及实用性探讨.煤炭工程师,1994,(2):48-49
[128] 赵永贵,李勤,郭鸿等.频率域地震吸收CT方法及其工程应用.中国科学,2002,30(1):107-112
[129] 李丕武.黄河物探技术应用与发展.郑州:黄河水利出版社,1996.485-541
[130] 肖柏勋.工程地球物理学进展.武汉:武汉水利电力大学出版社,2000.1-8
[131] 夏可凡等.水利水电地基与基础工程新技术.天津:天津科学技术出版社,2002.526-528
[132] 曹俊兴,严忠琼.地表波跨孔旅行时层析成像分辨率的估计.成都理工学院学报,1995,22(4):290-293
[133] 陈国金.井间观测系统的讨论.石油物探,1996,35(3):345-34
[134] 裴正林,余钦范,狄帮让.井间地震层析成像分辨率研究.物探与化探,2002,26(3):276-278
[135] 骆循.ART算法及其在地球物理反演中的作用.物探化探计算技术,1987,9(1):36-38
[136] Gordon. R.,Bender, R. and Herman, G.T. (1970).Algebraic Reconstruction Techniques(ART) for three-dimensional electron microscopy and X-ray photography, J. Theor.Bio1.1970,29:471-481
[137] Clement, W.P. and M.D. Knoll. Tomographic inverseon of crosshole radar Data: Confidence in results,.Proceedings of the Symposium on the Application of Geophysics to Environmental and Engineering Problems. Arlington, VA, 2000,20:553-562
[138] GROUST NOLET. Seismic Tomography with applications in Global Seismology and Exploration Geophysics. D. REIDEL PUBLISHING COMPANY, 1987
[139] [美]G.T.赫尔曼.由投影重建图像,CT理论基础.北京:科学出版社,1985.1-74
[140] 王春镛,李幼铭.地震层析成像及其应用.北京学术书刊出版社,1989.1-56
[141] 王运生.弯曲射线地表波透射层析成像的一种实现方法.河海大学学报,1993(40):227-230
[142] R.T. Langan. Tracing of rays through heterogencous media: An accurate and efficient procedure. GEOPHYSICS, 1985,50(9):116-127
[143] Moser. T. J. Shortest path ealculation of seismic rays. GEOPHYSICS, 1991, 56(1):59-67