天然气模型_天然气动态仿真真实数据分析

1.现货天然气行情怎么分析?

2.天然气水合物热动力学模型

3.石油天然气关键参数研究与获取

4.开展石油天然气动态评价

（1）、重点技术线分析

天然气投资与投机只有一字之差但是含义却大不相同，投机是短线的操作，讲究的是对风险的控制与人品，而投资则是指中长期的行为。

作为投资者来说一定要学习并掌握一些重点技术线的分析，不能盲目做单，必须要根据现货天然气的行情与重点技术线的分析来做单。如果，有一个专业的老师或团队指导你长期做单，那也许会更容易上手。

（2）、对基本面的关注

天然气这种投资产品跟其他一般的投资产品还不一样，它的特殊性在于，它不仅仅是一种商品也能说是一种政治品种。

影响天然气的原因有很多，天然气的产量，天然气主要产于中东地区，导致中东地区战争局势也会影响天然气的价格变动，还有可以代替天然气的新能源的发展情况，以及每周三公布的EIA数据，美国的非农数据等，都可以影响天然气价格。所以做天然气投资关注时事也是非常必要的。

现货天然气行情怎么分析?

（一）无机与有机天然气类型划分

天然气成因类型的判识主要依赖于天然气的组分和碳、氢同位素组成，并以天然气伴生的轻质油、凝析油、原油的轻烃地球化学特征以及稀有气体同位素组成为辅。腰英台地区的甲烷碳同位素明显偏重，其δ13C1＞30‰。据戴金星（1992），除高成熟和过成熟的煤型气外，δ13C1＞-30%。的均为无机成因的甲烷，因此利用CH4（%）与δ13C1（‰）图可知（图3-33），腰英台构造带主要分布煤型气区内，ChaS1井与YS1井（3466m）登娄库组可能为无机成因甲烷气或者少量的无机气混入的有机气，另外ChaSl井区的个别样品介于无机气与有机气之间，从而表明此研究区有深部的无机气混入，达尔罕构造带以及双坨子地区主要分布有机成因煤成气，煤型气与油型气需要进一步的判识（张枝焕、童亨茂等，2008）。

图3-33 无机与有机天然气类型划分

1—YS1（K1d）；2—YS1（K1yc）；3—YP1（K1yc）；4—YP7（K1yc）；5—YS2（K1yc）；6—DB11（K1yc）；7—D2（K1yc）；8—DB33井区；9—ChaS1井区；10—双—坨子地区

（二）有机烷烃气体进一步鉴别

在有机成因的烷烃气中，生物气和裂解气均具有高甲烷含量、低重烃含量的特点，它们的区别之一是生物气甲烷碳同位素较低，而裂解气的甲烷碳同位素值偏重，根据生物气的一个良好鉴别标志δ13C1＜-55%来看，长岭断陷天然气均属于裂解气。从δ13C1—1gC1/C2+3关系图来看（图3-34），腰英台构造带与ChaS1井区的天然气均属于煤型气，ChaS1井个别样品明显有无机气的混入，为煤成气与无机气的混合气。双坨子地区与腰英台地区的天然气组成特征明显存在差别，主要为原油伴生气以及凝析油与原油伴生气的混合气，由此表明两研究区的天然气的气源是不一致的，腰英台与达尔罕构造带的天然气主要为腐殖型干酪根裂解气，而非原油裂解气（张枝焕、童亨茂等，2008）。

苏联学者Гуцадо（1981）从CH4与CO2共生体系碳同位素热平衡原理出发，以世界上已有CH4与CO2共生体系中测得的δ13C.和δ13Cco2为依据，将自然界不同成因类型的CH4与CO2共生体系划分为三个区，即Ⅰ区为无机成因区，Ⅱ区为生物化学气区，Ⅲ区为有机质热裂解气区。根据图3-35不难看出，研究区腰英台构造带主要分布有机质热裂解气，YS1井与YS2井营城组天然气个别样品分布在无机气的成因区域，大部分样品介于有机质热裂解气区与无机成因气区，达尔罕构造带的天然气主要为有机质裂解气，因此腰英台构造区块的天然气极有可能存在混源特征，可能有无机气的混入，其混源单元还需要进一步的鉴别。

图3-34 天然气δ13C1—lg（C2+（C3）关系图

1—ChaS1井区；2—双坨子地区；3—YS1（K1yc）；4—YS1（K1d）；5—YP1（K1yc）；6—YP1（K1yc）；7—YS2（K1yc）；8—DB11（K1yc）；9—DB33井区

图3-35 CH4与CO2共生体系碳同位素分布图

1—YS1（K1d）；2—YS1（K1yc）；3—YP1（K1yc）；4—YP7（K1yc）；5—YS2（K1yc）；6—DB11（K1yc）；7—D2（K1yc）；8—DB33井区

（三）无机成因甲烷气及识别标志

自然界烃类的大规模形成是有机－无机物质相互作用的结果，而现今油气勘探都是在有机烃源发育的盆地中进行，有机和无机烷烃气混合成藏使无机烷烃气不如非烃气易于识别。尽管如此，目前在许多裂谷盆地中发现了一系列可能的无机成因天然气的聚集，说明无机成因油气仍有一定的发展前景。

到目前为止，对无机成因烃类气体的判断主要依据有烃类气体的组分、碳同位素、烷烃碳同位素系列、与烃类气体伴生的非烃气体、稀有气体的含量及同位素以及地质背景综合分析等方法。松辽盆地有无机成因CH4的一些重要判别依据：

1.该区与无机CO2气藏等伴生的CH4气藏，有特高甲烷碳同位素及负碳同位素系列

在松辽盆地送的与无机CO2气藏等伴生的甲烷碳同位素分析样品，碳同位素值出现了大量的δ13C1值大于－30‰，其中还有大量大于－20‰的样品，并出现了大量负碳同位素系列样品，且上述两种特征还同时出现在同一气田（藏），显示了无机成因烃气的存在。

碳同位素是判识无机成因天然气最直接的证据。我国许多地区如云南腾冲县澡塘河、四川甘孜县拖坝、吉林长白山天池、内蒙古克什克腾旗热水镇以及国外许多地区如新西兰地热区、东太平洋热液喷出口、俄罗斯希比尼地块岩浆岩、美国黄石公园等都发现了无机CH4。这些地区的甲烷碳同位素虽然变化较大，但一般都大于30‰。

许多学者亦提出了鉴定无机成因CH4的下限值，有的为大于－20‰，有的为－30‰。但必须指出的是不论哪一个值都不是划分无机甲烷的绝对值，因为某些高（过）成熟的煤型CH4也有显示重碳同位素特征的特点，因此在确定其成因时还需综合考虑其他资料，如烷烃气碳同位素系列、地质构造背景等。其中碳同位素系列是识别有机、无机烷烃气最有效的手段之一。

有机成因的天然气主要源于沉积物中分散有机质的分解。在生烃母质干酪根热降解生成烷烃气的过程中，由于12C—12C键的键能低于12C—13C键，因此生物成因天然气中CH4及其同系物的碳同位素组成具有随碳数的增大而变重的分布特征，即δ13C1＜δ13C2＜δ13C3＜δ13C4正碳同位素系列。这种分布特征几乎存在于所有有机成因的天然气藏，并被有机质热解成烃的模拟实验和理论推导所证实。而对于无机成因的烷烃气来说，重烃气含量很少，而且主要是由甲烷通过放电作用聚合形成的。在由CH4聚合形成高分子烃类或CO加氢合成烃类的过程中，由于12C—13C键的键能低于12C—12C键，使12C随分子量的增加而逐渐富集，从而形成甲烷同系物的碳同位素组成与有机成因的同位素系列正好相反，即形成δ13C1＞δ13C2＞δ13C3负碳同位素系列。如前面提到的俄罗斯希比尼地块与岩浆岩有关的天然气中δ13C1为3.2‰，δ13C2为9.1‰，δ13C3为16.2‰；美国黄石公园泥火山气的δ13C1为21.5‰，δ13C2为26.5‰。

徐家围子断陷在昌德、汪家屯、肇州以及朝阳沟等地区及腰英台气田均发现了甲烷碳同位素异常和负碳同位素系列，表明该区有无机烃类气体存在。汪家屯地区W a903井甲烷碳同位素最重达12.22‰，而乙烷的碳同位素为22.99‰；昌德地区表现的最为明显，FaS1、FaS2等井多个气样显示负碳同位素系列，且甲烷碳同位素偏重。从这些气样组分来看，干燥系数 一般都在0.98以上，显得很干，也与无机成因烷烃气的特征相似。

