静态酸化实验装置

『壹』 　沉积盆地流体-岩石相互作用研究方法和手段

在盆地沉积物埋藏后所经历的成岩过程中，会发生复杂的微生物、有机质、水、岩之间的相互作用过程。若烃类发生侵位，还涉及烃类参与的反应。传统上往往将它们单独地分别研究。流体-岩石相互作用研究力图将烃源岩、储集岩矿物和孔隙流体（油、气、水）及其中的微生物作为一个完整的地球化学系统来研究其相互作用，这就要求进行沉积学、水文地质学、同位素地球化学、微生物学等多学科交叉研究，将地质观察、实验模拟、计算机模拟结合在一起，解决一些单一学科的问题。下面介绍实验地球化学测试、实验室模拟、热力学理论计算等方面的研究方法。计算机软件模拟将专门分章讨论。

一、实验地球化学测试

沉积盆地流体-岩石相互作用研究需要对储层中油、气、水、岩进行全面的分析。所分析的项目及数量取决于研究的内容和目标，不能一概而论。

1.分析测试内容

岩石分析岩石的矿物成分、化学组成和储层物性；碳酸盐胶结物的碳、氧、锶同位素组成；硫酸盐和硫化物的产状、矿物习性、硫同位素组成；粘土矿物的X射线衍射分析和氧同位素分析。

流体包裹体分析流体包裹体包括液相和气相包裹体，液相又包括水相和烃类。均一化温度是各类流体包裹体常分析的内容，用以确定胶结物形成时期、油气注入时间。对于水相包裹体，需测定Na、K、Ca、Cl组成及盐度，用激光拉曼光谱测定溶解的CH₄、H₂S、CO₂气体质量分数，H₂S硫同位素和CO₂的碳同位素。对烃类包裹体则可进行全烃色谱分析，以确定是否发生蚀变。

油田水分析用毛细管等速电泳或高效液相色谱（HPLC）分析有机酸中甲酸、乙酸、丙酸、丁酸、苯甲酸等的浓度及总量。利用等离子发射光谱（ICP）分析微量元素K、Sr、Mn、Al、Fe、Zn、B、Li、Cs、Cd等。用钼-硅法分析其中二氧化硅的含量。用质谱仪分析碳、氢、氧、硫、锶、硼的同位素组成。

烃类分析分析稠油或沥青的物性和族组成、气相色谱特征、生物标志物和硫同位素，并与正常原油对比，以研究其成因机制。分析伴生气的气体组分和碳、硫同位素。

2.分析测试技术

国内众多的实验室已建立起了成熟的方法，来分析上述岩石学、流体包裹体及烃类分析的项目。唯粘土矿物（高岭石、蒙脱石和伊利石）的氧同位素分析国内尚未开展，但国外已有报道。油田水有机组分、微量元素及同位素分析，尚未为人熟知，有必要简要介绍。

1）有机酸分析技术

（1）等速电泳法（ITP）该法采用在中空的毛细管内进行恒流电泳的独特的分离分析方法。油水样经水相蒸发预处理，除去大量无机盐类后，即可直接进样进行有机酸分离。所用仪器为瑞典LKB-2127等速电泳仪及岛津IP-2A型等速电泳仪，检测器为电导检测器、紫外检测器及电位梯度检测器，配以200mm×0.5mm聚四氟乙烯毛细管（LKB-2127）及50cm×1mm、100cm×0.5mm两级聚四氟乙烯毛细管（IP-2A）。采用电解质溶液及尾随电解质溶液分别为组氨酸盐＋组氨酸溶液及2-N吗啉代乙磺酸溶液，或为HCl＋β-丙氨酸溶液及正己酸溶液。水相蒸发处理过程为：取水样低温蒸发，调至酸性，然后以丙酮洗涤过滤，再调节至碱性，浓缩定容。方法的回收率及相对标准偏差分别为96%～105%和2.4%～7.6%。

