深水探井转开发井全过程环空压力预测与干预研究

南海深水油气资源极为丰富，占南海油气总资源量超过70%。目前我国在南海深水区域的勘探已获得了较大成果[1-3]，但深水探井单井作业成本可达数亿人民币，而深水开发井作业成本更高，高经济门限成为限制深水油气资源开发的关键因素，深水勘探开发一体化策略则是缓解这一问题的有效途径之一[4-5]。深水探井转开发井技术是深水勘探开发一体化的重要一环，探井转开发井过程中环空压力会经历多个井筒状态，包括探井钻井阶段、探井弃井后长期闲置阶段、转开发井完井重入阶段、生产阶段，井筒在各阶段所处的工况条件和温压环境不同，环空压力变化较大。

L-X气田开发是我国首个自营深水开发项目，为降低开发成本，将该区域若干口前期探井转化为开发井，针对深水气井井筒和油藏特征，考虑深水探井转开发井各阶段的热量传递模式，建立井筒环空温度和压力计算模型，实现深水探井转开发井全过程的环空压力精确预测，基于预测结果制定从探井设计到后期生产全寿命的环空压力管理措施，不仅为保障目标气田安全高效开发提供参考，也对于降低深水油气开发成本，推动深水油气资源勘探开发进程具有积极作用。

1 井筒温度压力分布模型

深水探井转开发井，除了油藏布井方面的考量，对于钻完井工程而言，需要在钻井前期研究阶段即对包括环空压力管理在内的一系列工程因素采用后期可转开发井的设计。深水气井的泥线温度和井底温度相差较大，较强的温度效应将引起密闭环空内流体升温膨胀，产生环空带压现象，严重时会造成套管挤毁、变形或上顶井口等事故。在陆地油田和浅海油田的勘探开发实践中，可以通过打开套管头侧翼阀释放环空压力，但在深水气田开发中，水下井口和生产系统设计的特殊结构使得环空压力难以监测控制和调节释放。目前对于深水井筒环空压力已有一些研究，但多针对某单一时段，如生产阶段，且井筒类型为单纯探井或者开发井，针对深水探井转开发井全阶段的环空压力研究极少。

阿里巴巴内部人士告诉记者，集团对于大文娱的表现并不满意，这也是大文娱成立两年多以来三度换帅的重要原因。从俞永福到杨伟东，再到如今的樊路远。

1.1 井筒温度分布模型

井筒温度是影响环空压力的核心因素，众多学者对于海上油气井井筒温度分布已有较多研究，建立了一系列相关模型[6-10]，但大多针对深水油气井测试或生产阶段[11-12]，热量传递过程稳定单一，而深水探井转开发井各阶段存在不同的热量传递模式(如图1所示)。

图1 深水探井转开发井各阶段井筒热量传递模式

探井钻井阶段(如图1a所示)，地层热量沿轴向、径向2个方向，在地层中深段、地层浅段、海水段3个层面进行叠加传导：在地层中深段，热量沿地层—水泥环/套管复合层—A环空传导，随着钻井泥浆正循环路径沿A环空流向井筒上部，这一过程对井筒中深段环空进行降温；热量随泥浆行进至某一地层深度(传热转换节点)，泥浆温度高于地层温度，热量径向传递路径转换，由A环空泥浆通过水泥环/套管复合层散到地层中，这一过程对井筒浅段进行升温；与此同时，钻井泥浆继续往上流动，进入海水段后，环空泥浆与环境温差骤增，热量快速大量散失。探井弃井后长期闲置阶段(如图1b所示)，井筒温度逐渐稳定，最终趋于接近地层原始温度，在地层浅部井段，微量的地层热量沿井筒传递至水下井口，对稳定后的井筒温度场几乎没有影响，在模型构建中忽略。完井重入与生产阶段(如图1c、图1d所示)，热量随地层流体沿生产油管向上传递，沿途径向上通过水泥环/套管复合层向中浅部地层传递，封隔器以上环空不存在流体流动热量，轴向上热量传递忽略。

基于以上热量传递过程，主要研究并建立了泥线以下地层段井筒和海水段井筒热量传递模型。 1.1.1 泥线以下地层段井筒传热模型

根据经典传热学理论，相较于对流换热系数、水泥环导热系数以及地层导热系数，钢质套管的导热系数是非常小的，可以忽略。因此，在泥线以下传热模型中，认为热传导的热量损失为以下几部分的叠加：

