[发明专利]一种新的处理CSAMT资料的方法

权利要求书

1.一种新的处理CSAMT资料的方法，其特征在于：包括：在CSAMT的电场和磁场分量的数学表达式中，电场幅值的表达式是包含有地下介质的电阻率信息的：确定电偶源一维层状介质CSAMT电场分量数学表达式：

式中： 表示与源同方向的水平电场分量， 为偶极源的电偶极矩，，为供电电流，为电偶极长度，为自由空间磁导率，为角频率，，，为第一层电阻率，为电偶源中点到接收点的矢径和电偶源方向的夹角，为收发距，和分别为零阶和一阶贝塞尔函数，和分别为联系各层电性的函数。

2.根据权利要求1所述的一种新的处理CSAMT资料的方法其特征在于：确定基于光滑模型约束的总目标函数：

式中：为总目标函数，为模型水平粗糙度，为模型横向粗糙度矩阵，定义为模型参数的一阶导数，为正则化因子，为期望的拟合差，为数据残差，为模型修改量，为Ex对电导率一阶偏导数矩阵，为数据加权矩阵。

3.根据权利要求1所述的一种新的处理CSAMT资料的方法其特征在于：对模型参数进行约束：

对模型进行分层剖分，假设模型有M层，测量得到的CSAMT单分量电场幅值数据个数为N个，为了有效地减少多解性和确保参数的非负性，取：

式中：为实测Ex数据，为模型正演响应数据，为模型第层的预测的电导率值，为第层的当前电导率值，为第个观测数据的标准差。

4.根据权利要求1所述的一种新的处理CSAMT资料的方法其特征在于：求偏导数矩阵：

偏导数矩阵的求解使用差商的方法，即：

。

5.根据权利要求1所述的一种新的处理CSAMT资料的方法其特征在于：对总目标函数求极小，并对初始模型进行修改迭代反演得到反演结果：

得到：

解上述方程组，得到，对初始模型进行修改迭代反演得到反演结果。