《一种气井标况下体积泄漏速率获取方法及装置.pdf》由会员分享，可在线阅读，更多相关《一种气井标况下体积泄漏速率获取方法及装置.pdf（17页珍藏版）》请在专利查询网上搜索。
1、10申请公布号CN104131811A43申请公布日20141105CN104131811A21申请号201410374225322申请日20140731E21B47/1020120171申请人中国石油天然气股份有限公司地址100007北京市东城区东直门北大街9号72发明人王云李文魁李隽曹光强74专利代理机构北京三友知识产权代理有限公司11127代理人汤在彦54发明名称一种气井标况下体积泄漏速率获取方法及装置57摘要本发明涉及一种气井标况下体积泄漏速率获取方法及装置，其中，该方法包括根据井下泄漏点处环空压力与油管内压力平衡，获得泄漏点深度、油管内泄漏点处气体压力以及对应温度；根据油管内泄漏点处。
2、气体压力以及不同时刻泄漏点处环空压力判断泄漏的流动状态；对所述泄漏点处质量泄漏速率表达式转换得到泄漏点处体积泄漏速率表达式，对泄漏点处体积泄漏速率表达式进行关于时间的积分，得到对应时间内泄漏点处的体积泄漏量表达式；根据泄漏点处体积泄漏量确定任意泄压时间段内泄漏点处的泄漏面积；根据任意泄压时间段内泄漏点处的泄漏面积获得气井油套环空泄漏气体在标况下的体积泄漏速率。51INTCL权利要求书3页说明书10页附图3页19中华人民共和国国家知识产权局12发明专利申请权利要求书3页说明书10页附图3页10申请公布号CN104131811ACN104131811A1/3页21一种气井标况下体积泄漏速率获取方法。
3、，其特征在于，该方法包括根据井下泄漏点处环空压力与油管内压力平衡，获得泄漏点深度、油管内泄漏点处气体压力以及对应温度；根据油管内泄漏点处气体压力以及不同时刻泄漏点处环空压力判断泄漏的流动状态；若泄漏的流动状态为亚临界流动状态，则确定泄漏点处质量泄漏速率表达式为若泄漏的流动状态为临界流动状态，则确定泄漏点处质量泄漏速率表达式为其中，QM表示泄漏点处质量泄漏速率，KG/S；P1表示油管内泄漏点气体压力，PA；PO表示泄漏点处环空压力，PA；TA表示泄漏点处温度，K；A表示泄漏面积，M2；表示流速系数，表示实际流速与理论流速之比，无因次；K表示绝热指数，即定压热容与定容热容之比，无因次；CD表示流量。
4、系数，与泄漏孔的形状有关，无因次；M表示天然气摩尔质量，KG/MOL；R表示摩尔气体常数，8314JMOLK；Z表示气体体积偏差因子，无因次；对所述泄漏点处质量泄漏速率表达式转换得到泄漏点处体积泄漏速率表达式，对泄漏点处体积泄漏速率表达式进行关于时间的积分，得到对应时间内泄漏点处的体积泄漏量表达式；根据泄漏点处体积泄漏量确定任意泄压时间段内泄漏点处的泄漏面积；根据任意泄压时间段内泄漏点处的泄漏面积获得气井油套环空泄漏气体在标况下的体积泄漏速率。2如权利要求1所述的方法，其特征在于，所述油管内泄漏点处气体压力以及不同时刻泄漏点处环空压力判断泄漏的流动状态的步骤包括当时，泄漏的流动状态为亚临界流动。
5、状态；当时，泄漏的流动状态为临界流动状态；其中，3如权利要求1所述的方法，其特征在于，对所述泄漏点处质量泄漏速率表达式转换得到泄漏点处的体积泄漏速率表达式的步骤具体为所述泄漏点处质量泄漏量表达式按照式转换得到泄漏点处的体积泄漏速率表达式；其中，QMI表示为I时刻泄漏点处的质量泄漏速率，KG/S；QVI表示为I时刻泄漏点处的体积泄漏速率，M3/S；GI表示为I时刻泄漏点处的天然气密度，KG/M3。4如权利要求3所述的方法，其特征在于，所述对所述泄漏点处体积泄漏速率表达式进行关于时间的积分的步骤具体为对所述泄漏点处体积泄漏速率表达式按照式进行关于时间的积分，得到对应时间内泄漏点处的体积泄漏量表达式。
6、；其中，Q表示为任意时间段内泄漏点处累计权利要求书CN104131811A2/3页3体积泄漏量，M3；表示等效积分量，M。5如权利要求4所述的方法，其特征在于，所述根据泄漏点处的体积泄漏速率确定任意泄压时间段内泄漏点处的泄漏面积的步骤具体包括当泄压平衡时，气体在井下泄漏与井口泄漏达到平衡，即VSEGVEGA；其中，V表示任意时间段T内井下体积泄漏量，M3，VS表示任意时间段T内井口体积泄漏量，M3，天然气井下和井口的体积相差一个天然气体积膨胀系数EG，M3/M3；Z表示天然气体积偏差因子，无因次；若现场测得任意时间段T内井口体积泄漏量VS，则直接根据VSEGVEGA确定任意泄压时间段内泄漏点处。
7、的泄漏面积；否则，根据所述泄漏点处体积泄漏量表达式获取任意时间段内井口累计体积泄漏量，根据VSEGVEGA确定任意泄压时间段T内泄漏点处对应的泄漏面积。6如权利要求5所述的方法，其特征在于，所述根据任意泄压时间段内泄漏点处的泄漏面积获得气井油套环空泄漏气体在标况下的体积泄漏速率的具体步骤包括根据所述泄漏点处的体积泄漏速率表达式以及任意泄压时间段内泄漏点处的泄漏面积获得标况下的体积泄漏速率；其中，若泄漏的流动状态为亚临界流动状态，则标况下的体积泄漏速率为若泄漏的流动状态为临界流动状态，则标况下的体积泄漏速率为7一种气井标况下体积泄漏速率获取装置，其特征在于，该装置包括第一处理单元，用于根据井下泄。
8、漏点处环空压力与油管内压力平衡，获得泄漏点深度、油管内泄漏点处气体压力以及对应温度；判断单元，用于根据油管内泄漏点处气体压力以及不同时刻泄漏点处环空压力判断泄漏的流动状态；若泄漏的流动状态为亚临界流动状态，则确定泄漏点处质量泄漏速率表达式为若泄漏的流动状态为临界流动状态，则确定泄漏点处质量泄漏速率表达式为其中，QM表示泄漏点处质量泄漏速率，KG/S；P1表示油管内泄漏点气体压力，PA；PO表示泄漏点处环空压力，PA；TA表示泄漏点处温度，K；A表示泄漏面积，M2；表示流速系数，表示实际流速与理论流速之比，无因次；K表示绝热指数，即定压热容与定容热容之比，无因次；CD表示流量系数，与泄漏孔的形状。
9、有关，无因次；M表示天然气摩尔质量，KG/MOL；R表示摩尔气体常数，8314JMOLK；Z表示气体体积偏差因子，无因次；第二处理单元，包括转换模块、积分模块和确定模块；所述转换模块用于对所述泄漏点处质量泄漏速率表达式转换得到泄漏点处的体积泄漏速率表达式；所述积分模块用于对泄漏点处体积泄漏速率表达式进行关于时间的积分，得到对应时间内泄漏点处的体积泄漏量权利要求书CN104131811A3/3页4表达式；所述确定模块用于根据泄漏点处体积泄漏速率确定任意泄压时间段内泄漏点处的泄漏面积；标况下体积泄漏速率获取单元，用于根据任意泄压时间段内泄漏点处的泄漏面积获得气井油套环空泄漏气体在标况下的体积泄漏速。
10、率。8如权利要求7所述的装置，其特征在于，所述判断单元具体用于当时，泄漏的流动状态为亚临界流动状态；当时，泄漏的流动状态为临界流动状态；其中，9如权利要求7所述的装置，其特征在于，所述转换模块具体用于所述泄漏点处质量泄漏速率表达式按照式转换得到泄漏点处的体积泄漏速率表达式；其中，QMI表示为I时刻泄漏点处的质量泄漏速率，KG/S；QVI表示为I时刻泄漏点处的体积泄漏速率，M3/S；GI表示为I时刻泄漏点处的天然气密度，KG/M3。10如权利要求9所述的装置，其特征在于，所述积分模块具体用于对所述泄漏点处体积泄漏速率表达式按照式进行关于时间的积分，得到对应时间内泄漏点处的体积泄漏量表达式；其中，。
11、Q表示为任意时间段内泄漏点处累计体积泄漏量，M3；表示等效积分量，M。权利要求书CN104131811A1/10页5一种气井标况下体积泄漏速率获取方法及装置技术领域0001本发明涉及采气工程技术领域，特别涉及一种气井标况下体积泄漏速率获取方法及装置。背景技术0002目前，随着天然气勘探开发的不断发展，一大批高温、高压、含酸性介质气井试产或投产，随之而来的气井完整性问题也愈发突出。气井的完整性要求是在整个寿命期间，通过管理和技术手段确保地下流体不发生无控制的溢流。如果发生完整性失效问题，轻则导致井下天然气进入油套环空，即环空带压，严重者可导致环空异常带压，甚至导致井下天然气泄漏至地面、发生爆炸，。
12、给人民生命财产安全造成巨大损失，并引发严重的社会负面影响。0003气井完整性失效最根本的原因是井屏障组件失效导致天然气发生泄漏，彻底解决该问题就是进行修井作业，如果不加以区分，对所有环空带压井全部进行修井作业，一方面作业成本太高，以塔里木油田为例，一口气井的修井成本在2000万元以上；另一方面高温、高压、含酸性介质气井的修井作业风险很大。因此，非常有必要评价气井的完整性，按照不同完整性等级区别对待。0004目前，气井完整性评价的主要标准挪威NORSOKD010油气井钻井与作业时的完整性要求和挪威石油联合会OLF117标准油气井完整性推荐做法指南都提出了以井屏障组件的标况下体积泄漏速率作为评价气。
13、井完整性最基本的参数气体042M3/MIN，液体040L/MIN。0005国内外对于油气井标况下体积泄漏速率计算研究较少。在国内外公开发表的文献资料中，还未有采用数学理论方法可同时确定标况下体积泄漏速率泄漏面积的报道。目前已知挪威SCANWELL公司研发一套设备，可在井口测量井下泄漏速率及泄漏点深度，该技术在国外已经现场应用，在国内还未见现场实验。发明内容0006为解决现有技术无法确定气井标况下体积泄漏速率的问题，本发明提出一种气井标况下体积泄漏速率获取方法及装置，利用气井常规基础数据，以地面输气管线小孔泄漏模型为基础，采用积分原理，确定气井井下环空泄漏面积，实现标况下体积泄漏速率和泄漏面积。。
14、0007为实现上述目的，本发明提供了一种气井标况下体积泄漏速率获取方法，该方法包括0008根据井下泄漏点处环空压力与油管内压力平衡，获得泄漏点深度、油管内泄漏点处气体压力以及对应温度；0009根据油管内泄漏点处气体压力以及不同时刻泄漏点处环空压力判断泄漏的流动状态；若泄漏的流动状态为亚临界流动状态，则确定泄漏点处质量泄漏速率表达式为说明书CN104131811A2/10页6若泄漏的流动状态为临界流动状态，则确定泄漏点处质量泄漏速率表达式为其中，QM表示泄漏点处质量泄漏速率，KG/S；P1表示油管内泄漏点气体压力，PA；PO表示泄漏点处环空压力，PA；TA表示泄漏点处温度，K；A表示泄漏面积，M。
15、2；表示流速系数，表示实际流速与理论流速之比，无因次；K表示绝热指数，即定压热容与定容热容之比，无因次；CD表示流量系数，与泄漏孔的形状有关，无因次；M表示天然气摩尔质量，KG/MOL；R表示摩尔气体常数，8314JMOLK；Z表示气体体积偏差因子，无因次；0010对所述泄漏点处质量泄漏速率表达式转换得到泄漏点处体积泄漏速率表达式，对泄漏点处体积泄漏速率表达式进行关于时间的积分，得到对应时间内泄漏点处的体积泄漏量表达式；根据泄漏点处体积泄漏量确定任意泄压时间段内泄漏点处的泄漏面积；0011根据任意泄压时间段内泄漏点处的泄漏面积获得气井油套环空泄漏气体在标况下的体积泄漏速率。0012优选地，所述。
16、油管内泄漏点处气体压力以及不同时刻泄漏点处环空压力判断泄漏的流动状态的步骤包括0013当时，泄漏的流动状态为亚临界流动状态；当时，泄漏的流动状态为临界流动状态；其中，0014优选地，述对所述泄漏点处质量泄漏速率表达式转换得到泄漏点处的体积泄漏速率的步骤具体为0015所述泄漏点处质量泄漏量表达式按照式转换得到泄漏点处的体积泄漏速率表达式；其中，QMI表示为I时刻泄漏点处的质量泄漏速率，KG/S；QVI表示为I时刻泄漏点处的体积泄漏速率，M3/S；GI表示为I时刻泄漏点处的天然气密度，KG/M3。0016优选地，所述对所述泄漏点处体积泄漏速率表达式进行关于时间的积分的步骤具体为0017对所述泄漏点。
17、处体积泄漏速率表达式按照式进行关于时间的积分，得到对应时间内泄漏点处的体积泄漏量表达式；其中，Q表示为任意时间段内泄漏点处累计体积泄漏量，M3；表示等效积分量，M。0018优选地，所述根据泄漏点处的体积泄漏速率确定任意泄压时间段内泄漏点处的泄漏面积的步骤具体包括0019当泄压平衡时，井下泄漏与井口泄漏达到平衡，即VSEGVEGA；其中，V表示任意时间段T内井下体积泄漏量，M3，VS表示任意时间段T内井口体积泄漏量，M3，天然气井说明书CN104131811A3/10页7下和井口的体积相差一个天然气体积膨胀系数EG，M3/M3；Z表示天然气体积偏差因子，无因次；0020若现场测得任意时间段T内井。
18、口体积泄漏量VS，则直接根据VSEGVEGA确定任意泄压时间段内泄漏点处的泄漏面积；否则，根据所述泄漏点处体积泄漏量表达式获取任意时间段内井口累计体积泄漏量，根据VSEGVEGA确定任意泄压时间段T内泄漏点处对应的泄漏面积。0021优选地，所述根据任意泄压时间段内泄漏点处的泄漏面积获得气井油套环空泄漏气体在标况下的体积泄漏速率的具体步骤包括0022根据所述泄漏点处的体积泄漏速率表达式以及任意泄压时间段内泄漏点处的泄漏面积获得标况下的体积泄漏速率；其中，若泄漏的流动状态为亚临界流动状态，则标况下的体积泄漏速率为若泄漏的流动状态为临界流动状态，则标况下的体积泄漏速率为0023为实现上述目的，本发明。
19、还提供了一种气井标况下体积泄漏速率获取装置，该装置包括0024第一处理单元，用于根据井下泄漏点处环空压力与油管内压力平衡，获得泄漏点深度、油管内泄漏点处气体压力以及对应温度；0025判断单元，用于根据油管内泄漏点处气体压力以及不同时刻泄漏点处环空压力判断泄漏的流动状态；若泄漏的流动状态为亚临界流动状态，则确定泄漏点处质量泄漏速率表达式为若泄漏的流动状态为临界流动状态，则确定泄漏点处质量泄漏速率表达式为其中，QM表示泄漏点处质量泄漏速率，KG/S；P1表示油管内泄漏点气体压力，PA；PO表示泄漏点处环空压力，PA；TA表示泄漏点处温度，K；A表示泄漏面积，M2；表示流速系数，表示实际流速与理论流。
20、速之比，无因次；K表示绝热指数，即定压热容与定容热容之比，无因次；CD表示流量系数，与泄漏孔的形状有关，无因次；M表示天然气摩尔质量，KG/MOL；R表示摩尔气体常数，8314JMOLK；Z表示气体体积偏差因子，无因次；0026第二处理单元，包括转换模块、积分模块和确定模块；所述转换模块用于对所述泄漏点处质量泄漏速率表达式转换得到泄漏点处的体积泄漏速率表达式；所述积分模块用于对泄漏点处体积泄漏速率表达式进行关于时间的积分，得到对应时间内泄漏点处的体积泄漏量表达式；所述确定模块用于根据泄漏点处体积泄漏速率确定任意泄压时间段内泄漏点处的泄漏面积；0027标况下体积泄漏速率获取单元，用于根据任意泄压。
21、时间段内泄漏点处的泄漏面积获得气井油套环空泄漏气体在标况下的体积泄漏速率。说明书CN104131811A4/10页80028优选地，所述判断单元具体用于当时，泄漏的流动状态为亚临界流动状态；当时，泄漏的流动状态为临界流动状态；其中，0029优选地，所述转换模块具体用于所述泄漏点处质量泄漏速率表达式按照式转换得到泄漏点处的体积泄漏速率；其中，QMI表示为I时刻泄漏点处的质量泄漏速率，KG/S；QVI表示为I时刻泄漏点处的体积泄漏速率，M3/S；GI表示为I时刻泄漏点处的天然气密度，KG/M3。0030优选地，所述积分模块具体用于对所述泄漏点处体积泄漏速率表达式按照式进行关于时间的积分，得到对应时。
22、间内泄漏点处的体积泄漏量表达式；其中，Q表示为任意时间段内泄漏点处累计体积泄漏量，M3；表示等效积分量，M。0031上述技术方案具有如下有益效果本技术方案应用气井日常生产数据和油套环空泄压参数，通过求任意时间段内泄漏点处体积泄漏速率关于时间的积分，确定泄漏点处的泄漏面积，最终确定气井井油套环空泄漏气体在标况下体积泄漏速率模型，首次通过数学理论计算方法确定气井泄漏点处的泄漏面积，利用获得的气井泄漏点处的泄漏面积建立标况下体积泄漏速率计算表达式，减少了现有利用仪器设备现场测量带来的成本过高，以及测量过程风险大的弊端。附图说明0032为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例。
23、或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。0033图1为本发明提出的一种气井标况下体积泄漏速率获取方法流程图；0034图2为本发明提出的一种气井标况下体积泄漏速率获取装置框图；0035图3为本实施例的气井泄漏模型示意图；0036图4为本实施例的气井环空泄压曲线。具体实施方式0037下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本。
24、领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。0038根据文献资料可知，关于天然气地面输气管线的泄漏速率模型研究较多，也较为成熟，主要有小孔模型、大孔模型和管道模型。小孔模型假设气体泄漏是绝热过程，天然气说明书CN104131811A5/10页9为理想气体考虑，泄漏过程，管道内压力维持不变，根据伯努利方程和绝热方程，可以得出气体泄漏速度计算公式。计算天然气地面输气管线小孔泄漏，输气管线内外压力、泄漏面积都是已知的。0039借鉴小孔泄漏速率模型确定气井标况下体积泄漏速率模型，实现气井标况下体积泄漏速率和泄漏面积的数学计算，不需要额外的测试设备和测试作业。
25、量，依靠气井自身的参数和油套环空泄压数据就可以确定气井标况下体积泄漏速率和泄漏面积。0040经上述分析可知，本发明技术方案的工作原理为利用高温、高压气井自身基础参数和油套环空泄压参数，确定井下泄漏点深度；通过求任意时间段内泄漏点处体积泄漏速率关于时间的积分，确定泄漏点处泄漏面积，最终得出气井标况下体积泄漏速率方程，包括以下步骤1确定泄漏点深度；2确定泄漏点处体积泄漏速率表达式；3确定任意泄压时间段内泄漏点处的泄漏面积；4确定标况下体积泄漏速率表达式。这样能够实现气井标况下体积泄漏速率和泄漏面积的计算，弥补目前气井标况下体积泄漏速率和泄漏面积无法计算的不足，为气井完整性评价提供必需的技术手段与参。
26、数，避免现场测量与安全隐患。0041如图1所示，为本发明提出的一种气井标况下体积泄漏速率获取方法流程图。该方法包括0042步骤101根据井下泄漏点处环空压力与油管内压力平衡，获得泄漏点深度、油管内泄漏点处气体压力以及对应温度；0043步骤102根据油管内泄漏点处气体压力以及不同时刻泄漏点处环空压力判断泄漏的流动状态；若泄漏的流动状态为亚临界流动状态，则确定泄漏点处质量泄漏速率表达式为若泄漏的流动状态为临界流动状态，则确定泄漏点处质量泄漏速率表达式为其中，QM表示泄漏点处质量泄漏速率，KG/S；P1表示油管内泄漏点处气体压力，PA；PO表示泄漏点处环空压力，PA；TA表示泄漏点处温度，K；A表示。
27、泄漏面积，M2；表示流速系数，表示实际流速与理论流速之比，无因次；K表示绝热指数，即定压热容与定容热容之比，无因次；CD表示流量系数，与泄漏孔的形状有关，无因次；M表示天然气摩尔质量，KG/MOL；R表示摩尔气体常数，8314JMOLK；Z表示气体体积偏差因子，无因次；0044步骤103对所述泄漏点处质量泄漏速率表达式转换得到泄漏点处体积泄漏速率表达式，对泄漏点处体积泄漏速率表达式进行关于时间的积分，得到对应时间内泄漏点处的体积泄漏量表达式；根据泄漏点处体积泄漏量确定任意泄压时间段内泄漏点处的泄漏面积；0045步骤104根据任意泄压时间段内泄漏点处的泄漏面积获得气井油套环空泄漏气体在标况下的体。
28、积泄漏速率。0046如图2所示，为本发明提供的一种气井标况下体积泄漏速率获取装置框图。该装置包括0047第一处理单元201，用于根据井下泄漏点处环空压力与油管内压力平衡，获得泄漏说明书CN104131811A6/10页10点深度、油管内泄漏点处气体压力以及对应温度；0048判断单元202，用于根据油管内泄漏点处气体压力以及不同时刻泄漏点处环空压力判断泄漏的流动状态；若泄漏的流动状态为亚临界流动状态，则确定泄漏点处质量泄漏速率表达式为若泄漏的流动状态为临界流动状态，则确定泄漏点处质量泄漏速率表达式为其中，QM表示泄漏点处质量泄漏速率，KG/S；P1表示油管内泄漏点气体压力，PA；PO表示泄漏点处。
29、环空压力，PA；TA表示泄漏点处温度，K；A表示泄漏面积，M2；表示流速系数，表示实际流速与理论流速之比，无因次；K表示绝热指数，即定压热容与定容热容之比，无因次；CD表示流量系数，与泄漏孔的形状有关，无因次；M表示天然气摩尔质量，KG/MOL；R表示摩尔气体常数，8314JMOLK；Z表示气体体积偏差因子，无因次；0049第二处理单元203，包括转换模块2031、积分模块2032和确定模块2033；所述转换模块2031用于对所述泄漏点处质量泄漏速率表达式转换得到泄漏点处体积泄漏速率表达式；所述积分模块2032用于对泄漏点处体积泄漏速率表达式进行关于时间的积分，得到对应时间内泄漏点处的体积泄漏。
30、量表达式；所述确定模块2033用于根据泄漏点处体积泄漏量确定任意泄压时间段内泄漏点处的泄漏面积；0050标况下体积泄漏速率获取单元204，用于根据任意泄压时间段内泄漏点处的泄漏面积获得气井油套环空泄漏气体在标况下的体积泄漏速率。00511确定泄漏点深度0052根据气井基本数据以及已有的压力数据，通过拟合对比，优选适合本气井的多相管流计算模型，计算得出气井的温度、压力剖面，并回归温度、压力剖面的计算公式。根据平衡原理，在井口环空压力稳定平衡时，认为在井下泄漏点处环空压力与油管内压力平衡，即式1。0053P1POPAAGH10054其中，P1为油管内泄漏点气体压力，PA；PO为泄漏点处环空压力，P。
31、A；PA为井口处环空压力，PA；G98M/S2；A为环空保护液的密度，KG/M3；H为泄漏点深度，M；0055本计算方法假设00561油套环空仅存在一处泄漏点；00572油套环空中气柱的高度HL可以忽略；现场气井环空充满保护液，即使未完全充满，HL也很低，由HL引起的压差也很小，可以忽略不计。0058根据平衡环空压力值，可得泄漏点深度。00592确定泄漏点处体积泄漏速率表达式0060可根据地面输气管线小孔泄漏模型，即式2和式3，得出泄漏点处质量泄漏速率的计算公式0061说明书CN104131811A107/10页1100620063其中，QM为泄漏点处质量泄漏速率，KG/S；P1为油管内泄漏点。
32、气体压力，PA；PO为泄漏点处环空压力，PA；TA为泄漏点处温度，K；A为泄漏面积，M2；为流速系数，表示实际流速与理论流速之比，一般为097098，无因次；K为绝热指数，即定压热容与定容热容之比，对于天然气一般取K128，无因次；CD为流量系数，取值范围在0610，与泄漏孔的形状有关，圆孔取10，无因次；M为天然气摩尔质量，KG/MOL；R为摩尔气体常数，8314JMOLK；Z为气体体积偏差因子，无因次。0064式2和式3分别是输气管道小孔质量泄漏速率计算公式，其中，式2是临界流动计算公式，式3为亚临界流动计算公式。0065临界流动与非临界流动是由临界压力比决定，见式400660067其中，。
33、RCP为临界压力比，无因次。0068当时，为亚临界流动，当时，为临界流动。0069由于泄漏点深度已知，则可知泄漏处出油管压力与温度，且认为其在泄漏阶段保持不变；泄漏点处环空压力可根据井口环空压力与保护液液柱压力计算得出，只要有一个井口环空压力值，则可求出对应的泄漏点处环空压力；此外，式2和式3中除泄漏面积外，其余均可容易获取，部分为系数，部分为常量，部分为气井基本参数，仅泄漏面积未知。00703确定任意泄压时间段内泄漏面积0071假设在单次环空泄压过程中，泄漏点处泄漏面积为定值，则可根据式2和式3分别得出任意时刻I泄漏点处的质量泄漏速率，即式5和式6；007200730074其中，QMI为I时。
34、刻泄漏点处的质量泄漏速率，KG/S；P0I为I时刻泄漏点处环空压力，PA；根据式5和式6，采用环空泄压数据以及气井基础数据可得出各个时刻的泄漏点处质量泄漏速率表达式，且其是泄漏面积A的函数，必须确定一个边界条件，得出泄漏面积A。0075以上得出的都是质量泄漏速率，可通过式7，将质量泄漏速率转换成体积泄漏速率，00760077其中，QVI为I时刻泄漏点处的体积泄漏速率，M3/S；GI为I时刻泄漏点处天然气密度，KG/M3。说明书CN104131811A118/10页120078根据积分原理，对泄漏点处体积泄漏速率求关于时间的积分，即可得到井下任意时间段T内的体积泄漏量，即式8，00790080其。
35、中，Q为任意时间段内井下累计体积泄漏量，M3；为等效积分量，M。表示式8中积分式中除泄漏面积A外的积分数值，在任意给定的时间段内，为定值。0081当泄压平衡时，井下泄漏与井口泄漏标况下达到平衡，由于天然气在井下和井口标况下的性质不同，两者的体积相差一个天然气体积膨胀系数，EG，即式9。00820083其中，EG为天然气体积膨胀系数，M3/M3；Z为天然气体积偏差因子，无因次。0084由此可得平衡等式10。0085VSEGVEGA100086其中，V表示为任意时间段内泄漏点处累计体积泄漏量，M3；VS表示为任意时间段内标况下累计体积泄漏量，M3。0087如果现场可以测得任意时间段T内井口累计体积。
36、泄漏量VS，则可以确定泄漏面积A；如果没有井口累计体积泄漏量VS数据，可以根据泄压数据中的平衡压力，选用式2或式3对应的质量泄漏速率模型确定质量泄漏速率表达式，并转换成体积泄漏速率表达式；由于泄压平衡时，井口环空压力与大气压均为定值，泄压管道横截面积已知，可以计算任意时间段T标况下累计的体积泄漏量，由此确定泄漏面积A。0088此外，现场每次环空泄压后，恢复的环空压力值并不是保持在泄压前的水平，认为是A发生变化，导致泄漏速率发生变，进而反过来影响环空压力。认为在每次泄压时，A保持恒定，可以通过在下次泄压曲线，得出变化后的A。00894确定环空泄漏速率计算表达式0090根据得出的泄漏面积，则可以得。
37、出气井标况下体积泄漏速率计算模型，即式11和式12。009100920093可以根据式11和式12得出的任意时间段内气井标况下的体积泄漏速率，判断气井的完整性状态，给予现场的气井生产运行及管理提出可靠的建议。0094实施例0095以XX气田XX1井为例为例，如图3所示，为本实施例的气井泄漏模型示意图，其中标注了各参数的代表的意义。00961确定泄漏点深度；说明书CN104131811A129/10页130097如图4所示，为本实施例的XX1井环空泄压压力恢复曲线图。需要根据该曲线确定环空压力平衡稳定阶段，进而可以确定环空泄漏点深度。根据气井的基础数据以及油套环空泄压数据，确定平衡环空压力，优选。
38、HAGEDORNBROWN方法计算气井的温度、压力剖面。XX1井泄压平衡时，井口环空压力为34710000MPA；当压力恢复平衡时，即泄漏点处油管内外压力相等，此时井口环空压力为64420000MPA。根据计算结果，回归温度、压力剖面的计算公式0098T107H20018H4504130099P12103H24100H86775000140100P0644200001200GH150101根据式1、式14和式15可得泄漏点深度为292063M，泄漏点处油管内压力为98766088PA，温度为37380K。01022确定泄漏点处体积泄漏速率表达式0103根据油管内泄漏点处气体压力以及式1可得出不。
39、同时刻泄漏点处环空压力，根据式4判断泄漏的流动状态。由于泄漏点深度已知，则可知泄漏点出油管内压力与温度，且认为其在泄漏阶段保持不变；泄漏点处环空压力可根据井口环空压力与保护液液柱压力计算得出，只要有一个井口环空压力值，则可求出对应的泄漏点处环空压力；此外，式2和式3中除泄漏面积外，其余均可容易获取，部分为系数，部分为常量，部分为气井基本参数，仅泄漏面积未知。实际运算过程中可进行编程计算。01043确定任意泄压时间段内泄漏点处的泄漏面积0105根据上一步确定的泄漏点处体积泄漏速率表达式，按照式7进行转换可以得出泄漏点处的体积泄漏速率，按照式8对泄漏点处体积泄漏速率进行积分，可得出对应时间内的体积。
40、泄漏量表达式，该表达式式关于泄漏面积的函数；按照式10转换可得标况下体积泄漏量。如果对应时间内体积泄漏量已知，则可得出泄漏面积。0106本例中，未取得对应时间内的井口体积泄漏量，可采用井口泄压的泄漏流量计算体积泄漏量。由图4可知，当泄压平衡时，环空压力基本保持不变，为34710000PA，可求取T15H至T18H之间的总泄漏量。本例中采用1/4“管线泄压，壁厚为225MM，则内径为185MM，面积为268106M2；取098；K取128；CD取10；R取8314JMOLK；PA取泄压平衡时井口环空压力，为34710000PA；TA取井口处温度，为3182K；PA和TA可根据泄压数据直接得出；M。
41、可根据天然气组分的体积百分比以及各组分的摩尔质量求出，本例中为0017628KG/MOL；Z可根据采用目前成熟的图版法和计算法得出，也可以采用气田的体积偏差因子，本例中取201。由于此时为临界流动，采用式5可得质量泄漏速率。当T15H至T18H时，井口泄压泄漏为临界流动，体积泄漏量见式16，01070108此时，EG为26396M3/M3，天然气密度为0803987KG/M3。0109由此可得泄漏点处泄漏面积为3877109M3，等效泄漏半径为3514105M。01104确定标况下体积泄漏速率表达式0111根据得出的泄漏点处泄漏面积和式11及式12，可得XX1井标况下体积泄漏速率表达式。最终得。
42、出XX1的最大标况下体积泄漏速率为824M3/MIN。该方法为气田现说明书CN104131811A1310/10页14场评价判断气井的完整性状态，为气井生产运行及管理提出可靠的建议。0112综上所述，本发明实施例在高温、高压、含酸性介质气井中的完整性评价中发挥重要作用。该方法仅依靠现场已有数据经过理论模型计算，可确定气井油套环空泄漏面积，最终得出油套环空标况下体积泄漏速率计算模型，实现基于气井基本数据的标况下体积泄漏速率确定。0113此外，本发明实施例还创新性提出利用气井现场最基本的已有数据确定油套环空泄漏面积和标况下体积泄漏速率的计算方法，对气井的完整性评价及风险评估有重要指导意义。0114以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。说明书CN104131811A141/3页15图1说明书附图CN104131811A152/3页16图2说明书附图CN104131811A163/3页17图3图4说明书附图CN104131811A17。