基于ANSYS的复杂地形条件下电法数据地形改正研究.pdf

第3 7卷 第3期 V o l.3 7,N o.32 0 2 3年6月 M I N E R A L R E S OUR C E S AN D G E O L OG YJ u n.,2 0 2 3D O I:1 0.1 9 8 5 6/j.c n k i.i s s n.1 0 0 1-5 6 6 3.2 0 2 3.0 3.0 2 1基于AN S Y S的复杂地形条件下电法数据地形改正研究尚建阁1,2,3,李晨晖1,4,苏丽娟1,4,李 冰1,2,3,谢 崇1,4,邢尚鑫1,5,6(1.河南省有色金属矿产探测工程技术研究中心,河南 郑州 4 5 0 0 1 6;2.河南省有色金属地质矿产局第二地质大队,河南 郑州 4 5 0 0 1 6;3.河南省自然资源科技创新中心(矿山生态环境保护修复研究),河南 郑州 4 5 0 0 1 6;4.河南有色地质矿产集团有限公司,河南 郑州 4 5 0 0 1 6;5.河南省有色金属地质勘查总院,河南 郑州 4 5 0 0 5 2;6.河南省有色金属深部找矿勘查技术研究重点实验室,河南 郑州 4 5 0 0 5 2)摘 要:直流激发极化法是现代地质勘查中应用较为广泛的地球物理探测方法之一,该方法能快速获取地质体的极化率和电阻率等电性参数。但受地形起伏的影响,通过该方法所获取的电阻率参数往往不是地质体的真实反映。运用AN S Y S有限元分析软件能快速准确的对复杂地形条件下不同装置类型的电阻率开展二维、三维地形改正,使改正后的电阻率参数能反映地质体的真实情况。通过在研究区运用AN S Y S软件对野外采集的电阻率数据进行地形改正,取得了良好效果,可以为广大地球物理工作者提供技术参考。关键词:AN S Y S软件;复杂地形;电法数据;有限元分析;地形改正中图分类号:P 6 3 1.3;T P 3 1 9 文献标识码:A 文章编号:1 0 0 1-5 6 6 3(2 0 2 3)0 3-0 5 9 7-0 80 引言随着全球能源矿产的需求量不断增加,找矿工作的程度也越来越深入。直流激发极化法是一种常规的物探方法,也是目前应用最为广泛的物探方法之一,通过该方法可以快速的获取地球内部介质的极化率、电阻率等物理参数,但地球表面并不是理论上的一个光滑的平面,而是一个地形高低不平、连绵起伏的复杂曲面,特别是在山区,地形变化更加复杂,正是由于这些复杂地形的存在,在使用直流激发极化法开展找矿过程中,发现的异常会产生畸变,可能导致淹没地质体的有用信息,呈现假异常,使研究者不能做出正确 的推断解释,降低了物探 综合参数 利 用 水平1-3。因此,复杂地形条件下电阻率地形改正方法的研究,成为物探工作者急需解决的关键问题之一。AN S Y S软件是一款大型的融电磁学、结构力学、流体力学、仿声学等于 一体的大型 通用有限元 分 析 软件4。AN S Y S具有丰富的单元库,提供了对各种物理场量的分析功能5-6,为地球物理工作者在进行数值模拟分析时提供了强大的数学分析功能。前人运用AN S Y S开展了一系列的地球物理电法数据的数值模拟,对典型的地电模型全空间介质影响进行模拟,验证了AN S Y S开展电法数据模拟的正确性,均取得了较为理想的理论效果7-9,但是对复杂地形条件下的模拟应用较少。文章通过介绍运用AN S Y S开展电法地形改正数值模拟,总结经验并运用到实际生产找矿工作中,取得了良好应用效果,为广大地球物理工作者在复杂地形条件下开展直流激发极化法数据解释推断提供技术参考。收稿日期:2 0 2 2-0 5-1 7;修回日期:2 0 2 2-0 6-2 4基金项目:河南省有色金属地质矿产局科技创新项目(编号:Y S D K 2 0 1 1-1 2)与河南省2 0 2 2年度财政地质勘查项目(豫自然资发2 0 2 2-2 2)共同资助。第一作者:尚建阁(1 9 8 4),女,工程师,主要从事地质及地球物理科研工作。E-m a i l:2 7 0 9 1 8 8 7 2q q.c o m通信作者:李冰(1 9 8 4),男,硕士,高级工程师,主要从事地球物理科研工作。E-m a i l:5 4 0 0 3 5 8 0 1q q.c o m引文格式:尚建阁,李晨晖,苏丽娟,等.基于AN S Y S的复杂地形条件下电法数据地形改正研究J.矿产与地质,2 0 2 3,3 7(3):5 9 7-6 0 3,6 1 0.1 方法原理地形影响问题是从2 0世纪5 0年代伴随着直流电法在矿产资源勘查等方面的大规模应用开始受到重视的1 0。复杂的地形对直流激发极化法中电阻率参数的影响是致命的,因此,正确消除起伏地形对电阻率参数的影响,对后期地球物理成果解释推断的真实性起到了至关重要的作用1 1。研究地形对视电阻率的影响,目前主要使用的有三种方法:物理模拟法、解析法和数值模拟法。所谓直流电法的地形改正,就是通过上述三种方法将在起伏地形上测量得到的电阻率数据转换为平地形影响下数据的过程,从而消除不同地形所引起的误差。由于物理模拟法和解析法不适用于场源和物性分布复杂的正演情况,因此,数值模拟法则成为电阻率法正演最主要的方法1 2。随着计算机技术的发展和物探解释方法的进步,数值模拟技术得到了飞跃式发展。常用的电阻率数值模拟方法有有限差分法、积分方程法、边界单元法和有限元法,其中有限元法具有较高的求解精度且理论较为成熟,能够适应任意模型形态。因此,有限元法是目前处理电法正演问题的主流方法。AN S Y S是一款能实现多场及多场耦合分析、实现前后处理、求解及多场分析统一数据库的一体化大型有限元分析软件,一般而言,AN S Y S的基本分析过程可以分为三步,即创建有限元模型、施加载荷进行求解和查看分析结果。对应软件结构的三个程序模块:前 处 理 模 块(P R E P 7)、分 析 求 解 模 块(S O L U-T I ON)和后处理模块(P O S T 1和P O S T 2 6)。前处理模块(P R E P 7)的主要任务是建立结构分析的有限单元模型,它是结构分析的开始,一般步骤为分析准备、设置单元类型、设定实常数、定义材料属性、创建模型、划分网格。分析求解模块(S O L UT I ON)的主要任务是进行结构的计算,一般步骤:定义自由度、施加荷载、设定分析类型、给定边界条件及求解等。后处理模块(P O S T 1和P O S T 2 6)的主要任务是进行结构的分析,它包括通用后处理器(P O S T 1)和时程后处理器(P S OT 2 6)两种,前者主要用于静态分析的结果后处理,而后者主要用于动态分析结果的后处理。2 模拟实验利用AN S Y S软件进行地形改正,其一般思路:根据电法工作的实际情况,依照测点高程建立物理模型,并要根据计算精度,对模型进行网格剖分和边界条件的设定;根据工作精度的要求,合理选择单元类型,利用AN S Y S程序进行纯地形条件下各测点的视电阻率值计算;通过比值法,对实测视电阻率数据进行地形影响的压制。计算公式采用s改=(s实/s地形)0,通过AN S Y S有限元分析得到的纯地形视电阻率数据,对实测数据进行比值压制,最终得到改正后的s改数据。总的来说,利用AN S Y S软件进行地形改正需要正确记录实际地形高程数据,利用软件所提供的有限元分析软件计算出纯地形条件下影响的数值,再利用实际测量获得的电阻率数值与之进行比值法进行地形影响消除得到改正后的真实数值,基本处理流程见图1。图1 AN S Y S有限元数据分析流程图F i g.1 AN S Y S f i n i t e e l e m e n t d a t a a n a l y s i s f l o w c h a r t下面利用AN S Y S编写命令流加以说明。(1)!进入前处理器,建模/p r e p 7/P NUM,a r e a,1e t,1,p l a n e 2 3 0 !定义单元类型m p,1,r e s i s t i v i t y !定义材料电阻率参数*a s k,K P_N,8,!定义创建地形模型的点数*d i m,K P_N X,A R R AY,K P_N,895矿 产 与 地 质2 0 2 3年*d i m,K P_NY,A R R AY,K P_N,!定义地形节点的X Y数值k,I,K P_N X(I),K P_NY(I)!定义关健节点a,1,2,3,4,5,6,7,8 !由点创建面(2)进入求解器!将模型网格划分,定义边界条件和施加荷载a s e l,s,1,2a a t t,1 a s e l,s,4a a t t,2a s e l,s,1,4e s i z e,2 0 !定义网络划分尺度m s h a p e,1,2 d !定义剖分类型a m e s h,a l l !开始网络划分l s s c a l e,a l ln s e l,s,e x t !选择边界上的节点n s e l,u,l o c,y*g e t,N n o d,n o d e,c o u n t*g e t,Nm i n,n o d e,n u m,m i n*d o,J,1,N n o d,1D i s t 1=S Q R T(N X(Nm i n)-NA X)*2+(NY(Nm i n)-NAY(Nm i n)*2)D i s t 2=S Q R T(N X(Nm i n)-N B X)*2+(NY(Nm i n)-N B Y(Nm i n)*2)V 0=C u r r t*R e s 1*(1/D i s t 1-1/D i s t 2)/(2*P i)D,Nm i n,VO L T,V 0 !给边界上的节点赋值Nm i n=n d n e x t(Nm i n)*e n d d of k,2,a m p s,1f k,1 1,a m p s,-1 !施加电流荷载f i n i s h/s o la n t y p e,s t a t i c!/s t a t u s,s o l u !定义分析类型和求解s o l v eF I N I S H(3)!进入后处理器,查看结果/p o s t*S E T,D J,0.2*S E T,MN,0.2*D I M,V 1,A R R AY,2 0,1,*D I M,R a,T A B L E,2 0,2,*D O,I,1,2 0,1V 1(I,1)=VO L T(NO D E(2+D J*(I-1),6,0)-VO L T(NO D E(2+D J*I,6,0)!求模拟电位差AM=2+D J*(I-1)AN=2+D J*IBM=6-D J*(I-1)B N=6-D J*IKK=2*3.1 4 1 5 6/(1/AM-1/AN-1/BM+1/B N)R a(I,1)=2+I*D JR a(I,2)=V 1(I,1)*k k/1 0*e n d d oF I N I S H/G C O L,1,d i s/G C O L,2,R e s i s t i v i t y*V P L OT,R a(1,1),R a(1,2),2F I N I S H995 第3 7卷 第3期 尚建阁等:基于AN S Y S的复杂地形条件下电法数据地形改正研究 根据AN S Y S软件 进行地形改 正的思 路 和 分析流程,分别开展了不 同坡度角的 山谷 地 形 和 山脊地形对 直 流 激 发 极 化 法 电 阻 率 影 响 的 数 值 模拟(图2)。图2 a为 山 脊 模 型,为 坡 度 角,材 料电阻率为1 0 0 m,为 了 研 究 山 脊 纯 地 形 对 视电阻率的影响程度,分别进行了坡度为0、1 0、2 0、3 0、4 0 的模型计算。由图2 a可见,在山脊纯地形上得到的视电阻率曲线呈现为低阻特征,并且随着坡度角的增大,其低阻特征更加明显。图2 b为山谷模型,为坡度角,材料电阻率为1 0 0 m,同样也进行了坡度为0、1 0、2 0、3 0、4 0 的模型计算。由图2 b可见,在山谷纯地形上得到的视电阻率曲线呈现为高阻特征,并且随着坡度角的增大,其高阻特征更加明显。由图2可见,纯地形异常主要发生在角域顶点附近,也就是说主要发生在坡度变化的地方。异常的大小与坡度角的大小有关,即与地形起伏的陡缓程度关系密切。形态相同的角域,异常幅度随坡度的增大而增大。图2 不同地形条件影响下模拟电阻率曲线图F i g.2 S i m u l a t e d r e s i s t i v i t y c u r v e u n d e r t h e i n f l u e n c e o f d i f f e r e n t t e r r a i n c o n d i t i o n s(a)不同坡度角山脊纯地形影响模拟电阻率曲线(b)不同坡度角山谷纯地形影响模拟电阻率曲线(c)不同坡度角山脊地形下高阻直立薄板模拟电阻率曲线(d)不同坡度角山谷地形下低阻水平薄板模拟电阻率曲线(e)高阻直立薄板山脊地形改正后模拟电阻率曲线(f)低阻水平薄板山谷地形改正后模拟电阻率曲线 为了更好地研究地形对异常体电阻率曲线的影响,在数值模拟器中加入异常体进行数值模拟,图2 c为山脊地形下高阻直立薄板电阻率曲线图,为坡度角,背 景电阻率为1 0 0 m,直 立 高 阻 体电阻率为1 0 0 0 m。由图2 c可见,在水平地面时,直立薄板视电阻率异常表现为高阻特征,随着地形坡度的变大,异常曲线表现为低阻特征,尤其当坡度角为3 0、4 0 时,基本为纯地形引起的异常曲线特征。图2 d为不同山谷地形下低阻水平薄板电阻率曲线,为坡度角,背景电阻率为1 0 0 m,水平薄板的电阻率为1 0 m。由图2 d可见,在水平地面时,水平薄板电阻率曲线表现为低阻特征,随着地形坡度的变大,异常曲线表现为高阻特征,尤其当坡度角为3 0、4 0 时,基本为纯地形引起的异常曲线特征。由此可以看出,地形起伏对异常体引起的电阻率曲线影响很大,特别是当地形切割严006矿 产 与 地 质2 0 2 3年 重时,对电阻率资料的影响更大,使异常体电阻率异常发现较 大畸变,甚 至出现相 反 的 结 果。图2 e为高阻直立薄板山脊地形改正前后的电阻率曲线。由图2 e可见,蓝色曲线为地改前的电阻率曲线,当存在地形影响时,高阻体也表现为低阻特征;红色曲线为地改后的电阻率曲线,通过比值法改正后,电阻率曲线表现出高阻特征,与实际地下介质情况相符。图2 f为高阻直立薄板山脊地形改正前后的电阻率曲线。由图2 f可见,蓝色曲线为 地改前的电阻率 曲 线,当 存 在 地 形 影 响 时,由 于 地 形 的 影响,低阻体在电阻率曲线上也表现为高阻特征;红色曲线为 地 改 后 的 电 阻 率 曲 线,通 过 比 值 法 改 正后,电阻率曲线表现出低阻特征,与实际地断面介质情况相符。通过电流源有限元分析计算,其结果充分表明,地形对电阻率的影响很大,特别是地形起伏较大地区,影响程度更大。通过地形模拟实验可以总结得出地形对电阻率能产生较大影响,具有正地形展现出低阻特性,负地形展现出高阻特性,这样的影响随着地形的复杂程度不断提高显现的越发明显,甚至能抵消地下目标体所产生的真实电阻率影响,使得目标体的真实特征不能被探获。利用AN S Y S软件进 行快速地形 改 正 有 效 地 解 决 了地形起伏对电阻率的影响,使反映的介质电阻率回归为真实。在使用AN S Y S软件进行数值分析的过程中,网格剖分的大小直接影响着计算结果的精度,所以网格剖分越细,计算精度也就有所提高。但是,网格划分的大小也影响着计算工作量的大小,直接影响着计算速度。因此,在应用中应权衡这两种因素,做到综合考虑。图3为不同网格剖分密度下的山脊地形电阻率曲线。由图3可见,随着剖分尺度的大小变化,曲线的光滑程度也发生变化,且剖分尺度越大,曲线振荡程度越小,说明计算精度越高。在AN S Y S稳态电流有限元分析中,程序没有提供远边界单元类型,在实际应用时,须构建一个远边界模型,对研究区和边界区进行不同的网格划分。3 应用实例河南省西部山区是我国重要的有色金属成矿区带之一,该区域分布有多个大大小小的有色贵金属矿床,为河南省矿产资源需求提供了保障1 3-1 4。研究区位于河南省西部洛宁县,熊耳山北坡,其大地构造位置处于中朝准地台南缘、华熊台缘凹陷、崤山鲁山图3 二维纯山脊地形电阻率曲线图(不同剖分尺度)F i g.3 T w o-d i m e n s i o n a l t o p o g r a p h y r e s i s t i v i t y c u r v e o fp u r e r i d g e(d i f f e r e n t s u b d i v i s i o n s c a l e s)拱褶断束中部(图4)。地层出露主要为太古宇太华群的角闪斜长片麻岩、熊耳群的安山岩和第四系。区内岩浆活动频繁,岩浆岩种类较多,其中以燕山期酸性岩最为发育,并与本区内生金属矿产的形成关系极为密切。断裂构造十分发育,并以N ENN E向断裂为主,且规模较大。研究区围岩蚀变主要为硅化、绢云母化、绿泥石化、黄铁矿化、泥化等。硅化表现为石英沿破碎带充填,形成巨厚层的硅化碎裂岩带出现。构造蚀变带内矿化与蚀变主要有黄铁矿化、方铅矿化、碳酸盐化、高岭土化、褐铁矿化、硅化、绢云母化等,金属矿物表现形式为星点状、浸染状、细脉状、细脉侵染状、团块状。表1为研究区电性参数统计结果,由表1可见,研究区矿石(含矿化体)与围岩存在着较大的物性差异,目标体电性表现特征为低阻高极化特征,为在该研究区开展物探工作提供物性前提。本次物探研究工作是在试验区K 9号矿脉上进行,采用时间域激电偶极-偶极测深装置,点距2 0 m,隔离系数为5,MN极距4 0 m。仪器为V 8,最大供电电流为1 4 A,6个电通道方式。图5为试验剖面测量结果,由图5 a可见,极化率断面图,极化率高值异常与K 9号矿脉基本对应,反映出地形对极化率测量结果影响较小,而图5 b中所反映电阻率异常与矿脉不对应,并且展现出了高电阻率的特征,这一高电阻率特征是地下矿脉真实电阻率与地形影响产生的电阻率叠加影响产生的综合反映,其低阻异常向左侧地形相对低缓的一侧进行偏离。这是由于地形切割较大,106 第3 7卷 第3期 尚建阁等:基于AN S Y S的复杂地形条件下电法数据地形改正研究图4 区域地质矿产图(a)与研究区地质简图(b)F i g.4 R e g i o n a l g e o l o g y m a p o f m i n e r a l r e s o u r c e s(a)a n d g e o l o g i c a l s k e t c h m a p o f t h e s t u d y a r e a(b)1第四系 2新生界砂砾岩 3熊耳群安山岩 4官道口群石英砂岩及白云岩 5太华群变质岩 6中生代花岗岩类 7拆离断层 8断层 9地质界线 1 0平行不整合地质界线 1 1银(铅)矿床 1 2金矿床 1 3矿体及编号 1 4试验剖面表1 研究区电性参数统计结果T a b l e 1 S t a t i s t i c a l r e s u l t o f e l e c t r i c a l p a r a m e t e r s i n t h e s u r v e y a r e a岩(矿)石名称标本量/块极化率/%电阻率/(m)最小值最大值平均值最小值最大值平均值片麻岩4 20.2 71.5 81.0 5 51 8 4 73 2 7 6 28 5 9 1安山岩4 90.9 14.2 22.8 98 7 3 91 1 5 2 4 35 6 3 6 7蚀变岩3 50.1 41.0 50.4 9 93 2 7 52 0 1 3 91 1 3 1 2矿石(矿化)3 43.8 98.7 46.0 31 24 2 71 6 4山谷地形纯电阻率异常反映为高阻,致使矿体引起的低电阻异常被掩 盖,发生偏 离现象。图5 c为采用AN S Y S软件计 算出该地形 所 产 生 的 响 应 值 与实际探测的电阻率值的比值进行地形改正计算,其结果可见低电阻异常位置向右侧偏移,并与K 9号矿脉基本对应,彻底消除了复杂地形对探测结果的影响效应。结合图5 a和图5 c结果显示 与岩矿石物性结果一致,表现为低阻高极率特征,证明采用AN S Y S软件进行电阻率地形改正使测量结果回归真实反映,效果显著。4 结论(1)激发极化法是地质找矿中应用最为广泛的地球物理方法之一,但是其电阻率数据收复杂地形影响也最为严重,负地形将形成高电阻率反映,而正地形将产生低电阻率反映,从而实际物探测量所获得电阻率带有地形加载的影响,使得大地介质的真实电阻率得不到有效体现。(2)AN S Y S软件是一款能实现多场及多场耦合分析、实现前后处理、求解及多场分析统一数据库的一体化大型有限元分析软件,能快速准确进行多种场景数值模拟计算,能够对复杂地形条件下不同装置类型的电阻率开展二维、三维地形改正,利用该软件有限元分析功能模拟计算出纯地形下产生的电阻率响应,通过比值法将野外所获取的电阻率值中地形影响消除,获得大地介质真实电阻率,实现复杂地形对电阻率值的零影响效应,为物探解释的可靠性提供保障。(3)应用AN S Y S软件在研究区复杂地形条件下开展电阻率地形改正实际应用,对比结果显示,运用该软件进行地形改正方法可行,效果显著。206矿 产 与 地 质2 0 2 3年 图5 试验剖面电法资料二维反演断面图F i g.5 T w o-d i m e n s i o n a l i n v e r s i o n p r o f i l e o f t e s t s e c t i o n e l e c t r i c a l d a t a(a)试验剖面视极化率二维反演断面图(b)试验剖面视电阻率二维反演断面图(地形改正前)(c)试验剖面视电阻率二维反演断面图(地形改正后)参考文献:1 何继善,周正秀.电阻率法地形改正中的点电源场与线电源场问题J.地质与勘探,1 9 7 5(8):4 5-5 1.2 葛欢,卢君实,刘志强.高密度电法在场地地质综合评价中的应用J.矿产勘查,2 0 1 8,9(7):1 3 9 4-1 3 9 9.3 王烨,杨晓弘,戴亦军,等.基于有限单元法的激电参数正演数值模拟J.物探化探计算技术,2 0 0 9,3 1(1):4 4-4 7.4 孙明礼,胡仁喜,崔海蓉,等.AN S Y S 1 0.0电磁学有限元分析实例指导教程M.北京:机械工业出版社,2 0 0 7.5 赵斌.基于AN S Y S的激发极化法中地形影响问题的数值模拟D.长沙:中南大学,2 0 0 7.6 汤井田,肖晓,杜华坤,等.AN S Y S在直流电法正演中的应用J.地球物理学进展,2 0 0 6(3):9 8 7-9 9 2.7 阮百尧,熊彬.电导率连续变化的三维电阻率测深有限元模拟J.地球物理学报,2 0 0 2,4 5(1):1 3 1-1 3 8.8 高炳魁.基于AN S Y S直流电测深法的边坡勘探正演D.长沙:中南大学,2 0 0 8.9 张军,赵莹,马炳镇.基于AN S Y S的矿井直流电场研究J.地球物理学进展,2 0 1 2,2 7(6):2 6 0 9-2 6 1 6.1 0 徐自生,张丽,唐伟,等.复电阻率法(C R)在内蒙古那吉河铅锌矿探测中的应用J.矿产勘查,2 0 1 5,6(2):1 6 5-1 7 0.1 1 侯俊富,索朗顿珠,张革利,等.激发极化法在西藏昂青银多金属矿勘查中的应用J.矿产勘查,2 0 1 9,1 0(2):3 0 2-3 0 9.1 2 刘强,吴家祥,赵志成,等.在龙门店银铅矿区利用AN S Y S进行直流电法地形改正的研究J.矿产勘查,2 0 1 9,1 0(8):2 0 2 9-2 0 3 5.1 3 范少军,何文平.河南省洛宁县龙门店银铅矿地质特征及找矿前景J.矿产与地质,2 0 1 5,2 9(3):3 0 6-3 1 0.1 4 王杰,王小高,程宝成,等.豫西龙门店银铅锌矿碳氧同位素特征及其地质意义J.地质找矿论丛,2 0 2 1,3 6(2):1 6 5-1 7 0.(下转第6 1 0页)306 第3 7卷 第3期 尚建阁等:基于AN S Y S的复杂地形条件下电法数据地形改正研究r a d a r r e f l e c t e d w a v e a t t h e i n t e r f a c e b e t w e e n g r o u n d a n d u n d e r g r o u n d h i d d e n d a n g e r s t h r o u g h i n d o o r n u m e r i-c a l s i m u l a t i o n a n d a c t u a l c a s e c o m p a r i s o n.T h r o u g h t h e a n a l y s i s,i t c a n b e s e e n t h a t t h e a m p l i t u d e o f g r o u n d r e f l e c t e d w a v e i s s t r o n g,a n d t h e p h a s e i s o p p o s i t e t o t h a t o f r a d a r t r a n s m i t t i n g s o u r c e.T h e a m p l i t u d e o f t h e r e f l e c t e d w a v e o n t h e u p p e r i n t e r f a c e o f t h e c a v i t y a n d v o i d i s a l s o s t r o n g,a n d i t s p h a s e i s t h e s a m e a s t h a t o f t h e r a d a r t r a n s m i t t i n g s o u r c e a n d o p p o s i t e t o t h a t o f t h e g r o u n d r e f l e c t e d w a v e.T h e r e f o r e,i n p r a c t i c a l w o r k,t h e i d e n t i f i c a t i o n a c c u r a c y o f u n d e r g r o u n d h i d d e n d a n g e r s c a n b e i m p r o v e d b y c o m b i n i n g t h e a m p l i t u d e i n t e n s i t y a n d p h a s e r e l a t i o n s h i p o f g r o u n d r e f l e c t e d w a v e s a n d s u s p e c t e d a b n o r m a l r e f l e c t e d w a v e s w i t h o t h e r c h a r a c t e r i s t i c s o f u n d e r g r o u n d h i d d e n d a n g e r s,s u c h a s t i p d i f f r a c t i o n a n d b o t t o m r e f l e c t e d w a v e s.K e y w o r d s:g r o u n d p e n e t r a t i n g r a d a r,u n d e r g r o u n d c a v i t y,r e f l e c t e d w a v e,p h a s e,a m p l i t u d e,n u m e r i c a l s i m-u l a t i o n(上接第6 0 3页)R e s e a r c h o n t e r r a i n c o r r e c t i o n o f e l e c t r i c a l m e t h o d d a t au n d e r c o m p l e x t e r r a i n c o n d i t i o n b a s e d o n A N S Y SS HANG J i a n g e1,2,3,L I C h e n h u i1,4,S U L i j u a n1,4,L I B i n g1,2,3,X I E C h o n g1,4,X I NG S h a n g x i n1,5,6(1.N o n-F e r r o u s M i n e r a l E x p l o r a t i o n E n g i n e e r i n g T e c h n o l o g y R e s e a r c h C e n t e r o f H e n a n P r o v i n c e,Z h e n g z h o u 4 5 0 0 1 6,H e n a n,C h i n a;2.G e o l o g i c a l T e a m N o.2,B u r e a u o f N o n f e r r o u s M e t a l G e o l o g y a n d M i n e r a l R e s o u r c e s o f H e n a n P r o v i n c e,Z h e n g z h o u 4 5 0 0 1 6,H e n a n,C h i n a;3.N a t u r a l R e s o u r c e s S c i e n c e a n d T e c h n o l o g y I n n o v a t i o n C e n t e r o f H e n a n P r o v i n c e(R e s e a r c h o n M i n e E c o l o g i c a l E n v i r o n m e n t P r o t e c t i o n a n d R e s t o r a t i o n),Z h e n g z h o u 4 5 0 0 1 6,H e n a n,C h i n a;4.H e n a n N o n-f e r r o u s G e o l o g y a n d M i n e r a l R e s o u r c e s G r o u p C o.,L t d.,Z h e n g z h o u 4 5 0 0 1 6,H e n a n,C h i n a;5.G e n e r a l I n s t i t u t e o f N o n-f e r r o u s M e t a l s G e o l o g i c a l E x p l o r a t i o n o f H e n a n P r o v i n c e,Z h e n g z h o u 4 5 0 0 5 2,H e n a n,C h i n a;6.K e y L a b o r a t o r y o f D e e p O r e-p r o s p e c t i n g T e c h n o l o g y R e s e a r c h f o r N o n-f e r r o u s M e t a l s o fH e n a n P r o v i n c e,Z h e n g z h o u 4 5 0 0 5 2,H e n a n,C h i n a)A b s t r a c t:D C p o l a r i z a t i o n m e t h o d i s o n e o f t h e m o s t w i d e l y u s e d g e o p h y s i c a l e x p l o r a t i o n m e t h o d s i n m o d e r n g e o l o g i c a l e x p l o r a t i o n.I t c a n q u i c k l y o b t a i n t h e e l e c t r i c a l p a r a m e t e r s s u c h a s p o l a r i z a t i o n a n d r e s i s t i v i t y o f g e o l o g i c a l b o d i e s.H o w e v e r,d u e t o t h e i n f l u e n c e o f t o p o g r a p h i c r e l i e f,t h e r e s i s t i v i t y p a r a m e t e r s o b t a i n e d b y t h i s m e t h o d a r e o f t e n n o t t h e t r u e r e f l e c t i o n o f g e o l o g i c a l b o d y.AN S Y S f i n i t e e l e m e n t a n a l y s i s s o f t w a r e c a n q u i c k l y a n d a c c u r a t e l y c a r r y o u t t w o-d i m e n s i o n a l a n d t h r e e-d i m e n s i o n a l t e r r a i n c o r r e c t i o n f o r r e s i s t i v i t y o f d i f f e r e n t d e v i c e t y p e s u n d e r t h e c o m p l e x t e r r a i n c o n d i t i o n,s o t h a t,t h e r e s i s t i v i t y p a r a m e t e r s a f t e r c o r r e c t i o n c a n r e f l e c t t h e r e a l s i t u a t i o n o f g e o l o g i c a l b o d y.B y u s i n g AN S Y S s o f t w a r e i n t h e s t u d y a r e a,t e r r a i n c o r r e c-t i o n o f r e s i s t i v i t y d a t a c o l l e c t e d i n t h e f i e l d h a s a c h i e v e d t h e g o o d r e s u l t,w h i c h c a n p r o v i d e a t e c h n i c a l r e f e r-e n c e f o r t h e m a j o r i t y o f g e o p h y s i c i s t s.K e y w o r d s:AN S Y S s o f t w a r e,c o m p l e x t e r r a i n,e l e c t