地电学实验.docx
《地电学实验.docx》由会员分享，可在线阅读，更多相关《地电学实验.docx（19页珍藏版）》请在冰豆网上搜索。
地电学实验
地电学实验报告
姓名,,,,,
学号//////
地电学实验报告
实验一：
偶极子不同角度电场和电场强度正演
偶极子纵向放置，纵向为x，横向为y。
对电位而言，偶极子正极和负极距地面近的一面电势绝对值大，正极则对应高电势，负极则低电势。
对场强而言可通过电场线来判断，向右为正，左则负。
大小根据电场线密度确定。
或通过水平和垂直偶极子的情况旋转而得。
1：
a=0；
电场强度电位
2a=30
EU
3a=45
EU
4a=60
EU
5a=90
EU
6a=120
EU
7=135
EU
8a=180
EU
实验二：
球体中梯电阻率正演
1，u=p2/p1=0.1
2u=10
X为纵向，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
#include
voidmain（）
{
floatx[50],y[50],u[50],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"）;
inti,j;
for（i=0;i<50;i++）
{
x[i]=-50+i*2;
y[i]=-50+i*2;
}
for（j=0;j<50;j++）
{ra=0;rb=0;
for（i=0;i<50;i++）
{ra=sqrt（（x[i]+l*cos（a））*（x[i]+l*cos（a））+（y[j]）*（y[j]）+（h-l*sin（a））*（h-l*sin（a）））;
rb=sqrt（（x[i]-l*cos（a））*（x[i]-l*cos（a））+（y[j]）*（y[j]）+（h+l*sin（a））*（h+l*sin（a）））;
u[i]=I*p/（2*pi）*（（x[i]-l*cos（a））/（rb*rb*rb）-（x[i]+l*cos（a））/（ra*ra*ra））;
fprintf（fp2,"\n%f%f%f\n",y[j],x[i],u[i]）;}
}
fclose（fp2）;
}
实验二：
#include"stdio.h"
#include"math.h"
voidmain（）
{
floatx[100],y[100],u[100],h,r,v,m,n,k;
FILE*fp;
inti,j;
h=5;r=1；v=0.1;
fp=fopen（"r0.1.txt","w+"）;
for（i=0;i<100;i++）
{
x[i]=（-5+i*0.1）*h;
y[i]=（-5+i*0.1）*h;
u[i]=0;
}
for（i=0;i<100;i++）
{m=0;n=0；k=0;
for（j=0;j<100;j++）
{m=h*h+y[i]*y[i]-2*x[j]*x[j];
k=h*h+y[i]*y[i]+x[j]*x[j];
n=sqrt（k*k*k*k*k）;u[j]=1+2*（v-1）/（2*v+1）*（r*r*r）*（m/n）;
fprintf（fp,"\n%f%f%f\n",y[i],x[j],u[j]）;}}
fclose（fp）;}
实验三：
#include"stdio.h"
#include"math.h"
doubleg（doublea,intb）
{inti;
doubles=1;
for（i=0;i{s=s*a;
}
returns;
}
//d=2h/r
doublesum（doubleu,doubled）
{
inti;
doubles,m;
m=（u-1）/（u+1）;
s=0;
for（i=1;i<200;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））;
}
voidmain（）
{intj;
FILE*fp,*fp1;
fp=fopen（"w.txt","w"）;
fp1=fopen（"v.txt","w"）;
floatr[14]={1.5,2,3,4,5,6,9,12,15,20,30,45,60,90};
for（j=0;j<14;j++）
{
printf（"%f%f\n",sum（0.1,4/r[j]）,r[j]/4）;
fprintf（fp,"%f%f\n",sum（1.25,4/r[j]）,log10（r[j]/4））;
fprintf（fp1,"%f%f\n",sum（0.8,4/r[j]）,log10（r[j]/4））;
}
fclose（fp）;
fclose（fp1）;}
野外试验数据
1：
电测深
2.5
1
0
19.38
2.5
1
2
21.07
2.5
1
4
19.36
2.5
1
6
19.21
2.5
1
8
15.69
2.5
1
10
19.41
2.5
1
12
21.54
2.5
1
14
19.77
2.5
1
16
16.59
2.5
1
18
20.43
5
1
0
16.23
5
1
2
17.12
5
1
4
16.08
5
1
6
15.21
5
1
8
14.3
5
1
10
15.83
5
1
12
15.42
5
1
14
16.58
5
1
16
15.35
5
1
18
15.49
10
1
0
14.16
10
1
2
14.54
10
1
4
13.56
10
1
6
12.92
10
1
8
10.35
10
1
10
15.66
10
1
12
13.18
10
1
14
13.76
10
1
16
12.59
10
1
18
12.33
15
1
0
13.97
15
1
2
14.9
15
1
4
13.92
15
1
6
13.51
15
1
8
14.14
15
1
10
13.56
15
1
12
13.31
15
1
14
12.01
15
1
16
12.17
15
1
18
12.66
20
2
0
31.08
20
2
2
29.82
20
2
4
15.06
20
2
6
14.31
20
2
8
29.18
20
2
10
14.18
20
2
12
13.28
20
2
14
12.98
20
2
16
13.44
20
2
18
13.28
2：
中间梯度法
0
12
1146.5
0
10
2.06
0
8
8.89
0
6
6.01
0
4
45.36
0
2
21.4
0
0
2262.8
0
-2
229.4
0
-4
20.59
0
-6
57.13
0
-8
12.94
0
-10
5718.2
0
-12
1421.3
3
12
607.11
3
10
499.71
3
8
121.44
3
6
191.08
3
4
223.69
3
2
1943.3
3
0
1325.9
3
-2
108.06
3
-4
1811.3
3
-6
141.63
3
-8
366.56
3
-10
368.89
3
-12
84
升级会员 