成都理工大學(xué)電法勘探及地電學(xué)實(shí)習(xí)報告

《電法勘探及地電學(xué)》基礎(chǔ)理論實(shí)驗(yàn)報告學(xué)號：姓名：實(shí)習(xí)項(xiàng)目一中梯法的電阻率正演模擬一、實(shí)習(xí)目的：掌握中梯法球體正演場特征二、要求：設(shè)計球體的空間位置、物性參數(shù)和計算空間根據(jù)中梯法球體的正演公式及設(shè)計空間編寫程序；程序如下：#include<stdio.h>#include<math.h>voidmain(){ FILE*fp; fp=fopen("中間梯度電位.txt","w"); floatL,h,I,p,a,u1,u2,r0; floatx[101][101]={0},y[101][101]={0},z[101][101]={0},z1[101][101]={0}; inti,j; printf("請輸入中間梯度法參數(shù)，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得下面電位圖：實(shí)習(xí)項(xiàng)目二二層電測深電阻率正演模擬一、實(shí)習(xí)目的：掌握電測深電阻率正演方法及二層電測深曲線特征二、要求：根據(jù)二層電測深正演公式及設(shè)計參數(shù)編寫程序；設(shè)計參數(shù)表：第一層電阻率506060100250403020505第二層電阻率403020505506060100250第一層厚度246810246810r：{1.5,2,3,4,5,6,9,12,15,20,30,45,60,90}編寫程序如下：#include<stdio.h>#include<math.h>main(){ FILE*fp; fp=fopen("測深數(shù)據(jù).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成圖：橫坐標(biāo)為log10(r/h),縱坐標(biāo)為log10(ps/p1)(r=AB/2，h為深度，ps為視電阻率，p1為第一層點(diǎn))實(shí)習(xí)項(xiàng)目3電偶極子電位正演一、實(shí)習(xí)目的：掌握電偶極子正演場特征二、要求：設(shè)計電偶極子的空間位置、物性參數(shù)和計算空間根據(jù)電偶極子正演公式及設(shè)計空間編寫程序；編寫程序如下：#include<stdio.h>#include<math.h>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("此程序用于求偶極子電流在地面的電位和電位梯度，設(shè)偶極子中心為原點(diǎn)，長度為L，中心到地面距離為h，電流大小為I，電阻率為p,偶極子傾斜角度a\n "); printf("請分別給出參數(shù)，長度，距離，電流大小，電阻率，角度：\n"); scanf("%f%f%f%f%f",&L,&h,&I,&p,&a); for(i=0;i<11;i++) for(j=0;j<11;j++) { x[i][j]=-10.0+j*2; y[i][j]=-10.0+i*2; 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]*