（2）区带电泳法（CZE）由于油田水中Cl^-干扰测定结果，等速电泳法需对样品进行水相蒸发预处理，采用区带电泳法则避免了上述预处理。所用仪器为惠普HP3PCE高效毛细管电泳仪，毛细管为50cm×50μm内径熔融石英毛细管（有效长度48.5cm），检测器为二极管阵列检测器。电解质体系为：①邻苯二甲酸氢钾＋十六烷基三甲基溴化铵，pH=6.0；②3,5-二硝基苯甲酸＋十六烷基三甲基溴化铵＋5%甲醇，pH=9.0。检测波长为254nm及210nm，间接检测，压力样进，油田水样过滤后，即可直接进样进行有机酸分离。方法的相对标准偏差为1.1%～3.5%。

（3）毛细管气相色谱法（GC）利用AT1000大口径极性毛细管柱，对油田水中C₂—C₅一元羧酸进行分离分析。对油田水以水相蒸发除去大量无机盐类后，经浓缩再直接进样，无需酸化和萃取。方法回收率和相对标准偏差分别为79.6%～100%及1.9%～6.4%。

2）同步辐射X射线荧光分析

利用北京正负电子对撞机国家实验室同步辐射装置，在专用模式下进行工作。实验测试时，样品受同步辐射X射线激发，发生电离，被电离的原子产生次级特征X射线。每种元素有其固有的特征X射线能量及相应的特征波长，用Si（Li）探测器测定这些特征X射线的能量可判断元素的类别；根据测得的待测元素的特征X射线荧光计数与相同实验条件下标样所测的该元素的计数比较，可得出元素的含量。

由于同步辐射具有高亮度、高准直、线偏振及宽频可调等优异特性，因而用于样品的微量元素分析时灵敏度高，对制样要求简单，可在保持样品原始状态下进行测定，并能在相同的实验条件下同时测定一个油田水样品中的20多种微量元素，检测下限可达10^-6量级。

3）δD、δ¹⁸O、δ³⁴S和⁸⁷Sr/⁸⁶Sr的测定

δD的测试采用的是高纯锌（Zn）还原法，即将2μL水样在390℃下经过锌还原出氢气，然后用MAT251质谱仪测定氢气的D/H值。δ¹⁸O的测量采用CO₂-H₂O平衡法，即将一定量的CO₂高纯钢瓶二氧化碳与2mL水样平衡，用MAT251型质谱仪测定平衡后CO₂的¹⁸O/¹⁶O。δD、δ¹⁸O测试结果均以SMOW（标准平均大洋水）为标准给出，其标准偏差分别为1‰～2‰和0.20‰～0.30‰。

δ³⁴S硫化物硫同位素分析方法是，将硫化物与一定比例CuO混合，在1100℃下真空燃烧制备纯的SO₂气体。硫酸盐、自然硫或岩石中微量硫，均采用埃斯卡试剂处理，转化为氧同位素基本纯的硫酸钡。制样时，称取一定量的BaSO₄、V₂O₅、SiO₂（比例为1∶3.5∶3.5），混合均匀后放入瓷瓶内，并在其上覆盖一层铜丝，在980℃的真空热解下，制备纯的SO₂气，然后用MAT251型质谱计测定³⁴S/³²S值。δ³⁴S值以CDT（为迪亚布洛峡谷陨石中的陨硫铁）标准给出。其标准偏差为±0.10‰～0.30‰。

