Euler反褶积法的模型实验.pptVIP

Euler反褶积法的模型实验 基本原理 实验目的：学习Euler反演方法，不同构造指数和不同窗口对反演结果的影响。 实验的主要内容：编制简单模型的正、反演计算程序，计算出反演结果，并对反演结果进行解释，写出实验报告。 实验要求： 任选一个简单地质模型，进行实验； 不同构造指数；N=0.5，1.0，1.5，2.0，3.0 不同大小的窗口(最小是3*3) 提交程序（MATLAB或者其他语言代码）和实验报告 实验步骤： 选择一个简单模型（具体公式见教科书），计算正演数据（总场和梯度值），注意:如果选择球体模型，需要做三维的反演，其余地质体只要求做二维反演即可； 选择子窗口如3*3，或者更大的子窗口；对于二维情况，选择至少七个点的数据做窗口； 在子窗口内，根据坐标、总场值和梯度值，带入欧拉方程，建立超定方程组（如果是3*3则有9个线性方程组成），然后利用最小二乘法解此方程组。进而得到源体位置（x0，y0，z0 ），即背景值B。 4. 按照一定间隔移动窗口，重复(2)、(3)步骤，直至覆盖全区。 将各种构造指数反演结果以不同的符号绘制在同一张图内。从图中查看反演场源点的聚散情况，从而选择合适的结果最会最终的结果。（一般采用统计法挑选合适的值，比如筛选法，这里不做要求，只是均值和均方差） 注意：在Matlab环境中，左除法就是最小二乘意义下 的解。“\” * * 重磁场满足欧拉方程: （x0,y0,z0）是场源中心位置坐标，B是背景值。在已知总场强度及其梯度的情况下，对于不同的构造指数N，就可以得到不同的线性方程组，其中x0，y0，z0，B是未知数。反演场源，就是解上面的方程组。