申请/专利权人：中国地质科学院

申请日：2020-03-30

公开（公告）日：2020-06-26

公开（公告）号：CN111337993A

主分类号：G01V7/06(20060101)

分类号：G01V7/06(20060101)

优先权：

专利状态码：在审-实质审查的生效

法律状态：2020.07.21#实质审查的生效;2020.06.26#公开

摘要：本发明公开了一种基于变密度变深度约束的重力密度界面反演方法，包括：对获取的研究区内重力异常数据通过二维傅里叶变换转换，在频率域中进行正反演计算，得出修正后的重力异常场反演计算递推公式；基于全球地壳模型和实测地震资料信息获取研究区内密度界面的密度、深度信息情况；以得到的已知密度、深度信息数据作为约束条件，利用得出的修正后重力异常场反演计算递推公式，在频率域中进行密度界面反演计算。本发明通过同时考虑密度和深度变化的频率域重力密度界面反演方法，计算速度快，还能方便快捷的对大量数据进行处理。且本发明以已知的先验信息为约束条件，进行密度界面反演，能极大的提高重力方法反演的准确性，与实际情况相贴近。

主权项：1.一种基于变密度变深度约束的重力密度界面反演方法，其特征在于，所述方法包括如下步骤：1获取研究区范围内的重力数据，并对重力数据信息进行有效提取，得出由密度界面所反映出的重力异常信息，再将所述重力异常信息的数据分布通过二维傅里叶变换进行转换，在频率域中进行正反演计算，得出修正后的重力异常场的反演计算递推公式： 式中，F[·]表示对括号内对应量的傅里叶变换，Δh为界面的修正量，为扰动源坐标r在xy平面上的投影，上角标b、b-1、b+1为迭代次数序号，G为万有引力常数，ρ为密度差，ζ为界面S下某点的横坐标值；Δginv为实测的布格重力异常与第b次模型正演获得的重力异常之差，即为模型的改正量的重力异常，a为波数场的模值，n为导数阶数；为第b次迭代计算出深度正演出的重力异常，为第b次的迭代向下延拓算子，Ne为迭代次数；2基于全球地壳模型Crust1.0和实测地震资料信息获取研究区范围内密度界面的密度、深度信息分布情况；3以步骤2得到的研究区密度界面已知的密度、深度信息分布数据作为约束条件，利用步骤1得出的修正后重力异常场的反演计算递推公式，在频率域中进行密度界面反演计算。

全文数据：

权利要求：

百度查询： 中国地质科学院 一种基于变密度变深度约束的重力密度界面反演方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD