WinGlink处理MT数据说明书

1、WinGlink处理MT数据说明书WinGlink是一款意大利GEOSYSTE邮司开发的综合解释平台，主要侧重丁 大地电磁的建模和数据处理，能结合重磁资料和大地电磁的联合处理解释，并实现了钻孔等资料的约束反演。该系统实现了大地电磁的一、二维反演，并可进行 三维大地电磁的模拟，其正演采用了 Randy Mackie的二维交错采样有限差分迭 代算法程序。本文针对WinGlink处理大地电磁测深数据作详细的讲解。1、使用WinGlink软件前，需将电脑转换至英文系统下：开始菜单”一控制面板”一时钟、语言和区域”一区域和语言”一管理点击 更改系统区域设置”后，出现下图所示的语言选择框，选择任意英语国

2、家，并重新启动计算机。2、重新启动后，点击图标运行软件，出现下图中的对话框：WinGUnk其中1号框图表示创建一个新的数据库；2号框图表示打开一个已有的数 据库，3号框图则表示使用已购买的工具箱。对丁初学者或该软件为破解版的情况下，一般选择1号或者2号框图按钮。本文中点击“creat newdatabase”，然后 点击“0椭钮。3、点击“。顾钮后，会弹出一存储对话框，选择存储的文件夹并编辑数据 库名称。¥ New WinGLink DatabaseW头计律机汁Hi曜I卜4、新建数据库后，则需填写新建的数据库信息（Database propertied ,其 中共三项，分别是： Ge

3、nerak Projected Coords Geographic Coords第一项：General （常规信息），其中包含地区名称、所处半球，平面坐标单位、垂直坐标的单位。Propertieseri Projected Doords Geosrapliic Coordw |Area NameXr1 Ccordnetes Unit?|Tit»n 刑区名称XY坐标的单位,般使用公里（Kilometer）* & I Elevation UnitsL I加亦垂直方向（高程）上的单位，默汰为米（meters ）CancelNext» 影，处理MT数据一般默认使用该投影，故

4、无需进行改变。中间的True Origin（坐 标原点），Latitude框无法改变，而Longitude框的数值则需改成测区范围的任意 经度，并以度分秒的形式，例如改写成：107°00'00.0000'。'下方的Datum （基准）， 默认采用 WGS 1984即 World Geodetic System 1984是GPE球定位系统使用 而建立的坐标系统，无需进行改变。第二项：ProjectedCoords （投影坐标），对话框中的UTM即Universal Transverse Mercator表示通用横轴墨卡托投0第三项：Geographic Coor

5、ds（地理坐标），地理坐标让选用 WGS 1984坐标系 统，无需进行更改5、点击上图的“OKB,出现以下的界面:Rle Project Settings HelpD artabase PrcpertieAreaHarne:TfctclanHemispfiere:Not thXYUrilxKig metersElev-aton lints: McteisFTojeCIm1 ransvBtse Mere*rGrid:UTMrrMfeOiigriiLahbjde:cnomMwrLonabude:lOTWoamoErGnd Ggd顽es d Tiue 0«ignN a thing0.000

6、 m湖i®SOQOQaOOQmFoe0.999SM elnc C&ordnalesSpheroid:WG8 1994Ddum:WGS1984GBCtgrvhic CoordnaABS!Sphefori:WGS 1994DdurncVA3S19046、建立数据库之后，若需处理数据，则还需先建立一个工程，点击“Project ”一 “ Ne讨一山Snpleje'ct,的对话框后填写相应的信息:上图中，Summary选项中，从上至下，分别为日期、工程名、数据类型、station legend、位置、公司、承包商。其中的 Datatype 一栏相当重要，处理大地电磁 数据则需

7、选择“MT magnetotellurics项。”而对丁 “Summa纺的“Inf节U“Attachments*填可不填。信息填写完成后，点击 “OK7、创建好工程后，主界面中则会增加相应的功能图标:上图中，中间添加了一个新的工程，右侧所标注的图标则是进行 MT数据处 理和反演的功能。第1项“MapS功能为划分剖面，选择测点。第2项“Sounding” 则是修改测线上测点的数据。第3项“-Sections,”功能是绘制拟断面。第4项“2D Inversion贝是进行二维反演。第5项"3D modelin州是建立三维模型。由丁仅仅创立了一个工程，工程里却没有数据，故需从文件中导入数据。

8、8、导入数据。点击主界面左上角“File " Import。”弹出“Import对话框。选择“External Files从'外部导入数据文件。然后点击Next,”进入下一界面Inan Existing Project。若是之前未建立工程，则可选择第一个选项“Ia NewProject。”选择后，点击"Next”弹出下图所示的界面:上图对话框中，两个选项，由丁之前已经建立了一个工程，则该选择第二个,上图中，即在数据文件中导入数据，红色标记表示，选取全部的数据。若选取其 中的部分数据，可按住键盘上的“Ctrl,'然后点击鼠标左键选取。选择好后，点击“NextK