此外，也有学者提出负碳同位素系列并不是判断无机成因烃类气体最可靠的标志，由两种不同成因天然气混合，或由天然气的扩散引起同位素分馏均可造成这种现象的出现。以往的研究认为混合作用形成甲烷至丁烷碳同位素的完全反序排列可能性不大，但最近的同位素数值模拟研究结果表明，两种碳同位素正序排列的天然气，混合后可以得到碳同位素完全反序排列的天然气，但要求混合的两个端元的天然气必须具有不同的成因或来源，或它们是明显不同演化阶段的产物。从徐家围子地区的地质条件和同位素特征来看，很难用两种有机成因的气混合加以解释，因为要得到FaS1、FaS2那样重的甲烷负碳同位素系列，要求具有有机成因天然气甲、乙、丙碳同位素为15‰，－14‰，13‰相当的天然气存在，而这种天然气无法与有机质演化的任一阶段相对应，在徐家围子地区也未发现具这种特征的天然气。因此，混合作用不能合理解释该区存在的负碳同位素系列。

2.在该区火山岩的原生流体包裹体中发现CH4

地球深部流体的性质和成分是当前国内外学术界争论的热点课题。火山喷发物中含有大量的非烃气体、少量烃类气体、稀有气体以及沿一些深大断裂带及地震期前后有烃类气体、CO2和稀有气体释放已是公认的事实。近年来对火山岩及其地幔岩流体包裹体的研究进一步揭示其流体相主要为H2O、CO2、CH4、N2、H2、H2S及一些稀有气体。地幔物质及其所含流体在横向和纵向上分布也是极不均匀的，如河北坪尖晶石二辉橄榄岩幔源岩气体包裹体中还原性气体含量高达68.0%～93.4%，而山东栖霞大方山二辉橄榄岩样品中还原性气体为8.5%～39.3%。有学者研究了我国华北地区地幔岩的分布，认为地球深部由上到下依次为尖晶石二辉橄榄岩、尖晶石－石榴石二辉橄榄岩和石榴石二辉橄榄岩，分别代表岩石圈地幔和软流圈地幔。其中石榴石二辉橄榄岩中的H2和CH4的含量最高，而尖晶石二辉橄榄岩含H2和CH4相对较低，因而认为地球深部不同圈层可能孕育有不同性质和类型的天然气，由浅至深有H2O→CO2→CH4、H2富集的趋势，其中莫霍面附近可能是CO2的聚集带，岩石圈与软流圈界面附近可能是烃气的富集带，而H2可能有更深的来源。

在该区非气层段火山岩中集的火山岩流体包裹体，普遍有较高含量的无机烃气，证实无机成因烃类气体对该区气藏的贡献不容忽视。从徐家围子地区岩浆火山岩流体包裹体气液相成分来看，岩浆成分由基性变为酸性时，CO2有从少变多的趋势，CH4的变化趋势正好相反，因此上述研究成果及推断可能是正确的。在长岭达尔罕及腰南构造，在DB11 井的4017～4120m井段的基性岩中发现大量含CH4的气液相包裹体，其中CH4的最高含量可达到31.9%，该层测试产纯CH4，而在相邻的DS2井3670～3780m的酸性流纹岩中，产出以CO2为主的气藏，在该层中发育大量含CO2的气液相包裹本。

3.在该区发现大量示指深部低氧逸度环境的伴生气体

在松辽盆地，已发现部分高含H2及CO、H2S气的气藏，反映该区地壳深部存在低氧逸度环境，有利于甲烷的生成。无机成因气中低氧逸度组分往往构成共生组合，如DB11井营城组玄武岩段，H2含量达6%,H2S含量达（30～50）×10-6，与CH4共生。其各项同位素指标均反映这些组分源自无机成因，证实深部存在低氧逸度的大地构造环境。

4.从地质背景综合分析方法证实应当存在无机成因甲烷

一般认为，某些高（过）成熟的煤型甲烷也有显示重碳同位素特征的特点，并经不同成因天然气混合，或由天然气的扩散引起同位素分馏可造成负碳同位素系列。因此，在一些不含煤系的地区，如部分烃类气藏的δ13C1出现明显偏重，且出现负碳同位素系列，但周缘未发现明显的煤系烃源岩，可以确定存在较大规模的无机甲烷供给。

无机CO2与甲烷的共生，在各类有机烃类成藏条件差别不大的情况下，在局部地区出现特高、特大的气藏，或在有机烃类气体供给很少的区带，在圈闭中发现大量甲烷，揭示存在无机成因甲烷的供给。

以腰英台—达尔罕断凸带为例，该带已钻达基岩顶面的D2、DBIl井揭示，经二维、三维地震资料标定，该区周邻不存在煤系源岩，其它方向有机烃源的运移供给路线也很长。但在腰英台深层气田，发现富含CO2（含量15%～24%），以CH4为主（76%～85%）的气藏，也存在甲烷重碳同位素和碳同位素反向序列。在YS1、YS101、YS102、ChaS1、ChaS1-1、ChaS1-2、ChaS1-3井揭示大型腰英台气田，探明天然气地质储量达（600～700）×108m3的情况下，周围的ChaS2、D2、YN1井却仅发现了CO2气，未发现烃类聚集。这些表明腰英台深层气田有天然成因甲烷的混人。

由于岩石圈地幔及地壳深处广泛存在C、H、O、N等元素，无机成因天然气的主要组成是CO2，其次是CH4及N2等，无机成因气藏也是以CO2为主，含部分CH4、H2、N2、CO2等组分。在无机成因的甲烷气苗中，甲烷含量一般在5%～30%，但即使是这种较低含量，无机成因甲烷供给量也远大于有机成因甲烷供给量。19年Welham等指出，东太平洋北纬21°处中脊喷出的热液（400℃）中，含氢气、甲烷的氦，δ13C1值为17.6‰～－15‰，R/Ra约为8，说明这些气体是幔源的。该处喷出的H2的体积浓度为10%，每年喷出H2和CH4分别为12×108m3和1.6×108m3，如果以此喷出速度，即使仅按照与火山热的地质历史100万年来计算，该处喷出的H2和CH4即可达到1200×1012m3、160×1012m3，也远远大于有机物的生烃量。由此也可见，CO2的供给量是何等惊人。

同时在沉积盖层的深埋压实条件下，CO2易于与地壳中碳酸盐岩、碱性岩类发生反应，并大量溶解于水中，而产生大量的损耗。而在地壳沉积盖层的温度、压力条件下，CH4则有相对的化学稳定性，在CO2逃逸和散失量很大的条件下，无机成因CH4常可以形成相对富集，甚至形成无机成因甲烷为主的天然气藏。

（四）煤型气与油型气的鉴别

确认天然气属于煤型气还是属于油型气，对于追溯、对比烃源岩起着重要作用，目前最为常用的参数是乙烷或丙烷碳同位素。YS1井登娄库组天然气δ13C2为－24.7‰，为典型的煤型气，YS1井营城组天然气δ13C2为－26.4‰～－26.5‰，DBIl-1井与DBl1-2井营城组天然气δ13C2为－26.1‰～－28.7‰，均为煤型气和油型气混合气区，DB33-9-3井天然气的δ13C2为－29.3‰，也接近煤型气和油型气混合气区，按照δ13C2值－29%。为界限，长岭断陷天然气为高成熟的煤型气。

1.“V”型鉴别图（δ13C1-δ13C2-δ13C3）

考虑到甲烷、乙烷与丙烷三者碳同位素的综合信息，在δ13C1—δ13C2δ13C3相关图上（图3-36），利用烷烃成因天然气碳同位素系列数据，能够鉴别不同成因的有机天然气。其中Ⅰ区为煤型气，Ⅱ区为油型气，Ⅲ区为混合型气，Ⅳ区为深层混合气（戴金星，1992；顾忆等，1998）。从图3-36可以看出，腰英台构造带与达尔罕构造带的天然气主要分布在碳同位素倒转区以及煤型气和油型气或者深层气的混合气区，而且天然气的成熟度明显偏高，DBll井的天然气可能有少量的油型气混入，双坨子地区的天然气主要为煤型气与油型气，由此表明，双坨子构造带的天然气的特征明显不同于上述两个构造带，腰英台与达尔罕构造带的天然气明显具有多源的性质，而且可能混有深部的无机气，造成其甲烷的同位素明显偏重，导致其烃类组分的同位素发生倒转。

2.δ13C2-δ13C1图

通过利用δ13C2值的大小将天然气划分为煤型气、油型气以及煤型气与油型气的混合气区，再通过δ13C1受热演化程度的差异将天然气划分为未熟、低熟，成熟、高熟以及过成熟五个阶段，可以很好地将天然气中煤型气与油型气类型分开，从图3-37可以看出，腰英台与达尔罕构造带的DB33-9-3、DB33-5-5、DB11井以及ChaS1井的个别样品可能为高过成熟的煤型气与油型气混合气，而其余样品天然气均为高过成熟的煤型气，双坨子地区的天然气成熟度略低，分布油型气或煤型气，不同于腰英台与达尔罕构造带的天然气的特征。

图3-36 天然气δ13C2-δ13C1不同成因类型图

1—ChaS1井区；2—双坨子地区；3—YS1（K1d）;4—YS1（K1yc）；5—YP1（K1yc）；6—YS7（K1yc）;7—YS2（K1yc）;8—D2（K1yc）;9—DB11（K1yc）；10—DB33井区

图3-37 天然气δ13C2—δ13C1不同成因类型图

