资源描述：
《【最新精品版】直流电测深反演报告》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、电阻率测深反演报告姓名：黄鋆学号：20081001550班级：061084-11学院:地空学院指导教师：师学明、张文波2011.5.2514目录一、实验要求…………………………………………………3二、实验原理…………………………………………………3三、程序设计流程以及程序源代码…………………………3四、实验过程以及程序运行说明……………………………10五、实验不足与感想…………………………………………13六、参考文献…………………………………………………1414一、实验要求1.编制一个可视化反演程序，有可视化的对话窗口
2、。2.可用对话窗口实测深数据，然后读取数据马上显示实测数据3.能有可视化窗口修改地电模型参数，计算实测值，计算相对误差。二、实验原理1.原理：正演拟合法原理是根据实测视电阻率测深曲线特征确定预测岩层的层数，并试探地给一组层参数，用这组参数计算一条理论视电阻率测深曲线，将理论曲线与实测曲线进行对比，当两者偏差较大时，根据解正问题时所掌握的测深曲线特征，对模型进行修改，重算其理论曲线，再次进行对比，……，如此反复进行，直到两者偏差达到精度要求范围为止，最后的理论模型就可作为所求的反演问题的近似解。2.实现步骤：（1）输入实
3、测电阻率曲线数据；（2）根据实测电阻率测深曲线确定预测模型的岩层层数n；（3）根据测深曲线的特征，初步确定预测模型的层参数,；（4）递推计算电阻率转换函数Ti；（5）用数值滤波计算模型的视电阻率曲线；（6）对比和曲线，计算拟合误差；（7）当拟合误差满足精度要求（<5%），输出预测模型参数并保存，反演结束。否则，继续修改预测模型的层参数，知道满足精度要求。三、程序设计流程以及程序源代码141.程序设计流程：按照原理中的实现步骤，程序流程图如下：输入观测数据输出反演结果ε=
4、
5、-
6、
7、2<<修改初始模型，记为正演计算初始模型
8、否是2.程序源代码：functionvarargout=dianfa(varargin)gui_Singleton=1;gui_State=struct('gui_Name',mfilename,...'gui_Singleton',gui_Singleton,...'gui_OpeningFcn',@dianfa_OpeningFcn,...514'gui_OutputFcn',@dianfa_OutputFcn,...'gui_LayoutFcn',[],...'gui_Callback',[]);ifnargin
9、&&ischar(varargin{1})gui_State.gui_Callback=str2func(varargin{1});endifnargout[varargout{1:nargout}]=gui_mainfcn(gui_State,varargin{:});elsegui_mainfcn(gui_State,varargin{:});endfunctiondianfa_OpeningFcn(hObject,eventdata,handles,varargin)handles.output=hObject;
10、guidata(hObject,handles);functionvarargout=dianfa_OutputFcn(hObject,eventdata,handles)varargout{1}=handles.output;functionpushbutton1_Callback(hObject,eventdata,handles)globalpsceglobalabab=textread('ab.txt');%读入极距AB/2[filename,pathname]=uigetfile({'*.txt'},'请导入
11、实测数据');str=fullfile(pathname,filename);psce=load(str);functionpushbutton2_Callback(hObject,eventdata,handles)globalpsce%实测电阻率globalabglobalflag%是否反演的标志flag=1;axes(handles.axes1);loglog(ab,psce,'r-o','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','g','Mark
12、erSize',5);gridonxlabel('AB/2(m)');ylabel('ρs(Ωm)');legend('实测数据');holdonfunctionpushbutton3_Callback(hObject,eventdata,handles)globalpsglobalpsceglobalabglobalrn14gl