Cpa为环空流体的比热容，J/(kg·K)；

在密封的环空内加入可压缩的流体，以此来吸收因流体热膨胀而产生的高压。可压缩液体的种类很多，主要是指充有不同气体的各种类型泥浆：在工程上最常用的混入气体是氮气；泥浆的类型主要有复合油基泥浆、水基泥浆、盐水泥浆、或者淡水泥浆等。

(2)水泥环导热系数；

(3)地层的导热传热，为无限大介质中的一维非稳态传热问题。

因此，地层段传热数学模型归纳为：

(1)

式中：TD为无量纲温度；Kg为套管层导热系数，W/(m2·℃)；Tf为最外层水泥环温度，℃；Te为套管内温度，℃；Uto1为地层段总传热系数，W/(m2·℃)。

段文瀚解释说，第一个“绿色”之所以叫“绿色矿山”，是因为云天化采用最先进、效率最高的开采技术，极大提高了对磷矿石的利用率，做到了对资源的充分利用和循环利用。除此之外，云南的磷矿大都是胶磷矿，杂质含量较多，这就需要采用先进的浮选技术对胶磷矿进行加工，使其更干净、有害因子更少。同时，矿产开采完之后，云天化还会对土壤进行回填，利用先进技术进行浮土植被的种植。浮土植被的种植听起来简单，实际上对技术有很高的要求，前期投入巨大。据段文瀚介绍，上个世纪80年代至今，云天化单是复垦植被，恢复绿色的投入就将近10亿元，目的就是希望从源头支撑整个绿色发展理念的运行。

(2)

式中：hc为环空流体对流换热系数，W/(m2·℃)；hr为环空流体辐射传热系数，W/(m2·℃)；rto为油管外半径，m；rh为水泥环外半径，m；rco为套管外半径，m；Kcem为水泥环导热系数，W/(m2·℃)。

1.1.2 海水段井筒传热模型

能量守恒方程：

(3)

自然热对流中环空流体沿管壁的质量流量：

m=Vzρπ(δ+ri)2-ri2

(4)

“原则”的捍卫者——理论家们决不允许你随随便便地“搔痒”。 既然你想要以“瘙痒”的事实(或绝对正确性)对抗理论、击破理论，理论的事实(或绝对权威性)就一定想办法治疗你的“瘙痒”。 但理论治疗“瘙痒”的方法非常笨拙，那就是监狱和镣铐。 罗扎诺夫和陀思妥耶夫斯基一样，反对用监狱和镣铐对付地下人及其“瘙痒”：

δ=3.93Pr-0.5(0.952+Pr)0.25Grx-0.25X

(5)

格拉晓夫数Gr：

当夜，百里香去了小龙潭，将琵琶仙的话告诉了铁头大哥。第二天，又依琵琶仙之计，特意做了两笼香肠，抄僻巷小径送到北大营，果然畅通无阻，顺利接上了庄翻译，私下将琵琶仙的话传给了庄槐。

(6)

普朗特数Pr：

(7)

式中，

mf为生产流体质量流量，lb/h；

Cpf为生产流体比热容，J/(kg·K)；

Tbh为井底温度，℉；

Ts为井口温度，℉；

ma为由热对流引起的环空流体沿油管壁的质量流量，lb/h；

(1)井筒流体与套管的对流换热系数；

当前，我国对旅游行业的监管在一些方面存在不足之处。一些“有心人”钻空子，利用法律或制度漏洞，采取超低价竞争诱骗消费者跟团并强制购物等违规现象时有发生[1]。这些问题都反映出当前旅游业迫切需要通过完善相关法律制度来解决监管体系的漏洞。一些由于监管漏洞导致的问题，已经使部分游客的权益受到了侵害。这些漏洞和弊端不能得到解决，就无法有效地维护游客的权益。

Tba为环空流体在底部的平均温度，℉；

由出风口距端头不同距离L下的综掘工作面整体风流场分布图可以看到，风流从出风口射出后经过贴附射流区到达巷道壁面发生冲击，折返运动形成回流区；回风侧人行处最大风速可达到 4.625 m/s，超出《煤矿安全规程》规范的最大风速4 m/s[17]，掘进司机处受涡流影响，最小风速范围为0～0.05 m/s，低于《煤矿安全规程》规定的最小风速 0.25 m/s。风流分布不合理导致司机位置处及回风侧人行处粉尘积聚如图3所示。