1—ChaS1井区；2—双坨子地区；3—YS1（K1d）；4—YS1（K1yc）；5—YP1（K1yc）;6—YP7（K1yc）,7—YS2（K1yc）;8—D2（K1yc）;9—DB11（K1yc）;10—DB33井区

3.C1/C1-5与δ13C1图

利用干燥系数（C1/C1-5）与δ13C1同样也可以判识天然气类型.对于煤型气与油型气在不同的演化阶段过程中，其干燥系数与δ13C1存在一定的对应关系，对于成熟度高的油型气与煤型气，其干燥系数与δ13C1必然很高，图3-38中A1、B1、C1、D1、E1为煤型气演化阶段，界限由虚线表示，A2、B2、C2、D2、E2为油型气演化阶段，界限为由实线表示。通过图3-38可以看出，腰英台构造带与达尔罕构造带的营城组与登娄库组的天然气主要分布在高成熟的煤型气与油型气区，双坨子地区天然气具有煤型气与油型气的混合特征，明显不同于两构造带的天然气特征。

图3-38 利用C1/C1-5与δ13C1图判别不同类型烷烃气体

1—ChaS1井区；2—双坨子地区；3—YS1（K1d）;4—Ys1（K1cy）；5—D2（K1cy）；6—YP1（K1yc）；7—YP7（K1yc）；8—YS2（K1yc）;9—DB11（K1yc）;10—DB3井区

（五）天然气同位素倒转现象分析

长岭断陷腰英台与达尔罕构造带天然气碳同位素系列数据分析表明，碳同位素倒转系列和负碳同位素系列是其主体，并且碳同位素明显偏重。导致碳同位素异常的原因有很多，研究天然气碳同位素倒转的原因，对天然气的成因或其经受的次生变化作出判断，可以作为天然气运移途径和气源对比的一种间接方法。戴金星（1993）曾对烷烃气碳同位素系列倒转问题作过详细研究，认为引起碳同位素系列倒转的主要原因有：1）有机气与无机气的混合，二者分别属于正碳同位素系列与负碳同位素系列的典型，当二者混合时，很容易发生同位素分布的倒转现象；2）煤型气与油型气的混合，这是造成碳同位素系列倒转的主要原因；3）同型不同源或同源不同期天然气的混合，同源的早期形成的低成熟度的天然气散失一部分后的剩余气，与晚期较高成熟度形成的天然气形成混合天然气，可导致烷烃气同位素倒转；4）生物降解作用，细菌选择降解某些组分致使剩余组分变重；5）地温增高也可使碳同位素倒转，在碳同位素交换平衡下，若地温高于100℃，则出现正碳同位素系列；当温度高于200℃时，则正碳同位素系列改变成为负碳同位素系列（戴金星，1990）；6）源岩性质控制，在中国陆相河湖交替发育的含油气盆地，烃源岩有机质的分布是不均一的，同一套烃源岩中I型和Ⅲ型有机质可能同时存在，因此其产生的烃类烷烃气可能发生倒转，松辽盆地北部深层烃源岩就有混源的特点。

此外，盖层微渗漏造成的蒸发分馏作用也是许多天然气藏同位素出现倒转的重要原因，Prinzhofer等（1995）在对Jenden的资料进行重新解释时，认为微渗漏作用更能合理地解释Appalachian盆地天然气同位素的倒转现象，他们按Jenden等提出的混合模式计算后发现有些样品点并不符合混合模式，提出了一种新的微渗漏模式。黄海平（2000）利用微渗漏模式较好地解释了徐家围子断陷深层天然气同位素倒转的现象。从图3-39看出，腰英台构造带的ChaS1井区、达尔罕构造带的DB11-1、DB11-2、DB33-9-3、DB33-5-5等井天然气样品同位素发生倒转，是受到盖层微渗漏作用的影响。

导致天然气碳同位素倒转可能是上述因素之一，也可能是两种或两种以上的因素引起的。长岭断陷深层天然气普遍被认为主要来源于沙河子组和营城组，经历了较复杂的构造变形和较高的成熟演化阶段，可能存在多源气的混合，主力烃源岩发育于盆地断陷晚期和坳陷早期，火山活动频繁，烃源岩除正常的热演化外，还受到因火山活动引起的异常热，主力烃源岩沙河子组和火石岭组在盆地分布不均一，有机质具有非均质性，因生气层上下部位和层内成熟度及有机质性质不一样，也会使同层同时生成的天然气同位素发生混合而倒转。盆地基底发育深大断裂，无机成因的CO2、N2普遍存在，并且丰度较高，在腰英台地区CO2含量平均值为20%以上，因此天然气中可能有无机成因烷烃气加入，天然气藏产层主要在登娄库组与营城组，成藏模式比较复杂，天然气可能以垂直运移为主，运移路径较长，因而可以引起多期次的天然气碳同位素动力分馏效应。

图3-39 天然气同位素反转解释模式

1—ChaS1井区；2－双坨子地区；3—YS1（K1d）；4—YS1（K1yc）;5—D2（K1yc）;6—YP1（K1yc）;7—YP7（K1yc）;8—YS2（K1yc）;9—DB11（K1yc）；10—DB33井区

据此按照通常的天然气同位素的划分，结合长岭断陷腰英台地区天然气各种分析数据可知，YS1井登娄库组以及ChaS1井个别样品表现出无机成因气的特点，而腰英台构造带大部分井区的样品，如YS1、YS2、YP7井以及达尔罕构造带的DB33井区、DB1I井主要分布有机成因的烷烃气（张枝焕、童亨茂等，2008）。

天然气水合物热动力学模型

现货天然K线图包含了最高价、收盘价、开盘价、最低价这四个最基本的数据，所以就可以根据K线的形态判断出交易时间内的做多、做空的情况。

当开盘价和最低价相等，收盘价和最高价相等的时候，K线被称为大阳线，表示涨势强烈；当开盘价和最高价相等，收盘价和最低价相等的时候，K线被称为大阴线，表示大跌。

如果利用裸K判断走势，那么则需要划趋势线，如果利用技术指标MA均线以及boll看盘，那么则寻找支撑点位，针对MA均线更适合投资者学习

石油天然气关键参数研究与获取

一、天然气水合物热力学模型

1.理论基础

随着各种热力学研究的开展，现已有大量有关天然气水合物相平衡的数据和方法，可用来预测天然气水合物的形成。这些研究结果也有利于开发抑制天然气水合物形成的化学添加剂。一般说来，能影响溶液黏度性质的物质通常能抑制天然气水合物的形成。在工业应用上，甲醇是一种常用的阻凝剂。

Van der Waals和Platteeuw(1959)提出的热力学理论，一直是预测天然气水合物平衡模型的理论基础。Sloan(1990)指出，利用这些模型对Lw-H-V系统平衡压力的预测，误差应该不超过10%，而对温度的预测误差在2K左右。多年来，各国学者在Van der Waals和Platteeuw理论的基础上，提出了一些新的观点和天然气水合物相平衡预测的计算方法，对天然气水合物热力学模型的发展作出了贡献。

2.热力学模型

要描述各种天然气水合物相及其可能的多种共存流体相，需要使用一种以上的模型。状态方程是描述天然气水合物平衡的一种方法。为描述富水的流体相，Saito(1964)等使用了理想溶液方法(Raoult定律)，其基本前提是，设水中储存气体的溶解度在常规条件下可以忽略不计，尽管有事实证明这种设的有效性令人质疑，但这种方法在过去一段时间内仍为大多数学者所偏爱。当需要进行天然气水合物抑制计算时，要根据Menten(1981)提出的计算方法，用活度系数对Rao-ult定律进行校正。虽然该方法的可靠性要优于Hammerschmidt方程(1939)，但它不能用于评估阻凝剂(如甲醇)在共存相中的分布。为校正这个问题，Anderson等(1986)结合使用了Uniquac方程和用于超临界组分计算的亨利定律，计算液相中所有可凝聚组分的逸度。因此，要进行简单的天然气水合物抑制计算，有必要使用上述4个模型。由于这种内在的复杂性，对于现实中更复杂的系统，上述这些方法用处并不大。同时，这些方法都存在着收敛困难，不能作为进一步精确计算(如复合系统的稳定性分析)的基础。

Englezos等(1991)和Avlonitis等(1991)根据一个单一的状态方程，建立了全部流体相的模型。前者使用了有4个参数的立方状态方程，后者使用了有3个参数的立方状态方程，并开发了针对非对称相互作用的专用混合规则。根据目前的研究趋势看，对全部流体相使用单一的状态方程是最有发展潜力的方法。

3.模型的完善和发展

