一种考虑微观孔喉分布的孔渗模型建立方法及可读存储介质技术

本发明专利技术公开考虑微观孔喉分布的孔渗模型建立方法，包括统计第一岩心样品的高压压汞或核磁共振实验结果，绘制平均孔喉半径与常压气测渗透率的第一关系图版；统计第二岩心样品的覆压孔渗实验结果，绘制常压气测渗透率与覆压渗透率的第二关系图版及覆压渗透率与覆压孔隙度的第三关系图版；利用第一关系图版和第二关系图版，将靶区的孔喉半径分布关系转为覆压渗透率分布关系；基于靶区地质模型，建立服从0～1上平均分布的中间随机属性模型，利用覆压渗透率分布关系获得覆压渗透率累积分布，建立渗透率属性模型；利用第三关系图版，建立孔隙度属性模型；基于渗透率和孔隙度属性模型，构建数学模型进行数值模拟研究，量化不同类型油藏开发效果的差异。藏开发效果的差异。藏开发效果的差异。

【技术实现步骤摘要】
一种考虑微观孔喉分布的孔渗模型建立方法及可读存储介质


[0001]本专利技术涉及油田地质建模
，更具体地，涉及一种考虑微观孔喉分布的孔渗模型建立方法及可读存储介质。

技术介绍

[0002]地质模型建立技术是海上油田开发方案编制中的关键一环，是后续油藏数值模拟研究和合理开发指标预测的基础。不同类型储层微观孔喉差异明显，即使储层物性(孔隙度、渗透率、平均孔喉半径等)相同或相近，其微观孔喉分布各异，导致流体在储层中的渗流能力和油水/油气两相渗流特征存在较大的差异，严重影响了油藏合理开发方案编制和开发指标预测。此外，海上油田开发初期探井较少，以断块油藏为例每个区块普遍仅有1口探井，对应的常规孔渗模型建立方法包括基于探井资料的均质建模、结合周边邻近区块探井资料进行插值建模等，这些方法无法有效刻画不同类型储层孔喉分布差异，导致建立的孔渗模型难以有效考虑储层物性分布差异的影响。中国专利申请，公开了一种确定储层渗透率的方法，所述方法按照以下步骤进行操作：步骤(1)对岩心进行核磁T2实验及压汞实验，并获取核磁T2谱和压汞实验数据；步骤(2)利用压汞实验数据建立孔喉划分界限；步骤(3)根据孔喉划分标准建立地区核磁T2谱孔隙分量划分界限；步骤(4)根据步骤(3)所划分的界限，计算不同孔喉对应的弛豫时间划分区间的包络面积与总面积的比值，即孔隙分量值；步骤(5)根据步骤(4)所得孔隙分量值计算储层的渗透率。但是该技术方案仍无法有效刻画不同类型储层孔喉分布差异，导致建立的孔渗模型难以有效考虑储层物性分布差异的影响。

