珠江口盆地恩平凹陷断层封闭性评价

许新明，姜　建，陈胜红，刘丽华，刘贤来

(中海石油(中国)有限公司 深圳分公司，广东 广州　510240)



珠江口盆地恩平凹陷断层封闭性评价

许新明，姜建，陈胜红，刘丽华，刘贤来

(中海石油(中国)有限公司 深圳分公司，广东 广州510240)

为了综合评价断层的封闭性，以珠江口盆地恩平凹陷新近系为例，通过实钻断层圈闭油藏的封堵规律研究，标定断裂带泥质含量、泥岩涂抹系数和断面正压力临界值，建立了适合勘探实践的断层封闭性评价的新方法。结果表明：断层油藏与断层涂抹量呈现出泥岩涂抹因子随深度逐渐递增的趋势；断裂带泥质含量与泥岩发育情况是断层能够封堵的先决条件，断层的封闭性受构造图精度、构造调整和控圈断裂断层活动影响。应用该方法能有效评价未钻断层圈闭油气保存条件，有助于提高断层圈闭的勘探成功率，可推广至具有相似构造成因和沉积背景的勘探新区进行断层封闭性分析，是勘探早期评价断层封闭性有效的方法。

珠江口盆地；恩平凹陷；断层封闭；评价体系；新近系

0　引　言

断层封闭性的评价方法很多，比较常用的方法归纳有以下5类：(1)分析岩性侧向并置的方法，如Allen图解、Knipe图解、断移地层砂泥比值法和砂泥对接概率法等[1-4]；(2)分析断裂带侧向泥质涂抹能力的方法，如SSF、SGR、CSP法等[5-7]；(3)分析断层岩与围岩压力差的方法，如断-储排替压力差法、两侧流体相关参数对比法等[8-9]；(4)分析断裂垂向封闭性的方法，如断面正压力法[10-14]；(5)多因素综合判识方法，如断层岩破坏强度法、断裂开启指数法、模糊综合评判法[9,15]。但每种方法都有其局限性，单一的方法难以客观地评价不同地区断层的封堵风险。为了综合评价断层的封闭性，以恩平凹陷为例，建立了适合该区勘探实践的断层封闭性评价方法，即：首先对于勘探程度较低、探井少或无探井的地区，借助地震资料来预测断裂带的泥质含量，运用速度-岩性转换模版预测未钻圈闭的封闭性；其次，通过实钻资料厘定恩平凹陷泥岩涂抹系数SSF的门限值，运用SSF法可预测断层的侧向性；再次，分析并总结出恩平凹陷泥岩发生塑性变形的门限值，计算断层现今断面的正压力，通过比较该值与泥岩发生塑性变形的门限值，判别断层的垂向封闭性。

1　区域地质背景

图1　恩平凹陷新近系断裂分布特征Fig.1　Distributing features of the Neogene structure in Enping sag

恩平凹陷位于珠江口盆地珠一坳陷西南部，长期处在南东-近南北向拉张应力场作用下，新近系以来形成了北西西、近东西向2组张扭性断层[16-18]，其中以正断层为主，正、反向正断层共存，并发育了大量断层圈闭(图1)。目前该区仅在新近系钻遇了一定规模的油气，发现了A8、A9等油田，新近系主要目的层为韩江组下段、珠江组上段和珠江组下段，砂岩质量分数介于30%～50%，对应地震顶面标志层分别为T35、T40 和T50。该区石油地质探明率约13%，油气勘探仍具有较大潜力[17]。

研究区新近系断层主要以生长断层为主，断层圈闭发育，多数断层在油气大规模排烃期(10 Ma)以后仍然活动，在油气运聚与成藏中具有通道和封堵的双重作用[16]，此类断层封闭性与断裂带泥质含量、断裂带泥岩涂抹能力以及断面的紧闭程度密切相关，断层封堵与否以及封堵能力的大小是制约该地区新近系断层圈闭成藏的关键。

2　断层封闭性评价方法

断层不是一个简单的面，而是一个三维地质体，具有复杂的内部结构[19-21]。分析断层是否封堵，需半定量地计算断裂带内部泥质含量，并分析泥岩涂抹或岩性对接所形成的侧向封堵程度，还需分析断层面的紧闭程度所形成的垂向封堵能力。

2.1建立速度-岩性转换模版预测断层侧向封闭性的方法

