CN107796457A - 一种边远井高压槽车充装时间的计算方法 - Google Patents
一种边远井高压槽车充装时间的计算方法 Download PDFInfo
- Publication number
- CN107796457A CN107796457A CN201710972291.4A CN201710972291A CN107796457A CN 107796457 A CN107796457 A CN 107796457A CN 201710972291 A CN201710972291 A CN 201710972291A CN 107796457 A CN107796457 A CN 107796457A
- Authority
- CN
- China
- Prior art keywords
- pressure
- tank car
- natural gas
- delta
- pipeline
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
- G01F1/76—Devices for measuring mass flow of a fluid or a fluent solid material
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Fluid Mechanics (AREA)
- Filling Or Discharging Of Gas Storage Vessels (AREA)
Abstract
一种边远井高压槽车充装时间的计算方法，包括以下步骤：1)收集高压槽车充装系统的基础数据；2)取初始化循环次数n＝1，单位时长Δt，充装结束条件ΔP；3)判定管道中天然气的流态，并计算相应流态下的质量流量；4)计算单位时长后高压槽车的充装质量、井口压力和槽车压力；5)判定单位时长后井口压力与槽车压力之差是否大于ΔP；若是，则返回步骤3)进行循环，并按照n＝n+1累计循环次数，判定第n个Δt时长后井口压力与槽车压力之差是否大于ΔP，若否，则输出充装时间t＝nΔt，迭代结束。本发明的目的在于针对边远井难以实现高效低成本的回收提供一种高压槽车充装时间的计算方法，该方法为槽车选型及数量配置提供可靠依据，有助于提高回收效率及经济效益。
Description
技术领域
本发明涉及一种边远井高压槽车充装时间的计算方法，属于高压天然气回收工艺方案设计的技术领域。
背景技术
边远井内的天然气作为一种极具开发价值的非常规天然气资源，在世界天然气资源总量中所占比例正逐年攀升，受到世界各国的高度关注。但由于边远井具有产气含水量高、井间距大，距附近集气站距离一般也都在20km以上，交通不便，缺少外部供水、供电和天然气外输管道等非常规特征。若采用常规管道集气工艺，需铺设较长距离集气管道，导致经济效益差，且集气半径、流动保障等也难以满足集输工艺要求，故必须实现低成本策略才能实现盈利。目前，针对这类天然气的开发回收方式，主要有LNG、CNG、天然气就地发电、油田套管气利用、天然气吸附储存及天然气水合物固化存储等回收方式。其中，LNG技术因其液化装置及配套工艺投资高昂，要求气井产量高并能保持较长稳产期，故对于典型的低压低渗致密型气田内的边远井来讲是完全没有可行性的。因此，为了充分挖掘气田产能潜力、提高经济效益，通常采用高压槽车充装模式对偏远井进行低成本的开发利用。
针对边远井内天然气的回收，槽车数量测算是其方案设计的重要内容，而槽车数量测算与槽车充装时间又有很大的关联。若能够较为准确计算槽车充装时间，以此为依据，测算槽车数量，并进行合理配置，势必有助于提高回收效率与经济效益。然而，目前国内外关于边远井高压槽车充装时间的计算方法尚未提出。因此，本文提出了一种偏远井高压槽车充装时间的计算方法，该方法为槽车选型及数量配置提供可靠的依据，有助于提高回收效率与经济效益。
发明内容
本发明的目的是针对边远井难以实现高效回收的技术难题，提供一种边远井高压槽车充装时间的计算方法，该方法为槽车选型及数量配置提供可靠依据，有助于提高回收效率及经济效益。
为达上述目的，本发明技术方案的具体步骤如下：
步骤a、根据边远井高压槽车充装系统的设计资料，回收方案和现场勘测，收集脱水装置出口与槽车间的管道内径、管长、管道水力摩阻系数、井筒水容积、槽车水容积、井口初始压力、槽车内初始储存压力及天然气平均温度等基础数据；
步骤b、取初始循环次数n＝1，单位时长为Δt，充装结束条件为ΔP；
步骤c、判定管道中天然气的流态，并计算相应流态下的质量流量；
步骤d、计算第1个Δt时长后高压槽车的充装质量、井口压力和槽车压力；
步骤e、判定第1个Δt时长后井口压力与槽车压力之差是否大于ΔP；若是，则返回步骤c进行循环，并按照n＝n+1累计循环次数，判定第n个Δt时长后井口压力与槽车压力之差是否大于ΔP，若否，则输出高压槽车的充装时间t＝nΔt，迭代结束。
步骤c中所述的判定管道中天然气的流态，并计算相应流态下的质量流量；
若PQ0/PC0≥1.841，则管道中的天然气处于临界状态，质量流量的计算公式为：
式中，Mc—临界状态下天然气质量流量，kg/s；
PQ0—脱水装置出口的初始压力，Pa；
D—管道内径，m；
k—天然气的绝热系数，取1.314，无量纲；
MN—天然气平均摩尔质量，取17.62kg/kmol；
Z—压缩因子，无量纲；
R—气体常数，取8314.3J/(kmol·K)；
T—天然气的平均温度，K；
若PQ0/PC0<1.841，则管道中的天然气处于亚临界状态，质量流量的计算公式为：
式中，Ms—亚临界状态下天然气质量流量，kg/s；
PQ0—脱水装置出口的初始压力，Pa；
PC0—高压槽车内的初始储存压力，Pa；
D—管道内径，m；
λ—管道水力摩阻系数，无量纲；
L—脱水装置出口与高压槽车间的管长，m；
Z—压缩因子，无量纲；
R—气体常数，取8314.3/气体摩尔质量，m2/(s2·K)；
T—天然气的平均温度，K。
步骤d中所述第1个Δt时长内高压槽车充装质量的计算公式为：
式中，Δm1—第1个Δt时长内高压槽车的充装质量，kg；
Mc—临界状态下天然气质量流量，kg/s；
Ms—亚临界状态下天然气质量流量，kg/s；
Δt—单位时长，s；
所述第1个Δt时长后井口压力的计算公式为：
式中，P0—井口的初始压力，Pa；
V—井筒容积，m3；
m0—井筒内气体的初始储存质量，kg；
P1—第1个Δt时长后井口的压力，Pa；
Δm1—第1个Δt时长内高压槽车充装的质量，kg；
Z—压缩因子，无量纲；
R—气体常数，取8314.3/气体摩尔质量，m2/(s2·K)；
T—天然气的平均温度，K；
所述第1个Δt时长后槽车压力的计算公式为：
式中，PC0—高压槽车的初始储存压力，Pa；
VC—高压槽车容积，m3；
mC0—高压槽车内初始储存质量，kg；
PC1—第1个Δt时长后槽车的压力，Pa；
Δm1—第1个Δt时长内高压槽车的充装质量，kg；
Z—压缩因子，无量纲；
R—气体常数，取8314.3/气体摩尔质量，m2/(s2·K)；
T—天然气的平均温度，K。
本发明的有益效果是：针对边远井难以实现高效回收的技术难题，提供一种边远井高压槽车充装时间的计算方法，从而为槽车选型及数量配置提供可靠依据，有助于回收效率及经济效益的提高。该方法可用于国内煤层气、页岩气、致密砂岩气等非常规天然气田的回收开发。
附图说明
图1为边远井高压槽车充装时间计算方法的流程图；
图2为高压槽车的充装模型。
具体实施方式
下面结合实施例对本发明的具体实施方式作进一步详细描述，来说明本发明具体实施步骤和效果。以下实施例用于说明本发明，但不用来限制本发明的范围。
本发明提供了一种边远井高压槽车充装时间的计算方法，如图1所示，该方法包括以下步骤：
步骤a、根据边远井高压槽车充装系统的设计资料，回收方案和现场勘测，收集脱水装置出口与槽车间的管道内径、管长、管道水力摩阻系数、井筒水容积、槽车水容积、井口初始压力、槽车内初始储存压力及天然气平均温度等基础数据；
步骤b、取初始循环次数n＝1，单位时长为Δt，充装结束条件为ΔP；
步骤c、判定管道中天然气的流态，并计算相应流态下的质量流量；
若PQ0/PC0≥1.841，则管道中的天然气处于临界状态，质量流量的计算公式为：
式中，Mc—临界状态下天然气质量流量，kg/s；
PQ0—脱水装置出口的初始压力，Pa；
D—管道内径，m；
k—天然气的绝热系数，取1.314，无量纲；
MN—天然气平均摩尔质量，取17.62kg/kmol；
Z—压缩因子，无量纲；
R—气体常数，取8314.3J/(kmol·K)；
T—天然气的平均温度，K；
若PQ0/PC0<1.841，则管道中的天然气处于亚临界状态，质量流量的计算公式为：
式中，Ms—亚临界状态下天然气质量流量，kg/s；
PQ0—脱水装置出口的初始压力，Pa；
PC0—高压槽车内的初始储存压力，Pa；
D—管道内径，m；
λ—管道水力摩阻系数，无量纲；
L—脱水装置出口与高压槽车间的管长，m；
Z—压缩因子，无量纲；
R—气体常数，取8314.3/气体摩尔质量，m2/(s2·K)；
T—天然气的平均温度，K；
步骤d、第1个Δt时长内高压槽车的充装质量计算公式为：
式中，Δm1—第1个Δt时长内高压槽车的充装质量，kg；
Mc—临界状态下天然气质量流量，kg/s；
Ms—亚临界状态下天然气质量流量，kg/s；
Δt—单位时长，s；
第1个Δt时长后井口压力的计算公式为：
式中，P0—井口的初始压力，Pa；
V—井筒容积，m3；
m0—井筒内气体的初始储存质量，kg；
P1—第1个Δt时长后井口的压力，Pa；
Δm1—第1个Δt时长内高压槽车的充装质量，kg；
Z—压缩因子，无量纲；
R—气体常数，取8314.3/气体摩尔质量，m2/(s2·K)；
T—天然气的平均温度，K；
第1个Δt时长后槽车压力的计算公式为：
式中，PC0—高压槽车的初始储存压力，Pa；
VC—高压槽车容积，m3；
mC0—高压槽车内初始储存质量，kg；
PC1—第1个Δt时长后槽车的压力，Pa；
Δm1—第1个Δt时长内高压槽车的充装质量，kg；
Z—压缩因子，无量纲；
R—气体常数，取8314.3/气体摩尔质量，m2/(s2·K)；
T—天然气的平均温度，K；
步骤e、判定第1个Δt时长后井口压力与槽车压力之差是否大于ΔP；若是，则返回步骤c进行循环，并按照n＝n+1累计循环次数，判定第n个Δt时长后井口压力与槽车压力之差是否大于ΔP，若否，则输出高压槽车的充装时间t＝nΔt，迭代结束。
实施例1：以国内某低产井的天然气回收为例，并结合附图来说明本发明的内容及实现原理：
步骤a、根据高压槽车充装系统的设计资料，回收方案和现场勘测，收集高压槽车充装管道系统的基础数据；
如图2所示，脱水装置出口与高压槽车间的管道内径D为25mm、管长L为30m、管道水力摩阻系数λ为0.032、井筒水容积为36.79m3、槽车水容积为22m3、井口初始压力P0为25MPa、槽车内初始储存压力PC0为1.5MPa及天然气平均温度T为293.15K，脱水装置出口压力PQ随充装时间t的变化规律为：PQ＝0.0001t2+0.0470t+5.5655；
步骤b、取初始化循环次数n＝1，单位时长Δt＝5min，充装结束条件ΔP＝1MPa；
步骤c、判定管道中天然气的流态；
PQ0/PC0＝5.57/1.5＝3.71＞1.841，故管道中的天然气处于临界状态；
计算临界状态下天然气质量流量：
步骤d、计算第1个5min内高压槽车的充装质量：
Δm1＝Mc·Δt＝5.28×5×60＝1584kg
计算第1个5min后井口压力：
由P0V＝m0ZRT，得
由P1V＝(m0-Δm1)ZRT，得
计算第1个5min后高压槽车压力：
由PC0VC＝mC0ZRT，得
由PC1VC＝(mC0+Δm1)ZRT，得
步骤e、判定第1个5min后井口压力和槽车压力之差是否大于ΔP，即P1－PC1＝24.95－2.14＝22.81MPa>1MPa，则返回步骤c进行第2轮循环；
步骤c、判定管道中天然气的流态；
PQ1/PC1＝5.80/2.14＝2.71＞1.841，故管道中的天然气处于临界状态；
计算临界状态下天然气质量流量：
步骤d、计算第2个5min内高压槽车的充装质量：
Δm2＝Mc·Δt＝5.52×5×60＝1656kg
计算第2个5min后气井井口压力：
由P1V＝m1ZRT，得
由P2V＝(m1-Δm2)ZRT，得
计算第2个5min后高压槽车压力：
由PC1VC＝mC1ZRT，得
由PC2VC＝(mC1+Δm2)ZRT，得
步骤e、判定第2个5min后井口压力与槽车压力之差是否大于ΔP，即P2－PC2＝24.88－2.60＝22.28MPa>1MPa，则返回步骤c继续循环，并按照n＝n+1累计循环次数，判定第n个5min后井口压力与槽车压力之差是否大于ΔP，若否，则输出高压槽车的充装时间t＝5n min，迭代结束。
实施例2：以国内某高产井的天然气回收为例，并结合附图来说明本发明的内容及实现原理：步骤a、根据高压槽车充装系统的设计资料，回收方案和现场勘测，收集高压槽车充装管道系统的基础数据；
如图2所示，脱水装置出口与高压槽车间的管道内径D为25mm、管长L为30m、管道水力摩阻系数λ为0.032、井筒水容积为36.79m3、槽车水容积为22m3、井口压力P为25MPa且保持不变，槽车内初始储存压力PC0为1.5MPa及天然气平均温度T为293.15K，脱水装置出口压力PQ为6MPa；
步骤b、取初始化循环次数n＝1，单位时长Δt＝5min，充装结束条件ΔP＝1MPa；
步骤c、判定管道中天然气的流态；
PQ/PC0＝6/1.5＝4＞1.841，故管道中的天然气处于临界状态；
计算临界状态下天然气质量流量：
步骤d、计算第1个5min内高压槽车的充装质量：
Δm1＝Mc·Δt＝5.71×5×60＝1713kg
计算第1个5min后高压槽车压力：
由PC0VC＝mC0ZRT，得
由PC1VC＝(mC0+Δm1)ZRT，得
步骤e、判定第1个5min后井口压力和槽车压力之差是否大于ΔP，即P－PC1＝25－2.19＝22.81MPa>1MPa，则返回步骤c进行第2轮循环；
步骤c、判定管道中天然气的流态；
PQ/PC1＝6/2.19＝2.74＞1.841，故管道中的天然气处于临界状态；
计算临界状态下天然气质量流量：
步骤d、计算第2个5min内高压槽车的充装质量：
Δm2＝Mc·Δt＝5.71×5×60＝1713kg
计算第2个5min后高压槽车压力：
由PC1VC＝mC1ZRT，得
由PC2VC＝(mC1+Δm2)ZRT，得
步骤e、判定第2个5min后井口压力与槽车压力之差是否大于ΔP，即P－PC2＝25－2.65＝22.35MPa>1MPa，则返回步骤c继续循环，并按照n＝n+1累计循环次数，判定第n个5min后井口压力与槽车压力之差是否大于ΔP，若否，输出高压槽车的充装时间t＝5n min，迭代结束。从上述实例可知，基于本发明的一种偏远井高压槽车充装时间的计算方法，可以准确有效的计算出高压槽车的充装时间，进而为槽车选型及数量配置提供可靠依据，同时有助于提高回收效率及经济效益。该方法也可用于国内煤层气、页岩气、致密砂岩气等非常规天然气田的回收开发。
Claims (3)
1.一种边远井高压槽车充装时间的计算方法，其特征包含以下步骤：
步骤a、根据边远井高压槽车充装系统的设计资料，回收方案和现场勘测，收集脱水装置出口与槽车间的管道内径、管长、管道水力摩阻系数、井筒水容积、槽车水容积、井口初始压力、槽车内初始储存压力及天然气平均温度等基础数据；
步骤b、取初始循环次数n＝1，单位时长为Δt，充装结束条件为ΔP；
步骤c、判定管道中天然气的流态，并计算相应流态下的质量流量；
步骤d、计算第1个Δt时长后高压槽车的充装质量、井口压力和槽车压力；
步骤e、判定第1个Δt时长后井口压力与槽车压力之差是否大于ΔP；若是，则返回步骤c进行循环，并按照n＝n+1累计循环次数，判定第n个Δt时长后井口压力与槽车压力之差是否大于ΔP，若否，则输出高压槽车的充装时间t＝nΔt，迭代结束。
2.根据权利要求书1所述的一种边远井高压槽车充装时间的计算方法，其特征在于，步骤c中所述的判定管道中天然气的流态，并计算相应流态下的质量流量；
若PQ0/PC0≥1.841，则管道中的天然气处于临界状态，质量流量的计算公式为：
<mrow>
<msub>
<mi>M</mi>
<mi>c</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>Q</mi>
<mn>0</mn>
</mrow>
</msub>
<msup>
<mi>&pi;D</mi>
<mn>2</mn>
</msup>
</mrow>
<mn>4</mn>
</mfrac>
<msqrt>
<mrow>
<mfrac>
<mrow>
<msub>
<mi>kM</mi>
<mi>N</mi>
</msub>
</mrow>
<mrow>
<mi>Z</mi>
<mi>R</mi>
<mi>T</mi>
</mrow>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mn>2</mn>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式中：Mc—临界状态下天然气质量流量，kg/s；
PQ0—脱水装置出口的初始压力，Pa；
D—管道内径，m；
k—天然气的绝热系数，取1.314，无量纲；
MN—天然气平均摩尔质量，取17.62kg/kmol；
Z—压缩因子，无量纲；
R—气体常数，取8314.3J/(kmol·K)；
T—天然气的平均温度，K；
若PQ0/PC0<1.841，则管道中的天然气处于亚临界状态，质量流量的计算公式为：
<mrow>
<msub>
<mi>M</mi>
<mi>s</mi>
</msub>
<mo>=</mo>
<mfrac>
<mi>&pi;</mi>
<mn>4</mn>
</mfrac>
<msqrt>
<mfrac>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>P</mi>
<mrow>
<mi>Q</mi>
<mn>0</mn>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<msub>
<mi>P</mi>
<mrow>
<mi>C</mi>
<mn>0</mn>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>)</mo>
<msup>
<mi>D</mi>
<mn>4</mn>
</msup>
</mrow>
<mrow>
<mi>Z</mi>
<mi>R</mi>
<mi>T</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>&lambda;</mi>
<mi>L</mi>
</mrow>
<mi>D</mi>
</mfrac>
<mo>+</mo>
<mn>2</mn>
<mi>ln</mi>
<mfrac>
<msub>
<mi>P</mi>
<mrow>
<mi>Q</mi>
<mn>0</mn>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mrow>
<mi>C</mi>
<mn>0</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中：Ms—亚临界状态下天然气质量流量，kg/s；
PQ0—脱水装置出口的初始压力，Pa；
PC0—高压槽车内的初始储存压力，Pa；
D—管道内径，m；
λ—管道水力摩阻系数，无量纲；
L—脱水装置出口与高压槽车间的管长，m；
Z—压缩因子，无量纲；
R—气体常数，取8314.3/气体摩尔质量，m2/(s2·K)；
T—天然气的平均温度，K。
3.根据权利要求书1所述的一种边远井高压槽车充装时间的计算方法，其特征在于，步骤d中所述第1个Δt时长内高压槽车的充装质量的计算公式为：
<mrow>
<msub>
<mi>&Delta;m</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>M</mi>
<mi>c</mi>
</msub>
<mo>&CenterDot;</mo>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>M</mi>
<mi>s</mi>
</msub>
<mo>&CenterDot;</mo>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
式中：Δm1—第1个Δt时长内高压槽车的充装质量，kg；
Mc—临界状态下天然气质量流量，kg/s；
Ms—亚临界状态下天然气质量流量，kg/s；
Δt—单位时长，s；
所述第1个Δt时长后井口压力的计算公式为：
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mn>0</mn>
</msub>
<mi>V</mi>
<mo>=</mo>
<msub>
<mi>m</mi>
<mn>0</mn>
</msub>
<mi>Z</mi>
<mi>R</mi>
<mi>T</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mi>V</mi>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>m</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;m</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mi>Z</mi>
<mi>R</mi>
<mi>T</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
式中：P0—井口的初始压力，Pa；
V—井筒容积，m3；
m0—井筒内气体的初始储存质量，kg；
P1—第1个Δt时长后井口的压力，Pa；
Δm1—第1个Δt时长内高压槽车的充装质量，kg；
Z—压缩因子，无量纲；
R—气体常数，取8314.3/气体摩尔质量，m2/(s2·K)；
T—天然气的平均温度，K；
所述第1个Δt时长后槽车压力的计算公式为：
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>C</mi>
<mn>0</mn>
</mrow>
</msub>
<msub>
<mi>V</mi>
<mi>C</mi>
</msub>
<mo>=</mo>
<msub>
<mi>m</mi>
<mrow>
<mi>C</mi>
<mn>0</mn>
</mrow>
</msub>
<mi>Z</mi>
<mi>R</mi>
<mi>T</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>C</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>V</mi>
<mi>C</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>m</mi>
<mrow>
<mi>C</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&Delta;m</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mi>Z</mi>
<mi>R</mi>
<mi>T</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
式中：PC0—高压槽车的初始储存压力，Pa；
VC—高压槽车容积，m3；
mC0—高压槽车内初始储存质量，kg；
PC1—第1个Δt时长后槽车的压力，Pa
Δm1—第1个Δt时长内高压槽车的充装质量，kg；
Z—压缩因子，无量纲；
R—气体常数，取8314.3/气体摩尔质量，m2/(s2·K)；
T—天然气的平均温度，K。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710972291.4A CN107796457A (zh) | 2017-10-18 | 2017-10-18 | 一种边远井高压槽车充装时间的计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710972291.4A CN107796457A (zh) | 2017-10-18 | 2017-10-18 | 一种边远井高压槽车充装时间的计算方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107796457A true CN107796457A (zh) | 2018-03-13 |
Family
ID=61534236
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710972291.4A Pending CN107796457A (zh) | 2017-10-18 | 2017-10-18 | 一种边远井高压槽车充装时间的计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107796457A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110847883A (zh) * | 2018-07-24 | 2020-02-28 | 中国石油天然气股份有限公司 | 一种井口装置抬升高度预测方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1504401A (zh) * | 2002-12-02 | 2004-06-16 | 安阳钢铁股份有限公司 | 一种低温液体槽车充装方法 |
| CN202938010U (zh) * | 2012-11-27 | 2013-05-15 | 武汉钢铁(集团)公司 | 一种自动充装系统 |
| CN204005206U (zh) * | 2014-05-06 | 2014-12-10 | 武汉柯特达科技发展有限责任公司 | 低温槽车充装系统 |
| CN204201465U (zh) * | 2014-09-18 | 2015-03-11 | 济宁金威煤电有限公司 | 液氯槽车定量充装装置 |
| CN105864630A (zh) * | 2016-04-08 | 2016-08-17 | 海宁市英德赛电子有限公司 | 高纯氨充装系统及利用其充装高纯氨的充装方法 |
-
2017
- 2017-10-18 CN CN201710972291.4A patent/CN107796457A/zh active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1504401A (zh) * | 2002-12-02 | 2004-06-16 | 安阳钢铁股份有限公司 | 一种低温液体槽车充装方法 |
| CN202938010U (zh) * | 2012-11-27 | 2013-05-15 | 武汉钢铁(集团)公司 | 一种自动充装系统 |
| CN204005206U (zh) * | 2014-05-06 | 2014-12-10 | 武汉柯特达科技发展有限责任公司 | 低温槽车充装系统 |
| CN204201465U (zh) * | 2014-09-18 | 2015-03-11 | 济宁金威煤电有限公司 | 液氯槽车定量充装装置 |
| CN105864630A (zh) * | 2016-04-08 | 2016-08-17 | 海宁市英德赛电子有限公司 | 高纯氨充装系统及利用其充装高纯氨的充装方法 |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110847883A (zh) * | 2018-07-24 | 2020-02-28 | 中国石油天然气股份有限公司 | 一种井口装置抬升高度预测方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11414273B2 (en) | System and method for compressed air energy storage | |
| CN204062478U (zh) | 一种bog重新液化装置 | |
| Liu et al. | The role of underground salt caverns for large-scale energy storage: A review and prospects | |
| CN107145696A (zh) | 一种煤层气地上地下耦合求解的模拟方法 | |
| CN107542440A (zh) | 一种u型井微波加热开采天然气水合物的方法 | |
| CN104832131A (zh) | 一种基于机采和集油系统整体的集油参数优选方法 | |
| CN104975827B (zh) | 预测二氧化碳驱油藏指标的物质平衡方法 | |
| CN107542441A (zh) | 一种平行水平井注热水开采天然气水合物的方法 | |
| CN107676061A (zh) | 一种平行水平井开采海上天然气水合物的方法 | |
| CN107676060A (zh) | 一种天然气水合物u型井微波加热开采的方法 | |
| CN103968252A (zh) | 一种中低压天然气集气的注醇工艺方法 | |
| CN107514245A (zh) | 一种天然气水合物排式水平井开采的方法 | |
| CN103114828A (zh) | 一种井下节流工艺参数的设置方法 | |
| Abd et al. | A comprehensive thermodynamic analysis of an integrated solar enhanced oil recovery system for applications in heavy oil fields | |
| Zhang et al. | Research on the hydrate formation in the process of gas phase CO2 pipeline transportation | |
| CN103726815A (zh) | 一种co2驱采出井井筒流态确定及参数优化方法 | |
| CN105426666A (zh) | 天然气水合物分解气体释放速率计算方法及其装置 | |
| CN107575201A (zh) | 一种天然气水合物平行水平井微波加热开采的方法 | |
| CN107796457A (zh) | 一种边远井高压槽车充装时间的计算方法 | |
| CN107701150A (zh) | 一种海上天然气水合物平行水平井开采的方法 | |
| Zhou et al. | A mixed integer nonlinear programming model for optimal design of natural gas storage surface double-pipe network | |
| CN107575200A (zh) | 一种天然气水合物平行水平井开采的方法 | |
| CN211448630U (zh) | 一种降压和双管注入改性流体开采天然气水合物的装置 | |
| CN107558977A (zh) | 一种排式水平井微波加热开采天然气水合物的方法 | |
| CN207568572U (zh) | 一种模拟注入co2与生石灰粉开采天然气水合物实验装置 |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180313 |
|
| WD01 | Invention patent application deemed withdrawn after publication |