可控源音频大地电磁测深法应用实例_图文

1、第28卷第4期物探化探计算技术2006年11月 收稿日期:2005-04-29文章编号:10011749(200604033705可控源音频大地电磁测深法应用实例黄力军,陆桂福,刘瑞德(中国地质科学院地球物理地球化学勘查研究所,廊坊065000摘要:目前对电法的探测深度要求越来越高,常规电法有时难以满足实际工作需要。与常规电法相比,可控源音频大地电磁测深法具有探测深度大、设备相对轻便和横向分辨率高等特点,已广泛应用于地热资源、水文地质、有色金属矿产和工程地质勘查工作之中,并取得了令人十分满意的实际应用效果。这里介绍的是一些具有代表性可控源音频大地电磁测深应用实例和应用成果,这些实例基本函概了目

2、前电法的主要应用范畴,常规电法难以胜任的地质调查工作。实例中包括断陷盆地型地热资源和深部基岩裂隙水调查,危机矿山深部找矿和复杂地形条件下铁路隧道工程勘查。关键词:电磁测深;地热资源;水文地质;有色金属矿产;工程地质;调查中图分类号:P 63113+25文献标识码:A0前言可控源音频大地电磁测深是在音频大地电磁测深(AMT 基础上,发展起来的一种人工源频率测深方法。上世纪五十年代,在卡尼亚(L Cagn 2iard 论文1的基础上,发展形成了基于观测天然场大地电场和磁场正交分量,计算视电阻率的大地电磁测深方法。在音频(n ×10-1Hz n ×103Hz 范围内,大地电磁场相

3、对较弱,人文干扰较大。为了克服上述困难,七十年代初,D W Strang way 教授和他的学生M A Goldstein 2,3提出沿用AMT 的测量方式,观测人工供电产生的音频电磁场。由于所观测的电磁场频率、场强和方向均由人工控制,其观测的方式与AMT 相同,所以称其为可控源音频大地电磁测深(简称CS AMT 法。CS AMT 法可以采用磁性源或电性源二种人工场源,目前主要应用电性源可控源音频大地电磁测深法。可控源音频大地电磁测深是以有限长接地导线为场源,在距场源中心一定距离处同时观测电、磁场参数的一种电磁测深方法。目前应用大多采用赤道偶极装置进行标量测量,同时观测与场源平行的电场水平分量

4、E x 和与场源正交的磁场水平分量H y (测量装置布置见图1。利用电场振幅E x和磁场振幅H y 计算卡尼亚阻抗电阻率s ,电场相位E p 和磁场H p 计算卡尼亚阻抗相位s 。阻抗电阻率和阻抗相位联合反演计算反演电阻率参数,利用反演电阻率进行地质推断解释。目前,可控源音频大地电磁测深反演主要采用一维电阻率圆滑反演方法4。图1标量CS AMT 测量装置示意图Fig .1Scalar CS AM T survey array十几年来,我们一直从事可控源音频大地电磁测深应用研究和生产工作,在各个应用领域中均做过很多工作。在本文中,作者以典型实例介绍了可控源音频大地电磁测深方法的实际应用及应用效果

5、。1地热资源调查作为绿色能源,地热资源的开发利用日益受到人们的广泛重视。一般情况下,地热资源埋藏较深,开采风险很大。为了提高效率,降低成本,减少投资风险,在开发地热资源前,大多进行地质调查,地球物理勘查就是地热资源调查的重要手段之一。地下热水能够使热储层位电阻率降低,低阻异常是识别地热资源的重要标志,所以在以往地热资源调查中,直流电阻率法一直占主导地位。但直流电阻率法的工作效率低,勘探深度相对较浅,在使用中受较大限制。目前地热资源开发的深度越来越大,大多开采深度已超过2000m。随着深度加大,地表观测到的由地下热水引起的电阻率差异越来越小,以至难以观测到由地热变化引起的电阻率异常,所以,深部地

6、热资源调查的主要任务是勘查热储地层及地质构造分布情况。近几年,我们采用可控源音频大地电磁测深法,分别在北京、辽宁、河北、江西、江苏和浙江等地的深部地热资源调查中,均取得了较好的地质效果。断陷盆地型热水埋深相对较大,以中低温为主,分布范围较广,是目前利用最广的地热资源,在北京已经得到广泛地开发和利用。北京市及周边地区蓟县系雾迷山组白云岩厚度较大,埋深一般0mn×1000m,岩溶裂隙较为发育,热导率值较高,是北京市主要热储地层。该区地热勘查主要任务是了解蓟县系白云岩埋藏深度、上覆地层和构造分布情况。以北京市远郊某区为勘查实例,该测区大部为第四系所覆盖,第四系(Q由粘土、粉砂和砂砾石组成;

7、第四系下部为白垩系(K砂岩和侏罗系(J砂岩,侏罗系内穿插有火山岩;侏罗系下部为蓟县系(Jx,岩性主要为白云岩。在测区内地层间电性差异较为明显,第四系电阻率较低,约20m左右;白垩系电阻率20m 50m;侏罗系砂岩电阻率为50m 100m;蓟县系白云岩电阻率大于500m。各时代地层岩性之间电阻率存在明显差异,为开展电法勘探工作提供了充分的地球物理前提。图2是该测区某地可控源音频大地电磁测深反演电阻率和地质解释综合断面图。由图2可见,断面内除x=1500m附近出现断裂破碎带外,从剖面起点到端点断面电阻率逐渐升高,符合实际地层呈单斜变化趋势。断面内0m-500m为第四系,河滩局部出现高阻外,电阻率均

8、很低;除第四系和白垩系间电阻率差异较小,层间电阻率呈渐变关系,除两者间界面不太明显外,白垩系侏罗系蓟县系界面非常清楚。实际验证(地热井结果, 1296m见蓟县系白云岩,利用CS AMT反演电阻率推断蓟县系白云岩顶界面深度与实际地热井验证深度结果仅差几十米(见图2。该地热井已开发出温度71°C,出水量2000m3/d以上的地下热水。由于地热开发井终孔于断裂破碎内,是北京市目前温度和出水量最高的地热开发井 。图2北京市某地CS AM T地热勘查综合断面图Fig.2The comp rehensive cr oss secti on ofCS AM T exp l orati on,Bei

9、jing2水文地质调查直流电阻率测深法是浅部基岩裂隙水调查的主要手段。但其勘探深度相对较小,施工过程相对较为繁琐。当地表障碍物较多时,施工非常困难,尤其随着深度加大,直流电阻率测深法更加难以满足实际工作需要。下面是北京郊区水文地质调查实例。目前,北京应用的地下水主要分二个方面:一个是灰岩或白云岩岩溶裂隙水;另一个是断裂构造水。由于奥陶系灰岩岩溶破碎相对发育,与其它相应岩层相比埋深较浅,因而是目前北京市地下水开发的主要来源。勘探区为第四系覆盖,附近地表出露有石炭二迭系(C-P、奥陶系(O和寒武系(。其中,石炭二迭系主要为砂岩、页岩和泥岩,电阻率100m左右;奥陶系主要为灰岩,电阻率大于500m;

10、寒武系主要为砂岩、页岩和灰岩,电阻率100m300m。剖面附近出露有花岗岩体。833物探化探计算技术28卷可控源音频大地电磁测深的勘查目的是,利用剖面测量方法确定实测剖面内地层、花岗岩体和地质构造的分布情况,勘探深度2000m 。图3是区内实测剖面CS AMT 勘查综合断面图。由图3可见,除表层外,反演电阻率曲线可大致分为四个值区,其中从剖面起点到x =1300m 处为高值区,平均电阻率大于1000m ,推断为花岗岩体产生;从x =1300m 处到剖面末端的区域内,由上至下又分为三个值区,反演电阻率曲线为K 型,根据反演电阻率曲线及对应岩石电性,推断由上至下分别为石炭二迭系、奥陶系和寒武系(地

11、质推断结果见图3 。图3北京市某乡CS AM T 水文地质勘查综合断面图Fig .3The comp rehensive cr oss secti on of CS AMThydr ogeol ogical exp l orati on,Beijing图3内水文井为地下水开发钻孔,钻孔揭露上部见有砂页岩,下部为花岗岩和蚀变灰岩,地层深度与解释结果基本一致,终孔于断裂破碎带中(井深1000m ,出水量1000m 3/d 以上。3危机矿山深部找矿电法在金属矿产物探勘查工作中一直占有重要地位,随着找矿工作的进一步深入,浅表矿日益减少,寻找深部隐伏矿是今后找矿的主要目标。另一方面,找矿工作重点已经逐渐

12、转向我国西部地区,西部地区地表和地形条件更加复杂,一些常规电磁法已无法满足深部及复杂地形条件工作需要。可控源音频大地电磁测深方法采用定源观测,具有轻便、快速和勘探深度大等优点,其在复杂地形条件和深部找矿工作中,确定成矿有利地段和隐伏矿定位预测等方面起到了其它方法无法替代的重要作用。这里的定位预测实际上是指在详查找矿阶段,采用以电测深为主的方法手段对矿(化体进行二维或三维空间定位。根据其产状、规模和延伸(深情况,结合已知地质和其它资料进行合理地质解释的一种找矿方法技术。在实际工作中,我们一般采用中梯装置的时间域激发极化方法,寻找勘查区内成矿有利地段,使用可控源音频大地电磁测深法,在成矿有利地段对

13、矿(化体进行定位预测,并结合地质、地球化学和其它资料进行综合研究的找矿方法。事实表明,上述方法在实际找矿工作中可取得明显的地质效果。下面是甘肃白银有色金属公司小铁山矿深部找矿实例。白银有色金属公司是五十年代初兴建的我国最大的铜有色金属生产基地,近十年来,后备矿源出现危机。1992年我所参加白银矿区地质和物化探综合找矿研究工作,经过充分研究前人在中、大比例尺地质和物化探资料基础上,在小铁山矿区选择成矿有利地段进行100m ×20m 网度的可控源音频大地电磁面积性测量及重点剖面的岩石地球化学、重力测量工作。小铁山矿床是一个隐伏的海相火山岩Pb 、Zn 、Cu 多金属矿床。地表仅有一铁帽露

14、头,矿体全部埋藏于深部。区内出露的地层为白银群变质火山岩,可分为上、中、下三层,即上部石英角斑岩夹角斑凝灰岩;中部为硅质千枚岩夹石英角斑凝灰岩层;下部为石英角斑凝灰岩层。小铁山含矿带富存于一背斜轴部的石英角斑凝灰质碎裂岩中,矿化带上盘为硅质千枚岩,下盘为石英钠长斑岩体,已知矿的性质和分布与断裂构造关系密切。资料处理后分别在相距100m 的12线和13线可控源音频大地电磁测深反演电阻率断面内,出现相互对应的二处低阻异常(1和2异常,其中1(1-1和1-2异常异常对应已知矿体,2异常位于深部探矿工程的空白地段的175号点附近,12线的2异常位于海拔1100m 1800m 处,为单一低阻异常;13线

15、的2异常由二部分组成(21和22异常,其中21异常位于海拔1400m 1700m 处,22异常位于21异常下方海拔850m 1400m 处(见下页图4。根据地质资料和综合物化探异常,白银有色金属公司地质队1994年采用钻探验证(钻孔编号bzk112线的2异常。钻孔揭露484135m 终孔(724102m 之间均处于蚀变矿化岩石和近矿围岩内,其中铜矿化有向下增高趋势(Cu 平均含量0104%,最高含量达01095%,高于矿体近矿围岩含量(01036%和本区背景值(010022%。由于工程故障,钻孔未到设计深度(800m 。9334期黄力军等:可控源音频大地电磁测深法应用实例图4小铁山矿区13线C

16、S AM T反演电阻率异常和钻孔验证结果Fig.4CS AM T exp l orati on comp rehensive crs oss secti on,Xiaotieshan M ine在距12线100m处的13线布置22异常验证孔(钻孔编号13B2,见厚1717m,综合品位(Cu、Pb和Zn达17%的工业矿体;由于矿体埋深较大,后来加拿大明科公司在1547m的巷道内连续施工三个定向孔,每个钻孔均见工业矿体5,矿体的总金属量为450000t。这是目前我国矿山深部物化探找矿成功的实例。在小铁山深部找矿工作中CS AMT方法起到了关键的,其它方法难以做到的定位预测作用。4工程地质勘查地震勘

17、探一直是工程地质勘查的主要手段,在地震勘探难以施工地区,多采用直流电阻率测深或瞬变电磁测深方法等常规电法。在通常情况下,常规电法勘探深度较小,同时受到地表条件影响,在复杂地形条件下,直流电阻率测深和瞬变电磁测深方法等常规电法很难完成相应的勘查工作。以CS AMT法铁路隧道勘查工作为实例。寺铺尖隧道位于太行山深处,是朔黄铁路上最长的一个隧道。隧道全长6400m,隧道通过的地段地表条件非常复杂,相对高差超过1000m,大多地段行人难以通行。勘查目的是探测隧道通过的地段地层、可能富水地段和地质构造分布等情况,探测深度为1000m。在实际工作中,我们采用隧道上部和隧道二侧相互结合的测量方法。即尽量在设

18、计隧道上部测量,当设计隧道上部难以通过时,采用在隧道二侧布设测线,用二侧剖面测量结果来推断隧道通过地段地质情况的勘查方法。寺铺尖隧道区除沟谷中分布有第四系洪积层外,大部基岩裸露;隧道出口段出露有太古界五台山群混合花岗岩;隧道进口段主要出露有古生界中、上寒武系(2+3灰岩,其中,中寒武系下部为页岩;山顶处出露奥陶系下统(O1灰岩。出口处半坡见有下寒武系(1黄紫色页岩直接角度不正合于太古界地层之上。统计结果表明:奥陶系下统电阻率800m1200m;寒武系中、上统电阻率500m800m;寒武系下、中统电阻率500m800m;太古界五台山群电阻率大于1500m。图5(a(见下页是其中一条旁测剖面的CS

19、 AMT反演电阻率断面图。我们根据CS AMT反演电阻率横向变化(电阻率横向间断推断出剖面通过地段断层展布情况,这里推断出五条断层(分别将其编号为F2、F3、F4、F3+1和F5断层;结合岩石物性及电阻率随着深度变化推断出地层分布情况。图5(b是根据二条旁测剖面测量结果推断的设计隧道地质推断解释断面图(局部。施工单位利用我们预测结果施工,实际结果证明我们推断的地层地质构造均较准确6。这是目前CS AMT法在复杂地形条件下工程地质勘查工作中比较成功的应用实例。043物探化探计算技术28卷图5朔黄铁路寺铺尖隧道CS AMT勘查综合断面图Fig.5CS AMT exp l orati on comp

20、 rehensive crs oss secti on,Si pujian Tunnel5结语作为一种人工源电磁法,在可控源音频大地电磁测深方法测量中不可避免地会存在场源效应。所以在实际工作中,应根据工作任务,结合实际地质和地层电性分布情况,合理地选择场源方向、收发距和采样频率范围。在保证最佳信噪比情况下,要尽量使所采集的数据处在远区。在数据处理中,要注意识别场源效应和静态效应等影响引起的数据畸变,慎重处理这些效应的影响。无论资源开发还是工程施工进行前期调查均可以提高效率、降低成本、减少风险。实际应用结果表明,可控源音频大地电磁测深法能够适应复杂地形条件,满足探测深度较大的地质勘查工作要求。需

21、要说明的是,随着勘探深度加大,其勘查手段必须提高观测和解释精度,以加大勘探深度,提高方法的实际应用效果。另一方面,引起地下电阻率变化的因素很多,利用可控源音频大地电磁测探法勘查,需要结合已知地质、地球物理特征和其它有关资料进行综合研究,合理地进行地质解释。上述实例表明,在一些地热资源、水文地质、有色金属矿产资源和工程地质等地质勘查工作中,可控源音频大地电磁测深法能够起到其它方法无法起到的作用。参考文献:1C AG N I A RD L.Basic theory of the magnet otelluric meth2od of geophysical p r os pectingJ.Geop

22、hysics,1953,18:605.2STRANG WAY D W,S W I FT C M,HOLMER R C.Theapp licati on of audi o2frequency magnet otellurics(AM Tt o m ineral exp l orati onJ.Geophysics,1973,38:1159.3G OLDSTE I N M A,STRANG WAY D W.Audi o2frequen2cy magnet otellurics with a grangded electric di poles ourceJ.Geophysics,1975,40:

23、1159.4C ONST ABLE S C,P ARKER R L,C ONST ABLE C G.Occa ms inversi on:A p ractical algorith m for generatings mooth models fr om electr omagnetic s ounding dataJ.Geophysics,1987,52(3:289.5胡达骧,徐叶兵,陈群,译.明科矿业及金属公司在白银厂地区找矿新进展J.地质与勘探,1999,35(2:1.6杜彦军,王皓,李玮.可控源音频大地电磁测深法在寺铺尖隧道勘察中的应用J.铁路地质与地基,1998,37(1:19.作者

24、简介:黄力军(1954-,男,1982年毕业于长春地质学院物探系,1992年获中国地质大学(北京硕士学位,教授级高级工程师,主要从事电磁法勘查研究和生产工作。1434期黄力军等:可控源音频大地电磁测深法应用实例 Nov. 2006 COMPUTI G TECHN I UES FOR GEOPHYSICAL AND GEOCHEM ICAL EXPLORATI N N Q O 3 WANG Zhi2gang , HE Zhan 2xiang , W E I W en 2 3 ( 1. China University of Geosciences, Beijing bo 100083, Chin

25、a; 2. Departm ent of East Geophysical Company, Zhuozhou Heibei Province 072751, Chi2 na; 3. Geo - detection laboratory,M inistry of Educa2 tion, B eijing 100083, China . COM PU T IN G TECH 2 1, 2, 3 2 study the results of p rocession. Key words: high density resistivity m ethod; defect error; accide

26、ntal error; p rep rocess APPL ICAT I N O F PO TENT IAL F IELD NO R2 O M AL IZED DOW NW ARD EXTRAPOLAT I N IN O SP PRO BAB I ITY TOMO GRAPHY L RUAN Shuai, WANG Xu 2ben, GAO Yong2cai, et al . ( College of Information Engineering, University of Science and Technology of Chengdu, Chengdu Si2 chuan 61005

27、9, China . COM PU T IN G TECHN IQU ES RA T ION , 2006, 28 ( 4 : 332 N IQU ES FOR GEO PHYS ICAL AND GEOCHEM ICAL EX PLORA T ION , 2006, 28 ( 4 : 322 THE AR ITHM ET IC O F THE DATA PREPRO 2 CESS ING IN THE H IGH D ENS ITY RES IS IT IV2 ITY M ETHOD L I Yi2nian, LUO Yan 2zhong ( Institute of Geophysics

28、and Geom atics, China University of Geosciences, W uhan 430074, China . COM PU T IN G TECH 2 N IQU ES FOR GEO PHYS ICAL AND GEOCHEM ICAL EX PLORA T ION , 2006, 28 ( 4 : 328 The high density resisitivity m ethod through elec2 trode power supp ly, surveys electric current and volt2 age, gathers appare

29、nt resistivity data. B ecause of fault in operation and the p recision of instrum ent, the defect error and accidental error are involved in the p rocession of data collection, the defect error that are brought by m istakes of peop le and influence of outside is large, the contingency error that are

30、 inevitably and random icity is little. The m ission of computer is to lim it the defect error, supp ress the accidental error and m ini izes the influence of these errors to the in2 m version result This Paper analyzes the reason of error . s appearance and the characters of the m ethod of re2 sist

31、ivity p rocessing, compares the manual work w ith computer calculate p rogram , gives out several data perp rocessing p rogram s which base on different arith2 metics for the p rocession of error - lim itation, and The paper discusses the calculation method of field induced by vertical line sources

32、by segments and field stacking The method realizes 3D modeling of . field induced by DC in cased well The paper analy2 . zes comp ress storage of large coefficient m atrix which realizes 1D chain comp ress storage of coefficient ma2 trix by combination of array and structure. During so2 lution of su

33、pper - large sparse system of linear equa2 tions, incomp lete Cholesky factorization of stabilized biconjugate gradient ( I CBG method is app lied. Comparing of the analytic solutions and numerical ones of line sources in homogeneous half - space indi2 cates that the m ethod is accurate and reliable

34、. Key words: DC borehole 2ground m ethod; comp ress storage; 3D finite difference method; biconjugate gra2 dient m ethod; incomp lete cholesky factorization CO NTROLL ED SO URE AUD I - FREQUENCY O M AGNETO TELL UR ICS HUANG L i2jun, LU Gui2fu, L I Rui2de ( Institute of U Geophysical and Geochem ical

35、 Exp loration, CAGS Langfang 065000, China . COM PU T IN G TECH 2 N IQU ES FOR GEO PHYS ICAL AND GEOCHEM ICAL EX PLORA T ION , 2006, 28 ( 4 : 337 B ased on the normalized downward extrapolation theory, we have i p roved vertical p recision of tradi2 m tional surface self - potential ( SP tomography

36、meth2 od, of which p rocessing only uses the surface ob2 served voltages After that, external lim itations that . cause the shortages in rep resenting the whole under2 ground information can be volum inously avoided. During model analysis, the new tomography schem e can p rovide more exact results t

37、o fit the true original model. Towards the p ractical data exp lanations, it ac2 curately gives both horizontal abnor ity corresponding m w ith p ractical survey and vertical infor ation of whole m underground space. Key words: Self - potential, Surface SP Tomo2 graphy, Normalized Downward Extrapola

38、tion W ith the dep th increasing of the underground ob2 jects to be p rospected, the conventional electric meth2 ods are difficult to m eet the p rospecting requirem ents . Compared w ith the conventional electric m ethods, the Controlled source audio - frequency magnetotellurics m ethod is featured

39、 of large p rospecting dep th, porta2 ble equipm ent and high lateral resolution. It has been w idely used in geothermal exp loration, hydro - geo2 logical exp loration, non - ferrous m ineral exp loration and engineering geological investigation, and satisfac2 tory field results have been obtained

40、in these area. However, few papers concerning the app lication of the m ethod have been published on publicly circula2 ted journals until now, and even fewer are the papers FOR GEO PHYS ICAL AND GEOCHEM ICAL EX PLO 2 4 ABSTRACTSVol 28 No. 4 . THE I PRO VEM ENT AND D EVELO PM ENT M O F THE M UL T I -

41、 FUNCT I N ASTAT I SO FT2 O C W ARE L I J ian 2xin, WANG Gui2cai, L I Chun 2 ing U U m ( School of Info - Physics Engineering, Center South University , Changsha 410083, China . COM PU T IN G TECHN IQU ES FOR GEO PHYS ICAL AND GEO 2 CHEM ICAL EX PLORA T ION , 2006, 28 ( 4 : 346 concerning the system

42、atic app lication of the m ethod. This paper describes some rep resentative field cases and results of the Controlled source audio - frequency magnetotellurics method. These cases basically cover the main app lication areas of the electric method and the works to which the conventional electric meth

43、od can not competent, which include fault basin geother2 mal resource and deep bedrock fracture water p ros2 pecting, deep m ineral exp loration in the near dep lete m ine, and geotechnical investigation over the comp lex terrain for railw ay tunnels . Key words: eletromagnetic soundings; geother al

44、 m resources; hydrogeology; non 2ferrous m ine; engineer2 ing geology; investigation THE PARALL EL COM PUTAT I N O F LA 2 O GRANGE M UL T IPL IER SEARCH ING IN M AG2 NETO TELL UR IC OCCAM INVERS I N O 1, 2 ( 11China University of Geoscience,W uhan L I Yu U Hubei 430074, China; 2. Guilin Institute of

45、 Technolo2 gy, Guilin Guangxi 541004, Chian . COM PU T IN G Occam inversion has been w idely used in magne2 totelluric data p rocessing for its convergence stability and independent of startup model, but the Lagrange multip lier ( searching result in too m any forward calculations and reduce the inv

46、ersion speed greatly . Based on discussing the calculation method of , the paper put forward using parallel computation to solve this p roblem and realized scanning and searching parallel computation w ith m aster - slave p rogramm ing mode. To reduce the comm unication and achieve load balancing, c

47、omputation tasks were arranged in group s and each grope were dispatched to different computation node. The examp le show s that parallel computation can speed up the p rocessing obviously . Key words: M T; occam inversion; one di ensional m search; parallel computation; COW The ASTATIC is an accord

48、ing p rocessing p ro2 gram of GDP system , which is used for the static TECHN IQU ES FOR GEO PHYS ICAL AND GEO 2 CHEM ICAL EX PLORA T ION , 2006, 28 ( 4 : 342 THE STUDY O F S IGNAL EST I AT I N FROM M O GRO UND PENETRATE RADAR DATA IN W AVEL ET DOM A IN INTERFERED BY W H ITE NO ISE Q I N ing, CHANG

49、Yan 2jun, CAO Zhong2lin ( Institu2 U te of Geophysics and Geom atics, China University of Geosciences, W uhan, 430074, China . COM PU T IN G Considering that the observed signals in ground penetrate radar data are som etim es polluted by noise, we introduce, in this paper, the app lication of smoot2

50、 hing factor in wavelet domain to esti ate available m signal from GPR data interfered by white noise, ac2 cording to wavelet analysis and fractal theory both havng identical p roperties in m ulti - scale analysis and self sim ilarity First, the signal w ith white noise . is transformed via the wavelet analysis filter into multi - scale subsystem s in tim e scale dom ain. Some non - station