CN103136701A - Method for quantitatively evaluating pipe condensation probability of crude oil pipeline - Google Patents

Method for quantitatively evaluating pipe condensation probability of crude oil pipeline Download PDF

Info

Publication number
CN103136701A
CN103136701A CN2011103802161A CN201110380216A CN103136701A CN 103136701 A CN103136701 A CN 103136701A CN 2011103802161 A CN2011103802161 A CN 2011103802161A CN 201110380216 A CN201110380216 A CN 201110380216A CN 103136701 A CN103136701 A CN 103136701A
Authority
CN
China
Prior art keywords
pipeline
partiald
crude oil
temperature
pipe
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.)
Granted
Application number
CN2011103802161A
Other languages
Chinese (zh)
Other versions
CN103136701B (en
Inventor
苗青
姜保良
徐波
兰浩
王龙
闫峰
李其抚
薛向东
康阳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Oil and Gas Pipeline Network Corp
Original Assignee
China Petroleum and Natural Gas Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petroleum and Natural Gas Co Ltd filed Critical China Petroleum and Natural Gas Co Ltd
Priority to CN201110380216.1A priority Critical patent/CN103136701B/en
Publication of CN103136701A publication Critical patent/CN103136701A/en
Application granted granted Critical
Publication of CN103136701B publication Critical patent/CN103136701B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention discloses a method for quantitatively evaluating the pipe condensation probability of a crude oil pipeline, and relates to the technical field of pipeline systems. The method comprises the following steps of: 1) automatically extracting and analyzing the safety influence factor; 2) counting a data distribution rule of the influence factor; 3) calculating certainty shutdown and restart-up pressure; 4) setting an evaluation frequency and randomly sampling parameters; 5) calculating uncertainty shutdown and restart-up pressure; 6) counting a calculation frequency of restart-up failure; 7) obtaining a calculation result (the pipe condensation probability of the pipeline); 8) judging whether ineffectiveness judgment can be accepted; 9) if the ineffectiveness judgment cannot be accepted, adjusting a flowing guarantee scheme, turning to the step 2) of counting the data distribution rule of the influence factor, and if the ineffectiveness judgment can be accepted, submitting the flowing guarantee scheme. By the method, the crude oil pipeline flowing safety can be evaluated; and technical support can be supplied to safety and high-efficiency running of the pipeline.

Description

