CN105426692A - Ocean platform multi-stage task system reliability estimation method based on data drive - Google Patents
Ocean platform multi-stage task system reliability estimation method based on data drive Download PDFInfo
- Publication number
- CN105426692A CN105426692A CN201510925359.4A CN201510925359A CN105426692A CN 105426692 A CN105426692 A CN 105426692A CN 201510925359 A CN201510925359 A CN 201510925359A CN 105426692 A CN105426692 A CN 105426692A
- Authority
- CN
- China
- Prior art keywords
- phi
- sigma
- time
- phased mission
- stage
- 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
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Testing Or Calibration Of Command Recording Devices (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
The invention relates to an ocean platform multi-stage task system reliability estimation method based on data drive. The estimation method mainly includes the steps that a multi-stage task process model is established, and a probability density function of task time is estimated; a system degenerative process model is established, and the residual life is estimated; according to the reliability definition of the multi-stage task process model and the system degenerative model, system reliability parameter values are obtained through real-time calculation through the task time estimated value and the system service life estimated value. Through measured sensor data, the relation between historical data of an individual multi-stage task system and current real-time data information is established, reliability parameters of the multi-stage task system under the dynamic condition are accurately estimated in real time, and the problem that an existing method is only suitable for static hypothesis and a specific system is solved. The reliability parameters of the ocean platform multi-stage task system are accurately estimated in real time, the system is effectively maintained in real time, the maintenance cost of the system is reduced, and catastrophic faults are effectively avoided.
Description
Technical field
The invention belongs to reliability engineering technique field, specifically, relate to a kind of ocean platform phased mission systems reliability estimation methods.
Background technology
Ocean platform is the complex apparatus system for completing offshore and gas development work, usually, this kind of complication system of ocean platform needs through multiple stage when completing particular task, and the degradation characteristics in each stage is different, in existing literature research achievement, this type systematic is called as phased mission systems.The basic characteristics of phased mission systems are: continuous, each stage non-overlapping in multiple, stage.Within different consecutive hours intervals (i.e. stage), the phased mission that system completes is different, so the structural reliability in system each stage is also different.For ocean platform system, it completes offshore and gas development work need experience several stages such as offshore platform is built, drilling well, pressure break, running casing, oil recovery, only has the Mission Success in all stages reliably to complete finally completing of guarantee system task.
In complication system, reliability is one of of paramount importance index parameter in system, operation and maintaining process, be one for weighing the Important Parameters of engineering system performance.Therefore, the accurate estimation of ocean platform phased mission systems reliability for the maintenance of system and maintenance very important, system maintenance cost can be reduced to a certain extent, effectively the generation of averting a calamity property fault.
For the reliability estimation problem of this kind of phased mission systems of ocean platform, main employing Fault Tree Analysis and the method based on binary decision diagrams (bdds) in existing achievement in research, but the structure that these methods too rely on phased mission systems is known, but in fact, in engineer applied phased mission systems usually due to too complicated and be difficult to the complete structural information of Obtaining Accurate.For ocean platform, ocean platform at sea works, the interference such as random wave, tide, ocean current, floating ice, earthquake, corrosion is many, the non-intellectual in down-hole is strong, working environment is extremely complicated and uncertain, its mechanism dynamic model is difficult to accurately set up in advance, which dictates that its complete structural information can not Obtaining Accurate.Thus, above-mentioned two class main method are very limited in the application, thus are difficult to the reliability effectively accurately estimating system.On the other hand, under existing method is only applicable to static assumed condition, be a kind of off-line method in essence, contacting between reliability and history observation data and current Real-time Monitoring Data cannot be set up, be not useable for the dynamic estimation of system reliability.So, existing method is adopted to be difficult to estimate exactly the reliability of phased mission systems, thus cannot the fault of prognoses system effectively, and any little potential faults all may develop into large accident, cause the accidents such as the blowout of ocean platform, fire, blast, leakage of oil to occur, cause platform to damage, tilt or topple.Therefore, solution is estimated in the urgent need to seeking a kind of new phased mission systems reliability prediction.
Along with the development of infotech, ocean platform great majority are all equipped with data acquisition and transmission system, can physical quantity on Real-time Collection ocean platform, as voltage, electric current, crude oil temperature, the data such as polished rod load and displacement of motor.But these data only for simply detecting and form, lacking the data mining that manufacturing-oriented security monitoring and reliability prediction etc. are profound, thus causing the waste of resource at present.For this reason, these a large amount of sensor data information are that the reliability prediction of ocean platform multistage negotiation task system is estimated to provide a potential solution route.
Summary of the invention
The object of the invention is to the above-mentioned defect and the deficiency that overcome prior art existence, provide a kind of ocean platform phased mission systems reliability estimation methods based on data-driven, the method adopts the sensing data observed, build the relation between historical data information and current real-time information, accurately estimate in real time phased mission systems reliability parameter under dynamic condition, solve the problem that existing method is only applicable to static assumed condition.
According to one embodiment of the invention, provide a kind of ocean platform phased mission systems reliability estimation methods based on data-driven, the steps include:
Detected by the environment of sensor to system, the system environments Monitoring Data that comprehensive utilization sensor is observed, set up contacting between the historical data of individual phased mission systems and current real-time data information, stochastic filtering model is adopted to carry out modeling to stage interval, obtain phased mission system process model, estimate the probability density function of task time, obtain estimated value task time of phased mission systems.
Fusion application historical data information, sets up phased mission systems degenerative process model, estimates the life-span of phased mission systems, obtains the life estimation value of phased mission systems.
According to the reliability definition of phased mission system process model and phased mission systems degenerative process model, comprehensive utilization estimated value task time and life estimation value calculate the reliability reference value of phased mission systems in real time.
According in the method for estimation of the embodiment of the present invention, set up phased mission system process model and estimate that the concrete steps of probability density function of task time are: suppose that ocean platform offshore and gas development task has N number of stage, make stochastic variable X
nrepresent the n-th phase duration, its span is
make stochastic variable T
mrepresent the T.T. of completion system overall tasks, and have
order
for each stage X
n, the probability density function of n=1, L, N, wherein 2≤n≤N; Make Φ
i,nrepresent cut-off current time t
ithe environmental monitoring data information aggregate corresponding with the n-th stage; Adopt probability density function task time in environmental monitoring data estimation of the order 1 ~ N each stage.
According in the method for estimation of the embodiment of the present invention, the n-th phase process modeling and task time PDF estimation step be: suppose to get current time t by sensor
ithe environmental monitoring data information φ in lower n-th stage
i,n, make L
i,nrepresent t
ithe residue continuous time in stage in moment n-th, then through type L
i+1, n=L
i,n-(t
i+1-t
i), (if L
i,n> t
i+1-t
i) residue L continuous time of the n-th stage subsequent time can be obtained
i+1, n, L
0, nrepresent the residue continuous time of first stage; The consecutive hours area of a room considering task time always on the occasion of, utilize transform Z
i,n=lnL
i,n, obtain a process variable Z
i,n, to ensure L
i,n> 0, adopts a sliding scales parameter to set up φ
i,nand L
i,nstochastic relation model between the two, thus by φ
i,nestimate L
i,n, make φ
i,nbe expressed as Z
i,nfunction, set up the n-th phase process model as follows:
Z
i,n=lnL
i,n(1)
φ
i,n=g
n(Z
i,n)+η
i,n
Wherein, g
n(Z
i,n) be a function undetermined, it is described that the relation between task process and the environmental monitoring data corresponding with the stage; η
i,nthe measuring error of a Normal Distribution, and
Adopt extended Kalman filter to estimate or upgrade Z
i,nconditional probability density function, and ask for further residue continuous time L
i,n; Definition Z
i,nand Z
i+1, nrenewal amount and one-step prediction conditional probability density function amount be respectively Z
i,n| Φ
i,n: N (z
i|i, n, P
i|i, n) and Z
i+1, n| Φ
i,n: N (z
i+1|i, n, P
i+1|i, n), parameter z here
i|i, n, P
i|i, n, z
i+1|i, nand P
i+1|i, ncan be calculated by extended Kalman filter algorithm, its computing formula is as follows:
z
i|i,n=z
i|i-1,n+K
i,n[φ
i,n-g
n(z
i|i-1,n)](2)
Wherein, K
i,nrepresent the gain battle array of Kalman filtering, obtained by following formula:
Wherein
Correspondingly, the renewal equation of state estimation variance matrix can be obtained by following formula:
P
i|i,n=P
i|i-1,n(1-K
i,ng′
n(z
i|i-1,n))(4)
The initial parameter z of extended Kalman filter algorithm
0|0, nand P
0|0, nestimated to obtain by historical data.
Calculation expectation z is needed in the renewal equation of above-mentioned state estimation variance matrix
i|i-1, nwith variance P
i|i-1, na step estimated value, consider Z
i,n| Φ
i,n: N (z
i|i, n, P
i|i, n), and Z
i,n=lnL
i,n, obtain these two amounts by lower two formulas:
Consider the transformational relation of normal distribution and lognormal distribution, have L
i, 1| Φ
i, 1: lnN (z
i|i, 1, P
i|i, 1), then based on model expression, obtain and expect z
i|i-1,1with amount of variation P
i|i-1,1expression formula:
P
i|i-1,n=P
i-1|i-1,n(8)
By above-mentioned renewal expression formula, obtain the estimated probability density function that the first stage remains continuous time
for
Again because of
Then directly obtain X by correspondent transform
nat t
ithe distribution in moment
Wherein,
calculated by formula (9).
Thus, corresponding to t
imoment and X
1subordinate phase duration probability density function expression formula under condition is
Similarly, corresponding to t
imoment and X
npostorder under condition each phase duration probability density function expression formula is
Wherein s=n+1 ..., N.
Further, have
then task time T
mdistribution function be
At moment t
mprobability density function task time formula (14) differential process being obtained to the n-th stage is
So far, the environmental monitoring data information of being correlated with by the current generation obtains estimator task time in this stage.
According in the method for estimation of the embodiment of the present invention, Wiener-Hopf equation is used to set up phased mission systems degenerative process model, the concrete steps setting up phased mission systems degenerative process model are: assuming that the initial reading of degenerative process is Y (0)=0, if initial time is t
i(i=1,2 ..., N), then the monitored parameters evolution process model of time correlation is
Y(t)=y
i+λ(t-t
i)+σB(t-t
i)(16)
Further, consider fusion application historical data information, in model, λ develops into dynamic parameter, uses λ
i=λ
i-1+ η replaces, and wherein η ~ N (0, Q), then degradation model is reconstructed into
Wherein ε
i~ N (0, t
i-t
i-1), λ
iestimated value can be obtained by Kalman filtering algorithm.
λ is obtained by Gaussian distribution hypothesis and Bayesian filter principle
iprobability density function be
The life-span of equipment is defined by the concept of first-hitting time, when the degradation values Y (t) that formula (16) is determined reaches the failure threshold w preset first, just thinks equipment failure; By the time that the timing definition that equipment life stops is degenerative process Y (t) First excursion failure threshold w, then equipment t
ithe remaining life S in moment
ibe defined as:
S
i=inf{s
i:Y(s
i+t
i)≥w|Y
i}(19)
Order
for the cumulative distribution function of remaining life, predict t by degradation model formula (17)
ithe degeneration probability density function in moment
consideration model profile is normal distribution, has Y (s
i+ t
i) | Y
i:
then t
iprobability density function and the cumulative distribution function of the life expectancy of moment phased mission systems are respectively:
According in the method for estimation of the embodiment of the present invention, obtain λ by Kalman filtering algorithm
ithe step of estimated value be:
Initialization
p
0;
T
ithe state estimation in moment is
P
i|i-1=P
i-1|i-1+Q
K
i=(t
i-t
i-1)
2P
i|i-1+σ
2(t
i-t
i-1);
Upgrade
According in the method for estimation of the embodiment of the present invention, obtain the life-span probability density function of task system
with probability density function task time
after, according to the reliability definition of phased mission system process model and phased mission systems degenerative process model, estimate t
ithe reliability of moment phased mission systems, the probability obtaining the reliability parameter of the n-th stage fill is
The estimation of the reliability parameter value of ocean platform complex phased mission systems can be completed based on formula (21) and (22).
The ocean platform phased mission systems reliability estimation methods based on data-driven that the present invention proposes, adopt observe sensing data, by setting up contacting between the historical data of individual phased mission systems and current real-time data information, accurately estimate in real time the reliability parameter of phased mission systems under dynamic condition, not only solve the reliability estimation problem of individual phased mission systems, the reliability also completing system under dynamic condition is estimated.Existing method cannot utilize real-time information, be only applicable to static assumed condition and particular system, compared with the conventional method, by the ocean platform phased mission systems reliability estimation methods based on data-driven according to the embodiment of the present invention, accurately can estimate the reliability parameter of this kind of phased mission systems of dynamic condition Offshore Platform Under, solve the problem that existing method is only applicable to static assumed condition, thus can effectively safeguard system in time, for the maintenance of system and the determination of repair part ordering strategy provide theoretical foundation and the technical support of power, reduce the maintenance cost of system, the effectively generation of averting a calamity property fault, engineer applied is worth high.
Embodiment
Embodiments of the present invention is further illustrated below.
In order to without loss of generality, for ocean platform three stage fill, according to a kind of ocean platform reliability estimation methods based on data-driven of the embodiment of the present invention, the steps include:
Step one: detected by the environment of sensor to system, the system environments Monitoring Data that comprehensive utilization sensor is observed, set up contacting between the historical data of individual three stage fill and current real-time data information, stochastic filtering model is adopted to carry out modeling to stage interval, obtain three phased mission process models, estimate the probability density function of task time, obtain estimated value task time of three stage fill further.
Suppose that ocean platform offshore and gas development task has 3 stages, make stochastic variable X
nrepresent (n=1,2, the 3) duration in the n-th stage, its span is
make stochastic variable T
mrepresent the T.T. of completion system overall tasks, and have
order
for each phase duration X
n, n=1, the probability density function of 2,3; Make Φ
i,nrepresent cut-off current time t
ithe environmental monitoring data information aggregate corresponding with the n-th stage, adopts environmental monitoring data to estimate probability density function task time in each stage.
(1) first stage process model building and task time PDF estimation
Suppose to get current time t by sensor
ithe environmental monitoring data information φ of lower first stage
i, 1, make L
i, 1represent t
itime inscribe residue continuous time of first stage, then through type L
i+1,1=L
i, 1-(t
i+1-t
i), (if L
i+1,1> t
i+1-t
i) residue L continuous time of first stage subsequent time can be obtained
i+1,1.Further, the consecutive hours area of a room considering task time always on the occasion of, thus, transform Z can be utilized
i, 1=lnL
i, 1, obtain a process variable Z
i, 1, to ensure L
i, 1> 0, adopts a sliding scales parameter to set up φ
i, 1and L
i, 1stochastic relation model between the two, thus by φ
i, 1estimate L
i, 1, make φ
i, 1be expressed as Z
i, 1function, set up first stage process model as follows:
Z
i,1=lnL
i,1(23)
φ
i,1=g
1(Z
i,1)+η
i,1
Wherein, g
1(Z
i, 1) be a function undetermined, it is described that the relation between task process and the environmental monitoring data corresponding with the stage; η
i, 1the measuring error of a Normal Distribution, and
Extended Kalman filter (hereinafter referred to as EKF) is adopted to estimate/upgrade Z
i, 1conditional probability density function, and ask for further residue continuous time L
i, 1.Definition Z
i, 1renewal amount and one-step prediction conditional probability density function amount be respectively Z
i, 1| Φ
i, 1: N (z
i|i, 1, P
i|i, 1) and Z
i+1,1| Φ
i, 1: N (z
i+1|i, 1, P
i+1|i, 1), parameter z here
i|i, 1, P
i|i, 1, z
i+1|i, 1and P
i+1|i, 1can be calculated by EKF algorithm, its computing formula is as follows:
z
i|i,1=z
i|i-1,1+K
i,1[φ
i,1-g
1(z
i|i-1,1)](24)
Wherein, K
i, 1represent the gain battle array of Kalman filtering, obtained by following formula
Wherein
Correspondingly, the renewal equation of state estimation variance matrix can be obtained by following formula:
P
i|i,1=P
i|i-1,1(1-K
i,1g′
1(z
i|i-1,1))(26)
The initial parameter z of EKF algorithm
0|0,1and P
0|0,1estimated to obtain by historical data, in the renewal equation of above-mentioned state estimation variance matrix, need calculation expectation z
i|i-1,1with variance P
i|i-1,1a step estimated value, consider Z
i, 1| Φ
i, 1: N (z
i|i, 1, P
i|i, 1), and Z
i, 1=lnL
i, 1, these two amounts can be obtained by lower two formulas:
Consider the transformational relation of normal distribution and lognormal distribution, have L
i, 1| Φ
i, 1: lnN (z
i|i, 1, P
i|i, 1), then based on model expression, obtain and expect z
i|i-1,1with variance P
i|i-1,1expression formula:
P
i|i-1,1=P
i-1|i-1,1(30)
By above-mentioned renewal expression formula, obtain the estimated probability density function that the first stage remains continuous time
for
Again because of
X
1|Φ
i,1=L
i,1+t
i(32)
Then by X that correspondent transform directly obtains
1probability density function
Here
calculated by formula (31).
Corresponding to t
imoment and X
1the probability density function expression formula of the subordinate phase duration under condition is
Similarly, corresponding to t
imoment and X
1duration phase III probability density function expression formula under condition is
Further, if make T
m| Φ
i, 1=(X
1+ X
2+ X
3) | Φ
i, 1, then task time T
mdistribution function be
At moment t
mthe probability density function task time process of above formula differential being obtained to the first stage is
(2) subordinate phase process model building and task time PDF estimation
Subordinate phase process model is as follows:
Z
i,2=lnL
i,2(38)
φ
i,2=g
2(Z
i,2)+η
i,2
Wherein, g
2(Z
i, 2) be a function undetermined; η
i, 2the measuring error of a Normal Distribution, and
Obtain t
ithe environmental monitoring data information φ in moment
i, 2after, based on Monitoring Data set Φ
i, 2, adopt EKF estimation/update condition L
i, 2probability density function, make Z
i, 2| Φ
i, 2: N (z
i|i, 2, P
i|i, 2), Z
i+1,2| Φ
i, 2: N (z
i+1|i, 2, P
i+1|i, 2), then Z
i, 2have
z
i|i,2=z
i|i-1,2+K
i,2[φ
i,2-g
2(z
i|i-1,2)](39)
P
i|i,2=P
i|i-1,2(1-K
i,2g′
2(z
i|i-1,2))(41)
Wherein
Further, z is expected
i|i-1,2with amount of variation P
i|i-1,2expression formula:
P
i|i-1,2=P
i-1|i-1,2(43)
By above-mentioned renewal expression formula, obtain the estimated probability density function of subordinate phase residue continuous time
for
Based on X
2with L
i, 2relation, have
X
2|Φ
i,2=L
i,2+(t
i-x
1)(45)
Wherein, x
1obtained by first stage task process according to model specification.Thus, by variable transitions, directly X can be write out
2distribution and expression formula
Wherein
can be obtained by formula (44).
Corresponding to t
itime the duration phase III probability density function expression formula of inscribing be
Then task time T
mdistribution function be
Further, at moment t
mthe probability density function task time process of above formula differential being obtained to subordinate phase is
(3) phase III process model building and task time PDF estimation
Phase III process model is as follows:
Z
i,3=lnL
i,3(50)
φ
i,3=g
3(Z
i,3)+η
i,3
Wherein, g
3(Z
i, 3) be a function undetermined; η
i, 3the measuring error of a Normal Distribution, and
Obtain t
ithe environmental monitoring data information φ in moment
i, 3after, based on Monitoring Data set Φ
i, 3, adopt EKF estimation/update condition L
i, 3probability density function, then Z
i, 3have
z
i|i,3=z
i|i-1,3+K
i,3[φ
i,3-g
3(z
i|i-1,3)](51)
P
i|i,3=P
i|i-1,3(1-K
i,3g′
3(z
i|i-1,3))(53)
Wherein
Further, z is expected
i|i-1,3with amount of variation P
i|i-1,3expression formula be:
P
i|i-1,3=P
i-1|i-1,3(55)
By above-mentioned renewal expression formula, obtain the estimated probability density function that the phase III remains continuous time
for
Based on X
3with L
i, 3relation, have
X
3|Φ
i,3=L
i,3+(t
i-x
1-x
2)(57)
Wherein, x
1and x
2the amount obtained in the first two phase process, based on this, directly
Therefore, task time T
mdistribution and expression formula be
Further, at moment t
mthe process of above formula differential is obtained to probability density function task time of phase III:
So far, complete based on relevant environment Monitoring Data information task time probability density function estimation.
Step 2: fusion application historical data information, sets up phased mission systems degenerative process model, estimates the life-span of phased mission systems, obtains the life estimation value of phased mission systems.
(1) Wiener-Hopf equation is used to set up the degenerative process of phased mission systems
For without loss of generality, assuming that the initial reading of degenerative process is Y (0)=0, if initial time is t
i(i=1,2,3), then the monitored parameters evolution process model of time correlation is
Y(t)=y
i+λ(t-t
i)+σB(t-t
i)(60)
Further, consider fusion application historical data information, in model, λ develops into dynamic parameter, uses λ
i=λ
i-1+ η replaces, and wherein η ~ N (0, Q), then degradation model is reconstructed into
Wherein ε
i~ N (0, t
i-t
i-1), λ
iestimated value can be obtained by Kalman filter algorithm.
λ is obtained by Kalman filter algorithm
ithe step of estimated value be:
Step 1: initialization
p
0;
Step 2:t
ithe state estimation in moment is
P
i|i-1=P
i-1|i-1+Q
K
i=(t
i-t
i-1)
2P
i|i-1+σ
2(t
i-t
i-1);
Step 3: upgrade
λ is obtained by Gaussian distribution hypothesis and Bayesian filter principle
iprobability density function be
(2) estimate the life-span of phased mission systems, obtain the life estimation value of phased mission systems
The life-span of equipment is defined by the concept of first-hitting time, when the degradation values Y (t) that formula (60) is determined reaches the failure threshold w preset first, just thinks equipment failure.By the time that the timing definition that equipment life stops is degenerative process Y (t) First excursion failure threshold w, then equipment t
ithe remaining life S in moment
ibe defined as:
S
i=inf{s
i:Y(s
i+t
i)≥w|Y
i}(63)
Order
for the cumulative distribution function of remaining life, predict t by degradation model formula (61)
ithe degeneration probability density function in moment
consideration model profile is normal distribution, has Y (s
i+ t
i) | Y
i:
then t
iprobability density function and the cumulative distribution function of the life expectancy of moment phased mission systems are respectively:
Step 3: according to the reliability definition of phased mission system process model and phased mission systems degenerative process model, comprehensive utilization estimated value task time and life estimation value calculate the reliability reference value of phased mission systems in real time.
Obtain the life-span probability density function of task system
with probability density function task time
after, according to the reliability definition of phased mission system process model and phased mission systems degenerative process model, estimate t
ithe reliability of moment phased mission systems, the probability obtaining the reliability parameter of the n-th stage fill is
The estimation of the reliability parameter value of this kind of complicated phased mission systems of ocean platform can be completed based on formula (65) and (66) two formulas.
Above-described embodiment is used for explaining the present invention, instead of limits the invention, and in the protection domain of spirit of the present invention and claim, any amendment make the present invention and change, all fall into protection scope of the present invention.
Claims (6)
1. the ocean platform phased mission systems reliability estimation methods based on data-driven, it is characterized in that: the steps include: to be detected by the environment of sensor to system, the system environments Monitoring Data that comprehensive utilization sensor is observed, set up contacting between the historical data of individual phased mission systems and current real-time data information, stochastic filtering model is adopted to carry out modeling to stage interval, obtain phased mission system process model, estimate the probability density function of task time, obtain estimated value task time of phased mission systems; Fusion application historical data information, sets up phased mission systems degenerative process model, estimates the life-span of phased mission systems, obtains the life estimation value of phased mission systems;
According to the reliability definition of phased mission system process model and phased mission systems degenerative process model, comprehensive utilization estimated value task time and life estimation value calculate the reliability reference value of phased mission systems in real time.
2. the ocean platform phased mission systems reliability estimation methods based on data-driven according to claim 1, it is characterized in that: set up phased mission system process model and estimate that the concrete steps of probability density function of task time are: suppose that ocean platform offshore and gas development task has N number of stage, make stochastic variable X
nrepresent the n-th phase duration, its span is
make stochastic variable T
mrepresent the T.T. of completion system overall tasks, and have
order
for each phase duration X
n, the probability density function of n=1, L, N, wherein 2≤n≤N; Make Φ
i,nrepresent cut-off current time t
ithe environmental monitoring data information aggregate corresponding with the n-th stage; Adopt probability density function task time in environmental monitoring data estimation of the order 1 ~ N each stage.
3. the ocean platform phased mission systems reliability estimation methods based on data-driven according to claim 2, is characterized in that: the n-th phase process modeling and task time PDF estimation step be: suppose to get current time t by sensor
ithe environmental monitoring data information φ in lower n-th stage
i,n, make L
i,nrepresent t
ithe residue continuous time in stage in moment n-th, then through type L
i+1, n=L
i,n-(t
i+1-t
i), (if L
i,n> t
i+1-t
i) residue L continuous time of the n-th stage subsequent time can be obtained
i+1, n, L
0, nrepresent the residue continuous time of first stage; The consecutive hours area of a room considering task time always on the occasion of, utilize transform Z
i,n=lnL
i,n, obtain a process variable Z
i,n, to ensure L
i,n> 0, adopts a sliding scales parameter to set up φ
i,nand L
i,nstochastic relation model between the two, thus by φ
i,nestimate L
i,n, make φ
i,nbe expressed as Z
i,nfunction, set up the n-th phase process model as follows:
Z
i,n=lnL
i,n(1)
φ
i,n=g
n(Z
i,n)+η
i,n
Wherein, g
n(Z
i,n) be a function undetermined, it is described that the relation between task process and the environmental monitoring data corresponding with the stage; η
i,nthe measuring error of a Normal Distribution, and
Adopt extended Kalman filter to estimate or upgrade Z
i,nconditional probability density function, and ask for further residue continuous time L
i,n; Definition Z
i,nand Z
i+1, nrenewal amount and one-step prediction conditional probability density function amount be respectively Z
i,n| Φ
i,n: N (z
i|i, n, P
i|i, n) and Z
i+1, n| Φ
i,n: N (z
i+1|i, n, P
i+1|i, n), parameter z here
i|i, n, P
i|i, n, z
i+1|i, nand P
i+1|i, ncan be calculated by extended Kalman filter algorithm, its computing formula is as follows:
z
i|i,n=z
i|i-1,n+K
i,n[φ
i,n-g
n(z
i|i-1,n)](2)
Wherein, K
i,nrepresent the gain battle array of Kalman filtering, obtained by following formula:
Wherein
Correspondingly, the renewal equation of state estimation variance matrix can be obtained by following formula:
P
i|i,n=P
i|i-1,n(1-K
i,ng′
n(
zi|i-1,n))(4)
The initial parameter z of extended Kalman filter algorithm
0|0, nand P
0|0, nestimated to obtain by historical data; Calculation expectation z is needed in the renewal equation of above-mentioned state estimation variance matrix
i|i-1, nwith variance P
i|i-1, na step estimated value, consider Z
i,n| Φ
i,n: N (z
i|i, n, P
i|i, n), and Z
i,n=lnL
i,n, obtain these two amounts by lower two formulas:
Consider the transformational relation of normal distribution and lognormal distribution, have L
i, 1| Φ
i, 1: lnN (z
i|i, 1, P
i|i, 1), then based on model expression, obtain and expect z
i|i-1,1with variance P
i|i-1,1expression formula:
P
i|i-1,n=P
i-1|i-1,n(8)
By above-mentioned renewal expression formula, obtain the estimated probability density function that the first stage remains continuous time
for
Again because of
Then directly obtain X by correspondent transform
nat t
ithe distribution in moment
Wherein,
Calculated by formula (9);
Thus, corresponding to t
imoment and X
1subordinate phase duration probability density function expression formula under condition is
Similarly, corresponding to t
imoment and X
npostorder under condition each phase duration probability density function expression formula is
Wherein s=n+1 ..., N;
Further, have
then task time T
mprobability distribution function be
At moment t
mprobability density function task time formula (14) differential process being obtained to the n-th stage is
So far, the environmental monitoring data information of being correlated with by the current generation obtains estimator task time in this stage.
4. the ocean platform phased mission systems reliability estimation methods based on data-driven according to claim 1 or 3, it is characterized in that: use Wiener-Hopf equation to set up phased mission systems degenerative process model, the concrete steps setting up phased mission systems degenerative process model are: assuming that the initial reading of degenerative process is Y (0)=0, if initial time is t
i(i=1,2 ..., N), then the monitored parameters evolution process model of time correlation is
Y(t)=y
i+λ(t-t
i)+σB(t-t
i)(16)
Further, consider fusion application historical data information, in model, λ develops into dynamic parameter, uses λ
i=λ
i-1+ η replaces, and wherein η ~ N (0, Q), then degradation model is reconstructed into
Wherein ε
i~ N (0, t
i-t
i-1), λ
iestimated value obtained by Kalman filtering algorithm,
λ is obtained by Gaussian distribution hypothesis and Bayesian filter principle
iprobability density function be
The life-span of equipment is defined by the concept of first-hitting time, when the degradation values Y (t) that formula (16) is determined reaches the failure threshold w preset first, just thinks equipment failure; By the time that the timing definition that equipment life stops is degenerative process Y (t) First excursion failure threshold w, then equipment t
ithe remaining life S in moment
ibe defined as:
S
i=inf{s
i:Y(s
i+t
i)≥w|Y
i}(19)
Order
for the cumulative distribution function of remaining life, predict t by degradation model formula (17)
ithe degeneration probability density function in moment
consideration model profile is normal distribution, has
Then t
iprobability density function and the cumulative distribution function of the life expectancy of moment phased mission systems are respectively:
5. the ocean platform phased mission systems reliability estimation methods based on data-driven according to claim 4, is characterized in that: obtain λ by Kalman filtering algorithm
i(i=1,2 ..., N) the step of estimated value be:
Initialization
p
0;
T
ithe state estimation in moment is
P
i|i-1=P
i-1|i-1+Q
K
i=(t
i-t
i-1)
2P
i|i-1+σ
2(t
i-t
i-1);
Upgrade
6. the ocean platform phased mission systems reliability estimation methods based on data-driven according to claim 4, is characterized in that: the life-span probability density function obtaining task system
with probability density function task time
after, according to the reliability definition of phased mission system process model and phased mission systems degenerative process model, estimate t
ithe reliability of moment phased mission systems, the probability obtaining the reliability parameter of the n-th stage fill is
The estimation of the reliability parameter value of ocean platform complex phased mission systems can be completed based on formula (21) and (22).
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510925359.4A CN105426692B (en) | 2015-12-10 | 2015-12-10 | Ocean platform phased mission systems reliability estimation methods based on data-driven |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510925359.4A CN105426692B (en) | 2015-12-10 | 2015-12-10 | Ocean platform phased mission systems reliability estimation methods based on data-driven |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105426692A true CN105426692A (en) | 2016-03-23 |
CN105426692B CN105426692B (en) | 2018-06-05 |
Family
ID=55504900
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510925359.4A Active CN105426692B (en) | 2015-12-10 | 2015-12-10 | Ocean platform phased mission systems reliability estimation methods based on data-driven |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105426692B (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107480440A (en) * | 2017-08-04 | 2017-12-15 | 山东科技大学 | A kind of method for predicting residual useful life for modeling of being degenerated at random based on two benches |
CN110134675A (en) * | 2019-05-23 | 2019-08-16 | 大连海事大学 | A kind of data cleaning method and system towards oceanographic data stream |
CN111144606A (en) * | 2019-05-17 | 2020-05-12 | 深圳市德塔防爆电动汽车有限公司 | Safety failure risk prediction method for electric vehicle and electric vehicle |
CN111815107A (en) * | 2020-05-22 | 2020-10-23 | 中国人民解放军92942部队 | Task reliability modeling method for representing time elements |
CN113221372A (en) * | 2021-05-26 | 2021-08-06 | 河北工业大学 | BDD-based industrial robot PMS system reliability analysis method |
CN113435660A (en) * | 2021-07-13 | 2021-09-24 | 中国石油大学(华东) | Underwater Christmas tree re-prediction method integrating Kalman filtering and Bayesian network |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080082194A1 (en) * | 2006-09-29 | 2008-04-03 | Fisher-Rosemount Systems, Inc. | On-line multivariate analysis in a distributed process control system |
CN104166787A (en) * | 2014-07-17 | 2014-11-26 | 南京航空航天大学 | Aero-engine remaining life prediction method based on multi-stage information fusion |
CN104615866A (en) * | 2015-01-21 | 2015-05-13 | 北京航空航天大学 | Service life prediction method based on physical statistic model |
-
2015
- 2015-12-10 CN CN201510925359.4A patent/CN105426692B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080082194A1 (en) * | 2006-09-29 | 2008-04-03 | Fisher-Rosemount Systems, Inc. | On-line multivariate analysis in a distributed process control system |
CN104166787A (en) * | 2014-07-17 | 2014-11-26 | 南京航空航天大学 | Aero-engine remaining life prediction method based on multi-stage information fusion |
CN104615866A (en) * | 2015-01-21 | 2015-05-13 | 北京航空航天大学 | Service life prediction method based on physical statistic model |
Non-Patent Citations (2)
Title |
---|
彭宇,等: "数据驱动故障预测和健康管理综述", 《仪器仪表学报》 * |
谢红卫,等: "时变结构多阶段任务型***的可靠度研究", 《国防科技大学学报》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107480440A (en) * | 2017-08-04 | 2017-12-15 | 山东科技大学 | A kind of method for predicting residual useful life for modeling of being degenerated at random based on two benches |
CN107480440B (en) * | 2017-08-04 | 2020-01-21 | 山东科技大学 | Residual life prediction method based on two-stage random degradation modeling |
CN111144606A (en) * | 2019-05-17 | 2020-05-12 | 深圳市德塔防爆电动汽车有限公司 | Safety failure risk prediction method for electric vehicle and electric vehicle |
CN110134675A (en) * | 2019-05-23 | 2019-08-16 | 大连海事大学 | A kind of data cleaning method and system towards oceanographic data stream |
CN110134675B (en) * | 2019-05-23 | 2022-11-01 | 大连海事大学 | Data cleaning method and system for ocean data stream |
CN111815107A (en) * | 2020-05-22 | 2020-10-23 | 中国人民解放军92942部队 | Task reliability modeling method for representing time elements |
CN111815107B (en) * | 2020-05-22 | 2022-11-01 | 中国人民解放军92942部队 | Task reliability modeling method for representing time elements |
CN113221372A (en) * | 2021-05-26 | 2021-08-06 | 河北工业大学 | BDD-based industrial robot PMS system reliability analysis method |
CN113221372B (en) * | 2021-05-26 | 2023-02-24 | 河北工业大学 | BDD-based industrial robot PMS system reliability analysis method |
CN113435660A (en) * | 2021-07-13 | 2021-09-24 | 中国石油大学(华东) | Underwater Christmas tree re-prediction method integrating Kalman filtering and Bayesian network |
CN113435660B (en) * | 2021-07-13 | 2022-07-08 | 中国石油大学(华东) | Underwater Christmas tree re-prediction method integrating Kalman filtering and Bayesian network |
US11768976B2 (en) * | 2021-07-13 | 2023-09-26 | China University Of Petroleum (East China) | Subsea Christmas tree re-prediction method integrating Kalman filter and Bayesian network |
Also Published As
Publication number | Publication date |
---|---|
CN105426692B (en) | 2018-06-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105426692A (en) | Ocean platform multi-stage task system reliability estimation method based on data drive | |
Prakash et al. | Initial service life data towards structural health monitoring of a concrete arch dam | |
CN101789039B (en) | Calculation method for availability ratio and optimal repair cycle of relay protection device | |
CN103413023B (en) | A kind of multi-state system dynamic reliability assessment method | |
CN106372299B (en) | A kind of pipeline method for predicting residual useful life based on inverse Gaussian process and EM-PF | |
CN101382473B (en) | EWMA control chart method for bridge structure safety alarm | |
Xia et al. | Long‐term performance monitoring and assessment of concrete beam bridges using neutral axis indicator | |
Titova et al. | Lifetime of earth dams | |
CN102695943A (en) | Universal sensor self-diagnosis device and diagnosis method therefor | |
CN105046445A (en) | Flood prevention dispatching decision support system | |
Cai et al. | A dynamic-Bayesian-networks-based resilience assessment approach of structure systems: Subsea oil and gas pipelines as A case study | |
CN113435660B (en) | Underwater Christmas tree re-prediction method integrating Kalman filtering and Bayesian network | |
Huang et al. | An integrated risk sensing system for geo-structural safety | |
Kumar et al. | Methods for risk management of mining excavator through FMEA and FMECA | |
Huang et al. | Quantitative geotechnical risk management for tunneling projects in China | |
Iliopoulos et al. | Continuous fatigue assessment of an offshore wind turbine using a limited number of vibration sensors | |
Loughney et al. | Applications of industrial iot and wsns in o&m programmes for offshore wind farms | |
May et al. | Integrating structural health and condition monitoring: a cost benefit analysis for offshore wind energy | |
CN109977563B (en) | Submarine oil and gas pipeline residual life prediction method based on TV-HSMM | |
CN105354615A (en) | Variable-cost real-time prediction maintenance and spare part ordering method and system for blowout prevention valve | |
Phoon et al. | Special Issue on “Data-Centric Geotechnics for Practice” | |
KR100397191B1 (en) | Method of Evaluating Abnormal Groundwater Head of Observation Holes in Rock Surrounding Underground Highly Pressurized Gas Storage Caverns. | |
CN112363476A (en) | Underground drainage monitoring system and method based on positioning technology | |
Desnitsky et al. | Monitoring the state of materials in cyberphysical systems: water supply case study | |
Xiu et al. | Lifetime prognosis of blowout preventer valve for offshore platform based on Wiener process |
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 |