Tta为环空流体在顶部的平均温度，℉；

L为井深，ft；

将环空温度计算结果导入环空压力计算模型，环空A和环空B全寿命过程环空压力变化如图4所示。

Ti为油管壁温度，℉；

To为套管壁温度，℉；

(1)实体经济经营风险加大。互联网金融的健康发展离不开自身所处的实体经济环境。2007年，我国出现了最早的互联网金融平台。2008年～2010年，在全球金融危机的影响下，我国实体经济萧条，互联网金融平台发展艰难。2010年～2012年，我国实体经济开始复苏，但互联网金融停滞不前。2012年～2015年，我国实体经济逐渐进入到上升阶段，互联网金融企业也开始快速发展。2015年至今，我国实体经济进入波动阶段，在实体经济经营风险加大的影响下，不少互联网金融机构平台无法适应外部环境变化纷纷倒闭破产。

ri为油管外径，ft；

“昆南”阳平声“年”的唱调（《牡丹亭·寻梦》【懒画眉】“最撩人春色是今年”，711）。因阳平声字腔的音势是平后略升，故即为其字腔。此后的十六分休止符即为“年”字腔的结点，其后的音调即为过腔。

ro为套管内径，ft；

Vz为沿热板向上流动的速度，ft/s；

B为两板之间流体的热膨胀系数；

ΔT为两板之间的温度差，℉；

μ为流体的粘度，mPa·s；

δ为边界层的厚度，ft；

X为热传递开始至结束的一个特定的长度，ft。

新建模型专门考虑了深水海水井段的特殊性，根据钻井或完井重入不同阶段，海水井段隔水管与海水、环空流体间存在受迫热对流、自然对流以及热传导等多种热量传递模式，模型计算过程中，重点调整了热量散失程度大的海水段环空流体的切力、比热容、粘度、导热系数以及密度等影响因素对井筒传热性能和非凝胶对流层的影响。

8月以来，南方已进入秋季备肥阶段，北方市场也在逐步启动，但复合肥市场的行情并没有往年同季“热闹”。较高的产品及原材料价格与积极性较低的农民用肥需求，导致复合肥市场处于谨慎状态，产品价格和产需量处于上下游矛盾态势下的平稳状态。预计后期市场不会出现较大变动，还需要对上下游多方因素保持关注。

1.2 环空压力计算

目前密闭环空压力的计算方法有刚性空间法和基于PVT状态方程的相关方法[13-15]，但以上方法难以计算多环空情况下的压力。张波、管志川等[16]认为井筒受热温度升高以后，环空和环空内的液体体积同时发生改变，液体与套管之间的热物性差异导致环空的有限体积难以容纳受热膨胀以后的液体。为满足体积相容性，环空压力上升对液体体积产生压缩效应，从而使液体实际体积与环空体积相等，整个环空体积变化值可以表述为：

ΔVa=πΔrao2+2rao-Δrai2-2raiΔrai)dz

(8)

式中，ΔVa为环空改变的体积，z1为环空起始位置，m；z2为井口位置，m；rao为环空外径，m；Δrai为环空内径变化值，m；rai为环空内径，m；Δrai为环空内径变化值，m。

2 全过程环空温度压力预测

L-X气田是南海琼东南盆地一个深水气田，为降低开发成本，将该气田若干口前期探井转为开发井，L-X-1井即为其中一口，基本参数见表1。

表1 L-X-1井基本参数

水深/m1347井深/m341236＂导管下深/m9220＂套管下深/m46313-3/8＂套管下深/m25249-5/8＂套管下深/m3410地层温度/℃857泥线温度/℃21地层压力/MPa39

2.1 全过程环空温度预测

探井钻井是井筒全寿命过程的起点，钻井过程，包括各层次套管的构建过程中，各井段钻井、固井、循环等环节热量传递模式变化大，热量传递影响因素繁多，作用机理复杂，而钻井过程中环空温度的改变都是短时间的瞬态变化，对后续过程井筒温度状态无过多影响，因此对钻井过程简化处理，仅计算最下部12-1/4＂井段钻井过程温度变化，将最内层9-5/8＂套管固井完成(不考虑弃井)后的井筒瞬态温度作为起始节点温度，导入下一过程，钻井终点泥线温度计算结果见表2，环空A和环空B较为接近，环空C和环空D远小于前者，这是由于该井12-1/4＂井段长度较短，9-5/8＂套管与上层13-3/8＂套管管鞋距离较小，热量在环空A和环空B的传递路径涵盖大段相同深度，表现出较小的温度差异性；而20＂套管和36＂导管下深较浅，管鞋位置已在传热转换节点上部，热量在环空C和环空D传递路径逆转，且需要通过2(3)层套管和1(2)层水泥环，热量传递损耗增加，传递量下降，表现出较大的温度差异性。