⁸⁷Sr/⁸⁶Sr测定方法是，取一定量地层水，用超纯HCl酸化，经过标准离子交换技术分离后，在MTA261型多接收器质谱仪上进行测定。溶解碳酸盐全岩、胶结物是用超纯的HCl，溶解页岩采用超纯HF和HClO₄试剂。分析精度0.00003～0.00007。其中，地层水样来自中途测试或完井测试。但是，这类样品不可能有足够的采样覆盖面，尤其在井内更是如此。最有效的弥补方法是使用岩心样品，这就涉及岩心的保护及其水的离心分离。在应用了低浸染取心技术（即最大限度地减少泥浆对水的污染）以后，这种方法非常实用。还有一种是RSA法，即残余盐分析法。在实验室中用超纯水浸滤未经保护的常规岩心，以溶解孔隙中的盐。这种盐是岩心在储藏期间从蒸发的地层水中沉淀出来的。由于不可能浸滤出100%的盐类物质，所以浸滤出的盐不保留原始地层水总体化学性质。但是通过对RSA法的有效性严格检验后，发现锶同位素⁸⁷Sr/⁸⁶Sr比值却不受影响。在取样过程中必须避免在岩心边缘、裂隙面和含有高渗透性岩石的部位取样，筛去具有污染特征的数据（取决于渗透性与⁸⁷Sr/⁸⁶Sr之间的关系），还要沿一些岩样的半径方向测定RSA法的数据特征，以此来校验岩心中央未被污染水的稳定比值。与多种钻井泥浆渗透液相比，地层水中的高Sr含量意味着水中⁸⁷Sr/⁸⁶Sr比对污染作用相对地不太敏感。比较而言，地层水的⁸⁷Sr/⁸⁶Sr比值为0.705～0.730，砂岩中矿物的⁸⁷Sr/⁸⁶Sr比值变化范围更大：斜长石或碳酸盐小于0.710，钾长石大于0.730，而云母大于0.800。可见，用RSA法可以将油田水⁸⁷Sr/⁸⁶Sr比值十分精确地测定出来（Smalley,1987）。

二、实验室模拟

模拟实验是在实验室中通过控制实验条件来模拟自然条件下流体-岩石相互作用的过程。模拟实验包括动力学和热力学两种模拟方法。中国地质科学院张荣华研究员一直在模拟研究开放体系中方解石、萤石等矿物-水的反应动力学。而沉积盆地水-岩反应更常发生在半封闭-半开放体系中。模拟的内容包括：有机酸、CO₂的生成；有机组分（原油、有机酸等）参与的水-岩相互作用；金属有机配位化合物稳定性的实验测量等。常用的模拟实验方法是流动或动态实验装置（Barth等，1988；杨俊杰等，1995）。该方法是将反应溶液从一端注入，并在控制的温度、流速下与反应容器中涂有环氧树脂的岩心发生作用。反应溶液可以是各种合成地层水，可含有机酸或原油。在不同的持续时间里从另一端收集反应后的溶液，观测水化学的变化。另一方法采用间歇反应器（静态装置），反应容器可用不锈钢、钛制成。采集并分析经不同时间反应后的溶液，对比实验前后岩石的显微特征、物性或原油性质的变化，以达到模拟研究流体-岩石相互作用的目的。

三、热力学理论计算

热力学理论计算方法是运用热力学定律，对地球化学反应和过程进行理论计算来推断和解释各种地球化学现象（梅廉夫等，1994），可为实验结果的延拓、解释和检验提供理论依据。倪师军等（1993）根据流体包裹体温度、压力、成分及E_h-pH值，计算了成岩流体与矿物相互作用的趋势。而自由能更广泛应用于化学反应趋势的预测上。McBride（1987）、罗明高（1995）以反应的自由能模拟计算了成岩作用的序列；Meshri（1990）对比研究了碳酸和有机酸的热力学反应能力，计算了碳酸盐矿物方解石和铝硅酸盐矿物长石的溶解趋势和向粘土矿物转化趋势。Giles（1990）利用质量传递方程研究了矿物溶解-沉淀、离子迁移能力对次生孔隙和总孔隙度变化的影响。可见，热力学理论计算已用于地质现象的解释和预测上，是计算机软件模拟的基础。但相对而言，考虑的因素较为单一。

『贰』 污水处理絮凝沉淀时，PAM，PAC的用量怎么衡算

污水处理絮凝沉淀时，通过以下几个公式进行运算：

1、加药量mg/L=加药质量/处理水量/配药浓度

