一种针对RMT数据的双参数反演算法制造技术

本发明专利技术公开了一种针对RMT数据的双参数反演算法，包括以下步骤：步骤1、读取测点的位置、测量频率以及各个测点的观测数据；构建N维观测数据向量d；步骤2、为待反演的模型的两个参数赋初值；步骤3、根据测点的位置和测量频率，计算由两个参数共同作用得到的正演响应；步骤4、计算数据拟合差；步骤5、若数据拟合差大于目标拟合差，则进行第k次双参数反演迭代；否则反演结束；进行第k次双参数反演迭代的步骤为：首先计算正演响应对两个模型参数的灵敏度矩阵，然后求解双参数反演迭代方程组，得到两个模型参数的修正量；最后计算新的模型参数；步骤6、令迭代次数加1，返回步骤3。本发明专利技术能够降低数据拟合偏差，提高反演准确性。

A double parameter inversion algorithm for RMT data

The invention discloses a method for double parameter inversion algorithm of RMT data, which comprises the following steps: 1, read data points, and measuring the frequency of each measuring point; to construct N dimensional data vector D; step 2, the two parameter inversion model for given initial value; step 3 according to the position of the measuring points and measuring the frequency response, is calculated by the two parameters of speech interaction; step 4, the calculated data fitting difference; step 5, if the difference is greater than the target data fitting fitting difference, then a K double parameter inversion iteration; otherwise the end of the K inversion; iteration inversion two parameter steps: first calculate the forward response sensitivity matrix of two model parameters, then solving iterative double parameter inversion equations, correction two model parameters; the final calculation of new model parameters; step Suddenly 6, the number of iterations plus 1, return to step 3. The invention can reduce the deviation of data fitting and improve the accuracy of inversion.

【技术实现步骤摘要】
一种针对RMT数据的双参数反演算法
本专利技术涉及一种应用地球物理领域的高频电磁数据反演算法，特别涉及射频大地电磁法数据的反演。
技术介绍
目前，各类电磁资料的反演均是以寻求满足目标拟合差的地下介质电阻率分布为目的，忽略了地下介电常数差异对观测数据的影响。然而，实际中观测数据是由地下电阻率和介电常数共同决定的，因此采用传统的单参数电阻率反演来拟合观测数据是存在一定偏差的。尤其对于射频大地电磁法(RadioMagnetotelluric,RMT)数据，由介电常数所引起的波动场在总场中的占比可达20％以上，这种情况下，仅通过电阻率参数的反演来拟合目标函数势必引起一定偏差，降低资料解释的准确性。因此，有必要开发一种综合考虑地下电阻率和介电常数的双参数同步反演算法。
技术实现思路
本专利技术所解决的技术问题是，针对现有技术的不足，提供一种双参数同步反演算法，本专利技术能够同时得到地下电阻率和介电常数分布，降低了数据拟合偏差，提高了反演准确性。本专利技术的技术方案为：一种针对RMT数据的双参数反演算法，包括以下步骤：步骤1、读取测点的位置、测量频率以及各个测点的观测数据；各个测点的观测数据构成N维本文档来自技高网...

【技术保护点】
一种针对RMT数据的双参数反演算法，其特征在于，包括以下步骤：步骤1、读取测点的位置、测量频率以及各个测点的观测数据；各个测点的观测数据构成N维观测数据向量d＝(d1,d2,…,dN)

【技术特征摘要】
1.一种针对RMT数据的双参数反演算法，其特征在于，包括以下步骤：步骤1、读取测点的位置、测量频率以及各个测点的观测数据；各个测点的观测数据构成N维观测数据向量d＝(d1,d2,…,dN)T，其中di为第i个测点的观测数据，N为测点的个数；步骤2、为待反演的模型的两个参数m1(k)和m2(k)赋初值；k为迭代次数，其初始值为1；步骤3、根据测点的位置和测量频率，计算由m1(k)和m2(k)共同作用得到的正演响应F(m1(k),m2(k))，F(m1(k),m2(k))为N维向量，F(m1(k),m2(k))＝(f1,f2,…,fN)T，其中fi为第i个测点的正演响应值；步骤4、计算数据拟合差RMSk，计算公式为步骤5、判断RMSk是否大于目标拟合差(目标拟合差一般取为1)，若是，则进行第k次双参数反演迭代；否则退出迭代，反演结束；进行第k次双参数反演迭代的步骤为：(a)计算第k次迭代时的正演响应对两个模型参数的灵敏度矩阵J1(k)和J2(k)：(b)求解如下双参数反演迭代方程组，得到第k次迭代的两个模型参数的修正量Δm1(k)和Δm2(k)；

【专利技术属性】
技术研发人员：原源，汤井田，周聪，任政勇，
申请(专利权)人：中南大学，
类型：发明
国别省市：湖南,43