CN112241843A - 一种带腐蚀缺陷油套管剩余强度评价方法 - Google Patents
一种带腐蚀缺陷油套管剩余强度评价方法 Download PDFInfo
- Publication number
- CN112241843A CN112241843A CN202011124087.5A CN202011124087A CN112241843A CN 112241843 A CN112241843 A CN 112241843A CN 202011124087 A CN202011124087 A CN 202011124087A CN 112241843 A CN112241843 A CN 112241843A
- Authority
- CN
- China
- Prior art keywords
- oil
- residual
- oil casing
- formula
- corrosion
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Development Economics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Operations Research (AREA)
- Strategic Management (AREA)
- Pure & Applied Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Mathematical Physics (AREA)
- Economics (AREA)
- Educational Administration (AREA)
- Computational Mathematics (AREA)
- Geometry (AREA)
- Algebra (AREA)
- Evolutionary Computation (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Computer Hardware Design (AREA)
- Game Theory and Decision Science (AREA)
- Probability & Statistics with Applications (AREA)
- Marketing (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
Abstract
本发明提供一种带腐蚀缺陷油套管剩余强度评价方法,收集油套管实际参数,推导建立剩余强度计算方程,建立剩余强度公式;结合蒙特卡罗数学统计方法,建立腐蚀油套管状态函数;建立腐蚀油管失效模型,对随机变量分别进行一次模拟可以得到一个二维矩阵,把每次模拟得到的二维矩阵进行模拟计算,根据多次模拟结果可知,该油套管剩余强度可靠度。本发明基于常规评价标准,推导油套管剩余强度公式,结合蒙特卡罗数学统计方法,建立多元评价函数,可对多个影响指标同时进行评价,同时,通过数理统计的方法,大大提高了评价的准确度与高效性,适用性较好。
Description
技术领域
本发明涉及一种带腐蚀缺陷油套管剩余强度评价方法,属于钻完井工艺领域。
背景技术
随着我国很多油田进入后期开采阶段,油管腐蚀问题越来越成为影响油气生产安全的主要因素之一,对腐蚀后的油管进行全面系统可靠性评价变得越来越迫切。对于油管的腐蚀评价非常困难,难以像输油管道和其它压力容器那样得出明确的直接可以应用的评价标准。找到一种准确、高效的腐蚀缺陷油套管剩余强度评价方法是油气井长期安全生产的重点。
经过大量调研,现有的油套管剩余强度评价方法均通过行业标准进行,其步骤繁琐,流程复杂,且判别参数及标准参差不齐,大多评价方法中用到的评价参数为1至2个,评价指标较少,评价结果误差较大,使用局限性较大。
发明内容
本发明目的是:为了解决现今油套管剩余强度评价指标单一,评价结果误差较大,使用局限性大等问题,本发明蒙结合特卡罗数学统计方法,可同时实现多元指标的协同判别,评价结果较准确,可推广性强。
为实现上述目的,本发明提供了一种带腐蚀缺陷油套管剩余强度评价方法,包括以下步骤:收集油套管实际参数,建立剩余强度公式;结合蒙特卡罗数学统计方法,建立腐蚀油套管状态函数;建立腐蚀油管失效模型,求出概率密度函数及其概率分布函数;重复进行N次模拟,得到失效次数M;根据N次模拟结果可知,该油套管剩余强度可靠度为R=(N-M)/N。
具体的,该方法包括下列步骤:
S1,收集油套管实际参数,推导建立剩余强度计算方程,结合腐蚀深度d,油套管钢材流变应力σf,建立剩余强度公式;
S2,根据建立的剩余强度公式,结合蒙特卡罗数学统计方法,建立腐蚀油套管状态函数式中:P为油套管材许用应力,单位为MPa;d为腐蚀深度,单位为mm;t为实测壁厚,单位为mm;σf油套管钢材流变应力,单位为MPa;D为油套管外径,单位为mm;Mt为傅里叶因子,无量纲量;
S3,建立腐蚀油管失效模型Pf=∫z<0…∫f(d,L,t,D,σz,σb,P)dddLdDdσsdσbdP,分别求出d,σs,P的概率密度函数f(xi)及其概率分布函数F(xi),式中f(d,L,t,D,σs,σb,P)是d,L,t,D,σs,σb,P综合的概率密度函数,d为腐蚀深度,单位为mm;L为均厚长度,单位为mm;t为实测壁厚,单位为mm;D为油套管外径,单位为mm;σs为油管钢材的屈服强度,单位为MPa;σb为油管钢材的抗拉强度,单位为MPa;P为油套管材许用应力,单位为MPa;
S4,利用关系式生成一个均匀分布在(0.1)区间中的随机数xij,对给定的F(xij),对n个随机变量分别进行一次模拟可以得到一个1×n的二维矩阵(x1j,x2j,L,xnj),式中:i为随机变量个数,j为模拟的次数;
S6,根据N次模拟结果可知,该油套管剩余强度可靠度为R=(N-M)/N。
其中,步骤S1中的所述建立剩余强度公式包括以下步骤:
获取油套管基础参数,包括:实测壁厚t、最小实测厚度tmm,腐蚀余量FCA,最小许用壁厚tmin,油套管外径D,油套管内径Dt,缺陷深度d,油套管钢材流变应力σf,油管钢材的屈服强度σs;
根据剩余厚度比计算公式壳体系数计算公式傅里叶因子计算公式Mt=(1+0.48λ2)0.5,分别计算剩余厚度比、壳体系数以及傅里叶因子;式中:tmm为最小实测厚度,单位为mm;Rt为剩余厚度比,无量纲量;FCA为腐蚀余量,单位为mm;tmin为最小许用厚度,单位为mm;λ为壳体系数,无量纲量;L为均厚长度,单位为mm;D为油套管外径,单位为mm;t为实测壁厚,单位为mm;Mt为傅里叶因子,无量纲量;
本发明提供的一种带腐蚀缺陷油套管剩余强度评价方法,主要针对油套管均匀腐蚀情况,但针对局部腐蚀情况仍具有适用性。
本发明基于常规评价标准,推导油套管剩余强度公式,结合蒙特卡罗数学统计方法,建立多元评价函数,可对多个影响指标同时进行评价,同时,通过数理统计的方法,大大提高了评价的准确度与高效性,适用性较好。
与现有技术相比,本发明具有以下有益效果:(1)实现多元指标协同评价,评价结果较准确;(2)结合蒙特卡罗数学统计方法,使用范围广,可操作性高;(3)可推广性强。
附图说明
图1是本发明的技术路线图。
具体实施方式
下面结合实施方式和附图对本发明做进一步说明。
为实现上述目的,本发明提供了一种带腐蚀缺陷油套管剩余强度评价方法,其实施步骤如图1,该方法包括下列步骤:
S1,根据油套管实际参数,推导建立剩余强度计算方程,结合腐蚀深度d,油套管钢材流变应力σf,建立剩余强度公式;
S2,根据建立的剩余强度公式,结合蒙特卡罗数学统计方法,建立腐蚀油套管状态函数式中:P为油套管材许用应力,单位为MPa;d为腐蚀深度,单位为mm;t为实测壁厚,单位为mm;σf油套管钢材流变应力,单位为MPa;D为油套管外径,单位为mm;Mt为傅里叶因子,无量纲量;
S3,建立腐蚀油管失效模型Pf=∫z<0…∫f(d,L,t,D,σz,σb,P)dddLdDdσsdσbdP,分别求出d,σs,P的概率密度函数f(xi)及其概率分布函数F(xi),式中f(d,L,t,D,σs,σb,P)是d,L,t,D,σs,σb,P综合的概率密度函数,d为腐蚀深度,单位为mm;L为均厚长度,单位为mm;t为实测壁厚,单位为mm;D为油套管外径,单位为mm;σs为油管钢材的屈服强度,单位为MPa;σb为油管钢材的抗拉强度,单位为MPa;P为油套管材许用应力,单位为MPa;
S4,利用关系式生成一个均匀分布在(0.1)区间中的随机数xij,对给定的F(xij),对n个随机变量分别进行一次模拟可以得到一个1×n的二维矩阵(x1j,x2j,L,xnj),式中:i为随机变量个数,j为模拟的次数;
S6,根据N次模拟结果可知,该油套管剩余强度可靠度为R=(N-M)/N。
进一步的,所述建立剩余强度公式包括以下步骤:
获取油套管基础参数,包括:实测壁厚t、最小实测厚度tmm,腐蚀余量FCA,最小许用壁厚tmin,油套管外径D,油套管内径Dt,缺陷深度d,油套管钢材流变应力σf,油管钢材的屈服强度σs;
根据剩余厚度比计算公式壳体系数计算公式傅里叶因子计算公式Mt=(1+0.48λ2)0.5,分别计算剩余厚度比、壳体系数以及傅里叶因子;式中:tmm为最小实测厚度,单位为mm;Rt为剩余厚度比,无量纲量;FCA为腐蚀余量,单位为mm;tmin为最小许用厚度,单位为mm;λ为壳体系数,无量纲量;L为均厚长度,单位为mm;D为油套管外径,单位为mm;t为实测壁厚,单位为mm;Mt为傅里叶因子,无量纲量;
进一步的,所述方法主要针对油套管均匀腐蚀情况,但针对局部腐蚀情况仍具有适用性。
基于一种带腐蚀缺陷油套管剩余强度评价方法,接下来以某气藏A井套管剩余强度评价为例,进行该方法的进一步说明。
建立腐蚀油管失效模型Pf=∫z<0…∫f(d,L,t,D,σz,σb,P)dddLdDdσsdσbdP,分别求出d,σs,P的概率密度函数f(xi)及其概率分布函数F(xi),式中f(d,L,t,D,σs,σb,P)是d,L,t,D,σs,σb,P综合的概率密度函数;
获取油套管某一点基础参数,详见表1,包括:实测壁厚t、最小实测厚度tmm,腐蚀余量FCA,最小许用壁厚tmin,油套管外径D,油套管内径Dt,缺陷深度d,油套管钢材流变应力σf,油管钢材的屈服强度σs;
表1油套管某一点基础参数
实测壁厚t | 4.68 | 最小许用壁厚t<sub>min</sub> | 4 |
最小实测厚度t<sub>mm</sub> | 4.38 | 油套管外径D | 139 |
缺陷深度d | 0.84 | 许用应力P | 20 |
腐蚀余量FCA | 0.68 | 油套管内径D<sub>t</sub> | 127 |
钢材流变应力σ<sub>f</sub> | 216 | 钢材屈服强度σ<sub>s</sub> | 426 |
管道均匀减薄厚度LOSS | 0.53 | 剩余壁厚系数Q | 0.9 |
根据上述数据,计算得到壳体系数、傅里叶因子,均厚长度、剩余强度如表2所示:
壳体系数λ | 0.95 | 傅里叶因子M<sub>t</sub> | 1.198 |
均厚长度L | 18.89 | 剩余套管强度P<sub>e</sub> | 14.04 |
Z=Pe-P=-5.96(<0,套管失效)
重复上述步骤,选取第二点所测数据(x1j,x2j,L,x7j)进行判别,共重复1000次上述取点及计算步骤,失效次数为737次,可得该套管的剩余强度可靠度为:
R=(1000-737)/1000=26.3%
综上所述,结合蒙特卡罗统计学方法,得到该套管的剩余强度可靠度为26.3%,可靠度较低,失效风险高,需及时采取相应修井措施。
最后所应说明的是:以上实施例仅用以说明而非限制本发明的技术方案,尽管参照上述实施例对本发明进行了详细说明,本领域的普通技术人员应该理解:依然可以对本发明进行修改或者等同替换,而不脱离本发明的精神和范围的任何修改或局部替换,其均应涵盖在本发明的权利要求范围当中。
Claims (3)
1.一种带腐蚀缺陷油套管剩余强度评价方法,其特征在于,包括以下步骤:收集油套管实际参数,建立剩余强度公式;结合蒙特卡罗数学统计方法,建立腐蚀油套管状态函数;建立腐蚀油管失效模型,求出概率密度函数及其概率分布函数;重复进行N次模拟,得到失效次数M;根据N次模拟结果可知,该油套管剩余强度可靠度为R=(N-M)/N。
2.根据权利要求1所述的一种带腐蚀缺陷油套管剩余强度评价方法,其特征在于:具体包括以下步骤:
S1,收集油套管实际参数,推导建立剩余强度计算方程,结合腐蚀深度d,油套管钢材流变应力σf,建立剩余强度公式;
S2,根据建立的剩余强度公式,结合蒙特卡罗数学统计方法,建立腐蚀油套管状态函数Z:
式中:P为油套管材许用应力,单位为MPa;d为腐蚀深度,单位为mm;t为实测壁厚,单位为mm;σf油套管钢材流变应力,单位为MPa;D为油套管外径,单位为mm;Mt为傅里叶因子,无量纲量;
S3,建立腐蚀油管失效模型:
Pf=∫z<0…∫f(d,L,t,D,σz,σb,P)dddLdDdσsdσbdP
分别求出d、σs,P的概率密度函数f(xi)及其概率分布函数F(xi),式中f(d,L,t,D,σs,σb,P)是d、L、t、D、σs、σb、P综合的概率密度函数,d为腐蚀深度,单位为mm;L为均厚长度,单位为mm;t为实测壁厚,单位为mm;D为油套管外径,单位为mm;σs为油管钢材的屈服强度,单位为MPa;σb为油管钢材的抗拉强度,单位为MPa;P为油套管材许用应力,单位为MPa;
S4,利用关系式生成一个均匀分布在(0.1)区间中的随机数xij,对给定的F(xij),对n个随机变量分别进行一次模拟可以得到一个1×n的二维矩阵(x1j,x2j,L,xnj),式中:i为随机变量个数,j为模拟的次数;
S6,根据N次模拟结果可知,该油套管剩余强度可靠度为R=(N-M)/N。
3.根据权利要求2所述的一种带腐蚀缺陷油套管剩余强度评价方法,其特征在于:步骤S1中所述建立剩余强度公式,包括以下步骤:
获取油套管基础参数,包括:实测壁厚t、最小实测厚度tmm,腐蚀余量FCA,最小许用壁厚tmin,油套管外径D,油套管内径Dt,缺陷深度d,油套管钢材流变应力σf,油管钢材的屈服强度σs;
根据剩余厚度比计算公式壳体系数计算公式傅里叶因子计算公式Mt=(1+0.48λ2)0.5,分别计算剩余厚度比、壳体系数以及傅里叶因子;式中:tmm为最小实测厚度,单位为mm;Rt为剩余厚度比,无量纲量;FCA为腐蚀余量,单位为mm;tmin为最小许用厚度,单位为mm;λ为壳体系数,无量纲量;L为均厚长度,单位为mm;D为油套管外径,单位为mm;t为实测壁厚,单位为mm;Mt为傅里叶因子,无量纲量;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011124087.5A CN112241843A (zh) | 2020-10-20 | 2020-10-20 | 一种带腐蚀缺陷油套管剩余强度评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011124087.5A CN112241843A (zh) | 2020-10-20 | 2020-10-20 | 一种带腐蚀缺陷油套管剩余强度评价方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112241843A true CN112241843A (zh) | 2021-01-19 |
Family
ID=74169235
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011124087.5A Pending CN112241843A (zh) | 2020-10-20 | 2020-10-20 | 一种带腐蚀缺陷油套管剩余强度评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112241843A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114492001A (zh) * | 2022-01-14 | 2022-05-13 | 国家石油天然气管网集团有限公司 | 一种山地管道凹陷划痕组合缺陷的评价方法 |
CN114528700A (zh) * | 2022-01-26 | 2022-05-24 | 西安三维应力工程技术有限公司 | 一种含腐蚀坑油管剩余强度的确定方法 |
-
2020
- 2020-10-20 CN CN202011124087.5A patent/CN112241843A/zh active Pending
Non-Patent Citations (4)
Title |
---|
吴迪: "基于马尔可夫链蒙特卡洛法的天然气管道腐蚀可靠性研究", 《中国优秀硕士学位论文全文数据库 工程科技I辑》 * |
张占奎: "油气管道腐蚀失效预测及安全可靠性评估研究", 《中国博士学位论文全文数据库 工程科技I辑》 * |
张川: "含腐蚀损伤金属材料剩余寿命与剩余强度研究", 《中国博士学位论文全文数据库 工程科技I辑》 * |
蔡文军 等: "腐蚀管线剩余强度评估的研究进展", 《石油机械》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114492001A (zh) * | 2022-01-14 | 2022-05-13 | 国家石油天然气管网集团有限公司 | 一种山地管道凹陷划痕组合缺陷的评价方法 |
CN114492001B (zh) * | 2022-01-14 | 2024-08-02 | 国家石油天然气管网集团有限公司 | 一种山地管道凹陷划痕组合缺陷的评价方法 |
CN114528700A (zh) * | 2022-01-26 | 2022-05-24 | 西安三维应力工程技术有限公司 | 一种含腐蚀坑油管剩余强度的确定方法 |
CN114528700B (zh) * | 2022-01-26 | 2024-04-09 | 西安三维应力工程技术有限公司 | 一种含腐蚀坑油管剩余强度的确定方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112241843A (zh) | 一种带腐蚀缺陷油套管剩余强度评价方法 | |
CN108166963B (zh) | 一种海上油气井压裂效果评价方法 | |
CN105404776B (zh) | 一种含腐蚀缺陷的管道的可靠度评价方法 | |
CN113468757B (zh) | 一种基于模糊随机理论评估腐蚀天然气管道可靠性的方法 | |
CN107966545B (zh) | 一种任意成熟度烃源岩的有机碳恢复方法及其装置 | |
CN108920421A (zh) | 一种点腐蚀作用下油气管道系统的随机失效概率评估方法 | |
CN109736720B (zh) | 一种基于改进Kriging模型的深海连接器密封结构优化方法 | |
CN113806964B (zh) | 一种考虑多因素耦合作用的腐蚀与结垢速率预测方法 | |
CN111539142B (zh) | 一种管道裂纹扩展驱动力计算方法及系统 | |
CN109459302A (zh) | 一种腐蚀和高温联合作用下套管强度校核与优化设计方法 | |
CN112347414A (zh) | 一种单缺陷弯管极限内压计算方法 | |
CN109977563B (zh) | 基于tv-hsmm的海底油气管道剩余寿命预测方法 | |
CN111881619A (zh) | 基于matlab工具箱的bp神经网络实现管件冲蚀缺陷的预测方法 | |
CN114528700B (zh) | 一种含腐蚀坑油管剩余强度的确定方法 | |
CN115222053B (zh) | 一种基于蒙特卡罗方法的地层可压性概率分析方法 | |
CN111413558A (zh) | 一种基于dbscan的变压器故障诊断方法 | |
CN109784590B (zh) | 一种基于cagm(1,1)-bpnn的在役油气管道腐蚀预测方法 | |
CN112949190A (zh) | 基于r6-fad和bp-mcs的含裂纹缺陷管道的模糊随机可靠度评估方法 | |
CN112347583B (zh) | 一种增压站含双缺陷弯管极限内压的计算方法 | |
CN112364540A (zh) | 一种含双缺陷同心变径管极限内压的计算方法 | |
CN117688732A (zh) | 一种考虑多失效竞争与变量相关性的管道失效概率预测方法 | |
CN112966023B (zh) | 一种井筒的完整性预判方法 | |
CN102609602A (zh) | 多消防流量测试工况约束下给水管网管道粗糙度校核方法 | |
Wang et al. | A Novel Approach for Assessing the Accuracy of Prediction Models | |
CN118070458A (zh) | 数字孪生驱动的燃气管道剩余强度评估方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210119 |
|
RJ01 | Rejection of invention patent application after publication |