2、处理水量投加药量=处理水量m3/h*加药量g/m3

3、干泥量=处理水量*【（1-污泥含水率）/（1-泥饼含水率）】

4、每吨干泥的药剂消耗g/m3=加药量/干泥量

(2)静态酸化实验装置扩展阅读

用途

一、PAM用途

1、用于污泥脱水根据污泥性质可选用本产品的相应型号，可有效在污泥进入压滤之前进行污泥脱水，脱水时，产生絮团大，不粘滤布，压滤时不散，流泥饼较厚，脱水效率高，泥饼含水率在80%以下。

2、用于生活污水和有机废水的处理，本产品在配性或碱性介质中均呈现阳电性，这样对污水中悬浮颗粒带阴电荷的污水进行絮凝沉淀，澄清很有效。

3、用于以江河水作水源的自来水的处理絮凝剂，用量少，效果好，成本低，特别是和无机絮凝剂复合使用效果更好，它将成为治长江、黄河及其它流域的自来水厂的高效絮凝剂。

4、造纸用增强剂及其它助剂。提高填料、颜料等存留率、纸张的强度。

5、用于油田助剂，如粘土防膨剂，油田酸化用稠化剂。

6、用于纺织上浆剂、浆液性能稳定、落浆少、织物断头率低、布面光洁。

二、PAC用途：

1、城市给排水净化：河流水、水库水、地下水。

2、工业给水净化。

3、城市污水处理。

4、工业废水和废渣中有用物质的回收、促进洗煤废水中煤粉的沉降、淀粉制造业中淀粉的回收。

5、各种工业废水处理：印染废水、皮革废水、含氟废水、重金属废水、含油废水、造纸废水、洗煤废水、矿山废水、酿造废水、冶金废水、肉类加工废水、污水处理。

6、造纸施胶。

7、糖液精制。

8、铸造成型。

9、布匹防皱。

10、催化剂载体。

11、医药精制

12、水泥速凝。

13、化妆品原料。

『叁』 国内外研究现状

1.2.1 元素硫溶解度及沉积运移实验研究现状

（1）元素硫溶解度研究现状

对高含硫天然气中元素硫溶解度的认识是该类气藏开发过程中重要的环节之一。国内外对该问题进行了深入的研究。硫溶解度的研究主要包括实验和理论两个部分，以下为实验部分。

1960年，Kennedy^［7］等人研究了硫在不同含量的CH₄、CO₂和H₂S三种气体中的平衡溶解问题。并且首次说明了硫的溶解性能与气体压力、温度和组分有关。在一定温度压力的条件下，其溶解能力大小依次为H₂S、CO₂、CH₄。

1971年，Roof^[8]通过实验研究了低温低压条件下硫在硫化氢气体中的溶解度（压力6.8～30.6MPa，温度43.3℃～114℃）。

1976年，为了更好地研究深层气藏的高温高压条件下硫在酸性气体中的溶解度，Swift^［9］进行了溶解度实验研究（压力34.5 MPa～138 MPa，温度121℃～204℃）。

1980年和1988年，E.Brunner^{［10～11］}等人将Kennedy等人研究进行推广（压力6.6MPa ～155MPa，温度116℃～213℃），研究了硫在不同比例的CO₂、H₂S、C₁～C₄的14个合成酸性气体混合物中的溶解度。

1992年和1993年，P.M.Davis^［12］等人将E.Brunner等人的研究成果进行了深入研究（压力7 MPa～55MPa，温度60℃～150℃），将硫在简单多组分中的溶解扩展到实际的酸气组分中。

1993年，谷明星^{［13～14］}等人建立了静态法测定难挥发溶质（固体或液体）在超临界、近临界流体中溶解度的实验装置，针对硫化氢大于50％的富含H₂S酸性流体溶解度进行了测试。

2003年，C.Y.Sun^［15］在谷明星实验研究的基础上，在室内利用静态实验测试装置完成了元素硫在7个高含硫混合气体（H₂S CO₂、CH₄）中溶解度测定，并建立了能预测和关联硫在高含硫天然气中溶解度的气固热力学模型。

