[发明专利]一种基于紧凑约束的重力数据向下延拓方法及系统

说明书

本发明涉及一种基于紧凑约束的重力数据向下延拓方法及系统，包括：获取研究区域各测点的重力数据并进行预处理；确定研究区域的等效源的几何参数；基于预处理后的重力数据结合紧凑约束项和平滑策略得到等效源的密度分布；紧凑约束项根据等效源密度和紧凑因子确定；根据重力等效源的几何参数和等效源的密度得到等效源模型；基于等效源模型进行正演计算得到重力数据的向下延拓结果。通过引入紧凑约束项使得等效源物性集中分布，同时考虑平滑策略压制或消除等效源物性中的发散信号，显著改善了数据形态的畸变，提高了数据幅值的恢复精度，确保大深度向下延拓的稳定性。

技术领域

本发明涉及重力学及重力勘探领域，特别是涉及一种基于紧凑约束的重力数据向下延拓方法及系统。

背景技术

向下延拓是重力数据处理中的重要内容之一，将观测数据从观测面换算至观测面之下的无源空间中，可以提高其对场源的分辨率、增强弱信号、并满足对不同高度重力场信息的需求。例如将航空重力数据下延至地面、将船测重力数据下延至数千米深的海底、或将数百公里高度上的卫星重力数据下延至近地表高度。通过向下延拓获得的高精度与高分辨率重力数据在地学、海洋科学、环境科学、资源勘查和国防安全等领域发挥着重要的作用。

稳定、准确和大深度的向下延拓一直是重力数据处理中的难点和热点问题，其中稳定性是限制向下延拓最大距离和精度的主要因素。向下延拓分为频率域和空间域两类方法，如基于泰勒级数展开的方法、基于傅里叶变换的方法、解析延拓方法、积分-迭代方法、和基于位场中值理论的解算方法等。以上方法或因为求解高阶导数，或因为延拓因子的放大作用，使得高频信号被放大，尤其是噪声(或虚假干扰)信号，会引起数据形态的畸变，甚至“淹没”有效信号，造成下延过程的不稳定。为提高稳定性，通常采取的措施是正则化、滤波器、增加迭代次数、以及多源数据约束。正则化思想的引入，在压制噪声信号方面起到了积极的作用，但也会牺牲信号的分辨率和幅值；滤波器的作用与应用效果同正则化类似。增加迭代次数可以改善高频信号的过度放大，但仍未能根本性地解决含噪声的问题。多源数据间的相互约束，是目前较为流行的一种做法，然而在很多情况下会面临无法获取高质量多源数据的问题。因此，前人的研究工作虽然积极地推动了向下延拓方法的发展，但下延的不稳定问题依然严重制约着延拓距离和精度，需要持续开展深入的研究工作。

等效源方法是重力数据处理中较为常用的技术，该方法利用一组人工构造的虚拟场源(等效源)，代替真实场源在无源空间中产生重力数据；即通过拟合观测数据建立等效源模型，将重力数据处理转换为等效源模型正演计算。该方法具有直接处理原始点位观测数据和顾及数据间同源性的优点，在重力数据处理中表现出了良好的性能；但利用该方法进行大深度向下延拓时，同样会出现数据形态畸变和幅值亏损的问题；其根本原因是等效源物性中存在的“发散”信号。所以，获取高精度及高分辨率的下延结果，对等效源模型提出了更高的要求，需要探索新的技术方法以确保大深度向下延拓的稳定性。

发明内容

本发明的目的是提供一种基于紧凑约束的重力数据向下延拓方法及系统，通过引入紧凑约束项使得等效源物性集中分布，同时考虑平滑策略压制或消除等效源物性中的“发散”信号，显著改善了数据形态畸变的问题，提高了数据幅值的恢复精度，确保大深度向下延拓的稳定性。

为实现上述目的，本发明提供了如下方案：

一种基于紧凑约束的重力数据向下延拓方法，包括：

获取研究区域各测点的重力数据，并对所述重力数据进行预处理；

确定研究区域的等效源的几何参数；

基于预处理后的重力数据结合紧凑约束项和平滑策略得到等效源的密度分布；所述紧凑约束项根据等效源单元密度和紧凑因子确定；

根据所述等效源的几何参数和所述等效源的密度分布得到所述研究区域的等效源模型；

基于所述等效源模型进行正演计算得到所述重力数据的向下延拓结果。