The solidifying pipe of crude oil pipeline probability method for quantitatively evaluating
Technical field
The present invention is the solidifying pipe of a kind of crude oil pipeline probability method for quantitatively evaluating, relates to the piping system technical field.
Background technology
One of final ultimate limit state of wax-containing oil pipe is exactly the solidifying tubulose attitude of knowing, and this is to need in conduit running to be avoided as possible, and the production loss that brings is also maximum.Generally, the main cause that this ultimate limit state occurs is after pipeline plan, especially accident stopping transportation, the solidifying pipe that start-up course unsuccessfully causes.Judge that at present common practice that whether pipeline coagulates pipe is the current inlet temperature according to pipeline, according to the surrounding enviroment situation, considers that in given shutdown time, whether crude oil temperature is greater than its condensation point.Though this way is simple and easy to use, science is not strong.At first, the method is not considered the uncertainty of conduit running technological parameter, and for example ground temperature is actual fluctuates; Secondly, the index of simply using condensation point is science not, and whether required pressure is obviously more objective greater than the maximum allowable operating pressure of pipeline when starting.
Summary of the invention
The objective of the invention is to invent a kind of based on the certainty and uncertainty method reliably, the solidifying pipe of crude oil pipeline probability method for quantitatively evaluating accurately.
The pipeline actual production is in service, its throughput rate, out-station temperature, steady-state operation overall heat transfer coefficient, buried depth of pipeline, buried depth of pipeline place ground temperature, soil thermal conductivity etc. all have uncertainty, evaluation has been considered the uncertainty of these factors take the actual production form as the basis.Estimate the automatic extraction and analysis of involved calculating parameter and completed by software, as throughput rate, pressure, inlet temperature, out-station temperature etc., and the accuracy of these data is checked.
The solidifying pipe of the crude oil pipeline that the present invention proposes method for quantitatively evaluating, mainly flow process is as shown in Figure 1:
1) robotization of safety effects factor is extracted and is analyzed;
2) the data distributivity rule of influence factor statistics;
3) determinacy stopping transportation restart pressure calculates;
4) estimate the setting of number of times and the grab sample of parameter;
5) calculating of uncertain stopping transportation restart pressure;
6) the statistics pipeline is restarted the calculation times of failure;
7) the drawing of result of calculation (the solidifying pipe of pipeline probability);
8) whether satisfy and to accept the judgement of losing efficacy;
9) if do not satisfy, after adjusting the flow assurance scheme, turn 2) the data distributivity rule statistics of influence factor;
If satisfy, submit the flow assurance overall plan to.
The pipeline actual production is in service, its throughput rate, out-station temperature, steady-state operation overall heat transfer coefficient, buried depth of pipeline, buried depth of pipeline place ground temperature, soil thermal conductivity etc. all have uncertainty, evaluation has been considered the uncertainty of these factors take the actual production form as the basis.Estimate the automatic extraction and analysis of involved calculating parameter and completed by software, as throughput rate, pressure, inlet temperature, out-station temperature etc., and the accuracy of these data is checked.Solidifying pipe definition of probability is as follows:
For a given ultimate limit state, probability analysis comprises the model of broad sense random load S and the random drag R of broad sense.Corresponding limit state function can be expressed as following form in general manner:
g(x)=R-S (1)
Obviously, indicate inefficacy as g (x)<0.
The definition failure probability is as follows:
P f = ∫ g ( x ) ≤ 0 f x ( x ) dx - - - ( 2 )
In formula, x is the stochastic variable vector; f x(x) be joint probability density function.
Restart problem for the crude oil pipeline stopping transportation, definition load S is the minimum starting pressure of pipeline, and getting drag R is pipeline maximum allowable operating pressure (MAOP), and Pf is solidifying pipe probability.
Described data distributivity probability statistics and check are:
According to ISO 16708, the production run data of quantitative evaluation pipeline collection are carried out uncertainty and value analysis.To manage in defeated process the pipeline flow data regularity of distribution as example, chosen the data on flows of three months flowmeter statistics for certain east pipeline.Result is as follows:
The statistical law of table 1 pipeline throughput rate data
Figure BDA0000112386140000022
The result of table 1 shows, the kurtosis of throughput rate data is close to 3, and deviation is close to 0, can think the Normal Distribution rule.The uncertainty analysis process of other mobile Safety Influence Factors as above.
Described determinacy stopping transportation restart pressure calculates:
As mentioned before, this evaluation method is estimated framework based on secondary, and first order framework is not in the situation that consider the fluctuation of production run situation, for satisfying on-the-spot demand, with the security of the existing production and operation condition of Deterministic Methods assessment, if dangerous, will take the flow assurance measure according to result.The determinacy operating mode is generally on-the-spot given production data.
The setting of described evaluation number of times and the grab sample of parameter are:
After having completed the deterministic parsing in system one-level framework, just can carry out large-scale probability sampling, the sampling number of times is defaulted as 10000 times, generates corresponding numerical value by its regularity of distribution function during the every sub-sampling of production run parameter, but must set the bound of production run parameter.What reason was the generation of regularity of distribution function is the numerical value that meets the mathematical distribution rule, but must not have physical significance.In the situation that there is no the production scene appointment, the bound of system default is as follows, and throughput rate is between 3 times of variances of its average plus-minus, and temperature out of the station is between 0 ℃-90 ℃, and ground temperature is between-20 ℃-50 ℃.
Described uncertain stopping transportation restart pressure calculates and starts frequency of failure calculating:
It is evaluation number of times according to having set that uncertain stopping transportation restart pressure calculates, and whether each estimates production run parameter that number of times generates and analyze solidifying pipe and occur.As start unsuccessfully and to occur to be recorded as 1, do not occur to record 0.After completing, evaluation can record the number of times that total stopping transportation is restarted failure.
Described interpretation of result is:
Solidifying tube failure probability P f based on calculating is defined as pipeline shutdown with (1-Pf) and restarts safety indexes (SRSI).Under constant condition of pipeline shutdown time, the SRSI graded index is larger, and the possibility of pipeline shutdown failure is less restarting just, and the safe coefficient of pipeline restart procedure is just higher; Otherwise the possibility that pipeline is restarted failure is just larger, according to evaluation result, to the security of operation proposition safeguard of pipeline.
The defeated technology Calculation of evaluation procedure pipe
Waterpower, heating power during normal the conveying calculate
(1) normally moving the axial temperature drop of hot oil pipeline calculates
Stablizing under transport condition, between the temperature drop rate that in hot oil pipeline, oil flows and the coefficient of heat conductivity of conduit running parameter (out-station temperature, throughput rate etc.), caliber, buried depth of pipeline and surrounding medium thereof, temperature, heating station, length and crude oil property (specific heat capacity, density, viscosity etc.) etc. are closely related.In pipe, oil stream can be similar to Л e й б e н з о Н formula along the journey Temperature Distribution and describe:
T L=(T 0+b)+(T R-T 0-b)e -aL (3)
In formula, a=K π D/Gc, b=ig/c pa; G is the mass rate of oil product, kg/s; C pBe the specific heat capacity of crude oil, J/ (kg ℃); D is outer diameter tube, m; L is the distance apart from heating station, m; K is the pipeline overall heat transfer coefficient, W/ (m 2℃); T RBe heating station departures oil temperature, ℃; T LBe the oil temperature of distance heating station distance for the L place, ℃; T 0Be ambient temperature, buried pipeline is got the natural ground temperature in tube hub buried depth place, ℃.
(2) frictional resistance of hot oil pipeline calculates
The pressure gradient that fluid flows in horizontal pipe is:
i = λ ρ V 2 2 D = 4 f ρ V 2 2 D - - - ( 4 )
In formula, i is pressure gradient, Pa/m; V is the mean flow rate of pipeline inner fluid, m/s; ρ is the density of fluid, kg/m 3D is round tube inside diameter, m; λ is the darcy coefficient of friction resistance; F is the peaceful coefficient of friction resistance of model.
The calculating of stopping transportation temperature drop
In the evaluation procedure of initial set probability, according to evaluation requirements, can take Lumped-Capacity method or two kinds of solutions of numerical method to the stopping transportation temperature drop, Lumped-Capacity method has computing velocity and the relatively conservative characteristics of result, and numerical solution calculates accurately, but computing velocity is slower.
(1) Lumped-Capacity method
After hot oil pipeline shutdown, apart from starting point any distance l place, in stopping transportation t managed after the time, the calculating formula of oil temperature is as follows:
T = T 0 + ( T R - T 0 ) exp ( - Kπ D out l Gc y - bt ) - - - ( 5 )
b = KD D 1 2 4 c y ρ y + ( D 2 2 - D 1 2 ) 4 c g ρ g - - - ( 6 )
In formula, T is apart from pipeline starting point l place, the oil temperature of stopping transportation t after the time, ℃; T 0Be ambient temperature (being the ground temperature at buried depth place to buried pipeline), ℃; T RBe the oil temperature at beginning pipeline starting point place during stopping transportation, ℃; Overall heat transfer coefficient when K is the pipeline steady-state operation, W/ (m ℃); D outBe pipe external diameter, m; G is the pipeline throughput rate, kg/s; L is that calculation level is apart from the distance of pipeline starting point, m; T is shutdown time, s; B is coefficient, s -1, value is determined by formula (2-30); D 1, D 2Be respectively pipe diameter and external diameter, m; c y, c gBe respectively the specific heat capacity of crude oil and steel pipe, kJ/ (kg ℃); ρ y, ρ gBe respectively specific heat capacity and the density of steel pipe, kg/m 3
(2) numerical solution
Temperature drop process after pipe inner wax crude oil shutdown is a three-dimensional unsteady heat transfer problem of following phase transformation, natural convection and moving boundary.According to the difference of heat transfer type, after stopping transportation, the heat transfer of the interior crude oil of pipe can be divided into three phases again: free convection heat transfer stage, natural convection and heat conduction co-controlling stage and pure heat conduction stage.
Based on the complicacy of temperature drop process after the waxy crude oil stopping transportation, introduce following hypothesis in this model and realize simplifying:
1) introducing equivalent heat conductivity processes the Natural Convection Problem of liquid crude oil as Heat Conduction Problems;
2) in the process of stopping transportation temperature drop, along with reducing the wax crystalline substance, oil temperature separates out gradually, form solidifying oil reservoir, due to the complicacy of engineering problem, suppose that growth and the pipeline of solidifying oil reservoir is concentric;
3) utilize the specific heat capacity variation with temperature to describe the wax partial crystallization and go out to discharge latent heat to the impact of stopping transportation temperature drop process;
4) in the subordinate phase of stopping transportation temperature drop, the concept of quoting stagnant point is divided into natural convection district and thermal conductive zone with pipe interior.
Be a pure Heat Conduction Problems by above supposition with stopping transportation temperature drop problem reduction.
In pipe, the heat transfer equation of crude oil is:
ρC ∂ T ∂ τ = 1 r ∂ ∂ r ( λ e r ∂ T ∂ r ) - - - ( 7 )
Wherein:
λ e = - α 1 ( T L - T w ) ( ∂ T ∂ r ) w - - - ( 8 )
During program was calculated, during higher than stagnant point, coefficient of heat conductivity adopted following formula to calculate when oil temperature, and oil temperature is chosen the coefficient of heat conductivity of crude oil during lower than stagnant point.
&lambda; = &lambda; e t > T zl &lambda; o t < T zl - - - ( 9 )
The heat transfer equation of wax deposition layer, tube wall and anticorrosive coat
&rho; n C n &PartialD; T n &PartialD; &tau; = 1 r &PartialD; &PartialD; r ( &lambda; n r &PartialD; T n &PartialD; r ) - - - ( 10 )
In formula, n=1,2,3 represent respectively wax deposition layer, tube wall and anticorrosive coat.
The thermal conduction of soil equation
&rho; s C s &PartialD; T s &PartialD; &tau; = 1 r &PartialD; &PartialD; r ( &lambda; s r &PartialD; T s &PartialD; r ) - - - ( 11 )
The periphery of soil Domain is First Boundary Condition: temperature is the natural ground temperature at buried depth of pipeline place.
Correlation Criteria is satisfied in the heat transfer of wax deposition layer, tube wall, anticorrosive coat and soil, that is:
T k | r = R k = T k + 1 | r = R k k=1,2,3
&lambda; 3 &PartialD; T 3 &PartialD; r | r = R 3 = &lambda; s &PartialD; T s &PartialD; r | r = R 3
The stopping transportation restart pressure calculates
Adopt Liu Tian to help model and calculate starting pressure, in model, starting pressure is divided into three parts: liquid is by the static inertial force that accelerates to toggle speed, the friction force between liquid and tube wall and destroy the thixotropic additional force of non-newtonian fluid structure.Suppose to start throughput rate and remain unchanged, the maximal value of the starting pressure of derivation is:
P max = 4 Q&rho;a &pi; d 2 + 4 &tau; &infin; L d + 4 a ( &tau; 0 - &tau; &infin; ) dB ( 1 - e - BL a s ) - - - ( 12 )
Tube wall shear rate and the relation that starts throughput rate Newtonian fluid stable-state flow formula approximate representation:
Figure BDA0000112386140000064
In formula: τ 0, τ , B: oil transportation temperature down cut speed is The time initial shearing stress, balance shear stress and thixotropy index, L is duct length, a sBe pressure-wave emission speed.
First, formula (12) the right is the inertia pressure drop, makes crude oil produce acceleration, by the static toggle speed that reaches; Second is the pressure drop of balance frictional resistance, the viscous friction drag when overcoming mobile the shearing between fluid and tube wall; The 3rd is the thixotropy pressure drop, is used for destroying near the structure of the crude oil of tube wall.
The present invention can whether solidifying pipe be quantitatively described to pipeline, and the suggestion of proposition reliability is optimized in security and the operation of pipeline.The present invention has considered the uncertainty of influence factor first in solidifying pipe probability assessment; Determined concrete acceptable failure probability index first in solidifying pipe probability assessment; Thixotropy of Crude Oil model, pressure wave transmission speed have been considered in the computation process of solidifying pipe probability.The present invention can instruct the crude oil pipeline safety evaluatio that flows, and provides technical support to safe, the efficient operation of pipeline.
Description of drawings
The solidifying pipe of Fig. 1 crude oil pipeline probability quantitative evaluation process flow diagram
Embodiment
Embodiment. the solidifying pipe of crude oil pipeline probability method for quantitatively evaluating by the method that this example proposes, has carried out the evaluation of solidifying pipe probability to domestic certain crude oil pipeline as shown in Figure 1, draws the conclusion that is adapted to produced on-site.
The statistical study of field data
In the evaluation of solidifying pipe probability, when thinking pipeline throughput rate, out-station temperature, steady-state operation, overall heat transfer coefficient, wax deposition thickness, buried depth of pipeline, the buried depth of pipeline place gentle soil thermal conductivity in ground etc. all have uncertainty, this pipeline data that Field Research obtains in February, 09 are carried out automatic statistical analysis, the results are shown in Table 2 and table 3.
The Safety Influence Factors uncertainty analysis of flowing of certain crude oil pipeline of table 2
Figure BDA0000112386140000071
The uncertain value of the mobile Safety Influence Factors of certain crude oil pipeline of table 3
Figure BDA0000112386140000072
Solidifying pipe probability calculation result and analysis
Solidifying pipe probability when (1) stopping transportation reaches the longest permission shutdown time
The longest permission shutdown time of target pipeline according to stipulating in the conduit running rules has carried out solidifying pipe probability calculation.The longest permission shutdown time of this pipeline in different months of rules regulation sees Table 4.
Different months of table 4 the longest permission shutdown time
Figure BDA0000112386140000081
Based on conduit running parameter during in May, 2008~2009 year April, calculate the solidifying pipe probability between this pipeline 6 stations, the results are shown in Table 5.
The solidifying pipe probability calculation result of pipeline when table 5 stopping transportation reaches the longest permission shutdown time
Figure BDA0000112386140000082
By table 5, in conjunction with estimating basic data as can be known, with respect to commingled crude (29 ℃ of condensation points), No. 2 former crude oil fluidities relatively poor (32 ℃ of condensation points) of carrying of through-station, carry the inlet temperature of initial station-No. 1 through-station of this crude oil, No. 1-No. 2 through-station pipeline sections generally at 34~35 ℃, higher 2~3 ℃ than condensation point of crude oil.Therefore-No. 1 through-station in initial station is in November, and No. 1-No. 2 through-station is bigger than normal at the solidifying pipe probability in February and November.And the inlet temperature between all the other stations of conveying commingled crude is higher 4~5 ℃ than condensation point of crude oil generally at 33~34 ℃, and therefore, solidifying pipe probability is lower, can suitably reduce inlet temperature again.
(2) the fixing solidifying pipe probability under shutdown time
Calculate the solidifying pipe probability under stopping transportation 24h condition, the results are shown in Table 6.
The solidifying pipe of pipeline probability under table 6 stopping transportation 24h
As shown in Table 6, at first stop-No. 1 through-station pipeline section is in November, No. 1-No. 2 through-station pipeline sections February and November solidifying pipe probability higher.
(3) the solidifying pipe probability under different shutdown times
Choose ground temperature and carry out computational analysis minimum February.Solidifying pipe probability calculation between the different shutdown time conditions of target pipeline lower 6 stations be the results are shown in Table 7.
Solidifying pipe probability calculation result under table 7 different shutdown time conditions in February
Figure BDA0000112386140000091
Can be found out by result in table 7, along with the increase of shutdown time, between each station solidifying pipe probability gather way and inconsistent, therefore the solidifying manage-style between each station is distinguished by inches to some extent.
(4) the solidifying pipe probability under different inlet temperatures
Calculate respectively February at 1 ℃, 2 ℃, the solidifying pipe probability of 24h after stopping transportation under 3 ℃ of conditions that enter the station and under the actual fed operating mode more than condensation point more than condensation point more than condensation point, see Table 8.
Solidifying pipe probability under the different inlet temperature conditions of table 8
Figure BDA0000112386140000092
As shown in Table 8, when between each station during all according to same index regulation inlet temperature, for example, when entering the station, the solidifying pipe probability at each station is also incomplete same according to " condensation point+3 ℃ ".To same pipeline (region is identical) under the same time, even result is like this, for the pipeline of carrying different crude oils, in Various Seasonal, this difference can be larger, therefore, searching determines that than " condensation point+3 ℃ " index that more reasonably enters the station or inlet temperature mobile Safety Assessment System is expected to address this problem comprehensively according to being possible.Result based on above-mentioned evaluation draws can instruct produced on-site and procedure making, can find out thus the significance of carrying out the quantitatively evaluating safety that flows.

