CN106437630B - 一种基于扩散系数的分支水平井开采方法 - Google Patents
一种基于扩散系数的分支水平井开采方法 Download PDFInfo
- Publication number
- CN106437630B CN106437630B CN201610793237.9A CN201610793237A CN106437630B CN 106437630 B CN106437630 B CN 106437630B CN 201610793237 A CN201610793237 A CN 201610793237A CN 106437630 B CN106437630 B CN 106437630B
- Authority
- CN
- China
- Prior art keywords
- horizontal well
- branch horizontal
- branch
- indicate
- diffusion coefficient
- 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.)
- Expired - Fee Related
Links
- 238000009792 diffusion process Methods 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000011084 recovery Methods 0.000 title claims abstract description 24
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 68
- 238000004519 manufacturing process Methods 0.000 claims abstract description 51
- 239000011435 rock Substances 0.000 claims description 18
- 230000015572 biosynthetic process Effects 0.000 claims description 16
- 230000035699 permeability Effects 0.000 claims description 15
- 238000005553 drilling Methods 0.000 claims description 5
- 239000012530 fluid Substances 0.000 claims description 5
- 238000005065 mining Methods 0.000 claims description 5
- 230000006835 compression Effects 0.000 claims description 4
- 238000007906 compression Methods 0.000 claims description 4
- 238000000205 computational method Methods 0.000 claims description 3
- 238000005325 percolation Methods 0.000 claims description 3
- 238000012423 maintenance Methods 0.000 abstract description 6
- 230000006378 damage Effects 0.000 abstract description 5
- 230000007704 transition Effects 0.000 abstract description 5
- 230000002776 aggregation Effects 0.000 abstract description 3
- 238000004220 aggregation Methods 0.000 abstract description 3
- 239000000203 mixture Substances 0.000 abstract description 3
- 230000001186 cumulative effect Effects 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000001556 precipitation Methods 0.000 description 2
- 201000004569 Blindness Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000009467 reduction 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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开了一种基于扩散系数的分支水平井开采方法,具体按照以下步骤实施:首先计算分支水平井的页岩气产量的理论值,其次按照理论值进行分支水平井日实际产量的投产并绘制分支水平井日实际产量和时间的曲线图,当曲线图出现跳变点时对其进行调整直至曲线图不出现跳变点为止,当累计产水量和临界产水量相等时,分支水平井报废。一种基于扩散系数的分支水平井开采方法,在跃变点产生后进行调整,有利于页岩气及时携带水平井筒内的产出水,降低气水混合物返还井底的危害,减少产出水在井底的聚集,从而降低了维护费用,效率高,节省了维护成本。
Description
技术领域
本发明属于页岩气开采作业技术领域,涉及一种基于扩散系数的分支水平井开采方法。
背景技术
在页岩气开发中,在水平井基础上钻取分支井可以增加目的层长度、增大泄油面积,提高油气产量,因此,自2010年页岩气开采成本降低以后,我国分支水平井钻井技术发展迅速,并进入规模化应用阶段。目前,分支水平井技术是开采页岩气储层比较有效的手段。在页岩气井的开采中,为了获得最大的累计产量,一般地,一种方法是控制气井水平井井底的压力不变,另一种方法是控制气井的井口流量不变。
但是,现有的页岩气分支水平井开采方法通常存在不足之处:保持井底压力不变而随着时间的推移,气井的产量逐渐降低,气井开采的经济效益也随之降低;若是保持井口流量不变而随着时间的推移,地层压力差导致水平井井底产出地层水,地层水堵塞页岩气的气流通道,降低气井的开发寿命。这两种开采方法都没有考虑页岩气在储集层的扩散能力,即页岩气的分子扩散系数,这是页岩气的物理性质之一。扩散系数是沿扩散方向,在单位时间每单位浓度梯度的条件下,垂直通过单位面积所扩散某物质的质量或摩尔数,一般由实验测定。因此现有技术中亟需一种考虑页岩气扩散系数以便提高页岩气采收率的分支水平井开采方法。
发明内容
本发明的目的是提供一种基于扩散系数的分支水平井开采方法,解决了现有页岩气分支水平井开采中存在的随着开采时间的延长,水平井井底地层水堵塞页岩气的气流通道,降低气井的开发寿命的问题。
本发明所采用的技术方案是,一种基于扩散系数的分支水平井开采方法,具体按照以下步骤实施:
步骤1:确定目标储层的岩心绝对渗透率k0、岩心扩散系数Dk和页岩气的粘度μg;
步骤2:建立扩散系数影响下的气水两相渗流公式:
式中,krw表示水相相对渗透率;μw表示水相粘度;表示储集层的压力梯度;qw表示分支水平井的理论日产水量;qg表示分支水平井的日产气量;krg表示页岩气相对渗透率;Dk表示岩心扩散系数;p表示储集层任意位置的压力;
步骤3:计算分支水平井页岩气产量的理论值q0;
步骤4:根据步骤3中的q0进行分支水平井日实际产量qs的投产,即使qs=q0,同时记录生产时间,并绘制分支水平井日实际产量和时间的曲线图A,以及绘制分支水平井日实际产水量和时间的曲线图B;
步骤5:当所述步骤4的曲线图A出现跳变点时,调整分支水平井日实际产量qs,记为q1;
步骤6:根据步骤5中的q1进行分支水平井日实际产量qs的投产,即使qs=q1,并记录生产时间,并绘制分支水平井日实际产量和时间的曲线图A,以及绘制分支水平井日实际产水量和时间的曲线图B;
步骤7:当所述步骤6的曲线图A出现跳变点时,调整分支水平井日实际产量qs,记为q2;
步骤8:重复步骤6和步骤7直至分支水平井日实际产量和时间的曲线图A不出现跳变点为止,此时分支水平井日实际产量将不再进行调整;
步骤9:在分支水平井日实际产水量和时间的曲线图B上对日产水量求和得到累计产水量Q水总,当累计产水量和临界产水量相等时,分支水平井报废。
本发明的特点还在于:
步骤3中计算分支水平井页岩气产量的理论值q0的具体方法为:
按照稳定渗流理论,利用拟函数法对公式(1)进行积分得出页岩气产量公式为:
式中,l0表示气层水平长度;Zsc表示标准状态下气体压缩因子;Tsc表示标准状态下温度;ρgsc表示标准状态下气体密度;psc表示标准状态下压力;T表示地层温度;Z表示气体压缩因子;pe表示原始地层压力;pwn表示水平井分支井底压力,R(t)表示水平井分支影响下的动边界;rwn表示分支井井筒半径。
步骤5中调整分支水平井日实际产量qs的具体方法为:
步骤5.1:根据分支水平井的钻井和射孔资料确定分支数m和各相邻分支井射孔位置的间距l,并建立基于扩散系数的调整函数如下:
其中
式中,β表示页岩地层压力传播系数;v表示地层压力降传播速度;f(t)表示调整函数;t表示开采时间;INT()表示取整函数;p表示地层压力;
步骤5.2:按照步骤5.1的调整函数对分支水平井日实际产量qs进行调整,
式中,q1表示气井井口调整日产量。
步骤9中累计产水量Q水总的计算方法为:
式中,q日产水表示页岩气井日实际产水量;
步骤9中临界产水量通过目标储集层的探井和试井资料来确定。
步骤1中岩心绝对渗透率k0通过岩心渗透率自动测定仪测定。
步骤1中岩心扩散系数Dk通过实验室测定获得。
步骤1中页岩气的粘度μg通过曲线拟合法获得。
本发明的有益效果是:
1、一种基于扩散系数的分支水平井开采方法,按照分支水平井页岩气开采的特点,建立了考虑页岩气扩散系数的气水两相渗流公式,可以充分开采储集层的页岩气从而提高页岩气的采收率;
2、建立了基于扩散系数的调整函数,可以合理的调整气井的产量,避免了盲目的调整,减小了对气层的破坏;
3、建立了按照基于扩散系数的调整函数分级调整气井井口流量的公式,通过选择合理的分级调整页岩气的井口流量,可以降低地层生产压差,从而减缓地层水水窜,延长分支水平井的开采寿命;
4、按照基于扩散系数的开采方法在生产曲线图上出现跃变点时,选择合适的调整时机,可以延缓压降漏斗波及到的区域的压力降落速度,延长解吸气析出时间和析出量,从而整体上增加累计产量;在跃变点产生后进行调整,有利于页岩气及时携带水平井筒内的产出水,降低气水混合物返还井底的危害,减少产出水在井底的聚集,从而降低了维护费用,效率高,节省了维护成本。
具体实施方式
本发明提供一种基于扩散系数的分支水平井开采方法,具体包括以下步骤:
步骤1:确定目标储层的岩心绝对渗透率k0、岩心扩散系数Dk和页岩气的粘度μg,其中岩心绝对渗透率k0通过岩心渗透率自动测定仪测定,岩心扩散系数Dk通过实验室测定获得,页岩气的粘度μg通过曲线拟合法获得;
步骤2:建立扩散系数影响下的气水两相渗流公式:
式中,krw表示水相相对渗透率;μw表示水相粘度;表示储集层的压力梯度;qw表示分支水平井的理论日产水量;qg表示分支水平井的日产气量;krg表示页岩气相对渗透率;Dk表示岩心扩散系数;p表示储集层任意位置的压力;
步骤3:计算分支水平井页岩气产量的理论值q0,具体方法为:
按照稳定渗流理论,利用拟函数法对公式(1)进行积分得出页岩气产量公式为:
式中,l0表示气层水平长度;Zsc表示标准状态下气体压缩因子;Tsc表示标准状态下温度;ρgsc表示标准状态下气体密度;psc表示标准状态下压力;T表示地层温度;Z表示气体压缩因子;pe表示原始地层压力;pwn表示水平井分支井底压力;R(t)表示水平井分支影响下的动边界;rwn表示分支井井筒半径;
步骤4:根据步骤3中的q0进行分支水平井日实际产量qs的投产,即使qs=q0,并记录生产时间,并绘制分支水平井日实际产量和时间的曲线图A,以及绘制分支水平井日实际产水量和时间的曲线图B;
步骤5:当所述步骤4的曲线图A出现跳变点时,调整分支水平井日实际产量qs,记为q1;
调整分支水平井日实际产量qs的具体方法为:
步骤5.1:根据分支水平井的钻井和射孔资料确定分支数m和各相邻分支井射孔位置的间距l,并建立基于扩散系数的调整函数如下:
其中
式中,β表示页岩地层压力传播系数;v表示地层压力降传播速度;f(t)表示调整函数;t表示开采时间;INT()表示取整函数;p表示地层压力;
步骤5.2:按照步骤5.1的调整函数对分支水平井日实际产量qs进行调整,
式中,q1表示气井井口调整日产量;
步骤6:根据步骤5中的q1进行分支水平井日实际产量qs的投产,即使qs=q1,并记录生产时间,并绘制分支水平井日实际产量和时间的曲线图A,以及绘制分支水平井日实际产水量和时间的曲线图B;
步骤7:当所述步骤6的曲线图A出现跳变点时,调整分支水平井日实际产量qs,记为q2;
步骤8:重复步骤6和步骤7直至分支水平井日实际产量和时间的曲线图A不出现跳变点为止,此时分支水平井日实际产量将不再进行调整;
步骤9:在分支水平井日实际产水量和时间的曲线图B上对日产水量求和得到累计产水量Q水总,所述累计产水量Q水总的计算方法为:
式中,q日产水表示页岩气井日实际产水量;
当累计产水量和临界产水量相等时,分支水平井报废。
显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。由于页岩气井的总的产气量包含气层压力下降导致的产量和页岩气解吸气导致的产量两部分,因此通过调整页岩气的井口产量可以达到延缓地层水产出和延长气井寿命以便增加解吸气产量,从而最终达到增加页岩气井总产量的目的。
一种基于扩散系数的分支水平井开采方法,按照分支水平井页岩气开采的特点,建立了考虑页岩气扩散系数的气水两相渗流公式,可以充分开采储集层的页岩气从而提高页岩气的采收率;建立了基于扩散系数的调整函数,可以合理的调整气井的产量,避免了盲目的调整,减小了对气层的破坏;建立了按照基于扩散系数的调整函数分级调整气井井口流量的公式,通过选择合理的分级调整页岩气的井口流量,可以降低地层生产压差,从而减缓地层水水窜,延长分支水平井的开采寿命;按照基于扩散系数的开采方法在生产曲线图上出现跃变点时,选择合适的调整时机,可以延缓压降漏斗波及到的区域的压力降落速度,延长解吸气析出时间和析出量,从而整体上增加累计产量;在跃变点产生后进行调整,有利于页岩气及时携带水平井筒内的产出水,降低气水混合物返还井底的危害,减少产出水在井底的聚集,从而降低了维护费用,效率高,节省了维护成本。
Claims (8)
1.一种基于扩散系数的分支水平井开采方法,其特征在于,具体包括以下步骤:
步骤1:确定目标储层的岩心绝对渗透率k0、岩心扩散系数Dk和页岩气的粘度μg;
步骤2:建立扩散系数影响下的气水两相渗流公式:
式中,krw表示水相相对渗透率;μw表示水相粘度;表示储集层的压力梯度;qw表示分支水平井的理论日产水量;qg表示分支水平井的日产气量;krg表示页岩气相对渗透率;Dk表示岩心扩散系数;p表示储集层任意位置的压力;
步骤3:计算分支水平井页岩气产量的理论值q0;
步骤4:根据步骤3中的q0进行分支水平井日实际产量qs的投产,即使qs=q0,同时记录生产时间,并绘制分支水平井日实际产量和时间的曲线图A,以及绘制分支水平井日实际产水量和时间的曲线图B;
步骤5:当所述步骤4的曲线图A出现跳变点时,调整分支水平井日实际产量qs,记为q1;
步骤6:根据步骤5中的q1进行分支水平井日实际产量qs的投产,即使qs=q1,同时记录生产时间,并绘制分支水平井日实际产量和时间的曲线图A,以及绘制分支水平井日实际产水量和时间的曲线图B;
步骤7:当所述步骤6的曲线图A出现跳变点时,调整分支水平井日实际产量qs,记为q2;
步骤8:重复步骤6和步骤7直至分支水平井日实际产量和时间的曲线图A不出现跳变点为止,此时分支水平井日实际产量将不再进行调整;
步骤9:在分支水平井日实际产水量和时间的曲线图B上对日产水量求和得到累计产水量Q水总,当累计产水量和临界产水量相等时,分支水平井报废。
2.根据权利要求1所述的一种基于扩散系数的分支水平井开采方法,其特征在于,所述步骤1中岩心绝对渗透率k0通过岩心渗透率自动测定仪测定。
3.根据权利要求1所述的一种基于扩散系数的分支水平井开采方法,其特征在于,所述步骤1中岩心扩散系数Dk通过实验室测定获得。
4.根据权利要求1所述的一种基于扩散系数的分支水平井开采方法,其特征在于,所述步骤1中页岩气的粘度μg通过曲线拟合法获得。
5.根据权利要求1所述的一种基于扩散系数的分支水平井开采方法,其特征在于,所述步骤3中计算分支水平井页岩气产量的理论值q0的具体方法为:
按照稳定渗流理论,利用拟函数法对公式(1)进行积分得出页岩气产量公式为:
式中,l0表示气层水平长度;Zsc表示标准状态下气体压缩因子;Tsc表示标准状态下温度;ρgsc表示标准状态下气体密度;psc表示标准状态下压力;T表示地层温度;Z表示气体压缩因子;pe表示原始地层压力;pwn表示水平井分支井底压力;R(t)表示水平井分支影响下的动边界;rwn表示分支井井筒半径,m为根据分支水平井的钻井和射孔资料确定的分支数。
6.根据权利要求5所述的一种基于扩散系数的分支水平井开采方法,其特征在于,所述步骤5中调整分支水平井日实际产量qs的具体方法为:
步骤5.1:根据分支水平井的钻井和射孔资料确定分支数m和各相邻分支井射孔位置的间距l,并建立基于扩散系数的调整函数如下:
其中,式中,β表示页岩地层压力传播系数;v表示地层压力降传播速度;f(t)表示调整函数;t表示开采时间;INT()表示取整函数;
步骤5.2:按照步骤5.1的调整函数对分支水平井日实际产量qs进行调整,
式中,q1表示气井井口调整日产量。
7.根据权利要求1所述的一种基于扩散系数的分支水平井开采方法,其特征在于,所述步骤9中累计产水量Q水总的计算方法为:
式中,q日产水表示页岩气井日实际产水量,t表示开采时间。
8.根据权利要求1所述的一种基于扩散系数的分支水平井开采方法,其特征在于,所述步骤9中临界产水量通过目标储集层的探井和试井资料来确定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610793237.9A CN106437630B (zh) | 2016-10-31 | 2016-10-31 | 一种基于扩散系数的分支水平井开采方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610793237.9A CN106437630B (zh) | 2016-10-31 | 2016-10-31 | 一种基于扩散系数的分支水平井开采方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106437630A CN106437630A (zh) | 2017-02-22 |
CN106437630B true CN106437630B (zh) | 2018-08-21 |
Family
ID=58165055
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610793237.9A Expired - Fee Related CN106437630B (zh) | 2016-10-31 | 2016-10-31 | 一种基于扩散系数的分支水平井开采方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106437630B (zh) |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2783069B1 (en) * | 2011-11-22 | 2015-09-09 | Saudi Arabian Oil Company | Coupled pipe network - reservoir modeling for multi-branch oil wells |
CN103266874B (zh) * | 2013-05-17 | 2015-08-26 | 北京科技大学 | 一种开采页岩气的方法 |
CN104948163B (zh) * | 2014-03-24 | 2017-09-15 | 中国石油化工股份有限公司 | 一种页岩气井产能测定方法 |
CN104389594B (zh) * | 2014-10-13 | 2017-06-16 | 成都创源油气技术开发有限公司 | 页岩气井产能评价预测方法 |
CN105808913B (zh) * | 2014-12-31 | 2018-07-31 | 中国石油天然气股份有限公司 | 一种页岩气产量递减分析方法 |
-
2016
- 2016-10-31 CN CN201610793237.9A patent/CN106437630B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN106437630A (zh) | 2017-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106437629B (zh) | 一种基于启动压力梯度的分支水平井开采方法 | |
CN106651610B (zh) | 一种浅层超低渗砂岩油藏注水开发动态分析方法 | |
CN103527163B (zh) | 一种致密储层水平井体积压裂工艺 | |
CN105626036B (zh) | 一种确定油藏合理产液量油藏工程计算方法 | |
CN106484933A (zh) | 一种用于确定页岩气井井控动态储量的方法及系统 | |
CN107630686B (zh) | 水平井段间交错驱替和渗吸驱油的致密油能量补充方法 | |
CN108331574B (zh) | 一种水平井水平段相对出砂亏空剖面预测及防砂分段分级方法 | |
CN107705215A (zh) | 一种页岩储层重复压裂选井选段方法 | |
CN107832540A (zh) | 一种致密油藏技术极限井距确定方法 | |
CN103605874B (zh) | 无量纲系数压裂优化设计方法 | |
CN104234677A (zh) | 一种注气垂直驱替提高凝析气藏凝析油采收率方法 | |
CN110259421A (zh) | 一种裂缝性的致密油藏注水补充能量方法 | |
CN106651158B (zh) | 一种超低渗致密油藏水平井注水开发见效程度定量评价方法 | |
CN108518218A (zh) | 一种非常规油气藏多段压裂水平井单井动态储量确定方法 | |
CN106014372A (zh) | 一种基于砂体结构的水平井布井方法 | |
CN108518219A (zh) | 一种基于天然能量开发油田规模提液选井的评价方法 | |
CN111027789B (zh) | 一种煤层气井工作制度定量优化的计算方法 | |
CN104389566A (zh) | 一种判定气体窜逸时间的方法 | |
CN102364041B (zh) | 水平井压裂缝填充孚盛砂建立透油阻水筛的采油方法 | |
CN106285568B (zh) | 一种基于格林函数的分支水平井开采方法 | |
CN104100237A (zh) | 一种煤层气井的排采方法 | |
CN106437630B (zh) | 一种基于扩散系数的分支水平井开采方法 | |
CN106351612B (zh) | 一种基于分形维数的分支水平井开采方法 | |
CN106285569B (zh) | 一种基于滑脱系数的分支水平井开采方法 | |
Denbina et al. | Modelling cold production for heavy oil reservoirs |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180821 Termination date: 20191031 |