CN110162906B - 一种致密油藏渗流等值渗流阻力法和水电模拟*** - Google Patents
一种致密油藏渗流等值渗流阻力法和水电模拟*** Download PDFInfo
- Publication number
- CN110162906B CN110162906B CN201910457295.8A CN201910457295A CN110162906B CN 110162906 B CN110162906 B CN 110162906B CN 201910457295 A CN201910457295 A CN 201910457295A CN 110162906 B CN110162906 B CN 110162906B
- Authority
- CN
- China
- Prior art keywords
- well
- seepage
- equivalent
- pressure
- resistance
- 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
- 238000000034 method Methods 0.000 title claims abstract description 77
- 238000004088 simulation Methods 0.000 title claims description 55
- 238000004519 manufacturing process Methods 0.000 claims abstract description 69
- 230000001360 synchronised effect Effects 0.000 claims abstract description 15
- 239000012530 fluid Substances 0.000 claims abstract description 10
- 238000010586 diagram Methods 0.000 claims description 38
- 239000003921 oil Substances 0.000 claims description 38
- 238000002347 injection Methods 0.000 claims description 33
- 239000007924 injection Substances 0.000 claims description 33
- 230000035699 permeability Effects 0.000 claims description 25
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 claims description 16
- 229910052802 copper Inorganic materials 0.000 claims description 16
- 239000010949 copper Substances 0.000 claims description 16
- 239000000523 sample Substances 0.000 claims description 14
- 230000007423 decrease Effects 0.000 claims description 12
- 239000010779 crude oil Substances 0.000 claims description 7
- 230000009977 dual effect Effects 0.000 claims description 7
- 230000000694 effects Effects 0.000 claims description 6
- 230000005611 electricity Effects 0.000 claims description 6
- 239000008151 electrolyte solution Substances 0.000 claims description 6
- 230000000087 stabilizing effect Effects 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 230000015572 biosynthetic process Effects 0.000 claims description 4
- 230000035515 penetration Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 8
- 238000002474 experimental method Methods 0.000 abstract description 7
- 238000011160 research Methods 0.000 description 5
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000002608 ionic liquid Substances 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000002480 mineral oil Substances 0.000 description 2
- 235000010446 mineral oil Nutrition 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 239000011780 sodium chloride Substances 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000004401 flow injection analysis Methods 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 239000012088 reference solution Substances 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 230000005514 two-phase flow Effects 0.000 description 1
Images
Classifications
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Mining & Mineral Resources (AREA)
- Health & Medical Sciences (AREA)
- Geometry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种定产量条件下的致密油藏渗流等值渗流阻力法和水电模拟***，通过电压电流转换器，将稳压电源转换为稳流电源，用来模拟定产量生产，产量控制在较低的值，近井地带的流体缓慢采出，给远井地带流体足够的时间流向井筒，更有利于储层的动用和产量的稳定，因此实验结果更准确。本发明给出了考虑非均质性、非达西渗流、立方定律的渗流阻力公式，可更好地适用于分段压裂同步注采水平井等致密油藏开发新工艺。通过叠加原理和并行计算，研究了不同井之间的互相影响，提高了计算效率，对现场生产有较大的指导意义。通过戴维南定理、诺顿定理、函数发生器和示波器，用较简单的实验和计算步骤，就可求得不同条件下井筒的流动规律。
Description
技术领域
本发明涉及一种定产量条件下的致密油藏渗流等值渗流阻力法和水电模拟***，属于油气田开发研究模拟方法的技术领域。
背景技术
目前，油藏渗流模拟研究方法主要分为数值模拟、岩心实验和水电模拟等。多分支井、分段压裂水平井、鱼骨井、同井注采井等复杂结构井及三次采油技术如注气是致密油藏经济有效开发的必然选择，而这些手段的渗流模型、数值模拟方法较复杂，较难实现。岩心实验可采用实际油藏岩心进行，与实际地质情况、注采工艺较为接近，可直接测得较准确的生产数据，但是岩心、复杂结构井、裂缝模型等制作难度大，实验周期一般较长，费用较高，都无法满足越来越高的对致密油藏渗流模拟的要求。因此，水电模拟实验和等值渗流阻力法越来越受到人们的重视。根据文献《Generalized Dual Flow in Mesopore Medium》、《Analytical Solutions for 1-D Waterflood Problems Including CapillaryEffects》可根据水电相似原理，将电路定律与油藏渗流定律互相对偶，而研制出水电模拟实验和等值渗流阻力法，用导电介质模拟地层，在介质上施加一定的电势差产生的电场来模拟地层中的稳定渗流场是研究致密油藏稳态产能的一种简便而有效地方法，利用电路定律得到电流、电压数据，利用相似比例关系即可换算成生产数据，可适用于各种复杂井筒和裂缝模型，易实现，周期短，费用低，近年来这两种方法日益受到人们的重视。
因为稳流电源断路时容易被烧坏，不够安全，传统的水电模拟实验和等值渗流阻力法仅采用了稳压电源。因此只适用于定压生产的情况。稳压和稳流电源模型互换定律、戴维南定理、诺顿定理并没应用，因此求解模型需要求解大量方程组，对于复杂井网难以适用。
为此，本发明采用稳压和稳流电源模型互换定律对渗流等效电路图变换，可将井网不断化简，而戴维南定理、诺顿定理，不用重新求解其他井的方程，就可得到这口井不同压力下的产量或不同产量下的压力。采用叠加原理，可得到各井之间的互相影响。
发明内容
针对现有技术的不足，本发明公开一种定产量条件下的致密油藏渗流等值渗流阻力法。
本发明还公开实现上述方法的水电模拟***。
发明概述：
本发明是首次稳流电源引入到水电模拟实验中，在传统水电模拟装置的基础上，通过电压电流转换器，将稳压电源转换为稳流电源，用来模拟定产量生产，产量控制在较低的值，近井地带的流体缓慢采出，给远井地带流体足够的时间流向井筒，更有利于储层的动用和产量的稳定，因此实验结果更准确。
本发明提供一种新的揭示水平井渗流规律的电模拟***和等值渗流阻力法，以解决现有技术的如下问题：第一，现有的电模拟装置和等值渗流阻力法主要适用于定压生产的情况，但是对于致密油藏，定压生产时井筒附近流体流向井筒，远井地带难以动用，产量递减较快，结果不够精确。第二，现有的电模拟实验装置和等值渗流阻力法因为尚未采用戴维南定理、诺顿定理和示波器，要想研究井口压力半径和裂缝宽度对生产的影响需要重复实验和计算，步骤复杂。第三，现有的等值渗流阻力法并未将与油藏对偶电路化简，需要求解大量方程组。第四，现有的等值渗流阻力法并未采用叠加定理，无法研究各井对渗流场的影响。第五，现有的等值渗流阻力法并未考虑储层非均质性、非达西渗流，无法适用于致密油藏。第六，提供一种可以研究分段压裂同步注采水平井渗流规律的简便的实验和计算方法。
本发明详细的技术方案如下：
一种定产量条件下的致密油藏渗流等值渗流阻力法，包括利用水电模拟原理进行以下模拟实验，其特征在于：
在水电模拟电路中，稳压电源与电阻串联等效为稳流电源与电阻并联，对偶到油藏中，一口定压井或边界等效为两口井：一口内阻等于这个井而井口压力为0的注入井和一口定流量生产而没有内阻的井，流量等于井口压力除以井的内阻；
短路对偶为井口压力为0；
开路对偶为产量为0；
该方法的具体流程为：
1)画出油藏渗流的等效电路图，求出所述等效电流源的内阻和外阻；标出要求参数为点C；
2)判定所述等效电路图是否为最简电路图：如果是，则直接求出点C对应的参数；如果否，则固定点C，并将稳压电源化为稳流电源；
3)利用电阻、电源的串并联化简电路图至最简电路图。
根据本发明优选的，所述定产量条件下的致密油藏渗流等值渗流阻力法在应用于一个井网中时：
井组或井A单独作用表示除井A之外，将其他的定产量井开路，其他的定压井短路；
将叠加定律对偶到油藏中，其中，油藏压力场包括多个由不同的井单独作用产生的压力场；
该方法的具体流程为：
1)画出油藏渗流的等效电路图，求出所述等效电流源的内阻和外阻；
2)判断是否所有井对点C的影响都研究过了：如果是，则所有的井对点C的影响之和为点C的参数；如果否，取一口或一组井A；
3)其余定流量井开路，定压井短路；
4)利用电阻、电源的串并联化简电路图至最简电路图，重复步骤2)。
根据本发明优选的，通过考虑非均质性、非达西渗流、立方定律，用于模拟分段压裂同步注采井：
当油藏储层改造区存在分形结构时，渗透率呈线性、指数、幂律模式分布，表示为
k(x，y)＝k0fx(x)fy(y) (1)
其中，k为储层渗透率，m2；k0为近井近水力裂缝透率，m2；x为沿水力裂缝的方向，m；y为沿水平井的方向，m；D为分形维数；fx(x)、fy(y)分别表示渗透率沿x、y方向递减的规律。
渗透率以线性模式沿x、y方向递减时
k(x，y)＝k0(1-Dx)(1-Dy) (2)
渗透率以指数模式沿x、y方向递减时
k(x，y)＝k0e-Dxe-Dy (3)
渗透率以幂律模式沿x、y方向递减时
k(x，y)＝k0x-Dy-D (4)
储层改造区外阻表示为
其中，1为水力裂缝长度，m；μ为原油粘度，Pa·s；R为地层流体的渗流阻力；h为油藏高度，m。
将公式(2)(3)(4)代入公式(5)中得
未改造区原油存在低速非达西渗流，流速与压力梯度呈幂律关系时，由Ikuko-Ramey方法，未改造区的外阻为
其中，G为未改造区非达西流动系数，m为非达西流动指数。m＝1时，未改造区中的流动为达西流，G为未改造区的渗透率，lu为未改造区的长度，m。
根据本发明优选的，所述定产量条件下的致密油藏渗流等值渗流阻力法还包括：
将戴维南定理、诺顿定理对偶到油藏中，其中，从一个井网中取出任一口或一组井A，井网的其余部分对井A而言等效为：
一口定压生产井B，定压生产井B的压力等于井A开路且内阻无穷大时的压力，定压生产井B的内阻等于井A单独作用时，其余部分产生的流动阻力；
或者等效为生产井C和注入井D，等效为生产井C内阻等于井A单独作用时，其余部分产生的流动阻力，井口压力为0，注入井D没有内阻，注入量等于井A短路且内阻为0时的产量；
该方法的具体流程为：
1)画出油藏渗流的等效电路图，求出所述等效电流源的内阻和外阻；
2)假设井A之外的其余定压井短路、定流量井开路，求井网其余部分产生的阻力Rb，分别包括如下两种：
2-1)井A短路，且内阻为0时的压力pb，然后其余部分等效为井B，井B的压力等于pb,内阻等于Rb；
2-2)井A开路，且内阻无穷大时的流量Qb，然后其余部分等效为生产井C和注入井D，井C内阻等于Rb，井口压力为0，井D无内阻，注入量等于Qb；
3)将等效电路图化成最简电路图；
4)得到不同工作制度，井口尺寸下的生产数据。
一种水电模拟***，包括：
电解槽：用于盛放电解质溶液，可根据实际需要选择盐水、离子液体、矿物油；
电压电流转换器：将稳压电源转换为稳流电源，用来模拟定产量生产，其输出电流可调；
电压表：连在所述模拟井或探针和位置传感器与所述模拟裂缝之间，以获取模拟井底压力，并与数据采集器相连；
电流表：连在电压电流转换器和所述的模拟井之间，以获取产量，并与数据采集器相连；
数据采集器：用于将水电模拟***中自动记录的电压值和电流值传输给计算机。
根据本发明优选的，所述水电模拟***在用于模拟分段压裂同步注采井时，利用示波器代替电压表和电流表实现对所述电压值和电流值的数据采集；利用函数发生器作为稳压电源，该稳压电源产生锯齿波，得到不同模拟压力下的产量值或不同模拟产量下的压力值。
根据本发明优选的，在水电模拟***中，所述模拟井为分段压裂同步注采井水平井铜制模型。
根据本发明优选的，所述水电模拟***还包括：探针和位置传感器；
所述探针和位置传感器用于：利用探针在电解槽中进行模拟三维方向自由移动和模拟测量。
本发明的技术优势在于：
第一，本发明首次将稳流电源引入到水电模拟***和等值渗流阻力法，即可以考虑定产量生产，通过定产量生产，可缓解产量递减，动用远井地带流体，以适应致密油藏的来发，还可以简化渗流方程组求解步骤，使得结果更加精确。
第二，本发明将电压电流转换器将稳压电源转换为稳流电源，避免了稳流电源因为开路而被烧坏，极大地保护了实验人员的生命安全。
第三，本发明给出了考虑非均质性、非达西渗流、立方定律的渗流阻力公式，可更好地适用于分段压裂同步注采水平井等致密油藏开发新工艺。
第四，通过叠加原理和并行计算，研究了不同井之间的互相影响，提高了计算效率，对现场生产有较大的指导意义。
第五，通过戴维南定理、诺顿定理和示波器，用较简单的实验和计算步骤，就可求得不同条件下井筒的流动规律。
附图说明
图1为本发明所提供的分段压裂同步注采水平井渗流规律水电模拟***结构示意图；
图2为本发明所提供的研究裂缝尺寸、工作制度对生产参数影响的水电模拟***结构示意图；
在图1、2中，该水电模拟***包括：1、稳压电源；2、电压电流转换器；3、电流表；4、电压表；5、电解槽；6、用来模拟边界的铜片；7、用来模拟注入裂缝的铜片；8、用来模拟生产裂缝的铜片；9、探针和位置传感器；10、数据采集器；11、计算机；12、函数发生器；13、示波器；
图3为本发明所述等值渗流阻力法流程图；
图4为本发明采用叠加原理计算不同井之间的影响的流程图；
图5为本发明采用戴维南定理、诺顿定理计算裂缝尺寸、工作制度对生产参数影响的流程图；
图6为应用例1、2的油藏参数表；
图7为应用例1中五条裂缝流量，按从左到右的顺序列出；
图8为应用例2中五条裂缝流量、压力，按从左到右的顺序列出；
图9为应用例3的油藏参数；
图10为应用例3中五条裂缝压力，按从左到右的顺序列出。
具体实施方式
结合说明书附图来详细讲述本发明实施方式，但不限于此。
一种水电模拟***，包括：电解槽5：用于盛放电解质溶液，可根据实际需要选择盐水、离子液体、矿物油；
电压电流转换器2：将稳压电源转换为稳流电源，用来模拟定产量生产，其输出电流可调；
电压表4：连在电压电流转换器的两个输出端之间，以获取模拟井底压力，并与数据采集器相连；
电流表3：连在电压电流转换器和所述的模拟井之间，以获取产量，并与数据采集器相连；
数据采集器10：用于将水电模拟***中自动记录的电压值和电流值传输给计算机11。
所述水电模拟***在用于模拟分段压裂同步注采井时，利用示波器13代替电压表4和电流表3实现对所述电压值和电流值的数据采集；利用函数发生器12作为稳压电源，该稳压电源产生锯齿波，得到不同模拟压力下的产量值或不同模拟产量下的压力值。
在水电模拟***中，所述模拟井为分段压裂同步注采井水平井铜制模型，如图1、2中的用来模拟生产裂缝的铜片8。
所述水电模拟***还包括：探针和位置传感器9；
所述探针和位置传感器9用于：利用探针在电解槽5中进行模拟三维方向自由移动和模拟测量。
本发明水电模拟***工作流程包括：将所述电解槽5置于工作台上，将用来模拟边界、裂缝的铜片置于电解槽5中。电压电流转换器将稳压电源转换为稳流电源，连到用来模拟生产、注入裂缝的铜片上就可以模拟定流量注采。用来模拟边界的铜片连到稳压电源1，可用来模拟定压边界。电压表4连在模拟裂缝、边界的铜片之间，或连在模拟边界的铜片和探针之间，电流表3连在模拟裂缝的铜片和电压电流转换器2之间，调节电压电流转换器2输出电流，数据采集器10可实时采集探针位置、电流、电压，进而通过计算机11绘制压力场。
实施例1、
一种定产量条件下的致密油藏渗流等值渗流阻力法，包括利用水电模拟原理进行以下模拟实验，在水电模拟电路中，稳压电源与电阻串联等效为稳流电源与电阻并联，对偶到油藏中，一口定压井或边界等效为两口井：一口内阻等于这个井而井口压力为0的注入井和一口定流量生产而没有内阻的井，流量等于井口压力除以井的内阻；短路对偶为井口压力为0；开路对偶为产量为0；
该方法的具体流程为：
1)画出油藏渗流的等效电路图，求出所述等效电流源的内阻和外阻；标出要求参数为点C；
2)判定所述等效电路图是否为最简电路图：如果是，则直接求出点C对应的参数；如果否，则固定点C，并将稳压电源化为稳流电源；
3)利用电阻、电源的串并联化简电路图至最简电路图。
实施例2、
如实施例1所述的一种定产量条件下的致密油藏渗流等值渗流阻力法，在应用于一个井网中时：
井组或井A单独作用表示除井A之外，将其他的定产量井开路，其他的定压井短路；
将叠加定律对偶到油藏中，其中，油藏压力场包括多个由不同的井单独作用产生的压力场；
该方法的具体流程为：
1)画出油藏渗流的等效电路图，求出所述等效电流源的内阻和外阻；
2)判断是否所有井对点C的影响都研究过了：如果是，则所有的井对点C的影响之和为点C的参数；如果否，取一口或一组井A；
3)其余定流量井开路，定压井短路；
4)利用电阻、电源的串并联化简电路图至最简电路图，重复步骤2)。
实施例3、
如实施例1所述的一种定产量条件下的致密油藏渗流等值渗流阻力法，将戴维南定理、诺顿定理对偶到油藏中，其中，从一个井网中取出任一口或一组井A，井网的其余部分对井A而言等效为：
一口定压生产井B，定压生产井B的压力等于井A开路且内阻无穷大时的压力，定压生产井B的内阻等于井A单独作用时，其余部分产生的流动阻力；
或者等效为生产井C和注入井D，等效为生产井C内阻等于井A单独作用时，其余部分产生的流动阻力，井口压力为0，注入井D没有内阻，注入量等于井A短路且内阻为0时的产量；
该方法的具体流程为：
1)画出油藏渗流的等效电路图，求出所述等效电流源的内阻和外阻；
2)假设井A之外的其余定压井短路、定流量井开路，求井网其余部分产生的阻力Rb，分别包括如下两种：
2-1)井A短路，且内阻为0时的压力pb，然后其余部分等效为井B，井B的压力等于pb,内阻等于Rb；
2-2)井A开路，且内阻无穷大时的流量Qb，然后其余部分等效为生产井C和注入井D，井C内阻等于Rb，井口压力为0，井D无内阻，注入量等于Qb；
3)将等效电路图化成最简电路图；
4)得到不同工作制度，井口尺寸下的生产数据。
如图2所示，为研究定流量生产时，某裂缝不同压力下的产量或不同产量下的压力，采用函数发生器产生锯齿波，锯齿波输入到示波器X通道和电压电流转换器，得到随时间周期变化的电流，输入到模拟待研究裂缝的铜片，将该铜片或探针与模拟边界的铜片之间的电压输入到示波器Y通道，可得到该裂缝产量与压力场之间的关系。
本发明等值渗流阻力法工作流程：画出油藏渗流的等效电路图，求出内阻、外阻，标出要求参数所在的那一点，并设为点C。
当油藏储层改造区存在分形结构时，渗透率呈线性、指数、幂律模式分布，表示为
k(x，y)＝k0fx(x)fy(y) (1)
其中，k为储层渗透率，m2；k0为近井近水力裂缝透率，m2；x为沿水力裂缝的方向，m；y为沿水平井的方向，m；D为分形维数；fx(x)、fy(y)分别表示渗透率沿x、y方向递减的规律。
渗透率以线性模式沿x、y方向递减时
k(x，y)＝k0(1-Dx)(1-Dy) (2)
渗透率以指数模式沿x、y方向递减时
k(x，y)＝k0e-Dxe-Dy (3)
渗透率以幂律模式沿x、y方向递减时
k(x，y)＝k0x-Dy-D (4)
储层改造区外阻表示为
其中，1为水力裂缝长度，m；μ为原油粘度，Pa·s；R为地层流体的渗流阻力；h为油藏高度，m。
将公式(2)(3)(4)代入公式(5)中得
未改造区原油存在低速非达西渗流，流速与压力梯度呈幂律关系时，由Ikuko-Ramey方法，未改造区的外阻为
其中，G为未改造区非达西流动系数，m为非达西流动指数。m＝1时，未改造区中的流动为达西流，G为未改造区的渗透率，lu为未改造区的长度，m。
水力裂缝中流动符合立方定律，水力裂缝内阻为
其中，w为水力裂缝宽度，m；
如图3所示，固定点C，利用稳压电源与稳流电源互换定律和串并联定律，化简等效电路图。一口压力为p，内阻为R的井或边界，可等效为两口井：一口内阻为R，压力为0的井和一口流量为Q的井，反之亦然。p、Q、R之间的关系如下
图4所示的等值渗流阻力法可采用并行计算，每个计算节点对应一个井A，将除A之外的定流量井开路，定压井短路，可用图3所示的方法，求出井A对所有井流动的影响，填入一个矩阵，然后每个节点从矩阵中取出所有井对井A的影响，叠加得到井A的生产数据。
如图5所示，从井网中取井组或井A，井网的其余部分对井A而言可等效为一口压力为pb，内阻为Rb的井B，假设井A开路且内阻无穷大，可用图3图4所示的方法求出井A的压力，就是pb，假设其余部分的定流量井开路，定压井短路，其余部分产生的流动阻力就是Rb；也可等效为生产井C和注入井D，井C内阻就是Rb，井口压力为0，井D内阻为0，流量为Qb，假设井A短路且内阻为0，可用图3图4所示的方法求出井A的流量，就是Qb。改变井A的条件，可求得不同条件下井A的流动规律。
应用例1、
基于论文《Temporal scale analysis of two phase flow in fracturedwell》，专利《US 2015/0007996》和《US 2014/015603》，研究如图1所示分段压裂同步注采水平井模型，即同一水平井同时存在注入裂缝、生产裂缝，参数如图6所示。分别采用本渗流阻力法、中国专利CN201610822268.2中的实验装置得到水力裂缝产量、注入量，采用所述叠加原理、稳压电源稳流电源互换定理，得到的五条裂缝流量从左到右依次如图7所示。以实验结果作为参考解，验证了本方法的正确性。误差主要源于水力裂缝的边界效应。
应用例2、
研究数据如图6的分段压裂同步注采水平井模型，所述的生产、注入裂缝定流量生产，五条裂缝产量、注入量从左到右依次如图8所示。采用本发明所述的等值渗流阻力法和水电模拟***得到各裂缝压力，取一侧边界，另一侧边界和所有裂缝等效为一口定压井，得出井网总注入量等于总产量时，离边界最近的裂缝压力等于边界压力。采用所述稳压电源稳流电源互换定理，得到的五条裂缝压力从左到右依次如图8所示，进一步验证了本方法的正确性。误差主要源于计算的舍入误差。
应用例3、
研究如图1所示分段压裂同步注采水平井模型，所述裂缝之间区域包括改造区、未改造区，所述改造区渗透率k不仅沿x轴递减到10％，也沿y轴递减到10％，所述未改造区存在低速非达西渗流，参数如图9所示，采用所述等值渗流法求得的五条裂缝压力从左到右依次如图10所示。储层存在非达西渗流、非均质性时，生产井压力减小。
以上对本发明的原理及实施方式进行了阐述，并提供了使用案例，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。
Claims (7)
1.一种定产量条件下的致密油藏渗流等值渗流阻力法，包括利用水电模拟原理进行以下模拟实验，其特征在于：
在水电模拟电路中，稳压电源与电阻串联等效为稳流电源与电阻并联，对偶到油藏中，一口定压井或边界等效为两口井：井口压力为0的注入井，称为井A，一口定流量生产而没有内阻的井，称为井B，所述两口井可等效为定压井C，井C压力等于井B产量乘以井A的内阻；井C或定压边界等效为A、B两口井：注入井A内阻等于井C内阻；井C的井口压力为0；井B定产量生产而没有内阻，产量等于井口压力除以井的内阻；
短路对偶为井口压力为0；
开路对偶为产量为0；
该方法的具体流程为：
1)画出油藏渗流的等效电路图，求出所述等效电流源的内阻和外阻；标出要求参数为点C；所述要求参数就是所述待求产量或压力的井C；
2)判定所述等效电路图是否为最简电路图：如果是，则直接求出点C对应的参数；如果否，则固定点C，并将稳压电源化为稳流电源；
3)利用电阻、电源的串并联化简电路图至最简电路图。
2.根据权利要求1所述的一种定产量条件下的致密油藏渗流等值渗流阻力法，其特征在于，通过考虑非均质性、非达西渗流、立方定律，用于模拟分段压裂同步注采井：
当油藏储层改造区存在分形结构时，渗透率呈线性、指数、幂律模式分布，表示为
k(x，y)＝k0fx(x)fy(y) (1)
其中，k为储层渗透率，m2；k0为近井近水力裂缝透率，m2；x为沿水力裂缝的方向，m；y为沿水平井的方向，m；D为分形维数；fx(x)、fy(y)分别表示渗透率沿x、y方向递减的规律；
渗透率以线性模式沿x、y方向递减时
k(x，y)＝k0(1-Dx)(1-Dy) (2)
渗透率以指数模式沿x、y方向递减时
k(x，y)＝k0e-Dxe-Dy (3)
渗透率以幂律模式沿x、y方向递减时
k(x，y)＝k0x-Dy-D (4)
储层改造区外阻表示为
其中，1为水力裂缝长度，m；μ为原油粘度，Pa·s；R为地层流体的渗流阻力；h为油藏高度，m；
将公式(2)(3)(4)代入公式(5)中得
其中，l为水力裂缝长度，m；μ为原油粘度，Pa·s；R为地层流体的渗流阻力；h为油藏高度，m；
未改造区原油存在低速非达西渗流，流速与压力梯度呈幂律关系时，由Ikuko-Ramey方法，未改造区的外阻为
其中，G为未改造区非达西流动系数，m为非达西流动指数；m＝1时，未改造区中的流动为达西流，G为未改造区的渗透率，lu为未改造区的长度，m；
水力裂缝中流动符合立方定律，水力裂缝内阻为
其中，w为水力裂缝宽度，m。
3.根据权利要求1所述的一种定产量条件下的致密油藏渗流等值渗流阻力法，其特征在于，所述定产量条件下的致密油藏渗流等值渗流阻力法在应用于一个井网中时：
井组单独作用或井A单独作用，表示除井A之外，将其他的定产量井开路，其他的定压井短路；
将叠加定律对偶到油藏中，其中，油藏压力场包括多个由不同的井单独作用产生的压力场；
该方法的具体流程为：
1)画出油藏渗流的等效电路图，求出所述等效电流源的内阻和外阻；
2)判断是否所有井对点C的影响都研究过了：如果是，则所有的井对点C的影响之和为点C的参数；如果否，取一口或一组井A；
3)其余定流量井开路，定压井短路；
4)利用电阻、电源的串并联化简电路图至最简电路图，重复步骤2)。
4.根据权利要求1所述的一种定产量条件下的致密油藏渗流等值渗流阻力法，其特征在于，所述定产量条件下的致密油藏渗流等值渗流阻力法还包括：
将戴维南定理、诺顿定理对偶到油藏中，其中，从一个井网中取出任一口或一组井A，井网的其余部分对井A而言等效为：
一口定压生产井B，定压生产井B的压力等于井A开路且内阻无穷大时的压力，定压生产井B的内阻等于井A单独作用时，其余部分产生的流动阻力；
或者等效为生产井C和注入井D，等效的生产井C内阻等于井A单独作用时，其余部分产生的流动阻力，井口压力为0，注入井D没有内阻，注入量等于井A短路且内阻为0时的产量；
该方法的具体流程为：
1)画出油藏渗流的等效电路图，求出所述等效电流源的内阻和外阻；
2)假设井A之外的其余定压井短路、定流量井开路，求井网其余部分产生的阻力Rb，分别包括如下两种：
2-1)井A短路，且内阻为0时的压力pb，然后其余部分等效为井B，井B的压力等于pb,内阻等于Rb；
2-2)井A开路，且内阻无穷大时的流量Qb，然后其余部分等效为生产井C和注入井D，井C内阻等于Rb，井口压力为0，井D无内阻，注入量等于Qb；
3)将等效电路图化成最简电路图；
4)得到不同工作制度，井口尺寸下的生产数据。
5.一种如权利要求1所述的致密油藏渗流等值渗流阻力法的水电模拟***，其特征在于，包括：
电解槽：用于盛放电解质溶液，
电压电流转换器：将稳压电源转换为稳流电源，用来模拟定产量生产，其输出电流可调；
电压表：连在所述模拟井或探针和位置传感器与所述模拟裂缝之间，以获取模拟井底压力，并与数据采集器相连；
电流表：连在电压电流转换器和所述的模拟井之间，以获取产量，并与数据采集器相连；
数据采集器：用于将水电模拟***中自动记录的电压值和电流值传输给计算机。
6.如权利要求5所述的致密油藏渗流等值渗流阻力法的水电模拟***，其特征在于，所述水电模拟***在用于模拟分段压裂同步注采井时，利用示波器代替电压表和电流表实现对所述电压值和电流值的数据采集；利用函数发生器作为稳压电源，该稳压电源产生锯齿波，得到不同模拟压力下的产量值或不同模拟产量下的压力值。
7.如权利要求6所述的致密油藏渗流等值渗流阻力法的水电模拟***，其特征在于，在水电模拟***中，所述模拟井为分段压裂同步注采井水平井铜制模型。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910457295.8A CN110162906B (zh) | 2019-05-29 | 2019-05-29 | 一种致密油藏渗流等值渗流阻力法和水电模拟*** |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910457295.8A CN110162906B (zh) | 2019-05-29 | 2019-05-29 | 一种致密油藏渗流等值渗流阻力法和水电模拟*** |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110162906A CN110162906A (zh) | 2019-08-23 |
| CN110162906B true CN110162906B (zh) | 2020-06-23 |
Family
ID=67629655
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910457295.8A Active CN110162906B (zh) | 2019-05-29 | 2019-05-29 | 一种致密油藏渗流等值渗流阻力法和水电模拟*** |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110162906B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113919132B (zh) * | 2021-09-08 | 2024-07-09 | 中海石油(中国)有限公司海南分公司 | 一种适用于非均质潜山气藏的产能预测方法及装置 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9319005B2 (en) * | 2012-07-13 | 2016-04-19 | Rf Micro Devices, Inc. | Multi-band/multi-mode power amplifier with signal path hardware sharing |
| CA2820742A1 (en) * | 2013-07-04 | 2013-09-20 | IOR Canada Ltd. | Improved hydrocarbon recovery process exploiting multiple induced fractures |
| CN106442257A (zh) * | 2016-09-14 | 2017-02-22 | 中国石油大学(华东) | 基于离子液体的新型复杂结构井渗流电模拟*** |
| CN106401550B (zh) * | 2016-10-17 | 2022-12-23 | 长江大学 | 致密油体积压裂启动压力渗流规律电模拟装置及测试方法 |
-
2019
- 2019-05-29 CN CN201910457295.8A patent/CN110162906B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN110162906A (zh) | 2019-08-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101769151B (zh) | 基于油水前缘的电阻率数据采集处理方法 | |
| CN102930345B (zh) | 一种基于梯度算法的自适应井网优化方法 | |
| CN105089582B (zh) | 基于井下流量控制设备的油藏数值模拟方法及装置 | |
| Zhang et al. | Electricity generation and heating potential from enhanced geothermal system in Songliao Basin, China: Different reservoir stimulation strategies for tight rock and naturally fractured formations | |
| CN103104254B (zh) | 一种多功能油藏模拟实验装置及其实验方法 | |
| Aronofsky et al. | Mobility ratio-its influence on injection or production histories in five-spot water flood | |
| CN107145671B (zh) | 一种油藏数值模拟方法及*** | |
| CN206960673U (zh) | 一种探查河道溶洞的多维探测装置 | |
| Guo et al. | Parameter sensitivity analysis and optimization strategy research of enhanced geothermal system: A case study in Guide Basin, Northwestern China | |
| CN105631529A (zh) | 一种边水气藏见水时间预测方法 | |
| CN110162906B (zh) | 一种致密油藏渗流等值渗流阻力法和水电模拟*** | |
| CN106503407A (zh) | 存在部分连通断层的线性水侵油藏的试井分析方法及装置 | |
| CN109958413A (zh) | 一种特高含水期油藏动态流动单元划分方法 | |
| CN104712328A (zh) | 快速评价复杂油藏中单个流动单元动用状况的方法 | |
| CN106401550A (zh) | 致密油体积压裂启动压力梯度渗流规律电模拟装置及测试方法 | |
| CN209369793U (zh) | 一种致密油体积压裂启动压力渗流规律电模拟装置 | |
| CN102410015B (zh) | 一种复杂结构井渗流规律的电模拟***及方法 | |
| CN115587674A (zh) | 油藏改建储气库扩容达产过程气井动态产能预测方法 | |
| CN106014365B (zh) | 一种预测水驱开发油田产量递减率的方法 | |
| Wang et al. | Comprehensive evaluation of waterflooding front in low‐permeability reservoir | |
| CN103615239B (zh) | 一种测井地层成像***及其方法 | |
| CN109557280A (zh) | 一种储层应力敏感性的确定方法 | |
| CN104533401B (zh) | 一种适用于稠油油藏渗流规律的电模拟装置 | |
| CN103558109A (zh) | 一种测试电阻-含油饱和度标准关系曲线的方法 | |
| Xiao et al. | Laboratory study heterogeneity impact on microscopic residual oil distribution in tight sandstone cores during CO2 immiscible flooding |
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 |