某海上气田井筒完整性失效概率计算.pdf

五油管材与仪器68PETROLEUM TUBULAR GOODS&INSTRUMENTS2023年10 月完整性技术某海上气田井筒完整性失效概率计算*侯冠中,付建民,廖前华,孙连坡,张智,赵苑瑾（1.中海石油（中国）有限公司天津分公司天津30 0 459；2.西南石油大学油气藏地质及开发工程国家重点实验室四川成都6 10 50 0）摘要：某海上气田属于“中等高含二氧化碳、微含硫化氢”凝析气藏，长期生产时套管易发生腐蚀失效，破坏井筒完整性，诱发环空带压，不仅严重影响经济效益，还会危及到海洋的生态环境和工作人员的人身安全。因此，识别该气田潜在的安全隐患，建立完备的完整性评价方法具有重要意义。为此，系统调研了井筒完整性相关标准的发展现状，结合渤中该气田的特点，基于故障树分析法识别了井筒完整性失效的3个主要模式和34个因素，并通过故障树分析法与贝叶斯网络的映射关系，得到了气田井筒完整性失效的贝叶斯网络模型。最后以井X进行了实例计算，得到了其井筒完整性失效概率，并通过敏感因素分析，反向推理了该井并筒完整性失效的关键因素，并通过对受时间影响的参数进行调整，最终实现了气并筒完整性失效概率的动态预测。研究结果表明：X并服役第3年的井筒完整性失效概率为0.0 51，在第8 年达到0.16 7,主要影响因素为腐蚀、环空带压，在生产中应当加强管控。关键词：海上气田；井筒完整性；故障树；贝叶斯网络模型；失效概率中图法分类号：TE58D0I:10.19459/ki.61-1500/te.2023.05.014Wellbore Integrity Failure Probability Calculation of an Offshore Gas Field2.State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation,Southwest Petroleum University,Chengdu,Sichuan 610500,China)Abstract:An offshore gas field is a condensate gas reservoir with medium of high carbon dioxide content and slight hydrogen sulfide con-tent.During long-term production,casing is prone to corrosion failure,damaging wellbore integrity and inducing annulus pressure,whichnot only seriously affects economic benefits,but also endangers marine ecological environment and personal safety of workers.Therefore,itis of great significance to identify the potential safety hazards of this gas field and establish a complete integrity evaluation method.There-fore,this paper systematically investigated the development status of wellbore integrity related standards,combined with the characteristicsof the gas field,identified three main modes and 34 factors of wellbore integrity failure based on FTA(Failure Tree Analysis),and obtainedthe Bayesian network model of wellbore integrity failure of the gas field through the mapping relationship between FTA and Bayesian net-work.Finally,a case study of well X was conducted to obtain its wellbore integrity failure probability.Through the analysis of sensitive fac-tors,the key factors of wellbore integrity failure of the well were inversely deduced,and the dynamic prediction of wellbore integrity failureprobability was finally realized by adjusting the parameters affected by time.The research results show that the wellbore integrity failureprobability of Well X in the third year of service is 0.051,and it reaches 0.167 in the eighth year.The main influencing factors are corro-文献标识码：AHOU Guanzhong,FU Jianmin,LIAO Qianhua,SUN Lianpo,ZHANG Zhi?,ZHAO Yuanjin?(1.Tianjin Branch,CN00C（Ch i n a)Co.Lt d.,T i a n j i n 30 0 459,Ch i n a;文章编号：2 0 96-0 0 7 7(2 0 2 3)0 5-0 0 6 8-0 6初投稿收稿日期：2 0 2 3-0 3-0 4；修改稿收稿日期：2 0 2 3-0 5-17*基金项目：国家自然科学基金“四川盆地深层超深层气井环空带压预防与管控基础研究”（编号U22A20164）；国家自然科学基金“深层页岩气长位移水平井套管完整性研究”（编号52 0 7 42 34）；四川省青年科技创新研究团队专项计划“油气井建井安全四川省青年科技创新研究团队”（编号2 0 2 0 JDTD0016）；中国石油西南石油大学创新联合体科技合作项目“深井/超深井/水平井安全高效建井关键基础理论与技术研究”（编号：2 0 2 0 CX040100）第一作者简介：侯冠中，男，198 2 年生，工程师，2 0 0 5年毕业于大庆石油学院资源勘察专业，主要从事海洋石油钻井技术监督与管理工作。E-mail:hougz 通迅作者简介：张智，男，197 6 年生，教授，博士生导师，2 0 0 5年毕业于西南石油大学油气工程专业，获博士学位，主要从事油气井工程、井筒完整性与环空带压管控、材料腐蚀与防腐等方面的研究工作。E-mail：w is e z h 12 6.c o m2023年第9 卷第5期sion and annulus pressure,and management and control should be strengthened in production.Key words:offshore gas field;wellbore integrity;fault tree;Bayesian network;failure probability0引言为实现稳产30 0 0 万吨、上产40 0 0 万吨的生产目标，我国加大了对海上气田的油气勘探开发力度，但某海上气田钻遇多套地层，目的层属高温地层，压力体系复杂，钻井过程溢流漏失频发，地层流体含有腐蚀性气体 ,如此严苛的自然环境给井筒完整性安全带来了巨大的挑战，同时由于该气田缺乏类似区块的井筒完整性管理经验，针对井筒完整性失效概率的研究尤为必要。随着技术的发展、安全法规的制定、大众对环境问题的重视，井筒完整性管理已经发生了重大变化。最初，井筒完整性主要通过常规检查和压力测试进行泄漏识别。然而，日益复杂的勘探开发环境对井筒完整性管理提出了更高的要求,包括使用先进的监测系统在问题发生之前进行检测并采取相应措施预防事故发生，或实施定期检查和屏障维护计划，以及开发用于保障井筒完整性的新材料等。现阶段,国内外学者越来越重视对井筒完整性的研究。2 0 16 年,Loizzo M2建立了保障井筒完整性的地层屏障设计方法;2 0 17 年,Kiran R等3 提出井筒完整性受到水泥质量、油套管腐蚀、生产设计、钻井或生产过程温压动态变化等的影响；2 0 19年，周松民等4 将海洋弃置井风险因素归为井口装置、井筒、管柱及其他4类;2 0 2 0 年，张智等5 建立了涵盖井屏障设计、腐蚀控制、寿命预测、作业参数控制等的注水井井筒完整性设计方法;2 0 2 1年，时维才等6 对管柱渗漏引起的油套环空压力升高开展了研究;2 0 2 1年，Yousuf 等7 研究了水泥环引起的井筒完整性失效问题。可见现阶段井筒完整性管理转向考虑因素更多、针对井型更加具体的方向发展，而且定量评价是完整性管理的重要基础。然而针对这类含腐蚀性气体的海上凝析气田，国内外专家学者尚未形成一套井筒完整性定量评价方法，无法准确计算其失效概率并提出有效的管理手段。近几年,故障树一贝叶斯网络在井筒完整性评价方面开展了较多应用：2 0 15年，李新宏等8 利用贝叶斯模型实现了对海底管道泄漏事故动态分析；2 0 2 0 年，张智等9 识别了井下安全阀失效类型,并利用贝叶斯网络和Sintef数据库,计算了风险因素的权重和失效概率；袁晓兵等10 1利用贝叶斯网络研究了水下采油树系统的剩余寿命;李中等1 建立了深水探井井筒完整性评价模型；刘毅播等12 1结合贝叶斯法则分析了储气库井环空压力静态失效风险。由此可见,故障树贝叶斯网络模型在侯冠中等：某海上气田井筒完整性失效概率计算油气风险识别与分析方面具备可行性。为此,本文通过国内外井筒完整性相关标准，结合现场实际情况，建立海上气田井筒完整性失效的贝叶斯网络模型,形成井筒完整性失效概率的计算方法。最后以X井进行了实例计算,预测了该井井筒完整性未来的失效概率，并利用贝叶斯网络模型的反向推理能力识别出井筒完整性失效的关键因素。1井筒完整性相关标准的发展历程挪威石油工业协会NORSOKD-O10标准将井筒完整性定义为：采用有效的技术、管理手段来降低开采风险，保证油气井在成功废弃前的整个开采期间的安全13。目前国际上广泛接受的井简完整性概念是综合运用技术、操作和组织管理的解决方案来降低井筒在全生命周期内地层流体不可控泄漏的风险14。随着油气勘探开发逐渐向深海、极地等极端环境发展，井筒完整性也越来越受到石油行业重视，图1为近年来不同机构组织发布井筒完整性相关标准的统计数据。由图1可知,现有的井筒完整性标准涵盖了钻井、建井、深水井等多方面内容，但我国井筒完整性相关的标准较少。12钻井与建井过程完整性10口海上井及深水井完整性80口陆上井特殊作业完整性高温高压油气井完整性6全生命周期井简完整性420美国API挪威OLF国际ISO标准标准图1井筒完整性相关标准统计数据根据现有海上井井筒完整性相关标准,结合以往海上油气井井筒完整性失效的案例,可将该气田井筒完整性失效划分为3种主要类型。1)井口装置完整性失效，常见的井口装置完整性失效形式包括导向基座失效、低压井口头失效、高压井口头失效、套管悬挂器失效、防磨补芯装置失效、环形密封总成失效，一旦油气通过井口装置泄漏至海底，将导致不可估量的损失152)套管完整性失效，由于该气田高含腐蚀性气体，且井底温度压力较高，套管容易发生腐蚀失效。3)水泥环完整性失效,水泥环完整性失效主要是由于水泥环切向拉伸破坏、塑性变形和剪切失效或胶结面失效导致的水泥环密封性失效，水泥环一旦丧失密封性,将会形成多种渗漏通道,如图2 所示16 69国内企标准业标准70水泥环图2 水泥环完整性失效后形成的渗漏通道五油管材与仪器套管地层2故障树模型(a)水泥内部缝隙本文应用故障树分析法（FailureTree Analysis，(b)套管外壁与水泥环之间缝隙FTA)以井筒完整性失效作为顶事件T,井口装置完整性失效、套管完整性失效、水泥环完整性失效分别为中间(c）水泥环与地层间缝隙事件AIvA2A3,结合现场实际,将中间事件的影响因素作为故障树的底事件，根据中间事件失效影响因素的逻(d)套管外壁与地层之间通缝辑关系建立其故障树模型,如图3所示，故障树中的符号各自对应的事件如表1所示。2023年10 月BB.BB。D(X3图3井筒完整性失效故障树表1井筒完整性失效故障树中的符号和相应事件符号TA1水下井口装置完整性失效A2套管完整性失效A3水泥环完整性失效BI法兰泄漏B2仪表泄漏B阀门泄漏B4油套管头泄漏Bs外部压力改变B6套管抗挤强度不够B7水泥环切向拉伸破坏B8水泥环塑性变形B,水泥环剪切失效C套管抗挤强度降低3贝叶斯网络模型贝叶斯网络17-18 （Bayesian Network,BN）是一种基于概率分析的现代不确定性推理方法，适用于缺乏统计事件井筒完整性失效符号C2套管自身抗挤强度低C3轴向拉力大C4抗拉强度小Cs横向错断力大C6抗剪强度低DI套管腐蚀X法兰腐蚀X2密封不严X操作失误X4压力表破裂X考克密封失效X6闸阀内漏X7节流阀节流失效Xs阀体与阀盖连接处泄漏事件符号XX10XIX12X13X14X15X16X17X18Xi9X20X21X22数据的概率计算,其数学表达式为B=（G,P）,通常包括n个节点组成的有向无环图（DirectedAcyclicGraph，DAG)及表示节点间影响程度的条件概率表CPT(Condi-tionProbabilityTable）。图4为5个节点组成的简单贝叶事件盘根漏气注脂单流阀泄漏油套管顶丝泄漏油套管连接处泄漏地壳运动地层断裂地震地层强度差地层出水开采时间过长套管偏磨材料抗腐蚀性能差腐蚀环境防腐措施失效符号X23X24X25X26X27X28X29X30X31X32X33X34事件固井质量差射孔不合理套管材质差套管设计不合理环空带压套管内泄漏地层温度高水泥设计压裂施工蠕变地应力作用地层倾角大海底滑坡2023年第9 卷第5期斯网络,箭头代表了节点间的因果关系，通常由父节点（原因)指向子节点（结果），没有父节点的节点通常称为根节点。父节点与子节点之间不是相互独立的,可由CPT表示两者间存在的影响程度,如在to时刻,系统中节点B发生引起A发生的概率为0.1节点B2发生引起A发生的概率为0.5;B1、B,同时发生引起A发生的概率为0.7;B1、B,均不发生情况下A发生的概率为0.0 5；则A节点的CPT如表2 所示。随着时间变化,BI、B,对A的影响也会发生变化,此时只需要更改节点A的CPT中相应数值即可，因此贝叶斯网络适用于动态失效概率分析。图4贝叶斯模型示例图表2 节点A的CPT条件概率BI0B20A=10.05A=00.95贝叶斯网络中所有节点状态值取值为0 或1,取0时代表一定不发生,取1时代表一定发生,则根节点外的节点可以由不同状态值的根节点组合表示，如节点A可以表示为(X=a1,X,=a2,.X,=a;,.X,=an),a,为 0或1,则可由式(1)计算某事件A的发生概率19。P(A)=P(XI,X2,X,)=IIP(X,I(X,)(1)式中:P(A)为事件A发生的概率,n为贝叶斯网络的根节点数，（X）为节点X的父节点集合。节点C发生的概率为：P(C=1I(C)sub)=1-,I式中：（C)为节点C的父节点集合;（C)s u b 为节点C的父节点集合的子集;B,为不同的B节点事件;P,为B,的发生概率。由于表1中只是列举了井筒完整性失效的主要影响因素,还存在一些未知因素的影响没有考虑,此时可以采用Leakage Noisy-ORgate模型,将这些未知因素视作一个父节点BL,其概率设为 PL,加人模型中,则 LeakageNoisy-ORgate模型中节点C取真值的概率为：P(C=1I(C)sub)=1-(1-PL),II、(1-P.)i:B,ET(C)sl(3)侯冠中等：某海上气田井筒完整性失效概率计算式中:PL为未知因素B,的发生概率。将子节点C的父节点分为B,和除B,之外的全部父节点集合B。,根据式(2），可得：P(C=1IB,)=1-(1-P)(1-P.)P(C=1IB,)=Pa因此,父节点B,的概率为：P(CIB,)-P(CIB,)P(i)=1-P(CIB,)式中：P.为B.发生概率。式(1)将联合概率分布的计算转化为条件概率乘积的形式，大大降低了贝叶斯网络模型的计算和推理过程。通过式(3)可以得到考虑未知因素下父节点的发生概率。在已知基础事件发生概率条件下,贝叶斯网络可B以计算出其他事件的发生概率,从而计算井筒完整性失效概率，同时可反向推理导致井筒完整性失效的主要因素。根据故障树与BN模型的映射关系,可以将图3的故障树模型转化为BN模型。通过现场保存的历史数据1和专家评估可以得到根节点的发生概率，和其他节点的100.50.10.50.9=1.II(1-P.)(2)i:B,eCn.71(4)(5)(6)1CPT,即可计算顶事件井筒完整性失效的发生概率,通过0.7不断更新根节点发生概率和其他节点CPT,可以对井筒0.3完整性失效概率进行动态分析，同时可以将顶事件发生概率定为1,利用贝叶斯网络的反向推理功能,识别影响顶事件的关键因素。4现场应用该气田X井太古界潜山压力系数1.115,温度梯度为3.5/10 0 m，属正常压力和温度系统。截止2 0 2 1年3月15日,X井累计生产了2 7 0 d,平均日产气1.4110 m/d，平均日产水1.91m/d,平均日产油14.2 8 m/d,累计产油3.810m，累计产水515.49m，累计产气3.8 10 m。该井完整性失效因素由图3中的故障树模型表示,包含了34个根节点，投产2 年后井筒完整性BN模型的根节点概率分布如表3所示。表3根节点的发生概率值符号发生概率X1.90 10-4X23.19 10-4X32.90 10-4X44.28 10-6Xs5.04 10-4X8.51 10-3X2.02 10-4符号XXX10oXIX12X13X14发生概率1.02 10-41.96 10-43.62 10-4 5.02 10-4 4.22 10-45.30 10-14.96 10-472符号Xi5X16X17XI8X19X20X21X22X23X24根据上文建立的贝叶斯网络模型可以由根节点发生概率逐级向上计算，最终得到其他节点的条件概率表，由于该模型中节点较多,文中只展示节点A3、BI、D,等部分节点的CPT,如表4表6 所示。表 4节点 A;的 CPT条件概率B70B80B,0A3=00.950.4170.3750.2680.520.215 0.2870.187A3=10.050.5830.6250.7320.480.7850.7130.813表5节点B,的CPT条件概率X0X20X0Bi=00.9910.8550.7630.650.9660.7690.6990.48B;=10.0090.1450.2370.350.0340.2310.3010.52表6节点D,的CPT条件概率X200X210D,=10.998D;=00.002在利用本文建立的故障树贝叶斯网络模型得到各节点条件概率表CPT后,将数据代人匹兹堡大学开发的GeNie2.0软件,计算出其他中间事件的失效概率如表7所示，最终顶事件发生概率为0.0 51。若将顶事件发生概率设为1,反向推理主要影响因素,可以识别影响X井五油管材与仪器续表A的5个关键影响因素：X21、X2 7、X2 s、X2 9、X31,其失效概发生概率符号2.62 10-4 X255.13 10-2X26X278.13 10-3X287.17 10-42.18 10-45.26 10-44.96 10-17.21 10-41.90 10-33.90 10-3101102023年10 月发生概率率分别为：0.496、0.6 2 6、0.46 1、0.531、0.0 0 5。因此，在1.80 10-2生产中加强对腐蚀、环空带压的防控以及优化压裂施工4.60 10-3方案是降低X井井筒完整性风险的有效途径。表7 其他中间事件的发生概率计算结果6.25 10-1节点发生概率节点发生概率节点发生概率4.60 10-1A1X295.30 10-1X302.60 10-3X315.26 10-3X321.02 10-4 X332.02 10-4X341.62 10-41101001010.5750.4250.163A20.185A30.209BI0.009B20.001B30.469随着该井生产的进行,某些根节点的发生概率不再变化，比如X23、X2 4等，但大部分根节点的发生概率及非根节点的CPT都会发生变化，因此通过更新这两类参数，并利用贝叶斯网络计算受这两类参数影响的其他节点条件概率表后，可以对井筒完整性失效风险进行动态预测,根据该井目前的条件和该区块的历史数据可预测未来5年内根节点的发生概率,输入建立的BN模型中，预测该井未来5年内的井筒完整性失效概率,如图5所1示。可知随着生产进行，该井井筒完整性失效概率先加101110100.5210.479B4BsB6B7B:Bg1速增加，但在第5年后，由于产气量下降和工作制度等的改变导致部分根节点失效概率减小，X井井筒完整性失效概率增加幅度降低。0.180.160.140.120.10.080.060.0510.040.02103图5井X的井筒完整性失效动态预测结果5 结 论为保障海上气田的安全生产，本文建立了适用于海1上气田井筒完整性评价的模型,并得出了以下结论。1)传统井筒完整性的分析过程需要大量的数据，难0.177以定量计算，而通过建立井筒完整性失效的FTA-BN模0.823型可以克服数据的不确定性，既可以通过正向推理计算该井井筒完整性失效概率，又能反向推理识别该井井筒完整性失效的关键影响因素，同时还可以通过更新相关数据，对失效概率进行动态预测，有效指导现场安全经济生产。2)实例X井服役第3年发生井筒完整性失效的概0.0510.3090.1980.2010.0540.0980.1650.1670.1520.1390.0874CiC2C3C4CsC656生产时间/a0.2570.070.6120.0540.0500.056782023年第9 卷第5期率为0.0 51,随着服役时间的增长，失效概率逐渐增大，但在第5年后，由于产气量下降和工作制度等的改变，导致部分根节点失效概率减小，井筒完整性失效概率增大的幅度减缓。生产中,应当加强对腐蚀、环空带压和压裂施工等因素的监管，以实现气井安全、经济生产。3)由于本文中基于故障树贝叶斯模型的失效概率评价方法计算精度较大程度取决于根节点的失效概率,在实际应用过程中，需要可靠的现场数据做支撑,并结合专家经验进一步提高评价结果与现场实际的相符性。参考文献1周长所,耿亚楠，雷杨,等.渤中19-6 复杂断块构造地应力分布规律试验研究J.石油化工应用,2 0 2 0,39(1)：8 4-8 9.2 LOIZZO M,LECAMPION B,MOGILEVSKAYA S.The role ofgeological barriers in achieving robust well integrityJ.Ener-gy Procedia,2017,114:5193 5205.3 KIRAN R,TEODORIU C,DADMOHAMMADI Y,et al.Iden-tification and evaluation of well integrity and causes of failureof well integrity barriers(a review)J.Journal of NaturalGas Science and Engineering,2017,45:511-526.4周松民,李军伟,李瑾,等.基于权重分析法的海洋弃置井井筒完整性评价方法建立J.钻采工艺,2 0 19,42(2)：10-13.5张智,杨昆,刘和兴,等注水井井筒完整性设计方法J.石油钻采工艺,2 0 2 0,42(1)：7 6-8 5.6时维才,姚峰,金勇,等.二氧化碳驱注采井管柱失效原因分析及防护对策J.石油管材与仪器,2 0 2 1,7(6)：7 4-7 8.7 YOUSUF N,OLAYIWOLA O,GUO B Y,et al.A comprehen-sive review on the loss of wellbore integrity due to cement fail-ure and available remedial methods J.Journal of PetroleumScience and Engineering,2021,207:109123.8李新宏，朱红卫,陈国明，等.海底油气管道泄漏事故风险+“侯冠中等：某海上气田井筒完整性失效概率计算分析的贝叶斯动态模型J.中国安全科学学报,2 0 15,2 5(4):75 80.9张智,赵苑瑾,吴优，等.页岩气井井下安全阀完整性研究J.石油管材与仪器,2 0 2 0,6(4):56-6 2.10袁晓兵,陈国明，范红艳，等.基于贝叶斯网络的水下采油树系统剩余寿命预测J.中国海上油气,2 0 2 0,32（4）：171-178.11李中,吴建树,闫伟，等.基于贝叶斯网络的深水探井井筒完整性失效风险评估J.中国海上油气，2 0 2 0,32（5）：120-128.12刘毅播,闫怡飞，宋运存，等.储气库套管柱环空压力影响因素动态评价J.石油管材与仪器,2 0 2 2,8（5）：51-56.13 NORSOK.Well integrity in drilling and well operations:NORSOK D-010:2021 S.No r w a y:No r s k So k k e lsKonkuranseposisjon,2021.14张绍槐.井筒完整性的定义、功能、应用及进展J.石油钻采工艺,2 0 18,40(1):1-8.15韦红术,闫晴，曹波波,等.基于FTA定性的简易水下井口设计可靠性分析J.石油矿场机械,2 0 2 2,51（2）：8-14.16 FORD E P,MOEINIKIA F,LOHNE H P,et al.Leakage cal-culator for plugged and abandoned wells C /SPE BergenOne Day Seminar,Bergen,Norway:SPE,2017:SPE-185890-MS.17骆正山,赵乐新，王小完.基于动态贝叶斯网络的海底管道点蚀疲劳损伤失效模型研究J.表面技术，2 0 2 0,49(1):269-275.18王海朋,段富海.复杂不确定系统可靠性分析的贝叶斯网络方法J.兵工学报,2 0 2 0,41(1):17 1-18 2.19陈坤,陈序,魏鑫，等.基于贝叶斯网络的页岩气井喷风险分析J.安全与环境学报，2 0 19,19（6）：18 7 6-18 8 3.（编辑：屈忆欣）.73(上接第 6 7 页)8戴联双,考青鹏,杨辉,等.高强度钢管道环焊缝隐患治理措施研究J.石油管材与仪器,2 0 2 0,6(2）：32-37.9帅健，张银辉.高钢级管道环焊缝的应变集中特性研究J.石油管材与仪器,2 0 2 0,6(2）：8-14.10】沙胜义.输油管道环焊缝缺陷疲劳寿命评估J.管道技术与设备,2 0 17(2):2 8-32.11任俊杰，马卫锋，惠文颖，等.高钢级管道环焊缝断裂行为研究现状及探讨J.石油工程建设,2 0 19,45（1）：1-5.12】杨锋平,许彦，曹国飞，等.X80管线钢环焊缝力学性能及缺陷成因分析J.压力容器，2 0 2 0,37(1)：6 8-7 3.13】包志刚，马夜，李超，等.焊材匹配对低温弯管焊缝组织和低温冲击韧性的影响J.焊管,2 0 2 0,43（12）：7-11.14张振永.高钢级大口径天然气管道环焊缝安全提升设计关键J.油气储运，2 0 2 0,39（7）：7 40-7 48.15苏怀，张劲军.天然气管网大数据分析方法及发展建议J.油气储运,2 0 2 0,39(10):10 8 1-10 9 5.16于听雷,张学珍,郑景云.中亚地区过去2 0 0 0 年气候变化研究的数据挖掘J.第四纪研究,2 0 2 1,41(2）：56 4-57 8.17 SONG T M,SONG J.Prediction of risk factors of cyberbul-lying-related words in Korea:application of data mining u-sing social big data J.Telematics and Informatics,2021,58:101524.18 HE B,YIN L J.Prediction modelling of cold chain logisticsdemand based on data mining algorithm J.MathematicalProblems in Engineering,2021,2021:3421478.(编辑：姜婷）