基于动态更新加权矩阵的密度模型求取方法技术

技术编号:41479682 阅读:23 留言:0更新日期:2024-05-30 14:30
本发明专利技术公开了基于动态更新加权矩阵的密度模型求取方法,解决了现有技术中密度模型求取不准确和分辨率不高的问题,具体步骤包括:利用互相关系数为每个单元体分配初始权重,即构建初始加权矩阵;通过求解第一最小二乘问题,得出每个单元体的密度与权重之比;根据每个单元体的密度与权重之比,求出第k次迭代的密度向量;通过求解第二最小二乘问题,得出每个单元体的新权重,并构建新加权矩阵;以新加权矩阵代替初始加权矩阵,循环执行迭代。在迭代过程中,实时监测相邻两次迭代的密度向量之差的2范数,若该范数小于预设的阈值,即终止迭代,并输出最后一次迭代的密度向量作为最终的密度模型。本发明专利技术具有逻辑简单、高效可靠和高分辨率等优点。

【技术实现步骤摘要】

本专利技术涉及密度模型求解,尤其是基于动态更新加权矩阵的密度模型求取方法


技术介绍

1、在地球物理勘探与地球科学研究领域,重力反演是一种常用且重要的技术,它通过分析重力场的变化来推断地下密度异常体的特征。其中,经典的反演方法是将地下空间离散为一系列紧密排列的单元体,进而求解离散不适定问题来估算每个单元体的密度。在这一过程中,为每个单元体分配适当的权重,即引入合理的加权矩阵,无疑是提高反演分辨率的关键。在此,需要解释的是,地下单元体是指将目标地区的地下空间划分为许多个紧密排列的小区间。

2、然而,传统方法在为每个单元体分配权重时往往依赖于某种粗略的估计,并且在整个反演过程中这些权重保持恒定不变,这在很大程度上限制了反演分辨率的进一步提升。例如,li和oldenburg(1998)提出基于单元体深度的权重分配策略,portniaguine和zhdanov(2002)通过灵敏度分析来确定权重等。这些方法的本质都是为了克服重力的“屈肤效应”而施加的深度加权,旨在提高对深部异常的反映能力。然而,仅利用深度进行加权得出的反演结果分辨率较低,并且很难恢复出本文档来自技高网...

【技术保护点】

1.基于动态更新加权矩阵的密度模型求取方法,其特征在于,包括以下步骤:

2.根据权利要求1所述的基于动态更新加权矩阵的密度模型求取方法,其特征在于,所述互相关系数的值域为(-1,1),且所述初始加权矩阵的表达式为:

3.根据权利要求2所述的基于动态更新加权矩阵的密度模型求取方法,其特征在于,结合第一加权矩阵,求解第一最小二乘问题,得出每个单元体当前次迭代的密度与权重之比,包括以下步骤:

4. 根据权利要求3所述的基于动态更新加权矩阵的密度模型求取方法,其特征在于,根据每个单元体当前次迭代的密度与权重之比,求出当前次迭代的密度向量,其表达为:

...

【技术特征摘要】

1.基于动态更新加权矩阵的密度模型求取方法,其特征在于,包括以下步骤:

2.根据权利要求1所述的基于动态更新加权矩阵的密度模型求取方法,其特征在于,所述互相关系数的值域为(-1,1),且所述初始加权矩阵的表达式为:

3.根据权利要求2所述的基于动态更新加权矩阵的密度模型求取方法,其特征在于,结合第一加权矩阵,求解第一最小二乘问题,得出每个单元体当前次迭代的密度与权重...

【专利技术属性】
技术研发人员:梁生贤代堰锫王昌南唐高林李同柱廖国忠李华郭镜焦彦杰杨剑李富夏时斌李德伟
申请(专利权)人:中国地质调查局成都地质调查中心西南地质科技创新中心
类型:发明
国别省市:

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1