电测深正演模拟

1、本科生实验报告实验课程 电法与电磁法勘探 学院名称 地球物理学院 专业名称 勘查技术与工程 学生姓名 学生学号 指导教师 实验地点 地球物理学院实验室5417 实验成绩 二 年 月 二 年 月填写说明1、 适用于本科生所有的实验报告（印制实验报告册除外）；2、 专业填写为专业全称，有专业方向的用小括号标明；3、 格式要求： 用A4纸双面打印（封面双面打印）或在A4大小纸上用蓝黑色水笔书写。 打印排版：正文用宋体小四号，1.5倍行距，页边距采取默认形式（上下2.54cm，左右2.54cm，页眉1.5cm，页脚1.75cm）。字符间距为默认值（缩放100%，间距：标准）；页码用小五号字底端居中。

2、具体要求：题目（二号黑体居中）；摘要（“摘要”二字用小二号黑体居中，隔行书写摘要的文字部分，小4号宋体）；关键词（隔行顶格书写“关键词”三字，提炼3-5个关键词，用分号隔开，小4号黑体)； 正文部分采用三级标题；第1章 ××(小二号黑体居中，段前0.5行)1.1 ×××××小三号黑体×××××（段前、段后0.5行）1.1.1小四号黑体（段前、段后0.5行）参考文献（黑体小二号居中，段前0.5行），参考文献用五号宋体，参照参考文献著录规则（GB/T 77142005）。13成都

3、理工大学本科生课程实验报告水平层状介质视电阻率测深正演模拟摘 要电测深曲线在实际工作中有很大的用处，掌握理解曲线变化的本质对时间勘探工作有指导意义；本次实验用C语言编程计算了在地下为三层介质的情况下，地面同一点电阻率随深度增加而变化的视电阻率，并用绘图软件绘制了电测深图件分别的到A型、Q型、H型和K型电阻率测深曲线，最后分析了第一层电阻率变化及第一层深度变化所对应的曲线有哪些不同，得出在第一层深度先相同的情况下，电阻率增大导致曲线起止点的纵坐标增大，在电阻率相同的情况下，深度增大会使曲线首支近似水平的那一段长度增大；并且在实验过程中发现，第二层深度越大曲线中间的拐点越突出；这些曲线变化能帮助我

4、们判断地下介质的层数、电阻率变化、第一层的电阻率及深度等； 关键词：层状介质；视电阻率测深；正演；第3章 水平层状介质视电阻率测深正演模拟3.1 实验目的根据水平层状介质条件下电场理论导出的电测深视电阻率计算公式，设计程序计算水平层状介质视电阻率测深曲线，分析不同类型地电断面对应的视电阻率测深曲线特征，以及厚度变化对视电阻的影响规律。3.2 实验内容（1）根据水平层状介质电测深视电阻率计算公式，设计计算方法，设置相应的计算参数，完成二层或多层介质的视电阻率测深正演计算。（2）根据不同地电断面的计算结果绘制视电阻率测深曲线图，分析并总结不同电性及厚度变化时地面视电阻曲线的变化特征。3.3 实验设

5、备本次实验为理论计算，用到的主要设备为个人计算机。需要的工具软件分别是程序设计平台，及成果图件绘制软件Golden Software Grapher 7.0，完成测深曲线的绘制。3.4 实验步骤3.4.1 实验原理水平层状地层理论模型示意图如图3-1所示。根据水平层状介质电场理论，计算在地面某点视电阻率测深曲线。图3-1 水平层状地层理论模型示意图计算方法：根据多层层状介质的理论，地面观测点的电位可记为 （3-9）利用电场强度和视电阻之间的关系式： （3-10）上述视电阻率是供电极距的函数（=AB/2），因此上式可改写为： （3-11）令，则（3-11）变为 （3-12）即为电阻率转换函数。对

6、于N层介质，且具有以下递推关系：当n层以上全去除后，在第n层以下为均匀空间，则有第n层顶界面上 （3-13）当存在第n-1层和第n层这两层时，第n-1层顶界面上有 （3-14）上式（3-14）可改写为： （3-15）式（3-15）称为层状介质电阻率转换函数的递推公式。直流电测深计算视电阻率的公式（3-12）的积分在计算机求解时采用数字滤波法计算，通常采用20点滤波器实现，具体形式（3-16）为 （3-16）其中，为供电极距AB/2，一般取；，为位移系数，计算中取，为滤波系数，取值如下表3-1。表3-1 视电阻率测深正演采用的20点滤波系数表kCkkCkkCkkCk10.00304260.237

7、40112.7044016-0.025212-0.00119870.6194012-1.13240170.0112530.0128481.18170130.3930018-0.00497840.0235090.4248014-0.14360190.00207250.0868810-3.45070150.0581220-0.000318采用本方法计算层状介质视电阻率测深曲线时，只需要根据（3-16）式编写程序，通过输入不同的供电极距，即可获得不同供电极距的视电阻率。但应注意供电极距不是任意设定的，需要考虑采样值与供电极距的关系。计算流程为：（1）输入层参数，包括层数N，各层的层厚度和电阻率，存入