2005年，曾平^{［16～17］}对元素硫在天然气中的溶解度进行了实验研究，并对其机理进行了说明，分析了不同组分对元素硫溶解度的影响，提出混合物中含碳原子数目较多的烃类组分对硫溶解度有着重要的影响。

2009年，杨学锋^［18］通过自主设计的元素硫溶解度实验设备，针对Chrasnti^［19］和Roberts^［20］常系数模型进行了关联性研究，发现Chrasnti l溶解度计算模型更加科学可靠； 而Roberts溶解度模型，由于是根据有限特定的几组数据拟合得到，具有一定的局限性。

由于硫在含硫混合气中溶解度测试具有一定的危险性，故为了更好的得到硫在含硫混合气中的溶解度，国内外学者在理论模型方面也做了很多深入的研究。

1980年和1983年，J.B.Hyne^{［21～22］}等人研究发现随着温度压力的升高，元素硫和硫化氢会生成多硫化氢。反之，随着温度压力的降低，多硫化氢又会分解成为元素硫和硫化氢，从而导致硫沉积。

1982年，Chrastil^［19］基于理想溶液理论，提出了一个简化的热力学方程来计算硫的溶解度。该经验公式已经广泛用于超临界流体溶质溶解度的计算。

1989年，R.A.Tamxej^［23］等人在对大量实验数据进行拟合的基础上，得到了元素硫在含硫气体中溶解度的预测模型。

1997年，E.Bruce^［20］等人利用Brunner^［10］和Woll的实验数据，对Chrastil经验公式进行了回归拟合，建立了元素硫在酸性气体中的溶解度经验公式，该公式考虑了温度、压力和气体组分对元素硫溶解度的影响，因为方便应用，故一直被用于预测元素硫在含硫天然气中的溶解度。

1998年，Kunal Karan^［24］等人建立一个热动态模型，可用于预测酸气混合气体中硫溶解度，并利用该模型计算了元素硫在硫化氢和高含硫气体混合物中的溶解度。

2003年，C.Y.Sun^［15］等人采用与谷明星类似的方法，建立了能够预测和关联元素硫在高含硫天然混合气中溶解度的气固热力学模型。

2006年，杨学锋^［25］引入了超临界流体的压缩气体模型，建立了元素硫和高含硫天然气达到气固相平衡时定量计算元素硫溶解度的关联和预测模型。

（2）元素硫沉积运移实验研究现状

随着温度压力的降低，元素硫会从含硫天然气中析出，部分硫颗粒将会沉降，部分硫颗粒则会随储层流体运移。

目前，元素硫沉积实验主要集中在油藏方面^{［27～29］}，由于硫化氢的剧毒性，开展高含硫元素硫沉积储层伤害的实验极少。

2000年，Jamal H.Abou-Kassem^［30］利用氮气携带升华的元素硫进入碳酸岩岩心，观察和测定了元素硫对岩心的伤害。提出了一种简易的方法来模拟实际高含硫气藏元素硫对储层的伤害，但由于元素硫升华的温度极高，对其实验及数据的可行性值得深入探讨。

2008年，西南油气田分公司勘探开发研究院^{［31 ］}自主研制了模拟实际储层高温高压的条件下，元素硫沉积对储层伤害驱替实验仪器，完成了不同初始压力、温度下元素硫对天然碳酸盐岩岩心渗透率和孔隙度的伤害。

1.2.2 含硫气藏储层改造铁离子伤害研究现状

储层改造作为低渗透油气藏重要的增产措施已经得到了广泛的认可，目前含硫气藏也通常进行酸压改造增产作业。由于含硫气藏涉及元素硫沉积和酸性气体等因素，对其储层改造必要性的探讨还存在空白。