对天然气水合物相的理想固溶体，在设被圈闭的分子之间没有相互作用的前提下，Vander Waals等(1959)认为能够用一种Langmuir型吸附等温线描述固体天然气水合物相。他们利用这个设，证明天然气水合物相中水的化学势能与形成天然气水合物的气体性质无关，仅取决于天然气水合物相中两种不同类型空穴中气体的总浓度，天然气水合物与理想稀溶液具有相同的行为。在这个理论基础上，Parrish等(12)将用于计算分解压的天然气水合物模型延伸到多组分系统中。理想固溶体理论忽视了非理想状态所带来的影响，如“主”分子晶格的伸展或变形、被圈闭气体分子运动所受的限制，都有可能增加“主”分子和“客”分子的化学势。Hwang(1993)与他的同事们在分子动力学模拟的基础上，研究了“客”气体分子的大小对天然气水合物稳定性的影响。Avlanitis(1994)指出:这种方法的主要缺陷在于选取了不正确的势能参数，特别是乙烷的势能参数。为弥补这个缺陷，Avlonitis用一种折中方法优化了理想固溶体模型参数，在含甲醇或不含甲醇情况下，在Lw-H-V框架内，对天然的或合成的气体混合物都获得了令人满意的预测结果。

二、天然气水合物动力学模型

天然气水合物动力学是水合物领域的研究重点。通常以方程M+nH2O＜=＞［M·nH2O］表示水合物生成，这是一个气-液-固三相或气-固两相的多相反应过程，同时也是一个包含传热、传质和生成水合物反应机理的复杂反应方程，影响反应的条件很多，也很复杂。相对于天然气水合物热力学而言，对天然气水合物动力学的研究较少。天然气水合物动力学可以大略分为天然气水合物生成动力学和天然气水合物分解动力学两类。

1.天然气水合物生成动力学

针对目前研究亟待解决的水合物生成速率和效率问题，主要有以下两种研究方法(赵义等，2004):①热力学方法，即向反应体系中加入其他气体添加剂，让气体添加剂占据水合物结构中没有被占据的空腔，来降低水合物簇之间的转换活化能，提高水合物的晶体空腔填充率，从而达到促进水合物生成和提高水合物稳定性的目的，如向甲烷的水合物生成体系中加入少量的丙烷，就可以大大降低甲烷水合物的生成条件，并且生成的结构更稳定;②动力学方法，仅限于表面活性剂及助溶剂(hy-drotropes)的研究。对此有两种说:一是Sloan的观点，认为表面活性剂之所以促进水合物的生成，主要是因为它降低了气-液界面张力，增大了扩散传质速率，使气体更容易进入液相;二是Zhong等(2000)的观点，提出了一个4步骤的反应历程来解释观察到的现象，尚未得到充分的验证。以下对这4个步骤进行说明:

(1)气-水簇的形成

天然气水合物的成核过程是天然气水合物核向临界大小的靠近且生长的过程。气-水生长簇是天然气水合物形成的先兆。如果生长的核小于临界大小，核是不稳定的，可能在水溶液中生长或破裂。一个生长着的天然气水合物核，如已具有临界大小，就是稳定的，可以立即导致结晶天然气水合物的形成。

认识影响气-水簇形成的因素，有助于理解天然气水合物的成核过程。特别是水分子结构，它是指通过氢键相互联结水分子所形成的结构，在成核过程中起重要作用。冰是一种高度结构化的水，其水分子固定在一个呈四边形氢键结构的位置上。当温度升高到零点以上时，结构开始变得更加松散，与高度有序的冰的结构相比更加无序。

Sloan等提出了一种天然气水合物成核过程的分子机制，设想气-水簇开始形成临时结构，随后这些结构再生长成稳定的天然气水合物晶核。他们通过使用化学动力学方程，针对机制中设的每一种情况对成核过程进行了模拟。Lekvam和Ruoff也提出了反应作用的动力学机制。这种方法使用一种动力学速率模型，研究成核和生长过程，但他们的这种模型并没有强调天然气水合物核的稳定性。

Vysniauskas和Bishnoi在实验中使用不同来源的水进行了实验。结果显示，随着水的来源不同，平均成核开始期也不同。在实验中，来自于融化的冰水与实验中使用热自来水相比，前者的平均开始期较低;同样，使用来自于天然气水合物分解的水与使用热自来水相比，前者的平均开始期也较低，这就是所谓的“记忆效应”。这种现象在其他学者的研究中也出现过。研究发现，在已溶解的气体分子周围，水结构被强化了;这种作用于溶解气体分子周围的水结构强化现象，被认为是“疏水水合作用”现象。Frank等也提出了同样的观点。Glew在对甲烷天然气水合物和甲烷水溶液的热动力学性质进行研究时，发现了类似现象。Glew对甲烷-水系统分子模拟的研究显示，甲烷分子周围的水分子平均配位数对于Ⅰ型结构的小空穴来说，接近于21。Rahman和Stillinger认为，在溶解的溶质分子周围，水的网架与天然气水合物类型的孔型相似。另外，热力学分析显示，溶液具有很大的负熵，这正是水体内一种结构形成的标志。

气-水簇在天然气水合物成核过程中起了很重要作用。当溶液在过冷或过饱和状态下时，成核过程就发生了，学者们通常使用过冷或过饱和方法来研究成核作用。Bishoni等在研究时就用了过饱和方法，Kobayashi、Sloan等则用了过冷方法。

图10-5 典型的气体消耗简图

Bishnoi等在恒定温压下进行了天然气水合物形成实验。在实验温度下，实验压力比三相天然气水合物平衡压力要高，图10-5是实验过程中因气体溶解和天然气水合物形成，而导致的气体消耗的累积摩尔量随时间变化的曲线。

图10-5中A点的气体消耗摩尔量表示已溶解气体量，与三相天然气水合物平衡压力对应。A点与B点之间的准稳区域，代表着天然气水合物的成核过程。B点表示以突变方式出现的稳定临界大小天然气水合物核的出现点。Englezos和Bishnoi发现，在成核点B之前的溶解气体摩尔量，实际上要高于估算的二相(气-液)准稳定平衡状态下的摩尔量，估算来源于稳定区域的外推;气-水簇的形成能够耗尽在团块流体相中的天然气水合物形成的气体，从而导致超过两相值的气体溶解。Englezos等提出了计算天然气水合物核临界大小的方法，天然气水合物生长过程开始于图10-5中B点，并沿着线BC进行。根据Kobayashi和Sloan的实验结果，在容积不变的情况下，天然气水合物形成过程中的压力和温度轨迹如图10-6所示。图10-5中点A等同于图10-6中的点A。图10-6中点B也等同于图10-5中的点B，在B点，以突变方式形成的稳定天然气水合物颗粒的出现，导致了压力的突然下降。在图10-6中，点A与点B之间区域表示成核过程中的准稳定状态。过冷却方法和过饱和方法的相似性在图10-5与图10-6之间体现得相当明显。在图10-5中，与三相天然气水合物平衡相应，点A与点B之间，是处于准稳定状态的天然气水合物成核区域中气-水簇的生长区域。天然气水合物在点B的出现是突然的，Kobayashi描述它为“灾变性的”。尽管天然气水合物颗粒很小，但它们的数量很多，足以使溶液变得混浊。Kobayashi和Sloan认为，天然气水合物的突然出现使溶液不再处于过饱和状态，这样便导致了压力的剧降。

(2)天然气水合物的成核和生长过程

图10-6 天然气水合物形成实验温度-压力轨迹简图

从上面讨论可以看出，过冷方法与过饱和方法是等价的，对于天然气水合物成核过程来说都很重要。许多研究者建立了开始期和过冷之间的函数关系，过饱和同样也可以根据过冷却度进行转换。溶解中任何点的过饱和，都是在这点超过饱和浓度值的过量溶解气体浓度，可以用溶液中某一点的过饱和来判断稳定天然气水合物核最先出现在哪个地方。对于不流动系统，溶解气体浓度在分界面附近可能最高，天然气水合物的形成可能最先发生在气-液分界面上。对于搅拌系统来说，在溶液中最先形成天然气水合物的地方，取决于这点上溶解气体的浓度。溶液的水动力条件和气体溶解率可以影响天然气水合物成核的开始期。

Bishnoi等认为，天然气水合物成核作用开始期与过饱和作用相联系，根据对甲烷、乙烷以及二氧化碳天然气水合物的实验数据分析，揭示了成核开始期与过饱和的关系。当过饱和度减小时，成核开始期增大;当过饱和时，开始期增加到一个很大的值;相反，当过饱和度增加时，开始期减少到一个很小的值;当过饱和度很低时，开始期数据的分散程度很高，当过饱和度增加时，开始期数据的分散程度减小。天然气水合物成核过程，本质上是一个内在的随机过程，但高的过饱和度能够掩盖成核现象的随机本质，从而使观察到的开始期看起来像是早已被决定了一样。另外，天然气水合物成核的随机本质，也能够被实验系统中用来进行成核研究的其他因素所掩盖。在天然气水合物成核研究中，Parent和Bishnoi在原始实验状态下又观察到了开始期数据的随机性。

关于天然气水合物成核的研究还处于宏观层次上。对在溶液中的亚临界情况，还知之甚少。在建立基于分子级的模型之前，须通过实验研究揭示天然气水合物的成核机制。

