一种计算致密储层自发渗吸两相流体相对渗透率的方法技术

本发明专利技术涉及一种计算致密储层自发渗吸过程两相流体相对渗透率的方法，是一种评价致密储层自发渗吸两相流体渗流动态的有效手段。计算步骤依次包括：(1)岩心准备，用于自发渗吸核磁共振实验；(2)获取样品饱和模拟地层水条件核磁共振T2谱；(3)样品核磁共振自发渗吸实验，获取不同自发渗吸时间T2谱；(4)计算不同自发渗吸时间样品含水饱和度和T2谱对数平均值；(5)计算自发渗吸过程两相流体相对渗透率。本发明专利技术结合饱和水状态和自发渗吸过程核磁共振T2谱，准确计算得到致密储层自发渗吸过程两相流体相对渗透率，为致密油气藏开发动态评价评价奠定基础。

【技术实现步骤摘要】
一种计算致密储层自发渗吸两相流体相对渗透率的方法
本专利技术涉及一种储层自发渗吸两相流体相对渗透率的计算方法，尤其针对致密储层自发渗吸两相流体相对渗透率的计算。主要结合致密储层岩心自发渗吸核磁共振实验，可对致密油对致密储层自发渗吸过程气水(或油水)两相渗流相对渗透率进行定量计算，是一种准确有效的致密储层开发过程自发渗吸动态评价方法。
技术介绍
随着全球能源需求量的剧增和常规油气资源的日益枯竭，以致密油气为代表的非常规资源勘探开发日益受到重视。不同于常规储层，致密油气储层主要发育微纳米孔隙结构，储层物性极差，非均质强，因此，水平井多级压裂技术被用于致密储层开发，致密储层极差的物性使得常规注水驱油开发方式极难实现，而以自发渗吸驱油的开发方式可以取得一定的效果，因此，准确评价致密储层自发渗吸气水(或油水)两相动态，对于指导致密储层高效开发，具有极其重要的现实意义。目前，实验室内测试两相相对渗透率的实验方法主要有稳态法和非稳态法。参照中华人民共和国石油行业标准《SY/T5843-1997(2005)气水相对渗透率测定》和《SY/T5345-1999油水相对渗透率测定》，稳态法和非稳态法对样品物性均有严格要求，如针对气水两相相对渗透率测定，稳态法适用于样品渗透率范围为(0.5～1000)×10-3μm2，而非稳态法适用于样品渗透率大于0.01×10-3μm2。同时，在实验过程中需要准确计量气水、油水两相流体的流量、驱替压力和驱替时间，最终基于达西渗流理论求取对应各相的相渗透率，进而得到各相的相对渗透率。而在自发渗吸过程中，两相渗流过程的作用力主要为毛管压力和重力，并不存在外在驱替压力，因此亦无法准确计量得到，这就使得常规的相对渗透率测试及计算方法并不适用于致密油气藏开发流体自发渗吸过程相对渗透率的计算，从而对自发渗吸过程两相流体的动态表征带来一系列困难。研究表明，基于室内岩心核磁共振实验分析获取的T2谱，可有效指示两相流体渗流过程不同含水饱和度条件流体分布，从而可以建立起相对渗透率与含水饱和度的函数关系。因此，通过监测致密储层自发渗吸过程核磁共振T2谱分布，结合核磁共振渗透率计算SDR模型可推导得出不同含水饱和度条件流体相渗透率，进而获取不同含水饱和度条件各相流体的相对渗透率，为致密储层开发过程自发渗吸两相流体渗流动态评价提供支撑。
技术实现思路
本专利技术目的在于提供一种计算致密储层自发渗吸两相流体相对渗透率的方法，由此得到的计算模型、认识和结论，丰富了致密油气藏开发过程储层自发渗吸流体渗流动态评价方法，克服了传统实验技术对致密油气藏开发过程自发渗吸动态评价的局限。为达到以上技术目的，本专利技术提供以下技术方案。致密储自发渗吸两相流体相对渗透率计算公式如下：据文献(KenyonWE.Nuclearmagneticresonanceasapetrophysicalmeasurement.NuclearGeophysics,1992,6,153-171.)(CoatesGR,XiaoL,PrammerMG.NMRloggingprinciplesandapplications.HalliburtonEnergyServices:Houston,TX,1999.)得到一种基于核磁共振T2谱计算储层渗透率的方法，即SDR模型。可由下式进行表示：式中，k—核磁计算渗透率，单位为毫达西(mD)；φ—孔隙度，单位为小数(f)；C—与岩石矿物成分相关的系数，可通过实验数据拟合得，无量纲；m，n—与岩石物性和微观结构相关的系数，分别等于4和2，无量纲；T2lm—横向弛豫时间T2对数平均值，单位为毫秒(ms)。对于饱和水岩心样品，T2谱对数平均值可以通过以下公式计算：式中，T2i—核磁共振T2谱对应的横向弛豫时间，单位为毫秒(ms)；m(T2i)—不同横向弛豫时间T2i对应的幅度值，单位为幅度(a.u.)。从公式(1)可知，对于两相流体渗流(油水两相，油为非润湿相，水为润湿相；气水两相，气为非润湿相，水为润湿相；核磁共振T2谱监测中，水有核磁共振信号，模拟油(氟油)气均无核磁共振信号)，则有核磁共振信号润湿相流体(w)的相渗透率可以表示为，同理可得，对于两相流体渗流状态，无核磁共振信号非润湿相流体(nw)的相渗透率可以表示为，式中，kw—两相流体共渗时润湿相流体(水)相渗透率，单位为毫达西(mD)；knw—两相流体共渗时非润湿相流体(气或油)相渗透率，单位为毫达西(mD)；Sw—两相流体共渗时润湿相流体(水)饱和度，单位为小数(f)；Snw—两相流体共渗时非润湿相流体(气或油)饱和度，Snw＝1-Sw，单位为小数(f)；T2lmw—两相流体渗流润湿相流体横向弛豫时间T2谱对数平均值，单位为毫秒(ms)；T2lmnw—两相流体渗流非润湿相流体横向弛豫时间T2谱对数平均值，单位为毫秒(ms)。自发渗吸一段时间后，样品润湿相流体饱和度(Swt)可由下式表示：对于部分饱和水岩心样品，假定其含水饱和度为Swt，则核磁共振T2谱对数平均值可以通过以下公式计算：对于无核磁共振信号相(油或气)流体，对应的核磁共振T2谱对数平均值可以通过以下公式计算：式中，Sw—样品含水饱和度，单位为百分之(％)；Swt—自发渗吸一段时间后样品含水饱和度，单位为百分之(％)；m(T2i)Sw＝Swt—自发渗吸一段时间后，样品含水饱和度为Swt时核磁共振信号幅度值，单位为幅度(a.u.)；T2lmw—自发渗吸一段时间后，样品含水饱和度为Swt时有核磁共振信号相(水)T2谱对数平均值，单位为毫秒(ms)；T2lmnw—自发渗吸一段时间后，样品含水饱和度为Swt时无核磁共振信号相(气或油)T2谱对数平均值，单位为毫秒(ms)；t2—横向弛豫时间，单位为毫秒(ms)；T2max—横向弛豫时间最大值，单位为毫秒(ms)。则联立公式(1)和(3)可知，致密储层自发渗吸两相流体渗流润湿相流体(水)相对渗透率表示为：则联立公式(1)和(4)可知，致密储层自发渗吸两相流体渗流非润湿相流体(气或油)相对渗透率表示为：式中，krw—致密储层自发渗吸过程润湿相流体相对渗透率，无量纲；krnw—致密储层自发渗吸过程非润湿相流体相对渗透率，无量纲；kw—致密储层自发渗吸过程润湿相流体相渗透率，mD；knw—致密储层自发渗吸过程非润湿相流体相渗透率，mD。从公式(8)和(9)可以看出，要确定致密储层自发渗吸两相流体的相对渗透率，首先必须准确求取自发渗吸过程储层流体的饱和度(Sw，Snw)和不同饱和度条件横向弛豫时间T2谱对数平均值(T2lmw，T2lmnw)。一种计算致密储层自发渗吸两相流体相对渗透率的方法，依次包括以下步骤：(1)实验岩心准备；(2)样品饱和模拟地层水及对应核磁共振T2谱测试；(3)自发渗吸核磁共振实验；(4)自发渗吸过程T2谱变化分析，确定对应流体饱和度和核磁共振T2谱对数平均值；(5)计算致密储层岩心自发渗吸过程两相流体相对渗透率。以下对各步骤进行具体说明：(1)实验岩心准备按照中华人民共和国国家标准GB/T29172-2012《岩心分析方法》，对岩样进行选取、制备、清洗。岩样两端断面切割平整，用于自发渗吸核磁共振实验。用高精度电子天平称量样品干重m0，游标卡尺测量长度L本文档来自技高网...