对于勘探程度较低、探井较少的地区，可利用已钻井的声波时差和地震速度谱资料来预测断裂带的泥质含量，并通过已证实油藏所断移地层的砂泥岩比值来厘定达到封堵门限的断裂带泥质含量。具体原理如下[22-23]：(1)根据研究区已钻测井资料的声波时差曲线，获取大量的砂岩、泥岩的速度数据，对速度数据统计分析得出纯泥岩和纯砂岩的速度随深度的变化规律，根据纯泥岩和纯砂岩的速度随深度的变化规律做出该区砂泥岩的压实模型；(2)从速度谱上读取断层断开层段的地震反射时间t0对应的叠加速度Vs，依据Dix公式求得相应层段的地震层速度Vint；对其相邻的多个速度谱点的层速度进行平滑处理，将计算所得的地震层速度与测井声波时差计算的速度值进行比较校正，建立井震速度(测井声波时差计算的速度值与地震层速度)的可比性模型；(3)结合声波时差计算的速度所建立的砂泥岩的压实模型，建立“纯泥岩和纯砂岩的速度随深度的变化规律-地层埋深-地震层速度”的关系图版，通过岩石体积物理模型可以建立岩石指数转换模型：

(1)

从而得到：

(2)

式中：Vint为地震层速度；Vs和Vm分别为纯砂岩和纯泥岩的速度，此二者可由压实模型提供；Ps为单层砂岩百分比。由此可计算出单层i的泥岩百分含量Rmi=1-Psi，将获得的单层砂泥岩百分含量数据点(Rmi和Psi) 投到已建立的“砂泥岩百分含量—地震层速度—地层埋深”的关系图版上，可获取断层断开该层的地层埋深，通过计算断层断开各地层的埋深可求得断层断开地层的厚度，从而求取Σhi，由该公式可求得Rm。(4)通过上述方法，得到研究区恩平凹陷速度-岩性转换模版(图2)。通过计算实钻油藏(断层封堵成藏)所在深度的层速度，由井震速度的可比性模型，可以投影至该区砂泥岩的压实模型坐标系中，即恩平凹陷速度-岩性转换模版中。然后，分析实钻油藏计算的断裂带泥质含量分布状况，厘定断裂能够封堵时断裂带泥质含量临界线。分析认为：当地震速度谱预测目的层段断裂带的泥质含量大于该临界线时，认为断裂侧向封闭性良好，反之则反。

图2　恩平凹陷速度-岩性转换模版Fig.2　Velocity and lithology translation template of Enping sag

2.2建立泥岩涂抹临界值预测断层侧向封闭性的方法

泥岩或非储集层物质滑动注入断裂带内，形成对砂层“涂抹”的不渗透泥质带，泥质含量越高，涂抹带在空间上连续性越好，易封堵。在珠江口盆地，目前主要运用LINDSAY提出的泥岩涂抹因子(SSF)指标[6,23-24]来评价断层侧向封闭性，断层的断距越小，断开泥岩的厚度越大，泥岩涂抹层在空间上的连续性越好，泥岩涂抹因子的计算公式为：

(3)

断层泥质涂抹量属性仅是断层面形成的泥岩涂抹的相对可能估算值；断层能否封闭，门限值如何，需根据钻后已证实封堵的断层圈闭油藏反推的泥岩涂抹因子门限值来标定。根据恩平地区勘探实践，计算了恩平南部翘倾断裂构造带大量实钻断层圈闭油藏达到泥岩涂抹效果时，对应的泥质涂抹系数数据集，结果表明：断层成藏与断层涂抹量存在一定的统计相关性，呈现出泥岩涂抹因子随深度逐渐递增的趋势，与成岩作用使岩石的孔隙度总体趋势减小、深层比浅层易封堵有着良好的对应关系(图3)。

据此，建立了恩平凹陷新近系目的层段预测断层封闭性泥岩涂抹因子的取值规律(见公式(4))，分析认为：对于具有相似区域地质条件的断层圈闭，预测其不同层段断层的侧向封闭性，SSF的取值规律具有较好的应用价值。

(4)

式中：SSF为断层的泥岩涂抹因子，无量纲；H为断面的地层埋深，m。