天然气水合物的生长过程，是指作为固态天然气水合物的稳定天然气水合物核的生长，自20世纪60年代以来，许多学者就已对此进行了研究。在研究丙烷天然气水合物生长过程时，Knox认为晶体的近似大小取决于过冷度(指使液体冷到凝固点以下而不凝结)，较高的过冷度主要产生较小的颗粒，并导致明显的晶体生长。Pinder通过研究天然气水合物形成动力学，提出天然气水合物形成的反应速率随渗滤作用而定。Barrer和Esge在研究天然气水合物动力学时发现，对氪形成的天然气水合物来说，其晶体生长有一个明显的开始期。Falabella使用类似于Barrer和Esge的实验装置进行了研究，也得到了相似的结论。Falabella还发现，对于甲烷来说，其天然气水合物生长也有一个开始期，他根据冰的动力学数据，通过进行等温压换算，提出了一个次级动力学模型。Sloan和Fleyfel通过实验，研究了环丙烷天然气水合物的生长动力学。针对在纯水中的各种气体和气体混合物，Bishnoi等一直进行着天然气水合物形成动力学的系统性研究，在实验中使用一个搅拌反应器，其中装有电解质和表面活化剂溶液。他们认为，在稳压条件下，全部气体消耗量是时间的函数。

(3)天然气水合物生长动力学模型

在研究早期，Vysniauskas和Bishnoi提出了一个关于气体消耗速率的半经验模型。后来，Engl-ezos等把只有一个可调节参数的天然气水合物生长动力学模型公式化，这个模型是一个以结晶化和团块传递理论为基础的模型;它设固体天然气水合物颗粒被一个吸附“反应”层所包围，吸附反应层外是一层不流动的液体扩散层，溶解的气体从围绕在不流动液中向天然气水合物颗粒-水分界面扩散;然后，气体分子由于吸附作用而进入结构化的水分子构架并结合在一起。当水分子过量时，分界面被认为是气体最易集中的地方(反应速率用已溶解气体的逸度替代其浓度)。

在三相天然气水合物平衡压力和颗粒表面温度下，在扩散层中，溶解气体逸度值从fb变化到fs;在吸附层中，逸度值直降至feq，围绕颗粒的扩散动力等于fb-fs;但是对于“反应”阶段来说，这个值是fs-feq。在稳定状态下，扩散阶段和“反应”阶段的速率相等，fs能够从单个速率表达式中消去，可得到每一个颗粒的生长速率如下:

非常规油气地质学

式中:R*是扩散和吸附反应过程的组合速率常数;Ap是每个颗粒的表面积。在溶解气体的逸度中，fb-feq值不同于三相平衡逸度中的fb-feq值，它指的是全部动力。当在良好的搅拌系统中时，R*值表示反应的内在速率常数，R*值由甲烷和乙烷天然气水合物形成动力学的实验数据决定。在没有任何附加参数的情况下，这个模型可成功地扩展到甲烷和乙烷混合物的形成动力学;在这个模型中，纯水中甲烷天然气水合物形成时获得的R*值，可以应用到电解质溶液中的天然气水合物形成模型中，两者的R*值是相同的。

在液态二氧化碳和水的分界面上，Shindo等提出了二氧化碳天然气水合物形成模型;他们设天然气水合物主要发生在液态二氧化碳中，而不是在水中。最近，Skovborg和Rasmussen使用实验的气体消耗数据(数据来源于Bishnoi等)，提出了一种天然气水合物生成动力学模型;认为天然气水合物的形成，能够影响液体一侧的气-液团块传递系数。

(4)气-水体系中水合物的生成机理

天然气水合物结构和性质类似于冰(陈孝彦等，2004)，气-水体系中天然气水合物生成时，气体分子首先要溶解到水中，一部分气体分子与水一起形成水合物骨架，类似于冰的碎片(周公度等，1995)，形成了水合物结构中的第一种空穴。这些框架是一种亚稳定结构，相互结合形成更大的框架。在结合过程中，为保持水分子的4个氢键处于饱和状态，不可能做到紧密堆积，缔合过程中必然形成空的包腔，就形成了水合物结构中的另外一种空穴。另一部分溶解的气体分子通过扩散渗入到这些空穴中，并进行有选择的吸附;在吸附过程中满足Langmuir吸附定律，小气体分子进入小空穴，同时也能进入大空穴，大气体分子只能进入大空穴，即并不是每一个空穴都能被气体分子占据，这就解释了水合物平均只有三分之一左右的空穴被客体分子占据的机理。

陈孝彦等(2004)总结提出了气-水体系中水合物的生成机理，分为4步:①气体分子溶解过程，即气体分子溶解到水中;②水合物骨架形成过程，即气体分子的初始成核过程，溶解到水中的气体分子和水，形成一种类似冰碎片的天然气水合物基本骨架(一种空腔)，这种骨架通过结合形成另一种不同大小的空腔;③气体分子扩散过程，即气体分子扩散到水合物基本骨架中;④气体分子被吸附过程，即天然气气体分子在水合物骨架中进行有选择的吸附，从而使水合物晶体增长。

2.天然气水合物分解动力学

(1)理论基础

人们提出了许多基于相平衡的热力学模型来预测一定条件下水合物的生成条件及其抑制途径(赵义等，2004)，如通过改变其生成条件，来达到抑制目的的物理方法，包括干燥脱除法、加热保温法、降压法和加入非水合物形成气体法等，还包括通过加入添加剂的化学方法。

化学抑制法主要有热力学抑制剂和动力学抑制剂两种(赵义等，2004)。前者普遍取在生产设备和运输管线中注入甲醇、乙醇、乙二醇和氯化钠、氯化钙等，改变水合物热力学稳定条件，抑制或避免水合物生成;后者从降低水合物生成速度，以抑制水合物晶粒聚结和堵塞出发，通过加入一定量化学添加剂来改变水合物形成的热力学条件，显著降低水合物成核速率，延缓乃至阻止临界晶核生成，干扰水合物晶体的优先生长方向，影响水合物晶体的定向稳定性，具有用量少、效率高等优点，已成为了研究热点(吴德娟等，2000)。根据分子作用的不同机理，动力学抑制剂分为水合物生长抑制剂、水合物聚集抑制剂和具有双重功能的抑制剂，主要包括酰胺类聚合物、酮类聚合物、亚胺类聚合物、二胺类聚合物、共聚物类等，其中酰胺类聚合物是最主要的一类。

Holder等(1987)研究了在天然气水合物分解过程中的热传递过程，得出与成核沸腾现象相似的结论。Kamath等(1987)根据这种相似性，提出在丙烷分解期间，热传递率是ΔT的幂函数，其中天然气水合物表面的ΔT值与团块流体中的ΔT值是不相同的。后来，Kamath和Holder总结了它们的关系性，并用到甲烷天然气水合物分解过程中。

Selim等(1989)研究了甲烷水合物的热分解，认为水合物的分解是一个动态界面消融问题，并运用一维半无限长平壁的导热规律，建立了甲烷水合物的热分解动力学模型，Kamath等(1987)研究了甲烷和丙烷的热分解问题，认为水合物的分解主要受传热控制，其分解可与液体的泡核沸腾相比拟，而流体主体与水合物表面的温差ΔT是过程的推动力(Kamath et al.，1987)。

(2)实验研究

对天然气水合物分解动力学的基础研究是在带搅拌的大容积反应器中进行的，水合物以固体颗粒状分散于液体中，这用来研究分解本征动力学是可以的(周锡堂等，2006)。但用于研究与天然气生产有关，特别是天然气水合物分解的反应工程动力学，则缺乏实际意义(周锡堂等，2006)。自然赋存的天然气水合物可能是大块状的，更多的存在于多孔介质中。Sloan等报道过砂岩中的甲烷水合物生成和分解的一些实验数据，但没有仔细地研究水在孔隙里的分布情况;Circone等报道过以冰粒形成的水合物在272.5K的分解速率数据(Circone et al.，2000)，但也没有提供相应的动力学方程。存在于冻土带或海底沉积物中的天然气水合物，与人工合成的、仅仅存在于自由水中的水合物颗粒是大不相同的。因此从工程实际来考虑，研究多孔介质中水合物的分解动力学行为更有意义。Yousif等第一次将水合物分解动力学的研究与天然气的生产结合起来(Yousif et al.，1991)，不过其模型在估算水合物面积时却是经验性的。Goel等研究了天然气水合物的分解行为(Goel et al.，2001)，运用发散状扩散方程，分别得出了关于大块状和多孔介质中的天然气水合物的分解动力学解析模型。然而该模型忽略了分解水的流动和分解气出速率的变化，严重影响了其有效性。Hisashi等研究了多孔介质中水合物的形成和降压分解问题(Hisashi et al.，2002)。在其实验中，分别用了多种粒度的玻璃珠和合成陶粒来模拟多孔介质。最终结果表明，不同介质中水合物分解的表观反应速率常数不同，所得回归方程也不一样(周锡堂等，2006)。因此，在确定自然存在天然气水合物的分解速率时，有必要研究当地介质的孔隙性质及其粒度分布。