Claims (8)

1. a crude oil pipeline coagulates pipe probability method for quantitatively evaluating, it is characterized in that its main flow process is:
1) robotization of safety effects factor is extracted and is analyzed;
2) the data distributivity rule of influence factor statistics;
3) determinacy stopping transportation restart pressure calculates;
4) estimate the setting of number of times and the grab sample of parameter;
5) calculating of uncertain stopping transportation restart pressure;
6) the statistics pipeline is restarted the calculation times of failure;
7) pipeline solidifying pipe probability calculation result draws;
8) whether satisfy and to accept the judgement of losing efficacy;
9) if do not satisfy, after adjusting the flow assurance scheme, turn 2) the data distributivity rule statistics of influence factor;
If satisfy, submit the flow assurance overall plan to.
2. crude oil pipeline according to claim 1 coagulates pipe probability method for quantitatively evaluating, it is characterized in that the solidifying pipe of described pipeline probability is:
For a given ultimate limit state, probability analysis comprises the model of broad sense random load S and the random drag R of broad sense; Corresponding limit state function can be expressed as following form in general manner:
g(x)=R-S (1)
Obviously, indicate inefficacy as g (x)<0;
The definition failure probability is as follows:
P f = &Integral; g ( x ) &le; 0 f x ( x ) dx - - - ( 2 )
In formula, x is the stochastic variable vector; f x(x) be joint probability density function;
Restart problem for the crude oil pipeline stopping transportation, definition load S is the minimum starting pressure of pipeline, and getting drag R is pipeline maximum allowable operating pressure (MAOP), and Pf is solidifying pipe probability.
3. crude oil pipeline according to claim 1 coagulates pipe probability method for quantitatively evaluating, it is characterized in that described data distributivity probability statistics and check are:
According to ISO 16708, the production run data of quantitative evaluation pipeline collection are carried out uncertainty and value analysis.
4. crude oil pipeline according to claim 1 coagulates pipe probability method for quantitatively evaluating, it is characterized in that described determinacy stopping transportation restart pressure calculating is:
This evaluation method is estimated framework based on secondary, first order framework in the situation that do not consider production run situation fluctuation, is to satisfy on-the-spot demand, with the security of the existing production and operation condition of Deterministic Methods assessment, if dangerous, will take the flow assurance measure according to result; The determinacy operating mode is on-the-spot given production data.
5. crude oil pipeline according to claim 1 coagulates pipe probability method for quantitatively evaluating, it is characterized in that the setting of described evaluation number of times and the grab sample of parameter are:
After having completed the deterministic parsing in system one-level framework, just can carry out large-scale probability sampling, the sampling number of times is defaulted as 10000 times, generates corresponding numerical value by its regularity of distribution function during the every sub-sampling of production run parameter, but must set the bound of production run parameter; In the situation that there is no the production scene appointment, the bound of system default is as follows, and throughput rate is between 3 times of variances of its average plus-minus, and temperature out of the station is between 0 ℃-90 ℃, and ground temperature is between-20 ℃-50 ℃.
6. crude oil pipeline according to claim 1 coagulates pipe probability method for quantitatively evaluating, it is characterized in that described uncertain stopping transportation restart pressure calculates and the calculating of the startup frequency of failure is:
Whether uncertain stopping transportation restart pressure calculates is to analyze solidifying pipe according to the production run parameter that the evaluation number of times of having set and each evaluation number of times generate to occur; As start unsuccessfully and to occur to be recorded as 1, do not occur to record 0; After completing, evaluation can record the number of times that total stopping transportation is restarted failure.
7. crude oil pipeline according to claim 1 coagulates pipe probability method for quantitatively evaluating, it is characterized in that described interpretation of result is:
Solidifying tube failure probability P f based on calculating is defined as pipeline shutdown with (1-Pf) and restarts safety indexes (SRSI); Under constant condition of pipeline shutdown time, the SRSI graded index is larger, and the possibility of pipeline shutdown failure is less restarting just, and the safe coefficient of pipeline restart procedure is just higher; Otherwise the possibility that pipeline is restarted failure is just larger, according to evaluation result, to the security of operation proposition safeguard of pipeline.
8. crude oil pipeline according to claim 1 coagulates pipe probability method for quantitatively evaluating, it is characterized in that the defeated technology Calculation of evaluation procedure pipe is:
Waterpower, heating power during normal the conveying calculate
(1) normally moving the axial temperature drop of hot oil pipeline calculates
Stablizing under transport condition, between the temperature drop rate that in hot oil pipeline, oil flows and the coefficient of heat conductivity of conduit running parameter, caliber, buried depth of pipeline and surrounding medium thereof, temperature, heating station, length and crude oil property are closely related; In pipe, oil stream can be similar to Л e й б e н з о Н formula along the journey Temperature Distribution and describe:
T L=(T 0+b)+(T R-T 0-b)e -aL (3)
In formula, a=K π D/Gc, b=ig/c pa; G is the mass rate of oil product, kg/s; C pBe the specific heat capacity of crude oil, J/ (kg ℃); D is outer diameter tube, m; L is the distance apart from heating station, m; K is the pipeline overall heat transfer coefficient, W/ (m 2℃); T RBe heating station departures oil temperature, ℃; T LBe the oil temperature of distance heating station distance for the L place, ℃; T0 is ambient temperature, and buried pipeline is got the natural ground temperature in tube hub buried depth place, ℃;
(2) frictional resistance of hot oil pipeline calculates
The pressure gradient that fluid flows in horizontal pipe is:
i = &lambda; &rho; V 2 2 D = 4 f &rho; V 2 2 D - - - ( 4 )
In formula, i is pressure gradient, Pa/m; V is the mean flow rate of pipeline inner fluid, m/s; ρ is the density of fluid, kg/m 3D is round tube inside diameter, m; λ is the darcy coefficient of friction resistance; F is the peaceful coefficient of friction resistance of model;
The calculating of stopping transportation temperature drop
In the evaluation procedure of initial set probability, according to evaluation requirements, can take Lumped-Capacity method or two kinds of solutions of numerical method to the stopping transportation temperature drop, Lumped-Capacity method has computing velocity and the relatively conservative characteristics of result, and numerical solution calculates accurately, but computing velocity is slower;
(1) Lumped-Capacity method
After hot oil pipeline shutdown, apart from starting point any distance l place, in stopping transportation t managed after the time, the calculating formula of oil temperature is as follows:
T = T 0 + ( T R - T 0 ) exp ( - K&pi; D out l Gc y - bt ) - - - ( 5 )
b = KD D 1 2 4 c y &rho; y + ( D 2 2 - D 1 2 ) 4 c g &rho; g - - - ( 6 )
In formula, T is apart from pipeline starting point l place, the oil temperature of stopping transportation t after the time, ℃; T 0Be ambient temperature, ℃; T RBe the oil temperature at beginning pipeline starting point place during stopping transportation, ℃; Overall heat transfer coefficient when K is the pipeline steady-state operation, W/ (m ℃); D outBe pipe external diameter, m; G is the pipeline throughput rate, kg/s; L is that calculation level is apart from the distance of pipeline starting point, m; T is shutdown time, s; B is coefficient, s -1, value is determined by formula (2-30); D 1, D 2Be respectively pipe diameter and external diameter, m; c y, c gBe respectively the specific heat capacity of crude oil and steel pipe, kJ/ (kg ℃); ρ y, ρ gBe respectively specific heat capacity and the density of steel pipe, kg/m 3
(2) numerical solution
Temperature drop process after pipe inner wax crude oil shutdown is a three-dimensional unsteady heat transfer problem of following phase transformation, natural convection and moving boundary; According to the difference of heat transfer type, after stopping transportation, the heat transfer of the interior crude oil of pipe can be divided into three phases again: free convection heat transfer stage, natural convection and heat conduction co-controlling stage and pure heat conduction stage;
Based on the complicacy of temperature drop process after the waxy crude oil stopping transportation, introduce following hypothesis in this model and realize simplifying:
1) introducing equivalent heat conductivity processes the Natural Convection Problem of liquid crude oil as Heat Conduction Problems;
2) in the process of stopping transportation temperature drop, along with reducing the wax crystalline substance, oil temperature separates out gradually, form solidifying oil reservoir, due to the complicacy of engineering problem, suppose that growth and the pipeline of solidifying oil reservoir is concentric;
3) utilize the specific heat capacity variation with temperature to describe the wax partial crystallization and go out to discharge latent heat to the impact of stopping transportation temperature drop process;
4) in the subordinate phase of stopping transportation temperature drop, the concept of quoting stagnant point is divided into natural convection district and thermal conductive zone with pipe interior;
Be a pure Heat Conduction Problems by above supposition with stopping transportation temperature drop problem reduction;
In pipe, the heat transfer equation of crude oil is:
&rho;C &PartialD; T &PartialD; &tau; = 1 r &PartialD; &PartialD; r ( &lambda; e r &PartialD; T &PartialD; r ) - - - ( 7 )
Wherein:
&lambda; e = - &alpha; 1 ( T L - T w ) ( &PartialD; T &PartialD; r ) w - - - ( 8 )
During program was calculated, during higher than stagnant point, coefficient of heat conductivity adopted following formula to calculate when oil temperature, and oil temperature is chosen the coefficient of heat conductivity of crude oil during lower than stagnant point;
&lambda; = &lambda; e t > T zl &lambda; o t < T zl - - - ( 9 )
The heat transfer equation of wax deposition layer, tube wall and anticorrosive coat
&rho; n C n &PartialD; T n &PartialD; &tau; = 1 r &PartialD; &PartialD; r ( &lambda; n r &PartialD; T n &PartialD; r ) - - - ( 10 )
In formula, n=1,2,3 represent respectively wax deposition layer, tube wall and anticorrosive coat;
The thermal conduction of soil equation
&rho; s C s &PartialD; T s &PartialD; &tau; = 1 r &PartialD; &PartialD; r ( &lambda; s r &PartialD; T s &PartialD; r ) - - - ( 11 )
The periphery of soil Domain is First Boundary Condition: temperature is the natural ground temperature at buried depth of pipeline place; Correlation Criteria is satisfied in the heat transfer of wax deposition layer, tube wall, anticorrosive coat and soil, that is:
T k | r = R k = T k + 1 | r = R k k=1,2,3
&lambda; 3 &PartialD; T 3 &PartialD; r | r = R 3 = &lambda; s &PartialD; T s &PartialD; r | r = R 3
The stopping transportation restart pressure calculates
Adopt Liu Tian to help model and calculate starting pressure, in model, starting pressure is divided into three parts: liquid is by the static inertial force that accelerates to toggle speed, the friction force between liquid and tube wall and destroy the thixotropic additional force of non-newtonian fluid structure; Suppose to start throughput rate and remain unchanged, the maximal value of the starting pressure of derivation is:
P max = 4 Q&rho;a &pi; d 2 + 4 &tau; &infin; L d + 4 a ( &tau; 0 - &tau; &infin; ) dB ( 1 - e - BL a s ) - - - ( 12 )
Tube wall shear rate and the relation that starts throughput rate Newtonian fluid stable-state flow formula approximate representation:
Figure FDA0000112386130000054
In formula: τ 0, τ , B: oil transportation temperature down cut speed is
Figure FDA0000112386130000055
The time initial shearing stress, balance shear stress and thixotropy index, L is duct length, a sBe pressure-wave emission speed;
First, formula (12) the right is the inertia pressure drop, makes crude oil produce acceleration, by the static toggle speed that reaches; Second is the pressure drop of balance frictional resistance, the viscous friction drag when overcoming mobile the shearing between fluid and tube wall; The 3rd is the thixotropy pressure drop, is used for destroying near the structure of the crude oil of tube wall.
CN201110380216.1A 2011-11-25 2011-11-25 Pipe condensation probability of crude oil pipeline method for quantitatively evaluating Active CN103136701B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110380216.1A CN103136701B (en) 2011-11-25 2011-11-25 Pipe condensation probability of crude oil pipeline method for quantitatively evaluating

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110380216.1A CN103136701B (en) 2011-11-25 2011-11-25 Pipe condensation probability of crude oil pipeline method for quantitatively evaluating