图3　A9油田钻遇断层油藏所需泥岩涂抹因子门限值随深度的变化关系Fig.3　Threshold value of shale smear factor in drilled faulted reservoirs changing with depth in A9 Oilfield

2.3建立现今断面正压力临界值曲线并预测断层垂向封闭性的方法

影响断层垂向封闭性好坏的最关键因素是断层面的紧闭程度。当断裂活动停止后，断层面受上覆地层的重力作用，发生紧闭愈合，断层面所受压力越大，断层紧闭愈合程度越好，易形成垂向封闭[9-13]。

2.3.1现今断面正压力计算原理

断层面的紧闭程度可用断层面所受正压力大小来衡量,其计算公式[11]为：

P=0.009 876H(ρr-ρw)cosθ

(5)

式中：P为断面所受上覆地层闭合正压力，MPa；H为断面埋深，m；ρr为上覆地层的平均密度，取2.33 t/m3；ρw为地层水的密度，取1.03 t/m3；θ为断层倾角，(°)。

当断面所受正压力大于泥岩的变形强度时，泥质成分发生塑性变形而堵塞断裂带的渗漏空间，断面必将合拢封闭。泥岩塑性变形极限值可根据断层所处的埋深由下列公式[10]来确定：

Pshale=248.17H-207.47

(6)

式中：Pshale为泥岩塑性变形极限值，MPa；H为断层面埋深，m。

2.3.2断层垂向封闭临界值的确定

图4　已知含油构造主要目的层段断面正压力Fig.4　Normal pressure on fault plane of the objective interval in oil-bearing structure

本次研究采用与断层有关的测井资料、地震资料和油藏资料来计算4个已钻含油构造T35、T40、T50、T60对应的断面正压力值(T60为珠海组地震顶面标志层)，厘定预测全区断层面正压力临界值。该方法的前提是已钻构造的控圈断层具有相同的力学形成机制，且地震反射界面断面正压力的平均值作为该地震反射界面能够封闭的临界值。通过对凹陷内已知含油构造多个地震反射界面断面正压力的计算与统计(图4)，对断块油藏而言，其地震反射界面所计算的断面正压力是大于泥岩发生塑性变形所需的最大压力值，泥岩发生塑性变形使断裂带垂向封闭，断块构造油气能够成藏。因此，可以应用各地震反射界面计算的断面正压力值，取其平均值来预测恩平凹陷各层段的断层垂向封闭性，在韩江组下段、珠江组上段、珠江组下段、珠海组分别为6.1 MPa、7.6 MPa、9.9 MPa、11.1 MPa。当未钻构造目的层位计算的断面正压力小于该临界值时，认为断面不闭合，存在油气泄漏点，断层垂向不封堵。

3　评价标准与适用性

在珠江口盆地恩平地区断层圈闭风险评估的勘探实践中，证实该断层封闭评价体系具有运用基础资料多、量化后预测精度较高的特点，但对断层封堵的定量评价是相对的封堵，常受资料精度的制约(表1)。

速度-岩性转换模版预测的Rm可以反映断裂填充物中泥质含量的大小，但由于速度谱的选择存在一定的误差，因此它不是断裂填充物中真正的泥质百分含量，只是一个大致相当的预测值[22-23]，且该值往往介于围岩地层泥地比和断层错断地层泥质含量理论值之间[23]，反映了断层岩内泥质含量不仅受错断围岩砂地比影响，并且不是按照错断地层岩石成分均匀地充填于断裂带内部。

泥岩涂抹法的精度受断层圈闭的构造图、邻井的岩性数据以及SSF门限值的制约。实践证实，构造图对断层上下盘的地层在断面上的相对位置、形态有显著的影响，造成钻后在技术上的对接封堵与实钻完全不封堵的差异[23]。断层泥质涂抹量属性仅是断层面形成的泥岩涂抹的相对可能估算值，断层能否封闭，门限值如何，需根据钻后已证实封堵油藏反推的SSF门槛值与断裂带泥质含量综合标定。同时，要注意到，该定量计算方法具备诸多的不确定性，譬如符合封闭门限的层段也不一定封堵。因此，SSF的门限值在定量评价过程中可在相似构造条件、沉积环境下的地区推广使用。