Bishnoi等开展了对甲烷天然气水合物分解的实验研究，实验是在一个搅拌良好的反应器中进行的;天然气水合物在三相平衡压力以上存在;然后，在保持温度不变的条件下，把压力降低到低于三相平衡压力，这时，天然气水合物分解就开始了;实验在快速搅拌中进行，以保证避免团块传递的影响。他们提出，天然气水合物分解可能分为两个阶段:颗粒表面原结晶“主”格子破坏和随后的“客”分子从表面的解吸过程。Kim等提出了天然气水合物分解原内在动力学模型，他们设天然气水合物的颗粒为球形，并且被云雾状气体所包围，如图10-7所示。在图中，正在分解的颗粒被解吸“反应”层所围绕，再外层是排放出的气体云，天然气水合物颗粒分解速率公式如下:

非常规油气地质学

式中:kd为分解速率常数;Ap为颗粒表面积;feq为气体三相平衡逸度;fvg为气体分解策动力，定义为feq与fvg之差，即feq-fvg。

(3)研究进展和意义

与前文提到的对天然气水合物生长的研究一样，对天然气水合物分解的研究，应该包括对决定分解颗粒大小分布因素的研究。

图10-7 天然气水合物分解图

对天然气水合物分解和形成动力学的研究，给我们提出了大的挑战。天然气水合物形成被认为是一种包括成核过程和生长过程的结晶化过程。成核作用是一种内在的随机过程，它涉及气-水簇向具临界大小的稳定天然气水合物核的形成和生长问题。因较大的成核策动力和多相性的存在，成核作用随机性质不易被察觉。目前，对天然气水合物成核过程仍没有在分子级别上的测试方法。

天然气水合物生长包括作为固态天然气水合物的稳定水合核的生长，正在生长的天然气水合物颗粒表面积，强烈影响着生长速率。天然气水合物分解是一系列晶格的破坏和气体解吸过程，在分解时的热传递率与成核沸腾现象是相似的。应该深入研究天然气水合物颗粒在分解和生长过程中的大小分布，并应用于这些过程的模型化中。

尽管有多个天然气水合物形成模型已经被提出，但天然气水合物形成核的过程并没有完全被揭示。目前，科学家通过研究气体-水的接触面，已取得了一些实验上的进展，但是这些实验都是最近做的，并且至今没有充足的信息来提供一个确切的描述。这些实验通过研究熔点附近的热力学状态范围，来揭示与接近天然气水合物形成条件相联系的界面结构特征。在实验中，科学家把分子动力学模拟，应用到Ⅰ型甲烷天然气水合物和甲烷气体的接触面，发现接触面在270K以下是稳定的，在300K时发生熔解，同时发现了导致接触面稳定的压力条件。在伴随着表面层的无序化过程中，预熔现象是明显的。动力学性质显示了水平面格子振动的各向异性，这被认为是与在Ⅰ型天然气水合物(001)面上存在着晶轴相联系。这个意想不到的结果还有待于进一步研究。

在研究天然气水合物形成模型的同时，由于天然气水合物有时能对高纬度地区石油和天然气的运输造成意想不到的麻烦(如形成管塞)，有的学者(Monte Carlo)也开始了怎样抑制天然气水合物形成的研究。通过实验研究发现，可以使用一种无毒的、能溶解于水的聚合物———科利当(PVP)，来抑制天然气水合物的形成。Monte Carlo通过不同条件下PVP对单体、二聚物、四聚物、八聚物吸附性的研究，发现吸附作用主要在吡硌烷酮氧(pyrrolidone oxygen)和水面之间两个氢键的形成过程中出现。这种研究结果表明，通过在天然气水合物生长点上PVP的吸附，来抑制天然气水合物的形成是可行的，并且影响吸附的主要因素具有内在的统计性。

开展石油天然气动态评价

评价参数直接影响评价方法的有效性，不同类型的参数作用不同。有效烃源岩有机碳下限、产烃率图版、运聚系数是成因法的关键参数；最小油气田规模对统计法计算结果有较大影响；油气丰度是应用类比法的依据，由已知区带的油气丰度评价未知区带的丰度；可系数是将地质量转化成可量的关键参数。

（一）刻度区解剖

1.刻度区的定义

刻度区解剖是本次评价的特色之一，也是油气评价的重要组成部分。刻度区解剖的目的是通过对地质条件和潜力认识较清楚的地区的分析，总结地质条件与潜力的关系，建立两者之间的参数纽带，进而为潜力的类析提供参照依据。

刻度区是为取准评价关键参数，以保证评价的客观性而选择的满足“勘探程度高、探明率高、地质认识程度高”三高要求的三维地质单元。刻度区可以是一个盆地（凹陷）、一个油气运聚单元、一个区带、一个成藏组合、一个层系或一个二级构造带等。为了正确和客观认识地质条件和潜力，刻度区的选取在考虑“三高”条件的基础上，应尽量考虑不同地质类型的综合，这样可以更充分体现油气丰度与地质因素之间的关系。

2.刻度区解剖内容与方法

刻度区解剖主要围绕油气成藏条件、量及参数三个核心展开，剖析三者之间的关联规律和定量关系。

（1）成藏特征和成藏主控因素分析。成藏特征和成藏主控因素分析实质上是对选择的刻度区进行成藏特征总结，精细刻画出成藏的定性、定量的主控因素与参数，便于评价区确定类比对象。在一个含油气盆地、含油气系统、坳陷、凹陷的成藏规律刻画中，其成藏特征差异大，故一般最好选择以含油气系统（或坳陷）及其间的运聚单元作为对象，更便于有效的类比应用。油气运聚单元是盆地（凹陷）中具有相似油气聚集特征的独立的和完整的石油地质系统，是以盆地（凹陷）的油气聚集带为核心，并包含为该油气聚集带提供油气源的有效烃源岩。油气运聚单元是有效烃源岩、油气运移通道、有效储集层、有效盖层、有效的圈闭等要素在时间和空间上的有机组合。一个油气运聚单元可以有多个有效烃源岩体和烃源岩区为其供烃，但同一个油气运聚单元的油气聚集特征是相似的。一个油气运聚单元可以只包含一个油气成藏组合，也可以包含在纵向上叠置的多个油气成藏组合。因此刻度区地质条件的评价与定量刻画就是按照运聚单元→成藏组合→油气藏的层次路线综合分析烃源条件、储层条件、圈闭条件、保存条件以及配套条件等油气成藏条件。盆地模拟是地质评价流程中的一个重要组成部分，其作用主要体现在三个方面：其一是通过盆地模拟反映流体势特征，进而确定油气运聚单元的边界；其二是提供烃源参数，如生烃强度、生烃量、有效烃源岩面积等；其三是通过关键时刻的获取来反映油气成藏的动态作用过程。

（2）油气量确定。刻度区量计算与一般意义上的量计算稍有不同，正是由于刻度区的“三高”背景，特别是选定的刻度区探明程度越高越好，计算出的量更准确有利于求准各类评价参数。在本次刻度区解剖研究中，主要用了统计法来计算刻度区的量，统计法中包括油藏规模序列法、油藏发现序列法、年发现率法、探井发现率法、进尺发现率法以及老油田储量增长法，不同方法估算出的量用特尔菲加权综合。盆地模拟在计算生烃量方面技术已经比较成熟，因此刻度区（运聚单元）的生烃量仍由盆地模拟方法计算。

（3）油气参数研究。通过刻度区解剖，建立了参数评价体系和预测模型，获得了地质条件定量描述参数、量计算参数和经济评价参数，如运聚系数、丰度等关键参数。从刻度区获得的量与生油量之比可计算出运聚系数，刻度区的量与面积之比可获得单位面积的丰度，还可得到其他参数等。由于盆地内坳陷（凹陷）内各单元成藏条件差异，求得的参数是不同的，故细分若干运聚单元，求取不同单元的参数，这样用于类比区会更符合实际。

3.刻度区研究成果与应用

通过刻度区解剖研究，系统地获得运聚系数、油气丰度等多项关键参数，为油气评价提供各类评价单元类比参数选取的标准，保证评价结果科学合理。如中国石油解剖的辽河坳陷大民屯凹陷级刻度区，通过对其烃源条件、储层条件、圈闭条件、保存条件以及配套条件五方面精细研究，获得了22项量化的成藏条件的系统参数。根据大民屯凹陷内划分的六个运聚单元，分别计算各单元的生油量和量，直接获得六个单元的运聚系数。同时计算出各运聚单元单位面积的量，获得不同成藏条件下的丰度参数（表4-5）。

表4-5 大民屯凹陷刻度区解剖参数汇总表

在中国石油128个刻度区的基础上，各单位根据评价需要，又解剖了一定数量的刻度区。其中，中国石油利用已有刻度区128个，新解剖刻度区4个，共应用132个；中石化新解剖42个；中海油新解剖4个；延长油矿新解剖3个。各项目共应用了181刻度区，这些刻度区涵盖了我国主要含油气盆地中的大部分不同类型的坳陷、凹陷、运聚单元和区带，基本满足了不同评价区的需要。各种类型刻度区统计见表4-6。