表2 钻井结束泥线环空温度

环空A环空B环空C环空D环空泥线位置温度/℃7364288

完井重入过程基于相似原因进行简化处理，以上一过程即弃井后闲置阶段终点温度作为完井重入起点温度，计算完井管柱构建完成，开始清井放喷时的温度，衔接生产过程。

在生产过程中，考虑后期见水对井筒温度的影响，生产10年周期出水情况见表3。

表3 生产10年内出水情况

生产时间/年12345678910含水率，‰0000010020030050101502

以泥线位置作为评估节点，深水探井转开发井井筒全寿命过程环空温度变化如图2所示，

图2 深水探井转开发井全过程环空温度变化

探井井筒在钻井结束弃井后进入闲置阶段，泥线位置井筒环空温度迅速下降，环空A温度在10天内降幅达到88.15%，之后下降趋势放缓，在47天降至4.73℃，之后基本稳定；环空B温度变化与环空A较为接近，在10天内下降86.42%，之后下降趋势放缓，在33天降至4.62℃，之后基本稳定；环空C和环空D钻井终点温度远小于环空A和环空B，温度更快降至接近泥线温度的稳定状态，分别在第9天和第1.5天降至3.49℃和3.72℃。

选取我院2017年3月-2018年3月期间收治的急性阑尾炎患者60例为研究对象,均分为两组,每组30例,其中,对照组男17例,女13例,年龄在22-75岁,平均年龄为(37.21±1.25)岁;观察组男18例,女12例,年龄在23-76岁,平均年龄为(37.26±1.24)岁。两组在(年龄、性别)等方面,统计学无意义(P>0.05)。纳入标准:依据WHO中关于急性阑尾炎的诊断标准[2],确诊为急性阑尾炎患者;所有患者均签署知情同意书。排除不签署知情同意书者;排除精神疾病者。

图3 放喷与生产6个月后全井段环空温度

钻井阶段与完井重入阶段环空温度变化最大，但钻井阶段时间较短且随后立即转入长时间的弃井闲置阶段，温度迅速下行恢复，针对环空压力预测，重点考虑完井重入放喷与生产阶段温度的影响，放喷与生产6个月后全井段环空温度计算结果如图3所示，4个环空温度在全井段表现出的差异性与全过程表现出的温度趋势一致；在井深1 681 m位置温度出现一个台阶变化，这一深度是9-5/8＂套管固井水泥返高深度，水泥与上部留空是完全不同的热量传递介质，热量在介质界面上下具有不同的传递属性，最终在环空温度表现出这一台阶差异。

2.2 环空压力预测

环空D水泥返至泥线，不会发生环空带压；环空C下深浅且水泥返高大，在此情况下温度对环空压力的影响非常小，因此重点关注环空A和环空B的压力变化。

k为环空流体的导热系数，W/(m2·℃)；

图4 全过程环空压力变化

随着井筒闲置阶段后温度迅速降低，环空A和环空B压力也随着下降，但压力降落曲线形态与温度降落曲线形态存在明显区别：温度的下降非常急剧，在弃井后的10天内就降至一个相对稳定的低温状态，而环空压力由初始33.55 MPa下降至闲置阶段终点的10.32 MPa，压力降落过程缓慢且均匀的分布在整个闲置阶段(90天)，环空压力与温度呈非线性相关性，表现出温度对压力影响的延时性。

这种延时性同样体现在完井重入后的环空压力变化过程中：以环空A为例，在完井重入的20天时间内，环空A温度就由闲置阶段终点的4.73℃急剧增至64.12℃，之后趋于相对稳定，而环空A压力在完井重入结束进入生产阶段时，仅由闲置阶段终点的15.31 MPa增至18.10 MPa，随着生产进行，在189天才逐渐增至30.02 MPa，之后趋于缓慢稳定增长。在生产至第4年开始见水后，环空压力没有出现明显的变化。

3 环空压力干预

