一种基于结构约束的归一化重磁电震联合反演方法技术

本发明专利技术公开了一种基于结构约束的归一化重磁电震联合反演方法，包括如下步骤：S1：建立归一化联合反演目标函数；S2：建立初始模型；S3：正演计算和雅克比矩阵求取；S4：计算归一化交叉梯度偏导数和拉格朗日算子；S5：计算获得迭代模型；S6：对获得的迭代模型进行正演计算，并与观测数据进行拟合差求取，如果拟合差增大、达到最大迭代次数或者拟合差达到阈值，则停止循环反演，输出归一化联合反演模型，否则执行S3。该基于结构约束的归一化重磁电震联合反演方法克服了单一方法的局限性、减小了地球物理反演多解性、不依赖于岩石物性关系，同时还可以有效的避免由于数量级相差较大的物性参数直接耦合反演所导致的错误反演结果。

【技术实现步骤摘要】
一种基于结构约束的归一化重磁电震联合反演方法
本专利技术涉及重力、磁法、大地电磁和地震初至波走时四种地球物理方法同时处理的领域，特别提供了一种基于结构约束的归一化重磁电震联合反演方法。
技术介绍
根据岩石物性参数的不同，地球物理方法可分为重力、磁法、地震和电法等，重力和磁法具有很好的横向分辨能力，对地质构造分布划分起到了重要的作用，但纵向分辨能力较差，且勘探深度浅；大地电磁测深法(MT)的勘探深度与频率有关，由于频率范围丰富，其勘探深度达到几十或者上百公里，同时具有不受高阻层屏蔽，对良导体反映灵敏等优点，因此在研究地球深部构造等方面具有特殊的优势，但是横向分辨能力较差；对于地震初至波走时方法，具有较高的纵向分辨能力，在划分地下结构时分辨能力敏感，但是在盐丘、高角度断裂存在时容易形成盲区，并且初至波传播深度较浅，导致勘探深度较浅。以上每种单一方法都存在一定的缺陷，每种方法都只能从单一角度来评价地下介质的岩石物理特性，单凭一种地球物理方法很难准确勘探地下结构，因此，为了得到更准确的地下信息，利用多种地球物理参数对同一地下地质体进行综合解释已经成为当今发展的趋势。联合反演是综合地球物理方法本文档来自技高网...

【技术保护点】
1.一种基于结构约束的归一化重磁电震联合反演方法，其特征在于，包括如下步骤：S1：建立归一化联合反演目标函数，所述联合反演目标函数如式(1)：

【技术特征摘要】
1.一种基于结构约束的归一化重磁电震联合反演方法，其特征在于，包括如下步骤：S1：建立归一化联合反演目标函数，所述联合反演目标函数如式(1)：交叉梯度约束条件：τ(m)＝0其中：Cd＝diag[Cd1,Cd2,Cd3,Cd4]，Cm＝diag[Cm1,Cm2,Cm3,Cm4]，α＝[α1,α2,α3,α4]，其中，m为模型参数，m1,m2,m3,m4分别为电阻率、速度、密度和磁化率模型参数；m0为先验模型参数，m01,m02,m03,m04分别为电阻率、速度、密度和磁化率先验模型参数；d为观测数据，d1,d2,d3,d4分别为视电阻率、地震走时、重力异常和磁异常；Cd为观测数据d的数据协方差矩阵，Cd1,Cd2,Cd3,Cd4分别为视电阻率、地震走时、重力异常和磁异常的数据协方差矩阵；Cm为模型参数m的模型协方差矩阵，Cm1,Cm2,Cm3,Cm4分别为电阻率、速度、密度和磁化率的模型协方差矩阵；α为阻尼因子，α1,α2,α3,α4分别为大地电磁、地震初至波走时、重力和磁法的阻尼因子，f(m)表示正演响应，f1(m),f2(m),f3(m),f4(m)分别为大地电磁、地震初至波走时、重力和磁法的正演响应，为梯度，τ为交叉梯度，τij为大地电磁、地震初至波走时、重力和磁法两两之间的交叉梯度；为数据拟合项、为模型平滑约束项；κ1,κ2,κ3,κ4分别为电阻率、速度、密度和磁化率归一化因子；T和-1分别为矩阵转置和矩阵求逆；S2：建立电阻率、速度、密度和磁化率初始模型：在笛卡尔坐标系下沿x，z二个坐标轴分别将初始模型空间划分成为Nx，Nz个小的矩形网格单元，间距为Δx(i...

【专利技术属性】
技术研发人员：张镕哲，李桐林，邓馨卉，
申请(专利权)人：吉林大学，
类型：发明
国别省市：吉林,22