CN113916747B - 一种应用于黏弹性流体测试启动压力梯度的修订方法 - Google Patents
一种应用于黏弹性流体测试启动压力梯度的修订方法 Download PDFInfo
- Publication number
- CN113916747B CN113916747B CN202111195761.3A CN202111195761A CN113916747B CN 113916747 B CN113916747 B CN 113916747B CN 202111195761 A CN202111195761 A CN 202111195761A CN 113916747 B CN113916747 B CN 113916747B
- Authority
- CN
- China
- Prior art keywords
- viscoelastic fluid
- pressure drop
- inner diameter
- flow
- flow channel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000012530 fluid Substances 0.000 title claims abstract description 117
- 238000012360 testing method Methods 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000013178 mathematical model Methods 0.000 claims abstract description 21
- 238000006073 displacement reaction Methods 0.000 claims abstract description 14
- 230000001052 transient effect Effects 0.000 claims abstract description 6
- 239000011148 porous material Substances 0.000 claims description 21
- 239000000243 solution Substances 0.000 claims description 11
- 230000008602 contraction Effects 0.000 claims description 10
- 230000000977 initiatory effect Effects 0.000 claims description 10
- 238000010008 shearing Methods 0.000 claims description 6
- 239000008154 viscoelastic solution Substances 0.000 claims description 5
- 238000004088 simulation Methods 0.000 claims description 4
- 238000002474 experimental method Methods 0.000 abstract description 6
- 230000035699 permeability Effects 0.000 description 13
- 230000035882 stress Effects 0.000 description 11
- 239000007788 liquid Substances 0.000 description 9
- 239000003921 oil Substances 0.000 description 8
- 239000011435 rock Substances 0.000 description 7
- 239000004973 liquid crystal related substance Substances 0.000 description 6
- 229920000642 polymer Polymers 0.000 description 6
- 238000002347 injection Methods 0.000 description 5
- 239000007924 injection Substances 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 239000010779 crude oil Substances 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 230000000704 physical effect Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 230000010339 dilation Effects 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 230000032683 aging Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 239000003795 chemical substances by application Substances 0.000 description 2
- 239000000295 fuel oil Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 238000001179 sorption measurement Methods 0.000 description 2
- 206010043528 Throat tightness Diseases 0.000 description 1
- 239000000853 adhesive Substances 0.000 description 1
- 230000001070 adhesive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000036571 hydration Effects 0.000 description 1
- 238000006703 hydration reaction Methods 0.000 description 1
- 230000002209 hydrophobic effect Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000005325 percolation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000000518 rheometry Methods 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000002904 solvent Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 238000012108 two-stage analysis Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/082—Investigating permeability by forcing a fluid through a sample
- G01N15/0826—Investigating permeability by forcing a fluid through a sample and measuring fluid flow rate, i.e. permeation rate or pressure change
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
Landscapes
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Dispersion Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Fluid Mechanics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种应用于黏弹性流体测试启动压力梯度的修订方法,具体步骤如下:S1、配制目标溶液,测试黏弹性流体溶液的流变曲线并进行流变规律拟合;S2、构建黏弹性流体通过多孔介质时所产生流动阻力的数学模型;步骤3、用瞬态动用法进行黏弹性流体驱替稠油的启动压力梯度测试,获取启动压力实验数据;应用数学模型,分析计算黏弹性流体注入多孔介质中产生的附加阻力;S4、将启动压力实验数据减去附加阻力,然后获得的压差值除以岩心长度,即获得启动压力梯度。这样的修订方法为黏弹性流体作为驱替剂开展启动压力梯度测试实验提供了基础。
Description
技术领域
本发明涉及一种应用于黏弹性流体测试启动压力梯度的修订方法,属于石油技术领域。
背景技术
在根据经典Darcy定律,渗流速度与压力梯度之间表现为一条过原点的线性关系,只要偏离这种线性关系的渗流过程,都称为非Darcy渗流,又称为非线性渗流。“启动压力梯度”就是这种非线性渗流中的一个关键名词:1)对于低渗透油藏,流体渗流时必须有一个附加的压力梯度克服岩石表面吸附膜或水化膜引起的阻力才能流动;2)对于稠油油藏,由于黏度高,渗流阻力大,液固界面及液液界面的相互作用力大,导致只有当驱动压力梯度超过某一初始压力梯度时,稠油才开始流动(王晓琴,吴聚,冉艳,等.非线性渗流对异常高压气藏产能的影响[J].岩性油气藏,2012,24(4):125-128)。目前确定启动压力梯度的方法主要有室内物理实验模拟、数值实验和试井解释3种方法,并且目前国内外没有统一的测定标准和方法(陈元千.线性流的启动压力梯度不能用于平面径向流方程[J].石油学报,2011,32(6):1088-1091;贝尔J.多孔介质流体动力学[M].李竞生,陈崇希,译.北京:中国建筑出版社,1983:95-97.)。
对于真实的启动压力梯度的测试方法中,“瞬间动用法”被广大学者所应用。瞬间动用法是当孔隙中充满流体时,在入口端加压情况下,驱替压差从低逐渐升高,岩心中压力梯度越来越大,当压力梯度超过某一值后,注入的流体在克服注入端面各种阻力后进入岩心孔道,如果忽略注入过程中流体发生微小压缩的体积,那么由于压力的传递作用,孔隙中的流体就会发生微流动并传递到出口端面,而在出口端面会有流体渗出,实验室就是通过测定该流体被动用瞬间的压力作为最小启动压力(李南星,刘林玉,郑锐,等.鄂尔多斯盆地镇泾地区超低渗透储层评价[J].岩性油气藏,2011,23(2):41-45;张代燕,王子强,王殿生,等.低渗透油藏最小启动压力梯度实验研究[J].新疆地质,2011,29(1):106-109;谢晓庆,张贤松,张凤久,等.薄层低品位油藏孔隙结构及渗流特征[J].成都理工大学学报(自然科学版),2013,40(1):34-39)。
因此,对于启动压力梯度测试过程中,一般采用了注气或者注气作为驱替介质,这两者几乎可以看成是理想流体;而且不同流体作为驱替剂时,测试出来的启动压力梯度存在差异(朱维耀,田巍,朱华银,等.一种测试油气藏岩石启动压力梯度的优化方法[J].石油与天然气地质,2017,38(1):175-180)。然后如果采用黏弹性流体作为驱替介质,再应用“瞬态动用法”测试,则存在着明显的压力误差,原因就是来自于高黏物性和弹性作用特征使其在驱动多孔介质时,依然产生了很高的注入压力值,从而使测试出来的启动压力梯度不够准确。
现目前的所有文献中,并没有阐述用黏弹性流体作为驱动介质进行实验的研究报道,应该就是受到黏弹性流体自身固有属性所限制。因此,如果采用黏弹性流体进行实验,有必要针对其黏弹性特征产生的影响进行修订。本发明就是针对上述问题,设计的一种应用于黏弹性流体测试启动压力梯度的修订方法。
发明内容
为了克服现有技术中的缺点,本发明提供一种应用于黏弹性流体测试启动压力梯度的修订方法,本发明有效修订了黏弹性流体流动过程中由于黏弹性作用产生较大的阻力,影响测试启动压力梯度的结果,实现了采用黏弹性流体作为驱替剂进行启动压力梯度的实验测试,为研究多孔介质中非线性渗流的基础研究提供了支持。
本发明解决上述技术问题所提供的技术方案是:一种应用于黏弹性流体测试启动压力梯度的修订方法,具体包括以下步骤:
S1、提供黏弹性流体,测得所述黏弹性流体溶液的流变数据;根据Carreau-Yasuda模型和所述流变数据,获取所述黏弹性流体的目标参数;
S2、构建黏弹性流体通过多孔介质时所产生流动阻力的数学模型;
S3、提供所述黏弹性流体的启动压力实验数据,根据所述启动压力实验数据、所述目标参数、所述数学模型,得到附加阻力;
S4、根据所述附加压力和启动压力实验数据获得所述黏弹性流体的有效启动的压力梯度。
可选的,在一些实施例中,在S3中,提供所述黏弹性流体的启动压力实验数据,根据所述启动压力实验数据、所述目标参数,得到附加阻力包括:
开展黏弹性流体的启动压力实验,获得整体启动压力梯度;
分析目标多孔介质的物性条件,得到多孔介质模型参数,根据多孔介质模型参数和目标参数,计算黏弹性流体通过多孔介质时的附加阻力。
可选的,在一些实施例中,所述S1中,所述黏弹性流体的流变数据通过流变仪测试得到;
所述流变数据为剪切区间范围不小于0.01-10000s-1的测试数据。
可选的,在一些实施例中,所述S1中,所述Carreau-Yasuda模型为:
其中,μ为不同剪切速率下的表观黏度,mPa·s;μ0为零剪切黏度,mPa·s;μinf为极限剪切黏度,mPa·s;λ为松弛时间,n为幂率指数,a为常数,x为剪切速率。
可选的,在一些实施例中,所述S2中,所述数学模型为:
其中,△P为黏弹性溶液通过m个收缩-扩张流道组成孔喉结构时的总压降;△Pc为黏性压降;△Pd为弹性压降;q为流体通过微元段的流量;m为收缩-扩张流道组成个数;μe为弹性黏度;μv为黏性黏度;α为流道的壁面与轴线的角度;n为幂率指数;r2为收缩道半径;r1为扩张道半径。
可选的,在一些实施例中,所述S2中,所述数学模型通过以下步骤得到:
S21:建立包括m个收缩流道和m个扩张流道的流道模型,进行黏弹性流体在所述流道模型中的流通模拟;其中m为大于或等于1的自然数;
S22:假设仅存在剪切流动,建立所述黏弹性流体在所述收缩流道中的黏性压降和扩张流道中的剪切压降;
假设仅存在拉伸流动,建立所述黏弹性流体在所述收缩流道中的弹性压降和在所述扩张流道中的弹性压降;
S23:根据所述黏弹性流体在所述收缩流道中的黏性压降、在所述扩张流道中的剪切压降、在所述收缩流道中的弹性压降和在所述扩张流道中的弹性压降,得到所述数学模型。
可选的,在一些实施例中,所述黏弹性流体在所述收缩流道中的黏性压降为:
所述黏弹性流体在所述收缩流道中的剪切压降为:
所述黏弹性流体在所述扩张流道中的弹性压降为:
其中,β为收缩流道中的外径与内径之间的夹角;
所述黏弹性流体在所述扩张流道中的弹性压降为:
可选的,在一些实施例中,所述收缩流道为内径逐渐变小的变径管道;
所述扩张流道为内径逐渐变大的变径管道;
所述收缩流道的内径较小一端与所述扩张流道的内径较小一端连接,且所述收缩流道的内径较小的一端的内径等于所述扩张流道的内径较小一端的内径、所述收缩流道的内径较大的一端的内径等于所述扩张流道的内径较大一段的内径。
可选的,在一些实施例中,所述S3中,所述启动压力实验数据用岩心样品通过瞬态动用法进行黏弹性流体驱替稠油的启动压力梯度测试得到。
可选的,在一些实施例中,所述S3中,所述启动压力实验数据包括实验结果数据和实验参数;
所述实验参数包括岩心样品的长度。
可选的,在一些实施例中,所述S4中的多孔介质物性条件,根据目标油藏测试的渗透率和微观孔喉结构分析,明确渗透率与孔喉结构的对应关系,在气测渗透率的基础上,通过渗透率明确该条件下的多孔介质孔喉结构组成,即孔喉长度、宽度、夹角、数目范围等。
可选的,在一些实施例中,所述S4中,所述黏弹性流体的有效启动的压力梯度的值等于所述实验结果数据与附加阻力的差值与所述岩心样品的长度的比值。
本发明具有以下有益效果:本发明解决难以测定现有黏弹性流体动态吸附量的客观问题,为更为准确构建数学模型奠定基础,为油田聚合物驱数值化技术提供数据和方法支撑。
附图说明
图1为本发明的孔喉扩张-收缩通道示意图;
图2为本发明的孔喉收缩-扩张通道示意图;
图3为本发明实施例中的聚合物流变曲线图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
此外,下面所描述的本发明不同实施方式中所涉及的技术特征只要彼此之间未构成冲突就可以相互结合。
本发明提供一种应用于黏弹性流体测试启动压力梯度的修订方法,具体包括以下步骤:
S1、提供黏弹性流体,测得所述黏弹性流体溶液的流变数据;根据Carreau-Yasuda模型和所述流变数据,获取所述黏弹性流体的目标参数;
S2、构建黏弹性流体通过多孔介质时所产生流动阻力的数学模型;
S3、开展黏弹性流体的启动压力实验,获得整体启动压力梯度;
分析目标多孔介质的物性条件,得到多孔介质模型参数,根据多孔介质模型参数和目标参数,计算黏弹性流体通过多孔介质时的附加阻力;
S4、根据所述附加压力和启动压力实验数据获得所述黏弹性流体的有效启动的压力梯度。
需要说明的是,所述S1中,所述黏弹性流体的流变数据通过流变仪测试得到。根据所述测试数据获得流变曲线,所述流变数据为剪切区间范围不小于0.01-10000s-1的测试数据;使得到的流变曲线较为完整。将得到的流变曲线,得到目标参数。
需要说明的是,所述目标参数包括:Carreau-Yasuda中的所有未知量,包括零剪切黏度μ0,极限剪切黏度μinf,常数a,幂率指数n,松弛时间λ。其中,所述零剪切黏度和极限剪切黏度通过流变曲线得到,即流变曲线中剪切速率的最小值为零剪切黏度μ0,流变曲线中剪切速率的最大值为极限剪切黏度μinf。将流变数据带入Carreau-Yasuda进行非线性回归拟合,得到拟合曲线常数a、幂率指数n,松弛时间λ。
需要说明的是,所述S1中,所述Carreau-Yasuda模型为:
在一些实施例中,所述S2中,所述数学模型为:
其中,△P为黏弹性溶液通过m个收缩-扩张流道组成孔喉结构时的总压降;△Pc为黏性压降;△Pd为弹性压降;q为流体通过微元段的流量;m为收缩-扩张流道组成个数;μe为弹性黏度;μv为黏性黏度;α为流道的壁面与轴线的角度;n为幂率指数;r2为收缩道半径;r1为扩张道半径。
在一些实施例中,所述S2中,所述数学模型通过以下步骤得到:
S21:建立包括m个收缩流道和m个扩张流道的流道模型,进行黏弹性流体在所述流道模型中的流通模拟;其中m为大于或等于1的自然数;
S22:假设仅存在剪切流动,建立所述黏弹性流体在所述收缩流道中的黏性压降和扩张流道中的剪切压降;
假设仅存在拉伸流动,建立所述黏弹性流体在所述收缩流道中的弹性压降和在所述扩张流道中的弹性压降;
S23:根据所述黏弹性流体在所述收缩流道中的黏性压降、在所述扩张流道中的剪切压降、在所述收缩流道中的弹性压降和在所述扩张流道中的弹性压降,得到所述数学模型。
在一些实施例中,所述黏弹性流体在所述收缩流道中的黏性压降为:
所述黏弹性流体在所述收缩流道中的剪切压降为:
所述黏弹性流体在所述扩张流道中的弹性压降为:
其中,β为扩张流道中的外径与内径之间的夹角;
所述黏弹性流体在所述扩张流道中的弹性压降为:
需要说明的是,所述流道模型为孔喉收缩流道特征模型,其为采用CFD建立的物理模型。其中,收缩流长度为了l1,收缩流道壁面与轴线的角度为a,扩张流道长度l2,扩张流道壁面与轴线的角度为β,左端截面半径为r1,孔喉处截面半径为r2,右端截面半径为r3。
所述数学模型的建立基于以下假设:黏弹性流体在收缩流道中的压降等于剪切压降和弹性压降的和值。
需要说明的是,所述数学模型分收缩流道和扩张流道两段分析,在分析过程中,在收缩流道和扩张流道沿轴线方向分别任取一截长度为dx、半径为r的流体微元段,其具体如图1和图2所示。
所述收缩流道为内径逐渐变小的变径管道;
所述扩张流道为内径逐渐变大的变径管道;
所述收缩流道的内径较小一端与所述扩张流道的内径较小一端连接,且所述收缩流道的内径较小的一端的内径等于所述扩张流道的内径较小一端的内径、所述收缩流道的内径较大的一端的内径等于所述扩张流道的内径较大一段的内径。
需要说明的是,式(1)通过以下步骤推导获得:
A1:假设仅存在剪切流动,不考虑弹性作用,则对于收缩流道作用在流体微元段上的黏性阻力和压降的平衡方程式为:
πr2dPcs=-2πrdxsecα(τccosα) (5)
其中,τc为收缩流道中的剪切应力,dPcs为收缩流道微元段两端的剪切压降;x为流体运动长度;sec为正割函数,secα=1/cosα。
A2:根据聚合物溶液的剪切应力与剪切速率符合是幂律关系的,那么流体微元段的剪切速率可近似用圆管中的剪切速率,则剪切应力为:
其中,K为稠度系数;n为幂率指数;为剪切速率;μv为剪切黏度;q为流体通过流体微元段的流量。
A3:建立dPcs的计算式,其计算式为:
其中,dPcs通过以下方法建立:
建立式(8)、式(9):
组合将式(6)、式(8)、式(9)带入式(5)即可得到式(7)。
A4:在区间[r1,r2]上对式(7)两端分别进行积分即可得到式(1).
需要说明的是,式(2)通过以下步骤推导获得:
B1:假设仅存在拉伸流动,则对于收缩流道作用在流体微元段上的平均拉伸阻力和弹性压降为:
πr2dPce=-σcd(πr2) (10)
其中,σc为收缩流道中与x轴平行的平均拉伸应力,dPce为收缩流道微元段两端的弹性压降。
B2:平均拉伸应力与流速剖面位置无关,则平均拉伸应力为:
其中,为平均拉伸速率。
B3:在流体微元段取一个厚度为h,半径为的环面a,收缩角度为α1,其计算式为:
则对于这个收缩环面,其拉伸速度为:
其中,
其中,va为与环面中心线平行的速度。
B4:对于幂律流体:
将式(15)带入式(17),即可得到:
B5:因此,作用在整个收缩环面上的应力为:
而平均拉伸应力:
B6:将式(10)、式(16)带入式(19),即可得到:
B7:在区间[r1,r2]上对式(20)两端分别进行积分即可得到式(2)。
需要说明的是,式(3)通过以下步骤推导获得:
C1:假设仅存在剪切流动,不考虑弹性作用,则对于扩张流道作用在流体微元段上的黏性阻力和压降的平衡方程式为:
πr2dpds=-2πrdxsecα(τdcosβ) (21)
其中,τd为扩张流道中的剪切应力,dpds为扩张流道两端的剪切压降;sec为为正割函数,secα=1/cosα;
其中,
C2:假设剪切应力为剪切速率符合幂律关系,流体微元段的剪切速率r可近似用圆管中的剪切速率,则剪切应力可表示为:
C3:式(22)两端分别对x求导,可得式:
其中,将式(23)、(24)带入式(21):即可得到:
C4:在区间[r2,r3]上对式(25)两端分别进行积分即可得到式(3).
需要说明的是,式(4)通过以下步骤推导获得:
D1:假设仅存在拉伸流动,则对于扩张流道作用在流体微元段上的平均拉伸阻力和弹性压降为:
πr2dPde=-σd(πr2) (26)
其中,σd为收缩流道中与x轴平行的平均拉伸应力,dPde为扩张流道微元段两端的弹性压降。
D2:将式(19)、式(23)带入式(26)可得:
D3:在区间[r2,r3]上对式(27)两端分别进行积分即可得到式(4).
需要说明的是,所述数学模型通过以下步骤得到:
假设l1=l2,r1=r3,α=β,那么聚合物黏弹性流体在一个整体流动单元中的压降为:
则:由m个缩扩流道单元串联组成的毛细管两端的总压降,见式6-37。
在一些实施例中,所述S3中,所述启动压力实验数据用岩心样品通过瞬态动用法进行黏弹性流体驱替稠油的启动压力梯度测试得到。
需要说明的是,所述启动压力实验数据通过以下步骤获得:
S31:提供实验岩心;
S32:采用所述实验岩心进行饱和原油,束缚水建立后,静置,使所述实验岩心内的原油充分老化,并保证所述实验岩心的出口端管线内充满液体;
S33:设定驱替流量,进行岩心样品的驱替,观察岩心样品的出口端液体移动情况,当所述岩心样品的出口端液体开始移动时,记录压差,即为实验结果数据。
在一些实施例中,所述S3中,所述启动压力实验数据包括实验结果数据和实验参数;
所述实验参数包括岩心样品的长度。
在一些实施例中,所述S3中,多孔介质物性条件,根据目标油藏测试的渗透率和微观孔喉结构分析,明确渗透率与孔喉结构的对应关系,在气测渗透率的基础上,通过渗透率明确该条件下的多孔介质孔喉结构组成,即孔喉长度、宽度、夹角、数目范围等。
在一些实施例中,所述S4中,所述黏弹性流体的有效启动的压力梯度的值等于所述实验结果数据与附加阻力的差值与所述岩心样品的长度的比值。
以下以疏水缔合聚合物(AP-P4)溶液作为黏弹性流体为例,进行启动压力梯度的修订以获得有效启动的压力梯度。
S1:采用疏水缔合聚合物(AP-P4)溶液作为溶质,以目标油藏地层水为溶剂,配制浓度为2500mg/L的AP-P4溶液,静置后测试黏弹性流体溶液的流变数据,并得到流变曲线(如图3所示);根据流变数据进行流变规律拟合,确定应用流变Carreau-Yasuda模型(见式1),获取目标参数,见表1;
表1流变Carreau-Yasuda模型下的参数
S2:按照图1和图2所示的示意图,构建黏弹性流体通过多孔介质时所产生流动阻力的数学模型,具体步骤如下。
S21)不考虑弹性作用,仅考虑剪切流动的情况下,黏弹性流体在收缩流道中的黏性压降△Pcs,如式(2)。
S22)在收缩流道中仅考虑拉伸流动的情况下,黏弹性流体在收缩流道中的弹性压降△Pce,如式(3)。
S23)在扩张流道中仅考虑剪切流动的情况下,黏弹性流体在扩张流道中由剪切流动产生的压降△Pds,如式(4)。
S24)在扩张流道中仅考虑拉伸流动的情况下,黏弹性流体在扩张流道中的弹性压降△Pde,如式(5)。
最终形成(6)所示的总压降公式,m个收缩-扩张流道组成。
式中:△Ps是黏弹性溶液通过m个收缩-扩张流道组成孔喉结构时的总压降;△Pc是黏性压降;△Pd是弹性压降;q为流体通过微元段的流量;m是收缩-扩张流道组成个数;μe是弹性黏度;μv是黏性黏度;α是壁面与轴线的角度;n是幂率指数;r2是收缩道半径;r1是扩张道半径。
S3:用瞬态动用法进行黏弹性流体驱替稠油的启动压力梯度测试,具体实验步骤如下:首先提供长为8cm,直径为25mm,孔隙度为35%,渗透率为1805mD的实验岩心,并用黏度为6429mPa·s,温度25℃原油进行饱和,建立好束缚水后,静置24h,使岩心内的原油充分老化,并保证出口端管线内充满液体;然后通过泵设定微小的流量驱替,在岩心入口端以较为缓慢的速度建立压差,并利用液柱高度计量压差,仔细观察出口端液体移动情况,当液体开始移动时,记录此时的液柱高度;获取启动压力实验数据为0.00438MPa;
然后,根据多孔介质物性参数,获得到的对应孔喉条件尺寸、长=0.00091mm、宽=0.000023mm、夹角α=42°、夹角β=37°、m=8791,分析计算聚合物溶液注入岩心中产生的附加阻力为0.0031MPa。
S4:将S3中索取的启动压力实验数据(0.00438MPa)减去附加阻力(0.0031MPa),然后获得的压差值除以S3中的岩心长度,即获得启动压力梯度为0.016MPa/m。这样获得启动压力梯度值即为稠油在多孔介质中被黏弹性流体有效启动的压力梯度。
以上所述,并非对本发明作任何形式上的限制,虽然本发明已通过上述实施例揭示,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,可利用上述揭示的技术内容作出些变动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (9)
1.一种应用于黏弹性流体测试启动压力梯度的修订方法,其特征在于,包括以下步骤:
S1、提供黏弹性流体,测得所述黏弹性流体溶液的流变数据;根据Carreau-Yasuda模型和所述流变数据,获取所述黏弹性流体的目标参数;
S2、构建黏弹性流体通过多孔介质时所产生流动阻力的数学模型;
所述S2中,所述数学模型为:
其中,△P为黏弹性溶液通过m个收缩-扩张流道组成孔喉结构时的总压降;△Pc为黏性压降;△Pd为弹性压降;q为流体通过微元段的流量;m为收缩-扩张流道组成个数;μe为弹性黏度;μv为黏性黏度;α为流道的壁面与轴线的角度;n为幂率指数;r2为收缩道半径;r1为扩张道半径;
S3、提供所述黏弹性流体的启动压力实验数据,根据所述启动压力实验数据、所述目标参数和和所述数学模型,得到附加阻力;
S4、根据所述附加压力和启动压力实验数据获得所述黏弹性流体的有效启动的压力梯度。
2.根据权利要求1所述的一种应用于黏弹性流体测试启动压力梯度的修订方法,其特征在于,所述S1中,所述黏弹性流体的流变数据通过流变仪测试得到;
所述流变数据为剪切区间范围不小于0.01-10000s-1的测试数据。
3.根据权利要求1所述的一种应用于黏弹性流体测试启动压力梯度的修订方法,其特征在于,所述S1中,所述Carreau-Yasuda模型为:
其中,μ为不同剪切速率下的表观黏度,mPa·s;μ0为零剪切黏度,mPa·s;μinf为极限剪切黏度,mPa·s;λ为松弛时间,n为幂率指数,a为常数,x为剪切速率。
4.根据权利要求1所述的一种应用于黏弹性流体测试启动压力梯度的修订方法,其特征在于,所述S2中,所述数学模型通过以下步骤得到:
S21:建立包括m个收缩流道和m个扩张流道的流道模型,进行黏弹性流体在所述流道模型中的流通模拟;其中m为大于或等于1的自然数;
S22:假设仅存在剪切流动,建立所述黏弹性流体在所述收缩流道中的黏性压降和扩张流道中的剪切压降;
假设仅存在拉伸流动,建立所述黏弹性流体在所述收缩流道中的弹性压降和在所述扩张流道中的弹性压降;
S23:根据所述黏弹性流体在所述收缩流道中的黏性压降、在所述扩张流道中的剪切压降、在所述收缩流道中的弹性压降和在所述扩张流道中的弹性压降,得到所述数学模型。
5.根据权利要求4所述的一种应用于黏弹性流体测试启动压力梯度的修订方法,其特征在于,所述黏弹性流体在所述收缩流道中的黏性压降为:
所述黏弹性流体在所述扩张流道中的剪切压降为:
所述黏弹性流体在所述收缩流道中的弹性压降为:
其中,β为收缩流道中的外径与内径之间的夹角;
所述黏弹性流体在所述扩张流道中的弹性压降为:
6.根据权利要求5所述的一种应用于黏弹性流体测试启动压力梯度的修订方法,其特征在于,所述收缩流道为内径逐渐变小的变径管道;
所述扩张流道为内径逐渐变大的变径管道;
所述收缩流道的内径较小一端与所述扩张流道的内径较小一端连接,且所述收缩流道的内径较小的一端的内径等于所述扩张流道的内径较小一端的内径、所述收缩流道的内径较大的一端的内径等于所述扩张流道的内径较大一段的内径。
7.根据权利要求1所述的一种应用于黏弹性流体测试启动压力梯度的修订方法,其特征在于,所述S3中,所述启动压力实验数据用岩心样品通过瞬态动用法进行黏弹性流体驱替稠油的启动压力梯度测试得到。
8.根据权利要求7所述的一种应用于黏弹性流体测试启动压力梯度的修订方法,其特征在于,所述S3中,所述启动压力实验数据包括实验结果数据和实验参数;
所述实验参数包括岩心样品的长度。
9.根据权利要求8所述的一种应用于黏弹性流体测试启动压力梯度的修订方法,其特征在于,所述S4中,所述黏弹性流体的有效启动的压力梯度的值等于所述实验结果数据与附加阻力的差值与所述岩心样品的长度的比值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111195761.3A CN113916747B (zh) | 2021-10-14 | 2021-10-14 | 一种应用于黏弹性流体测试启动压力梯度的修订方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111195761.3A CN113916747B (zh) | 2021-10-14 | 2021-10-14 | 一种应用于黏弹性流体测试启动压力梯度的修订方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113916747A CN113916747A (zh) | 2022-01-11 |
CN113916747B true CN113916747B (zh) | 2023-09-15 |
Family
ID=79240254
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111195761.3A Active CN113916747B (zh) | 2021-10-14 | 2021-10-14 | 一种应用于黏弹性流体测试启动压力梯度的修订方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113916747B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116087053B (zh) * | 2022-12-15 | 2024-05-07 | 长江大学 | 一种研究稠油启动压力梯度对剩余油分布影响的方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08116945A (ja) * | 1994-10-27 | 1996-05-14 | Rheon Autom Mach Co Ltd | 食品材料の間欠定量吐出装置 |
EP0868976A2 (en) * | 1997-03-06 | 1998-10-07 | Keltech Engineering | Lapping apparatus and method for high speed lapping with a rotatable abrasive platen |
EP2088281A1 (en) * | 2008-02-08 | 2009-08-12 | Bp Exploration Operating Company Limited | Method for determining the viscosity of a waxy fluid in a shut-in pipeline |
CN106596336A (zh) * | 2016-12-13 | 2017-04-26 | 东北石油大学 | 黏弹性聚合物在多孔介质中弹性湍流特性测定装置及方法 |
KR20170088166A (ko) * | 2016-01-22 | 2017-08-01 | 한국과학기술원 | 점탄성유체와 점성유체의 병행층류를 이용한 미소 입자 분리 및 세정 방법 |
CN107503721A (zh) * | 2017-08-18 | 2017-12-22 | 西南石油大学 | 一种有效改善高渗高含水油藏聚合物驱效果的方法 |
CN109339774A (zh) * | 2018-10-18 | 2019-02-15 | 中国石油化工股份有限公司 | 基于边界层的低渗透变形介质非线性渗流数值模拟方法 |
DE102018209386B3 (de) * | 2018-06-13 | 2019-10-10 | Conti Temic Microelectronic Gmbh | Verfahren zum Bestimmen eines Fluiddrucks |
CN111929219A (zh) * | 2020-08-12 | 2020-11-13 | 西南石油大学 | 一种页岩油藏油水两相相对渗透率计算方法 |
CN112129682A (zh) * | 2020-09-04 | 2020-12-25 | 中国石油天然气股份有限公司 | 一种基于边界层动态变化的低渗透油藏驱替介质优选方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102010008146B4 (de) * | 2010-02-12 | 2022-03-31 | Carl Zeiss Meditec Ag | Messsystem und Verfahren zum Ermitteln des Innendrucks eines Auges sowie Verfahren und System zum Einstellen des Augeninnendrucks |
-
2021
- 2021-10-14 CN CN202111195761.3A patent/CN113916747B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08116945A (ja) * | 1994-10-27 | 1996-05-14 | Rheon Autom Mach Co Ltd | 食品材料の間欠定量吐出装置 |
EP0868976A2 (en) * | 1997-03-06 | 1998-10-07 | Keltech Engineering | Lapping apparatus and method for high speed lapping with a rotatable abrasive platen |
EP2088281A1 (en) * | 2008-02-08 | 2009-08-12 | Bp Exploration Operating Company Limited | Method for determining the viscosity of a waxy fluid in a shut-in pipeline |
KR20170088166A (ko) * | 2016-01-22 | 2017-08-01 | 한국과학기술원 | 점탄성유체와 점성유체의 병행층류를 이용한 미소 입자 분리 및 세정 방법 |
CN106596336A (zh) * | 2016-12-13 | 2017-04-26 | 东北石油大学 | 黏弹性聚合物在多孔介质中弹性湍流特性测定装置及方法 |
CN107503721A (zh) * | 2017-08-18 | 2017-12-22 | 西南石油大学 | 一种有效改善高渗高含水油藏聚合物驱效果的方法 |
DE102018209386B3 (de) * | 2018-06-13 | 2019-10-10 | Conti Temic Microelectronic Gmbh | Verfahren zum Bestimmen eines Fluiddrucks |
CN109339774A (zh) * | 2018-10-18 | 2019-02-15 | 中国石油化工股份有限公司 | 基于边界层的低渗透变形介质非线性渗流数值模拟方法 |
CN111929219A (zh) * | 2020-08-12 | 2020-11-13 | 西南石油大学 | 一种页岩油藏油水两相相对渗透率计算方法 |
CN112129682A (zh) * | 2020-09-04 | 2020-12-25 | 中国石油天然气股份有限公司 | 一种基于边界层动态变化的低渗透油藏驱替介质优选方法 |
Non-Patent Citations (5)
Title |
---|
Linear and Nonlinear Viscoelastic Modeling of Aorta and Carotid Pressure-Area Dynamics Under in Vivo and Ex Vivo Conditions;Valdez-Jasso 等;《ANNALS OF BIOMEDICAL ENGINEERING》;第39卷(第05期);第1438-1456页 * |
Numerical simulation of a bubble rising in shear-thinning fluids;Zhang L 等;《Journal of Non-newtonian Fluid Mechanics》;第165卷(第11期);第555-567页 * |
扩展流道内修正上随体MaxweII粘弹性流体的流动;尹洪军 等;《特种油气藏》(第04期);第36-39+105-106页 * |
聚集行为对聚合物溶液性能的影响研究;朱诗杰 等;《应用化工》;第49卷(第03期);第597-601页 * |
聚驱井流入动态特性研究;代兆国;《中国优秀硕士学位论文全文数据库工程科技I辑》(第07期);第B019-174页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113916747A (zh) | 2022-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110907334B (zh) | 一种砾岩全直径岩心径向流油水相对渗透率测量装置及方法 | |
Zeng et al. | Low velocity non-linear flow in ultra-low permeability reservoir | |
CN109557010B (zh) | 一种测定渗流启动压力梯度的实验方法 | |
CN104237099B (zh) | 测定致密岩心径向渗透率的装置及方法 | |
Shubao et al. | Dynamic effect of capillary pressure in low permeability reservoirs | |
CN101968423B (zh) | 低渗透储层启动压力测试方法 | |
CN102866082B (zh) | 一种测量泡沫油流变性的装置及方法 | |
CN105804726B (zh) | 一种泡点压力测试装置及方法 | |
Dong et al. | Effect of dynamic pseudo threshold pressure gradient on well production performance in low-permeability and tight oil reservoirs | |
CN107525746A (zh) | 一种表征致密油藏低速非达西渗流特征的方法及装置 | |
CN113866069B (zh) | 一种页岩岩心渗透率实验装置和方法 | |
CN113916747B (zh) | 一种应用于黏弹性流体测试启动压力梯度的修订方法 | |
CN112557276B (zh) | 一种同时测量多孔介质渗透率和孔隙率的方法 | |
CN107725046A (zh) | 一种评价油藏注水过程中毛管力的设备和方法 | |
Fuquan et al. | Single-and two-phase flow model in low-permeability reservoir | |
An et al. | Non-Klinkenberg slippage phenomenon at high pressure for tight core floods using a novel high pressure gas permeability measurement system | |
CN111595731B (zh) | 一种驱油用非均相粘弹性颗粒溶液阻力系数测试系统及测试方法 | |
CN114577677B (zh) | 基于储层非均质性的致密油藏非线性渗流表征方法 | |
CN202735216U (zh) | 一种钻井液半透膜评价设备 | |
Yang et al. | A new measurement method for radial permeability and porosity of shale | |
Fang et al. | Experimental study of gas flow characteristics in micro-/nano-pores in tight and shale reservoirs using microtubes under high pressure and low pressure gradients | |
CN109556996B (zh) | 油水两相干扰压力梯度的测量方法 | |
Haobo et al. | A comprehensive hydraulic calculation method of non-newtonian fluids used fourparameter model | |
Chen et al. | A fractal model of low-velocity non-Darcy flow considering viscosity distribution and boundary layer effect | |
CN115248177A (zh) | 一种基于光纤传感的低渗透岩石突破压力的测量方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |