一维电测深反演实验

《一维电测深反演实验》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、煤田电法勘探数字解释基础实验报告姓名：龙斌学号：05102344班级：地物10-2班指导老师：姜志海、刘志新2013年11月中国矿业大学实验一一维电测深反演实验一、实验目的通过对课堂学习的一维电测深反演的基础理论知识，对一维电测深反演的方法（拟合核函数法、直接拟合视电阻率法、DZ法）进一步的加强学习并练习，通过对程序的编制更深层次的理解一维电测深反演的基本思想与基本步骤。二、实验内容利用软件编制一维电测深反演程序，并对最终的成果图进行理解与解释，同时对于数据采集、处理及解释的程序进一步的理解与体会。三、实验原理1.拟合核函数法基本思想：确定初始

2、模型：（层参数初值）；由实测曲线视电阻率转换函数曲线；由初值计算理论曲线；按照一定法则进行曲线对比，判断与拟合程度是否达到要求。根据拟合差确定是否进行层参数调查；调整，以改善与的拟合程度；重复③～⑤步骤，直到与曲线达到拟合要求为止。2.视电阻率转化函数基本思想与拟合核函数法相似。3.DZ法自动反演基本思想：①根据实测曲线给出初值②由计算理论曲线，计算DZ曲线（只计算DZ点坐标）①比较与，修正DZ曲线②由DZ曲线②计算新的模型参数⑤把当作，重复②③④步骤，直到与之间拟合差满足精度要求。三、实验成果图1、模型参数设置：2、运行曲线3、反演结果五、实

3、验体会通过对本次实验的学习，对一维电测深反演有了更深刻的，对VB的编写也有了一定的提升，但是由于各方面原因，在实验编写的过程当中也遇到了很多不可避免的问题，比如说，对反演运算过程的一些认识和程序的编写，最后通过请教其他同学才得以完成，总的来说，这次实验的收获还是不小。附加程序：一维电测深反演VB程序PubliczbxAsIntegerPublicnAsIntegerPublicxps()AsDoublePublicnj1AsIntegerPublicnj2AsIntegerPublicppsz()AsDoublePublicppsf()AsDo

4、ubleForm1部分Dimx11AsDoubleDimx22AsDoubleDimy11AsDoubleDimy22AsDoublePrivateSubCommand1_Click()DimnnAsIntegerDimhh()AsDoubleDimpp()AsDoubleDimzpp()AsDoubleDimzhh()AsDoublenn=Val(Text1.Text)ReDimhh(nn-1)AsDoubleReDimpp(nn)AsDoubleReDimzpp(nn)AsDoubleReDimzhh(nn-1)AsDoubless1=Te

5、xt2.Textss1=ss1+""ss2=Text3.Textss2=ss2+""Fori=1Tonn-1ss1=LTrim(ss1)hh(i)=Val(Left(ss1,InStr(ss1,"")-1))ss1=Mid(ss1,InStr(ss1,"")+1)NextiFori=1Tonnss2=LTrim(ss2)pp(i)=Val(Left(ss2,InStr(ss2,"")-1))ss2=Mid(ss2,InStr(ss2,"")+1)Nextiznxs=Val(Text4.Text)ddcshu=Val(Text5.Text)yyy

6、=9Text6.Text=""DimssAsStringDimccf(9)AsDoubleDimccz(18)AsDoublefname="滤波系数施伦贝尔反演滤波系数"&".txt"OpenApp.Path&fnameForInputAs#1Fori=9To1Step-1Input#1,ccf(i)ccf(i)=ccf(i)*10^(-4)NextiFori=0To18Input#1,ccz(i)ccz(i)=ccz(i)*10^(-4)NextiClose#1Dimztt()AsDoubleReDimztt(nj1Tonj2)AsDou

7、bleForj=nj1Tonj2Fori=-9To18If(j-i)>=0Thenls1=ppsz(j-i)Elsels1=ppsf(Abs(j-i))EndIfIfi>=0Thenls2=ccz(i)Elsels2=ccf(Abs(i))EndIfztt(j)=ztt(j)+ls1*ls2NextiNextjFori=nj1Tonj2-1Picture1.Line(js(Exp(i*Log(10)/6)),js(ztt(i)))-(js(Exp((i+1)*Log(10)/6)),js(ztt(i+1))),RGB(0,0,10)NextiI

8、fddcshu=0ThenDimhtt()AsDoubleReDimhtt(nj1Tonj2)AsDoubleFori=nj1Tonj2nmta=Ex