基于微扰理论的伽马-伽马密度测井快速正演方法  

A fast forward method of gamma-gamma density logging based on perturbation theory

作　　者：陈前 张锋[1,2] 梁启轩 刁书渊 田立立 范继林 王才志 刘英明[4] 遆永周[5] CHEN Qian;ZHANG Feng;LIANG QiXuan;DIAO ShuYuan;TIAN LiLi;FAN JiLin;WANG CaiZhi;LIU YingMing;TI YongZhou(School of Geosciences,China University of Petroleum(East China),Shandong Qingdao 266580,China;Pilot National Laboratory for Marine Science and Technology,Shandong Qingdao 266580,China;Karamay Hongshan Oilfield Co.,Ltd.,Xinjiang Branch of China National Petroleum Corporation,Xinjiang Karamay 834000,China;Research Institute of Petroleum Exploration and Development,Beijing 100000,China;Institute of Isotope Research,Henan Academy of Sciences,Zhengzhou 450000,China)

机构地区：[1]中国石油大学(华东)地球科学与技术学院,山东青岛266580 [2]青岛海洋科学与技术试点国家实验室,山东青岛266580 [3]中国石油天然气集团总公司新疆分公司克拉玛依红山油田有限责任公司,新疆克拉玛依834000 [4]中国石油勘探开发研究院,北京100000 [5]河南省科学院同位素研究所有限责任公司,郑州450000

出　　处：《地球物理学报》2022年第4期1528-1535,共8页Chinese Journal of Geophysics

基　　金：国家自然科学基金(42174147,41974127);山东省专项经费(2021QNLM020001);国家科技重大专项(2017ZX05019005-005)联合资助.

摘　　要：伽马-伽马密度快速正演结果作为中子和密度联合反演模型的输入,对实时反演确定地层几何结构及属性参数具有重要的意义.基于时间独立的玻耳兹曼输运方程,结合微扰理论与蒙特卡罗(Monte Carlo N Particle Transport Code,MCNP)模拟,形成了伽马-伽马密度快速计算方法.以双源距补偿密度仪器参数为基准,建立了不同基准密度地层条件下的扰动灵敏度函数数据库.分析对比了快速计算和蒙特卡罗模拟在不同岩性密度、水平层状地层及倾斜地层条件下的密度结果与误差.结果表明:不同地层条件下的伽马-伽马密度快速计算结果与MCNP模拟结果吻合,三种岩性(砂岩、碳酸岩和白云岩)密度计算的平均误差小于0.01 g·cm^(-3)且与MCNP模拟速度相比,快速计算效率提高10^(6)倍.The fast forward modeling result of gamma-gamma density is used as the input of the neutron and density joint inversion model,which is of great significance for real-time inversion to determine the formation geometry and attribute parameters.Based on the time-independent Boltzmann transport equation,combined with the perturbation theory and Monte Carlo simulation,a fast calculation method of gamma-gamma density is proposed.The compensation density logging is adopted to establish the database of the perturbation sensitivity function under different reference formation conditions.The density responses and errors of fast calculation and MCNP simulation under different lithological densities,horizontally layered formations and inclined formations are analyzed and compared.The results show that the fast calculation results of gamma-gamma density logging under different formation conditions are consistent with the MCNP simulations.The average error of the three lithological(sandstone,carbonatite and dolomite)density calculations is less than 0.01 g·cm^(-3) and compared with the MCNP simulation speed,the fast calculation efficiency is increased by 10^(6) times.

关 键 词：密度测井 快速正演 微扰理论 蒙特卡罗模拟 

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