高含硫气藏物质平衡方程通式的新推导

维普资讯 http://www.cqvip.com 断　块　油　气　田　第１５卷第３期　ＦＡＵＩ　ＩＩ—ＢＩ』ＯＣＫ　０ＩＬ＆ＧＡＳ　ＦＩＥＬＤ　２００８年５月　文章编号：１００５－８９０７（２００８）０３—０８１—０３　高含硫气藏物质平衡方程通式的新推导　卞小强　杜建芬　李明军　杨廷玉　邓生辉　（１、西南石油大学油气藏地质及开发工程国家重点实验室，四川成都６１０５００；　２、川气东送建设工程指挥部工程技术部，四川达州６３５０００；３．吐哈油田公司钻采院，鄯善８３８２０２）　摘　要　高含硫气藏是一种非常规的气藏，其物质平衡方程式与常规气藏物质平衡方程式有着一定的差异。高含硫气　藏存在元素硫的析出与沉积，而常规气藏物质平衡方程通式并没有考虑硫沉积现象，所以不能完整地描述和预测高含硫气　藏的开发动态。在考虑硫沉积等多种因素的基础上，从物质的量守恒角度推导出了高含硫气藏的物质平衡方程通式及压降　方程通式。最后，利用新建的高含硫气藏物质平衡方程通式计算出某一含硫气藏的地质储量及压降变化，分析引起误差的　原因，验证了模型的正确性和可靠性，为科学、安全、高效开发高含硫气田提供了一定的理论依据。　关键词　气藏　高含硫硫沉积　物质平衡压降方程　中图分类号：ＴＥ３７５　文献标识码：Ａ　Ｎｅｗ　ｄｅｒｉｖｉｎｇ　ｍｅｔｈｏｄ　ｆｏｒ　ｇｅｎｅｒａｌ　ｅｘｐｒｅｓｓｉｏｎ　ｏｆ　ｍａｔｅｒｉａｌ　ｂａｌａｎｃｅ　ｅｑｕａｔｉｏｎ　ｏｆ　ｇａｓ　ｒｅｓｅｒｖｏｉｒｓ　ｗｉｔｈ　ｈｉｇｌｌ　ｓｕｌｆｕｒ　ｃｏｎｔｅｎｔ　Ｂｉａｎ　Ｘｉａｏｑｉａｎｇ（Ｓｔａｔｅ　Ｋｅｙ　Ｌａｂｏｒａｔｏｒｙ　ｏｆ　Ｏｉｌ　ａｎｄ　Ｇａｓ　Ｒｅｓｅｒｖｏｉｒ　Ｇｅｏｌｏｇｙ　ａｎｄ　Ｅｘｐｌｏｉｔａｔｉｏｎ，Ｓｏｕｔｈｗｅｓｔ　Ｐｅｔｒｏｌｅｕｍ　Ｕｎｉｖｅｒｓｉｔｙ，　Ｃｈｅｎｇｄｕ　６１０５００，Ｃｈｉｎａ），Ｄｕ　Ｊｉａｎｆｅｎ，Ｌｉ　Ｍｉｎｅｊｕｎ，ｅｔ　ａ１．　Ｔｈｅ　ｇａｓ　ｒｅｓｅｒｖｏｉｒ　ｗｉｔｈ　ｈｉｇｈ　ｓｕｌｆｕｒ　ｃｏｎｔｅｎｔ　ｉｓ　ａ　ｔｙｐｅ　ｏｆ　ｕｎｃｏｎｖｅｎｔｉｏｎａｌ　ｇａｓ　ｒｅｓｅｒｖｏｉｒｓ．Ｔｈｅｒｅ　ａｒｅ　ｓｏｍｅ　ｄｉｆｅｒｅｎｃｅｓ　ｂｅｔｗｅｅｎ　ｔｈｅ　ｇａｓ　ｒｅｓｅｒｖｏｉｒ　ｗｉｔｈ　ｈｉｇｈ　ｓｕｌｆｕｒ　ｃｏｎｔｅｎｔ　ａｎｄ　ｔｈａｔ　ｗｉｔｈｏｕｔ　ｓｕｌｆｕｒ．Ｔｈｅｒｅ　ａｒｅ　ｔｈｅ　ｓｅｐａｒａｔｉｏｎ　ａｎｄ　ｄｅｐｏｓｉｔｉｏｎ　ｏｆ　ｓｕｌｆｕｒ　ｉｎ　ｔｈｉｓ　ｔｙｐｅ　ｏｆ　ｇａｓ　ｒｅｓｅｒｖｏｉｒ．Ｂｕｔ　ｔｈｅ　ｓｕｌｆｕｒ　ｄｅｐｏｓｉｔｉｏｎ　ｄｏ　ｎｏｔ　ｂｅ　ｃｏｎｓｉｄｅｒｅｄ　ｉｎ　ｔｈｅ　ｇｅｎａｒａｌ　ｅｘｐｒｅｓｓｉｏｎ　ｏｆ　ｍａｔｅｒｉａｌ　ｂａｌａｎｃｅ　ｅｑｕａｔｉｏｎ　ｏｆ　ｃｏｎｖｅｎｔｉｏｎａｌ　ｇａｓ　ｒｅｓｅｒｖｏｉｒｓ，ｗｈｉｃｈ　ｃａｎ　ｎｏｔ　ｂｅ　ｕｓｅｄ　ｔｏ　ｄｅｓｃｒｉｂｅ　ａｎｄ　ｐｒｅｄｉｃｔ　ｔｈｅ　ｐｒｏｄｕｃｔｉｏｎ　ｐｅｒｆｏｒｍａｎｃｅ　ｏｆ　ｇａｓ　ｒｅｓｅｒｖｏｉｒ　ｏｆ　ｈｉｇｈ　ｓｕｌｆｕｒ　ｃｏｎｔｅｎｔ．　Ｗｈｅｒｅｆｏｒｅ　ｉｔ　ｉｓ　ｏｆ　ａｎ　ｉｍｐｏｒｔａｎｃｅ　ｓｉｇｎｉｉｆｃａｎｃｅ　ｔｏ　ａｃｃｕｒａｔｅｌｙ　ｄｅｒｉｖｅ　ｔｈｅ　ｇｅｎｅｒａｌ　ｅｘｐｒｅｓｓｉｏｎ　ｏｆ　ｍａｔｅｒｉｌａ　ｂａｌａｎｃｅ　ｅｑｕａｔｉｏｎ　ｏｆ　ｈｉｇｈ　ｓｕｌｆｅｒ　ｃｏｎｔｅｎｔ．Ｏｎ　ｔｈｅ　ｂａｓｉｓ　ｏｆ　ｃｏｎｓｉｄｅｒｉｎｇ　ｓｕｌｆｕｒ　ｄｅｐｏｓｉｔｉｏｎ，ｔｈｅ　ｇｅｎｅｒａｌ　ｅｘｐｒｅｓｓｉｏｎ　ｏｆ　ｍａｔｅｒｉｌａ　ｂａｌａｎｃｅ　ｅｑｕａｔｉｏｎ　ｏｆ　ｇａｓ　ｒｅｓｅｒｖｏｉｒｓ　ｗｉｔｈ　ｈｉｇｈ　ｓｕｌｆｕｒ　ｃｏｎｔｅｎｔ　ａｎｄ　ｐｒｅｓｓｕｒｅ　ｄｅｃｌｉｎｅ　ｆｏｒｍｕｌａ　ａｒｅ　ｄｅｒｉｖｅｄ　ｆｒｏｍ　ｍｏｌｅ　ｃｏｎｓｅｒｖａｔｉｏｎ．Ｔｈｅｎ　ｇｅｏｌｏｇｉｃａｌ　ｒｅｓｅｒｖｅｓ　ａｎｄ　ｐｒｅｓｓｕｒｅ　ｄｅｃｌｉｎｅ　ｆｏｒ　ｔｈｅ　ｇａｓ　ｒｅｓｅｒｖｏｉｒｓ　ｗｉｔｈ　ｈｉｇｈ　ｓｕｌｆｕｒ　ｃｏｎｔｅｎｔ　ａｒｅ　ｃａｌｃｕｌａｔｅｄ　ｂｙ　ｍｅａｎｓ　ｏｆ　ｔｈｅ　ｎｅｗ　ｅｑｕａｔｉｏｎｓ　ａｎｄ　ｔｈｅ　ｃａｕｓｅ　ｏｆ　ｅｒｒｏｒ　ｉｓ　ａｎａｌｙｚｅｄ．Ｆｉｎａｌｌｙ，　ｔｈｅ　ａｃｃｕａｒａｃｙ　ａｎｄ　ｒｅｌｉａｂｉｌｉｔｙ　ｏｆ　ｔｈｉｓ　ｍｅｔｈｏｄ　ｈａｖｅ　ｂｅｅｎ　ｖｅｒｉｉｆｄｅ，ｗｈｉｃｈ　ｐｒｏｖｉｄｅｓ　ｓｏｍｅ　ｔｈｅｏｒｅｔｉｃａｌ　ｒｅｆｅｒｅｎｃｅ　ｆｏｒ　ｓｃｉｅｎｔｉｉｆｃ　ｄｅｖｅｌｏｐｍｅｎｔ　ｏｆ　ｔｈｅ　ｇａｓ　ｒｅｓｅｒｖｏｉｒｓ　ｏｆ　ｈｉｇｈ　ｓｕｌｆｅｒ　ｃｏｎｔｅｎｔ．　Ｋｅｙ　ｗｏｒｄｓ：ｇａｓ　ｒｅｓｅｒｖｏｉｒｓ，ｈｉｇｈ　ｓｕｌｆｕｒ　ｃｏｎｔｅｎｔ，ｓｕｌｆｕｒ　ｄｅｐｏｓｉｔｉｏｎ，ｍａｔｅｒｉａｌ　ｂａｌａｎｃｅ，ｐｒｅｓｓｕｒｅ　ｄｅｃｌｉｎｅ　ｆｏｒｍｕｌａ．　近年来，四川东北部发现了越来越多的高含硫天　基本假设如下：①气藏温度在天然气开发过程中不发　然气藏。由于气藏温度在天然气开发过程中基本上不　生变化；②高含硫天然气在初始压力（　）下饱和溶解　发生变化，而物质平衡方法只需要高压物性资料和实　元素硫；⑧气藏温度低于元素硫在气藏压力下的凝固　际生产数据　，因此，物质平衡方法已成为气藏工程　点；（垒）析出的硫为固态且不随气流运移，即析出就沉积　动态分析的方法之一。由于硫的沉积，常规气藏物质平　在孑Ｌ隙中；⑤固态元素硫的密度不随压力变化，即为常　衡方程已经不能完整地描述高含硫气藏复杂的开发动　量；⑥忽略开发过程中储层的组分变化。　态。近年来，有学者从定容气藏出发，研究了高含硫气　当气藏投人开采后，地层压力将不断降低，硫在酸　藏的物质平衡方程　］，但没有实例说明方程的实用性。　气中的溶解度（Ｒ　）也相应地降低，并析出固态硫。沉积　１假设条件　下来的硫将占据一定的孑Ｌ隙空间。因此，高含硫气藏的　高含硫气藏是非常规气藏，随着开发的不断进行，　收稿日期：２００７一ｌ１—３Ｏ；改回日期：２００８一Ｏ１—３１。　气藏压力将不断地下降，溶解在酸气中的元素硫也将　作者简介：卞小强，男，１９７９年生，在读硕士研究生，主要从事高　不断地析出并沉积，这种现象在近井区域尤其明显。　含硫气藏开发研究工作。Ｅ－ｍａｉｈｂｘｑｉａｎｇ３２１０＿８８＠１６３．ｔｏｍ。　８１　维普资讯 http://www.cqvip.com

