超深高含硫气藏气—液硫两相渗流实验

超深高含硫气藏气—液硫两相渗流实验
顾少华1,2　石志良1,2　胡向阳1,2　史云清1,2　秦世江3　郭肖3
1.中国石化海相油气田开发重点实验室
2.中国石化石油勘探开发研究院
3.“油气藏地质及开发工程”国家重点实验室•西南石油大学
摘　要　超深高含硫气藏开发过程中会在储层中出现气—液硫同流的现象，其对气井产能的影响目前还缺乏实验数据的验证。

为此，研制了一套适用于高温高压条件下气—液硫两相驱替实时测试的装置，并制订了相应的测试流程，选取四川盆地元坝气田的取样岩心开展气—液硫两相驱替实验，并采用非稳态法计算气、液硫两相相对渗透率，得到气—液硫相对渗透率曲线，进而开展气—液硫两相渗流规律的定量化研究。

结果表明：①气、液硫两相共渗区较窄，当液硫临界饱和度高于40%时，井筒附近的液硫饱和度达到液硫临界流动饱和度，从而阻碍井筒附近气体的流动；②围压的变化会引起气—液硫相对渗透率曲线的变化，当围压增大时，气相相对渗透率及液硫相对渗透率均下降；③随着驱替压差增大，气体流速加快，携硫能力增强，气相相对渗透率及液硫相对渗透率均有所上升。

结论认为，气—液硫两相相对渗透率曲线的获得，实现了对气—液硫两相渗流规律的定量化研究，可用于超深含硫气井的产能评价。

