一种自适应的震源区重力密度变化实时反演计算方法技术

本发明专利技术提供了一种自适应的震源区重力密度变化实时反演计算方法，通过在震源区野外实时观测到的重力动态变化结果数据，自动、实时的反演震源区重力密度变化分布特征，得到初步的密度反演结果；根据结果修正或补充测量点位，较好的得到重力密度结构分布，实现了野外数据标准化、重力数据自动化改正处理和自适应的重力密度反演的功能。本发明专利技术提高了重力观测数据的处理效率，在野外边观测边动态反演，实时了解地下结构变化；本发明专利技术采用了高度集成的数据处理方法，使重力数据处理从分散、独立变为统一化；本发明专利技术自动分片提取需要的海量地形数据参与计算，计算过程不需要人工干预，降低了开发和维护成本。

【技术实现步骤摘要】
一种自适应的震源区重力密度变化实时反演计算方法
本专利技术属于重力密度反演
，具体涉及一种自适应的震源区重力密度变化实时反演计算方法。
技术介绍
地震孕育、发生以及发展历程，是现阶段地球物理学科当中一个重要的研究领域。而地震的孕育、发生和发展的整个过程和地球内部构造有着极为密切且重要的联系。研究地壳密度结构是了解地球内部构造的有效方法之一。通过研究地壳密度结构，可探索研究地下岩体物质的密度分布或异常岩体的空间位置、大小、形状以及断裂构造等地震发生的主要构造环境。重力反演是获取地壳密度结构较为有效的方法之一。重力反演是以地壳中不同的密度体之间的密度差为物性基础，通过重力异常的变化来研究地下空间的地址构造特征及其变化，主要用于地球深部结构的探测，区域地质构造单元的划分，沉积盆地圈定以及油气资源勘查等领域。通常，密度变化反演需要野外采集完数据后，经过室内处理，得到改正后的重力变化结果，再根据具体区域设置参数，采用相关反演算法进行反演，得到密度分布结果。这种传统的处理方式不便于在野外实时优化观测网络，如果在某些重要区域观测数据不够本文档来自技高网...

【技术保护点】
1.一种自适应的震源区重力密度变化实时反演计算方法，其特征在于：包括以下步骤：/nS1：建立地形数据库和重力场数据库；/nS2：调用读取函数读取并拼接数据库中的数据，生成地形数据文件；/nS3：根据地形数据和重力场数据进行重力改正，计算重力异常场；/nS4：根据重力异常场进行密度结构反演，具体步骤为：/nS41：将地下块体进行3D网格划分为m个棱柱体单元，得到初始的空间密度分布M0；/nS42：选择第k个棱柱体单元填充密度差值，根据重力异常场设置密度差范围，通过反演得到由k个棱柱体单元组成的异常体模型的空间密度分布Mk；/nS43：根据观测点的重力影响正演得到观测异常；/nS44：根据重力异常...

【技术特征摘要】
1.一种自适应的震源区重力密度变化实时反演计算方法，其特征在于：包括以下步骤：
S1：建立地形数据库和重力场数据库；
S2：调用读取函数读取并拼接数据库中的数据，生成地形数据文件；
S3：根据地形数据和重力场数据进行重力改正，计算重力异常场；
S4：根据重力异常场进行密度结构反演，具体步骤为：
S41：将地下块体进行3D网格划分为m个棱柱体单元，得到初始的空间密度分布M0；
S42：选择第k个棱柱体单元填充密度差值，根据重力异常场设置密度差范围，通过反演得到由k个棱柱体单元组成的异常体模型的空间密度分布Mk；
S43：根据观测点的重力影响正演得到观测异常；
S44：根据重力异常场设置密度差范围，通过重力异常场和观测异常反演得到异常体模型的空间密度分布M(k+1)；
S45：根据密度差和残差判断空间密度分布M(k+1)的误差是否满足预设要求；若不满足则合并空间密度分布M(k+1)与空间密度分布Mk，执行步骤S42；若满足要求则将当前空间密度分布作为最终的密度异常体M。


2.根据权利要求1所述的一种自适应的震源区重力密度变化实时反演计算方法，其特征在于：所述的步骤S3中，重力异常场包括自由空气异常、布格重力异常和均衡重力异常。


3.根据权利要求1所述的一种自适应的震源区重力密度变化实时反演计算方法，其特征在于：所述的步骤S42中，具体步骤为：
S421：设单位密度下重力观测点pi处产生的重力影响为Δgij，密度差值为Δρj，则搜索到第k步时由k个棱柱体单元组成的异常体模型的重力异常为：



S422：对第k+1个棱柱体单元j(j≠1,2,…,k)填充密度差值Δρj并聚合到模型构成新模型。


4.根据权利要求3所述的一种自适应的震源区重力密度变化实时反演计算方法，其特征在于：所述的步骤S43中，具体步骤为：
S431：设万有引力常量为G，重力观测点的pi坐标为(xi,yi,zi)，棱柱体单元的坐标为则单位密度下重力观测点pi处产生的重力影响Δgij为：



S432：单位密度下所有棱柱体单元在重力观测点pi处产生的重力值为：





5.根据权利要求4所述的一种自适应的震源区重力密度变化实时反演计算方法，其特征在于：所述的步骤S44中，具体步骤为：
调整未知的比例因子fk拟合新模型的重力异常与观测异常Δgi：






设模型的协方差矩阵QM由非空元素的对角归一化矩阵给出，协...

【专利技术属性】
技术研发人员：杨光亮，谈洪波，王嘉沛，周怀斌，申重阳，
申请(专利权)人：湖北省地震局中国地震局地震研究所，
类型：发明
国别省市：湖北;42