２００８年５月　断　块　油　气　田　第ｌ５卷第３期　物质平衡方程应该考虑元素硫沉积的影响。　其物质平　衡方程的基本形式为　ｎｓｉ＝ｎｗｐ＋ｎ　＋ｎｓ　（１）　２方程通式推导　２．１　气藏原始地下储量　由气体非理想状态方程可得，气藏原始地下储量　物质的量为　器＝器　（２）　将气藏容积换算到地面标准条件下的体积为　Ｇ＝Ｖｏ／Ｂ　（３）　２．２开发过程中某一压力下的累积采出量　高含硫气藏累积采出井流物物质的量为　㈩　２．３沉积固态元素硫量　１９８２年，Ｃｈｒａｓｔｉｌ提出了固体在超临界气体中的　溶解度公式［１０］：　Ｒｓ＝ｐｓ　ｋｅｘｐ（ａ／Ｔ＋ｂ）　（５）　而Ｒｏｂｅｒｔｓ　Ｅ“　利用Ｃｈｒａｓｔｉｌ的热力学模型，并结合　Ｂｒｕｎｎｅｒ和Ｗｏｌｌ的实验数据ｎ　，建立了硫在酸气中　溶解度的关联式：　Ｒｓ＝ｐｓ４ｅｘｐ（一０．４６６６／Ｔ一４．５７１１）　（６）　又Ｐ　＝　／ＺＲＴ　（７）　所以Ｒｓ＝（　／ＺＲＴ）　ｅｘｐ（－０．４６６６／Ｔ一４．５７１１）（８）　已知：　Ａｍｓ＝Ｖｇ　ｓｉ—　ｇＲｓ　（９）　ｇｓ＝（ＶｇＲｓｌ—ＶｇＲｓ）ｐｓ　（１０）　忽略固体硫的压缩性，元素硫在气藏目前条件下　的体积为　Ｇｓ＝Ｖｓ＝（ＶｓｉＲｓｌ—Ｖ．Ｒｓ）Ｐｓ　（１１）　沉积固态元素硫物质的量为　ｎｓ＝　＝——　Ｍ　：　：—　Ｍ　【ｌＺ　、　‘Ｊ　　２．４气藏剩余气量　当气藏采出一定气量后，气藏的压力从原始地层　压力下降到某个地层压力，气藏压降为△Ｐ（△Ｐ＝　一　Ｐ）。对于高含硫气藏，气藏的孔隙体积会因压力的下　降而收缩，气藏中的束缚水体积会因压力的下降而膨　胀，水体中的水会因压力的下降而侵入气藏，固态硫也　会因压力的下降而析出沉积并且膨胀，以上诸因素均　减小了气藏的容积。　孔隙体积的减小量为　８２　ｃｐ△Ｐ＝　（１３）　其中　ｃｐ＝啬　（１４）　束缚水体积的膨胀量为　ＡＶ　：　Ｃ　△Ｐ：一Ｖｘ，Ｓ，￣Ｃ　，，ＡＰ．　．（１５）　ｌ一．）　由于气藏一般不采用注水开发，因此，气藏被水侵　量占据的孔隙体积即气藏的存水量为　△Ｗ。＝Ｗ。一Ｗ　（１６）　考虑以上几个因素，剩余气量的地下体积为　Ｖｇ＝Ｖｇｉ—ＡＶＰ—ＡＶ　一ＡＷ。一Ｇｓ　（１７）　即　一　ＡＰ一　）＿Ｇｓ（１８）　４－ｃｏ＝竿　告测式（１８　写成　Ｖ　＝Ｖ　ｉ（１一Ｃ　△　一∞）一Ｇｓ　（１９）　从而得到气藏剩余气物质的量为　－，ｇ　（２０）【ｚｕ）　２．５高含硫气藏物质平衡方程通式　联立式（１）一（４）、（１２）、（１９）和（２０），得　虬一　＋　［　ｉ　二　垒　二　２二　！］＋　Ｒ　　。ＺＲＴ　　。ｇ　二　！　（２１）　ｓ　将式（３）代入式（２１）并经整理后，得高含硫气藏物　质平衡方程通式为　一　兹　Ｂｇ　咫　（　１　－Ｃ　ＡＰ　一∞）　奇警＝轰’等　开发过程中，任一压力下天然气的体积系数为　轰‘　ＺＴ　（２４）　将式（２２）－－（２４）变形，可得高含硫气藏物质平衡　方程通式的另一种表达式：　ＣｐＢｓ　－Ｇｓ　＋ｐｓＺＲＴＧｓ　Ｇ：　（２５）５）　维普资讯 http://www.cqvip.com 第１５卷第３期　卞小强，等：高含硫气藏物质平衡方程通式的新推导　２００８年５月　Ｐ一互ＰｓｃＴＧ　ｐ＋　ＧｓｐｓＲ　Ｔ　：　Ｐｉ　ＧＢ￣　ｚｉ　量，ｍ３；Ｗ　为气藏采出的水在地面条件下的体积，ｉｎ。；Ｚ　ｒ，２６　，　、　。　』。　ｌ一乙。‘　ｆ，　』　８　Ｚ　Ｇｓ—Ｇ　ｉ（１一Ｃ　△Ｐ一∞）　、～　为气藏原始条件下天然气的偏差系数；Ｚ为气藏目前　条件下天然气的偏差系数；　为天然气在标准条件下　的偏差系数；Ｔ为气藏的温度，Ｋ；　为标准条件下的　温度，Ｋ；　ｇｉ为天然气在原始压力下的体积系数；　为　天然气目前压力下的体积系数；　为气藏侵入水的体　积系数；Ｃｐ为岩石孔隙体积压缩系数，ＭＰａ～；Ｃｗ为气　３实例分析　以四川Ｉ东北部飞仙关某高含硫干气气藏Ｘ为例　进行研究，气藏的组分及基本参数见表１。　表１　高含硫气藏×的组成及基本参数　组分　浓度／（ｍｏｌ・Ｌ－’）　基本参数　数据　Ｈ２Ｓ　１１．３　／Ｋ　３７６．１　ＣＯ２　６．９８　只／ＭＰａ　４０．０２　ＣＨ－８０．６１　（实测）　１．０１８　４　Ｃ２Ｈ６　０．０４　Ｃ。／ＭＰａ－’　１．２×１０　Ｎ２　１．０６　０　该气藏从２００２年开始投产，至２００６年底地层压　力下降到３０　ＭＰａ，实测的Ｚ＝０．９２８　６，累积产气量为　１６．０２×１０８ｍ。。利用推导的物质平衡方程及文献［１４］　中的数据（ｐｓ＝１　９６０　ｋｇ・ｍ　；Ｍｓ＝２５６　ｋｇ・ｋｍｏｌ　），可以　得到计算的地质储量为８９．７５×１０。ｉｎ。。　综合地质资料、试井资料和测井资料等，该气藏飞　仙关实际地质储量是８６．３１×１０。ｍ３，计算值与实际值　相对误差为３．９８％。可能的原因为：气藏开采后，随着压　力的下降可能有外来水进入地层，但计算时未考虑；气　藏原始条件下，束缚水饱和度很小，在计算时将其忽略。　通过以上分析可以得出，新建的高含硫气藏物质　平衡方程通式可以计算含硫气藏的地质储量及压降的　变化，并分析引起误差的原因，给气藏工程师们提供了　又一新的手段。　４符号说明　／７，　ｉ为气藏原始物质的量，ｋｍｏｌ；凡　为累积产出井　流物物质的量，ｋｍｏｌ；ｒｔ　为气藏剩余物质的量，ｋｍｏｌ；ｒｔｓ　为气藏中沉积元素硫的物质的量，ｋｍｏｌ；只为气藏的原　始地层压力，ＭＰａ；Ｐ为气藏的目前地层压力，ＭＰａ；尸ｓ。　为地面标准条件下的压力；Ｖｇｉ为气藏原始容积，ｉｎ。；Ｖ　为目前条件下气藏的体积，ｉｎ。；Ｖ。为岩石的孔隙体积，　ｍ３；Ｖ　为气藏原始束缚水的体积，ｉｎ。；Ｖｓ为沉积的固　态元素硫在地层条件下的体积，ｉｎ。；Ｇ为气藏地质储　量，ｍ３；Ｇ。为气藏累积采出气量，ｉｎ。；Ｇｓ为沉积的元素　硫折算到标准条件下的体积，ｍ３；ＡＶ。为孔隙体积的减　小量，ｍ３；ＡＶ　为束缚水体积的膨胀量，ｉｎ。；△　。为气　藏的存水量，ｍ３；Ｗｅ为侵入气藏的边底水体积即水侵　藏中束缚水的弹性压缩系数，ＭＰａ　；Ｃ。为气藏容积压　缩系数，ＭＰａ～；尺　为气藏原始条件下硫在酸气中的溶　解度，ｇ・ｉｎ　；尺。为气藏目前条件下硫在酸气中的溶解　度，ｇ・ｉｎ　；Ｐ　为高含硫气藏天然气的密度，ｋｇ・ｉｎ　；Ｐｓ为　固态元素硫的密度，ｋｇ・ｉｎ　；　为天然气的分子量，ｋｇ・　ｋｍｏｌ　；ｓＭ为固态元素硫的分子量，ｋｇ・ｋｍｏｌ　；Ａｒｎｓ为地　层沉积固态元素硫的质量，ｋｇ；　为地层束缚水的饱和　度；∞为气藏存水体积系数；尺为通用体积常数，Ｒ＝　０．００８　３１４　ＭＰａ・ｉｎ。・（Ｋ・ｋｍｏ１）～；　、口、ｂ、ｄ为通过实验确　定的经验系数。　参　考　文　献　［１］　陈元千．异常高压气藏物质平衡方程的推导及应用［Ｊ］　石油学报，　１９８３。４（１）：４５—５３．　［２］　马永祥．对凝析气藏物质平衡方程的研讨［Ｊ］．石油勘探与开发，　１９９７，２４（６）：４５－５０．　［３］陈元千，董宁宇．气藏和凝析气藏物质平衡方程式的新推导［Ｊ］．断　块油气田，１９９９，６（３）：２４—２８．　［４］　陈元千．预测油田天然气可采储量的方法［Ｊ］．断块油气田，１９９９，　６（１）：１８－２５．　［５］李传亮　油藏工程原理［Ｍ］．北京：石油工业出版社，２００５：１２０－１３４　［６］　陈元千，李鋈．现代油藏工程［Ｍ］．北京：石油工业出版社，２００４：　１　１２—１２７．　［７］李中华，鲁章成，张庆生．应用物质平衡法测试抽油井的环空拟液　面［Ｊ］．断块油气田，１９９７，５（１）：３２—３６．　［８］付亚荣，马永忠，童立强．利用ＢＰ神经网络预测高含硫油井的硫　化物应力腐蚀［Ｊ］．断块油气田，２０００，７（１）：４７—５０．　［９］张勇，陈定明．高含硫气藏物质平衡方程的推导［Ｊ］．石油地　质，２００６，２７（６）：７２４—７２５．　［１０］Ｃｈｒａｓｔｉｌ　Ｊ．Ｓｏｌｕｂｉｌｉｔｙ　ｏｆ　ｓｏｌｉｄｓ　ａｎｄ　ｌｉｑｕｉｄｓ　ｉｎ　ｓｕｐｅｒｃｒｉｔｉｅａｌ　ｇａｓｅｓ［Ｊ］．　Ｔｈｅ　Ｊｏｕｒｎａｌ　ｏｆＰｈｙｓｉｃａｌ　Ｃｈｅｍｉｓｔｒｙ，１９８２，８６（１５）：３０１６—３０２１．　［１１］　Ｒｏｂｅａｓ　Ｂ　Ｅ．Ｔｈｅ　ｅｆｆｅｃｔ　ｏｆ　ｓｕｌｆｕｒ　ｄｅｐｏｓｉｔｉｏｎ　ＯＵ　ｇａｓｗｅｌｌ　ｉｎｆｌｏｗ　ｐｅｒｆｏｒｍａｎｃｅ［Ａ］．ＳＰＥ　３６７０７，１９９６．　［１２］Ｂｒｕｎｎｅｒ　Ｅ，Ｗｏｌｌ　Ｗ．Ｓｏｌｕｂｉｌｉｔｙ　ｏｆ　ｓｕｌｆｕｒ　ｉｎ　ｈｙｄｒｏｇｅｎ　ｓｕｌｆｉｄｅ　ａｎｄ　ｓｏｕｒ　￣ｓｅｓ［Ｊ］．ＳＰＥＪ，１９８０，２０（５）：３７７—３８４．　［１３］Ｂｒｕｎｎｅｒ　Ｅ，Ｐｌａｃｅ　Ｊｒ　Ｍ　Ｃ，ｗ０ｕ　Ｗ　Ｈ．Ｓｕｌｆｕｒ　ｓｏｌｕｂｉｌｉｔｙ　ｉｎ　ｓｏｕｒ　ｇａｓ［Ｊ］．　ＪＰＴ，１９８８，４０（１２）：１５８７—１５９２．　［１４］Ｓｈｕａｉ　Ｘ　Ｓ，Ｍｅｉｓｅｎ　Ａ．Ｎｅｗ　ｃｏｒｒｅｌａｔｉｏｎｓ　ｐｒｅｄｉｃｔ　ｐｈｙｓｉｃａｌ　ｐｒｏｐｅｒｔｉｅｓ　ｏｆ　ｅｌｅｍｅｎｔｌａ　ｓｕｌｆｕｒ［Ｊ］．Ｏｉｌ＆ＧａｓＪｏｕｒｎａｌ。１９９５，９３（４２）：５０－５５．　（编辑滕春鸣Ｊ　８３