中梯法的电阻率正演模拟

2016年7月9日 477 次阅读 0 条评论 0 人点赞

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

一、实习目的:

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

二、要求:

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


四、实验结论及认识

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

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

菜鸟

文章评论(0)