二层电测深电阻率正演模拟

实习项目三 二层电测深电阻率正演模拟

一、实习目的：

巩固电测深理论，掌握电测深电阻率正演方法及二层电测深曲线特征

二、要求：

1、熟练掌握电测深正演参数意义；

2、熟练掌握电测深曲线特征；

3、会使用相关软件绘制电测深视电阻率双对数图。

三、实验过程

1、二层电测深正演计算原理

2、二层电测深计算实验参数表

表1 二层电测深参数表

r：{1.5,2,3,4,5,6,9,12,15,20,30,45,60,90}

3、程序设计及代码（可附打印源代码）

#include <stdio.h>

#include <stdlib.h>

#include<math.h>

/* run this program using the console pauser or add your own getch, system("pause") or input loop */

double jisuan(double rr,double pp1,double pp2,double hh);

int main(int argc, char *argv[])

{

    double ps[10][14];int i,j;

    double r[]={1.5,2,3,4,5,6,9,12,15,20,30,45,60,90};// AB/2

    double p1[]={50,60,60,100,250,40,30,20,50,5};//一层电阻率

    double p2[]={40,30,20,50,5,50,60,60,100,250};//二层电阻率

    double h[]={2,4,6,8,10,2,4,6,8,10};//层厚度

    FILE *fp;

    fp=fopen("xxx.txt","w");

    for (i=0;i<10;i++)

    {

     for(j=0;j<14;j++)

     {

     ps[i][j]=jisuan(r[j],p1[i],p2[i],h[i]);

     //fprintf(fp,"%f\t%f\t",r[j]/2/h[i],ps[i][j]/p1[i]);

     //printf("%f\t%f\t",r[j]/2/h[i],ps[i][j]/p1[i]);

         }

         //fprintf(fp,"\n");

        //printf("\n");

    }

    for (j=0;j<14;j++)

{

     for(i=0;i<10;i++)

     {

         printf("%f\t%f\t",r[j]/2/h[i],ps[i][j]/p1[i]);

         fprintf(fp,"%f\t%f\t",r[j]/2/h[i],ps[i][j]/p1[i]);

     }

    fprintf(fp,"\n");

        printf("\n");

}

return 0;

}

double jisuan(double rr,double pp1,double pp2,double hh)

{

    double K12,pps,temp;int n;

    K12=(pp2-pp1)/(pp2+pp1);

    //printf("%f",K12);

    pps=0.0;

    temp=0.0;

    for (n=1;n<=30;n++)

    {

        temp+=pow(K12,n)*rr/pow((rr*rr+pow(2*n*hh,2)) ,1.5);

    }

    pps=pp1*(1+2*rr*rr*temp);

    //printf("%f",pps);

    return pps;

}

1. 双对数坐标视电阻率异常图（给出每条曲线参数，可打印后粘贴空白处）
   

四、实验结论及认识

两层水平地层的上层电阻率为，厚度为，下层电阻率为，厚度为无限大。两层电阻率之比，当大于1时，下层电阻率高，称为G型曲线；当小于1时，下层电阻率低，称为D型曲线。

五、试验数据