浅水域堆积体精准探测技术应用研究
|
摘要
对水下堆积体的物质组成及其分布特征等内容的准确了解,往往是水利工程建设与维护的重要基础工作,现多采用综合地球物理勘查(简称物探)方法,确定堆积体的分布范围与深度,再结合相应的钻探,对物探成果进行验证并进一步查明堆积体的物质组成。在浅水区,受水体深度、目标体的结构特征、波浪等外界因素的干扰,使各种物探方法的有效性受到明显影响；不同时期堆积体其组份差异往往较大,由此所构成的堆积体,其物理性质的复杂多变,给物探成果的准确解译形成很大的障碍。
以淮河中游干流护岸抛石和大型水库除险抛石为主要研究对象,针对高密度电阻率法、探地雷达法和高密度地震映像法这三种目前主要的工程物探方法,利用不同水深、不同厚度堆积体的数值模拟和现场实测成果,进行有效性对比研究；针对堆积体物性复杂多变的特点,对物探解译方法进行相应的改进；研究过程中,结合实际工作中所常遇的复杂工况,对仪器设备进行适用性改良。所获得的创新成果：(1)对高密度地震映像方法,研制了一种轻型气压式浅水域全自动触发器设备,在经济合理的前提下,有效克服了大能量震源的气泡效应。(2)在综合研究高密度电法、探地雷达和高密度地震映像法在浅水域抛石体探测中所表现的特征和规律的基础上,建立了浅水域抛石体探测的地电模型,并进行了正反演模拟,为浅水域堆积体精准探测应用提供理论支撑；(3)根据浅水域堆积体上覆水体的不同深度,堆积体的结构特征和赋存状态,探索精准探测方法,明确了对应的有效物探方法组合；上述工作,将为不同物探方法的有效选择与野外工作的高效实现提供技术支撑。
The exploration of accumulation bodies' composition and its distribution characteristics is always the basis for hydraulic engineering construction and maintenance. The integrated geophysical survey is applied mostly to deduce the distribution range and depth of accumulation bodies. Then, the results of geophysics is verified and the objects' materials are discovered according to drilling. However, almost all of the geophysical methods are affected by external factors, such as water depth, structure characteristics of the target body and wave. And it is believed that the material difference in the accumulation bodies during various periods is often larger, leading to difficulty of processing geophysical interpretation.
The paper takes dumped riprap in the middle reaches of Huaihe River and danger elimination riprap as main research objects. In view of the high density resistivity method, ground penetrating radar and seismic imaging as the three major engineering geophysical exploration methods, the targets' density, velocity and dielectric constant the forwarding computing is conducted where the physical parameters of water depths change,.
The geophysical interpretation technology according to the complicated property of the accumulation bodies is improved for the complex accumulation body physical characteristics. And the geophysical equipment's applicability is also improved in order to adapt to the complex condition during surveying. The innovations are:(1) It improves a light pressure type automatic trigger device which could be used in shallow waters during high density seismic image surveying, and the improved device will overcome the bubble effect of the large energy source.(2) On the basis of the characteristics of high density resistivity method, ground penetrating radar and high density seismic image method, a shallow waters riprap model is established, which will provide theoretical support for accurately detecting accumulation bodies under the premise of economic rationality.(3) A precise survey method according to different water depth, structure characteristics and occurrence status of accumulation bodies is probed, the corresponding effective combination geophysical surveying is also clear out. The above results will be effective for different geophysical exploration method selection and field work efficiently to provide technical support.
以淮河中游干流护岸抛石和大型水库除险抛石为主要研究对象,针对高密度电阻率法、探地雷达法和高密度地震映像法这三种目前主要的工程物探方法,利用不同水深、不同厚度堆积体的数值模拟和现场实测成果,进行有效性对比研究；针对堆积体物性复杂多变的特点,对物探解译方法进行相应的改进；研究过程中,结合实际工作中所常遇的复杂工况,对仪器设备进行适用性改良。所获得的创新成果：(1)对高密度地震映像方法,研制了一种轻型气压式浅水域全自动触发器设备,在经济合理的前提下,有效克服了大能量震源的气泡效应。(2)在综合研究高密度电法、探地雷达和高密度地震映像法在浅水域抛石体探测中所表现的特征和规律的基础上,建立了浅水域抛石体探测的地电模型,并进行了正反演模拟,为浅水域堆积体精准探测应用提供理论支撑；(3)根据浅水域堆积体上覆水体的不同深度,堆积体的结构特征和赋存状态,探索精准探测方法,明确了对应的有效物探方法组合；上述工作,将为不同物探方法的有效选择与野外工作的高效实现提供技术支撑。
The exploration of accumulation bodies' composition and its distribution characteristics is always the basis for hydraulic engineering construction and maintenance. The integrated geophysical survey is applied mostly to deduce the distribution range and depth of accumulation bodies. Then, the results of geophysics is verified and the objects' materials are discovered according to drilling. However, almost all of the geophysical methods are affected by external factors, such as water depth, structure characteristics of the target body and wave. And it is believed that the material difference in the accumulation bodies during various periods is often larger, leading to difficulty of processing geophysical interpretation.
The paper takes dumped riprap in the middle reaches of Huaihe River and danger elimination riprap as main research objects. In view of the high density resistivity method, ground penetrating radar and seismic imaging as the three major engineering geophysical exploration methods, the targets' density, velocity and dielectric constant the forwarding computing is conducted where the physical parameters of water depths change,.
The geophysical interpretation technology according to the complicated property of the accumulation bodies is improved for the complex accumulation body physical characteristics. And the geophysical equipment's applicability is also improved in order to adapt to the complex condition during surveying. The innovations are:(1) It improves a light pressure type automatic trigger device which could be used in shallow waters during high density seismic image surveying, and the improved device will overcome the bubble effect of the large energy source.(2) On the basis of the characteristics of high density resistivity method, ground penetrating radar and high density seismic image method, a shallow waters riprap model is established, which will provide theoretical support for accurately detecting accumulation bodies under the premise of economic rationality.(3) A precise survey method according to different water depth, structure characteristics and occurrence status of accumulation bodies is probed, the corresponding effective combination geophysical surveying is also clear out. The above results will be effective for different geophysical exploration method selection and field work efficiently to provide technical support.
引文
[1]中国工程院“21世纪中国可持续发展水资源战略研究”项目组.中国可持续发展水资源战略研究综合报告[R].中国工程科学,2000,2(8)：1-17.
[2]钱正英.中国水资源战略研究中几个问题的认识[J].河海大学学报(自然科学版),2001,29(3)：1-7.
[3]程晓陶.新时期大规模的治水活动迫切需要科学理论的指导—论有中国特色的洪水风险管理[J].水利发展研究,2001,1(4)：1-6.
[4]董哲仁.堤防除险加固实用技术[M].北京：中国水利水电出版社,1998.
[5]中华人民共和国水利部和国家统计局.第一次全国水利普查公报[R].2013.
[6]冷元宝,黄建通,张震夏等.堤坝隐患探测技术研究进展[J].地球物理学进展.2003,18(3)：370-379.
[7]徐晓东.安徽省长江2011年度河势演变分析及崩岸监测预警报告[R].安徽长江河道设计研究院,2011.
[8]马洪蛟,李宇,蔡辉等.探地雷达探测抛填厚度研究[J].南京林业大学学报(自然科学版),2004,28(2)：51-53.
[9]葛双成,张孚专,孙伯永等.洞头县杨文围涂东围堤抛石层雷达检测与分析[J].浙江水利科技,2006,(3)：35-36.
[10]孔宪立.工程地质学[M].北京：中国建筑工业出版社,2003.
[11]张玉,朱维申.小湾水电站左岸坝前堆积体在自然状态下稳定性的平面离散元与有限元分析[J].云南水力发电,2000,16(1)：36-39.
[12]唐川,黄楚兴,万哗.云南省丽江大地震及其诱发的崩塌滑坡灾害特征[J].自然灾害学报,1997,6(3)：76-84.
[13]邵景林,祁建忠,高杰.西气东输管道通过长江工程地质勘察水上钻探工艺合理性分析[J].工程勘察,2002(2)：52-53.
[14]高建波.浅谈水上地质钻探方法选择与应用[J].铁道勘察,2012,3：51-54.
[15]季喜平.轻型静力触探技术及其在中小型水利工程中的应用[J].水利水电科技进展,1998,18(4)：4749.
[16]Waite A D.王德石,译.实用声纳工程[M].北京：电子工业出版社,2004.
[17]吴自银,郑玉龙,初凤友,等.海底浅表层信息声探测技术研究现状及发展[J].地球科学进展,2005,20(11)：1210-1217.
[18]Jun Chen, Hong Jie Wan, Xue Wei Wang. A design of underwater acoustic responder based on spread-spectrum communication[J]. Procedia Engineering,2012,29(1):128-132.
[19]Dario Pompili, Tommaso Melodia, Ian F. Akyildiz. Three-dimensional and two-dimensional deployment analysis for underwater acoustic sensor networks[J]. Ad Hoc Networks,2009,7(4): 778-790.
[20]John J. Schultz, Carrie A. Healy, Kenneth Parker, Bim Lowers. Detecting submerged objects: The application of side scan sonar to forensic contexts[J]. Forensic Science International,2013, 231(1-3):306-316.
[21]B. Lafferty, R. Quinn, C. Breen. A side-scan sonar and high-resolution Chirp sub-bottom profile study of the natural and anthropogenic sedimentary record of Lower Lough Erne, northwestern Ireland[J]. Journal of Archaeological Science,2006,33(6):756-766.
22]褚宏宪,白大鹏,史慧杰等.海底抛石声学探测方法-以渤海湾人工岛抛石检测为例[J].海洋地质前沿,2011,27(4)：65-70.
[23]Vladimir Shtivelman. Shallow water seismic surveys for site investigation in the Haifa Port Extension area, Israel[J]. Journal of Applied Geophysics,2001,46(2):143-158.
[24]T. Hussein, A.G. Mula, H. Schneider. Geological and seismic investigations with regard to shallow ground water exploration in Eastern Sudan Republic[J]. Journal of African Earth Sciences (and the Middle East),1989,8(1):75-78.
[25]Shtivelman, V.. Using surface waves for estimating shear wave velocities in the shallow subsurface onshore and offshore area[J]. Eur. J. Environ. Eng. Geophys.,19994,15-35.
[26]Shtivelman, V., Frieslander, U., Zilbennan, R., Amit, R.. Mapping shallow faults at the Evrona playa site using high resolution reflection method[J]. Geophysics,1998,63:1257-1264.
[27]王传雷,曲赞,祁明松等.高精度磁测在长江马当要塞沉船探测中的应用[J].地质科技情报,2000,19(3)：98-102.
[28]喻忠鸿,王传雷,吴文贤.磁法探测炸弹有效深度的物理模拟试验及分析[J].工程地球物理学报,2007,4(2)：118-122.
[29]Dwain K. Butler. Potential fields methods of location of unexploded ordnance[J]. US Army Engineer Research and Development Center, The leading EDGE,2001,20(8):890-895.
[30]吕邦来,王传雷,曲赞.磁测在港航工程勘察中的几个实例[J].工程地球物理学报,2004,1(6)：494498.
[31]Chuanlei Wang. The Magnetic Detection of Sunken Ships in the Madang Section of the Yangtze River[J]. Journal of Environmental & Engineering Geophysics (USA),2006,11(2):123-131.
[32]Zeev Zalevsky, Yuri Bregman, Nizan Salomonski, Hovav Zafrir. Resolution Enhanced Magnetic Sensing System for Wide Coverage Real Time UXO Detection[J]. Journal of Applied Geophysics,2012,84(1):70-76.
[33]SL326-2005,水利水电工程物探规程[S].
[34]赵铁虎,张志珣、许枫.浅水区浅地层剖面测量典型问题分析[J].物探化探计算技术,2002,24(3)：215-216.
[35]滕惠忠,严晓明,李胜全.侧扫声纳图像增强技术[J].海洋测绘,2004,24(2)：47-48.
[36]李金铭,罗延钟.电法勘探新进展[M].北京：地质出版社,1996.
[37]严加永,孟贵祥,吕庆田等.高密度电法的进展与展望[J].物探与化探,2012,36(4)：576-583.
[38]罗延钟,王传雷,董浩斌.府密度电阻率法的电极装置选择[J].地质与勘探,2005,增刊：174-178.
[39]Edwards R N. On the resource evaluation of marine gas hydrate deposits using sea-floor transient electric dipole-dipole method [J]. Geophysics,1997,62:63-74.
[40]王 剑,张赛民,马晓成.水上高密度电法在油气管道工程中应用[J].吉林地质,2012,31(2)：76-79.
[41]吴长盛.北大港水库堤坝裂缝检测与评定技术研究[J].水利水电技术,2001(5)：61-63.
[42]王士鹏.高密度电法在水文地质和工程地质中的应用[J].水文地质工程地质,2000(1)：52-56.
[43]王传雷,董浩斌.科学技术部重点项目、武汉市配套项目“堤坝隐患探测技术研究”(1999—2001年)[R].中国地质大学,2001.
[44]Cheng Q, Tuo XG, Ge B, Li HL, Wang SD. The application of the high-density electrical method to the survey of Chayuangou landslide in Gaochuan, Sichuan province[J]. Geophysical and Geochemical exploration,2012,36:69-72.
[45]Osazuwa IB, Abdullahi NK.2-D Electrical Resistivity and Induced Polarisation Investigation at an Open Solid Waste Dumpsite:Case Study from Kaduna, North Central Nigeria[J]. Journal of Environmental Hydrology.2008,16:1-11.
[46]Al-Tarazi E, Abu Rajab J, Al-Naqa A, El-Waheidi M. Detecting leachate plumes and groundwater pollution at Ruseifa municipal landfill utilizing VLF-EM method[J]. Journal of Applied Geophysics,2008,65:121-131.
[47]Jegede SI, Osazuwa IB, Abdullahi NK. Geoenvironmental Study of groundwater contamination in a dual aquifer environment using earth resistivity imaging[J]. Journal of American Science, 2011,7:367-377.
[48]Fanny Genelle, Colette Sirieix, Joelle Riss, Veronique Naudet. Monitoring landfill cover by electrical resistivity tomography on an experimental site[J]. Engineering Geology,2012, 145-146:18-29.
[49]Luciana Orlando. Some considerations on electrical resistivity imaging for characterization of waterbed sediments[J]. Journal of Applied Geophysics,2013,95:77-89.
[50]A. Kampke. Focused imaging of electrical resistivity data in archaeological prospecting[J]. Journal of Applied Geophysics,1999,41(2-3):215-227.
[51]A. Samouelian, I. Cousin, A. Tabbagh, A. Bruand, G. Richard. Electrical resistivity survey in soil science:a review[J]. Soil and Tillage Research,2005,83(2):173-193.
[52]Dale F. Rucker, Marc T. Levitt, William J. Greenwood. Three-dimensional electrical resistivity model of a nuclear waste disposal site[J]. Journal of Applied Geophysics,2009,69(3-4): 150-164.
[53]王兴泰,李晓芹.电阻率图像重建的佐迪CZo-hdy)反演及其应用效果[J].物探与化探,1996,20(3)：228-233.
[54]张大海,王兴泰.二维视电阻率断面的快速最小二乘反演[J].物探化探计算技术,1999,21(1)：2-8.
[55]Dw groothedlinc. constable s. Occam inversion to generate smooth, two-dimensional models form magnetotelluric data [J]. Geophysics,1990,55:1613-1624.
[56]LOKE M H, BARKER R D. Rapid least-squares inversion of apparent resistivity pseudo sections by a quasi-Newton method[J]. Geophysical Prospecting,1996,44:131-152.
[57]吕邦来,匡翠萍.浅地层剖面探测与地震映像法探测的应用对比和分析[J].港工技术,2009,46(B07)：86-89.
[58]Ricker. The form of laws of propagation of seismic wavelets[J]. Geophysics,1953,18(1): 10-40.
[59]A.T. Dennison. An introduction to synthetic seismogram techniques [J]. Geophysics Prospecting, 1960,8(2):231-241.
[60]P. C. Wuenschel. Seismogram synthesis including multiples and transmission coefficients [J]. Geophysics,1960,25(1):106-129.
[61]付燕.人工地震信号去噪方法研究[D].西安：西北工业大学,2002.
[62]杨农合,徐小林,吴朝俊等.应用瞬变电磁法小线框大电流探测水下构造[J].工程地球物理学报,2012，2,9(3)：279-284.
[63]李大心等.探地雷达探测方法与应用[M].北京：地质出版社,1994.
[64]杜树春.探地雷达及其在环境地质中的应用[J].物探与化探,1996,20(5)：384-392.
[65]李大心等.江河堤防隐患的探地雷达调查[J].中国地质灾害与防治学报,1996,7(1)：21-25.
[66]武可贵,程学标.探地雷达在工程地基检测中的应用研究[J].水运工程,1997(9)：5-8.
[67]邓世坤.探地雷达用于深圳湾软基处理工程中的质量检测[J].岩土工程学报,2000,22(1)：33-35.
[68]雷林源.探地雷达应用中的几个基本问题[J].物探与化探,1998,22(6)：408-414.
[69]王水强,万明浩,谢雄耀等.不同探地雷达测量参数对数据采集效果的影响[J].物探与化探,1999,23(3)：211-215.
[70]Timo Saarenketo, Tom Scullion. Road evaluation with ground penetrating radar[J]. Journal of Applied Geophysics,2000,43(2-4) 119-138.
[71]J.A. Doolittle, B. Jenkinson, D. Hopkins, M. Ulmer, W. Tuttle. Hydropedological investigations with ground-penetrating radar (GPR):Estimating water-table depths and local ground-water flow pattern in areas of coarse-textured soils[J]. Geoderma,2006,131(3-4):317-329.
[72]Jianghai Xia, Evan K. Franseen, Richard D. Miller, Thomas V. Weis. Application of deterministic deconvolution of ground-penetrating radar data in a study of carbonate strata[J]. Journal of Applied Geophysics,2004,56(3):213-229.
[73]D. Gomez-Ortiz, T. Martin-Crespo, S. Martin-Velazquez, P. Martinez-Pagan, H. Higueras, M. Manzano. Application of ground penetrating radar (GPR) to delineate clay layers in wetlands. A case study in the Soto Grande and Soto Chico watercourses, Doflana (SW Spain) [J]. Journal of Applied Geophysics,2010,72(2):107-113.
[74]Christopher S. Lowry, Dante Fratta, Mary P. Anderson. Ground penetrating radar and spring formation in a groundwater dominated peat wetland[J]. Journal of Hydrology,2009,373(1-2): 68-79.
[75]王复明,尚向阳,钟燕辉.探地雷达在河流水深探测中的应用研究[J].人民黄河,2008,30(4)：17-18.
[76]郭成超,王复明,谢晓莉等.探地雷达在河道探测中的应用[J].物探与化探,2012,36(5)：887-892.
[77]尚向阳,张汝印,苏茂荣等.探地雷达在黄河枣树沟水下探测中的应用[J].人民黄河,2011,32(3)：10-12.
[78]丁凯,查恩来,周紧东等.应用地质雷达进行水下抛石探测的试验研究[J].水利学报,2002,12(S1)：560-565.
[79]Willianm P, Anderson Jr, David G et al. The effects of Holocene barrier-island evolution on water-table elevations, Hatteras Island, North Carolina, USA[J]. Hydrogeology Journal, 2000(8):390-404.
[80]Guy Serbin, Dani Or. Ground Penetrating radar measurement of crop and surface water content dynamics Slsevier [J]. Remote Sensing Environment,2005,96:119-134.
[81]Mohammad H Makkawi. Integrating GPR and geostatistical techniques to map the spatial extent of a shallow groundwater system[J]. Journal of Geophysics And Engineering,2004(1): 56-62.
[82]邓世坤,王惠濂.探地雷达图像的正演合成与偏移处理[J].地球物理学报,1993,36(4)：528-535.
[83]A A Lestari, A G Ynrovoy, L P Ligthart. Electric Field Integral Equation Model For Transient GPR Antenns [C]. Ninth International Conference on Ground Penetrating Radar,2002.
[84]陈义群,肖柏勋.论探地雷达现状与发展[J].工程地球物理学报,2005,2(2)：149-155.
[85]傅良魁.应用地球物理教程-电法勘探[M].北京：地质出版社,1991.
[86]杨发杰,巨妙兰.高密度电阻率探测方法及其应用[J].矿产与地质,2004,18(4)：356-360.
[87]吕玉增,阮百尧.高密度电法工作中的几个问题研究[J].工程地球物理学报,2005,2(4)：264-269.
[88]严加永,孟贵祥,吕庆田等.高密度电法的进展与展望[J].物探与化探,2012,36(4)：576-577.
[89]Dr. M. H. Loke. Electrical imaging surveys for environmental and engineering studies[M]. Malaysia,1999.
[90]肖宏跃,雷宛.地电学教程[M].北京：地质出版社,2008.
[91]雷宛,肖宏跃,邓一谦.工程与环境物探[M].北京：地质出版社,2006.
[92]陈宁,雷宛,肖宏跃等.边界系数法在高密度电阻率法正反演解释中的应用效果研究[J].成都理工大学学报(自然科学版),2011,38(3)：284-290.
[93]肖宏跃,雷宛,王永庆.实验室高密度电法微测系统的模拟研究[J].物探化探计算技术,2008,30(5)：419-424.
[94]Pedro Martinez-Paganl, Angel Faz Cano.2-D electrical resistivity imaging to assess slurry pond subsoil pollution in the southeastern region of murcia, Spain[J]. Journal of Environmental and Engineering Geophysics,2010,5:29-47.
[95]张大海,王兴泰.二维视电阻率断面的快速最小二乘反演[J].物探化探计算技术,1999,21(1)：2-8.
[96]周熙襄，钟本善.电法勘探数值模拟技术[M].成都：四川科学出版社,1986.
[97]SASAKI Y. Resolution of resistivity tomography inferred from numerical simulation[J]. Geophysical Prospecting,1992,40:453-464.
[98]房纯纲,姚成林,贾永梅.堤坝隐患及渗漏无损检测技术与仪器[M].北京：中国水利水 电出版社,2010.
[99]王爱国,马巍,王大雁.高密度电法不同电极排列方式的探测效果对比[J].工程勘察,2007(1)：72-75.
[100]Beiki, M.. Analytic signals of gravity gradient tensor and their application to estimate source location[J]. Geophysics,2010,75(6):159-174.
[101]Blakely, R. J.. Potential theory in gravity and Magnetic Applications[M]. Cambridge University Press.1995.
[102]Pawlowski, R. S., Hansen, R. O.. Gravity anomaly separation by Wiener filtering[J]. Geophysics,1990,55(5),539-548.
[103]Syberg, F. J. R.. A Fourier method for the regional residual problem of potential fields[J]. Geophysical Prospecting,1972,20(1),47-75.
[104]Wang, W. Y., Pan, Y., Qiu, Z. Y.. A new edge recognition technology based on the normalized vertical derivative of the total horizontal derivative for potential field data[J]. Applied Geophysics,2009,6(3),226-233.
[105]Wang, W. Y, Zhang, G C., and Liang, J. S.. Spatial variation law of vertical derivative zero points for potential field data[J]. Applied Geophysics,2010,7(3),197-209.
[106]肖锋,吴燕冈,孟令顺.位场数据处理中小子域滤波的改进[J].石油地球物理勘探,2010,45(1)：136-139.
[107]Yang, GY. A new technique for potential field data processing:small sub-domain filtering[J]. Oil Geophysical Prospecting,1995,30(1) 240-244.
[108]Jiang, F. Y, Gao, L. K.. Edge enhancement of gravity anomalies and gravity gradient tensors using an improved small sub-domain filtering method[J]. Applied Geophysics,2012,9(2), 119-130.
[109]Mickus, K. L. Hinojosa, J. H.. The complete gravity gradient tensor derived from the vertical component of gravity:a Fourier transform technique[J]. Journal of Applied Geophysics, 2001,46(3),159-174.
[110]王萍.探地雷达检测土壤紧实性的实验研究和信号反演[D].北京：中国矿业大学(北京),2010.
[111]李大心.探地雷达方法与应用[M].北京：地质出版社,1994.
[112]冯德山,戴前伟.用Visual C++开发探地雷达正演模拟软件[J].物探化探计算技术,2004,26(3)：231-234.
[113]冯德山.探地雷达二维时域有限差分正演[D].长沙：中南大学,2003.
[114]Casper D.A. Kung K.J.S. Simulation of ground-penetrating radar wavesin a 2D soil model[J]. Geophysic.1996,61:1034-1049.
[115]葛德彪,闫玉波.电磁波时域有限差分方法[M].西安：西安电子科技大学出版社,2002.
[116]Yee K S. Numerical solution of initial boundary value problems involving Maxwell equations in isotropic media[J]. IEEE Trans:Antennas Propagate,1966,14(3):302-307.
[117]Taylor C D, Lam D H, Shumpert T H. EM pulse scattering in time varying inhomogeneous media[J]. IEEE Trans:Antennas Propagate,1969,17950:585-589.
[118]Taflove A. Review of the formulation and applications of the finite-difference time-domain method for numerical modeling of electromagnetic wave interactions with arbitrary structures[J]. Wave Motion,1998,10:547-582.
[119]Kunz K S, Luebbers R J. The Finite Difference Time Domain Method for Electromagnetic[J]. Boca Raton, FL:CRC Press,1993.
[120]邓艳,于长杰,吴有亮等.水域地震映像法在油气管道穿越工程勘察中的应用[J].中国石油和化工标准与质量,2011,31(8)：104-105.
[121]肖宏跃,雷行建,雷宛.环境物探技术在岩溶勘察中的应用与效果[J].灾害学,2007,22(3)：58-60.
[122]韩连生,范广海,贾建军等.地震映像法在铁矿隐伏空区探测中的应用研究[J].工程地球物理学报,2012,9(4)：469-472.
[123]林文太,卞坚耐.地震映象法在港口码头工程地质勘察中的应用[J].港工技术,2001,1(3)：45-46.
[124]李瑞强.煤矿采空区的地震波场特征研究[D].太原理工大学,2011.
[125]Cerveny V., Molotkovl A. Ray method in seismology [M]. Kari ova University:Prada,1977.
[126]张中杰,何樵登.含裂隙介质中地震波运动学问题的正演模拟[J].石油地球物理勘探,1989,24(3)：290-300.
[127]符力耘,牟永光.弹性波边界元法正演模拟[J].地球物理学报,1994,37(4)：521-529.
[128]闫贫,何樵登.二维横向各向同性介质中的格林函数积分法合成记录[J].物化探计算技术,1995,17(2)：1-6.
[129]Fu L Y.. Numerical study of generalized Lippmann-Schwinger integral equation including surface topography[J]. Geophysics,2003.68:665-671.
[130]Dong Z, McMechan G A..3-D viscoelasitc anisotropic modeling of data from a multi-component, multi-azimuth experiment in northeast Texas[J]. Geophysics,1995,60(4): 1128-1150.
[131]Ramos-martinez J, Ortega A A.. Shear-wave splitting at vertical incidence in media containing intersecting fracture systems[C].68th Ann. Internat. Mtg, Soc. Expl. Geophys., Expanded Abstracts,1998.
[132]Ramos-Martinez J, Ortega AA, McMechan G A.3-D seismic modeling for cracked media: shear-wave splitting at zero offset[J]. Geophysics,2000,65(1):211-221.
[133]侯安宁,何樵登.各向异性介质中弹性波动方程高阶差分及其稳定性研究[J].地球物理学报,1995,38(4)：519-527.
[134]裴正林.三维各向异性介质中弹性波方程交错网格高阶有限差分法数值模拟[J].石油大学学报(自然科学版),2004,25(5)：23-29.
[135]裴正林.三维双相各向异性介质弹性波方程交错网格高阶有限差分法模拟[J].中国石油大学学报(自然科学版),2006,30(2)：608-615.
[136]裴正林.层状各向异性介质中横波分裂和再分裂数值模拟[J].石油地球物理勘探,2006,41(1)：17-25.
[137]Nordstrom J.. Carpenter M H. High-order finite difference methods multidimensional linear problems and curvilinear coordinates[J]. Journal of Computational physics,2001,173(1): 149-174.
[138]Pitarka A.3D Elastic Finite-Difference Modeling of Seismic Motion Using Staggered Grids with Nonuniform Spacing [J]. Bulletin of Seismological Soeiety of America,1999,89(1): 54-68.
[139]Oprsal I, Zahradnik J. Elastic finite-difference method for irregular grids[J]. Geophysics,1999, 64(1):240-250.
[2]钱正英.中国水资源战略研究中几个问题的认识[J].河海大学学报(自然科学版),2001,29(3)：1-7.
[3]程晓陶.新时期大规模的治水活动迫切需要科学理论的指导—论有中国特色的洪水风险管理[J].水利发展研究,2001,1(4)：1-6.
[4]董哲仁.堤防除险加固实用技术[M].北京：中国水利水电出版社,1998.
[5]中华人民共和国水利部和国家统计局.第一次全国水利普查公报[R].2013.
[6]冷元宝,黄建通,张震夏等.堤坝隐患探测技术研究进展[J].地球物理学进展.2003,18(3)：370-379.
[7]徐晓东.安徽省长江2011年度河势演变分析及崩岸监测预警报告[R].安徽长江河道设计研究院,2011.
[8]马洪蛟,李宇,蔡辉等.探地雷达探测抛填厚度研究[J].南京林业大学学报(自然科学版),2004,28(2)：51-53.
[9]葛双成,张孚专,孙伯永等.洞头县杨文围涂东围堤抛石层雷达检测与分析[J].浙江水利科技,2006,(3)：35-36.
[10]孔宪立.工程地质学[M].北京：中国建筑工业出版社,2003.
[11]张玉,朱维申.小湾水电站左岸坝前堆积体在自然状态下稳定性的平面离散元与有限元分析[J].云南水力发电,2000,16(1)：36-39.
[12]唐川,黄楚兴,万哗.云南省丽江大地震及其诱发的崩塌滑坡灾害特征[J].自然灾害学报,1997,6(3)：76-84.
[13]邵景林,祁建忠,高杰.西气东输管道通过长江工程地质勘察水上钻探工艺合理性分析[J].工程勘察,2002(2)：52-53.
[14]高建波.浅谈水上地质钻探方法选择与应用[J].铁道勘察,2012,3：51-54.
[15]季喜平.轻型静力触探技术及其在中小型水利工程中的应用[J].水利水电科技进展,1998,18(4)：4749.
[16]Waite A D.王德石,译.实用声纳工程[M].北京：电子工业出版社,2004.
[17]吴自银,郑玉龙,初凤友,等.海底浅表层信息声探测技术研究现状及发展[J].地球科学进展,2005,20(11)：1210-1217.
[18]Jun Chen, Hong Jie Wan, Xue Wei Wang. A design of underwater acoustic responder based on spread-spectrum communication[J]. Procedia Engineering,2012,29(1):128-132.
[19]Dario Pompili, Tommaso Melodia, Ian F. Akyildiz. Three-dimensional and two-dimensional deployment analysis for underwater acoustic sensor networks[J]. Ad Hoc Networks,2009,7(4): 778-790.
[20]John J. Schultz, Carrie A. Healy, Kenneth Parker, Bim Lowers. Detecting submerged objects: The application of side scan sonar to forensic contexts[J]. Forensic Science International,2013, 231(1-3):306-316.
[21]B. Lafferty, R. Quinn, C. Breen. A side-scan sonar and high-resolution Chirp sub-bottom profile study of the natural and anthropogenic sedimentary record of Lower Lough Erne, northwestern Ireland[J]. Journal of Archaeological Science,2006,33(6):756-766.
22]褚宏宪,白大鹏,史慧杰等.海底抛石声学探测方法-以渤海湾人工岛抛石检测为例[J].海洋地质前沿,2011,27(4)：65-70.
[23]Vladimir Shtivelman. Shallow water seismic surveys for site investigation in the Haifa Port Extension area, Israel[J]. Journal of Applied Geophysics,2001,46(2):143-158.
[24]T. Hussein, A.G. Mula, H. Schneider. Geological and seismic investigations with regard to shallow ground water exploration in Eastern Sudan Republic[J]. Journal of African Earth Sciences (and the Middle East),1989,8(1):75-78.
[25]Shtivelman, V.. Using surface waves for estimating shear wave velocities in the shallow subsurface onshore and offshore area[J]. Eur. J. Environ. Eng. Geophys.,19994,15-35.
[26]Shtivelman, V., Frieslander, U., Zilbennan, R., Amit, R.. Mapping shallow faults at the Evrona playa site using high resolution reflection method[J]. Geophysics,1998,63:1257-1264.
[27]王传雷,曲赞,祁明松等.高精度磁测在长江马当要塞沉船探测中的应用[J].地质科技情报,2000,19(3)：98-102.
[28]喻忠鸿,王传雷,吴文贤.磁法探测炸弹有效深度的物理模拟试验及分析[J].工程地球物理学报,2007,4(2)：118-122.
[29]Dwain K. Butler. Potential fields methods of location of unexploded ordnance[J]. US Army Engineer Research and Development Center, The leading EDGE,2001,20(8):890-895.
[30]吕邦来,王传雷,曲赞.磁测在港航工程勘察中的几个实例[J].工程地球物理学报,2004,1(6)：494498.
[31]Chuanlei Wang. The Magnetic Detection of Sunken Ships in the Madang Section of the Yangtze River[J]. Journal of Environmental & Engineering Geophysics (USA),2006,11(2):123-131.
[32]Zeev Zalevsky, Yuri Bregman, Nizan Salomonski, Hovav Zafrir. Resolution Enhanced Magnetic Sensing System for Wide Coverage Real Time UXO Detection[J]. Journal of Applied Geophysics,2012,84(1):70-76.
[33]SL326-2005,水利水电工程物探规程[S].
[34]赵铁虎,张志珣、许枫.浅水区浅地层剖面测量典型问题分析[J].物探化探计算技术,2002,24(3)：215-216.
[35]滕惠忠,严晓明,李胜全.侧扫声纳图像增强技术[J].海洋测绘,2004,24(2)：47-48.
[36]李金铭,罗延钟.电法勘探新进展[M].北京：地质出版社,1996.
[37]严加永,孟贵祥,吕庆田等.高密度电法的进展与展望[J].物探与化探,2012,36(4)：576-583.
[38]罗延钟,王传雷,董浩斌.府密度电阻率法的电极装置选择[J].地质与勘探,2005,增刊：174-178.
[39]Edwards R N. On the resource evaluation of marine gas hydrate deposits using sea-floor transient electric dipole-dipole method [J]. Geophysics,1997,62:63-74.
[40]王 剑,张赛民,马晓成.水上高密度电法在油气管道工程中应用[J].吉林地质,2012,31(2)：76-79.
[41]吴长盛.北大港水库堤坝裂缝检测与评定技术研究[J].水利水电技术,2001(5)：61-63.
[42]王士鹏.高密度电法在水文地质和工程地质中的应用[J].水文地质工程地质,2000(1)：52-56.
[43]王传雷,董浩斌.科学技术部重点项目、武汉市配套项目“堤坝隐患探测技术研究”(1999—2001年)[R].中国地质大学,2001.
[44]Cheng Q, Tuo XG, Ge B, Li HL, Wang SD. The application of the high-density electrical method to the survey of Chayuangou landslide in Gaochuan, Sichuan province[J]. Geophysical and Geochemical exploration,2012,36:69-72.
[45]Osazuwa IB, Abdullahi NK.2-D Electrical Resistivity and Induced Polarisation Investigation at an Open Solid Waste Dumpsite:Case Study from Kaduna, North Central Nigeria[J]. Journal of Environmental Hydrology.2008,16:1-11.
[46]Al-Tarazi E, Abu Rajab J, Al-Naqa A, El-Waheidi M. Detecting leachate plumes and groundwater pollution at Ruseifa municipal landfill utilizing VLF-EM method[J]. Journal of Applied Geophysics,2008,65:121-131.
[47]Jegede SI, Osazuwa IB, Abdullahi NK. Geoenvironmental Study of groundwater contamination in a dual aquifer environment using earth resistivity imaging[J]. Journal of American Science, 2011,7:367-377.
[48]Fanny Genelle, Colette Sirieix, Joelle Riss, Veronique Naudet. Monitoring landfill cover by electrical resistivity tomography on an experimental site[J]. Engineering Geology,2012, 145-146:18-29.
[49]Luciana Orlando. Some considerations on electrical resistivity imaging for characterization of waterbed sediments[J]. Journal of Applied Geophysics,2013,95:77-89.
[50]A. Kampke. Focused imaging of electrical resistivity data in archaeological prospecting[J]. Journal of Applied Geophysics,1999,41(2-3):215-227.
[51]A. Samouelian, I. Cousin, A. Tabbagh, A. Bruand, G. Richard. Electrical resistivity survey in soil science:a review[J]. Soil and Tillage Research,2005,83(2):173-193.
[52]Dale F. Rucker, Marc T. Levitt, William J. Greenwood. Three-dimensional electrical resistivity model of a nuclear waste disposal site[J]. Journal of Applied Geophysics,2009,69(3-4): 150-164.
[53]王兴泰,李晓芹.电阻率图像重建的佐迪CZo-hdy)反演及其应用效果[J].物探与化探,1996,20(3)：228-233.
[54]张大海,王兴泰.二维视电阻率断面的快速最小二乘反演[J].物探化探计算技术,1999,21(1)：2-8.
[55]Dw groothedlinc. constable s. Occam inversion to generate smooth, two-dimensional models form magnetotelluric data [J]. Geophysics,1990,55:1613-1624.
[56]LOKE M H, BARKER R D. Rapid least-squares inversion of apparent resistivity pseudo sections by a quasi-Newton method[J]. Geophysical Prospecting,1996,44:131-152.
[57]吕邦来,匡翠萍.浅地层剖面探测与地震映像法探测的应用对比和分析[J].港工技术,2009,46(B07)：86-89.
[58]Ricker. The form of laws of propagation of seismic wavelets[J]. Geophysics,1953,18(1): 10-40.
[59]A.T. Dennison. An introduction to synthetic seismogram techniques [J]. Geophysics Prospecting, 1960,8(2):231-241.
[60]P. C. Wuenschel. Seismogram synthesis including multiples and transmission coefficients [J]. Geophysics,1960,25(1):106-129.
[61]付燕.人工地震信号去噪方法研究[D].西安：西北工业大学,2002.
[62]杨农合,徐小林,吴朝俊等.应用瞬变电磁法小线框大电流探测水下构造[J].工程地球物理学报,2012，2,9(3)：279-284.
[63]李大心等.探地雷达探测方法与应用[M].北京：地质出版社,1994.
[64]杜树春.探地雷达及其在环境地质中的应用[J].物探与化探,1996,20(5)：384-392.
[65]李大心等.江河堤防隐患的探地雷达调查[J].中国地质灾害与防治学报,1996,7(1)：21-25.
[66]武可贵,程学标.探地雷达在工程地基检测中的应用研究[J].水运工程,1997(9)：5-8.
[67]邓世坤.探地雷达用于深圳湾软基处理工程中的质量检测[J].岩土工程学报,2000,22(1)：33-35.
[68]雷林源.探地雷达应用中的几个基本问题[J].物探与化探,1998,22(6)：408-414.
[69]王水强,万明浩,谢雄耀等.不同探地雷达测量参数对数据采集效果的影响[J].物探与化探,1999,23(3)：211-215.
[70]Timo Saarenketo, Tom Scullion. Road evaluation with ground penetrating radar[J]. Journal of Applied Geophysics,2000,43(2-4) 119-138.
[71]J.A. Doolittle, B. Jenkinson, D. Hopkins, M. Ulmer, W. Tuttle. Hydropedological investigations with ground-penetrating radar (GPR):Estimating water-table depths and local ground-water flow pattern in areas of coarse-textured soils[J]. Geoderma,2006,131(3-4):317-329.
[72]Jianghai Xia, Evan K. Franseen, Richard D. Miller, Thomas V. Weis. Application of deterministic deconvolution of ground-penetrating radar data in a study of carbonate strata[J]. Journal of Applied Geophysics,2004,56(3):213-229.
[73]D. Gomez-Ortiz, T. Martin-Crespo, S. Martin-Velazquez, P. Martinez-Pagan, H. Higueras, M. Manzano. Application of ground penetrating radar (GPR) to delineate clay layers in wetlands. A case study in the Soto Grande and Soto Chico watercourses, Doflana (SW Spain) [J]. Journal of Applied Geophysics,2010,72(2):107-113.
[74]Christopher S. Lowry, Dante Fratta, Mary P. Anderson. Ground penetrating radar and spring formation in a groundwater dominated peat wetland[J]. Journal of Hydrology,2009,373(1-2): 68-79.
[75]王复明,尚向阳,钟燕辉.探地雷达在河流水深探测中的应用研究[J].人民黄河,2008,30(4)：17-18.
[76]郭成超,王复明,谢晓莉等.探地雷达在河道探测中的应用[J].物探与化探,2012,36(5)：887-892.
[77]尚向阳,张汝印,苏茂荣等.探地雷达在黄河枣树沟水下探测中的应用[J].人民黄河,2011,32(3)：10-12.
[78]丁凯,查恩来,周紧东等.应用地质雷达进行水下抛石探测的试验研究[J].水利学报,2002,12(S1)：560-565.
[79]Willianm P, Anderson Jr, David G et al. The effects of Holocene barrier-island evolution on water-table elevations, Hatteras Island, North Carolina, USA[J]. Hydrogeology Journal, 2000(8):390-404.
[80]Guy Serbin, Dani Or. Ground Penetrating radar measurement of crop and surface water content dynamics Slsevier [J]. Remote Sensing Environment,2005,96:119-134.
[81]Mohammad H Makkawi. Integrating GPR and geostatistical techniques to map the spatial extent of a shallow groundwater system[J]. Journal of Geophysics And Engineering,2004(1): 56-62.
[82]邓世坤,王惠濂.探地雷达图像的正演合成与偏移处理[J].地球物理学报,1993,36(4)：528-535.
[83]A A Lestari, A G Ynrovoy, L P Ligthart. Electric Field Integral Equation Model For Transient GPR Antenns [C]. Ninth International Conference on Ground Penetrating Radar,2002.
[84]陈义群,肖柏勋.论探地雷达现状与发展[J].工程地球物理学报,2005,2(2)：149-155.
[85]傅良魁.应用地球物理教程-电法勘探[M].北京：地质出版社,1991.
[86]杨发杰,巨妙兰.高密度电阻率探测方法及其应用[J].矿产与地质,2004,18(4)：356-360.
[87]吕玉增,阮百尧.高密度电法工作中的几个问题研究[J].工程地球物理学报,2005,2(4)：264-269.
[88]严加永,孟贵祥,吕庆田等.高密度电法的进展与展望[J].物探与化探,2012,36(4)：576-577.
[89]Dr. M. H. Loke. Electrical imaging surveys for environmental and engineering studies[M]. Malaysia,1999.
[90]肖宏跃,雷宛.地电学教程[M].北京：地质出版社,2008.
[91]雷宛,肖宏跃,邓一谦.工程与环境物探[M].北京：地质出版社,2006.
[92]陈宁,雷宛,肖宏跃等.边界系数法在高密度电阻率法正反演解释中的应用效果研究[J].成都理工大学学报(自然科学版),2011,38(3)：284-290.
[93]肖宏跃,雷宛,王永庆.实验室高密度电法微测系统的模拟研究[J].物探化探计算技术,2008,30(5)：419-424.
[94]Pedro Martinez-Paganl, Angel Faz Cano.2-D electrical resistivity imaging to assess slurry pond subsoil pollution in the southeastern region of murcia, Spain[J]. Journal of Environmental and Engineering Geophysics,2010,5:29-47.
[95]张大海,王兴泰.二维视电阻率断面的快速最小二乘反演[J].物探化探计算技术,1999,21(1)：2-8.
[96]周熙襄，钟本善.电法勘探数值模拟技术[M].成都：四川科学出版社,1986.
[97]SASAKI Y. Resolution of resistivity tomography inferred from numerical simulation[J]. Geophysical Prospecting,1992,40:453-464.
[98]房纯纲,姚成林,贾永梅.堤坝隐患及渗漏无损检测技术与仪器[M].北京：中国水利水 电出版社,2010.
[99]王爱国,马巍,王大雁.高密度电法不同电极排列方式的探测效果对比[J].工程勘察,2007(1)：72-75.
[100]Beiki, M.. Analytic signals of gravity gradient tensor and their application to estimate source location[J]. Geophysics,2010,75(6):159-174.
[101]Blakely, R. J.. Potential theory in gravity and Magnetic Applications[M]. Cambridge University Press.1995.
[102]Pawlowski, R. S., Hansen, R. O.. Gravity anomaly separation by Wiener filtering[J]. Geophysics,1990,55(5),539-548.
[103]Syberg, F. J. R.. A Fourier method for the regional residual problem of potential fields[J]. Geophysical Prospecting,1972,20(1),47-75.
[104]Wang, W. Y., Pan, Y., Qiu, Z. Y.. A new edge recognition technology based on the normalized vertical derivative of the total horizontal derivative for potential field data[J]. Applied Geophysics,2009,6(3),226-233.
[105]Wang, W. Y, Zhang, G C., and Liang, J. S.. Spatial variation law of vertical derivative zero points for potential field data[J]. Applied Geophysics,2010,7(3),197-209.
[106]肖锋,吴燕冈,孟令顺.位场数据处理中小子域滤波的改进[J].石油地球物理勘探,2010,45(1)：136-139.
[107]Yang, GY. A new technique for potential field data processing:small sub-domain filtering[J]. Oil Geophysical Prospecting,1995,30(1) 240-244.
[108]Jiang, F. Y, Gao, L. K.. Edge enhancement of gravity anomalies and gravity gradient tensors using an improved small sub-domain filtering method[J]. Applied Geophysics,2012,9(2), 119-130.
[109]Mickus, K. L. Hinojosa, J. H.. The complete gravity gradient tensor derived from the vertical component of gravity:a Fourier transform technique[J]. Journal of Applied Geophysics, 2001,46(3),159-174.
[110]王萍.探地雷达检测土壤紧实性的实验研究和信号反演[D].北京：中国矿业大学(北京),2010.
[111]李大心.探地雷达方法与应用[M].北京：地质出版社,1994.
[112]冯德山,戴前伟.用Visual C++开发探地雷达正演模拟软件[J].物探化探计算技术,2004,26(3)：231-234.
[113]冯德山.探地雷达二维时域有限差分正演[D].长沙：中南大学,2003.
[114]Casper D.A. Kung K.J.S. Simulation of ground-penetrating radar wavesin a 2D soil model[J]. Geophysic.1996,61:1034-1049.
[115]葛德彪,闫玉波.电磁波时域有限差分方法[M].西安：西安电子科技大学出版社,2002.
[116]Yee K S. Numerical solution of initial boundary value problems involving Maxwell equations in isotropic media[J]. IEEE Trans:Antennas Propagate,1966,14(3):302-307.
[117]Taylor C D, Lam D H, Shumpert T H. EM pulse scattering in time varying inhomogeneous media[J]. IEEE Trans:Antennas Propagate,1969,17950:585-589.
[118]Taflove A. Review of the formulation and applications of the finite-difference time-domain method for numerical modeling of electromagnetic wave interactions with arbitrary structures[J]. Wave Motion,1998,10:547-582.
[119]Kunz K S, Luebbers R J. The Finite Difference Time Domain Method for Electromagnetic[J]. Boca Raton, FL:CRC Press,1993.
[120]邓艳,于长杰,吴有亮等.水域地震映像法在油气管道穿越工程勘察中的应用[J].中国石油和化工标准与质量,2011,31(8)：104-105.
[121]肖宏跃,雷行建,雷宛.环境物探技术在岩溶勘察中的应用与效果[J].灾害学,2007,22(3)：58-60.
[122]韩连生,范广海,贾建军等.地震映像法在铁矿隐伏空区探测中的应用研究[J].工程地球物理学报,2012,9(4)：469-472.
[123]林文太,卞坚耐.地震映象法在港口码头工程地质勘察中的应用[J].港工技术,2001,1(3)：45-46.
[124]李瑞强.煤矿采空区的地震波场特征研究[D].太原理工大学,2011.
[125]Cerveny V., Molotkovl A. Ray method in seismology [M]. Kari ova University:Prada,1977.
[126]张中杰,何樵登.含裂隙介质中地震波运动学问题的正演模拟[J].石油地球物理勘探,1989,24(3)：290-300.
[127]符力耘,牟永光.弹性波边界元法正演模拟[J].地球物理学报,1994,37(4)：521-529.
[128]闫贫,何樵登.二维横向各向同性介质中的格林函数积分法合成记录[J].物化探计算技术,1995,17(2)：1-6.
[129]Fu L Y.. Numerical study of generalized Lippmann-Schwinger integral equation including surface topography[J]. Geophysics,2003.68:665-671.
[130]Dong Z, McMechan G A..3-D viscoelasitc anisotropic modeling of data from a multi-component, multi-azimuth experiment in northeast Texas[J]. Geophysics,1995,60(4): 1128-1150.
[131]Ramos-martinez J, Ortega A A.. Shear-wave splitting at vertical incidence in media containing intersecting fracture systems[C].68th Ann. Internat. Mtg, Soc. Expl. Geophys., Expanded Abstracts,1998.
[132]Ramos-Martinez J, Ortega AA, McMechan G A.3-D seismic modeling for cracked media: shear-wave splitting at zero offset[J]. Geophysics,2000,65(1):211-221.
[133]侯安宁,何樵登.各向异性介质中弹性波动方程高阶差分及其稳定性研究[J].地球物理学报,1995,38(4)：519-527.
[134]裴正林.三维各向异性介质中弹性波方程交错网格高阶有限差分法数值模拟[J].石油大学学报(自然科学版),2004,25(5)：23-29.
[135]裴正林.三维双相各向异性介质弹性波方程交错网格高阶有限差分法模拟[J].中国石油大学学报(自然科学版),2006,30(2)：608-615.
[136]裴正林.层状各向异性介质中横波分裂和再分裂数值模拟[J].石油地球物理勘探,2006,41(1)：17-25.
[137]Nordstrom J.. Carpenter M H. High-order finite difference methods multidimensional linear problems and curvilinear coordinates[J]. Journal of Computational physics,2001,173(1): 149-174.
[138]Pitarka A.3D Elastic Finite-Difference Modeling of Seismic Motion Using Staggered Grids with Nonuniform Spacing [J]. Bulletin of Seismological Soeiety of America,1999,89(1): 54-68.
[139]Oprsal I, Zahradnik J. Elastic finite-difference method for irregular grids[J]. Geophysics,1999, 64(1):240-250.