考虑到元素硫沉积和酸性气体的影响，含硫气藏储层改造的核心就是控硫控铁。在处理含硫化氢气井的储层改造问题上，国内外主要集中在控制铁沉积上^{［32～37］}。在酸压作业中，对于控制铁离子沉淀，通常有三种方法：

一是对主体工作液进行研究，采用弱酸体系来控制残酸液的pH值，使得残酸pH值处于一个相对较低的位置，以便于抑制残液中析出含铁的硫化物。

二是采用铁离子络合剂。由于络合剂对高价的金属离子具有较强的亲和力，从而使得溶液中铁离子浓度低于析出沉淀的浓度，从而抑制铁离子沉积的产生。

三是采用还原剂，将溶液中的三价铁离子还原成为二价铁离子，从而达到避免沉淀析出的目的。

2004年，陈红军^［38］等人对于含硫化氢气井酸化过程中，硫化铁沉淀预测及抑制剂研究进行了详细的调研和研究，并提出了一套适应含硫气井酸压作业且与之匹配的添加剂，优化了酸液体系的整体性能，其具体表现为铁离子稳定剂、硫化氢吸收剂和控硫剂。

2007年，Jairo Leal^［39］等人在分析了在对解除硫化铁沉积过程中可能会出现的问题，提出了一系统有序的方法来对硫化铁沉淀进行移除。

2009年，Tao Chen^［40］等人建立了一套新的硫化铁测试方法来评价硫化铁抑制剂的性能。在此基础上，研制了一种新的抑制剂并对硫化铁抑制剂机理进行了说明。

1.2.3 元素硫沉积对储层伤害研究现状

为了研究地层条件下元素硫沉积对储层的伤害，国内外学者分别建立了考虑元素硫伤害的含硫气藏伤害模型，分析元素硫沉积对储层参数及产能的影响。

1966年，C.H.Kuo^［41］建立流体流动数学模型，该模型能够描述多孔介质中固相沉积。该模型假设初始状态含硫天然气饱和溶解元素硫。

1972年，C.H.Kuo^［42］将硫沉积模型引入，在黑油模型的基础上，建立元素硫沉积的储层伤害数学模型，该模型考虑了硫溶解度的变化和硫沉积对渗透率和孔隙度的影响。该模型能够模拟均质气藏一维径向流动情况下，采气速度、井距和井筒半径对硫沉积的影响。

1980年，J.B.Hyne^［21］等人通过统计学原理，分析了100多口含硫气井的元素硫沉积问题，分析了混合物中不同碳原子数、CO₂、硫化氢含量对元素硫沉积的影响。

1997年，E.Boberts^［20］在等温稳态理想流动的条件下，研究了酸性气井中元素硫沉积对流人动态的影响，建立了考虑元素硫沉积储层伤害模型，分析了不同时间，不同径向距离处元素硫饱和度的分布。发现硫的聚集速度与径向距离平方成反比，径向距离小，元素硫沉降距离的越快。同时还考虑表皮的影响，表皮越小，硫的聚集速度越小，但该模型假设元素硫析出就地沉降，没有考虑元素硫运移。

1997年，王琛^［43］在Roberts建立的理论基础上，研究了硫沉积对气井产能的影响及各因素对硫沉积的影响。

2001年，Faruk Civan^［44］将延迟效应引入到元素硫沉积里面，考虑元素硫动态沉积，即元素硫析出后不会就地沉降，而是运移一段时间或位移后再沉降。但并没有说明元素硫何时沉降，运移多长时间和位移。

2002年，Nicholas Hands^［45］等建立了天然裂缝性含硫气藏硫沉积预测解析模型，该模型考虑了温度和近井地带的气流临界流速的影响，对元素硫在近井地带的分布进行了分区和详细地研究，并给出了相应的井底除硫时间，但对于元素硫颗粒临界流速计算并没有给出具体计算方法。

2004年，杨满平^［46］考虑非达西渗流的影响，建立了高含硫气藏元素硫沉积模型。该模型在完善硫沉积伤害模型基础上，对比了考虑非达西和达西流动下，不同径向距离，不同时间元素硫饱和度随时间的变化关系，同时还分析了产能对硫颗粒沉积堵塞的影响。

