刘宇-作业5-自动化反演-2016.3.31探究.docx

  1. 1、本文档共17页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
中国地质大学(武汉)地空学院 姓 名: 刘 宇 班 级: 061132班 学 号: 20131003811 方法原理 一、正演计算 水平层状介质对称四极法电侧深的积分关系式(当MN-0时)可表示为 在上式中, 对理论转换电阻率 T1(λ) 进行积分, 便可以算出理论电阻率ρs 。 但是直接对该式积分较为困难, 所以采用数字滤波方法运算。经过代换、 褶积等一系列运算后, 可得到离散化的数字逮波形式如下 (2) C为滤波系数。I为取点样号。K为序号,s是位移 理论转换电阻率的递推公式 由下而上、如下 (3) 滤波系数C通过查找资料给出,观测值也能给出,就能根据公式算出视电阻率的值,通过MATLAB程序实现的结果如下: clear close all clc format short g nlayer=3 %层数 res(1)=10 res(2)=100 res(3)=5 thick(1)=2 thick(2)=8 ab2=[2 3 4.5 6 9 12 15 20 30 45 60 90 120 150 200 300 450 600 900 1200 1500 2000 3000 4500] nab2=length(ab2); for i=1:nab2 r=ab2(i); [rhos]=sdcs1dford(nlayer,res,thick,r); rhoscal(i)=rhos; end ab2=ab2(:); rhoscal=rhoscal(:); clear a a= [ab2 rhoscal]; dlmwrite(d01obs.dat,a); %画图 loglog(ab2,rhoscal,r-o,... linewidth,2,... MarkerEdgeColor,k,... MarkerFaceColor,g,... MarkerSize,10) axis equal grid on xlabel(AB/2 / m) ylabel(RHO / ohm) 副程序 function [rhos]=sdcs1dford(nlayer,res,thick,r) c=[0.003042 -0.001198 0.01284 0.0235 0.08688 0.2374 0.6194 1.1817 0.4248 -3.4507 2.7044 -1.1324 0.393 -0.1436 0.05812 -0.02521 0.01125 -0.004978 0.002072 -0.000318]; % rhos=10; % whos; for k=1:20 mk=exp(k*log(10)/6-2.1719)/r; %m的值 T(nlayer)=res(nlayer); for i=nlayer-1:-1:1 t1=T(i+1)+res(i); t2=T(i+1)-res(i); t3=exp(-2*mk*thick(i)); t4=t1+t2*t3; t5=t1-t2*t3; T(i)=res(i)*t4/t5; %求T的值 end T1(k)=T(1); %限定值的范围 end rhos=0; for k=1:20 rhos=rhos+T1(k)*c(k); end 反演拟合的实现 对于任一条电侧深曲线来说, 总会有某一段曲线对某一层层参数的改变反应特别明显, 而对其他层层参数的改变反应较为迟缓, 也就是说, 一个电性层对应着一段测深曲线段。 由此就可以建立电性参数和侧深曲线上的各曲线段一一对应的关系。 在反演拟合中,如果理论曲线和实测曲线的某一段重合不好、

文档评论(0)

希望之星 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档