断面压力法的精度受地震资料品质和断裂属性的影响，具有相似区域应力和力学形成机制的断裂才能由已知断面压力临界值来预测断面是否紧闭。通常断层平面展布平直的断层，如A10、A12的控圈断层，断面紧闭性较好；而具有扭动或雁行式发育的断层如A8、A9的控圈断层(图1)，断面挠曲或分叉处断面不紧闭而存在垂向封闭漏点的概率要大得多，实际工作中，应加密求取该区域断面正压力的数值，谨慎评估该类断层的垂向封闭风险。

表1　断层封闭性评价体系的评价标准

图5　恩平凹陷A14构造主要目的层段断裂带泥质含量预测Fig.5　Predicting the shale content of fault rock in the objective interval of A14 trap in Enping sag

4　实例应用与效果

4.1断层圈闭钻前风险评估实例

A14构造位于恩平南部翘倾断裂构造带西段(图1)，为一单条断层控圈的断鼻构造，断层以中晚期活动为主且断穿粤海组、韩江组、珠江组、珠海组和恩平组地层，在韩江组下段至珠海组地层(对应地震反射层位为T35—T60)均发育断层圈闭，勘探主要目的层位是韩江组下段和珠江组。断层是否封堵是该构造油气成藏的关键因素，因此钻前应用断层封闭性评价体系对该构造断层封堵风险进行系统的分析，并综合评估断层圈闭的勘探风险。

图6　A14构造目的层段岩性对接与泥岩涂抹评价结果Fig.6　Evaluation results of lithology juxtaposition and shale smearing in the objective intervals of A14 trap

首先，分析A14断裂带的速度谱资料，选取合理的叠加速度并计算为层速度，将各目的层段断裂带部位的层速度投影至恩平凹陷速度-岩性转换模版中(图5)，得出以下认识：(1)A14构造在韩江组下段预测的断裂带泥质含量介于56%～63%,断层在相应地层埋深处断裂带泥质含量都大于临界线,判断该构造在韩江组下段(介于地震反射层T35—T40之间)具备较好的侧向封闭性；(2)A14构造在预测的断裂带泥质含量介于49%～56%，断层在相应地层埋深处断裂带泥质含量基本略小于临界线,但已钻含油构造A10、A12、A13(图2)部分断裂带泥质含量也存在类似现象，仍然有少量油藏为断层封堵成藏，判断该构造在珠江组(介于地震反射层T40—T60之间)地层断层封堵存在一定风险。

其次，应用泥岩涂抹因子算法计算断层泥岩涂抹能力。此时，为了便于展示该断层各目的层段断层侧向封闭性，选取地震标志层位T35、T40、SB17.5、T50进行断层封堵分析，应用公式(4)计算SSF取值分别为1.67，1.97，2.1，2.5，运用断层封堵分析评价系统泥岩涂抹因子算法软件，得到A14构造各目的层段断层两盘岩性对接封堵评价结果和岩性与泥岩涂抹评价结果(图6)。结果表明：韩江组下段、珠江组地层泥岩涂抹效果较理想，预测断层在韩江组下段、珠江组都有较好的断层封堵能力，油气保存条件良好。

最后，根据公式(5)计算出A14构造控圈断层现今条件下不同位置、不同层段的断面正压力，并将计算值与现今断面正压力临界值进行对比，判断该断层断面是否存在垂向封闭“漏点”，从而预测A14构造控圈断层断面的垂向封堵能力(图7)。结果表明：(1)断层在同一层位不同位置所受断面正压力具有较大差异；在联络测线上，珠海组和珠江组下段呈现出由东南往西北现今断面正压力逐渐增大的趋势，而珠江组上段、韩江组下段呈现出由东南往西北现今断面正压力逐渐减小的趋势，由公式(5)可知断面正压力与断层倾角成反比，说明A14构造控圈断层断面在空间上呈现出中浅层断面扭转的形态，这与该区长期处在张扭应力背景吻合，而张扭断层的断面易于闭合，有利于断层垂向封堵。(2)该断层在主要目的层段计算的现今断面正压力均大于恩平凹陷现今断面正压力临界值，反映了该断层断面闭合程度较好，预测不存在垂向“漏点”，断层垂向封堵效果很好。

4.2应用效果