关键词　超深　液硫　酸性气体　相渗曲线　应力敏感　速敏效应　元坝气田　多相渗流
DOI: 10.3787/j.issn.1000-0976.2018.10.010
An experimental study on gas-liquid sulfur two-phase flow in ultradeep high-sulfur
gas reservoirs
Gu Shaohua1,2, Shi Zhiliang1,2, Hu Xiangyang1,2, Shi Yunqing1,2, Qin Shijiang3 & Guo Xiao3
(1. Sinopec Key Laboratory for Marine Oil & Gas Field Development, Beijing 100083, China; 2. Sinopec Exploration & Production Research Institute, Beijing 100083, China; 3. State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation//Southwest Petroleum University, Chengdu, Sichuan 610500, China)
NATUR. GAS IND. VOLUME 38, ISSUE 10, pp.70-75 10/25/2018. (ISSN 1000-0976; In Chinese)
Abstract:During the development of ultradeep high-sulfur gas reservoirs, gas-liquid sulfur co-existence occurs in reservoirs, but its ef-fect on gas well productivity has not been verified by experimental data. In this paper, a set of real-time test devices for gas–liquid sulfur two-phase flow displacement under high temperature and high pressure were developed, and the corresponding test process was formu-lated. Then, the sampling cores of Yuanba Gas Field in the Sichuan Basin were selected for gas-liquid sulfur two-phase displacement experiments. The relative permeability of gas and liquid sulfur was calculated by using non-steady state method and the gas–liquid sulfur relative permeability curve was plotted. Finally, the laws of gas-liquid sulfur two-phase flow were studied quantitatively. And the fol-lowing research results were obtained. First, the co-flow zone of gas and liquid sulfur is relatively narrow. When the critical liquid sulfur saturation is higher than 40%, the liquid sulfur saturation near the wellbore reaches the critical flowing saturation of liquid sulfur, so as to hinder the flowing of gas near the wellbore. Second, the gas–liquid sulfur relative permeability curve varies with the confining pressure. With the increase of confining pressure, both gas relative permeability and liquid sulfur relative permeability decrease. Third, with the increase of displacement pressure difference, the gas flowing speed increases, its sulfur carrying capacity increases and both gas relative permeability and liquid sulfur relative permeability increase slightly. In conclusion, the quantitative research on gas-liquid sulfur two-phase flow is realized based on the establishment of gas-liquid sulfur two-phase relative permeability curve, which can be used to evalu-ate the productivity of ultradeep sulfur gas wells.
Keywords: Ultra deep; Liquid sulfur; Sour gas; Relative permeability curve; Stress sensitivity; Velocity sensitivity effect; Yuanba Gas Field; Multiphase flow
基金项目：国家科技重大专项“高含硫气藏安全高效开发技术”（编号：2016ZX05017）。

作者简介：顾少华，1984年生，高级工程师，博士；主要从事气藏工程方面的研究工作。

地址：（100083）北京市海淀区北四环中路267号。

ORCID: 0000-0001-5203-8113。

E-mail: cc0012@
0　引言
我国含硫天然气资源丰富且拥有数量较多的高含硫气藏，如四川盆地东北地区的普光气田下三叠统飞仙关组及上二叠统长兴组气藏、渡口河气田飞仙关组气藏、罗家寨气田飞仙关组鲕滩气藏[1-3]，以及川北地区的元坝长兴组气藏[4]。

根据前人的研究成果，物理/化学作用后硫可溶解于酸性气体中。

高含硫酸性气藏在开采过程中，随压力和温度变化，酸气溶解硫的能力逐渐下降。

当气体中硫含量达到过饱和时，元素硫会逐渐析出[5-7]。

酸性气藏存在硫析出现象，析出的硫又会因温度的高低分别呈液态和固态，发生在多孔介质中则会影响气体流动。

因此，高含硫气藏在开采过程中存在着多相渗流问题，导致其流动机理异常复杂，对气井生产的影响难以判别[8-10]。

之前酸气中硫析出研究主要针对硫以固态形式析出[11-13]。

这是因为先期开发的酸性气藏普遍埋深较浅，地层温度难以达到硫熔点（115.21 ℃），析出的硫在地层中以固态形式出现。

随着勘探开发技术的不断进步，埋深超过6 000 m 的超深气藏（如四川盆地元坝长兴组气藏等）相继投入开发，地层温度达到硫熔点，导致析出的硫以液态形式存在[14-15]。

随着析出的液硫不断聚集，将对储层渗透性造成伤害，从而降低气井产能，影响气藏的开发效果[16-17]。

目前针对气—液硫两相渗流的研究多从理论方面考虑，而缺乏实验数据的验证。

由于气—液硫同流现象在近期发现，而目前还没有能够完成气—液硫两相驱替实验的现成装置和测试工艺，要进行气—液硫两相驱替实验，还面临一系列技术难点。

首先，为使硫以液态形式流动，液硫
从制备到驱替全过程的温度都需高于硫熔点；其次，液硫易与空气反应，一旦在高温环境下与空气接触会剧烈氧化生成二氧化硫，导致硫质量减少而产生较大测试误差，同时又会对实验人员的健康造成严重危害，因此必须实现全程密封；最后，实验过程中液硫易发生凝固，从而严重损毁设备和管材。

由此可见，气—液硫两相驱替实验的难度和危险并存，导致该技术的研究门槛极高，实验成本高昂，相关研究成果稀少，对气—液硫两相渗流特征缺乏规律性认识。

鉴于上述状况，笔者研制了一套适用于高温高压条件下气—液硫两相驱替实时测试装置，并制订了相应的测试流程，选取元坝气田的取样岩心，开展了气—液硫两相驱替实验，并采用非稳态法对相对渗透率实验数据进行处理，得到了气—液硫相对渗透率曲线，实现了对气—液硫两相渗流规律的定量化研究，可用于超深含硫气井的产能评价。

1　实验装置
实验装置主要由驱替系统（岩心夹持器）、增压系统（增压泵、回压控制器）、储集系统（储氮气罐、储硫中间容器）、围压控制系统（围压泵）、环境模拟系统（恒温箱）、气—液硫分离及收集系统（液硫收集容量瓶）、气—液硫自动计量系统（高精度天平、气体计量器、计算机）、压力监测系统（回压压力表、增压泵/围压泵自带压力表）以及控制软件和数据处理软件组成（图1）。

该装置的主要性能和技术指标如下：驱替压力为25 MPa ，流量为0.000 1～25 mL/min ，围压为80 MPa ，出口回压为20 MPa ，工作温度为150 ℃。

图1　
实验装置示意图
该实验装置针对气—液硫相渗测试中的难点进行了全面考虑。

首先，为解决液硫堵塞问题，在开始实验时将储硫中间容器和岩心夹持器置于恒温箱中，将温度加热至150℃，以保证硫在管线中保持液态形式。

而由于高精度天平难以承受高温，必须将其置于恒温箱外。

为使连接岩心夹持器与液硫收集容量瓶的管线内的硫不发生凝固，可在管线外裹上绝缘电热丝并加热至150℃，以妥善解决因液硫凝固产生的堵塞问题。

其次，由于液硫易与空气反应产生二氧化硫，进而对测试精度及实验人员的健康产生不利影响，在本实验中将实验气体替换为氮气，并持续使用强排气扇以保持实验室通风，并在测试过程中采用自动计量及远程监测技术，以有效减少实验人员与测试装置的接触，防止实验人员二氧化硫中毒的发生。

最后，针对液硫易发生凝固而损耗设备和管材的问题，也采取了以下手段：①针对岩心夹持器内橡胶筒与高温液硫反应导致橡胶筒性能改变的问题，批量采购大量备用密封橡胶筒，每做完一组岩心驱替实验应立即更换橡胶筒；②针对液硫凝固在设备中的问题，准备二硫化碳溶硫剂进行清洗，由于该溶硫剂有恶臭气味，清洗时需佩戴隔绝式防毒面具。

通过对上述系列技术的改进，有效解决了气—液硫实验面临的困难，使气—液硫驱替实验的顺利开展得以实现。

2　实验设计
2.1　测试流程
基于改进的实验装置，需制订相应的实验流程与之配套，参照国家标准GB/T 28912—2012[18]，同时针对某些步骤进行了修改。

2.1.1　岩心的选取与处理
制备直径为2.50 cm的岩心，其长度不小于直径的1.5倍，按照上述标准对岩心样本进行抽提、清洗及烘干，然后测量岩心样本的长度、直径、孔隙度及渗透率。

2.1.2　液硫的制备和实验气体的选取
将硫粉装满储硫中间容器（以下简称中间容器），开动恒温箱对其加热，将硫粉制备成液态硫。

由于硫粉变成液态硫后体积变小，故应选用大型中间容器以便尽可能多装入硫粉，以保证制备出充足的液态硫。

因为液硫在高温高压的空气环境中可自燃，由此选取氮气作为驱替气体，利用增压泵将其注入储氮气罐。

在进行液硫驱替时，将氮气泵入岩心夹持器以驱出岩心中饱和的液硫。

2.1.3　岩心饱和液硫
待装液硫的中间容器冷却后，将其移至恒温箱内，同时将岩心夹持器、回压阀也置于恒温箱内，需使用管线将回压阀接至恒温箱外的液硫收集容量瓶，并利用绝缘电热丝对处于恒温箱外的管线进行加热，防止液硫遇冷凝固而堵塞管路；将恒温箱内部温度升至150℃，回压阀以及相关管路的电加热丝温度亦升至150℃，然后启动驱替泵将中间容器中的液硫泵入岩心中，使岩心充分饱和液态硫。

2.1.4　气—液硫两相物理驱替实验
为保证驱替过程中液硫不会凝固而堵塞管线及岩心，整个驱替管线应保持150℃的高温。

开启驱替泵，使用氮气驱替岩心中的液态硫，直至液硫收集瓶内的液体出口端不再出现液态硫，此时驱替过程结束。

将驱替出的气体和液体进行气液分离，计量气体流量及液硫质量，并记录驱替时岩心两端的进、出口压力。

2.2　测试方案
选取四川盆地元坝气田的取样岩心，采用超低渗气体渗透率仪对所取岩心进行渗透率、孔隙度测试，如表1所示。

表1　岩心测试数据表
岩心编号长度/cm截面积/cm2渗透率/mD孔隙度
13.854.971.653.27%
24.514.988.446.60%
34.614.984.314.24%
44.534.972.913.54%
为研究气—液硫两相相渗曲线特征及其影响因素，针对不同岩心开展气—液硫两相物理驱替实验（各组实验的温度皆全程保持在150℃），实验条件如表2所示。

2.3　非稳态法相对渗透率实验数据处理方法
非稳态法依据Buckley和Leverett[19]提出的一维两相驱替前沿推进理论，认为在两相驱替过程中，气、液硫饱和度是驱替时间和距离的函数，由于气、液硫的相对渗透率随其饱和度变化而变化，气、液硫在岩心某横截面上的流量也随时间不同而变化。

因此，只要在驱替过程中准确记录压力及气、液硫流量的变化，便可计算出两相相对渗透率，具体计算方法如下。

气体通过岩心，进口处压力为p1，出口处压力为p2，气体体积发生变化，必须采用平均体积流量。

表2　气—
液硫两相物理驱替实验条件数据表
岩心编号
围压/MPa
驱替压力/MPa
回压/MPa 4
2018179820
18
17.97
注：各实验温度全程保持在150 ℃
将岩样出口压力（p 2）下测量的流体总产量修正为岩
样平均压力下的数值，即
（1）
式中V i 表示i 时刻的流体（含液硫和气）产量，
mL ；V i －1表示（i －1）时刻的流体产量，mL ；ΔV s i
表示（i －1）时刻到i 时刻的液硫增量，mL ；
p a 表示大气压，MPa ；Δp 表示驱替压差（即p 1和p 2的差值），
MPa ；ΔV g i 表示大气压下测得的某时间间隔对应的气增量，mL 。

将流体总产量修正后，采用式（2）、（3）、（4）、（5）计算非稳态气、液硫相对渗透率，其中气体为驱替相，
液硫为被驱替相。

（2
）
（3）
其中
（4）
式中f s (S g )表示含液硫率，无量纲；
V s (t )表示累积产出液硫量与孔隙体积的比值，无量纲；
V (t )表示液硫和气的总产出量与孔隙体积的比值，无量纲；
K rs 表示液硫相对渗透率，无量纲；
I 表示相对注入能力（即流动能力比），无量纲；
Q (t )表示t 时刻岩心出口端液硫流量，cm 3/s ；
Q s 表示初始时刻岩心出口端液硫流量，cm 3/s ；
Δp 0表示初始时刻驱替压差，MPa ；Δp (t )表示t 时刻驱替压差，MPa ；K rg 表示气相相对
渗透率，无量纲；
μg 表示气体黏度，mPa ·s ；μs 表示液硫黏度，mPa ·s 。

3　测试结果
3.1　气—液硫相渗曲线形态特征
1号岩心的气—液硫相渗曲线说明液硫产生后气相相对渗透率大幅下降，而液硫相对渗透率缓慢上升，液硫最大相对渗透率仅0.17（图2）。

气、液硫两相共渗区较窄，液硫临界流动饱和度较大，共渗区液硫饱和度介于56%～87%。

由于液硫临界流动饱和度较大，因此含硫气藏内绝大多数区域液硫产生后很难流动。

由于在井筒附近，压降最明显，液硫析出量最多，同时气流汇聚至井筒，大量的液硫被携带至井筒附近并聚集，造成井筒附近的液硫饱和度容易达到液硫临界流动饱和度，进而阻碍井筒附近气体的流动。

图2　1号岩心气—液硫相对渗透率曲线图（围压为20 MPa ）
3.2　不同围压下气—液硫相渗曲线
通过改变围压，测得1号、2号岩心的相渗曲线，
随着围压的增大，气相相对渗透率及液硫相对渗透率均下降（图3）。

由于围压升高，增加了作用在岩石骨架上的有效应力，同时由于气体具有较强可压缩性，使岩石孔隙及吼道变小。

根据微管流动理论[20-21]，由于液体边界层的存在，使流体（含气相和液相）可通过的空间变小，从而导致气相及液硫的相对渗透率下降。

图3　1号、2号岩心在不同围压下气—液硫相对渗透率曲线图
3.3　驱替压差对气—液硫相渗曲线的影响
由3号、4号岩心的气—液硫相对渗透率曲线，可见，随着驱替压差增大，气相相对渗透率及液硫相对渗透率均有一定程度地上升（图4）。

这是由于驱替压差增大，气体流速增加，部分孔道壁面处被束缚的液硫受到气体的携带而脱离孔道壁面，使得液硫的相对渗透率上升。

该现象称为Henderson 效应[21]，是Henderson 等针对凝析气通过驱替实验研究气—凝析液相对渗透率曲线时发现的，可见气—液硫两相渗流过程中同样存在Henderson 效应，且气体流速越快，携液硫能力越强。

随着孔隙中的液硫被部分携
图4　3号、4号岩心在不同驱替压差下气—
液硫相对渗透率曲线图
带出，液硫量减少，使气—液硫两相共渗时流动能力均有所提升。

因此，在进行气—液硫渗流研究时，需适当考虑气体流速的影响。

4　结论和认识
1）气、液硫两相共渗区较窄，液硫临界饱和度高于40%，井筒附近的液硫饱和度容易达到液硫临界流动饱和度，从而阻碍井筒附近气体的流动。

2）围压的变化会引起气—液硫相对渗透率曲线
的变化。

当围压增大时，气相相对渗透率及液硫相
对渗透率均下降。

3）气—液硫两相渗流过程中同样存在Hender-son效应，随着驱替压差增大，气体流速加快，携硫能力增强，气相相对渗透率及液硫相对渗透率均有所上升，从而使气—液硫两相共渗时流动能力均有所提升。

参　考　文　献
[ 1 ] 段金宝. 普光与元坝礁滩气田天然气成藏特征对比[J]. 西南石油大学学报(自然科学版), 2016, 38(4): 9-18.
Duan Jinbao. The comparative study of natural gas accumula-tion characteristics between Puguang and Yuanba reef-bank gas field[J]. Journal of Southwest Petroleum University (Science & Technology), 2016, 38(4): 9-18.
[ 2 ] 曾大乾, 彭鑫岭, 刘志远, 宿亚仙. 普光气田礁滩相储层表征方法[J]. 天然气工业, 2011, 31(3): 9-13.
Zeng Daqian, Peng Xinling, Liu Zhiyuan & Su Yaxian. Charac-terization methods of reef-beach facies reservoirs in the Puguang Gas Field[J]. Natural Gas Industry, 2011, 31(3): 9-13.
[ 3 ] 王海, 林然, 张晨阳, 黄炳光, 缪长生. 串珠状缝洞型碳酸盐岩储层压力变化特征研究[J]. 西南石油大学学报(自然科学版), 2017, 39(1): 124-132.
Wang Hai, Lin Ran, Zhang Chenyang, Huang Bingguang & Miao Changsheng. Pressure variation characteristics in bead-shaped fractured vuggy carbonate gas reservoirs[J]. Journal of Southwest Petroleum University (Science & Technology), 2017, 39(1): 124-132.
[ 4 ] 武恒志, 李忠平, 柯光明. 元坝气田长兴组生物礁气藏特征及开发对策[J]. 天然气工业, 2016, 36(9): 11-19.
Wu Hengzhi, Li Zhongping & Ke Guangming. Characteristics of the Changxing Fm biohermal gas reservoir in the Yuanba Gas-field, Sichuan Basin and development countermeasures[J]. Natu-ral Gas Industry, 2016, 36(9): 11-19.
[ 5 ] Sun CY & Chen GJ. Experimental and modeling studies on sulfur solubility in sour gas[J]. Fluid Phase Equilibria, 2003, 214(2): 187-195.
[ 6 ] Robert BE. The effect of sulfur deposition on gas well inflow per-formance[J]. SPE Reservoir Engineering, 1997, 12(2): 118-123. [ 7 ] Gu Minxin, Li Qun, Zhou Shanyan, Chen Weidong & Guo Tian-min. Experimental and modeling studies on the phase behavior of high H2S-content natural gas mixtures[J]. Fluid Phase Equilibria, 1993, 82(82): 173-182.
[ 8 ] Hyne JB. Study aids prediction of sulfur deposition in sour-gas wells[J]. Oil and Gas Journal, 1968, 25(11): 107-113.
[ 9 ] Al-awadhy F, Kocabas I, Abou-Kassem JH & Islam MR. Experi-mental and numerical modeling of sulfur plugging in a carbonate oil reservoir[J]. Energy Sources, 2005, 27 (1-2): 3-18.
[10] 胡书勇, 张烈辉. 高含硫气藏硫沉积的研究进展[J]. 钻采工艺,
2009, 32(1): 71-74.
Hu Shuyong & Zhang Liehui. Research on elemental sulfur dep-osition in gas reservoir with high H2S[J]. Drilling & Production Technology, 2009, 32(1): 71-74.
[11] 曾平, 赵金洲, 周洪彬. 高含硫气藏元素硫沉积对储集层的伤
害[J]. 石油勘探与开发, 2005, 32(6): 113-115.
Zeng Ping, Zhao Jinzhou & Zhou Hongbin. Study on gas per-colation behavior for gas reservoirs with high sulfur content[J].
Petroleum Exploration and Development, 2005, 32(6): 113-115.
[12] 杨学锋, 黄先平, 杜志敏, 李传亮. 考虑非平衡过程元素硫
沉积对高含硫气藏储层伤害研究[J]. 大庆石油地质与开发, 2007, 26(6): 67-70.
Yang Xuefeng, Huang Xianping, Du Zhimin & Li Chuanliang.
Reservoir damage caused by elemental sulfur deposition to high sulfur content gas reservoirs with the consideration of non-equi-librium process[J]. Petroleum Geology & Oilfield Development in Daqing, 2007, 26(6): 67-70.
[13] 翟广福, 李晓平. 气藏中硫的沉积机理及其对生产的影响[J].
新疆石油地质,2005,26(3): 301-303.
Zhai Guangfu & Li Xiaoping. Study on the sulfide precipitation mechanism and its affect on production and exploitation of gas reservoir[J]. Xinjiang Petroleum Geology, 2005, 26(3): 301-303.
[14] 杜志敏. 高含硫气藏流体相态实验和硫沉积数值模拟[J]. 天
然气工业, 2008, 28(4): 78-81.
Du Zhimin. High sulfur gas reservoir fluid phase state experiment and numerical simulation of sulfur deposition[J] Natural Gas In-dustry, 2008, 28(4): 78-81.
[15] 张勇, 杜志敏, 杨学锋, 游建国. 高含硫气藏气—固两相多组
分数值模型[J]. 天然气工业, 2006, 26(8): 93-95.
Zhang Yong, Du Zhimin, Yang Xuefeng & You Jianguo. Gas-sol-id two-phase multicomponent numerical model of acid gas reser-voirs[J]. Natural Gas Industry, 2006, 26(8): 93-95.
[16] Mahmoud M. Effect of elemental-sulfur deposition on the rock
petrophysical properties in sour-gas reservoirs[J]. SPE Journal, 2014, 19(4): 703-715.
[17] Guo Xiao & Wang Qing. A new prediction model of elemental
sulfur solubility in sour gas mixtures[J]. Journal of Natural Gas Science & Engineering, 2016, 31: 98-107.
[18] 中华人民共和国国家质量监督检验检疫总局, 中华人民共和
国国家标准化管理委员会. 岩石中两相流体相对渗透率测定方法: GB/T 28912—2012[S]. 北京: 中国标准出版社, 2012.
General Administration of Quality Supervision, Inspection and Quarantine of the PRC & Standarization Administration of the PRC. Test method for two phase relative permeability in rock: GB/T 28912-2012[S]. Beijing: Standards Press of China, 2012. [19] Buckley SE & Leverett MC. Mechanism of fluid displacement in
sands[J]. Transactions of the AIME, 1942, 146(1): 107-116. [20] 朱维耀, 田英爱, 于明旭, 韩宏彦, 张雪龄. 微圆管中流体的
微观流动机制[J]. 科技导报, 2014, 32(27):23-27.
Zhu Weiyao, Tian Ying'ai, Yu Mingxu, Han Hongyan & Zhang Xueling. Mechanism of microscopic fluid flow in microtubes[J].
Science & Technology Review, 2014, 32(27): 23-27.
[21] Henderson GD, Danesh A, Tehrani DH & Al-Kharusi B. The
relative significance of positive coupling and inertial effects on gas condensate relative permeabilities at high velocity[C]//SPE Annual Technical Conference and Exhibition, 1-4 October 2000, Dallas, Texas, USA. DOI: https:///10.2118/62933-MS.
（修改回稿日期　2018-08-02　编　辑　孔　玲）。