8、相应的数组中；（2）输入要计算的个供电极距值，存入对应的数组中；（3）读取要计算的第个供电极距值，（4）根据滤波系数序列计算第个值，即； （5）用电阻率转换函数递推公式，循环计算对应的（6）将计算得到的第个与第个滤波系数相乘，重复步骤（4）-（5）；（6）将得到的20个值求和即可得到供电极距时的视电阻率值； （7）重复步骤（3）-（6），即可获得所有供电极距对应的视电阻率值。（8）输出供电极距与的值，即可获得电测生正演结果。3.4.2 电测深视电阻率计算过程（1）计算参数设计 供电极距的计算与设置计算视电阻率采用如表3-2的供电极距参数，表中的参数系根据计算得到，计算结果如表3-2，共35个供

9、电极距。表3-2 视电阻率测深采用的供电极距参数表序号AB/2（m）序号AB/2（m）序号AB/2（m）序号AB/2（m）11.211528106.8129211938.31186728215.44345521.467799118.2540422046.41588629261.01570431.7782791210.0000002156.23413030316.22774442.1544351312.1152762268.12920331383.11865852.6101571414.6779922382.54041432464.15884963.1622781517.7827932499.9

10、9999433562.34128273.8311871621.5443462512129201684.6415891726.10157126146.77991835825.40411995.6234131831.62277527177.82793036999.999917 模型参数设置：本次试验中需要计算的模型参数如下表3-3，共12个模型，分为四组。表3-3 水平层状三层介质模型层参数表模型组号模型序号第一层电阻率第二层电阻率第三层电率第一层厚度第二层厚度第一组120.00100.00200.004.050.0220.00100.00200.008.050.0330

11、.00100.00200.0012.050.0440.00100.00200.0012.050.0第二组5200.00100.00204.050.06200.00100.00208.050.07220.00100.002012.050.08240.00100.002012.050.0第三组9100.0020.00200.004.016.010100.0020.00200.008.016.011120.0020.00200.0012.016.012140.0020.00200.0012.016.0第四组13100.00200.00304.016.014100.00200.00308.016.01

12、5120.00200.003012.016.016140.00200.003012.016.0（2）计算程序设计与代码 程序设计语言C语言程序 程序中的变量说明程序中使用的主要变量名、数据类型及其数值含义说明如表3-4。表3-4 程序中使用的主要变量及其说明序号变量名数据类型功能含义备 注1N整型地层层数输入参数2i、j、k、n整型循环参数3H浮点型数组名存放输入的地层深度数据输入参数3Rho浮点型数组名存放输入的地层电阻率数据输入参数5T浮点型数组名存放计算中得出的电阻率转换函数数据6file_name字符型数组名存放输入索要输出的文件名输入参数 程序源代码及相关注释#include<

13、stdlib.h>#include<stdio.h>#include<math.h>#include<string.h>void main()void dian_ce_shen(int N,float*h,float*rho,float*T);/声明计算电测深函数int N;float*h,*rho,*T;printf("输入地层层数（至少2层）："); scanf("%d",&N);h=(float*)malloc(N-1)*sizeof(float); /h,rho,T均为定义动态数组，分别表示每层的

