CN104863547A - 一种基于压力传播的煤层气开发初期排采速度确定方法 - Google Patents
一种基于压力传播的煤层气开发初期排采速度确定方法 Download PDFInfo
- Publication number
- CN104863547A CN104863547A CN201510174027.7A CN201510174027A CN104863547A CN 104863547 A CN104863547 A CN 104863547A CN 201510174027 A CN201510174027 A CN 201510174027A CN 104863547 A CN104863547 A CN 104863547A
- Authority
- CN
- China
- Prior art keywords
- pressure
- well
- inter
- propagation
- initial stage
- 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
Landscapes
- Measuring Fluid Pressure (AREA)
Abstract
本发明涉及一种基于压力传播的煤层气开发初期排采速度确定方法,包括以下步骤:1)根据压力传播距离与时间关系,确定给定井网井距条件下井间压力发生干扰所需时间;2)根据步骤1)得到的井间压力发生干扰所需时间,按照均衡降压的理念,确定降压速率;3)根据步骤2)得到的降压速率,确定井筒液面下降速率。本发明具有严格的理论基础,实现了煤层气开发初期排采速度的定量确定。本发明可以广泛应用于不同储层物性条件煤层气田的煤层气开发初期排采速度的定量确定。
Description
技术领域
本发明涉及一种排采速度确定方法,特别是关于一种基于压力传播的煤层气开发初期排采速度确定方法。
背景技术
目前煤层气主要通过排水降压采气进行开发,煤层一般都是饱和水的,要使煤层中吸附气解吸产出,开发初期需进行排水降压。排水有快慢的问题,从缩短开发周期角度来说,初期排采速度越快,气井见气越快,但排采速度过快,一方面会造成煤层激动,使得煤粉产出,堵塞渗流通道;另一方面,排采速度过快会使得煤层气井过早出现气液两相流,使得压降漏斗小,泄气面积小。目前煤层气排采制度的确定主要考虑控制煤粉产出和减小应力敏感效应,排采速度均来自现场经验积累,缺乏理论基础,在不同煤层气田适用性差。对煤层气开发过程中相态变化对压降漏斗及井间压力干扰的影响方面,目前还停留在定性认识上,缺乏定量确定方法。
发明内容
针对上述问题,本发明的目的是提供一种可以定量确定煤层气开发初期排采速度且适用范围广的基于压力传播的煤层气开发初期排采速度确定方法。
为实现上述目的,本发明采取以下技术方案:一种基于压力传播的煤层气开发初期排采速度确定方法,包括以下步骤:1)根据压力传播距离与时间的关系,确定给定井网井距条件下井间压力发生干扰所需时间;2)根据步骤1)得到的井间压力发生干扰所需时间,按照均衡降压的理念,确定降压速率:
△p/△t≤[pi-pd]/tint (1)
其中:△p/△t为降压速率,MPa/d;pi为原始煤层压力,MPa;pd为临界解吸压力,MPa;tint为井间压力发生干扰所需时间,d;3)根据步骤2)得到的降压速率,确定煤层气开发初期井筒液面降低速率:
其中:△h/△t为井筒液面降低速率,m/d;ρw为水的密度,kg/m3;g为重力加速度,等于9.8m/s2。
所述步骤1)中井间压力发生干扰所需时间tint依据直井单相水压力传播距离公式确定;对于直井未压裂情况,压力传播在平面上以圆形向外扩展,压力传播半径满足:
其中:rinv为圆形传播半径,m;t为生产时间,d;k为割理渗透率,mD;φ为煤层割理孔隙度;μ为地层水粘度,mPa·s;Ct为割理系统综合压缩系数,1/MPa;当rinv等于井距一半时,求得井间压力发生干扰所需时间tint;对于压裂直井,压力传播在平面上以椭圆形态向外扩展,椭圆的长短半轴分别满足:
其中:ai、bi分别为椭圆的长半轴和短半轴,m;ξi为中间参数,m;Lf为压裂直井裂缝半长,m;根据实际井距分别求得椭圆长、短半轴方向的井间压力干扰时间,取二者间的最大值为最终井间压力发生干扰所需时间tint。
本发明由于采取以上技术方案,其具有以下优点:1、本发明由于依据压力传播模型确定出给定井距条件下井间压力发生干扰所需时间,并由该时间确定出合理液面下降速率,因此实现了定量确定煤层气开发初期排采速度。2、本发明由于是基于不同开发阶段压力传播规律而建立起来的,因此具有理论基础。3、本发明由于建立了合理的液面下降速率模型,因此适合不同储层物性条件时排采速度的确定。本发明可以广泛应用于不同储层物性条件煤层气田的煤层气开发初期排采速度的定量确定。
附图说明
图1是煤层气初期定排水而后定压生产时压降漏斗曲线示意图
图2是初期排采速度对压力传播距离影响示意图
图3为本发明的流程图
图4为井间压力干扰示意图
图5为井间干扰时间确定示意图
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
在描述本发明的技术方案之前,首先简单说明一下本发明基于的原理,主要包括以下两点:
第一,煤层气开发过程气水两相压力传播速度显著小于初期单相水阶段:由图1可以看出,初期定排水时,压降速度快;随着压力降到解吸压力,形成气液两相流,压降主要消耗在气水两相区,两相区向外扩展速度很慢。
第二,排采速度过快会导致气井过早见气,一旦见气,压降主要消耗在两相区,压力传播速度慢,导致地层远处割理中的水无法排出,井间干扰慢:由图2可以看出,初期采用小排水量生产时,压降漏斗较平缓且延伸远;大排水量生产时,气井见气过早,气液两相渗流阻力大导致压力传播速度显著降低,压降漏斗曲线陡且延伸距离短。为了使煤层气压降漏斗传播的尽可能远,使得井与井之间达到压力干扰,要求在两井压力干扰之前煤层保持单相水流状态。换言之,就是保证井间压力干扰前井底压力尚未达到临界解吸压力。通过确定井间压力干扰所需时间,可制定确定排采速度的方案。
基于以上原理,本发明提出了一种基于压力传播的煤层气开发初期排采速度确定方法,如图3所示,其包括以下步骤:
1)收集煤层物性及流体资料,包括煤层割理孔隙度、地层水粘度和割理系统综合压缩系数。
2)收集井参数,包括井距、压裂井的裂缝半长和方向。
3)根据压力传播距离与时间的关系,确定给定井网井距条件下井间压力发生干扰所需时间tint。由于该时间tint内煤层气保持单相水流状态,因此临界状态是压力开始干扰时刻的井底压力等于临界解吸压力。
井间压力发生干扰所需时间tint依据直井单相水压力传播距离公式确定:
①对于直井未压裂情况,压力传播在平面上以圆形向外扩展,压力传播半径满足:
其中:rinv为圆形传播半径,m;t为生产时间,d;k为煤层割理渗透率,mD;φ为煤层割理孔隙度;μ为地层水粘度,mPa·s;Ct为割理系统综合压缩系数,1/MPa。当rinv等于井距一半时,求得井间压力发生干扰所需时间tint。
②对于压裂直井,压力传播在平面上以椭圆形态向外扩展,椭圆的长短半轴分别满足:
其中:ai、bi分别为椭圆的长半轴和短半轴,m;ξi为中间参数,m;Lf为压裂直井裂缝半长,m。根据实际井距可分别求得椭圆长、短半轴方向的井间压力干扰时间,取二者间的最大值为最终井间压力发生干扰所需时间tint,以保证两方向井间压力均达到干扰。
4)收集原始地层压力和临界解吸压力值,并根据得到的井间压力发生干扰所需时间tint,按照均衡降压的理念,合理降压速率应控制在:
△p/△t≤[pi-pd]/tint (4)
其中:△p/△t为降压速率,MPa/d;pi为原始煤层压力,MPa;pd为临界解吸压力,MPa。
5)降压速率是通过排水降低液面高度来实现的,因此,煤层气开发初期井筒液面降低速率(煤层气开发初期排采速度可通过煤层气开发初期井筒液面降低速率来表达)为:
其中:△h/△t为井筒液面降低速率,m/d;ρw为水的密度,kg/m3;g为重力加速度,等于9.8m/s2。
下面通过一个实施例来具体说明本发明的工作流程:
1)收集煤层物性及流体资料,包括煤层割理孔隙度、地层水粘度和割理系统综合压缩系数。举例:k=1~10mD,φ=0.02,μ=1mpa·s;Ct=0.0367MPa-1。
2)收集井参数,包括井距、压裂井的裂缝半长和方向。以裂缝井为例,井距为400m×300m,裂缝半长Lf=60m,裂缝方向沿着大井距方向,井点中心到两井距离分别为ai=200m和bi=150m。
3)计算井间压力发生干扰所需时间,该时间内煤层气保持单相水流状态,临界状态是压力开始干扰时刻的井底压力等于临界解吸压力,如图4所示。对于无压裂直井,采用公式(1)求取井间压力发生干扰所需时间;对于压裂直井,采用公式(2)和(3)求取椭圆长半轴方向和短半轴方向压力传播与时间的关系。以裂缝井为例,计算的压力传播距离与时间关系如图5所示,其中(a)为椭圆长半轴方向压力传播距离与时间关系示意图,(b)为椭圆短半轴方向压力传播距离与时间关系示意图,根据ai、bi可确定不同渗透率情况下井间压力干扰时间tint1、tint2、tint3和tint4,当k=2mD时,椭圆长、短半轴方向干扰时间分别为38d和23d,取最大干扰时间38d为最终井间压力发生干扰所需时间。
4)收集原始煤层压力和临界解吸压力值,根据公式(4)可计算合理降压速率,根据公式(5)计算煤层气开发初期液面降低速率。以裂缝井为例,假定原始煤层压力为3MPa,临界解吸压力为1.5MPa,则
上述各实施例仅用于说明本发明,其中的一些方法和步骤等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。
Claims (2)
1.一种基于压力传播的煤层气开发初期排采速度确定方法,包括以下步骤:
1)根据压力传播距离与时间的关系,确定给定井网井距条件下井间压力发生干扰所需时间;
2)根据步骤1)得到的井间压力发生干扰所需时间,按照均衡降压的理念,确定降压速率:
△p/△t≤[pi-pd]/tint (1)
其中:△p/△t为降压速率,MPa/d;pi为原始煤层压力,MPa;pd为临界解吸压力,MPa;tint为井间压力发生干扰所需时间,d;
3)根据步骤2)得到的降压速率,确定煤层气开发初期井筒液面降低速率:
其中:△h/△t为井筒液面降低速率,m/d;ρw为水的密度,kg/m3;g为重力加速度,等于9.8m/s2。
2.如权利要求1所述的一种基于压力传播的煤层气开发初期排采速度确定方法,其特征在于:所述步骤1)中井间压力发生干扰所需时间tint依据直井单相水压力传播距离公式确定;对于直井未压裂情况,压力传播在平面上以圆形向外扩展,压力传播半径满足:
其中:rinv为圆形传播半径,m;t为生产时间,d;k为割理渗透率,mD;φ为煤层割理孔隙度;μ为地层水粘度,mPa·s;Ct为割理系统综合压缩系数,1/MPa;当rinv等于井距一半时,求得井间压力发生干扰所需时间tint;
对于压裂直井,压力传播在平面上以椭圆形态向外扩展,椭圆的长短半轴分别满足:
其中:ai、bi分别为椭圆的长半轴和短半轴,m;ξi为中间参数,m;Lf为压裂直井裂缝半长,m;根据实际井距分别求得椭圆长、短半轴方向的井间压力干扰时间,取二者间的最大值为最终井间压力发生干扰所需时间tint。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510174027.7A CN104863547B (zh) | 2015-04-13 | 2015-04-13 | 一种基于压力传播的煤层气开发初期排采速度确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510174027.7A CN104863547B (zh) | 2015-04-13 | 2015-04-13 | 一种基于压力传播的煤层气开发初期排采速度确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104863547A true CN104863547A (zh) | 2015-08-26 |
CN104863547B CN104863547B (zh) | 2017-08-11 |
Family
ID=53909670
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510174027.7A Active CN104863547B (zh) | 2015-04-13 | 2015-04-13 | 一种基于压力传播的煤层气开发初期排采速度确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104863547B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107816335A (zh) * | 2016-09-13 | 2018-03-20 | 中国石油天然气股份有限公司 | 一种煤层气含气量开发下限的确定方法及装置 |
CN108468540A (zh) * | 2018-03-29 | 2018-08-31 | 西南石油大学 | 一种低渗-特低渗油气藏单井泄流面积的计算方法 |
CN109488284A (zh) * | 2018-11-23 | 2019-03-19 | 中海油能源发展股份有限公司 | 煤层气直井排采过程水压传播距离模拟测试装置 |
CN110413951A (zh) * | 2018-04-28 | 2019-11-05 | 中国石油天然气股份有限公司 | 一种煤层气井排采速度的确定方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102493788A (zh) * | 2011-12-15 | 2012-06-13 | 中国石油天然气股份有限公司 | 高煤阶煤层气井排水采气的方法 |
-
2015
- 2015-04-13 CN CN201510174027.7A patent/CN104863547B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102493788A (zh) * | 2011-12-15 | 2012-06-13 | 中国石油天然气股份有限公司 | 高煤阶煤层气井排水采气的方法 |
Non-Patent Citations (2)
Title |
---|
李清等: "延川南工区煤层气排采速率定量分析", 《石油天然气学报》 * |
段品佳等: "煤层气排采初期阶段合理降压速率的研究", 《煤炭学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107816335A (zh) * | 2016-09-13 | 2018-03-20 | 中国石油天然气股份有限公司 | 一种煤层气含气量开发下限的确定方法及装置 |
CN107816335B (zh) * | 2016-09-13 | 2019-12-10 | 中国石油天然气股份有限公司 | 一种煤层气含气量开发下限的确定方法及装置 |
CN108468540A (zh) * | 2018-03-29 | 2018-08-31 | 西南石油大学 | 一种低渗-特低渗油气藏单井泄流面积的计算方法 |
CN108468540B (zh) * | 2018-03-29 | 2021-04-09 | 西南石油大学 | 一种低渗-特低渗油气藏单井泄流面积的计算方法 |
CN110413951A (zh) * | 2018-04-28 | 2019-11-05 | 中国石油天然气股份有限公司 | 一种煤层气井排采速度的确定方法 |
CN110413951B (zh) * | 2018-04-28 | 2023-04-25 | 中国石油天然气股份有限公司 | 一种煤层气井排采速度的确定方法 |
CN109488284A (zh) * | 2018-11-23 | 2019-03-19 | 中海油能源发展股份有限公司 | 煤层气直井排采过程水压传播距离模拟测试装置 |
Also Published As
Publication number | Publication date |
---|---|
CN104863547B (zh) | 2017-08-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105587298B (zh) | 多流管模式的水驱油井含水率反演方法 | |
CN109209333B (zh) | 页岩气多井群高效开采间距优化方法 | |
CN104863547A (zh) | 一种基于压力传播的煤层气开发初期排采速度确定方法 | |
CN103939066B (zh) | 一种一注多采井组定注水量确定油井产液量的方法 | |
CN104834807B (zh) | 一种基于分形理论的应力敏感储层相对渗透率计算方法 | |
CN110359899B (zh) | 一种页岩气水平井重复压裂提高有效改造体积的方法 | |
CN105569613B (zh) | 一种中高阶煤煤层气排采方法 | |
CN106545336A (zh) | 考虑致密气藏渗流机理的产能计算方法 | |
CN107630686B (zh) | 水平井段间交错驱替和渗吸驱油的致密油能量补充方法 | |
CN105952432A (zh) | 一种超低渗致密油藏体积压裂水平井准自然能量开发交错布井方法 | |
CN104975827B (zh) | 预测二氧化碳驱油藏指标的物质平衡方法 | |
CN111749658B (zh) | 二氧化碳吞吐采油方法及装置 | |
CN103643928B (zh) | 一种基于压力场和流速场分布的逐级深部调剖方法 | |
CN109933845B (zh) | 一种暂堵压裂时天然裂缝的开启范围计算方法 | |
CN105625990A (zh) | 一种体积压裂水平井开采页岩油的方法 | |
CN103334725A (zh) | 评价低渗透油藏驱替有效性的方法及装置 | |
Christman et al. | Comparison of laboratory-and field-observed CO2 tertiary injectivity | |
CN104594854A (zh) | 底水油藏水平井变参数射孔控水完井方法及装置 | |
CN111963163B (zh) | 一种基于气水两相流动页岩气藏双孔双渗模型构建及压力动态预测方法 | |
CN106958437A (zh) | 一种油井压裂提高采收率新方法 | |
CN109918769A (zh) | 利用瞬时方程计算缝洞型油藏非稳态水侵水侵量的方法 | |
CN108830410A (zh) | 低渗透底水气藏的见水时间预测方法及装置 | |
CN211448630U (zh) | 一种降压和双管注入改性流体开采天然气水合物的装置 | |
CN104573215A (zh) | 一种调剖剂颗粒浓度及直径的计算方法 | |
CN113486537B (zh) | 一种基于气水交替驱油的井距设计方法、装置及设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP01 | Change in the name or title of a patent holder |
Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Co-patentee after: CNOOC research institute limited liability company Patentee after: China Offshore Oil Group Co., Ltd. Address before: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Co-patentee before: CNOOC Research Institute Patentee before: China National Offshore Oil Corporation |
|
CP01 | Change in the name or title of a patent holder |