中梯装置脉状体ρ_s异常的计算

中梯装置脉状体ρs异常的计算姓名：scnjlybd学号：xxxxxxx专业：地球物理学【摘要】：有一位于均匀场中，电阻率为𝜌2的脉状体（椭球体），其长半轴为c、短半轴为a、埋深为h，求在围岩介质为均匀各向同性，电阻率为𝜌1的条件下，主剖面上椭球体的视电阻率分布。【关键词】：中梯装置；脉状体；视电阻率异常CalculationofApparentResistivityAnomalyofIntermediateGradientDevice【ABSTRACT】：Anellipsoidsituatedintheuniformfieldisprovidedofwhichresistivityis𝜌2,semi-majoraxisisc,semi-minoraxisisa,andburialdepthish.Nowweknowtheresistivityofthemediuminthehomogeneousandisotropicsurroundingrockis𝜌1,andthenseekfortheapparentresistivitydistributiononmainsectionoftheellipsoid.【KEYWORDS】：intermediategradientdevice;ellipsoid;apparentresistivityanomaly1方法原理脉状体是野外实际条件下，常见的一种地质体。对于这类形状的矿体，我们可以通过对三轴椭球体的理论计算方法或通过模型实验的方法研究ρs异常的分布特征及规律。对寻找有一定走向长度的脉状矿体而言，当AB垂直矿体走向供电时，称为纵向中梯或常用中梯；当AB沿着（或平行）矿体走向供电时便称为横向中梯。这里主要简要说明前者。1.1均匀电流场中椭球体的一次场解答及ρ表达式设在电阻率为ρ1的均匀各向同性无限岩石中有一电阻率为ρ2，3个半轴分别为a、b、c的椭球体。我们取直角坐标系，原点位于椭球体中心，坐标轴X,Y，Z分别与a、b、c重合。均匀外电场的方向与X轴一致，见图1。在椭球坐标系中解拉普拉斯方程并引用极限条件，可求解椭球内部和外部一次场的电位表达式Fig1均匀电流场中的椭球体（何裕盛等,1978）为：121𝜌1(1)111𝜌1(2)式中：π∫2√222π∫2√222L中的积分下限为以下方程的最大实根：222222以上为全空间情况，当有平行于XOY平面的地面（z=）存在时，地面一次场电位表达式，可用（2）式中异常部分加倍的办法近似地写出：1ρ1（3）沿X方向的一次场电场强度为11ρ1(4)式中：,为椭球体的几何系数。此时，L中的积分下限是方程222222之最大实根，其中为椭球中心深度。而222√222222[]由（4）式可写出地面上沿X方向的视电阻率表达式为𝜌1ρ1（5）对于情况，（5）式中的1，各量可以简化，以简化后的各量代入（5）式，便得到走向无限、柱轴平行行地面的椭球柱体𝜌表达式𝜌𝜌1{21[2⁄√√]}(6)1.2脉状体上的中梯𝜌异常可以理解，脉状体矿体上中梯𝜌异常的量值大小和分布特征，有仅与其埋藏深度和相对围岩电阻率的大小有关，而且还与脉状矿体的产状有关。当令（6）式中的ac时，可得到相当于直立薄脉的𝜌表达式。如令ac,则得相当于水平薄脉的𝜌表达式。对于倾斜薄脉的情况也可导出与其相应的𝜌表达式。根据这些𝜌表达式即可算出各种𝜌剖面曲线。考虑到理论计算和模型实验结果是一致的，为简单起见，这里主要用模型实验结果，对不同产状、不同电阻率的脉状矿体上之中梯𝜌异常进行讨论。1.2.1低阻脉状体上的中梯𝜌异常对于直立铜板来说，在适当条件下，由于均匀电流场的方向与模型走向方向垂直，故在矿顶上视电阻率的变化不明显，近于和围岩（水）电阻率值（约为20Ω·m）相等。因此，用中间梯度法寻找直立良导薄矿脉是不利的。当矿体倾斜时，𝜌则出现了明显异常。其特点是在倾斜方向上𝜌明显降底，在反倾斜方向上则𝜌有所升高，曲线呈不对称状。根据𝜌有极小值的那一边可指明矿体的倾向，并能说明矿体为低阻体。这是因为矿体倾斜时，均匀外电流场的方向与矿体斜交，因此吸引电流的能力加强，致使倾斜方向的近地表处电流密度𝑀𝑁减小（小于）,根据𝜌𝑗𝑀𝑁𝑗ρMN的关系式，应有𝜌𝜌1的低阻异常。1.2.2高阻脉状体上的中梯𝜌异常在高阻脉状体上，不论产状如何均有大于背景值（𝜌水）的𝜌高阻异常，但以直立高阻脉上的异常为最大。这是因为对高阻脉来说，当脉的产状为直立时，均匀外电流场的方向与脉的走向方向垂直，因而排斥电流的能力最强。当高阻脉为倾斜产状时，由于排斥电流的能力减弱，所以𝜌异常变小。曲线不对称且在倾斜方向上𝜌下降较快，有明显低于𝜌1的极小值。依此可判断高阻脉的倾斜方向，并能说明异常源为高阻体。当高阻脉的产状为水平时，由于排斥电流的能力最差，所以𝜌异常最小。若脉的厚度很薄，甚至观测不到𝜌异常。因此，用中间梯度法寻找高阻直立脉最有利，而寻找高阻水平脉最不利。对高阻情况，视电阻率极大值点约位于矿顶上，两侧有次级极小点（倾向一侧较明显）；而低阻体情况则不同，矿顶附近对应的视电阻率异常为零，两侧有量值相近的极小（倾向一侧）和极大（把倾向一侧）约呈倒像对应。视电阻率异常形态随𝜌2的变化是由于矿内、外异常电流的分布为装置类型和矿体的几何特征所决定的。虽然表达公式较为繁琐，但通过计算可以明确其间之相应关系，这无论对理论研究和野外资料解释均有实际意义。2实现过程2.1建立数学模型电阻率为𝜌2的椭球体矿体，其长半轴为c、短半轴为a、埋深为h，其围岩介质为均匀各向同性，电阻率为𝜌1，如图2所示。观测点位于主剖面上的视电阻率表达式：𝜌ρ1{2[√12⁄⁄]}(3-1)其中[2222222222]12⁄{2222[222222222]12⁄}2.2编程计算2.2.1计算要求a.求主剖面上ρ值，（x由-30.0m到+30.0m，步长为1.0m）；b.h=10,30,50m,ρ1Ω,ρ2Ω；c.c=3.0,5.0,8.0(m)a=0.005,0.01,0.05(m),要求按顺序对应搭配；d.将u与的表达式编成二个函数子程序言e.以ρ（或ρρ1⁄）为纵轴，x为横轴，打印出ρ曲线。2.2.1C程序源代码#includestdio.h#includemath.h#definer11.0//ρ1#definer2100.0//ρ2#defineh10//h=10,30,50m#definec3.0//c=3.0,5.0,8.0(m)#definea0.005//a=0.005,0.01,0.05(m)voidmain(){intx;doublers[61][2]={0};//ρsdoubleu[61]={0},px[61]={0};FILE*fp;fp=fopen(rs.txt,w);double*fun_u(doubleu1[]);double*fun_px(doublepx1[]);fun_u(u);fun_px(px);for(x=-30;x31;x++){rs[x+30][1]=r1*(1+(2*(r1-r2)*a*c)/((a-c)*(r2*a+r1*c))*(1-sqrt((c*c+u[x+30])/(a*a+u[x+30]))-(1.0/2)*(a*a-c*c)/(sqrt(c*c+u[x+30])*pow(a*a+u[x+30],1.5))*px[x+30]*x));rs[x+30][1]/=r1;rs[x+30][0]=x;fprintf(fp,%f%f\n,rs[x+30][0],rs[x+30][1]);printf(%f\n,rs[x+30][1]);}fclose(fp);}double*fun_u(doubleu1[]){intx;for(x=-30;x31;x++)u1[x+30]=(1.0/2)*(h*h+x*x-a*a-c*c)+(1.0/2)*sqrt(pow(x*x+h*h,2)+pow(a*a-c*c,2)-2*(a*a-c*c)*(x*x-h*h));returnu1;}double*fun_px(doublepx1[]){intx;for(x=-30;x31;x++)px1[x+30]=x*(1+(x*x+h*h-a*a+c*c)/sqrt(pow(x*x+h*h,2)+(a*a-c*c)-2*(a*a-c*c)*(x*x-h*h)));returnpx1;}2.2.1计算结果a.当H=10m,c=3.0m,a=0.005m时，所得结果如图3所示；b.当H=30m,c=5.0m,a=0.01m时，所得结果如图4所示；c.当H=50m,c=8.0m,a=0.05m时，所得结果如图5所示。3结论及几点认识a.当ρ2ρ1，即脉状体为低阻时，ρs或ρρ1在其主剖面的中心上方有极小值。因此，根据ρ曲线主极值点的坐标，可以确定球心在地面的投影位置。b.从探测原理上，中间梯度法非常适合直立的高阻异常体。c.利用电法勘探寻找矿产资源，是一种简便、快速、经济的方法，经过资料解释和推断，能够大致确定工区的资源分布情况，应综合其他物探资料和当地的地质资料，作出更加准确的解释。参考文献：[1]程荃等.地球物理与勘查技术与工程·课程设计指导讲义[M].成都理工大学信息工程学院，2005，25～26.[2]肖宏跃，雷宛.地电学教程[M].北京：地质出版社，2008，44～46.