地电学实验.docx

1、地电学实验 地电学实验报告 姓名 , 学号 / 地电学实验报告实验一：偶极子不同角度电场和电场强度正演偶极子纵向放置，纵向为x，横向为y。对电位而言，偶极子正极和负极距地面近的一面电势绝对值大，正极则对应高电势，负极则低电势。对场强而言可通过电场线来判断，向右为正，左则负。大小根据电场线密度确定。或通过水平和垂直偶极子的情况旋转而得。 1 ：a=0； 电场强度 电位 2 a=30 E U 3 a=45 E U 4 a=60 E U 5 a=90 E U 6 a=120 E U 7=135 E U 8a=180 E U 实验二：球体中梯电阻率正演 1 ， u=p2/p1=0.1 2 u=10 X

2、为纵向，y为横向，可见高阻球体视电阻率较背景值高，低阻球体则低，同时低阻球体的异常比高阻明显实验三：电测深二层模板正演理解：对不同的u当u1*u2=1时，曲线关于x轴对称，曲线右只近与x轴平行，p2很大是，一斜率为1的直线上升，反之p2-0，以斜率-1向下。野外实验 1电测深：AB电极固定在中点两侧对称放MN，然后改变AB极距，测量完后，依次向右平移两米测下一个点，测点为MN中点，不同AB代表深度的大小。2中间梯度法 固定AB电极，MN在其横向中间2/3;纵向1/3;测点为MN中点副实验程序及野外试验数据实验一：/电流I ，电阻率 p，深度 h，角度 a，电极半长度l；/#include #i

3、nclude void main () float x50,y50,u50,ra,rb,I,l,h,pi,p,a; pi=3.14;I=10;l=5;h=10;p=0.01;a=pi;FILE *fp2 ,*fp; if(fp=fopen(正1.txt,w+)=NULL)return ; fp2=fopen(正1.txt ,w); int i,j; for(i=0;i50;i+) xi=-50+i*2; yi=-50+i*2; for(j=0;j50;j+) ra=0; rb=0; for(i=0;i50;i+) ra=sqrt(xi+l*cos(a)*(xi+l*cos(a)+(yj)*(y

4、j)+(h-l*sin(a)*(h-l*sin(a);rb=sqrt(xi-l*cos(a)*(xi-l*cos(a)+(yj)*(yj)+(h+l*sin(a)*(h+l*sin(a); ui=I*p/(2*pi)*(xi-l*cos(a)/(rb*rb*rb)-(xi+l*cos(a)/(ra*ra*ra); fprintf(fp2,n%f %f %fn,yj,xi,ui); fclose(fp2);实验二：#includestdio.h#includemath.hvoid main()float x100,y100,u100,h,r,v,m,n,k;FILE *fp;int i,j;h=

5、5;r=1；v=0.1;fp=fopen(r0.1.txt,w+);for(i=0;i100;i+) xi=(-5+i*0.1)*h; yi=(-5+i*0.1)*h; ui=0;for(i=0;i100;i+)m=0;n=0；k=0; for(j=0;j100;j+) m=h*h+yi*yi-2*xj*xj; k=h*h+yi*yi+xj*xj;n=sqrt(k*k*k*k*k);uj=1+2*(v-1)/(2*v+1)*(r*r*r)*(m/n); fprintf(fp,n %f %f %fn,yi,xj,uj);fclose(fp);实验三：#includestdio.h#include

6、math.hdouble g(double a,int b) int i;double s=1;for(i=0;ib;i+) s=s*a; return s;/d=2h/rdouble sum(double u,double d) int i;double s,m; m=(u-1)/(u+1); s=0;for(i=1;i200;i+) s=s+(g(m,i)/(sqrt(1+i*i*d*d)*sqrt(1+i*i*d*d)*sqrt(1+i*i*d*d) ) ) ;return (log10(2*s+1);void main() int j; FILE *fp,*fp1; fp=fopen(

7、w.txt,w); fp1=fopen(v.txt,w); float r14=1.5,2,3,4,5,6,9,12,15,20,30,45,60,90; for(j=0;j14;j+) printf(%f %fn,sum(0.1,4/rj),rj/4); fprintf(fp,%f %fn,sum(1.25,4/rj),log10(rj/4); fprintf(fp1,%f %fn,sum(0.8,4/rj),log10(rj/4); fclose(fp);fclose(fp1);野外试验数据1：电测深2.51019.382.51221.072.51419.362.51619.212.518

8、15.692.511019.412.511221.542.511419.772.511616.592.511820.4351016.2351217.1251416.0851615.2151814.3511015.83511215.42511416.58511615.35511815.49101014.16101214.54101413.56101612.92101810.351011015.661011213.181011413.761011612.591011812.33151013.97151214.9151413.92151613.51151814.141511013.561511213

9、.311511412.011511612.171511812.66202031.08202229.82202415.06202614.31202829.182021014.182021213.282021412.982021613.442021813.282：中间梯度法0121146.50102.06088.89066.010445.360221.4002262.80-2229.40-420.590-657.130-812.940-105718.20-121421.3312607.11310499.7138121.4436191.0834223.69321943.3301325.93-2108.063-41811.33-6141.633-8366.563-10368.893-1284