表4-6 各种类型刻度区统计表

（二）有效烃源岩有机碳下限

有效烃源岩有机碳下限是指烃源岩中有机碳含量的最小值，小于该值的烃源岩生成的烃量不能形成有规模的油气聚集。有效烃源岩有机碳下限是确定烃源岩体积的主要参数，直接影响生烃量的计算结果。

在大量烃源岩样品分析化验和有关地质资料研究基础上，明确了不同岩类有效烃源岩有机碳下限标准。陆相泥岩有效烃源岩有机碳下限为0.8%，海相泥岩为0.5%，碳酸盐岩为0.2%～0.5%，煤系源岩为1.5%。例如，陆相泥岩TO C与S1+S2关系表明，S1+S2在TO C为0.8%时出现拐点，有效烃源岩有机碳下限定为0.8%；碳酸盐岩气源岩残余吸附气量与有机碳关系表明，残余吸附气量在有机碳为0.2%处出现拐点，有效烃源岩有机碳下限定为0.2%（图4-1、图4-2）。

图4-1 陆相泥岩TOC与S1+S2关系图

图4-2 碳酸盐岩气源岩残余吸附气量与有机碳关系图

对于勘探实践中已经发现油气藏，但烃源岩有机碳含量未达统一下限的盆地，根据实际情况可进行适当调整。如柴达木盆地柴西地区，在分析了大量烃源岩有机碳和S1+S2指标资料后，明确该区有机碳含量下限为0.4%时，即达到有效烃源岩标准，并被发现亿吨级尕斯库勒大油田的勘探实践所证实。在渤海湾盆地评价过程中，建立起相对统一的有效烃源岩丰度取值下限标准：碳酸盐岩气源岩丰度下限取0.2%，碳酸盐岩油源岩丰度下限取0.5%，湖相泥岩丰度下限取1.0%。

有效烃源岩有机碳下限的基本统一，保证了生烃量计算标准的相对一致和全国范围内的可比。

（三）产烃率图版

烃源岩产烃率图版是用盆地模拟方法计算烃源岩生烃量和量的关键参数。产烃率图版一般用烃源岩热模拟实验方法获得。

1.液态烃产率图版

利用密闭容器加水热模拟实验方法，对中国陆相盆地不同类型烃源岩进行了热模拟实验。模拟实验所用样品取自松辽、渤海湾等10个盆地，包括侏罗系、白垩系和古近系的湖相泥岩、煤系泥岩和煤3大类烃源岩。其中湖相泥岩烃源岩的有机质类型包括Ⅰ型、Ⅱ1型、Ⅱ2型和Ⅲ型，煤系泥岩烃源岩的有机质类型包括Ⅱ2型和Ⅲ型，煤烃源岩的有机质包括Ⅱ1型、Ⅱ2型和Ⅲ型。根据模拟实验结果，编制了不同类型烃源岩的液态烃产率图版（图4-3、图4-4、图4-5）。

图4-3 湖相泥岩烃源岩液态烃产率图版

图4-4 煤系泥岩烃源岩液态烃产率图版

图4-5 煤烃源岩液态烃产率图版

2.产气率图版

由于生物气生气机制与干酪根成气和原油热裂解气的生气机制不同，因此，其产气率与干酪根和原油裂解气产气率求取方式不同。

（1）生物气产气率。对生物气源岩样品在25℃～75℃的条件下进行细菌培养产生生物气，由此得到不同温阶下各类有机质的生物气产率。在模拟实验结果的基础上，结合前人的研究结果，分别建立了淡水环境、滨海环境和盐湖环境中不同类型有机质的生物气产气率图版及演化模式。

（2）干酪根和原油裂解气产气率。对于不同类型气源岩油产气率，国内外学者及一、二轮评价中已做过大量的工作。较多的实验是应用热压模拟方法对各种类型烃源岩进行产油及产气率实验，这种方法所计算的产气率包括了原油全部裂解成气的产率，亦即常说的封闭体系下源岩的产气率，所得到的天然气产率是气源岩的最大产气率。另一种求取气源岩产气率的方法是在开放体系下对源岩进行热模拟实验，各阶段生成的天然气和原油均全部排出源岩，原油不能在源岩中进一步裂解为天然气。这两种情况都是地质中的极端情况。但是实际的地质条件大多是半开放体系，在这种情况下，源岩生成的油既不能全部排出烃源岩，也不能完全滞留于源岩中。不同地质条件下亦即开放程度不同情况下源岩产气率如何计算？具体方法为：求得封闭和开放体系下相同类型源岩的产气率，将上述两种体系下的产气率图版（中值曲线）输入盆地模拟软件中，得出烃源岩层在不同渗透条件下产气率图版。

（四）运聚系数

运聚系数是油气聚集量占生烃量的比例，是成因法计算量的一个关键参数，直接影响量计算结果。运聚系数的确定方法包括运聚系数模型建立法和运聚单元成藏条件分析法。

1.运聚系数模型建立法

通过刻度区解剖，确定影响运聚系数的主要地质因素及其与运聚系数的相关关系。刻度区解剖研究表明，烃源岩的年龄、成熟度、上覆地层区域不整合的个数和运聚单元的圈闭面积系数等地质因素与石油运聚系数之间存在相关关系。依此建立地质因素与石油运聚系数之间关系的统计模型，包括双因素模型和多因素模型。双因素模型（相关系数为0.922）的地质因素选用烃源岩年龄和圈闭面积系数：

lny=1.62-0.0032x1+0.01696x4

多因素模型（相关系数为0.934）的地质因素选用烃源岩年龄、烃源岩的成熟度、区域不整合个数和圈闭面积系数：

lny=1.487-0.00318x1+0.186x2-0.112x3+0.02118x4

式中：y——运聚单元的石油运聚系数，%;

x1——烃源岩年龄，Ma;

x2——烃源岩成熟度（Ro）,%;

x3——不整合面个数；

x4——圈闭面积系数，%。

2.运聚单元成藏条件分析法

依据刻度区提供的大量运聚系数，依盆地类型和影响运聚系数的主要地质因素，分类建立运聚系数取值标准与应用条件。在评价中，根据刻度区解剖结果，确定了油气运聚系数分级取值标准（表4-7）。在评价中得到了推广应用，取得了良好的效果。

表4-7 石油运聚系数分级评价表

（五）最小油气田规模

最小油气田规模是指在现有工艺技术和经济条件下开地下，当预测达到盈亏平衡点时的油气田可储量。最小油气田规模对统计法计算的量结果有较大影响。为此，中国石油天然气集团公司等三大石油公司和延长油矿管理局对最小油田规模进行了专门研究。

通过对不同油价、不同开发方式和未来可能技术条件下最小油气田规模研究，确定了不同地区的最小油气田规模的取值。在地理环境相对较好的东部地区，其勘探开发成本较低，最小油气田规模一般在10×104～30×104t，在地理环境相对较差的西部地区，其勘探开发成本高，最小油气田规模一般在50×104t以上，对于海域来说，油气勘探开发成本更高，最小油气田规模更大，一般在150×104～500×104t。

（六）丰度

油气丰度是指每平方公里内的油气量，是类比法计算量的关键参数。通过统计分析，建立了丰度模型和取值标准。

1.丰度模型

通过刻度区解剖，建立刻度区内评价单元油气丰度和相关地质要素之间的统计预测模型：

新一轮全国油气评价

式中：y——运聚单元的石油丰度，104t/km2;

x1——烃源岩生烃强度，104t/km2;

x2——储集层厚度/沉积岩厚度，小数；

x3——圈闭面积系数，%;

x4——不整合面个数。

2.丰度取值标准

通过统计不同含油气单元丰度的分布特点，结合地质成藏条件，总结出各类刻度区丰度的取值标准。

（1）不同层系丰度：古近系凹陷由于成藏条件优越，成藏时间晚，石油地质丰度一般大于20×104t/km2；中生代凹陷成藏时间相对较长，石油地质丰度相对较低，一般约为10×104t/km2；古生代凹陷由于生、储层时代老，多期成藏多期改造、破坏，预计其丰度更低。

（2）不同类型运聚单元丰度：中新生代断陷或坳陷盆地长垣型、潜山型和断陷型中央背斜构造型，石油地质丰度高，一般大于40×104t/km2；中新生代裂陷盆地、坳陷盆地边缘构造型和古近系缓坡构造型石油丰度次之，一般为10×104～30×104t/km2；中生代盆地岩性型和古生代压陷盆地的构造型石油丰度相对较低，一般小于10×104t/km2。

（3）不同区块或区带级丰度：区块或区带级石油丰度差异更大，从小于1×104t/km2到大于200×104t/km2。其中潜山型、岩性—构造型、披覆背斜区块丰度较高，一般大于50×104t/km2，最大可大于200×104t/km2。构造—岩性型、断裂构造型丰度一般为30×104～50×104t/km2。地层—岩性型、断鼻型以及裂缝型区块、丰度较低，一般小于30×104t/km2。

