CN106640021A - 压后放喷参数的计算方法及装置 - Google Patents
压后放喷参数的计算方法及装置 Download PDFInfo
- Publication number
- CN106640021A CN106640021A CN201611092861.2A CN201611092861A CN106640021A CN 106640021 A CN106640021 A CN 106640021A CN 201611092861 A CN201611092861 A CN 201611092861A CN 106640021 A CN106640021 A CN 106640021A
- Authority
- CN
- China
- Prior art keywords
- pressure
- model
- parameters
- fluid
- bottom hole
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 45
- 239000012530 fluid Substances 0.000 claims abstract description 84
- 239000007788 liquid Substances 0.000 claims abstract description 51
- 238000004519 manufacturing process Methods 0.000 claims abstract description 44
- 238000000034 method Methods 0.000 claims abstract description 40
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 22
- 230000008859 change Effects 0.000 claims abstract description 20
- 238000013461 design Methods 0.000 claims abstract description 15
- 230000008569 process Effects 0.000 claims abstract description 14
- 238000001914 filtration Methods 0.000 claims abstract description 12
- 238000012360 testing method Methods 0.000 claims abstract description 9
- 239000003129 oil well Substances 0.000 claims abstract description 8
- 239000011435 rock Substances 0.000 claims description 14
- 239000010410 layer Substances 0.000 claims description 12
- 238000009533 lab test Methods 0.000 claims description 11
- 230000035699 permeability Effects 0.000 claims description 9
- 230000015556 catabolic process Effects 0.000 claims description 8
- 238000006731 degradation reaction Methods 0.000 claims description 8
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 8
- 239000011229 interlayer Substances 0.000 claims description 6
- 238000003825 pressing Methods 0.000 claims description 5
- 238000007906 compression Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 4
- 230000035515 penetration Effects 0.000 claims description 3
- 230000006835 compression Effects 0.000 claims description 2
- 208000010392 Bone Fractures Diseases 0.000 claims 14
- 206010017076 Fracture Diseases 0.000 claims 14
- 206010010214 Compression fracture Diseases 0.000 claims 2
- 238000012937 correction Methods 0.000 abstract description 4
- 239000003921 oil Substances 0.000 description 14
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000010779 crude oil Substances 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000012065 filter cake Substances 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000005094 computer simulation Methods 0.000 description 1
- 238000010924 continuous production Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000008161 low-grade oil Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Classifications
-
- 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
-
- 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
- E21B47/00—Survey of boreholes or wells
-
- 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
- E21B47/00—Survey of boreholes or wells
- E21B47/06—Measuring temperature or pressure
-
- 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)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Theoretical Computer Science (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Business, Economics & Management (AREA)
- General Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Evolutionary Computation (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Geometry (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Computer Hardware Design (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种压后放喷参数的计算方法及装置,其用于计算压裂后油井的井口压力参数以及产液量参数,压后放喷参数的计算方法包括:获取包括压裂设计数据、测井数据以及室内试验数据在内的模型参数;根据所述模型参数建立包括裂缝体积变化模型、产出液模型、压裂液滤失模型、地层供液量模型在内的多个基础模型,并计算所述各个基础模型的代表参数;根据所述各个基础模型的代表参数、体积守恒定律和压力守恒定律,建立压力计算模型,并确定所述井口压力参数;根据所述产出液模型,确定所述产液量参数。本发明的压后放喷参数的计算方法及装置考虑了地层供液对返排过程压力和产量的影响,并结合多种计算模型进行比对修正,计算结果更为准确。
Description
技术领域
本发明涉及石油开采领域,尤其涉及压裂后放喷生产领域,具体涉及一种压后放喷参数的计算方法及装置。
背景技术
油田一般采用压裂作为主要的增产措施。而随着油田勘探开发的不断深入,地层石油资源品位越来越低,主要表现为孔隙度低、渗透率低、原始地层压力低、原油物性差、储层丰度低等特点,如何衔接好压裂和压后生产成为解决这些低品位油藏的关键之一。压后放喷就是这个重要的衔接环节,精细优化压后放喷生产制度,一方面可以评价压裂效果,指导压裂工艺技术的改进;另一方面可以保证压裂效果不被破坏,保证压后高效生产。
发明内容
为了解决压裂和压后生产的衔接问题以及精细优化压后放喷生产制度,本发明提出了一种压后放喷参数的计算方法及装置。
为了达到上述目的,本发明实施例提出了一种压后放喷参数的计算方法,用于计算压裂后油井的井口压力参数以及产液量参数,包括:
获取包括压裂设计数据、测井数据以及室内试验数据在内的模型参数;
根据所述模型参数建立包括裂缝体积变化模型、产出液模型、压裂液滤失模型、地层供液量模型在内的多个基础模型,并计算所述各个基础模型的代表参数;
根据所述各个基础模型的代表参数、体积守恒定律和压力守恒定律,建立压力计算模型,并确定所述井口压力参数;
根据所述产出液模型,确定所述产液量参数。
为了达到上述目的,相应的提出了一种压后放喷参数的计算装置,用于计算压裂后油井的井口压力参数以及产液量参数,包括:
参数获取模块,用于获取包括压裂设计数据、测井数据以及室内试验数据在内的模型参数;
代表参数获取模块,用于根据所述模型参数建立包括裂缝体积变化模型、产出液模型、压裂液滤失模型、地层供液量模型在内的多个基础模型,并计算所述各个基础模型的代表参数;
井口压力参数计算模块,用于根据所述各个基础模型的代表参数、体积守恒定律和压力守恒定律,建立压力计算模型,并确定所述井口压力参数;
产液量参数计算模块,用于根据所述产出液模型,确定所述产液量参数。
本发明实施例的压后放喷参数的计算方法及装置,通过结合分段压裂和体积压裂的特点,将压裂裂缝体积模型由以往的椭圆体修正为长圆柱体,更为接近真实体积;传统计算压力-产液量-油嘴尺寸的计算结果偏小,本发明的压后放喷参数的计算方法及装置考虑了地层供液对返排过程压力和产量的影响,并结合多种计算模型进行比对修正,计算结果更为准确;并且,本发明给出的滤失量计算模型与以往的模型相比,更具可操作性,计算过程也表直观,方便。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例的压后放喷参数的计算方法的流程图。
图2为本发明实施例的裂缝体积模型侧视图。
图3为本发明实施例的裂缝体积模型俯视图。
图4为本发明实施例的压裂液滤失过程示意图。
图5为本发明实施例的压后放喷参数的计算装置的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为了对本发明的表述方便,现将本发明实施例中出现的参数进行约定,具体参数与其表达的含义如下:
Hw表示最大裂缝高度,单位为m;Lp表示裂缝的长度,单位为m;Hp表示涉及压裂底层的厚度,单位为m;Ln表示侵入厚度,单位为m;Wmax表示裂缝的最大宽度,单位为m;Pw表示井底压力,单位为MPa;Pn表示地层内压力,单位为MPa;v表示岩石泊松比,无因次;E表示岩石岩性模量,Mpa;σmin 2表示上下隔层最小主应力,单位为MPa;σmin 1表示压裂层的最小主应力,单位为MPa;a表示压裂液降解指数,无因次;n表示压裂液流变指数,无因次;μe n表示压裂液稠度系数;Ct表示综合滤失系数,单位为m/min-0.5;Z表示流体的压缩系数,无因次;θ表示局部阻力系数,无因次;r为油嘴半径,单位为m;R为油管半径,单位为m;ρw表示流体密度,单位为kg/m3;kd表示压后裂缝平均渗透率,单位m2。
油田采用压裂作为主要的增产措施,其中压后放喷生产经历如下过程:首先,压后放喷初期由于压裂裂缝内存有较高压力,远高于井口压力和地层深部压力,并形成两个压差,即生产压差△Pw和滤失压差△Pc;生产压差△Pw使压裂裂缝中的流体可以通过井筒排出地面,滤失压差△Pc使压裂裂缝中的流体向地层中滤失。然后,随着压裂裂缝内的流体不断被排出和滤失,当井底压力Pw<裂缝压力Pt<地层压力Pn时,地层停止滤失,并开始向压裂裂缝内供液,其供液体积为V供。最后,当井底压力Pw=裂缝压力Pt=地层压力Pn时,整个系统停止供液,放喷结束。
图1为本发明实施例的压后放喷参数的计算方法的流程图。如图1所示,本发明实施例的压后放喷参数的计算方法,用于计算压裂后油井的井口压力参数以及产液量参数,包括:
S100,获取包括压裂设计数据、测井数据以及室内试验数据在内的模型参数;
S200,根据所述模型参数建立包括裂缝体积变化模型、产出液模型、压裂液滤失模型、地层供液量模型在内的多个基础模型,并计算所述各个基础模型的代表参数;
S300,根据所述各个基础模型的代表参数、体积守恒定律和压力守恒定律,建立压力计算模型,并确定所述井口压力参数;
S400,根据所述产出液模型,确定所述产液量参数。
在本实施例的S100步骤中,获取包括压裂设计数据、测井数据以及室内试验数据在内的模型参数。其中,压裂设计数据包括:最大裂缝高度Hw、裂缝的长度Lp等数据参数;测井数据包括:涉及压裂底层的厚度Hp、岩石泊松比v、岩石岩性模量E、上下隔层最小主应力σmin 2、压裂层的最小主应力σmin 1等参数;室内试验数据包括:压裂液降解指数a、压裂液流变指数n、压裂液稠度系数μe n、综合滤失系数Ct、流体的压缩系数Z、局部阻力系数θ、流体密度ρw、压后裂缝平均渗透率kd等参数。以上参数的具体获取方式不限于传感器测量、实验或计算机模拟计算,本发明对此不作限定。
在本实施例的S200步骤中,根据所述模型参数建立包括裂缝体积变化模型、产出液模型、压裂液滤失模型、地层供液量模型在内的多个基础模型,并计算所述各个基础模型的代表参数;
S200步骤中的裂缝体积变化模型,其代表参数为裂缝体积变化量ΔV缝。随着放喷过程的进行,裂缝体积随压力不断变化,并同时受井井底压力Pw和地层内压力Pn大小的影响。主要表现为裂缝的最大宽度Wmax和井底压力Pw的对应关系,且
当Hw≤Hp时,
当Hw>Hp时,
其中,Pw=Pt+Ph,其中Pt为井口压力,通过井口压力表读取,Ph为井筒液柱压力吗,Ph=ρw gh,ρw为产出液密度,h为从油层中深到井口的垂直深度,单位为m。
在此说明的是,现有的椭球形模型针对直井单层单缝压裂,而对于目前常见的直井多层多缝压裂和水平井压裂,不能考虑缝间的干扰性,而采用椭圆柱体模型模拟是将多条压裂裂缝通盘考虑,在准确性上更具有优势。如图2至图4所示。则从t1时刻到t2时刻的裂缝体积变化量可根据椭圆体体积变化量计算,且:
当Hw≤Hp时,将式1带入式3中,得:
其中,βs为施工时缝内平均压力与井底压力的比值,无因次,计算式为式中n为压裂液流变指数,无因次,a为压裂液降解指数,无因次;M为中间变量,无因次,计算式为式中Wmax(t1)为t1时刻的最大井底缝宽,Wmax(t2)为t2时刻的最大井底缝宽。
当Hw>Hp时,同样将式2带入式3中,即可得裂缝体积变化量ΔV缝的表达式,在此不做赘述。
S200步骤中的产出液模型,其代表参数为产出液体积变化量ΔVout。根据伯努利基本方程,考虑油嘴修正系数,从t1时刻到t2时刻的产出液体积变化量为:
其中,Pw(t1)为t1时刻时井底压力,单位为MPa;Pw(t2)为t2时刻时井底压力,单位为MPa;为返液过程中流体连续系数,无因次,该数据为采油现场经验数据,一般连续生产无气的常规稀油取0.8~1,稠油根据井口温度和原油粘温曲线综合考虑放喷取0.3~0.6,含气油井根据含气量取0.4~0.9。
S200步骤中的压裂液滤失模型,其代表参数为滤失量ΔV滤。如图4所示,压裂裂缝内流体的滤失过程为发生滤失后,根据地层岩性的特性,首先在压裂裂缝壁面上形成滤饼区,滤饼区由许多微小颗粒和地层原岩组成,渗透性差,可以阻止缝内流体向地层内滤失,但并不完全;然后滤失继续发生,形成侵入区,直到压力平衡。在压裂裂缝长度方向上,每一点的滤失量受该点滤失压差作用,并不相同,但符合一定的数学规律,从t1时刻到t2时刻的滤失量为这些点的滤失量之和,且:
其中,Pn’(t)为滤失后的实际底层压力,单位为Mpa。
S200步骤中的地层供液量模型,其代表参数为地层供液体积ΔV供。根据瞬态径向点源扩散原理,地层任意点的注入或排出的体积,与压力和时间呈指数关系,且:
对该等式左边进行展开,可得:
进一步化简变形,得:
其中,表示压裂层平均孔隙度。
该计算地层供液体积ΔV供的方法所需要的参数较多,而大部分参数对于同一层位或区块来说,一般变化不大,比较相似,在油田勘探开发初期,地质资料不全的情况下,可以根据邻井放喷生产数据,将地层供液体积ΔV供可以简化为:
其中,a和b可通过临井生产数据拟合得到,为经验系数,是对以往复杂产能模型的简化,需要根据邻井生产数据进行回归拟合求取。
在本实施例的S300步骤中,根据所述各个基础模型的代表参数、体积守恒定律和压力守恒定律,建立压力计算模型,并确定所述井口压力参数;
该压力计算模型为
将式4、式5、式6和式10同时代入式11的体积守恒等式ΔVout=ΔV缝-ΔV供+ΔV滤中,可得:
在此,为了方便对式12的简化,令: 所以式12可以简化为:
整理得:
由此可知,
根据对上式结果取其正,则:
将式11中的压力守恒等式Pw(t)=Pt(t)+Ph(t)代入式16中,得:
将 代入至式17中,可得:
所以在给定压后初始放喷压力的情况下,通过式18即可计算出某一时间段内的井口压力。
在本实施例的S400步骤中,根据所述产出液模型,确定所述产液量参数。通过式5可计算出对应时间段内的产液量,同时根据实际产液量控制的需要,采用凑数法,可计算对应需要的油嘴尺寸。
在介绍了本发明实施例的压后放喷参数的计算方法之后,接下来,参考图5对本发明实施例的压后放喷参数的计算装置进行介绍。该装置的实施可以参见上述方法的实施,重复之处不再赘述。以下所使用的术语“模块”、“单元”,可以是实现预定功能的软件和/或硬件。
图5为本发明实施例的压后放喷参数的计算装置的结构示意图,如图5所示。本发明实施例的压后放喷参数的计算装置,用于计算压裂后油井的井口压力参数以及产液量参数,包括:
参数获取模块100,用于获取包括压裂设计数据、测井数据以及室内试验数据在内的模型参数;
代表参数获取模块200,用于根据所述模型参数建立包括裂缝体积变化模型、产出液模型、压裂液滤失模型、地层供液量模型在内的多个基础模型,并计算所述各个基础模型的代表参数;
井口压力参数计算模块300,用于根据所述各个基础模型的代表参数、体积守恒定律和压力守恒定律,建立压力计算模型,并确定所述井口压力参数;
产液量参数计算模块400,用于根据所述产出液模型,确定所述产液量参数。
在本实施例中,所述裂缝体积变化模型的代表参数为ΔV缝,且:
其中,Hw为井底最大缝高,单位为m;Lp为裂缝长度,单位m;βs为施工时缝内平均压力与井底压力的比值,无因次,计算式为式中n为压裂液流变指数,无因次,a为压裂液降解指数,无因次;M为中间变量,无因次,计算式为式中Hp为压裂地层厚度,单位为m;Wmax为最大井底缝宽,单位为m,Wmax(t1)为t1时刻的最大井底缝宽,Wmax(t2)为t2时刻的最大井底缝宽;
当Hw≤Hp时,
当Hw>Hp时,
式中,v为岩石泊松比,无因次;E为岩石岩性模量,单位为MPa;Hp为压裂地层厚度,单位为m;Pw为井底压力,单位为MPa;σmin 2为所述上下隔层最小主应力、σmin 1为所述压裂层的最小主应力,单位为MPa。
以上参数中,Hw、Lp来源于所述压裂设计数据;Hp、v、E、σmin 2、σmin 1来源于所述测井数据;n、a来源于所述室内试验数据。
在本实施例中,所述产出液模型的代表参数为ΔVout,且:
其中,Z为流体的压缩系数,无因次;为返液过程中流体连续系数,无因次;r为油嘴半径,单位为m;R为油管半径,单位为m;ρw为液体密度,单位为kg/m3;θ为局部阻力系数,无因次;Pw(t1)为t1时刻时井底压力,单位为MPa;Pw(t2)为t2时刻时井底压力,单位为MPa;
以上参数中,Z、ρw、θ来源于所述室内试验数据;r、R来源于现场实际情况。
在本实施例中,所述压裂液滤失模型的代表参数为ΔV滤,且:
其中,kd为压后裂缝平均渗透率,单位为m2;μe n为压裂液稠度系数,n为流态指数;Ln=Ct(t2-t1)0.5为侵入厚度,单位为m,Ct为综合滤失系数;Pn’(t)为t时刻滤失后的实际底层压力,单位为MPa;Pw(t)为t时刻井底压力,单位为MPa;Hp为压裂地层厚度,单位为m;Lp为裂缝长度,单位m;
以上参数中,Lp来源于所述压裂设计数据;Hp、Pn’(t)来源于所述测井数据;kd、μe n、n来源于所述室内试验数据。
在本实施例中,所述地层供液量模型的代表参数为ΔV供,且:
其中,a、b为邻井生产数据回归拟合求得的参数;Pw为井底压力,Pw(t1)为t1时刻井底压力,Pw(t2)为t2时刻井底压力,单位为MPa;
以上参数中,a来源于所述室内试验数据。
在本实施例中,所述压力计算模型为
其中,ΔVout为所述产出液模型的代表参数;ΔV缝为所述裂缝体积变化模型的代表参数;ΔV供为所述地层供液量模型的代表参数;ΔV滤为所述压裂液滤失模型的代表参数;Pw(t)为t时刻的井底压力,Pt(t)为t时刻的井口压力,Ph(t)为t时刻的液柱压力。
在本实施例中,还包括:油嘴尺寸确定模块,用于根据实际产液量控制的需要,确定油嘴尺寸参数。
在此说明的是,虽然上述提及了压后放喷参数的计算装置的若干模块,但是所属领域技术人员应该知道以上模块并非强制限定的,在实际的实施过程中多个模块的特征和功能可以在一个模块中集成实现。
本发明实施例的压后放喷参数的计算方法及装置,通过结合分段压裂和体积压裂的特点,将压裂裂缝体积模型由以往的椭圆体修正为长圆柱体,更为接近真实体积;传统计算压力-产液量-油嘴尺寸的计算结果偏小,本发明的压后放喷参数的计算方法及装置考虑了地层供液对返排过程压力和产量的影响,并结合多种计算模型进行比对修正,计算结果更为准确;并且,本发明给出的滤失量计算模型与以往的模型相比,更具可操作性,计算过程也表直观,方便。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (14)
1.一种压后放喷参数的计算方法,用于计算压裂后油井的井口压力参数以及产液量参数,其特征在于,包括:
获取包括压裂设计数据、测井数据以及室内试验数据在内的模型参数;
根据所述模型参数建立包括裂缝体积变化模型、产出液模型、压裂液滤失模型、地层供液量模型在内的多个基础模型,并计算所述各个基础模型的代表参数;
根据所述各个基础模型的代表参数、体积守恒定律和压力守恒定律,建立压力计算模型,并确定所述井口压力参数;
根据所述产出液模型,确定所述产液量参数。
2.根据权利要求1所述的压后放喷参数的计算方法,其特征在于,
所述压裂设计数据包括:井底最大缝高、裂缝长度;
所述测井数据包括:压裂地层厚度、岩石泊松比、岩石岩性模量、上下隔层最小主应力、压裂层的最小主应力;
所述室内试验数据包括:压裂液流变指数、压裂液降解指数;
所述裂缝体积变化模型的代表参数为ΔV缝,且:
其中,Hw为所述井底最大缝高,单位为m;Lp为所述裂缝长度,单位m;βs为施工时缝内平均压力与井底压力的比值,无因次,计算式为式中n为压裂液流变指数,无因次,a为压裂液降解指数,无因次;M为中间变量,无因次,计算式为式中Hp为所述压裂地层厚度,单位为m;Wmax为最大井底缝宽,单位为m,Wmax(t1)为t1时刻的最大井底缝宽,Wmax(t2)为t2时刻的最大井底缝宽;
当Hw≤Hp时,
当Hw>Hp时,
式中,v为所述岩石泊松比,无因次;E为所述岩石岩性模量,单位为MPa;Hp为所述压裂地层厚度,单位为m;Pw为井底压力,单位为MPa;σmin2为所述上下隔层最小主应力、σmin1为所述压裂层的最小主应力,单位为MPa。
3.根据权利要求1所述的压后放喷参数的计算方法,其特征在于,
所述室内试验数据包括:流体的压缩系数、液体密度、局部阻力系数;
所述产出液模型的代表参数为ΔVout,且:
其中,Z为流体的压缩系数,无因次;为返液过程中流体连续系数,无因次;r为油嘴半径,单位为m;R为油管半径,单位为m;ρw为所述液体密度,单位为kg/m3;θ为所述局部阻力系数,无因次;Pw(t1)为t1时刻时井底压力,单位为MPa;Pw(t2)为t2时刻时井底压力,单位为MPa。
4.根据权利要求1所述的压后放喷参数的计算方法,其特征在于,
所述压裂设计数据包括裂缝长度;
所述测井数据包括滤失后的实际底层压力、压裂地层厚度;
所述室内试验数据包括压后裂缝平均渗透率、压裂液稠度系数、流态指数以及综合滤失系数;
所述压裂液滤失模型的代表参数为ΔV滤,且:
其中,kd为所述压后裂缝平均渗透率,单位为m2;μe n为所述压裂液稠度系数,n为所述流态指数;Ln=Ct(t2-t1)0.5为侵入厚度,单位为m,Ct为综合滤失系数;Pn’(t)为t时刻所述滤失后的实际底层压力,单位为MPa;Pw(t)为t时刻井底压力,单位为MPa;Hp为所述压裂地层厚度,单位为m;Lp为裂缝长度,单位m。
5.根据权利要求1所述的压后放喷参数的计算方法,其特征在于,所述地层供液量模型的代表参数为ΔV供,且:
其中,a、b为邻井生产数据回归拟合求得的参数;Pw为井底压力,Pw(t1)为t1时刻井底压力,Pw(t2)为t2时刻井底压力,单位为MPa。
6.根据权利要求1所述的压后放喷参数的计算方法,其特征在于,所述压力计算模型为
其中,ΔVout为所述产出液模型的代表参数;ΔV缝为所述裂缝体积变化模型的代表参数;ΔV供为所述地层供液量模型的代表参数;ΔV滤为所述压裂液滤失模型的代表参数;Pw(t)为t时刻的井底压力,Pt(t)为t时刻的井口压力,Ph(t)为t时刻的液柱压力。
7.根据权利要求1所述的压后放喷参数的计算方法,其特征在于,在所述确定所述产液量参数之后还包括:根据实际产液量控制的需要,确定油嘴尺寸参数。
8.一种压后放喷参数的计算装置,用于计算压裂后油井的井口压力参数以及产液量参数,其特征在于,包括:
参数获取模块,用于获取包括压裂设计数据、测井数据以及室内试验数据在内的模型参数;
代表参数计算模块,用于根据所述模型参数建立包括裂缝体积变化模型、产出液模型、压裂液滤失模型、地层供液量模型在内的多个基础模型,并计算所述各个基础模型的代表参数;
井口压力参数计算模块,用于根据所述各个基础模型的代表参数、体积守恒定律和压力守恒定律,建立压力计算模型,并确定所述井口压力参数;
产液量参数计算模块,用于根据所述产出液模型,确定所述产液量参数。
9.根据权利要求8所述的压后放喷参数的计算装置,其特征在于,
所述压裂设计数据包括:井底最大缝高、裂缝长度;
所述测井数据包括:压裂地层厚度、岩石泊松比、岩石岩性模量、上下隔层最小主应力、压裂层的最小主应力;
所述室内试验数据包括:压裂液流变指数、压裂液降解指数;
所述裂缝体积变化模型的代表参数为ΔV缝,且:
其中,Hw为所述井底最大缝高,单位为m;Lp为所述裂缝长度,单位m;βs为施工时缝内平均压力与井底压力的比值,无因次,计算式为式中n为所述压裂液流变指数,无因次,a为所述压裂液降解指数,无因次;M为中间变量,无因次,计算式为式中Hp为所述压裂地层厚度,单位为m;Wmax为最大井底缝宽,单位为m,Wmax(t1)为t1时刻的最大井底缝宽,Wmax(t2)为t2时刻的最大井底缝宽;
当Hw≤Hp时,
当Hw>Hp时,
式中,v为所述岩石泊松比,无因次;E为所述岩石岩性模量,单位为MPa;Hp为所述压裂地层厚度,单位为m;Pw为井底压力,单位为MPa;σmin2为所述上下隔层最小主应力、σmin1为所述压裂层的最小主应力,单位为MPa。
10.根据权利要求8所述的压后放喷参数的计算装置,其特征在于,
所述室内试验数据包括:流体的压缩系数、液体密度、局部阻力系数;
所述产出液模型的代表参数为ΔVout,且:
其中,Z为所述流体的压缩系数,无因次;为返液过程中流体连续系数,无因次;r为油嘴半径,单位为m;R为油管半径,单位为m;ρw为所述液体密度,单位为kg/m3;θ为所述局部阻力系数,无因次;Pw(t1)为t1时刻时井底压力,单位为MPa;Pw(t2)为t2时刻时井底压力,单位为MPa。
11.根据权利要求8所述的压后放喷参数的计算装置,其特征在于,
所述压裂设计数据包括裂缝长度;
所述测井数据包括滤失后的实际底层压力、压裂地层厚度;
所述室内试验数据包括压后裂缝平均渗透率、压裂液稠度系数、流态指数以及综合滤失系数;
所述压裂液滤失模型的代表参数为ΔV滤,且:
其中,kd为所述压后裂缝平均渗透率,单位为m2;μe n为所述压裂液稠度系数,n为流态指数;Ln=Ct(t2-t1)0.5为侵入厚度,单位为m,Ct为综合滤失系数;Pn’(t)为t时刻所述滤失后的实际底层压力,单位为MPa;Pw(t)为t时刻井底压力,单位为MPa;Hp为所述压裂地层厚度,单位为m;Lp为裂缝长度,单位m。
12.根据权利要求8所述的压后放喷参数的计算装置,其特征在于,所述地层供液量模型的代表参数为ΔV供,且:
其中,a、b为邻井生产数据回归拟合求得的参数;Pw为井底压力,Pw(t1)为t1时刻井底压力,Pw(t2)为t2时刻井底压力,单位为MPa;
以上参数中,a来源于所述室内试验数据。
13.根据权利要求8所述的压后放喷参数的计算装置,其特征在于,所述压力计算模型为
其中,ΔVout为所述产出液模型的代表参数;ΔV缝为所述裂缝体积变化模型的代表参数;ΔV供为所述地层供液量模型的代表参数;ΔV滤为所述压裂液滤失模型的代表参数;Pw(t)为t时刻的井底压力,Pt(t)为t时刻的井口压力,Ph(t)为t时刻的液柱压力。
14.根据权利要求8所述的压后放喷参数的计算装置,其特征在于,还包括:油嘴尺寸确定模块,用于根据实际产液量控制的需要,确定油嘴尺寸参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611092861.2A CN106640021B (zh) | 2016-12-01 | 2016-12-01 | 压后放喷参数的计算方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611092861.2A CN106640021B (zh) | 2016-12-01 | 2016-12-01 | 压后放喷参数的计算方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106640021A true CN106640021A (zh) | 2017-05-10 |
CN106640021B CN106640021B (zh) | 2019-03-15 |
Family
ID=58814079
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611092861.2A Active CN106640021B (zh) | 2016-12-01 | 2016-12-01 | 压后放喷参数的计算方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106640021B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108166963A (zh) * | 2017-12-13 | 2018-06-15 | 中国海洋石油集团有限公司 | 一种海上油气井压裂效果评价方法 |
CN110656923A (zh) * | 2018-06-28 | 2020-01-07 | 中国石油化工股份有限公司 | 一种碳酸盐岩缝洞型油藏生产井产量预测方法 |
CN110805436A (zh) * | 2019-10-09 | 2020-02-18 | 中国石油大学(北京) | 基于压降数据的单段裂缝产液贡献率评价方法及设备 |
CN111608647A (zh) * | 2020-04-14 | 2020-09-01 | 中国石油化工股份有限公司 | 一种注水注聚井区地层压力预测方法 |
CN111963149A (zh) * | 2020-08-26 | 2020-11-20 | 中国石油天然气集团有限公司 | 一种考虑滞地液量增压的压裂后地层压力求取方法 |
CN112112609A (zh) * | 2020-09-14 | 2020-12-22 | 西南石油大学 | 一种实现气藏压后返排油嘴尺寸调控的方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103577886A (zh) * | 2012-08-06 | 2014-02-12 | 中国石油化工股份有限公司 | 一种低渗气藏水平井分段压裂产能预测方法 |
CN104727798A (zh) * | 2015-03-30 | 2015-06-24 | 中国石油集团川庆钻探工程有限公司长庆井下技术作业公司 | 一种低渗透气藏转向重复压裂工艺方法 |
CN104863560A (zh) * | 2015-03-09 | 2015-08-26 | 东方宝麟科技发展(北京)有限公司 | 用于页岩气开采的宽网压裂方法 |
CN105134159A (zh) * | 2015-07-02 | 2015-12-09 | 西南石油大学 | 一种适用于缝洞型碳酸盐岩储层的高效酸压方法 |
CN105298479A (zh) * | 2015-09-09 | 2016-02-03 | 中国石油大学(北京) | 压裂直井产油(气)位置的诊断方法及其系统 |
CN105574283A (zh) * | 2015-12-24 | 2016-05-11 | 中国石油化工股份有限公司 | 一种压裂后放喷油嘴尺寸的确定方法 |
-
2016
- 2016-12-01 CN CN201611092861.2A patent/CN106640021B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103577886A (zh) * | 2012-08-06 | 2014-02-12 | 中国石油化工股份有限公司 | 一种低渗气藏水平井分段压裂产能预测方法 |
CN104863560A (zh) * | 2015-03-09 | 2015-08-26 | 东方宝麟科技发展(北京)有限公司 | 用于页岩气开采的宽网压裂方法 |
CN104727798A (zh) * | 2015-03-30 | 2015-06-24 | 中国石油集团川庆钻探工程有限公司长庆井下技术作业公司 | 一种低渗透气藏转向重复压裂工艺方法 |
CN105134159A (zh) * | 2015-07-02 | 2015-12-09 | 西南石油大学 | 一种适用于缝洞型碳酸盐岩储层的高效酸压方法 |
CN105298479A (zh) * | 2015-09-09 | 2016-02-03 | 中国石油大学(北京) | 压裂直井产油(气)位置的诊断方法及其系统 |
CN105574283A (zh) * | 2015-12-24 | 2016-05-11 | 中国石油化工股份有限公司 | 一种压裂后放喷油嘴尺寸的确定方法 |
Non-Patent Citations (3)
Title |
---|
王才等: "压裂直井压后返排油嘴直径优选方法", 《科学技术与工程》 * |
胡景宏等: "压裂液强制返排及支撑剂回流模型研究", 《西安石油大学学报(自然科学版)》 * |
胡景宏等: "压裂液返排率的理论计算", 《钻采工艺》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108166963A (zh) * | 2017-12-13 | 2018-06-15 | 中国海洋石油集团有限公司 | 一种海上油气井压裂效果评价方法 |
CN108166963B (zh) * | 2017-12-13 | 2020-02-14 | 中国海洋石油集团有限公司 | 一种海上油气井压裂效果评价方法 |
CN110656923A (zh) * | 2018-06-28 | 2020-01-07 | 中国石油化工股份有限公司 | 一种碳酸盐岩缝洞型油藏生产井产量预测方法 |
CN110656923B (zh) * | 2018-06-28 | 2022-09-20 | 中国石油化工股份有限公司 | 一种碳酸盐岩缝洞型油藏生产井产量预测方法 |
CN110805436A (zh) * | 2019-10-09 | 2020-02-18 | 中国石油大学(北京) | 基于压降数据的单段裂缝产液贡献率评价方法及设备 |
CN111608647A (zh) * | 2020-04-14 | 2020-09-01 | 中国石油化工股份有限公司 | 一种注水注聚井区地层压力预测方法 |
CN111963149A (zh) * | 2020-08-26 | 2020-11-20 | 中国石油天然气集团有限公司 | 一种考虑滞地液量增压的压裂后地层压力求取方法 |
CN111963149B (zh) * | 2020-08-26 | 2024-03-01 | 中国石油天然气集团有限公司 | 一种考虑滞地液量增压的压裂后地层压力求取方法 |
CN112112609A (zh) * | 2020-09-14 | 2020-12-22 | 西南石油大学 | 一种实现气藏压后返排油嘴尺寸调控的方法 |
CN112112609B (zh) * | 2020-09-14 | 2022-04-05 | 西南石油大学 | 一种实现气藏压后返排油嘴尺寸调控的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106640021B (zh) | 2019-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106640021B (zh) | 压后放喷参数的计算方法及装置 | |
CN104899411B (zh) | 一种储层产能预测模型建立方法和系统 | |
CN103902827B (zh) | 一种碳酸盐岩水平井流动单元划分方法 | |
CN106150477A (zh) | 一种确定单井控制储量的方法 | |
CN108518218B (zh) | 一种非常规油气藏多段压裂水平井单井动态储量确定方法 | |
CN108979612B (zh) | 一种致密油气藏压裂酸化复杂裂缝流动能力优化方法 | |
CN103590824A (zh) | 经过多段压裂改造后的致密气藏水平井的产能计算方法 | |
CN110984973A (zh) | 缝洞型碳酸盐岩气藏单井控制储量的确定方法 | |
CN104612635A (zh) | 煤层群瓦斯联合抽采达标预评判方法 | |
Zhang et al. | An evaluation method of volume fracturing effects for vertical wells in low permeability reservoirs | |
Li et al. | Effect of perforation friction on 3D In-stage multiple fracture propagation: A numerical study | |
CN107462936A (zh) | 利用压力监测资料反演低渗透储层非达西渗流规律的方法 | |
Liu et al. | A rigorous hydraulic-fracture equilibrium-height model for multilayer formations | |
CN109441415B (zh) | 基于邻井干扰的聚合物驱油藏测试井的试井解释方法 | |
CN112541287A (zh) | 疏松砂岩压裂充填防砂增产调剖一体化设计方法 | |
CN110439519A (zh) | 一种基于极限限流设计的压裂方法及系统 | |
CN108150160A (zh) | 一种地层中欠压实作用超压的求取方法 | |
CN105443093A (zh) | 用于注聚井的井口组合测试装置及其方法 | |
CN108763751A (zh) | 一种石油地质录井中气测全烃数据的校正方法 | |
Meng et al. | Variation of induced stress field during hydraulic fracture closure and its influence on subsequent fracture | |
CN109522579B (zh) | 水平井压裂施工破裂压力预测方法 | |
CN111950112A (zh) | 一种适用于底部封闭的碳酸盐岩储层动态分析方法 | |
CN111950111B (zh) | 一种适用于底部开放的碳酸盐岩储层动态分析方法 | |
CN112177604A (zh) | 确定压裂井间干扰程度的量化评价方法 | |
CN112035993B (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 |