电测深正演模拟(共15页).docx
《电测深正演模拟(共15页).docx》由会员分享,可在线阅读,更多相关《电测深正演模拟(共15页).docx(15页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
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)。专心-专注-专业水平层状介质视电阻率测深正演模拟摘 要电测深曲线在实际工作中有很大的用处,掌握理解曲线变化的本质对时间勘探工作
3、有指导意义;本次实验用C语言编程计算了在地下为三层介质的情况下,地面同一点电阻率随深度增加而变化的视电阻率,并用绘图软件绘制了电测深图件分别的到A型、Q型、H型和K型电阻率测深曲线,最后分析了第一层电阻率变化及第一层深度变化所对应的曲线有哪些不同,得出在第一层深度先相同的情况下,电阻率增大导致曲线起止点的纵坐标增大,在电阻率相同的情况下,深度增大会使曲线首支近似水平的那一段长度增大;并且在实验过程中发现,第二层深度越大曲线中间的拐点越突出;这些曲线变化能帮助我们判断地下介质的层数、电阻率变化、第一层的电阻率及深度等; 关键词:层状介质;视电阻率测深;正演;第3章 水平层状介质视电阻率测深正演模
4、拟3.1 实验目的根据水平层状介质条件下电场理论导出的电测深视电阻率计算公式,设计程序计算水平层状介质视电阻率测深曲线,分析不同类型地电断面对应的视电阻率测深曲线特征,以及厚度变化对视电阻的影响规律。3.2 实验内容(1)根据水平层状介质电测深视电阻率计算公式,设计计算方法,设置相应的计算参数,完成二层或多层介质的视电阻率测深正演计算。(2)根据不同地电断面的计算结果绘制视电阻率测深曲线图,分析并总结不同电性及厚度变化时地面视电阻曲线的变化特征。3.3 实验设备本次实验为理论计算,用到的主要设备为个人计算机。需要的工具软件分别是程序设计平台,及成果图件绘制软件Golden Software G
5、rapher 7.0,完成测深曲线的绘制。3.4 实验步骤3.4.1 实验原理水平层状地层理论模型示意图如图3-1所示。根据水平层状介质电场理论,计算在地面某点视电阻率测深曲线。图3-1 水平层状地层理论模型示意图计算方法:根据多层层状介质的理论,地面观测点的电位可记为 (3-9)利用电场强度和视电阻之间的关系式: (3-10)上述视电阻率是供电极距的函数(=AB/2),因此上式可改写为: (3-11)令,则(3-11)变为 (3-12)即为电阻率转换函数。对于N层介质,且具有以下递推关系:当n层以上全去除后,在第n层以下为均匀空间,则有第n层顶界面上 (3-13)当存在第n-1层和第n层这两
6、层时,第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.60.23740112.7044016-0.025212-0.70.6194012-1.13240170.0112530.0128481.18170130.39
7、30018-0.40.0235090.4248014-0.14360190.50.0868810-3.45070150.0581220-0.采用本方法计算层状介质视电阻率测深曲线时,只需要根据(3-16)式编写程序,通过输入不同的供电极距,即可获得不同供电极距的视电阻率。但应注意供电极距不是任意设定的,需要考虑采样值与供电极距的关系。计算流程为:(1)输入层参数,包括层数N,各层的层厚度和电阻率,存入相应的数组中;(2)输入要计算的个供电极距值,存入对应的数组中;(3)读取要计算的第个供电极距值,(4)根据滤波系数序列计算第个值,即; (5)用电阻率转换函数递推公式,循环计算对应的(6)将计算
8、得到的第个与第个滤波系数相乘,重复步骤(4)-(5);(6)将得到的20个值求和即可得到供电极距时的视电阻率值; (7)重复步骤(3)-(6),即可获得所有供电极距对应的视电阻率值。(8)输出供电极距与的值,即可获得电测生正演结果。3.4.2 电测深视电阻率计算过程(1)计算参数设计 供电极距的计算与设置计算视电阻率采用如表3-2的供电极距参数,表中的参数系根据计算得到,计算结果如表3-2,共35个供电极距。表3-2 视电阻率测深采用的供电极距参数表序号AB/2(m)序号AB/2(m)序号AB/2(m)序号AB/2(m)11.106.1938.28215.21.118.2046.29261.3
9、1.1210.2156.30316.42.1312.2268.31383.52.1414.2382.32464.63.1517.2499.33562.73.1621.25121.34681.84.1726.26146.35825.95.1831.27177.36999. 模型参数设置:本次试验中需要计算的模型参数如下表3-3,共12个模型,分为四组。表3-3 水平层状三层介质模型层参数表模型组号模型序号第一层电阻率第二层电阻率第三层电率第一层厚度第二层厚度第一组120.00100.00200.004.050.0220.00100.00200.008.050.0330.00100.00200.0
10、012.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.015120.00200.003
11、012.016.016140.00200.003012.016.0(2)计算程序设计与代码 程序设计语言C语言程序 程序中的变量说明程序中使用的主要变量名、数据类型及其数值含义说明如表3-4。表3-4 程序中使用的主要变量及其说明序号变量名数据类型功能含义备 注1N整型地层层数输入参数2i、j、k、n整型循环参数3H浮点型数组名存放输入的地层深度数据输入参数3Rho浮点型数组名存放输入的地层电阻率数据输入参数5T浮点型数组名存放计算中得出的电阻率转换函数数据6file_name字符型数组名存放输入索要输出的文件名输入参数 程序源代码及相关注释#include#include#include#i
12、ncludevoid 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均为定义动态数组,分别表示每层的深度、电阻率、电阻率转换函数rho=(float*)malloc(N*sizeof(float);T=(float*)malloc(N*sizeof(float);dian_ce_shen(N,h,rho,
13、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为采样间隔,s表示位移系数float C20=0.0; /定义固定长度的数组,存放滤波系数char file_name256=0; /该数组用于存放指定的文件名printf(输入要建立的文件名:);scanf(%s
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 测深 模拟 15
限制150内