固井水泥设计是影响环空压力的关键因素，但影响机理复杂，包括水泥浆性能、水泥返高设计、水泥环强度等都会在生产过程中与套管、留空气体等介质产生复杂的动态耦合变化，需要构建专门模型研究，这里不作赘述。

除固井水泥设计外，环空压力干预途径主要有2种：破裂盘和可压缩材料。

3.1 破裂盘

破裂盘是一个用来限制井筒或套管环空压力的一次性压力释放保护装置。通常把向外破裂时设计的装置称为破裂盘；而把向内破裂的装置称为坍塌盘。根据额定压力的大小划分为单向破裂盘和双向破裂盘。在套管柱中安装破裂盘(Burst Disk)：在外层套管上安装一到二个破裂盘，当密封环空内压力达到破裂盘的破裂压力时，外层套管上的破裂盘破裂，从而保护内、外套管不被挤毁或压裂，同时保证了内层套管串的完整性。

使用破裂盘，通常需要在外层套管留出地层沟通通道作为破裂盘激活后的泻压通道，会导致套管管柱存在薄弱点且破裂后内外层空间连通，存在安全隐患。

3.2可压缩材料

可压缩材料通过多种形式应用于环空压力干预中，其主要原理是通过到达临界压力后材料发生的压缩变形，为环空提供一定的流体膨胀空间和压力释放空间，进而降低环空压力。

(1)可压缩泡沫材料

针对库区消落带所引发的一系列环境生态问题，大量专家和学者对三峡库区消落带的环境治理以及生态修复进行研究和实践尝试。植被作为生态功能的主体，在消落带生态修复中至关重要，目前消落带生态修复的问题主要是在两个方面:一是适合消落环境植被筛选; 二是提供适应植被生长的环境，对三峡库区消落带生态修复的研究也大多集中在适合水淹植物物种的筛选以及适宜植物生长的生境构筑两方面。

在套管外安装可压缩泡沫，是目前国外深水井中常用的一种减小由温度引起的套管附加载荷的方法。原理是在内层套管上安装一定数量的可压缩的泡沫材料，当环空压力增加到一定程度时，可压缩泡沫材料开始变形，产生一定流体膨胀的空间，降低环空压力。

(2)可破裂泡沫球

(1)深水探井转开发井经历探井钻井段、弃井闲置、完井重入和生产4个阶段，各阶段热量传递模式有所区别，建立针对性温度计算模型，精确预测探井转开发井全过程温度变化，为环空压力预测提供基础参数。

(3)可压缩液体

“终极关怀”作为一个明确的问题域是由蒂利希提出的，但是事实上，它作为哲学的核心理论旨趣却是与生俱来的。马克思哲学终极关怀维度作为马克思哲学的应有之义，集中体现了马克思哲学对人的意义世界的积极探寻和构建，由于马克思对实践和自由的本真性阐释和践行，从而使人安身立命问题的解决获得了比以往思想理论更大的现实性和科学性。因此，对于马克思哲学终极关怀问题的研究就成为一项具有重要理论意义和现实意义的工作。

3.3 目标井环空压力干预

无论使用哪种可压缩材料，最终目的都是在环空压力达到临界值时激活材料压缩变形，为环空提供一个冗余空间，释放环空压力。根据设定的临界安全压力和前文所述模型，反算环空膨胀体积，进而获得所需冗余体积，即最小可压缩体积。如设定25 MPa和20 MPa 为L-X-1井环空A和环空B临界压力，反算体积结果见表4。

传热边界厚度δ：

表4 环空所需冗余及设计可压缩体积 /m3

环空所需冗余体积设计可压缩体积环空A03205环空B02103

环空压力在生产阶段变化如图5所示。

图5 环空压力干预

环空压力达到设定的临界压力时，压缩材料激活，开始提供冗余体积泻压，环空A泻压至21.87 MPa，环空B泻压至17.83 MPa，材料不再压缩，之后趋于相对稳定。

4 结论

该预防措施是在套管环空内放入一定数量的合成可破裂的泡沫球，这些合成的泡沫球内部充满空气，当密封环空内压力达到某一数值时，泡沫球就会破裂，小球体积减小，释放一定的空间，从而降低环空压力。可破裂泡沫球由复合材料制成，最常见的是空心玻璃球，其内充满标准大气压的空气。小球本身有一定的强度，而且有很小的加工公差，因此能保证在某一确定的压力下小球破裂。小球的直径根据环空空间的大小而不同，一般在19.05～38.10 mm。同时由于环空液体的压缩性不大，因此加入的可破裂泡沫球的体积长度一般占整个环空体积长度的2%～8%，大约相当于2～20个套管单根的长度。

