- 1、本文档共17页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 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)