CN110067537A - 用于煤层气排采的小变温吸附曲线制作方法 - Google Patents
用于煤层气排采的小变温吸附曲线制作方法 Download PDFInfo
- Publication number
- CN110067537A CN110067537A CN201910161492.5A CN201910161492A CN110067537A CN 110067537 A CN110067537 A CN 110067537A CN 201910161492 A CN201910161492 A CN 201910161492A CN 110067537 A CN110067537 A CN 110067537A
- Authority
- CN
- China
- Prior art keywords
- temperature
- pressure
- coal
- bed gas
- lower limit
- 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.)
- Pending
Links
- 239000003245 coal Substances 0.000 title claims abstract description 198
- 238000001179 sorption measurement Methods 0.000 title claims abstract description 84
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000000605 extraction Methods 0.000 title claims abstract description 20
- 230000000274 adsorptive effect Effects 0.000 claims abstract description 76
- 230000009466 transformation Effects 0.000 claims abstract description 14
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000010521 absorption reaction Methods 0.000 description 23
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 14
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 9
- 238000012512 characterization method Methods 0.000 description 7
- 238000003795 desorption Methods 0.000 description 7
- 238000005065 mining Methods 0.000 description 7
- 238000002474 experimental method Methods 0.000 description 6
- 241001269238 Data Species 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 239000007788 liquid Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000003337 fertilizer Substances 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 238000003860 storage Methods 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/12—Methods or apparatus for controlling the flow of the obtained fluid to or in wells
- E21B43/121—Lifting well fluids
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)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明提供一种用于煤层气排采的小变温吸附曲线制作方法,包括:设定煤储层温度为上限温度;设定煤层气井筒温度为下限温度;设定上限吸附压力为M兆帕;设定下限吸附压力为N兆帕;计算小变温吸附的变温变压比;根据下限吸附压力、下限温度和变温变压比,计算每增加一个单位吸附压力后的摄氏温度的温度值,得到M‑N组温度值;将计算所得的温度值由摄氏温度℃换算成热力学温度K;将M‑N组热力学温度以及对应的吸附压力值分别代入该种煤压力温度吸附方程,吸附量;作小变温吸附曲线。与相关技术相比,本发明提供的用于煤层气排采的小变温吸附曲线制作方法制作出的小变温吸附曲线所表征的煤层气更切合实际情况,实用性强且准确度高。
Description
技术领域
本发明涉及小变温吸附曲线制作方法领域,尤其涉及一种用于煤层气排采的小变温吸附曲线制作方法。
背景技术
煤层气是赋存于煤储层中以甲烷为主要成分,并吸附在煤基质颗粒表面为主的烃类气体。煤层气的地面抽采是指在煤储层的地面往下钻一井直至预定深度(穿过煤储层)用于排采煤层气的操作。煤储层的含水为承压水。煤储层与煤层气井筒之间有许多裂隙、节理可以作为承压水的渗流通道。煤层气的生产是通过抽排煤层气井筒的水,以降低煤储层压力,并使煤层气解吸、扩散、渗流至煤层气井中而实现。煤层气井排水前,井筒中的液面高度即为煤储层中承压水的水头高度,此时没有压力差,不会产生地下水流动。当煤层气井开始排水后,井筒中的液面不断下降,在煤储层和煤层气井筒之间形成逐渐下降的压力差。煤储层的承压水从压力高的煤储层流向压力低的煤层气井筒,使得煤储层中的压力不断下降,并逐渐向外扩展,形成一个以煤层气井筒为水头最低点中心的压降漏斗。当煤层气井筒的压力小于煤储层临界解吸压力以后形成煤层气解吸压降漏斗,压降漏斗内的吸附煤层气从煤的表面解吸出来。根据所形成的煤层气压降漏斗体积,结合计算吸附体积随吸附条件(煤层温度与煤层压力)的吸附方程,求出该井所能产出的煤层气总量。尽管所有地下煤层都可能蕴藏煤层气,但不是所有地下煤层的煤层气都具有商业开采价值。因此在煤层气商业开采潜力评价中,煤的吸附能力、解吸能力就是不可或缺的指标。
吸附量是吸附温度和吸附压力的函数,同时甲烷在煤中的解吸属于吸热反应,即随着甲烷的解吸煤层中的温度会局部下降,所以煤层气井筒与煤储层之间既有正压力梯度(煤储层又有高于煤层气井筒的压力),又有正温度梯度(煤储层有高于煤层气井筒的温度)。现在采用的Langmuir等温吸附方程,也称为兰氏等温吸附方程,是以完全不考虑煤层气井筒与煤储层之间的温度梯度为代价,仅计算在某一特定温度下,吸附量与吸附压力的关系。这与煤层气井筒与煤储层之间既有正压力梯度,又有正温度梯度的实际情况不符,完全不考虑煤层气井筒与煤储层之间的温度梯度的Langmuir等温吸附方程必然产生一定的误差。
等温吸附实验是在实验室里按选定的一个或数个温度做的实验测试。在实测的温度上,用Langmuir等温吸附方程计算吸附量与压力有精确的函数关系,得到两个参数(兰氏体积VL、兰氏压力PL)。兰氏体积VL是煤所有表面积都已吸附时(最大)的量,而兰氏压力PL是吸附速率与解吸速率之比。尽管做再多个温度的等温吸附实验,但Langmuir等温吸附方程的温度函数的连续性是不存在。这必然产生煤储层的实际温度与选择实验室的测试温度并不匹配,也会产生一定的误差。
当煤储层的实际温度与选择实验室的多个测试温度不同时,也有根据多个测试温度下兰氏参数变化的大小和变化方向用差分法计算煤储层的实际温度下的兰氏参数,从而间接得到在煤储层的实际温度下吸附量与吸附压力的关系。使用差分法精确计算煤储层的实际温度下Langmuir等温吸附方程的兰氏体积VL和兰氏压力PL的前提是温度对这两个参数的影响是线性的。但是大量的实验证实吸附温度对兰氏体积VL和兰氏压力PL的影响不是线性的,因此用差分法计算煤储层的实际温度下的兰氏参数必然也会产生一定的误差,进而对煤层气商业开采潜力评价会产生不同程度的误判。
因此,有必要提供一种新的用于煤层气排采的小变温吸附曲线制作方法以解决上述问题。
发明内容
本发明需要解决的技术问题是提供一种表征的煤层气更切合实际情况,实用性强且准确度高的用于煤层气排采的小变温吸附曲线制作方法。
为解决上述问题,本发明提供了一种用于煤层气排采的小变温吸附曲线制作方法,该方法包括如下步骤:
步骤S1、设定煤储层温度为上限温度。
步骤S2、设定煤层气井筒温度为下限温度。
步骤S3、设定上限吸附压力为M兆帕。
步骤S4、设定下限吸附压力为N兆帕。
步骤S5、计算小变温吸附的变温变压比,单位℃/兆帕,
步骤S6、根据所述下限吸附压力、所述下限温度和所述变温变压比,计算由所述下限吸附压力直至上限吸附压力范围内,每增加一个单位吸附压力后的摄氏温度的温度值,得到M-N组所述温度值。
步骤S7、将计算所得的所述温度值由摄氏温度℃换算成热力学温度K:
K=273+℃。
步骤S8、将M-N组所述热力学温度以及对应的吸附压力值分别代入该种煤压力温度吸附方程,得M-N组相应小变温压力的吸附量。
步骤S9、将对应的M-N组吸附压力值、所述吸附量以及在所述下限吸附压力时的吸附量为零,共M-N+1组对应数据作X煤下限温度℃-上限温度℃的小变温吸附曲线。
优选的,所述上限吸附压力为M=14。
优选的,所述上限吸附压力为N=0。
与相关技术比较,本发明的用于煤层气排采的小变温吸附曲线制作方法制作出的小变温吸附曲线所表征的煤层气更切合实际情况,煤层气井筒与煤储层之间既有正压力梯度,又有正温度梯度;由排采工程技术人员自己确定煤层气井筒温度、煤储层温度、和吸附压力上限,实用性更好,从而达到准确度高的目的。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图,其中:
图1为本发明的用于煤层气排采的小变温吸附曲线制作方法的流程框图。
图2为本发明的用于煤层气排采的小变温吸附曲线制作方法实施例一的小变温吸附曲线图;
图3为本发明的用于煤层气排采的小变温吸附曲线制作方法实施例二的小变温吸附曲线图;
图4为本发明的用于煤层气排采的小变温吸附曲线制作方法实施例三的小变温吸附曲线图;
图5为本发明的用于煤层气排采的小变温吸附曲线制作方法实施例四的小变温吸附曲线图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
实施例一
请参阅图1-2,以下结合具体实例进行说明:
以山西南峪矿肥煤(NY)为例:利用已有的压力温度吸附方程,根据已测的山西南峪矿肥煤(NY)所测30℃和50℃的等温吸附数据回归计算出表征肥煤的几何形体常数A,表征肥煤吸附表面的吸附流量系数B,衡量温度对肥煤吸附量影响的参数,衡量压力对肥煤吸附量影响的参数,如表1-表2所示:
表1肥煤的基本性质和平衡水条件下煤等温吸附结果
表2肥煤的基本性质和在30℃-50℃区间的压力温度吸附方程参数
序号 | 采样地点 | 煤级 | A | B | ||
NY | 山西南峪矿 | 肥煤 | 0.168 | 0.002031 | 953.3 | 0.5342 |
表1列出肥煤的基本性质和平衡水条件下煤等温吸附结果;
表2列出肥煤在30℃至50℃表征煤的几何形体常数A,表征煤吸附表面的吸附流量系数B,衡量温度对煤吸附量影响的参数,衡量压力对煤吸附量影响的参数。
如果有一煤层气井在山西南峪矿的煤田里排采操作。虽然知道煤是肥煤,并且有30℃等温吸附的兰氏体积VL和兰氏压力PL和50℃等温吸附的兰氏体积Vl和兰氏压力Pl,其基本性质和等温吸附如表1中所列。但是现场测得煤储层的温度的确是28℃,而煤层气井筒的温度却是28℃。按本发明的方法步骤计算并制作28℃-30℃小变温吸附曲线,该方法包括如下步骤:
步骤S1、设定煤储层温度为上限温度。
根据实际情况,因现场测得煤储层的温度是30℃,即所述上限温度为30℃。
步骤S2、设定煤层气井筒温度为下限温度。
根据实际情况,因现场测得煤层气井筒的温度是28℃,即所述下限温度为28℃。
步骤S3、设定上限吸附压力为M兆帕。
具体的,本实施方式中,根据文献建议设定上限吸附压力为14兆帕,即M=14。
步骤S4、设定下限吸附压力为N兆帕。
具体的,本实施方式中,根据实际情况,因在煤层气井筒内不可能出现负压,则设定下限吸附压力为0兆帕,即N=0。
步骤S5、计算小变温吸附的变温变压比,单位℃/兆帕,
其中,(上限温度-下限温度)为T,(上限吸附压力-下限吸附压力)为P,本实施方式中,T/P=0.1428℃/兆帕。
步骤S6、根据所述下限吸附压力、所述下限温度和所述变温变压比,计算由所述下限吸附压力直至上限吸附压力范围内,每增加一个单位吸附压力后的摄氏温度的温度值,得到M-N组所述温度值。
即由下限吸附压力开始,计算每增加一个单位吸附压力(兆帕)后的摄氏温度的温度值(℃)。本实施方式中则得到14组所述温度值。而该14组所述温度值以及对应的吸附压力具有正压力梯度和正温度梯度的特征。
步骤S7、将计算所得的所述温度值由摄氏温度℃换算成热力学温度K:
K=273+℃。
步骤S8、将14组所述热力学温度以及对应的吸附压力值分别代入该种煤压力温度吸附方程:
其中,压力温度吸附过程中,V是吸附量,单位为cm3/g;M是分子量,甲烷的分子量为16;T是热力学温度,单位为k;P是压力,单位MPa;其几何形体常数A,表征肥煤吸附表面的吸附流量系数B,衡量温度对肥煤吸附量影响的参数,衡量压力对肥煤吸附量影响的参数列于上述表2。
由此得到14组相应小变温压力的吸附量。
步骤S9、将对应的14组压力值、所述吸附量以及在所述下限吸附压力时的吸附量为零,共15组对应数据作“肥煤28℃-30℃的小变温吸附曲线,如图2所示。图2中所包括的30℃和50℃等温吸附线都是按表1数据所作,用于参考比较。
实施例二
以山西南峪矿肥煤(NY)为例,在山西南峪矿的煤田里排采操作。虽然知道煤是肥煤,并且有30℃等温吸附的兰氏体积VL和兰氏压力PL和50℃等温吸附的兰氏体积Vl和兰氏压力Pl,其基本性质和等温吸附如上述实施例一中表1所列。但是现场测得煤储层的温度的确是40℃,而煤层气井筒的温度却是37℃。按本发明的方法步骤计算并制作37℃-40℃小变温吸附曲线,该方法包括如下步骤:
步骤S1、设定煤储层温度为上限温度。
根据实际情况,因现场测得煤储层的温度是40℃,即所述上限温度为40℃。
步骤S2、设定煤层气井筒温度为下限温度。
根据实际情况,因现场测得煤层气井筒的温度是37℃,即所述下限温度为37℃。
步骤S3、设定上限吸附压力为M兆帕。
具体的,本实施方式中,根据文献建议设定上限吸附压力为14兆帕,即M=14。
步骤S4、设定下限吸附压力为N兆帕。
具体的,本实施方式中,根据实际情况,因在煤层气井筒内不可能出现负压,则设定下限吸附压力为0兆帕,即N=0。
步骤S5、计算小变温吸附的变温变压比,单位℃/兆帕,
其中,(上限温度-下限温度)为T,(上限吸附压力-下限吸附压力)为P,本实施方式中,T/P=0.2143℃/兆帕。
步骤S6、根据所述下限吸附压力、所述下限温度和所述变温变压比,计算由所述下限吸附压力直至上限吸附压力范围内,每增加一个单位吸附压力后的摄氏温度的温度值,得到M-N组所述温度值。
即由下限吸附压力开始,计算每增加一个单位吸附压力(兆帕)后的摄氏温度的温度值(℃)。本实施方式中则得到14组所述温度值。而该14组所述温度值以及对应的吸附压力具有正压力梯度和正温度梯度的特征。
步骤S7、将计算所得的所述温度值由摄氏温度℃换算成热力学温度K:
K=273+℃。
步骤S8、将14组所述热力学温度以及对应的吸附压力值分别代入该种煤压力温度吸附方程:
其中,压力温度吸附过程中,V是吸附量,单位为cm3/g;M是分子量,甲烷的分子量为16;T是热力学温度,单位为k;P是压力,单位MPa;其几何形体常数A,表征肥煤吸附表面的吸附流量系数B,衡量温度对肥煤吸附量影响的参数,衡量压力对肥煤吸附量影响的参数列于上述实施例一中的表2。
由此得到14组相应小变温压力的吸附量。
步骤S9、将对应的14组压力值、所述吸附量以及在所述下限吸附压力时的吸附量为零,共15组对应数据作“肥煤37℃-40℃的小变温吸附曲线,如图3所示。图3中所包括的30℃和50℃等温吸附线都是按上述实施例一中表1数据所作,用于参考比较。
实施例三
以陕西象山矿贫煤(XS)为例:利用已有的压力温度吸附方程,根据已测的陕西象山矿贫煤(XS)所测30℃和50℃的等温吸附数据回归计算出表征贫煤的几何形体常数A,表征贫煤吸附表面的吸附流量系数B,衡量温度对贫煤吸附量影响的参数,衡量压力对贫煤吸附量影响的参数,如表3-表4所示:
表3贫煤的基本性质和平衡水条件下煤等温吸附结果
表4贫煤的基本性质和在30℃-50℃区间的压力温度吸附方程参数
序号 | 采样地点 | 煤级 | A | B | ||
XS | 陕西象山矿 | 贫煤 | 0.168 | 0.002059 | 1149.5 | 0.4644 |
表3列出贫煤的基本性质和平衡水条件下煤等温吸附结果;
表4列出贫煤在30℃至50℃表征煤的几何形体常数A,表征煤吸附表面的吸附流量系数B,衡量温度对煤吸附量影响的参数,衡量压力对煤吸附量影响的参数。
如果有一煤层气井在陕西象山矿的煤田里排采操作。虽然知道煤是贫煤,并且有30℃等温吸附的兰氏体积VL和兰氏压力PL和50℃等温吸附的兰氏体积Vl和兰氏压力Pl,其基本性质和等温吸附如表3中所列。但是现场测得煤储层的温度的确是30℃,而煤层气井筒的温度却是28℃。按本发明的方法步骤计算并制作28℃-30℃小变温吸附曲线,该方法包括如下步骤:
步骤S1、设定煤储层温度为上限温度。
根据实际情况,因现场测得煤储层的温度是30℃,即所述上限温度为30℃。
步骤S2、设定煤层气井筒温度为下限温度。
根据实际情况,因现场测得煤层气井筒的温度是28℃,即所述下限温度为28℃。
步骤S3、设定上限吸附压力为M兆帕。
具体的,本实施方式中,根据文献建议设定上限吸附压力为14兆帕,即M=14。
步骤S4、设定下限吸附压力为N兆帕。
具体的,本实施方式中,根据实际情况,因在煤层气井筒内不可能出现负压,则设定下限吸附压力为0兆帕,即N=0。
步骤S5、计算小变温吸附的变温变压比,单位℃/兆帕,
其中,(上限温度-下限温度)为T,(上限吸附压力-下限吸附压力)为P,本实施方式中,T/P=0.1428℃/兆帕。
步骤S6、根据所述下限吸附压力、所述下限温度和所述变温变压比,计算由所述下限吸附压力直至上限吸附压力范围内,每增加一个单位吸附压力后的摄氏温度的温度值,得到M-N组所述温度值。
即由下限吸附压力开始,计算每增加一个单位吸附压力(兆帕)后的摄氏温度的温度值(℃)。本实施方式中则得到14组所述温度值。而该14组所述温度值以及对应的吸附压力具有正压力梯度和正温度梯度的特征。
步骤S7、将计算所得的所述温度值由摄氏温度℃换算成热力学温度K:
K=273+℃。
步骤S8、将14组所述热力学温度以及对应的吸附压力值分别代入该种煤压力温度吸附方程:
其中,压力温度吸附过程中,V是吸附量,单位为cm3/g;M是分子量,甲烷的分子量为16;T是热力学温度,单位为k;P是压力,单位MPa;其几何形体常数A,表征贫煤吸附表面的吸附流量系数B,衡量温度对贫煤吸附量影响的参数,衡量压力对贫煤吸附量影响的参数列于上述表4。
由此得到14组相应小变温压力的吸附量。
步骤S9、将对应的14组压力值、所述吸附量以及在所述下限吸附压力时的吸附量为零,共15组对应数据作“贫煤28℃-30℃的小变温吸附曲线,如图4所示。图4中所包括的30℃和50℃等温吸附线都是按上述表3数据所作,用于参考比较。
实施例四
以陕西象山矿贫煤(XS)为例,在陕西象山矿的煤田里排采操作。虽然知道煤是贫煤,并且有30℃等温吸附的兰氏体积VL和兰氏压力PL和50℃等温吸附的兰氏体积Vl和兰氏压力Pl,其基本性质和等温吸附如上述实施例三中表3所列。但是现场测得煤储层的温度的确是40℃,而煤层气井筒的温度却是37℃。按本发明的方法步骤计算并制作37℃-40℃小变温吸附曲线,该方法包括如下步骤:
步骤S1、设定煤储层温度为上限温度。
根据实际情况,因现场测得煤储层的温度是40℃,即所述上限温度为40℃。
步骤S2、设定煤层气井筒温度为下限温度。
根据实际情况,因现场测得煤层气井筒的温度是37℃,即所述下限温度为37℃。
步骤S3、设定上限吸附压力为M兆帕。
具体的,本实施方式中,根据文献建议设定上限吸附压力为14兆帕,即M=14。
步骤S4、设定下限吸附压力为N兆帕。
具体的,本实施方式中,根据实际情况,因在煤层气井筒内不可能出现负压,则设定下限吸附压力为0兆帕,即N=0。
步骤S5、计算小变温吸附的变温变压比,单位℃/兆帕,
其中,(上限温度-下限温度)为T,(上限吸附压力-下限吸附压力)为P,本实施方式中,T/P=0.2143℃/兆帕。
步骤S6、根据所述下限吸附压力、所述下限温度和所述变温变压比,计算由所述下限吸附压力直至上限吸附压力范围内,每增加一个单位吸附压力后的摄氏温度的温度值,得到M-N组所述温度值。
即由下限吸附压力开始,计算每增加一个单位吸附压力(兆帕)后的摄氏温度的温度值(℃)。本实施方式中则得到14组所述温度值。而该14组所述温度值以及对应的吸附压力具有正压力梯度和正温度梯度的特征。
步骤S7、将计算所得的所述温度值由摄氏温度℃换算成热力学温度K:
K=273+℃。
步骤S8、将14组所述热力学温度以及对应的吸附压力值分别代入该种煤压力温度吸附方程:
其中,压力温度吸附过程中,V是吸附量,单位为cm3/g;M是分子量,甲烷的分子量为16;T是热力学温度,单位为k;P是压力,单位MPa;其几何形体常数A,表征贫煤吸附表面的吸附流量系数B,衡量温度对贫煤吸附量影响的参数,衡量压力对贫煤吸附量影响的参数列于上述实施例三中的表4。
由此得到14组相应小变温压力的吸附量。
步骤S9、将对应的14组压力值、所述吸附量以及在所述下限吸附压力时的吸附量为零,共15组对应数据作“贫煤37℃-40℃的小变温吸附曲线,如图5所示。图5中所包括的30℃和50℃等温吸附线都是依据实施例三中的表3数据所作,用于参考比较。
与相关技术比较,本发明的用于煤层气排采的小变温吸附曲线制作方法制作出的小变温吸附曲线所表征的煤层气更切合实际情况,煤层气井筒与煤储层之间既有正压力梯度,又有正温度梯度;由排采工程技术人员自己确定煤层气井筒温度、煤储层温度、和吸附压力上限,实用性更好,从而达到准确度高的目的。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其它相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (3)
1.一种用于煤层气排采的小变温吸附曲线制作方法,该方法包括如下步骤:
步骤S1、设定煤储层温度为上限温度;
步骤S2、设定煤层气井筒温度为下限温度;
步骤S3、设定上限吸附压力为M兆帕;
步骤S4、设定下限吸附压力为N兆帕;
步骤S5、计算小变温吸附的变温变压比,单位℃/兆帕,
步骤S6、根据所述下限吸附压力、所述下限温度和所述变温变压比,计算由所述下限吸附压力直至上限吸附压力范围内,每增加一个单位吸附压力后的摄氏温度的温度值,得到M-N组所述温度值;
步骤S7、将计算所得的所述温度值由摄氏温度℃换算成热力学温度K;
K=273+℃;
步骤S8、将M-N组所述热力学温度以及对应的吸附压力值分别代入该种煤压力温度吸附方程,得M-N组相应小变温压力的吸附量;
步骤S9、将对应的M-N组吸附压力值、所述吸附量以及在所述下限吸附压力时的吸附量为零,共M-N+1组对应数据作X煤下限温度℃-上限温度℃的小变温吸附曲线。
2.根据权利要求1所述的用于煤层气排采的小变温吸附曲线制作方法,其特征在于,所述上限吸附压力为M=14。
3.根据权利要求2所述的用于煤层气排采的小变温吸附曲线制作方法,其特征在于,所述下限吸附压力为N=0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910161492.5A CN110067537A (zh) | 2019-03-04 | 2019-03-04 | 用于煤层气排采的小变温吸附曲线制作方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910161492.5A CN110067537A (zh) | 2019-03-04 | 2019-03-04 | 用于煤层气排采的小变温吸附曲线制作方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110067537A true CN110067537A (zh) | 2019-07-30 |
Family
ID=67366033
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910161492.5A Pending CN110067537A (zh) | 2019-03-04 | 2019-03-04 | 用于煤层气排采的小变温吸附曲线制作方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110067537A (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5470823A (en) * | 1993-05-03 | 1995-11-28 | Exxon Chemical Patents Inc. | Stimulation of coalbed methane production |
CN103592687A (zh) * | 2013-11-22 | 2014-02-19 | 中国石油化工集团公司 | 煤岩吸附气含量的定量计算方法 |
CN105930667A (zh) * | 2016-04-27 | 2016-09-07 | 中国矿业大学银川学院 | 一种计算气井煤层气吸附量的方法 |
-
2019
- 2019-03-04 CN CN201910161492.5A patent/CN110067537A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5470823A (en) * | 1993-05-03 | 1995-11-28 | Exxon Chemical Patents Inc. | Stimulation of coalbed methane production |
CN103592687A (zh) * | 2013-11-22 | 2014-02-19 | 中国石油化工集团公司 | 煤岩吸附气含量的定量计算方法 |
CN105930667A (zh) * | 2016-04-27 | 2016-09-07 | 中国矿业大学银川学院 | 一种计算气井煤层气吸附量的方法 |
Non-Patent Citations (1)
Title |
---|
张学梅: "变温变压下煤层气吸附量的预测", 《变温变压下煤层气吸附量的预测》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104389594B (zh) | 页岩气井产能评价预测方法 | |
CN111353205B (zh) | 用于致密气藏产水气井地层压力和动态产能的计算方法 | |
Gamadi et al. | Compositional simulation evaluation of EOR potential in shale oil reservoirs by cyclic natural gas injection | |
CN114372352B (zh) | 渗流-温度双场耦合数值模拟预测复杂断块油藏储气库调峰能力方法 | |
CN103161435B (zh) | 一种稠油热采直井试井解释方法 | |
CN103902827B (zh) | 一种碳酸盐岩水平井流动单元划分方法 | |
CN108133080B (zh) | 一种考虑非稳态吸附的非均质裂缝页岩气藏数值模拟方法 | |
CN106198338A (zh) | 页岩储层压裂裂缝应力敏感性测试装置及使用其的方法 | |
CN108959767B (zh) | 一种窄河道型气藏不同井型凝析油伤害数值模拟方法 | |
CN107526891B (zh) | 一种聚合物驱大孔道油藏试井分析方法 | |
CN105893679B (zh) | 低产水平井续流修正试井解释方法 | |
CN106522928B (zh) | 一种酸化压裂后停泵测井口压降不稳定试井方法 | |
CN109242364B (zh) | 一种高温高压气井模拟井筒体积置换产能评价法 | |
CN110984973A (zh) | 缝洞型碳酸盐岩气藏单井控制储量的确定方法 | |
CN102809526A (zh) | 一种测量二氧化碳在饱和油岩心中扩散系数的方法 | |
CN104330344A (zh) | 岩心气水两相渗流动态测试方法与装置 | |
CN106761613B (zh) | 一种二氧化碳驱替前缘的试井确定方法 | |
CN105651665A (zh) | 一种钻完井液对岩心油水渗透率影响评价方法 | |
CN103161436B (zh) | 一种稠油热采水平井试井解释方法 | |
CN109209307A (zh) | 一种定量分析低渗透油藏注水开发效果的方法 | |
CN105804713B (zh) | 一种快速确定注水井各小层井口注水启动压力的方法 | |
CN110067537A (zh) | 用于煤层气排采的小变温吸附曲线制作方法 | |
CN106321076A (zh) | 一种注水井启动压力测试方法 | |
CN110965992B (zh) | 一种确定地层含气原油粘度的方法 | |
CN107975356B (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190730 |