Publications (2)

Publication Number Publication Date
CN103136701A true CN103136701A (en) 2013-06-05
CN103136701B CN103136701B (en) 2017-07-14

Family

ID=48496500

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110380216.1A Active CN103136701B (en) 2011-11-25 2011-11-25 Pipe condensation probability of crude oil pipeline method for quantitatively evaluating

Country Status (1)

Country Link
CN (1) CN103136701B (en)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103953854A (en) * 2014-04-02 2014-07-30 中国石油大学(北京) Method and device for determining running parameters of waxy crude oil pipeline
CN104123425A (en) * 2014-08-07 2014-10-29 中国石油大学(北京) Crude oil pipeline preheating production thermal process simulation method
CN104318479A (en) * 2014-10-17 2015-01-28 中国石油天然气股份有限公司 Method for determining critical state and safety allowance when crude oil pipeline runs
CN105042627A (en) * 2015-07-24 2015-11-11 中国石油化工股份有限公司 Closed-loop regulation method and system for gas consumption of heating furnace of transfer station
CN105404776A (en) * 2015-11-13 2016-03-16 中国石油天然气股份有限公司 Reliability evaluating method for pipe with corrosion defect
WO2017024699A1 (en) * 2015-08-11 2017-02-16 深圳朝伟达科技有限公司 Safety margin determination method
CN104199804B (en) * 2014-08-15 2017-04-05 中国石油天然气股份有限公司 High water cut soil buried heated pipeline flowing failure mode decision method
CN108119761A (en) * 2017-12-14 2018-06-05 中国石油天然气股份有限公司 A kind of gas-liquid two-phase flows down natural gas line extent of corrosion Forecasting Methodology
CN109723977A (en) * 2018-12-27 2019-05-07 中国石油大学(北京) A kind of stopping transportation leak judgement method, equipment and the system of processed oil pipeline
CN111542678A (en) * 2018-01-05 2020-08-14 罗克斯软件解决方案有限公司 Well flow simulation system
CN111859295A (en) * 2020-07-15 2020-10-30 中国安全生产科学研究院 Quantitative risk evaluation method for oil pipeline
CN112749469A (en) * 2019-10-31 2021-05-04 中国石油天然气股份有限公司 Crude oil pipeline monitoring method and device and computer readable storage medium

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6527015B2 (en) * 1999-07-02 2003-03-04 F. Glenn Lively Insulated pipe
CN101256026A (en) * 2007-03-01 2008-09-03 中国石油天然气股份有限公司 Primary oil long-distance pipeline electric heating membrane safety heating system and construction method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6527015B2 (en) * 1999-07-02 2003-03-04 F. Glenn Lively Insulated pipe
CN101256026A (en) * 2007-03-01 2008-09-03 中国石油天然气股份有限公司 Primary oil long-distance pipeline electric heating membrane safety heating system and construction method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
帅健: "油气管道可靠性的极限状态设计方法", 《石油规划设计》 *
苗青等: "秦京输油管道停输再启动安全性评价与建议", 《油气储运》 *

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103953854A (en) * 2014-04-02 2014-07-30 中国石油大学(北京) Method and device for determining running parameters of waxy crude oil pipeline
CN104123425B (en) * 2014-08-07 2017-12-08 中国石油大学(北京) Crude oil pipeline preheats operation thermal procession emulation mode
CN104123425A (en) * 2014-08-07 2014-10-29 中国石油大学(北京) Crude oil pipeline preheating production thermal process simulation method
CN104199804B (en) * 2014-08-15 2017-04-05 中国石油天然气股份有限公司 High water cut soil buried heated pipeline flowing failure mode decision method
CN104318479A (en) * 2014-10-17 2015-01-28 中国石油天然气股份有限公司 Method for determining critical state and safety allowance when crude oil pipeline runs
CN105042627A (en) * 2015-07-24 2015-11-11 中国石油化工股份有限公司 Closed-loop regulation method and system for gas consumption of heating furnace of transfer station
WO2017024699A1 (en) * 2015-08-11 2017-02-16 深圳朝伟达科技有限公司 Safety margin determination method
WO2017024534A1 (en) * 2015-08-11 2017-02-16 深圳朝伟达科技有限公司 Method for determining critical state during operation of crude oil pipeline
CN105404776A (en) * 2015-11-13 2016-03-16 中国石油天然气股份有限公司 Reliability evaluating method for pipe with corrosion defect
CN105404776B (en) * 2015-11-13 2018-07-13 中国石油天然气股份有限公司 A kind of reliability evaluation method of the pipeline containing corrosion default
CN108119761A (en) * 2017-12-14 2018-06-05 中国石油天然气股份有限公司 A kind of gas-liquid two-phase flows down natural gas line extent of corrosion Forecasting Methodology
CN108119761B (en) * 2017-12-14 2020-02-14 中国石油天然气股份有限公司 Method for predicting corrosion degree of natural gas pipeline under gas-liquid two-phase flow
CN111542678A (en) * 2018-01-05 2020-08-14 罗克斯软件解决方案有限公司 Well flow simulation system
US11809793B2 (en) 2018-01-05 2023-11-07 Roxar Software Solutions As Well flow simulation system
CN109723977A (en) * 2018-12-27 2019-05-07 中国石油大学(北京) A kind of stopping transportation leak judgement method, equipment and the system of processed oil pipeline
CN112749469A (en) * 2019-10-31 2021-05-04 中国石油天然气股份有限公司 Crude oil pipeline monitoring method and device and computer readable storage medium
CN112749469B (en) * 2019-10-31 2022-11-04 中国石油天然气股份有限公司 Crude oil pipeline monitoring method and device and computer readable storage medium
CN111859295A (en) * 2020-07-15 2020-10-30 中国安全生产科学研究院 Quantitative risk evaluation method for oil pipeline
CN111859295B (en) * 2020-07-15 2023-06-02 中国安全生产科学研究院 Quantitative risk evaluation method for oil pipeline