通过刻度区解剖标定多种成藏因素下评价单元的丰度，不但为广泛应用类比法计算量提供了可靠的参数，同时也摆脱了过去以盆地总量为基础，利用地质评价系数类比将量分配到各评价单元的做法，使类比法预测的油气量在空间位置上更准确，提高了油气空间分布的预测水平。

（七）可系数

国外主要用建立在类比基础上的统计法计算油气可量，而我国第一轮、第二轮全国油气评价没有计算油气可量。本轮评价开展的油气可系数研究，通过可系数将地质量转化为可量，这在国内外油气评价中尚属首次。可系数是指地质中可出的量占地质量的比例，是从地质量计算可量的关键参数。

可系数研究与应用是常规油气评价的重要组成部分，主要目的是通过重点解剖、统计和类析方法，对我国油气可系数进行研究，为科学合理地计算油气可量提供依据，进而对重点盆地和全国油气可潜力进行评价。

1.评价单元类型划分

为使可系数研究成果与评价单元划分体系有机结合，遵循分类科学性、概括性和实用性三个基本原则，以油气类型、盆地类型、圈闭类型、储层岩性、储层物性等地质因素为依据，对评价单元进行了分析和分类，将国内石油评价单元分为中生代坳陷高渗、古近纪与新近纪断陷盆地复杂断块高渗等24种类型，天然气评价单元分为克拉通盆地古隆起、前陆盆地冲断带等16种类型（表4-8、表4-9）。

表4-8 不同类型评价单元石油可系数取值标准

表4-9 不同类型评价单元天然气可系数取值标准

2.刻度油气藏数据库的建立

已发现油气赋存在油气藏中，建立刻度油气藏数据库是统计已发现油气收率、分析影响收率主控因素、预测油气可系数的基础。刻度油气藏是油气可系数研究中作为类比标准的，地质认识清楚、开发程度高、已实施二次油或三次油技术的油气藏。

刻度油气藏选择原则：①典型性——能代表国内外主要的油气藏类型，保证类比法应用基础的广泛性；②针对性和实用性——针对油气评价，有效地指导相应类型评价单元油气可系数的确定；③开发程度高——油气藏开发程度高，地质参数和开发参数基本齐全；④三次油技术应用具有代表性——尽量选择已实施三次油技术的油藏，保证技术可系数的可靠性。

对国内43个油藏、30个气藏，国外59个油藏、22个气藏进行了剖析：收集整理每个油气藏的主要地质和开发参数；每个油气藏的地质条件主要包括储层特征、圈闭条件、流体性质等，开发条件主要包括开方式、开速度、增产措施等；研究不同因素对收率的影响程度，进而确定该油气藏收率的主控因素；针对开方式的不同，油藏的收率可分为一次、二次或三次收率；气藏主要是一次收率。通过对每个油气藏的地质条件、开发条件和收率进行分析，建立起国内外刻度油气藏数据库。

3.可系数主控因素分析

对影响可系数的地质条件、开发条件和经济条件进行了分析，建立起可系数主控因素的评价模型。

（1）在大量统计和重点解剖的基础上，对油气地质条件中的因素逐一进行分析，并提炼出15项油气收率的主控因素，即盆地类型、储层时代、圈闭类型、沉积相类型、储层岩性、储层厚度、储集空间类型、孔隙度、渗透率、埋深、含油饱和度、原油粘度、原油密度、变异系数、原始气油比。

（2）在诸多开发条件中，提高收率技术是极为重要的因素，不同提高收率技术适用条件不同，其提高收率的潜力也差距很大。通过综合分析，主要技术对不同类型油藏的提高收率潜力为：最小5%，中间值10%，最大值15%。

（3）利用石油公司提高收率模拟研究成果，建立了大型背斜油藏、复杂背斜油藏、断块油藏、岩性油藏、复杂储层油藏等在税后内部收益率为12%、油田开发到含水95%时聚合物驱和化学复合驱油时的油价与油田收率之间的关系，若这五类油藏要达到相同的收率，条件好的如大型背斜油藏、复杂背斜油藏所需的油价低于条件差的如岩性油藏、复杂储层油藏。

4.可系数取值标准的建立

在研究中，解剖了国内43个油藏、30个气藏，国外59个油藏、22个气藏，统计分析了大量油气田收率数据，给出了不同类型评价单元油气技术可系数和经济可系数取值范围，建立了不同类型评价单元油气可系数取值标准（表4-8、表4-9）。

（1）不同类型评价单元石油可系数相差较大，以技术可系数为例：中生代坳陷高渗和古近纪与新近纪断陷盆地复杂断块高渗评价单元可系数最大，其中间值大于40%；中生代坳陷中渗、古近纪与新近纪断陷盆地复杂断块中渗、中生代断陷、中新生代前陆、古生界潜山、古生界碎屑岩、古近纪残留型断陷、陆缘裂谷断陷古近纪与新近纪海相轻质油、陆缘弧后古近纪与新近纪海陆交互相轻质油等评价单元可系数为30%～40%；中生代坳陷低渗、古近纪与新近纪断陷盆地复杂断块低渗、古生界缝洞、南方古近纪与新近纪中小盆地、低渗碎屑岩、重（稠）油中高渗、变质岩、砾岩、陆内裂谷断陷新近纪重质油、陆内裂谷断陷古近纪复杂断块等评价单元可系数为20%～30%；低渗碳酸盐岩、重（稠）油低渗、火山岩等评价单元可系数为15%～20%。

（2）不同类型评价单元天然气可系数相差也较大：克拉通碳酸盐缝洞、礁滩和前陆冲断带等评价单元可系数最大，其平均值大于70%；克拉通古隆起、克拉通碎屑岩、前陆前渊、南方中小盆地、陆缘断陷、火山岩、变质岩和海域古近纪与新近纪砂岩等评价单元可系数为60%～70%；前陆斜坡、生物气、中生代坳陷、古近纪与新近纪断陷盆地复杂断块、残留断陷、砾岩等评价单元可系数为50%～60%；致密砂岩等评价单元可系数最小，其平均值小于50%。

5.可系数计算方法的建立

可系数计算方法包括可系数标准表法和刻度区类比法两种方法。

（1）标准表取值法。利用可系数标准表求取不同评价单元可系数的步骤如下：在不同类型评价单元可系数取值标准表中找到已知评价单元的所属类型；明确评价单元与可系数相关因素（宏观、微观）的定性、定量资料；对照可系数的类比评分标准表和类比评分计算方法，对评价单元进行类比打分；根据类比评价结果求取可系数。

（2）刻度区类比法。以建立的国内外刻度油气藏数据库为基础，利用刻度区类比法来求取不同评价单元的可系数。具体步骤如下：根据评价单元分类标准，将具体评价单元归类，并分析整理该评价单元的油气地质条件和开发条件；根据评价单元的类型及其地质条件和开发条件，从国内外刻度油气藏数据库选择适合的类比对象；对照可系数的类比评分标准表和类比评分计算方法，对该评价单元及其类比对象进行打分并计算它们的得分差值；根据得分差值求取该评价单元的可系数。

通过油气可系数标准和计算方法在全国129个盆地中的推广应用，既检验了可系数取值标准和所用基础数据的可靠性、可行性和适用性，保证了油气可量计算的客观性，又获得了全国油气可量。

（一）动态评价的必要性

随着地质认识和勘探开发形势的不断变化，对油气的认识也会不断更新，需要开展经常性评价，实现评价系统化、制度化、动态化，为制定能源政策和编制国家能源中长期发展规划，提供重要的科学依据。

（二）动态评价的成果要求

油气动态评价仍取组织、专家技术把关、石油公司和大学及研究院所具体承担实施任务的工作思路，继续实行产、学、研相结合的工作方式，集中优势力量完成评价工作。

全国油气动态评价工作是我国第一次由部门组织的以国家利益为主、突出国家需要的公益性和基础性的油气动态评价工作。评价成果还要体现及时性和前瞻性，分层次进行，既有及时跟踪勘探新进展、地质新认识，预测勘探新趋势，客观选择评价盆地和评价领域的年度动态评价；又包括对新区、新领域油气的调查评价。

油气年度动态评价项目为国家层面经常性的油气调查评价专项，每年滚动进行，每个五年的前一年进行一次全面评价，评价周期为每年的7月1日到下年的6月30日。第一轮初步为5年，即2006～2010年。

（三）动态评价完成的目标

进行动态评价的油气类型为常规和非常规油气，需要完成的主要目标有：

（1）及时跟踪盆地勘探进展，判断勘探开发趋势，探索新区、新领域、新层系油气前景，总结地质新认识，确定每年动态评价重点目标；开展油气年度动态评价，获得重点评价目标的油气地质量、可量数据，分析油气品质状况。

（2）总结油气分布规律，分析油气储量、产量增长趋势和开发利用前景。

（3）完善油气动态评价体系，包括油气评价方法、参数体系、评价规范、评价流程等。完善国家油气评价系统。对油气进行动态管理，实现油气评价工作的信息化、规范化、制度化，最终实现全国油气的动态评价。