【技术实现步骤摘要】
一种基于牛顿迭代法的密度测井六参数反演计算方法
本专利技术属于石油天然气开发领域，具体涉及一种基于牛顿迭代法的密度测井六参数反演计算方法。
技术介绍
密度测井作为裸眼井测井的重要组成部分，在地层岩性、孔隙度、泥质含量等参数评价方面具有不可或缺的作用。密度测井方法利用多探头的伽马能谱信息，采用精确的多参数反演方法准确评价地层，对于指导储层开发具有重要意义。目前，密度测井方法利用双探测器或三探测器仪器进行地层评价。该方法通过双源距“脊肋图”的方法消除泥饼的影响，评价地层参数，但是由于岩性和密度测量间存在一定的相关性，导致该方法误差较大，且无法评价泥饼及仪器井眼间隙；斯伦贝谢公司曾利用多元回归的方法求解多参数与能窗计数拟合的经验公式，但也存在一定的局限性，并且也无法解决实际测量过程中仪器与地层间存在间隙的问题。因此，亟需提出一种密度测井多参数反演方法，弥补多探测器密度测井仪器在多参数反演方法上存在的不足，解决密度测井多参数反演问题。
技术实现思路
为解决上述技术问题，本专利技术提出了一种基于牛顿...
【技术保护点】
1.一种基于牛顿迭代法的密度测井六参数反演计算方法，采用由一个
【技术特征摘要】
1.一种基于牛顿迭代法的密度测井六参数反演计算方法，采用由一个137Cs伽马源和多个伽马探测器组成的密度测井仪器，其特征在于，具体包括以下步骤：
步骤一：根据密度测井仪器参数及标准刻度井地层参数，利用MCNP模拟软件，建立蒙特卡罗数值模拟MCNP计算模型，模拟密度测井仪器在标准刻度井中的测量过程，得到标准刻度井的模拟能谱，并将刻度井的模拟能谱定义为标准条件能谱；基于密度测井仪器在标准刻度井中测量得到的实测能谱或模拟能谱，确定密度测井仪器中各探测器的能窗范围及伽马计数率；
步骤二：定义地层密度ρb、地层体积光电吸收指数Ub、泥饼密度ρmc、泥饼体积光电吸收指数Umc、泥饼厚度hmc和仪器井眼间隙hs为密度测井六参数，利用MCNP模拟软件，分别改变蒙特卡罗数值模拟MCNP计算模型中的密度测井六参数，得到不同密度测井参数条件下的测井响应及各能窗内的伽马计数率，确定能窗计数率与各密度测井参数的响应关系，建立密度测井六参数响应数据库；
步骤三：密度测井仪器在井中测量得到实测能谱，计算实测能谱与标准条件能谱中对应能窗计数率的差值，结合密度测井六参数响应数据库，利用牛顿迭代法求解非线性超定方程组，计算密度测井六参数的最优解，确定实测能谱对应的地层密度ρb、地层体积光电吸收指数Ub、泥饼密度ρmc、泥饼体积光电吸收指数Umc、泥饼厚度hmc和仪器井眼间隙hs的取值。
2.根据权利要求1所述的一种基于牛顿迭代法的密度测井六参数反演计算方法，其特征在于，所述密度测井仪器中至少设置有三个伽马探测器。
3.根据权利要求1所述的一种基于牛顿迭代法的密度测井六参数反演计算方法，其特征在于，所述步骤一中，蒙特卡罗数值模拟MCNP计算模型中，将地层划分为0.5cm×0.5cm的环形网格栅元。
4.根据权利要求1所述的一种基于牛顿迭代法的密度测井六参数反演计算方法，其特征在于，所述步骤二中，能窗计数率与各密度测井参数的响应关系如下所示：
能窗计数率与地层密度ρb和地层体积光...
【专利技术属性】
技术研发人员：张锋，范继林，陈前，宋灏辰，
申请(专利权)人：中国石油大学华东，
类型：发明
国别省市：山东;37
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。