2005年，曾平^［47］就高含硫气藏渗流规律进行了研究，得到孔隙度，渗透率随时间的变化关系，进一步完善了考虑非达西影响的元素硫沉积伤害模型。

2006年，杨学锋^［48］在Faruk Civan建立的模型基础上，考虑元素硫沉积的延迟效应，完善了元素硫动态沉积预测模型。

2006年，H.Mei^［49］等人在Roberts建立伤害模型基础上，根据实际井参数，建立了无阻流量与渗透率和储层厚度之间的关系。

2006年，Du Zhi-Ming^［50］等人建立了裂缝性气藏气液固三相耦合数学模型，并利用Roberts实例井数据进行计算，同时进行了结果对比分析。

2006年，Guo Xiao^［51］等人将气液固三相耦合模型与硫沉积实验相对比，分析了流速，初始硫浓度和岩心渗透率对元素硫沉积的影响。

2007年，Guo Xiao^［52］等人基于组分模型和相平衡原理建立了气液固三相数学模型，该模型可用于预测元素硫沉积，并提出需要进行储层解堵时间。

1.2.4 考虑元素硫沉积的产能方程及物质平衡方程研究现状

由于压力降最快的地方在近井地带，导致元素硫析出最快的地方聚集在近井地带，从而使得常规的产能方程需要进一步考虑元素硫沉积的影响。含硫气藏开发过程中元素硫沉积而导致试井曲线发生变化，对此学者们也进行了相应的研究。

2005年，李成勇^［53］等人进行了高含硫气藏解释方法研究，建立了高含硫气藏两区复合试井模型，并用Stehfest反演算法对井底压力响应典型曲线进行了计算，分析了污染半径和流度比对井底压力动态的影响。

2007年，段永刚^［54］等人建立了基于含硫气藏与井筒耦合的非稳态产能预测新方法，该方法为没有试采资料的气井合理配产提供了一种方法。

2008年，张烈辉^［55］等人基于渗流力学相关理论，对高含硫气藏的渗流模式进行了分析，建立了考虑附加表皮的复合渗流模型与产能试井解释数学模型。

2009年，晏中平^［56］等人在现代试井解释方法和油气渗流理论基础上，建立了考虑含硫气井硫污染区和未污染区两区双孔介质复合试井解释数学模型，并利用Stehfest反演算法对井底压力响应典型曲线进行了计算，同时完成了多参数对井底压力的敏感性分析。

随着高含硫气藏的开发，储层压力会不断降低，析出的元素硫将会占据储层部分孔隙空间，使得在建立含硫气藏物质平衡方程的时候，体积平衡方程发生了变化。

1936年，R.J.Schilthuis^［57］根据物质平衡原理首先建立了油藏的物质平衡方程式，因为该方法需要的相关地质及流体生产数据较少，同时计算方法相对简单，故一直在油藏工程中得到广泛使用。

国内的陈元千^{［58～60］}等人在物质平衡原理的基础上建立了气藏的物质平衡方程，并完善了不同类型的气藏物质平衡方程式。

在凝析气藏物质平衡方程式的问题上，国内的马永详^{［61～62］}利用摩尔平衡原理对凝析气藏物质平衡方程进行了研讨。

2006年，张勇^［63］等人给出了高含硫气藏物质平衡方程的推导，该模型考虑了元素硫沉积的影响，但仅仅是基于体积平衡原理，没有考虑元素硫的析出会导致混合天然气密度发生变化。

2008年，卞小强^［64］考虑了元素硫析出后，会使得天然气密度发生变化，必须使用质量平衡原理来建立含硫气藏物质平衡方程，故其利用摩尔平衡原理建立了气藏物质平衡方程，并进行了实例计算，但在建立物质平衡时，由于对元素硫产生的机理认识不足，使得摩尔平衡原理建立的方程求解具有一定难度。