直流电测深反演报告.docVIP

PAGE 电阻率测深反演报告 姓 名： 黄鋆 学 号： 20081001550 班 级：061084-11 学 院: 地空学院 指导教师： 师学明、张文波 2011.5.25 目录 实验要求…………………………………………………3 实验原理…………………………………………………3 程序设计流程以及程序源代码…………………………3 实验过程以及程序运行说明……………………………10 实验不足与感想…………………………………………13 参考文献…………………………………………………14 实验要求 1.编制一个可视化反演程序，有可视化的对话窗口。 2.可用对话窗口实测深数据，然后读取数据马上显示实测数据 3.能有可视化窗口修改地电模型参数，计算实测值，计算相对误差。 实验原理 原理： 正演拟合法原理是根据实测视电阻率测深曲线特征确定预测岩层的层数，并试探地给一组层参数，用这组参数计算一条理论视电阻率测深曲线，将理论曲线与实测曲线进行对比，当两者偏差较大时，根据解正问题时所掌握的测深曲线特征，对模型进行修改，重算其理论曲线，再次进行对比，……，如此反复进行，直到两者偏差达到精度要求范围为止，最后的理论模型就可作为所求的反演问题的近似解。 实现步骤： （1）输入实测电阻率曲线数据； （2）根据实测电阻率测深曲线确定预测模型的岩层层数n； （3）根据测深曲线的特征，初步确定预测模型的层参数,； （4）递推计算电阻率转换函数Ti； （5）用数值滤波计算模型的视电阻率曲线； （6）对比和曲线，计算拟合误差； （7）当拟合误差满足精度要求（5%），输出预测模型参数并保存，反演结束。否则，继续修改预测模型的层参数，知道满足精度要求。 程序设计流程以及程序源代码 程序设计流程： 按照原理中的实现步骤，程序流程图如下： 输入观测 输入观测数据 输出反演结果ε=||-||2修改初始模型，记为正演计算初始模型 输出反演结果 ε=||-||2 修改初始模型，记为 正演计算 初始模型 否 是 程序源代码： function varargout = dianfa(varargin) gui_Singleton = 1; gui_State = struct(gui_Name, mfilename, ... gui_Singleton, gui_Singleton, ... gui_OpeningFcn, @dianfa_OpeningFcn, ...5 gui_OutputFcn, @dianfa_OutputFcn, ... gui_LayoutFcn, [] , ... gui_Callback, []); if nargin ischar(varargin{1}) gui_State.gui_Callback = str2func(varargin{1}); end if nargout [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:}); else gui_mainfcn(gui_State, varargin{:}); end function dianfa_OpeningFcn(hObject, eventdata, handles, varargin) handles.output = hObject; guidata(hObject, handles); function varargout = dianfa_OutputFcn(hObject, eventdata, handles) varargout{1} = handles.output; function pushbutton1_Callback(hObject, eventdata, handles) global psce global ab ab=textread(ab.txt);%读入极距AB/2