CN112012728A - 一种聚合物在储层隙间流速的确定方法 - Google Patents
一种聚合物在储层隙间流速的确定方法 Download PDFInfo
- Publication number
- CN112012728A CN112012728A CN202010835220.1A CN202010835220A CN112012728A CN 112012728 A CN112012728 A CN 112012728A CN 202010835220 A CN202010835220 A CN 202010835220A CN 112012728 A CN112012728 A CN 112012728A
- Authority
- CN
- China
- Prior art keywords
- viscosity
- polymer
- shear
- flow rate
- polymer solution
- 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.)
- Granted
Links
- 229920000642 polymer Polymers 0.000 title claims abstract description 133
- 238000000034 method Methods 0.000 title claims abstract description 46
- 230000003715 interstitial flow Effects 0.000 claims abstract description 22
- 238000012937 correction Methods 0.000 claims abstract description 18
- 238000004088 simulation Methods 0.000 claims abstract description 14
- 238000002474 experimental method Methods 0.000 claims abstract description 12
- 238000013461 design Methods 0.000 claims abstract description 9
- 238000012360 testing method Methods 0.000 claims abstract description 8
- 238000010586 diagram Methods 0.000 claims abstract description 7
- 239000012530 fluid Substances 0.000 claims description 16
- 230000035699 permeability Effects 0.000 claims description 14
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 13
- 238000010008 shearing Methods 0.000 claims description 12
- 239000011148 porous material Substances 0.000 claims description 11
- 239000008398 formation water Substances 0.000 claims description 10
- HJUFTIJOISQSKQ-UHFFFAOYSA-N fenoxycarb Chemical compound C1=CC(OCCNC(=O)OCC)=CC=C1OC1=CC=CC=C1 HJUFTIJOISQSKQ-UHFFFAOYSA-N 0.000 claims description 8
- 239000000126 substance Substances 0.000 claims description 8
- 238000000518 rheometry Methods 0.000 claims description 6
- 238000013178 mathematical model Methods 0.000 claims description 3
- 230000000704 physical effect Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 abstract description 3
- 239000000243 solution Substances 0.000 description 52
- 238000004364 calculation method Methods 0.000 description 4
- 230000002349 favourable effect Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000033558 biomineral tissue development Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000008719 thickening Effects 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Geochemistry & Mineralogy (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Theoretical Computer Science (AREA)
- Geophysics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种聚合物在储层隙间流速的确定方法，具体步骤为：步骤一，对任意开展聚合物驱的油藏，采集数据开展聚合物流变性实验，测试不同浓度、不同剪切速率下聚合物溶液的粘度；步骤二，给定聚合物溶液浓度条件下，建立剪切速率与聚合物溶液粘度的散点图，采用幂率模式拟合粘度和剪切速率数据，求出幂律指数；步骤三，通过步骤二得到的幂律指数确定Rabinowitsch校正系数，将Rabinowitsch校正系数代入剪切速率与隙间流速的关系式，生成给定浓度的聚合物溶液的隙间流速公式；本发明建立了聚合物溶液粘度与储层隙间流速的转换方法，聚合物溶液粘度与储层隙间流速关系数据是聚合物驱工程设计和数值模拟动态预测必不可少的重要依据。
Description
技术领域
本发明属于油气田开发领域，具体涉及一种聚合物在储层隙间流速的确定方法。
背景技术
聚合物溶液的流变性是聚合物驱所涉及的一个重要基础领域。由于聚合物溶液在油层孔隙介质中的流动行为及其传输性质受其本身流变性质制约，因而聚合物溶液的流变参数是矿场聚合物驱工程设计和数值模拟动态预测必不可少的重要依据。
在低流速下,聚合物的相对分子质量越大，溶液浓度越高,聚合物驱阻力系数越大；在高流速下,受相对分子质量、聚合物溶液浓度、注入速度的综合影响,且不同的聚合物溶液表现出不同的特征。
通常，聚合物流变性实验得到的是剪切速率与粘度之间的关系，部分ECLIPSE数值模拟软件需要输入流速与粘度之间的关系，即使目前的关键字支持剪切速率直接输入，也没有零剪切粘度计算和剪切变稠方面的考虑。
目前，CMG化学驱数值模拟部分可实现流变性实验数据的直接输入，但是在不同油藏，特别是低渗透油藏，用不同岩心驱替实验得出的阻力系数和残余阻力系数，不可及孔隙体积等数据需要同剪切速率和流速转化数据相统一。
数值模拟软件提供了相对简单的实验数据输入方式，但仍然需要找到一种便捷的聚合物在储层中隙间流速方法，有助于给数值模拟软件更加精确的数据，从而使得结果更加准确。
发明内容
为了克服现有不能精确的确定聚合物在储层中隙间流速，从而导致后续矿场聚合物驱工程设计和数值模拟动态预测出现偏差的问题，本发明提供一种聚合物在储层中隙间流速的确定方法，本发明便捷且易掌握，将剪切速率转化为地层中的流速；确定不同浓度聚合物的零剪切粘度，流速与粘度之间的关系数据，可直接用于输入ECLIPSE数值模拟软件进行化学驱数值模拟，有助于给数值模拟软件更加精确的数据，从而使得结果更加准确。
本发明采用的技术方案是：
一种聚合物在储层隙间流速的确定方法，具体步骤为：
步骤一，对任意开展聚合物驱的油藏，采集数据开展聚合物流变性实验，测试不同浓度、不同剪切速率下聚合物溶液的粘度；
步骤二，给定聚合物溶液浓度条件下，建立剪切速率与聚合物溶液粘度的散点图，采用幂率模式拟合粘度和剪切速率数据，求出幂律指数；
步骤三，通过步骤二得到的幂律指数确定Rabinowitsch校正系数，将Rabinowitsch校正系数代入剪切速率与隙间流速的关系式，生成给定浓度的聚合物溶液的隙间流速公式；
步骤四，在步骤三中生成的给定浓度的聚合物溶液的隙间流速公式中输入开展聚合物驱的油藏的平均孔隙度、渗透率，确定给定浓度的聚合物溶液在不同剪切速率条件下对应的隙间流速；
步骤五，通过规划求解方法对步骤四中的实验数据进行拟合，确定给定浓度聚合物在零剪切速率条件下所对应的粘度；
步骤六，重复步骤五，得出不同浓度下聚合物溶液在油藏物性条件下的隙间流速与粘度的对应关系数据；
步骤七，将步骤六得到的隙间流速与粘度的对应关系，输入化学驱数值模型，进行后续矿场聚合物驱工程设计和数值模拟动态预测。
所述的步骤二中，采用表格软件中建立剪切速率与粘度数据的散点图，采用幂率模式拟合粘度和剪切速率数据，求出幂律指数；表格软件可为excel或WPS。
所述的步骤二中，幂率指数的确定，具体方法为：
聚合物水溶液倾向于牛顿流体行为时，即溶液粘度与剪切速率无关；聚合物溶液呈假塑性或剪切变稀行为时，这种状态下采用幂律模型来描述：
式(1)中：
γ——剪切应力,mPa
K——稠度系数,mPa.sn
n——幂律指数或流动行为指数，无因次由于粘度是剪切应力除以剪切速率，那么幂律流体的粘度则为：
式(2)中：η——粘度，mPa.s或cP
通过测定不同剪切速率下聚合物溶液的粘度，采用上述指数函数拟合曲线可求得幂律指数；
实际操作中，用目标油藏地层水配制给定浓度的聚合物溶液，在固定温度条件下，测定不同剪切速率下聚合物溶液的粘度，用幂律模式拟合粘剪曲线，可求得幂律指数n。
所述的步骤四中，储层隙间流速的确定方法为：
为了估算非牛顿流体在孔隙介质中流动时的平均剪切速率，在最为通用的表达式中，把孔隙体系模拟为具有相同半径和长度的毛细管束，因而:
式(3)中:r为等效毛细管半径，cm；φ是孔隙度；K为聚合物流动前水的渗透率,cm2；
把毛细管壁处的剪切速率取作
将式(3)代人替代r, 并对幂律流体应用非牛顿校正，将孔隙介质内的剪切速率可改写为：
式中:
参照式(4)可将剪切速率转变为ECLIPSE数值模拟软件中要求输入的流速
根据式(4)要求，我们需用通过聚合物溶液流变性实验，得出剪切速率与粘度之间的关系，拟合求出幂率指数，进而计算 Rabinowitsch校正系数；通过在式(4)中代入储层孔隙度和渗透率，确定出该聚合物在储层的隙间流速。
所述的步骤五中，零剪切粘度的确定方法为：
聚合物溶液的流变性试验所对应的数理方程是Meter模式：
式中：
μp-------剪切速率为
时聚合物溶液的粘度，mPa.s
μw-------剪切速率为
时水的粘度，mPa.s
μ0-------零剪切粘度，mPa.s
μ∞-------无限剪切粘度，mPa.s
pα-------指数系数
令：
Meter模式可转化为一元线性方程y＝a+bx的形式，对于该形式的数学模型依照最小二乘法原理应有：
根据目标油藏流体性质，近似取值无限剪切粘度和水的粘度，应用表格软件中的规划求解方法对实验数据进行拟合，设置零剪切粘度为自变量，设置式(8)中的Q为目标值，在目标值Q尽可能接近于0的条件下，求得零剪切粘度。
所述的步骤一中，采集的数据至少包括地层水矿化度、油藏温度、平均孔隙度和渗透率。
所述的步骤一中，在地层水矿化度和油藏温度条件下配制聚合物溶液。
本发明的有益效果为：
本发明便捷且易掌握，将剪切速率转化为地层中的流速；确定不同浓度聚合物的零剪切粘度，流速与粘度之间的关系数据，可直接用于输入ECLIPSE数值模拟软件进行化学驱数值模拟，有助于给数值模拟软件更加精确的数据，从而使得结果更加准确。
以下将结合附图进行进一步的说明。
附图说明
图1为给定浓度聚合物粘度与剪切速率的关系曲线及幂率指数的拟合求取。
具体实施方式
实施例1：
为了克服现有不能精确的确定聚合物在储层中隙间流速，从而导致后续矿场聚合物驱工程设计和数值模拟动态预测出现偏差的问题，本发明提供如图1所示的一种聚合物在储层中隙间流速的确定方法，本发明便捷且易掌握，将剪切速率转化为地层中的流速；确定不同浓度聚合物的零剪切粘度，流速与粘度之间的关系数据，可直接用于输入ECLIPSE数值模拟软件进行化学驱数值模拟，有助于给数值模拟软件更加精确的数据，从而使得结果更加准确。
一种聚合物在储层隙间流速的确定方法，具体步骤为：
步骤一，对任意开展聚合物驱的油藏，采集数据开展聚合物流变性实验，测试不同浓度、不同剪切速率下聚合物溶液的粘度；
步骤二，给定聚合物溶液浓度条件下，建立剪切速率与聚合物溶液粘度的散点图，采用幂率模式拟合粘度和剪切速率数据，求出幂律指数；
步骤三，通过步骤二得到的幂律指数确定Rabinowitsch校正系数，将Rabinowitsch校正系数代入剪切速率与隙间流速的关系式，生成给定浓度的聚合物溶液的隙间流速公式；
步骤四，在步骤三中生成的给定浓度的聚合物溶液的隙间流速公式中输入开展聚合物驱的油藏的平均孔隙度、渗透率，确定给定浓度的聚合物溶液在不同剪切速率条件下对应的隙间流速；
步骤五，通过规划求解方法对步骤四中的实验数据进行拟合，确定给定浓度聚合物在零剪切速率条件下所对应的粘度；
步骤六，重复步骤五，得出不同浓度下聚合物溶液在油藏物性条件下的隙间流速与粘度的对应关系数据；
步骤七，将步骤六得到的隙间流速与粘度的对应关系，输入化学驱数值模型，进行后续矿场聚合物驱工程设计和数值模拟动态预测。
本发明的目的是提供一种聚合物在储层中隙间流速的确定方法，针对于目前室内试验给出一般是剪切速率与粘度之间的关系，将剪切速率转化为地层中的流速，同时计算不同浓度聚合物的零剪切粘度，可直接用于输入数值模拟软件进行模拟，同时便于油藏工程计算。
本发明建立了聚合物溶液粘度与储层隙间流速的转换方法，聚合物溶液粘度与储层隙间流速关系数据是聚合物驱工程设计和数值模拟动态预测必不可少的重要依据。
实施例2：
基于实施例1的基础上，本实施例中，优选地，所述的步骤二中，采用表格软件中建立剪切速率与粘度数据的散点图，采用幂率模式拟合粘度和剪切速率数据，求出幂律指数；表格软件可为excel或WPS。
优选地，所述的步骤二中，幂率指数的确定，具体方法为：
聚合物溶液的流变性是相当复杂的，聚合物水溶液倾向于牛顿流体行为时，即溶液粘度与剪切速率无关；聚合物溶液呈假塑性或剪切变稀行为时，这种状态下采用幂律模型来描述：
式(1)中：
γ——剪切应力,mPa
K——稠度系数,mPa.sn
n——幂律指数或流动行为指数，无因次由于粘度是剪切应力除以剪切速率，那么幂律流体的粘度则为：
式(2)中：η——粘度，mPa.s或cP
通过测定不同剪切速率下聚合物溶液的粘度，采用上述指数函数拟合曲线可求得幂律指数；
实际操作中，用目标油藏地层水配制给定浓度的聚合物溶液，在固定温度条件下，测定不同剪切速率下聚合物溶液的粘度，用幂律模式拟合粘剪曲线，可求得幂律指数n。
优选地，所述的步骤四中，储层隙间流速的确定方法为：
为了估算非牛顿流体在孔隙介质中流动时的平均剪切速率，在最为通用的表达式中，把孔隙体系模拟为具有相同半径和长度的毛细管束，因而:
式(3)中:r为等效毛细管半径，cm；φ是孔隙度；K为聚合物流动前水的渗透率,cm2(1cm2＝1011mD)；
把毛细管壁处的剪切速率取作
将式(3)代人替代r, 并对幂律流体应用非牛顿校正，Jenning等人给出的孔隙介质内的剪切速率可改写为：
式中:
参照式(4)可将剪切速率转变为ECLIPSE数值模拟软件中要求输入的流速
根据式(4)要求，我们需用通过聚合物溶液流变性实验，得出剪切速率与粘度之间的关系，拟合求出幂率指数，进而计算 Rabinowitsch校正系数；通过在式(4)中代入储层孔隙度和渗透率，确定出该聚合物在储层的隙间流速。
优选地，所述的步骤五中，零剪切粘度的确定方法为：
数值模拟计算需要的剪切数据需要从零开始，零剪切粘度实验数据无法得到，必须通过拟合才能得到。
聚合物溶液的流变性试验所对应的数理方程是Meter模式：
式中：
μp-------剪切速率为
时聚合物溶液的粘度，mPa.s
μw-------剪切速率为
时水的粘度，mPa.s
μ0-------零剪切粘度，mPa.s
μ∞-------无限剪切粘度，mPa.s
pα-------指数系数(实验确定)
令：
Meter模式可转化为一元线性方程y＝a+bx的形式，对于该形式的数学模型依照最小二乘法原理应有：
根据目标油藏流体性质，近似取值无限剪切粘度和水的粘度，应用表格软件中的规划求解方法对实验数据进行拟合，设置零剪切粘度为自变量，设置式(8)中的Q为目标值，在目标值Q尽可能接近于0的条件下，求得零剪切粘度。
优选地，所述的步骤一中，采集的数据至少包括地层水矿化度、油藏温度、平均孔隙度和渗透率。
优选地，所述的步骤一中，在地层水矿化度和油藏温度条件下配制聚合物溶液。
本发明可为在低渗透油藏化学驱时不同粘度聚合物在储层中隙间流速的确定方法。利用建立剪切速率与储层隙间流速的关系式，利用Meter幂率模式中的幂率指数，通过Rabinowitsch校正系数校正非牛顿流体，将Jenning给出的孔隙介质内的剪切速率可改写为隙间流速，得出不同浓度、不同隙间流速下聚合物溶液的粘度，直接输入 ECLIPSE数值模拟软件进行模拟计算。
实施例3：
基于实施例1或2的基础上，本实施例中对油藏A进行一种聚合物在储层隙间流速的确定方法，具体过程为：
步骤1、流变性实验
用模拟目标油藏A矿化度的地层水配制不同浓度的聚合物溶液，完全溶解后用布氏粘度计分别在目标油藏A温度下测定不同剪切速率下的溶液粘度。转子：0号；转速：0.3～30rpm，参见附表1。
步骤2、计算幂率指数
选取配制浓度为1500mg/L聚合物溶液在油藏A温度条件下测定的不同剪切速率下聚合物溶液的粘度，用幂律模式的式(2)公式拟合粘剪曲线，求得幂率指数n为0.6564，参见图1。
步骤3、计算隙间流速
步骤2求得油藏A条件下，聚合物浓度为1500mg/L时，幂律流动行为指数n为0.6564，实验用岩心φ＝13.9％，水测渗透率k＝3.56 ×10-3μm2,参数代入式(3)公式，得到不同剪切速率对应的隙间流速值，参见附表2。
步骤4、零剪切粘度
油藏A条件下，近似取无限剪切粘度值为地层水的粘度0.5mPa.s，忽略剪切速率对地层水的粘度影响，近似取地层水水在不同剪切速率的下的粘度为0.5mPa.s。选择1500mg/L浓度条件下实验数据，初设零剪切粘度为35mPa.s，应用“规划求解”方法对实验数据进行拟合，设置零剪切粘度为自变量，设置式(8)中的Q为目标值，求得最终的1500浓度条件下的零剪切粘度为53.632mPa.s，进而得到1500mg/L聚合物隙间流速与粘度关系，参见附表3。
步骤5、不同流速与粘度数据
重复步骤4，计算出不同浓度聚合物溶液的零剪切粘度，参见附表4。进而得到在计算出不同浓度聚合物溶液隙间流速与粘度的关系。
表1油藏A条件下不同剪切速率下聚合物粘浓关系
表2 1500mg/L浓度聚合物剪切速率与隙间流速对照表
表3 1500mg/L浓度聚合物隙间流速与粘度对照表
表4油藏A条件下不同剪切速率下聚合物粘浓关系
以上举例仅是对本发明的说明，并不构成保护范围限制，凡是与本发明相同或似的设计均属于发明保护范围之内。实施例没有详细叙述的方法步骤及实验过程均属于本行业公知常用或手段，这里将不再一一叙述。
Claims (7)
1.一种聚合物在储层隙间流速的确定方法，其特征在于：具体步骤为：
步骤一，对任意开展聚合物驱的油藏，采集数据开展聚合物流变性实验，测试不同浓度、不同剪切速率下聚合物溶液的粘度；
步骤二，给定聚合物溶液浓度条件下，建立剪切速率与聚合物溶液粘度的散点图，采用幂率模式拟合粘度和剪切速率数据，求出幂律指数；
步骤三，通过步骤二得到的幂律指数确定Rabinowitsch校正系数，将Rabinowitsch校正系数代入剪切速率与隙间流速的关系式，生成给定浓度的聚合物溶液的隙间流速公式；
步骤四，在步骤三中生成的给定浓度的聚合物溶液的隙间流速公式中输入开展聚合物驱的油藏的平均孔隙度、渗透率，确定给定浓度的聚合物溶液在不同剪切速率条件下对应的隙间流速；
步骤五，通过规划求解方法对步骤四中的实验数据进行拟合，确定给定浓度聚合物在零剪切速率条件下所对应的粘度；
步骤六，重复步骤五，得出不同浓度下聚合物溶液在油藏物性条件下的隙间流速与粘度的对应关系数据；
步骤七，将步骤六得到的隙间流速与粘度的对应关系，输入化学驱数值模型，进行后续矿场聚合物驱工程设计和数值模拟动态预测。
2.根据权利要求1所述的一种聚合物在储层隙间流速的确定方法，其特征在于：所述的步骤二中，采用表格软件中建立剪切速率与粘度数据的散点图，采用幂率模式拟合粘度和剪切速率数据，求出幂律指数；表格软件可为excel或WPS。
3.根据权利要求1所述的一种聚合物在储层隙间流速的确定方法，其特征在于：所述的步骤二中，幂率指数的确定，具体方法为：
聚合物水溶液倾向于牛顿流体行为时，即溶液粘度与剪切速率无关；聚合物溶液呈假塑性或剪切变稀行为时，这种状态下采用幂律模型来描述：
式(1)中：
γ——剪切应力,mPa
K——稠度系数,mPa.sn
n——幂律指数或流动行为指数，无因次
由于粘度是剪切应力除以剪切速率，那么幂律流体的粘度则为：
式(2)中：η——粘度，mPa.s或cP
通过测定不同剪切速率下聚合物溶液的粘度，采用上述指数函数拟合曲线可求得幂律指数；
实际操作中，用目标油藏地层水配制给定浓度的聚合物溶液，在固定温度条件下，测定不同剪切速率下聚合物溶液的粘度，用幂律模式拟合粘剪曲线，可求得幂律指数n。
4.根据权利要求1所述的一种聚合物在储层隙间流速的确定方法，其特征在于：所述的步骤一中，采集的数据至少包括地层水矿化度、油藏温度、平均孔隙度和渗透率。
5.根据权利要求1所述的一种聚合物在储层隙间流速的确定方法，其特征在于：所述的步骤四中，储层隙间流速的确定方法为：
为了估算非牛顿流体在孔隙介质中流动时的平均剪切速率，在最为通用的表达式中，把孔隙体系模拟为具有相同半径和长度的毛细管束，因而:
式(3)中:r为等效毛细管半径，cm；φ是孔隙度；K为聚合物流动前水的渗透率,cm2；
把毛细管壁处的剪切速率取作
将式(3)代人替代r,并对幂律流体应用非牛顿校正，将孔隙介质内的剪切速率可改写为：
式中:
参照式(4)可将剪切速率转变为ECLIPSE数值模拟软件中要求输入的流速
根据式(4)要求，我们需用通过聚合物溶液流变性实验，得出剪切速率与粘度之间的关系，拟合求出幂率指数，进而计算Rabinowitsch校正系数；通过在式(4)中代入储层孔隙度和渗透率，确定出该聚合物在储层的隙间流速。
6.根据权利要求1所述的一种聚合物在储层隙间流速的确定方法，其特征在于：所述的步骤五中，零剪切粘度的确定方法为：
聚合物溶液的流变性试验所对应的数理方程是Meter模式：
式中：
μp-------剪切速率为
时聚合物溶液的粘度，mPa.s
μw-------剪切速率为
时水的粘度，mPa.s
μ0-------零剪切粘度，mPa.s
μ∞-------无限剪切粘度，mPa.s
pα-------指数系数
令：
Meter模式可转化为一元线性方程y＝a+bx的形式，对于该形式的数学模型依照最小二乘法原理应有：
根据目标油藏流体性质，近似取值无限剪切粘度和水的粘度，应用表格软件中的规划求解方法对实验数据进行拟合，设置零剪切粘度为自变量，设置式(8)中的Q为目标值，在目标值Q尽可能接近于0的条件下，求得零剪切粘度。
7.根据权利要求1所述的一种聚合物在储层隙间流速的确定方法，其特征在于：所述的步骤一中，在地层水矿化度和油藏温度条件下配制聚合物溶液。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010835220.1A CN112012728B (zh) | 2020-08-19 | 2020-08-19 | 一种聚合物在储层隙间流速的确定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010835220.1A CN112012728B (zh) | 2020-08-19 | 2020-08-19 | 一种聚合物在储层隙间流速的确定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112012728A true CN112012728A (zh) | 2020-12-01 |
| CN112012728B CN112012728B (zh) | 2024-03-26 |
Family
ID=73505016
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010835220.1A Active CN112012728B (zh) | 2020-08-19 | 2020-08-19 | 一种聚合物在储层隙间流速的确定方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112012728B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115600052A (zh) * | 2022-12-15 | 2023-01-13 | 昆明理工大学（Cn） | 一种确定幂律水泥浆液稠度系数的方法 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
| US20090293594A1 (en) * | 2008-05-30 | 2009-12-03 | Lincoln Industrial Corporation | System and Method for Estimating Apparent Viscosity of a Non-Newtonian Fluid |
| JP2011069754A (ja) * | 2009-09-28 | 2011-04-07 | Mitsubishi Rayon Co Ltd | せん断粘度推算方法および流動曲線の作成方法 |
| US20120158378A1 (en) * | 2010-12-15 | 2012-06-21 | Guillaume Enchery | Method of developing a petroleum reservoir from a flow model calibrated through pressure and saturation map scaling |
| CN103775078A (zh) * | 2014-02-14 | 2014-05-07 | 武汉科技大学 | 幂律流体在多孔介质中的有效粘度的预测方法 |
| US20180023374A1 (en) * | 2015-02-03 | 2018-01-25 | Schlumberger Technology Corporation | Multi-Phase Polymer Shear Viscosity Calculation in Polymer Coreflood Simulation Study Workflow |
| CN107869338A (zh) * | 2017-11-03 | 2018-04-03 | 中国石油化工股份有限公司 | 化学驱普通稠油动态降粘评价方法 |
| CN109308358A (zh) * | 2017-07-26 | 2019-02-05 | 中国石油化工股份有限公司 | 一种稠油油藏化学驱数值模拟方法 |
-
2020
- 2020-08-19 CN CN202010835220.1A patent/CN112012728B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
| US20090293594A1 (en) * | 2008-05-30 | 2009-12-03 | Lincoln Industrial Corporation | System and Method for Estimating Apparent Viscosity of a Non-Newtonian Fluid |
| JP2011069754A (ja) * | 2009-09-28 | 2011-04-07 | Mitsubishi Rayon Co Ltd | せん断粘度推算方法および流動曲線の作成方法 |
| US20120158378A1 (en) * | 2010-12-15 | 2012-06-21 | Guillaume Enchery | Method of developing a petroleum reservoir from a flow model calibrated through pressure and saturation map scaling |
| CN103775078A (zh) * | 2014-02-14 | 2014-05-07 | 武汉科技大学 | 幂律流体在多孔介质中的有效粘度的预测方法 |
| US20180023374A1 (en) * | 2015-02-03 | 2018-01-25 | Schlumberger Technology Corporation | Multi-Phase Polymer Shear Viscosity Calculation in Polymer Coreflood Simulation Study Workflow |
| CN109308358A (zh) * | 2017-07-26 | 2019-02-05 | 中国石油化工股份有限公司 | 一种稠油油藏化学驱数值模拟方法 |
| CN107869338A (zh) * | 2017-11-03 | 2018-04-03 | 中国石油化工股份有限公司 | 化学驱普通稠油动态降粘评价方法 |
Non-Patent Citations (1)
| Title |
|---|
| 张星;李兆敏;孙仁远;苏成祥;: "聚合物流变性实验研究", 新疆石油地质, no. 02, pages 197 - 199 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115600052A (zh) * | 2022-12-15 | 2023-01-13 | 昆明理工大学（Cn） | 一种确定幂律水泥浆液稠度系数的方法 |
| CN115600052B (zh) * | 2022-12-15 | 2023-09-12 | 昆明理工大学 | 一种确定幂律水泥浆液稠度系数的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112012728B (zh) | 2024-03-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Sanghani et al. | Rheology of foam and its implications in drilling and cleanout operations | |
| CN105840187A (zh) | 致密性油藏水平井分段压裂产能计算方法 | |
| CN108920849A (zh) | 一种缝洞型储层试井确定地层中洞流量方法 | |
| CN107291667B (zh) | 一种井间连通程度确定方法及系统 | |
| CN105738252A (zh) | 一种裂缝内稠油可流动开度界限的测量方法 | |
| CN107869338A (zh) | 化学驱普通稠油动态降粘评价方法 | |
| CN111062129A (zh) | 页岩油复杂缝网离散裂缝连续介质混合数值模拟方法 | |
| CN113034003B (zh) | 一种页岩气井产能快速评价方法 | |
| CN104060985A (zh) | 一种层状油藏调剖堵水堵剂进入深度测试方法及系统 | |
| CN114427445A (zh) | 缝洞型油藏无限大地层动态产能计算方法和系统 | |
| CN111502652A (zh) | 一种三孔介质气藏水平井产量递减及生产动态预测方法 | |
| CN110797092A (zh) | 一种考虑变粘酸动态过程的酸压模拟方法 | |
| Shen et al. | Spontaneous imbibition in asymmetric branch-like throat structures in unconventional reservoirs | |
| CN112012728B (zh) | 一种聚合物在储层隙间流速的确定方法 | |
| CN103743657A (zh) | 聚合物凝胶颗粒在多孔介质中通过能力的实验测定方法 | |
| CN112257349B (zh) | 一种判断致密砂岩可动水气藏气井是否具有开发价值的方法 | |
| CN113445988A (zh) | 一种低渗碳酸盐岩气藏气井产能评价方法 | |
| CN114580100B (zh) | 压裂水平井全井筒压力计算方法、设备和计算机可读储存介质 | |
| Deawwanich | Flow and displacement of viscoplastic fluids in eccentric annuli. | |
| CN105160056B (zh) | 高温高压油气直井两相流射孔完井参数与产能优化方法 | |
| Sonier et al. | A numerical model of multiphase flow around a well | |
| CN110644975B (zh) | 一种缝洞型油藏示踪剂曲线定量解释方法 | |
| CN109403942A (zh) | 一种确定页岩气井压裂设计中滑溜水与胶液混合比的方法 | |
| CN107194067A (zh) | 一种低渗透油藏数值模拟中启动压力梯度等效表征方法 | |
| Fan et al. | A generalized hydraulic calculation model for non-Newtonian fluid pipe flow and its application evaluation |
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 |