【技术保护点】
1.一种计算致密储层自发渗吸两相流体相对渗透率的方法，依次包括以下步骤：(1)岩心准备，样品两端断面磨平用于自发渗吸核磁共振实验；(2)样品饱和模拟地层水及对应核磁共振T2谱测试，采用抽真空加压饱和48h、再以0.01ml/min驱替模拟地层水24h的方法饱和致密储层岩心样品，后采用低场核磁共振实验装置检测饱和水状态样品核磁共振T2谱；(3)将岩心样品洗盐烘干，样品用无核磁共振信号的热缩管包裹，采用端面自发渗吸的方式，在实验室环境下恒温箱内进行自发渗吸实验，检测自发渗吸过程不同时间核磁共振T2谱，直至检测核磁共振T2谱趋于重叠；(4)自发渗吸过程核磁共振T2谱变化分析，确定自发渗吸过程样品含水饱和度和两相流体核磁共振T2谱对数平均值；(5)计算得到致密储层岩心自发渗吸过程两相流体相对渗透率。

【技术特征摘要】
1.一种计算致密储层自发渗吸两相流体相对渗透率的方法，依次包括以下步骤：(1)岩心准备，样品两端断面磨平用于自发渗吸核磁共振实验；(2)样品饱和模拟地层水及对应核磁共振T2谱测试，采用抽真空加压饱和48h、再以0.01ml/min驱替模拟地层水24h的方法饱和致密储层岩心样品，后采用低场核磁共振实验装置检测饱和水状态样品核磁共振T2谱；(3)将岩心样品洗盐烘干，样品用无核磁共振信号的热缩管包裹，采用端面自发渗吸的方式，在实验室环境下恒温箱内进行自发渗吸实验，检测自发渗吸过程不同时间核磁共振T2谱，直至检测核磁共振T2谱趋于重叠；(4)自发渗吸过程核磁共振T2谱变化分析，确定自发渗吸过程样品含水饱和度和两相流体核磁共振T2谱对数平均值；(5)计算得到致密储层岩心自发渗吸过程两相流体相对渗透率。2.如权利要求1所述一种计算致密储层自发渗吸两相流体相对渗透率的方法，其特征在于，步骤(4)采用如下公式计算自发渗吸过程样品含水饱和度：式中，Swt—自发渗吸一段时间后样品含水饱和度，％；T2max—横向弛豫时间最大值，ms；m(T2i)Sw＝Swt—自发渗吸一段时间后，样品含水饱和度为Swt时核磁共振信号幅度值，a.u.；m(T2i)Sw＝100％—样品...

【专利技术属性】
技术研发人员：陈猛，戴家才，刘向君，况晏，秦民君，裴阳，王中涛，苟顺超，
申请(专利权)人：西南石油大学，
类型：发明
国别省市：四川,51