先验岩石物性约束的地球物理正则化反演  

Geophysical regularization inversion constrained by priori petrophysical properties

作　　者：李曼 张志勇 刘佳成 杨军 鄢毛毛 周小宝 LI Man;ZHANG Zhiyong;LIU Jiacheng;YANG Jun;YAN Maomao;ZHOU Xiaobao(School of Geophysics and Measurement-control Technology,East China University of Technology,Nanchang,Jiangxi 330013,China;Jiangxi Provincial Institute of Cultural Relics and Archaeology,Nanchang,J iangxi 330095,China;Jiangxi Provincial Institute of Geophysical Exploration Technology,Xinyu,J iangxi 338000,China)

机构地区：[1]东华理工大学地球物理与测控技术学院,江西南昌330013 [2]江西省文物考古研究院,江西南昌330013 [3]江西省地球物理勘察技术研究院,江西新余338000

出　　处：《石油地球物理勘探》2024年第6期1420-1432,共13页Oil Geophysical Prospecting

基　　金：国家自然科学基金项目“浅地表可控场源电磁法多参数勘探研究”(42164008);“基于A-φ势三维CSEM高阶自适应有限元正演”(42004061);核资源与环境国家重点实验室联合创新基金项目“地浸砂岩铀矿溶浸液运移井地电法实时监控技术”(2022NRE-LH-08);江西省交通运输厅科技项目“孔底岩溶探测装备与判读技术研究”(2022H0025)联合资助。

摘　　要：随着勘探程度的深入,可用于正则化反演的先验信息越来越丰富和精确。在分析参数变换函数、模糊C均值聚类模型约束、交替方向乘子优化等引入岩石先验物性约束方法的基础上,提出了一种基于反演单元物性互异原则的模型约束方法。基于单元互异原则的模型约束从地质单元数量有限、每个离散反演单元只能属于某一种岩性的本质出发,通过限定离散反演单元的物性取值实现物性引入。将上述先验物性引入与经典正则化反演相结合,构建了统一的目标函数并通过高斯—牛顿法进行优化求解;对四种方案实现岩石物性引入的作用机理和数学本质进行了分析,提出了同时利用多种方案提高反演效果的组合策略;对比了L1和L2范数最小结构模型约束与先验岩石物性引入策略结合的反演效果。模型和实测数据反演结果均表明,充分利用岩石物性测量的先验信息可有效提高反演效果。With the deepening of exploration,the a priori information available for regularization inversion becomes increasingly rich and accurate.A model constraint method is proposed based on the principle of mutual differences in the physical properties of inversion units,after an analysis of methods that introduce a priori rock physical property constraint,including parameter transformation functions,the fuzzy C‑means(FCM)clustering model constraint and the alternating direction optimization of multipliers.The physical properties are introduced via the imposition of constraints on the values assigned to discrete inversion units,given that the essence of the model constraint based on the principle of mutual differences among inversion units is the finite number of geological units and that each discrete inversion unit can only belong to one lithology.A unified objective function is constructed through the combination of the aforementioned introduction of a priori physical properties with classical regularization inversion and is then optimized through the Gauss‑Newton method.The mechanisms and mathematics essence underlying the introduction of rock physical properties through four distinct schemes are examined,and a strategy of employing multiple schemes to improve the inversion effect is proposed.Finally,the inversion effect of combining L1 norm minimum structure model constraints and the introduction strategy of a priori physical property constraints is compared with that of combining L2 norm minimum structure model constraints and the same strategy.The inversion results of the model and field data all demonstrate that leveraging a priori rock physical property information can effectively improve the inversion effect.

关 键 词：先验岩石物性 正则化反演 变换函数 模糊C均值聚类 交替方向乘子算法 

分 类 号：P631[天文地球—地质矿产勘探]