Also Published As

Publication number Publication date
CN103136701B (en) 2017-07-14

Similar Documents

Publication Publication Date Title
CN103136701A (en) Method for quantitatively evaluating pipe condensation probability of crude oil pipeline
Li et al. Advances and future challenges of wax removal in pipeline pigging operations on crude oil transportation systems
US20230204162A1 (en) Automated Re-Melt Control Systems
Rittirong et al. Experimental study of paraffin deposition under two-phase gas/oil slug flow in horizontal pipes
Kjolaas et al. Experiments for low liquid loading with liquid holdup discontinuities in two-and three-phase flows
CN102778418A (en) Crude oil pipeline initial condensation probability quantitative evaluation method
Abdul-Majeed et al. Prediction of slug frequency for medium liquid viscosity two-phase flow in vertical, horizontal, and inclined pipes
CN107859880A (en) A kind of flow investigation analysis method that behavior is stolen for metering process after the pressure regulation of combustion gas elder generation
Njoku et al. Comparison of energy, exergy and entropy generation-based criteria for evaluating stratified thermal store performances
CN103741692B (en) A kind of intelligent water communication method of mass concrete and the system of use the method
Kutyłowska Application of MARSplines method for failure rate prediction
Zhang et al. Numerical investigation on the longitudinal distribution characteristics of falling film flow outside a vertical tube
JP6201640B2 (en) Steam pipe loss measurement system and measurement method
Jiang et al. The effects of oil property and inclination angle on oil–water core annular flow through U-Bends
Hilgefort Big data analysis using Bayesian network modeling: a case study with WG-ICDA of a gas storage field
Liu Variable water flow pumping for central chilled water systems
US20180334621A1 (en) Crude hydrocarbon fluids demulsification system
CN104102852A (en) Computing method of preheating commissioning thermodynamic process of crude oil pipeline
Liu et al. Study on the segmental pigging process for the wet gas gathering and transportation pipelines of the Puguang gas field
Zhao et al. Application of internal corrosion direct assessment in CO2 slug flow submarine pipelines
Hamad et al. Development of kerosene–water two‐phase up‐flow in a vertical pipe downstream of A 90° bend
Yap et al. Investigation of Two-Phase Liquid-Droplet Flow With Particle Deposition in the Heat Exchanger
Ryu Hydrodynamic Cause Analysis for Roof Damage on Water Storage Tank
Zhang et al. A new approach to calculate shutdown and restart pressure of waxy crude pipeline
Jiang et al. Robust pressure prediction of oil and gas pipeline networks based on equipment embedding neural network

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20211105

Address after: Room 08-10, 6 / F, block a, No. 5, Dongtucheng Road, Chaoyang District, Beijing 100013

Patentee after: National Petroleum and natural gas pipeline network Group Co.,Ltd.

Address before: 100007 Oil Mansion, 9 North Avenue, Dongcheng District, Beijing, Dongzhimen

Patentee before: PETROCHINA Co.,Ltd.