杨博，曹书锦,，毛雅静，周勇，朱自强，陈新跃

1.湖南科技大学 页岩气资源利用湖南重点实验室，湖南 湘潭 411201；2.湖南科技大学 资源环境与安全工程学院，湖南 湘潭 411201；3.中南大学 地球科学与信息物理学院，长沙 410083；4.湖南文理学院 资源环境与旅游学院，湖南 常德 415000

0 引言

欧拉反褶积法能使用少量先验信息且不需要准确的解释模型，即可快速准确地圈定异常源的基本轮廓或标识异常源中心，特别适合大尺度重磁位场数据的快速解释。欧拉反褶积法一般通过一定大小的滑动窗口在整个勘探区域移动，根据人工经验或枚举构造指数，在滑动窗口内基于欧拉齐次方程构建线性方程组[1--7]。该线性方程组由人为设定大小的滑动窗口内的观测数据构建，其条件数极大，这导致欧拉反褶积易于受到观测数据噪声、异常源形态及其空间位置和滑动窗口尺寸等因素的影响。

在实际地质情况中，地质构造一般为多个不同形态异常源的组合，传统欧拉反褶积解释工作流程中仅使用单一构造指数易于导致欧拉解发散，而通过枚举构造指数易于导致欧拉解集过于庞大[8--14]，单一过滤标准又将忽视其他异常源[1]。张季生等[15]针对欧拉反褶积方法易受高频干扰的问题，提出低阶的欧拉反褶积与解析信号相结合的位场反演方法，将高频噪声的干扰减少到可以控制的水平[16]。此外，将重力梯度张量解析信号引入欧拉反褶积，或将重力梯度张量数据进行联合欧拉三维反演，以提高欧拉反褶积解的收敛性[17--18]。……