钻前预测在韩江组的断裂带泥质含量均大于临界线(图5)，实钻结果发现油层集中在韩江组，而在珠江组，预测的断裂带泥质含量低于临界值，虽然珠江组钻遇了大量荧光显示段，测井解释多为含油水层,未发现一个油层。实钻结果与断裂带泥质含量的高低有较好的对应关系，说明断裂带泥质含量与泥岩的发育情况是断层能够封堵的先决条件。以韩江组下段为例(T35-T40)，钻前和钻后均选取泥岩涂抹系数计算的参考取值为1.67，计算的涂抹效果仍然高估了实际钻探的封堵效果，在泥岩涂抹SSF法预测断层的侧向封闭性上，这种状况非常普遍。即使在断面紧闭,断层垂向封堵的情况下(图7),运用已知油藏反推较合理的SSF值，计算得到的预测含油砂层数仍然可能多于实钻钻遇的油层数，表明油气侧向封闭预测效果受构造图精度、古构造调整对古油藏破坏程度和控圈断裂晚期活动对油气调整作用等的影响(表2)。因此，断层封闭性评价体系对未钻断层圈闭油气保存评价具有较好的预测效果，但仍需结合多重因素评估断层封闭性的风险。

图7　A14构造目的层段断层现今垂向封闭性Fig.7　The evaluation results of vertical seal in the objective interval of A14 trap

Table 2Comparative analysis of fault sealing between before and after drilled in the fault block traps

层段钻前预测封堵效果速度-岩性转换模板法泥岩涂抹法 断面压力法 断层封闭性评价体系综合评定 钻探效果韩江组下段 较好一般较好较好4个油层珠江组上段 较差较好较好一般多个含油水层,大量气测显示珠江组下段 较差较好较好一般若干含油水层,大量气测显示

5　结　论

(1)对于勘探程度较低的地区而言，可以充分利用速度谱资料，运用速度-岩性转换模版预测断裂带的泥质含量，进而对断层封闭性进行初步的评价，方法可行性好。

(2)恩平凹陷断层成藏与断层涂抹量存在一定的统计相关性，呈现出泥岩涂抹因子随深度逐渐递增的趋势，与成岩作用使岩石的孔隙度总体趋势减小、深层比浅层易封堵有着良好的对应关系。

(3)断裂带泥质含量与泥岩的发育情况是断层能够封堵的先决条件，泥岩涂抹法在预测断层的侧向封闭性上，常高估实际钻探的封堵效果，应充分考虑构造图精度、古构造调整、控圈断裂晚期活动等因素对油气保存的影响。

(4)运用该方法系统评价未钻圈闭的断层封堵风险，能够从垂向和侧向封闭性的角度，对圈闭的断层封闭风险进行综合分析，有助于提高断层圈闭的勘探成功率；对具有相似构造成因和沉积背景的断层进行封闭性分析，具有良好的推广价值。

[1]ALLAN U S.Model for hydrocarbon migration and entrapment within faulted structures[J]. AAPG Bulletin,1989,73( 7):803-811.

[2]KNIPE R J.Juxtaposition and seal diagrams to help analyze fault seals in hydrocarbon reservoirs [J].AAPG Bulletin,1997,81(2):187-195.

[3]赵蜜福，李阳，张煜，等.断层两盘岩性配置关系及断层的封闭性[J].中国石油大学学报(自然科学版)，2006，30(1)：7-11.

[4]孙龙涛，陈长民，詹文欢，等.珠江口盆地断层封堵特征及其影响因素[J]. 石油学报，2007，28(4)：36-40.

[5]付广，孙同文，吕延防.南堡凹陷断-砂配置侧向输导油气能力评价方法[J].中国矿业大学学报，2014，43(1):79-87.

[6]LINDSAY N G, MURPHY F C, WALSH J J. Outcrop studies of shale mear on fault surface[J]. International Association of Sedimentologists Special Publication, 1993, 15(1): 113-123.

[7]王蛟，束青林.孤岛油田西南缘馆下段断层封闭性与油气成藏[J].西南石油大学学报(自然科学版)，2010，32(5)：47-53.

[8]刘哲，吕延防，付晓飞，等.贝尔凹陷断层侧向封堵能力定量研究[J].吉林大学学报(地球科学版)，2012，42(2)：353-361.

[9]付广，史集建，吕延防，等.断层侧向封闭性定量研究方法的改进[J].石油学报，2012，33(3)：414-418.

[10]王秋菊.断层垂向封闭所需断层面正压力下限的确定方法[J].大庆石油地质与开发，2008，27(4)：30-34.