9、钮，进入下一界面。这是会出现一个“Warning的对话框，直接点击“NexfP 可。Next!”下一个界面则比较重要点击“Next!”出现下图所示的界面上图中，无需更改，直接点击Plot Parameters选项中，要选取 “No Recalculation不需要重新计算。而 在“Rotate tO勺选项中，由丁没有进行电性主轴分析，不知道旋转角度。故在这 里将"Principal Axi核为” “User Defined角度仍未0.0° 。设置好后，点击按钮 “Star开始读入数据。Import MT DataImport done!72 Statiom read72 S

10、tations appendedFinish< Back.数据读入完成!9、选取测线与测点。之前录入 MT数据的格式为.edi ”，其中包含了该点 的GPS坐标，高程等一系列信息。点击主界面右端“Maps®标。点击图标后，会弹出一对话框，如下:点击上图的Open选项后，弹出下图的对话框:设置完成后，弹出Maps的功能界面:划分测线：对MT数据进行反演处理需划分测线并选取测点，故点击“Profile ”一 "Add profile tfei南在图中选择一起点，点击左点，确定终点后, 点击右键！紧接着会弹出一对话框，需填写该测线的信息。点击OK后，在界面上则会显示一条标注

11、了名称的测线。划分测线后，还需选择测线周围的测点进行反演。点击左上角“Profile ” oto station model。图形界面则会发生变化。（注意：划分测线时要 注意测线的方向，以上图为例，划分测线时，是以 A为起点，B为终点。因此在 反演时，左端为A,右端为B）选取测点：选取测点的方式有三种，如下图所表示:选取全部测点，点击"profile- "'save”然后关闭该界面。10、处理测点数据划分好测线，并选取测点后，则需对数据进行相应的处理，如剔除测点中的 飞点”、删除数据质量差的测点、平滑曲线等。点击下图中的“Soundings”弹出下图中的对话框本文中

12、选择“Open stations from List其余圭选项较为简单，学生可自行摸 索。点击OK选择要处理测测点。（没有全选功能，只能一个个点击）。选择完后，点击Operf界面则会显示出选中测点的曲线图： 回 S3f f n£ i <13I urtrn店E3SESJ®蚌g Ua£ mfr-1> t2 j.f *岫 a aanaq .b u uq 申ift5 »"P*w« |bK：KilMfl> =M-X产渤1 .即的一 * f钟T Mm-E-ET-1 0£_选中测点的曲线= y可由L卢击墨大化，开始编辑处

13、理数据点击“Edi暗，界面右侧增加相应编辑功能按钮!"J!t I |hI夏或牛还；Tzy倾子信息1?'T-1 11'Sm 000 貌1声1 陶 Nc»e日di g阳中上图右侧以上自下分别是：(1) Mask X¥删除/添加界面中红色点。(2) MaskYx删除/添加界面中蓝色点。(3) Mask XY-YX删除/添加界面上的点。(相对应红蓝点同时删除或添加)(4) Mask Tipper:编辑倾子信息。(5) Swap:(6) Strip:(7) S-Shift:(8) SS Reset(9) Shift：整体上下移动曲线(10) Smooth:平

14、'滑，可设定平滑参数(11) + 360°:(12) - 360 :。(13) Cancel:取消(14) HF mod(15) Close：关闭该点数据编辑编辑好测点数据后，点击“Close”选择保存后退出11、二维反演:点击后，出现的对话框中，一般不带地形反演。Profile and Model SelectionProfilesLiiel选择的测线Models<Nsw »建一个新模型!r Uss Topograph带地形反演F" Deale Marina Model K.I Cancel学1-1 时 nifff nY YW=Z .ULU LJ

15、标*嗣。宇户号 岸叩n i in1 uu nil 4 j umin m =fr* *+* 岸 h * 言 W" +?+国 ,岸日 哗 '叩 点击OK!7淼-Qi Hh就J££ i4 S t. lElSffiJl P 35t对丁剖分网格，一般的要求是：(1) 两相邻测点至少有一条垂直线，若两点问有多条，则垂直网格线间距 相同！测线两侧的的垂直网格线，间距则按 1-2倍递增。测线两长度 一般为测线长度2-3倍。(2) 水平网格线一般按照趋肤深度来剖分，间距则可按1-2倍递增，剖分 的深度一般为4倍的反演深度。功能：(1) 添加网格线：添加行：点击 “Edit