技术实现思路
<br/>[0003]本专利技术为克服上述现有技术中孔渗模型建立方法无法有效刻画不同类型储层孔喉分布差异，导致建立的孔渗模型难以有效考虑储层物性分布差异的问题，提供一种考虑微观孔喉分布的孔渗模型建立方法。
[0004]为解决上述技术问题，本专利技术采用的技术方案是：一种考虑微观孔喉分布的孔渗模型建立方法，包括以下步骤：
[0005]S1：统计研究区域大量第一岩心样品的高压压汞或核磁共振实验结果，根据该高压压汞或核磁共振实验结果绘制平均孔喉半径与常压气测渗透率的第一关系图版；统计与第一岩心样品同一区域内相同井段的第二岩心样品的覆压孔渗实验结果，根据该覆压孔渗实验结果绘制常压气测渗透率与覆压渗透率的第二关系图版以及覆压渗透率与覆压孔隙度的第三关系图版；
[0006]S2：利用步骤S1的第一关系图版和第二关系图版，将目标靶区获得的孔喉半径分布关系转化为覆压渗透率分布关系；
[0007]S3：基于目标靶区地质模型，建立服从0～1上平均分布的中间随机属性模型，并利用步骤S2中覆压渗透率分布关系计算获得覆压渗透率累积分布，将中间随机属性模型中每个网格数值与覆压渗透率累积分布结果匹配，建立渗透率属性模型；同时，利用步骤S1中覆
压渗透率与覆压孔隙度的第三关系图版，建立孔隙度属性模型；
[0008]S4：基于步骤S3中的渗透率属性模型和孔隙度属性模型，构建相应的数学模型进行数值模拟研究，量化不同类型油藏开发效果的差异。
[0009]优选地，在所述步骤S1中，根据第一关系图版拟合获得平均孔喉半径与常压气测渗透率关系式为式中，为平均孔喉半径，k
g
为常压气测渗透率，a、b为拟合系数。
[0010]优选地，在所述步骤S1中，根据第二关系图版拟合获得常压气测渗透率与覆压渗透率关系式为式中，k
op
为覆压渗透率，c、d为拟合系数。
[0011]优选地，在所述步骤S1中，根据第三关系图版拟合获得覆压渗透率与覆压孔隙度关系式为式中，φ
op
为覆压孔隙度,e、f为拟合系数。
[0012]优选地，在所述步骤S2中，根据S1中的平均孔喉半径与常压气测渗透率关系式和常压气测渗透率与覆压渗透率关系式，将目标靶区得到的孔喉半径分布关系转化为覆压渗透率分布关系。
[0013]优选地，所述步骤S2具体为：利用平均孔喉半径与常压气测渗透率关系式将孔喉半径分布关系转化为常压气测渗透率分布关系，然后利用常压气测渗透率与覆压渗透率关系式将常压气测渗透率分布关系转化为覆压渗透率分布关系。
[0014]优选地，在所述步骤S3中，将中间随机属性模型中每个网格数值与覆压渗透率累积分布结果进行匹配，如某一网格中数值不超过某一覆压渗透率对应的累积概率值，将该覆压渗透率数值赋予该网格，利用该方法对中间随机属性模型中所有网格的数值进行逐一匹配，获得模型全部网格的覆压渗透率值，建立渗透率属性模型。
[0015]优选地，在所述步骤S3中，利用步骤S1中的覆压渗透率与覆压孔隙度关系式计算获得模型所有网格的覆压孔隙度值，建立孔隙度属性模型。
[0016]优选地，在所述步骤S2中，利用高压压汞或核磁共振实验获得的孔喉半径分布关系。
[0017]本专利技术另一方面提供一种存储有程序指令的可读存储介质，当所述程序指令被计算设备读取并执行时，使得所述计算设备执行如上所述的考虑微观孔喉分布的孔渗模型建立方法。
[0018]与现有技术相比，本专利技术的有益效果是：本专利技术所针对常规孔渗模型建立方法无法有效刻画不同类型储层微观孔喉分布差异的问题，基于区域大量高压压汞、核磁共振和覆压孔渗实验，建立孔喉半径分布关系转化覆压孔渗分布关系的方法，然后将目标靶区高压压汞或核磁共振实验获得的储层微观孔喉分布结果引入到孔渗模型中，建立考虑储层微观孔喉分布的地质模型。在此基础上，构建数学模型进行数值模拟研究，量化不同类型储层微观孔喉分布差异对油藏开发效果的影响。该方法拓展了海上油田早期开发的地质建模方法，解决了现有地质建模方法难以考虑储层微观孔喉分布差异对开发方案与指标影响的难题，为不同类型油藏数值模拟研究和开发方案编制提供了一种表征合理、可靠性强、操作性强的孔渗模型建立方法。
附图说明
[0019]图1是本专利技术孔渗模型建立方法的流程图；
[0020]图2是本专利技术实施例2中平均孔喉半径与常压气测渗透率的第一关系图版；
[0021]图3是本专利技术实施例2中常压气测渗透率与覆压渗透率的第二关系图版；
[0022]图4是本专利技术实施例2中覆压渗透率与覆压孔隙度的第三关系图版；
[0023]图5是本专利技术实施例2中细砂岩的平均孔喉半径分布图；
[0024]图6是本专利技术实施例2中砂砾岩的平均孔喉半径分布图；
[0025]图7是本专利技术实施例2中细砂岩覆压孔隙度分布图；
[0026]图8是本专利技术实施例2中细砂岩覆压渗透率分布图；
[0027]图9是本专利技术实施例2中砂砾岩覆压孔隙度分布图；
[0028]图10是本专利技术实施例2中砂砾岩覆压渗透率分布图；
[0029]图11是本专利技术实施例2中细砂岩渗透率属性模型图；
[0030]图12是本专利技术实施例2中砂砾岩渗透率属性模型图；
[0031]图13是本专利技术实施例2中细砂岩孔隙度属性模型图；
[0032]图14是本专利技术实施例2中砂砾岩孔隙度属性模型图本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种考虑微观孔喉分布的孔渗模型建立方法，其特征在于，包括以下步骤：S1：统计研究区域大量第一岩心样品的高压压汞或核磁共振实验结果，根据该高压压汞或核磁共振实验结果绘制平均孔喉半径与常压气测渗透率的第一关系图版；统计与第一岩心样品同一区域内相同井段的第二岩心样品的覆压孔渗实验结果，根据该覆压孔渗实验结果绘制常压气测渗透率与覆压渗透率的第二关系图版以及覆压渗透率与覆压孔隙度的第三关系图版；S2：利用步骤S1的第一关系图版和第二关系图版，将目标靶区获得的孔喉半径分布关系转化为覆压渗透率分布关系；S3：基于目标靶区地质模型，建立服从0～1上平均分布的中间随机属性模型，并利用步骤S2中覆压渗透率分布关系计算获得覆压渗透率累积分布，将中间随机属性模型中每个网格数值与覆压渗透率累积分布结果匹配，建立渗透率属性模型；同时，利用步骤S1中覆压渗透率与覆压孔隙度的第三关系图版，建立孔隙度属性模型；S4：基于步骤S3中的渗透率属性模型和孔隙度属性模型，构建相应的数学模型进行数值模拟研究，量化不同类型油藏开发效果的差异。2.根据权利要求1所述的考虑微观孔喉分布的孔渗模型建立方法，其特征在于，在所述步骤S1中，根据第一关系图版拟合获得平均孔喉半径与常压气测渗透率关系式为式中，为平均孔喉半径，k
g
为常压气测渗透率，a、b为拟合系数。3.根据权利要求2所述的考虑微观孔喉分布的孔渗模型建立方法，其特征在于，在所述步骤S1中，根据第二关系图版拟合获得常压气测渗透率与覆压渗透率关系式为式中，k
op
为覆压渗透率，c、d为拟合系数。4.根据权利要求3所述的考虑微观孔喉分布的孔渗模型建立方法，其特征在于，在所述步骤S1中，根据第三关系图版拟合获得覆压渗透率...

【专利技术属性】
技术研发人员：王磊，唐宋军，雷霄，王雯娟，洪楚侨，隋波，岳绍飞，彭小东，
申请(专利权)人：中海石油中国有限公司湛江分公司，
类型：发明
国别省市：