[11]刘震，谭卓，蔡东升，等.用断层面正压力法分析北部湾盆地涠西南凹陷断层垂向封闭性及其演化[J].地质科学，2008，43(4)：695-711.

[12]梁全胜，刘震，何小胡，等.断层垂向封闭性定量研究方法及其在准噶尔盆地白家海凸起东道海子断裂带应用[J]. 现代地质，2008，22(5)：803-809.

[13]侯亚伟，杨庆红，黄凯，等.基于构造稳定性的断层垂向封闭能力评价 [J]. 断块油气田，2013，20(2)：166-169.

[14]姜大朋，代一丁，刘丽华，等.断裂输导油气的机制及侧向分流控制因素探讨 [J].现代地质，2014，28(5)：1023-1031.

[15]吕延防，王有功，付广，等.珠江口盆地珠一坳陷断层圈闭钻探风险性评价[J].石油学报，2011，32(1)：95-100.

[16]许新明，陈胜红，王福国，等.珠江口盆地恩平凹陷断层特征及其对新近系油气成藏的影响[J].现代地质，2014，28(3)：543-550.

[17]许新明，刘丽华，陈胜红，等.珠江口盆地恩平凹陷新近系油气成藏主控因素分析[J].地质科技情报，2015，34(1)：100-106.

[18]许新明，刘贤来，陈胜红，等.张扭性断陷盆地构造样式与油气成藏的关系——以珠江口盆地恩平凹陷新生界为例[J].海洋地质前沿，2014，31(1)：31-36.

[19]孟令东，付晓飞，王雅春，等.徐家围子断陷火山岩断层带内部结构与封闭性[J].石油勘探与开发，2014，41(2)：150-157.

[20]罗胜元，何生，王浩.断层内部结构及其对封闭性的影响[J].地球科学进展，2012，27(2)：154-164.

[21]BRETAN P，YIELDING G，JONES H. Using calibrate shale gouge ratio to estimate hydrocarbon column heights[J]. AAPG Bulletin, 2003, 87(3):397-413.

[22]马洪，柳广弟，朱立林，等.地震资料在断层封闭性评价中的应用与探讨[J].勘探地球物理进展，2008，31(1)：48-52.

[23]许新明.基于断层岩泥质质量分数预测的断层侧向封闭性评价[J].断块油气田，2014，21(5)：574-578.

[24]付广，夏云清.断层对接型和断层岩型侧向封闭的差异性 [J].天然气工业，2013，33(10)：11-17.

Study on the Evaluation Method of Fault Sealing and Its Application:An Example of Neogene Layer in Enping Sag, Pearl River Mouth Basin

XU Xin-ming, JIANG Jian, CHEN Sheng-hong, LIU Li-hua, LIU Xian-lai

(ShenzhenBranchofCNOOCLtd.,Guangzhou,Guangdong510240,China)

Taking the Neogene layers of the Enping sag as an example, a new way to evaluate fault sealing of structural closure is introduced under the current exploration practice to calibrate critical values, that is, obtaining clay content in fault zone, shale smear value and normal pressure of fault surface from sealing rules of faulted reservoir. The results show that the threshold values of shale smear factor in drilled faulted reservoirs gradually increases with depth. Development of shale and its content of fault zone are prerequisites for fault sealing, and the evaluation of fault sealing is affected by structure map accuracy, structural adjustment and fault activity. The risk of fault sealing, including fault’s vertical and lateral sealing, could be analyzed by the evaluation method, so as to improve the success ratio of fault traps exploration. And the evaluation method could be applied to the new exploration area with the similar tectonic genesis and sedimentary background as Enping sag, which is an effective methods to evaluate fault sealing during the early exploration periods.

Pearl River Mouth Basin; Enping sag; fault sealing; evaluation method; Neogene

2014-12-20；改回日期：2015-04-10；责任编辑：孙义梅。

国家“十二五”科技重大专项“大型油气田及煤层气开发”子课题“南海东部海域已证实的富生烃凹陷再评价及新领域勘探方向”(2011ZX05023-001-006)。

许新明，男，硕士，工程师，1984年出生, 海洋地质学专业，主要从事珠江口盆地(东部)地区油气勘探与成藏方面的研究工作。Email：xuxm3@cnooc.com.cn。

TE122.3

A

1000-8527(2016)01-0122-08