16、” 一 "AddRumn ”一 “row”添加歹0:点击 “Edit” 一 "Add Row/column” 一 "column”(2) 删除网格线：删除行：点击 “Edit” 一 "Delete Row/ColumrowV “删除歹0:点击 “Edit” 一 "Delete Row/Column ” 一 "Column(3) 移动网格线：移动行：“Edit” 一 "Cells Height/Width ”一 "Row Height移动列：“Edit” 一 S Celght/Width ”一 "Colum

17、n Width”(4) 设置色棒：点击 “RangeS'Range Editor- Profile: BLinel" - Model： « New»T ablesB1522W1CI24312Ranges|3 deoimalsSpitDeleteReset to defeukNormatzeRainbow-T emplatesJ LoadJ SaveJ DeleteBound ValueBlue Cyan Red Green Yellow Magenla 囱剪 Light Gfisy rangB R* Red and Black Green end Blue

18、 Water, Blue. Red Blu&, Orange Blue. While. Red Blue, Magenla RanG心nEdit |匚岫屈-UnOKCancela'-Jj Reverse .使得幻色太作低阻蓝色袁小局阻:点击Automatic,弹出另一界面:设置好后，点击OK!(5)设置边界(6)保存模型：点击：“File Fave model ”弹出对话框:其中命名以按参数命名。以上图为例，tm表示反演模式，tau为光滑 系数，a为圆滑因子。(7)反演:上述工作完成之后，可开始进行反演。先退出编辑：点击“Edit”-“Exit Edit Mesh。退出之后才能

19、在菜单栏里出现“Inversio施项。设置反演系数：点击"Inversion ”一 "Setting由如下对话框:该参数设置非常重要!MainParms本文中，仅选择TM模式的反演。反演误差可自主设定！建议值为Parameters SettingData SelectSlatic ShiftMain FarmsFixed Farms Cgata Select |D ala Err or $Smooth Inv.Error FlooiSharp Inv.Min Frequency最小的频率TE data|0.01000TMdataf0.01000HZ data|a 01000

20、TE: ttTH: tt DecadesHZ: tt Decadec4 CO4W级数ITooUse Stetian DataUse Interpolated Datatt Freqj in Decade |5Use Smoothed Curves if existing非常重嬖，对于M【处理，仅仅点要Use Station DataM，取消,TUse Interpolate da :a 以及，"Uwp Smoothed Curves if existing"OK |CancelSmooth InvOKCancelParameters SettingStatic Shift

21、Fixed Parm | Data Errors | Enor FloorMain Parnms Data Select. Tfiooth Sharp Inv.f* Solve for the smoothes! model r 俄r厂 Solve fa the smoothest variationsaway frorn the apriori modelRegularization Laplacian'* Standard grid Laplacian operator'Uriform grid Laplacian operator(the uniform grid Lap

22、lacian Operator Mil likelyresult ir a much smoothei model 1R egulaiLiioM*F Minimize integral of Igrad(m)r2=_Minimize integral of |Laplacian(m)|weighting function smooth.par file)l?TWeightirig eriabieO Horizortal alpha factor Beta factor for weighting Min. block dirnenjianj (m|黄向上的光|纵向上的*滑500.0 H | 5

23、00.0 VFocussing叩| 0.00010 exponent |1.250For L2 smoolh inversion： epsilon = U, exponent = 2A!Jaufor Smoothing operatgp 3.00Oiutput sensitivity map 厂 S ave intermediate models光滑因子，越大越平滑。Parameters SettingError FloorMain Farms | Data Select. | Smooth Inv. | Sharp Inv. | Static Shift | Fined Farms | Da

24、taEnorc (£ror FbalError Floor 阂TE Rho |5.000TM Rho |5W0TE Phase |5.000 TM Phase)5.000两种模式下视阻率和相位的误差lyz Data, absolute |0.10000Note: Phase errors are entered in % equivalent rho error. A1 泛 rho error equivalent to a 0.29 degree absolute error in phaseCancel根据反演要求自行设定误差级数Parameters SettingData Er

25、rorsOKCancelMain Farms | Data Select. | Smooth Inv. Sharp Inv.Static Shift | Fixed Farms Erro| Error Floor |10.000TM Rho 10.000TM P嵌e5.000Default Standard Deviation Error IE Rho IE Phase 5.000For equal weighting of rho and phase in the inversion, the errors for phase input here should be half those

26、input for apparent resistivities.Tyz Data | IQDOOOONote: Phase Birors are cntefBd n % eqjivaleri rho error. A 1X rho error i$ equivalent to e 0.25 degree absolute error in pha沿/ Use Data Errors if existingFixedParmsMain Farms | Data Select. | Smooth Inv. Sharp Inv Static Shift ijieJ Par?j| DataEnor Error Floor'Fixed Paramelers10000000000T孙 for clamping fixed paramete