电法实验报告.docx

电法实验报告.docx电法实验报告
实习项目一
中梯法的电阻率正演模拟
一、实习冃的: 掌握中梯法球体正演场特征
二、实验要求： 设计球体的空间位置、物性参数和计算空间
根据中梯法球体的止演公式及设计空间编写程序;
作出工演成果的平面等值线图和中心剖面图 比较和分析不同电性比值时屮心剖面图的变化规律
三、球体模型
在地下半空间条件下，可按镜像法原理将半空间映射成全空间，用地
lz
而上方的一个镜像球代替地面影响。若以球心在 地面投影点0为原点，Z轴垂直向下，地面观测 点坐标为M (x, y, 0),球心坐标为(0,0, h0), 如右图所示。利用全空间球外电位表达式中异常 部分加倍的方法，可求得地面电位，进一步求得 任一点电场强度，再根据视电阻率微分表达式得 中梯法球体的止演公式:
Ps
3 + h2 - 2x2
% n 2 , 2~.~~2x5/2
(九 +y +X )
Ps 二］| ? “2 _ 1 3 嶠 _ 2/
P~ 2“2+1°(硝 + 疋严
四、程序代码：
〃======平面等值线=======11
#include<>
#include<>
void main()
{
FILE *fp;
fp=fopen(n ”);
float pl=l ,p2= ,ps,r0=2,h0= 10,x,y;
for(y=-20;y<20;y+4-)
for(x=-20;x<20;x++)
{
ps=pl*(l+2*(p2-pl)/(p2+pl)*r0*r0* r0*((h0*h0+y*y-2*x*x)/pow(h0 *h0+y*y+x*x,)));〃 计算电阻率//
fprintf(fp,n%f\t%f\t%f\nn,x,y,ps);
}
fclose(fp);
//===========中心剖面==========//
#include<>
#include<>
void main()
{
FILE *fp;
fp=fopen(n ；,wH);
float pl=1 ,p2= ,r0=2,h0= 10,x,s;//用 s 表示 ps/pl//
for(x=-20;x<20;x+4-)
{
S=l+2*(p2・pl)/(p2+pl)*r0*r0*r0*((h0*h0-2*x*x)/pow(h0*h0+x*x,2 ・5));〃计算电阻率〃
fprintf(fp,n%f\t%f\nn,x,s);
fclose(fp);
四、结果图件
15
-5
低阻平面等值线图
C C1014
C C1013
C C1012

C C101
C C10D9
C C10G8
C C10IJ7
C C1005
C C10C5
CC1004
C C1003
CC1002
C C1001
C C1
C C0999
CC0998
C C0997
高阻平面等值线图
电阻率剖血图
实习项目二 二层电测深电阻率正演模拟
一、 实习目的：
掌握电测深电阻率正演方法及二层电测深曲线特征
二、 要求：
根据二层电测深正演公式及设计参数编写程序；
30 r
or（1 + 2 述 —）
心（厂2 + （2必］）2 /
设计参数表:
第一层电阻率
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<math・ h> void main()
{
FILE *fp;
fp=fopen(，", "w");
double s,u2,k , ul;
double r[] ={,2, 3, 4, 5, 6, 9, 12, 15, 20, 30, 45, 60, 90};
double pl [] = {50,