14、深度、电阻率、电阻率转换函数rho=(float*)malloc(N*sizeof(float);T=(float*)malloc(N*sizeof(float);dian_ce_shen(N,h,rho,T);free(h);free(rho);free(T);void dian_ce_shen(int N,float*h,float*rho,float*T)int i,j,k,n=1; /循环变量float r=0.0,rho_s=0.0,mj=0.0; /r为电极距、rho_s为视电阻率，均初始为零， float delta=log(10)/6,s=-2.1719; /delta为采样间

15、隔，s表示位移系数float C20=0.0; /定义固定长度的数组，存放滤波系数char file_name256=0; /该数组用于存放指定的文件名printf("输入要建立的文件名:");scanf("%s",file_name); /输入文件名strcat(file_name,".txt");FILE*fp_LBXS,*fp_rhos; fp_rhos=fopen(file_name,"w");fp_LBXS=fopen("滤波系数.txt","r");for(i=0

16、;i<20;i+) /将滤波系数读入到一个数组fscanf(fp_LBXS,"%f",&Ci);fclose(fp_LBXS);for(i=0;i<N-1;i+) /输入各层的深度和电阻率printf("第%d层的深度和电阻率：",i+1);scanf("%f %f",&hi,&rhoi);printf("第%d层电阻率：",N);scanf("%f",&rhoN-1);TN-1=rhoN-1;while(1) rho_s=0.0;r=exp(0.5*

17、n*delta);printf("%fn",r); /电极距的大小，可以改变0.5来改变电极距的密度n+;if(r>1000)break; /我们只需要电极距在1000米以内for(j=1;j<21;j+)mj=exp(j*delta+s)/r;for(k=0;k<N-1;k+) /用循环体计算T1TN-k-2=rhoN-k-2*(rhoN-k-2*(1-exp(-2*mj*hN-k-2)+TN-k-1*(1+exp(-2*mj*hN-k-2)/(rhoN-k-2*(1+exp(-2*mj*hN-k-2)+TN-k-1*(1-exp(-2*mj*hN-k-

18、2);/电阻率转换函数的地推公式rho_s=T0*Cj-1+rho_s;fprintf(fp_rhos,"%f %fn",r,rho_s);fclose(fp_rhos); 输出数据文件名称及格式输出文件名称：根据运行程序时输入的文件名而把计算的数据输出到响应的文件中；成果数据文件格式：二列，第1列为AB/2值，第2列为供电极距对应的视电阻率值。（3）二层地电模型电测深视电阻率计算结果计算获得的不同模型电测深视电阻率结果（列出前8个模型计算结果）如表3-5。表3-5 不同模型电测深视电阻正演计算结果供电极距模型序号对应的视电阻率AB/2模型1模型2模型3模型4模型5模型6模

19、型7模型81.21152820.09820.01230.00440.004199.57199.94219.98239.981.46779920.17320.02230.00740.008199.25199.9219.96239.961.77827920.30220.03930.01440.014198.69199.82219.93239.922.15443520.52320.0730.02540.026197.74199.68219.88239.862.61015720.89320.12430.04540.047196.16199.44219.78239.753.16227821.49820.

20、21830.0840.085193.61199.03219.62239.553.83118722.45320.3830.14240.15189.65198.31219.33239.214.64158923.920.65530.2540.265183.8197.1218.83238.625.62341325.97721.11230.43640.461175.71195.1217.96237.66.81292128.7921.85130.75240.796165.36191.9216.49235.878.25404232.36923.00331.2841.355153.23186.99214.06

21、2331036.67324.71932.13642.261140.25179.83210.13228.3912.1152841.61927.14433.47943.681127.46170.05204.05221.2614.6779947.13530.37435.49345.811115.55157.62195.08210.7817.7827953.18634.4438.36148.842104.58142.98182.64196.3121.5443559.78439.32242.21652.91593.983126.89166.55177.7326.1015766.97744.98747.1

22、1458.08583.055110.14147.25155.6731.6227874.83451.43253.03964.33171.40893.294125.79131.4838.3118783.40958.68159.94371.59659.3576.831103.55106.8746.4158992.70766.76367.78379.81947.84361.3882.07283.59956.23413102.6575.67576.51388.92238.05147.86962.84563.23768.1292113.0685.35486.05698.77830.75337.21247.

23、16446.99982.54041123.795.67896.276109.1926.01429.80635.75835.43999.99999134.26106.46106.98119.8923.29925.30528.48328.214121.1528144.47117.47117.93130.621.87222.87324.42124.256146.7799154.04128.46128.85141.0221.13321.64222.37822.292177.8279162.75139.13139.47150.8820.72921.00921.38321.341215.4435170.4

24、5149.23149.52159.9320.48920.65720.86820.847261.0157177.05158.52158.7716820.33820.44420.57520.565316.2277182.53166.82167.03174.9920.2420.3120.39520.39383.1187186.95174.02174.18180.8420.17520.22120.27920.277464.1588190.42180.07180.2185.6120.13220.16320.20220.203562.3413193.08185.01185.11189.3920.10220

25、.12420.15220.154681.292195.06188.92189192.2920.08220.09820.11820.121825.4041196.51191.95192194.4820.06720.0820.09520.099999.9999197.55194.22194.27196.0920.05720.06820.0820.084（4）四组模型视电阻率测深曲线图及其特征分析通过计算，获得的视电阻率测深数据绘制的测深曲线图如下图3-2、图3-3、图3-4、图3-5，共四幅图，每个图为一组模型的计算结果，采用的坐标为双对数坐标图3-2 第一组模型电测深曲线图 图3-3 第二组模型

26、电测深曲线图图3-4 第三组模型电测深曲线图 图3-5 第四组模型电测深曲线图视电阻率测深曲线图特征分析：对于第一组曲线，曲线类型为A类型，曲线整体增大，最终趋于第三层电阻率值，曲线首支起始点的纵坐标随着第一层的电阻率增大而增大；在第一层电阻率相同的情况下，第一层深度小的曲线相对深度更大的曲线在增大到第二层电阻率时斜率更显增大；曲线中段出现两个拐点，在第一层电阻率与第二层电阻率之间和第二层与第三层之间各一个， 曲线尾支趋于第三层电阻率；对于第二组曲线，曲线类型为Q类型，曲线整体逐渐减小，最终趋于第三层电阻率值，曲线起始点纵坐标随着第一层的电阻率增大而增大，在第一层电阻率相同的情况下，第一层深度小的曲线相对深度更大的曲线在增大到第二层电阻率时斜率更显增大；曲线中段出现两个拐点，在第