中梯法的电阻率正演模拟

实习项目二 中梯法的电阻率正演模拟

一、实习目的：

掌握中梯法球体正演视电阻率异常特征

二、要求：

1、设计球体的空间位置、物性参数和计算空间；

2、根据中梯法球体的正演公式及设计空间编写程序；

3、熟悉不同电性差异、不同埋深、以及不同球体半径对异常特征影响。

三、实验过程：

1、试验计算原理：

1）平面视电阻率分布计算公式

2）剖面视电阻率分布计算公式

1. 空间球体分布结构图
   

图2 地下球体结构图

1. 实验参数
   

1）计算区域坐标范围___(-800，800)_______________；

2）球体几何参数：半径_____50___________；埋深_______500_________。

3) 物理参数：球体电阻率_____100___________________________（可取两个或以上）；围岩介质电阻率_______1000_________________________________（可取两个或以上，如取两者相对互为高阻或低阻）。

1. 程序设计及代码（可附打印源代码），通过其它计算方法本项可略。
   

#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）平面视电阻率异常图

四、实验结论及认识

当，即球体为低阻时，在球心正下方有极小值，两侧有的极大值；当，即球体为高阻时，在球心正上方有极大值，两侧则有的极小值，可确定球心在地面的投影的位置。

低阻球和高阻球的异常特征，均可用水平电流偶极子的电场分布规律加以解释。因此此时地下总的电流场乃是正常的均匀电流场与球体异常电流场（即水平电流偶极子的电流场）的叠加。