(2)计算了环空A和环空B的从弃井闲置到生产10年后的全过程环空压力，结果表明环空压力与温度呈非线性相关性，表现出温度对压力影响的延时性。

(3)使用压缩材料进行环空压力干预，设定25 MPa和20 MPa 为目标井环空A和环空B临界压力，反算环空膨胀体积，进而获得所需冗余体积，根据所需冗余体积设计可压缩体积，在环空压力达到临界值时激活材料压缩变形，为环空提供一个冗余空间，有效降低环空压力。

综上所述，语文这一门学科在与语言的应用上具有综合性的同时还有实践性。因此，在语文教学活动中需要做到的不再是以往空洞式的说教，现在开始重视对词语的品读、对句子的赏析，更注重课外拓展，使语文教学回归本真。本文对部编版三年级语文“语用型教学”措施展开分析，指出“语用型教学”的教学策略，就是需要贴合学生，找到语用的起点、品读语言，对语用的基础加以夯实，对写作手法加以琢磨，学会语用的方法。与学生的实习生活相结合，站在学生的视角下对语文教学内容进行解读。在对课文进行学习的重点就是对语言加以品读。对词语意思的理解，就是对语言的品读，对课文加以学习，指引着学生对写作的手法加以琢磨。

(4)深水探井转开发井是降低深水油气开发成本的有效途径，这一过程中的环空压力变化和干预是影响深水油气田安全生产的重要因素。解决深水探井转开发井全过程环空压力预测与干预问题，对于今后我国深水油气资源经济高效开发具有积极的推动作用。

参考文献:

[1] 吴时国，袁圣强.世界深水油气勘探进展与我国南海深水油气前景[J].天然气地球科学，2005，16(6):63-69.

[2] 江怀友，赵文智，裘怿楠，等.世界海洋油气资源现状和勘探特点及方法[J].中国石油勘探，2008，5(3):27-29.

[3] 李中，方满宗，李磊.南海西部深水钻井实践[J].石油钻采工艺，2015，37(1):92-95.

[4] 胡光义，孙福街，范廷恩，等. 海上油气田勘探开发一体化理念、基本思路和对策[J].中国海上油气，2013，25(6)：61-65.

[5] 赵会卿，张新建，孙建津，等.海上边际油气田勘探开发一体化设想[J].中国海上油气，2008，20(6)：421-425.

[6] 高永海，孙宝江，王志远，等.深水钻探井筒温度场的计算与分析[J].中国石油大学学报(自然科学版)，2008，32(2):58-62.

[7] 钟兵，方铎，施太和.井内温度影响因素的敏感性分析[J].天然气工业.2000，20(2):57-60.

[8] 何世明，何平，尹成，等.井下循环温度模型及其敏感性分析[J].西南石油学院学报，2002，24(1):57-60.

[9] HARRIS O O，OSISANYA O S .evaluation of equivalent circulating density of drilling fluids under high-pressure/high-temperature conditions[J].SPE 97018 ，2005(10):9-12.

[10] 宋洵成，管志川.深水钻井井筒全瞬态传热特征[J].石油学报.2011，32(4):704-708.

[11] 张崇，任冠龙，董钊，等.深水测试井筒温度场预测模型的建立及应用[J].中国海上油气，2016，28(5)：78-84.

[12] 董钊，余意，任冠龙，等. 南海超深水浅层低温气井测试温度场干预[J].石油钻采工艺，2016，38(5)：661-666.

[13] NOJABAEI B，HASAN AR，KABIR，C S．Modelling wellbore transient flu-id temperature and pressure during diagnostic fracture-injection tes-ting in unconventional reservoirs[J]．Journal of Canadian Petroleum Technology，2014，53(3):161-167．

[14] 杨进，唐海雄，刘正礼，等．深水油气井套管环空压力预测模型[J]．石油勘探与开发，2013，40(5):616-619．

[15] 车争安，张智，施太和，等．高温高压含硫气井环空流体热膨胀带压机理[J]．天然气工业，2010，30(2):88-90．

[16] 张波，管志川，张琦.深水油气井开采过程环空压力预测与分析[J]. 石油学报，2015，36(8)：1012-1017.