CN110797092A - 一种考虑变粘酸动态过程的酸压模拟方法 - Google Patents
一种考虑变粘酸动态过程的酸压模拟方法 Download PDFInfo
- Publication number
- CN110797092A CN110797092A CN201911060761.5A CN201911060761A CN110797092A CN 110797092 A CN110797092 A CN 110797092A CN 201911060761 A CN201911060761 A CN 201911060761A CN 110797092 A CN110797092 A CN 110797092A
- Authority
- CN
- China
- Prior art keywords
- acid
- viscosity
- variable
- concentration
- calculating
- 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
- 239000002253 acid Substances 0.000 title claims abstract description 142
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000004088 simulation Methods 0.000 title claims abstract description 22
- 230000008569 process Effects 0.000 title claims abstract description 17
- DSLZVSRJTYRBFB-UHFFFAOYSA-N Galactaric acid Natural products OC(=O)C(O)C(O)C(O)C(O)C(O)=O DSLZVSRJTYRBFB-UHFFFAOYSA-N 0.000 title claims abstract description 11
- DSLZVSRJTYRBFB-DUHBMQHGSA-N galactaric acid Chemical compound OC(=O)[C@H](O)[C@@H](O)[C@@H](O)[C@H](O)C(O)=O DSLZVSRJTYRBFB-DUHBMQHGSA-N 0.000 title claims abstract description 11
- 239000004094 surface-active agent Substances 0.000 claims abstract description 37
- 229910001424 calcium ion Inorganic materials 0.000 claims abstract description 36
- BHPQYMZQTOCNFJ-UHFFFAOYSA-N Calcium cation Chemical compound [Ca+2] BHPQYMZQTOCNFJ-UHFFFAOYSA-N 0.000 claims abstract description 35
- 238000006243 chemical reaction Methods 0.000 claims abstract description 23
- 239000011159 matrix material Substances 0.000 claims abstract description 18
- 239000007788 liquid Substances 0.000 claims abstract description 15
- 238000001914 filtration Methods 0.000 claims abstract description 12
- 230000008859 change Effects 0.000 claims abstract description 7
- 239000011435 rock Substances 0.000 claims description 27
- 238000002474 experimental method Methods 0.000 claims description 12
- 239000000243 solution Substances 0.000 claims description 11
- 239000011148 porous material Substances 0.000 claims description 7
- 238000012546 transfer Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000009792 diffusion process Methods 0.000 claims description 6
- 230000035699 permeability Effects 0.000 claims description 5
- 239000012530 fluid Substances 0.000 claims description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000002347 injection Methods 0.000 claims description 3
- 239000007924 injection Substances 0.000 claims description 3
- 238000012417 linear regression Methods 0.000 claims description 3
- 229910052739 hydrogen Inorganic materials 0.000 claims description 2
- GPRLSGONYQIRFK-UHFFFAOYSA-N hydron Chemical compound [H+] GPRLSGONYQIRFK-UHFFFAOYSA-N 0.000 claims description 2
- 238000000518 rheometry Methods 0.000 claims 1
- 238000005516 engineering process Methods 0.000 abstract description 3
- 230000009467 reduction Effects 0.000 abstract description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000020477 pH reduction Effects 0.000 description 2
- CYIDZMCFTVVTJO-UHFFFAOYSA-N pyromellitic acid Chemical compound OC(=O)C1=CC(C(O)=O)=C(C(O)=O)C=C1C(O)=O CYIDZMCFTVVTJO-UHFFFAOYSA-N 0.000 description 2
- 239000011575 calcium Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000000265 homogenisation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/10—Analysis or design of chemical reactions, syntheses or processes
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Chemical & Material Sciences (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Geochemistry & Mineralogy (AREA)
- Analytical Chemistry (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Crystallography & Structural Chemistry (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Environmental & Geological Engineering (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种考虑变粘酸动态过程的酸压模拟方法：分别测定变粘酸粘度与剪切速率、温度、钙离子浓度、pH值、表活剂浓度关系，建立变粘本构方程；模拟变粘酸酸液在裂缝系统中的流动及反应过程，得到裂缝系统内压力剖面；以裂缝系统内压力剖面为内边界条件，地面压力为外边界条件，模拟计算变粘酸在基质中的反应过程；模拟酸液在基质系统中的粘度变化情况，并计算得到滤失速度；将得到的滤失速度带入至步骤的裂缝系统内压力剖面。本发明用以解决现有技术中酸压模拟技术未考虑变粘酸的动态滤失过程，导致变粘酸酸压模拟缝长与实际情况符合率不高的问题，实现能够模拟得到考虑变粘酸动态降率后的酸压裂缝形态的目的。
Description
技术领域
本发明涉及压裂领域，具体涉及一种考虑变粘酸动态过程的酸压模拟方法。
背景技术
变粘酸是由斯伦贝谢公司在1997率先合成的一种粘度可变的酸液，用于通过水力压裂、酸化方式进行油气增产的压裂液中，在国外有大量使老油田恢复产能、新油田提高可采储量的成功运用。现有的变粘酸的初始黏度在170-1的剪切速率下可为30mpa·s左右，随着酸液与地层基质的反应，pH值升高，当pH值升至2左右时，酸液体系的粘度急剧上升，可达到1000mpa·s。之后随着酸岩反应的继续进行，pH值达到4左右时，激发另一化学反应，将残酸的粘度降低到20mpa·s或更低。酸液消耗过程中的高粘状态使得后注入的酸液转入渗透率相对较低的地层，并阻止了孔洞的发育、减缓了活性酸向天然裂缝的滤失，从而实现了均质及深度酸化的目的。而最终残酸粘度的下降，又有利于压裂液或酸液的返排。
在压裂增产领域内，酸液滤失是影响酸压缝长的主要因素之一，滤失速度越小，酸压缝长越长。变粘酸在反应过程中粘度不断发生变化，因而滤失速度也不断发生变化，最终影响酸压裂缝形态。现有的酸压模拟技术是基于常规酸液所得出的，未考虑变粘酸的动态滤失过程，导致变粘酸酸压模拟缝长与实际情况符合率不高。
发明内容
本发明的目的在于提供一种考虑变粘酸动态过程的酸压模拟方法，以解决现有技术中酸压模拟技术未考虑变粘酸的动态滤失过程，导致变粘酸酸压模拟缝长与实际情况符合率不高的问题，实现能够模拟得到考虑变粘酸动态降率后的酸压裂缝形态的目的。
本发明通过下述技术方案实现：
一种考虑变粘酸动态过程的酸压模拟方法，包括以下步骤：
(a)分别测定变粘酸粘度与剪切速率、温度、钙离子浓度、pH值、表活剂浓度关系，建立变粘本构方程；
(b)模拟变粘酸酸液在裂缝系统中的流动及反应过程，得到裂缝系统内压力剖面；
(c)以裂缝系统内压力剖面为内边界条件，地面压力为外边界条件，模拟计算变粘酸在基质中的反应过程；
(d)模拟酸液在基质系统中的粘度变化情况，并计算得到滤失速度；
(e)将得到的滤失速度带入至步骤(b)的裂缝系统内压力剖面。
进一步的，所述步骤(a)中变粘本构方程的建立方法包括：
①实验测定变粘酸在不同剪切速率下的剪切应力，通过线性回归得到变粘酸的流变指数；
②实验测定变粘酸在不同pH值下的粘度，拟合得到粘度与pH值的关系；
③实验测定变粘酸在不同钙离子浓度下的粘度，拟合得到粘度与钙离子浓度的关系；
④实验测定变粘酸在不同表面活性剂浓度下的粘度，拟合得到粘度与表面活性剂浓度的关系；
⑤建立受剪切速率、温度、钙离子浓度、pH值、表活剂浓度共同影响的变粘本构方程。
进一步的，步骤①中流变指数n的回归公式为：η＝Hγn-1，其中η为剪切应力，Pa；H为流体粘度，mPa·s；γ为剪切速率，s-1。
步骤②中拟合得到的粘度μ与pH值的关系公式为：μ(pH)＝μ0+μmax(pH)erf(b·pH+c)，其中μmax为最大黏度，μ0为初始黏度，mPa·s；pH为体系pH值，无量纲；b、c均为拟合系数，无量纲；
步骤③中拟合得到的粘度与钙离子浓度的关系公式为：
步骤④中拟合得到粘度与表面活性剂浓度的关系公式为：
步骤⑤所建立的本构方程为：
其中，α为与变粘酸性质有关的常数，无量纲；T为地层温度，K；T0为实验温度，K。
进一步的，所述步骤(b)中得到的裂缝系统内压力剖面方程为：
其中，wa为裂缝宽度，m；vx、vy分别为变粘酸在x、y方向的流速，m/s；μa为粘度，mPa·s；k为基质渗透率，Md；p为基质压力，MPa；pe为储层压力，MPa；vlm为变粘酸在裂缝壁面上不同位置的动态滤失速度，m/s。
进一步的，步骤(c)中模拟计算变粘酸在基质中的反应过程包括：计算剪切速率及流场、计算pH值及酸浓度场、计算钙离子浓度场、计算表面活性剂浓度场、计算温度场、计算酸与岩石反应后的孔隙度场、边界条件及迭代计算。
进一步的，剪切速率及流场通过流动方程计算，所述流动方程包括：
pH值及酸浓度场通过酸液浓度方程计算，所述酸液浓度方程包括：
钙离子浓度场通过钙离子平衡方程计算，所述钙离子平衡方程包括：
表面活性剂浓度场通过表面活性剂平衡方程计算，所述表面活性剂平衡方程包括：
温度场通过酸岩反应放热的热传递模型计算，所述酸岩反应放热的热传递模型包括：
其中，Tf—酸液温度，K；Ts—岩石温度，K；CPf—酸液比热容，J/kg K；CPs—岩石比热容，J/kg K；kef—酸液导热系数，W/m K；kes—岩石导热系数，W/m K；hc—酸液与岩石对流换热系数，W/m2 K；△Hr—反应焓，J/mol；ρf—酸液密度；ρs—岩石密度；
酸与岩石反应后的孔隙度场通过如下公式计算：
进一步的，所述步骤(d)中滤失速度的计算方法为：
(1)将步骤(c)中得到的剪切速率、温度、钙离子浓度、pH值、表活剂浓度带入步骤⑤所建立的本构方程中，得到基质系统中的粘度ueff分布情况；
(2)将ueff带入所述流动方程，得到不同位置的酸液流动速度，其中内边界上的流速即是裂缝壁面上的滤失速度。
进一步的，所述步骤(e)中，将得到的滤失速度带入至步骤(b)的裂缝系统内压力剖面，直至注酸循环结束。
本发明与现有技术相比，具有如下的优点和有益效果：
本发明一种考虑变粘酸动态过程的酸压模拟方法，克服了现有技术未考虑变粘酸的动态滤失过程、导致变粘酸酸压模拟缝长与实际情况符合率不高的缺陷，充分考虑了变粘酸与岩石反应的动态过程，通过模拟不同时刻下地层中的剪切速率、温度、钙离子浓度、pH值、表活剂浓度场，从而计算瞬态下的酸液粘度变化情况，并且得到酸液的动态滤失速度，最终模拟得到考虑变粘酸动态降率后的酸压裂缝形态，显著提高了变粘酸酸压模拟时的缝长准确性。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解，构成本申请的一部分，并不构成对本发明实施例的限定。在附图中：
图1为本发明具体实施例中剪切速率与剪切应力的关系图；
图2为本发明具体实施例中变粘酸在不同pH值下的粘度示意图；
图3为本发明具体实施例中变粘酸在不同钙离子浓度下的粘度示意图；
图4为本发明具体实施例中变粘酸在不同表面活性剂浓度下的粘度示意图；
图5为本发明具体实施例中变粘酸在裂缝系统内的流速剖面示意图；
图6为本发明具体实施例中变粘酸在裂缝系统内的压力剖面示意图；
图7为本发明具体实施例中裂缝系统内酸浓度场示意图；
图8为本发明具体实施例中裂缝系统内钙离子浓度场示意图；
图9为本发明具体实施例中裂缝系统内表面活性剂浓度场示意图；
图10为本发明具体实施例中变粘酸与岩石反应后的孔隙度场示意图；
图11为本发明具体实施例中基质系统中的粘度ueff分布示意图；
图12为本发明具体实施例中变粘酸在裂缝壁面上的滤失速度示意图；
图13为本发明具体实施例模拟完成后的裂缝形态图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明作进一步的详细说明，本发明的示意性实施方式及其说明仅用于解释本发明，并不作为对本发明的限定。
一种考虑变粘酸动态过程的酸压模拟方法，具体方法如下：
(一)、分别测定粘度与剪切速率、温度、钙离子浓度、pH值、表活剂浓度关系，建立变粘本构方程：
①实验测定变粘酸在不同剪切速率下的剪切应力，实验结果如图1所示，通过线性回归得到变粘酸的流变指数n：σ＝4.669·γ0.449。
②实验测定变粘酸在不同pH值下的粘度，实验结果如图2所示，通过以下公式拟合得到粘度与pH值的关系公式：
μ(pH)＝μ0+μmax(pH)erf(b·pH+c)，其中b—拟合系数，取0.609；c—拟合系数，取0.136。
③实验测定变粘酸在不同钙离子浓度下的粘度，实验结果如图3所示，通过以下公式拟合得到粘度与钙离子浓度的关系公式：
其中，d—拟合系数，取2.302；e—拟合系数，取15.203；W1——拟合系数，取18.012。
④实验测定变粘酸在不同表活剂浓度下的粘度，实验结果如图4所示，通过以下公式拟合得到粘度与表活剂浓度的关系公式：
其中，f—拟合系数，取1.000；g—拟合系数，取6.702；W2—拟合系数，取4.334。
⑤在上述实验结果的基础上，联袂变粘酸在不同环境下的粘度结果，建立受剪切速率、温度、钙离子浓度、pH值、表活剂浓度共同影响的变粘本构方程：
(二)模拟计算酸液在裂缝系统中的流动及反应过程，如图5所示，并计算得到如图6所示的裂缝系统内压力剖面。
(三)以裂缝系统内压力剖面为内边界条件，地面压力为外边界条件，计算转向酸在基质系统中的物理化学反应过程：
①通过流动方程计算剪切速率及流场，其中流动方程如下：
②通过酸液浓度方程计算pH值及酸浓度场，计算结果如图7所示，其中酸液浓度方程如下：
③钙离子平衡方程计算钙离子浓度场，计算结果如图8所示，其中钙离子平衡方程如下：
式中CCa 2+—钙离子摩尔浓度(mol/l)；De,Ca 2+—钙离子的有效扩散系数(m2/s)。
④通过表面活性剂平衡方程计算表面活性剂浓度场，计算结果如图9所示，其中表面活性剂平衡方程如下：
CSDVA—表活剂质量浓度，wt％；De,SDVA—表活剂的有效扩散系数，m2/s。
⑤考虑酸岩反应放热的热传递模型计算温度场。所述酸岩反应放热的热传递模型包括：
其中，Tf—酸液温度，K；Ts—岩石温度，K；CPf—酸液比热容，J/kg K；CPs—岩石比热容，J/kg K；kef—酸液导热系数，W/m K；kes—岩石导热系数，W/m K；hc—酸液与岩石对流换热系数，W/m2 K；△Hr—反应焓，J/mol；ρf—酸液密度；ρs—岩石密度；
⑥计算酸与岩石反应后的孔隙度场，计算结果如图10所示。
⑦边界条件(以裂缝系统内压力剖面为内边解条件，地面压力为外边界条件，等等)及迭代计算步骤，均为现有技术，在此不做赘述。
(四)、模拟酸液在基质系统中的粘度变化情况，并计算得到滤失速度：
将上述步骤(三)中得到的剪切速率、温度、钙离子浓度、pH值、表活剂浓度带入步骤(一)中第⑤小步，计算得到基质系统中的粘度分布情况ueff。
再将ueff带入步骤(三)的流动方程中，得到不同位置的酸液流动速度，其中内边界上的流速即是如图12所示的裂缝壁面上的滤失速度。
(五)将滤失速度带入第二步，循环结束，直至注酸结束，得到的模拟结束后的裂缝形态如图13所示。
以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (8)
1.一种考虑变粘酸动态过程的酸压模拟方法，其特征在于，包括以下步骤：
(a)分别测定变粘酸粘度与剪切速率、温度、钙离子浓度、pH值、表活剂浓度关系，建立变粘本构方程；
(b)模拟变粘酸酸液在裂缝系统中的流动及反应过程，得到裂缝系统内压力剖面；
(c)以裂缝系统内压力剖面为内边界条件，地面压力为外边界条件，模拟计算变粘酸在基质中的反应过程；
(d)模拟酸液在基质系统中的粘度变化情况，并计算得到滤失速度；
(e)将得到的滤失速度带入至步骤(b)的裂缝系统内压力剖面。
2.根据权利要求1所述的一种考虑变粘酸动态过程的酸压模拟方法，其特征在于，所述步骤(a)包括：
①实验测定变粘酸在不同剪切速率下的剪切应力，通过线性回归得到变粘酸的流变指数；
②实验测定变粘酸在不同pH值下的粘度，拟合得到粘度与pH值的关系；
③实验测定变粘酸在不同钙离子浓度下的粘度，拟合得到粘度与钙离子浓度的关系；
④实验测定变粘酸在不同表面活性剂浓度下的粘度，拟合得到粘度与表面活性剂浓度的关系；
⑤建立受剪切速率、温度、钙离子浓度、pH值、表活剂浓度共同影响的变粘本构方程。
3.根据权利要求2所述的一种考虑变粘酸动态过程的酸压模拟方法，其特征在于，
步骤①中流变指数n的回归公式为：η＝Hγn-1，其中η为剪切应力，Pa；H为流体粘度，mPa·s；γ为剪切速率，s-1；
步骤②中拟合得到的粘度μ与pH值的关系公式为：μ(pH)＝μ0+μmax(pH)erf(b·pH+c)，其中μmax为最大黏度，mPa·s；μ0为初始黏度，mPa·s；pH为体系pH值，无量纲；b、c均为拟合系数，无量纲；
步骤③中拟合得到的粘度与钙离子浓度的关系公式为：
步骤④中拟合得到粘度与表面活性剂浓度的关系公式为：
步骤⑤所建立的本构方程为：
其中，α为与变粘酸性质有关的常数，无量纲；T为地层温度，K；T0为实验温度，K。
4.根据权利要求3所述的一种考虑变粘酸动态过程的酸压模拟方法，其特征在于，所述步骤(b)中得到的裂缝系统内压力剖面方程为：
其中，wa为裂缝宽度，m；vx、vy分别为变粘酸在x、y方向的流速，m/s；μa为粘度，mPa·s；k为基质渗透率，Md；p为基质压力，MPa；pe为储层压力，MPa；vlm为变粘酸在裂缝壁面上不同位置的动态滤失速度，m/s。
5.根据权利要求4所述的一种考虑变粘酸动态过程的酸压模拟方法，其特征在于，步骤(c)中模拟计算变粘酸在基质中的反应过程包括：计算剪切速率及流场、计算pH值及酸浓度场、计算钙离子浓度场、计算表面活性剂浓度场、计算温度场、计算酸与岩石反应后的孔隙度场、边界条件及迭代计算。
6.根据权利要求5所述的一种考虑变粘酸动态过程的酸压模拟方法，其特征在于：
剪切速率及流场通过流动方程计算，所述流动方程包括：
pH值及酸浓度场通过酸液浓度方程计算，所述酸液浓度方程包括：
De为氢离子扩散系数，m2/s；Cs为酸液在孔隙表面的摩尔浓度mol/l；av为孔隙比表面，m2/m3；
钙离子浓度场通过钙离子平衡方程计算，所述钙离子平衡方程包括：
表面活性剂浓度场通过表面活性剂平衡方程计算，所述表面活性剂平衡方程包括：
其中，CSDVA—表活剂质量浓度，wt％；
De,SDVA—表活剂的有效扩散系数，m2/s；
温度场通过酸岩反应放热的热传递模型计算，所述酸岩反应放热的热传递模型包括：
其中，Tf—酸液温度，K；Ts—岩石温度，K；CPf—酸液比热容，J/kg K；CPs—岩石比热容，J/kg K；kef—酸液导热系数，W/m K；kes—岩石导热系数，W/m K；hc—酸液与岩石对流换热系数，W/m2 K；△Hr—反应焓，J/mol；ρf—酸液密度；ρs—岩石密度；
酸与岩石反应后的孔隙度场通过如下公式计算：
7.根据权利要求6所述的一种考虑变粘酸动态过程的酸压模拟方法，其特征在于，所述步骤(d)中滤失速度的计算方法为：
(1)将步骤(c)中得到的剪切速率、温度、钙离子浓度、pH值、表活剂浓度带入步骤⑤所建立的本构方程中，得到基质系统中的粘度ueff分布情况；
(2)将ueff带入所述流动方程，得到不同位置的酸液流动速度，其中内边界上的流速即是裂缝壁面上的滤失速度。
8.根据权利要求1所述的一种考虑变粘酸动态过程的酸压模拟方法，其特征在于，所述步骤(e)中，将得到的滤失速度带入至步骤(b)的裂缝系统内压力剖面，直至注酸循环结束。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911060761.5A CN110797092B (zh) | 2019-11-01 | 2019-11-01 | 一种考虑变粘酸动态过程的酸压模拟方法 |
| PCT/CN2019/120505 WO2021082130A1 (zh) | 2019-11-01 | 2019-11-25 | 考虑变粘酸动态过程的酸压模拟方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911060761.5A CN110797092B (zh) | 2019-11-01 | 2019-11-01 | 一种考虑变粘酸动态过程的酸压模拟方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110797092A true CN110797092A (zh) | 2020-02-14 |
| CN110797092B CN110797092B (zh) | 2022-07-19 |
Family
ID=69442496
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911060761.5A Active CN110797092B (zh) | 2019-11-01 | 2019-11-01 | 一种考虑变粘酸动态过程的酸压模拟方法 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN110797092B (zh) |
| WO (1) | WO2021082130A1 (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112179809A (zh) * | 2020-09-18 | 2021-01-05 | 清华大学 | 基于平行板流变测试得到壁面滑移曲线的方法及系统 |
| CN113671148A (zh) * | 2020-05-14 | 2021-11-19 | 中国石油天然气股份有限公司 | 一种酸蚀裂缝有效长度的实验测试方法 |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113833459A (zh) * | 2021-09-18 | 2021-12-24 | 中国科学院武汉岩土力学研究所 | 一种用于室内多段水力压裂物理模拟试验的封孔方法 |
| CN117054284B (zh) * | 2023-10-12 | 2023-12-22 | 西南石油大学 | 一种酸岩反应速率预测装置和方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102339326A (zh) * | 2010-07-16 | 2012-02-01 | 中国石油化工股份有限公司 | 一种分析模拟缝洞型油藏流体流动的方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5704386B2 (ja) * | 2010-10-12 | 2015-04-22 | 日本電気株式会社 | タンパク質分子のアミノ酸置換部位選択装置、置換アミノ酸選択装置、アミノ酸置換部位選択方法、置換アミノ酸選択方法、プログラムおよび記録媒体 |
| CN102928556B (zh) * | 2012-10-29 | 2014-05-14 | 西南石油大学 | 酸压裂缝内酸液指进大尺寸可视化物理模拟装置及方法 |
| JP6737022B2 (ja) * | 2016-07-15 | 2020-08-05 | 日立化成株式会社 | シミュレーション方法及びシミュレーション装置 |
| CN106050211A (zh) * | 2016-07-19 | 2016-10-26 | 中国石油大学(北京) | 一种缝洞型碳酸盐岩油气藏体积压裂模拟试件制备方法 |
| CN108225998B (zh) * | 2016-12-15 | 2020-10-13 | 中国石油化工股份有限公司 | 7厘米直径岩心地层温压下酸化增产模拟实验方法 |
| CN107725042B (zh) * | 2017-09-13 | 2020-06-19 | 中国石油大学(北京) | 一种高温高压大型碳酸盐岩缝洞型油藏物理模拟实验装置和方法 |
-
2019
- 2019-11-01 CN CN201911060761.5A patent/CN110797092B/zh active Active
- 2019-11-25 WO PCT/CN2019/120505 patent/WO2021082130A1/zh active Application Filing
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102339326A (zh) * | 2010-07-16 | 2012-02-01 | 中国石油化工股份有限公司 | 一种分析模拟缝洞型油藏流体流动的方法 |
Non-Patent Citations (2)
| Title |
|---|
| 薛衡: ""复杂碳酸盐岩储层水平井高效酸化数值模拟研究"", 《中国优秀博硕士学位论文全文数据库(博士) 工程科技Ⅰ辑》 * |
| 薛衡: ""考虑岩矿非均质性的前置液酸压模拟研究"", 《天然气工业》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113671148A (zh) * | 2020-05-14 | 2021-11-19 | 中国石油天然气股份有限公司 | 一种酸蚀裂缝有效长度的实验测试方法 |
| CN113671148B (zh) * | 2020-05-14 | 2024-05-28 | 中国石油天然气股份有限公司 | 一种酸蚀裂缝有效长度的实验测试方法 |
| CN112179809A (zh) * | 2020-09-18 | 2021-01-05 | 清华大学 | 基于平行板流变测试得到壁面滑移曲线的方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2021082130A1 (zh) | 2021-05-06 |
| CN110797092B (zh) | 2022-07-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110797092B (zh) | 一种考虑变粘酸动态过程的酸压模拟方法 | |
| CN102796498B (zh) | 泡沫钻井液 | |
| WO2016000090A1 (zh) | 两级封窜抑制低渗透裂缝型油藏co2驱过程中发生窜逸的采油方法 | |
| CN105699196A (zh) | 岩石渗流-应力-温度-化学耦合流变测试装置及其方法 | |
| CN109611067B (zh) | 深层石灰岩储层酸压酸液有效作用距离的数值计算方法 | |
| CN106469333B (zh) | 一种低渗稠油油藏水平井热采压力分布预测方法 | |
| CN108843303B (zh) | 一种基于泥岩蠕变模型的油水井套管损坏预测方法 | |
| CN104514531A (zh) | 一种三低砂岩油藏用双液法调剖堵水方法 | |
| CN107338033A (zh) | 一种低渗砂岩油藏复配聚合物型调剖剂及其使用方法 | |
| CN204789530U (zh) | 一种清水压裂液携砂能力测试装置 | |
| CN103539902A (zh) | 一种促进地层粘土改性的耐高温缩膨体系的制备方法 | |
| CN111353205A (zh) | 用于致密气藏产水气井地层压力和动态产能的计算方法 | |
| CN103726815B (zh) | 一种co2驱采出井井筒流态确定及参数优化方法 | |
| CN107402286A (zh) | 一种聚合物微球岩心内缓膨性能测试实验方法 | |
| CN110006788B (zh) | 在多孔介质气水界面测定堵水剂铺展性能的装置及方法 | |
| CN110298057B (zh) | 一种超临界二氧化碳压裂裂缝延伸计算方法 | |
| CN104357033A (zh) | 堵剂增效剂、含有该堵剂增效剂的凝胶及其制备方法与应用 | |
| CN107118751A (zh) | 一种适用热力采油的无机凝胶调剖剂及其使用方法 | |
| CN112980407A (zh) | 一种温度可控的凝胶堵漏剂及其制备方法和应用 | |
| CN104212418B (zh) | 高密度压井液 | |
| CN108587587B (zh) | 一种高强度可降解的油气井暂堵球及其制备方法与应用 | |
| CN107218009A (zh) | 一种低渗砂岩油藏铬铝离子复配调剖剂及其使用方法 | |
| CN104140801B (zh) | 一种铝凝胶调驱剂及其制备方法 | |
| CN107247816B (zh) | 一种稠油油藏开展冷采化学降粘最远施工半径的判定方法 | |
| CN209841613U (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 |