《电法勘探及地电学》.docx
- 文档编号:15601182
- 上传时间:2023-07-06
- 格式:DOCX
- 页数:9
- 大小:607.51KB
《电法勘探及地电学》.docx
《《电法勘探及地电学》.docx》由会员分享,可在线阅读,更多相关《《电法勘探及地电学》.docx(9页珍藏版)》请在冰点文库上搜索。
《电法勘探及地电学》
《电法勘探及地电学》
基础理论实验报告
实习项目一中梯法的电阻率正演模拟
一、实习目的:
掌握中梯法球体正演场特征
二、要求:
设计球体的空间位置、物性参数和计算空间
根据中梯法球体的正演公式及设计空间编写程序;
’「2二r3代厂2x2
_2\[°(h2y2x2)5/2
ho-2x2
程序如下:
#include
#include
voidmain()
{
FILE*fp;
fp=fopen(”中间梯度电位.txt","w");
floatL,h,l,p,a,u1,u2,r0;
float
x[101][101]={0},y[101][101]={0},z[101][101]={0},z1[101][101]={0};inti,j;
printf("请输入中间梯度法参数,r0,u2:
");
scanf("%f%f",&r0,&u2);
h=2*r0;
for(i=0;i<101;i++)
for(j=0;j<101;j++)
{
x[i][j]=-15.0+j*0.3;
y[i][j]=-15.0+i*0.3;
z1[i][j]=1+2*(u2-1)/(2*u2+1)*r0*r0*r0*(h*h+y[i][j]*y[i][j]-2*x[i][j]*x[i][j])/pow((h*h+y[i][j]*y[i][j]+x[i][j]*x[i][j]),2.5);
fprintf(fp,"%f%f%f\n",x[i][j],y[i][j],z1[i][j]);
}fclose(fp);
}
球的空间位置(0,0,10),球体半径r0=3,
情况一:
u2=ps/p1=10
得下面电位图:
15
111
11051J1.0951.091.0&51.0B
1.075
1.07
1.055
1.06
1.055
1.05
1.045
1.04
1X36
1.03
1.025
102
1.015
1.01
1.005
1
Q.9950990.9&5093
0.975
实习项目二二层电测深电阻率正演模拟
一、实习目的:
掌握电测深电阻率正演方法及二层电测深曲线特征
二、要求:
根据二层电测深正演公式及设计参数编写程序;
K;2r
(r2(2nhj2)2
30
i(r2宀
nW
设计参数表:
第一层电阻率
50
60
60
100
250
40
30
20
50
5
第二层电阻率
40
30
20
50
5
50
60
60
100
250
第一层厚度
2
4
\6
\8
10
2
4
6
8
10
r:
{1.5,2,3,4,5,6,9,12,15,20,30,45,60,90}
编写程序如下:
#include
#include
main()
{
FILE*fp;
fp=fopen("测深数据.txt","w");
floatr[14]={1.5,2,3,4,5,6,9,12,15,20,30,45,60,90};
floatp1[10]={50,60,60,100,250,40,30,20,50,5};
floatp2[10]={40,30,20,50,5,50,60,60,100,250};
floath[10]={2,4,6,8,10,2,4,6,8,10};
floatps[14]={0};
inti,j,a,n=30;
floats,sun=0,k;
for(i=0;i<10;i++)
{
k=(p2[i]-p1[i])/(p2[i]+p1[i]);
for(a=0;a<14;a++)
{
for(j=1;j<=30;j++)
{
s=pow(k,j)*r[a]/pow((r[a]*r[a]+2*j*h[i]*2*j*h[i]),1.5);
sun=sun+s;
}
ps[a]=p1[i]*(1+2*r[a]*r[a]*sun);
sun=0;
printf("%f%f\n",log10(r[a]/h[i]),log10(ps[a]/p1[i]));
fprintf(fp,"%f%f\n",log10(r[a]/h[i]),log10(ps[a]/p1[i]));
}
printf("\n”);
fprintf(fp,"\n");
}
fclose(fp);
}
最终在EXCEL成图:
横坐标为log10(r/h),纵坐标为log10(ps/p1)(r=AB/2,h为深度,ps为视电阻率,p1为第
层点)
-2
实习项目3电偶极子电位正演
一、实习目的:
掌握电偶极子正演场特征
二、要求:
设计电偶极子的空间位置、物性参数和计算空间
根据电偶极子正演公式及设计空间编写程序;编写程序如下:
#include
#include
voidmain()
{
FILE*fp;
fp=fopen(”偶极子电位.txt","w");
floatL,h,I,p,a;
floatx[11][11]={0},y[11][11]={0},z[11][11]={0};
inti,j;
printf("此程序用于求偶极子电流在地面的电位和电位梯度,设偶极子中心
为原点,长度为L,中心到地面距离为h,电流大小为I,电阻率为p,偶极子倾斜角度a\n");
printf("请分别给出参数,长度,距离,电流大小,电阻率,角度:
\n");scanf("%f%f%f%f%f",&L,&h,&l,&p,&a);
for(i=0;i<11;i++)
for(j=0;j<11;j++)
z[i][j]=(l*p/(2*3.1415))*(1/pow(((x[i][j]+L*cos(a))*(x[i][j]+L*cos(a))+y[i][j]*y[i][j]+(h+L*sin(a))*(h+L*sin(a))),0.5)-1/pow(((x[i][j]-L*cos(a))*(x[i][j]-L*cos(a))+y[i][j]*y[i][j]+(h-L*sin(a))*(h-L*sin(a))),0.5));
fprintf(fp,"%f%f%f\n",x[i][j],y[i][j],z[i][j]);
}fclose(fp);
}
最终缩成电位图如图:
z[i][j]=(I*p/(2*3.1415))*(1/pow(((x[i][j]+L*cos(a))*(x[i][j]+L*cos(a))+y[i][j]*y[i][j]+(h+L*sin(a))*(h+L*sin(a))),0.5)-1/pow(((x[i][j]-L*cos(a))*(x[i][j]-L*cos(a))+y[i][j]*y[i][j]+(h-L*sin(a))*(h-L*sin(a))),0.5));
fprintf(fp,"%f%f%f\n",x[i][j],y[i][j],z[i][j]);
}
fclose(fp);
}
z[i][j]=(I*p/(2*3.1415))*(1/pow(((x[i][j]+L*cos(a))*(x[i][j]+L*cos(a))+y[i][j]*y[i][j]+(h+L*sin(a))*(h+L*sin(a))),0.5)-1/pow(((x[i][j]-L*cos(a))*(x[i][j]-L*cos(a))+y[i][j]*y[i][j]+(h-L*sin(a))*(h-L*sin(a))),0.5));
fprintf(fp,"%f%f%f\n",x[i][j],y[i][j],z[i][j]);
}
fclose(fp);
}
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 电法勘探及地电学 电法勘探 地电