实习项目二 中梯法的电阻率正演模拟
一、实习目的:
掌握中梯法球体正演视电阻率异常特征
二、要求:
1、设计球体的空间位置、物性参数和计算空间;
2、根据中梯法球体的正演公式及设计空间编写程序;
3、熟悉不同电性差异、不同埋深、以及不同球体半径对异常特征影响。
三、实验过程:
1、试验计算原理:
1)平面视电阻率分布计算公式
2)剖面视电阻率分布计算公式
-
空间球体分布结构图
图2 地下球体结构图
-
实验参数
1)计算区域坐标范围___(-800,800)_______________;
2)球体几何参数:半径_____50___________;埋深_______500_________。
3) 物理参数:球体电阻率_____100___________________________(可取两个或以上);围岩介质电阻率_______1000_________________________________(可取两个或以上,如取两者相对互为高阻或低阻)。
-
程序设计及代码(可附打印源代码),通过其它计算方法本项可略。
#include<stdio.h>
#include<math.h>
#include<stdlib.h>
int main()
{ int x,y;float pps,ps,p1=100,p2=1000,r0=50,h0=500;
FILE *fp,*fpp;
fp=fopen("xxx.txt","r+");
fpp=fopen("www.txt","r+");
for(x=-800;x<=800;x+=10)
{
for(y=-800;y<=800;y++)
{
ps=p1*(1+2*(p2-p1)/(2*p2+p1)*r0*r0*r0*(h0*h0+y*y-2*x*x)/pow((h0*h0+y*y+x*x),2.5));
printf("%d\t%d\t%f\n",x,y,ps);
fprintf(fp,"%d\t%d\t%.10f\n",x,y,ps);
}
}
for(x=-800;x<=800;x+=10)
{
pps=p1*(1+2*(p2/p1-1)/(p2/p1+1)*r0*r0*r0*(h0*h0-2*x*x)/pow((h0*h0+x*x),2.5));
printf("%d\t%.10f\n",x,pps);
fprintf(fpp,"%d\t%.10f\n",x,pps);
}
return 0;
}5、视电阻率异常图(可打印后粘贴空白处)
1)主剖面视电阻率异常图
2)平面视电阻率异常图
四、实验结论及认识
当,即球体为低阻时,在球心正下方有极小值,两侧有的极大值;当,即球体为高阻时,在球心正上方有极大值,两侧则有的极小值,可确定球心在地面的投影的位置。
低阻球和高阻球的异常特征,均可用水平电流偶极子的电场分布规律加以解释。因此此时地下总的电流场乃是正常的均匀电流场与球体异常电流场(即水平电流偶极子的电流场)的叠加。
文章评论